A Novel Analytic Method to Calculate the Equivalent Stray Capacitance of the Low-Speed Maglev Train’s Suspension Electromagnet

: The suspension electromagnet (SEM) is one of the most critical components of the low-speed maglev train to achieve the stable suspension and non-contact operation with the track. Therefore, it is valuable for optimizing the maglev system to study the characteristics of the SEM. This paper proposes a novel analytic method to calculate the equivalent stray capacitance of the SEM efficiently. It considers the capacitances inside the winding (turn to turn and layer to layer) and between the winding and core. Firstly, utilizing the compensation analytic method (CAM) calculates the static capacitances to reduce the calculation error, and the results of the CAM have significant improvement comparing with that of the traditional analytic method (TAM). Secondly, the analytical formula of the core floating potential is derived based on the partial capacitance theory, which has a 2.1% relative error to the finite element method (FEM). Finally, the experiment results of the SEM of a test rig prove the accuracy and efficiency of the proposed method.


Introduction
Suspension electromagnet is one key component of the low-speed maglev train, which makes the maglev train stable, suspending through conversing electrical energy to magnetic energy [1,2]. Essentially, the SEM is a time-varying inductor, whose impedance would become capacitive at the natural frequency, owing to the existence of the stray capacitance [3]. The large stray capacitance would lead to the impulse current for the fast-switching occasion [4][5][6] and contribute to the electromagnetic interference (EMI) noise [7,8]. In Shanghai Lin'gang Maglev Test Line, a capacitancefree differential mode inductance is installed between the maglev chopper and the suspension electromagnet winding to reduce the inrush current of the chopper switches and suppress the EMI interference of the suspension controller. So far, a few pieces of literature have investigated the utilization of stray capacitance in resonant converters [9,10] and EMI filter [11,12]. On the contrary, a significant amount of research has been devoted to reducing the stray capacitance by optimizing the winding's structure; preserving its inductive impedance in a high-frequency range also prevents the EMI noise [4,13,14]. Playing merit or avoiding malpractice, in any case, we should understand its formation mechanism and parameter extracting method firstly.
At present, the methods of extracting the stray capacitance mentioned in the literature are the measurement method [6,15,16], the finite element method (FEM) [17,18], and the analytic method [19][20][21]. The measurement methods are black-box techniques, which can help to obtain the equivalent stray capacitance of the suspension electromagnet but have no help to optimize the stray capacitance due to vague physical meaning. The FEM can extract the stray capacitance and is more accurate compared to the analytic methods, but it is a time-consuming method, especially for the threedimensional FEM model. Therefore, an analytical method, which is time-saving and helpful to design and optimize the winding structure, is requisite. L. F. F. Gutierres et al. derived an analytic method for calculating stray capacitances of an air-core inductor [19], which is not suitable for the calculation of stray capacitance of the suspension electromagnet due to the existence of an iron core. Majid Ezati Mosleh et al. proposed an analytic method to calculate a slide-conductive rotor magneto cumulative generator's stray capacitance [20], which neglects the effects of the core potential on the electric field energy storage. Zhan Shen et al. [21] systematically took the influence of the iron core potential on the electric field energy into account for the first time, furthermore given several approximate methods to estimate the iron core floating potential. At the same time, Zhan Shen's team gave the equivalent circuit method [21], which is suitable for the case of including a few capacitors, whose calculation error increases with the increase in the number of the winding layers. This paper derives an analytical formula to calculate the equivalent stray capacitance of the SEM, considering the effect of the iron core on electric field energy storage. Besides, the core potential is derived utilizing the partial capacitance theory, which makes the proposed method applicable to the situation of core floating. The stray capacitance can effectively reduce by increasing the thickness of the insulation between the winding and core. Finally, the FEM numerical simulation and experiment results verify the derivation.
The rest of this paper is organized as follows. Section 2 applies the TAM and CAM, respectively, to calculate the static capacitances. Section 3 derives the analytical expressions of the iron core floating potential and the equivalent stray capacitance. In Section 4, we have utilized the FEM in different SEM models to prove the correctness and efficiency of the proposed method. Section 5 delivers the experiment results of a test rig to verify the derivation and discussion. Finally, several conclusions are drawn in Section 6. Figure 1 presents the suspension electromagnet's structure (without the track), which mainly includes a winding and an iron core. Regarding the SEM as a multi-conductor system (MCS), the stray capacitance of SEM is the result of the potential difference between any two units (iron core and per-turn winding) in the MCS. The SEM winding of the low-speed maglev train is manufactured with the rectangular section wire, and the standard winding method [9,10] is adopted. The capacitance is defined, which is determined by winding shape, size, space position, and insulation parameters, as the static capacitance (inherent capacitance) [9,21]. Figure 2 shows the winding arrangement method and the composition of static capacitances. Where the x-axis represents the marker of the turns of the winding, the y-axis represents the marker of the layers of the winding. Nt is the number of turns per-layer, Nc is the number of the layers of the winding, CtL is the static capacitance of two adjacent turns of the same layer, CtLL is the static capacitance of two adjacent turns of different layers, Ctn is the static capacitance of non-adjacent turns. Ccd is the static capacitance between 2th~(Nt-1)th turns of the first layer and the iron core, Ccu is the static capacitance between 1th/Ntth of the 2th~Ncth layers and the iron core, and Ccdu is the static capacitance between the 1th/Ntth turn of the first layer and the iron core. The equivalent stray capacitance of the winding is calculated based on electric field energy equivalence; that is, the electrical field energy stored in the equivalent capacitance is equal to that of stored in the electromagnet MCS. Therefore, the static capacitances and voltage distribution are the basis for calculating the equivalent capacitance.

Static Capacitance Calculation by TAM
The TAM is based on the parallel-plate capacitance model. Assuming that the winding is in order and the faceplate is smooth, it is appropriate to regard all static capacitances in the SEM as parallel-plate capacitors. The most excellent parallel-plate capacitor model is [9,10,21] where ε0 and εeff are the vacuum permittivity and effective relative permittivity, respectively. Leff and weff are the length and width of the polar plates, respectively. deef is the effective distance between two plates. For the general case of three-layer insulation in series, the relative permittivity of insulation is ( ) ε ε ε δ δ δ ε ε δ ε ε δ ε ε δ ε where ε1, ε2, and ε3 (δ1, δ2, and δ3) are the relative permittivity (thickness) of the three series connection insulating layers.  where a1~a9 are the markers of the bare wires, wc and hc are the width and thickness of bare wire, respectively, wa is the thickness of iron core insulation, wt is the thickness of conductor insulation, wL is the thickness of filled insulation between adjacent turns of the same layer, wLL is the thickness of insulation between adjacent turns of different layers, and wcc is the thickness of filled insulation between the winding and the core. R1 and R2 are the radius of the arc part of the innermost winding and the outermost winding, respectively. Lc is the length of the linear part of the winding. Let the relative permittivity of the filled insulation (F-insulation), the wire insulation (W-insulation), and the iron core insulation (C-insulation) be ε f , ε w , and ε c , respectively. Combined with Figure 3 and  π ε ε ε ε ε Afterward, based on the geometrical relationship between the bare wire a9 and the iron core in Figure 3a, Ccdu is equivalent to that of Ccd parallel and connects with a part of Ccu. Therefore, Ccdu is calculated by Equation (7).

Static Capacitance Calculation by CAM
Compared with the TAM, the CAM compensates for the edge effect by constructing the virtual capacitor. Part 2.1 derives the TAM to calculate the static capacitances of the SEM, which ignores the influence of the edge effect. However, there are four types of special zones in the SEM shown in Figure 4. The A-type zone shown in Figure 4a is surrounded by four adjacent turns winding. The Btype zone shown in Figure 4b is surrounded by two adjacent turns winding of the same layer and the iron core's side plate. The C-type zone shown in Figure 4c is surrounded by the iron core's side and bottom plates and one turn winding. The D-type zone shown in Figure 4d is surrounded by two adjacent turns winding of the first layer and the iron core's bottom plate. This section considers the influence of the electric field energy stored in the zones above-mentioned to improve the calculation accuracy of the static capacitances. The specific technique is to construct the virtual capacitors in these areas to compensate for the edge effect. The A-type zone in Figure 4a is divided into areas A1~A4. In A1, a virtual parallel plate capacitor is constructed, Cv-tL (the subscript v represents the virtual capacitance, the same below), whose two pole plates are a1a6 and a2a7. The electric field line's route in Cv-tL is parallel to (a1-o-a2). Therefore, the Submitting Equation (8) into Equation (1), Cv-tL can be obtained. By symmetry, the virtual capacitance of A3 is the same as that of A1. Besides, A1 accounts for half of the virtual capacitor space, and the compensated capacitance of CtL should include both A1 and A3 at the same time. Therefore, the static capacitance Cr-tL (the subscript r represents the capacitance after recompensing the edge effect, the same below) is calculated by A virtual parallel plate capacitor (Cr-tLL), corresponding to A2, has the plates a2a8 and a3a9, whose electric field line's route is coordinated to (a2-o-a3). So, the parameters to calculate Cr-tLL are ( ) ( ) Submitting Equation (10) into Equation (1), Cv-tLL is obtained. By symmetry, the virtual capacitance of A4 is the same as that of A2. In addition, A2 accounts for half of the virtual capacitor space, and the compensated capacitance of CtLL includes both A2 and A4 at the same time. So, Cr-tLL is calculated by The B-type zone in Figure 4b is divided into areas B1~B3. The compensation capacitance of B1 is the same as that of A2. In A3, a virtual parallel plate capacitor is constructed, Cv-cu, whose two plates are b1b6 and b4b5, respectively. The electric field line's route in Cv-cu is parallel to (b1-o-b5). Therefore, the parameters to calculate Cv-cu are ( ) ( ) Submitting Equation (12) into Equation (1), Cv-cu can be obtained. According to symmetry, the virtual capacitance of B2 is the same as that of B3. In addition, A2 accounts for three-quarters of the virtual capacitor space, and the compensated capacitance of Ccu includes both B2 and B3 at the same time. So, Cr-cu is calculated by ( ) The C-type zone in Figure 4c is divided into areas C1~C3. The compensation capacitance corresponding to C1 is the same as that of A1. In C2, a virtual capacitor Cv-cd is calculated, whose two plates are c2c6 and c5c3, respectively. The electric field line's route is parallel to (c2-o-c5). Therefore, the parameters to calculate Cv-cd are ( ) Submitting Equation (14) into Equation (1), Cv-cd can be obtained. According to symmetry, the virtual capacitance of C3 is the same as that of C3. In addition, C2 accounts for three-quarters of the virtual capacitor space, and the compensated capacitance of cd C includes both C2 and C3 at the same time. So, Cr-cd is calculated by ( ) The D-type area in Figure 4d is divided into areas D1 and D2. In D1, a virtual parallel plate capacitor is constructed, Cv-cdu, whose two plates are d1d5 and d3d2. The electric field line's route in Cvcdu is coordinated to (d1-d3). So, the parameters to calculate Cv-cdu are Submitting Equation (16) into Equation (1), Cv-cdu can be obtained. Because D1 accounts for half of the virtual capacitor space, the ratio of the effective length of D2 to that of D1 is In addition, the compensated capacitance of Ccd includes D1、D2 、B2, and C3 at the same time. So, Cr-cdu is calculated by Through the above derivation, the static capacitances of the SEM with recompensing the edge effect are obtained by Equation (1) and Equations (8) ~ (17).

Calculation of Electric Field Energy and Equivalent Stray Capacitance
Section 2 of this paper calculates the static capacitances at different positions of the SEM by TAM and CAM. This section calculates the electric field energy stored in the SEM, then derives the equivalent stray capacitance based on the principle of energy storage equivalent [21][22][23][24].

Electric Field Energy Inside the Winding
The electric field energy inside the winding is mainly stored in the capacitances CtL and CtLL, and the potential and potential difference distribution of the winding influence the amount of accumulation energy. The effective capacitance changes according to the charging state of the static capacitance, but the charge distribution depends on the potential distribution. For the Nc×Nt turns suspension electromagnet wingding, assuming that each turn of the winding bears the same voltage, Figure 5 gives the potential and the potential difference distribution of the winding. In Figure 5a, the x-axis is the length of the winding, Lt is the length of per turn winding, u-axis is the potential and potential difference, ut1 and ut2 are the potential distribution of the first turn winding and the second turn winding, respectively, utd is the potential difference between ut1 and ut2, Ut is the voltage borne by per turn winding. In Figure 5b, the x-axis is the number of turns winding, the u-axis is the same as Figure 5a, Nt is the turns number of per layer winding, uc1 and uc2 are the potential distribution of the first layer winding and the second layer winding, respectively, ucd is the potential difference between uc1 and uc2, and Uc is the voltage borne by per layer winding.
According to Figure 5a,b, the electric field energy stored in CtL and CtLL is calculated by Equations (18) and (19), respectively.
where tL,all W and tLL,all W are the electric field energy stored in all tL C and all tLL C , respectively.

Electric Field Energy between Winding and Iron Core
The iron core of the SEM of the low-speed maglev train is a solid-core conductor, which is insulated with the winding and surrounded by an electrostatic field produced by windings. So, the core is a floating potential conductor. With the hypothesis that the SEM winding's arrangement is according to the structure shown in Figure 2, the lowest potential of the winding is U0, the potential difference between the first and the last turn windings is U, the iron core potential is Ue. In case the potential distribution is from low to high, the electric field energy stored in Ccdu, Ccd, and Ccu is obtained by Equation (20).
According to the above analysis, the electrostatic induction floating potential of the iron core should be solved firstly for calculating the electric field energy stored between the winding and the iron core. Many achievements have been achieved in the numerical solution and simulation of floating potential [25][26][27][28]. But a few pieces of literature have mentioned the analytical calculation method of floating potential. Figure 6 illustrates the relationship between the winding and the iron core potential distribution. The iron core floating potential is between U1 and UN (Range 1 in Figure  6). Because the distance between the inner side winding and the iron core is closer than that of the outer side winding, Ue is limited between U1 and ( ) Figure 6). Consequently, an assumption is ( ) Figure 6), which is suitable for approximate calculation. Furthermore, utilizing the equivalent circuit method based on KCL, Zhan Shen et al. derived the analytical formula to calculate the iron core floating potential. The iron core floating potential calculation error of the equivalent circuit method is expanding with the increase in winding layers number [21]. For example, for an 8-layers winding, the calculation error reaches 15%. To improve the accuracy of the floating potential calculation result, the partial capacitance theory is applied to deduce the analytical expression, which is more flexible than the equivalent circuit method.
For an N-conductor system, the − Q V equation is written as Equation (22).
where 0 j C is the capacitance between the j conductor and ground. Suppose that the Nth conductor in the conductors' system is floating potential, then For the SEM, the iron core floating potential is calculated by Equation (24), with the known winding potential. The capacitances between the iron core and the adjacent conductor are much larger than the capacitances between the iron core and ground and the non-adjacent conductor. So, it is reasonable to ignore the influence of the latter to simplify the calculation. When the potential distribution of the winding is from low to high, according to Equation (24), the iron core floating potential is calculated by Equation (25). (25) where ′ k is the voltage division factor, which is calculated by Equation (26).
When the potential distribution of the winding is from high to low, the iron core floating potential is calculated by Equation (27).

Equivalent Stray Capacitance
According to the principle of energy storage equivalent, Equation (31) gives the formula to calculate the equivalent stray capacitance.
,all tL,tLL,cdu,cd,cd Where Wall is the total electric field energy stored in the SEM. Submitting Equations (18) where tL k , tLL k , cdu k , cd k , and cu k are the coefficients, which are determined by the inherent structure of the winding. Replacing tL C , tLL C , cdu C , cd C , and cu C of Equations (18) ~ (29) with r-tL C , r-tLL C , r-cdu C , r-cd C , and r-cu C , respectively, the parameters compensating the edge effect are obtained.

Numerical Simulation and Analysis by FEM
The FEM is an effective and accurate electromagnetic field analysis method, which is widely recognized by academics and engineering, and it can extract the capacitances of the suspension electromagnet. Two alternatives-two-dimensional (2D) and three-dimensional (3D) model of the researched object can verify the correctness of the theory and analysis. The 2D model has the advantages of simple geometric modeling, few meshes, and low computational cost; the 3D model can obtain more accurate results but with the price of a large amount of calculation and timeconsuming, especially for the large size and complex model. Considering that the winding turns' number is large, and the bare conductor and insulation dimensions are very small relative to that of the iron core, the 2D finite element model is used to study the stray capacitance of the SEM. The following works are delivered to verify the correctness of the proposed analytic method: firstly, the FEM, TAM, and CAM are applied to solve the same electromagnet model, then the results obtained by three methods are compared; secondly, taking the thickness of the insulation between the iron core and winding (wcc) as the independent variable, the applicability of TAM and CAM to different parameter conditions is studied.

Results Comparison of the Different Methods
The general steps of extracting the stray capacitance between two conductors by the Maxwell software are modeling geometric, assigning excitations for the conductors, solving the scalar potential, calculating the electric field strength, and computing the stray capacitance. Figure 7 shows the geometric model of the SEM. where Wcu, Hcd, Wwin, and Hwin represent the width of the iron core, the thickness of the core leg, the width of the window, and the height of the window, respectively. The constraint relationship shown in Equation (30) should be satisfied among the dimensions' parameters of the iron core (Wcu, Hcd, Wwin, and Hwin), winding (N, Nt, and Nc), insulations (wt, wL, wLL, wa, and wcc), and bare wires (wc and hc, they are the width and the thickness of the bare wire, respectively). Table 1 gives the parameters of the SEM. The 2D finite element model of the SEM is built with the parameters given in Table 1. The materials of the iron core and winding are A3 steel and aluminum, respectively. The 330 V voltage is equally divided into 360 and is assigned to each turn according to the winding mode shown in Figure  1. The natural boundary condition and the adaptive mesh generation method are adopted.
The default depth of the Maxwell software 2D electrostatic field model is 1 m (cannot be modified), and the 2D model cannot deal with the arc-shaped part at both ends of the electromagnet winding. Therefore, the parameters 1 2 0 = = R R and 1 = c m L are set to ensure the consistency between the analytical model and the finite element simulation model. Figure 8 illustrates the potential distribution of the finite element simulation of the above model. In Figure 8a, the potential distribution of the winding is from low to high, and the iron core's floating potential is 54.55 V; in Figure 8b, the winding's potential distribution is from high to low, and the iron core's floating potential is 275.45 V. It proves that the relationship between Equation (25) and Equation (27) is correct. Therefore, the potential difference between the iron core and the adjacent conductors makes part of the electric field energy stored in the insulation between them, which affects the equivalent stray capacitance of the winding. Figure 9 shows the energy density cloud chart of typical zones of the finite element simulation, in which (a), (b), (c), and (d) correspond to the A-type zone, B-type zone, C-type zone, and D-type zone mentioned in Section 2.2, respectively. There is some electric field energy in the above areas, so it is necessary to compensate for these areas when calculating the equivalent stray capacitance of the winding. The calculation results of the static capacitances, the electric field energy, the equivalent stray capacitance, and the iron core's floating potential, which are by the FEM, TAM, and CAM, are given in Table 2.   where RErr is the relative error and is defined as Equation (33).

( )
where FEM P , TAM P , and CAM P are the parameters' value calculated by the FEM, TAM, and CAM, respectively.
From Figure 10, we can know that the errors of the static capacitances, including tL tLL cd cu , , , , C C C C and cdu C , calculated by the CAM, are less than those of the TAM. Among them, the error of cu C is the largest, which is caused by the smaller ratio of the plate width to the distance between plates of the corresponding parallel plate capacitor, which leads to the edge effect obvious.
With the proposed edge effect compensating algorithm, the calculation error of cu C reduces from 9.77% to 1.16%, which is improved. Compared to the FEM, the relative errors of the iron core floating potential calculated by the TAM and the CAM are 2.1% and 2.11%, respectively, which is caused by neglecting the capacitance between the core and non-adjacent windings (seeing the simplification of Equation (24) to Equation (25)). The errors of the equivalent stray capacitance obtained by the TAM and the CAM are 6.77% and 1.53%, respectively, which are enough to meet the requirements of engineering application.

Effects of the Insulation Thickness
To adequately prove the validity of the proposed analytic method, different models are solved by the FEM, TAM, and CAM, respectively. Subsequently, we deliver out a comparison between different results. The dimensions of the iron core remain unchanged, and the parameters  Figure 11 gives the wc-wcc and hc-wcc curves, in which wcc varies from 0.2 mm to 4 mm at intervals of 0.2 mm. Figure 11. wc-wcc and wc-wcc curves.
Three methods (FEM, TAM, and CAM) are utilized to solve 20 models of the SEM (parameters are shown in Figure 11). Figure 12 shows the calculation results. From Figure 12a,b, we can know that the values of CtL and CtLL decrease linearly with the increase in wcc, which is because wc and hc decrease linearly with the increase in wcc, so that the plate area of the parallel plate capacitors corresponding to CtL and CtLL decreases linearly. From Figure 12c-e, it is known that Ccd, Ccu, and Ccdu are inversely proportional to wcc. This is because, under the constraints of Equation (32), the area of plates and the distance between two plates of the parallel plate capacitors corresponding to Ccd, Ccu, and Ccdu are increased, and the distance between two plates is the main factor affecting the capacitance parameters. It can be seen from Figure 12f that Ue obtained by FEM increases with the increase in wcc, while the trend of the variation of Ue obtained by TAM is opposite to that of FEM, so the calculation error increases with the increase in wcc; the changing trend of Ue obtained by CAM is the same as that of FEM, and the maximum difference between them is less than 1.8 V. Figure 12g,h present the trend of variation of the electric field energy Wall and the equivalent stray capacitance Ceq, respectively: both of them decrease with the increase in wcc and gradually converge to a constant. When wcc is small, more electric field energy is stored in the insulation between the winding and iron core, and the capacitances Ccd, Ccu, and Ccdu have a significant effect on the equivalent stray capacitance; conversely, the electric field energy is mainly stored inside the winding, so the equivalent stray capacitance is mainly determined by Ctl and Ctll.
With the above analysis, both TAM and CAM can solve the equivalent stray capacitance of the electromagnet winding, but the latter is more accurate. Taking the model of 0 2 = cc . mm w as an example, the errors of equivalent stray capacitance calculated by TAM and CAM are 6.77% and 1.53%, respectively.

Analytic Results of the Test Rig
Substituting the real parameters of the test rig's electromagnet  Table 3.

Experiments
This section carries out experiments to verify the accuracy of the analytical model. H.Y. Lu et al. proposed the methods of measuring the stray capacitance, including the low-frequency impedance and forced resonant frequency method [22], the external series inductance resonance method, the external parallel capacitor resonance method, and the step response method. M. Zdanowski et al. introduced the frequency scanning method [5], self-oscillation excitation method, parasitic charge measurement, and the impedance analyzer method. Yet the parasitic charge method is the most effective in measuring the stray capacitance of the SEM due to the non-linearity of the iron core. This method measures the current peak during the chopper operation to calculate the equivalent stray capacitance. At the moment the transistors turn-on or turn-off, the stray capacitance is charging by the shock current, and the charge is the integration of the current waveform. Figure 13 gives the schematic diagram utilized to measure the equivalent stray capacitance of the SEM, where G1 and G2 are the field-effect transistors; T1 and T2 are the fast recovery diodes. G1, G2, T1, and T2 constitute the maglev chopper together. In Figure 13, HCS, HPF, and LPF present the hall current sensor, the high-pass filter, and the low-pass filter, respectively; L0, R0, and Ceq are the equivalent inductance, the equivalent resistance, and the equivalent stray capacitance, respectively. Udc is the supply voltage, u(t) is the terminal voltage, icL is the terminal current (the current flowing the HCS), ic is the charging current of the equivalent stray capacitance, iL is the exciting current. P1, P2, and P3 are three data-acquisition points.
Considering that the shock current caused by the equivalent stray capacitance occurs at the moment of switches shifting, and the duration is far less than the switching period, we design the LPF and the HPF to separate two signal components. The signal sampled from P1 is the superposition of iL and ic. The signals sampled from P1 and P3 are, respectively, approximate to iL and ic.
The PI current controller is designed based on the Digital Signal Processor TMS320F28335. The supply voltage is 50 V, the switching frequency is 10 kHz, and the target current is 10 A. The stable current signal waveforms are shown in Figure 14.  Figure 14a shows the signal waveform of the HCS, which proves that the shock current is generated at the moment the chopper switches action. Figure 14b shows the signal waveform of the low-pass filter, in which the average current value is stable near the target current, and the ripple current amplitude is 0.05 A. Figure 14c,d show the signal waveform of the high-pass filter, among which Figure 14d shows the current transient response waveform of once switching action. From Figure 14d, we can know that the shock current is a process of damping oscillation, and the oscillation period is 0.25 us, which is about 0.25% of the switching period. The current oscillation process completes the charging of the equivalent stray capacitance of the winding. According to the definition of capacitance, Equation (32) is obtained.
where Q is the electric charge on the plate of the capacitor, which can be obtained by time-domain integration of ic. Regarding each half period of the oscillation waveform as a triangle, Equation (32) is simplified to Equation (33).
where k i is the peak value or valley value of the current waveform. Submitting the data in Figure   14d into Equation (33), the value of equivalent stray capacitance is 3 nF. The measurement result of the equivalent stray capacitance is larger (17.35% and 11.57%) than that calculated by TAM and CAM, respectively. The reasons for the above error may be as follows: firstly, ignoring a part of the electric field energy stored in the electromagnet conductor system in the analytical method; secondly, due to the limitation of the manufacturing process, the analytical model is not the same from the actual model; thirdly, the measurement error in the experiment process.

Conclusions
This paper gives a general procedure for calculating the equivalent stray capacitance of the suspension electromagnet winding. Firstly, the TAM and CAM methods to calculate the static capacitance of suspension electromagnet winding are derived. The static capacitance parameters calculated by the CAM are closer to FEM simulation results than that of the TAM, which proves the effectiveness of the proposed edge compensation method that is based on the virtual capacitor construction. Then, according to the electrostatic field partial capacitance theory, the analytical method to calculate the floating potential of the electromagnet's iron core is derived. The error of the iron core's floating potential calculated by the TAM increases with the increase in the insulation thickness wcc (between the winding and iron core). However, the iron core's floating potential error calculated by the CAM is less than 1.8 V with the variation of wcc. Finally, based on the equivalent of electric field energy storage, the analytical calculation formula of the winding port equivalent capacitance is derived. Compared with FEM simulation results, the calculation errors of TAM and CAM are, respectively, 6.77% and 1.53%, and the relative error of the latter with experimental results is 11.57%. Overall, the analytical method proposed in this paper quickly estimates the static capacitances of the winding, the floating potential of the iron core, the electrostatic field energy stored in the winding, and the equivalent capacitance of the winding port, under the conditions of known suspension electromagnet structure parameters and the winding's potential distribution. The above work provides a theoretical basis for the winding structure optimization and winding port differential mode inductance selection, which is for suppressing the EMI interference of the maglev controller.
Meanwhile, many problems still need to be solved. In the future, more works will be delivered to verify the above results. The EMI caused by the current shock, owing to the equivalent stray capacitance, needs to be further studied and measured. The optimization method of winding equivalent stray capacitance needs to be researched based on the results of this paper. Moreover, the iron core of the high-speed maglev train is composed of mutually insulated laminations. In this case, the calculation method of core potential and stray capacitance needs to be further studied.