Next Article in Journal
DM-MQTT: An Efficient MQTT Based on SDN Multicast for Massive IoT Communications
Next Article in Special Issue
Development of a Multiwavelength Visible-Range-Supported Opto–Ultrasound Instrument Using a Light-Emitting Diode and Ultrasound Transducer
Previous Article in Journal
Intelligent Ball Bearing Fault Diagnosis Using Fractional Lorenz Chaos Extension Detection
Previous Article in Special Issue
Equivalent Circuit Modeling for a Valveless Piezoelectric Pump
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Key Laboratory of Metallurgical Equipment and Control Technology of Ministry of Education, Wuhan University of Science and Technology, Wuhan 430081, China
2
Hubei Key Laboratory of Mechanical Transmission and Manufacturing Engineering, Wuhan University of Science and Technology, Wuhan 430081, China
3
Hubei Wuhan Tobacco Company Cigarette Logistics Distribution Center, Wuhan 430000, China
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(9), 3070; https://doi.org/10.3390/s18093070
Submission received: 29 July 2018 / Revised: 6 September 2018 / Accepted: 10 September 2018 / Published: 12 September 2018
(This article belongs to the Special Issue Recent Advances of Piezoelectric Transducers and Applications)

Abstract

:
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.

1. Introduction

At present, the functional materials used in the field of micro-actuation technology mainly include piezoelectric materials, shape memory alloys, and giant magnetostrictive materials [1,2,3,4,5,6].
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 method for disturbing force self-sensing in the application of giant magnetostrictive actuators, a Kelvin bridge combined with an active device is constructed instead of a conventional Wheatstone bridge for extraction of the induced voltage, the method for solving the nonlinear problems in GMA self-sensing signal extraction has been demonstrated. Mingzhang Luo [10] attempts to develop a portable, non-destructive evaluation method for assessing the length of installed rock bolts for inspection purposes, we proposed a portable device for the non-destructive evaluation of rock bolt conditions based on a giant magnetostrictive (GMS) actuator. The GMS actuator generates enough energy to ensure multiple reflections of the stress waves along the rock bolt and a lead zirconate titantate (PZT) sensor is used to detect the reflected waves. This paper can also be adopted for property assessment of others to accurately determine the rock bolt length. At present, GMA has reported on the application of precision fluid transmission and control, precision machining, precision peristaltic mechanism, sonar, vibration suppression, measurement and control, and so on.
Research [8,11,12,13] shows the following points: There is a nonlinear eigen relation between the output signal and the input signal of the GMA. It is a major problem in the practical application of intelligent materials for its intrinsic nonlinearity. Because of the strong coupling property between the material itself and the electromagnetic heat, it is very complicated for the nonlinear problem, which in turn makes modeling very difficult. In the research of a multi-physics dynamic magnetic-mechanical coupling model, scholars at home and abroad have also done many analyses and much research.
In the research by Yongxin Guo [14], a new Hammerstein model is proposed for modeling the rate-dependent hysteresis nonlinearity of a GMA, the proposed Hammerstein model can describe the hysteresis loops well within a frequency range of 0–100 Hz. However, based on this model, only a 2 DOF (DOF is an abbreviation for degree of freedom) control system is designed.
A novel adaptive filter is proposed to model the rate-dependent hysteresis nonlinearity in a giant magnetostrictive actuator [15]. In the proposed filter, generalized play operators are combined with linear delayed adaptive transversal filter to compose a new serial structure of adaptive filter model. The proposed adaptive filter is applied to model the rate dependent hysteresis of giant magnetostrictive actuator. The experimental results show that the proposed generalized play operator adaptive filter can describe the rate-dependent hysteresis behaviors including different single frequency input signal and multi-frequency composite input signal.
A.E. Clark [16] has proposed a linear piezomagnetic constitutive equation for a giant magnetostrictive actuator based on a large number of studies, this equation accurately describes the magnetic–mechanical coupling relationship of giant magnetostrictive materials, making it a magnetic–machine coupling basic equations in the study of giant magnetostrictive actuators.
As early as 1995, M.E.H. Benbouzid et al. [17,18] used the finite element simulation model to analyze the nonlinear dynamic characteristics of the giant magnetostrictive rod, which can guide the optimization of the magnetic circuit design of the giant magnetostrictive actuator. However, the model is a two-dimensional planar model, and its dynamic characteristics are obtained from the static characteristics of giant magnetostrictive materials. Because static modeling has many parameters, the calculation is complex and susceptible to other field-related parameters; therefore, the model itself has certain limitations.
Azoum [19] established a three-dimensional generalized finite element model based on the coupled constitutive equation of the giant magnetostrictive actuator, and simulated the solenoid coil. Simultaneously, Benatar [20] established a three-dimensional electro-magnetic-mechanical coupling model of GMM transducers used the multiphysics coupling software FEMLAB (the predecessor of COMSOL software which eveloped by the COMSOL group who founded in Stockholm, Sweden). Zhao Zhangrong and others from Zhejiang University [21] also used the weak solution of COMSOL software to simulate and analyze the electro–magnetic–mechanical coupling model of intelligent mast components of GMM, but the coefficient matrix substituted is a constant coefficient matrix, Therefore, the simulation result obtained is linear, and there is a certain error with the experimental results.
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.

2. The Theoretical Basis of the Magnetic Machine Coupling Model for the Giant Magnetostrictive Actuator (GMA)

2.1. 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:
× H = J S + J E
In this equation, H denotes the magnetic field intensity vector.
Electric field intensity E can be expressed by magnetic vector potential:
E = A t ϕ
where ϕ is a scalar potential function.
Then the applied current density caused by eddy current effect is expressed as (3):
J E = σ E = σ ( A t ϕ ) = σ ( A t 0 ) = σ A t
where σ is electric conductivity.
Then the magnetic flux density B is expressed by magnetic vector potential:
B = × A = | i j k x y z A x A y A z | = ( A z y A y z ) i + ( A x z A z x ) j + ( A y x A x y ) k

2.2. 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:
· T + f B = m d 2 u d t 2
where the displacement u = [ u x u y u z ] T , the stress T = [ T 1 T 12 T 13 T 21 T 2 T 23 T 31 T 32 T 3 ] , T 12 = T 21 , 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):
· T + f B = m d 2 u d t 2 + c u t
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:
S = s u
where s = [ x 0 0 0 y 0 0 0 z y 0 z x z 0 0 y x ] T . If it is substituted into Equation (7), Formula (8) is obtained as follows:
[ S 1 S 2 S 3 S 12 S 23 S 13 ] = s u = [ x 0 0 0 y 0 0 0 z y 0 z x z 0 0 y x ] T [ u x u y u z ]

2.3. The Conversion of the Weak Form Solution

(1)   Einstein notation
According to Einstein Notation [24], tensor i j k is a Levi–Civita signal, which is defined as follows:
i j k = { + 1 ,       i f ( i , j , k ) = ( 1 , 2 , 3 ) , ( 2 , 3 , 1 ) o r ( 3 , 1 , 2 ) ; 1 ,       i f ( i , j , k ) = ( 3 , 2 , 1 ) , ( 2 , 1 , 3 ) o r ( 1 , 3 , 2 ) ; 0 ,             i = j , j = k , o r , k = i ;
The cross product of the two variables can be expressed as follows:
( a × b ) i = i j k a j b k
The curl of the variables is shown as follows:
( × a ) i = i j k a k x j
Use Einstein Notation to mark Equations (1) and (3):
( × H ) i = ( J S ) i + ( J E ) i ( J E ) i = ( σ A t )
Substituting Equation (11) into the above two equations:
i j k H k x j = ( J S ) i + ( σ A t ) i = ( J S ) i σ A i t
Use Einstein Notation to mark Equation (6) and expand the formula:
( · T ) i + ( f B ) i = ( m d 2 u d t 2 ) i + ( c u t ) i
T i j x j + ( f B ) i = m d 2 u i d t 2 + c u i t
Use weighted residual method to integrate Equations (13) and (15):
V B i j k H k x j ψ i d V + V B σ A i t ψ i d V = V B ( J S ) i ψ i d V V u T i j x j ϕ i d V + V u ( f B ) i ϕ i d V = V u m d 2 u i d t 2 ϕ i d V + V u c u i t ϕ i d V
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:
V B i j k H k x j ψ i d V = V B i j k ( H k ψ i ) x j d V V B i j k H k ψ i x j d V V u T i j x j ϕ i d V = V u ( T i j ϕ i ) x j d V V u T i j ϕ i x j d V
Define the divergence V · F d V = S F · n d S , according to divergence theorem, the first term on the right of the Equation (17) can be expressed as follows:
V B i j k ( H k ψ i ) x j d V = V B i j k H k ψ i n j d V V u ( T i j ϕ i ) x j d V = V u T i j ϕ i n j d V
Substitute the above Equation into equation (17):
V B i j k H k x j ψ i d V = V B i j k H k ψ i n j d V V B i j k H k ψ i x j d V V u T i j x j ϕ i d V = V u T i j ϕ i n j d V V u T i j ϕ i x j d V
Then the Equation (19) can be rearranged as follows:
V B i j k H k x j ψ i d V = V B i j k i j k H j n k ψ i d V + V B i j k H i ψ k x j d V V u T i j x j ϕ i d V = V u T i j ϕ i n j d V V u T i j ϕ i x j d V
Substitute the above equation into Equation (16):
V B i j k H i ψ k x j d V + V B σ A i t ψ i d V = V B i j k H j n k ψ i d V + V B ( J S ) i ψ i d V V u m d 2 u i d t 2 ϕ i d V + V u c u i t ϕ i d V + V u T i j ϕ i x j d V = V u T i j ϕ i n j d V + V u ( f B ) i ϕ i d V
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:
ψ i = δ A i ϕ i = δ u
Equation (21) can be expressed as follows:
V B i j k H i δ A k x j d V + V B σ A i t δ A i d V = V B i j k H j n k δ A i d V + V B ( J S ) i δ A i d V V u m d 2 u i d t 2 δ u i d V + V u c u i t δ u i d V + V u T i j δ u i x j d V = V u T i j δ u i n j d V + V u ( f B ) i δ u i d V
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):
V B H · ( × δ A ) d V + V B σ A t δ A d V = V B J S δ A d V + V B ( H × n ) · δ A d V V u m d 2 u d t 2 δ u d V + V u c u t δ u d V + V u T · δ u d V = V u ( T · n ) δ u d V + V u f B δ u d V
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 = T n , 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:
V B H · δ B d V + V B σ A t δ A d V = V B J S δ A d V + V B H T · δ A d V V u m d 2 u d t 2 δ u d V + V u c u t δ u d V + V u T · δ S d V = V u t · δ u d V + V u f B δ u d V

2.4. 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 d x and d ξ through the Jacobian matrix 𝒥 e , where d x = 𝒥 e d ξ , d V = det ( 𝒥 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 e A and q e u :
A e = N A ( ξ ) q e A u e = N u ( ξ ) q e u
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:
δ A e = N A ( ξ ) δ q e A δ u e = N u ( ξ ) δ q e u
Because A is a three-dimensional vector, the interpolating matrix N A presents three lines, N n A nods in each element ( N n A 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 q A columns, where N q A = 3 N n A . The dimension of the vector q e A is N q A , 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 q u , which depended on the number of node N n u . From the above information, it is observed that the total degree of freedom of each element is N q = N q u + N q A .
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:
N 1 = ξ 1 N 2 = ξ 2 N 3 = ξ 3 N 4 = 1 ξ 1 ξ 2 ξ 3 ,
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 .
x 1 = N 1 x 1 , 1 + N 2 x 1 , 2 + N 3 x 1 , 3 + N 4 x 1 , 4
Because the four-node tetrahedral element in the finite elements is used, the dimension of vector q e A 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).
N ( ξ ) = [ N 1 0 0 0 N 1 0 0 0 N 1 N 2 0 0 0 N 2 0 0 0 N 2 N 3 0 0 0 N 3 0 0 0 N 3 N 4 0 0 0 N 4 0 0 0 N 4 ]
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.
B e = × A e = × ( N A q e A ) = C e q e A S e = u e = ( N u q e u ) = G e q e u
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 𝒥 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:
e = 1 N A ( Δ H · C e δ q e A J e d Δ + Δ σ e N A q e A t · N A δ q e A J e d Δ )                                                               = b = 1 N S A Δ H T · N A δ q b A J e , S d Δ S + e = 1 N A Δ J s , e N A δ q e A J e d Δ e = 1 N u ( Δ T · G e δ q e u J e d Δ + Δ ρ e N u 2 q e u t 2 · N u δ q e u J e d Δ + Δ c e N u q e u t · N u δ q e u J e d Δ )                                                               = b = 1 N S u Δ S t b · N u δ q b u J e , S d Δ S
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 S A and N S u 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).
Δ J e d Δ = 1 1 1 1 1 1 det ( x ξ ) d ξ 1 d ξ 2 d ξ 3 = V e
The integral Δ S J S , b d Δ S means boundary element integration in the global coordinate, which represents the element area.
Δ S J S , b d Δ S = 1 1 1 1 det [ x i ξ i x i ξ j x j ξ i x j ξ j ] d ξ i d ξ j = A b

2.5. Magneto Mechanical Coupling Nonlinear Model

Magneto mechanical coupling relationship [26] of the magnetostrictive actuator system can be described by the piezomagnetic equation:
B = μ T H + d T S = d T H + S H T
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 d T 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:
Δ H = a Δ B b Δ S Δ T = b T Δ B + c Δ S
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:
H e = a e C e q e A b e G e q e u T e = b e T C e q e A + c e G e q e u
The above equation can be substituted into Equation (32), and the symbols of the integral variables in the equation are defined as follows:
k e u = Δ G e T c e G e J e d Δ k e A = Δ C e T a e C e J e d Δ k e u , A = Δ C e T b e G e J e d Δ d e A = Δ ( N A ) T σ e N A J e d Δ d e u = Δ ( N u ) T c e N u J e d Δ m e = Δ ( N u ) T ρ e N u J e d Δ f b u = Δ S ( N u ) T t J e , S d Δ S f b A = Δ S ( N A ) T H T J e , S d Δ S f e J = Δ ( N A ) T J s , e J e , S d Δ
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:
e = 1 N A ( d e A q e A t + k e A q e A k e u , A q e u ) · δ q e A = e = 1 N A f e J · δ q e A + e = 1 N S A f b A · δ q b A e = 1 N u ( m e 2 q e u t 2 + d e u q e u t + k e u q e u ( k e u , A ) T q e A ) · δ q e u = b = 1 N S u f b u · δ q b u
Reorganizing the above equation with variational principle, it can be written into a matrix form:
[ 0 0 0 m e ] [ q ¨ e A q ¨ e u ] + [ d e A 0 0 d e u ] [ q ˙ e A q ˙ e u ] + [ k e A k e u , A ( k e u , A ) T k e u ] [ q e A q e u ] = [ f e J + f b A f b u ]
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 e u and the eddy current damping d e A (the eddy current is depended on the conductivity). In the stiffness matrix, k e A depends mainly on the magnetic permeability and it characterizes the capability of the magnetic excitation system. The coupling matrix k e u , A characterizes the capability of the magnetic field converted into mechanical energy by the mechanical traction vector f b u under the surface traction and capability of the mechanical energy converted into the magnetic energy at the magnetic field load vector f e J + f b A .
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:
a = B H ( H 0 , T 0 ) b = B T ( H 0 , T 0 ) b T = S H ( H 0 , T 0 ) c = S T ( H 0 , T 0 )
Coefficients in above formula can be written into matrix form shown as (42):
[ a b b T c ] = [ B H B T S H S T ]
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 .

3. 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.

3.1. 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.

3.2. 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.

3.3. 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.
Figure 6 presents GMM rod’s λ H curve provided by Gansu Star Company (pre-pressure is 10 MPa), as shown in the figure, the linear working area of GMM is located between the AB segments. According to the design requirements, the intensity of the magnetic field distributed on the GMM rod should be not less than 40 kA/m, taking into consideration of flux leakage and other factors, we know that the giant magnetostrictive actuator designed in this paper can fully meet the design requirements.
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%.

4. 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.

4.1. 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.
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 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.

4.2. 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.
The hysteresis is defined as follows:
e = max | y 2 i y 1 i | max ( y 1 ) × 100 %
where y 1 is the displacement during the lift range, y 1 i denotes the displacement during the lift stage (when the current is i ), and y 2 i 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).
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.

5. 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.

Author Contributions

Conceptualization, Z.Y. and T.W.; Methodology, Z.Y.; Validation, Z.Y. and M.Z.; Formal Analysis, M.Z.; Investigation, M.Z.; Resources, Z.Y.; Data Curation, M.Z.; Writing—Original Draft Preparation, Z.Y.; Writing—Review & Editing, Z.Y.; Supervision, T.W.; Project Administration, Z.Y.; Funding Acquisition, Z.Y.

Funding

The research is funded by Natural Science Fund of Hubei Province [Grant No. 2016CFB581] and Key Laboratory Open Fund of Ministry of Education of Metallurgical Equipment and Control of Wuhan University of Science and Technology [2013A07].

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Park, G.; Bement, M.T.; Hartman, D.A.; Smith, R.E.; Farrar, C.R. The use of active materials for machining processs: A review. Int. J. Mach. Tools Manuf. 2007, 47, 2189–2206. [Google Scholar] [CrossRef]
  2. Schellekens, P.; Rosielle, N.; Vermeulen, H.; Vermeulen, M.; Wetzels, S.; Pril, W. Design for precision: Current status and trends. Ann. CIRP 1997, 47, 557–586. [Google Scholar] [CrossRef]
  3. Liang, S.Y.; Hecker, R.L.; Landers, R.G. Machining process monitoring and control: The state of the art. ASME J. Manuf. Sci. Eng. 2004, 126, 297–310. [Google Scholar] [CrossRef]
  4. Srinivasan, S.; McFarland, M. Smart Structure: Analysis and Design; Cambridge University Press: Cambridge, UK, 2001. [Google Scholar]
  5. Chopram, I. Review of state of art of smart structures and integrated systems. AIAA J. 2002, 40, 2145–2187. [Google Scholar] [CrossRef]
  6. Inman, D.J. Smart materials in damage detection and prognosis. In Proceedings of the Fifth International Conference on Damage Assessment of Structures, Southampton, UK, 1–3 July 2003; pp. 3–16. [Google Scholar]
  7. Luo, M.; Li, W.; Wang, J.; Wang, N.; Chen, X.; Song, G. Development of a Novel GuidedWave Generation System Using a Giant Magnetostrictive Actuator for Nondestructive Evaluation. Sensors 2018, 18, 779. [Google Scholar] [CrossRef] [PubMed]
  8. Jia, Z.; Yang, X.; Guo, D.; Hou, L. Theories and Methods of Designing Microdisplacement Actuator Based on Giant Magnetostrictive Materials. Chin. J. Mech. Eng. 2001, 37, 46–49. [Google Scholar] [CrossRef]
  9. Yang, Y.; Wang, L.; Tan, J.; Zhao, B. Induced Voltage Linear Extraction Method Using an Active Kelvin Bridge for Disturbing Force Self-Sensing. Sensors 2016, 16, 739. [Google Scholar] [CrossRef] [PubMed]
  10. Luo, M.; Li, W.; Wang, B.; Fu, Q.; Song, G. Measurement of the Length of Installed Rock Bolt Based on StressWave Reflection by Using a Giant Magnetostrictive (GMS) Actuator and a PZT Sensor. Sensors 2017, 17, 444. [Google Scholar] [CrossRef] [PubMed]
  11. Jenner, A.G.; Smith, R.J.E.; Wilkinson, A.J. Actuation and transduction giant magnetostrictive alloys. Mechatronics 2000, 10, 457–466. [Google Scholar] [CrossRef]
  12. Lacheisserie, E.D.T.D. Magnetostriction Theory and Applications of Magnetoelasticity; CRC Press, Inc.: Boca Raton, FL, USA, 1993. [Google Scholar]
  13. Yang, Z.; He, Z.; Yang, F.; Rong, C.; Cui, X. Design and analysis of a voltage driving method for electro-hydraulic servo valve based on giant magnetostrictive actuator. Int. J. Appl. Electromagnet. Mech. 2018, 57, 439–456. [Google Scholar] [CrossRef]
  14. Guo, Y.; Mao, J.; Zhou, K. Rate-Dependent Modeling and H∞ Robust Control of GMA Based on Hammerstein Model With Preisach Operator. IEEE Trans. Control Syst. Technol. 2015, 23, 2432–2439. [Google Scholar] [CrossRef]
  15. Zhang, Z.; Ma, Y.; Guo, Y. A Novel Nonlinear Adaptive Filter for Modeling of Rate-Dependent Hysteresis in Giant Magnetostrictive Actuators. In Proceedings of the 2015 IEEE International Conference on Mechatronics and Automation, Beijing, China, 2–5 August 2015; pp. 670–675. [Google Scholar]
  16. Clark, A.E. Magnetostrictive Rare Earth-Fe2 Compounds; Wohlfarh, E.P., Ed.; North-Holland Publishing Company: New York, NY, USA, 1980. [Google Scholar]
  17. Benbouzid, M.E.H.; Reyne, G.; Meunier, G.; Kvarnsjo, L.; Engdahl, G. Dynamic modelling of giant magnetostriction in Terfenol-D rods by the finite element method. IEEE Trans. Magnet. 1995, 31, 1821–1824. [Google Scholar] [CrossRef]
  18. Pawel, I.; Krzysztof, K.; Lech, N.; Ƚukasz, K. FE transient analysis of the magnetostrictive actuator. Int. J. Appl. Electromagnet. Mech. 2016, 51, S81–S87. [Google Scholar]
  19. Azoum, K.; Besbes, M.; Bouillault, F. 3D FEM of magnetostriction phenomena using coupled constitutive laws. Int. J. Appl. Electromagnet. Mech. 2004, 19, 367–371. [Google Scholar]
  20. Benatar, J.G. Fem Implementations of Magnetostrictive-Based Applications. Master’s Thesis, University of Maryland, College Park, MD, USA, 2005. [Google Scholar]
  21. Zhao, Z.; Wu, Y.; Gu, X.; Xu, J.; Ge, R. Three-dimensional nonlinear dynamic finite element model for giant magnetostrictive actuators. J. Zhejiang Univ. 2008, 42, 203–208. [Google Scholar]
  22. Zhong, W. Ferromagnetism, 2nd ed.; Science Press: Beijing, China, 1992. [Google Scholar]
  23. Yuan, L.; Hu, Q. A Plane Wave Discontinuous Petrov-Galerkin Method for Helmholtz Equation and Time-harmonic Maxwell Equations with Complex Wave Numbers. J. Numer. Methods Comput. Appl. 2015, 36, 185–196. [Google Scholar]
  24. Hou, X.; Fan, Z.; Gu, Y.; Yang, H. A Wavelet Interpolation Galerkin Algorithm for Static Electromagnetic Field Analysis in Irregular Regions. Chin. J. Comput. Phys. 2005, 22, 539–548. [Google Scholar]
  25. Kannan, K.S.; Dasgupta, A. A nonlinear Galerkin finite-element theory for modeling magnetostrictive smart structure. Smart Mater. Struct. 1997, 6, 341–350. [Google Scholar] [CrossRef]
  26. Jia, Z.; Guo, D. Principle and Application of Micro-Displacement Actuator for Giant Magnetostrictive Materials; Science Press: Beijing, China, 2008. [Google Scholar]
  27. Claeyssen, F.; Lhermet, N.; Le Letty, R.; Bouchilloux, P. Actuators transducers and motors based on giant magnetostrictive materials. J. Alloys Compd. 1997, 258, 61–73. [Google Scholar] [CrossRef]
Figure 1. Four node tetrahedral element and its coordinate system.
Figure 1. Four node tetrahedral element and its coordinate system.
Sensors 18 03070 g001
Figure 2. Structural diagram of giant magnetostrictive actuator with non-uniform coils. 1. shell; 2. coil; 3. inlet and outlet; 4. output shaft sleeve; 5. preloading spring; 6. front-end cover; 7. output shaft; 8. fixed plate; 9. backend cover; 10. small end cover; 11. giant magnetostrictive material (GMM) rod; 12. coil bobbin outer tube.
Figure 2. Structural diagram of giant magnetostrictive actuator with non-uniform coils. 1. shell; 2. coil; 3. inlet and outlet; 4. output shaft sleeve; 5. preloading spring; 6. front-end cover; 7. output shaft; 8. fixed plate; 9. backend cover; 10. small end cover; 11. giant magnetostrictive material (GMM) rod; 12. coil bobbin outer tube.
Sensors 18 03070 g002
Figure 3. Flowchart of the model solving.
Figure 3. Flowchart of the model solving.
Sensors 18 03070 g003
Figure 4. Meshing of a giant magnetostrictive actuator.
Figure 4. Meshing of a giant magnetostrictive actuator.
Sensors 18 03070 g004
Figure 5. The distribution of magnetic field in the center axis of GMM rod.
Figure 5. The distribution of magnetic field in the center axis of GMM rod.
Sensors 18 03070 g005
Figure 6. The λ H curve of the GMM rod.
Figure 6. The λ H curve of the GMM rod.
Sensors 18 03070 g006
Figure 7. Axial output displacement curve of giant magnetostrictive actuator.
Figure 7. Axial output displacement curve of giant magnetostrictive actuator.
Sensors 18 03070 g007
Figure 8. Hysteretic loop of a giant magnetostrictive actuator.
Figure 8. Hysteretic loop of a giant magnetostrictive actuator.
Sensors 18 03070 g008
Figure 9. Overall experimental system platform. GMA—giant magnetostrictive actuator.
Figure 9. Overall experimental system platform. GMA—giant magnetostrictive actuator.
Sensors 18 03070 g009
Figure 10. Experimental prototype of a giant magnetostrictive actuator.
Figure 10. Experimental prototype of a giant magnetostrictive actuator.
Sensors 18 03070 g010
Figure 11. Hysteresis loop under the action of current 5 A and preloading 10 MPa.
Figure 11. Hysteresis loop under the action of current 5 A and preloading 10 MPa.
Sensors 18 03070 g011
Figure 12. Contrast diagram of hysteretic loop between experiment and finite element analysis.
Figure 12. Contrast diagram of hysteretic loop between experiment and finite element analysis.
Sensors 18 03070 g012
Table 1. Comparison between the finite element model analysis and experimental results.
Table 1. Comparison between the finite element model analysis and experimental results.
Electric Current (A)Output Displacement (μm)Error (%)
Finite Element Simulation ResultsExperimental Results
0000
0.527.9414.6147.71
192.6567.5227.12
1.5205.88182.6511.28
2282.35293.993.96
2.5348.53358.522.79
3398.53397.40.28
3.5432.35430.270.48
4466.18455.422.31
4.5493.12475.993.47
5498.09491.631.30

Share and Cite

MDPI and ACS Style

Yu, Z.; Wang, T.; Zhou, M. Study on the Magnetic-machine Coupling Characteristics of Giant Magnetostrictive Actuator Based on the Free Energy Hysteresis Characteristics. Sensors 2018, 18, 3070. https://doi.org/10.3390/s18093070

AMA Style

Yu Z, Wang T, Zhou M. Study on the Magnetic-machine Coupling Characteristics of Giant Magnetostrictive Actuator Based on the Free Energy Hysteresis Characteristics. Sensors. 2018; 18(9):3070. https://doi.org/10.3390/s18093070

Chicago/Turabian Style

Yu, Zhen, Tao Wang, and Min Zhou. 2018. "Study on the Magnetic-machine Coupling Characteristics of Giant Magnetostrictive Actuator Based on the Free Energy Hysteresis Characteristics" Sensors 18, no. 9: 3070. https://doi.org/10.3390/s18093070

APA Style

Yu, Z., Wang, T., & Zhou, M. (2018). Study on the Magnetic-machine Coupling Characteristics of Giant Magnetostrictive Actuator Based on the Free Energy Hysteresis Characteristics. Sensors, 18(9), 3070. https://doi.org/10.3390/s18093070

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop