Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
3703.pdf
Скачиваний:
9
Добавлен:
15.11.2022
Размер:
11.37 Mб
Скачать

156

Danial Senejohnny, Mohammadreza Faieghi and Hadi Delavari

 

 

The stability is guaranteed if the derivative of the Lyapunov function is negative definite, also known as the reaching condition. Taking the first time derivative from (18) one gets:

 

 

 

V S

S

k

(1

sb k D )

k D S

f

k

(x) k Db k s

S

k

 

 

 

(19)

 

 

 

 

k

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

k k k

k k

 

 

 

 

 

 

 

 

 

 

k k k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

 

(x)

 

 

Assuming that

f

k

(x)

are bounded, namely

k

 

sup

x

k

, and s 1 k Db , then

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

k

k

one gets 1 s k Db

0 . As a result, the expression (19) can be rewritten as

 

 

k

k k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Sk

 

 

 

1 ksbk kkD kkD k

 

Sk

 

 

kkDbk kks

 

 

 

 

 

 

Vk

Sk

 

 

 

 

 

Sk

 

 

(20)

 

 

 

 

 

 

Sk

 

 

 

Sk

 

 

1 ksbk kkD kkDbk kks kkD k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

Sk

 

ks kkDbk 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Therefore,

provided

k s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

to

dominate the model

 

matched

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

bk

 

 

 

 

 

kkD

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

uncertainties [33], then Vk will be negative which implies asymptotic stability of the system and Sk (t) and Sk (t) tend to zero as t .

4. Adaptation Law Synthesis

In this section an adaptive scheme based on Lyapunov function for switching feedback control gains is derived. The adaptation law guarantees the system stability regardless of model uncertainty bound. Take the following Lyapunov function

 

 

 

 

 

Vk

 

1

Sk

Sk

2k 1

(k

)

2

2k

(kk )

2

 

(21)

 

 

 

 

 

 

 

 

2

2

1

s

 

1

s

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

s

ˆs

s

 

s

ˆs

 

s

 

 

 

 

 

 

 

 

 

 

 

where k

k

k

,

kk

kk

kk

denote the estimation error. Taking time derivative from

(21) with respect to time and with the sense of (20), one has:

Vk

 

Sk

 

 

Sk

 

s

D

D

 

s

D

 

Sk

 

1

 

s

s

1 s s

(22)

 

 

 

 

 

 

 

 

 

1 k bk kk

kk

bk kk

kk k

2k 1 k

k

2k kk kk

adding and subtracting some terms, Eq. (22) can be rewritten as follows

 

 

 

 

 

 

 

Vk

 

Sk

 

 

 

Sk

 

 

ˆs D

 

ˆs

D

 

D

 

s

D

 

Sk

 

 

s

D

 

Sk

 

Sk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 k kk bk

kk bk kk

kk k

kk bk kk

 

 

k bk kk

 

 

 

(23)

1

s s

1k s k s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2k 1

 

 

k k

 

 

 

2k k

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆs

 

 

 

 

ˆs

can be chosen such that

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

It is clear that kk

and k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Complimentary Contributor Copy

Adaptive Second-Order Fractional Sliding Mode Control …

157

 

 

 

 

 

 

 

 

 

 

Sk

Sk 1

k bk kk

kk bk kk

kk k 0 ,

 

 

 

 

 

 

ˆs D

ˆs D

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

so the adaptation law can be derived directly from (23) as follow

sb k D

 

S

k

 

S

k

1

s s 0 s

2k 1

k Db

 

S

k

 

S

k

 

 

 

 

 

 

k k k

 

 

 

 

2k 1

k k

k

 

 

k k

 

 

 

 

 

 

(24)

 

 

 

 

1k s k s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k sb k D

 

S

k

0

k s

 

2k

k Db

S

k

 

 

 

 

 

 

k k k

 

 

 

 

 

2k k k

 

k

 

k k

 

 

 

 

 

 

 

The switching controller parameters ks and kks in switching control (13) are evaluated

through the adaptive algorithm derived in (24).

As the bound of uncertainty is not available, choosing switching controller gains manually can make chattering effect worst. Therefore, by virtue of adaptation law (24) switching controller parameters assume their gains by making a compromise between chattering suppression and reaching condition.

5. Numerical Studies

The twin-tanks system comprised of two small tanks mounted above a reservoir which provides storage for the water, Figure 2. Water is pumped into the bottom of each tank by two independent pumps. The pump only increases the liquid level and is not responsible for pumping the water out of the tank [8].

Figure 2. The plant of a liquid-level control system.

The twin-connected tanks system is a nonlinear dynamical system and the governing dynamical equations can be written as [8, 34]

Complimentary Contributor Copy

158

Danial Senejohnny, Mohammadreza Faieghi and Hadi Delavari

 

 

 

k1

 

 

 

 

 

sign(h h )

1

 

 

 

 

 

h

 

 

h h

q

 

 

 

 

 

 

 

 

1

 

 

A1

 

 

1

2

 

1

2

 

A1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(25)

 

 

k1

 

 

 

 

 

 

 

 

 

 

k2

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h2

 

 

h1 h2

sign(h1

h2 )

h2

q2

A2

 

A2

A2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

where h1 and h2 are the total water heads in Tank 1 and Tank 2 respectively, which are the two outputs of interest, and q1 , q2 are the two inflows into the tanks. It is assumed that A1 and A2 are the capacities of Tank-1 and Tank-2 respectively. Both tanks assumed to have the same cross sectional area i.e. A1 A2 A .

It has been seen in [35] that in chemical plants, selecting the flow rate as an input is more effective than using flow as the input. Thus, if the flow rates are considered as the inputs, i.e., q1 u1 and q2 u2 , the system dynamics can be written as

h1

k1

 

 

h1 h2

 

 

 

 

1

u1

 

 

 

 

2A

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

h1 h2

 

 

 

 

 

(26)

 

k1

 

 

 

h1 h2

 

 

 

 

 

 

k2

 

 

h2

 

 

 

 

 

h2

1

u2

2A

 

 

 

 

 

 

 

 

 

A

 

 

h1 h2

 

 

 

2A h2

It seems from (25) and (26) that there is a discontinuity at h1 h2 , so desired water level is selected as hd1 hd 2 .

A mathematical model of Twin Tank can be expressed as follows:

x1 (t) x2 (t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k1

 

x2 x4

1

 

 

 

 

 

x2 (t)

 

 

 

 

 

 

 

 

 

 

 

 

 

u1

 

 

 

 

2 A

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

x1 x3

 

 

 

 

(27)

x3 (t) x4 (t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k1

x2 x4

 

k2

 

 

1

 

x4 (t)

 

 

 

 

 

 

 

 

 

x4

 

 

u2

2A

 

 

x1 x3

 

 

2A x3

A

 

 

 

 

System constraints: The flowed fluid into

the tanks

( q1 and q2 ) cannot be negative

because the pumps can only pump water into the tanks. Therefore constraints on the inflow

are given by

q1 , q2 0 . In the steady state, for constant water level set points, the respective

derivatives

must be zero separately i.e. h1 0,

h2 0 . Therefore, the inequality

 

 

k 2

 

h

1 must hold [8].

 

 

 

1

2

 

 

k 2

k 2

h

 

 

 

 

 

1

2

 

1

 

 

Complimentary Contributor Copy

 

Adaptive Second-Order Fractional Sliding Mode Control …

 

 

 

159

 

 

 

 

 

 

Both tanks have the same area of 70 cm2 ,

k

14 and

k

2

10 , the set point water level is

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

selected as

hd ,1 xd ,1

8 cm and

hd ,3 xd ,3

6 cm , with

initial

water

level

h1 (0) 4 ,

h2 (0) 2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

In Eq.

(9) sliding

surface

parameters

are

selected

to

be as

k P 3 ,

k D

7

, k FD 1,

 

 

 

 

 

 

 

 

 

 

 

 

1

1

 

 

1

k I 0.5 ,

1,

0.8 ,

k P

3 , k D

7 , k FD 1,

 

2

1, k I 0.5

,

2

0.8 . The

1

1

1

 

2

2

 

 

2

 

2

 

 

 

initial conditions are chosen as

S1 (0) 0.5

and

S2 (0) 0.5

as well. In order to smooth the

control action, instead of sgn(.) function tanh(.) is used.

 

 

 

 

 

 

 

 

The simulation

results with

20% variations in

system’s

nominal

parameters (i.e

A 84 cm2 , k1 16.8 and k2 12 ) shows good stabilization of the water levels in both tanks

Figs. 3, dashed-line and solid-line represents system’s response and desired trajectory, respectively. The value of inflow in both tanks settles to the steady state values shown in Figure 3, the pumps only increases the liquid level, so their value is always positive.

h1(m)

inlet flow, q1

inlet flow, q2

9

 

 

 

 

7

 

 

 

8

 

 

 

 

6

 

 

 

7

 

 

 

 

5

 

 

 

6

 

 

 

(m)

4

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

h

 

 

 

 

5

 

 

 

 

3

 

 

 

4

 

 

 

 

2

 

 

 

30

10

20

30

 

10

10

20

30

 

time (sec)

 

 

 

 

time (sec)

 

300

 

 

 

 

 

 

 

 

200

 

 

 

 

 

 

 

 

100

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

-1000

5

10

15

 

20

25

30

35

 

 

 

 

time (sec)

 

 

 

 

300

 

 

 

 

 

 

 

 

200

100

0

-1000

2

10

15

20

25

30

35

 

 

 

 

time (sec)

 

 

 

Figure 3. 2-SMC with PIDDα surface and +20% variation in parameters of Twin Tank system.

Complimentary Contributor Copy

160

Danial Senejohnny, Mohammadreza Faieghi and Hadi Delavari

 

 

8

1

8

1

k

λ

8

2

8

2

k

λ

Figure 4. Adaptive parameters versus time.

1

 

2

 

 

 

 

 

 

 

S1

S2

Figure 5. Phase plot.

Adaptive parameters 1s , k 1s , 2s , k 2s in (24) are estimated as follow with zero initial conditions shown in Figure 4.

s k Db

 

S

 

 

S

, k s k Db

S

, s

k Db

 

S

2

 

S

2

k s

k Db

S

2

 

 

 

 

1

1

1

1

 

1

 

 

1

1

2

1

1

1

2

3

2

2

 

 

 

 

2

4

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Complimentary Contributor Copy

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]