Improvement of Ride Comfort of Truck Semi-Active Cab Suspension System Based on Accurate Modeling

Preprint

Article

Improvement of Ride Comfort of Truck Semi-Active Cab Suspension System Based on Accurate Modeling

Altmetrics

Downloads

102

Views

Comments

A peer-reviewed article of this preprint also exists.

Qihao Sun,

Changcheng Yin^*,Baohua Wang

Qihao Sun,

Changcheng Yin^*,Baohua Wang

This version is not peer-reviewed

Submitted:

06 August 2024

Posted:

08 August 2024

You are already at the latest version

Alerts

Abstract

The semi-active cab suspension system for trucks is gaining increasing importance due to its economic advantages, low energy consumption, and significant enhancement of ride comfort. This paper investigates the effects of three control methods on improving ride comfort of semi-active cab suspension systems under random and bump road conditions: Proportional-Integral-Derivative (PID) control, fuzzy PID control, and Model Predictive Control (MPC). Initially, an accurate multi-degree-of-freedom truck cab suspension model was developed and validated through actual road tests. Based on this model, three control strategies were designed and implemented. Finally, the effectiveness of each control strategy was evaluated under various road conditions, including random and bump road scenarios. The results demonstrate that these control strategies can effectively reduce vibrations and impacts, significantly enhancing ride comfort, which is crucial for reducing driver fatigue and improving driving safety.

Keywords:

Subject: Engineering - Automotive Engineering

1. Introduction

Trucks are composed of multiple intricate components, including the cab, frame, and wheels. During operation, uneven road surfaces can induce vibrations and shocks that are transmitted to the cab through the cab suspension , suspension, and tires, consequently impacting the driver’s comfort and driving safety. The primary function of the cab suspension system is to isolate road-induced vibrations, thereby reducing vibrations within the cab to enhance driver comfort and ensure driving safety. Therefore, investigating cab suspension systems is crucial for enhancing driving comfort and safety [1,2,3,4,5].

For the cab suspension system, research primarily focuses on three categories: passive suspension, semi-active suspension, and active suspension. Among these, semi-active suspension has garnered attention due to its combination of the advantages of both passive and active suspension systems. Compared to passive suspension, semi-active suspension demonstrates superior performance in enhancing driving comfort and vehicle stability [6],Compared to active suspension systems, semi-active suspension systems require only a small amount of external energy input, resulting in lower energy consumption and cost [7]. Control technology plays a crucial role in semi-active cab suspension systems. Superior control technology can mitigate vehicle jolts on uneven surfaces, effectively manage body posture, and enhance driving safety. Conversely, if the control technology is poorly implemented, its performance may be inferior to that of traditional passive suspension systems. Inappropriate control logic may lead to excessive adjustments, which not only reduces driving comfort but also increases the risk of driving. Therefore, developing high-performance control technologies is of utmost importance. Common control methods for semi-active cab suspensions include PID [8,9], sliding mode control [10,11], optimal control [5,12], and fuzzy control [13,14].

Optimal control techniques are widely used in the semi-active control systems of cab suspensions. Two primary methods, the Linear Quadratic Regulator (LQR) and Linear Quadratic Gaussian (LQG), achieve optimal performance by determining the weighted matrices for system states and control variables [15]. Xia et al. developed an LQG control strategy for truck active suspension systems to improve cargo integrity and other performance aspects. Their four-degree-of-freedom active suspension model showed significant enhancements in cargo integrity, ride comfort, and handling stability over passive suspension systems through time-domain and frequency-domain simulations [16]. Gandhi et al. evaluated the LQR’s effectiveness in a four-degree-of-freedom half-car suspension model using MATLAB-Simulink simulations, which indicated that LQR significantly improves ride quality and passenger comfort [17]. Yan and Zhang explored LQG control for heavy truck cab active suspension systems to boost ride comfort and safety. Their study, involving a half-truck cab active suspension model and a linear quadratic optimal control strategy, demonstrated marked reductions in vertical cab acceleration, roll angle, and pitch angle acceleration, enhancing overall ride comfort and driving safety [18]. Model Predictive Control (MPC) is an advanced strategy that guides current control decisions by predicting future system behavior. Unlike LQR and PID control, which only constrain input and output variables, MPC can handle various constraints on state variables, providing exceptional performance for complex systems. Lu et al. proposed a coordinated semi-active suspension system for trucks using MPC to improve ride comfort by addressing vertical, roll, and pitch vibrations. Co-simulations between ADAMS/Car and MATLAB/Simulink showed that MPC outperforms Skyhook control, optimizing truck cab dynamics and significantly enhancing driver comfort during long-haul transportation [19]. Yang et al. utilized an MPC algorithm that enhances vehicle suspension systems through road surface recognition. By tracking road parameters and employing a Bayesian approach to optimize suspension control, the MPC algorithm significantly reduces acceleration, suspension displacement, and tire displacement compared to passive suspension systems, thereby improving passenger comfort across various speeds and road conditions [20].

Fuzzy control is extensively utilized in semi-active cab suspension systems, primarily due to its ability to manage complex and nonlinear dynamic behaviors. This method enhances driving comfort and safety across diverse road conditions. Furthermore, fuzzy control is highly robust in addressing uncertainties and disturbances, which improves the system’s stability and reliability in challenging driving environments. Ning et al. introduced a semi-active control approach for active seat suspension using a disturbance observer-based Takagi-Sugeno (TS) fuzzy controller. This approach incorporates acceleration feedback to boost the controller’s performance and stability, compensating for disturbances like friction and model simplification. It emphasizes minimizing low-frequency vibrations to enhance ride comfort and reduce driver fatigue [21]. Nguyen et al. devised a PID-fuzzy control strategy for cab isolation systems, merging fuzzy and PID control to mitigate vibrations in the cab isolation frame. Research has shown significant reductions in the power spectral density (PSD) and weighted root mean square (RMS) values of vertical driver seat and pitch cab angle vibrations, especially under high-density elastoplastic soil conditions [22]. Liem et al. explored the use of fuzzy logic controllers in heavy truck cab isolation systems, aiming to improve ride comfort by reducing the RMS acceleration of the driver’s seat and the pitch and roll angles of the cab. Their findings demonstrated that RMS acceleration is considerably decreased across various road conditions, highlighting the effectiveness of fuzzy logic control in enhancing ride comfort and reducing vibrations compared to passive isolation systems [23].

In semi-active cab suspension control research, simplified multi-degree-of-freedom models are commonly employed. Nguyen et al. studied the ride comfort of an 8-degree-of-freedom heavy truck cab. Their study utilized an 8-degree-of-freedom dynamic model with Matlab/Simulink and fuzzy logic control, demonstrating significant improvements in ride comfort on both step and random road surfaces[24]. Akçay et al. developed a three-degree-of-freedom cab model using stochastic optimal control to analyze and improve truck driving dynamics. Their research compared LQG design with reduced-order controllers, showing that both effectively attenuated cab vibrations. HIFOO optimization presented a promising alternative, excelling in minimizing heave, pitch, and roll accelerations under random road excitation [25]. Basaran et al. proposed a seven-degree-of-freedom truck cab suspension model using an adaptive vector backstepping control method. The proposed control system aimed to improve the vibration control of the truck cab, significantly suppressing vibrations under various road conditions, thereby enhancing driver comfort and safety [26]. The truck cab suspension system is a complex system, and establishing an accurate dynamic model is crucial for designing suspension control strategies. However, current research on semi-active suspension models mostly lacks model validation. Therefore, the development of suspension control strategies for vehicle models based on real vehicle verification remains to be further investigated.

In this study, we developed an accurate ten-degree-of-freedom truck cab suspension dynamics model, validated through real vehicle testing. To enhance ride comfort, three control strategies were devised and their performances under different driving conditions were compared. The key contributions of this research are summarized as follows:

An accurate four-wheel road model was constructed, with actual road surface roughness of the test road measured;
The accuracy of the full vehicle cab suspension model was verified experimentally;
Three control strategies for semi-active suspension control were proposed;
The performance of the control strategies was evaluated on both random and bumpy roads.

The structure of this paper is as follows: Section 2 details the establishment and validation of the accurate semi-active cab suspension and road models. Section 3 describes the three proposed control strategies. Section 4 examines the performance of the various control strategies under different scenarios. Finally, Section 5 provides a summary.

2. Model and Verification

2.1. The Truck Cab Suspension Model

The truck cab suspension model discussed in this paper encompasses the cab, frame, wheels, cab suspension, and suspension system. As depicted in Figure 1, the model includes various degrees of freedom: the wheels possess four vertical degrees of freedom, the frame has three degrees of freedom (vertical motion, roll, and pitch), and the cab has three degrees of freedom (vertical motion, roll, and pitch). The model also takes into account the uneven effects of road disturbances on the front and rear wheels, as well as on the left and right sides. Compared to the traditional quarter-car model, this more comprehensive multi-degree-of-freedom vehicle model can more accurately represent the dynamic behavior of the vehicle under actual driving conditions [27,28].

According to Newton’s second law, the vertical vibration differential equation of the wheel assembly is illustrated by Equation (1):

\{\begin{matrix} m_{1} {\ddot{z}}_{w 1} = F_{w 1} - F_{k s 1} - F_{c s 1} \\ m_{2} {\ddot{z}}_{w 2} = F_{w 2} - F_{k s 2} - F_{c s 2} \\ m_{3} {\ddot{z}}_{w 3} = F_{w 3} - F_{k s 3} - F_{c s 3} \\ m_{4} {\ddot{z}}_{w 4} = F_{w 4} - F_{k s 4} - F_{c s 4} \end{matrix}

(1)

The frame is supported by four suspensions and, in turn, the cab is supported by four cab mounts. The differential equations governing the vertical, pitch, and roll motions of the frame are presented in Equation (3):

\{\begin{matrix} m_{f} {\ddot{z}}_{f} & = F_{k s 1} + F_{c s 1} + F_{k s 2} + F_{c s 2} \\ + F_{k s 3} + F_{c s 3} + F_{k s 4} + F_{c s 4} \\ - F_{k c s 1} - F_{c c s 1} - F_{k c s 2} - F_{c c s 2} \\ - F_{k c s 3} - F_{c c s 3} - F_{k c s 4} - F_{c c s 4} \\ l_{f x} {\ddot{θ}}_{f} & = (F_{k s 4} + F_{c s 4} + F_{k s 3} + F_{c s 3}) d_{2} \\ - (F_{k s 1} + F_{c s 1} + F_{k s 2} + F_{c s 2}) d_{1} \\ + (F_{k c s 1} + F_{c c s 1} + F_{k c s 2} + F_{c c s 2}) b_{1} \\ + (F_{k c s 3} + F_{c c s 3} + F_{k c s 4} + F_{c c s 4}) b_{2} \\ I_{f y} {\ddot{φ}}_{f} & = (F_{k s 1} + F_{c s 1} + F_{k s 4} + F_{c s 4}) d_{4} \\ - (F_{k s 2} + F_{c s 2} + F_{k s 3} + F_{c s 3}) d_{3} \\ - (F_{k c s 1} + F_{c c s 1} + F_{k c s 4} + F_{c c s 4}) d_{4} \\ + (F_{k c s 2} + F_{c c s 2} + F_{k c s 3} + F_{c c s 3}) b_{2} \end{matrix}

(2)

The cab, supported by four cab mounts, exhibits vertical, pitch, and roll motions, the control differential equations of which are demonstrated in Equation (3):

\{\begin{matrix} m_{f} {\ddot{z}}_{f} & = F_{k c s 1} + F_{c c s 1} + F_{k c s 2} + F_{c c s 2} \\ + F_{k c s 3} + F_{c c s 3} + F_{k c s 4} + F_{c c s 4} \\ I_{c x} {\ddot{θ}}_{c} & = (F_{k c s 3} + F_{c c s 3} + F_{k c s 4} + F_{c c s 4}) a_{2} \\ - (F_{k c s 1} + F_{c c s 1} + F_{k c s 2} + F_{c c s 2}) a_{1} \\ I_{c y} {\ddot{φ}}_{c} & = (F_{k c s 1} + F_{c c s 1} + F_{k c s 4} + F_{c c s 4}) a_{4} \\ - (F_{k c s 2} + F_{c c s 2} + F_{k c s 3} + F_{c c s 3}) a_{3} \end{matrix}

(3)

where

m_{i}

represents the mass of the four tire components, while

m_{f}

and

m_{c}

denote the masses of the vehicle frame and cab, respectively.

I_{f x}

and

I_{c x}

signify the moments of inertia of the frame and cab in the x-direction, and

I_{f y}

and

I_{c y}

denote the moments of inertia in the y-direction for the frame and cab, respectively. Parameters

a_{i}

b_{i}

, and

d_{i}

are geometric dimensions of the vehicle, with their specific meanings illustrated in Figure 1. The forces

F_{w i}

F_{k s i}

F_{c s i}

F_{k c s i}

F_{c c s i}

, and

F_{i}

correspond to the tire force, suspension stiffness force, suspension damping force, cab suspension stiffness force, cab suspension damping force, and semi-active suspension control force, respectively. Table 1 parameters enumerates the specific values for the vehicle’s mass, moments of inertia, and geometric dimensions.

The absolute displacement of the suspension is described by Equation (4):

\{\begin{matrix} z_{s t 1} = z_{f} + d_{1} θ_{f} + d_{4} φ_{f}, z_{c b 1} = z_{f} + b_{1} θ_{f} + d_{4} φ_{f}, z_{c t 1} = z_{c} + a_{1} θ_{c} + a_{4} φ_{c} \\ z_{s t 2} = z_{f} + d_{1} θ_{f} - d_{3} φ_{f}, z_{c b 2} = z_{f} + b_{1} θ_{f} - d_{3} φ_{f}, z_{c t 2} = z_{c} + a_{1} θ_{c} - a_{3} φ_{c} \\ z_{s t 3} = z_{f} - d_{2} θ_{f} - d_{3} φ_{f}, z_{c b 3} = z_{f} + b_{2} θ_{f} - d_{3} φ_{f}, z_{c t 3} = z_{c} - a_{2} θ_{c} - a_{3} φ_{c} \\ z_{s t 4} = z_{f} - d_{2} θ_{f} + d_{4} φ_{f}, z_{c b 4} = z_{f} + b_{2} θ_{f} + d_{4} φ_{f}, z_{c t 4} = z_{c} - a_{2} θ_{c} + a_{4} φ_{c} \end{matrix}

(4)

where

z_{s t i}

is the distance from the suspension to the vehicle’s center of mass,

z_{c b i}

is the distance from the lower support point of the cab suspension to the vehicle’s center of mass, representing geometric parameters, and

z_{c t i}

is the distance from the upper support point of the cab suspension to the center of mass of the cab. The variable q denotes the road excitation,

z_{w i}

is the vertical displacement of the suspension spring,

z_{f}

is the vertical displacement of the frame, and

z_{c}

is the vertical displacement of the cab. The relative displacements and velocities are provided in Section 6.

The mechanical characteristics of the damper exhibit significant nonlinearity, which is crucial to the overall vehicle performance. The damping force

F_{c}

as a function of relative velocity v is obtained by fitting experimental data using a cubic polynomial, as shown in Equation (5):

F_{c} = c_{1} v^{3} + c_{2} v^{2} + c_{3} v

(5)

where,

F_{c}

represents the damping force of the shock absorber. The coefficient

c_{i}

denotes the damping coefficients obtained through polynomial fitting. The specific parameters

c_{i}

are listed in Table 2.

To accurately describe the nonlinear characteristics of tires, this study employs a quadratic curve fitting method to analyze experimental data obtained under different tire pressure conditions. A quadratic polynomial is chosen for fitting, as shown in Equation (6):

F_{w} = k_{1} x^{2} + k_{2} x

(6)

where

F_{w}

represents the vertical tire force,

k_{i}

denotes the stiffness coefficients, and x represents the tire displacement.

The stiffness coefficients

k_{i}

exhibit a linear relationship with the tire pressure p, fitted using a linear polynomial as shown in Equation (7):

\{\begin{matrix} k_{1} = α_{1} p + α_{2} \\ k_{2} = α_{3} p + α_{4} \end{matrix}

(7)

The specific values of the coefficients

α_{i}

are detailed in Table 3.

2.2. Road Surface Model

2.2.1. Single Wheel Road Surface Model

The primary factor contributing to vibrations caused by vehicles is the roughness of the road surface, which can be characterized by its spatial power spectral density. This spatial frequency correlates with temporal frequency and varies with the vehicle’s speed [29,30]. When a vehicle travels at a speed u on a road surface with a spatial frequency n, the resulting temporal frequency f can be determined, as shown in Equation (8):

f = n u

(8)

The relationship between the spatial power spectral density and the temporal power spectral density is described by Equation (9):

G_{q} (f) = \frac{1}{u} G_{q} (n)

(9)

where

G_{q} (f)

is the temporal power spectral density, and

G_{q} (n)

is the spatial power spectral density.

The road’s power spectral density

G_{q} (n)

is defined by Equation (10):

G_{q} (n) = G_{q} (n_{0}) {(\frac{n}{n_{0}})}^{- w}

(10)

where

G_{q} (n_{0})

represents the roughness coefficient of the road surface, n is the spatial frequency,

n_{0}

is the reference spatial frequency (commonly set to 0.1), and w is the slope of the power spectral density, typically valued at -2.

By substituting Equation (9) into Equation (10), the temporal power spectral density can be expressed as Equation (11):

G_{q} (f) = \frac{1}{u} G_{q} (n_{0}) {(\frac{n}{n_{0}})}^{- 2} = G_{q} (n_{0}) n_{0}^{2} \frac{u}{f^{2}}

(11)

Expressing the vehicle speed u (in m/s) in terms of angular frequency

ω

(in rad/s), where

ω = 2 π f

, and substituting this into Equation (11) gives us Equation (12):

G_{q} (ω) = \frac{{(2 π)}^{2} G_{q} (n_{0}) n_{0}^{2} u}{ω^{2}}

(12)

When

ω = 0

G_{q} (ω)

approaches infinity, which is not feasible in practical engineering contexts. To accurately depict the road conditions at low frequencies, a spatial cutoff frequency

n_{m i n}

is introduced. Including this cutoff frequency into Equation (12) results in Equation (13):

G_{q} (ω) = {(2 π)}^{2} G_{q} (n_{0}) n_{0}^{2} \frac{u}{ω^{2} + ω_{0}^{2}}

(13)

Using the principles of random vibration, Equation (14) can be derived:

G_{q} (ω) = {| H (ω) |}^{2} S_{w}

(14)

where

H (ω)

is the frequency response function of the system, and

S_{w}

is the power spectral density of white noise, with

S_{w} = 1

. By combining Equations (13) and (14), we get the system’s frequency response function as expressed in Equation (15):

G_{q} (ω) = {| H (ω) |}^{2} S_{w}

(15)

Therefore, the time domain model for the road surface roughness, modeled as filtered white noise, is given by Equation (16):

\dot{q} (t) = - 2 π f_{0} \cdot q (t) + 2 π \sqrt{G_{q} (n_{0}) u} \cdot w (t)

(16)

2.2.2. Correlation of Road Inputs for Front and Rear Wheels

Assuming the front and rear wheel tracks of the vehicle are the same and the vehicle is moving in a straight line at a constant speed, the rear wheel road input is a time-delayed version of the front wheel road input. The relationship between the front and rear wheel road white noise inputs is given by Equation (17):

q_{4} (t) = q_{1} (t - t_{d})

(17)

The time delay

t_{d}

is defined by Equation (18):

t_{d} = \frac{L}{u}

(18)

where

q_{4} (t)

represents the rear wheel input,

q_{1} (t)

represents the front wheel input, L is the wheelbase, and u is the vehicle speed.

2.2.3. Correlation of Road Inputs for Left and Right Wheels

The correlation function of the road inputs for the left and right wheel tracks,

c o h (s)

, is given by Equation (19):

c o h (s) = \frac{q_{2} (s)}{q_{1} (s)}

(19)

Reference [52] presents a relatively simple and realistic mathematical model for four-wheel road input, described by Equation (20):

c o h (s) = e^{- \frac{B}{u} s}

(20)

where B is the vehicle’s wheel track and s is a fitting parameter.

By combining equ:coh and equ:cohh and using a first-order Pade approximation, we obtain Equation (21):

c o h (s) = \frac{q_{2} (s)}{q_{1} (s)} \approx \frac{1 - \frac{B}{2 v} s}{1 + \frac{B}{2 v} s}

(21)

Applying the inverse Laplace transform to the above equation, we get Equation (22):

{\dot{q}}_{2} (t) = - \frac{2 v}{B} q_{2} (t) + \frac{2 v}{B} q_{1} (t) - {\dot{q}}_{1} (t)

(22)

Thus, the trajectory equation for the left and right wheels can be expressed by Equation (23):

{\dot{q}}_{2} (t) = (2 π f_{0} + \frac{2 ν}{B}) q_{1} (t) - 2 π n_{0} \sqrt{G_{q} (n_{0}) ν} \cdot w (t) - \frac{2 ν}{B} q_{2} (t)

(23)

2.2.4. Four-Wheel Pavement Simulation Results

Based on the pavement model formula established in Section 2.2, the primary variables for different pavements are the pavement roughness coefficient

G_{q} (n_{0})

and vehicle speed u. According to the ISO 8608 standard, the range of roughness coefficients for various pavement grades is shown in Table 4.

Here, we utilized the MATLAB/Simulink platform to simulate vehicle travel on Class A pavement. According to Table 1, the power spectral density function for Class A pavement is represented as

G_{q} (n_{0}) = 16 \times 10^{- 6} m^{3}

. The vehicle speed u was set to 10 m/s.The time domain plot of the Class A pavement, as obtained from the simulation, is shown in Figure 2.

We conducted spectral analysis using the pwelch function in MATLAB. This function calculated the temporal power spectral density, which we then converted to spatial power spectral density using Equation (8), as illustrated in Figure 3.

From Figure 2 and Figure 3, it is evident that the constructed four-wheel pavement model effectively reflects the lag phenomenon between the front and rear wheels, as well as the correlation between the left and right wheels.

2.3. Measurement of Road Surface Roughness

In this study, the three-meter ruler method was employed to measure road surface roughness. The primary instruments used were a three-meter road inspection ruler and a feeler gauge, as illustrated in Figure 4. The feeler gauge with height markings was inserted into the gap to measure the maximum gap height, accurate to 0.2 millimeters.

The experimental road surface is located on Sanxia Road in Shiyan City, Hubei Province. The calculation method for single-rod road surface flatness involves taking the maximum gap value between the three-meter ruler and the road surface as the measurement result. Continuous measurements were taken 100 times along a 300-meter road section, and the average value was obtained, the measured value d =3 mm.

The Power Spectral Density (PSD) values were calculated following the method described in the literature [31]. By fitting the d values and the International Roughness Index (IRI) data, a conversion formula between the two was derived, as shown in Equation (24) below:

I_{I R I} = 0.9043 d - 0.0606

(24)

Based on the aforementioned relationship, the formula for calculating the Power Spectral Density

G_{q} (n_{0})

of the tested pavement is as follows:

G_{q} (n_{0}) = {(\frac{I_{IRI}}{0.6})}^{2}

(25)

Finally, the measured

G_{q} (n_{0})

of the tested road surface, rounded to the nearest whole number, is

20 \times 10^{- 6} m^{3}

2.4. Road Test Experiment

2.4.1. Experimental Equipment

In this experiment, the equipment used is illustrated in Figure 5. The data acquisition system utilizes a Siemens LMS data collector, operated in conjunction with Test.Lab software. The sensors employed are PCB Model 356A16 accelerometers from the United States, characterized by a ±50 g peak range and a 0.0001 g rms resolution. Calibration of the sensors follows the sensitivity calibration data provided by the manufacturer. The test vehicle is a Dongfeng Tianjin commercial truck. During testing, the driver ensures that the speed measurement error is within 0.5 m/s to meet the precision requirements of the experiment.

The accelerometers are mounted on the seat floor of the cab, and the data acquisition system is connected to a computer to verify the integrity and normalcy of data collection. After installation, it is necessary to calibrate the sensor’s measurement direction to align with the actual spatial coordinate system and adjust the measurement range to avoid data loss due to exceeding the range. The specific installation location is shown in Figure 6.

2.5. Data Filtering

During vehicle operation, road surface irregularities induce vertical vibrations in the vehicle, typically within the frequency range of 0.5 to 15 Hz. In contrast, high-frequency vibrations caused by short-wave stimuli such as engine vibrations fall within the range of 50 to 250 Hz. To effectively filter out these high-frequency vibrations, a second-order Infinite Impulse Response (IIR) digital band-pass filter can be employed, with cutoff frequencies set between 0.5 and 15 Hz. Figure 7 illustrates the effect of filtering on the acceleration signal within the cab of a vehicle traveling at a speed of 20 km/h.

From Figure 7, it is evident that the filter successfully eliminates high-frequency vibration noise (e.g., engine noise) while preserving the low-frequency vibration characteristics induced by road surface irregularities. This significantly enhances the usability and accuracy of the signal.

2.6. Model Validation

Ensuring the accuracy of the model is crucial in the development of a semi-active suspension system. This study validates the model’s accuracy by comparing the acceleration data generated from the simulation with the actual measured data. The specific method is illustrated in Figure 8, with the following steps:

Measure the actual road surface irregularities of the test road using a three-meter ruler.
Input the measured road surface irregularities into the road model.
Compare the acceleration signal obtained from the road model with the filtered measured acceleration signal.

The test vehicle was driven straight for 30 seconds at speeds of 20 km/h, 40 km/h, and 60 km/h on Sanxia Road in Shiyan City, Hubei Province. The data measured and simulated at different speeds are shown in Figure 9, Figure 10, and Figure 11:

To evaluate the model’s accuracy, the following formula was used to calculate the accuracy rate:

A_{c} = (1 - \frac{| M_{r m s} - S_{r m s} |}{M_{r m s}}) \times 100 %

(26)

where

A_{c}

represents the model accuracy,

M_{r m s}

denotes the measured RMS value of the acceleration signal, and

S_{r m s}

signifies the simulated RMS value of the acceleration signal. Table 5 illustrates the measured and simulated results at various speeds, along with their respective accuracies.

As demonstrated in Table 5, the proposed model exhibits high accuracy across different speeds, ranging from 84.66% to 96.99%. This confirms the robust accuracy of the proposed model.

3. Semi-Active Cab Suspension Control Strategies

In this section, we delve into the control strategies and constraints faced by the semi-active cab suspension system of commercial vehicles. To enhance ride comfort, we have designed and compared three control strategies: Proportional-Integral-Derivative (PID) control, Fuzzy PID control, and Model Predictive Control (MPC). These control strategies are constructed to optimize the vibration performance of the cab, addressing dynamic responses under varying road conditions.

3.1. Constraints of Semi-Active Cab Suspension Systems

When designing the controller for an active suspension system, multiple performance-limiting factors must be comprehensively considered to ensure stable operation under various conditions and to achieve the desired outcomes.

Mechanical Constraints on Suspension Movement: The suspension’s mechanical structure imposes limits on its movement. To avoid vehicle damage and ensure safety, the system must be designed to prevent the suspension from exceeding its physical limits. Specifically, the suspension deflection should be controlled within permissible limits to avoid excessive bottoming out, which can degrade ride comfort and potentially cause structural damage [32]. The limitation on suspension movement can be expressed by Equation (27):

$| z_{c t i} - z_{c b i} | \leq S_{m a x}$

(27)

where $S_{m a x}$ is the maximum allowable deflection of the suspension.
Actuator Force Limits: The active suspension system relies on actuators to adjust the suspension’s movement in real-time. However, the force that these actuators can exert is limited. To avoid damage and maintain system performance, the control force exerted by the actuators must not exceed their maximum limit [33]. This constraint is specified by Equation (28):

$| F_{i} | \leq F_{m a x}$

(28)

where $F_{m a x}$ is the maximum control force that the suspension actuators can provide.
Performance Objectives of the Active Suspension System: The primary performance objectives of the active suspension system include improving ride comfort by isolating the cab from road-induced shocks and maintaining road grip by ensuring continuous contact between the tires and the road. Achieving these objectives requires balancing conflicting demands. The static load limit of the tires is described by Equation (29):

$| F_{w i} | \leq F_{g}$

(29)

where $F_{g}$ is the static load of the tires.

3.2. PID Semi-active Cab Suspension Control Strategy

3.2.1. PID Control Optimization Objective

The primary objective of optimizing the semi-active suspension system is to minimize the vibrational acceleration response at the cab’s center of mass. To achieve this goal, the active controller must provide appropriate forces to adjust the cab’s posture. The vibration posture of the cab can be controlled in three directions: vertical, roll, and pitch, as illustrated in Figure 12.

Vertical Vibration Control: When the cab experiences vertical vibrations, the four dampers provide damping forces opposite to the direction of the vibrations. The vertical controller calculates the control force $F_{z}$ based on the acceleration on the cab.
Roll Vibration Control: When the cab undergoes roll vibrations, the two left-side dampers and two right-side dampers provide damping forces opposite to the direction of the vibrations. The roll controller calculates the control force $F_{θ}$ based on the roll angular acceleration on the cab.
Pitch Vibration Control: When the cab experiences pitch vibrations, the two front dampers and two rear dampers provide damping forces opposite to the direction of the vibrations. The pitch controller calculates the control force $F_{φ}$ based on pitch angular acceleration on the cab.

By performing logical force weighting calculations, the control forces from the three controllers are integrated to derive the optimal control force at each suspension point. The optimal control forces for the four suspension points can be expressed by Equation (30):

F_{i} = β_{1} F_{z} + β_{2} F_{θ} + β_{3} F_{φ} (i = 1, 2, 3, 4)

(30)

where

β_{i}

are the weighting coefficients.

3.2.2. PID Control Strategy Design

The PID controller is a classic control system algorithm widely applied in semi-active cab suspension systems. The PID controller calculates the control error through three components: proportional, integral, and derivative, and adjusts the control output to achieve the desired control objective. Equation (31) is given by:

F (t) = K_{p} e (t) + K_{i} \int_{0}^{t} e (t) d τ + K_{d} \frac{d e (t)}{d t}

(31)

where

F (t)

denotes the controller output, and

e (t)

represents the error at the current time. The system’s response characteristics are adjusted through three parameters: proportional gain

K_{p}

, integral gain

K_{i}

, and derivative gain

K_{d}

3.3. Fuzzy PID Semi-active Control Strategy for Cab Suspension System

Fuzzy control is a method based on fuzzy logic, suitable for handling complex and uncertain systems [34]. It formulates rules through expert experience and reasoning, simulating the fuzzy reasoning and decision-making processes of humans. In the design of the fuzzy control strategy, we select speed error e and acceleration error

e_{c}

as input variables, with output variables being

Δ k_{p}

Δ k_{i}

, and

Δ k_{d}

The value ranges for the variables

e_{\ddot{z}}

e_{c \ddot{z}}

e_{\ddot{θ}}

e_{c \ddot{θ}}

e_{\ddot{φ}}

, and

e_{c \ddot{φ}}

are

[- 1, 1]

[- 0.1, 0.1]

[- 0.6, 0.6]

[- 0.06, 0.06]

[- 0.6, 0.6]

, and

[- 0.06, 0.06]

respectively. The ranges for

Δ k_{p}

Δ k_{i}

, and

Δ k_{d}

are

[- 400, 400]

[- 20, 20]

, and

[- 4, 4]

respectively.

Input and output variables are fuzzified into seven membership levels: Positive Big (PB), Positive Medium (PM), Positive Small (PS), Zero (ZE), Negative Small (NS), Negative Medium (NM), and Negative Big (NB). These levels are represented by triangular membership functions. Next, fuzzy rules are designed, with some rules shown in Table 4. The fuzzy inference uses the Mamdani lookup table (LUT) method, and defuzzification is performed using the centroid method. The expression for the centroid method is as follows:

y^{*} = \frac{\sum_{i = 1}^{n} μ_{b_{i}} (y_{i}) \cdot y_{i}}{\sum_{i = 1}^{n} μ_{b_{i}} (y_{i})}

(32)

where

y^{*}

is the precise value after defuzzification,

μ_{b_{i}} (y_{i})

is the membership degree of the output fuzzy set, and

y_{i}

is the value of the output variable. The selected fuzzy rules are shown in Table 6.

The correction values of the PID parameters are calculated by Equation (33):

\{\begin{matrix} k_{p}^{'} = k_{p 0} + Δ k_{p} \\ k_{i}^{'} = k_{i 0} + Δ k_{i} \\ k_{d}^{'} = k_{d 0} + Δ k_{d} \end{matrix}

(33)

where

k_{p 0}

k_{i 0}

, and

k_{d 0}

are the initial PID parameters;

Δ k_{p}

Δ k_{i}

, and

Δ k_{d}

are the increments after fuzzification and fuzzy inference of e and

e_{c}

3.4. Model Predictive Control

Model Predictive Control (MPC) employs dynamic models to forecast the future behavior of controlled objects and accomplishes the control process through three steps [35]: predictive modeling, rolling optimization, and feedback correction. MPC is particularly advantageous in handling multi-input multi-output systems. The fundamental components of MPC consist of three main stages: prediction, optimization, and implementation.

Prediction Phase: Based on the current system state and predetermined control laws, the system’s output over a future horizon is predicted.
Optimization Phase: An optimal control sequence is derived by solving a finite-horizon optimization problem.
Implementation Phase: The first control input of the optimized sequence is applied to the system, and the entire process is repeated at the next sampling instant using updated measurements.

In MPC, the state-space equations of the system are first required. By linearizing the full-vehicle cab suspension system described in Section 2, the state-space representation of the system can be expressed as follows:

\{\begin{matrix} \dot{x} (t) = A x (t) + B_{u} u (t) + B_{w} ω (t) \\ y (t) = C x (t) + D_{u} u (t) + D_{w} ω (t) \end{matrix}

(34)

where, The matrices A and C are the system matrices;

B_{u}

is the input matrix;

D_{u}

is the output matrix;

B_{w}

and

D_{w}

are the disturbance matrices.

x (t) = {[z_{w 1}, z_{w 2}, z_{w 3}, z_{w 4}, z_{f}, θ_{f}, φ_{f}, z_{c}, θ_{c}, φ_{c}, {\dot{z}}_{w 1}, {\dot{z}}_{w 2}, {\dot{z}}_{w 3}, {\dot{z}}_{w 4}, {\dot{z}}_{f}, {\dot{θ}}_{f}, {\dot{φ}}_{f}, {\dot{z}}_{c}, {\dot{θ}}_{c}, {\dot{φ}}_{c}]}^{T}

(35)

u (t) = {[F_{1}, F_{2}, F_{3}, F_{4}]}^{T}

(36)

ω (t) = {[q_{1}, q_{2}, q_{3}, q_{4}]}^{T}

(37)

where

x (t)

represents the system states,

u (t)

denotes the control inputs, which are the control forces calculated by the controller to suppress vibrations of the suspension and wheel hub motors,

ω (t)

represents the road disturbances.

The system output is defined as the vertical acceleration, pitch angular acceleration, and roll angular acceleration:

y (t) = {[{\ddot{z}}_{f}, {\ddot{θ}}_{f}, {\ddot{φ}}_{f}]}^{T}

(38)

Since the optimization problem in MPC must be performed in the discrete-time domain, the continuous-time model is discretized using the zero-order hold method with a sampling time

T_{s}

, yielding the discrete-time model represented as:

\{\begin{matrix} x_{k + 1} = A_{d} x_{k} + B_{d w} w_{k} + B_{d u} u_{k} \\ y_{k} = C_{d} x_{k} + D_{d w} w_{k} + D_{d u} u_{k} \end{matrix}

(39)

where

A_{d}

B_{d u}

B_{d w}

C_{d}

D_{d u}

, and

D_{d w}

are the matrices of the discrete state-space representation.

The essence of the optimal control in MPC lies in the design of the cost function. The optimization objective function is given by Equation (40):

J = \sum_{i = 1}^{N} {(y (k + i) - y_{r e f} (k + i))}^{T} Q_{i} (y (k + i) - y_{r e f} (k + i)) + \sum_{i = 0}^{N - 1} u {(k + i)}^{T} R_{i} u (k + i)

(40)

where

y_{r e f}

denotes the reference signal, and

Q_{i}

and

R_{i}

are symmetric positive definite weighting matrices.

4. Experiments and Results

This section evaluates the effectiveness of the three control strategies described in Section 3 in terms of ride comfort for a semi-active cab suspension system. Section 3 provided detailed control design, while Section 2 offered vehicle modeling and parameter information. To assess the proposed controllers’ effectiveness, we selected vertical acceleration, pitch acceleration, and roll acceleration as comprehensive evaluation metrics. All analyses were conducted using MATLAB.

4.1. Random Road Surface

To verify the controllers’ performance, a road model was established as described in Section 2.2. Common road types and speeds were selected for evaluation: Class A road (general highway) and Class B road (rural concrete road), with speeds of 10 m/s and 5 m/s, respectively.

4.1.1. Class A Road at 10 m/s

In practice, typical highways are classified as Class A roads, tested at a speed of 10 m/s. Figure 13 respectively shows the performance of the three different controllers in terms of vertical acceleration (

\ddot{z}

), pitch acceleration (

\ddot{θ}

), and roll acceleration (

\ddot{φ}

4.1.2. Class B Road at 5 m/s

For Class B roads, representing rural concrete roads, the evaluation was conducted at a speed of 5 m/s. Figure 14 summarizes the test results of each controller concerning vertical acceleration (

\ddot{z}

), pitch acceleration (

\ddot{θ}

), and roll acceleration (

\ddot{φ}

)..

4.1.3. Performance Comparison

To quantitatively compare the performance of the control strategies, the RMS values of the accelerations were calculated. Table 7 and Table 8 present the RMS values for different road types and speeds.

Based on the analysis of Table 7, we can conclude that on Class A roads, MPC control significantly outperforms other control methods in reducing

\ddot{z}

\ddot{θ}

\ddot{φ}

. Specifically, compared to passive control, MPC control reduces vertical acceleration by 47%, by 19% compared to PID control, and by 13% compared to fuzzy PID control. Regarding pitch acceleration, MPC control reduces it by 39% compared to passive control, by 25% compared to PID control, and by 6% compared to fuzzy PID control. For roll acceleration, MPC control reduces it by 38% compared to passive control, by 19% compared to PID control, and by 18% compared to fuzzy PID control.

The analysis results of Table 8 indicate that on Class B roads, MPC control also performs best in reducing vehicle longitudinal acceleration, pitch acceleration, and roll acceleration. Specifically, MPC control reduces longitudinal acceleration by 45% compared to passive control, by 12% compared to PID control, and by 9% compared to fuzzy PID control. In terms of pitch acceleration, MPC control reduces it by 53% compared to passive control, by 31% compared to PID control, and by 9% compared to fuzzy PID control. Regarding roll acceleration, MPC control reduces it by 49% compared to passive control, by 34% compared to PID control, and by 25% compared to fuzzy PID control.

The analysis of Table 8 indicates that on Class B roads, MPC control also performs best in reducing vehicle longitudinal acceleration, pitch acceleration, and roll acceleration. Specifically, MPC control reduces vertical acceleration by 45% compared to passive control, by 12% compared to PID control, and by 9% compared to fuzzy PID control. In terms of pitch acceleration, MPC control reduces it by 53% compared to passive control, by 31% compared to PID control, and by 9% compared to fuzzy PID control. Regarding roll acceleration, MPC control reduces it by 49% compared to passive control, by 34% compared to PID control, and by 25% compared to fuzzy PID control.

On random roads, MPC control consistently excels in reducing vehicle acceleration, significantly outperforming passive control, PID control, and fuzzy PID control.

4.2. Bump response

The bumpy road is described by Equation (1) [36]:

Figure 15 shows

\ddot{z}

\ddot{θ}

\ddot{φ}

test results of each controller on the bumpy road. Table 9 lists the RMS values of

\ddot{z}

\ddot{θ}

and

\ddot{φ}

on the bumpy road.

The analysis of Table 9 indicates that on bumpy roads, MPC control also performs best in reducing vehicle vertical acceleration (

\ddot{z}

), pitch acceleration (

\ddot{θ}

), and roll acceleration (

\ddot{φ}

). Specifically, compared to passive control, MPC control reduces vertical acceleration by 45%, by 12% compared to PID control, and by 9% compared to fuzzy PID control. In terms of pitch acceleration, MPC control reduces it by 53% compared to passive control, by 31% compared to PID control, and by 9% compared to fuzzy PID control. Regarding roll acceleration, MPC control reduces it by 49% compared to passive control, by 34% compared to PID control, and by 25% compared to fuzzy PID control.

It is noteworthy that in the control of roll acceleration, a phenomenon of negative optimization was observed, indicating that in some cases, PID and fuzzy PID control might perform worse than passive control. However, MPC control still exhibits significant robustness in this aspect, effectively reducing roll acceleration even in the face of negative optimization.

5. Conclusions

This study investigates the effectiveness of three control methods in enhancing ride comfort for truck semi-active cab suspension system under random and bumpy road conditions: Proportional-Integral-Derivative control, Fuzzy PID control, and Model Predictive Control .

The main conclusions are as follows:

A single-wheel and four-wheel road model based on the power spectral density of road roughness was established. The simulations verified the lag between the front and rear wheels and the correlation between the left and right wheels. The actual road roughness coefficient of the test road was obtained using the three-meter straightedge method, with the road surface irregularity $G_{q} (n_{0})$ is rounded to $20 \times 10^{- 6} m^{3}$ .;
The accuracy of the cab suspension system model was experimentally validated by comparing the measured and simulated vertical acceleration data of the cab at different speeds. The model accuracy ranged from 84.66% to 96.99%, demonstrating high precision;
Three control strategies for the truck’s semi-active cab suspension were designed based on the model: PID control, Fuzzy PID control, and MPC control;
The validation results on both random and bumpy roads indicate that the proposed control strategy significantly enhances ride comfort. Among the various methods, Model Predictive Control (MPC) demonstrated superior performance and robustness, reducing vibrations by at least 31% without causing any adverse lift effects. This performance is notably better than passive control, PID control, and fuzzy PID control.

In conclusion, the proposed semi-active cab suspension control strategies significantly enhance ride comfort, with MPC showing the best performance. Future research will focus on developing more intelligent and adaptive control systems to address the uncertainties in vehicle parameters.

Author Contributions

Conceptualization, Methodology, Supervision, Writing—original draft, Software, Q.S.; Data curation, Writing—review and editing, C.Y.; Investigation, Validation, B.W. All authors have read and agreed to the published version of the manuscript.

Funding

Supported by the National Natural Science Foundation of China (52072116) and Science and Technology Program of Hubei Province (2023BEB047).

Data Availability Statement

The data are contained within this article.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

The relative displacement of the cab suspension (

x_{c s}

)is expressed as Equation (A1):

x_{c s} = z_{w i} - z_{s t i}

(A1)

The relative displacement of the suspension (

x_{s}

) is defined as Equation (A2):

x_{s} = z_{c b i} - z_{c t i}

(A2)

The absolute velocity is given by Equation (A3):

\{\begin{matrix} {\dot{z}}_{s t 1} = {\dot{z}}_{f} + d_{1} {\dot{θ}}_{f} + d_{4} {\dot{φ}}_{f}, {\dot{z}}_{c b 1} = {\dot{z}}_{f} + b_{1} {\dot{θ}}_{f} + d_{4} {\dot{φ}}_{f}, {\dot{z}}_{c t 1} = {\dot{z}}_{c} + a_{1} {\dot{θ}}_{c} + a_{4} {\dot{φ}}_{c} \\ {\dot{z}}_{s t 2} = {\dot{z}}_{f} + d_{1} {\dot{θ}}_{f} - d_{3} {\dot{φ}}_{f}, {\dot{z}}_{c b 2} = {\dot{z}}_{f} + b_{1} {\dot{θ}}_{f} - d_{3} {\dot{φ}}_{f}, {\dot{z}}_{c t 2} = {\dot{z}}_{c} + a_{1} {\dot{θ}}_{c} - a_{3} {\dot{φ}}_{c} \\ {\dot{z}}_{s t 3} = {\dot{z}}_{f} - d_{2} {\dot{θ}}_{f} - d_{3} {\dot{φ}}_{f}, {\dot{z}}_{c b 3} = {\dot{z}}_{f} + b_{2} {\dot{θ}}_{f} - d_{3} {\dot{φ}}_{f}, {\dot{z}}_{c t 3} = {\dot{z}}_{c} - a_{2} {\dot{θ}}_{c} - a_{3} {\dot{φ}}_{c} \\ {\dot{z}}_{s t 4} = {\dot{z}}_{f} - d_{2} {\dot{θ}}_{f} - d_{4} {\dot{φ}}_{f}, {\dot{z}}_{c b 4} = {\dot{z}}_{f} + b_{2} {\dot{θ}}_{f} + d_{4} {\dot{φ}}_{f}, {\dot{z}}_{c t 4} = {\dot{z}}_{c} - a_{2} {\dot{θ}}_{c} + a_{4} {\dot{φ}}_{c} \end{matrix}

(A3)

The relative velocity of the cab suspension(

v_{c s}

)is given by Equation (A4):

v_{c s} = {\dot{z}}_{w i} - {\dot{z}}_{s t i}

(A4)

The relative velocity of the suspension(

v_{s}

)is defined as Equation (A5):

v_{s} = {\dot{z}}_{c b i} - {\dot{z}}_{c t i}

(A5)

References

Zhang, J.; Guo, P.; Lin, J.; Wang, K. A Mathematical Model for Coupled Vibration System of Road Vehicle and Coupling Effect Analysis. Appl. Math. Model. 2016, 40, 1199–1217. [Google Scholar] [CrossRef]
Xu, P.; Bernardo, B.; Tan, K. Optimal Mounting Design for Cab Vibration Isolation. Int. J. Veh. Des. 2011, 57, 292. [Google Scholar] [CrossRef]
ElMadany, M. The Performance of Passive Cab Suspension Systems in Tractor-Semitrailer Vehicles. J. King Saud Univ. - Eng. Sci. 1990, 2, 131–151. [Google Scholar] [CrossRef]
Nguyen, V.; Zhang, J.; Le, V.; Jiao, R. Vibration Analysis and Modeling of an Off-Road Vibratory Roller Equipped with Three Different Cab’s Isolation Mounts. Shock Vib. 2018, 2018, 1–17. [Google Scholar] [CrossRef]
Ahn, D.V.; Kim, K.; Oh, J.; Seo, J.; Lee, J.W.; Park, Y.J. Optimal Control of Semi-Active Suspension for Agricultural Tractors Using Linear Quadratic Gaussian Control. Sens. 2023, 23, 6474. [Google Scholar] [CrossRef] [PubMed]
Bashir, A.O.; Rui, X.; Zhang, J. Ride Comfort Improvement of a Semi-active Vehicle Suspension Based on Hybrid Fuzzy and Fuzzy-PID Controller. Stud. Inform. Control 2019, 28, 421–430. [Google Scholar] [CrossRef]
Gui, L.m.; Shi, W.k.; Liu, W. A Semi-Active Suspension Design for off-Road Vehicle Base on Magneto-rheological Technology. 2012 9th International Conference on Fuzzy Systems and Knowledge Discovery; IEEE: Chongqing, Sichuan, China, 2012; pp. 2565–2568. [Google Scholar] [CrossRef]
Tian, M.; Nguyen, V. Control Performance of Suspension System of Cars with PID Control Based on 3D Dynamic Model. J. Mech. Eng. Autom. Control Syst. 2020, 1, 1–10. [Google Scholar] [CrossRef]
Ni, S.; Nguyen, V. Performance of Semi-Active Cab Suspension System with Different Control Methods. J. Mechatron. Artif. Intell. Eng. 2023, 4, 8–17. [Google Scholar] [CrossRef]
Maciejewski, I.; Blazejewski, A.; Pecolt, S.; Krzyzynski, T. A Sliding Mode Control Strategy for Active Horizontal Seat Suspension under Realistic Input Vibration. J. Vib. Control 2023, 29, 2539–2551. [Google Scholar] [CrossRef]
Zhao, Z.; Wang, C.; Zhao, J.; Du, W. LQR Force Command Planning–Based Sliding Mode Control for Active Suspension System. Proc. Inst. Mech. Eng. I: J. Syst. Control Eng. 2024, 238, 373–385. [Google Scholar] [CrossRef]
Lu, Y.; Khajepour, A.; Soltani, A.; Li, R.; Zhen, R.; Liu, Y.; Wang, M. Gain-Adaptive Skyhook-LQR: A Coordinated Controller for Improving Truck Cabin Dynamics. Control Eng. Pract. 2023, 130, 105365. [Google Scholar] [CrossRef]
Gheibollahi, H.; Tehrani, M.M.; Najafi, A. Improving Ride Comfort Approach by Fuzzy and Genetic-Based PID Controller in Active Seat Suspension. Int. J. Autom. Control 2024, 18, 184–213. [Google Scholar] [CrossRef]
Zhang, B.; Liu, M.; Wang, K.; Tan, B.; Deng, Y.; Qin, A.; Liu, J. Takagi–Sugeno Fuzzy Model-Based Control for Semi-Active Cab Suspension Equipped with an Electromagnetic Damper and an Air Spring. Machines 2023, 11, 226. [Google Scholar] [CrossRef]
He, Q.; Li, H.; Zhang, R.; Jia, T. Study on the LQR Control of High-speed Elevator Car Horizontal Vibration Based on the Jumping Inertia Weight Particle Swarm Optimization. Int. J. Acoust. Vib. 2022, 27, 122–137. [Google Scholar] [CrossRef]
Xia, R.x.; Li, J.h.; He, J.; Shi, D.f.; Zhang, Y. Linear-Quadratic-Gaussian Controller for Truck Active Suspension Based on Cargo Integrity. Adv. Mech. Eng. 2015, 7, 168781401562032. [Google Scholar] [CrossRef]
Gandhi, P.; Adarsh, S.; Ramachandran, K. Performance Analysis of Half Car Suspension Model with 4 DOF Using PID, LQR, FUZZY and ANFIS Controllers. Procedia Comput. Sci. 2017, 115, 2–13. [Google Scholar] [CrossRef]
Yan, G.H.; Zhang, S. Research on Modeling and Optimization Control of Heavy Truck Cab Active Suspension System. Appl. Mech. Mater. 2014, 687–691, 359–362. [Google Scholar] [CrossRef]
Lu, Y.; Khajepour, A.; Soltani, A. A Coordinated Control System for Truck Cabin Suspension Based on Model Predictive Control. Int. J. Heavy Veh. Syst. 2022, 29, 518. [Google Scholar] [CrossRef]
Yang, T. Active Suspension Control Strategy for Vehicles Based on Road Surface Recognition. Nonlinear Dyn 2024, 112, 11043–11065. [Google Scholar] [CrossRef]
Ning, D.; Sun, S.; Zhang, F.; Du, H.; Li, W.; Zhang, B. Disturbance Observer Based Takagi-Sugeno Fuzzy Control for an Active Seat Suspension. Mech. Syst. Sig. Process. 2017, 93, 515–530. [Google Scholar] [CrossRef]
Nguyen, V.; Jiao, R.; Le, V.; Hoang, A. Performance of PID-Fuzzy Control for Cab Isolation Mounts of Soil Compactors. Math. Models Eng. 2019, 5, 137–145. [Google Scholar] [CrossRef]
Liem, N.V.; Jianrun, Z.; Quynh, L.V.; Renqiang, J. Study of Fuzzy Control for Cab’s Isolation System of Heavy Truck. Vibroengineering Procedia 2016, 10, 309–314. [Google Scholar]
Nguyen, V.L.; Le, V.Q. Ride Comfort Performance of Heavy Truck with Three Control Cases of Semi-Active Isolation Systems. Vibroengineering Procedia 2019, 22, 93–98. [Google Scholar] [CrossRef]
Akçay, H.; Türkay, S. Stochastic Optimal Control of Truck Cabin with Active Suspension. Int. J. Heavy Veh. Syst. 2014. [Google Scholar] [CrossRef]
Basaran, S.; Basaran, M. Vibration Control of Truck Cabins With the Adaptive Vectorial Backstepping Design of Electromagnetic Active Suspension System. IEEE Access 2020, 8, 173056–173067. [Google Scholar] [CrossRef]
Yan, S.; Sun, W.; Gao, H. Vibration Suppression for Motor-Driven Uncertain Active Suspensions With Hard Constraints and Analysis of Energy Consumption. IEEE Trans. Syst., Man, Cybern., Syst. 2022, 52, 618–630. [Google Scholar] [CrossRef]
Yu, M.; Evangelou, S.A.; Dini, D. Advances in Active Suspension Systems for Road Vehicles. Engineering-london. 2024, 33, 160–177. [Google Scholar] [CrossRef]
Wang, C.; Cui, X.; Zhao, S.; Zhou, X.; Song, Y.; Wang, Y.; Guo, K. A Deep Reinforcement Learning-Based Active Suspension Control Algorithm Considering Deterministic Experience Tracing for Autonomous Vehicle. Appl. Soft Comput. 2024, 153, 111259. [Google Scholar] [CrossRef]
Qin, A.; Zhang, B.; Ning, D.; Tan, B.; Du, H. A Self-Sensing Approach for Estimating Suspension Displacement and Velocity in Semi-Active Electromagnetic Dampers. Mech. Syst. Sig. Process. 2024, 208, 111049. [Google Scholar] [CrossRef]
Li, P.; Hu, G.; Xia, H.; Guo, R. Efficient Method Based on Recurrent Neural Networks for Pavement Evenness Detection. Measurement 2023, 212, 112676. [Google Scholar] [CrossRef]
Zhao, R.; Xie, H.; Gong, X.; Sun, X.; Cao, C. Neural Network-Based Adaptive Height Tracking Control of Active Air Suspension System with Magnetorheological Fluid Damper Subject to Uncertain Mass and Input Delay. Sens. 2023, 24, 156. [Google Scholar] [CrossRef] [PubMed]
Krauze, P. Identification of Control-Related Signal Path for Semi-Active Vehicle Suspension with Magnetorheological Dampers. Sens. 2023, 23, 5770. [Google Scholar] [CrossRef] [PubMed]
Dai, D.; Zhang, J.; Zhang, B.; Li, P.; Hu, W. Adaptive Hierarchical Optimization Control for Electrohydraulic Suspension with Resistor-Capacitor Operator. Appl. Math. Model. 2024, 126, 606–624. [Google Scholar] [CrossRef]
Zhang, N.; Yang, S.; Wu, G.; Ding, H.; Zhang, Z.; Guo, K. Fast Distributed Model Predictive Control Method for Active Suspension Systems. Sens. 2023, 23, 3357. [Google Scholar] [CrossRef]
Wong, P.K.; Li, W.; Ma, X.; Yang, Z.; Wang, X.; Zhao, J. Adaptive Event-Triggered Dynamic Output Feedback Control for Nonlinear Active Suspension Systems Based on Interval Type-2 Fuzzy Method. Mech. Syst. Sig. Process. 2024, 212, 111280. [Google Scholar] [CrossRef]

Figure 1. Schematic diagram of the truck cab suspension model.

Figure 2. The Class A road.

Figure 3. The spatial power spectral density of the Class A road.

Figure 4. The three-meter road inspection ruler.

Figure 5. Experimental equipment: (a) Dongfeng Tianjin experimental truck. (b) LMS data acquisition equipment. (c) Triaxial accelerometer.

Figure 6. Sensor installation position.

Figure 7. Experimental data filtering.

Figure 8. Model validation methods.

Figure 9. Simulation Data Comparison with Actual Measurement Data at a Speed of 20 km/h on Grade A Road Surface.

Figure 10. Simulation Data Comparison with Actual Measurement Data at a Speed of 40 km/h on Grade A Road Surface.

Figure 11. Simulation Data Comparison with Actual Measurement Data at a Speed of 60 km/h on Grade A Road Surface.

Figure 12. Analysis of vertical, roll, and pitch vibration of the cab.

Figure 13. Performance Comparison of Three Control Strategies on Grade A Road Surface: (a) Vertical Acceleration. (b) Pitch Angle Acceleration. (c) Roll Angle Acceleration.

Figure 14. Performance Comparison of Three Control Strategies on Grade B Road Surface: (a) Vertical Acceleration. (b) Pitch Angle Acceleration. (c) Roll Angle Acceleration.

Figure 15. Comparison of Control Methods on the bumpy road.

Table 1. Vehicle parameters.

Parameter Name	Value	Unit
Geometric Parameter $d_{1}$	2.4	m
Geometric Parameter $d_{2}$	2.6	m
Geometric Parameter $d_{3}$	0.6	m
Geometric Parameter $d_{4}$	0.6	m
Geometric Parameter $b_{1}$	3.2	m
Geometric Parameter $b_{2}$	2	m
Geometric Parameter $a_{1}$	0.7	m
Geometric Parameter $a_{2}$	0.8	m
Geometric Parameter $a_{3}$	0.5	m
Front Tire Mass $m_{2}$	60	kg
Rear Tire Mass $m_{3}$	90	kg
Frame Mass $m_{f}$	3150	kg
Cabin Mass $m_{c}$	980	kg
Frame Moment of Inertia in x-direction $I_{f x}$	813	m²·kg
Frame Moment of Inertia in y-direction $I_{f y}$	764	m²·kg
Cabin Moment of Inertia in x-direction $I_{c x}$	346	m²·kg
Cabin Moment of Inertia in y-direction $I_{c y}$	2264	m²·kg

Table 2. Shock absorber parameters.

Shock Absorber Position and Condition	$c_{1}$	$c_{2}$	$c_{3}$
Cabin Front Suspension - Compression	0.000031	0.027476	11.389424
Cabin Front Suspension - Tension	-0.000017	0.011704	6.975855
Cabin Rear Suspension - Compression	-0.000023	-0.007452	4.482131
Cabin Rear Suspension - Tension	0.000032	-0.022581	6.182971
Front Suspension - Compression	0.000022	0.021163	7.101091
Front Suspension - Tension	0.000073	-0.080829	29.347936

Table 3. The coefficients

α_{i}

Table 3. The coefficients

α_{i}

Parameter	Value
$α_{1}$	0.3857
$α_{2}$	372.18
$α_{3}$	0.0138
$α_{4}$	-4.8361

Table 4. Road Surface Levels.

Road Surface Level	$G_{q} (n_{0})$ Upper Limit	$G_{q} (n_{0})$ Lower Limit	$G_{q} (n_{0})$ Geometric Mean
A	32	8	16
B	128	32	64
C	512	128	256
D	2048	512	1024
E	8192	2048	4096
F	32768	8192	16384
G	131072	32768	65536
H	524288	131072	262144

Table 5. Comparison of Simulation Data with Actual Measurement Data.

Speed	$M_{r m s}$	$S_{r m s}$	$A_{c}$
20 km/h	261.8304	221.6586	84.66%
40 km/h	408.7856	385.043	94.19%
60 km/h	457.6561	471.4536	96.99%

Table 6. Fuzzy Logic Table.

e/ ec	NB	NM	NS	ZE	PS	PM	PB
NB	PB	PB	PM	PM	PS	ZE	ZE
NM	PB	PB	PM	PS	PS	ZE	NS
NS	PM	PM	PM	PS	PS	ZE	NS
ZE	PM	PM	PS	ZE	NS	NS	NM
PS	PS	PS	ZE	NS	NS	NM	NM
PM	PS	ZE	NS	NS	NM	NM	NS
PB	ZE	ZE	NM	NM	NB	NB	NB

Table 7. Comparison of Control Methods on Grade A Road Surface.

Control Method	Passive	PID	Fuzzy PID	MPC
$\ddot{z} (R M S)$	0.27379	0.17853	0.16679	0.14448
$\ddot{θ} (R M S)$	0.095454	0.077903	0.062371	0.058613
$\ddot{φ} (R M S)$	0.097932	0.075125	0.073781	0.060532

Table 8. Comparison of Control Methods on Grade B Road Surface.

Control Method	Passive	PID	Fuzzy PID	MPC
$\ddot{z} (R M S)$	0.31288	0.19592	0.19078	0.17313
$\ddot{θ} (R M S)$	0.23141	0.15958	0.1206	0.10971
$\ddot{φ} (R M S)$	0.31039	0.23677	0.20873	0.15711

Table 9. Comparison of Control Methods at 15m/s.

Control Method	Passive	PID	Fuzzy PID	MPC
$\ddot{z} (R M S)$	0.25207	0.13749	0.131	0.11558
$\ddot{θ} (R M S)$	0.05446	0.04205	0.035943	0.032332
$\ddot{φ} (R M S)$	0.011419	0.012213	0.01207	0.0078793

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

MDPI Initiatives

Important Links

Choose an area of interest and we will send you notifications of new preprints at your preferred frequency.

Disclaimer

e/ ec	NB	NM	NS	ZE	PS	PM	PB
NB	PB	PB	PM	PM	PS	ZE	ZE
NM	PB	PB	PM	PS	PS	ZE	NS
NS	PM	PM	PM	PS	PS	ZE	NS
ZE	PM	PM	PS	ZE	NS	NS	NM
PS	PS	PS	ZE	NS	NS	NM	NM
PM	PS	ZE	NS	NS	NM	NM	NS
PB	ZE	ZE	NM	NM	NB	NB	NB

e/ ec	NB	NM	NS	ZE	PS	PM	PB
NB	PB	PB	PM	PM	PS	ZE	ZE
NM	PB	PB	PM	PS	PS	ZE	NS
NS	PM	PM	PM	PS	PS	ZE	NS
ZE	PM	PM	PS	ZE	NS	NS	NM
PS	PS	PS	ZE	NS	NS	NM	NM
PM	PS	ZE	NS	NS	NM	NM	NS
PB	ZE	ZE	NM	NM	NB	NB	NB

Improvement of Ride Comfort of Truck Semi-Active Cab Suspension System Based on Accurate Modeling

Abstract

1. Introduction

2. Model and Verification

2.1. The Truck Cab Suspension Model

2.2. Road Surface Model

2.2.1. Single Wheel Road Surface Model

2.2.2. Correlation of Road Inputs for Front and Rear Wheels

2.2.3. Correlation of Road Inputs for Left and Right Wheels

2.2.4. Four-Wheel Pavement Simulation Results

2.3. Measurement of Road Surface Roughness

2.4. Road Test Experiment

2.4.1. Experimental Equipment

2.5. Data Filtering

2.6. Model Validation

3. Semi-Active Cab Suspension Control Strategies

3.1. Constraints of Semi-Active Cab Suspension Systems

3.2. PID Semi-active Cab Suspension Control Strategy

3.2.1. PID Control Optimization Objective

3.2.2. PID Control Strategy Design

3.3. Fuzzy PID Semi-active Control Strategy for Cab Suspension System

3.4. Model Predictive Control

4. Experiments and Results

4.1. Random Road Surface

4.1.1. Class A Road at 10 m/s

4.1.2. Class B Road at 5 m/s

4.1.3. Performance Comparison

4.2. Bump response

5. Conclusions

Author Contributions

Funding

Data Availability Statement

Conflicts of Interest

Appendix A

References

MDPI Initiatives

Important Links

Subscribe

e/ ec	NB	NM	NS	ZE	PS	PM	PB
NB	PB	PB	PM	PM	PS	ZE	ZE
NM	PB	PB	PM	PS	PS	ZE	NS
NS	PM	PM	PM	PS	PS	ZE	NS
ZE	PM	PM	PS	ZE	NS	NS	NM
PS	PS	PS	ZE	NS	NS	NM	NM
PM	PS	ZE	NS	NS	NM	NM	NS
PB	ZE	ZE	NM	NM	NB	NB	NB