# Investigation into the Physical Mechanism and Influencing Factors of Sympathetic Inrush for Transformers in Series

Authors: Fang Peng; Houlei Gao
Journal: Mathematical Problems in Engineering
 (2017)

---

## Abstract

Sympathetic inrush is a magnetizing inrush phenomenon generated in an in-service transformer, when a nearby transformer in series or in parallel is energized. This phenomenon lasts longer than the general inrush and can cause protection maloperation. The study focuses on sympathetic inrush in series transformers. Flux linkage expressions for the series transformers are derived using the Laplace transform. The residual flux, switching angle, and system impedance factors are included in the expressions. The physical mechanism of the sympathetic inrush for series transformers is explained based on the characteristics of the flux linkages. The interaction between the saturated transformers during sympathetic inrush is also investigated. Moreover, the influence of the switching angle, residual flux, system impedance, and load on the inrush currents is analyzed and validated.


---

## Body

## 1. Introduction

The magnetizing inrush phenomenon has been well known since the transformer was first used in energy conversion. This phenomenon may cause maloperation of the current differential protection of the transformer. Moreover, the large electromagnetic forces in the transformer windings can damage the insulation of the transformer [1], and the voltage decreases can lead to poor power quality [2, 3] during magnetizing inrush. Magnetizing inrush can be classified into energization inrush, sympathetic inrush, and recovery inrush [4]. Energization inrush occurs in the transformer being energized, as a result of the reapplication of the system voltage to the transformer. When one transformer is in service and the other (in series or in parallel) is energized, sympathetic inrush occurs in the in-service transformer. Energization inrush also occurs in the energized transformer and is calledinitial inrush hereinafter. Additionally, the recovery inrush is generated during the restoration of the voltage after clearance. The present study focuses on sympathetic inrush.Several maloperation cases of current differential protection due to the sympathetic inrush have been reported [5]. Generally, the second harmonic component of the sympathetic inrush is high [6]. Thus, the second harmonic restraint algorithm is effective in preventing the protection of the in-service transformer from maloperation. However, the initial inrush and sympathetic inrush decay more slowly than the energization inrush in the single transformer due to the sympathetic interaction between the transformers [7]. Current transformer saturation caused by the long-lasting inrush current can decrease the second harmonic content, potentially causing the second harmonic restraint algorithm to fail [8–10]. Then, maloperation of the current differential protection of the in-service transformer will occur. Moreover, the sympathetic inrush may cause maloperation of other types of protection, such as the differential protection of the nearby generator and overcurrent protection, which leads to significant concern [8, 11].The mechanism of sympathetic inrush for transformers in parallel has been investigated previously [12–14]. However, no study has been performed on transformers in series (calledseries transformers below). In this context, the present study aims to perform an analysis of sympathetic inrush for series transformers. Their flux linkage expressions are derived; then, the mechanism of sympathetic inrush and the interaction between the series transformers are explained using the characteristics of the flux linkages. Moreover, the impacts of several factors including the switching angle, residual flux, system impedance, and load are also investigated. A practical model is used in the case studies to demonstrate the theoretical analysis.The rest of the paper is organized as follows. In Section2, the equivalent model of series transformers is presented, and the flux linkage expressions of the series transformers are derived based on the Laplace transform. In Section 3, the physical mechanism of the sympathetic inrush for series transformers is explained. In Section 4, the characteristics of the flux linkages under various levels of magnetizing inductance and the effect of transformer saturation on the inrush currents are investigated. In Section 5, the influence of the switching angle, residual flux, system impedance, and load on the sympathetic inrush is analyzed and verified. Finally, conclusions are drawn in Section 6.


## 2. Equivalent Model of Series Transformers during Sympathetic Inrush

### 2.1. Series Transformers and Equivalent Model

Series transformers are generally used in a step-down distribution system connected to large and special power consumers, such as traction systems, steel mills, and electric boilers. The series transformers are shown in Figure1. us is the power system equivalent source, and ZS=RS+jωLS is the equivalent impedance, where ω is the angular frequency. The circuit breaker (CB) is closed and T2 is energized. Initial inrush in T2 and sympathetic inrush in T1 are generated.Figure 1
Transformers connected in series.An equivalent T circuit is used to represent the transformer [15]. Thus, the equivalent circuit of series transformers is shown in Figure 2, where the symbols represent the following:Figure 2
Equivalent circuit of series transformers.R
c
1 and Lm1: the core loss resistance and magnetizing inductance of the iron core of T1,R
c
2 and Lm2: the core loss resistance and magnetizing inductance of the iron core of T2,R
11 and R12: the primary and secondary winding resistance of T1,R
21 and R22: the primary and secondary winding resistance of T2,L
σ
11 and Lσ12: the primary and secondary winding leakage inductance of T1,L
σ
21 and Lσ22: the primary and secondary winding leakage inductance of T2.All parameters are normalized to the primary side of T1.


### 2.2. Flux Linkage Deviation

The equivalent circuit is simplified as in Figure3 to investigate the transformer flux linkages. First, the current in the iron core consists of the core loss component and magnetizing component. The core loss resistance can be ignored since the core loss component is much smaller than the magnetizing component [13]. Second, the leakage inductance of the windings is far smaller than the magnetizing inductance of the iron core when the transformer is not saturated; therefore, Lσ can be neglected. Third, the secondary branch of T2 is inactive due to the lack of load. In addition, let R0=RS+R11, R1=R12, R2=R21, L0=LS, L1=Lm1, and L2=Lm2. Assume the source voltage varies sinusoidally with time; that is, us=Umsin⁡(ωt+φ), where Um and φ represent the amplitude and phase angle, respectively. The CB is closed at t=0 s.Figure 3
Simplified circuit of series transformers.The nonlinear characteristics of the iron core can be modeled using a piecewise linear inductance with two different slopes [16, 17]. However, the magnetizing inductance is assumed to be linear to simplify the analysis [13]. The purpose of the study is to qualitatively understand the relationship between the flux linkages and inrush currents for series transformers. Equation (1) can be obtained by Kirchhoff’s voltage law and Kirchhoff’s current law:(1)Umsin⁡ωt+φ=R0i0t+L0di0tdt+dψ1tdtdψ1tdt=dψ2tdt+R1+R2i2i0t=i1t+i2t,where i0(t), i1(t), and i2(t) represent the branch currents in L0, L1, and L2, respectively. ψ1(t)=L1i1(t) and ψ2(t)=L2i2(t) represent the flux linkages of T1 and T2, respectively.The Laplace transform of (1) is used to calculate the flux linkages in the frequency domain. Let R3=R1+R2. Then the flux linkage expressions in the time domain are obtained using the inverse Laplace transform, as presented in(2)ψ1t=-J2CFL2eτ1t1-AB2/4-ACB2A+AKJ+eτ2t1+AB2/4-ACB2A+AKJ+R32+ωL22L2F-Gsin⁡ωt+φ+θ1+BDωsin⁡ωt+φ-θ1ψ2t=H2F·eτ1t1-AB2/4-ACB2A-AIH+eτ2t1+1B2/4-ACB2A-AIH+1FωBDω2+G2sin⁡ωt+φ+θ2,where(3)A=L0L1+L1L2+L0L2B=R0L1+R0L2+R3L0+R3L1C=R0R3D=L1L2UmE=-R0L2ψ10--ψ20-F=A2ω4-2ACω2+B2ω2+C2G=ADω2-CDH=ψ20-F+ωGI=EF+CGsin⁡φ-BDωcos⁡φJ=-CGR3-BDω2L2sin⁡φ+GωL2+BDωR3cos⁡φ+EFR3-CFL2ψ20-K=ACFR3ψ20-+EFCL2-BR3-GL2+BDR3sin⁡φ+AGR3+BDBR3-CL2Cωcos⁡φτ1=-B+B2-4AC2Aτ2=-B-B2-4AC2Aθ1=arctan⁡ωL2R3θ2=arctan⁡-GBDω.Here,τ1 and τ2 represent the decay coefficients of the DC components. ψ1(0-) and ψ2(0-) are the initial flux linkages of T1 and T2 at t=0 s, respectively. ψ2(0-) represents the residual flux linkage of the energized transformer T2. Therefore, the effect of residual flux can be included in (2).Definea1, a2, a3, b1, b2, c1, and c2 as in (4)a1=-GR32+ωL22L2Fa2=BDR32+ωL22L2Fa3=ωBDω2+G2Fb1=-J2CFL21-AB2/4-ACB2A+AKJb2=H2F1-AB2/4-ACB2A-AIHc1=-J2CFL21+AB2/4-ACB2A+AKJc2=H2F1+AB2/4-ACB2A-AIH.a
1, a2, and a3 represent the amplitude of the AC components, while b1, b2, c1, and c2 represent the amplitude of the DC components. Then, the transformer flux linkages are simplified as follows: (5)ψ1t=a1sin⁡ωt+φ+θ1+a2sin⁡ωt+φ-θ1+b1eτ1t+c1eτ2tψ2t=a3sin⁡ωt+φ+θ2+b2eτ1t+c2eτ2t.Bothψ1(t) and ψ2(t) consist of the AC component and decaying DC component. The two AC components of ψ1(t) can be simplified as an AC component, which is explained in detail in Section 3. The DC components consist of the τ1 component and τ2 component for both transformers.


## 2.1. Series Transformers and Equivalent Model

Series transformers are generally used in a step-down distribution system connected to large and special power consumers, such as traction systems, steel mills, and electric boilers. The series transformers are shown in Figure1. us is the power system equivalent source, and ZS=RS+jωLS is the equivalent impedance, where ω is the angular frequency. The circuit breaker (CB) is closed and T2 is energized. Initial inrush in T2 and sympathetic inrush in T1 are generated.Figure 1
Transformers connected in series.An equivalent T circuit is used to represent the transformer [15]. Thus, the equivalent circuit of series transformers is shown in Figure 2, where the symbols represent the following:Figure 2
Equivalent circuit of series transformers.R
c
1 and Lm1: the core loss resistance and magnetizing inductance of the iron core of T1,R
c
2 and Lm2: the core loss resistance and magnetizing inductance of the iron core of T2,R
11 and R12: the primary and secondary winding resistance of T1,R
21 and R22: the primary and secondary winding resistance of T2,L
σ
11 and Lσ12: the primary and secondary winding leakage inductance of T1,L
σ
21 and Lσ22: the primary and secondary winding leakage inductance of T2.All parameters are normalized to the primary side of T1.


## 2.2. Flux Linkage Deviation

The equivalent circuit is simplified as in Figure3 to investigate the transformer flux linkages. First, the current in the iron core consists of the core loss component and magnetizing component. The core loss resistance can be ignored since the core loss component is much smaller than the magnetizing component [13]. Second, the leakage inductance of the windings is far smaller than the magnetizing inductance of the iron core when the transformer is not saturated; therefore, Lσ can be neglected. Third, the secondary branch of T2 is inactive due to the lack of load. In addition, let R0=RS+R11, R1=R12, R2=R21, L0=LS, L1=Lm1, and L2=Lm2. Assume the source voltage varies sinusoidally with time; that is, us=Umsin⁡(ωt+φ), where Um and φ represent the amplitude and phase angle, respectively. The CB is closed at t=0 s.Figure 3
Simplified circuit of series transformers.The nonlinear characteristics of the iron core can be modeled using a piecewise linear inductance with two different slopes [16, 17]. However, the magnetizing inductance is assumed to be linear to simplify the analysis [13]. The purpose of the study is to qualitatively understand the relationship between the flux linkages and inrush currents for series transformers. Equation (1) can be obtained by Kirchhoff’s voltage law and Kirchhoff’s current law:(1)Umsin⁡ωt+φ=R0i0t+L0di0tdt+dψ1tdtdψ1tdt=dψ2tdt+R1+R2i2i0t=i1t+i2t,where i0(t), i1(t), and i2(t) represent the branch currents in L0, L1, and L2, respectively. ψ1(t)=L1i1(t) and ψ2(t)=L2i2(t) represent the flux linkages of T1 and T2, respectively.The Laplace transform of (1) is used to calculate the flux linkages in the frequency domain. Let R3=R1+R2. Then the flux linkage expressions in the time domain are obtained using the inverse Laplace transform, as presented in(2)ψ1t=-J2CFL2eτ1t1-AB2/4-ACB2A+AKJ+eτ2t1+AB2/4-ACB2A+AKJ+R32+ωL22L2F-Gsin⁡ωt+φ+θ1+BDωsin⁡ωt+φ-θ1ψ2t=H2F·eτ1t1-AB2/4-ACB2A-AIH+eτ2t1+1B2/4-ACB2A-AIH+1FωBDω2+G2sin⁡ωt+φ+θ2,where(3)A=L0L1+L1L2+L0L2B=R0L1+R0L2+R3L0+R3L1C=R0R3D=L1L2UmE=-R0L2ψ10--ψ20-F=A2ω4-2ACω2+B2ω2+C2G=ADω2-CDH=ψ20-F+ωGI=EF+CGsin⁡φ-BDωcos⁡φJ=-CGR3-BDω2L2sin⁡φ+GωL2+BDωR3cos⁡φ+EFR3-CFL2ψ20-K=ACFR3ψ20-+EFCL2-BR3-GL2+BDR3sin⁡φ+AGR3+BDBR3-CL2Cωcos⁡φτ1=-B+B2-4AC2Aτ2=-B-B2-4AC2Aθ1=arctan⁡ωL2R3θ2=arctan⁡-GBDω.Here,τ1 and τ2 represent the decay coefficients of the DC components. ψ1(0-) and ψ2(0-) are the initial flux linkages of T1 and T2 at t=0 s, respectively. ψ2(0-) represents the residual flux linkage of the energized transformer T2. Therefore, the effect of residual flux can be included in (2).Definea1, a2, a3, b1, b2, c1, and c2 as in (4)a1=-GR32+ωL22L2Fa2=BDR32+ωL22L2Fa3=ωBDω2+G2Fb1=-J2CFL21-AB2/4-ACB2A+AKJb2=H2F1-AB2/4-ACB2A-AIHc1=-J2CFL21+AB2/4-ACB2A+AKJc2=H2F1+AB2/4-ACB2A-AIH.a
1, a2, and a3 represent the amplitude of the AC components, while b1, b2, c1, and c2 represent the amplitude of the DC components. Then, the transformer flux linkages are simplified as follows: (5)ψ1t=a1sin⁡ωt+φ+θ1+a2sin⁡ωt+φ-θ1+b1eτ1t+c1eτ2tψ2t=a3sin⁡ωt+φ+θ2+b2eτ1t+c2eτ2t.Bothψ1(t) and ψ2(t) consist of the AC component and decaying DC component. The two AC components of ψ1(t) can be simplified as an AC component, which is explained in detail in Section 3. The DC components consist of the τ1 component and τ2 component for both transformers.


## 3. Transformer Flux Linkage Analysis

### 3.1. AC Component of the Flux Linkages

Before the CB is closed, the transformer T1 is not saturated. Since the impedance of the magnetizing branch is much larger than the system impedance in this situation, the former divides most of the source voltage. Thus, the flux linkage of T1 is approximately equal to the integral of the source voltage, as expressed in(6)ψ1t=∫ustdt=Umω·sin⁡ωt+φ-π2.Only the AC component is included in the flux linkage of T1 before the CB is closed. ψ1(0-)=-Um/ω can be obtained from (6).After the CB is closed, the AC components of T1 and T2 can be obtained from (2) and (5). The magnetizing inductance changes from several henries to several hundred henries for different types of transformers with different rated voltages, capacity, and so on. In comparison, the maximum winding resistance is several ohms, which means that ωL2≫R3 and θ1≈π/2. Thus, ψ1AC(t)=a1sin⁡ωt+φ+θ1+a2sin⁡(ωt+φ-θ1) can be written as(7)ψ1ACt=DR32+ωL22Aω2-C+BDωL2Aω2-C2+Bω2sin⁡ωt+φ-π2=a4sin⁡ωt+φ-π2.As forψ2AC=a3sin⁡(ωt+φ+θ2), θ2≈-π/2 can be obtained because Aω2-C≫B. Similarly, it can be verified that a4 and a2 are approximately equal to Um/ω. Therefore, both ψ1AC and ψ2AC are approximately equal to the integral of the source voltage, that is, Um/ω·sin⁡(ωt+φ-π/2). The variation of R and L will not change this feature.Series transformers typically differ in their capacity and voltage level. After being normalized to the primary side of T1, a group of equivalent parameters is taken as an example:(8)RS=10ΩLS=0.3ΗR11=2ΩLm1=67ΗR12=2ΩLm2=130ΗR21=3Ωω=100πrad/s.The AC components of flux linkages for T1 and T2 can be obtained from (2) and (8), as shown in Figure 4. The base value is taken as the RMS value. The waveform is consistent with the analysis above.Figure 4
AC components of the flux linkages.


### 3.2. Decaying DC Component

As shown in Figure5, the decaying DC components of the flux linkages are also obtained from (2) and (8), where ψ1DC=b1eτ1t+c1eτ2t and ψ2DC=b2eτ1t+c2eτ2t.Figure 5
DC components of the flux linkages.The DC components of the two flux linkages have opposite polarities. The amplitude ofψ2DC continues decaying from the initial maximum (|ψ2DC_max|=1.405), while ψ1DC increases first and then decreases (|ψ1DC_max|=0.364). The former has a peak value at the switching instant (t=0 s), and the latter has a peak value after several seconds (t=9.4 s). As shown in Figure 6, each flux linkage’s DC component consists of two DC components with different decay coefficients (τ1 component and τ2 component, resp.).Figure 6
τ
1 and τ2 components of the flux linkages.Although the decay coefficients ofψ1DC and ψ2DC are identical, their polarities and amplitudes are different. As a result, ψ1DC and ψ2DC change in different ways. Two positive DC components (b2eτ1t and c2eτ2t) exist in ψ2DC. The superimposition of the two components is still positive and continually decays over time. In comparison, a positive component c1eτ2t and a negative component b1eτ1t exist in ψ1DC. Since the negative component initially decays more slowly than the positive component, the amplitude of the superimposition component ψ1DC increases in the negative polarity. After several seconds, the negative component decays faster than the positive component; therefore, the amplitude of the superimposition component ψ1DC decreases in the negative polarity. Furthermore, the initial value of ψ1DC is approximately zero because the two components in opposite polarities have approximately equal amplitude at t=0 s.


### 3.3. Sympathetic Inrush Case

A transformer typically operates near the knee point of the saturation curve. Only a single AC component is included in the transformer flux linkage when the transformer operates in the steady state. However, the superimposition of the DC component makes the transformer operate far above the knee point. Then, the transformer is driven into saturation and the magnetizing inrush is generated. If a positive DC component is added to the AC component in the positive half cycle, then positive inrush is generated. If a negative DC component is added to the AC component in the negative half cycle, negative inrush will be generated. During sympathetic inrush, the DC components of the two transformers have opposite polarities, while the AC components have the same polarity. Therefore, the initial inrush and sympathetic inrush are alternately generated in the positive and negative half cycle of the AC component. Further explanation of the process can be found in Section4.1. The physical mechanism of the sympathetic inrush for series transformers is similar to that for parallel transformers [13].A sympathetic inrush simulation model of series transformers is constructed as Figure7 in Power Systems Computer Aided Design (PSCAD) according to a real maloperation case in China. The system parameters are US=235 kV (RMS value of us), ω=100π rad/s, and ZS=30+j283 Ω. The parameters of T1 and T2 are listed in Table 1, where the magnetizing reactance represents the reactance of the iron core that is fully saturated.Table 1
Parameters of T1 and T2.
Rated capacity (MVA)
Voltage ratio (kV/kV)
Winding connection
Knee voltage (p.u.)
Magnetizing reactance (p.u.)
T1
31.5
220/6
Y/Y
1.12
0.14
T2
26
6.3/13.8
Δ/Y
1.15
0.3Figure 7
Series transformers simulation model.In the simulation, the CB is closed att=0.5 s. T2 is energized when T1 is in service. Taking Phase A as an example, the switching angle of Phase A is 0 (φA = 0°). The magnetizing inrush currents corresponding to the inrush phenomenon can be obtained from the differential currents of T1 and T2. The differential currents IT1dA and IT2dA are shown in Figure 8, which are normalized to the 6 kV side. Positive initial inrush is generated in T2 and negative sympathetic inrush is generated in T1. Only the Phase A currents are shown because the inrush currents of Phase A are the largest. The variation of the flux linkages is shown in Figure 9. ψ1 consists of the AC component and negative DC component, while ψ2 consists of the AC component and positive DC component.Figure 8
Initial inrush and sympathetic inrush of series transformers.Figure 9
Flux linkages of series transformers.


### 3.4. Interaction between the Series Transformers

If there is a single in-service transformer in the power system, only an AC component will be included in the transformer flux linkage.ψ1DC is generated due to the other energized transformer. Thus, it is the interaction between the series transformers that generates the sympathetic inrush.The flux linkage of a single energized transformerψ2tsg can also be derived in a similar way to (2): (9)ψ2tsg=L2Umsin⁡ωt+φ-θ3R42+ωL32+ψ20--L2Umsin⁡φ-θ3R42+ωL32e-R4/L3t,where R4=R0+R2, L3=L0+L2, and θ3=arctan⁡(ωL3/R4).An AC componentψ2ACsg and a DC component ψ2DCsg are included in the flux linkage. θ3≈-π/2 can be obtained, and ψ2ACsg is from the sinusoidal source voltage. It can be seen from (2) and (9) that the in-service transformer changes the decay coefficient of the DC component and increases a DC component to the flux linkage of the energized transformer. ψ2DCsg can be obtained from (8) and (9), as shown in Figure 10.Figure 10
DC component of the energized transformer.ψ
2
D
C for series transformers is also shown in Figure 10 as a comparison. ψ2DC decays more slowly than ψ2DCsg. As a result, the initial inrush decays more slowly due to the in-service transformer.


## 3.1. AC Component of the Flux Linkages

Before the CB is closed, the transformer T1 is not saturated. Since the impedance of the magnetizing branch is much larger than the system impedance in this situation, the former divides most of the source voltage. Thus, the flux linkage of T1 is approximately equal to the integral of the source voltage, as expressed in(6)ψ1t=∫ustdt=Umω·sin⁡ωt+φ-π2.Only the AC component is included in the flux linkage of T1 before the CB is closed. ψ1(0-)=-Um/ω can be obtained from (6).After the CB is closed, the AC components of T1 and T2 can be obtained from (2) and (5). The magnetizing inductance changes from several henries to several hundred henries for different types of transformers with different rated voltages, capacity, and so on. In comparison, the maximum winding resistance is several ohms, which means that ωL2≫R3 and θ1≈π/2. Thus, ψ1AC(t)=a1sin⁡ωt+φ+θ1+a2sin⁡(ωt+φ-θ1) can be written as(7)ψ1ACt=DR32+ωL22Aω2-C+BDωL2Aω2-C2+Bω2sin⁡ωt+φ-π2=a4sin⁡ωt+φ-π2.As forψ2AC=a3sin⁡(ωt+φ+θ2), θ2≈-π/2 can be obtained because Aω2-C≫B. Similarly, it can be verified that a4 and a2 are approximately equal to Um/ω. Therefore, both ψ1AC and ψ2AC are approximately equal to the integral of the source voltage, that is, Um/ω·sin⁡(ωt+φ-π/2). The variation of R and L will not change this feature.Series transformers typically differ in their capacity and voltage level. After being normalized to the primary side of T1, a group of equivalent parameters is taken as an example:(8)RS=10ΩLS=0.3ΗR11=2ΩLm1=67ΗR12=2ΩLm2=130ΗR21=3Ωω=100πrad/s.The AC components of flux linkages for T1 and T2 can be obtained from (2) and (8), as shown in Figure 4. The base value is taken as the RMS value. The waveform is consistent with the analysis above.Figure 4
AC components of the flux linkages.


## 3.2. Decaying DC Component

As shown in Figure5, the decaying DC components of the flux linkages are also obtained from (2) and (8), where ψ1DC=b1eτ1t+c1eτ2t and ψ2DC=b2eτ1t+c2eτ2t.Figure 5
DC components of the flux linkages.The DC components of the two flux linkages have opposite polarities. The amplitude ofψ2DC continues decaying from the initial maximum (|ψ2DC_max|=1.405), while ψ1DC increases first and then decreases (|ψ1DC_max|=0.364). The former has a peak value at the switching instant (t=0 s), and the latter has a peak value after several seconds (t=9.4 s). As shown in Figure 6, each flux linkage’s DC component consists of two DC components with different decay coefficients (τ1 component and τ2 component, resp.).Figure 6
τ
1 and τ2 components of the flux linkages.Although the decay coefficients ofψ1DC and ψ2DC are identical, their polarities and amplitudes are different. As a result, ψ1DC and ψ2DC change in different ways. Two positive DC components (b2eτ1t and c2eτ2t) exist in ψ2DC. The superimposition of the two components is still positive and continually decays over time. In comparison, a positive component c1eτ2t and a negative component b1eτ1t exist in ψ1DC. Since the negative component initially decays more slowly than the positive component, the amplitude of the superimposition component ψ1DC increases in the negative polarity. After several seconds, the negative component decays faster than the positive component; therefore, the amplitude of the superimposition component ψ1DC decreases in the negative polarity. Furthermore, the initial value of ψ1DC is approximately zero because the two components in opposite polarities have approximately equal amplitude at t=0 s.


## 3.3. Sympathetic Inrush Case

A transformer typically operates near the knee point of the saturation curve. Only a single AC component is included in the transformer flux linkage when the transformer operates in the steady state. However, the superimposition of the DC component makes the transformer operate far above the knee point. Then, the transformer is driven into saturation and the magnetizing inrush is generated. If a positive DC component is added to the AC component in the positive half cycle, then positive inrush is generated. If a negative DC component is added to the AC component in the negative half cycle, negative inrush will be generated. During sympathetic inrush, the DC components of the two transformers have opposite polarities, while the AC components have the same polarity. Therefore, the initial inrush and sympathetic inrush are alternately generated in the positive and negative half cycle of the AC component. Further explanation of the process can be found in Section4.1. The physical mechanism of the sympathetic inrush for series transformers is similar to that for parallel transformers [13].A sympathetic inrush simulation model of series transformers is constructed as Figure7 in Power Systems Computer Aided Design (PSCAD) according to a real maloperation case in China. The system parameters are US=235 kV (RMS value of us), ω=100π rad/s, and ZS=30+j283 Ω. The parameters of T1 and T2 are listed in Table 1, where the magnetizing reactance represents the reactance of the iron core that is fully saturated.Table 1
Parameters of T1 and T2.
Rated capacity (MVA)
Voltage ratio (kV/kV)
Winding connection
Knee voltage (p.u.)
Magnetizing reactance (p.u.)
T1
31.5
220/6
Y/Y
1.12
0.14
T2
26
6.3/13.8
Δ/Y
1.15
0.3Figure 7
Series transformers simulation model.In the simulation, the CB is closed att=0.5 s. T2 is energized when T1 is in service. Taking Phase A as an example, the switching angle of Phase A is 0 (φA = 0°). The magnetizing inrush currents corresponding to the inrush phenomenon can be obtained from the differential currents of T1 and T2. The differential currents IT1dA and IT2dA are shown in Figure 8, which are normalized to the 6 kV side. Positive initial inrush is generated in T2 and negative sympathetic inrush is generated in T1. Only the Phase A currents are shown because the inrush currents of Phase A are the largest. The variation of the flux linkages is shown in Figure 9. ψ1 consists of the AC component and negative DC component, while ψ2 consists of the AC component and positive DC component.Figure 8
Initial inrush and sympathetic inrush of series transformers.Figure 9
Flux linkages of series transformers.


## 3.4. Interaction between the Series Transformers

If there is a single in-service transformer in the power system, only an AC component will be included in the transformer flux linkage.ψ1DC is generated due to the other energized transformer. Thus, it is the interaction between the series transformers that generates the sympathetic inrush.The flux linkage of a single energized transformerψ2tsg can also be derived in a similar way to (2): (9)ψ2tsg=L2Umsin⁡ωt+φ-θ3R42+ωL32+ψ20--L2Umsin⁡φ-θ3R42+ωL32e-R4/L3t,where R4=R0+R2, L3=L0+L2, and θ3=arctan⁡(ωL3/R4).An AC componentψ2ACsg and a DC component ψ2DCsg are included in the flux linkage. θ3≈-π/2 can be obtained, and ψ2ACsg is from the sinusoidal source voltage. It can be seen from (2) and (9) that the in-service transformer changes the decay coefficient of the DC component and increases a DC component to the flux linkage of the energized transformer. ψ2DCsg can be obtained from (8) and (9), as shown in Figure 10.Figure 10
DC component of the energized transformer.ψ
2
D
C for series transformers is also shown in Figure 10 as a comparison. ψ2DC decays more slowly than ψ2DCsg. As a result, the initial inrush decays more slowly due to the in-service transformer.


## 4. Effect of Transformer Saturation on Inrush Currents

Although the linear magnetizing inductance assumption is made to simplify the transformer flux linkage analysis above, the magnetizing inductance of the saturated transformer is actually nonlinear. The magnetic characteristics of the transformer are typically represented using the piecewise lines. The magnetizing inductance will significantly decrease if the transformer is driven into saturation. In this section, the flux linkages of saturated transformers are analyzed to investigate the effect of transformer saturation.


### 4.1. Variation of Magnetizing Inductance during Sympathetic Inrush

The piecewise lines representing the magnetizing characteristics are illustrated in Figure11(a). The magnetizing inductance below the knee point (nonsaturation zone) is Lnon, and that above the knee point (saturation zone) is Lsat. Since the equivalent inductance around the knee point varies dramatically, Lnon is much larger than Lsat. As shown in Figure 11(a), the addition of the DC components to the AC components during one cycle of the flux linkages can be used to explain the variation of the magnetizing inductance further. As analyzed in Section 3, the DC components of the two transformers have opposite polarities, and the AC components are approximately equal.Figure 11
(a) The flux linkages and piecewise linear magnetization curve. (b) The equivalent magnetizing inductance during sympathetic inrush.
(a)
(b)(1) Beforet1, positive ψ2AC and ψ2DC are added together, but ψ2AC is small, and thus ψ2=ψ2AC+ψ2DC is below the knee point. Positive ψ1AC and negative ψ1DC also cause ψ1=ψ1AC+ψ1DC to be below the knee point. Therefore, T1 and T2 are not saturated and the magnetizing inductance of each transformer is Lnon. (2) As the instantaneous value of ψAC increases to the peak value during t1-t2, ψ2 is above the knee point. T2 is driven into saturation while the positive ψ1AC and negative ψ1DC keep ψ1 below the knee point. Therefore, the magnetizing inductance of T1 is Lnon, while that of T2 is Lsat. (3) As the instantaneous value of ψAC decreases from t2, T2 enters the nonsaturation zone. As ψAC decreases further, T1 is driven into saturation in the negative half cycle in a similar way.Therefore, T1 and T2 (T1/T2) go through four stages sequentially in each cycle during sympathetic inrush: nonsaturation/nonsaturation, nonsaturation/saturation, nonsaturation/nonsaturation, and saturation/nonsaturation. The corresponding four stages are for the periods0-t1, t1-t2, t2–t4, and t4-t5. In the period t5-t6, T1 and T2 keep unsaturated to the next cycle. The magnetizing inductance values will change with different stages, as shown in Figure 11(b). When the sympathetic inrush is generated, the equivalent magnetizing inductance of the in-service transformer is alternately presented at Lnon and Lsat. The same effect occurs for the magnetizing inductance of the energized transformer when the initial inrush is generated. In addition, Lnon and Lsat of the two transformers are actually different. Lnon1 and Lsat1 are the magnetizing inductance of T1 while Lnon2 and Lsat2 are the magnetizing inductance of T2.For the AC components, the decrease of the magnetizing inductance will not change the feature ofψ1AC and ψ2AC because the conditions ωL2≫R3 and Aω2-C≫B still hold true, even if the magnetizing inductance decreases. Therefore, ψ1AC and ψ2AC are both approximately equal to Um/ω·sin⁡(ωt+φ-π/2). The saturation of the transformers has a slight impact on the AC components.For the DC components, the magnetizing inductance is assumed to be linear in the DC component analysis in Section3. This assumption can be explained further as follows. When the CB is closed, the DC component of the flux linkage is generated to prevent the flux linkage from changing discretely. At the switching instant, the magnetizing inductance is Lnon, in the period 0-t1. Therefore, the variation feature of the DC component is determined predominantly by Lnon at the switching instant. The main characteristic of the flux linkage will not change even if the transformer is driven into saturated saturation. The reason is that the circuit is no longer forced to change and there is no new DC component generated to prevent a discrete change in the flux linkage. The effect of transformer saturation is to change the amplitude and decay rate to some degree.Further research is implemented to analyze the influence of the decrease of magnetizing inductance on the DC components.


### 4.2. DC Components of Flux Linkages When T2 Is Saturated

When T2 is saturated at thet1-t2 stage, Lm2 is reduced to Lsat2, while Lm1=Lnon1. The DC components of the flux linkages are obtained by taking a smaller Lm2 into (2), where Lsat2=60 H or 20 H. Although the stage of Lm2=Lsat2 is only a short period in a cycle, the effect is analyzed qualitatively using Lsat2 in the entire decay process. The DC components of the flux linkages are shown in Figure 12. ψ2DC decays more quickly in the main decay process and ψ1DC has an earlier and larger peak value with the decrease in Lm2. The reason is that the faster decaying component of ψ1DC (τ2 component) decays much faster and the amplitude of each component of ψ1DC (τ1 and τ2 components) increases. The τ1 and τ2 components under different Lm2 values are shown in Figure 13, and the decay coefficients are listed in Table 2. |b1|, |c1|, and |c2| increase, while |b2| decreases with decreasing Lm2. |τ1| and |τ2| increase with decreasing Lm2.Table 2
Decay coefficients varying withLm2.
L
m
2 (H)
130
60
20
τ
1
−0.0241
−0.0349
−0.0455
τ
2
−0.2839
−0.4236
−0.9650Figure 12
DC components of the flux linkages under differentLm2.Figure 13
τ
1 and τ2 components of the flux linkages under different Lm2.As a result, the initial inrush and sympathetic inrush are influenced. The initial inrush will decay more quickly, and the sympathetic inrush will be strengthened.


### 4.3. DC Components of Flux Linkages When T1 Is Saturated

When T1 is saturated at thet4-t5 stage, Lm1 is reduced to Lsat1, while Lm2=Lnon2. The DC components of the flux linkages are then obtained by taking a smaller Lm1 into (2), where Lsat1 = 30 H or 10 H. The DC components are shown in Figure 14. ψ2DC decays more slowly, and ψ1DC decreases with decreasing Lm1. When Lm1 decreases to 10 H, ψ1DC is too small to lead to the sympathetic inrush. τ1 and τ2 components of the flux linkages are shown in Figure 15. |b1|, |c1|, and |c2| decrease while |b2| increases with decreasing Lm1. As listed in Table 3, |τ1| and |τ2| increase with decreasing Lm1. As a result, the decay rate of each component increases. ψ2DC decays more slowly with decreasing Lm1 because the more slowly decaying component of ψ2DC (τ1 component) has a larger amplitude.Table 3
Decay coefficients varying withLm1.
L
m
1 (H)
67
30
10
τ
1
−0.0241
−0.0308
−0.0356
τ
2
−0.2839
−0.4939
−1.2546Figure 14
DC components of the flux linkages under differentLm1.Figure 15
τ
1 and τ2 components of the flux linkages under different Lm1.In regard to the inrush currents of the transformers, the sympathetic inrush is weakened and the initial inrush decays more slowly due to the saturation of T1.In conclusion, the series transformers interact with each other due to the transformer saturation. The initial inrush decays more quickly due to the saturation of the energized transformer, and the sympathetic inrush is weakened due to saturation of the in-service transformer. The saturation of the transformer resists further saturation, which can be thought of as a self-regulating mechanism. On the other hand, the initial inrush decays more slowly due to saturation of the in-service transformer, and the sympathetic inrush is strengthened due to saturation of the energized transformer. The saturation of one transformer strengthens the other transformer’s saturation, which can be explained in terms of the energy balance of the two transformers.


## 4.1. Variation of Magnetizing Inductance during Sympathetic Inrush

The piecewise lines representing the magnetizing characteristics are illustrated in Figure11(a). The magnetizing inductance below the knee point (nonsaturation zone) is Lnon, and that above the knee point (saturation zone) is Lsat. Since the equivalent inductance around the knee point varies dramatically, Lnon is much larger than Lsat. As shown in Figure 11(a), the addition of the DC components to the AC components during one cycle of the flux linkages can be used to explain the variation of the magnetizing inductance further. As analyzed in Section 3, the DC components of the two transformers have opposite polarities, and the AC components are approximately equal.Figure 11
(a) The flux linkages and piecewise linear magnetization curve. (b) The equivalent magnetizing inductance during sympathetic inrush.
(a)
(b)(1) Beforet1, positive ψ2AC and ψ2DC are added together, but ψ2AC is small, and thus ψ2=ψ2AC+ψ2DC is below the knee point. Positive ψ1AC and negative ψ1DC also cause ψ1=ψ1AC+ψ1DC to be below the knee point. Therefore, T1 and T2 are not saturated and the magnetizing inductance of each transformer is Lnon. (2) As the instantaneous value of ψAC increases to the peak value during t1-t2, ψ2 is above the knee point. T2 is driven into saturation while the positive ψ1AC and negative ψ1DC keep ψ1 below the knee point. Therefore, the magnetizing inductance of T1 is Lnon, while that of T2 is Lsat. (3) As the instantaneous value of ψAC decreases from t2, T2 enters the nonsaturation zone. As ψAC decreases further, T1 is driven into saturation in the negative half cycle in a similar way.Therefore, T1 and T2 (T1/T2) go through four stages sequentially in each cycle during sympathetic inrush: nonsaturation/nonsaturation, nonsaturation/saturation, nonsaturation/nonsaturation, and saturation/nonsaturation. The corresponding four stages are for the periods0-t1, t1-t2, t2–t4, and t4-t5. In the period t5-t6, T1 and T2 keep unsaturated to the next cycle. The magnetizing inductance values will change with different stages, as shown in Figure 11(b). When the sympathetic inrush is generated, the equivalent magnetizing inductance of the in-service transformer is alternately presented at Lnon and Lsat. The same effect occurs for the magnetizing inductance of the energized transformer when the initial inrush is generated. In addition, Lnon and Lsat of the two transformers are actually different. Lnon1 and Lsat1 are the magnetizing inductance of T1 while Lnon2 and Lsat2 are the magnetizing inductance of T2.For the AC components, the decrease of the magnetizing inductance will not change the feature ofψ1AC and ψ2AC because the conditions ωL2≫R3 and Aω2-C≫B still hold true, even if the magnetizing inductance decreases. Therefore, ψ1AC and ψ2AC are both approximately equal to Um/ω·sin⁡(ωt+φ-π/2). The saturation of the transformers has a slight impact on the AC components.For the DC components, the magnetizing inductance is assumed to be linear in the DC component analysis in Section3. This assumption can be explained further as follows. When the CB is closed, the DC component of the flux linkage is generated to prevent the flux linkage from changing discretely. At the switching instant, the magnetizing inductance is Lnon, in the period 0-t1. Therefore, the variation feature of the DC component is determined predominantly by Lnon at the switching instant. The main characteristic of the flux linkage will not change even if the transformer is driven into saturated saturation. The reason is that the circuit is no longer forced to change and there is no new DC component generated to prevent a discrete change in the flux linkage. The effect of transformer saturation is to change the amplitude and decay rate to some degree.Further research is implemented to analyze the influence of the decrease of magnetizing inductance on the DC components.


## 4.2. DC Components of Flux Linkages When T2 Is Saturated

When T2 is saturated at thet1-t2 stage, Lm2 is reduced to Lsat2, while Lm1=Lnon1. The DC components of the flux linkages are obtained by taking a smaller Lm2 into (2), where Lsat2=60 H or 20 H. Although the stage of Lm2=Lsat2 is only a short period in a cycle, the effect is analyzed qualitatively using Lsat2 in the entire decay process. The DC components of the flux linkages are shown in Figure 12. ψ2DC decays more quickly in the main decay process and ψ1DC has an earlier and larger peak value with the decrease in Lm2. The reason is that the faster decaying component of ψ1DC (τ2 component) decays much faster and the amplitude of each component of ψ1DC (τ1 and τ2 components) increases. The τ1 and τ2 components under different Lm2 values are shown in Figure 13, and the decay coefficients are listed in Table 2. |b1|, |c1|, and |c2| increase, while |b2| decreases with decreasing Lm2. |τ1| and |τ2| increase with decreasing Lm2.Table 2
Decay coefficients varying withLm2.
L
m
2 (H)
130
60
20
τ
1
−0.0241
−0.0349
−0.0455
τ
2
−0.2839
−0.4236
−0.9650Figure 12
DC components of the flux linkages under differentLm2.Figure 13
τ
1 and τ2 components of the flux linkages under different Lm2.As a result, the initial inrush and sympathetic inrush are influenced. The initial inrush will decay more quickly, and the sympathetic inrush will be strengthened.


## 4.3. DC Components of Flux Linkages When T1 Is Saturated

When T1 is saturated at thet4-t5 stage, Lm1 is reduced to Lsat1, while Lm2=Lnon2. The DC components of the flux linkages are then obtained by taking a smaller Lm1 into (2), where Lsat1 = 30 H or 10 H. The DC components are shown in Figure 14. ψ2DC decays more slowly, and ψ1DC decreases with decreasing Lm1. When Lm1 decreases to 10 H, ψ1DC is too small to lead to the sympathetic inrush. τ1 and τ2 components of the flux linkages are shown in Figure 15. |b1|, |c1|, and |c2| decrease while |b2| increases with decreasing Lm1. As listed in Table 3, |τ1| and |τ2| increase with decreasing Lm1. As a result, the decay rate of each component increases. ψ2DC decays more slowly with decreasing Lm1 because the more slowly decaying component of ψ2DC (τ1 component) has a larger amplitude.Table 3
Decay coefficients varying withLm1.
L
m
1 (H)
67
30
10
τ
1
−0.0241
−0.0308
−0.0356
τ
2
−0.2839
−0.4939
−1.2546Figure 14
DC components of the flux linkages under differentLm1.Figure 15
τ
1 and τ2 components of the flux linkages under different Lm1.In regard to the inrush currents of the transformers, the sympathetic inrush is weakened and the initial inrush decays more slowly due to the saturation of T1.In conclusion, the series transformers interact with each other due to the transformer saturation. The initial inrush decays more quickly due to the saturation of the energized transformer, and the sympathetic inrush is weakened due to saturation of the in-service transformer. The saturation of the transformer resists further saturation, which can be thought of as a self-regulating mechanism. On the other hand, the initial inrush decays more slowly due to saturation of the in-service transformer, and the sympathetic inrush is strengthened due to saturation of the energized transformer. The saturation of one transformer strengthens the other transformer’s saturation, which can be explained in terms of the energy balance of the two transformers.


## 5. Influencing Factors of Sympathetic Inrush

External conditions for the series transformers that have been connected to the power system vary greatly during sympathetic inrush. The switching angle and residual flux of the energized transformer vary widely, and their values depend on the switching time of the CB. Moreover, the system impedance and the load vary greatly under different operating conditions of the power system. Therefore, their impacts on the sympathetic inrush are investigated in this section.


### 5.1. Switching Angle

For a single energized transformer without residual flux, the flux linkage has a maximum DC component when the CB is closed atφ = 0° because the prospective flux linkage is at the maximum, and the existing flux linkage cannot be changed suddenly. The flux linkage has a minimum DC component when the CB is closed atφ = 90°.Similarly, the largest initial and sympathetic inrush currents are generated atφ = 0° and the smallest are generated atφ = 90° without residual flux. By taking differentφ in (2), the DC components of the flux linkages are obtained and are shown in Figure 16. |b1|, |b2|, |c1|, and |c2| all decrease whenφ varies from 0° to 90°. Therefore, the magnetizing inrush and sympathetic inrush are both weakened asφ increases and will cease at φ = 90°.Figure 16
Variation of the DC components of the flux linkages with the switching angle.A simulation model (shown in Figure7) is constructed in PSCAD. The inrush currents can be represented by the differential currents, as shown in Figure 17, which are normalized to the 6 kV voltage side. A Y/Y energized transformer is used in this case instead of the Δ/Y transformer because the variation feature under different switching angles can be observed more clearly using a Y/Y transformer. It is observed that the variation of the inrush currents is the same as the variation of the DC components of the flux linkages.Figure 17
Variation of the initial inrush and sympathetic inrush with the switching angle.


### 5.2. Residual Flux

The residual flux of the energized transformer is another important factor that influences the inrush phenomenon. The difference between the residual flux linkage and prospective flux linkage of the energized transformer determines the amplitude of the initial inrush and sympathetic inrush. If the difference is large, the initial inrush and sympathetic inrush are also large. As shown in Figure18, the DC components of flux linkages vary with the residual flux (φ = 0°). When the CB is closed atφ = 0°, the prospective flux linkage is approximately equal to -Um/ω. The DC components are the largest when the residual flux is Um/ω and are zero when the residual flux is -Um/ω. A simulation case is conducted in Real Time Digital Simulator (RTDS) to further demonstrate the analysis. The inrush currents are shown in Figure 19.Figure 18
Variation of the DC components of the flux linkages with residual flux.Figure 19
Variation of the initial inrush and sympathetic inrush with residual flux.


### 5.3. System Impedance

The system impedanceZS=|RS+jωLS| also has a significant effect on the inrush currents. The DC components of the flux linkages are obtained from (2) and are shown in Figure 20. ψ1DC has larger and earlier peak values as ZS increases. A larger ZS is more beneficial for the generation of the sympathetic inrush. ψ2DC is smaller in the main decay process as ZS increases, particularly before t = 8 s. The system impedance has a greater effect on the sympathetic inrush than the initial inrush.Figure 20
Variation of the DC components of the flux linkages with system impedance.The inrush currents varying with the system impedance by simulation conducted in PSCAD are shown in Figure21. The larger ZS weakens the initial inrush and causes the sympathetic inrush to have a larger and earlier peak value, validating the theoretical analysis.Figure 21
Variation of the initial inrush and sympathetic inrush with system impedance.


### 5.4. Load of the In-Service Transformer

The load of the in-service transformer is also an important factor, especially for series transformers. The load adds an impedance branch to the right side ofL1σ. The added branch restrains the influence of the initial inrush on the in-service transformer. Although the flux linkage expressions including the load parameter have not been obtained, a simulation can be conducted to analyze the effect of the load. As shown in Figure 22, the initial inrush and sympathetic inrush decrease with increasing load, and the sympathetic inrush may disappear under a heavy load.Figure 22
Variation of the initial inrush and sympathetic inrush with the load.In contrast, the maximum amplitude of the sympathetic inrush generated in the parallel transformers atPL = 30 MW exhibits a 15% decrease relative to that at PL = 0 MW. Here, the type of parallel transformers in the simulation is the same as T1. The results show that the influence of the load on the series transformers is more significant than that on the parallel transformers.


## 5.1. Switching Angle

For a single energized transformer without residual flux, the flux linkage has a maximum DC component when the CB is closed atφ = 0° because the prospective flux linkage is at the maximum, and the existing flux linkage cannot be changed suddenly. The flux linkage has a minimum DC component when the CB is closed atφ = 90°.Similarly, the largest initial and sympathetic inrush currents are generated atφ = 0° and the smallest are generated atφ = 90° without residual flux. By taking differentφ in (2), the DC components of the flux linkages are obtained and are shown in Figure 16. |b1|, |b2|, |c1|, and |c2| all decrease whenφ varies from 0° to 90°. Therefore, the magnetizing inrush and sympathetic inrush are both weakened asφ increases and will cease at φ = 90°.Figure 16
Variation of the DC components of the flux linkages with the switching angle.A simulation model (shown in Figure7) is constructed in PSCAD. The inrush currents can be represented by the differential currents, as shown in Figure 17, which are normalized to the 6 kV voltage side. A Y/Y energized transformer is used in this case instead of the Δ/Y transformer because the variation feature under different switching angles can be observed more clearly using a Y/Y transformer. It is observed that the variation of the inrush currents is the same as the variation of the DC components of the flux linkages.Figure 17
Variation of the initial inrush and sympathetic inrush with the switching angle.


## 5.2. Residual Flux

The residual flux of the energized transformer is another important factor that influences the inrush phenomenon. The difference between the residual flux linkage and prospective flux linkage of the energized transformer determines the amplitude of the initial inrush and sympathetic inrush. If the difference is large, the initial inrush and sympathetic inrush are also large. As shown in Figure18, the DC components of flux linkages vary with the residual flux (φ = 0°). When the CB is closed atφ = 0°, the prospective flux linkage is approximately equal to -Um/ω. The DC components are the largest when the residual flux is Um/ω and are zero when the residual flux is -Um/ω. A simulation case is conducted in Real Time Digital Simulator (RTDS) to further demonstrate the analysis. The inrush currents are shown in Figure 19.Figure 18
Variation of the DC components of the flux linkages with residual flux.Figure 19
Variation of the initial inrush and sympathetic inrush with residual flux.


## 5.3. System Impedance

The system impedanceZS=|RS+jωLS| also has a significant effect on the inrush currents. The DC components of the flux linkages are obtained from (2) and are shown in Figure 20. ψ1DC has larger and earlier peak values as ZS increases. A larger ZS is more beneficial for the generation of the sympathetic inrush. ψ2DC is smaller in the main decay process as ZS increases, particularly before t = 8 s. The system impedance has a greater effect on the sympathetic inrush than the initial inrush.Figure 20
Variation of the DC components of the flux linkages with system impedance.The inrush currents varying with the system impedance by simulation conducted in PSCAD are shown in Figure21. The larger ZS weakens the initial inrush and causes the sympathetic inrush to have a larger and earlier peak value, validating the theoretical analysis.Figure 21
Variation of the initial inrush and sympathetic inrush with system impedance.


## 5.4. Load of the In-Service Transformer

The load of the in-service transformer is also an important factor, especially for series transformers. The load adds an impedance branch to the right side ofL1σ. The added branch restrains the influence of the initial inrush on the in-service transformer. Although the flux linkage expressions including the load parameter have not been obtained, a simulation can be conducted to analyze the effect of the load. As shown in Figure 22, the initial inrush and sympathetic inrush decrease with increasing load, and the sympathetic inrush may disappear under a heavy load.Figure 22
Variation of the initial inrush and sympathetic inrush with the load.In contrast, the maximum amplitude of the sympathetic inrush generated in the parallel transformers atPL = 30 MW exhibits a 15% decrease relative to that at PL = 0 MW. Here, the type of parallel transformers in the simulation is the same as T1. The results show that the influence of the load on the series transformers is more significant than that on the parallel transformers.


## 6. Conclusion

This paper presents a thorough study of the sympathetic inrush for series transformers. Flux linkage expressions of series transformers, including the switching angle, residual flux, and system impedance factors, are derived using the Laplace transform. Flux linkages of series transformers consist of the AC and DC components when the sympathetic inrush occurs. The series transformers interact with each other during the sympathetic inrush. Saturation of each transformer resists further self-saturation and strengthens the other transformer’s saturation.The influencing factors of sympathetic inrush for series transformers, including the switching angle, residual flux, and system impedance, are analyzed and validated. The smaller switching angle (0°–90°), larger impedance, and lighter load are more beneficial for the generation of the sympathetic inrush. If the difference between the residual flux in the energized transformer and the integral of the energizing voltage is large, the initial inrush and sympathetic inrush are also large.



---

*Source: 9104707-2017-12-11.xml*