Torque and Battery Distribution Strategy for Saving Energy of an Electric Vehicle with Three Traction Motors

: A torque and battery distribution (TBD) strategy is proposed for saving energy for an electric vehicle (EV) that is driven by three traction motors. Each traction motor is driven by an independent inverter and a battery pack. When the vehicle is accelerating or cruising, its vehicle control unit determines the optimal torque distribution of the three motors by particle swarm optimization (PSO) theory to minimize energy consumption on the basis of their torque–speed–e ﬃ ciency maps. Simultaneously, the states of charge (SOC) of the three battery packs are controlled in balance for improving the driving range and for avoiding unexpected battery depletion. The proposed TBD strategy can increase 7.7% driving range in the circular New European Driving Cycle (NEDC) of radius 100 m and 28% in the straight-line NEDC. All the battery energy can be e ﬀ ectively distributed and utilized for extending the driving range with an improved energy consumption e ﬃ ciency. the proposed TBD shows better travel distance and the higher energy consumption e ﬃ ciency than a pure PSO method or PTR strategies. Similar results were also proved on a real vehicle for hardware-in-the-loop experiments on dynamometer and road tests. The road test proved that the TBD strategy extended 11% and 23.5%


Introduction
Hybrid and pure electric vehicles (EVs) have been commercially available for many years. A lot of research that focused on multiple propulsion and energy storage systems and the related power split and energy economy strategies has become imperative issue in EVs. For EVs with only one battery pack, it is important to balance the capacity of battery cells for improving its lifetime. Li et al. [1] presented a real-time state-of-charge (SOC) calculation method for a pure EV, where the lithium battery was simulated with a second-order resistance-capacitance (RC) model and the remaining capacity in battery cells was balanced by fuzzy control through a set of bi-directional fly-back direct current-direct current (DC-DC) converters. Gallardo-Lozano et al. [2] introduced a shunting transistor method to balance battery cells during the recharging and driving modes. Huang and Abu Qahouq [3] proposed an energy sharing control scheme to regulate the DC bus voltage, and simultaneously, to balance the SOC of battery cells with micro DC-DC converters. Pham et al. [4] addressed a fast-balancing topology for lithium-ion batteries in an EV by transferring the power in high-voltage cells directly to low-voltage cells through DC-DC converters. During battery charging, Dung et al. [5] eliminated the racing phenomenon by a pulse width modulation (PWM) based equalization process among battery packs so that the charging time was reduced by 48%.
For EVs with hybrid energy storage systems (HESSs), Jin et al. [6] and Akar et al. [7], respectively, proposed for their HESS of batteries and ultracapacitors a fuzzy control-based power management strategy to reduce battery degradation. A PWM technique was introduced by Menon et al. [8] for balancing the SOC of independent battery packs by continuously regulating the power flow from Figure 2a,b provide the TNE maps from the experiments for the driving modes of the three traction motors. The braking modes were estimated from the mirror image of 75% efficiency of the driving mode. The maximum torque was 150 Nm and the maximum speed was 2400 rpm for the 15-kW front motor, and they were 122 Nm and 1200 rpm for the 7-kW rear motors. Three motor control units were responsible for driving the three motors, and three lithium-ion batteries were deployed for providing the power: one pack of 144 V, 72 Ah, and 10.45 kWh for the front drive and two packs of 72 V, 72 Ah, and 5.2 kWh for the two rear drives. Figure 2a,b provide the TNE maps from the experiments for the driving modes of the three traction motors. The braking modes were estimated from the mirror image of 75% efficiency of the driving mode. The maximum torque was 150 Nm and the maximum speed was 2400 rpm for the 15-kW front motor, and they were 122 Nm and 1200 rpm for the 7-kW rear motors. Three motor control units were responsible for driving the three motors, and three lithium-ion batteries were deployed for providing the power: one pack of 144 V, 72 Ah, and 10.45 kWh for the front drive and two packs of 72 V, 72 Ah, and 5.2 kWh for the two rear drives.

Longitudinal Vehicle Dynamics Model
When the vehicle travels in a straight line, the longitudinal traction forces are usually simplified as and . The force with subscript 1 is the force exerted on the left tire, while 2 on the right tire. The tractive force and the normal forces and of the front and rear wheels were obtained by the following equations: (1) (2) (3) (4) (5) where is the total vehicle mass; is the vehicle mass excluding tire and wheel mass; is the front tire and wheel mass; is the rear tire and wheel mass; is the gravity acceleration; is the slope angle in degrees; is the air density; is the longitudinal velocity of vehicle; and are the traction forces exerted on the front and rear tires; is the distance from mass center to front tire; is the distance from mass center to rear tire; and is the tire radius. Other vehicle specifications used in this paper are described in Table 1.

Longitudinal Vehicle Dynamics Model
When the vehicle travels in a straight line, the longitudinal traction forces are usually simplified as F xr = F xr1 = F xr2 and F x f = F x f 1 = F x f 2 . The force with subscript 1 is the force exerted on the left tire, while 2 on the right tire. The tractive force F x and the normal forces N f and N r of the front and rear wheels were obtained by the following equations: where M is the total vehicle mass; M v is the vehicle mass excluding tire and wheel mass; m w f is the front tire and wheel mass; m wr is the rear tire and wheel mass; g is the gravity acceleration; θ is the slope angle in degrees; ρ is the air density; V x is the longitudinal velocity of vehicle; F x f and F xr are the traction forces exerted on the front and rear tires; L f is the distance from mass center to front tire; L r is the distance from mass center to rear tire; and r t is the tire radius. Other vehicle specifications used in this paper are described in Table 1.

Lateral Vehicle Dynamics Model
As shown in Figure 1b, the lateral vehicle dynamics are described by a four-wheel model when the vehicle is cornering. The longitudinal, lateral, and yaw vehicle dynamic equations are expressed as: ( Appl. Sci. 2020, 10, 2653 where δ is the steer angle; F x f 1 and F x f 2 are longitudinal traction forces on the left and right front tires and these forces are assumed equal for δ = 0; F y f 1 and F y f 2 are lateral traction forces on the left and right front tires; F xr1 and F xr2 are longitudinal traction forces on the left and right rear tires; F yr1 and F yr2 are lateral traction forces on the left and right rear tires; γ is the yaw velocity; β v is the vehicle sideslip angle; I z is the mass moment of inertia in the yaw direction; and M z is the yaw moment for cornering. These equations were used to determine torque distributions when the vehicle is cornering. For simplicity, the hill climbing resistance, aerodynamic drag, and rolling resistance of the last terms in (1) are omitted, and the roll and pitch motions are neglected. However, the vehicle in CarSim for the real-time HIL simulation was modelled with self-contained yaw, roll, and pitch dynamics.

Tire Model
The CarSim tire lookup table that was obtained directly from the laboratory measurements was used to model the tire characteristics. The friction coefficient between the tire and road surface was chosen at 0.85, the longitudinal and lateral traction (or friction) forces on the tire were expressed as a function of normal force and tire slip ratio. Once the vehicle velocity, acceleration or deceleration, the normal forces N f and N r are known, the slip ratio of each wheel can be obtained. The front and rear wheel speeds ω f and ω r are therefore determined by the definition of tire slip ratio, as follows: where the tire slip ratio λ represents λ f and λ r that correspond to the speeds ω f and ω r of the front and rear wheels.

Transmission Model
After the wheel speeds and accelerations are obtained, the output torque T m f of the front traction motor, the output torque T mr provided by the right rear in-wheel motor, and the output torque T ml provided by the left rear in-wheel motor can be calculated under the following assumptions:

1.
The rotor mass is so small compared with the vehicle mass that the rotational inertias of the three traction motors are neglected.

2.
The viscous and Coulomb frictions of motors and differentials in the transmission are all neglected.
Therefore, when the vehicle accelerates: Rear wheels : where n g is the reduction ratio of gearbox; I w is the wheel inertia; ω f is the speed of the front traction motor; ω r and ω l represent, respectively, the speed of the right and left rear in-wheel motors. In the steady state for a small steer angle and by neglecting frictions on wheel motors, the yaw moment Equation (9) can be simplified as where L w is the distance between two rear wheels. This equation will be used for determining the real-time torque distribution in the PSO process when the vehicle is cornering. For a straight-line driving, M z = 0 and T mr = T ml .

Battery State of Charge Model
The SOC of each battery pack on the EV is estimated by the following equation: where SOC i is the initial state of SOC, I b is the battery current, and Q b is the battery maximum capacity. By a simple internal resistance model, the output power of battery is estimated by: where R b is the internal resistance of the battery and was set at 0.17 Ω for all the battery packs in the MIL simulations in Section 4. The open circuit voltage (OCV) V oc is a function of SOC and is obtained from experiments for each battery pack. Figure 3 illustrates the OCV curve of the front battery pack. Accordingly, the battery current I b is easily calculated from (16). Therefore, when the vehicle accelerates: where is the reduction ratio of gearbox; is the wheel inertia; is the speed of the front traction motor; and represent, respectively, the speed of the right and left rear in-wheel motors. In the steady state for a small steer angle and by neglecting frictions on wheel motors, the yaw moment Equation (9) can be simplified as (14) where is the distance between two rear wheels. This equation will be used for determining the realtime torque distribution in the PSO process when the vehicle is cornering. For a straight-line driving, and .

Battery State of Charge Model
The SOC of each battery pack on the EV is estimated by the following equation: (15) where SOCi is the initial state of SOC, is the battery current, and is the battery maximum capacity. By a simple internal resistance model, the output power of battery is estimated by: (16) where is the internal resistance of the battery and was set at for all the battery packs in the MIL simulations in Section 4. The open circuit voltage (OCV) is a function of SOC and is obtained from experiments for each battery pack. Figure 3 illustrates the OCV curve of the front battery pack. Accordingly, the battery current is easily calculated from (16).

Torque Distribution Strategy: Particle Swarm Optimization
The total torque, T t , used to accelerate the vehicle was From (14), the corresponding yaw moment, M z , in the steady state of a small steer angle was provided by the differential torque of the left and right in-wheel motors. The maximum and minimum torque ranges of the three traction motors were first determined on the basis of the TNE maps of the three traction motors, wheel angular speeds, SOC of the batteries, and current limits. PSO was then used to determine the best torque distribution for the three traction motors under the constraints of their operation ranges and Equations (14) and (17).
The PSO theory originated by observing the hunting behavior of a swarm of birds or fish. In the process of torque distribution, a particle is a point on a search space of TNE map. As shown in Figure 2a,b, three swarms of particles were initially distributed at random on the three TNE maps of the traction motors, and each swarm had N particles. Each particle had a position state, which was defined as the torque where the particle was located. The pedal command given by the EV driver was their common target at a specific vehicle speed.
Each particle at its initial position in the search space determined its best direction, and all the particles approached their common target with the minimal global effort.
During the PSO process, the N particles were renewed through J generations and reached the target of minimal energy consumption in the end. After the least energy consumption converged in each generation, the particles were updated according to their own best and swarm best solutions as where the sub-index i stands for the i th particle; j stands for the generation number; c L1 and c L2 are learning factors; G represents the swarm's best known solution of all the particles; P i represents its own best known solution of the i th particle; rand 1 and rand 2 are random values between 0 and 1; and w is the inertia weight. At the final generation, the output torques of the rear left and right motors were calculated: Details of the original work of the real-time torque-distribution strategy by PSO for a pure EV with three traction motors were described in [16].

Torque Distribution Strategy: Priority Torque Ratio in Front and Rear Motors
For comparison, a priority torque ratio (PTR), P r , can be assigned from 0 to 1 to the front motor with respect to the rear motors. For example, P r = 1 means that the front motor takes the first priority for delivering the torque within its feasible range for driving the EV, while P r = 0 expresses that the Appl. Sci. 2020, 10, 2653 8 of 18 rear motors have the top priority for delivering the torque over the front motor. Therefore, the torque given by the front motor is calculated as where T m f ,max and T m f ,min are, respectively, the upper and lower limits of the front motor at a certain speed. If T m f is calculated higher than or equal to the total torque T t required for acceleration, the front motor will provide T t as demanded. If T m f is calculated lower than T t even though P r = 1, the rear motors must be responsible for delivering the rest portion of torque according to (21) and (22).

Battery Energy Consumption
The energy consumption and balance of the three battery packs were investigated when the torque distribution was executed by the PSO or PTR strategy. In this study, the urban driving cycle (UDC) of the New European driving cycle (NEDC) was used because of the limited motor speeds and battery voltage. Table 2 presents two simulation results for the energy consumption for the three battery packs after the EV drove (1) on a straight road and (2) clockwise on a circular path with a 100 m radius using the proposed PSO and PTR torque distribution strategies. For the cases of P r = 1 and 0.75, the front traction motor had the highest priority by delivering power over the rear motors, and the rear battery restored more power from regenerative braking than that consumed for driving. This was an example of negative battery energy consumption. For the three battery packs, energy consumption was found to be unbalanced, thereby causing an unbalanced SOC. If any of the three battery packs was depleted, the torque distribution strategy would fail. This could cause a serious deterioration in vehicle maneuverability and stability. Figure 4 shows the flowchart for the torque and battery distribution (TBD) strategy. Before the EV started from rest, the SOC of the front, rear right, and rear left battery packs was respectively measured as SOC F , SOC L , and SOC R . The SOC gap and SOC ratio were calculated:
SOC gap (SOC g ): The SOC gap was defined as the difference between the SOC of the front battery pack and the lower SOC of the two rear battery packs. It was negative when the front battery had less power remaining than the rear batteries, and it was positive when the front battery had more power than the rear batteries. In applications, a default value SOC g < −κ% (0 < κ < 2) can be assigned to determine the torque distribution mode (Figure 4).

2.
SOC ratio (SOC r ): The SOC ratio was defined as the ratio of energy consumption in terms of the SOC between the front battery pack and the two rear battery packs. On the basis of the simulation of straight road driving under the PSO strategy, as indicated in Table 2, the SOC r converged to 2.94 for the long-term operation of the UDCs. In applications, a default value, ρ (< 2.94), was assigned to determine the torque distribution mode under the TBD strategy.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 9 of 18 Through the use of SOCg, the state of balance (SOB) of the three battery packs was investigated. The three states of battery balance are described: 1. State 1: The SOC of the front battery pack was higher than that of any of the rear battery packs. 2. State 2: The SOC of the front battery pack was equal to that of any of the two rear battery packs. 3. State 3: The SOC of the front battery pack was lower than that of the rear battery pack.
Three torque distribution modes were then determined by the SOC gap (SOCg): 1. Mode 1: was proposed when the front motor had the first priority for delivering the torque under the condition of SOCg and SOCr .
At Mode 1, the SOC of the front battery pack was much higher than that of the rear battery packs. When the SOCr was larger than 2.94, the SOCg could increase continuously, and the battery balance could worsen even though the PSO strategy was executed. It was better for the front battery pack to reach a balance between the front and rear batteries.
2. Mode 2: The PSO strategy was prescribed when there was not much difference in the SOC of the three battery packs under the following conditions: SOCg or (SOCg and SOCr ) or ( SOCg − % and |SOCg| was increasing).
Because the PSO strategy is superior to the PTR strategy for distributing the torque among the traction motors, it should be used whenever the difference in the SOC of the battery packs is negligible. For example, SOCg = 0 is an ideal case in which all of the batteries are in balance, and the PSO strategy can save more energy than the other PTR strategies when distributing the torque. With SOCg and SOCr , the amount of power stored by the front and rear batteries was similar and sufficient. It was also an appropriate situation for torque distribution under the PSO strategy.
When SOCg %, there was not much difference in the energy storage of the battery packs, and it was still safe to execute PSO even though the SOC gap was increasing.

Mode 3:
was proposed when the rear motors took top priority for delivering more the Through the use of SOC g , the state of balance (SOB) of the three battery packs was investigated. The three states of battery balance are described: 1.
State 1: The SOC of the front battery pack was higher than that of any of the rear battery packs.

2.
State 2: The SOC of the front battery pack was equal to that of any of the two rear battery packs. 3.
State 3: The SOC of the front battery pack was lower than that of the rear battery pack.
Three torque distribution modes were then determined by the SOC gap (SOC g ): 1. Mode 1: P r = 1 was proposed when the front motor had the first priority for delivering the torque under the condition of SOC g > 0 and SOC r ≥ ρ.
At Mode 1, the SOC of the front battery pack was much higher than that of the rear battery packs. When the SOC r was larger than 2.94, the SOC g could increase continuously, and the battery balance could worsen even though the PSO strategy was executed. It was better for the front battery pack to reach a balance between the front and rear batteries.

2.
Mode 2: The PSO strategy was prescribed when there was not much difference in the SOC of the three battery packs under the following conditions: SOC g = 0 or (SOC g > 0 and SOC r < ρ) or (0 > SOC g ≥ −κ% and |SOC g | was increasing).
Because the PSO strategy is superior to the PTR strategy for distributing the torque among the traction motors, it should be used whenever the difference in the SOC of the battery packs is negligible. For example, SOC g = 0 is an ideal case in which all of the batteries are in balance, and the PSO strategy can save more energy than the other PTR strategies when distributing the torque. With SOC g > 0 and SOC r < ρ, the amount of power stored by the front and rear batteries was similar and sufficient. It was also an appropriate situation for torque distribution under the PSO strategy.
When 0 >SOC g ≥ −κ%, there was not much difference in the energy storage of the battery packs, and it was still safe to execute PSO even though the SOC gap was increasing.

3.
Mode 3: P r = 0 was proposed when the rear motors took top priority for delivering more the torque than the front motor under the following conditions: SOC g < −κ% or (0 >SOC g > −κ% and |SOC g | was decreasing).
At Mode 3, the SOC of the rear battery packs was much higher than that of the front battery. The consumption in the rear batteries had top priority so that the balance in the batteries could be restored. After the SOC g was restored within [0, −κ%], Mode 3 remained in operation because the SOC gap continued to decrease until the best balance state 2 was achieved. This avoided frequent shifts between Modes 2 and 3.

Model-in-the-Loop Simulations
In Figure 5, the TBD strategy was simulated on a model-in-the-loop (MIL) platform. MATLAB Simulink was applied to model the TBD strategy, battery model, driver model, slip ratio control (SRC), and direct yaw moment control (DYC), while CarSim provided the vehicle dynamics. The driver model simulated human driver behavior by a proportional-integral (PI) controller. The SRC was responsible for stabilizing vehicle motion through the tractive control system (TCS) during acceleration and the anti-lock brake system (ABS) during deceleration. The TBD strategy was performed after vehicle safety was confirmed. Because of the limited motor speeds and battery voltages, only the UDC part of NEDC was used for MIL simulations.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 10 of 18 restored. After the SOCg was restored within [0, %], Mode 3 remained in operation because the SOC gap continued to decrease until the best balance state 2 was achieved. This avoided frequent shifts between Modes 2 and 3.

Model-in-the-Loop Simulations
In Figure 5, the TBD strategy was simulated on a model-in-the-loop (MIL) platform. MATLAB Simulink was applied to model the TBD strategy, battery model, driver model, slip ratio control (SRC), and direct yaw moment control (DYC), while CarSim provided the vehicle dynamics. The driver model simulated human driver behavior by a proportional-integral (PI) controller. The SRC was responsible for stabilizing vehicle motion through the tractive control system (TCS) during acceleration and the anti-lock brake system (ABS) during deceleration. The TBD strategy was performed after vehicle safety was confirmed. Because of the limited motor speeds and battery voltages, only the UDC part of NEDC was used for MIL simulations.
Both the straight and circular road tests were simulated. Figure 6a shows the SOC records of the three battery packs during the TBD process for the EV driving a clockwise cornering on a circular path of radius 100 m for 3.5 UDCs of the NEDC. The indices ρ and were assigned at 1.0056 and 1, respectively.
At the beginning, the SOC of the front battery (91%) was higher than that of both rear batteries (90%). In addition, SOCr , the torque distribution of Mode 1, was executed until the SOCr reduced to 1.0056 at approximately 170 s, where the PSO of Mode 2 was executed for torque distribution.
The torque distribution mode shifted from Mode 2 to Mode 3 when the SOCg was less than 1% at approximately 1375 s. The torque distribution mode shifted back to Mode 2 about 1575 s when the SOC gap was reduced. The SOC of the front and rear batteries remained within a SOB by shifting the torque distribution mode during the driving cycle. The torque distribution histories from three traction motors are presented in Figure 6b.  Both the straight and circular road tests were simulated. Figure 6a shows the SOC records of the three battery packs during the TBD process for the EV driving a clockwise cornering on a circular path of radius 100 m for 3.5 UDCs of the NEDC. The indices ρ and κ were assigned at 1.0056 and 1, respectively.
At the beginning, the SOC of the front battery (91%) was higher than that of both rear batteries (90%). In addition, SOC r > 1.0056, the torque distribution of Mode 1, was executed until the SOC r reduced to 1.0056 at approximately 170 s, where the PSO of Mode 2 was executed for torque distribution.
The torque distribution mode shifted from Mode 2 to Mode 3 when the SOC g was less than −1% at approximately 1375 s. The torque distribution mode shifted back to Mode 2 about 1575 s when the SOC gap was reduced. The SOC of the front and rear batteries remained within a SOB by shifting the torque distribution mode during the driving cycle. The torque distribution histories from three traction motors are presented in Figure 6b. It was also interesting to compare the differences of the energy economy of the proposed TBD strategy by having the battery energy storage in balance and using other torque distribution strategies without balancing the batteries. In these simulations, the initial SOC of each of the three battery packs was 90%, and their corresponding amounts of energy are 10.45, 5.2, and 5.2 kWh for the front, rear right, and rear left battery packs, respectively. The EV stopped when any of the batteries was depleted. It was found that both the rear batteries were exhausted soon for , 0.25, and 0.5 when the rear motors took higher priority for delivering more torque than the front motor; while the front battery was depleted soon for and 1 when the front motor had the first priority for delivering the torque to the EV.
Energy consumption efficiency was defined as the ratio between the total energy consumption and the initial battery energy. The energy consumption rate was defined as the total energy consumption per travel distance. Table 3 shows the energy consumption results for both the circular and straight road simulations. For clockwise cornering on a circular path of radius 100 m during the UDCs, the proposed TBD It was also interesting to compare the differences of the energy economy of the proposed TBD strategy by having the battery energy storage in balance and using other torque distribution strategies without balancing the batteries. In these simulations, the initial SOC of each of the three battery packs was 90%, and their corresponding amounts of energy are 10.45, 5.2, and 5.2 kWh for the front, rear right, and rear left battery packs, respectively. The EV stopped when any of the batteries was depleted. It was found that both the rear batteries were exhausted soon for P r = 0, 0.25, and 0.5 when the rear motors took higher priority for delivering more torque than the front motor; while the front battery was depleted soon for P r = 0.75 and 1 when the front motor had the first priority for delivering the torque to the EV.
Energy consumption efficiency was defined as the ratio between the total energy consumption and the initial battery energy. The energy consumption rate was defined as the total energy consumption per travel distance. Table 3 shows the energy consumption results for both the circular and straight road simulations. For clockwise cornering on a circular path of radius 100 m during the UDCs, the proposed TBD strategy had a travel distance: 142.6 km, which was 7.7% higher than the PSO strategy without the balancing of the SOC of the battery packs. The torque distribution strategy under PSO without battery balancing had a better energy consumption rate at 104.5 Wh/km than the TBD strategy, but the front battery was delpeted after 132.4 km, and the energy consumption efficiency was 73.6%. The rear battery packs had 26.4% energy remaining when the EV stopped. On the straight road, the torque distribution strategy under PSO without battery balancing exhibited the best energy consumption rate: 75.21 Wh/km. However, the front battery was depleted after 167.6 km, and energy consumption efficiency was 67.5%. Thus, the rear battery packs had only 32.5% energy remaining when the EV stopped. The proposed TBD strategy of partly using torque distribution by PSO had the highest driving range: 214.4 km, i.e., approximately 248 UDCs. This was attributed to the SOC of the three battery packs being kept in balance to avoid unexpected battery depletion. Therefore, approximately 28% more driving range was extended by the TBD strategy than by the PSO strategy without battery balancing.
It was also found in the simulation that the battery energy was fully utilized by the TBD strategy for the EV on a straight road of UDC. The 99.9% energy consumption efficiency for the TBD strategy was calculated in Table 4. The energy of three battery packs can be effectively distributed and utilized to extend the driving range, when the SOC gap of the three battery packs remains within a prescribed limit during vehicle operation. Table 4. Energy consumption efficiency for the EV on a straight road of UDC by the TBD strategy.

Battery Front Rear Left Rear Right
Energy capacity (kWh) @ 100% SOC 10.45 5. Other strategies with a PTR, P r from 0 to 1, presented lower travel distances and less energy consumption efficiency than observed for the proposed TBD strategy. Figure 7 presents the architecture of hardware-in-the-loop experiment. A Mitsubishi Colt-Plus was retrofitted with a 15-kW radial-flux PMSM and a 144-V battery pack for front wheels and two 7-kW axial-flux PMSMs and two 72-V battery packs for rear wheels. This EV was set up on a Horiba MAHA-AIP ECDM-48 emission chassis dynamometer. The maximum test speed was 200 km/h, the maximum power absorbing was 150 kW at 100 km/h, the maximum tractive force was 5400 N for light duty and 6750 N for heavy duty, and the maximum vehicle inertia simulation was 4540 kg for 4-wheel drives.

Hardware-in-the-Loop Simulations
Appl. Sci. 2020, 10, x FOR PEER REVIEW 13 of 18 Figure 7 presents the architecture of hardware-in-the-loop experiment. A Mitsubishi Colt-Plus was retrofitted with a 15-kW radial-flux PMSM and a 144-V battery pack for front wheels and two 7-kW axial-flux PMSMs and two 72-V battery packs for rear wheels. This EV was set up on a Horiba MAHA-AIP ECDM-48 emission chassis dynamometer. The maximum test speed was 200 km/h, the maximum power absorbing was 150 kW at 100 km/h, the maximum tractive force was 5400 N for light duty and 6750 N for heavy duty, and the maximum vehicle inertia simulation was 4540 kg for 4-wheel drives. In experiments, the battery SOC information and throttle command were received by a real-time rapid control prototyping unit dSPACE MicroAutoBox II, in which the TBD strategy was built to determine the torque command for each motor. This MicroAutoBox took a role of vehicle control unit with an 800-MHz processor, 18-MB main memory, 16-MB flash memory, and dual CAN interfaces.

Hardware-in-the-Loop Simulations
In the HIL experiment, the vehicle followed the complete NEDC but the maximum speed at EUDC was restricted at 40 km/h in the HIL experiment. The initial SOCs of the front, rear left, and rear right batteries were 94.8%, 93.6%, and 95.5%, respectively. In order to make the experiment efficient, driving modes were shifted at SOCg % and SOCr , according to the flowchart of the TBD strategy in Figure 4. Figure 8 illustrates the time history of mode, state, SOCr, and SOCg and the corresponding torque distribution histories of the front, rear left, and rear right motors during the TBD process in the hardware-in-the-loop experiment. During the first 95 s, the SOC of the front battery pack was higher than the SOC of anyone of the rear battery packs. Therefore, State 1 (SOCg ) was identified and SOCr was larger than 1.006, the front motor took the first priority of delivering torque and Mode 1 ( ) was executed. Between 95 and 200 s, it was still at State 1 (SOCg ) but SOCr was less than 1.006, Mode 2 was executed with the PSO strategy.
Between 200 and 300 s, State 2 (SOCg ) was identified, i.e., the SOC of the front battery pack was equal to any one of the two rear battery packs. Mode 2 with the PSO strategy was working. Between 300 and 485 s, the SOC of the front battery pack was lower than the SOC of the rear battery packs, and State 3 (SOCg ) was identified. Because SOCg −1% and |SOCg| was increasing, Mode 2 remained.
When SOCg was less than 1% at State 3 between 485 and 585 s, the rear ba ery packs and motors started to take their top priority for delivering more power than the front battery and motor, and Mode 3 ( ) was executed. Between 585 and 680 s, the SOC gap was between 1% and 0%, but |SOCg| was decreasing, Mode 3 remained. In experiments, the battery SOC information and throttle command were received by a real-time rapid control prototyping unit dSPACE MicroAutoBox II, in which the TBD strategy was built to determine the torque command for each motor. This MicroAutoBox took a role of vehicle control unit with an 800-MHz processor, 18-MB main memory, 16-MB flash memory, and dual CAN interfaces.
In the HIL experiment, the vehicle followed the complete NEDC but the maximum speed at EUDC was restricted at 40 km/h in the HIL experiment. The initial SOCs of the front, rear left, and rear right batteries were 94.8%, 93.6%, and 95.5%, respectively. In order to make the experiment efficient, driving modes were shifted at SOC g = −κ% = −1% and SOC r = ρ = 1.006, according to the flowchart of the TBD strategy in Figure 4. Figure 8 illustrates the time history of mode, state, SOC r , and SOC g and the corresponding torque distribution histories of the front, rear left, and rear right motors during the TBD process in the hardware-in-the-loop experiment. During the first 95 s, the SOC of the front battery pack was higher than the SOC of anyone of the rear battery packs. Therefore, State 1 (SOC g > 0) was identified and SOC r was larger than 1.006, the front motor took the first priority of delivering torque and Mode 1 (P r = 1) was executed. Between 95 and 200 s, it was still at State 1 (SOC g > 0) but SOC r was less than 1.006, Mode 2 was executed with the PSO strategy.
Between 200 and 300 s, State 2 (SOC g = 0) was identified, i.e., the SOC of the front battery pack was equal to any one of the two rear battery packs. Mode 2 with the PSO strategy was working. Between 300 and 485 s, the SOC of the front battery pack was lower than the SOC of the rear battery packs, and State 3 (SOC g < 0) was identified. Because 0 > SOC g ≥ −1% and |SOC g | was increasing, Mode 2 remained.
When SOC g was less than −1% at State 3 between 485 and 585 s, the rear battery packs and motors started to take their top priority for delivering more power than the front battery and motor, and Mode 3 (P r = 0) was executed. Between 585 and 680 s, the SOC gap was between −1% and 0%, but |SOC g | was decreasing, Mode 3 remained. The corresponding SOCs of the three battery packs is shown in Figure 9. In the first 500 s, the front battery provided all power to the vehicle. After 500 s, the rear and front batteries powered the vehicle alternatively, so that the SOC gap would always remain in 1% as expected.  The corresponding SOCs of the three battery packs is shown in Figure 9. In the first 500 s, the front battery provided all power to the vehicle. After 500 s, the rear and front batteries powered the vehicle alternatively, so that the SOC gap would always remain in 1% as expected. The corresponding SOCs of the three battery packs is shown in Figure 9. In the first 500 s, the front battery provided all power to the vehicle. After 500 s, the rear and front batteries powered the vehicle alternatively, so that the SOC gap would always remain in 1% as expected.

Road Tests
The road test was executed on the Industrial Technology Research Institute (ITRI) campus. The total driving distance was about 5.4 km, during which the road slope varied and the maximum speed was 30 km/h. The driving curve is shown in Figure 10. The initial SOCs of the front, rear left, and rear right batteries were 72.8%, 70%, and 71.2%, respectively. Driving modes were shifted at SOC g = −κ% = −1% and SOC r = ρ = 1.02, according to the flowchart of the TBD strategy in Figure 4. Figure 11 illustrates the time history of mode, state, SOC r , SOC g , and the corresponding torque distributions.

Road Tests
The road test was executed on the Industrial Technology Research Institute (ITRI) campus. The total driving distance was about 5.4 km, during which the road slope varied and the maximum speed was 30 km/h. The driving curve is shown in Figure 10. The initial SOCs of the front, rear left, and rear right batteries were 72.8%, 70%, and 71.2%, respectively. Driving modes were shifted at SOCg and SOCr , according to the flowchart of the TBD strategy in Figure 4. Figure  11 illustrates the time history of mode, state, SOCr, SOCg, and the corresponding torque distributions.

Road Tests
The road test was executed on the Industrial Technology Research Institute (ITRI) campus. The total driving distance was about 5.4 km, during which the road slope varied and the maximum speed was 30 km/h. The driving curve is shown in Figure 10. The initial SOCs of the front, rear left, and rear right batteries were 72.8%, 70%, and 71.2%, respectively. Driving modes were shifted at SOCg and SOCr , according to the flowchart of the TBD strategy in Figure 4. Figure  11 illustrates the time history of mode, state, SOCr, SOCg, and the corresponding torque distributions.  Similar to the result of HIL test, the front battery pack took the first priority to supply power to drive the EV at Mode 1 (P r = 1) during the first 180 s. Because the EV moved upslope approximately at 150 s, two rear motors delivered extra torque. Between 180 and 480 s, Mode 2 with the PSO strategy worked when the SOC of the front battery was lower than that of the rear battery packs. The front motor still provided the major torque until the SOC gap was less than −1% (SOC g < −1%) while Mode 3 (P r = 0) was executed. Then, Modes 2 and 3 shifted alternatively. The SOC difference of the three battery packs was finally kept within 1%, as shown in Figure 12.
Similar to the result of HIL test, the front battery pack took the first priority to supply power to drive the EV at Mode 1 ( ) during the first 180 s. Because the EV moved upslope approximately at 150 s, two rear motors delivered extra torque. Between 180 and 480 s, Mode 2 with the PSO strategy worked when the SOC of the front battery was lower than that of the rear battery packs. The front motor still provided the major torque until the SOC gap was less than (SOCg ) while Mode 3 ( ) was executed. Then, Modes 2 and 3 shifted alternatively. The SOC difference of the three battery packs was finally kept within 1%, as shown in Figure 12.
The two strategies with a PTR, and , were also executed in the road test. The total energy consumption and the consumption of each battery pack are given in Table 5. It shows that the energy consumption rate of the TBD strategy was lower than those of and . It means that with the same battery energy, the EV has a longer driving range by the proposed TBD strategy than that by other PTR strategies. For example, the TBD strategy extended 11% more driving range than the PTR strategy when and the TBD strategy extended 23.5% more driving range than the PTR strategy when .

Conclusions
A novel TBD strategy has been proposed for the EV with three independent traction motors and battery packs. Upon acceleration of the EV, the demanded torque was provided by all three traction motors together at their highest efficiency under the PSO strategy for saving battery energy. Simultaneously, the SOC of the three battery packs had to be kept in balance to avoid any unexpected battery depletion and to improve the EV's driving range. Thus, a combination of PSO and the PTR for torque distribution strategy was applied to compromise between energy saving and energy balance. On the basis of the model-in-the-loop simulations, the proposed TBD strategy shows better travel distance and the higher energy consumption efficiency than a pure PSO method or PTR The two strategies with a PTR, P r = 0 and P r = 1, were also executed in the road test. The total energy consumption and the consumption of each battery pack are given in Table 5. It shows that the energy consumption rate of the TBD strategy was lower than those of P r = 0 and P r = 1. It means that with the same battery energy, the EV has a longer driving range by the proposed TBD strategy than that by other PTR strategies. For example, the TBD strategy extended 11% more driving range than the PTR strategy when P r = 0, and the TBD strategy extended 23.5% more driving range than the PTR strategy when P r = 1.

Conclusions
A novel TBD strategy has been proposed for the EV with three independent traction motors and battery packs. Upon acceleration of the EV, the demanded torque was provided by all three traction motors together at their highest efficiency under the PSO strategy for saving battery energy. Simultaneously, the SOC of the three battery packs had to be kept in balance to avoid any unexpected battery depletion and to improve the EV's driving range. Thus, a combination of PSO and the PTR for torque distribution strategy was applied to compromise between energy saving and energy balance. On the basis of the model-in-the-loop simulations, the proposed TBD strategy shows better travel distance and the higher energy consumption efficiency than a pure PSO method or PTR strategies. Similar results were also proved on a real vehicle for hardware-in-the-loop experiments on dynamometer and road tests. The road test proved that the TBD strategy extended 11% and 23.5% more driving range than other PTR strategies, when P r = 0 and P r = 1, respectively. The proposed TBD strategy is promising for extending the driving range of an EV with multiple traction motors and battery packs with an improved energy consumption efficiency.