Optimal Velocity Control for a Battery Electric Vehicle Driven by Permanent Magnet Synchronous Motors

The permanent magnet synchronous motor (PMSM) has high efficiency and high torque density. Field oriented control (FOC) is usually used in the motor to achieve maximum efficiency control. In the electric vehicle (EV) application, the PMSM efficiency model, combined with the EV and road load system model, is used to study the optimal energy-saving control strategy, which is significant for the economic operation of EVs.With the help of GPS, IMU, and other information technologies, the road conditions can be measured in advance. Based on this information, the optimal velocity of the EV driven by PMSM can be obtained through the analytical algorithm according to the efficiency model of PMSM and the vehicle dynamic model in simple road conditions. In complex road conditions, considering the dynamic characteristics, the economic operating velocity trajectory of the EV can be obtained through the dynamic programming (DP) algorithm. Simulation and experimental results show that the minimum energy consumption and global energy optimization can be achieved when the EV operates in the economic operation area.


Introduction
The permanent magnet synchronous motor (PMSM) has high efficiency and high torque density and is widely used in servo, wind power, and electric vehicle applications.Field oriented control (FOC) or vector control is used in PMSM, and the maximum toque/current control can be achieved.The optimal efficiency of PMSM has been previously studied [1][2][3][4].In many applications, the PMSM and its load constitute a system.Considering the load characteristics, how to achieve the optimal control performance of the system has rarely been reported.In this paper, the optimal control algorithm for the system is proposed and used in a battery electric vehicle (EV), especially in automatic cruise conditions.
The battery is the only energy source in the EV.The study of the minimum energy consumption control strategy is of great significance due to the limitations of battery energy.With the development of information technology and sensor technology, road condition information can be obtained more accurately, which provides the basic conditions for the global energy optimization of the EV.Global optimization is needed to minimize the power loss based on the characteristics of the entire driving cycle.The dynamic process of the system state in the optimization algorithm is taken into account, and the greatest optimal potential is obtained.There are several methods for solving the global optimization problems, such as the dynamic programming (DP) method [5,6], the offline optimization algorithm based on neural network and genetic algorithm [7], and the iterative algorithm based on optimal control theory [8].Among them, the DP algorithm has been applied in hybrid electric vehicles because it is simple and easy to understand [5].A lookahead control for heavy trucks to minimize trip time and fuel consumption is presented in [9].A dynamic programming algorithm is devised and used in a predictive control scheme by constantly feeding the conventional cruise controller with new set points.The velocity trajectory optimization in hybrid electric trucks has been proposed in [10].Taking advantage of satellite navigation, together with the vehicle's current mass and road load parameters, an optimization velocity is formulated and solved for a driver-defined time constraint.In this paper, the analytical and DP algorithms of minimum energy consumption are explored to solve the velocity optimization of an EV driven by PMSM.The dynamic characteristics in the minimum energy consumption control strategy for EV can be accurately acquired by the DP algorithm, but the solving process is complex and time-consuming.An analytical algorithm to solve the minimum energy consumption problem under different road loads is proposed and compared with the DP algorithm.It proved that the analytical algorithm can satisfy the engineering applications.The economic operating velocities for the EV in simple road conditions and complex road conditions solved by the DP algorithm are studied.Then, the economy operating states for the EV under different speed and time constraints in complex road conditions are analyzed.Finally, the energy-optimal driving velocity reference is provided for the driver to minimize the battery energy and maximize vehicle mileage.

Mathematical Problems in Engineering
The paper is organized as follows.Section 2 presents the vehicle and electric drive system model.The optimal velocity curves in simple road conditions calculated by the analytical model are described in Section 3.Then, the optimal velocity curves without constraint, with velocity constraint, and with time constraint in complex road conditions are studied and compared in Section 4. Simulation results and experiment validation analysis are discussed in Section 5. Section 6 presents the conclusions.

Vehicle Resistance Equation.
The force needed when driving on a sloped road in an EV can be expressed as [11] where   is the motor output torque;   is the ratio of the transmission;  0 is the ratio of the reducer;   is the efficiency of the transmission system;  is the wheel radius;  is the mass of the vehicle and payload;  is the rolling resistance coefficient;  is the slope angle;   is the dimensionless coefficient;  is the cross-sectional area exposed to flow; V is the vehicle velocity;  is the vehicle rotating mass conversion factor, which is involved in the determination of the added mass moment of the rotating part in the acceleration resistance; dV/d is the vehicle acceleration; and  is the air density and is equal to 1.2258 N⋅s 2 ⋅m −4 .

Efficiency Model of PMSM.
The PMSM is a rotating electric machine where the stator is a classic three-phase stator like that of an induction motor and the rotor has surface-mounted or interior permanent magnets.A PMSM is driven by sine wave voltage coupled with the given rotor position.The generated stator flux together with the rotor flux, which is generated by a rotor magnet, defines the torque and thus speed of the motor.When the motor is running, it has four kinds of watts losses: the copper loss, the iron loss, the mechanical loss, and the stray loss.
Considering the iron loss, the and -axis equivalent circuit model of PMSM [12,13] is shown in Figure 1.The axis stator current   and -axis stator current   are divided by the iron loss currents   ,   and the torque currents   ,   .In the steady-state case, the voltage balance equation can be expressed as The flux equation can be expressed as The electromagnetic torque equation is where   and   are the and -axis terminal voltages, respectively;   and   are the and -axis armature currents, respectively;   and   are the equivalent and -axis iron loss currents, respectively;   and   are the equivalent and -axis torque currents, respectively;   and   are the and -axis stator flux-linkages, respectively;   is the magnet fluxlinkage; and   and   are the and -axis inductances, respectively. is the electrical angular velocity,   is the armature resistance, and  is the number of pole pairs.The surface-mounted permanent magnet synchronous motor (SPMSM) will be studied in this paper, and the method is also adapted to study the interior PMSM (IPMSM) and induction motor (IM).In the SPMSM, the -axis inductance is approximately equal to the -axis inductance.Assuming   =   = , the torque equation can be simplified as (5).The torque generated by the SPMSM is proportional to the axis current and has no relationship with the -axis current: Figure 2 is the power flow diagram of the PMSM in driving mode.The input power of the PMSM can be deduced from ( 2) and ( 3) and the efficiency model in Figure 1 and is expressed as where  in is the input power and   is the equivalent iron consumption resistance, which can be obtained by experiment [13].In the right side of the equation, the first part is the copper loss  Cu , the second part is the iron loss  Fe , and the third part is the electromagnetic power   .The electromagnetic power, which is the sum of the mechanical loss   , stray loss   , and output power  out , can be expressed as The maximum torque/current control can be realized by controlling the -axis current   = 0 for the SPMSM, which is applied to the traction motor of this paper.
According to the equivalent circuit and formulas (2)-( 6), the relations of power loss, electromagnetic torque, and the electric angular velocity can be expressed as The input power can be deduced as The maximum electromagnetic torque of the PMSM in driving mode needs to satisfy the following conditions: The available maximum electromagnetic torque in regenerative braking mode must satisfy the conditions According to the model of the PMSM, field oriented control (FOC) is used in the motor to achieve maximum control efficiency.In field oriented control, the stator phase currents are measured and converted into a corresponding complex vector.This current vector is then transformed to a coordinate system rotating with the rotor of the machine.Now the real -axis component of the stator current vector (  ) in this rotor flux oriented coordinate system can be used to control the rotor flux-linkage and the imaginary -axis component (  ) can be used to control the motor torque.

Efficiency Model of the Inverter.
To accurately evaluate the inverter efficiency, one must have well-developed component models and an inverter switching algorithm.Otherwise, an accurate evaluation can only be obtained by actual tests.However, by making some assumptions and simplifying the device models, the analytical model can be applied to the inverter efficiency evaluation [14][15][16].For conduction loss evaluation, a simplified device model is employed: a pure resistor for power MOSFETs and a voltage source in series with a resistor for insulated gate bipolar transistors (IGBTs) and diodes.In this paper, the loss of a three-phase MOSFET-based full-bridge inverter is discussed, and the space vector pulse width modulation (SVPWM) is adopted by this inverter.
The conduction losses for the MOSFET and the antiparallel diode are shown in (12), respectively.The detailed derivation of the equation can be found in [17]: where   and   represent the on-state voltage drops of MOSFET and diode, respectively;   is the diode voltage drop at the zero-current condition;   and   are the resistive elements of MOSFET and the diode;  is the device current; and   is the per unit fundamental stator voltage,  is the load power factor angle, and   is the peak value of the sinusoidal wave.
Considering the dead time effect, the conduction loss of the MOSFETs and the antiparallel diodes in the three-phase inverter can be simplified as where   is the dead time and   is the period time of the SVPWM.
The conduction loss of the inverter is The switching loss of the three-phase inverter is where  sw-on and  sw-off are the turn-on loss and the turn-off loss, respectively.The total loss of the three-phase inverter, which is the sum of the conduction loss and the switching loss, can be expressed as where the relationship between   and   ,   is as follows: The loss of the inverter can be expressed by the electromagnetic torque and electrical angular velocity of the PMSM as 2.4.Battery Model.The electric model of the battery, which is described by a circuit that is composed of the basic elements, such as a resistor and capacitor, and is called an ECM, is widely used to analyze the dynamic properties of the battery voltage and current [18][19][20].The error in analyzing the battery charge and discharge properties does not exceed 5% [20], which satisfies the engineering requirement.Typically, the architecture of the circuit is composed of a fundamental ohmic resistor and one or more RC networks connected in series to simulate both the transient and steady responses of the battery.One of the ECMs used to simulate cell performance is illustrated in Figure 3, where an ohmic resistor with resistance   , an RC network (  //  ), and a DC source with voltage V oc (OCV) that is a function of SOC are connected in series.V  is defined as the battery terminal voltage and   is the outflow current.Using Kirchhoff 's law, the dynamics of the ECM shown in Figure 3 can be expressed as where V  is defined as the voltage across the RC network, as observed in Figure 3.
The energy consumption in a battery can also be expressed by the electromagnetic torque and electrical angular velocity of the PMSM as where  in and  inv are functions of the electromagnetic torque and electrical angular velocity of the PMSM, respectively, as shown in ( 9) and ( 18).

Optimal Velocity Curve in Simple Road Conditions
The simple road conditions proposed in this paper mean that the rolling resistance coefficient, the air density of fluid, and the road grade angle are constant.The economic velocity of the EV is defined as the best mileage that can be reached with the constant battery energy if the EV is traveling at this speed.Determination of the energy-time optimal velocity trajectories requires a combination of route information and vehicle conditions.If the sum of the rolling resistance force and the climbing resistance force is greater than zero, the motor works in driving mode.The optimal mileage is described by where  bat is the sum of  in ,  inv , and   and   is the available discharged energy in the battery.
For simple road conditions, the mileage equation in driving mode is given by where   is the velocity in km⋅h −1 and In ( 23), if the required electromagnetic torque   calculated by the driving conditions is larger than the maximum electromagnetic torque  edm in (10), the required electromagnetic torque is equal to the maximum electromagnetic torque.
If the sum of the rolling resistance force and the climbing resistance force is less than zero, the motor works in braking mode.The optimal mileage is defined as the battery being fully charged in the shortest mileage and is described by where   is the available charged energy in the battery.
For simple road conditions, the mileage equation in braking mode is given by In ( 27), if the required electromagnetic torque   calculated by the driving conditions is less than the negative maximum electromagnetic torque  ebm in (11), the required electromagnetic torque equals the negative maximum electromagnetic torque.
In driving mode, the optimal velocity means that the vehicle drives the longest distance at this speed while consuming the same amount of battery energy.In braking mode, the optimal velocity means that the vehicle drives the shortest distance at this speed while recharging the same amount of battery energy.
The numerical method is used to obtain the problem solution in the analytical approach.For every velocity from 0 to the maximum one, the mileage can be calculated when the battery energy is a certain value.Comparing the mileages for every velocity, we can obtain the maximum one in driving mode and the minimum one in braking mode, which are the optimal velocities.
In driving conditions, generally, the rolling resistance coefficient and air density are constants, but the slope angle often changes.Therefore, it is important to solve the optimal velocity at different slope angles.When the battery energy is a certain value, the mileage versus the velocity curve at different slope angles in a four-wheel drive electric vehicle is shown in Figure 4.At a slope angle, the EV can drive different miles for different constant velocities.The specifications of the EV and motors used in this and the following figures can be found in Section 5.
In Figure 4, the EV consumes energy at slope angles of 0%, 5%, and 10%.The percentages are the tangent values of slope angles which equal the vertical height divided by the horizontal movement distance and then multiplied by 100%.An optimal velocity can be found, and the mileage is monotonically decreasing on both sides of the optimal velocity at slope angles of 0%, 5%, and 10%.The motor recharged energy at the slope angles −10% and −5%.An optimal velocity can also be found, and the mileage is monotonically increasing on both sides of the optimal velocity.Therefore, an optimal velocity can be found for each slope angle, as shown in Figure 5.The figure shows that the greater the slope angle, the greater the velocity in downhill braking mode and the greater the slope angle, the greater the velocity in driving mode.

Optimal Velocity Curve in Complex Road Conditions
The complex road conditions proposed in this paper mean that the road grade angle changes.The road parameters can be measured with the help of the GPS navigation system and vehicle sensors.To compare the analytical approach and dynamic programming (DP) algorithm, this paper uses longdistance segment-unchanging road conditions.The road consists of five segments, whose slope angles are 0% ( 1 ), 5% ( 2 ), 0% ( 3 ), −5% ( 4 ), and 0% ( 5 ).

Analytical Approach for the Optimal Velocity Curve.
In the above-mentioned long-distance segment-unchanging road, the optimal velocity curve can be obtained by an analytical approach.For more complex road conditions, a guide for solving the optimal velocity can be provided by an analytical approach.The optimal velocity is equivalent to the least energy consumption on a particular road when travelling at a particular velocity and is described by where  bat is the battery power required in each segment,   is the discharged energy of the battery on the road, and   is the recharged energy of the battery on the road.
According to the analysis in Section 3, an optimal velocity can be found for each slope angle.The EV consumes the least energy when travelling on each slope angle road with the corresponding optimal velocity.For each optimal velocity, the motor operating states and the energy consumption can be obtained according to ( 21)-(28), as shown in Figure 6.It shows that the optimal operating velocities are different for different slope angle sections.
For different slope angles, an optimal operating point (an optimal torque and speed point) can be found in the motor efficiency map, as shown in Figure 7.The motor does not operate at the optimal efficiency point at the same speed and is under the optimal point.The reason is that the road load is considered, and high driving torque means high energy consumption.
The analytical method can accurately obtain the optimal velocity in the steady state, but it is difficult to solve the dynamic process because of the nonlinearity.the dynamic programming (DP) algorithm [21] can accurately calculate the optimal velocity at each moment.The classic DP algorithm is used in this section to solve the global optimization problems.The distance of the road  is divided into  steps, and the step length is Δ.Thus,  = Δ.

Dynamic Programming Algorithm for the
The cost function is defined as where   is the -step energy consumption and is given by where V  is the -step velocity of EV in units of m⋅s −1 and is the system state transition variable. bat is the -step battery power.
The system constraints are as follows: where   is the system control variable.
The system discrete state transition model is as follows: where   is the acceleration of the EV and is described by The solution to the dynamic programming problem can be expressed as follows.
(1) Obtain the optimal solution of the last step  * V,−1 . * V,−1 is the optimal energy consumption at the last step ( − 1 step), which is the function of the system state variable V −1 : Meanwhile, the optimal system control variable  * ,−1 corresponding to the system state variable can also be obtained: (2) Derive in reverse and obtain the optimal solution  * V, (V  ) and optimal control variable  * , (V  ) at  =  − 2, . . ., 0: where V +1 is deduced by V  according to (33).
After the backwards solution procedure is complete, the optimal objective function  * V, (V  ) and the optimal control variable  * , (V  ) for each system state can be obtained.Then, according to the initial state V 0 , combined with the system state transition equation, the optimal control variable and the corresponding system state at all steps can be solved.
According to the above equation, the optimal velocity calculated by the DP algorithm is shown in Figure 8  the slope angle is increasing, the motor torque increases in advance, and the EV begins to accelerate.When the slope angle is decreasing, the motor torque decreases in advance, and the EV begins to decelerate.The energy consumption calculated by the DP algorithm is close to the ideal result calculated by the analytical approach.
Figure 8 shows that the DP algorithm can accurately solve the system dynamic process compared with the analytical approach.The optimal velocities solved by both the analytical approach and the DP algorithm in the steady state are the same, which indicates that the analytical approach can be used to solve the optimal velocity in the steady state both rapidly and accurately.
The motor operating points calculated by DP are shown in Figure 9.In the steady state, the operating area is concentrated in three points calculated by the analytical method.In the dynamic programming (DP) algorithm, the operating area is concentrated in a few lines between the three points.The reason is that the DP algorithm can handle the dynamic process, but the analytical method does not work in this condition.Considering the road load, the solving motor operating points are also under the optimal motor efficiency lines.Thus, the optimal motor efficient operation cannot achieve the energy efficient vehicle travelling.

Dynamic Programming Algorithm for the Optimal Velocity
Curve with a Velocity Constraint.In the optimal velocity analysis, if the driver has a minimum velocity requirement, a minimum velocity constraint is added to the above considerations.The system constraint expression is changed to where V  = 20 km⋅h −1 is given as the minimum velocity limit.The solution is similar to that in Section 4.2, and the optimal velocity curve with velocity constraint solved by the DP algorithm is shown in Figure 10.
In Figure 10, the optimal velocity is the same as that calculated in Figure 8 if the velocity is higher than V  .If the   velocity is lower than V  in Figure 8, it will equal V  when the velocity constraint is added.This finding coincides with the conclusion deduced in Section 3: on both sides of the optimal velocity, the energy consumption is increasing.
The motor operating points calculated by DP with the minimum velocity constraint are shown in Figure 11.The motor operating area is concentrated in a few lines at the right side of the minimum velocity.

Dynamic Programming Algorithm for the Optimal Velocity
Curve with a Time Constraint.In the optimal velocity analysis, if the driver has a minimum time requirement, a time constraint is added to the above considerations.The system constraint expression is changed to where  is given as the travelling time.There are two system discrete states: velocity and time.It is a two-dimensional dynamic programming problem.
The system discrete state transition model is as follows: where V  is the -step velocity of EV in units of m⋅s −1 and   is the -step time.
The solution of the two-dimensional dynamic programming problem can be expressed as follows.
(1) Obtain the optimal solution of the last step  * V,−1 . * V,−1 is the optimal energy consumption at the last step ( − 1 step), which is the function of the system state variables V −1 and  −1 : Meanwhile, the optimal system control variable  * ,−1 corresponding to the system state variables can also be obtained: (2) Derive in reverse and obtain the optimal solution  * V, (V  ,   ) and optimal control variable  * , (V  ,   ) at  =  − 2, . . ., 0: where [V +1 ,  +1 ] is deduced by [V  ,   ] according to (40).
At each step of the DP algorithm, the objective function for all states [V  ,   ] under the control variable  , is calculated, which produces an unreasonable amount of system states and control variables and requires a lot of computing time.The amount of computation can be reduced by eliminating the unreasonable system states and control variables.To compare with the results in Section 4.3, the time constraint is defined as 840 s ≤  ≤ 860 s.The optimal velocity curve with the time constraint solved by the DP algorithm is shown in Figure 12, the shape of which is similar to that in Figure 8.
The motor operating points calculated by DP with the time constraint are shown in Figure 13.The operation area of the motor has more freedom than in Figure 11.Thus, it will be more energy efficient.The operating lines are in a higher velocity area compared to the results in Section 4.2 because of the shorter time constraint.

Simulation and Experimental Results
The EV used in the simulation and experiment is a four-wheel direct drive electric vehicle, as shown in Figure 14.The vehicle dynamics parameters of the EV are shown in Table 1.
The four hub motors are surface permanent magnet synchronous motors.The motor parameters are shown in Table 2.
The inverter is a three-phase full-bridge circuit, which is composed of MOSFETs and antiparallel diodes.The inverter parameters are shown in Table 3.

Simulation Analysis Based on MATLAB/Simulink.
The EV and SPMSM model is built in MATLAB/Simulink.The road conditions in the simulation are the same as those in  Section 4. The DP algorithm, 10 km⋅h −1 constant velocity, 20 km⋅h −1 constant velocity, 30 km⋅h −1 constant velocity, and 40 km⋅h −1 constant velocity are compared for motor torque and energy consumption, as shown in Figure 15.
As shown in Figure 15, the lower the velocity is, the less the electromagnetic torque is demanded.The energy consumption is not less with decreased speed, which is mainly related to the motor operating points.In the five conditions, the energy consumption when traveling the distance  is shown in Table 4.As shown in Table 4, the energy consumption is the lowest when the vehicle drives in accordance with the optimal velocity obtained by the DP algorithm.Compared to the 20 km⋅h −1 constant speed driving, the energy consumption for the DP algorithm is a 4.6% reduction.Therefore, the optimal mileage can be realized if the vehicle drives in accordance with the optimal velocity obtained by the DP algorithm.
Based on the above road conditions, the optimal velocity controls with and without constraints, including a velocity constraint (20 km⋅h −1 ) and a time constraint (840-860 s), are shown in Figure 16.The corresponding energy consumptions with a velocity constraint and a time constraint are shown in Table 5.In the case of a minimum velocity constraint of 20 km⋅h −1 , the energy consumption obtained by the DP algorithm is 791.6 kJ, which is less than that of the 10 km⋅h −1 constant speed and the 20 km⋅h −1 constant speed driving.The time consumption is 854 s, which is also less than that of the 10 km⋅h −1 constant speed and the 20 km⋅h −1 constant speed driving.This result indicates that the velocity obtained by the DP algorithm can save energy and time simultaneously.It also achieves the minimum energy consumption control.
In the case of a time constraint, the energy consumption is 766.5 kJ, and the time consumption is 848 s.Compared with that in the velocity constraint, the energy consumption in the time constraint is a 3.17% reduction because the motor operation area selection has more freedom.The test road consists of a level section (200 m) and an uphill section (150 m) with a 5% slope angle.The dynamometer is used to simulate the vehicle and the road load torque according to the vehicle dynamics parameters, and the traction motor, the inverter, and the battery are to simulate the vehicle driving system.The motor is controlled by the inverter with a reference curve of speed versus time.The dynamometer is controlled with a reference curve of torque versus time.The motor speed curve is obtained by the DP algorithm or constant velocities.The dynamometer torque curve is obtained by the road load, the vehicle dynamics parameters, and the motor or vehicle speed.They are both offline calculated.The DP algorithm, 10 km⋅h −1 constant velocity, 20 km⋅h −1 constant velocity, 30 km⋅h −1 constant velocity, and 40 km⋅h −1 constant velocity are compared for energy consumption, as shown in Table 6.
As shown in Table 6, the velocity calculated by the DP algorithm is the minimum energy consumption compared with the other four conditions.Therefore, the optimal mileage can be achieved if the vehicle is driven according to the velocity calculated by the DP algorithm.
The DP algorithm is also verified by a set of real road condition tests.The test road is on the edge of the Automobile Research Institute in Tsinghua University, as shown in Figure 18.The length of the road is 800 m.There is a crossroads at 268 m from the starting point and there is a traffic light here.The traffic light is green in [60 s, 90 s], during which the vehicle can turn right and the vehicle's maximum speed limit through the crossroads is 10 km⋅h −1 .There are two deceleration zones at a distance of 412 m and 515 m from the starting point.In order to ensure the driving comfort, the vehicle speed limit through the deceleration zone is also set to 10 km⋅h −1 .The road contains an uphill section (268 m) before the traffic lights and an uphill section (144 m), a downhill section (103 m), and an uphill section (285 m) after the traffic lights.
The optimal velocity calculated by DP algorithm and the velocities operated by drivers in the real road conditions are shown in Figure 19.The test results show that the vehicle operated by the optimal velocity reaches the crossroads with traffic lights at 61.2 s and at the speed of 10 km⋅h −1 .However, the vehicle operated by a driver is faster and arrives at the crossroads before 60 s and then stops to wait.After the traffic lights change to green, the vehicle starts to accelerate again.
In the remaining section of the road, the torque calculated by DP algorithm is smoother than that operated by a driver.The total energy consumption operated by the optimal velocity is 262.5 kJ, while the average of the total energy consumption operated by drivers is 290.9 kJ.The energy consumption is reduced about 9.8%.

Conclusion
It is significant to study the optimal velocity of the EV.
Because the efficiency model of the electric drive system is accurate, combined with the vehicle resistance equation, the driving range formula corresponding to the vehicle torque and velocity can be obtained.To solve the analytical expression, the optimal velocity is obtained, which provides the driver with the optimal velocity reference.However, in complex road conditions, the dynamic characteristics of the road load need to be considered.For known road conditions, the optimal velocity curve of the drive cycle is solved by the DP algorithm.The global energy optimization can be achieved.Compared with constant velocity driving, optimal velocity travelling not only reduces energy consumption but also saves time.The conclusion is verified by comparing the DP solved velocity and other constant velocity conditions in a four-wheel drive electric vehicle.

𝐴:
Cross-sectional area exposed to flow (m 2 )   : Dimensionless coefficient   : Equivalent capacitance of RC network in battery (F) dV/d: Vehicle acceleration (m⋅s −2 )   : Regenerative energy (J)   : D i s c h a r g e de n e r g y(J)   : Kinetic energy (J) : Rolling resistance coefficient   : S w i t c h i n gf r e q u e n c y(Hz) : Acceleration of gravity, equal to 9.8 (m⋅s −2 )   : B a t t e r yo u t p u tc u r r e n t(A)   : -axis stator current (A)   : E q u i v a l e n t-axis iron consumption current (A)   : E q u i v a l e n t-axis torque current (A)   : Ratio of transmission   : -axis stator current (A)   : E q u i v a l e n t-axis iron consumption current (A)   : E q u i v a l e n t-axis torque current (A)  0 : R a t i oo fr e d u c e r : D e v i c ec u r r e n t(A)   : Peak value of the sinusoidal wave (A)  on : Turn-on energy coefficient  off : Turn-off energy coefficient   : -axis inductance (H)   : -axis inductance (H) : M a s so fv e h i c l e(kg)   : Per unit fundamental stator voltage : Motor speed (r⋅min

Figure 2 :
Figure 2: Power flow diagram of the PMSM.

Figure 4 :Figure 5 :
Figure 4: Relations between velocity and mileage at different slope angles.

5 Figure 6 :
Figure 6: Optimal operating state for different slope angle sections.

Figure 7 :
Figure 7: The motor operating points for different slope angle sections.

Figure 8 :
Figure 8: Optimal operating conditions solved by DP.

Figure 9 :
Figure 9: The motor operating points solved by DP.

Figure 10 :
Figure 10: Optimal operating conditions with a velocity constraint.

Figure 11 :
Figure 11: The motor operating points with a velocity constraint.

Figure 12 :Figure 13 :
Figure 12: Optimal operating conditions with a time constraint.

Table 3 :Figure 15 :
Figure 15: Electromagnetic torque and energy consumption in different conditions.

Figure 16 :
Figure 16: Optimal operating conditions with and without constraints.

Figure 17 :
Figure 17: Motor and inverter efficiency map in driving and braking modes.(a) Driving mode.(b) Braking mode.

Figure 19 :
Figure 19: Test results in the real road conditions.
b e ro fp o l ep a i r s  bat : Ba t teryin p u tpo w er(W)  bat : Battery input power in time   (W)   : P o w e rc o n s u m p ti o ni nba t t e ry(W)   : C o n d u c t i o nl o s so ft h et h r e e -p h a s e inverter (W)  - : Conduction loss of antiparallel diodes in the three-phase inverter (W)  -MOSFETs : Conduction loss of the MOSFETs in the three-phase inverter (W)  Cu : Co p perlos s(W)   : E l e c t r o m a g n e t i cp o w e r(W)   : Rated electromagnetic torque (W)  Fe : I r o nl o s s(W) Optimal Velocity Curve.For the above-mentioned long-distance road, n (rpm) . The torque curve solved by the DP algorithm is very smooth.When

Table 1 :
Parameters of the electric vehicle and battery.

Table 2 :
Specifications of the surface permanent magnet synchronous motor.

Table 4 :
Consumption energy in distance s.

Table 5 :
Consumption energy with velocity and time constraints.

Table 6 :
Energy consumption on the test road.
Figure 18: Test road in Tsinghua University.