# Tracking Control of a Leg Rehabilitation Machine Driven by Pneumatic Artificial Muscles Using Composite Fuzzy Theory

Authors: Ming-Kun Chang
Journal: The Scientific World Journal
 (2014)

---

## Abstract

It is difficult to achieve excellent tracking performance for a two-joint leg rehabilitation machine driven by pneumatic artificial muscles (PAMs) because the system has a coupling effect, highly nonlinear and time-varying behavior associated with gas compression, and the nonlinear elasticity of bladder containers. This paper therefore proposes a T-S fuzzy theory with supervisory control in order to overcome the above problems. The T-S fuzzy theory decomposes the model of a nonlinear system into a set of linear subsystems. In this manner, the controller in the T-S fuzzy model is able to use simple linear control techniques to provide a systematic framework for the design of a state feedback controller. Then the LMI Toolbox of MATLAB can be employed to solve linear matrix inequalities (LMIs) in order to determine controller gains based on the Lyapunov direct method. Moreover, the supervisory control can overcome the coupling effect for a leg rehabilitation machine. Experimental results show that the proposed controller can achieve excellent tracking performance, and guarantee robustness to system parameter uncertainties.


---

## Body

## 1. Introduction

In cases of traumatic brain injury, bone injury, amputation, or spinal cord injury caused by misfortunes such as traffic accidents and cerebral apoplexy, lower limb rehabilitation machine can help patients recover extremity functions by means of continuous passive motion (CPM). Traditionally, physical therapy for achieving functional rehabilitation is carried out by medical therapists on a person-to-person basis. However, recently many automatic rehabilitation devices have been gradually applied in physical therapy programs. Rehabilitation machines are usually driven by electric motors, which are typically rigid in nature. Because of this, actuators can generate discomfort or pain when interfacing with humans. For this reason, current electromechanical actuation systems should be replaced to ensure adaptability, conformity, and safety. An adequate actuator for a rehabilitation device must provide physically adjustable compliance and safety and ensure soft contact with the patient, similar to the behavior of human muscles. It has been suggested that pneumatic artificial muscles (PAMs) can contribute towards achieving more comfortable devices for interfacing with human limb segments.PAMs behave in a manner very similar to the muscles that move the skeletons of animals and have many advantages, such as high power to weight ratio [1], high power to volume ratio [2], low maintenance, negligible mechanical wear, low cost, cleanliness, high reliability, flexibility, and compliance for use with humans. For these reasons, PAMs are commonly employed in rehabilitation engineering, nursing, and human-friendly therapeutic machine.However, PAMs exhibit highly nonlinear and time-varying behavior due to the compression of air and the nonlinear elasticity of bladder containers. This makes it difficult for classical controllers to achieve excellent control performance. In recent years, researchers have developed a wide variety of approaches to overcome these problems. Noritsugu and Tanaka [3] developed four modes of linear motion with impedance control to control force during movement and used an adaptive identification method to estimate the system model. Lilly and Yang [4] applied a sliding mode controller to a planar arm actuated by two PMA groups; simulation results were consistent with theoretical findings for two different masses. Ahn and Anh [5] adopted an ARNN controller in a PAM manipulator for reducing tracking errors. Shen [6] developed a full nonlinear model that encompassed all the major existing nonlinearities. Based on this model, the standard sliding mode control approach was applied to obtain robust control, even in the event of model uncertainties and disturbances.Since the inception of fuzzy set theory by Zadeh [7] in 1965, a great deal of research has been focused on fuzzy control systems. Takagi and Sugeno [8] proposed the T-S fuzzy model-based controller in 1985, and the T-S fuzzy model-based system subsequently emerged as one of the most active and fruitful areas of fuzzy control. Using a T-S fuzzy model-based controller, a complex dynamic model can be decomposed into a set of local linear subsystems via fuzzy inference. Stability analysis is carried out using the Lyapunov direct method, where the control problem is formulated into linear matrix inequalities (LMIs). Based on this approach, Ahn and Anh [9] also developed an inverse double nonlinear autoregressive model with exogenous control based on the T-S fuzzy model applied in a PAM robot. A novel H∞ control structure based on a Takagi-Sugeno model [10] was proposed to track the desired trajectories, and simulation results illustrated the efficiency of the proposed approach for the new rehabilitation device.The leg rehabilitation machine driven by PAMs is a two-input, two-output system. This paper proposes composite fuzzy theory, which includes T-S fuzzy tracking control and supervisor control in order to improve tracking performance. The proposed approach decomposes the model of a nonlinear system into a set of linear subsystems with associated nonlinear weighting functions, enabling the use of simple linear control techniques without the need for complicated nonlinear control strategies, and also provides a systematic framework for the design of a state feedback controller [11]. It has been shown that a composite fuzzy control system can be guaranteed to be asymptotically stable if a common positive definite solution exists for a set of Lyapunov inequalities. In addition, the supervisory control can overcome the coupling effect due to two-joint motion. In view of the above advantages, the proposed controller was applied to the output tracking control of this system, and experimental results verified that the proposed controller is capable of achieving excellent tracking performance.The remainder of the paper is organized as follows. Section2 describes the control strategies. Section 3 describes the system. In Section 4, the dynamics of the model are derived. Experimental results for output tracking are shown in Section 5. Finally, conclusions are presented in Section 6.


## 2. Control Strategies

### 2.1. Takagi-Sugeno Fuzzy Tracking Controller

Consider a general nonlinear dynamic equation(1)x˙(t)=f(x(t))+g(x(t))u(t)y(t)=q(x(t)),
where x∈Rn is the state vector, y∈Rm is the controlled output, u∈Rm is the control input vector, and f(x), g(x), and q(x) are nonlinear functions with appropriate dimensions. The nonlinear system (1) can then be expressed by the fuzzy system.Rulei:
(2)IFz1(t)isF1iand⋯andzg(t)isFgiTHENx˙(t)=Aix(t)+Biu(t),i=1,2,…,r,
where z(t)1~zg(t) are the premise variables including system states, Fij denotes the fuzzy sets, r is the number of fuzzy rules, and Ai and Bi are system matrices with appropriate dimensions. For simplicity, this study assumed that the membership functions had been normalized; that is, ∑i=1rΠj=1gFji(zj)=1. As in (1), using the singleton fuzzier, product inferred, and weighted defuzzier, the fuzzy system is inferred as
(3)x˙(t)=∑i=1rhi[Aix(t)+Biu(t)],
where hi(z(t))=Πj=1gFji(zj(t)). Note that ∑i=1rhi(z(t))=1 for all t, where ∑i=1rhi(z(t))≥0 for i=1,2,…,r are regarded as grade functions.For output tracking control, the control objective is required to satisfy(4)y-(t)-r(t)⟶0ast⟶∞,
where r(t) denotes the desired trajectory or reference signal. To convert the output tracking problem into a stabilization problem, a set of virtual desired variables xd(t) was introduced, to be tracked by the state variable x. Let x~(t)=x(t)-xd(t) denote the tracking error for the state variables. The time derivative of x~(t) yields
(5)x~˙(t)=x˙-x˙d=∑irhi[Aix(t)+Biu(t)]-x˙d(t).If the control inputu(t) is assumed to satisfy the following equation:
(6)∑i=1rhiBiτ(t)=∑i=1rhiBiu(t)+∑i=1rhiAixd(t)-x˙d(t),
where τ(t) is a new control to be designed, then the tracking error system (5) results in the following form:
(7)x~˙(t)=∑i=1rhiAix~(t)+∑i=1rhiBiτ(t).The design of the new controlτ(t) is similar to solving a stabilization problem. The purpose is to steer x~(t) to zero, which means that state x(t) tracks xd(t). The new fuzzy controller τ(t) is designed on the basis of parallel distributed compensation (PDC) and is represented as follows:Rulei:
(8)IFz1(t)isF1iand⋯andzg(t)isFgi:THENτ(t)=-Kix~(t),
where Ki represents feedback gain. The inferred output of the PDC controller is expressed in the following form:
(9)τ(t)=-∑i=1rhiKix~(t).Substituting (9) into (7) yields
(10)x~˙(t)=∑i=1r∑j=1rhihj(Ai-BiKj)x~(t).
The stability analysis of this tracking system (10) is carried out using the Lyapunov direct method, and the Lyapunov function is defined as
(11)V(x~(t))=x~T(t)Px~(t)>0,
where P is a positive symmetric matrix. Taking the derivative of V with respect to time yields
(12)V˙(x~(t))=x~˙T(t)Px~+x~T(t)Px~˙(t)V˙(x~(t))=∑i=1r∑j=1rhihj(x~T(AiT-KjTBiT)Px~)V˙(x~(t))=+∑i=1r∑j=1rhihjx~TP(Ai-BiKj)x~V˙(x~(t))=∑i=1r∑j=1rhihjx~T[(Ai-BiKi)TP+P(Ai-BiKi)]x~.
The controller is stable if V˙<0. Hence, the LMI form is expressed as follows:
(13)(Ai-BiKi)TP+P(Ai-BiKi)<0fori=1,2,…,r,GijP+PGij<01≤i<j≤r,
where Gij=(Ai-BiKj+Aj-BjKi)/2 and Gii=Ai-BiKi.The controller gainKi is obtained using the LMI toolbox of MATLAB. If there exists a common positive definite matrix P that satisfies inequalities (13), it can be guaranteed that the tracking error will approach zero.


### 2.2. Composite Fuzzy Tracking Controller

Because the leg rehabilitation machine has a coupling effect due to mechanism interaction, many fuzzy model controllers in the related literature exhibit restrictive tracking control in application. The proposed approach introduces supervisory control in order to overcome the coupling effect. Theith rule of the proposed controller is defined as follows.Rulei:
(14)IFz1(t)isF1iand⋯andzg(t)isFgiTHENu(t)=Kix~(t)+usi=1,2,…,r,
where us(t)∈Rm. The proposed controller consists of a local state feedback Kix~ and a supervisory control us. Therefore, the output of the proposed controller is
(15)u(t)=∑irhiKix~(t)+us.The closed-loop system is given by(16)x˙(t)=∑i=1r∑j=1rhihj[Ai-BiKj]x~(t)+∑i=1rhiBius(t)x˙(t)=∑i=1rhi2Giix~(t)+2∑i<jrhihjGijx~(t)+Bus(t).
Suppose that there exist a symmetric and positive definite matrix P and some matrices Ki so that the following reduced stability condition holds:
(17)(Ai-BiKi)TP+P(Ai-BiKi)≤-Qi,i=1,…,r,
where Qi is a positive definite matrix. Based on this assumption, each subsystem is locally controllable, and a stable feedback gain is obtainable. Intuitively, a common matrix P that satisfies (17) can be obtained more easily than can one that fulfills the basic stabilization conditions. When the LMI method is applied, conditions (17) can be efficiently verified. If a feasible solution is obtained, the design proceeds to exploit the supervisory control in order to deal with the coupling effects.Choose the Lyapunov function candidate,V1(x)=x~TPx~. The time derivative of V1(x) is as follows:
(18)V˙1(x)=∑i=1rhi2x~T(GiiTP+PGii)x~V˙1(x)=+2∑i<jrhihjx~T(GijTP+PGij)x~+2x~TPBusV˙1(x)≤-∑i=1rhi2x~TQix~V˙1(x)=+2∑i<jrhihjx~T(GijTP+PGij)x~+2x~TPBus.
Given the matrix property, clearly,
(19)λmin⁡(GijTP+PGij)∥x~∥2≤x~T(GijTP+PGij)x~λmin⁡(GijTP+PGij)∥x~∥2≤λmin⁡(GijTP+PGij)∥x~∥2,
where λmin⁡⁡(max⁡) denotes the smallest (largest) eigenvalue of the matrix. Define
(20)α=max⁡i,j⁡λmax⁡(GijTP+PGij)for1≤i<j≤r.
A relaxed condition concerning the coupling effect is expressed as
(21)∑i<jrhihjx~T(GijTP+PGij)x~≤k1∥x~∥2,k1=r(r-1)2α.
Finding the maximum value of ∑i<jrhihjxT(GijTP+PGij)x~ is equivalent to determining the maximum value of ∑i<jrhihjλmax⁡(GijTP+PGij). This can be presented as a nonlinear programming. The optimal algorithms are employed to seek the best solution. Moreover, the MATLAB Optimization Toolbox consists of functions that minimize or maximize general nonlinear functions. By using the toolbox, the nonlinear programming is expressed in the following form:
(22)max⁡i,j∑i<jrhihjλmax⁡(GijTP+PGij)1≤i<j≤rSubjectto∑i=1rμi=1μi≥0Subjectto∑j=1rμj=1μj≥0.
The largest eigenvalue of (GijTP+PGij) can be obtained in advance, so the maximum value is determined to be
(23)k2=max⁡i,j∑i<jhihjλmax⁡(GijTP+PGij).
The following supervisory control is chosen:
(24)us={-BTPx∥x~TPB∥2k∥x~∥2,if∥x~TPB∥≠00if∥x~TPB∥=0,
where k>kj, j=1 or 2. If ∥x~TPB∥≠0, then substituting (24) into (18) gives
(25)V˙1(x)≤-∑i=1rhi2x~TQix~+2kj∥x~∥2-2k∥x~∥2V˙1(x)≤-∑i=1rhi2x~TQix~=-V2(x),
where V2(x) is a positive definite function. When ∥x~TPB∥=0 can give the following form:
(26)x~T[(Ai-BiKi)TP+P(Ai-BiKi)]x=x~T(AiTP+PAi)x~≤-x~TQix~,i=1,…,r,∑i<jrhihjx~T(GijTP+PGij)x~=∑i<jrhihj[x~T(AiP-PAi)x~+x~T(AjTP-PAj)x~]≤-∑i<jhihjx~T(Qi+Qf)x~,1≤i<j≤r,the time derivative ofV1(x) becomes
(27)V˙1(x)≤-∑i=1rhi2x~TQix~-∑i<jhihjxT(Qi+Qf)x~=-V3(x),
where V3(x) is a positive definite function. Thus, the closed-loop fuzzy system is asymptotically stable.


## 2.1. Takagi-Sugeno Fuzzy Tracking Controller

Consider a general nonlinear dynamic equation(1)x˙(t)=f(x(t))+g(x(t))u(t)y(t)=q(x(t)),
where x∈Rn is the state vector, y∈Rm is the controlled output, u∈Rm is the control input vector, and f(x), g(x), and q(x) are nonlinear functions with appropriate dimensions. The nonlinear system (1) can then be expressed by the fuzzy system.Rulei:
(2)IFz1(t)isF1iand⋯andzg(t)isFgiTHENx˙(t)=Aix(t)+Biu(t),i=1,2,…,r,
where z(t)1~zg(t) are the premise variables including system states, Fij denotes the fuzzy sets, r is the number of fuzzy rules, and Ai and Bi are system matrices with appropriate dimensions. For simplicity, this study assumed that the membership functions had been normalized; that is, ∑i=1rΠj=1gFji(zj)=1. As in (1), using the singleton fuzzier, product inferred, and weighted defuzzier, the fuzzy system is inferred as
(3)x˙(t)=∑i=1rhi[Aix(t)+Biu(t)],
where hi(z(t))=Πj=1gFji(zj(t)). Note that ∑i=1rhi(z(t))=1 for all t, where ∑i=1rhi(z(t))≥0 for i=1,2,…,r are regarded as grade functions.For output tracking control, the control objective is required to satisfy(4)y-(t)-r(t)⟶0ast⟶∞,
where r(t) denotes the desired trajectory or reference signal. To convert the output tracking problem into a stabilization problem, a set of virtual desired variables xd(t) was introduced, to be tracked by the state variable x. Let x~(t)=x(t)-xd(t) denote the tracking error for the state variables. The time derivative of x~(t) yields
(5)x~˙(t)=x˙-x˙d=∑irhi[Aix(t)+Biu(t)]-x˙d(t).If the control inputu(t) is assumed to satisfy the following equation:
(6)∑i=1rhiBiτ(t)=∑i=1rhiBiu(t)+∑i=1rhiAixd(t)-x˙d(t),
where τ(t) is a new control to be designed, then the tracking error system (5) results in the following form:
(7)x~˙(t)=∑i=1rhiAix~(t)+∑i=1rhiBiτ(t).The design of the new controlτ(t) is similar to solving a stabilization problem. The purpose is to steer x~(t) to zero, which means that state x(t) tracks xd(t). The new fuzzy controller τ(t) is designed on the basis of parallel distributed compensation (PDC) and is represented as follows:Rulei:
(8)IFz1(t)isF1iand⋯andzg(t)isFgi:THENτ(t)=-Kix~(t),
where Ki represents feedback gain. The inferred output of the PDC controller is expressed in the following form:
(9)τ(t)=-∑i=1rhiKix~(t).Substituting (9) into (7) yields
(10)x~˙(t)=∑i=1r∑j=1rhihj(Ai-BiKj)x~(t).
The stability analysis of this tracking system (10) is carried out using the Lyapunov direct method, and the Lyapunov function is defined as
(11)V(x~(t))=x~T(t)Px~(t)>0,
where P is a positive symmetric matrix. Taking the derivative of V with respect to time yields
(12)V˙(x~(t))=x~˙T(t)Px~+x~T(t)Px~˙(t)V˙(x~(t))=∑i=1r∑j=1rhihj(x~T(AiT-KjTBiT)Px~)V˙(x~(t))=+∑i=1r∑j=1rhihjx~TP(Ai-BiKj)x~V˙(x~(t))=∑i=1r∑j=1rhihjx~T[(Ai-BiKi)TP+P(Ai-BiKi)]x~.
The controller is stable if V˙<0. Hence, the LMI form is expressed as follows:
(13)(Ai-BiKi)TP+P(Ai-BiKi)<0fori=1,2,…,r,GijP+PGij<01≤i<j≤r,
where Gij=(Ai-BiKj+Aj-BjKi)/2 and Gii=Ai-BiKi.The controller gainKi is obtained using the LMI toolbox of MATLAB. If there exists a common positive definite matrix P that satisfies inequalities (13), it can be guaranteed that the tracking error will approach zero.


## 2.2. Composite Fuzzy Tracking Controller

Because the leg rehabilitation machine has a coupling effect due to mechanism interaction, many fuzzy model controllers in the related literature exhibit restrictive tracking control in application. The proposed approach introduces supervisory control in order to overcome the coupling effect. Theith rule of the proposed controller is defined as follows.Rulei:
(14)IFz1(t)isF1iand⋯andzg(t)isFgiTHENu(t)=Kix~(t)+usi=1,2,…,r,
where us(t)∈Rm. The proposed controller consists of a local state feedback Kix~ and a supervisory control us. Therefore, the output of the proposed controller is
(15)u(t)=∑irhiKix~(t)+us.The closed-loop system is given by(16)x˙(t)=∑i=1r∑j=1rhihj[Ai-BiKj]x~(t)+∑i=1rhiBius(t)x˙(t)=∑i=1rhi2Giix~(t)+2∑i<jrhihjGijx~(t)+Bus(t).
Suppose that there exist a symmetric and positive definite matrix P and some matrices Ki so that the following reduced stability condition holds:
(17)(Ai-BiKi)TP+P(Ai-BiKi)≤-Qi,i=1,…,r,
where Qi is a positive definite matrix. Based on this assumption, each subsystem is locally controllable, and a stable feedback gain is obtainable. Intuitively, a common matrix P that satisfies (17) can be obtained more easily than can one that fulfills the basic stabilization conditions. When the LMI method is applied, conditions (17) can be efficiently verified. If a feasible solution is obtained, the design proceeds to exploit the supervisory control in order to deal with the coupling effects.Choose the Lyapunov function candidate,V1(x)=x~TPx~. The time derivative of V1(x) is as follows:
(18)V˙1(x)=∑i=1rhi2x~T(GiiTP+PGii)x~V˙1(x)=+2∑i<jrhihjx~T(GijTP+PGij)x~+2x~TPBusV˙1(x)≤-∑i=1rhi2x~TQix~V˙1(x)=+2∑i<jrhihjx~T(GijTP+PGij)x~+2x~TPBus.
Given the matrix property, clearly,
(19)λmin⁡(GijTP+PGij)∥x~∥2≤x~T(GijTP+PGij)x~λmin⁡(GijTP+PGij)∥x~∥2≤λmin⁡(GijTP+PGij)∥x~∥2,
where λmin⁡⁡(max⁡) denotes the smallest (largest) eigenvalue of the matrix. Define
(20)α=max⁡i,j⁡λmax⁡(GijTP+PGij)for1≤i<j≤r.
A relaxed condition concerning the coupling effect is expressed as
(21)∑i<jrhihjx~T(GijTP+PGij)x~≤k1∥x~∥2,k1=r(r-1)2α.
Finding the maximum value of ∑i<jrhihjxT(GijTP+PGij)x~ is equivalent to determining the maximum value of ∑i<jrhihjλmax⁡(GijTP+PGij). This can be presented as a nonlinear programming. The optimal algorithms are employed to seek the best solution. Moreover, the MATLAB Optimization Toolbox consists of functions that minimize or maximize general nonlinear functions. By using the toolbox, the nonlinear programming is expressed in the following form:
(22)max⁡i,j∑i<jrhihjλmax⁡(GijTP+PGij)1≤i<j≤rSubjectto∑i=1rμi=1μi≥0Subjectto∑j=1rμj=1μj≥0.
The largest eigenvalue of (GijTP+PGij) can be obtained in advance, so the maximum value is determined to be
(23)k2=max⁡i,j∑i<jhihjλmax⁡(GijTP+PGij).
The following supervisory control is chosen:
(24)us={-BTPx∥x~TPB∥2k∥x~∥2,if∥x~TPB∥≠00if∥x~TPB∥=0,
where k>kj, j=1 or 2. If ∥x~TPB∥≠0, then substituting (24) into (18) gives
(25)V˙1(x)≤-∑i=1rhi2x~TQix~+2kj∥x~∥2-2k∥x~∥2V˙1(x)≤-∑i=1rhi2x~TQix~=-V2(x),
where V2(x) is a positive definite function. When ∥x~TPB∥=0 can give the following form:
(26)x~T[(Ai-BiKi)TP+P(Ai-BiKi)]x=x~T(AiTP+PAi)x~≤-x~TQix~,i=1,…,r,∑i<jrhihjx~T(GijTP+PGij)x~=∑i<jrhihj[x~T(AiP-PAi)x~+x~T(AjTP-PAj)x~]≤-∑i<jhihjx~T(Qi+Qf)x~,1≤i<j≤r,the time derivative ofV1(x) becomes
(27)V˙1(x)≤-∑i=1rhi2x~TQix~-∑i<jhihjxT(Qi+Qf)x~=-V3(x),
where V3(x) is a positive definite function. Thus, the closed-loop fuzzy system is asymptotically stable.


## 3. System Descriptions

Figure1 shows the experimental setup, including four PAMs, two rotary potentiometers, four pressure proportional valves, and four pressure transducers. The hardware includes an IBM-compatible personal computer to calculate the control signal, which controls the pressure proportional valve through a D/A card. The angles of the joints are detected using rotary potentiometers, the air pressure of each PAM is measured using pressure transducers, and the measurements are then fed back to the computer through an A/D card. These specifications are listed in Table 1.Table 1
Component specifications.
Number
Component
Specifications
1, 2, 3, 4
PAM
Festo, MAS-20-150N
7, 8, 9, 10
Pressure proportional valve
Mac, PPC5C
5, 6
Rotary potentiometer
Keen Engineering, KRT2050
11, 12, 13, 14
Pressure transducer
Jihsense, SN: 911166, 0–10 kgf/cm2
16
IBM-compatible PC
Pentium 1.8 GHz
DAC, ADC
Automation, AIO3321Figure 1
The experimental setup.Figure2 presents the operation principle of the leg rehabilitation machine, depicting a two-joint leg. The behavior of the leg manipulated by the rehabilitation machine is similar to that of a human leg. Output angles θ1 and θ2 simulate the knee and ankle joints, and the ranges of the rotary angles θ1 and θ2 are from -45° to 45° and from -50° to 50°, respectively. The link mass m1=2.7kg, m2=0.81kg, and the link length l1=0.5m, l2=0.26m. The rotating torque τ is generated by the difference in pressure Δp between the two opposing PAMs. That is, when pa>pb, as in Figure 2, the torque exerted on the joint is counterclockwise and the rotation of the joint is also counterclockwise.Figure 2
Operation principle of the leg rehabilitation machine driven by PAMs.So, a pair of such PAMs is tied together around a pulley with a radiusri, as in Figure 2. Then, the torque values imparted to the pulley by the PAM pair are [12]
(28)τ1=(ϕ1a-ϕ1b)r1τ2=(ϕ2a-ϕ2b)r2,
where(29a)ϕ1a=F1a(p1a)-K1a(p1a)r1θ1-B1a(p1a)r1θ˙1(29b)ϕ1b=F1b(p1b)-K1b(p1b)r1θ1-B1b(p1b)r1θ˙1(29c)ϕ2a=F2a(p2a)-K2a(p2a)r2θ2-B2a(p2a)r2θ˙2(29d)ϕ2b=F2b(p2b)-K2b(p2b)r2θ2-B2b(p2b)r2θ˙2,where the spring coefficient K(p) and the damping coefficient B(p) are given by Reynolds et al. [13].The desired input pressuresPa=[p1ap2a]T and Pb=[p1bp2b]T for each PAM are generated by the following equation:
(30)Pa(t)=P0+ΔP(t),Pb(t)=P0-ΔP(t),
where P0=[p10p20]T is a nominal constant input PAM pressure and ΔP(t)=[Δp1Δp2]T is the control pressure input with an arbitrary function of time. Because the pressure input ΔP(t)=[Δp1Δp2]T and output θ=[θ1θ2]T, the system can be written as a two-input, two-output (TITO) control system. The control signal u=[u1u2]T is proportional to ΔP based on the pressure proportional valve’s characteristics. That is, ΔP can be used instead of u as a control input.


## 4. Dynamic Model of a Two-Joint Leg Rehabilitation Machine Driven by PAMs

Figure2 shows a two-joint leg rehabilitation machine driven by PAMs, and the dynamic equation is given as follows [14]:
(31)M(θ)θ¨+C(θ,θ˙)θ˙+G(θ)=τ,
where
(32)M(θ)=[(m1+m2)l2m2l1l2(s1s2+c1c2)m2l1l2(s1s2+c1c2)m2l22]C(θ,θ˙)=[0-m2l1l2(c1s2-s1c2)θ˙2-m2l1l2(c1s2-s1c2)θ˙10]G(θ)=[(m1+m2)l1gs1-m2l2gs2]
and M(θ) is the moment of inertia, C(θ,θ˙) includes Coriolis and centripetal force, and G(θ) is the gravitational force. Notation s1=sin(θ1), s2=sin(θ2), c1=cos⁡(θ1), and c2=cos⁡(θ2). Let x1=θ1, x2=θ˙1, x3=θ2, and x4=θ˙2; then (31) can be written as the following state-space form [14]:
(33)x˙1=x2x˙2=f1(x)+g11(x)τ1+g12τ2x˙3=x4x˙4=f2(x)+g21(x)+g22τ2,
where
(34)f1(x)=(s1c2-c1s2)[m2l1l2(s1s2+c1c2)x22-m2l22x42]l1l2[(m1+m2)-m2(s1s2+c1c2)2]+[(m1-m2)l2gs1-m2l2gs2(s1s2+c1c2)]l1l2[(m1+m2)-m2(s1s2+c1c2)2]f2(x)=(s1c2-c1s2)[-(m1+m2)l12x22+m2l1l2(s1s2+c1c2)x42]l1l2[(m1+m2)-m2(s1s2+c1c2)2]+[-(m1+m2)l1gs1(s1s2+c1c2)+(m1+m2)l1gs2]l1l2[(m1+m2)-m2(s1s2+c1c2)2]g11(x)=m2l22m2l12l22[(m1+m2)-m2(s1s2+c1c2)2]g12(x)=-m2l1l2(s1s2+c1c2)m2l12l22[(m1+m2)-m2(s1s2+c1c2)2]g21(x)=-m2l1l2(s1s2+c1c2)m2l12l22[(m1+m2)-m2(s1s2+c1c2)2]g22(x)=(m1+m2)l12m2l12l22[(m1+m2)-m2(s1s2+c1c2)2].


## 5. Experimental Studies

Figure3 shows the leg rehabilitation machine using an actual human loading with a 65 kg weight. The automatic device can help patients to recover lower limb motion function by means of continuous passive motion, such as a sinusoidal wave command, an irregular curve command, and an end-effect tracking command. The experiments include both the proposed approach and PDC for comparison in order to evaluate efficacy and control performance. The controllers were implemented on an Intel Pentium 1.8 GHz PC with a sampling time of 5 ms, and the entire control software was coded in C++.Figure 3
The two-joint leg rehabilitation device with actual human loading.This study attempts to use as few rules as possible in order to minimize design effort and complexity. The T-S fuzzy model of the system is thus given the following four-rule fuzzy model:(35)R1:IFx1isabout(π/4)andx3isabout(π/4)R1h:THENx˙=A1x1+B1uR2:IFx1isabout(π/4)andx3isabout(-π/4)R2h:THENx˙=A2x1+B2uR3:IFx1isabout(-π/4)andx3isabout(π/4)R3h:THENx˙=A3x1+B3uR4:IFx1isabout(-π/4)andx3isabout(-π/4)R4h:THENx˙=A4x1+B4u,
where
(36)A1=[01004.40820.00590.6742-0.00020001-1.45720.0002-3.97220.0002]A2=[01004.67340.00490.5320.00010001-1.2145-0.0001-3.5630.0001]A3=[01005.47820.00210.6876-0.000200011.45250.00024.67230.0002]A4=[01004.9821-0.00440.65430.00020001-1.1034-0.00023.56310.0003]B1=B4=[001.4811-2.84900-2.84923.7417]B2=B3=[001.19651.0297001.029719.1501]P=[16.2602-0.66401.0734-0.0095-0.66400.3007-0.34550.28711.0734-0.34550.4730-0.3458-0.00950.2871-0.34580.3658]K1=[-0.45090.1745-0.71820.0527-0.28070.0034-0.04710.2639]K2=[-0.46190.1238-0.78120.0351-0.28270.0068-0.08510.1401]K3=[-0.59020.2132-0.89100.0531-0.41070.0117-0.07300.2989]K4=[-0.48910.2109-0.67340.04145-0.31800.0085-0.06310.3893]
which guarantee the stability condition (17). MATLAB Toolbox is used to obtain parameters as k1=0.0251 and k2=0.01539. For comparison with the proposed controller, the PDC feedback gains are designed to be
(37)K^1=[-0.32930.1023-0.25380.02317-0.17830.0021-0.026720.4732]K2=[-0.42780.1845-0.24510.0731-0.21850.0026-0.07520.0923]K3=[-0.60130.3461-0.54820.0231-0.44210.0093-0.06510.3529]K4=[-0.37560.3150-0.53910.0421-0.42500.0023-0.05310.3597].


### 5.1. Sinusoidal Wave Response

Continuous reciprocation is required in order to foster the recovery of extremity function. The sinusoidal wave responses of the proposed approach and PDC for both knee and ankle joints are shown in Figure4. It is evident that angle trajectories of the proposed approach are close to the command. Figure 5 shows that the proposed approach exhibits less tracking errors than does PDC. The peak-peak error and phase lag are listed in Table 2. Because of the interaction of the two joints, PDC has significant angle errors for θ1, which will degrade the rehabilitation effect. However, supervisory control can overcome the coupling effect of the two joints to achieve excellent rehabilitation function for patients.Table 2
Peak-peak error and phase lag for Figure5.
The proposed approach
PDC
Peak-peak error
Phase lag
Peak-peak error
Phase lag
θ
1
θ
2
θ
1
θ
2
θ
1
θ
2
θ
1
θ
2
0.7%
0.35%
4.9°
4.5°
1.3%
0.65%
16.2°
10.8°Sinusoidal wave response for knee and ankle joints.
(a)
The proposed approach
(b)
PDCAngle tracking errors for both knee and ankle.
(a)
The proposed approach
(b)
PDC


### 5.2. Irregular Curve Response

In practical applications, it could be expected that the reference command will change with different input frequencies. The desired trajectories for both knee and ankle joints are(38)θ1=20*0.33(sin⁡(2πf1t)+sin⁡(2πf2t)+sin⁡(2πf3t))θ2=15*0.33(sin⁡(2πf1t)+sin⁡(2πf2t)+sin⁡(2πf3t))
with f1=0.05Hz, f2=0.1Hz, and f3=0.066Hz.Figure6 shows the tracking responses of irregular curves obtained using both the proposed approach and PDC. Tracking errors for the knee and ankle joints are shown in Figure 7. Clearly, the angle error of the proposed approach is average maintained within 2°. However, the proposed approach is capable of adapting to different frequencies.Irregular curve response for both knee and ankle joints.
(a)
The proposed approach
(b)
PDCAngle tracking errors for both the proposed approach and PDC.
(a)
The proposed approach
(b)
PDC


### 5.3. Elliptic Response

The desired end-effect or trajectory is given by(39)xd(t)=0.614-0.015·cos⁡(0.2π·t-π)yd(t)=-0.1·sin(0.2π·t-2π),
where 0≤t≤20 seconds.The end-effect tracking responses in thex, y coordinate for both the proposed approach and PDC are shown in Figure 8, and the end-effect position tracking errors are displayed in Figure 9. It is evident that tracking behavior of the proposed approach is better than that of the PDC. As can be seen, the tracking errors of the proposed approach are within 0.03 m. On the other hand, angle tracking errors of knee and ankle joints are shown in Figure 10.Elliptic response of the proposed approach and PDC.
(a)
The proposed approach
(b)
PDCThe end-effect position tracking errors.
(a)
The proposed approach
(b)
PDCAngle tracking errors of knee and ankle joints.
(a)
The proposed approach
(b)
PDCMoreover, it is difficult to enhance end-effector tracking performance using the PDC algorithm because the PDC cannot overcome the nonlinearity of PAMs and the structural interaction. However, the proposed approach overcomes successfully the coupling effect and parameter uncertainties of the system. As seen in the experimental results, the proposed approach can attain excellent end-effector tracking performance in rehabilitation function.


## 5.1. Sinusoidal Wave Response

Continuous reciprocation is required in order to foster the recovery of extremity function. The sinusoidal wave responses of the proposed approach and PDC for both knee and ankle joints are shown in Figure4. It is evident that angle trajectories of the proposed approach are close to the command. Figure 5 shows that the proposed approach exhibits less tracking errors than does PDC. The peak-peak error and phase lag are listed in Table 2. Because of the interaction of the two joints, PDC has significant angle errors for θ1, which will degrade the rehabilitation effect. However, supervisory control can overcome the coupling effect of the two joints to achieve excellent rehabilitation function for patients.Table 2
Peak-peak error and phase lag for Figure5.
The proposed approach
PDC
Peak-peak error
Phase lag
Peak-peak error
Phase lag
θ
1
θ
2
θ
1
θ
2
θ
1
θ
2
θ
1
θ
2
0.7%
0.35%
4.9°
4.5°
1.3%
0.65%
16.2°
10.8°Sinusoidal wave response for knee and ankle joints.
(a)
The proposed approach
(b)
PDCAngle tracking errors for both knee and ankle.
(a)
The proposed approach
(b)
PDC


## 5.2. Irregular Curve Response

In practical applications, it could be expected that the reference command will change with different input frequencies. The desired trajectories for both knee and ankle joints are(38)θ1=20*0.33(sin⁡(2πf1t)+sin⁡(2πf2t)+sin⁡(2πf3t))θ2=15*0.33(sin⁡(2πf1t)+sin⁡(2πf2t)+sin⁡(2πf3t))
with f1=0.05Hz, f2=0.1Hz, and f3=0.066Hz.Figure6 shows the tracking responses of irregular curves obtained using both the proposed approach and PDC. Tracking errors for the knee and ankle joints are shown in Figure 7. Clearly, the angle error of the proposed approach is average maintained within 2°. However, the proposed approach is capable of adapting to different frequencies.Irregular curve response for both knee and ankle joints.
(a)
The proposed approach
(b)
PDCAngle tracking errors for both the proposed approach and PDC.
(a)
The proposed approach
(b)
PDC


## 5.3. Elliptic Response

The desired end-effect or trajectory is given by(39)xd(t)=0.614-0.015·cos⁡(0.2π·t-π)yd(t)=-0.1·sin(0.2π·t-2π),
where 0≤t≤20 seconds.The end-effect tracking responses in thex, y coordinate for both the proposed approach and PDC are shown in Figure 8, and the end-effect position tracking errors are displayed in Figure 9. It is evident that tracking behavior of the proposed approach is better than that of the PDC. As can be seen, the tracking errors of the proposed approach are within 0.03 m. On the other hand, angle tracking errors of knee and ankle joints are shown in Figure 10.Elliptic response of the proposed approach and PDC.
(a)
The proposed approach
(b)
PDCThe end-effect position tracking errors.
(a)
The proposed approach
(b)
PDCAngle tracking errors of knee and ankle joints.
(a)
The proposed approach
(b)
PDCMoreover, it is difficult to enhance end-effector tracking performance using the PDC algorithm because the PDC cannot overcome the nonlinearity of PAMs and the structural interaction. However, the proposed approach overcomes successfully the coupling effect and parameter uncertainties of the system. As seen in the experimental results, the proposed approach can attain excellent end-effector tracking performance in rehabilitation function.


## 6. Conclusions

In this study, a novel composite fuzzy control is proposed and applied in the two-joint leg rehabilitation device driven by PAMs. The proposed controller is not only capable of decomposing nonlinear systems into a set of linear subsystems, but is also capable of simplifying a complex nonlinear system using linear control techniques, with the control gains determined using MATLAB’s LMI Toolbox based on the Lyapunov stability theorem. Moreover, the supervisory control can overcome the coupling effect for a leg rehabilitation machine. Experimental results show that the system response of the proposed approach was in good agreement with that of the reference input and guarantee robustness to system parameter uncertainties.



---

*Source: 464276-2014-03-18.xml*