Next Article in Journal
Anti-Yeast Synergistic Effects and Mode of Action of Australian Native Plant Essential Oils
Next Article in Special Issue
Modified Infinite-Time State-Dependent Riccati Equation Method for Nonlinear Affine Systems: Quadrotor Control
Previous Article in Journal
Design and Analysis of the Rolling and Jumping Compound Motion Robot
Previous Article in Special Issue
Development of Task-Space Nonholonomic Motion Planning Algorithm Based on Lie-Algebraic Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Mechanical Feedback Classification of Linear Mechanical Control Systems

by
Marcin Nowicki
1,* and
Witold Respondek
2,*
1
Institute of Automatic Control and Robotics, Poznan University of Technology, Piotrowo 3a, 61-138 Poznań, Poland
2
Laboratoire de Mathématiques, INSA de Rouen Normandie, 685 Av. de l’Université, 76801 Saint-Etienne-du-Rouvray, France
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2021, 11(22), 10669; https://doi.org/10.3390/app112210669
Submission received: 5 October 2021 / Revised: 29 October 2021 / Accepted: 8 November 2021 / Published: 12 November 2021

Abstract

:
We give a classification of linear nondissipative mechanical control system under mechanical change of coordinates and feedback. First, we consider a controllable case that is somehow a mechanical counterpart of Brunovský classification, then we extend the result to all linear nondissipative mechanical systems (not necessarily controllable) which leads to a mechanical canonical decomposition. The classification of Lagrangian systems is given afterwards. Next, we show an application of the classification results to the stability and stabilization problem and illustrate them with several examples. All presented results in this paper are expressed in terms of objects on the configuration space R n only, while the state-space of a mechanical control system is R n × R n consisting of configurations and velocities.

1. Introduction

In this paper we consider the problem of classification of linear mechanical control systems under mechanical feedback transformations. Therefore, by providing a solution to the classification problem, we answer several questions:
(a)
Given two linear mechanical systems, how to determine whether they are equivalent?
(b)
Is there a set of complete invariants that are computable in terms of objects on the configuration space only?
(c)
Is there a distinguished normal (or canonical) form?
We consider the above-defined classification problem in three important cases, namely for controllable and uncontrollable mechanical systems, and for the subclass of Lagrangian systems. A classification of controllable linear (first-order) systems x ˙ = A x + B u under general linear transformations and general linear feedback has been solved in the celebrated Brunovský classification [1], see also [2]. A general classification (including not necessarily controllable systems) leads to the canonical decomposition [3,4] and a short yet concise note [5]. Here, we consider the novel problem of classification of linear mechanical systems under linear transformations that respect splitting into configurations and velocities, and linear mechanical feedback. It turns out that the special form of considered (second-order) systems and the special form of (mechanical) transformations yield the counterpart of the above-mentioned classical results. In other words, we deal with a smaller class of control systems (than general linear systems) and we use more subtle mechanical feedback transformations (than general feedback transformations) and yet the invariants are perfectly analogous to those of the general case. Moreover, they can be computed on a half of the state space, namely using objects defined on the configurations space only. What is more, we show that mechanical feedback transformations are perfectly adapted to the class of mechanical systems, namely the classification of mechanical systems is the same if, instead of mechanical feedback transformations, we use all linear feedback transformations. Our analysis of the classification problem implies a series of results for stability and stabilization of linear mechanical systems that we present in Section 6. The problem of stability of motion of linear mechanical systems has been studied extensively, see e.g., [6,7], Section 5.2 in [8] and the reference therein. In those papers, certain structural assumptions (symmetry, positive definiteness, etc.) on matrices of the system are assumed, and what is crucial, dissipative forces are allowed an thus asymptotic stability is concerned. We present new results in this field. We do not assume, a priori, any structure on the mechanical system (see Section 2). Moreover, we consider mechanical systems without dissipative terms. This leads to important conclusions that, first, such systems are never asymptotically stable nor can be asymptotically stabilized by feedback. Second, only the Lagrangian subclass of mechanical control systems can be stable, therefore the problem of stabilization of (general) mechanical system reduces, actually, to the problem of finding mechanical transformations that make the system Lagrangian, which fits well to our classification problem.
The paper is organized as follows. In Section 2, we define the class of linear mechanical control systems and its special subclass, namely Lagrangian control systems, and we introduce linear mechanical feedback transformations. In Section 3 and Section 4, the complete classification under mechanical feedback transformations of, respectively, controllable (Theorem 1) and uncontrollable (Theorem 2) linear mechanical control systems is presented. Then, the classification of linear Lagrangian control systems is given in Section 5. In Section 6, we formulate results about stability and stabilization and show in Theorem 3 that, within the class of linear mechanical non-dissipative control systems, only Lagrangian ones can be stabilized (not asymptotically). Finally, we illustrate our theory by two multibody mass-spring systems. We conclude the paper in Section 8, where we summarize our results.

2. Problem Statement

2.1. Linear Mechanical Control Systems

In this subsection, we introduce the object of our study, namely linear mechanical control systems. This class is larger than linear Lagrangian control systems that will also be introduced in this subsection and form an important subclass in our study.
Consider the linear mechanical control system with n degrees of freedom and m controls
x ¨ = E x + B u ,
where x R n are the configurations (the generalized coordinates). The matrix E is an n × n constant real matrix corresponding to an uncontrolled (depending on configurations only but, possibly, non potential) force in the system and the input matrix B is an n × m constant real matrix describing external forces controlled by the controls u R m .
Equivalently, system (1) can be represented as a first-order system on the state-space R n × R n , equipped with coordinates ( x , y ) denoting, respectively, configurations and velocities
LMS x ˙ = y y ˙ = E x + B u ,
or as a linear control system of dimension 2 n , with coordinates z = x , y R n × R n (to be precise z is a “stacked” vector z = ( x T , y T ) T , however we decide to skip this notation for a clarity sake), given by, compare [9]
z ˙ = A ^ z + B ^ u ,
where
A ^ = 0 I n E 0 , B ^ = 0 B .
Notice that the equation y ˙ = E x + B u is not of the more general form y ˙ = E x + L y + B u , that is, contains neither gyroscopic nor dissipative terms L y , see e.g., [6,7,9,10] for a discussion of both classes.
Our obvious inspiration are linear Lagrangian (conservative) systems which constitute a subclass of (2) whose configuration space is the real vector space R n , equipped with an inner product on the space R n of velocities given by a real valued quadratic form 1 2 x ˙ T M x ˙ (describing the kinetic energy of the system), where M is a constant real symmetric and positive definite matrix ( M T = M > 0 ) . Moreover, we consider the potential energy given by a quadratic form V = 1 2 x T P x , where P is a symmetric potential matrix ( P T = P ) .
The corresponding quadratic Lagrangian reads L = 1 2 x ˙ T M x ˙ 1 2 x T P x . The derivation of the Euler-Lagrange equations yields the second-order system:
M x ¨ + P x = K u ,
where K is an n × m real matrix whose columns are vectors corresponding to the external controlled forces. Straightforward calculations show that any Lagrangian linear mechanical control system can be represented by a particular form of (2), namely
LLS x ˙ = y y ˙ = E L x + B u , where E L = M 1 P and B = M 1 K .
Notice that for Lagrangian systems (6) the matrix B can be any (since K is arbitrary) but the matrix E L is special, namely the product of two symmetric matrices (the first being invertible), which we call the Lagrangian structure.
Another characterization of (6) can be formulated as follows [11].
Proposition 1.
The linear mechanical control system (2) is a Lagrangian system of the form (6) if and only if there exists a real invertible map T : R n R n that diagonalizes the matrix E, that is, T E T 1 = diag ( λ 1 , , λ n ) .
For a proof see [11]. The above proposition can be rephrased as follows. All distinct eigenvalues λ j of the matrix E, 1 j q , where q is the number of distinct eigenvalues, are real, i.e., λ j R , and there are no Jordan blocks. More precisely, the algebraic multiplicity μ j of λ j (indicating how many times λ j appears as a root of the characteristic polynomial of E) is equal to its geometric multiplicity γ j = dim ker ( λ i I n E ) , i.e., the dimension of the eigenspace associated with λ j .
Remark 1.
Obviously, Lagrangian control systems can be represented using the Hamiltonian formulation, as Hamiltonian control systems. There exists a widely studied branch of control theory for this class of systems. For a relation between Lagrangian systems and Hamiltonian systems, including nonlinear case, see [12]. For a survey of port-Hamiltonian systems see [13].

2.2. Linear Mechanical Feedback Transformations

In this subsection, we introduce mechanical feedback transformations under which we classify linear mechanical systems, namely mechanical changes of coordinates and mechanical feedback.
A linear mechanical transformation is given by a linear transformation of the following form z ˜ = T z , where T = T 0 0 T , z = ( x , y ) and z ˜ = ( x ˜ , y ˜ ) , i.e.,
x ˜ = T x , y ˜ = T y ,
where T is an invertible n × n real matrix. This linear transformation preserves configurations, i.e., maps x-coordinates into x ˜ -coordinates. Moreover, since the derivatives of configurations are velocities, it induces the linear transformation y ˜ = T y (given by the same T) on velocities that maps the equation x ˙ = y into x ˜ ˙ = y ˜ .
The linear mechanical feedback is
u = F x + G u ˜ ,
where F is an m × n matrix and G is an m × m invertible matrix. The linear mechanical system (2) transformed by the transformations (7) and (8) reads
x ˜ ˙ = T y = y ˜ y ˜ ˙ = T E + B F T 1 x ˜ + T B G u ˜ = E ˜ x ˜ + B ˜ u ˜ .
Definition 1.
Two systems (2), given by x ˙ = y , y ˙ = E x + B u , and ( LMS ˜ ) , given by x ˜ ˙ = y ˜ , y ˜ ˙ = E ˜ x ˜ + B ˜ u ˜ , are called linear mechanical feedback equivalent, briefly LMF-equivalent, if there exist mechanical transformations (7) and (8) such that E ˜ = T ( E + B F ) T 1 and B ˜ = T B G .
It is natural to consider mechanical transformations (7) and (8), since the class of linear mechanical control systems (2) is closed under those transformations, i.e., the transformed system (9) is linear and mechanical. However, in case of the subclass of Lagrangian systems (6), the mechanical transformations need not preserve the Lagrangian structure (see Section 5 for a detailed analysis).
The group of linear mechanical transformations L M F , consisting of triplets T , F , G , preserves trajectories, that is, any element T , F , G of that group maps the trajectories of (2) into those of its LMF-equivalent system ( LMS ˜ ) , given by (9). Indeed, if z t , z 0 , u ( t ) is a trajectory of (2), passing through z 0   =   ( x 0 , y 0 ) and corresponding to a control u ( t ) , then z ˜ t , z ˜ 0 , u ˜ ( t ) = T z t , z 0 , u ( t ) is a trajectory of (9), passing through z ˜ 0 = T z 0 = ( T x 0 , T y 0 ) and corresponding to u ˜ ( t ) , where u ( t ) = F x ( t ) + G u ˜ ( t ) . Moreover, via T : R n R n , it establishes a correspondence between trajectories in the configuration space R n , i.e., x ˜ t , z ˜ 0 , u ˜ ( t )   =   T x ( t , z 0 , u ( t ) ) , making the following diagram commutative (notice, however, that π z ( t , z 0 , u )   =   x ( t , z 0 , u ) depends on z 0 = ( x 0 , y 0 ) consisting of an initial configuration x 0 and initial velocity y 0 ):
z ( t , z 0 , u ) π T , F , G z ˜ ( t , z ˜ 0 , u ˜ ) π x ( t , z 0 , u ) T , F , G x ˜ ( t , z ˜ 0 , u ˜ )
where π : R n × R n R n is the canonical projection π ( x , y ) = x , which assigns to the pair ( x , y ) the point x at which the velocity y is attached.

3. Classification of Controllable Systems (2)

A linear mechanical control system (2) is controllable if for any t 0 , any initial state z 0 = ( x 0 , y 0 ) , and any final state z f = ( x f , y f ) there exist t f > t 0 and a control u : [ t 0 , t f ] R m , such that z ( t 0 ) = ( x 0 , y 0 ) and z ( t f ) = ( x f , y f ) . The following result is a straightforward generalization of [10], see also [9] and cf. a classical work on modal controllability for Lagrangian class including dissipative forces [7,14].
Lemma 1 
(Controllability of (2)). For (2) the following statements are equivalent:
(i) 
LMS is controllable
(ii) 
rank B ^ , A ^ B ^ , , A ^ 2 n 1 B ^ = 2 n (Kalman Rank Condition)
(iii) 
rank B ^ , A ^ 2 B ^ , , A ^ 2 ( n 1 ) B ^ = n
(iv) 
rank B , E B , , E n 1 B = n (Mechanical Kalman Rank Condition)
Proof. 
From the Kalman controllability result, we have ( i ) ( i i ) . The rest of the proof follows from a direct computation of the Kalman controllability matrix
B ^ , A ^ B ^ , , A ^ 2 n 1 B ^ = 0 B 0 E B 0 E n 1 B B 0 E B 0 E n 1 B 0 .
Therefore we see that we can take only even powers A ^ 2 i B ^ in ( i i i ) or the lower part of the matrix (10) as in ( i v ) . □
Remark 2.
Note that the Machanical Kalman Rank Condition, i.e., item ( i v ) of Lemma 1, uses objects on the configuration space R n only, while the state-space of (2) is R n × R n . Apart of a pure mathematical value, that sort of reduction is practically motivated since computations are simpler. All of our further results share that property.
Now attach to the system (2) an n-tuple of indices ( r ¯ 0 , , r ¯ n 1 )
r ¯ 0 = rank B , r ¯ i = rank B , E B , . . . , E i B rank B , E B , . . . , E i 1 B ,
for 1 i n 1 . Furthermore, define the dual indices
ρ ¯ j = card r ¯ i j : 0 i n 1 for 1 j m .
These integers, which we will call mechanical half-indices, are mechanical analogues of the controllability (Brunovský, Kronecker) indices r i ’s and ρ j ’s of system (3) defined by (11) and (12), respectively, with E replaced by A ^ and B by B ^ . Note that the indices ρ ¯ i are invariant under (7) and (8), therefore they form a set of invariants attached to (2). Actually, they form a set of complete invariants, as we will show in Theorem 1 below. We denote the above sequences as R ¯ E , B = r ¯ 0 , , r ¯ n 1 and P ¯ E , B = ρ ¯ 1 , , ρ ¯ m .
Proposition 2.
For the mechanical control system LMS
(i) 
the sequence of indices R A ^ , B ^ = r 0 , r 1 , , r 2 n 1 is the doubled sequence of R ¯ E , B , i.e., r 0 , r 1 , , r 2 n 1 = r ¯ 0 , r ¯ 0 , r ¯ 1 , r ¯ 1 , , r ¯ n 1 , r ¯ n 1 ;
(ii) 
the mechanical half-indices are half of the controllability indices, i.e., ρ j = 2 ρ ¯ j , for 1 j m .
Proof. 
Let us invoke the controllability matrix of LMS given by (10) and calculate r 0 = rank B ^ and r i = rank B ^ , , A ^ i B ^ rank B ^ , , A ^ i 1 B ^ , and compare them with r ¯ i given by (11). The crucial observation is that we can calculate the ranks of the lower and upper submatrices separately and then add them.
r 0 = rank B ^ = rank B = r ¯ 0 r 1 = rank B ^ , A ^ B ^ rank B ^ = rank B + rank B rank B = r ¯ 0 .
Assume that ( r 0 , r 1 , , r 2 i 2 , r 2 i 1 ) = r ¯ 0 , r ¯ 0 , , r ¯ i 1 , r ¯ i 1 . Then
r 2 i = rank B ^ , , A ^ 2 i B ^ rank B ^ , , A ^ 2 i 1 B ^ = rank B , , E i B + rank B , , E i 1 B 2 rank B , , E i 1 B = rank B , , E i B rank B , , E i 1 B = r ¯ i , r 2 i + 1 = rank B ^ , , A ^ 2 i + 1 B ^ rank B ^ , , A ^ 2 i B ^ = 2 rank B , , E i B rank B , , E i B + rank B , , E i 1 B = rank B , , E i B rank B , , E i 1 B = r ¯ i .
By an induction argument, the sequence of n integers R ¯ E , B = r ¯ 0 , r ¯ 1 , , r ¯ n 1 and that of 2 n integers R A ^ , B ^ = r 0 , , r 2 n 1 satisfy the desired relation ( i ) . Using (12) and ( i ) , calculate, for 1 j m ,
ρ j = card r i j : 0 i 2 n 1 = 2 card r ¯ i j : 0 i n 1 = 2 ρ ¯ j ,
which proves ( i i ) . □
The following theorem asserts that mechanical half-indices form a set of complete invariants ρ ¯ j of linear controllable mechanical systems (2), as do controllability indices ρ j and indices r ¯ i .
Theorem 1.
The following statements are equivalent, for fixed n and m:
(i) 
Two controllable systems (2) and ( LMS ˜ ) , represented by pairs ( E , B ) and ( E ˜ , B ˜ ) , respectively, are LMF-equivalent,
(ii) 
R ¯ ( E , B ) = R ¯ ( E ˜ , B ˜ ) ,
(iii) 
P ¯ ( E , B ) = P ¯ ( E ˜ , B ˜ ) , i.e., the mechanical half-indices coincide,
(iv) 
P ( A ^ , B ^ ) = P ( A ˜ ^ , B ˜ ^ ) , i.e., the controllability indices coincide,
where A ^ , B ^ and A ˜ ^ , B ˜ ^ are of the form (4).
Proof. 
Equivalence of ( i i ) and ( i i i ) follows from the definition.
( i ) ( i i i ) . We associate with (2), given by the pair ( E , B ) , a virtual linear (first-order) control system
x ˙ = E x + B v ,
and similarly with ( LMS ˜ ) , given by ( E ˜ , B ˜ ) , we associate
x ˜ ˙ = E ˜ x ˜ + B ˜ v ˜ .
Now we directly use the Brunovský classification theorem [1] to prove that (13) and (14) are equivalent under a transformation x ˜ = T x and feedback v = F x + G v ˜ , if and only if their controllability indices coincide. Note that the controllability indices of (13), respectively of (14), coincide with the mechanical half-indices of associated (2) (respectively ( LMS ˜ ) ). Now notice that x ˜ = T x and v = F x + G v ˜ establish feedback equivalence between (13) and (14) if and only if x ˜ = T x , y ˜ = T y , and u = F x + G u ˜ establish LMF-equivalence between (2) and ( LMS ˜ ) . Therefore ( i ) is equivalent to ( i i i ) . Equivalence of ( i i i ) and ( i v ) follows immediately from Proposition 2. □
Remark 3.
Notice that the general feedback group acting on systems of the form (3) by A ^ S A ^ + B ^ F ^ S 1 , B ^ S B ^ G ^ (where S : R 2 n R 2 n is any, not necessarily of the form (7)) is much bigger than the mechanical feedback group (7) and (8). Nevertheless both group actions have exactly the same orbits when acting on linear mechanical systems (2) and thus, the same sets of complete invariants implying that if two linear mechanical systems are feedback equivalent they are also mechanical feedback equivalent.
We can formulate the following important corollary.
Corollary 1.
Any linear mechanical controllable system (2) is LMF-equivalent to the mechanical canonical form
x ˙ i j = y i j y ˙ i j = x i j + 1 y ˙ i ρ ¯ i = u i , , 1 i m , 1 j ρ ¯ i ,
where ( ρ ¯ 1 , , ρ ¯ m ) are mechanical half-indices.
The above corollary follows from the fact (see the proof of the equivalence ( i ) ( i i i ) of Theorem 1) that if a transformation x ˜ = T x , v = F x + G v ˜ brings (13) into the Brunovský form, then x ˜ = T x , y ˜ = T y , and u = F x + G u ˜ brings (2) into the above mechanical canonical form (15).
The mechanical canonical form (15) consists of m chains of even number 2 ρ ¯ i of integrators and can also be represented in the matrix form:
x ˙ = y y ˙ = E F x + B F u ,
where the pair ( E F , B F ) is in the Brunovský form, i.e., E F , B F , are block diagonal matrices, of dimension n × n and n × m , respectively, of the following forms:
E F = N 1 N 2 N m , B F = b 1 b 2 b m ,
where N i , 1 i m , is a superdiagonal nilpotent matrix of dimension ρ ¯ i (mechanical half-index) and b i is a ρ ¯ i × 1 vector:
N i = 0 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 , b i = 0 0 0 1 .

4. Classification of Uncontrollable Systems (2)

In this section, we assume that (2) is not controllable (cf. Lemma 1), i.e.,
rank B , E B , , E n 1 B = k < n .
We will use the same class of mechanical transformations (7) and (8) to establish a canonical form for systems (2) that are not controllable and thus we will classify them (for the controllable case, Theorem 1 and form (15) provide a complete classification).
Theorem 2.
For given n , m , k , the linear mechanical system (2) is LMF-equivalent to the following canonical form:
x ˙ c = y c y ˙ c = E F x c + B F u x ˙ d = y d y ˙ d = E d J x d ,
where dim x c = k , and E F , B F are block-diagonal matrices of dimension k × k , k × m , respectively, of the form (16) and (17), given by the mechanical half-indices ρ ¯ 1 , , ρ ¯ m satisfying j = 1 m ρ ¯ j = k , while E d J is in the Jordan form, that is, a block-diagonal ( n k ) × ( n k ) matrix whose diagonal blocks are of four possible forms:
(i)
for a real eigenvalue λ j of E d J
either D j R = λ j λ j λ j or J j R = λ j 1 λ j 1 λ j ,
(ii)
for a complex eigenvalue λ j = α j ± i β j of E d J
either D j C = C j C j C j or J j C = C j I 2 C j I 2 C j ,
where C j = α j β j β j α j and I 2 = 1 0 0 1 .
Therefore any mechanical system (2) can be decomposed under LMF into two independent mechanical subsystems:
(i)
a 2 k -dimensional controllable mechanical system (2) represented in ( x c , y c ) -coordinates by the pair ( E F , B F ) , which is in the canonical form (15), with mechanical half-indices ρ ¯ 1 , , ρ ¯ m ; the subindex “c” stands for controllable,
(ii)
a dynamical linear mechanical system (a system without controls) represented in ( x d , y d ) -coordinates by the second-order differential equation x ¨ = E d x d , where the matrix E d can always be transformed into the Jordan form E d J ; the subindex “d” stands for dynamical.
Proof. 
The system (2) is uniquely defined by the pair ( E , B ) on which an element T , F , G of the group L M F acts according to the rule:
E , B ( T , F , G ) T ( E + B F ) T 1 , T B G ,
which is the same equivalence transformation as in [5], where a slightly different notation ( A = E , P = T 1 , K = F T 1 , Q = G ) is used. By the result of [5], we bring the pair ( E , B ) into the form:
E = E F 0 0 E d , B = B F 0 ,
where E d is an ( n k ) × ( n k ) matrix and ( E F , B F ) is in Brunovský form. What remains to prove is to bring E d into its canonical form E d J . In order to do that, we use a change of coordinates x d T d x d , y d T d y d , where T d : R n k R n k brings E d into its real Jordan form E d J , see e.g., [15], and apply T = I k 0 0 T d . □
Denote B = Im B and X c = B + E B + + E k 1 B . We have E X c X c (since rank ( B , E B , , E k 1 B ) = k ), so E gives rise to a well defined map E ^ : R n / X c R n / X c between the quotient spaces (i.e., factor spaces). Choose a subspace X d R n such that X c X d = R n and any linear coordinates ( x c , x d ) on R n such that X c = x c 0 : x c R k and X d = 0 x d : x d R n k . Thus the matrix E takes in the ( x c , x d ) -coordinates the block-triangular form E = E c 1 E c 2 0 E d and the map E ^ in these coordinates is given by the matrix E d . Choosing another subspace X ˜ d (completing X c to R n ) will lead to another system of coordinates ( x ˜ c , x ˜ d ) and a matrix E ˜ d related to E d by E ˜ d = T E d T 1 . So the eigenvalues of E ^ are well defined as those of E d (being, obviously, the same as those of E ˜ d ).
Let λ 1 , , λ q be q mutually distinct eigenvalues of the matrix E d . For each λ j , 1 j q , let δ j be the dimension of the diagonal block of λ j and let κ j be the number of Jordan blocks of λ j , whose dimensions are ε j 1 , , ε j κ j , respectively. We define the eigenstructure Λ j of λ j by Λ j = ( λ j , δ j , ε j 1 , , ε j κ j ) . Clearly, the same eigenvalue λ j may appear in more than one block and the algebraic multiplicity of λ j is μ j = δ j + ε j 1 + + ε j κ j .
Proposition 3.
Two mechanical control systems (2) and ( LMS ˜ ) given by ( E , B ) and ( E ˜ , B ˜ ) , respectively, are LMF-equivalent if and only if their mechanical half-indices coincide, that is P ¯ ( E , B ) = P ¯ ( E ˜ , B ˜ ) , and the eigenstructures of E d and E ˜ d coincide ( Λ 1 , , Λ q ) = ( Λ ˜ 1 , , Λ ˜ q ) , up to a permutation of Λ j ’s, where q, the number of distinct eigenvalues, is the same for both systems.
Proof. 
By Theorem 2, both systems can be brought, by LMF-transformations, into their canonical form (18). Then, both controllable subsystems are LMF-equivalent if and only if their mechanical half-indices coincide (see Theorem 1). Finally, if the eigenstructures of E d J and E ˜ d J coincide, then they are the same matrices up to a permutation of blocks (note that both are in the Jordan form) and vice versa. Thus the composition of both LMF-transformations converts (2) into ( LMS ˜ ) . □
Example 1.
Consider two systems (2) and ( LMS ˜ ) , with n = 3 , m = 1 represented, respectively, by the following pairs ( E , B ) and ( E ˜ , B ˜ ) :
E = 1 0 0 0 1 0 0 0 1 , E ˜ = 1 1 0 0 1 0 0 0 1 , B = B ˜ = 0 0 1 .
It can be easily checked that the two systems are not LMF-equivalent. For both systems k = 1 and the mechanical half-indices are ρ ¯ 1 = 1 . The dynamical (uncontrolled) part of both systems has λ = 1 as an eigenvalue of algebraic multiplicity 2 but the eigenstructure of the first system is Λ = ( λ , δ , ε ) = ( 1 , 2 , 0 ) , meaning that λ = 1 defines a diagonal block of dimension 2, and that of the second system is Λ = ( λ ˜ , δ ˜ , ε ˜ ) = ( 1 , 0 , 2 ) , meaning that there is a single Jordan block of dimension 2.

5. Classification of Lagrangian Systems (6)

We start with a counterpart of Proposition 3. To formulate it, recall that the matrix E L = M 1 P , associated with the Lagrangian system (6), has real eigenvalues and is diagonalizable. It follows that, each eigenvalue λ j of E d L is real and its eigenstructure is Λ j = ( λ j , δ j , 0 ) , for 1 j q , where q is the number of mutually distinct eigenvalues of E d L and μ j = δ j are their algebraic multiplicities. In this case, we will simply denote Λ j = ( λ j , μ j ) . Recall that, see Section 4, the matrix E d represents the quotient map E ^ : R n / X c R n / X c , induced by E : R n R n , where X c = B + + E k 1 B .
Proposition 4.
Two Lagrangian systems (6) and ( LLS ˜ ) , given by ( E L , B ) and ( E ˜ L , B ˜ ) , respectively, are LMF-equivalent if and only if their mechanical half-indices coincide, that is P ¯ ( E L , B ) = P ¯ ( E ˜ L , B ˜ ) and, up to permutations, ( λ j , μ j ) = ( λ ˜ j , μ ˜ j ) , for 1 j q , where λ j ’s and λ ˜ j ’s are distinct eigenvalues of E d L and E ˜ d L , respectively, and q is their number (the same for both systems).
A proof follows directly from Proposition 3 applied to two Lagrangian systems. The statements of Propositions 3 and 4 are formally the same; the only difference is that the eigenstructure of E d of (2) can have any elements Λ j (real or complex eigenvalues λ j and diagonal or Jordan blocks), while Λ j of E d L of (6) consists of real eigenvalues λ j and diagonal blocks only. Contrary to Proposition 4, that is a Lagrangian counterpart of Proposition 3, the classification Theorem 2 does not apply to the Lagrangian systems (6) because, in general, the feedback transformation u = F x + G u ˜ does not perserve the Lagrangian structure. Actually, the mechanical canonical form, given by ( E F , B F ) , is never Lagrangian because E F has all eigenvalues λ j = 0 , and thus never can be written as a product of two symmetric matrices E F = M 1 P , for M being invertible. Indeed, a symmetric matrix, whose all eigenvalues are zero, is P = 0 but then M 1 P = 0 E F .
Instead of Theorem 2, we have the following result.
Proposition 5.
(i)
Any Lagrangian mechanical system (6) is LMF-equivalent to the following Lagrangian system:
x ˙ c = y c y ˙ c = E F L x c + B F u x ˙ d = y d y ˙ d = E d L x d ,
where E d L = diag ( λ 1 d , , λ l d ) , l = n k , with λ j d arbitrary real, and E F L is of the same form as E F , given by (16) and (17), with N i replaced by L i = 0 1 0 0 0 1 a 1 i a 2 i a ρ ¯ i i , such that the eigenvalues of L i are all real and mutually distinct.
(ii)
Any two Lagrangian systems of the form (19), with the same E d L but arbitrary terms a j i and a ˜ j i (such that, for each block 1 i m , the eigenvalues λ j i of L i are mutually distinct and so are λ ˜ j i of L ˜ i ), are LMF-equivalent.
Proof. 
( i ) . By Theorem 2, any (6) is LMF-equivalent to (18), whose x c -subsystem is not Lagrangian, for which E d J is E d L (since all eigenvalues of (6) are real). Then we change N i into L i using the feedback transformation that replaces u i by j = 1 ρ ¯ i a j i x c j + u i , such that the eigenvalues of L i are real and mutually distinct, and obtain a Lagrangian system of the form (19). To prove ( i i ) , take any system of the form (19) and using the feedback transformation that replaces u i by j = 1 ρ ¯ i ( a ˜ j i a j i ) x c j + u i , we get a system of the form (19) in which a ˜ j i take place of a j i . □
The above proposition suggests that there is no a Brunovský-like canonical form for Lagrangian systems. On one hand, the matrix L i with all a j i = 0 does not give a Lagrangian system, on the other hand, there is no a privileged choice of non-zero eigenvalues, and all of them are feedback equivalent as asserted by item ( i i ) . This suggests not to change the original eigenvalues and to apply a change of coordinates x ˜ = T x , y ˜ = T y only, which we treat in the next proposition.
Proposition 6.
Consider a single-input, m = 1 , Lagrangian mechanical system. It is equivalent via x ˜ = T x , y ˜ = T y to
x ˙ c i = y c i y ˙ c i = λ i c x c i + u 1 i k x ˙ d i = y d i y ˙ d i = λ i d x d i , 1 i l = n k
where all λ i c , λ i d R and, moreover, λ i c are mutually distinct.
The above form is clear. The matrix E is diagonal and the system (6) decouples into two independent subsystems: the uncontrollable one consisting of l = n k independent second order dynamical systems x ¨ d i = λ i d x d i (with real eigenvalues λ i d that can be any) and a completely controllable one consisting of k second order control systems x ¨ c i = λ i c x c i + u , whose eigenvalues λ i c are distinct, all controlled by the same control u. A similar form can be obtained for the multi-input case m > 1 (with a more complicated form of the matrix B) but we will not present it here because of lack of space.

6. Stability and Stabilization

Now we come back to (not necessarily Lagrangian) mechanical systems (2) and we turn our attention to the relation between the eigenvalues of E, denoted λ j , and the eigenvalues of A ^ , denoted σ j , since the latter are responsible for the stability of the system.
Lemma 2.
Let A ^ = 0 I n E 0 . If λ is an eigenvalue of E, then σ = ± λ are eigenvalues of A ^ .
Proof. 
By a direct calculation, we see that the characteristic polynomial of A ^ is
det ( σ I 2 n A ^ ) = det σ I n 0 0 σ I n 0 I n E 0 = det ( σ 2 I n E ) = ( σ 2 λ 1 ) ( σ 2 λ 2 ) ( σ 2 λ n ) ,
since all blocks of A ^ commute, and where λ 1 λ n are the eigenvalues of E. The above polynomial has the roots σ j = ± λ j . □
If λ R , then σ = ± a or σ = ± a i , where a = | λ | R + . If λ = α + β i C (and its conjugate α β i C ), then σ = ± ( a + b i ) and ± ( a b i ) , where α = a 2 b 2 and β = 2 a b . Denote R 0 = λ R : λ 0 . Let us visualize the relation between the eigenvalues with the following table.
λ ( E ) σ ( A ^ ) Sketch
R + ± a Applsci 11 10669 i001
R 0 ± a i Applsci 11 10669 i002
C ± ( a + b i ) and ± ( a b i ) Applsci 11 10669 i003
The above analysis of the eigenvalues leads to the following simple but important observation. We will say that a dynamical system z ˙ = A ^ z , for A ^ given by (4) and z = ( x , y ) , is stable if all its equilibrium points are stable. Recall that μ j denotes the algebraic multiplicity of an eigenvalue λ j and γ j = dim ker ( λ j I n E ) its geometric multiplicity.
Proposition 7.
Consider the dynamical system z ˙ = A ^ z , where A ^ = 0 I n E 0 or, equivalently, x ¨ = E x .
(AS) 
The system z ˙ = A ^ z is never asymptotically stable.
(S) 
The following conditions are equivalent:
(i) 
the system z ˙ = A ^ z is stable,
(ii) 
all eigenvalues λ j of E satisfy λ j R 0 and, moreover, their algebraic and geometric multiplicities coincide, i.e., μ j = γ j , for 1 j q , where q is the number of distinct eigenvalues of Ed,
(iii) 
E has a Lagrangian structure, i.e., E = M 1 P for some symmetric matrices M (invertible) and P and the eigenvalues of E satisfy λ j R 0 .
Proof. 
It is immediate to see ( A S ) , since at least for one eigenvalue σ j of A ^ we have Re σ j 0 . Equivalence of ( i ) ( i i ) follows immediately from Lemma 2 and the table above. For a proof of necessity of μ j = γ j (which is equivalent to rank ( λ j I n E ) = n μ j ), see the proof of Theorem 4.5 in [16]. Equivalence ( i i ) ( i i i ) is given by Proposition 1 and the comments below it. □
Of course, the class of second-order differential equations (or dynamical systems) x ¨ = E x is bigger than the class of Lagrangian systems x ¨ = E L x because the matrix E can be any, while E L has to be R -diagonalizable. It is obvious that x ¨ = E L x cannot be asymptotically stable because it preserves the energy 1 2 x ˙ T M x ˙ + 1 2 x T P x . It turns out that x ¨ = E x is never asymptotically stable either (although it may have, contrary to x ¨ = E L x , complex eigenvalues with non-zero real part) and, moreover, it is stable if and only if it is Lagrangian. So there are no stable second order differential equations (in other words, dynamical systems) x ¨ = E x others than Lagrangian ones.
A mechanical control system (2) is called asymptotically stabilizable if there exists a mechanical feedback of the form u = F x such that the closed loop
x ˙ = y y ˙ = ( E + B F ) x
is asymptotically stable and is called stabilizable if all equilibria of the above closed loop system are stable. Recall that the map E : R n R n induces the map E ^ : R n / X c R n / X c , where X c R k , that is represented by the matrix E d R l × l , where l = n k .
Proposition 7, applied to the uncontrolled system x ˙ d = y d , y ˙ d = E d x d , leads to the following result describing stabizability of mechanical control systems (2).
Theorem 3.
Consider a mechanical control system (2).
(AS) 
The system (2) is never asymptotically stabilizable.
(S) 
The following conditions are equivalent:
(i) 
(2) is stabilizable,
(ii) 
the matrix E d has all eigenvalues λ j R 0 and, moreover, is diagonalizable, i.e., its Jordan form E d J consists of diagonal blocks D j R only, corresponding to λ j R 0 ,
(iii) 
all eigenvalues λ j of E d satisfy λ j R 0 and, moreover, their algebraic and geometric multiplicities coincide, i.e., μ j = γ j , for 1 j q , where q is the number of distinct eigenvalues of E d ,
(iv) 
the uncontrolled subsystem x ˙ d = y d , y ˙ d = E d x d is T-equivalent to a Lagrangian system,
(v) 
(2) is LMF-equivalent to (6) and, moreover, the matrix E d satisfies the conditions of item ( i i ) above (or of the equivalent items ( i i i ) or ( i v ) ).
Proof. 
Item ( A S ) follows directly from Proposition 7. For ( S ) , the crucial observation is that the controllable system is stabilizable, thus we deal with the uncontrollable system only, which must already be stable for the whole system to be stabilizable. Therefore, the equivalence ( i ) ( i i i ) follows directly from Proposition 7. Equivalence of ( i i ) ( i i i ) is given by Proposition 1 and the comments below it. For ( i i ) ( i v ) , note that the matrix E d J of ( i i ) coincides with E d L of (19), so the transformation T that diagonalizes E d renders the uncontrolled subsystem Lagrangian. Finally, assume that (2) is LMF-equivalent to (6) and E d satisfies ( i i ) , or equivalent conditions ( i i i ) or ( i v ) , i.e., it is stable. Then, by Proposition 5 ( i i ) it is also LMF-equivalent to a Lagrangian system that is stable. Hence, LMF-transformation stabilizes the system (2). The inverse follows from the previous arguments. □
Corollary 2.
Any stabilizing feedback u = F x , for a mechanical control system (2), renders the system x ˙ = y , y ˙ = ( E + B F ) x Lagrangian.
In other words, for (2) to be stabilizable, all unstable modes x ¨ j = λ j x j , where Re λ j > 0 must be controllable, i.e., must be contained in the controllable x c -subsystem, or, equivalently, all uncontrollable modes must be stable. Checking that requires, however, a decomposition of (2) into controllable and uncontrollable subsystems. Therefore we provide below invariant conditions that can be verified for any (2). This is analogous to [7], where Lagrangian systems with dissipative forces were considered.
Proposition 8.
The system (2), given by ( E , B ) , is stabilizable if and only if for any eigenvalue λ j of E such that rank ( λ j I n E , B ) < n , we have λ j R 0 and μ j = γ j , where μ j and γ j are, respectively, the algebraic and geometric multiplicity of λ j .
Proof. 
Note that rank ( λ j I n E ) < n if and only if λ j is an eigenvalue of E. By Hautus lemma (see e.g., [7,15,17]), if rank ( λ j I n E , B ) = n , then the mode x ¨ j = λ j x j is controllable. Thus, by assuming rank ( λ j I n E , B ) < n we identify all uncontrollable modes so the corresponding subsystems have to be stable, i.e., λ j R 0 and μ j = γ j , see Theorem 3 ( i i ) . □

7. Examples

Classical examples of linear mechanical control systems are mass-spring systems. We present the equations of motions of n-coupled mass-spring system, which consists of n bodies, where the position of i-th body is denoted x i , and m i is the mass of i-th body. The bodies are connected by n + 1 springs with k i being the spring constant of i-th spring, as depicted in Figure 1. The external forces (controls) u i may, á priori, be applied to each body.
The dynamics of i-th body is given by the balance of forces acting on the body
m i x ¨ i = k i x i x i 1 + k i + 1 x i + 1 x i + u i ,
where x 0 x n + 1 0 . The equations can be formulated in the form of (5), where -4.6cm0cm
M = m 1 0 0 0 0 m 2 0 0 0 0 m 3 0 0 0 0 0 0 0 0 m n P = k 1 + k 2 k 2 0 0 k 2 k 2 + k 3 k 3 0 0 k 3 k 3 + k 4 0 0 0 0 k n 0 0 k n k n + k n + 1
or as a Lagrangian system of the form (6), where
E = E L = k 1 k 2 m 1 k 2 m 1 0 0 k 2 m 2 k 2 k 3 m 2 k 3 m 2 0 0 k 3 m 3 k 3 k 4 m 3 0 0 0 0 k n m n 1 0 0 k n m n k n k n + 1 m n B = 1 m 1 0 0 0 0 1 m 2 0 0 0 0 1 m 3 0 0 0 0 0 0 0 0 1 m n .
If all n controls are present, then the system is fully actuated. However, it is enough to apply only the control u n in order for the system to be controllable.
Example 2
(The mass-spring system with one control). Consider the n-mass-spring system (21) with only one control u : = u n , i.e u i = 0 for 1 i n 1 . The system is (2), where E is given by (22) and b = 0 , , 1 m n T . It is straightforward to show that rank b , E b , , E n 1 b = n . What is more, introduce c = ( m ¯ k ¯ , 0 , 0 , , 0 ) , where m ¯ = i = 1 n m i and k ¯ = i = 2 n k i , and apply the transformation T = ( c , c E , , c E n 1 ) T . The transformed system is a Lagrangian system of the form (19), given by:
E F L = T E T 1 = 0 1 0 0 0 1 a 1 a 2 a n and b F = T b = 0 0 1
and E d L is nonexistent since the system is controllable. Note that, by applying the feedback u = i = 1 n a i x ˜ i + u ˜ , we obtain the (non Lagrangian) mechanical canonical form (15) or we substitute a i ’s with any other a ˜ i ’s, and those a ˜ i ’s, for which the corresponding eigenvalues λ ˜ j satisfy λ ˜ j R 0 and are mutually distinct, give stable Lagrangian systems.
Example 3
(The uncontrollable mass-spring system). Consider the mass-spring system with 3 equal masses m and 4 equal springs with the spring constant k. The external force is applied to the second mass only. The equations of motion are of the form (6) with:
E = E L = 2 k m k m 0 k m 2 k m k m 0 k m 2 k m and b = 0 1 m 0 .
A direct calculation shows that rank ( b , E b , E 2 b ) = 2 < n = 3 . In order to decompose the system, we take T 1 = col ( b , E b , v ) , where v = ( 1 , 0 , 0 ) T and set:
E ˜ = T E T 1 = 0 2 k 2 m 2 k 1 4 k m 0 0 0 2 k m and b ˜ = T b = 1 0 0 .
Then, by another transformation T ¯ = 1 k m ( 2 + 2 ) m 2 1 k m ( 2 2 ) m 2 0 0 1 applied to the pair ( E ˜ , b ˜ ) , we can establish a Lagrangian system in the form (20) with coordinates ( x c 1 , x c 2 , x d 1 )
x ˙ c 1 = y c 1 y ˙ c 1 = λ 1 c x c 1 + u x ˙ c 2 = y c 2 y ˙ c 2 = λ 2 c x c 2 + u x ˙ d 1 = y d 1 y ˙ d 1 = λ 1 d x d 1 ,
where λ 1 c = k m ( 2 + 2 ) , λ 2 c = k m ( 2 2 ) and λ 1 d = 2 k m . It is immediate to see that the system is stable since λ 1 d , λ 1 c , λ 2 c R 0 and all eigenvalues σ j of A ^ are pure imaginary giving oscillations. Moreover, the frequency | λ 1 d | is invariantly related to the system, while λ 1 c , λ 2 c can be set freely (by an appropriate choice of feedback u = F x c ).

8. Conclusions

In this paper, we have studied a classification of linear mechanical control systems (2) under mechanical change of coordinates and feedback. In Section 3, Section 4 and Section 5, which are the heart of our paper, we completely solved the classification problem for both controllable and uncontrollable cases, which enabled us to establish the corresponding canonical forms. To our best knowledge this complete solution, expressing complete invariants in terms of objects on the configuration space only, have not been publish before and constitutes a novelty of our work. Our obvious inspirations are Lagrangian control systems, which apart from applicational importance, turn out to be crucial in Section 6, where we have discussed stability and stabilization of (2) (not asymptotic stabilization). The conclusion is that there are no other stabilizable mechanical control systems (2) than Lagrangian ones. Finally, we illustrated our results by a classical representative of linear mechanical control system, i.e., by the mass-spring system, which is considered in two special cases: one that is controllable and another that is not.

Author Contributions

Conceptualization, methodology and formal analysis, M.N. and W.R.; writing—original draft preparation, M.N.; writing—review and editing, W.R. and M.N. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partially supported by the statutory grant No. 0211/SBAD/0520.

Acknowledgments

The authors dedicate this work to the memory of Krzysztof Kozłowski.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Brunovský, P. A Classification of Linear Controllable Systems. Kybernetika 1970, 6, 173–188. [Google Scholar]
  2. Rosenbrock, H.H. State-Space and Multivariable Theory, 1st ed.; Nelsons: London, UK, 1970. [Google Scholar]
  3. Kalman, R.E. Mathematical description of linear dynamical systems. SIAM J. Control Optim. 1963, 1, 152–192. [Google Scholar] [CrossRef]
  4. Brockett, R.W. Finite Dimensional Linear Systems; Wiley: New York, NY, USA, 1970. [Google Scholar]
  5. Liang, B.C. A Classification of Linear Control Systems. Linear Multilinear Algebra 1980, 8, 261–264. [Google Scholar] [CrossRef]
  6. Bulatovic, R.M. A stability theorem for gyroscopic systems. Acta Mechanica 1999, 136, 119–124. [Google Scholar] [CrossRef]
  7. Laub, A.J.; Arnold, W.F. Controllability and Observability Criteria for Multivariable Linear Second-Order Models. IEEE Trans. Automat. Contr. 1984, 29, 163–165. [Google Scholar] [CrossRef]
  8. Müller, P.C. Stability and Optimal Control of Nonlinear Descriptor Systems: A Survey. Appl. Math. Comp. Sci. 1998, 8, 269–286. [Google Scholar]
  9. Bullo, F.; Lewis, A.D. Geometric Control of Mechanical Systems Modeling, Analysis and Design for Simple Mechanical Control Systems, 1st ed.; Springer: New York, NY, USA; Heidelberg/Berlin, Germany, 2004. [Google Scholar]
  10. Hughes, P.C.; Skelton, R.E. Controllability and Observability of Linear Matrix-Second-Order Systems. J. Appl. Mech. 1980, 47, 415–420. [Google Scholar] [CrossRef]
  11. Nowicki, M.; Respondek, W. Mechanical state-space linearization of mechanical control systems and symmetric product of vector fields. In Proceedings of the 7th IFAC LHMNC 2021, Berlin, Germany, 11–13 October 2021. [Google Scholar]
  12. van der Schaft, A. Controlled Invariance for Hamiltonian Systems. Math. Syst. Theory 1985, 18, 257–291. [Google Scholar] [CrossRef] [Green Version]
  13. van der Schaft, A. Port-Hamiltonian systems: An introductory survey. In Proceedings of the International Congress of Mathematicians, Madrid, Spain, 22–30 August 2006; Volume III, pp. 1339–1365. [Google Scholar]
  14. Belotti, R.; Richiedei, D.; Tamellin, I.; Trevisani, A. Pole Assignment for Active Vibration Control of Linear Vibrating Systems through Linear Matrix Inequalities. Appl. Sci. 2020, 10, 5494. [Google Scholar] [CrossRef]
  15. Zabczyk, J. Mathematical Control Theory, 2nd ed.; Springer Nature: Cham, Switzerland, 2020. [Google Scholar]
  16. Khalil, H.K. Nonlinear Systems, 3rd ed.; Prentice Hall: Hoboken, NJ, USA, 2002. [Google Scholar]
  17. Hautus, M.L.J. Controllability and observability conditions of linear autonomous systems. Proc. Kon. Ned. Akad. Wetensch. 1969, A, 443–448. [Google Scholar]
Figure 1. The n-coupled mass-spring system.
Figure 1. The n-coupled mass-spring system.
Applsci 11 10669 g001
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nowicki, M.; Respondek, W. A Mechanical Feedback Classification of Linear Mechanical Control Systems. Appl. Sci. 2021, 11, 10669. https://doi.org/10.3390/app112210669

AMA Style

Nowicki M, Respondek W. A Mechanical Feedback Classification of Linear Mechanical Control Systems. Applied Sciences. 2021; 11(22):10669. https://doi.org/10.3390/app112210669

Chicago/Turabian Style

Nowicki, Marcin, and Witold Respondek. 2021. "A Mechanical Feedback Classification of Linear Mechanical Control Systems" Applied Sciences 11, no. 22: 10669. https://doi.org/10.3390/app112210669

APA Style

Nowicki, M., & Respondek, W. (2021). A Mechanical Feedback Classification of Linear Mechanical Control Systems. Applied Sciences, 11(22), 10669. https://doi.org/10.3390/app112210669

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop