Published: 30 June 2015

Composite state variable based nonlinear backstepping design for the underactuated TORA system

Weiping Guo1
Diantong Liu2
1, 2School of Computer and Control Engineering, Yantai University, Yantai City, People’s Republic of China
Corresponding Author:
Diantong Liu
Views 67
Reads 39
Downloads 1077

Abstract

A nonlinear vibration controller is proposed for the translational oscillators with rotating actuator (TORA) system with the recursive technology. A composite state variable (CSV) is defined for the TORA system to start the recursive process. The design procedure treats the some state variables as virtual control inputs to design the virtual controllers step by step until the nonlinear vibration controller is obtained. The system stability is studied via a stability theorem and simulation results show the validity of the proposed controller.

1. Introduction

The TORA is composed of a non-actuated translational oscillators and an actuated pendulum so that it has been seen as a classical underactuated mechanical system. It was firstly developed as a simplified model for a dual spin spacecraft [1]. The TORA system has been a research focus with two reasons. One research motivation is the possible use to suppress the translational vibration with a rotational actuation and another is an independent interest to control the TORA as a benchmark problem to design or validate a nonlinear control system or an underactuated control system.

The control problems on TORA have a considerable amount of attention from the all world researchers. There are two focuses on the methods to design a controller, i.e. intelligent design methods and nonlinear design methods. Among the intelligent design methods, some fuzzy based control methods are proposed. The self-tuning fuzzy sliding mode control methods are designed [2, 3] but there are some fuzzy rules and parameters to be preset. A decoupled method is designed to divide the system into two subsystems [4] and two fuzzy controllers corresponding to the subsystems are designed to stabilize the TORA system. However, the fuzzy controllers must be predesigned, and stability cannot be guaranteed. Moreover, some neural network based controllers are proposed. An adaptive neural network controller with SPSA Algorithm is designed [5]. In the literature [6], a method was proposed that is based on system approximation using n neural network. However, a time-consuming learning process is needed and the disturbance rejection cannot be considered. Among the nonlinear design methods, the energy-based design [7-9], the recursive design [10, 11] and sliding mode based design [12], the design processes and the obtained controllers are relatively cumbersome. A comparison study of nonlinear control techniques for the TORA is performed [7]. The energy-based hybrid, entropy-based hybrid, and sliding mode control techniques are experimentally validated and compared on the example of the RTAC testbed. The experimental results indicate that the energy- and entropy-based hybrid controllers use less control effort than the sliding mode controller while stabilizing the zero equilibrium state in less time.

The backstepping design provides a systematic framework and recursive design methodology for nonlinear control [13]. The design procedure treats the state variables as virtual control inputs to design the virtual controllers step by step, and illustrates the stability via the Lyapunov stability theorem. However, a nonlinear feedback cascade model is needed in the traditional backstepping design process and the obtained controllers are relatively cumbersome. In this paper, we propose a backsteppping process that is not based on a nonlinear feedback cascade model but begins from a predefined CSV, and a simple control is expected.

This paper is organized as follows. The TORA system is transformed into a nonlinear system through a collocated partial feedback linearization. Then a nonlinear recursive design method is adopted and a CSV is defined to start the recursive process. The design procedure treats the some state variables as virtual control inputs to design the virtual controllers step by step, and the stability is studied via a stability theorem.

2. The TORA system model

The TORA is shown schematically in Fig. 1. The oscillator consists of a cart of mass m1 connected to a fixed wall by a linear spring of stiffness k. The cart is constrained to have one-dimensional travel and its position is q1. The rotational actuator with the output torque is attached to the cart and it has a pendulum with the equivalent mass m2, the rotate radius r, the moment of inertia I and the rotate angle q2. Since the rotating torque can not drive the cart directly, it is a benchmark example of the underactuated mechanical system, which has one control input τ and two configuration variables (q1,q2), and its Euler-Lagrange equations of motion can be obtained as:

1
(m1+m2)q¨1+m2rcos(q2)q¨2-m2rsinq2q˙22+kq1=0,m2rcosq2q¨1+m2r2+Iq¨2+m2grsinq2=τ.

Fig. 1TORA system configuration

TORA system configuration

Through the following collocated partial feedback linearization [14]:

2
τ=(m2r2+I)-(m2rcos(q2))2m1+m2u+(m2r)2sin(q2)cos(q2)m1+m2q˙12
-km2rq1cosq2m1+m2+m2grsinq2.

The dynamics can be reduced to:

3
q¨1=m2rsinq2q˙22-kq1-m2rcosq2um1+m2,q¨2=u.

3. The nonlinear controller design through backstepping

The controller is designed with the recursive technology and the design process is:

Step 1.

A composite state variable for the nonlinear underactuated TORA system is defined as:

4
x1=k1q1+k2q2+k3m1+m2q˙1+m2rcosq2q˙2,

where, k1, k2 and k3 are some positive constants for the system stability to be proved. The derivative of the CSV can be obtained:

5
x˙1=k1q˙1+k2q˙2+k3m1+m2q¨1+m2rcosq2q¨2-m2rsinq2q˙22.

According to Eq. (1), Eq. (5) yields:

6
x˙1=k1q˙1+k2q˙2-k3kq1.

The k2q˙2 is looked as the virtual input to control the above equation and define a reference trajectory q˙2r for q˙2 to follow as:

k2q˙2r=-k1q˙1+k3kq1-b1x1,

where b1 is a design constant, it leads to an error e1 defined as:

e1=k2q˙2-k2q˙2r=k1q˙1+k2q˙2-k3kq1+b1x1=x˙1+b1x1,

i.e.:

7
x˙1=e1-b1x1.

Its derivative is:

8
e˙1=k1q¨1+k2q¨2-k3kq˙1-b1x˙1.

According to Eq. (2), Eq. (8) yields:

9
e˙1=k1m2rsinq2q˙22-kq1-m2rcosq2um1+m2+k2u-k3kq˙1-b1x˙1.

Consider a scalar positive definite Lyapunov function given by:

10
V=12x12+12e12.

Using Eq. (7) and (9), the time derivative V˙ is given by:

11
V˙=x1x˙1+e1e˙1=-b1x12+(x1+e˙1)e1
=-b1x12+k1m2rsinq2q˙22-kq1-m2rcosq2um1+m2+k2u-k3kq˙1-b1x˙1+e˙1e1.

Step 2.

In step 1, the time derivative of the Lyapunov function V is obtained in Eq. (11). We note that the variable u enters the right hand side of Eq. (11).

In order to make the V˙ be negative definite, we can make the following equation hold:

12
k1m2rsinq2q˙22-kq1-m2rcosq2um1+m2+k2u-k3kq˙1-b1x˙1+x1=-b2e1.

Such that:

13
V˙=-b1x12-b1e12.

Therefore, the control law can be obtained from Eq. (12) as:

14
u=k2-k1m2rcosq2m1+m2-1k3kq˙1-b1x˙1-b2x2+k1m2rsinq2q˙22-kq1m1+m2-x1.

Since the system parameters m1, m2 and r are known and cosq2[-1, 1], the choices of k1 and k2 can make Eq. (14) away from the singularity.

Theorem 1. The TORA system described by Eq. (1) is asymptotically stable to its equilibrium point (q1,q˙1,q2,q˙2)=(0,0,0,0)under the control input Eq. (2) with Eq. (14) and positive constants k1, k2, k3, b1, b2.

Proof:

The backstepping design process has proved: the time derivative V˙ of the chosen positive definite Lyapunov function V is negative definite.

Integrating both sides of Eq. (13) yields:

00V˙dτ=00(-b1x12-b1e12)dτ.

Then:

Vt-V0=00-b1x12-b1e12dτ.

We find that:

Vt=12x12+12e12=V0+00-b1x12-b1e12dτV0<.

Therefore it can be obtained that x1L and e1L, i.e. supx1=x1<,supe1=e1<.

At the same time, from Eq. (13):

V˙=x1x˙1+e1e˙1=-b1x12-b1e12<.

It is obvious that x˙1L and e˙1L, i.e. supx˙1=x˙1<, supe˙1=e˙1<.

Now let:

s1=q1+k3(m1+m2)k1q˙1, s2=q2+k3m2rk2cosq2q˙2.

Thus:

x1=k1s1+k2s2.

From our research result (Theorem 1 in [15]), it can be known that s1L2, s2L2 and s1L, s˙1L, s2L, s˙2L. According to Barbalat’s lemma, limts1=0, limts2=0, i.e. when t:

15
q1+k3(m1+m2)k1q˙1=0,
16
q2+k3m2rk2cosq2q˙2=0.

For Eq. (15), a positive definite Lyapunov function V1=q12/2 is selected and V˙1=q1q˙1=-k3(m1+m2)q˙12/k1 is negative definite. For Eq. (16), a positive definite Lyapunov function V2=q22/2 is selected and V˙2=q2q˙2=-k3m2rcosq2q˙22/k2is negative definite when q2(-π/2,π/2). Therefore, the TORA system is asymptotically stable to the equilibrium point (q1,q˙1,q2,q˙2)=(0, 0, 0, 0).

Remark 1. As we know, a nonlinear feedback cascade model is needed in the traditional backstepping design process. The proposed design do not need a nonlinear feedback cascade model and it begins from a predefined CSV and treats some state variables as virtual control inputs to design the nonlinear controllers step by step. Therefore, it may be a choice to realize a backstepping design for those systems that can not be transformed into a nonlinear feedback cascade model.

Remark 2. The TORA is a system that can be transformed into a nonlinear feedback cascade model and the backstepping design can be realized [16]. However, the design process and the obtained controller are more cumbersome because the design process is composed of a second-order system design and a first-order system design.

Remark 3. From the process to prove the Theorem 1, it can be seen that the Eq. (15) is convergent only when k1 and k3 have same sign and the Eq. (16) is convergent only when k2 and k3 have same sign, i.e. only when k1, k2 and k3 have same sign, the system is stable. Therefore, some positive constants are chosen for k1, k2 and k3 in Eq. (4).

4. Simulation studies

In order to test the proposed nonlinear control algorithm, the following system parameters [13] are used: m1=10, m2=1, k=5, r=1, I=1. The parameters of the nonlinear controller are chosen as k1=1, k2=3, k3=3, b1=5, b2=5. The simulation results are shown in Fig. 2 and Fig. 3.

The simulation results in Fig. 2 are obtained under the initial state q1,q˙1,q2,q˙2=(0.1, 0, 0, 0) and the simulation results in Fig. 3 are obtained under the initial state q1,q˙1,q2,q˙2=(0, 0, 180°, 0). It can be seen form the simulation results and many other simulations that the TORA system is asymptotically stable under any initial states with the proposed control algorithm. At the same time, it is shown that it is possible to suppress the translational vibration with a rotational actuation. On the other hand, the control performance can be improved through adjusting the parameters of the proposed controller. Lots of simulation experiments show that the parameters k1, k2 respectively correspond to the system states q1, q2, therefore it is easy to adjust the parameters for a improved system performance. Moreover, compared with the traditional backstepping design [16], it can be seen that the control performance of the proposed method is improved from the simulation results in [16] with the same system parameters or more difficult-to-control initial state. If only the simulation or experimental results are compared with the other control algorithms [2-12], a shorter stabilization time, a better states response and a lower number of oscillations are realized. Therefore, the proposed controller is simpler and more efficient.

Fig. 2Simulation results under system initial state q1,q˙1,q2,q˙2=(0.1, 0, 0, 0)

Simulation results under system initial state q1,q˙1,q2,q˙2=(0.1, 0, 0, 0)

a) Position and velocity of the cart

Simulation results under system initial state q1,q˙1,q2,q˙2=(0.1, 0, 0, 0)

b) Angle and angle velocity of the pendulum

Simulation results under system initial state q1,q˙1,q2,q˙2=(0.1, 0, 0, 0)

c) Control torque

Simulation results under system initial state q1,q˙1,q2,q˙2=(0.1, 0, 0, 0)

d) Responses of CSV x1 and error e1

Fig. 3Simulation results under system initial state q1,q˙1,q2,q˙2=(0, 0, 180°, 0)

Simulation results under system initial state q1,q˙1,q2,q˙2=(0, 0, 180°, 0)

a) Position and velocity of the cart

Simulation results under system initial state q1,q˙1,q2,q˙2=(0, 0, 180°, 0)

b) Angle and angle velocity of the pendulum

Simulation results under system initial state q1,q˙1,q2,q˙2=(0, 0, 180°, 0)

c) Control torque

Simulation results under system initial state q1,q˙1,q2,q˙2=(0, 0, 180°, 0)

d) Responses of CSV x1 and error e1

5. Conclusions

Using a collocated partial feedback linearization, the underactuated TORA is transformed to a nonlinear system that is used in the recursive process. A CSV is defined and a nonlinear control algorithm is proposed with the recursive technology. The design procedure treats the state variables as virtual control inputs to design the virtual controllers step by step, and the system stability is studied via a stability theorem.

References

  • Robert T. B., Bernstein D. S., Coppola V. T. A benchmark problem for nonlinear control design. International Journal of Robust and Nonlinear Control, Vol. 8, Issue 4-5, 1998, p. 307-310.
  • Hung L. C., Lin H. P., Chung H. Y. Design of self-tuning fuzzy sliding mode control for TORA system. Expert Systems with Applications, Vol. 32, Issue 1, 2007, p. 201-212.
  • Li C. D., Yi J. Q., Zhao D. B. Control of the TORA system using SIRMs based type-2 fuzzy logic. Proceedings of IEEE International Conference on Fuzzy System, Jeju Island, 2009, p. 694-699.
  • Lee C. H., Pan H. Y., Chang H. H., Wang B. H. Decoupled adaptive type-2 fuzzy controller (DAT2FC) design for nonlinear TORA systems. IEEE International Conference on Fuzzy Systems, Vancouver, Canada, 2006, p. 506-512.
  • Lee C. H., Yu T. M., Chien J. C. Adaptive neural network controller design for a class of nonlinear systems using SPSA algorithm. International MultiConference of Engineers and Computer Scientists, Hong Kong, 2011, p. 868-873.
  • Lee C. H., Chang S. K. Experimental implementation of nonlinear TORA system and adaptive backstepping controller design. Neural Computing and Applications, Vol. 21, Issue 4, 2012, p. 785-800.
  • Avis J. M., Nersesov S. G., Nathanb R., et al. A comparison study of nonlinear control techniques for the RTAC system. Nonlinear Analysis: Real World Applications, Vol. 11, Issue 4, 2010, p. 2647-2658.
  • Jankovic M., Fontaine D., Kokotovic P. V. TORA example: cascade and passivity based control designs. IEEE Transactoins on Control System Technology, Vol. 4, Issue 3, 1996, p. 292-297.
  • Gao B. T., Zhang X. H., Chen H. J., et al. Energy-based control design of an underactuated 2-dimensional TORA system. Proceedings of the 2009 IEEE/RSJ International Conference on Intelligent Robots and Systems, St. Louis, 2009, p. 1296-1301.
  • Bupp R., Bernstein D. S., Coppola V. T. Experimental implementation of integrator backstepping and passive nonlinear controllers on the RTAC testbed. International Journal of Robust and Nonlinear Control, Vol. 8, Issue 4-5, 1998, p. 435-457.
  • Wang Y. Q., Li S., Chen Q. W. A recursive design method of controller for the stabilization of the translational oscillator with a rotational actuator system. Proceedings of the IEEE International Conference on Computer Science and Automation Engineering, Shanghai, 2011, p. 558-561.
  • Chen Y. F., Huang A. C. Controller design for a class of underactuated mechanical systems. IET Control Theory Application, Vol. 6, Issue 1, 2012, p. 103-110.
  • Khalil H. K. Nonlinear Systems, 3rd Edition. Prentice Hall, NJ, 2000.
  • Spong M. W. Underactuated Mechanical Systems. Control Problems in Robotics and Automation. Springer Verlag, London, UK, 1997.
  • Wang W., Yi J., Zhao D., Liu D. Design of a stable sliding-mode controller for a class of second-order underactuated systems. IEEE Proceedings on Control Theory and Applications, Vol. 151, Issue 6, 2004, p. 683-690.
  • Liu D. T., Guo W. P. Nonlinear backstepping design for the underactuated TORA system. Journal of Vibroengineering, Vol. 16, Issue 2, 2014, p. 552-559.

About this article

Received
01 November 2014
Accepted
13 January 2015
Published
30 June 2015
SUBJECTS
Vibration generation and control
Keywords
composite state variable
TORA
backstepping design
vibration control
Acknowledgements

This work is supported by the National Natural Science Foundation of China (61175086) and a Project of Shandong Province Higher Educational Science and Technology Program (No. J13LN25).