Abstract
This article introduces the basic theory of modal analysis. The modal testing system is established for a developing truck body. Using a pseudo-random excitation signal along the x, y, and z directions, the modal test is carried out to obtain the dynamic performance by multi-point excitations. The transfer function set is obtained by averaging all transfer functions, and the set is processed by order selection to get 28 order modes. Fitting modal parameters and normalizing modal mass can gain natural frequencies and vibration modes of cab. Then, the modal analysis of truck cab is calculated based on the finite element method, and the results are compared with the test results. The improved measures are put forward to enhance the local stiffness, avoid modal coupling, and reduce vibration noise. So, this article supplies a reference for the dynamic test for the large body.
Introduction
With the rapid growth of buildings, long-range transportation, and logistics, the heavy lorry has a huge development in the automotive industry 1 because of the advantage of dead-weight. In the running process, extrinsic and intrinsic motivations, induced by the uneven surface of road, the change of speed and direction, the vibration coming from wheel, engine and transmission system, the impact of gear, and so on easily cause a strong vibration of cab. The vibration brings many disadvantages, such as making the driver and passengers feel uncomfortable, bringing noise and fatigue failure of components, damaging the protective layer and sealing capacity, and weakening the corrosion resistance of the cab. 2 Therefore, the modal analysis of cab body-in-white (BIW) is important. Especially, the low-order elastic modal reflects the stiffness performance of the whole vehicle body and is the key parameter to control the vibration of vehicle.3,4
There are also many studies on the technique of modal parameter identification. Experimental modal analysis (EMA) or modal testing is an important engineering technique to determine the structural modal parameters. And EMA can be combined with the computer-aided engineering (CAE) that provides theoretical simulation to verify the structural theoretical model, and hence, the virtual testing can be implemented to assist product design and development.5–7 In the late 1960s, Cole and Henry 8 proposed a random decrement technique (RDT) based on the single-point excitation and single-order modal testing and applied it to the modal parameter identification of space shuttle. The random decrement method is a single-channel signal, which cannot be directly used for modal parameter identification. 9 Using the response signal in time domain, Ibrahim 10 put this technology into the field of multi-channel signals. With the development of the modal analysis technology, Ibrahim time domain (ITD) identification was put forward. This method can identify the modal parameters only using response without inputs; therefore, it can be only applied to Gaussian distribution whose mean value of response signal is zero. And it needs many tests, so it is easy to be interfered by noise and has high signal-to-noise ratio. 11 Box et al. 12 proposed a time-series analysis method used in modal analysis, and this approach did not produce energy leak, had high resolution, and can be used for real-time online modal analysis, but it can only identify local modal and cannot describe the overall modes under white noise excitation.13,14
Brown et al. 15 put forward a least-squares complex exponent (LSCE) method to identify the modal parameters of multiple responses, which is a single–reference point complex exponential (SRCE) method. Yu 16 adopted the single-input multi-output analysis method to obtain the modal parameters of white body, compared with the finite element simulation results, in order to optimize the design of the body. Because of uneven distribution of energy in the system, the method is applied to the local mode identification. Van der Auweraer et al. 17 and Verboven 18 developed a least-squares complex frequency domain (LSCF) method to identify modal parameters in frequency domain. This method is more widely used, but there are many defects, such as power consumption, truncation error, frequency aliasing, and off-line analysis. In time domain, the modal parameter identification can be carried out only using response signal; therefore, it can make the real-time online modal analysis for cars and other structures to reflect the actual dynamic performance.
Based on SRCE method, Leuridan and Vold 19 and Vold et al. 20 presented a poly–reference point complex exponential (PRCE) method, which extracted the natural frequencies according to incentive and response signals. Compared with the SRCE, it expanded the amount of data information, and the energy distribution was more even in the system. The real and complex modes can be identified, and the extracted modal parameters are more complete. In addition, this method has the strong ability to recognize the intensive and multiple root mode, and the identification accuracy is greatly improved. Ma et al. 21 used PRCE to get modal parameter identification of a car BIW, and the test results agreed with the simulations.
Due to a large number of truck under long-distance freight work, the drivers have to drive in a long time with heavy load, so the high request is put forward for the heavy truck comfort.22,23 The multi-point excitations are widely adopted in order to get more accurate modal parameters in modal test for vehicle, but the most modal tests are based on the assumption of stationary white noise excitation. In this article, considering the pseudo-random signal as the incentive in order to simulate the actual environment accurately, the dynamics parameters of a heavy truck are solved using the least squares method as parameter identification and adopting multi-point pulsing and multiple points receiving approach. Then, using the finite element method (FEM) to solve kinetic parameters under modal coordinates,24,25 this article describes the contrast analysis between the simulation results and the test to check the validity of the experiment. So, this article supplies a reference for the dynamic test for the large body.
Theory of modal analysis
The modal test theory is to test the system by the dynamic measuring technology.26,27 By processing input and output signals, the testing transfer function is identified to get the inherent characteristics using modal parameters. The motion equation can be expressed as
where [M] is the mass matrix, [C] is the damping matrix, [K] is the stiffness matrix, x is the displacement, f is the exciting force, and t is the time.
Equation (2) was derived from equation (1) based on Laplace transform (s is a variable), considering the initial displacement and velocity as zero
or
where [Z(s)] is the generalized impedance of system.
The modal analysis by FEM is based on the “Reticulation” material characteristics to get the mass matrix [M] and stiffness matrix [K]. For the linear or weak damping system, their damping matrix [C] is considered as proportional damping:
Transfer function matrix
So, equation (6) is expressed as
Equation (6) is converted to
EMA is the inverse process of equation (4). To the known input and output signals, the mass matrix, damping matrix, and stiffness matrix can be solved using formula (6).25,28
Modal test of truck cab
Theory of modal test
Using LSCE method,17,18 the modal parameters are identified based on equation (3). In modal test, LSCE is a time domain poly-reference method, which can identify the complete modes by residue. Formula (3) can be expanded in terms of residue at pole using inverse Laplace transformation
After sampling,
where n is the number of sample points, k is the degree of freedom system, and N is the residue needing to fit (modal number). The regression equation is defined as
Based on regression equation (10), formula (9) can be deduced to
All sample points can constitute Hankel matrix, 29 which can be written in the form of row vector
where the first subscript of h is the sample point, from 1 to
The overdetermined equation (12) is solved by using the least squares method to obtain all the residues. Then, the system poles can be obtained based on the above residues and equation (10). Modal shape, mass, stiffness, damping ratio, and so on are identified according to the relationship among the modal parameters.
Establishment of the test system
Composition of the test system
Experimental system consists of three parts: excitation system, response pick-up system, and modal process analysis systems.30–32 (1) Excitation system: Quad Digital-to-Analog Converter (QDAC) signal generation module, the power amplifier, and exciter of LMS SCADAS III SC316W; (2) response pick-up system: integrated circuit piezoelectric (ICP) acceleration sensor, ICP force sensor, signal amplification of LMS SCADAS III SC316W, and smart data acquisition system; (3) modal process analysis system: LMS Test.Lab. 33 The test system is described in Figure 1.

The test system.
Experiment scheme
In the modal test, the truck cab was in a free state and suspended with rubber rope. The truck body was hung on the front and rear ends of the left and right longerons using a rubber rope, and it was in four-point suspension status. The cab body was kept balanced as the front hitch point was located in the installation position of bumper, and the rear hitch point was situated in the holes of cab stringer. The lower end of the rubber rope was connected to the cab by a hook, and the upper end was connected to honing car through a pulley.34,35 The suspension is shown in Figure 2.

The hitch positions of truck cab in modal test.
The appropriate excitation signal can provide enough energy and excite all the modes in test, so pseudo-random excitation signal was chosen under the frequency range of 0–200 Hz. Multi-point excitations and multi-point responses (x, y, and z three-way responses) can prevent the loss of important information and avoid the incomplete mode.
In experiment, the excitation points were measured by triaxial sensor, and the pick-up points were tested in seven batches. Each batch included 30 moving sensors to test along the three, x, y, and z, directions. The method of collecting signal increased the experiment workload but improved the test accuracy. The front end of longeron was defined as X exciting direction, which was shown in Figure 3. The front end of right threshold beam was considered as Y exciting direction, shown in Figure 4. The rear end of longeron was Z exciting direction, shown in Figure 5.

The exciting point in X direction.

The exciting point in Y direction.

The exciting point in Z direction.
A total of 213 collection points were identified as shown in Figure 6, and their distribution was as follows: 40 points were located in firewall, 64 points in side plates, 43 points in floorboard, 29 points in roof panel, and 37 points in rear quarter panel.

The layout of test points.
Collection of experiment data
During collection, the transfer and coherence functions of test data were simultaneously processed. The signal was valid when coherent coefficient was above 0.8.36,37 Collecting signals while processing them can improve the precision and avoid repetitive trials. Coherence in x, y, and z directions was shown in Figures 7–9, respectively. The coherent coefficient of few frequencies was below 0.8, and the rest was close to 1. Because of the good correlation between measuring points and excitation points, the signal-to-noise ratio was high and excitation was effective.

Coherence in x direction.

Coherence in y direction.

Coherence in z direction.
The response acceleration was tested and analyzed in three, x, y, and z, directions. All transfer functions were averaged to obtain the transfer function set, and the set was processed by order selection to get 28 order modes, is shown in Figure 10.

The transfer function set.
Data processing
Response frequency and vibration modes
Based on the EMA theory, using PRCE method, the process of modal parameter identification is shown in Figure 11. During modal trials, error was inevitable. Exciting force should inspire each mode, but too high exciting force would cause local nonlinear. Due to the large volume of truck cab, the high frequencies above 100 Hz can be detected rarely. The average of all the transfer functions was computed to get the transfer function set. And the set was processed by the determination of the order to identify the number of modes. The result was shown in Figure 10 with adapting the curve of figure in order to show the peak clearly. Finally, 28 step modes were collected. Using the method of complex modal, single degree to fit the natural modal, and adopting mass normalization method 38 (modal mass is defined as 1) to obtain the vibration mode, all results were described in Table 1 and Figure 12.

Identification flow chart for modal parameter.
The first eight natural frequencies.

The first eight mode shapes: (a) the first torsional mode along X axis (F = 16.499 Hz), (b) swing and torsion mode along X axis (F = 17.483 Hz), (c) local mode in rear side plate (F = 25.030 Hz), (d) local mode in rear roof panel (F = 28.076 Hz), (e) torsional mode along Y axis and local modes in rear roof panel and the left side plate (F = 29.106 Hz), (f) torsional mode along the Y axis and local mode in roof panel (F = 30.773 Hz), (g) first-order bending vibration around Z axis (F = 30.773 Hz), and (h) mode of rear end of floorboard (F = 35.680 Hz).
In the modal test, there were many local modes described in Table 2.
Local modes of cab.
Establishment of the finite element model
Building the finite element model of the truck cab to do modal simulation, the simulating results are compared with the test. In the process, we should pay attention to the following points:39,40
To save time and cost of analysis, the model should be simplified with deleting the fillets less than 10 mm in diameter, small abdicating steps, and transition fillets.
The truck body is made of sheet metal by stamping and spot welding, and its thickness is relatively thin, so shell element in size of 20–30 mm is chosen to discrete the cab. In meshing process, we can control the element quality by aspect ratio, warpage, skew, and Jacobian. The finite model consists of quadrilateral elements in 90%.
The truck cab is assembled together by welding, bolts, riveting, and bonding, using the node coupling method to simulate welding spot, beam elements to simulate the bolts, rigid elements to simulate welding seam, and the spring elements to simulate the bonding.
The cab material is low-carbon steel, the elasticity modulus is 2.1 × 105 MPa, the ratio of Poisson is 0.3, and density is 7.8 × 103 Kg/m3, adopting isotropic material in simulation.
The finite element model of truck cab shown in Figure 13 includes 211,026 elements and 223,851 nodes. There is 203,212 shell elements, which contain 192,345 quadrilateral elements and 10,867 triangular elements. There are 2152 connecting elements, among them 1389 are rigid elements, 413 are beam elements, and 123 are spring elements.

Finite element model of truck cab.
Comparison analysis between simulation and test results
Results contrast
The simulation and test results were analyzed, as shown in Table 3.
Comparison between simulation and test modal results.
FEM: finite element method
From Table 3, the relative error between simulation and experiment results is less than 5%, so the simulation modes are consistent with the test measured. Therefore, the results of simulation and test are reliable.
Improvements of cab
The low-order modes less than 100 Hz of the truck cab are gained in test and simulation. Combined with the actual condition, some suggestions for improvement are proposed as follows.
Figure 14 shows the local mode in left center of rear side along x direction. We suggest that position 1 pointed by arrow in rear part should be stamped into position 2 pointed by arrow, or glued, to prevent vibration and noise.

The rear wall of cab.
From Figure 15, the vertical beams in both sides of roof also exist the local mode, so the vertical beams and the roof sides should be bonded together using foam.

Side panels of cab.
Owing to substandard quality of welding or lack of local stiffness, the local mode exists also in the right rear part of the floorboard (location 2 pointed by arrow in Figure 16). Position 1 in left front floor has an elliptical hole compared with the right symmetrical position, and it will result in insufficient stiffness and local modal, so the floorboard around hole should be stamped into flute type.

Floorboard of cab.
According to Figure 17, the local mode is intensive in rear roof and roof beam, and its frequency is relatively low, so it is easy to generate modal coupling with passenger compartment and cause the resonate noise. The roof pointed by arrow should be stamped into wavy stripes along X axis or glue. In addition, the local mode is relatively obvious in the rear of floor, and it may be related to the large opening at the bottom, hence more attention must be paid to improvements.

Roof of cab.
Conclusion
The dynamic characteristics of a truck cab are identified using experiments and simulation methods in this article. Modal test scheme is constructed, and the finite element model of the cab is established. Below 100 Hz, 28 step modes are extracted and the test and FEM simulation results are consistent according to comparison. The following observation can be made from the study:
The natural frequencies and modal shapes obtained from experiment and simulation are almost close to each other, so the modal test plan and simulation model are reliable, and this article provides the important design basis to improve truck cab.
The first frequency of truck body should not be equal or close to the external excitation frequencies. The external excitations mainly come from road, engine, transmission, and wheel, and among these incentives, the effect of engine is greatest. Based on the formula of excitation frequency, 41 the engine idling rotational frequency of truck is about 12.5 Hz. To avoid resonance, the first modal frequency should be different from the excitation frequency at least 2 Hz, so it should be larger than 14.5 Hz. As shown in Table 1, the first-order frequency is 16.5 Hz, so the truck will not resonate and the truck body structure is reasonable.
The first torsional mode and the first bending mode have the greatest impact on the truck body; therefore, in order to avoid coupling, the difference between them must be more than 3 Hz. From Table 1, the first torsional vibration frequency is 16.5 Hz, and the first bending frequency is 32.9 Hz, without coupling.
In the process of test and simulation, there are many local modes, and they are most likely caused by unreasonable local structure, low-quality welding, or insufficient local stiffness, such as in left center of rear side along x direction, in the vertical beams, on both sides of roof, and on the right rear part of the floorboard. Those local modes tend to generate modal coupling with passenger compartment and cause the resonate noise, so those positions should be paid more attention in the future improvements.
Footnotes
Acknowledgements
The authors owe a special thank you to Song Shaofang for her suggestions and insights on test.
Academic Editor: Sandra Velarde-Suáre
Declaration of conflicting interests
The authors declare that there is no conflict of interest.
Funding
This work was supported by the Project of National Natural Science Foundation of China (Grant no. 51175320), and partly supported by the Program for Professor of Special Appointment (Eastern Scholar) at Shanghai Institutions of Higher Learning and the Fund for Talents Development by the Shanghai Municipality, China.
