Analysis of the Electromechanical Characteristics of Power Transformer under Different Residual Fluxes

: When the electromagnetic transients occur in a power transformer, an inrush current is generated in its winding. The inrush current not only affects the performance of the transformer windings, but also affects the lifetime of the transformer. Many factors affect the inrush current, the most inﬂuential ones among which are the closing phase angle and the residual ﬂux. In this paper, a dry-type transformer simulation model is built to analyze the inﬂuence of the inrush current on the performance of transformer windings during no-load reclosing conditions. Firstly, the inrush current was generated in the transformer windings during the no-load reclosing operation under different residual ﬂuxes. Secondly, the ﬁeld-circuit coupling 3d ﬁnite element method is used to analyze the electromagnetic force at different locations of the transformer windings under the inﬂuence of different inrush currents. The results of winding structural parameter variations are obtained through electromagnetic-structural coupling simulation, and the electromagnetic forces are used as the input parameter for the structural analysis. Finally, the residual ﬂux is generated by controlling the opening and closing angle of the transformer through the phase-controlled switch, and the winding electromechanical characteristics are tested under different residual ﬂuxes. Finally, comparisons of the test and simulation results are drawn to verify the impact of the closing angle and residual ﬂux on inrush current and the winding deformation during the no-load reclosing conditions.


Introduction
When a power transformer is cut out of the system, a certain amount of residual flux will be generated due to the hysteresis characteristics of the core material. The residual flux will have a large impact on the stable operation of the transformer when the transformer is reclosed. The inrush current can be 6-10 times the rated current when electromagnetic transients occur in the reclosing process of the power transformer. Therefore, accurate measurement of the residual flux of the core is of great significance for analyzing the electromechanical characteristics of the winding under the inrush current.
To detect the residual flux in the transformer core, many algorithms and models have been proposed by domestic and foreign researchers. For example, Ref. [1] used the change in the permeability of local hysteresis lines to calculate the residual flux, but its calculation method is too complicated for application in engineering practice. Ref. [2] proposed a method for calculating the residual flux considering the transient magnetization current, and derived an empirical formula for the calculation of the residual flux magnetic density, the error of which can always be less than 10%. However, this method cannot determine the direction of the residual flux. Refs. [3,4] proposed a residual flux measurement method by establishing a known relationship between the residual flux and the current of a power transformer. They built a core detection circuit for ring transformers, processes the experiments to obtain transient measurement current signals, and substituted the results into the residual flux and current relationship equation. This method was able to obtain the residual flux in the core of power transformers more accurately.
For more results, many researchers have started their research from the perspective considering the residual flux's direction. For example, Ref. [5] proposed a residual flux measurement method considering the residual flux direction. The transient current of the coil under different voltages and residual flux densities was analyzed by means of both the principle calculation and simulation analysis, and an empirical formula for the residual flux was obtained. Finally, the results were verified by experiments. A method for interpreting the residual flux of a single-phase transformer by sequentially applying two Direct Current (DC) voltages of opposite polarity to the transformer windings and then detecting their corresponding transient response currents was proposed in [6]. The equation between the residual flux and the transient current difference was fitted using the data analysis method. Finally, the effectiveness of the method was verified through experiment. Although the above methods are able calculate the residual flux's value in the transformer core, they are difficult to implement in practice.
Excessive inrush current in the windings can seriously affect the electromechanical characteristics of the transformer, which may even in serious cases lead to the destruction of the transformer. Refs. [7,8] analyzed the electromagnetic forces generated by inrush currents and the impact they have on the transformer structure. Mechanical failures of the transformer usually manifest as the deformation and the displacement of the windings. When the deformation and displacement of the windings exceeds the threshold value, plastic deformation occurs, and the insulation capacity is gradually reduced. Furthermore, the failure resistance of the transformer is significantly reduced.
To analyze the effects of excessive inrush currents on transformers, structural optimization and simulation methods for analyzing magnetic fields of transformers have emerged. One of the most widely used computational analysis methods for electromagnetic as well as the structural simulation is the finite element method. Residual flux in the core could produce an inrush current which amplitude can far exceed the rated current. For example, in [9], a numerical analysis method was proposed to accurately calculate the inrush current of a three-phase three-arm core transformer considering the nonlinear magnetization curve of the iron core. Ref. [10] used Octave software to calculate the maximum value of the inrush current generated by a single-phase transformer. Meanwhile, the magnetic density, air-gap leakage, and winding losses of the transformer core under the maximum inrush current were analyzed using finite element simulation analysis. In Ref. [11], a method for calculating the inrush current of superconducting transformers was proposed, and the accuracy of the calculation results was verified by experiments. Refs. [12][13][14][15] analyzed the transient electromagnetic forces of large power transformer windings under shortcircuit currents. A three-dimensional transformer electromagnetic-structural coupling finite element model was established to analyze the dynamic electromagnetic force and winding mechanical structure interaction during short-circuit condition, and the simulation analysis results were compared with the experimental results. Ref. [16] qualitatively and quantitatively analyzed the mechanical power generated in the transformer winding under the action of inrush current by means of analytical and numerical methods. The analytical and computational results were compared and analyzed with the experimental results, which proved the validity of the method and the correctness of the conclusions. A method based on digital image processing of transformer winding deformation was proposed in Ref. [17]. Ref. [18] detected and classified different operating conditions of transformers by building a neural network classification model. The above studies focus on the deformation of transformer windings under excitation current and the ability of windings to resist fault current, with fewer studies examining the combined variation of the transformers under inrush current, with no clear conclusions being given.
Both the short-circuit current and the inrush current have a great impact on the transformer, and in some cases, the inrush current on the transformer winding can even produce a greater impact than the short-circuit current. Ref. [8] analyzed the short-circuit forces on the windings of a 50 KVA, dry-type transformer under short-circuit current using finite elements and verified the correctness of the simulation analysis by experiments. In Ref. [19], a single-phase, 50 MVA transformer model was established, and the model was analyzed using the finite element method to analyze the axial and radial forces on the transformer windings under the inrush current. In Ref. [20], the electromechanical stresses generated in the transformer windings by the inrush current and short-circuit current were studied by means of both numerical calculations and simulation analysis. In Ref. [21], the radial short-circuit forces generated on the transformer's low-voltage (LV) winding under short-circuit were studied. The relationship between the leakage flux in the transformer windings and the short-circuit electromagnetic force was analyzed. The radial short-circuit forces were calculated, and the ability of the LV winding to withstand short-circuit currents was evaluated according to International Electrotechnical Commission (IEC) standards and short-circuit experiments. Ref. [22] used the electromagnetic-thermal-fluid coupling method to analyze the effects of inrush currents on transformer windings. Ref. [23] used the field-circuit coupling method to calculate the leakage magnetic field and the short-circuit electromagnetic force under short-circuit conditions. According to the actual parameters, the flexural of the windings was analyzed to study the axial stability of the windings. In Ref. [24], the field-circuit coupling method was used to compare the leakage and the electromagnetic forces of transformer low-voltage windings under different short-circuit conditions. In Ref. [25], a simulation model of transformer winding loosening, winding deformation and winding insulation failure was developed using a simulation model of a transformer in normal operation, and the vibration characteristics of winding faults were analyzed. The above-mentioned studies mainly focus on the analysis of the impact of shortcircuit current on transformer windings from a structural perspective. Studies on the effect of the electromagnetic transients generated by the inrush current on the electromechanical characteristics of the transformer are fewer, and most of the conclusions are vague, while at the same time lacking experimental verification.
With increasing transformer capacity, the impact of the inrush current on the transformer becomes more and more significant. In particular, in the case of no-load reclosing, the impact of inrush current on the transformer is more prominent. In this case, the analysis of the inrush current during the no-load reclosing transient, as well as the electromechanical characteristics of the windings under the inrush current, are of great importance, and may extend the lifetime of the transformer and allow the diagnosis of transformer fault. Therefore, this paper provides a comprehensive analysis of the transformer magnetic field changes and leakage flux distribution at the moment of transformer reclosing. At the same time, the electromagnetic-structural field coupling method is used to analyze the impact of the inrush current by the transformer re-closing condition.
Inrush current is generated at the instant of no-load reclosing of the transformer, and is mainly a result of the closing phase angle and residual flux. Currently, most of the research on residual flux is focused on measurement and suppression, and most of the methods are difficult to implement in operating transformers. Therefore, this paper uses the phase-controlled switching technique to control and calculate the residual fluxes of the transformer core under different closing phase angles, as well as to analyze the peak excitation inrush currents and their effect on the transformer windings generated under different residual fluxes.
In summary, this paper presents a comprehensive analysis of the electromechanical characteristic of power transformer windings with different residual fluxes and reclosing angles. Firstly, a simulation circuit of the power transformer is established to analyze and obtain the inrush current under the no-load reclosing condition of the transformer with different residual fluxes. Secondly, the "field-path" coupling method and finite element method are combined to analyze the electromagnetic force of the transformer winding under the inrush current corresponding to different residual fluxes. After that, the power transformer electromagnetic-structural coupling analysis model is established to analyze the deformation of the winding structure. Then, based on the simulation results, the parametric relationships between residual flux and closing phase angle and electromagnetic force and winding deformation are obtained. Finally, the validity of the simulation method and the correctness of the results are verified by experiments.

Generation of Residual Flux Magnetism
Power transformer cores are mainly composed of silicon steel sheets with hysteresis characteristics, and the hysteresis characteristics exhibited by transformer cores differ depending on the core material. The hysteresis characteristics of the transformer core are mainly manifested when the external magnetization conditions are increased, i.e., after the transformer is energized, there will be a magnetization curve and a hysteresis return, i.e., the magnetization and demagnetization paths are different. Figure 1 shows the partitioning of the magnetization curve of the core material after the transformer is energized, where region I denotes the starting region, region II denotes the Rayleigh region, region III denotes the nonlinear steep region, and region VI denotes the converging saturation region. As can be seen from Figure 1, once the core starts to be magnetized, even if the magnetic field strength H decreases to zero, the magnetic induction strength B in the core does not drop to zero, and the magnetic induction strength value in the core at this time is the residual flux B r . obtain the inrush current under the no-load reclosing condition of the transformer with different residual fluxes. Secondly, the "field-path" coupling method and finite element method are combined to analyze the electromagnetic force of the transformer winding under the inrush current corresponding to different residual fluxes. After that, the power transformer electromagnetic-structural coupling analysis model is established to analyze the deformation of the winding structure. Then, based on the simulation results, the parametric relationships between residual flux and closing phase angle and electromagnetic force and winding deformation are obtained. Finally, the validity of the simulation method and the correctness of the results are verified by experiments.

Generation of Residual Flux Magnetism
Power transformer cores are mainly composed of silicon steel sheets with hysteresis characteristics, and the hysteresis characteristics exhibited by transformer cores differ depending on the core material. The hysteresis characteristics of the transformer core are mainly manifested when the external magnetization conditions are increased, i.e., after the transformer is energized, there will be a magnetization curve and a hysteresis return, i.e., the magnetization and demagnetization paths are different. Figure 1 shows the partitioning of the magnetization curve of the core material after the transformer is energized, where region I denotes the starting region, region II denotes the Rayleigh region, region III denotes the nonlinear steep region, and region VI denotes the converging saturation region. As can be seen from Figure 1, once the core starts to be magnetized, even if the magnetic field strength H decreases to zero, the magnetic induction strength B in the core does not drop to zero, and the magnetic induction strength value in the core at this time is the residual flux Br. Usually, the magnetic field strength H when the transformer is powered off varies, the resulting hysteresis return will be different, and the residual flux Br will be different. The degree of saturation of the core during power failure will also affect the residual flux Br. Therefore, it is possible to estimate the magnitude and direction of the residual flux Br in the core from the magnetic field strength H and the hysteresis return at the moment of power failure of the transformer. Usually, the magnetic field strength H when the transformer is powered off varies, the resulting hysteresis return will be different, and the residual flux B r will be different. The degree of saturation of the core during power failure will also affect the residual flux B r . Therefore, it is possible to estimate the magnitude and direction of the residual flux B r in the core from the magnetic field strength H and the hysteresis return at the moment of power failure of the transformer.

Calculation of Residual Flux Magnetism
The phase angle, amplitude, and direction of the residual flux at the moment of transformer closing are the main factors affecting the size of the inrush current amplitude and its duration under no-load reclosing. However, the peak of the inrush current is mainly affected by the phase angle θ and the remaining flux φ r . Under the unsaturated state of the core flux φ, the magnetization curve of the core basically conforms to the linear variation. When φ approaches the saturation flux φ s , the magnetization curve gradually becomes nonlinear, and the excitation current increases sharply as the core saturation increases, at which time the excitation current can even reach thousands of times the rated no-load current, i.e., the inrush current.
According to the full current law, it can be obtained: where l is the transformer magnetic circuit length (m); N is the number of turns of the excitation winding; H is the magnetic field strength (A/m); i(t) is the current in the excitation winding (A). Because the excitation winding direct current resistance R is very small and the leakage resistance is very small, the calculation ignores the direct current resistance voltage drop as well as the leakage resistance voltage drop, and according to the law of electromagnetic induction, the following can be obtained: where S is the effective cross-sectional area of the transformer (m 2 ); B(t) is the average magnetic density in the core (T); φ is the main magnetic flux (Wb). Based on the voltage and current data at the moment of the controller control break, the corresponding B and H values are calculated to obtain the B-H curve corresponding to the iron core. Additionally, the ratio of the maximum magnetic density B m formed by the hysteresis return to the residual flux magnetic density B r is determined according to the different operating voltages, and then B r is obtained from the measured B m .

Generation and Calculation of Excitation Inrush Current
The transformer needs to go through a transient process after reclosing before it can reach a steady state, i.e., the decay process of the excitation current. Equation (3) shows the relationship between voltage, current and flux in a transformer.
where V m is the peak voltage (V), θ is the voltage switching angle when the transformer is closed ( • ), i 0 (t) is the instantaneous current generated at the time of reclosing (A), φ m (t) is the instantaneous flux which contains the main flux and the leakage flux (Wb), R is the resistance of the windings in the transformer (Ω). Equation (4) is the value of the instantaneous flux when the transformer is energized instantaneously.
where the peak value of the steady-state flux of φ p (Wb), φ r is the value of the residual flux (Wb), L is the inductance of the excitation coil (H). When θ = 0 • and the φ r of the residual flux is positive, the flux waveform is shown in Figure 2.
From Figure 2, the magnetic flux φ in the transformer core is kept at a low level when the residual flux equals to zero. When residual flux is present, the initial value of the magnetic flux in the core will be increased, which makes it easier to saturate the core. Figure 3 shows a typical inrush current waveform at the moment of transformer energization.  From Figure 2, the magnetic flux ϕ in the transformer core is kept at a low lev the residual flux equals to zero. When residual flux is present, the initial value of th netic flux in the core will be increased, which makes it easier to saturate the core. Figure 3 shows a typical inrush current waveform at the moment of transfor ergization. When calculating the peak inrush current, it is first necessary to determine t ration flux in the core, which can be determined by Equation (5).
where S is the cross-sectional area of the iron core (cm 2 ). The leakage flux ϕair in the air gap can be calculated by Equation (6)    From Figure 2, the magnetic flux ϕ in the transformer core is kept at a low level when the residual flux equals to zero. When residual flux is present, the initial value of the magnetic flux in the core will be increased, which makes it easier to saturate the core. Figure 3 shows a typical inrush current waveform at the moment of transformer energization. When calculating the peak inrush current, it is first necessary to determine the saturation flux in the core, which can be determined by Equation (5).
where S is the cross-sectional area of the iron core (cm 2 ). The leakage flux ϕair in the air gap can be calculated by Equation (6) where H is the magnetic field strength (A/m). Equation (7) is the average cross-sectional area of the coil winding Sm.
where Dm average winding diameter (m).
In the case of magnetic saturation of the iron core, the system reactance value Xs can be calculated by Equation (8). When calculating the peak inrush current, it is first necessary to determine the saturation flux in the core, which can be determined by Equation (5).
where S is the cross-sectional area of the iron core (cm 2 ). The leakage flux φ air in the air gap can be calculated by Equation (6) (Wb).
where H is the magnetic field strength (A/m). Equation (7) is the average cross-sectional area of the coil winding S m .
where D m average winding diameter (m).
In the case of magnetic saturation of the iron core, the system reactance value X s can be calculated by Equation (8). where µ 0 is the air magnetic permeability (H/m). f is the power supply frequency (Hz). Equation (9) represents the switching angle.
where K 1 is the saturation angle correction factor. B s is the saturation magnetic flux density (T). B r is the residual flux magnetic flux density (T). K 1 , B s, and B r are equal to 0.9, 1.25 B mp , and 0.8 B mp . The peak value of magnetic induction (B mp ) in the transformer core can be calculated by Equation (10).
where V is the voltage value when the transformer is operating (V). S core is the crosssectional area of the iron core. According to the formula derived from our previous study [26], Equation (11) is used for calculating the peak inrush current.
where i 0max is the peak inrush current generated by the single-phase transformer (A). K 2 is the peak inrush current correction factor, whose value is 1.15.

Calculation of Winding Electromagnetic Force
The moment the transformer is energized, the inrush current and the leakage field interact to produce a large instantaneous electromagnetic force, and Equation (12) is the formula for calculating the electromagnetic force.
where F is the bulk density vector of electromagnetic force (N/m 3 ), J is the current density vector (A/m 2 ), B is the leakage flux density (T).
To analyze the impact of electromagnetic force on the transformer winding more accurately, the electromagnetic force is usually decomposed into radial force F r and axial force F a to analyze its impact on the transformer winding in detail.

Radial Force Analysis and Calculation
The radial force F r generated by the inrush current is a major source of mechanical stress in transformer high-voltage winding during excitation. The main reason for this situation is that the electromagnetic force generated by the inrush current is not generated in the primary side winding and secondary side winding at the same time as the shortcircuit current, but rather is mostly concentrated in the high-voltage winding. Therefore, the analysis of the electromechanical characteristics of the transformer winding under no-load reclosing is mainly focused on high-voltage winding. Figure 4 shows the radial electromagnetic force generated on the windings of the transformer.
Formula (13) is the formula for calculating the radial force F r on the transformer winding.  Formula (13) is the formula for calculating the radial force Fr on the transforme winding.
Equation (14) shows the average value of the electromagnetic force on the winding.
where Fr is the radial electromagnetic force (N). Frmed is the radial average electromagneti force (N). Iinr is the peak inrush current (A). h is the winding height (m). Dm is the averag diameter of the winding (m). N is the number of turns of the winding coil. Figure 5 shows two forms of deformation-forced and free flexion-when the trans former winding is subjected to excessive radial forces. The transformer winding is sub jected to inward pressure under the influence of radial force. When the pressure valu exceeds the mechanical strength of the low-voltage winding support structure, the wind ing will produce two different forms of deformation: free flexion and forced flexion Forced flexion refers to the deformation of transformer windings along the direction o electromagnetic force under the influence of strong electromagnetic force, which is mainl related to the strength of electromagnetic force and winding support structure. Free flex ion means that the transformer winding has a deformation along the direction of the non electromagnetic force under the influence of two electromagnetic forces with an angle Factors such as the strength of the electromagnetic force, the angle of the two electromag netic forces and the structure of the support bar may have an effect on the free flexion.  Equation (14) shows the average value of the electromagnetic force on the winding.
where F r is the radial electromagnetic force (N). F rmed is the radial average electromagnetic force (N). I inr is the peak inrush current (A). h is the winding height (m). D m is the average diameter of the winding (m). N is the number of turns of the winding coil. Figure 5 shows two forms of deformation-forced and free flexion-when the transformer winding is subjected to excessive radial forces. The transformer winding is subjected to inward pressure under the influence of radial force. When the pressure value exceeds the mechanical strength of the low-voltage winding support structure, the winding will produce two different forms of deformation: free flexion and forced flexion. Forced flexion refers to the deformation of transformer windings along the direction of electromagnetic force under the influence of strong electromagnetic force, which is mainly related to the strength of electromagnetic force and winding support structure. Free flexion means that the transformer winding has a deformation along the direction of the non-electromagnetic force under the influence of two electromagnetic forces with an angle. Factors such as the strength of the electromagnetic force, the angle of the two electromagnetic forces and the structure of the support bar may have an effect on the free flexion. Formula (13) is the formula for calculating the radial force Fr on the transformer winding.
Equation (14) shows the average value of the electromagnetic force on the winding.
where Fr is the radial electromagnetic force (N). Frmed is the radial average electromagnetic force (N). Iinr is the peak inrush current (A). h is the winding height (m). Dm is the average diameter of the winding (m). N is the number of turns of the winding coil. Figure 5 shows two forms of deformation-forced and free flexion-when the transformer winding is subjected to excessive radial forces. The transformer winding is subjected to inward pressure under the influence of radial force. When the pressure value exceeds the mechanical strength of the low-voltage winding support structure, the winding will produce two different forms of deformation: free flexion and forced flexion. Forced flexion refers to the deformation of transformer windings along the direction of electromagnetic force under the influence of strong electromagnetic force, which is mainly related to the strength of electromagnetic force and winding support structure. Free flexion means that the transformer winding has a deformation along the direction of the nonelectromagnetic force under the influence of two electromagnetic forces with an angle. Factors such as the strength of the electromagnetic force, the angle of the two electromagnetic forces and the structure of the support bar may have an effect on the free flexion.

Analysis and Calculation of Axial Forces
At the moment of no-load reclosing of the transformer, there will be both radial electromagnetic force component F r and axial electromagnetic force component F a in the winding. Figure 6 shows the axial electromagnetic force component F a generated on high-voltage winding under the inrush current. At the moment of no-load reclosing of the transformer, there will be b tromagnetic force component Fr and axial electromagnetic force component ing. Figure 6 shows the axial electromagnetic force component Fa generated age winding under the inrush current.
The radial leakage generated at the end of the winding interacts with rent to generate axial forces. The total axial force generated at the center o can be found by Equation (15) [19][20][21].
where Fa is the axial compression force (F/m). d is the radial thickness of the (m).

Matlab/Simulink Circuit Simulation Analysis
In this paper, a single-phase transformer circuit model is built, as show

Analysis and Calculation of Axial Forces
At the moment of no-load reclosing of the transformer, there will be both radial electromagnetic force component Fr and axial electromagnetic force component Fa in the winding. Figure 6 shows the axial electromagnetic force component Fa generated on high-voltage winding under the inrush current.  The radial leakage generated at the end of the winding interacts with the inrush current to generate axial forces. The total axial force generated at the center of the winding can be found by Equation (15) [19][20][21].
where Fa is the axial compression force (F/m). d is the radial thickness of the energized coil (m).

Matlab/Simulink Circuit Simulation Analysis
In this paper, a single-phase transformer circuit model is built, as shown in Figure 8.
The radial leakage generated at the end of the winding interacts with the inrush current to generate axial forces. The total axial force generated at the center of the winding can be found by Equation (15) [19][20][21].
where F a is the axial compression force (F/m). d is the radial thickness of the energized coil (m).

Matlab/Simulink Circuit Simulation Analysis
In this paper, a single-phase transformer circuit model is built, as shown in Figure 8. The research presented in this paper is a continuation of our previous work [26]. The parameters of the test transformer are shown in Table 1. To ensure that the simulation results are as close as possible to the actual values, the parameters of the prototype need to be transformed into those of the specific module, and the hysteresis characteristics of the core material, i.e., the i-φ curve, must be taken into account. s 2021, 14, x FOR PEER REVIEW 10 of 22 The research presented in this paper is a continuation of our previous work [26]. The parameters of the test transformer are shown in Table 1. To ensure that the simulation results are as close as possible to the actual values, the parameters of the prototype need to be transformed into those of the specific module, and the hysteresis characteristics of the core material, i.e., the i-ϕ curve, must be taken into account. To consider the hysteresis characteristics of the core, the magnetization curve of the transformer core needs to be set. In this paper, the hysteresis design function of the Powergui module is used to design and modify the hysteresis return of the core material to complete the setting of the hysteresis characteristics of the transformer core. Figure 9 shows the magnetization curve of the transformer core after modifying the setting parameters.   To consider the hysteresis characteristics of the core, the magnetization curve of the transformer core needs to be set. In this paper, the hysteresis design function of the Powergui module is used to design and modify the hysteresis return of the core material to complete the setting of the hysteresis characteristics of the transformer core. Figure 9 shows the magnetization curve of the transformer core after modifying the setting parameters.
parameters of the test transformer are shown in Table 1. To ens results are as close as possible to the actual values, the parameter to be transformed into those of the specific module, and the hys the core material, i.e., the i-ϕ curve, must be taken into account. To consider the hysteresis characteristics of the core, the ma transformer core needs to be set. In this paper, the hysteresis desig gui module is used to design and modify the hysteresis return of t plete the setting of the hysteresis characteristics of the transform the magnetization curve of the transformer core after modifying t    Figure 11 shows the corresponding inrush current amplitude under different residual fluxes and closing angles.  Figure 11 shows that the residual flux in the transformer core as well as the angle during no-load reclosing have an effect on the inrush current. The maximum amplitude of the inrush current can be obtained when other factors remain constant and the closing phase angle θ = 0°. When the closing angle is θ = 90°, the minimum amplitude of the inrush current is generated. While the reclosing angles are the same, with the increase of residual flux, the amplitude of the inrush current generated in the transformer winding is also gradually increasing, and its amplitude can reach up to 5165 A. Therefore, the inrush current in the windings is generated by the combined effect of residual flux and closing angle during transformer operation.

Windings Electromagnetic Simulation Analysis
A single-phase transformer is used in order to analyze its leakage, electromagnetic force and winding deformation. Simulation analysis shows that under the combined effect of inrush current and leakage, a large electromagnetic force is generated on the winding. Therefore, the accurate measurement of magnetic leakage is very important. To ensure the accuracy of the simulation results, the measurement points of magnetic leakage are   Figure 11 shows the corresponding inrush current amplitude under different r ual fluxes and closing angles.  Figure 11 shows that the residual flux in the transformer core as well as the a during no-load reclosing have an effect on the inrush current. The maximum ampli of the inrush current can be obtained when other factors remain constant and the clo phase angle θ = 0°. When the closing angle is θ = 90°, the minimum amplitude of the in current is generated. While the reclosing angles are the same, with the increase of resi flux, the amplitude of the inrush current generated in the transformer winding is gradually increasing, and its amplitude can reach up to 5165 A. Therefore, the inrush rent in the windings is generated by the combined effect of residual flux and closing a during transformer operation.

Windings Electromagnetic Simulation Analysis
A single-phase transformer is used in order to analyze its leakage, electromag force and winding deformation. Simulation analysis shows that under the combined e of inrush current and leakage, a large electromagnetic force is generated on the wind Therefore, the accurate measurement of magnetic leakage is very important. To en the accuracy of the simulation results, the measurement points of magnetic leakage  Figure 11 shows that the residual flux in the transformer core as well as the angle during no-load reclosing have an effect on the inrush current. The maximum amplitude of the inrush current can be obtained when other factors remain constant and the closing phase angle θ = 0 • . When the closing angle is θ = 90 • , the minimum amplitude of the inrush current is generated. While the reclosing angles are the same, with the increase of residual flux, the amplitude of the inrush current generated in the transformer winding is also gradually increasing, and its amplitude can reach up to 5165 A. Therefore, the inrush current in the windings is generated by the combined effect of residual flux and closing angle during transformer operation.

Windings Electromagnetic Simulation Analysis
A single-phase transformer is used in order to analyze its leakage, electromagnetic force and winding deformation. Simulation analysis shows that under the combined effect of inrush current and leakage, a large electromagnetic force is generated on the winding. Therefore, the accurate measurement of magnetic leakage is very important. To ensure the accuracy of the simulation results, the measurement points of magnetic leakage are selected as comprehensively as possible. Figure 12 shows a schematic diagram of the location of the measurement point selection for air gap leakage of the transformer winding.
Energies 2021, 14, x FOR PEER REVIEW 12 of 2 selected as comprehensively as possible. Figure 12 shows a schematic diagram of the lo cation of the measurement point selection for air gap leakage of the transformer winding  From Figure 13, the amplitude of radial leakage at measurement points A and D, and E, C and F are the same, but in opposite directions, while the direction and amplitud  selected as comprehensively as possible. Figure 12 shows a schematic diagram of the location of the measurement point selection for air gap leakage of the transformer winding.  From Figure 13, the amplitude of radial leakage at measurement points A and D, B and E, C and F are the same, but in opposite directions, while the direction and amplitude winding. The maximum axial magnetic density B a = 0.1157 T and the maximum radiation magnetic density B r = 0.0623 T.

Simulation Analysis of Winding Structure
When the external conditions are basically the same, the inrush current amplitude generated during no-load reclosing is higher than the inrush current amplitude generated under other circumstances. Meanwhile, most of the inrush current generated by no-load reclosing is concentrated in the high-voltage winding. Therefore, the paper focuses on the structural characteristics of transformer primary windings under different residual fluxes no-load reclosing conditions. Under normal operating conditions, the electromagnetic forces and leakage fields generated in the transformer windings are relatively small, so the winding support structure is fully capable of withstanding electromagnetic stresses. However, under the transient phenomenon of a sudden increase in inrush current, the leakage field due to high currents also reaches considerable values, so that the transformer is damaged when the electromagnetic force is greater than the stress that the winding can withstand.
When the winding height is the same, the magneto-motive force in the transformer is uniformly distributed, and the radial leakage at both ends interacts with the energized winding to produce axial force, with the maximum value at both ends of the winding, squeezing the winding. The radial force is generated by the interaction between the axial leakage and the energized winding, with the maximum value in the middle of the winding and decreasing toward the ends of the winding, stretching the winding, and usually generating a circumferential stress on the primary side windings. Figure 14 shows a schematic diagram of the positions of the primary windings in the transformer model. of axial leakage remain the same. The maximum value of axial leakage occurs in the middle of the winding, while the maximum value of radial leakage occurs at two ends of the winding. The maximum axial magnetic density Ba = 0.1157 T and the maximum radiation magnetic density Br = 0.0623 T.

Simulation Analysis of Winding Structure
When the external conditions are basically the same, the inrush current amplitude generated during no-load reclosing is higher than the inrush current amplitude generated under other circumstances. Meanwhile, most of the inrush current generated by no-load reclosing is concentrated in the high-voltage winding. Therefore, the paper focuses on the structural characteristics of transformer primary windings under different residual fluxes no-load reclosing conditions. Under normal operating conditions, the electromagnetic forces and leakage fields generated in the transformer windings are relatively small, so the winding support structure is fully capable of withstanding electromagnetic stresses. However, under the transient phenomenon of a sudden increase in inrush current, the leakage field due to high currents also reaches considerable values, so that the transformer is damaged when the electromagnetic force is greater than the stress that the winding can withstand.
When the winding height is the same, the magneto-motive force in the transformer is uniformly distributed, and the radial leakage at both ends interacts with the energized winding to produce axial force, with the maximum value at both ends of the winding, squeezing the winding. The radial force is generated by the interaction between the axial leakage and the energized winding, with the maximum value in the middle of the winding and decreasing toward the ends of the winding, stretching the winding, and usually generating a circumferential stress on the primary side windings. Figure 14 shows a schematic diagram of the positions of the primary windings in the transformer model.

Simulation Analysis of Winding Structure
When the external conditions are basically the same, the inrush current am generated during no-load reclosing is higher than the inrush current amplitude ge under other circumstances. Meanwhile, most of the inrush current generated by reclosing is concentrated in the high-voltage winding. Therefore, the paper focuse structural characteristics of transformer primary windings under different residua no-load reclosing conditions. Under normal operating conditions, the electromagnetic forces and leakag generated in the transformer windings are relatively small, so the winding suppo ture is fully capable of withstanding electromagnetic stresses. However, under th ent phenomenon of a sudden increase in inrush current, the leakage field due currents also reaches considerable values, so that the transformer is damaged w electromagnetic force is greater than the stress that the winding can withstand.
When the winding height is the same, the magneto-motive force in the tran is uniformly distributed, and the radial leakage at both ends interacts with the en winding to produce axial force, with the maximum value at both ends of the w squeezing the winding. The radial force is generated by the interaction between t leakage and the energized winding, with the maximum value in the middle of th ing and decreasing toward the ends of the winding, stretching the winding, and generating a circumferential stress on the primary side windings. Figure 14 shows a schematic diagram of the positions of the primary winding transformer model.     Figure 15, it can be obtained that the electromagnetic force at both ends of the winding is greater than that at the middle, and the maximum value of the electromagnetic force reaches 2.5707 × 10 6 N/m 3 . Figure 16 shows the stress and displacement clouds of the primary windings when the residual flux magnetization B r = 1.24 T and the closing angle θ = 0 • . From Figure 15, it can be obtained that the electromagnetic force at both ends of the winding is greater than that at the middle, and the maximum value of the electromagnetic force reaches 2.5707 × 10 6 N/m 3 . Figure 16 shows the stress and displacement clouds of the primary windings when the residual flux magnetization Br = 1.24T and the closing angle θ = 0°. The axial force is due to the radial magnetic field. The maximum stress of the winding reaches up to 1.7334 × 10 7 Pa under the excitation of the inrush current, as can be seen from Figure 16a. Figure 16b shows that the most significant deformation is in the middle of the winding, with a maximum deformation of 4.6187 × 10 −2 mm. Based on Figure 7, the analysis of the impact of the axial force Fa on the deformation of the winding, combined with Figure 16, that the axial electromagnetic force generated under the influence of the excitation surge will squeeze the winding, result in an intensification of the winding amp-turn imbalance, which will lead to the collapse of the transformer.
To ensure the accuracy of the winding structure analysis as much as possible, the transformer windings are partitioned as shown in Figure 17. In addition, the electromagnetic force at the selected measurement points on each partition winding is measured under the conditions of residual flux Br = 1.24T and closing angle θ = 0°.   The axial force is due to the radial magnetic field. The maximum stress of the winding reaches up to 1.7334 × 10 7 Pa under the excitation of the inrush current, as can be seen from Figure 16a. Figure 16b shows that the most significant deformation is in the middle of the winding, with a maximum deformation of 4.6187 × 10 −2 mm. Based on Figure 7, the analysis of the impact of the axial force F a on the deformation of the winding, combined with Figure 16, that the axial electromagnetic force generated under the influence of the excitation surge will squeeze the winding, result in an intensification of the winding amp-turn imbalance, which will lead to the collapse of the transformer.
To ensure the accuracy of the winding structure analysis as much as possible, the transformer windings are partitioned as shown in Figure 17. In addition, the electromagnetic force at the selected measurement points on each partition winding is measured under the conditions of residual flux B r = 1.24 T and closing angle θ = 0 • . From Figure 15, it can be obtained that the electromagnetic force at both ends of the winding is greater than that at the middle, and the maximum value of the electromagnetic force reaches 2.5707 × 10 6 N/m 3 . Figure 16 shows the stress and displacement clouds of the primary windings when the residual flux magnetization Br = 1.24T and the closing angle θ = 0°. The axial force is due to the radial magnetic field. The maximum stress of the winding reaches up to 1.7334 × 10 7 Pa under the excitation of the inrush current, as can be seen from Figure 16a. Figure 16b shows that the most significant deformation is in the middle of the winding, with a maximum deformation of 4.6187 × 10 −2 mm. Based on Figure 7, the analysis of the impact of the axial force Fa on the deformation of the winding, combined with Figure 16, that the axial electromagnetic force generated under the influence of the excitation surge will squeeze the winding, result in an intensification of the winding amp-turn imbalance, which will lead to the collapse of the transformer.
To ensure the accuracy of the winding structure analysis as much as possible, the transformer windings are partitioned as shown in Figure 17. In addition, the electromagnetic force at the selected measurement points on each partition winding is measured under the conditions of residual flux Br = 1.24T and closing angle θ = 0°.     Figure 18 shows that the value of the radial electromagnetic force is significantly higher in the middle of the primary windings than both ends, and the maximum value of the axial electromagnetic force at both ends is higher than the middle position of the primary windings, under inrush current excitation. The primary side winding is subjected to excessive radial electromagnetic force, Fr, which will not only lead to an increase in the air gap between the primary side and secondary side windings, but may also damage the insulation of the primary side winding, resulting in a short circuit between turns of the transformer. At the same time, the radial electromagnetic force in the middle of the winding is greater than that at the ends, thus making the structure of the primary windings of a transformer operating for a long time close to an ellipsoidal shape. Under the influence of the inrush current, the axial electromagnetic force generated at both ends of the primary windings will squeeze the winding to the middle of the winding, making the unbalance between the turns of the winding intensify, thus increasing the electromagnetic force on the winding and reducing the transformer life. In general, the axial electromagnetic force Fa generally has a more pronounced effect on the end of the winding, and the middle of the winding is more influenced by the radial electromagnetic force Fr.

Phase Control Switch Principle
To analyze the electromechanical characteristics of the windings under different residual fluxes, the transformer is controlled by a compound switch composed of thyristors to turn on and off. Figure 19 shows a schematic diagram of the phase-controlled switch.  Figure 18 shows that the value of the radial electromagnetic force is significantly higher in the middle of the primary windings than both ends, and the maximum value of the axial electromagnetic force at both ends is higher than the middle position of the primary windings, under inrush current excitation. The primary side winding is subjected to excessive radial electromagnetic force, F r , which will not only lead to an increase in the air gap between the primary side and secondary side windings, but may also damage the insulation of the primary side winding, resulting in a short circuit between turns of the transformer. At the same time, the radial electromagnetic force in the middle of the winding is greater than that at the ends, thus making the structure of the primary windings of a transformer operating for a long time close to an ellipsoidal shape. Under the influence of the inrush current, the axial electromagnetic force generated at both ends of the primary windings will squeeze the winding to the middle of the winding, making the unbalance between the turns of the winding intensify, thus increasing the electromagnetic force on the winding and reducing the transformer life. In general, the axial electromagnetic force F a generally has a more pronounced effect on the end of the winding, and the middle of the winding is more influenced by the radial electromagnetic force F r .

Phase Control Switch Principle
To analyze the electromechanical characteristics of the windings under different residual fluxes, the transformer is controlled by a compound switch composed of thyristors to turn on and off. Figure 19 shows a schematic diagram of the phase-controlled switch. The voltage signals of the circuit are obtained by the multi-channel data acquisition function of the hardware. After the Alternating Current (AC) excitation voltage crosses zero, the hardware controls the forward thyristor trigger angle to generate different phase signals, controls the reverse thyristor renewal, and turns on the AC contactor connecting circuit. Thus, the hardware realizes the communication between the master chip and the Personal Computer (PC) and sets the closing phase angle. Then the PC sends out the specified values to control the loading excitation to form different residual fluxes.

Windings Electromagnetic Characteristics Test
The test process (Figure 20) uses the recorder DL850 (Yokogawa, Tokyo, Japan) to record and store the final current and voltage signals, and uses the oscilloscope TEK DPO3054 (Tektronix, Beaverton, America) to display the voltage and current signals. Using the laser Doppler vibrometer LV-S01 (Sunny Optical Technology (Group) Co.,LTD, Zhejiang, China), the winding deformation displacement is measured under different residual fluxes of magnetic influence. The simultaneous recording of data using a recorder and an oscilloscope can mutually verify the correctness of experimental waveforms and experimental data. The leakage density between the winding air gaps is measured using a measuring coil and the leakage waveform is displayed by a wave recorder DL850. Figure 21 shows the peak inrush current under different residual fluxes.  The voltage signals of the circuit are obtained by the multi-channel data acquisition function of the hardware. After the Alternating Current (AC) excitation voltage crosses zero, the hardware controls the forward thyristor trigger angle to generate different phase signals, controls the reverse thyristor renewal, and turns on the AC contactor connecting circuit. Thus, the hardware realizes the communication between the master chip and the Personal Computer (PC) and sets the closing phase angle. Then the PC sends out the specified values to control the loading excitation to form different residual fluxes.

Windings Electromagnetic Characteristics Test
The test process (Figure 20) uses the recorder DL850 (Yokogawa, Tokyo, Japan) to record and store the final current and voltage signals, and uses the oscilloscope TEK DPO3054 (Tektronix, Beaverton, America) to display the voltage and current signals. Using the laser Doppler vibrometer LV-S01 (Sunny Optical Technology (Group) Co., LTD, Zhejiang, China), the winding deformation displacement is measured under different residual fluxes of magnetic influence. The simultaneous recording of data using a recorder and an oscilloscope can mutually verify the correctness of experimental waveforms and experimental data. The leakage density between the winding air gaps is measured using a measuring coil and the leakage waveform is displayed by a wave recorder DL850. Figure 21 shows the peak inrush current under different residual fluxes. The voltage signals of the circuit are obtained by the multi-channel data acquisition function of the hardware. After the Alternating Current (AC) excitation voltage crosses zero, the hardware controls the forward thyristor trigger angle to generate different phase signals, controls the reverse thyristor renewal, and turns on the AC contactor connecting circuit. Thus, the hardware realizes the communication between the master chip and the Personal Computer (PC) and sets the closing phase angle. Then the PC sends out the specified values to control the loading excitation to form different residual fluxes.

Windings Electromagnetic Characteristics Test
The test process (Figure 20) uses the recorder DL850 (Yokogawa, Tokyo, Japan) to record and store the final current and voltage signals, and uses the oscilloscope TEK DPO3054 (Tektronix, Beaverton, America) to display the voltage and current signals. Using the laser Doppler vibrometer LV-S01 (Sunny Optical Technology (Group) Co.,LTD, Zhejiang, China), the winding deformation displacement is measured under different residual fluxes of magnetic influence. The simultaneous recording of data using a recorder and an oscilloscope can mutually verify the correctness of experimental waveforms and experimental data. The leakage density between the winding air gaps is measured using a measuring coil and the leakage waveform is displayed by a wave recorder DL850. Figure 21 shows the peak inrush current under different residual fluxes.  From the above Figure 21, it can be seen that when the transformer is closed at noload, the inrush current generated is related to the residual flux. As the residual flux increases, the peak of the inrush current gradually increases. The main reason is that the larger the residual flux magnetic density Br in the core, the faster the core reaches saturation, the higher the degree of saturation, and the greater the value of the resulting inrush current. Figure 22 shows that when the direction of the residual flux magnetism changes, the direction of the leakage magnetism will also change, but the size of the residual flux magnetism has little effect on the leakage magnetism. The trend of air gap leakage at both ends of the winding is basically the same, and the difference is not significant.

Windings Mechanical Characteristics Test
To measure the mechanical characteristics of the transformer prototype winding under different residual flux magnetism, an acceleration sensor is used to measure the acceleration response, and a laser vibrometer is used to measure the winding displacement. By analyzing the acceleration frequency response during the winding vibration, important parameters such as the system intrinsic frequency, damping ratio and modal vibration pattern are obtained. Data acquisition of vibration signals is performed using the DH5902 (Donghua Test Technology Co., Jingjiang, Jiangsu, China) vibration test system. A general flow chart of the experiment is shown in Figure 23.  From the above Figure 21, it can be seen that when the transformer is closed at no-load, the inrush current generated is related to the residual flux. As the residual flux increases, the peak of the inrush current gradually increases. The main reason is that the larger the residual flux magnetic density B r in the core, the faster the core reaches saturation, the higher the degree of saturation, and the greater the value of the resulting inrush current. Figure 22 shows that when the direction of the residual flux magnetism changes, the direction of the leakage magnetism will also change, but the size of the residual flux magnetism has little effect on the leakage magnetism. The trend of air gap leakage at both ends of the winding is basically the same, and the difference is not significant. From the above Figure 21, it can be seen that when the transformer is c load, the inrush current generated is related to the residual flux. As the resi creases, the peak of the inrush current gradually increases. The main reaso larger the residual flux magnetic density Br in the core, the faster the core rea tion, the higher the degree of saturation, and the greater the value of the resu current. Figure 22 shows that when the direction of the residual flux magnetism direction of the leakage magnetism will also change, but the size of the residu netism has little effect on the leakage magnetism. The trend of air gap leakage of the winding is basically the same, and the difference is not significant.

Windings Mechanical Characteristics Test
To measure the mechanical characteristics of the transformer prototype der different residual flux magnetism, an acceleration sensor is used to measu eration response, and a laser vibrometer is used to measure the winding displ analyzing the acceleration frequency response during the winding vibratio parameters such as the system intrinsic frequency, damping ratio and mod pattern are obtained. Data acquisition of vibration signals is performed using (Donghua Test Technology Co., Jingjiang, Jiangsu, China) vibration test system flow chart of the experiment is shown in Figure 23.

Windings Mechanical Characteristics Test
To measure the mechanical characteristics of the transformer prototype winding under different residual flux magnetism, an acceleration sensor is used to measure the acceleration response, and a laser vibrometer is used to measure the winding displacement. By analyzing the acceleration frequency response during the winding vibration, important parameters such as the system intrinsic frequency, damping ratio and modal vibration pattern are obtained. Data acquisition of vibration signals is performed using the DH5902 (Donghua Test Technology Co., Jingjiang, Jiangsu, China) vibration test system. A general flow chart of the experiment is shown in Figure 23. The test focused on measuring the vibration information of the winding biscuit at the moment of transformer reclosing, based on the structural characteristics of the transformer winding biscuit and combined with the finite element simulation analysis results, the specific measurement acquisition points are shown in Figure 24 below. The arrangement of measurement points mainly follows the principle of avoiding the omission of the main location but also as convenient as possible, each transformer winding line cake on a total of seven points of measurement. The data of the seven vibration measurement points shown in Figure 24 above were collected at the residual flux Br = 0.349 and the closing angle θ = 0° and analyzed by spectrum. Figure 25 shows the peak winding acceleration for measurement points 1-7.   Figure 26 shows the peak winding acceleration at measurement points 1-7 with residual flux Br = 0.349 and closing angle θ = 0°. The test focused on measuring the vibration information of the winding biscuit at the moment of transformer reclosing, based on the structural characteristics of the transformer winding biscuit and combined with the finite element simulation analysis results, the specific measurement acquisition points are shown in Figure 24 below. The arrangement of measurement points mainly follows the principle of avoiding the omission of the main location but also as convenient as possible, each transformer winding line cake on a total of seven points of measurement. The test focused on measuring the vibration information of the winding biscuit at the moment of transformer reclosing, based on the structural characteristics of the transformer winding biscuit and combined with the finite element simulation analysis results, the specific measurement acquisition points are shown in Figure 24 below. The arrangement of measurement points mainly follows the principle of avoiding the omission of the main location but also as convenient as possible, each transformer winding line cake on a total of seven points of measurement. The data of the seven vibration measurement points shown in Figure 24 above were collected at the residual flux Br = 0.349 and the closing angle θ = 0° and analyzed by spectrum. Figure 25 shows the peak winding acceleration for measurement points 1-7.   Figure 26 shows the peak winding acceleration at measurement points 1-7 with residual flux Br = 0.349 and closing angle θ = 0°. The data of the seven vibration measurement points shown in Figure 24 above were collected at the residual flux B r = 0.349 and the closing angle θ = 0 • and analyzed by spectrum. Figure 25 shows the peak winding acceleration for measurement points 1-7. The test focused on measuring the vibration information of the winding biscuit at moment of transformer reclosing, based on the structural characteristics of the tra former winding biscuit and combined with the finite element simulation analysis resu the specific measurement acquisition points are shown in Figure 24 below. The arran ment of measurement points mainly follows the principle of avoiding the omission of main location but also as convenient as possible, each transformer winding line cake o total of seven points of measurement. The data of the seven vibration measurement points shown in Figure 24 above w collected at the residual flux Br = 0.349 and the closing angle θ = 0° and analyzed by sp trum. Figure 25 shows the peak winding acceleration for measurement points 1-7.   Figure 26 shows the peak winding acceleration at measurement points 1-7 with sidual flux Br = 0.349 and closing angle θ = 0°.  From Figure 26, it can be seen that the vibration characteristics ar ferent measurement points at different frequencies. Among them, me has the most serious vibration and the largest vibration acceleration a 97.7 Hz. Figure 27 shows the deformation of Line 1 cake and Line 5 cake u sidual fluxes.  Figure 27, it can be seen that the winding deformation gradu the increase of residual flux magnetic Br. The main reason is that whe magnetization Br in the core is larger, the peak of the inrush current wi air gap leakage will also increase, which in turn leads to an increase mation. In the middle of the winding, due to the increase of radial force, From Figure 26, it can be seen that the vibration characteristics are different for different measurement points at different frequencies. Among them, measurement point 3 has the most serious vibration and the largest vibration acceleration at the frequency of 97.7 Hz. Figure 27 shows the deformation of Line 1 cake and Line 5 cake under different residual fluxes. From Figure 26, it can be seen that the vibration characteristics are different for different measurement points at different frequencies. Among them, measurement point 3 has the most serious vibration and the largest vibration acceleration at the frequency of 97.7 Hz. Figure 27 shows the deformation of Line 1 cake and Line 5 cake under different residual fluxes. From Figure 27, it can be seen that the winding deformation gradually increases with the increase of residual flux magnetic Br. The main reason is that when the residual flux magnetization Br in the core is larger, the peak of the inrush current will increase and the air gap leakage will also increase, which in turn leads to an increase in winding deformation. In the middle of the winding, due to the increase of radial force, there will be more serious deformation, making the overall structure of the transformer an oval ball shape. Table 2 shows the comparison of the simulated value Imax1 and the experimental value Imax2 of the peak inrush current generated by the transformer under different residual fluxes Br when the closing phase angle θ = 0°. From Figure 27, it can be seen that the winding deformation gradually increases with the increase of residual flux magnetic B r . The main reason is that when the residual flux magnetization B r in the core is larger, the peak of the inrush current will increase and the air gap leakage will also increase, which in turn leads to an increase in winding deformation. In the middle of the winding, due to the increase of radial force, there will be more serious deformation, making the overall structure of the transformer an oval ball shape. Table 2 shows the comparison of the simulated value I max1 and the experimental value I max2 of the peak inrush current generated by the transformer under different residual fluxes B r when the closing phase angle θ = 0 • . From Table 2, it can be seen that the maximum error between the simulated value and the experimental value obtained by simulation is less than 0.3%, and the minimum error is less than 0.1%, demonstrating the effectiveness of the simulation analysis method. At the same time, with the increase of the residual flux B r , the amplitude of the inrush current gradually increases, compared with zero residual flux, the residual flux B r increases to 1.24 T when the inrush current amplitude I max increases by about 2.5 times. Table 3 shows the comparison between the simulated value ∆x 1 and the tested value ∆x 2 under different residual fluxes B r of line 5 cake.  Table 3, it can be seen that the maximum error between the simulated value ∆x 1 and the experimental value ∆x 2 obtained by simulation is less than 1%, and the minimum error is less than 0.1%, demonstrating the effectiveness of the simulation. Meanwhile, as the residual flux B r increases, the winding deformation gradually increases. Compared with zero residual flux, the winding deformation ∆x increases about 3 times when the residual flux B r increases to 1.24 T.

Discussion
Comparison and analysis of the experimental results and the simulation results proved that the different residual fluxes could lead to differences in the peak magnitude of the inrush current, the leakage magnetic field, the winding electromagnetic force, and the deformation of the winding. As the residual flux of the transformer core increases, the inrush current, as well as the winding deformation, increases significantly.

Conclusions
In this paper, the electromechanical characteristics of the no-load reclosing winding of a dry-type transformer under different residual fluxes are comprehensively analyzed. First, a simulation circuit was built to obtain the peak inrush current generated during no-load reclosing under different residual fluxes of the transformer core. At the same time, these current peaks are loaded as the excitation into the transformer model. Moreover, the electromechanical characteristics of the transformer winding under the influence of the inrush current generated during no-load reclosing are analyzed using the field-circuit coupling and the finite element method. Finally, for the verification tests, the closing operation of the transformer is controlled by the phase-controlled switch, and the residual flux at the closing moment is calculated based on the analysis of the breaking angle. The laser Doppler vibrometer and the measuring coil were used to measure the winding deformation and air gap leakage. On the basis of a comparative analysis of the experimental and simulation results, it can be seen that as the residual flux in the core increases, under the combined effect of the inrush current and the air gap leakage field, a large electromagnetic force is generated on the transformer winding. According to the simulation and experimental results, the axial electromagnetic force generally affects the winding end more obviously, and the middle position of the winding is more affected by the radial electromagnetic force.
The validity and reliability of the simulation method and the principle analysis are proved in this paper by their comparison. The comprehensive analysis and experimental verification of the electromechanical characteristics of the transformer windings in the case of no-load reclosing were performed in this paper, making it possible to effectively assess the performance of transformer winding parameters under the influence of inrush current. This research can be applied to the simulation analysis of different power transformers, and is a guide for extending the service life of transformers and the in-depth study of transformer fault analysis. With the increase of grid capacity, it is of great significance to the power system to ensure the safe and stable operation of transformers.
However, when transformers are subjected to multiple inrush current shocks in a relatively short period, their electromechanical characteristics are bound to be greatly affected. Therefore, in the next step, it is necessary to analyze the electromechanical characteristics of a transformer when it experiences multiple inrush currents.  Data Availability Statement: All data, models generated or used during the study appear in the submitted article.

Conflicts of Interest:
The authors declare no conflict of interest.