Vibration Information Acquisition and Representation Control Strategy for Construction Vehicle Driver Seat

Vibration information acquisition and representation control strategy in the experiment of construction of vehicle driver seat are investigated in this study. An improved adaptive algorithm of multisensor attitude data fusion is proposed. The nonlinear equations of vehicle attitude quaternions are constructed by using the test information of accelerometers and magnetometers. Then, a complementary ﬁlter is used to fuse the obtained attitude with gyroscope information to obtain attitude information during vehicle travel. Weight matrix and ﬁlter coeﬃcients in the algorithm are adjusted with vehicle motions to ensure accuracy of attitude measurement under static, low-, and high-acceleration conditions. The vehicle vibration displacement information is reconstructed by the low-frequency attenuation integration method on the basis of attitude measurement to inhibit integral trend items and reduce vibration displacement reconstruction errors. The three-state controller is designed to control a six-degree-of-freedom hydraulic simulated vibration test bed, and the obtained vibration signal is used as an excitation signal to perform a vibration representation control test. Finally, results verify that the proposed method demonstrates favourable data fusion and representation accuracy.


Introduction
When construction vehicles drive on complex roads, various parts of a driver seat are constantly affected by vibration and impact of uneven road, which may be a huge risk for vehicle safety and cause damage in various degrees [1][2][3]. Accelerating the development of vehicle seat products and ensuring the reliability of the product are effective approaches to avoid these damages by using a real vehicle to collect vehicle vibration information of different road conditions and then adopt a vibration test system to simulate a complex vibration environment [4][5][6][7]. At present, a single degree-offreedom (DOF) seat vibration fatigue test bed has a singlevibration form, making it difficult to simulate complex vibration environments [8,9]. Although the multi-DOF vibration test bed can simulate a complex vibration environment, researches on the multiaxis vibration test bed mostly focus on the vibration control algorithm [10][11][12][13][14][15][16][17]. Considerable researches on road excitation information acquisition have been conducted. However, most of these studies use road roughness information generated by a power spectral density function of a standard road level [18,19] or a rigid road [20,21], and studies on obtaining actual driving vibration information are limited. e basis of vibration environment simulation is to acquire vibration information during vehicle driving. Vibration information includes vehicle attitude and vibration displacement [22][23][24]. e calculating algorithm for an attitude angle is the core to obtain vehicle attitude, an important factor that affects the record accuracy of vibration information. e design and use of a reasonable attitudecalculating algorithm are the most important issues in an attitude measurement system [25]. e vehicle attitude measurement system in this study is composed of gyroscope, accelerometer, and magnetometer. ese components are suitable for various complex environments [26]. In accordance with the respective characteristics of sensors, the collected data are fused to enable a complementary sensor performance and obtain highprecision attitude measurement data. Researchers propose two kinds of filtering algorithms for the multisensor information fusion problem. One is Kalman filtering and various algorithms, such as Kalman filter (KF), extended KF, unscented KF, and federated KF, are derived on the basis of Kalman filtering [27][28][29]. e other is complementary filtering (CF), and various improved filtering algorithms, such as classical CF, linear CF, and CF based on the gradient descent method, are derived on the basis of CF [30,31]. e vehicle vibration displacement information is difficult to be detected directly, and an indirect measurement method is used to obtain the actual test. At present, the indirect measurement method of the accelerometer is extensively used [32,33]. e signals of velocity and displacement are obtained by single and double integrals of acceleration signals, respectively. e use of the acceleration signal to reconstruct velocity and displacement information is an effective method for obtaining a displacement signal.
ree main methods are previously mentioned. e first method considers the integral of the acceleration signal in a time domain.
is method is simple and has a minimal amount of calculation. It has been applied to the field of aircraft line planning and navigation [34]. However, this method has limited processing capability for trend terms, and the obtained displacement signal has serious trend-term errors occasionally. e second method is decomposing a signal in accordance with its characteristics. e empirical mode decomposition method is used to correct the displacement resulted from an integral of the acceleration signal [35]. e third method is using Fourier transform to convert acceleration, velocity, and displacement [36,37].
is method is effective in handling trend items; however, it has numerous calculations and is unsuitable for real-time data processing.
To improve the control accuracy of the multisensor data fusion, the adaptive control algorithms are developed to guarantee the boundedness of adaptive estimation [38,39] and have integrated adaptive control to achieve the online identification of component parameters of the PV system. Moreover, fuzzy control is used to approximate the unmodeled part of system. In [40], the fault is accommodated by reconfiguring the formation structure through the novel notion of the weighted absolute measurement formation digraph and then a robust controller for the partially low-level recovered vehicle is activated.
is study aims to increase accuracy of vehicle attitude measurement at the static, low-, and high-acceleration states and propose an improved adaptive algorithm of multisensor attitude fusion. Vehicle attitude quaternion nonlinear equations by using the test information of the accelerometer and magnetometer are developed in this paper. en, the complementary filter is used to fuse the attitude and gyroscope information to obtain the attitude information during vehicle driving. Moreover, according to the weight matrix and filtering coefficient of the vehicle motion state adjustment algorithm, the vehicle attitude measurement under the static, low-, and high-acceleration states achieves good results. Finally, a six-degree-of-freedom hydraulic simulation vibration test bed is designed, and the collected vibration signals are used as the excitation signals to carry out the vibration reproduction control experiment.

Vibration Representation Test Bed.
e vibration test bed is based on the Stewart structure and consists of an upper platform, a lower platform, and six branched chains. Each chain is composed of a hydraulic cylinder, an upper base, a lower base, and a universal joint. e top of the upper base is fixed to the lower surface of the upper platform. e lower end is fixed to the upper end of the universal joint, and the lower end of the lower base is fixed to the upper surface of the lower platform. e driver seat of a vehicle is fixed to the upper end of the upper platform. e upper platform has six DOFs, and its structure is illustrated in Figure 1. Table 1 displays the technical index of the vibration test bed.
In Table 1, α denotes the roll angle around the x-axis, β denotes the pitch angle around the y-axis, c denotes the yaw angle around the z-axis, a denotes a displacement range along the x-direction, b denotes a displacement range along the y-direction, and c denotes a displacement range along the z-direction.

Electrohydraulic Servo Control System.
A control system of a hydraulic cylinder comprises a computer processing unit, a controller, an A/D converter, a D/A converter, an electrohydraulic servo valve amplifier, and a displacement sensor. e collected acceleration signal is transmitted to the computer processing unit after being processed using a filter. A digital signal that controls the expected displacement of the hydraulic cylinder is compared with the actual displacement signal of the hydraulic cylinder fed back by a vibration representation device. Subsequently, the difference is transmitted to the central processing unit to be calculated by the control algorithm and then outputs a control result. e current signal in the electromagnetic coil of the electrohydraulic servo valve is converted by using a D/A converter to control the displacement of the electrohydraulic servo valve spool and the extension of the hydraulic cylinder and drives the motion representation of the 6-DOF test bed. e displacement sensor in the extended end of the hydraulic cylinder detects the displacement information and transmits this information to the A/D converter to convert it into a digital signal. e displacement signal is compared with the displacement signal calculated by using the obtained data to determine the difference which forms a closed-loop control of the hydraulic cylinder displacement. Figure 2 depicts the structure of the hydraulic servo control system. In Figure 2, u i (i � 1-6) is the control quantity which is the output of the D/A converter to the servo amplifier and y i (i � 1-6) is the piston displacement of the hydraulic cylinder. Moreover, x is the displacement of the servo spool; p b and q b are pressure and flow of the rod cavity of the hydraulic cylinder, respectively; p a and q a correspond to pressure and flow of a rodless cavity of the hydraulic cylinder; and p s is the working pressure of the system.

Attitude Measurement Algorithm.
A gyroscope has a remarkable dynamic performance, but an accumulative error during the calculating process results in a time lapse. Accelerometer and magnetometer use gravity and geomagnetic elds for attitude measurement, respectively. No accumulative error has been observed. However, the dynamic performance is poor, and the measurement accuracy is a ected by motion acceleration and electromagnetic interference. Di erent sensor characteristics are used in the frequency domain to obtain accurate attitude information, and the CF scheme is used to fuse the measurement data obtained by using the gyroscope, accelerometer, and magnetometer. Figure 3 demonstrates the principle of the fusion algorithm. In Figure 3   is the parameter of the complementary lter, and W is the weight matrix. Firstly, the vehicle motion state is judged by the accelerometer-measured values to determine the weight matrix W, and the interference of vehicle motion acceleration is limited by adjusting the weight matrix. Secondly, nonlinear attitude quaternion equations are constructed by using the accelerometer and magnetometer information and solved by using the Gauss-Newton method to fuse the accelerometer-and magnetometer-measured values.
irdly, lter parameter β f is determined by the vehicle motion state to achieve the adaptive cut-o frequency of the complementary lter. Finally, the measurement data obtained by using the gyroscope, accelerometer, and magnetometer are fused to ensure the accuracy of attitude measurement under static and dynamic conditions.
(1) Attitude Update Algorithm of the Gyroscope. Runge-Kutta method is used to solve quaternion di erential equations numerically in accordance with the measured angular velocity of a body. e discrete quaternion update equation is

Motion base
Hydraulic cylinder

Electrohydraulic servo system
Servo valve Platform Controller Figure 2: Principle of a closed-loop control of the hydraulic servo control system.  Shock and Vibration where Q g denotes the output attitude of the gyroscope, Q k−1 indicates the last estimated attitude of the complementary filter, ΔT represents the sampling interval, and K 1 , K 2 , K 3 , and K 4 are obtained by the following equation: where Ω k is the output angular velocity matrix of the gyroscope at the current time and Ω k is determined by the following equation: where w b xk , w b yk , and w b zk denote the output values of the three axes of the gyroscope at the current time.
(2) Attitude Information Fusion of the Accelerometer and Magnetometer.
e estimated gravitational acceleration vector in the vehicle coordinate system is e gravitational acceleration vector a b measured by the accelerometer should be equal to the estimated gravitational acceleration vector g b when a vehicle is in a static state. us, equation (5) can be obtained: e geomagnetic vector measured by the magnetometer in the initial attitude of the vehicle is [m x0 , m y0 , m z0 ] T . During the motion of vehicle, the estimated geomagnetic field vector in the vehicle coordinate system can be obtained using the coordinate transformation matrix b g R from the geodetic coordinate system to the vehicle coordinate system: e geomagnetic field vector m b measured by the magnetometer must be equal to the estimated geomagnetic field vector m b . us, equation (7) can be obtained: Normalised quaternion is used to represent an attitude, that is, ‖Q‖ � 1.
us, the quaternion also satisfies the following equation: Let f a (q 0 , q 1 , q 2 , q 3 ) � [f a1 , f a2 , f a3 , f a4 , f a5 , f a6 , f a7 ] T ; the overdetermined nonlinear equations for q 0 , q 1 , q 2 , and q 3 are obtained from equations (5), (7), and (8): Quaternion attitude information can be obtained by solving equation (9) using the accelerometer-and magnetometer-measured values at the current moment. However, the equations disregard the weight of each equation in the data fusion and can obtain accurate results under static and uniform motions. Acceleration interference information is obtained from the accelerometermeasured values when the vehicle has a moving acceleration, and errors are identified in the test result. To limit the influence of the vehicle moving acceleration on data fusion results, this study introduces the weight matrix W to modify equations. W is a diagonal matrix that contains three parameters, namely, λ a , λ m , and λ Q , which is defined as follows: Shock and Vibration 5 where I 3×3 is the unit matrix, λ a denotes the weight of the accelerometer-measured values, λ m denotes the weight of the magnetometer-measured values, and λ Q denotes the weight of the normalised quaternion equation. After considering the weight of each equation, the nonlinear equation for q 0 , q 1 , q 2 , and q 3 is where x � [q 0 , q 1 , q 2 , q 3 ]. e solutions of equations are transformed into a weighted least square problem: e aforementioned equation is solved to obtain the attitude quaternion Q a from the accelerometer-and magnetometer-measured values through the Gauss-Newton iteration method.
where Q is the optimal estimated attitude which is the output of the complementary filter and β f is the parameter of the complementary filter. e relationship between β f and filter cut-off frequency f T is e attitude data are dominated by the gyroscope when the signal frequency is greater than f T , whereas the attitude fusion results of the accelerometer and magnetometer and the low-and high-frequency interferences are separated by the complementary filter effectively.

(4) Adaptive Adjustment Scheme of Filter Parameters and Weight Matrix.
e accelerometer-measured value is the vector sum of the gravitational and moving acceleration when a vehicle has a moving acceleration. At this time, the mode of accelerometer measurement is no longer the magnitude of the gravitational acceleration. erefore, a vehicle motion state can be judged by the modulus a b of the accelerometer-measured value. e error of the upper limit σ 0 is given (σ 0 is a relatively small positive real number) because of the measurement error. e vehicle is considered in a static or a uniform motion state when the magnitude of accelerometermeasured value satisfies the following equation: ����������������� e error is small when a vehicle has low moving acceleration although the accelerometer introduces an attitude measurement error. At this time, the accelerometermeasured value still plays a beneficial role in the attitude fusion and cannot be completely ignored. Owing to a positive real number σ 1 , the vehicle is considered in a low acceleration state when the accelerometer-measured value satisfies the following equation: A vehicle is considered in a high-acceleration state when the accelerometer-measured value satisfies equation (17). At this time, the accelerometer-measured value is irrelevant to attitude measurement: ����������������� e adjustment scheme of the weight matrix and filter parameters is formulated by the motion state of a vehicle. e parameters are listed in Table 2.

Experiment for the Attitude Measurement.
e comparison experiment for attitude measurement is conducted by using the VG440 model gyroscope of Memsic Semiconductor Company ® to verify the effectiveness of the proposed algorithm. Vertical gyroscope and sensing device are fixed on the test carrier plate to ensure that they perform a random rotation motion in each direction. Subsequently, the attitude information obtained from these materials is recorded and compared. Figures 4-6 exhibit the comparison of the measured data and the error of roll, pitch, and yaw angles, correspondingly. In the random rotation, the attitude angle error is less than 1°, thereby satisfying the design requirements.

Vibration Displacement Reconstruction Algorithm.
e implementation of the vibration displacement reconstruction algorithm is as follows: Firstly, the accelerometer is compensated by the gravitational acceleration to obtain the moving acceleration that is indicated in its own coordinate system. Secondly, the moving acceleration of the vehicle coordinate system is transformed into a geodetic coordinate system by using an attitude transformation matrix, and the time history of vehicle acceleration is obtained in the geodetic coordinate system. Finally, lowfrequency attenuation integration algorithm is used to reconstruct the vibration displacement information from acceleration information. e principle of the vehicle vibration displacement reconstruction algorithm is displayed in Figure 7, where ϕ, θ, and ψ denote the roll, pitch, and yaw angles of the vehicle.
(1) Gravitational Acceleration Compensation. Whether a vehicle is still or moving, the output of the accelerometer contains components of the gravitational acceleration on each axis. In fact, the output value of accelerometer is the vector sum of the moving and gravitational accelerations. erefore, gravitational acceleration compensation must be performed on the output value of the accelerometer to obtain the moving acceleration of the vehicle. e specific method calculates as follows: 6 Shock and Vibration where a b m is the moving acceleration of the vehicle, a b is the output value of the accelerometer, b g R is the attitude transformation matrix from the vehicle coordinate system to the geodetic coordinate system, and g 0 0 G T represents the local gravitational acceleration vector in the geodetic coordinate system.
(2) Conversion of Motion Acceleration to the Geodetic Coordinate System. e obtained moving acceleration a b m after gravity compensation is expressed in the vehicle coordinate     Shock and Vibration 7 system and must be converted into the geodetic coordinate system to obtain the moving acceleration information in space. e representation of the moving acceleration a g m of a vehicle in the geodetic coordinate system is as follows: where g b R is the attitude transformation matrix from the vehicle coordinate system to the geodetic coordinate system.
(3) Displacement Is Reconstructed by Acceleration. e threeaxis acceleration is processed by using the low-frequency attenuation algorithm to determine displacement after obtaining the moving acceleration information in the geodetic coordinate system. e integral displacement S g m is obtained as follows: where F and F −1 represent the Fourier transform and inverse Fourier transform, correspondingly. A vibration of the test carrier is implemented in all directions, and the vibration attitude of the test carrier is reconstructed by using this algorithm. Figure 8 presents the 6-DOF motion information.

Design of a Representation System Controller.
e threestate feedbacks are position, velocity, and acceleration feedback. e three-state control includes three-state feedback and feedforward control. e introduction of acceleration feedback can improve the damping ratio of system under the stable condition of the inner loop. e introduction of velocity feedback reduces the position closed-loop damping ratio slightly but extends the bandwidth of the closed-loop system which is bene cial. e three-state feedforward control can o set the pole which is close to the virtual axis in the closed-loop system to extend the system bandwidth further.

Design of the ree-State Feedback Controller.
A mathematical model of the single hydraulic cylinder electrohydraulic servo system is established, and the open-loop transfer function of system is where k k is the open-loop gain of a system, ω k is the natural frequency of the hydraulic power mechanism, and ξ k is the damping ratio of the hydraulic power mechanism. Figure 9 illustrates the block diagram of a system after adding three-state feedback. In Figure 9, k df , k vf , and k af are the calculated parameters and represent the position, velocity, and acceleration feedback coe cients of the controller.
In accordance with the system block diagram, the closedloop transfer function of the three-state feedback system can be obtained as follows: e closed-loop transfer function after setting the expected three-state feedback is where ω r is the expected nal bandwidth of a system (system bandwidth after correction of three-state feedback and feedforward system) to achieve the nal dynamic characteristic requirements. ω r 157 rad/s (approximately 25 Hz) is selected. e value of ω c is generally 1.05-1.20 times that of the hydraulic natural frequency. In this study, ω c 1.2ω k and ξ c 0.7 are selected to ensure that the system has a su cient stability margin after the three-state feedback. e system transfer function G f (s) is the same as the expected transfer function W(s) after the three-state ϕ, θ, ψ   (22) and (23) can be used to solve the three-state feedback coe cients:

Design of the ree-State Feedforward
Controller. e three-state feedforward controller is added before a system completes the three-state feedback regulation, thereby eliminating the pole which is close to the virtual axis in the closed-loop transfer function. us, the system bandwidth can be extended. Let the transfer function of the three-state feedforward controller be B(s) k ap s 2 + k vp s + k dp , where k ap , k vp , and k dp are the determined controller parameters.
In equation (23), the three-state feedforward controller satis es equation (26) to o set the pole of the three-state feedback control system, which pole is close to the virtual axis: Equation (27) can be obtained by integrating equations (25) and (26). e parameters of the three-state feedback controller can be obtained by substituting the selected parameter values of ξ c and ω c into the following equation:

Experiment
When the construction vehicle drives on random roads, a test system collects the vibration information data of the vehicle and uses the attitude calculation algorithm in this study to obtain the vehicle vibration information. e experimental seat is installed on a 6-DOF vibration representation platform. Furthermore, the counterweight is fixed to the seat, and the vibration representation test bed is controlled by using real vehicle acquisition data, as depicted in Figure 10. Table 3 summarises the main parameters of the electrohydraulic servo system for the vehicle vibration representation test bed. e recorded attitude information after completing the test is compared with the attitude information of the moving platform. Because the measured value is relatively close to the target attitude, the motion curve within a short period is compared when the vibration is relatively intense. e attitude following curve of the system is demonstrated in Figure 11. e attitude following error of the vibration representation system is calculated, and the following error of the 6-DOF test platform is analysed. e attitude following error curve of the vibration representation system is exhibited in Figure 12. e position and attitude following curve and its error curve indicate the response lag of the vibration representation system. is phenomenon is mainly due to the inertia of the mechanical and hydraulic systems; such inertia causes a certain response delay to the input signal of the test platform. e response lag of the system is less than 10 ms, which satisfies the design requirements. A certain error exists between the position, attitude of the vibration platform, and the target signal. e attitude angle error is less than 0.5°, and the platform translation error is less than 2 mm. e errors are caused by two main reasons as follows: (1) further optimisation of control parameters and (2) assembly gap and error of each branched joint of the parallel platform. However, the errors of experimental results satisfy the design requirements, and it verifies the effectiveness of vibration information acquisition and representation control strategy.

Conclusions
(1) To increase accuracy of vehicle attitude measurement, an improved adaptive algorithm of multisensor attitude fusion is proposed. e algorithm uses weighted least squares method and complementary filtering fusion algorithm to fuse the information of accelerometer, magnetometer, and gyroscope and calculate the vehicle attitude. Experiments show that the algorithm can effectively fuse multisensor information and accurately obtain vehicle attitude.    (2) In order to obtain the vibration displacement of vehicle body, a displacement reconstruction algorithm based on acceleration measurement is proposed. e algorithm converts the acceleration measurement signal to the geodetic coordinate system according to vehicle attitude and then reconstructs the vibration displacement information from the acceleration signal by using the lowfrequency attenuation integral method. e algorithm can e ectively remove the trend term in the integration process, and the integration results are more accurate.
(3) A three-state controller is designed for the sixdegree-of-freedom vibration test bed. e damping ratio of system is improved by three-state feedback, and the dynamic characteristics of system are improved. en, the dominant poles of the closed-loop system are cancelled by three-state feedforward, which further expands the frequency range of the vibration test bed to reproduce road information. (4) e collected vibration signal of the test carrier is used as the excitation signal to test the vibration reproduction system. e experimental results show that the platform attitude angle error is less than 0.5 degrees and the displacement error is less than 2 mm. e system can meet the design requirements.
Data Availability e data used to support the ndings of this study are included within the article.

Conflicts of Interest
e authors declare that there are no con icts of interest in the publication of this paper.