Next Article in Journal
ECG Signal Features Classification for the Mental Fatigue Recognition
Previous Article in Journal
When Is a Graded Free Complex Exact?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Integrated Control of Spray System and Active Suspension Systems Based on Model-Assisted Active Disturbance Rejection Control Algorithm

1
School of Mechanical Engineering, Yanshan University, Qinhuangdao 066004, China
2
School of Electrical Engineering, Yanshan University, Qinhuangdao 066004, China
*
Author to whom correspondence should be addressed.
Submission received: 3 August 2022 / Revised: 10 September 2022 / Accepted: 14 September 2022 / Published: 19 September 2022

Abstract

:
Due to the lack of body stability of emergency rescue vehicles, their attitude stability is insufficient and they are unable to realize working while driving, resulting in low rescue efficiency. Aiming at the water tower fire truck, which is equipped with an active suspension system, the vehicle attitude stability is studied. First, combined with the active suspension system and spray system, a 13-DOF integrated dynamic model for the water tower fire truck is established. Using the model-assisted active disturbance rejection control method, the controllers are designed for the vertical displacement, pitch angle, and roll angle of the vehicle attitude. Then, the computer simulation is carried out to verify the effectiveness of this control method. Finally, the water spray obstacle crossing experiment is carried out with a JP32G water tower fire truck. The results show that when the vehicle runs over the triangular obstacle on one side and two sides in the integrated spray-active suspension mode, the peak–peak values of body pitch angle and roll angle are reduced by 10.9% and 23.2%, and 23.7% and 16.3%, respectively, compared with the passive hydro pneumatic suspension.

1. Introduction

With the acceleration of urbanization and the improvement of people’s living standards, disaster prevention and mitigation, public safety, and emergency hedging have attracted increasingly more attention. Emergency rescue is characterized by uncertainty, unexpectedness, and complexity. The chassis of existing emergency rescue vehicles is restricted by their poor mobility and insufficient vehicle attitude stability, which leads to low rescue efficiency. For example, when a water tower fire truck is working, the traditional suspension system cannot satisfy the requirement of vehicle posture stability. When working, the whole vehicle needs to be supported at a fixed place by outriggers, which makes it impossible for the vehicle to work while driving. The scope of work is limited and the rescue efficiency is not high. If the water tower fire truck is operated while driving, its operation efficiency can be improved, and the rescue time can be saved, which is of great significance to ensure the safety of personnel and property.
The water tower fire truck is equipped with a high-power fire pump and a fire monitor. The fire monitor will produce a jet reaction force, resulting in a large roll and pitch moment acting on the body. At the same time, the boom of the water tower fire truck will be deployed, resulting in the increase in the center of gravity of the whole vehicle. These factors seriously affect the pose stability of the vehicle. In addition, when the vehicle is driving, the road surface is uneven and sometimes full of explosives such as gravel. When working on uneven terrain, the vehicle may roll over [1]. Therefore, to realize the operation while driving, the water tower fire truck requires the body to have better posture stability, which puts forward higher requirements for the regulation performance of the vehicle suspension system. The suspension system is an important part of the vehicle chassis, and its performance directly determines the ride comfort, operating stability, and driving safety of the vehicle [2,3,4].
The active suspension system can adjust the output of the actuator in real-time according to the changes in road input and vehicle state [5,6,7]. Obtaining a better shock absorption and body posture control effect provides the possibility to realize the operation while driving the water tower fire truck.
At present, much research effort has been devoted to the control synthesis of active suspension systems, such as model predictive control [8], preview control [9], adaptive control [10,11], fault-tolerant control [12], and event-triggered control [13]. The above control methods have improved the performance of vehicles to varying degrees, and many scholars have studied vehicle attitude stability. Kim and Lee [14] proposed a new nonlinear controller to adjust the height of the vehicle sprung mass and to regulate the roll and pitch angles of the vehicle body by an air suspension system. Youn et al. [15] designed an attitude tracking controller with zero acceleration as the target and simulated the feasibility of the attitude tracking controller by using a 4-DOF half-car model. Westhuizen et al. [16] proposed the possibility of using slow active suspension control to reduce the body roll and, thus, reduce the rollover tendency. Tchamna [17] studied a variable-stiffness semi-active suspension system to control vehicle attitude under cornering conditions. Chen [18] proposed an attitude control strategy for heavy emergency rescue vehicles based on the road level. For different road grades, the control parameters of the active suspension system are changed to reduce the root-mean-square values of the vertical acceleration, pitch angle, and roll angle of the vehicle body. Youn [19] studied the attitude control method of the vehicle body and calculated the ideal driving posture of the vehicle to offset the lateral and longitudinal forces acting on the passengers. The simulation results showed that this method has great potential. He [1] analyzed the impact of fire monitors under different working conditions on the static rollover angle of fire engines. The results showed that when the fire monitor works under high pressure, the influence on the static rollover angle of the vehicle can reach 15.8°. By establishing a virtual prototype model, Sun et al. [20] analyzed the influence of jet reaction force on the lateral stability of urban main fire trucks. The results showed that when the water flow velocity of the fire monitor is 30 L/s to 120 L/s, the rollover angle of the fire engine ranges from 30.9° to 20.57°. Gong et al. [21] proposed a semi-active suspension variable damping control strategy for heavy vehicles using a 9-DOF vehicle mode with improved vehicle ride comfort and handling stability. Wang et al. [22] presented a vehicle attitude compensation algorithm based on the state observer for a vehicle semi-active suspension system, improved the ride comfort of magneto-rheological suspension, and optimized the vehicle body attitude under braking and steering states. Li et al. [23] developed a centralized-distributed control strategy with attitude information obtained by multi-sensor fusion. The proposed attitude control system improved the obstacle performance, mobility, and flexibility of the vehicles. Guo [24] made the multi-axle emergency rescue vehicle equivalent to a 3-degrees-of-freedom parallel mechanism; the action of the servo actuator was driven by solving the inverse solution of the parallel mechanism position. Zeng et al. [25] presented an aperiodic adaptive intermittent control scheme to stabilize the attitude of the vehicle by using the backstepping technique.
The above research on vehicle attitude mainly focused on body attitude control under steering or stationary conditions. Few papers have been published on the body stability of the water tower fire truck while driving. Most of the models used in the current literature have been quarter-car or half-car suspension models, and only the vertical attitude or the pitch angle has been considered, which is not practical in industrial applications. Although the whole vehicle suspension system model has been used in the literature [14,15,18,22,23,24], the dynamic model of the working device mounted on the vehicle has not been considered. These models are not accurate enough for the water tower fire truck. Although the literature [1,20] has studied the effects of the reaction force on the vehicle, they have not realized the control of the vehicle attitude. Therefore, in order to improve the stability of the high-pressure jet fire truck, it is necessary to study a new vehicle dynamics model and algorithm.
In this paper, the dynamic model and control method of the vehicle active suspension system is studied. The main contribution in this work is to explore a 13-DOF dynamic model of the spray-active suspensions integrated system considering factors such as the fire monitor reaction force and the boom posture, rather than the conventional active suspension. In addition, unlike the linear state observer without model reference information in the literature [26,27], we consider the factors such as fire monitor reaction force and boom posture, which is embedded into a nonlinear ESO. This can reduce the estimation and bandwidth burden of the observer, and the nonlinear fal function can accelerate the convergence rate of the observer. Thus, the anti-interference ability of the system is improved. The spray-active suspensions integrated system of the water tower fire truck in Figure 1 is different from conventional suspension; thus, the theoretical results-derived MADRC controller is particular in this study. Therefore, this paper intends to make the water tower fire truck have a high posture stability when moving, to realize fire truck operation while moving on an uneven road.
The rest of the paper is organized in the following structure. Section 2 presents a Dynamic model and a state–space representation of a spray-active suspensions integrated system. Section 3 shows the equations of the MADRC design method of the Vehicle Attitude Stability Controller. Section 4 describes the results obtained after simulating the attitude of the vehicle in types of conventional models and integrated models and in different control methods. Section 5 shows the Experiment and Results Analysis, and the discussion and conclusions are shown in Section 6 and Section 7, respectively.

2. Dynamic Model of Spray-Active Suspensions Integrated System

2.1. Reaction Force Calculation of Fire Monitor

The high pressure and high flow working environment of the fire monitor will certainly result in a large reaction force, which will affect the operation stability of the fire truck. Therefore, it is necessary to calculate the reaction force. The direction of the reaction force acting on the fire monitor is opposite to that of the water flow. According to the Navier–Stokes equation of the fluid, the ANSI Standard Specification gives the following general expressions for calculating the jet reaction force.
F = C . V . d d t ( ρ u ) d V + C . S . ρ u ( u d A ) + A i n p d A i n + A o u t p d A o u t A p i p e p a d A p i p e
where F represents the reaction force, ρ represents density of the fluids, u represents the velocity of the fluid, p represents the pressure at the center of unit area, pa represents the pressure of the environment around the pipe, V represents the volume of the control volume, A represents the area of the control volume, Ain represents the area that flows into the control volume, Aout represents the area that flows out of the control volume, Apipe represents the external area of the pipe, C.V. represents the control volume, and C.S. represents the control surface.
In the formula, there is no assumption or simplified calculation, and it is based on the most basic principles of the fluid force, so it is accurate and universal. There are many integral operations in the general formula, which are difficult to apply in the actual project. It needs to be further simplified.
Assume that the flow of water in the fire monitor is one-dimensional.
When the fluid reaches a steady state, the above equation can be simplified as [28]
F l = ρ e u e 2 A e A e ( p e p a )
where ρe represents the density of the fluids, Fl represents the steady-state jet reaction force, ue represents the fluid velocity, pe represents fluid pressure, Ae represents the area of the fire monitor outlet, and pa represents the surrounding pressure around the fire monitor.

2.2. Integrated Dynamic Model

The 13-DOF integrated dynamic model of the active suspension and the spray system is shown in Figure 1.
The 13 degrees of freedom in this model are 6 unsprung mass vertical movements; the vertical, roll, and pitch movements of the vehicle; the luffing, rotation of the boom, and the swing motion of the fire monitor in two directions. There are six suspension actuating units in the model; zbi represents vertical displacement at the hinge point of the suspension and chassis (i = 1,2,3,4,5,6), mwi represents each unsprung mass, zwi represents the vertical displacement of each unsprung mass, csi represents the damping coefficient, ksi represents the stiffness of the suspension, l represents the wheel tread, and kti represents the stiffness of the tire.
Assume that the mass of the boom is concentrated at the center of mass; according to Newton’s second law, the following equations of motion for the heave, pitch, and rolling of a vehicle can be obtained.
m b Z ¨ = k s 1 ( z w 1 z b 1 ) + c s 1 ( z ˙ w 1 z ˙ b 1 ) + U 1 + k s 2 ( z w 2 z b 2 ) + c s 2 ( z ˙ w 2 z ˙ b 2 ) + U 2 + k s 3 ( z w 3 z b 3 ) + c s 2 ( z ˙ w 3 z ˙ b 3 ) + U 3 + k s 4 ( z w 4 z b 4 ) + c s 2 ( z ˙ w 4 z ˙ b 4 ) + U 4 + k s 5 ( z w 5 z b 5 ) + c s 2 ( z ˙ w 5 z ˙ b 5 ) + U 5 + k s 6 ( z w 6 z b 6 ) + c s 2 ( z ˙ w 6 z ˙ b 6 ) + U 6 F l sin α
θ ¨ I θ = a [ k s 1 ( z w 1 z b 1 ) + c s 1 ( z ˙ w 1 z ˙ b 1 ) + U 1 + k s 2 ( z w 2 z b 2 ) + c s 2 ( z ˙ w 2 z ˙ b 2 ) + U 2 ] + b [ k s 3 ( z w 3 z b 3 ) + c s 2 ( z ˙ w 3 z ˙ b 3 ) + U 3 + k s 4 ( z w 4 z b 4 ) + c s 2 ( z ˙ w 4 z ˙ b 4 ) + U 4 ] + c [ k s 5 ( z w 5 z b 5 ) + c s 2 ( z ˙ w 5 z ˙ b 5 ) + U 5 + k s 6 ( z w 6 z b 6 ) + c s 2 ( z ˙ w 6 z ˙ b 6 ) + U 6 ] + F l [ cos α sin ( β + γ ) sin δ B + sin α ( d B cos δ cos γ ) ]
φ ¨ I φ = l 2 [ k s 1 ( z w 1 z b 1 ) + c s 1 ( z ˙ w 1 z ˙ b 1 ) + U 1 + k s 3 ( z w 3 z b 3 ) + c s 2 ( z ˙ w 3 z ˙ b 3 ) + U 3 + k s 5 ( z w 5 z b 5 ) + c s 2 ( z ˙ w 5 z ˙ b 5 ) + U 5 ] l 2   [ k s 2 ( z w 2 z b 2 ) + c s 2 ( z ˙ w 2 z ˙ b 2 ) + U 2 + k s 4 ( z w 4 z b 4 ) + c s 2 ( z ˙ w 4 z ˙ b 4 ) + U 4 + k s 6 ( z w 6 z b 6 ) + c s 2 ( z ˙ w 6 z ˙ b 6 ) + U 6 ] F l [ cos α cos ( β + γ ) B sin δ s B in α sin δ sin γ ]
where Z, θ, and φ represent vertical displacement, pitch angle, and roll angle of the vehicle, respectively; mb represents a sprung mass; mb = ma1 + ma2; ma1 represents the chassis part of the sprung mass; ma2 represents the boom part of the sprung mass; α represents the angle between the jet reaction force and horizontal direction; β represents the angle between the jet reaction force and vertical direction; a represents the distance between the front axle and mass center of the chassis; b represents the distance between the middle axle and mass center of the chassis; c represents the distance between the rear axle and the mass center of the chassis; d represents the distance between the hinge of the boom and chassis and the mass center of the chassis; δ represents the luffing angle of the boom; γ represents the rotating angle of the boom; B represents the length of the boom; B1 represents distance from the hinge point of the boom and chassis to the mass center of the boom; Iθ1 represents the pitching moment of inertia of the chassis part; Iφ1 represents the rolling moment of inertia of the chassis part; Iθ represents the equivalent pitching moment of inertia of the vehicle; Iφ represents the equivalent rolling moment of inertia of the vehicle body.
The calculation method of the equivalent pitch and roll moment of inertia of the vehicle body is as follows:
I θ = I θ 1 + m a ( B 1 sin δ ) 2 + ( d B 1 cos δ cos γ ) 2
I φ = I φ 1 + m a ( B 1 sin δ ) 2 + ( B 1 cos δ sin γ ) 2
The vertical dynamic equation of the six wheels is
m w i Z ¨ w i = k s i ( z b i z w i ) k t i ( z w i z r i ) + c s i ( z ˙ b i z ˙ w i ) U i
where zri represents the road input of the ith wheel.
As the pitch and roll angles of the vehicle are small, it can be considered that sinθθ, sinφφ. Suppose that the vehicle is a rigid body; according to the motion law of a rigid body, the vertical displacement of the hinge joint between the suspension and chassis can be obtained:
Z b 1 = Z b a sin θ l / 2 sin φ Z b a θ l / 2 φ
Z b 2 = Z b a sin θ + l / 2 sin φ Z b a θ + l / 2 φ
Z b 3 = Z b + a sin θ l / 2 sin φ Z b + a θ l / 2 φ
Z b 4 = Z b + a sin θ + l / 2 sin φ Z b + a θ + l / 2 φ
Z b 5 = Z b + a sin θ l / 2 sin φ Z b + a θ l / 2 φ
Z b 6 = Z b + a sin θ + l / 2 sin φ Z b + a θ + l / 2 φ
Select the vehicle system state vector X as follows:
X = [ Z   θ   φ   z b 1   z b 2   z b 3   z b 4   z b 5   z b 6   z w 1   z w 2   z w 3   z w 4   z w 5   z w 6   Z ˙   θ ˙   φ ˙   z ˙ w 1   z ˙ w 2   z ˙ w 3   z ˙ w 4   z ˙ w 5 ˙ z w 6 ]
Then, the system equation is expressed as
X ˙ = A X + B U + W
where A is the coefficient matrix of vector X, B is the coefficient matrix of vector U, W represents external disturbance, and U is the road profile of each wheel and the active output force of each actuator.
U = [ z r 1   z r 2   z r 3   z r 4   z r 5   z r 6   U 1   U 2   U 3   U 4   U 5   U 6 ]
W = [ w 1   w 2   w 3 ]
w 1 = F l sin α
w 2 = F l [ cos α sin β sin δ B + sin α ]
w 3 = F l cos α cos β sin δ B  
Select the vertical displacement, pitch angle, roll angle, and the corresponding speed, acceleration of the vehicle body, the suspension, and the tire displacement to form the output vector Y:
Y = [ Z   θ   φ   Z ˙   θ ˙   φ ˙   Z ¨   θ ¨   φ ¨   z b 1   z b 2   z b 3   z b 4   z b 5   z b 6   z w 1   z w 2   z w 3   z w 4   z w 5   z w 6 ]
The output equation of the system can be expressed as
Y = C X + D U
where C is a coefficient matrix of vector X and D is a coefficient matrix of vector U.

3. Design of Vehicle Attitude Stability Controller

The ADRC method is used to design the vehicle attitude controller. ADRC consists of three main parts [29,30]: extended state observer (ESO), feedback controller, and disturbance rejection law.
The available system model information can be used to reduce the estimation time of the ESO and improve the ability of disturbance rejection. Accordingly, the jet action force is regarded as a model-assisted part, which is embedded in the ESO and integrated into the control law to achieve the desired output [31,32].
According to the derived vehicle dynamics equations, let x1 = Z, x2 = θ, and x3 = φ. Thus, Equations (3)–(5) can be rewritten as follows:
x ¨ 1 = Z ¨ = f 1 ( Z   θ   φ   Z ˙   θ ˙   φ ˙   z r 1 z r 2 z r 6 ) F l sin α m b + U 1 m b + U 2 m b + U 3 m b + U 4 m b + U 5 m b + U 6 m b
x ¨ 2 = θ ¨ = f 1 ( Z   θ   φ   Z ˙   θ ˙   φ ˙   z r 1 z r 2 z r 6 ) + F l [ cos α sin β sin δ B + sin α ( d B cos δ ) ] I θ a U 1 I θ a U 2 I θ + a U 3 I θ + a U 4 I θ + a U 5 I θ + a U 6 I θ
x ¨ 3 = φ ¨ = f 1 ( Z   θ   φ   Z ˙   θ ˙   φ ˙   z r 1 z r 2 z r 6 ) F l cos α cos β sin δ B I φ + l U 1 2 I φ l U 2 2 I φ + l U 3 2 I φ l U 4 2 I φ + l U 5 2 I φ l U 6 2 I φ
Let X = [x1 x2 x3]T, take u = [ U1 U2 U3 U4 U5 U6]T as the input of the suspension actuator, and take Y = [x1 x2 x3]T as system output; the state space expression is
X ¨ = F ( Z   θ   φ   Z ˙   θ ˙   φ ˙   z r 1 z r 2 z r 6 ) + F l H + B u Y = X
where
F ( Z   θ   φ   Z ˙   θ ˙   φ ˙   z r 1 z r 2 z r 6 ) = f 1 ( Z   θ   φ   Z ˙   θ ˙   φ ˙   z r 1 z r 2 z r 6 ) f 2 ( Z   θ   φ   Z ˙   θ ˙   φ ˙   z r 1 z r 2 z r 6 ) f 3 ( Z   θ   φ   Z ˙   θ ˙   φ ˙   z r 1 z r 2 z r 6 )
B = 1 m b 1 m b 1 m b 1 m b 1 m b 1 m b a I θ a I θ b I θ b I θ c I θ c I θ l 2 I φ l 2 I φ l 2 I φ l 2 I φ l 2 I φ l 2 I φ
H = sin α m b cos α sin β sin δ B + sin α ( d B cos δ ) I θ cos α cos β sin δ B I φ
The coupling between each suspension output force Bu in the above system is called “Static coupling”, and the coupling that exists in the Vehicle model is called the “Dynamic coupling”. Let S = Bu be the virtual output of the system; then, Equation (27) can be rewritten as
X ¨ = F ( Z   θ   φ   Z ˙   θ ˙   φ ˙   z r 1   z r 2 z r 6 ) + F l H + S Y = X
Therefore, the input–output relationship of the ith channel of the system is expressed as
x ¨ i = f i ( Z   θ   φ   Z ˙   θ ˙   φ ˙   z r 1 z r 2 z r 6 ) + F l H i + S i y i = x i
where i = (1,2,3), respectively, represents the vertical displacement, the pitch angle, and the roll angle of the vehicle. In each channel, the relationship between the input and output is SISO, and the static coupling in the system is decoupled. Dynamic coupling and other disturbances are regarded as system disturbances, and a second-order extended state observer (ESO) is designed to estimate the system disturbance as follows:
e i 1 = Z i 1 y i 1 f e i 1 = f a l ( e i 1 , ζ , h i ) Z ˙ i 1 = Z i 2 β i 1 e i 1 Z ˙ i 2 = Z i 3 β i 2 f e i 1 + F l H i + b i 1 S i Z ˙ i 3 = β i 3 f e i 1
where ei1 is the output error of each channel; Zi3 is the estimated disturbance of the system; hi is the length of the sampling step; Zi1 and Zi2 are system state observations; bi1 is the compensation factor; fal is the nonlinear feedback function; βi1, βi2, and βi3 are adjustable parameters of the extended state observer. The parameters are adjusted by the Fibonacci series and engineering experience [33].
f a l ( e , ρ , h ) = e τ s i g n ( e ) e h e h 1 τ e h
Design the state error control law as
e i 1 = v i 1 Z i 1 e i 2 = v i 2 Z i 2 F 0 i = β 1 f a l ( e i 1 , α 1 , ρ ) + β 2 f a l ( e i 2 , α 2 , ρ ) F i = F 0 i ( Z i 3 + F l H i ) / b i 1
Where β1 and β2 are gains of the controller.
The principle of a single-channel MADRC is shown in Figure 2.
As the control objective of each channel is to stabilize vertical displacement, the pitch angle, and the roll angle, i.e., v1i = 0, a tracking differential does not need to be designed in the controller to extract the input.
Three virtual outputs Si are calculated by the three ADRCs designed above, and the suspension system studied in this paper is a 6-input 3-output system. Therefore, the 3 outputs need to be decoupled. As the input coefficient matrix B of the system is irreversible, the output of the system is decoupled by solving the pseudo-inverse matrix of B. Let T = BT(BBT)−1; then, the input force of each actuator is Uj =TSi, (i = 1,2,3; j = 1,2,…,6). The structure of the MADRC system of the vehicle is shown in Figure 3.
Vertical, pitch, and roll controllers are designed using the MADRC integrated control strategy. The parameters of the spray system are added into the ESO, which reduces the estimation burden of the ESO, improves the system control performance, and realizes the adjustment and control of vehicle attitude.

4. Simulation Analysis

To verify the effectiveness of the three designed controllers, the control effect of the active suspension system and passive suspension system is simulated. The required parameters in the vehicle simulation model are shown in Table 1.
Random white noise as shown in Figure 4 is taken as the road excitation, the driving speed of the vehicle is set to 10 km/h, and the instantaneous flow rate of the fire monitor is set to a sinusoidal form with an amplitude of 0–120 L/min, as shown in Figure 5a. The simulation time is set to 10 s, and the simulation results of the output response of the vehicle attitude are shown in Figure 5b and Figure 6.
Figure 5b and Figure 6 show the comparison curves of vertical displacement, pitch angle, and roll angle of the mass center of the chassis. It can be seen from the figure that the vertical displacement, pitch angle, and roll angle of the integrated spray-active suspension system under MADRC control are significantly optimized compared with the conventional active suspension system under ADRC control and passive suspension.
The quantitative analysis of the curves is carried out, and the root-mean-square values of the vertical displacement, pitch angle, and roll angle of the vehicle are calculated. The results are shown in Table 2.
It can be seen from Table 2 that compared with the passive suspension system, the vertical displacement, pitch angle, and roll angle of the traditional active suspension system under ADRC control decrease by 56.78%, 37.84%, and 56.92%, respectively, and those of the integrated spray-active suspension system under MADRC control are further reduced by 13.42%, 41.38%, and 20.65%, respectively, compared with the passive suspension system. Each evaluation index is reduced by 70.20%, 79.22%, and 77.57% respectively.
In order to verify the advantages of the integrated model compared with the traditional suspension model in the literature [21], sinusoidal wave road excitation is implemented in the simulation. The road excitation of the two sides of the vehicle is shown in Figure 7. Other required parameters in the simulation model are the same as before. Figure 8 and Figure 9 show the comparison curves of vertical displacement, pitch angle, and roll angle of the mass center of the chassis, and the peak–peak values are presented in Table 3.
As shown in the previous figures, the proposed spray-active suspensions integrated model presents a better performance when compared with the conventional active suspension model. The vertical displacement, pitch angle, and roll angle of the vehicle are decreased by 19.8%, 20.7%, and 45.6%, respectively.

5. Experiments and Results Analysis

In order to further verify the effectiveness of the controller for improving the vehicle attitude, a barrier-crossing experiment was carried out on a water tower fire truck JP32G. The suspension system of the vehicle has active and passive working modes. The active mode is driven by a hydraulic servo and the passive mode is the hydro-pneumatic suspension system. The structure of the integrated control system is shown in Figure 10.
The controller uses a SCM9022 main board, which meets the standard PC/104 bus standard. The data acquisition system uses a ADT882-AT expansion board and provides 32 16-bit analog channels at a sampling speed of 200 kHz. The inertia measuring unit adopts an MTi 300AHRS IMU and the pressure sensor uses a PM95B series pressure transmitter with a measuring range of 0–40 MPa and accuracy of <0.25%. The actuator is a hydraulic cylinder; p1i and p2i (i = 1, 2, …, 6) are pressures of the rod and rodless cavity, respectively. Basic parameters of the tested vehicle are shown in Table 4.
The water spray obstacle crossing test was carried out on the water tower fire truck working in the integrated spray-active suspension system mode and the passive suspension mode. The test pavement was cement pavement and tire excitation was achieved by setting triangular obstacles on the pavement. Figure 11 shows the water spray obstacle crossing test vehicle and site. The triangular obstacle is shown in Figure 12.
During the experiment, triangular obstacles were placed on one or two sides in front of the vehicle. The control program was started, the vehicle was moved over the triangular obstacles, and changes in body pitch and roll angle were recorded. The control procedure was switched off, the suspension system was switched to hydro-pneumatic passive suspension, the vehicle was moved over the triangular obstacles at the same speed and manner, and the pitch and roll angle of the vehicle were recorded.
Figure 13 shows the pitch and roll angle curves of the vehicle with the wheels on one side crossing a triangular obstacle at a speed of 3.2 km/h. Figure 14 shows the pitch and roll angle curves of the vehicle with the wheels of two sides crossing a triangular obstacle at a speed of 3.2 km/h.
From the test comparison curves, it can be seen that the integrated spray-active suspension under MADRC control significantly improves the pitch and roll stability of the vehicle compared with the passive suspension system when the vehicle passes through uneven road surfaces.
Quantitative analysis of Figure 13 and Figure 14 is carried out, and the peak–peak value of vehicle attitude is shown in Table 5.
It can be seen from Table 5 that compared with the passive suspension system, the peak–peak values of pitch angle and roll angle of the integrated spray-active suspension system under MADRC control decrease by 10.9% and 23.2%, respectively, when the vehicle runs over triangular obstacles on one side, and decrease by 23.7% and 16.3%, respectively, when the vehicle runs over triangular obstacles on two sides.

6. Discussion

In the 13-DOF integrated dynamic model studied in this paper, it is assumed that the boom is a rigid body, but in fact, the boom will produce large deformation under the torque due to the clearance and compliance of the boom. This directly affects the spray angle of the fire monitor at the top of the boom and the action point of the spray reaction force, so the elastic deformation of the boom is a critical factor to be considered in future work.
A model-assisted nonlinear disturbance rejection controller is proposed and applied to water tower fire trucks. The innovation lies in introducing the known vehicle dynamic information into the extended state observer to reduce the bandwidth burden of the observer and improve the overall control performance. Due to the complicated theoretical analysis of the nonlinear ADRC (the nonlinearity of the fal function), the stability analysis of the disturbance rejection controller is still in the development stage, and there is relatively little related theoretical research. Analysis methods include the self-stable region approach [34,35], frequency domain method [36,37], and Popov hyperstability theory [38,39]. The difficulty of stability analysis theory leads to the limited application of nonlinear ADRC. However, it has been proved in practice that the nonlinear ADRC can stabilize the system and has higher convergence efficiency than the linear form. From the perspective of engineering application, this paper preliminarily explores its application in the field of active suspension. However, the stability analysis of the model-assisted expanding state observer is still under further study, and this is also work to be carried out in the future.

7. Conclusions

In this paper, a 13-DOF integrated spray-active suspension dynamic model of the water tower fire truck is established. The MADRC control method is used to design the vehicle attitude controller, and computer simulation and real vehicle tests are carried out.
The simulation results show that when the vehicle speed is 10 km/h, the vertical displacement, pitching angle, and roll angle of the integrated spray-active suspension system under the MADRC control decrease by 70.20%, 79.22%, and 77.57%, respectively, compared with the passive suspension system, and decrease by 19.8%, 20.7%, and 45.6%, respectively, compared with the conventional active suspension model under sinusoidal wave road excitation.
The experimental results show that when the vehicle runs over the triangular obstacle on one side and two sides in the integrated spray-active suspension mode, the peak–peak values of pitch angle and roll angle are reduced by 10.9% and 23.2%, and 23.7% and 16.3%, respectively, compared with the passive hydro pneumatic suspension.

Author Contributions

Conceptualization, D.Z. and J.Z.; methodology, D.Z. and J.Z.; software, S.L.; validation, S.L. and Z.Z.; investigation, D.Z. and J.Z.; resources, D.Z. and S.L; formal analysis, J.C.; data curation, Z.Z.; writing—original draft preparation, G.L. and J.C.; writing—review and editing, D.Z. and J.Z.; funding acquisition, D.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Regional Innovation and Development Funds of the National Natural Science Foundation of China (U20A20332), and Innovative research group of the Natural Science Foundation of Hebei Province, China (E202003174).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. He, M.; Jia, C. Analysis on roll stability of urban main battle fire truck under working conditions. Fire Sci. Technol. 2018, 37, 3. [Google Scholar]
  2. Zirkohi, M.M.; Lin, T. Interval type-2 fuzzy-neural network indirect adaptive sliding mode control for an active suspension system. Nonlinear Dyn. 2015, 79, 513–526. [Google Scholar] [CrossRef]
  3. Gohrle, C.; Schindler, A.; Wagner, A.; Sawodny, O. Road Profile Estimation and Preview Control for Low-Bandwidth Active Suspension Systems. IEEE/ASME Trans. Mechatron. 2015, 20, 2299–2310. [Google Scholar] [CrossRef]
  4. Pan, H.; Jing, X.; Sun, W. Robust finite-time tracking control for nonlinear suspension systems via disturbance compensation. Mech. Syst. Signal Process. 2017, 88, 49–61. [Google Scholar] [CrossRef]
  5. Du, M.; Zhao, D.; Yang, B.; Wang, L. Terminal sliding mode control for full vehicle active suspension systems. J. Mech. Sci. Technol. 2018, 32, 2851–2866. [Google Scholar] [CrossRef]
  6. Mustafa, G.I.Y.; Wang, H.P.; Tian, Y. Vibration control of an active vehicle suspension systems using optimized model-free fuzzy logic controller based on time delay estimation. Adv. Eng. Softw. 2019, 127, 141–149. [Google Scholar] [CrossRef]
  7. Zhang, Y.; Liu, Y.; Liu, L. Minimal learning parameters-based adaptive neural control for vehicle active suspensions with input saturation. Neurocomputing 2020, 396, 153–161. [Google Scholar] [CrossRef]
  8. Yao, J.; Wang, M.; Li, Z.; Jia, Y. Research on Model Predictive Control for Automobile Active Tilt Based on Active Suspension. Energies 2021, 14, 671. [Google Scholar] [CrossRef]
  9. Li, P.; Lam, J.; Cheung, K.C. Multi-objective control for active vehicle suspension with wheelbase preview. J. Sound Vib. 2014, 333, 5269–5282. [Google Scholar] [CrossRef]
  10. Liu, Y.; Zeng, Q.; Tong, S.; Chen, C.L.P.; Liu, L. Adaptive Neural Network Control for Active Suspension Systems with Time-Varying Vertical Displacement and Speed Constraints. IEEE Trans. Ind. Electron. 2019, 66, 9458–9466. [Google Scholar] [CrossRef]
  11. Liu, Y.; Zeng, Q.; Tong, S.; Chen, C.L.P.; Liu, L. Actuator Failure Compensation-Based Adaptive Control of Active Suspension Systems with Prescribed Performance. IEEE Trans. Ind. Electron. 2020, 67, 7044–7053. [Google Scholar] [CrossRef]
  12. Li, H.; Liu, H.; Gao, H.; Shi, P. Reliable Fuzzy Control for Active Suspension Systems with Actuator Delay and Fault. IEEE Trans. Fuzzy Syst. 2012, 20, 342–357. [Google Scholar] [CrossRef]
  13. Li, H.; Zhang, Z.; Yan, H.; Xie, X. Adaptive Event-Triggered Fuzzy Control for Uncertain Active Suspension Systems. IEEE Trans. Cybern. 2019, 49, 4388–4397. [Google Scholar] [CrossRef] [PubMed]
  14. Kim, H.; Lee, H. Height and Leveling Control of Automotive Air Suspension System Using Sliding Mode Approach. IEEE Trans. Veh. Technol. 2011, 60, 2027–2041. [Google Scholar]
  15. Youn, I.; Im, J.; Tomizuka, M. Level and attitude control of the active suspension system with integral and derivative action. Veh. Syst. Dyn. 2006, 44, 659–674. [Google Scholar] [CrossRef]
  16. Van der Westhuizen, S.F.; Els, P.S. Slow active suspension control for rollover prevention. J. Terramech. 2013, 50, 29–36. [Google Scholar] [CrossRef]
  17. Tchamna, R.; Lee, M.; Youn, I. Attitude control of full vehicle using variable stiffness suspension control. Optim. Control Appl. Methods 2015, 36, 936–952. [Google Scholar] [CrossRef]
  18. Chen, H.; Gong, M.; Zhao, D.; Zhu, J. Body attitude control strategy based on road level for heavy rescue vehicles. Proc. Inst. Mech. Eng. Part D 2021, 235, 1351–1363. [Google Scholar] [CrossRef]
  19. Youn, I.; Wu, L.; Youn, E.; Tomizuka, M. Attitude motion control of the active suspension system with tracking controller. Int. J. Automot. Technol. 2015, 16, 593–601. [Google Scholar] [CrossRef]
  20. Sun, J.; Li, W.; He, M. Analysis of Fire Water Monitor Jet Reaction Forces and Their Influences on the Roll Stabilities of Urban Firefighting Vehicle. Fire Technol. 2019, 55, 2547–2566. [Google Scholar] [CrossRef]
  21. Gong, M.; Chen, H. Variable damping control strategy of a semi-active suspension based on the actuator motion state. J. Low Freq. Noise Vib. Act. Control 2020, 39, 787–802. [Google Scholar] [CrossRef]
  22. Wang, R.; Sheng, F.; Ding, R.; Meng, X.; Sun, Z. Vehicle attitude compensation control of magneto-rheological semi-active suspension based on state observer. Proc. Inst. Mech. Eng. Part D 2021, 235, 3299–3313. [Google Scholar] [CrossRef]
  23. Li, B.; Zheng, G.; Wang, Z. Attitude Control of the Vehicle with Six In-Wheel Drive and Adaptive Hydro Pneumatic Suspensions. In Proceedings of the WCX SAE World Congress Experience, Detroit, MI, USA, 9–11 April 2019. [Google Scholar]
  24. Guo, Q.; Zhao, D.; Zhao, X.; Li, Z.; Shi, X. Active Suspension Control Strategy of Multi-Axle Emergency Rescue Vehicle Based on Inertial Measurement Unit. Sensors 2021, 21, 6877. [Google Scholar] [CrossRef] [PubMed]
  25. Zeng, Q.; Ma, D.; Zhao, J. Adaptive Aperiodic Intermittent Control for Nonlinear Active Suspension Systems. In Proceedings of the Chinese Control Conference, Shanghai, China, 26–28 July 2021; pp. 1437–1441. [Google Scholar]
  26. Khodamipour, G.; Khorashadizadeh, S.; Farshad, M. Observer-based adaptive control of robot manipulators using reinforcement learning and the Fourier series expansion. Trans. Inst. Meas. Control 2021, 43, 2307–2320. [Google Scholar] [CrossRef]
  27. Izadbakhsh, A.; Khorashadizadeh, S.; Kheirkhahan, P. Tracking Control of Electrically Driven Robots Using a Model-free Observer. Robotica 2018, 37, 729–755. [Google Scholar] [CrossRef]
  28. Liu, B. Calculation of Blowdown Thrust Force due to The Rupture of Main Pipe in Nuclear Power Plant. Ph.D. Thesis, Harbin Institute of Technology, Harbin, China, 2014. [Google Scholar]
  29. Hasbullah, F.; Faris, W.F. Simulation of disturbance rejection control of half-car active suspension system using active disturbance rejection control with decoupling transformation. J. Phys. Conf. Ser. 2017, 949, 12025. [Google Scholar] [CrossRef]
  30. Herbst, G. A Simulative Study on Active Disturbance Rejection Control (ADRC) as a Control Tool for Practitioners. Electronics 2013, 2, 246–279. [Google Scholar] [CrossRef]
  31. Li, S.; Cao, M.; Li, J.; Cao, J.; Lin, Z. Sensorless-Based Active Disturbance Rejection Control for a Wind Energy Conversion System with Permanent Magnet Synchronous Generator. IEEE Access 2019, 7, 122663–122674. [Google Scholar] [CrossRef]
  32. Sobhy, A.; Lei, D. Model-assisted active disturbance rejection controller for maximum efficiency schemes of DFIG-based wind turbines. Int. Trans. Electr. Energy Syst. 2021, 31, e13107. [Google Scholar] [CrossRef]
  33. Wang, L.; Zhao, D.; Liu, F.; Liu, Q.; Meng, F. Linear active disturbance rejection control for electro-hydraulic proportional position synchronous. Control Theory Appl. 2018, 35, 1618–1625. [Google Scholar]
  34. Huang, Y. A new synthesis method for uncertain systems the Self-Stable Region approach. Int. J. Syst. Sci. 1999, 30, 33–38. [Google Scholar] [CrossRef]
  35. Zhao, L.; Yang, Y.; Xia, Y.; Liu, Z. Active Disturbance Rejection Position Control for a Magnetic Rodless Pneumatic Cylinder. IRE Trans. Ind. Electron. 2015, 62, 5838–5846. [Google Scholar] [CrossRef]
  36. Li, J.; Qi, X.; Xia, Y.; Pu, F.; Chang, K. Frequency domain stability analysis of nonlinear active disturbance rejection control system. ISA Trans. 2015, 56, 188–195. [Google Scholar] [CrossRef] [PubMed]
  37. Wu, D.; Chen, K. Limit cycle analysis of active disturbance rejection control system with two nonlinearities. ISA Trans. 2014, 53, 947–954. [Google Scholar] [CrossRef] [PubMed]
  38. Li, J.; Xia, Y.; Qi, X.; Gao, Z. On the necessity, Scheme, and Basis of the Linear-Nonlinear Switching in Active Disturbance Rejection Control. IRE Trans. Ind. Electron. 2017, 64, 1425–1435. [Google Scholar] [CrossRef]
  39. Li, J.; Xia, Y.; Qi, X.; Zhao, P. Robust absolute stability analysis for interval nonlinear active disturbance rejection based control system. ISA Trans. 2017, 69, 122–130. [Google Scholar] [CrossRef] [PubMed]
Figure 1. 13-DOF Integrated dynamic model.
Figure 1. 13-DOF Integrated dynamic model.
Mathematics 10 03391 g001
Figure 2. Principle of MADRC.
Figure 2. Principle of MADRC.
Mathematics 10 03391 g002
Figure 3. Structure of MADRC system.
Figure 3. Structure of MADRC system.
Mathematics 10 03391 g003
Figure 4. (a) Road profile of left side under stochastic signals; (b) road profile of right side under stochastic signals.
Figure 4. (a) Road profile of left side under stochastic signals; (b) road profile of right side under stochastic signals.
Mathematics 10 03391 g004
Figure 5. (a) Spray flow curve; (b) comparison curve of vehicle vertical displacement under stochastic signals.
Figure 5. (a) Spray flow curve; (b) comparison curve of vehicle vertical displacement under stochastic signals.
Mathematics 10 03391 g005
Figure 6. (a) Comparison curve of vehicle pitch angle under stochastic signals; (b) comparison curve of vehicle roll angle under stochastic signals.
Figure 6. (a) Comparison curve of vehicle pitch angle under stochastic signals; (b) comparison curve of vehicle roll angle under stochastic signals.
Mathematics 10 03391 g006
Figure 7. (a) Road profile of left side under sinusoidal signals; (b) road profile of right side under sinusoidal signals.
Figure 7. (a) Road profile of left side under sinusoidal signals; (b) road profile of right side under sinusoidal signals.
Mathematics 10 03391 g007
Figure 8. (a) Comparison curve of vehicle vertical displacement under sinusoidal signals; (b) comparison curve of vehicle pitch angle under sinusoidal signals.
Figure 8. (a) Comparison curve of vehicle vertical displacement under sinusoidal signals; (b) comparison curve of vehicle pitch angle under sinusoidal signals.
Mathematics 10 03391 g008
Figure 9. Comparison curve of vehicle roll angle under sinusoidal signals.
Figure 9. Comparison curve of vehicle roll angle under sinusoidal signals.
Mathematics 10 03391 g009
Figure 10. Structure of the integrated control system.
Figure 10. Structure of the integrated control system.
Mathematics 10 03391 g010
Figure 11. Test vehicle and site.
Figure 11. Test vehicle and site.
Mathematics 10 03391 g011
Figure 12. Triangle obstacle.
Figure 12. Triangle obstacle.
Mathematics 10 03391 g012
Figure 13. (a) Pitch angle comparison curve of wheel-crossing obstacle on one side; (b) roll angle comparison curve of wheel-crossing obstacle on one side.
Figure 13. (a) Pitch angle comparison curve of wheel-crossing obstacle on one side; (b) roll angle comparison curve of wheel-crossing obstacle on one side.
Mathematics 10 03391 g013
Figure 14. (a) Pitch angle comparison curve of wheel-crossing obstacle on two sides; (b) roll angle comparison curve of wheel-crossing obstacle on two sides.
Figure 14. (a) Pitch angle comparison curve of wheel-crossing obstacle on two sides; (b) roll angle comparison curve of wheel-crossing obstacle on two sides.
Mathematics 10 03391 g014
Table 1. Main parameters of the vehicle.
Table 1. Main parameters of the vehicle.
ParameterSymbolValueUnit
Sprung massma125,649Kg
Boom massma24249Kg
Unsprung massmw650Kg
Stiffness of suspension springks130,000N/m
Damping of suspensioncs38,500N/m
Stiffness coefficient of tirekt1.9 × 106N/m
Equivalent pitching moment of inertiaIφ78,593Kg m2
Equivalent rolling moment of inertiaIφ26,506Kg m2
Wheel treadl2.15m
Distance of front axle and mass center of chassisa3.65m
Distance of middle axle and mass center of chassisb1.05m
Distance of rear axle and mass center of chassisc2.70m
Distance of hinge of boom and mass center of chassisd3.12m
Luffing angle of boomδ85°
Rotating angle of boomγ0°
Length of boomB18m
Distance of hinge point of boom and mass center of boomB18.4m
Angle between jet reaction force and horizontal directionα20°
Angle between jet reaction force and vertical directionβ45°
Area of fire monitor outlets3.11 × 10−3m2
Table 2. Root-mean-square value of vehicle attitude under stochastic signals.
Table 2. Root-mean-square value of vehicle attitude under stochastic signals.
PassiveADRCMADRC
vertical displacement (mm)2.6381.140(↓56.78%)0.786(↓70.20%)
pitch angle0.0770.048(↓37.84%)0.016(↓79.22%)
roll angle0.4770.205(↓56.92%)0.107(↓77.57%)
Table 3. Peak–peak value of vehicle attitude under sinusoidal signals.
Table 3. Peak–peak value of vehicle attitude under sinusoidal signals.
Vertical Displacement/mmPitch Angle/°Roll Angle/°
Conventional model11.170.290.68
Integrated model8.960.230.37
Difference2.210.060.31
Performance improvement19.8%20.7%45.6%
Table 4. Main parameters of the tested vehicle.
Table 4. Main parameters of the tested vehicle.
ParametersValue
Vehicle typeJP32G
Weight/kg29,898
Chassis typeXPD36
Number of axles3
Number of wheels6
Number of drive wheels6
Tire specification445/95R25
Wheelbase of first and second axle /mm4700
Wheelbase of second and third axle /mm1650
Wheel tread /mm2150
Table 5. Peak–peak value of vehicle attitude.
Table 5. Peak–peak value of vehicle attitude.
Pitch Angle/°Roll Angle/°
One sideTwo sidesOne sideTwo sides
Passive suspension0.731.851.680.55
Integrated spray-active suspension0.651.411.290.46
Difference0.080.440.390.08
Performance improvement10.9%23.7%23.2%16.3%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhu, J.; Zhao, D.; Liu, S.; Zhang, Z.; Liu, G.; Chang, J. Integrated Control of Spray System and Active Suspension Systems Based on Model-Assisted Active Disturbance Rejection Control Algorithm. Mathematics 2022, 10, 3391. https://0-doi-org.brum.beds.ac.uk/10.3390/math10183391

AMA Style

Zhu J, Zhao D, Liu S, Zhang Z, Liu G, Chang J. Integrated Control of Spray System and Active Suspension Systems Based on Model-Assisted Active Disturbance Rejection Control Algorithm. Mathematics. 2022; 10(18):3391. https://0-doi-org.brum.beds.ac.uk/10.3390/math10183391

Chicago/Turabian Style

Zhu, Jianxu, Dingxuan Zhao, Shuang Liu, Zilong Zhang, Guangyu Liu, and Jinming Chang. 2022. "Integrated Control of Spray System and Active Suspension Systems Based on Model-Assisted Active Disturbance Rejection Control Algorithm" Mathematics 10, no. 18: 3391. https://0-doi-org.brum.beds.ac.uk/10.3390/math10183391

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