Dynamic behavior of a vehicle with rear axle compliance steering
Xiao Mei Xu^{1} , Yi Ping Jiang^{2} , Ning Chen^{3} , Heow Pueh Lee^{4}
^{1, 3}College of Automobile and Traffic Engineering, Nanjing Forestry University, Nanjing, China
^{1, 4}Department of Mechanical Engineering, National University of Singapore, Singapore, Singapore
^{2}College of Engineering, Nanjing Agricultural University, Nanjing, China
^{2}Corresponding author
Journal of Vibroengineering, Vol. 19, Issue 6, 2017, p. 44834497.
https://doi.org/10.21595/jve.2017.17580
Received 19 August 2016; received in revised form 15 June 2017; accepted 18 June 2017; published 30 September 2017
Rear axle compliance steering (RACS) is a technology of passive fourwheel steering, which is designed to improve the vehicle handling and stability at medium or high speed. This paper focuses on the dynamic behavior of the vehicle with RACS. Firstly, the compliance steering principle for different rear suspensions is illustrated. Then, the viscoelastic members with fractional order derivative properties are introduced into RACS, and the fractional order model of RACS is formulated. Next, the dynamic model of the vehicle with RACS is established, the adjusting rules for the compliance steering stiffness are derived, and the vehicle stability is investigated. Finally, numerical experiments are performed to illustrate the effects on the vehicle dynamic behavior caused by the compliance steering stiffness, the viscoelastic members and the vehicle longitudinal velocity. Research results show that, the vehicle with RACS has better dynamic characteristics than that without RACS at medium or high speed; and the compliance steering stiffness, the viscoelastic members and the vehicle longitudinal velocity have different impacts on the vehicle lateral dynamic behavior.
Keywords: dynamic behavior, compliance steering, fractional order, sideslip angle, yaw rate.
1. Introduction
As one of effective rear wheel steering modes for passenger cars, rear axle compliance steering (RACS) is popularly adopted to improve the vehicle stability at high speed [1, 2]. RACS utilizes the tire lateral force and the lateral elasticity of the rear suspension rather than any special steering mechanism. This can make the rear axle turn in the same direction as that of the front wheels at medium or high speed, reduce the vehicle sideslip angle at the center of gravity (CG), and thus improve the vehicle lateral stability. These features make RACS applied broadly and widely in various kinds of cars.
From the academic perspective, the research problem for RACS is how to rationally use the compliance characteristics of the rear suspension. Previous studies on this problem can be classified as two streams in terms of the research method: one is experimental method, and the other is numerical method. In the former stream, Momiyama and Miyazaki [2] investigated the compliance under steer technology of the rigid axle rear suspension through calculation and the actual vehicle test; Mou [3] studied the steering characteristics of the compliance suspension used in CitroënZX series cars through experiments; Pan et al. [4] analyzed the rear axle compliance steering performance of the vehicle with fivelink nonindependent rear suspension and verified the results through experiments. In the latter stream, Guo and Yin [1] illustrated the rear axle compliance structure for torsion beam suspension and revealed its compliance steering principle; Li and Chen [5] performed stability and robustness of compliance steering system; Ticã et al. [6] studied the influence of compliance on an elastokinematic model of a proposed rear suspension; Wang et al. [7] carried out a preliminary research on the relationship between compliance steering stiffness and vehicle yaw rate. Compared with experimental method, numerical method is advantageous to simulate and explore the parameter relationships as well as sensitivity analyses from dynamic and evolution perspective.
Generally, the research on RACS is still in development in the academic community. Several previous studies investigated the rear suspension design problem to improve vehicle handling and stability [8, 9], but, to the best of our knowledge, very rare literature investigated the relationships among compliance steering stiffness and dynamic behavior and vehicle working conditions.
The motivation of this study is to investigate the relationships among compliance steering stiffness and vehicle dynamic behavior and vehicle working conditions, by using numerical method. The contributions contain three aspects: (1) incorporate viscoelastic members with fractional order derivative properties into RACS, and then formulate it as a fractional order model; (2) reveal the adjusting rule of the compliance steering stiffness varying with the vehicle longitudinal velocity; (3) perform the influence of the viscoelastic member’s characteristic parameters on the vehicle dynamic behavior.
The structure of this paper is organized as follows. In Section 2, the operating principle of RACS is illustrated. In Section 3, the dynamic model of the vehicle with RACS is formulated and discussed. In Section 4, the numerical experiments are performed and analyzed. Finally, the concluding remarks are provided in Section 5.
2. Operating principle of RACS
This section mainly illustrates the operating principle of RACS for different rear suspensions. Even though RACS has different implementation ways due to different types of suspensions, it has the same technical principle. Through ingenious design of the connection structure between the rear axle and the vehicle body, RACS makes the rotation center of the rear axle pertain to the vehicle body move backward, and then by means of the lateral tire force compels the rear axle to rotate through an angle.
The technology of RACS was firstly used in those cars with torsion beam rear suspensions, such as Citroën Fukang, Citroën Elysée, Volkswagen Golf IV, the Audi A3 and Bora. Fig. 1 shows the arrangement of the rubber bushings on the rear axle frame of the torsion beam suspension [1]. Suppose that the vehicle turns left, the lateral force ${F}_{r}$ suffered by the rear wheels points to the left. Because the two front bushings are elastic in $A$$A$ direction and rigid in BB direction, they can only move in $A$$A$ direction and their instantaneous velocities intersect at the point $cc$. The distance from point $cc$ to the acting point of ${F}_{r}$ is denoted by $e$. With the action of ${F}_{r}$ the whole rear axle frame rotates counterclockwise around the point $cc$, namely the rear axle turns in the same direction as that of the front wheels, which intensifies the trend of understeer and reduces the equivalent sideslip angles of the rear axle.
Fig. 1. Top view of RACS schematic diagram for torsion beam suspension [1]
For the double wishbone suspension, the compliance steering of rear wheels is actualized by using different stiffness rubber bushings which connect the $A$shaped control arm to the vehicle body [10].
Fig. 2 shows the compliance steering principle of the rear axle for fivelink suspension [4]. The lower drag links $EF$ and $GH$ are parallel to the vehicle longitudinal axis. The intersect point $cc$ of the upper drag links $AC$ and $BD$ is the yaw rotation center of the rear axle. The transversal push rod $JK$ lies behind the rear axle, and the two ends of rod $JK$ are connected to the vehicle body through rubber bushings. The radial stiffness of rubber bushings used in rod $JK$ is almost five times as those bushings used in drag links $AC$ and $BD$. When the vehicle turns a corner, the lateral force of the rear axle is mainly undertaken by the rubber bushings used in rod $JK$. And then the whole rear axle turns in the same direction as the front wheels. Obviously, the positions of rods $AC$, $BD$ and $JK$ are crucial for determining the steering characteristics of the vehicle.
Fig. 2. Top view of RACS schematic diagram for fivelink suspension [4]
3. Theoretical formulation
In this section, the fractional order model of RACS is formulated, and the dynamic model of the vehicle with RACS is established and discussed.
3.1. Fractional order model of RACS
From Section 2 it can be seen that, elastic bushings are key components for actualizing rear axle compliance steering. In the present study, in order to investigate the feasibility of adjusting the compliance steering angle through changing the viscoelastic material parameters, the viscoelastic members replacing elastic bushings are introduced into RACS. Fig. 3 shows the schematic diagram of RACS with viscoelastic members. As shown in Fig. 3, $l$ is the axial length of the viscoelastic member; $c$ is the distance from the axis of the viscoelastic member to the rear axle; $cc$ is the yaw rotation center of the rear axle and $e$ is the distance from the point $cc$ to the acting point of lateral force ${F}_{r}$.
Fig. 3. Schematic diagram of RACS with viscoelastic members
A generalized standard linear solid model with two independent fractional parameters is adopted to describe the behavior of the viscoelastic member [11], as given in Eq. (1):
where $\sigma \left(t\right)$ is the stress, $\epsilon \left(t\right)$ is the strain, ${\tau}_{\epsilon}$ is the relaxation time, ${\tau}_{\sigma}$ is the retardation (creep) time, ${E}_{0}$ is the relaxed magnitude of the elastic modulus, and ${\tau}_{\epsilon}$, ${\tau}_{\sigma}$ and ${E}_{0}$ are parameters dependent on the material. Parameters $\alpha $ and $\gamma $$(0<\alpha <1$, $0<\gamma <1)$ are fraction, ${D}^{\alpha}(\cdot )$ and ${D}^{\gamma}(\cdot )$ respectively represent RiemannLiouville fractional derivatives with the order $\alpha $ and $\gamma $, which are defined in Eq. (2) [12]:
In which, $\xi $ is the fractional order of the time derivative, ${D}^{\xi}$ represents the operator of fractional derivation of $\xi $th order, $\mathrm{\Gamma}$ denotes the Euler gamma function [13].
Assume that the viscoelastic member is homogeneous and let the two sides of Eq. (1) multiply by the cross section area $A$ of the viscoelastic member, then the Eq. (3) can be obtained:
In which, ${F}_{v}\left(t\right)$ is the axial force of the viscoelastic member caused by the compliance steering.
According to Fig. 3, the strain $\epsilon \left(t\right)$ of the viscoelastic member can be written as:
where ${\delta}_{c}\left(t\right)$ is the compliance steering angle of the rear axle.
Ignoring the rotational inertia of the RACS, and considering the equilibrium relation between the rear wheel lateral force ${F}_{r}$ and the axial force of the viscoelastic member ${F}_{v}\left(t\right)$, the Eq. (3) can be derived as:
where ${C}_{s}$ is the compliance steering stiffness. Implementing Laplace transform of Eq. (5), the compliance steering angle ${\delta}_{c}\left(s\right)$ of the rear axle can be derived as:
3.2. Dynamic model of the vehicle with RACS
Fig. 4 shows a twodegreeoffreedom model including yaw and lateral motion dynamics of the vehicle with RACS.
Fig. 4. Dynamic model of the vehicle with RACS
In this model the coordinate system is fixed on the vehicle center of gravity which is denoted as $CG$. $x$ and $y$ denote the vehicle longitudinal and lateral directions, respectively; ${F}_{f}$, ${F}_{r}$ represent the front and rear tire lateral forces, respectively; ${u}_{f}$, ${u}_{r}$ represent the velocities of the front and rear axles, respectively; ${\alpha}_{f}$, ${\alpha}_{r}$ represent the sideslip angles of the front and rear tires, respectively; ${\delta}_{f}$ is the steering angle of the front wheel and ${\delta}_{c}$ is the compliance steering angle of the rear axle; ${\delta}_{e}$ is the angle between the $x$ axis and the velocity direction of ${u}_{r}$, which is usually called equivalent sideslip angle and its numerical value equals to $({\alpha}_{r}{\delta}_{c})$; $u$ and $v$ denote the longitudinal and lateral velocities of the vehicle $CG$, respectively, and $u$ is supposed to be constant in the present study; $\beta $ is the vehicle sideslip angle at $CG$ and $\beta =v/u$; ${\omega}_{r}$ is the yaw rate of the vehicle.
Assume that the sideslip angle $\beta $ is small and $\left\beta \right\ll 1$, then the sideslip angels of front and rear tires can be written as:
where $a$ and $b$ are the distance from $CG$ to the front and rear axles respectively, and $L=a+b$ is the wheelbase.
The tire lateral force of the front and rear wheels are considered as linear functions of their sideslip angles, which can be written as:
where ${k}_{f}$ and ${k}_{r}$ represent the cornering stiffness of the front and rear tires, respectively.
According to Fig. 4 and considering the previous force, the vehicle model including lateral and yaw motions can be expressed as:
where $m$ is the vehicle mass, ${I}_{z}$ denotes the vehicle moment of inertia about the yaw axis.
3.3. Adjusting rules for compliance steering stiffness ${\mathit{C}}_{\mathit{s}}$
Implementing Laplace transform of Eq. (9), two transfer functions ${G}_{\beta}\left(s\right)$ and ${G}_{{\omega}_{r}}\left(s\right)$ from the front wheel step steering input ${\delta}_{f}$ to the sideslip angle $\beta $ and the yaw rate ${\omega}_{r}$ are yielded as Eq. (10) and Eq. (11):
where:
$T=K=\frac{2au{k}_{f}}{{I}_{z}},Z=\frac{4b{k}_{f}{k}_{r}(a+b)+2a{k}_{f}m{u}^{2}}{mu{I}_{z}}\cdot \frac{{C}_{s}}{{k}_{r}}.$
Sideslip angle $\beta $ of the vehicle $CG$ is an important parameter for evaluating the vehicle lateral stability. Zero sideslip angle ($\beta =$ 0) is always the control object for the vehicle with active or semiactive steering system. $\beta =$ 0 is helpful to improve the vehicle lateral stability. According to Eq. (10), Eq. (12) can be obtained and which indicates that the steadystate value of sideslip angle $\beta $ for the vehicle with RACS is influenced by the vehicle mass, the position of $CG$, the cornering stiffness of front and rear tires, the vehicle longitudinal velocity and the compliance steering stiffness of RACS. Therefore, when vehicle working conditions change, the steadystate value of the sideslip angle $\beta $ can be kept unchanged through adjusting the compliance steering stiffness ${C}_{s}$:
According to Eq. (12) and Table 1, when $u=$ 20 m/s Fig. 5 can be obtained, which depicts the relationship between the steadystate value of ${G}_{\beta}\left(s\right)$ and compliance steering stiffness ${C}_{s}.$ In Fig. 5, the steadystate value of ${G}_{\beta}\left(s\right)$ changes nonlinearly with ${C}_{s}.$ When ${C}_{s}$ equals to some value, the steadystate value of ${G}_{\beta}\left(s\right)$ equals to zero. That means changing ${C}_{s}$ can make the steadystate value of the sideslip angle equal to zero, which is helpful to improve the lateral stability of the vehicle with RACS.
Fig. 5. Relationship between steadystate value of ${G}_{\beta}\left(s\right)$ and ${C}_{s}$
Let the numerator of Eq. (12) equal to zero, the relationship between the longitudinal velocity $u$ and the compliance steering stiffness ${C}_{s}$ can be obtained as:
Fig. 6 shows the adjusting rule of the compliance steering stiffness ${C}_{s}$ changing with the vehicle longitudinal velocity $u$, where ${C}_{s}$ decreases nonlinearly with $u$ increasing. Especially, when $u$ lies in [18, 40] m/s (i.e., [64.8, 144] km/h), ${C}_{s}$ has better adjusting effect; when $u$ is larger than 40 m/s, ${C}_{s}$ has very little changes; when $u$ is lower than 18 m/s, ${C}_{s}$ approaches to infinite. Therefore, at the case of lower speed, RACS should stop working, and the vehicle should switch into the common working mode which turns a corner with two front wheels.
Fig. 6. Adjusting rule for the compliance steering stiffness ${C}_{s}$
3.4. Investigation of the vehicle stability
In order to assess the ability of the vehicle dynamic system withstanding the disturbance, the stability of the vehicle dynamic system is investigated in this subsection.
According to Eq. (5)Eq. (9), the state space model of the vehicle motion can be expressed as:
where $X$ is the matrix of the state variables consisting of the yaw rate ${\omega}_{r}$ and the sideslip angle $\beta $. $X$ can be expressed as: $X={\left[\begin{array}{ll}{\omega}_{r}& \beta \end{array}\right]}^{T}.$
The coefficient matrices in Eq. (14) are listed as follows:
In which, $R=1/(PW+1)$ and $W=(1+{\tau}_{\sigma}^{\gamma}{s}^{\gamma})/(1+{\tau}_{\epsilon}^{\alpha}{s}^{\alpha})$.
According to Lyapunov’s method [14], the Lyapunove function for verifying the stability of Eq. (14) can be given as:
where $P$ is a 2×2 matrix, $I$ is the unit matrix. If the eigenvalue of $P$ are all positive, the system defined in Eq. (14) is stable; otherwise, the system is unstable.
The key parameters of the vehicle and tires used in this study are listed in Table 1.
Table 1. Vehicle and tire parameters [15]
Parameters

Values

Unit

Parameters

Values

Unit

$m$

1740

kg

${k}_{f}$

–35000

N/rad

$a$

1.035

m

${k}_{r}$

–37500

N/rad

$b$

1.655

m

${I}_{z}$

3048

kg·m^{2}

Based on the above Lyapunove stability criteria, some numerical experiments under different working conditions and sensitive parameters given in Table 2 are carried out to verify the system stability against the perturbation and disturbance. Research results show that the vehicle dynamic systems are all stable under different numerical experiments. In other words, the vehicle with RACS has strong ability to resist the perturbation and disturbance.
Table 2. Working conditions and sensitive parameters
Working conditions

Sensitive parameters


$u=$ 30 (m/s)
$\alpha =$ 0.3
$\gamma =$ 0.7

$a$ (m)


0.935

1.035

1.135

1.235

1.335

1.435

1.535

1.635


$a=$ 1.035 (m)
$\alpha =$ 0.3
$\gamma =$ 0.7

$u$ (m/s)


5

10

15

20

25

30

35

40


$a=$ 1.035 (m)
$u=$ 30 (m/s)
$\gamma =$ 0.5

$\alpha $


0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8


$a=$ 1.035 (m)
$u=$ 30 (m/s)
$\alpha =$ 0.5

$\gamma $


0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

4. Numerical experiments
In this section, numerical experiments are carried out to investigate the dynamic behavior of the vehicle with RACS, and the effects on the vehicle dynamic behavior caused by the compliance steering stiffness ${C}_{s}$, the viscoelastic members used in RACS and the vehicle longitudinal velocity $u$. All performed numerical experiments are conducted with the same step steering input ${\delta}_{f}$. The output variables of the vehicle dynamic system are yaw rate ${\omega}_{r}$, sideslip angle $\beta $ and compliance steering angle ${\delta}_{c}$.
4.1. Effects of compliance steering stiffness
In this part, the effects of compliance steering stiffness on the dynamic behavior of the vehicle with RACS are investigated. For all numerical cases in this part, the related parameters of the vehicle are set as $u=$ 20 m/s, $\alpha =$ 0.3, $\gamma =$ 0.7 and the values of ${\tau}_{\epsilon}^{\alpha}$ and ${\tau}_{\sigma}^{\gamma}$ are referred to Eldred et al. [16]. And ${C}_{s}/{k}_{r}$ is used to measure the compliance steering stiffness and the values of ${C}_{s}/{k}_{r}$ are set as –1.5, –3, –4.5, respectively. Note that, the compliance steering stiffness remarkably grows with the increasing of ${C}_{s}/{k}_{r}$ absolute value.
Fig. 7 shows the effects of ${C}_{s}/{k}_{r}$ on the sideslip angle $\beta $ of the vehicle $CG$. The steadystate value of sideslip angle for the vehicle with RACS is obviously smaller than the vehicle without RACS, which means the vehicle with RACS has better vehicle body posture when it turns a corner. Moreover, the steadystate value of $\beta $ is more and more close to zero with the decreasing of $\left{C}_{s}/{k}_{r}\right.$ In other words, reducing the compliance steering stiffness can dramatically make the sideslip angle $\beta $ decreased to zero.
Fig. 8 shows the effects of ${C}_{s}/{k}_{r}$ on the yaw rate ${\omega}_{r}$ of the vehicles. It can be seen that the steadystate value of yaw rate for the vehicle with RACS is obviously smaller than the vehicle without RACS. The steadystate value of yaw rate is obviously decreasing with the reduction of the compliance steering stiffness. It should be noted that the smaller yaw rate means the driver needs to turn more steering wheel angle for the same turning radius.
Actually, the vehicle without RACS is a special case of the proposed model of vehicles with RACS. As the curves demonstrated in Fig. 7 and Fig. 8, when ${C}_{s}/{k}_{r}$ changes from –1.5 to –4.5, the curves of vehicle with RACS are approaching to those of vehicle without RACS. At the extreme case, let compliance steering stiffness equal to a large enough value, the curves would have overlapped each other between vehicles with RACS and without RACS. Furthermore, in Fig. 7 and Fig.8, the curves of vehicle without RACS are consistent with the results given in [17, 18]. Therefore, it could be concluded that, the proposed model of vehicle with RACS in this study is more generalized and includes the dynamic behavior of the vehicle without RACS.
Fig. 9 shows the effects of ${C}_{s}/{k}_{r}$ on the compliance steering angle ${\delta}_{c}$. The compliance steering angle ${\delta}_{c}$increases with $\left{C}_{s}/{k}_{r}\right$ decreasing, namely reducing the compliance steering stiffness can increase the compliance steering angle. Furthermore, the compliance steering angles are all small and within the deformation range for the viscoelastic members used in RACS. Therefore, it is completely feasible to change the compliance steering angle through adjusting the compliance steering stiffness.
Fig. 7. Effects of ${C}_{s}/{k}_{r}$ on the sideslip angle
Fig. 8. Effects of ${C}_{s}/{k}_{r}$ on the yaw rate
Fig. 9. Effects of ${C}_{s}/{k}_{r}$ on the compliance steering angle
Fig. 10 and Fig. 11 depict the effects of ${C}_{s}/{k}_{r}$ on the bode diagrams of sideslip angle and yaw rate transfer functions. Bode diagram includes magnitudefrequency characteristic and phasefrequency characteristic. In the present study, magnitudefrequency characteristic reflects the ability of the vehicle to accurately perform the steering instructions at different frequencies. Phasefrequency characteristic reflects the phase difference between the vehicle dynamic response and the front wheel steering angle input varying with the frequency. Considering the vehicle steering process at medium speed or high speed belongs to a low frequency motion and the working frequency for the driver turns the steering wheel is always within 1 rad/s [17], this study focuses on the steering frequency between 0.1 rad/s and 10 rad/s.
In Fig. 10, the magnitudefrequency characteristic shows the varying of sideslip angle gain with the steering frequency, where the magnitudefrequency curve of the vehicle without RACS is obviously higher than those of the vehicle with RACS. With the increasing of steering frequency, the phase difference between the sideslip angle and the steering angle input is decreasing. Besides, the phase difference of the vehicle without RACS is generally a little larger than those of the vehicle with RACS.
In Fig. 11, the magnitudefrequency characteristic shows the curves almost have the same steadystate gain of the yaw rate. For the curve of the vehicle without RACS, there is a peak point, which means at that frequency a little distortion appears for the vehicle performing the driver’s steering instructions. For the curves of ${C}_{s}/{k}_{r}=3$ and ${C}_{s}/{k}_{r}=$ –45, they both have good performance in carrying out the steering motion; for the curve of ${C}_{s}/{k}_{r}=$ –1.5, the magnitude of the yaw rate is decreasing with the steering frequency increasing. Summarily, the compliance steering stiffness has little impact on the phasefrequency characteristic.
Fig. 10. Effects of ${C}_{s}/{k}_{r}$ on the bode diagram of ${G}_{\beta}\left(s\right)$
Fig. 11. Effects of ${C}_{s}/{k}_{r}$ on the bode diagram of ${G}_{{\omega}_{r}}\left(s\right)$
4.2. Effects of viscoelastic member
This subsection investigates effects of the viscoelastic member used in RACS on vehicle dynamic behavior. In the numerical cases, the vehicle longitudinal velocity $u=$ 20 m/s, the fractional orders $\alpha $ and $\gamma $ are set as $\alpha =$ 0.2, 0.4, and $\gamma =$ 0.6, 0.8, respectively. The parameters $\alpha $ and $\gamma $ are related to the material properties of the viscoelastic member. The relationship between ${C}_{s}$ and $u$ is in accordance with Eq. (13).
Fig. 12 and Fig. 13 show the effects of fractional orders $\alpha $ and $\gamma $ on the sideslip angle and yaw rate of the vehicle with RACS. It can be observed that $\alpha $ has little effect on the sideslip angle and nearly no effect on the yaw rate. The fractional order $\gamma $ has obvious effect on the sideslip angle and a little effect on the yaw rate. Furthermore, the larger the order $\gamma $ is, the closer to zero the sideslip angle is, and the smaller the yaw rate is. Fig. 14 depicts the effects of fractional orders $\alpha $ and $\gamma $ on the compliance steering angle. It is easy to see that the order $\alpha $ has little effect and the order $\gamma $ has obvious effect on the compliance steering angle. The larger the order $\gamma $ is, the larger the compliance steering angle is.
Fig. 12. Effects of $\alpha $ and $\gamma $ on the sideslip angle
Fig. 13. Effects of $\alpha $ and $\gamma $ on the yaw rate
Fig. 14. Effects of $\alpha $ and $\gamma $ on the compliance steering angle
Obviously, the two fractional parameters $\alpha $ and $\gamma $ of the viscoelastic member have distinct effects on the vehicle dynamic behavior, which provides an important guidance for the choice of viscoelastic member used in RACS. With the increasing of $\gamma $ the sideslip angle is more and more close to zero, which is more and more beneficial for improving the vehicle lateral stability. The effects of increasing order $\gamma $ on the vehicle dynamic behavior are similar to those of reducing compliance steering stiffness.
4.3. Effects of vehicle longitudinal velocity
In this subsection, the effects of vehicle longitudinal velocity $u$ on the dynamic behavior of the vehicle with RACS are investigated. For all numerical cases in this part, the fractional orders are set as $\alpha =$ 0.3 and $\gamma =$ 0.7, $u$ are set as 20 m/s, 30 m/s, 40 m/s and 50 m/s, respectively. The relationship between ${C}_{s}$ and $u$ is in accordance with Eq. (13).
Fig. 15 and Fig. 16 show the effects of $u$ on the sideslip angle and the yaw rate. From the two figures, it can be seen that the curves present obvious fluctuation phenomenon before reaching the steadystate values when $u$ is larger than 40 m/s. This phenomenon is related to the structure parameters of the vehicle. Fig. 15 indicates that for different vehicle velocity, if the compliance steering stiffness is adjusted according to Eq. (13), the sideslip angle of vehicle $CG$ can always tend to be zero. Fig. 16 shows that steadystate yaw rate decreases obviously with the increasing of $u$, which means at higher speed the driver needs to steer more steering wheel angle for the same turning radius.
Fig. 15. Effects of $u$ on the sideslip angle
Fig. 16. Effects of $u$ on the yaw rate
Effects of $u$ on the compliance steering angle are shown in Fig. 17. Before reaching the steadystate value, the curves also experience slight fluctuation at high speed, and the larger the velocity is, the more obvious the fluctuation is. The smaller compliance steering stiffness will make the compliance steering angle appear unstable before reaching its steadystate value.
In Fig. 18 and Fig. 19, the curves demonstrate how $u$ impacts the bode diagrams of ${G}_{\beta}\left(s\right)$ and ${G}_{{\omega}_{r}}\left(s\right)$. For the cases of $u$ equals to 20 m/s, 30 m/s, 40 m/s, 50 m/s, the adjusting rules of ${C}_{s}$ are obtained by Eq. (13). In Fig. 18, the magnitudefrequency curves show the vehicle with RACS has different steadystate gains of sideslip angles at different $u$. Along with the decreasing of $u$, the magnitude of ${G}_{\beta}\left(s\right)$ is also decreasing within the working frequency, which means the sideslip angle is gradually close to zero. The phasefrequency curves show that within the working frequency the phase difference between the input and the output has little variation for different $u$. In Fig. 19, the magnitude of yaw rate is inversely proportional to $u$, which is compatible with the time domain results in Fig. 16. The magnitudefrequency curves are slightly decreasing with the increasing frequency. In other words, at the same $u$, the vehicle with RACS has good steering performance within the working frequency. Moreover, the phasefrequency curves have little difference at different $u$ within the working frequency. When the steering frequency equals to 0.1 rad/s the phase delay is small, which means the vehicle has good steering response performance when the driver turns the steering wheel at low speed.
Fig. 17. Effects of $u$ on the compliance steering angle
Fig. 18. Effects of $u$ on the bode diagram of ${G}_{\beta}\left(s\right)$
Fig. 19. Effects of $u$ on the bode diagram of ${G}_{{\omega}_{r}}\left(s\right)$
5. Conclusions
In this paper, the dynamic behavior of the vehicle with RACS has been investigated. Firstly, the operating principle of RACS has been illustrated. Secondly, the fractional order model of RACS with viscoelastic members has been formulated, and the dynamic model of the vehicle with RACS has been proposed. Thirdly, the relationship between the steadystate value of sideslip angle transfer function and the compliance steering stiffness ${C}_{s}$, the adjusting rule of ${C}_{s}$ has been analyzed, and the vehicle stability has been investigated. Fourthly, the impacts on the vehicle dynamic behavior have been investigated from the three aspects, i.e., the compliance steering stiffness, the viscoelastic members and the vehicle longitudinal velocity. Finally, many numerical experiments have been conducted in terms of time domain and frequency domain, respectively.
Some findings can be summarized as follows: (1) if the vehicle runs at lower speed, RACS should stop working and the vehicle should switch into the common working mode which turns with two front wheels; (2) the compliance steering angle is inversely proportional to the compliance steering stiffness, and the sideslip angle and the yaw rate are proportional to the compliance steering stiffness; (3) the vehicle with RACS has strong ability to resist the perturbation and disturbance. (4) compared with the vehicle without RACS, the vehicle with RACS has lower sideslip angle, yaw rate and a little larger response time; (5) for the dynamic behavior of the vehicle with RACS, the viscoelastic parameter $\gamma $ has obvious impact, while the parameter $\alpha $ has tiny impact.
The most direct natural extension involves carrying out physical experiments to verify the analytical results and numerical simulation in this study. Other promising directions of future research are the investigation on control strategies and control algorithms of the compliance steering stiffness.
Acknowledgements
Authors of the paper would like to express appreciation of the supports from National Natural Science Foundation of China (51605228), Jiangsu Natural Science Foundation (BK20160742), Youth Science and Technology Innovation Project of Nanjing Forestry University (CX2016028), Jiangsu Postdoctoral Research Foundation (1402084C).
References
 Guo K. H., Yin C. L. The steering principle and handling stability analysis of compliance suspension. Automotive Engineering, Vol. 19, Issue 6, 1997, p. 321324, (in Chinese). [Search CrossRef]
 Momiyama F., Miyazaki K. Compliance steer and road holding of rigid rear axle for enhancing the running straightness of large sized vehicles. International Truck and Bus Meeting and Exposition, 1993. [Publisher]
 Mou X. D. The Preliminary Study of Side Leans Steering the Characteristic to Passenger Vehicle Compliance Suspension. SAE Technical Paper, 2007013639. [Search CrossRef]
 Pan X., Lin Y., Chen X. K., Wang D. C. Rear axle compliance steering performance analysis of fivelink dependent rear suspension. Journal of Southwest Jiaotong University, Vol. 47, Issue 3, 2012, p. 477483., (in Chinese). [Search CrossRef]
 Li B., Chen S. H. A research on the analysis and synthesis method of compliance steering. Automotive Engineering, Vol. 28, Issue 9, 2006, p. 829833. [Search CrossRef]
 Ticã M., Dobre G., Mateescu V., VirziMariotti G. Influence of compliance for an elastokinematic model of a proposed rear suspension. International Journal of Automotive Technology, Vol. 15, Issue 6, 2014, p. 885891. [Publisher]
 Wang Y., Chen N., Tian J., Xu X. M. Handing performances of vehicle with a fractional compliant rearwheel steering system. IEEE International Conference on Fractional Differentiation and Its Applications (ICFDA), 2014, p. 16. [Publisher]
 Blundell M. V. The influence of rubber bush compliance on vehicle suspension movement. Materials and Design, Vol. 19, Issues 12, 1998, p. 2937. [Publisher]
 Gao X., Jiang D., Wang Y. The research of effects of suspension bushing stiffness on vehicle yaw response performance. International Conference on Electric Information and Control Engineering, 2011, p. 62836286. [Search CrossRef]
 Afkar A., MahmoodiKaleibar M., Paykani A. Geometry optimization of double wishbone suspension system via genetic algorithm for handling improvement. Journal of Vibroengineering, Vol. 14, Issue 2, 2012, p. 827837. [Search CrossRef]
 Rossikhin Y. A., Shitikova M. V. Analysis of rheological equations involving more than one fractional parameters by the use of the simplest mechanical systems based on these equations. Mechanics of TimeDependent Materials, Vol. 5, Issue 2, 2001, p. 131175. [Publisher]
 Zhou Y. Basic Theory of Fractional Differential Equations. World Scientific, Singapore, 2014. [Publisher]
 Ross B. A Brief History and Exposition of the Fundamental Theory of Fractional Calculus Fractional Calculus and its Applications. Springer Berlin Heidelberg, 1975. [Publisher]
 Khalil, H. K. Nonlinear Systems. Third Edition, Upper Saddle River, Prenticehall, 2002. [Search CrossRef]
 Yin G. D., Chen N., Wang J. X., Wu L.Y. A study on μsynthesis control for fourwheel steering system to enhance vehicle lateral stability. Journal of Dynamic Systems Measurement and Control, Vol. 133, Issue 1, 2011, p. 1100211006. [Publisher]
 Eldred L. B., Baker W. P., Palazotto A. N. KelvinVoigt versus fractional derivative model as constitutive relations for viscoelastic materials. AIAA Journal, Vol. 33, Issue 3, 2015, p. 547550. [Search CrossRef]
 Guiggiani M. The Science of Vehicle Dynamics. Springer, Netherlands, 2014. [Publisher]
 Bai Y., Ren X., Xu Y. F., Wu F. W. Factors affecting dynamic characteristics of driver steering control. Journal of Chang’an University, Vol. 35, Issue 4, 2015, p. 112118. [Search CrossRef]
Cited By
Applied Mathematical Modelling
Yusheng Zhou, Kwokwai Chung

2022

Science and Engineering of Composite Materials
Sun Songsong, Wan Maosong

2021

Shock and Vibration
Xiaomei Xu, Dong Chen, Lei Zhang, Ning Chen

2019

Shock and Vibration
Maohua Xiao, Kai Wen, Cunyi Zhang, Xiao Zhao, Weihua Wei, Dan Wu

2018
