Next Article in Journal
A Machine Learning App for Monitoring Physical Therapy at Home
Previous Article in Journal
Quasi-Zero Stiffness Vibration Sensing and Energy Harvesting Integration Based on Buckled Piezoelectric Euler Beam
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Neural Network-Based Adaptive Height Tracking Control of Active Air Suspension System with Magnetorheological Fluid Damper Subject to Uncertain Mass and Input Delay

1
School of Mechanical and Electrical Engineering, Guizhou Normal University, Guiyang 550001, China
2
School of Vehicle and Mobility, Tsinghua University, Beijing 10084, China
3
Automotive Engineering Research Institute, Jiangsu University, Zhenjiang 212013, China
4
Research Institute of Highway, Ministry of Transport, Beijing 100088, China
*
Author to whom correspondence should be addressed.
Submission received: 2 November 2023 / Revised: 14 December 2023 / Accepted: 20 December 2023 / Published: 27 December 2023
(This article belongs to the Section Vehicular Sensing)

Abstract

:
In this paper, we present a novel robust adaptive neural network-based control framework to address the ride height tracking control problem of active air suspension systems with magnetorheological fluid damper (MRD-AAS) subject to uncertain mass and time-varying input delay. First, a radial basis function neural network (RBFNN) approximator is designed to compensate for unmodeled dynamics of the MRD. Then, a projector-based estimator is developed to estimate uncertain parameter variation (sprung mass). Additionally, to deal with the effect of input delay, a time-delay compensator is integrated in the adaptive control law to enhance the transient response of MRD-AAS system. By introducing a Lyapunov–Krasovskii (LK) functional, both ride height tracking and estimator errors can robustly converge towards the neighborhood of the desired values, achieving uniform ultimate boundness. Finally, comparative simulation results based on a dynamic co-simulator built in AMESim 2021.2 and Matlab/Simulink 2019(b) are given to illustrate the validity of the proposed control framework, showing its effectiveness to operate ride height regulation with MRD-AAS systems accurately and reliably under random road excitations.

1. Introduction

The growing need to isolate vibrations caused by uneven road surfaces has led to more concentration on the vehicle suspension system, which is considered a vital component of the chassis in automotive manufacturing companies, such as Li Auto and Tesla. The electronically controlled active air suspension (AAS) systems adjust the ride height and posture movements by inflating and deflating the air spring so as to protect the vehicle body on bumpy roads and reduce air drag. In connection with the magnetorheological fluid damper, the AAS systems are more effective in enhancing ride comfort and vehicle maneuverability during complex driving conditions as compared with conventional passive and semi-active suspensions. The active air suspension systems with magnetorheological fluid damper (MRD-AAS) provide an alternative solution to meet the stringent requirements of chassis functions.
Practically, there exist some challenges in maintaining the ride height control under random road excitations robustly and accurately with MRD-AAS systems, since it is difficult to establish a mathematical model to depict the highly nonlinear dynamic characteristics of the air spring and magnetorheological fluid damper (MRD) simultaneously, especially considering the variable stiffness and hysteresis of the adjustable air spring and MRD during the ride height regulation process. Moreover, the uncertain mass caused by the omnipresent changes in payload [1,2] and the effect of input delay induced by pneumatic and hydraulic actuators [3] impose complexities in control system design. These nonlinearities, parametric uncertainties and input delay could severely worsen the ride height tracking control performance or even lead to instability if they are not well addressed.
Presently, a series of models have been established to describe the dynamic characteristics of AAS systems for ride height control. For instance, the authors of [4] analyzed the relationship between spring force and air spring volume and built a mapping model to track the spring force using a back propagation neural network method. However, the modeling accuracy depended on a large amount of experiment tests. The authors of [5] established a nonlinear model of an AAS system on the basis of thermodynamics to represent the dynamic characteristics of the air spring, which are classified as unmodeled dynamics and parametric uncertainties. However, the nonlinear and hysteretic characteristics of the MRD are not fully considered, which would degrade ride height control performance. A number of MRD dynamic models have been constructed to describe the nonlinear and hysteretic characteristics under random road excitations, such as the normalized phenomenological model [6], the hyperbolic tangent model [7], the restructured phenomenological model [8], and the adjustable sigmoid model [9]. However, these models rely largely on the shape functions and a series of parameter identification with optimization techniques. However, by considering the omnipresent parameter variations due to the uncertain sprung mass, it is difficult to maintain the accuracy of the optimized phenomenological model. Therefore, it is important to design a real-time approximator of the phenomenological MRD force to ensure effective ride height control with the MRD-AAS system. Due to the good generalization ability and simplified form, an RBFNN will be considered to design the damper force approximator to make up for the defect of the phenomenological damper model [6,7,8,9].
Despite the challenge caused by the MRD-AAS system modeling for depicting the intrinsic strong nonlinearities and hysteretic characteristics, taking into account the designing control methods for dealing with these characteristics, sensitive parameter variations (such as sprung mass) and actuator input delay are also necessary during the ride height control procedure. In order to handle the highly nonlinear and hysteretic characteristics of the MRD-AAS system, numerous control methods have been proposed, such as sliding mode control, backstepping control, and fuzzy logic and neural network-based control techniques. The authors of [10] proposed a compensated backstepping controller combining traditional backstepping techniques with an adaptive radial basis function. In order to estimate unmodeled dynamic characteristics, fuzzy logic systems and neural networks are employed as the approximators for improving tracking control performance [11,12,13]. However, to our best knowledge, the number of required rules or neurons remains uncertain for making the approximation error stay within a small specific range. In particular, the estimation errors of the uncertain sprung mass were not proofed to converge to zero, theoretically.
Moreover, the MRD-AAS system is composed of the hydraulic and pneumatic actuators, which induce input delay during the ride height controller design procedure. It is also a challenge caused by designing controllers for handling the effect of the input delay, avoiding system instability. Although the delay may be quite short, it can restrict the performance of the controller. To conquer this problem, a Lyapunov–Krasovskii functional is introduced to handle the effect of unknown sufficiently slow time-varying input time lags. The authors of [14] used linear matrix inequalities (LMIs) to solve for the time delay problem of magnetorheological semi-active suspensions and obtained conclusions such as a maximum critical delay of 32.4 (ms) for MRD. The paper [15] investigated the nonlinear dynamics of a constant coefficient time lag vibration control system in detail, and obtained a series of important criteria such as the ability of the system to remain stable in a smaller interval of time delay feedback. The authors of [16] used a locally optimal hybrid PSO/LMI algorithm for perturbations and time delays and tuned the continuous saturation controller by adjusting the feedback gain, but they did not consider the effect of sprung mass variations from the controller. Thus, inspired by the above, a predictor-based compensation scheme is used in this work to address the effects of input delay and uncertainty.
Inspired by the above discussion and considerations, this paper presents a solution to address the height tracking control problem for nonlinear MRD-AAS systems with dynamic characteristics, uncertain parameter variation, and time-varying input delay. The novelties and contributions of this paper are summarized as follows:
1.
A robust adaptive ride height tracking control framework is proposed, where a MRD force approximator is designed to improve the accuracy of MRD-AAS system modeling, achieving uniform ultimate boundedness in the presence of uncertain sprung mass and time-varying input delay.
2.
A projector-based nonlinear estimator is developed to estimate the uncertain sprung mass for enhancing the adaptive performance of the proposed control system.
3.
A time-delay compensator is introduced to deal with the time-varying input delay induced by the pneumatic and hydraulic actuators for improving the robustness of the proposed control system.
With respect to the existing literature, the main merits of the proposed adaptive robust ride height control strategy are as follows. Different from the phenomenological damper model presented in [6,7,8,9], the proposed control system is designed based on a nonlinear mathematical model with the MRD force approximator so as to improve the modeling accuracy of the MRD-AAS system. Compared with the mass estimation method presented in [11,12,13], the developed projector-based mass estimator can ensure that the estimation error is bounded by a certain value, achieving adaptive performance. To improve the anti-interference performance, a time-delay compensator is integrated in the proposed control strategy that allows us to compensate for time-varying input delay. The remainder of this paper is structured as follows.
The organization of this paper is as follows: In Section 2, the notation used throughout this paper is presented. Section 3 involves the dynamic modeling of the MRD-AAS. Section 4 describes the design method for the MRD-AAS system. Section 5 illustrates the effectiveness of proposed control strategy in co-simulation. Finally, Section 6 summarizes the paper.

2. Notation

R n denotes the n-dimensional Euclidean space and is used throughout the paper. A function f is defined as R R n is of class C n if the derivation exists and f , f , , f n are continuous. For your reference, Table 1 lists the main symbols and descriptions of the model, controller, and estimator parameters as follows.

3. Problem Formulation

The objective of this section is to formulate the problem of ride height tracking control with application to the MRD-AAS suspension system. As displayed in Figure 1, there are several components among the sprung mass and unsprung mass of a quarter vehicle MRD-AAS suspension, such as an adjustable air spring, an MRD, and other accessories. We start by describing the modeling of a quarter vehicle with MRD-AAS by considering the variable passenger or/and payload, time lags induced by the pneumatic and hydraulic actuators, and disturbance. Then, the issue of maintaining ride height tracking control is formulated by taking into account the uncertain sprung mass, time-varying input delay, and disturbance.

3.1. MRD-AAS System Modeling

For addressing the problem of the ride height tracking control, based on the modified Binghan mathematical model for magnetorheological dampers, a mathematical model of a quarter vehicle with MRD-AAS is established. The schematic diagram of this quarter vehicle with MRD-AAS is displayed in Figure 1, and based on Newton’s second law, the dynamical equations of the masses are expressed by
m s z ¨ s = F a s F d F g m u z ¨ u = F a s + F d F k D t ,
where F a s denotes the air spring force; F d represents MRD damper force; F g is the gravitational force; z s and z u are the sprung mass and unsprung mass displacement of a quarter vehicle; and m s and m u are the sprung mass and unsprung mass (the wheel assembly). F t and D t are the elasticity force and damping force of the tire, respectively. Forces produced by the adjustable air spring and the tire are given as
F a s = A ¯ a s ( p a s p a t m ) F k = k t ( z u z r ) D t = c t ( z ˙ u z ˙ r ) ,
where p a t m denotes the atmospheric pressure, s p = 10 5 , and A ¯ a s = A a s s p is the effective area of the adjustable air spring. z r denotes the random road excitation, and k t and c t are the stiffness and damping coefficient of the tire, respectively. According to the performance tests of the magnetorheological hysteresis suspension presented in [17,18], a hyperbolic tangent model is employed to effectively describe the nonlinear hysteretic characteristics of the MRD current force, expressed by
F d = f t a n h ( α z ˙ s + β s i g n ( z s ) ) + c z ˙ s ,
where
f = f 1 I + f 2 α = α 1 I + α 2 β = β 1 I + β 2 c = c 1 I + c 2 ,
where f , f 1 , f 2 denote the hysteresis loop scaling factors; α , α 1 , α 2 are the hysteresis loop slope scaling factors, β , β 1 , β 2 are the half-width factors of the hysteresis loop, I is the current, and c , c 1 , c 2 are the MRD damping coefficients.
Due to the inherent nonlinearities and uncertainties in the MRD-AAS system, the modeling lumped error, time-varying input delay, and disturbance should be considered in the employed model for ride height control. Following [5], the MRD-AAS is expressed as
x ˙ 1 = x 2 x ˙ 2 = m s 1 A ¯ a s ( x 3 p a t m ) f t a n h ( α x 2 + β s i g n ( x 1 ) ) c x 2 F δ g x ˙ 3 = u τ x 2 x 3 A a s γ v a s 1 + ( γ 1 ) ( s p v a s ) 1 Q ˙ + d n ( t ) ,
where x 1 = z s , x 2 = z ˙ s , x 3 = p a s / s p , v a s is air volume, Q ˙ is the heat transfer rate, and γ denotes the specific heat ratio. u τ denotes the delayed control input, d n ( t ) is external time-varying disturbance, satisfying d n ( t ) d ¯ n with d ¯ n denoting a known positive constant, and F δ represents the concentrated modeling error of the MRD. Details of the nonlinear active air suspension system modeling can be found in [5].
Remark 1. 
By considering the practical operating scenarios of ride height tracking control with MRD-AAS systems, the vertical displacement of sprung mass x 1 , velocity x 2 , and x 3 could be measured by the configured signal acquisition devices, e.g., the vehicle-equipped displacement and pressure sensors.
Assumption 1. 
Under the same vertical excitation, the air spring force F a s , the MRD force F d , and their derivatives are bounded and satisfy
| F a s | F ¯ a s , | F d | F ¯ d , | F ˙ a s | F ¯ a t , | F ˙ d | F ¯ d t ,
where F ¯ a s , F ¯ d , F ¯ a t , and F ¯ d t are known positive constants [11].
Assumption 2. 
By considering the omnipresent changes of passenger or/and payload, the sprung mass m s is uncertain, constant and satisfied
m m i n m s m m a x ,
where m m i n and m m a x are known minimum and maximum constants of the vehicle sprung mass.
Assumption 3. 
The time delay τ is unknown, positive, time-varying, and satisfies τ τ ¯ , where τ ¯ is a known positive constant. There exists a sufficiently accurate constant τ k , satisfying
| τ τ k | τ u , | u ( t τ ) u ( t τ k ) | u e ,
where τ u and u e are known constants, and τ k τ m a x . The MRD-AAS system (5) stays finite during the period of [ t 0 , t 0 + τ ¯ ] .

3.2. RBFNN Approximator

Since the normalized phenomenological models, such as the hyperbolic tangent model, are not accurate enough to capture the nonlinear and hysteretic characteristics of MRD (3), a RBFNN approximator is designed to approximate F δ online, which caused by the unmodeled dynamics of the MRD. Following [19], the hidden Gaussian function h j for the neural net j is adopted as
h j = exp ( ( χ ϕ j ) 2 / 2 b j 2 ) , j = 1 , 2 , . . . , l ,
where χ = [ χ 1 , χ 2 , , χ 3 ] T is the input vector, ϕ j = [ ϕ 1 , ϕ 2 , , ϕ n ] is the center vector of neural net j and b j is the width of the Gaussian function that are chosen based on a range of input values [20]. Then, F δ can be expressed as
F δ ( t ) = W * T h + ϵ ,
where W * denotes the desired weights of neural net j, h = [ h 1 , h 2 , , h l ] T denotes the hidden layer function, and ϵ is the approximation error. Therefore, approximate smooth function F ^ δ can be estimated by the RBFNN as follows
F ^ δ ( t ) = W ^ T h ,
where W ^ T the weights of the RBFNN [21]. To evaluate the performance of the RBFNN approximator, a performance index function is adopted as
E ( t ) = F d ( t ) + F ^ δ ( t ) F m ( t ) 2 / 2 ,
where F m ( t ) is the training data obtained from experimental tests of the MRD under the current conditions of I = 0 3 ( A ) . By using the gradient descent learning algorithm for effectively adjusting the weights Δ w j ( t ) , centers Δ ϕ j , i ( t ) , and widths Δ b j ( t ) numerically to minimize the squared error function E [20], the following equations are given as
Δ w j ( t ) = E / w j = ( F d ( t ) + F δ ( t ) F m ( t ) ) h j Δ ϕ j , i ( t ) = E / Δ ϕ j , i = ( F d ( t ) + F δ ( t ) F m ( t ) ) w j h j ( χ i ϕ j i ) / b j 2 ) Δ b j ( t ) = E / b j = ( F d ( t ) + F δ ( t ) F m ( t ) ) w j h j ( ( χ ϕ j ) / b j 3 ) .
The momentum factor ı domain is 0 ı 1 , the updated parameters w j , ϕ j i , and b j are obtained as
w j ( t ) = w j ( t 1 ) + Δ w j + ı ( w j ( t 1 ) w j ( t 2 ) ) ϕ j i ( t ) = ϕ j i ( t 1 ) + Δ ϕ j i + ı ( ϕ j i ( t 1 ) ϕ j i ( t 2 ) ) b j ( t ) = b j ( t 1 ) + Δ b j + ı ( b j ( t 1 ) b j ( t 2 ) ) .
To verify the effectiveness of the constructed MRD model, a MRD test bench is fabricated as shown in Figure 2a, which consists of an upper fixture, DC power supply, hydraulic exciter, and MRD. As Figure 2b shows, the measured external characteristics of the employed MRD at the current from 0 to 3 ( A ) are collected and given. Moreover, Figure 3a–h show the comparison of a hyperbolic tangent model with and without the designed RBFNN approximator. It can be found that the designed RBFNN approximator could compensate for the concentrated modeling errors and keep the approximated error in a bound of 20 ( N ) . Therefore, the constructed MRD model combined the hyperbolic tangent model with the designed RBFNN approximator, which was then applied in the following control system design.

3.3. Problem Statement

Taking into account nonlinearities, parameter variations, and time-varying input delays, we tackle the ride height tracking control problem of MRD-AAS systems, stating that the desired height x d should be located into a curve of class C 3 or above, with limited temporal derivatives. An arbitrary small neighborhood of the desired height x d may be stabilized at all times by designing a control input u τ in the presence of unmodeled dynamics, uncertain sprung mass, time-varying input delays, and disturbance such that the ride height z s can always be stabilized in an arbitrarily small neighborhood of its desired height x d .

4. Adaptive Controller Design

As mentioned in the Introduction, the uncertain sprung mass, time-varying input delays, and highly nonlinear dynamic characteristics of the MRD-AAS system impose difficulties on the height tracking control system design. Based on the established quarter vehicle model with the MRD-AAS (5), a delay-free ride height control law u τ is synthesized for guaranteeing ultimate global uniform ultimate boundedness. Figure 4 shows the control block diagram of the proposed control strategy. The on-board sensors can measure the states of the MRD-AAS system x 1 , x 2 , and x 3 , as stated in Remark 1. To compensate for the hyperbolic tangent modeling error of MRD, a RBFNN approximator is developed with an adaptive gradient descent learning algorithm. A projector-based estimator is designed to estimate the uncertain sprung mass m s mostly caused by the payload. We start the controller design by defining the height tracking error, given as
e 1 = x 1 x d ,
and an auxiliary error term e n as
e n = 0 t e 1 ( θ ) d θ .
Our first Lyapunov candidate function is defined as
V 1 = 1 2 e 1 2 + 1 2 ξ 2 ,
where ξ = e n + α e 1 , whose time derivative yields
V ˙ 1 = W 1 ( ξ , η ) + η ( x ˙ 1 x ˙ d + k 1 η + ξ ) ,
where η = α ξ + e 1 , W 1 ( ξ , η ) = k 1 η 2 + α ξ 2 , and α , k 1 are positive constants.
Following the backstepping technique procedure, we define a new error e 2 as
e 2 = x ˙ 1 x ˙ d + k 1 η + ξ .
Constructing a new Lyapunov function candidate by incorporating e 2 , we define
V 2 = V 1 + 1 2 e 2 2 .
Computing the time derivative of V 2 , we have
V ˙ 2 = W 2 ( ξ , η , e 2 ) + e 2 k 2 e 2 + η + ϖ ( F a s F d F ^ δ ) g x ¨ d + k 1 η ˙ + ξ ˙ e 2 ϖ F ˜ δ ,
where W 2 ( ξ , η , e 2 ) = k 1 η 2 + α ξ 2 + k 2 e 2 2 , ϖ = 1 / m s , and k 2 is a positive control gain, F ˜ δ is the estimation error of MRD force, given as
F ˜ δ = F δ F ^ δ = W * T h + ϵ W ^ T h = W ˜ T h + ϵ ,
where W ˜ T = W ^ T W * T . Notice that ϖ [ 1 / m m a x , 1 / m m i n ] is unknown. Then, we introduce the estimated ϖ ^ , which is bonded by make the following projection function, given as
p ( ϖ ^ ) = ϖ m a x + ϱ 1 exp ( ϖ m a x ϖ ^ ϱ ) , ϖ ^ ϖ m a x ϖ ^ , e l s e ϖ m i n + ϱ 1 e x p ( ϖ ^ ϖ m i n ϱ ) , ϖ ^ ϖ m i n
where ϖ m i n , and ϱ is an arbitrarily small positive number. Then, (21) is rewritten as
V ˙ 2 = W 2 ( ξ , η , e 2 ) + e 2 k 2 e 2 + η + p ( ϖ ^ ) ( F a s F d F ^ δ ) x ¨ d g + k 1 η ˙ + ξ ˙ + z a e 2 z a + e 2 ϖ p ( ϖ ^ ) ( F a s F d F ^ δ ) e 2 ϖ ( W ˜ T h + ϵ ) ,
and z a is defined as
z a = t τ k t u ( θ ) d θ .
The term ϖ p ( ϖ ^ ) and ( W ˜ T h + ϵ ) in the last two terms of (24) result from the uncertain sprung mass and the concentrated modeling error of the MRD force. A new Lyapunov candidate function is then defined to obtain the estimation law for ϖ ^ ˙ and W ^ ˙ , given as
V 2 b = V 2 + 1 γ w W ˜ T W ˜ + 1 λ 0 ϖ ˜ ( p ( ϖ + τ ) ϖ ) d τ ,
where ϖ ˜ = ϖ ^ ϖ denotes the estimation error for sprung mass. Computing the time derivative of V 2 b , we obtain
V ˙ 2 b = W 2 ( ξ , η , e 2 ) + e 2 k 2 e 2 + η + p ( ϖ ^ ) ( F a s F d F ^ δ ) x ¨ d g + k 1 η ˙ + ξ ˙ + z a e 2 z a + ϖ p ( ϖ ^ ) e 2 ( F a s F d F ^ δ ) 1 λ ϖ ^ ˙ e 2 ϖ ϵ W ˜ T e 2 ϖ h + 1 γ w W ^ ˙ .
We choose the updating law for ϖ ^ ˙ and W ^ ˙ as
ϖ ^ ˙ = λ e 2 ( F a s F d F ^ δ ) ,
W ^ ˙ = γ w e 2 m m i n h .
Substituting (28) and (29) into (24), we obtain
V ˙ 2 b = W 2 ( ξ , η , e 2 ) + e 2 k 2 e 2 + η + p ( ϖ ^ ) ( F a s F d F ^ δ ) x ¨ d g + k 1 η ˙ + ξ ˙ + z a e 2 z a e 2 ϖ ϵ W ˜ T h e 2 ( m m i n ϖ ) W 2 ( ξ , η , e 2 ) + e 2 k 2 e 2 + η + p ( ϖ ^ ) ( F a s F d F ^ δ ) x ¨ d g + k 1 η ˙ + ξ ˙ + z a e 2 z a e 2 ϵ W ˜ T h e 2 ( m m i n ϖ ) W 2 ( ξ , η , e 2 ) + e 2 k 2 e 2 + η + p ( ϖ ^ ) ( F a s F d F ^ δ ) x ¨ d g + k 1 η ˙ + ξ ˙ + z a e 2 z a e 2 ϵ .
Remark 2. 
The introduced functions 1 γ w W ˜ T W ˜ and 1 λ 0 ϖ ˜ ( p ( ϖ + τ ) w ) d τ in (26) are positive.
Continuing with the backstepping procedure, we define the next error term as
e 3 = k 2 e 2 + η + p ( ϖ ^ ) ( F a s F d F ^ δ ) x ¨ d g + ξ ˙ + k 1 η ˙ + z a .
Defining the next Lyapunov candidate function by augmenting e 3 as
V 3 = V 2 b + 1 2 e 3 2 .
and computing the time derivative of V 3 , we have
V ˙ 3 = W 3 ( ξ , η , e 2 , e 3 ) e 2 z a + e 3 z ˙ a e 2 ϵ + e 3 [ k 3 e 3 + e 2 + k 2 e ˙ ^ 2 + η ˙ + p ( ϖ ^ ) [ A ¯ a s ( u τ γ A a s x 2 x 3 v a s 1 + Q ˙ ( γ 1 ) ( s p v a s ) 1 + d n ) α f x ˙ ^ 2 ( 1 t a n h 2 ( α x 2 + β s i g n ( x 1 ) ) ) c x ˙ ^ 2 F ^ ˙ δ ] + ξ ¨ + p ˙ ( ϖ ^ ) ( F a s F d F ^ δ ) x d + k 1 η ¨ ] ,
where z ˙ a = u ( t ) u ( t τ k ) , W 3 ( ξ , η , e 2 , e 3 ) = k 1 η 2 + α ξ 2 + k 2 e 2 2 + k 3 e 3 2 , where k 3 is a positive control gain, and
x ˙ ^ 2 = p ( ϖ ^ ) ( F a s F d F ^ δ ) g e ¨ ^ 1 = x ˙ ^ 2 x ¨ d e ˙ ^ 2 = e ¨ ^ 1 + k 1 η ˙ + e 1 + α e ˙ 1 .
We note that d n is an unknown perturbation and cannot be included in the control law, so it is replaced by the following inequality [22], using
| Υ | Υ t a n h ( Υ ϵ r ) 0.2875 ϵ r ,
where ϵ r is a positive number. Rewriting (33), we have
V ˙ 3 W 3 ( ξ , η , e 2 , e 3 ) e 2 z a + e 3 z ˙ a e 2 ϵ + 0.2758 ϵ r + e 3 [ k 3 e 3 + e 2 + k 2 e ˙ ^ 2 + η ˙ + p ( ϖ ^ ) [ A ¯ a s ( u γ A a s x 2 x 3 v a s 1 + d m a x t a n h ( e 3 A ¯ a s p ( ϖ ^ ) d m a x ϵ r ) + Q ˙ ( γ 1 ) ( s p v a s ) 1 ) F ^ ˙ δ α f x ˙ ^ 2 1 t a n h 2 ( α x 2 + β s i g n ( x 1 ) ) c x ˙ ^ 2 ] + ξ ¨ + p ˙ ( ϖ ^ ) ( F a s F d F ^ δ ) x d + k 1 η ¨ ] .
Now, we choose the delay-free control input as
u = A ¯ a s 1 [ p ( ϖ ^ ) 1 k 3 e 3 e 2 p ˙ ( ϖ ^ ) ( F a s F d F ^ δ ) + x d k 1 η ¨ ξ ¨ k 2 e ˙ ^ 2 η ˙ + α f x ˙ ^ 2 1 t a n h 2 ( α x 2 + β s i g n ( x 1 ) ) + c x ˙ ^ 2 + F ^ ˙ δ ] d m a x t a n h ( e 3 A ¯ a s d m a x p ( ϖ ^ ) ϵ r ) + γ A a s x 2 x 3 v a s 1 Q ˙ ( γ 1 ) ( s p v a s ) 1 .
Substituting (37) into (36), we obtain
V ˙ 3 k 1 η 2 α ξ 2 k 2 e 2 2 k 3 e 3 2 + 0.2785 ϵ r e 2 z a + e 3 z ˙ a .
where z ˙ a = u ( t ) u ( t τ k ) .
Remark 3. 
The introduced projection function (23) is of class C 1 , and ensures that p ( ϖ ^ ) satisfies p ( ϖ ^ ) [ ϖ m i n ϖ , ϖ m a x ϖ ] . Setting m m i n = 1 / ( ϖ m a x + ϖ ) and m m a x = 1 / ( ϖ m i n ϖ ) , we can conclude that the designed delayed-free input (37) is always well defined since p ( ϖ ^ ) is always nonzero.
In summary, the main result of the height tracking problem is summarized in the following theorem.
Theorem 1. 
Let x d R 3 in (4) be the proposed height whose time derivatives are bounded and continuous. By considering the closed-loop system obtained by the following control law (37), the RBFNN approximator online compensates for MRD force in (11), and the estimated sprung mass in (28). Then, the tracking errors e i = [ e 1 , e 2 , e 3 ] , given by (4), (19), and (31) can eventually converge to a small neighborhood of zero in the presence of modeling error d n and input delay τ, achieving uniform ultimate boundedness.
We start the proof by using mean-value theorem to (38), leading to
V ˙ 3 k 1 η 2 α ξ 2 θ 2 e 2 2 θ 3 e 3 2 1 4 ( z a 2 z ˙ a 2 ) + 0.2785 ϵ r ,
where θ 2 , θ 3 are positive constants. Furthermore, we define a differentiable positive-definite functional as
V = V 3 b + κ V L K ,
where λ is a positive constant, V L K is a Lyapunov–Krasovskii function, given as
V L K = t τ k t s t | u ( θ ) | 2 d θ ) d s ,
whose derivation derivative yields
V ˙ = k 1 η 2 α ξ 2 θ 2 e 2 2 θ 3 e 3 2 1 4 ( z a 2 z ˙ a 2 ) + 0.2785 ϵ r + κ τ k | u | 2 κ t τ k t | u ( θ ) | 2 d θ .
If the sufficient conditions in (41) are satisfied, then the following inequality can be defined
| u | ϕ + k 3 | e 3 | ,
where ϕ is a known positive constant, and using the fundamental inequality it gives us
κ τ k | u | 2 = κ τ k ϕ 2 + 2 k 3 κ τ k e 3 ϕ + κ τ k k 3 2 e 3 2 ( k 3 + 1 ) κ τ k ϕ 2 + κ τ k k 3 2 e 3 2 .
Utilizing the Cauchy–Schwarz inequality, the integral in (41) can be upper bounded as
κ t τ k t | u ( θ ) | 2 d θ κ 2 τ k | z a | 2 κ 2 t τ k t | u ( θ ) | 2 d θ κ 2 t τ k t | u ( θ ) | 2 d θ κ 2 τ k t τ k t s t | u ( θ ) | 2 d θ d s .
Substituting (44), and (45) into (41), V ˙ is upper bounded as
V ˙ = k 1 η 2 α ξ 2 Φ 2 e 2 2 Φ 3 e 3 2 Φ 4 z a 2 Φ L K V L K + Φ e ,
where Φ 2 = k 2 + 1 4 , Φ 3 = k 3 1 4 κ τ k k 3 2 , Φ 4 = 1 4 + κ 2 τ k , Φ L K = κ 2 τ k , Φ e = 0.2785 ϵ r + ( k 3 + 1 ) κ τ k ϕ 2 , and then setting z = [ η , ξ , e 2 , e 3 ] T , V ˙ 3 can be further upper bounded by
V ˙ k m i n z 2 + Φ e ,
where the auxiliary constant k m i n R 3 is defined as
k m i n m i n { k 1 , α , Φ 2 , Φ 3 , Φ 4 } .
Consequently, the above equation is ultimately negative for
z > Z = Φ e k m i n .
Remark 4. 
According to Equation (48), the ultimate boundary size can be made very small by choosing a larger value of gain k 3 . For arbitrary intervals of delay or arbitrarily large uncertain masses, the control gain required to satisfy the sufficient conditions in (28) and (44) depends in its size on the external level of road excitation and the inherent error ϕ resulting from u. Thus, to obtain a smaller error, a smaller α and a larger arbitrary control gain k 1 , k 2 , k 3 need to be used. However, it is necessary to find the appropriate balance between tracking accuracy and oscillator strength to avoid the controller causing unwanted oscillations, divergences, etc.
Remark 5. 
This work stabilises the conservativeness and computational complexity of the system by constructing a Lyapunov function that satisfies uncertain quality information and input delay information (as in (26) and (41)), and uses some suitable boundary techniques that can conservatively reduce the desired result to a minor value.

5. Simulation Validation

To verify the effectiveness of proposed control strategy, in this section, by using the toolbox provided by AMESim(2021.2), a virtual plant of quarter vehicle with the AAS system is established. To obtain more realistic simulation results, we use a built-in secondary development component, Submodel Editor, to create and update the mathematical logic embedded in the MRDs and the mass blocks. In addition, the mathematical model and the proposed controller for the MRD-AAS system are programmed in Matlab/Simulink(2019(b)). Figure 5 shows the control block diagram for co-simulation. Major parameters used in the co-simulation are given in Table 2. Following [5], the control law u τ is converted into the required changes of air mass m ˙ d e s , given by
m ˙ d e s = s p v a s u τ γ R T s , if u τ > 0 s p v a s u τ γ R T s , otherwise .
Remark 6. 
The secondary development programme in AMESim, mentioned in [23,24] of this work to satisfy specific simulation demands, creates the required application libraries and software (as in (3) and (21)) flexibly through secondary development using some suitable mathematical logic and modifying the output and input port signals for specific demands and computational complexity.
The corresponding PWM duty cycle ( D c ) of the chosen virtual control input u can be viewed as the control signal which can be determined as
Pro . I : D c = m i n ( s p v a s u τ m ˙ i n γ R T s , 1 ) Pro . II : D c = m i n ( s p v a s u τ m ˙ i n γ R T s , 1 ) Pro . III : D c = 0 .
In order to demonstrate the height tracking performance of the proposed control system in the presence of sprung mass variations and time-varying input delay simultaneously, the uncertain sprung mass is divided into five stages: (i) 0 to t 1 : m = 300 ( kg ) ; (ii) t 1 to t 2 : m = 330 ( kg ) ; (iii) t 2 to t 3 : m = 360 ( kg ) ; (iv) t 3 to t 4 : m = 330 ( kg ) ; and (v) t 4 to end: m = 300 ( kg ) . Meanwhile, with a current of I = 0 ( A ) , the input time delay is divided into three groups: τ 1 = τ k 1 + 5 s i n ( π t ) ( ms ) , τ 2 = τ k 2 + 5 c o s ( 0.5 π t ) ( ms ) , τ 3 = τ k 3 + 5 c o s ( π t ) s i n ( 0.5 π t ) ( ms ) , τ 4 = τ k 2 + 5 c o s ( 0.5 π t ) ( ms ) , τ 5 = τ k 1 + 5 s i n ( π t ) ( ms ) , with τ k 1 = 10 ( ms ) , τ k 2 = 20 ( ms ) , τ k 3 = 30 ( ms ) , and the disturbance occurring on the system is set as d n = 0.018 c o s ( 0.01 t ) . The co-simulation is carried out on a rough road corresponding to class B of an ISO road profile with a driving speed of 40 ( km / h ) .
The time evolution of height z s and tracking error e 1 with uncertain mass m s and time-varying delays of τ 1 , τ 2 and τ 3 under the driving speeds of 40 ( km / h ) and 100 ( km / h ) are displayed in Figure 6a,b, respectively. The proposed delay-free control law enables the height to converge towards the desired value within 2 ( s ) . The estimated error of sprung mass m ˜ is bounded by the range of 0.5 ( kg ) , as shown in Figure 7. Correspondingly, the tracking error e 1 undergoes new transients at t = t i , i = 1 , 2 , 3 , 4 , but always converges to a neighborhood of zero as time increases and is ultimately bounded by a tolerant value 1 ( mm ) . Moreover, the designed RBFNN approximator could compensate for the non-simulated dynamics of the MRD within the lagged time range 30 ( ms ) , as shown in Figure 8.
In order to further demonstrate the advantages of the proposed control strategy without consideration of uncertain mass and time-varying input delay, Figure 9 shows that the proposed controller and the method presented in [25] could drive the height tracking error close to zero under the same co-simulation parameters and road disturbances. The comparison results obtained from both two controllers are given in Table 3. It is noted that the root mean square (RMS) and response time (RT) of the ride height and the adjustment time are decreased by using the proposed controller. Simulation results in Figure 9 and the performance comparison in Table 3 indicate that the proposed control technique outperforms the HMPC given in [25].

6. Conclusions

This paper provided an innovative neural network-driven approach to tackle the issue of adjusting vehicle height with an MRD-AAS system in the presence of unmodeled dynamics of MRD, uncertain mass, and time-varying input delays. A RBFNN approximator was designed with the adaptive gradient descent learning algorithm to compensate for the concentrated modeling error of the MRD caused by the unmodeled dynamics. A delay-free control strategy was synthesized to enable the ride height converge to an arbitrarily small neighborhood of the preset desired height. The designed RBFNN approximator and mass law were integrated into the delay-free control input, achieving uniform ultimate boundedness. Co-simulation results were given to validate the performance and effectiveness of the proposed control strategy.
With respect to future work, a deep reinforcement learning-based optimal MRD-AAS control system will be designed to maintain leveling and posture motion control while improving performance in terms of isolating vibrations holding capacity and handling stability.

Author Contributions

Conceptualization, R.Z. and H.X.; data curation, X.G. and X.S.; formal analysis, R.Z. and H.X.; funding acquisition, R.Z. and X.S.; investigation, R.Z. and H.X.; methodology, H.X. and C.C.; project administration, R.Z. and X.S.; resources, H.X. and C.C.; Software, R.Z. and H.X.; supervision, X.S. and C.C.; validation, X.G. and C.C.; visualization, H.X. and X.G.; writing—original draft preparation, R.Z. and H.X.; writing—review and editing, X.S. and C.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded in part by the Science and Technology Foundation of Guizhou Province under grant numbers ZK [2021] 320 and ZK [2023] 260; in part by the National Natural Science Foundation of China under grant numbers 52072161 and 52175127; and in part by the International Science and Technology Cooperation Program of Jiangsu Province (Project No. BZ2022055).

Institutional Review Board Statement

Our research did not involve human participants or animals.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The datasets generated during and/or analyzed during the current study are not publicly available but are available from the corresponding author upon reasonable request.

Acknowledgments

All the authors are deeply grateful to the editors and anonymous referees for their valuable suggestions for improving the quality of this paper. The authors also gratefully acknowledge contributions from the School of Mechanical and Electrical Engineering, Guizhou Normal University, and Automotive Engineering Research Institute, Jiangsu University.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Li, W.; Xie, Z.; Wong, P.K.; Hu, Y.; Guo, G.; Zhao, J. Event-Triggered Asynchronous Fuzzy Filtering for Vehicle Sideslip Angle Estimation with Data Quantization and Dropouts. IEEE Trans. Fuzzy Syst. 2022, 30, 2822–2836. [Google Scholar] [CrossRef]
  2. Yang, J.; Wang, X.M.; Baldi, S.; Singh, S.; Fari, S. A Software-in-the-Loop Implementation of Adaptive Formation Control for Fixed-Wing UAVs. IEEE/CAA J. Autom. Sin. 2019, 6, 1230–1239. [Google Scholar] [CrossRef]
  3. Zhang, Z.Y.; Wang, J.B.; Wu, W.G.; Huang, C.X. Semi-active control of air suspension with auxiliary chamber subject to parameter uncertainties and time-delay. Int. J. Robust Nonlinear Control 2020, 30, 7130–7149. [Google Scholar] [CrossRef]
  4. Sun, L.Q.; Wang, Y.; Li, Z.X.; Geng, G.Q.; Liao, G. H∞ Robust Control of Interconnected Air Suspension Based on Mode Switching. IEEE Access 2022, 10, 62377–62390. [Google Scholar] [CrossRef]
  5. Zhao, R.C.; Xie, W.; Wong, P.K.; Cabecinhas, D.; Silvestre, C. Robust Ride Height Control for Active Air Suspension Systems with Multiple Unmodeled Dynamics and Parametric Uncertainties. IEEE Access 2019, 7, 59185–59199. [Google Scholar] [CrossRef]
  6. Spencer, B.F.; Dyke, S.J.; Sain, M.K.; Carlson, J.D. Phenomenological model for magnetorheological dampers. J. Eng. Mech. 1997, 123, 230–238. [Google Scholar] [CrossRef]
  7. Kwok, N.M.; Ha, Q.P.; Nguyen, T.H.; Li, J.; Samaln, B. A novel hysteretic model for magnetorheological fluid dampers and parameter identification using particle swarm optimization. Sens. Actuators A 2006, 132, 441–451. [Google Scholar] [CrossRef]
  8. Bai, X.X.; Chen, P.; Qian, L.J. Principle and validation of modified hysteretic models for magnetorheological dampers. Smart Mater. Struct. 2015, 24, 085014. [Google Scholar] [CrossRef]
  9. Hu, G.; Liu, Q.; Li, G.; Xu, M. Simulation and analysis of adjustable Sigmoid model for a typical magnetorheological damper. Mach. Tool Hydraul. 2016, 44, 9–15. [Google Scholar]
  10. Gao, Z.; Wong, P.K.; Zhao, J.; Hua, X.Q.; Ma, X.B.; Xie, Z.C. Design of compensatory backstepping controller for nonlinear magnetorheological dampers. Appl. Math. Modell. 2023, 114, 318–337. [Google Scholar] [CrossRef]
  11. Na, J.; Huang, Y.; Wu, X.; Su, S.F.; Li, G. Adaptive Finite-Time Fuzzy Control of Nonlinear Active Suspension Systems with Input Delay. IEEE Trans. Cybern. 2020, 50, 2639–2650. [Google Scholar] [CrossRef] [PubMed]
  12. Min, X.; Li, Y.; Tong, S. Adaptive fuzzy output feedback inverse optimal control for vehicle active suspension systems. Neurocomputing 2020, 403, 257–267. [Google Scholar] [CrossRef]
  13. Khan, L.; Qamar, S.; Khan, M.U. Comparative analysis of adaptive NeuroFuzzy control techniques for full car active suspension system. Arab. J. Sci. Eng. 2014, 39, 2045–2069. [Google Scholar] [CrossRef]
  14. Zhu, M.F.; Lv, G.; Zhang, C.P.; Jiang, J.M.; Wang, H.R. Delay-Dependent Sliding Mode Variable Structure Control of Vehicle Magneto-Rheological Semi-Active Suspension. IEEE Access 2022, 10, 51128–51141. [Google Scholar] [CrossRef]
  15. Yan, G.; Fang, M.X.; Xu, J. Analysis and experiment of time-delayed optimal control for vehicle suspension system. J. Sound Vib. 2019, 446, 144–158. [Google Scholar] [CrossRef]
  16. Wu, H.; Zheng, L.; Li, Y.N.; Zhang, Z.D.; Liang, Y.X.; Hu, Y.M. Comprehensive analysis for influence of complex coupling effect and controllable suspension time delay on hub-driving electric vehicle performance. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 2022, 236, 40–58. [Google Scholar] [CrossRef]
  17. Feng, Z.Q.; Yu, D.; Chen, Z.B.; Xing, X.D.; Yan, H. A New Hysteresis Model with Resistor–Capacitor Operator for Magnetorheological Dampers. J. Vib. Eng. Technol. 2023, 11, 3229–3244. [Google Scholar] [CrossRef]
  18. Hu, G.L.; Liu, Q.j.; Li, R.Q.; Li, G. Vibration control of semi-active suspension system with magnetorheological damper based on hyperbolic tangent model. Adv. Mech. Eng. 2017, 9, 1–15. [Google Scholar] [CrossRef]
  19. Corns, S.M.; James, K.; Derong, L.; David, F. Fundamentals of Computational Intelligence: Neural Networks, Fuzzy Systems, and Evolutionary Computation. Genet. Program. Evolvable Mach. 2017, 18, 119–120. [Google Scholar] [CrossRef]
  20. Liu, J.K. Radial Basis Function (RBF) Neural Network Control for Mechanical Systems; Tsinghua University Press: Beijing, China, 2013; pp. 34–39. [Google Scholar]
  21. Ge, S.S.; Wang, J. Robust adaptive neural control for a class of perturbed strict feedback nonlinear systems. IEEE Trans. Neural Netw. 2002, 13, 1409–1419. [Google Scholar] [CrossRef]
  22. Polycarpou, M.M. Stable adaptive neural control scheme for nonlinear systems. IEEE Trans. Autom. Control 1996, 41, 447–451. [Google Scholar] [CrossRef]
  23. LMS International Corporation. AMESim Reference Manual Rev 13: LMS Image; Siemens AG Lab: Munich, Germany, 2013. [Google Scholar]
  24. Fu, Y.; Qi, X. The Manual of System Modeling & Simulation Reference for LMS Imagine. LabAMESim; Beijing University of Aeronautics and Astronautics Press: Beijing, China, 2011; pp. 272–283. [Google Scholar]
  25. Sun, X.Q.; Cai, Y.F.; Chen, L.; Liu, Y.L.; Wang, S.H. Vehicle height and posture control of the electronic air suspension system using the hybrid system approach. Veh. Syst. Dyn. 2016, 54, 328–352. [Google Scholar] [CrossRef]
Figure 1. Schematic of quarter vehicle with MRD-AAS.
Figure 1. Schematic of quarter vehicle with MRD-AAS.
Sensors 24 00156 g001
Figure 2. MRD external dynamic tests: (a) Test bench; (b) Measured dynamic characteristics.
Figure 2. MRD external dynamic tests: (a) Test bench; (b) Measured dynamic characteristics.
Sensors 24 00156 g002
Figure 3. RBFNN approximator performance: (a) MRD force at 0 A; (b) RBFNN approximated error at 0 A; (c) MRD force at 1 A; (d) RBFNN approximated error at 1 A; (e) MRD force at 2 A; (f) RBFNN approximated error at 2 A; (g) MRD force at 3 A; (h) RBFNN approximated error at 3 A.
Figure 3. RBFNN approximator performance: (a) MRD force at 0 A; (b) RBFNN approximated error at 0 A; (c) MRD force at 1 A; (d) RBFNN approximated error at 1 A; (e) MRD force at 2 A; (f) RBFNN approximated error at 2 A; (g) MRD force at 3 A; (h) RBFNN approximated error at 3 A.
Sensors 24 00156 g003
Figure 4. Block diagram of proposed control strategy.
Figure 4. Block diagram of proposed control strategy.
Sensors 24 00156 g004
Figure 5. Control block diagram of co-simulation.
Figure 5. Control block diagram of co-simulation.
Sensors 24 00156 g005
Figure 6. Height tracking performance of quarter vehicle MRD-AAS in co-simulation: (a) At speed of 40 km/h; (b) At speed of 100 km/h.
Figure 6. Height tracking performance of quarter vehicle MRD-AAS in co-simulation: (a) At speed of 40 km/h; (b) At speed of 100 km/h.
Sensors 24 00156 g006
Figure 7. Time evolution of the estimated mass m ^ and mass estimation error m ˜ with input delays τ 1 , τ 2 , τ 3 .
Figure 7. Time evolution of the estimated mass m ^ and mass estimation error m ˜ with input delays τ 1 , τ 2 , τ 3 .
Sensors 24 00156 g007
Figure 8. MRD force comparison with time-varying input delay τ 3 in co-simulation.
Figure 8. MRD force comparison with time-varying input delay τ 3 in co-simulation.
Sensors 24 00156 g008
Figure 9. Performance comparison and height error on a quarter vehicle in co-simulation. HMPC in [25].
Figure 9. Performance comparison and height error on a quarter vehicle in co-simulation. HMPC in [25].
Sensors 24 00156 g009
Table 1. Nomenclature.
Table 1. Nomenclature.
SymbolDescriptionSymbolDescription
z s ( m ) sprung mass displacement D c duty cycle of PWM signal
z u ( m ) unsprung mass displacement T p ( t ) PWM signal period
z r ( m ) road disturbance m ˙ i n ( kg · s 1 ) inflating air mass rate
m s ( kg ) sprung mass m ˙ o u t ( kg · s 1 ) deflating air mass rate
m u ( kg ) unsprung mass m ˙ d e s ( kg · s 1 ) desired air mass change
A a s ( m 2 )surface area of air spring f ( M R D ) hysteresis loop scaling factor
A ¯ a s ( m 2 ) effective area of air spring α ( M R D ) hysteresis loop scaling factor
v a s ( m 3 ) air spring volume β ( M R D ) hysteresis loop half-width factor
F d ( N ) MRD force of hyperbolic model c ( M R D ) MRD damping coefficient
γ adiabatic index of ideal gases τ k ( ms ) delay-time
Rideal-gases constant x d ( m ) desired height
Q ˙ ( J · s 1 ) heat transfer rateuactuator command
P u ( Pa ) upstream pressure k 1 , k 2 , k 3 positive constants
P d n ( Pa ) downstream pressure m ˜ s ( kg ) mass estimation error
p a t m ( MPa ) atmosphere pressure m ^ s ( kg ) mass estimation
P c r ( Pa ) critical pressure ratio γ w constant coefficient
c q deflating air mass rate λ constant coefficient
T ( K ) air temperature F δ ( N ) modeling error of MRD
κ constant coefficient F ˜ δ ( N ) estimation error of MRD force
ϱ constant coefficient F ^ δ ( N ) estimation of MRD force
Table 2. Parameter used in simulation.
Table 2. Parameter used in simulation.
ParameterValueParameterValue
A a s 0.0072 ( m 2 ) α 1 1.8506
A h e a t 500 ( J / K s m 2 ) α 2 100.48
g 9.807 ( m / s 2 ) β 1 0.0819
γ 1.4 β 2 0.5245
P a t m 1.0133 ( bar ) c 1 651.23
s p 10 5 c 2 438.3
z 0 0.2047 ( m ) k 1 1
p 0 5.4 ( bar ) k 2 5
S x 10 5 ( m 2 ) k 3 10
C q 8 × 10 2 ϵ 1 1 × 10 2
R 287.1 z 0 i 0.05 ( m )
v 40 ( k m / h ) I0 A
G 0 0.64 × 10 4 d m a x 0.02
m u 30 ( kg ) t o l e r a n c e 1 × 10 5
f 1 160.88 α 1
f 2 32.8981 190
Table 3. Comparison of controller tracking performance.
Table 3. Comparison of controller tracking performance.
RMSE of z 1 RT of z 1
Times Range (s)Method in [25]Proposed MethodImprovement 1Method in [25]Proposed MethodImprovement 1
0 to 5 s 0.6273 0.1521 75.75 % 0.2549 0.1459 42.76 %
5 to 10 s 7.3058 2.7969 61.72 % 7.6781 5.7588 25.00 %
10 to 13 s 0.2974 0.1858 37.53 % ///
13 to 15 s 17.0576 6.8637 59.76 % 16.134814.317311.26%
15 to 20 s 1.0783 0.1626 84.92 % ///
Improvement 1 = (Proposed method—Method in [25])/Method in [25].
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, 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. Sensors 2024, 24, 156. https://0-doi-org.brum.beds.ac.uk/10.3390/s24010156

AMA Style

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. Sensors. 2024; 24(1):156. https://0-doi-org.brum.beds.ac.uk/10.3390/s24010156

Chicago/Turabian Style

Zhao, Rongchen, Haifeng Xie, Xinle Gong, Xiaoqiang Sun, and Chen Cao. 2024. "Neural Network-Based Adaptive Height Tracking Control of Active Air Suspension System with Magnetorheological Fluid Damper Subject to Uncertain Mass and Input Delay" Sensors 24, no. 1: 156. https://0-doi-org.brum.beds.ac.uk/10.3390/s24010156

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