Pitch Controller Design of Wind Turbine Based on Nonlinear PI / PD Control

/is paper analyzed and discussed the nonlinear PI/PD control for wind energy conversion system (WECS) pitch controller design. /e nonlinear PI/PD controller consists of a classical PI/PD controller and a nonlinear gain table. In this paper, the detailed design procedures of the controller including not only the selection of control parameters but also the formulation of nonlinear gain table were given. In order to verify the effectiveness and correctness of the algorithm, a detailed 2MW simulation test platform based on MATLAB/SIMULINK environment was established. Meanwhile, some valuable conclusions were also obtained. /e presented contents of this article have the reference value and guiding significance for practical engineering application.


Introduction
With the global climate problem increasingly prominent and nonrenewable energy sources being largely consumed, the exploration and utilization of renewable energy sources such as wind energy and solar energy have been paid to a high attention.As the fastest growing renewable energy sources in the recent years [1][2][3][4], wind generation is the most prevalent in coastal regions spanning temperate and boreal climates.Countries such as China, USA, Denmark, and Canada possess significant wind power potential due to their high average wind velocities [5][6][7][8].Hence, the research and development about wind power have much practical value nowadays.Compared with the constant speed-constant frequency power system, the greatest advantage of variable speed-constant frequency wind power system widely used is the ability to get access to the maximum energy conversion.Mainstream variable speed-constant frequency power generation models include doubly-fed induction generator (DFIG) and permanent magnet synchronous generator (PMSG).Compared with DFIG, PMSG has many superior characteristics such as more efficient performance, higher reliability, and wider speed control range and is gradually becoming the first choice [9].erefore, the WECS based on PMSG was selected for research in this paper.
As Figure 1 shows, according to the wind speed, the working region of WECS typically can be divided into two regions, namely, partial-load region that has wind speed below the rated wind speed and full-load region that has wind speed above the rated wind speed.In the partial-load region, the control goal of the wind power system is generally to capture the maximum wind energy in order to achieve the maximum economic benefit.At this time, the partial-load region can also be called the maximum power point tracking (MPPT) region.In the full-load region, the variable propeller is particularly important.At this time, the generator-output power is limited at the rated value by the pitch control since the capacity of the generator and converter is limited [10,11].Furthermore, pitch control is also needed in some special working conditions such as the limited power control under rated wind speed, the low-voltage crossing in power grid, power grid needs to be injected into the inertia, and lighten the wind turbine loads.Meanwhile, it is pointed out that pitch control can also smoothen the power [12]. is paper focuses on the pitch control above the rated wind speed.
Figure 2 illustrates the pitch-controlled system.Obviously, Figure 2(a) depicts the blades action under pitchcontrolled, and Figure 2(b) shows the variable-pitch drive system in wind turbine hub.e variable-pitch drive system consists of pitch drive cabinet, variable-pitch drive, emergency power supply, pitch bearing and variable-pitch lubrication systems.e pitch control system and communication system are in the pitch drive cabinet.In general, the pitch controller and master controller are connected by the optical ber communication.
e PROFIBUS or CANopen protocols are generally used in this communication.Usually, the master controller sends the pitch instructions to the pitch controller and gets some important information from it at the same time.When the instructions are obtained, the pitch controller starts to control the variable-pitch drive-operating system.
At present, the research on variable-pitch control mainly includes two main categories: linear control and nonlinear control.
e research of the linear controller mainly has PI/PD control, H∞ controller with a linear matrix inequality approach [13], linear quadratic Gaussian LQR control [14], and generalized predictive control (GPC) [15].e nonlinear control mainly contains LPV control [16] and some intelligent control methods such as fuzzy control [17].Because wind turbine is a multi-input multi-output system containing strongly nonlinear dynamics, the e ect of nonlinear control is better than the linear control.Taking into account the complexity of the algorithm, the most applied controllers in engineering is still PI/PD control.erefore, the design of the variable propeller with the nonlinear PI/PD controller considering the nonlinear characteristics of the wind turbine has become the key.So, this article takes this as the focus of the study and summarizes some valuable conclusions.And this study can also provide the necessary reference and guidance for the follow-up engineering practice.e structure of this paper is as follows: the dynamic model and structure of WECS is brie y introduced in Section 2. Section 3 explains the controller design in details, including the controller parameter tuning and the calculation of nonlinear gain.In Section 4, the validation of the nonlinear PI/PD controller is surveyed using the MATLAB/Simulink simulator.Finally, Section 5 provides some useful conclusions.

Mathematical Model of WECS
WECS is mainly constitutive of wind turbine, transmission chain, permanent magnet synchronous generator, generatorside converter, DC bus, and the grid-side inverter.Its structure is shown in Figure 3.
where ρ is the air density, R is the radius of wind wheel, v is the wind speed, β is the pitch angle, ω t is the speed of wind turbine, λ is the tip speed ratio (TSR), and C P is defined as the wind energy conversion coefficient.A 2 MW baseline wind turbine-power coefficient is shown in Figure 4. Obviously, C P is the function of β and λ. e tip-speed ratio λ is given by 2.2.Dynamic Model of PMSG.If PMSG was supposed to be an ideal motor, it should meet the following criteria: (1) three-phase stators are symmetry.
(3) there is no damping winding on the rotor (4) magnetic saturation of iron core, vortex, and hysteresis loss can be neglected.(5) electronic conductivity of permanent magnet material is zero.
Under this assumption, the mathematical model of PMSG [23,24] is (3) Electromagnetic torque is given by where L d and L q are the d-axis and q-axis stator reluctance, R s is the stator resistance, Ψ is the permanent flux, n p is pole pairs, ω g is the speed of PMSG, i d and i q are d-axis and q-axis currents, respectively, and u d and u q are d-axis and q-axis voltages, respectively.Because nonsalient PMSG meets L d � L q � L, the mathematical model of PMSG has coupling terms such as n p ω gen i q , −n p ω gen (i d + Ψ/L), according to Equation (3). Given e mathematical model of PMSG could be rewritten as e above formula also can be rewritten into a transfer function G 0 (s).
Figure 5 shows the inner control of current.e closed loop transfer function of the current or torque can be obtained by Figure 5.And the closed loop transfer function also could be taken as a first-order inertial system.
where τ T is the inertial time constant and T g_ref is the command signal of T g .Equation (8) shows that the dynamic model of PMSG could be equivalent to a first-order inertial system by current feedforward decoupling control and first-order tuning for PI parameters.Usually, when τ T ≤ 0.01, we can leave out the inertial time delay of PMSG.

Dynamics of Blade Pitch Actuator.
ere are two kinds of drives in the blade pitch actuator at present.ey are the motor drive and hydraulic drive.Compared with the hydraulic drive, motor drive has the characteristics of lower cost and higher reliability and has been widely used.Actual motor drive is shown in Figure 6.
e collecting variable-pitch strategy was used.According to Figure 7, the blade dynamics is expressed as By leaving out the inertia time delay of the motor drive, the drive torque of motor is assumed as e blade pitch actuator dynamic can be given by Figure 8 where β ref is the command of β and the PD parameters of the blade pitch actuator are K βP and K βD .
If the delay coefficient meets the blade pitch actuator dynamic model can be simplified as

Shock and Vibration
Obviously, the inertia delay coe cient, τ β in Equations ( 12) and ( 13) can be reduced by increasing PD parameters.When K βP and K βD take the appropriate value, the inertia delay can also be left out.Actually, because the value of the parameters J Blade , μ, and f is unknown, the value of the control parameters K βP and K βD can only depend on experience.

Dynamic Model of Transmission Chain.
Figure 9 shows the generator cabin internal components of high-speed permanent magnet wind turbine.It is clear that its transmission chain mainly consists of spindle (or low-speed shaft), gearbox, and high-speed shaft.
e dynamics of the transmission chain is characterized by Figure 10.
According to Figure 10, the dynamics of the transmission chain is derived as e gearbox ratio is de ned as From Equations ( 14) and ( 15), the dynamic model of the transmission chain can be simpli ed as where   Shock and Vibration

Design of Controller
When the wind speed v is above the rated wind v rate , the system needs pitch control to unload the power and enable the speed ω t to approach the rated speed ω rate .If the inertia delay τ T and τ β are neglected by adjusting control parameters, the electromagnetic torque T g and pitch angle β can be approximated as their command values T g_ref and β ref .
where s is the Laplace factor.
By linearizing Equations ( 1) and ( 18), the small signal values are where Superscripts ∼ and -refer to the small signal value and mean value, respectively.o is high-order in nitesimal.Meanwhile, Equation ( 16) can also be rewritten as  Shock and Vibration 5 By eliminating the mean value of the upper equation and Equations ( 19) and ( 22) can be obtained.
Due to a < 0, k D can increase the inertia of the system.In general, k D and k I are not used at the same time because they have the opposite e ect on the system, namely, either PI control or PD control can be used in pitch control.

Design of Nonlinear PI Controller.
When k D 0, Equation ( 22) could be rewritten as By comparing d(s) with the typical second-order system, we get Equation ( 22) can be obtained.
Figure 11 shows the torque characteristics of WECS.Obviously, there are two equilibrium points A and B in the system.For equilibrium point B, if there is a perturbation increasing the rotor speed, then PMSG torque will be over wind turbine torque.After the perturbation is eliminated, rotor speed can decrease.By the same method, we can also con rm the system will be back to point B by the assumption that a little perturbation reduced the rotor speed.Hence, the working point B is stable.It could be con rmed that the equilibrium point A is unstable by this perturbation observation method.erefore, only the point B is a stable equilibrium point, and the range of stability is ω t ≥ ω A . e slope of wind turbine torque characteristic, T t is b and the slope of PMSG torque characteristic, nT g is −P rate /ω 2 t .When the system works at the point B, we can get Equation (26) by Figure 11.
Generally, the damping D eq is unknown.From Equation (26), the damping ratio ζ in Equation (25) meets Given (−ak P )/(2 −ak I J eq ) ε ≥ 1 and the damping ratio ζ and PI parameters will meet ζ > ε and erefore, the parameters K P and K I are given by K P 2J eq ω n ε, By the assumption of K P η I K I , Equation (29) could be rewritten as

Design of Nonlinear PD Controller.
Inertia has the ability to prevent frequency mutation in the system.Low frequency oscillation and rotor oscillation of the system can be suppressed by increasing inertia.erefore, PD controlproviding the virtual inertia for the system, is often applied to the small inertia systems and situations that need to be injected into inertia.Similarly, when k I 0, Equation ( 22) also could be rewritten as If D eq − b − P rate /ω 2 t in Equation ( 31) is neglected, the transfer function G vω is given by Given and K D η D J eq , the parameter K P is where BW is the bandwidth of G 0 (s).
In general, we keep K D and K P (or K I and K P ) xed and make −1/a as a nonlinear gain table.Meanwhile, the limit of the variable pitch rate and angle is taken into account, and pitch control based on nonlinear PI or PD controller is shown in Figure 12. .Obviously, all of the wind speed v, rotor speed of wind turbine ω t , and pitch angle β can have e ect on the nonlinear gain a by Equation (20).But, the sensitivity of the system to the pitch angle is far greater than other factors generally.When the system works at the rated working point ω t ω rate and T t T rate , the relationship between v and β is

Nonlinear Gain Table
In general, the nonlinear gain a only depends on the pitch angle β. e calculation process of nonlinear gain is shown in Figure 13, and the parameters of WECS are shown in Table 1.
In practice, the function C P (β, λ) can only be obtained by data tting, and these data for tting are generated by Bladed Software.In order to verify the correctness of the algorithm shown in Figure 8, an empirical formula of C P is used.
e truth value of a could be calculated by Equations ( 20) and (37).Table 2 shows the calculated and truth values of nonlinear gain at different wind speed.From Table 2, the calculated value of a is very close to the truth value of a. erefore, the algorithm shown in Figure 13 is valid.e minimum of β is 2.21 °above the rated wind, and the relationship between −1/a and β is shown in Figure 14

Validation of Algorithm by Simulation
In this section, simulations are carried out with MATLAB/Simulink environment.And the block diagram of WECS is shown in Figure 15.In order to verify the correctness and effectiveness of nonlinear PI and PD pitch control, a detailed 2 MW simulation test platform of WECS was constructed.e simulation parameters of the system are shown in Table 1, and the function C P (β, λ) is the empirical formula.e step-change wind shown in Figure 16(a) is used to simulate the actual wind condition.And the response of PI control is shown in Figure 17.Given η I � ε � 2 × 10 3 , the PI parameters K P and K I were calculated by Equation (30).Meanwhile, simulation results of PD control under stepchange wind shown in Figure 16 Figures 17(a) and 18(a) show the electromagnetic power curves.It is clear that the output power of WECS is maintained at 2 MW, and the power fluctuation with nonlinear gain control is obviously less than the one with fixed parameters control.e generator rotor speed curves are shown in Figures 17(b) and 18(b).Obviously, the speed steady-state errors of PI control is less than the speed steadystate errors of PD control.Figures 17(c) and 18(c) describe the pitch angle curve.Compared with the PI control, the actual value of the pitch angle under PD control failed to achieve the value in Table 2.However, the pitch control is generally used as a coarse control.erefore, the small pitch angle errors and rotor speed errors of PD control are acceptable.By Figures 17(d) and 18(d), the stator three-phase currents could be observed.e simulation results also indicate that the power of WECS is more sensitive to pitch control than the rotor speed, when the wind speed changes.
In Figure 19 2. Under the nonlinear PI/PD control, the output power uctuate is small.And the output power is very sensitive to the change in pitch angle under di erent wind conditions.Furthermore, the merits of nonlinear PI/PD control algorithms have been clearly revealed in Figures 19 and 20.

Conclusion
In this research, the design procedure of a pitch controller with nonlinear PI/PD control is given in detail, including the selection of PI/PD parameters and the calculation of Shock and Vibration nonlinear gain.After that, a simulation test platform for a two-megawatt WECS was built.Finally, the e ectiveness and correctness of the algorithm was veri ed by this test platform, and some useful conclusions were summed up: (i) e output power of WECS is the most sensitive to the change in the pitch angle.
(ii) Compared with the xed parameters control, the nonlinear control with nonlinear gain can make the power uctuate smaller.(iii) e precision of speed tracking is high, and the pitch angle can change into the right place under the PI control, compared with the PD control.

Figure 5 :
Figure 5: Inner control of current (remark: PI parameters are denoted as k iP and k iI ).(a) Current feedforward decoupling control.(b) First-order tuning for PI parameters.

Figure 10 :
Figure 10: Dynamic model of transmission chain (remark: Sun and Sate represent the sun wheel and planetary gear of gearbox).

Figure 13 :
Figure 13: e calculation process of nonlinear gain a (remark: Δβ is the precision value of β, and the initial value of Δβ is 0.01 °).
(b) are illustrated in Figure 18.Given η D � 0.7 and BW � 50, the PD parameters K P and K D were calculated by Equation (34).
(a), the wind speed decays at step change.Under this wind speed, the responses of WECS are shown in Figures 19(b)-19(f).Meanwhile, simulation results under gradual change wind speed shown in Figure 20(a) are given by Figures 20(b)-20(f).Di erent from other responses curves such power, torque, and rotor speed, the changes in the pitch angle are not smooth.is is because the rate of pitch angle change is limited at ±7 °/s.By thesis comparisons among classical PI control algorithm, classical PD control algorithm, nonlinear PI control algorithm, and nonlinear PD control algorithm, the similar conclusions could be summed up.Compared with the PD control, PI control has high tracking accuracy and the pitch angle has changed into the right place in Table

Figure 18 :
Figure 18: Responses under PD control.(a) Electromagnetic power curve.(b) Generator rotor speed curve.(c) Pitch angle curve.(d) Generator stator three-phase current under PI nonlinear gain.

Table 2 :
Nonlinear gain a.