A Nonlocal Strain Gradient Approach for Out-of-Plane Vibration of Axially Moving Functionally Graded Nanoplates in a Hygrothermal Environment

School of Rail Transportation, Soochow University, Suzhou 215131, China Engineering Research & Development Centre for Underground Technology of Jiangxi Province, School of Civil Engineering and Architecture, East China Jiaotong University, Nanchang, Jiangxi 330013, China Guangxi Key Laboratory of Cryptography and Information Security, Guilin University of Electronic Technology, Guilin 541004, China


Introduction
With the rise of micro/nanotechnologies, the development and application of micro/nano-electromechanical systems (MEMS/NEMS) have developed rapidly, such as micro/nanosensors, micro/nanoresonators, smart wearable devices, and medical nanorobots [1,2]. e systematic and in-depth understandings of the mechanical properties of these devices and their key micro/nanocomponents are needed to provide theoretical guidance for the mechanical-based design and regulation [3,4]. So far, there have been numerous investigations on this topic, especially the studies on the mechanical behavior of uniform nanomaterials and structures based on the molecular dynamic simulation or nonclassical continuum approaches [5][6][7][8][9]. ese results play an important role in perfecting performances and promoting the industrialization of MEMS/NEMS [10] and cellular mechanics [11]. However, while the MEMS/NEMS work in some harsh environments, such as high temperature, high pressure, humidity and heat, high-strain rate, and highspeed impact, the nanostructure composed of a single uniform material gradually fails to meet the engineering requirements [12]. Introducing the concept of functionally graded materials at macroscale into micro/ nanoscale and constructing functionally graded nanostructures is one of the effective ways to overcome the problems in extreme working environments [13][14][15]. Functionally graded nanostructures, as one of the components in nanodevices or nanosystems, may produce an unwanted vibration caused by external excitations or their own factors during work. Hence, the studies on the vibration characteristics of functionally graded nanostructures are of great significance.
To study the dynamic behaviors of micro/nanodevices composed of functionally graded nanostructures, there are generally three kinds of theoretical approaches including the discrete methods [16], modified continuum theories [17], and atom-continuum multiscale models [18]. e discrete methods refer to that the micro/nanomaterials and structures are equivalent to a number of molecules or atoms, and the dynamics of materials and structures can be simulated by calculating the interaction between molecules or atoms. For a special micro/nanosystem with a large number of molecules/atoms, the application of discrete methods may take plenty of time with a higher requirement of computing hardware. erefore, the modified continuum theories that regard the nanomaterials and structures as generalized continuous structures are more suitable for this case. Using the traditional analysis thought and modeling procedure (e.g., Hamilton's principle) of the macroscopic continuum, the mechanical properties of nanomaterials and structures can be predicted by establishing the nonclassical continuum constitutive relations containing internal characteristic scale parameters of the nanomaterials. Although the modified continuum theories still adopt the basic ideas of a classical macroscopic continuum, the constitutive relations of macroscopic continuum mechanics have not been duplicated directly but improve them to include new phenomena at the nanoscale, such as a long-range interaction between atoms and strains with higher-order gradient. e atom-continuum multiscale models aim to establish an atomic calculation framework on the micro/nanoscale and a continuum framework on the macroscale according to different geometric sizes in different dimensions of micro/nanomaterials and structures, so as to develop the spatial multiscale model. It may also introduce temporal multiscale into the dynamic analysis of atom-continuum multiscale models. In the present study, the dynamics of two-dimensional functionally graded nanostructures are investigated by using the modified continuum theories. Hence, the nonclassical continuum model is applied to characterize the inherent features at the nanoscale. Since functionally graded nanostructures are composed of inhomogeneous materials and there are a huge number of molecules/atoms contained in two-dimensional functionally graded nanostructures, it is relatively difficult to deal with inhomogeneous materials and structures by molecular dynamic simulation, atomic simulation, or atomcontinuum multiscale models. e modified continuum theories mainly include the micropolar and micromorphic theory [19], the couple stress theory [20], the strain gradient theory [21], and the nonlocal theory [22], which have been widely applied in the mechanical analysis of micro/nano-structures. For example, Li et al. explored the size dependence of nanorods [23] and nanobeams [24] based on the nonlocal integral model, and the results show that the nonlocal effect was more sensitive in the thickness direction, which implies that the thickness effect plays a dominant role in the contribution of size dependence. Jalaei and Civalek [25] investigated the dynamic instability of viscoelastic porous functionally graded nanobeams embedded on the visco-Pasternak medium according to Eringen's differential law, and the model is exposed to the magnetic field as well as an axial oscillating load. Lu et al. [26] performed the factors influencing the postbuckling properties of multilayer microtubules reinforced by functionally graded graphene platelets containing initial geometrical imperfection via the modified couple stress theory.
In recent years, with the continuous deepening of the research, these modified continuum theories have shown some disadvantages [27]. For example, the internal characteristic scale parameters of the same nanomaterial are different in different studies due to the inconsistency in the determination methods. For another example, there have been different kinds of nonlocal effects predicted. e equivalent stiffness softening and stiffness hardening models were established, and even some paradoxes were found in the literature. On the basis of these theoretical models, Lim et al. [28] proposed a nonlocal strain gradient theory. By introducing three internal characteristic scale parameters, the nonlocal effects of classical strain and strain gradients, as well as the gradient effects of traditional Eringen's nonlocal stress and high-order nonlocal stresses are comprehensively considered in the theory. erefore, the theoretical predictions by the nonlocal strain gradient theory were satisfactorily consistent with the experimental results [29]. Recently, the nonlocal strain gradient theory has become one of the most widely applied nonclassical continuum approaches [30]. For instance, Saffaria et al. [31] studied the dynamics of a cantilever nanotube resting on a nonlinear viscous Pasternak foundation in the context of the nonlocal strain gradient theory, in which the viscoelastic fluid was conveyed, and the magnetic field and harmonic loads were also taken into account. Based on the nonlocal strain gradient theory, Li et al. [32] examined the effect of size-dependent parameters on the stiffness of functionally gradient beams using the Timoshenko beam model. To sum up, the strain gradient theory is applied to analyze the free vibration of functionally graded nanoplates so as to obtain accurate mechanical results.
Generally, the functionally graded nanostructures can be modeled as one-dimensional beam/tube/rod structures, two-dimensional plate/shell/film structures, and three-dimensional solid structures according to their size characteristics. Among them, the one-dimensional functionally graded structures have been fully studied [33]. In terms of functionally graded nanoplates plate/shell/film, there is also some research literature. For example, Lu et al. [34] developed a functionally graded cylindrical shell model in the framework of the nonlocal strain gradient theory and the Gurtin-Murdoch surface elasticity theory. ree typical size effects including nonlocal effect, strain gradient effect, and surface energy effect on the free vibration were carried out. Vinh and Tounsi [35] analyzed the influence of the spatial variation of the nonlocal parameter on the free vibration of functionally graded sandwich nanoplates by modifying the classical nonlocal theory, and the motion equation was established based on the new inverse hyperbolic shear deformation theory and Hamilton's principle. e rigid body motion of the model itself is not taken into account in most literature on nanoplates. However, there are some phenomena similar to the motion at the macroscale [36] in nanosystems. For instance, medical nanorobots often need to move in the human blood as a whole to complete the auxiliary diagnosis and treatment of diseases through structural vibration. At this moment, kinematics and dynamics are coupled with each other. Of course, the out-ofplane vibration of a plate is more obvious than the in-plane vibration, and the nanoplate is no exception.
As mentioned above, functionally graded nanostructures tend to be designed for extreme working environments, and taking NEMS (e.g., nanosensors/resonators/nanorobots) as an example, they often need to work in some humid and hot environments. erefore, it is of great significance to investigate the vibration characteristics of functionally graded nanoplates in the hygrothermal environment. In the investigation of mechanical behaviors of functionally gradient nanoplates, the hygrothermal effects have been attended by many researchers. Refrafi et al. [37] analyzed the buckling responses of simply supported functionally graded sandwich plates rested on the Winkler-Pasternak elastic foundation subjected to hygrothermal and mechanical loads. Tounsi et al. [38] proposed a simple four-variable trigonometric integral shear deformation model to perform the mechanical properties of advanced functionally graded ceramic-metal plates seated on a two-parameter elastic foundation, which exposed to a nonlinear hygro-thermomechanical load.
ere are some performed studies on the axially moving nanoplates. In [39,40], the nonlocal theory was adopted to analyze the free vibration of axially moving nanoplates, but the strain gradient effect was not taken into account. In [41,42], the nonlocal strain gradient theory was employed to study the mechanical behaviors of axially moving nanoplates. However, sometimes, it is necessary to consider the hygrothermal effects because of the extreme working environments of functionally graded nanomaterials. For this purpose, the hygrothermal effect on the axially moving functionally graded nanoplates is considered herein. e novelty of this research is to explore the effects of the axial speed, gradient index, and internal and external characteristic parameters on the vibration of axially moving functionally graded nanoplates exposed to hygrothermal environments. In addition, the coupling effect between the internal characteristic parameters is also considered. e remaining parts of the present work are arranged as follows: In Section 2, the theoretical model of the axially moving functionally graded nanoplate subjected to hygrothermal loads is established, and the effect of the gradient index on the bending stiffness is demonstrated. en, the partial differential governing equations describing the out-of-plane vibration are developed based on Hamilton's principle in the context of the nonlocal strain gradient theory. In Section 3, in order to calculate the natural frequencies, the composite beam function method in conjunction with the complex mode approach is employed according to the boundary conditions, and then the correctness of the approach is verified by comparisons. Afterwards, the effects of temperature variation, moisture concentration, axial speed, gradient index, nonlocal parameter, material characteristic scale parameter, and aspect ratio on the out-of-plane vibration frequency are analyzed via numerical examples. Finally, the summary is concluded in Section 4. e results of the present work are of reference significance for the design, optimization, and control of axially moving MEMS/ NEMS based on functionally gradient nanoplates, such as nanorobots and nanosensors subjected to hygrothermal environments.

Theoretical Constitutive and
Governing Equations e core constitutive equations of the nonlocal strain gradient theory were expressed in an integral form initially [28]. However, the spatial integration of strains and its gradients is difficult to be solved. Hence, an equivalent differential constitutive equation was transformed for a two-dimensional nanostructure as where t is a total stress tensor, ε ′ is a classical strain tensor, C is an elastic tensor, ∇ 2 is a Laplace operator, a is a material internal scale, l is a strain gradient characteristic scale, and e is a nonlocal material parameter. As depicted in Figure 1, given an axially moving functionally graded nanoplate with length L a , width L b , and thickness h, the x-y plane coordinate system is established in the geometric middle plane and z coordinate is along the thickness of the nanoplate. e axial speed is denoted by υ.
e functionally graded nanoplate is exposed to the hygrothermal environment with uniform temperature changing ΔT (K) and moisture concentration rise ΔH (%), where ΔT � T− T 0 , in which T is the applied temperature and T 0 � 300K is the reference temperature; ΔH � H− H 0 , in which H is the applied moisture concentration and H 0 � 0 is the reference moisture. In the present work, we assume that the nanoplate is composed of the ceramic and metal, namely, the inner compositions of the two-dimensional thin nanostructure are the mixture of two materials by the gradient law so that the characteristics vary regularly according to the gradient index. erefore, Young's modulus E, mass density ρ, thermal expansion coefficient α, and moisture expansion coefficient β of the functionally graded nanoplate can be respectively expressed as [43] where p is the gradient index and two subscripts c and m denote the material properties of ceramic and metal phases, respectively.

Shock and Vibration
Because the changing trend of the above four types of material properties is identical, only taking Young's modulus of functionally graded nanoplates as an example, the variation along the thickness of the nanoplate is shown in Figure 2 under different gradient indexes, where E c � 380 GPa and E m � 70 GPa are adopted for ceramic and metal materials, respectively.
Let u, v, and w represent the displacements of the nanoplate along x, y, and z axes, respectively. Because the physical neutral layer is not equal to the geometric middle plane (x-y plane coordinate established herein) for an inhomogeneous functionally graded nanoplate, the position coordinate of the physical neutral layer is en, by applying the classic plate theory [44] and the concept of physical neutral plane [45], the displacement fields can be expressed as According to the relationship between the classic Kirchhoff thin plate theory and the nonlocal strain gradient differential constitutive equation, the physical equations of the functionally graded nanoplate can be acquired as where the functionally graded nanoplate is regarded as a plane stress problem, B � E(z)/(1 − μ 2 ), and μ is Poisson's ratio.
It is noted that the definition of the bending moment and its relationship with stress are valid in both macroscale and microscale. erefore, the defining expression of the bending moment in the nonlocal strain gradient theory is consistent with that in the classical continuum theory. e difference is that the stress in the definition is nonlocal stress and classical stress, respectively. en, by applying the classic plate theory [44], the bending moment can be expressed as According to equation (3), the relationship between the bending moment and the displacement can be derived from equation (4) as where w represents the out-of-plane displacement and D � is the equivalent bending stiffness. e effect of the gradient index on the equivalent bending stiffness is shown in Figure 3.
As clearly visible in Figure 3, when the gradient index increases, the equivalent bending stiffness decreases quickly at first and then gradually slows down, which indicates that a higher gradient index reduces the equivalent bending stiffness of functionally graded nanoplates. e observations can provide reference for the design of functionally graded nanoplates subjected to different working conditions.
To analyze the out-of-plane vibration of axially moving functionally graded nanoplates in hygrothermal environments, Hamilton's principle is utilized to derive the governing equation of motion, which requires where δU, δV, and δK are the variations of the strain potential energy, external potential energy, and kinetic energy, respectively, and t 1 and t 2 represent two arbitrary moments. Firstly, the strain potential energy for the nanoplate may be calculated as where A is the region of the middle plane of the nanoplate. en, we assume that the temperature changing results in the external potential energy which is reverse to work done by external loads, and it can be written as where N x � N y � N T + N H are the hygrothermal loads, in which N T is the temperature load and N H is the humid environment load, respectively. And, the hygrothermal resultant can be expressed by in which Finally, the kinetic energy of the nanoplate is given by Shock and Vibration 5 [45]. Substituting equations (7), (8), and (11) into equation (6), and then integrating by parts, the partial differential equation that governs the out-of-plane vibration can be written as According to the variational principle, the essential/ natural boundary conditions of the model along x � 0, L a can be obtained as And the boundary conditions along y � 0, L b are acquired as Substituting equation (5) into (12) and the governing equation in terms of out-of-plane displacement can be expressed as In order to obtain the dimensionless governing equation and analyze the out-of-plane vibration characteristics of axially moving functionally graded nanoplates qualitatively, the following variables and parameters are necessary: Accordingly, the dimensionless form of the governing equation (15) can be rewritten as ζ 2 z 6 w zx 6 + 3ζ 2 λ 2 z 6 w zx 4 zy 2 + 3ζ 2 λ 4 z 6 w zx 2 zy 4 + ζ 2 λ 6 z 6 w zy 6 Shock and Vibration e spatial variables x and y and time variable t are coupled in equation (17). In the linear elasticity and small deformation, one can separate these variables. e solution of out-of-plane displacement is set to w(x, y, t) � W(x, y)e iω mn t , (18) where ω mn is the vibration circular frequency, usually a complex frequency, of which the real and imaginary components can be used to describe the vibration frequency and related stability, respectively. Substituting equation (18) into (17), the nondimensional frequency-domain governing equation can be obtained as

Numerical Results and Discussion
Unlike most research studies applying Galerkin method to investigate the moving nanosystems [39,41,46], the present work employs the composite beam function method [47] to develop the dynamic stability of axially moving nanoplates. e approach assumes the displacement function according to boundary conditions, so that the governing equation can be solved numerically directly. Compared with other accurate solution methods such as Galerkin and Navier's solution, the advantage of the present method is that it can analyze systems with various boundary conditions.
In addition to the governing equations of axially moving functionally graded nanoplates established in the previous section, the boundary conditions are needed to solve the outof-plane vibration theoretical model. Taking the edge constraint of simply supported on four sides as an example, according to equations (13) and (14), the corresponding natural boundary conditions at x � 0, L a are chosen as M xx � 0, and the natural boundary conditions at y � 0, L b are chosen as M yy � 0. en, considering the relationship between bending moment and displacement in equation (5), the details of the basic boundary conditions of the nanoplates can be obtained as 8 Shock and Vibration erefore, the dimensionless boundary conditions can be written as Utilizing the composite beam function method [47], the dimensionless out-of-plane displacement can be assumed as where ϕ m (x) and φ n (y) are the adopted mode functions of x and y directions, respectively, and correspondingly, ψ mn (x, y) is the product of two mode functions.
For the nanoplate simply supported on four sides, the ydirection mode function can be defined as a sine function that satisfies the boundary conditions in equation (21). Hence, the product of mode functions can be expressed as Substituting equation (23) into (22), and further into equation (19) yields e solution to equation (24) can be written in form as ϕ m (x) � Ce iβx , and the following equation can be obtained by substituting the solution into equation (24): As a sixth-order ordinary differential equation, the solution to equation in (25) can be assumed as

Shock and Vibration
where C jm (j � 1, 2, 3, . . .. . ., 6) denotes six unknown coefficients introduced by integral, and β jm (j � 1, 2, 3, . . .. . ., 6) are corresponding six roots of algebraic equation (25). According to equation (21), the boundary conditions can be written as As a result, a set of algebraic equations are acquired by substituting equation (26) into the simplified boundary condition equation (27). To ensure the existence of the nontrivial solutions of the coefficients or mode function in equation (26), the determinant of the coefficient matrix in the algebraic equations should be set to zero.
Combining equations (25) and (28), the dimensionless vibration frequency can be determined numerically. In the calculation process, the other parameters are set as constant values firstly; then, an interval of the remaining parameter is provided, so that the variation of vibration frequencies with respect to the remaining parameter can be achieved via a numerical calculation program. Consequently, the relationships between the dimensionless vibration frequency and each parameter can be obtained.
In order to demonstrate the accuracy and validity of the present analyses, a benchmark result is given by the comparison with previous results available in the literature. First, by ignoring the strain gradient, axial velocity, and hygrothermal environment effects, the present results are compared with those reported by Natarajan et al. [48] and Jandaghian and Rahmani [49] in Table 1. e results agree well with those from the literature, which proves that the model and method proposed in this paper are reliable in vibration analyses. In addition, the influence of uniform thermal effects on nanoplates is also compared in Table 2 [50]. e verification example reveals that the proposed approach can accurately predict the behavior of nanoplates under the influence of external environmental loads.
In the present examples, we take n � 1, and the material properties of the functionally graded nanoplate are shown as follows [51].
As for the magnitude of the nonlocal parameter ea and the strain gradient parameter l, some literature [52,53] has made a special discussion, but there are no experimental methods to determine it. In the physical sense, these two parameters represent the continuous range of nonlocal longrange action and the strength of strain gradient influence. In the present work, following the method used by Lu et al. [54], the two length parameters are normalized by length L a to obtain numerical results, and it is used that both ea/L a and l/L a are in the range of 0 to 0.2.
Firstly, the influences of the hygrothermal environment on vibrating functionally graded nanoplates are examined.
Taking L a /h � 10, τ � 0.1, ζ � 0.1, λ � 1, ξ � 1, p � 1, and ΔH � 0, 1%, and 2%, respectively, we reveal the variation of dimensionless first-and second-order vibration frequencies with temperature changing. It can be seen from Figures 4(a) and 4(b) that, the smaller the moisture concentration is, the greater the vibration frequency is. In addition, when the moisture concentration is fixed, the dimensionless vibration frequency decreases with the increase of temperature changing. For example, when the temperature changing increases from 0 to 200 K and the moisture concentration remains unchanged compared with H 0 , the first-order vibration frequency decreases by 9.3% and the second-order decreases by 3.56%. is is because the high temperature leads to expansion in nanoplate structures, which reduces the structural stiffness. Consequently, it is believed that the influence of temperature changing should be paid more attention to in lower-frequency vibration stage.
Subsequently, vibration frequencies versus continuous moisture concentrations are presented in Figures 5(a) and 5(b). Let ΔT be 0, 50 K, and100 K, respectively, and other fixed parameters remain the same as Figures 4(a) and 4(b). e change of vibration frequencies with increasing moisture concentration can be observed. It can be revealed that the dimensionless vibration frequency is smaller when the temperature is higher. For example, if the temperature is kept as a constant at T 0 , when the moisture concentration changes from 0 to 2%, the first-order vibration frequency decreases by 12.25% and the second-order decreases by 4.63%. e influence of hygrothermal environmental factors on fundamental frequencies of vibrating nanoplates is indispensable in the design of nanoplates working in special environments. e decrease of vibration frequency means that the equivalent stiffness of the nanoplate is weakened. erefore, an increase in the moisture concentration may reduce the stiffness of functionally graded nanoplates. is is because the nanoplate structure degrades and swells with the adsorption of water molecules. Hence, in the actual design of NEMS, the damage effect of the high temperature and moisture on the structure should be considered.
In addition, the effect of the axial speed of functionally graded nanoplates on the vibration is considered. Let other parameters remain unchanged. For example, take ΔT �100 K, ΔH � 1%, and p � 1. en, let τ � 0.01, ζ � 0.05; τ � 0.05, ζ � 0.05; and τ � 0.05, ζ � 0.1, respectively. e trend of the dimensionless vibration frequency with the increasing axial speed is presented in Figures 6(a) and 6(b). It can be seen that the vibration frequency increases when the nonlocal parameter becomes smaller or the material  Shock and Vibration 11 characteristic scale parameter becomes larger. At the same time, it is apparent that with the increase of the axial speed, the vibration frequency decreases faster and faster. For example, choosing τ � 0.05 and ζ � 0.05, when the dimensionless speed increases from 0 to 2, the first-order vibration frequency decreases by 9.36%, while the dimensionless speed increases from 2 to 4, the first-order frequency decreases by 59.99%, and so the decline rate increases significantly.
Hence, the increase of axial speed leads to the instability of functionally graded nanoplates, and such a conclusion is a common sense which also holds at the macroscale. In the actual NEMS design, the appropriate speed should be chosen to ensure the stability of the system by considering the reduction effect of speed on the structure stiffness. Moreover, it can be observed that, when the nonlocal parameter is a constant, there is a significant difference between the   erefore, there is a remarkable effect of the material characteristic scale parameter on the stability of the functionally graded nanoplates. Unlike the axial speed, the material characteristic scale parameter belongs to an internal influence factor, which is particularly significant at the nanoscale but can be ignored at the macroscale. It can also be seen from the Figure 6 that, for different nonlocal parameters and the material characteristic scale parameters, the curve has basically the same change trend. is shows that internal influence factors basically have no influence on the effect of axial speed changes on the system stability, which also indicates that the internal characteristic parameters affect little on the system at the macro scale.
As the gradient index increases, the properties of nanoplates gradually change from ceramic to metal. e effect of the gradient index on the dimensionless natural frequency under different material characteristic scale parameters is described in Figures 7(a) and 7(b). It can be seen that the natural frequency decreases as the gradient index increases, and for smaller material characteristic scale parameters, the rate of decrease is greater. In other words, when the material characteristic parameter is small, the change of the gradient index has a stronger influence on the vibration behaviors of the nanoplate. In the range of the gradient index from 0 to 4, the vibration frequency is constant and the first-order frequency is always lower than the second-order frequency, indicating that the nanoplate system is stable in this range. Similarly, for both the first-and second-order frequencies, the frequencies are significantly different as the material characteristic scale parameters increase, which proves the size dependence of the nanoplates. For different d material characteristic scale parameters, the slope of the curve in Figure 7 is different when it drops, which indicates that the material characteristic scale parameters will affect the stiffness effect of the gradient index on the system. en, the influence of the nonlocal parameter on the vibration stability is analyzed. e other parameters are kept unchanged and then take ξ � 1, ζ � 0.05; ξ � 2.5, ζ � 0.05; and ξ � 2.5, ζ � 0.1, respectively. e changes of the dimensionless vibration frequency with increasing the nonlocal parameter under different dimensionless axial speeds and material characteristic scale parameters are plotted in Figures 8(a) and 8(b). It can be observed that the vibration frequency is larger with a smaller speed or a larger material characteristic scale parameter. Besides, it is clear that the dimensionless vibration frequency decreases with increasing the nonlocal parameter. Taking the case of ξ � 1 and ζ � 0.05 as an example, when τ increases from 0 to 0.2, the dimensionless first-order vibration frequency decreases by 34.99%, and the secondorder decreases by 48.24%, which reflects the softening behavior in the nonlocal effect. Meanwhile, we observe that, when ξ � 2.5 and ζ � 0.1, the calculated results are in the middle in Figure 8(a), but they are at the top in Figure 8(b). is is because the increase of stiffness caused by the increase of material characteristic scale parameter is greater than the stiffness weakening effect by the increase of the axial speed.
Furthermore, the influence of the nonlocal parameter on the vibration of functionally graded nanoplates with different rectangular shapes is discussed. e aspect ratio of the nanoplate λ is chosen as 1, 1.05, and 1.1, respectively. As clearly visible in Figures 9(a) and 9(b), the dimensionless vibration frequency changes with the nonlocal parameter under different aspect ratios. While the nonlocal parameter changes from 0 to 0.2, if the aspect ratio is 1, the first-order frequency decreases by 33.2% and the second-order frequency decreases by 46.26%. If the aspect ratio is 1.05, the first-order frequency decreases by 33.74% and the second-order frequency decreases by 46.58%. If the aspect ratio is 1.1, the first-order frequency decreases by 34.29% and the second-order frequency decreases by 46.90%. It is obvious that, in the case of higher-order vibration, when the aspect ratio increases, the vibration frequency decreases more and   more significantly. e vibration of the nanoplate can be controlled by changing the shape of the nanoplate to avoid system resonance and ensure the stability of the structure. en, the influence of the material characteristic scale parameter is investigated. Keep the aspect ratio λ � 1, and then select different axial speeds and nonlocal parameters as comparison; that is, taking ξ � 1, τ � 0.05; ξ � 2, τ � 0.05; and ξ � 2, τ � 0.1, respectively. Figures 10(a) and 10(b) present the variation of vibration frequencies with the material characteristic scale parameter. Similar to previous conclusions, a lower axial speed and a smaller nonlocal parameter result in a higher vibration frequency. In addition, it can be seen that vibration frequencies increase with an increase in the material characteristic scale parameter. For instance, if ξ � 1 and τ � 0.05, when the material characteristic scale parameter increases to 0.2, then the first-order vibration frequency increases by 30.27% and the second-order increases by 79.42%. is indicates that the stiffness of the nanoplate becomes larger while increasing the material characteristic scale parameter, which reflects the stiffness hardening effect. Besides, it is found that the influence of the material characteristic scale parameter on higher-order frequencies is greater than that on lower-order ones with different nonlocal parameters, reflecting a coupling effect between the material characteristic scale parameter and nonlocal parameter.
Finally, based on the coupling effect of the nonlocal parameter and the material characteristic scale parameter mentioned above, keeping ξ � 1 and λ � 1 invariant, the different material characteristic scale parameters are taken to observe the changes of vibration frequencies with the internal scale parameter ratio in Figure 11. It is concluded that vibration frequencies become greater in case of increasing the material characteristic scale parameter when the ratio τ/ζ is less than 1 (i.e., τ < ζ). Herein, the material characteristic scale parameter representing the strain gradient effect plays a dominant role, and the increase of vibration frequencies means the stiffness hardening. On the other hand, vibration frequencies diminish with an increase in the material characteristic scale parameter in the case of τ > ζ. Herein, the nonlocal parameter representing the nonlocal effect dominates the mechanism of stiffness variation and thus exhibits a softening behavior. In particular, the vibration frequency is restored to that of the classical thin plate vibration theory when the ratio τ/ζ is equal to 1 (i.e., τ � ζ). e stiffness softening and hardening behaviors were commonly seen in the nanostructures. In fact, the coupling effect of the nonlocal parameter and material characteristic scale parameter results in the behaviors of stiffness softening and hardening. In the nonlocal strain gradient theory, two nanoscaled internal scale parameters including the nonlocal parameter and the material characteristic scale parameter are in the same magnitude, and the nonlocal effect and strain gradient effect respectively represented by the two parameters are equivalent. Moreover, the two kinds of internal scale effects are coupled in out-of-plane vibrations of axially moving functionally graded nanoplates in hygrothermal environments as analyzed before. In this coupling interaction, the effect of a nanoscaled nonlocal parameter is almost equal to that of the nanoscaled material characteristic scale parameter. is explains why the boundary between the softening behavior and hardening behavior is determined by the magnitude relationship between the two internal parameters.
In this study, the competition relationship between different internal characteristic parameters under the  external elastic deformation and rigid motion is revealed for the first time, which can provide a reference for the design of MEMS/NEMS with a macroscopic motion in hygrothermal environments. Due to the coupling effect, the nonlocal effect and strain gradient effect should be considered comprehensively to control the stability of the system in an actual situation. It should also be mentioned that although the coupling phenomenon in the nonlocal and strain gradient effects is observed, the influences of internal characteristic parameters and external physical parameters on axially moving hygrothermal functionally graded nanoplates are independent.

Conclusions
Based on the nonlocal strain gradient theory, the functionally graded nanoplate with an axial motion in hygrothermal environments is investigated. e governing equation of motion for a vibrating and moving functionally graded nanoplate is derived via Hamilton's principle. e out-of-plane vibration characteristics of functionally graded nanoplates with boundary constraints simply supported on four sides are presented by solving the frequency-domain equation numerically. e effects of the dimensionless parameters on the vibration frequencies such as the temperature changing, moisture concentration, axial speed, gradient index, nonlocal parameter, material characteristic scale parameter, and aspect ratio are discussed. e following conclusions are obtained: (i) e effect of hygrothermal environmental factors on fundamental frequencies of vibrating nanoplates is obvious and the increase of moisture concentration and temperature change reduces the vibration frequency.
(ii) e variation of axial speed plays an important role in the stability of axially moving nanoplates. With the increase of axial speed, the vibration frequency decreases faster and faster. Hence, the increase of speed leads to the instability of the system.  (iii) e shape of the nanoplate and the mixing mode of its internal components have effects on the dynamic stability of the system. e greater the aspect ratio is, the higher the vibration frequency is. When the gradient index becomes larger, the natural frequency gets decreasing. (iv) e vibration frequencies decrease with the increase of the nonlocal parameter and increase with the increase of the material characteristic scale parameter. e two internal characteristic parameters influence and restrict each other. If the nonlocal parameter is greater than the material characteristic scale parameter, the nonlocal effect plays a dominant role, and the stiffness softening is observed. If not, the strain gradient effect dominates, and the stiffness hardening is observed. Hence, the competition relationship existing in internal characteristic parameters at the micro/nanoscale is revealed.
Data Availability e data generated or analyzed during this study are included within the article.

Conflicts of Interest
e authors declare that there are no conflicts of interest.

Authors' Contributions
Chengxiu Zhu and Jianwei Yan contributed equally to this work.