Vibration Characteristics of Magnetostrictive Composite Cantilever Resonator with Nonlocal Effect

Taking the nonlocal effect into account, the vibration governing differential equation and boundary conditions of a magnetostrictive composite cantilever resonator were established based on the Euler magnetoelastic beam theory. The frequency equation and vibration mode function of the composite cantilever were obtained by means of the separation of variables method and the analytic solution of ordinary differential equations. The lateral deflection, vibration governing equations, and boundary conditions were nondimensionalized. Furthermore, the natural frequency and modal function of the composite beam were quantitatively analyzed with different nonlocal parameters and transverse geometry dimensions using numerical examples. Compared with the results without considering the nonlocal effect, the influence of the nonlocal effect on the vibration characteristics was analyzed. The numerical results show that the frequency shift and frequency band narrowing of the magnetostrictive cantilever resonator are induced by nonlocal effects. In particular, the high-frequency vibration characteristics, such as vibration amplitude and modal node of the composite beam, are significantly affected. These analysis results can provide a reference for the functional design and optimization of magnetostrictive resonators.


Introduction
Smart materials are materials that have both intrinsic and extrinsic capabilities to respond to stimuli and environmental changes and activate their functions in response to these changes.Magnetostrictive materials with good magneto-mechanical coupling properties undergo dimensional changes when exposed to a magnetic field, and conversely, they can also effectively transform the alternating magnetic energy in the environment into mechanical vibration through this magnetostrictive effect.As an important component device in micro-electro-mechanical systems (MEMS), magnetostrictive resonators are widely used in actuators, sensors, energy collectors, frequency controllers, and logic switches [1][2][3][4].High-performance MEMS require the resonator to have the advantages of high sensitivity, high precision, and low energy consumption when operating.Due to the complexity and miniaturization of resonators, higher requirements are put forward for their vibration performance.For example, the operating bandwidth is wider, the energy density is higher, and the driving voltage is lower [5][6][7][8].
Therefore, it is important to accurately predict the vibration characteristics of microsized magnetostrictive resonators to improve the design and use of resonator devices.
Magnetostrictive resonators are divided in structure into two major categories, rodtype and beam-type, which are mostly made into cantilevers.For the longitudinal vibration Sensors 2024, 24, 5390 2 of 14 of rod-type resonators, relative experiments, numerical simulations, and theoretical analyses are available in the references [9][10][11].Recently, beam-type magnetostrictive resonators with layered composite structures have attracted much attention because of their multiinterface physical field effects, such as magnetomechanics and magnetoelectric.Yu et al. [12] used Kirchhoff theory to simulate the surface effect by applying the spring force on the boundary and analyzed the influence of electromechanical resonance characteristics of the nanocantilever beam device in MEMSs.Jia et al. [13] proposed a nonlinear magnetomechanical coupling model at low magnetic fields using Gibbs free energy method, and the predicted results were in good agreement with the experimental data under given prestress conditions.Reddy et al. [14] obtained an analytical solution to the vibration problem of giant magnetostrictive laminated simply supported beams based on Hamilton's principle.Xu and Shang [15] used high-order shear deformation theory and linear coupled magnetoelastic constitutive relations to construct the vibration governing equation of a magnetostrictive composite cantilever actuator and presented the influence of geometrical, material parameters, and excitation frequency on the structure's magneto-mechanical coupling coefficient in detail.Xu and Shang [16] applied the nonlinear constitutive relation of magnetostrictive materials to obtain the vibration differential equation of a double-layer cantilever beam using Hamilton's principle and analyzed frequency-doubling dynamic characteristics of the amplitude-frequency curve and deflection response under periodic magnetic field excitation.At present, the theoretical research of magnetostrictive resonator model building is almost entirely based on classical beam or plate theory, but the mathematical model building of micro-devices rarely considers or ignores the influence of microscale on the structural vibration characteristics.
As early as 1994, Quandt et al. [17] studied and designed a double-layer thin-film cantilever-type magnetostrictive actuator and obtained the relationship curve between deflection and magnetic field strength through experiments.In recent years, researchers [8,18,19] designed magnetoelectric composite resonators with different structural forms and observed the effects of frequency doubling, mixing, and resonance frequency shifts of resonators under a magnetic field through experiments.Using the finite element method, the sensing characteristics of the cantilever resonator under different vibration modes [20] and the design optimization of the device were analyzed [21], as well as the influence of the magnetic field on the resonant frequency [22].In addition, magnetostrictive resonators have also become a research focus of scholars in biosensing technology due to their advantages in wireless detection, simple operation, high sensitivity, and fast detection speed [23,24].
Experimental results showed the size-dependent effect of the magnetostrictive material on a small scale (i.e., micro/nano-scale), which plays an important role in designing miniaturized smart devices [25,26].Therefore, it is necessary to understand the vibration behavior of magnetostrictive micro resonators by considering microscale theoretical models.So far, continuum mechanics is commonly used to describe the size-dependent effect in microstructures [27,28].Based on previous studies [15,16], this paper introduced microscale parameters into the material constitutive equation via nonlocal elasticity theory and established a vibration control model of a magnetostrictive composite cantilever beam.The effect of microscale parameters on the vibration characteristics of the resonator is further studied by numerical calculation.

Theoretical Foundation and Modeling
We consider the vibration of a rectangular double-layer composite cantilever beam generated by the uniformly distributed horizontal magnetic field H, as depicted in Figure 1.The length of the beam is L, the width is b, and the thickness is h.Layer one and layer two refer to the magnetostrictive layer and the elastic substrate layer, respectively.A Cartesian coordinate system Oxz is attached to the beam with z along the beam thickness direction and x-axis on the interface of the magnetostrictive layer and the elastic substrate layer.The elastic modulus is denoted as E i , the shear elastic modulus as G i , Poisson's ratio as υ i , the density as ρ i , the thickness as h i , the cross-sectional area A i , and the moment of inertia to the y i axis as I yi = bh 3 i /12 for the i-th layer of the composite cantilever beam (i = 1, 2), where i = 1 and i = 2 represent the magnetostrictive layer and the base layer, respectively.as υi, the density as ρi, the thickness as hi, the cross-sectional area Ai, and the moment of inertia to the yi axis as  Different from classical elasticity theory, the nonlocal continuum theory is capable of capturing the size-dependent effect in magnetostrictive microstructures.In nonlocal theory, the stress at an arbitrary point of the magnetostrictive microstructure is a function of the strain and magnetic field at all points in the structure.Thus, based on the nonlocal elasticity theory of Eringen [29], the magnetoelastic constitutive relation of a one-dimensional composite beam can be written in the following form: where ε and ( ) H t are the stress, strain, and magnetic field, respectively.It is emphasized that ea is the nonlocal parameter, which represents the influence of small scale on the behavior of microstructure, where e is a given material constant and can be determined by molecular dynamics simulations, experiments, or other computational methods [30] and a depends on the internal microstructure size of the material such as C-C bond length, lattice spacing, and particle distance.Furthermore, d is the magnetostrictive coefficient, as the elastic substrate layer is not influenced by the magnetic field, where the parameter d2 = 0. Based on the Euler-Bernoulli beam theory, the displacements of the composite beam can be assumed to be as follows: ( ) ( , , ) ( , ) where u(x, t) and w(x, t) are the longitudinal displacement and the deflection at the interface (z = 0), respectively.
The axial force, shear force and bending moment of the layer ith cross section are denoted by M M M = + are the total transverse shear force and the total bending moment acting on the whole cross section, and the torque of the axial force on the y-axis is omitted.
Ignoring the horizontal inertial force and inertial couple, we can obtain the internal force dynamic balance equation: Different from classical elasticity theory, the nonlocal continuum theory is capable of capturing the size-dependent effect in magnetostrictive microstructures.In nonlocal theory, the stress at an arbitrary point of the magnetostrictive microstructure is a function of the strain and magnetic field at all points in the structure.Thus, based on the nonlocal elasticity theory of Eringen [29], the magnetoelastic constitutive relation of a one-dimensional composite beam can be written in the following form: where σ x , ε x and H(t) are the stress, strain, and magnetic field, respectively.It is emphasized that ea is the nonlocal parameter, which represents the influence of small scale on the behavior of microstructure, where e is a given material constant and can be determined by molecular dynamics simulations, experiments, or other computational methods [30] and a depends on the internal microstructure size of the material such as C-C bond length, lattice spacing, and particle distance.Furthermore, d is the magnetostrictive coefficient, as the elastic substrate layer is not influenced by the magnetic field, where the parameter d 2 = 0. Based on the Euler-Bernoulli beam theory, the displacements of the composite beam can be assumed to be as follows: where u(x, t) and w(x, t) are the longitudinal displacement and the deflection at the interface (z = 0), respectively.The axial force, shear force and bending moment of the layer ith cross section are denoted by S and M (i) (i = 1, 2), respectively.
N and M = M (1) + M (2)  are the total transverse shear force and the total bending moment acting on the whole cross section, and the torque of the axial force on the y-axis is omitted.Ignoring the horizontal inertial force and inertial couple, we can obtain the internal force dynamic balance equation: It is assumed that the nonlocal length parameters of the upper and lower layers of the composite beam are the same, namely µ = (e 1 a 1 ) 2 = (e 2 a 2 ) 2 .Substituting the geometric relation ε x /∂x into Equation (1), and using Equation (2), the axial stress of each layer can be expressed as follows: (1) Integrate both sides of Equations ( 6) and ( 7) with respect to z in the interval [0, h 1 ] and [−h 2 , 0], using the axial force definition F (1) x dz and F (2) x dz, we obtain (1) Similarly, multiply both sides of Equations ( 6) and ( 7) by z, and integrate with respect to z in the interval [0, h 1 ] and [−h 2 , 0].Using the bending moment definition x dz and the equilibrium Equations ( 4) and ( 5), the total bending moment and total shear force can be obtained: The axial force expression ( 8) and ( 9) are derived with respect to x and then added together; using the equilibrium Equation (3), we obtain: Substituting Equation ( 11) into equilibrium Equation (4) yields: Taking the derivative of both sides of Equation ( 12) with respect to x and substituting it into Equation ( 13), we obtain the vibration control equation of the magnetostrictive composite cantilever beam considering nonlocal effects: where EI is the equivalent flexural stiffness, ρA is the equivalent inertial mass, and EA is equivalent tensile stiffness of the composite beam; their expressions are as follows: It should be noted that in the expression of the equivalent flexural stiffness EI, in addition to the direct superposition of the flexural stiffness of the two layers of the beam, there is also an additional contribution from the equivalent tensile stiffness EA.This is Sensors 2024, 24, 5390 5 of 14 because the total deformation of the double-layer composite beam under the action of the magnetic field is the coupling deformation of axial stretching and transverse bending, which accords with the deformation mechanism of the magnetostrictive cantilever beam resonator.In addition, the nonlocal parameter µ does not affect the stiffness and inertial mass of the beam [27].
At the fixed end of the cantilever beam (x = 0, −h 2 ≤ z ≤ h 1 ), the following boundary condition of zero displacement should be satisfied: At the free-end (x = L), the boundary condition for free stress is usually not exactly satisfied; since we propose the Saint-Venant boundary condition, i.e., the internal force is zero, we have: By integrating both sides of Equation ( 12) with respect to x, and using the fixed end boundary conditions (15) and the free end boundary conditions ( 16), the axial displacement of the interface can be written as follows: Substituting Equation ( 17) into Equations ( 10) and ( 11), the bending moment and shear force of the composite beam section can be expressed as follows: From Equations ( 17)- (19), we obtain the boundary conditions of the magnetostrictive composite cantilever beam as follows: For the practical application of the resonators, it is usually necessary to consider the steady-state periodic vibration of cantilever beam elements caused by periodic magnetic field excitation, so that their displacement response satisfies the following periodic conditions: where 0 ≤ x ≤ L, T = 2π/ω and ω are the vibration period and relative circular frequency respectively.Up to now, the vibration problem of the magnetostrictive composite cantilever resonator considering the nonlocal effect is mathematically reduced to solving the governing Equation ( 14) associated with the boundary conditions (20) and periodic conditions (21).Different from Euler's classical beam theory, the vibration equation and boundary condi-tions contain the nonlocal parameter µ that describes nonlocal effects, and when µ = 0, the vibration governing equation degenerates to the classical Euler beam vibration governing equation without considering nonlocal effects.It is worth noting that the vibrational excitation term does not appear in the vibration equation, but in the boundary conditions.

Solution to the Vibration Problem
The magnetic field H is assumed to be a harmonic function of time t; that is, H(t) = H 0 e iωt (i = √ −1).Denote the dimensionless variables and parameters as follows: We can rewrite the governing Equation ( 14) as follows: Accordingly, the dimensionless form of the definite solution condition (20) can be rewritten as follows: where H(τ) = H 0 e i2πτ .The dimensionless form of the periodic condition ( 21) is as follows: Based on the periodic condition (25), the magnetoelastic vibration response of the composite beam can be expressed as the following separate variable form: where W(ξ) is the displacement amplitude.Substituting Equation (26) into Equation (23) and eliminating the time variable τ, the following ordinary differential equation with constant coefficients is obtained: The general solution of Equation ( 27) is given in the form: where Substituting the solution expression Equation ( 28) into the boundary Equations (24a)-(24c), we obtain the system of linear algebraic equations for determining the coefficients B j (j = 1, 2, 3, 4) as follows: Sensors 2024, 24, 5390 Denote the coefficient matrix of (29) as A(ω).
For the case of excited vibration H(τ) ̸ = 0, that is, the closed circuit boundary conditions, we consider the case of periodic vibration: the magnetostrictive layer is subjected to a uniform magnetic field of sinusoidal variation with time, that is H(τ) = H 0 sin 2πτ.By substituting the known excitation frequency ω into Equation ( 29), the coefficient B j can be determined, and then using Equation ( 26), we finally obtain the vibration response of the composite beam.
For the case of free vibration, the external excited magnetic field H(τ) ≡ 0, that is the open-circuit boundary conditions, the magnetic induction on the surfaces of the beam is zero; hence the boundary condition (29) becomes homogeneous.The condition B j ̸ = 0 from Equation (29) gives det(A(ω)) = 0. Thus, the natural frequency equation of the composite cantilever beam can be derived as follows: Substituting the expression of λ i (ω) (i = 1, 2) into Equation ( 30), yields The natural frequencies of each order B 2 j = 1, the co- efficients B j (j = 1, 2, 3, 4) corresponding to λ i (ω) can be obtained from Equation (29).Consequently, the vibration mode function W n (ξ) is given by Equation (28).

Numerical Examples and Discussion
We now apply our solution to a double-layer beam composed of the magnetostrictive (Tb 0.3 Dy 0.7 ) 0.42 Fe 0.58 and silicon to investigate the influence of the nonlocal parameter on natural frequency and mode shapes.Take the magnetostrictive composite cantilever beam device mentioned in reference [17] as an example, which has a length L = 20 mm, width b = 5 mm, and a total thickness h = 60 µm.For the magnetostriction layer, the thickness h 1 = 10 µm, the elastic modulus E 1 = 50 Gpa, the density ρ 1 = 9.25 × 10 3 kg/m 3 , the Poisson's ratio ν 1 = 0.3, and the magnetostriction coefficient d 1 = 5.027 × 10 −8 m/A.For the substrate layer, the thickness h 2 = 10 µm, the elastic modulus E 2 = 130 Gpa, the density ρ 2 = 2.33 × 10 3 kg/m 3 and Poisson's ratio ν 2 = 0.28.
By referring to the value of the dimensionless nonlocal parameter l in reference [27], we consider the parameter l to be 0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.35 and 0.4 in this paper.The first three natural frequencies for the composite beam with different parameters l calculated by Equation (30) are shown in Table 1 and compared with the results without nonlocal effects (when l = 0, the nonlocal composite beam is reduced to the classical one).It is found that, with increasing nonlocal parameter l and the order of modes, the natural frequency differences between, with, and without nonlocal effects increase monotonically.For example, when the parameter l = 0 and l = 0.2, the relative deviations of natural frequencies are 1.86% for f 1 , 6.14% for f 2 and 40.3% for f 3 .Obviously, the nonlocal effect has a greater impact on the higher order frequencies than the lower-order ones.In the practical application of micro/nano electromechanical systems, magnetostrictive resonators are more likely to exhibit high-frequency vibration characteristics, so the influence of nonlocal effects on their frequency cannot not be ignored.Note: Data marked with asterisk '*' are the results of l = 0, which are consistent with the natural frequency results of classic Euler beam theory in reference [16].
Figure 2 shows the variation curves of the first three orders of natural frequencies for the composite beam with the dimensionless nonlocal parameter l.The fundamental frequency is the first natural frequency of the structure.For sensors and actuators, the fundamental frequency is the most easily excited vibration frequency of the structure, and so it is one of the important vibration indicators.It is interesting that with the increasing parameter l, the fundamental frequency of the composite beam increases, while the second and third natural frequencies decrease.That is, the frequency bandwidths of the composite beam become narrower, and the narrowing speed accelerates with the increase in the parameter l.This demonstrates that the nonlocal effect weakens the effective structural stiffness of the composite beam and thus causes the frequency shift effect, which affects the vibration characteristics of the resonator.Note: Data marked with asterisk '*' are the results of l = 0, which are consistent with the natural frequency results of classic Euler beam theory in reference [16].
Figure 2 shows the variation curves of the first three orders of natural frequencies for the composite beam with the dimensionless nonlocal parameter l.The fundamental frequency is the first natural frequency of the structure.For sensors and actuators, the fundamental frequency is the most easily excited vibration frequency of the structure, and so it is one of the important vibration indicators.It is interesting that with the increasing parameter l, the fundamental frequency of the composite beam increases, while the second and third natural frequencies decrease.That is, the frequency bandwidths of the composite beam become narrower, and the narrowing speed accelerates with the increase in the parameter l.This demonstrates that the nonlocal effect weakens the effective structural stiffness of the composite beam and thus causes the frequency shift effect, which affects the vibration characteristics of the resonator.
Figure 3 gives the first three orders of mode shape curves of the composite beam with different dimensionless nonlocal parameters l.It is observed that the mode shapes roughly follow three groups.The solid line of each group curve represents the mode shape when the parameter l = 0, and the dashed line represents the corresponding mode curve when the parameter l is 0.1, 0.2 and 0.3.It can be seen that for the same-order mode, the larger the nonlocal parameter l, the greater the change in beam amplitude and the more obvious the vibration effect.It is worth noting that nonlocal effects could also cause changes in the positions where the amplitude of the beam is zero (i.e., the modal node).Taking the second-order mode (blue curves) as an example, the modal node position moved from 0 .7 5 ξ = for l = 0 towards the free end to about 0 .8ξ = for l = 0.3.Therefore, in the practical application of resonators, the position errors of modal nodes and antinodes (i.e., where amplitude is at its peak) would affect the function of high-performance devices.In other words, both modal node and anti-node positions could be tuned by the nonlocal parameter, which could be very beneficial to the functional design of composite micro beam devices.Figure 3 gives the first three orders of mode shape curves of the composite beam with different dimensionless nonlocal parameters l.It is observed that the mode shapes roughly follow three groups.The solid line of each group curve represents the mode shape when the parameter l = 0, and the dashed line represents the corresponding mode curve when the parameter l is 0.1, 0.2 and 0.3.It can be seen that for the same-order mode, the larger the nonlocal parameter l, the greater the change in beam amplitude and the more obvious the vibration effect.It is worth noting that nonlocal effects could also cause changes in the positions where the amplitude of the beam is zero (i.e., the modal node).Taking the second-order mode (blue curves) as an example, the modal node position moved from ξ = 0.75 for l = 0 towards the free end to about ξ = 0.8 for l = 0.3.Therefore, in the practical application of resonators, the position errors of modal nodes and antinodes (i.e., where amplitude is at its peak) would affect the function of high-performance devices.In other words, both modal node and anti-node positions could be tuned by the nonlocal parameter, which could be very beneficial to the functional design of composite micro beam devices.In order to further demonstrate the effect of the nonlocal parameter on the natural frequency of a composite cantilever beam with different transverse dimensions, the first three orders of natural frequencies with different layer thickness ratios c ( ) and nonlocal parameters l are listed in Table 2 for a fixed total thickness of the composite beam h.Here, not only is the natural frequency determined, but the shift in the natural frequency is also studied.We introduce the frequency difference and the relative frequency shift, noted as ), where j f , 0j f are frequen- cies with and without nonlocal effects, respectively.It is observed clearly that, as the thickness of the magnetostrictive layer increases (that is, c increases), the natural frequencies decrease, and the frequency difference j f Δ also decreases, but the relative frequency shift j R almost does not change with the thickness ratio c.
Figure 4 shows the variation curves of the natural frequency of the composite cantilever beam with thickness ratio c when the nonlocal parameter l = 0, 0.1, 0.2.It can be seen that, when the magnetostrictive layer is thicker (the larger c is), the overall stiffness of the composite beam is smaller, and the natural frequency of the structure is lower.In addition, it is clear that the difference in frequency shift caused by nonlocal effects is larger for higher vibration modes.
Figures 5 and 6 show the first three natural frequency difference and relative frequency shift curves of the composite beam with different transverse dimensions as a function of nonlocal parameters.With the decrease in the interlayer thickness ratio and the increase in the microscale parameters, the frequency difference j f Δ with and without nonlocal effects presents a more nonlinear difference, while the thickness ratio has little influence on the relative frequency shift j R for the same mode and almost no change.This feature indicates that, for a given composite beam thickness, a slight adjustment to the nonlocal parameter l would not influence the relative frequency shift of the structure, regardless of the layer thickness of the two materials.In order to further demonstrate the effect of the nonlocal parameter on the natural frequency of a composite cantilever beam with different transverse dimensions, the first three orders of natural frequencies with different layer thickness ratios c (c = h 1 /h) and nonlocal parameters l are listed in Table 2 for a fixed total thickness of the composite beam h.Here, not only is the natural frequency determined, but the shift in the natural frequency is also studied.We introduce the frequency difference and the relative frequency shift, noted as ∆ f j = f j − f 0j and R j = ∆ f j / f 0j • 100% (j = 1, 2, 3), where f j , f 0j are frequencies with and without nonlocal effects, respectively.It is observed clearly that, as the thickness of the magnetostrictive layer increases (that is, c increases), the natural frequencies decrease, and the frequency difference ∆ f j also decreases, but the relative frequency shift R j almost does not change with the thickness ratio c.
Figure 4 shows the variation curves of the natural frequency of the composite cantilever beam with thickness ratio c when the nonlocal parameter l = 0, 0.1, 0.2.It can be seen that, when the magnetostrictive layer is thicker (the larger c is), the overall stiffness of the composite beam is smaller, and the natural frequency of the structure is lower.In addition, it is clear that the difference in frequency shift caused by nonlocal effects is larger for higher vibration modes.
Figures 5 and 6 show the first three natural frequency difference and relative frequency shift curves of the composite beam with different transverse dimensions as a function of nonlocal parameters.With the decrease in the interlayer thickness ratio and the increase in the microscale parameters, the frequency difference ∆ f j with and without nonlocal effects presents a more nonlinear difference, while the thickness ratio has little influence on the relative frequency shift R j for the same mode and almost no change.This feature indicates that, for a given composite beam thickness, a slight adjustment to the nonlocal parameter l would not influence the relative frequency shift of the structure, regardless of the layer thickness of the two materials.In order to further demonstrate the effect of the nonlocal parameter on the output of the magnetostrictive resonator, Figure 7 shows the deflection response curves (at one period T) on the free end of the magnetostrictive thin-film cantilever device [17] with different nonlocal parameters subject to the external excitation frequency H0 = 500 Hz.The curves of the solid line group represent the magnetostrictive coefficient   In order to further demonstrate the effect of the nonlocal parameter on the output of the magnetostrictive resonator, Figure 7 shows the deflection response curves (at one period T) on the free end of the magnetostrictive thin-film cantilever device [17] with different nonlocal parameters subject to the external excitation frequency H0 = 500 Hz.The curves of the solid line group represent the magnetostrictive coefficient  In order to further demonstrate the effect of the nonlocal parameter on the output of the magnetostrictive resonator, Figure 7 shows the deflection response curves (at one period T) on the free end of the magnetostrictive thin-film cantilever device [17] with different nonlocal parameters subject to the external excitation frequency H 0 = 500 Hz.The curves of the solid line group represent the magnetostrictive coefficient d 1 = 5.027 × 10 −8 m/A, and the curves of the dashed line group represent the magnetostrictive coefficient d 2 = 8.437 × 10 −8 m/A in Figure 7.It can be seen from the figure that the magnetostrictive coefficient d has a significant influence on the output of deflection for the free end, while the nonlocal parameter l has a slight effect on the deflection of the free end.

Conclusions
The vibration problem of the magnetostrictive composite cantilever beam resonator is investigated using nonlocal elastic and Euler beam theory.The frequency equation and vibration mode function of the composite magnetostrictive cantilever are obtained by means of the separation of variables method and the analytic solution of ordinary differential equations.The effects of nonlocal parameters, interlayer thickness, and magnetostriction coefficient on the vibration characteristics of the magnetostrictive resonator are investigated numerically.1. Different from the classical Euler beam theory, the vibration equation and boundary conditions of the nonlocal theory contain the parameter µ to describe the nonlocal effects, and when µ = 0, the vibration equation and boundary conditions degenerate to the classical theoretical form (without considering nonlocal effects).2. The coupling deformation of axial stretching and transverse bending occurs in a double-layer composite beam subjected to the magnetic field, and the vibration governing equation is a fourth-order ordinary differential equation.In particular, the equivalent flexural stiffness of the beam includes not only the flexural stiffness but also the tensile stiffness, while the nonlocal parameter does not affect the stiffness and inertial mass of the beam.3. The natural frequencies and vibration modal shapes of magnetostrictive composite beams depend strongly on the nonlocal parameters.Consequently, the nonlocal effects reinforce the vibration characteristics of the micro/nano resonator, especially on the high-order vibration modes.4. Compared to the vibration modes, the nonlocal parameter l has only a slight effect on the deflections (i.e., elastic displacements) of the free end.The significant influence of l on natural frequency and vibration modes in different thickness laminated beams could be utilized in designing better laminated micro/nano magnetostrictive resonators with nonlocal effects.
Author Contributions: Conceptualization, Y.X. and X.S.; methodology, Y.X. and X.S.; formal analysis, Y.X.; investigation, Y.X.; writing-original draft preparation, Y.X.; Writing-review & editing, X.S.; funding acquisition, Y.X. and K.X.All authors have read and agreed to the published version of the manuscript.

Conclusions
The vibration problem of the magnetostrictive composite cantilever beam resonator is investigated using nonlocal elastic and Euler beam theory.The frequency equation and vibration mode function of the composite magnetostrictive cantilever are obtained by means of the separation of variables method and the analytic solution of ordinary differential equations.The effects of nonlocal parameters, interlayer thickness, and magnetostriction coefficient on the vibration characteristics of the magnetostrictive resonator are investigated numerically.

1.
Different from the classical Euler beam theory, the vibration equation and boundary conditions of the nonlocal theory contain the parameter µ to describe the nonlocal effects, and when µ = 0, the vibration equation and boundary conditions degenerate to the classical theoretical form (without considering nonlocal effects).

2.
The coupling deformation of axial stretching and transverse bending occurs in a double-layer composite beam subjected to the magnetic field, and the vibration governing equation is a fourth-order ordinary differential equation.In particular, the equivalent flexural stiffness of the beam includes not only the flexural stiffness but also the tensile stiffness, while the nonlocal parameter does not affect the stiffness and inertial mass of the beam.

3.
The natural frequencies and vibration modal shapes of magnetostrictive composite beams depend strongly on the nonlocal parameters.Consequently, the nonlocal effects reinforce the vibration characteristics of the micro/nano resonator, especially on the high-order vibration modes.4.
Compared to the vibration modes, the nonlocal parameter l has only a slight effect on the deflections (i.e., elastic displacements) of the free end.The significant influence of l on natural frequency and vibration modes in different thickness laminated beams could be utilized in designing better laminated micro/nano magnetostrictive resonators with nonlocal effects.
= / for the i-th layer of the composite cantilever beam (i = 1, 2), where i = 1 and i = 2 represent the magnetostrictive layer and the base layer, respectively.

Figure 1 .
Figure 1.Schematic diagram of geometric dimensions and coordinates for the composite cantilever beam.

Figure 1 .
Figure 1.Schematic diagram of geometric dimensions and coordinates for the composite cantilever beam.

3 Figure 2 .
Figure 2. The first three orders of nature frequencies of the composite cantilever beam vary with the dimensionless nonlocal parameters.

Figure 2 .
Figure 2. The first three orders of nature frequencies of the composite cantilever beam vary with the dimensionless nonlocal parameters.

Figure 3 .
Figure 3. Curves of first three orders of vibration modes of composite beam with different nonlocal parameters (h1/h = 0.167).

Figure 3 .
Figure 3. Curves of first three orders of vibration modes of composite beam with different nonlocal parameters (h 1 /h = 0.167).

Figure 4 .
Figure 4. Nature frequencies f versus interlayer thickness ratios c with different microscale parameters l.

Figure 4 .Figure 5 .
Figure 4. Nature frequencies f versus interlayer thickness ratios c with different microscale parameters l.

Figure 6 .
Figure 6.Relative frequency shifts j R versus nonlocal parameters l with different transverse di-

,in Figure 7 .
and the curves of the dashed line group represent the magnetostric-It can be seen from the figure that the magnetostrictive coefficient d has a significant influence on the output of deflection for the free end, while the nonlocal parameter l has a slight effect on the deflection of the free end.

Figure 5 .Figure 5 .
Figure 5. Frequency differences ∆ f j versus nonlocal parameters l with different transverse dimensions.

Figure 6 .
Figure 6.Relative frequency shifts j R versus nonlocal parameters l with different transverse di-

,in Figure 7 .
and the curves of the dashed line group represent the magnetostric-It can be seen from the figure that the magnetostrictive coefficient d has a significant influence on the output of deflection for the free end, while the nonlocal parameter l has a slight effect on the deflection of the free end.

Figure 6 .
Figure 6.Relative frequency shifts R j versus nonlocal parameters l with different transverse dimensions.

Figure 7 .
Figure 7.The deflection response curves on free end of composite cantilever beam with different magnetostrictive coefficient d varies with microscale parameter l (c = 0.167).

Figure 7 .
Figure 7.The deflection response curves on free end of composite cantilever beam with different magnetostrictive coefficient d varies with microscale parameter l (c = 0.167).

Table 1 .
The first three order nature frequencies of the composite cantilever beam with different dimensionless nonlocal parameters l (Hz).

Table 2 .
First three orders of nature frequencies and relative frequency shift rates of composite beam with different interlayer thickness ratios c.

Table 2 .
First three orders of nature frequencies and relative frequency shift rates of composite beam with different interlayer thickness ratios c.