Motor Vector Control Based on Speed-Torque-Current Map

: In order to e ﬀ ectively extend the mileage of pure electric vehicles, the inﬂuence of the electromechanical energy conversion principle and the dynamic control of a permanent magnet synchronous motor (PMSM) on the performance of pure electric vehicles are studied with a dual-motor drive system as a carrier in this paper. A vector control strategy based on the speed-torque-current map of a motor is proposed, considering the bus voltage ﬂuctuation inﬂuenced by battery charge and discharge. By constructing a complete vehicle model including the dynamic control model of the motor, the power distribution control strategy of a dual-motor coupling mode considering the dynamic characteristics of the motor is developed and simulated on the MATLAB platform. The results show that the dynamic model of the motor is closer to the actual running conditions. The proposed control strategy e ﬀ ectively reduces the demand power, decreases the energy consumption of the electric drive system, and extends the mileage of the vehicle. coupling condition. They have been used to study the dynamic characteristics and energy consumption characteristics of the double motor coupling mode. The simulation results show that compared to the static control model, the dynamic control model can better reﬂect the actual working process of the motor. In the torque coupling mode, the minimum demand power of the system, corresponding to the optimal MG1 target torque, is obtained from the dynamic characteristics of the motor. This power is 1.27% lower than that obtained from the static model and 20.6% lower than the maximum demand power of the system. In the speed coupling mode, the minimum demand power of the system, corresponding to the optimal MG1 target speed, is also obtained from the dynamic characteristics of the motor. This power is 6.16% lower than the demand power corresponding to the MG1 target speed obtained by the static characteristics and 22.74% lower than the maximum demand power of the system. The above results show that the proposed vector control strategy based on the dynamic model of the speed–torque–current map achieves optimal power distribution of the electric drive system, and e ﬀ ectively reduces the energy consumption of the electric drive system and extends the mileage of the car. This research results in this article lay a theoretical foundation for further improving the dynamic performance of electric vehicles.


Introduction
As the only source of propulsion, the dynamic characteristics and the performance of the motor are very important to vehicles. In order to meet the needs of vehicle research, the study of motors from static control strategy to dynamic control strategy can make the vehicle run more efficiently. The efficiency and torque characteristics of the motor directly determine the effect of the control strategy. It also has a direct impact on the power loss of the motor and the mileage of the vehicle.
In recent years, the permanent magnet synchronous motor (PMSM) is increasingly used as the power source in pure electric vehicles. Multi-power source systems are also taken seriously because of their wide efficiency range and multi-mode flexibility [1]. At present, the drive control strategy of a multi-power driven vehicle is mainly designed based on the static characteristics of the transmission system components [2]. In [3], a rule-based control strategy for a dual-motor coupling driven pure electric vehicle is designed to obtain the optimal acceleration performance. In [4], the dynamic programming algorithm is used to optimize the mode switching threshold value of a dual-motor coupling electric drive system with rotational speed coupling mode, and the optimized control strategy significantly improves the economic efficiency of the vehicle. Wang [5] comprehensively considered the problem of motor temperature rise, which leads to its efficiency and torque decline in the process of designing the coordinated control strategy of the whole vehicle powertrain. It effectively reduces the frequency of motor failure due to the temperature rising. Li [6] investigated the PMSM control system through combining the method of rotor field-oriented control and the space vector control strategy. Zhou [7] established the dynamic model of the PMSM and the efficiency model of the inverter

C2
Final Drive Output shaft   The main parameters of the system are shown in Table 2.

Motor Model and Operation Constraints
The motor MG1 is taken as an example for modeling. The PMSM with three-phase and four-pair poles is selected, whose parameters are shown in Table 3. There is always power loss in the working process of a motor, which leads to the fluctuation of speed and torque. The iron loss model can not only reflect the changes of motor voltage, current, flux, and torque, but can also express the changes of the controllable loss including copper loss and iron loss caused by the winding current, which is widely used in the PMSM loss modeling [17]. The equivalent circuit diagram of the PMSM dynamic model employed in this paper is shown in Figure 3. Considering the loss caused by the stator eddy current and hysteresis effect, an iron loss resistance R c (ω) is added to the equivalent circuit in parallel.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 4 of 22 The main parameters of the system are shown in Table 2.

Motor Model and Operation Constraints
The motor MG1 is taken as an example for modeling. The PMSM with three-phase and four-pair poles is selected, whose parameters are shown in Table 3.

Symbol
Value There is always power loss in the working process of a motor, which leads to the fluctuation of speed and torque. The iron loss model can not only reflect the changes of motor voltage, current, flux, and torque, but can also express the changes of the controllable loss including copper loss and iron loss caused by the winding current, which is widely used in the PMSM loss modeling [17]. The equivalent circuit diagram of the PMSM dynamic model employed in this paper is shown in Figure  3. Considering the loss caused by the stator eddy current and hysteresis effect, an iron loss resistance R (ω) is added to the equivalent circuit in parallel.  The voltage equivalent equations of the dynamic model are as follows.
The iron loss resistance R c is calculated by [18]: Due to the limits of battery voltage, temperature, the bearing capacity of components, and PWM (pulse width modulation) modulation methods, the current and voltage supply capacity of the inverter is limited respectively to i lim and U dc_lim . The motor armature current i s and terminal voltage u s expression and constraints are as follows: The expressions of electromagnetic torque T e and the motion equation of the PMSM are as follows: The relationship between ω r and ω is as follows:

Speed-Torque-Current Map Integrated Current Control
At present, the control strategies of the motor current mainly include the maximum torque current ratio (MTPA) control, maximum current voltage (MCMV) control, and maximum torque voltage ratio (MTPV) control.
MTPA control is a control method aiming to obtain the maximum output torque with the minimum input of a stator current vector within the limits of voltage and current [19,20]. The control equation of its boundary line is as follows: MCMV control is one of the weak magnetic control methods. Its characteristic is that the current vector and voltage vector are the limits of the current and the voltage of the inverter bus, Appl. Sci. 2020, 10, 78 6 of 21 respectively [21]. When the motor speed is adjusted to a certain speed of ω m , the control equation of (i q , i d ) in this region is as follows: MTPV control enables the motor to achieve the maximum speed by making full use of the voltage vector supplied by the inverter [22][23][24]. For the tab-mounted PMSM, the reference current vector and the voltage control equation of MTPV is as follows: To sum up, the running trajectories of the boundaries of the three control strategies in the d-q coordinate system are shown in Figure 4. OA, AB, and BC represent the running trajectories of the three control strategies (MTPA, MCMV, and MTPV) in the d-q coordinate system, respectively.
vector and voltage vector are the limits of the current and the voltage of the inverter bus, respectively [21]. When the motor speed is adjusted to a certain speed of ω , the control equation of (i , i ) in this region is as follows: MTPV control enables the motor to achieve the maximum speed by making full use of the voltage vector supplied by the inverter [22][23][24]. For the tab-mounted PMSM, the reference current vector and the voltage control equation of MTPV is as follows: To sum up, the running trajectories of the boundaries of the three control strategies in the d-q coordinate system are shown in Figure 4. OA, AB, and BC represent the running trajectories of the three control strategies (MTPA, MCMV, and MTPV) in the d-q coordinate system, respectively.  To analyze the influence of different current control strategies on the motor speed-torque characteristics, the external characteristics under different strategies are obtained, as shown in Figure 5. To analyze the influence of different current control strategies on the motor speed-torque characteristics, the external characteristics under different strategies are obtained, as shown in Figure 5. Appl. Sci. 2020, 10,   Aiming to combine the advantages of the above three control strategies, a comprehensive control strategy is proposed in this paper to expand the motor working range to its maximum within the scope of the inverter bus capacity. The principles of the comprehensive control strategy are as follows: (a) The MTPA control strategy is adopted in the peak torque area to minimize the copper loss and iron loss. (b) In the high-torque region of the weak magnetic field, the MCMV control is adopted to transit to the MTPV control range. (c) MTPV control is used to extend the speed range of the motor in the high-speed region of the weak magnetic field.
With the above principles, the external characteristic boundary of current control is obtained, as shown in Figure 6. Compared with the single current control strategy, the designed comprehensive control strategy has a wider speed range under the same torque, which improves the utilization ratio of the input current vector and voltage vector and maximizes the overall working range of the motor.  Aiming to combine the advantages of the above three control strategies, a comprehensive control strategy is proposed in this paper to expand the motor working range to its maximum within the scope of the inverter bus capacity. The principles of the comprehensive control strategy are as follows: (a) The MTPA control strategy is adopted in the peak torque area to minimize the copper loss and iron loss. (b) In the high-torque region of the weak magnetic field, the MCMV control is adopted to transit to the MTPV control range. (c) MTPV control is used to extend the speed range of the motor in the high-speed region of the weak magnetic field.
With the above principles, the external characteristic boundary of current control is obtained, as shown in Figure 6. Compared with the single current control strategy, the designed comprehensive control strategy has a wider speed range under the same torque, which improves the utilization ratio of the input current vector and voltage vector and maximizes the overall working range of the motor.  Aiming to combine the advantages of the above three control strategies, a comprehensive control strategy is proposed in this paper to expand the motor working range to its maximum within the scope of the inverter bus capacity. The principles of the comprehensive control strategy are as follows: (a) The MTPA control strategy is adopted in the peak torque area to minimize the copper loss and iron loss. (b) In the high-torque region of the weak magnetic field, the MCMV control is adopted to transit to the MTPV control range. (c) MTPV control is used to extend the speed range of the motor in the high-speed region of the weak magnetic field.
With the above principles, the external characteristic boundary of current control is obtained, as shown in Figure 6. Compared with the single current control strategy, the designed comprehensive control strategy has a wider speed range under the same torque, which improves the utilization ratio of the input current vector and voltage vector and maximizes the overall working range of the motor.  Based on the designed current control boundary curve, the speed and torque are traversed and substituted into Equations (8) and (9) to obtain the current of the d axis and q axis of all the work points. The resulting speed-torque-current map of the d-axis and q-axis are shown in Figure 7.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 8 of 22 Based on the designed current control boundary curve, the speed and torque are traversed and substituted into Equations (8) and (9) to obtain the current of the d axis and q axis of all the work points. The resulting speed-torque-current map of the d-axis and q-axis are shown in Figure 7.

Motor Vector Control Strategy
In order to solve the nonlinear problem of the PMSM current in the weak magnetic region, a vector control strategy is designed for the motor based on the speed-torque-current map, as shown in Figure 8. When the motor runs in the weak magnetic region, the d and q-axis currents are obtained using the lookup table method. This method avoids the control problems caused by the nonlinear current, which in turn extends the weak magnetic working region of the motor. It is convenient to obtain the current in the deep weak magnetic region by using the lookup table. Besides, the following problems can also be solved: it is difficult to obtain the current simply by using the weak magnetic PI controller, and the dynamic response of the motor is slow.

Motor Vector Control Strategy
In order to solve the nonlinear problem of the PMSM current in the weak magnetic region, a vector control strategy is designed for the motor based on the speed-torque-current map, as shown in Figure 8. When the motor runs in the weak magnetic region, the d and q-axis currents are obtained using the lookup table method. This method avoids the control problems caused by the nonlinear current, which in turn extends the weak magnetic working region of the motor. It is convenient to obtain the current in the deep weak magnetic region by using the lookup table. Besides, the following problems can also be solved: it is difficult to obtain the current simply by using the weak magnetic PI controller, and the dynamic response of the motor is slow. Appl. Sci. 2020, 10, x FOR PEER REVIEW 9 of 22  As shown in Figure 8, a speed-torque-current map link is added to the vector control strategy between the speed control link and the current control link. The output reference torque T _ of the speed loop and the actual motor speed are the inputs to this link. The reference currents of the d and q axes corresponding to the current speed and the target torque can be obtained by a lookup table.
Since the speed-torque-current map contains the information of the weak magnetic region (including the deep weak magnetic region), the currents of the d and q axes are constrained in the voltage limit circle and current limit circle. It solves the transition problems of current control in different working ranges and the problem of torque fluctuation caused by current or voltage over-saturation.
Charging and discharging are always present in the driving process of pure electric vehicles, which leads to the change of battery SOC (state of charge). The relationship between the battery voltage and SOC is shown in Figure 9. The change of the battery voltage causes the dynamic change of bus voltage amplitude [25], which in turn affects the d-axis current and q-axis current, as well as the speed and torque characteristics of the motor. The external characteristics of the motor under different bus voltages are shown in Figure 10. As shown in Figure 8, a speed-torque-current map link is added to the vector control strategy between the speed control link and the current control link. The output reference torque T _ref of the speed loop and the actual motor speed are the inputs to this link. The reference currents of the d and q axes corresponding to the current speed and the target torque can be obtained by a lookup table. Since the speed-torque-current map contains the information of the weak magnetic region (including the deep weak magnetic region), the currents of the d and q axes are constrained in the voltage limit circle and current limit circle. It solves the transition problems of current control in different working ranges and the problem of torque fluctuation caused by current or voltage over-saturation.
Charging and discharging are always present in the driving process of pure electric vehicles, which leads to the change of battery SOC (state of charge). The relationship between the battery voltage and SOC is shown in Figure 9. The change of the battery voltage causes the dynamic change of bus voltage amplitude [25], which in turn affects the d-axis current and q-axis current, as well as the speed and torque characteristics of the motor.  As shown in Figure 8, a speed-torque-current map link is added to the vector control strategy between the speed control link and the current control link. The output reference torque T _ of the speed loop and the actual motor speed are the inputs to this link. The reference currents of the d and q axes corresponding to the current speed and the target torque can be obtained by a lookup table.
Since the speed-torque-current map contains the information of the weak magnetic region (including the deep weak magnetic region), the currents of the d and q axes are constrained in the voltage limit circle and current limit circle. It solves the transition problems of current control in different working ranges and the problem of torque fluctuation caused by current or voltage over-saturation.
Charging and discharging are always present in the driving process of pure electric vehicles, which leads to the change of battery SOC (state of charge). The relationship between the battery voltage and SOC is shown in Figure 9. The change of the battery voltage causes the dynamic change of bus voltage amplitude [25], which in turn affects the d-axis current and q-axis current, as well as the speed and torque characteristics of the motor. The external characteristics of the motor under different bus voltages are shown in Figure 10. The external characteristics of the motor under different bus voltages are shown in Figure 10. It is shown in Figure 10 that when the bus voltage decreases, the control boundary and the maximum torque boundary of the motor move to the left, which reduces the motor working range. If the influence of bus voltage variation is ignored, the operation of these areas cannot be controlled. When the voltage increases, the bus voltage cannot be fully utilized. When the voltage decreases, the working point may be outside the allowable motor working range, which leads to the over-current or over-voltage protection. The influence of bus current on the motor can also be seen from the perspective of current control of d and q axes, which is shown in Figure 11. When the bus voltage is U1, the motor works at points A1 and A2 to realize torques T1 and T2, respectively. If the bus voltage increases from U1 to U2, then the motor working points are B1 and B2, which improves the utilization ratio of the bus voltage. On the contrary, if T1 and T2 are the required torques, the motor works at points B1 and B2. If the bus voltage decreases from U2 to U1 and the control strategy still takes B1 and B2 as working points. In this situation, the current regulator saturates and loses its ability to regulate the current. It is shown in Figure 10 that when the bus voltage decreases, the control boundary and the maximum torque boundary of the motor move to the left, which reduces the motor working range. If the influence of bus voltage variation is ignored, the operation of these areas cannot be controlled. When the voltage increases, the bus voltage cannot be fully utilized. When the voltage decreases, the working point may be outside the allowable motor working range, which leads to the over-current or over-voltage protection. The influence of bus current on the motor can also be seen from the perspective of current control of d and q axes, which is shown in Figure 11. When the bus voltage is U1, the motor works at points A1 and A2 to realize torques T1 and T2, respectively. If the bus voltage increases from U1 to U2, then the motor working points are B1 and B2, which improves the utilization ratio of the bus voltage. On the contrary, if T1 and T2 are the required torques, the motor works at points B1 and B2. If the bus voltage decreases from U2 to U1 and the control strategy still takes B1 and B2 as working points. In this situation, the current regulator saturates and loses its ability to regulate the current.  It is shown in Figure 10 that when the bus voltage decreases, the control boundary and the maximum torque boundary of the motor move to the left, which reduces the motor working range. If the influence of bus voltage variation is ignored, the operation of these areas cannot be controlled. When the voltage increases, the bus voltage cannot be fully utilized. When the voltage decreases, the working point may be outside the allowable motor working range, which leads to the over-current or over-voltage protection. The influence of bus current on the motor can also be seen from the perspective of current control of d and q axes, which is shown in Figure 11. When the bus voltage is U1, the motor works at points A1 and A2 to realize torques T1 and T2, respectively. If the bus voltage increases from U1 to U2, then the motor working points are B1 and B2, which improves the utilization ratio of the bus voltage. On the contrary, if T1 and T2 are the required torques, the motor works at points B1 and B2. If the bus voltage decreases from U2 to U1 and the control strategy still takes B1 and B2 as working points. In this situation, the current regulator saturates and loses its ability to regulate the current.   Considering the fluctuation of bus voltage, the relationship between the actual bus voltage U dc_act and the currents of the d and q axes is as follows: Equation (15) indicates that the bus voltage fluctuation can be corrected by the motor speed. Specifically, when the bus voltage changes, the motor speed can be appropriately increased or decreased by µ times to correct the voltage, as shown in the following equation: Based on Equation (16), a speed modifier considering the influence of bus voltage is added to the speed-torque-current map link to develop a vector control strategy, as shown in Figure 12. The real-time bus voltage U _cur obtained from the voltage sensor is used as the input to the speed compensator, and variable N _cps represents the speed after the compensation.
Considering the fluctuation of bus voltage, the relationship between the actual bus voltage U _ and the currents of the d and q axes is as follows: Equation (15) indicates that the bus voltage fluctuation can be corrected by the motor speed. Specifically, when the bus voltage changes, the motor speed can be appropriately increased or decreased by μ times to correct the voltage, as shown in the following equation: Based on Equation (16), a speed modifier considering the influence of bus voltage is added to the speed-torque-current map link to develop a vector control strategy, as shown in Figure 12. The real-time bus voltage U _ obtained from the voltage sensor is used as the input to the speed compensator, and variable N _ represents the speed after the compensation. To validate the control effect of the proposed vector control strategy, a working condition is designed in which the bus voltage and load torque are 650 V and 20 Nm, respectively. Besides, the motor speed first accelerates from 0 to 1000 rpm, and then again accelerates from 1000 to 5000 rpm after the speed stabilizes for a short period (around 100 ms). The simulation results are shown in Figures 13 and 14. To validate the control effect of the proposed vector control strategy, a working condition is designed in which the bus voltage and load torque are 650 V and 20 Nm, respectively. Besides, the motor speed first accelerates from 0 to 1000 rpm, and then again accelerates from 1000 to 5000 rpm after the speed stabilizes for a short period (around 100 ms). The simulation results are shown in Figures 13 and 14.   It is seen from Figure 13 that in the process of speed regulation, torque and current fluctuate greatly due to the effect of the bus voltage changes. However, in Figure 14, the current and torque fluctuations are greatly reduced during the speed regulation period (1000-5000 rpm), when the bus voltage changes are taken into consideration. The results in these two figures indicate that the speed It is seen from Figure 13 that in the process of speed regulation, torque and current fluctuate greatly due to the effect of the bus voltage changes. However, in Figure 14, the current and torque fluctuations are greatly reduced during the speed regulation period (1000-5000 rpm), when the bus voltage changes are taken into consideration. The results in these two figures indicate that the speed modifier with bus voltage regulation well reduces the torque and current fluctuation in the speed regulation process, which is beneficial to the output speed-torque characteristics of the motor.

System Power Distribution Strategy
In the actual operation process of the motor, when the output demand is constant, the position of the working point has an impact on the speed regulation process and the energy consumption of the motor. However in the static model, the working point is well determined. Based on vehicle model and PMSM dynamic model established in the previous section, the speed-torque characteristics and the energy consumption characteristics of system are analyzed when the system demand changes. In the coupling modes, the torque distribution strategy and speed distribution strategy are shown in Figures 15 and 16, respectively.
Due to the difference between the static model and the actual operation, the minimum energy consumption point calculated by the static model may not be accurate for the motor. Therefore, the power distribution strategy based on the dynamic control enables the vehicle to operate at the working point with the minimum energy demand, thereby reducing the vehicle energy consumption.
In the actual operation process of the motor, when the output demand is constant, the position of the working point has an impact on the speed regulation process and the energy consumption of the motor. However in the static model, the working point is well determined. Based on vehicle model and PMSM dynamic model established in the previous section, the speed-torque characteristics and the energy consumption characteristics of system are analyzed when the system demand changes. In the coupling modes, the torque distribution strategy and speed distribution strategy are shown in Figures 15 and 16, respectively.   Due to the difference between the static model and the actual operation, the minimum energy consumption point calculated by the static model may not be accurate for the motor. Therefore, the power distribution strategy based on the dynamic control enables the vehicle to operate at the working point with the minimum energy demand, thereby reducing the vehicle energy consumption.

Simulation Analyses
Based on the configuration considered in this paper, two vehicle models which are configured with the static and dynamic models are established for simulation studies. The following three working conditions are employed for verification purposes: 0-100 km acceleration, torque coupling mode, and speed coupling mode.

0-100 km/h Acceleration Simulation
In the simulation of 0-100 km/h acceleration, when the vehicle speed is lower than 60 km/h, the torque coupling mode is adopted; otherwise, the speed coupling mode is adopted. The results are shown in Figure 17.

Torque Coupling Mode Performance Simulation
According to the flow shown in Figure 15, the designed simulation condition of the torque coupling mode is as follows: during the period from 0 to 0.05 s, the vehicle speed is adjusted to 18 km/h. In this process, the torque is adjusted to the required torque. After 0.05 s, the vehicle accelerates to 19.7 km/h with 1.8889 m/s 2 acceleration during 0.25 s. When the vehicle speed is 18 km/h and the acceleration is 1.8891 m/s 2 , the required torque and speed of the main reducer are 224 Nm and 720 rpm, respectively. According to the dynamic equation of the torque couple mode, the speeds of MG1 and MG2 are 1494 rpm and 2237 rpm, respectively. Under the corresponding speed constraint, the maximum electromagnetic torques of MG1 and MG2 can achieve 153 Nm and 75 Nm, respectively. The torque intervals of MG1 and MG2 are [−153 Nm, 153 Nm] and [−75 Nm, 75 Nm], respectively. Moreover, the torques of MG1 and MG2 (i.e., T and T ) are also limited by the demand torque T of the main reducer, namely: The motor electromagnetic torque is the sum of the load torque and the speed regulating torque. If the load torque is too large, the motor will lose the ability to regulate the speed. In order to meet the requirements of the output characteristics of the main reducer, the initial and target torques of MG1 are limited in the range of [−2 Nm, 148 Nm] based on Equation (17). Response delay and torque fluctuation always exist in the actual process of motor speed regulation. Figure 17 shows that due to the response delay, the 0-100 km/h acceleration time resulting from the dynamic model is 14.22 s, which is 0.12 s more than that of the static model (14.10 s). The reason is that the speed and torque response characteristics of the dynamic model are closer to the actual motor operating characteristics. Different from the step response in the mode switching process of the static model, it can be seen from the speed diagram that the time required for the motor to complete mode switching is almost 0.28 s in the dynamic model, which is not shown in the static model. Torque fluctuation always exists in the dynamic model but not in the static model. The fluctuation in the process of mode switching is more apparent, which is shown in the enlarged red circle areas for the torques of MG1 and MG2. In this process, the torque of MG1 is greater than that of the static model, and the torque of MG2 is less than that of the static model. The lines of the torque aligned with the static model after 5.57 s because in the mode switching process, the speed regulation of the two motors caused the change of motor torque. This situation is similar to the actual motor operation, indicating that the dynamic model can well reflect the actual motor speed regulation process.

Torque Coupling Mode Performance Simulation
According to the flow shown in Figure 15, the designed simulation condition of the torque coupling mode is as follows: during the period from 0 to 0.05 s, the vehicle speed is adjusted to 18 km/h. In this process, the torque is adjusted to the required torque. After 0.05 s, the vehicle accelerates to 19.7 km/h with 1.8889 m/s 2 acceleration during 0.25 s. When the vehicle speed is 18 km/h and the acceleration is 1.8891 m/s 2 , the required torque and speed of the main reducer are 224 Nm and 720 rpm, respectively. According to the dynamic equation of the torque couple mode, the speeds of MG1 and MG2 are 1494 rpm and 2237 rpm, respectively. Under the corresponding speed constraint, the maximum electromagnetic torques of MG1 and MG2 can achieve 153 Nm and 75 Nm, respectively. The torque intervals of MG1 and MG2 are [−153 Nm, 153 Nm] and [−75 Nm, 75 Nm], respectively. Moreover, the torques of MG1 and MG2 (i.e., T m1 and T m2 ) are also limited by the demand torque T o of the main reducer, namely: The motor electromagnetic torque is the sum of the load torque and the speed regulating torque. If the load torque is too large, the motor will lose the ability to regulate the speed. In order to meet the requirements of the output characteristics of the main reducer, the initial and target torques of MG1 are limited in the range of [−2 Nm, 148 Nm] based on Equation (17).
Under the same working condition, the torque of MG1 is 88 Nm according to the motor static model control strategy, which is set as the initial MG1 torque. That is, T1 = 88 Nm. In this section, the MG1 target load torque value (T1* = −2 Nm) is taken as an example to analyze the speed-torque characteristics during the transition process of two motor working points under this working condition. The dynamic response process of the two motors is shown in Figure 18. characteristics during the transition process of two motor working points under this working condition. The dynamic response process of the two motors is shown in Figure 18.  In the torque coupling mode, the torques of MG1 and MG2 can be coupled arbitrarily, as long as the output torque of the vehicle meets requirements. This means that countless combinations of the torques of MG1 and MG2 exist. The economy of each combination can be measured by the demand power of the electric drive system. According to the flow shown in Figure 16, the MG1 target torque  In the torque coupling mode, the torques of MG1 and MG2 can be coupled arbitrarily, as long as the output torque of the vehicle meets requirements. This means that countless combinations of the torques of MG1 and MG2 exist. The economy of each combination can be measured by the demand power of the electric drive system. According to the flow shown in Figure 16, the MG1 target torque is traversed to obtain the corresponding demand power of the electric drive, as shown in Figure 19.  Figure 18 demonstrates the change of the two motor speeds, as the vehicle speed is adjusted to 18 km/hat, no load condition during the period of 0-0.045 s. In the period of 0.045-0.05 s, the torques of the two motors are adjusted to the demand torque corresponding to the condition (i.e., 18 km/h vehicle speed and 1.8991 m/s2 acceleration. At the moment of 0.05 s, the motor MG1 is loaded by the target torque of −2 Nm, and the torque of MG2 is calculated from the real-time load torque and the load torque of MG1. In the torque coupling mode, the torques of MG1 and MG2 can be coupled arbitrarily, as long as the output torque of the vehicle meets requirements. This means that countless combinations of the torques of MG1 and MG2 exist. The economy of each combination can be measured by the demand power of the electric drive system. According to the flow shown in Figure 16, the MG1 target torque is traversed to obtain the corresponding demand power of the electric drive, as shown in Figure 19.  When the demand torque is constant in the torque coupling mode, the torque of MG1 is 88 Nm on the point of the allocated power obtained from the static model. In the dynamic model, the demand power of the electric drive system first decreases and then increases with the increase of the MG1 target torque, as shown in Figure 19. Note that the starting point T 1 * = −2 Nm, the final point T 1 * = 148 Nm, and the trough point T 1 * = 68 Nm are employed in the simulation studies, as shown in Table 4. The motor static model cannot reflect the loss during the transition process of the motor working point, so the optimal power distribution scheme of the two motors cannot be obtained when the power demand of the system is fixed. If the vehicle demand torque is constant, MG1's target torque is too high, which results in MG1 working near the peak torque. If the target torque of MG1 is too small, then MG2 will work near its peak torque, and vice versa. These situations will lead to an increase in system losses and lower efficiency. The transformation law of demand power with MG1 target torque is shown in Figure 19, when the demand torque is constant. Table 4 shows that when the MG1 target torque is 148 Nm, the demand power of the electric power system reaches a maximum of 27.140 kW. When the target torque is 88 Nm in the static model, the average demand power of the system is 21.822 kW. When the target torque is 68 Nm (obtained from the dynamic model), the demand power of the system has a minimum of 21.549 kW, which is 1.27% lower than the average power consumption obtained by the static model, and 20.6% lower than the maximum demand power of the electric drive system.

Speed Coupling Mode Performance Simulation
In the speed coupling mode, when the required torque is constant, the speeds of both motors can change within a certain range, and the system power consumption corresponding to different speed combinations is also different. According to the flow shown in Figure 16, the simulation condition is designed as follows: from 0 to 0.05 s, the vehicle speed is adjusted to 60 km/h. In this process, the torque is adjusted to the demand torque. Then, the vehicle travels from the initial speed of 60 km/h to the target speed of 60.43 km/h, with an acceleration of 0.48 m/s 2 for 0.25 s.
Under the above conditions, the demand torque and speed of the main reducer are 75 Nm and 2430 rpm, respectively. The corresponding torques of MG1 and MG2 are 51 Nm and 24 Nm, respectively. Under the restriction of corresponding speed, the maximum speeds of MG1 and MG2 are 5600 rpm and 8500 rpm, respectively. The speed ranges of MG1 and MG2 are [−5600 rpm, 5600 rpm] and [−8500 rpm, 8500 rpm], respectively. In addition, the speeds of MG1 and MG2 are also constrained by the demand speed of final drive, which is described as follows: To meet the requirements of the main reducer output characteristics, the speed range of MG1 under this condition is [−295 rpm, 5600 rpm] based on Equation (18).
The target speed n 1 * = 3590 rpm of MG1 is taken as an example to analyze the speed-torque characteristics during the working point transfer process. The dynamic responses of MG1 and MG2 are shown in Figure 20.  As shown in Figure 20, in the period of 0-0.045 s, the speed of the two motors are adjusted to the speed when the vehicle speed is 60 km/h. In the period of 0.045-0.05 s, the torques of the two motors are adjusted to the demand torque corresponding to the condition (i.e., 60 km/h speed and 0.48 m/s 2 acceleration). At the moment of 0.05 s, the MG1 speed is adjusted to the target speed of 3590 rpm; meanwhile, the torque of MG1 increases and that of MG2 decreases as these two motors need extra torque to overcome the rotational inertia force when the speed changes.

Time (s)
According to the power distribution flow in Figure 16, all of the energy consumption corresponding to the target speed n1* of MG1 can be obtained. When the vehicle speed is 60.43 km/h, the range of n1* is [−295 rpm, 5600 rpm]. When the initial speed of MG1 is set to 3242 rpm, the target speed of MG1 traverses in the interval of [−295 rpm, 5600 rpm] to obtain the demand power of the electric drive system. The demand power of the system, corresponding to different target speeds of MG1, is shown in Figure 21. As shown in Figure 20, in the period of 0-0.045 s, the speed of the two motors are adjusted to the speed when the vehicle speed is 60 km/h. In the period of 0.045-0.05 s, the torques of the two motors are adjusted to the demand torque corresponding to the condition (i.e., 60 km/h speed and 0.48 m/s 2 acceleration). At the moment of 0.05 s, the MG1 speed is adjusted to the target speed of 3590 rpm; meanwhile, the torque of MG1 increases and that of MG2 decreases as these two motors need extra torque to overcome the rotational inertia force when the speed changes.
According to the power distribution flow in Figure 16, all of the energy consumption corresponding to the target speed n1* of MG1 can be obtained. When the vehicle speed is 60.43 km/h, the range of n1* is [−295 rpm, 5600 rpm]. When the initial speed of MG1 is set to 3242 rpm, the target speed of MG1 traverses in the interval of [−295 rpm, 5600 rpm] to obtain the demand power of the electric drive system. The demand power of the system, corresponding to different target speeds of MG1, is shown in Figure 21.
0.48 m/s acceleration). At the moment of 0.05 s, the MG1 speed is adjusted to the target speed of 3590 rpm; meanwhile, the torque of MG1 increases and that of MG2 decreases as these two motors need extra torque to overcome the rotational inertia force when the speed changes.
According to the power distribution flow in Figure 16, all of the energy consumption corresponding to the target speed n1* of MG1 can be obtained. When the vehicle speed is 60.43 km/h, the range of n1* is [−295 rpm, 5600 rpm]. When the initial speed of MG1 is set to 3242 rpm, the target speed of MG1 traverses in the interval of [−295 rpm, 5600 rpm] to obtain the demand power of the electric drive system. The demand power of the system, corresponding to different target speeds of MG1, is shown in Figure 21.

MG1 target speed (rpm)
Demand power (kW) The optimal speed of MG1 n1* obtained from the static model of the motor is 3590 rpm. In the dynamic model, it is seen from Figure 21 that with the increase of target speed n1* of MG1, the system demand power first decreases and then increases. The following points of the MG1 target speed n1* are selected in the simulation studies: the starting point of the decreasing function is −295 rpm, the trough point is 1867 rpm, and the final point of the increasing function is 5600 rpm. Table 5 shows the demand power at different MG1 target speeds. The optimal speed of MG1 n1* obtained from the static model of the motor is 3590 rpm. In the dynamic model, it is seen from Figure 21 that with the increase of target speed n1* of MG1, the system demand power first decreases and then increases. The following points of the MG1 target speed n1* are selected in the simulation studies: the starting point of the decreasing function is −295 rpm, the trough point is 1867 rpm, and the final point of the increasing function is 5600 rpm. Table 5 shows the demand power at different MG1 target speeds.  Table 5 shows that the power of MG1 increases and that of MG2 decreases, with the increase of the MG1 target speed n1*. When the optimal speed (3590 rpm) obtained by the static model is taken as the target speed of MG1, the average demand power of the drive system is 23.543 kW, which is not the minimum demand power of the electric drive system under this condition. When the MG1 target speed n1* is 5600 rpm, the maximum power demand of the electric drive system is 28.595 kW. This is because when the target speed of MG1 is 5600 rpm, the speeds of motor MG1 and MG2 need to be adjusted from 3242 rpm to 5600 rpm, and from 733 rpm to −4165 rpm, respectively. This is a large speed regulation process, which leads to power loss. When the target speed of MG1 is 1867 rpm, the minimum demand power of the electric drive system can be obtained, i.e., 22.092 kW. This minimum value is 6.16% lower than the average demand power obtained from the static model and 22.74% lower than the maximum demand power of the electric drive system.

Conclusions
In this paper, the mechanical and electrical energy conversion mechanism and dynamic control of PMSM are studied in order to analyze the impact on the performance of pure electric vehicle. Considering the influence of bus voltage fluctuation on motor performance, a vector control strategy based on the speed-torque-current map is designed. In this study, three working conditions were formulated, including the 0-100 km/h acceleration condition, speed coupling condition, and torque coupling condition. They have been used to study the dynamic characteristics and energy consumption characteristics of the double motor coupling mode. The simulation results show that compared to the static control model, the dynamic control model can better reflect the actual working process of the motor. In the torque coupling mode, the minimum demand power of the system, corresponding to the optimal MG1 target torque, is obtained from the dynamic characteristics of the motor. This power is 1.27% lower than that obtained from the static model and 20.6% lower than the maximum demand power of the system. In the speed coupling mode, the minimum demand power of the system, corresponding to the optimal MG1 target speed, is also obtained from the dynamic characteristics of the motor. This power is 6.16% lower than the demand power corresponding to the MG1 target speed obtained by the static characteristics and 22.74% lower than the maximum demand power of the system. The above results show that the proposed vector control strategy based on the dynamic model of the speed-torque-current map achieves optimal power distribution of the electric drive system, and effectively reduces the energy consumption of the electric drive system and extends the mileage of the car. This research results in this article lay a theoretical foundation for further improving the dynamic performance of electric vehicles.