Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Awakening City: Traces of the Circadian Rhythm within the Mobile Phone Network Data
Previous Article in Journal
The Impact of Student Characteristics for Working with AR Technologies in Higher Education—Findings from an Exploratory Study with Microsoft HoloLens
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Observer-Based Active Fault-Tolerant Control of an Asymmetric Twin Wind Turbine

1
Laboratory of Mechanics, Modelling and Manufacturing (LA2MP), National Engineering School of Sfax, Sfax 3038, Tunisia
2
Quartz Laboratory, EA-7393, ENSEA, 95014 Cergy-Pontoise, France
3
LS2N, UMR CNRS 6004, ECN, 44300 Nantes, France
4
Higher Institute of Applied Science and Technology of Sousse, Sousse University, Sousse 4003, Tunisia
*
Author to whom correspondence should be addressed.
Information 2022, 13(3), 113; https://doi.org/10.3390/info13030113
Submission received: 30 December 2021 / Revised: 30 January 2022 / Accepted: 11 February 2022 / Published: 25 February 2022

Abstract

:
This paper deals with the design of an active fault-tolerant control based on observers for a twin wind turbine, consisting of two wind turbines mounted on the same tower. We consider an asymmetric conditions case, when only one turbine is affected by an inter-turn short circuit fault of a permanent magnet synchronous machine. A diagnosis design is developed which combines the fault estimation method together with an active fault-tolerant control. The main advantage of the proposed method is to detect and correct the considered fault in a short time in order that the twin wind turbine behaves as it does in the healthy case.

1. Introduction

The huge amounts of electricity consumption and limited energy supply have spearheaded the necessity of renewable energy as solar photovoltaics and wind turbines. In this work, a new concept of a twin wind turbine (TWT), which is composed of two turbines mounted on the same tower, is recalled [1]. It is intended to be used on offshore, where severe conditions lead to the appearance of faults. The topic of fault diagnosis of wind turbines during operation has been addressed increasingly in recent years. In fact, because of the non-stationary operating conditions and the variable load, wind turbines are subject to many types of faults, such as blade faults, mechanical faults [2] and electrical faults. Here we consider an electrical fault which affect only one turbine of the TWT. The inter-turn short circuit is the most common fault in permanent magnet synchronous machines (PMSM). The latter is usually related to insulation degradation, that grows and leads to devastating consequences in a short period of time, even if an early detection and diagnosis is not established. The considered fault, which influences only one turbine in our case, leads to consider the problem of control of the TWT under asymmetric conditions. For that, an active fault-tolerant control based on observers is developed. The principle of the method is analyzed at two levels: the existence of fault and its severity are determined (level 1), and the specific fault is then compensated to stabilize the system around a reference trajectory (level 2). Various works have picked up on the real-time detection of this type of failure in the literature [3,4].
A motor current signature analysis (MCSA) technique for inter-turn stator fault detection have been reported in [5]. However, the MCSA is usually used on stationary analysis, which leads to unsatisfactory results when they are applied under non-stationary conditions [6]. To analyze non-stationary conditions, the discrete wavelet transforms (DWT) and the continuous wavelet transforms (CWT) are used in [7,8] to prove the time–frequency features in the detection of a stator winding fault and a rotor imbalance fault in the wind turbine. The DWT method has a powerful filtering capability for wind turbine signals. The application of the DWT strengthens the viability of the proposed technique in detecting the changes of wind turbine running condition. Indeed, the feasibility of detecting a wind turbine mechanical fault through analyzing the generator power signal using the CWT technique has been demonstrated. The wavelet analysis is capable of detecting the stator turns fault under both stationary and non-stationary conditions. Nevertheless, it is still not able to locate the fault position and to predict the fault severity.
In [9], a diagnosis method based on fuzzy logic is applied to detect inter-turn short circuit and open phase faults in the PMSM of a wind turbine.
Apart from analyzing conventional vibration studies that are articulated on current, voltage, torque and power signals [10], model-based methods, such as parity space techniques, parameters estimation and observers methods [11] represent some of the most successful fault diagnosis approaches. In fact, observers constitute an important tool which can estimate outputs of the system even under faulty condition case. Specifically, sliding mode observers which have been vastly used in many dynamic systems, such as robotics, vehicles and motors. This method allows the identification of the fault and, furthermore, its isolation. This foremost step is called a fault detection and isolation (FDI). For a linear system, two sliding mode observers in cascade are considered in [12]. Observers-based approaches, based on sliding mode techniques, are also reviewed in [13,14] to estimate the system state variables asymptotically, in the presence of faults. For wind turbines, in [15], an observer-based fault detection and isolation FDI approach using a Kalman filter is developed for benchmark model. The same approach is used in [16,17] to detect sensor and actuator faults and to estimate blade root moment sensor faults and pitch actuator stuck faults, respectively. In [18], an observer-based fault detection approach for stator inter-turn short circuit fault in doubly fed electric machine, based on d q -frame, is established.
This paper proposes an observer-based approach for a TWT under asymmetric fault. Indeed, an active fault-tolerant control scheme based on a b c -frame is developed in order to compensate the imbalance caused by the fault and to be stable around the nominal equilibrium. As a main feature of this method, the equipment may be protected, and subsequently the TWT continues operating under nominal conditions with acceptable performances.
The paper is organized as follows. In Section 2, the TWT in both healthy and faulty cases is described. Section 3 presents our method of active fault-tolerant control. In Section 4, simulation results highlight the interest of the proposed method. Conclusion and open problems for future research are drawn in Section 5.

2. Problem Statement

The model of the TWT, patented by [1], is presented on this section in both healthy and faulty cases. The principal notations used in this paper are given by Table 1.

2.1. Healthy Twin Wind Turbine Description

As shown in Figure 1, the TWT is composed of two turbines mounted on the same tower. The advantage of this structure is that its orientation in front of the wind is ensured without needing an actuator. In fact, to keep this orientation while maintaining optimal power production requires robust control. More details on the wind turbine description can be found in [1]. The innovative concept, which is intended to be used offshore, is exposed to variable wind speeds. The mechanical power computed from the kinetic energy and the angular speed and aerodynamic torque applied to the rotor by the wind are, respectively, given by Equations (1)–(3) :
P i = π ρ 2 C p i λ i , β i R p 2 W w cos ( ψ α ) 3 ,
Ω i = λ i W w cos ψ α R p ,
Γ a i = P i Ω i = π ρ 2 λ i C p i λ i , β i R p 3 W w cos ( ψ α ) 2 ,
Remark 1.
In this paper, the index, i = 1 , 2 , denotes the i-th turbine (i.e, the first and second turbines, respectively).
The angle α represents the angle between the wind direction and the fixed direction and the angle ψ is defined as the angle between the horizontal axis of structure and the fixed direction. C p i is the power coefficient [19], characterized by nonlinear functions depend on the pitch angle of blade β i and the tip speed ratio λ i , which is proportional to the rotor angular speed (see Equation (2)).
To extract the maximum of wind energy, the TWT must be oriented face to wind and this is realized when the angle of orientation, ψ , attains its reference, α . It can be argued that this orientation is the common input of the two turbines.
As it was mentioned, the orientation of the system in front of the wind is obtained without need of an actuator. In contrast, it is made by the torque extracted from the difference between forces generated by the two turbines, which enables system rotation around its vertical axis. Thus, considering that ψ ˙ is the derivative of ψ with respect to time, dynamics of system rotation is defined by:
d r ψ ¨ = f r ψ ˙ + F 1 F 2 ,
In Equation (4), d r is the inertia moment, f r is the friction coefficient associated with yaw motion, is the distance between the horizontal and vertical axis of structure, and F 1 F 2 is the difference between the two drag forces. The expression of F i is given by the following equation:
F i = π ρ 2 C d i λ i , β i R p 2 W w cos ( ψ α ) 2 ,
where C d i is the drag coefficient, whose expression is given in the Appendix A by (A3).
To achieve optimal power production, power coefficients have to reach their optimum C p i o p t , which correspond to β i o p t and λ i o p t [20]. Based on healthy assumptions, the two turbines have the same power coefficient and the same optimal tip speed ratio,
C p 1 λ 1 o p t , β 1 o p t = C p 2 λ 2 o p t , β 2 o p t and λ 1 o p t = λ 2 o p t = λ o p t ,
So, the same optimal pitch angle of blades β 1 o p t = β 2 o p t = β o p t is adopted. Furthermore, in this concept, the pitch angles are adjusted well to maximize power and to allow the rotation of the system (drag force of each turbine depends on its pitch blade (1)). The dynamics of the pitch angles can be given by the first-order system:
β ˙ i = 1 t β ( β i r e f β i ) , i = 1 , 2 ,
where the constant t β represents the time constant of blades actuation system and β i r e f is the pitch blade reference
β 1 r e f = β o p t + Δ β and β 2 r e f = β o p t Δ β .
The electrical part of the TWT, is equipped with two salient permanent magnet synchronous machines. This type is widely used in such applications because of many features as its high efficiency and compactness, simplification of construction and maintenance and excellent dynamic response characteristics. For the healthy case, the three-phase permanent magnet synchronous machine model in a b c -frame [21] for each turbine, is represented by the following equations:
V i = L i h I i + L i h d d t I i + d d t E m i h , i = 1 , 2 ,
where L i h is given by:
L i h = r s + d L a ( θ e i ) d t d M a b ( θ e i ) d t d M a c ( θ e i ) d t d M b a ( θ e i ) d t r s + d L b ( θ e i ) d t d M b c ( θ e i ) d t d M c a ( θ e i ) d t d M c b ( θ e i ) d t r s + d L c ( θ e i ) d t .
I i = i a i i b i i c i T and V i = v a n i v b n i v c n i T are, respectively, the current and voltage vectors of the three phases related to the i-th turbine.
Assumption 1.
For safety reasons, the three-phase voltages are limited by V m a x , for some V m a x > 0 .
The electromotive force vector, E m i h = e m a i e m b i e m c i T , depends on the electrical angular position, θ e i ( d θ e i d t = Ω i ) , and on the flux linkage, ϕ , produced by the permanent magnet, where:
e m a i = ϕ cos ( θ e i ) , e m b i = ϕ cos θ e i 2 π 3 , e m c i = ϕ cos θ e i + 2 π 3 ,
The inductance matrix L i h also depends on the rotor electrical angular position, due to saliency:
L i h = L a ( θ e i ) M a b ( θ e i ) M a c ( θ e i ) M b a ( θ e i ) L b ( θ e i ) M b c ( θ e i ) M c a ( θ e i ) M c b ( θ e i ) L c ( θ e i ) ,
The inductance matrix, L i h , contains self inductance of each phase in its diagonal L a ( θ e i ) , L b ( θ e i ) and L c ( θ e i ) , and the mutual inductance, M a b ( θ e i ) , M a c ( θ e i ) and M b c ( θ e i ) . Details of inductance matrix terms are carried out by the following expressions:
L a ( θ e i ) = L + L 0 + L 1 cos ( 2 θ e i ) , M b c ( θ e i ) = M c b ( θ e i ) = M 0 + L 1 cos ( 2 θ e i ) ,
L b ( θ e i ) , = L + L 0 + L 1 cos ( 2 θ e i 2 π 3 ) , M a b ( θ e i ) = M b a ( θ e i ) = M 0 + L 1 cos ( 2 θ e i 2 π 3 ) ,
L c ( θ e i ) = L + L 0 + L 1 cos ( 2 θ e i + 2 π 3 ) , M a c ( θ e i ) = M c a ( θ e i ) = M 0 + L 1 cos ( 2 θ e i + 2 π 3 ) .
where L , L 0 , M 0 , and L 1 are constants denoting, stator leakage inductance, average value of the winding magnetizing inductance, and the magnitude of the inductance variation due to the non-uniformity of the air gap, respectively.
Remark 2.
The inductance matrix is uniformly bounded over time (see Equations (11)–(13)). In addition, it admits three eigenvalues which are different from zero and not dependent on θ e . Otherwise, the determinant is a non-zero constant and then the matrix is even invertible.
Considering that i d i , i q i v d i , v q i are, respectively, the currents and voltages of the stator for each turbine, L d and L q are the dq-axis inductance and p is the pole-pair number, in d q -frame, we have:
d i d i d t = r s L d i d i + p L q L d Ω i i q i + 1 L d v d i ,
d i q i d t = r s L q i q i p L d L q Ω i i d i p ϕ L q Ω i + 1 L q v q i .
The interconnection between the input torque (aerodynamic torque) and the output (electromagnetic torque) is ensured by the shaft of the electrical machine, and is given by the following mechanical model:
d Ω i d t = 1 J ( Γ a i Γ e m i f v Ω i ) ,
where J is the total inertia, f v is the viscous friction and Γ e m i is the electromagnetic torque, given by:
Γ e m i = p L d L q i d i i q i + p ϕ i q i .
Equations (14), (16) and (17) are used later in the health-control design. As previously explained, the rotation of the TWT is ensured without need of an actuator. Consequently, a control strategy must be developed in order to achieve three important objectives. To produce the maximum of electrical production, the first objective is to keep the whole structure in front of the wind to extract the maximum of energy. Secondly, the rotational speeds of the two wind turbines are controlled at a reference Ω i r e f , by keeping their tip–speed ratios at their optimal values. Furthermore, finally, as the fatigues loads in the mechanical shaft caused by the ripples of the electromagnetic torque can impact the produced electrical power, one solution consists of forcing the direct current i d i to zero. To sum up, objectives on the healthy case can be reaching by forcing the output, y, to zero, where Y is given by:
Y = ψ α i d 1 Ω 1 Ω r e f i d 2 Ω 2 Ω r e f T .

2.2. Faulty Twin Wind Turbine Description

Taking into account electrical fault, inter-turn short circuit failure in the a-phase has been studied in [20] to test the robustness of the control, based on the d q -frame. The same fault has also been examined with a different severity in [22] to develop an active fault-tolerant control, based on the a b c -frame.
In this paper, we assume that this fault persists on the three stator phases in only one turbine. For the healthy case, where the three stator phases are symmetric and have the same total number of winding N t , the severity can be defined as the factor of the short-circuited turns N a d , N b d and N c d of the a b c -phases devising by the total number of winding. With respect to the severities, δ a , δ b and δ c of a-phase, b-phase and c-phase, respectively, the model can be written in the a b c -frame as:
V i = L i f I i + L i f d d t I i + d d t E m i f ,
where the matrix L i f is given by:
L i f = 1 δ a r s + d L a ( θ e i ) d t 1 δ a + δ b d M a b ( θ e i ) d t 1 δ a + δ c d M a c ( θ e i ) d t 1 δ b + δ a d M b a ( θ e i ) d t 1 δ b r s + d L b ( θ e i ) d t 1 δ b + δ c d M b c ( θ e i ) d t 1 δ c + δ a d M c a ( θ e i ) d t 1 δ c + δ b d M c b ( θ e i ) d t 1 δ c r s + d L c ( θ e i ) d t ,
The electromagnetic force vector is given by
E m i f = 1 δ a e m a i 1 δ b e m b i 1 δ c e m c i T ,
and the faulty inductance matrix L i f is given by:
L i f = L i h δ a L a ( θ e i ) M a b ( θ e i ) M a c ( θ e i ) M b a ( θ e i ) 0 0 M c a ( θ e i ) 0 0 δ b 0 M a b ( θ e i ) 0 M b a ( θ e i ) L b ( θ e i ) M b c ( θ e i ) 0 M c b ( θ e i ) 0 δ c 0 0 M a c ( θ e i ) 0 0 M b c ( θ e i ) M c a ( θ e i ) M c b ( θ e i ) L c ( θ e i ) .
Remark 3.
For a healthy machine, the inductance has the same expression as the fault inductance just by replacing the three severities of faults with zero δ a = 0 , δ b = 0 , δ c = 0 .
Assumption 2.
The number of shorted turns is smaller than 50 % of the total number of turns per phase. This implies that δ a , δ b and δ c 0 , 0.5 .
As the fault cannot exceed the 50 % (if not, then the faulty electrical machine will be drastically damaged and thus the whole system will be influenced by the defect), the time-varying inductance matrix L i f is always invertible. In the considered case, the inverse of the matrix is carried out numerically by reason of complicated equations.

2.3. Control Strategy

Concerning control strategies, one is used in [23], in the healthy case, based on the high order sliding mode (HOSM) and modified in [20], when considering the faulty case (an insulation degradation fault persists on only one turbine), based on a passive fault-tolerant control. The limit of this latter control is that it can be exploited only for incipient faults. Thus, if the defect exceeds 8 % , the objectives of the control are not achieved.
Nevertheless, the appearance of the fault break reflects not only the symmetry of the phases but also the symmetry of the two turbines. Consequently, if the severity of the fault is important, the stability of the TWT is disturbed, which highlights the requirement of an efficient fault detection method. Hence, an active fault-tolerant control is developed.
The idea is once the fault exceeds 8 % , the system switches to the proposed observer-based active fault-tolerant control method. This latter will be able to pinpoint the inter turn short-circuit fault quickly and precisely. Concerning control strategies, one is used in [23], in the healthy case, based on high order sliding mode and modified in [20], when considering the faulty case (an insulation degradation fault persists on only one turbine), based on a passive fault-tolerant control. The limit of this latter control is that it can be exploited only for incipient faults. Thus, if the defect exceeds 8 % , then the objectives of the control are not achieved.
Nevertheless, the appearance of the fault breaks not only the symmetry of phases but also the symmetry of the two turbines. Consequently, if the severity of the fault is important, the stability of the TWT is disturbed, which highlights the requirement of an efficient fault detection method. Hence, an active fault-tolerant control is developed.
The idea is once the fault exceeds 8 % , the system switches to the proposed observer-based active fault-tolerant control method. This latter will be able to pinpoint the inter-turn short circuit fault quickly and precisely.

3. Active Fault-Tolerant Control and Diagnosis Method

3.1. Active Fault-Tolerant Control

In this section, the active fault-tolerant control, “when the fault is known”, is developed. Considering all faults equal to zero, the proposed control is also valid for the healthy system. Taking into consideration an inter-turn short circuit on the three phases of PMSM in only one turbine, the control based on d q -rotating frame will not be adequate. So, a new control strategy based on a b c -frame is adapted by depending on the direct current—indeed the homopolar component ( i h i = 1 3 ( i a i + i b i + i c i ) )—which must be assigned to zero in order to respect the Kirchhoff laws. This approach leads to two additional states in the state vector ( i a 1 , i b 1 , i c 1 , i a 2 , i b 2 , i c 2 three-phase currents on each turbine are considered instead of two-phase currents i d 1 , i q 1 , i d 2 , i q 2 ).
In the same way, two extra inputs are introduced ( v a n 1 , v b n 1 , v c n 1 , v a n 2 , v b n 2 , v c n 2 instead of v d 1 , v q 1 , v d 2 , v q 2 ). Therefore, an extended nonlinear system based on a b c -frame is given by:
x ˙ = f ( x , d ) + g ( x , d ) u ,
where expressions of f ( x , d ) and g ( x , d ) are given in the Appendix A (Equations (A1) and (A2)), the state vector x and the input vector u are given by:
x = β 1 β 2 ψ ψ ˙ i a 1 i b 1 i c 1 Ω 1 i a 2 i b 2 i c 2 Ω 2 T ,
u = Δ β v a n 1 v b n 1 v c n 1 v a n 2 v b n 2 v c n 2 T .
Remark 4.
As Δ β represents the difference between the two blade pitch angles, which is limited by 2 π , and the three-phase voltage, V i , which is also limited by V m a x , the input vector u in our case is bounded.
Comparing to [22], the active fault-tolerant control has 7 objectives to achieve instead of 5. In fact, the two additional outputs i h 1 = 1 3 ( x 5 + x 6 + x 7 ) and i h 2 = 1 3 ( x 9 + x 10 + x 11 ) are used to correct the asymmetry of phases and thus the expression of the new outputs vector:
y = h ( x ) = x 3 α i d 1 x 8 Ω 1 r e f i h 1 i d 2 x 12 Ω 2 r e f i h 2 T ,
Note that i d 1 and i d 2 depend, respectively, on ( x 5 , x 6 , x 7 ) and ( x 9 , x 10 , x 11 ) , with respect to Park’s transformation. The successive derivatives of the outputs lead to the following equation:
y ( ε ) = Λ ( x , d ) + Θ ( x , d ) u ( x , d ) ,
where y ( ε ) = ( y 1 ( 3 ) , y 2 ( 1 ) , y 3 ( 2 ) , y 4 ( 1 ) , y 5 ( 1 ) , y 6 ( 2 ) , y 7 ( 1 ) ) T , with y ( i ) = d i y d t i , for i 1 . The control is developed based on high order sliding mode, for that each output has its relative degree ε = ( ε 1 , ε 2 , ε 3 , ε 4 , ε 5 , ε 6 , ε 7 ) = ( 3 , 1 , 2 , 1 , 1 , 2 , 1 ) related to ( y 1 , y 2 , y 3 , y 4 , y 5 , y 6 , y 7 ) , respectively.
The vector field Λ is given in the Appendix A (A5) and the decoupling matrix Θ is given by:
Θ ( x , d ) = Θ 1 0 Θ 2 0 Θ 3 Θ 1 = 2 ρ π W w cos x 3 α 2 d r T β B λ Θ 2 = P i L i f 1 Θ 3 = P i L i h 1 ,
where P i is the Park’s transformation matrix. The decoupling matrix, Θ , is a square block diagonal matrix in which the diagonal elements are three square matrices of ( 1 × 1 , 3 × 3 and 3 × 3 ) and the off-diagonal elements are 0. The second block is the product of the Park transformation matrix and the inverse of inductance matrix of the faulty turbine L i f and the third block is the product of the Park transformation matrix and the inverse of inductance matrix of the healthy one L i h . As a result, this matrix is regular, which leads to a following decoupling control based on a new control ϑ (which will be detailed later):
u ( x , d ) = Θ ( x , d ) 1 ϑ Λ ( x , d ) .
Diffeomorphism
From the previously presented control scheme, we define the diffeomorphism Z = h ( x ) , which given as follows z 1 = y 1 , z 2 = y ˙ 1 , z 3 = y ¨ 1 , z 4 = y 2 , z 5 = y 3 , z 6 = y ¨ 3 , z 7 = y 4 , z 8 = y 5 , z 9 = y 6 , z 10 = y ˙ 6 , z 11 = y 7 and z 12 = β 1 + β 2 .
z 1 = x 3 α z 2 = x 4 z 3 = f r d r x 4 + 2 π l ρ R p 2 d r B ( λ ) W w 2 cos ( x 3 α ) 2 ( x 1 x 2 ) z 4 = 2 3 cos θ e x 5 + cos θ e 2 π 3 x 6 + cos θ e + 2 π 3 x 7 z 5 = x 8 Ω 1 r e f z 6 = 1 J Γ a 1 Γ e m 1 f v x 8 z 7 = 1 3 x 5 + x 6 + x 7 z 8 = 2 3 cos θ e x 9 + cos θ e 2 π 3 x 10 + cos θ e + 2 π 3 x 11 z 9 = x 12 Ω 2 r e f z 10 = z 6 = 1 J Γ a 1 Γ e m 1 f v x 12 z 11 = 1 3 x 9 + x 10 + x 11 z 12 = x 1 + x 2
Knowing that k = 1 7 ε k = 11 < 12 , there is a zero dynamics of dimension 1, which is expressed as:
z ˙ 12 = 2 β r e f t β 1 t β z 12 .
The zero dynamics depends on the external input which is not directly influenced by the defect. For that, the dynamic of zero, z 12 , is always input-to-state stable (input: β r e f ) [24].
Definition 1.
Consider an integer, ϵ 2 , and for j = 2 , ϵ , define the vector, Z j = [ z 1 , z j ] T , the decreasing sequence of positive real numbers, r j = r 1 ( j 1 ) , for some r 1 2 , the non decreasing sequence of positive real numbers α j , so that α j α j 1 r 1 . Defining also, for j = 2 , ϵ , the C 1 r-homogeneous function:
σ j = z j α j r j + k j 1 α j r j σ j 1 α i α j 1
where c being a function defined as . c . s g n ( . ) , the index σ 1 = z 1 α 1 r 1 and k j 0 .
Based on the definition and the diffeomorphism, the proposed controller in this paper reads as:
ϑ = K y σ y μ j , k = K ψ σ ψ μ 1 , 3 K i d z 4 μ 4 K Ω σ Ω 1 μ 5 , 6 K i h z 7 μ 7 K i d z 8 μ 8 K Ω σ Ω 2 μ 9 , 10 K i h z 11 μ 11 ,
where σ p s i , σ Ω 1 , σ Ω 2 can be calculated based on Definition 1 and the parameter μ reads as the following adaptive law:
μ j , k = max { 1 γ j k | z j | | z j | + c , 0 } . 1 j , k 11 , c > 0 , γ > 1
The advantage of the active fault-tolerant control, that it makes the machine behave as in the case where the three phases are symmetric and that way it enables the fault to be compensated. In this section, we assume that the fault is known. Because this is not a realistic case, we will treat the global problem in the next sections.

3.2. Fault Estimation

In this section, a fault detection method is developed in order to estimate the defect, d. Coupled with the control designed in Section 3.2, the fault detection method gives the possibility to design an active fault-tolerant control in order to achieve the reference trajectory. The global fault vector, d, that contains the three-phase defects is carried out by:
d = δ a δ b δ c T .
The monitoring sensors in the faulty model are x 5 , x 6 and x 7 (which correspond, respectively, to the three-phase currents i a , i b and i c ).
The idea here is to reorganize terms which are expressed linearly in function of the defect in the right hand and which are independent on the fault in the left. Otherwise, from Equations (7) and (19), and knowing that the defect is modeled in a linear way under the following form:
L i f = L i h Δ L ; Δ L = δ a r s + d L a ( θ e i ) d t δ a + δ b d M a b ( θ e i ) d t δ a + δ c d M a c ( θ e i ) d t δ b + δ a d M b a ( θ e i ) d t δ b r s + d L b ( θ e i ) d t δ b + δ c d M b c ( θ e i ) d t δ c + δ a d M c a ( θ e i ) d t δ c + δ b d M c b ( θ e i ) d t δ c r s + d L c ( θ e i ) d t L i f = L i h Δ L ; Δ L = δ a L a ( θ e i ) δ a + δ b M a b ( θ e i ) δ a + δ c M a c ( θ e i ) δ b + δ a M b a ( θ e i ) δ b L b ( θ e i ) δ b + δ c M b c ( θ e i ) δ c + δ a M c a ( θ e i ) δ c + δ b M c b ( θ e i ) δ c L c ( θ e i ) E m i f = E m i h Δ E m ; Δ E m = δ a ϕ cos θ e i δ b ϕ cos θ e i 2 π 3 δ c ϕ cos θ e i + 2 π 3 T
We can obtain relation (34) between the measurement vector, γ , and the fault vector, d, which takes the form γ = D ( x ) d .
V i L i h I i + L i h d d t I i + d d t E m i h γ = Δ L I i + Δ L d d t I i + Δ E m , D ( x ) d
where D(x) is given by:
D x = D 11 R L i h ( 1 , 2 ) i b i + L i h ( 1 , 2 ) d i b i d t R L i h ( 1 , 3 ) i c i + L i h ( 1 , 3 ) d i c i d t R L i h ( 2 , 1 ) i a i + L i h ( 2 , 1 ) d i a i d t D 22 R L i h ( 2 , 3 ) i c i + L i h ( 2 , 3 ) d i c i d t R L i h ( 3 , 1 ) i a i + L i h ( 3 , 1 ) d i a i d t R L i h ( 3 , 2 ) i b i + L i h ( 3 , 2 ) d i b i d t D 33 ,
For D k k = j = 1 3 R L i h ( k , j ) I i ( j ) + L i h ( k , j ) d I i ( j ) d t + E m i ( k ) . f o r k = 1 , 2 , 3
Note that R L i h ( k , j ) refers to the element of the row, k, and the column, j, is related to the matrix, R L i h .
Remark 5.
If the matrix D ( x ) is regular, it is possible at the same time to solve the problem of detection, isolation and estimation of the faults.
Assumption 3.
The matrix D ( x ) is always invertible.
It is proved, numerically (because of the equation’s complexity), that D ( x ) is invertible. In our case, the number of measurements is equal to the number of unknown uncertainties, which can solve the left inverse problem by inverting directly the matrix D ( x ) in (34).
In this work, the estimation method is stable and the active fault-tolerant control is also stable. Even so, what can we say about the global stability in the case of the coupling of the two tools?

3.3. Observer-Based Active Fault-Tolerant Control

The diagram in Figure 2 provides a comprehensive summary of inputs and outputs of the system, the concept of the proposed technique and the entire system, which is made up of two turbines: a healthy turbine and a faulty turbine. After verifying the stability of both control strategy and the fault estimation method, the two strategies are coupled in order to detect the fault and to correct it at the same time.
The diagnosis method focus only on one part of the global system, the dynamics of observers—indeed the global model—are carried out by the following system:
x ˙ = f ( x , d ) + g ( x , d ) u ( x , d ^ ) d ^ ˙ = d d ^ 0.5 s i g n d d ^
The power coefficient 0.5 is used in order to reduce sensitivity to noise [25]. Nevertheless, if it is equal to 0, then the dynamical system becomes a first order sliding mode, and the chattering phenomenon occurs at the equilibrium point. Furthermore, if this coefficient is equal to 1, then the dynamical system becomes linear.
Remark 6.
In (37), d comes from the algebraic resolution of (34), if we consider that this resolution is carried out without computation and measurement noise, then we can use directly d in the active fault-tolerant control, and there is no more stability problem due to the fact that d ^ is not always equal to d in the control.
Theorem 1.
Under Assumptions 1, A2 and A3, the closed-loop system (37) is exponentially stable.
Proof of Theorem 1. 
Note that, since the zero dynamic (i.e., (30)) is only exponentially stable and does not depend on the fault d, the stability is proved only exponentially even if all the other states are finite time stable.
Now, considering d ^ = d , the first equation of the system (37) becomes:
x ˙ = f ( x , d ) + g ( x , d ) u ( x , d ) ,
and this gives y ( ε ) = ϑ , which is finitely time stable, see [26]. Moreover, the dynamic d ^ ˙ = d d ^ 0.5 s i g n d d ^ is also finitely time stable.
Consequently, it is only sufficient to verify that during the time d ^ is different to d, the state x remains bounded. This is equivalent to prove that the system (39) is locally ISS with respect to Δ u .
x ˙ = f ( x , d ) + g ( x , d ) u ( x , d ^ ) = f ( x , d ) + g ( x , d ) u ( x , d ) + g ( x , d ) Δ u .
where Δ u = u ( x , d ^ ) u ( x , d ) .
Note that Δ u is bounded by Assumption 1 and g ( x , d ) is bounded under the Assumption 2 (which is also detailed in the Appendix A). The local ISS property of the previous equation follows from the fact that, when d ^ = d ( Δ u = 0 ), the system (38) (or equivalently (29) in Z) is finitely time stable (and therefore at least locally exponentially stable) together with the fact that g ( x , d ) Δ u is bounded (see [27] for more details).
This ends the proof of the theorem. □

4. Simulation Results

The faulty TWT model extracted in the problem statement has been applied to evaluate the performances of the proposed diagnosis method. In this section, simulation results, using Ode1 solver in Matlab Simulink (Mathworks Inc., Natick, MA, USA) with a fixed time step of 10 4 s, are performed and commented.
To quantify the influence of the faults on the outputs of the system, an inter-turn short circuit fault is considered on one phase with a severity of 4 % , which occurs after 7 s. A passive control based on HOSM is used in this case [20]. As can be seen in Figure 3a, once the fault appears, the three-phase currents of the faulty machine increase and the direct current is not well controlled to be assigned to zero (Figure 3b). This is because the fault breaks the symmetry of the three phases of the machine. In Figure 4, the angular speed of the healthy machine follows its reference (the blue and the green curves are confused), whereas the angular speed of the faulty machine does not reach its reference. Therefore, if the fault exceeds 8 % , then a new control method is indispensable to compensate the fault and to keep the TWT in front of the wind. For that, an active fault-tolerant control is proposed.
To prove the effectiveness of the proposed method, three faults are considered, respectively, in the a-phase, b-phase and c-phase of the PMSM of only one turbine, with important severities which are greater than 8 % ( 30 % , 16 % and 10 % of turns are short-circuited, respectively, in the a-phase, b-phase and c-phase).
From Figure 5, the first fault is activated at the time 3 s, which affects the a-phase, while the other faults are set to zero.
After 3 s, 30 % of the totality of turns is defected. At time 5 s, the second fault is started and then after 8 s the third fault is setting off, which leads to the existence of three faults simultaneously. For the second fault, 16 % of turns are defected, whereas only 10 % of turns are defected on the c-phase (see Figure 6).
In this case, a diagnosis method has been carried out: firstly, we are able to detect the presence of faults, secondly, to identify and distinguish them from the obtained signals and finally, the active fault-tolerant control has been implemented in order to ensure a well-functioning of the system. With the presence of faults, the three-phase currents are subject to increase and generate an uncontrolled direct current, i d i . With the AFTC, these currents are well guided to attain their references, as shown in Figure 7.

5. Conclusions

In this paper, an active fault-tolerant control based on fault estimation was proposed for a new concept of a TWT. Numerical simulations performed in MATLAB/Simulink proved that the proposed AFTC design can ensure the stability of the electrical machine currents output with and without a fault. In addition, the proposed approach makes it possible to keep the whole system oriented to face toward the wind, despite the existence of a fault on only one turbine. Moreover, the fault detection approach detailed in this paper can be used for many other faults affecting wind turbines. For that, further work will be required to apply the sparse recovery diagnosis method, while considering other types of faults related to the operating conditions which the turbine is connected to.

Author Contributions

Methodology, I.H., J.-P.B. and F.P.; Supervision, I.H. and J.-P.B.; Validation, I.H. and J.-P.B.; Visualization, N.F. and M.S.A.; Writing—original draft, M.M.; Writing—review & editing, I.H. and J.-P.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors decalre no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
TWTTwin Wind Turbine
PMSMPermanent Magnet Synchronous Machine
MCSAMotor Current Signature Analysis
DWTDiscrete Wavelet Transforms
CWTContinuous Wavelet Transforms
FDIFault Detection and Isolation
AFTCActive Fault-Tolerant Control
ISSInput to-State Stable
HOSMHigh Order Sliding Mode

Appendix A

The vector f ( x , d ) and the matrix g ( x , d ) presented on the nonlinear system (22) are expressed as:
f ( x , d ) = β 1 o p t x 1 T β β 2 o p t x 2 T β x 4 1 d r f r x 4 + F 1 x 1 , x 3 F 2 x 2 , x 3 L L 1 f 1 L 1 f x 5 x 6 x 7 T + d d t E m 1 f 1 J Γ a 1 x 1 , x 3 Γ e m 1 x 5 , x 6 , x 7 f v x 8 L 2 h 1 L 2 h x 9 x 10 x 11 T + d d t E m 2 h 1 J Γ a 1 x 2 , x 3 Γ e m 1 x 9 , x 10 , x 11 f v x 12 ,
g ( x , d ) = 1 T β 0 0 0 0 0 0 1 T β 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 L 1 f 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 L 2 h 1 0 0 0 0 0 0 0 0 0 0 0 .
As T β is a constant, the two matrix L 1 f and L 2 h are also bounded and always invertible. This implies that the matrix g ( x , d ) is bounded by the inverse of their maximum max L 1 f , L 2 h = L l + L 0 + L 1 .
The expression of the drag coefficient c d i is given by:
c d i λ i , β i = A i λ i + B i λ i β i ,
where
A i λ i = a 0 + a 1 λ i + a 2 λ i 2 + a 3 λ i 3 and B i λ i = b 0 + b 1 λ i + b 2 λ i 2 + b 3 λ i 3 .
The constants a k and b k for k = 1 , 2 , 3 are equal, respectively, to: a 0 = 0.25382 , a 1 = 0.1369 , a 2 = 0.04345 , a 3 = 0.00263 , b 0 = 0.008608 , b 1 = 0.0063 , b 2 = 0.0015 and b 3 = 0.000118 .
Considering that C = ρ π 2 R p 2 W w 2 cos 2 ψ α , and with repect to the fault, the vector field Λ ( x , d ) is given by:
Λ ( x , d ) = f r d r ψ ¨ + B λ C d r T β l β 1 β 2 + B λ C ˙ d r l β 1 β 2 + B ˙ λ C d r l β 1 β 2 Λ 2 1 J Γ ˙ a 1 Γ ˙ e m 1 1 3 f 5 + f 6 + f 7 Λ 5 1 J Γ ˙ a 2 Γ ˙ e m 2 1 3 f 9 + f 10 + f 11 ,
with
Λ 2 = 2 3 cos θ e 1 cos θ e 1 2 π 3 cos θ e 1 + 2 π 3 f 5 f 6 f 7 T p x 8 sin θ e 1 sin θ e 1 2 π 3 sin θ e 1 + 2 π 3 x 5 x 6 x 7 T , Λ 5 = 2 3 cos θ e 2 cos θ e 2 2 π 3 cos θ e 2 + 2 π 3 f 9 f 10 f 11 T p x 12 sin θ e 2 sin θ e 2 2 π 3 sin θ e 2 + 2 π 3 x 9 x 10 x 11 T .
The derivative of the two electromagnetic torques are given, respectively, by these expressions:
Γ ˙ e m 1 = σ 21 Λ 2 + σ 22 Λ 21 and Γ ˙ e m 2 = σ 51 Λ 5 + σ 52 Λ 51 . with
Λ 21 = 2 3 sin θ e 1 sin θ e 1 2 π 3 sin θ e 1 + 2 π 3 f 5 f 6 f 7 T p x 8 cos θ e 1 cos θ e 1 2 π 3 cos θ e 1 + 2 π 3 x 5 x 6 x 7 T ,
Λ 51 = 2 3 sin θ e 2 sin θ e 2 2 π 3 sin θ e 2 + 2 π 3 f 9 f 10 f 11 T p x 12 cos θ e 2 cos θ e 2 2 π 3 cos θ e 2 + 2 π 3 x 9 x 10 x 11 T ,
σ 21 = p ϕ f + 2 3 p L d L q cos θ e 1 cos θ e 1 2 π 3 cos θ e 1 + 2 π 3 x 5 x 6 x 7 T ,
σ 22 = 2 3 p L d L q sin θ e 1 sin θ e 1 2 π 3 sin θ e 1 + 2 π 3 x 5 x 6 x 7 T ,
σ 51 = p ϕ f + 2 3 p L d L q cos θ e 2 cos θ e 2 2 π 3 cos θ e 2 + 2 π 3 x 9 x 10 x 11 T ,
σ 52 = 2 3 p L d L q sin θ e 2 sin θ e 2 2 π 3 sin θ e 2 + 2 π 3 x 9 x 10 x 11 T .
The park transformation matrix P i is given by:
P i = 2 3 cos θ e i cos θ e i 2 π 3 cos θ e i + 2 π 3 sin θ e i sin θ e i 2 π 3 sin θ e i + 2 π 3 1 2 1 2 1 2 .

References

  1. Herskovits, A.; Laffitte, O.; Thome, P.; Tobie, A. V-Shaped, bi-Rotor Wind Generator on a Spar Floating Structure. French Patent WO2014060420 A, 24 April 2014. Volume 1. [Google Scholar]
  2. Ghane, M.; Nejad, A.R.; Blanke, M.; Gao, Z.; Moan, T. Statistical fault diagnosis of wind turbine drivetrain applied to a 5 MW floating wind turbine. J. Phys. Conf. Ser. 2016, 753, 052017. [Google Scholar] [CrossRef] [Green Version]
  3. Liu, L. Robust Fault Detection and Diagnosis for Permanent Magnet Synchronous Motors. Ph.D. Thesis, The Florida State University, Tallahassee, FL, USA, 2006. [Google Scholar]
  4. Otava, L.; Graf, M.; Buchta, L. Interior Permanent Magnet Synchronous Motor Stator Winding Fault Modelling. IFAC-PapersOnLine 2015, 48, 324–329. [Google Scholar] [CrossRef]
  5. Nandi, S.; Toliyat, H.A.; Li, X. Condition monitoring and fault diagnosis of electrical motors—A review. IEEE Trans. Energy Convers. 2005, 20, 719–729. [Google Scholar] [CrossRef]
  6. Henao, H.; Capolino, G.A.; Fernandez-Cabanas, M.; Filippetti, F.; Bruzzese, C.; Strangas, E.; Pusca, R.; Estima, J.; Riera-Guasp, M.; Hedayati-Kia, S. Trends in fault diagnosis for electrical machines: A review of diagnostic techniques. IEEE Ind. Electron. Mag. 2014, 8, 31–42. [Google Scholar] [CrossRef]
  7. Yang, W.; Tavner, P.; Wilkinson, M. Condition monitoring and fault diagnosis of a wind turbine with a synchronous generator using wavelet transform. In Proceedings of the Conference on Power Electronics, Machines and Drives, York, UK, 2–4 April 2008. [Google Scholar]
  8. Yang, W.; Tavner, P.; Wilkinson, M. Condition monitoring and fault diagnosis of a wind turbine synchronous generator drive train. IET Renew. Power Gener. 2008, 3, 1–11. [Google Scholar] [CrossRef]
  9. Bechkaoui, A.; Ameur, A.; Bouras, S.; Hadjadj, A. Open-circuit and inter-turn short-circuit detection in PMSG for wind turbine applications using fuzzy logic. Energy Procedia 2015, 74, 1323–1336. [Google Scholar] [CrossRef] [Green Version]
  10. Basak, D.; Tiwari, A.; Das, S. Fault diagnosis and condition monitoring of electrical machines—A Review. In Proceedings of the 2006 IEEE International Conference on Industrial Technology, Singapore, 16–18 August 2006; pp. 3061–3066. [Google Scholar]
  11. Garcia, E.A.; Frank, P.M. Deterministic nonlinear observer-based approaches to fault diagnosis: A survey. Control Eng. Pract. 1997, 5, 663–670. [Google Scholar] [CrossRef]
  12. Ng, K.Y.; Tan, C.P.; Edwards, C.; Kuang, Y.C. New results in robust actuator fault reconstruction for linear uncertain systems using sliding mode observers. Int. J. Robust Nonlinear Control. IFAC-Affil. J. 2007, 17, 1294–1319. [Google Scholar] [CrossRef]
  13. Yan, X.G.; Edwards, C. Adaptive sliding-mode-observer-based fault reconstruction for nonlinear systems with parametric uncertainties. IEEE Trans. Ind. Electron. 2008, 55, 4029–4036. [Google Scholar]
  14. Alwi, H.; Edwards, C.; Tan, C.P. Fault Detection and Fault-Tolerant Control Using Sliding Modes; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  15. Chen, W.; Ding, S.X.; Haghani, A.; Naik, A.; Khan, A.Q.; Yin, S. Observer-based FDI schemes for wind turbine benchmark. IFAC Proc. Vol. 2011, 44, 7073–7078. [Google Scholar] [CrossRef]
  16. Wei, X.; Verhaegen, M. Fault detection of large scale wind turbine systems: A mixed H∞/H- index observer approach. In Proceedings of the 16th Mediterranean Conference on Control and Automation, Ajaccio, France, 25–27 June 2008; pp. 1675–1680. [Google Scholar]
  17. Wei, X.; Verhaegen, M. Sensor and actuator fault diagnosis for wind turbine systems by using robust observer and filter. Wind Energy 2011, 14, 491–516. [Google Scholar] [CrossRef]
  18. Odgaard, P.F.; Stoustrup, J. Unknown input observer based scheme for detecting faults in a wind turbine converter. IFAC Proc. Vol. 2009, 42, 161–166. [Google Scholar] [CrossRef] [Green Version]
  19. Guenoune, I.; Plestan, F.; Chermitti, A.; Evangelista, C. Modeling and robust control of a twin wind turbines structure. Control Eng. Pract. 2017, 69, 23–35. [Google Scholar] [CrossRef]
  20. Makni, M.; Haidar, I.; Barbot, J.P.; Plestan, F.; Feki, N.; Abbes, M.S. Analysis and control of Twin Wind Turbine subject to asymmetric fault. In Proceedings of the IEEE Conference on Control Technology and Applications (CCTA), Montreal, QC, Canada, 24–26 August 2020; pp. 1–6. [Google Scholar]
  21. Pillay, P.; Krishnan, R. Modeling of permanent magnet motor drives. In Motor Control and Power Electronics, Proceedings of the IECON’87, Cambridge, MA, USA, 4–6 November 1987; International Society for Optics and Photonics: Bellingham, WA, USA, 1987; Volume 854, pp. 289–293. [Google Scholar]
  22. Makni, M.; Haidar, I.; Barbot, J.P.; Plestan, F.; Feki, N.; Abbes, M.S. Active fault tolerant control for twin wind turbine subject to asymmetric fault. arXiv 2021, arXiv:2108.07238. [Google Scholar]
  23. Tahoumi, E.; Evangelista, C.; Plestan, F.; Ghanes, M.; Barbot, J.; Puleston, P. Energy efficient control derived from homogeneous algorithm—Application to a wind system. Control Eng. Pract. 2020, 103, 104586. [Google Scholar] [CrossRef]
  24. Khalil, H.K. High-Gain Observers in Nonlinear Feedback Control; SIAM: Philadelphia, PA, USA, 2017. [Google Scholar]
  25. Yu, L.; Zheng, G.; Barbot, J.P. Dynamical sparse recovery with finite-time convergence. IEEE Trans. Signal Process. 2017, 65, 6146–6157. [Google Scholar] [CrossRef] [Green Version]
  26. Cruz-Zavala, E.; Moreno, J.A. Homogeneous high order sliding mode design: A Lyapunov approach. Automatica 2017, 80, 232–238. [Google Scholar] [CrossRef]
  27. Bernuau, E.; Polyakov, A.; Efimov, D.; Perruquetti, W. Verification of ISS, iISS and IOSS properties applying weighted homogeneity. Syst. Control Lett. 2013, 62, 1159–1167. [Google Scholar] [CrossRef] [Green Version]
Figure 1. SEREO concept—face (left) and top (right) views.
Figure 1. SEREO concept—face (left) and top (right) views.
Information 13 00113 g001
Figure 2. Diagnosis and fault-tolerant control diagram.
Figure 2. Diagnosis and fault-tolerant control diagram.
Information 13 00113 g002
Figure 3. The effect of an inter-turn short circuit fault of severity 4 % at time 7 s before the AFTC: (a) the three-phase currents of the faulty machine and (b) the direct current of the faulty machine.
Figure 3. The effect of an inter-turn short circuit fault of severity 4 % at time 7 s before the AFTC: (a) the three-phase currents of the faulty machine and (b) the direct current of the faulty machine.
Information 13 00113 g003
Figure 4. Comparison between angular speeds of the healthy and the faulty machines before the AFTC.
Figure 4. Comparison between angular speeds of the healthy and the faulty machines before the AFTC.
Information 13 00113 g004
Figure 5. Diagnosis of inter-turn short circuit fault in the electrical machine of TWT: 30 % of a-phase turns are defected at t = 3 s.
Figure 5. Diagnosis of inter-turn short circuit fault in the electrical machine of TWT: 30 % of a-phase turns are defected at t = 3 s.
Information 13 00113 g005
Figure 6. Diagnosis of inter-turn short circuit fault in the electrical machine of TWT: (a) severity = 16 % of: b-phase at t = 5 sand (b) severity = 10 % of: c-phase at t = 8 s.
Figure 6. Diagnosis of inter-turn short circuit fault in the electrical machine of TWT: (a) severity = 16 % of: b-phase at t = 5 sand (b) severity = 10 % of: c-phase at t = 8 s.
Information 13 00113 g006
Figure 7. The effectiveness of the active fault-tolerant control on the three-phase currents of faulty PMSM.
Figure 7. The effectiveness of the active fault-tolerant control on the three-phase currents of faulty PMSM.
Information 13 00113 g007
Table 1. Principal notations.
Table 1. Principal notations.
NameDescriptionNameDescription
W w wind speed L h inductance matrix in healthy case
R p blade radius L f inductance matrix in faulty case
ρ air density E m i electromagnetic force
β pitch blade Ω i r e f reference of the i-th angular speed
λ tip speed ratio δ severity of the fault
r s stator resistancedfault vector
L f h is the Lie derivative of h with respect to f [ f , g ] = L f g L g f is the Lie bracket of f and g
ref is used for references.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Makni, M.; Haidar, I.; Barbot, J.-P.; Plestan, F.; Feki, N.; Abbes, M.S. Observer-Based Active Fault-Tolerant Control of an Asymmetric Twin Wind Turbine. Information 2022, 13, 113. https://doi.org/10.3390/info13030113

AMA Style

Makni M, Haidar I, Barbot J-P, Plestan F, Feki N, Abbes MS. Observer-Based Active Fault-Tolerant Control of an Asymmetric Twin Wind Turbine. Information. 2022; 13(3):113. https://doi.org/10.3390/info13030113

Chicago/Turabian Style

Makni, Mariem, Ihab Haidar, Jean-Pierre Barbot, Franck Plestan, Nabih Feki, and Mohamed Slim Abbes. 2022. "Observer-Based Active Fault-Tolerant Control of an Asymmetric Twin Wind Turbine" Information 13, no. 3: 113. https://doi.org/10.3390/info13030113

APA Style

Makni, M., Haidar, I., Barbot, J. -P., Plestan, F., Feki, N., & Abbes, M. S. (2022). Observer-Based Active Fault-Tolerant Control of an Asymmetric Twin Wind Turbine. Information, 13(3), 113. https://doi.org/10.3390/info13030113

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop