Research and control of a new dual-modulation magnetic gear compound motor for electric vehicles based on a mathematical model and FEA co-simulation

. According to the development of the electric vehicle motor drives, a magnetic gear compound motor with small size, lightweight, non-contact, and high-power has a good development prospect in the new energy electric vehicle industry. A new dual-modulation magnetic gear compound motor (DMFMCM) with high torque, low torque ripple, and high mechanical strength is proposed in the paper. The topology and driving principle of DMFMCM are analyzed in this paper, and the ﬁ nite elements analysis (FEA) is used to compare and analyze the DMFMCM and the conventional magnetic gear compound motor (CMGCM). Furthermore, the mathematical model of DMFMCM is established to achieve the decoupling of the three-phase voltage or current, and the mathematical model and FEA co-simulation is established to ensure the accuracy of the mathematical model. Finally, according to the principle of SVPWM controlled by i d =0, the paper simulates a PI-adjusted three-phase DMFMCM control system model. The result shows that the DMFMCM controlled by SVPWM has high stability, strong anti-interference ability and good speed regulation performance, thus meeting the development of electric vehicles.


Introduction
In order to respond to the sustainable social development trend, electric vehicles have gradually replaced traditional fuel vehicles in the automotive industry. Electric vehicles convert electrical energy into mechanical energy by driving motors, which are environmentally friendly, stable, efficient and more convenient to operate than traditional fuel vehicles [1][2][3]. In order to respond to the trend of sustainable social development, the drive motor should have the characteristics such as low speed, high torque, high stability, and good safety and comfort [4][5].
The power system of electric vehicles is generally composed of a drive motor and a gear transmission system. The driving motors of electric vehicles mainly include DC motors, induction motors, switched reluctance motors, and permanent magnet motors [6]. The principle of the DC motor is simple, and the torque control is simple. However, its volume is usually large, the energy density and power density are low, and the commutation device has a life limit. Therefore, there are certain restrictions on the use of electric vehicles [7]. The induction motor has a simple structure and high reliability. However, with its electrically exciting system, the copper and eddy current losses can lead to lower efficiency and higher control costs. [8]. Switched reluctance motors have the characteristics of a simple rotor structure, high mechanical strength and reliability. However, the switched reluctance motor has high torque fluctuations at low speeds, which will seriously affect the comfort of electric vehicles during starting, braking and low-speed driving [9]. In contrast, permanent magnet synchronous motors use permanent magnets instead of excitation winding. Permanent magnet synchronous motors reduce the copper and eddy current losses caused by the field current, so the motor efficiency is higher [10]. The permanent magnet motor is combined with the magnetic gear to form a compound motor with integrated drive and transmission. The transmission system of the magnetic gear compound motor has a transmission ratio, which can realize variable speed transmission. In addition, the magnetic gear compound motor has the advantages of small size, lightweight, low friction loss, convenient control, and high efficiency, which can meet the driving performance requirements of electric vehicles.
The magnetic gear compound motor, a new directdrive motor with low speed and high torque characteristics, was proposed by Professor K.T. Chau of the University of Hong Kong in 2009. This compound motor makes the magnetic gear and the permanent magnet motor share a rotor to form a compound structure, which improves the efficiency and power density of the device [11]. However, the compound motor has three air gaps, which are not conducive to control, and the rotor surface attaches to three layers of permanent magnets, which will fall off when the motor rotates at high speed. In 2014, Fan Ying from Southeast University in China proposed a unipolar permanent magnet arrangement magnetic gear compound motor with a self-deceleration function [12]. This structure adds a modulation ring between the rotor and the stator to achieve variable-speed transmission. In addition, the two-layer air gap reduces the magnetic leakage phenomenon, which improves motor torque. In 2017, K. Atallah from the University of Sheffield in the UK proposed a pseudo-direct drive magnetic gear compound motor [13]. The motor adopts an outer stator structure and has permanent magnets attached to the stator teeth, so there is one less air gap than the traditional magnetic gear compound motor. According to the analysis and experiment, the feasibility of the variable winding structure is verified, and the magnetic leakage phenomenon is reduced. In the same year, S.X. Niu from the Hong Kong Polytechnic University proposed a double-pole magnetic gear compound motor [14]. The outer rotor and the modulation ring have equal permanent magnets, so the two layers of permanent magnets are combined with iron sheets to form a modulation ring with modulation. Compared with the single-pole magnetic gear compound motor, this structure has a higher torque density and is more suitable for direct-drive systems with low speed and high torque. In addition, Hang Zhao from the City University of Hong Kong proposed a dualmodulator magnetic gear reduction motor in 2020 [15]. The motor structure adds a modulator outside the spoketype permanent magnet motor, which makes it have a magnetic flux focusing effect and increases the flux density on the air gap to increase the torque. In summary, it can be found that the magnetic gear compound motor will strive to simplify the mechanical structure of the compound motor, reduce the number of air gap layers, and increase the torque density and reduce torque ripple.
In order to improve the performance of the drive motor, it is not only necessary to optimize the electromagnetic structure of the motor but also to realize the control algorithm. The controller is an important part of the electric vehicle, and the control algorithm is its core elements. In recent years, motor control technology is keeping maturing, providing reliable technical support for magnetic gear compound motors. The magnetic gear compound motor is a multi-variable and strongly coupled control object. Therefore, decoupling its control parameters is the most effective way to improve the accuracy of the control algorithm, and vector control and direct torque control are the most common decoupling control strategies [16][17][18]. Although direct torque control has a low degree of model dependence, it has a longer control cycle and a smaller motor inductance. During startup and load changes, the torque ripple of the motor is relatively large. In addition, it is difficult to accurately observe the flux linkage and torque when the motor is at low speeds, and it is also hard to achieve high-performance control of the motor torque and flux linkage [19]. Vector control is also called field-oriented control, and the control effect is based on the accurate mathematical model of the controlled object. By decoupling the stator current, the vector control decomposes the current into electromagnetic torque current and excitation current, and then the two currents are independently controlled. The vector control system can ensure the torque of the actual load object, produce stable electromagnetic torque, and have a wide range of speed regulation. During the starting and braking, the armature current of the motor generates electromagnetic torque, which can make full of the overload capacity of the motor and improve the starting and braking performance [20,21]. Therefore, although the vector control relies on the accuracy of the mathematical model, the drive motor using the control method has better system stability and faster response speed. Commonly used vector control methods include i d = 0 control, cos ' = 1 control, maximum torque/current ratio control, field weakening control, and maximum efficiency control. Among them, the effect of i d = 0 control is comparable to that of DC motors, especially in the low-speed control range, the control effect is more prominent [22][23][24]. Therefore, the i d = 0 control is the most used in the vector control method.
According to the development trend of the magnetic gear compound motor, a new dual-modulation magnetic gear compound motor (DMFMCM) applied to electric vehicles is proposed. The DMFMC is to improve the torque density of the magnetic gear compound motor, suppress the torque ripple, enhance the mechanical structure, and improve the reliability and stability of the motor during operation. The magnetic gear compound motor proposed in this paper combines the modulation pieces with the stator teeth to reduce an air gap, which is convenient for the assembly of the motor. Furthermore, it can realize the secondary speed modulation and increase the gear ratio. Then, the space vector pulse width modulation (SVPWM) method with i d = 0 is used to control the speed of the DMFMCM so that the conversion of the motor speed is more stable. Based on the research of motor performance, this paper establishes mathematical model and finite elements analysis (FEA) co-simulation of DMFMCM to verify the accuracy of the mathematical model. In addition, a Simulink simulation model based on the mathematical model was established to verify the motor parameters during motor control. After ensuring the accuracy of the mathematical model, the SVPWM with i d = 0 is adopted for DMFMCM, which makes the motor have better performance to meet the needs of the electric vehicles, and build the SVPWM control module simulation model.

Analysis and research of DMFMCM
The magnetic gear compound motor is a structure that can realize low-speed and high-torque transmission, which has a good application prospect in the field of new energy such as wind power generation and electric vehicles. Figure 1 shows the topology of DMFMCM. The outer rotor of DMFMCM adopts a salient pole structure. When the convex pole coefficient d e = 0, it is a uniform permanent magnet; when the convex pole coefficient d e > 0, the permanent magnet is convex; when d e < 0, the permanent magnet is concave. The DMFMCM uses a convex pole factor of 0.5, i.e. the ratio of the protrusion of the permanent magnet to the width of the air gap is 0.5. The permanent magnets of the inner rotor have a tangential magnetized spoke-type structure with pole shoes embedded between the permanent magnets, which can form a flux-magnetic structure. In addition, the magnetic gear compound motor combines the modulated pieces and stator teethes, which makes the stator have the effect of magnetic field modulation. Therefore, the magneto-gear compound motor reduces an air gap and achieves a twice modulating effect of the stator teeth and the modulation ring. Furthermore, this salient pole of the permanent magnets and the flux-magnetic structure effectively enhance the torque density, increase the utilization of the permanent magnets and reduce the magnetic leakage. Therefore, the DMFMCM can effectively increase torque and suppress torque ripple [25,26].
The parameters of DMFMCM are shown in Table 1.

Topology of DMFMCM drive system
The driving system of DMFMCM studied in this paper is a three-phase 18-slot/28-pole field modulation permanent magnet motor (FMPMM). The three-phase alternating current is applied to the winding coil of DMFMCM, which makes the energized winding generate a rotating magnetic field to couple with the FMPMM motor rotor and generate torque. The torque generated by the rotor of FMPMM is coupled with the outer permanent magnet rotor through the outer modulation ring so that the outer rotor generates torque to work. The topological structure of the drive system FMPMM is shown in Figure 2. The rotor of FMPMM adopts a tangentially magnetized spoke structure, and the pole shoes are embedded between the permanent magnets to increase the torque. The stator adopts a slotted structure, and the stator teeth and the modulation pieces are combined to form a permanent magnet motor with field modulation. Due to the combination of the modulated block with the stator teeth, the two air gaps that would have been required with the addition of the modulating ring become one air gap. Therefore, the motor reduces the air gap and suppresses the magnetic leakage. With the introduction of the modulation ring, the air gap space in the magnetic gear section is divided into two layers of the air gap. Due to the magnetic permeability of the modulation pieces is much greater than the air permeability, the magnetic path in the radial direction of the air gap is uneven, resulting in a change in the magnetic flux density in the air gap. The magnetic field modulation of modulation ring can be expressed by the air gap permeability function where j represents the harmonics number, v p represents the modulation ring speed, ' up and ' rp respectively represent the initial radial and tangential angles of modulation ring, l 0 up , l 0 rp , l j up and l j rp represent Fourier coefficients. The permanent magnets of the inner rotor are embedded in pole-shoes and distributed along the circumference. Therefore, the pole-shoes and the modulation ring have a similar modulation effect. The specific expression is as follows where k represents the harmonics number, n s represents the pole-shoes number, v s is the pole-shoes speed, ' us and ' rs respectively represent the initial radial and tangential angles of the pole-shoes, l 0 us , l 0 rs , l k us and l k rs represent Fourier coefficients.
According to the modulating effect of the modulating commutation and pole shoe, The radial flux density distribution produced by the permanent magnet rotor is expressed as follows.
where b r j,k,a represents Fourier coefficients. The tangential magnetic flux distribution produced by any rotor is calculated as follows.
where b u j,k,a represents Fourier coefficients. The spatial harmonic pole logarithm contained in the air gap magnetic field of the magnetic gear proposed in this paper is expressed as follows.
According to the operation principle of CMG in [8], the modulated pole-pieces number and PM pole-pairs should be satisfied as (8), From the above equation, the pole-shoes and outer rotor will generate p i + n s pole-pairs harmonics. Similarly, the modulated ring and outer rotor will produce p o + n p polepairs harmonics The Fourier expansion of the magnetic induction intensity at a certain position is where B r represents the radial magnetic field component, B u represents the tangential magnetic field component. B rk and B uk respectively represent the kth Fourier coefficients of the radial and tangential component. a rk and a uk are the kth phase angles at that position. The harmonic torque of the k th order is obtained according to Maxwell tensor method as where r represents the air-gap radius, m 0 represents the vacuum permeability, and l ef represents the axis length.

The working principle of FMPMM
The electromechanical energy conversion of the magnetic gear compound motor is carried out by the magnetic field generated by the rotor and the armature winding. Based on the equivalent magnetic circuit method, the magnetic circuit can be analogous to an electric circuit to simplify the calculation and analysis process of the motor. It can be regarded as a loop composed of permanent magnets, magneto-motive force source and air gap. In order to simplify the electromagnetic torque derivation process of FMPMM, it is assumed that the magnetic field changes only in the cross-section, and the local magnetic saturation and flux leakage of the iron core are ignored. In the process of analyzing the change of air gap permeability, the equivalent air gap length d g 0 is as follows.
where d g is the length of the air gap, k c is the Carter coefficient.
Due to the permanent magnet rotor using highperformance NdFeB, the demagnetization curve changes linearly, and the characteristics of the second quadrant can be approximated as a straight line [27]. Therefore, the permanent magnet can be equivalent to a magnetomotive force source with constant outer reluctance. The magnetic circuit analysis model of FMPMM and various parameters are shown in Figure 3a. Where u m is the thickness of the permanent magnet, r g is the air gap radius, l m is the radial length of the permanent magnet, and u p is the thickness of the rotor core. According to the known parameters, the permanent magnet magnetomotive force F m can be calculated as follows.
where B r is the remanence, m 0 is the vacuum permeability, and m rm is the relative permeability. The outer magnetic resistance R m of the permanent magnet is calculated as (3).
where L a is the axial length of FMPMM, and the equivalent air gap magnetoresistance R g can be expressed as (14).
According to the magnetic flux low, it can be found that the magnetic flux of the rotor is calculated as follows.
The amplitude of air gap magnetomotive force can be expressed as follows.
Based on the calculation and analysis, the air gap magnetomotive force waveform of FMPMM can be equivalent to a square wave is shown in Figure 4. Considering the rotation of the rotor, and according to  the Fourier transform, the calculation formula for the change of the air gap magnetomotive force F ag with the circumferential position u can be obtained as follows.
where p r is the pole pairs number of the rotor, and the jth component amplitude F agj is calculated as (18).
The above is a simplified analysis of the change in permeability caused by the change in the stator tooth slot with the aid of an equivalent air gap. In the following, the effect of the stator tooth slot change on the air gap permeability is considered, and the air gap permeability is calculated. Figure 5 shows the simplified air gap permeability waveform considering the stator tooth slot alternation. The permeability l t and l s at the corresponding air gap of the stator tooth slot can be expressed as follows.
where u t is the width of the stator teeth, and u s represents the width of the slots. According to the permeance waveform and using Fourier transform, it can be seen that the calculation for the change of the air gap permeance with the position angle u is as follows.
where l 0 is the DC component, N st is the number of stator teeth, and the expression is as (22).
The harmonic component l i is calculated as follows.
In summary, it can be found that the air gap flux density of the DMFMCM is calculated as follows where B h is the high-order magnetic density harmonic component, which occupies a small proportion and it can be ignored. The formula shows that the fundamental component with the same number of rotor pole pairs, there are also N stÀpr and N st+pr pairs of harmonic magnetic fields. In a three-phase permanent magnet motor, the formula of the phase flux linkage is as follows.
where k d is the fundamental distribution coefficient of the winding, N p is the number of turns in series, and s is the coil span. Furthermore, the stator slot pitch u t can be expressed as (26).
The expression of the phase flux linkage of the DMFMCM can be obtained as shown in (27) by combining (24) and (25).
where q is the number of stator slots per pole. The induced electromotive force of FMPMM is as follows.
The FMPMM is an 18-slot/28-pole permanent magnet synchronous machine. The winding pole pair p s is 4, the gear ratio G r between the inner rotor and the winding is 3.5, the number of stator slots per pole is 2.25, and the coil span s is 2. Substituting the parameters into equation (18) above gives the phase induction potential amplitude of the FMPMM as follows.
In order to analyze the electromagnetic torque of FMPMM, using the brushless AC control method with i d = 0, the electromagnetic torque T e of the permanent magnet motor is expressed as follows.
In order to verify the accuracy of the above theory, an FEA analysis was performed on the FMPMM. The 2D model of FMPMM can be constructed using AutoCAD software and then imported into Maxwell FEA simulation software, or it can be created directly in Maxwell software. By feeding three-phase currents into the windings and simulating the model, the FMPMM air-gap flux density distribution and harmonic analysis are obtained in Figure 6. It can be found from the figure that the stator teeth will cause the air-gap flux density to generate some spatial harmonics. In addition to the fundamental components with the same number of rotor pole pairs, N stÀpr and N st+pr pairs harmonic component with larger amplitudes are generated. Therefore, it can be seen that the FEA analysis is consistent with the theoretically derived flux density expression (24). According to the equivalent circuit method of FMPMM and FEA, the working principle of the torque generated by the drive system of DMFMCM is determined to drive the entire compound motor to operate.

Comparison and analysis of DMFMCM and the conventional magnetic gear compound motor (CMGCM)
DMFMCM generates a rotating magnetic field by passing an alternating current, which couples with the permanent magnets on the spoke rotor through the stator teeth and drives the rotor to rotate. Then the rotating magnetic field is coupled with the outer rotor through the modulation ring to make the outer rotor rotate and work. In order to verify the superiority of DMFMCM, an FEA simulation model of the traditional 18-slot magnetic gear compound motor (CMGCM) was established for comparison and analysis with DMFMCM. The FEA simulation results of the outer rotor torque of the two structures are compared to analyze their performance differences.
Based on the principle of quantitative comparison, the number of permanent magnets used to the two motors is the same, and the size of the motor are the same. The FEA is used to obtain the magnetic field distribution of two compound motors as shown in Figure 7.
As shown in Figure 8, the torque of DMFMCM is significantly higher than that of CMGCM, and the torque ripple of DMFMCM is significantly lower than that of CMGCM. According to (31), the torque ripple can be obtained in Table 2. It can be found that the static torque of DMFMCM is 321.02Nm, and it is 1.4 times that of CMGCM. Furthermore, the torque ripple of DMFMCM is 0.19%, which is 29.3% of CMGCM.
where T ripple is the torque ripple, T max is the maximum torque, T min is the minimum torque, and T avg is the average electromagnetic torque.

Mathematical model of DMFMCM
The mathematical model of the drive motor is the basis of the simulation modelling and the analysis and design of the control system. Therefore, according to the electromagnetic characteristics of DMFMCM, a mathematical model under the two-phase rotating coordinate system is established. The FEA method is used to verify the results of the mathematical model to ensure its accuracy. In addition, a simulation model of DMFMCM is built by Simulink components to verify the motor parameters of the control system. The mathematical model of the three-phase static coordinate system (abc coordinate system) based on the stator armature windings of the motor is shown in Figure 9. There are three-phase symmetrical windings A, B, and C on the stator, which differ by 120°in electrical angle. d axis is on the direction of the permanent magnet pole axis of the rotor, and q is advanced by 90°in the counterclockwise direction of the rotor. Take a reference axis, the electrical angle between the d and the reference axis is u. Based on the analysis method of the mathematical model, the flux linkage equation, voltage equation and torque equation of DMFMCM are established in abc coordinate system and the two-phase rotating coordinate system (d-q coordinate system).

Mathematical model in abc coordinate system
The total flux linkage of each winding of DMFMCM includes permanent magnet flux linkage and armature reaction flux linkage. Therefore, the flux equation expression of    When the brushless AC control mode with i d = 0 is adopted, the three-phase current of the motor is expressed as follows.
where I m is the amplitude of the phase current. The current flow direction is selected as the positive direction of the current, and the positive direction of the winding flux-linkage and winding current conforms to the right-handed spiral rule. The three-phase voltage balance equation of DMFMCM is as follows.
The resistance value of each phase is equal, The electromagnetic torque formula under the electromagnetic equation is as follows.
The general formula of the mechanical motion equation of the motor is as follows.
where Tm is the total load torque of the system, J is the moment of inertia, and v is the air gap speed. The load torque is the torque of the entire transmission system, which can be interpreted as the mechanical torque obtained after the torque added to the outer rotor passes through the transmission system. Therefore, equation (36) is the mechanical equation for calculating the electromagnetic torque, which is the dynamics of a structure.
According to the mathematical model of DMFMCM in the abc coordinate system, the mathematical model in the d-q coordinate system can be established, which lays the foundation for the simulation model of DMFMCM and the research of control strategy. In addition, based on the established mathematical model, a simulation model of the DMFMCM motor is built using Simulink simulation software.

Flux linkage equation in d-q coordinate system
The inductance matrix of DMFMCM in the d-q coordinate system is as follows.
Using the FEA method to calculate inductance can be divided into two steps. Firstly, the phase flux linkage c m of the armature winding without current should be calculated; Secondly, the direct current I is applied to one phase of the armature winding, and the total flux linkage c of the phase winding is obtained. The total flux linkage is mainly composed of permanent magnet flux linkage and armature flux linkage, which can be expressed as follows.
The phase winding inductance can be calculated as follows.
The calculated inductance is self-inductance when the measured phase and the energized phase are the same phases, and the calculated inductance is mutual inductance when the measured phase and the energized phase are not the same phases. According to the method, the selfinductance and mutual inductance waveforms of phase A are shown in Figure 10.
The phase winding inductance of DMFMCM will periodically change with the change of the rotor position. According to the phase inductance calculated by the FEA method shown in Figure 10, when the rotor d-axis advances the A-phase winding axis u electrical angle, the three-phase winding inductance can be expressed as follows.
where L D is the average self-inductance of the phase, L S is the amplitude of the phase self-inductance fundamental component, and M D is the average of mutual inductance. The calculation of inductance component is as follows.
According to (41), it can be seen that there is a small cosine component in L d and L q . Because DMFMCM adopts a spoke-type inner rotor and a surface-mount outer rotor, it has a "salient pole effect". Furthermore, due to the nonuniform circumferential reluctance caused by the alternating stator slots and the salient pole rotor, a part of the flux linkages on the d-axis and the q-axis turns to each other. Therefore, the mutual inductance L dq between the d and q axes is not 0 but a sine value. Figure 11 shows the phase self-inductance and mutual inductance of DMFMCM in the abc coordinate system obtained by FEA. The inductance components L d , L q , L 0 and L dq calculated by FEA can be obtained according to (37). And the inductance component L d 0 , L q 0 , L 0 0 and L dq 0 in the d-q coordinate system derived by the mathematical model can also be calculated according to (37)-(39). The inductance waveforms calculated by FEA and the mathematical model are shown in Figure 12, and the average inductance is shown in Table 3. As shown in this table, the FEA shows that the average values of L d and L q are approximately equal, but the waveforms are different. It can be found that the DMFMCM has a salient pole effect. A comparison of the inductance waveforms in the d-q coordinate system shows that the two waveforms are similar and have small errors, indicating that the inductance results derived from the mathematical model are accurate and valid. According to the calculation of FEA and mathematical model, the results of the motor inductance in the d-q coordinate system can be obtained as shown in the following Table 3.
It can be seen from (41) that the inductance components L d and L q of DMFMCM on the d-q axis are not constant values, and the mutual inductances L dq and L qd are not zero. Therefore, the flux linkage in the d-q coordinate system can be obtained as follows.
Substituting (40) into (42), the flux linkage in the d-q coordinate system can be calculated as follows.
where c s is the fundamental wave amplitude of the phase permanent magnet chain under no-load in the abc coordinate system, where c md , c mq , c 0 is calculated as  Figure 13 shows the no-load permanent magnet flux waveform of DMFMCM. The values of c md , c mq , and c m0 can be calculated, and the flux linkages in the d-q coordinate system can also be calculated by substituting into (43) and (44).
In order to verify the accuracy of the flux linkages of the DMFMCM motor in the d-q coordinate system, the paper adopts the id = 0 control method. The flux linkage waveforms on the d-axis and q-axis are calculated by FEA and mathematical model, as shown in Figure 14. The figure shows that the waveforms calculated by FEA c d and c q are similar to the waveforms obtained by mathematical model c d 0 and c q 0 , and the average flux linkage is almost equal. Therefore, it can be found that the flux linkage of DMFMCM calculated by the mathematical model can be verified, and the accuracy of the mathematical model can be ensured. Table 4 shows the average flux linkage of DMFMCM in the d-q coordinate system calculated by the two methods. It can be found that the average flux linkages obtained by the two methods are almost equal to determining the d-axis and q-axis flux linkage components. And through co-simulation, the calculation results of the mathematical model can be verified to ensure the accuracy of the mathematical model.

Voltage equation in coordinate system
Due to the torque control adopting alternating current, the current in the d-q coordinate system obtained by Park transformation in the abc coordinate system is as follows.
where a is the current factor angle, I m is the current amplitude, and when a = 0, it is i d = 0 AC control.
where v is the rotational angular velocity of the air-gap magnetic field.

Torque equation in d-q coordinate system
According to the voltage and current equations in the d-q coordinate system, the electromagnetic torque in the d-q coordinate system can be calculated as follows.  The torque waveform of the DMFMCM obtained after using the mathematical model with i d = 0 and finite element analysis described in this paper is shown in Figure 15. By analyzing and calculating the torque using FEA and mathematical model co-simulation, the torque and torque ripple of DMFMCM can be determined more accurately. As shown in Table 5, the torque waveforms calculated by the two methods are similar, and the torque ripple is also low.
According to the torque obtained in the figure, combined with the (31), the torque ripple calculated by the two methods can be calculated. The calculation results are shown in Table 5.
Based on the mathematical model, the stator winding of the DMFMCM is modelled using Simulink. According to the needs of SVPWM, the whole model adopts a modular design. The stator winding of the DMFMCM requires output current, electromagnetic torque, and rotational angular velocity. Each module is independent and communicated through shared variables. The relationship of the modules can be intuitively seen from the Simulink simulation model of stator windings for DMFMCM and determine the relationship between the parameters. The specific simulation model is shown in Figure 16.
According to the mathematical model, Simulink is used to model the DMFMCM motor. As shown in Figure 16b, three-phase alternating current and total load torque T m are applied to the motor. After a series of motor coupling actions, three-phase AC I abc , electromagnetic torque T e and mechanical angular velocity v m can be obtained.

Judgment of voltage vector
In order to simplify the judgment process, the voltage is equivalently judged. Figure 17.
Define intermediate variables A, B, and C, with the following judgments.
-If U c > 0, then C = 1, otherwise C = 0. The sector calculation formula is as follows.
The sector of the reference voltage vector is determined by N, and the corresponding relationship is shown in Table 6.

Determination of voltage working time
According to the sector where the obtained voltage vector is located, the action time of the space vector can be calculated.
The U dc is the DC bus voltage, and the assignment relationship corresponding of T 1 , T 2 , and T 3 is obtained as shown in Table 7.
Saturation judgment is performed on the copied T 1 and T 2 . If T 1 +T 2 < T, then T 1 and T 2 keep unchanged; if T 1 +T 2 > T, the relationship is as follows.

Determination of sector vector switching point
According to SVPWM algorithm, it is defined as follows.
The relationship between the sector vector switching points T cm1 , T cm2 , T cm3 and N can be determined as shown in Table 8.

Modeling and simulation of SVPWM
According to the results obtained by the SVPWM algorithm, the triangular carrier signal is compared with the vector switching points of each sector to output a PWM pulse signal. The SVPWM is mainly to control the rotor magnetic field. In the rotating d-q coordinate system, the magnetic current and torque current are independently controlled to improve the control performance of DMFMCM. In addition, the mathematical model of DMFMCM in the d-q coordinate system is established. When i d = 0, the electromagnetic torque of DMFMCM is only related to the variable i q , thereby improving the speed regulation performance of the system. The functional block diagram of SVPWM is shown in Figure 18.
As shown in Figure 18, the control strategy includes the current loop, speed loop and SVPWM algorithm, and the SVPWM control of DMFMCM is mainly the control of the current component on the d-q coordinate system. The difference value between the reference speed and the actual speed after power-on is adjusted by the PI regulator, and output is the given value of the q-axis current i q 0 . The threephase currents i a , i b and i c from DMFMCM are transformed by Clark and Park to obtain the current components i d and i q in the d-q coordinate system. In addition, the given value i q 0 is compared with the actual value i q , and then the q-axis voltage U q 0 adjusted by the PI regulator can be output. In the same     T B  T A  T A  T C  T C  T B   T cm2  T A  T C  T B  T B  T A  T C   T cm3  T C  T B  T C  T A  T B  T A way, the d-axis goes through the same process, and the d-axis voltage U d 0 is output. U q 0 and U d 0 are transformed by Park to obtain U a and U b , and then the PWM signal is obtained by the SVPWM algorithm to control the inverter. The simulation model of SVPWM established by Simulink is shown in the Figure 19.
The conditions of the simulation model are set as follows: reference speed n 0 = 700 r/min, PWM switching frequency f = 50 Hz, simulation time is 0.75 s, variable step size is ode23tb algorithm, relative error is 1e-4. The simulation parameters of DMFMCM are shown in Table 9, the results are shown in Figure 20.    Figure 18, the three-phase DMFMCM starts without load. The speed adjusted by the PI regulator will rapidly increase from 0 and exceed the reference speed within 0-0.05 s and will quickly stabilize at the reference speed. At the same time, the torque suddenly increased and stabilized at 16.13 Nm within 0.05 s. At 0.2 s, the rotational speed changes abruptly and immediately stabilize at the reference speed. And the torque also jumped to 25.34 Nm and stabilized. In addition, changes in current are also synchronized with changes in speed and torque. It can be seen that when the electric vehicle starts and shifts, it can output stable torque in a short time, and when shifting and adjusting the speed, the output torque is quickly stable. Therefore, the simulation results show that the SVPWM control is based on PI regulation, which makes DMFMCM have the advantages of start-up and shifting time short, high stability, strong anti-interference ability, and good speed regulation performance.

Conclusion
The paper analyzes and studies a new DMFMCM used in electric vehicles. Compared with the conventional magnetic gear compound motor (CMGCM), the DMFMCM has the advantages of strong mechanical strength, high torque density, and low torque ripple. The torque of the outer rotor is 321.02Nm, and the output torque of the outer rotor is 1.4 times that of CMGCM. Furthermore, the torque ripple of the outer rotor is 0.19%, which is only 29.3% of CMGCM. Therefore, the electromagnetic performance of DMFMCM is much better than CMGCM.
The mathematical model of DMFMCM is established, and it is also modelled by Simulink software. The flux linkages equation, electromagnetic torque and voltage equation obtained based on the mathematical model are compared with the FEA simulation results of DMFMCM. According to the FEA and mathematical model cosimulation, the accuracy of the mathematical model is verified through the FEA and mathematical model cosimulation, which lays the foundation for SVPWM control.
This paper establishes a Simulink simulation model of DMFMCM and outputs the results of speed, current and electromagnetic torque based on the input parameters of DMFMCM. It provides a basis for the SVPWM control and verifies the parameters of DMFMCM. By analyzing the principle of SVPWM control, the control strategy mainly includes the speed control loop, current control loop, and SVPWM algorithm. In addition, the system adopts the AC control mode with id = 0, builds a three-phase DMFMCM motor control system model based on PI regulation, and outputs three-phase current, electromagnetic torque and speed. The simulation results show that the SVPWM control strategy makes DMFMCM have the advantages of short start-up shift time, high stability, strong antiinterference ability, and good speed regulation performance, which meets the development needs of electric vehicle drive motors.