# A Robust Control Method for Lateral Stability Control of In-Wheel Motored Electric Vehicle Based on Sideslip Angle Observer

Authors: Yaxiong Wang; Feng Kang; Taipeng Wang; Hongbin Ren
Journal: Shock and Vibration
 (2018)

---

## Abstract

In-wheel motored powertrain on electric vehicles has more potential in maneuverability and active safety control. This paper investigates the longitudinal and lateral integrated control through the active front steering and yaw moment control systems considering the saturation characteristics of tire forces. To obtain the vehicle sideslip angle of mass center, the virtual lateral tire force sensors are designed based on the unscented Kalman filtering (UKF). And the sideslip angle is estimated by using the dynamics-based approaches. Moreover, based on the estimated vehicle state information, an upper level control system by using robust control theory is proposed to specify a desired yaw moment and correction front steering angle to work on the electric vehicles. The robustness of proposed algorithm is also analyzed. The wheel torques are distributed optimally by the wheel torque distribution control algorithm. Numerical simulation is carried out in Matlab/Simulink-Carsim cosimulation environment to demonstrate the effectiveness of the designed robust control algorithm for lateral stability control of in-wheel motored vehicle.


---

## Body

## 1. Introduction

In recent years, active safety systems have been developed and commercialized aiming at improving vehicle high speed safety and making drivers more aware of the situation around them. The objects of these systems are to keep the vehicle stable under extreme driving conditions and follow the driver’s intention. All body motion control and disturbance forces as well as aerodynamic forces applied to the vehicle are generated in the contact patches between the tire and road surfaces. The maximum tire friction forces are determined by the road surface conditions and tire vertical forces. The summations of longitudinal and lateral tire forces cannot exceed the boundary of the tire friction ellipse. In-wheel motor EVs have many advantages as a platform for vehicle motion control in the viewpoint of vehicle stability control [1–4]. Motors are mounted inside each wheel and driven independently, and this brings more potential in the improvement of vehicle handling and stability performance because of its fast response with precise control and high efficiency. The traction and brake forces of each individual wheel can be independently controlled by means of controlling the corresponding electric drive motor [5]. Active safety systems can significantly prevent vehicle accident or relieve drivers’ workload. And with the in-wheel motor powertrain system, the advanced driving assistance system (ABS, braking assistance, electronic stability program (ESP), etc.) could be realized by software without additional mechanical or hydraulic components with greater reliability and increased safety [6]. In-wheel motored cars are capable of offering more passenger space than traditional engine-powered cars [7]. It is necessary to optimize the tire forces distribution and make the forces exerted by all tires efficiently, so as to keep them inside the tire friction ellipse. Additionally, with stringent vehicle emission regulation and fast development of motor and battery technologies, the in-wheel motor powertrain has become one of the most popular configurations in the area of alternative energy vehicles (AEVs) and are currently being developed by various car manufacturers [8–10]. AEVs have better mobility and fuel economy performances compared to the conventional internal combustion engine vehicles [11–13].Extensive research has been conducted by using different control methodologies for in-wheel motored vehicle control systems in recent years [14–16]. Zhao et al. [17] proposed the model predictive control allocation for yaw stability improvement of four-wheel drive electric vehicles in critical driving conditions. The proposed allocation strategy obtained a desired yaw moment to apply on the EVs by the upper level control system, while the control allocation is used to determine control inputs for four driving motors by commanding appropriate wheel slips. A correctional linear quadratic regulator combined the feedback and feed forward control algorithm is introduced by Li et al. to deduce the object of the stability yaw moment in order to guarantee the yaw rate and sideslip angle stability. Her et al. [18] presented an integrated chassis control of the differential braking, the front and rear traction torques, and the active roll moment for optimized tire force coordination to enhance the limit handling performance. Zhang and Wang [19] presented a robust gain-scheduling approach for vehicle lateral dynamics control through active front steering and direct yaw moment control. van der Sande et al. [20] proposed robust control analysis for a steer-by-wire vehicle considering the nonlinear characteristics of tire forces. Li studied the active front steering and direct yaw control coordination control problem based on model predictive control theory. Wu et al. [21] proposed an improved robust internal model control algorithm blending model tracking and internal model control for active steering system in order to reach high performance of yaw rate tracking with certain robustness.In order to improve the robust performance of vehicle lateral stability especially for high speed condition, a robust controller is designed in this paper, which is insensitive to external and internal disturbances and parameters variations, such as tire stiffness variation when tire forces are in saturation. Combining the sideslip angle observer, the robust controller is designed to generate the yaw moment to improve the longitudinal performances and lateral stability of the in-wheel motored vehicle.The paper is organized as follows: vehicle dynamic model and control reference model are presented in Section2, while in Section 3, the sideslip angle estimation is designed for the controller, and system controller is designed and analyzed based on robust control theory. Tire forces distribution and in-wheel motor control are illustrated in Sections 4 and 5, respectively; the efficiency of the proposed control is demonstrated by simulation results in Section 6. At the end of this paper, the conclusion and future work are given.


## 2. Vehicle Dynamic Modeling and Reference Model

The single-track model is widely used in the study of vehicle lateral dynamics. The single-track model is also called bicycle model, as shown in Figure1, which considers lateral and yaw motions under the assumption of negligible lateral load transfer, roll, and compliance steer while traveling on smooth road with constant speed. One of the main effects taken into account is the tire lateral slip phenomenon. However, the other main effects regarding vehicle behavior (such as axle deformation, axle kinematics, and suspension elastokinematics effect) are not modeled. The lateral tire force can be approximated as linear relationship with tire slip angle when vehicle lateral acceleration is less than 0.3 g. So it is fairly accurate for the linear bicycle model, and the lateral tire force generation falls into the linear range. As tire slip angle increases due to a large steering maneuver which results in large lateral acceleration, the model will lose the fidelity compared to those responses of actual vehicle [22]. In this paper, a single-track vehicle model considering the tire forces saturation is used for state estimation and stability control.Figure 1
Dynamic bicycle model.The bicycle model can be expressed as(1)v˙=Fyfcosδf+Fyrm-ru,r˙=aFyfcosδf+bFyr+ΔTIz, where v is lateral velocity; r is yaw rate; δf is steering angle input; ΔT is yaw control torque.(2)Fy=kiαii=f,r.With front tire stiffness, kf∈[kfmin,kfmax], and with rear tire stiffness, kr∈[krmin,krmax], where the subscripts f and r denote front axle and rear axle and α is tire lateral slip angle.The sideslip angle of mass center is(3)β=vu,ax=u˙-rv,ay=v˙+ru, where ax,ay are longitudinal and lateral acceleration of vehicle mass center, respectively.For simplicity, a neutral steering model (i.e.,bkr-akf=0) is chosen as reference model. The state space equations of linear bicycle model can be expressed as(4)X˙ref=ArefXref+BrefU,v˙r˙=-kf+krmu-10-a2kf+b2krIzuvr+krm0akfIz1IzδΔT.The transfer functions can be deduced as(5)Gδ→rsref=akfuIzus+a2kf+b2kr,GΔT→rsref=ΔTIzus+a2kf+b2kr.So, the reference yaw rate can be defined as(6)rref=Gδ→rsrefδs+GΔT→rsrefΔTs.Because the yaw rate is mainly affected by the steering angle input and the contribution ofΔT is relatively small, (6) can be further simplified as(7)rref=KGδ→rsrefδs,where K is correction coefficient which considers the effects from ΔT and the maximum lateral tire forces can be generated.


## 3. Robust Controller Design

Optimal control algorithms are widely used in implementation of the industry control, but it is not always tolerant to changes in the environment disturbances or the control system. A control system is robust if it is insensitive to differences between the actual system and the nominal model of the system which is used in the controller design [23]. These differences are considered as model uncertainty or perturbations. Application of the robust control theory is important to build the reliable and precise systems [24, 25]. A robust control algorithm is designed to improve the lateral stability for in-wheel EV combined with online sideslip angle observer, as shown in Figure 2.Figure 2
The diagram of the control logic.


### 3.1. Sideslip Angle Observer Based on Unscented Kalman Filter

From a practical point of view, sideslip angle is considered as a significant signal in assessment of the transient response of vehicle stability. In this section, the estimation method of sideslip angle of vehicle mass center is presented based on the virtual tire forces sensor, as shown in Figure3. The vehicle state information is estimated by using the UKF, which was proposed by Julier and Uhlman [26]. Detailed explanations of the UKF theory can be found in Ref [27]. The flowchart of the sideslip angle estimator is presented in Figure 4.Figure 3
The flowchart of the sideslip angle estimator.Figure 4
Lateral tire forces under different road friction coefficient.Based on (1)–(3), the vehicle dynamic system can be rewritten in the discrete-time nonlinear state transition equation:(8)xk+1=fxk,uk,wk,yk=hxk,uk,vk,where f is the vehicle system dynamics; x(k) is the state at the sampling instant k,x=[uk,vk,r]; u(k) is the input to the system at the sampling instant k; w(k) is the process noise which is assumed to be drawn from a zero mean with covariance Q; y(k) is a set of noisy measurements, y=[ax,ay]; v(k) is the observation noise which is assumed to be zero mean Gaussian white noise with covariance R. The process of UKF can be summarized as follows:(a) Sigma points’ calculation: in this section, the sigma points are set by symmetrical sampling method. The mean of the vehicle state vectorx is x- and the covariance is Px. And the elements of sigma matrix X=[x1,x2,…,xk,…,x2n+1] can be generate by(9)xk=x-k=0,x-+n+λPxkk=1,2,…,n,x--n+λPxk-nk=n+1,…,2n, where λ=α2n+κ-n is the scaling parameter. The constant α determines the spread of the sigma points around x- and is usually set as a small value 10-4≤α≤1; κ≥0; it make sure that the covariance matrix is positively definite.(b) These sigma vectors are updated by the system function(10)Xk+1=Xk+1+fXk+1,uk,wk·ΔT.(c) Measurement vector update is as follows:(11)Yk+1=hxk+1,uk+1,vk+1. The updated state vector xk+1|k and the measurement value can be approximated by using a weighted sample of the sigma matrix,(12)xk+1∣k≈∑i=02nWimXk+1i,yk+1≈∑i=02nWimYk+1i.(d) Covariance update: the covariance of the state vector and the measurement value can be calculated by(13)Pk+1∣kx=∑i=02nWicXk+1i-xk+1Xk+1i-xk+1′+Q,Pk+1y=∑i=02nWicYk+1i-yk+1Yk+1i-yk+1′+R,Pk+1xy=∑i=02nWicXk+1i-xk+1Yk+1i-yk+1′.The weight can be acquired by (14)Wim=λn+λ,1/2n+λ,…,1/2n+λ︸2n,Wic=λn+λ+1-α2+γ,1/2n+λ,…,1/2n+λ︸2n.Assume that the system noise and the measurement noise are white Gaussian noise and covariance isQ andR, respectively. γ considers the high order moment of the prior distribution; for Gaussian distribution, here, γ=2 is optimal [27].(e) The correction gain, system state vector, and covariance can be updated by(15)Kk+1=Pk+1xyPk+1y-1,xk+1∣k+1=xk+1∣k+Kk+1zk+1-yk+1,Pk+1∣k+1x=Pk+1∣kx-Kk+1Pk+1yKk+1′, where zk+1 is the measurement value from the vehicle sensor.


### 3.2. Robust Controller Design

The linear bicycle model can be rewritten as(16)X˙=AX+BU,(17)v˙r˙=-kf+krmu-akf+bkrmu-1-akf+bkrIzu-a2kf+b2krIzuvr+krm0akfIz1IzδfΔT,(18)y=r=Cx,(19)Gpsδ→r=kfua+bkr+amusa+b2kfkr+mbkr-akfu2+kf+krIz+ma2kf+b2krus+mIzu2s2,(20)GpsΔT→r=umus+kf+kra+b2kfkr+mbkr-akfu2+kf+krIz+ma2kf+b2krus+mIzu2s2.The zero poles can be deduced as(21)s0δ→r=-a+bkramu;s0ΔT→r=-kf+krmu.And the system frequency and damping ratio are(22)f0=12πa+b2kfkr+mbkr-akfu2mIzu2,β0=kf+krIz+ma2kf+b2kr2mIza+b2kfkr+mbkr-akfu2.For the unit step response, the steady state gain of yaw rate response is as follows:(23)Gpsδ→rlims→0=a+bukrkfa+b2kfkr+bkr-akfmu2,GpsΔT→rlims→0=ukf+kra+b2kfkr+mbkr-akfu2.The controller can be designed as(24)C1s=kδ→r1Gpsδ→rlims→011-s/s0δ→rs2/2πf02+β0/πf0s+111s2/2πfc2+βc/πfcs+1,(25)C2s=kΔT→r1GpsΔT→rlims→011-s/s0ΔT→rs2/2πf02+β0/πf0s+111s2/2πfc2+βc/πfcs+1,where fc and βc can be chosen freely for system pole assignment in order to improve the stability of the system. To improve the yaw rate response, we should increase the bandwidth of the controller, which means frequency fc>f0. To suppress the vibration of second-order system, the damping coefficient should be larger than the controlled system, as βc>β0.The open loop transfer function is(26)Lns=C1sC2sGpsδ→rGpsΔT→r=C1sGpsδ→rC2sGpsΔT→r=kδ→rs2/2πfc2+βc/πfc+1kΔT→rs2/2πfc2+βc/πfc+1=Lδ→rLΔT→r.The close loop transfer function can be deduced as(27)Tns=Lδ→rLδ→r+1LΔT→rLΔT→r+1Τ,Tnslims→0=kδ→rkδ→r+1kΔT→rkΔT→r+1Τ.When the coefficient[kδ→rkΔT→r] is chosen sufficiently large, the steady state of close loop transfer function will go to 1; that is, the yaw rate can follow the reference yaw rate precisely.


### 3.3. Robust Analysis

Vehicle stability heavily depends on the front and rear tire forces. It is appropriate to consider the effect of parameter-based uncertainty by using theH∞ analysis and design the controller to ensure that the system remains stable for all perturbations in the uncertain conditions.The lateral tire force is increasing with the increase of tire slip angle, and when the slip angle is larger than a certain value, the tire force will be saturate, as shown in Figure4. And the tire lateral stiffness varies with the increase of tire slip angle and tire road friction conditions. So we assume that there is multiplicative uncertain of magnitude wI(s). With the uncertainty, the loop transfer function becomes(28)Gps=Gns1+wIsΔIΔI∞≤1,where the subscript I denotes input; Gn(s) is nominal system plant and always assumes that the close loop system is stable; ΔI is system perturbations; wI is system weighting function. It is always difficult to obtain the system perturbations, and ΔI will extend the perturbation range. In special case, ΔI∞=1 is the system perturbation is the worst case; we choose the weighting function wI(s) to avoid the perturbations from deviating the reality of system. Here wI(s) is any stable transfer function.The open loop transfer function becomes(29)Lps=HsGps=GnsHs1+wIsΔI=Ln+wILnΔI, where ΔI is the worst case of the system perturbation ΔI≤1 and wI is a weighting function of perturbation, as shown in Figure 5.Figure 5
Feedback system with multiplicative uncertainty.Considering the Nyquist plot ofLp, at each frequency find the smallest radius lp(w) which includes all the possible perturbation plant:(30)lps=maxGpjw-GnjwGnjw=maxwIsΔI≤wIs,∀w.According to the Nyquist criterion [28], a feedback system is stable if and only if the contour in the Lp(s)-plant does not encircle the (−1, 0) point. The distance from nominal plant Ln to point (−1, 0) must be bigger than the maximum perturbations wILn; that is, (31)wILn<Ln+1,∀w,wILnLn+1<1⟺wITn∞<1,∀w.Combined with (27), the control system robust stability conditions can be deduced as(32)kδ→rkδ→r+1kΔT→rkΔT→r+1Τ<1wI<1lps, where lps=[lp0δ→rlp0ΔT→r].(33)lp0δ→r=lpsδ→rlims→0=mu2akfkfmaxkr-krmin+bkrkrminkfmax-kfkrkfa+b2kfmaxkrmin+bkrmin-akfmaxmu2,lp0ΔT→r=lpsΔT→rlims→0=a+b2kfkr+mbkr-akfu2a+b2kfmaxkrmin+mbkrmin-akfmaxu2kfmax+krminkf+kr-1,lpsδ→rlims→0=kfmaxkf-1,lpsΔT→rlims→0=0.


## 3.1. Sideslip Angle Observer Based on Unscented Kalman Filter

From a practical point of view, sideslip angle is considered as a significant signal in assessment of the transient response of vehicle stability. In this section, the estimation method of sideslip angle of vehicle mass center is presented based on the virtual tire forces sensor, as shown in Figure3. The vehicle state information is estimated by using the UKF, which was proposed by Julier and Uhlman [26]. Detailed explanations of the UKF theory can be found in Ref [27]. The flowchart of the sideslip angle estimator is presented in Figure 4.Figure 3
The flowchart of the sideslip angle estimator.Figure 4
Lateral tire forces under different road friction coefficient.Based on (1)–(3), the vehicle dynamic system can be rewritten in the discrete-time nonlinear state transition equation:(8)xk+1=fxk,uk,wk,yk=hxk,uk,vk,where f is the vehicle system dynamics; x(k) is the state at the sampling instant k,x=[uk,vk,r]; u(k) is the input to the system at the sampling instant k; w(k) is the process noise which is assumed to be drawn from a zero mean with covariance Q; y(k) is a set of noisy measurements, y=[ax,ay]; v(k) is the observation noise which is assumed to be zero mean Gaussian white noise with covariance R. The process of UKF can be summarized as follows:(a) Sigma points’ calculation: in this section, the sigma points are set by symmetrical sampling method. The mean of the vehicle state vectorx is x- and the covariance is Px. And the elements of sigma matrix X=[x1,x2,…,xk,…,x2n+1] can be generate by(9)xk=x-k=0,x-+n+λPxkk=1,2,…,n,x--n+λPxk-nk=n+1,…,2n, where λ=α2n+κ-n is the scaling parameter. The constant α determines the spread of the sigma points around x- and is usually set as a small value 10-4≤α≤1; κ≥0; it make sure that the covariance matrix is positively definite.(b) These sigma vectors are updated by the system function(10)Xk+1=Xk+1+fXk+1,uk,wk·ΔT.(c) Measurement vector update is as follows:(11)Yk+1=hxk+1,uk+1,vk+1. The updated state vector xk+1|k and the measurement value can be approximated by using a weighted sample of the sigma matrix,(12)xk+1∣k≈∑i=02nWimXk+1i,yk+1≈∑i=02nWimYk+1i.(d) Covariance update: the covariance of the state vector and the measurement value can be calculated by(13)Pk+1∣kx=∑i=02nWicXk+1i-xk+1Xk+1i-xk+1′+Q,Pk+1y=∑i=02nWicYk+1i-yk+1Yk+1i-yk+1′+R,Pk+1xy=∑i=02nWicXk+1i-xk+1Yk+1i-yk+1′.The weight can be acquired by (14)Wim=λn+λ,1/2n+λ,…,1/2n+λ︸2n,Wic=λn+λ+1-α2+γ,1/2n+λ,…,1/2n+λ︸2n.Assume that the system noise and the measurement noise are white Gaussian noise and covariance isQ andR, respectively. γ considers the high order moment of the prior distribution; for Gaussian distribution, here, γ=2 is optimal [27].(e) The correction gain, system state vector, and covariance can be updated by(15)Kk+1=Pk+1xyPk+1y-1,xk+1∣k+1=xk+1∣k+Kk+1zk+1-yk+1,Pk+1∣k+1x=Pk+1∣kx-Kk+1Pk+1yKk+1′, where zk+1 is the measurement value from the vehicle sensor.


## 3.2. Robust Controller Design

The linear bicycle model can be rewritten as(16)X˙=AX+BU,(17)v˙r˙=-kf+krmu-akf+bkrmu-1-akf+bkrIzu-a2kf+b2krIzuvr+krm0akfIz1IzδfΔT,(18)y=r=Cx,(19)Gpsδ→r=kfua+bkr+amusa+b2kfkr+mbkr-akfu2+kf+krIz+ma2kf+b2krus+mIzu2s2,(20)GpsΔT→r=umus+kf+kra+b2kfkr+mbkr-akfu2+kf+krIz+ma2kf+b2krus+mIzu2s2.The zero poles can be deduced as(21)s0δ→r=-a+bkramu;s0ΔT→r=-kf+krmu.And the system frequency and damping ratio are(22)f0=12πa+b2kfkr+mbkr-akfu2mIzu2,β0=kf+krIz+ma2kf+b2kr2mIza+b2kfkr+mbkr-akfu2.For the unit step response, the steady state gain of yaw rate response is as follows:(23)Gpsδ→rlims→0=a+bukrkfa+b2kfkr+bkr-akfmu2,GpsΔT→rlims→0=ukf+kra+b2kfkr+mbkr-akfu2.The controller can be designed as(24)C1s=kδ→r1Gpsδ→rlims→011-s/s0δ→rs2/2πf02+β0/πf0s+111s2/2πfc2+βc/πfcs+1,(25)C2s=kΔT→r1GpsΔT→rlims→011-s/s0ΔT→rs2/2πf02+β0/πf0s+111s2/2πfc2+βc/πfcs+1,where fc and βc can be chosen freely for system pole assignment in order to improve the stability of the system. To improve the yaw rate response, we should increase the bandwidth of the controller, which means frequency fc>f0. To suppress the vibration of second-order system, the damping coefficient should be larger than the controlled system, as βc>β0.The open loop transfer function is(26)Lns=C1sC2sGpsδ→rGpsΔT→r=C1sGpsδ→rC2sGpsΔT→r=kδ→rs2/2πfc2+βc/πfc+1kΔT→rs2/2πfc2+βc/πfc+1=Lδ→rLΔT→r.The close loop transfer function can be deduced as(27)Tns=Lδ→rLδ→r+1LΔT→rLΔT→r+1Τ,Tnslims→0=kδ→rkδ→r+1kΔT→rkΔT→r+1Τ.When the coefficient[kδ→rkΔT→r] is chosen sufficiently large, the steady state of close loop transfer function will go to 1; that is, the yaw rate can follow the reference yaw rate precisely.


## 3.3. Robust Analysis

Vehicle stability heavily depends on the front and rear tire forces. It is appropriate to consider the effect of parameter-based uncertainty by using theH∞ analysis and design the controller to ensure that the system remains stable for all perturbations in the uncertain conditions.The lateral tire force is increasing with the increase of tire slip angle, and when the slip angle is larger than a certain value, the tire force will be saturate, as shown in Figure4. And the tire lateral stiffness varies with the increase of tire slip angle and tire road friction conditions. So we assume that there is multiplicative uncertain of magnitude wI(s). With the uncertainty, the loop transfer function becomes(28)Gps=Gns1+wIsΔIΔI∞≤1,where the subscript I denotes input; Gn(s) is nominal system plant and always assumes that the close loop system is stable; ΔI is system perturbations; wI is system weighting function. It is always difficult to obtain the system perturbations, and ΔI will extend the perturbation range. In special case, ΔI∞=1 is the system perturbation is the worst case; we choose the weighting function wI(s) to avoid the perturbations from deviating the reality of system. Here wI(s) is any stable transfer function.The open loop transfer function becomes(29)Lps=HsGps=GnsHs1+wIsΔI=Ln+wILnΔI, where ΔI is the worst case of the system perturbation ΔI≤1 and wI is a weighting function of perturbation, as shown in Figure 5.Figure 5
Feedback system with multiplicative uncertainty.Considering the Nyquist plot ofLp, at each frequency find the smallest radius lp(w) which includes all the possible perturbation plant:(30)lps=maxGpjw-GnjwGnjw=maxwIsΔI≤wIs,∀w.According to the Nyquist criterion [28], a feedback system is stable if and only if the contour in the Lp(s)-plant does not encircle the (−1, 0) point. The distance from nominal plant Ln to point (−1, 0) must be bigger than the maximum perturbations wILn; that is, (31)wILn<Ln+1,∀w,wILnLn+1<1⟺wITn∞<1,∀w.Combined with (27), the control system robust stability conditions can be deduced as(32)kδ→rkδ→r+1kΔT→rkΔT→r+1Τ<1wI<1lps, where lps=[lp0δ→rlp0ΔT→r].(33)lp0δ→r=lpsδ→rlims→0=mu2akfkfmaxkr-krmin+bkrkrminkfmax-kfkrkfa+b2kfmaxkrmin+bkrmin-akfmaxmu2,lp0ΔT→r=lpsΔT→rlims→0=a+b2kfkr+mbkr-akfu2a+b2kfmaxkrmin+mbkrmin-akfmaxu2kfmax+krminkf+kr-1,lpsδ→rlims→0=kfmaxkf-1,lpsΔT→rlims→0=0.


## 4. Tire Force Distribution

The yaw stability mainly depends on the steering angle input and yaw moment control. And the active front steering control is directly input as yaw moment generator. The tire longitudinal tire forces can be calculated by(34)ΔTc=-Bfcosδ2+asinδFxfl-BrFxrl2+Bfcosδ2+asinδFxfr+BrFxrr2,ΔTr=∑Txi,ΔTxi=Fxi·Ri, where ΔTc is the corrected yaw moment generated from differential driving from motors; ΔTr is total required driving torque from driver; Fxi is the ith longitudinal tire force; Ri is the ith tire roll radius.When the slip angle is larger than the upper boundary, the tire forces will be in the saturation region, according to the friction circle theory [29]; the tire longitudinal and lateral forces follow the nonlinear constraint,(35)Fxi_sat2+Fyi_sat2≤μFz2, where μ is the road friction coefficient.So,(36)0≤Fxi_sat≤μFz2-Fyi_sat2.Assuming that the torque split is 40 to 60 percent between front and rear axle, the four wheels driving torque can be obtained according to (34)~(36) with consideration of the tire force friction circle constrain.In order to improve the yaw stability of vehicle, the active front steering system and motor differential driving system work coordinately. The rule of the control coordination can be described as follows: when vehicle lateral acceleration is relative large, it means that tire tends to work close to the saturation zone. The values of yaw control torque will increase by differential longitudinal tire forces, so as to make a good use of tire forces. Likewise, in order to improve the longitudinal stability and traction ability of vehicle, when lateral acceleration is small, the active front steering will work.


## 5. In-Wheel Motor Model and Control

The function of the permanent-magnetic motor is as follows:(37)Te=1τs+1Tm,Tm=Kmotori, where τ is the closed-loop response time which is a control characteristic of motor torque controller; Tm is the motor driving torque command from control algorithm; Kmotor is the motor torque coefficient.(38)V=Ri+Li˙+Kcemfw, where V is the armature voltage, R is the armature resistance, i is the armature current, w is the motor rotational speed, and L is the armature inductance, Kcemf is the counter-electromotive force coefficient. And the control logic of the permanent magnet synchronous motor (PMSM) is shown in Figure 6.Figure 6
Control logic of permanent magnet synchronous motor.


## 6. Simulation and Discussion

The sine-with-dwell maneuver is widely used to evaluate the performance of electronic stability control system [30]. In this section, in order to evaluate the performance of the stability control algorithm, the numerical simulation was implemented in the Matlab/Simulink environment and cosimulated with the Carsim software. Carsim is commercial vehicle dynamics software which is widely used in automotive industry. The performance of robust control is validated under sine-with-dwell maneuver. Figure 7 shows the steering wheel input angle of the sine-with-dwell maneuver. From this simulation scenario, all the dynamic characteristics of the vehicle can be adequately described, and the yaw stability of motor in-wheel EV can be fully verified. In the simulation, a vehicle runs at a constant speed of 50 km/h. The vehicle parameters are listed in Table 1. And the simulation results are plotted in Figures 7– 16.Table 1
Vehicle parameters.
Vehicle mass
m
1351 kg
Yaw moment of inertia
I
z
2031.41 kgm2
C.g. distance to front wheels
a
1.04 m
C.g. distance to rear wheels
b
1.56 m
Sprung c.g. height
h
c
g
0.542 m
Front track width
B
f
1481 m
Rear track width
B
r
1.486 m
The wheel radius
R
w
0.305 mFigure 7
Sine-with-dwell steering wheel input.Figure 8
Estimation results comparison of vertical forces.Figure 9
Estimation results comparison of tire lateral forces.Figure 10
Estimation results comparison of sideslip angle.Figure 11
Estimation results comparison of yaw rate.Figure 12
Sideslip angle in mass center comparison.Figure 13
Lateral acceleration in mass center comparison.Figure 14
Active front steering correction angle.Figure 15
The total yaw moment requirement from controller.Figure 16
Four-wheel torque control input for in-wheel motors.Under the sine-with-dwell simulation scenario, the robust control for vehicle stability combined with the wheel torque distribution strategy should not only meet the demand of vehicle steering maneuverability improvement but also improve the vehicle body stability under high speed condition. Yaw rate response is the maneuverability index, while the sideslip angle of mass center is an important index for vehicle stability [31].In Figure11, it can be clearly observed that compared with no control, the vehicle with the robust control strategy could follow the driver’s intention perfectly, which means the maneuverability of controlled vehicle is improved significantly. Figure 12 compared the sideslip angle of mass center; we can find that the designed estimation strategy can constrain the sideslip angle effectively. This means that the robust controller could improve the vehicle stability especially when vehicle is in high speed or sharp curve turning condition.Figure13 shows the lateral acceleration of vehicle mass center comparison. It is well known that large lateral acceleration brings more load transfer and could lead to inadequate using of the tire forces. With the proposed control algorithm, the lateral acceleration is smaller compared with no control. So the proposed controller has potentials in the improvement of vehicle stability.The proposed systems help improve stability and response at high speed and make the vehicle have a tendency to understeer. An understeering vehicle is self-stabilizing through its tendency to return to straight ahead motion after any external disturbances. Figures14 and 15 plotted the front wheel active steering angle and the total added yaw moment from the controller.Figure16 presents tire force allocation results from the total added yaw moment. It can be found that, combined with the active front steering and motor in-wheel actuators, the vehicle stability is improved significantly. And the controller integrates the longitudinal and lateral dynamic to make good using of tire forces so as to avoid the saturation of tire forces.


## 7. Conclusion and Future Works

This paper proposed theH∞ robust control strategy for motor in-wheel four wheels driving electric vehicles. The proposed control strategy integrates the longitudinal and lateral dynamic control through the active front steering and yaw moment control systems considering the saturation characteristics of tire forces. An upper level control system is designed by using H∞ robust control strategy combined with the online sideslip angle observer to specify a desired yaw moment and correction front steering angle to work on the EVs. The four-wheel torque is obtained optimally by the wheel torque distribution control algorithm. According to the simulation results, it can be asserted that the robust control algorithm could improve vehicle maneuverability and high speed stability. This is very important for vehicle safety. The future works will be focusing on how to distribute the torque smoothly among the axles according to different road friction coefficients or driving situations.



---

*Source: 8197941-2018-02-07.xml*