Mathematical Problems in Engineering

Volume 2014 (2014), Article ID 409268, 8 pages

http://dx.doi.org/10.1155/2014/409268

Research Article

Dynamic Model of a Wind Turbine for the Electric Energy Generation

José de Jesús Rubio,¹ Luis Arturo Soriano,² and Wen Yu²

¹Sección de Estudios de Posgrado e Investigación, ESIME Azcapotzalco, Instituto Politécnico Nacional, Avenida de las Granjas No. 682, Colonia Santa Catarina, 02250 México, DF, Mexico

²Departamento de Control Automatico, CINVESTAV, Instituto Politécnico Nacional, 07360 México, DF, Mexico

Received 5 November 2013; Accepted 7 January 2014; Published 26 February 2014

Academic Editor: Feng Ding

Copyright © 2014 José de Jesús Rubio 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.

Abstract.

A novel dynamic model is introduced for the modeling of the wind turbine behavior. The objective of the wind turbine is the electric energy generation. The analytic model has the characteristic that considers a rotatory tower. Experiments show the validity of the proposed method.

1. Introduction

Researchers are often trying to improve the total power of a wind turbine. The dynamic model of a wind turbine plays an important role in some applications as the control, classification, or prediction.

Some authors have proposed the equations to model the dynamic behavior of the wind turbine as shown in [1–5].

This paper presents a novel dynamic model of a wind turbine; the first dynamic model is for the wind turbine, the second is for the tower, and both are related. Because the rotatory tower can turn, it may help the wind turbine to increase the air intake. Some companies propose wind turbines that consider rotatory blades; however, in this study, a rotatory tower is proposed instead of the rotatory blades because the first dynamic model is easier to obtain than the second.

The proposed model is linear in the states as the works analyzed by [6–10] show. The technique provides an acceptable approximation of the wind turbine behavior.

The paper is structured as follows. In Section 2, the dynamic model of a windward wind turbine of three blades with a rotatory tower is introduced. In Section 3, the simulations of the dynamic model are compared with the real data obtained by a real wind turbine prototype. Finally, in Section 4, the conclusion and future research are detailed.

2. Dynamic Model of the Wind Turbine with a Rotatory Tower

This section is divided into four parts: the first is the description of the mechanic model, the second is the description of the aerodynamic model, the third is the description of the electric model, and, finally, the fourth is the combination of the aforementioned models to obtain the final dynamic model.

The dynamic model of the wind turbine is, first, the equations that represent the change between the wind energy and mechanic energy and, second, the equations that represent the change between the mechanic energy and electric energy.

2.1. The Mechanic Model

A windward wind turbine of three blades with a rotatory tower is considered. First, the Euler Lagrangian method [11, 12] is used to obtain the model that represents the change from the wind energy to mechanic energy for the wind turbine and the change from the electric energy to mechanic energy for the tower. The masses are concentrated at the center of mass. Consider the lateral view of Figure 1 and upper view of Figure 2.

Figure 1: Lateral view of the wind turbine.

Figure 2: Upper view of the wind turbine.

From Figures 1 and 2, it can be seen that

where , , , , is the angular position of the tower motor in rad, is the angular position of a wind turbine blade in rad, and is the length to the center of the wind turbine blade in . Consequently, the kinetic energy , and potential energy , are given as

where is the tower mass in kg, is the blade mass in kg, is the acceleration gravity in m/s², is the constant length of the tower in m, and is the length to the center of the tower in m. The torques and are given as

where is the torque of the generator moved by the blade in kg m² rad/s², is the torque of the motor used to move the tower in kg m² rad/s², and are the spring effects presented when the blade is near to a stop in kg m²/s², and and are the shock absorbers in kg m² rad/s. Using the Euler Lagrange method in [11, 12] gives the following equations:

The angular position of the blade 1 is related to the angular position of the blades 2 and 3 as follows:

where and are the angular positions for the blades 2 and 3, respectively. Then, using (4), it yields the equations for blades 2 and 3 as a function of as follows:

where , , , and are the torques applied to move the blades 2 and 3, respectively. Adding , , and , respectively, it gives

Now, adding the three equations of (4), (6), and (7) and using (8), (9), and (10), it gives

where and . From Figures 1 and 2, of (11) is defined as follows:

where and and are the force of the air received by the three blades. Equation (12) describes the assumption that the air goes in one direction; if , then the maximum air intake moves the blades of the wind turbine, but if the tower turns to the left or to the right and changes, then the wind turbine turns and the air intake decreases. From Figures 1 and 2 and [12], of (11) is defined as follows:

where is a motor magnetic flux constant of the tower in Wb and is the motor armature current of the tower in A. Equations (11), (12), and (13) are the main equations of the mechanic model that represents the wind turbine and tower.

2.2. Aerodynamic Model

The aerodynamic model is the dynamic model of the torque applied to the blades. The mechanic power captured by the wind turbine is given by [2–5]

where is the air density in Kg/m³, is the area swept by the rotor blades in m² with radius in m, is the wind speed in m/s, and is the performance coefficient of the wind turbine, whose value is a function of the tip speed ratio , defined as [2–5]

For the purpose of simulation, the following model of is presented [2–5]:

where

and the coefficients are , , , , , , and is the blade pitch angle in rad. Using (12), (14), (15), (16), and (17) gives the dynamic model of the torque applied to the wind turbine blades as follows:

where is defined in (12), is defined in (15), is defined in (16), and is defined in (17).

2.3. The Electric Model

Now, analyze the change from the mechanic to electric energy for the wind turbine generator and the electric energy to mechanic energy for the tower motor. Figure 3 shows the wind turbine generator and Figure 4 shows the tower motor.

Figure 3: The wind turbine generator.

Figure 4: The tower motor.

In [12, 13], they presented the model of a motor, and similarly, from Figure 3, a model for the generator is obtained by using the Kirchhoff voltage law. Therefore, the dynamic models of the motor and generator are as follows:

where is the motor back emf constant in V s/rad, is the generator back emf constant in V s/rad, is the motor armature resistance in , is the generator armature resistance in , is the motor armature inductance in H, is the generator armature inductance in H, is the motor armature voltage in V, is the generator armature voltage in V, is the motor armature current in A, and is the generator armature current in A. For the generator of this paper, . Thus, (19) becomes

2.4. The Final Dynamic Model

Thus, (11), (13), and (18) that represent the change between the wind energy and mechanic energy and (20) that represents the change between the mechanic energy and electric energy are considered as the wind turbine dynamic model with a rotatory tower.

Define the state variables as , , , , , , the inputs as and , and the output as . Consequently, the dynamic model of (11), (13), (15), (18), and (20) becomes

where is defined in (16) and is defined in (17).

Remark 1. The model of this research considers a direct current motor and a direct current generator. This wind turbine is proposed to be used on the house roof for the daily usage to feed a light as a renewable source. If more than one wind turbine is used, more electricity is generated. The proposed model could be extended to other kinds of machines by changing (20).

Remark 2. Note that the dynamic model states of the wind turbine with a rotatory tower (21) is a linear system [6–10].

3. Experiments to Validate the Dynamic Model

Figure 5 shows the prototype of a wind turbine with a rotatory tower which is considered for the simulations of the dynamic model. This prototype has three blades with a rotatory tower which does not use a gear box. Table 1 shows the parameters of the prototype. The parameters and are obtained from the wind turbine blades. The parameters , , and are obtained from the tower motor. The parameters , , , and are obtained from the wind turbine generator. The parameters , , , and are obtained from [2–5].

Table 1: Parameters of the prototype.

Parameter
Value

0.5 m

0.5 kg

0.45 Vs/rad

0.5 rad

Figure 5: Prototype of a wind turbine with rotatory tower.

The dynamic model of the wind turbine with a rotatory tower is given by (21) with the parameters of Table 1. are considered as the initial conditions for the plant states , , , , , and . The root mean square error (RMSE) is used, and it is given as [13–15]

where , , or and , , and are the real data of , , and , respectively, , .

3.1. Example 1: The First Behavior

Figures 6 and 7 show the wind turbine inputs, outputs, and states. Table 2 shows the RMSE for the errors.

Table 2: The RMSE for Example 1.

RMSE

Figure 6: Inputs and output of the wind turbine of Example 1.

Figure 7: States of the wind turbine of Example 1.

From Figures 6 and 7, it can be seen that the dynamic model has good behavior described as follows: (1) from  s to  s, both inputs are fed; consequently, the tower moves far of the maximum air intake, the generator current is decreased, and the wind turbine blades stop moving; (2) from  s to  s, both inputs are not fed; consequently, current is not generated and both the tower and wind turbine blades do not move; (3) from  s to  s, both inputs are fed, but the air intake is positive and tower voltage is negative; consequently, the tower returns to the maximum air intake, the generator current is increased, and the wind turbine blades move; (4) from  s to  s, both inputs are not fed; consequently, current is not generated and the tower and wind turbine blades do not move. The dynamic model is a good approximation of the wind turbine behavior because the signals of the first are near to the signals of the second. From Table 2, it is shown that the dynamic model is a good approximation of the real process because the RMSE is near to zero.

3.2. Example 2: The Second Behavior

Figures 8 and 9 show the wind turbine inputs, outputs, and states. Table 3 shows the RMSE for the errors.

Table 3: The RMSE for Example 2.

RMSE

Figure 8: Inputs and output of the wind turbine of Example 2.

Figure 9: States of the wind turbine of Example 2.

From Figures 8 and 9, it can be seen that the dynamic model has good behavior described as follows: (1) from  s to  s, the input air is fed and the tower input is not fed; consequently, the tower remains in the maximum air intake, the generator current is maximum, and the wind turbine blades have motion; (2) from  s to  s, the air is not fed and the tower input is fed; consequently, current is not generated, the tower moves far of the maximum air intake, and the wind turbine blades do not have motion; (3) from  s to  s, the air is fed and the tower input is not fed; consequently, the tower does not move, the generator current is minimum, and the wind turbine blades almost do not move; (4) from  s to  s, the air is not fed and the tower input is fed with a negative voltage; consequently, current is not generated, the tower returns to the maximum air intake, and the wind turbine blades do not have motion. The dynamic model is a good approximation of the wind turbine behavior because the signals of the first are near to the signals of the second. From Table 2, it is shown that the dynamic model is a good approximation of the real process because the RMSE is near to zero.

4. Conclusion

In this paper, a dynamic model of a wind turbine with a rotatory tower was introduced; the simulations using the parameters of a prototype showed that the proposed dynamic model has an acceptable approximation of the wind turbine behavior. The proposed dynamic model could be used on control, prediction, or classification. As a future research, the modeling will be improved using other interesting methods as the least squares in [16–26], neural networks in [14, 27–32], or fuzzy systems in [15, 33–35].

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors are grateful to the editor and the reviewers for their valuable comments and insightful suggestions, which helped to improve this paper significantly. The authors thank the Secretaría de Investigación y Posgrado, Comisión de Operación y Fomento de Actividades Académicas del IPN, and Consejo Nacional de Ciencia y Tecnología for their help in this paper.

References

L.-L. Fan and Y.-D. Song, “Neuro-adaptive model-reference fault-tolerant control with application to wind turbines,” IET Control Theory & Applications, vol. 6, no. 4, pp. 475–486, 2012.

A. Žertek, G. Verbič, and M. Pantoš, “Optimised control approach for frequency-control contribution of variable speed wind turbines,” IET Renewable Power Generation, vol. 6, no. 1, pp. 17–23, 2012.

E. B. Muhando, T. Senjyu, A. Yona, H. Kinjo, and T. Funabashi, “Disturbance rejection by dual pitch control and self-tuning regulator for wind turbine generator parametric uncertainty compensation,” IET Control Theory and Applications, vol. 1, no. 5, pp. 1431–1440, 2007.

C. Y. Tang, Y. Guo, and J. N. Jiang, “Nonlinear dual-mode control of variable-speed wind turbines with doubly fed induction generators,” IEEE Transactions on Control Systems Technology, vol. 19, no. 4, pp. 744–756, 2011.

R. Vepa, “Nonlinear, optimal control of a wind turbine generator,” IEEE Transactions on Energy Conversion, vol. 26, no. 2, pp. 468–478, 2011.

M. Dehghan and M. Hajarian, “Modied AOR iterative methods to solve linear systems,” Journal of Vibration and Control, 2012.

M. Hajarian, “Matrix form of the Bi-CGSTAB method for solving the coupled Sylvester matrix equations,” IET Control Theory & Applications, vol. 7, no. 14, pp. 1828–1833, 2013.

M. Hajarian, “Effcient iterative solutions to general coupled matrix equations,” International Journal of Automation and Computing, vol. 10, no. 5, pp. 481–486, 2013.

M. Hajarian, “Matrix iterative methods for solving the Sylvester-transpose and periodic Sylvester matrix equations,” Journal of the Franklin Institute, vol. 350, no. 10, pp. 3328–3341, 2013.

J. J. Rubio, G. Ordaz, M. Jiménez-Lizárraga, and R. I. Cabrera, “General solution of the Navier-Stokes equation to describe the dynamics of a homogeneous viscous fluid in an open tube,” Revista Mexicana de Física, vol. 59, no. 3, pp. 217–223, 2013.

F. L. Lewis, D. M. Dawson, and C. T. Abdallah, Control of Robot Manipulators, Theory and Practice, 2004.

M. W. Spong, S. Hutchinson, and M. Vidyasagar, Robot Modeling and Control, John Wiley & Sons, New York, NY, USA, 2006.

J. J. Rubio, J. Pacheco, J. H. Pérez-Cruz, and F. Torres, “Mathematical model with sensor and actuator for a transelevator,” Neural Computing and Applications, vol. 24, no. 2, pp. 277–285, 2014.

J. J. Rubio and J. H. Perez-Cruz, “Evolving intelligent system for the modelling of nonlinear systemswith dead-zone input,” Applied Soft Computing B, vol. 14, pp. 289–304, 2014.

J. J. Rubio, “Evolving intelligent algorithms for the modelling of brain and eye signals,” Applied Soft Computing B, vol. 14, pp. 259–268, 2014.

S. Ding, R. Ding, and E. Yang, “A filtering based recursive least squares estimation algorithm for pseudo-linear auto-regressive systems,” Journal of the Franklin Institute, 2013.

F. Ding, “Combined state and least squares parameter estimation algorithms for dynamic systems,” Applied Mathematical Modelling, vol. 38, no. 1, pp. 403–412, 2014.

F. Ding, “Coupled-least-squares identification for multivariable systems,” IET Control Theory & Applications, vol. 7, no. 1, pp. 68–79, 2013.

F. Ding, “Decomposition based fast least squares algorithm for output error systems,” Signal Processing, vol. 93, no. 5, pp. 1235–1242, 2013.

F. Ding, “Two-stage least squares based iterative estimation algorithm for CARARMA system modeling,” Applied Mathematical Modelling, vol. 37, no. 7, pp. 4798–4808, 2013.

L. Han, F. Wu, J. Sheng, and F. Ding, “Two recursive least squares parameter estimation algorithms for multirate multiple-input systems by using the auxiliary model,” Mathematics and Computers in Simulation, vol. 82, no. 5, pp. 777–789, 2012.

H. Hu and F. Ding, “An iterative least squares estimation algorithm for controlled moving average systems based on matrix decomposition,” Applied Mathematics Letters, vol. 25, no. 12, pp. 2332–2338, 2012.

X. L. Li, L. C. Zhou, R. F. Ding, and J. Sheng, “Recursive least-squares estimation for Hammerstein nonlinear systems with nonuniform sampling,” Mathematical Problems in Engineering, vol. 2013, Article ID 240929, 8 pages, 2013.

S. Pan and J.-S. Chen, “A least-square semismooth Newton method for the second-order cone complementarity problem,” Optimization Methods & Software, vol. 26, no. 1, pp. 1–22, 2011.

D. Q. Wang, F. Ding, and Y. Y. Chu, “Data filtering based recursive least squares algorithm for Hammerstein systems using the key-term separation principle,” Information Sciences, vol. 222, pp. 203–212, 2013.

D. Wang, F. Ding, and L. Ximei, “Least squares algorithm for an input nonlinear system with a dynamic subspace state space model,” Nonlinear Dynamics, vol. 75, no. 1-2, pp. 49–61, 2014.

J. Chen and R. F. Ding, “An auxiliary-model-based stochastic gradient algorithm for dual-rate sampled-data Box-Jenkins systems,” Circuits, Systems, and Signal Processing, vol. 32, no. 5, pp. 2475–2485, 2013.

J. Chen, “Several gradient parameter estimation algorithms for dual-rate sampled systems,” Journal of the Franklin Institute, vol. 351, no. 1, pp. 543–554, 2014.

F. Ding, X. Liu, H. Chen, and G. Yao, “Hierarchical gradient based and hierarchical least squares based iterative parameter identification for CARARMA systems,” Signal Processing, vol. 97, pp. 31–39, 2014.

M. Pratama, S. G. Anavatti, P. P. Angelov, and E. Lughofer, “PANFIS: a novel incremental learning machine,” IEEE Transactions on Neural Networks and Learning Systems, vol. 25, no. 1, pp. 55–68, 2014.

J. H. Pérez-Cruz, J. J. Rubio, J. Pacheco, and E. Soriano, “State estimation in MIMO non- linear systems subject to unknown deadzones using recurrent neural networks,” Neural Computing and Applications, 2013.

J. Sun, J.-S. Chen, and C.-H. Ko, “Neural networks for solving second-order cone constrained variational inequality problem,” Computational Optimization and Applications, vol. 51, no. 2, pp. 623–648, 2012.

A. Buchachia, “Dynamic clustering,” Evolving Systems, vol. 3, no. 3, pp. 133–134, 2012.

E. Lughofer, “Sigle pass active learning with conflict and ignorance,” Evolving Systems, vol. 3, no. 4, pp. 251–271, 2012.

L. Maciel, A. Lemos, F. Gomide, and R. Ballini, “Evolving fuzzy systems for pricing fixed income options,” Evolving Systems, vol. 3, no. 1, pp. 5–18, 2012.