Mathematical Problems in Engineering

Volume 2014 (2014), Article ID 193960, 14 pages

http://dx.doi.org/10.1155/2014/193960

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

State Key Laboratory of Automotive Safety and Energy, Tsinghua University, Beijing 100084, China

Received 3 June 2014; Revised 23 August 2014; Accepted 25 September 2014; Published 3 November 2014

Academic Editor: Kui Fu Chen

Copyright © 2014 Dongbin Lu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

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.

#### 1. 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–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 look-ahead 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.

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.

#### 2. Dynamics Modeling

##### 2.1. 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; 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; 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; is the vehicle acceleration; and is the air density and is equal to 1.2258 N·s^{2}·m^{−4}.

##### 2.2. 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 flux-linkage; 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 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 , the second part is the iron loss , and the third part is the electromagnetic power . The electromagnetic power, which is the sum of the mechanical loss , stray loss , and output power , can be expressed as

The maximum torque/current control can be realized by controlling the -axis current 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.

##### 2.3. 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–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 and 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–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 (OCV) that is a function of SOC are connected in series. 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 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 and are functions of the electromagnetic torque and electrical angular velocity of the PMSM, respectively, as shown in (9) and (18).

#### 3. 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 is the sum of , , 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 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 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.

#### 4. 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 long-distance segment-unchanging road conditions. The road consists of five segments, whose slope angles are 0% (), 5% (), 0% (), −5% (), and 0% ().

##### 4.1. 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 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.

##### 4.2. Dynamic Programming Algorithm for the Optimal Velocity Curve

For the above-mentioned long-distance road, 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, .

The cost function is defined as
where is the -step energy consumption and is given by
where is the -step velocity of EV in units of m·s^{−1} and is the system state transition variable. 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 . is the optimal energy consumption at the last step ( step), which is the function of the system state variable : Meanwhile, the optimal system control variable corresponding to the system state variable can also be obtained: (2)Derive in reverse and obtain the optimal solution and optimal control variable at : where is deduced by according to (33).

After the backwards solution procedure is complete, the optimal objective function and the optimal control variable for each system state can be obtained. Then, according to the initial state , 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 torque curve solved by the DP algorithm is very smooth. When 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.

##### 4.3. 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 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 . If the velocity is lower than in Figure 8, it will equal 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.

##### 4.4. 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 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 . is the optimal energy consumption at the last step ( step), which is the function of the system state variables and : Meanwhile, the optimal system control variable corresponding to the system state variables can also be obtained: (2)Derive in reverse and obtain the optimal solution , and optimal control variable , at : where [, ] is deduced by [, ] according to (40).

At each step of the DP algorithm, the objective function for all states [, ] 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. 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.

#### 5. 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.

##### 5.1. 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.

##### 5.2. Experimental Results

The motor and inverter efficiency model is verified first. Figures 17(a) and 17(b) show the motor and inverter efficiency calculated by the analytical model and the experimental data in driving and braking modes, respectively. The figure indicates that the motor and inverter efficiency calculated by the analytical model can reflect the actual motor efficiency.

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%.

#### 6. 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.

#### Nomenclature

Cross-sectional area exposed to flow | |

: | Dimensionless coefficient |

: | Equivalent capacitance of RC network in battery |

: | Vehicle acceleration |

: | Regenerative energy |

: | Discharged energy |

: | Kinetic energy |

: | Rolling resistance coefficient |

: | Switching frequency |

: | Acceleration of gravity, equal to 9.8 |

: | Battery output current |

: | -axis stator current |

: | Equivalent -axis iron consumption current |

: | Equivalent -axis torque current |

: | Ratio of transmission |

: | -axis stator current |

: | Equivalent -axis iron consumption current |

: | Equivalent -axis torque current |

: | Ratio of reducer |

: | Device current |

: | Peak value of the sinusoidal wave |

: | Turn-on energy coefficient |

: | Turn-off energy coefficient |

: | -axis inductance |

: | -axis inductance |

: | Mass of vehicle |

: | Per unit fundamental stator voltage |

: | Motor speed |

: | Number of pole pairs |

: | Battery input power |

: | Battery input power in time |

: | Power consumption in battery |

: | Conduction loss of the three-phase inverter |

: | Conduction loss of antiparallel diodes in the three-phase inverter |

: | Conduction loss of the MOSFETs in the three-phase inverter |

: | Copper loss |

: | Electromagnetic power |

: | Rated electromagnetic torque |

: | Iron loss |

: | Input power |

: | Total loss of the three-phase inverter |

: | Mechanical loss |

: | Output power |

: | Stray loss |

: | Turn-on loss |

: | Turn-off loss |

: | Switching loss of the three-phase inverter |

: | Armature resistance |

: | Resistive element of diode |

: | Resistive element of MOSFET |

: | Ohmic resistor in battery |

: | Ohmic resistor of RC network in battery |

: | Distance |

: | Period time of the SVPWM |

: | Dead time |

: | Time in every speed step |

: | Motor electromagnetic torque |

: | Optimum electromagnetic torque in regenerative braking |

: | Maximum braking electromagnetic torque |

: | Maximum driving electromagnetic torque |

: | Maximum electromagnetic torque |

: | Motor mechanical torque |

: | Motor output torque |

: | -axis terminal voltage |

: | -axis terminal voltage |

: | Vehicle velocity |

: | Vehicle velocity |

: | On-state voltage drops of diode |

: | Battery output voltage |

: | Voltage across the RC network in battery |

: | Final velocity |

: | Initial velocity |

: | On-state voltage drops of MOSFET |

: | Diode voltage drops at the zero-current condition |

: | Slope angle |

: | Vehicle rotating mass conversion factor |

: | Efficiency of the transmission system |

: | Air density |

: | Load power factor angle |

: | -axis stator flux-linkage |

: | Magnet flux-linkage |

: | -axis stator flux-linkage |

: | Electrical angular velocity . |

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgment

This research is supported by the US-China Clean Energy Research Collaboration: Collaboration on Cutting-Edge Technology Development of Electric Vehicle (Program of International S&T Cooperation) (no. 2010DFA72760).

#### References

- R. S. Colby and D. W. Novotny, “Efficient operation of surface mounted PM sychronous motors,”
*IEEE Transactions on Industry Applications*, vol. 23, no. 6, pp. 1048–1054, 1987. View at Publisher · View at Google Scholar · View at Scopus - N. Urasaki, T. Senjyu, and K. Uezato, “Investigation of influences of various losses on electromagnetic torque for surface-mounted permanent magnet synchronous motors,”
*IEEE Transactions on Power Electronics*, vol. 18, no. 1, pp. 131–139, 2003. View at Publisher · View at Google Scholar · View at Scopus - C. Mademlis, J. Xypteras, and N. Margaris, “Loss minimization in surface permanent-magnet synchronous motor drives,”
*IEEE Transactions on Industrial Electronics*, vol. 47, no. 1, pp. 115–122, 2000. View at Publisher · View at Google Scholar · View at Scopus - H. Kim, J. Hartwig, and R. D. Lorenz, “Using on-line parameter estimation to improve efficiency of IPM machine drives,” in
*Proceedings of the IEEE 33rd Annual Power Electronics Specialists Conference (PESC '02)*, pp. 815–820, June 2002. View at Scopus - C.-C. Lin, H. Peng, J. W. Grizzle, and J.-M. Kang, “Power management strategy for a parallel hybrid electric truck,”
*IEEE Transactions on Control Systems Technology*, vol. 11, no. 6, pp. 839–849, 2003. View at Publisher · View at Google Scholar · View at Scopus - M. Koot, J. T. B. A. Kessels, B. de Jager, W. P. M. H. Heemels, P. P. J. van den Bosch, and M. Steinbuch, “Energy management strategies for vehicular electric power systems,”
*IEEE Transactions on Vehicular Technology*, vol. 54, no. 3, pp. 771–782, 2005. View at Publisher · View at Google Scholar · View at Scopus - V. Paladini, T. Donateo, A. de Risi, and D. Laforgia, “Super-capacitors fuel-cell hybrid electric vehicle optimization and control strategy development,”
*Energy Conversion and Management*, vol. 48, no. 11, pp. 3001–3008, 2007. View at Publisher · View at Google Scholar · View at Scopus - S. Delprat, T. M. Guerra, G. Paganelli, J. Lauber, and M. Delhom, “Control strategy optimization for an hybrid parallel powertrain,” in
*Proceedings of the American Control Conference*, pp. 1315–1320, Arlington, Va, USA, June 2001. View at Scopus - E. Hellström, M. Ivarsson, J. Åslund, and L. Nielsen, “Look-ahead control for heavy trucks to minimize trip time and fuel consumption,”
*Control Engineering Practice*, vol. 17, no. 2, pp. 245–254, 2009. View at Publisher · View at Google Scholar · View at Scopus - T. van Keulen, B. de Jager, D. Foster, and M. Steinbuch, “Velocity trajectory optimization in hybrid electric trucks,” in
*Proceedings of the American Control Conference (ACC '10)*, pp. 5074–5079, Baltimore, Md, USA, July 2010. View at Scopus - H. Wallentowitz,
*Automotive Engineering I: Longitudinal Dynamics of Vehicles*, China Machine Press, 2009. - N. Urasaki, T. Senjyu, and K. Uezato, “Accurate modeling for permanent magnet synchronous motor drives,” in
*Proceedings of the 15th Annual IEEE Applied Power Electronics Conference and Exposition (APEC '00)*, pp. 387–392, February 2000. View at Scopus - S. Morimoto, Y. Tong, Y. Takeda, and T. Hirasa, “Loss minimization control of permanent magnet synchronous motor drives,”
*IEEE Transactions on Industrial Electronics*, vol. 41, no. 5, pp. 511–517, 1994. View at Publisher · View at Google Scholar · View at Scopus - J. S. Lai, R. W. Young, and J. W. McKeever, “Efficiency consideration of DC link soft-switching inverters for motor drive applications,” in
*Proceedings of the 25th Annual IEEE Power Electronics Specialists Conference (PESC '94)*, vol. 2, pp. 1003–1010, June 1994. View at Publisher · View at Google Scholar · View at Scopus - J. S. Lai, R. W. Young Sr., G. W. Ott, and J. W. McKeever, “Efficiency modeling and evaluation of a resonant snubber based soft-switching inverter for motor drive applications,” in
*Proceedings of the 26th Annual IEEE Power Electronics Specialists Conference*, pp. 943–949, Atlanta, Ga, USA, June 1995. View at Scopus - W. Dong, J.-Y. Choi, Y. Li et al., “Efficiency considerations of load side soft-switching inverters for electric vehicle applications,” in
*Proceedings of the 15th Annual IEEE Applied Power Electronics Conference and Exposition (APEC '00)*, pp. 1049–1055, New Orleans, La, USA, February 2000. View at Scopus - D. Lu and M. Ouyang, “Torque-based optimal acceleration control for electric vehicle,”
*Chinese Journal of Mechanical Engineering*, vol. 27, no. 2, pp. 319–330, 2014. View at Publisher · View at Google Scholar · View at Scopus - L. Benini, G. Castelli, A. Macii, E. Macii, M. Poncino, and R. Scarsi, “Discrete-time battery models for system-level low-power design,”
*IEEE Transactions on Very Large Scale Integration Systems*, vol. 9, no. 5, pp. 630–640, 2001. View at Publisher · View at Google Scholar · View at Scopus - L. Gao, S. Liu, and R. A. Dougal, “Dynamic lithium-ion battery model for system simulation,”
*IEEE Transactions on Components and Packaging Technologies*, vol. 25, no. 3, pp. 495–505, 2002. View at Publisher · View at Google Scholar · View at Scopus - Y.-H. Chiang, W.-Y. Sean, and J.-C. Ke, “Online estimation of internal resistance and open-circuit voltage of lithium-ion batteries in electric vehicles,”
*Journal of Power Sources*, vol. 196, no. 8, pp. 3921–3932, 2011. View at Publisher · View at Google Scholar · View at Scopus - R. Bellman,
*Dynamic Programming*, Princeton University Press, New Jersey, NJ, USA, 1957. View at MathSciNet