Next Article in Journal
Comparative Analysis of Active Bonded Piezoelectric Repair Systems for Damaged Structures under Mechanical and Thermo-Mechanical Loads
Next Article in Special Issue
An Analysis of and Improvements in the Gear Conditions of the Automated Mechanical Transmission of a Battery Electric Vehicle Considering Energy Consumption and Power Performance
Previous Article in Journal
Fixed-Time Congestion Control for a Class of Uncertain Multi-Bottleneck TCP/AWM Networks
Previous Article in Special Issue
Design of a Suspension Controller with an Adaptive Feedforward Algorithm for Ride Comfort Enhancement and Motion Sickness Mitigation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Research on Coordinated Control of Vehicle Inertial Suspension Using the Dynamic Surface Control Theory

1
School of Mechanical Science & Engineering, Huazhong University of Science and Technology, Wuhan 430074, China
2
Chinese Scholartree Ridge State Key Laboratory, China North Vehicle Research Institute, Beijing 100072, China
3
School of Mechanical Engineering, Beijing Institute of Technology, Beijing 100081, China
4
School of Automotive and Traffic Engineering, Jiangsu University, Zhenjiang 212013, China
5
State Key Laboratory of Advanced Design and Manufacturing Technology for Vehicle, Hunan University, Changsha 410082, China
6
Faculty of Vehicle and Energy Engineering, Thai Nguyen University of Technology, Thai Nguyen 250000, Vietnam
*
Author to whom correspondence should be addressed.
Actuators 2024, 13(10), 389; https://doi.org/10.3390/act13100389
Submission received: 27 August 2024 / Revised: 24 September 2024 / Accepted: 25 September 2024 / Published: 2 October 2024
(This article belongs to the Special Issue Modeling and Control for Chassis Devices in Electric Vehicles)

Abstract

:
In the process of driving, the steering, braking, and driving conditions and different road conditions affect the vibration characteristics of the vehicle in the vertical, roll, and pitch directions. These factors greatly impact the riding comfort of the vehicle. Among them, the uneven distribution of vertical load between the left and right or the front and rear suspension is one of the important factors affecting the performance indicators of the vehicle’s roll angle acceleration and pitch angle acceleration. In order to improve the ride comfort of the vehicle in vertical, roll, and pitch motion, the inerter is introduced in this paper to form a new type of suspension structure with the “spring-damping” base element, inertial suspension. It breaks away from the traditional “spring-damping” base element of the inherent suspension structure. In this paper, the mechatronic inerter is taken as the actual controlled object, and the inertial suspension structure is considered as the controlled model based on the dynamic surface control theory and the pseudo-inverse matrix principle. Thus, the coordinated control of the inertial suspension can be achieved. Under random road input, compared with passive suspension, the ride comfort performance indicators of the vehicle with inertial suspension based on dynamic surface control are significantly improved. Finally, a Hardware-in-the-Loop (HiL) test of the controller based on dynamic surface control is carried out to verify that the performance of the vehicle inertial suspension using the dynamic surface control algorithm had improved in terms of vehicle ride comfort. The error between the experimental results and the simulation results is about 8%, which verifies the real-time performance and effectiveness of the dynamic surface controller in the real controller.

1. Introduction

The automobile suspension system [1] is responsible for the transmission of force and moment between the body and the wheel, weakening the body vibration [2,3,4] caused by the ground impact, which plays a direct or indirect role in the pitch, roll, and vertical motion of the vehicle, affecting the driving safety of the car and the ride comfort of the occupant. Extreme conditions such as high-speed vehicle driving or slippery road steering have great impact on vehicle roll motion, which is prone to rollover and other accidents [5]. Frequent braking and acceleration conditions of vehicles are very unfavorable to vehicle pitching motion, which easily causes discomfort to drivers and passengers [6]. Jiang Zhidong [7] built a half car model including braking, analyzed the influence of hydraulic interconnected suspension parameters on the anti-pitch performance of the vehicle, and significantly improved the ride comfort and anti-pitch performance of the vehicle by designing ceiling and PID controllers. Therefore, it is very important to study the vehicle suspension system to improve the vehicle ride comfort.
As for vehicle coordinated control, research experts have carried extensive research on vehicle roll and pitch motion. The coordinated control of roll motion mainly encompasses the following key aspects: control of steering angle, distribution of wheel torque, control of vertical load transfer, etc., typically achieved through differential braking, steering-by-wire, active suspension or active lateral stabilizer bar, and other ways to achieve relevant functions [8]. Benyamin Anafjeh [9] designed a model predictive control based on Kalman filter and realized the active control of anti-roll bar by using a brushless motor as the actuator, which improved the driving safety and ride comfort of high-speed railway vehicles. In China, Jiang Heyao [10] designed a model predictive four-wheel hub motor controller considering ride comfort and safety. The fishhook working condition proved that this strategy effectively enhanced the anti-roll performance of the vehicle. Li Yushan [11] compared vehicle coordination control strategies using the Pareto optimal method and linear quadratic optimal control, and the results showed that the strategy using the Pareto optimal method significantly enhanced the roll stability of commercial vehicles while maintaining the path tracking accuracy. Active suspension design aims to reduce vehicle roll vibration caused by uneven vertical load distribution of left and right wheels due to uneven road surface. Zhang Liangxiu [12] studied the three-free-roll model that only considers the driver’s steering input and the seven-degrees-of-freedom model that considers both the driver’s input and uneven road surface input, as well as controlled variable damping by designing a sliding mode variable structure. The anti-roll performance of vehicles has been significantly improved. Pan Gongyu [13] proposed a vehicle roll evaluation index based on the principle of zero moment point, and designed an H∞ robust controller to adjust the stiffness of the lateral stabilizer bar. Compared with the lateral stabilizer bar controlled by passive and fuzzy PID, the lateral stabilizer bar controlled by H robust has a lower roll evaluation index, effectively enhancing the anti-roll performance of the vehicle. Li Shengqin [14] established a parameterized model of a minibus on the CarSim platform, set the lateral load transfer rate as the rollover threshold, and designed the active suspension control based on fuzzy PI and the differential braking control based on fuzzy PID. Simulation results show that the control strategy of integrated active suspension control and differential braking can improve the anti-roll performance and handling stability of the vehicle. Based on the vehicle’s 14-DOF model, Zhao Shun [15] established hierarchical control of multiple sub-systems including active suspension control, active front steering, and anti-lock braking, and the simulation shows that the anti-roll performance of the vehicle has been significantly improved.
Vehicle pitch motion is generally caused by frequent acceleration, deceleration, and uneven road input [16], and the pitch vibration resulting from uneven road input is usually suppressed by suspension control. Aiming at the coordinated control in the pitch direction, Chen Guanpeng [17] designed a hydraulic-pneumatic suspension with active wheelbase adjustment for six-wheel independently driven unmanned vehicles, and proposed a layered control strategy in which the upper PID sliding mode estimated the pitching torque, the middle-layer particle swarm optimization algorithm determined the accumulator pressure, and the lower layer used the particle swarm optimization algorithm again to achieve wheel torque distribution. Multi-scenario simulation verified the effectiveness of this layered control strategy. Liang Jinhao [18] proposed a decentralized cooperative control framework to integrate the active front steering system and the active suspension system, enhancing the vehicle lateral and vertical stability during path tracking. Shen Yujie [19] investigated a novel type of power-driven damper integrated with an inertial suspension system and designed a fractional derivative sliding mode controller to trace the optimal suspension structure. The findings indicate that the proposed control approach significantly enhances the vibration suppression and control efficiency. Igor Škrjanc [20] designed a control method based on decomposition model prediction to meet the needs of aircraft pitch angle control with fast response and small vibration, which can reduce the instantaneous vibration amplitude compared with the classical compensator control method. Donald L. Margolis [21] analyzed the performance of semi-active suspension with isolated active suspension actuator power supply and verified that the performance of semi-active suspension with isolated power supply was similar to that of full active suspension, which could effectively improve the anti-pitch performance of vehicles. Cao Jiangtao [22] proposed an improved semi-vehicle model by introducing vehicle pitch angle state, which can verify the effectiveness of linear quadratic control. Based on the linear parameter variational theory, C. Poussot Vassal [23] proposed a semi-active suspension control strategy that meets prior brake constraints, and further studies revealed the law of the influence of semi-active suspension control on the anti-pitch performance of the vehicle. Pierpaolo De Filippi [24] improved the anti-pitching performance of vehicles by adjusting the control strategy of braking and traction torque.
In this paper, based on inertial suspension structure, the NSGA-II algorithm was used to identify the inertial suspension structure with superior performance, the dynamic surface control theory was combined with the vehicle inertial suspension using the mechatronic inerter, and the vehicle coordination controller was designed based on the dynamic surface control theory with the optimal inertial suspension structure as the basic structure, which effectively reduced the root-mean-square (RMS) value of the body vibration, roll angle and pitch angle acceleration. It provides a new scheme and direction for improving vehicle ride comfort.
The content is structured as follows: Section 2 presents the vehicle inertial suspension model. In Section 3, the NSGA-II algorithm is briefly introduced, and the objective optimization function is constructed. Section 4 presents the designed coordination controller based on dynamic surface control theory. Section 5 analyzes the performance of vehicle inertial suspension based on dynamic surface control. Verification analyses using the HiL test are presented Section 6. In Section 7, the conclusions of this paper are summarized.

2. Vehicle Inertial Suspension Model

We consider the full-vehicle model depicted in Figure 1 to study the dynamic behavior of the vehicle suspension system. In this study, L4 is used as the inertial suspension structure [25].
As shown in Figure 2. The mechatronic inerter is connected to the spring and damper in series, and the load spring is connected to the mechatronic inerter, spring, and damper in parallel. The blue box shows the structure of the mechatronic inerter. The mechatronic inerter consisting of motor, ball screw and flywheel, When the mechatronic inerter is stimulated by vibration, the linear motion is converted into the rotating motion of the flywheel and motor through the ball screw, and the electromagnetic force generated by the motor is output in parallel with the inertial force.
The kinematic equation of the controlled model is
m s z ¨ s = F f l + F f r + F r l + F r r I x φ ¨ = m t ( v ˙ + u ψ ˙ ) h 2 + m s g h 1 sin ( φ ) + ( F f l F f r ) w f 2 + ( F r l F r r ) w r 2 I y θ ¨ = m t ( u ˙ v ψ ˙ ) h 3 + m s g h 1 sin ( θ ) + l r ( F r l + F r r ) l f ( F f l + F f r ) I z ψ ¨ = l f [ ( F y f l + F y f r ) cos ( δ ) + ( F x f l + F x f r ) sin ( δ ) ] l r ( F y r l + F y r r ) + w f 2 [ ( F x f r F x f l ) cos ( δ ) + ( F y f l F y f r ) sin ( δ ) ] + w r 2 ( F x r r F x r l ) + ( M z f l + M z f r + M z r l + M z r r ) m t ( v ˙ + u ψ ˙ ) = ( F x f l + F x f r ) sin ( δ ) + ( F y f l + F y f r ) cos ( δ ) + F y r l + F y r r m t ( u ˙ v ψ ˙ ) = ( F x f l + F x f r ) cos ( δ ) ( F y f l + F y f r ) sin ( δ ) + F x r l + F x r r
The vertical displacement equation for the four corners of the body is
Z s f l = z s l f θ + w f 2 φ Z s f r = z s l f θ w f 2 φ Z s r l = z s + l r θ + w r 2 φ Z s r r = z s + l r θ w r 2 φ
The equation of motion of the unsprung mass is
m u f l Z ¨ u f l = k t ( Z r f l Z u f l ) F f l m u f r Z ¨ u f r = k t ( Z r f r Z u f r ) F f r m u r l Z ¨ u r l = k t ( Z r r l Z u r l ) F r l m u r r Z ¨ u r r = k t ( Z r r r Z u r r ) F r r
The suspension force equation is
F f l = K f ( Z u f l Z s f l ) + F b f l k 1 2 w f ( φ Z u f l Z s f l 2 w f ) F f r = K f ( Z u f r Z s f r ) + F b f r k 1 2 w f ( φ Z u f r Z s f r 2 w f ) F r l = K r ( Z u r l Z s r l ) + F b r l k 2 2 w r ( φ Z u r l Z s r l 2 w r ) F r r = K r ( Z u r r Z s r r ) + F b r r k 2 2 w r ( φ Z u r r Z s r r 2 w r )
Among
F b f l = b f ( Z ¨ b f l Z ¨ s f l ) + F e f l = k f ( Z u f l Z b f l ) + c f ( Z ˙ u f l Z ˙ b f l ) F b f r = b f ( Z ¨ b f r Z ¨ s f r ) + F e f r = k f ( Z u f r Z b f r ) + c f ( Z ˙ u f r Z ˙ b f r ) F b r l = b r ( Z ¨ b r l Z ¨ s r l ) + F e r l = k r ( Z u r l Z b r l ) + c r ( Z ˙ u r l Z ˙ b r l ) F b r r = b r ( Z ¨ b r r Z ¨ s r r ) + F e r r = k r ( Z u r r Z b r r ) + c r ( Z ˙ u r r Z ˙ b r r )
The equation of the tire vertical load change is
F z f l = F f l + m u f l g + m s g l r 2 ( l r + l f ) F z f r = F f r + m u f r g + m s g l r 2 ( l r + l f ) F z r l = F r l + m u r l g + m s g l f 2 ( l r + l f ) F z r r = F r r + m u r r g + m s g l f 2 ( l r + l f )
where, m t represents the vehicle mass; m s indicates the sprung mass; m u f l , m u f r , m u r l , and m u r r represent the four wheel spring mass; I x , I y , and I z represent the moments of inertia of vehicle body roll, pitch, and yaw, respectively; u represents the lateral speed of the vehicle; v represents the lateral speed of the vehicle; ψ represents the yaw angle; θ indicates the pitch angle; φ indicates the roll angle; k t denotes the radial stiffness of the tire; z s represents the vertical displacement of the body; Z s f l , Z s f r , Z s r l , and Z s r r represent the vertical displacement of the sprung mass at the four corners of the body; Z u f l , Z u f r , Z u r l , and Z u r r represent the vertical displacement of unsprung mass; Z r f l , Z r f r , Z r r l , and Z r r r represent road excitation; k 1 and k 2 represent the roll stiffness of front and rear suspension; l f and l r represent the distance from the center of mass to the front and rear wheels; w f and w r denote the front and rear wheel base; h 1 represents the distance between the body’s center of mass and the ground; h 2 indicates the distance from the roll center to the ground; h 3 represents the distance from the pitch center to the ground; F f l , F f r , F r l , and F r r indicate the vertical load of suspension; F z f l , F z f r , F z r l , and F z r r represent the vertical supporting force of the ground facing the tire; F x f l , F x f r , F x r l , and F x r r indicate the longitudinal force of the tire; F y f l , F y f r , F y r l , and F y r r indicate tire lateral force; M z f l , M z f r , M z r l , and M z r r represent the tire righting torque; b f and b r are, respectively, the inertial coefficient of the front and rear wheel hydraulic and electric coupling inerter; F b f l , F b f r , F b r l , and F b r r represent the output force of the left front wheel, right front wheel, left rear wheel, and right rear wheel suspension using the adjustable current of the inerter; F e f l , F e f r , F e r l , and F e r r represent the force of the mechatronic inerter in the left front wheel, right front wheel, left rear wheel, and right rear wheel suspension, respectively; c f and c r are the damping coefficient of the front and rear wheel suspension, respectively; k f and k r are the secondary spring stiffness of the front and rear wheel suspension, respectively; and Z b f l , Z b f r , Z b r l , and Z b r r represent the vertical displacement of the mechatronic inerter of the left front wheel, right front wheel, left rear wheel, and right rear wheel suspension, respectively.
Aiming at the control of ride comfort, this paper takes the body vibration acceleration, pitch angle acceleration, and roll angle acceleration as the control targets, while taking into account the vertical, pitch, and roll movements of the vehicle. Therefore, the dynamics equations of vertical, pitch, and roll movements in the vehicle model are rewritten as follows:
z ¨ s = f z s ( z s , t ) + g z s ( z s , t ) u z s θ ¨ = f θ ( θ , t ) + g θ ( θ , t ) u θ φ ¨ = f φ ( φ , t ) + g φ ( φ , t ) u φ
From the above formula, f z s ( z s , t ) , f θ ( θ , t ) , and f φ ( φ , t ) can be found in Appendix A points 1, 2, and 3.
g z s ( z s , t ) = 1 / m s g θ ( θ , t ) = 1 / I y g φ ( φ , t ) = 1 / I x
u z s = F e f l + F e f r + F e r l + F e r r u θ = l f ( F e f l + F e f r ) + l r ( F e r l + F e r r ) u φ = w f 2 ( F e f l F e f r ) + w r 2 ( F e r l F e r r )
The parameters of the whole car model are from a mature passenger car, and the detailed parameters are shown in Table 1.
The single wheel pavement input model with white noise as the excitation source of random pavement input is mathematically expressed as [26]:
Z ˙ r ( t ) = 2 π f 0 Z r ( t ) + 2 π G 0 u w ( t )
where G0 is the road roughness coefficient, set to 256 × 10−6 m−3. u is the speed and w(t) is the integral white noise. Assuming the vehicle is traveling at a speed of 20 m/s, the difference in excitation between the front and rear wheels is a delay determined by the speed and the length of the front to rear wheels. The unevenness displacement of the left front wheel and the right front wheel is shown in Figure 3.

3. Optimization of Vehicle Inertial Suspension Parameters

In this paper, the NSGA-II algorithm is used to optimize inertial suspension parameters, mainly considering vehicle ride comfort. The algorithm flowchart is shown in Figure 4. The optimal individual coefficient is set at 0.3, the population size amounts to 200, the number of iterations is 100, the stopping criterion is reached after 200 generations, and the deviation of the fitness function is 10−10.
To evaluate vehicle ride comfort, the RMS values of body vertical acceleration, lateral camber acceleration, and pitch angle acceleration were selected as the optimization objectives. In the optimization process, in order to reflect the effect of the inerter, we keep the stiffness of the support spring unchanged, and other parameters to be optimized are the following:
X = [ k f , k r , c f , c r , b f , b r ]
The objective function J 1 , J 2 , J 3 is obtained:
min J 1 = B A ( X ) B A p a s s min J 2 = R A ( X ) R A p a s s min J 3 = P A ( X ) P A p a s s
In the above formula, J 1 , J 2 , and J 3 , respectively, represent the objective function of body acceleration, pitch angle acceleration, and roll angle acceleration; B A ( X ) , P A ( X ) , and R A ( X ) , respectively, represent the RMS values of body, pitch angle, and roll angle acceleration of the inertial suspension to be optimized; B A p a s s , P A p a s s , and R A p a s s , respectively, represent the RMS values of body, pitch angle, and roll angle acceleration of passive suspension.
Firstly, the performance indicators for evaluating vehicle ride comfort in vertical, pitch, and roll directions should not be worse than those of passive suspension. Then, in the design of suspension parameters, considering that the working space between the wheel and the body is limited to a certain range to ensure better grip when the wheel is in contact with the road, the working space of the designed inertial suspension and the dynamic load of the tire are smaller than those of the passive suspension. Therefore, the optimization conditions are as follows:
s t B A ( X ) < B A p a s s R A ( X ) < R A p a s s P A ( X ) < P A p a s s D T L ( X ) D T L p a s s S W S ( X ) S W S p s s s U B < X < L B
In the above formula: D T L ( X ) is the RMS value of the dynamic tire load of the left front wheel suspension; S W S ( X ) represents the RMS value of the suspension working space of the left front wheel; D T L p a s s represents the RMS value of the dynamic tire load of the left front wheel passive suspension; S W S p s s s represents the RMS value of suspension working space of the left front wheel passive suspension; U B represents the lower limit of the optimization target parameter, the values for [10,000, 10,000, 500, 500, 1, 1]; and L B represents the upper limit of the optimization target parameter, the values for [40,000, 40,000, 4000, 4000, 1000, 1000]. When the vehicle is driving in a straight line at a speed of 20 m/s on a Class C random road surface, the performance indicators of the passive suspension are shown in Table 2.
After optimization by the NSGA-II algorithm, the component parameters of four kinds of inertial suspension structures are obtained. The results of parameter optimization are shown in Table 3.

4. Design of Coordination Controller Based on Dynamic Surface Control Theory

4.1. Controller Design Procedure

In this paper, the controller actually needs the ideal displacement z d , the ideal pitch angle θ d , and the ideal roll angle φ d of the body.
First, the ideal displacement of the vertical motion of the body is defined as z d , so the position error is defined as
e z 1 = z s z d
To ensure that u z s causes e z 1 to converge to 0, the derivative of e z 1 yields
e ˙ z 1 = z ˙ s z ˙ d
We design the Lyapunov function
V 1 = 1 2 e z 1 2
Then
V ˙ 1 = e ˙ z 1 e z 1 = ( z ˙ s z ˙ d ) e z 1
Assuming an intermediate dummy variable α z , to make V ˙ 1 less than or equal to 0, we define e z 2 = z ˙ s α z , then
V ˙ 1 = ( e z 2 + α z z ˙ d ) e z 1
In the Backstepping design, the intermediate virtual variable is designed as α z = c z 1 e z 1 + z ˙ d , where c z 1 is a positive real number, which can satisfy the stability of Lyapunov function. However, differentiation explosion will occur when differentiating the intermediate virtual variable, which increases the workload of the computer and the simulation time. Therefore, dynamic surface control adopts a low-pass filter to avoid this defect.
The control law of the intermediate virtual variable in the Backstepping design is defined as z ˙ ̄ s = c z 1 e z 1 + z ˙ d , α z is output through the low-pass filter 1 τ s + 1 , and intermediate virtual variable α z satisfies
τ α ˙ z + α z = z ˙ ̄ s α z ( 0 ) = z ˙ ̄ s ( 0 )
Equation (6) can obtain α ˙ z = z ˙ ̄ s α z τ , and its filtering error is y 2 = α z z ˙ ̄ s .
The design of the Lyapunov function satisfying position, intermediate dummy variable, and filtering error is as follows:
V 2 = 1 2 e z 1 2 + 1 2 e z 2 2 + 1 2 y 2 2
Taking the derivative of this gives us V ˙ 2 :
V ˙ 2 = e z 1 e ˙ z 1 + e z 2 e ˙ z 2 + y 2 y ˙ 2
From e ˙ z 2 = z ¨ s α ˙ z = f z s ( z s , t ) + g z s ( z s , t ) u z s z ˙ ̄ s α z τ , y ˙ 2 = α ˙ z z ˙ ̄ ˙ s = y 2 τ + c z 1 e ˙ z 1 z ¨ d , The converted V ˙ 2 are presented in Appendix A point 4.
The dynamic surface control law is set as
u z s = 1 g z s ( z s , t ) ( f z s ( z s , t ) + α ˙ z c z 2 e z 2 )
where c z 2 is a positive real number, and the same can be obtained:
u θ = 1 g θ ( θ , t ) ( f θ ( θ , t ) + α ˙ θ c θ 2 e θ 2 )
u φ = 1 g φ ( φ , t ) ( f φ ( φ , t ) + α ˙ φ c φ 2 e φ 2 )
According to the control force and control moment obtained by vertical, roll, and pitch motion, the force of mechatronic inerter in four suspensions can be obtained, which can be sorted as
A f X = b f
b f = u z s u θ u φ ,   A f = 1 1 1 1 l f l f l r l r w f 2 w f 2   w r 2 w r 2 , X = F e f l F e f r F e r l F e r r
When there is a solution to Equation (25), it is equivalent to r a n k ( A f ) = r a n k ( A f , b f ) , that is, the number of linearly independent columns that need to satisfy the coefficient matrix is the same as the number of linearly independent columns of the augmented matrix. From Formula (26), it can be calculated that
r a n k ( A f ) = 3 , r a n k ( A f , b f ) = 3
Since A f is 3 × 4 dimensional matrix, that is, A f is not a square matrix, it is impossible to solve the linear equations in the form of computational matrices by finding the inverse matrix of A f . Therefore, this paper constructs a generalized inverse matrix solution for A f , then the pseudo-inverse matrix C * of A f is C * = ( A f T A f ) 1 A f T , and C * A f = I can be obtained, I is the identity matrix, and its expression is C * = p i n v ( A f ) . Thus, the force of mechatronic inerter can be obtained:
F e f l F e f r F e r l F e r r = C * u z s u θ u φ
Therefore, the schematic diagram of the dynamic surface control system designed in this paper is shown in Figure 5. The road input offers external excitation. Vertical displacement and pitch and roll angle are derived from the dynamic inertial suspension system based on the hydraulic and electric coupling inerter, and the error equation is acquired with the predefined ideal vertical displacement and pitch and roll angle. Through intermediate virtual variable and first-order filter, the dynamic surface control law of vertical, pitch, and roll motion is obtained, and the pseudo-inverse matrix is constructed to address the issue. The electromagnetic force of the four suspensions can be obtained.

4.2. Stability Proof

Then, the stability of the designed dynamic surface control law is analyzed by Lyapunov stability criterion.
If V 2 ( 0 ) p , p > 0 is taken, then all signals in the closed loop system are in a bounded closed set, that is, in a compact set.
Proof. 
If V 2 = p , V 2 = 1 2 e z 1 2 + 1 2 e z 2 2 + 1 2 y 2 2 = p , and since A is a function of e z 1 , e z 2 , y 2 , z ¨ d , then A is bounded, denoted B , then A 2 B 2 1 0 .
V ˙ 2 = e z 1 e z 2 + y 2 c z 1 e z 1 2 c z 2 e z 2 2 + y 2 y 2 τ + A           e z 1 e z 2 + e z 1 y 2 c z 1 e z 1 2 c z 2 e z 2 2 1 τ y 2 2 + A y 2           1 2 e z 1 2 + e z 2 2 + 1 2 e z 1 2 + y 2 2 c z 1 e z 1 2 c z 2 e z 2 2 1 τ y 2 2 + 1 2 y 2 2 A 2 + 1 2           = ( 1 c z 1 ) e z 1 2 + ( 1 2 c z 2 ) e z 2 2 + ( 1 2 A 2 + 1 2 1 τ ) y 2 2 + 1 2
When c z 1 1 + r , c z 2 1 2 + r , 1 τ 1 2 B + 1 2 + r ,   r > 0 , then
V ˙ 2 r e z 1 2 r e z 2 2 + 1 2 A 2 B 2 2 r y 2 2 + 1 2           = 2 r V 2 + B 2 2 B 2 A 2 B 2 2 y 2 2 + 1 2 = 2 r V 2 + A 2 B 2 1 B 2 y 2 2 2 + 1 2           2 r V 2 + 1 2
V 2 = p is known from the condition, so Equation (30) can be changed to V ˙ 2 2 r p + 1 2 .
Therefore, in order to ensure that V ˙ 2 is negative definite, that is, to ensure the stability of the control system, 2 r p + 1 2 0 is r p 4 .
That is, it has been proven that when r p 4 , all the information of the closed-loop control system is in the compact set, that is, if V 2 ( 0 ) p , then V ˙ 2 0 , and thus V 2 ( t ) p . □
Lemma 1.
For any V : [ 0 , ) R , the solution of the inequality equation  V ˙ α V + f , t t 0 0  is 
V ( t ) e α ( t t 0 ) V ( t 0 ) + t 0 t e α ( t τ ) f ( τ ) d τ , t t 0 0
where  α  is a real number.
According to the conclusion of Lemma 1, Formula (32): V ˙ 2 2 r V 2 + 1 2 is solved and its convergence is analyzed, then Formula (32) can be rewritten as:
V 2 t e 2 r t t 0 V 2 t 0 + 0.5 e 2 r t t 0 t e η 1 τ d τ                                                                                     = e 2 r t t 0 V 2 t 0 + 0.5 e 2 r t 2 r e 2 r t e 2 r t 0                                                                                     = e 2 r ( t t 0 ) V 2 ( t 0 ) + 1 4 r ( 1 e 2 r ( t t 0 ) )
Obtain
l i m t V 2 ( t ) 1 4 r
Note V 2 ( t ) converges progressively and is determined by r .
Then, when 1 τ 1 2 B + 1 2 + r , take τ 0 , then r . That is what τ is based on.
The unknown parameter can be set according to the stability condition. The value of τ is 0.01. Since r p 4 and p are related to the initial value, set to 2, then r can be 1. Because of c z 1 1 + r ,   c z 2 1 2 + r , then c z 1 ,   c z 2 take 16 and 2.5, respectively. Similarly, c θ 1 ,   c θ 2 and c φ 1 ,   c φ 2 are both 2 and 2.5. The linear driving has a greater influence on the acceleration of the vehicle body, so the gain value of c z 1 ,   c z 2 is larger.

5. Performance Analysis of Inertial Suspension Based on Dynamic Surface Control

5.1. Time Domain Performance Analysis

In order to analyze the mechanical properties of suspension, a speed of 20 m/s was selected for simulation, and the road input was set as C-class random road input. Figure 6 shows the performance response comparison of passive suspension and Inertial suspension, while Table 4 shows the RMS values comparison of performance indicators of passive suspension and inertial suspension.
As can be observed from the time domain comparison presented in Figure 6, the acceleration of the vehicle body and the acceleration of the side angle are conspicuously suppressed throughout the entire time domain. The suspension working space of the four-wheel is significantly smaller than that of the passive suspension. The acceleration of the pitch angle and the dynamic tire load of the four-wheel are marginally less than those of the passive suspension. As can be seen from Table 4, the performance of the inertial suspension based on dynamic surface control is significantly better than that of the passive suspension. Compared with the passive suspension, the RMS value of the body acceleration of the inertial suspension decreases by 23.73%. Compared with that of the passive suspension, the RMS of the performance indicators of the roll angle acceleration of the inertial suspension decreased by 15.28%. For the RMS of pitch angle acceleration, the performance of inertial suspension improved by 8.44%. In addition, the RMS of dynamic tire load of the controllable inertial suspension also decreased to a certain extent. In terms of the RMS of the suspension working space, the performance of the inertial suspension with the RMS of the suspension working space of the right front wheel suspension and the right rear wheel suspension was significantly improved by 23.21%, and the improvement effect of the left front wheel suspension and the left rear wheel suspension was better than that of the passive suspension. A comprehensive analysis of Figure 6 shows that the inertial suspension based on dynamic surface control achieves the goal of improving ride comfort, and its performance indicators are superior to that of the passive suspension.

5.2. Frequency Performance Analysis

In order to analyze the performance of the inertial suspension in the frequency domain, the time domain performance indicators diagram under random road input at 20 m/s is obtained by Fourier transform, and the power spectral density comparison diagram at 1–15 Hz is obtained, as shown in Figure 7.
As can be seen from Figure 7, compared with passive suspension, the performance of the inertial suspension based on dynamic surface control is better, especially in the low frequency band 0–3 Hz, where the peak power spectral density of the inertial suspension decreases significantly. It can be concluded that the inertial suspension based on dynamic surface control can significantly improve the vibration isolation performance of the vehicle suspension in the frequency domain and effectively improve the ride comfort of the vehicle.

6. Verification Analysis of HiL Test

In this section, the overall layout scheme of Hybrid Control Unit (HCU) HiL test is designed, the dynamic surface control algorithm is written into the physical HCU, and the HiL test is carried out to obtain the control effect of the dynamic surface control algorithm in the physical controller HCU. The test and simulation results are compared to analyze the error size of the RMS value of its performance indicators. It is proven that the dynamic surface controller designed in this paper can run normally and reliably in the physical controller HCU.
The overall scheme of vehicle inertial suspension HiL test based on dynamic surface control is shown in Figure 8.
In order to increase the persuasiveness of the inertial suspension control algorithm based on dynamic surface control and make it run in the real vehicle controller, the test condition is set to run in a straight line under the C-class random road input of 20 m/s, so as to ensure that the test results of HiL are similar to the simulation results. Figure 8 shows the comparison of simulation results and HiL test results of performance indicators of the inertial suspension based on dynamic surface control.
Figure 9 and Table 5 show the time domain curves and RMS values of body acceleration, body roll angle acceleration, body pitch angle acceleration, suspension working space, and dynamic tire load of the inertial suspension based on dynamic surface control in simulation and HiL test, respectively. It can be seen from Table 5 that in the HiL test, the RMS values of the body acceleration, roll angle acceleration, and pitch angle acceleration of the inertial suspension based on dynamic surface control are 5.76%, 4.91%, and 3.12% compared with the simulation results, which are shown in Figure 9. The performance indicators of the suspension based on the HiL test are consistent with the simulation results, and the error is controlled within 8%. In summary, the experimental results show that the controller based on dynamic surface control designed in this paper can meet the needs of vehicle controllers.

7. Results

The 14-DOF vehicle suspension model was constructed, the multi-objective function for vehicle riding comfort was formulated, and the parameters were optimized through the NSGA-II algorithm. By solving the position error equations in vertical, roll, and pitch directions, and designing the intermediate dummy variable under dynamic surface control based on the Lyapunov function, the dynamic surface control law in vertical, roll, and pitch directions is obtained. Then, by solving the pseudo-inverse matrix in vertical, roll, and pitch directions, the thrust force of mechatronic inerter is obtained. Under random road input, compared with the passive suspension, the RMS value of the body acceleration, the RMS value of the body roll angle acceleration, and the RMS value of the body pitch angle acceleration of the inertial suspension based on the dynamic surface control are reduced by 23.73%, 15.28%, and 15.19%, respectively, and ride comfort of the vehicle are effectively improved.
In the simulation of the connected controller, the feasibility and practicability of the control algorithm running in the real HCU are verified. The error between the test results and the simulation results under the random road surface is about 8%.

Author Contributions

Conceptualization, Y.Z.; methodology, F.D.; software, D.S.; validation, X.W.; formal analysis, H.W.; investigation, X.Y.; writing—original draft preparation, T.Z.; writing—review and editing, V.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China, grant number 52072157.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

1.
f z s ( z s , t ) = [ K f ( Z u f l Z s f l ) + F R f l k 1 2 w f ( φ Z u f l Z s f l 2 w f ) + K f ( Z u f r Z s f r ) + F R f r k 1 2 w f ( φ Z u f r Z s f r 2 w f ) + K r ( Z u r l Z s r l ) + F R r l k 2 2 w r ( φ Z u r l Z s r l 2 w r ) + K r ( Z u r r Z s r r ) + F R r r k 2 2 w r ( φ Z u r r Z s r r 2 w r ) ] / m s
2.
f θ ( θ , t ) = { l f [ K f ( Z u f l Z s f l ) + F R f l k 1 2 w f ( φ Z u f l Z s f l 2 w f ) + K f ( Z u f r Z s f r ) + F R f r k 1 2 w f ( φ Z u f r Z s f r 2 w f ) ] + l r [ K r ( Z u r l Z s r l ) + F R r l k 2 2 w r ( φ Z u r l Z s r l 2 w r ) + K r ( Z u r r Z s r r ) + F R r r k 2 2 w r ( φ Z u r r Z s r r 2 w r ) ] + m t ( u ˙ v ψ ˙ ) h 3 + m s g h 1 sin ( θ ) } / I y
3.
f φ ( φ , t ) = { w f 2 [ K f ( Z u f l Z s f l ) + F R f l k 1 2 w f ( φ Z u f l Z s f l 2 w f ) K f ( Z u f r Z s f r ) F R f r + k 1 2 w f ( φ Z u f r Z s f r 2 w f ) ] + w r 2 [ K r Z u r l Z s r l + F R r l k 2 2 w r φ Z u r l Z s r l 2 w r K r ( Z u r r Z s r r ) F R r r + k 2 2 w r ( φ Z u r r Z s r r 2 w r ) ] + m t ( v ˙ + u ψ ˙ ) h 2 + m s g h 1 sin ( φ ) } / I x
4.
V ˙ 2 = e z 1 e z 2 + y 2 + z ˙ ̄ s z ˙ d + e z 2 f z s z s , t + g z s z s , t u z s z ˙ ̄ s α z τ + y 2 ( y 2 τ + c z 1 e ˙ z 1 z ¨ d )               = e z 1 e z 2 + y 2 + z ˙ ̄ s z ˙ d + e z 2 f z s z s , t + g z s z s , t u z s z ˙ ̄ s α z τ + y 2 ( y 2 τ + A )
5.
A = c z 1 e ˙ z 1 z ¨ d = c z 1 ( z ˙ s z ˙ d ) z ¨ d = c z 1 ( e z 2 + α z z ˙ d ) z ¨ d = c z 1 ( e z 2 + y 2 + z ˙ ̄ s z ˙ d ) z ¨ d = c z 1 ( e z 2 + y 2 c z 1 e z 1 ) z ¨ d

References

  1. Shuping, C.; Bin, Z. Automobile Construction; China Machine Press: Beijing, China, 2007. [Google Scholar]
  2. Qin, Y.C.; Zhao, Z.; Wang, Z.F.; Li, G. Study of longitudinal–vertical dynamics for in-wheel motor-driven electric vehicles. Automot. Innov. 2021, 4, 227–237. [Google Scholar] [CrossRef]
  3. Yang, X.; Zhang, T.; Shen, Y.; Liu, Y.; Bui, V.; Qiu, D. Tradeoff analysis of the energy-harvesting vehicle suspension system employing inerter element. Energy 2024, 308, 132841. [Google Scholar] [CrossRef]
  4. Shen, Y.; Chen, A.; Du, F.; Yang, X.; Liu, Y.; Chen, L. Performance enhancements of semi-active vehicle air ISD suspension. Proc. Inst. Mech. Eng. Part D-J. Automob. Eng. 2024. [Google Scholar] [CrossRef]
  5. Zhao, D.; Yin, Y.; Ni, T.; Zhang, W. Research review on chassis integrated control of heavy duty truck. J. Yanshan Univ. 2020, 44, 189–197. [Google Scholar]
  6. Pan, G.; Ding, C.; Wang, W.; Li, A. Application of extended zero moment point in vehicle pitch control evaluation and Active control. J. Jiangsu Univ. (Nat. Sci. Ed.) 2021, 42, 85–91. [Google Scholar]
  7. Jiang, Z.; Zheng, M.; Zhang, N. Research on pitch dynamics of semi-active anti-pitch hydraulic interconnected suspension. J. Vib. Shock 2020, 39, 272–278. [Google Scholar]
  8. Xia, P. Simulation study on effect of multi-axle vehicle suspension system on roll stability. Master’s Thesis, Wuhan University of Technology: Wuhan, China, 2015. [Google Scholar]
  9. Anafjeh, B.; Moosavi, H.; Danesh, M. Active optimal roll control of railway vehicles in curved tracks using an electrically actuated anti-roll bar System. Int. J. Control. Autom. Syst. 2023, 21, 1127–1142. [Google Scholar] [CrossRef]
  10. Jiang, H.; Wang, Y.; Wu, Y.; Wang, P. Cooperative control strategy for lateral stability and roll prevention of four-wheel drive electric vehicle. J. Jilin Univ. (Eng. Technol. Ed.) 2024, 54, 540–549. [Google Scholar]
  11. Li, Y.; Chi, Y.; Ji, X.; Liu, Y.; Wu, J. Research on path tracking and anti-roll cooperative control of commercial vehicle based on pareto optimal equilibrium theory. Automot. Eng. 2019, 41, 654–661+687. [Google Scholar]
  12. Zhang, L.; Zhou, J.; Cheng, Q.; Wu, G. Estimation and active control of vehicle roll state considering road uneven input. J. Vib. Shock. 2018, 37, 133–141+154. [Google Scholar]
  13. Pan, G.; Ding, C.; Li, Y. Research on vehicle roll evaluation index and roll control based on zero moment point. Chin. J. Mech. Eng. 2023, 59, 175–187. [Google Scholar]
  14. Li, S.; Du, P.; Feng, X. Integrated control strategy of anti-rollover chassis for minibus. J. Jiangsu Univ. (Nat. Sci. Ed.) 2022, 43, 131–138. [Google Scholar]
  15. Zhao, S.; Li, Y.; Yu, Q. Vehicle stability control based on coordinated control of chassis sub-systems. J. Traffic Transp. Eng. 2015, 15, 77–85. [Google Scholar]
  16. Chen, G.; Lv, S.; Li, G.; Dai, J.; Yang, Z. Research on semi-vehicle semi-active fuzzy PI control suspension based on particle swarm optimization. J. Henan Univ. Technol. (Nat. Sci. Ed.) 2020, 39, 69–79. [Google Scholar]
  17. Chen, G.; Xu, X.; Tang, Y.; Jiang, G. A Three-Axis Unmanned Vehicle Hydraulic Drive System and Its Comprehensive Energy-Saving Control. Strategy. Patent 202211532845.6, 18 April 2023. [Google Scholar]
  18. Liang, J.; Lu, Y.; Pi, D.; Yin, G.; Zhuang, W.; Wang, F.; Feng, J.; Zhou, C. A decentralized cooperative control framework for active steering and active suspension: Multi-agent approach. IEEE Trans. Transp. Electrif. 2021, 8, 1414–1429. [Google Scholar] [CrossRef]
  19. Shen, Y.; Jia, M.; Yang, X.; Liu, Y.; Chen, L. Vibration suppression using a mechatronic PDD-ISD-combined vehicle suspension system. Int. J. Mech. Sci. 2023, 250, 108277. [Google Scholar] [CrossRef]
  20. Škrjanc, I. A decomposed-model predictive functional control approach to air-vehicle pitch-angle control. J. Intell. Robot. Syst. Theory Appl. 2007, 48, 115–127. [Google Scholar] [CrossRef]
  21. Margolis, D.L. Semi-active heave and pitch control for ground vehicles. Veh. Syst. Dyn. 1982, 11, 31–42. [Google Scholar] [CrossRef]
  22. Cao, J.T.; Liu, H.H.; Li, P.; Brown, D.J.; Dimirovski, G. A study of electric vehicle suspension control system based on an improved half-vehicle model. Int. J. Autom. Comput. 2007, 4, 236–242. [Google Scholar] [CrossRef]
  23. Poussot, V.C.; Sename, O.; Dugard, L.; Gáspár, P.; Szabó, Z.; Bokor, J. A new semi-active suspension control strategy through LPV technique. Control. Eng. Pract. 2008, 16, 1519–1534. [Google Scholar] [CrossRef]
  24. De, F.P.; Tanelli, M.; Corno, M.; Savaresi, S.M.; Santucci, M.D. Electronic stability control for powered two-wheelers. IEEE Trans. Control. Syst. Technol. 2014, 22, 265–272. [Google Scholar]
  25. Liu, Y.; Shi, D.; Du, F.; Yang, X.; Zhu, K. Topological optimization of vehicle ISD suspension under steering braking condition. World Electr. Veh. J. 2023, 14, 297. [Google Scholar] [CrossRef]
  26. Zhang, Y.; Zhang, J. Numerical simulation of stochastic road process using white noise filtration. Mech. Syst. Signal Process. 2005, 20, 363–372. [Google Scholar]
Figure 1. Vehicle control model.
Figure 1. Vehicle control model.
Actuators 13 00389 g001
Figure 2. Structure of vehicle inertial suspension.
Figure 2. Structure of vehicle inertial suspension.
Actuators 13 00389 g002
Figure 3. Pavement excitation displacement of the left front wheel and the right front wheel.
Figure 3. Pavement excitation displacement of the left front wheel and the right front wheel.
Actuators 13 00389 g003
Figure 4. Flowchart of the NSGA-II algorithm.
Figure 4. Flowchart of the NSGA-II algorithm.
Actuators 13 00389 g004
Figure 5. Schematic diagram of dynamic surface control system of inertial suspension based on hydraulic and electric coupling inerter.
Figure 5. Schematic diagram of dynamic surface control system of inertial suspension based on hydraulic and electric coupling inerter.
Actuators 13 00389 g005
Figure 6. Comparison of performance response of inertial suspension. (a) Body acceleration. (b) Roll angle acceleration. (c) Pitch angle acceleration. (d) Left front wheel suspension working space. (e) Left front dynamic tire load. (f) Right front wheel suspension working space. (g) Right front dynamic tire load. (h) Left rear wheel suspension working space. (i) Left rear dynamic tire load. (j) Right rear wheel suspension working space. (k) Right rear dynamic tire load.
Figure 6. Comparison of performance response of inertial suspension. (a) Body acceleration. (b) Roll angle acceleration. (c) Pitch angle acceleration. (d) Left front wheel suspension working space. (e) Left front dynamic tire load. (f) Right front wheel suspension working space. (g) Right front dynamic tire load. (h) Left rear wheel suspension working space. (i) Left rear dynamic tire load. (j) Right rear wheel suspension working space. (k) Right rear dynamic tire load.
Actuators 13 00389 g006aActuators 13 00389 g006bActuators 13 00389 g006c
Figure 7. Power spectral density analysis of inertial suspension performance parameters. (a) Body acceleration. (b) Roll angle acceleration. (c) Pitch angle acceleration. (d) Left front wheel suspension working space. (e) Left front dynamic tire load. (f) Right front wheel suspension working space. (g) Right front dynamic tire load. (h) Left rear wheel suspension working space. (i) Left rear dynamic tire load. (j) Right rear wheel suspension working space. (k) Right rear dynamic tire load.
Figure 7. Power spectral density analysis of inertial suspension performance parameters. (a) Body acceleration. (b) Roll angle acceleration. (c) Pitch angle acceleration. (d) Left front wheel suspension working space. (e) Left front dynamic tire load. (f) Right front wheel suspension working space. (g) Right front dynamic tire load. (h) Left rear wheel suspension working space. (i) Left rear dynamic tire load. (j) Right rear wheel suspension working space. (k) Right rear dynamic tire load.
Actuators 13 00389 g007aActuators 13 00389 g007bActuators 13 00389 g007c
Figure 8. Design of HiL test scheme of vehicle inertial suspension based on dynamic surface control theory.
Figure 8. Design of HiL test scheme of vehicle inertial suspension based on dynamic surface control theory.
Actuators 13 00389 g008
Figure 9. HiL and simulation test results. (a) Body acceleration. (b) Roll angle acceleration. (c) Pitch angle acceleration. (d) Left front wheel suspension working space. (e) Left front dynamic tire load. (f) Right front wheel suspension working space. (g) Right front dynamic tire load. (h) Left rear wheel suspension working space. (i) Left rear dynamic tire load. (j) Right rear wheel suspension working space. (k) Right rear dynamic tire load.
Figure 9. HiL and simulation test results. (a) Body acceleration. (b) Roll angle acceleration. (c) Pitch angle acceleration. (d) Left front wheel suspension working space. (e) Left front dynamic tire load. (f) Right front wheel suspension working space. (g) Right front dynamic tire load. (h) Left rear wheel suspension working space. (i) Left rear dynamic tire load. (j) Right rear wheel suspension working space. (k) Right rear dynamic tire load.
Actuators 13 00389 g009aActuators 13 00389 g009b
Table 1. Vehicle model parameters.
Table 1. Vehicle model parameters.
ParameterValue
Vehicle mass mt/(kg)1659
Sprung mass ms/(kg)1410
Front wheel unsprung mass mufl, mufr/(kg)26.5
Rear wheel unsprung mass murl, murr/(kg)24.4
Front wheel base wf/(m)1.574
Rear wheel base wr/(m)1.593
Distance from front axle to center of mass lf/(m)1.278
Distance from rear axis to center of mass lr/(m)1.430
Height of center of mass h1/(m)0.50
Roll height h2/(m)0.40
Distance from center of roll to center of mass h3/(m)0.25
Body roll moment of inertia Ix/(kg·m2)925
Body pitch moment of inertia Iy/(kg·m2)2577
Body yaw moment of inertia Iz/(kg·m2)2603
Wheel inertia Iw/(kg·m2)0.99
Wheel radius Rw/(m)0.345
Front suspension roll stiffness k1/(N·m·rad−1)47,298
Rear suspension roll stiffness k2/(N·m·rad−1)37,311
Tire equivalent spring stiffness kt/(N/m)192,000
Front suspension spring stiffness of the original model Kf/(N·m−1)25,000
Rear suspension spring stiffness of the original model Kr/(N·m−1)22,000
Front suspension spring stiffness kf/(N·m−1)9325
Rear suspension spring stiffness kr/(N·m−1)7897
Front suspension damping coefficient cf/(N·s·m−1)1579
Rear suspension damping coefficient cr/(N·s·m−1)1407
Front suspension inertial coefficient bf/(kg)330
Rear suspension inertial coefficient br/(kg)265
Table 2. Performance indicators of passive suspension.
Table 2. Performance indicators of passive suspension.
Performance IndicatorsRMS Value
Body acceleration/(m·s−2)1.3644
Roll angle acceleration/(rad·s−2)0.8326
Pitch angle acceleration/(rad·s−2)0.7094
Left front wheel suspension working space/(m)0.0225
Left front dynamic tire load/(kN)1.1633
Right front wheel suspension working space/(m)0.0151
Right front wheel dynamic tire load/(kN)0.8796
Left rear wheel suspension working space/(m)0.0180
Left rear wheel dynamic tire load/(kN)1.0851
Right rear wheel suspension working space/(m)0.0157
Right rear wheel dynamic tire load/(kN)0.8787
Table 3. Parameter optimization results.
Table 3. Parameter optimization results.
Optimization ParameterValue
Front suspension spring stiffness/(N·m−1)9325
Rear suspension spring stiffness/(N·m−1)7897
Front suspension damping coefficient/(N·s·m−1)1579
Rear suspension damping coefficient/(N·s·m−1)1407
Front suspension inertial coefficient/(kg)330
Rear suspension inertial coefficient/(kg)265
Table 4. Performance comparison table of inertial suspension.
Table 4. Performance comparison table of inertial suspension.
The RMS ValuesPassive SuspensionInertial SuspensionPerformance Improvement
Body acceleration/(m·s−2)1.36441.040623.73%
Roll angle acceleration/(rad·s−2)0.83260.705415.28%
Pitch angle acceleration/(rad·s−2)0.70940.64958.44%
Left front wheel suspension working space/(m)0.02250.019115.27%
Left front dynamic tire load/(kN)1.16331.09995.45%
Right front wheel suspension working space/(m)0.01510.011623.21%
Right front dynamic tire load/(kN)0.87960.82725.95%
Left rear wheel suspension working space/(m)0.01800.015812.05%
Left rear dynamic tire load/(kN)1.08511.04693.52%
Right rear wheel suspension working space/(m)0.01570.011228.57%
Right rear dynamic tire load/(kN)0.87870.81667.06%
Table 5. RMS and errors of HiL and simulation test performance indicators.
Table 5. RMS and errors of HiL and simulation test performance indicators.
Performance IndicatorsSimulationHiLError
Body acceleration/(m·s−2)1.04061.10065.76%
Roll angle acceleration/(rad·s−2)0.70540.74014.91%
Pitch angle acceleration/(rad·s−2)0.64950.66973.12%
Left front wheel suspension working space/(m)0.01910.02036.35%
Left front dynamic tire load/(kN)1.09991.15505.01%
Right front wheel suspension working space/(m)0.01160.01247.42%
Right front wheel dynamic tire load/(kN)0.82720.86885.03%
Left rear suspension working space/(m)0.01580.01644.05%
Left rear dynamic tire load/(kN)1.04691.09104.21%
Right rear wheel suspension working space/(m)0.01120.01174.13%
Right rear dynamic tire load/(kN)0.81660.84964.03%
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.

Share and Cite

MDPI and ACS Style

Zhao, Y.; Du, F.; Wang, H.; Wang, X.; Yang, X.; Shi, D.; Bui, V.; Zhang, T. Research on Coordinated Control of Vehicle Inertial Suspension Using the Dynamic Surface Control Theory. Actuators 2024, 13, 389. https://doi.org/10.3390/act13100389

AMA Style

Zhao Y, Du F, Wang H, Wang X, Yang X, Shi D, Bui V, Zhang T. Research on Coordinated Control of Vehicle Inertial Suspension Using the Dynamic Surface Control Theory. Actuators. 2024; 13(10):389. https://doi.org/10.3390/act13100389

Chicago/Turabian Style

Zhao, Yanhui, Fu Du, Hujiang Wang, Xuelin Wang, Xiaofeng Yang, Dongyin Shi, Vancuong Bui, and Tianyi Zhang. 2024. "Research on Coordinated Control of Vehicle Inertial Suspension Using the Dynamic Surface Control Theory" Actuators 13, no. 10: 389. https://doi.org/10.3390/act13100389

APA Style

Zhao, Y., Du, F., Wang, H., Wang, X., Yang, X., Shi, D., Bui, V., & Zhang, T. (2024). Research on Coordinated Control of Vehicle Inertial Suspension Using the Dynamic Surface Control Theory. Actuators, 13(10), 389. https://doi.org/10.3390/act13100389

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop