Modeling of Dynamic Hysteresis Based on Takagi-Sugeno Fuzzy Duhem Model

International Journal of Fuzzy Logic and Intelligent Systems.
2013.
Dec,
13(4):
277-283

This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which permits unrestricted noncommercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

- Received : November 21, 2013
- Accepted : December 25, 2013
- Published : December 30, 2013

Download

PDF

e-PUB

PubReader

PPT

Export by style

Article

Metrics

Cited by

TagCloud

In this study, we propose a novel method for modeling dynamic hysteresis. Hysteresis is a widespread phenomenon that is observed in many physical systems. Many different models have been developed for representing a hysteretic system. Among them, the Duhem model is a classical nonlinear dynamic hysteresis model satisfying the properties of hysteresis. The purpose of this work is to develop a novel method that expresses the local dynamics of the Duhem model by a linear system model. Our approach utilizes a certain type of fuzzy system that is based on Takagi-Sugeno (T-S) fuzzy models. The proposed T-S fuzzy Duhem model is achieved by fuzzy blending of the linear system model. A simulated example applied to shape memory alloy actuators, which have typical hysteretic properties, illustrates the applicability of our proposed scheme.
Hysteresis nonlinearity.
model
[9]
is a generalization of the Duhem model.
The hysteresis phenomenon is a typical nonlinear system. Hence, it is impossible to apply various linear system theories, and it is difficult to analyze a hysteresis system, even though the Duhem model expresses the dynamics of hysteresis well. To solve this problem, we propose a dynamic hysteresis model based on a Takagi-Sugeno (T-S) fuzzy model to express the local dynamics of a Duhem model with a linear system model. The proposed T-S fuzzy Duhem model is achieved by fuzzy blending of the linear system model. By using the proposed model, we can apply various linear system theories to a dynamic hysteresis model.
The remainder of this paper is organized as follows. In Section 2, we explain hysteresis and its properties. In Section 3, the dynamics of the Duhem model are introduced. Section 4 presents the proposed T-S fuzzy Duhem model. Simulated results verifying the applicability of the proposed model are provided in Section 5. Concluding remarks are presented in the final section.
Hysteretic memory.
as inherent to hysteresis or as necessary conditions for a system to display hysteresis. Two of the most important properties of hysteresis are memory and rate-independence.
Mayergoyz has said that hysteresis can be regarded as nonlinearity with memory. Memory implies that the output may depend not only on the input but also on the previous evolution of the input. According to Mayergoyz, all hysteresis nonlinearities can be classified into two categories: hysteresis nonlinearities with local memories and hysteresis nonlinearities with nonlocal memories. For the hysteresis nonlinearities with local memories, the future output depends uniquely upon the future input for any given output.
Figure 2
(a) shows that the graph can follow only one path if the input increases and only one path if the input decreases. However, for the hysteresis nonlinearities with nonlocal memories, the future output depends not only upon the current output and future input but also on the past history of input extremum values. Hence, the graph can trace any number of paths, depending on the previous evolution of input, as shown in
Figure 2
(b).
The decisive characteristic differentiating hysteretic nonlinearities from all other systems with memory is the rateindependence property. Hysteresis is said to be rate-independent if the rate of change of input has no influence on branching, while variation in the branches of the hysteresis nonlinearity are determined only by the past extremum of the input. These properties often determine whether a hysteresis model is suitable for representing a system or not
[5]
.
v
(
t
) and any initial value
w
_{0}
of the hysteresis output, the hysteresis output
w
is defined to be the solution of the following initial value problem:
The slopes of the hysteresis curves can be determined by slope functions
g
_{+}
and
g
_{－}
. The subscripts + and － represent increasing and decreasing curves, respectively. For (
ύ
(
t
))
^{±}
, the slope can be defined as
By using the Duhem model with the appropriate slope functions for both increasing and decreasing conditions, we can model the hysteresis. Consequently, how to choose suitable slope functions is key for successful modeling using the Duhem model. The literature
[1]
proposes the use of Gaussian probability density functions (PDFs) as slope functions for the Duhem model, which is given by
where
v
is the input,
k
is the constant parameter,
σ
^{2}
is the variance, and
μ
is the mean of the Gaussian PDF. We assume that hysteresis curves have zero slopes at both the beginning and end of each curve. Hence, we can also choose Gaussian PDFs for the slope functions of the Duhem model. From Eqs. (1-3), we can observe that the major loop slope function is given by the differential equation
An approach for representing the differential model of a minor loop was proposed by Likhachev. According to Likhachev, the major loop slope functions are multiplied by a scaling constant to obtain the minor loop slope functions. Consequently, the complete hysteresis model including both the major loop
Hysteresis curves based on the Duhem model.
and the minor loop can be expressed as
where
Figure 3
shows an example of hysteresis curves modeled by the Duhem model with Gaussian PDF slope functions. These curves consist of a major loop and a minor loop.
i
rules of the T-S fuzzy models are of the following form:
Here,
M_{ij}
is the fuzzy set, and
r
is the number of model rules:
x
(
t
)∈
ℜ
^{n}
is the state vector,
u
(
t
)∈
ℜ
^{m}
is the input vector,
y
(
t
)∈
ℜ
^{q}
is the output vector, and
A_{i}
∈
ℜ
^{n×n}
,
B_{i}
∈
ℜ
^{n×m}
, and
C_{i}
∈
ℜ
^{q×n}
;
z
_{1}
(
t
),…,
z_{p}
(
t
) are known premise variables that may be functions of the state variables and time.
By modifying the differential equation of the Duhem model, we can apply a T-S fuzzy model to the Duhem model. Eq. (5) can be modified as
We define
v
and
ύ
as premise variables
z
_{1}
(
t
) and
z
_{2}
(
t
), respectively. Then, the
i
rules of the proposed T-S fuzzy Duhem model are derived as
where
x
(
t
) =
w
(
t
) and
C_{i}
= 1.
There are two approaches for constructing a fuzzy model. The identification approach, which uses input-output data, is suitable for plants that are too difficult to be expressed as analytical models. On the other hand, when a nonlinear dynamic model is obtainable, a fuzzy model can be modeled by a derivation from the given nonlinear system equations. Indeed, if we know the nominal values of the parameters
μ
_{±}
and
σ
_{±}
, then we can use the second approach. In general, however, choosing the parameter values of the slope functions may be difficult. In such cases, the identification approach is more appropriate.
Input signal for major loop.
Input signal for mixed loop.
obtained from the literature
[16]
.
As mentioned above, we can model the hysteresis of an SMA actuator by using the classical Duhem model. Under the above parameters and equations, hysteresis curves modeled by the classical Duhem model are shown in
Figures 6
and
7
. Because an identification approach is chosen to construct T-S fuzzy systems, input-output pairs should be obtained. We can obtain input-output pairs from the simulation results for the classical Duhem model and construct a T-S fuzzy system with the data. Eleven fuzzy sets are defined as fuzzy sets of the first input
v
, as shown in
Figure 8
. To reflect two conditions of Eq. (9), two fuzzy sets are defined as fuzzy sets of the second input
ύ
, as shown in Figure 9. These fuzzy sets reflect whether the SMA actuator is in a cooling state or heating state. By using the input-output pairs and the previous membership functions, we can determine the values of
A_{i}
and
B_{i}
in Eq. (11) for each rule. Consequently, we can plot the hysteresis curves modeled by the T-S fuzzy Duhem model, as in
Figures 10
and
11
.
The results in
Figures 10
and
11
exhibit the hysteresis nonlinearity well. More analytically, we computed modeling errors between the results for the classical Duhem model and the T-S fuzzy Duhem model. The modeling errors are shown in
Figures 12
and
13
.
Simulation parameters
Hysteresis curves for the Duhem model with major loop.
Hysteresis curves for the Duhem model with mixed loop.
Membership functions for v .
The simulation results show that the proposed model approximates the Duhem model within a bounded range of [－1.2 1]. Considering the scaling constant
k
= 40, our proposed model has a modeling error within 3%. We conclude that this bounded range is acceptable.
Membership functions for ύ .
Hysteresis curves for the T-S fuzzy Duhem model with major loop.
Hysteresis curves for the T-S fuzzy Duhem model with mixed loop.
Modeling error for major loop case.
Modeling error for mixed loop case.
of hysteresis with a linear system. Additionally, the proposed model was verified to have applicability to simulations. Considering the wide range of hysteresis modeling, the proposed model is useful.

1. Introduction

Hysteresis is a phenomenon that has long been observed in a wide variety of physical systems, including mechanical hysteresis, magnetic hysteresis, and material hysteresis. For instance, hysteresis displayed by shape memory alloy (SMA) during the process of a phase transition is an example of hysteresis in a mechanical system
[1]
. There is no precise agreement in the literature regarding the definition of hysteresis. In this paper, we follow the definition proposed by Mayergoyz
[2]
. According to this definition, a system is said to be hysteretic if its input-output relationship is a multi-branch nonlinearity.
Figure 1
shows a typical hysteresis nonlinearity.
Modeling of hysteresis nonlinearity is an important goal of research in the field of hysteresis. For several decades, various models for hysteresis have been developed, such as the Preisach model
[3]
and Duhem model
[4]
. The Duhem model is a classical nonlinear dynamic model, whereas the Preisach model is not a dynamic model. Additionally, the Duhem model satisfies the properties of hysteresis, such as memory and rate-independence. Moreover, the Duhem model is able to tune the minor loop shapes independently from the major loop shapes. This increases the usefulness of the Duhem model in describing complicated minor loops. Hence, this model has been widely used to describe hysteretic systems, and many other models of hysteresis are related to the Duhem model
[5]
. For example, the Bouc-Wen model
[6
,
7]
as well as the Coleman-Hodgdon model
[8]
are special cases of the Duhem model. The Chua-Stromsmoe
Lager Image

2. Hysteresis and Its Properties

- 2.1 Hysteresis

Hysteresis arises in diverse physical systems, including mechanical hysteresis, magnetic hysteresis, and material hysteresis. The word “hysteresis” has been derived from the Greek word “hysterein,” which connotes lag, and hysteretic systems are generally described as having memory. Although there is no precise definition of hysteresis, our work follows the definition proposed by Mayergoyz. According to Mayergoyz, a system is said to be hysteretic if its input-output relationship is a multi-branch nonlinearity for which transitions from a branch to another branch occur at input extrema.
- 2.2 Properties of Hysteresis

Hysteresis has several properties. These properties are regarded
Lager Image

3. Classical Duhem Model

Many different mathematical models have been developed to express hysteresis nonlinearity, such as the Preisach model, Duhem model, Bouc-Wen model, and Chua-Stromsmoe model. Considering the usefulness of the dynamics provided by models and the capability of satisfying the above-mentioned properties, we chose a classical Duhem model to design the hysteresis model based on T-S fuzzy models.
The Duhem model is a differential-equation-based hysteresis model. For any differentiable hysteresis input
Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

4. T-S fuzzy Duhem Model

In Eqs. (5) and (6), the Duhem model is expressed in the form of differential dynamic equations. Because these dynamics consist of nonlinear functions, various linear system theories and analysis techniques cannot be applied to the hysteresis system. To solve this problem, we propose a T-S fuzzy Duhem model, which expresses the local dynamics of the Duhem model by a linear system model. As expressed in
[10]
, the
Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

5. Simulations

This section presents a simulation of the T-S fuzzy Duhem model with models of SMA actuators. SMAs are metallic alloys that display the shape memory effect. When SMAs are cooled or heated to a different temperature, phase transitions can occurred. During the process of the phase transition, SMAs exhibit a hysteretic transformation. To demonstrate the effectiveness of the proposed method, we compare the simulation results for modeling using the Duhem model with those using the T-S fuzzy Duhem model. Simulations for two cases wereexecuted to verify the applicability of the T-S fuzzy Duhem model. The first case is that the hysteresis curve displays only a major loop. The second case is that the hysteresis curve displays both a major loop and minor loop.
Figures 4
and
5
show the temperature input signal used in the simulations. The range of the input is [－300 300]. Table 1 shows the parameter values used in this modeling simulation. These parameters are
Lager Image

Lager Image

Simulation parameters

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

6. Conclusion

Many hysteresis models have been proposed and studied. However, the existing hysteresis models include nonlinear functions. Hence, it is impossible to apply various linear system theories and difficult to analyze a hysteretic system. To solve this problem, we propose a system based on a T-S fuzzy model. To design the T-S fuzzy system, the Duhem model, which is expressed by dynamic equations, can be used. Unlike the previous model, the proposed model can express the local dynamics
Lager Image

Lager Image

- Conflict of Interest

No potential conflict of interest relevant to this article was reported.
Dutta S. M.
,
Ghorbel F. H.
,
Dabney J. B.
“Modeling and control of a shape memory alloy actuator,”
Proceedings of the 20th IEEE International Symposium on Intelligent Control
Limassol, Cyprus
June 27-29, 2005
1007 -
1012
** DOI : 10.1109/.2005.1467151**

Mayergoyz I. D.
1991
“The classical preisach model of hysteresis,” in Mathematical Models of Hysteresis, I. D. Mayergoyz, Ed.
New York, NY
New York, NY
1 -
63
** DOI : 10.1007/978-1-4612-3028-1_1**

Mayergoyz I.
,
Friedman G.
1988
“Generalized preisach model of hysteresis,”
IEEE Transactions on Magnetics
24
(1)
212 -
217
** DOI : 10.1109/20.43892**

Macki J.W.
,
Nistri P.
,
Zecca P.
1993
“Mathematical models for hysteresis,”
SIAM Review
35
(1)
94 -
123
** DOI : 10.1137/1035005**

Dutta S. M.
,
M.S. thesis
2004
“Dynamic hysteresis modeling and applications,”
Rice University
Houston, TX
M.S. thesis

Bouc R.
1967
“Forced vibration of mechanical systems with hysteresis,”
Preceedings of the 4th International Conference on Nonlinear Oscillations
Prague, Czechoslovakia

Wen Y. K.
1976
“Method for random vibration of hysteretic systems,”
ASCE Journal of the Engineering Mechanics Division
102
(2)
249 -
263

Coleman B. D.
,
Hodgdon M. L.
1987
“On a class of constitutive relations for ferromagnetic hysteresis,”
Archive for Rational Mechanics and Analysis
99
(4)
375 -
396
** DOI : 10.1007/BF00282052**

Chua L. O.
,
Stromsmoe K. A.
1971
“Mathematical model for dynamic hysteresis loops,”
International Journal of Engineering Science
9
(5)
435 -
450
** DOI : 10.1016/0020-7225(71)90046-2**

Tanaka K.
,
Wang H. O.
2001
Fuzzy Control Systems Design and Analysis: a Linear Matrix Inequality Approach
Wiley
New York, NY

Wang L. X.
1997
A Course in Fuzzy Systems and Control
Prentice Hall PTR
Upper Saddle River, NJ

Jun S. W.
,
Kim D. W.
,
Lee H. J.
2011
“Design of T-S fuzzy-model-based controller for control of autonomous underwater vehicles,”
Journal of Korean Institute of Intelligent Systems
21
(3)
302 -
306
** DOI : 10.5391/JKIIS.2011.21.3.302**

Kim H. J.
,
Joo Y. H.
,
Park J. B.
202
“Controller design for continuous-time Takagi-Sugeno fuzzy systems with fuzzy lyapunov functions: LMI approach,”
International Journal of Fuzzy Logic and Intelligent Systems
12
(3)
187 -
192
** DOI : 10.5391/IJFIS.2012.12.3.187**

Song M. K.
,
Park J. B.
,
Kim J. K.
,
Joo Y. H.
2009
“Delayrange-dependent stability analysis and stabilization for nonlinear systems: T-S fuzzy model approach,”
Journal of Korean Institute of Intelligent Systems
19
(3)
337 -
342
** DOI : 10.5391/JKIIS.2009.19.3.337**

Lee H. J.
,
Joo Y. H.
,
Lee S. Y.
,
Park J. B.
2001
“Stochastic stabilization of TS fuzzy system with Markovian input delay,”
Journal of Korean Institute of Intelligent Systems
11
(6)
459 -
464

Feng Y.
,
Rabbath C. A.
,
Su C. Y.
“Inverse Duhem model based robust adaptive control for flap positioning system with SMA actuators,”
Proceedings of the 18th IFAC World Congress
Milano, Italy
August 28-September 2, 2011
8126 -
8131
** DOI : 10.3182/20110828-6-IT-1002.01744**

Citing 'Modeling of Dynamic Hysteresis Based on Takagi-Sugeno Fuzzy Duhem Model
'

@article{ E1FLA5_2013_v13n4_277}
,title={Modeling of Dynamic Hysteresis Based on Takagi-Sugeno Fuzzy Duhem Model}
,volume={4}
, url={http://dx.doi.org/10.5391/IJFIS.2013.13.4.277}, DOI={10.5391/IJFIS.2013.13.4.277}
, number= {4}
, journal={International Journal of Fuzzy Logic and Intelligent Systems}
, publisher={Korean Institute of Intelligent Systems}
, author={Lee, Sang-Yun
and
Park, Mignon
and
Baek, Jaeho}
, year={2013}
, month={Dec}