No Chattering and Adaptive Sliding Mode Control of a Fractional-Order Phase Converter with Disturbances and Parameter Uncertainties

Discussing the dynamical properties of various power system models is of significant importance in order to understand its complete behavior. Even though there are many literatures discussing about the chaotic behavior shown by phase converter circuits, none of them have reported the hazardous phenomenon of multistability. In this paper, we derive the fractional-order model of a phase converter circuit and investigate its dynamics. Bifurcation of the system with the parameters and fractional order are investigated. A forward and backward continuation scheme is adopted to display various coexisting attractors; the property of multistability is also discussed. Using forward and backward continuation, various coexisting attractors and the property of multistability are discussed. Two different sliding mode controllers for controlling chaotic oscillations with model disturbances and parameter uncertainties are derived, and the effectiveness of the controllers is discussed with numerical simulations.


Introduction
The nonlinear dynamical systems are described using nonlinear differential equations and represented in state (phase) space.The system state in steady state can be defined by equilibrium point and by limit cycles in state space; when the system is subjected to aperiodic oscillations, the system can be in quasiperiodic or in chaotic state.For different parameter values of a nonlinear system, the location and the number of equilibrium points can change.According to Lyapunov, regular attractors (steady states, limit cycles) are fully stable or fully unstable, but this is not true for chaotic attractors.Deterministic chaos is a phenomenon when small discrepancies in the initial conditions lead to an unpredictable behavior.The average rate of expansion along the principle axes are known as Lyapunov exponents.The presence of at least one positive Lyapunov exponent confirms chaos in the system.The chaotic behavior in electrical drive systems are noted as undesirable.
Most of the power systems are composed of ordinary differential equations; the system state changes (drastic changes in current or voltage decreasing or increasing rate) lead to the nonlinear dynamical behavior and makes the circuit modeling complex [1].The last two decades have witnessed multiple literatures discussing about the chaotic behavior of various power system models.To be specific, three cases, (i) high-frequency, time-sharing inverter, (ii) hysteresis current-controlled three-phase voltage source converter, and (iii) dual-channel resonant converter, were studied for existence and effectiveness of chaotic behavior in power electronics [2].Nonlinear behaviors, especially chaotic phenomena of the three-phase voltage source inverter (VSI) with hysteresis current comparator, are studied [3].Chaotic behavior of various types of switching converters is discussed in [2,[4][5][6][7]; the analytical model of switching converters with piecewise switched circuits is developed during the 1970s.The changes in the switching state in power converters results in nonsmooth system trajectories, and hence, it demands nonlinear methods for analysis.Basically, two modeling strategies, namely, continuous-time averaging approach (CTAA) and discrete-time iterative map approach (DTIA), are used to formulate converter models.The CTAA is simple and more suitable to characterize low-frequency phenomena whereas DTIA provides detailed information on the dynamical behavior of the system.Investigation on these approaches reveals that the CTAA is suitable to study the behavior when regular periodic motion collapsed into a quasiperiodic orbit of a much longer period [8].The interest on investigating switching converters increases since the system shows unusual bifurcation phenomena such as border collisions and grazings with usual properties of chaotic and quasiperiodic bifurcation.
Considerable works are identified on stabilizing such systems.Notable stabilization techniques that have been tested on electrical systems are time delayed approach [9,10], artificial intelligence techniques, nonlinear feedback, sliding mode, and adaptive backstepping.In 2000, Chen et al. [11] studied the time-delayed feedback approach to stabilize the chaotic dynamics in an experimental DC drive.Asakura et al. [12] opted neural networks to stabilize the chaotic dynamics of an induction drive system.Nonlinear feedback control [13], sliding mode control [14], and adaptive backstepping control [15] have all been applied towards stabilization of chaos in electrical systems.After the contribution of Wang on the universal approximation theorem, the effective utilization of fuzzy logic is recognized in adaptive controllers [16].A fuzzy logic-based adaptive controller is introduced in a MEMS triaxial gyroscope to approximate the model uncertainties and external disturbances [17].Recently, Fang et al. [18] investigated vibration control of a cantilever beam with an adaptive backstepping fuzzy sliding mode controller and ascribed that the sliding mode control complements an additional compensator for achieving system stability.In order to reduce the influence of chattering effect, a nested dynamic sliding mode controller is formulated and robustness is proved for a three-phase active power filter [19].An adaptive sliding mode control system using a double-loop recurrent neural network structure is designed and provided better approximation performance compared with the regular neural network approach [20].In [21], it is shown that an adaptive sliding mode control with the fuzzy neural network based on a radial basis function can promise that the system could reach the sliding surface and converge to equilibrium point asymptotically.
Fractional calculus is a generalization of ordinary differentiation and integration to an arbitrary (noninteger) order.Exponential laws are a classical approach to study dynamical changes in systems, but there are many systems where dynamics undergo faster or slower changes than exponential laws.Fractional-order systems play a significant role to describe the irregular changes in dynamics.Naturally, dynamical systems are composed in the form of fractional order, for example, the voltage-current relation of a semi-infinite lossy transmission line [22] or diffusion of the heat through a semi-infinite solid, where heat flow is equal to the halfderivative of the temperature [23].Recent studies reveal that fractional-order treatment allows us to describe and model a real object more accurately and precisely.Hartley et al. introduced the fractional-order Chua system with three single differential equations containing the noninteger (fractional) order derivatives [24].The fractional term can generate a new degree of freedom which can significantly improve the control performance as well as the parameter identification [25].
Chaos suppression research findings on the permanent magnet synchronous generator (PMSG) in wind turbine [26], permanent magnet synchronous motor (PMSM) [27], brushless DC motor [28,29], and smart grid [17] expose the unidentified chaotic regions while treating them as fractional orders.To the best of our knowledge, there is no power converter model with fractional order reported.In this paper, we formulated the fractional-order three-phase converter with saturable inductors from [28,29].

Phase Converter Circuit (PCC)
In this paper, we consider the phase converter circuit with saturable inductors as described in [28,29].The equivalent circuit of the single phase to the three-phase converter circuit of Ferro resonance type as given in [28] is shown in Figure 1.
The magneto motive force f φ is related to the flux ϕ by a nonlinear relation f φ = a 1 φ + a 3 φ 2 , and the other elements in the circuit shown in Figure 1 are assumed to be linear.The circuit equations are derived as magnetic and The equivalent circuit of the single-to three-phase converter circuit.
2 Complexity electrical subcircuits.The magnetic subcircuit equations can be derived as and the electrical subcircuit equations can be derived as We also assume that the variables are restricted to For simplicity, let us assume (symmetrical case) Using the relations discussed in ( 1), (2), and (3) and using the assumptions, the two state variables are taken as x = φ a − φ b = φ, y = φ c and τ = ωt; the state equations of Duffing type can be derived as To further derive the phase converter state equations, we use c 0 = c 1 = 0, c 3 = 1 and k = 0 1 as detailed in [28].The value of k is taken less than 0 3 to avoid nonperiod doubling bifurcations.Let a = k, b = c 3 /4, c = 3c 3 /4, d = B, e = B 0 and the state equations for the autonomous case can be derived as The PCC system described in (5) shows chaotic oscillations for the parameter values a = 0 1, b = 0 25, c = 0 75, d = 3, e = 0 2 as discussed in [28,29].

Fractional-Order Phase Converter Circuit (FOPCC)
This paper is aimed at analyzing the PCC system in its fractional-order form and also proposing sliding mode controllers to suppress chaotic oscillations with external disturbances and parameter uncertainties.To derive the fractional-order PCC (FOPCC), we have three commonly used definition of the fractional-order differential operator: Grünwald-Letnikov (GL), Riemann-Liouville, and Caputo [30][31][32].Compared to the other two methods, Grünwald-Letnikov gives a smooth derivative because of the recursively defined binomial coefficients [33][34][35][36][37].The continuous Riemann-Liouville approach is discretized using the GL method as seen in many literatures [23,38,39] and involves a discrete convolution between a "weight" or binomial coefficient function and the function of interest for differentiation.
As can be seen in the literatures, the Caputo method has been widely used for numerical solutions of the fractionalorder systems.But the GL method has benefits over the other methods of solving fractional orders due to the smoothness of the resultant approximations [40].Hence, we use the GL method to derive the FOPCC.The GL derivative can be defined as For numerical calculations, the above equation is modified as and the binomial coefficients required for the numerical simulation are calculated as Let the general form of the 3D fractional-order system be defined as In order to simulate system (9) using the GL method, we use the discretization method discussed [40]: where β is the binomial coefficient calculated using (8).The value of N is taken as the truncation window size and L is the memory element available and k when all the available memory elements are used.
Let us define the FOPCC oscillator as Using (10) in (11), the discrete form of the FOPCC system can be derived as For the parameter values a = 0 1, b = 0 25, c = 0 75, e = 0 2, fractional order q = 0 995, and initial conditions 1, 1, 1, 1, 1 , the FOPCC system (11) 2D phase portrait is shown in Figure 2 and chaotic oscillations also shown in Figure 2. In the entire numerical analysis of the FOPCC, the step size is h = 0 001.

Dynamical Properties of the FOPCC System
4.1.Bifurcation.We discuss the bifurcation of the FOPCC with respect to parameters and fractional order.The autonomous FOPCC system shown in (5) has five parameters.Our interest falls on parameter d for the reason that it represents the amplitude of the external excitation of the Duffing type equation.In [28,29], numerical bifurcation analysis was done, but system behavior for various values of the parameters is not sufficiently analyzed.Even though numerical bifurcation analysis was partly done in [28,29], fewer have been investigated on the system behavior for various values of the parameters.We hence derive the bifurcation plots with the other parameters taken as a = 0 1, b = 0 25, c = 0 75, e = 0 2, fractional order q = 0 995, and using d as the control parameter and the initial conditions taken as 1, 1, 1, 1, 1 .Figure 3 shows the bifurcation of the FOPCC with d.As can be seen from Figure 3, the FOPCC system shows period-1 limit cycles for 2 ≤ d ≤ 2 6, and at d = 2 6, the first period doubling is seen.For region 2 694 ≤ d ≤ 3 988, the FOPCC system shows the first chaotic zone and exits chaos and shows period-3 limit cycles for 3 989 ≤ d ≤ 7 694.We could see the second chaotic region for 7 8 ≤ d ≤ 8, and the FOPCC system exits chaos and shows period-3 oscillations for 8 001 ≤ d ≤ 8 756 and period-6 oscillations for 8 757 ≤ d ≤ 8 932.The third chaotic region of the system is seen for 8 98 ≤ d ≤ 9 56 and exits chaos with period-2 limit cycles.
Figure 3(b) shows the bifurcation of the FOPCC system with the fractional order q which is varied between 0 98, 1 while keeping the other parameters as a = 0 1, b = 0 25, c = 0 75, d = 3, e = 0 2 and the initial conditions taken as 1, 1, 1, 1, 1 .The FOPCC shows chaotic oscillations for q > 0 989, and we could see transient spirals in the region 0 98 ≤ q ≤ 0 986 which is mainly due to the effect of the truncation in the calculation of binomial coefficients.To remove such spiral transients, we could select a larger memory size which would result in longer simulation times.

Multistability and Coexisting
Attractors.Multistability and coexisting attractors in physical systems play a significant role as such oscillations may create hazardous effects in system behavior.There are few discussions on the multistable properties of mechanical systems such as the controllability of multistability in a quasiperiodically driven system [41][42][43] and multistability and coexisting attractors in a horizontal platform system [44].But there are no such investigations of coexisting attractors in an electrical system and especially in fractional-order physical systems.Hence, we are interested in investigating the multistability feature of the FOPCC system.

Complexity
To derive the multistability plots of the FOPCC, the bifurcation diagram is obtained by plotting local maxima of the coordinate x in terms of the parameter that is increased (or decreased) in tiny steps in the range of 2 ≤ d ≤ 10 with the parameters a = 0 1, b = 0 25, c = 0 75, e = 0 2 and the fractional order q = 0 995.The initial condition for the first iteration is taken as 1, 1, 1, 1, 1 and is reinitialized to the end values of the state variables in every iteration.This strategy, known as forward and backward continuation, represents a simple way to localize the window in which the system develops multistability.The existence of multistability can be confirmed by comparing the forward (Figure 4

Chaos Suppression in the Phase Converter Circuit
Recent literatures have shown that fractional-order controls are effective in suppressing chaos compared to their integer-order counterparts [1-8, 22, 24].Engineering systems which change under conditions demanding closed loops, various control schemes were used to get desirable output.PID controllers are a widely used control algorithm in general engineering applications.But systems with sensitivity in initial conditions and with drastic changes in parameters such as chaotic and hyperchaotic systems need an advanced control algorithm.An adaptive control scheme is introduced and shows effectiveness in suppressing chaos [1,6].The extended backstepping control algorithm provided a successful output in [2].Three famous control schemes named as sliding mode control, robust control, and extended backstepping control are implemented in fractional-order brushless DC motor to suppress chaos [4].The sliding mode control is free from chattering effects, so it is more suitable for control of complex nonlinear systems like the fractional-order Willamowski-Rössler chemical system [5].For systems suffering from external disturbances and/or parameter uncertainties, the sliding mode controller adapts to the sliding condition first so the adaptive sliding mode control algorithm is introduced [7].Nonlinear phenomena in complex systems such as voltage collapse and oscillatory phenomenon in power systems can be studied using the chaos theory.Chaos control methodologies have been of greater interest in nonlinearity suppression of power systems.Chaos control problems in fractional-order complex systems with higher harmonics like the induction motor drive system and smart grid systems are discussed where the control objective is achieved with a genetically optimized fractional-order PID controller [3,8].We discuss the chaos suppression problem of FOPCC with two different scenarios.In the first one, we assume that the FOPCC system has model uncertainties and external disturbances and design a fractional-order no-chattering sliding mode controller to suppress chaotic oscillations.In the second scenario, we assume that the FOPCC system has parameter uncertainties and design a fractional-order adaptive sliding mode controller to suppress chaos.In both cases, we use the Lyapunov approach to show that the designed controllers are effective in suppressing chaos in finite time t.

Chaos Suppression of the FOPCC System with
Disturbances.In this section, we derive the chaos control for the FOPCC system considered with external disturbances.We use the no-chattering sliding mode control defined in [24].Let us define the nonautonomous FOPCC system with the controllers as Let the nonlinear generalized equation of the FOPCC system (13) be defined as where A is the system matrix, C is the control matrix, d is the disturbance, b is the nonlinearity in the FOPCC system, and Δp is the parameter.Let us assume that the disturbances and the parameters are bounded to d ≤ α 1 and Δp ≤ α 2 , respectively.The desired signal for control can be defined as The error dynamics of the system defined in ( 14) is The proportional integral sliding surface can be defined as  6 Complexity where λ ∈ R 4×4 such that det λC ≠ 0 and L ∈ R 4×4 such that the eigenvalues of the matrix A − CL are negative and hence stable.
In the sliding surface, we know the condition that S = D q S = 0 and hence D q S = λD q e − λ A − CL e 17 Using ( 14) and ( 15) in ( 16), The controller satisfying the sliding condition can be defined as 7 Complexity where and h is a positive constant.

Stability of the Controller.
We use the Lyapunov approach to derive the stability of the controller and to show that the controller can make the FOPCC system reach the sliding condition.Let us define the Lyapunov function as The Lyapunov first derivative can be defined as The direct numerical solution of ( 22) is more complex, and hence, we use the modified Lyapunov method for the fractional order used in [4] to simplify (23) as Using ( 18) in ( 23), we get Using the designed sliding mode controller (19) in (24), Further approximating (25), we can derive (26) as 8 Complexity As per the boundedness of the disturbance and the parameters, If we assume S > h, then S i G S i /h = S i and ( 27) modify to If the positive constant μ > α 1 + α 2 + β and with β > 0, the controller (19) will make system (13) reach the sliding condition.
The Lyapunov derivative is negative semidefinite, confirming that there exists a finite time T where S ≤ h.For simplicity, the desired signal is taken as 0 5 sin t , 1 2 sin 5t , 1 5 sin 2t , 0 2 sin 1 5t .Figure 6 shows the time history of the controlled states with the controllers switched on at t = 20 s.

Chaos Suppression of FOPCC System with Parameter
Uncertainties.In this section, we assume that the FOPCC system has uncertain parameters and hence we derive an adaptive sliding mode controller which could estimate the parameters and suppress chaotic oscillations.Let us define the FOPCC system with adaptive sliding mode controllers as The uncertainties in the parameters are estimated with the parameter estimates â, b, ĉ, d, ê, and the parameter estimation errors are defined as The dynamics of the parameter estimates are derived as Let us define the proportional integral sliding surface as The fractional-order sliding dynamics can be derived from (35): D q y s y = D q y y + λ y y = −ay − bx 3 − cxz 2 + u y + λ y y, D q w s w = D q w w + λ w w = −aw − bx 2 z − cz 3 + d cos t + e + u w + λ w w Let us define the adaptive controllers which can make the FOPCC system with parameter uncertainties reach the sliding condition as u y = ây + bx 3 + ĉxz 2 − λ y y − η y sgn s y − ρ y s y , 37 where η y , η w , ρ y , ρ w , λ y , λ w are all positive constants.

Stability of the Controllers.
The stability of the controllers is established using the Lyapunov stability approach.The Lyapunov candidate function can be defined as The dynamics of the Lyapunov candidate function is From the definition of fractional calculus, V = s y D q y D 1−q y s y + s w D q w D 1−q w s w + e i D q i D 1−q i e i 41 The direct numerical simulation of ( 41) is complex, and hence, we use the method as that used to solve ( 22) and (41) can be simplified as V = s y −ay − bx 3 − cxz 2 + u y + λ y y + s w −aw − bx 2 z − cz 3 + d cos t + e + u w + λ w w + e a D q â + e b D q b + e c D q ĉ + e d D q d + e e D q ê , 42 where q is the commensurate fractional order of the system.
Using the controller (37) in ( 42 The Lyapunov first derivative in ( 45) is negative semidefinite confirming that the controller will make the FOPCC system reach the sliding condition in a finite time.
For numerical simulations, we use the initial conditions of the state variables as 0, 1, 0, 1 and the parameter estimates as â = 0, b = 0, ĉ = 0, d = 0, ê = 0.The commensurate fractional order is q = 0 995.Figures 7 and 8 show the time history of the parameter estimates and controlled states of the FOPCC system, respectively.

Conclusion
A fractional-order phase converter circuit is derived from its integer-order mathematical model, and various Figure 6: The time history of the states with the controller in effect when t ≥ 20 s and the initial conditions for the simulation is 0, 1, 0, 1 , and the commensurate fractional order is q = 0 995. 10 Complexity dynamical properties are discussed.Bifurcation analysis of the fractional-order system with forward and backward continuation shows the existence of multistability.Such coexisting attractors in a power system model were not explored earlier, and multistability existence is hazardous in any physical nonlinear system and needs to be further investigated.Coexisting attractors shown by the fractional-order system are shown but not limited to only those discussed as there may be other initial conditions and parameters that can contribute to coexistence.To suppress chaotic oscillations, two different scenarios are discussed.In the first scenario, we use a chattering free sliding mode controller to suppress chaotic oscillations with disturbances, and in the second scenario, we use an adaptive sliding mode controller to control chaos with parameter uncertainties.Numerical simulations are conducted to validate the effectiveness of the controllers.
(a) blue) and backward (Figure 4(a) red) bifurcation diagrams.As can be seen from Figure 4(a), a chaotic attractor coexists with a period-2 limit cycle in the regions 7 ≤ d ≤ 8, 8 001 ≤ d ≤ 8 756, and 9 89 ≤ d ≤ 10.The corresponding Lyapunov exponents (LEs) for forward (Figure 4(b)) and backward (Figure 4(c)) are also presented to confirm the findings on coexisting attractors.The LEs are calculated using the wellknown Wolfs algorithm [45] by replacing the ODE integrator

Figure 3 :
Figure 3: (a, b) The bifurcation of the FOPCC system with d.

Figure 4 :
Figure 4: (a) Bifurcation plots of the FOPCC system with parameter d where the forward continuation is shown with blue dots and backward continuation is shown with red dots.The corresponding LEs are shown in (b) and (c) for forward and backward continuation, respectively.

5. 1 . 2 .−bx 3 − cxz 2 0−bx 2 z 30 For
Numerical Analysis.For the FOPCC system, − cz 3 + d cos t , numerical simulations, the initial conditions of the nonautonomous FOPCC system are taken as 1, 1, 1, 1 and the disturbances are 0 8 sin 0 5t ; 1 2 sin 0 4t ; 0 4 sin t ; sin 2t .The gain matrix is taken as λ = diag 1, 1, 1, 1 , and the eigenvalues of A − CL are assumed to be −3, ), V = s y e a y + e b x 3 + e c xz 2 − η y sgn s y − ρ y s y + s w e 1 w + e b x 2 z + e c z 3 − e d cos t − e e − η w sgn s w − ρ w s w + e a â + e b b + e c ĉ + e d d + e e ê 43 Let us define the parameter update laws as D q â = −ys y − ws w , D q b = −x 3 s y − x 2 zs w , D q ĉ = −x 2 zs y − z 3 s w , D q d = cos t s w , D q d = s w 44 Using the Parameter update laws in (43), V = −η y s y − ρ y s 2 y − η w s w − ρ w s 2 w 45

Figure 7 :
Figure 7: The time history of the parameter estimates of the FOPCC system.

Figure 8 :
Figure 8: The time history of the controlled states of the FOPCC system.