Research Article

Advances in Mechanical Engineering

2021, Vol. 13(5) 1–25
Ó The Author(s) 2021
Analysis and optimization of the DOI: 10.1177/16878140211015523
journals.sagepub.com/home/ade
vehicle handling stability with
considering suspension kinematics and
compliance characteristics

Jin Gao and Fuquan Wu

Abstract
The dynamic model of the front double wishbone suspension and the rear multi-link suspension of the vehicle are estab-
lished. On the basis of detailed analysis of suspension kinematics, calculation method of wheel alignment angle and force
calculation of suspension bushing, the influence mechanism of suspension bushing on the vehicle transient state is clari-
fied, and the vehicle transient characteristic index is derived from the vehicle three-free dynamic model. The sensitivity
analysis of the suspension bushing is carried out, and the bushing stiffness which has a great influence on the transient
state of the vehicle is obtained. The bushing stiffness scale factor is used as the optimization variable, the vehicle transi-
ent characteristic index is used as the optimization target, and the NSGA-II optimization algorithm is used for multi-
objective optimization. After optimization, one Pareto solution is selected to compare with the original vehicle, the com-
parison results show that the yaw rate gain, resonance frequency and delay time of yaw rate in the vehicle transient char-
acteristic index are all improved, other optimization targets change less. In the steady-state comparison, the understeer
tendency of the vehicle increases, and the roll angle of the vehicle increases but is within an acceptable range.

Keywords
Suspension kinematics, suspension compliance, bushing deformation, vehicle transient characteristics, multi-objective
optimization

Date received: 7 December 2020; accepted: 14 April 2021

Handling Editor: James Baldwin

Introduction bushing. At present, many experts and scholars have

carried out a lot of research on suspension bushing,
The K (Kinematics) characteristic and C (Compliance) suspension movement and suspension force.
characteristic of the suspension are important perfor- Garcı́a et al. studied the frequency and amplitude
mances of the suspension, which largely determines the dynamics of a commercial bushing with carbon black
handling stability of the vehicle. During the movement
of the vehicle, the suspension and the suspension bush-
ing play an important role. The deformation of the Faculty of Transportation Engineering, Kunming University of Science and
bushing will affect the wheel alignment angle of the Technology, Kunming, China
vehicle’s suspension, and will also affect the force of
the vehicle’s tires, and ultimately affect the vehicle’s Corresponding author:
Jin Gao, Faculty of Transportation Engineering, Kunming University of
dynamic characteristics. It is of great significance to Science and Technology, No. 727 South Jingming Road, Chenggong
carry out research on suspension and bushing, espe- District, Kunming 650500, China.
cially the optimization of the stiffness of the suspension Email: [email protected]

Creative Commons CC BY: This article is distributed under the terms of the Creative Commons Attribution 4.0 License
(https://creativecommons.org/licenses/by/4.0/) which permits any use, reproduction and distribution of the work
without further permission provided the original work is attributed as specified on the SAGE and Open Access pages
(https://us.sagepub.com/en-us/nam/open-access-at-sage).
2 Advances in Mechanical Engineering

filled rubber bushings in the axial and radial directions. model. The model can achieve a reasonable compro-
Based on the measured data, the axial and radial mise between prediction accuracy, identification diffi-
dynamic stiffness models of the rubber bushing are culty and calculation amount, and is an ideal tool for
established. The accuracy of the new engineering model vehicle chassis dynamics simulation and analysis.6
describing the axial and radial dynamic characteristics Sancibrian et al. use a spatial synthesis method to
of carbon black filled rubber bushing is verified by design the kinematics of the vehicle’s double wishbone
experiments. This work allows the bushing characteris- suspension system, using a multi-angle synthesis tech-
tics in a complex system to be calculated quickly and nology to establish the RSSR-SS spatial link, and
easily.1 Gibanica proposed a new method of calculating express it through a constraint equation. Numerical
the stiffness of the bushing. In the experiment, four examples show that this method can effectively improve
bushings were used for comparative analysis to quantify the design of double wishbone suspension system and
the uncertainty of the bushing parameters, and an has good robustness. This study can achieve high accu-
uncertainty quantization method based on bootstrap- racy in the definition of functional parameters during
ping linear alternative model within parameters was suspension movement, and provides a better reference
proposed. For three identical rear subframes, finally a for suspension design.7 Balike et al. proposed a compre-
model update program using damping equalization has hensive power quarter vehicle model to study the kine-
been successfully used to update all bushes of the rear matics and dynamics of the link suspension and the
subframe and find the optimal parameters of the influence of link geometry on selected performance
bushes.2 Plagge et al. proposed a model that is based on measures. Combined with the kinematics of the double
physical ideas and reasonable assumptions about the wishbone suspension, a plane two-degree-of-freedom
microstructure of materials and is designed for high model was established.8
efficiency and robustness in finite element applications.3 The above experts combined different bushing struc-
Sedlaczek et al., proposed an advanced phenomenologi- tures and materials to put forward an equivalent model
cal Daimler model of rubber bushings. The modeling of bushing stiffness, and analyzed the characteristics of
complexity of this model can be easily adapted to the the bushing. From different perspective, use the differ-
specific problem requirements of various simulation ent methods to study vehicle suspension. However, the
tasks. Examples of two typical chassis bushings under above experts still have some problems that have not
durable load conditions demonstrate the practical been solved. First of all, experts did not analysis the
applicability and superior performance of the modular relationship between bushing deformation and force in
and adaptable Daimler model. The modular approach detail, and did not calculate the amount of bushing
allows the basic building blocks of a specific combina- deformation and force during suspension movement.
tion of problems, representing stiffness, static friction Secondly, the experts did not analysis the impact of
and frequency-dependent dynamics. They can be easily suspension kinematics characteristics on the transient
parameterized using Excel’s graphical preprocessing, and steady-state characteristics of the vehicle, and the
and the achievable model accuracy can be directly impact of the suspension compliance characteristics
checked and adjusted to current requirements. In order caused by the deformation of the suspension bush on
to overcome the problem of low accuracy due to the the transient and steady-state characteristics of the
limitation of expressing the frequency-dependent vehicle. The experts did not conduct sensitivity analysis
dynamic characteristics of the vibration isolator.4 on the stiffness of the suspension bushings, and did not
Aydemir and Sendur5 developed and proposed a simpli- find the bushes that have a greater impact on the tran-
fied transfer function model that varies with frequency, sient characteristics of the vehicle. The third is that the
by generating complex stiffness and damping varying experts did not match and optimize the bushing stiff-
with frequency from the horizontal measurement of the ness to improve the transient performance of the
vehicle. The proposed method has been verified on a vehicle.
heavy commercial truck, and the accuracy of the pro- Aiming at the problems that experts have not solved
posed model is better than that of Voigt modeling yet, this article will start our research from three
method. Zenglian et al. first derives the mathematical aspects. The first will analysis the impact of suspension
equations of the Berg model and Dzierzek model, and movement on the wheel alignment angle and calculate
discusses the theoretical methods of parameter identifi- the value of the wheel alignment angle, and analysis the
cation. Secondly, the static and dynamic characteristics displacement and force of the specific suspension, espe-
of the rubber sleeve were tested, and the parameters of cially the displacement and force analysis of the multi-
the two models were determined based on the experi- link suspension. In addition, the deformation and force
mental data, and the prediction accuracy of the two of the bushing during suspension motion will be calcu-
models was analyzed. On the basis of detailed analysis lated, and the influence of bushing deformation on the
of dynamic stiffness and damping, a new empirical compliance of the suspension will be analyzed in detail.
model is established using Berg model and Dzierzek The second aspect will analysis the impact of the
Gao and Wu 3

Figure 2. Schematic diagram of rubber bushing’s structure.

Figure 1. Multi-body dynamic model of vehicle’s suspension:

(a) front double wishbone suspension and (b) rear multi-link
suspension.

bushing on the transient characteristics of the vehicle,

focusing on explaining the influence mechanism of the
bushing on the transient and steady-state characteristics
of the vehicle. Suspension bushings will be analyzed for
sensitivity, and suspension bushings that have a greater
impact on the vehicle’s transient conditions will be
selected. In the third aspect, the selected bushing stiff-
ness is used as the optimization variable, and the fre- Figure 3. Schematic diagram of bushing mechanical model.
quency characteristic index of the vehicle under the sine
swept steering input is used as the optimization target diagram of the cylindrical rubber bushing is shown in
to match and optimize the suspension bushing, thereby Figure 2.
improving the dynamic characteristics of the vehicle. The rubber bushing is fixed on the metal inner sleeve,
and a radial (y or z axis direction) load W is applied to
the metal outer sleeve to cause the outer sleeve produc-
Establish vehicle suspension model and ing a displacement S. In theory, the dimensionless
calculation of rubber bushing’s stiffness radial stiffness coefficient n is often used to characterize
the ability of the bush to resist radial deformation. The
In order to analysis the kinematics characteristics and
mechanical model of the rubber bushing is shown in
compliance characteristics of the vehicle suspension and
Figure 3 and the stiffness curve of the bush is shown in
the influence of the suspension bushings on the suspen-
Figure 4.
sion performance, the multi-body dynamic model of the
The radial stiffness coefficient of the bushing is
vehicle suspension was established according to a real
defined as:
vehicle,9 as shown in Figure 1.
The rubber bushing is a rubber metal piece com- W
posed of metal and rubber, the outer cylinder and the n= ð1Þ
GSL
inner cylinder are composed of metal, and the rubber
material is in the middle, the rubber bushing is inte- Among them, G is the shear modulus of rubber. The
grated after vulcanization. In order to accurately analy- radial stiffness is:
sis the dynamic characteristics of the vehicle suspension
system, the mechanical properties of the rubber compo- W
Kty = = nGL ð2Þ
nents need to be analyzed accurately. The structure S
4 Advances in Mechanical Engineering

Figure 4. Diagram of rubber bush’s stiffness.

The approximate calculation formula of radial stiffness

coefficient is:
p
n = h   r 2 r2 i ð3Þ
1 r2 2
4 ln r1  r2 2 + r2 +
1   12 
1
L
10 c1 + c2
r1 Figure 5. Schematic diagram of double wishbone suspension.

among them, influence on the wheel alignment angle. Controlling the

8 change of the wheel alignment angle within a reason-
> 1 able range can maintain good dynamic characteristics
>
> c1 = 60s1
>
>
>
>
>
> s2
>
> c2 = 2
>
> s1
>
<    
r2 2 ln2 rr21 r2 2 ln rr21 3ðr2 2  r1 2 Þ ð4Þ
>
> s1 = + 
>
> 2 2 2 2 2
>
> 2ðr2  r1 Þ 4 ðr2 + r1 Þ 16r1
>
>      
>
> 3 r 2 r
>
> r2 4 ln r21 r2 4 ln r21 11r2 2 ln rr21 ðr2 2  r1 2 Þð7r2 2  2r1 2 r2 2 + 7r1 4 Þ
>
> s2 = +  +
:
4ðr2 2  r1 2 Þ2 8ðr2 4  r1 4 Þ 96r1 2 576r1 4 ðr2 2 + r1 2 Þ

When the rubber bush is subjected to an axial (x axial

direction) load, the elastic theory can be used to derive
the axial stiffness calculation formula:
 2
2pLG 1  ln rr21
Ktx =        ð5Þ
3 1 1 r2 r2 r2 2 1 r2
2  ln  1 + ln + ln2  1 + ln2
r1 r1 r2 r1 r1 r1 r2 r2 r1

The calculation of torsional stiffness (around the x-axis) is:

of the vehicle during the movement. After reasonable
 1
1 1 simplification of the suspension, the wheel alignment
Krx = 4pLG 2
 2 ð6Þ angle of the suspension can be obtained through geo-
r1 r2
metric relation and vector calculation of suspension.

Kinematics analysis of vehicle’s suspension Analysis and calculation of vector ~

ek
The influence from suspension’s movement on the vehi- In order to calculate the wheel alignment angle of the
cle’s handling stability is specifically reflected in the ek and ~
suspension, the vectors ~ ehub in Figure 5 need to
Gao and Wu 5

be calculated.10 In order to calculate the vector ~ ek , the The position vector of the wheel knuckle can be defined
*
upper control arm, lower control arm and wheel hub of *
as: BD = lk ek , the unit vector of the wheel knuckle is
the suspension need to be taken into consideration.
defined as:
Because the mechanical structure of the double wish-
bone suspension can be seen as consisting of RSSR- 2 3
d + lu cos x cos j + lu sin x sin F sin j  ll cos u
SSP (R: revolute, S: spherical, P: prismatic) links,11 the * 14 5
ek = h + lu cos x sin j  lu sin x sin F cos j
suspension in Figure 5 can perform two independent lk
q + lu sin xcosF  ll sin u
complete loop analysis of OBDAO and OBNHJO.
ð12Þ
First consider the RSSR closed-loop connection rod
composed of OABCO. The ball hinge position at J is The unit vector of the lower control arm can be defined
considered to be fixed in a specific position. In order to as:
analysis the RSSR connecting rod, the required struc-
tural parameters are: ll , lu , lk , j, and F. When these para- *
e1 = ½cos u 0 sin uT ð13Þ
meters are determined and given a set of lower control
arm angle u, the RSSR connection rod can be used for By formulas (10)–(13), a point position can be obtained
kinematic analysis. The vector equation of this loop is: when the u value of a set of lower control arm is given.
 *  * *
  
BD = OD  OB. The vector on the right side of the
equation Analysis and calculation of vector ~
ehub
 * can be written as:
  * *
In order to calculate vector ~ehub , now discuss the SSP
OB = ll cos u i + ll sin uk
 ! connection rod system. Due to the needs of analyzing,
  *
 OD  = ðd + lu cos x cos j + lu sin x sin F sin jÞ i the OBNHJO connecting rod system was considered.
ð7Þ
* For this analysis, the required parameters are: the ini-
+ ðh + lu cos x sin j  lu sin x sin F cos jÞ j tial coordinates of the spherical hinge points H and J.
 *
*
 
+ ðq + lu sin xcosFÞk The length of the steering rod is tr = HJ , the effective
 * *  *
Usually the displacement equation of the RSSR con- 
steering arm length is sa = HB 3 HD=DB. The
necting rod is:
change in toe angle of the double wishbone suspension
due to the input angle u can be determined by con-
k1 + ðk2  k3 cos uÞ cos x  ðk4 + k5 sin u  k6 cos uÞ sin x
straining OBNHJO. The first constraint is that for each
 k7 cos u  k8 sin u = 0 position of the mechanism, the length of the tie rod
ð8Þ remains unchanged, so there is:
among them,
tr 2 = ðHx  Jx0 Þ2 + ðHv  Jv0 Þ2 + ðHz  Jz0 Þ2 ð14Þ
8
>
> k1 = h2 + d 2 + q2 + ll2 + lu2  lk2 The second constraint is that the length of the steering
>
>
>
> k2 = 2h lu sin j + 2dlu cos j arm remains constant for each position of the mechan-
>
>
>
> k = 2ll lu cos j ism, so there is:
< 3
k4 = 2h lu cos j sin F  2qlu sin j sin F
ð9Þ  2
>
> k5 = 2ll lu cos F sa 2 = ðNx  Hx Þ2 + Ny  Hy + ðNz  Hz Þ2 ð15Þ
>
>
>
> k 6 = 2ll lu sin j sin F
>
>
>
> k = 2dll among them,
: 7
k8 = 2qll rffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
ffi

For a given set of lower control arm angle u and com-
* * *  *2
ON = ll e1 + l a ek , la = BH   s2a ð16Þ
bine formula (8), use the half-tangent formula can get
the value of x: The third constraint is that for each position of the
pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi! mechanism, the steering arm vector and the steering
1 b 6 b2  4ac knuckle vector are perpendicular to each other
x = 2tan ð10Þ  * *
2a HN ? BN . Because the scalar product of two vertical
* *
among them, vectors is equal to zero, so: HN  BN = 0.
8  *  * *  
< a = k1 + k2  k3 cos u + k7 cos u + k8 sin u * * *
l l e 1 + l a ek  H x i + H y j + H z k  l a ek = 0
b = 2k4 + 2k5 sin u + 2k6 cos u ð11Þ
:
c = k1  k2  k3 cos u + k7 cos u + k8 sin u ð17Þ
6 Advances in Mechanical Engineering

Figure 6. Vector transformation: (a) schematic diagram of wheel static surface position vector, (b) schematic diagram of Angle
change of vector ek, and (c) schematic diagram of the rotation of wheel knuckle relative to ek.
*
In equation (17), it can get: is not considered. Since the unit matrix ek of any posi-
tion of the mechanism can be calculated by formula
*
c3 H x + c4 H y + c5 H z + c6 = 0 ð18Þ (13), the rotation amount of ek can be determined by
the outer product as:
among them,
h * *
 *i
* * * l = sin1 ek0 3 ek  ekl ð22Þ
c3 = la ekx , c4 = , c5 =
la eky ,
la ekz
* * * * * *
 ð19Þ *
c6 = la 2 + ll la ekx e1x + eky e1y + ekz e1z ekl is the vector on the rotation axis perpendicular to
* *
the plane formed by vectors ek and ek0 ,as shown in
In equations (17)–(19), the unknown parameters Hx , Figure 6(b):
Hy , and Hz are the new coordinates of the steering rod *  * 
* *  *
ball hinge point H. These three equations can solve ekl = ek0 3 ek =ek0 3 ek  ð23Þ
three unknowns, at this stage, the positions of all joints
are determined. In the static position of the suspension The second rotation analysis requires a complete analy-
system, the unit vector GW of the hub can be deter- sis of the rotation of wheel knuckle with respect to the
*
mined as: unit vector ek , as shown in Figure 6(c):
* h *  *i
ehub0 = ½cos e0 0  sin e0 T ð20Þ lp = sin1
*
ea0l 3 ea  ek ð24Þ
e0 is the position of the positive static surface of the
The matrix with a rotation angle of u around a certain
wheel, as shown in Figure 6(a). The static position of *
axis in the direction of the unit vector u is:
2 3
*  cos u + u2x ð1  cos uÞ ux uy ð1  cos uÞ  uz sin u ux uz ð1  cos uÞ + uy sin u
6 u u ð1  cos uÞ + u sin u cos u + u2y ð1  cos uÞ uy uz ð1  cos uÞ  ux sin u 7
R u ,u =4 y x z 5 ð25Þ
uz ux ð1  cos uÞ  uy sin u uz uy ð1  cos uÞ + ux sin u cos u + u2z ð1  cos uÞ

Because of the related rotation matrix has obtained,

*
the wheel center can be determined: the unit vector ehub from the initial position to the new
position can be determined as:
* *
lhub0 = c0 ½0 1 0T + lhub ehub0 ð21Þ *
*  *  *
ehub = R ek , lp R ekl , l ehub0 ð26Þ
In order to get the any position of the vector on the
wheel knuckle, the rotation matrix can be used and take
the following steps: starting from the initial position of
the steering arm, determine the rotation matrix of the
Calculation of wheel alignment angle
*
unit vector of the steering arm from ea0 to any position In order to analysis and calculate the wheel alignment
*
ea . Among them, the rotation of the wheel knuckle is angle of the double wishbone suspension, the double
considered as the superposition of two different rota- wishbone suspension is simplified to obtain a schematic
tions. Initially, the rotation of the joint around its axis diagram as shown in Figure 7. The camber angle, caster
Gao and Wu 7

Figure 8. Variation diagram of toe angle with wheel travel.

Figure 7. Schematic diagram of wheel alignment angle of
double wishbone suspension.
* *

s = cos1 ek  ehubp  p=2 ð31Þ
angle, kingpin angle, and toe angle in Figure 7 can all
be calculated by vectors.12 among them,
In Figure 7, ~k represents the unit vector in the z-axis * T .qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
ffi
direction, ~j represents the unit vector in the y-axis direc- ehubp = ehubpx ehubpy 0 e2hubpx + e2hubpy ð32Þ
ek represents the unit vector of the wheel hub, ~
tion, ~ ehub
represents the unit vector of the wheel axis, ~ ehubp repre- After analyzing the influence of the suspension motion
sents the unit vector obtained by projecting the vector on the wheel alignment angle, the opposite travel simu-
ehub in the XY plane, ~
~ ecaster is the unit vector perpendi- lation on the front suspension and the rear suspension
cular to vector ~ ehubp in the XY plane, a is the toe angle was carried out, the changes of toe angle and camber
of the wheel, u is the camber angle, t is caster angle, s angle with wheel travel as shown in Figures 8 and 9
is kingpin angle. The toe angle of the suspension system were obtained.
can be calculated from the vectors ~ ehubp and ~j: In Figure 8, during the wheel travel downward, the
toe angle of the front suspension wheel (the inner wheel
 
a = p=2  cos1 e ! ~
hubp  j ð27Þ when turning) increases in the positive direction, and
the toe angle of the rear suspension wheel (the inner
The camber angle shown in Figure 7 can be calculated wheel when turning) increases in the negative direction.
ehub and ~
by vectors ~ k: The toe angle of the corresponding outer wheel of the
 * * front suspension increases negatively, and the toe angle
u = cos1 ehub k  p=2 ð28Þ of the outer wheel of the rear suspension increases posi-
tively. For the front suspension, the toe angle of the
The caster angle shown in Figure 7 can be calculated by inner wheel increases positively, and the toe angle of
ek and ~
vectors ~ ecaster : the outer wheel increases negatively, which will increase
*  the understeer characteristics of the vehicle. For the
*
t = cos1 ek  ecaster  p=2 ð29Þ rear suspension, the toe angle of the inner wheel
increases negatively, and the toe angle of the outer
among them, wheel increases positively, which will also increase the
understeer characteristics of the vehicle.
ecaster * = In Figure 9, when the wheel travels downward, the
2 3T camber angle of the front and rear suspension wheel
vffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi vffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
6 ehuby u u 1 u 1 7 (the inner wheel when turning) increases in the positive
47 t 2 6u
t 2 05 direction. When the wheel travels upward, the camber
ehubx ehuby ehuby
ehubx +1 ehubx +1 angle of the wheel (the outer wheel when turning) of the
ð30Þ front suspension and the rear suspension both increases
in the negative camber direction. In terms of adhesion,
The kingpin angle of the suspension system can be cal- the front suspension’s and the rear suspension’s outer
culated by vector ~
ek and vector ~
ehubp : wheel camber angle increase in the negative camber
direction, the tire’s adhesion will be improved. For the
8 Advances in Mechanical Engineering

Figure 9. Variation diagram of camber angle with wheel travel.

perspective of understeer characteristics of the vehicle,

when the vehicle is turning, the camber angle of the
front suspension’s outer wheel increases to the negative
camber direction, the understeer characteristic of the
vehicle will decrease, and the camber angle of the rear
Figure 10. Schematic diagram of double wishbone suspension
suspension’s outer wheel increases to the negative cam- system.
ber direction, the understeer characteristics of the vehi-
cle will increase.
The two ends of the coil spring and the shock absorber
Statics analysis of double wishbone are installed between the lower control arm and the
vehicle’s body, which are installed between the ball
suspension joints C and E, respectively.
When the wheels of the double wishbone suspension If an object rotates from one position to another
are excited by the road and jump, the hinge point of without displacement, the coordinates of the object in
the suspension link and the link rod will move with the the global coordinate system o  xyz can be determined
travel of the wheel, so the position of each hinge point by the coordinate system fixed on the rigid body. Take
will change accordingly. Through the theory of coordi- the upper control arms A1D and A2D of the double
nate transformation, after selecting a reference point, wishbone suspension as an example, assuming that the
the coordinates of other hinge points can be calculated rotation angles around the x, y, and z axisare h1, r1,
0 0
accordingly. The force balance and torque balance and z1 respectively, the new coordinates A1 , A2
0
equations can be used to calculate the force of the vehi- and D of points A1, A2 and D on the rigid body
cle suspension hinge points and bushing. A1D and A2D in the coordinate system fixed to itself
can be obtained.
n o  
Displacement analysis of double wishbone 0
A1 = ½X ðh1 Þ½Y ðr1 Þ½Z ðz1 Þ A01 ð33Þ
suspension’s connecting point n 0o  
The schematic diagram of the double wishbone suspen- A2 = ½X ðh1 Þ½Y ðr1 Þ½Z ðz1 Þ A02 ð34Þ
sion is shown in Figure 10. In order to analysis the force n 0o  
and moment of suspension’s ball joint and bushing, D = ½X ðh1 Þ½Y ðr1 Þ½Z ðz1 Þ D0 ð35Þ
firstly, the position transformation of each connection
point of the double wishbone suspension is analyzed.13 ½X ðh1 Þ, ½Y ðr1 Þ, and ½Z ðz1 Þ are the ð3 3 3Þ rotation
The double wishbone suspension is composed of upper matrices of the upper control arm around the x, y, and
control arm, lower control arm, wheel knuckle, springs, z axis on their own coordinate  system
  andtherotation
shock absorbers and tie rods. The four bushes A1 and angles are h1 , r1 , and z1 . A01 , A02 , and D0 are the
A2, A3 and A4 in Figure 10 connect the upper and lower initial coordinates of points A1 , A2 , and D.
control arms to the vehicle body respectively. The ball When the rigid body rotates and moves at the same
joint D and the ball joint G respectively connect the time, select point
 0  D as the
 0 reference
 point, then the new
wheel knuckle with the upper and lower control arms. coordinates A1 and A2 of points A1 and A2 in the
Gao and Wu 9

global coordinate system can be obtained by formula H 0 and P01 are the initial coordinates of points H and
(33) subtracting formula (35) and formula (34) subtract- P1 , and the matrix ½Rðh3 , r3 , z3 Þ can be calculated as
ing formula (35): follows:
( 0 )  
A1 ½Rðh1 , r1 , z1 Þ, 0 ½Rðh3 , r3 , z3 Þ = ½X ðh3 Þ½Y ðr3 Þ½Z ðz3 Þ ð41Þ
 0  =
A2 0, ½Rðh1 , r1 , z1 Þ ½X ðh3 Þ, ½Y ðr3 Þ, and ½Z ðz3 Þ are the 3 3 3 rotation
(     ) (  0  ) ð36Þ
A 0  D0 D matrices of the wheel knuckle around the x, y, and z
 10   0  +  0  axis on its own coordinate system and the rotation
A2  D D angles are h3 , r3 , and z3 . According to equations (36),
(38), (40), when external force is applied to the suspen-
The
 0 0 in
 formula
 (36)represents a zero matrix of 3 3 3, sion, the new position of each connection point and
A1 , A02 and D0 are the initial coordinates of the bushing can be obtained. The displacement of the bush-
connection
 0 points
 A1 , A2 and D in the global coordinate ing and spring can also be easily estimated.
system; D is the new coordinate value of hinge point
D, matrix ½Rðh1 , r1 , z1 Þ is a rotation matrix:
Calculate the forces and moments acting on the
½Rðh1 , r1 , z1 Þ = ½X ðh1 Þ½Y ðr1 Þ½Z ðz1 Þ ð37Þ bushing
 g 
In the same way, taking point If a static external force FP1 is applied to the contact
n o  0G as thereference
0
point,
the new coordinates A3’ , A4 , and C of the con- point between the tire and the ground, the bushing will
nection points A3 , A4 , and C on the lower control arm rotate and translate. The displacement of the bushing S
can be obtained: (S means A1 , A2 ,A3 , and A4 ) can be calculated:
8 0 9 2 3   0    g 
  S  fS0 g DS
< A3  >
> = ½Rðh2 , r2 , z2 Þ, 0 0 QgS = +  0g  ð42Þ
0 6 7 ffS g Df S0
A4 = 4 0, ½Rðh2 , r2 , z2 Þ, 0 5
:  0 >
> ;
C 0 0½Rðh2 , r2 , z2 Þ Superscript
 0 g means in the global coordinate system.
8  0  0 9 8  0  9 ð38Þ S represents the translational position of the bush
< A3   G  >
> = > < G  > = after deformation. ffS g is the rotation angle of the
0 0 0
A4  G + G bush S in the global
:  0  0 >
> ; > :  0 > ;  coordinate
  system when the exter-
C  G G nal force is input. DS0g and DfgS0 in equation (42)
are the initial displacement and initial angle of the bush-
   
A03 , A04 , fC0 g, and fG0 g are the ing in the global coordinate system under the action of
 initial
0
positions of
the hinge points A3 , A4 , C, and G. G is the new coor- the vehicle’s own weight. In the global coordinate sys-
dinate of the hinge point G; Matrix ½Rðh2 , r2 , z2 Þ is a tem, the force and moment applied to the bush due to
rotation matrix: the displacement of bush S are:
  g   g  
½Rðh1 , r1 , z1 Þ = ½X ðh2 Þ½Y ðr2 Þ½Z ðz2 Þ ð39Þ FSg  = K St g, 0 fQg g ð43Þ
MS 0, KSr S

½X ðh2 Þ, ½Y ðr2 Þ, and ½Z ðz2 Þ are the 3 3 3 rotation  g  

matrices of the lower control arm around the x, y, and FS and MSg are the force matrix and moment
z axis on its own coordinate system and the rotation matrix of 3 3 1, respectively.
angles are h2 , r2 , z2 . In the same way, with point D as (  
0  0
the reference
 0 point, the new coordinates G , H , KStg = ½RS T KStl ½RS 
g l  ð44Þ
and P1 of the hinge points G, H, and P1 can be KSr = ½RS T KSr ½R S 
obtained: 
8  0 9 2 3 T represents the transpose of the matrix. KStg and
g
< G  >
> = ½Rðh3 , r3 , z3 Þ, 0, 0 KSr are the translational movement stiffness matrix
0 6 7 and the rotation stiffness matrix respectively of the rub-
H = 4 0, ½Rðh3 , r3 , z3 Þ, 0 5
>
: 0 ; > ber
P1 0, 0, ½Rðh3 , r3 , z3 Þ  bushing  S in the  global
 coordinate
 system.
8  0   0  9 8  0  9 ð40Þ KSt =diag kStx ,kSty ,kStz and KSr =diag kSrx ,kSry ,kSrz
< G   D  >
> = > < D  > = are the translational stiffness matrix and rotational
0 0 0
H  D + D stiffness matrix of the rubber bushing S in the local
>
:  0   0  > ; > :  0 > ;
P1  D D coordinate system. ½RS  is the 333 transformation
matrix of the bushing S from the local coordinate sys-
tem to the global coordinate system. In the global
10 Advances in Mechanical Engineering

 
coordinate system, the force and moment applied to FGg is the reaction force appliedto the
 lower control
the upper control arm by the bush Su (Su refers to the arm from the wheel knuckle G. FCg is the reaction
bush A1 and the bush A2) are: force of the spring acting on the lower control arm.
8     Under the action of the external load of the upper and
> fRFSgu g I,
 0 fFSgu g lower control arms, tie rods and the contact point
>
< fRM g g = ~
R ,I fMSgu g between the tire and ground, the wheel knuckle is in a
 Su   Su D   ð45Þ
>
>
g
fRFSl g Ih i0 fFSgl g balanced state,14 and the static balance equation is
: g = R ~ g
fRMSl g Sl G , I fMSl g obtained:
  g    g   g    g
I is the unit matrix of ð3 3 3Þ, R ~ Su D is an oblique sym-
FP1g   FDg   FHg   FG
= f0g
metric matrix of the position vector of the bush Su in MP1 MD MH 0
the global coordinate
 system, and point D is the refer- ð49Þ
ence point. R ~ Sl G is an oblique symmetric matrix of
 g 
the position vector of the bush Sl in the global coordi-
 FP1 is gthe
 external force input
 g to
 the tire from the
nate system. The oblique symmetric matrix ~ Su D and
R road, MP1 isthemoment of FP1 about the connec-

~ Sl G are:
R tion point D. FHg is the reaction force acting on the
8 2 3 wheel knuckle from the tie rod, expressed as:
> 0  zgSu D ygSu D   
>
>
>
>  6 g 7 FHg = ns1 RJH ð50Þ
>
>
>
~ Su D = 6 z S D
R 4 u 0  xgSu D 7
5
>
>   
< yg xgSu D ns1 is the scale
2 Su D
0
3  factor,
 RJH is the vector  from
 point H
ð46Þ to point J . MHg is the moment of force FHg relative to
>
> 0  zgSl D ygSl D
>
>  6 7 point G. The length, displacement and initial length of the
>
> ~ Sl D = 6 zgS D
>
> R 4 l 0  xgSl D 7
5 spring after deformation satisfy the following formula:
>
>
:
ygSl D xgSl D 0 0
LEC + DLEC  LEC = 0 ð51Þ
xgSu D , ygSu D , zgSu D are the displacement of the bush Su along 0
the x, y, and z directions in the global coordinate sys- LEC is the new length of the spring after deformation,
tem. xgSl D , ygSl D , zgSl D are the displacement of the bush Sl LEC is the initial length of the spring, DLEC is the defor-
along the x, y, and z directions in the global coordinate mation of the spring:
system.  0
.
The upper control arm is in a balanced state under DLEC = ns2 LEC  ns20 LEC Kd ð52Þ
the action of the reaction force and the reaction
moment of the bushing, the spring, and the ball joint of ns20 is the initial scale factor of the spring force, Kd is
the wheel knuckle. Therefore, the sum of each force the stiffness of the spring. Due to the length of the tie
component of all forces, as well as the sum of each rod does not change during the movement of the sus-
force and each moment component of the connection pension, there is:
point D in the global coordinate system is zero.
0
According to this principle, the static balance equation LJH  LJH = 0 ð53Þ
of the upper control arm can be obtained and expressed 0
as: LJH and LJH are the length of the tie rod after and
before the suspension change, respectively. Through
 g    g    g
the above analysis, the overall balance equation can be
RFA1g  + RFA2g  + FD
= f0g ð47Þ
RMA1 RMA2 0 obtained as:
 g 8 g    g    g 9
FD is the combined force applied to the upper control > RFA1 RFA2 F >
> RM g  + RM g  + 0 D
>
>
>
>
>
arm from the wheel knuckle and the connection point >
> A1 A2 >
>
>
>                 >
>
>
> RF g
RF g
F g
F g >
>
D. Similar to obtaining the static balance equation of >
>  A3  +  A4  +  C  + G >
>
< RMA3g
RMA4g
MC g
0 =
the upper control arm, the static balance equation of  g    g   g    g = f0g
> >
the lower control arm can be obtained, expressed as: > FP1g   FDg   FHg  
>
>
FG >
>
>
>
> M MD MH 0 >
>
>
>  0 P1 >
>
( g 
) ( g  ) >
>
>
 >
>
>
RFA3 RF >
: L EC + DL EC  L EC >
;
+  A4g 
0
 g 
LJH  LJH
RMA3 RMA4
( g ) ( ) ð48Þ ð54Þ
F FGg
+  Cg  + = f0g Take the lateral force in the interval ½3000, 3000, take
MC 0 a value into the formula (54), the equation (54) can be
Gao and Wu 11

Figure 11. The force in the x direction of the rubber bushing.

Figure 13. Schematic diagram of multi-link suspension.

For the wheel knuckle, when point B5 isselected as

0
the reference point, the new position Bi for other
points Bi (i = 5;8) can be calculated:
n 0o     n 0 o
Bi = ½X ðh4 Þ½Y ðr4 Þ½Z ðz4 Þ B0i  P02 + P2

Figure 12. The moment of the rubber bushing around the x- ð55Þ
axis.  0
  0

P2 and Bi are the initial coordinates of points P2
0 
and fBi g respectively; P2 is the new coordinate of
point P2 ;½X ðh4 Þ, ½Y ðr4 Þ, and ½Z ðz4 Þ are the ð3 3 3Þ
solved by using the Newton iteration method, the force
rotation matrices when the wheel knuckle rotates
and moment situations shown in Figures 11 and 12 can
around the x, y, and z axis respectively, and the rotation
be obtained.
angles are h4 , r4 , and z4 .
For connection rod Ai Bi (i = 5;8), when point Ai is
Statics analysis of multi-link suspension selected as the reference point, the new coordinates of
point Bi in the global coordinates system is:
When an external load acts on the contact point n 0o     n 0 o
between the tire and the ground, the wheels, wheel Bi = ½X ðhi Þ½Y ðri Þ½Z ðzi Þ B0i  A0i + Ai
knuckle and connection rods will move to new position
due to the elastic deformation of the bushes and springs ð56Þ
in the suspension. In this way, a new static equilibrium    
Where A0i and B0i are  the initial coordinates of
is reached. The new position of each bush in this load 0
points Ai and Bi respectively; Ai is the new position
can then be calculated.
of point Ai ;½X ðhi Þ, ½Y ðri Þ, and ½Z ðzi Þ are the ð3 3 3Þ
rotation matrix when connecting rod Ai Bi rotates
Displacement analysis of multi-link suspension’s around x, y, and z axis respectively, and the rotation
angles are hi , ri , and zi .
connection point According to equations (55) and (56), when external
The schematic diagram of the multi-link suspension is force acts on the suspension, the new position and dis-
shown in Figure 13. placement of each bush can be obtained.
12 Advances in Mechanical Engineering

Calculate the forces and moments acting on the knuckle is in a balanced state. Therefore, the sum of
rubber bushing each reaction force (RF), reaction moment (RM) and
external force applied to the wheel knuckle is zero in
When an external force is applied to the wheel knuckle,
the global coordinate system. The static balance equa-
the bushing will rotate and translate.15 For the bushing
tion of the wheel knuckle is expressed as:
w (w means A5 , A6 ,A7 , and A8 ) connectingthe connec-
tion rod and the body, the displacement Qgw in the 8 9
8 n
P o
>
> g >
>
global coordinate system is:   >
< RFBj >
=
 g 
g
FB0 j=5
  0    g  FP2 + g
MB0
+
> 8 n
P o > = f0g ð61Þ
  w  fw0 g Dw >
> g >
>
g
Qw = +  0g  ð57Þ : RMBj ;
ffw g Dfw0 j=5
 g 
The superscript g in equation (57) means in the global FP2is the external force applied to the wheel at point
coordinate system; ffw g is the rotation angle matrix
 gof P2 .The rod is in equilibrium under the reaction force
the bushing
 g  w in the global coordinate system; Dw0 and moment of bushing Ai and hinge point Bi , there-
and Dfw0 are the initial displacement and rotation fore, in the global coordinate system, the sum of each
angle of bushing w in the global coordinate  gsystem
 force generated by the elastic deformation of the bush
under the weight
 of the vehicle. The force Fw and and each moment generated by the force relative to
moment Mwg due to the displacement of the bushing point Ai is zero. Therefore, the static balance equation
are: of rod Ai Bi (i = 5–8) can be obtained and expressed as:
  g   g    g    g  
 
Fwg  = ½Kwt , 0g  Qgw ð58Þ RFwg  + RFBig  = f0g ð62Þ
Mw 0, Kwr RMw RMBi

0 is the zero matrix of ð3 3 3Þ. The initial length, deformed length and displacement of
( the spring satisfy the following formula:
g l 
½Kwt  = ½Rw T Kwt ½R w 
g T l  ð59Þ
Kwr = ½Rw  Kwr ½Rw  LA‘ 0 B0 + DLA0 B0  LA0 B0 = 0 ð63Þ
0
the transpose of the matrix. ½Kwt g
 and LA0 B0 is the length of the spring after deformation; LA0 B0
grepresents
T  is the initial length of the spring; DLA0 B0 is the deforma-
Kwr are the translational movement stiffness matrix
and the rotation stiffness matrix respectively of the tion of the spring:
l  .
bush w in the global coordinate
  system. Kwt = 0
l
diag kwtx , kwty , kwtz and Kwr = diag kwrx , kwry , kwrz are DLA0 B0 = ns3 LA0 B0  ns30 LA0 B0 Ks ð64Þ
the translational stiffness matrix and rotational stiff-
ness matrix of the bush w in the local coordinate sys- ns30 is the initial scale factor of the spring force, Ks is
tem. ½Rw  is the 3 3 3 transformation matrix of the the stiffness of the spring.
bushing w from the local coordinate system to the glo- On the basis of the above analysis, the balance equa-
bal coordinate system. According to Newton’s third tions of the suspension system can be obtained. By
law and equation (58), the reaction force and moment combining the balance equation of the wheel knuckle
applied by the rubber bushing to the control arm in the (equation (61)), the balance equation of the rod (equa-
global coordinate system are: tion (62) and the constraint equation of spring (equa-
  g    g    tion (63)), the equation of the multi-link suspension is
g
RFwg  = Fwg  = ½Kwt , 0g  fQgw g expressed as:
RMw Mw 0,  Kwr 8
8 8 n
P o 99
>
> >
> g >
> >
ð60Þ >
>
>  
 g  >
< RFBj =>
> >
>
>
>
> g F B0 j=5 >
>
Mark the forces and moments on the hinge points B5 , >
> F + g + n o >
>
>
<
P2 M >
> P
8 >
> >
;=
B0 g
B6 , B7 and B8 as fFBi g and fMBi g, then for the rigid >
: RM > Bj
  g     g j =5 = f0g
body A5 B5 , A6 B6 , A7 B7 , and A8 B8 , the reaction force >
> >
>
>
> RFwg  + RFBig  >
>
and moment at the Bi ði = 5;8Þ end are denoted as >
> >
>
>
> n RMw RMBi o >
>
fRFBi g and fRMBi g. >
> >
>
: L0 + DL L ;
A0 B0 A0 B0 A0 B0

Force balance and moment balance analysis ð65Þ

Under the action of the reaction force and moment of Figure 14 shows the change of the force on the rear
the hinge point Bi and the external load, the wheel multi-link suspension bushing with the lateral force on
Gao and Wu 13

Figure 14. The change of the bush’s force in the x-direction Figure 16. The toe angle and camber angle of the front double
with the lateral force of the wheel. wishbone suspension vary with the lateral force of the wheel.

8 and 9. The characteristic of the tire force and the

change of the wheel alignment angle caused by the
deformation of the elastic elements such as the bushing
in the suspension is called the compliance characteristic
of the suspension. The compliance characteristic of the
suspension includes the characteristic of the tires sub-
jected the same direction lateral force, reverse lateral
force, longitudinal force, same direction align torque,
and reverse align torque. The impact of the bushing on
the dynamic characteristic of the vehicle is mainly in
two aspects, one is that due to the deformation of the
bushing, when the tire is subjected to lateral force, long-
itudinal force and align torque, it will cause the wheel
Figure 15. The torque of the bushing around the x-axis to produce additional steering angle, which will affect
changes with the lateral force of the wheel. the steering characteristic of the vehicle. At the same
time, the bushing will affect the force on the vehicle
wheels, that is, it affects the dynamic characteristic of
the wheels. Figure 15 shows the change of the torque the vehicle by influencing the tire force.
on the rear multi-link suspension bushing with the lat-
eral force on the wheels.
The influence of suspension motion and bushing on
the change of toe Angle and camber Angle
Analysis of the influence of suspension From the analysis of the front two parts, it can be seen
kinematics and compliance characteristic that the force of the vehicle wheel will cause the defor-
on vehicle dynamics mation of the suspension bushing, and the deformation
of the bushing will cause the change of the wheel align-
The actual deformation of the rubber bush during the ment angle. Based on the force analysis about front
movement of the suspension is more complicated, double wishbone suspension and the rear multi-link
including the axial movement of the inner ring (metal suspension, the compliance characteristic simulation
inner surface) relative to the outer ring (metal outer sur- analysis of the double wishbone suspension and the
face), the movement of the inner and outer rings in dif- rear multi-link suspension of the vehicle is carried out,
ferent radial directions, and the rotation of the inner the changes of the toe angle and the camber angle of
and outer rings around the axis etc. The suspension the wheel with the lateral force and the longitudinal
kinematics characteristic is the change characteristic of force as shown in Figures 16 to 19 are obtained.
the wheel alignment angles caused by the geometric Figure 16 shows the change of toe angle and camber
movement of the suspension. The kinematics character- angle of the front double wishbone suspension with the
istic of the suspension include wheel bouncing in the lateral force of the wheel. In Figure 16, both the toe
same direction and opposite direction. The kinematics angle and the camber angle increase in the positive
characteristics of the suspension are analyzed in Figures direction as the reverse lateral force increases. The toe
14 Advances in Mechanical Engineering

Figure 17. The toe angle and camber angle of the front double
wishbone suspension vary with the longitudinal force of the Figure 19. Changes in toe angle and camber angle of rear
wheel. multi-link suspension with wheel’s longitudinal force.

Figure 18. Changes in toe angle and camber angle of rear Figure 20. The effect of toe angle on wheel’s side slip angle.
multi-link suspension with wheel’s lateral force.
by the vehicle is a positive lateral force, the toe angle of
angle and camber angle both increase in the negative the wheels changes less, and the camber angle increases
direction as the positive lateral force increases. Figure in the negative direction. Figure 19 shows the change
17 shows the change of toe angle and camber angle of of the toe angle and camber angle of the rear multi-link
the front double wishbone suspension with the longitu- suspension wheels when subjected to longitudinal force.
dinal force of the wheel. In Figure 17, when the longitu- In Figure 19, whether the wheel receives the driving
dinal force received by the wheels is the braking force, force or the braking force, the toe angle changes little.
the toe angle increases in the negative direction with the For the camber angle, when the wheel receives the
increase of the braking force, and the camber angle braking force, the camber angle increases in the nega-
increases in the positive direction with the increase of tive direction, when the wheel receives the driving force,
the braking force. When the longitudinal force received the camber angle of the wheel increases in the positive
by the wheels is the driving force, the toe angle increases direction.
in the positive direction as the driving force increases,
and the camber angle decreases in the negative direction
as the driving force increases. Effect of toe angle change caused by suspension
Figure 18 shows the change of the toe angle and movement and bushing deformation on tire side slip
camber angle of the rear multi-link suspension with the
stiffness
lateral force received by the wheels. In Figure 18, when
the lateral force received is the reverse lateral force, the When there is toe angle of the wheel, the side slip angle
toe angle of the wheels changes less, and the camber of the outer wheel is aa = a + d0 , and the side slip
angle increases in the positive direction as the reverse angle of the inner wheel is aa = a  d0 , there will be a
lateral force increases. When the lateral force received deviation as shown in Figure 20. At this time, the sum
Gao and Wu 15

of the side slip stiffness at point 1 is too large. The rea-

son is: due to the higher wheel load at point 2, the side
slip stiffness increases more than at point 3. The result
is the toe angle reduces the axle’s side slip angle when
driving in a curve. This means that the toe angle of the
front wheels and the toe angle change caused by the
driving force of the front wheels reduce the side slip
angle of the front axle, the toe of the rear wheels and
the change of the toe angle caused by the driving force
of the rear wheels will increase the side slip angle of the
rear axle. As a result, the vehicle tend to oversteer.
The average side slip stiffness of the inner and outer
wheels is denoted as Kf , without considering the toe
angle of the wheel, the side slip stiffness of the inner
wheel of the front wheel is denoted as Kyi , and the side
slip stiffness of the outer wheel of the front wheel is
denoted as Kya . When considering the toe angle of the Figure 21. Effect of camber on understeer/oversteer
wheel, the side slip stiffness of the inner wheel of the characteristic.
0
front wheel is denoted as Kyi , and the side slip stiffness
0
of the outer wheel of the front wheel is denoted as Kya . The effect of suspension motion and the change of
Regardless of toe angle, the total side slip stiffness camber caused by bushing on vehicle’s steering
Kyi + Kya of the inner and outer wheels is greater than
characteristic
the average side slip stiffness Kf , when considering the
toe angle of the wheel, the total side slip stiffness The variation characteristics of wheel’s side slip Angle
0 0
Kyi + Kya of the inner and outer wheels is greater than with relative centripetal acceleration are shown in
the side slip stiffness Kyi + Kya of the inner and outer Figure 21. When the camber is positive, the side slip
wheels without considering the toe angle, that is angle is greater than when there is no camber, and
when the camber is negative, it is smaller than when
0 0
Kyi + Kya Kyi + Kya there is no camber. The understeer/oversteer trend of
\ \Kf ð66Þ the vehicle will be affected by changes in camber as the
2 2
function of centripetal acceleration (or the function of
The toe angle of the wheel itself and the change of the roll or the function of lateral force). When the centripe-
toe angle will affect the tire side slip angle. The change tal acceleration is low, the front wheel’s camber changes
of the wheel’s toe angle caused by the wheel lateral in positive direction/rear wheel no change, or front
force, longitudinal
 force and align torque is expressed wheel’s camber no change/rear wheel’s camber change
as Da Fy , DaðFx Þ, and DaðT Þ, denote the change of in negative direction, the understeer tendency of the
toe angle caused by the suspension roll as DaðkÞ, there- vehicle will increase; otherwise, the understeer will
fore, the toe angle of the vehicle’s wheel can be decrease. When the lateral acceleration is high, the
expressed as: front wheels have negative camber/the rear wheels have
   no camber, there is a slight oversteer of the vehicle;
av = 6 aV 0 + DaV ðkÞ + DaV FyV + DaV ðFxV Þ + DaV ðT Þ while the front wheels without camber/the rear wheels
aH = 6 aH 0 + DaH ðkÞ + DaH FyH + DaH ðFxH Þ + DaH ðT Þ
have a positive camber, there is a strong oversteer char-
ð67Þ acteristic of the vehicle.
The camber of the wheel will have an effect on the side
aV is the toe angle of the front wheels, aH is the toe slip angle of the wheel. The change of camber angle
angle of the rear wheels, aV 0 is the initial toe angle of caused by the lateral force, longitudinal force
the front wheels of the vehicle, and aH0 is the initial toe   and the
align torque of the wheel is expressed as Du Fy , DuðFx Þ,
angle of the rear wheels of the vehicle. and DuðT Þ, and the change of camber angle caused by
For the toe angle, when designing the suspension, it the suspension roll is expressed as DuðkÞ, so the camber
requires to increase the side slip angle of the wheel dur- angle of the vehicle’s wheel can be expressed as:
ing the change of the toe angle of the front axle caused
8  
by the deformation of the bushing, the toe angle of the < uv = 6 uV 0 + DuV ðdVL Þ + DuV ðkÞ + DuV FyV + DuV ðFxV Þ
rear axle wheels requires to reduce the wheel’s side slip :
+ DuV ðT Þ  
angle during the change of the toe angle of the rear axle, uH = 6 uH 0 + DuH ðkÞ + DuH FyH + DuH ðFxH Þ + DuH ðT Þ
so that the vehicle tends to understeer. ð68Þ
16 Advances in Mechanical Engineering

X
uV is the camber angle of the front wheels, uH is the Fy = 2Yf + 2Yr + 2Ycf + 2Ycr ð71Þ
camber angle of the rear wheels, uV 0 is the initial cam-
ber angle of the front wheels of the vehicle, and uH0 is
the initial camber angle of the rear wheels of the Multi-objective optimization of vehicle’s
vehicle. handling stability
When designing the suspension, for the wheel’s cam-
ber angle, when the centripetal acceleration is low, the There are multiple evaluation parameters for the transi-
deformation of the bushing is required to make the ent characteristics of vehicle. In order to explain these
front wheel camber change in positive direction/rear evaluation parameters in the form of formulas, the
wheel no change, or front wheel camber no change/rear three degree of freedom vehicle dynamics model is used.
wheel’s camber changes in negative direction, which will Since there are multiple optimization targets and multi-
increase the understeer tendency of the vehicle. When ple optimization variables in the optimization, there-
the lateral acceleration is high, let the front wheels have fore, optimizing the bushing of the suspension based on
negative camber/the rear wheels have no camber, so the sine swept input to improve the transient dynamics
that the vehicle has a slight oversteer to avoid strong of the vehicle requires multi-objective optimization.
oversteer characteristic. Therefore, on the basis of the three degree of freedom
vehicle’s dynamics, the optimization targets of vehicle’s
transient characteristic are derived, the multi-objective
Lateral force on vehicle’s tires optimization program are established.
During the vehicle’s movement, the lateral force
received by a vehicle’s tire can be calculated by two Optimize objective function
parts, one is the product of the tire’s side slip stiffness
For parameter-oriented vehicle’s dynamics models, the
and the cornering angle, and the other is the product of
three degree of freedom vehicle dynamics model is
the lateral camber thrust coefficient and the camber
mostly used at present, because the three degrees of
angle. The side slip angle and camber angle of the tire
freedom vehicle’s dynamic model is simpler than other
are obtained by the suspension movement and suspen-
dynamic models, but it can describe the main motion
sion bushing deformation. The force obtained by the
characteristic and force characteristic of the vehicle.
product of the tire side slip stiffness and the tire side
The vehicle’s three degree of freedom model considers
slip angle is:
the vehicle’s lateral, yaw and roll these three degrees of
 freedom. For the transient characteristic of the vehicle,
Yf = Yf 1 + Yf 2 = kfL afL + kfR afR
ð69Þ the three degree of freedom vehicle dynamics model is
Yr = Yr1 + Yr2 = krL arL + krR arR
also easy to describe the frequency characteristic. The
The lateral force obtained by product of the tire’s lat- main parameters such as gain and delay time can also
eral thrust coefficient and the camber angle is: be derived with the help of the three degree of freedom
vehicle dynamics model. Because the sine swept steer-

2Ycf = 2Kcf ff ing input simulation is used in this article, the vehicle’s
ð70Þ roll, lateral and yaw these three degrees of freedom are
2Ycr = 2Kcr fr
also mainly considered. The three degree of freedom
The total lateral force produced by the suspension vehicle dynamics model is:
movement and the bushing’s deformation is:
8      
>
> db _ d2k ∂af a _ ∂ar b _
>
> mV + c  ms hs 2 = 2Kf d + k  b  c + 2Kr k  b+ c
>
> dt dt  ∂k V ∂k V
>
>
>
> ∂k f ∂kr
>
>  2 Kcf f + Kcr k
>
> ∂k ∂k 
>
< dr 2
  
d k ∂af a ∂ar b
I  Ixz 2 = 2Kf d + k  b  c_ a  2Kr k b + c_ b ð72Þ
> dt
> dt  ∂k V  ∂k V
>
>
>
> ∂kf ∂kr
>
>  2 aKcf k  bKcr k
>
> ∂k ∂k
>
>  
>
> d2k dr db dk
>
: Ik 2  Ixz  ms hs V + c_ = ðKk + ms ghs Þk  Ck
dt dt dt dt
Gao and Wu 17

After Laplace transform of equation (72), the charac-

teristic equation is:
2 3
2(aKf  bKr )
mVs + 2(Kf + Kr ) mV +  ms hs s2  2Yk
6 V 7
6 7
6 2a2 Kf + b2 Kr 2 7 = 0 ð73Þ
4 2(aKf  bKr ) Is +  Ixz s  2Nk 5
V
2
 ms hs Vs  Ixz s  ms hs V Ik s + Ck s + (Kk  ms ghs )

Based on the formula (73), the objective function based

on the sine swept steering input simulation is:
 
 2(aKf  bKr ) 
 2Kf
 mV +  ms hs s2  2Yk 

 V 
2 2
 2a Kf + b Kr 
 2aKf Is +  Ixz s2  2Nk 
 V 
b(s)  0  Ixz s  ms hs V Ik s + Ck s + (Kk  ms ghs ) 
2
=  ð74Þ
d(s)  2(aKf  bKr ) 2 
 mVs + 2(Kf + Kr ) mV +  m h
s s s  2Y k 
 V 
 
 2a2 Kf + b2 Kr 2 
 2(aKf  bKr ) Is +  Ixz s  2Nk 
 V 
  ms hs Vs  Ixz s  ms hs V Ik s + Ck s + (Kk  ms ghs ) 
2

 
 mVs + 2(Kf + Kr ) 2Kf  ms hs s2  2Yk 
 
 2(aKf  bKr ) 2aKf 2
 Ixz s  2Nk 
 
_c(s)   ms hs Vs 0 Ik s + Ck s + (Kk  ms ghs ) 
2
=  ð75Þ
d(s)  2(aKf  bKr ) 2 
 mVs + 2(Kf + Kr ) mV +  ms hs s  2Yk 
 V 
 
 2a2 Kf + b2 Kr 2 
 2(aKf  bKr ) Is +  Ixz s  2Nk 
 V 
  ms hs Vs  Ixz s  ms hs V Ik s2 + Ck s + (Kk  ms ghs ) 
 
 2(l K  lr Kr ) 
 mVs + 2(Kf + Kr ) mV + f f 2K f 
 V
 
 2lf2 Kf + lr2 Kr 
 2(lf Kf  lr Kr ) Is + 2lf f 
K
 V
k(s)   ms hs Vs  Ixz s  ms hs V 0
=  ð76Þ
d(s)  2(lf Kf  lr Kr ) 
 mVs + 2(Kf + Kr ) mV +  ms hs s2  2Yk 
 V 
 2 2
2lf Kf + lr Kr 
 2 
 2(lf Kf  lr Kr ) Is +  I xz s  2N k 
 V 
 2 
 ms hs Vs  Ixz s  ms hs V Ik s + Ck s + (Kk  ms ghs )

In equations (74)–(76), Yk and Nk are:  b _


side slip angle ∂a ∂kr
∂k k  b + V c and camber angle ∂k k,
r

8     only consider the additional corners generated by the

>
> ∂af ∂ar ∂kf ∂kr
>
< Y k = K f + K r  K cf + K cr movement of the suspension, and do not take into
∂k ∂k ∂k ∂k account the impact of the deformation of the suspen-
   
>
> ∂a f ∂a r ∂kf ∂kr sion bushing on the wheel side slip angle and camber
> Nk =
: aKf  bKr  aKcf  bKcr
∂k ∂k ∂k ∂k angle. When considering the deformation of the sus-
ð77Þ pension bushing, the deformation of the bushing will
cause the compliance characteristic of the suspension to
In equation (72), the side slip stiffness Kf and Kr of the change. It can be seen from equation (69) that the sum
vehicle’s tires do not consider the influence of suspen- of the side slip stiffness of the left and right tires will
sion movement and bushing deformation. The front change accordingly, and the side slip stiffness of the
∂a
wheels’ side slip angle d + ∂kf k  b  Va c_ and cam- front axle is no longer Kf , the side slip stiffness of the
∂k
ber angle ∂kf k of the vehicle, as well as the rear wheels’ rear axle is no longer Kr . From equations (70) and (71),
18 Advances in Mechanical Engineering

Table 1. Vehicle parameters.

Parameter Value Parameter Value

Sprung mass (kg) 1067 Height from center of mass to ground (mm) 330
Unsprung mass (kg) 595 Distance between rear wheels (mm) 1487.8
Center of mass to front wheelbase (mm) 1452.6 Distance between front wheels (mm) 1508.3
Center of mass to rear wheelbase (mm) 1249.5 Wheelbase (mm) 2702.1

it can be seen that the side slip angle and camber angle There are eight bushings connecting the connection
will also change further, which will cause the transient rods and the vehicle body in Figures 10 and 13, four for
characteristics of the vehicle to change. pffiffiffiffiffiffiffi the front suspension and four for the rear suspension,
In equations (74)–(76), make s = jv, where j = 1, and each bushing has stiffness in six directions. Taking
the frequency response characteristic function of the the upper rear control arm on the front suspension as
vehicle can be obtained, including the amplitude- an example, the bushing is named K2 , then the displace-
frequency characteristic and
 the  phase-frequency
  char-
 ment stiffness of the bushing in the three directions x, y,
_ 
acteristic. Among them, cd ðvÞ, kd ðvÞ, and bd ðvÞ are z, and the rotational stiffness around the x, y, and z axis
the gains of the yaw rate, roll angle and the side slip are denoted as K2x , K2y , K2z ,K2tx , K2ty , and K2tz .
angle relative to the steering wheel angle respectively. There are eight bushings in Figures 10 and 12, so
The frequency corresponding to the maximum yaw rate there are 48 stiffness values, but some stiffness values
gain is the resonance frequency of the vehicle. The reso- have little effect on the handling stability of the vehicle,
nance frequency of the vehicle is set as one of the objec- and some stiffness values have a greater influence on
tive functions. uc ðvÞ, uk ðvÞ, and ub ðvÞ are the the handling stability of the vehicle. In order to find out
corresponding phase lag angles. Since the general driver the bushing stiffness values that have a greater impact
turns the steering wheel at a frequency of about 0.5 Hz, on the handling stability of the vehicle, the sensitivity
the gain of the yaw rate, roll angle, and side slip angle analysis of these 48 stiffness values is carried out first,
at 0.5 Hz relative to the steering wheel angle and the and then several stiffness values with higher sensitivity
corresponding delay time at 0.5 Hz are taken as opti- are selected for optimization analysis. After the sensitiv-
mize targets. The phase lag angle can be converted to ity analysis, five stiffness values with greater sensitivity
delay time by the following formula: are selected as shown in Table 2.

u + 180 u  180
t= or ð78Þ
360  f 360  f Optimization algorithm
The overall optimization objective functions are: By optimizing the stiffness of the suspension bushing to
improve the handling stability of the vehicle, the multi-
8  
> c_  objective optimization method and algorithm are
>
> f =  ð v Þ 
>
> 1  d  adopted. The NSGA-II algorithm is a relatively suc-
>
> v = 0:5Hz
>
>   cessful algorithm used in multi-objective optimization
>
> k 
> f2 = d ðvÞv = 0:5Hz
> algorithms.16 In non-dominated sorting, individuals
>
>
>
>   close to the Pareto frontier are selected, which enhances
>
> b 
>
> f =  ð v Þ  Pareto’s forward ability. Combined with the analysis of
< 3 d 
 _
v = 0:5Hz
  ð79Þ section 7.2, the mathematical expression of multi-
>
> c  objective optimization is as follows:
>
>
>

f4 = v  ðvÞ
>
> d
>
>
max
>
> u + 180 uc  180 F ð X Þ = fmin F1 ff1 ð X Þ, f2 ð X Þ, f3 ð X Þ, f5 ð X Þ, f6 ð X Þg
>
> f5 =
c
or
>
> 360  f 360  f + max F2 ff4 gg
>
>
>
>
>
> u + 180 uk  180 ð80Þ
: f6 = k or
360  f 360  f 
X = fx1 , x2 , . . . , x5 g
ð81Þ
0 ł xi ł 2:5, i = 1, 2, . . . , 5
Sensitivity analysis of the influence of suspension
F ð X Þ is the objective function, fk ð X Þ is the sub-objective
bushing’s stiffness on vehicle handling stability function, X is the optimized variable. In the NSGA-II
The vehicle’s parameters used in the optimization are algorithm, the calculation mechanism of crossover and
shown in Table 1: sudden mutation is called SBX method. The crossover
Gao and Wu 19

Table 2. Bushing stiffness with larger sensitivity value.

Stiffness Resonance Yaw rate Max yaw Roll angle Delay time of Delay time of
frequency (f4 ) gain (f1 ) rate gain (f3 ) gain (f2 ) lateral acceleration (f6 ) yaw rate (f5 )
z
KA1 0.97 1.69 20.19 1.62 2.08 2.16
tz
KA3 2.22 1.99 21.13 1.77 1.24 1.60
y
KA4 1.86 2.38 20.44 2.24 0.83 1.45
z
KA7 22.39 21.96 1.43 22.53 21.50 2.00
tz
KA8 1.63 1.69 21.18 1.96 2.05 2.28

operation in generating sub-individuals according to

the SBX method is:
8
>
> ð1, t + 1Þ 1 + bqi ð1, tÞ 1 + bqi ð2, tÞ
< xi = xi + xi
2 2 ð82Þ
>
> xð2, t + 1Þ = 1 + bqi xð1, tÞ + 1 + bqi xð2, tÞ
: i i i
2 2
The mutation operation in generating sub-individuals
according to the SBX method is:

ð1, t = 1Þ ð1, tÞ  
xi = xi + dq xUB LB
i  xi ð83Þ
8h i1
>
< 2u + ð1  2uÞð1  dÞhm + 1 1, u ł 0:5
hm + 1
dq = h i1 Figure 22. Iterative diagram of the maximum gain value of the
>
: 1  2ð1  uÞ + 2ðu  0:5Þð1  dÞhm + 1 , u.0:5
hm + 1
yaw rate relative to the steering wheel angle.

ð84Þ
  UB 
d = min xi  xLB UB
i , xi  xi xi  xLB
i , u 2 ½1, 0 ð85Þ

At present, the optimization objective function used in

the optimization of vehicle handling stability is mainly
realized by vehicle dynamic parameter modeling. In the
process of establishing parameter model, some second-
ary factors and degrees of freedom of the vehicle are
often ignored, which leads to the optimized model com-
pared with the actual vehicle, the error is larger. In
order to make the vehicle’s model as close as possible
to the real vehicle model, this paper adopts the multi-
body dynamics model, that is, the multi-objective opti-
mization of vehicle’s handling stability is carried out
for the multi-body dynamics model. Figure 23. Iterative graph of yaw rate gain at 0.5 Hz relative to
Since the sine swept steering input test is used in this the steering wheel angle.
paper to evaluate the handling stability of a vehicle, the
parameters of the sine swept steering input simulation
to evaluate the vehicle’s handling stability are taken as 0.5 Hz relative to the steering wheel angle. In the per-
the optimization objectives, and the suspension bushing formance of the frequency characteristics of the vehicle,
stiffness with higher sensitivity are taken as the optimi- the smaller the maximum gain value of the yaw rate
zation variables to perform multi objective optimiza- relative to the steering wheel angle, the more stable the
tion. After the multi objective optimization of the vehicle. Therefore, it is hoped that the maximum gain
whole vehicle, the optimization iterative diagram of the value of the yaw rate relative to the steering wheel angle
optimization targets is shown in Figures 22 to 27. will be minimized during optimization. The change
Figures 22 and 23 are iterative graphs of yaw rate trend of the iterative graphs in Figures 22 and 23 meets
gain with the steering wheel angle. Figure 22 is the max- the requirements of the optimization goal.
imum value of the yaw rate gain relative to the steering Figure 24 is an iterative diagram of the resonance
wheel angle, and Figure 23 is the yaw rate gain value at frequency. The larger the resonance frequency, the
20 Advances in Mechanical Engineering

Figure 24. Iterative diagram of resonance frequency. Figure 26. Iterative diagram of lateral acceleration delay time.

Figure 25. Iteration diagram of roll angle gain relative to the Figure 27. Iterative diagram of yaw rate delay time.
lateral acceleration.

more stable the vehicle. The resonance frequency in

Figure 24 is continuously increasing as the optimization
progresses, which meets the requirements of the optimi-
zation goal. Figure 25 shows the roll angle gain relative
to the lateral acceleration. In Figure 25, the roll angle
gain relative to the lateral acceleration is increased com-
pared to the original vehicle, but the increased value is
small, so it has little effect on the handling stability of
the vehicle.
Figures 26 and 27 show the delay time of lateral
acceleration relative to the steering wheel angle and the
delay time of yaw rate relative to lateral acceleration.
The smaller the delay time, the faster the response of
the vehicle. Figures 26 and 27 reflect that the delay time Figure 28. Pareto solution sets of yaw rate gain at 0.5 Hz with
is decreasing, which is beneficial to improve the resonance frequency.
response speed of the vehicle, indicating that the opti-
mization is carried out in accordance with the required
target. solution sets are obtained, as shown in Figures 28 to
33. Among them, the black point represents the general
solution, the blue point represents the Pareto solution,
Optimization results and the green point represents the relatively optimal
After the multi-objective optimization for the transient solution. Delay time 1 refers to the delay time of the
characteristic of the vehicle, the corresponding Pareto lateral acceleration relatives to the steering wheel angle,
Gao and Wu 21

Figure 32. Pareto solution sets of delay time 1 relative to

Figure 29. Pareto solution sets of roll angle gain at 0.5 Hz
delay time 2.
relative to yaw rate gain at 0.5 Hz.

Figure 30. Pareto solution sets of delay time 1 with yaw rate Figure 33. Pareto solution sets of delay time 2 relative to yaw
gain at 0.5 Hz. rate gain.

comparison of parameters before and after optimiza-

tion of the vehicle are shown in Table 3.

Comparative analysis of optimization

results of vehicle handling and stability
In order to analysis the optimization effect, a Pareto
solution is selected from the optimization results and
compared with original vehicle. In order to ensure that
the handling stability of the vehicle is optimized and
improved, comparisons are made from two aspects of
the vehicle’s transient characteristics and steady-state
characteristics.
Figure 31. Pareto solution sets of delay time 1 relative to
resonance frequency. Comparative analysis of vehicle’s transient
characteristics
and delay time 2 refers to the delay time of yaw rate Bring the selected Pareto solution into the original vehi-
relatives to the lateral acceleration. cle for simulation analysis, and compare the results
Through optimization analysis, Pareto solution sets with the original vehicle, the result comparison chart in
and relative optimal solution are finally obtained. The Figures 34 to 39 are obtained. Figure 34 shows the yaw
22 Advances in Mechanical Engineering

Table 3. Comparison of parameters before and after optimization.

f1 (s21) f2 (s21) f3 (s21) f4 (Hz) f5 (s) f6 (s) x1 x2 x3 x4 x5

Original vehicle 0.47 0.46 4.40 0.63 0.20 0.14 1 1 1 1 1

Optimized vehicle 0.44 0.42 4.52 0.83 0.17 0.14 1.31 0.001 0.03 0.18 0.21
Percentage change (%) 26.4 28.7 2.7 31.7 215 0 + 31 299.9 297 282 279

Figure 34. Yaw rate gain relative to the steering wheel angle.
Figure 36. Roll angle gain relative to the lateral acceleration.

difference to the delay time of yaw rate between opti-

mized and original vehicle is obvious, and it is reduced
after optimization compared with the original vehicle,
which meets the requirements of optimization goals.
Between 0 and 1 Hz, the optimized delay time of yaw
rate is basically smaller than the original vehicle. When
the frequency exceeds 1 Hz, the difference of yaw rate
delay time between the optimized vehicle and the origi-
nal vehicle is small. Since the frequency rarely reaches
above 1 Hz during the operation of the vehicle, the situ-
ation above 1 Hz will not be discussed too much.
Figure 36 shows the roll angle gain relative to the
lateral acceleration. In Figure 36, from low frequency
to high frequency, as the frequency increases, the differ-
Figure 35. Delay time of yaw rate relatives to the lateral ence of roll angle gain between the optimized vehicle
acceleration.
and original vehicle is relatively small. After enlarging
the curve at 0.5 Hz, the situation shown in Figure 37 is
obtained. In Figure 37, the roll angle gain relative to
rate gain relative to the steering wheel angle. In Figure
the lateral acceleration after optimization is increased
34, in the range of 0 to 1 Hz, especially at 0.5 Hz, the
than original vehicle, but the increased value is rela-
curve of the optimized yaw rate gain relative to the tively small and within an acceptable range. Because it
steering wheel angle is smaller than that of the original is a multi-objective optimization, it is normal for indi-
vehicle, which meets the requirements of the optimiza- vidual optimization objectives to fail to reach the
tion goal. When the frequency exceeds 1 Hz, the differ- optimum.
ence of the yaw rate gain before and after optimization Figure 38 shows the delay time of the lateral accel-
is not obvious. Since the frequency of the driver in the eration relative to the steering wheel angle, and after
process of driving the vehicle rarely exceeds the range enlarging the range around 0.5 Hz in Figure 38, the
of 1 Hz, the high frequency range is not considered. delay time of the lateral acceleration around 0.5 Hz
Figure 35 is the delay time of yaw rate relatives to relative to the steering wheel angle as shown in Figure
the lateral acceleration. In Figure 35, at 0.5 Hz, the 39 is obtained.
Gao and Wu 23

Figure 37. An enlarged view of roll angle gain relative to the Figure 39. The enlarged view of the lateral acceleration’s delay
lateral acceleration. time relative to the steering wheel angle.

Figure 40. The curve of steering wheel angle with the lateral
acceleration.
Figure 38. Delay time of lateral acceleration relative to the
steering wheel angle.
angle with the lateral acceleration shown in Figure 40
and the curve of the roll angle with the lateral accelera-
In Figures 38 and 39, the absolute value of the lat-
tion shown in Figure 41 are obtained.
eral acceleration’s delay time of the optimized vehicle is
In Figure 40, at the same lateral acceleration, the
reduced, but in general, the difference of the lateral
steering wheel angle of the optimized vehicle is larger,
acceleration’s delay time between the optimized vehicle
indicating that the understeer characteristic of the vehi-
and original vehicle in the entire frequency range is
cle is increased, and the vehicle tends to be more stable.
smaller. Because it is a multi-objective optimization, it
In Figure 41, under the same lateral acceleration, the
is difficult to achieve all the optimization objectives to
optimized vehicle’s roll angle has increased, but the
reach the optimal state, the lateral acceleration’s delay
increase value is not much and within in the acceptable
time after optimizing didn’t significantly reduced but it
range.
is acceptable.

Comparison and analysis of steady state with Conclusion

constant radius cornering simulation The stiffness calculation of vehicle suspension bushing
In order to further verify the improvement of the vehi- is analyzed and the multi-body dynamic model of vehi-
cle’s dynamics characteristics after the vehicle optimiza- cle’s front and rear suspension is established according
tion, the constant radius cornering simulation of the to the actual vehicle. The influence mechanism of the
whole vehicle before and after the optimization is com- bushing on the suspension wheel alignment angle and
pared and analyzed. The curve of the steering wheel the suspension force is analyzed in detail and
24 Advances in Mechanical Engineering

upper front control arm bushing in the z direc-

tion is reduced by 97%. For the rear suspension,
when the stiffness of the lower front control arm
bushing along the z direction is reduced by 82%
and the torsional stiffness of the lower rear con-
trol arm bushing around the z axis is reduced by
79%, the yaw rate gain at 0.5 Hz is reduced
by 6.4%, and the roll angle gain is increased by
2.7%, the delay time of yaw rate is reduced by
15%, the delay time of roll angle is unchanged,
and the resonance frequency is increased by
31.7%. From the perspective of steady-state
conditions, the understeer trend of the vehicle
increases, and the roll angle increases slightly,
Figure 41. The curve of roll angle with the lateral acceleration.
but it is within an acceptable range.

systematically. Through sensitivity analysis, five bushes Declaration of conflicting interests

that have a greater impact on the transient characteris- The author(s) declared no potential conflicts of interest with
tics of the vehicle are selected, and multi-objective opti- respect to the research, authorship, and/or publication of this
mization is performed for the selected bushings. The article.
specific conclusions obtained are:
Funding
1. Through the calculation of the vector, the link
The author(s) disclosed receipt of the following financial sup-
position and link angle during the movement of port for the research, authorship, and/or publication of this
the double wishbone suspension are analyzed in article: This project was supported by National Natural
detail, and the wheel alignment angle of the dou- Science Foundation of China (NSFC) (No.51965026),
ble wishbone suspension is finally calculated. Yunnan Province Applied Basic Research Foundation (No.
Through the simulation calculation of the multi- 2018FB097) and Scientific Research Fund Project of Yunnan
body dynamics model, the toe angle and camber Provincial Education Department (No. 2018JS022). The
angle change with the wheel travel are obtained. authors are greatly appreciated for the financial support.
2. The bushing deformation and bushing force of
double wishbone suspension and multi-link sus-
pension are analyzed, and the relationship ORCID iD
between bushing deformation and suspension Jin Gao https://orcid.org/0000-0001-6375-2865
compliance characteristics is explained.
3. Sensitivity analysis shows that for the double
wishbone suspension, the torsional stiffness of References
the lower front control arm bushing around the 1. Garcı́a Tárrago MJ, Kari L, Vinolas J, et al. Frequency
z axis, the translational stiffness of the lower and amplitude dependence of the axial and radial stiff-
rear control arm bushing along the y direction, ness of vehiclebon-black filled rubber bushings. Polymer
and the translational stiffness of upper front Test 2007; 26: 629–638.
control arm bushing along the z direction have a 2. Gibanica M, Abrahamsson TJS and Olsson M. Model
greater impact on the vehicle’s transient charac- updating of multiple nominally identical vehicle compo-
teristic. For the rear multi-link suspension, the nents. SAE Technical Paper 2018-01-0554, 2018.
3. Plagge J, Ricker A, Kröger NH, et al. Efficient modeling
stiffness of the lower front control arm bushing
of filled rubber assuming stress-induced microscopic
along the z direction and the torsional stiffness
restructurization. Int J Eng Sci 2020; 151: 1021–1028.
of the lower rear control arm bushing around 4. Sedlaczek K, Dronka S and Rauh J. Advanced modular
the z-axis have a greater impact on the vehicle’s modelling of rubber bushings for vehicle simulations.
transient characteristics. Veh Syst Dyn 2011; 49: 741–759.
4. For the front suspension, when the torsional 5. Aydemir E and Sendur P. Simplified transfer function
stiffness of the lower front control arm bushing approach for modeling frequency dependency of damp-
around the z axis increases by 31%, the transla- ing characteristics of rubber bushings. Proc Inst Mech
tional stiffness of the lower rear control arm Eng D J Automobile Eng 2019; 233: 2518–2531.
bushing along the y direction is reduced by 6. Zenglian Y, Lljun Z and Zhuoping Y. A novel empirical
99.9%, and the translational stiffness of the model of rubber bushing in automotive suspension
Gao and Wu 25

system. Int Conf Motion Vib Control 2010; 2010: _3A12 K Stiffness, N =m
to 1_. S Bushing displacement, m
7. Sancibrian R, Garcia P, Viadero F, et al. Kinematic r The radius of the bushing, m
design of double-wishbone suspension systems using a c Parameter
multiobjective optimisation approach. Veh Syst Dyn l The length of the rod, m
2010; 48: 793–813. p Parameter
8. Balike KP, Rakheja S and Stiharu I. Development of k Parameter
kineto-dynamic quarter-vehicle model for synthesis of a e unit vector
double wishbone suspension. Veh Syst Dyn 2011; 49:
R Rotation matrix
107–128.
9. Attia HA. Dynamic modelling of the double wishbone
F Force,N
motor-vehicle suspension system. Eur J Mech 2002; 21: Q Bushing displacement
167–174. M Torque, N  m
10. Tanik E and Parlaktasx V. On the analysis of double wish- m Mass,Kg
bone suspension. J Adv Mech Des 2016; 10: 1–11. a the distance between the front axle and
11. Shen Q, Lee W-T and Russell K. A small-scale optimiza- the vehicle’s center of mass, m
tion model for RSSR-SS motion generation with branch b the distance between the rear axle and the
and order defect elimination in Matlab. J Adv Mech Des vehicle’s center of mass, m
Syst Manuf 2014; 8: 1–8. I Moment of inertia, m4
12. Reddy KV, Kodati M, Chatra K, et al. A comprehensive V Velocity,m=s
kinematic analysis of the double wishbone and MacPher-
f objective function
son strut suspension systems. Mech Mach Theory 2016;
105: 441–470.
s Laplace operator
13. Tang L, Shangguan W-B and Dai L. A calculation
method of joint forces for a suspension considering non- Geek symbols
linear elasticity of bushings. Proc Inst Mech Eng K J
Multi-body Dyn 2012; 226: 281–297. x Angle of upper control arm, degree
14. Nabawy AE, Abdelrahman AA, Abdalla WS, et al. Anal- j Upper control arm rotation Angle, degree
ysis of the dynamic behavior of the double wishbone sus- F Upper control arm bushing rotation
pension system. Int J Appl Mech 2019; 11: 1–19. angle, degree
15. Wang X-L, Dai L and Shangguan W-B. Calculation of e Wheel hub rotation angle, degree
joint forces of a multi-link suspension for strength and l Rotation angle of vector, degree
fatigue analysis of bushings and control arms. Int J Veh m unit vector
Des 2014; 66: 217–234. a Toe angle, degree
16. Zhao Z, Liu B, Zhang C, et al. An improved adaptive u Rotation angle, degree
NSGA-II with multi-population algorithm. Appl Intell
u Camber angle, degree
2019; 49: 569–5801.
t Caster angle, degree
s Kingpin angle, degree
Appendix h Rotation angle, degree
r Rotation angle, degree
Notation
z Rotation angle, degree
n Bushing radial stiffness coefficient f Bushing rotation angle, degree
W load, N n Scale factor
G Rubber shear modulus, Pa c Yaw rate angle, degree
L Length,m