# The Study of Optimization and Matching to Spring and Antiroll Bar Stiffness of Suspension for Multiresponse Target of Whole Vehicle under Sine-Swept Steering Input

Authors: Jin Gao; Fuquan Wu
Journal: Mathematical Problems in Engineering
 (2020)

---

## Abstract

This article first leads from the specific double-wishbone suspension and multilink suspension structure form. And then systematically and detailedly analyse the change of spring's stiffness, and antiroll bar's stiffness causes the change of the side slip stiffness and rotation angle of tire, which will lead to the change of tire force, and then affect the dynamic characteristics of the whole vehicle. Based on this, the vehicle dynamics model considering the suspension is established, and the transfer function of the vehicle’s response index to steering wheel angle with coupling spring stiffness and antiroll bar stiffness is derived. Based on the dynamic theory analysis of the suspension and the whole vehicle, the multibody dynamics model of the whole vehicle with front double-wishbone suspension and rear multilink suspension was established. By calculating the frequency response characteristics of the vehicle under the sine-swept input, the frequency response index at the normal steering wheel operating frequency of 0.5 Hz was obtained. In addition, these frequency response indexes at 0.5 Hz were taken as optimization objectives, and the spring stiffness and antiroll bar stiffness of the front and rear suspension were taken as optimization variables, which were optimized by the NSGA-II algorithm. The results show that at 0.5 Hz, the gain value in the frequency response index is reduced, and the delay time is not significantly different from other group schemes, but it is not the worst; the value is within an acceptable range, and the dynamic characteristics of the car in the low frequency range have been improved.


---

## Body

## 1. Introduction

The characteristic of suspension is very important to the vehicle’s handling stability. The matching design of suspension’s stiffness is one of the important methods to improve vehicle’s handling stability. The matching of the suspension system’s roll stiffness has direct influence on the handling stability and safety of the vehicle. For cars with independent suspension, the roll stiffness of the suspension is mainly composed of the spring’s and the antiroll bar’s stiffness.Up to now, many experts and scholars have done a lot of research on the springs and antiroll bars of vehicle’s suspension, so as to improve the vehicle’s handling stability. Kazemi et al. studied the improvement of the operation handling stability of the whole vehicle whose front suspension is McPherson suspension and rear suspension is semitrailing arm suspension. A 9-DOF automotive dynamic model with McPherson suspension as the front suspension and semitrailing arm suspension as the rear suspension was established. After defining the offline optimization objective function and the online optimization objective function, the bee algorithm was used to optimize the vehicle’s handling stability by optimizing the geometric parameters of the suspension. Finally, the comparison results show that the offline objective function optimization is better than the online objective function optimization. The handing stability has been improved [1]. Mastinu et al. proposed a formula to describe the dynamic of suspension based on the 1/4 vehicle model and gave the standard deviations of the analytical forms of vehicle’s body acceleration, the relative displacement of sprung and unsprung mass, and the ground force, respectively. The invariant points of the frequency response functions of active suspension and passive suspension are derived. The Pareto formula for selecting appropriate suspension parameters is given. The analytical formula they presented is useful to understand the suspension system [2–5]. Šagi et al. developed a new multiobjective optimization model to determine the optimal parameters of the suspension system. The new optimization model developed is to integrate a rapid simulation tool for suspension kinematics analysis and vehicle dynamics analysis with appropriate accuracy into a multiobjective optimization environment. On the basis of identifying the parameters that affect the suspension and defining the criteria for evaluating the dynamic characteristics of the vehicle, the multiobjective optimization of the suspension is carried out. Finally, the FMOGA-II algorithm is used to obtain the best results in terms of convergence, number of solutions, calculation time, and Pareto Frontier [6]. Javanshir et al. built a model based on dynamics software Trucksim, and the geometric parameters of the suspension system were optimized with antiroll bars and coil springs to improve ride comfort and handling stability of vehicles. The optimized vehicle is able to pass through the target direction with the minimum possible deviation, the minimum lateral acceleration, and minimal lateral slip [7]. Fossati et al. proposed an effective method to optimize the design of the vehicle’s passive suspension system under random excitation and developed a numerical calculation program for vehicle state analysis. The objective function considers comfort and safety and combines with the NSGA-II algorithm for multiobjective optimization. The dynamic analysis results of the vehicle model are compared with the optimized and unoptimized suspension systems, and it is verified that the optimization can reduce the weighted RMS value of the vertical acceleration of the driver’s seat by up to 21.14%, while improving the safety of the car [8]. Gobbi et al. proposed an optimization algorithm based on the local approximation of objective function and constraint function, and the suspension system of ground vehicle is optimized to achieve the optimal balance through grip, comfort, working space, and turning performance. Numerical results show that this algorithm has better precision and higher efficiency than some widely used optimization methods [9]. Gadhvi et al. use the most extensive multiobjective optimization algorithms NSGA-II, SPEA2, and PESA-II to optimize the passive suspension. The comparison results show that the Pareto Frontier of NSGA-II algorithm achieves an extreme trade-off advantage, and the optimal value of the target vector is slightly better than SPEA2 and PESA-II. SPEA2 and PESA-II were superior in maintaining Pareto optimal solution diversity [10, 11]. Shi et al. by parameterizing the geometric dimensions of the antiroll bar and optimizing the design of the antiroll bar for a certain car model, it is concluded that the antiroll bar with large lateral dip stiffness is helpful to improve the roll stability of the car [12].The study of the above scholars improved the vehicle’s handling stability from the suspension system, the stiffness of suspension’s spring, and the antiroll bar. However, their researches are mainly based on the parameter-oriented model to improve the vehicle’s handling stability and rarely involve the specific suspension structure form. Firstly, there is no in-depth discussion on the suspension characteristics of specific suspension structure and no in-depth explanation on the influence mechanism of spring and antiroll bar stiffness on suspension characteristics and vehicle’s handling stability; it mainly includes three aspects: first, the above experts and scholars have not analyzed and clarified the relationship between the stiffness of spring and antiroll bar and the roll stiffness of suspension in detail; second, there is no detailed analysis of the roll center and instantaneous rotation center of the suspension, especially the analysis of instantaneous rotation center of the multilink suspension, which has a direct impact on the lateral load transfer; third, it did not analyze in detail the influence mechanism of suspension’s roll stiffness and pitch stiffness on load transfer and the influence mechanism of load transfer on tire side slip stiffness, nor did it discuss in detail the influence mechanism of suspension’s roll stiffness on wheel’s toe angle and deformation angle. In addition, the above scholars did not propose an effective and feasible optimization matching method for the stiffness of spring and antiroll bar to make the vehicle’s frequency characteristic indexes reach a better state for the specific suspension structure form.For problems that experts have not solved, this article will be divided into four parts to write a statement of the research work done: the first part, first of all, the relationship between the stiffness of spring and the antiroll bar and the roll stiffness of the suspension is analyzed in detail, and then analysis the roll center and instantaneous rotation center of the suspension is analyzed in detail, especially the instantaneous rotation center of the multilink suspension. The second part will analyze in detail the influence mechanism of suspension’s roll stiffness on tire’s side slip stiffness from the angle of lateral and longitudinal load transfer, and the influence mechanism of suspension roll stiffness and pitch stiffness on toe angle, deformation angle, camber angle, and vehicle’s lateral dynamics. In the third part, on the basis of the analysis in the previous two parts, the total lateral force generated by suspension on the tire was mainly described. After ignoring some factors, the multibody dynamic model of the whole vehicle with the front double-wishbone suspension and rear multilink suspension coupling spring stiffness and antiroll bar stiffness was established. In the fourth part, the spring stiffness and antiroll bar stiffness of the front and rear suspension are taken as optimization variables, and the frequency response indexes corresponding to 0.5 Hz of the whole vehicle under the of sine-swept input are taken as optimization target for multiobjective optimization, so as to carry out matching design for the suspension spring stiffness and antiroll bar stiffness. Finally, Pareto solution set and relative optimal solution of frequency response indexes under sine-swept input are obtained. The relative optimal solution is compared with the original vehicle’s simulation data to verify whether the frequency response performance of the car has been improved.


## 2. Analysis of the Force of Suspension’s Spring and the Antiroll Bar

The stiffness of the suspension’s spring and the antiroll bar has important effect on the roll characteristic of the vehicle. And it greatly affects and determines the roll stiffness of the suspension. The installation of coil spring and the force of suspension are shown in Figure1 [13].Figure 1
The installation of coil spring and the force of suspension.In Figure1, F is the support force of the suspension from the ground; FL is the force on the suspension’s transverse arm from the vehicle’s body; F′ is the force on the shock absorber from the vehicle’s body; and F′ can be decomposed into the force FF that acts on the axis of shock absorber and the force FQ that acts on the axis that perpendicular to the axis of shock absorber. The force FF acting on the shock absorber’s axle is supported by the shock absorber’s spring, and the force FQ acting on the axle that perpendicular to the shock absorber’s axle is borne by the shock absorber’s upper end. η is the angle between the force F′ that acts on the shock absorber and the shock absorber’s axle.The force of the antiroll bar is shown in Figure2. When the vertical displacement of the right and left wheel is different, the antiroll bar will produce a torsional stiffness Kφw. And then the rolling force FA that acts on the end of the antiroll bar will be produced.Figure 2
The force diagram of the car’s antiroll bar.According to the force diagram of the antiroll bar, in the moment of the car’s rolling, the torque acting on the sprung mass and unsprung mass is(1)MAR=Kφwlcl0le2⋅ϕ−Kφwlc2le2⋅ϕu,where lc is the length of the antiroll bar’s straight rod part; l0 is the distance between the two rubber bushes; and le is the vertical displacement between the antiroll bar’s endpoint and the straight rod part. ϕ is the roll angle of sprung mass; ϕu is the roll angle of unspung mass.


## 3. Analysis of Suspension’s Roll Stiffness and Roll Center

The roll stiffness and roll center are important parameters of suspension which can affect the loading transfer between left and right wheels and make tire’s side slip stiffness change in the course of the vehicle’s motion. For the vehicle in this paper, the front suspension is a double-wishbone suspension and the rear suspension is multilink suspension.


### 3.1. Roll Stiffness Analysis of the Suspension

The total roll stiffness of the car is mainly composed of the roll stiffness of the front suspension and the rear suspension, that is,(2)Kϕ=Kϕf+Kϕr,where Kϕ is the total roll stiffness of the car; Kϕf is the roll stiffness of the front suspension; and Kϕr is the roll stiffness of the rear suspension. The roll stiffness of front suspension and rear suspension are combined with the stiffness that the spring contributed and the stiffness that the antiroll bar contributed, that is,(3)Kϕf=Kϕfs+Kϕfw,(4)Kϕr=Kϕrs+Kϕrw,where Kϕfs and Kϕrs are the roll stiffness that the front and rear suspension’s spring contributed, respectively, and Kφfw and Kϕrw are the roll stiffness that the front and rear suspension’s antiroll bar contributed, respectively.Figure3 shows the mechanical model where the spring stiffness of the double-wishbone independent suspension contributes to the roll stiffness of the suspension. The sprung mass is fixed, and the ground rotates an angle Φ around the intersection E of the car’s center line and the ground.Figure 3
The dynamic model of double-wishbone suspension’s roll stiffness.It can be obtained from Figure3 that the roll stiffness is(5)Kϕf=dTdΦ=12⋅Ks⋅c⋅n1⋅tfp⋅n2⋅cosγ2,where T is the torque that the ground acts on the suspension to prevent it from rotating. ΔWf is the increment of the force of the wheel to the ground. tf is the front wheel center distance. Ks is the spring’s linear stiffness. γ is the intersection angle. The point M′ is the instantaneous rotation center. p,, c, n1, and n2 are the lever’s length.In order to analysis the roll stiffness and the change of the roll center of double-wishbone suspension and multilink suspension during the driving of the car. Dynamic models of the car’s front double-wishbone suspension and rear multilink suspension shown in Figures4 and 5 are built in ADAMS/CAR software [14].Figure 4
The dynamic model of vehicle’s double-wishbone suspension.Figure 5
The dynamic model of vehicle’s multilink suspension.Figures6 and 7 show the influence trend of antiroll bar’s stiffness and spring’s stiffness on the suspension’s roll stiffness. In Figure 6, the front and rear suspensions increase in a nonlinear manner as the stiffness of the antiroll bar increases. When the antiroll bar stiffness gradually increases in the range of 40 N/m and 60 N/m, the roll stiffness of the front suspension and rear suspension increases slowly. When the antiroll bar stiffness increases to 60 N/m, the roll stiffness of the front suspension and rear suspension will increase at a large rate as the stiffness of the antiroll bar increases. In Figure 7, the roll stiffness of front suspension and rear suspension increases linearly with the increase of spring stiffness of front and rear suspensions.Figure 6
The influence trend of antiroll bar’s stiffness on suspension’s roll stiffness.Figure 7
The influence trend of spring’s stiffness on suspension’s roll stiffness.


### 3.2. Analysis of the Multilink Suspension’s Roll Center

Figure8 shows the solid model of multilink suspension. The motion of the multilink suspension is complicated. In order to find the multilink suspension’s instantaneous rotation axis, some necessary simplifications are needed [15].Figure 8
The solid model of multilink suspension.Figure9 shows a schematic diagram after equivalent transformation of the actual multilink suspension. Plane a is determined through three points ①, ③, and ② of the lower suspension arm. The extension line of extroverted arm ④ and ⑤ intersects plane a at point A. Making the plane d perpendicular to extroverted arm through point 5, the axle 1 can be obtained through points ① and ②, and the axle 2 can be obtained through points ③ and ⑧. V3 and V8 can be obtained as point ③, and point ⑧ rotates around axle 1. It is easy to get V53 as point 5 exercises relative to point ③. Make the plane e that parallel to V3 through point 5 and point V53. The alternating line of plane d and plane e is the motion direction of point 5 shown as V5. Make plane f perpendicular to V5 through point 5. The alternating line of plane f and plane a is the instantaneous roll axle during multilink suspension’s motion.Figure 9
Schematic diagram of the instantaneous roll axle of the multilink suspension.The instantaneous axis of rotation of the suspension intersects the cross section through the axle of the car at a point, which is the instantaneous center of rotation. The point where the connection line between the instantaneous rotation center and the bottom of the suspension’s tire intersects the longitudinal section of the car is the roll center. As shown in Figure10, it is the change of the roll center’s height of the suspension with the tire jumping. In Figure 10, the roll center’s height of the rear suspension is higher than that of the front suspension, and the difference of the roll center’s height of the front and rear suspensions becomes larger and larger as the wheel is jumping.Figure 10
The change of the roll center’s height of the suspension with the tire jumping.


## 3.1. Roll Stiffness Analysis of the Suspension

The total roll stiffness of the car is mainly composed of the roll stiffness of the front suspension and the rear suspension, that is,(2)Kϕ=Kϕf+Kϕr,where Kϕ is the total roll stiffness of the car; Kϕf is the roll stiffness of the front suspension; and Kϕr is the roll stiffness of the rear suspension. The roll stiffness of front suspension and rear suspension are combined with the stiffness that the spring contributed and the stiffness that the antiroll bar contributed, that is,(3)Kϕf=Kϕfs+Kϕfw,(4)Kϕr=Kϕrs+Kϕrw,where Kϕfs and Kϕrs are the roll stiffness that the front and rear suspension’s spring contributed, respectively, and Kφfw and Kϕrw are the roll stiffness that the front and rear suspension’s antiroll bar contributed, respectively.Figure3 shows the mechanical model where the spring stiffness of the double-wishbone independent suspension contributes to the roll stiffness of the suspension. The sprung mass is fixed, and the ground rotates an angle Φ around the intersection E of the car’s center line and the ground.Figure 3
The dynamic model of double-wishbone suspension’s roll stiffness.It can be obtained from Figure3 that the roll stiffness is(5)Kϕf=dTdΦ=12⋅Ks⋅c⋅n1⋅tfp⋅n2⋅cosγ2,where T is the torque that the ground acts on the suspension to prevent it from rotating. ΔWf is the increment of the force of the wheel to the ground. tf is the front wheel center distance. Ks is the spring’s linear stiffness. γ is the intersection angle. The point M′ is the instantaneous rotation center. p,, c, n1, and n2 are the lever’s length.In order to analysis the roll stiffness and the change of the roll center of double-wishbone suspension and multilink suspension during the driving of the car. Dynamic models of the car’s front double-wishbone suspension and rear multilink suspension shown in Figures4 and 5 are built in ADAMS/CAR software [14].Figure 4
The dynamic model of vehicle’s double-wishbone suspension.Figure 5
The dynamic model of vehicle’s multilink suspension.Figures6 and 7 show the influence trend of antiroll bar’s stiffness and spring’s stiffness on the suspension’s roll stiffness. In Figure 6, the front and rear suspensions increase in a nonlinear manner as the stiffness of the antiroll bar increases. When the antiroll bar stiffness gradually increases in the range of 40 N/m and 60 N/m, the roll stiffness of the front suspension and rear suspension increases slowly. When the antiroll bar stiffness increases to 60 N/m, the roll stiffness of the front suspension and rear suspension will increase at a large rate as the stiffness of the antiroll bar increases. In Figure 7, the roll stiffness of front suspension and rear suspension increases linearly with the increase of spring stiffness of front and rear suspensions.Figure 6
The influence trend of antiroll bar’s stiffness on suspension’s roll stiffness.Figure 7
The influence trend of spring’s stiffness on suspension’s roll stiffness.


## 3.2. Analysis of the Multilink Suspension’s Roll Center

Figure8 shows the solid model of multilink suspension. The motion of the multilink suspension is complicated. In order to find the multilink suspension’s instantaneous rotation axis, some necessary simplifications are needed [15].Figure 8
The solid model of multilink suspension.Figure9 shows a schematic diagram after equivalent transformation of the actual multilink suspension. Plane a is determined through three points ①, ③, and ② of the lower suspension arm. The extension line of extroverted arm ④ and ⑤ intersects plane a at point A. Making the plane d perpendicular to extroverted arm through point 5, the axle 1 can be obtained through points ① and ②, and the axle 2 can be obtained through points ③ and ⑧. V3 and V8 can be obtained as point ③, and point ⑧ rotates around axle 1. It is easy to get V53 as point 5 exercises relative to point ③. Make the plane e that parallel to V3 through point 5 and point V53. The alternating line of plane d and plane e is the motion direction of point 5 shown as V5. Make plane f perpendicular to V5 through point 5. The alternating line of plane f and plane a is the instantaneous roll axle during multilink suspension’s motion.Figure 9
Schematic diagram of the instantaneous roll axle of the multilink suspension.The instantaneous axis of rotation of the suspension intersects the cross section through the axle of the car at a point, which is the instantaneous center of rotation. The point where the connection line between the instantaneous rotation center and the bottom of the suspension’s tire intersects the longitudinal section of the car is the roll center. As shown in Figure10, it is the change of the roll center’s height of the suspension with the tire jumping. In Figure 10, the roll center’s height of the rear suspension is higher than that of the front suspension, and the difference of the roll center’s height of the front and rear suspensions becomes larger and larger as the wheel is jumping.Figure 10
The change of the roll center’s height of the suspension with the tire jumping.


## 4. Analysis of Lateral and Longitudinal Load Transfer during Vehicle’s Motion

Car’s body will roll when the car turns and changes line, and the roll angle should be controlled in a reasonable range in car’s movement. Therefore, it is necessary to design reasonable suspension’s roll stiffness. The rolling of the car’s body is shown in Figure11. Figure 12 shows the rolling model of the front axle of the car, and this model is obtained by the simplification of the front suspension of this paper.Figure 11
Vehicle’s rolling model.Figure 12
Vehicle’s front axle’s rolling model.In Figure11, ms is the sprung mass. mu is the unsprung mass. Aj is the lateral acceleration of sprung mass. ay is the lateral acceleration of unsprung mass. hs is the distance between the mass center of sprung mass and the roll center os. h0∗ is the distance between the suspension’s roll center and the road surface. r is the wheel’s radius. ϕ is the roll angle of the vehicle’s body. The torque caused by the suspension centrifugal force is M1=msAjhs. The roll torque caused by the mass gravity of the sprung mass is M2=msghsϕ. For independent suspension, the roll torque caused by the centrifugal force of unspung mass is M3=−muayh0∗−r.In Figure11, according to the torque equilibrium condition, the following can be obtained:(6)Kϕϕ=msAjhs+msghsϕ−muayh0∗−r.According to equation (6), the total roll stiffness and steady-state roll angle of the suspension can be obtained as follows:(7)Kϕ=msAjhs+muayh0∗−rϕ+msgh,(8)ϕ=msAjhs−muayh0∗−rKϕf+Kϕr−msghs.For the simplified suspension model shown in Figure12, the transformation between the linear stiffness and angular stiffness of the front and rear suspension’s spring is(9)Kϕfs=12Kfstfl2fl1f2,(10)Kϕrs=12Krstrl2rl1r2,where Kfs and Krs are the linear stiffness of the front and rear suspension’s springs, respectively. l1f is the length of the fore arm. l1r is the length of the rear arm. l2f is the distance between the spring’s installation position of the front suspension and the hinge point of the arm. l2r is the distance between the spring’s installation position of the rear suspension and the hinge point of the arm. tf and tr are the wheel center distance of the front and rear suspension, respectively.Substituting equations (3) and (4) into equation (2), the total roll stiffness Kϕ of vehicle can be obtained:(11)Kϕ=Kϕf+Kϕr=Kϕfs+Kϕfw+Kϕrs+Kϕrw.For the inertial forces acting on the vehicle, there must be an equilibrium tire force. The inertial force acting on the center of mass decomposes to the front and rear wheels:(12)msAj+muayaa+b,(13)msAj+muayba+b,where a is the distance from the vehicle’s center of mass to the front axle. b is the distance from the vehicle’s center of mass to the rear axle.If the vehicle’s body rolls, load transfer will occur between the left and right wheels of the front and rear axles. Let the load transfer quantities of the front axle and rear axle of the vehicle beΔWf and ΔWr, respectively. In the plane of the front and rear wheels perpendicular to the longitudinal direction of the vehicle, it can be obtained from the torque balance around the corresponding roll center:(14)Kϕfϕ=ΔWftf−msAj+muayba+bhf,(15)Kϕrϕ=ΔWrtf−msAj+muayaa+bhr.Substitute equation (8) into equations (14) and (15), ΔWf and ΔWr can be obtained as follows:(16)ΔWf=KϕftfmsAjhs−muayhof∗−rKϕf+Kϕr−msghs+msAj+muayba+bhf,(17)ΔWr=KϕrtrmsAjhs−muayhor∗−rKϕf+Kϕr−msghs+msAj+muayaa+bhr.The load transfer of the car not only occurs between the left and right wheels but also happens between the front and rear axles of the car. When the car is accelerating or braking, load transfer will occur between the front and rear axles. The load transfer between the front and rear axles of the vehicle is(18)ΔW=hms+mux¨a+b,where h is the distance between the vehicle’s center of mass and the ground, h=hs+h0∗.


## 5. Effect of Load Transfer and Longitudinal Force on Tire’s Side Slip Stiffness

If the vehicle’s body rolls, load transfer will occur between the left and right wheels. When the car’s body pitches, the load transfer will occur between the front and rear axles. The tire’s side slip stiffness will be affected when the load transfer occurs, and vehicle’s handling stability will be affected too.


### 5.1. Change of Tire’s Side Slip Stiffness While Only considering Load Transfer between Left and Right Wheels

According to the analysis in the above part, it is known that load transfer will occur between the left and right wheels when the vehicle is rolling, and the load transfer will cause the change of the tires’ side slip stiffness. Taking the front axle as an example, when no lateral force is applied to the vehicle, the vertical load of the left and right wheels of the axle isWf0, and the side slip stiffness of each tire is Kf0. When there is lateral force acting on the car and the ground, there is a corresponding reaction force Ffy acting on the front axle’s two tires. If the vertical load of the left and right wheels does not change, the corresponding side slip angle is(19)λf0=Ffy2Kf0.In fact, the vertical load of the left and right wheels will change when the lateral force acts on the vehicle. The inside wheel reducesΔWf, and the outboard wheel increases ΔWf; the side slip stiffness of the two tires then becomes Kf1 and Kf2. Due to left and right tire’s side slip angle is equal,(20)Ffy=Kf1λf+Kf2λf.MakeKf0′=Kf1+Kf2/2, where Kf0′ is the average side slip stiffness of each tire after vertical load redistribution, so the side slip angle of the left and right wheels is(21)λf=Ffy2Kf0′.Due to the change curve of the side slip stiffness of the left and right tires is a smooth curve protruding upward, from the geometrical point of view,Kf0>Kf0′ and λf>λf0. Further analysis shows that the greater the vertical load difference between the left and right wheels, the smaller the average side slip stiffness. The side slip stiffness of the rear tires is the same as the front tires.Thus, under the action of lateral force, if the vertical load of the front axle’s left and right wheels of the car changes greatly, cars tend to increase understeer. If the vertical load of the left and right wheels of the rear axle changes greatly, the car tends to reduce the understeer. The load variation of the left and right wheels of the front and rear axles of a car is determined by a series of parameters such as the roll stiffness of the front and rear suspensions, the sprung mass, the unspung mass, the center of mass’ position, and the roll center of the front and rear suspensions.


### 5.2. Tires’ Side Slip Stiffness While considering Longitudinal Load Transfer and Longitudinal Forces

The previous analysis is based on the conclusion that longitudinal load transfer and longitudinal force are ignored. However, there is a longitudinal load transfer during the actual motion of the car. If the influence of load transfer on roll stiffness is small, it can be regarded as the first order. If the driving/braking force is a small relative to the tire load, the side slip stiffness of the tire can be simulated by a simple parabolic function. For small longitudinal and lateral accelerations,(22)∂Kf∂WΔWfKf0,∂Kf∂WΔWKf0,2XfμWf2.They are treated as tiny quantities of the same order of magnitude. Based on these simplifications, the side slip stiffness can be expressed by the following formula at a small side slip angle:(23)Kf1=Kf01−∂Kf∂WΔWfKf0−∂Kf∂WΔW2Kf0−122XfμWf2,Kf2=Kf01+∂Kf∂WΔWfKf0−∂Kf∂WΔW2Kf0−122XfμWf2,Kr1=Kr01−∂Kr∂WΔWrKr0−∂Kr∂WΔW2Kr0−122XrμWr2,Kr2=Kr01+∂Kr∂WΔWrKr0−∂Kr∂WΔW2Kr0−122XrμWr2,where W is the mass of the car, W=ms+mu. ΔWf is the load transfer amount between the left and right wheels of the front axle when the car turns. ΔW is the amount of load transfer between the front and rear axles caused by the driving force or braking force, and Wf is the vertical load of the front axle Wf=bW/a+b. μ is the coefficient of friction between the tire and the road. Xf and Xr are the longitudinal force on the front and rear tires of the vehicle; Xf=αcWx¨/2, and Xr=1−αcWx¨/2. αc is the distribution coefficient of driving force or braking force. l is the wheelbase of the car, l=a+b. x¨ is the longitudinal acceleration caused by driving force or braking force. h is the height between the vehicle’s center of mass and the ground. ΔWr is the load transfer amount of the left and right wheels of the rear axle. Wr is the vertical load of the rear axle, Wr=aW/a+b.From the above analysis, it can be known that the equivalent lateral stiffness of the front and rear axles can be obtained by superposing the lateral stiffness of the left and right tires, respectively. Ifx¨ is used to represent the equivalent side slip stiffness of the front and rear axles of the car, then(24)2Kf∗≈2Kf01+hW2lKf0∂Kf∂Wx¨+12αclub2x¨2,2Kr∗≈2Kr01−hW2lKr0∂Kf∂Wx¨+121−αclua2x¨2.


## 5.1. Change of Tire’s Side Slip Stiffness While Only considering Load Transfer between Left and Right Wheels

According to the analysis in the above part, it is known that load transfer will occur between the left and right wheels when the vehicle is rolling, and the load transfer will cause the change of the tires’ side slip stiffness. Taking the front axle as an example, when no lateral force is applied to the vehicle, the vertical load of the left and right wheels of the axle isWf0, and the side slip stiffness of each tire is Kf0. When there is lateral force acting on the car and the ground, there is a corresponding reaction force Ffy acting on the front axle’s two tires. If the vertical load of the left and right wheels does not change, the corresponding side slip angle is(19)λf0=Ffy2Kf0.In fact, the vertical load of the left and right wheels will change when the lateral force acts on the vehicle. The inside wheel reducesΔWf, and the outboard wheel increases ΔWf; the side slip stiffness of the two tires then becomes Kf1 and Kf2. Due to left and right tire’s side slip angle is equal,(20)Ffy=Kf1λf+Kf2λf.MakeKf0′=Kf1+Kf2/2, where Kf0′ is the average side slip stiffness of each tire after vertical load redistribution, so the side slip angle of the left and right wheels is(21)λf=Ffy2Kf0′.Due to the change curve of the side slip stiffness of the left and right tires is a smooth curve protruding upward, from the geometrical point of view,Kf0>Kf0′ and λf>λf0. Further analysis shows that the greater the vertical load difference between the left and right wheels, the smaller the average side slip stiffness. The side slip stiffness of the rear tires is the same as the front tires.Thus, under the action of lateral force, if the vertical load of the front axle’s left and right wheels of the car changes greatly, cars tend to increase understeer. If the vertical load of the left and right wheels of the rear axle changes greatly, the car tends to reduce the understeer. The load variation of the left and right wheels of the front and rear axles of a car is determined by a series of parameters such as the roll stiffness of the front and rear suspensions, the sprung mass, the unspung mass, the center of mass’ position, and the roll center of the front and rear suspensions.


## 5.2. Tires’ Side Slip Stiffness While considering Longitudinal Load Transfer and Longitudinal Forces

The previous analysis is based on the conclusion that longitudinal load transfer and longitudinal force are ignored. However, there is a longitudinal load transfer during the actual motion of the car. If the influence of load transfer on roll stiffness is small, it can be regarded as the first order. If the driving/braking force is a small relative to the tire load, the side slip stiffness of the tire can be simulated by a simple parabolic function. For small longitudinal and lateral accelerations,(22)∂Kf∂WΔWfKf0,∂Kf∂WΔWKf0,2XfμWf2.They are treated as tiny quantities of the same order of magnitude. Based on these simplifications, the side slip stiffness can be expressed by the following formula at a small side slip angle:(23)Kf1=Kf01−∂Kf∂WΔWfKf0−∂Kf∂WΔW2Kf0−122XfμWf2,Kf2=Kf01+∂Kf∂WΔWfKf0−∂Kf∂WΔW2Kf0−122XfμWf2,Kr1=Kr01−∂Kr∂WΔWrKr0−∂Kr∂WΔW2Kr0−122XrμWr2,Kr2=Kr01+∂Kr∂WΔWrKr0−∂Kr∂WΔW2Kr0−122XrμWr2,where W is the mass of the car, W=ms+mu. ΔWf is the load transfer amount between the left and right wheels of the front axle when the car turns. ΔW is the amount of load transfer between the front and rear axles caused by the driving force or braking force, and Wf is the vertical load of the front axle Wf=bW/a+b. μ is the coefficient of friction between the tire and the road. Xf and Xr are the longitudinal force on the front and rear tires of the vehicle; Xf=αcWx¨/2, and Xr=1−αcWx¨/2. αc is the distribution coefficient of driving force or braking force. l is the wheelbase of the car, l=a+b. x¨ is the longitudinal acceleration caused by driving force or braking force. h is the height between the vehicle’s center of mass and the ground. ΔWr is the load transfer amount of the left and right wheels of the rear axle. Wr is the vertical load of the rear axle, Wr=aW/a+b.From the above analysis, it can be known that the equivalent lateral stiffness of the front and rear axles can be obtained by superposing the lateral stiffness of the left and right tires, respectively. Ifx¨ is used to represent the equivalent side slip stiffness of the front and rear axles of the car, then(24)2Kf∗≈2Kf01+hW2lKf0∂Kf∂Wx¨+12αclub2x¨2,2Kr∗≈2Kr01−hW2lKr0∂Kf∂Wx¨+121−αclua2x¨2.


## 6. Change of Wheels’ Toe Angle and Deformation Steering Angle

If the vehicle has longitudinal accelerationx¨ and lateral acceleration y¨ at meantime, the body will have pitch and roll motion. This will cause the suspension’s stroke to change, resulting in the change of the toe angle and the deformation steering angle.


### 6.1. Change of Wheels’ Toe Angle

Suppose the vehicle’s pitch and roll axes are on the ground, and the toe angle’s change of the front wheel caused by the pitching motion is∂αf/∂zahWx¨/Kθ. Similarly, the toe angle’s change caused by the roll motion is ∂αf/∂ztfhWy¨/2Kϕ. ∂αf/∂z is the toe angle’s change in the unit suspension’s stroke. Kθ is the pitch stiffness. Kϕ is the roll stiffness.Figure13 shows the curve of the front and rear suspension’s toe angle with the wheel’s displacement when the wheels’ bouncing in the reverse direction [16].Figure 13
The curve of the toe angle of the front and rear wheels with the wheel’s displacement.In Figure13, during the downward bouncing of the wheels, the toe angle of the front suspension wheels (the inner wheel when turning) increases in the positive direction, and the toe angle of the rear suspension wheels (the inner wheel when turning) increases in the negative direction. The toe angle of the corresponding outer wheels of the front suspension increases negatively, and the toe angle of the outer wheels of the rear suspension increases positively. For the front suspension, the toe angle of the inner wheels increases positively, and the toe angle of the outer wheels increases negatively, which will increase the understeer characteristics of the car. For the rear suspension, the toe angle of the inner wheel increases negatively, and the toe angle of the outer wheel increases positively, which will also increase the understeer characteristics of the car.


### 6.2. Change of Wheels’ Deformation Steering Angle

In addition to of the toe angle’ change, the torqueTs applied to the steering system due to the tire’s lateral force 2Fy=b/lWy¨ and longitudinal force Xf=αcWy¨ acting on the front wheels will cause deformation steering angle to change, where Ts is(25)Ts=2ξFy+FyKy2Xf=ξ+αcW2Kyx¨bWly¨,where Ky is the side slip stiffness of the tire. ξ is the sum of the tire drag and the rearward drag of the kingpin. For each wheel, the sum of the toe angle’s change, and the deformation steering angle can be expressed as follows:(26)αf1=∂αf∂zahWKθx¨+∂αf∂ztfhW2Kϕy¨−∂αf∂XαcW2x¨−∂αf∂Tξ+αcW2Kyx¨bWly¨,(27)αf2=−∂αf∂zahWKθx¨+∂αf∂ztfhW2Kϕy¨−∂αf∂XαcW2x¨−∂αf∂Tξ+αcW2Kyx¨bWly¨,(28)αr1=−∂αr∂zahWKθx¨+∂αr∂ztrhW2Kϕy¨−∂αr∂X1−αcW2x¨,(29)αr2=∂αr∂zbhWKθx¨+∂αr∂ztrhW2Kϕy¨−∂αr∂X1−αcW2x¨,where ∂αf/∂X and ∂αr/∂X are the deformation steering angles of the unit longitudinal force of the front and rear suspension systems, respectively. ∂αf/∂X is the deformation steering angle of the unit steering torque.Figure14 shows the change of vehicle’s toe angle with longitudinal force. In Figure 14, when the longitudinal force is positive, it is the driving force, and when the longitudinal force is negative, it is the braking force. As the braking force increases, the toe angle of the front suspension increases in the negative direction, and the toe angle of the rear suspension increases in the positive direction. As the driving force increases, the toe angle of the front suspension increases positively, and the toe angle of the rear suspension increases negatively.Figure 14
The curve of vehicle’s toe angle with longitudinal force.Figure15 shows the curve of the front and rear suspensions’ toe angle with the aligning torque. In Figure 15, as the wheel’s aligning force increases, the toe angle of the front suspension and the rear suspension decreases correspondingly.Figure 15
The change of the car’s toe angle with the aligning force.


## 6.1. Change of Wheels’ Toe Angle

Suppose the vehicle’s pitch and roll axes are on the ground, and the toe angle’s change of the front wheel caused by the pitching motion is∂αf/∂zahWx¨/Kθ. Similarly, the toe angle’s change caused by the roll motion is ∂αf/∂ztfhWy¨/2Kϕ. ∂αf/∂z is the toe angle’s change in the unit suspension’s stroke. Kθ is the pitch stiffness. Kϕ is the roll stiffness.Figure13 shows the curve of the front and rear suspension’s toe angle with the wheel’s displacement when the wheels’ bouncing in the reverse direction [16].Figure 13
The curve of the toe angle of the front and rear wheels with the wheel’s displacement.In Figure13, during the downward bouncing of the wheels, the toe angle of the front suspension wheels (the inner wheel when turning) increases in the positive direction, and the toe angle of the rear suspension wheels (the inner wheel when turning) increases in the negative direction. The toe angle of the corresponding outer wheels of the front suspension increases negatively, and the toe angle of the outer wheels of the rear suspension increases positively. For the front suspension, the toe angle of the inner wheels increases positively, and the toe angle of the outer wheels increases negatively, which will increase the understeer characteristics of the car. For the rear suspension, the toe angle of the inner wheel increases negatively, and the toe angle of the outer wheel increases positively, which will also increase the understeer characteristics of the car.


## 6.2. Change of Wheels’ Deformation Steering Angle

In addition to of the toe angle’ change, the torqueTs applied to the steering system due to the tire’s lateral force 2Fy=b/lWy¨ and longitudinal force Xf=αcWy¨ acting on the front wheels will cause deformation steering angle to change, where Ts is(25)Ts=2ξFy+FyKy2Xf=ξ+αcW2Kyx¨bWly¨,where Ky is the side slip stiffness of the tire. ξ is the sum of the tire drag and the rearward drag of the kingpin. For each wheel, the sum of the toe angle’s change, and the deformation steering angle can be expressed as follows:(26)αf1=∂αf∂zahWKθx¨+∂αf∂ztfhW2Kϕy¨−∂αf∂XαcW2x¨−∂αf∂Tξ+αcW2Kyx¨bWly¨,(27)αf2=−∂αf∂zahWKθx¨+∂αf∂ztfhW2Kϕy¨−∂αf∂XαcW2x¨−∂αf∂Tξ+αcW2Kyx¨bWly¨,(28)αr1=−∂αr∂zahWKθx¨+∂αr∂ztrhW2Kϕy¨−∂αr∂X1−αcW2x¨,(29)αr2=∂αr∂zbhWKθx¨+∂αr∂ztrhW2Kϕy¨−∂αr∂X1−αcW2x¨,where ∂αf/∂X and ∂αr/∂X are the deformation steering angles of the unit longitudinal force of the front and rear suspension systems, respectively. ∂αf/∂X is the deformation steering angle of the unit steering torque.Figure14 shows the change of vehicle’s toe angle with longitudinal force. In Figure 14, when the longitudinal force is positive, it is the driving force, and when the longitudinal force is negative, it is the braking force. As the braking force increases, the toe angle of the front suspension increases in the negative direction, and the toe angle of the rear suspension increases in the positive direction. As the driving force increases, the toe angle of the front suspension increases positively, and the toe angle of the rear suspension increases negatively.Figure 14
The curve of vehicle’s toe angle with longitudinal force.Figure15 shows the curve of the front and rear suspensions’ toe angle with the aligning torque. In Figure 15, as the wheel’s aligning force increases, the toe angle of the front suspension and the rear suspension decreases correspondingly.Figure 15
The change of the car’s toe angle with the aligning force.


## 7. Lateral Force of the Car’s Tires Generated by Suspension

For vehicles in motion, the front and rear wheels are subjected to lateral forces resulting from the deformation of the tire. As mentioned earlier, the suspension can cause changes in the side slip stiffness and steering angle of the tire, which can change the tire force.The roll of the body will cause the wheel’s camber angle to change. The camber angle of the wheel will produce lateral force acting on the tire. At the same time, for a moving vehicle, the front and rear wheels are subjected to lateral forces caused by tire deformation. Taking into account the roll of the vehicle body, the forceFy1 caused by suspension that acts on the side of the tires is expressed as follows:(30)Fy1=Yf1+Yf2+Yr1+Yr2+2Ycf+2Ycr,(31)2Ycf=−2Zcf∂σf∂ϕϕ,2Ycr=−2Zcr∂σr∂ϕϕ,where Yf1, Yf2, Yr1, and Yr2 are the lateral force caused by the deformation of the vehicle’s tires. Ycf and Ycr are the lateral camber thrust on the front and rear tires of the car. Zcf and Zcr are the lateral camber thrust coefficients of the front and rear tires, respectively. σf and σr are the camber angles of front and rear tires.After performing the opposite travel simulation on the vehicle’s front and rear suspension’s multibody dynamics models established in Figures4 and 5, the curve of the wheel’s camber angle with wheel’s displacement is shown in Figure 16 is obtained [17]. In Figure 16, when the wheels bounce down, the camber angles of the front and rear suspension wheels (the inner wheel when turning) change to the positive camber direction. When the wheels jump up, the camber angles of the wheels (outer wheels when turning) of the front suspension and the rear suspension both change in the direction of negative camber. In terms of adhesion, whether it is the front suspension’s or the rear suspension’s outer wheel camber angle changing to the negative camber direction, the tire adhesion will be improved. From the perspective of understeer characteristics of the car, when the car is turning, the camber angle of the outer wheels of the front suspension changes to the negative camber direction, which will reduce the understeer characteristic of the car, and the camber angle of the outer wheels of the rear suspension changes to the negative camber direction, which will increase the understeer characteristics of the car.Figure 16
Curves of the suspension’s camber angle with wheel’s displacement.From equations (23), (26), and (27), the front wheel’s steering force caused by the change of toe angle is(32)Yf1+Yf2=Kf1αf1+Kf2αf2=2Kf01−∂Kf∂WΔW2Kf0−122XfμWf2×∂αf∂ztfhW2Kϕy¨−∂αf∂Tξ+αcW2Kyx¨bWly¨−2Kf0∂αf∂zahWKθx¨−∂αf∂XαcW2x¨∂Kf∂WΔWfKf0≈2Kf∗∂αf∂ztfhW2Kϕ−∂αf∂Tξ+αcW2Kyx¨bWly¨−∂αf∂zlfhWKθx¨−∂αf∂XαcW2x¨∂Kf∂WΔWfKf0=2Kf∗af.Similarly, for the rear wheels,(33)Yr1+Yr2=Kr1αr1+Kr1αr2≈2Kr∗∂αr∂zbhWKθy¨+∂αr∂zbhW2Kθx¨+∂αr∂X1−αc2Wx¨∂Kr∂WΔWrKr0=2Kr∗αr.In the above equations, tiny quantities of order 2 and above are omitted for small longitudinal and lateral accelerations.In the above derivation process of equations (32) and (33), the changes of toe angles αf1,αf2 and deformation steering angle αr1,αr2 of the front and rear wheels are replaced by a single equivalent angle αf,αr. Sum up, the lateral force caused by suspension of vehicle is(34)Fy1=2Kf∗αf+2Kr∗αr+2Ycf+2Ycr.


## 8. Establishment of Simplified Car Dynamics Model

The car itself is a complex nonlinear dynamic system, it is difficult to describe it accurately with formulas. However, in order to reflect the performance of the car under some special working conditions, after omitting some minor parts of the complex automobile system, a simplified 5-DOF vehicle dynamic model is established. The model’s degrees of freedom include three angular motions of the sprung mass (roll, pitch, and yaw) and two planar motions of the unsprung mass (longitudinal and lateral) [1, 18].


### 8.1. Kinematics Analysis of Vehicle

In the vehicle’s course of motion, the lateral acceleration and longitudinal acceleration will change, and the sprung mass will move relative to unsprung mass. The sprung mass centerOs and unsprung mass center Ou of the car are shown in Figure 17.Figure 17
The center of mass of body’s different parts.In Figure17, hp is the height from pitch center Op to the mass center of sprung mass Os; e is the height from the unsprung’s mass center to the ground; C1 is the distance from the unsprung’s mass center to the pitch center; and C2 is the distance from the mass center of sprung to the pitch center.As shown in Figure18, α is side slip angle of front wheels. ψ is the yaw angle of the car’s unsprung mass. δf is the steering angle of the car’s front wheels. αf is the side slip angle of the car’s front tires. β is the side slip angle of the car’s center of mass.Figure 18
Side slip angle of front wheels.The longitudinal speedu and lateral speed v of the unsprung mass can be calculated by(35)X˙=ucosψ−vsinψ,Y˙=usinψ+vcosψ.The longitudinal and lateral accelerations of the unsprung mass can be calculated by(36)aOux=X¨cosψ+Y¨sinψ=u˙−ψ˙v,aOuy=−X¨sinψ+Y¨cosψ=v˙+ψ˙u.The acceleration vector of vehicle’s sprung massa⟶Os is obtained by summing the vector of the vehicle’s unsprung mass acceleration a⟶Ou, the relative acceleration a⟶Os/Ou, the Coriolis acceleration a⟶C, and the implicated acceleration a⇀Trac:(37)a⟶Os=a⇀Ou+a⇀Os/Ou+a⇀C+a⇀Trac=Aii^+Ajj^+Akk^.Among them,(38)a⟶C=2ψ˙k^×v⟶Os/Ou;a⟶Trac=a⟶Tracn+a⟶Tract,a⟶Tracn=ψ˙k^×ψ˙k^×r⟶Os/Ou;a⟶Tract=ψ¨k^×r⟶Os/Ou,(39)Ai=−C2θ˙2cosθ−C2θ¨sinθ−ψ˙2C1+C2cosθ+hpcosϕsinθ+ψ¨hpsinϕ+2ψ˙hpϕ˙cosϕ−hpϕ˙2cosϕ+ϕ¨sinϕsinθ+2θ˙ϕ˙sinϕcosθ+θ˙2cosϕsinθ−θ¨cosθcosϕ+u˙−vψ˙,Aj=hpϕ˙2sinϕ−ϕ¨cosϕ+ψ˙2sinϕ+ψ˙C1+C2cosθ+hpcosϕsinθ−2ψ˙C2θ˙sinθ+hpϕ˙sinϕsinθ−hpθ˙cosϕcosθ+v˙+uψ˙,Ak=−hpθ˙2cosθcosϕ+θ¨sinθcosϕ−2θϕ˙sinθsinϕ+ϕ¨cosθsinϕC2θ˙2cosθ+C2θ¨sinθ.


### 8.2. Dynamics Analysis of Vehicle

On the basis of the lateral stiffness, lateral force, and automobile’s kinematics analysis, the force of the vehicle can be obtained.In equation (39), Ai,Aj,Ak are the accelerations of the sprung mass in the three directions x,y,z. r⇀Os/Ou is the distance between the sprung mass center and the unsprung mass center. θ is the angle of the car’s sprung mass rotates relative to the unsprung mass around the x-axis, also known as the pitch angle. ϕ is the angle by which the sprung mass of the car rotates relative to the unsprung mass about the y-axis, also known as the roll angle. From these acceleration components, the inertial force of the vehicle’s body and the inertial force of the unsprung mass can be expressed as(40)Fs=−msAii+Ajj+Akk,Fu=−muu˙−ψ˙vi+v˙+ψ˙uj,where ms is the sprung mass of the car. mu is the unsprung mass of the car. The longitudinal force, lateral force, and inertia force of the wheel are shown in Figure 19.Figure 19
Model of vehicle subjected to external and inertial forces.In Figure19, Fiz is the vertical force on the tires. Fy is the lateral force, and Xit is the longitudinal force on the tires, where i=fR,fL,rR,rL. Op is pitch center, and Os is the mass center of sprung mass. Considering the external force of the entire vehicle, the equation of motion can be summarized as [19](41)muvdβdt−dψdt+msAi+FfLy+FfRysinδf+XfLt+XfRt−XrLt+XrRt=0,muudβdt+dψdt+msAj−FfLy+FfRycosδf−FrLy−FrRy+XfLt+XfRtsinδf=0,Ixd2ϕdt2−Ixzd2ψdt2+msAkhpsinϕ−mshpgsinϕ+12tfFfLz−FfRz+12trFrLz−FrRz+msAje−C2cosθ+hcosθcosϕ=0,Iyd2θdt2+msAkC1+C2cosθ+hpcosϕsinθ−msAie−C2cosθ+hpcosθcosϕ−mshpgsinθsinϕ+aFfLz+FfRz−bFrLz+FrRz=0,Izd2ψdt2−Ixzd2ϕdt2=msAihpsinϕ+msAjC1+C2cosθ+hpcosϕsinθ−12tfFfRy−FfLysinδf−bFrRy+FrLy−12tfXfRt−XfLtcosδf−aXfRt+XfLtsinδf+aFfRy+FfLycosδf+12trXrRt−XrLt.


### 8.3. Simplified Vehicle Dynamics Model

Because this paper uses the sine-swept steering input simulation, the longitudinal speed of the car remains the same in the whole process of the vehicle’s movement. So vehicle’s pitch degree of freedom is omitted.If the vehicle’s speed remains constant, the vehicle’s longitudinal forcesXfLt, XfRt, XrLt, and XrRt are ignored, and the car’s roll steering and roll camber are only considered. So the longitudinal acceleration x¨ is zero, and then the model can be simplified into a 3-DOF vehicle dynamics model including roll degrees of freedom, lateral degrees of freedom, and yaw degrees of freedom. Combining equations (32) and (33), the sum of the lateral force Fy1 caused by the suspension and the lateral force Fy2 caused by the tire is the total lateral force acting on the vehicle’s tires:(42)∑Fy=Fy1+Fy2=2Yf+2Yr+2Ycf+2Ycr+Fy2=2Kf∗∂αf∂ztfhW2Kϕ−∂αf∂TbWlξy¨−∂αf∂TbWlαcW2Kyy¨x¨−∂αf∂zlfhWKθx¨−∂αf∂XαcW2x¨∂Kf∂WΔWfKf0+2Kr∗∂αr∂zbhWKθy¨+∂αr∂zbhW2Kθx¨+∂αr∂X1−αc2Wx¨∂Kr∂WΔWrKr0−2Kcf∂σf∂ϕϕ−2Kcr∂σr∂ϕϕ+2Kf∗δ−β−auψ˙+2Kr∗β+buψ˙.Due to only the car’s roll steering and roll camber are considered, the car’s aligning torque is ignored, combining to equation (24), and then equation (42) becomes:(43)∑Fy=Fy1+Fy2=2Yf+2Yr+2Ycf+2Ycr+Fy2=2Kf0∂αf∂ϕϕs+2Kr0∂αr∂ϕϕ−2Kcf∂σf∂ϕϕ−2Kcr∂σr∂ϕϕ+2Kf0δ−β−auψ˙+2Kr0β+buψ˙.Combining equations (24) and (43), the dynamic equation of the model is as follows:(44)ms+muudβdt+2Kf0+Kr0β+ms+muu+2aKf0−bKr0uψ˙−mshsd2ϕdt2−2Yϕϕ=2Kfδ,2aKf0−bKr0β+Izdψ˙dt+2a2Kf0−b2Kr0vψ˙−Ixzd2ϕdt2−2Nϕϕ=2aKf0δ,−mshsudβdt−Ixzdψ˙dt−mshsuψ˙+Ixd2ϕdt2+Cϕdϕdt+Kϕ−msghsϕ=0.Among them,(45)Yϕ=∂αf∂ϕKf+∂αr∂ϕKr−∂σf∂ϕKcf+∂σr∂ϕKcr,Nϕ=∂αf∂ϕaKf−∂αr∂ϕbKr−∂σf∂ϕaKcf−∂σr∂ϕbKcr.Performing Laplace transformation on equations (44), the following equations represented by a matrix are obtained:(46)ms+muus+2Kf0+Kr0ms+muu+2aKf0−bKr0u−mshss2−2Yϕ2aKf0−bKr0Izs+2a2Kf0+b2Kr0u−Ixzs2−2Nϕ−mshsus−Ixzs−mshsuIxs2+Cϕs+Kϕ−msghsβψ˙ϕ=2Kf0δs2aKf0δs0.The transfer function can be obtained from equation (46):(47)βsδs=2Kf0ms+muu+2aKf0−bKr0/u−mshss2−2Yϕ2aKf0Izs+2a2Kf0+b2Kr0/u−Ixzs2−2Nϕ0−Ixzs−mshsvIxs2+Cϕs+Kϕ−msghsms+muus+2Kf0+Kr0ms+muu+2aKf0−bKr0/u−mshss2−2Yϕ2aKf0−bKr0Izs+2a2Kf0+b2Kr0/u−Ixzs2−2Nϕ−mshsus−Ixzs−mshsuIxs2+Cϕs+Kϕ−msghs,(48)ψ˙sδs=ms+muus+2Kf0+Kr02Kf0−mshss2−2Yϕ2aKf0−bKr02aKf0−Ixzs2−2Nϕ−mshsus0Ixs2+Cϕs+Kϕ−msghsms+muus+2Kf0+Kr0ms+muu+2aKf0−bKr0/u−mshss2−2Yϕ2aKf0−bKr0Izs+2a2Kf0+b2Kr0/u−Ixzs2−2Nϕ−mshsus−Ixzs−mshsuIxs2+Cϕs+Kϕ−msghs,(49)ϕsδs=ms+muus+2Kf0+Kr0ms+muu+2aKf0−bKr0/u2Kf02aKf0−bKr0Izs+2a2Kf0+b2Kr/u2aKf0−mshsus−Ixzs−mshsu0ms+muus+2Kf0+Kr0ms+muu+2aKf0−bKr0/u−mshss2−2Yϕ2aKf0−bKr0Izs+2a2Kf0+b2Kr0/u−Ixzs2−2Nϕ−mshsus−Ixzs−mshsuIxs2+Cϕs+Kϕ−msghs,where Kϕf is the roll stiffness of the front suspension and Kϕr is the roll stiffness of the rear suspension. Kϕf includes the front suspension’s spring stiffness Kϕfs and the antiroll bar’s stiffness Kϕfw. The stiffness Kϕr includes the rear suspension’s spring stiffness Kϕrs and the stiffness Kϕrw of the rear antiroll bar. Substitute equations (3) and (4) into equations (47)–(49), the transfer function of response indicators that couple stiffness of the spring and the antiroll bar relative to the steering wheel angle is obtained.


## 8.1. Kinematics Analysis of Vehicle

In the vehicle’s course of motion, the lateral acceleration and longitudinal acceleration will change, and the sprung mass will move relative to unsprung mass. The sprung mass centerOs and unsprung mass center Ou of the car are shown in Figure 17.Figure 17
The center of mass of body’s different parts.In Figure17, hp is the height from pitch center Op to the mass center of sprung mass Os; e is the height from the unsprung’s mass center to the ground; C1 is the distance from the unsprung’s mass center to the pitch center; and C2 is the distance from the mass center of sprung to the pitch center.As shown in Figure18, α is side slip angle of front wheels. ψ is the yaw angle of the car’s unsprung mass. δf is the steering angle of the car’s front wheels. αf is the side slip angle of the car’s front tires. β is the side slip angle of the car’s center of mass.Figure 18
Side slip angle of front wheels.The longitudinal speedu and lateral speed v of the unsprung mass can be calculated by(35)X˙=ucosψ−vsinψ,Y˙=usinψ+vcosψ.The longitudinal and lateral accelerations of the unsprung mass can be calculated by(36)aOux=X¨cosψ+Y¨sinψ=u˙−ψ˙v,aOuy=−X¨sinψ+Y¨cosψ=v˙+ψ˙u.The acceleration vector of vehicle’s sprung massa⟶Os is obtained by summing the vector of the vehicle’s unsprung mass acceleration a⟶Ou, the relative acceleration a⟶Os/Ou, the Coriolis acceleration a⟶C, and the implicated acceleration a⇀Trac:(37)a⟶Os=a⇀Ou+a⇀Os/Ou+a⇀C+a⇀Trac=Aii^+Ajj^+Akk^.Among them,(38)a⟶C=2ψ˙k^×v⟶Os/Ou;a⟶Trac=a⟶Tracn+a⟶Tract,a⟶Tracn=ψ˙k^×ψ˙k^×r⟶Os/Ou;a⟶Tract=ψ¨k^×r⟶Os/Ou,(39)Ai=−C2θ˙2cosθ−C2θ¨sinθ−ψ˙2C1+C2cosθ+hpcosϕsinθ+ψ¨hpsinϕ+2ψ˙hpϕ˙cosϕ−hpϕ˙2cosϕ+ϕ¨sinϕsinθ+2θ˙ϕ˙sinϕcosθ+θ˙2cosϕsinθ−θ¨cosθcosϕ+u˙−vψ˙,Aj=hpϕ˙2sinϕ−ϕ¨cosϕ+ψ˙2sinϕ+ψ˙C1+C2cosθ+hpcosϕsinθ−2ψ˙C2θ˙sinθ+hpϕ˙sinϕsinθ−hpθ˙cosϕcosθ+v˙+uψ˙,Ak=−hpθ˙2cosθcosϕ+θ¨sinθcosϕ−2θϕ˙sinθsinϕ+ϕ¨cosθsinϕC2θ˙2cosθ+C2θ¨sinθ.


## 8.2. Dynamics Analysis of Vehicle

On the basis of the lateral stiffness, lateral force, and automobile’s kinematics analysis, the force of the vehicle can be obtained.In equation (39), Ai,Aj,Ak are the accelerations of the sprung mass in the three directions x,y,z. r⇀Os/Ou is the distance between the sprung mass center and the unsprung mass center. θ is the angle of the car’s sprung mass rotates relative to the unsprung mass around the x-axis, also known as the pitch angle. ϕ is the angle by which the sprung mass of the car rotates relative to the unsprung mass about the y-axis, also known as the roll angle. From these acceleration components, the inertial force of the vehicle’s body and the inertial force of the unsprung mass can be expressed as(40)Fs=−msAii+Ajj+Akk,Fu=−muu˙−ψ˙vi+v˙+ψ˙uj,where ms is the sprung mass of the car. mu is the unsprung mass of the car. The longitudinal force, lateral force, and inertia force of the wheel are shown in Figure 19.Figure 19
Model of vehicle subjected to external and inertial forces.In Figure19, Fiz is the vertical force on the tires. Fy is the lateral force, and Xit is the longitudinal force on the tires, where i=fR,fL,rR,rL. Op is pitch center, and Os is the mass center of sprung mass. Considering the external force of the entire vehicle, the equation of motion can be summarized as [19](41)muvdβdt−dψdt+msAi+FfLy+FfRysinδf+XfLt+XfRt−XrLt+XrRt=0,muudβdt+dψdt+msAj−FfLy+FfRycosδf−FrLy−FrRy+XfLt+XfRtsinδf=0,Ixd2ϕdt2−Ixzd2ψdt2+msAkhpsinϕ−mshpgsinϕ+12tfFfLz−FfRz+12trFrLz−FrRz+msAje−C2cosθ+hcosθcosϕ=0,Iyd2θdt2+msAkC1+C2cosθ+hpcosϕsinθ−msAie−C2cosθ+hpcosθcosϕ−mshpgsinθsinϕ+aFfLz+FfRz−bFrLz+FrRz=0,Izd2ψdt2−Ixzd2ϕdt2=msAihpsinϕ+msAjC1+C2cosθ+hpcosϕsinθ−12tfFfRy−FfLysinδf−bFrRy+FrLy−12tfXfRt−XfLtcosδf−aXfRt+XfLtsinδf+aFfRy+FfLycosδf+12trXrRt−XrLt.


## 8.3. Simplified Vehicle Dynamics Model

Because this paper uses the sine-swept steering input simulation, the longitudinal speed of the car remains the same in the whole process of the vehicle’s movement. So vehicle’s pitch degree of freedom is omitted.If the vehicle’s speed remains constant, the vehicle’s longitudinal forcesXfLt, XfRt, XrLt, and XrRt are ignored, and the car’s roll steering and roll camber are only considered. So the longitudinal acceleration x¨ is zero, and then the model can be simplified into a 3-DOF vehicle dynamics model including roll degrees of freedom, lateral degrees of freedom, and yaw degrees of freedom. Combining equations (32) and (33), the sum of the lateral force Fy1 caused by the suspension and the lateral force Fy2 caused by the tire is the total lateral force acting on the vehicle’s tires:(42)∑Fy=Fy1+Fy2=2Yf+2Yr+2Ycf+2Ycr+Fy2=2Kf∗∂αf∂ztfhW2Kϕ−∂αf∂TbWlξy¨−∂αf∂TbWlαcW2Kyy¨x¨−∂αf∂zlfhWKθx¨−∂αf∂XαcW2x¨∂Kf∂WΔWfKf0+2Kr∗∂αr∂zbhWKθy¨+∂αr∂zbhW2Kθx¨+∂αr∂X1−αc2Wx¨∂Kr∂WΔWrKr0−2Kcf∂σf∂ϕϕ−2Kcr∂σr∂ϕϕ+2Kf∗δ−β−auψ˙+2Kr∗β+buψ˙.Due to only the car’s roll steering and roll camber are considered, the car’s aligning torque is ignored, combining to equation (24), and then equation (42) becomes:(43)∑Fy=Fy1+Fy2=2Yf+2Yr+2Ycf+2Ycr+Fy2=2Kf0∂αf∂ϕϕs+2Kr0∂αr∂ϕϕ−2Kcf∂σf∂ϕϕ−2Kcr∂σr∂ϕϕ+2Kf0δ−β−auψ˙+2Kr0β+buψ˙.Combining equations (24) and (43), the dynamic equation of the model is as follows:(44)ms+muudβdt+2Kf0+Kr0β+ms+muu+2aKf0−bKr0uψ˙−mshsd2ϕdt2−2Yϕϕ=2Kfδ,2aKf0−bKr0β+Izdψ˙dt+2a2Kf0−b2Kr0vψ˙−Ixzd2ϕdt2−2Nϕϕ=2aKf0δ,−mshsudβdt−Ixzdψ˙dt−mshsuψ˙+Ixd2ϕdt2+Cϕdϕdt+Kϕ−msghsϕ=0.Among them,(45)Yϕ=∂αf∂ϕKf+∂αr∂ϕKr−∂σf∂ϕKcf+∂σr∂ϕKcr,Nϕ=∂αf∂ϕaKf−∂αr∂ϕbKr−∂σf∂ϕaKcf−∂σr∂ϕbKcr.Performing Laplace transformation on equations (44), the following equations represented by a matrix are obtained:(46)ms+muus+2Kf0+Kr0ms+muu+2aKf0−bKr0u−mshss2−2Yϕ2aKf0−bKr0Izs+2a2Kf0+b2Kr0u−Ixzs2−2Nϕ−mshsus−Ixzs−mshsuIxs2+Cϕs+Kϕ−msghsβψ˙ϕ=2Kf0δs2aKf0δs0.The transfer function can be obtained from equation (46):(47)βsδs=2Kf0ms+muu+2aKf0−bKr0/u−mshss2−2Yϕ2aKf0Izs+2a2Kf0+b2Kr0/u−Ixzs2−2Nϕ0−Ixzs−mshsvIxs2+Cϕs+Kϕ−msghsms+muus+2Kf0+Kr0ms+muu+2aKf0−bKr0/u−mshss2−2Yϕ2aKf0−bKr0Izs+2a2Kf0+b2Kr0/u−Ixzs2−2Nϕ−mshsus−Ixzs−mshsuIxs2+Cϕs+Kϕ−msghs,(48)ψ˙sδs=ms+muus+2Kf0+Kr02Kf0−mshss2−2Yϕ2aKf0−bKr02aKf0−Ixzs2−2Nϕ−mshsus0Ixs2+Cϕs+Kϕ−msghsms+muus+2Kf0+Kr0ms+muu+2aKf0−bKr0/u−mshss2−2Yϕ2aKf0−bKr0Izs+2a2Kf0+b2Kr0/u−Ixzs2−2Nϕ−mshsus−Ixzs−mshsuIxs2+Cϕs+Kϕ−msghs,(49)ϕsδs=ms+muus+2Kf0+Kr0ms+muu+2aKf0−bKr0/u2Kf02aKf0−bKr0Izs+2a2Kf0+b2Kr/u2aKf0−mshsus−Ixzs−mshsu0ms+muus+2Kf0+Kr0ms+muu+2aKf0−bKr0/u−mshss2−2Yϕ2aKf0−bKr0Izs+2a2Kf0+b2Kr0/u−Ixzs2−2Nϕ−mshsus−Ixzs−mshsuIxs2+Cϕs+Kϕ−msghs,where Kϕf is the roll stiffness of the front suspension and Kϕr is the roll stiffness of the rear suspension. Kϕf includes the front suspension’s spring stiffness Kϕfs and the antiroll bar’s stiffness Kϕfw. The stiffness Kϕr includes the rear suspension’s spring stiffness Kϕrs and the stiffness Kϕrw of the rear antiroll bar. Substitute equations (3) and (4) into equations (47)–(49), the transfer function of response indicators that couple stiffness of the spring and the antiroll bar relative to the steering wheel angle is obtained.


## 9. Build the Vehicle’s Dynamics Model and Multitarget Optimization

The impact of the suspension on the tire’s force can be known from the previous analysis. The load transfer will cause the tire’s side slip stiffness to change, and the longitudinal force and lateral force will cause the toe angle change and deformation steering of the tire. The influence of suspension on the side slip stiffness and the tire’s steering angle can be equivalent to the lateral force generated by the deformation of the automobile’s tire.


### 9.1. Build the Car’s Dynamics Model

In order to analyze and simulate the vehicle’s handling stability, based on the above analysis of vehicle’s dynamics, a vehicle’s dynamic model was established in ADAMS/CAR software. The established vehicle’s dynamic model includes a front suspension system, a rear suspension system, a steering system, and a braking system.The tire’s Magic formula was used in building the tire’s dynamics model. The lateral force of a car’s tire can be obtained by the Magic formula:(50)Fijy=FijzsinDttan−1Bitαij−EtBitαij−tan−1Bitαij+γijε1Fijz+ε2,where i and j are the labels for the wheels, i=f,r, j=R,L. The coefficients Dt, Bit, and Et need to be measured through actual experiment and have no practical physical significance. The side slip angle αfj of the front wheels and the side slip angle αrj of the rear wheels can be calculated by the following formulas:(51)αfj=tan−12v+2aψ˙2u+β1tfψ˙−δf,αrj=tan−12v−2bψ˙2u+β1trψ˙,j=R,L,β1=+1,j=L,−1,j=R.The longitudinal force of a car’s tires depends on the vehicle’s forward acceleration and whether the tire is a driving wheel. Considering the longitudinal force of the car and the force of the tire, the equation of motion of the car’s tire can be obtained.The motion equation of the driving wheel is(52)Irτ˙rj=Td−rXrjt−εFrjz,j=R,L.The motion equation of the slave wheel is(53)Ifτ˙fj=rXrjt−εFfjz,j=R,L.In equations (52) and (53), Fiz is the vertical force on the wheels, and Xit is the longitudinal force on the wheels, i=fR,fL,rR,rL. If and Ir are the moment inertia of front and rear wheels. τ is the angular velocity of wheels. r is wheels’ radius. Td is the driving torque. ε is rolling resistance of wheels. The vertical displacement dit between the connection point of the spring and the damper caused by the roll and pitch movement of the car and the vertical speed vij of the sprung mass can be calculated by the following formulas:(54)dij=−p1sinθ+12tip2sinϕcosθ+hzcosθcosϕ−1,vij=12tip2ϕ˙cosθ−p2p1θ˙,where(55)p1=a,i=f,b,i=r,p2=−1,j=R,1,j=L,The vertical load on each wheel is(56)Fijz=ms+mugξ1a+ξ2ba+b+dijki+vijci,where(57)ξ1=0,i=f,1,i=r,ξ2=1,j=R,0,i=r.For pure rotation of each wheel, its angular velocity can be expressed as(58)τij=1Ru+β2ψ˙ti22+v+β1aψ˙2,where(59)β1=1,i=f,−1,i=r,β2=1,j=L,−1,j=R.The parameters of the vehicle’s dynamics model finally established are shown in Table1.Table 1
Vehicle’s parameters.
ParameterValuems (kg)1067mu (kg)595A (mm)1452.6B (mm)1249.5tf (mm)1508.3tr (mm)1487.8mf (kg)891.89ktf (N/mm)210kϕfs (N/mm)43kϕrs (N/mm32kϕwf (N/mm)29400kϕwr (N/mm)12800hs (mm)330l (mm)2702.1mr (kg)770.11ktr (N/mm)240This article uses the sine-swept steering input simulation method to evaluate the handling stability of the car; that is, the steering wheel angle’s input waveform is changed to a sine wave. Using sine wave input with different frequencies, the amplitude ratio of the output to the input at different frequencies can be obtained, that is, the amplitude-frequency characteristics of the car. At the same time, the phase difference between the output and the input at different frequencies can be obtained, that is, the phase frequency characteristics of the car. The simulation data of the sine-swept steering input simulation in this paper are shown in Table2.Table 2
The simulation parameters.
ParameterValueTime (s)20Steps2048Velocity (km/h)100Max steer wheel angle (°)24.5Initial frequency (Hz)0.2Max frequency (Hz)3.8Frequency rate (Hz/s)0.2Start time (s)2


### 9.2. Set Up Multiple Target Optimization Functions

In equations (48)–(50), make s=jω and j=−1, the frequency response characteristic function of the car can be obtained, including amplitude-frequency characteristics and phase frequency characteristics. ψ˙/δω, ϕ/δω, and β/δω are the gain of yaw rate, roll angle, and side slip angle relative to the steering wheel angle in the amplitude-frequency characteristic function. ρψω, ρϕω, and ρβω are the corresponding phase lag angles. The frequency corresponding to the maximum gain of the yaw rate relative to the steering wheel angle is the resonance frequency of the car. The phase lag angle can be converted into the delay time by the following formula:(60)t=ρ+180360⋅ωorρ−180360⋅ω.Because the frequency of driver rotates the steering wheel is about 0.5 Hz. Based on the theoretical analysis of the above part, in the sine-swept steering input simulation, the yaw rate gain relative to the steering wheel angle at 0.5 Hz, the resonance frequencyfψ, roll angle gain relative to the steering wheel angle, the delay time t1 of lateral acceleration relative to the steering wheel angle, and the delay time t2 of yaw rate relative to lateral acceleration are used as optimized target. The stiffness of suspension’s spring and the antiroll bar are used as the optimized variables to do multiobjective optimization [20, 21].In order to achieve a better handling stability of the car, the yaw rate gain relative to the steering wheel angle, the roll angle gain relative to the steering wheel angle, the side slip angle gain relative to the steering wheel angle, the delay time of roll angle relative to the steering wheel angle, and the delay time of yaw rate relative to the lateral acceleration should be minimum, and at the same time, it is hoped that the resonance frequency of the yaw rate relative to the steering wheel angle is the maximum; that is, the objective function can be expressed as(61)f1=minψ˙δωω=0.5,f2=minϕδωω=0.5,f3=minβδωω=0.5,f4=maxfψ,f5=mint1,f6=mint2.


### 9.3. Multiobjective Optimization and Optimization Results

The Pareto solution is also called nondominated solution. When there are multiple goals, due to conflicts between goals, one solution is the best on a goal and maybe the worst on the other. While improving any objective function, these solutions will necessarily weaken at least one other objective function, which is called a nondominated solution or a Pareto solution. The set of optimal solutions of a set of objective functions is called the Pareto solution set. In 1986, Pareto introduced the concept of a multiobjective solution, which is defined as follows: Assume that any two solutions S1 and S2 for all targets, if S1 is better than S2, then we call S1 dominate S2; if S1 is not dominated by other solutions, S1 is called a nondominated solution, also called the Pareto solution. This paper adopts multiobjective genetic algorithm to optimize the handling stability of the car.After optimization calculation, the Pareto solution set of the vehicle’s resonance frequency and yaw rate gain value at 0.5 Hz is obtained as shown in Figure20. Figure 21 shows the Pareto solution set of resonance frequency and the delay time of lateral acceleration relative to the steering wheel angle [22].Figure 20
Pareto solution set of resonance frequency and yaw rate gain at 0.5 Hz.Figure 21
Pareto solution set of the delay time of lateral acceleration relative to the steering wheel angle and resonance frequency.Figure22 is a Pareto solution set of the delay time of yaw rate relative to the lateral acceleration and the resonance frequency. Figure 23 shows the Pareto solution set of the delay time of yaw rate relative to lateral acceleration (delay time 1) and the delay time of lateral acceleration relative to steering wheel angle (delay time 2).Figure 22
Pareto solution set of the delay time of yaw rate relative to lateral acceleration and resonance frequency.Figure 23
Pareto solution set of side slip angle gain relative to the yaw rate gain.Figure24 shows the Pareto solution set of the delay time of the lateral acceleration relative to the steering wheel angle and the yaw rate gain value. Figure 25 shows the Pareto solution set of side slip angle gain relative to the yaw rate gain.Figure 24
Pareto solution set of delay time of the lateral acceleration relative to the steering wheel angle and the yaw rate gain.Figure 25
Pareto solution set of side slip angle gain relative to the yaw rate gain.As shown in Figure26, the blue dots are the Pareto solution set of the side slip angle gain relative to the delay time of yaw rate relative to the lateral acceleration.Figure 26
Pareto solution set of the side slip angle gain relative to the delay time.In the Pareto solution set of Figures20–26, a relatively optimal solution is obtained after the optimization operation. The relatively optimal solution and two sets of solutions selected from the Pareto solution set are, respectively, brought into the vehicle model for sine-swept steering input simulation and then compared with the original vehicle’s simulation data.The simulation results obtained is shown in Table3. Compared with the original car, the increase or decrease percentages of spring stiffness, antiroll bar stiffness, and optimization targets at 0.5 Hz are shown in Table 4. In Table 3, solution 1 is a relatively optimal solution, and solution 2 and solution 3 are two sets of solutions selected from the Pareto solution set.Table 3
The results before and after optimization.
kϕfs (N/mm)kϕrs (N/mm)kϕfw (N/mm)kϕrw (N/mm)f1s−1f2s−1f3s−1f4Hzf5sf6sOriginal car432029400128000.480.090.060.610.200.14Solution 186322648791950.430.120.050.590.180.11Solution 251.638.435280153600.450.080.080.640.220.16Solution 334.425.623520102400.490.100.090.690.210.13Table 4
The increase or decrease percentage of the optimized parameters compared with the original car.
kϕfs (N/mm)kϕrs (N/mm)kϕfw (N/mm)kϕrw (N/mm)f1s−1f2s−1f3s−1f4Hzf5sf6sPercentage50%37.5%−10%−28%−10.4%−25%−17%−3.3%−10%−21.4%As shown in Table4, compared with the original car, the optimized car’s front suspension spring stiffness and rear suspension spring stiffness increased by 50% and 37.5%, respectively. And the front suspension antiroll bar stiffness and rear suspension antiroll bar stiffness decreased by 10% and 28%, respectively. The yaw rate gain is reduced by 10.4%, the roll angle gain is reduced by 25%, the side slip angle gain is reduced by 17%, the resonance frequency is reduced by 3.3%, the delay time of lateral acceleration is reduced by 10%, and the delay time of yaw is reduced by 21.4%.After simulating the 4 sets of data in Table3 with sine-swept steering input simulation, the comparison figures shown in Figures 27–33 were obtained. In Figures 27–33, Solution 1 is the relatively optimal solution, and Solution 2 and Solution 3 are the two sets of solutions selected from the Pareto solution set.Figure 27
Yaw rate gain relative to the steering wheel angle.Figure 28
Phase lag of yaw rate relative to the lateral acceleration.Figure 29
The delay time of yaw rate relative to lateral acceleration.Figure 30
Roll angle gain relative to the lateral acceleration.Figure 31
Delay time of lateral acceleration relative to steering wheel angle.Figure 32
Phase lag of yaw rate relative to the steering angle.Figure 33
Side slip angle gain relative to the steering wheel angle.Figure27 shows the curve of yaw rate gain relative to the steering wheel angle obtained after the four sets of data are brought into the model for simulation. In Figure 27, within 1 Hz, the difference between the curves is obvious. After optimization, the gain of yaw rate relative to the steering wheel at 0.5 Hz is smaller than the other three solutions, which meets the requirements of the optimization goal. Between 0 and 1 Hz, the gain value of the yaw rate relative to the steering wheel angle corresponding to the optimized scheme is smaller than the other three groups of schemes. The curve in Figure 27 has a small difference between 1 Hz and 4 Hz. Because the driver rarely exceeds the range of 1 Hz during the handling of the car, the change of the yaw rate gain value at high frequency is not important to be considered.Figure28 shows the phase lag of yaw rate relative to the lateral acceleration. Convert the phase lag in Figure 28 to the delay time of yaw rate relative to the lateral acceleration according to equation (60), and the delay time shown in Figure 29 can be obtained.Figure29 shows the delay time of the yaw rate relative to the lateral acceleration. In Figure 29, because it is a multiobjective optimization, it is difficult to make all optimization objectives reach the most ideal state. The delay time of yaw rate relative to the lateral acceleration at 0.5 Hz obtained after optimization is not much different from the other two sets of solutions. But its value is within the acceptable range. Between 0 and 1 Hz, the delay time of the yaw rate relative to the lateral acceleration corresponding to the optimized scheme was compared with the other three groups of schemes, and the result was the same as that at 0.5 Hz.Figure30 shows the roll angle gain relative to the lateral acceleration. In Figure 30, as the frequency increases, the curve of the roll angle gain relative to the lateral acceleration corresponding to the relatively optimal solution is smaller than the curve corresponding to the other three solutions. It shows that the obtained relatively optimal solution is beneficial to improve the handling stability of the car.Figure31 shows the delay time of lateral acceleration relative to steering wheel angle. In Figure 31, for the optimized curve, the delay time of lateral acceleration relative to the steering wheel angle at 0.5 Hz is slightly different from the curves of the other two sets of solutions, but its value is not the worst among the four sets of curves and within the acceptable range. Between 0 and 1 Hz, the delay time of the lateral acceleration relative to the steering wheel angle corresponding to the optimized scheme is compared with the other three schemes, and the result is the same as that at 0.5 Hz. As the frequency increases, the absolute value of the delay time gradually decreases. Since the driver’s manipulation of the car rarely reaches more than 1 Hz, the vehicle’s frequency characteristics at high frequencies are not considered.Figure32 shows the phase lag of yaw rate relative to the steering wheel angle. Obviously, the absolute value of the phase lag of the yaw rate relative to the steering wheel angle after the optimization is significantly reduced. This shows that the relative optimal solution of the car will make the car respond faster. This has an advantageous effect on the handling stability of the vehicle.In Figure33, within 1 Hz, the difference between the curves of the side slip angle gain relative to the steering wheel angle is relatively obvious, but as the frequency increases, the difference between the curves becomes smaller and smaller. The driving frequency of the general driver is generally maintained in the low frequency range, that is, within 1 Hz, so generally the car’s high frequency range is not considered. Moreover, the optimized vehicle’s side slip angle gain relative to the steering wheel angle at 0.5 Hz is smaller than that of the other three solutions, which is beneficial to improve the handling stability of the vehicle and meets the requirements of the optimization goal.


## 9.1. Build the Car’s Dynamics Model

In order to analyze and simulate the vehicle’s handling stability, based on the above analysis of vehicle’s dynamics, a vehicle’s dynamic model was established in ADAMS/CAR software. The established vehicle’s dynamic model includes a front suspension system, a rear suspension system, a steering system, and a braking system.The tire’s Magic formula was used in building the tire’s dynamics model. The lateral force of a car’s tire can be obtained by the Magic formula:(50)Fijy=FijzsinDttan−1Bitαij−EtBitαij−tan−1Bitαij+γijε1Fijz+ε2,where i and j are the labels for the wheels, i=f,r, j=R,L. The coefficients Dt, Bit, and Et need to be measured through actual experiment and have no practical physical significance. The side slip angle αfj of the front wheels and the side slip angle αrj of the rear wheels can be calculated by the following formulas:(51)αfj=tan−12v+2aψ˙2u+β1tfψ˙−δf,αrj=tan−12v−2bψ˙2u+β1trψ˙,j=R,L,β1=+1,j=L,−1,j=R.The longitudinal force of a car’s tires depends on the vehicle’s forward acceleration and whether the tire is a driving wheel. Considering the longitudinal force of the car and the force of the tire, the equation of motion of the car’s tire can be obtained.The motion equation of the driving wheel is(52)Irτ˙rj=Td−rXrjt−εFrjz,j=R,L.The motion equation of the slave wheel is(53)Ifτ˙fj=rXrjt−εFfjz,j=R,L.In equations (52) and (53), Fiz is the vertical force on the wheels, and Xit is the longitudinal force on the wheels, i=fR,fL,rR,rL. If and Ir are the moment inertia of front and rear wheels. τ is the angular velocity of wheels. r is wheels’ radius. Td is the driving torque. ε is rolling resistance of wheels. The vertical displacement dit between the connection point of the spring and the damper caused by the roll and pitch movement of the car and the vertical speed vij of the sprung mass can be calculated by the following formulas:(54)dij=−p1sinθ+12tip2sinϕcosθ+hzcosθcosϕ−1,vij=12tip2ϕ˙cosθ−p2p1θ˙,where(55)p1=a,i=f,b,i=r,p2=−1,j=R,1,j=L,The vertical load on each wheel is(56)Fijz=ms+mugξ1a+ξ2ba+b+dijki+vijci,where(57)ξ1=0,i=f,1,i=r,ξ2=1,j=R,0,i=r.For pure rotation of each wheel, its angular velocity can be expressed as(58)τij=1Ru+β2ψ˙ti22+v+β1aψ˙2,where(59)β1=1,i=f,−1,i=r,β2=1,j=L,−1,j=R.The parameters of the vehicle’s dynamics model finally established are shown in Table1.Table 1
Vehicle’s parameters.
ParameterValuems (kg)1067mu (kg)595A (mm)1452.6B (mm)1249.5tf (mm)1508.3tr (mm)1487.8mf (kg)891.89ktf (N/mm)210kϕfs (N/mm)43kϕrs (N/mm32kϕwf (N/mm)29400kϕwr (N/mm)12800hs (mm)330l (mm)2702.1mr (kg)770.11ktr (N/mm)240This article uses the sine-swept steering input simulation method to evaluate the handling stability of the car; that is, the steering wheel angle’s input waveform is changed to a sine wave. Using sine wave input with different frequencies, the amplitude ratio of the output to the input at different frequencies can be obtained, that is, the amplitude-frequency characteristics of the car. At the same time, the phase difference between the output and the input at different frequencies can be obtained, that is, the phase frequency characteristics of the car. The simulation data of the sine-swept steering input simulation in this paper are shown in Table2.Table 2
The simulation parameters.
ParameterValueTime (s)20Steps2048Velocity (km/h)100Max steer wheel angle (°)24.5Initial frequency (Hz)0.2Max frequency (Hz)3.8Frequency rate (Hz/s)0.2Start time (s)2


## 9.2. Set Up Multiple Target Optimization Functions

In equations (48)–(50), make s=jω and j=−1, the frequency response characteristic function of the car can be obtained, including amplitude-frequency characteristics and phase frequency characteristics. ψ˙/δω, ϕ/δω, and β/δω are the gain of yaw rate, roll angle, and side slip angle relative to the steering wheel angle in the amplitude-frequency characteristic function. ρψω, ρϕω, and ρβω are the corresponding phase lag angles. The frequency corresponding to the maximum gain of the yaw rate relative to the steering wheel angle is the resonance frequency of the car. The phase lag angle can be converted into the delay time by the following formula:(60)t=ρ+180360⋅ωorρ−180360⋅ω.Because the frequency of driver rotates the steering wheel is about 0.5 Hz. Based on the theoretical analysis of the above part, in the sine-swept steering input simulation, the yaw rate gain relative to the steering wheel angle at 0.5 Hz, the resonance frequencyfψ, roll angle gain relative to the steering wheel angle, the delay time t1 of lateral acceleration relative to the steering wheel angle, and the delay time t2 of yaw rate relative to lateral acceleration are used as optimized target. The stiffness of suspension’s spring and the antiroll bar are used as the optimized variables to do multiobjective optimization [20, 21].In order to achieve a better handling stability of the car, the yaw rate gain relative to the steering wheel angle, the roll angle gain relative to the steering wheel angle, the side slip angle gain relative to the steering wheel angle, the delay time of roll angle relative to the steering wheel angle, and the delay time of yaw rate relative to the lateral acceleration should be minimum, and at the same time, it is hoped that the resonance frequency of the yaw rate relative to the steering wheel angle is the maximum; that is, the objective function can be expressed as(61)f1=minψ˙δωω=0.5,f2=minϕδωω=0.5,f3=minβδωω=0.5,f4=maxfψ,f5=mint1,f6=mint2.


## 9.3. Multiobjective Optimization and Optimization Results

The Pareto solution is also called nondominated solution. When there are multiple goals, due to conflicts between goals, one solution is the best on a goal and maybe the worst on the other. While improving any objective function, these solutions will necessarily weaken at least one other objective function, which is called a nondominated solution or a Pareto solution. The set of optimal solutions of a set of objective functions is called the Pareto solution set. In 1986, Pareto introduced the concept of a multiobjective solution, which is defined as follows: Assume that any two solutions S1 and S2 for all targets, if S1 is better than S2, then we call S1 dominate S2; if S1 is not dominated by other solutions, S1 is called a nondominated solution, also called the Pareto solution. This paper adopts multiobjective genetic algorithm to optimize the handling stability of the car.After optimization calculation, the Pareto solution set of the vehicle’s resonance frequency and yaw rate gain value at 0.5 Hz is obtained as shown in Figure20. Figure 21 shows the Pareto solution set of resonance frequency and the delay time of lateral acceleration relative to the steering wheel angle [22].Figure 20
Pareto solution set of resonance frequency and yaw rate gain at 0.5 Hz.Figure 21
Pareto solution set of the delay time of lateral acceleration relative to the steering wheel angle and resonance frequency.Figure22 is a Pareto solution set of the delay time of yaw rate relative to the lateral acceleration and the resonance frequency. Figure 23 shows the Pareto solution set of the delay time of yaw rate relative to lateral acceleration (delay time 1) and the delay time of lateral acceleration relative to steering wheel angle (delay time 2).Figure 22
Pareto solution set of the delay time of yaw rate relative to lateral acceleration and resonance frequency.Figure 23
Pareto solution set of side slip angle gain relative to the yaw rate gain.Figure24 shows the Pareto solution set of the delay time of the lateral acceleration relative to the steering wheel angle and the yaw rate gain value. Figure 25 shows the Pareto solution set of side slip angle gain relative to the yaw rate gain.Figure 24
Pareto solution set of delay time of the lateral acceleration relative to the steering wheel angle and the yaw rate gain.Figure 25
Pareto solution set of side slip angle gain relative to the yaw rate gain.As shown in Figure26, the blue dots are the Pareto solution set of the side slip angle gain relative to the delay time of yaw rate relative to the lateral acceleration.Figure 26
Pareto solution set of the side slip angle gain relative to the delay time.In the Pareto solution set of Figures20–26, a relatively optimal solution is obtained after the optimization operation. The relatively optimal solution and two sets of solutions selected from the Pareto solution set are, respectively, brought into the vehicle model for sine-swept steering input simulation and then compared with the original vehicle’s simulation data.The simulation results obtained is shown in Table3. Compared with the original car, the increase or decrease percentages of spring stiffness, antiroll bar stiffness, and optimization targets at 0.5 Hz are shown in Table 4. In Table 3, solution 1 is a relatively optimal solution, and solution 2 and solution 3 are two sets of solutions selected from the Pareto solution set.Table 3
The results before and after optimization.
kϕfs (N/mm)kϕrs (N/mm)kϕfw (N/mm)kϕrw (N/mm)f1s−1f2s−1f3s−1f4Hzf5sf6sOriginal car432029400128000.480.090.060.610.200.14Solution 186322648791950.430.120.050.590.180.11Solution 251.638.435280153600.450.080.080.640.220.16Solution 334.425.623520102400.490.100.090.690.210.13Table 4
The increase or decrease percentage of the optimized parameters compared with the original car.
kϕfs (N/mm)kϕrs (N/mm)kϕfw (N/mm)kϕrw (N/mm)f1s−1f2s−1f3s−1f4Hzf5sf6sPercentage50%37.5%−10%−28%−10.4%−25%−17%−3.3%−10%−21.4%As shown in Table4, compared with the original car, the optimized car’s front suspension spring stiffness and rear suspension spring stiffness increased by 50% and 37.5%, respectively. And the front suspension antiroll bar stiffness and rear suspension antiroll bar stiffness decreased by 10% and 28%, respectively. The yaw rate gain is reduced by 10.4%, the roll angle gain is reduced by 25%, the side slip angle gain is reduced by 17%, the resonance frequency is reduced by 3.3%, the delay time of lateral acceleration is reduced by 10%, and the delay time of yaw is reduced by 21.4%.After simulating the 4 sets of data in Table3 with sine-swept steering input simulation, the comparison figures shown in Figures 27–33 were obtained. In Figures 27–33, Solution 1 is the relatively optimal solution, and Solution 2 and Solution 3 are the two sets of solutions selected from the Pareto solution set.Figure 27
Yaw rate gain relative to the steering wheel angle.Figure 28
Phase lag of yaw rate relative to the lateral acceleration.Figure 29
The delay time of yaw rate relative to lateral acceleration.Figure 30
Roll angle gain relative to the lateral acceleration.Figure 31
Delay time of lateral acceleration relative to steering wheel angle.Figure 32
Phase lag of yaw rate relative to the steering angle.Figure 33
Side slip angle gain relative to the steering wheel angle.Figure27 shows the curve of yaw rate gain relative to the steering wheel angle obtained after the four sets of data are brought into the model for simulation. In Figure 27, within 1 Hz, the difference between the curves is obvious. After optimization, the gain of yaw rate relative to the steering wheel at 0.5 Hz is smaller than the other three solutions, which meets the requirements of the optimization goal. Between 0 and 1 Hz, the gain value of the yaw rate relative to the steering wheel angle corresponding to the optimized scheme is smaller than the other three groups of schemes. The curve in Figure 27 has a small difference between 1 Hz and 4 Hz. Because the driver rarely exceeds the range of 1 Hz during the handling of the car, the change of the yaw rate gain value at high frequency is not important to be considered.Figure28 shows the phase lag of yaw rate relative to the lateral acceleration. Convert the phase lag in Figure 28 to the delay time of yaw rate relative to the lateral acceleration according to equation (60), and the delay time shown in Figure 29 can be obtained.Figure29 shows the delay time of the yaw rate relative to the lateral acceleration. In Figure 29, because it is a multiobjective optimization, it is difficult to make all optimization objectives reach the most ideal state. The delay time of yaw rate relative to the lateral acceleration at 0.5 Hz obtained after optimization is not much different from the other two sets of solutions. But its value is within the acceptable range. Between 0 and 1 Hz, the delay time of the yaw rate relative to the lateral acceleration corresponding to the optimized scheme was compared with the other three groups of schemes, and the result was the same as that at 0.5 Hz.Figure30 shows the roll angle gain relative to the lateral acceleration. In Figure 30, as the frequency increases, the curve of the roll angle gain relative to the lateral acceleration corresponding to the relatively optimal solution is smaller than the curve corresponding to the other three solutions. It shows that the obtained relatively optimal solution is beneficial to improve the handling stability of the car.Figure31 shows the delay time of lateral acceleration relative to steering wheel angle. In Figure 31, for the optimized curve, the delay time of lateral acceleration relative to the steering wheel angle at 0.5 Hz is slightly different from the curves of the other two sets of solutions, but its value is not the worst among the four sets of curves and within the acceptable range. Between 0 and 1 Hz, the delay time of the lateral acceleration relative to the steering wheel angle corresponding to the optimized scheme is compared with the other three schemes, and the result is the same as that at 0.5 Hz. As the frequency increases, the absolute value of the delay time gradually decreases. Since the driver’s manipulation of the car rarely reaches more than 1 Hz, the vehicle’s frequency characteristics at high frequencies are not considered.Figure32 shows the phase lag of yaw rate relative to the steering wheel angle. Obviously, the absolute value of the phase lag of the yaw rate relative to the steering wheel angle after the optimization is significantly reduced. This shows that the relative optimal solution of the car will make the car respond faster. This has an advantageous effect on the handling stability of the vehicle.In Figure33, within 1 Hz, the difference between the curves of the side slip angle gain relative to the steering wheel angle is relatively obvious, but as the frequency increases, the difference between the curves becomes smaller and smaller. The driving frequency of the general driver is generally maintained in the low frequency range, that is, within 1 Hz, so generally the car’s high frequency range is not considered. Moreover, the optimized vehicle’s side slip angle gain relative to the steering wheel angle at 0.5 Hz is smaller than that of the other three solutions, which is beneficial to improve the handling stability of the vehicle and meets the requirements of the optimization goal.


## 10. Conclusions

In this paper, the influence mechanism of spring stiffness and antiroll bar stiffness on suspension and vehicle performance is analyzed, the simplified vehicle dynamics model is established, and the transfer function coupling spring and antiroll bar is obtained. Finally, a multibody dynamics model considering the suspension is established, and a relative optimal matching method of spring stiffness and antiroll bar stiffness is obtained through NSGA-II optimization algorithm. The specific conclusions are as follows:(1)
Suspension roll stiffness and roll center will affect the load transfer, thus changing the tire’s side slip stiffness.According to a real vehicle object studied in this paper, the law of the roll stiffness of the front and rear suspension with the spring stiffness and the stiffness of the antiroll bar is obtained. For the rear multilink suspension of the car, the instantaneous axis of motion and roll center were found through complex geometric analysis, and the law of the front and rear suspension roll center changing with wheel jump was obtained.(2)
The change of tire side slip stiffness, tire rotation angle, and camber angle caused by suspension will lead to the change of tire force.Consider the two factors of the suspension causing the change of the tire rotation angle. One is the change of the suspension stroke caused by the vehicle’s roll and pitch movement, which leads to the change of the toe angle, and this is related to the suspension’s roll stiffness and pitch stiffness. On the contrary, the force received by the tire will cause a change in the rotation angle of the tire, that is deformation steering.The influence of the suspension on the side slip stiffness and tire’s rotation angle is equivalent to the lateral force generated by tire deformation plus the camber thrust generated by the camber angle to obtain the total lateral force experienced by the tire.(3)
The corresponding Pareto solution set and the relative optimal solution are obtained by using the NSGA-II algorithm to optimize. By comparison, the dynamic performance of the car is improved.When the spring stiffness of the front suspension is increased by 50%, the stiffness of the rear suspension antiroll bar is increased by 37.5%, the stiffness of the front antiroll bar is decreased by 10%, and the spring stiffness of the rear suspension is decreased by 28%, the car’s corresponding yaw rate gain, the roll angle gain and the side slip angle gain at 0.5 Hz are reduced, and the delay time of lateral acceleration and yaw rate at 0.5 Hz are within a reasonable and acceptable range. At other frequency points between 0 and 1 Hz, the results of the optimized scheme compared with the other three groups of schemes were the same as those at 0.5 Hz.(4)
The changing law of the response index under sine-swept input at low frequency and high frequency is as follows:When the frequency is below 1 Hz, the gain and delay time do not change according to a certain law with the increase of frequency. Between 1 Hz and 3 Hz, the yaw rate gain, side slip angle gain, delay time of yaw rate, and lateral acceleration decrease as the frequency increases, and the roll angle gain increases as the frequency increases.(5)
In this paper, the optimization is carried out for the frequency point of 0.5 Hz. Next, the optimization target will be increased to a frequency range, so as to make the optimization more comprehensive and broader and improve the optimization effect.



---

*Source: 8820108-2020-12-21.xml*