Modeling the Converse Magnetoelectric Effect in the Low-Frequency Range

This article is devoted to the theory of the converse magnetoelectric (CME) effect for the longitudinal, bending, longitudinal-shear, and torsional resonance modes and its quasi-static regime. In contrast to the direct ME effect (DME), these issues have not been studied in sufficient detail in the literature. However, in a number of cases, in particular in the study of low-frequency ME antennas, the results obtained are of interest. Detailed calculations with examples were carried out for the longitudinal mode on the symmetric and asymmetric structures based on Metglas/PZT (LN); the bending mode was considered for the asymmetric free structure and structure with rigidly fixed left-end Metglas/PZT (LN); the longitudinal-shear and torsional modes were investigated for the symmetric and asymmetric free structures based on Metglas/GaAs. For the identification of the torsion mode, it was suggested to perform an experiment on the ME structure based on Metglas/bimorphic LN. All calculation results are presented in the form of graphs for the CME coefficients.


Introduction
The continuous development of technology requires constant involvement in the process of materials with new multifunctional properties.One of these materials is a composite material with a magnetoelectric (ME) effect [1][2][3].A large number of studies are already known on these materials in connection with the prospects for their use in the development of new electronic devices [4].There is the direct ME effect (DME) and converse ME effect (CME) [5], where the direct effect is characterized by the induction of an output electrical voltage in the ME composite under the influence of a magnetic field, and the converse effect is the excitation of magnetization under the influence of an external electric field.It should be noted that the converse effect is often referred to in the literature as inverse.The DME effect has already been studied in sufficient detail [5,6], and on its basis, a number of promising electronic devices [7][8][9][10] have been proposed.Relatively few works have been devoted to the CME effect.The first studies of the CME effect were associated with the microwave ME effect, which consists of a shift in the ferromagnetic resonance (FMR) line in the ME structure under the influence of an external electric field [11][12][13][14][15] and the further use of this effect to create tunable microwave devices [5,6].In recent years, research interest in the CME effect in the low-frequency region has increased.This is due to the great interest in the design of magnetic field sensors [16][17][18] and low-frequency ME antennas [19][20][21][22] and the possibility of their use in spintronics [23][24][25][26][27][28] and ME memory devices [29][30][31].The purpose of this article is a brief analysis of the works on the CME effect Sensors 2024, 24, 151 2 of 32 and to obtain expressions for CME coefficients for the low-frequency region, including the electromechanical resonance (EMR) modes: longitudinal, bending, longitudinal-shear, and torsional.
Below is a brief overview of the work on the low-frequency CME effect in composites.In one of the first works, Jia et al. [32] theoretically and experimentally studied the CME effect in the Terfenol-D/PMN-PT structure and obtained an ME coefficient of 105 mG/V at a bias field of 170 Oe.A linear relationship was found between the measured magnetic induction and the applied alternating electrical voltage in the range of 50-160 V.In [33][34][35][36][37][38][39][40][41][42][43][44], the authors examined the features of the CME effect in various structures.In [33], the measurements of CME effects were conducted in the trilayers of PZT/Ni/PZT.A pickup coil wound around the sample was used to measure the CME effect due to the change in the magnetic induction in Ni.The measured static magnetic-field dependence of the ME voltage has been attributed to the variation in the piezomagnetic coefficient for Ni.The frequency dependence of the CME effect showed a resonance character due to radial acoustic modes in PZT.Chen et al. [34] investigated the CME effect in a transition metalsbased ferromagnetic-shape memory alloy/piezoelectric ceramic-laminated composite.The strong CME effect was observed at room temperature over a broad bandwidth, under weak magnetic bias and an electric field.The authors of [35] studied the CME effect in a two-layer composite consisting of piezoelectric PZT and ferrite MZF plates.Dependences of the magnetic induction variation Bm on the amplitude and frequency of the AC electric field, as well as on magnetic field strength, temperature, and composition for layered PZT/MZF composite, have been determined.It is shown in [36] that the voltage transformation coefficient depends not only on the CME coefficient and the number of turns of the inductor but also on the relative orientation of the magnetic field and the geometric dimensions of the sample.The CME effect was studied [37] in a layered composite Metglas/PMN-PT by the induction method at various frequencies.A large CME coefficient of 3.05 G/V was observed at a resonant frequency of 76.5 kHz in a weak bias magnetic field of 50 Oe.The CME coefficient of the heterostructure is practically constant and has a relatively high value in a wide frequency range of 1-64 kHz.Xuan et al. [38] investigated the CME effect in Metglas/PZT, and the highest value of the CME coefficient was 4.18 G/V at the resonant frequency of 103 kHz in a low bias field of 45 Oe.In addition, the authors presented data on the known CME coefficients, where the highest value was 18.6 G/V at a frequency of 47 kHz, obtained for a three-layer Ni-Mn-Ga/PMN-PT/Ni-Mn-Ga structure.In [39], using the design of a Rosen-type piezoelectric transducer and a Terfenol-D plate with a Terfenol-D thickness of 0.7 mm, a gigantic CME coefficient of 35.7 G/V was obtained.Lu et al. [40] demonstrated the equivalence of the DME effects and the CME ones in two-phase systems consisting of piezoelectric and magnetostrictive materials.This was achieved by reformulating Maxwell's relation in terms of the effective electric and magnetic dipole moments of the system and comparing the coupling forces at the same electric and magnetic DC biases.In [41], a dispersed composite consisting of 65 mol.% Na 0.5 Bi 0.5 TiO 3 and 35 mol.% CoFe 2 O 4 and its structure, microstructure, ferroelectric, magnetostrictive, magnetic, and direct/converse magnetoelectric properties were studied.The composite showed different magnetization behaviors under electrically polarized and unpolarized conditions.The percentage change in magnetization due to polarization is approximately 15% at a magnetic field of 500 Oe.Bilayer modeling of DME and CME coefficients for longitudinal modes in the EMR region has been considered and presented in [45].As an example, specific cases of bilayers of cobalt ferrite/PZT and nickel/PZT bilayers were studied.In [45], the DME and CME in Metglas/LN and Metglas/PMN-PT trilayers have been studied in the range from 20 Hz to 0.4 MHz.It is shown that, as a result of the difference in piezoelectric constants, the CME in the structure with PMN-PT is significantly greater (~2 orders of magnitude) than with LN.The article by [46] is devoted to the study of the linearity of the CME effect and improving sensitivity in the low-frequency range.Kalgin et al. [43] researched the CME effect in the ME structures consisting of a layer of Terfenol-D ferromagnetic powder, epoxy adhesive, and a polarized piezoelectric layer of PZT.The dependences of the CME effect on the frequency and strength of the electric field, the strength of the constant magnetic field, the thickness of the ferromagnetic layer, the average size of Terfenol-D granules in the ferromagnetic layer, and temperature have been established.The conditions for obtaining the maximum CME effect are analyzed.In [44], the authors proposed a new bilayer-laminated ME composite consisting of magnetostrictive Ni and Terfenol-D plates and a piezoelectric PZT plate.The bias magnetic field and the electric field frequency dependences of the CME coefficient were investigated.It was shown that the Ni/Terfenol-D/Ni/PZT exhibits a large CME coefficient of 6.2 × 10 −7 s/m at the electric field frequency of 42 kHz under a low bias magnetic field of 230 Oe.Along with the listed works, which present a measurement technique and provide estimates of the CME effect, there are known, as noted earlier, works that describe the first practical applications of the CME effect.Chu et al. [16] studied the potential of an electrically driven bulk magnetic field sensor based on the CME effect.It was experimentally found that a limit of detection of 115 pT for a magnetic field of 10 Hz and 300 pT for a magnetic field of 1 Hz was achieved by exciting the ME laminate by 1 V without any bias field.In this case, the power consumption for the ME laminate is only 0.56 mW, which is much lower compared to tens of milliwatts (10-100 mW) for optically pumped or fluxgate sensors and also shows advantages over conventional ME magnetic field sensors based on the DME effect with a current pump.In [17], the authors presented the equivalent magnetic noise of a multi-push-pull magneto (Elasto) electric sensor, implemented as a magnetometer with phase modulation circuitry.Based on the constitutive equations and the Nyquist theorem, they have determined the expected sensitivities and the equivalent magnetic noise for the phase modulation and provided a comparison with experimental results.For measurement applications, the authors [18] investigated the CME effect on a thin-film silicon cantilever of a size of 25 mm × 2.45 mm × 0.35 mm using a piezoelectric AlN film and a magnetostrictive FeCoSiB film, each 2 µm thick.The measured sensitivity at the resonant frequency of 515 kHz increases nearly linearly with a carrier signal amplitude, reaching 64 kV/T at 160 mV.The DC detection limit is 210 pT/Hz1/2, and also about 70 pT/Hz1/2 at 10 Hz without magnetic field bias.The obtained high characteristics are of great interest for biomagnetic applications.Recently, research into ME transmitting antennas based on the CME effect in the low-frequency range has been very active [19][20][21] since their creation will significantly reduce the size and ensure effective underwater and underground communications.Yao et al. [19] proposed a bulk acoustic wave-mediated multiferroic antenna structure.Its potential for efficient radiation of electromagnetic waves is evaluated by analytically deriving the lower bound of its Q factor.A one-dimensional multiscale finite-difference time-domain technique was developed to predict the bilateral, dynamic coupling between acoustic waves and electromagnetic waves.Dong et al. [20] presented a novel, very low-frequency communication system using one pair of ME antennas.With 80 V driving voltage, the power consumption of the ME transmitter has been measured as 400 mW at the maximum communication distance of 120 m.In [21], Y. Wang et al. investigated a ME antenna based on mechanical resonance made by MEMS technology.The ME resonant disk structure with a diameter of 70 µm and a thickness of 1 µm with SiO 2 /Cr/Au/AlN/Cr/Au/FeGaB-stacked layers was prepared on a 300 µm silicon wafer and showed a giant ME coefficient of 2.928 kV/cm/Oe in resonance at 224.1 kHz.The article by [22] is devoted to assessing the CME effect within the framework of a two-dimensional multifield communication model for an acoustically controlled antenna.To compare the results obtained in our paper, it is necessary to clarify the details of the calculation in [22].In addition to the use of the CME effect in sensors and antennas, this effect, as noted, is beginning to be widely explored for applications in spintronics [23][24][25][26][27][28] and ME memory devices [29][30][31].Zavaliche et al. [23] used piezoelectric force microscopy and magnetic force microscopy to locally image the coupled piezoelectric-magnetic switching in epitaxial ferroelectric BiFeO 3 -ferrimagnetic CoFe 2 O 4 columnar nanostructures.They presented a perpendicular CME susceptibility of 1.0 × 10 −2 G cm/V induced by an electric field.In [24], the authors showed that the magnetization of the NiFe/CoFe film on PZT substrate changes by around 6% at the low voltage of 60 V, applied to PZT substrate and is up to two times larger at 800 V in an external magnetic field of 50 Oe.The use of voltage also allows the reversible adjustment of the magnetization orientation in ferromagnetic layers.Chopdekar et al. [25] investigated a (011)-oriented ferroelectric PMN-PT substrate with an epitaxial ferromagnetic (La,Sr)MnO 3 film and showed that the electric field pulses of 6 kV/cm induce large, reversible, and bistable remanent strains.The authors demonstrated that an electric field pulse can be used to 'set' and 'reset' the magnetic anisotropy orientation and resistive state in the film for the electric-field manipulation of nanostructures at room temperature.Spaldin and Ramesh in [26] discussed that the change in magnetic properties by an electric field in ME multiferroic materials has driven significant research activity, with the goal of realizing their transformative technological potential.Ghidini et al. [27,28] investigated the Ni films coupled via strain to ferroelectric substrates of BTO and PMN-PT.The authors studied the shear strain-induced CME effects in Ni thin films on PMN-PT and the reversible switching of magnetization in Ni films on barium titanate.Eerenstein et al. [29] showed electrically induced giant, sharp, and persistent magnetic changes (up to 2.3 × 10 −7 sm −1 ) in ferromagnetic 40 nm LSMO films on 0.5 mm ferroelectric BaTiO 3 substrate.Observed X-ray diffraction confirms strain coupling via ferroelastic BaTiO 3 domains.Hu and Nan [30] presented a phenomenological model for the electric-field-induced reorientation of the magnetic easy axis in ferromagnetic Fe, CFO, Ni, Fe 3 O 4 /ferroelectric PZN-PT, BTO, and PZT-layered heterostructures.The authors showed that as the applied electric field increases, the easy axes tend to switch the directions in CFO and Ni films and highlight that this effect has great potential in prototype electrically controlled magnetic recording devices.Wang et al., in their review [31], examined the theoretical and experimental possibilities of an electric field for controlling magnetism through strain coupling in CoFeB/PMN-PT and discussed the problems and prospects of using this method to create ME memory.The authors showed that investigating the effect of low-frequency CME is useful for future research.
The structure of the article is as follows.The purpose of the article is indicated in the introduction.Sections 2 and 3 are devoted to the consideration of the longitudinal and bending modes in symmetric and asymmetric ME structures.Sections 5 and 6 describe the longitudinal-shear and torsional modes-in such structures.Note that in Sections 4 and 8, as a special case, the quasi-static regime is also considered.A discussion of the calculation method specifics used is included in Section 9. Section 10 presents the conclusions of this paper.

EMR Longitudinal Mode
We accept the alternating magnetic field intensity in the ME composite's magnetostrictive phase to be zero according to [32]: The full thickness of the ME composite is where p t and m t are the piezoelectric and magnetostrictive phase thicknesses.The magnetostrictive and piezoelectric volume fractions are Figure 1 shows the ME composite's structure for calculating the longitudinal mode of the CME effect: Figure 1 shows the ME composite's structure for calculating the longitudinal mode of the CME effect: The strain tensor longitudinal component for magnetostrictive and piezoelectric phases is where x U is the longitudinal displacement.
The material equation of the piezoelectric phase is where 3 E is the alternating electric field intensity applied to the piezoelectric phase of the ME composite, 11 p s is the piezoelectric phase's elastic compliance coefficient, and The material equations of the magnetostrictive phase are ( ) where m Y is the Young's modulus of the magnetostrictive phase, 11 q is the piezomag- netic modulus, and μ is the magnetostrictive phase's permeability.Equations ( 5) and ( 6) take the form for 1 0 h = : Substitute Equation (7) into Equation ( 8): The longitudinal component of the ME composite stress tensor is: The strain tensor longitudinal component for magnetostrictive and piezoelectric phases is where U x is the longitudinal displacement.
The material equation of the piezoelectric phase is where E 3 is the alternating electric field intensity applied to the piezoelectric phase of the ME composite, p s 11 is the piezoelectric phase's elastic compliance coefficient, and d 31 is piezoelectric coefficient.The material equations of the magnetostrictive phase are where m Y is the Young's modulus of the magnetostrictive phase, q 11 is the piezomagnetic modulus, and µ is the magnetostrictive phase's permeability.Equations ( 5) and ( 6) take the form for h 1 = 0: Substitute Equation (7) into Equation (8): The longitudinal component of the ME composite stress tensor is: where Sensors 2024, 24, 151 6 of 32 The motion equation for deformations of ME composite is [47] Substitute Equation (10) into Equation (12): where ω is the angular frequency of the ME composite's deformations, and the composite effective density is where p ρ and m ρ are the piezoelectric and magnetostrictive phase densities.
The solution to the motion equation for deformations of ME composite is where the wave number is Then, the longitudinal component of the strain tensor has the form: The equilibrium conditions for a free ME composite are We obtained a linear inhomogeneous equation system for two unknowns, A and B, by substituting Equations ( 17) and (10) into Equation (18), and found two unknowns A and B from this equation system and substituted them into Equation (17).Then, the obtained expression for the longitudinal component of the strain tensor is substituted into Equation ( 9), and we find the alternating magnetic induction occurring in the magnetostrictive phase during the CME effect: where η = kL 2 .The alternating magnetic induction in the ME composite averaged over the magnetostrictive phase volume during the CME effect is as follows: The CME coefficient is The fundamental resonant frequency for this case is Sensors 2024, 24, 151 7 of 32

Rigid End of the Plate
Boundary conditions for the case when the left end of the ME composite is rigidly fixed, and the right one is free: Substitute Equations ( 10) and ( 15) into Equation ( 23): Then, we obtain Completely analogous to the case with a free ME composite, we find The alternating magnetic induction in the ME composite averaged over the magnetostrictive phase volume during the CME effect is The CME coefficient is The fundamental resonant frequency for this case is ( ) The fundamental resonant frequency for this case is (29)

Free Plate
Figure 2 shows the asymmetric ME composite's structure for calculating the longitudinal and bending mode of the CME effect: The CME coefficient in the resonant regime of the longitudinal mode can be found in Equation ( 21) for an asymmetric free ME composite, and only in Equations ( 1) and ( 2) is it necessary to remove the number 2 before m t .Also, it is necessary to remove the number 2 when integrating the thickness of the magnetostrictive phase in Equation ( 20).The CME coefficient in the resonant regime of the longitudinal mode can be found in Equation ( 21) for an asymmetric free ME composite, and only in Equations ( 1) and ( 2) is it necessary to remove the number 2 before m t.Also, it is necessary to remove the number 2 when integrating the thickness of the magnetostrictive phase in Equation (20).

Rigid End of the Plate
The CME coefficient in the resonant regime of the longitudinal mode can be found in Equation (28) for an asymmetric ME composite with a rigidly fixed left-end and the right one is free, and only in Equations ( 1) and ( 2) is it necessary to remove the number 2 before m t.Also, it is necessary to remove the number 2 when integrating the thickness of the magnetostrictive phase in Equation (27). Figure 3 below shows the dependence of the CME coefficient on the alternating electric field frequency applied to the ME composite's piezoelectric phase for the longitudinal mode of the ME effect.The following parameters were used in the calculations: the magnetostrictive phase is Metglas with a thickness of m t = 29 µm; piezoelectric is PZT (LN cut y + 128 • ) with a thickness of p t = 0.5 mm; the resonance quality factor Q = 130; the length of the ME composite is 10 mm.Piezoelectric PZT is ceramic that is polarized along the thickness direction (the Z-axis is index 3).The Metglas width is equal to the piezoelectric width, and the ME composite width is much greater than its thickness; the ME composite's width is less than its length.To take into account losses in the calculation, it is assumed ω = 2π 1 + 1 Q i f .The resonant frequency calculated by Equation ( 22) for the longitudinal mode of the ME effect for the symmetric free ME composite Metglas/PZT/Metglas is 165.8 kHz; for the symmetric free ME composite, Metglas/LN cut y + 128 • /Metglas is 267 kHz.The resonant frequency calculated by Equation ( 29) for the longitudinal mode of the ME effect for symmetric ME composite with rigidly fixed left-end Metglas/PZT/Metglas is 82.9 kHz, and for symmetric ME composite with rigidly fixed left-end Metglas/LN cut y + 128 • /Metglas, it is 133.5 kHz.The resonant frequency, calculated by Equation (22) for the longitudinal mode of the ME effect for the asymmetric free ME composite Metglas/PZT is 164.6 kHz, and for the asymmetric free ME composite Metglas/LN cut y + 128 • , it is 272.6 kHz.The resonant frequency calculated by Equation ( 29) for the longitudinal mode of the ME effect for asymmetric ME composite with rigidly fixed left-end Metglas/PZT is 82.3 kHz, and for asymmetric ME composite with rigidly fixed left-end Metglas/LN cut y + 128 • , it is 136.3 kHz.
The CME coefficient value is two times less for an asymmetric ME composite for the longitudinal mode of the CME effect in comparison with a symmetric ME composite because an asymmetric ME composite has the volume of the magnetostrictive phase and is two times less than for a symmetric one.The resonant frequency of the CME effect for the longitudinal mode is greater for the composite with LN cut y + 128 • as a piezoelectric than for a composite with PZT as a piezoelectric due to the lower effective density and greater effective stiffness of the composite with LN cut y + 128 • .The CME coefficient value is greater for a composite with PZT as a piezoelectric than for a composite with LN cut y + 128 • , which is due to a much larger piezoelectric coefficient PZT compared to LN cut y + 128 • and a greater elastic compliance coefficient for the PZT than for LN cut y + 128 • .The CME effect resonant frequency is two times less for the ME composite with a fixed left end than for the free ME composite, while the CME coefficient resonant value almost does not change.
resonant frequency, calculated by Equation ( 22) for the longitudinal mode of the ME effect for the asymmetric free ME composite Metglas/PZT is 164.6 kHz, and for the asymmetric free ME composite Metglas/LN cut y + 128°, it is 272.6 kHz.The resonant frequency calculated by Equation ( 29) for the longitudinal mode of the ME effect for asymmetric ME composite with rigidly fixed left-end Metglas/PZT is 82.3 kHz, and for asymmetric ME composite with rigidly fixed left-end Metglas/LN cut y + 128°, it is 136.3 kHz.The CME coefficient value is two times less for an asymmetric ME composite for the longitudinal mode of the CME effect in comparison with a symmetric ME composite because an asymmetric ME composite has the volume of the magnetostrictive phase and is

EMR Bending Mode
The calculation is performed for an asymmetric ME composite (see Figure 2) for the bending mode in the EMR region.

Free Plate
The full thickness of the ME composite is The longitudinal component of the strain tensor for magnetostrictive and piezoelectric phases is [47] where w is the transverse displacement.
The longitudinal component of the magnetostrictive and piezoelectric phases stress tensor, electric field intensity's third component applied to the piezoelectric phase of the ME composite, and the alternating magnetic induction in the magnetostrictive phase for where [47] ε T 33 is the permittivity tensor's component of the piezoelectric phase, and p D 3 is the piezoelectric phase's electrical displacement.
The alternating magnetic induction in the magnetostrictive phase is determined by Equation (9).
The ME composite's bending moment is where b is the composite width, z 0 is the position of the boundary between the piezoelectric and magnetostrictive phases relative to the ME composite's neutral line.
Let us find the electric voltage on the piezoelectric: Find the electrical displacement in piezoelectric p D 3 from Equation (37) and substitute into Equation ( 34) where z 0 is determined from the minimum condition c 11 : The shear force is [47]: The motion equation for bending vibrations of ME composite is [47] Substitute Equation ( 41) for shear force into Equation (42): The dependence of the displacement on time is harmonic w ∼ e iωτ , therefore: where the composite effective density is determined by Equation ( 14).
The general solution to the motion equation is Then, we obtain The equilibrium conditions for a free ME composite are We obtain a linear inhomogeneous equation system for four unknowns C 1 , C 2 , C 3 , C 4 by substituting Equations ( 38) and (41) into Equation (47).We find four unknowns C 1 , C 2 , C 3 , C 4 from this equation system.
Then, substitute the obtained constants C 1 , C 2 , C 3 , C 4 into Equation ( 46) and find the alternating magnetic induction in the ME composite averaged over the magnetostrictive phase volume during the CME effect: where Then, the CME coefficient is calculated by the equation The fundamental resonant frequency for this case is

Rigid End of the Plate
The boundary conditions for the case when the left end of the ME composite is rigidly fixed and the right one is free: The general solution to the motion equation corresponds to Equation (45).
After substituting Equations ( 45), ( 38) and ( 41) into Equation (52), we obtain a linear inhomogeneous system for determining unknown constants C 1 , C 2 , C 3 , C 4 .Having found them from this system, we substitute them into Equation (46) and then find the alternating magnetic induction averaged over the volume of the ME composite during the CME effect: Then, the CME coefficient is calculated by the equation The fundamental resonant frequency for this case is Figure 4 below shows the dependence of the CME coefficient on the alternating electric field frequency applied to the ME composite's piezoelectric phase for the bending mode of the ME effect.The ME composite's material parameters are the same as for calculating the dependence of the CME coefficient on the alternating electric field frequency for the longitudinal mode.The CME effect resonant frequency for the bending mode is greater for a composite with the LN section y + 128° as a piezoelectric phase than for a composite with PZT due to the lower effective density and greater effective stiffness of the composite with LN cut y + 128°.The value of the CME coefficient is greater for a composite with PZT as a piezoelectric phase than for a composite with LN cut y + 128°, which is due to a much higher value of the piezoelectric coefficient PZT compared to LN cut y + 128°-a greater elastic compliance coefficient for the PZT than for LN cut y + 128°.The CME effect resonant frequency is three times less for the case when the composite's left end is rigidly fixed than for the free ME composite, while the CME coefficient resonant value increases by more than three times for the bending mode of vibrations of the ME composite.The CME effect resonant frequency is much less for the bending mode than for the longitudinal mode for the free ME composite and the composite with the rigidly fixed left end, while the CME coefficient resonant values for the bending mode of vibrations are much less than for the longitudinal mode.

Quasi-Static Regime
When the bias field 0 H is oriented along the ME composite length, the CME coef- ficient in the quasi-static regime for a symmetric ME composite is determined by the contribution of only the longitudinal mode; for an asymmetric ME composite, the CME coefficient in the quasi-static regime is determined by the contribution of both the longitudinal and bending modes.

Symmetric Structure
Assuming in Equation ( 21) the frequency f is equal to zero, we obtain:  The CME effect resonant frequency for the bending mode is greater for a composite with the LN section y + 128 • as a piezoelectric phase than for a composite with PZT due to the lower effective density and greater effective stiffness of the composite with LN cut y + 128 • .The value of the CME coefficient is greater for a composite with PZT as a piezoelectric phase than for a composite with LN cut y + 128 • , which is due to a much higher value of the piezoelectric coefficient PZT compared to LN cut y + 128 • -a greater elastic compliance coefficient for the PZT than for LN cut y + 128 • .The CME effect resonant frequency is three times less for the case when the composite's left end is rigidly fixed than for the free ME composite, while the CME coefficient resonant value increases by more than three times for the bending mode of vibrations of the ME composite.The CME effect resonant frequency is much less for the bending mode than for the longitudinal mode for the free ME composite and the composite with the rigidly fixed left end, while the CME coefficient resonant values for the bending mode of vibrations are much less than for the longitudinal mode.

Quasi-Static Regime
When the bias field H 0 is oriented along the ME composite length, the CME coefficient in the quasi-static regime for a symmetric ME composite is determined by the contribution of only the longitudinal mode; for an asymmetric ME composite, the CME coefficient in the quasi-static regime is determined by the contribution of both the longitudinal and bending modes.

Symmetric Structure
Assuming in Equation ( 21) the frequency f is equal to zero, we obtain:

Asymmetric Structure
Consider the motion equation for deformations: Equation ( 57) has the form of the quasi-static regime: This means that T 1 should not depend on x.This means that S 1 should also not depend on x.Since both the longitudinal and bending modes are excited in the asymmetric ME composite for the quasi-static regime, then The magnetostrictive phase's material equations are Equations ( 60) and (61) have the form for h 1 = 0: The longitudinal component of the piezoelectric phase stress tensor and electric field intensity third component, applied to the piezoelectric phase of the ME composite, are The electrical voltage on the piezoelectric, taking into account that D 3 does not depend on z, is as follows: Let us express the electric displacement's third component in the piezoelectric phase from Equation (66): Substitute Equation (59) into Equation (62): The first condition for the static equilibrium of an ME composite is that the total longitudinal force is equal to zero: The second condition for the static equilibrium of the ME composite is the zero total moment, which is given by: We use the well-known condition to determine z 0 , which is that the longitudinal force does not depend on the coefficient B associated with bending vibrations, and the bending moment does not depend on the coefficient A associated with longitudinal vibrations.This condition is Then, we obtain Then, we find the alternating magnetic induction averaged over the volume of the ME composite m B 1 during the CME effect, substituting in Equation (63) the obtained constants A and B; then, the CME coefficient is as follows: Figure 5 below shows the dependence of the CME coefficient on the piezoelectric volume fraction for the quasi-static regime of the CME effect.The material parameters of the ME composites are the same as for calculating the dependence of the CME coefficient on the alternating electric field frequency for the longitudinal mode of the CME effect.
The optimal piezoelectric volume fraction of a symmetric ME composite is 0.55, at which the CME coefficient is maximal for a composite with PZT and 0.45 for a composite with LN cut y + 128 • .The optimal piezoelectric volume fraction is 0.8 for an asymmetric ME composite with PZT, and the optimal piezoelectric volume fraction is 0.75 for a composite with LN cut y + 128 • .The CME coefficient value is greater for a composite with PZT as a piezoelectric phase than for a composite with LN cut y + 128 • , which is due to a much larger piezoelectric coefficient PZT compared to LN cut y + 128 • .The optimal piezoelectric volume fraction of a symmetric ME composite is 0.55, at which the CME coefficient is maximal for a composite with PZT and 0.45 for a composite with LN cut y + 128°.The optimal piezoelectric volume fraction is 0.8 for an asymmetric ME composite with PZT, and the optimal piezoelectric volume fraction is 0.75 for a composite with LN cut y + 128°.The CME coefficient value is greater for a composite with PZT as a piezoelectric phase than for a composite with LN cut y + 128°, which is due to a much larger piezoelectric coefficient PZT compared to LN cut y + 128°.

EMR Longitudinal Shear Mode
The equations for calculating the CME coefficient during the longitudinal shear mode are the same as for the longitudinal mode, but the longitudinal mechanical stress tensors are replaced by с is on 66 с .
Figure 6 shows the symmetric ME composite's structure for calculating the longitudinal shear mode of the CME effect.

EMR Longitudinal Shear Mode
The equations for calculating the CME coefficient during the longitudinal shear mode are the same as for the longitudinal mode, but the longitudinal mechanical stress tensors are replaced by m T 1 , p T 1 on the shear mechanical stress tensors m T 6 , p T 6 , the piezoelectric and piezomagnetic modulus d 31 , q 11 are on d 36 , q 16 , the piezoelectric phase compliance coefficient p s 11 is on p s 66 , the magnetostrictive phases Young's modulus m Y is on shift modulus m G, and the composite's effective stiffness coefficient c 11 is on c 66 .
Figure 6 shows the symmetric ME composite's structure for calculating the longitudinal shear mode of the CME effect.

Symmetric Structure
The CME coefficient is calculated by the following equation:

Symmetric Structure
The CME coefficient is calculated by the following equation: where The fundamental resonant frequency for this case is [47]

Asymmetric Structure
The CME coefficient can be found in Equation ( 75) for an asymmetric ME composite in the resonant mode of the longitudinal-shear ME mode for a free composite, and it is necessary to remove the number 2 before m t in Equations ( 1), ( 2) and (75).

EMR Torsional Mode for ME Composite of Metglas/GaAs
The calculation is performed for an asymmetric ME composite (see Figure 7) for the torsional mode in the EMR region.The full thickness of the ME composite is determined by Equation (30).The ME composite's strain tensor shear components are as follows: The X-axis is drawn along the composite length in the corresponding symmetry plane of the composite, along the composite beam rotation axis during torsional vibrations in the direction of the composite length.The Y-axis will be directed along the composite width.The bias field H 0 is directed along the Y-axis, and the alternating electric field p E 3 is directed along the Z-axis.
The full thickness of the ME composite is determined by Equation (30).
The ME composite's strain tensor shear components are as follows: where θ is the twist angle.
The material equation for the piezoelectric phase is where p G is the piezoelectric phase's shift modules.We find the piezoelectric phase's stress tensor tangent components from Equation (82): The material equations of the magnetostrictive phase for h 1 = 0 as follows: We find the magnetostrictive phase's stress tensor tangent components from Equation (85): The ME composite's magnetic induction for h 1 = 0 is: The electric displacement in piezoelectric is: where ε is the piezoelectric phase's permittivity.We express p E 3 from Equation (88): where and substitute it into Equation (84).As a result, we obtain The torque is as follows: where where z 0 is the position of the boundary between the piezoelectric and magnetostrictive phases relative to the ME composite's rotation axis.
Let us find the electric voltage on the piezoelectric: We express the electric displacement in a piezoelectric from Equation (95) and substitute it into Equation (93): z 0 is determined from the condition of the minimum effective shear modulus of the ME composite G : The torsional vibrations are as follows: where the ME composite's moment of inertia per unit width is Sensors 2024, 24, 151 20 of 32 where the piezoelectric and magnetostrictive phases polar moments are Substitute Equation (97) into Equation (100): The dependence of the twist angle on time is harmonic θ ∼ e iωτ , therefore: where the wave number is The general solution to Equation ( 104) is as follows: The equilibrium conditions for a free ME composite are We obtain a linear inhomogeneous equation system for two unknowns A and B by substituting Equation (97) into Equation (107).Then, we find two unknowns A and B, from this equation system and substitute them into Equation (87).After that, we find the alternating magnetic induction occurring in the magnetostrictive phase m B 1 during the CME effect.And we will find the CME coefficient: where The fundamental resonant frequency for this case is Figure 8 below shows the dependence of the CME coefficient on the alternating electric field frequency applied to the ME composite's piezoelectric phase for the longitudinal-shear and torsional modes.The magnetostrictive phase is Metglas, with a thickness of m t = 29 µm.The piezoe- lectric is GaAs with a thickness of p t = 0.2 mm.The value of the quality factor Q = 300.
The ME composite length was 23 mm, and the width was 0.3 mm.The magnetostrictive phase is Metglas, with a thickness of m t= 29 µm.The piezoelectric is GaAs with a thickness of p t= 0.2 mm.The value of the quality factor Q= 300.The ME composite length was 23 mm, and the width was 0.3 mm.
The resonant frequency calculated by Equation (80) for the longitudinal shear mode of the CME effect for the symmetric free ME composite Metglas/GaAs/Metglas is 67.1 kHz, and for the asymmetric free ME composite Metglas/GaAs, it is 69.5 kHz.The resonant frequency, calculated by Equation (110) for the torsional mode of the CME effect for the asymmetric free ME composite Metglas/GaAs, is 67.8 kHz.
The resonant frequency of the CME effect is slightly less for the longitudinal shear mode for an asymmetric ME composite with GaAs as a piezoelectric phase than for the torsional mode, while the resonant value of the CME coefficient is more than ten times greater for the longitudinal shear mode than for the torsional mode.The value of the CME coefficient is almost two times less for an asymmetric ME composite for the longitudinal shear mode in comparison with a symmetric ME composite because the volume of the magnetostrictive phase is two times less for an asymmetric ME composite than for a symmetric one.It was determined in previous studies that in an asymmetric ME composite, with the orientation of the bias field shown in Figure 6, both longitudinal shear and torsional modes exist simultaneously during the CME effect.Moreover, the resonant frequency of the longitudinal shear mode and the torsional mode differ slightly, and the CME coefficient's value for the longitudinal shear mode is an order of magnitude greater than for the torsional mode.Because of this, the contribution of the torsional mode to the full CME coefficient is insignificant against the background of the longitudinal shear mode, and it is difficult to determine experimentally.One of the possible ways out of this situation may be the use of a bimorph LN Zy + 45 • as a piezoelectric phase in the ME composite.

EMR Torsional Mode in a ME Composite Based on Bimorph LN
The X-axis is drawn along the composite length in the corresponding symmetry plane of the composite, along the composite beam rotation axis during torsional vibrations in the direction of the composite length.The Y-axis will be directed along the composite width.The bias field H 0 is directed along the Y-axis, and the alternating electric field p E 3 is directed along the Z-axis.
Figure 9 shows the asymmetric ME composite's structure with bimorph LN Zy + 45 • as a piezoelectric phase for calculating the torsional mode of the CME effect: mode for an asymmetric ME composite with GaAs as a piezoelectric phase than for the torsional mode, while the resonant value of the CME coefficient is more than ten times greater for the longitudinal shear mode than for the torsional mode.The value of the CME coefficient is almost two times less for an asymmetric ME composite for the longitudinal shear mode in comparison with a symmetric ME composite because the volume of the magnetostrictive phase is two times less for an asymmetric ME composite than for a symmetric one.
It was determined in previous studies that in an asymmetric ME composite, with the orientation of the bias field shown in Figure 6, both longitudinal shear and torsional modes exist simultaneously during the CME effect.Moreover, the resonant frequency of the longitudinal shear mode and the torsional mode differ slightly, and the CME coefficient's value for the longitudinal shear mode is an order of magnitude greater than for the torsional mode.Because of this, the contribution of the torsional mode to the full CME coefficient is insignificant against the background of the longitudinal shear mode, and it is difficult to determine experimentally.One of the possible ways out of this situation may be the use of a bimorph LN Zy + 45° as a piezoelectric phase in the ME composite.

EMR Torsional Mode in a ME Composite Based on Bimorph LN
The X-axis is drawn along the composite length in the corresponding symmetry plane of the composite, along the composite beam rotation axis during torsional vibrations in the direction of the composite length.The Y-axis will be directed along the composite width.The bias field 0 H is directed along the Y-axis, and the alternating electric field 3 p E is directed along the Z-axis.
Figure 9 shows the asymmetric ME composite's structure with bimorph LN Zy + 45° as a piezoelectric phase for calculating the torsional mode of the CME effect:  The ME composite's strain tensor shear components are determined by Equation (81).
The material equations of the magnetostrictive phase for h 1 = 0 correspond to Equations ( 86) and (87). where The position of the interface between the piezoelectric and magnetostrictive phases relative to the composite beam rotation axis z 0 is determined from the condition of the minimum polar moments of the shear stiffness coefficient Q: The torsional vibrations are determined by Equation (100), the ME composite's moment of inertia per unit width is determined by Equation ( 101), where the piezoelectric and magnetostrictive phases polar moments are determined by Equation (102).
Substitute Equation (120) into Equation (100): The dependence of the twist angle on time is harmonic θ ∼ e iωτ , therefore: where the wave number is The general solution of Equation ( 125) is determined by Equation ( 106).The equilibrium conditions for a free ME composite correspond to Equation (107).
We obtain a linear inhomogeneous equation system for two unknowns A and B by substituting Equation (120) into equilibrium conditions for a free ME composite (107).Solving this system, we find A and B.Then, we substitute the found constants into Equation (87) and find the alternating magnetic induction occurring in the magnetostrictive phase m B 1 during the CME effect.After that, we will find the CME coefficient: where Sensors The fundamental resonant frequency for this case is   The value of the CME coefficient is much greater for the torsional mode for an asymmetric composite with bimorph LN Zy + 45 • as a piezoelectric phase than for a composite with GaAs.
The longitudinal shear mode of the CME effect in an asymmetric composite with a bimorph LN will not be excited.If some type of CME effect is observed in the experiment, then it can only occur from the torsional mode.Thus, the longitudinal shear mode interfering with the observation of the torsional mode will be eliminated.

ME Composite Based on Bimorph LN
Assuming in Equation (127) the frequency f is equal to zero, we obtain: Figure 11 below shows the dependence of the CME coefficient on the piezoelectric volume fraction for the quasi-static regime of the CME effect for longitudinal-shear and torsional modes.The material parameters of the ME composites are the same as for calculating the dependence of the CME coefficient on the alternating electric field frequency in the EMR region.
The optimal piezoelectric volume fraction of a symmetric ME composite is 0.54, at which the CME coefficient is maximal for a composite with GaAs as piezoelectric and 0.45 for an asymmetric ME composite with GaAs.The optimal piezoelectric volume fraction is 0.68 for an asymmetric ME composite with bimorph LN Zy + 45 • .The maximum value of the CME coefficient is approximately five times greater for an asymmetric ME composite with a bimorph LN Zy + 45 • than for an asymmetric ME composite with GaAs as a piezoelectric in the quasi-static regime.

ME Composite Based on Bimorph LN
Assuming in Equation (127) the frequency f is equal to zero, we obtain: Figure 11 below shows the dependence of the CME coefficient on the piezoelectric volume fraction for the quasi-static regime of the CME effect for longitudinal-shear and torsional modes.The material parameters of the ME composites are the same as for calculating the dependence of the CME coefficient on the alternating electric field frequency in the EMR region.

Discussion
We would first like to discuss the applicability limits of our theoretical calculation of the CME coefficient in this section.One of our main assumptions is that we consider the mechanical coupling between the magnetostrictive and piezoelectric phases to be ideal.Therefore, we consider the corresponding strain tensor components of the magnetostrictive and piezoelectric phases to be equal.In practice, the mechanical coupling between the magnetostrictive and piezoelectric phases may be quite far from ideal.This is due to the methods of connecting the phases: gluing or some other methods.At the same time, the mechanical deformation effect is transferred from the piezoelectric phase to the magne-tostrictive one and is not transferred as well as it is described by theory during the CME effect in the experiment.Therefore, the experimental CME coefficient is always less than the theoretical one.Another assumption we use in modeling is that we consider the alternating magnetic field inside the magnetostrictive phase to be zero.This assumption justifies itself well when the ME composite length is greater than its width and the width is significantly greater than its thickness.In fact, this is not exactly true.A more accurate picture could be given by finite element modeling in such a software package as Comsol Multiphysics 6.2 for the alternating magnetic field distribution inside and around the magnetostrictive phase.Another assumption close in meaning is that we consider the electric induction vector's component directed along the ME composite's thickness independent of the coordinate along the thickness when considering bending and torsional modes.Again, this assumption works well with a long, narrow, very thin ME composite.However, a more accurate calculation could be made by the method of finite element modeling in Comsol Multiphysics.It is also necessary to discuss the complexity of observing the CME effect's torsional mode against the background of the CME effect's longitudinal shear mode.These modes exist simultaneously in an asymmetric ME composite when the CME effect is excited by applying an alternating electric voltage to the piezoelectric phase.Calculations show that the CME effect's main resonant frequency for the longitudinal shear mode for an asymmetric ME composite with GaAs as a piezoelectric phase is slightly less than for the torsional mode, while the resonant value of the CME coefficient is more than ten times greater for the longitudinal shear mode than for the torsional mode.Therefore, it is difficult to observe such a small CME effect from the torsional mode against the background of a relatively large CME effect from the longitudinal shear mode in a small general frequency range.One of the possible ways to solve this problem is to use bimorphic LN Zy + 45 • with the same LN layer thickness with oppositely directed polarization along the ME composite's thickness instead of GaAs.The theory shows that in the ME composite based on the bimorphic LN Zy + 45 • , the longitudinal shear mode is not excited.Therefore, the observation of the CME effect's torsional mode will be greatly simplified.Also, we would like to explain why the article does not compare our theoretical calculations with the experimental data of other authors.This is because we consider relatively simple ME structures with the same length and width as the magnetostrictive and piezoelectric phases in modeling.Other authors use more complex ME structures in their experiments, in which the lengths or widths of the magnetostrictive and piezoelectric phases do not coincide [34,35,40,[48][49][50], and other additional layers are present in the design; for example, a permanent magnet [51].Such complex constructions of ME composites require separate, rather laborious calculations for comparison with our theory when observing the CME effect.

Conclusions
The general theory of the CME effect in composites in the low-frequency range, including the EMR region, is presented.The main EMR modes are considered in more detail, such as longitudinal, bending, longitudinal-shear, and torsional.To demonstrate the theory, expressions for the CME coefficients were obtained for symmetric and asymmetric GaAs/Metglas and LN/Metglas layered structures.To complete the analysis, a brief overview of the main works on the CME effect in the low-frequency range, including in the field of EMR, is given.
This article notes a number of limitations of the presented theory of the CME effect that must be taken into account.This is, first of all, the imperfection of the mechanical connection between the piezoelectric and magnetostrictive phases associated with the commonly used adhesive technology.The following restrictions are due to the presence of an alternating magnetic field inside the magnetostrictive phase and the constant direction of the electrical induction vector along the thickness of the composite in the case of bending and torsional modes.
Recommendations are given for the experimental study of the CME effect in the torsional mode in a bimorph structure based on Metglas/LN, which will eliminate the longitudinal shear mode and study only the spectrum of the torsional mode.To reliably calculate the CME effect, taking into account the complexity of the problem, it is proposed to use our analytical estimates presented in this article and compare them with calculations based on the Comsol Multiphysics package and the results of the experiments.

Figure 3 .
Figure 3. Dependence of the CME coefficient on the alternating electric field frequency applied to the ME composite's piezoelectric phase for the longitudinal mode of the ME effect.(А) Free symmetric ME composite: the blue line is composite Metglas/PZT/Metglas, and the red line is composite Metglas/LN cut y + 128°/Metglas.(В) Symmetric ME composite with a rigidly fixed left end: the blue line is composite Metglas/PZT/Metglas, and the red line is composite Metglas/LN cut y + 128°/Metglas.(С) Free asymmetric ME composite: the blue line is composite Metglas/PZT, and the red line is composite Metglas/LN cut y + 128°.(D) Asymmetric ME composite with a rigidly fixed left end: the blue line is composite Metglas/PZT, and the red line is composite Metglas/LN cut y + 128°.

Figure 3 .
Figure 3. Dependence of the CME coefficient on the alternating electric field frequency applied to the ME composite's piezoelectric phase for the longitudinal mode of the ME effect.(A) Free symmetric ME composite: the blue line is composite Metglas/PZT/Metglas, and the red line is composite Metglas/LN cut y + 128 • /Metglas.(B) Symmetric ME composite with a rigidly fixed left end: the blue line is composite Metglas/PZT/Metglas, and the red line is composite Metglas/LN cut y + 128 • /Metglas.(C) Free asymmetric ME composite: the blue line is composite Metglas/PZT, and the red line is composite Metglas/LN cut y + 128 • .(D) Asymmetric ME composite with a rigidly fixed left end: the blue line is composite Metglas/PZT, and the red line is composite Metglas/LN cut y + 128 • .

Figure 4 .
Figure 4. Dependence of the CME coefficient on the alternating electric field frequency applied to the ME composite's piezoelectric phase for the bending mode of the ME effect.(А) Free asymmetric ME composite: the blue line is composite Metglas/PZT, and the red line is composite Metglas/LN cut y + 128°.(В) Asymmetric ME composite with a rigidly fixed left end: the blue line is composite Metglas/PZT, and the red line is composite Metglas/LN cut y + 128°.

Figure 4 .
Figure 4. Dependence of the CME coefficient on the alternating electric field frequency applied to the ME composite's piezoelectric phase for the bending mode of the ME effect.(A) Free asymmetric ME composite: the blue line is composite Metglas/PZT, and the red line is composite Metglas/LN cut y + 128 • .(B) Asymmetric ME composite with a rigidly fixed left end: the blue line is composite Metglas/PZT, and the red line is composite Metglas/LN cut y + 128 • .The resonant frequency calculated by Equation (51) for the bending mode of the CME effect for the asymmetric free ME composite Metglas/PZT is 19.7 kHz, and for the asymmetric free ME composite Metglas/LN cut y + 128 • , it is 33.1 kHz.The resonant frequency calculated by Equation (55) for the asymmetric ME composite with rigidly fixed left-end Metglas/PZT is 3.1 kHz, and for the asymmetric ME composite with rigidly fixed left-end Metglas/LN cut y + 128 • , it is 5.2 kHz.The CME effect resonant frequency for the bending mode is greater for a composite with the LN section y + 128 • as a piezoelectric phase than for a composite with PZT due to the lower effective density and greater effective stiffness of the composite with LN cut y + 128 • .The value of the CME coefficient is greater for a composite with PZT as a piezoelectric phase than for a composite with LN cut y + 128 • , which is due to a much higher value of the piezoelectric coefficient PZT compared to LN cut y + 128 • -a greater elastic compliance coefficient for the PZT than for LN cut y + 128 • .The CME effect resonant frequency is three times less for the case when the composite's left end is rigidly fixed than for the free ME composite, while the CME coefficient resonant value increases by more than three times for the bending mode of vibrations of the ME composite.The CME effect resonant frequency is much less for the bending mode than for the longitudinal mode for the free ME composite and the composite with the rigidly fixed left end, while the CME coefficient resonant values for the bending mode of vibrations are much less than for the longitudinal mode.

Figure 5 .
Figure 5. Dependence of the CME coefficient on the piezoelectric volume fraction in the quasi-static regime.(А) The blue line is the symmetric composite Metglas/PZT/Metglas, and the red line is the symmetric composite Metglas/LN cut y + 128°/Metglas.(В) The blue line is asymmetric composite Metglas/PZT, and the red line is asymmetric composite Metglas/LN cut y + 128°.

Т
T , the piezoelectric and piezomagnetic modulus 31 d , 11q are on 36 d , 16 q , the piezoelectric phase compliance coefficient 11 p s is on 66 p s , the magnetostrictive phases Young's modulus m Y is on shift modulus m G , and the composite's effective stiffness coefficient 11

Figure 5 .
Figure 5. Dependence of the CME coefficient on the piezoelectric volume fraction in the quasi-static regime.(A) The blue line is the symmetric composite Metglas/PZT/Metglas, and the red line is the symmetric composite Metglas/LN cut y + 128 • /Metglas.(B) The blue line is asymmetric composite Metglas/PZT, and the red line is asymmetric composite Metglas/LN cut y + 128 • .

Figure 6 .
Figure 6.Symmetric ME composite scheme for the longitudinal shear mode.

Figure 6 .
Figure 6.Symmetric ME composite scheme for the longitudinal shear mode.

Figure 7 .H
Figure 7. Asymmetric ME composite scheme for the torsional mode.The orientations of the crystallographic axes are [100], [010], and [001], where [100] is the length (X-axis is index 1), [010] is the width (Y-axis is index 2), and [001] is a thickness (Z-axis is index 3), respectively, for GaAs as a ME composite's piezoelectric.The X-axis is drawn along the composite length in the corresponding symmetry plane of the composite, along the composite beam rotation axis during torsional vibrations in the direction of the composite length.The Y-axis will be directed along the composite width.The bias field

Figure 7 .
Figure 7. Asymmetric ME composite scheme for the torsional mode.The orientations of the crystallographic axes are [100], [010], and [001], where [100] is the length (X-axis is index 1), [010] is the width (Y-axis is index 2), and [001] is a thickness (Z-axis is index 3), respectively, for GaAs as a ME composite's piezoelectric.The X-axis is drawn along the composite length in the corresponding symmetry plane of the composite, along the composite beam rotation axis during torsional vibrations in the direction of the composite length.The Y-axis will be directed along the composite width.The bias field H 0 is directed along the Y-axis, and the alternating electric field p E 3 is directed along the Z-axis.The full thickness of the ME composite is determined by Equation(30).The ME composite's strain tensor shear components are as follows:

Figure 8
Figure 8 below shows the dependence of the CME coefficient on the alternating electric field frequency applied to the ME composite's piezoelectric phase for the longitudinalshear and torsional modes.

Figure 8 .
Figure 8. Dependence of the CME coefficient on the alternating electric field frequency applied to the free ME composite's piezoelectric phase for the longitudinal-shear and torsional modes.(А) Longitudinal shear mode, a symmetric free ME composite Metglas/GaAs/Metglas; (В) longitudinalshear mode, an asymmetric free ME composite Metglas/GaAs; (С) torsional mode, a symmetric free ME composite Metglas/GaAs.

Figure 8 .
Figure 8. Dependence of the CME coefficient on the alternating electric field frequency applied to the free ME composite's piezoelectric phase for the longitudinal-shear and torsional modes.(A) Longitudinal shear mode, a symmetric free ME composite Metglas/GaAs/Metglas; (B) longitudinal-shear mode, an asymmetric free ME composite Metglas/GaAs; (C) torsional mode, a symmetric free ME composite Metglas/GaAs.

Figure 9 .
Figure 9. Asymmetric ME composite scheme with bimorph LN Zy + 45° phase.A bimorph LN Zy + 45° consists of two LN Zy + 45° plates having the opposite polarization direction along the thickness.The crystallographic axis, corresponding to the thickness of the bimorph layer for one plate, is [101] (the Z-axis is index 3), and for the other plate, it is [−10−1].Crystallographic axes, directed along the length and width, are [10−1] (the X-axis is index 1) and [010] (the Y-axis is index 2), respectively.The ME composite's strain tensor shear components are determined by Equation (81).

Figure 9 .
Figure 9. Asymmetric ME composite scheme with bimorph LN Zy + 45 • phase.A bimorph LN Zy + 45 • consists of two LN Zy + 45 • plates having the opposite polarization direction along the thickness.The crystallographic axis, corresponding to the thickness of the bimorph layer for one plate, is [101] (the Z-axis is index 3), and for the other plate, it is [−10−1].Crystallographic axes, directed along the length and width, are [10−1] (the X-axis is index 1) and [010] (the Y-axis is index 2), respectively.The ME composite's strain tensor shear components are determined by Equation (81).The material equations of the magnetostrictive phase for h 1 = 0 correspond to Equations (86) and (87).

)Figure 10 Figure 10
Figure 10 below shows the dependence of the CME coefficient on the alternating electric field frequency applied to the ME composite's piezoelectric phase for the longitudinalshear and torsional modes for a composite with a bimorph LN Zy + 45 • .

Figure 10 .
Figure 10.Dependence of the CME coefficient on the alternating electric field frequency applied t the free ME composite's piezoelectric phase for the torsional mode.Asymmetric free ME composit bimorph LN Zy + 45°/Metglas.

Figure 10 .
Figure 10.Dependence of the CME coefficient on the alternating electric field frequency applied to the free ME composite's piezoelectric phase for the torsional mode.Asymmetric free ME composite bimorph LN Zy + 45 • /Metglas.The following material parameters were used in the calculation: The length of the ME composite was L = 23 mm, and the width was b = 0.5 mm.In the calculation, the following material parameters of the initial components were used: For Metglas: m ρ = 7180 kg/m 3 , m G = 3.85 × 10 10 Pa, q 16 = 1.0 × 10 −9 m/A, and m t = 29 µm.For LN Zy + 45 • : p ρ = 4647 kg/m 3 , p c E 55 = 6.75 × 10 10 Pa, p c E 56 = 7.5 × 10 9 Pa, p c E 66 = 6.75 × 10 10 Pa, ε 33 = 36.5,p1 e 35 = − p2 e 35 = 2.5 C/m 2 , p1 e 36 = − p2 e 36 = 2.5 C/m 2 , and p t = 0.4 mm.The value of the quality factor for this ME composite Q = 100.The resonant frequency calculated by Equation(129) for the torsional mode of the CME effect for the asymmetric free ME composite bimorph LN Zy + 45 • /Metglas is 85.2 kHz.The value of the CME coefficient is much greater for the torsional mode for an asymmetric composite with bimorph LN Zy + 45 • as a piezoelectric phase than for a composite with GaAs.The longitudinal shear mode of the CME effect in an asymmetric composite with a bimorph LN will not be excited.If some type of CME effect is observed in the experiment, then it can only occur from the torsional mode.Thus, the longitudinal shear mode interfering with the observation of the torsional mode will be eliminated.
D 3 = e 35 S 5 + e 36 S 6 + ε 33 ε 0 66 are the shear components at a constant electric field intensity of the piezoelectric phase's stiffness tensor, and e 35 , e 36 are the piezoelectric coefficients at a constant electric field intensity.