Evaluating the effect of the working conditions on the ride comfort and road friendliness of the heavy truck
Van Liem Nguyen^{1} , Khac Tuan Nguyen^{2}
^{1, 2}Hubei Polytechnic University, Huangshi City, China
^{1, 2}Thai Nguyen University of Technology, Thai Nguyen City, Vietnam
^{2}Corresponding author
Vibroengineering PROCEDIA, Vol. 21, 2018, p. 8388.
https://doi.org/10.21595/vp.2018.20309
Received 12 October 2018; accepted 31 October 2018; published 13 December 2018
This study proposes a nonlinear dynamic model of the heavy trucks to evaluate the influence of the various working conditions on the ride comfort and road friendliness of the heavy vehicles based on Matlab/Simulink software. The weighted root mean square (RMS) acceleration responses of the vertical driver’s seat, the cab’s pitch and roll angles, and the dynamic load coefficient (DLC) are chosen as objective functions. The influence of the vehicle velocity and loads on the different road surfaces are evaluated based on the standard ISO 26311. The results show that the vehicle’s ride comfort is clearly affected on the road surfaces of ISO level C and D at the range of the vehicle velocity. Especially, the vehicle’s ride comfort is very poor at the velocity of 1030 m.s^{1} on ISO level D. The road friendliness of the vehicle is greatly decreased in the condition of the halfloaded vehicle. Particularly the vertical dynamic wheel loads at 2nd axle strongly affect the road damage at the velocity of 1030 m.s^{1}.
 Heavy trucks, ride comfort, road friendliness, DLC, weighted RMS acceleration responses.
 Vibration isolation system, heavy truck's vibration, vehicle dynamic model.
Keywords: heavy truck vehicle, dynamic model, ride comfort, road friendliness.
1. Introduction
The market of the heavy trucks was competing very fierce and the features of the heavy trucks were claimed higher. One of the most important requirements is to improve the ride comfort and road friendless. In order to solve these problems, the heavy truck dynamic parameters were analyzed and optimized. The influence of the heavy truck dynamic parameters on the ride comfort was analyzed via 3D vehicle dynamic model, and the suspension system parameters of the cab and the vehicle of the heavy trucks were then designed and optimized through a halfvehicle model [13], the objective functions were to improve the driver’s seat comfort and controlling the cab’s pitch angle.
In recent years, the traditional suspension systems almost were studied and replaced by using air suspension system, active or semiactive suspension systems with its capabilities of consuming less power, low cost and providing better ride quality. A halfvehicle model with hybrid balance control of a heavy truck suspension system was studied [4], and the PIDfuzzy control and wheelbase preview were combined to control the semiactive air suspension of the halfheavy truck model [5]. The heavy truck’s ride comfort also was the goal of researches.
The performance of the heavy trucks is also evaluated via the road friendliness apart from the ride comfort. The effect of dynamic load coefficients (DLC) of type suspension systems from truck axles on the road friendliness and the design of the suspension systems to minimize road damage were studied [6, 7]. An optimal design of “roadfriendly” the vehicle suspension system was also done based on the rough pavement surfaces [8], a 3D heavy truck model with air suspension systems was optimized by semiactive fuzzy control [9]. However, all the above researches almost focus to study the vehicle dynamic parameters and optimize the structure of type suspension systems to improve the ride comfort as well as the road friendliness. Influence of the various working conditions of the heavy trucks received only a little attention.
This study proposes a 3D dynamic model of a typical heavy truck, the weighted RMS acceleration responses of the driver’s seat, cab’s pitch and roll angles, and the DLC’s value are chosen as the objective function. Matlab/Simulink software is used to simulate and calculate the objections. This major goal is to evaluate the influence of the ride comfort and road friendliness of heavy truck under the various working conditions.
2. Materials and methods
2.1. Heavy truck dynamic model
A threeaxle heavy truck with the dependent suspension for the steering axle (1st axle) and the trailing arm type suspensions for two rear vehicle axles (2nd axle and 3rd axle) is selected for the vehicle dynamic analysis. A 3D model of vehicle with 13DOF is established to analyze the effect of the various operating conditions, as modeled in Fig. 1, where ${z}_{j}$ and ${m}_{j}$ are the vertical displacements and masses of the driver’s seat, the cab, the vehicle body and the axles; ${\varphi}_{\mathrm{2,3}}$ and ${\theta}_{u}$ are the angular displacements of the cab, the vehicle body and the axles; ${c}_{s}$, ${c}_{ci}$, ${c}_{j}$, ${c}_{wj}$ and ${k}_{s}$, ${k}_{ci}$, ${k}_{j}$, ${k}_{wj}$ are damping and stiffness coefficients of the suspension systems of the seat, the cab, the vehicle and the wheels; ${q}_{wj}$ is the excitation of the road on the tire; ${l}_{v}$ and ${b}_{i}$ are the distances of the vehicle ($i=$ 14, $j=$ 16, $u=$ 26, $v=$ 17).
Fig. 1. 3D dynamic model of a heavy truck
a) Side view
b) Front view
The dynamic reaction forces ${F}_{s}$, ${F}_{ci}$, ${F}_{j}$, and ${F}_{wj}$ of the suspension systems of the seat, cab, vehicle, and wheels are respectively written by:
$+{c}_{ci}[{\dot{z}}_{2}{\dot{z}}_{3}+{(1)}^{n}{l}_{n+1}{\dot{\varphi}}_{2}+{(1)}^{i+1}{b}_{1}{\dot{\theta}}_{2}+{l}_{p}{\dot{\varphi}}_{3}+{(1)}^{i}{b}_{1}{\dot{\theta}}_{3}],$
$+{c}_{j}[{\dot{z}}_{3}{\dot{z}}_{n+3}+{(1)}^{m}{l}_{q}{\dot{\varphi}}_{3}+{(1)}^{j+1}{b}_{4}{\dot{\theta}}_{3}+{(1)}^{j}{b}_{4}{\dot{\theta}}_{n+3}],$
when $i=j=$ 12, then $m=$ 1, ${l}_{p}={l}_{2}+{l}_{3}+{l}_{4}$, ${l}_{q}={l}_{5}$; when $i=j=$ 34 then $n=m=$ 2, ${l}_{p}={l}_{4}$, ${l}_{q}={l}_{6}$; when $j=$ 56 then $m=$ 2, ${l}_{q}={l}_{6}+{l}_{7}$.
Based on the vehicle dynamic model in Fig. 1, and by applying Newton’s second law of motion, the motion equations of the vehicle can be represented in the matrix form as follows:
where [$\mathbf{M}$], [$\mathbf{C}$] and [$\mathbf{K}$] are ($a$×$a$) mass, damping and stiffness matrices, respectively, {$\mathbf{Z}$} is the ($a$×$1$) displacement vector, {$\mathbf{F}\left(t\right)$} is the ($a$×1) exciting force vector, and $a$ is the number of DOF, ($a=$ 13).
2.2. Random excitation of the road surface
The vehicle’s ride comfort is strongly influenced not only by the vehicle dynamic coefficients but also by the dynamic interaction between vehicles and roads. The random excitation of the road surface can be represented with a periodic modulated random process. Thus, the spectral density of the road terrain surface is written in accordance with the ISO proposal [10] over the different spatial frequency ranges by$S\left(\mathrm{\Omega}\right)=S\left({\mathrm{\Omega}}_{0}\right).(\mathrm{\Omega}/{\mathrm{\Omega}}_{0}{)}^{{w}_{0}}$. Herein, space frequency $\mathrm{\Omega}$ is the reciprocal of the wavelength and ${\mathrm{\Omega}}_{0}=$ 0.1 m^{1} is reference space frequency, ${w}_{0}=$ 2 is the frequency index, and the value $S\left({\mathrm{\Omega}}_{0}\right)$ provides a measure for the random terrain with the reference spatial frequency ${\mathrm{\Omega}}_{0}=$ 1/2$\pi $ cycle.m^{1}.
More specifically, road surface roughness is assumed to be a zeromean stationary Gaussian random process. It can be generated through an inverse Fourier transformation:
where $N$ is the number of intervals and ${\phi}_{i}$ is a random phase uniformly distributed from 0 to 2$\pi $.
In this study, four types of road surface roughness including ISO level AD are adopted according to the standard ISO 8068 [11], and the simulation results are shown in Fig. 2.
Fig. 2. Typical road surface roughness according to ISO 8068
a)
b)
3. Results and analysis
3.1. Evaluate criteria
In order to evaluate the effect of the various operating conditions on the ride comfort and road friendliness of heavy truck, two evaluate methods are carried out including:
Basic evaluation method: According to the international standard ISO 26311 [12], the effect of vibration on ride comfort is evaluated based on the weighted rootmeansquare (RMS) acceleration responses, which is defined as in Eq. (7). Herein, ${a}_{w}\left(t\right)$ is the acceleration (translational and rotational) that depends on the time of measurement $T$.
The weighted RMS value (${a}_{wz}$) of the vertical driver’s seat, the cab’s pitch and roll angles are calculated in Eq. (7), and the results are compared with ${a}_{wz}$ given in ISO 26311:
Dynamic load coefficient (DLC): The DLC is frequently used to characterize dynamic loading of axles and it is defined as in Eq. (8) [9]. Herein, ${F}_{w,RMS}$_{}is the RMS’s value of the vertical dynamic wheel load and ${F}_{st}$ is the average vertical wheel load.
The DLC’s value is a range of 0.05 to 0.3 under normal operating conditions. The DLC’s value can reduce to the zero when the wheels move on a special smooth road or increase up to 0.4 when the wheelsroad contact is broken [7]. In this study, the DLC’s value to analyze the vehicleroad interaction and evaluate the road friendliness of the vehicle.
3.2. Effect of the vehicle velocity and the random road surfaces
In order to evaluate the effect, the vehicle velocity and road surfaces on the ride comfort, simulations are carried out when the vehicle moves on various road surfaces in the vehicle velocities. The parameters of the heavy truck are given in Table 1. The simulation results of the acceleration responses of the vertical driver’s seat, cab’s pitch and roll angles when the vehicle moves on the road surface conditions of ISO level A and C at ${v}_{0}=$ 20 m.s^{1} forward speed are shown in Fig. 3. Besides, the weighted RMS values also are listed in Table 2.
Table 1. The heavy truck dynamic parameters
Parameters

Values

Parameters

Values

Parameters

Values

Parameters

Values

${m}_{1}$ / kg

120

${l}_{1}$; ${l}_{2}$ / m

0.2; 1.1

${K}_{s}$ / kN.m^{1}

20

${C}_{s}$ / kNm.s^{1}

0.2

${m}_{2}$ / kg

500

${l}_{3}$; ${l}_{4}$ / m

1; 4.68

${K}_{ci}$ / kN.m^{1}

100

${C}_{ci}$ / kNm.s^{1}

0.75

${m}_{3}$ / kg

19 000

${l}_{5}$; ${l}_{6}$ / m

5.18; 0.62

${K}_{\mathrm{1,2}}$ / kN.m^{1}

102

${C}_{\mathrm{1,2}}$ / kNm.s^{1}

7.029

${m}_{4}$ / kg

450

${l}_{7}$ / m

1.35

${K}_{36}$ / kN.m^{1}

545.4

${C}_{36}$ / kNm.s^{1}

24.09

${m}_{5}$ / kg

1 025

${b}_{1}$; ${b}_{2}$ / m

0.38; 0.8

${K}_{w\mathrm{1,2}}$ / kN.m^{1}

690

${C}_{w\mathrm{1,2}}$ / kNm.s^{1}

1.5

${m}_{6}$ / kg

1 025

${b}_{3}$; ${b}_{4}$ / m

1; 0.4

${K}_{w36}$ / kN.m^{1}

1 380

${C}_{w36}$ / kNm.s^{1}

3.0

Fig. 3. The acceleration responses on the ISO level A and ISO level C
a) Vertical driver’s seat
b) Cab’s pitch angle
c) Cab’s roll angle
Fig. 4. The weighted RMS acceleration responses on various roads in the vehicle velocity region
a) Vertical driver’s seat
b) Cab’s pitch angle
c) Cab’s roll angle
From results in Table 2 and comparison with the standard ISO 26311 [12], the driver feels “Not uncomfortable” and “A little uncomfortable” when the vehicle moves on road surfaces of ISO level A and B. The driver feels “Fairly uncomfortable” and “Very uncomfortable” when the vehicle moves on road surfaces of ISO level C and D. Therefore, the vehicle’ ride comfort is strongly affected on the very poor road surfaces of ISO level C and D, especially of ISO level D.
The simulation results of the weighted RMS acceleration responses on the various road surfaces in a range of the vehicle velocity from 2.530 m.s^{1} are shown in Fig. 4. All weighted RMS acceleration responses in Fig. 4(a), (b) and (c) are greatly affected by the vehicle velocity. When the vehicle velocity is increased, the weighted RMS values are also increased, thus the vehicle’s ride comfort is reduced. Especially the ride comfort is strongly decreased at velocities from 1530 m.s^{1}. However, the ride comfort is insignificantly decreased at velocities from 2.57.5 m.s^{1}. In addition, all the weighted RMS values are also clearly affected by the random road surfaces in velocity region. The weighted RMS values are the smallest when the vehicle moves on the road surface of ISO level A and the weighted RMS values are the biggest on the road surface of ISO level D. Therefore, the vehicle’s ride comfort is very poor on the road surface of ISO level D, meanwhile, good on the road surface of ISO level A.
Fig. 5. The vertical dynamic wheel load at the axles
Table 2. The weighted RMS values
Parameters

Road surface roughness


ISO A

ISO B

ISO C

ISO D


${a}_{wz}$ /m.s^{2}

0.26

0.55

0.99

1.76

${a}_{w\varphi c}$ / rad.s^{2}

0.04

0.10

0.18

0.32

${a}_{w\theta c}$ / rad.s^{2}

0.09

0.20

0.36

0.64

3.3. Effect of the loaded conditions
In order to evaluate the effect of the loaded condition on the ride comfort, three conditions of the halfloaded vehicle (Half), fullyloaded vehicle (Fully) and overloaded vehicle (Over) are respectively simulated on a road surface of ISO level B in the vehicle velocity. The result is plotted in Fig. 6. The weighted RMS values of the vertical driver’s seat, cab’s pitch and roll angles in Fig. 6(a), (b) and (c) shown that the vehicle’s ride comfort is significantly affected by the loaded conditions. All the weighted RMS values are reduced in the halfloaded condition and increased in the overloaded condition in comparison with the fullyloaded condition. Especially, the results are clearly seen at the velocity of 1030 m.s^{1}.
Fig. 6. The weighted RMS values under the loaded conditions in the vehicle velocity region
a) Vertical driver’s seat
b) Cab’s pitch angle
c) Cab’s roll angle
Besides, the road friendliness is concerned under the loaded conditions at ${v}_{0}=$ 20 m.s^{1} forward speed. The dynamic reaction forces between wheels and road at the axles are shown in Fig. 5. The results shown that the vertical dynamic wheel load at the 2nd axle is the biggest meanwhile at the 1st axle is the smallest. Fig. 7 also shown that the DLC’s values at the 2nd axle is higher than both 1st axle and 3rd axle when the vehicle moves on various road surfaces of ISO level AD. This implies that the vertical dynamic wheel loads at 2nd axle greatly affect the road damage as well as the durability of the suspension system. However, this DLC’s values are in the range of 0.05 to 0.3 under normal operation conditions.
The DLC’s values are clearly affected by the velocity and the loaded conditions. At the 2nd axle and on the road surface of ISO level B, the simulation results in Fig. 8 shown that when the velocity of vehicle increasing, the DLC’s values are also increased. Besides, the DLC’s values are strongly increased in the halfloaded condition and reduced in the overloaded condition in comparison with the fullyloaded condition. Especially at the velocity of 1030 m.s^{1} and in the condition of the halfloaded vehicle, the DLC’s value is over 0.3. It implies that the wheels on the 2nd axle and road contact are broken.
Fig. 7. The DLC at the axles on the various road roughness
Fig. 8. The DLC at 2nd axle on ISO level B
4. Conclusions
The vehicle’s ride comfort is poor when the vehicle moves on road surfaces of ISO level C and D at a range of the vehicle velocity. Especially, the vehicle’s ride comfort is very poor at the velocity of 1030 m.s^{1} on ISO level D.
In the condition of the overloaded vehicle, the DLC’s values are reduced but the weighted RMS values are increased, on the contrary, in the halfloaded condition, the DLC’s values are increased and the weighted RMS values are decreased. Therefore, the vehicle’s ride comfort is reduced in the overloaded condition, and the road friendliness of the vehicle on the road surface is greatly decreased in the condition of the halfloaded vehicle.
The vertical dynamic wheel load and the DLC’s value at the 2nd axle is the biggest. Therefore, the wheels at 2nd axle greatly affect the road damage, especially at the velocity of 1030 m.s^{1 }in the halfloaded condition.
References
 Li Peng Fei, Ma Li, He Tian Ming A simulation study on vibration isolation of cab mounts in a commercial vehicle. Automotive Engineering, Vol. 27, Issue 6, 2005, p. 741743. [Search CrossRef]
 Zhang Jun Feng, He Yan Song, Yang Hai Wei Modification of cab suspension system based on pitch angular acceleration. China Mechanical Engineering, Vol. 23, Issue 18, 2012, p. 22582262. [Search CrossRef]
 Mahmoodabadi M. J., et al. A novel combination of particle swarm optimization and GA for pareto optimal design of a vehicle vibration model. Applied Soft Computing, Vol. 13, 2013, p. 25772591. [Publisher]
 Tsampardoukas Georgios, Stammers Charles W. Hybrid balance control of a magnetorheo logical truck suspension. Journal of Sound and Vibration, Vol. 317, 2008, p. 514536. [Publisher]
 Xie Zheng Chao, et al. Noiseinsensitive semiactive air suspension for heavyduty vehicles with an integrated fuzzywheelbase preview control. Mathematical Problems in Engineering, Vol. 2013, 2013, p. 121953. [Publisher]
 Buhari R., Rohani M. M., Abdullah M. E. Dynamic load coefficient of type forces from truck axles. Applied Mechanics and Materials, Vol. 405, Issue 408, 2013, p. 19001911. [Publisher]
 Cole D. J., Cebon D. Truck suspension design to minimise road damage. Proceedings of the Institute of Mechanical Engineers, Part D: Journal of Automobile Engineering, Vol. 210, 1996, p. 95107. [Publisher]
 Lu Sun Optimum design of “roadfriendly” vehicle suspension systems subjected to rough pavement surfaces. Applied Mathematical Modelling, Vol. 26, Issue 5, 2002, p. 635652. [Publisher]
 Nguyen V., Zhang J., et al. Performance analysis of air suspension system of heavy truck with semiactive fuzzy control. Journal of Southeast University, Vol. 33, Issue 2, 2017, p. 159165. [Search CrossRef]
 ISO/TC108/SC2/WG4 N57. Reporting Vehicle Road Surface Irregularities, 1982. [Search CrossRef]
 ISO 8068. Mechanical VibrationRoad Surface ProfilesReporting of Measured Data, 1995. [Search CrossRef]
 12ISO 26311. Mechanical Vibration and ShockEvaluation of Human Exposure to Whole Body VibrationPart 1: General Requirements, 1997. [Search CrossRef]
Cited By
Vibroengineering PROCEDIA
Huaxiang Zhou, Xiaoyan Wu, Peiling Wang

2020

Journal of Mechanical Engineering, Automation and Control Systems
Yongzhu Hu, Vanliem Nguyen, Yong Ye

2020

Vibroengineering PROCEDIA
Beiping Zhang, Vanliem Nguyen, Yong Wang

2020

Vibroengineering PROCEDIA
Van Liem Nguyen, Van Quynh Le

2019
