Ride comfort analysis of a nonlinear vehicle excited by the consecutive speedcontrol humps
Shan Liang^{1} , Yongsheng Sun^{2} , Qin Zhu^{3} , Zhiyong Yang^{4} , Cheng Zeng^{5}
^{1, 2, 4, 5}College of Automation, Chongqing Univers, Chongqing 400044, China
^{3}Department of Mechanical Engineering, Oyama National College of Technology, Oyama 3230806, Japan
^{4}College of Information and Engineering, Chongqing Vocational Institute of Engineering, Chongqing 400037, China
^{1}Corresponding author
Journal of Vibroengineering, Vol. 15, Issue 4, 2013, p. 16561664.
Received 18 June 2013; accepted 5 December 2013; published 31 December 2013
JVE Conferences
The consecutive speedcontrol humps (SCHs) possess the function of controlling speed forcibly, but cause violent vibration of a vehicle inevitably. This paper tries to further explore the inherent link among parameters of the SCHs, velocity and vehicle vibration. A 4DOF nonlinear halfvehicle model with nonlinear springs and nonlinear dampers is established. The consecutive SCHsspeed coupling excitation function is presented by combination of trapezoidal and sine wave of constant amplitude and variable frequency. The nonlinear dynamics of halfvehicle model is investigated by numerical simulation. It reveals that various forms of vibrations, such as periodic, quasiperiodic and chaotic vibrations, could appear in the system with the change of the velocity. Further it is found that quasiperiodic motions will affect vehicle ride comfort most and can be avoided by changing parameters of the consecutive SCHs. Results are conductive to deep understanding of nonlinear vibration in vehicle and rational design of the consecutive SCHs.
Keywords: ride comfort, consecutive speedcontrol humps, nonlinear vibration, halfvehicle model.
1. Introduction
The speedcontrol hump (SCH) which is installed on highway is currently widely used as one of the compulsive speed control facilities. In general, the consecutive SCHs are local elevations of residential road surface of limited height, which are placed in series at hundreds or thousands of meters with a fixed interval. When the automobile passes over the consecutive SCHs at a certain velocity, the impact from the uneven road surface will inevitably bring to the vehicle unwanted vibrations and noise. Prolonged exposure to repeated vibrations and impacts of wholebody nature when vehicle travels on uneven road has been associated to occupational health disorders [1]. Therefore the effect of the consecutive SCHs on vehicle vibrations, and in particular on the ride comfort, is still needed to be further studied [25].
Since the automobile is a kind of highly nonlinear system, the vibrations caused by road surface will appear in nonlinear characters, such as chaos and bifurcation. To investigate nonlinear response of vehicle model excited by road surface, some productive works have been done by scholars. In these studies researchers often rely on the quartercar model for studying heave motion as this model is the simplest to analyze and yet can reasonably predict the response of a system [68]. However the actual vehicle is a more complex nonlinear system than a quartercar model system. The vehicle speed, pitch motion, roll motion and the impact of front and rear wheels are also important to the dynamics of actual vehicle. Zhu et al. [9, 10] have investigated the possibilities of chaotic vibrations which existed in halfcar and fullcar models under sine wave excitations. However the influence of vehicle speed on vehicle dynamics and how the nonlinear vibrations affect the vehicle ride comfort have rarely been investigated in the current studies.
In this paper the nonlinear vibration behaviours of the 4DOF halfvehicle model with nonlinear springs and damping elements through the excitations from the consecutive SCHs are investigated. The study begins by establishing a combined consecutive SCHsspeed dynamic excitation function and the nonlinear 4DOF halfvehicle mechanical model. Then the numerical simulation is carried out for investigating the dynamic responses of the simplified vehicle system. Results show that the complex behaviours, such as quasiperiodic and chaotic vibrations, may exist in the system. Quasiperiodic motion deteriorates vehicle ride comfort significantly, which can be avoided by changing the parameters of the consecutive SCHs.
2. The excitation model of the consecutive SCHs
In general, there are some shapes of the consecutive SCHs in the special sections on the highway. When a vehicle passes over a series of SCHs with a certain speed, the trajectory of the vehicle can be approximately considered as the periodic trapezoidal wave [11]. The consecutive speedcontrol humps and its model of trapezoidal excitation are illustrated as in Fig. 1a and Fig. 1b, respectively.
Fig. 1. Consecutive SCHs excitation
a) Consecutive SCHs road
b) Road profile of the SCH surface
In Fig. 1b, $h$ is the amplitude of consecutive SCHs excitation, ${s}_{1}$ and ${s}_{2}$ are the width of a hump and the gap between two adjacent humps, ${T}_{\text{1}}$ is the time to pass a hump for a vehicle, $[{T}_{1},T]$ is the time to pass the gap between two adjacent humps. If we take the roughness of the road surface into account, then a sinusoidal wave can be viewed as the intrinsic excitation of the road surface. The combination of a sinusoidal wave and a trapezoidal wave is used to describe the excitation generated by consecutive speedcontrol hump area on highway. Thus, the consecutive speed control humps can be approximated as:
where $A$ is the amplitude of sinusoidal excitation, ${f}_{2}$ is the intrinsic excitation frequency caused by road surface roughness, $trap\left(t\right)$ stands for the trapezoidal wave, which can be expressed as follows:
where ${T}_{1}=T/2$, ${t}_{a}={T}_{1}/4$ and ${t}_{b}={T}_{1}{T}_{1}/4$ and $T$ is the period of consecutive SCHs excitation. Define ${f}_{1}$ is the frequency of consecutive SCHs excitation, then:
Obviously ${f}_{1}<{f}_{2}$. Assume the following equation holds:
where $\sigma $ is defined as the ratio of ${f}_{1}$ and ${f}_{2}$, ${s}_{1}$ and ${s}_{2}$ are width of a hump and the gap between two adjacent humps, respectively.
For halfvehicle model, the excitations to the front and rear tire are defined respectively as:
where $\Delta t$ is the time delay between the forcing functions ${x}_{fd}$ and ${x}_{rd}$, which can be calculated by the vehicle speed and the distance between front and rear wheels.
3. FourDOF suspension model
Assuming the vehicle has a bilaterally symmetric structure, it is reasonable to simplify the vehicle as a fourDOF suspension model as in Fig. 2 [9].
Fig. 2. Nonlinear fourDOF halfvehicle model
The notations used in the model (Fig. 2) are defined as following:
${m}_{b}$: vehicle body mass; ${m}_{f}$, ${m}_{r}$: front and rear unsprung masses, respectively; $\theta $: angular displacement of vehicle body; ${x}_{b}$: vertical displacement of vehicle body; ${x}_{f}$: heave displacement of ${m}_{f}$; ${x}_{r}$: heave displacement of ${m}_{r}$; ${x}_{fd}$, ${x}_{rd}$: excitations to the front and rear tires; ${F}_{sf1}$, ${F}_{sr1}$: front and rear nonlinear tire spring forces;$\mathrm{}{F}_{sf2}$, ${F}_{sr2}$: front and rear nonlinear suspension spring forces; ${F}_{cf1}$, ${F}_{cr1}$: front and rear nonlinear tire damper forces; ${F}_{cf2}$, ${F}_{cr2}$: front and rear nonlinear suspension damper forces. By analyzing the mechanical model, the equation of ${x}_{rd}$, the equations of motion can be expressed as:
where $J$ is vehicle body inertia. The spring forces ${F}_{sij}$ are assumed to have the following characteristics [12, 13]:
where $i=f,r$ denotes front or rear wheel, $j=1,2$ denotes suspension or tire, the definitions of $i$ and $j$ are suitable for other formulas in the paper, ${k}_{ij}$ is the equivalent stiffness, ${\Delta}_{ij}$ is the deformation of the spring, $n$ is an exponent suspension spring. Therefore ${\Delta}_{f1}$, ${\Delta}_{f2}$, ${\Delta}_{r1}$ and ${\Delta}_{r2}$ can be given as follows:
${\Delta}_{f1}={x}_{f}{\Delta}_{sf1}{x}_{fd}$, ${\Delta}_{f2}={x}_{b}{\Delta}_{sf2}{x}_{f}{l}_{f}\mathrm{s}\mathrm{i}\mathrm{n}\theta $,
${\Delta}_{r1}={x}_{r}{\Delta}_{sr1}{x}_{rd}$, ${\Delta}_{r2}={x}_{b}{\Delta}_{sr2}{x}_{r}+{l}_{r}\mathrm{s}\mathrm{i}\mathrm{n}\theta $.
The dampers of tires and springs are assumed to be viscous, so the characteristics of damping force can be expressed as:
where ${F}_{cij}$ is the damping force and $\Delta {\dot{x}}_{ij}$ is the relative velocity of the extremes of the damper. The damping coefficients are defined as:
where ${c}_{iju}$ and ${c}_{ijd}$ are the damping coefficients for tension and compression, respectively.
Letting ${x}_{1}={x}_{b}$, ${x}_{2}={\dot{x}}_{b}$, ${x}_{3}=\theta $, ${x}_{4}=\dot{\theta}$, ${x}_{5}={x}_{f}$, ${x}_{6}={\dot{x}}_{f}$, ${x}_{7}={x}_{r}$, ${x}_{8}={\dot{x}}_{r}$, the state equations of the 4DOF suspension system can be expressed as:
where ${m}_{bf}$ and ${m}_{br}$ are the partial masses for front and rear axles of the system, and:
When the suspension system is in the relatively static state of equilibrium, there are no external excitations and relative speed, then the following equations are valid:
Therefore, the static deformations ${\Delta}_{sf1}$, ${\Delta}_{sf2}$, ${\Delta}_{sr1}$ and ${\Delta}_{sr2}$ can be evaluated by (12).
Assume that ${\delta}_{xb}$ is the initial heave displacement of ${m}_{b}$ and ${\Delta}_{\theta}$, then the initial angular displacement ${\delta}_{xb}$ and ${\Delta}_{\theta}$ can be calculated as:
Therefore the initial condition of state variables can be set as:
4. Numerical analysis
Owing to the nonlinearity of the system (10), the dynamic response of the vehicle model can be studied numerically by the fourthorder RungeKutta algorithm provided by MATALB.
Table 1 shows parameters for numerical simulation as an example [9].
Table 1. Parameters for numerical simulation
Parameter

Value

Parameter

Value

Parameter

Value

${m}_{b}$

1180.0 kg

$J$

633.615 kg·m^{2}

${c}_{f2u}$, ${c}_{r2u}$

500 N·m/s

${m}_{f}$

50.0 kg

$g$

9.81 N/kg

${c}_{f2d}$, ${c}_{r2d}$

359.7 N·m/s

${m}_{r}$

45.0 kg

${l}_{f}$

1.123 m

${c}_{f1u}$, ${c}_{r1u}$

10 N·m/s

${k}_{f2}$

36952.0 N/m

${l}_{r}$

1.377 m

${c}_{f1u}$, ${c}_{r1u}$

10 N·m/s

${k}_{r2}$

30130.0 N/m

${n}_{f2}$, ${n}_{r2}$

1.5

$h$

0.02 m

${k}_{f1}$, ${k}_{r1}$

14000.0 N/m

${n}_{f1}$, ${n}_{r1}$

1.25

${s}_{1}$, ${s}_{2}$

0.5 m

$A$

0.0015 m

$\sigma $

10

4.1. Nonlinear dynamics analysis of the system
In general, for nonlinear dynamics analysis of a system, the transient process and stationary process are two issues of crucial importance. However the kinetic characteristics of a mechanical system could be evident only when the system tends to be stable. In order to investigate the nonlinear dynamics of the studied system, in this paper the system kinetic characteristics are explored by investigating dynamics of stationary process and the corresponding speed ranges of the different kinetic performances are obtained through speed bifurcation diagrams, phase portraits, Poincaré maps and power spectrums.
In this section speed bifurcation diagrams are first introduced for analyzing the dynamics of the vehicle suspension system, which are obtained by plotting the Poincaré points of the state variables (${x}_{b}$ and $\theta $ etc. in this paper) at different velocities. For a certain speed, the complexity of the system is usually proportional to the number of Poincaré points. Singlepoint indicates periodic motion and multipoint may mean perioddoubling, quasiperiodic or chaotic motion. Speed bifurcation diagrams could roughly represent the system kinetic characteristics and intuitively investigate the effects of speed on the dynamics of the system.
Fig. 3. Speed bifurcation diagram of ${x}_{b}$
Fig. 4. Speed bifurcation diagram of $\theta $
Fig. 34 respectively represent the speed bifurcation diagrams of ${x}_{b}$ and $\theta $ when $v$ varies in 45.0 km/h ~ 70.0 km/h, where simulation time: $t=$20 s, simulation step: $\Delta v=$0.2 km/h. As can be shown in the Figures, there are visible differences in the vibration performance when the speed is in different scopes. Based on these, the considered speed range can be broadly classified into three regions: region A (45.0 km/h ~ 50.6 km/h), region B (50.6 km/h ~ 62.0 km/h), region C (62.0 km/h ~ 70.0 km/h). It’s obvious that the system kinetic characteristics are similar when $v$ is in a same region. To investigate specific kinetic characteristics in the three regions, the phase portrait, Poincaré map and power spectrum at $v=$47.6 km/h, 54.0 km/h, 64.0 km/h are investigated in Fig. 57, respectively.
As shown in Fig. 57, the characteristics of vertical vibration of vehicle body are separately periodic, quasiperiodic or chaotic when speeds are 47.6 km/h, 54.0 km/h or 64.0 km/h, which implies that the system performs stable periodic motion in region A; when $v$ varies in region B, the system appears in the quasiperiod motion; as $v$ increases to region C, the chaotic motion takes place in the system.
Fig. 5. Periodic motion ($v=$47.6 km/h)
Fig. 6. Quasiperiodic motion ($v=$54.0 km/h)
Fig. 7. Chaotic motion ($v=$64.0 km/h)
4.2. Ride comfort analysis of vehicle model
By the analysis above, the nonlinear vibrations of the vehicle may present periodic, quasiperiodic and chaotic features when system is in the stable regimes. But in fact the time to pass over the consecutive SCHs is relatively short for a real car. Therefore, in order to analyse vehicle model under the excitation of the consecutive SCHs, transient process is also important. In this paper vehicle Vibration Amplitude (VA) and Vibration Strength (VS) are introduced for ride comfort analysis of vehicle model.
Fig. 8. Timedomain responses of vibration amplitude (blue: $v=$ 47.6 km/h, red: $v=\mathrm{}$54.0 km/h, green: $v=\mathrm{}$64.0 km/h)
Fig. 9. Timedomain responses of vibration strength (blue: $v=$ 47.6 km/h, red: $v=\mathrm{}$54.0 km/h, green: $v=\mathrm{}$64.0 km/h)
Fig. 8 and Fig. 9 separately represent the timedomain responses of the VA and VS of the vehicle model under the consecutive SCHs excitation. In both of them it is supposed that the vehicle begins to pass into the consecutive SCHs zone when Time$=$2 s, and the three kinds of lines with different colors are the timedomain curves of the vehicle vibrations when the vehicle is at 47.6 km/h, 54.0 km/h and 64.0 km/h, which represent nonlinear vibrations in region A, B, and C, respectively.
As shown in Fig. 8 and Fig. 9, there are significant differences in the behaviours of the system when the speed is in region A, B and C. In region B, the VA and VS of the system diverges dramatically since the vehicle begins to pass into the consecutive SCHs zone and keeps quasiperiodic motions at high levels when stable, which means ride comfort gets worse sharply when passing over the consecutive SCHs. In the meantime, the VA and VS values of the system, speeds of which belong to region A or C, are convergent at the transient process and fluctuate at relatively low levels, which means that the effects of extra vibrations caused by consecutive SCHs on ride comfort are not obvious.
4.3. Relations between the consecutive SCH and vehicle ride comfort
From the analysis above, it is clear that vehicle model has a bad ride comfort when the system meets the conditions of quasiperiodic motions. For the sake of safety, vehicles are required to travel over a speed hump road in a restricted range of speed. Therefore it’s a feasible way to design the appropriate parameters of consecutive SCH and avoid such scope of speed becoming to the quasiperiod region. Table 2 represents the relations between the parameters of the consecutive SCHs and the range of speed in which results in vehicle excessive vibrations. $s$ is defined as the width of SCH, $s={s}_{1}+{s}_{2}$, $\delta $ is defined as the duty ratio of humps, $\delta /100={s}_{1}/({s}_{1}+{s}_{2})$.
As shown in Table 2, the values of speed range become larger when the width of SCH increases. However, the values decrease with the duty ratio of humps increasing. It means one can obtain a good ride comfort by changing $s$ or $\delta $ (or the both ones) so that the vehicle which is traveling in the restricted speed will not cause excessive vibrations.
Table 2. Relations between the parameters of consecutive SCH and vehicle speed of great vibrations
$\delta =50$, $s$ / (m)

The speed range / (km/h)

${s}_{1}=$0.05 m, $\delta $

The speed range / (km/h)

0.8

38.7 $<v<$ 49.8

45

56.0 $<v<$ 66.7

1

50.6 $<v<$ 62.0

50

50.6 $<v<$ 62.0

1.2

60.3 $<v<$ 73.2

55

46.4 $<v<$ 55.3

5. Conclusions
In the paper the nonlinear dynamics of a 4DOF vehicle suspension model under the consecutive SCHs excitation is studied through numerical simulation. Results show that the periodic motion, quasiperiodic motion and chaotic motion will occur in the system when the vehicle travels over the consecutive SCHs with the increase of vehicle velocity. Vehicle ride comfort becomes bad when quasiperiodic response occurs in the system and it can be improved effectively by changing the parameters of consecutive SCHs.
Although the mechanical model of vehicle and the consecutive SCHs are only simplified ones and the parameters used do not agree precisely with the practical data for actual vehicles, the results may still be useful in the design of a ground vehicle and the cognition of dynamics response of vehicle when driving in a special pavement structure.
Acknowledgements
The authors would like to thank anonymous referees for their helpful comments and suggestions, and the Doctoral Fund of Ministry of Education of China (20100191110037) for support.
References
 BengtOlov Wikström, et al. Health effects of longterm occupational exposure to wholebody vibration: a review. International Journal of Industrial Ergonomics, Vol. 14, 1994, p. 273292. [Search CrossRef]
 A. Sezgin, Yz. Arslan. Analysis of the vertical vibration effects on ride comfort of vehicle driver. Journal of Vibroengineering, Vol. 14, 2012, p. 559571. [Search CrossRef]
 Zifan Fang, et al. Mechatronics in Japanese rail vehicles: active and semiactive suspensions. Control Engineering Practice, Vol. 10, 2011, p. 428437. [Search CrossRef]
 WeiYen Wang, et al. Hierarchical fuzzyneural control of antilock braking system and active suspension in a vehicle. Automatica, Vol. 48, 2012, p. 16981706. [Search CrossRef]
 G. Priyandoko, M. Mailah, H. Jamaluddin. Vehicle active suspension system using skyhook adaptive neuro active force control. Mechanical Systems and Signal Processing, Vol. 23, 2009, p. 855868. [Search CrossRef]
 S. Liang, C. Li, Q. Zhu, Q. Xiong. The influence of parameters of consecutive speed control humps on the chaotic vibration of a 2DOF nonlinear vehicle model. Journal of Vibroengineering, Vol. 13, 2011, p. 406413. [Search CrossRef]
 S. Li, S. Yang, W. Guo. Investigation on chaotic motion in hysteretic nonlinear suspension system with multifrequency excitations. Mechanics Research Communications, Vol. 31, 2004, p. 229236. [Search CrossRef]
 Utz Von Wagner. On nonlinear stochastic dynamics of quarter car models. Int. J. NonLinear Mech., Vol. 39, 2004, p. 753765. [Search CrossRef]
 Q. Zhu, M. Ishitobi. Chaos and bifurcations in a nonlinear vehicle model. Journal of Sound and Vibration, Vol. 275, 2004, p. 11361146. [Search CrossRef]
 Q. Zhu, M. Ishitobi. Chaotic vibration of a nonlinear fullvehicle model. International Journal of Solids and Structures, Vol. 43, 2006, p. 747759. [Search CrossRef]
 F. Liu, S. Liang, Q. Zhu, Q. Y. Xiong. Effects of the consecutive speed humps on chaotic vibration of a nonlinear vehicle model. ICIC Express Letters, Vol. 4, 2010, p. 16571664. [Search CrossRef]
 A. Moran, M. Nagai. Optimal active control of nonlinear vehicle suspension using neural networks. JSME International Journal, Vol. 37, 1994, p. 707718. [Search CrossRef]
 J. C. Dion. Tires, Suspension and Handling. 2nd Edition, Society of Automotive Engineers, Warrendale, PA, 1996. [Search CrossRef]