# I. INTRODUCTION

Environmental problems such as global warming, the depletion of fossil fuels, and air pollution are becoming serious. Therefore, electric vehicles (EVs) have attracted a great deal of interest as zero-emission vehicles. EVs have the following advantages [1]:

Some applications in the field of electrical drives require the use of several electric machines and many static converters that have an important place among electromechanical systems. These systems are called multi-machine multi-converter systems (MMSs) [2]. MMSs can be considered as extensions of classical drives. In numerous applications, one motor is controlled by one converter known as a single-machine single-converter system. However, in high-power applications such as traction systems, conveyer lines, and steel processing, two or more machines are fed by one converter. This topology results in a light, highly compact, and low-cost system [3]. This system is called a multi-machine single-converter system. The control of multi-machine single-converter systems is the subject of this study. Several methods have been proposed for the control of multi-machine single-inverter systems. In this case, a master–slave based direct torque control (DTC) strategy is developed to generate the electric differential.

Electric differential is another technology that eliminates traditional mechanical differentials [4]-[6]. This system can reduce the weight and friction loss in a driving train by removing the mechanical differential components and consequently improving the overall reliability and efficiency during power transmission. The working principle of an electric differential system is similar to that of a traditional mechanical differential system. When an EV is cornering or making a turn, the outer wheels travel longer paths than the inner wheels. This implies that the rotational speed of the outer wheels will be faster than that of the inner wheels. The resulting speed differences between the inner and outer wheels are automatically tuned by an electric differential system.

An advantage of EVs is the possible increase in safety because of the DTC and speed during cornering and risky maneuvers. Electric differential is characterized by the following features: 1) no mechanical link exists between the driven wheels; 2) traction power is separately applied to each wheel by a speed controller, and 3) the speed controller applies minimal power to the inner wheel.

EVs equipped with in-wheel motors are envisioned to be future personal EVs on the basis of their energy efficiency and motion control. Over the past few years, research on motion control, including traction control or yaw stability control, has been conducted to utilize independently driven in-wheel motors [1], [7], [8]. In-wheel electrical motor control [9], [11] is a novel technology that provides fast and accurate torque actuation.

Permanent magnet synchronous motors (PMSMs) have been extensively analyzed as feasible candidates for variable speed EV traction applications [12]. PMSMs exhibit high efficiency when operating at constant speeds in the constant torque region because of their low (copper) losses. Vector control techniques are preferred to improve the dynamic performance of induction motor drives for EV propulsion. An innovative control method called DTC has gained attention for electric propulsion systems [13], [15] because, unlike vector control, it can produce fast torque control of PMSMs and does not require heavy on-line computation.

Sensorless control is an attractive feature in DTC for PMSMs without a speed control loop, in which no position or speed signal is needed for the motor and only an initial position is needed for the PMSMs. However, in high-performance applications, speed control is usually required. In this case, the rotor speed should be detected for feedback. As discussed earlier, the use of speed sensors outweighs the advantages of DTC systems. Many researchers have paid considerable attention to the sensorless method of DTC for PMSMs to eliminate the cost of the speed sensor and to improve the reliability of the driving system.

This paper is organized as follows: To improve the dynamic performance of permanent magnet synchronous in-wheel motors, a DTC based on a speed model reference adaptive system (MRAS) observer is proposed. Simulation results are presented to validate the effectiveness of the proposed sensorless control. Section II introduce a three degrees-of-freedom (3-DOF) yaw plan model to describe the lateral motion in EVs. Section III discusses the design of an electric differential for EVs with four independent in-wheel motors based on the new DTC algorithm. Finally, a series of Matlab/Simulink simulation tests is conducted to evaluate the performance of the proposed sensorless master–slave DTC system applied to electric differentials to improve the stability of in-wheel motor-driven EVs.

# II. SENSORLESS DTC WITH A MRAS OBSERVER

A block diagram of a sensorless DTC-PMSM that employs a MRAS as speed estimator is shown in Fig. 1. MRAS based estimators provide the desired state from two different models: the reference model and the adjustable model [16] [20]. The error between the two models is used to estimate the unknown parameter (speed in this case). In the MRAS, only the adjustable model depends on the unknown parameter [21], and the reference model is not dependent on speed. The error signal is fed into an adaptation mechanism to provide the estimated quantity, which is used to tune the adjustable model. A sensorless control algorithm (Fig. 2) is employed in this study. A PMSM is considered as the reference model, and the stator current equations are considered as the adjustable model. The estimated rotor speed is used to tune the adjustable model based on the current equations of the motor. The speed error is continuously monitored to ensure negative feedback and, hence, overall system stability. The current equations of the PMSM are given as:

The above current equations can be rewritten as:

where:

**Fig. 1.**diagram of sensorless DTC-PMSM drive with MRAS based speed estimator.

**Fig. 2.**Structure of the MRAS Block.

Eq. (4) has the speed as a variable which will be used for the adjustable model. A PMSM is used as the reference model and provides id and iq . In this algorithm, the

estimated value of the rotor speed is calculated as:

The performance of the MRAS observer is based on the controller used in the adaptation mechanism. In this case, a proportional-integral controller is used to tune the adaptive model based on the reference model. The controller forces the error to be zero. Fig. 3 shows the results of a test in which the sensorless DTC drive is subjected to a stepwise speed changes from 100 rad/s to zero speed in a series of 20 rad/s steps and then back to 100 rad/sec at no load. Fig. 3 indicates that the actual speed and the estimated speed follow the reference speed. Fig. 4 shows a waveform of the speed error between the estimated speed and the actual speed. Fig. 5 shows the sensorless performance during a speed change from 100 rad/s to zero speed at no load. Stable operation and no-oscillating-speed performance are obtained at around the zero speed.

**Fig. 3.**Sensorless performance for the staircase test.

**Fig. 4.**Error between estimated and actual speeds.

**Fig. 5.**Estimated and actual speeds at low speeds.

To test the load torque disturbance capability of the sensorless DTC-PMSM drive at low speed, different levels of load torque are applied at 20 rad/s. Fig. 11 shows the electromagnetic torque generated by the PMSM as a step change in the load torque is applied to the motor. The sensorless performance at low and zero speeds is also good.

**Fig. 6.**Error between estimated and actual speeds at low speed.

**Fig. 7.**Sensorless performance for speed change (a) from 60 rad/s to 40 rad/s, and (b) from 40 rad/s to 60 rad/s.

**Fig. 8.**Error between estimated and actual speeds at low speed.

**Fig. 9.**Sensorless performance for load torque application at 20 rad/s (50% rated load).

**Fig. 10.**Error between estimated and actual speeds for load torque application.

**Fig. 11.**Electromagnetic torque generated by PMSM.

# III. VEHICLE DYNAMICS

The structure of the EV in this study is composed of four in-wheel motors (i.e., PMSMs) mounted in each wheel, as shown in Fig. 12 [28]. Therefore, the wheel torque for each wheel can be controlled completely and independently for vehicle motion control. Regenerative braking is also possible.

**Fig. 12.**Configuration of an EV with four independent in-wheel motors.

In this section, a 3-DOF yaw plane model is introduced to describe the lateral motion in EVs. The yaw plane representation is shown in Fig. 14.

**Fig. 13.**In-wheel motors.

**Fig. 14.**3-DOF yaw plane vehicle model.

The vehicle dynamics are described by the longitudinal velocity, lateral velocity and yaw rate as follows [6], [29], [31] :

and:

With a linear tire model, the front and rear cornering forces can be expressed as a product of the cornering stiffness (Cf, Cr) and the sideslip angle (αf , αr).

The sideslip angles of the wheels can be easily expressed in terms of the longitudinal, lateral, and angular velocities, as well as the steering angle δ. Explicit expressions of the sideslip angles for the front and rear axles are represented by Eq. (10).

The longitudinal slip can be defined for the four wheels as:

where:

i = 1,2,3, and 4 correspond to the front left, front right, rear left, and rear right ( = lf , fr, rl, rr ) wheels, respectively; Rω is the wheel radius; ωi is the angular velocity of the in-wheel motor; and v is the linear speed at which the contact zone moves on the ground.

The interrelationships between the slip ratio λ and the traction coefficient μ can be described by various formulas. In this study, the widely adopted magic formula [31] [33] is applied to describe the relationship between the slip and traction forces and to build a vehicle model for the following simulations, as shown in Eq. (11).

The coefficient sets of c1 , c2 , c3 and c4 are defined in [31].

The longitudinal forces for the four in-wheel motors can be calculated using the following equation:

The drive system model can be described by the following mechanical equations:

where Tri is the resistive torque; Nf, Nr are the front and rear normal forces and have the following expressions:

# IV. KINEMATIC AND DYNAMIC MODELS FOR THE ELECTRIC DIFFERENTIAL

The electric differential operation supposes the solution of two technological problems: 1) the wheel synchronization and 2) the computation of the relative wheel speed as a function of the turn angle. The relationship between the inner two and outer two wheel velocities in a corner are usually described through the Ackerman steering principle. This principle computes the relative speed differences between the wheels by using the data of the turn (steering) angle. A new DTC algorithm is used to control a bi-machine traction system [22]- [27].

As shown in Fig. 15, the linear velocities of the four wheels are represented by Vfl , Vfr , Vrl , and Vrr . The steering angles of the front wheels are shown as δfl and δfr . The distances from the front wheels to the rear wheels of the EV [centre of gravity (CG)] are shown as If and Ir . Half of the vehicle track is represented by Is .

**Fig. 15.**Kinematic model of an EV in a turn maneuver.

When the EV is driven by all four wheels in a cornering maneuver, the inner two wheels and the outer two wheels will have different rotational speeds. On the basis of the kinematic model shown in Fig. 15, the linear velocities of the four wheels are described as follows [34], [35]:

The steering angles for the front right and front left wheels are assumed to be the same, namely, δfl = δfr = δ .

All of the tires are assumed to have the same effective radius Rω . The reference rotational speed of each wheel during a turn maneuver is described as follows:

The difference between the angular speeds of the wheel drives is expressed by Eq. (17). The signal of the steering angle indicates the curve direction [Eq. (18)].

When a vehicle arrives at the beginning of a curve, the driver applies a steering angle to the wheels. However, the electronic differential immediately acts on the two motors, reduces the speed of the inner wheel, and increases the speed of the outer wheels. The driving wheel angular speeds are:

# V. SIMULATION RESULTS

The objective of the stability control is to improve the vehicle steadiness and transient response properties, thereby enhancing the vehicle handling performance and maintaining stability with regard to maneuvers. The desired vehicle responses are defined on the basis of the concerning intention of the driver (e.g., the steering command and vehicle speed imposed by the driver). The vehicle speed starts from zero to a chosen reference speed of 30 km/h. At this operating point, a cycle variation is imposed on the vehicle chassis by the steering command of the driver, as illustrated in Fig. 17(a). The input hand-steering angle linearly decreases to - 70° within 1 s because the gear ratio from the hand-steering wheel to the front driving wheels is approximately 10 s, and the maximum steering angles for both of the front wheels (assumed to be the same) can reach - 7°. Fig. 17(a) shows the longitudinal velocity of the vehicle vx when the vehicle reaches a constant speed at 13 s. A delay in reaching this speed is allowed due to the acceleration of the vehicle mass. From Figs. 17(c) and (d), the lateral velocity vy and the yaw moment r depend on the steering command of the driver. These two speeds occur only during cornering and they vanish when the vehicle is traveling on a straight road.

**Fig. 16.**Block diagram of the proposed electric differential based on a new master–slave DTC technique.

**Fig. 17.**Traction system behavior in a turn maneuver.

When the EV is driven by all four wheels in a cornering maneuver, the inner two wheels and the outer two wheels will have different rotational speeds, as shown in Figs. 17(g), (h), (i). This turn is assumed to be carried out at a constant speed. Figs. 17(e) and (f) illustrate the linear speeds of the wheels. It can be seen that a slight decrease occurs when the driver applies a steering command. However, the driver makes a steering angle which will be the steering angle of the front wheels. Thus, the electric differential immediately acts on the four motors simultaneously. This in turn lowers the speed of the inner two wheels in comparison with the outer two wheels. At such times, the speeds of the wheels change values. This affects the electrical and mechanical quantities, as shown in Figs. 17(s), (t), (u). During the first steering manoeuvers, the outer motors (M1 and M2) rotate at higher speeds than the inner motors (M3 and M4). When the speeds of the inner wheels are decreasing in accordance with their new references, the torque tries to change its sign because of the sudden change in speed [Figs. 17(j) and (k)]. Hence, the motors work in a braking mode by developing a negative torque. This working phase can be exploited for energy recuperation to the battery. The motors (M1 and M2) during the second steering maneuvers are presented in Figs. 17(r) and (i). In this case, the vehicle is turning on the left side at a speed of 30 km/h [Fig. 17(e)]. Despite the fact that the driving wheels are following different paths, they turn in the same direction with different speeds. The left driving wheels turn with less speed than the right wheels. The behavior of these speeds is illustrated in Figs. 17(h) and (i). When the speeds of the right wheels increase, the torque motors associated with these wheels increase and try to reach the growth of the speeds. Figs. 17(o) and (p) illustrate the variation in the electromagnetic torque of the motors. In the transient state, the torque is limited to 145 N.m.

The variation in the phase currents for each motor is presented in Figs. 17(l)-(o). Two intervals for the vehicle dynamics exist. The first interval occurs when the vehicle starts moving. In this interval, each motor solicits a high current to attain the reference speed imposed by the driver. This current is due to the starting torque caused by the vehicle inertia. The second interval starts after 17 s. In this interval, each motor tries to develop an electromagnetic torque.

Fig. 17(p) illustrates the traction forces generated by the motors of the front axle (M1 and M3) and the motors of the rear axle (M2 and M4). High traction forces are provided by the motors to move the vehicle from a standstill. This phenomenon seems logical because these forces must overcome the overall forces resisting vehicle movement. A discrepancy exists between the traction forces during a movement through a turn. The motors located within the curvature of the turn produce lower traction forces than the motors that are outside the curvature of the turn. This scheme also applies to the load torques imposed on the motors, as shown in Figs. 17(q) and (r).

Fig. 17(p) shows the difference in traction forces of the motors during the completion of a turn. A difference is also observed for the speed of the motors [Fig. 17(g)] and the resistant torque of the motors [Fig. 17(q) and (r)]. When the speeds of the motors on the left side and right side diverge, the left wheels slow down, whereas the right wheels speed up to traverse the turn.

As shown in Fig. 17(s) and during the turn trajectory of the vehicle, the lateral forces applied to both the front and rear driving wheels show different values. The lateral side slip of the front and rear wheels are shown in Fig. 17(t). The behaviour of the lateral slip of the wheels is illustrated in Fig. 17(u). The longitudinal slips of the four wheels are maintained in the adhesive region. Fig. 17(v) shows that the adhesive coefficient of the outer driven wheels is larger than that of the inner wheels. Therefore, stability is maintained during vehicle turns. The simulation results validate the design of the electric differential under normal driving manoeuvres.

# VI. CONCLUSION

In this study, a novel multi-machine control based on an electric differential system for EV applications is proposed. A new DTC algorithm is used to control the bi-machine traction system based on a speed MRAS observer. Simulation results shown that a vehicle with the proposed multi-machine robust control, which is composed of four in-wheel PMSMs, can successfully follow the trajectories defined by the steering maneuvers of a driver.