1. Introduction
The equations of motion of a control mechanical system contain terms quadratically dependent on velocity that are usually called the Coriolis terms. Controller synthesis becomes tractable in the absence of these nonlinear terms [6], so it is useful to find a transformation that eliminates them from the equations of motion.
A control mechanical system is called quasilinearizable if there is a linear transformation of the velocity variables such that the Coriolis terms all vanish after the transformation. There has been active research on quasilinearization [2, 5, 7, 8], but the results were obtained by the zero curvature condition or by some complicated PDE conditions, producing restrictive outcomes. Then, very strong results were finally obtained in [4] where easily verifiable quasilinearizability conditions were derived.
In this paper we consider feedback transformations as well as state transformations, in order to increase possibility of removing the Coriolis terms from the dynamics. A control mechanical system is called feedback quasilinearizable if all Coriolis terms can be eliminated by a linear velocity transformation followed by a feedback transformation. We here obtain a sufficient condition for feedback quasilinearizability and apply it to prove the feedback quasilinearizability of the Acrobot system. We also derive a condition for partial quasilinearizability via a linear velocity transformation in the course of obtaining the result on feedback quasilinearizability.
2. Main Results
2.1 Review of quasilinearization theory
We review the theory of quasilinearization of mechanical systems in [4] from a slightly different viewpoint. We here use a linear bundle map from the tangent bundle TQ of a given configuration space Q to its cotangent bundle T*Q instead of a linear bundle map from TQ to itself. This different style of presentation, however, does not affect the validity of the results in [4].
Let Q be an n - dimensional manifold and q = (qi) a local coordinate system on Q ; refer to [1], [3] for manifolds theory. Let TQ and T*Q denote the tangent bundle and the cotangent bundle of Q , respectively. The natural pairing between TQ and T*Q is denoted by <, > . The natural local coordinate bases of TQ and T*Q are used.
TQ = span{∂1,⋯,∂n} T*Q = span{dq1,⋯,dqn} .
The symbol ∂i is also used as the operator of partial differentiation with respect to qi. We use the Einstein summation convention throughout this paper and the following convention for the ranges of various indices:
i, j, k, l, r, s = 1,⋯,n; a, b, c = 1,⋯,p.
Consider a control mechanical system on the configuration space Q with Lagrangian
and p-dimensional control bundle W ⊂ T*Q , where m = (mij) is the positive definite symmetric mass matrix and V(q) is the potential energy of the system. Since our results are all local, we assume that W is generated by p independent 1-forms as follows:
W = span{W1,⋯,Wp}
where each 1-form Wa, a = 1,⋯, p, is written in coordinates as
Wa = Wiadqi.
The equations of motion of this control mechanical system are given by
for i = 1,⋯, n, where u = (ua) ∈ ℜp is the control vector. Here, mij denotes the (i, j) entry of the inverse matrix of m = (mij) and Γijk are the Christoffel symbols defined by
The quadratic terms in the equations of motion are called Coriolis terms.
Consider an invertible linear bundle map A : TQ → T*Q given by
In coordinates,
where α = αidqi. Let (Bij) be the inverse matrix of (Aij ), i.e., , where is the Kronecker delta. In (x, α) coordinates on T*Q , the equations of motion (1) become
where i = 1, ⋯, n, or in vector form
where
Notice that all the Coriolis terms vanish in the equations in (4) if and only if
for all i, j, k, in which case the equations of motion become
or in vector form
Definition 2.1: A control mechanical system is said to be quasilinearizable if there is an invertible linear transformation of the form (2) that transforms the equations of motion of the system (1) to the form (6) and (7).
We can regard the configuration space Q of a mechanical system as a Riemannian manifold equipped with the metric m = (mij) that is induced from the kinetic energy of the system. A vector field X = Xi∂i on a Riemannian manifold (Q, m) is called a Killing vector field if it satisfies
for all 1 ≤ i ≤ j ≤ n. Letting α = mX = m jk Xk dqj, we can write (8) as
in terms of the 1-form α. A 1-form that satisfies (9) is called a Killing 1-form. Both (8) and (9) are called the Killing equation. Comparison of (5) and (9) implies that Eq. (5) is the Killing equations in (9) for the 1-form Ai := Aijdqj for each i. Hence, a quasilinearizing transformation consists of n pointwise independent Killing 1-forms, where each row of A is a Killing 1-form.
Let iso (Q, m) denote the set of all Killing vector fields on (Q, m). It is a Lie algebra over ℜ under the usual bracket operation on vector fields. Let Δ denote the distribution on Q that is generated by Killing vector fields, i.e.
for each q ∈ Q. The rank of Δq is, by definition, the dimension of Δq as a vector subspace of TqQ. Then the quasilinearizability can be geometrically stated as follows.
Theorem 2.2 ([4]): Let q be a point in (Q, m).The quasilinearization of the system (1) is possible around q if and only if Δq = TqQ , i.e., rank Δq = dim Q.
2.2 Partial quasilinearization and feedback quasilinearization
We now pose the following two main questions for control mechanical systems that are not quasilinearizable:
Q1. (Partial Quasilinearization) How many of the equations in (4) can be made free of the Coriolis terms via a transformation of the form (2)?
Q2. (Feedback Quasilinearization) If an affine feedback transformation of the form
with h : Q → ℜp and u ∈ ℜp , is allowed in addition to the linear transformation of the form (2), when can a given system be transformed to the form (6) and (7), i.e, to the following form
which is free of the Coriolis terms?
Definition 2.3: A control mechanical system is called feedback quasilinearizable if its equations of motion can be transformed to the form (12) and (13) via a transformation of the form (2) followed by a feedback transformation of the form (11).
Definition 2.4: A point q in (Q, m) is called regular if the rank of the distribution Δ defined in (10) is constant in a neighborhood of q .
We now provide an answer to the first question we posed in the beginning of this section.
Theorem 2.5 (Partial Quasilinearization): Let q0 be a regular point in (Q, m). Then, at least k -equations can be made free of the Coriolis terms via an invertible transformation of the form (2) around q0 if and only if rank Δ ≥ k in a neighborhood of q0 .
Proof: (⇒) By hypothesis there is a linear transformation such that the first k - equations can be written as
for i = 1,⋯, k in a neighborhood of q0 . In other words Eq. (5) holds for i = 1,⋯, k. Hence, the first k row vectors of A are pointwise independent Killing 1-forms, which implies that rank Δ ≥ k in a neighborhood of q0 .
(⇐) This direction can be proven similarly.
The above theorem can be also interpreted as follows: k is the maximum number of the equations that can be made free of Coriolis terms via a transformation of the form (2) around a regular point q0 if and only if rank Δ = k in a neighborhood of q0 .
We now answer the second question posed in the beginning of this section.
Theorem 2.6 (Feedback Quasilinearization): A control mechanical system is feedback-quasilinearizable around a regular point q0 if
for each q in a neighborhood of q0 , where Δ0 is the codistribution on Q that annihilates Δ , i.e., pointwise
Proof: Let k be the constant rank of Δ around q0 . Then there exist k Killing vector fields X1, ⋯, Xk that span Δ pointwise around q0 . Choose (n - k) more vector fields X k+1, ⋯, Xn such that the set of vector fields {X1, ⋯, Xn} span TQ around q0 . One can find (n - k) 1-forms βk+1,⋯, βn in Δ0 around q0 such that
for k + 1 ≤ i ≤ n. Since Δ0 ⊂ W by hypothesis, there exist vectors uk + 1,⋯, un in ℜp such that
βi = Wui
for k + 1 ≤ i ≤ n, where βi and ui are assumed to be in column vector form. Let
the first k rows of which are Killing 1-forms since X1, ⋯, Xk are Killing vector fields. Let
or in coordinates . Change coordinates from to α to transform (1) to (3) and (4), where the first k - equations in (4) become free of the Coriolis terms. Apply the following control u ∈ ℜp
where
for k + 1 ≤ i ≤ n . It is then easy to see that the system (3) and (4) is transformed via this feedback control to the system (12) and (13). Therefore, the system is feedback quasilinearizable around q0 .
3. Example
Consider the Acrobot system in Fig. 1, where there is an actuation u on the outer joint. Let M1 and M2 be the masses of the bobs and l1 and l2 the lengths of the massless rods. The gravitational acceleration is denoted by g. Let θ1 denote the angle of the first rod measured counter-clockwise from the upward vertical, and θ2 the angle measured counterclockwise from the ray containing the first rod to the second rod.
Fig. 1.The Acrobat system
The Lagrangian of the system is given by
where
The scalar curvature RS of the metric m = (mij) is computed as
which is not constant. Hence, the system is not quasilinearizable by Theorem III.1 in [4].
Let us now investigate feedback quasilinearizability of this system. The Acrobat has only one Killing vector field up to a scalar factor and it is given by X = ∂1 , which can be easily obtained using software Maple. Hence,
Δ = span{∂1}, Δ0 = span{dθ2}.
The control bundle of the Acrobot is given by
W = span{dθ2}.
Since Δ0 ⊂ W, the Acrobot is feedback quasilinearizable by Theorem 2.6.
References
- R. Abraham, J.E. Marsden and T. Ratiu, Manifolds, Tensor Analysis, and Applications, 2nd Edition, Springer, 1998.
- N.S. Bedrossian, “Linearizing coordinate transformation and Riemannian curvature,” in Proc. IEEE Conf. Dec. Control, Tucson, AZ, Dec., pp. 80-85, 1992.
- W.M. Boothby, An Introduction to Differentiable Manifolds and Riemannian Geometry, 2nd Ed., Academic Press, 2002.
- D.E. Chang and R.G. McLenaghan, “Geometric criteria for the quasilinearization of the equations of motion of mechanical systems,” IEEE Trans. Automatic Control, 58 (4), 1046-1050, 2013. https://doi.org/10.1109/TAC.2012.2218671
- D. Koditschek, “Robot kinematics and coordinate transformations,” in Proc. IEEE Conf. on Decision and Control, Ft. Lauderdale, FL, Dec. , pp. 1-4,1985.
- A. Loria and E. Panteley, “A separation principle for Euler-Lagrange systems”, in New Directions in Nonlinear Observer Design (H. Nijemeijer and T.I. Fossen eds.), London, Springer Verlag, 1999.
- M.W. Spong, “Remark on robot dynamics: Canonical transformations and Riemannian geometry,” in Proc. IEEE Int. Conf. Robot. Autom., Nice, France, May, pp. 554-559, 1992.
- A. Venkatraman, R. Ortega, I. Sarras, and A. van der Schaft, “Speed observation and position feedback stabilization of partially linearizable mechanical systems,” IEEE Trans. Automatic Control, 55 (5), pp. 1059-1074, 2010. https://doi.org/10.1109/TAC.2010.2042010
Cited by
- A globally exponentially stable speed observer for a class of mechanical systems: experimental and simulation comparison with high-gain and sliding mode designs 2017, https://doi.org/10.1080/00207179.2017.1404130