Abstract
The present article deals with the buckling response of functionally graded multilayer graphene platelet-reinforced composite (FG-GPL RC) rectangular plates with circular/elliptical cutouts resting on a Winkler-type elastic foundation under uniaxial and biaxial normal and shear loads. Rule of mixtures and the Halpin–Tsai approach are applied to obtain the effective Poisson’s ratio, mass density, and elastic modulus of the reinforced composite. The governing equations are developed by applying the third-order shear deformation plate theory. Then, the finite element procedure is used to solve the problem. Four different types of graphene platelet distributions, namely, UD, FG-X, FG-V, and FG-O, are considered. A broad range of factors such as plate aspect ratio, plate slenderness ratio, applying uniaxial and biaxial normal and shear loads to the plate, several Winkler elastic foundation stiffness parameters, different displacement boundary conditions, the effect of size of the circular cutout and orientation of the elliptical cutout, and the influence of GPL weight fraction are discussed in several tabular and graphical data in detail.
1 Introduction
The increasing usage of composite materials, such as carbon nanotube (CNT) and graphene-reinforced composite (GRC) structures, due to their countless applications and superior material properties, has led to several studies and articles investigating the different aspects and behaviors of these appreciable materials. GRC structures have several advantages, such as lightweight, flexibility, transparency, energy absorbency, and electrical and thermal conductivity, and their foreseeable wide application in different industries has attracted significant attention in recent years. Adding a low content of nanofiller reinforcements with a specific distribution pattern to the matrix can surprisingly improve the structural behavior. Hence, several research studies have been carried out on the structural analyses of shell- and plate-type structures reinforced with graphene platelets (GPLs). Those studies mentioned here are related to the behavior of structures that are reinforced by GPLs. In detail, by applying the first-order shear deformation theory (FSDT), showed that the weight fraction of nanofillers and GPL distribution pattern are two main factors affecting the bending and buckling behaviors of the FG-GPL RC plates. investigated FG-GRC-laminated cylindrical panels and noticed that the piece-wise FG pattern of graphene reinforcement may enhance the buckling forces and the post-buckling strength. The governing equations for the post-buckling responses of the panels were proposed based on Reddy’s higher-order shear deformation shell theory and von Kármán strain–displacement relationships. employed the higher-order shear deformation plate theory (HSDT) to present the thermal post-buckling of FG-GPL-laminated plates supporting on an elastic medium. They presented that the critical buckling temperature of the plates may be enhanced by the FG graphene reinforcement distribution pattern of the structure. examined the buckling and post-buckling responses of the uniaxially and biaxially compressed FG-GPL RC plates and revealed that by adding only 1.0 wt% GPLs, the FG-X structure achieved the highest critical buckling force up to 5.55 times that of the neat epoxy plate. In their study, theoretical formulations were proposed based on the FSDT and von Kármán-type nonlinear kinematics and included the effect of an initial geometric imperfection. A two-step perturbation technique was employed to determine the asymptotic post-buckling solutions, and the biaxial compressive post-buckling equilibrium paths of both perfect and imperfect plates simply supported on all edges. applied the Hamilton principle, FSDT, and finite element method (FEM) to study the influence of the volume fraction of CNTs on the displacements of annular sector plate and observed that enhancing the volume fraction of CNTs from 0.11 to 0.17 causes a reduction of more than 36% in the transverse displacements. used the Hamilton principle based on FSDT and FEM to derive the governing equations of motion of FG porous annular sector plates reinforced by GPLs for three various types of porosity patterns. They revealed that in the case of the GPLA pattern, the transverse displacement is decreased by 91% by adding only 1 wt% of GPLs to the metal matrix. applied FSDT of shells and Donnell kinematic relationships and Hamilton principle to derive the governing equations for the free vibrations and the boundary conditions of GPL-RC cylindrical shells. They showed that when the inner and outer layers of the GPL-RC cylindrical shells have the maximum content of reinforcements, maximum frequencies are obtained. proposed the NURBS procedure based on the four-variable refined plate theory to analyze FG-GPL RC plates and proved that the structures with FG-O and FG-X patterns lead to the minimum and maximum amounts of natural frequencies and buckling loads, respectively. presented the free vibration behavior of thick, multilayer composite plates with different distribution patterns of GPLs and obtained the same results as . They applied a quasi-3D plate model which captures the thickness stretching effects and non-uniform shear strains through the thickness. With the aid of the Navier solution, suitable for plates with all edges simply supported, Fourier expansions were implemented for the essential variables of the displacement field. Closed-form expressions were provided to obtain the natural frequencies of FG-GPL RC plates. investigated the buckling and post-buckling behaviors of functionally graded multilayer nanocomposite beams reinforced with a low content of GPLs resting on an elastic foundation. The nonlinear governing equations of the beam on an elastic foundation were derived within the framework of the first-order shear deformation beam theory and then converted into a nonlinear algebraic system by using the differential quadrature method. They showed that GPLs with a bigger surface area and comprising fewer single-GPL laminates may provide better reinforcing influence; however, for GPL aspect ratios and width-to-thickness ratios more than 4 and 1000, critical buckling force and post-buckling path will be almost identical. Furthermore, many studies have been implemented to show that functionally graded nanofiller distribution patterns also play a key role in the structural response of composite structures (; ; ; ; ; ; ; ; ; ; ; ; ; ; ).
Many studies have investigated the effect of an elastic foundation on the structural behavior of composite structures to show that the existence of an elastic foundation can exhibit better performance of the structure. For example, achieved highly accurate results by employing a three-dimensional elasticity-graded element to study the buckling response of FG annular sector plates. Their results revealed that the size and location of buckling deformations are mainly influenced by how the plate is partially rested on the elastic foundation. observed that applying an elastic medium as the foundation leads to higher buckling forces and causes narrower, local, and shorter buckling deformations. investigated the buckling response of a GPL sheet resting on an elastic foundation subjected to thermal conditions and shear loads using the differential quadrature method (DQM). studied the buckling behavior of GPL sheets resting on an elastic Pasternak and Winkler medium and revealed that critical axial buckling forces of the rested GPL sheet for the Pasternak foundation were greater than those obtained for the Winkler foundation.
Structural members, including plates, beams, and shells, are often subjected to different load conditions that may cause the buckling phenomenon. Several research studies have been conducted to investigate the buckling problem of composite structures under different loading conditions, including uniaxial and biaxial normal and shear loads. For example, presented the buckling of heterogeneous FGM plates subjected to different normal and shear forces by applying a 3D elasticity approach. applied a 3D elasticity-based graded FEM to study the shear buckling analysis of FGM annular sector plates. applied the Gurtin–Murdoch surface theory of elasticity to analyze the buckling behavior of FGM skew nanoplates subjected to shear loads.
Composite plates find application as a structural member in various fields, like automobiles, space vehicles, submarines, and many others, where a high strength and stiffness-to-weight ratio are desired. Sometimes, holes and cutouts of various shapes are unavoidable in these types of plates for their proper functioning. These structural members are usually subjected to in-plane loading. With geometrical discontinuities, these composite plates with holes are more vulnerable due to an increase in stress near discontinuities (holes), which results in the buckling phenomenon. Therefore, the buckling analysis of these structures for better safety and evaluation of the load-bearing capacity of these composite plates is essential. Limited studies have been conducted on the structural analyses of these composite structures to find the effect of cutout size and shape on their behavior. For example, conducted a numerical study using FEM to show the effects of different shapes of cutouts on the thermal buckling response of graphite/epoxy-laminated plates. They found that for a plate with an elliptical cutout, the values of the critical buckling temperatures are enhanced by increasing the cutout orientation from a horizontal case to a vertical case. analyzed the buckling of FG-GPL RC plates with a circular hole subjected to normal loads based on the FSDT and FE approaches for three different GPL distribution patterns. showed that stresses around the circular cutout of FGM plates under biaxial tensions may be considerably decreased by the FG distribution of the material properties. applied FEM to provide mechanical stress distribution of FGM plates around the circular cutout and natural frequencies of the plate. utilized computational modeling to study the effect of circular cutouts’ diameter and arrangements on the buckling response of FGM plates. The obtained results denoted that horizontal arrangements exhibit the highest critical buckling loads and best buckling performance.
Reviewing relevant articles showed that several research studies have been conducted on the behaviors of the composite structures reinforced by CNTs and GPLs. However, fewer studies have been conducted to study the effects of circular/elliptical cutouts on the buckling response of composite structures, especially FG-GRC plates. Due to the various industrial applications of graphene platelet-reinforced composite plates in mechanical and aerospace engineering, the need to fill the void of how circular and elliptical cutouts influence the buckling response of GPL RC plates is sensed. Particularly, the effect of shear loads and elliptical-shaped cutouts on the buckling response of plates is less investigated. Furthermore, plates on an elastic foundation are common structural elements employed in many civil engineering applications, such as foundations, storage tanks, swimming pools, floor systems of buildings and highways, and airfield pavements. Therefore, the influence of biaxial normal and shear loads and elastic foundation on the buckling response of FG-GPL RC plates with circular/elliptical cutouts is studied for the first time in the present work. In addition, most of the previous studies are based on FSDT and CLPT, which are suitable for moderately thick and thin plates, respectively. FSDT overestimates the buckling loads for thick plates rather than the third-order shear deformation plate theory (TSDT). Hence, in the present study, the governing equations are derived applying the TSDT to investigate the influence of a broad range of factors such as plate aspect ratio, plate width/thickness ratio, size and shape of the cutout, uniaxial and biaxial normal and shear loads, several Winkler elastic foundation stiffness parameters, different boundary conditions, and GPL nanofiller weight fraction on the buckling loads of FG-GPL RC plates. The influence of four major graphene platelet (GPL) distribution patterns, namely, UD, FG-X, FG-V, and FG-O, on the buckling response of plates is investigated. Rule of mixtures and the Halpin–Tsai approach are applied to obtain the effective Poisson’s ratio, mass density, and elastic modulus of the reinforced composite. The finite element procedure is employed to solve the problem.
2 Problem modeling
2.1 Geometry and material properties of an FG-GPL RC plate
We combine the advantages of both FGM and GPL using a functionally graded graphene platelet-reinforced composite (FG-GPL RC) structure. GPLs and isotropic epoxy are considered the reinforcement and the matrix of the composite media, respectively. An FG-GPL RC plate with length a, width b, total thickness h, and a circular cutout with radius r, as shown in Figure 1, is assumed in this paper. The plate is resting on a Winkler elastic foundation, where kw indicates the stiffness of the Winkler elastic foundation. The plate is constructed of layers with identical thickness t. GPLs in each individual layer are uniformly distributed, but their volume fraction changes from layer to layer, which leads to a functionally graded pattern. Therefore, each individual layer of the GPL RC plate is isotropic and homogeneous. Figure 2 shows the five uniform or functionally graded distribution patterns of GPL nanofillers. UD denotes the uniform distribution, and FG-X, FG-O, FG-V, and FG-Λ show the functionally graded patterns of GPLs. The UD distribution 1) denotes a homogeneous and isotropic multilayer plate in which nanofillers are uniformly dispersed across the thickness of the plate. The GPL volume fraction is dispersed symmetrically and linearly in FG-X 2) and FG-O 3) distributions, in which the GPL amount gradually changes from the top or bottom surfaces to the middle layer. The unsymmetrical patterns are FG-V 4) and FG-Λ 5) distributions. An adequately large number of layers, , leads to a gentle and connected distribution pattern along the plate thickness. The volume fraction of the kth layer for different GPL patterns is given as follows (:
FIGURE 1
FIGURE 2
where is total number of layers and assumed to be an even number. is the average magnitude of the volume fraction of GPLs and can be evaluated based on the weight fraction of GPLs , density of the matrix , and GPLs .
In the aforementioned equations, the subscripts (m) and (GPL) stand for the matrix and GPLs, respectively. The effective elastic modulus of the nanocomposite plate is obtained through the Halpin–Tsai model, which is dependent on the geometry of the nanoplatelets (; ). Hence, the elasticity modulus of the kth layer of the plate, which is denoted by , is regarded asIn Eq. 3, the elasticity modulus of the matrix is shown by . The two parameters and in Eq. 3 depend on the elastic properties of the constituents and the geometry of the platelets:
The other two parameters, and , in Eqs 4, 5 depend on the geometrical dimensions of the nanofillers:where , , and are assumed to be the length, width, and thickness of GPLs, respectively.
The simple rule of mixture estimation is used to evaluate the Poisson’s ratio and the mass density of the kth layer of the nanocomposite plate (; ; ):
wherewhere , and , are the Poisson’s ratio and mass density of the GPL and matrix, respectively.
In addition, the rigidity modulus of the kth layer is as follows:
2.2 Governing equations of an FG-GL RC plate
The governing equations of the plate resting on an elastic foundation are derived applying the third-order shear deformation plate theory. Hence, the displacement field is assumed aswhere , , and represent the displacements at (the midplane of a plate); and are the transverse normal rotations of the and axes. The finite element procedure is used to estimate the eigenvalue buckling forces of the FG-GPL RC plate. A four-node element with seven DOFs per node is employed for modeling the plate. The displacement field of an individual node is as follows:
Among the displacement components, is estimated by the Hermitian shape functions and the other components by the Lagrangian functions. Hence, the displacement components may be presented as follows:where is the Lagrangian shape function and is the Hermitian function. The components of and are presented in Appendix A.
The displacement field of an element is approximated by using its nodal values:where is the shape function matrix; and are expressed as follows:
The strains of the plate when it is subjected to in-plane pre-buckling loads may be expressed aswith
By substituting Eq. 15 into Eqs 24, 25, the strain components may be presented aswith
The stress–strain relationships based on Hooke’s law are as follows:with
The elastic strain energy of the plate and Winkler elastic foundation is written as follows:
or in the matrix formwhere the stiffness matrices corresponding to the strain energy of deformation of element and elastic foundation are given as
where
The geometric strain energy due to in-plane pre-buckling loads is obtained as
By substituting the geometric strain into Eq. 41, we obtainwhere
Then, Eq. 42 becomes
with
For the buckling problem, after assembly of element matrices, we derive the following equation:where , , and are the global stiffness matrix, stiffness matrix of the elastic foundation, and global geometric stiffness matrix, respectively. Eq. 47 should be solved to obtain the buckling mode shapes and buckling forces .
3 Parametric and comparison studies
A FG-GPL RC plate, as shown in Figure 1, with width is considered. As studied in , , and , results and values for free vibration and buckling analysis of a functionally graded multilayer GPL RC plate with a total number of layers ≥10 provide accurate and acceptable results. In the present study, to achieve more accurate results, the total number of layers is considered to be = 20. Other required data, geometric ratios, and boundary conditions are given in each table. GPL nanofillers considered in this paper have a length of = 2.5 μm, width of = 1.5 μm, and thickness of = 1.5 nm. The other required material properties such as elastic modulus and Poisson’s ratio of epoxy matrix and GPL reinforcement are provided in Table 1.
TABLE 1
| Property name | Epoxy matrix | GPL |
|---|---|---|
| Modulus of elasticity | 3 | 1,010 |
| Poisson’s ratio | 0.34 | 0.186 |
| Density | 1,200 | 1,062.5 |
Mechanical properties of the epoxy matrix and GPLs.
Three boundary conditions for different loadings are considered:
3.1 Verification of results
In order to validate the accuracy of the results of the present study, the critical buckling load ratios are recalculated and compared as shown in Table 2 with those from . A comparative study is carried out for the critical buckling load ratios of a simply supported FG-GPL RC plate without cutouts and with the same material properties listed in Table 1. Here, and are critical buckling loads of the plates with and without GPLs, respectively. Geometric dimensions of the plate are , , and . The number of layers, , and 1.0 wt% GPLs are used for comparing results. It can be seen that the present results are in good agreement with those obtained by .
TABLE 2
| Distribution pattern | Present | Difference % | |
|---|---|---|---|
| UD | 4.3236 | 4.2847 | 0.9077 |
| FG-X | 5.5544 | 5.5439 | 0.1893 |
| FG-O | 3.0476 | 2.9540 | 3.1686 |
Comparison of the critical buckling load ratio () for a simply supported FG-GPL RC plate under uniaxial compression.
3.2 Numerical results and discussion for buckling loads of an FG-GPL RC plate with circular/elliptical cutouts
In this section, several numerical examples are demonstrated to study the effects of different factors on the buckling behavior of an FG-GPL RC plate with a cutout.
The effect of different thickness ratios () and GPL patterns on the buckling response of a multilayered FG-GPL RC plate with a circular cutout under uniaxial compression is shown in Table 3. The cutout radius is assumed to be and the plate aspect ratio is , while the plate is fully clamped at all edges. The first four buckling loads are reported in Table 3. This table reveals that with the same volume fraction and same width-to-thickness ratio, the FG-X distribution pattern results in the highest buckling values, while the FG-O distribution pattern results in the lowest values. Hence, UD and FG-V distribution patterns result in intermediate values. In addition, results show that when the plate slenderness ratio increases, the buckling loads decrease, as expected. The influence of GPL weight fraction and distribution pattern of nanofillers on the buckling loads of the FG-GPL RC plate with a circular cutout is provided in Table 4. The buckling loads of the pure epoxy plate are tabulated for comparison purposes. As can be seen in the table, by dispersing only a low amount of GPL nanofillers, the buckling loads can remarkably improve. By dispersion of only 0.5 wt% and 1.0 wt% GPLs for the FG-X plate, the buckling load enhanced 2.97 and 4.84 times, respectively. The effects of two different boundary conditions and GPL patterns on the buckling loads of the plate under uniaxial compression are compared in Table 5. This table shows the first four buckling loads of the plate for , , and circular cutout radius . As can be seen in this table, buckling forces are higher for the boundary conditions that restrict more degrees of freedom (DOFs). Values tabulated in Table 6 are provided to indicate the effect of the ratio and GPL pattern on the buckling load results of the FG-GPL RC clamped plate with a circular cutout under uniaxial load. Results denote that by increasing the plate aspect ratio from 1 to 2, the buckling loads enhanced up to 99.6% for the FG-X plate, but by changing the aspect ratio from 2 to 4, the buckling loads decreased by 7.8%. This can be justified by investigating the pre-buckling stresses of the plates with different aspect ratios. In this case, for plates with a/b = 4, the magnitude of the compressive stresses around the hole is more than the plates with a/b = 2. Results of Tables 3–6 denote that the stiffness of the outer layers generally has the most remarkable effect on the stiffness of the composite structure. Therefore, as shown previously, the structure with the FG-X gradient pattern shows the highest stiffness and buckling capabilities, and UD, FG-V, and FG-O patterns in the sequence have lower values.
TABLE 3
| b/h | UD | FG-X | FG-O | FG-V | |
|---|---|---|---|---|---|
| 10 | Mode 1 | 426.790 | 555.260 | 269.750 | 348.840 |
| Mode 2 | 604.830 | 710.790 | 401.040 | 503.030 | |
| Mode 3 | 747.950 | 805.940 | 522.560 | 633.980 | |
| Mode 4 | 923.590 | 999.980 | 627.820 | 775.960 | |
| 15 | Mode 1 | 29.944 | 33.861 | 20.671 | 26.038 |
| Mode 2 | 31.556 | 34.682 | 22.236 | 27.738 | |
| Mode 3 | 32.177 | 35.561 | 22.589 | 28.151 | |
| Mode 4 | 39.842 | 40.566 | 29.257 | 35.412 | |
| 20 | Mode 1 | 13.433 | 16.394 | 8.926 | 11.414 |
| Mode 2 | 14.379 | 17.120 | 9.700 | 12.327 | |
| Mode 3 | 14.650 | 17.541 | 9.853 | 12.522 | |
| Mode 4 | 18.809 | 21.312 | 12.993 | 16.352 | |
| 30 | Mode 1 | 4.156 | 5.384 | 2.686 | 3.458 |
| Mode 2 | 4.516 | 5.755 | 2.945 | 3.779 | |
| Mode 3 | 4.589 | 5.874 | 2.986 | 3.835 | |
| Mode 4 | 6.041 | 7.517 | 3.980 | 5.094 |
Effect of the b/h ratio on the buckling load of the FG-GPL RC plate with a circular cutout (; , CCCC, uniaxial load, and 1.0 wt%).
TABLE 4
| Weight fraction (wt%) | UD | FG-X | FG-O | FG-V | |
|---|---|---|---|---|---|
| 0.0 | Mode 1 | 3.390 | 3.390 | 3.390 | 3.390 |
| Mode 2 | 3.629 | 3.629 | 3.629 | 3.629 | |
| Mode 3 | 3.698 | 3.698 | 3.698 | 3.698 | |
| Mode 4 | 4.747 | 4.747 | 4.747 | 4.747 | |
| 0.5 | Mode 1 | 8.408 | 10.081 | 6.192 | 7.635 |
| Mode 2 | 9.000 | 10.603 | 6.709 | 8.221 | |
| Mode 3 | 9.170 | 10.847 | 6.818 | 8.358 | |
| Mode 4 | 11.773 | 13.415 | 8.950 | 10.856 | |
| 1.0 | Mode 1 | 13.433 | 16.394 | 8.926 | 11.414 |
| Mode 2 | 14.379 | 17.120 | 9.700 | 12.327 | |
| Mode 3 | 14.650 | 17.541 | 9.853 | 12.522 | |
| Mode 4 | 18.809 | 21.312 | 12.993 | 16.352 |
Effect of the GPL weight fraction and GPL pattern on the buckling load of the FG-GPL RC plate with a circular cutout (; ; , CCCC, and uniaxial load).
TABLE 5
| BC | UD | FG-X | FG-O | FG-V | |
|---|---|---|---|---|---|
| SSSS | Mode 1 | 12.807 | 15.630 | 8.452 | 10.666 |
| Mode 2 | 13.740 | 16.346 | 9.215 | 11.564 | |
| Mode 3 | 13.896 | 16.637 | 9.279 | 11.623 | |
| Mode 4 | 17.854 | 20.259 | 12.234 | 15.182 | |
| CCCC | Mode 1 | 13.433 | 16.394 | 8.926 | 11.414 |
| Mode 2 | 14.379 | 17.120 | 9.700 | 12.327 | |
| Mode 3 | 14.650 | 17.541 | 9.853 | 12.522 | |
| Mode 4 | 18.809 | 21.312 | 12.993 | 16.352 |
Effect of different boundary conditions and GPL patterns on the buckling load of the FG-GPL RC plate with a circular cutout (; ; , uniaxial load, and 1.0 wt%).
TABLE 6
| a/b | UD | FG-X | FG-O | FG-V | |
|---|---|---|---|---|---|
| 1 | Mode 1 | 13.433 | 16.394 | 8.926 | 11.414 |
| Mode 2 | 14.379 | 17.120 | 9.700 | 12.327 | |
| Mode 3 | 14.650 | 17.541 | 9.853 | 12.522 | |
| Mode 4 | 18.809 | 21.312 | 12.993 | 16.352 | |
| 2 | Mode 1 | 26.083 | 32.723 | 17.075 | 21.900 |
| Mode 2 | 26.417 | 32.997 | 17.337 | 22.223 | |
| Mode 3 | 34.562 | 41.668 | 23.020 | 29.355 | |
| Mode 4 | 35.932 | 42.732 | 24.118 | 30.670 | |
| 4 | Mode 1 | 24.056 | 30.361 | 15.707 | 20.163 |
| Mode 2 | 24.057 | 30.361 | 15.707 | 20.163 | |
| Mode 3 | 30.854 | 37.645 | 20.446 | 26.086 | |
| Mode 4 | 30.854 | 37.646 | 20.447 | 26.086 |
Effect of the plate aspect ratio and GPL pattern on the buckling load of the FG-GPL RC plate with a circular cutout (; , CCCC, uniaxial load, and ).
Table 7 presents the buckling loads for three different radiuses of the circular cutout of the FG-GPL RC plate under uniaxial load. As expected, with the reduction in the cutout radius, the buckling loads are increased. For example, by increasing the cutout radius-to-width ratio ( from 0.1 to 0.3, the buckling loads are decreased by 22.3%. This is due to the fact that for larger cutouts, the pre-buckling compressive stresses are more than smaller cutouts. Table 8 lists the first four buckling loads for the FG-GPL RC plate under uniaxial load with a circular cutout resting on the Winkler-type elastic foundation. Values are provided for several elastic foundation stiffness parameters and a plate with CCCC supports. The comparison of the values in Table 8 with those without an elastic foundation (Table 3) denotes that applying an elastic foundation postpones the buckling phenomenon. For example, the buckling loads are increased by 42.6% for the FG-GPL RC plate when the elastic foundation stiffness is increased from 0 to 1 .
TABLE 7
| Cutout size | UD | FG-X | FG-O | FG-V | |
|---|---|---|---|---|---|
| ×0.5 (r/b = 0.1) | Mode 1 | 15.170 | 19.293 | 9.873 | 12.715 |
| Mode 2 | 17.610 | 21.667 | 11.664 | 14.881 | |
| Mode 3 | 26.881 | 30.675 | 18.462 | 23.207 | |
| Mode 4 | 30.533 | 36.491 | 20.428 | 26.123 | |
| ×1 (r/b = 0.2) | Mode 1 | 13.433 | 16.394 | 8.926 | 11.414 |
| Mode 2 | 14.379 | 17.120 | 9.700 | 12.327 | |
| Mode 3 | 14.650 | 17.541 | 9.853 | 12.522 | |
| Mode 4 | 18.809 | 21.312 | 12.993 | 16.352 | |
| ×1.5 (r/b = 0.3) | Mode 1 | 11.871 | 13.408 | 8.252 | 10.326 |
| Mode 2 | 12.027 | 13.556 | 8.355 | 10.488 | |
| Mode 3 | 12.112 | 13.571 | 8.413 | 10.573 | |
| Mode 4 | 12.122 | 13.678 | 8.471 | 10.577 |
Effect of the circular cutout radius and GPL pattern on the buckling load of the FG-GPL RC plate with a circular cutout (; , CCCC, uniaxial load, and 1.0 wt%).
TABLE 8
| Kw (N/m3) | UD | FG-X | FG-O | FG-V | |
|---|---|---|---|---|---|
| 0 | Mode 1 | 13.433 | 16.394 | 8.926 | 11.414 |
| Mode 2 | 14.379 | 17.120 | 9.700 | 12.327 | |
| Mode 3 | 14.650 | 17.541 | 9.853 | 12.522 | |
| Mode 4 | 18.809 | 21.312 | 12.993 | 16.352 | |
| 1.00E+08 | Mode 1 | 14.983 | 17.704 | 10.486 | 13.033 |
| Mode 2 | 15.420 | 18.059 | 10.761 | 13.415 | |
| Mode 3 | 16.157 | 18.860 | 11.400 | 14.095 | |
| Mode 4 | 19.477 | 21.861 | 13.714 | 17.074 | |
| 5.00E+08 | Mode 1 | 18.446 | 20.832 | 13.663 | 16.504 |
| Mode 2 | 18.622 | 20.980 | 13.819 | 16.679 | |
| Mode 3 | 20.249 | 22.427 | 15.160 | 18.215 | |
| Mode 4 | 21.706 | 23.701 | 15.985 | 19.426 | |
| 1.00E+09 | Mode 1 | 21.313 | 23.385 | 16.344 | 19.400 |
| Mode 2 | 21.544 | 23.586 | 16.533 | 19.628 | |
| Mode 3 | 23.203 | 24.948 | 17.739 | 21.130 | |
| Mode 4 | 23.870 | 25.475 | 18.103 | 21.669 |
Effect of the Winkler elastic foundation and GPL pattern on the buckling load of the FG-GPL RC plate with a circular cutout (; ; , CCCC, uniaxial load, and 1.0 wt%).
Buckling loads of the FG-GPL RC plate under uniaxial and biaxial normal loads are compared with those resting on an elastic foundation in Table 9 for two different plate aspect ratios and four GPL patterns. As can be seen in the table, for the plate aspect ratio of 1, the buckling loads of uniaxial and biaxial compressions are almost identical for the plate without an elastic foundation. For the higher plate aspect ratios (), the buckling load of the plate under biaxial load is considerably lower compared to those under uniaxial load. For all aspect ratios, by adding an elastic foundation to the plate, the buckling loads are increased for both uniaxial and biaxial compressions. Values provided in Table 10 are tabulated with the intent of comparing the critical buckling loads of the FG-GPL RC plate under uniaxial load with different shapes of cutout with/without an elastic foundation. The effect of the elliptical cutout orientation of the plate under uniaxial load is also investigated. The elliptical cutout is considered to have the same area as the circular cutout with . Here, c and d are the semi-major and semi-minor radius of the elliptical cutout, and their radius-to-width ratios are assumed to be 0.3 and , respectively. As seen in Table 10, if the major radius of the elliptical cutout is parallel to the direction of the applied load, the critical buckling load of the UD pattern is approximately 2.39 times lower compared to those when the semi-major radius of the elliptical cutout is perpendicular to the applied load. Furthermore, in this case, the buckling loads are lower than the case of the plate with the circular cutout. Results show that the influence of the elastic foundation on the plate with an elliptical cutout is more remarkable compared to the plate with a circular cutout. For example, when the plate with an elliptical cutout is embedded in an elastic foundation, the critical buckling load of the plate with the UD pattern is approximately 1.73 times greater compared with the plate without an elastic foundation, while for the plate with a circular cutout, it is approximately 1.37 times higher. Buckling loads of the FG-GPL RC plate under shear loads are compared with those resting on the elastic foundation, as shown in Table 11, for two different plate aspect ratios, different shapes of cutouts, and four GPL patterns. By increasing the aspect ratio from 1 to 2, the critical shear buckling load of the plate without an elastic foundation and for the UD pattern reduces by 1.83 times, while that of the plate resting on an elastic foundation reduces by 1.55 times. As can be seen from this table, the shear buckling loads of plates with circular and elliptical cutouts are almost identical. Figure 3, Figure 4, and Figure 5 show the first four buckling mode shapes of the full-clamped FG-GPL RC plate discussed in previous numerical tables. Mode shape results indicate that the elastic foundation with higher stiffness generates more slender and local buckling waves.
TABLE 9
| a/b | Load case | Kw (N/m3) | UD | FG-X | FG-O | FG-V | |
|---|---|---|---|---|---|---|---|
| 1 | Uniaxial | 0 | Mode 1 | 13.433 | 16.394 | 8.926 | 11.414 |
| Mode 2 | 14.379 | 17.120 | 9.700 | 12.327 | |||
| Mode 3 | 14.650 | 17.541 | 9.853 | 12.522 | |||
| Mode 4 | 18.809 | 21.312 | 12.993 | 16.352 | |||
| Biaxial | 0 | Mode 1 | 13.301 | 17.324 | 8.566 | 11.164 | |
| Mode 2 | 15.762 | 19.632 | 10.417 | 13.403 | |||
| Mode 3 | 15.762 | 19.632 | 10.417 | 13.403 | |||
| Mode 4 | 18.760 | 22.766 | 12.531 | 16.037 | |||
| Uniaxial | 5.00E+08 | Mode 1 | 18.446 | 20.832 | 13.663 | 16.504 | |
| Mode 2 | 18.622 | 20.980 | 13.819 | 16.679 | |||
| Mode 3 | 20.249 | 22.427 | 15.160 | 18.215 | |||
| Mode 4 | 21.706 | 23.701 | 15.985 | 19.426 | |||
| Biaxial | 5.00E+08 | Mode 1 | 25.996 | 29.847 | 19.786 | 23.486 | |
| Mode 2 | 26.829 | 30.724 | 20.072 | 24.124 | |||
| Mode 3 | 28.139 | 31.840 | 21.065 | 25.474 | |||
| Mode 4 | 28.139 | 31.840 | 21.065 | 25.474 | |||
| 2 | Uniaxial | 0 | Mode 1 | 26.083 | 32.723 | 17.075 | 21.900 |
| Mode 2 | 26.417 | 32.997 | 17.337 | 22.223 | |||
| Mode 3 | 34.562 | 41.668 | 23.020 | 29.355 | |||
| Mode 4 | 35.932 | 42.732 | 24.118 | 30.670 | |||
| Biaxial | 0 | Mode 1 | 5.468 | 7.230 | 3.497 | 4.524 | |
| Mode 2 | 6.304 | 8.191 | 4.072 | 5.236 | |||
| Mode 3 | 9.804 | 12.376 | 6.414 | 8.244 | |||
| Mode 4 | 10.911 | 13.725 | 7.140 | 9.169 | |||
| Uniaxial | 5.00E+08 | Mode 1 | 42.357 | 48.324 | 31.756 | 37.977 | |
| Mode 2 | 42.456 | 48.406 | 31.767 | 38.048 | |||
| Mode 3 | 49.374 | 54.496 | 37.003 | 44.390 | |||
| Mode 4 | 49.478 | 54.640 | 37.114 | 44.392 | |||
| Biaxial | 5.00E+08 | Mode 1 | 17.724 | 20.249 | 13.487 | 16.052 | |
| Mode 2 | 17.792 | 20.315 | 13.666 | 16.103 | |||
| Mode 3 | 17.900 | 20.547 | 13.749 | 16.205 | |||
| Mode 4 | 18.038 | 20.549 | 13.804 | 16.283 |
Influence of uniaxial and biaxial normal loads on the buckling load of the FG-GPL RC plate with a circular cutout resting on an elastic foundation (; , CCCC, and 1.0 wt%).
TABLE 10
| c/b | d/b | UD | FG-X | FG-O | FG-V | ||
|---|---|---|---|---|---|---|---|
| 0.2 | 0.2 | 0 | Mode 1 | 13.433 | 16.394 | 8.926 | 11.414 |
| Mode 2 | 14.379 | 17.120 | 9.700 | 12.327 | |||
| Mode 3 | 14.650 | 17.541 | 9.853 | 12.522 | |||
| Mode 4 | 18.809 | 21.312 | 12.993 | 16.352 | |||
| 0.3 | 0.133 | 0 | Mode 1 | 9.582 | 12.129 | 6.266 | 8.047 |
| Mode 2 | 10.156 | 12.626 | 6.715 | 8.585 | |||
| Mode 3 | 11.292 | 14.002 | 7.453 | 9.569 | |||
| Mode 4 | 13.293 | 15.966 | 8.930 | 11.380 | |||
| 0.133 | 0.3 | 0 | Mode 1 | 22.885 | 24.207 | 15.983 | 20.012 |
| Mode 2 | 23.261 | 24.402 | 16.368 | 20.372 | |||
| Mode 3 | 23.270 | 24.646 | 16.818 | 20.453 | |||
| Mode 4 | 26.699 | 27.676 | 19.704 | 23.540 | |||
| 0.2 | 0.2 | 5.00E+08 | Mode 1 | 18.446 | 20.832 | 13.663 | 16.504 |
| Mode 2 | 18.622 | 20.980 | 13.819 | 16.679 | |||
| Mode 3 | 20.249 | 22.427 | 15.160 | 18.215 | |||
| Mode 4 | 21.706 | 23.701 | 15.985 | 19.426 | |||
| 0.3 | 0.133 | 5.00E+08 | Mode 1 | 16.531 | 18.973 | 12.355 | 14.849 |
| Mode 2 | 16.551 | 18.976 | 12.459 | 14.879 | |||
| Mode 3 | 16.894 | 19.282 | 12.730 | 15.219 | |||
| Mode 4 | 17.207 | 19.648 | 12.764 | 15.407 |
Buckling loads of the FG-GPL RC plate with circular/elliptical cutouts under uniaxial load resting on an elastic foundation (; , CCCC, and 1.0 wt%).
TABLE 11
| a/b | c/b | d/b | UD | FG-X | FG-O | FG-V | ||
|---|---|---|---|---|---|---|---|---|
| 1 | 0.2 | 0.2 | 0 | Mode 1 | 32.959 | 39.384 | 22.042 | 28.184 |
| Mode 2 | 34.870 | 40.818 | 23.572 | 29.972 | ||||
| Mode 3 | 50.748 | 54.719 | 35.527 | 44.495 | ||||
| Mode 4 | 53.578 | 55.548 | 38.094 | 47.328 | ||||
| 5.00E+08 | Mode 1 | 46.028 | 48.628 | 34.709 | 41.620 | |||
| Mode 2 | 46.100 | 48.656 | 34.755 | 41.679 | ||||
| Mode 3 | 57.427 | 58.683 | 42.939 | 51.813 | ||||
| Mode 4 | 58.180 | 58.755 | 43.472 | 52.496 | ||||
| 0.3 | 0.133 | 0 | Mode 1 | 32.833 | 38.978 | 22.056 | 28.144 | |
| Mode 2 | 33.777 | 39.595 | 22.858 | 29.055 | ||||
| Mode 3 | 49.547 | 53.796 | 34.936 | 43.676 | ||||
| Mode 4 | 50.335 | 54.099 | 35.973 | 44.682 | ||||
| 5.00E+08 | Mode 1 | 45.735 | 48.080 | 34.480 | 41.375 | |||
| Mode 2 | 45.761 | 48.088 | 34.519 | 41.410 | ||||
| Mode 3 | 54.999 | 60.474 | 41.055 | 49.639 | ||||
| Mode 4 | 55.216 | 60.527 | 41.233 | 49.859 | ||||
| 2 | 0.2 | 0.2 | 0 | Mode 1 | 18.035 | 22.545 | 11.837 | 15.168 |
| Mode 2 | 19.434 | 24.030 | 12.828 | 16.382 | ||||
| Mode 3 | 20.871 | 25.600 | 13.840 | 17.635 | ||||
| Mode 4 | 24.820 | 30.118 | 16.455 | 20.959 | ||||
| 5.00E+08 | Mode 1 | 29.613 | 33.697 | 22.382 | 26.634 | |||
| Mode 2 | 29.914 | 34.026 | 22.499 | 26.853 | ||||
| Mode 3 | 31.346 | 35.463 | 23.484 | 28.026 | ||||
| Mode 4 | 32.069 | 36.132 | 23.830 | 28.579 | ||||
| 0.3 | 0.133 | 0 | Mode 1 | 16.976 | 21.185 | 11.166 | 14.275 | |
| Mode 2 | 18.002 | 22.235 | 11.905 | 15.180 | ||||
| Mode 3 | 19.536 | 24.162 | 12.898 | 16.470 | ||||
| Mode 4 | 23.134 | 28.387 | 15.292 | 19.504 | ||||
| 5.00E+08 | Mode 1 | 28.504 | 32.388 | 21.687 | 25.661 | |||
| Mode 2 | 28.659 | 32.576 | 21.752 | 25.773 | ||||
| Mode 3 | 30.581 | 34.632 | 22.695 | 27.309 | ||||
| Mode 4 | 31.174 | 35.226 | 22.932 | 27.725 |
Shear buckling loads of the FG-GPL RC plate with circular/elliptical cutouts resting on an elastic foundation (, CCCC, and 1.0 wt%).
FIGURE 3
FIGURE 4
FIGURE 5
4 Conclusion
A set of numerical studies has been conducted to describe the buckling behavior of a functionally graded multilayer GPL-reinforced composite plate with circular/elliptical cutouts resting on a Winkler-type elastic foundation under uniaxial and biaxial normal and shear loads. Formulations are proposed based on the third-order shear deformation plate theory and finite element procedure. A broad range of factors such as plate aspect ratio, plate width/thickness ratio, size of the circular cutout, applying uniaxial and biaxial normal and shear loads, several Winkler elastic foundation stiffness parameters, different displacement boundary conditions, circular and elliptical cutouts, orientation of the elliptical cutout, and influence of the GPL weight fraction have been investigated for four GPL distribution patterns to analyze the buckling behavior of the FG-GPL RC plates. The main conclusions of the present study are summarized as follows:
➢ In general, the stiffness of the outer layers has the most remarkable effect on the stiffness of the composite structure. Therefore, as described previously, the structure with the FG-X gradient pattern presented the highest stiffness and buckling capabilities, and UD, FG-V, and FG-O patterns in the sequence have lower values.
➢ By dispersing a low content of GPL nanofillers to the matrix, the buckling loads can remarkably improve. For example, by dispersion of only 0.5 wt% and 1.0 wt% GPLs in the epoxy matrix, the FG-X plate buckling load enhanced 2.97 and 4.84 times, respectively.
➢ By increasing the cutout radius-to-width ratio ( ) from to 0.3, the buckling loads decreased by 22.3%.
➢ Buckling loads increased by 42.6% for the FG-GPL RC plate when the elastic foundation stiffness increased from 0 to 1 .
➢ If the major axis of the elliptical cutout is parallel to the direction of the applied load, the critical buckling load for the plate with the UD pattern is approximately 2.39 times lower than that when the semi-major radius of the elliptical cutout is perpendicular to the applied load.
➢ The influence of the elastic foundation on the plate with an elliptical cutout is more remarkable compared to the plate with a circular cutout. When the plate with an elliptical cutout is embedded in an elastic foundation, the critical buckling load of the plate with the UD pattern is approximately 1.73 times greater compared with that of the plate without an elastic foundation, while for the plate with a circular cutout, it is approximately 1.37 times higher.
➢ By increasing the plate aspect ratio from 1 to 2, the critical shear buckling load of the plate without an elastic foundation for the UD pattern reduces by 1.83 times, while that of the plate resting on an elastic foundation reduces by 1.55 times.
➢ The shear buckling loads of the plates with circular and elliptical cutouts with an equal area are almost identical.
Statements
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Author contributions
AK: formal analysis, investigation, software, validation, writing–original draft, and writing–review and editing. MB: data curation, formal analysis, investigation, methodology, software, validation, writing–original draft, and writing–review and editing. KA: conceptualization, data curation, formal analysis, investigation, methodology, supervision, validation, writing–original draft, and writing–review and editing.
Funding
The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
Abedini BaghbadoraniA.KianiY. (2021). Free vibration analysis of functionally graded cylindrical shells reinforced with graphene platelets. Compos. Struct.276, 114546. 10.1016/j.compstruct.2021.114546
2
AlashkarA.ElkafrawyM.HawilehR.AlHamaydehM. (2022). Buckling analysis of functionally graded materials (FGM) thin plates with various circular cutout arrangements. J. Compos. Sci.6 (9), 277. 10.3390/jcs6090277
3
AnsariR.TorabiJ.HassaniR. (2018). In-plane and shear buckling analysis of FG-CNTRC annular sector plates based on the third-order shear deformation theory using a numerical approach. Comput. &Amp; Math. Appl.75 (2), 486–502. 10.1016/j.camwa.2017.09.022
4
AsemiK.AshrafiH.ShariyatM. (2016a). Three-dimensional stress and free vibration analyses of functionally graded plates with circular holes by the use of the graded finite element method. J. Appl. Mech. Tech. Phys.57 (4), 690–700. 10.1134/s0021894416040131
5
AsemiK.BabaeiM.KiarasiF. (2020). Static, natural frequency and dynamic analyses of functionally graded porous annular sector plates reinforced by graphene platelets. Mech. Based Des. Struct. Mach.50 (11), 3853–3881. 10.1080/15397734.2020.1822865
6
AsemiK.SalehiM.AkhlaghiM. (2014). Three dimensional biaxial buckling analysis of functionally graded annular sector plate fully or partially supported on Winkler elastic foundation. Aerosp. Sci. Technol.39, 426–441. 10.1016/j.ast.2014.04.011
7
AsemiK.SalehiM.AkhlaghiM. (2015). Three dimensional graded finite element elasticity shear buckling analysis of FGM annular sector plates. Aerosp. Sci. Technol.43, 1–13. 10.1016/j.ast.2015.02.009
8
AsemiK.ShariyatM.SalehiM.AshrafiH. (2013). A full compatible three-dimensional elasticity element for buckling analysis of FGM rectangular plates subjected to various combinations of biaxial normal and shear loads. Finite Elem. Analysis Des.74, 9–21. 10.1016/j.finel.2013.05.011
9
AshrafiH.AsemiK.ShariyatM. (2013). A three-dimensional boundary element stress and bending analysis of transversely/longitudinally graded plates with circular cutouts under biaxial loading. Eur. J. Mech. - a/Solids42, 344–357. 10.1016/j.euromechsol.2013.07.009
10
BabaeiM.AsemiK. (2020). Static, dynamic and natural frequency analyses of functionally graded carbon nanotube annular sector plates resting on viscoelastic foundation. SN Appl. Sci.2 (10), 1652. 10.1007/s42452-020-03421-7
11
BuiT. Q.DoT. V.TonL. H. T.DoanD. H.TanakaS.PhamD. T.et al (2016). On the high temperature mechanical behaviors analysis of heated functionally graded plates using FEM and a new third-order shear deformation plate theory. Compos. Part B Eng.92, 218–241. 10.1016/j.compositesb.2016.02.048
12
CongP. H.DucN. D. (2018). New approach to investigate the nonlinear dynamic response and vibration of a functionally graded multilayer graphene nanocomposite plate on a viscoelastic Pasternak medium in a thermal environment. Acta Mech.229 (9), 3651–3670. 10.1007/s00707-018-2178-3
13
EsmaeiliH. R.KianiY.BeniY. T. (2022). Vibration characteristics of composite doubly curved shells reinforced with graphene platelets with arbitrary edge supports. Acta Mech.233 (2), 665–683. 10.1007/s00707-021-03140-z
14
FanF.LeiB.SahmaniS.SafaeiB. (2020). On the surface elastic-based shear buckling characteristics of functionally graded composite skew nanoplates. Thin-Walled Struct.154, 106841. 10.1016/j.tws.2020.106841
15
FarsadiT.RahmanianM.KurtaranH. (2022). Nonlinear stability of multilayered graphene platelet-reinforced functionally graded wing-like plates. Acta Mech.233 (8), 3233–3252. 10.1007/s00707-022-03238-y
16
FattahiA. M.SafaeiB.AhmedN. A. (2019). A comparison for the non-classical plate model based on axial buckling of single-layered graphene sheets. Eur. Phys. J. Plus134 (11), 555. 10.1140/epjp/i2019-12912-7
17
García-MacíasE.Rodríguez-TemblequeL.SáezA. (2018). Bending and free vibration analysis of functionally graded graphene vs. carbon nanotube reinforced composite plates. Compos. Struct.186, 123–138. 10.1016/j.compstruct.2017.11.076
18
GengX.ZhaoL.ZhouW. (2020). Finite-element buckling analysis of functionally graded GPL-reinforced composite plates with a circular hole. Mech. Based Des. Struct. Mach.49 (7), 1028–1044. 10.1080/15397734.2019.1707688
19
HungD. X.TuT. M.HaoT. D. (2021). Free vibration analysis of laminated CNTRC plates using the pb2-ritz method. J. Mech. Eng.18 (1), 213–232. 10.24191/jmeche.v18i1.15181
20
JafariP.KianiY. (2021). Free vibration of functionally graded graphene platelet reinforced plates: a quasi 3D shear and normal deformable plate model. Compos. Struct.275, 114409. 10.1016/j.compstruct.2021.114409
21
Lakshmi narayanaA.Vijaya KumarR.Krishnamohana RaoG. (2018). Thermal buckling analysis of laminated composite plate with square/rectangular, elliptical/circular cutout. Mater. Today Proc.5 (2), 5354–5363. 10.1016/j.matpr.2017.12.121
22
LeiZ.ZhangL.LiewK. (2015). Free vibration analysis of laminated FG-CNT reinforced composite rectangular plates using the kp-Ritz method. Compos. Struct.127, 245–259. 10.1016/j.compstruct.2015.03.019
23
LinF.XiangY.ShenH. S. (2017). Temperature dependent mechanical properties of graphene reinforced polymer nanocomposites – a molecular dynamics simulation. Compos. Part B Eng.111, 261–269. 10.1016/j.compositesb.2016.12.004
24
LinH. G.CaoD. Q.XuY. Q. (2018). Vibration, buckling and aeroelastic analyses of functionally graded multilayer graphene-nanoplatelets-reinforced composite plates embedded in piezoelectric layers. Int. J. Appl. Mech.10 (03), 1850023. 10.1142/s1758825118500230
25
MaR.JinQ. (2022). Stability of functionally graded graphene-reinforced composite laminated thick plates in thermal environment. Acta Mech.233 (10), 3977–3996. 10.1007/s00707-022-03300-9
26
MohammadiM.FarajpourA.MoradiA.GhayourM. (2014). Shear buckling of orthotropic rectangular graphene sheet embedded in an elastic medium in thermal environment. Compos. Part B Eng.56, 629–637. 10.1016/j.compositesb.2013.08.060
27
MollaeiS.BabaeiM.AsemiK. (2023). Torsional buckling of functionally graded graphene reinforced composite laminated cylindrical panel. Archive Appl. Mech.93 (2), 427–435. 10.1007/s00419-022-02132-2
28
Muni Rami ReddyR.KarunasenaW.LokugeW. (2018). Free vibration of functionally graded-GPL reinforced composite plates with different boundary conditions. Aerosp. Sci. Technol.78, 147–156. 10.1016/j.ast.2018.04.019
29
NguyenH. N.TanT. C.LuatD. T.PhanV. D.ThomD. V.MinhP. V. (2019). Research on the buckling behavior of functionally graded plates with stiffeners based on the third-order shear deformation theory. Materials12 (8), 1262. 10.3390/ma12081262
30
PradhanS. C.PhadikarJ. K. (2010). Scale effect and buckling analysis of multilayered graphene sheets based on nonlocal continuum Mechanics. J. Comput. Theor. Nanosci.7 (10), 1948–1954. 10.1166/jctn.2010.1565
31
SaiahB.BacheneM.GuemanaM.ChikerY.AttafB. (2022). On the free vibration behavior of nanocomposite laminated plates contained piece-wise functionally graded graphene-reinforced composite plies. Eng. Struct.253, 113784. 10.1016/j.engstruct.2021.113784
32
ShahrjerdiA.YavariS. (2018). Free vibration analysis of functionally graded graphene-reinforced nanocomposite beams with temperature-dependent properties. J. Braz. Soc. Mech. Sci. Eng.40 (1), 25. 10.1007/s40430-017-0943-1
33
ShariyatM.AsemiK. (2014b). Three-dimensional non-linear elasticity-based 3D cubic B-spline finite element shear buckling analysis of rectangular orthotropic FGM plates surrounded by elastic foundations. Compos. Part B Eng.56, 934–947. 10.1016/j.compositesb.2013.09.027
34
ShenH. S.XiangY.FanY. (2018). Postbuckling of functionally graded graphene-reinforced composite laminated cylindrical panels under axial compression in thermal environments. Int. J. Mech. Sci.135, 398–409. 10.1016/j.ijmecsci.2017.11.031
35
ShenH. S.XiangY.LinF. (2017). Thermal buckling and postbuckling of functionally graded graphene-reinforced composite laminated plates resting on elastic foundations. Thin-Walled Struct.118, 229–237. 10.1016/j.tws.2017.05.006
36
SongM.YangJ.KitipornchaiS. (2018). Bending and buckling analyses of functionally graded polymer composite plates reinforced with graphene nanoplatelets. Compos. Part B Eng.134, 106–113. 10.1016/j.compositesb.2017.09.043
37
SongM.YangJ.KitipornchaiS.ZhuW. (2017). Buckling and postbuckling of biaxially compressed functionally graded multilayer graphene nanoplatelet-reinforced polymer composite plates. Int. J. Mech. Sci.131 (132), 345–355. 10.1016/j.ijmecsci.2017.07.017
38
ThaiC. H.FerreiraA.TranT.Phung-VanP. (2019). Free vibration, buckling and bending analyses of multilayer functionally graded graphene nanoplatelets reinforced composite plates using the NURBS formulation. Compos. Struct.220, 749–759. 10.1016/j.compstruct.2019.03.100
39
YangJ.ChenD.KitipornchaiS. (2018). Buckling and free vibration analyses of functionally graded graphene reinforced porous nanocomposite plates based on Chebyshev-Ritz method. Compos. Struct.193, 281–294. 10.1016/j.compstruct.2018.03.090
40
YangJ.WuH.KitipornchaiS. (2017a). Buckling and postbuckling of functionally graded multilayer graphene platelet-reinforced composite beams. Compos. Struct.161, 111–118. 10.1016/j.compstruct.2016.11.048
41
ZhangL.LeiZ.LiewK. (2015). Free vibration analysis of functionally graded carbon nanotube-reinforced composite triangular plates using the FSDT and element-free IMLS-Ritz method. Compos. Struct.120, 189–199. 10.1016/j.compstruct.2014.10.009
Appendix
A four-noded Lagrangian quadrangular element is used. The associated linear Lagrangian interpolation functions in terms of the natural coordinates are as follows:
In addition, the Hermitian interpolation functions in terms of the natural coordinates are as follows:
Summary
Keywords
graphene, buckling, multilayer plate, third-order shear deformation plate theory, Winkler
Citation
Kalhori A, Bayat MJ and Asemi K (2023) Buckling response of functionally graded multilayer graphene platelet-reinforced composite plates with circular/elliptical cutouts supporting on an elastic foundation under normal and shear loads. Front. Mech. Eng 9:1293713. doi: 10.3389/fmech.2023.1293713
Received
13 September 2023
Accepted
16 October 2023
Published
26 October 2023
Volume
9 - 2023
Edited by
Wei Zhao, Oklahoma State University, United States
Reviewed by
Pranav Borwankar, Virginia Tech, United States
Riccardo Vescovini, Polytechnic University of Milan, Italy
Updates
Copyright
© 2023 Kalhori, Bayat and Asemi.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Kamran Asemi, k.asemi@iau-tnb.ac.ir
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.