Open access peer-reviewed chapter

Robust Nonlinear Control Strategy for Small Wind Turbines: A Case Study

By Ridha Cheikh and Hocine Belmili

Submitted: January 28th 2019Reviewed: December 3rd 2019Published: April 15th 2020

DOI: 10.5772/intechopen.90750

Downloaded: 43

Abstract

This chapter presents a case study of robust nonlinear control strategy using nonlinear feedback control technique based on Lyapunov theory and associated with robust control laws. The proposed approach aims to enhance robustness of the wind turbine control scheme. In fact, we selected as a case study, the most used electrical generator in small-scale wind applications, the Permanent Magnet Synchronous Generator (PMSM). Indeed, the control strategy presented in this chapter allows an efficient operation of the wind turbine in the standalone operating mode, offers a nonlinear handling of the WECS(s) and guarantees maximum wind power harvesting and robustness against critical working conditions. Talking about stability, in several wind generator control schemes; a such classical PI controllers-based scheme can easily be disturbed by any uncertainty of the system parameters, thus, in this chapter, we focused on how to overcome this issue by proposing a robust control strategy based on nonlinear controller derived from the Lyapunov Theory. The chapter presents numerical simulations within Matlab/SIMULINK environment. These results proved the effectiveness and the benefits of the proposed approach.

Keywords

  • WECS
  • PMSG
  • eigenvalues
  • nonlinear feedback control
  • Lyapunov theory
  • robust control

1. Introduction

During the last decade, the world has known a significant increase in electricity production from renewable energy sources, due in part to the liberalization of the electricity market, which has attracted the greed of new producers who are oriented towards cheaper productions (cogeneration, biomass, wind, etc.) and not in overly heavy investments (thermal power plants), and, secondly, because the emerging ecological awareness on climate change due to the emission of greenhouse gases has resulted in a political will that results from international protocols, such as Kyoto 1997 and Paris 2015, to encourage the use of clean and renewable energies for electricity generation [1].

Currently, wind energy is an important part of renewable energy production. Based on statistics from the Global Wind Energy Council (GWEC), more than 50 GW of clean, emission-free wind power was installed in 2018, bringing total installations to 591 GW globally; those statistics forecast that the global installed capacity can reach more than 817 GW by 2021 [2]. However, as the wind energy is very different from conventional sources, it is an intermittent energy source that is not continuously available for conversion into electricity and outside direct control because the used primary energy cannot be stored, it may be predictable but cannot be dispatched to meet the demand of an electric power system. Thus, the rapid growth of wind energy penetration into power systems causes many problems regarding the power flow control flexibility, which has forced many countries to revise their grid codes to ensure stable and reliable network operation. Power systems can reach very high wind penetration levels, as is the case of Denmark (42.4%), Germany (14.8%) and Spain (14.1%) [3, 4].

Nowadays the major challenge of the wind energy industry is to respect the deferent electrical grid requirements of each country; this means that power systems doest accept except a pure and balanced electrical energy comes from wind farms and in case of grid disturbances; those wind farms have to remain connected and support the grid to avoid black outs [5, 6]. As a result, the successful installation of wind turbines subject to requirements imposed by any country’s energy market will not be possible without the contribution of advanced and robust control strategies [7].

Currently, wind turbine technology is almost based on two topologies: fixed speed-based topology (FST) which uses Induction Machines (IMs) and variable speed-based topology (VST) which generally uses Double Fed Induction Machines (DFIMs) or Permanent Magnet Synchronous Machines (PMSMs) [8].

However, the FST topology is recognized by many advantages such as costless and no complexity even more, nevertheless it has many disadvantages such as lower efficiency; power fluctuations in the grid due to wind gusts, short life cycle due to significant efforts undergone by its structure [9]. Thus, even the wind system complexity will increase, but VST, in which power fluctuations could be reduced and hence wind gusts can be stored as kinetic energy in large rotating masses, has been proposed to mitigate all drawbacks of FST.

In recent years, choosing the type of the electrical generator of WTs is a serious subject of several research works. This component represents the heart of a WT; it is used to produce electric power with lower cost and under variable nature of the wind speed. In the wind industry, two types of generators dominate the sector: (a) the Permanent Magnet Synchronous Generator (PMSG) with very low inertia, high volumetric torque, good efficiency and better controllability, which, furthermore, with the possibility to direct drive train, has become a serious competitor of (b) the Double Fed Induction Machines (DFIG). This latest (DFIG) is the most used in wind industry sector due to the ability to control powers flow with significantly improved yield. The DFIG can operate in different speed modes (sub-synchronous, synchronous and super-synchronous). Due to the progress of modern technologies, the use of the DFIG in large power scales seems to have become a more efficient solution for electrical energy generation whether for onshore or offshore installations [10, 11, 12].

Therefore, developing robust nonlinear control algorithms for electro-mechanical systems actuated by different types of electric machinery, emphasizing system stability and robustness, is a very essential issue [13]. Under this scope, this chapter proposes a robust nonlinear control strategy based on Lyapunov Theory of a standalone Permanent Magnet Synchronous Generator (PMSM). Therefore, the chapter is focused on the design of nonlinear control strategy using the nonlinear feedback control techniques and Lyapunov theory to guarantee robustness of the wind energy conversion systems.

Finally, based on simulation results, a general conclusion is presented in this chapter showing the performance of the proposed control strategy used for the studied wind turbine.

2. Lyapunov theory and robust control association

2.1 Lyapunov theory fundamentals

The theory introduced in the late 19th century by the Russian mathematician Alexandr Mikhailovich Lyapunov is generally the most useful approach for analyzing the stability and designing control of linear or nonlinear systems. In fact, Lyapunov’s work for stability problems includes two methods: the so-called linearization method and the direct method. The linearization method draws conclusions about a nonlinear system’s local stability around an equilibrium point from the stability properties of its linear approximation [14]. However, the direct method (or second method) is not restricted to local motion, and determines the stability properties of a nonlinear system by constructing a scalar “energy-like” function for the system and examining the function’s time variation [15]. In other words, the direct method is a generalization of the idea that if there is some “measure of energy” in a system, then we can study the rate of change of the system’s energy to ascertain stability [16]. In this chapter, we will focus on the direct method for the control scheme design.

2.1.1 Lyapunov’s direct method concept

Lyapunov’s direct method for stability is now referred to as the Lyapunov stability criterion and makes use of a scalar “energy-like” function V(x), which has an analogy to the potential function of classical dynamics.

It can be introduced for a system ẋ=fxhaving a point of equilibrium at x=0. Consider a function Vx: RnRsuch that

  • Vx=0if and only if x=0. Vx>0if and only if x0.

  • V̇x=ddtVx=i=1nVxifix=V.fx0for all values of x0.

For asymptotic stability, the condition V̇x<0for x0is required.

Then Vxis called a Lyapunov function and the system is stable in the sense of Lyapunov.

2.1.2 Stability of a linear system

Consider a linear system described by its state space model as follows:

Ẋ=AX(A is a finite matrix), we can say that A is asymptotically (or exponentially) stable if all its eigenvalues’ real parts are negative; in other words, all the system states will converge to the vicinity of their equilibrium points [16].

The above stability condition can be verified equivalently through Lyapunov concept, using the quadratic Lyapunov function as follows:

V=XtPXwith P is a symmetric positive definite matrix,

  • If X0XtPX>0

  • then V̇=ẊtPX+XtPẊ. Replacing the value of Ẋwe have the following:

V̇=XtAtP+PAX=XtQX,so,AtP+PA=Q

The stability condition is that V̇<0so the matrix Qshould be positive definite for some definite matrix P=Pt[17].

2.2 Modeling of parametric uncertainty for AC machines

2.2.1 Resistance variation

Heating of the machine windings due to temperature increase causes an increase in stator and rotor resistances. This variation can be generally presented in the following relation:

R=R01+αT

where αis the temperature coefficient value (≈4.1 × 10−3) and T is the temperature in degree Celsius.

2.2.2 Inductance variation

The variation of the inductances is caused by the saturation of the ferromagnetic materials of the machine. The dynamics of this variation and that of currents cannot be considered slow [18].

2.2.3 Skin effect

This concerns the induction machine; in fact, the rotor bars’ resistance increases with the value of the frequency of rotor currents. An increase in the frequency causes a backflow of the currents to the periphery of the bars. The variation of the resistance of the bars is calculated numerically using a finite element method [18].

In steady-state operation, that is to say in the absence of parametric disturbances, the machine parameters take nominal values.

In general, the machine parameters are subject to unknown variations resulting from the different situations in which it evolves. The variation of the resistances is directly related to the fluctuations of the temperature of the machine; whereas the variations of the inductances are related to the phenomenon of magnetic saturation. This type of disturbance is called structured uncertainty.

From above, let us express the electrical parametric variations of the machine, such as:

Rs,rRs,r0+ΔRs,r,Ls,rLs,r0+ΔLs,rσσ0+Δσ,TrTr0+ΔTr

2.2.4 Global uncertainty model

For an AC machine represented by the nonlinear system below:

ẋ=fx+gxuy=hx

To globally model uncertainty of the previous system, we can rewrite it as follows:

ẋ=fx+Δfx+gx+Δgxuy=hx

where the exact part of the plant is represented by the functions (f, g, h), while (∆f, ∆g) represent the uncertain part of the system [19].

A robust control of the uncertainty is not guaranteed by the linear control law. Thus, one should design a robust law based on the Lyapunov theory to achieve an input–output stability and robustness against uncertainty [19].

3. Case study: robust control of a standalone WECS-based PMSG

3.1 System modeling

The schema of the real studied system is presented in Figure 1(a). However, for simplification purposes, and due to the power electronics and local grid dynamics being much faster than other system dynamics, it will not be considered and the simplified system is presented in Figure 1(b) [14].

Figure 1.

Schema of the studied WECS-based PMSG: (a) real system; (b) simplified system.

From Figure 1(b), the control goal is to change the chopper equivalent resistance RL (load) at the generator terminals in order to adjust the generator speed to its optimum value, and then the maximum power capture will be achieved.

The global system model can be expressed part per part in the following equations:

vt=v¯+Δvt,Pt=ΓtωrΓt=12πρR3v2CΓλCΓλ=q2λ2+q1λ+q0λ=ωrRvCP=λ×CΓAeromechanicalpartE1
ddtid=Rs+RLLd+LLid+pLqLLLd+LLiqωhddtiq=Rs+RLLq+LLiqpLd+LLLq+LLidωh+pΦmLq+LLωhelectricalpartE2
Jhdωhdt=ηiΓtΓgωh=ωr×iΓg=pΦmiqElectromecanicalpartE3

where v(t) is the wind speed (average and turbulent components); ΓtΓgare respectively the turbine and generator torques; Ptis the mechanical power; ωrωhare turbine and generator speeds respectively; ρis the air density; Ris the swept surface radius; CΓλis the torque coefficient; q0q1q2are given in the appendix; CPis the power coefficient; λis the tip speed ratio (TSR); ηis the efficiency; iis the gearbox ration; id, iq, Ld, Lq are respectively the d/q axis currents and inductances; Rs is the stator resistance; p is the pair poles; RL and LL are respectively the resistance and inductance of the equivalent load.

From Eq. (1) to Eq. (3), choosing radius x=x1x2x3T=idiqωhTas state vector, u=u1u2=RLvas control signal and y=ωhas a desired output, the nonlinear state space model of the system can be defined as follows:

ẋ=1Ld+LLRsx1+pLqLLx2x31Lq+LLRsx2pLd+LLx1x3+pΦmx31Jhd3x32pΦmx2+1Ld+LLx1u11Lq+LLx2u21Jhd1u22+d2u2x3E4

where d1=ηρπR32iq0,d2=ηρπR42i2q1,d3=ηρπR52i3q2.

3.2 Lyapunov linearization method

In this subsection, we will apply the Lyapunov Linearization (first method) to the studied system (Eq. (4)) in order to deal with the system’s small-signal stability under wind speed variation.

Suppose that x¯u¯is an equilibrium point and input [20].

Using Taylor’s series expansion of the function f we have:

δx=xx¯=δx1δx2δx3Tδu=uu¯=δu1δu2TE5

The derivative variable within a unique function can be written as:

ẋ=FxtutE6
δ̇xtFx¯u¯+Fxx=x¯u=u¯δxt+Fux=x¯u=u¯δut+h.o.t0E7

where Fx¯u¯=0.

Then:

δ̇xtFxx=x¯u=u¯δxt+Fux=x¯u=u¯δutE8

The matrices, A, B of the LTI system are constant matrices:

AFxx=x¯u=u¯Rn×Rn,BFux=x¯u=u¯Rn×Rm

The linearized version around an equilibrium point characterized by the quintuple x¯u¯=i¯sdi¯sqω¯hR¯Lv¯can be written as follows:

δ̇x=a1+a2R¯La3x¯3a3x¯2b1x¯3b2+b3R¯Lb1x¯1+b40c4c2v¯+c3x¯3δx+a2x¯10b3x¯200c1v¯+c2x¯3δuE9

where,

a1=RsLd+LL;a2=1Ld+LL;a3=pLqLLLd+LL.b1=pLd+LLLq+LL;b2=RsLq+LL;b3=1Lq+LL;b4=pΦmLq+LL;c1=2d1Jh;c2=d2Jh;c3=2d3Jh;c4=pΦmJh.

3.2.1 Small-signal stability

In this section, the small-signal stability of the linearized system will be checked regarding wind speed variation; therefore, we will have a point of view about the impact of the wind speed variation on the system stability. However, in open loop, without control, the system is excited by a wide range of wind speed variation starting from the cut-in to cut-out wind speed. The eigenvalues of the matrix A (Jacobian) are calculated and presented in Table 1.

v¯(m/s)45671012
λ 1,2−8.56±15i−7.09±18i−6.29±20i−5.74±23i−5.07±32i−4.95±37i
λ 3−279−342−406−470−666−798

Table 1.

Linearized system eigenvalues under wind speed variation.

Figure 2 shows clearly that for all values of wind speed, all the eigenvalues still have negative real parts and take places in the left-half of the complex plane, which means the system is almost stable under small perturbation. However, with wind speed increase, the complex-conjugated eigenvalues start immigrating almost close to the imaginary axis and then exhibit bad damping characteristics (Figure 3). Hence, damping should be improved in closed loop-based linear controller.

Figure 2.

Eigenvalues (poles) immigrating under wind speed variation.

Figure 3.

System characteristics regarding damping/overshoot for two wind speed cases (4 m/s left, 12 m/s write).

3.3 System nonlinear feedback linearization

If we assume the wind speed dynamics to be very slow compared to the other system dynamics, the nonlinear system of Eq. (4) can be reformulated as follows:

ẋ=RsLd+LLx1+pLqLLLd+LLx2x3RsLq+LLx2pLd+LLLq+LLx1x3+pΦmLq+LLx31Jhd1v2+d2vx3+d3x32pΦmx2fx+1Ld+LLx11Lq+LLx2,0gxuy=hx=x3E10

where d1=ηρπR32iq0,d2=ηρπR42i2q1,d3=ηρπR52i3q2.

The system of Eq. (10) has smooth nonlinearities; therefore, the Feedback Linearization Technique can be applied to control the system.

From Eq. (10), taking derivatives of output y with respect to time up to r times will give:

yr=Lfrhx+LgLfr1hxu0E11

r is the relative degree (r < n).

Lfrhxis the Lie derivative of h(x) along the direction of the vector field f(x) up to r times, LgLfr1hxis the Lie derivative of h(x) along the direction of the vector field g(x).

If we haveαx=Lfrhxβx=LgLfr1hxE12

Eq. (11) becomes:

yr=αx+βxuE13

Then the feedback linearization control can be chosen as follows:

ux=1βxαx+υE14
yr=υE15

The relative degree of Eq. (20) is r = 2, it is a linear input–output double-integrator, using linear methods to guarantee the system control (tracking problem).

Note that, the field vectors αxand βxshould be completely known to apply the feedback linearization technique [21].

The computation of the field vectors αxand βxgives:

αx=pΦmJhb2x2+b1x1x3+b4x3+1Jhd2iv+2d3x3ẋ3βx=b4Jhx2E16

To harvest the maximum wind power, the generator should be running at optimal speed for each wind speed, thus the generator speed reference can be expressed as follows:

yrt=ωh=λivtRE17

where λis the optimal tip speed ratio.

Let us take an output tracking error as follows:

et=yrtytE18

and introduce the following error vector:

ξt=etėtTE19

Hence, the control objective can be achieved using the ideal control law:

ux=1βxαx+yr2tkTξE20

The polynomial s2+k1s+k2would have all its roots in the left-half of the complex plane if we choose an appropriate gain matrix k=k2k1T. This means that the tracking error asymptotically converges to zero (e¨t+k1ėt+k2et=0). Thus, the system poles are placed following the specifications concerning the overshoots and the settling times.

3.4 Robust control law design

As we mentioned previously, AC machines are subject to a large model uncertainty due to parameter variations, noises, measurement errors etc. In this chapter, we will take into account just electrical and mechanical parametric uncertainties. In the studied case, it is clear that a successful control strategy of the WECS-based PMSG hinges on the good computing of control u, which is strongly dependent on the machine parameters (Eq. (20)). Therefore, we propose a novel robust control law to handle the uncertainty issue and offer an accurate description of the system model [19].

From Section 2.2.1, Eq. (10) can be rewritten as follows:

ẋ=fx+Δfx+gx+Δgxuy=hxE21

However, the exact parts of the model are f(x), g(x) and h(x), while Δfand Δgpresent the uncertainty.

Consequently, we have:

y¨=υ+ΔυΔυ=Δαx+ΔβxuE22

Δυis the effect of uncertainty, it is an unknown term, which makes the linear feedback control law of Eq. (15) incapable of providing a robust tracking performance of the system. Thus, the following robust control law is proposed to handle this issue.

From Eq. (22), we rewrite the error dynamics as follows:

ξ̇=01k2k1Acξ+01Bυ+ΔυE23

We chose the Lyapunov function as follows:

V=12ξTE24

The derivation of V gives

V̇=12ξ̇+ξTPξ̇E25

P is a positive semi-definite matrix in which

AcTP+PAc=QE26

Q is an identity matrix.

After computation, V̇becomes

V̇=12ξT+ξTPBυ+ΔυE27

Vsatisfies the condition V̇<0along the solution trajectory of the system, which gives the following:

ξTPBυ+Δυ<0E28

υis chosen as

υ=FsgnξTPBE29

with F>Δυ.

The substitution of Eq. (29) in Eq. (28) gives

ξTPBFsgnξTPB+Δυ<0E30

Eq. (30) is always negative for any value of ξTPBprovided F>Δυ. The value of the gain F is chosen accordingly by trial and error.

We have

sgnξTPB=1ifξTPB>00ifξTPB=01ifξTPB<0E31

After computation, the P matrix element can be set as follows:

P=K12K2+1+K22K112K212K21+K22K1K2E32

Finally, the linear control signal is expressed as

υ=yr2tkTξFsgn12K21+K22K1K2ξE33

4. Simulation results and discussion

Based on the control scheme shown in Figure 4, we have simulated within Matlab/Simulink® the studied system shown in Figure 1. The system data are given in the appendix. Through simulation, we have considered that all system variable states (x 1, x 2, x 3) are available for feedback. The simulation aims to check two control performances, a wind maximum power tracking performance under a realistic wind speed, and robustness performance against both rapid wind speed variation (sharp) and parametric uncertainty of the generator inertia.

Figure 4.

Robust control simulation scheme.

4.1 Tracking performance check

For tracking performance under wind speed fluctuations (Figure 5(a)–(d)), the power coefficient Cp holds easily its maximum value (Cpmax ) and the same for tip speed ratio, TSR (λ). Consequently, the PMSG runs at optimal speeds for each wind speed value (see Figure 5(d)) and the maximum power extraction is indeed guaranteed (Figure 6(b)).

Figure 5.

Check of tracking performance under a realistic wind speed. (a) Wind speed; (b) TSR; (c) power coefficient; (d) generator speed.

Figure 6.

(a) Control signal; (b) turbine power; (c) optimal points neighborhood; (d) electromagnetic torque versus TSR.

However, from Figure 6(a), it should be mentioned that the control signal, which represents the equivalent chopper resistance (R L), has the chattering phenomenon due to discontinuous control effect (see Eq. (33)). From Figure 6(c) and (d), all system operating points are in the neighborhood of the optimal operating points.

Then, the next simulation is to check the performance of the proposed controller facing critical wind speed conditions. Hence, we have subjected the system to a sharp rise and drop variation of the wind speed (Figure 7(a)). The simulation results show that even under abrupt changes of the wind speed, the system kept its stability and continued working after some short transitory disturbances and small static errors, which can be observed through Figure 7(a)-(c) and Figure 8(a)-(c).

Figure 7.

Robustness check under sharp wind speed variation. (a) Wind speed profile; (b) wind turbine electrical power.

Figure 8.

(a) Power coefficient; (b) tip speed ratio; (c) generator speed.

4.2 Parametric uncertainty robustness check

Following Eq. (14) and Eq. (16), it is obvious that the control signal is calculated based on field vectors αxand βxwhose values are strongly dependent on the system parameters, especially the high speed shaft inertia J h. Thus, in order to check the proposed controller robustness against this parametric uncertainty, we have simulated the system with sharp variation of high speed shaft inertia as shown in Figure 9(a). From the obtained simulation results (Figures 10 and 11), except small static errors and very small disturbances, this parametric uncertainty would not have any notable effect on the stability of the control scheme. However, the field vectors αxand βxdynamics, the control signal and the generator speed can be shown in and Figure 12.

Figure 9.

(a) High speed shaft inertia; (b) average wind speed.

Figure 10.

Power coefficient Cp.

Figure 11.

Tip speed ratio.

Figure 12.

(a) Control signal; (b) generator speed; (c) α(x); (d) β(x).

5. Conclusion

In this chapter, we have discussed the robust control technique based on Lyapunov theory to ensure tracking performance and robustness of a nonlinear system. Mainly, the chapter presents this control strategy based on a case study of a standalone nonlinear WECS-based PMSG. In fact, we focused on tracking performance to ensure maximum wind power extraction and robustness performance against parametric uncertainty. Whereas, the control strategy uses the state feedback linearization technique associated with a control law derived from the Lyapunov theory. The obtained results through simulation have proved the efficiency of the proposed control technique.

  1. Turbine rotor:

  2. R = 2.5 m ρ = 1.25 kg/m3, CPmax = 0.477, λ* = 7

  3. Drive train:

  4. η = 1, i = 7; Jh = 0.0552 kg m2

  5. Torque coefficient parameters:

  6. q 0 = 0.0061, q 1 = −0.0013, q 2 = −9.7477 × 10−4

  7. PMSG:

  8. p = 3, Rs = 3.3 Ω

  9. Ld = 0.0416 H

  10. Lq = 0.0416 H

  11. LL = 0.08 H

  12. Φ m = 0.4382 Wb

  13. Vs = 380 V

  14. Control parameters

  15. Linear control: k 1 = 100. k 2 = 4000

  16. Robust control: F = 5000.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Ridha Cheikh and Hocine Belmili (April 15th 2020). Robust Nonlinear Control Strategy for Small Wind Turbines: A Case Study, Design Optimization of Wind Energy Conversion Systems with Applications, Karam Y. Maalawi, IntechOpen, DOI: 10.5772/intechopen.90750. Available from:

chapter statistics

43total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Design Optimization of Wind Energy Conversion Systems with Applications

Edited by Karam Maalawi

Next chapter

Optimal Energy Production Planning in Power Systems with Wind Energy Sources

By Cristian-Dragos Dumitru and Adrian Gligor

Related Book

First chapter

Introductory Chapter: An Introduction to the Optimization of Composite Structures

By Karam Maalawi

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us