Next Article in Journal
Investigations on the Processing of Ceramic Filled Inks for 3D InkJet Printing
Next Article in Special Issue
Declining Performance of Silicone-Based Magnetorheological Elastomers after Accelerated Weathering
Previous Article in Journal
Sensor-Assisted Assessment of the Tribological Behavior Patterns of AA7075 Hybrid MMC Reinforced with Multi-Wall Carbon Nanotubes and Pulverized Fuel Ash
Previous Article in Special Issue
Multidisciplinary Design Optimization of a Novel Sandwich Beam-Based Adaptive Tuned Vibration Absorber Featuring Magnetorheological Elastomer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Asymmetric Hysteresis Modeling Approach Featuring “Inertial System + Shape Function” for Magnetostrictive Actuators

1
Laboratory for Adaptive Structures and Intelligent Systems (LASIS), Department of Vehicle Engineering, Hefei University of Technology, Hefei 230009, China
2
College of Mechanical Engineering, Anhui Science and Technology University, Fengyang 233100, China
*
Author to whom correspondence should be addressed.
Materials 2020, 13(11), 2585; https://doi.org/10.3390/ma13112585
Submission received: 14 April 2020 / Revised: 5 May 2020 / Accepted: 7 May 2020 / Published: 5 June 2020

Abstract

:
Hysteresis of the actuators based on magnetostrictive materials influences the control performance of the application systems. It is of importance and significance to establish an effective hysteresis model for the magnetostrictive actuators for precision engineering. In this paper, based on the analysis of the Duhem model, a first-order inertial system with hysteresis characteristic under harmonic input is used to describe the hysteresis caused by the inertia of the magnetic domains of magnetostrictive materials. Shape function is employed to describe the pinning of domain walls, the interactions of different magnetic domains of magnetostrictive materials, and the saturation properties of the hysteresis. Specifically, under an architecture of “inertial system + shape function” (ISSF-Duhem model), firstly a new hysteresis model is proposed for magnetostrictive actuators. The formulation of the inertial system is constructed based on its general expression, which is capable of describing the hysteresis characteristics of magnetostrictive actuators. Then, the developed models with a Grompertz function-based shape function, a modified hyperbolic tangent function-based shape function employing an exponential function as an offset function, a one-sided dead-zone operator-based shape function are compared with each other, and further compared with the classic modified Prandtl–Ishlinskii model with a one-sided dead-zone operator. Sequentially, feasibility and capability of the proposed hysteresis model are verified and evaluated by describing and predicting the hysteresis characteristics of a commercial magnetostrictive actuator.

1. Introduction

Magnetostrictive materials, a smart material, produce mechanical stress and strain when exposed to external magnetic fields [1]. Because of high energy density, large output strain/force, and ultrafast response (within several microseconds), magnetostrictive materials have been widely developed/studied in various applications, such as active vibration control [2] and ultrasonic positioning [3]. However, similar to other smart material-based actuators [4,5,6,7,8], such as piezoelectric ceramic-based actuators and shape memory alloy-based actuators, the magnetostrictive actuators behave strong hysteresis nonlinearity [9,10]. Such hysteresis nonlinearity would decrease the actuating accuracy and degrade control system performance.
Physical mechanism of the magnetostriction is shown in Figure 1 (where H , H s , λ , and λ s are the magnetic field, saturation magnetic field, magnetostriction strain, and saturation magnetostriction strain, respectively). The magnetic domains of magnetostrictive materials generate motion when exposed to external magnetic fields, which involves the motion of the domain walls and the rotation of the magnetic domains [11]. Firstly, impedances to domain wall motion due to pinning sites encountered by the domain walls as they move are the reasons for the magnetostrictive hysteresis. Secondly, the magnetic domain itself has inertia. During the motion, the magnetic domain is inevitably affected by its own inertia. Inertia would be one reason for the magnetostrictive hysteresis. Because of the non-uniformity of magnetostrictive materials, the effect of inertia on different magnetic domains is also different with each other. Thirdly, during the motion, interactions between different magnetic domains of magnetostrictive materials occur, which also appear as hysteresis macroscopically. It can also be seen from Figure 1c that the magnetostrictive material has a maximal magnetostrictive strain. This can be used to explain the saturation characteristics of magnetostrictive hysteresis.
Investigations on the nonlinearity of magnetostrictive actuators and hysteresis models have been conducted for years. Jiles–Atherton (J–A) model [11] and its extended models [12,13,14,15,16,17] were proposed to capture the hysteresis nonlinearity via the domain wall theory. However, the J-A model shows disadvantages in describing minor loops [16]. The extended J–A models are actually mainly aiming at solving such issue. In reference [18], a modified Jiles–Atherton–Sablik model is constructed for asymmetry in magnetomechanical effect under tensile and compressive stress. Operator-based hysteresis models in the phenomenological hysteresis models describe the nonlinearity by weighted superposition of elementary hysteresis operators with simple mathematical forms. The typical operator-based hysteresis models are the relay operator of the Preisach model [10,19,20,21,22], the KP kernel of the Krasnoselskii–Pokrovskii model [7,23], and the play/stop operator of the Prandtl–Ishlinskii model [24,25,26,27,28]. The disadvantages of the operator-based hysteresis models are the characteristics of the model, including symmetry and saturation, depend on the elementary hysteresis operator. The more complex the shape is, the more the number of operators will be used. In reference [29], the modified Prandtl–Ishlinskii model is constructed with a one-sided dead-zone operator (PD-Model) and used to model complex hysteresis. Another kind of phenomenological hysteresis model—differential equation-based hysteresis models—need only one auxiliary nonlinear differential equation to describe the hysteretic behavior, which shows the advantage of higher efficiency [30]. In reference [31], the Bouc–Wen model was used to describe the hysteresis nonlinearity of magnetostrictive actuators. In references [32,33], a new Duhem model is established by incorporating a shape function to describe the hysteresis of magnetorheological fluid actuators. The Duhem model describing the hysteresis characteristics of ferromagnetically soft materials was given in reference [34] but has not been found for the hysteresis characteristics of magnetostrictive materials.
The Duhem model uses a differential equation to characterize the output change with the input direction [35]. The specific input–output relationship is
d w d t = f 1 u t , w t u ˙ + t f 2 u t , w t u ˙ t
u ˙ ± t = u ˙ t ± u ˙ t 2
where f 1 u t , w t and f 2 u t , w t C 0 R 2 ; w t and u t , the output and input functions, are continuous differentiable functions, respectively. When the input function u t increases, the output function w t increases along a path, and when u t decreases, w t decreases along another path.
Construction of the functions f 1 u t , w t and f 2 u t , w t is important for describing the system hysteresis when employing the Duhem model. The Duhem model has defined the properties of the functions f 1 u t , w t and f 2 u t , w t , but it is still challenging to construct a suitable function based on hysteresis characteristics in practice.
Based on the analysis and literature reviewed above, especially the analysis of the physical mechanism of the hysteresis of magnetostrictive actuators, the technical contribution of this paper is that a hysteresis modeling methodology based on a first-order inertial system and shape function (ISSF-Duhem model) is constructed for magnetostrictive actuators with properties of asymmetry and saturation hysteresis. Specifically, a differential equation of the inertial system is used to represent the hysteresis caused by the inertia in the motion of magnetic domains, and a general expression of the inertial system describing the hysteresis is established. As compared with the original Duhem model, the constructed expression reduces the number of unknown functions. For the hysteresis caused by the pinning of domain walls, the interactions between the magnetic domains and the saturation characteristics of the hysteresis, a shape function is incorporated to correct the hysteresis curve represented by the inertial system. The hysteresis model will be verified and evaluated based on the description and prediction of the hysteresis of a commercial magnetostrictive actuator through experimental tests. In this study, the obtained simulation results using an ISSF-Duhem model with different shape functions (Grompertz function-based shape function, modified hyperbolic tangent function-based shape function, and one-sided dead-zone operator-based shape function) will be compared with the results obtained by using a PD-model.

2. ISSF-Duhem Model

According to Equation (1), the Duhem model essentially represents the hysteresis as the input-output relationship of a first-order system. Based on this, this study starts from the analysis of the hysteresis phenomenon of the inertial system and then constructs a new hysteresis model.

2.1. Hysteresis in an Inertial System

The first-order inertial system can be expressed as
w ˙ T + w = u
where u , w , and T are the input, output, and the system parameter, respectively. According to Equation (3), Figure 2 presents the hysteresis performance of the inertial system. As seen from the figure, the hysteresis is related to the system parameter T .

2.2. Construction of Inertial System

As shown in Figure 1, the magnetic domains of magnetostrictive materials generate motion when exposed to external magnetic fields. The effect of inertia on different magnetic domains is also different with each other because of the non-uniformity of magnetostrictive materials. If the hysteresis system is considered to be a first-order inertial system, then the parameter T of the system is variable.
Based on reference [36], let
T = α g u ˙
where u ˙ is the input rate of the system; g u ˙ represents the function of the input rate, which is positively homogeneous; α is a constant related to the characteristics of the system.
Substituting Equation (4) into Equation (3), and to adapt to different hysteresis phenomena, expression is given as
w ˙ + α g u ˙ w = α γ g u ˙ f u
where f 1 u t , w t is a function of u and is related to the characteristics of the hysteresis itself; γ is a constant.
Let g u ˙ = u ˙ , then Equation (5) can be rewritten as
w ˙ + α u ˙ w = α γ u ˙ f u
Equation (6) represents the inertial system hysteresis model, and its input–output relationship can be further expressed as
w ˙ = α u ˙ β α f u w
where β = α γ .
Equation (7) can be further rewritten as
d w U d u = β f u α w U                       u ˙ > 0 d w D d u = β f u + α w D                   u ˙ < 0  
where   w U represents the value on the rising curve of the hysteresis loop; w D represents the value on the falling curve of the hysteresis loop. Assuming that the two endpoints of the hysteresis loop are ( u L , w L ) and ( u R , w R ), the function expression of the hysteresis loop rise and fall curves obtained by Equation (8) are
w U u = β α f u + w L β α f u L e α u u L β α e α u u L u f ε e α ε d ε u u L ,   u R
w D u = β α f u w R β α f u R e α u R u β α e α u u u R f ε e α ε d ε u u L ,   u R
In Equation (9), f ε e α ε is bounded within ( u L ,   u R ), then
w L β α f u L e α u u L β α e α u u L u f ε e α ε d ε w L β α f u L + β α f τ e α τ u R u L = ρ L
where τ u L , u R and is a constant; ρ L is the maximal value of the absolute value of the hysteretic part on the rising curve of the hysteresis loop.
Similarly, it can be obtained by Equation (10),
w R β α f u R e α u R u β α e α u u u R f ε e α ε d ε w R β α f u R + β α f ε e α φ u R u L = ρ R  
where φ u L ,   u R and is a constant; ρ R is the maximal value of the absolute value of the hysteretic part on the falling curve of the hysteresis loop.
Rewrite Equations (11) and (12) as
w u = β α f u + σ u
σ u = w L β α f u L e α u β α e α u u L u f ε e α ε d ε w R β α f u R e α u R u β α e α u u u R f ε e α ε d ε
where
σ u ρ = max ρ L , ρ R
In Equation (7), f u is a monotone increasing, piecewise smooth, odd function with a bounded derivative, that is, lim u f u is bounded. According to these conditions, the following function is selected
f u = u
Then, Equation (7) can be rewritten as
w ˙ = α u ˙ β α u w

2.3. Construction of Shape Function

The hysteresis curve of the inertial system represented by Equation (17) exhibits centrosymmetry. However, the hysteresis curve of the magnetostrictive actuators has the characteristics of asymmetry and saturation, due to the pinning of domain walls, the interactions between different magnetic domains of the magnetostrictive materials during motion, and the saturation characteristics of the magnetostrictive strain. A shape function is necessary to be incorporated to correct the hysteresis curve represented by the inertial system. In general, sigmoid functions and one-sided dead-zone operator are often used to describe the hysteresis of the magnetostrictive actuators. In this study, we construct shape functions based on these functions. It is noted that the Grompertz function, the hyperbolic tangent function, and the one-sided dead-zone operator are chosen as the basic functions for constructing the shape functions as examples; other similar functions can also be used.

2.3.1. Grompertz Function-Based Shape Function

Grompertz function is employed as the shape function in this study. It is expressed as
z 1 x = 1 + 2 e e μ 1 x μ 2
where μ 1 and μ 2 are both constants determining the shape of the Grompertz function output. The function z 1 x is within the range of [−1, 1].
When x = 0 ,
z 1 0 = 1 + 2 e e μ 2
It is seen that z 1 0 is dependent on the value of μ 2 . In order to describe the hysteresis of magnetostrictive actuators, Equation (18) can be given by
z 1 x = z 1 x z 1 0 = 2 e e μ 1 x μ 2 2 e e μ 2
Combining Equations (17) and (20) and using w in Equation (17) as the input of the shape function, considering that the hysteresis of magnetostrictive actuators generally contains a linear component, the ISSF-Duhem model with Grompertz function-based shape function (G-Model) describing the hysteresis of magnetostrictive actuators can be given by
w ˙ = α u ˙ β α u w   z 1 w = 2 e e μ 1 w μ 2 2 e e μ 2 y = k 1 u + k 2 z 1
where y represents the displacement of magnetostrictive actuators; k 1 and k 2 are both constants.
The schematic of the proposed ISSF-Duhem hysteresis model approach is presented in Figure 3. As compared with the Duhem model represented by Equation (1), the ISSF-Duhem model reduces the number of unknown functions. Only a suitable shape function is necessary to be selected according to the shape characteristics of the hysteresis curve to completely determine the model. As compared with continuous-time dynamic model of a class of backlash-like hysteresis in reference [37], formally Equation (7) of the ISSF-Duhem model is a specific case of Equation (4) of the mentioned research. However, Equation (7) has a clearer physical meaning since it is derived based on the analysis of the hysteresis of the magnetostrictive actuator, combined with characteristics of the first-order inertial system. At the same time, the hysteresis curve represented by Equation (7) is corrected by the shape function, making the ISSF-Duhem model more suitable for describing the hysteresis characteristics of magnetostrictive actuators.

2.3.2. Modified Hyperbolic Tangent Function-Based Shape Function

To compare, a shape function based on the hyperbolic tangent function [38] is constructed as
z 2 x = tanh δ x
where δ is a constant determining the slope of the hyperbolic tangent function.
The function z 2 x within the range of [−1, 1] cannot be directly used as a shape function because of its centrosymmetry. A modified shape function is constructed based on the hyperbolic tangent function:
z 2 x = 1 + 2 tanh e μ 3 x μ 4
where μ 3 and μ 4 are both constants.
The modified hyperbolic tangent function is not centrosymmetric, which is actually constructed through offsetting the hyperbolic tangent function by exponential function e x . In addition, the range of modified hyperbolic tangent function is still [−1, 1] by the coefficients.
Similar to the Grompertz function, the modified hyperbolic tangent function z 2 x depends on μ 4 , when x = 0 . To describe the hysteresis of magnetostrictive actuators, Equation (23) can be rewritten as
z 2 x = 2 tanh e μ 3 x μ 4 2 tanh e μ 4
Using Equation (24) with the input of w in Equation (17) instead of z 1 w in Equation (21), the ISSF-Duhem model with modified hyperbolic tangent function-based shape function (M-Model) can be given by
w ˙ = α u ˙ β α u w   z 2 w = 2 tanh e μ 3 w μ 4 2 tanh e μ 4 y = k 1 u + k 2 z 2

2.3.3. One-Sided Dead-Zone Operator-Based Shape Function

To further compare, the one-sided dead-zone operator is employed as the shape function. The one-sided dead-zone operator [29] is expressed as
Γ r i x = max x r i ,   0               for     r i > 0                   x                                               for     r i = 0 min x r i ,   0               for     r i < 0
where r i represents the threshold of the one-sided dead-zone operator; Γ r i represents the output of the one-sided dead-zone operator with the threshold r i .
Using Equation (26) with the input of w in Equation (17) instead of z 1 w in Equation (21), the ISSF-Duhem model with a one-sided dead-zone operator-based shape function (D-Model) can be given by
w ˙ = α u ˙ β α u w   z 3 w = i = n n g i Γ r i w y = k 1 u + k 2 z 3
where n represents the number of one-sided dead-zone operators; g i represents the weighting factor of the i-th one-sided dead-zone operator.
When r i 0 , the thresholds and weighting factors of one-sided dead-zone operators are selected as
r i = μ 5 i i = 0 , 1 , , n
g i = μ 6 e μ 7 i
When r i < 0 , the thresholds and weighting factors of one-sided dead-zone operators are selected as
r i = μ 8 i     i = n , , n + 1 , 1
g i = μ 9 e μ 10 i
where μ 5 , μ 6 , μ 7 , μ 8 , μ 9 , and μ 10 are all constants.

3. Experimental Validation

Aiming at verifying the performance of the ISSF-Duhem model for describing and predicting the hysteresis of magnetostrictive actuators, (i) the input–output characteristics of a commercial magnetostrictive actuator under different excitations are tested; and (ii) the ISSF-Duhem model is verified and evaluated with the measured experimental data.

3.1. Experimental Setup and Hysteresis Loop Test

Figure 4 shows the experimental setup for the hysteresis of the magnetostrictive actuator. The core component of the magnetostrictive actuator is a Terfenol-D rod with a length of 79.6 mm, a diameter of 20 mm, and the cut-off frequency of 76.07–380.34 Hz, which is prestressed about 3.023 MPa by a preload spring, and provided a bias magnetic field of about 3.2 kA/m by permanent magnets. The displacement of the magnetostrictive actuator is measured by an LVDT (Linear Variable Differential Transformer, Type: GA-2), which is a rebound type with a relatively small-stiffness spring. The excitation signal generated by LabVIEW software (version: 2014, National Instruments, Austin, TX, USA) and NI CompactRIO-9030, and amplified by a controllable current driver, is applied to the actuator. The input current and output displacement signals of the magnetostrictive actuator are acquired via the data acquisition system (Type: DEWE-501). The resonance of the mechanical setup is about 159.12 Hz.
Figure 5 presents the experimental tests in profile of the displacement versus current, when the magnetostrictive actuator is under harmonic input with different currents (1–6 A with an interval of 1 A). It can be seen that both the hysteresis and its asymmetry of the magnetostrictive actuator increases with the increase of excitation current.
Figure 6a,b present the triangle excitation current with various amplitudes of 6–1 A (with a decrement of 1 A) and output displacement of the magnetostrictive actuator under the excitation current, respectively. Observing Figure 6b, both the hysteresis and its asymmetry of the magnetostrictive actuator decrease with the decrease of the current.

3.2. Parameter Identification

Using the experimental data of the magnetostrictive actuator under harmonic input u(t) = Amsin(2πt), Am = 2, 3, 4 and 5 A, the multi-islands genetic algorithm (GA) is used to fit the established ISSF-Duhem model parameters [39]. Then, the performance of the ISSF-Duhem model is evaluated by comparing the model-based fitted results with the experimental results.
Because of the complexity and nonlinearity of the models as represented by Equations (21), (25), and (27), and the hysteresis itself as a relatively complex nonlinear phenomenon, parameter identification becomes a key issue in this study. There are many existing methods for parameter identification, such as a least square method, modified Gauss–Newton method, differential evolution method, particle swarm algorithm, and GA. GA originates from computer simulation study on the biological evolution systems and follows the principle of natural selection and survival of the fittest. The multi-islands GA in ISIGHT software is employed for parameter identification. The specific process of parameter identification given in Figure 7 consists of two parts. One is the fitting of experimental data and the other is the update of the model parameters. The fitting process is implemented in MATLAB software, and then the fitted result is transmitted to ISIGHT software (version: 5.9.2, Dassault Systèmes Simulia Corp., Providence, RI, USA). The multi-islands GA in the ISIGHT software adjusts the values of the parameters according to the received fitted results and transmits them to the MATLAB software (version: R2014b, the Mathworks, Inc., Natick, MA, USA). Finally, the parameter identification is realized by the joint calculation of the two software programs. The results of the parameter identification are listed in Table 1. The accuracy of the hysteresis of the magnetostrictive actuator described by the ISSF-Duhem model is evaluated with the mean square error (MSE), the maximal error (Emax), and relative error (Erelative). MSE, Emax, and Erelative are given by
MSE = 1 m i = 1 m ( D i s p sim , i D i s p exp , i ) 2
E max = max D i s p sim ,   i D i s p exp ,   i     ( i = 1 , , m )
E relative = E max max D i s p exp , i min D i s p exp , i
where m represents the number of data; D i s p sim , i and D i s p exp ,   i represent i-th fitted and experimental displacement value, respectively; E max is the maximal error; D i s p sim ,   i D i s p exp ,   i represents the absolute value.
Figure 8 presents the comparison between the fitted results of the ISSF-Duhem model with different shape functions as well as the PD-Model and the experimental tests in the profile of the displacement versus current, when the magnetostrictive actuator is under excitation u(t) = Am sin(2πt) with Am = 2, 3, 4, and 5 A. Observing Figure 8a, the ISSF-Duhem model can effectively describe the relationship between the displacement and current of the magnetostrictive actuator when the excitation current is 2 A. However, the fitted results of the ISSF-Duhem model with different shape functions are different. Comparing with the G-Model and M-Model, the fitted results of the D-Model are better for approximating the experimental results. At the same time, it can be seen that the fitted results of the D-Model and M-Model are better than those of the PD-Model for approximating the experimental results. Figure 8b (when Am = 3), 8c (when Am = 4) and 8d (when Am = 5) show the similar results.
Figure 9 presents the comparison between the fitted results of the ISSF-Duhem model with different shape functions as well as the PD-Model and the experimental tests in profile of the error versus time, when the magnetostrictive actuator is under excitation u(t) = Am sin(2πt) with Am = 2, 3, 4 and 5 A. As shown in Figure 9, all the G-Model, the M-Model, the D-Model, and the PD-Model can effectively track the displacement time history of the magnetostrictive actuator and effectively describe the hysteresis with less than 3 μm error. In other words, the ISSF-Duhem model can describe the hysteresis of magnetostrictive actuators as well as the PD-Model.
In order to analyze the effectiveness of the ISSF-Duhem model more in depth, the fitted results of the ISSF-Duhem model and the PD-Model are further analyzed with respect to experimental results. MSE, Emax, and Erelative between the fitted and experimental displacement under different excitations are calculated respectively, as shown in Figure 10. MSE reflects the overall condition of the fitted results of the ISSF-Duhem model and the PD-Model with respect to the experimental results. As seen from Figure 10a, MSE gradually increases with the increase of excitation current. When the excitation current is 5 A, MSE with a relatively small value is the largest, indicating that the ISSF-Duhem model is quite fair for the overall condition of the simulation of the hysteresis of the magnetostrictive actuator as well as the PD-Model. As shown in Figure 10a, MSE of the D-Model is smaller than that of the G-Model and M-Model. It can be concluded that, as compared with the Grompertz function-based and the modified hyperbolic tangent function-based shape function, the one-sided dead-zone operator-based shape function can enhance the performance of the ISSF-Duhem model. As seen from Figure 10a, the ISSF-Duhem model with the modified hyperbolic tangent function-based and the one-sided dead-zone operator-based shape function show better performance than the PD-Model. Emax and Erelative reflect the local condition of the fitted results of the ISSF-Duhem model and the PD-Model with respect to experimental results. Observing Figure 10b, with the increase of excitation current, Emax gradually increases, but Erelative decreases for the G-Model and the M-Model and increases for the D-Model slightly. In Figure 10b, Emax and Erelative of the D-Model are smaller than those of the G-Model, the M-Model, and the PD-model when the excitation current is less than 5 A. It indicates that the one-sided dead-zone operator-based shape function would enhance the performance of the ISSF-Duhem model. When the excitation current is 5A, Emax and Erelative of the M-Model is the smallest among the four models. Observing Figure 10a,b, the ISSF-Duhem model with the modified hyperbolic tangent function-based and the one-sided dead-zone operator-based shape function provide superior performance over the PD-Model.

3.3. Prediction Performance

Based on the identified parameters of the ISSF-Duhem model and the PD-Model in Section 3.2, the ISSF-Duhem model and the PD-Model are used to predict the hysteresis of the magnetostrictive actuator under two excitations: (i) u(t) = Amsin(2πt), Am = 1 and 6 A; (ii) triangle excitation current with various amplitudes of 6–1 A (with a decrement of 1 A). Then, the performance of the ISSF-Duhem model is further evaluated by comparing the model-based predicted results with those of the PD-Model and with the experimental results.
Figure 11a–f present the comparison between the predicted results of the ISSF-Duhem model with different shape functions as well as the PD-Model with the experimental tests in profiles of the displacement versus current, the displacement versus time and the error versus time, respectively, when the magnetostrictive actuator is under excitation u(t) = Amsin(2πt) with Am = 1 A and 6 A. As shown in Figure 11a–e, the ISSF-Duhem model can effectively predict the hysteresis of the magnetostrictive actuator when the excitation currents are 1 A and 6 A as well as the PD-Model. As shown in Figure 11c, when the excitation current is 1 A, the prediction errors of the M-Model, the D-Model, and the PD-Model are similar, and all smaller than the G-Model. However, when the excitation current increases to 6 A, the prediction error of the PD-Model is larger than that of the M-Model and the D-Model, as shown in Figure 11f.
Figure 12a–c presents the comparison between the predicted results of the ISSF-Duhem model with different shape functions as well as the PD-Model and the experimental tests in profiles of the displacement versus current, the displacement versus time, and the error versus time, respectively, when the magnetostrictive actuator is under the triangle excitation current with various amplitudes. Observing Figure 12, the ISSF-Duhem model can effectively predict the hysteresis of the magnetostrictive actuator under the triangle excitation and the PD-Model.
For the predicted results shown in Figure 11 and Figure 12, MSE, Emax, and Erelative between the predicted and the experimental displacements are calculated, as listed in Table 2. It is seen that the ISSF-Duhem model with different shape functions has different prediction accuracy for the hysteresis of the magnetostrictive actuator under different excitations. When the harmonic excitation current is 1 A, MSE of the D-Model of the predicted result is smaller than those of predicted results of the G-Model, M-Model, and PD-Model. Emax and Erelative of the predicted results of the D-Model, the M-Model, and the PD-Model are almost the same, but all lower than that of the G-Model. When the harmonic excitation current is 6 A, MSE of the D-Model predicted result is also smaller than those of the predicted results of the G-Model, the M-Model, and the PD-Model, but Emax and Erelative of the M-Model predicted result are the smallest. When the excitation current is altered to the triangle form as shown in Figure 6a, MSE, Emax, and Erelative of the D-Model predicted results are smaller than those of predicted results of the G-Model, the M-Model, and the PD-Model.
To sum up, MSE, Emax, and Erelative of the fitted and predicted results are comprehensively analyzed and shown in Figure 13. It can be seen from Figure 13a that MSE of the predicted result of the excitation current of 6 A is significantly larger than that of the fitted result of the excitation current of 5 A. The main reason is that the hysteresis increases with the excitation current increasing. As seen from Figure 13, the one-sided dead-zone operator-based shape function can enhance the performance of the ISSF-Duhem model as compared to the Grompertz function-based and the modified hyperbolic tangent function-based shape function. However, as shown in Figure 13b, the one-sided dead-zone operator-based shape function cannot make the performance of the ISSF-Duhem model as good as the modified hyperbolic tangent function-based shape function in some local condition. Based on the above comparison and analysis, the proposed modeling approach using the inertial system combining any shape function would be helpful for modeling asymmetric hysteresis of smart material-based actuators.

4. Discussion

The established ISSF-Duhem model focuses only on the asymmetric hysteresis of magnetostrictive actuators and the effect of stress and temperature on their hysteresis characteristics are not considered in this study. In fact, modeling the hysteresis with consideration of the effect of stress and temperature on hysteresis itself is an interesting and challenging task. It will be investigated in the future on the basis of this work.
As regarding the inverse problems of the nonlinear hysteresis models, it is actually to establish the mathematical relationship between the target output y d e s i r e and the required input u ( = f 3 y desire ). In order to inverse the ISSF-Duhem, firstly, the explicit functional relationship between w and u should be established according to the general expression of the first-order inertial system, w = f 4 u . Secondly, substituting the explicit functional relationship between w and u into the shape function to get the function relationship between z i i = 1 , 2 and u , z i i = 1 ,   2 = f 5 u . Lastly, the mathematical relationship between the target output y desire and the required input u , i.e., the inverse model of the ISSF-Duhem model is established. The general expression of the first-order inertial system shown as Equation (7) is a differential equation, which requires a discrete relationship in order to establish the functional relationship between w and u [40]. However, due to the coupling term between w and u in the ISSF-Duhem model, the explicit function relationship between w and u cannot be obtained, so the inverse model of the ISSF-Duhem model cannot be obtained. Big data and artificial intelligence methods might be used to describe nonlinear phenomena and to obtain the inverse model of the ISSF-Duhem model [41,42,43].

5. Conclusions

Based on the analyses of the Duhem model and the causes of hysteresis of magnetostrictive materials, introducing the concept of inertial system and shape function, the ISSF-Duhem model for describing the hysteresis of magnetostrictive actuators featured with saturation and asymmetry is proposed and investigated in this paper. Specifically, the hysteresis caused by the inertia in the motion of magnetic domains is described by a differential equation of the inertial system. The hysteresis curve represented by the inertial system is corrected by incorporation of a shape function for the hysteresis caused by the pinning of domain walls, the interactions between the magnetic domains and the saturation characteristics of the hysteresis.
The feasibility of the ISSF-Duhem model for describing the hysteresis is analyzed by solving the model. A multi-islands GA is applied for parameter identification for the ISSF-Duhem model. According to the identified parameters, the hysteresis of the magnetostrictive actuator is simulated and predicted by the ISSF-Duhem model. The comparison between the model simulated and predicted results and the experimental results show that the ISSF-Duhem model can effectively characterize the hysteresis of magnetostrictive actuators.
The shape function affects the performance of the ISSF-Duhem model. The comparative analysis shows that compared with the Grompertz function-based and the modified hyperbolic tangent function-based shape function, the one-sided dead-zone operator-based shape function could enhance the performance of the ISSF-Duhem model. In addition, the ISSF-Duhem model with the one-sided dead-zone operator-based and the modified hyperbolic tangent function-based shape function provide superior performance over the modified Prandtl–Ishlinskii model with a one-sided dead-zone operator.
Precise and simple hysteresis models of actuators are the foundation to realize the accurate, rapid, and effective control systems. The concept of the developed approach of the ISSF-Duhem model is of help and significance for hysteresis modeling of the magnetostrictive material-based actuators. Compared to the conventional Duhem model, the ISSF-Duhem model and the modeling approach with incorporation of the inertial system and shape function reduce the number of unknown functions. Only a suitable shape function is necessary to be selected according to the shape characteristics of the hysteresis curve to completely determine the model. As compared with the continuous-time dynamic model of a class of backlash-like hysteresis, the proposed ISSF-Duhem model has better flexibility in describing different hysteresis phenomena, and is more suitable for describing the hysteresis of magnetostrictive actuators featured with saturation and asymmetry, since it incorporates a shape function. Therefore, the modeling approach using an inertial system combining any shape function would be helpful for modeling asymmetric hysteresis of smart material-based actuators.

Author Contributions

Conceptualization, Z.-Y.S. and X.-X.’F’.B.; methodology, Z.-Y.S. and X.-X.’F’.B.; software, Z.-Y.S.; validation, Z.-Y.S., X.-X.’F’.B. and L.-J.Q.; formal analysis, Z.-Y.S.; investigation, Z.-Y.S. and X.-X.’F’.B.; data curation, Z.-Y.S.; writing—original draft preparation, Z.-Y.S.; writing—review and editing, X.-X.’F’.B.; supervision, L.-J.Q. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

References

  1. Claeyssen, F.; Lhermet, N.; Le Letty, R.; Bouchilloux, P. Actuators, transducers and motors based on giant magnetostrictive. J. Alloy. Compd. 1997, 258, 61–73. [Google Scholar] [CrossRef]
  2. Yi, S.; Yang, B.; Meng, G. Microvibration isolation by adaptive feedforward control with asymmetric hysteresis compensation. Mech. Syst. Signal Process. 2019, 114, 644–657. [Google Scholar] [CrossRef]
  3. Calkins, F.T.; Smith, R.C.; Flatau, A.B. Energy-based hysteresis model for magnetostrictive transducers. IEEE Trans. Magn. 2000, 36, 429–439. [Google Scholar] [CrossRef] [Green Version]
  4. Hassani, V.; Tjahjowidodo, T. Structural response investigation of a triangular-based piezoelectric drive mechanism to hysteresis effect of the piezoelectric actuator. Mech. Syst. Signal Process. 2013, 36, 210–223. [Google Scholar] [CrossRef]
  5. Minorowicz, B.; Leonetti, G.; Stefanski, F.; Binetti, G.; Naso, D. Design, modelling and control of a micro-positioning actuator based on magnetic shape memory alloys. Smart Mater. Struct. 2016, 25, 075005. [Google Scholar] [CrossRef]
  6. Zakerzadeh, M.R.; Sayyaadi, H. Precise position control of shape memory alloy actuator using inverse hysteresis model and model reference adaptive control system. Mechatronics 2013, 23, 1150–1162. [Google Scholar] [CrossRef]
  7. Zhou, M.; He, S.; Hu, B.; Zhang, Q. Modified KP model for hysteresis of magnetic shape memory alloy actuator. IETE Tech. Rev. 2015, 32, 29–36. [Google Scholar] [CrossRef]
  8. Zhu, W.; Yang, F.; Rui, X. Robust independent modal space control of a coupled nano-positioning piezo-stage. Mech. Syst. Signal Process. 2018, 106, 466–478. [Google Scholar] [CrossRef]
  9. Olabi, A.G.; Grunwald, A. Design and application of magnetostrictive materials. Mater. Des. 2008, 29, 469–483. [Google Scholar] [CrossRef] [Green Version]
  10. Tan, X.; Baras, J.S. Modeling and control of hysteresis in magnetostrictive actuators. Automatica 2004, 40, 1469–1480. [Google Scholar] [CrossRef] [Green Version]
  11. Jiles, D.C.; Atherton, D.L. Theory of ferromagnetic hysteresis. J. Magn. Magn. Mater. 1986, 61, 48–60. [Google Scholar] [CrossRef]
  12. Jiles, D.C. Frequency dependence of hysteresis curves in conducting magnetic materials. J. Appl. Phys. 1994, 76, 5849–5855. [Google Scholar] [CrossRef] [Green Version]
  13. Dapino, M.J.; Smith, R.C.; Faidley, L.E.; Flatau, A.B. A coupled structural-magnetic strain and stress model for magnetostrictive transducers. J. Intell. Mater. Syst. Struct. 2000, 11, 135–152. [Google Scholar] [CrossRef] [Green Version]
  14. Baghel, A.P.; Chwastek, K.; Kulkarni, S.V. Modeling of minor hysteresis loops in rolling and transverse directions of grain-oriented laminations. IET Electr. Power Appl. 2015, 9, 344–348. [Google Scholar] [CrossRef]
  15. Jiles, D.C. A self consistent generalized model for the calculation of minor loop excursions in the theory of hysteresis. IEEE Trans. Magn. 1992, 28, 2602–2604. [Google Scholar] [CrossRef]
  16. Valadkhan, S.; Morris, K.; Shum, A. A new load-dependent hysteresis model for magnetostrictive materials. Smart Mater. Struct. 2010, 19, 125003. [Google Scholar] [CrossRef]
  17. Zhu, Y.; Yang, X.; Wereley, N.M. Research on hysteresis loop considering the prestress effect and electrical input dynamics for a giant magnetostrictive actuator. Smart Mater. Struct. 2016, 25, 085030. [Google Scholar] [CrossRef]
  18. Li, J.W.; Xu, M.Q. Modified Jiles-Atherton-Sablik model for asymmetry in magnetomechanical effect under tensile and compressive stress. J. Appl. Phys. 2011, 110, 063918. [Google Scholar] [CrossRef]
  19. Iyer, R.V.; Tan, X.; Krishnaprasad, P.S. Approximate inversion of the Preisach hysteresis operator with application to control of smart actuators. IEEE Trans. Autom. Control 2005, 50, 798–810. [Google Scholar] [CrossRef]
  20. Mayergoyz, I. Mathematical models of hysteresis. IEEE Trans. Magn. 1986, 22, 603–608. [Google Scholar] [CrossRef] [Green Version]
  21. Kuhnen, K.; Krejci, P. Compensation of complex hysteresis and creep effects in piezoelectrically actuated systems—A new Preisach modeling approach. IEEE Trans. Autom. Control 2009, 54, 537–550. [Google Scholar] [CrossRef]
  22. Mayergoyz, I.D. Dynamic Preisach models of hysteresis. IEEE Trans. Magn. 1988, 24, 2925–2927. [Google Scholar] [CrossRef]
  23. Koutný, J.; Kružík, M.; Kurdila, A.J.; Roubíček, T. Identification of Preisach-Type Hysteresis Operators. Numer. Funct. Anal. Optim. 2008, 29, 149–160. [Google Scholar] [CrossRef]
  24. Al Janaideh, M.; Aljanaideh, O. Further results on open-loop compensation of rate-dependent hysteresis in a magnetostrictive actuator with the Prandtl-Ishlinskii model. Mech. Syst. Signal Process. 2018, 104, 835–850. [Google Scholar] [CrossRef]
  25. Feng, Y.; Li, Z.; Rakheja, S.; Jiang, H. A modified Prandtl-Ishlinskii hysteresis modeling method with load-dependent delay for characterizing magnetostrictive actuated systems. Mech. Sci. 2018, 9, 177–188. [Google Scholar] [CrossRef]
  26. Li, Z.; Su, C.Y.; Chen, X. Modeling and inverse adaptive control of asymmetric hysteresis systems with applications to magnetostrictive actuator. Control Eng. Pract. 2014, 33, 148–160. [Google Scholar] [CrossRef]
  27. Al Janaideh, M.; Su, C.Y.; Rakheja, S. Development of the rate-dependent Prandtl–Ishlinskii model for smart actuators. Smart Mater. Struct. 2008, 17, 035026. [Google Scholar] [CrossRef]
  28. Yang, M.J.; Li, C.X.; Gu, G.Y.; Zhu, L.M. A Modified Prandtl-Ishlinskii Model for Rate-dependent Hysteresis Nonlinearity Using m th-power Velocity Damping Mechanism. Int. J. Adv. Robot. Syst. 2014, 11, 163. [Google Scholar] [CrossRef] [Green Version]
  29. Kuhnen, K. Modelling, identification, and compensation of complex hysteretic and log (t)-type creep nonlinearities. Control Intell. Syst. 2005, 33, 134–147. [Google Scholar] [CrossRef] [Green Version]
  30. Lin, C.J.; Yang, S.R. Precise positioning of piezo-actuated stages using hysteresis-observer based control. Mechatronics 2006, 16, 417–426. [Google Scholar] [CrossRef] [Green Version]
  31. Zhu, W.; Bian, L.X.; Cheng, L.; Rui, X.T. Non-linear compensation and displacement control of the bias-rate-dependent hysteresis of a magnetostrictive actuator. Precis. Eng. 2017, 50, 107–113. [Google Scholar] [CrossRef]
  32. Bai, X.X.; Cai, F.L.; Chen, P. Resistor-capacitor (RC) operator-based hysteresis model for magnetorheological (MR) dampers. Mech. Syst. Signal Process. 2019, 117, 157–169. [Google Scholar] [CrossRef]
  33. Chen, P.; Bai, X.X.; Qian, L.J.; Choi, S.B. An approach for hysteresis modeling based on shape function and memory mechanism. IEEE/ASME Trans. Mechatron. 2018, 23, 1270–1278. [Google Scholar] [CrossRef]
  34. Coleman, B.D.; Hodgdon, M.L. On a class of constitutive relations for ferromagnetic hysteresis. Arch. Ration. Mech. Anal. 1987, 99, 375–396. [Google Scholar] [CrossRef]
  35. Macki, J.W.; Nistri, P.; Zecca, P. Mathematical models for hysteresis. Siam Rev. 1993, 35, 94–123. [Google Scholar] [CrossRef] [Green Version]
  36. Oh, J.; Bernstein, D.S. Semilinear Duhem model for rate-independent and rate-dependent hysteresis. IEEE Trans. Autom. Control 2005, 50, 631–645. [Google Scholar]
  37. Su, C.Y.; Stepanenko, Y.; Svoboda, J.; Leung, T.P. Robust adaptive control of a class of nonlinear systems with unknown backlash-like hysteresis. IEEE Trans. Autom. Control 2000, 45, 2427–2432. [Google Scholar] [CrossRef] [Green Version]
  38. Kciuk, M.; Chwastek, K.; Kluszczyński, K.; Szczygłowski, J. A study on hysteresis behaviour of SMA linear actuators based on unipolar sigmoid and hyperbolic tangent functions. Sens. Actuators A Phys. 2016, 243, 52–58. [Google Scholar] [CrossRef]
  39. Bai, X.X.; Chen, P.; Qian, L.J. Principle and validation of modified hysteretic models for magnetorheological dampers. Smart Mater. Struct. 2015, 24, 085014. [Google Scholar] [CrossRef]
  40. Li, W.; Chen, X.; Li, Z. Inverse compensation for hysteresis in piezoelectric actuator using an asymmetric rate-dependent model. Rev. Sci. Instrum. 2013, 84, 115003. [Google Scholar] [CrossRef]
  41. Nugroho, P.W.; Li, W.; Du, H.; Alici, G.; Yang, J. An adaptive neuro fuzzy hybrid control strategy for a semiactive suspension with magneto rheological damper. Adv. Mech. Eng. 2014, 6, 487312. [Google Scholar] [CrossRef]
  42. Weber, F.; Bhowmik, S.; Høgsberg, J. Extended neural network-based scheme for real-time force tracking with magnetorheological dampers. Struct. Control Health Monit. 2014, 21, 225–247. [Google Scholar] [CrossRef]
  43. Zong, L.H.; Gong, X.L.; Guo, C.Y.; Xuan, S.H. Inverse neuro-fuzzy MR damper model and its application in vibration control of vehicle suspension system. Veh. Syst. Dyn. 2012, 50, 1025–1041. [Google Scholar] [CrossRef]
Figure 1. Schematic of the physical mechanism of magnetostriction. (a) H = 0; (b) H ≠ 0 (c) H = Hs.
Figure 1. Schematic of the physical mechanism of magnetostriction. (a) H = 0; (b) H ≠ 0 (c) H = Hs.
Materials 13 02585 g001
Figure 2. The hysteresis characteristics of the inertial system with input: u = 1.5sin(2πt), when T = 2, 5 and 10.
Figure 2. The hysteresis characteristics of the inertial system with input: u = 1.5sin(2πt), when T = 2, 5 and 10.
Materials 13 02585 g002
Figure 3. The proposed ISSF-Duhem model.
Figure 3. The proposed ISSF-Duhem model.
Materials 13 02585 g003
Figure 4. Experimental setup for characteristic tests of the magnetostrictive actuator: (a) the photograph and (b) the schematic.
Figure 4. Experimental setup for characteristic tests of the magnetostrictive actuator: (a) the photograph and (b) the schematic.
Materials 13 02585 g004aMaterials 13 02585 g004b
Figure 5. Displacement–current relationship of the magnetostrictive actuator with different excitations: u(t) = Amsin(2πt), Am = 1, 2, 3, 4, 5, and 6 A.
Figure 5. Displacement–current relationship of the magnetostrictive actuator with different excitations: u(t) = Amsin(2πt), Am = 1, 2, 3, 4, 5, and 6 A.
Materials 13 02585 g005
Figure 6. Tests of the magnetostrictive actuator under triangle excitation current: (a) the excitation current with different amplitudes and (b) the output displacement.
Figure 6. Tests of the magnetostrictive actuator under triangle excitation current: (a) the excitation current with different amplitudes and (b) the output displacement.
Materials 13 02585 g006aMaterials 13 02585 g006b
Figure 7. The process of parameter identification for the models.
Figure 7. The process of parameter identification for the models.
Materials 13 02585 g007
Figure 8. Comparison of the experimental data and the fitted results of the established ISSF-Duhem model and the PD-Model when the magnetostrictive actuator under 1 Hz harmonic input with different currents: (a) 2 A, (b) 3 A, (c) 4 A, and (d) 5 A. G-Model: the ISSF-Duhem model with Grompertz function-based shape function; M-Model: the ISSF-Duhem model with modified hyperbolic tangent function-based shape function; D-Model: the ISSF-Duhem model with a one-sided dead-zone operator-based shape function; PD-Model: the modified Prandtl–Ishlinskii model with a one-sided dead-zone operator.
Figure 8. Comparison of the experimental data and the fitted results of the established ISSF-Duhem model and the PD-Model when the magnetostrictive actuator under 1 Hz harmonic input with different currents: (a) 2 A, (b) 3 A, (c) 4 A, and (d) 5 A. G-Model: the ISSF-Duhem model with Grompertz function-based shape function; M-Model: the ISSF-Duhem model with modified hyperbolic tangent function-based shape function; D-Model: the ISSF-Duhem model with a one-sided dead-zone operator-based shape function; PD-Model: the modified Prandtl–Ishlinskii model with a one-sided dead-zone operator.
Materials 13 02585 g008aMaterials 13 02585 g008b
Figure 9. Time history of the displacement difference between experimental results and the fitted results of the established ISSF-Duhem model and the PD-Model when the magnetostrictive actuator under 1 Hz harmonic input with different currents: (a) G-Model: the ISSF-Duhem model with Grompertz function-based shape function, (b) M-Model: the ISSF-Duhem model with modified hyperbolic tangent function-based shape function, (c) D-Model: the ISSF-Duhem model with a one-sided dead-zone operator-based shape function, and (d) PD-Model: the modified Prandtl–Ishlinskii model with a one-sided dead-zone operator.
Figure 9. Time history of the displacement difference between experimental results and the fitted results of the established ISSF-Duhem model and the PD-Model when the magnetostrictive actuator under 1 Hz harmonic input with different currents: (a) G-Model: the ISSF-Duhem model with Grompertz function-based shape function, (b) M-Model: the ISSF-Duhem model with modified hyperbolic tangent function-based shape function, (c) D-Model: the ISSF-Duhem model with a one-sided dead-zone operator-based shape function, and (d) PD-Model: the modified Prandtl–Ishlinskii model with a one-sided dead-zone operator.
Materials 13 02585 g009aMaterials 13 02585 g009b
Figure 10. Displacement difference between experimental data and the fitted results of the established ISSF-Duhem model and the PD-Model when the magnetostrictive actuator under 1 Hz harmonic input with different currents: (a) MSE versus current; (b) Emax and Erelative versus current. Emax: bars, Erelative: lines. G-Model: the ISSF-Duhem model with Grompertz function-based shape function; M-Model: the ISSF-Duhem model with modified hyperbolic tangent function-based shape function; D-Model: the ISSF-Duhem model with a one-sided dead-zone operator-based shape function; PD-Model: the modified Prandtl–Ishlinskii model with a one-sided dead-zone operator.
Figure 10. Displacement difference between experimental data and the fitted results of the established ISSF-Duhem model and the PD-Model when the magnetostrictive actuator under 1 Hz harmonic input with different currents: (a) MSE versus current; (b) Emax and Erelative versus current. Emax: bars, Erelative: lines. G-Model: the ISSF-Duhem model with Grompertz function-based shape function; M-Model: the ISSF-Duhem model with modified hyperbolic tangent function-based shape function; D-Model: the ISSF-Duhem model with a one-sided dead-zone operator-based shape function; PD-Model: the modified Prandtl–Ishlinskii model with a one-sided dead-zone operator.
Materials 13 02585 g010
Figure 11. Comparison of the experimental data and the predicted results of the established ISSF-Duhem model and the PD-Model when the magnetostrictive actuator under 1 Hz harmonic input with different currents: (a) the displacement versus current with the current of 1 A, (b) the displacement versus time with the current of 1 A, (c) the error versus time with the current of 1 A, (d) the displacement versus current with the current of 6 A, (e) the displacement versus time with the current of 6 A, and (f) the error versus time with the current of 6 A. G-Model: the ISSF-Duhem model with Grompertz function-based shape function; M-Model: the ISSF-Duhem model with modified hyperbolic tangent function-based shape function; D-Model: the ISSF-Duhem model with a one-sided dead-zone operator-based shape function; PD-Model: the modified Prandtl–Ishlinskii model with a one-sided dead-zone operator.
Figure 11. Comparison of the experimental data and the predicted results of the established ISSF-Duhem model and the PD-Model when the magnetostrictive actuator under 1 Hz harmonic input with different currents: (a) the displacement versus current with the current of 1 A, (b) the displacement versus time with the current of 1 A, (c) the error versus time with the current of 1 A, (d) the displacement versus current with the current of 6 A, (e) the displacement versus time with the current of 6 A, and (f) the error versus time with the current of 6 A. G-Model: the ISSF-Duhem model with Grompertz function-based shape function; M-Model: the ISSF-Duhem model with modified hyperbolic tangent function-based shape function; D-Model: the ISSF-Duhem model with a one-sided dead-zone operator-based shape function; PD-Model: the modified Prandtl–Ishlinskii model with a one-sided dead-zone operator.
Materials 13 02585 g011aMaterials 13 02585 g011bMaterials 13 02585 g011c
Figure 12. Comparison of experimental data and the predicted results of the established ISSF-Duhem model and the PD-Model when the magnetostrictive actuator under the triangle excitation current with various amplitudes: (a) the displacement versus current, (b) the displacement versus time, and (c) the error versus time. G-Model: the ISSF-Duhem model with Grompertz function-based shape function; M-Model: the ISSF-Duhem model with modified hyperbolic tangent function-based shape function; D-Model: the ISSF-Duhem model with a one-sided dead-zone operator-based shape function; PD-Model: the modified Prandtl–Ishlinskii model with a one-sided dead-zone operator.
Figure 12. Comparison of experimental data and the predicted results of the established ISSF-Duhem model and the PD-Model when the magnetostrictive actuator under the triangle excitation current with various amplitudes: (a) the displacement versus current, (b) the displacement versus time, and (c) the error versus time. G-Model: the ISSF-Duhem model with Grompertz function-based shape function; M-Model: the ISSF-Duhem model with modified hyperbolic tangent function-based shape function; D-Model: the ISSF-Duhem model with a one-sided dead-zone operator-based shape function; PD-Model: the modified Prandtl–Ishlinskii model with a one-sided dead-zone operator.
Materials 13 02585 g012aMaterials 13 02585 g012b
Figure 13. Displacement difference between the experimental results and the fitted, predicted results of the established ISSF-Duhem model and the PD-Model when the magnetostrictive actuator under 1 Hz harmonic input with different currents: (a) MSE versus current and (b) Emax and Erelative versus current. Emax: bars, Erelative: lines. G-Model: the ISSF-Duhem model with Grompertz function-based shape function; M-Model: the ISSF-Duhem model with modified hyperbolic tangent function-based shape function; D-Model: the ISSF-Duhem model with a one-sided dead-zone operator-based shape function; PD-Model: the modified Prandtl–Ishlinskii model with a one-sided dead-zone operator.
Figure 13. Displacement difference between the experimental results and the fitted, predicted results of the established ISSF-Duhem model and the PD-Model when the magnetostrictive actuator under 1 Hz harmonic input with different currents: (a) MSE versus current and (b) Emax and Erelative versus current. Emax: bars, Erelative: lines. G-Model: the ISSF-Duhem model with Grompertz function-based shape function; M-Model: the ISSF-Duhem model with modified hyperbolic tangent function-based shape function; D-Model: the ISSF-Duhem model with a one-sided dead-zone operator-based shape function; PD-Model: the modified Prandtl–Ishlinskii model with a one-sided dead-zone operator.
Materials 13 02585 g013aMaterials 13 02585 g013b
Table 1. Parameter identification results for ISSF-Duhem model with three shape functions under 1 Hz harmonic input with different currents. G-Model: the ISSF-Duhem model with Grompertz function-based shape function; M-Model: the ISSF-Duhem model with modified hyperbolic tangent function-based shape function; D-Model: the ISSF-Duhem model with a one-sided dead-zone operator-based shape function.
Table 1. Parameter identification results for ISSF-Duhem model with three shape functions under 1 Hz harmonic input with different currents. G-Model: the ISSF-Duhem model with Grompertz function-based shape function; M-Model: the ISSF-Duhem model with modified hyperbolic tangent function-based shape function; D-Model: the ISSF-Duhem model with a one-sided dead-zone operator-based shape function.
ParametersG-ModelM-ModelD-Model
k 1 3.0372.9502.979
k 2 34.92533.1020.778
α 0.5930.6120.575
β 3.4683.3221.859
μ 1 0.059
μ 2 0.988
μ 3 0.061
μ 4 1.251
μ 5 0.452
μ 6 2.378
μ 7 0.098
μ 8 0.587
μ 9 3.997
μ 10 0.214
n 6
Table 2. MSE, Emax and Erelative of the displacement difference between the predicted results of the ISSF-Duhem model and PD-Model with the identified parameters and the experimental tests when the magnetostrictive actuator is under 1 Hz harmonic input at different currents: 1 A and 6 A, and the triangle excitation current with various amplitudes. G-Model: the ISSF-Duhem model with Grompertz function-based shape function; M-Model: the ISSF-Duhem model with modified hyperbolic tangent function-based shape function; D-Model: the ISSF-Duhem model with a one-sided dead-zone operator-based shape function; PD-Model: the modified Prandtl–Ishlinskii model with a one-sided dead-zone operator. Tri: the triangle excitation current with various amplitudes.
Table 2. MSE, Emax and Erelative of the displacement difference between the predicted results of the ISSF-Duhem model and PD-Model with the identified parameters and the experimental tests when the magnetostrictive actuator is under 1 Hz harmonic input at different currents: 1 A and 6 A, and the triangle excitation current with various amplitudes. G-Model: the ISSF-Duhem model with Grompertz function-based shape function; M-Model: the ISSF-Duhem model with modified hyperbolic tangent function-based shape function; D-Model: the ISSF-Duhem model with a one-sided dead-zone operator-based shape function; PD-Model: the modified Prandtl–Ishlinskii model with a one-sided dead-zone operator. Tri: the triangle excitation current with various amplitudes.
ErrorMSEEmax (μm)Erelative (%)
Model 1 A6 ATri1 A6 ATri1 A6 ATri
G-Model0.1232.87832.82270.55832.93684.69297.203.465.50
M-Model0.01841.65291.12670.25041.97962.95263.232.343.46
D-Model0.00971.50470.92100.25842.50321.28613.242.951.53
PD-Model0.0142.18831.42720.26133.23832.37383.723.822.82

Share and Cite

MDPI and ACS Style

Si, Z.-Y.; Bai, X.-X.‘F.; Qian, L.-J. Asymmetric Hysteresis Modeling Approach Featuring “Inertial System + Shape Function” for Magnetostrictive Actuators. Materials 2020, 13, 2585. https://doi.org/10.3390/ma13112585

AMA Style

Si Z-Y, Bai X-X‘F, Qian L-J. Asymmetric Hysteresis Modeling Approach Featuring “Inertial System + Shape Function” for Magnetostrictive Actuators. Materials. 2020; 13(11):2585. https://doi.org/10.3390/ma13112585

Chicago/Turabian Style

Si, Zhi-Yuan, Xian-Xu ‘Frank’ Bai, and Li-Jun Qian. 2020. "Asymmetric Hysteresis Modeling Approach Featuring “Inertial System + Shape Function” for Magnetostrictive Actuators" Materials 13, no. 11: 2585. https://doi.org/10.3390/ma13112585

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