Study on the Magnetic-machine Coupling Characteristics of Giant Magnetostrictive Actuator Based on the Free Energy Hysteresis Characteristics

A giant magnetostrictive actuator presents advantages such as large strain, high precision, and quick response. It is a hotly debated research topic in the field of micro drivers; however, the nonlinear intrinsic relationship between its output and input signals make it difficult to construct its nonlinear eigen model in the process of its practical application. Therefore, the motivation of this paper is to study the nonlinear magnetic–mechanical coupling characteristics of the giant magnetostrictive actuator, which is driven by free energy hysteresis characteristics. The nonlinear magnetic–mechanical coupling model under the weak form solution is deduced from the basic electromagnetic and mechanical theories, based on the distribution law of the axial magnetic field simulation, carried out to analyze the output displacement characteristics of the giant magnetostrictive actuator under preload. Experimental characterization of the device is also studied in the built experiment setup. Research results show that the experimental results coincide well with the simulation results, which show that the designed magnetic circuit for the giant magnetostrictive actuator is correct, and the coupling model of magnetic and machine of the giant magnetostrictive actuator based on the free energy hysteresis characteristics is reasonable.

Based on the giant magnetostrictive (GMS) effect, a novel giant magnetostrictive actuator with high actuation power is designed and implemented by Mingzhang Luo [7]. This design enables the generation of stress waves with high energy, and the focusing of the generated stress waves on the test object, it can be used in the quality assessment of rock bolt-reinforced structures and other nondestructive testing and evaluation applications that require high-power stress wave generation. Zhenyuan Jia et al. [8] used a giant magnetostrictive actuator (GMA) made of giant magnetostrictive material (GMM) to apply a magnetic field to strain it by using its magnetostrictive effect and to produce output force or displacement using its longitudinal deformation, with large strain, high precision, fast response, high reliability, and so on. It quickly became a research hotspot in the field of micro-drive. Yuanyuan Yang and his college [9] present an induced voltage linear extraction None of the above methods establishes a model based on the characteristics of the giant magnetostrictive material itself, and the giant magnetostrictive material is a ferromagnetic material, its essential feature is that there are intrinsic nonlinearities and hysteresis characteristics, so it is a reasonable method to study the intrinsic relationship between the output and input of the actuator based on the intrinsic nonlinearity and hysteresis characteristics of the giant magnetostrictive material.
In this paper, the hysteresis characteristics and nonlinearity of giant magnetostrictive actuators are studied in depth by the free energy hysteresis model. The three-dimensional magnetic machine coupling model of the giant magnetostrictive actuator is established, and a three-dimensional nonlinear coupled model of the giant magnetostrictive actuator is established according to the free-energy hysteresis model that used the material eigen model combined with the three-dimensional coupled model, then the finite element simulation software is used for analysis and calculation to obtain the magnetic induction and strain increments. Following this, the total strain output is obtained, and the simulation results were verified by experiments.

The Fundamental Theory of the Magnetic Field Excited by the Winding for GMA
For the giant magnetostrictive actuator system, Maxwell's equations [22][23][24] are used to solve the distribution of the magnetic field excited by the winding. The solving variables are set as the vector magnetic potential (A x , A y , A z ) in the magnetic field. According to the differential Maxwell equations of the magnetic field and the principle of current continuity, the driving frequency of the magnetic field excited by the winding is lower than 30 MHz. The driving current source is chosen to be a constant current source so that the magnetic field is formed to be a constant magnetic field. Therefore, the electric displacement field vector D can be neglected, then ∂D ∂t = 0. The total current density consists of two parts: source current density J S and applied current density J E , which is caused by eddy current effect, then the divergence of magnetic field intensity is as follows: In this equation, H denotes the magnetic field intensity vector. Electric field intensity E can be expressed by magnetic vector potential: − ∇φ (2) where φ is a scalar potential function. Then the applied current density caused by eddy current effect is expressed as (3): where σ is electric conductivity. Then the magnetic flux density B is expressed by magnetic vector potential:

The Fundamental Theory of the Mechanical Field of GMA
According to Newton's Second Law, the motion of a giant magnetostrictive rod can be described as follows: where the displacement u = u x u y u z T , the stress T =    T 13 = T 31 , T 23 = T 32 , and f B is body force. When damping is considered, and Navier-stocks equation [23] is combined to derive Formula (6): (6) where c is the viscous damping coefficient, and m is quality of the system. The solving variable of the mechanical field is set as displacement vector (u x , u y , u z ), then the relationship between the elastic strain S and the displacement can be expressed as follows: If it is substituted into Equation (7), Formula (8) is obtained as follows:

The Conversion of the Weak Form Solution
(1) Einstein notation According to Einstein Notation [24], tensor ∈ ijk is a Levi-Civita signal, which is defined as follows: The cross product of the two variables can be expressed as follows: The curl of the variables is shown as follows: Use Einstein Notation to mark Equations (1) and (3): Substituting Equation (11) into the above two equations: Use Einstein Notation to mark Equation (6) and expand the formula: Use weighted residual method to integrate Equations (13) and (15): where ψ i and φ i are weight function, V B is magnetic field distribution region, and V u is action area of mechanical field. Take partial integration for the first term on the left of the Equation (16), then Formula (17) is derived as: Define the divergence V ∇·FdV = S F·ndS, according to divergence theorem, the first term on the right of the Equation (17) can be expressed as follows: Substitute the above Equation into equation (17): Then the Equation (19) can be rearranged as follows: Substitute the above equation into Equation (16): According to Galerkin Method [22,24], weight function is equal to basis function, where minimal variables A i , u i are used. The weight function is expanded to the following approximate expression: Equation (21) can be expressed as follows: This is the weak solution equation for solving the magnetic and mechanical field variables of GMA system under the Einstein Notation.
(2) Matrix notation Equation (23) can be expressed by matrix notation as Equation (24): The equations for the weak solution formally express the balance of virtual work both inside and outside the system. According to the kinematics relationship, the variation variables δB = ∇ × δA, δS = ∇δu are introduced to calculate H and T. Because the surface traction at the mechanical field boundaries is t = Tn, the tangential component of the magnetic field boundaries is H T = H × n, then the virtual work of the system can be expressed as follows:

Three-Dimensional Finite Element Discrimination of Virtual Work Model for the System
Using the finite element method, the solvability domain is discretized into a finite number of units (the number of element is assumed to be e), the weak form solution Equation (24) is solved on the element. The solving variables, the vector potential energy, offset variables of every element should be calculated by node value interpolation. Thus, their interpolation will be written into matrix form using the interpolation and the shape functions to calculate the variable node value.
Interpolation and integration are calculated at the local coordinate ξ, which is combined with dx and dξ through the Jacobian matrix J e , where dx = J e dξ, dV = det(J e )dξ 1 dξ 2 dξ 3 = J e dξ 1 dξ 2 dξ 3 . For the linear interpolation of geometry structure, the variables A e and u e of each element can be interpolated by node values q A e and q u e : According to the Galerkin method [23,24], weight function is equal to basis function, then the shape function of minimal variables is written as follows: Because A is a three-dimensional vector, the interpolating matrix N A presents three lines, N A n nods in each element (N A n is depended on the order of the element), and every nod can be described as a three-dimensional vector, which is relative with the nodal value A e . So the interpolating matrix N A consists of N A q columns, where N A q = 3N A n . The dimension of the vector q A e is N A q , and all three components in each node are A e . Because the displacement variables does not require the same shape function, the available number of columns of N u can be expressed by N u q , which depended on the number of node N u n . From the above information, it is observed that the total degree of freedom of each element is N q = N u q + N A q . In this paper, the typical shape functions are mainly tetrahedron elements of linear functions or quadratic Lagrange function [24]. Assuming that vector potential energy and displacement variables are expressed in four-node tetrahedral element (as shown in Figure 1), then N A = N q = N. Matrix shape function N is consisted of a Lagrange function: where, for node i, when j = i, N j = 1, otherwise, when j = i, N j = 0. According to the shape function (28), in the global coordinate system, x 1 can be interpolated with the node value x 1,n .
Sensors 2018, 18, x FOR PEER REVIEW 8 of 21  , the vector potential energy can be simplified to the value of node 2. It is the same for the other coordinates. For other points in the tetrahedron, generally, the linear interpolation of node values is used for the shape function; for instance, local coordinates is written as In the finite element model, the node values of the vector potential energy and the displacement variables are unknown. The virtual vector potential energy and the displacement variables can be set Because the four-node tetrahedral element in the finite elements is used, the dimension of vector q A e is 12, and the first three dimensions at node 1 are the three components of A e , which is followed by the three components at nodes 2, 3, and 4. The matrix of shape function is shown as Formula (30).
When local coordinate ξ = (1, 0, 0) is located at node 1, A e = A e,1 , the vector potential energy can be simplified to the value of node 1. Similarly, when local coordinate ξ = (0, 1, 0) is located at node 2, A e = A e,2 , the vector potential energy can be simplified to the value of node 2. It is the same for the other coordinates. For other points in the tetrahedron, generally, the linear interpolation of node values is used for the shape function; for instance, local coordinates is written as ξ = (1/2, 1/2, 0), A e = 1/2(A e,1 + A e,2 ).
In the finite element model, the node values of the vector potential energy and the displacement variables are unknown. The virtual vector potential energy and the displacement variables can be set to be any values. In order to show finite element discretization in the virtual work expression (25), the magnetic flux density and the strain can be expressed in terms of vector potential energy and the displacement variables.
The elements in the matrices C e and G e are the derivatives of the local coordinate system ξ about the global coordinate system x, which is assumed to be the discrete form of curl and gradient. For linear elements, Jacobean matrix J e is a constant matrix. Then, matrix C e and G e does not depend on ξ.
Assuming that magnetic field be divided into N A elements, the mechanical field is discretized into N u elements. The virtual work balance expression can be shown as follows: Because the influence of gravity and Lorentz force in the magnetostrictive device can be neglected, body force could be ignored [25]. Where subscript b denotes the number of boundary elements, N A S and N u S represents the number of boundary elements of the magnetic field and the mechanical field, respectively. Integral ∆ J e d∆ means element integration in the global coordinates. Then, the element volume can be described as (33).
The integral ∆ S J S,b d∆ S means boundary element integration in the global coordinate, which represents the element area.

Magneto Mechanical Coupling Nonlinear Model
Magneto mechanical coupling relationship [26] of the magnetostrictive actuator system can be described by the piezomagnetic equation: where µ T is permeability tensor at constant stress; S H is flexibility coefficient tensor at constant magnetic field intensity; and d and S are the piezomagnetic tensor and strain, respectively.
The above equations show that the total magnetic flux density consisted of two parts. The first part shown in the first term denotes the magnetic flux density µ T H that is induced by the applied magnetic field, and the second terms represents magnetic flux density dT that is caused by the application of mechanical stress. Similarly, the total strain is also formed with two terms, where the first term is the magnetic strain d T H that is caused by applied magnetic field, while the second one denotes the elastic strain S H T, which is caused by the applied stress. Inside the giant magnetostrictive materials, its strain is immediately relevant with the magnetic intensity, stress state, and its material properties. The actuator's output displacement and output force are the result of mutual coupling of the magnetic and elastic fields.
The above eigenvalue equations can be expressed in incremental form as follows: The increments of the magnetic field and the stress field are combined with the solution result calculated by the finite element through the Equation (36). Solving from the known initial solution, the pre-compression t, surface magnetic field strength H T , and source current density J S are all in incremental form. Replacing the discrete variables in the Equation (31) with Formula (36) and (37) can be derived as follows: The above equation can be substituted into Equation (32), and the symbols of the integral variables in the equation are defined as follows: According to the definition of the symbols shown above, the virtual work balance expression of the magnetic field and the mechanical field can be written as follows: Reorganizing the above equation with variational principle, it can be written into a matrix form: As it can be seen from the above equation, the mass matrix is singular and only the mass of mechanical field is included. As the giant magnetostrictive actuator is excited by a frequency lower than MHz under normal working conditions, the electric displacement vector can be neglected. As a result, the second derivative term of the vector potential energy will be zero. The damping matrix is mainly derived from the material's internal damping d u e and the eddy current damping d A e (the eddy current is depended on the conductivity). In the stiffness matrix, k A e depends mainly on the magnetic permeability and it characterizes the capability of the magnetic excitation system. The coupling matrix k u,A e characterizes the capability of the magnetic field converted into mechanical energy by the mechanical traction vector f u b under the surface traction and capability of the mechanical energy converted into the magnetic energy at the magnetic field load vector f J e + f A b . The incremental form of the coefficient relationship (36) can be solved by the free energy model. Its constitutive model can be written as a function of magnetic flux density B and the mechanical strain S about the magnetic field intensity H and the mechanical stress T, and Equation (35) is derivative. The coefficients can be expressed in the form of the binary function's derivative: Coefficients in above formula can be written into matrix form shown as (42): For the giant magnetostrictive materials, the differences of the lattice orientation, the driving magnetic field, and the stress will cause the changes of the parameter matrix values. It is assumed that the pre-pressure of the actuator is large enough to dominate the lattice anisotropy so that the changes of the parameter matrix caused by anisotropy ignoring can be ignored. Then, the magnetostrictive material is considered to be isotropous. According to equation of M = B/µ 0 − H and the scalar free energy model, the above matrix values can be solved respectively by M.

Calculation of the Coupling Model of Giant Magnetostrictive Actuator (GMA)
A giant magnetostrictive actuator is a new type of precision actuator, which is based on the magnetostriction deformation effect of magnetostrictive material rod when the input magnetic field changes. Giant magnetostrictive actuators is designed according to the properties of giant magnetostrictive materials. In this paper, we perform finite element simulation analysis on the established magnetomechanical coupling model.

The Structure and Working Principle of the Giant Magnetostrictive Actuator
The structure of giant magnetostrictive actuator is shown in Figure 2; it is mainly composed of giant magnetostrictive material rod, coil, coil bobbin outer tube, output shaft, front-end cover, rear end cover, shell, inlet and outlet, preloading spring, and fixing plate. In order to reduce the magnetic flux leakage and rise the magnetic field intensity of the center axis of the giant magnetostrictive material rod, the closed magnetic circuit consists of the giant magnetostrictive material rod, the output shaft, the shell, and the rear end cover. The output shaft, the front end cover, and the preloading spring form a preloading mechanism together, and pre-pressure is given to the giant magnetostrictive material rod in the pre-experiment period so that the giant magnetostrictive material rod obtains a relatively larger axial output strain under the action of the magnetic field. Then, greater output displacement and force are obtained. At the same time, the giant magnetostrictive rod is exerted with a certain pre-pressure to keep it working on a state of compression and prevent material's fragmentation caused by its high brittleness. The coil bobbin outer tube composes a water-cooled cavity. The water inlet and outlet are used as the lead-in and lead-out holes to allow the overall system of the actuator to carry out water circulation, thereby keeping the temperature of the GMM rod within a certain range to suppress thermal deformation. The working principle of the giant magnetostrictive actuator is that the driving coil generates a magnetic field under the effect of the direct current signal, which causes the giant magnetostrictive rod to induce magnetostrictive deformation to move the output shaft and realize the output of the force and the displacement. Using the method of adjusting the driving DC signal can change the amplitude of the magnetic field. By this way, we can obtain the output of the displacement and the force in different sizes.

Finite Element Simulation of the Coupling Model of Giant Magnetostrictive Actuator (GMA)
In this study, the MATLAB2009a software demo (which is a commercial mathematics software produced by MathWorks, MA, USA) is used to program the model, and a demo version of COMSOL 4.0a software demo (which eveloped by the COMSOL group who founded in Stockholm, Sweden) [17] is used to simulate the coupling model of giant magnetostrictive actuator. The analysis process includes the following: Firstly, COMSOL finite element simulation software is used to construct the geometric model and mesh the actuator, and the coefficient matrix is iterated through the MATLAB software. Secondly, both two are invoked to complete the coupled analysis of the material model and the dynamic model. In theory, the magnetic field of the coil fills the whole space. In the finite element simulation, a larger closed cuboid air field is built to substitute the space. The giant magnetostrictive actuator is placed in the air field, and the boundary of air field is set to be magnetic insulation.
The diameter of the GMM rod is 10 mm, and its length is 450 mm. The material of the shell is chosen to be 45 steel. The diameter of single wire of driving coil is 1 mm, and its height is 22 mm. The total number turns of winding coil is 9570 turns. The solution process of the model is shown in Figure 3: Firstly, the current values are inputted and the model is meshed. Secondly, the increment of the initial magnetic field strength can be obtained by the input of initial coefficient matrix. Thirdly, the coefficient matrix is updated by solving the model, where the increments of magnetic induction and strain should be calculated so that further new incremental magnetic field intensities can be obtained. The cycle is repeated until the current waveform function is terminated. The working principle of the giant magnetostrictive actuator is that the driving coil generates a magnetic field under the effect of the direct current signal, which causes the giant magnetostrictive rod to induce magnetostrictive deformation to move the output shaft and realize the output of the force and the displacement. Using the method of adjusting the driving DC signal can change the amplitude of the magnetic field. By this way, we can obtain the output of the displacement and the force in different sizes.

Finite Element Simulation of the Coupling Model of Giant Magnetostrictive Actuator (GMA)
In this study, the MATLAB2009a software demo (which is a commercial mathematics software produced by MathWorks, MA, USA) is used to program the model, and a demo version of COMSOL 4.0a software demo (which eveloped by the COMSOL group who founded in Stockholm, Sweden) [17] is used to simulate the coupling model of giant magnetostrictive actuator. The analysis process includes the following: Firstly, COMSOL finite element simulation software is used to construct the geometric model and mesh the actuator, and the coefficient matrix is iterated through the MATLAB software. Secondly, both two are invoked to complete the coupled analysis of the material model and the dynamic model. In theory, the magnetic field of the coil fills the whole space. In the finite element simulation, a larger closed cuboid air field is built to substitute the space. The giant magnetostrictive actuator is placed in the air field, and the boundary of air field is set to be magnetic insulation.
The diameter of the GMM rod is 10 mm, and its length is 450 mm. The material of the shell is chosen to be 45 steel. The diameter of single wire of driving coil is 1 mm, and its height is 22 mm. The total number turns of winding coil is 9570 turns. The solution process of the model is shown in Figure 3: Firstly, the current values are inputted and the model is meshed. Secondly, the increment of the initial magnetic field strength can be obtained by the input of initial coefficient matrix. Thirdly, the coefficient matrix is updated by solving the model, where the increments of magnetic induction and strain should be calculated so that further new incremental magnetic field intensities can be obtained. The cycle is repeated until the current waveform function is terminated.

Simulation Result of Magneto Mechanical Coupling Model of Giant Magnetostrictive Actuator (GMA)
A meshing strategy is chosen in the finite analysis as it is shown in the literature [19,21]. In order to reflect the output characteristics of giant magnetostrictive material rods in a better way, the giant magnetostrictive material rods and output shaft are meshed more intensively (Figure 4). Because the driving magnetic circuit and driving coil are of a symmetrical structure, and both of them have large sizes, the mesh is relatively sparse so that the number of elements and the number of degrees of freedom in the model solving process can be reduced. Then, the operation efficiency can be improved.

Simulation Result of Magneto Mechanical Coupling Model of Giant Magnetostrictive Actuator (GMA)
A meshing strategy is chosen in the finite analysis as it is shown in the literature [19,21]. In order to reflect the output characteristics of giant magnetostrictive material rods in a better way, the giant magnetostrictive material rods and output shaft are meshed more intensively (Figure 4). Because the driving magnetic circuit and driving coil are of a symmetrical structure, and both of them have large sizes, the mesh is relatively sparse so that the number of elements and the number of degrees of freedom in the model solving process can be reduced. Then, the operation efficiency can be improved.

Simulation Result of Magneto Mechanical Coupling Model of Giant Magnetostrictive Actuator (GMA)
A meshing strategy is chosen in the finite analysis as it is shown in the literature [19,21]. In order to reflect the output characteristics of giant magnetostrictive material rods in a better way, the giant magnetostrictive material rods and output shaft are meshed more intensively (Figure 4). Because the driving magnetic circuit and driving coil are of a symmetrical structure, and both of them have large sizes, the mesh is relatively sparse so that the number of elements and the number of degrees of freedom in the model solving process can be reduced. Then, the operation efficiency can be improved.  The maximum input electric current is set to be 5 A, and preloading is 10 MPa. The axial magnetic field distribution curve is shown in Figure 5, from which we can see that GMM rod axial magnetic field strength is around 250 kA/m. The maximum input electric current is set to be 5 A, and preloading is 10 MPa. The axial magnetic field distribution curve is shown in Figure 5, from which we can see that GMM rod axial magnetic field strength is around 250 kA/m.   When the input current increases monotonically from 0 A to 5 A, the axial output displacement curve of the giant magnetostrictive actuator is shown as Figure 7. When the input current is from 0 A to 5 A, then from 5 A to 0 A, the hysteretic loop of the giant magnetostrictive actuator is shown in Figure 8.   The maximum input electric current is set to be 5 A, and preloading is 10 MPa. The axial magnetic field distribution curve is shown in Figure 5, from which we can see that GMM rod axial magnetic field strength is around 250 kA/m.   When the input current increases monotonically from 0 A to 5 A, the axial output displacement curve of the giant magnetostrictive actuator is shown as Figure 7. When the input current is from 0 A to 5 A, then from 5 A to 0 A, the hysteretic loop of the giant magnetostrictive actuator is shown in Figure 8. When the input current increases monotonically from 0 A to 5 A, the axial output displacement curve of the giant magnetostrictive actuator is shown as Figure 7. When the input current is from 0 A to 5 A, then from 5 A to 0 A, the hysteretic loop of the giant magnetostrictive actuator is shown in Figure 8.  As shown in Figure 7, the maximum output displacement is 498 μm , and when the input current is changed from 0 A to 3 A, the linearity of giant magnetostrictive actuator is much better. When the input current is greater than 3 A, the output displacement presents worse linearity. From Figure 8, the worst hysteresis happens in situations when the current is around 1 A, and the hysteresis is 25.02%.

Testing and Experimental Research on Giant Magnetostrictive Actuator
Based on the finite element analysis of the magneto mechanical coupling model of the giant magnetostrictive actuator [27], a system platform for performance testing of the developed giant magnetostrictive actuator is set up, and the experiment research on the output displacement characteristics of giant magnetostrictive actuator is conducted. Through comparison of the experimental and simulation results, the output displacement and hysteresis of the two systems will be comprehensively analyzed.

Experimental Platform
The overall experimental system platform is shown in Figure 9, it consists of a giant  As shown in Figure 7, the maximum output displacement is 498 μm , and when the input current is changed from 0 A to 3 A, the linearity of giant magnetostrictive actuator is much better. When the input current is greater than 3 A, the output displacement presents worse linearity. From Figure 8, the worst hysteresis happens in situations when the current is around 1 A, and the hysteresis is 25.02%.

Testing and Experimental Research on Giant Magnetostrictive Actuator
Based on the finite element analysis of the magneto mechanical coupling model of the giant magnetostrictive actuator [27], a system platform for performance testing of the developed giant magnetostrictive actuator is set up, and the experiment research on the output displacement characteristics of giant magnetostrictive actuator is conducted. Through comparison of the experimental and simulation results, the output displacement and hysteresis of the two systems will be comprehensively analyzed.

Experimental Platform
The overall experimental system platform is shown in Figure 9, it consists of a giant As shown in Figure 7, the maximum output displacement is 498 µm, and when the input current is changed from 0 A to 3 A, the linearity of giant magnetostrictive actuator is much better. When the input current is greater than 3 A, the output displacement presents worse linearity. From Figure 8, the worst hysteresis happens in situations when the current is around 1 A, and the hysteresis is 25.02%.

Testing and Experimental Research on Giant Magnetostrictive Actuator
Based on the finite element analysis of the magneto mechanical coupling model of the giant magnetostrictive actuator [27], a system platform for performance testing of the developed giant magnetostrictive actuator is set up, and the experiment research on the output displacement characteristics of giant magnetostrictive actuator is conducted. Through comparison of the experimental and simulation results, the output displacement and hysteresis of the two systems will be comprehensively analyzed.

Experimental Platform
The overall experimental system platform is shown in Figure 9, it consists of a giant magnetostrictive actuator, vibration isolation platform, constant current source, control system, displacement detection system, and constant temperature water cooling system. magnetostrictive actuator, vibration isolation platform, constant current source, control system, displacement detection system, and constant temperature water cooling system. In this study, a prototype of a giant magnetostrictive actuator was designed and manufactured as shown in Figure 10. For the system platform, which is built for the prototype, industrial PC (PC is an abbreviation for personal computer) is taken as its control system. Firstly, the input voltage signal is entered from the industrial PC to controllable DC (DC is an abbreviation for direct-current) voltage constant-current power supply named DH1716A (Beijing Dahua Radio Instrument Factory, Beijing, China), the constant current source receives the signal and outputs the corresponding current. The driving coil generates a magnetic field under the effect of the input current. Because of the action of the magnetic field, GMM rod generates deformation and outputs displacement. The displacement signal is detected by the LVDT (LVDT is an abbreviation for linear variable differential transformer, which is a linear displacement sensor) micrometer (manufactured by Bosch Precision Measurement and Control Company, which was founded in Stuttgart, Germany, the version is MDS-L-0500-M6-1A), whose accuracy reached the micrometer-level, and which will finally be sent to industrial PC. The fluxmeter used in the experiment is the HT701 digital fluxmeter produced by Shanghai Hengtong Magnetism Technology Company, Ltd.,which was founded in Shanghai, China. Its measuring range can be divided into six steps (0.2, 0.5, 1, 2, 4, 8) × 20 wb, and the measurement accuracy can reach ±1% of the full scale. The HAP-100 series air-cushion vibration isolation platform manufactured by Shanghai Tiannuo Electromechanical Company, Ltd. (which was founded in Shanghai, China) was used to meet the vibration isolation requirements in the experimental process.
The application test software developed based on visual basic (VB) language runs directly on the industrial PC. It communicates with the constant current source and the LVDT micrometer through the data serial port line. By controlling the output current magnitude of the constant current source, the real-time displacement value displayed from the micrometer can be read, and the test In this study, a prototype of a giant magnetostrictive actuator was designed and manufactured as shown in Figure 10. For the system platform, which is built for the prototype, industrial PC (PC is an abbreviation for personal computer) is taken as its control system. Firstly, the input voltage signal is entered from the industrial PC to controllable DC (DC is an abbreviation for direct-current) voltage constant-current power supply named DH1716A (Beijing Dahua Radio Instrument Factory, Beijing, China), the constant current source receives the signal and outputs the corresponding current. The driving coil generates a magnetic field under the effect of the input current. Because of the action of the magnetic field, GMM rod generates deformation and outputs displacement. The displacement signal is detected by the LVDT (LVDT is an abbreviation for linear variable differential transformer, which is a linear displacement sensor) micrometer (manufactured by Bosch Precision Measurement and Control Company, which was founded in Stuttgart, Germany, the version is MDS-L-0500-M6-1A), whose accuracy reached the micrometer-level, and which will finally be sent to industrial PC. magnetostrictive actuator, vibration isolation platform, constant current source, control system, displacement detection system, and constant temperature water cooling system. In this study, a prototype of a giant magnetostrictive actuator was designed and manufactured as shown in Figure 10. For the system platform, which is built for the prototype, industrial PC (PC is an abbreviation for personal computer) is taken as its control system. Firstly, the input voltage signal is entered from the industrial PC to controllable DC (DC is an abbreviation for direct-current) voltage constant-current power supply named DH1716A (Beijing Dahua Radio Instrument Factory, Beijing, China), the constant current source receives the signal and outputs the corresponding current. The driving coil generates a magnetic field under the effect of the input current. Because of the action of the magnetic field, GMM rod generates deformation and outputs displacement. The displacement signal is detected by the LVDT (LVDT is an abbreviation for linear variable differential transformer, which is a linear displacement sensor) micrometer (manufactured by Bosch Precision Measurement and Control Company, which was founded in Stuttgart, Germany, the version is MDS-L-0500-M6-1A), whose accuracy reached the micrometer-level, and which will finally be sent to industrial PC. The fluxmeter used in the experiment is the HT701 digital fluxmeter produced by Shanghai Hengtong Magnetism Technology Company, Ltd.,which was founded in Shanghai, China. Its measuring range can be divided into six steps (0.2, 0.5, 1, 2, 4, 8) × 20 wb, and the measurement accuracy can reach ±1% of the full scale. The HAP-100 series air-cushion vibration isolation platform manufactured by Shanghai Tiannuo Electromechanical Company, Ltd. (which was founded in Shanghai, China) was used to meet the vibration isolation requirements in the experimental process.
The application test software developed based on visual basic (VB) language runs directly on the industrial PC. It communicates with the constant current source and the LVDT micrometer through the data serial port line. By controlling the output current magnitude of the constant current source, the real-time displacement value displayed from the micrometer can be read, and the test The fluxmeter used in the experiment is the HT701 digital fluxmeter produced by Shanghai Hengtong Magnetism Technology Company, Ltd.,which was founded in Shanghai, China. Its measuring range can be divided into six steps (0.2, 0.5, 1, 2, 4, 8) × 20 wb, and the measurement accuracy can reach ±1% of the full scale. The HAP-100 series air-cushion vibration isolation platform manufactured by Shanghai Tiannuo Electromechanical Company, Ltd. (which was founded in Shanghai, China) was used to meet the vibration isolation requirements in the experimental process.
The application test software developed based on visual basic (VB) language runs directly on the industrial PC. It communicates with the constant current source and the LVDT micrometer through the data serial port line. By controlling the output current magnitude of the constant current source, the real-time displacement value displayed from the micrometer can be read, and the test data will be output as TXT (TXT is an abbreviation for text document format) text document for archiving.
The constant temperature hydrocooling system consists of four parts: cooling water tank, temperature sensor, cooling fan, and heating rod. The temperature controller used in the experiment is an ATC-800 Microprocessor Temperature Controller made in Shanghai Jingchuang electric appliance manufacturing Company. Ltd., Shanghai, China, whose measuring accuracy is ±1 • C and whose indicating range is 0 • C-45 • C.

Experimental Research on Giant Magnetostrictive Actuator
In this study, a magnetostrictive experimental study of a giant magnetostrictive actuator is performed based on an experimental system of a giant magnetostrictive actuator. Through the acquired measurement data, and with the help of the analysis part of the test software we have already written, the useful information is extracted from the measured data and it will be analyzed.
(1) The experiment of output characteristic When the pre-pressure is 10 MPa and the input current is from 0 A to 5 A, then from 5 A to 0 A, the hysteretic loop of the giant magnetostrictive actuator is shown in Figure 11. When the input current increases from 0 A to 5 A monotonically, the relationship between output displacement and current is in the lift range. When the input current decreases from 5 A to 0 A monotonously, the relationship between the output displacement and the current is in the return stage, and the output displacement exhibits a significant hysteretic effect throughout the entire return stage. data will be output as TXT (TXT is an abbreviation for text document format) text document for archiving. The constant temperature hydrocooling system consists of four parts: cooling water tank, temperature sensor, cooling fan, and heating rod. The temperature controller used in the experiment is an ATC-800 Microprocessor Temperature Controller made in Shanghai Jingchuang electric appliance manufacturing Company. Ltd., Shanghai, China, whose measuring accuracy is 1 ± °C and whose indicating range is 0 °C-45 °C.

Experimental Research on Giant Magnetostrictive Actuator
In this study, a magnetostrictive experimental study of a giant magnetostrictive actuator is performed based on an experimental system of a giant magnetostrictive actuator. Through the acquired measurement data, and with the help of the analysis part of the test software we have already written, the useful information is extracted from the measured data and it will be analyzed.
(1) The experiment of output characteristic When the pre-pressure is 10 MPa and the input current is from 0 A to 5 A, then from 5 A to 0 A, the hysteretic loop of the giant magnetostrictive actuator is shown in Figure 11. When the input current increases from 0 to 5 A monotonically, the relationship between output displacement and current is in the lift range. When the input current decreases from 5 A to 0 monotonously, the relationship between the output displacement and the current is in the return stage, and the output displacement exhibits a significant hysteretic effect throughout the entire return stage.  The hysteresis is defined as follows: where y 1 is the displacement during the lift range, y 1i denotes the displacement during the lift stage (when the current is i), and y 2i represents the displacement during the return stage (when the current is i). From the equation above, the hysteresis is 29.74%.
(2) Comparison between experimental results and simulation results The experimental results are compared with the simulation results (as shown in Table 1). The results of finite element analysis and the results of output data are basically consistent, and most of the relative errors are stable at less than 15%, which means the established finite element model of magneto mechanical coupling is reasonable and the analysis result is also accurate. This model presents certain reference effects in the design process of giant magnetostrictive actuators. Figure 12 shows the comparison chart of the hysteretic loop between the experimental results and the finite element analysis results (the real line shows the experimental results, and the dotted line shows the finite element analysis results). (2) Comparison between experimental results and simulation results The experimental results are compared with the simulation results (as shown in Table 1). The results of finite element analysis and the results of output data are basically consistent, and most of the relative errors are stable at less than 15%, which means the established finite element model of magneto mechanical coupling is reasonable and the analysis result is also accurate. This model presents certain reference effects in the design process of giant magnetostrictive actuators. Figure 12 shows the comparison chart of the hysteretic loop between the experimental results and the finite element analysis results (the real line shows the experimental results, and the dotted line shows the finite element analysis results).   It can be seen from Figure 12 that the experimental results (the blue curve) are basically coincident with the hysteretic loops of the finite element analysis (the red dotted curve). When the current value is lower than 2.0 A, the deviation between the experimental curve and the simulation curve is large; when the current value exceeds 4 A, the deviation is also significantly larger, the current interval with the smallest deviation is between 2.0 A and 4.0 A. The reason for this phenomenon is that there is an optimal working current range for the giant magnetostrictive material, under the action of outside the optimal working current range, the hysteresis effect is obvious, and the displacement output error is large. In the actual application process, it should try to make it work in the best current range (the optimal working current range needs to be drawn by experiment). At the same time, there is a coincidence error between the finite element simulation curve and the experimental curve, which is because of the measurement error of the LVDT micrometer and the influence of the temperature change for the giant magnetostrictive material during the experiment. The giant magnetostrictive material has an optimum working temperature range, although the constant temperature treatment is adopted, the magnetic properties of the giant magnetostrictive material change when the temperature value changes because of the influence of heat generated by the current. At the beginning, the giant magnetostrictive material works at a constant temperature of 25 • C, when the current is entered for a period of time, the internal temperature of the rod of the giant magnetostrictive material increases, the hysteresis characteristic changes, and the linearization of the displacement output is improved. When the working time is slightly longer, the temperature rise is larger, the increase of the temperature affects the hysteresis effect, and curvature of the experimental curve of the displacement output is inconsistent with the simulation curve. The optimal working temperature and the optimum working current range require further post-experimental studies. At the same time, the multi-physics coupling characteristics of mechanical giant magnetostrictive materials consisting of the hysteresis characteristics, the mechanical properties, temperature, and others need to be further studied based on the research of this paper.
Comparing the experimental hysteresis with the finite element analysis hysteresis, the experimental result is 29.74%, the finite element analysis result is 25.02%, and the relative error is 15.87%. This shows that the proposed three-dimensional magneto mechanical coupling finite element model can describe the dynamic performance of the giant magnetostrictive actuator accurately and can play an auxiliary role in structural optimization.

Conclusions
For the giant magnetostrictive actuator, in this paper, at first, the eigen model of giant magnetostrictive material is constructed in combination with the free energy hysteresis model, and according to the fundamental equations of magnetic field and mechanical field of the giant magnetostrictive actuator, the three-dimensional nonlinear magnetic machine coupling model under weak solution is derived from MAXWELL equation and Newton's second law. Then, the axial magnetic field distribution of the giant magnetostrictive actuator is simulated using COMSOL finite element simulation software and MATLAB software, and the magnetic induction intensity value and output strain value with the current change controlled by the magnetic field are obtained, and the simulation results are verified experimentally. The results show that the experimental results of the output displacement and hysteresis of the giant magnetostrictive actuator are basically consistent with the simulation results, and the relative error is 15.87%, these show that the magnetic field strength on the central axis of the giant magnetostrictive actuator falls within the linear working area of the GMM under the action of the solenoid magnetic field (working current linear range is 2.0 A-4.0 A), and the magnetic circuit design is reasonable, the three-dimensional magnetic machine coupling model of the giant magnetostrictive actuator based on free energy hysteresis characteristics is designed to meet the input and output characteristics of the giant magnetostrictive actuator. The research conclusions have important guiding significance for the structural optimization of giant magnetostrictive actuator and the practical application of giant magnetostrictive materials.