Abstract

A modified control scheme based on the combination of online trained neural network and sliding mode techniques is proposed to enhance maximum power extraction for a grid connected permanent magnet synchronous generator (PMSG) wind turbine system. The proposed control method does not need the knowledge of the uncertainty bounds nor the exact model of the nonlinear system. Since the neural network is trained online, the time to estimate good weights can affect the dynamic performance of the process during the startup phase. Therefore an appropriate way to smoothly and explicitly accelerate the neural network rate of convergence during the startup phase is proposed. Furthermore, a flexible grid side voltage source converter control structure which can handle both grid connected and standalone modes based on conventional proportional integral (PI) control method is presented. Simulations are done in Matlab/Simulink environment to verify the effectiveness and assess the performance of the proposed controller. The results analysis shows the superiority of the proposed RBF neuro-sliding mode controller compared to a nonlinear controller based on sliding mode control method when the system undergoes parameter uncertainties.

1. Introduction

In recent years, a new challenge that seeks to reduce the effects of climate change has emerged in the world energy production domain. The exploitation of renewable energy sources has been intensified so as to tackle the growing energy demand with reduced environmental impact. Wind energy source, as a sustainable and renewable resource, has gained more interest for electrical power generation [1]. Wind energy conversion systems (WECSs) can be divided into two categories: fixed speed and variable speed. Fixed speed wind generators are directly connected to the electrical network but have limited controllability and range of operation [2]. Variable speed wind turbines (WTs) although equipped with additional components bring more advantages such as maximum power point tracking (MPPT) control strategy, increased power capture, and improved power quality [3]. The captured electrical energy can be injected into the utility grid or used to power up standalone loads. The electrical generator involved in the variable wind conversion process can be interfaced to the load by means of different topologies of controlled power electronic converters [4, 5].

Various investigations on wind turbine generators have shown that permanent magnet synchronous generator (PMSG) has many advantages compared to other kinds of electrical machines [6, 7] and therefore becomes an attractive choice for variable speed wind turbines. Meanwhile, the PMSG is a nonlinear dynamic and highly coupled system. In addition, the PMSG is subjected to many disturbances including change in wind speed and parameter uncertainties such as variations in stator resistance, stator inductance, rotor flux, and rotor inertia during operation. Hence, in order to tap efficiently power from PMSG based WT, it is important to take into consideration the above disturbances in the control design.

Linear controllers [814] have been studied and implemented. But the reduced range of operating points and poor robustness of these controllers with respect to machine uncertainties are their main limitations. Optimal tuning of control parameters for PMSG WECS using democratic joint operation algorithm has been reported in [15] to boost the performance of a linear controller. Meanwhile an alternative way to cope efficiently with the nonlinear nature and varying operation conditions of PMSG WECS is to exploit nonlinear control designs. Accordingly feedback [16, 17] and sliding mode [1825] techniques have been successfully processed. An association of passivity theory with sliding mode control was proposed in [26] to remove the PMSG WECS nonlinearities and cope with various uncertainties. However, the implementation of the above control algorithms required the knowledge of the mathematical model of the plant which is not accurate due to various assumptions considered. In addition, the knowledge of the upper limit of the system uncertainties, including parametric variations and external disturbances, is assumed. But the difficulty to obtain the exact system uncertainties bound in practical application may lead to unsatisfactory global control performance.

In searching appropriate solutions to the aforementioned shortcomings, artificial neural network (NN) control approaches have been explored. NNs are considered to be powerful tools for identification/nonlinear control by imitating cognitive expertise of the human [2730], thanks to their self-learning abilities and capabilities to approximate nonlinear functions. Some results from the application of NN based control schemes for PMSG WECSs can be found in [3134]. In [31], an online learning modified recurrent Elman NN control method is designed for a standalone PMSG WECS with practical implementation. An estimation algorithm is developed to observe the upper limit of the network approximation error. This estimation helped to improve the dynamic performance of the controller. In [32], online trained wavelet recurrent NN control method to improve the transient performance of a standalone PMSG WT has been reported and some experimental results were provided. The approaches in [31, 32] require extra memory to store previous information and high number of neurones which in turn induces large computation time. Another fact is that the total efficiency of the controlled PMSG WECS is limited since the MPPT algorithm, which is essential for variable speed WTs, was not addressed. In [33], an adaptive RBF NN control strategy is adopted to achieve tip speed ratio MPPT control of a PMSG WT with experimental results. The NN is trained offline. The parameters of the RBF NN are adjusted using a gradient method coupled to a modified particle swarm optimization algorithm. But the difficulty to obtain offline a training pattern which matches the exact behaviour of the closed loop system can hinder the performance of this scheme. In [34], sliding mode RBF NN control is applied to extract maximum power from PMSG WT. Simulation results demonstrated the robustness of the proposed method against wind speed variations. Nevertheless, the last two control methods assumed that the uncertainties of the system lie within a known boundaries. The main contributions of this paper can be summarized as follows:(i)A new control design of NN controller to tap maximum power from a grid connected PMSG wind turbine is investigated.(ii)Online trained (as opposed to offline trained methods, e.g., [33]) RBF neuro-sliding mode controller is selected to benefit from the simple structure and the properties of approximation capability of RBF NN and the robustness feature of sliding mode.(iii)An estimation of the upper limit of the system uncertainty bound which is not obtainable in practice is provided (as opposed to some existing works [33, 34]).(iv)A solution to reduce the control effort and accelerate the rate of convergence of the network during the startup phase (as opposed to the method in [30]) is proposed.(v)A PI control scheme for the grid side which is able to handle both grid connected and stand modes of operation is proposed.

The remainder of this paper is organized as follows. Section 2 focuses on the system modeling. In Section 3, the proposed hybrid control development is presented and the GSC control design is provided in Section 4. Computer simulations using a nonlinear controller based on sliding mode control method and the proposed sliding neural control approach are reported in Section 5. The conclusions of the paper are drawn up in Section 6.

2. PMSG WECS Modeling

Figure 1 presents an overview of the grid connected PMSG WECS under analysis. A wind turbine rotor is coupled to a PMSG shaft via a gearbox. The PMSG is controlled through a PWM machine side converter (MSC) to enhance MPPT. An electronic load controller (ELC), composed of a single switch PWM dc-to-dc converter feeding an electronic load, is used to provide a fixed dc voltage needed for PWM grid side converter (GSC) operation. The electronic load, working in a constant voltage mode, is designed such that when the duty cycle of the dc-to-dc converter is unity (during fault or over generation), it should consume the maximum output power of the generator. A LCL filter is introduced to improve the quality of the power injected into the local AC load or the grid by lessening the harmonics produced by the GSC switching effects. A circuit breaker is used in case of grid connected/standalone transition.

2.1. Turbine Modeling

The mechanical power extracted by the WT from the kinetic energy of the wind is expressed in terms of the air density , the area swept by the blades ( is the blade radius), the power coefficient , and the wind speed as follows [35]:The power coefficient is a function of the tip speed ratio and the blade pitch angle and can be stated asThe tip speed ratio is expressed aswhere represents the gearbox ratio. and are rotational speed and torque developed by the turbine at the blade side of the gearbox while and are the machine side shaft speed and mechanical torque, respectively. is computed by combining (1) and (3) asThe turbine equation of motion referred to as the generator side of the gearbox is given bywith and being the total rotating parts inertia and viscous factor, respectively. is the electromagnetic torque of the generator and is defined later.

When the wind speeds belong to region 2 [36, 37], i.e., between the cut-in and the rated values, the MPPT algorithm should be enabled to increase the efficiency of the turbine captured power. Therefore the power coefficient should be kept at its maximum value for each wind speed in this region. This is possible if the turbine is working at the optimum tip speed ratio which in turn leads to generator optimum speed as follows:Hence the optimum power caption is achieved if the generator speed is steered to the optimum speed . This MPPT scheme, known as tip speed ratio MPPT, is adopted in the present work. Thus the target optimum power of the generator is computed aswhere is the optimal torque constant.

2.2. Modeling of PMSG Associated with MSC and dc-Link

The nonlinear dynamics of a surface mounted PMSG coupled with PWM voltage source rectifier and PWM dc-to-dc converter (see Figure 1) in the synchronously rotating coordinates system is described as follows [3840]:where denote the stator current components and and are the stator windings resistance and inductance. is the number of pole pairs. are the stator flux components and is the rotor flux. is the electrical power delivered by the generator. are the stator voltage components while is the voltage across the dc-link capacitor . is the electronic load resistance. is the duty cycle of the dc-to-dc converter switch. and are the GSC current and active power, respectively. state the to reference frame transformation of the PWM duty ratios of the MSC switching devices.

2.3. GSC and LCL Filter Dynamics

The mathematical model of the GSC and LCL filter in the coordinates system is expressed as follows [41]:where , , and are inverter current, filter capacitor voltage, and grid side current components. and represent, respectively, the local AC load current and injected grid current components. and denote the inverter output voltage components. and are the duty ratios of the GSC. is the angular frequency of the grid voltage. , , and denote, respectively, the LCL filter inverter side inductance, capacitance, and grid side inductance. The filter inductors have internal resistances and .

3. Proposed RBF Neuro-Sliding Mode Control Design

The MSC, acting as voltage source rectifier, is controlled to efficiently extract the maximum power from the wind turbine while the dc-to-dc converter control task is to track the voltage across the dc-link capacitor to the desired level. It is important to regulate the dc voltage for the efficient working of the GSC. The -axis current reference is set to zero in order to obtain unity power factor of the PMSG while the generator speed reference is equal to the optimal rotor speed defined earlier for MPPT.

In [30] authors proposed a combined online trained RBF neural network and sliding mode control algorithm for a class of SISO nonlinear uncertain systems. In the present study the same approach is modified and exploited. Unless otherwise specified, the proposed control development is conducted assuming the following:(A1) The state variables as well as the system disturbances are bounded; i.e., they belong to a compact subset .(A2) The rate of change of the system parameters is slow compared to the fast dynamics of the system state variables.(A3) The upper bound of the system uncertainties is unknown.

Before processing the proposed RBF neuro-sliding mode control scheme, let us derive each state variable to be control (i.e., , and ) until its control input explicitly appears and thus rearrange the -axis current dynamic in (8) asThe term has been added to represent the unknown uncertainties of the system. Since the mechanical torque is unmeasurable, replace it by its optimal expression . Then consider the electromagnetic torque expression in (10) and derive (5) with respect to time as follows: represents the unknown uncertainties of the system.

From (14), rewrite the dc-link voltage dynamic as ( and state, respectively, the square dc voltage, and the system uncertainties):Now define the d-axis current, the rotor speed, and square dc voltage tracking errors as , , and , respectively (the reference signals are marked with ). Let also the associated sliding variables be , , and such that (, and )Derive the sliding surfaces with respect to time and obtainwith , and .

Consider systems (21), (23), and (25) with proposed respective controllers in (31), (32), and (33)where , , and are design parameters. Thus (28), (29), and (30) becomeFrom (34), (35), and (36) one can conclude that , , and converge exponentially to 0.

Remark 1. Considering the assumption that the functions , , , , , and along with the terms , , , , , and are unknown, the controllers in (31), (32), and (33) can not be implemented in real operation condition. Consequently three sets of RBF neuro-sliding mode controllers are employed to estimate the control signals , , and .

Remark 2 (the conditions , , and can be verified). Compute the proposed approximations of the desired control actions , , and as follows:where , , and are the NNs input vectors. , , and denote the RBF NNs. , , and are estimations of optimal weight vectors , , and . Bearing in mind Remark 2, the terms , , and , which help to improve the neural networks rate of convergence in the presence of uncertainties, are described as follows:, , and are estimates of the positive unknown limits of the NNs uncertainties and are computed aswhere , , and , , . For number of nodes in the hidden layer, the functions and are described as follows:where is a Gaussian activation function. and are, respectively, the center and the width of the th node. is the projection function on the compact . The weight vectors are initialized in a uniformly random manner while the centers and the widths are selected according to the clustering method as [42]where and state, respectively, the lower and the upper bounds of the input vector .
The proof of the convergence of the above NNs controllers is found in [30].

Remark 3. During the startup phase, given that the weights vector is randomly initialized, the network may take a relatively long time to converge. This implies that the network uncertainties are more significant during this phase. Unfortunately the estimations of the network uncertainties bounds given in (43), (44), and (45) start from 0 before increasing to sufficient values. An obvious solution is to take a too high estimation rate . But a too large rate may excite high frequency unmodeled dynamics, maintain higher control effort, and render the system unstable.
To address this issue, a transformation for the estimation of uncertainties bound is proposed. This transformation seeks to steer the tracking error to the origin and also allow the tracking process to be explicitly controllable so that the rate of convergence during the startup phase can be smoothly pre-specified and faster. In concrete terms, definewhere is the transformed and , a smooth decreasing function that satisfies is the maximum value of . For example, the functionis a good candidate and may be used in (49).
Considering the aforementioned transformation, the bias terms in (40), (41), and (42) becomewhere and , , are specified by the designer.
It is well known that online identification and adaptation of a time varying parameter in a controller contribute to relax the requirement of the system parameter knowledge and help to improve the control performance [43]. Such is the purpose of the following rotor flux identifier. Substitute by in (5) and rearrange it asBased on (57), a generator speed observer can be written as ( is a constant gain)Computing the dynamic of observation error yieldswhere and denote the estimated quantities while and are the error quantities. Assuming that the parameter identifiability condition is verified, the rotor flux estimation law can be processed as ( and is Laplace variable)The proof of convergence of the identification algorithm in ((61)–(63)) can be found in [24].

The schematic description of the proposed control algorithm is presented in Figure 2. In this figure, graphs (a), (b), and (c) show the equivalent control systems while graph (d) exhibits the interface of the control inputs with PWM units. is the generator shaft position.

4. Control Design for GSC

The GSC, acting as a voltage source inverter, is controlled to track the active and reactive power of the grid by using grid current regulator. Conventional PI control method is exploited in the present study to achieve the control objective. Details concerning the selection of PI gains are described in Appendix. Figure 3 illustrates the proposed GSC control structure. is the grid voltage angular position estimated by a PLL for grid synchronization.

The references of the grid current components are computed from the reference active power and reference reactive power as follows:In order to allow the GSC to operate at unity power factor, is set to zero. The reference active power is derived directly from the optimum power captured from the WT as shown in (65).

Remark 4. In standalone mode, the GSC is controlled to regulate the amplitude and frequency of the AC load voltage. From the grid current control structure in Figure 3, it is possible to switch to standalone mode by disconnecting the grid and setting manually the reference components ( and ) of the voltage across the filter capacitor and the value of the load voltage angular frequency for reference frame transformation. The switches in Figure 3 together with the circuit breaker are dedicated to this end. That is an interesting feature of the proposed inverter control structure.

5. Numerical Simulation Results

Numerical simulations of the studied WECS (see Figure 1 with control structures in Figures 2 and 3) are carried out in Matlab/Simulink software to assess the performance of the proposed RBF neuro-sliding mode controller (RBF-NSMC). Simulations based on a sliding mode control (SMC) method [24] are conducted and the results are presented to compare the capabilities of the proposed RBF-NSMC. The WECS parameters [44] and the numerical values of the control parameters are listed in Tables 1 and 2, respectively. It is worth mentioning that the PMSG can directly be connected to the WT. However, a gearbox ratio is selected to comply with realistic transmission condition. An inductive RL load ( and ) is used to materialize a general local AC load. More details on the selection criteria of the LCL filter parameters are given in Appendix. Each neural network is built up with four neurons in the hidden layer.

The analysis of the system response is conducted considering three operation conditions. Firstly the effectiveness of the proposed RBF-NSMC is investigated under nominal WECS parameters (see Figures 49). Secondly the robustness of the controller is examined under WECS parameters uncertainties (see Figures 1012). The wind speed is time varying in both cases and its behaviour is shown in Figure 4. Finally the performance of the control system is evaluated under transition from grid connected to standalone mode during low wind speeds (see Figures 13 and 14).

Case 1 (test under varying wind speed with WECS nominal parameters). It can be observed on Figures 47 that the control objective of the generator side converter and dc-link is achieved since the generator axis current is rapidly set to zero (see Figure 5), the dc voltage is regulated to the desired reference value, and the maximum power is accurately harnessed from the wind (see Figure 6). For this test, one can remark clearly on Figure 7 that the results provided by both proposed RBF-NSMC and SMC methods are practically similar. Figures 8 and 9 exhibit the behaviour of the grid side signals when the proposed control approach is maintained on the machine side and dc-link. It is noticeable from Figure 8(a) that the grid current components are steered to their estimated references with small response time. Figure 8(b) shows that the regulation of the voltage across the LCL filter capacitor is satisfactory. This voltage regulation is necessary for local AC load feeding and grid synchronization. The sinusoidal waveforms of the filter capacitor voltage and the total grid plus local load current are depicted in Figure 9. From the above observations one can undoubtedly conclude that the proposed control scheme is effective and the extracted optimum power is transferred to the load with unity power factor.

Case 2 (test under varying wind speed with WECS parameters uncertainties). The graphs on Figures 1012 show the results obtained when uncertainties are introduced on the process parameters at time . These perturbations consist of , , , and sudden changes on the nominal values of stator resistance, stator inductance, rotor flux, and rotor inertia, respectively. From Figure 10, one can note that the proposed controller presents fast dynamic response compared to that of the SMC one. The absolute values of tracking errors representing maximum deviations from the optimal curve, for different time periods, are picked up from Figure 10(b) and presented in Table 3. It comes out from this table that the deviations under the proposed RBF-NSMC method are always smaller than those provided by the SMC. This demonstrates that the proposed control scheme features higher robustness with respect to varying wind speed and parameter disturbances. It is shown in Figure 11(a) that the transformed estimation of the unknown uncertainties bound starts with a high value to accelerate the process and converges to a low value to reduce the energy consumed by the controller and also prevent from unwanted high frequency unmodeled dynamics. The resulting smooth control actions for both MSC and GSC are shown in Figure 11(b). These control actions are very significant since they are not saturated and therefore are favorable for interfacing with PWM units. This remark supports the fact that the proposed RBF-NSMC is feasible for real time implementation. The total harmonic distortions (THDs) of the LCL filter capacitor voltage and injected grid current are both less than which is the upper limit as seen in Figure 12. Although this study is focused on the grid connected mode, one may say that, for wind speeds lower than those depicted in the used wind profile, one can switch to standalone mode as stated in Remark 4. This helps not to exceed the THD maximum value on the injected grid current.

Case 3 (performance evaluation under grid connected to standalone mode transition). In this test, the system starts with the same conditions as in Case 2. But at time , the mean value of the wind profile changes to low value (from m·s−1 to m·s−1, see Figure 13). Simultaneously the system is switched to standalone mode (see Remark 4). It is observed on Figure 14(a) that the global stability of the system is preserved and better transient performance is kept by the proposed RBF-NSMC as compared to the SMC. It is seen on Figure 14(b) that the filter capacitor voltage is still well regulated by the proposed grid side PI controller despite the grid disconnection. It is also noted that, after this grid disconnection, the injected grid current vanishes and only the portion absorbed by the local load remains. These observations demonstrate the flexible feature of the proposed grid side control structure.

6. Conclusions

A control scheme of a grid connected variable speed PMSG WECS is investigated. The control approach based on online learning RBF-sliding network has been proposed to extract efficiently maximum power without the knowledge of neither the exact model of the nonlinear system nor the requirement of the upper limit of the system disturbances. A transformation of the estimated system uncertainties bound for explicitly accelerating the rate of convergence of the NN during startup phase has been proposed. Simulation results have shown the effectiveness and the feasibility of the proposed control scheme. The proposed RBF-NSMC exhibited higher robustness with fast dynamic response face to wind speed variations and WECS parameters changes. The proposed grid side inverter control structure is flexible so as to switch between grid connected mode and standalone mode in order to cover both high and low wind speeds.

Appendix

A. Selection of GSC PI Control Gains

In this reasoning and represent proportional and integral gains, respectively, and quotes the error signal, while is a signal with being the corresponding reference. Rearrange (19) and (20) in the following form: and are the transformed control inputs. And compute the open loop transfer functions of (A.1) and (A.3)where is the open loop time constant. Now select the correctors in the form of Proceed to the compensation of the open loop time constant by choosing asand calculate the corrected closed loop transfer functions in the form of first-order systems aswhere is the closed loop time constant. The practical settling time with criteria () of a first-order system is equal to three times its time constant. On the other hand the settling time of the closed loop system should be at least equal to the open loop settling time; that is,Hence the outputs of the closed loop systems (A.1) and (A.3) will exponentially converge to their reference values if conditions (A.8) and (A.10) are satisfied. Thus the reference values of components for the voltage across the LCL filter capacitor are expressed asAdopt the same methodology for the GSC currents control loops and obtainFor the control design of the LCL filter capacitor voltage, rewrite (17) and (18) in the form of The corresponding transfer functions are written asLet the correctors beThen calculate the closed loop transfer functions and arrange them in the form of second-order systems as and are, respectively, the system oscillation damping factor and undamped natural angular frequency. Thus the reference currents of the GSC are computed as

B. LCL Filter Design

Grid side inverter rated power , nominal grid line to line voltage , switching frequency , and grid frequency . In this study the calculation of LCL filter parameters is based on the methodology presented in [45].

Base Values Definition. The base impedance , inductance , capacitance , and nominal (rms) converter current are computed, respectively, as follows:where is the grid angular frequency.

Converter Side Inductance Selection. The converter side inductance is chosen according to the following expression:where is the maximum ripple current which must lie between and of the nominal converter current.

Filter Capacitance Selection. The filter capacitance should be less than of the base value; i.e., .

Grid Side Inductance Selection. The grid inductance is selected as is the ratio between converter side and grid side inductances and is the harmonics attenuation of the phase current. In this study the parameters have been chosen by considering attenuation of maximum converter ripple current, base capacitance, and . The resulting LCL filter parameters are , , and . Thus the filter resonant frequency is well between the desired range . The inductors internal resistances have been taken as and .

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

There are no conflicts of interest.