# Adaptive Fixed-Time Fast Terminal Sliding Mode Control for Chaotic Oscillation in Power System

Authors: Caoyuan Ma; Faxin Wang; Zhijie Li; Jianyu Wang; Chuangzhen Liu; Wenbei Wu; Yuzhou Cheng
Journal: Mathematical Problems in Engineering
 (2018)

---

## Abstract

The second-order chaotic oscillation system model is used to analyze the dynamic behavior of chaotic oscillations in power system. To suppress chaos and stabilize voltage within bounded time independent of initial condition, an adaptive fixed-time fast terminal sliding mode chaos control strategy is proposed. Compared with the conventional fast terminal sliding mode control strategy and finite-time control strategy, the proposed scheme has advantages in terms of convergence time and maximum deviation. Finally, simulation results are given to demonstrate the effectiveness of the proposed control scheme and the superior performance.


---

## Body

## 1. Introduction

As a typical multivariable and strongly coupled nonlinear system, power system exhibits a lot of nonlinear dynamic behaviors during operation, such as low-frequency oscillation, bifurcation, and chaos [1]. If the amplitude of the disturbance satisfies certain conditions, the power system will appear with continuous and random chaotic oscillations which may lead to instability of the system, voltage collapse, and even catastrophic blackouts [2, 3]. Therefore, it is imperative to study control schemes for chaos suppression in power system. In 1993, Chiang et al. [4] observed the chaotic behavior under different load conditions in the power system firstly. Since then, chaos control had become a hot issue of academia. The control methods of chaotic oscillation can be divided into two aspects: guiding the development of the chaotic oscillation to the expected orbit and suppressing the occurrence of system chaos [5, 6]. At present, the commonly used control methods include parameter perturbation method (OGY) [7–10], feedback control [11–14], adaptive control [15–17], fuzzy control [18–20], and backstepping control [21]. These approaches have important theoretical and practical significance to guarantee the stability and secure operation of power system while some defects did exist at the same time. For example, the backstepping control was proposed to systematize and structure the design process of the system's Lyapunov functions and controllers through the reverse design. But its structure is very complex, and the complexity of regression matrix would become stronger especially when the nonlinear damping existed for the system parameters uncertainty [22, 23]. Fuzzy control is difficult to adapt to the requirements of large-scale adjustment. It needs to constantly adjust the control rules and parameters [24–27]. Moreover, all the above control methods can only reach the asymptotic stability; that is, the convergence time cannot be allocated in advance. From the point of view of the operation of the power system, the oscillation is acceptable if it is damped within a limited time.Fast terminal sliding mode (FTSM) control can achieve system stability in finite time. It has the advantages of strong robustness to external disturbances and parameter disturbances [28, 29]. The fast terminal sliding mode (FTSM) control can converge the system to the equilibrium point on the sliding mode surface in finite time while the traditional sliding mode control can only ensure that it reaches the sliding mode surface during the same period [30]. Finite-time control can also keep the system stable within bounded time [31]. Because the finite-time control has stronger robustness, better anti-interface, and faster convergence, it has become a research hotspot of academia in recent years. Cai et al. [32] applied the finite-time theory for the first time to deal with the generalized synchronization problem of different order chaotic systems. Khaki et al. [33] proposed a novel fuzzy finite-time variable structure controller for double integrator power systems to damp the complicated chaotic oscillations of an interconnected power system, when such oscillations can be made by load perturbation of a power system working on its stability edges. The above two control methods can make the system stable in finite time but also subject to the limitless of finite-time stability. The system stability is affected by the initial condition while it is difficult to obtain accurate parameters in practical project.Fixed-time stability is an extension of finite-time stability. Compared with the above control approaches, fixed-time control can not only maintain stronger robustness and better anti-interference and ensure a clear upper bound of settling time, but also make the system globally uniformly ultimately bounded stable one with any initial conditions [34]. Polyakov firstly proposed the concept of fixed-time in 2012 [35]. At present, due to this attractive property, fixed-time control has begun to be widely used in hydraulic turbine governing systems [36], PMSM [37], multiagent systems [38], aircraft systems [39], chaotic systems [40], and so on. Ni et al. [40] firstly applied the STATCOM to suppress chaos of power system and authors designed a fixed-time dynamic surface high-order sliding mode control approach to achieve semiglobally fixed-timely uniformly ultimately bounded stabilization. In view of the discussion of the above control schemes, this paper presents an adaptive fixed-time fast terminal sliding mode control scheme which can accelerate the convergence, suppress chaos in the power system, and avoid voltage collapse. The main advantages of proposed controller are that it can guarantee the system stability in finite time independent of initial state and the settling time can be directly calculated.The structure of this paper is as follows: the second section gives some theorems and lemmas to facilitate the derivation of the system control method in the later text. In the third section, an adaptive fixed-time fast terminal sliding mode control scheme is proposed to suppress chaos in the power system. The fourth section provides the simulation results in this paper to demonstrate the effectiveness of the proposed control scheme. The fifth section has made some summary conclusion based on the above work.


## 2. Preliminaries

In this paper, for the convenience of analysis, we firstly introduce a necessary definition and some lemmas which play an important role in design process.Consider the following differential equation system:(1)x˙t=fxt,x0=x0where x∈Rn is the system state variable;f: R+×Rn→Rn is a smooth nonlinear function. Assume the origin is an equilibrium point of (1).Definition 1 (see [41–43]).
The origin of system (1) is a finite-time stable equilibrium if the origin is Lyapunov stable and there exists a functionT: Rn→R+, called the settling time function, such that, for every x0∈Rn, the solution x(t,x0) of system (1) is defined on [0,Tx0),x(t,x0)∈Rn, for all t∈[0,T(x0)), and limt→T(x0)⁡x(t,x0)=0.Definition 2 (see [35, 43]).
The origin of system (1) is said to be fixed-time stable equilibrium point if it is globally finite-time stable with bounded convergence time T(x0); that is, there exists a bounded positive constant Tmax such that T(x0)<Tmax satisfies.Lemma 3 (see [35, 42]).
If there exists a continuous functionV: Rn→R+U|0| such that(1)
V(x)=0⇔x=0,(2)
V(x) is radially unbounded,(3)
for someα, β,p,q,k>0, andpk<1
any solutionx(t) satisfied the inequality (2)D∗Vxt≤-αVpxt+βVqxtkwhere D∗V(x(t)) is the upper right-hand derivative of the function V(x(t)), then the origin is globally fixed-time stable and the following estimate holds:(3)Tx0≤1αk1-pk+1βkqk-1,∀x0∈RnLemma 4 (see [44]).
For any nonnegative real numbers, that is,x1,x2,  …,xN≥0, the following inequality holds:(4)∑i=1Nxiη≥∑i=1Nxiη,0<η≤1∑i=1Nxiς≥N1-ς∑i=1Nxiς,ς>1


## 3. Main Results

Ignoring the dynamic process of excitation loop and damping winding, it is assumed that the mechanical power of generator is always the same in transient process, and the transient salient effect of generator is not considered. The external factors of system, such as disturbance’s influence on the system, are mainly considered. The equivalent circuit of power grid second-order chaotic oscillation system is presented in Figure1, whereG1,G2 andT1,T2 are equivalent generators and main transformers of system, respectively.Figure 1
Equivalent circuit of power grid two-order chaotic oscillation system.The power system is composed of two generator buses and one load bus. Here, the second-order nonlinear mathematical model of the synchronous generators used in this paper is as follows:(5)δ˙=ωω˙=-1HPssin⁡δ-DHω+1HPm+1HPecos⁡λtwhere δ and ω are state variables that denote relative angle and angular frequency between equivalent generators, respectively;H andD refer to generator inertia and damping coefficient, respectively;Ps andPm are electromagnetic power and mechanical power of generator, respectively; Pe and λ represent amplitude and frequency of load disturbance. To simplify the model, we make the following transformation: (6)x˙1=x2x˙2=-asin⁡x1-bx2+c+Fcos⁡λt+u=fx+Fcos⁡λt+uwhere a=Ps/H, b=D/H, c=Pm/H, F=Pe/H, [x1,x2]=[δ,ω], fx=-asin⁡x1-bx2+c, andu is control input.Then the fast terminal sliding mode surface can be designed as(7)s=x2+α0x1+β0x1q0/p0where parameters satisfy α0, β0 >0, andq0,p0 are positive odd numbers. The derivative of the system (7) gives(8)s˙=x˙2+α0x˙1+q0p0β0x1q0/p0-1x2=fx+Fcos⁡λt+u+α0x2+q0p0β0x1q0/p0-1x2For system (6), we can see thatX0 = (z, 0) is a point of equilibrium of the system, wherez is a constant. To realize the global stability of system (6) based on the fixed-time stability theory in the case of the equilibriumX0, another global fast terminal sliding mode is taken as the following form: (9)s˙+ϕs+γsq/p+ksgnssα+ksgnssβ=0where parameters satisfy ϕ, γ>0, 0<α<1, β>1, andq<p.p andq are positive odd numbers. According to the above analysis, we propose our main results in the following theorem.Theorem 5.
System (6) can be stabilized in fixed-time under the following adaptive controller:(10)u=-fx-F^cos⁡λt-α0x2-q0p0β0x1q0/p0-1x2-ϕs-γsq/p-ksgnssα-ksgnssβwhere F^ is the estimated value of the uncertain parameter andk is the tuning parameter of the terminal attractor satisfying the following adaptive laws: (11)F^˙=scos⁡λt-F^-Fα-F^-Fβk˙=sα+1+sβ+1-k-gα-k-gβwhere g is the arbitrary constant.Proof.
Firstly, substitute (10) into (8) to get the following formula which will be used directly in the process of proof: (12)s˙=Fcos⁡λt-F^cos⁡λt-ϕs-γsq/p-ksgnssα-ksgnssβ
Consider the following Lyapunov candidate function as(13)V=12s2+12F^-F2+12k-g2
By using the designed controlleru and the appropriate tuning parameters, the time derivative ofV can be obtained as(14)V˙t=ss˙+F^-FF^˙+k-gk˙=sFcos⁡λt-F^cos⁡λt-ϕs-γsq/p-ksgnssα-ksgnssβ+F^-Fscos⁡λt-F^-Fα-F^-Fβ+k-gsα+1+sβ+1-k-gα-k-gβ=-ϕs2-γsp+q/p+sF-F^cos⁡λt+sF^-Fcos⁡λt-F^-Fα+1-F^-Fβ+1-ksα+1-ksβ+1+k-gsα+1+k-gsβ+1-k-gα+1-k-gβ+1=-ϕs2-γsp+q/p-gsα+1-gsβ+1-F^-Fα+1-F^-Fβ+1-k-gα+1-k-gβ+1
Here, owing to the fact thatp andq are positive odd numbers, the process of derivation continues:(15)
V
˙
t
≤
-
g
s
α
+
1
-
g
s
β
+
1
-
F
^
-
F
α
+
1
-
F
^
-
F
β
+
1
-
k
-
g
α
+
1
-
k
-
g
β
+
1
=
-
2
1
/
2
α
+
1
g
1
2
s
2
1
/
2
α
+
1
-
2
1
/
2
α
+
1
1
2
F
^
-
F
2
1
/
2
α
+
1
-
2
1
/
2
α
+
1
1
2
k
-
g
2
1
/
2
α
+
1
-
2
1
/
2
β
+
1
g
1
2
s
2
1
/
2
β
+
1
-
2
1
/
2
β
+
1
1
2
F
^
-
F
2
1
/
2
β
+
1
-
2
1
/
2
β
+
1
1
2
k
-
g
2
1
/
2
β
+
1
≤
-
m
1
2
s
2
1
/
2
α
+
1
+
1
2
F
^
-
F
2
1
/
2
α
+
1
+
1
2
k
-
g
2
1
/
2
α
+
1
-
n
1
2
s
2
1
/
2
β
+
1
+
1
2
F
^
-
F
2
1
/
2
β
+
1
+
1
2
k
-
g
2
1
/
2
β
+
1where m=min21/2(α+1)g,21/2(α+1),n=min21/2(β+1)g,21/2(β+1). Thus, it follows from Lemma 4 that(16)V˙t≤-m12s2+12F^-F2+12k-g21/2α+1-31-β/2n12s2+12F^-F2+12k-g21/2β+1=-mV1/2α+1-31-β/2nV1/2β+1
According to Lemma3, system (6) reaches the equilibrium within a bounded time and the bound of convergence time can be estimated by(17)T≤1m1-α+1/2+131-β/2nβ+1/2-1=2m1-α+231-β/2nβ-1which means that the state variables of the system can stabilize when t≥T.


## 4. Numerical Simulations

In this section, numerical simulations are performed to demonstrate the effectiveness and the superiority of the proposed control method. System (6) is selected for simulation. Without loss of generality, the system parameters area=1,b=0.02, andc =0.2 and controller parameters are α0=1, β0 =2,q0 =5,p0 =9, ϕ=20, γ=15, α=0.5, β=1.5,p=3,q=1, g=0.3, andk(0)=0.2. The initial values of the state variable areS(δ0, ω0) = (0.43,0.003). In the power system, the damping coefficient and the inertial coefficient of the generator are often constant, and the variation of the load disturbance can usually cause the chaotic oscillation of the power system. Further, when amplitude of periodic load disturbanceF varies to 0.2593, the power system exhibits chaotic oscillations which are displayed in Figures 2 and 3. As the parameters are selected, Lyapunov exponents can be calculated as (λ1,λ2,λ3) =(0.0174, 0, −0.0374). There is one positive Lyapunov exponent, which validates the existence of chaotic attractor. Figure 2 is the phase portrait of chaotic power system and it shows chaotic behavior of the power system clearly. Figure 3 shows the time responses of state variables in chaotic power system. As can be seen from Figure 3, the time responses of state variables are in an irregular and aperiodic oscillatory state and their trajectories are unpredictable after a long period of time.Figure 2
Phase portrait of chaotic power system.Figure 3
Time responses of state variables in chaotic power system.The chaotic oscillation state has great harm to the power system. Circuit voltage and current waveform distortion caused by the oscillation, particularly over voltage, will cause severe local instability with potentially serious impact and damage to the power system. Therefore, an immediate control action needs to be activated to suppress chaos. Simulations are conducted to examine the proposed controller performance in terms of suppressing chaos in power system and stabilizing the power system to its desired operating point. The time responses of state variables and the phase portrait of chaotic power system under proposed controller are presented in Figures4 and 5, respectively. Figures 6 and 7 are the time response of tuning parameter of the terminal attractor and the time response of the proposed controlleru. As Figure 7 shows, the controller was added at 50s. Further, the bound of convergence timeT can be calculated by taking the controller parameters we selected into (17). After calculation, we getT=14.59s. Therefore, the theoretical limit of stabilization time in the numerical simulations ist1=64.59s. Figure 4 shows that the control objectives have been controlled to the stability value. The system state variables δ and ω reached stable states at 59.2s and 53.3s, respectively, and they are all smaller thant1. Figure 5 depicts that the control objective has been stabilized to its desired operating point and the chaotic oscillation has been suppressed completely, which verifies the effectiveness of proposed controller.Figure 4
Time responses of state variables in chaotic power system with proposed controller.Figure 5
Phase portrait of chaotic power system with proposed controller.Figure 6
Tuning parameter of terminal attractork.Figure 7
Time response of the proposed controlleru.To explore the influence of different initial conditions on the control effect, the simulation has been done. Figure8 shows the time responses of state variables when the initial states of system areS1=(0.43,0.003),S2=(0.45,0.005), andS3=(0.47,0.007). As can be seen from Figure 8, the oscillation amplitude and overshoot are different, but the difference between the settling times is very small. The system state variable δ achieves stable state at 59.2s, 59.4s, and 59.7s, respectively, and the system state variableωachieves stable state at 53.3s, 53.6s, and 53.7s, respectively. They will always be less than the theoretical derivation timet1.Figure 8
The time responses of state variables with different S.Figures9(a)–9(d) show the variations of the system state variables δ and ω with different values of α and β. The simulated condition for the controllers is chosen to be the same. In Figures 9(a) and 9(b), where β = 1.5 and α = 0.6, 0.7, 0.8, 0.9, it is shown that, for both δ and ω, the smaller the value of α is, the faster the convergence rate is. In Figures 9(c) and 9(d), where α = 0.6 and β = 1.5, 1.6, 1.7, 1.8, it is shown that the law of β is the same as that of α. In other words, the convergence time is a decreasing function of the values of α and β for both δ and ω, which is consistent with the theoretical analysis of the maximum stable timeT of the system in the previous section. In addition, it is found from Figures 9(a) and 9(c) that the values of α and β affect the stability value of the system state variable δ. Therefore, we can get the stability value of the system state variable δ we need by adjusting the values of the controller parameters α and β.Figure 9
The time responses of state variables with different values ofα and β.
(a) (b) (c) (d)Figures10(a) and 10(b) are the comparisons of time responses of state variables in chaotic power system with the proposed controller, the fast terminal sliding mode (FTSM) controller, and finite-time controller. To make a fair comparison, the initial conditions, parameters, and the tuning parameters are identical. The simulation time is set to 20s, and the controllers are added at 10s. Figures 10(a) and 10(b) show that the settling time of proposed control scheme is smaller than that of the other two. Moreover, the time response curves of system state variables δ, ω under proposed control scheme have smaller maximum deviation than that under the other two, which demonstrates the proposed controller has better control effect. Table 1 records the exact numerical values of the settling time and maximum deviation in the simulation results, which verifies the conclusion more accurately.Table 1
Record of the numerical values of the settling time and the maximum deviation.
System state variable Settling time Maximum deviation δ Proposed control 12.06s -0.0236 FTSM control 12.76s -0.0480 Finite-time control 12.31s -0.0304  ω Proposed control 11.47s 0.0787 FTSM control 11.97s 0.1662 Finite-time control 11.81s 0.1353Figure 10
Time responses of state variables in chaotic power system with different controllers.
(a) (b)


## 5. Conclusion

In this paper, to investigate the problem of chaos suppression and voltage stabilization in chaotic power system, a new control scheme based on the fixed-time stability theory is proposed. An adaptive fixed-time fast terminal sliding mode chaos control strategy is presented to design controller. Simulation results illustrate the effectiveness and superiority of the proposed controller. Compared with the conventional fast terminal sliding mode control strategy and finite-time control strategy, the proposed controller has more advantages in the aspect of convergence time and maximum deviation. Moreover, the settling time is independent of the initial state and can be directly calculated. Note that the present study did not consider the effect of noise perturbation; future research will extend the proposed control strategy to high-order system with noise perturbation.



---

*Source: 5819428-2018-07-05.xml*