Research Article  Open Access
Robust Control of Wind Turbines by Using Singular Perturbation Method and Linear Parameter Varying Model
Abstract
The maximum power point tracking problem of variablespeed wind turbine systems is studied in this paper. The wind conversion systems contain both mechanical part and electromagnetic part, which means the systems have time scale property. The wind turbine systems are modeled using singular perturbation methodology. A linear parameter varying (LPV) model is developed to approximate the nonlinear singularly perturbed model. Then stability and robust properties of the openloop linear singularly perturbed system are analyzed using linear matrix inequalities (LMIs). An algorithm of designing a stabilizing statefeedback controller is proposed which can guarantee the robust property of the closedloop system. Two numerical examples are provided to demonstrate the effectiveness of the control scheme proposed.
1. Introduction
Wind energy has become one of the fastestgrowing energies during the last two decades [1–3]. Because the wind turbines are large, flexible structures operating in noisy environment can offer abundant energy to mankind, and control problems are necessary to be studied to improve the reliability and conversion efficiency of the wind energy conversion systems (WECS).
Generally, wind turbines can be classified into four categories by [4], namely, fixedspeed fixedpitch turbines, fixedspeed variablepitch turbines, variablespeed fixedpitch turbines, and variablespeed variablepitch turbines. Compared to fixedspeed wind turbines, variablespeed ones can capture more energy and achieve better dynamic loads alleviation and fewer grid connection power peaks [4, 5].
Even though a wind energy conversion system can operate over a wide range of wind speeds, basically it is only active in two regions including partial load and full load [2]. In the partial load region, the control objective is to track the desired rotor speed corresponding to varying wind speed and maintain optimal tipspeed ratio [2, 6]. In the full load region, the generation goal is to limit the generated power to avoid overloading [4, 7]. Many researchers focus their interest on the control problems in the partial load region. In [8], an adaptive control scheme was developed for partial load region control of a variablespeed wind turbine, and the stability properties of the adaptive controller and the rotor speed were analyzed. Adaptive control approaches for maximizing power capture were also studied in [9, 10]. In [11], a twomass model and a wind speed estimator were used to propose a nonlinear controller, which was aimed to optimize the wind power capture. A highorder sliding mode controller is designed based on a high gain observer to optimize the capture wind energy by tracking the optimal torque in [12, 13].
The LPV reformulation of gain scheduling has become very popular during the last decades [4]. Compared with switching linear time invariant (LTI) controllers, LPV method can achieve better stability property. The LPV model was firstly introduced by Shamma and Athans [14] and then has been widely applied to wind turbine systems [15–18].
As known to all, the wind conversion systems include mechanical part and electromagnetic part, and the time scales of mechanical dynamics (slow) and electrical dynamics (fast) are quite different from each other [2, 19]. So the wind conversion systems have twotimescale property, which means these kinds of systems have high dimension and stiffness problems. The methodology of singular perturbation is always considered as a powerful tool to achieve dimension reduction, stiffness relief, and better control precision of multitime scale systems in reality [20–23]. Therefore, the singular perturbation technique is introduced to reduce the stiffness and to reach better tracking precision in this paper.
Motivated by above, in this paper, the control problems of variablespeed variablepitch turbines in the partial load region are considered. The rotor speed tracking problem of wind energy conversion system is studied to maintain the tipspeed ratio at optimal value . Namely, our goal is to keep rotor speed tracking the reference while the wind speed varies with time , where is the radius of the wind rotor plane.
Firstly, the nonlinear model of the WECS is presented. By extracting a small singular perturbation parameter, the nonlinear singularly perturbed system is obtained from the original nonlinear system. Then, the nonlinear singularly perturbed system is approximated by a LPV model. The stability properties of the openloop singularly perturbed LPV model are analyzed. Also, it is proved that holds for a given constant , where , where is the wind speed, is the operating point, is the tracking error, and is the transfer function from to . Next, the controller design method is developed using LMI skills. Then, an algorithm is proposed to demonstrate the whole process of designing a feedback controller for the singularly perturbed nonlinear system. In the end, two numerical examples demonstrate the effectiveness of the algorithm obtained in this paper. The whole control scheme is depicted in Figure 1.
The main contributions of this paper are outlined as below:(1)Singular perturbation methodology and LPV model are combined, for the first time, to solve the maximum power point tracking problem of wind energy conversion system. Compared to switching LTI controllers, LPV model can achieve better stability property. And, singular perturbation methodology is a wonderful tool to reduce the stiffness and improve the tracking precision.(2)The stability and robust property of the singularly perturbed LPV systems are proved under certain conditions. Besides, an algorithm is presented to design a robust stabilizing control scheme for the singularly perturbed LPV system. Using LMI technology, the stability and robust property are proved for the closedloop system.
The rest of this paper is organized as follows. In Section 2, the nonlinear model of the WECS is presented which involves the dynamic characteristics of the mechanical part and electrical part. Then, a singular perturbation parameter is extracted which leads to a singularly perturbed nonlinear model. The nonlinear singularly perturbed model is reformed into a singularly perturbed LPV model that contains a family of LTI models. Later, Lyapunov theory is utilized to prove that the tracking error between the rotor speed and the desired rotor speed decays to zero in Section 3. Furthermore, for a given constant , is proved. In Section 4, the conditions in LMI forms are given to design a robust stabilizing parameterdependent controller. Section 5 presents an algorithm to design a parameterdependent controller. In Section 6, two numerical examples are given to illustrate the effectiveness of the results obtained. Finally, conclusions are drawn in Section 7.
2. System Description
In this section, the mathematical model of the wind turbine systems with permanent magnet synchronous generators is developed.
Because the electrical part of the system changes much faster than the mechanical part, namely, the states of the wind energy conversion systems have two different time scales. A singular perturbation parameter is extracted to obtain the singularly perturbed nonlinear model. The following LPV technique is used to linearize the singularly perturbed nonlinear system.
Our consideration is focused on partial load. In the partial load region, the wind rotor speed is adjusted to maintain the optimal tipspeed ratio as the wind speed changes. For this control purpose, only aerodynamics, drive train dynamics, and generator dynamics are taken into account.
Commonly, the aerodynamic torque is given as follows [24, 25]:where is the air density, is the wind speed, is the radius of the wind rotor plane, and power coefficient is approximated by a secondorder polynomial of tipspeed ratio [4] as below:where is the maximum power coefficient, is the optimal tipspeed ratio corresponding to , and is defined aswhere is the wind rotor speed.
The drive train block has the model as below [2, 26]:where is the generator speed, is the internal torque, is the wind rotor inertia, is the generator inertia, is the stiffness coefficient of the highspeed shaft (the generator shaft), is the damping coefficient of the highspeed shaft (the generator shaft), is the gearbox ratio, and is the gearbox efficiency.
Then, the generator dynamics are modeled as [2, 25]where is the generator electromagnetic torque, , , and , , are the and components of the stator current, inductance, and voltage, respectively, is the stator resistance, is the number of pole pairs, and is the flux.
By combining (4) and (5), the complete nonlinear model of the wind energy conversion system is obtained. Next, we will introduce the singular perturbation method.
Considering the order of magnitude of and , we define the socalled singular perturbation parameter and obtainwhere , .
Now the singularly perturbed nonlinear system is obtained by uniting (4) and (6).
And now we are ready to derive the linear model using LPV method. Choose an operating point and linearize the nonlinear parts in the singularly perturbed nonlinear system at point :where
By substituting (7) into (4) and (6), we have the following singularly perturbed linear system at the operating point :whereFurthermore, note , , and rewrite (15) aswhereThen, by appropriately choosing operating points ,is a convex polytope with being vertices, . Note that the LPV model (17), with and approximated by (11) and (12), is affine in the parameters. Namely, there exist scalars , , , and such that and . Hence, it is easy to see that, for any , there exists a set of positive numbers , , such thatwhere .
Therefore, for any , we have derived the singularly perturbed LPV model as below:
Remark 1. For the details of skills to choose operating points appropriately, please refer to [4].
Remark 2. The wind conversion system considered here contains both mechanical and electrical parts which are of different time scales. Singularly perturbed model is developed of this system and LPV method is used to reform the model. The combination of singular perturbation theory and LPV method is novel for the maximum power point tracking problem.
3. Stability Analysis of OpenLoop System
In this paper, the objective is to maintain the optimal tipspeed ratio by adjusting wind rotor speed as the wind speed changes. So the operating points are chosen such that the tipspeed ratio is optimal. If the states and in (21) decay to zero, it means that the errors between the actual states and the desired states tend to zero. In this case, the wind turbine runs to extract all the available power.
And, this section will analyze the stability of the singularly perturbed system (21) when the control input based on LMI technique.
When , from (21) we can have
Before the main result, Schur Complement Lemma is given below.
Lemma 3 (Schur Complement Lemma). The partitioned matrixif and only if
Theorem 4. For the singularly perturbed system (22), given a positive scalar , if there exist five positive matrices satisfying the following conditions:where , , andThen, the equilibrium point of system (22) at point is asymptotically stable, and is satisfied.
Proof. Construct a Lyapunov functionwhere , . Since and , we have , and obviously holds.
Firstly, the asymptotic stability of the system is proved under conditions in Theorem 4 when the disturbance is zero. Derive with the respect to along the trajectory of (22), with , and we obtainFrom the LMI (25), it is not difficult to have the following inequality:And, by adding the weight values (i.e., ) of to (29), we can getTherefore,Hence, the system of (22) at point is asymptotically stable when .
Next, the robust property of system (22) will be proved when the disturbance .
Because , using the LMI (25) it can be obtained thatwhich leads toThen by using Schur Complement Lemma twice, LMI (33) can be transformed into inequality below:Derive along the trajectory of (22) and use the inequality (34); becomesBased on the asymptotic stability proved at the first part of this proof, we have . Assume , and integrate both sides of (35) from to ; we can haveThen, (36) can be rearranged to haveTherefore, as can be seen from the inequality above,Now, it is easy to see that is satisfied. This completes the proof.
Remark 5. Even though Theorem 4 can guarantee the stability and robust property of system (22), the LMI (25) is dependent on small parameter , so it might be singular and difficult to solve.
The following result improves the illcondition inequality.
Theorem 6. For the singularly perturbed system (22), given a positive scalar , if there exist five matrices satisfying the following conditions,where , , andthen, the equilibrium point of system (22) at point is asymptotically stable, and holds.
Proof. Define a Lyapunov function as below:where and .
According to condition (39), we haveAs a consequence, is satisfied.
Then, the timederivative of along the solution of (21) is given byThe following part is similar to the proof of Theorem 4 and is consequently omitted here.
4. Controller Design
In this section, a robust statefeedback controller is designed for system (21), and the stability property of the closedloop system is analyzed.
Since the coefficients of system (21) depend on , it is reasonable to design a controller whose feedback gain matrix also depends on .
As can be seen from (20) and (22), at any point , the dynamic equation of the nonlinear singularly perturbed system can be expressed as a weightsum of the dynamic equations at the vertices of , . Therefore, we will design controllers for the LTI systems operating at the vertices of , , and use a weightsum of the controllers at vertices as the control input to the system at point .
At the vertex point , , design a robust statefeedback controller as below:The closedloop system of (22) at the vertex point with (45) as control input is as follows:where .
Then, for the nonlinear singularly perturbed system (21) at , the statefeedback controller is as follows:where and .
Applying (47) to system (21), the closedloop system is obtained as below:
The properties of closedloop system of (48) are analyzed in Theorem 8. To prove Theorem 8, the following lemma from [27] is needed.
Lemma 7 (see [27]). Let , , and be a scalar; then
Theorem 8. The closedloop system of (21) with (47) as input is asymptotically stable; namely, (48) system is asymptotically stable, and is satisfied for a given , if there exist matrices and of appropriate dimension, such that the following LMIs hold:where .
Proof. Similar to the proof of Theorem 4, andAccording to (51), we can haveLet . From (50) and (53) we can obtain , andBy applying Schur Complement Lemma four times to (54), it yieldsAccording to Lemma 7, it is easy to obtainHence, from (55) and (56), we can getwhich is equal toDefine a Lyapunov function as below:where . Then, the timederivative of along the solution of (21) is given byThe following part is similar to the proof of Theorem 4 and is consequently omitted here.
5. Algorithm of Synthesis
In order to clarify the whole process of designing a parameterdependent controller for the original nonlinear system (4) and (5), the following algorithm is presented.
Step 1. Choose five operating points .
Step 2. Abstract singular perturbation parameter , and get the nonlinear singularly perturbed systems (4) and (6).
Step 3. Linearize the nonlinear singularly perturbed systems (4) and (6) at , and obtain the linear parameterdependent coefficients and , where .
Step 4. For a given , at each operating point , LMIs (50) and (51) are solved to obtain control gain matrices , .
Step 5. At time , the variable is measured, and the weighting coefficients satisfyingare computed with and .
Step 6. The control gain matrix at time is obtained as below:
Step 7. Apply the controller (62) to the original nonlinear system (4) and (5).
Step 8. At time , repeat Step 5 to Step 8.
Remark 9. This algorithm can only be applied to the point , where is defined by (19). If , it can not be guaranteed that we can compute , , with satisfied.
Remark 10. Since the operating point involves wind speed , implies the range of the wind speed within which our algorithm is effective.
Remark 11. In [12, 13], a highorder sliding mode control strategy is proposed based on a high gain observer to optimize the maximum power point tracking problem of wind energy conversion system. Their method presents chatteringfree, behavior, finite reaching time and robustness. However, highorder controller and observer need more online computation. In this paper, the wind energy conversion system is modeled using singular perturbation theory which considers the generator dynamics as fast subsystem and the drive train block as slow subsystem. This method combines the singular perturbation methodology and LPV model for the first time to solve the MPPT problem.
6. Numerical Examples
In this paper, the goal is to track the desired wind rotor speed and maintain the optimal tipspeed ratio when the wind speed changes. The simulation study is performed to verify the effectiveness of the proposed control algorithm. The two examples consider the CART 3blades wind turbine taken from [20] as an objective. The parameters of the wind turbine are given in Table 1. The experiment is carried out on the MATLAB FAST® software which is developed by the American National Renewable Energy Laboratory.

Example 12. Consider a wind turbine system with the parameters depicted in Table 1. The wind speed is assumed as a constant in this example. The tracking result controlled by the algorithm developed in Section 5 is compared with that of the optimal torque (OT) method, and the compared rotor tracking results are shown in Figure 2.
It can be seen that the tracking error between actual rotor speed controlled by the algorithm presented here and the desired rotor speed decays to zero after a small overshoot. But the tracking error between the desired rotor speed and the rotor speed controlled by the optimal torque method can not decay to zero with time, which means a static error exists. So this example illustrates that the algorithm developed in this paper is more effective than the optimal torque method.
Example 13. In this example, a turbulence of 600 seconds produced by using TurbSim software is applied. This turbulence satisfies the IEC614001 standard and the turbulence is shown in Figure 3.
The rotor speeds controlled by the algorithm presented in this paper and the optimal torque method are shown in Figure 4. It is obvious that the rotor speed controlled by the algorithm can track the desired rotor speed much better than that of the optimal torque method.
Furthermore, the wind power capture efficiency and energy conversion efficiency controlled by this novel method and optimal torque method are compared in Table 2. It can be seen from Table 2 that the method developed can obtain better wind power capture efficiency and energy conversion efficiency.

7. Conclusions
This paper extended the continuoustime infinite horizon nonlinear quadratic optimal control problem of NSPSs to discretetime version with the weight matrices dependent on the states in the cost function. For a class of discretetime NSPSs in this paper, we used the theory of singular perturbations and time scales to decouple the original highorder NSPS into orderreduced slow and fast (boundary layer) subsystems. Then, via the statedependent Riccati equation, suboptimal controllers for the two subsystems are designed with the weight matrices varying with states in the cost functions. A composite controller consisting of two suboptimal controllers is developed for the original system. It is proved that the equilibrium point of the original closedloop system with a composite controller is locally asymptotically stable. In the end, an example is given to show the effectiveness of the results obtained.
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
This work was supported by the Natural Science Foundation of Jinling Institute of Technology (40620029) and the Natural Science Foundation of Jiangsu Province (15KJB520010).
References
 K. E. Johnson, L. Y. Pao, M. J. Balas, and L. J. Fingersh, “Control of variablespeed wind turbines: standard and adaptive techniques for maximizing energy capture,” IEEE Control Systems Magazine, vol. 26, no. 3, pp. 70–81, 2006. View at: Publisher Site  Google Scholar
 H. M. Nguyen and D. S. Naidu, “Time scale analysis and control of Wind Energy Conversion Systems,” in Proceedings of the 5th International Symposium on Resilient Control Systems (ISRCS '12), pp. 149–154, Salt Lake City, Utah, USA, August 2012. View at: Publisher Site  Google Scholar
 P. Musgrove, Wind Power, Cambridge University Press, New York, NY, USA, 2010.
 F. D. Bianchi, H. D. Battista, and R. J. Mantz, Wind Turbine Control Systems Principles, Modelling and Gain Scheduling Design, Springer, London, UK, 2007.
 B. Boukhezzar and H. Siguerdidjane, “Nonlinear control of variable speed wind turbines without wind speed measurement,” in Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference (CDCECC '05), pp. 3456–3461, Seville, Spain, December 2005. View at: Publisher Site  Google Scholar
 J. H. Laks, L. Y. Pao, and A. D. Wright, “Control of wind turbines: past, present, and future,” in Proceedings of the American Control Conference (ACC '09), pp. 2096–2103, IEEE, St. Louis, Mo, USA, June 2009. View at: Publisher Site  Google Scholar
 M. Liserre, R. Cárdenas, M. Molinas, and J. Rodríguez, “Overview of multiMW wind turbines and wind parks,” IEEE Transactions on Industrial Electronics, vol. 58, no. 4, pp. 1081–1095, 2011. View at: Publisher Site  Google Scholar
 K. E. Johnson, L. Y. Pao, M. J. Balas, V. Kulkarni, and J. L. Fingersh, “Stability analysis of an adaptive torque controller for variable speed wind turbines,” in Proceedings of the 2004 43rd IEEE Conference on Decision and Control (CDC '04), pp. 4087–4094, Paradise Island, Bahamas, December 2004. View at: Google Scholar
 K. E. Johnson and L. J. Fingersh, “Adaptive pitch control of variablespeed wind turbines,” Journal of Solar Energy Engineering, vol. 130, no. 3, Article ID 031012, 2008. View at: Publisher Site  Google Scholar
 Y. D. Song, B. Dhinakaran, and X. Y. Bao, “Variable speed control of wind turbines using nonlinear and adaptive algorithms,” Journal of Wind Engineering and Industrial Aerodynamics, vol. 85, no. 3, pp. 293–308, 2000. View at: Publisher Site  Google Scholar
 B. Boukhezzar and H. Siguerdidjane, “Nonlinear control of a variablespeed wind turbine using a twomass model,” IEEE Transactions on Energy Conversion, vol. 26, no. 1, pp. 149–162, 2011. View at: Publisher Site  Google Scholar
 M. Benbouzid, “Highorder sliding mode control of DFIGbased wind turbines,” in Wind Turbine Control and Monitoring, Advances in Industrial Control, pp. 23–48, Springer, Berlin, Germany, 2014. View at: Google Scholar
 B. Beltran, M. Benbouzid, T. AhmedAli, and H. Mangel, “DFIGbased wind turbine robust control using highorder sliding modes and a high gain observer,” International Review on Modelling and Simulations, vol. 4, no. 3, pp. 1148–1155, 2011. View at: Google Scholar
 J. S. Shamma and M. Athans, “Guaranteed properties of gain scheduled control for linear parametervarying plants,” Automatica, vol. 27, no. 3, pp. 559–564, 1991. View at: Publisher Site  Google Scholar  MathSciNet
 F. D. Bianchi, R. J. Mantz, and C. F. Christiansen, “Gain scheduling control of variablespeed wind energy conversion systems using quasiLPV models,” Control Engineering Practice, vol. 13, no. 2, pp. 247–255, 2005. View at: Publisher Site  Google Scholar
 F. Wu, X. H. Yang, A. Packard, and G. Becker, “Induced L2norm control for LPV system with bounded parameter variation rates,” in Proceedings of the American Control Conference, vol. 3, pp. 2379–2383, Seattle, Wash, USA, 1995. View at: Google Scholar
 F. Lescher, J. Y. Zhao, and P. Borne, “Switching LPV controllers for a variable speed pitch regulated wind turbine,” in Proceedings of the IMACS Multiconference on Computational Engineering in Systems Applications, pp. 1334–1340, October 2006. View at: Publisher Site  Google Scholar
 C. Sloth, T. Esbensen, and J. Stoustrup, “Active and passive faulttolerant LPV control of wind Turbines,” in Proceedings of the 2010 American Control Conference (ACC '10), pp. 4640–4646, Baltimore, Md, USA, July 2010. View at: Google Scholar
 Y. Zhang, D. S. Naidu, H. Nguyen, C. Cai, and Y. Zou, “Time scale analysis and synthesis for Model Predictive Control under stochastic environments,” in Proceedings of the 7th International Symposium on Resilient Control Systems (ISRCS '14), Denver, Colo, USA, August 2014. View at: Publisher Site  Google Scholar
 H. M. Nguyen, Advanced control strategies for wind energy conversion systems [Ph.D. thesis], Idaho State University, Pocatello, Idaho, 2013.
 C. Cai, Z. Wang, J. Xu, X. Liu, and F. E. Alsaadi, “An integrated approach to global synchronization and state estimation for nonlinear singularly perturbed complex networks,” IEEE Transactions on Cybernetics, vol. 45, no. 8, pp. 1597–1609, 2015. View at: Publisher Site  Google Scholar
 D. S. Naidu and A. J. Calise, “Singular perturbations and time scales in guidance and control of aerospace systems: a survey,” Journal of Guidance, Control, and Dynamics, vol. 24, no. 6, pp. 1057–1078, 2001. View at: Publisher Site  Google Scholar
 D. S. Naidu, “Singular perturbations and time scales in control theory and applications: an overview,” Dynamics of Continuous, Discrete and Impulsive Systems Series B: Applications and Algorithms, vol. 9, no. 2, pp. 233–278, 2002. View at: Google Scholar  MathSciNet
 F. A. Shirazi, K. M. Grigoriadis, and D. Viassolo, “Wind turbine integrated structural and LPV control design for improved closedloop performance,” International Journal of Control, vol. 85, no. 8, pp. 1178–1196, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 C.N. Wang, W.C. Lin, and X.K. Le, “Modelling of a PMSG wind turbine with autonomous control,” Mathematical Problems in Engineering, vol. 2014, Article ID 856173, 9 pages, 2014. View at: Publisher Site  Google Scholar
 I. Munteanu, A. Bratcu, N. Cutuluslis, and E. Ceanga, Optimal Control of Wind Energy Systems: Toward a Global Approach, Springer, 2008.
 L. E. Ghaoui, E. Feron, and V. Balakrishnan, Linear Matrix Inequalities in System and Control Theory, Society for Industrial and Applied Mathematics, Philadelphia, Pa, USA, 1994.
Copyright
Copyright © 2016 Yan Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.