Model-Free Adaptive Control for Tank Truck Rollover Stabilization

Influenced by lateral liquid sloshing in partially filled tanks, tank vehicles are apt to encounter with rollover accidents. Due to its strong nonlinearity and loading state uncertainty, it has great challenges in tank vehicle active control. Based on the model-free adaptive control (MFAC) theory, the roll stability control problem of tank trucks with different tank shapes and liquid fill percentages is explored. First, tank trucks equipped with cylinder or elliptical cylinder tanks are modelled, and vehicle dynamics is analyzed. ,is dynamic model is used to provide I/O data in the controlled system. Next, the control objective of tank vehicle rollover stabilization is analyzed and the controlled variable is selected. Subsequently, differential braking and active front steering controller are designed by MFAC algorithm. Finally, the effectiveness of the designed controllers is verified by simulation, and difference between the controllers is analyzed. ,e controller designed by MFAC algorithm is proven to be adaptive to vehicle loading and driving states. ,e controlled system has great robustness.


Introduction
Tank vehicles are widely used in road transportation of liquid cargoes. As affected by liquid sloshing in partially filled tank, the tank vehicle usually has poor driving performance, especially rollover stability [1,2]. e active control on it is of great significance.
Model based and data driven are the two main approaches proposed for the control-system design. In the model-based method, the system dynamic model is first derived, and then, a controller is designed based on that. e tank vehicle is complex fluid-solid coupling multibody system. Its dynamics is greatly influenced by liquid sloshing, which makes it difficult to construct an accurate dynamic model [3,4]. To simplify the modelling of the tank vehicle, the equivalent mechanical model is popularly used to approximately describe liquid sloshing in partially filled tank; then, the coupling mechanism of liquid sloshing and vehicle driving could be expressed easily and the vehicle dynamic modelling could be achieved. As an essential nonlinear system, the complex dynamic model renders it difficult for many control methods based on system dynamics to design the rollover stabilization controller. Furthermore, as liquid sloshing is greatly influenced by tank shapes, the equivalent mechanical model for liquid sloshing would change with it [5][6][7][8][9][10]. By far, only tank truck equipped with the horizontal cylinder or elliptical cylinder tank was modelled [6,[9][10][11][12][13][14]. Quite few studies on the modelling of tank truck with other tank shapes have been reported. Owing to this, the control algorithm developed by model-based design strategy is poorly adaptive. It is difficult for the model-based control method to achieve wide usage in practical application.
Two common data-driven controllers, PID and fuzzy logic controller, had been studied by many researchers to achieve tank vehicle rollover and yaw stabilization. Zhao et al. designed a PID controller to obtain the additional yaw moment and applied differential braking on tires to improve vehicle roll stability. e authors used the difference between the practical and theoretical roll angle to design the PID controller [15]. Hu and Zhao analyzed the yaw stability of tractor tanker semitrailer and then conducted a control on the tractor's yaw rate. Error between the practical tractor yaw rate and its desired value, as well as its derivation change rate, was used to design the fuzzy logic controller [16]. Li used the axle lateral load transfer ratio (LTR) as the controlled variable to design active antirollover bar by PID control algorithm [17]. Zeng used active front steering, active suspension, and their combination, respectively, to restrain liquid sloshing in tanks. Liquid oscillation angle was used as the controlled variable [18]. Besides the PID and fuzzy logic control method, Acarman andÖzgüner used the sliding mode control method to regulate liquid sloshing frequency, which would greatly reduce the settling time of the vehicle driving state [19]. e PID and fuzzy logic control are classic, typically offline data-driven control method. e controller design by those methods has a simple structure. However, tank vehicle dynamics is greatly affected by tank shape, liquid fill percentage, and cargo physical characters, thus rendering it difficult for those controllers to maintain a consistent control performance. Controller parameters need to be readjusted to stabilize the system. erefore, the model-free adaptive control (MFAC) method, an online data-driven control method with good robustness and adaptability, is proposed to stabilize tank vehicle. e MFAC algorithm was proposed by Hou in 1994 [20]. It is a data-driven control method based on system online data. e algorithm performs dynamic linearization at each working point of the controlled system by the aid of I/O data, thus describing system dynamics by this incremental dynamic linear model. Based on the dynamic linear model, the adaptive controller could be designed [21,22]. e pseudopartial derivative in the dynamic linear model is insensitive to the time-varying structure and time-varying phase of the system; thus, the controller exhibits strong robustness and adaptability. Wang Wang studied the fault-tolerant control problem for a class of discrete-time systems subjected to the sensor fault. For the purpose of obtaining better control performance, a model-free adaptive fault-tolerant controller is developed by employing more past control information [23]. Jiang et al. explored the heading control problem of unmanned surface vehicles with uncertainties based on MFAC theory. e conventional MFAC algorithm was developed to solve problems existing in USVs' heading control [24]. Jiang et al. studied the heading tracking problem of the six-wheel independent drive and four-wheel independent steering unmanned ground vehicle with variable wheelbase under the influence of uncertainties [25]. Li et al. addressed a pattern-moving-based modified modelfree adaptive control (PMFAC) scheme and illustrated the convergence of its tracking error for a class of nonlinear systems with the unknown model and multithreshold quantized observations [26]. erefore, the paper aims to investigate adaptive rollover stabilization for tank trucks by the MFAC method, and its outline is as follows. In Section 1, tank trucks equipped with cylinder or elliptical cylinder tanks are modelled, and vehicle dynamics is analyzed. e dynamic model of tank trucks is used to provide I/O data in the controlled system. In Section 2, control objective is analyzed and the controlled variable is selected; then, controllers based on differential braking and active front steering strategies are designed by the MFAC method. In Section 3, the effectiveness of the designed controllers is verified by simulation, and the difference between the two controllers is analyzed. In the end, a major conclusion about this research is given.  [6,7,11], as shown in Figure 1. e movement of the elliptical pendulum should be equivalent to lateral liquid sloshing in partially filled tanks.

Dynamic Modelling and Analysis of
us, the pendulum should have the same dynamics as liquid sloshing. By the principle of dynamic equivalence, parameters of the elliptical pendulum could be modelled as the function of liquid sloshing dynamics [14]: where b p is half length of the minor axis of the elliptical pendulum, a p is half length of the major axis of the elliptical pendulum, g is acceleration of gravity, ζ is length ratio of the major and minor axis of tank cross section, m p is the pendulum ball mass, m f is the static liquid mass, m is the liquid cargo mass, ω is liquid oscillation frequency, F s and a y are the peak liquid sloshing force in a sloshing cycle and the corresponding sloshing lateral acceleration at the same time, b cg is the vertical distance between the liquid center of mass and tank lowest point when liquid-free surface is leveled, b is half length of the major axis of the tank, and b f is the vertical distance between the static liquid and tank lowest point. After the sloshing frequency and force of lateral liquid sloshing have been obtained, parameters of the elliptical pendulum can be obtained according to equation (1). By curve fitting, they can be modelled as the function of the tank shape and liquid fill percentage, that is, 2 Mathematical Problems in Engineering where Δ is liquid fill percentage, which is the ratio of liquid level height to the tank height.

Modelling of Tank Truck.
As a component of the vehicle, the movement of the tank coincides with that of the vehicle. at is to say, the tank translates along vehicle Y-axis and rotates on vehicle roll axis. Consequently, the reference frame of liquid sloshing is a noninertial coordinate system with both lateral translation and roll rotation, and liquid sloshing is affected by the movement of the vehicle.
While the elliptical pendulum model is used to equivalently describe liquid sloshing in partially filled tanks, its dynamics in the noninertial coordinate system should be derived, that is, where θ is pendulum oscillation angle. When vehicle is stationary, the value is π/2, η is dimensionless damping coefficient of lateral liquid sloshing, V is vehicle driving speed, β is vehicle body slip angle, r is the vehicle yaw rate, e 2 is the coordinate of the liquid center of mass along vehicle Xaxis, φ is the vehicle roll angle, and H f is the height from the tank center to vehicle roll axis.
Since equation (3) contains nonlinear composite terms such as _ θ 2 , _ ϕ 2 , r 2 , and _ ϕr, liquid sloshing is an essential nonlinear subsystem that is difficult to linearize near its equilibrium point.
In lateral direction, vehicle dynamics along Y-axis can be expressed by where m s and m u are vehicle sprung and unsprung mass, respectively, c is the coordinate of the center of gravity of vehicle sprung mass along X-axis, e is the coordinate of the center of gravity of vehicle unsprung mass along X-axis, h s is the vertical distance from the center of gravity of vehicle sprung mass to vehicle roll axis, H 1 is the height from the static liquid to vehicle roll axis, H is the vertical distance from pendulum ball to vehicle roll axis, which is expressed by H � H f − b p sin θ, and F f and F r are lateral tire-road forces of front and rear axles, respectively. Also, in lateral direction, vehicle dynamics about Z-axis can be expressed by where I zs , I zu , and I zc are the yaw moment of inertia of vehicle sprung mass, vehicle unsprung mass and liquid cargo around their centroid coordinate system, respectively, I xzs is the product of inertia of vehicle sprung mass around its centroid coordinate system, and l f and l r are the distance from the center of gravity of the vehicle to front and rear axles, respectively.

Mathematical Problems in Engineering
Vehicle roll moment about X-axis can be expressed by where I xs and I xc are the roll moment of inertia of vehicle sprung mass and liquid cargo around their centroid coordinate system, respectively, k φ is suspension angular stiffness, and c φ is suspension angular damping. Equations (4)-(6) describe three degrees of freedom of tank trucks. Along with equation (3) that describes liquid sloshing dynamics, the tank truck is dynamically modelled. Owing to nonlinear terms existing in the model, the tank truck is an essential nonlinear system that is difficult to approximately linearize near its equilibrium point.
In the state space form, the tank truck's dynamic model is where vehicle state vector X is [θ, _ θ, β, r, ϕ, _ ϕ] T . Magic formula is used to describe tire cornering characteristics, that is, where D � −0.0004m 2 te + 8.9012m te + 163.94, m te is tire vertical load, B � 8.4; C � 1.59, and α is tire slip angle.
Considering lateral load transfer, vertical loads of the four tires are expressed by where m FL , m FR , m RL , and m RR are the vertical loads of left front tire, right front tire, left rear tire, and right rear tire, respectively, m sf , m sr , m uf, m ur , m ff , m fr , m pf , and m pr are the load distribution of vehicle sprung mass, vehicle unsprung mass, static liquid mass, and pendulum mass on the front and rear axles, respectively, d f and d r are track width of the front and rear axle, respectively, and h is the distance from tank lowest point to vehicle roll axis.

Tank Truck Dynamic Analysis.
A widely used 6 × 4 tank truck is selected to analyze its driving performance. e tank equipped on the vehicle has length of 5.8 m and crosssectional area of 2.5 m 2 . Liquid cargo density is assumed to be 1000 kg/m 3 . Other basic parameters of the vehicle are listed in Table 1.
ree elliptical cylinder tanks that have the same crosssectional area have been designed. Length ratio of major and minor axes of the tank changes from 1 to 2 with a 0.5 step. Tank dimensions are listed in Table 2. Under the same liquid fill percentage, cargo in the three tanks has the same weight.
A normal truck (labelled as NT) with the same basic parameters as the tank truck is used as a reference. An angular step test is conducted on the four vehicles. Vehicle driving speed is 25 km/h, the step steering wheel angle is 180 degrees (steering system gear ratio is about 25 : 1), and liquid fill percentage is 0.5 in the test. Dynamic response of LTAB10, LTAB15, LTAB20, and NT are shown in Figure 2. Figure 2 shows that tank trucks have slightly poor lateral stability and rather poor roll stability compared with normal truck. e yaw rate and lateral acceleration of tank trucks decrease about 4.5%∼8.2% compared to those of normal trucks, and the roll angle and LTR decrease about 40%∼ 200%. A simple conclusion can be obtained that tankers are prone to nontripping rollover accidents.

Control Objectives and Controlled
Variables. LTR and roll angel are variables that can directly illustrate vehicle roll motion. Compared with vehicle roll angle, LTR has a threshold limit. It reaches or almost reaches 1.0, which means a rollover accident would happen. us, LTR could be used as a warning signal. e rollover controller is inactive when vehicle LTR is below the set threshold. Otherwise, the rollover controller is activated. Given the consideration that rollover accident will happen as long as LTR reaches 1.0 and the active control is virtually useless, 0.8 is used as the set threshold.
LTR is a comprehensive variable that could not be measured directly, which makes it not a good choice as the controlled variable in the controlled system. It is known that lateral acceleration is the main reason which leads to vehicle roll motion. Under the action of lateral acceleration, the vehicle turns around roll axis. According to Figure 3, the roll dynamics under vehicle steady state can be written as where ΔF is axle lateral load transfer, ϕ s is vehicle roll angle under steady state, and Y ps is the coordinate of pendulum ball along vehicle Y-axis under vehicle steady state. Under vehicle steady state, vehicle sideslip angular velocity approximates to 0, lateral acceleration is mainly determined by the vehicle yaw rate. us, in equation (11), lateral acceleration of vehicle sprung mass, the static liquid, and the pendulum ball could be expressed by When LTR reaches the set threshold, equation (11) can be rewritten as Lateral acceleration obtained by equation (12) will be the threshold to keep vehicle roll stability. While the practical lateral acceleration exceeds the limitation, a rollover accident would happen.
In equation (12), the static liquid mass m f and its vertical position described by H 1 , the mass of pendulum ball m p and its position described by Y ps and H change along with tank shape and cargo fill percentage. erefore, lateral acceleration threshold for the tank truck with different tank shapes and cargo fill percentages differs. By equation (12), lateral acceleration threshold for LTAB10, LTAB15, and LTAB20 under different liquid fill percentages can be obtained, and the results are shown in Figure 4.

Mathematical Problems in Engineering
Since lateral acceleration under vehicle steady-stable state is mainly caused by yaw motion, the limitation on vehicle lateral acceleration turns into the limitation on the vehicle yaw rate. e target yaw rate of the tank trucks can be obtained by r � a y /V. e vehicle yaw rate could be measured or estimated directly, and it will be the controlled variable in the controlled system. Consequently, the architecture of tank trucks' rollover control is illustrated in Figure 5.

Control
Strategy. Analysis in Section 2.1 tells that tank vehicle roll stabilization can be achieved by limiting the vehicle yaw rate. Differential braking and active front steering are strategies widely used to control vehicle yaw motion. While differential braking is used to control vehicle yaw motion, the closed-loop system can be expressed by where ΔM z is the additional yaw moment achieved by differential braking on left or/and right tire. Active front steering applies an additional front-wheel steering angle on the basis of driver steering wheel input to achieve yaw rate control. e closed-loop system under this control strategy can be expressed by where ΔF f (δ) is tire lateral force of the front axle caused by the additional front-wheel steering angle. While the vehicle yaw rate is selected as the controlled variable, output of the controlled system is

Controller Design.
e basic principle of the MFAC method is establishing a dynamic linear model that is equivalent to a nonlinear system at every operating point and using the I/O data of the controlled system to estimate the pseudopartial derivative in the incremental dynamic linear model online. Subsequently, the weighted one-step forward controller is designed to realize the data-driven MFAC of nonlinear systems. e dynamic linearization techniques include the compact form dynamic linearization (CFDL), the partial form dynamic linearization (PFDL), and the full form dynamic linearization (FFDL). e CFDL and PFDL are used for the SISO system. ey assume that the partial derivative of the controlled system output is only related to its input.
e FFDL are suitable both for the SISO and MIMO system. It assumes that the partial derivative of the controlled system output is related to its input and its previous output.
In this research, control input of the tank truck is the additional yaw moment or the additional front-wheel steering angle, and control output of the tank truck is the vehicle yaw rate. Different from the model-driven controller design, the controller designed by the MFAC method only needs to describe the dynamic relationship between control input and output; other vehicle driving state variables do not need to be measured or estimated. is would get rid of the online estimation of pendulum swing state variables which are extremely difficult to estimate and greatly alleviate controller design difficulty.
Equations (13) and (14) show that both input and output of the controlled system have affection on the vehicle yaw rate. us, the controller is designed based on FFDL. It is assumed that the control input in time of [k − n u + 1, k] and control output in time of [k − n y + 1, k] have great influence on the vehicle yaw rate at k + 1. On this hypothesis, the linear dynamic model of the controlled system will be where ϕ n y ,n u (k) � [ϕ 1 (k), ϕ 2 (k), . . . , ϕ n y (k), ϕ n y +1 (k), . . . , ϕ n y +n u (k)] T is the pseudopartial derivative matrix of the controlled system, ΔH n y ,n u (k) � [Δr(k), . . . , Δr(k − n y + 1), Δu(k), . . . , Δu(k − n u + 1)] T is the column vector composed by input and output of the controlled system that has influence on system output at the next moment, And n y and n u are control input and output length constant of linearization, respectively. ϕ n y ,n u (k) in equation (16) is a time-varying matrix, each element in this matrix will be estimated online using the system I/O data. With the aid of this time-varying linear dynamic model, the relation between input and output of the controlled system can be obtained. Mathematical Problems in Engineering e estimation result of pseudopartial derivative matrix obtained by system I/O data is ϕ n y ,n u (k) � ϕ n y ,n u (k − 1) + η 1 ΔH n y ,n u (k − 1) μ + ΔH n y ,n u (k − 1) where η 1 is the step size factor, η 1 ⊆ (0, 1], and µ is the penalty factor, μ > 0. e controller designed on the basis of the linear dynamic model can be expressed as

Mathematical Problems in Engineering
where r * is the limit threshold of vehicle yaw rate, ϕ n y +1 is the estimated value of the (n y + 1)th submatrix of the pseudopartial derivative matrix ϕ n y ,n u , ϕ i is the estimated value of the ith submatrix of the pseudopartial derivative matrix ϕ n y ,n u , ρ is the step size factor, ρ ⊆ (0, 1], and λ is the weighting factor, λ > 0.

e Performance of Differential Braking Control Strategy. e actual additional yaw moment input under differential braking strategy is
where U is the actual additional yaw moment input, u is the theoretical control input calculated by equations (16)- (18), and M is the amplification coefficient. It has been studied in [6,7] that tankers are the most prone to encounter with rollover accident when liquid fill percentage is 0.6∼0.8. us, liquid fill percentage of 0.6 is chosen as the cargo loading state for LTAB10. With vehicle driving speed of 15 m/s, the rear axle LTR of LTAB10 quickly reaches 1.0 while the tanker is given a step input of the frontwheel steering angle of 0.07 rad (Figure 6), and the rollover accident happens. After applying the differential braking that contributes to an additional yaw moment U, tanker's rear axle LTR shall not reach to the rollover threshold of 1.0. Furthermore, by the aid of differential braking, LTAB10 can still keep roll stability while the front-wheel steering angle increased to 0.4 rad or vehicle driving speed increased to 25 m/s. Under the three vehicle driving conditions, the differential braking controller helps LTAB10 keep roll stability. Under active control, the peak and steady-state value of rear axle LTR are smaller than 0.89 and 0.75, respectively. Lateral liquid sloshing is also restrained by the controller. mentioned that the required additional yaw moment increases as vehicle instability aggravates.
Keep the initial pseudopartial derivative matrix, the step size factor, the weight factor, and the penalty factor of the controller constant, and only adjust the limit threshold of tanker yaw rate according to Figure 4. e liquid fill percentage of 0.7 and 0.8 are chosen as the cargo loading state for LTAB15 and LTAB20, respectively. Dynamic response of LTAB15 and LTAB20 with and without control is shown in Figures 7 and 8. In Figure 7, LTAB15 encounters with rollover accident under angular step test with vehicle driving speed of 15 m/s and the front-wheel steering angle of 0.05 rad. In Figure 8, LTAB20 encounters with rollover accident under the angular step test with vehicle driving speed of 15 m/s and the front-wheel steering angle of 0.04 rad. Figures 6-8 tell that the differential braking controller based on the MFAC method could work in a wide range of factors that have great influence on tanker roll stability, such as tank shape, liquid fill percentage, vehicle driving speed, and front-wheel steering angle. After applying differential braking control, the peak and steady-state value of rear axle LTR is limited below 0.96 and 0.75, respectively. Furthermore, lateral liquid sloshing is also restrained.

e Performance of Active Front Steering Control Strategy.
e actual additional front steering angle input under active front steering strategy can be obtained by equations (16)- (18). Dynamic response of tank vehicles LATB10, LTAB15, and LTAB20 with and without the active front steering controller is shown in Figures 9-11. Liquid fill percentage for these three tankers is 0.6, 0.7, and 0.8, respectively.
Under different test conditions, the initial pseudopartial derivative matrix, the step size factor, the weight factor, and the penalty factor of the controller keep constant, only the limit threshold of the vehicle yaw rate is adjusted. e simulation shows that the active front steering controller based on the MAFC method plays a positive role on tanker roll stability, and it is adaptive to tank shape, liquid cargo fill percentage, vehicle driving speed, and front-wheel steering angle. By the aid of the active controller, the peak and steadystate value of rear axle LTR are limited below 0.91 and 0.75, respectively. Pendulum oscillation is also restrained. differences. Transient characters of vehicle rear axle LTR, such as the settling time, overshoot, steady-state value, and damped oscillation frequency, are calculated for this work. e results are plotted in Figure 12. e damped oscillation frequency and the steady-state value of LTR under these two control strategies have quite slightly difference, but settling time and overshoot differ a lot. LTR under active front steering control has a much smaller settling time and overshoot than that under differential braking control, and the differences are about 54% and 37%, respectively. e comparison shows that active front steering controller is much more suitable for tanker roll stability.

Discussion and Conclusions
Tank trucks are apt to encounter with rollover accident due to lateral liquid sloshing in partially filled tanks, and active rollover control would bring positive influence on it. Given the consideration of the essential nonlinear characters of the tank truck, the model-driven controller design would not be a good choice. Not only is there a giant difficulty in the controller design, but the designed controller will be complex and fragile due to the nonlinear items. Literature review shows two common data-driven controllers, PID and fuzzy logic controllers, have been studied in tank vehicle rollover stabilization. ese two data-driven controllers are designed offline. Controller parameters keep constant and cannot be tuned adaptively. us, they are generally used in the linear controlled system. Tank vehicle dynamics reveal nonlinear characters, and it is greatly influenced by vehicle loading state (tank shape, liquid fill percentage, and physical characters of liquid cargo are included) and vehicle driving state (vehicle driving speed and the front-wheel steering angle are included). Once vehicle loading or driving state changes, the offline tuned PID or fuzzy logic controller will not be suitable. e changeable vehicle loading and driving states calls for an adaptive data-driven controller. erefore, MFAC, one of data-driven methods, is used to design the desired controller in this research. e MFAC method has the advantage of adaptability under system working condition, which meets the needs of the tanker controller design as tank trucks usually have different tank shapes, liquid cargo fill percentage, vehicle driving speed, and front-wheel steering angle. e differential braking and active front steering controller based on the MFAC method are considered in this research. In the controlled system, the vehicle yaw rate is used as the output of the controlled system, and the additional yaw moment or the additional front-wheel steering angle is used as the control input of the controlled system.
Simulation on the controlled system shows that both the differential braking and active front steering controller plays a positive role on tanker roll stability. LTR of vehicle rear axle is used as the variable to reflect vehicle roll stability. With the aid of controller, the peak and steady-state values of LTR are limited below 0.96 and 0.75, respectively. Pendulum oscillation is effectively restrained, too. e simulation also shows the controller's adaptability to tank shape, liquid fill percentage, vehicle driving speed, and front-wheel steering angle.
A comparison between the two controllers is carried out to figure out the differences. e result shows that the damped oscillation frequency and the steady-state value of LTR under the two control strategies have quite slight difference, but settling time and overshoot differ a lot. e active front steering strategy controller is much more suitable for tank roll stability, as settling time and overshoot of LTR under active front steering control is about 54% and 37% smaller than that under differential braking control.
Data Availability e tank vehicle information and modelling and controller design data used to support the findings of this study are included within the article.