Size-Dependent Buckling and Post-Buckling Analysis of the Functionally Graded Thin Plate Al–Cu Material Based on a Modified Couple Stress Theory

Size-dependent functionally graded material thin plate buckling and post-buckling problems are considered using the framework of the MCST (Modified Couple Stress Theory). Based on modified couple stress theory and power law, the post-buckling deflection and critical buckling load of simply supported functionally graded material thin plate are derived using Hamilton’s minimum potential energy principle. The analysis compares the simulation results of linear buckling and nonlinear buckling. Innovatively, a power-law distribution with scale effects is considered. The influences of scale effect parameters l and power-law index parameters k on buckling displacement, load, and strain energy of plates have been investigated. In this article, it is found that the critical buckling displacement, critical buckling load, and buckling strain energy increase with increases in the power-law index parameters k. The membrane energy decreases as the power-law index parameter increases. If the upper and lower layers are swapped, the opposite result is obtained. In comparison, the scale effect parameter is more influential than the power-law exponent. The critical buckling displacement in the x-direction is not affected by scale effects. The critical buckling load, the membrane energy, and buckling strain energy increase as the scale effect parameter increases. Scale effects increase material stiffness compared with traditional theory, and the power-law index parameters affect FGM properties such as elastic modulus, Poisson’s ratio, density, etc. Both scale effects parameters and power-law index parameters have important effects on the mechanical behavior of materials.


Introduction
MEMS/NEMS (micro-electro-mechanical systems/nano-electro-mechanical systems) devices play a very significant role in high-precision fields, showing unprecedented application prospects in aerospace, consumer electronics, medical equipment, and other fields [1][2][3]. The main structural units of MEMS/NEMS micro-sensor devices are microstructures such as micro-beams, micro-plates, and micro-shells. The size of the structure is generally in the micrometer or even nanometer scale. A large number of experiments [4][5][6] have shown that when the structural size of the material reaches the micro-nano level, the mechanical properties of the material are greatly different from those of the macro, and the classical theory will no longer be applicable. Functionally graded material is a new type of composite material composed of two or more materials with one-dimensional or multi-dimensional continuous component distribution. The combination of multiple components can achieve excellent performance [7][8][9] such as higher thermal resistance, higher fracture toughness, better corrosion resistance, etc. The application of functionally graded materials in MEMS/NEMS devices is expected to improve device performance. Compared with traditional composite materials, functionally graded materials have their unique advantages. Due to mutation in the material interface of traditional composite materials, the deformation difference caused by the different thermal expansion coefficients of the two materials in the thermal environment eventually generates thermal stress concentration on the device, which may bring failure risks such as buckling, post-buckling, delamination, and even fracture. Compared with traditional composite materials, functionally graded materials have no obvious interface, and their composition and structure show continuous gradient changes [10]. The continuous material distribution can greatly improve the bond strength [11], reduce the crack driving force [12], reduce the residual stress and thermal stress [13], and reduce the risk of stress concentration [14]. The schematic diagram of the cross-section of the functionally graded material distributed along the thickness direction is shown in Figure 1. Functionally graded materials contain micro-porous defects, and porosity is one of the most important parameters affecting the mechanical behavior of FGM (functionally graded materials). Souhir Zghal et al. [15] concluded that the porosity leads to a decrease in the elastic modulus, so the deflection of the FGMs beams increases. The coefficient and distribution of porosity directly affect the mechanical properties of the FGMs. Gao et al. [16] concluded that increased porosity reduces the critical buckling load, reducing the stiffness of the beam. The lowest critical buckling load is obtained with the porosity in the T-3 distribution (symmetric with weaker top and bottom surfaces). Static and dynamic responses of functionally graded micro/nano-plates are the basis for the development of functionally graded MEMS devices. Jung et al. [17] studied the bending and vibration of a FGM micro-plate on a Pasternak elastic foundation by coupling the first-order shear deformation theory and the modified couple stress theory. The changes in the bending deflection and natural frequency of the FGM micro-plate with the foundation parameters, material gradient parameters, and scale effect parameters under uniform load are analyzed. Eshraghi et al. [18] studied the free vibration response of the FGM micro-rectangular plate based on the higher-order shear deformation theory and the modified couple stress theory and obtained the analytical solution of the natural frequency. Song [19] et al. studied the vibration response of the FGM micro-beam in the medium under photothermal excitation, obtained the resonance frequency in different fluid environments, and determined the effect of the gradient factor on the resonance frequency. Recently, functionally graded (FG) grapheme platelets (GPLs) and reinforced composites (FG-GPLRC) have great application potential in MEMS/NEMS devices due to their superior mechanical properties such as high strength, good thermal conductivity, exceptionally high elastic modulus, and large specific surface area [20]. Yang's team [21][22][23] was the first to propose a multilayer functionally graded graphene platelets reinforced composite material with a layered variation along the thickness direction, and conducted pioneering work on its mechanical properties. The results show that the mechanical properties of FG-GPLRC are significantly better than those of GPLs. Furthermore, incorporating graphene origami into FGMs can also yield unique mechanical properties such as a negative Poisson's ratio [24], resulting in versatility and superior structural performance. Yang et al. [24] propose an efficient micro-mechanical model based on a genetic programming (GP) algorithm and molecular dynamics (MD) simulation to predict the mechanical properties of graphene origami metal functionally graded materials. Theoretical results show that the proposed functionally graded metamaterial beam achieves significantly improved bending properties.
Compared with macro-devices, micro-devices represent not only a reduction in size but also a substantial change in mechanical properties. The traditional theory is no longer applicable. The traditional couple stress theory [25][26][27] considers the displacement and rotation of the point particle at the same time, but there are many material parameters considered in the scale effect, making it unsuitable for practical engineering applications. In 2002, Yang et al. [28] proposed a modified couple stress theory that only requires one scale parameter, which reduces the difficulty of solving and also has the potential to guide engineering practice. Chen et al. [29] proposed a new modified couple stress theory for anisotropic materials in 2011. In recent years, Richa Priyanka [30] et al. studied the stability and free vibration characteristics of micro-lamination under varying axial loads based on the MCST (modified couple stress theory). The study found that beams with high buckling loads and high natural frequencies under C-C (Clamped-Clamped) support conditions have more significant scale effects. J. N. Reddy et al. [31] verified the applicability of the modified coupled stress theory for web-core sandwich panels. Li et al. [32] designed a standard experimental method to determine the scale effect parameters of the modified couple stress theory. The scale effect parameters were calculated by measuring the vibration frequency of the micro-cantilever beam. The couple stress theory is an effective method to describe the micro size scale effect, and its dimension is at the micron level. However, when the structural size of the material reaches the nanometer level, the specific surface area also needs to be considered, and the influence of the surface effect cannot be ignored. Non-local theory [33,34], strain gradient theory [35] and surface elasticity theory [36] can be used to describe the scale effect and surface effect of nanostructures. Song [37] et al. considered the refraction and reflection of thermal, elastic, and plasma waves in semiconductor microstructures based on nonlocal theory. The research indicates that wave properties show a significant difference between nonlocal theory and classical theory. Li et al. [38] studied the propagation properties of flexural waves in FGM beams based on the theory of nonlocal strain gradients and showed that the optical phase velocities and acoustical phase velocities increase with decreasing nonlocal or material length scale parameters. Y.S. Li et al. [39] studied the bending properties of piezoelectric nano-plates based on surface effects, concluding that surface effects enhance plate stiffness depending on the thickness of the plate. These theories describing micro-scale structures are derived from continuum mechanics as a complement to classical mechanics to describe the properties of micro-scales.  Compared with macro-devices, micro-devices represent not only a reduction in size but also a substantial change in mechanical properties. The traditional theory is no longer applicable. The traditional couple stress theory [25][26][27] considers the displacement and rotation of the point particle at the same time, but there are many material parameters considered in the scale effect, making it unsuitable for practical engineering applications. In 2002, Yang et al. [28] proposed a modified couple stress theory that only requires one scale parameter, which reduces the difficulty of solving and also has the potential to guide engineering practice. Chen et al. [29] proposed a new modified couple stress theory for anisotropic materials in 2011. In recent years, Richa Priyanka [30] et al. studied the stabil- The post-buckling problem is a geometric nonlinear problem [40], which is caused by large out-of-plane deflection. Von Karman was the first to derive nonlinear equilibrium equations for large deflections of plates. L. L. Ke et al. [41,42] studied the post-buckling behavior and nonlinear vibration characteristics of piezoelectric nano-beams under thermo-electro-mechanical loadings based on the Von Karman nonlinear equation. Based on the Von on Karman plate theory, Chen et al. [43] considered the influence of geometric nonlinearity and damage effects on the post-buckling characteristics of micro-plates. The study found that the larger the scale effect parameter under the same load condition, the less destructive the effect on the post-buckling properties. Amin Ghorbani Shenas et al. [44] analyzed the thermal buckling, post-buckling, and free vibration properties of rotating pre-twisted FGM beams. The effects of angular velocity, torsion angle, scale effect parameter, and power-law index on the post-buckling equilibrium path and linear free vibration behavior of FGM beams were investigated. Li et al. [45] studied the influence of scale effects on post-buckling behavior of FGM nano-beams with Von Karman geometric nonlinearity based on nonlocal strain gradient theory, which found that the values of scaling parameters determine stiffness-softening or stiffness-hardening.
According to the above, there are many studies on the mechanical behavior of FGM materials, and power-law functions [46][47][48] are often used to describe the material properties of FGM. Density, elastic modulus, Poisson's ratio, etc. are inherent properties of materials. When the size of the material reaches the micro-scale, the scale effects cannot be ignored, and the scale effect parameter is also an inherent property of the micro-scale material based on MCST. Therefore, considering the power-law distribution of scale effects is of great significance. However, there is no mathematical model study on the distribution of the power-law function of the scale effect of FGM. In this paper, the equilibrium method based on MCST combined with the power-law is adopted to obtain the elastic modulus and stiffness of a thin FGM plate. The power-law distribution function of the scale effect parameter is introduced into the equilibrium equation. Scale effect parameters of thin FGM plates vary with the thickness and proportions of materials and are uniquely determined. Coupling the modified couple stress theory into the S. Timoshenko et al. [49][50][51] plate theory, the equilibrium equation of the Al-Cu FGM thin plate is obtained. The critical buckling displacement and critical buckling load of the FGM thin plate are obtained according to Hamilton's principle [52,53]. The changes in bending energy and membrane energy with power-law index k and post-buckling stiffness are analyzed. The stability of MEMS/NEMS in multi-field and multi-scale environments is still a challenge [54][55][56][57][58]. The content of this paper is of great significance for MEMS/NEMS manufacturing and stability analysis.

Power-Law of FGM
Considering that the FGM thin plate is composed of Al and Cu, the material changes continuously from Cu on the lower surface to Al on the upper surface. Power-law functions [46][47][48] are often used to describe the material properties of FGM. The gradient change of material properties can be expressed as a power-law function as [46][47][48]: where P(z) are the material property parameters of the FGM thin plate at z, P U and P L are the material property parameters of the upper and lower surfaces, respectively, k is a power-law index representing the volume fraction, and h is FGM plate thickness. When k = 0, the FGM plate consists entirely of upper surface Al, and when k is infinite, the material consists entirely of lower surface Cu. Therefore, Young's modulus and Poisson's ratio of the FGM plate according to Equation (1) can be expressed as follows: Considering the scale effect parameter l varies with volume fraction, the power-law function of gradient variation introduced into the scale effect parameter l is:

FGM Plate Theory Based on MCST
The constitutive equation of FGM plates can be expressed as follows based on MCST: where σ ij , ε ij , χ ij , m ij are expressed FGM plate stress tensor, strain tensor, symmetric curvature tensor, and couple stress tensor, respectively, θ i is the FGM plate rotation vector and e ijk is the permutation symbol, µ and λ are Lame constants, E is the elastic modulus, υ is the Poisson ratio based on the power-law of FGM, and l is the scale parameter of FGM, which is related to the structure size and power-law index k.
Consider an FGM thin plate of length a, width b, and thickness h under the external lateral load q. The displacement field u i based on Kirchhoff thin plate theory is given as: u(x, y, z) = −z ∂w ∂x (12) v(x, y, z) = −z ∂w ∂y (13) w(x, y, z) = w(x, y, t) The bending moment and shear force according to power-law can be expressed as: The equilibrium equation can be expressed as follows: ∂M y ∂y By using Equations (20)- (22).The FGM thin plate governing equation can be expressed as: where where D(z) is the FGM plate flexural rigidity. Compared with classical governing equations [49][50][51], the flexural rigidity of FGM thin plate based on power-law and MCST follows as:

Post-Buckling of FGM Thin Plate
Consider the geometry and force diagram of the FGM thin plate as shown in Figure 2. The FGM thin plate is of length a, width b, thickness h, and force P. The displacement in the x-direction under the pressure of P is u x ; the displacement in the y-direction is u y , and the deflection is w 0 .
where ( ) D z is the FGM plate flexural rigidity. Compared with classical governing equations [49][50][51], the flexural rigidity of FGM thin plate based on power-law and MCST follows as:

Post-Buckling of FGM Thin Plate
Consider the geometry and force diagram of the FGM thin plate as shown in Figure  2. The FGM thin plate is of length a, width b, thickness h, and force P. The displacement in the x-direction under the pressure of P is ux; the displacement in the y-direction is uy, and the deflection is w0. The expression of displacement function of the FGM thin plate is as follows: Assume that the flexural deflection is a double harmonic function as: The expression of displacement function of the FGM thin plate is as follows: (27) Assume that the flexural deflection is a double harmonic function as: According to the Von Karman nonlinear equation, the strain of a thin plate is expressed as: ε yy = ∂u y ∂x The curvature components can be defined as follows: The total energy at the location (u 0 , w 0 ) is: where U b and U m are strain and membrane energies, respectively. According to references [49][50][51], the FGM thin plate strain energy can be seen as: By using Equations (32)- (34) in Equation (36), the FGM thin plate strain energy can be derived as: Reducing Equation (37) results in the following: More specifically, the strain energy expression considering the MCST of the FGM thin plate is as follows: The membrane energy expression is as follows: where Substituting Equations (29)-(31) into Equation (40) yields the solution for membrane energy.
A more specific expression is shown below: According to Hamilton's principle [52,53], among all possible displacement fields for deformation, the real displacement field minimizes the total potential energy.
The equilibrium at the x-directional displacement of u 0 and out-of-plane deflection of w 0 requires that the first-order energy partial derivatives of different directions are all zero. Thus, we can obtain the following results: Substituting Equations (38) and (40) into Equation (35), the total energy equation can be obtained: Therefore Equations (45) and (46) can be reduced to According to Equation (49), the relationship between out-plane and in-plane displacement can be written as: Considering the pre-buckling problem, when w 0 = 0, Equations (48) and (50) can be derived to From Equation (52), it can be seen that the scale effect parameter does not affect the critical buckling displacement.
Eliminating w 0 based on Equations (48) and (49), we can obtain an approximate linear post-buckling solution as follows: where (u 0 ) c is the critical displacement in the x-direction. The critical buckling force P cr can be obtained by substituting Equation (52) with Equation (51).
A more specific expression is as follows:

Numerical Results
This section investigates numerical examples of pre-buckling and post-buckling problems for FGM thin plates. The FGM is made of Al and Cu. The upper surface is Cu and the lower surface is Al. The material properties of the FGM are described as shown in Table 1.

Effects of the Scale Parameter
By considering various scale parameters, the relationship between force and x-direction displacement is plotted in Figure 3 based on CT and MCST. Equation (23) shows that the MCST can be converted into CT when scaling parameter l = 0. It can be seen from Figure 3 that when u 0 < u 0c , it represents the thin plate in the pre-buckling state. Conversely, u 0 > u 0c represents the post-buckling state. Figure 3 shows that with the same x-direction displacement, as the scale parameter l increases, the force P tends to increase. The increase in image slope with the increase in scale effect parameter indicates the increase in flexural stiffness of the thin plate. At the instant of buckling, the buckling load P decreases. The slope of the post-buckling stage is smaller than that of the pre-buckling stage, indicating that the flexural stiffness of the thin plate decreases after buckling. According to Equation (53), consider a 10 µm thick Cu plate and Al plate. The Cu plate flexural stiffness postbuckling is 0.52 times that pre-buckling based on CT, and 0.5204 times based on MCST. The Al plate flexural stiffness post-buckling is 0.52 times that pre-buckling based on CT, and 0.53 times based on MCST. The details are shown in Table 2, where D pre is the pre-buckling flexural stiffness and D post is post-buckling flexural stiffness.     For the critical buckling displacement results, the linear analysis in Figure 3 and the nonlinear simulation in Figure 4 are approximately the same. However, the critical load for nonlinear buckling is significantly smaller than that for linear buckling. Linear buckling obtains the critical buckling load based on the small deformation assumption. In the nonlinear buckling analysis, when the load is small, the analysis results are consistent with the linear analysis. As the load gradually increases, the deformation increases, the structure exhibits nonlinear response, and the force-displacement curve begins to deviate from the linear result. Therefore, the critical buckling load is less than the linear analysis result.    Figure 5 shows the linear buckling displacement response of Cu plates. The whole deformation is regular, and the maximum out-of-plane displacement is in the center of the plate.  Figure 6 shows the irregular stress distribution in post-buckling nonlinear analysis. It also exhibits irregular deformation compared with the linear post-buckling analysis. There is maximum stress at the center of the Cu plate in the nonlinear post-buckling analysis, which is also the most prone to failure. The minimum stress is in the middle of the edge of the Cu plate.  Figure 7 shows that with the same deflection, as the scale parameter l increases, the buckling strain energy tends to increase. The buckling strain energy increases as the deflection increases. Figure 8 shows the two possibilities for buckling: in-plane, and out-of-  Figure 6 shows the irregular stress distribution in post-buckling nonlinear analysis. It also exhibits irregular deformation compared with the linear post-buckling analysis. There is maximum stress at the center of the Cu plate in the nonlinear post-buckling analysis, which is also the most prone to failure. The minimum stress is in the middle of the edge of the Cu plate.  Figure 7 shows that with the same deflection, as the scale parameter l increases, the buckling strain energy tends to increase. The buckling strain energy increases as the deflection increases. Figure 8 shows the two possibilities for buckling: in-plane, and out-of-  Figure 5 shows the linear buckling displacement response of Cu plates. The whole deformation is regular, and the maximum out-of-plane displacement is in the center of the plate. Figure 6 shows the irregular stress distribution in post-buckling nonlinear analysis. It also exhibits irregular deformation compared with the linear post-buckling analysis. There is maximum stress at the center of the Cu plate in the nonlinear post-buckling analysis, which is also the most prone to failure. The minimum stress is in the middle of the edge of the Cu plate. Figure 7 shows that with the same deflection, as the scale parameter l increases, the buckling strain energy tends to increase. The buckling strain energy increases as the deflection increases. Figure 8 shows the two possibilities for buckling: in-plane, and out-ofplane buckling. Figure 9 is the contour line of membrane energy. Figures 8 and 9 indicate the variation of membrane energy with x-direction displacement and deflection. Nanomaterials 2022, 12, x FOR PEER REVIEW 14 of 24 plane buckling. Figure 9 is the contour line of membrane energy. Figures 8 and 9 indicate the variation of membrane energy with x-direction displacement and deflection.   plane buckling. Figure 9 is the contour line of membrane energy. Figures 8 and 9 indicate the variation of membrane energy with x-direction displacement and deflection.

Effects of Power-Law Index Parameter k
By considering various power-law index parameters k, the changes in Young's modulus, Poisson's ratio, and density through the thickness direction of the FGM plate are investigated and shown in Figures 10-12 for different material distributions. The powerlaw index parameter k has an important effect on the Young's modulus and Poisson's ratio of the FGM thin plate. Figure 13 shows that with the same material proportions, the critical displacement in the x-direction decreases with the decreasing power exponent k. At the same power exponent k, the critical displacement in the x-direction decreases as the proportion of the upper material increases. According to Equation (52), the scale effect does not change the critical displacement in the x-direction. However, power-law index parameter k has a significant effect on it.

Effects of Power-Law Index Parameter k
By considering various power-law index parameters k, the changes in Young's modulus, Poisson's ratio, and density through the thickness direction of the FGM plate are investigated and shown in Figures 10-12 for different material distributions. The powerlaw index parameter k has an important effect on the Young's modulus and Poisson's ratio of the FGM thin plate. Figure 13 shows that with the same material proportions, the critical displacement in the x-direction decreases with the decreasing power exponent k. At the same power exponent k, the critical displacement in the x-direction decreases as the proportion of the upper material increases. According to Equation (52), the scale effect does not change the critical displacement in the x-direction. However, power-law index parameter k has a significant effect on it.

Influence of Power-Law Distribution Considering the Scale Effect
According to Equations (39) and (50), the relationship of buckling strain energy, x-direction displacement and material proportion can be seen in Figure 14. The power-law index parameters have a significant effect on buckling strain energy. When the material ratio is determined, the value of the scale effect parameter is determined, and the buckling energy is determined accordingly. Figure 14 shows that the buckling strain energy increases with the increase in x-direction displacement. At the same proportion of material and x-direction displacement, the buckling strain energy increases with the increasing power-law index k. The membrane energy decreases as the power-law index parameter increases, as is shown in Figure 15.

Influence of Power-Law Distribution Considering the Scale Effect
According to Equations (39) and (50), the relationship of buckling strain energy, xdirection displacement and material proportion can be seen in Figure 14. The power-law index parameters have a significant effect on buckling strain energy. When the material ratio is determined, the value of the scale effect parameter is determined, and the buckling energy is determined accordingly. Figure 14 shows that the buckling strain energy increases with the increase in x-direction displacement. At the same proportion of material and x-direction displacement, the buckling strain energy increases with the increasing power-law index k. The membrane energy decreases as the power-law index parameter increases, as is shown in Figure 15.   The significance of the power-law index parameters k and scale effect parameters l is shown in Figures 16-18. However, compared with power-law index parameters k, the scale effect parameters have a greater influence on critical buckling load. From Figure 16, with the same scale parameter l and power exponent k, the critical buckling load varies little with the material ratio. Figure 15 shows that the higher the modulus of elasticity, the higher the critical buckling load. When l = 0, the maximum difference of critical buckling load is less than 0.5 N. However, it shows a large difference in Figure 18 when considering the power-law distribution of the scale effect. The maximum difference of the critical buckling load is more than 30 N, as shown in Figure 18. The scale effect parameter of the lower layer of Al is about 4.6 times that of the upper layer of Cu. Therefore, as the proportion of the lower layer Al material increases, the scale effect parameter of the FGM plate also increases, thus leading to a significant increase in the critical buckling load. The significance of the power-law index parameters k and scale effect parameters l is shown in Figures 16-18. However, compared with power-law index parameters k, the scale effect parameters have a greater influence on critical buckling load. From Figure 16, with the same scale parameter l and power exponent k, the critical buckling load varies little with the material ratio. Figure 15 shows that the higher the modulus of elasticity, the higher the critical buckling load. When l = 0, the maximum difference of critical buckling load is less than 0.5 N. However, it shows a large difference in Figure 18 when considering the power-law distribution of the scale effect. The maximum difference of the critical buckling load is more than 30 N, as shown in Figure 18. The scale effect parameter of the lower layer of Al is about 4.6 times that of the upper layer of Cu. Therefore, as the proportion of the lower layer Al material increases, the scale effect parameter of the FGM plate also increases, thus leading to a significant increase in the critical buckling load.        Present I represents the critical buckling load without considering the scale effect. The critical buckling load decreases as the power law index parameter increases. With the increase in the power law index, the volume fraction of SUS304 increases. The increase in the volume fraction of SUS304 results in a decrease in the hardness of the FGM plates, resulting in a decrease in the critical buckling load. Present II and Present III represent the critical buckling loads when considering a power-law distribution with scale effects. There are no experimental data and literature on the scale effect parameters of FGM, so we assume a scale effect parameter of 1 µm for SUS304 and a scale effect parameter of 2 µm for Si 3 N 4 in Present II, and assume a scale effect parameter of 2 µm for SUS304 and a scale effect parameter of 1 µm for Si 3 N 4 in Present III as a comparison. Present II shows that when the Si 3 N 4 volume fraction increases, the elastic modulus and scale effect parameters of the FGM increase, resulting in an increase in the critical buckling load compared with Present I. Present III shows that when the Si3N4 volume fraction increases, the elastic modulus of the FGM increases but the scale effect parameter decreases, resulting in less fluctuation in the critical buckling load.

Concluding Remarks
Based on the modified coupled stress theory, the post-buckling of an Al-Cu FGM thin plate is analyzed. According to Hamilton's principle, deflection solutions and stress solutions for post-buckling under-scale effects are derived. The influence of the scale effect on the buckling strain energy and membrane energy is investigated in detail. The most significant conclusions can be summarized as follows:

•
Under the same conditions, the critical buckling load, the membrane energy, and buckling strain energy increase as the scale effect parameter increases. The critical x-direction displacement is not affected by scale effects.

•
The critical buckling displacement for nonlinear buckling is consistent with the linear buckling analysis, but the critical buckling load is smaller considering the nonlinearity of the structure • Under the same circumstances, the modulus of elasticity, Poisson's ratio, critical buckling load, critical buckling displacement, and buckling strain energy all increase with increasing power-law index parameter k. The membrane energy decreases as the power-law index parameter increases. As the power-law index k increases, the volume fraction of bottom copper increases. The increased modulus of elasticity of the FGM plate results in an increase in stiffness and thus an increase in the critical buckling load.

•
Both the scale effect parameter l and the power-law exponent k have important effects on the FGM thin plate buckling and post-buckling problems, but in comparison, the scale effect parameter is more influential than the power-law exponent.

•
The scale effect parameter can greatly increase the critical buckling load, and the correct choice of material is important for the stable design of the device. The scale effect cannot be ignored, and the research on the scale effect is of great significance for MEMS/NEMS design, and manufacture.
It is of great significance to consider continuously distributed scale effect parameters for micro-scale FGM. Both scale effect parameters and power-law index parameters have important effects on the mechanical behavior of micro-scale FGM. The scale effect parameter affects the stiffness of the material, and the power-law index parameter affects the intrinsic properties of the FGM. The present results will be used for the design of micro-scale FGM MEMS/NEMS, and it is necessary to consider both the MSCT and Power-law in the design of micro-scale FGM.