Next Article in Journal
A New Consideration for Validating Battery Performance at Low Ambient Temperatures
Next Article in Special Issue
Simulation of the Oxygen Reduction Reaction (ORR) Inside the Cathode Catalyst Layer (CCL) of Proton Exchange Membrane Fuel Cells Using the Kinetic Monte Carlo Method
Previous Article in Journal
Mixed Kernel Function Support Vector Regression with Genetic Algorithm for Forecasting Dissolved Gas Content in Power Transformers
Previous Article in Special Issue
Analytical Calculation of Armature Reaction Field of the Interior Permanent Magnet Motor
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Steering Stability Control for a Four Hub-Motor Independent-Drive Electric Vehicle with Varying Adhesion Coefficient

1
National Engineering Laboratory for Electric Vehicles, Beijing Institute of Technology, Beijing 100081, China
2
Co-Innovation Center of Electric Vehicles in Beijing, Beijing Institute of Technology, Beijing 100081, China
3
BAIC BJEV Inc., Beijing 100021, China
*
Author to whom correspondence should be addressed.
Energies 2018, 11(9), 2438; https://doi.org/10.3390/en11092438
Submission received: 30 August 2018 / Revised: 10 September 2018 / Accepted: 11 September 2018 / Published: 14 September 2018
(This article belongs to the Collection Electric and Hybrid Vehicles Collection)

Abstract

:
In order to enhance the steering stability of a four hub-motor independent-drive electric vehicle (4MIDEV) on a road with varying adhesion coefficient, for example on a joint road, this paper proposes a hierarchical steering stability control strategy adapted to the road adhesion. The upper control level of the proposed strategy realizes the integrated control of the sideslip angle and yaw rate in the direct yaw moment control (DYC), where the influences of the road adhesion and sideslip angle are both studied by the fuzzy control. The lower control level employs a weight-based optimal torque distribution algorithm in which weight factors for each motor torque are designed to accommodate different adhesion of each wheel. The proposed stability control strategy was validated in a co-simulation of the Carsim and Matlab/Simulink platforms. The results of double-lane-change maneuver simulations under different conditions indicate that the proposed strategy can effectively achieve robustness to changes in the adhesion coefficient and improve the steering stability of the 4MIDEV.

1. Introduction

With the rapid development of motor technology, the four hub-motor independent-drive electric vehicle (4MIDEV) concept has emerged [1,2]. Different from the traditional centralized drive vehicle, the 4MIDEV is directly driven by four hub motors without the need for clutches, drive shafts and other mechanical components. The drive/brake torque of four hub motors not only responds quickly and accurately, but also can be independently controlled, which helps a lot to improve the stability and maneuverability of the vehicles [3]. Since traditional electronic stability program (ESP) cannot be directly applied to the 4MIDEV, its steering stability control has attracted much attention from researchers [4,5,6,7].
The stability control strategy with hierarchical structure for the 4MIDEV has been widely used and proved to be effective, including an upper control level to generate the yaw moment command, and a lower control level to change the yaw moment to the torque required for four drive wheels [5,6,7,8]. However, few studies have enabled the stability control strategy to adapt to the adhesion coefficient, but only focus on the constraints of road adhesion. Furthermore, relevant verifications are always neglected when the adhesion coefficient changes, which may bring more control errors and adversely affect the final steering stability control effect. As the driving conditions of vehicles become more complex, it is necessary to study a stability control strategy adaptive to road adhesion to ensure good steering stability of the 4MIDEV, which puts forward higher requirements for the upper control level and the lower control level.
In general, the upper control level is usually designed based on the direct yaw moment control (DYC) which controls the lateral motion of vehicle through a yaw moment. The sideslip angle and yaw rate corresponding to stability and maneuverability are always selected as the control targets in the DYC [9,10,11,12]. However, unlike the yaw rate, there is no direct relationship between the sideslip angle and the yaw moment, but they are related to each other by nonlinear tire forces. Most of the DYCs previously proposed only control the yaw rate and ignore the impact of the sideslip angle. For example, Tahami et al. [13] and Lin et al. [14] adopted fuzzy control and sliding mode control, respectively, in the DYC, and both successfully realized the control of the yaw rate without the consideration of the sideslip angle. However, the importance of the sideslip-angle control is much greater than that of the yaw-rate control when the adhesion coefficient is low. The lack of consideration for the sideslip angle will be detrimental to the steering stability control effect, especially when the road adhesion coefficient changes. Some studies have improved the DYCs on the basis of the previous ones in which the sideslip angle was only taken as an intermediate variable in the control of the yaw rate, and finally achieved good results. Zou et al. [15] and Tahami et al. [16] redesigned the reference yaw rate based on the sideslip angle respectively, and both finally achieved a good control of the sideslip angle and the yaw rate by tracking the yaw rate with the reference value. Zhai et al. [17] weighted the sideslip angle and the yaw rate to obtain a comprehensive control quantity and also achieved good control performance. Although the above methods have basically considered the influence of the sideslip angle, the impact of adhesion coefficient was neglected, and there was a lack of correlation verification when the adhesion coefficient changes. With different adhesion coefficients, the effects of the sideslip angle and the yaw rate on the stability are quite different. Therefore, a new DYC adapted to road adhesion is required to achieve the integrated control of sideslip angle and yaw rate.
The torque distribution algorithm is a key part in the lower control level. The optimal distribution algorithm has been widely adopted and proved to be effective compared with other torque distribution algorithms such as average distribution [18], etc., which makes the allocation to meet optimal objectives as much as possible. In current studies, the tire workload usage is always selected as the main optimal objective to reflect steering stability. Jin et al. [19] and Park et al. [20] developed an optimal distribution algorithm based on the tire workload usage objective respectively, and the steering stability of the vehicle was improved. In [21], Zhai et al. designed an energy consumption objective in addition to tire workload usage objective, and finally achieved energy savings while ensuring stability through multi-objective optimization. Furthermore, Li et al. [22] additionally considered the control objective and ensured the accuracy of stability control. These methods have been proven to be effective in the face of the road with constant adhesion coefficient, but they did not consider the influence of changes in adhesion coefficient. For the stability control, the workload usage of each tire was generally considered equally important to each other, i.e., different adhesion conditions at each wheel were ignored, which may result in a wheel with low adhesion experiencing extreme loads while other wheels maintain a relatively high stability margin.
To solve the above problems, this paper proposes a hierarchical stability control strategy that is adaptive to the road adhesion coefficient. The upper control level of the proposed strategy realizes the integrated control of the sideslip angle and yaw rate in which the influences of the sideslip angle and road adhesion are both studied by the fuzzy control. As the road adhesion coefficient is also took into account in the yaw moment control, the output of the DYC becomes more reasonable. The lower control level employs a weight-based optimal distribution algorithm in which weight factors for each tire workload usage are designed to accommodate different adhesion of each wheel. The different weight control of each wheel torque allows the 4MIDEV to better adapt to changes in adhesion coefficient.
The rest of this study is organized as follows: in Section 2, the 4MIDEV system model, including the vehicle dynamics model, driver model, etc., is presented. The proposed hierarchical strategy for the 4MIDEV is then designed in Section 3. In Section 4, the proposed strategy is evaluated by a co-simulation using CarSim and Matlab/Simulink under different double-lane-change maneuver conditions. Finally, the conclusions are summarized in Section 5.

2. System Overview and Modeling

The 4MIDEV, as shown in Figure 1, is directly driven by four hub motors without the need for clutches, drive shafts and other mechanical components. The drive/brake torque of four hub motors can be independently controlled by a main electronic controller using the CAN bus (a kind of area network), according to the driver’s operation and actual vehicle condition.

2.1. Vehicle Dynamics Model

Ignoring the pitch and roll motions, a simplified three-degree-of-freedom (3-DOF) dynamic model is elaborated to study the vehicle’ planar motion: yaw motion, longitudinal motion, and lateral motion, as shown in Figure 2a, where V ˙ x and V ˙ y denote the derivative of the longitudinal velocity V x and the lateral velocity V y , respectively, γ ˙ is the derivative of the yaw rate γ . m denotes the mass of the vehicle, I z is the yaw inertia of the vehicle. F x i j and F y i j denote the longitudinal and lateral forces of the respective wheels, respectively, where i { f , r } denotes the front or the rear and j { l , r } denotes the left or the right. δ f denotes the front wheel steering angle. a and b are the distance between the center of gravity and the front axle and rear axle, respectively. d is half of the tread.
The dynamic equations of the 3-DOF model can be expressed in matrix form as:
[ V ˙ x V ˙ y γ ˙ ] = [ V y γ V x γ 0 ] + B x [ F x f l F x f r F x r l F x r r ] + B y [ F y f l F y f r F y r l F y r r ]
where the coefficient matrices B x and B y are:
B x = [ cos δ f m cos δ f m 1 m 1 m sin δ f m sin δ f m 0 0 a sin δ f d cos δ f I z a sin δ f + d cos δ f I z d I z d I z ]
B y = [ sin δ f m sin δ f m 0 0 cos δ f m cos δ f m 1 m 1 m a cos δ f + d sin δ f I z a cos δ f d sin δ f I z b I z b I z ]

2.2. Tire Model

The tire model reflects the nonlinear characteristics of tire dynamics, which is a critical part for steering stability control. The “Magic Formula” tire model [23] and “Unified Index” tire model [24] have been widely used and proven to be effective, but too much characteristic parameters of the tire are required. In order to improve the real-time performance of the system, a simplified tire model in [21] is used as follows:
F y i j = K α i j 1 ( F x i j μ i j F z i j ) 2 μ i j k tan 1 ( k μ i j α i j ) k = K α i j π 2 F z i j
where μ i j denotes the road adhesion coefficient, K α i j and α i j denote the lateral stiffness and slip angle of each tire, respectively.

2.3. Hub Motor Model

The hub motors typically employ a permanent magnet synchronous motor (PMSM) that can be modeled using a space vector control method. Considering that the torque dynamic response of the PMSM is much faster than the vehicle dynamic response, the dynamic response of the motor can be simplified to a second-order system whose torque transfer function can be expressed as:
G ( s ) = T i j T i j * = 1 2 τ 2 s 2 + 2 τ s + 1
where τ denotes the damping ratio, and its value is initially set to 0.01. T i j * denotes the motor drive torque command, and T i j is the actual drive torque.

2.4. Driver Model

A “preview-follow” driver model [25] is adopted in this paper to simulate the driver’s operations on the acceleration/brake pedal and steering wheel based on the desired vehicle speed and the expected route. The driver’s preview time is set to 1 s for the closed-loop test. Then the front-wheel-steer angle can be obtained.

3. Design of Hierarchical Steering Stability Control Strategy for the 4MIDEV

In this paper, the stability control strategy proposed for the 4MIDEV is designed as a hierarchical structure, including a reference output level, an upper control level with the yaw moment control, and a lower control level based on the torque allocation, as shown in Figure 3. A reference dynamic model is established in the reference output level to obtain the desired vehicle dynamics parameters. The upper control level is composed of speed control and yaw moment control to realize adaptive control of the vehicle speed, yaw rate, and sideslip angle, respectively. The lower control level optimally distributes the virtual control forces generated by the upper control level, such as the yaw moment, to each hub-motor controller in the form of the driving/brake torque commands, so as to realize the steering stability control of the 4MIDEV.

3.1. Reference Output Level

To directly and clearly reflect the driver’s steering intention, the reference output level simplifies the above 3-DOF model to 2 DOFs as the reference dynamics model, as shown in Figure 2b, which makes the steady-state responses of the vehicle linear with the driver input. The kinetic equations can be formulated as follows:
{ m V x ( β ˙ + γ ) = F y f cos δ f + F y r I z γ ˙ = a F y f cos δ f b F y r
where γ ˙ and β ˙ are the derivative of the yaw rate γ and the sideslip angle β , respectively. F y f and F y r denote the lateral tire force of the front and rear wheels.
By ignoring the system transient characteristics, the desired target values for the sideslip angle and the yaw rate can be obtained:
{ γ = 1 1 + A V x 2 V x l δ f β = 1 ( m L ) ( a b K r ) V x 2 1 + A V x 2 b l δ f
where A = m l 2 a K f b K r K f K r , K f and K r are the front and rear tire cornering stiffness, respectively.
Considering the above target values appear in the linear working area of the tire, the target values of the yaw rate and the sideslip angle must be limited when the tire works on the nonlinear working area [26]:
{ γ limit = 0.85 V x μ g β limit = tan 1 ( 0.02   μ g )
Thus, the target values of the yaw rate and the sideslip angle obtained from the reference output level can be expressed as:
{ γ d e s = min { | γ | , γ lim } sgn ( δ f ) β d e s = min { | β | , β lim } sgn ( δ f )

3.2. Upper Control Level

The upper control level, as shown in Figure 3, consists of speed control and yaw moment control, which generates the virtual commands of traction force and yaw moment to realize the integrated control of vehicle speed, yaw rate, and sideslip angle. Furthermore, the road adhesion coefficient is also taken into account in the yaw moment control, and the control effect becomes better.

3.2.1. Speed Control

The sliding mode control is adopted in the speed control to compute the traction force required to maintain the desired longitudinal speed, as shown in Figure 3, which takes the control error of speed as the sliding surface:
S v x = V x V x d e s
where V x and V x d e s denote the actual and desired longitudinal speed, respectively. The following formula can be derived from Equation (10):
S ˙ v x = V ˙ x V ˙ x d e s = 1 m [ F x ( F y f l + F y f r ) sin δ ] + V y γ V ˙ x d e s
The reaching law is defined as:
S ˙ v x = ε v x s a t ( S v x ) k v x S v x
where ε v x > 0 , k v x > 0 , and the saturation function s a t ( S v x ) can be formulated as:
s a t ( S v x ) = { sgn ( x / ϕ ) | s | ϕ x / ϕ | s | < ϕ
The longitudinal-speed control law can be finally obtained by combining Equations (11) and (12):
F x d e s = m [ V ˙ x d e s V y γ ε v x s a t ( S v x ) k v x S v x ] + ( F y f l + F y f r ) sin δ

3.2.2. Yaw Moment Control

The sideslip angle and the yaw rate are both selected as the control variables in the yaw moment control, as shown in Figure 3. Considering that the sideslip angle and the yaw rate are mutually coupled, and the sideslip angle is inconvenient to measure directly and its estimation accuracy is limited, the yaw moment control takes the sideslip angle as an intermediate variable in the control of the yaw rate. Based on the estimation of the sideslip angle, a fuzzy controller is developed in the yaw moment control to compute the correction factor K , so as to adjust the nominal yaw rate γ 0 m , i.e., γ 0 m = K γ . The yaw moment is then calculated by a sliding mode controller based on the deviation of the nominal yaw rate γ 0 m from the desired yaw rate γ d e s . The structure of the yaw moment control is shown in Figure 4.
Based on the phase plane method [27] and the β-method [28], a comprehensive control quantity, i.e., β + B 1 β ˙ , is designed as an input of the fuzzy controller, in which the coefficient B 1 reflects the influence of the rate of sideslip angle and is related to the adhesion coefficient. In addition, since the influence of the sideslip angle on the steering stability is different under different road adhesion coefficients, the adhesion coefficient μ is also introduced into the design of the correction coefficient K as another input of the fuzzy controller. There are seven fuzzy subsets to describe the values of inputs and outputs, i.e., Positive-Big (PB), Negative-Big (NB), Positive-Medium (PM), Negative-Medium (NM), Positive-Small (PS), Negative-Small (NS), and Zero (Z). The fuzzy rules are shown in Table 1.
The sliding mode controller, as shown in Figure 4, computes the yaw moment so that the nominal yaw rate γ 0 m follows the desired value γ d e s . The sliding surface is defined as:
S γ = γ 0 m γ d e s
Then the following formula can be derived from (15):
S ˙ γ = γ ˙ 0 m γ ˙ d e s = M z I z + 1 I z [ ( F y f l F y f r ) sin δ d ( F y r l + F y r r ) b + ( F y f r + F y f l ) cos δ a ] γ ˙ d e s
The reaching law is:
S ˙ γ = ε γ s a t ( S γ ) k γ S γ
where ε γ > 0 , k γ > 0 . The yaw moment control law can be finally formulated as:
M z d e s = ( F y f l F y f r ) sin δ d + ( F y r l + F y r r ) b ( F y f r + F y f l ) cos δ a + I z [ γ ˙ d ε y s a t ( S y ) k y S y ]

3.3. Lower Control Level

The lower control level optimally distributes the virtual control forces generated by the upper control level, such as the traction force and the yaw moment, to each hub-motor controller in the form of the driving/brake torque commands. The relationship between the upper level and lower level can be formulated as:
v = Bu
where u = [ T f l T f r T r l T r r ] T , v = [ F x d e s M z d e s ] T , B = 1 R [ cos δ f cos δ f 1 1 d cos δ f + a sin δ f d cos δ f + a sin δ f d d ] .
In addition to the virtual control forces from the upper control level, the torque distribution also needs to meet the friction ellipse constraint:
F x i j 2 + F y i j 2 ( μ i j F z i j ) 2
where F z i j and μ i j denote the vertical load and the road adhesion coefficient at each wheel, respectively. The lateral force can be calculated from the simplified tire model.
The longitudinal force F x i j from the ground is related to the torque T i j applied to the wheels:
ω ˙ i j = 1 J c ( T i j F x i j R M f i j )
where ω ˙ i j , J c and R denote the angular acceleration, moment of inertia and rolling radius of the wheel, respectively. M f i j is the rolling resistance moment.
The lower control level employs a weight-based optimal distribution algorithm in which weight factors for each tire workload usage were designed to accommodate different adhesion of each wheel.
A weight-based optimal distribution algorithm is investigated in the lower control level with the tire workload usage Ω i j as the main optimization target, which can be formulated as:
Ω i j = F x i j 2 μ i j 2 F z i j 2
The smaller the sum of tire workload usages, the greater the stability margin of the wheel, and the better the steering stability of the vehicle [17,29]. Considering that the vehicle is facing a complex road surface, the adhesion conditions at each wheel are often different. If the tire workload usages are simply added as the overall optimization target, the different adhesion conditions will be ignored. There is a possibility that the wheels in low adhesion conditions are unstable, and the wheels in the high adhesion condition still have a large stability margin. The weight factors for each tire workload usage are designed to accommodate different adhesion of each wheel as follows:
ψ i j = χ i j Ω i j Ω i j
where χ i j reflects the influence of the adhesion coefficient at each wheel, the value of which is positively correlated with the adhesion coefficient.
In summary, the torque distribution algorithm of the lower control level can be described as:
min J 1 = Γ u 2
where u = [ T f l T f r T r l T r r ] T , Γ = d i a g ( ψ i j μ i j F z i j R ) . The constraints for the torque distribution are:
{ Bu = v ( I I ) u ( u ¯ u - )
where u ¯ and u - denote the upper and lower limits of u , respectively. Considering that the intensity of the equality constraint B u = v is too strong, it is converted to min W v ( Bu v ) 2 and introduced into (24) as a penalty item with a weight coefficient ξ as follows:
min J 2 = Γ u 2 + ξ W v ( Bu v ) 2
where W v denotes the distribution weight matrix, which is related to the control priority of the traction force and the yaw moment. For the low-adhesion road, the weights of the yaw moment and the traction force can be increased and reduced, respectively, so as to obtain good stability performance.
The (26) can be reformulated as:
J 2 = arg min u - u u ¯ ( ξ 1 2 W v B Γ 1 ) u ( ξ 1 2 W v v 0 ) 2 2 = Au C 2 2
The above optimization problem can be solved by the active set method [30]. For the kth iteration, at a given u k , the correction value p k needs to be determined:
min p k J = A ( u k + p k ) C 2 2 B p k = 0
If the u k + p k is a feasible solution, then u k + 1 = u k + p k , and the Lagrange multiplier ( ε λ ) T can be calculated by:
A T ( A u k + 1 C ) = ( B T D T ) ( ε λ )
where D = ( I I ) T . If λ 0 , then u k + 1 is the final optimal solution and the iteration ends. Otherwise, the constraint corresponding to the minimum λ is removed and the next iteration is performed.
If the u k + p k is not a feasible solution, then u k + 1 = u k + α k p k , and the next iteration is performed. The step α k is limited by:
α k = max { α k [ 0 , 1 ] ; u - u k + α k p k u ¯ }
The iteration continues in accordance with the above steps until an optimal solution is found.

4. Simulation Test and Analysis

The proposed adaptive stability control strategy was verified in a co-simulation of Carsim and Matlab/Simulink. The 4MIDEV model and driver model were developed in the Carsim, while the proposed strategy and the hub-motor model were implemented in the Simulink, as shown in the Figure 3. To further prove the improvement in steering stability, the proposed hierarchical stability control strategy, referred to as “adaptive control,” was compared with other two strategies, referred to as “common control” and “speed control”, respectively. The “common control” is proposed in [17], in which the DYC and torque distribution algorithm didn’t consider the change of the adhesion coefficient. In addition, the “speed control” has only developed a “preview-follow” driver model and a speed sliding mode controller to ensure the stability. The aforementioned stability strategies were evaluated at double-lane-change (DLC) maneuvers under the μ-split road and joint road. The path setting for the DLC maneuver is shown in Figure 5. The parameters of the 4MIDEV and hub motors used in this study are presented in Table 2.

4.1. μ-Split Road

The DLC maneuver was first conducted on the μ-split road, as shown in Figure 6. The desired vehicle speed was set to 60 km/h and the driver’s preview time was 1 s. The dynamic responses of the vehicle under different control strategies are shown in Figure 7.
The adhesion coefficient of μ-split road is shown in Figure 7a, which was initially set to 0.75 and suddenly became 0.75 on the half side and 0.1 on the other half side at a distance of 105 m from the starting point. Figure 7b,c show the vehicle’s track and speed responses, respectively. It can be seen that the adaptive control, compared with the other two controls, performed best in achieving the following of the trajectory and speed with the smallest errors.
Figure 7d,e show the yaw rate and the sideslip angle responses of the vehicle, respectively. Both the sideslip angle and the yaw rate of aforementioned controls appeared to deviate from the desired values when the adhesion coefficient was changed, i.e., about at 7 s. However, the adaptive control could make these two quantities basically follow the desired values. The control error of the sideslip angle and the yaw rate under the adaptive control, compared with the common control, are reduced by 58.8% and 55%, respectively, which indicates that the steering stability can be ensured even in the face of varying adhesion coefficient. Figure 7f shows the sideslip angle rate and the sideslip angle in the phase plane. The curve of the adaptive control, compared with the other two controls, is most concentrated on the origin, which indicates that the stability control effect is best [21,31].
Figure 7g–i show the torque of each hub motor under the speed control, common control, and adaptive control, respectively. It can be seen that the adaptive control, compared with the other two controls, can effectively adapt to the change of the adhesion coefficient and realize different weight control of each wheel torque. For example, the output torque of the 1th and 3th wheel shown in the Figure 7i, were both increased to ensure the desired yaw moment.

4.2. Joint Road

A joint road was also developed for the DLC maneuver, as shown in Figure 8, in which the vehicle maintained a constant speed of 60 km/h. The dynamic responses of the vehicle under different stability control strategies are shown in Figure 9.
Figure 9a shows the adhesion coefficient of joint road, which was initially 0.75 and then changed to 0.2 at the middle of DLC maneuver. Figure 9b,c show the vehicle’s track and speed responses, respectively. It can be seen that all the tracks of the aforementioned controls were obviously offset from the scheduled track after the change of road adhesion coefficient. In particular, the track under the speed control had a large fluctuation and tended to be unstable, whereas the adaptive control could suppress the disturbances with less deviation and fast response. In addition, both the adaptive control and the common control could basically ensure the desired speed, but the speed under the speed control was significantly deviated from the expected value.
Figure 9d,e show the yaw rate and the sideslip angle responses of the vehicle, respectively. Figure 9f shows the phase plane about the sideslip angle rate and the sideslip angle. Although the change of the adhesion coefficient brought disturbances to the control, the adaptive control can effectively reduce the control error of the yaw rate and the sideslip angle, showing the best stability control effect. In contrast, the vehicle under the speed control was already unstable. The maximum control error of the sideslip angle and yaw rate under the common control was about 3 times and 4.1 times larger than that under the adaptive control, respectively. Furthermore, the curves of the speed control, common control, and adaptive control in the Figure 9f, tended to be more focused on the origin, which means that the stability control effect of the aforementioned controls increased gradually.
Figure 9g–i show the torque of each hub motor under the speed control, common control, and adaptive control, respectively. After the change of adhesion coefficient, the torque under speed control had an unstable oscillation, whereas the adaptive control could realize different weight control of each wheel torque, which made the 4MIDEV better adapt to changes in adhesion coefficient.

5. Conclusions

In order to enhance the steering stability of a four hub-motor independent-drive electric vehicle, this paper proposes a hierarchical stability control strategy that is adaptive to the road adhesion coefficient. The upper control level of the proposed strategy is composed of speed control and yaw moment control, which realizes adaptive control of the vehicle speed, yaw rate, and sideslip angle, respectively. The influences of the sideslip angle and road adhesion are both studied by fuzzy control in the yaw moment control, and the output of the DYC becomes more reasonable. The lower control level employs a weight-based optimal distribution algorithm in which weight factors for each tire workload usage are designed to accommodate different adhesion of each wheel. The different weight control of each wheel torque allows the 4MIDEV to better adapt to changes in adhesion coefficient.
The proposed stability control strategy was validated in a co-simulation using Carsim and Matlab/Simulink, in which double-lane-change maneuvers under the μ-split road and joint road conditions were conducted. The simulation results indicate that the proposed strategy, compared with other methods, can effectively achieve robustness to changes in the adhesion coefficient and improve the steering stability of the 4MIDEV. Since the proposed method adopts sliding mode control, the chattering phenomenon is an important factor affecting the control effect. Due to the characteristics of the sliding mode control, the chattering cannot be eliminated, which may make the stability control effect become worse when the adhesion coefficient changes very rapidly and greatly. However, since the vehicle motion is slow enough compared to the sampling time, this situation hardly occurs in practice. Therefore, the good control effect of the proposed method can be basically guaranteed. In future work, the hierarchical stability control strategy will be tested through the hardware-in-the-loop test, and then the real-car test.

Author Contributions

Conceptualization, L.Z. and T.S.; Funding acquisition, L.Z.; Investigation, R.H.; Methodology, R.H. and T.S.; Validation, R.H.; Writing–original draft, R.H.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 51475045.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, X.D.; Göhlich, D. Integrated traction control strategy for distributed drive electric vehicles with improvement of economy and longitudinal driving stability. Energies 2017, 10, 126. [Google Scholar] [CrossRef]
  2. Liu, M.C.; Gu, F.H.; Zhang, Y.Z. Ride comfort optimization of in-wheel-motor electric vehicles with in-wheel vibration absorbers. Energies 2017, 10, 1647. [Google Scholar] [CrossRef]
  3. Fallah, S.; Khajepour, A.; Fidan, B.; Chen, S.K.; Litkouhi, B. Vehicle optimal torque vectoring using state-derivative feedback and linear matrix inequality. IEEE Trans. Veh. Technol. 2013, 62, 1540–1552. [Google Scholar] [CrossRef]
  4. Shuai, Z.B.; Zhang, H.; Wang, J.M.; Li, J.Q.; Ouyang, M.G. Combined AFS and DYC control of four-wheel-independent-drive electric vehicles over CAN network with time-varying delays. IEEE Trans. Veh. Technol. 2014, 63, 591–602. [Google Scholar] [CrossRef]
  5. Chen, Y.; Wang, J.M. Design and evaluation on electric differentials for overactuated electric ground vehicles with four independent in-wheel motors. IEEE Trans. Veh. Technol. 2012, 61, 1534–1542. [Google Scholar] [CrossRef]
  6. Kim, D.H.; Kim, C.J.; Kim, S.H.; Choi, J.Y.; Han, C.S. Development of adaptive direct yaw-moment control method for electric vehicle based on identification of yaw-rate model. In Proceedings of the 2011 IEEE Intelligent Vehicles Symposium (IV), Baden-Baden, Germany, 5–9 June 2011; pp. 1098–1103. [Google Scholar]
  7. Chen, J.; Yu, J.Z.; Zhang, K.X.; Ma, Y. Control of regenerative braking systems for four-wheel-independently-actuated electric vehicles. Mechatronics 2018, 50, 394–401. [Google Scholar] [CrossRef]
  8. Wang, Z.P.; Wang, Y.C.; Zhang, L.; Liu, M.C. Vehicle stability enhancement through hierarchical control for a four-wheel-independently-actuated electric vehicle. Energies 2017, 10, 947. [Google Scholar] [CrossRef]
  9. Ding, S.; Liu, L.; Zheng, W.X. Sliding mode direct yaw-moment control design for in-wheel electric vehicles. IEEE Trans. Ind. Electron. 2017, 64, 6752–6762. [Google Scholar] [CrossRef]
  10. Ren, B.T.; Chen, H.; Zhao, H.Y.; Yuan, L. MPC-based yaw stability control in in-wheel-motored EV via active front steering and motor torque distribution. Mechatronics 2016, 38, 103–114. [Google Scholar] [CrossRef]
  11. Alipour, H.; Sabahi, M.; Bannae Sharifian, M.B. Lateral stabilization of a four wheel independent drive electric vehicle on slippery roads. Mechatronics 2015, 30, 275–285. [Google Scholar] [CrossRef]
  12. Zhang, X.Z.; Wei, K.X.; Yuan, X.F.; Tang, Y.Q. Optimum torque distribution for stability improvement of four-wheel distributed driven electric vehicle using coordinated control. J. Comput. Nonlinear Dyn. 2016, 11. [Google Scholar] [CrossRef]
  13. Tahami, F.; Kazemi, R.; Farhanghi, S. A novel driver assist stability system for all-wheel-drive electric vehicles. IEEE Trans. Veh. Technol. 2003, 52, 683–692. [Google Scholar] [CrossRef]
  14. Lin, C.; Xu, Z.F. Wheel Torque Distribution of four-wheel-drive electric vehicles based on multi-objective optimization. Energies 2015, 8, 3815–3831. [Google Scholar] [CrossRef]
  15. Zou, G.C.; Luo, Y.G.; Lian, X.M.; Li, K.Q. A research of DYC for independent 4WD EV based on control target dynamic regulated. In Proceedings of the 2007 IEEE International Conference on Vehicular Electronics and Safety, Beijing, China, 13–15 December 2008; pp. 1–7. [Google Scholar]
  16. Tahami, F.; Farhangi, S.; Kazemi, R. A fuzzy logic direct yaw-moment control system for all-wheel-drive electric vehicles. Veh. Syst. Dyn. 2004, 41, 203–221. [Google Scholar] [CrossRef]
  17. Zhai, L.; Sun, T.M.; Wang, J. Electronic stability control based on motor driving and braking torque distribution for a four in-wheel motor drive electric vehicle. IEEE Trans. Veh. Technol. 2016, 65, 4726–4739. [Google Scholar] [CrossRef]
  18. De Novellis, L.; Sorniotti, A.; Gruber, P.; Orus, J.; Rodriguez Fortun, J.-M.; Theunissen, J.; De Smet, J. Direct yaw moment control actuated through electric drivetrains and friction brakes: Theoretical design and experimental assessment. Mechatronics 2015, 26, 1–15. [Google Scholar] [CrossRef] [Green Version]
  19. Jin, X.J.; Yin, G.D.; Chen, N. Gain-scheduled robust control for lateral stability of four-wheel-independent-drive electric vehicles via linear parameter-varying technique. Mechatronics 2015, 30, 286–296. [Google Scholar] [CrossRef]
  20. Park, J.; Park, Y. Optimal input design for fault identification of overactuated electric ground vehicles. IEEE Trans. Veh. Technol. 2016, 65, 1912–1923. [Google Scholar] [CrossRef]
  21. Zhai, L.; Hou, R.F.; Sun, T.M.; Kavuma, S. Continuous steering stability control based on an energy-saving torque distribution algorithm for a four in-wheel-motor independent-drive electric vehicle. Energies 2018, 11, 350. [Google Scholar] [CrossRef]
  22. Li, B.; Goodarzi, A.; Khajepour, A.; Chen, S.K.; Litkouhi, B. An optimal torque distribution control strategy for four-independent wheel drive electric vehicles. Veh. Syst. Dyn. 2015, 53, 1172–1189. [Google Scholar] [CrossRef]
  23. Bakker, E.; Pacejka, H.B.; Lidne, L. A New Tire Model with an Application in Vehicle Dynamics Studies. J. Passeng. Car. 1989, 98, 101–113. [Google Scholar]
  24. Li, S.Y.; Min, Y.J.; Wang, L.M.; An, L.H. Establishment and Simulation Analysis of Tire Dynamic Model. J. Nanjing Inst. Technol. (Nat. Sci. Ed.) 2009, 7, 34–38. [Google Scholar]
  25. Guo, K.H.; Pan, F.; Ma, F.J. Preview optimized artificial neural network driver model. Chin. J. Mech. Eng. 2003, 39, 26–28. [Google Scholar] [CrossRef]
  26. Song, P.; Tomizuka, M.; Zong, C.F. A novel integrated chassis controller for full drive-by-wire vehicles. Veh. Syst. Dyn. 2015, 53, 215–236. [Google Scholar] [CrossRef]
  27. Shimada, K.; Shibahata, Y. Comparison of Three Active Chassis Control Methods for Stabilizing Yaw Moments. J. Passeng. Car. 1994, 103, 1178–1187. [Google Scholar]
  28. Inagaki, S.; Kushiro, I.; Yamamoto, M. Analysis on vehicle stability in critical cornering using phase-plane method. In Proceedings of the International symposium on advanced vehicle control, Tsukuba, Japan, 24–28 October 2014. [Google Scholar]
  29. Mokhiamar, O.; Abe, M. Simultaneous optimal distribution of lateral and longitudinal tire forces for the model following control. J. Dyn. Syst. Meas. Control 2004, 126, 753–763. [Google Scholar] [CrossRef]
  30. Harkegard, O. Efficient active set algorithms for solving constrained least squares problems in aircraft control allocation. In Proceedings of the 41st IEEE Conference on Decision and Control, Las Vegas, NV, USA, 10–13 December 2002; pp. 1295–1300. [Google Scholar]
  31. Zhang, H.; Li, X.S.; Shi, S.M.; Liu, H.F.; Guan, R.; Liu, L. Phase Plane analysis for vehicle handling and stability. Int. J. Comput. Intell. Syst. 2011, 4, 1179–1186. [Google Scholar] [CrossRef]
Figure 1. System Configuration of the 4MIDEV.
Figure 1. System Configuration of the 4MIDEV.
Energies 11 02438 g001
Figure 2. Dynamic models of the 4MIDEV. (a) 3-DOF model. (b) 2-DOF model.
Figure 2. Dynamic models of the 4MIDEV. (a) 3-DOF model. (b) 2-DOF model.
Energies 11 02438 g002
Figure 3. The configuration of hierarchical stability control strategy.
Figure 3. The configuration of hierarchical stability control strategy.
Energies 11 02438 g003
Figure 4. Structure of the yaw moment control.
Figure 4. Structure of the yaw moment control.
Energies 11 02438 g004
Figure 5. The path of the DLC maneuver.
Figure 5. The path of the DLC maneuver.
Energies 11 02438 g005
Figure 6. The μ-split road. (a) top view. (b) driver view.
Figure 6. The μ-split road. (a) top view. (b) driver view.
Energies 11 02438 g006
Figure 7. Simulation results for DLC on μ-split road (μ = 0.75–0.2). (a) adhesion coefficient. (b) track. (c) speed. (d) yaw rate. (e) sideslip angle. (f) phase plane. (g) torque under the speed control. (h) torque under the common control. (i) torque under the adaptive control.
Figure 7. Simulation results for DLC on μ-split road (μ = 0.75–0.2). (a) adhesion coefficient. (b) track. (c) speed. (d) yaw rate. (e) sideslip angle. (f) phase plane. (g) torque under the speed control. (h) torque under the common control. (i) torque under the adaptive control.
Energies 11 02438 g007aEnergies 11 02438 g007b
Figure 8. The joint road. (a) top view. (b) driver view.
Figure 8. The joint road. (a) top view. (b) driver view.
Energies 11 02438 g008
Figure 9. Simulation results for DLC on joint road (μ = 0.75–0.2). (a) adhesion coefficient. (b) track. (c) speed. (d) yaw rate. (e) sideslip angle. (f) phase plane. (g) torque under the speed control. (h) torque under the common control. (i) torque under the adaptive control.
Figure 9. Simulation results for DLC on joint road (μ = 0.75–0.2). (a) adhesion coefficient. (b) track. (c) speed. (d) yaw rate. (e) sideslip angle. (f) phase plane. (g) torque under the speed control. (h) torque under the common control. (i) torque under the adaptive control.
Energies 11 02438 g009aEnergies 11 02438 g009b
Table 1. Fuzzy rules of correction coefficient K.
Table 1. Fuzzy rules of correction coefficient K.
K β + B 1 β ˙
NBNMNSZPSPMPB
μ NBPBPBPBPMPBPBPB
NMPBPBPMPMPMPBPB
NSPMPMPSPSPSPMPM
ZNMNSNSZPSPSPM
PSNMNMNMNSNSNSNM
PMNBNBNMNMNMNBNB
PBNBNBNBNMNBNBNB
Table 2. Parameters of the 4MIDEV and hub motors.
Table 2. Parameters of the 4MIDEV and hub motors.
SymbolQuantityValue
mvehicle mass1411 kg
Izmoment of inertia about the yaw axis2031.4 kg m3
hgheight of the center of mass0.54 m
blength from the center of gravity (CG) to the rear wheel axis1.56 m
alength from the CG to the front wheel axis1.04 m
rtire radius0.3 m
dtread width1.48 m
Pmmaximum power28 kW
Tmmaximum torque340 Nm
nmmaximum speed1200 rpm
Perated power14 kW
Terated torque170 Nm
nerated speed800 rpm

Share and Cite

MDPI and ACS Style

Hou, R.; Zhai, L.; Sun, T. Steering Stability Control for a Four Hub-Motor Independent-Drive Electric Vehicle with Varying Adhesion Coefficient. Energies 2018, 11, 2438. https://doi.org/10.3390/en11092438

AMA Style

Hou R, Zhai L, Sun T. Steering Stability Control for a Four Hub-Motor Independent-Drive Electric Vehicle with Varying Adhesion Coefficient. Energies. 2018; 11(9):2438. https://doi.org/10.3390/en11092438

Chicago/Turabian Style

Hou, Rufei, Li Zhai, and Tianmin Sun. 2018. "Steering Stability Control for a Four Hub-Motor Independent-Drive Electric Vehicle with Varying Adhesion Coefficient" Energies 11, no. 9: 2438. https://doi.org/10.3390/en11092438

APA Style

Hou, R., Zhai, L., & Sun, T. (2018). Steering Stability Control for a Four Hub-Motor Independent-Drive Electric Vehicle with Varying Adhesion Coefficient. Energies, 11(9), 2438. https://doi.org/10.3390/en11092438

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop