Suppression Research Regarding Low-Frequency Oscillation in the Vehicle-Grid Coupling System Using Model-Based Predictive Current Control

Abstract: Recently, low-frequency oscillation (LFO) has occurred many times in high-speed railways and has led to traction blockades. Some of the literature has found that the stability of the vehicle-grid coupling system could be improved by optimizing the control strategy of the traction line-side converter (LSC) to some extent. In this paper, a model-based predictive current control (MBPCC) approach based on continuous control set in the dq reference frame for the traction LSC for electric multiple units (EMUs) is proposed. First, the mathematical predictive model of one traction LSC is deduced by discretizing the state equation on the alternating current (AC) side. Then, the optimal control variables are calculated by solving the performance function, which involves the difference between the predicted and reference value of the current, as well as the variations of the control voltage. Finally, combined with bipolar sinusoidal pulse width modulation (SPWM), the whole control algorithm based on MBPCC is formed. The simulation models of EMUs’ dual traction LSCs are built in MATLAB/SIMULINK to verify the superior dynamic and static performance, by comparing them with traditional transient direct current control (TDCC). A whole dSPACE semi-physical platform is established to demonstrate the feasibility and effectiveness of MBPCC in real applications. In addition, the simulations of multi-EMUs accessed in the vehicle-grid coupling system are carried out to verify the suppressing effect on LFO. Finally, to find the impact of external parameters (the equivalent leakage inductance of vehicle transformer, the distance to the power supply, and load resistance) on MBPCC’s performance, the sensitivity analysis of these parameters is performed. Results indicate that these three parameters have a tiny impact on the proposed method but a significant influence on the performance of TDCC. Both oscillation pattern and oscillation peak under TDCC can be easily influenced when these parameters change.


Introduction
With the rapid development of high-speed railway, alternating current (AC)-direct current (DC)-AC drive electric multiple units (EMUs) and electric locomotives are increasingly put into operation. Meanwhile, low-frequency oscillation (LFO) accidents have happened in many countries, such as Norway, Germany, Switzerland, the United States, and France [1][2][3][4]. Since 2008, the phenomenon has frequently occurred in China's high-speed railway depots. The LFO in railway is characterized by the amplitude fluctuation of grid-side voltage, current, and DC-side voltage, and happens when the multiple vehicles are concentrated in one power supply district and get power from a traction network [5]. With the larger voltage oscillation peak of the traction network, The paper is organized as follows. In Section 2, the mathematical model of MBPCC based on the topology structure of one traction LSC of the CRH3 EMUs is deduced. In Section 3, the simulation model of EMUs' dual traction LSCs is built, and a comprehensive simulation verification compared with TDCC is performed in order to verify the superiority of the proposed control strategy from several aspects. In Section 4, a whole dSPACE semi-physical platform is established to certify the control strategy's feasibility in real applications. In Section 5, based on a reduced-order model of the traction network, the simulation model of seven vehicles accessed in the traction network is constructed to verify the suppression effect of LFO, and the sensibility analysis of parameters under the condition of seven EMUs accessed in the traction network is discussed, in order to further explore the impact on LFO occurrence and the performance of traction LSCs when system parameters change. The study idea is shown in Figure 1. model of EMUs' dual traction LSCs is built, and a comprehensive simulation verification compared with TDCC is performed in order to verify the superiority of the proposed control strategy from several aspects. In Section 4, a whole dSPACE semi-physical platform is established to certify the control strategy's feasibility in real applications. In Section 5, based on a reduced-order model of the traction network, the simulation model of seven vehicles accessed in the traction network is constructed to verify the suppression effect of LFO, and the sensibility analysis of parameters under the condition of seven EMUs accessed in the traction network is discussed, in order to further explore the impact on LFO occurrence and the performance of traction LSCs when system parameters change.
The study idea is shown in Figure 1.

Mathematical Model of One Traction Line-Side Converter
The vehicle-grid coupling system of high-speed railway is composed of the traction network and the EMUs traction drive system, as shown in Figure 2. In this paper, a CRH3 EMU is set as the study object; thus, the topology of the traction drive unit of EMUs is a single-phase two-level structure. When the LFO happens, the EMUs just start up and only supply power for the auxiliary facility by Verification of one traction drive unit

Section 5
The two-step prediction prediction

Mathematical Model of One Traction Line-Side Converter
The vehicle-grid coupling system of high-speed railway is composed of the traction network and the EMUs traction drive system, as shown in Figure 2. In this paper, a CRH3 EMU is set as the study object; thus, the topology of the traction drive unit of EMUs is a single-phase two-level structure. When the LFO happens, the EMUs just start up and only supply power for the auxiliary facility by DC-side voltage. Therefore, the inverter and motor can be regarded as a pure resistance, and only part of the rectifier is involved in the equivalent circuit of one traction LSC. DC-side voltage. Therefore, the inverter and motor can be regarded as a pure resistance, and only part of the rectifier is involved in the equivalent circuit of one traction LSC.
In the equivalent circuit of EMUs' traction LSC, N L and N R denote the leakage inductance and leakage resistance of the secondary winding of the vehicle transformer in EMUs, respectively. Four IGBTs ( 1 S , 2 S , 3 S , and 4 S ) are used to construct the two-level topology, d C is the DC-side capacitor, inductance 2 L and capacitor 2 C compose the second-order filter circuit, and the inverter and motor are replaced by a pure resistive load load R .
According to the Kirchhoff voltage law, the relationship between the voltage and current in the AC-side is listed below:    In the equivalent circuit of EMUs' traction LSC, L N and R N denote the leakage inductance and leakage resistance of the secondary winding of the vehicle transformer in EMUs, respectively. Four IGBTs (S 1 , S 2 , S 3 , and S 4 ) are used to construct the two-level topology, C d is the DC-side capacitor, inductance L 2 and capacitor C 2 compose the second-order filter circuit, and the inverter and motor are replaced by a pure resistive load R load .
According to the Kirchhoff voltage law, the relationship between the voltage and current in the AC-side is listed below: where u N , i N , and u ab represent the line voltage, line current, and input voltage of the converter, respectively. If the harmonics are neglected, u N , i N , and u ab in the dq reference frame are defined as follows:      u N = u Nd sin(ωt) + u Nq cos(ωt) i N = i Nd sin(ωt) + i Nq cos(ωt) u ab = u abd sin(ωt) + u abq cos(ωt) (2) www.manaraa.com where u Nd , i Nd , and u abd are the d-axis components of u N , i N and u ab . u Nq , i Nq , and u abq are the q-axis components of u N , i N , and u ab , respectively. Substituting Equation (2) into Equation (1), the mathematical model of one traction LSC can be described as: Applying the first-order discrete approximation to the mathematical model expressed in Equation (3), a discrete dynamic model of one traction LSC in the dq reference frame is depicted by Equation (4): where T s is the sampling interval. i Nd (k) and i Nq (k) are the discrete values of i Nd and i Nq , respectively; i Nd (k + 1) and i Nq (k + 1) are the one-step predictive discrete values of i Nd and i Nq , respectively; and u abd (k) and u abq (k) are the discrete values of the control voltages at the k-th sampling instant. Thus, i Nd (k + 1) and i Nq (k + 1) can be predicted at the k-th sampling instant: The relationship between u ab (k) and u ab (k − 1) is satisfied in the following way: where ∆u abd (k) represents the variation between u abd (k) and u abd (k − 1), and ∆u abq (k) represents the variation between u abq (k) and u abq (k − 1). Substituting Equation (6) into Equation (5), then: Equation (7) is the prediction model. It should be noted that prediction line current values are decided by present currents, voltages, and the LSC's system parameters.

The Two-Step Prediction
In the digital control system, both the computation time delay and sampling time delay surely exist [28]. Because of the restriction of the hardware and the digital control algorithm, the control voltages calculated by the sampling and predictive currents have to be adopted at the (k + 1)-th sampling instant. Therefore, there is one sampling period delay (T d ) between the calculated control voltage by the controller and the real adopted control voltage, as shown in Figure 3. To eliminate the error caused by T d , a two-step current prediction is adopted.
www.manaraa.com voltage by the controller and the real adopted control voltage, as shown in Figure 3. To eliminate the error caused by d T , a two-step current prediction is adopted.
One-step prediction Two-step prediction When two-step prediction is adopted, Nq i k+ , respectively, cannot be omitted in Equation (9):

The Design of Performance Function
The key to achieving MPC is to obtain the most effective control quantities by solving the performance function optimally. In this paper, the optimal control voltage would be obtained.
The performance function is composed of the predicted current components at the (k + 2)-th sampling instant and the control voltage variations, with corresponding weighting coefficients [29]. The function is defined as follows: When two-step prediction is adopted, i Nd (k + 1) and i Nq (k + 1) are first predicted at the k-th sampling instant according to Equation (8), and then i Nd (k + 2) and i Nq (k + 2) are predicted according to Equation (9). Thus, the control voltage u abd (k + 1) and u abq (k + 1) can be calculated at the k-th sampling instant, and then adopted to the rectifier at the k+1-th sampling instant. The delay error is compensated by two-step prediction.
It is worth noting that i Nd (k + 1) and i Nq (k + 1) are also predicted through two-step prediction at the (k − 1)-th sampling instant. Thus, i Nd (k + 1) and i Nq (k + 1), which are predicted at the k-th sampling instant, are considered equal to the predicted values at the k−1-th sampling instant. Thus, the control voltage variations ∆u abd (k) and ∆u abq (k), which exist in the Equation (7), can be set to zero. The prediction Equation (7) can be rewritten as Equation (8). While i Nd (k + 2) and i Nq (k + 2) are not predicted before the k-th sampling instant, it cannot be affirmed that the prediction current values are equal to the k-th actual current values. The control voltages ∆u abd (k + 1) and ∆u abq (k + 1) that determine i Nd (k + 2) and i Nq (k + 2), respectively, cannot be omitted in Equation (9):

The Design of Performance Function
The key to achieving MPC is to obtain the most effective control quantities by solving the performance function optimally. In this paper, the optimal control voltage would be obtained.
The performance function is composed of the predicted current components at the (k + 2)-th sampling instant and the control voltage variations, with corresponding weighting coefficients [29]. The function is defined as follows: where α 1 , α 2 , β 1 , and β 2 represent the weighting coefficients of the line current and voltage variations. There are no analytical or numerical methods or control design theories to adjust these parameters; currently, they are determined based on empirical procedures. In [30], an approach based on an empirical procedure is presented to obtain suitable weighting factors. When more objectives are www.manaraa.com considered, the weighting coefficients are usually obtained using trial and error procedures and running time-consuming simulations [31]. Since the front two components in Equation (10) are both the current variables, weighting factors can be considered as the same, as are the back two components. Thus, if two of weighting coefficients are decided, the other two can also be set, which can be clearly seen from Table A1.
In order to make i Nd (k + 2) and i Nq (k + 2) track their references i Nd (k) * and i Nq (k) * , the variations of the control voltage need to be kept as small as possible. To do this, take the derivative of Equation (10) to find the extreme point, as shown in Equation (11): Substituting i Nd (k + 2) and i Nq (k + 2) of Equation (9) and w(k) of Equation (10) into Equation (11), the optimal control variables ∆u abd (k + 1) and ∆u abq (k + 1) can be derived as Equation (12). In addition, the output control voltages u abd (k) * and u abq (k) * can be depicted by Equation (13): The optimal control voltages that minimize the performance function are calculated by Equation (13), and then fed to a modulator stage to generate the PWM drive signal g at the k-th instant. An MBPCC controller of LSCs can be designed based on the deducing process of prediction model, the delay compensation, and the solution of performance function. Figure 4 shows the block diagram of the MBPCC. A PI controller is adopted in the voltage loop to regulate the DC-side voltage, and its output is set as the current d-axis reference i Nd (k) * . In addition, the q-axis current reference i Nq (k) * generally needs to equal zero to achieve the unit power factor.

Simulations of One Traction Drive Unit of Electric Multiple Units
To validate the feasibility and effectiveness of the proposed method, simulation verifications based on the model of EMUs' dual traction LSCs, namely the traction drive unit, were carried out from several aspects, by comparing with the control performance of traditional TDCC in MATLAB /SIMULINK, respectively. The simulation model of EMUs' dual traction LSCs was built as shown in Figure 5. The adjustable parameters in the TDCC and MBPCC controllers were regulated into the most appropriate values, as listed in Tables A1 and A2. When LFO happens in high-speed railways, the fluctuation of the traction network's voltage will disturb the performance of the EMUs' LSC. Currents and voltages of the controller cannot trace well with the references, and lead to greater deterioration of fluctuation. MBPCC is a current-controlled strategy constituted by a prediction model and the performance function. By making the performance function minimal, the optimal control variables can be calculated to ensure the prediction current www.manaraa.com value tracking the references. Therefore, it is possible to adopt MBPCC to optimize the performance of the converter controller and suppress the LFO.

Simulations of One Traction Drive Unit of Electric Multiple Units
To validate the feasibility and effectiveness of the proposed method, simulation verifications based on the model of EMUs' dual traction LSCs, namely the traction drive unit, were carried out from several aspects, by comparing with the control performance of traditional TDCC in MATLAB /SIMULINK, respectively. The simulation model of EMUs' dual traction LSCs was built as shown in Figure 5. The adjustable parameters in the TDCC and MBPCC controllers were regulated into the most appropriate values, as listed in Tables A1 and A2.

Simulations of One Traction Drive Unit of Electric Multiple Units
To validate the feasibility and effectiveness of the proposed method, simulation verifications based on the model of EMUs' dual traction LSCs, namely the traction drive unit, were carried out from several aspects, by comparing with the control performance of traditional TDCC in MATLAB /SIMULINK, respectively. The simulation model of EMUs' dual traction LSCs was built as shown in Figure 5. The adjustable parameters in the TDCC and MBPCC controllers were regulated into the most appropriate values, as listed in Tables A1 and A2.

Recifier1
Powergui The block diagram of TDCC is shown in Figure 6. The outer voltage control loop uses a PI controller to keep the DC-side voltage equal to its reference value, and the PI output provides the reference of input current. The proportional controller of the inner current control loop makes the input current track its reference value. The block diagram of TDCC is shown in Figure 6. The outer voltage control loop uses a PI controller to keep the DC-side voltage equal to its reference value, and the PI output provides the reference of input current. The proportional controller of the inner current control loop makes the input current track its reference value.

Start-Operation Process
The start-operation process of EMUs' dual traction LSCs is divided into three periods-namely, the pre-charge period, the uncontrolled rectifier period, and the nominal load period. In the simulation, the initial time of the three stages was set to 0 s, 0.2 s, and 0.4 s, respectively. To compare MBPCC with TDCC, the waveforms of N u and N i , as well as d u are depicted in Figure 7. The performance indexes of the d u of EMUs' dual traction LSCs are shown in Table 1. Figure 8 shows

Start-Operation Process
The start-operation process of EMUs' dual traction LSCs is divided into three periods-namely, the pre-charge period, the uncontrolled rectifier period, and the nominal load period. In the simulation, the initial time of the three stages was set to 0 s, 0.2 s, and 0.4 s, respectively. To compare MBPCC with www.manaraa.com  Table 1. Figure 8 shows the fast Fourier transform (FFT) analysis results of the i N with MBPCC and TDCC.

Start-Operation Process
The start-operation process of EMUs' dual traction LSCs is divided into three periods-namely, the pre-charge period, the uncontrolled rectifier period, and the nominal load period. In the simulation, the initial time of the three stages was set to 0 s, 0.2 s, and 0.4 s, respectively. To compare MBPCC with TDCC, the waveforms of N u and N i , as well as d u are depicted in Figure 7. The performance indexes of the d u of EMUs' dual traction LSCs are shown in Table 1. Figure Figure 7, the pre-charge and uncontrolled rectifier periods of TDCC are almost the same as MBPCC's. When the nominal load is accessed, the u d based on TDCC achieves the reference value after 0.4 s, with a large overshoot. However, the u d based on MBPCC achieves stability after 0.25 s, and the overshoot is about 100 V, far less than that of the PI control. The voltage fluctuation range of MBPCC is only ±10 V, while the voltage fluctuation based on TDCC is larger. In Figure 8, the total harmonic distortion (THD) of i N based on MBPCC is 4.76%, apparently lower than that of TDCC. In the control system based on MBPCC, there are more high-order harmonics around the odd switching frequency, which is considered as a drawback of the control strategy. Overall, MBPCC presents a better performance and dynamic response, due to the smaller overshoot, shorter adjustment time, and tinier voltage fluctuation.
www.manaraa.com and the overshoot is about 100 V, far less than that of the PI control. The voltage fluctuation range of MBPCC is only ±10 V, while the voltage fluctuation based on TDCC is larger. In Figure 8, the total harmonic distortion (THD) of N i based on MBPCC is 4.76%, apparently lower than that of TDCC.
In the control system based on MBPCC, there are more high-order harmonics around the odd switching frequency, which is considered as a drawback of the control strategy. Overall, MBPCC presents a better performance and dynamic response, due to the smaller overshoot, shorter adjustment time, and tinier voltage fluctuation.

Sudden-Load-Change Process
To further validate the dynamic response performance, the simulations under the sudden-loadchange condition were carried out. Figure 9 shows

Sudden-Load-Change Process
To further validate the dynamic response performance, the simulations under the sudden-load-change condition were carried out. Figure 9 shows the simulation waveforms of the dynamic response of the u d under the condition of the load changing suddenly, based on MBPCC and TDCC, respectively. The load changes at 1 s, from 10 Ω to 0.01 Ω. Based on MBPCC, the u d has gone through a period of voltage decline from 3000 V to about 2500 V, and returns to 3000 V after 0.5 s, with a ±20 V voltage fluctuation later. In contrast, the u d based on TDCC drops by nearly 800 V, and fluctuates in the range of ±50 V after becoming stable again. Therefore, it can be concluded that the MBPCC has a better capacity for resisting disturbance than the transient PI control, because of smaller voltage decline and fluctuation when the load changes suddenly.

Track Performance
At first, we measured  www.manaraa.com

Track Performance
At first, we measured i Nd * when the outer voltage loop plays a role in obtaining the reference current value. i Nd * was about equal to 830 A, and i Nq * was set to zero. Based on the measured values, we replaced the outer voltage loop by a step signal, to test the track performance while the reference current steps up or down similarly, a step signal replaces the constant module, in order to detect the track performance while the reference current i Nq * steps up or down.
As shown in Figure 10a, i Nd decreased with i Nd * , varying from 830 A to 600 A at 0.3 s, and i Nq returned to zero after a prompt downward fluctuation. In Figure 10b, i Nd increases, with i Nd * varying from 830 A to 1000 A, and i Nq rapidly returns the setting value i Nq * after an upward shock. Similar conditions can be seen in Figure 11, when i Nq * varies. In two cases,

Semi-Physical Test of Electric Multiple Units' Dual Traction Line-Side Converter
To reflect the real application condition, a whole dSPACE semi-physical experimental platform was established. The dSPACE semi-physical experimental platform included the dSPACE simulator, the physical control circuit chassis, a power supply, and an external PC, as shown in Figure 12. The dSPACE simulator was used to simulate the circuit topology of EMUs' dual traction LSCs, and connected with the PC through a network line, so as to import the circuit simulation model from the PC to the dSPACE simulator. The physical control circuit chassis was connected to the PC through the data cable, in order to achieve the control algorithm program import and debugging. The power supply was supplied to the physical control circuit chassis.

Semi-Physical Test of Electric Multiple Units' Dual Traction Line-Side Converter
To reflect the real application condition, a whole dSPACE semi-physical experimental platform was established. The dSPACE semi-physical experimental platform included the dSPACE simulator, the physical control circuit chassis, a power supply, and an external PC, as shown in Figure 12. The dSPACE simulator was used to simulate the circuit topology of EMUs' dual traction LSCs, and connected with the PC through a network line, so as to import the circuit simulation model from the PC to the dSPACE simulator. The physical control circuit chassis was connected to the PC through the data cable, in order to achieve the control algorithm program import and debugging. The power supply was supplied to the physical control circuit chassis.   In conclusion, MBPCC presents a better control performance for EMUs' dual traction LSCs, due to the smaller overshoot, shorter adjustment time, and tinier voltage fluctuation in the start-operation process, as well as its greater capacity for resisting disturbance and better track performance between the actual current and the reference current.

Semi-Physical Test of Electric Multiple Units' Dual Traction Line-Side Converter
To reflect the real application condition, a whole dSPACE semi-physical experimental platform was established. The dSPACE semi-physical experimental platform included the dSPACE simulator, the physical control circuit chassis, a power supply, and an external PC, as shown in Figure 12. The dSPACE simulator was used to simulate the circuit topology of EMUs' dual traction LSCs, and connected with the PC through a network line, so as to import the circuit simulation model from the PC to the dSPACE simulator. The physical control circuit chassis was connected to the PC through the data cable, in order to achieve the control algorithm program import and debugging. The power supply was supplied to the physical control circuit chassis.  The most important part of the dSPACE semi-physical experimental platform was the physical control circuit chassis, which was composed of the main five modules, as shown in right-hand section of Figure 12: (1) Switching power supply module: the switching power supply module has the function of a power supply for the entire physical control circuit chassis. (2) SMC module: the SMC module realizes the data transmission of the physical control circuit chassis and the outside. After control program is compiled on the computer, the SMC module achieves the connection with the computer through the data cable, and control strategy is imported into the physical control circuit chassis.  The most important part of the dSPACE semi-physical experimental platform was the physical control circuit chassis, which was composed of the main five modules, as shown in right-hand section of Figure 12: (1) Switching power supply module: the switching power supply module has the function of a power supply for the entire physical control circuit chassis. (2) SMC module: the SMC module realizes the data transmission of the physical control circuit chassis and the outside. After control program is compiled on the computer, the SMC module achieves the connection with the computer through the data cable, and control strategy is imported into the physical control circuit chassis. www.manaraa.com The experimental results on the dSPACE semi-physical platform are almost in accordance with the simulation results in MATLAB. In Figure 13a, the d u based on MPBCC achieves the reference value without the overshoot. As seen in Figure 13b, N i can remain in phase with N u , which means that the dual traction LSCs operate with a unified power factor. Compared with TDCC, as shown in Figure 14, N i based on MBPCC has the smaller harmonic distortion.

System Verification
The dSPACE semi-physical experiment when multi-EMUs are accessed in the power supply system is still a challenge. To verify the suppression effect of LFO and perform further sensibility analyses, the simulations of seven EMUs accessed in traction network were performed in MATLAB/SIMULINK. The autotransformer power supply system is mostly adopted in high-speed railways of China. The structure of a traction network is very complicated because of a mass of multiconductor transmission lines that are distributed, as well as the mutual coupling effect [5]. Therefore, considering the practical factors, such as the skin effect of lines and external disturbances, it is more reasonable to adopt a reduced-order method to model the traction network [5,32] than the Theveninequivalent method used in [33]. Avoiding a duplication of effort, the modeling process of the vehiclegrid coupling system is no longer described. The experimental results on the dSPACE semi-physical platform are almost in accordance with the simulation results in MATLAB. In Figure 13a, the u d based on MPBCC achieves the reference value without the overshoot. As seen in Figure 13b, i N can remain in phase with u N , which means that the dual traction LSCs operate with a unified power factor. Compared with TDCC, as shown in Figure 14, i N based on MBPCC has the smaller harmonic distortion.

System Verification
The dSPACE semi-physical experiment when multi-EMUs are accessed in the power supply system is still a challenge. To verify the suppression effect of LFO and perform further sensibility analyses, www.manaraa.com the simulations of seven EMUs accessed in traction network were performed in MATLAB/SIMULINK. The autotransformer power supply system is mostly adopted in high-speed railways of China. The structure of a traction network is very complicated because of a mass of multi-conductor transmission lines that are distributed, as well as the mutual coupling effect [5]. Therefore, considering the practical factors, such as the skin effect of lines and external disturbances, it is more reasonable to adopt a reduced-order method to model the traction network [5,32] than the Thevenin-equivalent method used in [33]. Avoiding a duplication of effort, the modeling process of the vehicle-grid coupling system is no longer described.

The Effect of Suppressing Low-Frequency Oscillation
Based on the TDCC controller, when seven vehicles are accessed in the reduced-order model of traction network, LFO occurs in both voltage and current, as shown in Figure 15. The u d fluctuates between 2600 V and 3600 V. It can also be found that the traction network voltage arrives at the wave crest while the grid-side current reaches the wave trough. The amplitude fluctuation of voltages and currents will influence the performance of the EMUs' LSC, and even lead to a traction blockade.

The Effect of Suppressing Low-Frequency Oscillation
Based on the TDCC controller, when seven vehicles are accessed in the reduced-order model of traction network, LFO occurs in both voltage and current, as shown in Figure 15. The d u fluctuates between 2600 V and 3600 V. It can also be found that the traction network voltage arrives at the wave crest while the grid-side current reaches the wave trough. The amplitude fluctuation of voltages and currents will influence the performance of the EMUs' LSC, and even lead to a traction blockade. By adopting the proposed MBPCC strategy, the system can still achieve stability when the number of multi-EMUs accessed in the traction network reaches seven. Figure 16  By adopting the proposed MBPCC strategy, the system can still achieve stability when the number of multi-EMUs accessed in the traction network reaches seven. Figure 16 shows the simulation results when seven EMUs are accessed. It can be observed that the LFO does not appear even though vehicle numbers have reached their critical value. In a stable state, the u d remains stable at about 3000 V, and the voltage deviation is ±10 V. The fluctuations of grid-side voltage and current are small.
www.manaraa.com By adopting the proposed MBPCC strategy, the system can still achieve stability when the number of multi-EMUs accessed in the traction network reaches seven. Figure 16 shows the simulation results when seven EMUs are accessed. It can be observed that the LFO does not appear even though vehicle numbers have reached their critical value. In a stable state, the d u remains stable at about 3000 V, and the voltage deviation is ±10 V. The fluctuations of grid-side voltage and current are small.

Analysis of System Parameters
The analysis of parameter sensitivity is necessary to find out their influence on the stability of the vehicle-grid coupling system. In this paper, in order to find out the effects of different external parameters on the LFO and the performance of traction LSC, three parameters (the load resistance  can be concluded that the load resistance has fewer effects on the proposed method than TDCC.

Analysis of System Parameters
The analysis of parameter sensitivity is necessary to find out their influence on the stability of the vehicle-grid coupling system. In this paper, in order to find out the effects of different external parameters on the LFO and the performance of traction LSC, three parameters (the load resistance R load , equivalent leakage inductance, L N and the distance to power supply D) are discussed when seven vehicles are accessed in the traction network. Figure 17 shows the waveforms of the u d under TDCC and MBPCC, when the load resistance R load is set as 20, 50, 75, or 100 Ω, respectively. As seen in Table 2 and Figure 17, when the load resistance R load varies, both the overshoot and adjustment time of the u d under MBPCC have only some tiny variations.

Analysis of System Parameters
The analysis of parameter sensitivity is necessary to find out their influence on the stability of the vehicle-grid coupling system. In this paper, in order to find out the effects of different external parameters on the LFO and the performance of traction LSC, three parameters (the load resistance is set as 20, 50, 75, or 100 Ω, respectively. As seen in Table 2 and Figure      However, LFO happens under TDCC when the load resistance is small. The smaller R load is, the severer the oscillation is. LFO would not occur once R load exceeds a value about 75 Ω. Therefore, it can be concluded that the load resistance has fewer effects on the proposed method than TDCC.

Equivalent Leakage Inductance L N
In Figure 18, when the equivalent leakage inductance of vehicle transformer L N varies from 0.001 H to 0.006 H, the system under MBPCC maintains an initially good dynamic performance, with little change in the u d 's overshoot and regulation time. Hence, MBPCC can keep good performance and be unaffected by the change of L N , while TDCC is very sensitive. When L N is set as 0.001 H, the LFO does not happen. With the increase of L N , oscillation appears, and the oscillation peak enlarges gradually.

Distance to Power Supply D
D represents the distance to the power supply in the traction network. The larger D is, the larger the equivalent impedance on the line side will be. As seen in Figure 19 and Table 2, with the change of D, the influence on MBPCC's control performance is slight, because the overshoot and adjustment time change only slightly more than before, and the LFO does not happen.
As for TDCC, the oscillation phenomenon can be observed when D was chosen as 10, 20, 30, and 40 km. The oscillation peak value is about 3500 V at first, and then gradually diminishes in the former three cases; when D is 40 km, the oscillation peak will be kept at 3500 V. The details are listed in Table  2. It was found that the oscillation peak varies from 3530 V to 3240 V within 1 s when D is 10 km, and from 3520 V to 3300 V when D is 30 km. The damping capability of LFO decreases with the increase of D under TDCC.
Contrasting results show that MBPCC can be almost unaffected by the three parameters, and LFO does not happen when these parameters change. However, the performance of TDCC is sensitive to the three parameters. Both the oscillation pattern and oscillation peak can be easily

Distance to Power Supply D
D represents the distance to the power supply in the traction network. The larger D is, the larger the equivalent impedance on the line side will be. As seen in Figure 19 and Table 2, with the change of D, the influence on MBPCC's control performance is slight, because the overshoot and adjustment time change only slightly more than before, and the LFO does not happen.
www.manaraa.com Contrasting results show that MBPCC can be almost unaffected by the three parameters, and LFO does not happen when these parameters change. However, the performance of TDCC is sensitive to the three parameters. Both the oscillation pattern and oscillation peak can be easily influenced when the parameters change.   As for TDCC, the oscillation phenomenon can be observed when D was chosen as 10, 20, 30, and 40 km. The oscillation peak value is about 3500 V at first, and then gradually diminishes in the former three cases; when D is 40 km, the oscillation peak will be kept at 3500 V. The details are listed in Table 2. It was found that the oscillation peak varies from 3530 V to 3240 V within 1 s when D is 10 km, and from 3520 V to 3300 V when D is 30 km. The damping capability of LFO decreases with the increase of D under TDCC.
Contrasting results show that MBPCC can be almost unaffected by the three parameters, and LFO does not happen when these parameters change. However, the performance of TDCC is sensitive to the three parameters. Both the oscillation pattern and oscillation peak can be easily influenced when the parameters change.

Conclusions
To optimize the control performance of EMUs' traction drive units, namely dual traction LSCs, and suppress LFO in vehicle-grid coupling system of high-speed railways, MBPCC is proposed in order to apply the traction LSC of CRH3 EMUs. After theoretical analysis, simulation verifications, and semi-physical verifications, the performance of MBPCC is demonstrated in comparison to TDCC, and the advantages of each listed below: (1) Simulation verifications in MATLAB of EMUs' dual traction LSCs based on MBPCC and TDCC were implemented from three aspects. The results prove that MBPCC can obtain better dynamic and static performances, such as a smaller THD, tinier overshoot in start-operation process, greater capacity for resisting disturbance under load changing suddenly, and a better track performance. (2) Semi-physical verifications in the dSPACE semi-physical experimental platform were realized.
The results certified the effectiveness of MBPCC and its superiority in real applications, when compared with TDCC. (3) When the multi-EMUs were assessed in the reduced-order model of traction network, the results showed that MBPCC can ensure the system's stability and suppress LFO more efficiently compared with TDCC. (4) The influences of different external parameters R load , L N , and D in the vehicle-grid coupling system under MBPCC and TDCC have been discussed in detail. It could be concluded that these three parameters have a tiny impact on MBPCC, while they greatly influence the performance of TDCC. Both the oscillation pattern and oscillation peak under TDCC can be easily influenced when parameters change.
The proposed method can be applied for the control of EMUs' traction LSCs, and provide a good effect on the suppression of LFO. Moreover, MBPCC is insensitive to system parameters, which provides greater possibilities for its application. There is some advanced work regarding MBPCC that can be studied in the future. In the aspect of physical verification, multi-EMUs accessed in the model of the traction network should be realized in order to increase the reliability for suppressing LFO. With respect to algorithm optimization, combining predictive control with the disturbance observer can further advance the performance and robustness of the system, which can deal with the modeling errors and uncertainties, disturbances, and sensor noise.
Author Contributions: The individual contribution of each co-author with regards to the reported research and writing of the paper is as follows. Z.L. and Y.W. conceived the idea, Y.W. performed experiments and data analysis, and all authors wrote the paper. All authors have read and approved the final manuscript.