Modeling and Analysis of New Power Devices Based on Linear Phase-Shifting Transformer

: With the rapid development of new power systems, various new power devices have also been developed. It is very important to establish analytical models of new power devices to ensure or even improve the reliability and stability of the power system. A linear phase-shifting transformer (LPST) is a new type of power device that mainly relies on air gaps to transfer energy, so establishing an accurate air-gap magnetic ﬁeld model is very important for improving the efﬁciency of this system. In this paper, an analytical model of an unequal-pitch linear phase-shifting transformer (UP-LPST) was established by combining the distributed magnetic circuit method (DMCM) and Schwartz–Christopher transformation (SCT). Taking the magnetic ﬁeld strength as a variable, an accurate magnetic ﬁeld analysis model for a UP-LPST considering saturation, cogging, and edge was established. Taking a 1 kw UP-LPST as a prototype, the accuracy of the model was veriﬁed by the ﬁnite element method and experiments. This modeling method could also be used to establish magnetic ﬁeld models of other similar structures in new energy power systems, especially those with cogging structures.


Introduction
The operation of a new energy power system must have a high level of reliability.However, ensuring the stable operation of any line component of a power system is complex.Its stability and reliability are determined by different characteristic parameters of the system.The main influencing factors include different impedances of parallel lines in the system, the power factor, changes in the input power, and changes in the load [1][2][3].A variety of electrical devices have been developed, with the phase-shifting transformer being one of the most important technologies.This device can control the current distribution between the branches of a parallel power system through the adjustment of the phase angle.At the same time, it can solve the overload problem caused by the unbalanced impedance of parallel transmission lines and improve the stability and efficiency of the power system, so it has received extensive attention from scholars at home and abroad [4][5][6].
As a new type of phase-shifting transformer, the linear phase-shifting transformer (LPST) has many advantages.Compared with the traditional phase-shifting transformer, its core structure is simple, the air gap adjustment is easy, the phase-shifting angle is wide, and the volume is small [7][8][9].At the same time, it can effectively eliminate low-order harmonics, improve the quality of output waveforms, and reduce power grid harmonic pollution [10].The LPST is a power device that mainly relies on air gaps to transfer energy, and there are several rectangular slots on the primary and secondary sides of the LPST to hold winding coils.This slotted structure produces a degree of cogging, which results in a significant increase in the tooth harmonic amplitude and distorts the air-gap magnetic field.More seriously, it reduces the quality of the output waveform and increases the energy loss [11], thereby reducing the stability and efficiency of the power system.Furthermore,

Analytical Model of Slotless Magnetic Field
The structure of a UP-LPST is basically the same as that of a linear motor.A diagram of its main structure is shown in Figure 1.The difference is that the length of the core on the primary and secondary sides of the UP-LPST is the same, and its core is symmetric about the air gap.Four groups of three-phase windings are distributed longitudinally along the core on the primary side, and one group of three-phase windings on the secondary side.In contrast to conventional phase-shifting transformers, the energy conversion of a UP-LPST is mainly realized by an air-gap traveling-wave magnetic field.When the primary winding of the UP-LPST is energized, a linear traveling-wave magnetic field can be generated in the core, and then a three-phase electromotive force is induced in the secondary side.The UP-LPST studied in this paper comprised a pair of poles, as shown in Figure 1.The primary side was composed of four groups of three-phase bridge inverter circuits with a twelve-phase input (a1-x1, a2-x2, …, c3-z3, c4-z4), and the windings were distributed with full pitch.The secondary side comprised a three-phase output (A1-X1-A2-X2-A3-X3-A4-X4, …), and the windings were distributed in a combination of long-distance winding and short-distance winding.This winding distribution could effectively improve the three-phase asymmetry caused by the discontinuity of the core.The core structures of the primary and secondary sides were completely consistent and symmetric with respect to the air gap.
Firstly, in order to obtain the magnetic field without grooves, the following basic assumptions were made [26]: 1.The primary and secondary side end flux leakage is ignored.2. The flux lines in the virtual teeth are all arranged in the longitudinal direction, and the flux lines in the yokes are all arranged in the normal direction.3. The size of the primary-and secondary-side virtual teeth is the same.4. The interaction of adjacent slots (virtual slots) is ignored.
The overall model was divided into five regions along the longitudinal direction for magnetic circuit calculation, as shown in Figure 2, where I is the air-gap region; II and III are the primary-side virtual teeth and the yoke region, respectively; and IV and V are the secondary-side virtual teeth and the yoke region, respectively.For the loop shown in Figure 2: where Fδ_n (m) and Fδ_n (m + 1) are the air-gap normal magnetic pressure drops at node m and node m + 1, respectively; Fj1_l (m) and Fj2_l (m) are the longitudinal magnetic pressure drops of the yoke at the core node m on the primary and secondary sides, respectively; The UP-LPST studied in this paper comprised a pair of poles, as shown in Figure 1.The primary side was composed of four groups of three-phase bridge inverter circuits with a twelve-phase input (a1-x1, a2-x2, . . ., c3-z3, c4-z4), and the windings were distributed with full pitch.The secondary side comprised a three-phase output (A1-X1-A2-X2-A3-X3-A4-X4, . . .), and the windings were distributed in a combination of long-distance winding and short-distance winding.This winding distribution could effectively improve the three-phase asymmetry caused by the discontinuity of the core.The core structures of the primary and secondary sides were completely consistent and symmetric with respect to the air gap.
Firstly, in order to obtain the magnetic field without grooves, the following basic assumptions were made [26]: 1.
The primary and secondary side end flux leakage is ignored.

2.
The flux lines in the virtual teeth are all arranged in the longitudinal direction, and the flux lines in the yokes are all arranged in the normal direction.

3.
The size of the primary-and secondary-side virtual teeth is the same.4.
The interaction of adjacent slots (virtual slots) is ignored.
The overall model was divided into five regions along the longitudinal direction for magnetic circuit calculation, as shown in Figure 2, where I is the air-gap region; II and III are the primary-side virtual teeth and the yoke region, respectively; and IV and V are the secondary-side virtual teeth and the yoke region, respectively.The UP-LPST studied in this paper comprised a pair of poles, as shown in Figure 1.The primary side was composed of four groups of three-phase bridge inverter circuits with a twelve-phase input (a1-x1, a2-x2, …, c3-z3, c4-z4), and the windings were distributed with full pitch.The secondary side comprised a three-phase output (A1-X1-A2-X2-A3-X3-A4-X4, …), and the windings were distributed in a combination of long-distance winding and short-distance winding.This winding distribution could effectively improve the three-phase asymmetry caused by the discontinuity of the core.The core structures of the primary and secondary sides were completely consistent and symmetric with respect to the air gap.
Firstly, in order to obtain the magnetic field without grooves, the following basic assumptions were made [26]: 1.The primary and secondary side end flux leakage is ignored.2. The flux lines in the virtual teeth are all arranged in the longitudinal direction, and the flux lines in the yokes are all arranged in the normal direction.3. The size of the primary-and secondary-side virtual teeth is the same.4. The interaction of adjacent slots (virtual slots) is ignored.
The overall model was divided into five regions along the longitudinal direction for magnetic circuit calculation, as shown in Figure 2, where I is the air-gap region; II and III are the primary-side virtual teeth and the yoke region, respectively; and IV and V are the secondary-side virtual teeth and the yoke region, respectively.For the loop shown in Figure 2: where Fδ_n (m) and Fδ_n (m + 1) are the air-gap normal magnetic pressure drops at node m and node m + 1, respectively; Fj1_l (m) and Fj2_l (m) are the longitudinal magnetic pressure drops of the yoke at the core node m on the primary and secondary sides, respectively; For the loop shown in Figure 2: where F δ_n (m) and F δ_n (m + 1) are the air-gap normal magnetic pressure drops at node m and node m + 1, respectively; F j1_l (m) and F j2_l (m) are the longitudinal magnetic pressure drops of the yoke at the core node m on the primary and secondary sides, respectively; and F t1_n (m), F t1_n (m + 1), F t2_n (m), and F t2_n (m + 1) are the normal magnetic pressure drops of the virtual teeth at the core node m and node m + 1 on the primary and secondary sides, respectively.The positive direction of the coordinate axis is the positive direction of the magnetic pressure drop of each section.
Because the primary and secondary sides of the straight-line phase-shifting transformer were symmetric about the center line of the air gap, the model could be simplified to a one-sided model, as shown in Figure 3 [27].
Processes 2022, 10, x FOR PEER REVIEW 4 of 17 and Ft1_n (m), Ft1_n (m + 1), Ft2_n (m), and Ft2_n (m + 1) are the normal magnetic pressure drops of the virtual teeth at the core node m and node m + 1 on the primary and secondary sides, respectively.The positive direction of the coordinate axis is the positive direction of the magnetic pressure drop of each section.
Because the primary and secondary sides of the straight-line phase-shifting transformer were symmetric about the center line of the air gap, the model could be simplified to a one-sided model, as shown in Figure 3 [27].The total magnetic pressure drop of the circuit could be simplified as: Each MMF could be expressed as follows: where Hδ_n is the magnetic field intensity in the normal direction of the air gap; Bt1_n and Ht1_n denote the flux density and magnetic field intensity in the normal direction of the primary-side virtual teeth, respectively; Bj1_l and Hj1_l denote the flux density and magnetic field intensity in the longitudinal direction of the primary yoke, respectively; δ is the length of the air gap; ht1 is the virtual tooth height of the primary core; μ0 is the air permeability; and μj1 denotes the permeability of each node of the primary yoke.
When the tooth is not saturated, the main flux within a pitch can be considered to pass entirely through the tooth.At this point, the tooth magnetic flux densities at node m are as follows: where H and L are the height and length of the core of the UP-LPST, respectively; bt is the tooth width of the UP-LPST; t1 is the tooth pitch; and KFe is the superposition coefficient of the core.However, when the teeth are saturated, most of the main flux passes through the teeth, and the rest enters the yoke through the slot.At this time, the actual magnetic flux in the tooth becomes smaller, so it was necessary to revise (5) as follows: where B ꞌ t_n (m) is the normal apparent magnetic flux density of the tooth at node m, representing the magnetic flux density when all the flux enters the tooth; Bt_n (m) is the actual normal magnetic flux density of the tooth at node m; Ht_n (m) is the actual normal magnetic field intensity of the tooth at node m; and kδ is the slot coefficient, kδ = (H•bs)/(KFe•L•bt).The total magnetic pressure drop of the circuit could be simplified as: Each MMF could be expressed as follows: where H δ_n is the magnetic field intensity in the normal direction of the air gap; B t1_n and H t1_n denote the flux density and magnetic field intensity in the normal direction of the primary-side virtual teeth, respectively; B j1_l and H j1_l denote the flux density and magnetic field intensity in the longitudinal direction of the primary yoke, respectively; δ is the length of the air gap; h t1 is the virtual tooth height of the primary core; µ 0 is the air permeability; and µ j1 denotes the permeability of each node of the primary yoke.When the tooth is not saturated, the main flux within a pitch can be considered to pass entirely through the tooth.At this point, the tooth magnetic flux densities at node m are as follows: where H and L are the height and length of the core of the UP-LPST, respectively; b t is the tooth width of the UP-LPST; t 1 is the tooth pitch; and K Fe is the superposition coefficient of the core.However, when the teeth are saturated, most of the main flux passes through the teeth, and the rest enters the yoke through the slot.At this time, the actual magnetic flux in the tooth becomes smaller, so it was necessary to revise (5) as follows: where B t_n (m) is the normal apparent magnetic flux density of the tooth at node m, representing the magnetic flux density when all the flux enters the tooth; B t_n (m) is the actual normal magnetic flux density of the tooth at node m; H t_n (m) is the actual normal magnetic field intensity of the tooth at node m; and k δ is the slot coefficient, According to the continuity of magnetic flux, the longitudinal magnetic flux density of the yoke at node m is equal to the total normal magnetic flux density of the air gap from node 1 to node m.The magnetic flux density of the yoke was calculated as follows: where D is the thickness of the UP-LPST.Therefore, the magnetic flux density of the yoke of the primary core could be represented as follows: Meanwhile, it could be shown that the total air-gap flux in the whole length range of the UP-LPST is 0: For the UP-LPST, the harmonic frequency of the fundamental wave current is 24k ± 1, so the lowest is the 23rd harmonic.Due to the high harmonic order, only the fundamental MMF was considered to simplify the calculation.Thus, the MMF of each air-gap node was calculated as: where x(m) = (m − 1)dx (10) In this equation, N 1 is the number of turns on the primary side, I 1 is the primary measured current, k w1 is the winding coefficient, p is the polar logarithm, τ is the polar distance, and dx is the length of each block after segmentation.The normal magnetic flux density of the air gap at node m was calculated as follows: where K s is the preset saturation coefficient and µ 0 is the vacuum permeability.Finally, we judged whether the iteration precision value was satisfied: When the actual error was greater than the accuracy requirement, the air-gap flux density was corrected as follows [26]: where k s is the iterative coefficient.
The iteration flow chart is shown in Figure 4. Firstly, select the DC bus voltage and calculate the effective value of the primary current.Then, calculate the initial air-gap flux density B δ0_slotless_n (m) according to Equation (11).Calculate the magnetic density B j1_l (m) and B t1_n (m) of each node of the primary yoke and teeth according to Equations ( 6) and (8).Determine the magnetic permeability µ j1 (m) and µ t1 (m) of each node of the primary yoke and teeth using the B-H curve.Finally, carry out the iterative calculation until the judgment conditions are satisfied.
Processes 2022, 10, x FOR PEER REVIEW 6 of 17 The iteration flow chart is shown in Figure 4. Firstly, select the DC bus voltage and calculate the effective value of the primary current.Then, calculate the initial air-gap flux density Bδ0_slotless_n (m) according to Equation (11).Calculate the magnetic density Bj1_l (m) and Bt1_n (m) of each node of the primary yoke and teeth according to Equations ( 6) and (8).Determine the magnetic permeability μj1(m) and μt1(m) of each node of the primary yoke and teeth using the B-H curve.Finally, carry out the iterative calculation until the judgment conditions are satisfied.By extending the DMCM to the LPST, the magnetic field distribution of the air gap and iron core yoke could be obtained.This paper only describes the air gap in detail.The final result is shown in Figure 5: By extending the DMCM to the LPST, the magnetic field distribution of the air gap and iron core yoke could be obtained.This paper only describes the air gap in detail.The final result is shown in Figure 5:

Cogging Effect
According to [20], when the ratio of teeth width to air-gap length is greater than 2.44, the influence between adjacent slots can be ignored.However, the ratio of teeth width to air-gap length of the model used in this paper was much higher than 2.44.In order to facilitate the analysis, a single-slot model was adopted in this paper.Firstly, the following assumptions were made [20]: 1.The primary side is slotted, and the secondary side is a smooth plane.2. The magnetic conductivity of the core on the primary and secondary sides is infinite.3.Both the primary and secondary sides of the iron core are planes with equal magnetic potential, one of which is 0 and the other φ0.
When the slot depth of the UP-LPST had been determined, the polygon of the z plane could be obtained, as shown in Figure 6a.The relationship between the z plane and the w plane could be obtained from Table 1.below: where

Cogging Effect
According to [20], when the ratio of teeth width to air-gap length is greater than 2.44, the influence between adjacent slots can be ignored.However, the ratio of teeth width to air-gap length of the model used in this paper was much higher than 2.44.In order to facilitate the analysis, a single-slot model was adopted in this paper.Firstly, the following assumptions were made [20]: The primary side is slotted, and the secondary side is a smooth plane.

2.
The magnetic conductivity of the core on the primary and secondary sides is infinite.

3.
Both the primary and secondary sides of the iron core are planes with equal magnetic potential, one of which is 0 and the other ϕ 0 .
When the slot depth of the UP-LPST had been determined, the polygon of the z plane could be obtained, as shown in Figure 6a.

Cogging Effect
According to [20], when the ratio of teeth width to air-gap length is greater than 2.44, the influence between adjacent slots can be ignored.However, the ratio of teeth width to air-gap length of the model used in this paper was much higher than 2.44.In order to facilitate the analysis, a single-slot model was adopted in this paper.Firstly, the following assumptions were made [20]: 1.The primary side is slotted, and the secondary side is a smooth plane.2. The magnetic conductivity of the core on the primary and secondary sides is infinite.3.Both the primary and secondary sides of the iron core are planes with equal magnetic potential, one of which is 0 and the other φ0.
When the slot depth of the UP-LPST had been determined, the polygon of the z plane could be obtained, as shown in Figure 6a.The relationship between the z plane and the w plane could be obtained from Table 1.below: where The relationship between the z plane and the w plane could be obtained from Table 1.below: where * The interior angle of a polygon.
Introducing the Jacobian elliptic function, let w be the inverse function of the intermediate variable k, as follows: Then ( 14) can be expressed as: The transformation between the z-w plane is: where snα, cnα, and dnα are Jacobian elliptic functions; Π (k, α) is the elliptic integral of the third kind; and δ is the air-gap length.Through the corresponding relation of the z-w plane, the relation between a, α, b s , and δ could be obtained as follows: where K(1/a2) is the elliptic integral of the first kind; K (1/a2) the elliptic integral of the first kind of complementary modules; and Z(α) is the Jacobian zeta function.The expressions are: where F(x,1/a), E(x,1/a), and E(1/a 2 ) are the incomplete elliptic integral of the first kind, incomplete elliptic integral of the second kind, and complete elliptic integral of the second kind, respectively.
Processes 2022, 10, 1596 9 of 17 After logarithmic transformation, the corresponding relation of the t-w plane was written as: When the slot depth had been determined, the magnetic field density distribution of the air gap and slot in the UP-LPST was calculated as:

End Effect
Because the core of the UP-LPST is not continuous, the end effect as its inherent property also needed to be analyzed.The specific SCT analysis model is shown in Figure 7.
where F(x,1/a), E(x,1/a), and E(1/a2) are the incomplete elliptic integral of the first kind, incomplete elliptic integral of the second kind, and complete elliptic integral of the second kind, respectively.
After logarithmic transformation, the corresponding relation of the t-w plane was written as: When the slot depth had been determined, the magnetic field density distribution of the air gap and slot in the UP-LPST was calculated as:

End Effect
Because the core of the UP-LPST is not continuous, the end effect as its inherent property also needed to be analyzed.The specific SCT analysis model is shown in Figure 7.The relationship between the z plane and the w plane could be obtained from Table 2 below: The relationship between the z plane and the w plane could be obtained from Table 2 below: Similarly, the relationship of side z-w-t is: After considering the end effect, the magnetic flux density on the primary core was determined to be: The magnetic flux density reached its maximum at point z 2 in the z plane: With B max as the base value, the magnetic field density at any point is as follows: By inserting certain w values, the magnetic field density distribution curve could be drawn as shown in Figure 8.
After considering the end effect, the magnetic flux density on the primary c determined to be: The magnetic flux density reached its maximum at point z2 in the z plane:   The relative permeability distribution function of the primary side considering the end effect could be obtained by fitting the curve in Figure 8.
The relative permeability of the air gap obtained when the secondary side was slotted separately was similar to that of the primary side, because the core structures of the primary and secondary sides of the UP-LPST are exactly the same.According to the above analysis, the air-gap magnetic field density B(x) of the UP-LPST along the length direction could be calculated as follows:

Analytical Results
In this paper, an unequal-pitch linear phase-shifting transformer was taken as an example.The specific parameters are shown in Table 3.Based on the equations presented in Section 2, the air-gap relative permeance at different normal positions within a range of tooth pitch could be obtained.In order to analyze the variation trend of the air-gap relative permeance λ in a single slot, four different positions were selected: the center line of the air gap (y = 0 mm), the outer surface of the primary side (y = δ/2 mm), the center of the primary-side slot (y = δ/2 + h/3 mm), and the bottom of the primary-side slot (y = δ/2 + hmm).Figure 9 shows that the air-gap permeability distribution was different for different air-gap radii.The closer to the opening surface of the slot, the greater the influence of the opening slot on the air-gap magnetic field, that is, the deeper the pit.The air-gap relative permeance in a single slot was decomposed by Fourier transformation with the tooth pitch as the period.In this way, the air-gap relative permeance across the whole length range of the UP-LPST could be obtained, as shown in Figure 10a.Figure 10b shows that the magnetic field density distribution at the center of the air gap considering the effects of cogging and end could be obtained by multiplying the calcu- The air-gap relative permeance in a single slot was decomposed by Fourier transformation with the tooth pitch as the period.In this way, the air-gap relative permeance across the whole length range of the UP-LPST could be obtained, as shown in Figure 10a.Figure 10b shows that the magnetic field density distribution at the center of the air gap considering the effects of cogging and end could be obtained by multiplying the calculated air gap relative permeability by the slotless magnetic field.The air-gap magnetic field was not highly sinusoidal when considering the effect of slot and end and was not affected by slot openings that were opposite teeth.However, when the slots were opposite each other, the air-gap magnetic field was affected by the interaction of the iron cores on both sides, presenting a concave shape.At the same time, due to the influence of "out and in", the magnetic flux density distribution at both ends was different.The air-gap relative permeance in a single slot was decomposed by Fourier transformation with the tooth pitch as the period.In this way, the air-gap relative permeance across the whole length range of the UP-LPST could be obtained, as shown in Figure 10a.Figure 10b shows that the magnetic field density distribution at the center of the air gap considering the effects of cogging and end could be obtained by multiplying the calculated air gap relative permeability by the slotless magnetic field.The air-gap magnetic field was not highly sinusoidal when considering the effect of slot and end and was not affected by slot openings that were opposite teeth.However, when the slots were opposite each other, the air-gap magnetic field was affected by the interaction of the iron cores on both sides, presenting a concave shape.At the same time, due to the influence of "out and in", the magnetic flux density distribution at both ends was different.Figure 11 presents the schematic diagram of harmonic amplitude and flux density distribution under different slot spacing ratios.The main harmonic amplitudes could be obtained by the Fourier decomposition of the air-gap magnetic fields with different slot spacing ratios.When the slot opening increased, the harmonic amplitude of the basic teeth increased significantly, which reduced the performance of the straight-line phase-shifting transformer.In Figure 11b, the lengths of the slot openings are 10 mm, 11 mm, 12 mm, and 13 mm.The minimum flux density was located in the middle of the slot opening, while the maximum flux density was located in the tooth, and the large slot opening had a greater impact on the flux density.

FEM Verification
To verify the correctness of the analytical results, an FEM model of the UP-LPST was constructed according to Table 1, as shown in Figure 12.A 100V DC voltage was applied to the FEM model, and the magnetic field distribution data at the center line of the air gap were extracted and compared with the analytical results.The comparison results are shown in Figure 13.The results showed that the overall distribution fitted well, but the details were slightly different.The errors were mainly focused on the top and edge of the sawtooth wave.This was because the magnetic field lines at the teeth and air gap were assumed to have only normal components.However, near the opening of the slot, the effect of magnetic focusing on the air-gap magnetic field was great, which was one of the reasons for the error.At the same time, this method did not consider the magnetic flux leakage effect of the UP-LPST, so there was a 2% error between the analytical results and the FEM results.
obtained by the Fourier decomposition of the air-gap magnetic fields with different slot spacing ratios.When the slot opening increased, the harmonic amplitude of the basic teeth increased significantly, which reduced the performance of the straight-line phase-shifting transformer.In Figure 11b, the lengths of the slot openings are 10 mm, 11 mm, 12 mm, and 13 mm.The minimum flux density was located in the middle of the slot opening, while the maximum flux density was located in the tooth, and the large slot opening had a greater impact on the flux density.

FEM Verification
To verify the correctness of the analytical results, an FEM model of the UP-LPST was constructed according to Table 1, as shown in Figure 12.A 100V DC voltage was applied to the FEM model, and the magnetic field distribution data at the center line of the air gap were extracted and compared with the analytical results.The comparison results are shown in Figure 13.The results showed that the overall distribution fitted well, but the details were slightly different.The errors were mainly focused on the top and edge of the sawtooth wave.This was because the magnetic field lines at the teeth and air gap were assumed to have only normal components.However, near the opening of the slot, the effect of magnetic focusing on the air-gap magnetic field was great, which was one of the reasons for the error.At the same time, this method did not consider the magnetic flux leakage effect of the UP-LPST, so there was a 2% error between the analytical results and the FEM results.Compared with the exact subdomain method (ESM) in [28], the error of the DMCM-SCT method was smaller due to the consideration of the influence of core saturation and edge, and the specific performance results are shown in Table 4.   Compared with the exact subdomain method (ESM) in [28], the error of the DMCM SCT method was smaller due to the consideration of the influence of core saturation and edge, and the specific performance results are shown in Table 4. Compared with the exact subdomain method (ESM) in [28], the error of the DMCM-SCT method was smaller due to the consideration of the influence of core saturation and edge, and the specific performance results are shown in Table 4.According to Table 4, compared with ESM, the results obtained by DMCM-SCT were closer to those obtained by the finite element method.Although a certain calculation dimension was added, the overall effect was still better than ESM.However, the influence of end magnetic flux leakage was not considered in the calculation method, so there were still some errors in the results.

Test Verification
According to the data in Table 1, a prototype of the UP-LPST was manufactured, as shown in Figure 14.The test platform was mainly composed of a power supply driver module, a DSP, a signal amplification control circuit, a rectifier bridge, an inverter module, a terminal, and a UP-LPST.Since the air gap of the experimental prototype was only 0.3 mm, it was difficult to directly measure the actual air-gap flux density.Therefore, we carried out indirect verification by testing the no-load output voltage and current of the prototype.

Test Verification
According to the data in Table 1, a prototype of the UP-LPST was manufactured, as shown in Figure 14.The test platform was mainly composed of a power supply driver module, a DSP, a signal amplification control circuit, a rectifier bridge, an inverter module, a terminal, and a UP-LPST.Since the air gap of the experimental prototype was only 0.3 mm, it was difficult to directly measure the actual air-gap flux density.Therefore, we carried out indirect verification by testing the no-load output voltage and current of the prototype.Figure 15 shows the experimental waveform of the no-load output voltage of the prototype and the waveform of the FEM simulation calculation.By comparison, it could be found that the waveforms of the no-load output voltage calculated by FEM were basically the same as those measured by the experiment, but there were still slight differences.The main reason was that the manufacturing technology of the prototype was not ideal and the precision of core was not sufficient.At the same time, in the actual processing experiment, it was difficult to ensure that the length of the air gap was 0.3 mm, so the length of the air gap was not completely equal.However, the overall waveform direction was consistent, which indirectly verified the effectiveness of the analysis method in this paper.Figure 15 shows the experimental waveform of the no-load output voltage of the prototype and the waveform of the FEM simulation calculation.By comparison, it could be found that the waveforms of the no-load output voltage calculated by FEM were basically the same as those measured by the experiment, but there were still slight differences.The main reason was that the manufacturing technology of the prototype was not ideal and the precision of core was not sufficient.At the same time, in the actual processing experiment, it was difficult to ensure that the length of the air gap was 0.3 mm, so the length of the air gap was not completely equal.However, the overall waveform direction was consistent, which indirectly verified the effectiveness of the analysis method in this paper.
found that the waveforms of the no-load output voltage calculated by FEM were basically the same as those measured by the experiment, but there were still slight differences.The main reason was that the manufacturing technology of the prototype was not ideal and the precision of core was not sufficient.At the same time, in the actual processing experiment, it was difficult to ensure that the length of the air gap was 0.3 mm, so the length of the air gap was not completely equal.However, the overall waveform direction was consistent, which indirectly verified the effectiveness of the analysis method in this paper.

Conclusions
The impact of connecting new energy sources to a power system on the power quality of the grid is mainly reflected in the voltage and current.The new power system will consume more reactive power during operation, which will cause a serious drop in the grid voltage.Therefore, a large number of rectifier and inverter devices need to be used in the grid connection.The existence of these devices will inevitably inject harmonic currents into the grid, resulting in the distortion of grid voltage and current waveforms.While affecting the power supply quality, it will also cause additional load losses to the power equipment flowing through the distorted current.The special structure and phase-shifting method of the linear phase-shifting transformer can effectively eliminate low-order harmonics and improve the output waveform quality.For the LPST, which mainly relies on air gaps for energy transfer, the establishment of an accurate magnetic field model has an important influence on the calculation and even the improvement of the system efficiency.
By focusing on the specifics of LPST energy transfer and analyzing its structure, the DMCM and SCT were extended to the LPST, and a magnetic field analysis model for the LPST was established.The model can simultaneously consider the effects of saturation, cogging, and edges on the magnetic field.Therefore, the model is more in line with the actual situation of the linear phase-shifting transformer.Under the conditions of a given input DC bus voltage, the accuracy of the model was proven by means of direct verification using a finite element model and indirect verification via an experiment.There was a degree of error in the results, because the magnetic leakage at the end of the linear phase-shifting transformer was not considered.However, the model can still be applied to phase-shifting transformers with similar structures, and even other power devices.Our model could play an important role in improving the efficiency of transformers and new power systems and reducing losses.However, this paper only analyzed the model under no-load conditions.An analysis under load conditions is the next step.

Figure 2 .
Figure 2. Schematic diagram of the block model of the UP-LPST.

Figure 2 .
Figure 2. Schematic diagram of the block model of the UP-LPST.

Figure 2 .
Figure 2. Schematic diagram of the block model of the UP-LPST.

Figure 3 .
Figure 3. Diagram of the simplified model.

Figure 3 .
Figure 3. Diagram of the simplified model.

Figure 4 .
Figure 4.The iterative flowchart of the UP-LPST based on the DMCM.

Figure 4 .
Figure 4.The iterative flowchart of the UP-LPST based on the DMCM.

Processes 2022 , 17 Figure 5 .
Figure 5.The slotless air-gap magnetic field within the length range of the UP-LPST.

Figure 5 .
Figure 5.The slotless air-gap magnetic field within the length range of the UP-LPST.

Processes 2022 , 17 Figure 5 .
Figure 5.The slotless air-gap magnetic field within the length range of the UP-LPST.

With
Bmax as the base value, the magnetic field density at any point is as follo By inserting certain w values, the magnetic field density distribution curve c drawn as shown in Figure8.

Figure 8 .
Figure 8. Distribution of air-gap magnetic field density at the UP-LPST end.

Figure 8 .
Figure 8. Distribution of air-gap magnetic field density at the UP-LPST end.

Figure 9 .
Figure 9. Air-gap relative permeance at different positions in a single slot.

Figure 9 .
Figure 9. Air-gap relative permeance at different positions in a single slot.

Figure 9 .
Figure 9. Air-gap relative permeance at different positions in a single slot.

Figure 10 .
Figure 10.Distribution of air-gap relative permeance and flux density within the length range of the UP-LPST.(a) Air-gap relative permeance; (b) distribution of air-gap flux density under the influence of cogging and end.

Figure 10 .
Figure 10.Distribution of air-gap relative permeance and flux density within the length range of the UP-LPST.(a) Air-gap relative permeance; (b) distribution of air-gap flux density under the influence of cogging and end.

Figure 11 .
Figure 11.Distribution of harmonic amplitude and magnetic field density at different slot spacing ratios.(a) Harmonic amplitude value; (b) single-slot flux density distribution.

Figure 11 . 17 Figure 12 .
Figure 11.Distribution of harmonic amplitude and magnetic field density at different slot spacing ratios.(a) Harmonic amplitude value; (b) single-slot flux density distribution.Processes 2022, 10, x FOR PEER REVIEW 14 of 17

Figure 13 .
Figure 13.The air-gap flux density distribution of the UP-LPST.

Figure 13 .
Figure 13.The air-gap flux density distribution of the UP-LPST.

Figure 13 .
Figure 13.The air-gap flux density distribution of the UP-LPST.

Figure 15 .
Figure 15.No-load phase-A output voltage diagram of finite element model and experiment.(a) FEM waveform; (b) experimental waveform.Figure 15.No-load phase-A output voltage diagram of finite element model and experiment.(a) FEM waveform; (b) experimental waveform.

Figure 15 .
Figure 15.No-load phase-A output voltage diagram of finite element model and experiment.(a) FEM waveform; (b) experimental waveform.Figure 15.No-load phase-A output voltage diagram of finite element model and experiment.(a) FEM waveform; (b) experimental waveform.

Table 1 .
SCT table of z plane and w plane of finite slot depth model.

Table 2 .
SCT table of z plane and the w plane at the primary-side end.

Table 3 .
Parameters of the UP-LPST.

Table 4 .
Comparison of the three methods.

Table 4 .
Comparison of the three methods.