Abstract

This paper presents a stability analysis of train converters in order to evaluate how the controller parameters affect the absorbed current. The new dynamic model presented in this paper is capable of considering the time-variant nature of the system for the correct tuning of the feedback proportional-integral PI controller, applying a current controlled modulation technique never used in high-power traction converters. The reduction of the harmonic content of the current absorbed by a converter employed at the input stage onboard high-speed trains is really important, considering the interaction with the signaling system set up for traffic control. A computer model of the converter, considering both the power and the control structure, has also been implemented in order to deliver a validated tool for the developed theoretical analysis.

1. Introduction

The main purpose of public transport is to ensure, in every operational situation, the safety and regularity of the service provided by rolling stocks. To this purpose it is necessary that the harmonic disturbances generated by traction units are conformed to the EMC limits imposed by the compatibility with the signalling system set up for traffic control.

The propulsion motor drives and auxiliary services of the modern traction units are powered by microprocessor controlled bidirectional electronic converters. Therefore, there is a continuous energy exchange with the power supply system in large spectral contents, since the new traction units cannot be considered anymore as simple passive loads absorbing energy from the line.

Nowadays the input stage of the locomotives and high speed trains supplied in AC at 25 kV, 50 Hz and 15 kV, 16.7 Hz are constituted by more four-quadrant (4Q) converters. In fact, the high power requested by the traction motor drives and the DC/DC converters for auxiliary services request a stable input DC voltage in the range of 1500–1800 V to better exploit the modern IGBT switches. This DC link voltage has to be provided by 4Q converters that allow its stabilization and the power inversion during regenerative braking. The high power value of the input stage does not allow high switching frequencies for the 4Q converter. Consequently the absorbed current presents a high ripple value characterized by high harmonic current components that cannot be tolerated by the system. Indeed the track circuit used for signaling and communication for the traffic management and safety employs signal currents overlapped with the power ones. For this reason it is really important to control the current harmonic content in order to prevent the possible interference phenomena and to meet the requests of the new Technical Specification for the Interoperability [1]. First of all, it is necessary to perform a stability analysis of the 4Q converter. The paper presents a mathematical model for the single-phase four-quadrant (4Q) converter, introducing a method able to consider the time-variant nature of the system. The main goal of this algorithm is the analysis of the control parameter effects on the absorbed current in order to perform a correct tuning of the feedback PI controller. The developed method applies a particular modulation technique instead of the traditional PWM control, in order to better avail the converter characteristic, controlling directly the absorbed current and maintaining, at the same time, a constant switching frequency. Therefore, the main characteristic of the proposed algorithm is its application to a single-phase converter that employs a current controlled control strategy (Smart Modulation). This allows to linearize the system in order to apply the traditional methodology for the stability analysis. Other algorithms employ vectorial techniques, Park Transform, and so forth, but they are useful for three-phase converters or other methodologies are applied to the traditional PWM modulation.

Starting from real cases referred to European High-Speed Trains, a circuital model of the train input stage has been implemented in EMTP environment, considering both the power and the control structure. Many simulations have been performed to validate the theoretical analysis presented.

2. AC/DC Conversion Input Stage

The input stages of a High-Speed Train (HST) is composed of dedicated converters that realize the bidirectional AC/DC conversion. The HST here considered is constituted by two locomotives with two motor drives each. Every motor drive is supplied by double 4Q converters input stage, in order to realize redundancy and then guarantee the continuity of the service, through the main transformer as depicted in Figure 1.

The 4Q converter function is to absorb/deliver power from/to the AC network according to the traction motors and auxiliary services power demand [2]. First of all, the converter has to absorb current from the contact line at a voltage and frequency with a fundamental harmonic, , in phase with and with limited harmonic disturbance, that means a power factor and a distortion factor . Moreover, it has to absorb power from the contact line with a mean value, , pulsing at the frequency , as well as to supply a continuous power, , to the three-phase motor drive inverters connected to the direct current DC link [3].

The use of a forced switching converter allows for higher DC link voltages than those generated by diode or thyristor bridges. In this manner, it is possible to increase the flowing power at the same converter losses [4]. The maximum value of the DC link voltage is imposed by the presence of the switches and the voltage smoothing capacitors [5]. Moreover, the filtering on the AC side is simplified and is only related to the switching frequencies [6].

3. Mathematical Model of the 4Q Converter

The 4Q converters used in modern locomotives are two-level converters, because they employ the new high voltage IGBT. Therefore, their stability analysis needs to consider the four different circuit configurations as a function of the switch status. Figure 2 reports a representation of the two converter levels considered in which its equivalent circuit is composed of an ideal voltage generator (, where is the transformer ratio of the onboard transformer), a commutation inductance , an equivalent resistances , a DC link capacitance , and an equivalent current absorbed by the train motor drives .

The modeling process starts by representing the discrete modes of the converter mathematically. The discrete model of the converter describes each status of the converter by separate linear equations. Taking the state variables as and , the four-converter state equations are reported in Table 1.

Observing the matrices , , and , it is possible to note that the working conditions do not influence the elements on the diagonal. Therefore, the state equation in generic condition can be reduced to where depends on the mode of the converter and can be equal to 1, −1, or 0.

The 4Q converter function is to absorb or deliver power from/to the AC line depending on the power value requested by the traction motors and auxiliary services. This objective can be achieved using different control logic [7]. For this application PWM is often used, which is a well-known commutation technique. This control has the advantage of having a constant switching frequency, which allows for the sizing of dedicated filters tuned at that frequency; however, this control has a disadvantage in that it has indirect control of the current through the variation of the modulation index, , and the load angle, , causing possible transient overcurrents.

Assuming the use of PWM modulation, there are many commutations in a period at the main frequency, and, therefore, can be replaced with the time-varying function, , where is the modulation index and is the load angle. The time-varying nature of complicates the resolution of (1).

To eliminate the time variance in the equations, instead of transforming the reference frame and applying the Park Transform, a substitution method has been adopted [8]. The first step is to rewrite the state equations including : From (2) a solution for and can be where has one constant forcing term and another sinusoidal term at the primary angular frequency, , and has one constant component and a ripple at an angular frequency of . Both the constant term and phase are forced to zero by the converter control. The ripple component, , is quite low due to an LC filter in the DC link that is tuned at a frequency, , as well as a high capacitance and voltage in the DC link. As it will be presented in the next session, it reaches the 0.2% of the DC component, therefore, the voltage can be approximated with its average value .

Deriving (3) and substituting it into (2), a new system of three equations in three unknown terms can be obtained equating the cosine and sine terms: These equations are time invariant, but nonlinear. To be implemented in a control system design they must be linearized around their steady-state operating point. To linearize the system, the following substitutions are made, where capital letters represent the steady-state operating point and the prime symbol represents a small deviation from the operating point: (i) modulation index;(ii) DC link voltage;(iii) load angle;(iv) absorbed current.

The value of is stiff and imposed by the contact line connected to the main line, whereas is imposed by the power requested by the train.

The first-order approximations from the Taylor Series for the sine and cosine functions are and used to simplify the system of (4):

4. Stability Analysis of 4Q Converter with Current Controlled Strategy

Retaining only the first-order terms, it is possible to obtain the small-signal model, useful to find the frequency response of the converter [9]: The small-signal model is developed in terms of the Laplace operator. Applying the Laplace transformation to (6) and assuming zero initial conditions, we obtain Instead of the PWM control that generates an AC voltage trough the modulation index and its phase , in the following a new modulation technique that can directly control the current value is proposed.

Considering the variability of the switching frequency, , and the difficulty defining an adequate current ripple, the Smart Modulation, presented in [10, 11] and briefly described in the appendix, is applied.

This strategy allows direct control of the current value maintaining, at the same time, a constant value of the switching frequency. Therefore, the modulation index, , does not directly appear in the control algorithm, but, assuming a unitary power factor, the converter absorbed current is controlled by the equivalent conductance, . From the second equation of the system in (7) it is possible to determine the value of . Substituting into the third equation of (7) we obtain Equation (8) gives the value, , as a function of . It is now possible to determine the variations of the DC link voltage for small changes of the current absorbed by the converter The reference current is given by the product between the equivalent conductance, , and the supply voltage, (), where G comes from the PI regulator that maintains the DC link voltage constant: The block diagram of the DC voltage regulation loop is shown in Figure 3, where only one of the two-input stage 4Q converter is represented. In the block diagram shown in Figure 3 three control loops can be identified. The first control loop is the DC link voltage control loop comprised of the voltage measurement, the low-pass filter explained above, the comparison with the reference value, , and the PI controller. Its output is the value of the equivalent conductance, , that keeps the DC link voltage constant while varying the power requested or injected by the traction motors and auxiliary services. The second loop is the reference current generator, which considers the input voltage measurement followed by a filter dedicated to high-frequency disturbances. The obtained value from this loop multiplied with the equivalent conductance, , gives the reference current that the converter, through Smart Modulation, has to generate in order to balance the input and output powers. The third loop is related to the DC component compensation in the AC input current. In fact, its output value is a constant current that, when added to the reference current, allows the DC component to be cancelled out in order to avoid saturation of the input transformer.

The DC link voltage variation for small changes of given by the regulator is obtained multiplying (9) for the input voltage : The open-loop function given by the product of (10) and (11) is The above presented theory leads to a procedure for the converter design. In particular, the following steps can be identified:(1)converter parameters identification, such as rated and maximum powers, AC input and DC output voltages, and consequently input and output currents;(2)switching inductance and DC link capacitance calculation. These two parameters are determined considering the maximum input current ripple and the maximum acceptable DC link voltage oscillation;(3)calculation of the correct tuning parameters of the regulator (proportional and the integral ) to have good stability and at the same time a fast time response of the system;(4)calculation of the tuning parameters of the regulator (proportional and the integral ) at the limit of stability in order to perform the sensitivity analysis;(5)sensitivity analysis.

5. Analysis Application in a Real Case

The converter design steps previously formulated have been applied to a real case as example consisting in the converter studied in this paper (Figure 1) employed as input stage onboard High-Speed Railway (HSR) Trains operating on the new Italian HSR lines.

Converter Parameters Identification.
In transportation systems many typologies of power converters are employed that can have different applications (i.e., front-end converters, auxiliary services, UPS, high-speed trains, heavy traction, naval propulsion) and sizes. The typical size of the converter and its usual data for this application are reported in Table 2.

Switching Inductance and DC Link Capacitance Calculation.
Referring to Figure 2, the elements that more influence the control parameters are the switching inductance and the DC link capacitance . The inductance value is chosen to limit the input current ripple, while the capacitance influences the DC link oscillation. Therefore, the obtained values are and .

Calculation of the Correct Tuning Parameters of the Regulator (Proportional and the Integral ).
In order to study the system stability using the Bode criteria through module and phase diagrams, identified by the open loop transfer function (12), a dedicated code in MATLAB environment has been implemented. From this program it is possible to calculate the control parameters and , giving as input the specific data of the converter considered (Table 2), that have to stabilize the converter dynamics. Under these conditions and referring to the studied case, a good compromise to obtain a fast system response and, at the same time, a good phase margin is and . The Bode diagrams obtained from (12) using these values are shown in Figure 4(a). The cut-off frequency that guarantees a good response with a sufficient phase margin is about 200 Hz.

Calculation of the Tuning Parameters of the Regulator (Proportional and the Integral ) at the Limit of Stability.
The parameter values for the limit of the stability condition have been investigated using the same Matlab code, obtaining and . The relative Bode diagrams for these parameters are shown in Figure 4(b). As can be seen from this figure, the phase margin is close to zero.

Sensitivity Analysis.
Applying the same algorithm, it is possible to obtain the effect of changes of the parameters on the stability of the controller. As it can be shown in Figure 5, an increase of the proportional gain causes significant reduction of the gain margin leading to the loss of stability. On the contrary, the reduction of improves the stability but it degrades the time response performance. Once determined the power system structure, the traction power converters are mainly subjected to the line voltage variation. In fact, in traction systems the voltage can vary in the range of −33%–+20% of the nominal one. As it can be seen in (12), the line voltage (reported on the secondary side of the onboard transformer) is common to all the numerator terms, therefore, it can be considered as a proportional gain. The sensibility analysis carried out and presented in Figure 5 shows that around the optimal value of (−33%–+20%) there is enough gain margin to maintain the stability.

6. Harmonic Absorption due to Single-Phase Supply

One of the problems that characterize the 4Q converter and, in general, all of the AC/DC single-phase converters, is the absorption of variable instantaneous power. In fact, in contrast to a balanced three-phase system, the flowing power in a single-phase system is characterized by a mean value equal to the active power requested from the load and a variable term oscillating at a frequency twice that of the supply frequency. This oscillation, in the case of traction converters, causes the generation of a ripple in the DC output voltage. The 4Q converters are equipped with a filter that is tuned to a frequency of 100 Hz in 25 kV, 50 Hz train systems or 33.3 Hz in 15 kV/16.7 Hz train systems. However, it is not possible to completely cancel the ripple in the DC link voltage. Due to the regulator response time, the equivalent conductance at the end of the transient state condition coming out from the regulator can be defined as the sum of its mean theoretic value requested and a ripple at double frequency   . This phenomenon has not been considered in the stability analysis because it is not influential; however, it is the cause of harmonic current absorption from the converter. The reference current leaving the regulator, given as the product between and the sinusoidal input voltage , is The reference current reported in (13) can be defined as the sum of the following three terms:(i) is the theoretical value of the reference current proportional to the absorption power; it is in phase with the voltage signal and at the same frequency;(ii) is a term at the same net frequency, but shifted π/2 radians, which is the source of absorption of the reactive power;(iii) is a term at three times the frequency that determines the presence of a third harmonic current.

In order to reduce the component, it is necessary to filter the DC link voltage with a suitable low-pass filter. This filter tuning is critical because it simultaneously allows for a fast system response and a low cut-off frequency. It is impossible to achieve both of these objectives, but for the system stability it is necessary to focus on the system response, which implies a small ripple.

7. Model and Simulation Results

In order to validate the theoretical analysis described above, a suitable model of the system has been implemented using the EMTP-ATP dynamic simulation tool. The data employed for the model is obtained from a real HSR Train operating in Italy on a 25 kV/50 Hz line. Each 4Q converter is sized for a rated power equal to 900 kW and a maximum peak power equal to 1500 kW. These are relevant power values for a switching converter that has to be small and light enough to be installed onboard.

The 4Q converters consist of a power circuit and a control system. The power circuit has been modeled with the circuital elements already available in EMTP-ATP. The MODELS language has been used for implementing the proposed converter control.

In the following, some simulation results are reported in order to verify the outcomes of the performed analysis. The simulation cases considered for the theoretical analysis validation are the following:(A)nominal condition;(B)loading;(C)converters regulators working at the stability limit.

(A) Nominal Condition
The normal operating condition considers the traction phase at the nominal power of 3.6 MW due to four 4Q converters. The parameters and obtained for the stability condition are employed in the model. The current absorbed by one 4Q converter is shown in Figure 6.
From Figure 6 it is easy to see that the converter commutations are numerically limited because of the low switching frequency of the power GTO due to the high rated power of the converter. Therefore, many harmonic components are generated in the input current. The most significant harmonic components are reported in Table 3.
In particular, considerable amplitude is associated to the 10th harmonic order, corresponding to the switching frequency of the converter (17% of the fundamental component). Moreover, there is the third harmonic component (which has an amplitude that is 2.7% of the fundamental amplitude) and is caused by the interaction between the ripple at 100 Hz of the DC link voltage and the supply voltage at 50 Hz. The THD is equal to 19.07%.
With respect to the DC link voltage, its behavior at the end of the transient state condition is shown in Figure 7.
The mean value of the voltage is equal to 1800 V, and it is maintained at a constant value equal to the reference value by the PI regulator. Figure 7 shows the clear oscillation at 100 Hz due to the oscillation of the power absorbed by the single-phase converter.

(B) Loading
The loading case, obtained applying a step load to the system, has been simulated in order to stress the system with high-power variation for stability test. Indeed, the transient behaviour has been analyzed in the case of a load change on the DC side from 50 A (presuming that 10% is dedicated to the auxiliary services functioning) to 500 A. This means that a power variation of 3.24 MW (from 360 kW to 3.6 MW) has been simulated. Also in this case, the parameters and obtained for the stability condition are employed in the model.
The absorbed current behaviour in the loading transient at (instant at which the initial simulation transient is finished) is shown in Figure 8.
Moreover, this power variation causes a DC link voltage swell equal to 30–35 V (1.6%), as shown in Figure 9 in continues line. In fact, in the transient phase, the unbalance between requested/remised power from the traction motors and absorbed/injected power from the 4Q converter, caused by the delay of the PI regulator, is compensated by the energy stored in the smoothing capacitors. The speed of the regulator is an important design parameter because it affects the dynamic of the DC voltage during transient events. A lower speed can cause high overvoltages that can damage the electronic components or deep voltage sags that can lead to loss of converter control.
The DC link voltage ripple implies that the measured value is not suitable for the regulator because it will cause the absorption of a high 3° harmonic current as shown in (13). For this reason, a filter that limits the 100 Hz component, but is sufficiently fast to follow the DC link voltage transients, is inserted. The filtered input signal of the converter regulator is shown in Figure 9 by the dashed line. In order to guarantee the stability of the system, this filtered signal has to limit the 100 Hz component before and after the transient phase maintaining, at the same time, a fast response to follow the transient phenomena.

(C) Converters Regulators Working at the Stability Limit
This last case has been simulated in order to verify the efficacy of the choice of the regulator parameters. For this reason, different values of and have been inserted in the model. In particular, the following examples will demonstrate the results obtained using the same values obtained previously for the limit of the stability condition: and .
The DC link current and the conductance at full load are demonstrated in Figures 10 and 11. For both of these behaviours it is possible to note a clear increase of an irregular ripple due to the control operation at the limit of stability. These results confirm the theoretical analysis carried out.
The equivalent conductance at the limit of stability assumes, at the end of the transient state condition, a mean value equal to the value at the stability condition because it is correlated to the power transfer requested by the loads. However, in this last case, there is a highly irregular oscillation around this mean value that, combined with the input sinusoidal voltage, causes the absorption of important current harmonic components (Table 4).

8. Conclusions

In this paper the influence of controller parameters on the quality of the absorbed current has been investigated. Moreover the application of a current control method to four-quadrant (4Q) input stage converters employed in traction systems has been proposed. In particular, it has been chosen the Smart Modulation technique with good results, that has been originally developed for low power high switching frequency converters. The changing of modulation technique, from the traditional PWM to the innovative Smart Modulation, has requested a new and deep stability analysis due to the changing of the mathematical equations that describe the whole system. This analysis has also led to a methodology for setting the regulator parameters.

In order to validate the theoretical analysis above described, a circuital model of real cases referred to European High Speed Trains has been implemented in EMTP environment. The outputs of the numerous simulations carried out confirm the results of the theoretical analysis and furnish important indications for the performance of the train and the observance of the European Technical Specification for the Interoperability [1]. In particular, they show how the oscillations in the control loop can cause great current harmonic absorption that can interfere with the track signalling system.

Appendix

The Smart Modulation (SM) technique is a predictive control that imposes the switching frequency independently from operating conditions. In addition, the SM allows interlacing of switching instants in case of more converters connected in parallel. Referring to Figures 2 and 3, this algorithm considers the reference current and the inductor flux error to calculate the commutation instants. The expression that defines the ideal voltage at the 4Q AC terminals when the reference current is absorbed is

Considering that the trajectories of the real current are linear in a switching period , it is possible to calculate the commutation instants as the intersection of such trajectories:

A more detailed description of this algorithm can be found in [10, 11].