Research Article  Open Access
YuHsing Chen, Catalin Gabriel Dincan, Philip Kjær, Claus Leth Bak, Xiongfei Wang, Carlos Enrique Imbaquingo, Eduard Sarrà, Nicola Isernia, Alberto Tonellotto, "ModelBased Control Design of Series Resonant Converter Based on the Discrete Time Domain Modelling Approach for DC Wind Turbine", Journal of Renewable Energy, vol. 2018, Article ID 7898679, 18 pages, 2018. https://doi.org/10.1155/2018/7898679
ModelBased Control Design of Series Resonant Converter Based on the Discrete Time Domain Modelling Approach for DC Wind Turbine
Abstract
This paper focuses on the modelling of the series resonant converter proposed as a DC/DC converter for DC wind turbines. The closedloop control design based on the discrete time domain modelling technique for the converter (named SRC#) operated in continuousconduction mode (CCM) is investigated. To facilitate dynamic analysis and design of control structure, the design process includes derivation of linearized statespace equations, design of closedloop control structure, and design of gain scheduling controller. The analytical results of system are verified in zdomain by comparison of circuit simulator response (in PLECS™) to changes in pulse frequency and disturbances in input and output voltages and show a good agreement. Furthermore, the test results also give enough supporting arguments to proposed control design.
1. Introduction
MEDIUMvoltage DC (MVDC) collection of wind power is an attractive candidate to reduce overall losses and installation cost, especially within offshore HVDCconnected wind generation as illustrated in Figure 1 [1]. To connect DC wind turbine with MVDC network (±50), the series resonant converter (SRC) serves as a stepup solidstate transformer as shown in Figure 2. With the series resonant converter, the DC turbine converter can take advantages of high efficiency, high voltage transformation ratio, and galvanic fault isolation for different ratings of turbine generator [2–6].
Traditional closedloop control of SRC for the DC distribution system is easily implemented by detecting the zerocrossing of the resonant inductor current and controlling the length of transistor and diode conduction angle α without considering circuit parameters of SRC [7]. Additionally, the output power flow control of SRC for DC network is achieved by controlling the phaseshift angle and frequency between the two arms of Hbridge inverter [6, 8, 9].
Based on the discrete time domain modelling approach, the smallsignal model of an improved SRC (named SRC#) is proposed [9, 10]. This paper continues with the smallsignal plant model addressed in Section 3 and the Appendix and mainly focuses on the closedloop control design for the system. In the following sections, the mode of operation of SRC# and smallsignal plant model based on the discrete time domain modelling approach will be briefly introduced first. The structure of closedloop control based on the proposed smallsignal plant model and the improvement in the disturbance rejection capability will be revealed. To satisfy the power flow control with variable switching frequency, the gain scheduling technique will be given. Finally, the analytical solution of overall system is revealed and verified by comparing with timedomain trace in circuit simulation model implemented in PLECS™ under different operating points. Furthermore, the proposed control deign will be demonstrated by a scaleddown laboratory test bench.
2. Mode of Operation of Series Resonant Converter
The mode of operation of series resonant converter (SRC#) in Figure 2 is decided by the ratio between natural frequency of tank ( and ) and the switching frequency of Hbridge inverter: subresonant, resonant, and super resonant mode. In subresonant mode, the switching frequency of Hbridge inverter is lower than the natural frequency of tank. The resonant operating mode is selected when the switching frequency is equal to the natural frequency of tank. If converter’s switching frequency is higher than the natural frequency of tank, the converter is operated in the super resonant mode [9].
Contrasting with the constant frequency with phase shift control which is normally applied for operation in super resonant mode, to achieve ZVS at turnon, Figure 3 illustrates the concept of frequencydepended power flow control of SRC#. The converter leg of SRC# consisting of switches T_{1} and T_{2} is referred to as the leading leg and the one consisting of switches T_{3} and T_{4} is referred to as the lagging leg as indicated in Figure 2. Both converter legs operate at a 50% duty cycle [6, 9].
To achieve ZCS character at turnoff or minimize the turnoff current, the IGBTbased SRC# is designed to operate at subresonant continuousconduction mode (subresonant CCM). This control design can drive the implemented phase shift having the same length as the resonant pulse without sacrificing the advantage of linear relation to the number of resonant pulses, as depicted in Figure 3. Compared to a traditional SRC with frequency control design in subresonant mode, therefore, the medium frequency transformer in the SRC# addressed in this paper can be designed for a higher frequency and avoids saturation for lower frequencies.
3. Discrete Time Domain Modelling Approach for Series Resonant Converter
Considering the efficiency, subresonant mode is selected for the mode of operation of SRC# for the DC wind turbine [5, 6]. Based on the circuit topology shown in Figure 2, Figure 5 illustrates the steadystate voltage and current waveforms of SRC# in subresonant mode, where is the switching frequency (=2π·) of SRC#. To apply the linear control theory to the SRC# control design, deriving the plant model of SRC# with the discrete time domain modelling approach includes the derivation of largesignal equations based on the interesting interval shown in Figure 5, linearization of discrete state equations, and derivation of smallsignal transfer function. In the derivation, the voltages and are assumed to be discrete in nature, having the constant values and in interval of event, and then switch to next states and at the start of event. This procedure is only valid when the variation in or v_{g}(t) in the event is relatively smaller than its initial and final values [10].
With the discrete time domain modelling approach, (1) gives a linearized statespace model of SRC# in subresonant mode and the transfer functions between input state variables and the defined interesting states are shown in (3) and (5). To simplify the derivation, the output filter of SRC# (i.e., L_{f} and C_{f}) is neglected and only the DC component of output current diode rectifier i_{out,Rec} is selected as an output variable I_{o}.
To obtain the harmonic model of DC turbine converter, Figure 4 gives a complete flow chart of mathematical derivation of SRC# plant model, which describes how the SRC# plant model is obtained. First of all, the circuit topology and mode of operation are decided as shown in Figures 2 and 5 and then the equivalent circuit based on the switching sequence of transistors is generated in Figure 6. Based on the circuit topology shown in Figure 2, Figure 5 illustrates the voltage and current waveforms of SRC# n subresonant mode and the equivalent circuit for each event (switching interval) is given in Figure 6. According to Figure 6, the large signal model of converter is created (step 4) and then the interesting state variables (step 5) are defined to generate the smallsignal equation and the space model of converter as in steps 6 and 7, respectively (, , , and ). Eventually, the converter plant model (power stage of converter) is established based on the interesting transfer function (g_{1}, g_{2}, and g_{3}, in step 8). The correction of model (plant model) has been confirmed and the details of derivation are given in the Appendix.whereand the derivation of [A], [B], [C], and [D] matrixes is shown in the Appendix. Transfer functions between defined internal state variables and input states are given bywhereand transfer functions between converter output current and input state variables are where the transfer functions , , and can be obtained via The transfer function describes how the output current is influenced by the control input variable and the transfer functions and describe how the output current is affected if any disturbance occurs in input voltage () and the output voltage (). For example, the array network (MVDC grid) contains voltage harmonics. The transfer function can be used to evaluate the effect of voltage harmonics on the converter output current. Detailed derivation of the above linearized statespace model and the expression of elements in [A], [B], [C], and [D] matrix in (1) have been revealed in the Appendix.
(a)
(b)
(c) (T2, T3 ON)
(d) (D2, T4 ON)
4. ModelBased ClosedLoop Control Design
Figure 7 gives an overview of smallsignal control model of SRC# based on the average plant model in (5) and (6). The control design of SRC# includes derivation of smallsignal plant model and the design of the compensator . The smallsignal transfer functions of SRC# between converter output current and input state variables are given by (5) and (6), where the output current variation is the expression of linear combination of the three independence inputs. The relationship between α and f_{s} in subresonant CCM in largesignal model is whereBy substituting the perturbation terms of smallsignal analysis into expression in (8), the smallsignal expression of α and f_{s} can be obtainedwhere the AC component isEventually, the system transfer function in Figure 5 can be expressed asEquation (12) can be further expressed as the following:with a loop gain.where the loop gain is defined by the product of gains around forward and feedback paths [11].
5. Disturbance Rejection Capability
The closedloop control design of SRC# is implemented via the compensator , which is applied to shape the loop gain of the system (i.e., ). Considering the transfer function of output current given in (13), the relationship between and is shaped by closedloop control asThe variation in output current I_{o} caused by can be alleviated by increasing the magnitude of the loop gain when the closedloop control design is integrated with the SRC# plant model. The system transfer functions in (13) also show that the variation reduction of I_{o} due to variation in MVDC network will benefit from a high loop gain :Furthermore, consider the tracking performance of output current control in (17).Assume that a constant power reference is applied to the control loop with a constant MVDC source and a constant LVDC source. A large loop gain (i.e., ) can also make sure of a good DC current tracking performance as shown in (18)
Therefore, the objective of the compensator is to govern the system with a desired loop gain (i.e., ), where the deviation of desired loop transfer function can be found by simply evaluating the magnitude asymptote in Figure 8:Considering the desired loop gain illustrated in Figure 8, the disturbance rejection capability of the output current for a frequency range below the crossover frequency () can be improved with closedloop control. For example, at the low frequency range (), the output current is almost in direct proportion to the power reference signal . Furthermore, a high loop gain provides a good disturbance reduction to the variation on input voltage and output voltage by the factor 1/.Typically, the crossover frequency f_{c} should be less than approximately 10% of switching frequency of SRC# ( < 0.1) to limit the harmonics caused by PWM switching [11]. Based on (19), therefore, compensator under a certain operating point (OP) can be expressed by
Equations (23)–(27) summarize the parameters (i.e., Q, , , , , and ) which are used to shape the loop gain via the compensator . The crossover frequency f_{c} and the lowfrequency pole at are defined asThe lowfrequency zero at and highfrequency pole at can be chosen according to crossover frequency and required phase margin as follows:where the angle is a phase lead angle of compensator at f_{c}. The DC gain of target loop gain is
The Qfactor is used to characterize the transient response of closedloop system. Using a high Qfactor can increase the dynamic response during transient, but it can also cause overshoot and ringing on power devices. In practical application, the Qfactor must be sufficiently low to keep enough phase margins and alleviate voltage and current stress on power devices [11]. Additionally, since the power flow control of SRC# depends on the control of switching frequency f_{s}, the parameters of target curve and the coefficient of transfer function g_{c} have to be changed according to different operating points (different output powers). To make sure that the compensator g_{c} can match with different output power requirements, therefore, a gain scheduling approach is proposed which will be revealed in the next section.
6. Design of Digital Gain Scheduling Controller
Gain scheduling controller is designed to access the parameter of compensator in real time and then adjust it based on the different operating points. Figure 9 gives a complete digital controller of SRC# based on the smallsignal control model and the bilinear transformation. The digital controller of SRC# consists of a smallsignal controller, a gain scheduling controller, a feedforward control loop, and a DC component calculator ( calculator). The controller is implemented in zdomain with a variable interrupt frequency (∝ switching frequency ). With the bilinear transform, the general form of the discretetime representation of the compensator can be expressed aswhere coefficients and (n=0~5) are used to specify the coefficients of numerator and denominator.
To design the gain scheduling controller, coefficients a_{n} and b_{n} in (28) are evaluated under different operating points (i.e., different output power) with (22)–(27). A trend in the variation of each coefficient (i.e., P_{REF} vs. a_{n} and b_{n}) is recorded and then is formulated via the polynomial approximation as shown in Figures 12 and 13 which will be discussed in the next section (Section 7). Eventually, the coefficient of g_{c}(z) for SRC# in the subresonant CCM can be adjusted by a continuous function such as and in real time to avoid any potential turbulences caused by gainchanging.
7. Verification of ClosedLoop Control Design
With the SRC# topology in Figure 2 and the controller shown in Figure 9, Tables 1 and 2 give the parameters used in the statespace model and circuit simulation models (tools) for verifying the validation of overall system in zdomain. The control model in the subresonant CCM is verified to identify the accuracy of proposed smallsignal model, and then the results of coefficient assessment of with the gain scheduling controller are integrated with control loop and are tested by a ramppower reference.


By applying a +0.5% stepping perturbation to all input state variables, Figures 10 and 11 give the analytical solutions of smallsignal model of SRC# and the results obtained from the timedomain switching model implemented in PLECS™. The SRC# with closedloop control is commanded to deliver around 9.0MW DC power and 7.5MW DC power to MVDC network, respectively. Figures 10 and 11 show that both the steady state and transient state in the analytical model match with the results generated by switching model. Therefore, dynamics of SRC# switching model can be predictable and controlled with the proposed smallsignal model.
Figures 12 and 13 give the result of coefficient assessment of for the design of the gain scheduling controller. Based on (28), the trend in the variation of coefficients and in subresonant CCM from 5.75MW to 10MW (0.5MW/step) is identified and then the variation of each coefficient is approximated with a polynomial (i.e., and ). According to the variation in output power reference , the gain scheduling controller accesses the polynomial to regulate its coefficient in real time. To evaluate the adequacy of control design of overall system, finally, the timetrace simulation of output power flow control is given in Figure 14 with a ramppower reference from 0.1MW to 10MW, and vice versa. The results show that the output current/power () of the series resonant converter can be well controlled when magnitude output powers references are changed.
(a) Rampup power: ==0.1MW → 10MW, slop=4.95MW/sec
(b) Rampdown power: ==10MW → 0.1MW, slop=4.95MW/sec
8. Laboratory Test Results
To verify the control design, first the circuit simulation is carried out with circuit simulation tool of PLECS™, and then the controller is implemented in a scaleddown laboratory test bench. The circuit configuration of test bench and the corresponding parameters are shown in Figure 15 and Table 3, respectively, where the MVDC network is simulated by a unidirectional power flow DC power source with a controllable perturbation.

Figures 16(a) and 16(b) depict the system response when a positive and a negative step perturbation (0.01p.u) in MVDC network are applied, respectively. The control design exhibits a close behavior in either simulation or experimental test. There is some small tracking error during the transient between the simulation and test results. This usually is caused by the estimated error of components and stray inductance which is not considered in simulation model. Figure 16(c) represents how the output current behaves when a stepchange (0.26p.u) is applied in the power reference signal . Under the proposed control law for SRC#, both the simulation and test result show that the DC component of DC turbine output current () tracking performance can be guaranteed. However, a small oscillation (≈40Hz) during the transient of stepchange of power reference signal in the experimental test is observed due to the series diode (in Figure 15) which is reversebiased at this test occasion.
(a) 4.0 stepup disturbance (0.01p.u) in the MVDC voltage ()
(b) 4.0 stepdown disturbance (0.01p.u) in the MVDC voltage ()
(c) Stepup disturbance (0.26p.u) in the power reference signal (equivalent of current reference signal: 1.37A→1.73A)
9. Conclusion
A modelbased control design of SRC# for DC wind power plant based on smallsignal plant model in the discrete timedomain modelling is revealed. This paper continues with the modelling of SRC# given in the Appendix and mainly addresses the closedloop control design for the system. The control design process contains the derivation of statespace plant model, design of closedloop control structure, and design of gain scheduling controller. Compared with the traditional frequencydepended power flow control which relied on openloop structure, the SRC# with the closedloop structure can gain a better disturbance rejection capability for the output power control. The verification of proposed digital controller including plant model is addressed in both the analytical model and the timedomain circuit simulation implemented in PLECS™ in Section 7 by evaluating the SRC# with the steppingperturbation under the subresonant CCM. Furthermore, gain scheduling approach is implemented by the polynomial approximation and tested under different operating points (different output powers). Integrating the gain scheduling controller with closedloop structure enables the system to automatically adjust parameters of controller in real time to satisfy different output power requirements without sacrificing the control performance. Finally, Section 8 shows that all the test results give enough supporting arguments to the proposed control design.
Appendix
The objective of the study is to understand the harmonics distribution of offshore DC wind farm and how the DC wind turbines are affected by harmonics from MVDC gird. This section summarizes the derivation of plant model of DC wind turbine based on the discrete time domain modelling approach (discrete time domain modelling approach [10], steps 18) which can help the reader to reach the plant model of DC wind turbine (SRC#) and then conduct control deign of DC wind turbine. The following discussion will give a complete derivation process including the corresponding flow chart of the derivation of SRC# plant model given in Figure 4.
Steps 1 3: Decide the Circuit Topology of DC Turbine Converter, Resonant Tank Waveform, and Equivalent Circuit. Steps 13 describe the circuit topology of SRC# (DC wind turbine converter) and mode of operation, which is operated in subresonant CCM as in Figures 2 and 5. The corresponding equivalent circuit for the SRC# in subresonant CCM is given in Figure 6, where the waveform is divided by different time zone (different switching sequence) based on the discrete time domain modelling approach proposed by King, R. J. [10]. Those figures (Figures 2, 5, and 6) are used to generate the large signal model of SRC#.
Step 4: Large Signal Model. Based on Figure 6, the objective of derivation of largesignal model is to express the final value of interesting state variables in each switching interval with the initial values. The procedure is only valid when the variation in output voltage (MVDC grid voltage) or input voltage (LVDC voltage) in the event (switching) is relatively smaller than its initial and final values [10]. Equations (A.1) to (A.16) give the derivation of largesignal model of resonant inductor current i_{r}(t) and resonant capacitor voltage v_{Cr}(t) and their end values at event in terms of initial values of event.
For (T_{1}, T_{4} ON)whereThe resonant inductor current (t) and resonant capacitor voltage v_{Cr}(t) can be obtained by solving (A.1).whereAt time the tank current i_{r} makes a zero crossing, commutating T_{1} and T_{4} off and turning on D_{1} and T_{3}. Therefore,whereFor ≤ t ≤ (D_{1}, T_{3} ON)wherewhereEventually, the inductor current i_{r}(t) and capacitor voltage v_{Cr}(t) at time t= can be represented bywhereThe largesignal expression of resonant inductor current i_{r}(t) and resonant capacitor voltage v_{Cr}(t) in event ( and ) can be obtained with the same process as derivation of equations, as in (A.1)(A.16).
SteadyState Solution of LargeSignal Model. Equation (A.18) gives the conditions for calculating steadystate solution (operating points) of discrete state equation.By substituting (A.18) into (A.15) and (A.16), the steadystate solution of and can be expressed in terms of , , β_{k}, and α_{k}:where the overbar is used to indicated the steadystate value of interesting state variables.
To simplify the derivation, the output filter of SRC (i.e., L_{f} and C_{f}) is neglected due to very slow dynamics in voltage and current compared with the resonant inductor current and resonant capacitor and only the DC component of output current diode rectifier i_{out,Rec} is selected as an output variable i_{o}. Therefore, during the K_{th} event, the output current equation delivered by the SRC is expressed aswhereand is the initial value of Inductor current, is the initial value of capacitor voltage, is the initial value of rectifier output voltage, and is then initial value of input voltage of SRC#. Z_{r} ( =√(L_{r}/C_{r}) ) is characteristic impedance defined by parameter of resonant tanks, α_{k} (= ) is the transistor and diode conduction angle during the switching interval (event k), and θ_{s} (=) is represented by the switching frequency of converter.
The steadystate solution of discrete state equation for output variable i_{o} is obtained by substituting steadystate condition into (A.21) as
Step 5: Define State Variable. Since the discrete largesignal state equations in (A.15), (A.16), and (A.21) have a high nonlinearity, control design technique based on the linear control theory cannot directly be applied. To obtain a linear statespace model, therefore, the linearization of largesignal equation is necessary. Equation (A.24) gives the definitions of interesting state variables in both the k_{th} switching event ( ≤ t ≤ ) and the switching event ( ≤ t ≤ ). Finally, the equations of approximation of derivative in (A.24) and (A.25) are used to convert the discrete stateequation (largesignal model) into continuous time [10].whereBy replacing the state variables in (A.15) and (A.16) with the defined state variables in (A.24) and applying the approximation of (A.25) for derivative, the nonlinear statespace model is given bywhere the output equation is defined as
Step 6: Linearization and SmallSignal Model. Consider that all the interesting state variables in pervious steps are in the steadystate (near the certain operating point, OP) with a small perturbation; therefore, the nonlinear state equations can be formalized with Taylor Series Expansion in terms of the operating point (OP) and the perturbations:
(i) Resonant inductor current: whereand thenwhere the subscript OP indicates the steadystate point, where the derivatives are evaluated at that point.
(ii) Resonant capacitor voltage:whereand then
where
(iii) Output current equation:whereand thenwhere