# Disturbance Observer-Based Integral Backstepping Control for a Two-Tank Liquid Level System Subject to External Disturbances

Authors: Xiangxiang Meng; Haisheng Yu; Herong Wu; Tao Xu
Journal: Mathematical Problems in Engineering
 (2020)

---

## Abstract

A novel method of disturbance observer-based integral backstepping control is proposed for the two-tank liquid level system with external disturbances. The problem of external disturbances can be settled in this paper. Firstly, the mathematical model of the two-tank liquid level system is established based on fluid mechanics and principle of mass conservation. Secondly, an integral backstepping control strategy is designed in order to ensure liquid level tracking performance by making the tracking errors converge to zero in finite time. Thirdly, a disturbance observer is designed for the two-tank liquid level system with external disturbances. Finally, the validity of the proposed method is verified by simulation and experiment. By doing so, the simulation and experimental results prove that the scheme of disturbance observer-based integral backstepping control strategy can suppress external disturbances more effective than the disturbance observer-based sliding mode control method and has better dynamic and steady performance of the two-tank liquid level system.


---

## Body

## 1. Introduction

In the last decade or so, the two-tank liquid level control system has been widely used in industrial production. From the point of view of control, the liquid level system represented complex multivariable nonlinear problems, which is one of the important fields of nonlinear control theory application. The liquid level system is affected by external disturbances and the parameters are of great uncertainty, which makes the control task more complicated [1]. So, the precise control of the two-tank liquid level system is very difficult. Sliding mode control is a robust control method for uncertain systems, which maintains high robustness against various uncertainties, such as external disturbances and measurement errors. For a class of nonlinear systems with multiple input and multiple output, a second-order sliding mode control algorithm [2] was proposed, which takes the input and output of nonlinear systems as the standard form of controller design. In [3], feedback linearization and sliding mode algorithm are designed for a controller based on feedback linearization for a four-tank system. By inserting a boundary layer around the sliding surface, buffeting associated with sliding mode control can be reduced, which is far better than traditional proportional integral (PI) control. A fractional-order proportional integral-differential sliding surface sliding mode controller [4] for liquid level control in the two-tank system was proposed. In [5], for a class of uncertainty and unknown disturbance of multiple-input-multiple-output data-driven sliding mode control problem of nonlinear discrete systems, Weng and Gao used the nonparametric dynamic linearization technique and the second-order discrete sliding mode observer and proposed a discrete sliding mode control law based on the proportion integration differentiation (PID) sliding surface in order to obtain the faster transient response and the smaller steady-state tracking error. Benamor and Messaoud [6] studied the robust adaptive sliding mode control law for uncertain discrete systems with unknown time-varying time-delay input. Finally, the verification results of the control law on the actual system showed that although there was no overmodulation phenomenon, there was still chattering phenomenon. Although the advanced control algorithm can obtain better control performance in the simulation research, its structure is often complex and there are many undetermined parameters, so it is difficult to be popularized in practice, limiting its application in practice. Therefore, it is necessary to select a control strategy with relatively simple method and good control performance.For the two-tank liquid level system, a model-based backstepping controller and an adaptive backstepping controller were designed [7]. In order to illustrate the effectiveness of the adaptive inverse control design, a detailed experimental comparison is made with the proportional integral controller. A nonlinear generalized predictive control and a backstepping algorithm were applied and tested in [8, 9]. Khalili et al. [10] have developed a robust backstepping sliding mode controller for tracking control of a 2-DOF piezoelectric micro-operating system to eliminate chattering phenomenon in the control process. Combining the finite time current observer with the adaptive backstepping control scheme, a control mechanism with high cost performance and strong robustness was obtained. The results show that the control scheme can successfully estimate the unknown current, providing a possibility for the realization of the current-free sensor controller [11]. An active queue management scheme to control network congestion in combination with H∞ theory and integral backstepping technique to ensure better tracking performance and asymptotically stable of all signal probabilities in the closed-loop system is given in [12]. In [13, 14], the authors have applied nonlinear backstepping to ship control for quite a long time. And the strong coupling characteristic of underactuated system was resolved by a virtual controller [15].In practice, more and more attention has been paid to the influence of external interference on nonlinear systems and the suppression of interference. The distributed adaptive command-filtered backstepping scheme [16] based on the neural network was presented, which can ensure that the tracking error of the container reaches the desired origin neighborhood and all signals in the closed-loop system are bounded. An adaptive output feedback control problem [17] was studied for a class of uncertain nonlinear systems with input delay and disturbance. In [18–21], the general framework of the nonlinear system can be obtained by using the disturbance observer-based control technology, and the application of this method in the industrial field was illustrated. In [22], when the static feedback cannot guarantee the closed-loop stability, the disturbance observer was allowed to feedback as a dynamic system. Liu et al. [23, 24] studied the speed tracking control problem of the synchronous motor drive system under a matched and unmatched interference, a terminal sliding mode control method, and a port control Hamiltonian control method based on nonlinear disturbance observer, which was proposed to realize the speed and current tracking control of the permanent magnet synchronous motor drive system. In addition, an adaptive sliding mode control strategy based on interference observer [25] was proposed to solve the problems of multiple actuator faults, parameter uncertainty, and external disturbances in a quadrotor helicopter. Different mechanisms compensate actuator fault, parameter uncertainty, and external disturbances, respectively, adopt an adaptive scheme to adjust actuator fault and parameter uncertainty, and design a disturbance observer to attenuate external disturbance. Moreover, the multivariable disturbance observer was proposed to improve the antidisturbance performance of traditional advanced feedback control [26]. In [27], a robust control method based on the finite time disturbance observer was proposed to track the output of the three-tank system in the presence of mismatched uncertainties. Then, a robust control strategy with time delay compensation was designed for multiple-input-multiple-output processes [28] with matching uncertainty and process delay. Based on the mixed fuzzy reasoning system, artificial hydrocarbon network was used in the defuzzification stage, which was called fuzzy molecular control [29]. Gouta et al. [30] proposed an adaptive control and a generalized predictive control method of the coupled two-tank system, which are to minimize the multilevel cost function defined on the prediction layer. Aiming at the problem of four-tank control [31], a set of disturbance uncertainty suppression control laws were compared and proved, whose control requirements were usually expressed in the literature in the form of a set value sequence. The uncertainty class was defined as the union of four subclasses: unknown disturbance, parameter uncertainty, measurement error, and neglected dynamics. In addition to the above literature, the second-order sliding controller [32] was successfully used to adjust the liquid level of the two-tank coupling liquid level system and the computer simulation results show that the controller can adjust the liquid level with little difference in performance. A real-time implementation method of fuzzy coordinated classical proportional integral (PI) control scheme was proposed. Experimental results show that the controller structure track parameter changed quickly and had good performance when load disturbance and the set value changed [33]. The performance of an active disturbance rejection control method in coupling tank system control was studied [34]. The results show that compared with other controllers, the active disturbance rejection control method was effective in improving time domain measurement and suppressing interference, but it lacked experimental verification. Smida et al. [35] adopted the observation scheme combining the high-gain observer and sliding mode observer to improve the robustness of state estimation quality and to reconstruct the disturbance waveform in a better way.In this paper, mathematical model for the two-tank liquid level system is established based on the principle of hydromechanics and the principle of mass conservation. Then, an integral backstepping control method for the two-tank liquid level system and disturbance observer is developed. Furthermore, research studies of the integral backstepping control method and the system with disturbance have been carried out on the innovative experimental platform for the complex control system of four-tank NTC-I type in quantity. Simulation and experiment results prove this suggested control strategy and disturbance attenuation strategy highly effective compared with the disturbance observer-based sliding mode control strategy of the two-tank liquid level system. Moreover, the controller has fewer adjusting parameters, simpler structure, and easier implementation. And the steady-state and dynamic performance of the proposed controller are both far better than some complex algorithms listed in the references, such as [32]. The results show that the proposed method has high-blooded dynamic and steady-state performances. In practical industrial application, the factory workshop will change the given value irregularly according to production demand. However, it is not realistic to readjust the controller parameters every time when the given value is changed. Therefore, the control method proposed in this paper has a wide range of practical application prospects.


## 2. System Description and Modeling

The schematic drawing in Figure1 represents the model of a two-tank liquid level system. This system consists of two tanks, two level sensors (the level sensor is located at the top of each tank), two orifices (the orifice at the bottom of each tank), two pumps, a storage tank, and four manual valves. The two tanks have same cross sections, but the cross sections of the four manual valves are different. In this experimental device, the pump 1 feeds tank 1 and the outflow of tank 1 turns into partial input of tank 2. The pump 2 feeds tank 2, and the effluent from tank 2 is discharged into the storage tank.Figure 1
The schematic diagram of the two-tank process.It can be known from the law of conservation of mass and the time rate of change of liquid in each tank is given by(1)ddtρAihit=ρqint−ρqoutt,where ρ is density of liquid, Ai is the cross section of tank i (the unit is cm2), ai is the cross section of the outlet manual valve i, hit is the height of liquid inside tank i (the unit is cm), i∈1,2,3,4, the qint is the output flow rate of electric control valve, and qoutt is the output flow rate of the tank at the bottom (the unit is cm3/s).By the Bernoulli equation in hydrodynamics, the liquid flow velocityqoutt of flowing out of the tank at the bottom is related to the cross-sectional area ai of the manual control valve at the bottom of the tank and the liquid level height in the tank. So, the outflow velocity of the valve outlet the tank can be expressed as(2)qoutt=ai2ghit,where g is the gravitational acceleration.According to the conservation of mass principle, the differential equation of tank 1 and tank 2 can be written as follows:(3)h˙1t=a3A1q1−a1A12gh1t,h˙2t=a4A2q2+a1A22gh1t−a2A22gh2t,where q1 and q2 are the output flow of the pump 1 and pump 2, respectively.Define the system state and input as(4)xi=hi,ui=qi,where i∈1,2.Thus, the mathematical model of the two-tank liquid level system can be expressed as(5)x˙1=−a1A12gx1+a3A1u1,x˙2=a1A22gx1−a2A22gx2+a4A2u2.


## 3. Controller Design of the Two-Tank Liquid Level System

To facilitate the calculation and interpretation of the main contents in this section, the following constants are defined:(6)A=a1A12g,B=a3A1,C=a1A22g,D=a2A22g,E=a4A2.Based on the situation above, the system described by equation (5) can be rewritten as follows:(7)x˙1=−Ax1+Bu1,x˙2=Cx1−Dx2+Eu2.


### 3.1. The Integral Backstepping Controller Design and Stability Analysis

The accuracy of the control target will be quantified by the liquid level tracking errorse1 and e2 of tank 1 and tank 2, respectively. And tracking errors are as follows:(8)e1=x1−x1d,e2=x2−x2d.Then, first-order derivative of equation (8) can be written as(9)e˙1=x˙1−x˙1d,e˙2=x˙2−x˙2d.Define the first-order derivative of the errors as(10)e˙1=−k1e1−c1∫0te1dt,e˙2=−k2e2−c2∫0te2dt,where k1>0, k2>0, c1>0, and c2>0.Substitute equation (10) into equation (9) and end up with(11)x˙1=x˙1d−k1e1−c1∫0te1dt,x˙2=x˙2d−k2e2−c2∫0te2dt.By combining equations (7) and (11), the controller can be computed as(12)u1=x˙1d+Ax1−k1e1−c1∫0te1dtB,u2=x˙2d+Dx2−Cx1−k2e2−c2∫0te2dtE.Hence, the integral backstepping controller can be derived from equations (6), (8), and (12). The result is(13)u1=a12gx1/A1+x˙1d−k1x1−x1d−c1∫0tx1−x1ddta3/A1,u2=a22gx2/A2−a12gx1/A2+x˙2d−k2x2−x2d−c2∫0tx2−x2ddta4/A2.Select the Lyapunov functionV as(14)V1=12e12,V2=12e22,(15)V=V1+V2=12e12+12e22.Taking the first-order time derivative for equation (15), we can obtain(16)V˙=V˙1+V˙2=e1e˙1+e2e˙2.Then, we can apply equation (10) to equation (16) and the result of the equation (16) can be calculated as(17)V˙=e1−k1e1−c1∫0te1dt+e2−k2e2−c2∫0te2dt=−k1e12−c1∫0te12dt−k2e22−c2∫0te22dt.Obviously,V is positive definite and V˙ is negative definite, and both satisfy Lyapunov stability theorem, so system (5) is asymptotically stable.Remark 1.
Assume that there exists a scalar functionV of the state x, with continuous first-order derivatives such that(i)
Vx is positive definite(ii)
V˙x is negative definite(iii)
Vx⟶∞ as x⟶∞
Then, the equilibrium at thexd is globally asymptotically stable [36]. xd is the expected value of the liquid level.


### 3.2. Design and Stability Analysis of the Controller with the Disturbance Observer

This paper considers a class of nonlinear systems with exogenous disturbances, which can be described as(18)x˙1=−Ax1+Bu1+d1,x˙2=Cx1−Dx2+Eu2+d2,where d1 and d2 are exogenous disturbances, and they are bounded.(19)d1≤D1,d2≤D2,where D1 and D2 are known positive constants.Construct the disturbance observers as(20)x^˙1=−Ax^1+Bu1+d^1+h1x1−x^1,d^˙1=h2x1−x^1,(21)x^˙2=Cx1−Dx^2+Eu2+d^2+h3x2−x^2,d^˙2=h4x2−x^2,where h1, h2, h3, and h4 are the adjustable parameters.With the estimated disturbanced˜1 and d˜2, the disturbance estimated errors are defined as(22)x˜1=x1−x^1,d˜1=d1−d^1,x˜2=x2−x^2,d˜2=d2−d^2.By combining equations (19) and (21), the first-order derivative of the disturbance estimation errors can be calculate as (23)d˜˙1=d˙1−d^˙1=−λ1d˜1,d˜˙2=d˙2−d^˙2=−λ2d˜2,there λ1 and λ2 are both the nonlinear disturbance observer gains.Then, the observer error equation is given as(24)z˙i=Hizi.Here,(25)z1z2=x˜1d˜1,z3z4=x˜2d˜2,z˙1z˙2=x˜˙1d˜˙1,z˙3z˙4=x˜˙2d˜˙2,i=1,2.So, one can get(26)z˙1z˙2=−A2−h1001z1z2,z˙3z˙4=−D2−h2001z3z4,H1=−A2−h1001,H2=−D2−h2001.It is easy to prove that the observer is asymptotically stable with proper choice ofh1, h2, h3, and h4.Remark 2.
If the linearized system is strictly stable (i.e., if all eigenvalues ofH are strictly in the left-half complex plane), then the equilibrium point is asymptotically stable (for the actual nonlinear system) [36].
Finally, according to equations (13), (18), (20), and (21), the controller for the nonlinear system with disturbances can be obtained as follows:(27)u1=a12gx1/A1+x˙1d−k1x1−x1d−c1∫0tx1−x1ddt−d^1a3/A1,u2=a22gx2/A2−a12gx1/A2+x˙2d−k2x2−x2d−c2∫0tx2−x2ddt−d^2a4/A2.
In this case, tracking error equation (8) can be rewritten as(28)e1′=x1−x1d,e2′=x2−x2d.
Consider the rate of error change, the time derivative of (28) can be expressed as(29)e˙1′=x˙1−x˙1d,e˙2′=x˙2−x˙2d.
Choose a Lyapunov functionV′ for nonlinear system (18) as(30)V1′=12e1′2+12d˜12,V2′=12e2′2+12d˜22,(31)V′=V1′+V2′=12e1′2+12d˜12+12e2′2+12d˜22.
Then, the first-order time derivative of the selected Lyapunov function (31) is obtained as(32)V˙′=V˙1′+V˙2′=e1′e˙1′+d˜1′d˜˙1′+e2′e˙2′+d˜2′d˜˙2′.
By combining equations (10) and (32), (32) can be rewritten as(33)V˙′=e1′−k1e1′−c1∫0te1′dt+e2′−k2e2′−c2∫0te2′dt−λ1d˜12−λ2d˜22=−k1e1′2−c1∫0te1′2dt−k2e2′2−c2∫0te2′2dt−λ1d˜12−λ2d˜22.
Obviously,V′ is positive definite. In addition, from equation (33), one can know that V˙′ is negative definite. Hence, it satisfied Lyapunov stability theorem, which means system (18) is asymptotically stable.


## 3.1. The Integral Backstepping Controller Design and Stability Analysis

The accuracy of the control target will be quantified by the liquid level tracking errorse1 and e2 of tank 1 and tank 2, respectively. And tracking errors are as follows:(8)e1=x1−x1d,e2=x2−x2d.Then, first-order derivative of equation (8) can be written as(9)e˙1=x˙1−x˙1d,e˙2=x˙2−x˙2d.Define the first-order derivative of the errors as(10)e˙1=−k1e1−c1∫0te1dt,e˙2=−k2e2−c2∫0te2dt,where k1>0, k2>0, c1>0, and c2>0.Substitute equation (10) into equation (9) and end up with(11)x˙1=x˙1d−k1e1−c1∫0te1dt,x˙2=x˙2d−k2e2−c2∫0te2dt.By combining equations (7) and (11), the controller can be computed as(12)u1=x˙1d+Ax1−k1e1−c1∫0te1dtB,u2=x˙2d+Dx2−Cx1−k2e2−c2∫0te2dtE.Hence, the integral backstepping controller can be derived from equations (6), (8), and (12). The result is(13)u1=a12gx1/A1+x˙1d−k1x1−x1d−c1∫0tx1−x1ddta3/A1,u2=a22gx2/A2−a12gx1/A2+x˙2d−k2x2−x2d−c2∫0tx2−x2ddta4/A2.Select the Lyapunov functionV as(14)V1=12e12,V2=12e22,(15)V=V1+V2=12e12+12e22.Taking the first-order time derivative for equation (15), we can obtain(16)V˙=V˙1+V˙2=e1e˙1+e2e˙2.Then, we can apply equation (10) to equation (16) and the result of the equation (16) can be calculated as(17)V˙=e1−k1e1−c1∫0te1dt+e2−k2e2−c2∫0te2dt=−k1e12−c1∫0te12dt−k2e22−c2∫0te22dt.Obviously,V is positive definite and V˙ is negative definite, and both satisfy Lyapunov stability theorem, so system (5) is asymptotically stable.Remark 1.
Assume that there exists a scalar functionV of the state x, with continuous first-order derivatives such that(i)
Vx is positive definite(ii)
V˙x is negative definite(iii)
Vx⟶∞ as x⟶∞
Then, the equilibrium at thexd is globally asymptotically stable [36]. xd is the expected value of the liquid level.


## 3.2. Design and Stability Analysis of the Controller with the Disturbance Observer

This paper considers a class of nonlinear systems with exogenous disturbances, which can be described as(18)x˙1=−Ax1+Bu1+d1,x˙2=Cx1−Dx2+Eu2+d2,where d1 and d2 are exogenous disturbances, and they are bounded.(19)d1≤D1,d2≤D2,where D1 and D2 are known positive constants.Construct the disturbance observers as(20)x^˙1=−Ax^1+Bu1+d^1+h1x1−x^1,d^˙1=h2x1−x^1,(21)x^˙2=Cx1−Dx^2+Eu2+d^2+h3x2−x^2,d^˙2=h4x2−x^2,where h1, h2, h3, and h4 are the adjustable parameters.With the estimated disturbanced˜1 and d˜2, the disturbance estimated errors are defined as(22)x˜1=x1−x^1,d˜1=d1−d^1,x˜2=x2−x^2,d˜2=d2−d^2.By combining equations (19) and (21), the first-order derivative of the disturbance estimation errors can be calculate as (23)d˜˙1=d˙1−d^˙1=−λ1d˜1,d˜˙2=d˙2−d^˙2=−λ2d˜2,there λ1 and λ2 are both the nonlinear disturbance observer gains.Then, the observer error equation is given as(24)z˙i=Hizi.Here,(25)z1z2=x˜1d˜1,z3z4=x˜2d˜2,z˙1z˙2=x˜˙1d˜˙1,z˙3z˙4=x˜˙2d˜˙2,i=1,2.So, one can get(26)z˙1z˙2=−A2−h1001z1z2,z˙3z˙4=−D2−h2001z3z4,H1=−A2−h1001,H2=−D2−h2001.It is easy to prove that the observer is asymptotically stable with proper choice ofh1, h2, h3, and h4.Remark 2.
If the linearized system is strictly stable (i.e., if all eigenvalues ofH are strictly in the left-half complex plane), then the equilibrium point is asymptotically stable (for the actual nonlinear system) [36].
Finally, according to equations (13), (18), (20), and (21), the controller for the nonlinear system with disturbances can be obtained as follows:(27)u1=a12gx1/A1+x˙1d−k1x1−x1d−c1∫0tx1−x1ddt−d^1a3/A1,u2=a22gx2/A2−a12gx1/A2+x˙2d−k2x2−x2d−c2∫0tx2−x2ddt−d^2a4/A2.
In this case, tracking error equation (8) can be rewritten as(28)e1′=x1−x1d,e2′=x2−x2d.
Consider the rate of error change, the time derivative of (28) can be expressed as(29)e˙1′=x˙1−x˙1d,e˙2′=x˙2−x˙2d.
Choose a Lyapunov functionV′ for nonlinear system (18) as(30)V1′=12e1′2+12d˜12,V2′=12e2′2+12d˜22,(31)V′=V1′+V2′=12e1′2+12d˜12+12e2′2+12d˜22.
Then, the first-order time derivative of the selected Lyapunov function (31) is obtained as(32)V˙′=V˙1′+V˙2′=e1′e˙1′+d˜1′d˜˙1′+e2′e˙2′+d˜2′d˜˙2′.
By combining equations (10) and (32), (32) can be rewritten as(33)V˙′=e1′−k1e1′−c1∫0te1′dt+e2′−k2e2′−c2∫0te2′dt−λ1d˜12−λ2d˜22=−k1e1′2−c1∫0te1′2dt−k2e2′2−c2∫0te2′2dt−λ1d˜12−λ2d˜22.
Obviously,V′ is positive definite. In addition, from equation (33), one can know that V˙′ is negative definite. Hence, it satisfied Lyapunov stability theorem, which means system (18) is asymptotically stable.


## 4. Simulation Results and Analysis

In order to verify the effectiveness of the above control strategy, we compared it with the disturbance observer-based sliding mode control method.The disturbance observer-based sliding mode controller can be described as(34)u1=A1a3−m1sgns1+n1+b1x1−x1d+n1b1∫e1dt+x˙1d+a12gx1A1−d^1,u2=A2a4−m2sgns2+n2+b2x2−x2d+n2b2∫e2dt+x˙2d+a22gx2A2−a12gx1A2−d^2,where b1, b2, m1, m2, n1, and n2 are positive constants.The comparative study results of the two control strategies are carried out in Simulink environment. Then, the adjustable parameters of the two-tank coupled liquid level system are given in Table1.Table 1
Adjustable parameters of the system mode.
Parameters
Value
Unit
a
1
0.2
cm
2
a
2
0.3
cm
2
a
3
0.3
cm
2
a
4
0.3
cm
2
A
1
196
cm
2
A
2
196
cm
2The desired equilibrium points are(35)x1d=20cm,x2d=16cm.


### 4.1. The Undisturbed Simulation Results

#### 4.1.1. Sliding Mode Controller

The controlled liquid level trajectories of the two tanks are given in Figures2 and 3, respectively. The adjustable parameters are determined, as shown in Table 2.Figure 2
The liquid-level curve of tank 1.Figure 3
The liquid-level curve of tank 2.Table 2
The sliding mode controller adjustable parameters.
Parameters
Value
m
1
1
n
1
0.1
b
1
0.05
m
2
1
n
2
0.1
b
2
0.06As shown in Figures2 and 3, the dynamic response time of the liquid-level curve of tank 1 is short, its rise time is tr=10s, but overshoot occurs, its overshoot is σ%=0.0000025%, and there has always been a steady-state error of 0.0002 cm. Although the liquid-level curve of tank 2 did not overshoot, its dynamic stage rose slowly when approaching the equilibrium point, and the rise time is tr=100s.


#### 4.1.2. Integral Backstepping Controller

The controlled liquid level trajectories of the two tanks are shown in Figures2 and 3, respectively. The adjustable parameters are determined, as shown in Table 3.Table 3
The integral backstepping controller adjustable parameters.
Parameters
Value
k
1
1
c
1
0.0001
k
2
1
c
2
0.0001As shown in Figures2 and 3. The rise time of liquid-level curve of tank 1 is also tr=10s. Although overshoot occurs and overshoot is σ%=0.0000025%, the steady-state error is zero. The liquid-level curve of tank 2 did not overshoot and its dynamic response time is short, and its rise time is tr=10s.


### 4.2. Simulation Results with Disturbance

The simulation verification in this paper is obtained by injecting the disturbance in the form of step signal into the system without changing the controller parameters of the two-tank liquid level system. The adjustable parameters are determined, as shown in Table4.Table 4
The disturbance observer parameters.
Parameters
Value
h
1
10
h
2
−1
h
3
1
h
4
−0.01


#### 4.2.1. Sliding Mode Controller

The liquid-level curves and control inputs trajectories of two tanks are given in Figures4 and 5 and Figures 6 and 7, respectively. The disturbances are added to tank 1 and tank 2 at the simulation time of 100 seconds, respectively. Simulation results show that the injection of external disturbances into a single tank has no effect on the liquid level of another tank. It can be seen from Figures 8 and 9 that the recovery time is relatively long of liquid levels in tank 1 and tank 2 after injection of external disturbance, and the adjustment time is ts=50s.Figure 4
The liquid-level curve of tank 1 after adding disturbance.Figure 5
The liquid-level curve of tank 2 after adding disturbance.Figure 6
The input curves of controlleru1.Figure 7
The input curves of controlleru2.Figure 8
The liquid-level curve of tank 1 after adding disturbance.Figure 9
The liquid‐level curve of tank 2 after adding disturbance.


#### 4.2.2. Integral Backstepping Controller

The liquid-level curves and control inputs trajectories of two tanks are given in Figures10 and 11 and Figures 6 and 7, respectively. At the simulation time of 100 seconds, the disturbances are added to tank 1 and tank 2, respectively. Then, the simulation results show that the injection of external disturbances into a single tank can also achieve almost no effect on the liquid level of another tank. It can be seen from Figures 8 and 9 that the recovery time is faster of liquid levels in tank 1 and tank 2 after injection of external disturbance, and the adjustment time is ts=10s.Figure 10
The liquid-level curve of tank 1 after adding disturbance.Figure 11
The liquid-level curve of tank 2 after adding disturbance.


### 4.3. Comparison of Simulation Results

For the determination and optimization of the parameters, the Lyapunov stability theorem is used to determine the range of parameters, and then the specific parameter values are determined by the empirical trial and error method, such as in Figures12 and 13.Figure 12
The liquid-level curve of tank 1.Figure 13
The liquid-level curve of tank 2.Simulation results show that the two-tank liquid level system with sliding mode control algorithm has a slight liquid level overshoot, which not only has a certain amount of steady-state error but also takes a long time to restore stability after injection of disturbances. In the two-tank liquid level system, using the integral backstepping control algorithm based on the disturbance observer, the liquid level has a slight overshoot, but its steady-state error is zero, and it can recover to a stable value in a very short time after injecting the external disturbances. According to the above discussion, the integral backstepping control method based on the disturbance observer has good dynamic, stability, and disturbance suppression characteristics.


## 4.1. The Undisturbed Simulation Results

### 4.1.1. Sliding Mode Controller

The controlled liquid level trajectories of the two tanks are given in Figures2 and 3, respectively. The adjustable parameters are determined, as shown in Table 2.Figure 2
The liquid-level curve of tank 1.Figure 3
The liquid-level curve of tank 2.Table 2
The sliding mode controller adjustable parameters.
Parameters
Value
m
1
1
n
1
0.1
b
1
0.05
m
2
1
n
2
0.1
b
2
0.06As shown in Figures2 and 3, the dynamic response time of the liquid-level curve of tank 1 is short, its rise time is tr=10s, but overshoot occurs, its overshoot is σ%=0.0000025%, and there has always been a steady-state error of 0.0002 cm. Although the liquid-level curve of tank 2 did not overshoot, its dynamic stage rose slowly when approaching the equilibrium point, and the rise time is tr=100s.


### 4.1.2. Integral Backstepping Controller

The controlled liquid level trajectories of the two tanks are shown in Figures2 and 3, respectively. The adjustable parameters are determined, as shown in Table 3.Table 3
The integral backstepping controller adjustable parameters.
Parameters
Value
k
1
1
c
1
0.0001
k
2
1
c
2
0.0001As shown in Figures2 and 3. The rise time of liquid-level curve of tank 1 is also tr=10s. Although overshoot occurs and overshoot is σ%=0.0000025%, the steady-state error is zero. The liquid-level curve of tank 2 did not overshoot and its dynamic response time is short, and its rise time is tr=10s.


## 4.1.1. Sliding Mode Controller

The controlled liquid level trajectories of the two tanks are given in Figures2 and 3, respectively. The adjustable parameters are determined, as shown in Table 2.Figure 2
The liquid-level curve of tank 1.Figure 3
The liquid-level curve of tank 2.Table 2
The sliding mode controller adjustable parameters.
Parameters
Value
m
1
1
n
1
0.1
b
1
0.05
m
2
1
n
2
0.1
b
2
0.06As shown in Figures2 and 3, the dynamic response time of the liquid-level curve of tank 1 is short, its rise time is tr=10s, but overshoot occurs, its overshoot is σ%=0.0000025%, and there has always been a steady-state error of 0.0002 cm. Although the liquid-level curve of tank 2 did not overshoot, its dynamic stage rose slowly when approaching the equilibrium point, and the rise time is tr=100s.


## 4.1.2. Integral Backstepping Controller

The controlled liquid level trajectories of the two tanks are shown in Figures2 and 3, respectively. The adjustable parameters are determined, as shown in Table 3.Table 3
The integral backstepping controller adjustable parameters.
Parameters
Value
k
1
1
c
1
0.0001
k
2
1
c
2
0.0001As shown in Figures2 and 3. The rise time of liquid-level curve of tank 1 is also tr=10s. Although overshoot occurs and overshoot is σ%=0.0000025%, the steady-state error is zero. The liquid-level curve of tank 2 did not overshoot and its dynamic response time is short, and its rise time is tr=10s.


## 4.2. Simulation Results with Disturbance

The simulation verification in this paper is obtained by injecting the disturbance in the form of step signal into the system without changing the controller parameters of the two-tank liquid level system. The adjustable parameters are determined, as shown in Table4.Table 4
The disturbance observer parameters.
Parameters
Value
h
1
10
h
2
−1
h
3
1
h
4
−0.01


### 4.2.1. Sliding Mode Controller

The liquid-level curves and control inputs trajectories of two tanks are given in Figures4 and 5 and Figures 6 and 7, respectively. The disturbances are added to tank 1 and tank 2 at the simulation time of 100 seconds, respectively. Simulation results show that the injection of external disturbances into a single tank has no effect on the liquid level of another tank. It can be seen from Figures 8 and 9 that the recovery time is relatively long of liquid levels in tank 1 and tank 2 after injection of external disturbance, and the adjustment time is ts=50s.Figure 4
The liquid-level curve of tank 1 after adding disturbance.Figure 5
The liquid-level curve of tank 2 after adding disturbance.Figure 6
The input curves of controlleru1.Figure 7
The input curves of controlleru2.Figure 8
The liquid-level curve of tank 1 after adding disturbance.Figure 9
The liquid‐level curve of tank 2 after adding disturbance.


### 4.2.2. Integral Backstepping Controller

The liquid-level curves and control inputs trajectories of two tanks are given in Figures10 and 11 and Figures 6 and 7, respectively. At the simulation time of 100 seconds, the disturbances are added to tank 1 and tank 2, respectively. Then, the simulation results show that the injection of external disturbances into a single tank can also achieve almost no effect on the liquid level of another tank. It can be seen from Figures 8 and 9 that the recovery time is faster of liquid levels in tank 1 and tank 2 after injection of external disturbance, and the adjustment time is ts=10s.Figure 10
The liquid-level curve of tank 1 after adding disturbance.Figure 11
The liquid-level curve of tank 2 after adding disturbance.


## 4.2.1. Sliding Mode Controller

The liquid-level curves and control inputs trajectories of two tanks are given in Figures4 and 5 and Figures 6 and 7, respectively. The disturbances are added to tank 1 and tank 2 at the simulation time of 100 seconds, respectively. Simulation results show that the injection of external disturbances into a single tank has no effect on the liquid level of another tank. It can be seen from Figures 8 and 9 that the recovery time is relatively long of liquid levels in tank 1 and tank 2 after injection of external disturbance, and the adjustment time is ts=50s.Figure 4
The liquid-level curve of tank 1 after adding disturbance.Figure 5
The liquid-level curve of tank 2 after adding disturbance.Figure 6
The input curves of controlleru1.Figure 7
The input curves of controlleru2.Figure 8
The liquid-level curve of tank 1 after adding disturbance.Figure 9
The liquid‐level curve of tank 2 after adding disturbance.


## 4.2.2. Integral Backstepping Controller

The liquid-level curves and control inputs trajectories of two tanks are given in Figures10 and 11 and Figures 6 and 7, respectively. At the simulation time of 100 seconds, the disturbances are added to tank 1 and tank 2, respectively. Then, the simulation results show that the injection of external disturbances into a single tank can also achieve almost no effect on the liquid level of another tank. It can be seen from Figures 8 and 9 that the recovery time is faster of liquid levels in tank 1 and tank 2 after injection of external disturbance, and the adjustment time is ts=10s.Figure 10
The liquid-level curve of tank 1 after adding disturbance.Figure 11
The liquid-level curve of tank 2 after adding disturbance.


## 4.3. Comparison of Simulation Results

For the determination and optimization of the parameters, the Lyapunov stability theorem is used to determine the range of parameters, and then the specific parameter values are determined by the empirical trial and error method, such as in Figures12 and 13.Figure 12
The liquid-level curve of tank 1.Figure 13
The liquid-level curve of tank 2.Simulation results show that the two-tank liquid level system with sliding mode control algorithm has a slight liquid level overshoot, which not only has a certain amount of steady-state error but also takes a long time to restore stability after injection of disturbances. In the two-tank liquid level system, using the integral backstepping control algorithm based on the disturbance observer, the liquid level has a slight overshoot, but its steady-state error is zero, and it can recover to a stable value in a very short time after injecting the external disturbances. According to the above discussion, the integral backstepping control method based on the disturbance observer has good dynamic, stability, and disturbance suppression characteristics.


## 5. Experimental Results and Analysis

In order to verify the performances of the propose control strategy, experiments are developed on the Innovative experimental platform for the complex control system of the four-tank NTC-I type. As shown in Figures14 and 15, they are the experimental platform and Wincc liquid-level monitoring interface, respectively. The experimental apparatus is a typical process control object in industrial process-four tank. For the object, a feedback control system based on programmable logic controller (PLC) and Matlab/simulink are constructed. The PLC part uses the Siemens CPU s7-300 module to extend the special analog input. The module collects the actual water level of the four-capacity water tank and expands the dedicated analog output module feed water. The pump provides the actual required analog voltage. Communication between Matlab and PLC/HMl industrial control is realized through OLE for process control (OPC) communication technology. The control system is integrated with some advanced system identification techniques and advanced control methods that can be applied to the controlled object of the four-tank easily. The upper computer uses a common desktop computer, and the operating system uses Microsoft’s Windows XP system. Other development software used in this device includes WinCC V6.2, monitoring software of SIEMENS, down-bit STEP7 V 5.4 of PLC controller configuration programming and configuration communication mode, and OPC server configuration SIMATIC NET software. The adjustable parameters of the two-tank liquid level system are given in Table 1.Figure 14
Experimental platform.Figure 15
Wincc liquid-level monitoring interface.The desired equilibrium points are(36)x1d=20cm,x2d=16cm.


### 5.1. The Undisturbed Experimental Results

#### 5.1.1. Sliding Mode Controller

The controller adjustable parameters are given in Table5.Table 5
The sliding mode controller adjustable parameters.
Parameters
Value
m
1
100
n
1
10
b
1
0.01
m
2
100
n
2
10
b
2
0.01The actual controlled liquid-level trajectories of the two tanks are shown in Figure16. The dynamic response time of the experimental curves of tank 1 and tank 2 are relatively fast, rising time is tr=35s and tr=45s, respectively, but there is always obvious chattering phenomenon of the liquid level, and the amplitude is ±0.5cm.Figure 16
The liquid-level curves of sliding mode control.


#### 5.1.2. Integral Backstepping Controller

The integral backstepping controller adjustable parameters are given in Table6.Table 6
The integral backstepping controller adjustable parameters.
Parameters
Value
k
1
0.08
c
1
10
k
2
0.12
c
2
1The actual controlled liquid-level trajectories of the two tanks are given in Figure17. Under the same condition of valve opening, the dynamic response time of tank 1 and tank 2 liquid-level curves is still very fast and the rising time is tr=30s and tr=40s, respectively. Moreover, there is no chattering phenomenon in the liquid-level curve under this control method and the steady-state process is very good.Figure 17
The liquid-level curves of integral backstepping control.


### 5.2. Experimental Results with Disturbance

In this paper, disturbances are added to the two-tank liquid level system and the controller parameters are unchanged. The disturbance observer parameters are shown in Table7.Table 7
The disturbance observer parameters.
Parameters
Value
h
1
1
h
2
−0.01
h
3
1
h
4
−0.01


#### 5.2.1. Sliding Mode Controller

The actual controlled liquid-level curves and control input trajectories of the two tanks with disturbances are shown in Figures18–20. A certain amount of water is injected into tank 1 at 20:07:20. The liquid level of tank 1 has apparent fluctuation, and the adjustment time is ts=30s. At 17:26:40, a certain amount of water is injected into tank 2, the liquid level of tank 2 has apparent fluctuation, and the adjustment time is ts=30s.Figure 18
The liquid-level curves of tank 1 with disturbance on sliding mode control.Figure 19
The liquid-level curves of tank 1 with disturbance on integral backstepping control.Figure 20
The input curves of the sliding mode controller.


#### 5.2.2. Integral Backstepping Controller

The actual controlled liquid-level curves and control input trajectories of the two tanks with disturbances are given in Figures21–23. A certain amount of water is injected into tank 1 and tank 2 at 11:12:30 and 11:00:30, respectively. The liquid levels of tank 1 and tank 2 have hardly any fluctuation, and the adjustment time is ts=15s and ts=15s, respectively.Figure 21
The liquid-level curves of tank 2 with disturbance on sliding mode control.Figure 22
The liquid-level curves of tank 2 with disturbance on integral backstepping control.Figure 23
The input curves of the integral backstepping controller.As shown in Figures24 and 25, without changing the controller parameters of the two-tank liquid level system, It still has good dynamic and steady-state performance when the set value is changed.Figure 24
The liquid-level changed curves of sliding mode control.Figure 25
The liquid-level changed curves of integral backstepping control.


### 5.3. Comparison of Experimental Results

The sliding mode control algorithm is tested and verified on the equipment platform. It can be seen from Figures16, 18 and 19 that the two-tank liquid level system controlled by the sliding mode algorithm always has the phenomenon of chattering, especially when the liquid level is disturbed externally. As can be seen from Figures 17, 21 and 22, the disturbance observer-based integral backstepping control method not only has no chattering phenomenon but also exhibits good steady-state performance and disturbance-suppression characteristics in the case of disturbance.


## 5.1. The Undisturbed Experimental Results

### 5.1.1. Sliding Mode Controller

The controller adjustable parameters are given in Table5.Table 5
The sliding mode controller adjustable parameters.
Parameters
Value
m
1
100
n
1
10
b
1
0.01
m
2
100
n
2
10
b
2
0.01The actual controlled liquid-level trajectories of the two tanks are shown in Figure16. The dynamic response time of the experimental curves of tank 1 and tank 2 are relatively fast, rising time is tr=35s and tr=45s, respectively, but there is always obvious chattering phenomenon of the liquid level, and the amplitude is ±0.5cm.Figure 16
The liquid-level curves of sliding mode control.


### 5.1.2. Integral Backstepping Controller

The integral backstepping controller adjustable parameters are given in Table6.Table 6
The integral backstepping controller adjustable parameters.
Parameters
Value
k
1
0.08
c
1
10
k
2
0.12
c
2
1The actual controlled liquid-level trajectories of the two tanks are given in Figure17. Under the same condition of valve opening, the dynamic response time of tank 1 and tank 2 liquid-level curves is still very fast and the rising time is tr=30s and tr=40s, respectively. Moreover, there is no chattering phenomenon in the liquid-level curve under this control method and the steady-state process is very good.Figure 17
The liquid-level curves of integral backstepping control.


## 5.1.1. Sliding Mode Controller

The controller adjustable parameters are given in Table5.Table 5
The sliding mode controller adjustable parameters.
Parameters
Value
m
1
100
n
1
10
b
1
0.01
m
2
100
n
2
10
b
2
0.01The actual controlled liquid-level trajectories of the two tanks are shown in Figure16. The dynamic response time of the experimental curves of tank 1 and tank 2 are relatively fast, rising time is tr=35s and tr=45s, respectively, but there is always obvious chattering phenomenon of the liquid level, and the amplitude is ±0.5cm.Figure 16
The liquid-level curves of sliding mode control.


## 5.1.2. Integral Backstepping Controller

The integral backstepping controller adjustable parameters are given in Table6.Table 6
The integral backstepping controller adjustable parameters.
Parameters
Value
k
1
0.08
c
1
10
k
2
0.12
c
2
1The actual controlled liquid-level trajectories of the two tanks are given in Figure17. Under the same condition of valve opening, the dynamic response time of tank 1 and tank 2 liquid-level curves is still very fast and the rising time is tr=30s and tr=40s, respectively. Moreover, there is no chattering phenomenon in the liquid-level curve under this control method and the steady-state process is very good.Figure 17
The liquid-level curves of integral backstepping control.


## 5.2. Experimental Results with Disturbance

In this paper, disturbances are added to the two-tank liquid level system and the controller parameters are unchanged. The disturbance observer parameters are shown in Table7.Table 7
The disturbance observer parameters.
Parameters
Value
h
1
1
h
2
−0.01
h
3
1
h
4
−0.01


### 5.2.1. Sliding Mode Controller

The actual controlled liquid-level curves and control input trajectories of the two tanks with disturbances are shown in Figures18–20. A certain amount of water is injected into tank 1 at 20:07:20. The liquid level of tank 1 has apparent fluctuation, and the adjustment time is ts=30s. At 17:26:40, a certain amount of water is injected into tank 2, the liquid level of tank 2 has apparent fluctuation, and the adjustment time is ts=30s.Figure 18
The liquid-level curves of tank 1 with disturbance on sliding mode control.Figure 19
The liquid-level curves of tank 1 with disturbance on integral backstepping control.Figure 20
The input curves of the sliding mode controller.


### 5.2.2. Integral Backstepping Controller

The actual controlled liquid-level curves and control input trajectories of the two tanks with disturbances are given in Figures21–23. A certain amount of water is injected into tank 1 and tank 2 at 11:12:30 and 11:00:30, respectively. The liquid levels of tank 1 and tank 2 have hardly any fluctuation, and the adjustment time is ts=15s and ts=15s, respectively.Figure 21
The liquid-level curves of tank 2 with disturbance on sliding mode control.Figure 22
The liquid-level curves of tank 2 with disturbance on integral backstepping control.Figure 23
The input curves of the integral backstepping controller.As shown in Figures24 and 25, without changing the controller parameters of the two-tank liquid level system, It still has good dynamic and steady-state performance when the set value is changed.Figure 24
The liquid-level changed curves of sliding mode control.Figure 25
The liquid-level changed curves of integral backstepping control.


## 5.2.1. Sliding Mode Controller

The actual controlled liquid-level curves and control input trajectories of the two tanks with disturbances are shown in Figures18–20. A certain amount of water is injected into tank 1 at 20:07:20. The liquid level of tank 1 has apparent fluctuation, and the adjustment time is ts=30s. At 17:26:40, a certain amount of water is injected into tank 2, the liquid level of tank 2 has apparent fluctuation, and the adjustment time is ts=30s.Figure 18
The liquid-level curves of tank 1 with disturbance on sliding mode control.Figure 19
The liquid-level curves of tank 1 with disturbance on integral backstepping control.Figure 20
The input curves of the sliding mode controller.


## 5.2.2. Integral Backstepping Controller

The actual controlled liquid-level curves and control input trajectories of the two tanks with disturbances are given in Figures21–23. A certain amount of water is injected into tank 1 and tank 2 at 11:12:30 and 11:00:30, respectively. The liquid levels of tank 1 and tank 2 have hardly any fluctuation, and the adjustment time is ts=15s and ts=15s, respectively.Figure 21
The liquid-level curves of tank 2 with disturbance on sliding mode control.Figure 22
The liquid-level curves of tank 2 with disturbance on integral backstepping control.Figure 23
The input curves of the integral backstepping controller.As shown in Figures24 and 25, without changing the controller parameters of the two-tank liquid level system, It still has good dynamic and steady-state performance when the set value is changed.Figure 24
The liquid-level changed curves of sliding mode control.Figure 25
The liquid-level changed curves of integral backstepping control.


## 5.3. Comparison of Experimental Results

The sliding mode control algorithm is tested and verified on the equipment platform. It can be seen from Figures16, 18 and 19 that the two-tank liquid level system controlled by the sliding mode algorithm always has the phenomenon of chattering, especially when the liquid level is disturbed externally. As can be seen from Figures 17, 21 and 22, the disturbance observer-based integral backstepping control method not only has no chattering phenomenon but also exhibits good steady-state performance and disturbance-suppression characteristics in the case of disturbance.


## 6. Conclusions

In view of the two-tank liquid level system, a methodology of disturbance observer-based integral backstepping control is proposed in this paper. According to the principle of fluid mechanics and mass conservation, the mathematical model of the two-tank liquid level system is established. As can be seen from the simulation results, when the external disturbance is added to a single tank, the liquid-level value of another tank has no influence. Furthermore, it can be known from the experiment results that the given value can be set discretionarily and reached quickly. Then, the results of simulation and experiment verified that the disturbance observer-based integral backstepping control methodology has good dynamic and static performance compared with the disturbance observer-based sliding mode control method. Moreover, the steady-state and dynamic performance of the proposed controller are both far better than some complex algorithms listed in the references. Its structure is simpler, and adjustable parameters are fewer and easier to realize. Therefore, the control method proposed in this paper has a wide range of practical application prospects. However, the controller proposed in this paper does not consider the influence of time delay on the system, and it is only limited to the two-tank water level system. Whether it is suitable for the four-tank water level system needs to be further verified.Further, in addition to considering the influence of disturbance on the system, the influence of time delay on the system is also considered and corresponding solution is given, which is combined with the method proposed in this paper to achieve better control effect. Then, the composite control method is extended and applied to the four-tank liquid level system.



---

*Source: 6801205-2020-01-20.xml*