Next Article in Journal
Impact of the Taxes on Used Nuclear Fuel on the Fuel Cycle Economics in Spain
Next Article in Special Issue
Forces and Moments on Flat Plates of Small Aspect Ratio with Application to PV Wind Loads and Small Wind Turbine Blades
Previous Article in Journal
Harmonic Injection-Based Power Fluctuation Control of Three-Phase PV Systems under Unbalanced Grid Voltage Conditions
Previous Article in Special Issue
Simulating Dynamic Stall Effects for Vertical Axis Wind Turbines Applying a Double Multiple Streamtube Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Tuning of the PI Controller Parameters of a PMSG Wind Turbine to Improve Control Performance under Various Wind Speeds

1
School of Electrical and Computer Engineering, Seoul National University, Seoul 151-744, Korea
2
School of Electrical Engineering, Kookmin University, Seoul 136-702, Korea
*
Author to whom correspondence should be addressed.
Energies 2015, 8(2), 1406-1425; https://doi.org/10.3390/en8021406
Submission received: 2 January 2015 / Accepted: 2 February 2015 / Published: 13 February 2015
(This article belongs to the Collection Wind Turbines)

Abstract

:
This paper presents a method to seek the PI controller parameters of a PMSG wind turbine to improve control performance. Since operating conditions vary with the wind speed, therefore the PI controller parameters should be determined as a function of the wind speed. Small-signal modeling of a PMSG WT is implemented to analyze the stability under various operating conditions and with eigenvalues obtained from the small-signal model of the PMSG WT, which are coordinated by adjusting the PI controller parameters. The parameters to be tuned are chosen by investigating participation factors of state variables, which simplifies the problem by reducing the number of parameters to be tuned. The process of adjusting these PI controller parameters is carried out using particle swarm optimization (PSO). To characterize the improvements in the control method due to the PSO method of tuning the PI controller parameters, the PMSG WT is modeled using the MATLAB/SimPowerSystems libraries with the obtained PI controller parameters.

1. Introduction

The penetration of renewable energy resources has seen consistent growth in recent years [1,2]. In particular, the use of wind power has grown considerably, and as of the end of 2011 significant commercial installations have appeared in 80 countries, totaling approximately 240 GW; the average annual cumulative growth rate of wind power over the past 15 years has been approximately 28% [2,3]. Due to this increase in capacity, the specific technical requirements for wind turbines (WTs) have led to control issues [4].
In this paper we focus on improving the control of WTs by tuning the controller parameters depending on the wind speed. The wind speed is an important factor that determines the operating conditions of WTs, and may vary arbitrarily and in an unpredictable manner; controller parameters should therefore be tuned for various wind speeds. In this paper, the WTs considered are permanent magnet synchronous generator (PMSG) turbines. Along with doubly fed induction generator (DFIG) WTs [2,5], PMSGs are the most common type of WT; however, PMSG WTs are preferable, since the former have lower maintenance costs [2,6]. The maintenance cost is an important factor since most of the large-scale wind farms have been sited on off-shore recently. This phenomenon is due to its geographical limitation and environmental regulations for on-shore site. The controller is designed as a proportional-integral (PI) controller since it is the most common controller used in industrial applications today [7,8].
Many literatures on WT control performance issue have been published [4,9,10,11,12,13,14,15]. A comparison between a PMSG WT and a DFIG WT was studied from the stability perspective in [13], but no method for tuning controller parameters to improve control performance was proposed. In reference [4], two different PMSG WT control modes were compared with each other and analyzed from the power coefficient and stability perspectives. However, it did not contain any theoretical analysis such as small-signal analysis which gives more precise result of analysis. For precise analysis of the control performance, some early works have investigated using optimization algorithms and small-signal analysis. For instances, in reference [9], DFIG controller parameters were tuned by particle swarm optimization (PSO) and in reference [10], small-signal modeling and analysis was implemented and PI controller parameters were adjusted by investigating traces of eigenvalues. For other instances, a genetic algorithm (GA) was used in reference [11] to adjust the controller parameters of a PMSG WT and a torque compensation strategy for PMSG WTs to improve the system stability was presented in reference [12]. A major demerit of references [9,10,11,12], which is the main focus in the paper, is that they have only considered a fixed wind speed and that of references [10,12] is that they did not apply optimization algorithms, so they cannot guarantee that the designed controller parameters are optimized to their objective. A small-signal stability analysis of a DFIG WT was presented in references [14,15]. A controller parameter adjustment for different wind speeds was implemented in them, but the number of considered wind speed points for designing the controller parameters is insufficient to apply to various wind speeds. Overall, the small-signal modeling and analysis of a PMSG WT has been rarely researched [10].
In this paper, PI controller parameters of a PMSG WT for various wind speeds are designed to improve control performance. To this end, small-signal modeling is implemented for analyzing the control performance of a PMSG WT. Eigenvalues which can be obtained from a system matrix of the small-signal model are coordinated by adjusting the PI controller parameters. The parameters to be tuned are selected by using participation factors to make the problem simple. The process of adjusting the parameters is conducted by using a PSO algorithm which uses the heuristic-based swarm intelligence concept. It is known to effectively solve large-scale nonlinear optimization problems and has fast convergence characteristics [16]. The effectiveness of the PI controller parameters obtained by the proposed method is verified by a simulation which applies them to a PMSG WT modeled in MATLAB/SimPowerSystems.

2. Small-Signal Modeling

Figure 1 shows the system to be analyzed, which contains a PMSG WT integrated with an infinite bus. It is composed of the WT drive train, a PMSG, a machine-side converter (MSC), a DC-link capacitor, a GSC, an RL-filter, and a transformer. In this section, each component of the system will be formulated into a state-space model, which will then be used to find the eigenvalues.
Figure 1. Schematic diagram of the system.
Figure 1. Schematic diagram of the system.
Energies 08 01406 g001

2.1. Wind Turbine Model

The power that can be extracted from the wind by a WT can be expressed as [17]:
P w = 1 2 C p ( λ , β ) ρ π r 2 v w 3
λ = ω t r v w
where λ is the tip-speed ratio, β is the pitch angle, ρ is the air density, r is the blade radius, vw is the wind speed, ωt is the rotational speed of the WT, and Cp is the power coefficient of WT which is a function of λ and β. Cp is given as [17]:
C p = 0.73 ( 151 λ i 0.58 β 0.002 β 2.14 13.2 ) e 18.4 / λ i
λ i = ( 1 λ 0.02 β + 0.003 β 3 + 1 ) 1 .
The pitch angle β is controlled only when vw exceeds the rated wind speed to maintain the rated active power of the PMSG. Otherwise, it is kept constant at β = 0˚ which is the only considered value of β since only the operating conditions under the rated wind speed will be treated in this paper. Hence, it does not influence the small-signal modeling.
For a given vw and β, there is an optimal rotational speed of the WT ωt,opt which gives an optimal tip-speed ratio λopt [18,19]. It is essential to maintain ωt = ωt,opt to maximize Pw which is an objective of the maximum power point tracking (MPPT) control. This means that it is to make Cp to the maximum value Cp,max, so the maximum power that can be extracted from the wind is given as [17]
P w , m a x = K o p t ω t , o p t 3
K o p t = 1 2 C p , m a x ρ π r 2 ( r λ o p t ) 3 .
It can be implemented by controlling the electrical rotational speed of the generator rotor ωe, which has a following relationship with ωt [17]:
ω e = N p p ω m = N g r N p p ω t
where Npp is the number of pole pairs in the PMSG, Ngr is the gear ratio, and ωm is the mechanical rotational speed of the generator rotor. ωm accelerates or decelerates respect to the WT driving torque equation [20]:
T e T m = J e q d ω m d t
where Te and Tm are the electromagnetic and the mechanical torque of the generator, respectively, and Jeq is the total equivalent inertia referred to the generator which contains the turbine and the generator inertia.

2.2. Permanent Magnet Synchronous Generator Model

Transforming a three-phase stator current imabc into the machine-side dq-reference frame aligned with the electrical rotor position θe, the voltage and torque equations of the salient-pole PMSG can be expressed as [20]:
v m d = R s i m d + L d d i m d d t ω e L q i m q v m q = R s i m q + L q d i m q d t + ω e L d i m d + ω e λ p m
T e = 3 2 N p p { λ p m i m q + ( L d L q ) i m d i m q }
where Rs is the stator resistance, λpm is the permanent magnetic flux, and vmd, vmq, imd, imq, Ld, and Lq are the d- and q-component of stator voltage, current, and inductance, respectively.

2.3. Machine-Side Converter Controller Model

The MSC control structure is shown in Figure 2. It consists of three PI controllers each of which controls the current imd, imq, and the active power output Pout, respectively. Superscript * denotes the reference values of the PI controller inputs. From Equations (5)–(7), (9), we may determine expressions for the relevant MSC controller as follows:
i m d * = 0   A R s i m d + L d d i m d d t = K p 1 ( i m d * i m d ) + K i 1 φ 1 d φ 1 d t = i m d * i m d P o u t * = K o p t ( ω e N p p N g r ) 3 P o u t = 3 2 ( v s d i s d + v s q i s q ) = 3 2 ( v s d i g d + v s q i g q ) = 3 2 v s d i g d i m q * = K p 2 ( P o u t P o u t * ) + K i 2 φ 2 d φ 2 d t = P o u t P o u t * R s i m q + L q d i m q d t = K p 3 ( i m q * i m q ) + K i 3 φ 3 d φ 3 d t = i m q * i m q
where Kp1Kp3 and Ki1Ki3 are the p- and i-gains of the PI controllers, vsd, vsq, isd, and isq are the d- and q-axis voltage and current at the terminal of RL-filter, igd and igq are the GSC terminal d- and q-axis current, respectively, and φ1φ3 are the intermediate state variables introduced by the PI controllers. Imposing i m d * = 0 is an optimal approach to avoid the nonlinearity of Equation (10) and to minimize current coupling [21,22]. This reduces Equation (10) to [23]:
T e = 3 2 N p p λ p m i m q .
Regulating i*md = 0 guarantees minimum stator currents for a given torque, which enhances the mechanical efficiency [22]. When calculating Pout, since it passes through same current path, isd and isq can be set equal to igd and igq, respectively. Since the grid-side dq-reference frame is aligned with θs, which is the angular position of the RL-filter terminal voltage, vsq = 0 which allows to eliminate vsq term in Pout shown in Equation (11). The dynamic process of the converter can be ignored since it is much faster than electromagnetic and mechanical dynamic process, so it can be assumed that v*md = vmd, v m q * = vmq [10].
Figure 2. Control structure of MSC.
Figure 2. Control structure of MSC.
Energies 08 01406 g002

2.4. DC-Link Model

Figure 3 shows the power flow through the DC-link part from the PMSG to the power grid where Cdc is the capacitance of the DC-link capacitor, PP and PM are the active powers from the PMSG and the MSC, respectively, PC, PG, and PF are the active powers flow into the DC-link capacitor, the GSC, and the RL-filter, respectively, and Vdc and Idc are the DC-link voltage and current, respectively. Assuming that power losses in the MSC, the GSC, the RL-filter, and the transformer can be neglected, the active power transferred from the PMSG to the grid can be expressed as:
P P = 3 2 ( v m d i m d v m q i m q ) = P M P C = V d c I d c = V d c C d c d V d c d t = P M P G P G = P F = P o u t .
Figure 3. DC-link power flow.
Figure 3. DC-link power flow.
Energies 08 01406 g003

2.5. Grid-Side Converter Controller Model

The voltage balance across the RL-filter shown in Figure 1 can be expressed as [20]:
v g d = R f i g d + L f d i g d d t ω s L f i g q + v s d v g q = R f i g q + L f d i g q d t + ω s L f i g d + v s q
where vgd and vgq are the GSC terminal d- and q-axis voltage, respectively, Rf and Lf are the resistance and the inductance of the RL-filter, respectively, and ωs is the angular frequency of the RL-filter terminal voltage.
Figure 4 shows the structure of the GSC controller. It includes four PI controllers, which control Vdc, igd, igq, and the reactive power output Qout, respectively. As that of the MSC controller, superscript * denotes for the reference values of the PI controller inputs. Based on Equation (14) and the control structure shown in Figure 4, the equations relevant to the GSC controller can be expressed as follows:
V d c * = 5400  V i g d * = K p 4 ( V d c V d c * ) + K i 4 φ 4 d φ 4 d t = V d c V d c * R f i g d + L f d i g d d t = K p 5 ( i g d * i g d ) + K i 5 φ 5 d φ 5 d t = i g d * i g d Q o u t * = 0  VAR Q o u t = 3 2 ( v s q i s d v s d i s q ) = 3 2 ( v s q i g d v s d i g q ) = 3 2 v s d i g q i g q * = K p 6 ( Q o u t Q o u t * ) + K i 6 φ 6 d φ 6 d t = Q o u t Q o u t * R f i g q + L f d i g q d t = K p 7 ( i g q * i g q ) + K i 7 φ 7 d φ 7 d t = i g q * i g q
where Kp4Kp7 and Ki4Ki7 are the p- and i-gains of the PI controllers, respectively, and φ4φ7 are the intermediate state variables introduced by the PI controllers. Imposing Q o u t * = 0 is to maintain the unity power factor. As mentioned in the previous section, isd and isq are equal to igd and igq, respectively, and vqs = 0, so Qout can be simplified as shown in Equation (15).
Figure 4. Control structure of GSC.
Figure 4. Control structure of GSC.
Energies 08 01406 g004

2.6. Power Grid Model

As shown in Figure 1, the PMSG WT is connected to the infinite bus via the transformer and the transmission line. To establish the equations of the power grid, assuming that the resistance of the transformer is neglected, the circuit on the right side of node s in Figure 1 can be simplified as Figure 5 where xT and xL are the reactances of the transformer and the transmission line, respectively, igdq is the GSC terminal dq-axis current, vsdq is the RL-filter terminal dq-axis voltage, vidq is the infinite bus dq-axis voltage, |vsdq| and |vidq| are the magnitude of vsdq and vidq, respectively, and δ is the angle of the infinite bus voltage relative to the reference angle. Since the reference angle is aligned to that of vsd, the equation of the current flows into the power grid can be calculated as:
i g d q = | v s d q | sin 0 | v i d q | sin δ x T L + j | v s d q | cos 0 + | v i d q | cos δ x T L
where xTL is the summation of xT and xL. Since the power factor is maintained to unity, igdq is in phase with vsdq. And this makes the real part and the imaginary part of igdq to be igd and igq, respectively. Thus igd and igq can be written as:
i g d = | v i d q | sin δ x T L i g q = v s d + | v i d q | cos δ x T L
Figure 5. One-line diagram of the power grid.
Figure 5. One-line diagram of the power grid.
Energies 08 01406 g005

2.7. Complete Small-Signal Model

To analyze the system stability following a change in the wind speed, and to evaluate the effectiveness of the controller parameters for the whole system, it is necessary to establish a complete small-signal model of the PMSG WT. By linearizing Equations (1)–(17), small-signal state-space equations of the PMSG WT can be obtained as follows:
x ˙ = A x + B u y = C x + D u x = [ ω e   i m d   i m q   V d c   i g d   i g q   φ 1   φ 2   φ 3   φ 4   φ 5   φ 6   φ 7 ] T u = [ v w   i m d *   V d c *   Q o u t * ] T y = [ P o u t   Q o u t ] T
where the individual elements of matrix A, which are derived as explained in the Appendix, are given at the bottom of the page and the subscript 0 appears in the elements of A denotes for the initial value. The matrices B, C, and D are unnecessary to calculate since the main purpose of the small-signal modeling is to acquire eigenvalues of A. It is sufficient to analyze the system stability only with the eigenvalues of A since they are identical to the system poles before cancellation of common poles and zeroes in the transfer function [24].
A = [ a 1 , 1 a 1 , 2 a 1 , 3 0 0 0 0 0 0 0 0 0 0 0 a 2 , 2 0 0 0 0 a 2 , 7 0 0 0 0 0 0 a 3 , 1 0 a 3 , 3 0 a 3 , 5 a 3 , 6 0 a 3 , 8 a 3 , 9 0 0 0 0 a 4 , 1 a 4 , 2 a 4 , 3 a 4 , 4 a 4 , 5 a 4 , 6 a 4 , 7 a 4 , 8 a 4 , 9 0 0 0 0 0 0 0 a 5 , 4 a 5 , 5 0 0 0 0 a 5 , 10 a 5 , 11 0 0 0 0 0 0 a 6 , 5 a 6 , 6 0 0 0 0 0 a 6 , 12 a 6 , 13 0 a 7 , 2 0 0 0 0 0 0 0 0 0 0 0 a 8 , 1 0 0 0 a 8 , 5 a 8 , 6 0 0 0 0 0 0 0 a 9 , 1 0 a 9 , 3 0 a 9 , 5 a 9 , 6 0 a 9 , 8 0 0 0 0 0 0 0 0 a 10 , 4 0 0 0 0 0 0 0 0 0 0 0 0 a 11 , 4 a 11 , 5 0 0 0 0 a 11 , 10 0 0 0 0 0 0 0 a 12 , 5 a 12 , 6 0 0 0 0 0 0 0 0 0 0 0 a 13 , 5 a 13 , 6 0 0 0 0 0 a 13 , 12 0 ]
a 1 , 1 = 9.305 ρ π r 2 v w 0 3 N p p 2 2 J e q ω e 0 2 ( 217.965 v w 0 2 N p p 2 N g r 2 r 2 ω e 0 2 42.092 v w 0 N p p N g r r ω e 0 + 1 ) exp ( 18.4 v w 0 N p p N g r r ω e 0 0.0552 ) ,   a 1 , 2 = 3 N p p 2 2 J e q ( L d L q ) i m q 0 ,
a 1 , 3 = 3 N p p 2 2 J e q { λ p m + ( L d L q ) i m d 0 } ,   a 2 , 2 = K p 1 R s L d ,   a 2 , 7 = K i 1 L d ,   a 3 , 1 = 3 K o p t K p 2 K p 3 ω e 0 2 L q N p p 3 N g r 3 ,   a 3 , 3 = K p 3 R s L q ,   a 3 , 5 = 3 K p 2 K p 3 2 L q ( v s d 0 x T L 2 i g d 0 2 v s d 0 + x T L i g q 0 ) ,
a 3 , 6 = 3 K p 2 K p 3 x T L i g d 0 2 L q ,   a 3 , 8 = K i 2 K p 3 L q ,   a 3 , 9 = K i 3 L q ,   a 4 , 1 = 3 i m q 0 2 C d c V d c 0 { i m d 0 L q ( 3 K o p t K p 2 K p 3 ω e 0 2 N p p 3 N g r 3 + L d i m d 0 + λ p m ) } ,   a 4 , 2 = 3 2 C d c V d c 0 ( K p 1 i m d 0 v m d 0 ω e 0 L d i m q 0 ) ,
a 4 , 3 = 3 2 C d c V d c 0 ( ω e 0 L q i m d 0 v m q 0 + K p 3 i m q 0 ) ,   a 4 , 4 = 3 ( v m d 0 i m d 0 + v m q 0 i m q 0 + v s d 0 i g d 0 ) 2 C d c V d c 0 2 ,   a 4 , 5 = 3 2 C d c V d c 0 { 3 K p 2 K p 3 i m q 0 2 ( x T L 2 i g d 0 2 v s d 0 + x T L i g q 0 v s d 0 ) v s d 0 + x T L 2 i g d 0 2 v s d 0 + x T L i g q 0 } ,
a 4 , 6 = 3 x T L i g d 0 2 C d c V d c 0 ( 1 3 K p 2 K p 3 i m q 0 2 ) ,   a 4 , 7 = 3 K i 1 i m d 0 2 C d c V d c 0 ,   a 4 , 8 = 3 K i 2 K p 3 i m q 0 2 C d c V d c 0 ,   a 4 , 9 = 3 K i 3 i m q 0 2 C d c V d c 0 ,   a 5 , 4 = K p 4 K p 5 L f ,   a 5 , 5 = K p 5 R f L f ,   a 5 , 10 = K i 4 K p 5 L f ,   a 5 , 11 = K i 5 L f ,
a 6 , 5 = 3 K p 6 K p 7 x T L 2 i g d 0 i g q 0 2 L f ( v s d 0 + x T L i g q 0 ) ,   a 6 , 6 = 1.5 K p 6 K p 7 ( x T L i g q 0 v s d 0 ) R f K p 7 L f   a 6 , 12 = K i 6 K p 7 L f ,   a 6 , 13 = K i 7 L f ,   a 7 , 2 = 1 ,   a 8 , 1 = 3 K o p t ω e 0 2 N p p 3 N g r 3 ,   a 8 , 5 = 3 2 ( v s d 0 x T L 2 i g d 0 2 v s d 0 + x T L i g q 0 ) ,
a 8 , 6 = 3 x T L i g d 0 2 ,   a 9 , 1 = 3 K o p t K p 2 ω e 0 2 N p p 3 N g r 3 ,   a 9 , 3 = 1 ,   a 9 , 5 = 3 K p 2 2 ( v s d 0 x T L 2 i g d 0 2 v s d 0 + x T L i g q 0 ) ,   a 9 , 6 = 3 K p 2 x T L i g d 0 2 ,   a 9 , 8 = K i 2 ,   a 10 , 4 = 1 ,   a 11 , 4 = K p 4 ,   a 11 , 5 = 1 ,
a 11 , 10 = K i 4 ,   a 12 , 5 = 3 x T L 2 i g d 0 i g q 0 2 ( v s d 0 + x T L i g q 0 ) ,   a 12 , 6 = 3 ( x T L i g q 0 v s d 0 ) 2 ,   a 13 , 5 = 3 K p 6 x T L 2 i g d 0 i g q 0 2 ( v s d 0 + x T L i g q 0 ) ,   a 13 , 6 = 3 K p 6 ( x T L i g q 0 v s d 0 ) 2 1 ,   a 13 , 12 = K i 6 .

3. Optimization Algorithm

Here PSO is employed due to the effectiveness in solving optimization problems. A brief introduction to the PSO algorithm follows, together with the problem formulation.

3.1. Particle Swarm Optimization Algorithm

PSO is part of the swarm intelligence family of algorithms, which are known to effectively solve large-scale nonlinear optimization problems [16]. It exhibits excellent performance in searching for the global optimum because it can diversify the swarm using a stochastic velocity term [25]. Particles look for the optimal solution in search space based on position, which is learnt heuristically and iteratively. The movement of the particles is based on the following relations [25]:
V i k + 1 = w V i k + c 1 r a n d 1 ( X l b , i k X i k ) + c 2 r a n d 2 ( X g b k X i k )
X i k + 1 = X i k + V i k + 1
w = w m a x k ( w m a x w m i n ) N
where X is the position vector, V is the velocity vector, w is the inertia weight, I is the particle index, k is the iteration index, Xlb and Xgb are the local and global optimum positions of the swarm, and N is the total number of iterations. The parameters c1 and c2 are the two acceleration constants, and rand1 and rand2 are two random numbers with a uniform distribution in the range [0.0, 1.0]. Typically, the maximum value of c1 + c2 should be 4.0, and an effective initial condition is c1 = c2 = 2 [16,26,27]. In this paper, c1, c2, wmax, and wmin are set to 2, 2, 1, and 0.1, respectively, grounds for experimental results.

3.2. Problem Formulation

The goal of the optimization algorithm is to minimize the objective function J which is given as:
J = 1 | max j = 1 , 2 , ... , 13 { Re ( λ j ) } | + f p f p = { 1000 ,     if  max j = 1 , 2 , ... , 13 { Re ( λ j ) } 0 0 ,       otherwise
where λj is the jth mode eigenvalue of the system matrix A in Equation (18) and fp is the penalty function which is for determining whether the system is stable or not. Equation (22) is referred to in reference [25] with a slight modification of the penalty function. If the system is unstable or marginally stable which means the maximum value of the real part of eigenvalues is larger than or equal to 0, fp becomes larger to prevent unstable or marginally stable system. The minimum and the maximum values for p- and i-gains of the PI controllers are set to 0.01 p.u. and 20 p.u., respectively. The PI controller parameters will be set to the minimum if they are lesser than the minimum and will be set to the maximum if they are larger than the maximum.

4. Case Studies

Three cases are studied in this paper to investigate the effectiveness of the proposed method. In Case I, the PI controller parameters, which are obtained by an experimental simulation test, are used to test the control performance of the PMSG WT. A classical method of tuning the parameters is adopted in Case II. In this case, participation factors of state variables of Case I are acquired, and then the PI controller parameters to be tuned are chosen. With changing the selected parameters, traces of eigenvalues are obtained. By observing the trace of eigenvalues, parameters that place eigenvalue at the leftmost point of the complex plane can be selected. The last case is the proposed method which uses PSO to minimize the objective function J in Equation (22). This optimization algorithm is applied to various wind speed points to change the parameters if the wind speed is changed. At every wind speed, participation factors are obtained to select the PI controller parameters to be controlled.

4.1. Case I—Adjusting Experimentally

In Case I, the PI controller parameters are adjusted heuristically by simulation test. This was done by MATLAB/SimPowerSystems model of the test system. After the test model is established, PI controller parameters are acquired by running simulation with trial and error. Among a number of trial runs of simulation by changing the parameters, a set of them that gives the best simulation results are chosen. They are shown in Table 1. With these parameters, the eigenvalues are obtained at the wind speed of 8 m/s and they are shown in Table 2. As shown in Table 2 for Case I, all of the real parts of the eigenvalues are negative, so the system is stable. However, to improve control performance, it is essential to move the dominant system poles, which are the eigenvalues that have large value of real part relative to the others, to the left side of the complex plane. In this case, λ5 and λ6 are of interest to relocate. The classical method to carry out this procedure is investigated in Case II.
Table 1. PI controller parameters for three cases at vw0 = 8 m/s.
Table 1. PI controller parameters for three cases at vw0 = 8 m/s.
ParametersValue (in p.u.)
Case ICase IIProposed method
Kp1114.14
Ki1114.10
Kp20.10.010.01
Ki20.010.060.14
Kp3110.73
Ki30.20.22.57
Kp4112.91
Ki40.50.53.78
Kp51117.46
Ki50.20.217.15
Kp6114.61
Ki6113.82
Kp7112.12
Ki71.21.20.21
Table 2. Eigenvalues for three cases at vw0 = 8 m/s.
Table 2. Eigenvalues for three cases at vw0 = 8 m/s.
IndexEigenvalues
Case ICase IIProposed method
λ1−720.76−713.70−13945.12
λ2−333.29−341.44−165.88 + j627.30
λ3−13.98 + j80.97−7.67 + j78.26−25.21 + j202.28
λ4−13.98 − j80.97−7.67 − j78.26−25.21 − j202.28
λ5−2.36 + j80.59−7.01 + j6.53−15.01
λ6−2.36 − j80.59−7.01 − j6.53−15.03
λ7−80.56−80.32−379.52
λ8−95.76−98.24−165.88 − j627.30
λ9−261.31 + j357.16−261.31 + j357.16−1673.04
λ10−261.31 − j357.16−261.31 − j357.16−480.39
λ11−737.57 + j503.99−737.57 + j503.99−9524.57
λ12−737.57 − j503.99−737.57 − j503.99−37.22
λ13−176.62−176.62−264.98

4.2. Case II—Classical Method

Since there were 14 PI controller parameters in the PMSG WT system, selecting parameters that are strongly related to the dominant poles is helpful to make the problem simpler. This procedure can be carried out by searching the relationship between the eigenvalues and the state variables from participation factors of those state variables [28,29].
Table 3 lists the participation factors. It can be seen that the eigenvalues to be relocated are dependent on the two state variables ωe and φ2, which are related to the parameters Kp2 and Ki2, as can be seen from Figure 2. To coordinate dominant system poles at the leftmost point of the complex plane, the traces of the dominant system poles Kp2 and Ki2 are plotted, as shown in Figure 6. Here λ5 moves to the left-hand side of the complex plane, whereas λ6 moves to the right-hand side and, eventually, the real part of it becomes larger than 0, which indicates that the system became unstable.
Table 3. Magnitudes of participation factors of state variables for Case I.
Table 3. Magnitudes of participation factors of state variables for Case I.
λ1λ2λ3,4λ5,6λ7λ8λ9,10λ11,12λ13
ωe0000.5100000
imd0000000.6200
imq0.011.380.01000.39000
Vdc0.030.010.5400.030000
igd1.140.020.0500.080000
igq00000001.050.05
φ10000000.6200
φ2000.010.5200000
φ300.4000.011.41000
φ40.010.020.5300.050.01000
φ50.1300.0301.170.02000
φ600000000.270.86
φ700000000.820.09
Figure 6b illustrates that λ5 and λ6 initially move toward the left-hand side of the complex plane, until the real part becomes approximately ‒15; subsequently one of them moves in the opposite direction and the other continues to move in the same direction. The traces of λ3 and λ4 should be noted, since they move to the right-hand side as Ki2 increases. From the traces of eigenvalues, we may conclude that Kp2 should be decreased and Ki2 should be increased appropriately. Tuned values of the parameters and the corresponding eigenvalues for Case II are listed in Table 1 and Table 2, respectively. It can be seen that the real parts of λ5 and λ6 in Case II were smaller than in Case I; however, the real parts of λ3 and λ4 became larger and approached those of λ5 and λ6. In this case, other parameters that those related to λ3 and λ4 should also be adjusted to improve the control system; however, it is not straightforward to search the traces of eigenvalues. To solve this problem, we apply PSO, as described in the following subsection.
Figure 6. Eigenvalue analysis. (a) Trace of the dominant poles when Kp2 increases; (b) Trace of the dominant poles when Ki2 increases.
Figure 6. Eigenvalue analysis. (a) Trace of the dominant poles when Kp2 increases; (b) Trace of the dominant poles when Ki2 increases.
Energies 08 01406 g006

4.3. Case III—Proposed Method

Figure 7 shows a flow chart of the PSO method. At wind speed intervals of 0.1 m/s from the cut-in wind speed of 3 m/s to the rated wind speed 11 m/s, the PI controller parameters that minimize J were searched for. Initially, the wind speed was set to vw = 3 m/s and all of the parameters were adjusted to minimize J. Then vw was increased by 0.1 m/s, with the optimized parameters acquired using PSO at vw = 3 m/s used to initialize the problem. The dominant pole pd was determined, which had the largest real part, and the other poles with real parts that were of similar magnitude to that of pd (with a difference of less than a tolerance ε) were found, if they existed. This tolerance for locating additional significant pole was ε = 0.5. A set of PI controller parameters Sp, which are strongly dependent on the dominant poles, were searched for by using the participation factors. By setting these parameters Sp as variables to be tuned, rather than all 14 parameters, the minimization of J is simplified. The PSO algorithm was then implemented to obtain the parameter set Sp. This process was repeated in 0.1-m/s increments of vw, up to vw = 11 m/s.
Figure 7. The whole flow chart of the proposed method.
Figure 7. The whole flow chart of the proposed method.
Energies 08 01406 g007
The resulting PI controller parameters are shown in Figure 8, where Kp2, Ki2, Kp4, and Ki4 are plotted as a function of the wind speeds in the range 3 ≤ vw ≤ 11 m/s. The other parameters remained the same as those listed in Table 1, which were obtained at vw = 3 m/s. As shown in Table 4, the dominant poles at vw = 3 m/s were λ5 and λ6 and were strongly dependent on the state variables ωe and φ2, which follows from the participation factors. Those variables are related to Kp2 and Ki2, as mentioned previously, which is to be expected since they are the controller parameters that determine the active power, and react much more slowly than those that determine the currents. During the PSO optimization, Kp2 remained constant at the minimum, since it tended to decrease for the whole range of the wind speeds. The parameter Ki2 also tended to decrease as the wind speed increased. Only these two parameters were controlled until the real part value of λ3 and λ4 became closer to those of λ5 and λ6 at vw = 9.9 m/s (see Table 4). Thus at vw = 10 m/s, Kp4 and Ki4 (which are related to λ3 and λ4) were also selected as parameters to be modified. After adjusting Kp4 and Ki4 at vw = 10 m/s, the real parts of λ5 and λ6 became much smaller than those of λ3 and λ4 (see Table 4). Subsequently, only Kp2 and Ki2 were adjusted for wind speeds of vw > 10 m/s.
Figure 8. The parameters obtained by using PSO algorithm for various wind speeds. (a) Kp2 and Ki2 for vw = 3–11 m/s. (b) Kp4 and Ki4 for vw = 3–11 m/s.
Figure 8. The parameters obtained by using PSO algorithm for various wind speeds. (a) Kp2 and Ki2 for vw = 3–11 m/s. (b) Kp4 and Ki4 for vw = 3–11 m/s.
Energies 08 01406 g008
Table 4. Eigenvalues for the proposed method at particular wind speeds.
Table 4. Eigenvalues for the proposed method at particular wind speeds.
IndexEigenvalues
vw = 3 m/svw = 9.9 m/svw = 10 m/s
λ1−13944.43−13946.34−13888.08
λ2−160.40 + j626.99−168.52 + j627.74−173.50 + j623.71
λ3−39.60 + j208.09−18.75 + j197.66−42.21 + j205.74
λ4−39.60 − j208.09−18.75 − j197.66−42.21 − j205.74
λ5−5.68 + j0.002−18.45−18.65 + j0.006
λ6−5.68 − j0.002−18.47−18.65 − j0.006
λ7−379.47−379.55−380.57
λ8−160.40 − j626.99−168.52 − j627.74−173.50 − j623.71
λ9−1673.04−1673.04−1673.04
λ10−480.39−480.39−480.39
λ11−9524.57−9524.57−9524.57
λ12−37.22−37.22−37.22
λ13−264.98−264.98−264.98

5. Simulation Verification

The effectiveness of the proposed method was investigated using simulations. The parameters listed in Table 1 were applied to the PMSG WT modeled using the MATLAB/SimPowerSystems libraries. The parameters describing the PMSG WT system are given in the Appendix. Due to the large number of sets of PI controller parameters for the whole range of wind speeds considered in Case III, we describe a simulation test at vw = 8 m/s, since this can demonstrate the superiority of the PSO method.
Figure 9. Simulation test results for three cases. (a) Wind speed; (b) Active power output; (c) Reactive power output; (d) DC-link voltage; (e) MSC d-axis current; (f) MSC q-axis current; (g) GSC d-axis current; (h) GSC q-axis current.
Figure 9. Simulation test results for three cases. (a) Wind speed; (b) Active power output; (c) Reactive power output; (d) DC-link voltage; (e) MSC d-axis current; (f) MSC q-axis current; (g) GSC d-axis current; (h) GSC q-axis current.
Energies 08 01406 g009
Figure 9a shows the time dependence of the wind speed that was applied to the PMSG WT, and Figure 9b–h show the dynamics of the system state variables and output variables for the three control methods. The data were filtered using a low-pass filter to provide a better comparison, as the original signals contained harmonic distortion. From Figure 9b,f–h, it can be seen that the signals from the PSO method (Case III) stabilized the system to a steady-state more rapidly than those of the other control methods. Moreover, overshoots of the signals were smaller with the PSO method. From Figure 9c,d, we can see improved performance in terms of Qout and Vdc with the PSO method, which maintained their reference values well, whereas those of the other methods temporarily deviated from the reference values.

6. Conclusions

We have described a tuning method for the PI controller parameters of a PMSG WT. To acquire the dominant poles of the system and evaluate the stability, a small-signal analysis of the PMSG WT was carried out. The PI controller parameters were then adjusted to minimize the real part of the dominant poles. This process was implemented using a PSO algorithm, which is an efficient algorithm to solve nonlinear optimization problems.
Since the operating conditions of the PMSG WT change when the wind speed fluctuates, parameter sets were obtained for the wind speeds with 0.1 m/s increments in the range from the cut-in speed to the rated speed (i.e., 3–11 m/s). To simplify the problem, the PI controller parameters that were strongly related to the dominant poles were selected to be adjusted using the PSO algorithm using the participation factors; this was carried out for each wind speed investigated, except vw = 3 m/s, for which all parameters were optimized.
To verify the effectiveness of the proposed method, the PMSG WT shown in Figure 1 was modeled using the MATLAB/SimPowerSystems libraries. Based on these simulated data, we may conclude that the control performance of the proposed method gives faster rising and settling time and has smaller overshoot than those of the other cases.

Acknowledgments

This work was supported by the Human Resources Development program (No. 20124030200030) of the Korea Institute of Energy Technology Evaluation and Planning (KETEP) grant funded by the Korea government Ministry of Trade, Industry and Energy.
This work was supported by the New & Renewable Energy program (No. 2012T100100064) of the Korea Institute of Energy Technology Evaluation and Planning (KETEP) grant funded by the Korea government Ministry of Trade, Industry and Energy.

Author Contributions

Yun-Su Kim modeled and tested simulations and compiled the manuscript. Il-Yop Chung collected the data and investigated early works. Seung-Il Moon supervised this work.

Appendix

A1. Derivation of Matrix A

From Equations (1)–(4), (7) and (8), a derivative of ωe can be given as:
d ω e d t = N p p J e q ( T e T m ) = N p p J e q [ 3 2 N p p i m q { λ p m + ( L d L q ) i m d } P w ω m ]         = N p p 2 2 J e q [ 3 i m q { λ p m + ( L d L q ) i m d } + ρ π r 2 v w 3 ω e × ( 110.23 N p p N g r v w r ω e 9.30531 )                 × exp ( 18.4 N p p N g r v w r ω e 0.0552 ) ]
A linearization of Equation (23) around steady-state values gives:
d Δ ω e d t = a 1 , 1 Δ ω e + a 1 , 2 Δ i m d + a 1 , 3 Δ i m q
From Equations (9) and (11), derivatives of imd and imq can be calculated. The derivative of imd can be given as:
d i m d d t = 1 L d { ( 0 i m d ) K p 1 + φ 1 K i 1 R s i m d }
A linearization of Equation (25) around steady-state values gives:
d Δ i m d d t = a 2 , 2 Δ i m d + a 2 , 7 Δ φ 1
And the derivative of imq can be given as:
d i m q d t = 1 L q [ { ( 3 2 v s d i g d K o p t ω e 3 N p p 3 N g r 3 ) K p 2 + φ 2 K i 2 i m q } K p 3 + φ 3 K i 3 R s i m q ]
In Equation (27), the variable vsd, which is not defined as the state variable, is contained so a linearization of vsd should be expressed as terms of linearization of state variables. From Equation (17), the infinite bus voltage can be given as
| v i d q | cos δ + j | v i d q | sin δ = v s d + i g q x T L j i g d x T L
From Equation (28), square of the magnitude of the infinite bus voltage can be given as:
| v i d q | 2 = ( v s d + i g q x T L ) 2 + ( i g d x T L ) 2
Since |vidq| is a constant value, a linearization of Equation (29) gives
Δ v s d = i g d 0 x T L 2 v s d 0 + i g q 0 x T L Δ i g d x T L Δ i g q
which is expressed in terms of the state variables igd and igq. By using Equation (30), Equation (27) can be linearized as:
d Δ i m q d t = a 3 , 1 Δ ω e + a 3 , 3 Δ i m q + a 3 , 5 Δ i g d + a 3 , 6 Δ i g q + a 3 , 8 Δ φ 2 + a 3 , 9 Δ φ 3
From Equations (11) and (13), a derivative of Vdc can be given as:
d V d c d t = 3 2 C d c V d c ( v m d i m d + v m q i m q + v s d i g d )
Equation (32) contains variables vmd and vmq that are not defined as the state variables so they should be expressed in terms of the state variables when linearizing. From Equations (9) and (11), vmd and vmq can be linearized as:
Δ v m d = K p 1 Δ i m d + K i 1 Δ φ 1 L q i m q 0 Δ ω e L q ω e 0 Δ i m q Δ v m q = K p 3 Δ i m q + ( 3 K o p t ω e 0 2 K p 2 K p 3 N p p 3 N g r 3 + L d i m d 0 + λ p m ) Δ ω e + K i 2 K p 3 Δ φ 2 + K i 3 Δ φ 3 3 K p 2 K p 3 i g d 0 x T L 2 Δ i g q + L d ω e 0 Δ i m d + 3 K p 2 K p 3 2 ( v s d 0 i g d 0 2 x T L 2 v s d 0 + i g q 0 x T L ) Δ i g d
By using Equation (33), Equation (32) can be linearized as:
d Δ V d c d t = a 4 , 1 Δ ω e + a 4 , 2 Δ i m d + a 4 , 3 Δ i m q + a 4 , 4 Δ V d c + a 4 , 5 Δ i g d + a 4 , 6 Δ i g q + a 4 , 7 Δ φ 1 + a 4 , 8 Δ φ 2 + a 4 , 9 Δ φ 3
Derivatives of igd and igq can be obtained from Equation (15). The derivative of igd can be given as:
d i g d d t = 1 L f [ { ( V d c 5400 ) K p 4 + φ 4 K i 4 i g d } K p 5 + φ 5 K i 5 R f i g d ]
Subsequently, Equation (35) can be linearized as:
d Δ i g d d t = a 5 , 4 Δ V d c + a 5 , 5 Δ i g d + a 5 , 10 Δ φ 4 + a 5 , 11 Δ φ 5
And the derivative of igq can be given as:
d i g q d t = 1 L f [ { ( 3 2 v s d i g q 0 ) K p 6 + φ 6 K i 6 i g q } K p 7 + φ 7 K i 7 R f i g q ]
Then a linearization of Equation (37) can be given as
d Δ i g q d t = a 6 , 5 Δ i g d + a 6 , 6 Δ i g q + a 6 , 12 Δ φ 6 + a 6 , 13 Δ φ 7
Derivatives of φ1φ3 are given in Equation (11) and by using Equations (11) and (30), their linearization can be given as
d Δ φ 1 d t = a 7 , 2 Δ i m d
d Δ φ 2 d t = a 8 , 1 Δ ω e + a 8 , 5 Δ i g d + a 8 , 6 Δ i g q
d Δ φ 3 d t = a 9 , 1 Δ ω e + a 9 , 3 Δ i m q + a 9 , 5 Δ i g d + a 9 , 6 Δ i g q + a 9 , 8 Δ φ 2
Derivatives of φ4φ7 are given in Equation (15) and by using Equations (15) and (30), their linearization can be given as:
d Δ φ 4 d t = a 10 , 4 Δ V d c
d Δ φ 5 d t = a 11 , 4 Δ V d c + a 11 , 5 Δ i g d + a 11 , 10 Δ φ 4
d Δ φ 6 d t = a 12 , 5 Δ i g d + a 12 , 6 Δ i g q
d Δ φ 7 d t = a 13 , 5 Δ i g d + a 13 , 6 Δ i g q + a 13 , 12 Δ φ 6

A2. Parameters of the PMSG WT System

ρ = 1.225 (kg/m3); r = 83.5 (m); Cdc = 8 (mF); Cp,max = 0.4412; λopt = 7.2064; Ngr = 30; Npp = 9; Rs = 8.67 (mΩ); Ld = 2.86 (mH); Lq = 3.44 (mH); λpm = 7.15 (V∙s); Lf = 1.8 (mH); Rf = 0.65 (mΩ); xT = 0.2371 (Ω); xL = 0.0013 (Ω).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Thresher, R.; Robinson, M.; Veers, P. To capture the wind. IEEE Power Energy Mag. 2007, 5, 34–46. [Google Scholar] [CrossRef]
  2. Goudarzi, N.; Zhu, W.D. A review on the development of wind turbine generators across the world. Int. J. Dynam. Control. 2013, 1, 192–202. [Google Scholar] [CrossRef]
  3. Global Wind Energy Council. Global Wind Energy Outlook 2012; Global Wind Energy Council: Brussels, Belgium, 2012. Available online: http://www.gwec.net/wp-content/uploads/2012/11/GWEO_2012_lowRes.pdf (accessed on 5 June 2013).
  4. Kim, Y.-S.; Chung, I.-Y.; Moon, S.-I. An analysis of variable-speed wind turbine power-control methods with fluctuating wind speed. Energies 2013, 6, 3323–3338. [Google Scholar] [CrossRef]
  5. Iglesias, R.L.; Arantegui, R.L.; Alonso, M.A. Power electronics evolution in wind turbines—A market-based analysis. Renew. Sustain. Energy Rev. 2011, 15, 4982–4993. [Google Scholar] [CrossRef]
  6. Benelghali, S.; Benbouzid, M.E.H.; Charpentier, J.F. Comparison of PMSG and DFIG for Marine Current Turbine Applications. In Electrical Machines (ICEM); IEEE: Rome, Italy, 2010; pp. 1–6. [Google Scholar]
  7. Knospe, C. PID control. IEEE Control Systems Mag. 2006, 26, 30–31. [Google Scholar] [CrossRef]
  8. Åström, K.J.; Panagopoulos, H.; Hägglund, T. Design of PI controllers based on non-convex optimization. Automatica 1998, 34, 585–601. [Google Scholar] [CrossRef]
  9. Wu, F.; Zhang, X.-P.; Godfrey, K.; Ju, P. Small signal stability analysis and optimal control of a wind turbine with doubly fed induction generator. IET Gener. Transm. Distrib. 2007, 1, 751–760. [Google Scholar] [CrossRef]
  10. Huang, H.; Mao, C.; Lu, J.; Wang, D. Small-signal modelling and analysis of wind turbine with direct driven permanent magnet synchronous generator connected to power grid. IET Renew. Power Gener. 2012, 6, 48–58. [Google Scholar] [CrossRef]
  11. Hasanien, H.M.; Muyeen, S.M. Design optimization of controller parameters used in variable speed wind energy conversion system by genetic algorithms. IEEE Trans. Sustain. Energy 2012, 2, 200–208. [Google Scholar] [CrossRef]
  12. Geng, H.; Xu, D. Stability analysis and improvements for variable-speed multipole permanent magnet synchronous generator-based wind energy conversion system. IEEE Trans. Sustain. Energy 2011, 2, 459–467. [Google Scholar] [CrossRef]
  13. Wu, F.; Zhang, X.-P.; Ju, P. Small signal stability analysis and control of the wind turbine with the direct-drive permanent magnet generator integrated to the grid. Electr. Power Syst. Res. 2009, 79, 1661–1667. [Google Scholar] [CrossRef]
  14. Mishra, Y.; Mishra, S.; Li, F.; Dong, Z.Y.; Bansal, R.C. Small-signal stability analysis of a DFIG-based wind power system under different modes of operation. IEEE Trans. Energy Convs. 2009, 24, 972–982. [Google Scholar] [CrossRef]
  15. Yang, L.; Yang, G.Y.; Xu, Z.; Dong, Z.Y.; Wong, K.P.; Ma, X. Optimal controller design of a doubly-fed induction generator wind turbine system for small signal stability enhancement. IET Gener. Transm. Distrib. 2010, 4, 579–597. [Google Scholar] [CrossRef]
  16. Valle, Y.; Venayagamoorthy, G.K.; Mohagheghi, S.; Hernandez, J.-C.; Harley, R.G. Particle swarm optimization: Basic concepts, variants and applications in power systems. IEEE Trans. Evol. Comput. 2008, 12, 171–195. [Google Scholar] [CrossRef]
  17. Ackermann, T. Wind Power in Power Systems; John Wiley & Sons: Hoboken, NJ, USA, 2005. [Google Scholar]
  18. Lanzafame, R.; Messina, M. Horizontal axis wind turbine working at maximum power coefficient continuously. Renew. Energy 2010, 35, 301–306. [Google Scholar] [CrossRef]
  19. Abdullah, M.A.; Yatim, A.H.M.; Tan, C.W.; Saidur, R. A review of maximum power point tracking algorithms for wind energy systems. Renew. Sustain. Energy Rev. 2012, 16, 3220–3227. [Google Scholar] [CrossRef]
  20. Yazdani, A; Iravani, R. Voltage-Sourced Converters in Power Systems; John Wiley & Sons: Hoboken, NJ, USA, 2010. [Google Scholar]
  21. Novotny, D.W.; Lipo, T.A. Vector Control and Dynamics of AC Drives; Oxford University Press: Oxford, UK, 1996. [Google Scholar]
  22. Leonhard, W. Control of Electrical Drives; Springer-Verlag: New York, NY, USA, 2001. [Google Scholar]
  23. Yazdani, A.; Iravani, R. A neutral-point clamped converter system for direct-drive variable-speed wind power unit. IEEE Trans. Energy Convs. 2006, 21, 596–607. [Google Scholar] [CrossRef]
  24. Nise, N.S. Control Systems Engineering; John Wiley & Sons: Hoboken, NJ, USA, 2008. [Google Scholar]
  25. Chung, I.-Y.; Liu, W.; Cartes, D.A.; Collins, E.G.; Moon, S.-I. Control methods of inverter-interfaced distributed generators in a microgrid system. IEEE Trans. Ind. Appl. 2010, 46, 1078–1088. [Google Scholar] [CrossRef]
  26. Eberhart, R.; Shi, Y.; Kennedy, J. Swarm Intelligence; Morgan Kaufmann: San Francisco, CA, USA, 2001. [Google Scholar]
  27. Ozcan, E.; Mohan, C. Particle swarm optimization: Surfing the waves. Proc. IEEE Congr. Evol. Comput. 1999, 3, 1939–1944. [Google Scholar]
  28. Verghese, G.C.; Perez-Arriaga, I.J.; Schweppe, F.C. Selective modal analysis with application to electric power systems, Part I: Heuristic introduction. IEEE Trans. Power App. Syst. 1982, PAS-101, 3117–3125. [Google Scholar]
  29. Kundur, P. Power System Stability and Control; McGraw-Hill: New York, NY, USA, 1994. [Google Scholar]

Share and Cite

MDPI and ACS Style

Kim, Y.-S.; Chung, I.-Y.; Moon, S.-I. Tuning of the PI Controller Parameters of a PMSG Wind Turbine to Improve Control Performance under Various Wind Speeds. Energies 2015, 8, 1406-1425. https://doi.org/10.3390/en8021406

AMA Style

Kim Y-S, Chung I-Y, Moon S-I. Tuning of the PI Controller Parameters of a PMSG Wind Turbine to Improve Control Performance under Various Wind Speeds. Energies. 2015; 8(2):1406-1425. https://doi.org/10.3390/en8021406

Chicago/Turabian Style

Kim, Yun-Su, Il-Yop Chung, and Seung-Il Moon. 2015. "Tuning of the PI Controller Parameters of a PMSG Wind Turbine to Improve Control Performance under Various Wind Speeds" Energies 8, no. 2: 1406-1425. https://doi.org/10.3390/en8021406

Article Metrics

Back to TopTop