Next Article in Journal
Energy Flexibility Management Based on Predictive Dispatch Model of Domestic Energy Management System
Next Article in Special Issue
Techno-Economic Analysis of Solar Tower Aided Coal-Fired Power Generation System
Previous Article in Journal
Thermal Impact Assessment of Groundwater Heat Pumps (GWHPs): Rigorous vs. Simplified Models
Previous Article in Special Issue
Improved Reactive Current Detection Method of SVG
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Proportional-Type Performance Recovery DC-Link Voltage Tracking Algorithm for Permanent Magnet Synchronous Generators

Department of Creative Convergence Engineering, Hanbat National University, 125, Dongseo-daero, Yuseong-gu, Daejeon 34158, Korea
Energies 2017, 10(9), 1387; https://doi.org/10.3390/en10091387
Submission received: 14 August 2017 / Revised: 5 September 2017 / Accepted: 7 September 2017 / Published: 12 September 2017

Abstract

:
This study proposes a disturbance observer-based proportional-type DC-link voltage tracking algorithm for permanent magnet synchronous generators (PMSGs). The proposed technique feedbacks the only proportional term of the tracking errors, and it contains the nominal static and dynamic feed-forward compensators coming from the first-order disturbance observers. It is rigorously proved that the proposed method ensures the performance recovery and offset-free properties without the use of the integrators of the tracking errors. A wind power generation system has been simulated to verify the efficacy of the proposed method using the PSIM (PowerSIM) software with the DLL (Dynamic Link Library) block.

1. Introduction

Because of the high efficiency, high power, and simple structures, the use of permanent magnet synchronous machines (PMSMs) for various motoring and generating applications has been preferred. Eliminating the external rotor excitation removes the rotor losses, which improves the PMSM efficiency dramatically. Moreover, the absence of rotor winding eliminates not only the slip rings but also the brushes, which leads to a considerable reduction of the maintenance costs [1,2,3,4,5,6].
In power generating applications, the PMSM plays the role of a balanced three-phase AC power supply depending on the input mechanical power source, such as wind power, thermoelectric power, nuclear power, and so on. The DC-link voltage across the output capacitor should be controlled by the three-phase inverter with properly designed control algorithms. Because the corresponding DC-link voltage control problem of the PMSM-based power system is equivalent to the case of the AC/DC converter with a variable AC power source, the extant solutions for the AC/DC converter output voltage control problems can be applied. The cascade control strategy has mainly been adopted for controlling the output voltage of AC/DC converters, where the outer-loop controller produces a desired d-axis current reference for the inner-loop current controller so as to regulate the output voltage. Because of the simple structure, the proportional-integral (PI) controller has commonly been used for the both inner and outer loops [7]. In the AC/DC converter model, there exist inherent nonlinearities and disturbance terms in the output voltage and current dynamics, which could cause a severe closed-loop degradation or instability. Several novel control techniques have been applied for developing the inner-loop current controller to improve the output voltage tracking performance by dealing with the nonlinearities and disturbance terms accordingly. For instance, feedback-linearizing [8,9], passivity-based [10], and predictive [11] schemes remarkably enhanced the closed-loop performance by canceling the nonlinearity and disturbances using the true parameters of the AC/DC converter. The model predictive controllers (MPCs) in [12] were designed to try to enhance the transient performance through minimizing the cost function tracking errors, by using an exhaustive search method for each sampling time. These advanced methods, however, require the use of the true parameters of the converter, such as the resistance, inductance, and capacitance, to ensure closed-loop stability with the desired closed-loop performance. Thus, it is questionable that these advanced schemes can still provide the desired closed-loop performance with the stability guarantee, as the converter parameters and the load current can be dramatically changed in accordance with the operating conditions. An adaptive current controller in [13] has been devised in order to deal with the parameter uncertainties, but it is ambiguous whether the current convergence property is still valid without an offset error in the real implementation due to the absence of the integrators in the tracking errors. The recent sliding mode controller in [14] was developed in order to enhance the output voltage tracking performance in the transient periods by suppressing the disturbance terms caused by a plant-model mismatch. In the case of the proportional type output voltage controller in [15], there were no stability analysis and steady-state analysis with respect to the offset errors. Disturbance observer (DOB)-based output voltage controllers are used in [16,17] that include the integrators of the tracking errors, which may require the use of anti-windup algorithms.
There have been many multi-variable approaches in [18,19,20,21] to solve the output voltage tracking problem of the AC/DC converter. These methods have been made through two steps; first, a positive definite function with respect to the tracking errors was defined, and second, a control action was derived to make the positive definite function monotonically decreasing for all time. However, knowledge of the true parameters is required to ensure the closed-loop performances.
This paper offers a DOB-based proportional-type DC-link voltage tracking algorithm, considering the nonlinearity in the DC-link voltage and PMSM dynamics with the model-plant mismatches. The proposed method is derived through a multi-variable approach, which combines the simple proportional-type feedback linearizing (FL) controller with the first-order DOBs. The contribution of this article is summarized as follows. First, it incorporates the first-order DOBs in the proportional-type FL DC-link voltage control algorithm so as to stabilize the tracking errors. Second, it is rigorously proved that the closed-loop system ensures the performance recovery property as well as the offset-free property without the use of the integrators of the tracking errors. The realistic simulation results confirm the efficacy of the proposed method, in which the wind power system with the output capacitor has been emulated by PSIM (PowerSIM) software.

2. Permanent Magnet Synchronous Generator Model in Rotational d-q Axis

In the rotational d-q axis, the electrical and mechanical dynamics of the permanent magnet synchronous generator (PMSG) are given by [22]:
L d d i d ( t ) d t = R s i d ( t ) + L q ω r ( t ) i q ( t ) + u d ( t )
L q d i q ( t ) d t = R s i q ( t ) L d ω r ( t ) i d ( t ) λ P M ω r ( t ) + u q ( t )
J d ω ( t ) d t = B ω ( t ) + T m ( t ) T e ( i d ( t ) , i q ( t ) ) ,
t 0 , where i d ( t ) , i q ( t ) , and ω r ( t ) : = P ω ( t ) denote the d-q axis’ current and electrical speed, respectively, and ω ( t ) and P represent the mechanical speed and the number of pole pairs, respectively. The PMSG parameters R s , L d , L q , λ P M , B, and J represent the stator resistance, d-q inductances, magnet flux, viscous friction, and rotor moment of inertia, respectively, and T m ( t ) is the input mechanical torque. The input voltages u d ( t ) and u q ( t ) are treated as the control input, and the electrical torque T e ( i d ( t ) , i q ( t ) ) is given by
T e ( i d ( t ) , i q ( t ) ) : = 3 2 P Δ L d q i d ( t ) i q ( t ) + λ P M i q ( t ) , t 0
where Δ L d q : = L d L q .
Letting P c a p ( t ) , P i n ( t ) , and P g r i d ( t ) be the capacitor power, input power, and grid power, respectively, as shown in Figure 1, it holds that
P c a p ( t ) = C v d c ( t ) d v d c ( t ) d t = P i n ( t ) P g r i d ( t ) , t 0
and, because the input power of P i n ( t ) can be decomposed as P i n ( t ) = P G ( t ) + P l o s s , i n v ( t ) , t 0 , where P G ( t ) and P l o s s , i n v ( t ) are the generator output power and the inverter power loss, respectively, it follows from Equation (5) that
C d v d c ( t ) d t = ω ( t ) v d c ( t ) T e ( t ) + P l o s s , i n v ( t ) v d c ( t ) I g r i d ( t ) , t 0
where I g r i d ( t ) denotes the load current to the grid.
It should be noted that it is difficult to consistently identify the true values of the PMSG parameters, with the exception of the number of pole pairs P as well as the output capacitance value, because the parameters can dramatically vary with the operating conditions, such as the phase current, DC-link voltage, and PMSG temperature. Thus, it is reasonable for the dynamical Equations (1), (2), and (6) to be rewritten as
L d , 0 d i d ( t ) d t = R s , 0 i d ( t ) + L q , 0 ω r ( t ) i q ( t ) + u d ( t ) + d d , o ( t )
L q , 0 d i q ( t ) d t = R s , 0 i q ( t ) L d , 0 ω r ( t ) i d ( t ) λ P M , 0 ω r ( t ) + u q ( t ) + d q , o ( t )
C 0 d v d c ( t ) d t = ω ( t ) v d c ( t ) T e , 0 ( t ) + d v , o ( t ) , t 0
with the nominal parameters of R s , 0 , L d , 0 , L q , 0 , λ P M , 0 , and C 0 , where the nominal torque T e , 0 ( t ) is defined as
T e , 0 ( t ) : = T e ( t ) | L d = L d , 0 , L q = L q , 0 , λ P M = λ P M , 0 , t 0
and d d , o ( t ) , d q , o ( t ) , and d v , o ( t ) denote the unknown disturbances caused from the parameter mismatch, unmodeled dynamics, and load uncertainties.
The next section develops a DC-link voltage tracking algorithm based on the perturbed PMSG current dynamics of Equations (7) and (8) with the DC-link voltage dynamics of Equation (9). To this end, the following assumptions are made:
  • The disturbances d d , o ( t ) , d q , o ( t ) , and d v , o ( t ) are unknown, but their time-derivatives are bounded for all time.
  • The stationary a-b-c phase current and the PMSG speed ω ( t ) are measurable so that the rotational d-q axis currents are available for feedback.

3. DC-Link Voltage Tracking Algorithm Design

The DC-link voltage control objective is to force the closed-loop system to behave according to the transfer function given by
V d c * ( s ) V d c , r e f ( s ) = ω v c s + ω v c , s C
where ω v c is the desired cut-off frequency as a design parameter, and V d c * ( s ) and V d c , r e f ( s ) denote the Laplace transforms of the DC-link voltage v d c * ( t ) and its reference v d c , r e f ( t ) , respectively.
The time domain expression of the target dynamics of Equation (10) can be written as
v ˙ d c * ( t ) = ω v c v d c , r e f ( t ) v d c * ( t ) , t 0
and, then, with the target dynamics of Equation (11), defining the error as
v ˜ d c * ( t ) : = v d c * ( t ) v d c ( t ) , t 0
the DC-link voltage tracking error dynamics can be obtained as
C 0 v ˜ ˙ d c * ( t ) = C 0 v ˙ d c * ( t ) C 0 v ˙ d c ( t ) = ω ( t ) v d c ( t ) b i q ( t ) ω ( t ) v d c ( t ) 3 2 P Δ L d q , 0 i d ( t ) i q ( t ) + d v ( t ) = ω ( t ) v d c ( t ) b i q , r e f ( t ) + ω ( t ) v d c ( t ) b i ˜ q ( t ) ω ( t ) v d c ( t ) 3 2 P Δ L d q , 0 i d ( t ) i q ( t ) + d v ( t ) , t 0
where b : = 3 2 P λ P M , 0 , i q , r e f ( t ) denotes the q-axis current reference, the q-axis current tracking error i ˜ q ( t ) is defined as i ˜ q ( t ) : = i q , r e f ( t ) i q ( t ) , t 0 , L d q , 0 : = L d q | L d , L q = L d , 0 , L q , 0 , and d v ( t ) : = C 0 v ˙ d c * ( t ) d v , o ( t ) , t 0 .
The stabilization of the tracking error dynamics of Equation (12) can be established by the proposed q-axis current reference:
i q , r e f ( t ) : = v d c ( t ) b ω ( t ) C 0 λ v c v ˜ d c * ( t ) ω ( t ) v d c ( t ) 3 2 P Δ L d q , 0 i d ( t ) i q ( t ) + d ^ v ( t ) , t 0
where λ v c > 0 refers to the design parameter for adjusting the decay ratio of the DC-link voltage tracking error of v ˜ d c * ( t ) , and d ^ v ( t ) represents the estimated disturbance defined as
d ^ v ( t ) : = ζ v ( t ) + l v C 0 v ˜ d c * ( t ) , t 0
The state ζ v ( t ) satisfies the nonlinear observer dynamics:
ζ ˙ v ( t ) = l v ζ v ( t ) l v 2 C 0 v ˜ d c * ( t ) + l v ω ( t ) v d c ( t ) b i q ( t ) + 3 2 P Δ L d q , 0 i d ( t ) i q ( t ) , t 0
where l v > 0 is the observer gain as a design parameter. For the rest of this article, the nonlinear observer of Equation (15) with the output of Equation (14) is called the DC-link voltage DOB.
Substitution of the proposed q-axis current reference of Equation (13) to the DC-link voltage error dynamics of Equation (12) yields that
v ˜ ˙ d c * ( t ) = λ v c v ˜ d c * ( t ) + ω ( t ) C 0 v d c ( t ) b i ˜ q ( t ) + 1 C 0 d ˜ v ( t ) , t 0
where d ˜ v ( t ) : = d v ( t ) d ^ v ( t ) , t 0 . Lemma 1 presents the stability property of the closed-loop error dynamics of Equation (16), which helps the whole closed-loop stability analysis to be performed considerably easily.
Lemma 1.
For any specified designed parameters λ v c > 0 and l v > 0 , the proposed q-axis current reference of Equation (13) forces the DC-link voltage error dynamics of Equation (16) to be strictly passive with respect to the input–output mapping given by
ω ( t ) C 0 v d c ( t ) b i ˜ q ( t ) γ v d ˙ v ( t ) v ˜ d c * ( t ) d ˜ v ( t )
for some positive constant γ v .
The result of Lemma 1 implies that the DC-link voltage of v d c ( t ) exponentially converges to the target trajectory of v d c * ( t ) if the two signals of i ˜ q ( t ) and d ˙ v ( t ) exponentially vanish. The Appendix A proves Lemma 1 by showing the positive definite function, defined as
V v d c ( t ) : = 1 2 ( v ˜ d c * ( t ) ) 2 + γ v 2 d ˜ v 2 ( t ) , t 0
to be
V ˙ v d c ( t ) α v d c V v d c ( t ) + v ˜ d c * ( t ) ω ( t ) C 0 v d c ( t ) b i ˜ q ( t ) + γ v d ˜ v ( t ) d ˙ v ( t ) , t 0
for some constants α v c > 0 and γ v > 0 . For details, see the Appendix A.
The next step is to derive the d-q axis current error dynamics, which can be obtained as
L d , 0 i ˜ ˙ d ( t ) = L d , 0 i ˙ d , r e f ( t ) L d , 0 i ˙ d ( t ) = R s , 0 i d ( t ) L q , 0 ω r ( t ) i q ( t ) u d ( t ) + d d ( t )
L q , 0 i ˜ ˙ q ( t ) = L q , 0 i ˙ q , r e f ( t ) L q , 0 i ˙ q ( t ) = R s , 0 i q ( t ) + L d , 0 ω r ( t ) i d ( t ) + λ P M , 0 ω r ( t ) u q ( t ) + d q ( t )
where
d d ( t ) : = L d , 0 i ˙ d , r e f ( t ) d d , o ( t ) d q ( t ) : = L q , 0 i ˙ q , r e f ( t ) d q , o ( t ) , t 0
For stabilizing the error dynamics of Equations (16), (20) and (21), the control law is proposed as
u d ( t ) = R s , 0 i d ( t ) L q , 0 ω r ( t ) i q ( t ) + L d , 0 λ c c i ˜ d ( t ) + d ^ d ( t )
u q ( t ) = R s , 0 i q ( t ) + L d , 0 ω r ( t ) i d ( t ) + λ P M , 0 ω r ( t ) + L q , 0 ω ( t ) C 0 v d c ( t ) b v ˜ d c * + L q , 0 λ c c i ˜ q ( t ) + d ^ q ( t ) , t 0
where the design parameter λ c c > 0 is introduced to tune the decay ratio of the current tracking errors i ˜ x ( t ) , x = d , q ; the estimated disturbances, d ^ x ( t ) , x = d , q , are defined as
d ^ x ( t ) : = ζ x ( t ) + l x L x , 0 i ˜ x ( t ) , x = d , q , t 0
and ζ x ( t ) , x = d , q denote the states of the nonlinear observers:
ζ ˙ d ( t ) = l d ζ d ( t ) l d 2 L d , 0 i ˜ d ( t ) + l d R s , 0 i d ( t ) + L q , 0 ω r ( t ) i q ( t ) + u d ( t )
ζ ˙ q ( t ) = l q ζ q ( t ) l q 2 L q , 0 i ˜ q ( t ) + l q R s , 0 i q ( t ) L d , 0 ω r ( t ) i d ( t ) λ P M , 0 ω r ( t ) + u q ( t ) , t 0
where l x > 0 , x = d , q are the observer gains as design parameters. For the rest of this article, the nonlinear observers of Equations (25) and (26) with the outputs of Equation (24) are called the current DOBs.
Substituting the control laws of Equations (22) and (23) to the d-q current error dynamics of Equations (20) and (21), it is easy to see that
i ˜ ˙ d ( t ) = λ c c i ˜ d ( t ) + 1 L d , 0 d ˜ d ( t )
i ˜ ˙ q ( t ) = λ c c i ˜ q ( t ) ω ( t ) C 0 v d c ( t ) b v ˜ d c * ( t ) + 1 L q , 0 d ˜ q ( t ) , t 0
where d ˜ x ( t ) : = d x ( t ) d ^ x ( t ) and x = d , q , t 0 . Lemma 2 presents the stability property of the closed-loop error dynamics of Equations (27) and (28), which helps the whole closed-loop stability analysis to be performed considerably easily.
Lemma 2.
For any specified designed parameters λ c c > 0 , l x > 0 and x = d , q , the proposed control law of Equations (22) and (23) forces the d-q axis error dynamics of Equations (27) and (28) to be strictly passive with respect to the input–output mapping given by
0 ω ( t ) C 0 v d c ( t ) b v ˜ d c * ( t ) γ d d ˙ d ( t ) γ q d ˙ q ( t ) i ˜ d ( t ) i ˜ q ( t ) d ˜ d ( t ) d ˜ q ( t )
for some positive constants γ x , x = d , q .
The result of Lemma 2 implies that the d-q axis current tracking errors exponentially vanish if the signals of v ˜ d c * ( t ) and d ˙ x ( t ) , x = d , q also vanish. The Appendix A proves Lemma 2 by showing the positive definite function defined as
V d q ( t ) : = x = d , q 1 2 i ˜ x 2 ( t ) + x = d , q γ x 2 d ˜ x 2 ( t ) , t 0
to be
V ˙ d q ( t ) α d q V d q ( t ) i ˜ q ( t ) ω ( t ) C 0 v d c ( t ) b v ˜ d c * ( t ) + x = d , q γ x d ˜ x ( t ) d ˙ x ( t ) , t 0
for some constants α d q > 0 and γ x > 0 , x = d , q . For details, see the Appendix A.
Now, it is ready to analyze the closed-loop tracking error behaviors by using two inequalities of Equations (19) and (31), which are obtained as the results of Lemmas 1 and 2. See the Appendix A for the proof of Theorem 1.
Theorem 1.
For any given control parameters of λ v c > 0 , λ c c > 0 , l v > 0 and l x > 0 , x = d , q , the proposed control law of Equations (22) and (23), with the DOBs of Equations (14), (15), (24) – (26), makes the closed-loop system exponentially recover the target DC-link voltage tracking performance of Equation (11), that is,
| v ˜ d c * ( t ) | γ c l e α c l 2 t , | d ˜ x | d ̲ x , t 0
for some α c l > 0 , γ c l > 0 , and d ̲ x > 0 .
Note that it is not clear whether the closed-loop system suffers the offset errors or not in the real implementations, as the proposed control algorithms of Equations (22) and (23) with the d-axis current reference of Equation (13) do not include the integrators of the tracking errors. Interestingly, the proposed control law guarantees the offset-free property without the use of the integrators of the tracking errors, which means that the control algorithm can be further simplified by removing the corresponding anti-windup algorithms. This point can be considered a practical advantage of this article. For details, see Theorem 2; the proof is given in the Appendix A.
Theorem 2.
In the steady state, the closed-loop system always removes the DC-link voltage steady errors, that is,
v d c 0 = v d c , r e f 0
where v d c 0 and v d c , r e f 0 denote the steady-state values of v d c ( t ) and v d c , r e f ( t ) , respectively.
Remark 1.
The DC-link voltage and the current DOB dynamics can be written in the first-order low-pass filer form:
d ^ ˙ x ( t ) = l x d x ( t ) d ^ x ( t ) , x = v , d , q , t 0
whose transfer functions are given by
D ^ x ( s ) D x ( s ) = l x s + l x , x = v , d , q , s C
which shows that the DOB gains l x > 0 and x = v , d , q , can be tuned for the transfer functions of Equation (34) to have a desired cut-off frequency in r a d / s .

4. Simulations

This section discusses the simulations to show the efficacy of the proposed method, comparing this with the classical feedback linearizing method in [23]. PSIM software was utilized to perform the simulations in which the control algorithm was implemented using the DLL (Dynamic Link Library) block with C-language. The true parameters of the PMSG and the output capacitor were given by
R s = 0.099 Ω , L d = L q = 4.07 mH , λ P M = 0.3166 Wb , P = 40 , J = 0.02 kg · m 2 , B = 0.000425 Nm / rad / s , C = 2350 μ F
and it was assumed that the nominal parameters were given by
R s , 0 = 0.7 R s , L d , 0 = 1.5 L d , L q , 0 = 1.5 L q , λ P M , 0 = 1.2 λ P M , C 0 = 0.6 C
which were used for constructing the control algorithm. The wind turbine part was emulated by adopting the wind turbine block of the renewable energy package in the PSIM software with the nominal output power of 10 kW, a moment of inertia of 1 × 10 3 kg · m 2 , a base wind speed of 20 m/s, a base rotational speed of 50 rpm, and an initial rotational speed of 10 rpm. It was assumed that the wind came from the Weibull distribution-based wind model [24]. The control and the pulse-width modulation (PWM) period were set to be 0.1 ms. Figure 2 depicts the implemented closed-loop system, for which the load-side converter system with an AC load was assumed to be a single-phase passive load for convenience.
The cut-off frequencies were adjusted as f c c = 200 Hz and f v c = 5 Hz so that
ω c c = 2 π f c c = 1256 rad / s , ω v c = 2 π f v c = 31.4 rad / s ,
which was applied to both of the proposed and feedback linearizing methods. The DOB gains were selected as l v = l d = l q = 314 rad/s, for the cut-off frequency of the corresponding transfer functions of Equation (34). Note that the d-axis current reference of i d , r e f ( t ) was set to zero for the maximum torque per ampere (MTPA) operation of the surface-mounted-type PMSG used in this simulation. The decay ratios of the DC-link voltage tracking error and the d-q axis current tracking errors were tuned to be λ v c = 125.6 and λ c c = ω c c = 1256 .
The FL technique in [23] used for a comparison is described by
u d ( t ) = L d , 0 ω c c i ˜ d ( t ) + R s , 0 ω c c 0 t i ˜ d ( τ ) d τ L q , 0 ω r ( t ) i q ( t )
u q ( t ) = L q , 0 ω c c i ˜ q ( t ) + R s , 0 ω c c 0 t i ˜ q ( τ ) d τ + L d , 0 ω r ( t ) i d ( t ) + λ P M , 0 ω r ( t ) , t 0
Here, the q-axis current reference comes from
i q , r e f ( t ) = v d c ( t ) b ω ( t ) 2 C 0 ω v c v ˜ d c ( t ) + C 0 ω v c 2 0 t v ˜ d c ( τ ) d τ , t 0
where v ˜ d c ( t ) = v d c , r e f ( t ) v d c ( t ) , t 0 . Note that it is easy to verify that the FL method of Equations (37)–(39) approximately yields the current and DC-link voltage-loop transfer functions given by
G x ( s ) = ω c c s + ω c c , x = d , q , G v ( s ) = ω v c s + ω v c , s C
through the pole-zero cancellation, provided that the nominal parameters are equal to the true parameters, that is, R s , 0 = R s , L x , 0 = L x , λ P M , 0 = λ P M and C 0 = C , x = d , q . The cut-off frequencies ω c c and ω v c were selected to be the same as for the proposed method.
The first simulation was conducted to evaluate the DC-link voltage tracking performance with a resistive load of R L = 100 Ω ; the DC-link voltage reference was increased from 300 to 500 V, and it was then decreased to 300 V. The simulation results are depicted in Figure 3, and they imply that the proposed method successfully drives the DC-link voltage to rapidly follow its target dynamics of Equation (11) in the presence of the model-plant mismatches. The behaviors of the estimated disturbances and the wind pattern used in this simulation are presented in Figure 4.
The second simulation compared the DC-link voltage regulation performances by applying the pulse resistive load from R L = 100 Ω to R L = 28.6 Ω at the 300 V operation mode. Figure 5 depicts the comparison results, which imply that the proposed method considerably reduces the over/under shoots caused by the load variations, compared to the classical FL method. The correponding PMSG output power behaviors are presented in Figure 6.
The third simulation was performed to show the closed-loop robustness by investigating the DC-link voltage tracking performance variations for different load conditions, such as R L = 100 , 50 , 30 Ω . The DC-link voltage reference was increased from 300 to 600 V. The corresponding result is given in Figure 7, which indicates that the proposed method effectively suppresses the DC-link voltage tracking performance variations for several load conditions.
In the last simulation, the DC-link voltage tracking behaviors were compared at a resistive load of R L = 100 Ω with the same DC-link voltage reference by increasing the cut-off frequencies, such as 2, 5, 8 Hz. As can be clearly seen from the result presented in Figure 8, the proposed method efficiently assigns the desired DC-link voltage tracking performance to the closed-loop system, unlike the classical FL technique.
As shown in the simulation results, in contrast to the classical FL method, the proposed method successfully maintains a satisfactory closed-loop performance despite the operating mode changes. Hence, the proposed method offers almost the same closed-loop performance for various operating ranges without any additional gain scheduling method, which corresponds to a practical advantage of the proposed technique.

5. Conclusions

In this article, a robust DC-link voltage tracking algorithm was proposed by incorporating first-order DOBs in the proportional-type feedback-linearizing DC-link voltage control algorithm. It is rigorously shown that the proposed method guarantees the performance recovery property with the offset-free property. The efficacy of the proposed method has been confirmed by performing realistic simulations for wind power system applications.

Acknowledgments

This research was supported by the newly created professor research fund of Hanbat National University in 2017.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

This section proves Lemmas 1 and 2, and Theorems 1 and 2 in a sequential manner. First, Lemma 1 is proved as follows.
Proof. 
The nonlinear observer of Equation (15) with the estimated disturbance of Equation (14) can be written as
d ^ ˙ v l v C 0 v ˜ ˙ d c * = l v ( d ^ v l v C 0 v ˜ d c * ) l v 2 C 0 v ˜ d c * + l v ω v d c ( b i q + 3 2 P Δ L d q , 0 i d i q ) , t 0
which can be rearranged as
d ^ ˙ v = l v C 0 v ˜ ˙ d c * + ω v d c ( b i q + 3 2 P Δ L d q , 0 i d i q ) d ^ v = l v d ˜ v , t 0
where d ˜ v : = d v d ^ v , t 0 , because it follows from the error dynamics of Equation (12) that
d v = C 0 v ˜ ˙ d c * + ω v d c ( b i q + 3 2 P Δ L d q , 0 i d i q ) , t 0
Now, consider the positive definite function defined in Equation (30):
V v d c = 1 2 ( v ˜ d c * ) 2 + γ v 2 d ˜ v 2 , t 0
where the positive constant γ v will be determined later. The corresponding time derivative of the function of V x , x = d , q , along the trajectories of Equations (A1) and (16), is obtained as
V ˙ v d c = v ˜ d c * v ˜ ˙ d c * + γ v d ˜ v d ˜ ˙ v = v ˜ d c * λ v c v ˜ d c * + ω C 0 v d c b i ˜ q + 1 C 0 d ˜ v + γ v d ˜ v l v d ˜ v + d ˙ v λ v c 2 ( v ˜ d c * ) 2 γ v l v 1 2 λ v c C 0 2 d ˜ v 2 + v ˜ d c * ω C 0 v d c b i ˜ q + γ v d ˜ v d ˙ v , t 0
where the last inequality is verified by applying Young’s inequality:
x y ϵ 2 x 2 + 1 2 ϵ y 2 , x , y R , ϵ > 0
It is easy to see that the positive constant γ v , determined as
γ v : = 1 l v 1 2 λ v c C 0 2 + 1 2
makes the inequality of Equation (A2) as follows:
V ˙ v d c λ v c 2 ( v ˜ d c * ) 2 1 2 d ˜ v 2 + v ˜ d c * ω C 0 v d c b i ˜ q + γ v d ˜ v d ˙ v , α v d c V v d c + v ˜ d c * ω C 0 v d c b i ˜ q + γ v d ˜ v d ˙ v , t 0
where α v d c : = min { λ v c , 1 γ v } , which implies the strict passivity of the input–output mapping of Equation (17).  ☐
The proof of Lemma 2 is given as follows.
Proof. 
The nonlinear observers of Equations (25) and (26) with the estimated disturbances of Equation (24) can be written as
d ^ ˙ d l d L d , 0 i ˜ ˙ d = l d ( d ^ d l d L d , 0 i ˜ d ) l d 2 L d , 0 i ˜ d + l d ( R s , 0 i d + L q , 0 ω r i q + u d ) , d ^ ˙ q l q L q , 0 i ˜ ˙ q = l q ( d ^ q l q L q , 0 i ˜ q ) l q 2 L q , 0 i ˜ q + l q ( R s , 0 i q L d , 0 ω r i d λ P M , 0 ω r + u q ) , t 0
which can be rearranged as
d ^ ˙ d = l d ( L d , 0 i ˜ ˙ d R s , 0 i d + L q , 0 ω r i q + u d d ^ d ) = l d d ˜ d ,
d ^ ˙ q = l q ( L q , 0 i ˜ ˙ q R s , 0 i q L d , 0 ω r i d λ P M , 0 ω r + u q d ^ q ) = l q d ˜ q , t 0
because it follows from the error dynamics of Equations (20) and (21) that
d d = L d , 0 i ˜ ˙ d R s , 0 i d + L q , 0 ω r i q + u d , d q = L q , 0 i ˜ ˙ q R s , 0 i q L d , 0 ω r i d λ P M , 0 ω r + u q , t 0
Now, consider the positive definite function defined in Equation (30):
V d q = x = d , q 1 2 i ˜ x 2 + x = d , q γ x 2 d ˜ x 2 , t 0
where the positive constant γ x , x = d , q is determined later. The corresponding time derivative of the function of V d q , along the trajectories of Equations (A6), (A7), (27) and (28), is obtained as
V ˙ d q = x = d , q i ˜ x i ˜ ˙ x + x = d , q γ x d ˜ x d ˜ ˙ x = i ˜ d λ c c i ˜ d + 1 L d , 0 d ˜ d + i ˜ q λ c c i ˜ q ω C 0 v d c b v ˜ d c * + 1 L q , 0 d ˜ q + x = d , q γ x d ˜ x l x d ˜ x + d ˙ x x = d , q λ c c 2 i ˜ x 2 x = d , q γ x l x 1 2 λ c c L x , 0 2 d ˜ x 2 i ˜ q ω C 0 v d c b v ˜ d c * + x = d , q γ x d ˜ x d ˙ x , t 0
where the last inequality is verified by applying Young’s inequality of Equation (A3). It is easy to see that the positive constants γ x and x = d , q , determined as
γ x : = 1 l x 1 2 λ c c L x , 0 2 + 1 2 , x = d , q
make the inequality of Equation (A8) as follows:
V ˙ d q x = d , q λ c c 2 i ˜ x 2 x = d , q 1 2 d ˜ x 2 i ˜ q ω C 0 v d c b v ˜ d c * + x = d , q γ x d ˜ x d ˙ x α d q V d q i ˜ q ω C 0 v d c b v ˜ d c * + x = d , q γ x d ˜ x d ˙ x , t 0
where α d q : = min { λ c c , 1 γ d , 1 γ q } , which implies the strict passivity of the input–output mapping of Equation (29).  ☐
Theorem 1 is proved as follows.
Proof. 
Define the composite-type positive definite function as
V c l : = V v d c + V d q , t 0
where the two positive definite functions of V x , x = v d c , d q are defined in Equations (18) and (30). Then, the time-derivative of V along the trajectories of Equations (16), (27) and (28) is obtained as
V ˙ c l = V ˙ v d c + V ˙ d q α v d c V v d c + v ˜ d c * ω C 0 v d c b i ˜ q + γ v d ˜ v d ˙ v α d q V d q i ˜ q ω C 0 v d c b v ˜ d c * + x = d , q γ x d ˜ x d ˙ x = α v d c V v d c α d q V d q + x = v , d , q γ x d ˜ x d ˙ x , t 0
where the first inequality is obtained by the two inequalities of Equations (19) and (31). The application of Young’s inequality of Equation (A3) makes the inequality of Equation (A12) as follows:
V ˙ c l α v d c 2 V v d c α d q 2 V d q x = v , d , q c x d ˜ x 2 + x = v , d , q γ x d ˜ x d ˙ x α c l V c l , | d ˜ x | γ x d ˙ ¯ x c x , x = v , d , q , t 0
where α c l : = min { α v d c 2 , α d q 2 } , c v : = α v d c γ v 4 , c d : = α d q γ d 4 , c q : = α d q γ q 4 , and d ˙ ¯ x stands for the positive constant satisfying | d ˙ x | d ˙ ¯ x , x = v , d , q , which verifies that the inequality of Equation (32) holds by the comparison principle in [25].  ☐
The proof of Theorem 2 is given as follows.
Proof. 
As can be seen in the proof of Theorem 1 in the Appendix A, the DOB dynamics can be written as
d ^ ˙ x ( t ) = l x d ˜ x ( t ) , x = v , d , q , t 0
which gives the steady-state relationship:
d ˜ x 0 = 0 , l x > 0 , x = v , d , q
where ( ( · ) ) 0 refers to the steady-state value of ( · ) . On the other hand, the closed-loop error dynamics of Equations (16), (27) and (28) with the target dynamics of Equation (11) yield the steady-state equations:
0 = J R x ˜
0 = λ c c i ˜ d 0
0 = ω v c v d c , r e f 0 ( v d c * ) 0
where
x ˜ = ( v ˜ d c * ) 0 i ˜ q 0 , J = ω 0 b C 0 v d c 0 0 1 1 0 , R = λ v c 0 0 λ c c
Because the matrix ( J R ) is negative definite, that is,
x ˜ T ( J R ) x ˜ = x ˜ T J x ˜ x ˜ T R x ˜ = x ˜ T J T x ˜ x ˜ T R x ˜ = x ˜ T R x ˜ < 0 , x ˜ 0
Equation (A15) implies that
x ˜ = 0
which, together with Equation (A17), concludes that
( v ˜ d c * ) 0 = 0 v d c 0 = ( v d c * ) 0 = v d c , r e f 0
  ☐

References

  1. Chu, L.; Jia, Y.F.; Chen, D.-S.; Xu, N.; Wang, Y.-W.; Tang, X.; Xu, Z. Research on Control Strategies of an Open-End Winding Permanent Magnet Synchronous Driving Motor (OW-PMSM)-Equipped Dual Inverter with a Switchable Winding Mode for Electric Vehicles. Energies 2017, 10, 616. [Google Scholar] [CrossRef]
  2. Ding, X.; Du, M.; Cheng, J.; Chen, F.; Ren, S.; Guo, H. Impact of Silicon Carbide Devices on the Dynamic Performance of Permanent Magnet Synchronous Motor Drive Systems for Electric Vehicles. Energies 2017, 10, 364. [Google Scholar] [CrossRef]
  3. Liu, X.; Du, J.; Liang, D. Analysis and Speed Ripple Mitigation of a Space Vector Pulse Width Modulation-Based Permanent Magnet Synchronous Motor with a Particle Swarm Optimization Algorithm. Energies 2016, 9, 923. [Google Scholar] [CrossRef]
  4. Yan, H.; Xu, Y.; Zou, J. A Phase Current Reconstruction Approach for Three-Phase Permanent-Magnet Synchronous Motor Drive. Energies 2016, 9, 853. [Google Scholar] [CrossRef]
  5. Xu, B.; Mu, F.; Shi, G.; Ji, W.; Zhu, H. State Estimation of Permanent Magnet Synchronous Motor Using Improved Square Root UKF. Energies 2016, 9, 489. [Google Scholar] [CrossRef]
  6. Zhao, J.; Gao, X.; Li, B.; Liu, X.; Guan, X. Open-Phase Fault Tolerance Techniques of Five-Phase Dual-Rotor Permanent Magnet Synchronous Motor. Energies 2015, 8, 12810–12838. [Google Scholar] [CrossRef]
  7. Dixon, J.W.; Ooi, B.T. Indirect current control of a unity power factor sinusoidal current boost type three phase rectifier. IEEE Trans. Ind. Electron. 1988, 35, 508–515. [Google Scholar] [CrossRef]
  8. Lee, T.S. Input-output linearization and zero-dynamics control of three-phase AC/DC voltage-source converters. IEEE Trans. Power Electron. 2003, 31, 11–22. [Google Scholar]
  9. Kazmierkowski, M.P.; Krishnan, R.; Blaabjerg, F. Control in Power Electronics—Selected Problems; Academic Press: Cambridge, MA, USA, 2002. [Google Scholar]
  10. Lee, T.S. Lagrangian Modeling and Passivity Based Control of Three Phase AC to DC Voltage Source Converters. IEEE Trans. Ind. Electron. 2004, 51, 892–902. [Google Scholar] [CrossRef]
  11. Bouafia, A.; Gaubert, J.P.; Krim, F. Predictive Direct Power Control of Three-Phase Pulsewidth Modulation (PWM) Rectifier Using Space-Vector Modulation (SVM). IEEE Trans. Power Electron. 2010, 25, 228–236. [Google Scholar] [CrossRef]
  12. Rodriguez, J.; Cortes, P. Predictive Control of Power Converters and Electrical Drives; Wiley-IEEE Press: Hoboken, NJ, USA, 2012. [Google Scholar]
  13. Vazquez, S.; Sanchez, J.A.; Carrasco, J.M.; Leon, J.I.; Galvan, E. A model based direct power control for three-phase power converters. IEEE. Trans. Ind. Electron. 2008, 55, 1647–1657. [Google Scholar] [CrossRef]
  14. De Araujo Ribeiro, R.L.; de Oliveira Alves Rocha, T.; de Sousa, R.M.; dos Santos, E.C.; Lima, A.M.N. A Robust DC-Link Voltage Control Strategy to Enhance the Performance of Shunt Active Power Filters Without Harmonic Detection Schemes. IEEE Trans. Ind. Electron. 2015, 62, 803–813. [Google Scholar] [CrossRef]
  15. Ge, J.; Zhao, Z.; Yuan, L.; Lu, T.; He, F. Direct Power Control Based on Natural Switching Surface for Three-Phase PWM Rectifiers. IEEE Trans. Power Electron. 2015, 30, 2918–2922. [Google Scholar] [CrossRef]
  16. Salomonsson, D.; Sannino, A. Direct Power Control Based on Natural Switching Surface for Three-Phase PWM Rectifiers. In Proceedings of the Conference Record of the 2007 IEEE Industry Applications Conference, Forty-Second IAS Annual Meeting, New Orleans, LA, USA, 23–27 September 2007. [Google Scholar]
  17. Wang, C.; Li, X.; Guo, L.; Li, Y.W. A Nonlinear-Disturbance-Observer-Based DC-Bus Voltage Control for a Hybrid AC/DC Microgrid. IEEE Trans. Power Electron. 2014, 29, 6162–6177. [Google Scholar] [CrossRef]
  18. Lee, D.C.; Lee, G.M.; Lee, K.D. DC-bus voltage control of three-phase AC/DC PWM converters using feedback linearization. IEEE Trans. Ind. Appl. 2000, 36, 826–833. [Google Scholar]
  19. Komurcugil, H.; Kukrer, O. Lyapunov-based control for three-phase PWM AC/DC voltage-source converters. IEEE Trans. Power Electron. 1998, 13, 801–813. [Google Scholar] [CrossRef]
  20. Gomez, M.H.; Ortega, R.; Lagarrigue, F.L.; Escobar, G. Adaptive PI Stabilization of Switched Power Converters. IEEE. Trans. Control Syst. Technol. 2010, 18, 688–698. [Google Scholar] [CrossRef]
  21. Flores, D.D.P.; Scherpen, J.M.A.; Liserre, M.; de Vries, M.M.J.; Kransse, M.J.; Monopoli, V.G. Passivity-Based Control by Series/Parallel Damping of Single-Phase PWM Voltage Source Converter. IEEE. Trans. Control Syst. Technol. 2012, 3, 459–471. [Google Scholar]
  22. Krause, P.; Wasynczuk, O.; Sudhoff, S. Analysis of Electric Machinery; IEEE Press: Baltimore, MD, USA, 1995. [Google Scholar]
  23. Sul, S.K. Control of Electric Machine Drive Systems; Wiley: Hoboken, NJ, USA, 2011; Volume 88. [Google Scholar]
  24. Mathew, S. Wind Energy: Fundamentals, Resource Analysis and Economics; Springer: Berlin, Germany, 2006. [Google Scholar]
  25. Khalil, H.K. Nonlinear Systems; Prentice Hall: Upper Saddle River, NJ, USA, 2002. [Google Scholar]
Figure 1. The permanent magnet synchronous generator (PMSG) power system topology.
Figure 1. The permanent magnet synchronous generator (PMSG) power system topology.
Energies 10 01387 g001
Figure 2. The implementation of the closed-loop system
Figure 2. The implementation of the closed-loop system
Energies 10 01387 g002
Figure 3. The closed-loop performance comparison between the proposed method and the feedback linearizing (FL) method with a resistive load of R L = 100 Ω .
Figure 3. The closed-loop performance comparison between the proposed method and the feedback linearizing (FL) method with a resistive load of R L = 100 Ω .
Energies 10 01387 g003
Figure 4. The estimated disturbance behaviors and the wind pattern based on the Weibull distribution.
Figure 4. The estimated disturbance behaviors and the wind pattern based on the Weibull distribution.
Energies 10 01387 g004
Figure 5. The comparison result of DC-link voltage regulation performance between the proposed method and the feedback linearizing (FL) method at the DC-link voltage of 300 V with a pulse resistive load torque from R L = 100 Ω to R L = 28 . 6 Ω ; (a) the resulting output voltage wave forms, (b) a magnified version of (a), (c) d-axis current wave forms, (d) q-axis current wave forms
Figure 5. The comparison result of DC-link voltage regulation performance between the proposed method and the feedback linearizing (FL) method at the DC-link voltage of 300 V with a pulse resistive load torque from R L = 100 Ω to R L = 28 . 6 Ω ; (a) the resulting output voltage wave forms, (b) a magnified version of (a), (c) d-axis current wave forms, (d) q-axis current wave forms
Energies 10 01387 g005
Figure 6. The permanent magnet synchronous generator (PMSG) output power behaviors of the proposed method and the feedback linearizing (FL) method at the DC-link voltage of 300 V with a pulse resistive load torque from R L = 100 Ω to R L = 28 . 6 Ω .
Figure 6. The permanent magnet synchronous generator (PMSG) output power behaviors of the proposed method and the feedback linearizing (FL) method at the DC-link voltage of 300 V with a pulse resistive load torque from R L = 100 Ω to R L = 28 . 6 Ω .
Energies 10 01387 g006
Figure 7. The DC-link voltage tracking performance variation comparison between the proposed method (a) and the feedback linearizing (FL) method (b) at three resistive loads: R L = 100 , 50 , 30 Ω .
Figure 7. The DC-link voltage tracking performance variation comparison between the proposed method (a) and the feedback linearizing (FL) method (b) at three resistive loads: R L = 100 , 50 , 30 Ω .
Energies 10 01387 g007
Figure 8. The DC-link voltage tracking behavior comparison between the proposed method (a) and the feedback linearizing (FL) method (b) by increasing the cut-off frequencies, such as 2, 5, 8 Hz.
Figure 8. The DC-link voltage tracking behavior comparison between the proposed method (a) and the feedback linearizing (FL) method (b) by increasing the cut-off frequencies, such as 2, 5, 8 Hz.
Energies 10 01387 g008

Share and Cite

MDPI and ACS Style

Kim, S.-K. Proportional-Type Performance Recovery DC-Link Voltage Tracking Algorithm for Permanent Magnet Synchronous Generators. Energies 2017, 10, 1387. https://doi.org/10.3390/en10091387

AMA Style

Kim S-K. Proportional-Type Performance Recovery DC-Link Voltage Tracking Algorithm for Permanent Magnet Synchronous Generators. Energies. 2017; 10(9):1387. https://doi.org/10.3390/en10091387

Chicago/Turabian Style

Kim, Seok-Kyoon. 2017. "Proportional-Type Performance Recovery DC-Link Voltage Tracking Algorithm for Permanent Magnet Synchronous Generators" Energies 10, no. 9: 1387. https://doi.org/10.3390/en10091387

APA Style

Kim, S. -K. (2017). Proportional-Type Performance Recovery DC-Link Voltage Tracking Algorithm for Permanent Magnet Synchronous Generators. Energies, 10(9), 1387. https://doi.org/10.3390/en10091387

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