Abstract
This paper studies the transient process of synchronous generator after threephase short circuit, deduces the expression of stator current in the transient process of synchronous generator, and verifies the theoretical conclusion through simulation. After the short circuit occurs, the electromagnetic transient process is complex. In the electromagnetic transient process, the stator current and rotor current of synchronous generator appear attenuation oscillation due to directaxis reactance, crossaxis reactance, transient reactance and subtransient reactance. Based on the basic principle of synchronous generator, this paper analyzes the analytical expressions of stator straight axis current, stator cross axis current and total current. After that, the simulation model is built in Matlab/Simulink, and the simulation results are consistent with the theoretical analysis results. This simulation method provides convenience for analyzing the transient process of synchronous generator and has good engineering application value.
Highlights
 The research method of the transient process of threephase shortcircuit fault of synchronous generator is proposed.
 Analytical expression and Matlab/Simulink were used to simulate the transient process, and the characteristics of the transient process were discussed.
 The results show that the method proposed in this paper can accurately simulate the transient process of threephase short circuit of synchronous generator, which lays a foundation for the parameter d
1. Introduction
Synchronous generator is the core equipment of modern power system. In the event of threephase shortcircuit fault, the transient process of synchronous generator has an important influence on the safe operation and protection of synchronous generator [1]. After the large disturbance of the power system, the synchronous generator experiences the electromagnetic transient process and the electromechanical transient process. This paper studies the electromagnetic transient process of the synchronous generator [2]. The variation of terminal voltage, stator current, excitation current, power angle and rotational speed of synchronous generator need to be studied in electromagnetic transient process. A synchronous generator model is built in MATLAB/Simulink simulation environment. It is important to study the protection device of generator that studies the voltage and current change of synchronous generator in threephase shortcircuit fault.
2. Threephase short circuit analysis of synchronous generator
When discussing the threephase short circuit of synchronous generator, the normal operation and fault operation state of synchronous generator can be superimposed and solved [3, 4]. This is equivalent to adding a threephase voltage equal to and reverse to the generator terminal voltage in normal operation of the generator at the generator port. After adding threephase voltage at both ends of the stator winding, in order to keep the flux chain of the winding unchanged, an instantaneous current is produced in the stator winding. The decay of these currents is as follows, the DC and frequency doubling components of the stator attenuate to zero at the moment of short circuit, and the AC component of the stator gradually attenuates from the initial value and finally attenuates to the steady value. In addition, most of the rotor current will decay to zero, leaving only the rotor DC.
As time goes by, the free current that appears in the shortcircuit state will decay to zero due to the resistance of the stator and rotor. The time constant will determine the attenuation of the current components of each winding. The larger the time constant, the slower the attenuation, and the smaller the time constant, the faster the attenuation. in synchronous generator, there are windings on stator and rotor, and the time constant of winding is related to inductance, then the time constant is expressed as the ratio of winding inductance value to winding resistance. The transient differential equation of winding is explained by the negative number of the ratio of winding resistance to winding inductance, that is, the characteristic root of differential equation [5]. When magnetic coupling exists in adjacent windings, mutual inductance should be considered to calculate the time constant, otherwise it will affect the accuracy of the final result. Under the circumstance of ensuring the constant flux linkage, studying the free current generated during short circuit, it can be concluded that according to the rotor excitation time constant Td’, it can be concluded that it attenuates according to the time constant of excitation winding which the DC and frequency doubling components of stator windings correspond to the fundamental frequency components of rotor current [6]. The time constant 9 is determined by the ratio of inductance of stator winding to resistance considering the mutual inductance of excitation winding [7, 8]. When this magnetic flux passes through the longitudinal axis of the rotor, the equivalent reactance of the stator winding is ${X}_{d}$; when it passes through the horizontal axis, it is ${X}_{q}$. The calculation process of the stator and rotor current is as follows, the $d$, $q$ axis component of stator current is calculated first, and then the rotor current is calculated by the integration of formula.
The $d$axis component of the stator current is:
The $q$axis component of the stator current is:
Through the finishing transformation of the above formula, the current in the rotor winding can be obtained as show in Eq. (3):
a phase current of stator is:
where ${X}_{d}$ is the synchronous reactance of stator winding straight axis (longitudinal axis); ${X}_{q}$ is the synchronous reactance of stator winding cross axis (transverse axis); ${X}_{f}$ is the armature reaction reactance between the longitudinal axis windings; ${X}_{ad}$ is the reactance of the longitudinal axis of the generator rotor; ${X}_{aq}$ is the reactance of the horizontal axis of the generator rotor; ${X}_{d}$, ${X}_{q}$_{}represent the reactance of d and q damper windings; ${X}_{d}^{\text{'}}$ is the transient reactance of longitudinal axis; ${X}_{d}^{\text{'}\text{'}}$^{}is the secondary transient reactance on the longitudinal axis; ${X}_{q}^{\text{'}\text{'}}$ is the secondary transient reactance of the longitudinal axis; ${E}_{q}^{\text{'}}$ is the transient electromotive force of the transverse axis; ${E}_{q}^{\text{'}\text{'}}$ is the secondary transient electromotive force of the transverse axis; ${E}_{d}^{\text{'}\text{'}}$ is the transient electromotive force of longitudinal axis; ${E}_{q}\left[0\right]$ is the instantaneous noload electromotive force before short circuit; $V\left[0\right]$ is the terminal voltage before short circuit.
3. Threephase short circuit simulation of synchronous generator
3.1. Threephase shortcircuit model of synchronous generator
A threephase shortcircuit simulation model of synchronous generator is established in Matlab/Simulink simulation environment. The model is shown in Fig. 1.
Fig. 1the Simulation model of threephase shortcircuit of synchronous generator
3.2. Parameter setting of simulation model
The basic parameters of the synchronous generator are as follows: rated power is ${P}_{N}=$200 MW, ${U}_{N}=$ 13.8 kV, $f=$ 50 HZ, ${x}_{d}=$ 1.0, ${x}_{d}^{\text{'}}=$ 0.3, ${x}_{d}^{\text{'}\text{'}}=$ 0.21, ${x}_{q}=$ 0.6, ${x}_{q}^{\text{'}\text{'}}=$ 0.31, ${x}_{1}=$ 0.15, ${T}_{d}^{\text{'}}=$ 1.64 s, ${T}_{d}^{\text{'}\text{'}}=$ 0.34 s, ${T}_{q0}^{\text{'}\text{'}}=$ 1.4 s, ${R}_{S}=$ 0.005.
The simulation model is established and the simulation parameters are modified. Through the Powergui module, the fault time is set to 0.02025 s, the simulation step is 0.01 s, the simulation time is set to 1 s, and the ode23tb algorithm is selected to solve the problem.
3.3. Analysis of simulation results
The waveform of stator directaxis current ${i}_{d}$ of synchronous generator is shown in Fig. 2. It can be seen from Fig. 2 that the stator directaxis current increases sharply at the moment of threephase short circuit of synchronous generator and the maximum current amplitude is 9 p.u., After reaching the maximum current, the amplitude attenuates rapidly with time constant and finally attenuates to 0 p.u. The waveform of stator crossaxis current ${i}_{q}$ in Fig. 3 shows that when threephase short circuit occurs in threephase fault device, the stator crossaxis current increases rapidly to the current amplitude, the maximum current is 3.02 p.u., then rapidly attenuates up and down to zero value, and the stator crossaxis current attenuates to 0 p.u. due to resistance.
Fig. 2The waveform diagram of stator straight axis current
Fig. 3The waveform of stator crossaxis current
The waveform of the excitation current ${i}_{f}$ is shown in Fig. 4. At the moment of short circuit of synchronous generator, the excitation current does not decay for a period of time after the fault, and increases rapidly to the current standard value of 2.8, which is due to the attenuation of damping winding current. In order to maintain the flux conservation, the excitation winding increases the current, and the excitation current continues to decay after increasing to the amplitude, but it will not attenuate to 0.
Fig. 4The waveform diagram of excitation current if
The threephase current waveform of stator when synchronous generator is shortcircuited is shown in Fig. 5. The amplitude of the threephase waveform is equal and the phase difference is 120 degrees before the failure. In the event of a threephase short circuit fault, the a phase current waveform first moves to the coordinate axis $y$ the positive half axis, then gradually moves to the coordinate axis $y$ the negative half axis. The moving direction of b phase and c phase is opposite to that of a phase. The threephase amplitude is equal, the phase difference is still 120 degrees, and the stator threephase current finally attenuates to the steadystate value.
Fig. 5The simulation waveform of stator current ia, ib, ic during threephase short circuit of synchronous generator
It can be seen from the above analysis that when the synchronous generator is operating normally, there is no directaxis and quadratureaxis currents. At the moment of short circuit, the vertical axis and horizontal axis currents gradually increase, and decay after a period of time after short circuit. The direct axis current increases faster than the quadrature axis current, and the attenuation is also slower than the quadrature axis. At the same time, the excitation current not only does not attenuate but increases sharply after short circuit. The main result is that the DC current attenuates quickly. The excitation winding attenuates the DC current and increases the excitation current to offset the attenuation of the DC current, thus maintaining the flux conservation. At the moment of short circuit of synchronous generator, stator current contains three components: fundamental frequency, DC frequency and frequency doubling component. When a threephase shortcircuit fault occurs, the A phase current increases the waveform first and then gradually moves down because of attenuation, while the B, C phase current waveform is opposite to the A phase current. The amplitude of the threephase current is equal, the phase difference is invariant, The stator threephase voltage decays directly to zero from the value before the shortcircuit moment. In general, the threephase short circuit of synchronous generator has an effect on the voltage and current, which makes the voltage drop to zero directly. For the current, the process is more complicated, but the total current is increased.
4. Conclusions
The proposed simulation model is convenient to modify various reactance parameters and time constants of the synchronous generator, and can accurately simulate the electromagnetic transient process of the synchronous generator. The expressions of stator current and rotor current of synchronous generator are given, which lay a foundation for further research on the influence of various reactance on shortcircuit current.
References

Best R. J., Morrow D. J., Crossley P. A. Current transients in the small salientpole alternator during sudden shortcircuit and synchronisation events. IET Electric Power Applications, Vol. 4, Issue 9, 2010, p. 687700.

Zhao L., Feng Y. Simulation and analysis of short circuit in 1000 MW synchronous generator based on Simulink/SPS module. AsiaPacific Power and Energy Engineering Conference, 2010.

Spunei E., Piroi I. Comparative analysis between stationary and dynamic parameters of a synchronous generator, with the main variable of the air gap magnetic induction. International Conference on Applied and Theoretical Electricity, 2012.

Wang J. X., Zobaa A. F., Bie Z. H., et al. From mathematical analysis to experimental calculation: teaching threephase shortcircuits of a synchronous generator. International Journal of Electrical Engineering Education, Vol. 49, Issue 4, 2012, p. 444463.

Sohn J., Hong S., Sunwoo M. Alternator torque model based on equivalent circuit of synchronous generator for electric power management. IEEE Transactions on Vehicular Technology, Vol. 62, Issue 8, 2013, p. 35933602.

Piroi I., Elisabeta S., Piroi F. Experimental measurements of the rotor oscillations in an synchronous generator during the threephased sudden shortcircuit test. International Conference on Applied and Theoretical Electricity, 2014.

Zhu Y. N., Jin J. X. Simulation analysis of a synchronous generator system under fault conditions. International Conference on Applied Superconductivity and Electromagnetic Devices, 2013, p. 180184.

Hackbart M. Novel approach to calculate electrical currents in stator, field and damperwindings at threephase sudden shortcircuit for large synchronous generators. Elektrotechnik Und Informationstechnik, Vol. 133, Issue 2, 2016, p. 112120.
About this article
This paper was supported by the following research projects: by the Natural Science Fund of Fujian Province (2019 J01845).