ORIGINAL RESEARCH article

Front. Mech. Eng., 26 October 2023

Sec. Solid and Structural Mechanics

Volume 9 - 2023 | https://doi.org/10.3389/fmech.2023.1293713

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

  • Department of Mechanical Engineering, Islamic Azad University, North Tehran Branch, Tehran, Iran

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 nameEpoxy matrixGPL
Modulus of elasticity 31,010
Poisson’s ratio0.340.186
Density 1,2001,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 patternPresentDifference %
UD4.32364.28470.9077
FG-X5.55445.54390.1893
FG-O3.04762.95403.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 36 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/hUDFG-XFG-OFG-V
10Mode 1426.790555.260269.750348.840
Mode 2604.830710.790401.040503.030
Mode 3747.950805.940522.560633.980
Mode 4923.590999.980627.820775.960
15Mode 129.94433.86120.67126.038
Mode 231.55634.68222.23627.738
Mode 332.17735.56122.58928.151
Mode 439.84240.56629.25735.412
20Mode 113.43316.3948.92611.414
Mode 214.37917.1209.70012.327
Mode 314.65017.5419.85312.522
Mode 418.80921.31212.99316.352
30Mode 14.1565.3842.6863.458
Mode 24.5165.7552.9453.779
Mode 34.5895.8742.9863.835
Mode 46.0417.5173.9805.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%)UDFG-XFG-OFG-V
0.0Mode 13.3903.3903.3903.390
Mode 23.6293.6293.6293.629
Mode 33.6983.6983.6983.698
Mode 44.7474.7474.7474.747
0.5Mode 18.40810.0816.1927.635
Mode 29.00010.6036.7098.221
Mode 39.17010.8476.8188.358
Mode 411.77313.4158.95010.856
1.0Mode 113.43316.3948.92611.414
Mode 214.37917.1209.70012.327
Mode 314.65017.5419.85312.522
Mode 418.80921.31212.99316.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

BCUDFG-XFG-OFG-V
SSSSMode 112.80715.6308.45210.666
Mode 213.74016.3469.21511.564
Mode 313.89616.6379.27911.623
Mode 417.85420.25912.23415.182
CCCCMode 113.43316.3948.92611.414
Mode 214.37917.1209.70012.327
Mode 314.65017.5419.85312.522
Mode 418.80921.31212.99316.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/bUDFG-XFG-OFG-V
1Mode 113.43316.3948.92611.414
Mode 214.37917.1209.70012.327
Mode 314.65017.5419.85312.522
Mode 418.80921.31212.99316.352
2Mode 126.08332.72317.07521.900
Mode 226.41732.99717.33722.223
Mode 334.56241.66823.02029.355
Mode 435.93242.73224.11830.670
4Mode 124.05630.36115.70720.163
Mode 224.05730.36115.70720.163
Mode 330.85437.64520.44626.086
Mode 430.85437.64620.44726.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 sizeUDFG-XFG-OFG-V
×0.5 (r/b = 0.1)Mode 115.17019.2939.87312.715
Mode 217.61021.66711.66414.881
Mode 326.88130.67518.46223.207
Mode 430.53336.49120.42826.123
×1 (r/b = 0.2)Mode 113.43316.3948.92611.414
Mode 214.37917.1209.70012.327
Mode 314.65017.5419.85312.522
Mode 418.80921.31212.99316.352
×1.5 (r/b = 0.3)Mode 111.87113.4088.25210.326
Mode 212.02713.5568.35510.488
Mode 312.11213.5718.41310.573
Mode 412.12213.6788.47110.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)UDFG-XFG-OFG-V
0Mode 113.43316.3948.92611.414
Mode 214.37917.1209.70012.327
Mode 314.65017.5419.85312.522
Mode 418.80921.31212.99316.352
1.00E+08Mode 114.98317.70410.48613.033
Mode 215.42018.05910.76113.415
Mode 316.15718.86011.40014.095
Mode 419.47721.86113.71417.074
5.00E+08Mode 118.44620.83213.66316.504
Mode 218.62220.98013.81916.679
Mode 320.24922.42715.16018.215
Mode 421.70623.70115.98519.426
1.00E+09Mode 121.31323.38516.34419.400
Mode 221.54423.58616.53319.628
Mode 323.20324.94817.73921.130
Mode 423.87025.47518.10321.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/bLoad caseKw (N/m3)UDFG-XFG-OFG-V
1Uniaxial0Mode 113.43316.3948.92611.414
Mode 214.37917.1209.70012.327
Mode 314.65017.5419.85312.522
Mode 418.80921.31212.99316.352
Biaxial0Mode 113.30117.3248.56611.164
Mode 215.76219.63210.41713.403
Mode 315.76219.63210.41713.403
Mode 418.76022.76612.53116.037
Uniaxial5.00E+08Mode 118.44620.83213.66316.504
Mode 218.62220.98013.81916.679
Mode 320.24922.42715.16018.215
Mode 421.70623.70115.98519.426
Biaxial5.00E+08Mode 125.99629.84719.78623.486
Mode 226.82930.72420.07224.124
Mode 328.13931.84021.06525.474
Mode 428.13931.84021.06525.474
2Uniaxial0Mode 126.08332.72317.07521.900
Mode 226.41732.99717.33722.223
Mode 334.56241.66823.02029.355
Mode 435.93242.73224.11830.670
Biaxial0Mode 15.4687.2303.4974.524
Mode 26.3048.1914.0725.236
Mode 39.80412.3766.4148.244
Mode 410.91113.7257.1409.169
Uniaxial5.00E+08Mode 142.35748.32431.75637.977
Mode 242.45648.40631.76738.048
Mode 349.37454.49637.00344.390
Mode 449.47854.64037.11444.392
Biaxial5.00E+08Mode 117.72420.24913.48716.052
Mode 217.79220.31513.66616.103
Mode 317.90020.54713.74916.205
Mode 418.03820.54913.80416.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/bd/bUDFG-XFG-OFG-V
0.20.20Mode 113.43316.3948.92611.414
Mode 214.37917.1209.70012.327
Mode 314.65017.5419.85312.522
Mode 418.80921.31212.99316.352
0.30.1330Mode 19.58212.1296.2668.047
Mode 210.15612.6266.7158.585
Mode 311.29214.0027.4539.569
Mode 413.29315.9668.93011.380
0.1330.30Mode 122.88524.20715.98320.012
Mode 223.26124.40216.36820.372
Mode 323.27024.64616.81820.453
Mode 426.69927.67619.70423.540
0.20.25.00E+08Mode 118.44620.83213.66316.504
Mode 218.62220.98013.81916.679
Mode 320.24922.42715.16018.215
Mode 421.70623.70115.98519.426
0.30.1335.00E+08Mode 116.53118.97312.35514.849
Mode 216.55118.97612.45914.879
Mode 316.89419.28212.73015.219
Mode 417.20719.64812.76415.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/bc/bd/bUDFG-XFG-OFG-V
10.20.20Mode 132.95939.38422.04228.184
Mode 234.87040.81823.57229.972
Mode 350.74854.71935.52744.495
Mode 453.57855.54838.09447.328
5.00E+08Mode 146.02848.62834.70941.620
Mode 246.10048.65634.75541.679
Mode 357.42758.68342.93951.813
Mode 458.18058.75543.47252.496
0.30.1330Mode 132.83338.97822.05628.144
Mode 233.77739.59522.85829.055
Mode 349.54753.79634.93643.676
Mode 450.33554.09935.97344.682
5.00E+08Mode 145.73548.08034.48041.375
Mode 245.76148.08834.51941.410
Mode 354.99960.47441.05549.639
Mode 455.21660.52741.23349.859
20.20.20Mode 118.03522.54511.83715.168
Mode 219.43424.03012.82816.382
Mode 320.87125.60013.84017.635
Mode 424.82030.11816.45520.959
5.00E+08Mode 129.61333.69722.38226.634
Mode 229.91434.02622.49926.853
Mode 331.34635.46323.48428.026
Mode 432.06936.13223.83028.579
0.30.1330Mode 116.97621.18511.16614.275
Mode 218.00222.23511.90515.180
Mode 319.53624.16212.89816.470
Mode 423.13428.38715.29219.504
5.00E+08Mode 128.50432.38821.68725.661
Mode 228.65932.57621.75225.773
Mode 330.58134.63222.69527.309
Mode 431.17435.22622.93227.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), 486502. 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), 690700. 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), 38533881. 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, 426441. 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, 113. 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, 921. 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, 344357. 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, 218241. 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), 36513670. 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), 665683. 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), 32333252. 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, 123138. 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), 10281044. 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), 213232. 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), 53545363. 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, 245259. 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, 261269. 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), 39773996. 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, 629637. 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), 427435. 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, 147156. 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), 19481954. 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, 934947. 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, 398409. 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, 229237. 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, 106113. 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), 345355. 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, 749759. 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, 281294. 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, 111118. 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, 189199. 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

*Correspondence: Kamran Asemi,

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics