Adaptive Sliding Mode Control for PMSG Wind Turbine Systems

In this paper, variable speed PMSG wind turbine systems with unknown system parameters, such as vicious friction coefficient and total inertia, are considered. The errors and variations of wind speed are modeled as a disturbance in mechanical torque. In general, the optimum rotating speed is given based on the MPPT (Maximum Power Point Tracking) algorithm and the designed controller tracks the reference (optimum) rotating speed in spite of these parametric uncertainties and disturbances. In order to have a desired rotor speed, a sliding mode current controller is proposed to have robustly stabilizing torque input. From the robustly stabilizing q-axis current iq, q-axis voltage input uq is obtained. Additionally, the d-axis control input ud is designed to regulate the d-axis current id. The adaptive estimator, for the total inertia J and the viscous friction coefficient F, is designed by a backstepping control technique. The robust stability of the closed-loop system is shown using a Lyapunov function. The proposed controller is verified via a simulation using MATLAB/Simulink.


Introduction
Wind turbines in wind power generation systems can be classified as inductive or synchronous depending on the rotor of the generator [1].Initially, the wind power generation industry used doubly-fed induction generators, for which the operating speed of the generator can be controlled using a slip.However, permanent magnet synchronous generators are increasingly being adopted for scenarios requiring high generation capacity, easy maintenance management, and cost efficiency, such as offshore wind power generation complexes.Permanent magnet synchronous generators (PMSG) have the advantages of reduced production costs due to developments in processing technology, reduced operating noise due to the use of permanent magnets in the rotors, reduced generator weight, and reduced maintenance costs.They also have the advantage of being able to create high torque output compared to doubly-fed induction generators due to high magnetic flux density and efficiency [2][3][4][5][6][7].
Wind turbines can be designed to operate at fixed speed or at variable speed [2,8].Fixed-speed operating systems transform wind energy into electrical energy using induction generators directly connected to a three-phase power grid.The rotor of the wind turbine is coupled to the generator shaft with a fixed ratio gearbox and operates at constant speed.By contrast, variable speed operating systems are systems which vary the rotation speed of the generator's rotors, and they are mainly used in synchronous generators.Variable operating wind power generators do not connect directly with the grid, instead connecting with it through a converter, and the generator is operated through aerodynamic control to adjust the torque, speed, and power.This reduces the mechanical stress and aerodynamic noise, and the generator can be controlled such that the wind power turbine operates at the maximum output coefficient over a wider range of wind speeds [2,9].
Here, the control method for energy production at maximum efficiency is called maximum power point tracking (MPPT) [10].An MPPT controller calculates the optimal rotor rotation speed for varying wind speeds.Therefore, the location and accuracy of the sensors which measure wind speed can be considered as important factors.To extract maximum power, linear controllers are designed based on an approximated linear model, such as conventional vector control with proportional-integral loops [3,9].However, these control strategies may not provide satisfactory performances due to the system nonlinearity of the PMSG and wide-range operation points.To improve the performance, a feedback linearizing control based MPPT is proposed [11], where the mechanical rotation speed controller and current controllers are designed via linear control methods.However, this results in a complex control law and has weak robustness against parameter uncertainties and external disturbances.Generally, wind power turbine rotors have very large diameters, and wind speeds vary according to the location on the turbine from which wind speeds are measured, so wind speed measurements may not be a suitable value for maximum energy production, depending on the measurement location.Therefore, wind speed measurement is more suitable for systems with severe noise due to the environmental uncertainty caused by measurement location and by noise inherent in their sensors.Because wind power turbine systems have significant uncertainty in their mechanical elements, control techniques are required to overcome this problem.
Sliding mode control (SMC) has been frequently used as a robust control for disturbance and uncertainty among the various robust control methods.Some papers on the application of SMC to wind energy conversion system have been presented in recent decades [11][12][13][14][15][16][17][18][19][20][21][22].In the early studies of [11,12], SMC is applied for MPPT in the wind energy conversion system with uncertainties.In [12], an optimal torque SMC strategy for a variable speed wind turbine system is proposed to implement MPPT tasks.In other words, SMC strategy was applied for controlling electromagnetic torque in MPPT for PMSG.In reference [11], an input-output linearization technique and SMC are applied to a wind energy conversion system with unstructured uncertainties.In [13], an SMC based on Enhanced Exponential Reaching Law was proposed and investigated on a grid-connected PMSG wind turbine system.In [14,15], SMC was used to improve robustness under various operating conditions such as parameter changes or load variation, and to reduce errors between desired command values and actual current for d-and q-axis currents.In [16], an induction generator (IG) speed drive was designed with the application of a sliding mode controller and a proposed artificial neural network controller.In [17], optimum torque and rotor speed were studied to be obtained by defining the appropriate sliding surface.In [18], an adaptive SMC was studied, first by designing an SMC for reducing the rotor speed error in PMSG, followed by a design for controller gain and mechanical torque estimation.
For uncertain nonlinear systems with known bounds, classical sliding mode control provides robustness.Recently, the adaptive SMC designs adapt the switching gain online without any predefined knowledge of the bound of uncertainty [23,24].To avoid over-estimation of gain, super twisting controllers [25,26] and Time-Delayed controllers [27][28][29][30] are proposed.[31] tracked the control of a class of uncertain nonlinear systems where the upper bound of the system uncertainty has explicit dependency on the system states.
In this paper, an adaptive robust sliding mode controller is proposed for a wind turbine system with parametric uncertainties and external disturbance.In the real system, some parameters, such as, mechanical inertia, viscous friction, stator resistance, inductance, are affected by operating conditions and manufacturing tolerance managed within predefined bounds.To improve the robustness and remove the effect of the noise and uncertainty which exist in wind speed measurements, an adaptive sliding mode controller is proposed.The wind speed sensor error is translated to an external torque disturbance with amplitude bound.To reject the torque disturbance, an SMC is designed for the nominal plant using the q-axis current as a virtual input.Next, the real physical control inputs, q-axis voltage and d-axis voltage, are designed to track the designed q-axis current and to regulate the d-axis current, respectively.An adaptation law is designed based on back-stepping technique while accounting for unknown parameters of the wind turbine system.The proposed controller is a dynamic controller which is composed of an SMC controller and parameter adaptation.The rotation speed error of rotor exponentially converges to the reference speed in spite of the external disturbance and unknown parameters.
The remainder of this paper is organized as follows.Chapter 2 briefly describes the wind energy conversion system in wind power generation systems.Chapter 3 defines the problem of noise and uncertainty in controller design, and describes the design of a controller to resolve this problem.Chapter 4 verifies the controller through simulations using MATLAB/Simulink, and Chapter 5 concludes the paper.

Modeling of a Wind Energy Conversion System
The wind power generation system model is shown in Figure 1.It is a system that converts wind energy into mechanical energy, and creates electrical energy through a generator.This system can be broadly divided into two parts: the generator side and the electrical grid side.This study focuses on controlling the generator side.accounting for unknown parameters of the wind turbine system.The proposed controller is a dynamic controller which is composed of an SMC controller and parameter adaptation.The rotation speed error of rotor exponentially converges to the reference speed in spite of the external disturbance and unknown parameters.The remainder of this paper is organized as follows.Chapter 2 briefly describes the wind energy conversion system in wind power generation systems.Chapter 3 defines the problem of noise and uncertainty in controller design, and describes the design of a controller to resolve this problem.Chapter 4 verifies the controller through simulations using MATLAB/Simulink, and Chapter 5 concludes the paper.

Modeling of a Wind Energy Conversion System
The wind power generation system model is shown in Figure 1.It is a system that converts wind energy into mechanical energy, and creates electrical energy through a generator.This system can be broadly divided into two parts: the generator side and the electrical grid side.This study focuses on controlling the generator side.

Wind Turbine Modeling
The kinetic energy of wind and the power of the wind can be expressed by Equations ( 1) and (2), where E is the wind's kinetic energy,  is the air density, A is the area that the wind passes through, v is the velocity of the wind, t is the time, w P is the power of the wind which can potentially be used, and b R is the radius of the wind power turbine through which the wind passes.The power of the wind expressed in Equation ( 2) only shows the maximum potential power that the wind has.In actuality, only part of this potential power can be used as electricity via the wind power turbine.The ratio of the wind's power and the mechanical power which can be generated by the turbine is the power coefficient p C .The maximum value of the output coefficient is 59.26%, which is referred to as Betz's limit; however, it actually exists in the range of approximately 25 to 45%, which can be expressed as follows [32,33]:

Wind Turbine Modeling
The kinetic energy of wind and the power of the wind can be expressed by Equations ( 1) and ( 2), where E is the wind's kinetic energy, ρ is the air density, A is the area that the wind passes through, v is the velocity of the wind, t is the time, P w is the power of the wind which can potentially be used, and R b is the radius of the wind power turbine through which the wind passes.The power of the wind expressed in Equation ( 2) only shows the maximum potential power that the wind has.In actuality, only part of this potential power can be used as electricity via the wind power turbine.The ratio of the wind's power and the mechanical power which can be generated by the turbine is the power coefficient C p .The maximum value of the output coefficient is 59.26%, which is referred to as Betz's limit; however, it actually exists in the range of approximately 25 to 45%, which can be expressed as follows [32,33]: where P m is the output mechanical power, λ is the tip speed ratio, β is the blade pitch angle, c 1 = 0.5, c 2 = 116, c 3 = 0.4, c 4 = 5 and c 5 = 21.The wind power turbine mechanical energy that can be extracted from the wind can be expressed via Equations ( 2) and (3) as shown below: Output power according to the rotor rotation speed for each wind speed is shown in Figure 2.
Energies 2019, 12, x FOR PEER REVIEW 4 of 18 where m P is the output mechanical power,  is the tip speed ratio,  is the blade pitch angle, 1 0.5 c  , c  .The wind power turbine mechanical energy that can be extracted from the wind can be expressed via Equations ( 2) and (3) as shown below: Output power according to the rotor rotation speed for each wind speed is shown in Figure 2.

C
that corresponds to the optimal value of p C [34].Therefore, where opt  and opt m  are the optimal values of tip speed ratio and rotor speed, respectively.By choosing the optimal values of tip speed ratio for the maximum power coefficient, maximum power can be extracted from the optimal rotor speed.It can be seen that the maximum output energy in Figure 2 varies according to wind speed.The generator's rotation speed at the maximum output energy is referred to as the optimal rotor speed (ω opt = ω re f ), and the technique for maintaining this speed at each wind speed, obtaining maximum output energy, is MPPT (Maximum Power Point Tracking) [10,34].According to the variation of β, typical characteristics of aerodynamic power coefficient correspond to tip speed ratio are illustrated in Figure 3.When β is maintained as a constant, we can see the power coefficient C p has only one maximum value C pmax that corresponds to the optimal value of C p [34].Therefore, C pmax can be expressed as where λ opt and ω opt m are the optimal values of tip speed ratio and rotor speed, respectively.By choosing the optimal values of tip speed ratio for the maximum power coefficient, maximum power can be extracted from the optimal rotor speed.
The wind power turbine's mechanical energy can be expressed as a product of the torque and rotation speed: where T m is the wind power turbine's mechanical torque, and ω m is the wind power turbine rotor's rotation speed.From Equations ( 5) and ( 7), the turbine's torque is The wind power generator's basic dynamic equation is where F is the viscous friction coefficient, J is the total inertia, and T e is the electromagnetic torque.
Energies 2019, 12, x FOR PEER REVIEW 5 of 18 The wind power turbine's mechanical energy can be expressed as a product of the torque and rotation speed: where m T is the wind power turbine's mechanical torque, and m  is the wind power turbine rotor's rotation speed.From Equations ( 5) and ( 7), the turbine's torque is The wind power generator's basic dynamic equation is where F is the viscous friction coefficient, J is the total inertia, and e T is the electromagnetic torque.

PMSG Modeling
Normally, the currents and the voltages in 3-phase systems are represented in dq frame using dq transformation.In dq transformation, the d axis is set as the generator rotor's N pole, and the q axis is perpendicular to the electrical angle, and according to the stator's 3-phase current and the electrical elements in a synchronously rotating coordinate axis, are transformed as shown in Figure 4.The PMSG kinetic equation, having been dq transformed, is shown below [1,4,7,9]: where d i is the d-axis current flowing to the stator, q i is the q-axis current flowing to the stator, d u is the input voltage for the stator's d axis, q u is the input voltage for the stator's q axis, m  is the generator's rotor speed, R is the resistance, L is the inductance, P is the number of pole pairs, and m  is the magnetic flux of the PMSG.

PMSG Modeling
Normally, the currents and the voltages in 3-phase systems are represented in dq frame using dq transformation.In dq transformation, the d axis is set as the generator rotor's N pole, and the q axis is perpendicular to the electrical angle, and according to the stator's 3-phase current and the electrical elements in a synchronously rotating coordinate axis, are transformed as shown in Figure 4.The PMSG kinetic equation, having been dq transformed, is shown below [1,4,7,9]: where i d is the d-axis current flowing to the stator,i q is the q-axis current flowing to the stator, u d is the input voltage for the stator's d axis, u q is the input voltage for the stator's q axis, ω m is the generator's rotor speed, R is the resistance, L is the inductance, P is the number of pole pairs, and ψ m is the magnetic flux of the PMSG.Then, the electromagnetic torque is where e T is the electromagnetic torque of the generator.If the inductance of each axis has the same value ( L L L   ), then Equation ( 12) can be simplified as Substitute Equation ( 13) into Equation ( 9), then the dynamics of the wind turbine are represented as Then, the electromagnetic torque is where T e is the electromagnetic torque of the generator.If the inductance of each axis has the same value (L = L d = L q ), then Equation ( 12) can be simplified as Substitute Equation ( 13) into Equation ( 9), then the dynamics of the wind turbine are represented as From Equations ( 10), (11), and ( 14), the whole model is In this model, the viscous friction coefficient F and the total inertia J are known a priori.However, in practical cases, the true values are not the same as the designed ones, because of production variation and changes according to time.Furthermore, the mechanical torque T m includes disturbances such as wind speed changes.In the steady state, the mechanical torque T m is obtained using wind speed v and rotating speed ω m from Equation (8).Thus, ∆T m represents the time varying signal with bounded amplitude, which is induced by the measurement error of v and the rotation speed error.For the stable and robust operation, we should consider the unknown parameters F and J in the controller design stage.

Controller Design
For the nominal turbine's mechanical torque T m , we can obtain from the measured wind speed v and rotor rotation speed ω m given by MPPT algorithm.However, the measured wind speed has some errors due to sensor noise and may not represent the effective wind speed according to the measurement sensor's location.When the torque input error is ∆T m , the actual mechanical torque input is expressed as the sum of the known nominal value T m and the disturbance ∆T m .Under these conditions, the controller should track the reference rotation speed ω * m to generate the maximum power.To have a robustly stabilizing controller, we find a reference current input to stabilize the speed error based on the sliding mode control technique.The obtained control is the reference signal of q-axis current i q to control the rotation speed.Since the q-axis current i q is not the control input, the q-axis voltage u q is designed to make i q track the reference signal.The other control input u q , the d-axis voltage, regulates the d-axis current i d .For the unknown parameters, a parameter estimator is designed based on backstepping.

Sliding Mode Controller Design for the Rotation Speed Error Regulation
A robust torque controller is designed based on the sliding mode control technique to overcome torque disturbances.In general, sliding mode controllers achieve excellent tracking performances and robustness against modeling uncertainty and disturbances.If we define the sliding surface S as a rotating speed error and the system satisfies stable sliding mode condition, then the rotating error slides on the sliding surface, i.e., rotating speed error remains 0. To make a system satisfy the condition of achieving a stable sliding mode, S • .S < 0, where S is the sliding mode plane, and .S is the time derivative of S [35][36][37][38].
If the uncertainties in the mechanical torque are taken into account, the wind turbine PMSG system is represented as Define the rotor's rotation speed error (speed tracking error) where ω * m is the reference rotation speed which is obtained by MPPT algorithm based on the measured wind speed, which may include sensor noises and wind speed variation.To have a torque controller, take the time derivative of Equation (17) .
For the simple notations, the quantity −1.5Pψ m i q is defined as u τ (t), which can be interpreted as torque control input to track the reference rotation speed ω * m .To have a robustly stabilizing sliding mode controller, we define the sliding surface as the speed tracking error Assumption 1.For the mechanical torque disturbance ∆T m (t), there exists a positive real number δ that satisfies Equation (20).
The Assumption 1 implies that the influence on the torque by wind speed variation, sensor noises, and model uncertainties is limited by a bounded magnitude.Given the practical environment, this assumption is reasonable and acceptable.Now, we propose a torque control input u * τ (t) for a given δ, where c 1 > 0, η > 0, and γ = δ + η √ 2 .
Proposition 1. Suppose that a system (18) satisfies assumption 1 for a positive real number δ.Then, the error system ( 18) is robustly exponentially stable by the feedback input u * τ (t) given as (21).
Proof of Proposition 1. Define a Lyapunov function candidate as Take the time derivative and by Assumption 1, we have Since .
V ≤ 0, the rotation speed error z 1 converges to 0 exponentially as time goes to infinity.
As in the conventional SMC, if we choose for some positive constant α, then we can show that .
Remark 1. From Equation ( 21), we have the reference q-axis current i * q to control the torque of PMSG as in the following Equation.

Voltage Controller Design and Adaptative Estimation for Unknown Parameters
In practical cases, many parameters are different from the designed values or changes over time.In this subsection, we assume that the viscous friction coefficient F and the total inertia J in Equation ( 16) have unknown values due to manufacturing tolerance.If the coefficients F and J are substituted with estimated values, the torque control input Equation ( 20) is ûτ where F is the estimated viscous friction F, and Ĵ is the estimated total inertia J.
Define the input error z 2 as Substitute Equation ( 26) into Equation ( 18); we have .
where J = J − Ĵ, and F = F − F. Take derivative Equation (25), From Equations ( 27) and ( 28), the time derivative of z 2 is .  û * τ (t).However, because of this approximation, the states are absolutely bounded instead of displaying exponential convergence.Please refer to Appendix A. Assumption 2. The variation ranges for total inertia J and viscous friction coefficient F are known, i.e.J ∈ [J min , J max ] and F ∈ [F min , F max ].
As observed in Equation ( 27) and ( 29), the error system contains control input voltage u q and time derivatives .Ĵ and .F. For some Lyapunov function, if we can make the time derivative of Lyapunov function be negative by choosing appropriate functions for u q , .Ĵ, and .F, then the error system is stable in the sense of Lyapunov.Additionally, since PMSG uses a permanent magnet, the d-axis current may be zero.Therefore, the d-axis current is regulated by control input voltage u d .In other words, if d-axis current tracking error is z 3 such that z 3 = i d − i * d then z 3 goes to zero as time goes to infinity.Now we propose the d,q-axis control input voltage u d , u q , and the estimator of Ĵ, F as follows: . .
Proposition 2. Suppose that Assumption 1 and Assumption 2 are satisfied.By the control inputs ( 30) and ( 31) and the parameter estimator ( 32) and ( 33), the tracking error states, z 1 , z 2 , z 3 exponentially decay to the origin and the estimated parameters Ĵ, F converge to some bounded values.
Proof of Proposition 2. Define a Lyapunov function candidate as Equation ( 34) is differentiated as follows: . V 2 z 1 .
By Assumption 1, we have Hence, .
V 2 can be written as: .
By the Cauchy-Schwarz inequality, . ( For c 1 , c 2 and c 3 such that V 2 ≤ 0. Thus, the error states z 1 , z 2 and z 3 decay to zero exponentially.Furthermore, if the error state Ĵ and . F are also zero.By LaSalle's Invariance Principle, we can say the estimator output Ĵ, F converges to some bounded point.Hence, it can be concluded that the proposed controller ( 30)-( 33) of the model of the wind energy conversion system with uncertainties (10), ( 11) and ( 16) guarantees that the rotating speed converges to the reference rotating speed ω * m exponentially.However, the parameter estimation converges to some values which may not the true ones.

Simulation Results
In order to verify the proposed adaptive sliding mode controller, a PMSG wind turbine system is simulated using MATLAB/Simulink, where the system parameters and the designed controller gains are presented in Tables 1 and 2, respectively.The torque input used in the simulation is depicted in Figure 5, which starts at 1000[N • m] and drops to 900[N • m] at time t = 1[sec].The reference rotor speeds are 75[rad/s] and 70[rad/s] for each torque value.The sinusoidal input 5 sin(44t) + 5 sin(20t) + 5 sin(52t) is also applied as an external disturbance ∆T m which is bounded by δ = 17.The initial values of the controller states F and Ĵ are assumed to be zero.
Table 1.The parameters of the wind turbine and the PMSG.

Parameter Value Unit
The resistance of the PMSG (R) 0.15 Ω The inductance of the PMSG (L)  For the performance analysis, we use an SMC controller designed based on nominal values 90 nom J  and 9 nom F  without parameter adaption.As the simulation results, each of the system states and the estimated values are shown in Figure 6.In Figure 6a, the controller wind turbine rotor speeds The control inputs q-axis voltage q u and d-axis voltage d u are depicted in Figure 6d,e.Due to the 2 z , chattering is observed during a reaching phase.The result of parameter estimator is shown in Figure 6d,e.As mentioned before, the estimated parameters converge to some bounded values which may different from the true ones.
(a) For the performance analysis, we use an SMC controller designed based on nominal values J nom = 90 and F nom = 9 without parameter adaption.As the simulation results, each of the system states and the estimated values are shown in Figure 6.In Figure 6a, the controller wind turbine rotor speeds ω m and the reference rotor speed ω * m are shown as a solid line and a dotted line, respectively.The reference rotor speed is initially 75[rad/s], and at time instant t = 1 [sec], it changes to 70[rad/s].In the presence of a torque disturbance and unknown parameters, the rotor rotating speed ω m exponentially converges to the reference speed ω * m in a short time using adaptive SMC.However, in the case of generic SMC, the rotor speed converged to 74.8[rad/s] and 69.8[rad/s] with remaining steady state error.The q-and d-axis currents i q , i d and reference values i * q , i

Conclusions
An adaptive robust sliding mode controller is proposed for a PMSG wind turbine system which has parametric uncertainties and external torque disturbance.The torque disturbance includes the wind speed measurement error and mechanical vibrations, which influence the mechanical torque.The control purpose is to track the reference rotor speed which is given by the MPPT algorithm and to regulate the d-axis current.To obtain a robust controller, a q-axis current is used as virtual control input and q-axis voltage is designed for the q-axis current to track the designed virtual input using an adaptive SMC.Regarding the d-axis current, d-axis voltage is designed for d-axis current regulation.For the unknown parameters J and F, an adaptation law is designed to make the closed loop system be Lyapunov stable.The performance of the controller is verified through generator rotor speed tracking via a simulation using MATLAB/Simulink (R2017b, MathWorks).The error states converge to the origin exponentially, and estimated parameters converged to some bounded values.In addition, for a practical implementation sigmoid function, ultimate boundedness of error system is shown in Appendix A.
In the future, it will be necessary to extend to a complete wind power generation system including a grid side converter connected to a grid.

Conclusions
An adaptive robust sliding mode controller is proposed for a PMSG wind turbine system which has parametric uncertainties and external torque disturbance.The torque disturbance includes the wind speed measurement error and mechanical vibrations, which influence the mechanical torque.The control purpose is to track the reference rotor speed which is given by the MPPT algorithm and to regulate the d-axis current.To obtain a robust controller, a q-axis current is used as virtual control input and q-axis voltage is designed for the q-axis current to track the designed virtual input using an adaptive SMC.Regarding the d-axis current, d-axis voltage is designed for d-axis current regulation.For the unknown parameters J and F, an adaptation law is designed to make the closed loop system be Lyapunov stable.The performance of the controller is verified through generator rotor speed tracking via a simulation using MATLAB/Simulink (R2017b, MathWorks).The error states converge to the origin exponentially, and estimated parameters converged to some bounded values.In addition, for a practical implementation sigmoid function, ultimate boundedness of error system is shown in Appendix A.
In the future, it will be necessary to extend to a complete wind power generation system including a grid side converter connected to a grid.
From Equations (A5) and (A6), the time derivative of z 2 is .
Now we propose the d,q-axis control input voltage u d , u q , and the estimator of Ĵ, F as follows: .
where φ and θ are positive scalars to determine each boundary layer thickness.
Proof of Proposition A2.Define a Lyapunov function candidate as Equation (A12) is differentiated as follows: . F.
V 4 can be written as: .
Then we can rewrite Equation (A16) as  Therefore, the states are ultimately bounded.

Figure 2 .
Figure 2. Maximum Power Point Tracking with variable wind speeds.It can be seen that the maximum output energy in Figure2varies according to wind speed.The generator's rotation speed at the maximum output energy is referred to as the optimal rotor speed ( opt ref   

Figure 2 .
Figure 2. Maximum Power Point Tracking with variable wind speeds.

Figure 3 .
Figure 3. Characteristics of Power coefficient versus Tip speed ratio for different blade pitch angles.

Figure 3 .
Figure 3. Characteristics of Power coefficient versus Tip speed ratio for different blade pitch angles.

Figure 4 .
Figure 4. Schematic diagram of the Permanent Magnet Synchronous Generator.

Figure 4 .
Figure 4. Schematic diagram of the Permanent Magnet Synchronous Generator.

) Remark 2 .
In practical control systems, it is impossible to implement the signum function sgn(•) and a generalized function δ(•).By replacing the signum function with sigmoid function tanh(•), we obtain similar results.If we use a sigmoid function tanh(•), the Dirac delta function δ(•) will not appear in

Figure 5 .
Figure 5.The mechanical torque with external disturbance.

m
and the reference rotor speed * m  are shown as a solid line and a dotted line, respectively.The reference rotor speed is initially 75[rad/s], and at time instant   1 sec t  , it changes to 70[rad/s].In the presence of a torque disturbance and unknown parameters, the rotor rotating speed m  exponentially converges to the reference speed * m  in a short time using adaptive SMC.However, in the case of generic SMC, the rotor speed converged to 74.8[rad/s] and 69.8[rad/s] with remaining steady state error.The q-and d-axis currents , 6b,c.The actual currents are solid lines, and the reference values are dotted lines.The q-, d-axis currents converge to the reference values in exponential way.The conversion rate is determined by choosing proper control parameters 1 2 3 , , c c c , as shown in the proof of Proposition 2.

Figure 5 .
Figure 5.The mechanical torque with external disturbance.
re f d are plotted in Figure 6b,c.The actual currents are solid lines, and the reference values are dotted lines.The q-, d-axis currents converge to the reference values in exponential way.The conversion rate is determined by choosing proper control parameters c 1 , c 2 , c 3 , as shown in the proof of Proposition 2. The control inputs q-axis voltage u q and d-axis voltage u d are depicted in Figure 6d,e.Due to the z 2 , chattering is observed during a reaching phase.The result of parameter estimator is shown in Figure 6d,e.As mentioned before, the estimated parameters converge to some bounded values which may different from the true ones.

10 Figure 5 . and 9 nomF 2 zFigure 6 .
Figure 5.The mechanical torque with external disturbance.For the performance analysis, we use an SMC controller designed based on nominal values 90 nom J  and 9 nom F  without parameter adaption.As the simulation results, each of the system states and the estimated values are shown in Figure 6.In Figure 6a, the controller wind turbine rotor speeds m  and the reference rotor speed * m  are shown as a solid line and a dotted line, respectively.The reference rotor speed is initially 75[rad/s], and at time instant

Figure 6 .
Figure 6.The Simulation results of the proposed adaptive SMC and SMC: (a) The rotor speed ω m and the reference rotor speed ω * m ; (b) The q-axis current i q and reference value i * q ; (c) The d-axis current i d and reference value i re f d ; (d) The q-axis voltage input u q ; (e) The d-axis voltage input u d and reference value i re f d ; (f) The real value of the total inertia J and estimator Ĵ; (g) The real value of the viscous friction coefficient F and estimator F.

Table 2 .
The control parameters.