Next Article in Journal
Insight into the Influence of Grinding on the Extraction Efficiency of Selected Bioactive Compounds from Various Plant Leaves
Next Article in Special Issue
Power Flow Analysis on the Dual Input Transmission Mechanism of Small Wind Turbine Systems
Previous Article in Journal
An Open Source Framework Approach to Support Condition Monitoring and Maintenance
Previous Article in Special Issue
Stochastic Extreme Wind Speed Modeling and Bayes Estimation under the Inverse Rayleigh Distribution
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Finite-Time Fast Dynamic Terminal Sliding Mode Maximum Power Point Tracking Control Paradigm for Permanent Magnet Synchronous Generator-Based Wind Energy Conversion System

1
Department of Electrical and Computer Engineering, COMSATS University Islamabad, Abbottabad Campus, Abbottabad 22060, Pakistan
2
Department of Electrical and Computer Engineering, COMSATS University Islamabad, Islamabad 45550, Pakistan
3
Center for Advanced Studies in Telecommunication (CAST), COMSATS University Islamabad, Islamabad 45550, Pakistan
4
Department of Electrical Engineering, University of Engineering and Technology Peshawar, Bannu Campus, Bannu 28100, Pakistan
5
School of Electrical and Electronics Engineering, Chung-Ang University, Dongjak-gu, Seoul 06974, Korea
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(18), 6361; https://doi.org/10.3390/app10186361
Submission received: 23 August 2020 / Revised: 8 September 2020 / Accepted: 10 September 2020 / Published: 12 September 2020
(This article belongs to the Special Issue Wind Power Systems: Design, Operation, and Control)

Abstract

:
Due to the intermittent nature of wind, there exists a major disparity between the power generation from the wind and the demand of electricity. Hence, a sophisticated maximum power point tracking (MPPT) control paradigm must be formulated for maximizing the power extraction from the wind. This research article focuses on the formulation of a nonlinear fast dynamic terminal sliding mode control (FDTSMC)-based MPPT strategy for optimizing the power extraction from a 3kW, variable speed, fixed-pitch wind energy conversion system equipped with a permanent magnet synchronous generator. The proposed MPPT strategy is compared with the benchmark fast terminal sliding mode control, conventional sliding mode control, feedback linearization control and proportional integral derivative control-based MPPT strategies under a stochastic wind speed profile. The proposed paradigm has been found superior in its tracking performance by converging the output tracking error to zero in a finite time, realizing a high precision performance, offering fast dynamic response, reducing the chattering to a minute level and guaranteeing global robustness. The superior performance and effectiveness of the proposed FDTSMC-based MPPT control paradigm is tested and validated through extensive MATLAB/Simulink simulations.

1. Introduction

Wind energy possesses an enormous potential for simultaneously addressing the worldwide rapidly growing energy demand, contributing to sustainable development and alleviating global warming concerns. Over the past few decades, this technology has undergone a very fast development. At present, the variable speed wind turbines (VSWTs) have gained a lot of attention among their other counterparts for achieving the maximum efficiency and delivering better quality of power [1,2]. Generally, due to the inconsistency of the wind speed, the wind energy conversion system (WECS) exhibits variation in its electric power output. This makes the inherent characteristics of the WECS highly nonlinear; thus, the maximum power extraction becomes a challenging task. For maximizing the efficiency of a WECS, the application of a maximum power point tracking (MPPT) control strategy is essential [3].
Different versions of alternating current (AC) generators have been employed in WECS, such as the squirrel-cage induction generator (SCIG), doubly-fed induction generator (DFIG) and permanent magnet synchronous generator (PMSG). Among the stated generators, the PMSG is preferred for a variable speed WECS due to its higher power density, higher efficiency, lower maintenance cost, lack of excitation and a wide operating speed range (0–100%) [4,5,6].
Different MPPT strategies have been formulated for maximizing the power extraction from the WECS. Most of the MPPT control strategies, formulated in the domain of the WECS, have focused on the development of nonlinear techniques. Among these strategies, one of the most promising nonlinear control strategies is the sliding mode control (SMC). The SMC possesses several attractive attributes, including simple design, fast dynamic response, insensitivity to parameter variations and disturbance rejection [7]. Apparently, the SMC strategy seems to be an ideal solution for MPPT due to its positive attributes; however, the chattering phenomenon as well as the asymptotic convergence because of a linear sliding surface in the first-order SMC invites a high level of criticism. The chattering causes the system response to oscillate around the desired reference (sliding surface), thus leading to higher mechanical wear and tear, reduced tracking performance and overheating of the power circuits [8,9].
The SMC control signal basically consists of two parts: an equivalent (or a continuous control) signal and a discontinuous (or switching control) signal. The continuous control part controls the system states once they are on the sliding surface whereas the discontinuous control part deals with the uncertainties (disturbances) but also gives rise to the chattering. The greater the uncertainties in the system, the larger the amplitude of the discontinuous control signal, and hence the chattering phenomenon is more severe. A normal strategy for the chattering minimization is the insertion of a boundary layer near the sliding surface. In this strategy, the continuous control signal is replaced by the discontinuous one once the system response is in the vicinity of the sliding surface. Through this strategy, however, the chattering is minimized, but it also gives rise to a finite steady-state error. Hence, a trade-off is required between the tracking accuracy and chattering. Minimizing the amplitude of the discontinuous control signal can be employed as an alternate strategy for the chattering reduction [10]. However, this strategy degrades the robustness as well as the transient response of the controller.
In [11], the authors have proposed a nonlinear backstepping integral sliding mode control (BISMC)-based MPPT strategy for a standalone PMSG–WECS. The system Lie derivatives have been computed using a multilayer feed-forward artificial neural network (ANN). This technique exhibited a good performance in alleviating chattering. The authors in [12] have proposed a discrete-time integral SMC strategy, consisting of a new discrete-time reaching law, to accomplish the MPPT of a PMSG–WECS. The designed controller has been implemented using an Acorn RISC Machine-based LM4F120 microcontroller, and its performance has been examined experimentally on a WECS emulator prepared in the laboratory. It has been found to exhibit a better performance than the CSMC and P&O strategies in terms of the settling time and output voltages ripple. In [13], the authors have proposed a double integral SMC-based MPPT strategy for PMSG–WECS, where the reference signal (optimum DC-side current) has been generated through a fuzzy logic controller. The proposed strategy has been found to have more promising results than the CSMC and power signal feedback control strategies in terms of eliminating the generator-side current harmonics. The authors in [14] have proposed a voltage mode second-order SMC for improving the MPPT efficiency of a PMSG–WECS. In the proposed method, the second-order sliding surface has a proportional integral derivative (PID) type structure. The proposed technique has been found to exhibit a better performance than the CSMC.
Advanced SMC strategies offering finite-time convergence, such as the conventional terminal sliding mode control (TSMC) and fast terminal sliding mode control (FTSMC), have been proposed to effectively minimize chattering. Both of these stated strategies utilize nonlinear sliding surfaces, where fractional power is introduced into the sliding surface to guarantee a fast as well as finite-time convergence during the sliding phase. The convergence speed of FTSMC is higher than the conventional TSMC [15]. However, just like the conventional SMC, the chattering phenomenon by TSMC and FTSMC is not eliminated completely. In [10], a novel fast dynamic terminal sliding mode control (FDTSMC) strategy, based on Lyapunov theory is proposed for a class of single-input and single-output uncertain nonlinear systems, for complete elimination of chattering. However, to the best of the authors understanding and knowledge, the FDTSMC-based paradigm has never been applied to accomplish the MPPT of a variable speed PMSG–WECS.
In this article, a nonlinear MPPT control paradigm based on FDTSMC is proposed for a 3 kW, variable speed, fixed-pitch, PMSG–WECS-based standalone system. The proposed MPPT strategy converges the output tracking error to zero in a finite-time, realizes a high precision performance, offers fast dynamic response, reduces chattering to a minute level and guarantees global robustness. The performance of the proposed paradigm is tested and validated under a stochastic wind speed profile in Matlab/Simulink. A comparative analysis is performed among the proposed MPPT strategy, FTSMC, conventional sliding mode control (CSMC), feedback linearization control (FBLC) and proportional integral derivative (PID) control-based MPPT techniques, where the proposed MPPT strategy has been found superior to all the stated techniques in terms of the tracking performance.

2. Mathematical Modeling of PMSG–WECS

This section covers the mathematical modeling of the standalone, fixed-pitch, variable speed PMSG–WECS. A schematic of the overall variable speed PMSG–WECS is demonstrated in Figure 1. The principal components of a typical variable speed PMSG–WECS include: a variable speed wind turbine (VSWT), a gearbox, power electronic converters and a PMSG coupled with a VSWT.

2.1. Wind Turbine Mathematical Modeling

The mechanical power, P m , captured by the VSWT is expressed as follows [16]:
P m = 1 2 ρ π R t 2 v w 3 C p ( λ , β )
where ρ denotes the air density (generally, its value is 1.25 kg /m 3 , at 10 C , at sea level), R t stands for the turbine blade radius (m), v w represents the wind speed (m/s), and C p indicates the wind turbine rotor power conversion coefficient that defines the wind turbine rotor aerodynamic efficiency. Generally, C p varies with both the tip speed ratio, λ , and the blade pitch angle, β . However, β is assumed to be a constant ( β = 0 ) in this article.
C p is determined by numerical approximation; such as in [17], for β = 0 , it is given as follows:
C p ( λ ) = 0.0061 λ 0.0013 λ 2 + 0.0081 λ 3 9.7477 × 10 4 λ 4 6.5416 × 10 5 λ 5 + 1.3027 × 10 6 λ 6 4.54 × 10 7 λ 7
C p has a unique maximum value (say, C p m a x ) at a particular λ , called λ o p t . Thus, a VSWT tracks the C p m a x , while maintaining the tip speed ratio at its optimal value, λ o p t , for the maximum wind power extraction.
The ratio of the peripheral speed of the wind turbine blades to the wind speed indicates the tip speed ratio, as follows:
λ = Ω R t v w
where Ω represents the rotational speed of the wind turbine blades (rad/s, i.e., rotational speed of the low-speed shaft).
The wind turbine mechanical torque (aerodynamic torque), T m , can be obtained from Equations (1) and (3), as follows:
T m = P m Ω = 1 2 ρ π R t 3 v w 2 C t ( λ )
where C t ( λ ) = C p ( λ ) / λ represents the wind turbine torque coefficient. It can be found from a second-order polynomial expression, as a function of the tip speed ratio, as follows:
C t ( λ ) = α 0 + α 1 λ + α 2 λ 2
Now, substituting Equations (3) and (5) into Equation (4), T m can also be expressed as follows:
T m = d 1 v w 2 + d 2 v w Ω h g r + d 3 Ω h 2 g r 2
where Ω h is rotational speed of the high speed shaft (PMSG rotor) and g r = Ω h / Ω is the gear ratio (or transmission ratio). Moreover,
d 1 = 0.5 π ρ R t 3 α 0 , d 2 = 0.5 π ρ R t 4 α 1 , d 3 = 0.5 π ρ R t 5 α 2 α 0 = 0.1253 , α 1 = 0.0047 , α 2 = 0.0005
Knowing C p m a x and λ o p t , the optimal mechanical power, P m o p t , generated by the VSWT can be worked out from Equation (1) by substituting Equation (3) into it as follows:
P m o p t = ρ π R t 5 Ω 3 C p m a x ( λ ) 2 λ o p t 3
Figure 2 demonstrates the mechanical power output of the VSWT under different wind speeds. There exists a unique maximum power point, P m o p t , at each wind speed. All these maximum power points constitute an optimal regime characteristic (ORC). The ORC indicates the operating region where the maximum possible energy can be extracted from the fluctuating wind. The primary goal of the MPPT control strategy is to continuously operate the VSWT on the ORC, despite fluctuations in the wind speed.
The optimal mechanical torque can be obtained from the optimal mechanical power given in Equation (7) as follows:
T m o p t = P m o p t Ω = ρ π R t 5 Ω 2 C p m a x ( λ ) 2 λ o p t 3
Table 1 and Table 2 describe different significant parameters of the overall PMSG–WECS along with its control system.

2.2. Permanent Magnet Synchronous Generator Mathematical Modeling

The PMSG can be conveniently modeled in the d-q (or synchronously rotating) reference frame, where the zero components are discarded. The axis along the axis of the rotor is called the direct or the d-axis, while the axis perpendicular to the d-axis is known as the quadrature or q-axis.
For a standalone PMSG–WECS, the voltages on the d-q axes are considered as the output variables. Based on this concept, the PMSG model can be described in terms of the d-q axes voltages by the following group of equations [18,19]:
L d + L c h d i d d t = R s + R c h i d + p L q L c h i q Ω h L q + L c h d i q d t = R s + R c h i q p L d + L c h i q Ω h + p Φ m Ω h J h d Ω h d t = T m g r T e m = T m g r p L d L q i d i q + p Φ m i q = d 1 v w 2 g r + d 2 v w Ω h g r 2 + d 3 Ω h 2 g r 3 p Φ m i q }
where R s denotes the stator resistance, L d and L q represent the stator d- and q- a x i s inductances, respectively, i d and i q stand for the d- and q- a x i s stator currents, respectively, p indicates the number of pole pairs, Φ m indicates the maximum value of magnetic flux and T e m represents the electromagnetic torque of the PMSG, resulting from the interaction between the stator and rotor fluxes. Moreover, J h denotes the moment of inertia of the high-speed shaft (PMSG) and R c h , L c h indicate the resistance and inductance of the load, respectively. In Equation (9), the non-salient pole type of PMSG has been considered, for which L d = L q .
The WECS is inherently a highly nonlinear system. The objective of this work is to maximize the power extraction from the PMSG–WECS, which is accomplished by controlling the shaft speed. For this purpose, the nonlinear PMSG–WECS model can be written in state-space representation as follows:
x ˙ = f ( x ) + g ( x ) u y = h ( x )
where
x ˙ 1 x ˙ 2 x ˙ 3 x ˙ = a 2 x 1 + a 3 x 2 x 3 b 2 x 2 + b 3 x 1 x 3 + b 4 x 3 1 J h d 1 v w 2 g r + d 2 v w x 3 g r 2 + d 3 x 3 2 g r 3 p Φ m x 2 f ( x ) + x 1 a 1 x 2 b 1 0 g ( x ) . R c h u y = Ω h h ( x ) }
where f ( x ) = [ f 1 f 2 f 3 ] T , g ( x ) = [ g 1 g 2 g 3 ] T , x = [ x 1 x 2 x 3 ] T = [ i d i q Ω h ] T R n represents the state vector, f ( x ) and g ( x ) stand for the nonlinear smooth vector fields, and u and y indicate the control input and output, respectively. Moreover,
a 1 = L d + L c h ; a 2 = R s a 1 ; a 3 = p ( L q L c h ) a 1 b 1 = L q + L c h ; b 2 = R s b 1 ; b 3 = p ( L d + L c h ) b 1 ; b 4 = p Φ m b 1
Table 1 and Table 2 describe different significant parameters of the overall PMSG–WECS along with its control system.

3. Coordinates Transformation

To express the system in normal form, the coordinates transformation must be carried out, as follows:
z 1 = h ( x ) = x 3 = Ω h z 2 = L f h ( x ) = h ( x ) x f ( x ) = ξ 1 v w 2 + ξ 2 v w x 3 + ξ 3 x 3 2 ξ 4 x 2 z 3 = a 3 x 1 x 2 }
where ξ 1 = d 1 J h g r , ξ 2 = d 2 J h g r 2 , ξ 3 = d 3 J h g r 3 and ξ 4 = d 4 = p Φ m .
Since, the system under consideration possesses a relative degree, r < n (system order), where r = 2 and n = 3 , it implies that the system is only partially linearizable [20,21]. Moreover, the input-output form (normal form) can be expressed as follows:
z ˙ 1 = z 2 z ˙ 2 = L f 2 h ( x ) + L g L f h ( x ) u z ˙ 3 = L f 3 h ( x ) + L g L f 2 h ( x ) u ˙ }
where one of the transformed system states, i.e., z 3 represents the internal (or zero) dynamic state. Now, the Lie derivatives are given as follows:
L f 2 h ( x ) = d 4 ( b 2 x 1 + b 3 x 1 x 3 + b 4 x 3 ) + ( d 2 v w + 2 d 3 x 3 ) ( ξ 1 v w 2 + ξ 2 v w x 3 + ξ 3 x 3 2 ξ 4 x 2 ) L g L f h ( x ) = d 4 x 2 a 1 0 }

4. Lie Derivatives Estimation via Adaptive Neuro-Fuzzy Inference System

In practical contexts, one may not be available with the nonlinear terms and aerodynamic forces. Hence, in this work, an offline Takagi–Sugeno–Kang (TSK)-based ANFIS scheme is designed for the estimation of the nonlinear drift of the PMSG–WECS (i.e., L g L f h ( x ) , L g L f 2 h ( x ) , L f 2 h ( x ) and L f 3 h ( x ) ). The TSK-based ANFIS is a MISO (multiple-input-single-output) system [22]. Four ANFIS networks are trained for estimating the system Lie derivatives using the inputs ( z 1 , z 2 , z 3 and v w ) and outputs ( L g L f h ( x ) , L g L f 2 h ( x ) , L f 2 h ( x ) and L f 3 h ( x ) ). The input-output data are generated by simulating the system for a particular wind speed profile. The ANFIS network then mimics the relationships between the inputs and the corresponding outputs. The ANFIS estimation network comprises six layers. Layer 1 consists of all the inputs that are available during the estimation process. Layer 2 generates the membership functions based on triangular functions. Layer 3 represents the antecedent connective part of the extracted IF-THEN fuzzy rule. Layers 4, 5 and 6 represent the normalization layer, the consequent or regressor layer and the summation layer, respectively.
The complete ANFIS network structure for Lie derivatives estimation is demonstrated in Figure 3.

5. MPPT Control Design for PMSG–WECS

The main goal of the MPPT control design is to maximize the power output of the PMSG–WECS. This goal is accomplished by proposing an FDTSMC-based MPPT control strategy.
The conventional TSMC strategy may not be able to guarantee the same convergence speed for the system states located distant from the equilibrium point, O. However, it can offer a finite-time convergence for the system states located near the equilibrium point.
The conventional TSMC concept was first introduced in [23], for which the nonlinear sliding surface can be expressed as follows:
s 1 = e ˙ + β 1 e q 1 p 1 = 0
where e represents the tracking error and β 1 , p 1 , q 1 > 0 . Note that both p 1 and q 1 must be odd integers, such that p 1 > q 1 . It can easily be verified that for any given initial tracking error, e ( 0 ) 0 , the dynamics expressed in Equation (14) will converge to e = 0 in a finite time, as determined by the reaching time, t s 1 = [ p 1 / { β 1 p 1 q 1 } ] | e ( 0 ) | p 1 q 1 / p 1 , obtained from the analytical solution of Equation (14). The equilibrium point O is termed as the terminal attractor. The reaching time can be tuned by proper selection of the parameters β 1 , p 1 and q 1 .
The inclusion of the nonlinear term e q 1 / p 1 in Equation (14) improves the convergence towards the equilibrium point such that the closer the e q 1 / p 1 is to the equilibrium point, the faster the convergence will be, thus offering a finite-time convergence. However, for the system states located distant from the equilibrium point, the conventional TSMC strategy does not outperform its conventional SMC counterpart (where p 1 = q 1 ) because the e q 1 / p 1 term now tends to retard the convergence towards the equilibrium point. To counteract this problem, a solution has been proposed in [24] that is known as the FTSMC, for which the nonlinear sliding surface can be expressed as follows:
s 2 = e ˙ + c 1 e + β 2 e q 2 / p 2 = 0
where c 1 , β 2 , p 2 , q 2 > 0 . Note that both p 2 and q 2 must be odd integers, such that p 2 > q 2 . Thus, Equation (15) gives e ˙ = c 1 e β 2 e q 2 / p 2 . By properly choosing p 2 and q 2 , for a given initial tracking error e ( 0 ) 0 , the dynamics expressed in Equation (15) will converge to e = 0 in a finite time. Physically, it can be interpreted as: when e is located far away from the equilibrium point, the approximate dynamics can then be expressed as: e ˙ = c 1 e , which offers a fast and finite-time convergence. However, when e is located close to the equilibrium point, the approximate dynamics then become: e = β 2 e q 2 / p 2 , which is the terminal attractor offering fast and finite-time convergence. More precisely, the analytical solution of Equation (15) gives the exact time to reach the equilibrium point (reaching time) as follows:
t s 2 = p 2 c 1 p 2 q 2 ln { c 1 e ( 0 ) p 2 q 2 / p 2 + β 2 } ln β 2

FDTSMC-Based MPPT Control Design for PMSG–WECS

The FTSMC renders the advantages of fast dynamic response, high steady-state tracking accuracy and finite-time convergence. However, chattering of the control signal due to the switching action of the discontinuous control signal, and sensitivity to high-frequency noise, still remain as the main disadvantages of the FTSMC. To counteract these issues, a novel strategy has been proposed in [25], namely the global FDTSMC. In the stated strategy, the time derivative of the control input has been taken as the control variable for the augmented system such that the control input appears as chattering-free due to an integrator (acting as a low-pass filter) placed next to the system.
In this article, the FDTSMC offers a finite-time rotor speed tracking error convergence of the PMSG–WECS for ensuring the MPPT. This control strategy is developed using a function augmented nonlinear sliding surface that guarantees the tracking error convergence to the equilibrium point in a finite time. Moreover, the high-frequency chattering is completely eliminated by eliminating the reaching phase.
The input-output form for FDTSMC can be described as follows:
z ˙ 1 = z 2 z ˙ 2 = L ^ f 2 h ( x ) + L ^ g L f h ( x ) u z ˙ 3 = L ^ f 3 h ( x ) + L ^ g L f 2 h ( x ) u ˙ }
where L ^ f 3 h ( x ) and L ^ g L f 2 h ( x ) have been expressed in Equations (A1) and (A2), in Appendix A, respectively.
The nonlinear sliding surface for the proposed FDTSMC can be expressed as follows:
s 3 = σ ˙ + c 2 σ + β 3 σ q 3 / p 3
where c 2 , β 3 , p 3 , q 3 > 0 . Note that both p 3 and q 3 must be odd integers, such that ( p 3 > q 3 ). The term σ in Equation (18) denotes a sliding mode surface that can be expressed as follows:
σ = e ˙ + c 3 e = ( z 2 z ˙ r e f ) + c 3 ( z 1 z r e f )
Now, the successive differentiation of Equation (19) yields:
σ ˙ = ( z ˙ 2 z ¨ r e f ) + c 3 ( z ˙ 1 z ˙ r e f ) = ( z 3 z ¨ r e f ) + c 3 ( z 2 z ˙ r e f ) σ ¨ = ( z ˙ 3 z r e f ) + c 3 ( z ˙ 2 z ¨ r e f ) = L ^ f 3 h ( x ) + L ^ g L f 2 h ( x ) u ˙ z r e f + c 3 ( e ¨ ) }
Next, the FDTSMC sliding surface expressed in Equation (18) is differentiated and simplified using Equation (20) in order to derive the equivalent controller as follows:
s 3 ˙ = L ^ f 3 h ( x ) + L ^ g L f 2 h ( x ) u ˙ z r e f + c 3 e ¨ + c 2 σ ˙ + β 3 q 3 p 3 σ q 3 p 3 1 σ ˙
The equivalent controller, u ˙ e q D , can be obtained by setting Equation (21) to zero as follows:
u ˙ e q D = 1 L ^ g L f 2 h ( x ) L ^ f 3 h ( x ) z r e f + c 3 e ¨ + c 2 σ ˙ + β 3 q 3 p 3 σ q 3 p 3 1 σ ˙
While, the discontinuous controller, u ˙ d i s c D , can be given as follows:
u ˙ d i s c D = K 1 s 3 K 2 sign ( s 3 )
where K 1 , K 2 > 0 are the switching controller gains.
Finally, the FDTSMC-based MPPT law is obtained by combining Equations (22) and (23), and passing through an integrator (low-pass filter) as follows:
u F D T S M C = [ 1 L ^ g L f 2 h ( x ) L ^ f 3 h ( x ) z r e f + c 3 e ¨ + c 2 σ ˙ + β 3 q 3 p 3 σ q 3 p 3 1 σ ˙ u ˙ e q D K 1 s 3 K 2 sign ( s 3 ) u ˙ d i s c D ] . d t
The computational flow chart for implementation of the proposed MPPT strategy is depicted in Figure 4, while the closed-loop control system for PMSG–WECS is given in Figure 5.

6. Stability Analysis

To verify the PMSG–WECS closed-loop system stability under FDTSMC, Lyapunov theory is employed.

Stability Analysis for FDTSMC

First of all, a Lyapunov candidate function, V f 1 , and its derivative, V ˙ f 1 , are defined along the system trajectory as follows:
V f 1 = 1 2 s 3 2 V ˙ f 1 = s 3 s ˙ 3 }
Now, s ˙ 3 from Equation (21) is substituted into Equation (25) as follows:
V ˙ f 1 = s 3 L ^ f 3 h ( x ) + L ^ g L f 2 h ( x ) u ˙ F D T S M C z r e f + c 3 e ¨ + c 2 σ ˙ + β 3 q 3 p 3 σ q 3 p 3 1 σ ˙
Now, u F D T S M C from Equation (24) is substituted into Equation (26), giving:
V ˙ f 1 = K 1 s 3 2 K 2 s 3 sign ( s 3 ) K 1 s 3 2 K 2 | s 3 | K 1 , K 2 > 0 , V ˙ f 1 < 0

7. Simulation Results and Discussion

In this section, the proposed nonlinear MPPT control scheme, namely FDTSMC, is tested and evaluated in MATLAB/Simulink through simulations. The FTSMC [26], CSMC, FBLC [27] and PID have been used as benchmarks to manifest the superior tracking performance of the proposed MPPT paradigm. For performance evaluation, the wind speed is assumed to have a stochastic profile.
Figure 4 and Figure 5 demonstrate the application of the proposed nonlinear MPPT control paradigm to the standalone, variable speed, fixed-pitch, 3 k W PMSG–WECS. The MPPT strategy ensures the maximum power extraction from the VSWT while keeping the tip speed ratio, λ , at its optimal value, λ o p t . Table 1 describes different significant parameters of the proposed FDTSMC.

7.1. Performance Evaluation under Stochastic Wind Speed Profile

The performance of the proposed MPPT paradigm has been tested and evaluated, in this section, for a 3 k W PMSG–WECS under a stochastic wind speed profile that comprises rapid and frequent wind speed fluctuations. The maximum power conversion coefficient, C p m a x = 0.4762 , computed from Equation (2) occurs at an optimal tip speed ratio, λ o p t = 7 . The simulations are performed for a 100 s horizon at an average wind speed of, v w a v = 7 m / s (also depicted in Figure 2).
Figure 6 demonstrates the PMSG rotational speed tracking performance. It is evident from the zoomed-in segments of the figure that the FDTSMC is offering a superior tracking performance among its competitors. Moreover, both the FTSMC and FDTSMC offer a finite-time convergence. The superior performance of the FDTSMC to the FTSMC, CSMC, FBLC and PID, as well as the finite-time convergence of both the FTSMC and FDTSMC, are also evident from the speed tracking error of the PMSG, as depicted in Figure 7. The CSMC exhibits slightly better performance than both the FBLC and PID.
Figure 8 indicates the tip speed ratio performance of the VSWT. The PID has the worst performance, and it offers a lot of oscillations around the λ o p t . However, both the FTSMC and FDTSMC are accurately tracking the λ o p t with FDTSMC being the best MPPT candidate.
Figure 9 shows the power conversion coefficient of the VSWT. It is clear that both the FTSMC and FDTSMC strategies are tracking the optimal power conversion coefficient, C p m a x = 0.4762 , throughout the selected stochastic wind speed profile, thus guaranteeing the maximum wind power extraction. The zoomed-in segments of the figure clearly depict that, based on the C p m a x tracking, the FDTSMC is the best MPPT candidate, while the PID is the worst.
Figure 10 and Figure 11 indicate the low-speed shaft mechanical power and the high-speed shaft mechanical power versus the tip speed ratio, respectively. It is evident that both the FTSMC and the FDTSMC strategies maintain the low-speed as well as the high-speed shaft mechanical powers in the close neighborhood of the optimal tip speed ratio, λ o p t . On the other hand, both the FBLC and PID techniques have very poor performance on account of having a lot of power fluctuations around the optimal tip speed ratio.
Figure 12 represents the low-speed shaft mechanical power versus the low-speed shaft rotational speed. The proposed MPPT controller efficiently maintains the rotational speed of the VSWT in such a manner that the system operates on the ORC (optimal regimes characteristic)—that is, the region where the maximum wind power can be extracted. However, both the FBLC and PID schemes experience a lot of speed as well as mechanical power fluctuations.
Figure 13 illustrates the PMSG electromagnetic torque versus the tip speed ratio. Since both the FTSMC- and the FDTSMC-based MPPT candidates maintain the torque at the optimal tip speed ratio, λ o p t = 7 , it implies that both are capable of extracting the maximum wind power. The FDTSMC has the best, while the PID has the worst performance because of frequent torque fluctuations. Similarly, Figure 14 shows the VSWT mechanical torque versus the low-speed shaft rotational speed. Again, both the FTSMC- and the FDTSMC-based MPPT candidates maintain the torque around the ORC. It signifies that these two MPPT strategies are capable of extracting the maximum wind power, with FDTSMC being superior to the remaining MPPT candidates.

7.2. Performance Indices

To further validate the superior performance of the proposed MPPT strategy to the FTSMC-, CSMC-, FBLC- and PID-based MPPT strategies, the dynamic performance of each MPPT candidate is also evaluated through four different performance indices, given as follows [11,28,29]:
I S E = 0 t s i m e t 2 d t I T S E = 0 t s i m t e t 2 d t I A E = 0 t s i m | e t | d t I T A E = 0 t s i m t | e t | d t }
where t s i m indicates the total simulation time.
The performances indices, ISE (integral squared error), ITSE (integral of time squared error), IAE (integral absolute error) and ITAE (integral of time absolute error) have been computed for PMSG–WECS using e ( t ) = Ω h Ω r e f for each MPPT strategy. It is evident from Figure 15, Figure 16, Figure 17 and Figure 18 that as the time advances, the accumulative error of each MPPT controller also increases. However, the proposed FDTSMC renders the flattest profile and the smallest error in each case, thus guaranteeing its superior performance. Moreover, the FTSMC is the second best MPPT candidate, while the PID is the worst.

8. Conclusions

To extract the maximum possible power from the wind and to reduce chattering, in this research article, a nonlinear MPPT control strategy based on FDTSMC has been proposed for a variable speed, standalone, fixed-pitch, 3 k W PMSG–WECS. The proposed MPPT strategy has been simulated under a stochastic wind speed profile in MATLAB/Simulink. Its MPPT performance has been compared with and found superior to the FTSMC-, CSMC-, FBLC- and PID-based benchmark MPPT techniques in terms of offering a finite-time convergence, accurate MPPT, low steady-state error, fast dynamic response and minute chattering. The stability of the closed-loop system is demonstrated by the Lyapunov stability theory. The results are promising, and it has been observed that the FDTSMC has the best MPPT performance among the candidates, while the conventional PID has the worst.

Author Contributions

M.Z. developed the concept, conducted simulations and wrote the initial manuscript draft; L.K. and J.-S.R. supervised this research work; Q.K., S.U. and I.S. analyzed the data and reviewed/edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Basic Science Research Program through the National Research Foundation of Korea funded by the Ministry of Education (2016R1D1A1B01008058) and by Human Resources Development (No.20204030200090) of Korean Institute of Energy Technology Evaluation and Planning (KETEP) grant funded by the Korea government Ministry of Trade, Industry, and Energy.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

L ^ f 3 h ( x ) = d 4 b 2 x 2 + b 2 b 3 d 4 x 1 x 3 + b 2 b 4 d 4 x 3 2 b 2 d 4 x 2 u b 1 + a 2 b 3 d 4 x 1 x 3 + a 3 b 3 d 4 x 2 x 3 2 b 3 d 4 x 1 x 3 u a 1 + b 3 ξ 1 d 4 v w 2 + b 3 ξ 2 d 4 x 1 x 3 + b 3 ξ 3 d 4 x 1 x 3 2 + b 3 d 4 2 x 1 x 2 J h + b 3 d 4 x 1 x 3 u b 1 b 4 d 4 x 3 u b 1 + d 4 x 2 u 2 b 1 2 + d 2 v w g r 2 + 2 d 3 x 3 g r 3 [ ξ 2 v w + 2 ξ 3 x 3 ξ 1 v w 2 + ξ 2 v w x 3 + ξ 3 x 3 2 + d 4 x 2 + d 4 b 2 x 2 J h + b 3 d 4 x 1 x 3 d 4 x 2 u b 1 + d 4 x 3 ]
L ^ g L f 2 h ( x ) = d 4 x 2 J h b 1

References

  1. Soliman, M.A.; Hasanien, H.M.; Azazi, H.Z.; El-Kholy, E.E.; Mahmoud, S.A. Linear-quadratic regulator algorithm-based cascaded control scheme for performance enhancement of a variable-speed wind energy conversion system. Arab. J. Sci. Eng. 2019, 44, 2281–2293. [Google Scholar] [CrossRef]
  2. Worku, M.Y.; Abido, M.A. Fault Ride-Through and Power Smoothing Control of PMSG-Based Wind Generation Using Supercapacitor Energy Storage System. Arab. J. Sci. Eng. 2019, 44, 2067–2078. [Google Scholar] [CrossRef]
  3. Tahir, K.; Belfedal, C.; Allaoui, T.; Denai, M.; Doumi, M.H. A new sliding mode control strategy for variable-speed wind turbine power maximization. Int. Trans. Electr. Energy Syst. 2018, 28, e2513. [Google Scholar] [CrossRef] [Green Version]
  4. Housseini, B.; Okou, A.F.; Beguenane, R. Performance comparison of variable speed PMSG-based wind energy conversion system control algorithms. In Proceedings of the 2017 Twelfth International Conference on Ecological Vehicles and Renewable Energies (EVER), Monte-Carlo, Monaco, 11–13 April 2017; pp. 1–10. [Google Scholar]
  5. Yaramasu, V.; Dekka, A.; Durán, M.J.; Kouro, S.; Wu, B. PMSG-based wind energy conversion systems: Survey on power converters and controls. IET Electr. Power Appl. 2017, 11, 956–968. [Google Scholar] [CrossRef]
  6. Chinchilla, M.; Arnaltes, S.; Burgos, J.C. Control of permanent-magnet generators applied to variable-speed wind-energy systems connected to the grid. IET Electr. Power Appl. 2006, 21, 130–135. [Google Scholar] [CrossRef] [Green Version]
  7. Munteanu, I.; Bacha, S.; Bratcu, A.I.; Guiraud, J.; Roye, D. Energy-reliability optimization of wind energy conversion systems by sliding mode control. IEEE Trans. Energy Convers. 2008, 23, 975–985. [Google Scholar] [CrossRef]
  8. Shtessel, Y.; Edwards, C.; Fridman, L.; Levant, A. Sliding Mode Control and Observation; Springer: New York, NY, USA, 2014. [Google Scholar]
  9. Young, K.D.; Utkin, V.I.; Ozguner, U. A control engineer’s guide to sliding mode control. IEEE Trans. Control Syst. Technol. 1999, 7, 328–342. [Google Scholar] [CrossRef] [Green Version]
  10. Liu, J.; Sun, F. A novel dynamic terminal sliding mode control of uncertain nonlinear systems. J. Control Theoryand Appl. 2007, 5, 189–193. [Google Scholar] [CrossRef]
  11. Khan, I.U.; Khan, L.; Khan, Q.; Ullah, S.; Khan, U.; Ahmad, S. Neuro-Adaptive Backstepping Integral Sliding Mode Control Design for Nonlinear Wind Energy Conversion System. Turk. J. Electr. Eng. Comput. Sci. 2020, in press. [Google Scholar] [CrossRef]
  12. Yazici, İ.; Yaylaci, E.K. Maximum power point tracking for the permanent magnet synchronous generator-based WECS by using the discrete-time integral sliding mode controller with a chattering-free reaching law. IET Power Electron. 2017, 10, 1751–1758. [Google Scholar] [CrossRef]
  13. Yin, X.X.; Lin, Y.G.; Li, W.; Liu, H.W.; Gu, Y.J. Fuzzy-logic sliding-mode control strategy for extracting maximum wind power. IEEE Trans. Energy Convers. 2015, 30, 1267–1278. [Google Scholar] [CrossRef]
  14. Dursun, E.H.; Kulaksiz, A.A. Second-order sliding mode voltage-regulator for improving MPPT efficiency of PMSG-based WECS. Int. J. Electr. Power Energy Syst. 2020, 121, 106149. [Google Scholar] [CrossRef]
  15. Boukattaya, M.; Mezghani, N.; Damak, T. Adaptive nonsingular fast terminal sliding-mode control for the tracking problem of uncertain dynamical systems. ISA Trans. 2018, 77, 1–19. [Google Scholar] [CrossRef]
  16. El Mourabit, Y.; Derouich, A.; El Ghzizal, A.; El Ouanjli, N.; Zamzoum, O. Nonlinear backstepping control for PMSG wind turbine used on the real wind profile of the Dakhla-Morocco city. Int. Trans. Electr. Energy Syst. 2020, 30, e12297. [Google Scholar] [CrossRef]
  17. Munteanu, I.; Bratcu, A.I.; Cutululis, N.A.; Ceangǎ, E. Optimal Control of Wind Energy Systems: Towards a Global Approach; Springer: London, UK, 2008. [Google Scholar]
  18. Cutululis, N.A.; Ceangǎ, E.; Hansen, A.D.; Sørensen, P. Robust multi-model control of an autonomous wind power system. Wind Energy 2006, 9, 399–419. [Google Scholar] [CrossRef]
  19. Cheikh, R.; Menacer, A.; Alaoui, L.C.; Drid, S. Robust nonlinear control via feedback linearization and Lyapunov theory for permanent magnet synchronous generator-based wind energy conversion system. Front. Energy 2020, 14, 180–191. [Google Scholar] [CrossRef]
  20. Slotine, J.J.E.; Li, W. Applied Nonlinear Control; Prentice Hall: Englewood Cliffs, NJ, USA, 1991. [Google Scholar]
  21. Marquez, H.J. Nonlinear Control Systems; John Wiley & Sons: Hoboken, NJ, USA, 2003. [Google Scholar]
  22. Lau, C.K.; Ghosh, K.; Hussain, M.A.; Hassan, C.C. Fault diagnosis of Tennessee Eastman process with multi-scale PCA and ANFIS. Chemom. Intell. Lab. Syst. 2013, 120, 1–14. [Google Scholar] [CrossRef]
  23. Zhihong, M.; Paplinski, A.P.; Wu, H.R. A robust MIMO terminal sliding mode control scheme for rigid robotic manipulators. IEEE Trans. Autom. Control 1994, 39, 2464–2469. [Google Scholar] [CrossRef]
  24. Yu, X.; Zhihong, M. Fast terminal sliding-mode control design for nonlinear dynamical systems. IEEE Trans. Circuits Syst. I Fundam. Theory Appl. 2002, 49, 261–264. [Google Scholar]
  25. Ziong, J.J.; Zhang, G.B. Global fast dynamic terminal sliding mode control for a quadrotor UAV. ISA Trans. 2017, 66, 233–240. [Google Scholar]
  26. Dursun, E.H.; Kulaksiz, A.A. Second-Order Fast Terminal Sliding Mode Control for MPPT of PMSG-based Wind Energy Conversion System. Elektron. Elektrotech. 2020, 26, 39–45. [Google Scholar] [CrossRef]
  27. Soufi, Y.; Kahla, S.; Bechouatb, M. Feedback linearization control based particle swarm optimization for maximum power point tracking of wind turbine equipped by PMSG connected to the grid. Int. J. Hydrogen Energy 2016, 41, 20950–20955. [Google Scholar] [CrossRef]
  28. Sahib, M.A.; Ahmed, B.S. A new multiobjective performance criterion used in PID tuning optimization algorithms. J. Adv. Res. 2016, 7, 125–134. [Google Scholar] [CrossRef] [Green Version]
  29. Ali, K.; Khan, Q.; Ullah, S.; Khan, I.; Khan, L. Nonlinear robust integral backstepping based MPPT control for stand-alone photovoltaic system. PLoS ONE 2020, 15, e0231749. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Schematic of the overall variable speed permanent magnet synchronous generator–wind energy conversion system (PMSG–WECS).
Figure 1. Schematic of the overall variable speed permanent magnet synchronous generator–wind energy conversion system (PMSG–WECS).
Applsci 10 06361 g001
Figure 2. Wind turbine mechanical power output under different wind speeds.
Figure 2. Wind turbine mechanical power output under different wind speeds.
Applsci 10 06361 g002
Figure 3. ANFIS network structure for PMSG–WECS Lie-derivatives estimation.
Figure 3. ANFIS network structure for PMSG–WECS Lie-derivatives estimation.
Applsci 10 06361 g003
Figure 4. Computational flow chart for implementation of the proposed MPPT strategy.
Figure 4. Computational flow chart for implementation of the proposed MPPT strategy.
Applsci 10 06361 g004
Figure 5. Closed-loop control system for the PMSG–WECS.
Figure 5. Closed-loop control system for the PMSG–WECS.
Applsci 10 06361 g005
Figure 6. High-speed shaft rotational speed tracking.
Figure 6. High-speed shaft rotational speed tracking.
Applsci 10 06361 g006
Figure 7. High-speed shaft rotational speed tracking error.
Figure 7. High-speed shaft rotational speed tracking error.
Applsci 10 06361 g007
Figure 8. Tip speed ratio performance.
Figure 8. Tip speed ratio performance.
Applsci 10 06361 g008
Figure 9. Power conversion coefficient.
Figure 9. Power conversion coefficient.
Applsci 10 06361 g009
Figure 10. Low-speed shaft mechanical power versus the tip speed ratio.
Figure 10. Low-speed shaft mechanical power versus the tip speed ratio.
Applsci 10 06361 g010
Figure 11. High-speed shaft mechanical power versus the tip speed ratio.
Figure 11. High-speed shaft mechanical power versus the tip speed ratio.
Applsci 10 06361 g011
Figure 12. Mechanical power versus the low-speed shaft rotational speed.
Figure 12. Mechanical power versus the low-speed shaft rotational speed.
Applsci 10 06361 g012
Figure 13. Electromagnetic torque versus the tip speed ratio performance.
Figure 13. Electromagnetic torque versus the tip speed ratio performance.
Applsci 10 06361 g013
Figure 14. Mechanical torque versus the low-speed shaft rotational speed.
Figure 14. Mechanical torque versus the low-speed shaft rotational speed.
Applsci 10 06361 g014
Figure 15. Integral squared error.
Figure 15. Integral squared error.
Applsci 10 06361 g015
Figure 16. Integral of time squared error.
Figure 16. Integral of time squared error.
Applsci 10 06361 g016
Figure 17. Integral absolute error.
Figure 17. Integral absolute error.
Applsci 10 06361 g017
Figure 18. Integral of time absolute error.
Figure 18. Integral of time absolute error.
Applsci 10 06361 g018
Table 1. PMSG–WECS and the proposed maximum power point tracking (MPPT) controller parameters.
Table 1. PMSG–WECS and the proposed maximum power point tracking (MPPT) controller parameters.
NameQuantityValue
Wind TurbineAir density (at 10 C , at sea level), ρ 1.25 kg /m 3
Turbine blade radius, R t 2.50 m
Transmission (or gear) ratio, g r 7
Average wind speed, v w a v 7 m / s
Optimal tip speed ratio, λ o p t 7
Maximum power conversion coefficient, C p m a x 0.4762
PMSGStator resistance, R s 3.30 Ω
Stator d-axis inductance, L d 41.56 mH
Stator q-axis inductance, L q 41.56 mH
Magnetic flux, Φ m 438.20 mWb
Pole pairs, p3
High-speed shaft moment of inertia, J h 0.05522 kg m 2
Load inductance, L c h 0.008 H
FDTSMCGain, p 3 3
Gain, q 3 1
Gain, c 2 100
Gain, c 3 7000
Gain, β 3 200
Gain, K 1 4000
Gain, K 2 0.1
Table 2. Specification of various important parameters.
Table 2. Specification of various important parameters.
ParameterValueParameterValueParameterValue
a 1 0.1215 b 3 −6 d 4 2.6292
a 2 −27.1471 b 4 21.6288 ξ 1 9.9399
a 3 −0.9486 d 1 3.8422 ξ 2 −0.13316
b 1 0.1215 d 2 −0.3603 ξ 3 5.059 × 10 3
b 2 −27.1471 d 3 −0.0958 ξ 4 2.6292

Share and Cite

MDPI and ACS Style

Zafran, M.; Khan, L.; Khan, Q.; Ullah, S.; Sami, I.; Ro, J.-S. Finite-Time Fast Dynamic Terminal Sliding Mode Maximum Power Point Tracking Control Paradigm for Permanent Magnet Synchronous Generator-Based Wind Energy Conversion System. Appl. Sci. 2020, 10, 6361. https://doi.org/10.3390/app10186361

AMA Style

Zafran M, Khan L, Khan Q, Ullah S, Sami I, Ro J-S. Finite-Time Fast Dynamic Terminal Sliding Mode Maximum Power Point Tracking Control Paradigm for Permanent Magnet Synchronous Generator-Based Wind Energy Conversion System. Applied Sciences. 2020; 10(18):6361. https://doi.org/10.3390/app10186361

Chicago/Turabian Style

Zafran, Muhammad, Laiq Khan, Qudrat Khan, Shafaat Ullah, Irfan Sami, and Jong-Suk Ro. 2020. "Finite-Time Fast Dynamic Terminal Sliding Mode Maximum Power Point Tracking Control Paradigm for Permanent Magnet Synchronous Generator-Based Wind Energy Conversion System" Applied Sciences 10, no. 18: 6361. https://doi.org/10.3390/app10186361

APA Style

Zafran, M., Khan, L., Khan, Q., Ullah, S., Sami, I., & Ro, J. -S. (2020). Finite-Time Fast Dynamic Terminal Sliding Mode Maximum Power Point Tracking Control Paradigm for Permanent Magnet Synchronous Generator-Based Wind Energy Conversion System. Applied Sciences, 10(18), 6361. https://doi.org/10.3390/app10186361

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