Electric Vehicle Fast-Charging Station Unified Modeling and Stability Analysis in the dq Frame

The electric vehicle fast-charging station is an important guarantee for the popularity of electric vehicle. As the fast-charging piles are voltage source converters, stability issues will occur in the grid-connected fast-charging station. Since the dynamic input admittance of the fast-charging pile and the dynamic output impedance play an important role in the interaction system stability, the station and grid interaction system is regarded as load-side and source-side sub-systems to build the dynamic impedance model. The dynamic input admittance in matrix form is derived from the fast-charging pile current control loop considering the influence of the LC filter. Similarly, the dynamic output impedance can be obtained similarly by considering the regional power grid capacity, transformer capacity, and feed line length. On this basis, a modified forbidden region-based stability criterion is used for the fast-charging station stability analysis. The frequency-domain case studies and time-domain simulations are presented next to show the influence of factors from both the power grid side and fast-charging pile side. The simulation results validated the effectiveness of the dq frame impedance model and the stability analysis method.


Introduction
In recent years air pollution and haze have caused health issues around all the world.It is now becoming a more and more serious social problem that calls for less fossil fuels and more renewable energy.As a large consumer of fossil fuels, the transportation sector has emitted about 25% of the greenhouse gases produced by energy-related sectors, so the automobile's energy structure mode change is imperative [1].The demand of green transportation is growing rapidly due to these problems, and electric vehicles (EVs) are widely regarded as a proper solution for this problem.Some countries and regions, such as China and Europe, will start researching the exit schedule of traditional-fuel vehicles [2,3].This will contribute to the enhancement of the development and deployment of electric vehicles and lead to more and more EVs entering the market [4].
A large number of EVs require a large charging infrastructure to supply the electricity need.Generally, the EV power battery capacity varies from about 20 kWh to more than 100 kWh.According to the SAE J1772 standard [5], there are three charging standards for EVs that have been developed [6][7][8]: Level 1: The charging equipment uses single-phase AC power of 120 V voltage and 12 to 16 A current.It takes more than 10 h for a level 1 charger to charge a 20 kWh EV battery from 0% state of charge (SOC) to 100% SOC.
Level 2: The charging equipment uses 208-240 V single-phase AC power to provide up to 19.2 kW of charging power, needing at least 1 h for charging a 20 kWh EV battery to full SOC.
Level 3: This also known as fast-charging, the voltage used in this system is rated to 380 V in the AC side, 600 V in the DC side, and the maximum current is up to 400 A to provide a maximum 240 kW charge power.A 100 kWh EV battery can be fully charged within less than half an hour.
The common slow-charging piles are usually located in residential and office or shopping mall building parking lots, which cannot satisfy the sudden charging requirement due to the limitations of time and space.As a supplement to the slow charging, the EV fast-charging station is a necessary supplement and vital components for public acceptance of electric vehicles because it can greatly improve the ease of use for the EV owners.It makes the electricity energy refill as easy as refueling gasoline.
The EV fast-charging station generally consists of a distribution transformer and multiple fast-charging piles that are three-phase AC-DC voltage source converters (VSC) to generate the DC voltage for the EV battery.Since the fast-charging piles can provide a maximum charging power more than 200 kW, suppose that a EV fast-charging station contains 10 charging piles, when 10 are EVs charging simultaneously, the load can be up to 2000 kW, which is about the same as a residential or office building.Large numbers of EVs connecting to the urban power grid through EV fast-charging stations will change the dynamic characteristics of the power system since the fast-charging piles are controlled power electronic devices.EV fast-charging stations, together with the VSC-connected renewable generators, e.g., wind power and solar power, lead the power grid to power electronic features and make the AC grid relatively weak [9,10].Under such a situation, the interaction between the VSCs and the AC power grid will lead to complex oscillation problems.Such oscillation instability issues appeared in railway traction and supply systems where pantograph rising has been reported, a swing of the peak value of the grid side voltage appears on the sinusoidal wave as a low-frequency oscillation (LFO) [11].In Texas and the Hebei province of China, interaction between wind power plants and the power grid lead to a sub-synchronous oscillation that will influence the stability of the power plant [12,13].Sub-synchronous oscillation between the wind power plants even led to tripping of power plant units in Xinjiang.These instability issues are known to be caused by the power electronic interface loads and generators, such as electric vehicles, wind power plants, solar power plants, and electrified railway interaction with the power grid and interaction between these loads and generation.
Compared with the traditional LFO, the mechanism of interaction among VSCs and the AC grid is more complex and has not been fully understood.These kinds of power electronic interfaced oscillation and instability issues are generally analyzed from the viewpoint of small disturbances by the method of eigenvalues based on the state space model [14] and the impedance method based on frequency theory [15][16][17].The eigenvalue method is difficult to apply in the large-scale system with the high penetration of power electronic interfaced loads and generators due to its requirement of detailed models and parameters of the converters and the power system.Compared with the eigenvalue method, the impedance-based method regards the interaction system as two sub-systems, comprised of the load side converter and the source side AC grid at the point of common coupling (PCC).Only the output and input characteristics at the PCC of both sub-systems are concerned in this method to analyze the stability of the interaction system.In recent years the impedance method has been adopted in the stability analysis for grid-connected VSCs [18,19].Jian has proposed an impedance-based method in [19] by using the Nyquist criterion to analyze the stability of grid connected inverters.However, this method is developed for the single input and single output (SISO) system, and it is not suitable for stability analysis of the fast-charging station.
Generally, the three-phase AC-DC fast-charging piles are controlled in the dq frame.Thus, this paper builds the mathematical model of the three-phase fast-charging pile under the dq frame by using the transformation relationship between stationary frame and rotate frame.The fast-charging pile is considered as a multi-input multi-output (MIMO) system under the dq frame.In this mathematical model, the input variable and output variable are both two-dimensional vectors, the output impedance and input admittance are also a two-dimensional square matrix.The stability of the MIMO system can be analyzed based on the impedance-based method.That is based on the output impedance matrix of the grid side and the input admittance matrix of the fast-charging pile side.
In the dq frame dynamic model, the traditional impedance-based method cannot be used for stability analysis.For the MIMO system with the asymmetric impedance/admittance matrix, the generalized Nyquist criterion proposed by MacFarlane and Postlethwaite is widely used for stability assessment [20].Except the generalized Nyquist criterion, some other methods are proposed for the AC cascade system, such as singular value criterion [21,22], D channel criterion [23][24][25], and norm criterion [26].Compared with traditional methods, such as the Bode diagram, Nyquist criterion, dominate pole, etc., these methods are simpler in calculation.Singular value criterion and norm criterion have some conservatism in stability analysis, and the D channel criterion is obtained from experimental results of high power factor cascade systems, and it is not a sufficient condition, theoretically.The application of D channel criterion is limited since its validity in other AC cascade system cannot be guaranteed.A forbidden region-based stability analysis method in the dq frame that can be classified as a D channel criterion is proposed in [27].Since the high power factor condition cannot be met when instability occurs, the D channel criterion cannot be used for the EV fast-charging station stability analysis.Based on these methods, a novel stability criterion based on a different forbidden region from [28] is adopted for stability analysis of an EV fast-charging station in this paper.
This paper develops the dynamic impedance model of the EV fast-charging system in the dq frame at the point of common coupling, and stability of the fast-charging station is analyzed by using the novel forbidden region-based criterion.The rest of the paper is organized as follows: the detailed description of the EV fast-charging system, including the structure, connection, and control diagram is proposed in Section 2. Then the dynamic model of the EV fast-charging station and grid interaction system is presented in Section 3. In this part, the control system topology of the three-phase voltage source converter is presented, the dynamic model includes the matrix form input admittance of the fast-charging pile and the matrix output impedance of the AC grid side.The input admittance of the fast-charging pile is derived from the current control loop by using the frame transformation, and the dynamic model of the grid side is obtained through the output impedance of the grid system, including the power source, transformer, and distribution line, is presented.Frequency domain stability analysis is presented in Section 4. Furthermore, the time domain simulations are presented in Section 5 to validate the model proposed in this paper.

EV Fast-Charging System Description
Generally, the EV fast-charging station is supplied by the power grid through a distribution transformer via overhead lines or underground cables.The single-line electric diagram of the EV fast-charging station connected to the power grid through a substation transformer is shown in Figure 1.The source represents the power grid and one transformer is included in the substation, the AC bus represents the high voltage and low voltage of the substation, respectively.The fast-charging pile is the connection between the EV and the grid as Figure 1 shows.

Control System Introduction of the Fast-Charging Pile
The EV fast-charging pile in in the station is a three-phase AC/DC voltage source converter.The electrical tropology of the fast-charging pile is shown in Figure 2.  The LC-type filter is used to reduce harmonics and the DC link, that is, the EV battery side is represented by a constant DC voltage.The fast-charging pile is controlled by the voltage and current double loop control method based on the phase-locked loop (PLL).The grid-side voltage feedforward is included in the current control loop.The grid-side voltage at the common coupling point (PCC) is measured through the PLL to obtain the voltage angle.
For simplicity, the symbols used in this paper are clarified here first.The complex space vectors are represented by italic capital letters with subscript dq, e.g., I dq = I d + jI q or V dq = V d + jV q , here dq means the vectors are in the dq reference framework since the fast-charging pile is controlled in this reference frame.The corresponding real space vectors are represented by bold letters, such as I dq and T and V dq = V d , V q ] T ; this form represents the coordinates of the vectors in the corresponding frame.As the current and voltage in dq frame are two-dimensional, the admittance and impedance of the system in the dq frame are a 2 × 2 matrix, and they are denoted as Y 2m dq and Z 2m dq .In this figure, G c,dq (s) is the current loop PI controller in the dq frame, it is given by:

Current Control Loop Dynamic Modeling
G del,dq (s) depicts the effect of the digital computation delay and the PWM link delay, the computation delay is generally a switching period T s , and the PWM delay is generally half a switching period [29], thus: According to the first-order Padé approximation, G del,dq (s) is approximated by: Y V in i,s (s) and Y V g i,s (s) represent the influence of the LC filter and the grid-side impedance in the stationary frame.The subscript s means stationary frame; in other word, it is the αβ stationary frame.They can be derived from the following equations: where According to the frequency domain translation between the stationary frame to the rotation frame proposed in [30,31], as shown in Equation (8): As mentioned in [30], Y(S) is regarded as a combination of symbol s, which is actually a differential operator.Y V in i (s) and Y V g i (s) in dq frame can be obtained as follows: where θ is the angle between the stationary frame and the rotational frame, ω is the rated angular frequency of grid, and ω = dθ dt .The transfer function of the current control system can be expressed as follows: Here G cl,dq (s) and Y cl,dq (s) are the closed loop gain of the current loop and the dq frame input admittance of the fast-charging pile respectively, as shown in Equation (11).
T dq (s) is the loop gain of the current control loop, it is depicted as follows: Y o,dq (s) depicts the influence of grid side voltage feedforward and the LC filter, To obtain the model of the current control loop in matrix form in dq frame, all the variables should be transferred to matrix form by using the following transfer relationship: This transfer relationship is explained in Appendix A. According to the above transfer relationship, the current loop controller and the effect of the digital computation delay and the PWM link delay can be changed to matrix form as follows: The influence of the LC filter and the grid side impedance in the stationary frame Y V in i,dq (s) and Y V g i,dq (s) are transfer to matrix form by the following transformation: Then the transfer function of the current control system can be expressed as follows: Here G 2m cl,dq (s) and Y 2m cl,dq (s) are the closed loop gain of the current loop and the input admittance of the fast-charging pile in dq frame respectively, as shown in Equation (20): T 2m dq (s) is the loop gain of the current control loop, it is depicted as the following equation: Y 2m o,dq (s) depicts the influence of grid side voltage feedforward and the LC filter: I 2m is the 2 × 2 unit matrix, as follows: Thus, the dq frame dynamic function can be depicted as follows: where Y 2m cl,dq (s) is the input admittance matrix of the charging pile in dq frame.

Dynamic Model of the AC Power Grid
According to [19], the power grid side can be represented by the Thevenin equivalent model with an ideal voltage source and the output impedance.The equivalent impedance of the source side consists of the regional grid impedance, substation transformer impedance, and the charging station transformer impedance.

Regional Power Grid Impedance
The three-phase short circuit capacity of the regional 220 kV power grid and the distribution substation at the PCC is S SC , and the ratio of the inductance reactance X and resistance R is k X/R , thus, X and R can be obtained with: where U B is the base voltage.

Transformer Impedance
In the grid-side subsystem, there are two transformers, one is the 110 kV distribution transformer, the other one is the 10 kV distribution transformer.The rated capacity of the transformer is S T , the short circuit voltage percentage is U s %, the short circuit loss is P s , and the impedance of the transformer reduced to the higher voltage side can be obtained by:

Distribution Line Impedance
The impedance of the distribution line is determined by the line type and line length, the resistivity and reactance rate are ρ r and ρ x respectively.Then the impedance of the line is: where L is the length of the distribution line.

Equivalent Impedance of the Source Subsystem
The Thevenin equivalent model of source side subsystem is shown in Figure 5, the impedance can be obtained by: Furthermore, the inductance is: Thus, the impedance of the source side subsystem is Z s = R s + sL s .According to the control topology in Figure 5, the source side subsystem is consisting of an ideal AC source and a reactance, the linearized dynamic state equation of the reactance is: Thus, the output admittance of the grid side in stationary frame is: According to the frequency domain translation between the stationary frame to the rotation frame, the grid side output admittance is: The matrix form of the grid side output admittance is: (sL g +r g ) 2 +(ωL g ) 2 • sL g + r g ωL g −ωL g sL g + r g (33) Thus, the output impedance of the grid side is:

Stability Analysis Method for the MIMO System
Since both the input admittance matrix of the fast-charging pile and the output impedance matrix of the AC grid-side in the dq frame are derived in the previous section, the stability can be analyzed by using proper criteria.As it is known, the mathematical model of the fast-charging system in the dq frame is a MIMO system.The traditional stability analysis method is the Bode diagram and the generalized Nyquist criterion, but compared with the later-proposed singular value criterion, norm criterion, it is more complicated.However, the singular value criterion and norm criterion used for stability analysis are quite conservative.Thus, the forbidden region-based criterion proposed in [28] is used in this paper.

Generalized Nyquis, Singular Value, and Norm Criterion
The Nyquist criterion is extended for the scalar transfer functions by MacFarlane and Postlethwaite in 1980.It is developed to a generalized Nyquist criterion to address the stability by a matrix transfer function.The generalized Nyquist criterion for the system shown in Figure 6 can be stated as follows: The transfer function of this system is: where G(s)H(s) is defined as the return ratio of the system.
The assumption is made that the multivariable feedback system shown in Figure 6 have no open loop uncontrollable or unobservable modes whose corresponding characteristic frequencies lie in the right-half plane.Under this assumption, the system is stable if and only if the net sum of anti-clockwise encirclements of the critical point (−1,j0) by the set of characteristic loci of the return ratio L(s) = G(s)H(s) is equal to the total number of right-half plane poles of G(s) and H(s).For this paper, the fast-charging station and power grid are stable, respectively, so the number of right-half panel poles for both the grid-side output impedance Z 2m s,dq and the fast-charging pile-side input admittance Y 2m s,dq is 0, that is, the number of right-half panel poles for the return ratio is 0. Thus, the interaction system is stable if the number of anti-clockwise encirclements of the point (−1,j0) by all eigenvalues of the return ratio is 0. Hence the eigenvalues of the return ratio should be located within the region of a unit circle as shown in Figure 7.The region out of the unit circle is regard as the forbidden region to guarantee the stability of the system.
In the fast-charging station system, the output impedance of the grid side Z 2m go,dq (s) is written as: The fast-charging pile-side input admittance is written as: The return ratio is defined as: On the basis of generalized Nyquist criterion, singular value criterion and norm criterion for cascade system in dq frame are developed.
According to the singular value criterion, the cascade system is stable if the product of the maximum singular value of the output impedance matrix of the grid-side system and the maximum singular value of the input admittance matrix of the converter-side is less than 1 in the whole frequency domain.In mathematics it can be expressed as follows: where σ is the maximum module of the singular value for a matrix.As Equation ( 39) is satisfied, the eigenvalue of the return ratio L 2m dq (s) is surely located in the unit circle due to the module of eigenvalues being in the range determined by the module of the maximum and minimum singular values of a matrix.The mathematical explanation is proposed in Appendix B.
The norm criterion is easy for stability analysis of the cascade system due to its simplification, a modified norm criterion for the AC cascade system is proposed in [26] as follows: where the G norm and sum norm are defined as:

Modified Forbodden Region-Based Criterion
The modified forbidden region-based criterion focuses on reducing the conservatism of the previously-proposed criterion.As shown in Figure 7, the forbidden region of the eigenvalue of L 2m dq (s) is the outside region of the unit circle.In this modified criterion, the forbidden region is reduced to the left half panel of the critical point (−1,j0) as shown in Figure 8.In order to guarantee the eigenvalues are not located in the forbidden region shown in Figure 8, the Gershgorin circle cannot cross with the forbidden region, that is, the left most point of the Gershgorin circle should not less than −1, so the following conditions have to be met: According to [28], L 2m dq T (s) and L 2m dq (s) have the same eigenvalues, and all the eigenvalues of L 2m dq (s) are also located in the Gershgorin circle of L 2m dq T (s), that is: Thus, when any one of the Equations ( 43) and ( 44) is satisfied, the system is stable.The stable criterion for this system is: On the other hand, both Z 2m go,dq (s) and Y 2m cl,dq (s) have the same structure of: Thus, the return ratio L also have the same structure as mentioned above.It can be obtained: Thus the stability criterion can be simplified and expressed as follows: Re(L dd ) − L dq > −1 (48)

Parameters of the Fast-Charging Station System
The power grid-side subsystem in this paper generally contains four parts: the regional power grid, the substation distribution transformer, the distribution line (overhead line or underground cable), and the distribution transformer for the EV fast-charging station.The difference in any one of the four parts will affect the grid-side subsystem output impedance and then influence the fast-charging station stability.For the stability analysis, a standard condition is set as a contrast case.In this standard condition, the EV fast-charging station has five identical fast-charging piles.The related parameters of the main circuit are shown in Table 1.The control parameters of the current control loop are shown in Table 2 as follows: Table 2. Current loop control parameters.

Category Value
Proportional gain (K p ) 0.15 Integral gain (K i ) 1.5 The parameters of the power grid are shown in Table 3: Table 3. Power grid-side parameters.

Category Value
Regional Short Circuit Capacity(S) 300 MVA Ratio of inductance and reactance(k x/r ) 15 10 kV feeder cable(L feeder ) 2 km The parameters of the transformers are shown in Table 4:

Power Grid Side Subsystem Influence Factors
The relationship between the impedance of the grid side subsystem and the stability of the fast-charging station is investigated by regarding the four parts that constitute the grid-side subsystem as influence factors.Only one part is changed while the other parts are kept constant at each time to investigate the corresponding influence.The magnitude-frequency diagrams are shown in Figure 10.
Figure 10a shows the influence while the power grid short circuit capacity varies through 100 MVA-300 MVA-500 MVA-infinite.Figure 10b shows the influence while the 110 kV transformer capacity varies through 6.3 MVA-10 MVA-16 MVA-20 MVA. Figure 10c shows the influence while the 10 kV feeder length varies through 0.5 km-1 km-1.5 km-2 km. Figure 10d shows the influence, while the 10 kV transformer capacity varies through 0.5 MVA-0.63MVA-0.8MVA-1 MVA.According to the magnitude-frequency diagrams, it can be obtained that when any one factor is changed, the corresponding magnitude of the stability criterion changed.As shown in Figure 10a, along with the increasing of regional power grid short circuit capacity, the magnitude margin of the system is getting larger, which means the interaction system is getting more stable.Similarly, as shown in Figure 10b-d, along with the increasing of 110 kV and 10 kV distribution transformer capacity and decreasing of the 10 kV feeder length, the magnitude margin of the system is also getting larger.The law of change summarized from the figure is shown in the following Table 5.
The symbol ↑ in the table denotes the increase of influence factor variation, system stability, and grid-side output impedance.Similarly, the symbol ↓ denotes the decrease of these indices.
As shown in the table, no matter the increase of the power grid short circuit capacity, 110 kV transformer capacity and 10 kV transformer capacity, or the decrease of 10 kV feeder length will all lead to the improvement of the fast-charging system stability.Thus, it can be concluded that when the output impedance of the grid-side subsystem decreases, the system stability will be enhanced.On the other hand, when the grid-side output impedance decreases, it is generally regarded that the grid side is getting stronger.

Fast-Charging Pile-Side Influence Factors
The input admittance of the fast-charging pile-side is derived from the current control loop.According to the input admittance model, the current controller parameters and the filter parameters all have influence on the input admittance.Thus, this part investigates the influence of these parameters on the fast-charging station stability.The corresponding magnitude-frequency diagrams are shown in Figure 11.
Figure 11a shows the influence while the filter inductor varies through 1 mH-3 mH-5 mH-7 mH. Figure 11b shows the influence while the filter conductor varies through 50 µF-75 µF-100 µF-125 µF. Figure 11c shows the influence while the proportional gain varies through 0.1-0.15-0.2-0.5. Figure 11d shows the influence while the integral gain varies through 1-1.5-2-3.According to the magnitude-frequency diagrams, it can be obtained that when any one factor is changed, the corresponding magnitude of the stability criterion changed.The magnitude-frequency diagrams shown in Figure 11a,b depict the influence of filter parameters on the interaction system stability.According to Figure 11a, the magnitude margin increases along with the increasing of the filter inductance, while in Figure 11b, it is depicted that the stability margin decreases along with the increase in the filter conductance.The influence of the PI controller parameters are shown in Figure 11c,d, along with the increasing of the proportional factor and the integral factor, the stability margin of the interaction system decreases.The law of change summarized from the figure is shown in the following Table 6.
The symbols ↑ and ↓ denote the same meaning as in Table 5.
Along with the increase of the filter conductance, the proportional factor, and the integral factor, the stability margin of the fast-charging system decreases.The stability margin increases while the filter inductance increases.However, the charging pile-side input admittance variation tendency is difficult to determine since the expression of the input admittance is rather complex.The magnitude-frequency diagrams of the input admittance of the fast-charging system in these situations are shown in Figure 12.According to the diagrams shown in Figure 12, the influence of the filter parameters and current loop control parameters on the input admittance is uncertain.In some frequency ranges, the relationship is a positive correlation, while in other frequency ranges, the relationship is a negative correlation.The input admittance derived from the current control loop is actually a complex transfer function, and it has different expression with the common impedance that is expressed by r + s•L, so the relationship is uncertain.

Time Domain Simulation
In order to verify the modeling and stability analysis method, time domain simulation is used in this section.In the simulation, the equal output impedance of the grid-side subsystem and the LC filter inductor are changed to investigate their influence on the PCC voltage.The voltage waveforms of the PCC while the output impedance of the grid-side subsystem changes are shown in Figure 13.According to the waveform shown in Figure 13, as the output impedance increased to 10 times the original value, obvious distortion appears in the voltage waveform.From this point, the fast-charging system is becoming unstable when the output impedance increases.For more detailed information, the Fourier transformation is used to analysis the PCC voltage.The analysis result is shown in Figure 14.As shown in the spectrum, the harmonics increase significantly as the grid-side output impedance increases.Since the output impedance of the grid-side subsystem increases, the grid side subsystem is relatively weak; this change leads to instability.
The voltage waveforms of PCC while the LC filter inductance changes are shown in Figure 15.As the voltage waveform in the above picture shows, when the LC filter inductance decreases to 20% of the original value, the voltage waveform has serious distortion.When the LC filter inductance increases to 10 times the original value, the voltage waveform is very smooth.The time domain simulation depicts that the fast-charging system becomes stable along with the increase of the LC filter inductance.
For a more detailed and quantitative analysis of the system stability, Fourier transformation is used to analyze the PCC voltage.The analysis result is shown in Figure 16.
According to the spectrum shown in Figure 16, the harmonics decrease significantly as the fast-charging pile LC filter inductance increases.Since the fast-charging pile LC filter inductance increases, harmonics are filtered out by the inductor, so the system is more stable and the PCC voltage waveform is smooth.

Conclusions
This paper analyzes the electric vehicle fast-charging station and power grid interaction system stability by using the impedance-based modeling method.In this method, the interaction system is separated into two parts: the source-side subsystem and the load-side subsystem.The source side subsystem is represented by the Thevenin equivalent model and the load side subsystem is represented by the input admittance.
The dynamic modeling approach of the electric vehicle fast-charging station in dq frame is proposed.The input admittance of the fast-charging station in matrix form of dq frame is derived from the current control loop by using the complex vector and complex transfer function.On the basis of the dynamic model of the fast-charging system, the forbidden region-based stability criterion is adopted for stability analysis.Some conclusions can be drawn as follows: (1) The influence of the source-side and load-side parameters are analyzed by using the forbidden region-based stability criterion.Time domain simulation is provided to prove the proposed model and stability analysis method are reasonable.
(2) Parameters of both source-side and load-side can affect the stability of the interaction system.Different parameters have different influence tendencies.

Figure 1 .
Figure 1.Single-line sketch map of fast-charging station connection to the power grid.

Figure 3 .
Figure 3. Simplified block diagram of the three-phase fast-charging pile.

Figure 4
Figure 4 depicts the current control system block diagram with the PLL link and grid-side voltage feedforward, and the influence of grid-side voltage on the grid current is considered.

Figure 4 .
Figure 4. Block diagram of the inner current control loop.

Figure 6 .
Figure 6.Matrix transfer function feedback system.

Figure 7 .
Figure 7. Distribution of eigenvalues of the L dq determined by the generalized Nyquist criterion.

Figure 10 .
Figure 10.The magnitude-frequency of stability criterion under different situations of the grid-side subsystem.(a) Regional power grid with different short circuit capacity; (b) 110 kV distribution transformer with different capacity; (c) 10 kV feeder with different length; and (d) 10 kV charging station transformer with different capacity.

Figure 11 .
Figure 11.The magnitude-frequency of stability criterion under different situation of the fast-charging pile-side subsystem.(a) Different filter inductance of the fast-charging pile; (b) different filter capacitance of the fast-charging pile; (c) different proportional factor of the current control loop; and (d) different integral factor of the current control loop.

Figure 12 .
Figure 12.The admittance-frequency diagram of fast-charging pile-side input admittance.(a) Input admittance of different inductors; (b) input admittance of different capacitors; (c) input admittance of different proportion factors; and (d) input admittance of different integral factors.

Figure 13 .
Figure 13.The PCC voltage waveform of different grid-side output impedances.(a) Grid-side output impedance with the original value; and (b) the grid-side output impedance with 10 times the original value.

Figure 14 .
Figure 14.The PCC voltage Fourier analysis for different grid side output impedance.(a) FFT for the original output impedance; and (b) FFT for 10 times the output impedance.

Figure 15 .
Figure 15.The PCC voltage waveform of different LC filter inductance values.(a) Fast-charging pile LC filter inductance with 20% the original value; and (b) the fast-charging pile LC filter inductance with 10 times the original value.

Figure 16 .
Figure 16.The PCC voltage Fourier analysis of different LC filter inductor values.(a) FFT for 20% the LC filter inductance; and (b) FFT for 10 times the LC filter inductance.

Table 1 .
Fast-charging pile main circuit parameters.

Table 4 .
Parameters of the distribution transformers.

Table 5 .
Fast-charging pile main circuit parameters.

Table 6 .
Fast-charging pile main circuit parameters.