Abstract
Asymmetrical longwaved road unevenness can accelerate the fatigue failure of the dump truck frame. Due to the difficulty of the tire damping modelling and the complex pitch motion of balanced suspension, it is hard to reveal the relationship between the parameters of longwaved road and the fatigue life of dump truck frame. Firstly, on the premise of not introducing truncation error, by applying matrix operations, an 11 DOF vibration model of dump truck considering tire damping is established to simulate the dynamic responding loads of frame bearing points. Secondly, the proposed model is validated by comparing the predicted and measured vehicle response and subsequently used to predict the dynamic vehicle loads. Finally, using the simulated responding loadtime histories, the paper calculates and analyzes the fatigue life of dump truck under different road profile excitations. The results show that the increase of road amplitude can lead to fatigue weak area expanding to the back part of frame, and logarithmically shorten the fatigue life of frame.
1. Introduction
The excitation imposed on running vehicle is mainly road surface roughness [1]. While the interaction between the vehicle vibration and road roughness is a complex dynamic process. The dynamic responses depend on road roughness, suspension properties, vehicle speed and other factors simultaneously, which cannot be simply calculated via linear operation from road roughness.
Wheel force transducer (WFT) is a tool which can measure the wheel force effectively, which has been used in many vehicle test fields. However, due to the high price and poor versatility, the popularization and application of the WFT in the field of vehicle dynamics is limited [2]. In order to obtain more reliable dynamic responses, researchers have made several studies on relationship between road excitation and wheel loads.
H. M. Ngwangwa et al. [3] using an artificial neural network reconstructs the road surface profiles from measured vehicle accelerations. Jinshan Lin [4] applies a new cubic spline weight function neural network (CSWFNN) to the identification of road surface power spectrum density from measured accelerations. The introduction of these two neural network methods basically achieved the recognition of road roughness, and built the relationship between road roughness and acceleration signals respectively. However, these neural network methods cannot give the physical interpretation for the intrinsic relationship. A. González et al. [5] notes that the relationship between the power spectral densities of road surface and vehicle accelerations can be interpreted via a transfer function, but their vehicle model is oversimplified and without detailed study of dynamic behavior. M. A. Lak, G. Degrande and G. Lombaert [6] establish a multiDOF dynamic vehicle model to study the effect of road unevenness on the dynamic vehicle response.
There are two calculation methods of vehicle vibration response: timedomain method and frequencydomain method. For random response of linear vehicle system, frequencydomain method is often used. Wherein, the timedomain model considering tire damping parameter is rarer. On one hand, Yu Fan and Lin Yi [7] propose that vehicle system dynamics is mainly confined to rigid body dynamics of suspension, the research frequency of which limits to 0.25~15 Hz, meaning the shortest wavelength of road unevenness, $\lambda $, is 1/15 of vehicle driving distance per second. Generally, within the scope of normal driving speed, the shortest wavelength is largely greater than tire print, approximately 0.15~0.25 m. Therefore, vertical dynamic characteristics of rolling tire can be simplified to a spring ignoring damping. On the other hand, the reason to ignore tire damping in vehicle timedomain modeling is that differential equations containing derivative inputs cannot be easily converted to state space equations, which leads to a difficulty in solving analytical solutions of simultaneous ordinary differential equations (simultaneous ODE). For a linear system with one input and one output, i.e., singleinput singleoutput (SISO) system, the solving method for analytical solutions is complicated, let alone the conversion from differential equations containing derivative inputs into state space equation in multiinput multioutput system.
Compared to other commercial vehicles, the working speed of dump truck, especially mining dump truck, is relatively lower. The road conditions, especially near the loading and unloading sites, which are the beginning and ending of cycling mileage, are the toughest. The proportion of lowfrequency and longwaved undulations in these fields is higher than normal roads.
In order to avoid the wheel off the ground and uneven distribution of axle loads, balanced suspension has been widely used in threeaxle dump truck to connect the second and rear axle. Due to balanced suspension rotating, while one wheel is being lifted, the other wheel will be lowered. If two arms of traction bars have the same length, equal vertical loads in two axles will be guaranteed in any case. However, under the excitation of longwaved road, when the resultant force of longitudinal force components of front and rear axles does not equal to zero, the frame will suffer horizontal force from the balance shaft seat, causing different influences on the fatigue life of frame and balanced suspension.
This paper analyzes the fatigue life of 6×4 dump truck frame excited by low frequency longwaved road roughness. In order to obtain more actual loadtime histories of frame bearing points, an 11DOF vehicle vibration model has been established, considering pitch motion of balanced suspension. Considering the low coherence between the left and right road roughness, and excessive mutual disturbance caused by asymmetric road excitation, independent suspension model could not be used in traditional multiDOF modelling. Therefore, three integral axles of the dump truck are modelled as dependent suspension [8].
The tire damping of front wheel takes 10~35 % of the front suspension damping, and tire damping of latter wheels almost reached the 3~7 times of the balanced suspension damping [9]. Ignoring the tire damping in vehicle vibration modeling will cause obvious distortion [10, 11]. On the basis of conversion from SISO differential equations to state space, via introduction of matrix operations, the paper implements the solving for analytical solution of multivariate differential equations containing derivative inputs, and then completes the modelling of vehicle vibration system considering tire damping.
The paper builds a full dump truck vibration model considering tire damping to simulate the dynamic responding loads of frame bearing points in Section 2 and Section 3. Then the proposed model is validated by comparing the predicted and measured response results and subsequently used to predict the dynamic vehicle loads for the fatigue life calculation in Section 4. Section 5 calculates and analyzes the fatigue life of the dump truck frame under different road profile excitation. The article is concluded in Section 6.
2. The establishment of timedomain wheel excitation model
2.1. The timedomain modeling of random road profile
In the ISO classification, the relationship between the displacement power spectral density ${G}_{q}\left(n\right)$ and the wave number $n$ for different classes of road roughness may be approximated by (ISO 8608 [12]):
where, ${G}_{q}\left(n\right)$ is expressed in m^{3} cycle^{1}, $n$ is expressed in cycles m^{1}, ${G}_{q}\left({n}_{0}\right)$ is the power spectral density of road roughness at the datum or cutoff wave number ${n}_{0}$, which is equal to 1/2π cycles m^{1}. $W$ is the frequency index, which reflects the frequency structure of the pavement, usually, $W=\text{2}$. The ${n}_{1}$ and ${n}_{2}$ represent the upper and low limit of frequency bandwidth respectively. The bandwidth (${n}_{1}$, ${n}_{2}$) should ensure that it includes the main natural frequency of the vehicle when the dump truck with average speed on the road, here, ${n}_{1}=$0.005 m^{1}, ${n}_{2}=$1.422 m^{1}.
The displacement power spectral density as given in Eq. (1) is usually calculated from the measurement of surface roughness described by vertical ordinates at equally spaced points along the road. However, in the absence of such measurements, pseudorandom profiles can be generated to fit those spectral densities [3]. Ren Hongbin, Chen Sizhong and Wu Zhicheng [13] present a formula for generating a onedimensional random profile as:
where, ${q}_{r}$ is the discrete value of the random road profile in space domain. ${G}_{q}\left(n\right)$ is the discrete value of PSD function of the road roughness. $\Delta l$ is the distance interval between successive ordinates of the surface profile. In order to take pavement phase changes into consideration, ${\phi}_{n}$ was introduced into the equations. Only the real part was picked into the road model, ignoring the imaginary part, which will not cause complex function. ${\phi}_{n}$ is a set of independent random phase angles uniformly distributed between 0 and 2π. It is the key parameter of the coherence of left track and right track.
According to the coherence function of left track and right track, the inputs of left and right track has some characteristics: the correlation of their amplitude and phase decreases with the increase of the spatial frequency, and it fits the fact.
As suggested in [13], the empirical coherence functions of left track and right track is exponentially decreasing:
where, $\rho $ is empirical coherence parameter [13], here, $\rho =$1.5. ${d}_{v}$ is the width of two wheels on an axle.
We can calculate the phase of the right track ${\phi}_{n\_right}$ as follows:
where, ${\phi}_{n\_left}$ is the phase of left track. ${\phi}_{n\_new}$ is another new random number that obey the normal distribution of the interval of [0, 2π].
According to Eq. (4), phase information of random road roughness of the right road could be obtained, from which, it could be known that the coherence of ${\phi}_{n\_right}$ and ${\phi}_{n\_left}$ decreases as spatial frequency increases, which could guarantee the satisfaction of left and right wheel coherence.
The values of random road profile in time domain could be obtained by the following equation:
where, ${q}_{r}\left(l\right)$ is the road roughness along the travelling direction in the space domain. $v$ is the running speed of vehicle.
Suppose there is only spatial hysteresis caused by wheelbase between the front wheel and rear wheel road profiles. The random inputs of left second wheel and rear wheel [14] can be given by:
where, ${\tau}_{1}=L1/v$, ${\tau}_{2}=L2/v$. $L1$ is the distance of the wheelbase between front axle and second axle. $L2$ is the distance of the wheelbase between front axle and rear axle.
Based on the algorithm stated above, a sixwheel excitation model for dump truck is established via Matlab/Simulink. Take the class B road as an example. ${G}_{q}\left({n}_{0}\right)=\text{64\xd71}{\text{0}}^{\text{6}}\text{}{\text{m}}^{\text{3}}$, $B=$ 2 m, $L1=$ 3.8 m, $L2=$ 3.8+1.35 m, the vehicle’s velocity is $v=$2.78 m/s. The random road profiles for six wheels in space domain could be reached as shown in Fig. 1.
Fig. 1Random road profiles for six wheels of dump truck
2.2. The timedomain modelling of longwaved road profile
In order to quantitatively investigate the impacts of longwaved lowfrequency road topography on the vehicle vibration and fatigue life of frame, it is necessary to superimpose deterministic longwaved bosses on the standard random road profile, taking standard road profile as both background and reference, to build a fullband asymmetric boss test road. Thus, by superimposing a set of bosses, the crosssection geometric parameters of which were known, an asymmetric undulating road sample similar to the real test road is built.
As it is shown in Fig. 2, longwaved road is constituted by several bosses placing on a class B road. The geometrical parameters of the boss crosssection are shown in Fig. 2(a). The lift (between OA, which is 00.6 m) and return (between BC, which is 1.21.8 m) trip of undulating road are all circular curves, whose radius are $R=$1.275 m, which respectively tangent to point A and B. The upper surface of the boss is a horizontal straight line, with $H=$0.15 m. The mathematical expression of the boss road profile can be given as follows:
where $x$ represents the distance along the test road, ${q}_{B}$ is the boss’s vertical ordinate at distance $x$. As shown in Fig. 2(b), bosses are arranged alternately left and right on the road. The longitudinal separation between every boss and the adjacency one to the other side is 0.85 m.
Fig. 2The research sample of longwaved road profile
a) The geometrical parameters of the boss crosssection
b) The boss test road
Combined with the random road profile ${q}_{r}$ gained in section 2.4, the fullband longwaved road profile ${q}_{l}$ to be reconstructed was given by the summation of Eqs. (2, 68). As shown in Fig. 3.
Fig. 3The fullband longwaved road profile
3. The 11DOF vehicle vibration model of dump truck
The entire system was modelled by the linear 11degree of freedom (11DOF) model comprising 3DOFs associated with the vehicle body, 6DOFs associated with the bounce and roll motions of three integral axles, and 2DOF for the pitch motions of the balanced suspension.
This paper mainly studies the dynamic loads of frame bearing points. Therefore, the cab and other parts above the frame are simplified to the equivalent mass of the centroid point.
For most vehicles, resonances in roll occur at lower frequencies (typically between 0.5 and 1 Hz) than resonances in bounce (Gillespie [15]), which is easy to study the vibration characteristics of vehicle using the 1/2 model. However, the coherence of the left and right wheel tracks of asymmetric longwaved road tend to be smaller and the vehicle is more responsive to roll input. Therefore, it is necessary to establish the whole vehicle model rather than a pitchplane model.
Three axles in 6×4 dump truck are all integral axles, where the axle is a kind of dependent suspension. There is a direct geometric relationship between bounce motions of left and right axle ends, which could not be simplified as two independent vertical degrees of freedom [8]. The dependent model must be introduced, as shown in Fig. 4.
Fig. 411DOF vibration model of dump truck. Where, Zb is the bounce motion of vehicle body, θ is the pitch motion of vehicle body, and φ is the roll motion. Zi (i=1~6) represents the bounce motion of axle end. Zj (j=A, B, C, D) represents the bounce motion of frame bearing point. Zi and θi (i=f, m, r) are the bounce and roll motions of integral axle. Zj (j=E, F,G, H) represents the bounce of stabilizer rod end
The balanced suspension has been widely utilized in threeaxle dump truck to enhance the axle load; therefore, it is necessary to build a reasonable balanced suspension model and demonstrate its vibration characteristics [16]. As shown in Fig. 4, the leaf springs of the balanced suspension are regarded as a whole. The middle part of the leaf springs is connected with the sprung mass. Both ends of the leaf springs are separately linked with the front and rear axle through the rigid stabilizer rod. The front and rear axle's vibration are coupled together with the stabilizer rod.
The equation of motion of a linear vehicle model (with linear geometry and linear springs and dampers representing tire and suspension elements) can be written in matrix form as:
where $\mathbf{M}$ is the matrix representing vehicle body mass, $\mathbf{C}$ is a matrix of vehicle damping, $\mathbf{K}$ is the stiffness matrix of the vehicle system, ${\mathbf{C}}_{\mathbf{t}}$ is a matrix of tire damping, ${\mathbf{K}}_{\mathbf{t}}$ is the stiffness matrix of the tire, $\mathbf{Z}$ is the vector of vehicle dynamic responses and ${\mathbf{Q}}_{\mathbf{g}}$ is the vector of displacement inputs acting on the tires.
It was assumed that both the pitch angle and roll angle were so small that kinematical motions at the front axle, rear axle, second axle and balanced suspension mount might be approximated by the following equations, respectively:
Vertical bounce of axle ends could not be directly obtained because of the DOF definition of balanced suspension. Based on the small deflection hypothesis, the following equations could be approximately gained.
In order to simplify the solving for timedomain solution, to obtain a convenient and compact way to model and analyze systems with multiple inputs and outputs, the multivariable differential equations should be converted to the state space representation, which is a mathematical model of a physical system as a set of input, output and state variables related by firstorder differential equations. To abstract from the number of inputs, outputs and states, the variables of linear and time invariant system are expressed as vectors, and the differential and algebraic equations can be written in matrix form. It is essential that the matrix and vector modeling is very efficient from a computational standpoint for computer implementation.
Consideration of the tire damping in vehicle vibration model will put time derivatives of road roughness in the right hand side of Eq. (10). Since most of the control theory is done for systems without considering the input derivatives, it is necessary to lower the order of inputs, i.e. to remove the time derivatives of inputs [17].
The conventional order reduction approach is a direct numerical differential operation to road roughness, to obtain time derivatives of road roughness. Firstly, the road roughness input ${\mathbf{K}}_{\mathbf{t}}{\mathbf{Q}}_{\mathbf{g}}+{\mathbf{C}}_{\mathbf{t}}{\dot{\mathbf{Q}}}_{\mathbf{g}}$ was directly converted to force input acting on the vehicle $\mathbf{f}\left(t\right)$ through numerical differential operation. Then Eq. (10) was recast into multivariable state space equations given by:
where $\mathbf{I}$ and $0$ represent identity and null matrices with similar dimensions as the property matrices, $\widehat{0}$ is a null column vector with a length similar to the forcing vector $\mathbf{f}$ and $\mathbf{x}$ is a vector containing the states (displacement and velocity vectors). And the state variables $\mathbf{x}$ was selected as Eq. (26):
Although input order reduction has been realized in Eq. (25) via the conversion from ${\mathbf{K}}_{\mathbf{t}}{\mathbf{Q}}_{\mathbf{g}}+{\mathbf{C}}_{\mathbf{t}}{\dot{\mathbf{Q}}}_{\mathbf{g}}$ to $\mathbf{f}\left(t\right)$, the essence of this order reduction is an approximate computational process of numerical differentiation algorithm based on RungeKutta. Order reduction itself could cause truncation error input and reduce the reliability and accuracy of system output. In order to oversimplify input vector, therefore, it is essential to keep input derivative for solving reliability and accuracy.
At present, the transformation from SISO differential equation containing input derivative to state space representation is achieved mainly by selecting an appropriate intermediate variable. Based on the SISO transformation method, the input, output variable and status variable are expanded to vector, variable of each system is expanded to coefficient matrix and matrix operations is introduced to the conversion process. Then, input derivatives of multivariable equations, ${\dot{\mathbf{Q}}}_{\mathbf{g}}$, could be eliminated successfully. Statespace equation could be achieved, as shown in Eq. (27):
where the state variables ${x}_{2}$ was selected as Eq. (18):
Via a ss.m function in MATLAB, a 11DOF statespace model of dump truck was constructed based on Eq. (27).
Based on the statespace model, then the responses of the dynamic system excited by the input ${Q}_{g}$ assigned by the timedomain road roughness in section 2.4 and 2.5 were simulated via the lsim.m function in MATLAB.
4. Numerical simulation and validation
The acceleration and load signals of specific/arbitrary position could be obtained via numerical simulation by the proposed statespace model. Before further calculations, the 11DOF model needs to be validated and calibrated by comparing the simulated and measured vehicle responses.
4.1. Timedomain simulation of dynamic loads
The dynamic wheel load is defined as the tire load changes relative to static equilibrium. ${F}_{di}$ ($i=1,2,3,4,5,6$) is the dynamic tire load of each axle end in 11DOF dump truck. Therefore, the six tire loads can be given by:
Dynamic load of frame bearing point is the frame load changes imposed by bearing points. There are four bearing points for the 11DOF dump truck frame: 2 points on front spring supports (F_{SA} and F_{SB}) and each point on left and right of suspension (${F}_{SC}$ and ${F}_{SD}$). In Section 3.2, leaf spring of balanced suspension has been simplified to two springs, as shown in Fig. 4. The introduction of virtual stabilizer rod in the balanced suspension model results in the emergence of the intermediate variables, i.e., the hypothetical dynamic loads of stabilizer rod end, ${F}_{SE}$, ${F}_{SF}$, ${F}_{SG}$ and ${F}_{SH}$. The formula of front bearing points of frame and four dynamic loads of rod ends is shown in Eq. (30) and the relationship between dynamic loads of rod ends and that of balance shaft ends is shown in Eq. (31):
Generally speaking, the longitudinal force imposed by balanced suspension on vehicle frame is not significant when driving steadily on a class B road or above. When the road surface slightly undulates, which does not cause significant deflection of leaf spring, the longitudinal dynamic force generated by driving force fluctuations will be delivered to the 4th cross member and the bottom of balance shaft via traction bar.
The introduction of balanced suspension effectively improved load distribution between second and rear axles, moderated axle dynamic load shocks and enhanced axle service life. In the meantime, however, the structure principle of integral balanced suspension decides that, when the resultant force of longitudinal force components of front and rear axles does not equal to zero, the frame will suffer horizontal component force exerted by offroad road roughness. Computing conclusion shows that the root mean square of frame longitudinal dynamic load, under the excitation of the bosses road with $H=$0.15 m, can range up to 3 % to 5 % of that of vertical load, which cannot be ignored in analysis methods.
As it is shown in Fig. 5, longwaved road profile can cause the contact force between tire and ground deviating from vertical direction. The magnitude of horizontal component force depends on the angle between road profile and horizontal plane. The sum of vertical component forces of the second axle and rear axle equals to the burden of vehicle load on balance shaft, which is ${F}_{C}={F}_{EV}+{F}_{FV}={F}_{d3}+{F}_{d5}$. Whereas, the horizontal loads exerted by balance axle seat equal to the sum of component forces from two tires, which is ${F}_{Ch}={F}_{Eh}+{F}_{Fh}$. During the transfer process from tire contacting force to frame bearing force, the leaf spring and traction bar all played a delivering role. The resultant force of frame exerted by balance shaft always follows the mechanical analysis conclusion stated above.
Fig. 5Mechanical analysis of balanced suspension
Equivalently shift the horizontal force delivered by traction bar to the axis of the balance shaft. The horizontal loads on two sides of balance shaft, ${F}_{CH}$ and ${F}_{DH}$, could be written as:
Considering the envelope effect from tire to road, the tire grounding angle could not be directly solved via simple derivation operation from road displacement. Contacting print of tire and road profile can be equivalent to a straight line, $L$. Then by calculating angle ${\theta}_{i}$$(i=E,F,G,H)$, as shown in Fig. 6, the more reliable tire grounding angle, ${\theta}_{E}$, ${\theta}_{F}$, ${\theta}_{G}$, ${\theta}_{H}$, could be obtained.
Fig. 6Tire grounding angle
4.2. Road tests for model validation
In order to make the comparison between the predicted and measured vehicle responses, under excitations of class B road and boss test road, two 30meterlong test roads are chosen as the road samples. As shown in Fig. 3(b).
The parameters of the test vehicle are shown in Table 1.
Table 1Parameters of the 6×4 dump truck
Curb weight / (kg)  12815  Front spring stiffness / (N/m)  3.8×105  Front tire stiffness / (N/m)  1.1×106 
Wheelbase / (mm)  3800+1350  Rear spring stiffness / (N/m)  3.6×106  Rear tire stiffness / (N/m)  2.2×106 
Track / (mm)  1850  Front axle damping / (N s/m)  10000  Front tire damping / (N s/m)  3500 
Suspension damping / (N s/m)  2000  Rear tire damping / (N s/m)  7000 
As shown in Fig. 7, in order to acquire the dynamic loads of the front wheels, two wheel force transducers were installed on each side of front axle (Suzhou Brontek Measure & Control Technology Co., Ltd, Model No. SD823). For acquiring the acceleration signals of axle ends and frame bearing points, to verify the secondorder response output of the 11DOF model, acceleration sensors were installed near the leaf springs of front and rear axles, and on the main frame near the installing seat of balance shaft (Sinocera Piezotronics, Inc., Model No. X901).
Fig. 7Sensor arrangement
Fig. 8Power spectral density (PSD) of front wheel load (Fd1), acceleration of front axle (AccF) and acceleration of balance shaft seat (AccBS) under two road profiles
a) PSD of simulated Fd1 under class B road
b) PSD of measured Fd1 under class B road
c) PSD of simulated Fd1 under longwaved road
d) PSD of measured Fd1 under longwaved road
e) PSD of simulated AccF under class B road
f) PSD of measured AccF under class B road
g) PSD of simulated AccF under longwaved road
h) PSD of measured AccF under longwaved road
i) PSD of simulated AccBS under class B road
j) PSD of measured AccBS under class B road
k) PSD of simulated AccBS under longwaved road
l) PSD of measured AccBS under longwaved road
Fig. 8 shows that, from the perspective of energy distribution in the frequency domain, the simulated results are well consistent with the measured data, which proves that the proposed timedomain vehicle model can well describe the dynamic characteristics of dump truck.
Fig. 8(ab) and 8(cd) show that the PSD peak of Fd1 (front wheel load) near the 2.5 Hz, corresponding to sprung mass motion, has been distinctly restrained. Whereas, the low frequency component corresponding to frontaxle unsprung mass is released, and shows an obvious aggregation in low frequency band. It presents that under longwaved road excitation, front suspension is more susceptible to low frequency road, and the influence from sprung mass motion has been correspondingly weakened. The PSDs of acceleration signals show the similar conclusions, as shown in Fig. 8(ef) and 8(gh).
As shown in Fig. 8(ij) and 8(kl), the PSDs of acceleration on balance shaft seat under different road profiles show a higher consistency, which proves that the balanced suspension can well absorb the shock to frame caused by longwaved road profile, and eases the lowfrequency vibration of balance shaft seat.
The delicate difference lies in the aggregation of responding load between 25 Hz caused by asymmetrical longwaved road, which illustrates that roll motion caused by asymmetrical road could not be restrained by the balanced suspension. In other words, the influence of roll motion to balanced suspension vibration mainly embody at 2~5 Hz.
5. Fatigue life analysis
According to Eq. (2022), six loadtime histories of frame bearing points, loads of balanced suspension seat containing vertical load and horizontal load, could be gained. After a series of preprocessing, the loadtime histories were imported into the fatigue life calculation and damage analysis software, nCode, to anticipate and analyze fatigue life of frame. Several conclusions are drawn as follows.
Fig. 9The relationship among the fatigue life of weak area, wave amplitude and wavelength ratio in total distance
In order to make a quantitative research on influence of parameters of longwaved road to fatigue life of frame, take the wave amplitude, $H$ and the ratio of wavelength to total distance, $\lambda /L$, as variables separately and take responding loads under various road profile as the inputs for fatigue life calculation. The relationship among the fatigue life of weak area, wave amplitude and wavelength ratio can be drawn as Fig. 9.
Wherein, the red dashed line, between 10^{10} and 10^{12} km, is the frame fatigue life under class B road, which is treated as the life criterion for research on other working conditions.
As shown in Fig. 9, fatigue life of frame is shortened linearly with proportion of boss length in total distance increasing. Nevertheless, with the increase of boss amplitude, the shortening of fatigue life shows a logarithmic trend. Therefore, in Fig. 9, the logarithmic scale is adopted in vertical coordinate.
Take the road with amplitude = 0.12 m as an example, when the ratio of boss wavelength increases from 6 % to 24 %, fatigue life of fatigue weak area decreases from 1.391×106 km to 9.212×105 km. And the linear correlation coefficient of four ratios is 0. 9841. Under excitation of different boss amplitudes, $H=$0.15 m, 0.08 m and 0.04, the other three working conditions show the similar linear correlation conclusion, whose linear correlation coefficients are 0.9859, 0.9776 and 0.9639 respectively.
The logarithmic relationship, shown by fatigue life of weak area under different amplitudes, demonstrates that the impact of amplitude of longwaved road is more obvious than that of ratio on fatigue life of frame.
Fig. 10The fatigue life contour of frame under class B road excitation
Fig. 11The fatigue life contour of frame under longwaved road sample 1 (H= 0.04 m)
Fig. 12The fatigue life contour of frame under longwaved road sample 2 (H= 0.12 m)
It can be seen from Fig. 9 and Fig. 10, compared to the fatigue life contour under class B road, the fatigue weak area of frame under longwaved road ($H=$0.04 m, ratio of wavelength is 24 %) has become decentralized mainly to two parts. The one is on main frame side rail near the engine supporting beam and the other is between the third crossmember and front end of mainsub frame connecting plate. These fatigue analysis results are consistent with the statistical conclusions of frame fatigue cracking under mining and construction site [18].
In order to investigate the impact of wave amplitude on frame fatigue life, on the premise of same wavelength ratio, another longwaved road sample and its corresponding loadtime histories were created. Then the fatigue life of frame under the excitation of higher amplitude roads was calculated and analyzed.
As shown in Fig. 13, under the excitation of the second longwaved road sample, whose boss amplitude is 0.12 m, the fatigue life of frame is shortened dramatically and the fatigue weak area began to expand to rear of frame. It can be seen from the bottom view of balanced suspension in Fig. 10 and Fig. 11, with the increase of the wave amplitude, the fatigue weak area of fourth crossmember and balance shaft assemblies increased significantly. The fatigue weak area of other regions did not show much obvious positional deviation, which indicates that the impact of wave amplitude change on frame fatigue life is mainly focused on the main frame near the third crossmember and balance shaft assemblies.
6. Conclusions
This work was aimed at investigating the impacts of longwaved road roughness on fatigue life of dump truck frame by establishing 11DOF vibration model considering tire damping, dependent suspension and balanced suspension.
1. Introduction of matrix operations makes the transformation method of single input single output (SISO) differential equation containing input derivative to state space representation successfully expand to the multivariate equations (MIMO). Therefore, 11DOF dump truck model has been successfully established considering tire damping, on the premise of not introducing truncation error. The model has been validated and calibrated via dynamic wheel forces and acceleration signals of axle end collected in road tests. From the perspective of energy distribution of the signal in the frequency domain, the simulated results are basically consistent with measured data, which proves that the proposed timedomain vehicle model could well describe the dynamic characteristics of dump truck.
2. Under the longwaved road excitation, front suspension is more susceptible to low frequency road roughness, rather than the sprung mass motion. The balanced suspension can well absorb the shock to frame caused by longwaved road profile and ease the lowfrequency vibration of balance shaft seat. Whereas, the balanced suspension structure cannot well suppress the roll motion caused by asymmetrical road.
3. The calculation results show that fatigue life of frame is linearly shortened with the ratio of boss length in total distance increasing. Nevertheless, with the increase of boss amplitude, the shortening of fatigue life shows a logarithmic trend, which demonstrates that the impact of amplitude of longwaved road is more obvious than that of ratio on frame fatigue life.
4. Simulated analysis of frame fatigue life demonstrates that under longwaved road excitation, the fatigue weak area is more decentralized to two parts. The one is on main frame side rail near the engine supporting beam and the other is between the third crossmember and front end of mainsub frame connecting plate. The further analysis shows that the increase of longwaved road amplitude will lead to the fatigue weak area expanding to the back part of frame. The change of longwaved road amplitude has an influence on fatigue life, mainly focused near the third crossmember on the main frame and balance shaft assemblies.
The simulation conclusions and the statistical results of frame fatigue cracking in actual mining area and construction site have a higher consistency, which fully supports the reasonability of this paper.
References

Kim J.H., Song J.B. Control logic for an electric power steering system using assist motor. Mechatronics, Vol. 12, Issue 3, 2002, p. 447459.

Lin G., Zhang W., Yang F., Pang H., Wang D. An initial value calibration method for the wheel force transducer based on memetic optimization framework. Mathematical Problems in Engineering, Vol. 2013, 2013.

Ngwangwa H. M., Heyns P. S., Labuschagne F. J. J., Kululanga G. K. Reconstruction of road defects and road roughness classification using vehicle responses with artificial neural networks simulation. Journal of Terramechanics, Vol. 47, Issue 2, 2010, p. 97111.

Lin Jinshan Identification of road surface power spectrum density based on a new cubic spline weight neural network. Energy Procedia, Vol. 17A, 2012, p. 534539.

González A., O’brien E. J., Li Y.Y., Cashell K. The use of vehicle acceleration measurements to estimate road roughness. Vehicle System Dynamics, Vol. 46, Issue 6, 2008, p. 483499.

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 Dynamics and Earthquake Engineering, Vol. 31, Issue 10, 2011, p. 13571377.

Fan Yu, Yi Lin Vehicle System Dynamics, 1st ed. China Machine Press, 2005.

Hu Q., Wang W., Sun B., Lei X. Vibration property of dependent suspension based on Matlab. Journal of Chongqing Jiaotong University (Natural Science), Vol. 1, 2012, p. 133136.

Wang T., Shen J. Macro Design and Analysis of Truck, 1st ed. Standards Press of China, 2013.

Lu Y. Research on Dynamics of Heavy Vehicle and Road Surface Interaction. Ph.D. Theses, Beijing Jiaotong University, 2011.

Jin Y. Research on the Vibration Performance of a Threeaxle Heavyduty equipped with Balanced Suspension and the Road Damage it caused. M.S. Theses, Jilin University, 2009.

I.T. 108/SC 2, ISO 8608:1995, Mechanical vibration – Road surface profiles – Reporting of measured data. Multiple. Distributed through American National Standards Institute, 2007.

Ren Hongbin, Chen Sizhong, Wu Zhicheng Model of excitation of random road profile in time domain for a vehicle with four wheels. International Conference on Mechatronic Science, Electric Engineering and Computer (MEC), 2011, p. 23322335.

Zhang L., Zhang T. Study on general model of random inputs of the vehicle with four wheels correlated in time domain. Transactions of the Chinese Society of Agricultural Machinery, Vol. 12, 2005, p. 2931.

Gillespie T. D. Fundamentals of Vehicle Dynamics. SAE International, Warrendale, PA, SAE Technical Paper R114, 1992.

Lu Y., Yang S., Li S. Heavy vehicle dynamics with balanced suspension based on enveloping tire model. Frontiers of Mechanical Engineering in China, Vol. 5, Issue 4, 2010, p. 476482.

Mullari T., Kotta Ü. Simplification of the generalized state equations. Kybernetika, Vol. 42, Issue 5, 2006, p. 617628.

Li D., Lv J., Ying J. The analysis of static and dynamic characteristic for SX360 type vehicle frame. Mechanical Science and Technology, Vol. 21, Issue 3, 2002, p. 380386.
About this article
The support of the Hightech Industrialization Projects of Shanxi Province, China (No. 20112368) and The Natural Science Foundation of Shanxi Province, China (No. 20120110242), are gratefully acknowledged. Thanks are due to Dr. Hongbin Ren and Dr. Jinchao Zhao, of Beijing Institute of Technology, for their help in road roughness modeling.