Zhou 2019 IOP Conf. Ser. Mater. Sci. Eng. 538 012068
Zhou 2019 IOP Conf. Ser. Mater. Sci. Eng. 538 012068
Zhou 2019 IOP Conf. Ser. Mater. Sci. Eng. 538 012068
Quan Zhou1,a, Puguang Yuan1, Tim Tudor2, Xudong Wang1, Zhangwei Guo1,
Yihui Huang1, Jingfeng Sun1, Boyun Wang3
1
Institute of vehicle suspension, Automotive Engineering Department, Wuhan
University of Technology, China
2
School of Engineering, University of Wales Trinity Saint David, Swansea, UK
3
Bellevue College, Bellevue, Washington, US
Abstract. With the rapid development of automobile technology and the improvement of
people's quality of life, people put forward higher requirements for vehicle ride comfort. Shock
absorption has become an important direction for automobile research, manufacturing and
optimization. The traditional air suspension airbag delimits the initial air pressure range based
on the bearing capacity of the airbag material, so the range of the pressure is inaccurate.
Considering optimizing ride comfort of the vehicle, this paper proposes a method to determine
the accurate initial air pressure range. In order to determine the accurate initial air pressure
range of the airbag, this paper needs to know the vehicle ride comfort under any air pressure.
The airbag spring stiffness is calculated using the airbag model in MATLAB environment. A
shock absorber with appropriate damping coefficient is selected by using the single mass
system model and the airbag spring stiffness. Suspension stiffness is calculated by constructing
a single wishbone beam independent suspension model and using the airbag spring stiffness.
The damping coefficient and suspension stiffness are introduced into the Simulink five-degree-
freedom vehicle model to obtain the acceleration signals of every part including the seat whose
RMS value of the acceleration signal is used to judge the ride comfort of the vehicle. Under the
specific simulation conditions in this paper, there is an accurate range (14~16kPa) in the rough
range (12~20kPa) of the initial air pressure of the airbag to make the vehicle ride smoothly.
1. Introduction
Vehicle ride comfort is a significant performance indicator to compete for advantages in market
competition. With the improvement of people's living standards, the problem of structural vibration
comfort has attracted more and more researchers' attention. The vehicle ride comfort directly affect the
experience of the occupants, and indirectly affect the vehicle's power and economy. The suspension
system plays a significant role in the shock absorption performance of the vehicle. In practical
applications, the air suspension system can meet the various requirements of the suspension design
and save energy and material, thus driving many scholars and engineers extensively research and
gradually promote it. Air suspension has been widely used in various vehicles. However, for the vast
majority of air suspension airbags on the market, only a rough initial air pressure range is given, while
Content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution
of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.
Published under licence by IOP Publishing Ltd 1
ICMMME2019 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 538 (2019) 012068 doi:10.1088/1757-899X/538/1/012068
the accurate initial air pressure range of various airbags are not given. This paper predicts that there is
a specific pressure value or accurate pressure range in the rough initial air pressure range to maximize
the vehicle ride comfort. This paper intends to find the accurate pressure value or an accurate range of
the initial air pressure of the airbag by simulation.
The initial air pressure of the airbag is the main factor determining the natural frequency of the air
spring, which determines the air spring stiffness, the air suspension stiffness and the damping
coefficient of the shock absorber. Xiong established the vibration model of the air suspension system,
and proposed the method to calculate the parameters of the air spring and air suspension and the
method to select shock absorber with appropriate parameters [1]. Li proposed the formula to calculate
the stiffness of three different shaped air springs and the stiffness of air suspensions matching different
shaped springs [2].
Meanwhile, the rough road surface is the root cause of the vibration of the vehicle, so it is necessary to
establish a mathematical model to describe road roughness. Xie used computer simulation of random
road surface excitation as input, and established the road surface roughness model by
MATLAB/Simulink simulation software for computer simulation [3]. Liu built a Simulink simulation
model for Class B Road based on pseudo white noise method [4]. Zheng used the probability and
statistics theory and the central limit theorem to establish a mathematical model of the uneven road
surface spectrum based on the power spectral density of the principle of harmonic superposition [5].
Zhang's research summarizes the triangular series method and the Fourier inverse transformation
method for road roughness [6].
In general, whole body vibration is the main form of vibration that causes physical injury to the driver
so a standard should be established to evaluate the ride comfort. This paper decides to use the impact
of the body’s vibration to evaluate the ride comfort. In order to evaluate the impact of the body's
vibration, the International Organization for Standardization (ISO) introduced the international
standard "Guidelines for the Evaluation of Human Body's Resistance to Whole Body Vibration" (ISO
2631-1978) [9] in 1978. In order to accurately evaluate the ride comfort of the vehicle, Zhang uses the
RMS value of the vertical weighted acceleration of the passenger to evaluate the index of the flatness.
[7] Zhang summarized two internationally accepted evaluation criteria for ride comfort and three
quantitative calculation methods for passenger comfort [6].
The weighted RMS acceleration is related to seat acceleration signal. According to the relationship
between weighted RMS acceleration and riding comfort (Table 1), it is necessary to establish a
suitable mathematical model to analyze the vibration of every part of the vehicle to get seat
acceleration signal. In order to accurately evaluate the ride comfort of the vehicle, Zhang constructed a
five-degree-freedom vehicle model of human-vehicle-road, and analyzed the interaction between
human and vehicle roads by transfer matrix method [7]. Xie performed mechanical analysis to obtain
the system vibration differential equation and built the four-degree-freedom vehicle model by
Simulink. [3] Liu, through the dynamic analysis of vehicles, established the eight-degree-freedom
vehicle model, calculated the eight-degree- freedom coupled vibration equation of the vehicle four-
point random excitation, and used the block matrix method to derive the transfer function of the
2
ICMMME2019 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 538 (2019) 012068 doi:10.1088/1757-899X/538/1/012068
vehicle vibration system. [8] Zhang's research summarizes the two-degree-freedom, the four-degree-
freedom, five-degree-freedom, the seven-degree-freedom, seventeen-degree-freedom vehicle model
for vibration analysis of each part [6].
The above literatures provide a theoretical basis and modeling basis for studying the relationship between initial
air pressure and vehicle ride comfort of air suspension airbags, but none of these literatures study the initial air
pressure of air suspension airbags for the purpose of maximizing vehicle ride comfort. Therefore, the purpose of
this paper is to find a method for determining accurate initial air pressure range of air suspension airbag based on
vehicle ride comfort simulation analysis
In order to ensure the best ride comfort of a vehicle under a specific road roughness level, this paper proposes a
method to determine the accurate initial air pressure range of the air suspension airbag. For a certain type of
vehicle using air suspension, the control variable method is used to change the initial air pressure of the air
suspension airbag, thereby changing the suspension stiffness and suspension damping coefficient. Vehicle
parameters and the road roughness excitation derived from the Simulink road roughness model are introduced to
five-degree-freedom vehicle Simulink model. Using the transfer matrix method to solve the vibration equation,
the acceleration signal of the vehicle seat under each initial air pressure of airbag can be obtained. Then use the
acceleration signal of the vehicle seat to find the seat acceleration RMS value. The RMS value of the seat
acceleration at each initial air pressure is compared to find an accurate range of the initial airbag pressure that
optimizes ride comfort. The method can be used to determine the optimal initial air pressure range of the air
suspension airbag of the vehicle under various road conditions and the vehicle can adjust the initial air pressure
of the air suspension in real time according to the road condition; The method can also be used to initialize the
air suspension airbag gas pressure during maintenance.
Figure 1. System model of optimal initial air pressure optimization method for airbags.
2. System description
The working process of the optimal initial air pressure optimization system for airbags proposed in this paper is
shown in the figure 1. First, an initial air pressure is required to set the air suspension airbag. Next, the
suspension stiffness and suspension damping coefficient can be calculated by inputting the initial air pressure of
the airbag and the parameters related to airbag and air suspension into the single mass system model and the
airbag model. Then, the suspension stiffness damping coefficient, vehicle parameters and road roughness
excitation obtained from the road model are input into the five-degree-freedom vehicle Simulink model to obtain
the signals of the vibration displacement, velocity and acceleration of each part of the vehicle. Then, the RMS
3
ICMMME2019 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 538 (2019) 012068 doi:10.1088/1757-899X/538/1/012068
value of the existing seat acceleration signal is obtained to evaluate the ride comfort of the vehicle as an output
and recorded. Finally, change the initial air pressure of the airbag and repeat the above steps to find the initial air
pressure of the airbag with the optimal ride comfort.
3. Mathematical model
In this paper, a five-degree-freedom vehicle model is built in the Simulink environment, and the ride comfort
simulation analysis of the model's vibration differential equation is carried out. First, in order to obtain the input
parameters of the model, the road roughness excitation is obtained through the simulation of the road
roughness model. The second is to obtain the front air suspension stiffness and rear air
suspension and front suspension damping coefficient and rear suspension damping coefficient
through the airbag model, the single mass system model and the single wishbone beam independent
suspension model. The road roughness excitation , the suspension stiffness , the suspension
damping coefficient and the vehicle parameters are input to the five-degree-freedom vehicle model to
obtain the seat acceleration .
(1)
where: is the road surface roughness coefficient, which is the value of the road surface power
spectral density at the reference spatial frequency , in ; n represents the spatial frequency of
several wavelengths per m length ; is the reference spatial frequency.
According to the road power spectrum density, the road roughness is divided into 8 levels, as shown in
the following table, and the geometry mean of the corresponding root mean square value ( ) of
the road roughness in the range of 0.011 <n< 2.83 is listed.
Under normal circumstances, the numerical range of road roughness is: wavelength ,
amplitude road roughness classification as shown in the table 2.
4
ICMMME2019 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 538 (2019) 012068 doi:10.1088/1757-899X/538/1/012068
The method proposed in this paper takes B-level road as an example. For reference to the literature [4],
using the pseudo white noise method, the mathematical model of a single wheel subjected to road
roughness excitation in the time domain can be described by :
(2)
In equation , is the road roughness excitation; is the road surface roughness coefficient;
is the random simulated white noise; is the cut-off frequency; is the vehicle speed.
Because the probability of a passenger car driving on a highway or a road with good road conditions is
high, here is a B-level road as an example:
According to the road roughness model and the initially set parameters, the road roughness excitation
can be obtained in the Simulink environment.
The relationship between the load F on the air spring, the absolute pressure P of the air inside the
airbag, the atmospheric pressure , and the effective area is:
(3)
5
ICMMME2019 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 538 (2019) 012068 doi:10.1088/1757-899X/538/1/012068
When the load causes the height of the airbag to change, the volume and internal pressure of the airbag
also change, and the variation law satisfies the gas state equation:
(4)
where: V represents the volume within the airbag at any position; represents the absolute pressure
of the volume within the airbag at the static equilibrium position; represents the volume within the
airbag at the static equilibrium position; represents the polytropic index. When the vehicle speed is
slow, ; when the road condition is poor, the gas state changes close to the adiabatic process,
then ; in general, take , so in this paper . The vehicle simulation speed is
.
The vertical stiffness k of the air spring is:
(5)
where:
,
Consider when the airbag is in static equilibrium position. Spring stiffness at static equilibrium
position can be obtained:
(6)
From equation , the static equilibrium spring stiffness of the front suspension airbag is:
(7)
The static equilibrium spring stiffness of the rear suspension airbag is:
(8)
The static equilibrium spring stiffness of the front suspension airbag and the static equilibrium
spring stiffness of the rear suspension airbag can be obtained by the airbag model. With the single
mass system model, the front suspension damping coefficient and the rear suspension damping
coefficient can be calculated from the values of . The front suspension stiffness and the
rear suspension stiffness can be obtained by the single wishbone beam independent suspension
model.
6
ICMMME2019 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 538 (2019) 012068 doi:10.1088/1757-899X/538/1/012068
So:
ζ (11)
According to the equation , suspension damping coefficient can be obtained.
The front suspension damping coefficient is:
ζ (12)
The rear suspension damping coefficient is:
ζ (13)
where: is the front suspension mass. is the rear suspension mass. is the stiffness of front
airbag. is the stiffness of rear airbag.
7
ICMMME2019 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 538 (2019) 012068 doi:10.1088/1757-899X/538/1/012068
(16)
(17)
(18)
(19)
(20)
Road roughness excitation , the stiffness of vehicle front suspension and rear suspension ,
the damping coefficient of front and rear suspension and vehicle parameters are available.
These parameters are input together to the five-degree-freedom vehicle model and the seat acceleration
signal is obtained.
8
ICMMME2019 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 538 (2019) 012068 doi:10.1088/1757-899X/538/1/012068
For seat:
(21)
For front wheel:
(22)
For rear wheel:
(23)
For the vehicle body:
(24)
For the torque balance of vehicle body:
(25)
Using the transfer matrix method to solve the vibration equation of a five-degree-freedom vehicle
model, can be written as:
(26)
where:
9
ICMMME2019 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 538 (2019) 012068 doi:10.1088/1757-899X/538/1/012068
4. Simulation model
In order to determine the optimal initial air pressure of the suspension airbag under the B-level road
surface, the paper first gives an arbitrary initial air pressure, and calculates the suspension stiffness and
suspension damping coefficient through the initial air pressure in the MATLAB environment. Then,
the vehicle parameters and the road roughness excitation are introduced to five-degree-freedom
vehicle model to obtain the acceleration signal of seat whose RMS value can be an evaluation
index of the ride comfort of the vehicle. The simulation of this paper consists of two parts. The first
part is the road surface roughness excitation measured by the road roughness Simulink simulation
model. The second part is the five-degree-freedom Simulink model to measure the acceleration signal
of the seat and find its RMS value. The vehicle parameters are derived from the bus of the literature
[7]. The airbag parameters are the airbag parameters that are reasonably matched according to the
vehicle parameters of the bus. The most widely used vibration comfort standards in the world are BS
6841-1987 and ISO 2631-1-1997 [9], both of which use the weighted acceleration of seat RMS as the
10
ICMMME2019 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 538 (2019) 012068 doi:10.1088/1757-899X/538/1/012068
basic evaluation index. [7] At present, when studying the ride comfort of a vehicle, the human comfort
index is generally based on the weighted acceleration RMS , and the calculation formula is as
follows:
(27)
5. Simulation results
The parameters introduced in this paper are shown in the table 3.
11
ICMMME2019 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 538 (2019) 012068 doi:10.1088/1757-899X/538/1/012068
12
ICMMME2019 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 538 (2019) 012068 doi:10.1088/1757-899X/538/1/012068
13
ICMMME2019 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 538 (2019) 012068 doi:10.1088/1757-899X/538/1/012068
14
ICMMME2019 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 538 (2019) 012068 doi:10.1088/1757-899X/538/1/012068
Figure 10 to Figure 18 show the relationship between the acceleration of the seat with time
when takes different values. It can be seen that the difference between the maximum and
minimum values of the acceleration of the person and the seat increases as the initial air pressure of
the airbag increases.
/ Pa / (m/s2)
12000 0.08882
13000 0.06459
14000 0.03773
15000 0.009946
16000 0.01768
17000 0.04326
18000 0.06132
19000 0.0707
20000 0.09029
15
ICMMME2019 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 538 (2019) 012068 doi:10.1088/1757-899X/538/1/012068
Table 4 reflects the relationship between the seat acceleration RMS value and the initial air
pressure of airbag. Figure 19 is a line diagram drawn based on the data of Table 4. It can be seen
from the line diagram that as the initial air pressure of the airbag increases, the acceleration RMS
value of the seat decreases first and then increases while the ride comfort of the vehicle increases
first and then decreases. This reflects the existence of an accurate initial air pressure of airbag with
optimum vehicle ride comfort and the initial air pressure ranges between 14kPa and 16kPa.
6. Conclusion
Taking the B-level road as an example, in order to optimize the ride comfort of the vehicle, this paper
proposes a Method for determining accurate initial air pressure range of air suspension airbag based on
vehicle ride comfort simulation analysis. Through simulation and calculation, the main conclusions are
as follows:
The five-degree-freedom vehicle model established in this paper can fully reflect the vibration
characteristics of the vehicle. Using random road excitation, air suspension parameters and other
relevant vehicle parameters as inputs, the five-degree-freedom vehicle Simulink model can simulate
the road roughness excitation (Figure 9) and the vibration displacement signal, speed signal and
acceleration signal of each part of the vehicle. As is shown in Figure 8, the Simulink model can still
output displacement, velocity and acceleration signals of other parts such as body, wheel and centroid.
This paper only selects the person and the seat acceleration signal to analyse and research.
Displacement signal, velocity signal and acceleration signal of other parts can also be used as research
objects in other fields.
There is an optimal initial air pressure of the airbag to maximize the vehicle ride comfort. It can be
seen from Figure 19 that under the B-level road surface, as the initial air pressure of the airbag
increases, the ride comfort of the vehicle increases first and then decreases. Obviously, there is an
accurate initial air pressure maximizing the vehicle ride comfort, which is consistent with the expected
result of this paper.
The method proposed in this paper can obtain a more accurate initial air pressure range of the airbag.
The optimal initial air pressure range of the vehicle airbag is between 14~16kPa (Figure 19). The
rough range given in this paper is 12~20kPa. The length of the range is reduced from 8kPa to 2kPa,
and the length of the accurate range is reduced to 1/4.
Here are a few innovative points:
(1) Using the airbag model, the single mass system model and the single wishbone beam independent
suspension model can determine the specific values of the suspension stiffness and suspension
damping coefficient corresponding to any initial air pressure of airbag.
(2) The five-degree-freedom vehicle model established in this paper can comprehensively reflect the
vibration characteristics of every part of the vehicle under the condition of rough road surface. This
model is the basis for studying the ride comfort of vehicles and analysing the vibration state of
vehicles.
16
ICMMME2019 IOP Publishing
IOP Conf. Series: Materials Science and Engineering 538 (2019) 012068 doi:10.1088/1757-899X/538/1/012068
(3) When analysing the ride comfort of a rough road surface, it is not appropriate to use the vehicle
vibration acceleration to evaluate the ride comfort of the human body. Therefore, the vibration
condition of the seat is particularly important and can be directly used to measure the ride comfort of
the vehicle.
(4) The initial air pressure of the air suspension airbag has a great influence on the ride comfort of the
vehicle. Under the B-level road surface, as the initial air pressure of the airbag increases, the ride
comfort of the vehicle increases first and then decreases, so there is an initial air pressure. Therefore
there is an accurate initial air pressure makes the vehicle ride smoother. In this paper, the method of
finding the accurate range of the initial air pressure of the airbag can be used as the initial air pressure
reference value for air suspension production and can also be used as an important reference for airbag
maintenance.
7. Prospective
The five-degree-freedom vehicle model constructed in this paper can reflect the vibration
characteristics of every part of the vehicle more comprehensively, so it can be used to study the
influence of other vehicle parameter changes on vehicle ride comfort.
This paper provides a method to determine the optimal initial air pressure range of air suspension
airbags. It can be used as a reference for the study of shock absorption. It can also be used as a method
to determine the precise initial air pressure range of airbags during airbag maintenance. It can also be
used as the basis for producing variable stiffness air suspension.
The rough range given in this paper is 12~20kPa, and the accurate range found in this paper is
14~16kPa. In theory, if the selected pressure gradient is smaller, a more accurate initial pressure range
or even an accurate initial pressure value can be obtained. However, if the selected pressure gradient is
too small, the workload simulated and calculated using the method proposed in this paper will be too
large. So the models and method of this paper need to be optimized.
References
[1] S Xiong, W Bao, Calculation and Selection of Air Suspension System Parameters. Mechanical
Manufacturing, Page 8-10, 07(2008).
[2] X Li, Research and Simulation Analysis of Dynamic Characteristics of Single Air Chamber Air Spring,
Northeastern University, Page 89(2013).
[3] J Xie, Q Zhang, Simulation Analysis of Vehicle Ride Comfort Based on Simulink,
Electromechanical Technology, Page 14-18, 01(2013).
[4] H Liu, W Zhou, J Chen, Vehicle Ride Simulation and Optimization Based on Five-degree-
freedom Vehicle Model, 33rd China Control Conference, Page 5(2014)
[5] H Zheng, Standard Road Surface Reconstruction Based on Harmonic Superposition Method,
Mechanical Research and Application, 06(2014).
[6] B Zhang, Vehicle-Road System Coupling Vibration Analysis and Ride Comfort Evaluation,
Central South University, Page 156(2013).
[7] H Zhang, W Yang, Evaluation Method of Road Flatness Based on Human-vehicle-road Five-
degree-freedom Vibration Model, Journal of Transportation Engineering, Page 16-22,
04(2010).
[8] C Liu, Study on The Transfer Function of Eight-degree-freedom Vehicle Model, Journal of
Gansu Agricultural University, Page 113-117, 03(2006).
[9] Z Yu, Automotive Theory, Mechanical Industry Press(2000).
17