2018 International Conference on Sensor Network and Computer Engineering (ICSNCE 2018) 5 A New Non-Singular Terminal Sliding Mode Control and Its Application to Chaos Suppression in Interconnected Power System Wang Nian, Chen Hui, Ding Dawei school of Electronics and information engineering, Anhui university Hefei, 230601, China e-mail: wn_xlb@ahu.edu.cn Power Quality Engineering Research Center, Min of education Hefei, 230601, China e-mail: 1335965245@qq.com Abstract—Interconnected power system is a typical nonlinear dynamical system, it will cause great harm to the interconnected power system when the chaos occurred. This paper analyses the nonlinear dynamical behavior of the interconnected power system with a uncertain electromagnetic disturbance amplitude, and the influence of the electromagnetic disturbance amplitude on the stability of the system is obtained. Thus, this paper proposes a novel non-singular terminal sliding mode control to restrain the chaos, then the system will reach a stable state in a fixed time. According to the theoretical analysis, by introducing the saturation function, the control method can solve the singularity problem in the sliding mode control, and the interconnected power system will be stable in a short time when it’s in chaos. The simulation prove the correctness of the method. Keywords-Interconnected Power System; Saturation Function; Non-singularity; Sliding Mode Control; Fixed Time I. INTRODUCTION Power system is a kind of nonlinear dynamical system with multi-degree of freedom, strong coupling and multi-variable, which has rich dynamical behaviors. With the development of power grid system, grid interconnection has became a inevitable trend, and it can improve the quality of power. At the same time, it also provides convenience for the dispatching optimization of power system[1-3]. However, chaos often appears in the interconnected power system, it will bring great challenge to the stability of power system. For example, as the result of chaotic oscillations[4-6], several large area power outages appear in the United States, China and Canada in 1955. And it is difficult to suppress this phenomenon by using linear controller[7]. Therefore, it is necessary to study the mechanism of chaos in interconnected power systems and it is meaningful to design the nonlinear controller to suppress the chaos. Due to the high nonlinear characteristics of the interconnected power system, its stability is very sensitive to outside disturbance. If the perturbation is too large, its operating point will change obviously. There are also some tools to analyze its stability which include geometric method, energy function, bifurcation theory and numerical simulation[8]. Recently, there are a lot of scholars to study the stability of power grid system. For example, Nayfeh [9] uses the multi-scale perturbation method to study the stability of single machine power system and the bifurcation analysis of a single machine infinite power system is investigated by Duan[10]. Because coupling power angle exists in the interconnected power system, the inherent dynamical behavior of interconnected power system will be more abundant. Through the detailed numerical simulation, the influence of the conventional non-linearity index on the dynamic characteristics of the interconnected power system is expounded[11]. 2018 International Conference on Sensor Network and Computer Engineering (ICSNCE 2018) 6 In recent years, with the development of control method, the system is controlled from the single machine infinite system to multiple machine system, and the interconnected power system. In this paper, by proposing a fix-time non-singular terminal sliding mode control to realize stability when chaos occurred in interconnection power system. Some systems can be stabilized in a fixed time by using finite time control, this control method can be used in a lot of fields (for instance[12-13]). However, it is difficult to ensure the boundary convergence time when it is independent of the initial state. Therefore, in some practical system, it is not workable to use this control method when the initial condition is uncertain. Fortunately, this question was solved by Polyakov with the fixed-time stability theory[14]. Zuo[15] proposed a non-singular fixed-time terminal sliding mode controller for a class of second order nonlinear systems that can solve the singularity problem of terminal sliding mode controller in most instances. In this paper, by using the fast terminal sliding mode control can make the system convergence to steady state in finite time, the saturation function[16] and fast fixed time stability theory that can solve the singular problem through theoretical proof in this control method. This control method not only solve the singularity problem of the sliding mode controller, but the convergence speed is faster. Motivated by the above analysis, this paper investigated the dynamical characteristics and control of interconnected power system. Section 2 introduces the dynamical characteristics of the interconnected power system, by introducing the maximum Lyapunov index, power spectrum, phase diagram and timing diagram, the paper decribe the dynamic of the interconnected power system. when the amplitude of electromagnetic disturbance is v=1.3,the system is in chaos. a non-sigular sliding mode variable structure control method has been introduced in Section 3 and the effectiveness of the control method can be obtained by theoretical analysis, the advantages of this method are verified, and the convergence time is calculated in Section 4.And Section 5 gives the conclusion of this paper. II. THE ANALYSIS OF THE MODEL AND DYNAMIC CHARACTERISTICS OF THE INTERCONNECTED POWER SYSTEM There are two kinds of oscillation modes in interconnection power system, one is a single generator acts on other generators in the system with the frequency is between 0.5 and 2.0 Hz. The other oscillation mode is mainly expressed as a generator group in a region interacts with a generator group in another area, the frequency is between 0.1-1.0 Hz. This paper studies the interconnected power system model with two generators. Considering the influence of the amplitude of the electromagnetic disturbance power, the model is as follows: ( ) ( ), ( ) 1 [ sin( ( )) ( ) cos( ) sin( ( )) cos( )] s m k e d t t dt d t P t D t P P t t P t dt H                 Where, ( )t is the phase angle between the excitation potential and the terminal voltage between the two generators, and ( )t is the angular velocity of the two generators. D is the equivalent damping coefficient. , , , s m e k P P P P represents the amplitude of the electromagnetic power, the mechanical power, the load disturbance power, and the electromagnetic disturbance power, respectively. ,  is the electromagnetic power disturbance frequency and the load disturbance frequency. H is the equivalent moment of inertia. In order to analyses the dynamics of the interconnected power system, the simplified model is obtained: 1 2 2 1 2 1 ( ) ( ), ( ) sin( ( )) ( ) cos( ) sin( ( )) cos( ) dx x d dx x x v x d                     (2) Where 1 2 ( ) ( ), ( ) ( ) / , / , / , / , / , / , / , / s s m s k s e s s s s x t x t H P D HP P P v P P P P H P H P t P H                      Here, the parameters can be selected as 0.4, 0.02, 0.2, 0.8           . A. Lyapunov exponent The Lyapunov exponent is an important parameter of the system that can be used to measure and determine whether 2018 International Conference on Sensor Network and Computer Engineering (ICSNCE 2018) 7 the dynamic system is in the chaos. When the maximum Lyapunov exponent of the dynamic system is greater than 0, it can be concluded that the system is in chaotic oscillation state. Figure 1. The Lyapunov exponent B. Power spectrum To study the chaotic behavior of a system, it is effective to use power spectrum analysis methods. Actually, power spectrum analysis is through the time and space translate to the frequency space for the signal frequency structure. When the chaos occurs in the system, the power spectrum of the system behaves as a continuous irregular distribution, for example in fig.2(c). (a) (b) (c) Figure 2. Power spectrum Similarly, the phase space map also is a tool to determine whether the system is in chaos. (a) (b) (c) Figure 3. Phase plane plots 2018 International Conference on Sensor Network and Computer Engineering (ICSNCE 2018) 8 According to the maximum Lyapunov exponent, phase diagram, and power spectrum of interconnected power system, we can know that the interconnected power system is in an irregular non-periodic chaotic state when the parameter 1.3v  .This chaotic state will cause great harm to the stability of the system, it will lead to a large area of power outages. Therefore, it is necessary to study the method to restrain the chaos in interconnected power system. III. THE SLIDING MODE CONTROLLER DESIGN In this section, the paper proposes a control scheme that can restrain chaos in power system, this paper needs to add control law in the second item in control equation that can make the system output 1 x convergences to the control target, namely, 1 2 2 , ( ) ( , ) ( ) , x x x f x g x t d t u         (3) Where : 1 2 ( ) sinf x x x     : 1 ( , ) cos( ) sing x t v t x  , ( ) cos( )d t t  . By proposing a timing non-singular fast terminal control method, the system can be smoothly stabilized on the sliding surface s . The sliding surface s is designed as: 1 1 1 1 1 1 1 1 1 ( ) (| | 1) 2 2 2 2 1 1 1 1 1 m m p sign x n n q s x x x         (4) the paper can get the reach law : 2 2 2 2 2 2 1 1 ( ) (| | 1) 2 2 2 2 2 2 m m p sign s n n q s s s          Here, after reading the relevant literature [18 ] about how to overcome the singularity in terminal sliding mode control. This paper quotes saturation function that can solve the question. 1 1 1 1 1 1 1 2 2 2 2 2 2 1 1 1 1 1 ( ) (| 1|) 1 2 2 2 21 1 1 2 1 1 1 1 1 2 1 1 1 ( ) (| | 1) 2 2 2 2 2 2 1 1 [ ( ) ( 2 2 1 ( ) (| | 1)) ] 2 2 ( , ) m m sign x n n p q m m p sign s n n q m u f x dt n m sign x x x n p sat x x h q s s                         (6) In control input, the paper quotes saturation function to limit the amplitude of singularity term 1 1 1 1 2 p q x x  ,where the saturation function is , | | , ( , ) ( ), | | , x if x y sat x y y sign x if x y      (7) Theorem 1 In control law(6), if there is a positive number 1 2 1 2 , , ,    ,and also 1m , 1 n , 2 m , 2 n , 1 p , 1 q , 2 p , 2 q is odd positive integers satisfying 1 1m n , 2 2m n , 1 1p q , 2 2 p q , and 1 1 ( ) / 2m n , 2 2 ( ) / 2m n , 1 1 ( ) / 2p q , 2 2 ( )p q is odd positive integers. The system 11 will be stable in a fixed time. proof 2 1 1 2 V s  (8)  2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 1 ( ) (| | 1) 2 2 2 2 1 2 2 1 1 ( ) (| | 1) 1 1 2 2 2 2 2 2 1 1 ( ) (| | 1) 1 1 2 2 2 2 2 2 2 1 2 1 ( ) (2 ) (2 ) m m p sign s n n q m m p sign s n n q m m p sign s n n q V s s s s s s s V V                                     When | | 1s  2 2 2 2 1 1 2 2 1 2 1 2 1 (2 ) (2 ) m p n q V V V       (10) 2018 International Conference on Sensor Network and Computer Engineering (ICSNCE 2018) 9 When | | 1s  2 2 1 2 1 2 1 2 1 (2 ) (2 ) p q V V V      (11) Obviously, 1 0V  , so the control system will stable in expect target non-singularity. Theorem 2 The singularity item of the control input is restricted by the saturation function method, so that the system does not affect the stability analysis even if there is a singular region. Proof Defined inequality 1 1 1 1 1 1 2 1 | | p qp x x h q    as the singularity area. The state variable 1 x in first equation at system(3). 1 1 2 0 ( ) (0) ( ) t x t x x d   (12) When 2 ( ) 0x t  , 1 ( )x t will increase monotonically and leave the singularity. If 2 ( ) 0x t  , 1 ( )x t will decrease monotoniclly and leave the singularity. Therefore, the existence of the singular region does not affect the results of the stability analysis . Arrival time analysis Consider the following differential equation 1 1 ( ) (| | 1) 2 2 2 2 pm m sign y qn nx x x         (13) Where, assuming 0, 0   , and , , ,m n p q is odd positive integers, the system will stable in a fix time. Proof The above system can be written as follows: ,| | 1 ,| | 1 pm qn p q x x x x x x x x                  (14) A new variable is defined as 1 p q z x   , so, the first equation in the above system can be written 1 0 m p n q p qq p q p z z q q         (15) Let [( ) ] / [ ( )]m n q n q p    , can get 1 0 q p q p z z q q         (16) Similarly available, the second equation in system can be obtained 0 q p q p z z q q        (17) So, the maximum convergence time is 0 0 0 1 0 11 0 1 1 lim ( ) lim ( ) 1 1 1 1 ( ln(1 )) ln(1 ) z z z q T z dz dz q p z z q n q q p m n q p                                 (18) Without losing the general consideration of second order systems, the fix time for the controlled system to reach the slid surface is 2 21 1 2 2 2 2 2 2 2 1 2 2 2 2 2 2 2 2 21 1 ln(1 ) 2 2 2 m p n q n q T m n q p            (19) When the controlled system reaches the sliding surface 0s  , the target sliding mode of the system satisfies the following 1 1 1 1 1 1 1 1 1 ( ) (| | 1) 2 2 2 2 1 2 1 1 1 1 m m p sign x n n q x x x x          (20) The corresponding system 1 x will converge in a fixed time: 1 11 1 1 1 2 2 1 1 1 2 1 1 1 1 1 1 1 2 21 1 ln(1 ) 2 2 2 m p n q n q T m n q p            (21) The convergence time for system is 1 11 1 1 1 2 2 1 11 1 1 1 2 2 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 21 1 ln(1 ) 2 2 2 2 21 1 ln(1 ) 2 2 2 m p n q m p n q n q T T T m n q p n q m n q p                         (22) IV. SIMULATION EXPERIMENT The proposed control method is applied to suppress chaotic oscillation in studied power system. The parameters 2018 International Conference on Sensor Network and Computer Engineering (ICSNCE 2018) 10 of the controller is 1 10  , 1 5  , 2 10  , 2 5  , 1 9m  , 1 5n  , 2 9m  , 2 5n  , 100h  , 1 5p  , 1 9q  , 2 5p  , 2 9q  . The initial value of the controlled system is 1 2 [ , ] [0.5, 0.1]x x  .The dynamics of the system under this parameter have been obtained in the second part, and the system is in a chaotic state before it is controlled. As shown in figure 4, before being controlled, the system is in a chaos. Figure 4. Time domain waveform(uncontrolled) Figure 5. Time domain waveform after control Figure 6. Phase plane plots 2018 International Conference on Sensor Network and Computer Engineering (ICSNCE 2018) 11 Figure 7. Power spectrum From Fig.5, this paper can get that the system under chaotic oscillation is controlled by the control method proposed in this paper, then it will converge quickly to the desired target within a fixed time. By contrasting the Fig.6, Fig.7, this can find this result. The system will be in chaos when it’s uncontrolled, but the method which have proposed in this paper apply in the interconnected system, the system is stabilized. V. CONCLUSION In this paper, by plotting the Lyapunov exponent diagram, power spectrum and phase diagram of the interconnected power system, the influence of the amplitude of the electromagnetic disturbance for the system has been analyzed. According to the three stability criteria, when the system parameter 1.3v  , the interconnected power system will be in chaos. So, a non-singular terminal sliding mode control method with fixed time stability has been applied in the system when it’s in chaos. By comparing the system output, we can find that the control method proposed in this paper can restrain the chaotic oscillation, and the interconnection system is stabilized in a fixed time. The singularity problem in the terminal sliding mode control is eliminated by introducing the saturation function. Due to the timing convergence characteristics and non-singularity of the proposed method, it will be applied to the actual power equipment. VI. ACKNOWLEDGMENT The paper supported by National Natural Science Foundation of China (61201227) and Major science and technology projects in Anhui Province (16030901058). REFERENCES [1] Fuhong Min ,Yaoda wang, Guangya Peng, Enrong Wang,and Jane A.Auth,2016, “Bifurcations, chaos and adaptive backstepping sliding mode control of a power system with excitation limitation”.AIP Advances, 6, 085214 [2] Jyoti Ranjan Nayak, Tridipta Kumar Pati, Binod Kumar Sahu, Sanjeeb Kumar Kar, 2015, ICCPCT “Fuzzy-PID Controller Optimized TLBO Algorithm on Automatic Generation Control of a two-area Interconnected Power System” [3] Sat Sat Aung, Zaw Min Htike 2016 ASRJETS “Modeling and Simulation of Load Frequency Control for Three Area Power System Using Proportional Integral Derivative (PID) Controller” 26 2 301- 315 [4] Chen Ju-hua,Xu Nan.Application resrarch on small signal stability analysis of power systems”.Proceedings of the Csee, 2004, 24(12):52-57 [5] Du Quwei, Xiaoshu Luo.Passivity-based adaptive control of chaotic oscillations in power system”. chaos solitons&Fractals 2007. 31(3):665-671 [6] MM Zirkohi, Tkumbasar, T Lin, “Hybrid Adapive Type-2 Fuzzy Tracking Control of Chaoic Oscillation Damping of Power Systems” Asian Journal of Control 2017.19 [7] J Ni, L Liu, C Liu, X Hu.chattering-Free time scale separation sliding mode control design with application to power system chaos suppressin.Mathematical Problems in Engineering,2016:1-14. [8] Q Sun, Y Zhang, H He, D Ma, H Zhang.A Novel Energy Function-Based Stability Evaluation and Nonlinear Control Approach for Energy Internet.IEEE Transactions on Smart Grid,2017, PP(99) 1-16 [9] Ali Nayfeh DEAN T. MOOK Donald W.Lobitz 1974 12 9 Numerical Perturbation Method for the Nonlinear Analysis of Structural Vibrations 1222-1228 [10] Xiaodong Wang, Yushu Chen, Gang Han, Caiqin Song 2015 “Nonlinear dynamic analysis of a single-machine infinite bus power system Applied Mathematical Modelling “39 2951-2961 2018 International Conference on Sensor Network and Computer Engineering (ICSNCE 2018) 12 [11] Min Fu-Hong, Ma Mei-Ling Zhai Wei Wang En-Rong 2014” Chaotic control of the interconnected power system based on the relay characteristic function” Acta Phys. Sin. 63 050504 [12] X. Y. He, Q. Y. Wang, and W. W. Yu, “Finite-time containment control for second-order multiagent systems under directed topology,” IEEE Trans. Circuits Syst. II, Exp. Briefs, vol.61, no.8, pp.619-623, Aug.2014 [13] Y. Q. Wu, B. Wang, and G. D. Zong, “Finite-time tracking controller design for nonholonomic systems with extended chained form,” IEEE Trans. Circuits Syst. II, Exp. Briefs, vol.52, no.11, pp.798-802, Nov.2005 [14] A. Polyakov, “Nonlinear feedback design for fixed-time stabilization of linear control systems,” IEEE Trans. Autom. Control, vol.57, no.8, pp.2106-2110, Aug.2012 [15] E. Cruz-Zavala, J.A. Moreno, and L.M. Fridman, “Uniform robust exact differentiator,” IEEE Trans. Autom. Control, vol. 56, no. 11, pp.2727-2733, Nov.2011 [16] Z. Y. Zuo, “Non-singular fixed-time terminal sliding mode control of non-linear systems,” IET Contr. Theory Appl., vol.9, no.4, pp.545-552, Apr. 2015 [17] Junkang Ni, Ling Liu, Member, IEEE, Chongxin Liu, Xiaoyu Hu and Shilei Li 2016 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS-II: EXPRESS BRIEFS [18] B xu.Composite Learning Finite-Time Control With Application to Quadrotors.IEEE Transactions on Systems Man & Cybernetics Systems.2017, PP(99):1-9.