Investigating parametric homogenization models for natural frequency of FGM nano beams

: This research focuses on exploring the free vibration behavior of functionally graded (FG) nano-beams. To calculate the effective properties of the FG nano-beam, which varies solely in the thickness direction, the four homogenization schemes Mori-Tanaka, Tamura, Reuss and Voigt are employed. This study employs high-order shear deformation nano-beam theory and derives the governing equations of motion using nonlocal differential constitutive relations of Eringen. Hamilton’s principle is utilized in conjunction with the refined three variables beam theory. The consideration of a length scale parameter accounts for small-scale effects. Analytical solutions are obtained for a simply supported FG nano-beam and compared with existing literature solutions. The research also investigates the influence of different homogenization schemes, the nonlocal parameter, beam aspect ratio and various material compositions on the dynamic response of the FG nano-beam.


Introduction
Nano-beam materials are composites made up of nanotubes, nanowires and nanoparticles which are manipulated on a nanoscale level and used in a variety of applications.Researchers have utilized Eringen's nonlocal elasticity theory to investigate the static deformation, free vibration, bending, buckling and wave properties of nano-beams and other structures like nano-plates and curved nanoshells [1][2][3].
Recent research on nanobeam functionally graded beams has gained traction, with various studies conducted to explore the current state of this field.This discussion provides an in-depth analysis of the published articles, offering insight into potential applications of this technology.By examining current trends, potential improvements and applications, this research can help to further our understanding of nanobeam functionally graded beams and their utility.In their study, Vinh and Tounsi [4] investigated the free vibration of functionally graded doubly curved nanoshells using nonlocal first-order shear deformation theory.They compared the frequencies obtained to existing results and examined the effects of nonlocal parameters and the power-law index on the nanoshells.Additionally, in their later work [5] they determined the effect of a spatial variation in the nonlocal parameter on the free vibration of functionally graded sandwich (FGSW) nanoplates, finding that the natural frequencies were significantly reduced.Cuong-Le et al. investigated the linear and nonlinear solutions of a sigmoid functionally graded material nanoplate with porous effects, using isogeometric finite element formulation to examine the effects of power indexes, aspect ratios, nonlocal and strain gradient parameters [6].Liu et al. sought to understand the dynamic deflection response of an exponentially functionally graded material (E-FGM) nanoplate embedded in a visco-elastic foundation when subjected to a moving load, utilizing third-order shear deformation theory, Eringen nonlocal elasticity and a state-space method to analyze the influence of nonlocality, volume fraction index, porosity index, visco-elastic foundation coefficients and velocity and time span of a moving load on the forced vibration of the nanoplate [7].Faghidian and Tounsi studied wave dispersion and free vibrations of elastic nano-beams through the mixture unified gradient theory, providing insight into the augmented elasticity theory's effectiveness in softening and stiffening responses [8].John Peddieson utilized nonlocal elasticity theory to develop a nonlocal Benoulli/Euler beam model for cantilever beams [9], and Xu Mingtian [10] applied the same theory in combination with an integral equation approach to investigate the free transverse vibrations of nano-to-micron scale beams.Rebai Billel has studied the effects of homogenization models on the stress, deflection and vibration behavior of functionally graded plates subjected to thermal and mechanical loads [11][12][13].Reddy et al. have evaluated the static bending, vibration, and buckling responses of beams with various boundary conditions using nonlocal differential constitutive relations of Eringen.They have presented analytical solutions and numerical results to show the effect of nonlocal behavior on deflections, buckling loads and natural frequencies of beams [14,15].Zhang and Qing introduced well-posed two-phase nonlocal models for nanobeam vibrations based on higher-order shear deformation theory.Their work also encompassed the examination of refined shear deformation beams' buckling behavior under thermal loading [16].Additionally, Zhang, Schiavone and Qing extended their research to encompass hygro-thermal vibration analysis of nanobeams.This analysis employed stress-driven nonlocal theory and incorporated two-variable shear assumptions [17,18].Ebrahimi has examined the vibration characteristics of magneto-electro-thermo-elastic functionally graded nano-beams in the framework of third order shear deformation theory [19], while Thai et al. have presented a nonlocal sinusoidal shear deformation beam theory for the bending, buckling and vibration of nano-beams [20].Eltaher et al. have conducted a free vibration analysis of size-dependent functionally graded nano-beams using the finite element method based on the nonlocal continuum model of Eringen [21] and Nazemezhad et al. have investigated the nonlinear free vibration of functionally graded nano-beams with immovable ends using nonlocal elasticity and Euler-Bernoulli beam theory [22].Ebrahimi has studied the static stability and free vibration characteristics of functionally graded nano-beams using a nonlocal couple stress theory [23]; Hadji et al. have investigated the free vibration of porous functionally graded nano-beams using a new nonlocal hyperbolic shear deformation beam theory [24]; Youcef et al. have focused on the behavior of non-local shear deformation beam theory for the vibration of functionally graded (FG) nano-beams with porosities [25]; and Shariati et al. have examined the size-dependent vibrations and stability of axially functionally graded (AFG) nano-beams using a nonlocal elasticity theory [26].
Other research conducted by Cornacchia et al. delved into the analytical solutions pertaining to the linear vibrations and buckling of nano plates [27].Following suit, Tocci Monaco et al. conducted an in-depth investigation into the critical temperatures associated with magneto-electro-elastic nonlocal strain gradient plates [28].Luciano et al. contributed significantly by exploring free flexural vibrations in nanobeams with non-classical boundary conditions [29].Additionally, Fabbrocino et al. made substantial strides in the field with their research on dynamic crack growth utilizing a moving mesh method.These research endeavors collectively offer invaluable insights and research findings, advancing the knowledge base in these areas of study [30].
In recent findings by Wu et al. on dynamics of nanostructures including vibrations in graphene sheets and nanobeams, this work validates computational models for accurate prediction of natural frequencies across a wide parameter space [47,48].
B. Karami et al. investigated the free vibration behavior of a curved microbeam composed of a functionally graded material (FGM) using the modified strain gradient theory of elasticity and firstorder shear deformation theory.In the study, various homogenization models, such as Voigt, Reuss, Hashin-Shtrikman bounds and cubic local representative volume elements (LRVE) schemes were employed to determine the effective material properties of the FGM curved microbeam [49].
Compared to earlier investigations using simple models like Euler-Bernoulli beam theory, current research utilizes more advanced approaches like higher-order shear deformation theory and isogeometric finite element formulation to capture complex behaviors of functional nanobeams.Advantages include analyzing complex geometries, deflections, porosities and foundations for practicable insight.Challenges remain around computational efficiency for real applications.Overall trends point to a focus on novel materials, geometries and theories to progress functional nanobeam research toward utilization through solving practical engineering problems.
This research uses the Mori-Tanaka, Tamura, Reuss and Voigt homogenization schemes and the high-order shear deformation nano-beam theory to derive the governing equations of motion using Hamilton's principle and the nonlocal differential constitutive relations of Eringen.Analytical solutions are then obtained from this data and compared with other solutions in the literature.The effects of the various homogenization schemes, nonlocal parameter, beam aspect ratio and material compositions on the dynamic response of the FG nano-beam are also investigated.This study undertakes a first-of-its-kind comparative vibration analysis of FGM nano-beams using four homogenization techniques and comprehensive parametric study of nonlocal parameters, gradient indices, volume fractions and geometries.This work provides new insights into the natural frequencies and mode shapes.The comparative homogenization scheme analysis and extensive vibrational behavior investigation pave the way for optimized design and future stability, buckling and wave propagation studies of FGM nano-beams.The volume distribution fraction through the thickness has been identified as the following function, using Eq 1, where p is the material index.

Nano-beam with different homogenization schemes
The micromechanics models chosen for the comparison study are detailed in sections 2.1.to 2.4.

Mixture law (Voigt model)
Mixture law, also known as the Voigt model, is a mathematical model that describes the behavior of a composite material.It assumes that the properties of the composite material are obtained by combining the properties of the individual components in a linear manner.It is used to calculate the effective modulus of elasticity, strength and other properties of a composite material using Eq 2 [50].

Reuss model
The Reuss model is a mathematical model used to calculate the effective properties of a composite material.It assumes that the properties of the composite material are obtained by averaging the properties of the individual components.It is used to calculate the effective modulus of elasticity, strength and other properties of a composite material.Eq 3 calculates the proprieties of FG structures with the assumption that the stress is uniform through the thickness [51].

Tamura model
The Tamura model is a mathematical model used to calculate the effective properties of a composite material.It is based on the concept of strain energy density and assumes that the properties of the composite material are obtained by combining the properties of the individual components in a non-linear manner.It is used to calculate the effective modulus of elasticity, strength and other properties of a composite material using Eq 4 [52].
The method of Tamura is another way to express the linear law of Voigt where the empirical term q "stress-to-strain transfer" has been added in formulation [53].

Mori-Tanaka model
The Mori-Tanaka method is a homogenization scheme used to estimate the effective properties of FGM.This method was developed by Mori and Tanaka in 1973 and is based on the assumption that the material properties in the microstructure can be represented by a representative volume element (RVE).The method uses small strain theory and assumes that the RVE is composed of an assembly of small particles of different materials which interact through elasticity.The Mori-Tanaka method can also be used to calculate the effective thermal conductivity and effective electrical conductivity of the FGM using Eq 5 [54].
The Mori-Tanaka homogenization scheme formulation is given by [55][56][57][58]: where P(z) is the effective material property.Pm and Pc are the properties of the metal and ceramic faces of the beam respectively.

Kinematics
The displacement field of the current theory is based on the assumption that in-plane and transverse displacements are split into bending and shear components, that bending parts of in-plane displacement are similar to d'Euler Bernoulli theory and that shear parts of in-plane displacement induce hyperbolic variations of shear strains and shear stresses which vanish on the top and bottom surfaces of the beam.Those result in the following displacement field, as shown in Eq 6: Here,  0 (, ) denotes the displacements along the  coordinate directions of a point on the mid-plane of the beam;   (, )and   (, ) are the bending and shear components of the transverse displacements, respectively.The nonzero strains associated with the displacement field are shown in Eq 7: where:

Constitutive relations
The nonlocal constitutive relation for the macroscopic stress takes the following form, as shown in Eq 8: where E and G are Young's modulus and shear modulus, respectively.Moreover,  0 is a material constant, and  is the internal characteristic length.Once the nonlocal parameter  = ( 0 ) 2 is equal to zero, we obtain the constitutive relations of the local theories.

Equations of motion
Hamilton's principle is used herein to derive the equations of motion.We have obtained: Integrating by parts Eq 9, and collecting the coefficients of  0 ,   and   , the following equations of stability are obtained by Eq 10: where ,   ,   and  are stress resultants and they are defined by: (,   ,   ) = ∫  (1, , )   (11) And ( 0 ,  2 ) are mass inertias defined as: Substituting Eqs 11-13 into Eq 8 and integrating through the thickness of the beam, the stress resultants are related to the generalized displacements by the relations:

Exact solution for a simply-supported FG beam
Following the Navier solution procedure, we assume the solution form for  0 ,   ,  that satisfies simply-supported FG conditions.
where   ,  and   are arbitrary parameters to be determined,  is the eigen-frequency associated with (n) th eigen-mode.
Substituting Eq 16 into Eq 15, the analytical solutions can be obtained from the Eq 17:

Numerical results and discussion
This section compares the proposed homogenization model with the one published by Zemri et al. [59] in order to verify the accuracy of different homogenization models (Reuss, Tamura and Mori-Tanaka).It also provides a parametric analysis of the effects of the nonlocal parameter, the material index, the length-to-thickness ratio and the type of phases constituting the FGM nano-beam on its natural frequency.

Comparisons and validation
The Tanaka, Tamura, and Reuss model was compared to and validated with the Voigt model used by Zemri et al. [59] to ensure the accuracy of the results.The material properties of the FGM nanobeams, consisting of six different combinations of stainless steel (SUS304) and aluminum (Al) with either zirconia (ZrO2), aluminia (Al2O3) or silicon carbide (SiC), are shown in Table 1.These combinations are named FGM1-1 SUS304/ZrO2, FGM1-2 SUS304/Al2O3, FGM1-3 SUS304/SiC, FGM2-1 Al/ZrO2, FGM2-2 Al/Al2O3 and FGM2-3 Al/SiC FGM nano-beams, respectively.In all cases, we present the non-dimensional frequency defined as shown in Eq 18: Tables 2 and 3 provide a comprehensive comparison of the fundamental frequency of nano-beams for the four different models: Voigt, Tanka, Tamura and Reuss.The varying length-to-thickness ratios and material indices of the nonlocal parameter indicate that the results are reliable and accurate.The validation results show good agreement across the different models, demonstrating the trustworthiness of the models.

Effect of the length-to-thickness ratio the FGM nano beam on its natural frequency
Figure 2 demonstrates the variation in the natural frequency of an FGM nano beam when evaluated using four homogenization models (Voigt, Tamura, Tanak and Reuss) with a length-tothickness ratio of 10 to 100, material index of 1 and a nonlocal parameter of 4.
It can be seen from Figure 2 that the natural frequency of the FGM nano-beam increases with the length-to-thickness ratio l for all four homogenization models.As the l/h ratio increases, the beam becomes longer and thinner, resulting in an increase in natural frequency.It is also notable that the Reuss and Tamura models provide the same natural frequency values, while the Voigt and Tanak models give lower values.This discrepancy is due to the assumptions made in each model.Furthermore, the relative difference between the natural frequency evaluated using the four homogenization models decreases as the length-to-thickness ratio increases, albeit only slightly The relative difference = [ω(Reuss) -ω(Voigt, Tamura, Tanak)]/ω(Reuss) between the natural frequency of the FGM nano-beam evaluated using the four homogenization models is 0.059 when l is 10 and 0.038 when l is 100.

Effect of the nonlocal parameter the FGM nano-beam on its natural frequency
Figure 3 shows how the natural frequency of an FGM nano-beam changes as the nonlocal parameter ranges from 0 to 5, while keeping the material index at 1 and length-to-thickness ratio at 10, when evaluated using four different homogenization models (Voigt, Tamura, Tanak and Reuss).
The data in Figure 3 displays that, as the nonlocal parameter increases from 0 to 5, the natural frequency of an FGM nano-beam decreases and the Tanak model gives the lowest frequency value.Among the models, the Voigt model shows the most pronounced decrease in natural frequency, dropping from 6.881 to 5.630 (a decrease of 1.251).The Tanak model shows a decrease of 1.231 from 6.625 to 5.394.The Tamura and Reuss models show a decrease of 0.707 from 7.004 to 6.297.This decrease in frequency is due to the nonlocal parameter making the beam stiffer, which in turn reduces its natural frequency.The larger the nonlocal parameter, the more rigid the beam becomes and thus the lower the frequency.7.2.3.Effect of the material index of the FGM nano-beam on its natural frequency Figure 4 presents the natural frequency of an FGM nano-beam evaluated with four homogenization models (Voigt, Tamura, Tanak and Reuss) when the material index p is changed from 0 to 30, and length-to-thickness ratio and nonlocal parameter are both kept at 10 and 4, respectively.
The natural frequency decreases as the material index increases in Figure 4.The Tanak model produces the smallest values of natural frequency, which means that the nano-beam is the most flexible when evaluated with this model.

Effect of type of phases constituting the FGM nano-beam on its natural frequency
In this section, the impact of various metal and ceramic material combinations on the natural frequency of an FGM nano-beam is studied through a parametric analysis.The figure depicts the natural frequency of FGM nano-beams composed of six different combinations of stainless steel (SUS304) and aluminum (Al) with either zirconia (ZrO2), aluminia (Al2O3) or silicon carbide (SiC).The frequency is calculated using the Reuss model with different l/h ratios, ranging from 10 to 100, a material index of 1 and a nonlocal parameter of 4.
As shown in Figure 5, the length-to-thickness ratio (l/h) has a crucial impact on the natural frequency of FGM nano-beams.As l/h increases, the natural frequency of the FGM nano-beams also increases.In this figure, FGM12 SUS304/Al2O3 has the highest natural frequency among the FGM nano-beams due to its high elastic modulus and low density compared to other FGM nano-beams.Conversely, FGM2-3 Al/SiC has the lowest natural frequency due to its low elastic modulus and high density.It is evident that the natural frequency of FGM nano-beams is greatly influenced by both their length-to-thickness ratio (l/h) and material properties.

Conclusions
This comprehensive parametric study and validation of four homogenization models advances our understanding of predicting the natural frequencies of functionally graded nano-beams.The results for varying length-to-thickness ratios, material indices and nonlocal parameters show excellent agreement across the Voigt, Tanka, Tamura and Reuss models.The relative difference in predicted natural frequencies between models decreases slightly as the length-to-thickness ratio increases.Furthermore, there are significant effects of the nonlocal parameter, material index and constitutive phases on the natural frequency.These findings provide justification that the validated computational homogenization techniques presented can reliably determine the natural frequencies of FGM nanobeams across a range of geometric and material parameters.Researchers can confidently utilize these models to gain new insights into the dynamic characteristics of nano-beams.The comprehensive parametric study and model validation represents an important advancement to the field by establishing these homogenization techniques as trusted tools for the design and analysis of nano-beams.This work elucidates the complex interplay between nano-beam dimensions, materials and intrinsic size effects, moving nano-scale research forward.

Figure 1 shows
Figure 1 shows FG nano-beam of length L, width b and thickness h.The material on the top surface (z = +h/2) of the plate is ceramic and is graded to metal at the bottom surface of the plate (z = −h/2).

Figure 1 .
Figure 1.The geometry domain of the FG nano-beam.

Table 1 .
Material properties of metal and ceramic.

Table 2 .
Fundamental frequency of nano-beams for different scheme's models with varying length-to-thickness ratios, material index and nonlocal parameters.

Table 3 .
Fundamental frequency of nano-beams for different scheme's models with varying length-to-thickness ratios, material index and nonlocal parameters.