Fault-Tolerant Neuro Adaptive Constrained Control of Wind Turbines for Power Regulation with Uncertain Wind Speed Variation

: This paper presents a novel adaptive fault-tolerant neural-based control design for wind turbines with an unknown dynamic and unknown wind speed. By utilizing the barrier Lyapunov function in the analysis of the Lyapunov direct method, the constrained behavior of the system is provided in which the rotor speed, its variation, and generated power remain in the desired bounds. In addition, input saturation is also considered in terms of smooth pitch actuator bounding. Furthermore, by utilizing a Nussbaum-type function in designing the control algorithm, the unpredictable wind speed variation is captured without requiring accurate wind speed measurement, observation, or estimation. Moreover, with the proposed adaptive analytic algorithms, together with the use of radial basis function neural networks, a robust, adaptive, and fault-tolerant control scheme is developed without the need for precise information about the wind turbine model nor the pitch actuator faults. Additionally, the computational cost of the resultant control law is reduced by utilizing a dynamic surface control technique. The effectiveness of the developed design is verified using theoretical analysis tools and illustrated by numerical simulations on a high-fidelity wind turbine benchmark model with different fault scenarios. Comparison of the achieved results to the ones that can be obtained via an available industrial controller shows the advantages of the proposed scheme.


Introduction
The key factor for the enhancement of the efficiency of a wind turbine is how to develop the control structure.Specifically, the pitch control design is a vital step in variable pitch wind turbines working in high wind speed, i.e., so-called full load region, to avoid hazardous operation as well as to avoid conservative power generation, i.e., less than nominal power [1].This objective is often known as power regulation for nominal power generation [2].Accordingly, the control design of wind turbines in power regulation has gained significant importance during the last decades [3].Several industrial controllers for power regulation use the Proportional-Integral-Derivative (PID) type control, as the linear controller [4,5].However, as wind turbines are complex nonlinear dynamic processes, linear controllers may not accurately render the expected performance [6].Consequently, in the last decade, modern and advanced controller schemes have been adopted to regulate power generation accurately, e.g., linear parameter varying control [7], gain scheduling [8], adaptive nonlinear control [9], optimal control [10], evolutionary algorithms [11], robust control [12], and fuzzy logic systems [13].A detailed review of power regulation controllers designed for wind turbines can be found in [14].
Wind turbine operation in the presence of high wind speed variation may lead to pitch actuator faults, which in turn leads to poor power regulation and catastrophic operation [4,15].These faults can be considered as pitch actuator bias, effectiveness loss, and dynamic change [2].Also, the debris build-up and blade erosion change the aerodynamic characteristics of the blades [16].The presence of faults can increase the need for maintenance operations and downtimes, which may lead to a decrease of the power generation and increase of the cost, particularly for offshore wind farms, installed in remote places that are sometimes difficult to reach [17,18].Thus, it is desirable to integrate the fault tolerance capability into the pitch actuator controller to attenuate the fault effects and keep the performance at the desired level, especially for large rotor and offshore deployments [7].Therefore, in the development of the wind turbine control structure, the fault-tolerant pitch controller design has been considered, and different approaches have been proposed, such as fuzzy control [19], adaptive sliding mode control [4], and robust linear parameter varying control [16].
In reality, the wind turbine's aerodynamic torque is a nonlinear function of wind speed [20].On the other hand, the pitch angle variation adjusts the speed via regulating the aerodynamic torque.However, the wind speed is a highly stochastic variable.Accordingly, it can be stated that the control function, from the pitch angle to the aerodynamic torque, is not completely known.This leads to the unknown control direction problem.Therefore, considering uncertain wind speed variation in the pitch angle control design of wind turbines is a significant challenge [21].In [22], a nonstandard extended Kalman filter is developed to estimate the wind speed for maximum power extraction of variable speed wind turbines.In [23], a comparative study using soft computing methodologies for the estimation of wind speed was presented.A review of the effective wind speed estimation-based control of wind turbines can be found in [24].Even though numerous methodologies have addressed wind speed estimation of wind turbines, the presented structures are still found to be very complicated and ineffective in practice.
On the other hand, in the full load region, if the wind turbine speed increases and violates the predefined limits, the mechanical brakes, located on the rotor, are engaged [15].This leads to generated power reduction considerably lower than the nominal one.Also, excessive rotor speeds in wind turbine operation may lead to a hazardous situation.Thus, to ensure safe operation of wind turbines, the rotor speed and its variation must be constrained within the safe-to-operate bounds.By that means, the variation of the generated power around the nominal power can be constrained at some predefined bounds.These bounds are designed within which the engagement of the mechanical brake is avoided.However, there are very limited works available in the literature that consider the constrained power generation.To fill this gap in the past literature, the authors recently developed a new strategy for constrained power generation [5], which can be viewed as an extension of the direct Lyapunov method to constrained systems.The core of this method consists of developing a barrier Lyapunov function (BLF) to constrain the generator speed and the generated power.An essential advantage of the BLF is that it guarantees that the corresponding arguments are constrained [25].The algorithm proposed by the authors in [5] uses the logarithm-type BLF for nonlinear wind turbines to constrain the generator speed and thus to generate the constrained power.However, in the authors' study, constrained performance was not guaranteed in the presence of faults.Also, the uncertain wind speed variation was not considered in designing the previous constrained control scheme.
Motivated by the aforementioned considerations, the primary objective of this paper was to design the pitch actuator control of wind turbines under uncertain wind speed variation to constrain the rotor speed and the generated power within the safe-to-operate bounds.These bounds were defined in order to avoid the engagement of the mechanical brake.In this manner, the overspeeding as well as the conservative power generation problems are resolved.The main idea consists of developing the BLF-based control to provide constrained behavior and further utilizing a Nussbaumtype function to cope with the unknown wind speed variation.The former is utilized to keep the generated power within the given desirable constraints, provided by the designer, and the latter is exploited to regulate the power without requiring accurate wind speed measurement, observation, or estimation.The controller is further developed to tolerate the pitch actuator faults for nonlinear wind turbine models.This considerably increases the reliability and efficiency as well.Thanks to the dynamic surface control (DSC) technique, the computational complexity of the control scheme is reduced by avoiding the repetitive differentiation of virtual control in the controller structure.In addition, the pitch actuator effort is smoothly bounded to avoid pitch actuator saturation.The second objective of the paper was to design a fault identification scheme to investigate the estimated fault signal, and then accurately generate fault information, which consists of fault detection, isolation, and type/size estimation.The fault identification task may be employed, for example, for planning effective and timely maintenance of offshore deployments, thus reducing the downtime of plants and the cost.Also, the availability is enhanced.In addition, to estimate the uncertain aerodynamic torque, a radial basis function (RBF) neural network is used, whose weights are automatically tuned without requiring any early training scheme.Finally, the numerical simulation is conducted to evaluate the proposed controller performance.Also, a comparison between the proposed controller and the available industrial PID controller is made considering the numerical control criteria, which are used to quantify the performance of both controllers.Accordingly, the contribution of this paper can be summarized as follows.
1.With the adoption of BLF-based constrained control, the rotor speed and its variation are constrained, and consequently, the variation of generated power around the nominal power will not violate the predefined constraint.This guarantees safe desirable nominal power generation, and less mechanical brake engagement.2. The Nussbaum-type function is adopted to handle the unknown control direction problem, which stems from an uncertain wind speed and consequent uncertain aerodynamic torque variation.Accordingly, the need for accurate wind speed measurement is avoided.3. The pitch actuator fault effects, including effectiveness loss, pitch angle bias, hydraulic leak, high air content in the oil, and pump wear, are compensated for automatically via an adaptive faulttolerant controller design.Also, the effect of blade aerodynamic characteristic changes, due to debris build-up and erosion, is considered and mitigated.The fault information, including fault type, size, and time, is estimated, which can be used for maintenance operations.4. Smooth pitch actuator saturation is designed to avoid the harsh and fast pitch actuator saturation phenomenon, which may increase the structural load on the wind turbine and result in performance degradation.Also, a neural network estimator is adaptively augmented in the proposed controller to obtain the uncertain aerodynamic torque.5.The control design is fulfilled in the backstepping framework, utilizing the virtual control concept.In this regard, the repeated differentiation of virtual control is required, which increases the complexity of the designed controller order.The DSC technique is used to eliminate this problem by introducing a first-order filter [25].
The rest of this paper is organized as follows.In Section 2, the wind turbine model is summarized.In Section 3, pitch actuator saturation and faults are introduced.In Section 4, the desired operational mode and objectives are introduced.Accordingly, the proposed controller is designed in Section 5, and the fault identification scheme is described in Section 6.The numerical evaluation of the proposed controller is addressed in Section 7 and the results are discussed.Finally, concluding remarks and open problems are given in Section 8.

Nominal Wind Turbine Model
The wind kinetic energy is captured by the blades and transferred into the rotor, rotating at ω .The effective wind speed, V , causes an induced aerodynamic torque, T , and thrust, F , modelled as [2]: where ρ is the air density and R is the blade length.Also, C and C are the torque and thrust coefficients.These coefficients are functions of the blade pitch angle, β, and tip speed ratio, λ, defined as λ = Rω /V [2].Considering the long elastic tower, T causes a fore-aft oscillation of the nacelle, i.e., a bending oscillation of the tower.This motion is modelled as [7]: where B , K , and M represent the damping ratio, and the elasticity coefficient of the tower and nacelle mass, respectively.Also, the nacelle displacement from its equilibrium position is represented by x .The effective wind speed at the rotor plane is then obtained as V = V − x , where V is the free wind speed, which itself is the wind speed before encountering the blades [9].The captured aerodynamic power by the wind turbine is then written as: where, C represents the power coefficient.Also, considering P = T ω , the relation between the power and torque coefficients is C = C λ.The empirical equation of C is stated as [26]: where 1/λ = 1/(λ + 0.08β) − 0.035/(β + 1), C = 0.5176, C = 116, C = 0.4, C = 5, C = 21, and C = 0.0068.Also, C can be approximated as [27]: where A = 0.006, A = 0.095, A = −4.15,A = 2.75, A = 0.001, A = 7.8, A = −0.00016,and A = −8.88.The drive train is used to increase the rotor speed, ω , and transfer the kinetic energy into the generator shaft, rotating at ω .The drive train is modelled as a two-mass system.The rotor and generator shafts have an inertia of J and J , respectively.The elastic gear meshing is considered, with the inclusion of the torsion stiffness, K , and the torsion damping, B .This elastic gear meshing leads to a torsional angle of twist, θ , defined as: where θ and θ are the rotation angle of the rotor and generator shafts, respectively.Also, the viscous friction at the bearings of the rotor and generator shafts are modelled with coefficients B and B , respectively.The drive train efficiency is η .So, the drive train is modelled as [15]: The generator speed, rotor speed, and their time derivative sensors are modelled as, ω , = ω + ν , ω , = ω + ν , and ω , = ω + ν , where ν , ν , and ν are noise contents [2].
The generator shaft kinetic energy is converted into electrical energy in the generator.Also, a converter is located between the generator and the electrical grid to adjust the generated power frequency [7].The generator internal electronic controller is much faster than the mechanical dynamic behavior of wind turbines.So, it is assumed that the generator torque, T , is adjusted at the generator reference torque fast enough to ignore the generator dynamic response.Also, the generated electrical power, P , is approximated as a static function given by [7]: where η is the generator efficiency.
The hydraulic pitch actuator rotates the blades to regulate β at the desired one, β , tuned by the pitch controller.The pitch actuator is modelled as [4]: where ω and ξ are the natural frequency and the damping ratio of the pitch actuator, respectively.The pitch actuator operational ranges are limited as β ≤ β ≤ β , β ≤ β ≤ β .Note that in this paper, (•) and (•) stand for the maximum and minimum allowable value of the variable (•), respectively.The pitch angle and its time derivative sensors are modelled as β = β + ν , β = β + ν and β = β + ν , where ν , ν , and ν are the noise contents [15].The numeric values of the wind turbine benchmark model parameters are given in Table 1 [7,15].

Pitch Actuator Saturation, Faults, and Blade Aerodynamic Characteristics Change
In reality, the achievable pitch angle range is bounded.Hence, the practical operational range of the pitch actuator is limited.So, the high wind speed variation and the consequent high pitch angle variation may lead to pitch actuator saturation, which consequently causes violation of the constrained power regulation.So, the pitch angle saturation phenomenon should be considered in the pitch angle controller design to avoid any abrupt and long-lasting saturation, and smoothly pass any possible saturation period of the pitch actuator.For the given wind turbine model, the pitch actuator saturation, H(β ), can be considered as: where β = 30° and β = −2°, as in Table 1.H(β ) is illustrated in Figure 1.So, this actuation saturation function is integrated into the pitch actuator mechanism (9) as: Obviously, this saturation function of the pitch angle is non-smooth with sharp saturation behavior, which may cause pitch actuator failure [28].So, it is desirable to approximate this saturation behavior by a smooth function and to pass from each saturation period fluently.In this regard, the following smooth saturation function is proposed as: where ϱ = 2 , ϱ = 30 , Ρ = exp(ϵ + ηβ ) , ϵ = 0.5 ln ϱ/ϱ , and η is a positive constant to be selected.S(β ) is always in −ϱ, ϱ for all β ∈ ℝ. S(β ) is illustrated in Figure 1, for η = 0.1.Then, H(β ) can be expressed as: where D(β ) is the difference between S(β ) and H(β ).The bounded property of the function, S(β ), and saturation function, H(β ), yields the function, D(β ), to be bounded as, where D is a positive and unknown constant.For ease of pitch actuator controller design, the mean value theorem is employed on function S(β ) to get: where β = mβ + (1 − m)β and m ∈ (0, 1).By choosing β = 0 and using the fact that S(0) = 0, Equation ( 14) becomes: where S = 2η(ϱ + ϱ)/(Ρ + Ρ ) .S ∈ (0.2, 1.65) for β ∈ (−ϱ, ϱ ) .Then, S is a positive variable.Now, considering Equations ( 11), (13), and (15), the pitch actuator dynamic behavior with the smooth saturation function can be written as: Wind turbine operation in harsh offshore sites may lead to pitch actuator dynamic change, bias, and effectiveness loss.The dynamic change is because of the pressure drop due to hydraulic oil leakage, high air content in the oil, and pump wear.These dynamic change cases cause a slow pitch actuator response [7].Consequently, power regulation in full load operation is not satisfactorily achieved.The dynamic change is considered as the change of the natural frequency and damping ratio in the pitch actuator in Equation (16).The characteristics of these changes are summarized in Table 2 [4,7], where N, HL, PW, and HAC represent normal, hydraulic leaks, pump wear, and high air content situations, respectively.Also, ω , and ξ are the natural frequency and damping ratio, respectively, in the situation X.Also, α and α are fault indicators.The dynamic change case effects are illustrated in Figure 2, where the initial pitch angle is set to 5 ° and β = 0 °.It is obvious that the response for all dynamic change cases is slower than the normal one.The dynamic change is modelled as an uncertainty, which should be attenuated by the pitch angle controller.The dynamic changes are modelled as a convex function of the normal natural frequency and normal damping ratio [18].So, the pitch actuator (Equation ( 16)) is rewritten, including the dynamic change effect, i.e., added as an uncertainty in the model, as: where The pitch actuator output can be corrupted by an unanticipated fault, modelled as an additive bias and/or effectiveness loss.These faults deviate the pitch angle from the desired one [16].These faults are modelled as: where Φ(t) represents the unknown pitch actuator bias [29].Also, ρ(t) is the unknown effectiveness of the actuator, which is 0 < ρ(t) ≤ 1 , where ρ(t) = 1 indicates full effectiveness of the pitch actuator and ρ(t) = 0 is total actuator loss [13,29].β (t) is the reference pitch angle, which is designed by the pitch controller.It is obvious that in the case of full pitch actuator effectiveness and no pitch bias, β (t) = β (t).The pitch actuator dynamic behavior, Equation ( 17), associated with pitch actuator bias and effectiveness loss, can be rewritten as: The wind turbine operation in the presence of rain, snow, and dirt leads to blade erosion or debris build-up, which, in turn, leads to a blade aerodynamic efficiency reduction.Consequently, the captured aerodynamic power is decreased.On the other hand, the power regulation is not satisfactorily achieved due to the changed blade aerodynamic profile.The debris build-up effect is modelled as an aerodynamic change, ∆T ,∆ , due to a change in the power coefficient as C = C + ΔC .So, it is very important to foresee this potential change in the controller design.It is worth noting that debris build-up is challenging to detect, as it is hard to identify if the reason for the reduced power is the blade's debris or simply that the wind speed has decreased.As debris build-up occurs slowly, this change is mostly assumed to lie within the annual maintenance/inspection, in which the blades are cleaned/replaced.So, this paper aims to design a pitch controller that is insensitive to debris build-up that guarantees graceful degradation up to the next planned maintenance of the blades.

Desired Operational Mode and Control Objectives
As mentioned earlier in Section 1, in the full load region, it is desirable to keep P at the nominal value, P , , to avoid overspeeding and consequent brake engagement.To this aim, by taking Equation (8) into account, (i) T is to be kept at the nominal value, T , , and (ii) ω is to be kept at the nominal value, ω , , to have nominal power generation as, P = η T ω = η T , ω , = P , [7].Considering the fast response of the electrical generator, the objective (i) is achieved via tuning the generator reference torque at T , , which leads T to be set to T , quickly.The objective (ii) is achieved by adjusting β of the pitch actuator.This leads to tuning the aerodynamic torque, and consequently, the rotor speed and the generator speed [4].The wind speed is a highly stochastic variable.So, accurate nominal power generation is very challenging, and in the case of improper controller design, it may lead to overspeeding and braking.So, it is very beneficial to guarantee that the generated power and speed do not violate the given constraint, within which the mechanical brake is not engaged.It should be noted that for the power control purpose, the generator torque controller is not active.So, the faults in the generator are not considered in this paper and it is assumed that generator faults have already been accommodated for using the generator controller [4].
The controller is designed to adjust the reference pitch angle, β , and keep the rotor speed as close as possible to the nominal one, i.e., ω , , never violating the given constraint, in the presence of wind speed variation, disturbance, and pitch actuator faults and saturation.The primary objective of this paper consists of satisfying the above-mentioned requirements.
By considering structurally safe operation of the wind turbine, it is desirable to keep the drive train torsion angle variation, θ , as small as possible, which in turn leads to the drive train stress reduction.Accordingly, θ = 0 leads to N ω = ω [17].So, it is beneficial to keep the ratio between the rotor and generator speeds at the drive train ratio [14].As the generator speed has to follow the signal, ω , , then the rotor speed is kept at ω , = ω , /N [6].Moreover, the condition, θ = 0, with zero initial drive train torsion angle, leads to θ = 0, i.e., the reduced drive train stress trajectory [6].Accordingly, considering Equation ( 7), the desirable operational mode of the wind turbine with reduced drive train stress is given by [6,14]: where, 20) can be rewritten as, where c = a + a b , c = a a + a b , and c = a a , c = a b .Considering Equations ( 1) and ( 21), it is obvious that the rotor speed is controlled by regulating the pitch angle and the consequent aerodynamic torque.In this paper, it is assumed that at any operational point of the wind turbine, T is not a singular function.Also, for any pair of (V , ω ), there is a given pitch angle, i.e., β * , leading to the nominal power generation [9].So, in the presence of wind speed variations, β * will be set to the value that satisfies the control objective.For the considered benchmark model, the β * diagram is illustrated in Figure 3 [7].Note that as the wind speed is considered an uncertain disturbance, then β * is an unknown variable.It is obvious that the rotor dynamic behavior (Equation ( 21)) is a non-affine function of the pitch angle [9].Linearization is one obvious solution.However, it would lead to high inaccuracy.So, the paper solved this problem by using the mean value theorem in this paper.As stated earlier, T is not a singular function for any triple pair (V , ω , β) in the operational range of the wind turbine.So, according to the mean value theorem, for any given pair of (V , ω ), there exists Ξ ∈ (0, 1), such that [9]: where Using T in Equation ( 1) and C in Equation ( 4), for the considered wind turbine model, whose parameter values are given in Table 1, the diagram of ∂T / ∂β in the full load region is shown in Figure 4. Remark 1.In Figure 4, it is evident that − ≤  / ≤ − < 0, with 0 <  < .This means that as the wind speed increases, by increasing the pitch angle, the aerodynamic torque decreases.
Taking the time derivative of Equation ( 22) yields: where and ∂T / ∂β = T , .∆T ,∆ is due to ΔC , which is the result of the changes in the blade aerodynamic characteristics.The debris build-up and erosion occur slower than the mean time to the maintenance of the blades.So, all terms that are contributing to ∆T ,∆ are assumed to be bounded, then ∆T ,∆ is bounded as ∆T ,∆ ≤ k , where k is an unknown positive constant.Also, it should be noted that as the wind speed is not accurately measurable, T , in Equation ( 23) is an unknown variable.Substituting Equation (23) into Equation ( 21) leads to: Now, by considering Equation ( 19), the rotor dynamic behavior (Equation ( 24)) can be rewritten as: This expression describes the wind turbine rotor dynamic behavior in the desired operational mode, which takes into account possible pitch actuator dynamic changes.Also, smooth pitch angle saturation is included.

Constrained Fault-Tolerant Controller Design and Stability Analysis
In this section, the constrained fault-tolerant pitch controller is designed to guarantee that the generated power is kept within given constraints, in the presence of wind speed variation, disturbance, faults, and saturation.The robust stability of the wind turbine closed-loop system with the proposed controller is proved.First, some technical preliminaries are provided, which will be used for the controller design.

Technical Preliminaries
The wind speed is uncertain as the wind speed is measured with an anemometer, usually placed at the back of the nacelle.Therefore, its measurement is affected by the turbulence generated by the rotor.So, the wind speed is considered as an uncertain disturbance.Accordingly, the aerodynamic torque, T , is not accurately available.On the other hand, T is contributing to the rotor dynamic response (Equation ( 25)).So, T should be estimated to be used in the proposed controller structure.In this paper, an RBF neural network is designed to estimate the aerodynamic torque [30].To this end, T is approximated as [9,31]: where  * is the optimal weight vector, () = [h (), h (), … , h ()] ∈ ℝ is the known basis function vector, s > 1 is the number of neural network nodes,  = [T , ω , β] ∈ Ω  , ε ∈ ℝ is the approximation error, and h () is selected as a Gaussian function given by [9,21]: where  = [ϑ , , ϑ , , ϑ , ] is the i centre vector of the inputs, is the width vector of the Gaussian functions [32].T , is the approximation of T provided by the RBF, described as: The optimal weight vector, i.e.,  * , is defined as It should be noted that ε is bounded as |ε| ≤ ε , with an unknown bound, ε > 0. Now, the following definitions and lemmas are given, which will be used in the proposed controller design.
As T , is an unknown variable, contributing to the gain of β in Equation (25), this leads to the unknown control direction problem.To tackle this issue in the controller design, the Nussbaumtype function is utilized, which is defined as follows.
The BLF function is defined as follows, which is used in the constrained control construction.
Definition 2 [32].If the scalar function, (), is positive definite continuous with respect to the solution of the system,  = () ,on an open region,  , then () is a BLF with continuous first-order partial derivatives within all .As  approaches the boundary of the region, , then () approaches infinity.Finally, () satisfies () ≤ , ∀ ≥ 0 along the solution of  = () for (0) ∈ , and some positive constant, .
The following definition is given for the boundedness of the closed-loop system.

Controller Design Procedure
The main objective of the designed controller is to keep the rotor speed and acceleration within constraints, which in turn leads to bounded power generation around the nominal one.The proposed controller design requires the definition of the rotor speed tracking error and its time derivative as follows: where ω , is the desired rotor speed.As stated earlier, ω , in the full load region is ω , .z is a virtual control.Here, to avoid repetitive differentiation of z , which increases the implementation complexity, the DSC technique is utilized, which requires the filtering of z by means of a stabilizing function, α , to be determined.Let α pass through a first-order filter with a time constant, τ , defined as: The output error of the first-order filter is χ = z − α with its first-time derivative, z = −χ /τ .A Lyapunov function is chosen as: where Λ = πξ /2, ξ = e /k , whilst k represents a constraint on e .It should be noted that V is continuous in the set Ω = e : −k < e < k .V is positive definite and its first term captures the BLF characteristics of the modified tracking error, ξ , according to Definition 2. The first-time derivative of ξ is obtained as: On the other hand, the first-time derivative of Equation ( 31) can be obtained as: The virtual control, α , is designed as: where, γ is a positive design parameter.The substitution of Equation ( 34 Since α is a function of ω , ω , , and ω , , it can be shown that: Considering Equation (37), α is a continuous function.Then, given δ , and any positive number, δ , the set, Ω , : = {ω , ∈ ℝ: ω , + ω , + ω , < δ , }, for all initial conditions satisfying Ω : = {[e , χ ] : (k tanΛ )/π + χ /2 < δ } is compact [25].Thus, α has a maximum constant value, M , in the compact set, Ω , × Ω , for given initial conditions [28].So, based on Young's inequality: According to Equations (36) and (38), Equation (35) can be rewritten as: The parameter, τ , is selected as τ < 1 to satisfy (2/τ − 2) > 0. Also, considering Lemma 3, −γ e sec Λ < −(γ k tanΛ )/π holds true.Consequently, Equation (39) is rewritten as: where σ , = min {γ , (2/τ − 2)} and σ , = M /2.
Considering the measurement noise and the RBF neural network estimation of aerodynamic torque (Equation ( 26)), the rotor dynamic response (Equation ( 25)) can be rewritten as: where , and f = Δf /ω , + S ω , Φ. Considering the bounded achievable, β, β , and β , the boundedness of Δf , S , and Φ are concluded.This, in turn, leads to the bondedness of the fault signal, f, in Equation (41), i.e., |f| ≤ f ̅ , where f ̅ is an unknown positive constant.Also, the first time derivative of f is assumed to be bounded, i.e., f ≤ ρ , where ρ is an unknown positive constant [4].Indeed, it is assumed that the applied fault is a slowly varying function of time.Besides, the sensor noise contents are bounded, which is a reasonable assumption [2,17].Accordingly, by considering Figure 4 and |ε| ≤ ε , the disturbance, d, is bounded as |d| ≤ D , where D is a positive unknown constant [5].Finally, it is easily seen that G is unknown but bounded, due to the presence of T , .Now, the proposed pitch angle controller is designed as: ζ = e sec Λ υ , associated with the adaptive laws: to estimate the fault, the RBF neural network weights and disturbance, respectively, where, η , η , η , η , γ , σ , σ , and σ are positive design parameters.Also, Γ ∈ ℝ is a design matrix such that Γ = Γ > 0. Accordingly, the estimation errors are defined as: Considering the bounded applicable fault, it is practically reasonable to assume the estimation error is bounded as f ≤ ρ , where ρ is an unknown positive constant [4].It should be noted that this bound is only used to analyze the stability of the closed-loop system.Nevertheless, this will not be used in the designed control structure, as it is assumed to be unknown.So, the actual estimation will not be required in setting up and implementing the control scheme.Now, a Lyapunov function is selected as: where Λ = πξ /2, ξ = e /k and k is a considered constraint on e .It is worth noting that V is continuous in Ω = e : −k < e < k .V is positive definite and its first term captures the BLF characteristics of the modified tracking error, ξ , according to Definition 2. The time derivative of ξ is obtained as: The first-time derivative of Equation ( 46) can be obtained as: Substituting Equations ( 42)-(44) into Equation (48) leads to: where: Also, with the aid of Lemma 4, the following inequality is obtained: where λ (Γ ) is the maximum eigenvalue of Γ .Finally, considering Lemma 3, it leads to: Using Equations ( 50)-(53) in Equation (49), the following inequality is obtained: where σ , = min {γ , σ , σ , σ /λ (Γ )} and σ , = a ω , LKη /2ξ + a LKη /2ω , ξ + D Kη + σ D /2 + a Lf ̅ Kη /2ξ + σ f ̅ /2 + ρ ρ + σ ‖ * ‖ /2.Now, the main property of the designed pitch controller is proven by Theorem 1.
From the above-mentioned analysis, the objectives (i), (ii), and (iii) are achieved as follows: (i).
Consider the boundedness of V , V , e , and e .Therefore ω and ω are bounded.Now, from the boundedness of f , d , and  , the boundedness of α , ζ , υ , f , d ,  , and consequently β is proven.(ii).
As part of the closed-loop system analysis, it is shown that the tracking errors, e and e , always stay in the sets, Ω = e : |e | < k and Ω = e : |e | < k , respectively.(iii).
Multiplying both sides of Equation ( 59) by exp(σ , t) yields: Thus, the integration of Equation ( 60) over [0, t], becomes: where  = V (0) − c , /σ , exp(−σ , t) + c , /σ , .From the definition of V , it can be shown that: where  = k 2tan π/k /π .If V (0) = c , /σ , , then it holds that |e | < k 2tan πc , /k σ , /π .If V (0) ≠ c , /σ , , it can be concluded that given any  > k 2tan πc , /k σ , /π, there exists T, such that for any t > T, it has |e | < .As t → ∞, |e | < k 2tan πc , /k σ , /π, which implies that e can be made arbitrarily small by selecting the design parameters appropriately.Considering Definition 3, and the objectives (i), (ii), and (iii), it is guaranteed that the closed-loop system is UUB.This completes the proof.□ Remark 2. In Theorem 1, the objective (i) implies that the wind turbine equipped with the proposed pitch angle controller is stable.The objective (ii) states the constrained rotor speed and acceleration are guaranteed.Accordingly, the generator speed and generated power are retained in the given bounds.Considering Section 3, then the efficient power regulation requirements are met, satisfying the power grid demand.In this manner, both rotor overspeeding and mechanical brake engagement are avoided.The objective (iii) represents the expert's knowledge in the implementation stage of the proposed controller to satisfactorily make the tracking error small.These objectives are satisfied in the presence of pitch actuator faults, dynamic change, saturation, and blade aerodynamic characteristic change.

Fault Identification Scheme
In this section, a scheme is given to identify the pitch actuator fault, including the pitch bias, Φ, effectiveness loss, ρ, dynamic change, Δf , and aerodynamic characteristic change, ∆T ,∆ .To this aim, the auxiliary signals are calculated and compared to the estimated fault, f , to identify the case of the dynamic change.Considering Table 2 and the definition of Δf in Equation ( 17), the auxiliary signals are computed as follows: where in f , , the auxiliary signal is calculated for the fault case, X, which includes fault free, pump wear, hydraulic leakage, or high air.In order to finalize the pitch actuator fault identification scheme, first, the dynamic change case is considered, assuming no pitch actuator bias.Considering f = Δf /ω , + S ω , Φ, in the absence of Φ, it can be obtained that f = Δf /ω , .So, using the auxiliary signals, f , , reported in Equation ( 63), the most similar one to f is identified as the dynamic change case.Therefore, similarity indices are needed, which enhance the fault identification task.The indices adopted in this paper are the root mean squared error (RMSE) and variance accounted for (VAF), defined as follows: where T is the given operation period of the wind turbine.In the ideal fault identification case, the RMSE and VAF indices are zero and 100%, respectively.Accordingly, the dynamic change of X with the RMSE and VAF indices close to zero and 100%, respectively, is selected as the corresponding dynamic change case, which is indicated as X .The pitch bias, Φ, is considered as an added constant on f.So, having pitch bias occur with the dynamic change, the RMSE index only deviates significantly from zero.However, the VAF index still indicates the correct dynamic change properly.Now, after the identification of X , the pitch actuator bias is estimated.Considering f = f , /ω , + S ω , Φ, where f , is the calculated auxiliary signal using Equation (63) for the identified dynamic change, X , the estimation of pitch actuator bias, Φ , is computed as follows: Finally, if neither the dynamic change case nor pitch actuator bias are identified, and meanwhile, the fault-free case is not identified, then it can be concluded that the faulty case is either pitch actuator effectiveness loss or aerodynamic characteristic change.Considering Equation (41), it is clear that the aerodynamic characteristic change is considered as an additive disturbance, and attenuated by the proposed controller.On the other hand, the effectiveness loss is contributing in the control gain, i.e., G. So, the estimated fault, f , is affected by the effectiveness loss and is insensitive to the aerodynamic characteristic change.Accordingly, the given period is considered, in which no pitch actuator dynamic change, bias, or fault-free cases are identified.Then, if f considerably deviates from zero, this leads to the identification of the effectiveness loss.Otherwise, the aerodynamic characteristic change is identified.Therefore, the fault isolation task is accomplished.It should be mentioned that this fault identification scheme is robust against the disturbance, d, in Equation ( 41), as its effect is guaranteed to be attenuated using the proposed controller.

Numerical Evaluation and Comparison
In this section, numerical simulations are conducted to evaluate the features of the controller (Equation (42)).Moreover, the available industrial PID controller is briefly introduced, for comparison with the proposed controller performance.Suitable numerical metrics are introduced to quantify and compare the performance of the proposed and PID controllers.

Industrial Baseline PID Controller
The most commonly adopted industrial controller for power regulation of wind turbines in full load operation is the PID controller, due to its simplified implementation and effectiveness [4].The PID controller is thus used to regulate the pitch angle based on the generator speed tracking error, e , defined as: Accordingly, the PID controller, used for tracking the blade pitch angle given the reference pitch angle, has the form [4,7]: where K , K , and K are the proportional, integral, and derivative gains of the controller, respectively, to be set via traditional methods, in order to guarantee system stability as well as satisfying performance.K , K , and K are mostly chosen as constant gains for the whole operational region, although some works proposed the use of different gains for each operating condition of the plant [7].The values of the PID gains here settled as K = 1, K = 4, and K = 0 [4,26].
In the structure of the industrial controller (Equation ( 67)), the sensor noise, ω , , is not necessarily attenuated and may be amplified, even if a filter is used to remove noise content [7].Also, any possible loss of effectiveness, ρ(t), and blade aerodynamic characteristic change, ∆T ,∆ , are not analytically attenuated with this solution.Moreover, this controller does not guarantee that any pitch actuator bias, Φ(t), and dynamic change, Δf , is correctly managed.These remarks will help to highlight the advantages of the proposed controller compared to the PID controller, which are analyzed by means of a simulated example in the next section.

Performance Metrics
The comparative numerical performance metrics are defined in this section.The difference between the generator speed and the nominal one is considered as the first metric, defined as: where T is the given operation period of the wind turbine.Similarly, the difference between the generated power and the nominal one is considered as the second metric, defined as: Obviously, it is desirable to keep C1 and C2 as close to zero as possible.The maximum power deviation from nominal is calculated as: C3 indicates the instantaneous power deviation from the nominal, which may cause a sudden break down.In contrast, C2 accumulates all power deviation, which may lead to gradual failure.So, it is expected that C3 never violates the settled constraint.Also, the drive train torsion angle is calculated as: Which represents the applied drive train torsion stress due to variation in aerodynamic torque, as a result of pitch angle variation.It is desirable that the proposed controller maintains the value of C4 near the one provided by the PID regulator, which is accepted in industrial practice.Finally, in order to evaluate limited variations of β, C5 and C6 are defined as: (72)

Simulation Results
In this section, the numerical simulations are reported to evaluate the performance of the proposed controller (Equation ( 42)), both in fault-free and faulty situations.Also, a comparison is made to the industrial PID controller (Equation ( 67)) to illustrate the effectiveness of the proposed controller, considering suitable numerical metrics.It should be noted that different fault scenarios, including single and simultaneous occurrences, and wind speed variations are introduced to investigate the robustness of the proposed controller.
The parameter values of the proposed controller are summarized here.The constraints on the rotor speed and its time derivative are selected as: k = 0.02 rad/s, k = 0.04 rad/s .With these values, the inequalities, ω − ω , ≤ 0.02 rad/s and |ω | ≤ 0.04 rad/s , are satisfied.Consequently, considering the operational mode, the constraints on the generator shaft speed and generated power are ω − ω , ≤ 1.9 rad/s and P − P , ≤ 0.056 MW.As it is clearly highlighted in Equation ( 42 Also, the centres and width of the RBF neural network are selected as: and φ = 10ones(10,1), respectively.The other control parameters values are selected as: , L = 700000 , and Γ = I × .

Fault-Free Situation
Firstly, the performance of the proposed controller is analyzed for a simulation time of 1500 (s).Fault-free conditions are also considered, with wind speed with mean of 19.84 (m/s) and standard deviation of 1.94 (m/s), as shown in Figure 5.The design should lead to the following properties: (i) The considered constraints are not violated, (ii) the pitch angle saturation is smoothly avoided, and iii) the performance is improved compared to the conventional PID controller.The corresponding rotor speed, rotor acceleration, generator speed, and generated power, using the proposed controller, are shown in Figures 6-9, respectively.It can be verified that the considered constraints are not violated.On the other hand, with the same wind sequence, the PID controller results are given in Figures 6-9.The PID controller is not able to keep the corresponding outputs within the considered constraints, in the presence of the wind speed variation.Moreover, the obtained reference pitch angle using both controllers are compared in Figure 10, in which it is shown that the PID controller leads to pitch actuator saturation.In contrast, the proposed controller has smoothly avoided saturation.It should be noted that as the proposed controller maintains the rotor speed within the constraints, despite the high wind speed variation, faster pitch angle variations are generated.As mentioned in Section 3, the proposed controller is designed on the desired trajectory of the wind turbine, for which the drive train torsion angle is reduced.In order to analyze this issue, the induced drive train torsion angle using the proposed and PID controllers is depicted in Figure 11.It is shown that the induced drive train torsion angle using the proposed controller has values close to the ones achieved via the PID controller.This implies that the proposed controller has not considerably increased drive train torsion and, consequently, stress, despite the wind speed variation and more accurate nominal power tracking.The uncertain aerodynamic torque estimation is shown in Figure 12, in which the actual aerodynamic torque is reported to evaluate the estimation efficiency.It is highlighted that the aerodynamic torque is estimated quite accurately, and has been kept around the nominal one, the same as the actual aerodynamic torque.Finally, to accurately compare the results, the performance metrics using both controllers are summarized in Table 3.It can be noted that the performance metrics, C1 , C2 , and C3 , have been considerably reduced by using the proposed controller.These results correspond to Figures 6,8, and 9. Also, the metric C4 shows the same induced drive train torsion angle rate, as illustrated in Figure 11.The metric C5 shows the advantage of using the smooth pitch angle saturation, as depicted in Figure 10.As remarked above, the accurate nominal power tracking needs higher pitch angle change, in the presence of high wind speed variation.This aspect can be verified considering the index C6.So, it can be concluded that the proposed controller improves the wind turbine performance in the fault-free case compared to the industrial PID controller.The section evaluates the fault tolerance capabilities of the proposed controller in the presence of faults.It is expected that the mentioned constraints are not violated, whilst the fault effects are attenuated.Also, the estimated faults are analyzed with respect to the indices described in Equation (64).The results using the PID controller are also reported to study the effect of each fault as well as to highlight the benefit of the proposed controller.The considered fault scenario is defined in Table 4.The occurrence of single faults is considered, in order to accurately study their individual effects, as highlighted in Table 4.The same wind speed sequence shown in Figure 5 is considered here.Figures 13-16 illustrate the rotor speed, rotor acceleration, generator speed, and generated power with respect to the corresponding constraints, respectively, using both the proposed and PID controllers.Compared to Figures 6-9, it is obvious that the PID controller performances are degraded, while the proposed controller is able the attenuate the fault effects and maintain the considered outputs within the corresponding constraints.The designed reference pitch angle values using both controllers are compared in Figure 17.It highlights that the simple PID controller has led to pitch actuator saturation while saturation is smoothly avoided using the proposed controller.Moreover, it is worth noting that as the proposed controller is trying to counteract the pitch actuator dynamic change, the dynamic change has reduced the speed of the pitch actuator.Accordingly, it has led to slightly higher pitch angle variations compared to the fault-free case, during the dynamic change periods.Consequently, the drive train torsion angle induced by the proposed controller was increased, as illustrated in Figure 18.In order to analyze the changes of all variables, X, due to the fault effect with respect to the fault-free case, the following relation is defined: where δX is the change in the considered variable X, X is its fault-free value, and X represents the corresponding value in the faulty situation.In Figures 19-21, δP , δω , and δω , are illustrated, respectively, for the proposed and PID controllers.It is shown that the pitch actuator bias and the effectiveness loss have led to considerable limitations of the achievable PID controller performance.However, these performance degradations are significantly attenuated using the proposed controller.Also, the effect of the considered blade aerodynamic change, i.e., ∆T ,∆ = 5%, has led to smaller variations.It should be noted that the hydraulic leak has notable effects in the performance degradation.To accurately study the effect of faults on the pitch actuator response, δβ is depicted in Figure 22.Evidently, changes in the pitch actuator dynamic have been considerably attenuated using the proposed controller while the slower pitch actuator leads to worse performance with the PID controller.Also, the effect of pitch actuator bias is completely removed, since its effect is compensated by using the estimation of the bias itself.The situation is even worse when the PID controller is exploited.The same result is obtained considering the effectiveness loss fault.In the case of blade aerodynamic change, both controllers have led to the same trend in δβ while the variation using the PID controller is significantly higher.However, the fault is removed after 1300 (s), whilst its effect reduces the PID controller performance.Finally, to quantitatively compare the performance of the controllers, the values of the performance indices are summarized in Table 5.It is worth noting that the achieved performance of the proposed controller in the presence of faults are similar to the faultfree conditions.This represents the main point of the fault-tolerant control design.In this section, with the aid of the estimated fault, f , and the calculated auxiliary signals of Equation (63), the fault identification performance is analyzed using the indices in Equation (64).The identification task includes fault detection, isolation, and its reconstruction (i.e., its shape).Consequently, the estimated pitch actuator bias, Φ , is obtained using Equation (65), in which S is calculated using S = 2η(ϱ + ϱ)/(Ρ + Ρ ) . It should be noted that the indices in Equation ( 64) and the estimated fault, f , should be computed and compared in every time step of simulation in order to accurately identify the fault.However, since the overall performance of the proposed fault identification technique is analyzed, the comparisons are performed in each fault period of Table 4. Also, as the calculated auxiliary signal for the fault-free case, i.e., f , , is always zero, then the VAF in Equation ( 64) is calculated as var f × 100% .Now, to investigate the fault identification capability of the proposed controller, the estimated f signal is shown in Figure 23.Also, the auxiliary signals computed from Equation (63) are depicted in Figure 24.Moreover, by means of the fault identification indices in Equation (64), summarized in Table 6, the estimated pitch actuator bias is shown in Figure 25.As remarked above, to fulfil the fault identification task, the values of the RMSE and the VAF indices closer to zero and 100%, respectively, determine the fault case.Considering Table 6, it is clear that in each fault-free period, the indices are indicating the fault-free case.Also, the estimated pitch actuator bias is zero.It is worth noting that in the fault-free periods, the VAF indices calculated for the dynamic change cases, are negative.Indeed, this result is justified considering the form of Equation (64).Comparing Figures 23 and 24, it can be pointed out that in the fault-free periods, var f , − f is greater than var f , , which leads to negative VAF. Obviously, this is not the case considering the calculated VAF, using the fault-free auxiliary signal.For the case of the pitch actuator dynamic changes, it can be verified that the selected indices lead to accurate identification of the corresponding actual dynamic change case.In all dynamic change cases, the estimated pitch actuator bias is zero.So, the dynamic change cases are clearly distinguishable from the pitch actuator bias.Evidently, during the pitch actuator bias period, none of the indices satisfying the considered conditions can be selected.On the other hand, the pitch actuator bias is precisely estimated.The situation is different for the cases of the pitch actuator effectiveness loss and the aerodynamic characteristic change.Obviously, neither the dynamic change case nor the fault-free case is selected, as the corresponding indices do not satisfy the given conditions.On the other hand, the pitch actuator bias is estimated as zero.Therefore, it can be concluded that these two periods correspond to the pitch actuator effectiveness loss and/or aerodynamic characteristic change.However, in order to correctly identify these two cases, the estimated fault is first considered, i.e., Figure 23.From Equation (41), it can be noted that the aerodynamic characteristic change is described as an additive disturbance while the effectiveness loss affects the control gain, i.e., G in Equation (41).The proposed controller has shown robustness features with respect to the considered disturbances.Accordingly, the estimated fault is affected by the effectiveness loss and it is insensitive to the aerodynamic characteristic change.Therefore, between 1000 and 1100 (s), as the estimated fault is different from zero, the effectiveness loss case is identified.Also, between 1200 and 1300 (s), as the estimated fault is zero, the aerodynamic characteristic change case is identified.

Robustness Evaluation
In this section, the proposed controller is further evaluated in terms of robustness to different wind speed sequences and fault scenarios.The wind speed is shown in Figure 26, with a mean of 20.41 (m/s) and standard deviation of 3.01 (m/s) for 1100 (s).Compared to the former wind speed reported in Figure 5, the signal considered in Figure 26 presents more fluctuations.So, it is more challenging for the controller to satisfy the objectives.Also, for the different fault scenario reported in Table 7, the faults occur simultaneously for a longer period.Also, the pitch actuator bias and the aerodynamic characteristic change values are increased.On the other hand, the pitch actuator effectiveness is decreased.For the sake of brevity, only the generated power in the fault-free and fault cases are considered for both the controllers.Considering these conditions, the performance metrics are compared.The generated power is illustrated in Figures 27 and 28, in both the fault-free and faulty situations, respectively.It is obvious that in both situations, the generated power is maintained within the prescribed constraints.The performance metrics are summarized in Tables 8 and 9, which further confirm that the proposed controller can successfully maintain reliable performance under faulty conditions.Its performances are clearly better than the ones achievable with the PID controller.Considering the simulation results in the fault-free and first fault scenario of Table 4 using the wind profile of Figure 5, it can be pointed out that the proposed controller is able to reduce the degradation of the wind turbine performances.Firstly, it is shown that the generated power never violates the given bound.It should be noted that this bound is appropriately selected by the wind turbine manufacturer.Indeed, this bound represents the safe-to-operate bound in full load operation without mechanical brake engagement.On the other hand, mechanical braking increases the induced fatigue load on the drive train.Also, the rotor acceleration is constrained, thus leading to reduce the torque stress on the rotor shaft.So, it can be concluded that the proposed controller leads to induce lower fatigue load and stress to the drive train.One obvious benefit of the proposed controller is that harsh long-lasting pitch actuator saturation is avoided.In fact, using the smooth pitch actuator saturation function, the speed of the pitch actuator response increases before the saturation, as highlighted in Figure 1.This characteristic leads to improved power regulation.It illustrates that the induced drive train torsion angle rate was kept at the same level as the baseline controller one.So, if the behavior of the baseline PID regulator is universally accepted in industrial control, the proposed controller can be used by industries to reduce the applied torsional stress.It is evident that the estimated aerodynamic torque fluctuates around the nominal one due to the inherent features of the RBF neural network and the Gaussian basis functions.However, as mentioned in the controller design procedure, the estimation error is bounded.This is obvious when comparing the actual aerodynamic torque and the estimated one.As remarked earlier, the pitch actuator dynamic change leads to a slower response of the pitch actuator and consequently poor power regulation.This phenomenon is highlighted in Figures 19 and 22.Nevertheless, the proposed controller was able to attenuate this effect, which is the same as the pitch actuator bias, using the properly designed fault estimator.On the other hand, the effects of the pitch actuator effectiveness loss and debris build-up are mitigated appropriately, satisfying the performance objectives.Also, using the fault estimator information alongside the proposed fault identification scheme, different faults are identified.Similar results can be obtained for even more severe instantaneous faults and higher wind speed variation, as illustrated in Figures 27 and 28.

Conclusions and Open Problems
This paper proposed a novel adaptive constrained control methodology for wind turbine power regulation subject to actuation failures as well as unknown system dynamics.In contrast to previous works where an unknown wind speed observer/estimator was needed, using the Nussbaum-type function, the proposed method was able to handle the unpredictable wind speed variation effects in the control design without requiring accurate wind speed measurement.The constrained rotor speed and generated power were guaranteed while the pitch actuator remained within the desired bounds.Using the barrier Lyapunov function in conjunction with the concept of dynamic surface control, a constrained stable control structure with cheap computational cost was developed.In addition, utilizing a radial basis functions neural network together with a proper fault-tolerant scheme, a robust and adaptive scheme was developed without the need for precise information about either the wind turbine model or the pitch actuator faults.Numerical simulations were finally performed to validate the effectiveness of the reported theoretical developments, and comparisons with the available industrial controller performance were shown.
Finally, by considering the proposed controller in this paper and the investigated results, the future research direction of this paper can be outlined as the validation of the proposed controller by means of data acquired from real or experimental-scale wind turbines for the whole operational region, achieving the Industry 4.0 requirements.

Figure 2 .
Figure 2. Pitch actuator response in situations: normal (black line), hydraulic leak (red line), pump wear (blue line), and high air content (green line).

Figure 3 .
Figure 3. Diagram of β * (blue line), and lower and upper pitch angle bounds (red dashed lines).

Figure 4 .
Figure 4. ∂T / ∂β diagram in full load operation (blue line), and upper and lower bounds (red dashed lines).
), a Nussbaum-type function is needed.In this paper, the Nussbaum-type function, N(ζ ) = ζ cos(ζ ), is used, which fulfils Definition 1.The RBF neural network structure has s = 10.

Figure 6 .
Figure 6.Rotor speed using the proposed controller (dark blue line), PID controller (red line), nominal rotor speed (light blue line), and constraints (green line), in the fault-free situation.

Figure 7 .
Figure 7. Rotor acceleration using the proposed controller (dark blue line), PID controller (red line), and constraints (green line), in the fault-free situation.

Figure 8 .
Figure 8. Generator speed using the proposed controller (dark blue line), PID controller (red line), nominal generator speed (light blue line), and constraints (green line), in the fault-free situation.

Figure 9 .
Figure 9. Generated power using the proposed controller (dark blue line), PID controller (red line), nominal power (light blue line), and constraints (green line), in the fault-free situation.

Figure 10 .
Figure 10.Reference pitch angle using the proposed controller (dark blue line) and PID controller (red line), in the fault-free situation.

Figure 11 .
Figure 11.Induced drive train torsion angle rate using the proposed controller (dark blue line) and PID controller (red line), in the fault-free situation.

Figure 12 .
Figure 12.Actual aerodynamic torque (red line), estimated one (dark blue line), and nominal one (light blue line).

Figure 13 .
Figure 13.Rotor speed using the proposed controller (dark blue line), PID controller (red line), nominal rotor speed (light blue line), and constraints (green line), with the first fault scenario.

Figure 14 .
Figure 14.Rotor acceleration using the proposed controller (dark blue line), PID controller (red line), and constraints (green line), with the first fault scenario.

Figure 15 .
Figure 15.Generator speed using the proposed controller (dark blue line), PID controller (red line), nominal generator speed (light blue line), and constraints (green line), with the first fault scenario.

Figure 16 .
Figure 16.Generated power using the proposed controller (dark blue line), PID controller (red line), nominal power (light blue line), and constraints (green line), with the first fault scenario.

Figure 17 .
Figure 17.Reference pitch angle using the proposed controller (dark blue line) and PID controller (red line), with the first fault scenario.

Figure 18 .
Figure 18.Induced drive train torsion angle rate using the proposed controller (dark blue line) and PID controller (red line), with the first fault scenario.

Figure 19 .
Figure 19.δP using the proposed controller (dark blue line) and PID controller (red line), with the first fault scenario.

Figure 20 .
Figure 20.δω using the proposed controller (dark blue line) and PID controller (red line), with the first fault scenario.

Figure 21 .
Figure 21.δω using the proposed controller (dark blue line) and PID controller (red line), with the first fault scenario.

Figure 22 .
Figure 22. δβ using the proposed controller (dark blue line) and PID controller (red line), with the first fault scenario.

Figure 24 .
Figure 24.Auxiliary signal in the case of fault-free (red line), pump wear (dark blue line), high air content (green line), and hydraulic leak (light blue line).

Figure 26 .
Figure 26.Second free wind speed profile.

Figure 27 .
Figure 27.Generated power using the proposed controller (dark blue line), PID controller (red line), nominal power (light blue line), and constraints (green line), in the fault-free situation, with the second wind speed sequence.

Figure 28 .
Figure 28.Generated power using the proposed controller (dark blue line), PID controller (red line), nominal power (light blue line), and constraints (green line), under the second fault scenario, with the second wind speed sequence.

Table 1 .
Wind turbine benchmark model parameters.

Table 3 .
Performance metrics in the fault-free situation.

Table 5 .
Performance metrics with the first fault scenario.

Table 8 .
Performance metrics in the fault-free situation, with the second wind speed sequence.

Table 9 .
Performance metrics under the second fault scenario, with the second wind speed sequence.