The integration of internal leakage fault detection and tolerant control for single-rod hydraulic actuators is present in this paper. Fault detection is a potential technique to provide efficient condition monitoring and/or preventive maintenance, and fault tolerant control is a critical method to improve the safety and reliability of hydraulic servo systems. Based on quadratic Lyapunov functions, a performance-oriented fault detection method is proposed, which has a simple structure and is prone to implement in practice. The main feature is that, when a prescribed performance index is satisfied (even a slight fault has occurred), there is no fault alarmed; otherwise (i.e., a severe fault has occurred), the fault is detected and then a fault tolerant controller is activated. The proposed tolerant controller, which is based on the parameter adaptive methodology, is also prone to realize, and the learning mechanism is simple since only the internal leakage is considered in parameter adaptation and thus the persistent exciting (PE) condition is easily satisfied. After the activation of the fault tolerant controller, the control performance is gradually recovered. Simulation results on a hydraulic servo system with both abrupt and incipient internal leakage fault demonstrate the effectiveness of the proposed fault detection and tolerant control method.
1. Introduction
Hydraulic systems have been used in industry in a wide number of applications, including robotics and manipulators [1], machine tools [2], active suspension systems [3, 4], hydraulic positioning systems [5–8], and hydraulic load simulators [9–11], by virtue of their small size-to-power ratios, high response, high stiffness, and high load capability. Due to their ever increasing use in industry and elsewhere, it is essential that such systems should be made more reliable and safe in operation. The dynamic characteristics of hydraulic systems, however, suffer various faults [12] due to wear, battle damage, and/or unexpected failures of the system components. Therefore, fault detection and diagnosis (FDD) and fault tolerant control (FTC) of hydraulic systems have received more and more attentions.
FDD is a potential technique to provide efficient condition monitoring and/or preventive maintenance for hydraulic systems but is also a very challenging task, as it is extremely difficult to model such systems precisely owing to their nonlinear dynamics [13] and uncertainties [14]. Condition monitoring of hydraulic systems is therefore very useful in the early detection of component failure which would lead to better operational safety and economy. This has led to the increasing trend towards integrating elements of FDD as part of a control system design [15]. In order to develop FDD algorithms for hydraulic systems, a number of different approaches have been proposed in the literature. These include methods based on hardware redundancy [16], parameter estimation methods which utilize linearized models [17] and Volterra models [18], robust observer based methods using nonlinear system models [19, 20], and an informative statistical study [21]. In order to account for the presence of parametric uncertainties, Extended Kalman Filters [22] and adaptive observers [23] which estimate both the states and parameters of the system have found application in fault detection for hydraulic systems. To consider both the parametric uncertainties and the uncertain nonlinearities in hydraulic servo systems, Garimella and Yao [24] proposed a model based detection technique based on a nonlinear adaptive robust observer (ARO) [25]. State and parameter estimation are both utilized for fault detection in those papers with online PE condition monitoring.
On the other hand, fault tolerant control for electrohydraulic systems has received relatively less attention so far. But, with the increase in demand for robustness and reliability, it is fast becoming an active area of research. In [18], the authors present a fault accommodation (FA) technique based on the linear model of the system. Using the FDD results presented in [24], Gayaka and Yao [26] present a comprehensive scheme for FA which is based on adaptive robust philosophy [27], taking into account various uncertainties present in the hydraulic systems.
The internal leakage is a typical fault in hydraulic systems, which is difficult to measure directly, and will decrease the loop gain and increase the effective damping [28], thus degrading the performance of hydraulic systems. Fault detection methods for internal leakage based on Extended Kalman Filter [29] and wavelet analysis techniques [30] have been developed for double-rod hydraulic actuators. Robust passive fault tolerant control based on quantitative feedback theory (QFT) techniques [31, 32] has been widely investigated, but as pointed in [33], passive FTC is often used to save time until a more sophisticated active algorithm enters into the system. In addition, some performance has to be given up by the passive FTC scheme.
For single-rod hydraulic actuators, integrating fault detection and active fault tolerant control for internal leakage fault, which will be focused on in this paper, remains sparse. In contrast to the double-rod hydraulic actuator, the areas of the two chambers of a single-rod hydraulic actuator are different. As a result, the two dynamic equations relating the pressure changes in the two chambers to the servovalve opening cannot be combined into a single equation that relates the load pressure to the valve opening. This complicates the fault detection and tolerant controller design. In this paper, a normal nonlinear robust backstepping controller is firstly developed based on the nominal nonlinear system model, and a prescribed performance index is derived. In the fault detection module, the derived performance index is used as the fault threshold to detect the fault. That is to say, under normal condition, that is, no fault occur, the system performance is dominated by the derived performance index; if a slight leakage has occurred, the designed normal robust controller can cover this fault and the prescribed performance index is also respected. In this case, the control performance is still guaranteed though a leakage fault has occurred; thus the fault is not necessary to be alarmed and detected, the system continues to function under the normal robust controller; if a severe leakage has happened and the prescribed performance index is thus violated, then the fault is detected. Due to only internal leakage fault being considered in this paper, the fault diagnosis module is omitted, and then a fault tolerant controller is activated. The active fault tolerant controller is developed based on the adaptive control philosophy, as the internal leakage can be modeled as a parametric variation. The fault tolerant controller only needs to change a parameter in the normal robust controller and thus has a simple structure and can be easily implemented. The learning process is operated based on the PE condition which is monitored online to guarantee the excellent convergence of parameter adaptation. After the activation of the FTC controller, the system performance is gradually recovered and the leakage level can be provided for system maintenance.
Extensive simulation results have been obtained for the motion control of a hydraulic system. These results verify the effectiveness of the proposed fault detection and the fault tolerant controller.
The remainder of this paper is organized as follows. Problem formulation and system model are present in Section 2. Normal robust controller, fault detection module, and fault tolerant controller are designed in Section 3. The effectiveness of the proposed fault detection and tolerant control scheme is demonstrated in Section 4, by simulation results. Section 5 gives some conclusions.
2. Problem Formulation and Dynamic Models2.1. Nominal System Models
The system under consideration is depicted in Figure 1. The goal is to have the inertia load to track any specified motion trajectory as closely as possible in spite of various model uncertainties and/or internal leakage fault.
Schematic diagram of a single-rod hydraulic actuator system.
The dynamics of the inertia load can be described by
(1)my¨=P1A1-P2A2-F(y,y˙,t),F(y,y˙,t)=Ff+Fe,
where m and y represent the mass and the displacement of the load, respectively, P1 and P2 are the pressures inside the two chambers of the cylinder, A1 and A2 are the ram areas of the two chambers, and F is the lumped effect of uncertain nonlinearities such as friction Ff and external disturbance Fe (such as cutting force in machining, hinge moment in aircraft actuation systems). While there have been many friction models proposed [34], a simple and often adequate approach is to regard the friction force as a static nonlinear function of the velocity, which is given by
(2)Ff(y˙)=By˙+Fs(y˙),
where B represents the combined coefficient of the modeled damping and viscous friction forces and Fs is the nonlinear term that can be modeled as [34]
(3)Fs(y˙)=[fc+(fs-fc)e-|y˙/y˙v|σ]sgn(y˙),
where fc, fs represent the level of Coulomb friction and stiction friction, respectively, y˙v, σ are empirical parameters used to describe the Stribeck effect, and sgn(·) is the sign function.
The pressure dynamics in the two chambers can be written as [13]
(4)P˙1=β1V1(-A1y˙-qleak+Q1),P˙2=β2V2(A2y˙+qleak-Q2),
where V1=V01+A1y and V2=V02-A2y are the total control volume of the first and second chamber, respectively, V01 and V02 are the original total control volume of the two cylinder chambers, respectively, β1 and β2 are the effective bulk modulus of the two cylinder chambers, respectively, and qleak is the total internal leakage of the actuator due to pressure and can be modeled as follows in the normal case [13]:
(5)qleak=ClPL,
where Cl is the normal coefficient of the internal leakage, PL=P1-P2 is the load pressure of the hydraulic cylinder, Q1 is the supplied flow rate to the forward chamber, and Q2 is the return flow rate of the return chamber. Q1 and Q2 are related to the spool valve displacement of the servovalve, xv, by [13]
(6)Q1=kq1xv[s(xv)Ps-P1+s(-xv)P1-Pr],Q2=kq2xv[s(xv)P2-Pr+s(-xv)Ps-P2],
where
(7)kq1=Cdw12ρ,kq2=Cdw22ρ,s(*) is defined as
(8)s(*)={1,if*≥00,if*<0,
where Cd is the discharge coefficient, w1 and w2 are the spool valve area gradients, ρ is the density of oil, Ps is the supply pressure of the fluid, and Pr is the return pressure. As discussed in some works in the literature, the spool valve displacement xv can be related to the control input u by a first-order system [35] given by
(9)τvx˙v=-xv+kiu,
where τv and ki are the time constant and gain of the servovalve, respectively, or a second-order system [31] given by
(10)x¨v=-ωv2xv-2ζvωvx˙v+kiωv2u,
where ωv and ζv represent the servovalve natural frequency and damping ratio, respectively.
Thus, the nominal system models can be represented by the inertia load dynamics (1)~(3), pressure dynamics (4)~(8), and the servovalve dynamic (9) or (10).
2.2. Internal Leakage Fault Models
If an internal leakage fault has occurred (i.e., the actuator piston seal has been defective), as suggested by Thompson et al. [36] and widely used in [28, 31, 32], the rate of fluid flow across a faulty actuator piston seal can be modeled as a turbulent orifice flow. Combined with the normal internal leakage, the qleak in (5) can be expressed by
(11)qleak=ClPL+η(t-Tf)Ct|PL|sgn(PL),
where Ct is the leakage coefficient used to set the severity of the internal leakage fault in the following simulation and η(t-Tf) represents the time profile with Tf being the time of occurrence of the fault and is written in the following form [37]:
(12)η(t-Tf)={0ift<Tf1-e-μ(t-Tf)ift≥Tf,
where μ>0 denotes the rate of fault evolution. Small value of μ characterizes slowly developing fault, also known as incipient fault. For large value of μ, the time profile η approaches a step function, which models the abrupt fault.
Given the desired motion trajectory yd(t)=x1d(t), our objective is to synthesize a control input u such that the output y tracks yd(t) as closely as possible in spite of model uncertainties and/or internal leakage fault. Meanwhile, a fault detection module is designed to detect the severe internal leakage fault and a fault tolerant controller is synthesized to recover the control performance gradually in the presence of the severe fault.
For the desired motion trajectory, we have the following practical assumption.
Assumption 1.
The desired position yd(t) is C3 continuous and bounded.
3. Fault Detection and Controller Design3.1. Synthesizing Design Model and Issues to Be Addressed
For the controller design, a suitable design model can simplify the controller’s implementation; thus the dynamics of the inertia load can be rewritten as
(13)my¨=P1A1-P2A2-By˙-AfSf(y˙)-d~,
where d~≜Fe+Fs-AfSf(y˙) is the lumped unmodeled uncertain nonlinearity and AfSf is used to approximate the nonlinear Coulomb friction, in which Af is the amplitude of the Coulomb friction and the Sf is a continuous shape function to capture the switching characteristics of the Coulomb friction with respect to the system velocity.
Although the effects of servovalve dynamics have been included in some works in the literature, this requires an additional sensor to obtain the spool position and only minimal performance improvement can be achieved for positioning tracking, so many researchers neglect servovalve dynamics [6]. Since a high-response servovalve is used here, it is assumed that the control applied to the servovalve is directly proportional to the spool position; then the following equation is given: xv=kiu. Thus, from (8), s(xv)=s(u).
Therefore (6) can be transformed to
(14)Q1=g1uR1,Q2=g2uR2,
where g1=kq1ki, g2=kq2ki, and
(15)R1=s(u)Ps-P1+s(-u)P1-Pr,R2=s(u)P2-Pr+s(-u)Ps-P2.
Assumption 2.
For simplification, the servovalve used here is symmetric and then g=g1=g2=kq1ki=kq2ki. Since high pressure supply is used here and so the effective bulk modulus in the two chambers are large and change very small, assume βe=β1=β2. In a practical hydraulic system under normal working condition, P1 and P2 are both bounded by 0<Pr<P1<Ps, 0<Pr<P2<Ps [6].
Define the state variables as x=[x1,x2,x3,x4]T≜[y,y˙,P1,P2]T; the entire design model including (13), (4), (11), (12), (14), and (15) can be expressed in a state-space form as follows:
(16)x˙1=x2,x˙2=1m[A1x3-A2x4-Bx2-AfSf(x2)-d~],x˙3=βeV1[-A1x2-Cl(x3-x4)-η(t-Tf)Ct|x3-x4|FFFF×sgn(x3-x4)+gR1u-A1x2-Cl(x3-x4)-η(t-Tf)Ct|x3-x4|],x˙4=βeV2[A2x2+Cl(x3-x4)+η(t-Tf)Ct|x3-x4|DDD×sgn(x3-x4)-gR2uA2x2+Cl(x3-x4)+η(t-Tf)Ct|x3-x4|].
In order to make system (16) fall into the classic class of system known as strict feedback form to use the backstepping method [38] appropriately, define a new state variable as x-3=A1x3-A2x4 [5, 6, 14]; then, the state-space equation (16) is transformed as
(17)x˙1=x2,x˙2=1m[x-3-Bx2-AfSf(x2)-d~],x-˙3=f1(x)u-f2(x)-η(t-Tf)Ctf3(x),
where
(18)f1(x)=(A1R1V1+A2R2V2)βeg,f2(x)=A1V1βe[A1x2+Cl(x3-x4)]+A2V2βe[A2x2+Cl(x3-x4)],f3(x)=(A1V1+A2V2)βe|x3-x4|sgn(x3-x4).
In this paper, the full state feedback is employed and Ps, Pr are measured for nonlinear control. In general, the system is subjected to parametric uncertainties due to the variations of m, B, βe, Cl, Ct, and so on. Many adaptive controllers for hydraulic systems have been proposed to tackle these parameter uncertainties; however, we only consider the parametric uncertainty due to the internal leakage fault, that is, the parameter Ct for better fault detection and tolerant control. In many applications, for example, aircraft actuation systems, machine tools, testing instruments, and so on, this is not a strong assumption as the other parameters can be identified off-line, and the identification errors can be lumped into uncertain nonlinearities. Define the unknown parameter θ=Ct; thus the state-space equation (17) is transformed to
(19)x˙1=x2,x˙2=1m[x-3-Bx2-AfSf(x2)-d~],x-˙3=f1(x)u-f2(x)-η(t-Tf)θf3(x).
For most applications, the extent of the internal leakage fault and uncertain nonlinearities are known. Thus the following practical assumption is made [5, 14].
Assumption 3.
Parametric uncertainty and uncertain nonlinearities satisfy(20a)θ∈Ωθ≜{θ:θmin≤θ≤θmax},(20b)|d~|≤δd(t),where θmin, θmax, and δd(t) are known.
Since x1 is the displacement of the cylinder, noting (15) and Assumption 2, the following inequation always holds in hydraulic control systems
(21)f1(x)>0,∀x.
3.2. Discontinuous Projection Mapping
Let θ^ denote the estimate of unknown parameter θ and θ~ the estimation error (i.e., θ~=θ^-θ). Viewing (20a), a discontinuous projection can be defined as [14]
(22)Projθ^(•)={0ifθ^=θmaxand•>0,0ifθ^=θminand•<0,•otherwise
for an adaptation law given by θ^˙=Proj(Γτ) in which Γ>0 is an adaptation gain and τ is an adaptation function given later, and the projection mapping in (22) guarantees
(23)(P1)θ^∈Ωθ^≜{θ^:θmin≤θ^≤θmax},(P2)θ~T(Γ-1Projθ^(Γτ)-τ)≤0,∀τ.
Property (P1) implies that the parameter estimation is always within the known bounded set. Property (P2) enables one to show that the use of projection modification to the traditional adaptation law holds the perfect learning capability of the traditional one.
3.3. Normal Robust Controller Design and Results
The design parallels the backstepping procedure [38] with a robust design consideration.
Step 1.
Noting that the first equation of (19) does not have any uncertainties, a Lyapunov function can be constructed for the first two equations of (19) directly. Define a switching-function-like quantity as
(24)z2=z˙1+k1z1=x2-x2eq,x2eq≜x˙1d-k1z1,
where z1=x1-x1d(t) is the output tracking error and k1 is a positive feedback gain. Since Gs(s)=z1(s)/z2(s)=1/(s+k1) is a stable transfer function, making z1 small or converging to zero is equivalent to making z2 small or converging to zero. So the rest of the design is to make z2 as small as possible with a guaranteed transient performance. Differentiating (24) and noting (19), we have
(25)mz˙2=mx˙2-mx˙2eq=x-3-Bx2-AfSf(x2)-d~-mx˙2eq.
In this step, x-3 is treated as virtual control input. Then we can construct a robust control function α2(x1,x2,t) for the virtual control input x-3 such that output tracking error z1 converges to zero or a small value with a guaranteed transient performance. The resulting control function α2(x1,x2,t) is given by
(26)α2(x1,x2,t)=α2a+α2s,α2a(x1,x2,t)=Bx2+AfSf(x2)+mx˙2eq,α2s=α2s1+α2s2,α2s1=-k2s1z2,
where k2s1>0 is a constant gain.
In (26), α2a functions as a model compensation and α2s as a robust control law, in which α2s2 is chosen to satisfy the following conditions:
(27)Condition(i)z2[α2s2-d~]≤ε2,Condition(ii)z2α2s2≤0,
where ε2 is a positive design parameter which can be arbitrarily small. Essentially, Condition (i) of (27) represents the fact that α2s2 is synthesized to dominate the uncertain nonlinearities d~ with control accuracy measured by the design parameter ε2, and Condition (ii) is to make sure that α2s2 is dissipating in nature so that it does not interface with the functionality of the model compensation part α2a. How to choose α2s2 to satisfy constraints like (27) can be found in [27].
Remark 4.
One example of a smooth α2s2 satisfying (27) can be found in the following way. Let h2 be any smooth function satisfying
(28)h2≥δd2;
then α2s2 can be chosen as
(29)α2s2=-k2s2z2≜-h24ε2z2,
where k2s2>0 is a nonlinear gain. It can be found that (27) is satisfied.
Let z3=x-3-α2 denote the input discrepancy. For the positive-semi-definite (p.s.d.) Lyapunov function V2 defined by V2=mz22/2, noting (25), the time derivative of V2 is
(30)V˙2=z2[z3+α2-Bx2-AfSf(x2)-d~-mx˙2eq]
with the virtual control input (26); then
(31)V˙2=z2z3-k2s1z22+z2(α2s2-d~).
Step 2.
This step is to synthesize an actual control law for u such that x-3 tracks the virtual control function α2 with a guaranteed transient performance as follows. From (19), we can obtain
(32)z˙3=x-˙3-α˙2=f1(x)u-f2(x)-η(t-Tf)θf3(x)-α˙2,
where
(33)α˙2=∂α2∂t+∂α2∂x1x2+∂α2∂x2x˙2.
In (33), the acceleration x˙2 is assumed to be available to avoid the propagation of the uncertain nonlinearities d~, and thus the conservatism of the robust controller design can be reduced. In the lack of acceleration information case, one can refer to [14], the same general result form (the following Theorem 5) can also be obtained with an enhanced robust controller.
Based on (32), the final control law can be synthesized by
(34)u=ua+us,ua=1f1(x)[f2(x)+α˙2+θ^f3(x)-z2],us=1f1(x)(us1+us2)us1=-k3s1z3,
where k3s1>0 is a constant gain and θ^=0 in (34) for the case where there is no fault or only a slight internal leakage fault.
In (34), ua functions as a model compensation and us as a robust control law, in which us2 is chosen as
(35)us2=-k3s2z3≜-h34ε3z3,
where ε3 is a positive design parameter which can be arbitrarily small and h3 is a continuous function satisfying
(36)h3≥δq2(t)
in which δq(t) is a positive function for the tolerance level of the slight internal leakage fault; that is, the effect of the slight internal leakage qf=η(t-Tf)θf3(x) satisfies
(37)|qf|≤δq.
Then the nonlinear robust control law us2 (35) satisfies the following conditions:
(38)Condition(i)z3[us2-qf]≤ε3,Condition(ii)z3us2≤0.
Essentially, Condition (i) of (38) represents the fact that us2 is synthesized to tolerate the effect of the slight internal leakage qf with a control accuracy measured by the design parameter ε3, and Condition (ii) is to make sure that us2 is dissipating in nature so that it does not interface with the functionality of the model compensation part ua. Define a Lyapunov function as follows:
(39)V3=V2+12z32.
With the robust control law (34), the time derivative of V3 is
(40)V˙3=V˙2+z3z˙3=-k2s1z22+z2(α2s2-d~)-k3s1z32+z3(us2-qf).
Thus, we have the following performance theorem.
Theorem 5.
In the presence of uncertain nonlinearities satisfying (20b) and/or slight internal leakage fault satisfying (37), the robust control law (34) guarantees the following.
In general, all signals are bounded. Furthermore, the p.s.d. Lyapunov function V3 is bounded by
(41)V3≤Vr(t),
where
(42)Vr(t)=exp(-λt)V3(0)+ελ[1-exp(-λt)],
in which λ=2min{k2s1/m,k3s1} and ε=ε2+ε3.
If after a finite time t0, d~=0,qf=0, that is, in the absence of uncertain nonlinearities and internal leakage fault (i.e., d~=0,qf=0∀t≥t0), then, in addition to results in (a), asymptotic output tracking is also achieved; that is, z1→0 as t→∞.
Proof of Theorem 5.
From (40), and noting Condition (i) of (27) and (38), then
(43)V˙3≤-k2s1z22+ε2-k3s1z32+ε3≤-λV3+ε
which leads to (41). Thus z1, z2 and z3 are bounded. Noting Assumption 1 and (20a) and (20b), it follows that x2eq and the time derivative of x2eq are bounded; thus α2 is bounded, which will lead to the boundedness of x-3. Following Assumption 1, and noting (13), the acceleration x˙2 is bounded; thus α˙2 is bounded. With Assumption 2, we see that the state x is bounded. The control input u is thus bounded. This proves (a) of Theorem 5.
Now consider the situation in (b) of Theorem 5, from (40) and noting Condition (ii) of (27) and (38); then
(44)V˙3=-k2s1z22-k3s1z32≤-λV3.
Since all signals are bounded, following the standard proof procedure in ARC (more details can be found in [39, 40]), it is easy to check that z2→0 as t→∞, which leads to (b) of Theorem 5.
Remark 6.
Results of Theorem 5 indicate that the proposed normal robust controller has an exponentially converging transient performance with the exponentially converging rate λ and the final tracking error being able to be adjusted via certain controller parameters freely in a known form; it is seen from (42) that λ can be made arbitrarily large, and ε/λ, the bound of V3(∞) (an index for the final tracking errors), can be made arbitrarily small by increasing gains k1, k2s1, k3s1, and/or decreasing controller parameters ε2, ε3. In this sense, Vr(t) can be thought as the prescribed performance index and the actual performance index V3(t) is always governed by the prescribed performance index Vr(t) in the normal case and/or the presence of the slight internal leakage fault case. Such a guaranteed transient performance is especially important for the control of electrohydraulic systems since execute time of a run is very short. (b) of Theorem 5 implies that an improved performance can be obtained if there are no fault and modeling errors.
Remark 7.
Designing an extra robust controller us2 (35) is to emphasize the fact that, if the prescribed control performance Vr(t) is maintained even though in the presence of some internal leakage flow, it is not necessary to break off the function of hydraulic systems. And thus abundant unnecessary maintenance cost is saved and the economy is improved for many applications, for instance, machine tools and testing instruments. For aircraft system, this concept can improve the mission availability rate. On the other hand, the effect of the off-line identification errors for some parameters (i.e., βe, Cl, and so on) also can be though as a slight internal leakage fault, and the robust controller can guarantee the control performance with these identification errors.
3.4. Fault Detection and Tolerant Control Design
A performance-oriented fault detection scheme is proposed in this section which is used to detect severe internal leakage fault reliably for preventive maintenance and further fault accommodation (i.e., FTC). The severe internal leakage fault is defined as the violation of the inequality (37) for a period time and simultaneously the control performance is thus deteriorated heavily. The fault detection scheme consists of three components: (1) observed-error generation that indexes the control performance of the controlled hydraulic system; (2) the evaluation scheme that monitors the observed error to detect the presence of any off-nominal system behavior (i.e., decision whether a fault has occurred); and (3) residual generation that gives reliable information for the next fault diagnosis. These three components are detailed as follows.
(1) Observed-Error Generation. The selected observed error must be related to the control performance for the performance-oriented fault detection. It is clear from the normal robust controller design procedure that the p.s.d. Lyapunov function V3(t) is a suitable candidate for this purpose since it is a quadratic form with respect to those control performance indexes z1, z2, and z3. Thus, we use the performance index V3(t) defined in (39) as the observed error for the fault detection scheme.
(2) The Evaluation Scheme. In order to detect the severe internal leakage fault, we continuously monitor the observed-error V3(t) and therefore the threshold for V3(t) has to be so chosen such that, in the absence of the severe fault, V3(t) should be less than the selected threshold for all possible normal cases and/or the slight internal leakage cases; that is, the robustness has to be guaranteed to void any false alarms. From Theorem 5, we recommend the prescribed performance index Vr(t) to complete this mission and the robustness of this selected threshold is guaranteed by Theorem 5. The detection decision scheme is given as follows:
(45)Im={0ifV3(t)-Vr(t)≤0,1ifV3(t)-Vr(t)>0,
where Im is the flag whether a fault has occurred or not, Im=0 declares that the system is fault-free, and Im=1 declares that a fault has been detected.
(3) Residual Generation. The generation of the residual signal r(t) defined below is used for further fault diagnosis,
(46)r(t)={0fault-freeV3(t)otherwise.
Since only the internal leakage fault is considered here, the fault diagnosis module is omitted.
At the same time of the detection of a fault, the fault tolerant control designed later is activated immediately.
The overall block diagram of the fault detection scheme, fault tolerant control, and the flow of information in the system are shown in Figure 2. As shown in Figure 2, the proposed fault tolerant control is based on the parameter adaptation scheme. According to the value of Im, the learning mechanism is redefined as follows by using the discontinuous projection (22):
(47)θ^˙={0ifIm=0Proj(Γτ)ifIm=1.
The diagram of the proposed fault detection and tolerant control scheme.
From the properties of the discontinuous projection (22), that is, (23), the learning mechanism is controlled and the parameter estimation is thus retained in the prior known convex set Ωθ.
The tolerant controller is the same as the aforementioned normal robust controller (34), in which θ^ is updated by the parameter adaptation law (47) with the initial value θ^(0)=0. The schematic diagram of the proposed fault tolerant controller is shown in Figure 2.
After the activation of the fault tolerant controller, we have the following theorem.
Theorem 8.
With the projection type adaptation law (47) and adaptation function τ=φz3, in which φ=-f3(x), the proposed fault tolerant controller can gradually recover the control performance; that is, the results of Theorem 5 are restored as t→∞, if the following persistent excitation (PE) condition is satisfied:
(48)∫tt+Tφ2(υ)dυ≥κ,forsomeκ>0andT>0,∀t.
Proof of Theorem 8.
From the fault model (11) and (12), η(t-Tf)Ct→Ct as t-Tf→∞, thus the third equation of the system model (19) can be written as
(49)x-˙3=f1(x)u-f2(x)-θf3(x)
and the time derivative of z3 can be rewritten as
(50)z˙3=x-˙3-α˙2=f1(x)u-f2(x)-θf3(x)-α˙2.
After the activation of the fault tolerant controller, applying the control input u (34) with the adaptation law (47), then
(51)z˙3=-z2-k3s1z3+us2-θ~φ.
It is clear from (51) that there is only parametric uncertainty θ~ and no unmodeled disturbances. Therefore, it is well known that the parameter estimation converges to its true value (i.e., θ~→0 as t→∞) if the PE condition (48) is satisfied [41]. Then the effect of the fault is counteracted by the added term in the fault tolerant controller (i.e., the adaptive compensation term θ^f3(x)). That is to say, as t→∞, the time derivative of V3(t) satisfies
(52)V˙3≤-k2s1z22+ε2-k3s1z32+z3us2≤-λV3+ε
which will lead to the result (a) of Theorem 5. If there are no uncertain nonlinearities in the entire system, (i.e., d~=0), from (52) and the proof procedure of Theorem 5, we have
(53)V˙3≤-λV3
which will lead to the result (b) of Theorem 5. Thus the results of Theorem 5 are gradually recovered.
Remark 9.
It is clear from (48) that the PE condition is easily satisfied as the integrated term is a square of a scalar function φ. Thus the demand of the richness of the motion trajectory is greatly relaxed and the parameter adaptation is prone to implement. The proposed tolerant controller is also prone to realize as only a slight modification has to be made based on the normal robust controller (34). Furthermore, the fault leakage level can be caught by the excellent convergence of the parameter adaptation scheme.
Remark 10.
Based on the proposed fault detection scheme (45) and the fault tolerant controller (34) with (47), one can assign the required control performance based on the results of Theorem 5, and the tolerance level of the internal leakage (i.e., δq) via certain controller parameters freely in a known form. Small value of δq means the high sensitivity of the internal leakage fault. Large value of δq, on the contrary, means the high robustness of the internal leakage fault.
4. Simulation Results
To illustrate the above designs, simulation results are obtained for the hydraulic system discussed in Section 2, that is, (1)~(8), and (10). Since the double-rod hydraulic actuators are widely discussed [28, 31, 32], for more conveniently and easily comprehending the basis idea of the proposed fault detection and tolerant control scheme, in this section, the double-rod hydraulic actuator, which is a particular case of the single-rod hydraulic actuator, is used for the illustrative purpose. Another reason of choosing the double-rod hydraulic actuator as the plant to be discussed is that we only have a dual-vane hydraulic rotary actuator which is equivalent to the double-rod hydraulic actuator in the schematic sense in our laboratory. Moreover, it is clear from Section 3 that the type of the hydraulic actuator does not influence the verification of the effectiveness of the proposed fault detection and control scheme. The nominal equivalent parameters of the hydraulic actuator are summarized in Table 1 which have been identified in detail.
Nominal equivalent parameters of the considered hydraulic system.
Symbol
Parameter
Nominal value
m
The total mass of the load
40 kg
A1
The ram area of the cylinder
1.92×10-4 m2
A2
The ram area of the cylinder
1.92×10-4 m2
B
The coefficient of the viscous friction
45 N·s/m
fc
The Coulomb friction
52 N
fs
The stiction friction
78 N
y˙v
The Stribeck velocity
0.013 m/s
σ
The empirical parameter of the Stribeck effect
2
βe=β1=β2
The bulk modulus of the oil in the cylinder
200 Mpa
V01
The original total control volume of the cylinder
1.76×10-4 m3
V02
The original total control volume of the cylinder
1.76×10-4 m3
Cl
The normal coefficient of the internal leakage
8.9×10-12 m5/Ns
kq=kq1=kq2
The flow gain of the servovalve
9.92×10-6 m4/(s·A·N)
Ps
The supply pressure
7 Mpa
Pr
The return pressure
0 Mpa
ωv
The servovalve natural frequency
628 rad/s
ζv
The servovalve damping ratio
0.7
ki
The gain of the servovalve
4×10-3 A/V
The following two controllers are compared.
NRC: the nonlinear robust controller (34) without fault accommodation, that is, θ^≡0 no matter the severe leakage fault occurs or not. The controller parameters are listed in Table 2.
FTC: the proposed fault tolerant controller (34) with adaptation law (47), in which θmin=0, θmax=6×10-7, Γ=1×10-21. From the parameters listed in Table 2, it can be seen that the converging rate λ=100. Choosing δd=7200 means that about 1 L/min internal leakage fault flow can be tolerant in the following tests.
The parameters of the proposed control scheme.
Symbol
Parameter
Value
Af
The amplitude of the Coulomb friction
52 N
Sf
The continuous shape function
arctan (900y˙)
k1
The feedback gain for z1
10
k2s1
The feedback gain for z2
2000
δd
The bound of the uncertain nonlinearities d~
100 N
ε2
The controller parameter for α2s2
1
k3s1
The feedback gain for z3
50
δq
The tolerance level of the internal leakage fault
7200
ε3
The controller parameter for us2
1×105
In the following simulation tests, x1d=arctan(0.2sin(2πt))[1-exp(-0.1t3)]m is used as the motion trajectory which satisfies Assumption 1, and thus V3(0)=0.
To test the effectiveness of the proposed fault detection and tolerant control algorithms, three typical cases of the internal leakage fault are considered: Case 1 characterizes the abrupt internal leakage fault; Case 2 characterizes the incipient internal leakage fault; Case 3 characterizes the slight internal leakage fault. The parameters of the internal leakage fault are listed in Table 3.
The parameters of the internal leakage fault.
Case 1
Case 2
Case 3
Ct (m4/(s·N))
4.5×10-7
1×10-7
1.2×10-8
μ
10
0.5
5
Tf (s)
6
6
6
For Case 1, the fault detection is shown in Figure 3 (the index V3 is chopped off). It is clear that the fault is immediately detected (Im=1 when t=6 s) when it occurs (Tf=6 s) as the inequality (41) is immediately violated by the abrupt severe internal leakage fault. After the detection of the abrupt fault, the FTC controller is activated and the control performance is gradually restored which is shown in Figure 4. If one compares FTC with NRC, it is seen that, before the detection of the abrupt fault, the FTC and NRC are the same. Along with the activation of the parameter adaptation, the FTC can gradually counteract the effect of the fault and the serious performance deterioration of NRC is saved by FTC. The excellent convergence of the parameter estimation is given in Figure 5 which can indicate the level of the fault. The results of the control input are present in Figure 6, which show that the control input of FTC and NRC is close to the maximum authority of the control input (10 V) due to the severe internal leakage fault. In this case, the maximum flow of the internal leakage is 38 L/min, which is 60% of the load flow of the servovalve (63 L/min).
Case 1: fault detection.
Case 1: tracking errors of FTC and NRC.
Case 1: parameter estimation.
Case 1: the control input of FTC and NRC.
For Case 2, that is, the incipient internal leakage fault which develops very slowly, the fault detection is shown in Figure 7. The fault is detected (Im=1 when t=6.33 s) after a period of its occurrence (Tf= 6 s) as the fault is slowly developing and the effect of the fault between t=6 s and t=6.33 s is tolerant by the prior given bound δq. After the detection of the incipient fault, the FTC controller is activated and the control performance is gradually restored which is shown in Figure 8. The parameter estimation is given in Figure 9 and the level of the incipient fault is indicated. The results of the control input are present in Figure 10. Although the control effort of NRC and FTC is almost the same, the FTC controller makes a more accurate accommodation than that of the NRC controller. From the above results, it can be seen that the proposed scheme can effectively detect the fault which will destroy the control performance, that is, inequality (41). In the early developing period, the fault is small enough such that it does not influence the control performance; hence, it is not necessary to detect it. But once it develops largely enough and begins to deteriorate the tracking performance, the proposed scheme detected immediately and activated the FTC to compensate the fault and recover the control performance.
Case 2: fault detection.
Case 2: tracking errors of FTC and NRC.
Case 2: parameter estimation.
Case 2: control input of FTC and NRC.
For Case 3, the effect of the slight internal leakage fault does not exceed the threshold δq; that is, the inequality of (37) is respected which is shown in Figure 11. The control performance of NRC does not obviously deteriorate which is shown in Figure 12. In this case, the fault is very small such that the control performance is not affected. Hence, the fault is not necessary to be alarmed. So no fault is detected (Im=0 for all t) though the slight fault has occurred (Tf=6 s), which is given in Figure 13. In this case, the FTC controller is the same as the NRC controller. The control input is present in Figure 14. The results of Case 3 verify the prior design of the tolerance of the slight fault, that is, the bound δd. The unnecessary maintenance and the breakoff of hydraulic systems can be avoided for this slight fault. Following this design concept, the maintenance cost is reduced and the economy is improved; simultaneously the control performance is also guaranteed, that is, just the purpose of the proposed FDD scheme. The idea of the tolerance of the slight fault might make a better trade-off between the performance and the economy of the hydraulic control systems.
Case 3: the effect of the slight internal leakage fault.
Case 3: tracking errors of NRC.
Case 3: fault detection.
Case 3: the control input of NRC.
5. Conclusions
In this paper, a novel fault detection method is proposed which is a performance-oriented scheme. The main feature of the scheme is that a tolerance level of the slight fault is given, and the control performance of the hydraulic systems is guaranteed by a normal robust controller in the normal case and/or the presence of the slight fault case. For the severe fault, the detection scheme gives a reliable detection. The concept of the proposed fault detection may reduce the maintenance cost and the economy can be improved with the hydraulic control systems. Based on this fault detection scheme, the internal leakage fault is accommodated by the proposed fault tolerant controller via the parameter estimation. Since only the coefficient of the internal leakage fault is considered and other parameters are identified off-line and the identification errors can be lumped into the unmodeled nonlinearities, the PE condition is easily satisfied which gives an excellent convergence of the parameter estimation. Then the level of the fault can be provided. After the activation of the fault tolerant controller, the control performance is gradually restored. In addition, the fault tolerant controller has a simple structure which is based on the normal robust controller; thus the computing time is saved. Another contribution of the proposed fault detection and tolerant control scheme is that one can assign the required control performance and the tolerance level of the internal leakage via certain controller parameters freely in a known form. Three cases of the internal leakage fault, that is, the abrupt severe fault, the incipient fault, and the slight fault, are investigated to verify the effectiveness of the proposed fault detection and tolerant control strategy. From the model development, it can be seen that the synthesized model is particularly conservative, which might lead the proposed scheme to be inappropriate for physical hydraulic systems. In the future, a more general and comprehensive fault detection and tolerant control should be considered for hydraulic servo systems, with various parametric uncertainties and uncertain nonlinearities, to remove the prerequisites utilized in this paper.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The first author is grateful to Professor B. Yao for hosting his visit at Purdue University from October 2010 to October 2011 and the guidance in ARC design. This work was supported in part by the National Natural Science Foundation of China (Grant no. 51305203), in part by the Specialized Research Fund for the Doctoral Program of Higher Education (Grant no. 20133219120026), China, and in part by Jiangsu Planned Projects for Postdoctoral Research Funds (Grant no. 1302002A).
ZhuW. H.PiedboeufJ. C.Adaptive output force tracking control of hydraulic cylinders with applications to robot manipulators200512722062172-s2.0-2124443842410.1115/1.1898237YaoB.Al-MajedM.TomizukaM.High-performance robust motion control of machine tools: an adaptive robust control approach and comparative experiments19972263762-s2.0-003116363010.1109/3516.588624SunW.GaoH.YaoB.Adaptive robust vibration control of full-car active suspensions with electrohydraulic actuators201321624172422SunW.GaoH.KaynakO.Adaptive backstepping control for active suspension systems with hard constraints201318310721079YaoJ.JiaoZ.MaD.YanL.High-accuracy tracking control of hydraulic rotary actuators with modeling uncertainties2014192633641GuanC.PanS.Nonlinear adaptive robust control of single-rod electro-hydraulic actuator with unknown nonlinear parameters20081634344452-s2.0-4344911300110.1109/TCST.2007.908195YaoJ.JiaoZ.HanS.Friction compensation for low velocity control of hydraulic flight motion simulator: a simple adaptive robust approach201326381482210.1016/j.cja.2013.04.001YaoJ.JiaoZ.MaD.Extended-state-observer-based output feedback nonlinear robust control of hydraulic systems with backstepping201410.1109/TIE.2014.2304912YaoJ.JiaoZ.YaoB.Robust control for static loading of electro-hydraulic load simulator with friction compensation2012256954962KarpenkoM.SepehriN.Hardware-in-the-loop simulator for research on fault tolerant control of electrohydraulic actuators in a flight control application2009197106710772-s2.0-7034988397310.1016/j.mechatronics.2009.01.008YaoJ.JiaoZ.YaoB.ShangY.DongW.Nonlinear adaptive robust force control of hydraulic load simulator2012255766775SepehriN.KarpenkoM.AnL.KaramS.A test rig for experimentation on fault tolerant control and condition monitoring algorithms in fluid power systems: from design through implementation20052934414582-s2.0-27844510452MerrittH. E.1967New York, NY, USAWileyYaoB.BuF.ReedyJ.ChiuG. T.-C.Adaptive robust motion control of single-rod hydraulic actuators: theory and experiments20005179912-s2.0-003390388110.1109/3516.828592IsermannR.On the design and control of mechatronic systems-a survey19964314152-s2.0-003008449710.1109/41.481402WattonJ.1992Chichester, UKEllis HorwoodYuD.Fault diagnosis for a hydraulic drive system using a parameter-estimation method199759128312912-s2.0-003123826710.1016/S0967-0661(97)84367-5TanH. Z.SepehriN.Parametric fault diagnosis for electrohydraulic cylinder drive units2002491961062-s2.0-003647590710.1109/41.982253FrankP. M.Online fault detection in uncertain nonlinear systems using diagnostic observers: a survey199425122129215410.1080/00207729408949341MR1316312ZBL0822.93046HammouriH.KaboreP.OthmanS.BistonJ.Failure diagnosis and nonlinear observer. Application to a hydraulic process20023394-545547810.1016/S0016-0032(02)00027-3MR1931506ZBL1021.93011MunchhofM.2006Darmstadt, GermanyTechnischen Univestitat DarmstadtAnL.SepehriN.Hydraulic actuator circuit fault detection using extended Kalman filterProceedings of the American Control ConferenceJune 2003Denver, Colo, USA426142662-s2.0-0142138312RajamaniR.HedrickJ. K.Adaptive observers for active automotive suspensions: theory and experiment19953186932-s2.0-002926762210.1109/87.370713GarimellaP.YaoB.Model based fault detection of an electro-hydraulic cylinderProceedings of the American Control Conference (ACC'05)June 2005Portland, Ore, USA4844892-s2.0-23944497423GarimellaP.YaoB.Nonlinear adaptive robust observer design for a class of nonlinear systemsProceedings of the American Control ConferenceJune 2003Denver, Colo, USA439143962-s2.0-0142231041GayakaS.YaoB.Fault detection, identification and accommodation for an electro-hydraulic system: an adaptive robusProceedings of the 17th World Congress, International Federation of Automatic Control (IFAC'08)July 2008Seoul, Republic of Korea13815138202-s2.0-7996101779210.3182/20080706-5-KR-1001.2736YaoB.TomizukaM.Adaptive robust control of siso nonlinear systems in a semi-strict feedback form199733589390010.1016/S0005-1098(96)00222-1MR1451908ZBL0876.93083KarpenkoM.SepehriN.Robust position control of an electrohydraulic actuator with a faulty actuator piston seal200312534134232-s2.0-014198777710.1115/1.1592194AnL.SepehriN.Hydraulic actuator leakage fault detection using extended Kalman filter20056141512-s2.0-18344385701GoharriziA. Y.SepehriN.A wavelet-based approach to internal seal damage diagnosis in hydraulic actuators2010575175517632-s2.0-7795116159510.1109/TIE.2009.2032198KarpenkoM.SepehriN.Fault-tolerant control of a servohydraulic positioning system with crossport leakage20051311551612-s2.0-1234428512010.1109/TCST.2004.838570KarpenkoM.SepehriN.Quantitative fault tolerant control design for a leaking hydraulic actuator201013252372442-s2.0-7795614211710.1115/1.4001707TzafestasC. S.RigatosG. G.TzafestasS. G.Design of fault-tolerant control systems: passive and active approaches1998244528MR1693326Armstrong-HélouvryB.DupontP.De WitC. C.A survey of models, analysis tools and compensation methods for the control of machines with friction1994307108311382-s2.0-002846035510.1016/0005-1098(94)90209-7AlleyneA.LiuR.A simplified approach to force control for electro-hydraulic systems2000812134713562-s2.0-003454880110.1016/S0967-0661(00)00081-2ThompsonD. F.PruynJ. S.ShuklaA.Feedback design for robust tracking and robust stiffness in flight control actuators using a modified QFT technique199972161480149710.1080/002071799220128MR1723342ZBL0941.93547ZhangX.PolycarpouM. M.ParisiniT.A robust detection and isolation scheme for abrupt and incipient faults in nonlinear systems200247457659310.1109/9.995036MR1893516ZBL1125.93358KrsticM.KanellakopoulosI.KokotovicP. V.1995New York, NY, USAWileyYaoJ.YangG.JiaoZ.MaD.Adaptive robust motion control of direct-drive dc motors with continuous friction compensation201320131483754810.1155/2013/837548YaoJ.JiaoZ.MaD.Adaptive robust control of DC motors with extended state observer20146173630363710.1109/TIE.2013.2281165YaoB.JiangC.Advanced motion control: from classical PID to nonlinear adaptive robust controlProceedings of the 11th IEEE International Workshop on Advanced Motion Control (AMC'10)March 2010Nagaoka, Japan8158292-s2.0-7795373147110.1109/AMC.2010.5464026