Vibration analysis of a beam on a moving vehicle under the road excitation with different contact models
S. Q. Wu^{1} , H. X. Ai^{2} , Y. Z. Guo^{3} , Q. G. Fei^{4}
^{1, 2, 3, 4}Department of Engineering Mechanics, Southeast University, Nanjing, P. R. China
^{1, 2, 3, 4}Jiangsu Key Laboratory of Engineering Mechanics, Southeast University, Nanjing, P. R. China
Journal of Vibroengineering, Vol. 15, Issue 4, 2013, p. 16891700.
Received 11 July 2013; accepted 18 August 2013; published 31 December 2013
JVE Conferences
Dynamic analysis of a beam on a moving vehicle is presented in this paper. The vehicle is simulated by a four degreesoffreedom massspring system and the beam on top is supported by springdamping systems. Two contact models named the ‘point contact’ and the ‘patch contact’ respectively, are adopted to simulate the interaction between road surface and vehicular tyres. The equation of motion of the beamvehicle system is formulated and the dynamic response on the beam under the excitation of the irregular road surface is derived. Numerical simulations are conducted to demonstrate the influence of different factors, such as the length of the contact, the velocity of vehicle, the road condition and the bracing stiffness, etc. on the vibration level of the beam structure, which aims to provide references on the vibration problem in transporting a beamshaped package.
Keywords: beamvehicle system, moving, road excitation, dynamic response, contact model.
1. Introduction
Vibration problems generally existed during road transportation especially when packages are flexible with large slenderness ratio, such as missiles, prefabricated concrete beams, etc. These packages can be regarded as beamshaped structures which may suffer from long distance transportation before service and damages may initiate due to the vibration caused by the unevenness of the road profile. There are mainly two kinds of road excitations: (1) the stochastic excitation due to the road surface roughness; (2) the impact excitation due to an irregularity causing bouncing of vehicle. The impact excitation will result in a sudden increase in the dynamic response, which may cause damages in structure as well as the equipment inside; the amplitude of the vibration on packages under stochastic excitation is relatively small compared with that from the impact, but it may cause fatigue damage in the packages during long distance transportation. Therefore, to study the dynamic response of the beamshaped packages under the road excitation is very important for the condition assessment and vibration control during transportation.
Research work on the vibration analysis of vehicle as well as the packages on top due to the excitation from irregular road profile has been conducted for decades. Shi and Cai [1] investigated the dynamic effects of vehicle by a threedimensional vehiclepavement coupled model and concluded that the dynamic vehicular axle loads are significantly higher than the static ones under rough road condition. Michaltsos [2] proposed a mathematical model to determine the critical velocity when the vehicle loses touch with the road surface during the bouncing due to an irregularity. Lu et al. [3] investigated the acceleration on vehicle under both the stochastic and impact excitations at different traveling speeds with experimental and analyzing approaches, results indicate that the effect of truck speed on root mean square acceleration including shock and vibration was strong at lower speed, but slight at a higher speed. Lak et al. [4] studied the relationship among the road unevenness, the dynamic vehicle response and groundborne vibration with in situ measurements of road unevenness and a 3D twoaxle vehicle model, results show that the predicted RMS value of dynamic response on vehicle is generally larger than the measurement data. Similarly, based on the measurement of road surface roughness and a halfcar model, Barbosa [5] studied the vibration response of vehicle using a spectral method in frequency domain.
All the aforementioned research works focus on the vibration of vehicle body due to road excitation. During the road transportation of packages, more attention should be paid on the vibration level of the package and the design of its bracing system. Ghaith [6] proposed a nonlinear dynamic model to analyze the dynamic response of a BernoulliEuler beam fixed on a moving cart carrying lumped tip mass subjected to external periodic force, which aims to describe the motion of structures like forklift vehicles or ladder cars that carry heavy loads and military airplane wings. Ragulskis et al. [7] studied the dynamic response of a container with package during transportation and suggested that to avoid the overload, the velocity of vehicle must ensure the maximum difference between the excitation frequency and the nature, frequency of the system.
In most of the research works on vibration analysis of package during road transportation, the contact between vehicular tyres and road surface is assumed as a point in the dynamic model. In reality, pneumatic tyres will deform due to static/dynamic loads from the vehicle body which results in a footprint area called patch contact between tyres and road surface. Studies on modeling of pneumatic tyres [8] and the relationship among the loads, deformation of tyres and the inflation pressure have been conducted [9]. Yin et al. [10] introduced the patch contact model of pneumatic tyres in bridgevehicle interaction problem, results show that treating the interaction between tyres and road surface as point contacts may overestimate the dynamic deflection of the bridge. With patch contact models adopted, the accuracy of the dynamic response of the bridge under moving vehicular loads can significantly be improved compared with that analyzed with point contact simulations.
Since the patch contact model shows great advantage over the point contact model in bridgevehicle interaction problem, a study is conducted in this paper by introducing the patch contact model in the vibration analysis of a beamshaped package on a moving vehicle under road excitation. The package and vehicle are modeled by a BernoulliEuler beam and a four DegreesOfFreedom (DOFs) massspring system respectively, the bracing systems between the package and vehicle are simulated by springdamping systems. The road surface roughness is assumed to have a Power Spectral Density (PSD) defined according to the ISO standard [11]. A comparison of the dynamic response on the beam with the two different contact models adopted will be given and different factors which may affect the vibration level of the beam will also be discussed.
2. The beamvehicle model
The schematic diagram of the beamvehicle model with the point contact model for tyres is shown in Fig. 1a, where ${F}_{1}\left(t\right)$ and ${F}_{2}\left(t\right)$ are the two interaction forces at the braces between the vehicle body and the beam respectively. ${F}_{3}\left(t\right)$ and ${F}_{4}\left(t\right)$ are the two interaction forces between tyres and road surface. Parameters ${k}_{i}$ and ${c}_{i}(i=\mathrm{1,2})$ are the stiffness and the damping of the braces between the vehicle and beam respectively. ${l}_{i}\left(i=\mathrm{1,2}\right)$ is the distance from the left end of the beam to the $i$th bracing system. ${b}_{1}{L}_{0}$ and ${b}_{2}{L}_{0}$ represent the distance between the bracing system and the center of the vehicle body respectively.
Fig. 1. The beamvehicle model
a) The schematic diagram
b) The patch contact model
Similarly ${a}_{1}S$ and ${a}_{2}S$ are the distances between the wheels and the center of the vehicle body respectively. The symbols ${m}_{i}$, ${K}_{si}$, ${C}_{si}$ and ${y}_{i}\left(i=\mathrm{1,2}\right)$ represent the mass, stiffness, damping and displacment of the $i$th suspension system respectively. Variables ${y}_{i}$ and ${\widehat{y}}_{i}(i=\mathrm{3,4})$ are the deflection at the bracing systems on the beam and the vehicle body respectively. ${y}_{V}$ and ${\theta}_{V}$ are the vertical displacement and pitch displacement of the vehicle body respectively. Location #1, #2 and #3 are the right end, 1/4 point and midpoint of the beam respectively. ${K}_{ti}$, ${C}_{ti}$ are the stiffness and damping of the $i$th tyre modeled by point contact model. The patch contact model for tyre is shown in Fig. 1b and further explainations will be given in Section 2.2.
2.1. Finite element model of the beamshaped package
The beamshaped package is modeled by a BernoulliEuler beam. The equation of motion is given as:
where $\rho A$, $c$ and $EI$ are the mass per unit length, damping and the flexural rigidity respectively. $u(x,t)$ is the deflection on the beam at location $x$ and time $t$. The term in the right side of Eq. (1) represents the excitation forces from the bracing systems. ${N}_{p}$ is the number of the braces equal to 2 in this study. Though only two bracing systems are involved in the beamvehicle system, more forces can be included in Eq. (1) without additional changes of the expression. By adopting the finite element discretization, the finite element model can be obtained as:
where ${\mathbf{M}}_{s}$, ${\mathbf{C}}_{s}$ and ${\mathbf{K}}_{s}$ are the mass, damping and stiffness matrices of the beam model respectively; $\mathbf{X}$, $\dot{\mathbf{X}}$ and $\ddot{\mathbf{X}}$ are the vectors of nodal displacements, velocities and accelerations respectively; ${\mathbf{H}}_{s}$ is the location matrix for the dynamic interaction forces ${\mathbf{F}}_{d}\left(t\right)={\left\{{F}_{1}\left(t\right){F}_{2}\left(t\right)\right\}}^{T}$ at the braces which has the following expression:
where $N$ is the number of DOFs of the beam with boundary conditions. ${\mathbf{H}}_{i}$ represents the shape functions of the beam element on which the $i$th force is applied, and:
where ${\widehat{x}}_{i}$ is local coordinate on the beam element for the $i$th brace. ${\mathbf{F}}_{s}\left(t\right)={\left[{F}_{s1}{F}_{s2}\right]}^{T}$ is the vector of static components of the interaction force at the braces due to the weight of the beam and:
where $L$ is the length of the beam and $g$ is the acceleration of gravity.
2.2. Equation of motion of vehicle
The vehicle is modeled by a four DOFs massspring system. The equation of motion of the vehicle is derived using the Lagrange formulation as follows:
where $\mathbf{Y}$, $\dot{\mathbf{Y}}$ and $\ddot{\mathbf{Y}}$ are the vectors of displacement, velocity and acceleration on the vehicle model respectively, and $\mathbf{Y}={\left\{{y}_{V}{\theta}_{V}{y}_{1}{y}_{2}\right\}}^{T}$. Other submatrices and vectors in Eq. (6) have been listed in Appendix 1. It is noted that the modeling of vehicle with point contact model or patch contact model respectively, share the same form of equation of motion in Eq. (6). However submatrices related to the tyres and the interaction forces between pneumatic tyres and road surface included in the force vector $\mathbf{F}\left(t\right)$ will be different.
When the point contact model is adopted, the dynamic interaction forces between the pneumatic tyres and the road surface are:
where ${K}_{t1},{K}_{t2},{C}_{t1},{C}_{t2}$ are the stiffness and damping for the point contact model of the two tyres respectively; $r\left({x}_{1}\right)$, $r\left({x}_{2}\right)$ are the road surface roughness at the two contact points respectively. The overdot (∙) denotes the differentiation with respect to time $t$ and the right prime (′) denotes the differentiation with respect to space coordinate $x$ in this paper. Since the modeling of point contact in the vehicle model may result in an overestimation of the dynamic response on the beamshaped package [10], the patch contact model will also be introduced.
In Fig. 1b the patch contact model [10] is presented, where ${x}_{i}\left(i=\mathrm{1,2}\right)$ denote the location of the $i$th tyre on the road, ${l}_{txi}\left(i=\mathrm{1,2}\right)$ and ${\Delta}_{i}$ are the length of patch contact and the deformation at the center of the $i$th tyre respectively. ${U}_{txi}$ is the deformation at location $x$ on the $i$th tyre:
where $R$ is the radius of tyre, and:
The central deformation ${\Delta}_{i}$ and the corresponding length of patch contact ${l}_{txi}\left(i=\mathrm{1,2}\right)$ when the beamvehicle system is at rest can be obtained according to the practical inflation pressure of the $i$th tyre. The implementation procedure [12] is as follows:
(1) To determine the average deformation of all the tyres denoted as $\Delta $ with standard inflation pressure $P$ under the total weight of the beamvehicle system $G$:
where $n$ and ${l}_{b}$ denote the total number of the tyres on vehicle and the width of the tyres respectively. The standard inflation pressure $P$ is generally equal to $280\mathrm{K}\mathrm{P}\mathrm{a}$.
(2) To estimate the practical inflation pressure ${P}_{wi}$ according to the average deformation of tyre (in the case of a twoaxle vehicle):
where ${Q}_{i}\left(i=\mathrm{1,2}\right)$ are the static vehicular axle loads, $a=22.54$, $b=21.54$ are two constants for heavy vehicles.
(3) To calculate the parameter factor ${\alpha}_{i}$ and stiffness factor ${\beta}_{i}$ from the practical inflation pressure ${P}_{wi}\left(i=\mathrm{1,2}\right)$:
(4) Finally the central deformation ${\Delta}_{i}$ and the corresponding length of patch contact ${l}_{txi}\left(i=\mathrm{1,2}\right)$ of the $i$th tyre on the road surface can be obtained according to Eqs. (14) and (15) repectively:
where $\gamma $ is the ratio between the real contact area and the assumed rectangular contact area, which generally equals to $0.577$.
The axial deformation at location $x$ on the $i$th tyre denoted as ${d}_{i}$ can be calculated as:
When the stiffnesses of the patch contact model for the two tyres are represented by ${K}_{ti}\left(i=1,2\right)$, the relationship between the two static vehicular axle loads and the deformation of tyres can be expressed as:
Therefore when the patch contact model is adopted, the dynamic interaction forces between the pneumatic tyres and the road surface can be obtained as:
2.3. Equation of motion of the beamvehicle system
With the finite element model of the beamshaped package and the equation of motion of the transporation vehicle, the equation of motion of the beamvehicle system can be derived via the matching of the DOFs at the bracing systems as:
where $\mathbf{M}$, $\mathbf{C}$ and $\mathbf{K}$ are the mass, damping and stiffness matrices of the beamvehicle system respectively; $\mathbf{Z}$, $\dot{\mathbf{Z}}$ and $\ddot{\mathbf{Z}}$ are the vectors of displacements, velocities and accelerations respectively, and $\mathbf{Z}={\left\{\mathbf{X}\mathbf{}\mathbf{Y}\right\}}^{T}$; $\mathbf{P}\left(t\right)$ is the vector of forces applied on the beamvehicle system. Other submatrices and vectors in Eq. (21) have been listed in Appendix 2.
3. Numerical  simulations
To investigate the vibration level on the beam subjected to road excitation as well as the influence of the different factors on the results in dynamic analysis, numerical examples are conducted. The parameters of the beamvehicle system are shown in Table 1.
Modal analysis is conducted on the beam model with bracing systems. The damping ratio for each mode ${\zeta}_{i}$ is assumed to be $0.02$ and the constants in the Rayleigh damping model can be calculated from the results of the modal analysis. When the stiffnesses of the bracing systems are large, e. g. larger than ${10}^{8}\mathrm{N}/\mathrm{m}$, the vehicle and beam will have equivalent vibration levels at the connecting locations. Therefore ${10}^{8}\mathrm{N}/\mathrm{m}$ is adopted for the stiffness of the braces in numerical simulations except in Section 3.4 where the influence of the stiffness of the braces on the vibration level of the beam is discussed. When there are no dampers installed in the bracing systems of the beam, the damping may be relatively small and will be negleced in this study. By solving Eq. (21) with the Newmark$\beta $ method, the vibration response on the beamvehicle system can be obtained, in which the samples of the road surface roughness in the excitation are generated according to Eq. (3.2) as shown in Appendix 3. The beam is evenly divided into four beam elements and the sampling rate in numerical simulations is $200\mathrm{H}\mathrm{z}$ which covers the first five natural frequencies of the beamvehicle system.
Table 1. Parameters of the beamvehicle model
Beam model with braces

Vehicle model [13]


$L=8.66m$

${I}_{v}=1.47\times {10}^{5}kg/{m}^{2}$

${m}_{2}=1000kg$

$EI=3\times {10}^{7}{Nm}^{2}$

${m}_{v}=1.7735\times {10}^{4}kg$

${K}_{s1}=2.47\times {10}^{6}N/m$

$\rho A=230.95kg/m$

$S=4.27m$

${C}_{t1}=3.74\times {10}^{6}N/m$

${l}_{1}=7.05m$

${a}_{1}=0.519$

${C}_{s1}=3.0\times {10}^{4}N/m/s$

${l}_{2}=3.577m$

${a}_{2}=0.481$

${C}_{t1}=3.9\times {10}^{3}N/m/s$

${k}_{1}=1\times {10}^{8}N/m$

${L}_{0}=4m$

${K}_{s2}=4.23\times {10}^{6}N/m$

${k}_{2}=1\times {10}^{8}N/m$

${b}_{1}=0.519$

${K}_{t2}=4.6\times {10}^{6}N/m$

${c}_{1}=0$

${b}_{2}=0.481$

${C}_{s2}=4\times {10}^{4}N/m/s$

${c}_{2}=0$

${m}_{1}=1500kg$

${C}_{t2}=4.3\times {10}^{3}N/m/s$

${\zeta}_{i}=0.02$

$R=0.45m$

3.1. Length of the patch contact
The length of footprint between tyre and road surface in the patch contact model will affect the analyzed vibration response of the beamshaped package. The point contact model can be regarded as a special case of the patch contact model with the contact length equals to zero.
In this subsection, a case study is conducted to demonstrate the influence of the contact length on the analyzed vibration response of the beam. When the vehicle carrying a beamshaped package is traveling on a road with Grade B [11] at a speed of $10\mathrm{m}/\mathrm{s}$, the Root Mean Squre (RMS) of the accelerations at the three locations on the beam are varying with the contact length, which is shown in Fig. 2.
Fig. 2. The RMS of acceleration versus contact length
Results in Fig. 2 indicate that analyzed RMS of accelerations on beam decrease when the length of patch contact increases. Results from the bridgevehicle interaction problem [10] show that the calculated bridge response is only sightly larger than the measurement data by using the point contact model, and by adopting the patch contact model the accuracy of the vibration response of bridge can be improved. In contrast, the vibration response on the beam is greatly affected by the model for simulating the contact between tyre and road surface, e. g. when the patch contact model is adopted with a contact length equal to $16\mathrm{c}\mathrm{m}$, the RMS of acceleration on location #2 is equal to $0.172\mathrm{m}/{\mathrm{s}}^{2}$, which is only 1/16 of that in the case of point contact with the contact length equal to zero. Therefore to introduce the patch contact model for tyres in the vibration analysis of the beamvehicle system is of greater significance than in the vibration analysis of the bridge under moving vehicular axle loads.
In order to improve the accuracy of the analyzed vibration response of the beam, apropriate length for the patch contact should be selected in the modeling of interaction between tyres and road surface. Due to the deformation of the tyres, the contact length will vary with the dynamic interaction forces between tyres and road surface, i. e. the contact length is timedependent. The dynamic interaction forces are generally varying around a mean value due to the excitation nature of the irregular road profile. To simplify the problem, the contact length in this study is assumed to be a contant which can be calculated according to Eq. (15) from the deformation of the two tyres due to the static vehicular axle loads as shown in the right sides of Eqs. (17) and (18) respectively. In this study the length of the patch contact equals to 16 cm and this value will be adopted in further numerical simulations.
3.2. The velocity of vehicle
The velocity of the transportation vehicle will affect the vibration level on the beam structure. Assuming the road grade is B, the RMS of the vibration response at different locations on the beam varying with the vehicle’s velocity is shown in Fig. 3.
Fig. 3. RMS of response on beam versus the velocity of vehicle
a) The point contact model
b) The patch contact model
Results show that the analyzed RMS of the vibration response on the beam with point contact model is much larger than that from the patch contact model. In both cases the maximum and minimum RMS of dynamic displacement appears on Locations #1 and #2 respectively, the RMS of velocities on Locations #2 and #3 are equalvalent with different velocities of vehicle, the maximum and minimum RMS of accelerations appear on Location #1 and #3 respectively. The analyzed RMS of the response increases with fluctuations when the velocity of vehicle increases which can be explained as: when the vehicle moves faster, according to Eqs. (7) and (8) for the point contact model or (19) and (20) for the patch contact model, the upper limit of frequency in the excitations increases, which will result in a growth trend in the RMS of the response. The flucturations come from the resonance between the road excitation at a particular velocity of vehicle and the natrual frequencies of the beamvehicle system.
3.3. The road condition
Table 2. RMS of response at different locations with different road grades
Grade of road

A

B

C

D


Location #1

RMS of disp. ($\mathrm{m}\mathrm{m}$)

0.21

0.41

0.83

1.70

RMS of vel. ($\mathrm{m}/\mathrm{s}$)

0.0036

0.0070

0.0141

0.0283


RMS of acc. ($\mathrm{m}/{\mathrm{s}}^{2}$)

0.3573

0.7135

1.4263

2.8550


Location #2

RMS of disp. ($\mathrm{m}\mathrm{m}$)

0.12

0.23

0.47

0.93

RMS of vel. ($\mathrm{m}/\mathrm{s}$)

0.0018

0.0035

0.0071

0.0143


RMS of acc. ($\mathrm{m}/{\mathrm{s}}^{2}$)

0.0881

0.1760

0.3543

0.7087


Location #3

RMS of disp. ($\mathrm{m}\mathrm{m}$)

0.09

0.17

0.35

0.70

RMS of vel. ($\mathrm{m}/\mathrm{s}$)

0.0019

0.0038

0.0075

0.0150


RMS of acc. ($\mathrm{m}/{\mathrm{s}}^{2}$)

0.1324

0.2657

0.5306

1.0593

The unevenness of the road profile is the main source of the vibration problems during transportation. To investigate the influence of level of road surface roughness on the vibration level of the beam with the patch contact model adopted, the RMS of the vibration responses of the beam under different road conditions when the velocity of the vehicle is $10\mathrm{m}/\mathrm{s}$ are listed in Table 2. It is noted from the ISO standard that there are mainly four classes of roads graded by A, B, C and D with the roughness coefficient ${G}_{d}\left({n}_{0}\right)$ equal to $16\times {10}^{6}$, $64\times {10}^{6}$, $256\times {10}^{6}$ and $1024\times {10}^{6}{\mathrm{m}}^{3}$ which represent the road conditions as good, average, bad and worse respectively. Results in Table 2 show that the vibration level on the beam increases sigificantly when the road conditions become worse.
3.4. Stiffness of the bracing system
The moving vehicle and the beamshaped structure on top are connected by braces modeled by springdamping systems. The stiffness of the braces will affect the propogation of vibration on vehicle body to the beam. This study aims to provide the influence of the stiffness of the brace systems on the vibration level of the beam. To simplify the problem, the stiffnesses of the two braces are assumed to have the same values. When the vehicle travels on a road with Grade B at a speed of $10\mathrm{m}/\mathrm{s}$, the maximun displacements and the RMS of accelerations at different locations on the beam varying with stiffness of braces are shown in Fig. 4 when the patch contact model is adopted.
Fig. 4. Response on beam versus the stiffness of the braces
a) The maximum displacement
b) The RMS of acceleration
On the perspective of design, the maximum displacement on beam should be small in order to fix the beamshaped package which requires a relatively large stiffness of the braces. On the other hand, when the braces are very stiff, the vibration level on vehicle and beam at the location of braces tends to be identical, i. e. there is no vibration reduction when the vibration on the vehicle body propogates to the beam on top through the bracing systems, which should also be avoided in the design process. The design of the stiffnesses of bracing systems requires to assure the beam to have both small displacements and low vibration level. Results in Fig. 4 show that the maximum displacement on the beam has a sharp decrease to a converged value when the stiffness of bracing systems increases, the RMS of the acceleration on beam converges to a constant value before several fluctuations and the optimal stiffness can be obtained when the RMS reaches the minimum, which equals to $1\times {10}^{8}\mathrm{N}/\mathrm{m}$ in this study.
4. Conclusions
The equation of motion of a beamshaped package on moving vehicle under the road excitation is derived in this paper and the vibration analysis is further conducted, which aims to provide references on the condition assessment and vibration control in package transportation. Two different contact models, namely the point contact model and the patch contact model repectively, are adopted to simulate the interaction between tyre and road surface. Numerical simulations are conducted to investigate the influence of the length of patch contact, the velocity of vehicle, the road condition and the stiffnesses of the bracing systems on the vibration level of the beamshaped structure. Some conclusions are drawn:
(1) The vibration response on the beam is greatly affected by the modeling of the contact between tyre and road surface, the analyzed RMS of accelerations on the beam decreases when the length of patch contact increases. Since the point contact model will lead to an overestimation on the calculated response, to introduce the patch contact model for tyres in the vibration analysis of the beamvehicle system is of great significance.
(2) The length of the patch contact model should be carefully selected by the practical vehicular axle loads and the specification of tyres, to guarantee the accuracy of the analyzed vibration response on beamshaped package.
(3) The analyzed RMS of the vibration response on package increases with fluctuations when the velocity of vehicle increases or the road conditions become worse.
(4) When the stiffness of brace systems increases, the maximum displacement on the beam has a sharp decrease to a converging value, the RMS of the acceleration on beam converges to a constant value before several fluctuations. The opitimal value of the stiffnesses of bracing systems requires the beam to have both small displacements and low vibration level.
Acknowledgements
The work described in this paper was supported by a Research Grant from National Natural Science Foundation of China (No. 10902024), Jiangsu Natural Science Foundation (No. BK2010397) and also funded by the Priority Academic Program Development (PAPD) of Jiangsu Higher Education Institutions.
References
 Shi X. M., Cai C. S. Simulation of dynamic effects of vehicles on pavement using a 3D interaction model. Journal of Transportation Engineering, Vol. 135, Issue 10, 2009, p. 736744. [Search CrossRef]
 Michialtsos G. T. Bouncing of a vehicle on an irregularity: a mathematical model. Journal of Vibration and Control, Vol. 16, Issue 2, 2010, p. 181206. [Search CrossRef]
 Lu F., Ishikawa Y., Kitazawa H., Satake T. Effect of vehicle speed on shock and vibration levels in truck transport. Packaging Technology and Science, Vol. 23, 2010, p. 101109. [Search CrossRef]
 Lak M. A., Degrande G., Lombaert G. The effect of road unevenness on the dynamic vehicle response and groundborne vibrations due to road traffic. Soil Dynamic and Earthquake Engineering, Vol. 31, 2011, p. 13571377. [Search CrossRef]
 Barbosa R. S. Vehicle vibration response subjected to longwave measured pavement irregularity. Journal of Mechanical Engineering and Automation, Vol. 2, Issue 2, 2012, p. 1724. [Search CrossRef]
 Ghaith F. A. Nonlinear dynamic modeling of elastic beam fixed on a moving cart and carrying a lumped tip mass subjected to external periodic forces. Proceedings of the ASME 2009 International Design Engineering Technical Conferences & Computers and Information in Engineering Conference, San Diego, California, USA, DETC200986318, 2009, p. 17. [Search CrossRef]
 Ragulskis K., Gegeckienė L., Kibirkštis E., Miliūnas V., Zubrickaitė L., Pauliukas A., Ragulskis L. Dynamic study of transportation containers with packages. Journal of Vibroengineering, Vol. 14, Issue 4, 2012, p. 18761884. [Search CrossRef]
 Gwanghum G. Vehicle Dynamic Simulation with a Comprehensive Model for Pneumatic Tyres. Ph. D. Thesis, University of Arisona, 1988. [Search CrossRef]
 Lyasko M. I. The determination of deflection and contact characteristic of a pneumatic tyre on a rigid surface. Journal of Terramechanics, Vol. 31, Issue 4, 1994, p. 239246. [Search CrossRef]
 Yin X. F., Cai C. S., Fang Z., Deng L. Bridge vibration under vehicular loads: tyre patch contact verse point contact. International Journal of Structural Stability and Dynamics, Vol. 10, Issue 3, 2010, p. 529554. [Search CrossRef]
 ISO 8606:1995(E). Mechanical Vibration – Road Surface Profiles – Reporting of Measured Data. 1995. [Search CrossRef]
 Ge F., Qu X. C., Yang S. H., Chen D. K., Zhu Y. Study of selection and relation between tyre air pressure and automobile performance. Journal of Academy of Military Transportation, Vol. 11, Issue 1, p. 5355, (in Chinese). [Search CrossRef]
 Wu S. Q., Law S. S. Vehicle axle load identification on bridge deck with irregular road surface profile. Engineering Structures, Vol. 33, 2011, p. 591601. [Search CrossRef]