Modeling and Experimental Study of Oil-Cooled Stacked Giant Magnetostrictive Actuator for Servo Valve

: Giant magnetostrictive materials (GMMs) have broad application prospects in the ﬁeld of servo valves, but the giant magnetostrictive actuator (GMA) has problems such as large loss and severe heat generation, which a ﬀ ect the output e ﬀ ect and accuracy. To solve these problems, this paper designs a stacked giant magnetostrictive actuator (SGMA) and analyzes the magnetic circuit and magnetic ﬁeld distribution of the SGMA. Based on the magnetic ﬁeld analysis and the Jiles–Atherton model, we analyze the SGMA magnetization model, simplify the traditional model, and give a solution for the simpliﬁed model using the Runge–Kutta method. We analyze the eddy current loss of the SGMA, and according to Bessel’s equation and the Kelvin function, we calculate the relationship among eddy current loss, GMM rod radius, and magnetic ﬁeld frequency. By analyzing the inherent hysteresis of GMMs, a hysteresis loss model of the SGMA is established in this paper. We also calculate the coil impedance and obtain the coil loss model. Based on the loss model, the SGMA cooling system is designed. Based on the above analysis, we design a SGMA prototype, set-up the corresponding experimental platform, and conduct the necessary experiments. The experimental results show that the SGMA responds well to di ﬀ erent signals, but as frequency increases, attenuation, deformation, and hysteresis become more pronounced, which veriﬁes the amplitude and phase changes caused by various losses in the theoretical analysis. The experiment also observes the temperature rise of the oil-cooled SGMA at di ﬀ erent frequencies, indicating that the cooling system can e ﬀ ectively control the temperature change of the SGMA, which validates the foregoing analysis.


Introduction
In 1842, British physicist James Joule discovered that the shape and volume of nickel (Ni) would change under the action of an external magnetic field. This is magnetostriction. In 1974, Clark et al. discovered that the ternary rare earth alloy Tb 1-x Dy x Fe 2 could produce a large magnetostrictive strain at room temperature at x = 0.73 [1,2], so it is called a giant magnetostrictive material (giant magnetostrictive material, GMM). Based on a GMM, the giant magnetostrictive actuator (GMA) has a large output force and output displacement [3][4][5][6][7][8], so it can be used as a drive mechanism for direct drive electro-hydraulic servo valves. Goodfriend   A common GMA is mainly composed of a GMM rod, permanent magnets, coil, shell, ejector, and other parts. The alternating current will generate an alternating excitation magnetic field, which will generate an eddy current magnetic field inside the GMM rod. The excitation magnetic field, eddy current magnetic field, and bias magnetic field act together on the GMM rod, causing the magnetostrictive strain of the GMM rod and pushing the ejector rod to output a displacement. In addition, researchers find that the GMM rod will produce a larger output strain under a certain preload, so in a GMA, a disc spring is generally designed to provide the preload to the GMM rod.
It can be known from Figure 1 that when the material of the shell and the ejector of the GMA is selected with a low magnetic reluctance, the GMA can form an ideal closed magnetic circuit.
Most scholars [13] make the following assumptions during the magnetic field modeling phase: 1. The magnetic reluctance of the GMM is much higher than other materials such as the shell; 2. There is no air gap between all materials, or the air gap is very small. These two assumptions make the model ignore the magnetic circuit leakage and air gap leakage. In the forced oil-cooled GMA, the magnetic flux of the oil part is also ignored. This assumption is acceptable only when the GMA does not include an oil cooling system. However, when an oil cooling system is added to the GMA, a larger oil cavity will inevitably appear inside the GMA. As can be seen from reference [14], in GMA magnetic circuit modeling, the magnetic circuit law is similar to the circuit law: 1. In the magnetic circuit, the parallel magnetic segments have the same magnetomotive force; 2. In the magnetic circuit, the magnetic flux of the main circuit is the sum of the magnetic fluxes of the branches. Therefore, in the modeling of the magnetic circuit, the magnetic flux in the oil cavity should be taken into account, and the GMA cannot be simply regarded as a branchless magnetic circuit.
Because the structure of the GMA is axisymmetrical, its internal magnetic circuit is also symmetrical. To facilitate analysis, the GMA is cut along the axis and one side is taken for magnetic circuit analysis, as shown in Figure 2. Coils and permanent magnets work together to generate a magnetomotive force, and the magnetic flux will flow in the magnetic circuit along the sequence of the ejector ①-upper cover ⑥-shell ⑦ (oil cavity ⑧)-lower cover ⑨, where there are multiple places in series and parallel in the process. Because the structure of the GMA is axisymmetrical, its internal magnetic circuit is also symmetrical. To facilitate analysis, the GMA is cut along the axis and one side is taken for magnetic circuit analysis, as shown in Figure 2. Coils and permanent magnets work together to generate a magnetomotive force, and the magnetic flux will flow in the magnetic circuit along the sequence of the ejector 1 -upper cover 6 -shell 7 (oil cavity 8 )-lower cover 9 , where there are multiple places in series and parallel in the process.  In the SGMA, there are three sources of magnetomotive force: One is the excitation coil, the other is the permanent magnet, and the third is the eddy current magnetic field generated by the eddy current effect. When the three magnetic fields are superimposed, they act on the GMM rod together, causing the GMM rod to have a magnetostrictive effect.  In the SGMA, there are three sources of magnetomotive force: One is the excitation coil, the other is the permanent magnet, and the third is the eddy current magnetic field generated by the eddy current effect. When the three magnetic fields are superimposed, they act on the GMM rod together, causing the GMM rod to have a magnetostrictive effect.
The excitation magnetic field is generated by the current in the coil, and its magnetomotive force is where N is the number of coil turns and I is the excitation current. The magnetic field generated by the permanent magnet does not change, and its magnetomotive force U 2 is a certain value, which is determined by the nature of the permanent magnet itself. In the SGMA, as the permanent magnets and GMM rods are stacked alternately, the magnetic field generated by the permanent magnets is mainly absorbed by the GMM rods, and the magnetic leakage outside the magnetic circuit is small. Therefore, the magnetomotive force of a permanent magnet is approximately equal to the product of its equivalent coercive force H c and its length l 2 , that is, The eddy current magnetic field is generally manifested as hindering the change in the excitation magnetic field, so the eddy current magnetic field is the loss when the excitation magnetic field acts on the GMM rod, that is, the eddy current magnetomotive force U 3 and U 1 are related. In the SGMA, it is considered that U 3 only acts on GMM rods and does not affect other components in the magnetic circuit. Therefore, U 3 is regarded as the loss of U 1 and is not analyzed in the magnetic circuit. Thus, for the SGMA, In the magnetic circuit, the reluctance R Ωi of each component of the SGMA can be expressed as In Equation (4), l i is the length of each component (where l 2 is the sum of the lengths of the permanent magnets and l 3 is the sum of the lengths of the GMM rods), µ i is the permeability, and A i is the cross-sectional area.
Therefore, the total reluctance R Ω of the SGMA is Combining Equations (1), (2), and (5), it can be obtained that the average magnetic field strength H in the trunk road is Inside the coil, there is a gap between the GMM rod and the coil skeleton, so the average magnetic field strength H G distributed to the GMM rod and the permanent magnet can be expressed as Actuators 2020, 9, 37

of 19
The structure of the excitation coil and skeleton in the SGMA is shown in Figure 3. Set the coordinate axis with the center of the coil as the origin. The coil winding length is L C , the inner radius is R a , and the outer radius is R b .
Inside the coil, there is a gap between the GMM rod and the coil skeleton, so the average magnetic field strength HG distributed to the GMM rod and the permanent magnet can be expressed as ( ) The structure of the excitation coil and skeleton in the SGMA is shown in Figure 3. Set the coordinate axis with the center of the coil as the origin. The coil winding length is LC, the inner radius is Ra, and the outer radius is Rb. In Reference [14], the magnetic induction strength on the internal cross-section of the coil can be considered approximately equal everywhere. In a vacuum, the magnetic field strength excited at a point (0, Z) on the coil axis is Therefore, in the SGMA, the magnetic field intensity distributed to the stacked GMM rods is In Reference [14], the magnetic induction strength on the internal cross-section of the coil can be considered approximately equal everywhere. In a vacuum, the magnetic field strength excited at a point (0, Z) on the coil axis is Therefore, in the SGMA, the magnetic field intensity distributed to the stacked GMM rods is Through Equation (9), we establish the magnetic field model inside the SGMA and connect the input current with the magnetic field strength of the SGMA. In addition, the model describes the inhomogeneity of the internal magnetic field of the SGMA.

Magnetization Modeling of SGMA
When describing the relationship between the magnetic field strength and magnetization of a GMM, the Jiles-Atherton model is one of the currently used materialistic models [15]. The Jiles-Atherton model describes the obvious hysteresis and non-linear characteristics of the GMM magnetization process. From the perspective of physics, the GMM magnetization model is established in the form of calculus. In order to solve the problem where different scholars have different forms when quoting the Jiles-Atherton model, Qingyou et al. modified the Jiles-Atherton model of the plastic magnetization deformation of materials by discussing the physical meaning of different variables [16]. Guangming et al. unified the form of the Jiles-Atherton model based on the physical meaning and mathematical properties, and eliminated intermediate variables to simplify the Jiles-Atherton model [17,18].
According to Weiss molecular field theory, under the action of the external magnetic field H, the internal effective magnetic field H e of the GMM is In Equation (10), α is the molecular field parameter and M is the magnetization. In a GMM rod, the magnetization M comprises a reversible magnetization M rev and an irreversible magnetization M irr , that is, In Equation (11), the reversible magnetization M rev is related to the irreversible magnetization M irr and the hysteresis-free magnetization M an , that is, In Equation (12), c is a reversible coefficient, which is determined by the characteristics of the GMM, and represents the relationship between M rev and M irr .
From the law of conservation of energy, the static magnetization energy in the magnetization process is the difference between the non-hysteretic magnetization static magnetic energy and the hysteresis energy loss, that is, In Equation (13), k is a pinning coefficient, which indicates the strength of the hysteresis. δ and δ M are symbolic parameters, which are defined as From the theory of paramagnetic magnetization and the Langevin function, the non-hysteresis magnetization M an can be expressed as In Equation (15), M s is a saturation magnetization and a is a shape parameter. Equations (10)- (13), and (15) give the general expressions of the Jiles-Atherton model, that is Simplify Equation (16), remove the intermediate variables, and organize the Jiles-Atherton model into As the relationship between M and H is given by a nonlinear ordinary differential equation, an analytical solution cannot be obtained, so we use the Runge-Kutta method to solve its numerical solution. The fourth-order Runge-kutta formula is chosen here. Then, In Equation (19), h is the calculation step size, and M n and H n are the iteration values of M and H in the nth-step calculation, respectively.

Output Model of SGMA
According to the quadratic moment domain rotation model, we can establish an uncoupled magnetomechanical model of the GMM. The relationship between the magnetostrictive strain λ and the magnetization M of the GMM can be expressed as where λ S is the saturation magnetostrictive strain. The excitation coil generates an excitation magnetic field that is not uniformly distributed on the GMM rod, and the permanent magnet reduces the unevenness of the overall magnetic field of the SGMA. However, the permanent magnet increases the discontinuity of the local magnetic field intensity change of the GMM rod, that is, the magnetic field intensity change trend in the GMM rod on both sides of the permanent magnet is different. Therefore, in the dynamic analysis of the SGMA, it cannot be regarded as a whole, and the stacked GMM rods should be equivalent to n series mass-spring-damping systems, as shown in Figure 4. The excitation coil generates an excitation magnetic field that is not uniformly distributed on the GMM rod, and the permanent magnet reduces the unevenness of the overall magnetic field of the SGMA. However, the permanent magnet increases the discontinuity of the local magnetic field intensity change of the GMM rod, that is, the magnetic field intensity change trend in the GMM rod on both sides of the permanent magnet is different. Therefore, in the dynamic analysis of the SGMA, it cannot be regarded as a whole, and the stacked GMM rods should be equivalent to n series massspring-damping systems, as shown in As shown in Figure 4, m1 is the mass of the output jack, and m2-mn are the sum of the mass of a small piece of permanent magnet and the GMM rod, that is, In Equation (21), ρ1 and ρ2 are the mass densities of the permanent magnet and the GMM rod, respectively, and A2 is the cross-sectional area of the stacked GMM rod.
The multi-degree-of-freedom vibration equation is where M, C, and K are the mass matrix, damping matrix, and stiffness matrix, respectively, and Q is As shown in Figure 4, m 1 is the mass of the output jack, and m 2 -m n are the sum of the mass of a small piece of permanent magnet and the GMM rod, that is, In Equation (21), ρ 1 and ρ 2 are the mass densities of the permanent magnet and the GMM rod, respectively, and A 2 is the cross-sectional area of the stacked GMM rod.
The multi-degree-of-freedom vibration equation is where M, C, and K are the mass matrix, damping matrix, and stiffness matrix, respectively, and Q is the driving force vector.
Because the stacked GMM rods use homogeneous materials, and the performance of the same material in the axial direction is not much different, it can be considered that its elastic modulus and structural damping are uniform. Therefore, in a multi-degree-of-freedom model, In Figure 4, the displacement of each element of the multi-degree-of-freedom vibration system mainly manifests as the deformation of the spring and the damper. In essence, it is the deformation of the GMM rod.
Set the average strain of the i-th stage GMM to be λ i , then its magnetostrictive force is where E N is the GMM elastic modulus. The permanent magnet hardly generates magnetostrictive strain, so its force comes from the strain of the GMM rod on both its sides. Therefore, in a multi-degree-of-freedom vibration system, the driving force of m i is For a SGMA, both the output force F and the output displacement x ultimately appear as the displacement and force of the ejector m 1 , so Through the output model of the multi-degree-of-freedom system, we established the relationship between the magnetic field and the output displacement (and force). Although this relationship is established by a discrete model, because we only focus on the force and displacement of the ejector, a continuous output is obtained.

Loss Analysis of SGMA
Because of its low resistivity, a GMM has an eddy current effect in an alternating magnetic field, which will generate an eddy magnetic field. The eddy magnetic field is opposite to the excitation magnetic field, thus hindering the effect of the excitation magnetic field on the GMM rod. The eddy current magnetic field is not uniformly distributed inside the GMM rod and reaches the maximum at the center of the GMM rod, so the skin effect occurs. The effective magnetic field can reach the radial depth inside the material, and the skin depth δ can be expressed as In Equation (27), f is the frequency of the excitation magnetic field, σ is the conductivity of the GMM, µ 0 is the vacuum permeability, and µ r is the relative permeability of the material.
The external environment in which the GMM rod is located, mainly the magnetic field strength and the pre-stressing force, easily affects its relative permeability. The relationship between magnetostrictive strain and magnetic field strength under different prestresses is shown in Figure 5a, and the measured curve of the relative permeability of the GMM rod with the applied magnetic field strength under the preload of 10 MPa is shown in Figure 5b [19]. It can be seen from Figure 5 that when the magnetic field strength changes, the relative permeability of the GMM changes greatly. Taking μr = 7, 9, and 11, and the resistivity of the GMM 1/σ = 6 × 10 -7 Ω·m, the skin depth is shown in Figure 6. It can be seen from Figure 5 that when the magnetic field strength changes, the relative permeability of the GMM changes greatly. Taking µ r = 7, 9, and 11, and the resistivity of the GMM 1/σ = 6 × 10 -7 Ω·m, the skin depth is shown in Figure 6. (b) It can be seen from Figure 5 that when the magnetic field strength changes, the relative permeability of the GMM changes greatly. Taking μr = 7, 9, and 11, and the resistivity of the GMM 1/σ = 6 × 10 -7 Ω·m, the skin depth is shown in Figure 6. It can be seen from Figure 6 that the skin depth of the GMM decreases as the frequency of the applied magnetic field increases. When the frequency of the applied magnetic field reaches the order It can be seen from Figure 6 that the skin depth of the GMM decreases as the frequency of the applied magnetic field increases. When the frequency of the applied magnetic field reaches the order of 1000 Hz, the skin depth of the GMM rod is less than 5 mm, which seriously affects the magnetization effect of the GMM rod and the output accuracy of the SGMA.
Another key parameter of the GMM is the eddy current cutoff frequency, which is the minimum frequency at which the skin effect can be ignored. This frequency is expressed as [20] f c = 2σ πd 2 µ r µ 0 .
In Equation (28), d = 2r is the diameter of the GMM rod. Selecting the same parameter settings as in Figure 6, we obtain the relationship between the diameter of the GMM rod and the eddy current cutoff frequency, as shown in Figure 7.
Actuators 2020, 9, x FOR PEER REVIEW 11 of 20 of 1000 Hz, the skin depth of the GMM rod is less than 5 mm, which seriously affects the magnetization effect of the GMM rod and the output accuracy of the SGMA.
Another key parameter of the GMM is the eddy current cutoff frequency, which is the minimum frequency at which the skin effect can be ignored. This frequency is expressed as [20] In Equation (28), d = 2r is the diameter of the GMM rod. Selecting the same parameter settings as in Figure 6, we obtain the relationship between the diameter of the GMM rod and the eddy current cutoff frequency, as shown in Figure 7. Among different magnetic materials, when the operating frequency is the same, the lower the eddy current cut-off frequency fc, the greater the eddy current loss. It can be seen from Figure 7 that the diameter of the GMM has a significant negative correlation with the cutoff frequency of the Among different magnetic materials, when the operating frequency is the same, the lower the eddy current cut-off frequency f c , the greater the eddy current loss. It can be seen from Figure 7 that the diameter of the GMM has a significant negative correlation with the cutoff frequency of the material, and in the large-diameter GMM rod, the eddy current effect has a great influence. When the diameter of the GMM reaches 10 mm, the upper limit frequency of the applied magnetic field is less than 500 Hz, which severely limits the response speed of the SGMA.
Known from Maxwell's equations, In Equation (29), H is the magnetic field strength, J is the current density, E is the electric field strength, B is the magnetic induction strength, and σ is the conductivity.
From Equation (29), the propagation formula of H can be obtained For rod-shaped GMM materials, take any point inside it, and set the distance from this point to the axis of the GMM rod as x. At time t, when the external magnetic field H = H 0 e jωt , the internal magnetic field of the GMM rod can be expressed as where u x = x −jωσµ r µ 0 = xγ = x −jγ = x − j δ . Considering Equation (31) as the zero-order Bessel function of the imaginary argument, combined with the Kelvin function, the solution can be expressed as [21] H( where J 0 is the Bessel function of the first kind, u r = u x | x=r , and . Ber is the real part of the Kelvin function, Bei is the imaginary part of the Kelvin function, and J 0 − jγ x = Ber(γ x) + jBei(γ x).
It can be known from Figure 8 that when the frequency of the applied magnetic field increases, the relative magnetic field strength of the GMM rod decreases, that is, there is a part in the GMM that cannot be completely penetrated by the magnetic field, and the closer to the center of the rod, the greater the attenuation of the magnetic field. Moreover, GMM rods with different rod diameters have different blocking effects on high-frequency magnetic fields. The thicker the GMM rod (i.e., the larger the radius), the stronger the attenuation at the same frequency.
Take the parameters of the SGMA as follows: GMM resistivity 1/σ = 6 × 10 -7 Ω·m; relative permeability μr = 11; radius r = 4, 10 mm; frequencies are 100, 500, 1000, 2000, or 10,000 Hz. According to Equation (32), the internal magnetic field strength of the GMM rod is shown in Figure 8. It can be known from Figure 8 that when the frequency of the applied magnetic field increases, the relative magnetic field strength of the GMM rod decreases, that is, there is a part in the GMM that cannot be completely penetrated by the magnetic field, and the closer to the center of the rod, the greater the attenuation of the magnetic field. Moreover, GMM rods with different rod diameters have different blocking effects on high-frequency magnetic fields. The thicker the GMM rod (i.e., the larger the radius), the stronger the attenuation at the same frequency.
When the GMM rod is in an alternating magnetic field, due to the eddy current effect, the induced electromotive force generated in the rod can be expressed as The internal resistance of the GMM rod can be expressed as When the GMM rod is in an alternating magnetic field, due to the eddy current effect, the induced electromotive force generated in the rod can be expressed as The internal resistance of the GMM rod can be expressed as Therefore, the eddy current loss of the GMM is In the alternating magnetic field, due to the inherent hysteresis effect of the GMM, the magnetic induction intensity will be one phase ϕ behind the applied magnetic field, that is, Therefore, the permeability of the GMM in an alternating magnetic field can be expressed as where Due to the imaginary part of the relative permeability, the GMM also has hysteresis loss in the applied alternating magnetic field. The loss per unit time can be expressed as Simplified, Actuators 2020, 9, 37 13 of 19 From Equations (37) and (39), we can know that the hysteresis loss of the GMM is In the SGMA, the inductance of a coil can be expressed as [22] where k c is the size coefficient. The impedance of the coil is Then, under the alternating current, the heat loss of the coil is where Z = R 2 ct + L 2 ct is the coil impedance and R ct is the coil resistance. From Equations (35), (39), and (43), the total loss p in the SGMA is The output model describes the process of converting electrical energy into mechanical energy, and the loss model describes the process of converting electrical energy into internal energy. These two processes occur at the same time. Of course, changes in internal energy will also cause changes in material properties, so we need to reduce the conversion of electrical energy into internal energy or reduce the impact of internal energy on output.

Design of the Cooling System
There is a large loss in the SGMA, and a large amount of energy is converted into internal energy. Due to the compact structure of the SGMA itself, internal energy is generated inside the SGMA and the heat exchange with the outside is not smooth. Therefore, the internal energy generated by the loss cannot be lost in time. A large amount of internal energy can accumulate inside the SGMA, causing its temperature rise. In addition, the GMM is a heat-sensitive material, so when the temperature changes, the GMM magnetostrictive strain is as shown in Figure 9. The output model describes the process of converting electrical energy into mechanical energy, and the loss model describes the process of converting electrical energy into internal energy. These two processes occur at the same time. Of course, changes in internal energy will also cause changes in material properties, so we need to reduce the conversion of electrical energy into internal energy or reduce the impact of internal energy on output.

Design of the Cooling System
There is a large loss in the SGMA, and a large amount of energy is converted into internal energy. Due to the compact structure of the SGMA itself, internal energy is generated inside the SGMA and the heat exchange with the outside is not smooth. Therefore, the internal energy generated by the loss cannot be lost in time. A large amount of internal energy can accumulate inside the SGMA, causing its temperature rise. In addition, the GMM is a heat-sensitive material, so when the temperature changes, the GMM magnetostrictive strain is as shown in Figure 9. It can be seen from Figure 9 that when the temperature changes, the magnetostrictive strain value of the GMM changes greatly.
In addition, a sharp rise in the internal temperature of the SGMA will also cause thermal expansion of the GMM rod. When the temperature changes, the thermal deformation of the GMM is It can be seen from Figure 9 that when the temperature changes, the magnetostrictive strain value of the GMM changes greatly.
In addition, a sharp rise in the internal temperature of the SGMA will also cause thermal expansion of the GMM rod. When the temperature changes, the thermal deformation of the GMM is where λ t = 12 × 10 −6 K −1 is the linear expansion coefficient of the GMM, l is the length of the GMM rod, and ∆t is the temperature change. Therefore, temperature has a great influence on the SGMA. In the modeling of the SGMA, the influence of temperature must be considered. However, the heating source, heating mechanism, and heat exchange process in the SGMA are relatively complicated. When applying the non-steady-state heat conduction differential equation and Fourier's theorem, many parameters are not easy to obtain accurately, so we use a simplified design method [23].
From the SGMA itself, in order to suppress the effect of temperature changes on its output, known from Figure 9, it should be stabilized to work at about 40 • C. Therefore, in the design process, the parameters that need attention are only the size of the SGMA and the flow rate of the cooling oil. The size parameters of the SGMA are determined by its output characteristics, and it is only necessary to pay attention to the flow rate of the cooling oil in terms of temperature characteristics.
Assume that, in a SGMA, the heat transferred from the surface of each component is proportional to its surface area. According to the principle of pyrology, the flow rate of cooling oil is at least where c is the specific heat capacity of the cooling oil.

Experimental Verification of Oil-Cooled SGMA
According to the foregoing analysis, the prototype parameters of the SGMA are as Table 1: The prototype of the SGMA is shown in Figure 10.
Specific heat capacity of cooling oil c/kJ/(kg·K) 0.7-0.8 Density of cooling oil /g/cm 3 0.857 Thermal conductivity of cooling oil /W/mc 0.11 Flow rate of cooling oil q/mL/s 15 The diameter of the wire /mm 1.20 The prototype of the SGMA is shown in Figure 10. In order to test the relevant performance of the SGMA prototype, we set-up a SGMA test system, as shown in Figure 11  In order to test the relevant performance of the SGMA prototype, we set-up a SGMA test system, as shown in Figure 11.
Actuators 2020, 9, x FOR PEER REVIEW 16 of 20 Figure 11. Test system of SGMA prototype.
The test system uses the Rigol-DG1022U signal generator (Shenzhen Keruijie Technology Co., Ltd., Shenzhen, China) to generate the excitation signal, which is amplified by the GF800W power amplifier to drive the SGMA. The cooling oil is pumped into the SGMA by an oil pump to stabilize the temperature of the SGMA. The IT6932A programmable voltage source provides Microtrak TM 3-LTS-025-02 laser (MTI Instruments, Inc., Washington, USA) displacement sensor power, so that the laser displacement sensor can accurately measure the displacement of the SGMA. The system uses the TES-1310 temperature sensor (TES Electrical Electronic Corp., Taiwan, China) to monitor the actuator temperature, the XMK-7 temperature sensor (Mingxing Refrigeration Parts Factory, Yuyao, Figure 11. Test system of SGMA prototype. The test system uses the Rigol-DG1022U signal generator (Shenzhen Keruijie Technology Co., Ltd., Shenzhen, China) to generate the excitation signal, which is amplified by the GF800W power amplifier to drive the SGMA. The cooling oil is pumped into the SGMA by an oil pump to stabilize the temperature of the SGMA. The IT6932A programmable voltage source provides Microtrak TM 3-LTS-025-02 laser (MTI Instruments, Inc., Washington, USA) displacement sensor power, so that the laser displacement sensor can accurately measure the displacement of the SGMA. The system uses the TES-1310 temperature sensor (TES Electrical Electronic Corp., Taiwan, China) to monitor the actuator temperature, the XMK-7 temperature sensor (Mingxing Refrigeration Parts Factory, Yuyao, China) to monitor the fuel tank oil temperature, the PICO-TA189 current clamp (PicoTechnology Inc, Cambridgeshire, UK) to monitor the input current value, and the PICOscope digital oscilloscope to collect test data. The cooling system is shown in Figure 12.
Ltd., Shenzhen, China) to generate the excitation signal, which is amplified by the GF800W power amplifier to drive the SGMA. The cooling oil is pumped into the SGMA by an oil pump to stabilize the temperature of the SGMA. The IT6932A programmable voltage source provides Microtrak TM 3-LTS-025-02 laser (MTI Instruments, Inc., Washington, USA) displacement sensor power, so that the laser displacement sensor can accurately measure the displacement of the SGMA. The system uses the TES-1310 temperature sensor (TES Electrical Electronic Corp., Taiwan, China) to monitor the actuator temperature, the XMK-7 temperature sensor (Mingxing Refrigeration Parts Factory, Yuyao, China) to monitor the fuel tank oil temperature, the PICO-TA189 current clamp (PicoTechnology Inc, Cambridgeshire, UK) to monitor the input current value, and the PICOscope digital oscilloscope to collect test data. The cooling system is shown in Figure 12. Pass a 1, 10, 100, and 1000 Hz sine wave, square wave, and triangular wave through the SGMA. The corresponding experimental results are as follows: It can be seen from Figure 13 that the SGMA responds well to signals of different waveforms, but due to the influence of the response speed of the coil and the nature of the material itself, the displacement waveform has a deviation compared to the input waveform. Pass a 1, 10, 100, and 1000 Hz sine wave, square wave, and triangular wave through the SGMA. The corresponding experimental results are as follows: It can be seen from Figure 13 that the SGMA responds well to signals of different waveforms, but due to the influence of the response speed of the coil and the nature of the material itself, the displacement waveform has a deviation compared to the input waveform.
The output waveform has a lag relative to the input, and as the excitation frequency increases, the lag phase difference increases. In addition, the output displacement of the SGMA decreases with increasing frequency. On the waveform, the experimental results are similar to the theoretical calculation results, and the phase difference between them is small, indicating that the modeling theory can be used to describe the principle of the SGMA. As the frequency increases, the error between the theoretical and experimental results increases, and the loss from the actual output is significantly greater than the theoretical calculation, showing that the modeling theory does not accurately describe the high-frequency motion of the SGMA. This is because the calculation of high-frequency loss in the output modeling theory mainly comes from the J-A model, and does not consider the case of incomplete magnetization of the GMM rod due to the eddy current and skin effect under high-frequency action.
Sine waves of 50, 100, and 500 Hz are passed into the SGMA, and the internal temperature change of the SGMA is monitored by a temperature sensor. The results are shown in Figure 14.
It can be seen from Figure 14 that when cooling oil is introduced into the SGMA, the internal temperature of the SGMA slowly rises and remains stable after reaching a certain temperature, but the stable temperature is different at different frequencies, and the higher the frequency, the higher the temperature. The experimental results show that the cooling system can indeed keep the temperature of the SGMA stable and restrain the thermal deformation of the material and the change in the magnetostrictive coefficient caused by the temperature.   consider the case of incomplete magnetization of the GMM rod due to the eddy current and skin effect under high-frequency action.
Sine waves of 50, 100, and 500 Hz are passed into the SGMA, and the internal temperature change of the SGMA is monitored by a temperature sensor. The results are shown in Figure 14. It can be seen from Figure 14 that when cooling oil is introduced into the SGMA, the internal temperature of the SGMA slowly rises and remains stable after reaching a certain temperature, but the stable temperature is different at different frequencies, and the higher the frequency, the higher the temperature. The experimental results show that the cooling system can indeed keep the temperature of the SGMA stable and restrain the thermal deformation of the material and the change in the magnetostrictive coefficient caused by the temperature.

Conclusions
This paper analyzes the magnetic field distribution of a stacked giant magnetostrictive actuator, and it establishes the magnetization model and output model of a SGMA, analyzes various losses occurring in the work of the SGMA, and designs the cooling system of the SGMA. Based on theoretical analysis, we make a prototype of the SGMA, build a corresponding test system, conduct related experiments, and obtain the following conclusions: 1. Inside the SGMA, the magnetic field distribution is not uniform. Due to the existence of the cooling oil cavity, there are multiple series and parallel situations in the magnetic circuit, and the magnetic field generated by the excitation coil cannot be completely concentrated on the GMM rod.
2. On the GMM rod, the magnetic field distribution is also not uniform. In the axial direction, the distribution of the magnetic field on the GMM rod is related to its axial coordinates, and in the radial direction, it is related to its radius and the frequency of the applied magnetic field.
3. There are various losses in the SGMA, including eddy current loss, hysteresis loss, and coil loss. These three losses cause a large amount of energy to be dissipated as internal energy and affect the output accuracy of the SGMA.

Conclusions
This paper analyzes the magnetic field distribution of a stacked giant magnetostrictive actuator, and it establishes the magnetization model and output model of a SGMA, analyzes various losses occurring in the work of the SGMA, and designs the cooling system of the SGMA. Based on theoretical analysis, we make a prototype of the SGMA, build a corresponding test system, conduct related experiments, and obtain the following conclusions:

1.
Inside the SGMA, the magnetic field distribution is not uniform. Due to the existence of the cooling oil cavity, there are multiple series and parallel situations in the magnetic circuit, and the magnetic field generated by the excitation coil cannot be completely concentrated on the GMM rod.

2.
On the GMM rod, the magnetic field distribution is also not uniform. In the axial direction, the distribution of the magnetic field on the GMM rod is related to its axial coordinates, and in the radial direction, it is related to its radius and the frequency of the applied magnetic field.

3.
There are various losses in the SGMA, including eddy current loss, hysteresis loss, and coil loss. These three losses cause a large amount of energy to be dissipated as internal energy and affect the output accuracy of the SGMA.

4.
Considering the loss of heat generation, in order to reduce the influence of temperature on accuracy, this paper designs a SGMA cooling system. Through analysis, this paper uses a simplified calculation method and obtains the fact that the flow rate of the cooling oil is related to the size and loss of SGMA. 5.
In this paper, a complete modeling theory of SGMA is established, and all theoretical models from current to output and loss are obtained. Not only the output but also the loss are analyzed, and an active temperature control method is proposed. The SGMA designed in this paper has a fast response speed, stable output waveform, and small error under the conditions of laboratory and low-frequency excitation signals. When the cooling system is installed, the temperature rise of the SGMA can be effectively controlled, which can meet the application requirements of direct-drive electro-hydraulic servo valves.

Conflicts of Interest:
The authors declare no conflict of interest.