Fuzzy sliding-mode control with rule adaptation for nonlinear systems

Report 2 Downloads 167 Views
Article _____________________________ Fuzzy sliding-mode control with rule adaptation for nonlinear systems Lon-Chen Hung and Hung-Yuan Chung Department of Electrical Engineering, National Central University, Jhong-Li, Tao-Yuan 320, Taiwan, Republic of China E-mail: [email protected]

Abstract: A fuzzy sliding-mode control with rule adaptation design approach with decoupling method is proposed. It provides a simple way to achieve asymptotic stability by a decoupling method for a class of uncertain nonlinear systems. The adaptive fuzzy sliding-mode control system is composed of a fuzzy controller and a compensation controller. The fuzzy controller is the main rule regulation controller, which is used to approximate an ideal computational controller. The compensation controller is designed to compensate for the difference between the ideal computational controller and the adaptive fuzzy controller. Fuzzy regulation is used as an approximator to identify the uncertainty. The simulation results for two cart–pole systems and a ball– beam system are presented to demonstrate the effectiveness and robustness of the method. In addition, the experimental results for a tunnelling robot manipulator are given to demonstrate the effectiveness of the system.

Keywords: fuzzy, sliding-mode control, decoupling, adaptive, robot

1. Introduction Fuzzy logic control, as one of the most useful approaches for utilizing expert knowledge, has undergone extensive research in the past decade (Lim & Hiyama, 1991; Lau & Wong, 1998; Lee, 2002; Pham & Chen, 2002; Leung et al., 2003). Fuzzy logic control is generally applicable to plants that are mathematically poorly modelled and where experienced operators are available for providing qualitative guidance. Although achieving much practical success, fuzzy control has not been viewed as a rigorous science, due to a lack of formal synthesis techniques which guarantee the very basic requirements of global stability and acceptable performance (Wang, 1994; Kazemian, 2001). In stability analysis (Filev & Yager, 1994; Chai & Tong, 1999), it is commonly assumed that the mathematical 226

Expert Systems, September 2006, Vol. 23, No. 4

model of the plant is known; this assumption contradicts the very basic premise of fuzzy control systems, i.e. to control processes that are poorly modelled from a mathematical viewpoint. Based on fuzzy systems which are capable of approximating, with arbitrary accuracy, any real continuous function on a compact set, a globally stable adaptive controller is first synthesized from a collection of fuzzy IF–THEN rules (Lee, 1990; Wu et al., 2003). The fuzzy system, used to approximate an optimal controller, is adjusted by an adaptive law based on a Lyapunov function synthesis approach. In recent years, there have been attempts to design a fuzzy system based on the sliding-mode control law (Hwang & Lin, 1992; Palm, 1994; Kim & Lee, 1995; Glower & Munighan, 1997; Chen, 2001; Pham & Chen, 2002; Liang & Su, 2005). These authors have shown that the

 c 2006 The Authors. Journal Compilation  c 2006 Blackwell Publishing Ltd.

boundary layer can be reached in finite time and the ultimate boundedness of states is obtained asymptotically even though there exists some disturbance of dynamic uncertainties of the system. Palm (1994) showed the analogy between a simple and a sliding-mode controller with a boundary layer. Hwang and Li (1992) proposed a fuzzy sliding-mode controller and opened a way of designing for higher order nonlinear systems. Kim and Lee (1995) used some fuzzy control rules to construct reaching control under the assumption that equivalent control already exists. The Liang and Su (2005) design of a new two-input, one-output fuzzy sliding-mode controller is based on its distinct capability of improving the transient behaviour of the system during the reaching phase as well as the guaranteed steady-state tracking precision while in the boundary layer. In most studies, the fuzzy controller of secondorder systems is designed on a phase plane built by the error and change of error that are produced from the states and change of states. For example, in a cart–pole system only the pole subsystem is considered, ignoring the cart subsystem, and it is thus impossible to achieve a good control around the set point. In this study, a decoupling fuzzy controller design is proposed. This controller guarantees some properties, such as the robust performance and stability properties. Further, a class of fourth-order nonlinear systems is investigated. Lo and Kuo (1998) proposed a decoupled fuzzy sliding-mode control (DFSMC) to cope with the above issue. However, the method focuses on giving fuzzy rules, under the assumption that equivalent control exists. An adaptive decoupling fuzzy sliding-mode control (ADFSMC) design scheme is presented through width adaptation for a class of fourthorder uncertain under-actuated systems. Each subsystem, which decouples into two secondorder systems, is said to have main and subcontrol purposes. Two sliding surfaces are constructed through the state variables of the decoupling subsystem. We define the main and subtarget conditions for these sliding surfaces and introduce an intermediate variable from the sub-sliding surface condition. The proposed

adaptation law, which results from the direct adaptive approach, is used to appropriately determine the width of the unknown system variables. And the membership functions in the THEN part vary with the width adaptation of consequence. An adaptive law is then used to tune width in the THEN part to appropriately determine the distribution of each membership function. The proposed control is that the structure of the fuzzy controller does not need to be changed while using a common design procedure and the computing time may be reduced considerably. A tuning methodology is derived to tune the consequence parts of the fuzzy rules. The online tuning algorithm is derived in the Lyapunov sense; thus, the stability of the control system can be guaranteed. Furthermore, to relax the necessity for the uncertain bound of the compensation controller, an estimation mechanism is investigated. The uncertain bound is investigated so that the chattering phenomena of the control efforts can be decreased. To illustrate the effectiveness of the proposed design method, a comparison between a DFSMC and the proposed ADFSMC is made. The proposed ADFSMC has the following advantages. (1) With it we can control most complex systems without knowing the exact mathematical models. (2) The dynamic behaviour of the controlled system can be approximately represented by a decoupling sliding surface. (3) The proposed control approach has an advantage over the former model-based control schemes in that it does not require prior knowledge of the dynamic nonlinear system. (4) The ADFSMC not only can increase the robustness to system uncertainties but also can decrease the chattering phenomenon in the conventional sliding-mode controller. The paper is divided into five sections. In Section 2, the systems and the problem formulation are described. In Section 3, the ADFSMC is presented. In Section 4, the proposed controller is used to control a cart–single-pendulum system, a cart–double-pendulum system and a ball–beam system, and we show how the controller is synthesized as a Borland Cþþ Builder

 c 2006 The Authors. Journal Compilation  c 2006 Blackwell Publishing Ltd.

Expert Systems, September 2006, Vol. 23, No. 4

227

code and integrated into the software executed by the computer that controls a tunnelling robot manipulator system. Finally, we conclude with Section 5.

s1 ¼ c1 ðx1  zÞ þ x2

2. System description Consider a second-order nonlinear system, which can be represented by the following state-space model in canonical form:

x_ 2 ðtÞ ¼ f ðxÞ þ bðxÞu þ dðtÞ

¼ ½c1 1½x1 x2 T  c1 z ¼ cT x12  c1 z

ð3Þ

and then the other switching line as

x_ 1 ðtÞ ¼ x2 ðtÞ

s2 ¼ c 2 x 3 þ x 4

ð1Þ

yðtÞ ¼ x1 ðtÞ where x ¼ [x1 x2]T is the state vector, f(x) and b(x) are nonlinear functions, u is the control input and d(t) is an external disturbance. The disturbance is assumed to be bounded as follows: |d(t)|rD(t). For this kind of second-order system, we can use many kinds of control methods, such as proportional integral derivative control, sliding-mode control, fuzzy control etc. A control law u can be easily designed to make the secondorder system (1) arrive at our control goal. However, for such nonlinear models as an inverted pendulum system, the system dynamic representation is generally not in an exact canonical form. Rather, it has the following form: x_ 1 ðtÞ ¼ x2 ðtÞ

ð4Þ

In the design of a decoupling sliding-mode controller, an equivalent control is first given so that the states can stay on the sliding surface. Thus, in sliding motion, the system dynamic is independent of the original system and a stable equivalent control system is achieved. The equivalent control can be obtained by letting s_ 1 equal zero, i.e. _ þ x_ 2 s_1 ¼ c1 ðx_ 1  zÞ ¼ c1 x2  c1 z_ þ f1 þ b1 u þ d1 ¼ 0

ð5Þ

The decoupling sliding-mode control input is chosen as follows for a Lyapunov function candidate: 1 V ¼ s1 s21 2

ð6Þ

Taking the time derivative of (6), we have

x_ 2 ðtÞ ¼ f1 ðxÞ þ b1 ðxÞu1 þ d1 ðtÞ x_ 3 ðtÞ ¼ x4 ðtÞ

ð2Þ

V_ ¼ s1 s_ 1 ¼ s1 ðc1 x2  c1 z_ þ f1 þ b1 u þ d1 Þ

x_ 4 ðtÞ ¼ f2 ðxÞ þ b2 ðxÞu2 þ d2 ðtÞ T

where x ¼ [x1 x2 x3 x4] is the state vector, f1(x), f2(x) and b1(x), b2(x) are nonlinear functions, u1, u2 are the control inputs and d1(t), d2(t) are external disturbances. The disturbances are assumed to be bounded as follows: |d1(t)|rD1(t), |d2(t)| rD2(t). One can use (2) effectively to design u1 and u2; however, this approach is only utilized to control a subsystem of (2). For example, if the model is a cart–pole system, we only control either the pole or the cart of a system such as (2). Hence, the idea of decoupling 228

is employed to design a control u to govern the whole system. In equation (2), we first define one switching line as

Expert Systems, September 2006, Vol. 23, No. 4

ð7Þ

It can easily be shown from (7) that the decoupling sliding-mode controller u1 can be divided into an equivalent control input and a reaching control input which, if it has the following form, will be negative: u ¼ ueq  K sgnðs1 =F1 Þ

ð8Þ

where K > D1(t)=|b1(t)|. K is a positive constant. Then the system is controlled in such a way that the state always moves towards the sliding surface and hits it. Thus, the trajectory is always forced to move towards the sliding surface. But

 c 2006 The Authors. Journal Compilation  c 2006 Blackwell Publishing Ltd.

(8) will show high-frequency switching near the sliding surface (s1 ¼ 0) due to the sgn function involved. Thus, in order to reduce the chattering phenomenon, we replace sgn(s1) with sat(s1) as follows: u ¼ ueq  K satðs1 =F1 Þ

ð9Þ

where K > D1(t)=|b1(t)|. Hence, in the sliding motion, an equivalent controller will be ueq ¼

1 ðc1 x2 þ c1 z_ f1 þ z_ 1 þ ls1 Þ b1

ð10Þ

Substituting equation (10) into equation (5), we obtain s_ 1 þ ls1 ¼ 0

ð11Þ

where l is a positive value. The sliding surface on the phase plane can be defined as (3). The control objective is to drive the system state to the original equilibrium point. The switching line variables s1 and s2 are reduced to zero gradually at the same time by an intermediate variable z. In equation (3), z is a value transferred from s2; it has a value proportional to s2 and has the range proper to x1. Equation (3) denotes that the control objective of u1 is changed from x1 ¼ 0, x2 ¼ 0 to x1 ¼ z, x2 ¼ 0 (Lo & Kuo, 1998). Because the controller u ¼ u1 is used to govern the whole system, the bound of x1 can be guaranteed by letting jzjrZupper

0 < Zupper < 1

ð12Þ

where Zupper is the upper bound of abs(z). Equation (12) implies that the maximum absolute value of x1 will be limited. Summarizing what we have mentioned above, z can be defined as z ¼ satðs2 =Fz Þ Zupper

0 < Zupper < 1 ð13Þ

where Fz is the boundary layer of s2 to smooth z, Fz transfers s2 to the proper range of x1, and the definition of the sat(.) function is  sgnðjÞ if jjjZ1 ð14Þ satðjÞ ¼ j if jjj< 1

Notice that z is a decaying oscillation signal because Zupper is a factor less than one. Remark. Decouple the whole system into subsystems A and B. System A contains x1, x2 and its corresponding sliding surface s1 ¼ c1 x1 þ x2. System B contains x3, x4 and its corresponding sliding surface s2 ¼ c2x3 þ x4. By doing this, the main control objective is to keep the states of system A moving toward the surface s1 ¼ 0 and converging to x1 ¼ x2 ¼ 0 asymptotically, and a subtarget is to keep the states of system B moving toward the surface s2 ¼ 0 and sliding to x3 ¼ x4 ¼ 0 asymptotically. Because the main target is to stabilize system A, it is reasonable to consider the information from system B as secondary and this secondary information must be reflected through a mechanism to the primary. An intermediate variable z, which represents secondary information, is incorporated into s1 and therefore the sliding surface is modified to take the form c1(x1  z) þ x2. This modification reflects the fact that the main target is changed from x1 ¼ 0, x2 ¼ 0 to x1 ¼ z, x2 ¼ 0, where x1 is a function of z. Notice that the subtarget of system B, s2 ¼ 0, is embedded in the main target through the variable z. Both primary and secondary targets can be controlled simultaneously. z represents a decaying oscillation signal whose value is assigned through (13). Hence, consider equation (3). If s1 ¼ 0, then x1 ¼ z, x2 ¼ 0. Since z is a value transferred from s2, when s2 ! 0, then z ! 0 and x1 ! 0. From equation (4), if the condition s1 ! 0 holds, the control objective can be achieved. Moreover, the choice of c1 and c2 has a strong influence on the behaviour in the transient state of the system. Appropriate choice of c1 and c2 is necessary for achieving a favourable transient response. Hence, control implementation can be started with zero initial fuzzy rules to eliminate the trial-and-error process for designing fuzzy control rules and has been developed with rule bases that do not need to be defined a priori.

 c 2006 The Authors. Journal Compilation  c 2006 Blackwell Publishing Ltd.

Expert Systems, September 2006, Vol. 23, No. 4

229

If the perfect control law does not obtain, it is not possible to implement (9). To overcome such a problem, a novel approach of the rule adaptation of fuzzy control is proposed to estimate an equivalent control input ueq. This will be proposed in Section 3.

this term is to force the system state to slide on the sliding surface. Another is the corrective control uc that drives the states towards the sliding surface. Thus the control law can be represented as

3. Design of ADFSMC

where uADFSMC is the approximate equivalent control and the corrective control uc is designed to stabilize the states of the control system around a preselected uncertainty bound. Now, the rule base of ADFSMC is constructed as

In this section, we show how to develop an adaptive fuzzy sliding-mode controller for obtaining equivalent control through rule adaptation in the decoupling system. Then, we construct corrective control to guarantee the system’s stability. The proposed ADFSMC has online self-tuning fuzzy rules without the trial-and-error process for finding appropriate fuzzy rules. If the state trajectory can be forced to slide on the sliding surface, then a stable equivalent control system is achieved. However, if the function f1 is unknown (for simplicity, we assume b1 is known), there is no way to yield equivalent control ueq. In this paper, a set of fuzzy rule bases is applied to approximating (10). Motivated by the principle of decoupling sliding-mode control, the control law consists of two parts. One is the estimated decoupling sliding component uADFSMC that is constructed by an adaptive mechanism, as shown in Figure 1. The effect of

u ¼ uADFSMC þ uc

ð15Þ

(i; j)th rule: If s1 is Ai and s_1 is Bj then uADFSMC is uk where i 2 I ¼ fp; p þ 1; . . . ; 1; 0; 1; ð16Þ

. . . ; p  1; pg j 2 J ¼ fq; q þ 1; . . . ; 1; 0; 1;

ð17Þ

. . . ; q  1; qg k 2 K ¼ fr; r þ 1; . . . ; 1; 0; 1;

ð18Þ

. . . ; r  1; rg

The rule base of ADFSMC is constructed as in this study; triangular types and singletons,

x3 S2 SMC x4

Adaptive Law (39)

Z x1 S1 SMC

Adaptive Fuzzy Controller

x2

Corrective Controller (40),(41)

uADFSMC +

u Plant +

y

uc

Figure 1: Block diagram of the ADFSMC system. 230

Expert Systems, September 2006, Vol. 23, No. 4

 c 2006 The Authors. Journal Compilation  c 2006 Blackwell Publishing Ltd.

respectively, are used to define the membership functions of the IF part and the THEN part, which are shown in Figures 2–4. The membership functions of the IF part and the THEN part are arranged to have the same width and boundary in the universe of discourse [  p; p], [  q; q] and [  U; U] respectively. The rule is defined in the following analytical form (Chen, 2001): k ¼ ½si  ð1  sÞj

s 2 ½0; 1

ð19Þ

where k is an operation that takes an integer which is the nearest to k and s is a rule regulating factor. Obviously, by properly adjusting s, the value of k will be changed by (19), which indirectly determines which uk should be taken into account. So the controller is expected to provide different control actions corresponding to different s. According to (16)–(18), it is easily found that (19) is constrained by rrsi þ ð1  sÞjrr

ð20Þ

 A− p

Ap +1

A0

Ap − 1

u−r

ur

uADFSMC

−U

U

0

Figure 4: The membership functions of the THEN part.

Considering the extreme case, i ¼ p and j ¼ q, the following condition shows that s should be satisfied:

sr

rq pq

ð21Þ

It also shows that r ¼ p > q satisfies the condition s 2 [0, 1]. Define each membership function of the THEN part as (Chen, 2001) ð22Þ

where h ¼ U=r and kr is the modified function of k and is represented as

p

B0

kr ¼  si  ð1  sÞj Considering the case pirs1rpiþ1 qi r s_1 rqiþ1 , four rules are fired:

 Bq − 1

ð23Þ and

Bq

ði; jÞ

s1 0

ur−1

h

Figure 2: The membership functions of the input variable s1 of the IF part.

−q

u0

uk ¼ kr h

0

B− q B− q+ 1

u−r+1

Ap

s1 −p



ði; j þ 1Þ

ði þ 1; jÞ

ði þ 1; j þ 1Þ

ð24Þ

Therefore, we can get a crisp output through the defuzzification strategy

q

Figure 3: The membership functions of the input variable s_ 1 of the IF part.  c 2006 The Authors. Journal Compilation  c 2006 Blackwell Publishing Ltd.

uADFSMC ¼

O X

Expert Systems, September 2006, Vol. 23, No. 4

ð25Þ

231

where

where O ¼ w1 u1 þ w2 u2 þ w3 u3 þ w4 u4

y ¼ ½w1 ði þ jÞ þ w2 ði þ j  1Þ

¼ w1 ½si  ð1  sÞjh þ w2 ½sði þ 1Þ

þ w3 ði þ j  1Þ þ w4 ði þ jÞh

 ð1  sÞjh þ w3 ½si  ð1  sÞðj þ 1Þh

ð31Þ

The control law for the ADFSMC system is assumed to have the following form:

þ w4 ½sði þ 1Þ  ð1  sÞð1 þ jÞh ¼ fw1 ½sði þ jÞ  j þ w2 ½sði þ j  1Þ  j

u ¼ uADFSMC þ uc

ð32Þ

þ w3 ½sði þ j  1Þ  ðj þ 1Þ þ w4 ½sði þ jÞ  ð1 þ jÞgh

X ¼ w1 þ w2 þ w3 þ w4

ð26Þ

ð27Þ

and

where uADFSMC is the approximate equivalent control and the corrective control uc is designed to stabilize the states of the control system around a preselected uncertainty bound. Substituting equation (32) into equation (2), we obtain x_ 2 ¼ f1 ðxÞ þ b1 ðxÞu

w1 ¼ min½mAi ðs1 Þ; mBj ðs_1 Þ w2 ¼ min½mAiþ1 ðs1 Þ; mBj ðs_1 Þ

¼ f1 ðxÞ þ b1 ðxÞðuADFSMC þ uc Þ ð28Þ

¼ c1 x2 þ c1 z_ þ b1 ðxÞðuADFSMC þ uc  ueq Þ

w3 ¼ min½mA ðs1 Þ; mBjþ1 ðs_1 Þ

ð33Þ or, equivalently,

w4 ¼ min½mAiþ1 ðs1 Þ; mBjþ1 ðs_1 Þ The main aim is to derive an adaptive law to adjust the regulating factor s such that the estimated equivalent control uADFSMC can be optimally approximated to the equivalent control of the decoupling sliding-mode control under the situation of an unknown function f1. For simplicity, a triangular type membership function is chosen for the aforementioned fuzzy variables. The online tuning algorithm of parameters is proposed to adjust the consequent parameters for monitoring the system control performance. Suppose there exists an optimal regulating factor s* which is a constant such that the uADFSMC has minimum approximation error: e ¼ uADFSMC  uADFSMC

ð29Þ

Then from (19)–(21) we have uADFSMC  232

uADFSMC

1 ¼ ðs  s Þy X

Expert Systems, September 2006, Vol. 23, No. 4

x_ 12 ¼ Ac x12 þ Bc ðuADFSMC þ uc  ueq Þ þ Pz_ þ Qs1

ð34Þ

where " x12 ¼ ½ x1 Bc ¼ ½ 0 P¼½0

x2 

T

Ac ¼

0

1

0

c1

#

b1 ðxÞ T c1  T

Q¼½0

l T

Equation (5) can be rewritten s_D ¼ cT x12  c1 z ¼ cT Ac x12 þ cT Bc ðuADFSMC þ uc  ueq Þþ cT Pz þ cT Qs1  c1 z_ ¼ b1 ðuADFSMC þ uc  ueq Þ  lsD ð35Þ

ð30Þ

where c ¼ ½ c1

1 T and sD ¼ s1  F1 sat(s1=F1).

 c 2006 The Authors. Journal Compilation  c 2006 Blackwell Publishing Ltd.

The above properties of the boundary layer concept are exploited in the design of the ADFSMC system, our goal being rule adaptation as soon as the boundary layer is reached. This approach aims to avoid the possibility of unbounded growth. As described above we develop an adaptive fuzzy control law with rule regulation in an attempt to accomplish this objective. The adjustment method for the fuzzy rule is local, and therefore only the specific operation domain is adapted without detrimental effects on other operating areas. In this paper, a method is proposed to guarantee the closed-loop stability by the Lyapunov function. In order to analyse the overall system stability, we introduce assumptions as follows. Assumption 1: b1(x) in (2) is differentiable, i.e.

(32), in which the adaptation laws of the fuzzy controller are designed as (39) and the corrective controller is designed as (40) with the adaptive bound estimation shown in (41), then we can guarantee the asymptotic stability of the closedloop system. The fuzzy rule adaptive laws are given by

ð40Þ

K_^ ¼ Z2 jsD j

ð41Þ



ð36Þ

Assumption 2: There exist optimal values for the fuzzy rule such that      1 @b1 1  ¼ K u _ s  u þ ð37Þ x eq D  ADFSMC 2 @x 

1 2 1 2 1 2 C þ L s þ 2jb1 j D 2Z1 2Z2

ð42Þ

^  s , L ¼ K^ K  and Z1, Z2 are where C ¼ s positive constants. The variation of this function (42) with respect to time is sD s_D 1 2 @b1 1 _ þ 1 L L_ V_ ¼ þ sD sgnðb1 Þ 1 x_ þ C C Z1 Z2 jb1 j 2 @x

where the uncertainty bound K  is a positive constant.

¼

sD s_D 1 @b1 1 _ 1 ^ þ L K_^ þ s_D sgnðb1 Þ 1 x_ þ C s Z1 Z2 jb1 j 2 @x

¼

1 1 @b1 _ þ s2D sgnðb1 Þ 1 x_ sD ðcT x_ 12  c1 zÞ jb1 j 2 @x

This uncertainty bound cannot be measured for practical applications. Therefore, a bound estimation is developed to observe the bound of approximation error: ^  K X ¼ KðtÞ

uc ¼ K sgnðb1 Þ satðs1 =F1 Þ

Proof: Choose the Lyapunov function as

where  1  @b1 @b1 @b1 ; 1 ; :::; 1 @x1 @x2 @x4

ð39Þ

where Z1 and Z2 are positive constants. Moreover, the system states converge to the sliding surface asymptotically.

d 1 dx _ b ðxÞ ¼ Db1 ¼ Db1 1 ðxÞ 1 ðxÞðxÞ dt 1 dt

Db1 1 ðxÞ ¼

1 ^_ ¼ s Z sD sgnðb1 Þ y X 1

ð38Þ

^ is the estimated uncertainty bound. where KðtÞ Theorem 1: Considering the dynamic system described by (2) and the sliding surface (3) for the bounded continuous desired state trajectory, if the ADFSMC law is designed as  c 2006 The Authors. Journal Compilation  c 2006 Blackwell Publishing Ltd.

þ

¼

1 1 ^_ þ L K_^ Cs Z1 Z2

 l 2 sD þ sD sgnðb1 Þ uADFSMC þ uc  ueq jb1 j 1 @b1 þ sD 1 x_ þuADFSMC  uADFSMC 2 @x þ



1 1 ^_ þ L K_^ Cs Z1 Z2 Expert Systems, September 2006, Vol. 23, No. 4

233

¼

l 2 s þ sD sgnðb1 Þ jb1 j D   1 @b1  1  uADFSMC  ueq þ sD x_ 2 @x

þ sD sgnðb1 ÞðuADFSMC  uADFSMC Þ 1 _ 1 ^ þ L K_^ þ sD sgnðb1 Þ uc þ C s Z1 Z2 l 2 1 s þ jsD jK  þ sD sgnðb1ÞCyjsD jK^ r jb1 j D X 1 _ 1 ^ þ L K_^ þ Cs Z1 Z2   l 2 1 _ 1 ^ þ sD sgnðb1 Þ y s ¼ sD þ C T jb1 j Z1 X   1 _^ þL ð43Þ K jsD j Z2 By selecting appropriate values for F1, (43) implies that V_ is negative semidefinite: l 2 s V_ r jb1 j D

ð44Þ

If |sD|rF1, sD ¼ 0, then V ¼ 0 and V_ ¼ 0. If |sD| > F1, then sD ¼ s1  F1 sat(s1=F1) and s_D ¼ s_1 has the same sign as s1. From the algorithm, we have s1  s_1 < 0; therefore V_ ¼ s1  s_1 < 0. Then, for all t Z 0, V_ r 0 holds. This is thus a monotonic non-increasing function. Since V_ r 0, limt!1 V exists, i.e. V(1) exists. Thus sD is bounded, as is s. Since the continuous function of the closure set is bounded, then xi is bounded and so is s_D . As a consequence sD is uniformly continuous, as is V_ ¼RsD  s_D . Since V(t) is bounded and t limt!1 0 V_ dt ¼ Vð1Þ  Vð0Þ exists, then according to the Barbalat lemma (Slotine & Li, 1991) we have limt!1 V_ ¼ 0 and we obtain limt!1 s1 ¼ 0. To sum up the ADFSMC control law with ^ adjusted by (36), the method parameter vector s can find appropriate fuzzy rules in fuzzy control implementation and the online adjust rule also has the effect of improving the stability property. In particular, only one regulating factor is needed to tune the control rules. This will lead to less computational time and easier design than 234

Expert Systems, September 2006, Vol. 23, No. 4

conventional design approaches. The results show that the rule regulation of ADFSMC is stable in the sense of the Lyapunov function candidate. Compared to conventional design procedures of fuzzy sliding-mode controllers, the proposed controller design employs the expert knowledge of the sliding-mode to accomplish a control task that is easier and more efficient. Hence, the dead-zone concepts of nonlinear adaptive control are applied in this fuzzy rule regulation equation to improve the smoothness and robustness of the ADFSMC controller.

4. Computer simulations and experimental results In this section, we shall demonstrate that the ADFSMC is applicable to a cart–single-pendulum system, a cart–double-pendulum system, a ball–beam system (Lo & Kuo, 1998) and a tunnelling robot manipulator system to verify the theoretical development. 4.1. Case A: Cart–single-pendulum system The structure of a cart–single-pendulum system is illustrated in Figure 5. It consists of a straightline rail, a cart, a pendulum and a driving unit. The cart can move left or right on the rail. The pendulum is hinged at the centre of the cart and can rotate around the pivot in the same vertical



u

x

Figure 5: Structure of an inverted pendulum system.

 c 2006 The Authors. Journal Compilation  c 2006 Blackwell Publishing Ltd.

plane with the rail. Its dynamics are described below: x_ 1 ¼ x2 x_ 2 ¼

mt g sin x1  mp L sin x1 cos x1 x22 þ u cos x1 Lð43 mt  mp cos2 x1 Þ

simulation results. It is found that the pole and the cart can be stabilized to the equilibrium point. Further, the performance and robustness of the proposed control is better than those achieved by Lo and Kuo (1998) and Chen et al. (2002).

þd

x_ 4 ¼

4 3

ð45Þ

mp Lx22 sin x1 þ mp g sin x1 4 2 3 mt  mp cos x1

þ

3ð43

cos x1

4 uþd mt  mp cos2 x1 Þ

where x1 ¼ y is the angle of the pole with respect to the vertical axis; x2 ¼ y_ is the angular velocity of the pole with respect to the vertical axis; x3 ¼ x is the position of the cart; x4 ¼ x_ is the velocity of the cart; and mt ¼ mc þ mp. In what follows, we define the following variables: s1 ¼ c1 ðy  zÞ þ y_ ¼ c1 ðx1  zÞ þ x2 s2 ¼ c2 x þ x_ ¼ c2 x3 þ x4

4.2. Case B: Cart–double-pendulum system As shown in Figure 9, the cart–double-pendulum system considered here consists of a straight-line rail, a cart, a double pendulum and a driving unit. The cart can move left or right on the rail. The pendulum is hinged on the centre of the cart and can rotate around the pivot in the same vertical plane with the rail. Pole 1 is the pole

0.8 proposed control DFSM control SFSM control

0.6 0.4 0.2

ð46Þ ð47Þ

theta(rad)

x_ 3 ¼ x4

0 −0.2 −0.4 −0.6

and

−0.8

0 < Zu < 1

−1

ð48Þ

In the simulation, the following specifications are used: mp ¼ 0.05 kg, mc ¼ 1 kg, L ¼ 0.5 m, g ¼ 9.8 m=s2, c1 ¼ 5, c2 ¼ 0.5, F1 ¼ 5, Fz ¼ 15, Zu ¼ 0.9425, d ¼ 0.087sin(t), Z1 ¼ 2, Z2 ¼ 2. The initial values are y ¼  601, y_ ¼ 0 , x ¼ 0, x_ ¼ 0. In this simulation, it is impossible to regulate the cart position and pole angle simultaneously. If we want to regulate the cart position error, the pole will fall down. The sign of the pole angle is the same as that of the cart position; both the pole angle error and its velocity are zero, but the cart position is still in an undesired location; the cart position error and its velocity have the same sign. We produce the results for the ADFSMC assuming unknown nonlinear continuous functions f1 of the subsystem, with a DFSMC assuming known nonlinear continuous functions f1 of the subsystem. Figures 6–8 show the

−1.2

0

2

4

6

8

10 12 time(s)

14

16

18

20

Figure 6: Angular evolution of the pole.

1.2 proposed control DFSM control SFSM control

1 0.8 0.6 z

z ¼ satðs2 =Fz Þ Zu

0.4 0.2 0 −0.2

0

2

4

6

8

10 12 time(s)

14

16

18

20

Figure 7: System intermediate variable z of the cart–single-pendulum system.

 c 2006 The Authors. Journal Compilation  c 2006 Blackwell Publishing Ltd.

Expert Systems, September 2006, Vol. 23, No. 4

235

is the angle of pole 2 with respect to the vertical axis; x4 ¼ y_ 2 is the angular velocity of pole 2 with respect to the vertical axis; x5 ¼ x is the position of the cart; and x6 ¼ x_ is the velocity of the cart. In what follows, we define the following variables:

12 proposed control DFSM control SFSM control

10 8

x(m)

6 4 2 0

s1 ¼ c1 ðy  zÞ þ y_ ¼ c1 ðx1  zÞ þ x2

−2 −4

2

0

4

6

8

10 12 time(s)

14

16

18

20

s2 ¼ c2 x þ x_ ¼ c2 x3 þ x4

Figure 8: Position evolution of the cart.

sz ¼ cz x5 þ x6 z ¼ satðzp Þ Zu

2

 sz =Fz zp ¼ s2 =F2

1

ð51Þ

ð52Þ 0 < Zu < 1 if s2 r jsq j if s2 > jsq j

ð53Þ

where sq is the threshold value of s2. In the simulation, the following specifications are used: l1 ¼ 1 m, l2 ¼ 1 m, m1 ¼ 1 kg, m2 ¼ 1 kg, mc ¼ 1 kg, L ¼ 0.5 m, g ¼ 9.8 m=s2, c1 ¼ 5, c2 ¼ 1, cz ¼ 0.5, F1 ¼ 3, F2 ¼ 3, Fz ¼ 15, sq ¼ 0.01,

u

Zu ¼ 0.4712, d ¼ 0.08sin(t), Z1 ¼ 2, Z2 ¼ 2. The initial values are y1 ¼ 301, y2 ¼ 101, €y1 ¼ y_ 1 ¼ 0 ,

x

Figure 9: Structure of a cart–double-pendulum system. connected to the cart and pole 2 is above pole 1. The system’s dynamics are represented by x_ 1 ¼ x2 x_ 2 ¼ f1 þ b1 u þ d x_ 3 ¼ x4 x_ 4 ¼ f2 þ b2 u þ d

ð49Þ

x_ 6 ¼ f3 þ b3 u where x1 ¼ y1 is the angle of pole 1 with respect to the vertical axis; x2 ¼ y_ 1 is the angular velocity of pole 1 with respect to the vertical axis; x3 ¼ y2 Expert Systems, September 2006, Vol. 23, No. 4

€ y2 ¼ y_ 2 ¼ 0 , x€ ¼ x_ ¼ 0. We produce the results for the ADFSMC assuming unknown nonlinear continuous functions f1 of the subsystem, with a DFSMC assuming known nonlinear continuous functions f1 of the subsystem. Figures 10–12 show the simulation result. It is found that the pole and the cart can be stabilized to the equilibrium point. Further, the performance and robustness of the proposed control is better than those achieved by Lo and Kuo (1998) and Lin and Mon (2005). 4.3. Case C: Ball–beam system

x_ 5 ¼ x6

236

ð50Þ

The ball and beam system shown in Figure 13 consists of a beam pivoting at the centre point and a ball free to move along the beam. Let r be the distance of the ball from the origin and y be the angle that the beam makes with the horizontal line. The objective is to design an ADFSMC where the output converges to zero from

 c 2006 The Authors. Journal Compilation  c 2006 Blackwell Publishing Ltd.

arbitrary initial conditions in a certain region. Its dynamics are described below:

20 proposed control DFSM control SFSM control

15

x_ 1 ¼ x2

10

x(m)

x_ 2 ¼ u þ d ð54Þ

x_ 3 ¼ x4

5 0

x_ 4 ¼ Bðx3 x22  G sin x1 Þ −5

0

2

4

6

8

10 12 time(s)

14

16

18

20

Figure 12: Position evolution of the cart. 10 proposed control DFSM control SFSM control

8 6 r(m)

where x1 ¼ y is the angle of the pole with respect to the vertical axis; x2 ¼ y_ is the angular velocity of the pole with respect to the vertical axis; x3 ¼ r is the position of the cart; x4 ¼ r_ is the velocity of the cart; B ¼ MR2=(Jb þ MR2); Jb is the moment of inertia of the ball; M is the mass of the ball; R is the radius of the ball; and g is the acceleration due to gravity. The centre of rotation is assumed to be frictionless and the ball is free to roll along the

4 2

0.6

0

proposed control DFSM control SFSM control

0.5 0.4

−2

theta1(rad)

0.3

0

2

4

6

8

10 12 time(s)

14

16

18

20

Figure 13: Structure of a ball–beam system.

0.2 0.1 0

−0.1 −0.2 −0.3

0

2

4

6

8

10 12 time(s)

14

16

18

20

Figure 10: Angular evolution of pole 1. 1 proposed control DFSM control SFSM control

0.8 0.6 theta2(rad)

0.4 0.2 0

−0.2 −0.4 −0.6 −0.8

0

2

4

6

8

10 12 time(s)

14

16

Figure 11: Angular evolution of pole 2.

18

20

beam. It is required that the ball remains in contact with the beam and that rolling occurs without slipping. The objective is to keep the ball close to the centre of the beam close to the horizontal position. In the simulation, the following specifications are used: B ¼ 0.7143, Jb ¼ 2  10 6, M ¼ 0.05 kg, R ¼ 0.01 m, g ¼ 9.8 m=s2, c1 ¼ 5, c2 ¼ 0.5, F1 ¼ 5, Fz ¼ 15, Zu ¼ 0.9425, d ¼ 0.08 sin(t), Z1 ¼ 2, Z2 ¼ 2. The initial values are x1 ¼ y ¼ 601, x2 ¼ y_ ¼ 00 , x3 ¼ 10, x4 ¼ r_ ¼ 0. We produce the results for the ADFSMC assuming unknown nonlinear continuous functions f1 of the subsystem, with a DFSMC assuming known nonlinear continuous functions f1 of the subsystem. Figures 14–16 show the simulation result. It is found that the ball– beam system can be stabilized to the equilibrium point and that y and r both converge to zero.

 c 2006 The Authors. Journal Compilation  c 2006 Blackwell Publishing Ltd.

Expert Systems, September 2006, Vol. 23, No. 4

237

Further, the proposed control performance and robustness are better than those achieved by Lo and Kuo (1998) and Chen et al. (2002).

1 proposed control DFSM control SFSM control

0.8 0.6 0.4 z

4.4. Case D: Tunnelling robot manipulator system

0.2 0 −0.2 −0.4

0

2

4

6

8

10 12 time(s)

14

16

18

20

Figure 15: System intermediate variable z of the ball–beam system. 1.2 proposed control DFSM control SFSM control

1 0.8 theta(rad)

Robot manipulations in a tunnel for civil engineering (Figure 17) form the hardware of the system (Chou et al., 2003). The hardware of the robot arm (Figure 18) contains (1) the forward and back motion of the x axis, range 0–29.4 cm; (2) the swirl motion of the y axis, range þ 781 to –781; (3) the sliding motion of the R axis, range 0–20 cm. Each axis includes both a limited switch and a home-tested sensor switch. The home-tested sensor allows users to return to the initial position, write programs and measure exactly. Both the limited switch and the hometested sensor switch work by magnetic induction. A supersonic sensor is hung on the R axis to detect the surface of the tunnel. The wooden semicircular tunnel is used to simulate an unsmooth real tunnel for three-axis robot manipulation. There is a 10 cm  10 cm  3.2 cm concavity for the experiment. The aspects of the interfaces are as follows: (1) two encoder=DAC interfaces; (2) an analog input–output interface; (3) a digital input–output interface. The utilization of a D=A card makes u become an analog signal. The position control of the y axis and x axis means that the motion of the robot arm is from its initial position to the set point. A comparison of the proposed control with the fuzzy sliding-

0.6 0.4 0.2 0 − 0.2 − 0.4

0

2

4

6

8

10 12 time(s)

14

16

18

20

Figure 16: Angular evolution of the beam.

10 proposed control DFSM control SFSM control

8

Figure 17: A real robot manipulator for civil engineering in a tunnel.

r(m)

6 4 2 0 −2

0

2

4

6

8

10 12 time(s)

14

16

Figure 14: Position evolution of the ball. 238

Expert Systems, September 2006, Vol. 23, No. 4

18

20

mode control (Chou et al., 2003) for the same plant is shown in Figures 19 and 20, respectively. It is clear that the proposed method not only reduces overshoot to near zero but also maintains a small setting time and steady error. Experimental results showed that the proposed  c 2006 The Authors. Journal Compilation  c 2006 Blackwell Publishing Ltd.

20

proposed control FSM control

Translational motion R axis

x (cm)

15 Rotational motion  axis

10 5 0

Forward and Backward motion x axis

theta (degree)

proposed control FSM control

0

2

4

6

8

10 12 time (s)

14

16

18

2

4

6

8

10 12 time (s)

14

16

18

20

Figure 20: The proposed controller for the x axis is moved from 0 cm to 15 cm.

Figure 18: Diagram of the three-axis robot.

50 45 40 35 30 25 20 15 10 5 0

0

20

Figure 19: The response of the position controller for the y axis is moved from 01 to 401. control determines the robot position accurately and constructs tunnels precisely.

5. Conclusions An ADFSMC design scheme is presented through width adaptation for fourth-order nonlinear systems. The proposed adaptation law, which results from the direct adaptive approach, is used to appropriately determine the width of the unknown system variables. The membership functions in the THEN part vary with the width adaptation. An adaptive law is then used to tune width in the THEN part to appropriately determine the distribution of each membership function. The fuzzy rules of DFSMC should be preconstructed by trial-and-error tuning. By the ADFSMC, the fuzzy rules can be learned online by an adaptive law and the stability of the proposed ADFSMC system can be guaranteed. This control strategy establishes the appropriate fuzzy rules by continuous online learning in-

stead of a trial-and-error process. It simplifies the implementation of a fuzzy controller. The response of the system converges faster than results achieved in previous work. Simulation results have been provided to demonstrate the robustness and effectiveness of the proposed control system. Moreover, the experimental results have demonstrated that the proposed control scheme is effective and valid for real-world engineering applications.

References CHAI, T. and S. TONG (1999) Fuzzy direct adaptive control for a class of nonlinear systems, Fuzzy Sets and Systems, 103 (3), 379–387. CHEN, J.Y. (2001) Rule regulation of fuzzy slidingmode controller design: direct adaptive approach, Fuzzy Sets and Systems, 120 (1), 159–168. CHEN, S.Y., F.M. YU and H.Y. CHUNG (2002) Decoupled fuzzy controller design with single-input fuzzy logic, Fuzzy Sets and Systems, 129 (3), 335–342. CHOU, W.H., C.N. HUNG and H.Y. CHUNG (2003) A fuzzy sliding controller with gray predictor for the motion of a robot manipulator, Proceedings of the 25th National Symposium on Automatic Control, 805–810. FILEV, D.P. and R.R. YAGER (1994) On the analysis of fuzzy logic controllers, Fuzzy Sets and Systems, 68 (2), 39–66. GLOWER, J.S. and J. MUNIGHAN (1997) Design fuzzy controllers from a variable structures standpoint, IEEE Transactions on Fuzzy Systems, 5 (1), 138–144. HWANG, G.C. and S.C. LIN (1992) A stability approach fuzzy control design for nonlinear systems, Fuzzy Sets and Systems, 48 (3), 279–287.

 c 2006 The Authors. Journal Compilation  c 2006 Blackwell Publishing Ltd.

Expert Systems, September 2006, Vol. 23, No. 4

239

KAZEMIAN, H.B. (2001) Study of learning fuzzy controllers, Expert Systems, 18 (4), 186–193. KIM, S.W. and J.J. LEE (1995) Design of a fuzzy controller with fuzzy sliding surface, Fuzzy Sets and Systems, 71 (2), 359–367. LAU, H. and T.N. WONG (1998) A fuzzy expert system for complex closed-loop control: a non-mathematical approach, Expert Systems, 15 (2), 98–109. LEE, C.C. (1990) Fuzzy logic in control systems: fuzzy logic controller – Part I, IEEE Transactions on Systems, Man and Cybernetics, 20, 404–418. LEE, M.R. (2002) Expert system for nuclear power plant accident diagnosis using a fuzzy inference method, Expert Systems, 19 (4), 201–207. LEUNG, R.W.K., H.C.W. LAU and C.K. KWONG (2003) On a responsive replenishment system: a fuzzy logic approach, Expert Systems, 20 (1), 20–32. LIANG, C.Y. and J.P. SU (2005) A new approach to the design of a fuzzy sliding-mode controller, Fuzzy Sets and Systems, 139 (1), 111–124. LIM, C.M. and T. HIYAMA (1991) Application of fuzzy logic control to a manipulator, IEEE Transactions on Robotics and Automation, 7 (5), 688–691. LIN, C.M. and Y.J. MON (2005) Decoupling control by hierarchical fuzzy sliding-mode controller, IEEE Transactions on Control Systems Technology, 13 (4), 593–598. LO, J.C. and Y.H. KUO (1998) Decoupled fuzzy sliding-mode control, IEEE Transactions on Fuzzy Systems, 6 (3), 426–435. PALM, R. (1994) Robust control by fuzzy sliding-mode, Automatic, 30 (9), 1429–1437. PHAM, T.T. and G. CHEN (2002) Some applications of fuzzy logic in rule-based expert systems, Expert Systems, 19 (4), 208–223. SLOTINE, J.J.E. and W. LI (1991) Applied Nonlinear Control, Englewood Cliffs, NJ: Prentice Hall. WANG, L.X. (1994) Adaptive Fuzzy Systems and Control Design and Stability Analysis, Englewood Cliffs, NJ: Prentice Hall, 140–154.

240

Expert Systems, September 2006, Vol. 23, No. 4

WU, W.Z., W.X. ZHANG and H.Z. LI (2003) Knowledge acquisition in incomplete fuzzy information systems via the rough set approach, Expert Systems, 20 (5), 280–286.

The authors Lon-Chen Hung Lon-Chen Hung was born in Kaoshiung, Taiwan, Republic of China. He is currently pursuing a PhD degree in electrical engineering from the National Central University, Taiwan. His research interests include sliding-mode control, fuzzy theory and control, and neural network applications.

Hung-Yuan Chung Hung-Yuan Chung was born in Ping-Tung, Taiwan, Republic of China. He received a PhD in electrical engineering from the National Cheng Kung University, Tainan, Taiwan, in 1987. In August 1987 he joined the Department of Electrical Engineering at the National Central University, Chung-li, Taiwan, as an associate professor. He has been a professor since August 1992. His research and teaching interests include system theory and control, adaptive control, fuzzy control, neural network applications and microcomputer-based control applications.

 c 2006 The Authors. Journal Compilation  c 2006 Blackwell Publishing Ltd.