Study on the influence of bus suspension parameters on ride comfort

This study proposes a three-dimensional vibration model of bus with 10 DOF (degree of freedom) based on Dragan Sekulić model to analyze the suspension parameters directly influenced ride comfort. The suspension parameters which include the stiffness and damping parameters are analyzed based on the weighted r.m.s. (root-mean-square) acceleration responses of the space of a driver, passenger in the middle part of the bus and passenger in the rear overhang according to ISO 2631-1:1997. The results show that both stiffness and damping parameters of vehicle suspension have important influences on ride comfort. Especially, the stiffness and damping values of vehicle suspension are within the value ranges (0.5k ≤ k ≤ 0.75 k ) and (0.5c ≤ c ≤ 0.75c ) to improve the ride comfort of driver and passengers.


Study on the influence of bus suspension parameters on ride comfort
Le Xuan Long 1 , Le Van Quynh 2 , Bui Van Cuong 3   Faculty of Automotive and Power Machinery Engineering, Thai Nguyen University of Technology, Thai Nguyen, Vietnam 2 Corresponding author E-mail: 1 xuanlong_0307@yahoo.co.uk, 2 lequynhdl@yahoo.com, 3 buivancuong1301@gmail.com

Introduction
Suspension system is one of the important components to connect a vehicle body and its wheel, which not only affects the vehicle ride comfort, but also has the negative effects on the road surface as well as safe driving.Parameters matching, and optimization of suspension system are an important approach to improve vehicle ride comfort.In order to improve the ride comfort, road holding, workspace and preventing rollover, the method for the optimized design of the vehicle suspension system, considering a full vehicle model with 11-DOF based on the vertical acceleration values of seats was presented [1] and the optimal parameters of passive suspension system was found to improve ride comfort and rollover resistance including safety of the vehicle design.The MatLab/Simulink model of 7-DOF full vehicle model was developed and a genetic algorithm optimization technique was used to find the optimum values of spring stiffness and damping coefficient for front and rear passive suspension system at variable velocities which improve the vehicle ride performance [2].The optimal parameters of passive and semi-active suspension systems using the Pareto front optimization method in respect with ride comfort and road holding were presented [3]and the active control strategy was applied to the full-bus suspension system to improve ride comfort [4].The bus spatial vibration model with 10-DOF was developed to analyze the effects of vibrations on vehicle ride comfort and the analytical results have shown that the allowable vibration exposure time for the driver's body decreases as the spring stiffness of the driver's seat suspension system increases [5].
The evaluation method based on passenger crowdsourcing mode was proposed, which can fully perceive the vibration signals in different positions in the bus and indicate the vibration comfort in real time [6].The domestic vehicle was chosen as the research object which 2-DOF vibration model was established according to the vehicle dynamics theory [7].Then, the root mean square value of each vibration response was calculated by MatLab software when vehicle moves on grade C road surface and the study results indicate that the influence rule of each parameter of suspension on the vehicle ride comfort was obtained.
A full-vehicle vibration model with 10 DOF for 2-axle bus considering air suspension systems and seat suspension systems was established based on Dragan Sekulić model [5].Vibration differential equations and MatLab, Simulink model is established to simulate for full-vehicle vibration model with the random road surface profile according to the international standard ISO 8068 [8].The influences of vehicle suspension parameters such as stiffness and damping coefficients on seat's driver/passenger comfort is analyzed according to the international standard ISO 2631-1 (1997) [9].

Bus vibration model
A bus IK-301 with the dependent suspension systems including an air bag and two air bags for front and rear suspension systems was selected for establishing vehicle vibration model.And the vibration model with 10-DOF was established by Dragan Sekulić [5] to analyze the seat comfort of bus IK-301, as shown in Fig. 1.

Fig. 1. 3-D dynamic model of the bus IK-301
Fig. 1,  are the single air bag stiffness coefficients;  are the single shock-absorber damping coefficients;  are the single tire stiffness coefficients;  are the single tire damping coefficients; m are the elastic-suspended mass of the fully loaded bus;  and  are the front and rear axle mass;  are the mass of the driver and seat, mass of the passenger 1 and seat and mass of the passenger 2 and seat;  and  are the distances from the front and rear axles to the bus center of gravity;  ,  ,  and  are the distances, respectively; ,  and   are the angular displacements at the center of gravity of the elastic-suspended mass of the fully loaded bus and axles, respectively;  and  are the moments of inertia of the elastic-suspended mass of the fully loaded bus and axles, respectively;  are the random road surface profiles;  is the speed of vehicle ( = 1, 2, 3,  = 1, 2 and  = 1/6 and  = left, right).
The vehicle vibration equations can be formulated in different ways.One of the most popular methods is the Lagrange equation type II.For the vibration model showed in Fig. 1, the general vibration differential equation for the bus IK-301 is represented as the standard form of matrix equation: where, ,  and  indicate the mass matrix, stiffness matrix and damping matrix of vehicle and seats, respectively; the  ,  and  respectively refer to the acceleration, velocity, displacement vectors of vehicle and seats;  =   +   , is the vector of wheel-road contact forces, where  and  indicate the tire stiffness matrix and tire damping coefficient matrix of wheels, respectively,  is the random road excitation vector.

Random road surface excitation
There are two main methods for obtaining the random road surface excitation such as experimental measurement and numerical simulation using a time-domain model.In this study, the random road surface excitation is depicted by numerical simulation method [10].Road profile is typically assumed to be homogeneous and isotropic Gaussian random process and its statistical characteristics can be described by PSD.According to the International Standards Organization (ISO) 8608 [8], PSD of road roughness can be defined as Eq. ( 2): where,  is spatial frequency in m -1 ,  is reference spatial frequency with a value of 0.1 m -1 ,  ( ) is PSD value for reference spatial frequency in m 3 ,  is termed waviness, and reflects approximate frequency structure of the road profile, commonly taken as  = 2. Eq. ( 2) leads to estimation errors (overrated phenomenon) especially at low frequencies.To deal with this problem, PSD of road roughness based on rational white noise signal is proposed [10,11], and Eq. ( 2) can be modified as follows: where,  represents velocity of the vehicle in m/s,  represents variance of road roughness.Time domain representation of the road can be given as: where, () is the random road excitation; () is the white noise sequence.Simulink of MatLab can be used to solve the process described by Eq. ( 4) to generate random road profiles, the entry of the process is the standard deviation of the white noise process which depends of the road surface category, a 80 km/h-22,2 m/s speed and other parameter values [12] are set.Simulink model of random road generator is shown in Fig. 2

Vehicle ride comfort evaluation method
Currently there are many methods to evaluate the vehicle ride comfort such frequency-domain method, time-domain method.This study is based on ISO 2631-1 (1997), the vibration evaluation based on the basic evaluation method including measurements of the weighted root-mean-square (r.m.s.) acceleration is defined by: where,  () is the weighted acceleration (translational and rotational) as a function of time, m/s 2 ;  is the duration of the measurement, s.
In this way, a synthetic index-called vertical weighted rms acceleration,  can be calculated from formula Eq. ( 5) and the rms.value of the vertical acceleration in vehicle would be compared with the values in Table 1, for indications of likely reactions to various magnitudes of overall vibration in the public transport.

Simulation and analysis results
In this study, MatLab software and Simulink toolbox are selected to solve Eq. ( 1) with the vehicle simulation parameters [5] and the random road surface excitation.The simulations are carried out when vehicle moves on the ISO good road (grade B) condition and the bus speed of 22,2 m/s with full load.The vertical acceleration responses of the space of a driver, passenger in the middle part of the bus and passenger in the rear overhang is shown in Fig. 3.
The value of  is defined as 0.199 m.s -2 and this value of  satisfies the conditions for driver comfort (according to Table 1); the values of  and  are defined as 0.621 m.s -2 , 0.685 m.s -2 , respectively and this values satisfy uncomfortable conditions for passengers (according to Table 1).Both stiffness and damping parameters of vehicle suspension have important influences on vehicle ride comfort.In order to analyze the influence of bus suspension parameters such as stiffness and damping coefficients on ride comfort, discussions will present in the following section.Fig. 5 shows that vehicle suspension damping coefficient increases from 0.25 to 0.75 , the  values tends to decrease, which makes driver's ride comfort considerably improve.And then the  and  values tends to increase quickly, which has the negative effects on ride comfort of passengers.
According to Fig. 4 and Fig. 5, in order to decrease the vehicle vibrations and hence improve the ride comfort, the values of the vehicle suspension parameters are within the value ranges (0.5 ≤  ≤ 0.75 ) and (0.5 ≤  ≤ 0.75 ).

Conclusions
This study presents a methodology to analyze the effects of bus suspension parameters on ride comfort based on the full-vehicle vibration model using the random road surface profile and the analysis results as follow: 1) In the suspension system, the reasonable matching of suspension stiffness and damping can efficiently improve seat ride comfort, which keeping ride comfort keeping ride comfort within a safety limits and the safety of the road surface.
2) In order to decrease the vehicle vibrations and hence improve the ride comfort, the values of the vehicle suspension parameters are within the value ranges (0.5 ≤  ≤ 0.75 ) and (0.5 ≤  ≤ 0.75 ).

2 .
(a) and the simulation result of the typical grade B road surface is shown in Fig. 2(b).a) Simulink model of random road generator b) Grade B road profile Fig. Simulink model of random road generator and simulation result

3 . 4 . 5 .
a) The space of a driver b) Passenger in the middle part c) Passenger in the rear overhang Fig.The vertical acceleration responses of bus seats 5.1.Influence of vehicle suspension stiffness Vehicle suspension stiffness is an important factor that influences the ride comfort.To analyze its effect on the  ,  and  values, the values of suspension stiffness coefficients  = 0.25 0.5 0.75 1 1.25 1.5 1.75 2.0 ×  , where  =  ,  ,  ,  and the values of suspension damping coefficients  = 0.25 0.5 0.75 1 1.25 1.5 1.75 2.0 ×  , where  =  ,  ,  ,  were analyzed when vehicle moves on the ISO good road (grade B) condition and the bus speed of 22,2 m/s with full load, where  and  are used to designate the vehicle stiffness and damping coefficients in the reference document [5].The influence of vehicle suspension stiffness coefficients on  ,  and  values are shown in Fig. 4. From Fig. 4 shows that when the vehicle suspension stiffness coefficient increases, the  ,  and  values tend to increase, which makes the negative effects on ride comfort of driver and passengers.The  values increase quickly when the values of the damping coefficients are  = 0.25 and the  ,  and  values do not change much when the values of the damping coefficients are achieved by great values.When the values of the damping coefficients increase from 1.5 to 2 , the  ,  and  values decrease about from 3.6 to 3.76 %, which makes the ride comfort of driver and passengers negligible.a) The  values at vertical driver's seat b) The  values at passenger in the middle part c) The  values at passenger in the rear overhang Fig. Influence of vehicle suspension stiffness coefficients on ride comfort of bus seats 5.2.Influence of vehicle suspension damping Vehicle suspension damping is another important factor that influences the ride comfort.To analyze its effect on the  ,  and  values, the values of suspension damping coefficients  = 0.25 0.5 0.75 1 1.25 1.5 1.75 2.0 ×  , where  =  ,  ,  ,  and the values of suspension stiffness coefficients  = 0.25 0.5 0.75 1 1.25 1.5 1.75 2.0 ×  , where  =  ,  ,  ,  were analyzed when vehicle moves on the ISO good road (grade B) condition and the bus speed of 22,2 m/s with full load.The influence of vehicle suspension damping coefficient on  ,  and  values are shown in Fig. 5. a) The  values at vertical driver's seat b) The  values at passenger in the middle c) The  values at passenger in the rear overhang Fig. Influence of vehicle suspension damping coefficient on ride comfort of bus seats 0 Vehicle suspension damping/(N.s.m -1 0 0.5k 0 0.75k 0 1.0k 0 1.25k 0 1.5k 0 1.75k 0 2.0k 0 xc 0 STUDY ON THE INFLUENCE OF BUS SUSPENSION PARAMETERS ON RIDE COMFORT.LE XUAN LONG, LE VAN QUYNH, BUI VAN CUONG 82 VIBROENGINEERING PROCEDIA.DECEMBER 2018, VOLUME 21

Table 1 .
Comfort levels related to  threshold values