Nonlinear vibration of semitrailer suspension: wavelet and multiscale entropy-based approaches

We examine the dynamical response of a semitrailer suspension to a real road proﬁle under two different loads. The trailer is a low loader for transporting heavy and large loads, the conﬁguration of which can be adjusted as desired. Experiments are performed on both unladen and 57-ton laden trailers. The acceleration of unsprung and sprung parts of the vehicle is measured. The obtained acceleration time series are then analyzed using the Fourier and wavelet transforms as well as the multiscale entropy approach. The data analysis results reveal transmission of suspension vibration and complexity of vehicle vertical acceleration as a response to the road proﬁle.


Introduction
The effect of a rough road surface on vehicle vibrations is still researched by automotive manufacturers and research groups [1][2][3][4][5][6][7]. These studies usually focus on the drivers' and passengers' comfort during their travel. In the case of goods vehicles, conditions of load carrying are also sometimes taken into account [8,9]. For a better understanding of the problem, the behavior of different types of vehicle suspensions was tested under real road conditions. These tests were performed in order to identify and optimize the suspension parameters [10][11][12], and to minimize tire wear [13]. The relationship between anti-vibration performance and characteristics of a vehicle suspension has been analyzed intensively in recent years [16][17][18]. Based on real experimental conditions, several models have been proposed. These models take into account such variables as vehicle kinematics, road curvature, and internal vehicle excitations [12,14,15]. Despite the visible progresses in simplified vehicle modeling, semitrailer suspensions, especially heavier ones with a modular design, have not been studied extensively. The investigation of such vehicles is particularly important not only because of the structure overload and tire wear but also because of dynamic interaction between the wheels and road surface that causes surface degradation and exerts great impact on road infrastructure such as bridges and railway crossing. What is more, the total mass of heavy semitrailers can vary to a greater extent than it is the case with passenger cars or lorries; a fully laden semitrailer can be four or even more times heavier than an unladen one. The aim of this study was to analyze the vibration of a heavy semitrailer during its operation under usual road conditions. The results will give additional insight into performance of suspension of heavy modular vehicles. In addition, they can be used to estimate transmitted loads acting on the main supporting structure of a semitrailer. The knowledge of the dynamic overloads is vital to the design process because the assumed overload factor influences the accuracy of strength calculations of a new design [19,20].
The nonlinear dynamics of trailer suspension plays a main role in filtering the rough road excitation in a wide range of loading and in the presence of the coupling force from the tractor. Our study compares road excitations of the driving unsprung elements and the sprung platform in certain (natural) conditions concerning the two extreme cases: semitrailers unloaded and fully loaded. Both cases are important for producers and users of semitrailers due to transport safety, costs of operation, and maintenance of cargo transported in good condition. In the case of empty trailer (light transport), suspension system is more stiff, resulting in increased transmission of vibration into the platform. With a maximum load of energy, transmission can be caused by a failure of the suspension system with specific excitations that might occur under working conditions. The proposed studies would allow to evaluate the efficiency of the suspension system in both extreme cases, when operating in seminatural conditions.

Object and conditions of measurements
The examinations were conducted on a prototype of a modular semitrailer for carrying oversized and heavy loads up to 60 tons. Depending on the requirements, the length of the semitrailer platform can be extended, and the number of axles can be set to 2, 3, or 5. The number of axles can be changed owing to the modular structure of the trailer. The semitrailer can be used as a 2-or 3-axle bogie vehicle, or both. All of the axles are equipped with steered twin wheels and tires 245/70 R 17.5. The nominal load of the axles is 12 tons. The maximum construction weight of the semitrailer is 84 tons. The semitrailer is equipped with a typical air suspension consisting of a control arm, an air bag, and a hydraulic shock absorber. A schematic design of the semitrailer used in the experiments is shown in Fig. 1.
The vibration data were measured using PCB Piezotronic manufactured equipment. Single-axis accelerometers type ICP M338A34 with an average sensitivity of 1 mV/(m/s 2 ) were used. Their measurement range was +/−4900 m/s 2 , and their frequency range was from 0.7 to 3000 Hz. The acceleration signals were conditioned with the use of a four-channel amplifier type 482A16. A DAQPad 6070E from National Instruments was used for data acquisition. The maximum sampling frequency was set to 1.25 MHz, and resolution of analog-todigital converter was 12 bits. In order to utilize the whole scale of the converter, the voltage amplification of the signals was set to 100, and the measurement range was varied according to maximum signal amplitude. During the research, four vibration transducers mounted on the front axle of the trailer suspension were used. Two transducers were mounted at the flat spring rocker arms close to the right and left wheel axes, and the other two were mounted on the body above the right and left rocker arms mounting pins. Such configuration of the measurement points enabled measurement of both sprung (body) and unsprung (wheel axles) acceleration (Fig. 2).
The measurements were taken on a 5-axle semitrailer in two cases: (a) unladen semitrailer; (b) semitrailer loaded with 57 tons. During the experiments, the semitrailer was hauled by a 3-axle semitrailer truck at a speed of 9 km/h through the same section of the ground paved with concrete blocks. The vibration data were recorded simultaneously by four measurement transducers with sampling rate of 20 kHz per channel. At each measurement event, up to 200,000 samples were recorded.

Experimental results
The collected data for the unladen semitrailer are given in Fig. 3a-d and for the laden one in Fig. 4a-d. Note that the unsprung mass vibrations generally have a higher amplitude and, additionally, are characterized by   Table 1). (Color figure online) sudden increases in their value. This was caused by the external event surface irregularity. The results of vibration measurements of the semitrailer moving at a speed of 9 km/h on the surface made of cobblestones were analyzed. This interesting events could be investigated in a frequency domain. Consequently, the Fourier transform can be employed to provide the system response features including frequency of transmission of vibration through the suspension system. The results of the FFT were limited to 0-60 Hz band. The results for the unladen semitrailer are presented in Fig. 5a, b and for the laden one in Fig. 5c, d. Interestingly, we observed two peaks (for the frequency f of about 3 and 8 Hz) in the left-and right-hand side sensors of the unladen semitrailer. In laden trailer case, the peaks observed in the frequency spectrum of sprung masses decrease to about 1.5 and 4 Hz as an effect of a fairly larger load participating in sprung masses. The resulting signal transmittance Tr( f ) can be expressed as: where Acc u ( f ) and Acc s ( f ) are the unsprung and sprung accelerations, respectively. Note that in the most of the frequency region, the vibrations are damped effectively, and the resulting transmittance is less than unity. The damping of vibration is less effective in the range of 2-7 Hz and, especially in the case of unladen trailer, in the ranges of 40-45 and 55-58 Hz.

Wavelet analysis
For better clarity, we have also performed the wavelet analysis. Namely, we studied the time series Acc i (t) (i = s or u) given by the continuous wavelet transform (CWT) [21,22]. The corresponding continuous wavelet transform with respect to a wavelet function ψ(.) is defined as follows: where < Acc u/s > and σ Acc(u/s) are the averages and standard deviations for unsprung and sprung accelerations, respectively, and t i is the sampling time. Finally, the wavelet ψ(t) is referred to as the mother wavelet, and the letters s and n denote the scale and the time index, respectively. The wavelet power spectrum (WPS) of the Acc u/s time series is defined as the square modulus of the CWT: In our calculations, we have used a complex Morlet wavelet as the mother wavelet. The Morlet wavelet consists of a plane wave modulated by a Gaussian function and is described by where θ 0 is the center frequency, also referred to as the order of the wavelet, η is a renormalized time variable. The value of θ 0 defines the number of oscillations in the wavelet and thus controls the time-frequency resolutions. In our analysis, we used θ 0 = 6 This choice provides a good balance between the time and frequency resolutions. Also, for this choice, the scale is approximately equal to the period, and therefore, the terms scale and period can be interchanged for interpreting the results. For general discussion, see [21][22][23].
The corresponding results of CWT as the wavelet power spectra are presented in Fig. 6. One can see two characteristic periods 0.12 and 0.32 s of the unladen vehicle; 0.32 and 0.64 s of the laden vehicle. This is consistent with the corresponding frequency peaks in Fig. 5. Moreover, by using CWT, one can see the evolution of the response amplitude and frequency in time domain. The colors also help us to distinguish time-localized events (Fig. 6c) as well as the intermittent behavior mainly for the unsprung masses. Generally, the spectra in Fig. 6a, c are more uniform due to random wheel excitation by road surface corrugation.

Composite multiscale entropy approach
To gain more insight into complex systems and their nonlinear behavior, the multiscale entropy [24,25] analysis is widely used [7,26,27]. This method provides, for measured signals, a relative level of complexity of finite length time series. Unfortunately, there is no consensus over the complexity definition, but it is combined with "meaningful structural richness" [28] contained over multiple spatial and temporal correlations. The concept of multiscale entropy (MSE) [24,25] is based on the coarse-graining procedure that uses a coarse-grained time series, as an average of the original data points within nonoverlapping windows by increasing the scale factor τ according to the following formula (see Fig. 7): where A i is a raw one-dimensional acceleration time series x = {A 1 , A 2 , . . . , A N }. In this approach for each scale factor τ , the MSE calculation based on the time series of the coarse-grained A (τ ) j : where m = 2 is the pattern length, and r is the similarity criterion which is usually set to be r < σ [24]; here, σ is the standard deviation of the original time series and A (1) i = A i . To estimate SampEn (A (τ ) , m, r ) from Eq. 6 (see also Fig. 7), we count the number of vector pairs denoted by A (τ ) (i) and A (τ ) ( j) in the time series of length m and m + 1 having the distance d[A (τ ) (i), A (τ ) ( j)] < r . We denote them by P m and P m+1 , respectively. Finally, we define the sample entropy to be [29]: Strictly speaking, it is a shortcoming of the logarithm of the conditional probability that two sequences with a tolerance r form points that remain within r of each other at the next point. Continuing these research efforts, Wu et al. [26] introduced the concept of a composite multiscale entropy (CMSE), which for higher scale factor provides entropy more reliably than the usual multiscale entropy by including multiple combinations of neighbor points. The prescribed algorithm for CMSE calculations is the following formula: In the present paper, the CMSE is applied for real signals reporting a different behavior of the system characterized by the time series in Figs. 3 and 4. The voltage output results of composite multiscale entropy calculation are usually provided for several similarities r for different CMSE values with increasing the scale factor τ . The appearance of higher values of CMSE corresponds to the presence of higher complexity within the analyzed signal. It can be noticed that there is a variation in the irregularity of the system as a different τ is chosen.
The CMSE results, for the chain length of m = 2 and similarity criterion r = 0.2σ , are presented in Fig. 8. The CMSE for laden semitrailer at the unsprung point (green line) is higher than for the sprung point (blue line) as shown in Fig. 8a, b. This can be explained by the fact that the unsprung vehicle part undergoes significant in the vibration signal, which correspond to the signal carrying some information about the ground noise. Therefore, the entropy has a relatively low value. At the sprung part, these changes are selectively suppressed and complexity increases. Interestingly, the sprung response is close to the pink noise (more flat curves), while the unsprung response indicates rather the white noise (a hyperbola-like curves).
On the other hand, with the unladen semitrailer, the difference between the unsprung and sprung response entropy is fairly small as a result of larger vibration transmission (Fig. 8c, d). Both of responses (sprung and unsprung) indicate the white noise type of excitation (a hyperbolic shape of curves). Simply, the semitrailer suspension is designed to damp more effectively vibrations in the laden case.

Conclusions
We have analyzed the effect of load on vibration transmission through the heavy semitrailer suspension. The suspension operates satisfactorily, showing in a wide region of frequencies a fairly low value of vibration transmission Tr( f ) < 1. However, acceleration amplitudes of unsprung and sprung masses in a low-frequency range (up to 7 Hz) both for the laden and unladen vehicle are similar, which suggests that the suspension does not provide proper damping. In the range of higher frequencies, damping effectiveness increases both for the unladen and laden vehicles. Comparing the absolute values of the amplitudes for laden and unladen semitrailers, we found that they are higher for a relatively light vehicle. Additionally, the higher load provided more effective damping. We also observed a reflection of the resonance peaks of the sprung mass in the Fourier and Wavelet spectra where the transmission increased considerably. In particular, the wavelets enabled to investigate the transient and nonstationary behavior and the corresponding sprung mass responses including the time localizations of vibrations and intermittent behavior. On the other hand, the multiscaled entropy revealed a huge difference with respect to the sprung and unsprung parts in the laden vehicle case. At the same time, however, this difference does not hold for the unladen semitrailer case, which confirms that the semitrailer suspension works efficiently in the laden trailer case. To provide more information about the suspension system, it is necessary to perform a series of new tests that will take into account different road profiles and surfaces as well as various loading and vehicle speeds.