Next Article in Journal
A New Explicit Numerical Schemes for Time-Dependent PDEs with Application to Pressure Driven Fluid Flow in a Rectangular Duct
Next Article in Special Issue
Study and Analysis of Dynamics and Energy Efficiency of Arc Steelmaking Furnace Electrical Mode with a Fuzzy Control Algorithm
Previous Article in Journal
The Effect of Microencapsulated PCM Slurry Coolant on the Efficiency of a Shell and Tube Heat Exchanger
Previous Article in Special Issue
Exploiting the Moth–Flame Optimization Algorithm for Optimal Load Management of the University Campus: A Viable Approach in the Academia Sector
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sub-Optimal Stabilizers of the Pendubot Using Various State Space Representations

Institute of Automatic Control and Robotics, Poznan University of Technology, ul. Piotrowo 3a, 60-965 Poznan, Poland
*
Author to whom correspondence should be addressed.
Energies 2022, 15(14), 5146; https://doi.org/10.3390/en15145146
Submission received: 20 June 2022 / Revised: 12 July 2022 / Accepted: 13 July 2022 / Published: 15 July 2022
(This article belongs to the Special Issue Thermo-Mechanical and Electrical Measurements for Energy Systems)

Abstract

:
This paper considers the issue of linear-quadratic regulator (LQR) design for nonlinear systems with the use of smooth state and input transformations. The proposed design methodology is considered in the stabilisation task of the Pendubot, which is based on the concept of feedback equivalent control systems. It turns out that it is possible to find a controller that ensures comparable dynamics of the closed-loop system in the vicinity of the set point regardless of the state-space representation adopted. In addition, the synthesis of suboptimal controllers according to the LQR strategy ensuring equal dynamics at the equilibrium point is presented. The properties of the studied controllers were investigated in a simulation environment and using experimental tests. The detailed forms of transformations and linear approximations given can be regarded as ready-made procedures that can be applied to stabilise similar mechanical systems in robotics.

1. Introduction

The problem of designing optimal feedback to stabilize dynamic systems can be regarded as one of the fundamental issues of control theory. In particular, for a class of time-invariant linear (LTI) systems, this problem is already well recognized. One of the essential control design approaches to stabilize these systems is based on the minimization of an energetic-like quadratic performance index involving the state trajectory and input, which leads to the linear-quadratic regulator commonly known as LQR. This approach has a long and rich history, dating back to the early work of Kalman [1], who discussed the problem of optimal feedback, providing the design equations for LQR.
Feedback design becomes more challenging in the case of nonlinear systems, for which the search for an optimal controller can be a complex problem that cannot be solved by analytical methods. One important tool in control design is to resort to a linear approximation of the nonlinear dynamics around the equilibrium point under condition that the corresponding linear model is controllable in the sense of Kalman. A known limitation of this approach is the local nature of the solution, which imposes restrictions on the set of admissible initial conditions. In such circumstances, the LQR approach is capable of providing only suboptimal results. Despite that, LQR still constitutes an essential method for dealing with the control of nonlinear systems. In particular, it is able to guarantee a local exponential convergence, that provides robustness to some class of uncertainties [2]. For this reason, the method is of great importance in robotics, where the predominant group of mechanical systems exhibit strongly nonlinear properties. Recently, one can find publications that discuss applications of the LQR method, often in the context of complex hybrid techniques for the control of manipulators with rigid or compliant joints [3,4,5]. In particular, LQR-based control solutions have been proposed for a class of underactuated systems [6,7,8,9]. In some cases LQR approach is used as a tuning method for strictly nonlinear controllers, see for instance [10]. Modifications of LQR for nonlinear systems have been considered, such as the state-dependent Riccati equation (SDRE) approach [11] and its various implementations [12,13,14]. In addition, techniques inspired by differential dynamic programming [15,16] lead to iterative LQR (iLQR) [17,18], for which the dynamics are linearized in a sequence and the cost function about a given nominal trajectory is computed to find an LQR control policy. Based on a similar idea, the LQR-equivalent of Kalman smoothing has been reported in [18].
In this paper, we refer to nonlinear control theory and use the concept of feedback-equivalent control systems [19] in order to improve the performance of the LQR design by extending the set of feasible initial conditions that define the so-called basin of attraction. The considered methodology is to use a linear state feedback designed for the linear approximation of a feedback-equivalent control system, and determine a stabilizer taking advantage of state and input maps. A key ingredient is the transformation of the original dynamics to a form that exhibits better characteristics for synthesizing a linear regulator. In particular case, the equivalent dynamics can be even linear [20], which can considerably improve the controller performance in a given subset of the state space at which linearization is possible. In this context, it may be important to estimate the area of convergence. An important tool here is Lapunov analysis [21,22] and the use of nonlinear numerical and analytical methods such as the sum of squares (SOS) [23,24], which allow for a less restrictive approximation.
From the point of view of control objectives, the design of the LQR should take into account the quality index determined for the original system, i.e., taking into account the state of this system and the energy expenditure associated with the actual input. Therefore, in the case of the feedback design based on the transformed system, the question of the LQR tuning criterion seems to be important [25,26]. This issue is analyzed in this work. We show here a solution to determine the gains that ensure comparable dynamics of the closed-loop system in the vicinity of the set point regardless of the state and input representation adopted. In addition, we provide an equivalent quadratic form that describes the optimization criterion with respect to the transformed dynamics.
In this paper, the control methodology considered is used to design a Pendubot stabilizer. This system, along with Acrobot and the inverted pendulum, can be considered as a benchmark underactuated system in robotics [27,28,29]. We deal with the stabilization of the Pendubot at up-right position without taking into account the swing-up problem. Instead, the problem investigated here is focused on the design of a smooth state feedback at a neighborhood of the desired point taking advantage of the concept of feedback equivalent control systems. Although, the linear approximation of the Pendubot at an equilibrium point is controllable, the system is not fully linearizable using state and input transformations. Furthermore, as discussed precisely in [30], there are also significant obstructions in the application of the partial linearization approach, while this method cannot be directly employed for stabilization due to the presence of singular points [31].
To investigate the effect of the state-space representation on the characteristics of the closed-loop system, the original Langrange dynamics of the Pendubot is transformed into two alternative forms. These forms take advantage of the so-called quasi-velocities, which in analytical mechanics are understood as linear combinations of generalized velocities with coefficients that are functions of the generalized coordinates [32,33]. In the first form, the inertial normalized quasi-velocities (NQV) proposed by Jain and Rodriguez in [34,35] which comes from the factorization of the inertia matrix, are used. In particular, the transformations proposed for the Pendubot to design a swing-up controller [36], in this paper are used for the stabilization task. It turns out that quasi-velocities along with the nonlinear transformation of the coordinates can be used to represent the Pendubot dynamics in the so-called normal form [37]. Such a form, considered among others in classification problems, highlights in an organized way the essential features of a nonlinear dynamic system [30].
The new contributions to this paper include the following:
  • Comparison of Pendubot mathematical models using different representations, including application of quasi-velocities;
  • Synthesis of sub-optimal controllers according to the LQR strategy ensuring equal dynamics at the equilibrium point;
  • Simulation comparison of the controllers and determination of the convergence area under constrained input conditions;
  • Conducting experimental tests and obtaining results illustrating properties of controllers.
It is noteworthy that another purpose of the work is to adapt the control algorithm to the real system and to show reproducible results. To the authors’ knowledge, in many publications on stabilization of Acrobot and Pendubot-type mechanical structures, real-based models are not explicitly investigated and primarily the basic models proposed in Spong’s works are recalled. Often, other works do not contain sufficient information about the model and controller parameters, or are tested for a non-physical system. Our aim is to dispel these doubts through experimental verification and research on stability or the area of convergence of the algorithm. Therefore, the experiments have been carried out for a system that can be built from components of a commercially available system. For this reason, the results shown in the paper can be treated as a basis for future comparisons.
The paper is organized as follows. In Section 2 the nominal dynamics of the Pendubot are recalled. Section 3 deals with feedback-equivalent control systems and describes two equivalent models of the Pendubot taking into account the inertial normalized quasi-velocities (NQV) and transformation to the normal form (NF). In Section 4, the design of the controller based on the LQR approach and its stability issues are discussed. In Section 5 simulation and experimental results are presented. Section 6 discusses the results obtained, and Section 7 ends with general conclusions and plans for future research.

2. Model

Let us consider the mechanical system presented in Figure 1 with the state chosen as x = q T ω T T S 1 × S 1 × R 2 , where q = q 1 q 2 T stands for the system configuration and ω = ω 1 ω 2 T denote velocities, while q i and ω i describe the angular displacement and angular velocity of the i th link ( i = 1 , 2 ), respectively. The system input u = τ R corresponds to the torque exerted on the first joint. In the state-space representation, the dynamics of the Pendubot can be described as follows
x ˙ = ω D 1 q C q , q ˙ q ˙ G q + b u ,
where
D ( q ) = a 1 + a 2 + 2 a 3 cos q 2 a 2 + a 3 cos q 2 a 2 + a 3 cos q 2 a 2
is the inertia matrix,
C ( q , q ˙ ) = a 3 sin q 2 q ˙ 2 a 3 sin q 2 ( q ˙ 1 + q ˙ 2 ) a 3 sin q 2 q ˙ 1 0
is the centrifugal and Coriolis matrix,
G q = a 5 cos ( q 1 + q 2 ) + a 4 cos q 1 a 5 cos ( q 1 + q 2 ) ,
describes gravity torques and b = 1 0 T is the input matrix. The constant coefficients a j , j = 1 , 2 , , a 5 , satisfy: a 1 = m 1 l c 1 2 + m 2 l 1 2 + I 1 , a 2 = m 2 l c 2 2 + I 2 , a 3 = m 2 l 1 l c 2 , a 4 = g m 1 l c 1 + m 2 l 1 and a 5 = g m 2 l c 2 .

3. Equivalent Models of the Pendubot in Quasi-Velocities

3.1. Feedback Equivalent Control Systems

In this paper, we deal with feedback-equivalent control systems [19]. To describe the concept formally, let us introduce the following smooth dynamic systems.
Σ : x ˙ = f x , u ,
Σ * : ξ ˙ = f * ξ , ν ,
where x X R n , ξ Ξ R n denote states, while u , ν R m are inputs and f : X × R m R n , f * : Ξ × R m R n are smooth state functions. It is assumed that these systems are equivalent under the change of states,
ξ = ρ x
and inputs,
ν = φ x , u ,
with ρ : X Ξ and φ : X × R m R m being local diffeomorphisms. To facilitate further investigations, the inverse of the input map (8) is also introduced as:
u = φ 1 x , ν .
From there, system Σ will be considered the nominal dynamics of the Pendubot represented by (1), while system Σ * will be regarded as a transformed version of this dynamic model. In the following subsections, alternative representations of this model will be taken into account. The first representation is based on the decomposition of the mass matrix, while the second comes from the transformation to the so-called normal, which is characterized by a cascade-like structure.

3.2. Transformation Based on Inertial Normalized Quasi-Velocities (NQV)

The concept of inertial normalized quasi-velocities (NQV) proposed by Jain and Rodriguez in the work [34,35] comes from the factorization of the inertia matrix, which is positive definite, non-singular and symmetric. For this purpose, Cholesky factorization can be taken into account. Applying this method with respect to matrix (2) one can write that:
D ( q ) : = L ( q ) L T ( q ) ,
where
L ( q ) = d 2 q 2 a 2 ( 1 + a 32 cos q 2 ) 0 a 2 ,
with a 32 = a 3 a 2 and d 2 q 2 = a 1 a 3 2 a 2 cos 2 q 2 , denotes an upper triangular matrix with positive diagonal entries, [36]. The quasi-velocity NQV can be defined in terms of the matrix L, which contains inertial parameters of the mechanical system, as follows:
σ = L T ( q ) q ˙ .
The linear map (12) dependent on the configuration q can be seen as a part of the state transformation (7). Consequently, assuming that ξ = q T σ T T S 1 × S 1 × R 2 is a new state, and L is invertible for any q, the following global diffeomorphism can be obtained:
ξ = ρ x = q L T ( q ) ω .
In the new states, the dynamics (1) can be represented by (6) with:
f * ξ , ν = L T q 1 σ C σ q , σ σ G σ q + b ν ,
where
C σ ( q , ν ) = L 1 q C ( q , q ˙ ) L ˙ T q L T q 1 = q ˙ 1 + q ˙ 2 a 3 sin q 2 1 d 2 1 + a 32 cos q 2 1 d 2 1 a 2 0 L T q 1 = a 32 d 2 sin q 2 a 2 a 32 cos q 2 σ 1 + d 2 σ 2 0 1 1 0 ,
G σ ( q ) = L 1 ( q ) G ( q ) = 1 d 2 a 4 cos q 1 a 32 a 5 cos ( q 2 ) cos ( q 1 + q 2 ) a 5 a 2 cos ( q 1 + q 2 )
and ν is the transformed input, also known as a quasi-force, according to the general formula (8) with
φ x , u = 1 d 2 q 2 u ,
where q 2 S 1 , d 2 q 2 > 0 .
Remark 1.
The application of NQV makes it possible to transform the nominal system Σ into the form Σ N Q V * , which has a simpler structure. Recalling that b = 1 0 T , it can be concluded from (14) that the input signal ν only affects the variable σ 1 , while the evolution of the variable σ 2 is exclusively due to dynamic couplings and gravity. In addition, taking into account (15), one can state that the description of centrifugal and Coriolis forces is simplified compared to the nominal model.

3.3. Transformation to the Normal Form (NF)

Here, the transformation of the Pendubot based on the methodology investigated in [30,38] is taken into account. At first, the following feedback transformation according to [39] is employed,
φ x , u = ( d 11 + ψ 1 q 2 d 21 ) 1 C 1 * q ˙ q ˙ G 1 ψ 1 q 2 C 2 * q ˙ q ˙ + G 2 + u ,
where
ψ q 2 = a 2 a 2 + a 3 cos q 2 ,
C i * is i th row of matrix C in (3), G i denotes the i th row of vector G in (4) and ν R represents the new input. Using this transformation in the Pendubot dynamics (5), one can obtain:
x ˙ = ω 0 a 2 1 C 2 * ω ω + G 2 + 0 1 ψ 1 q 2 ν .
Next, introducing the new state ξ = θ 1 v 1 θ 2 v 2 T and applying the following state transformation:
ξ = ρ x = q 1 ϑ q 2 ω 1 ψ q 2 ω 2 q 1 ω 1 ,
where
ϑ q 2 = 0 q 2 ψ s ( s ) d s = 2 a 2 a 2 2 a 3 2 arctan a 2 a 3 a 2 + a 3 tan q 2 2 ,
one can obtain the equivalent dynamics with the state function represented as:
f * ξ , ν = v 1 α v 1 2 + β v 1 v 2 + γ v 2 2 + η v 2 0 + 0 0 0 1 ν
with
α ( q 2 ) = a 3 sin q 2 a 2 , β ( q 2 ) = 2 α ( q 2 ) , γ q 2 = a 3 2 sin q 2 cos q 2 a 2 ( a 2 + a 3 cos q 2 ) , η q 1 , q 2 = a 5 cos q 1 + q 2 a 2 + a 3 cos q 2
being scalar functions dependent on the configuration q, which, in view of the transformation (21), can also be expressed in terms of new configuration-like variables θ 1 and θ 2 . Furthermore, since v 1 is a linear combination of ω 1 and ω 2 , it can be considered as a quasi-velocity.
Remark 2.
It should be noted that the dynamic system Σ NF * represented by the state function (23) has a cascade-like structure, for which the control input directly affects only linear subsystem associated to states θ 2 and v 2 and all nonlinearities of the system are represented by the second component of the state function being a polynomial of the second degree with respect to v 1 and v 2 . At the same time, however, it can be seen that a clear interpretation of the state variables of the nominal system is lost as a result of the state transformation.

4. Design of Sub Optimal Stabilizers for the Pendubot

4.1. Equivalence of LQR Design

Since the most important control tool considered in this paper is based on the Linear Quadratic Regulator (LQR) strategy, let us recall the following theorem.
Theorem 1 (LQR feedback, based on [40]).
Consider the following LTI system
z ˙ = A z + B w ,
where z R n , w R m denote the state and the input, while A R n × n and B R m are the state and the input matrices, respectively. Next, consider the following linear feedback:
w = K z ,
where K R m × n denotes the gain matrix, and define the following integral performance index
J = 0 z T ( t ) w T ( t ) W z ( t ) w ( t ) d t ,
where W R ( n + m ) × ( n + m ) is a positive definite symmetric weight matrix, which can be decomposed as follows:
W = Q N N T R 0 ,
while Q R n × n 0 and R R m × m 0 , correspond to state (error) and input weight matrices, and N R n × m defines coupling terms.
Assuming that the feedback (26) makes the system (25) asymptotically stable while minimizing the performance index (27), the optimal gains of the feedback satisfies:
K = R 1 B T S + N T ,
while S is a matrix that solves the Algebraic Ricatti Equation in the form of
A T S + S A + Q S B + N R 1 B T S + N T = 0 .
The LQR can also be used to a nonlinear system Σ taking advantage of its linear approximation. For this purpose, let us assume that x = x 0 is an equilibrium point for system Σ defined at u = u 0 . Then the following approximation can be investigated:
Σ ˜ : x ˜ ˙ = A x ˜ + B u ˜ ,
where x ˜ = x x 0 , u ˜ = u u 0 , A = f x | u = u 0 x = x 0 and B = f u | u = u 0 x = x 0 are Jacobian matrices, while it is assumed that ( A , B ) is the controllable pair. Thus, it is possible to design an LQR-based feedback taking into account the performance index:
J = 0 x ˜ T ( t ) u ˜ T ( t ) W x ˜ ( t ) u ˜ ( t ) d t ,
which is based on (27). Feedback
u ˜ = K x ˜ ,
where K is chosen according to (29) and (30), can also be used for the nonlinear system Σ in a neighborhood of x 0 . In such a case one can consider the following controller
u = K x ˜ + u 0 ,
which can be regarded as a sub-optimal control solution due to the nonlinear nature of the system Σ . Namely, there exists a neighborhood B x 0 of point x 0 such that:
x 0 B x 0 , lim t x ( t ) = x 0 ,
and the criterion function (32) converges to a value close to the minimum established in the case of linear system Σ ¯ .
The stability of the closed-loop system can be conveniently proved using the Lyapunov-like analysis considered below.
Proof of the local exponential stability of the closed-loop system for the controller (34).
Let us consider the nonlinear system Σ . Taking advantage of its linear approximation (31) at x 0 and u 0 , one can write that:
x ˜ ˙ = A x ˜ + B u ˜ + r ˜ 1 ( x ˜ ) + r ˜ 2 ( x ˜ ) u ˜ ,
where r ˜ 1 ( x ˜ ) R n and r ˜ 2 ( x ˜ ) R n describes residual terms. Applying feedback (33) one obtains the following closed-loop dynamics:
x ˜ ˙ = H x ˜ + r ˜ 1 ( x ˜ ) r ˜ 2 ( x ˜ ) K x ˜ ,
with H = A B K . Assume now that the gains K are chosen according to the LQR paradigm which makes matrix H Hurwitz. Thus, the following Lyapunov equation is satisfied: H T P ¯ + P ¯ H = Q ¯ , where P ¯ , Q ¯ R n × n are symmetric positive definite matrices. The Lyapunov function candidate is chosen as:
V = x ˜ T P x ˜ ,
and satisfies the following bounds:
λ min P ¯ x ˜ 2 V λ max P ¯ x ˜ 2 ,
where · defines vector/matrix 2-norm. One can prove that the time derivative of V becomes:
V ˙ = x ˜ T Q ¯ x ˜ + 2 r ˜ 1 ( x ˜ ) r ˜ 2 ( x ˜ ) K x ˜ T P ¯ x ˜ .
To facilitate further analysis one can use the following bounds:
r ˜ 1 ( x ˜ ) r ˜ 2 ( x ˜ ) K x ˜ r ˜ 1 ( x ˜ ) + r ˜ 2 ( x ˜ ) K F x ˜ ,
where · F is the Frobenius matrix norm. Furthermore, for
x ˜ M ¯ ,
where M ¯ is a positive constant, one can consider:
r ˜ 1 ( x ˜ ) C 1 x ˜ 2 , r ˜ 2 ( x ˜ ) K F x ˜ C 2 x ˜ 2 ,
with C 1 , C 2 > 0 being constants. Consequently, the higher order terms can be represented by:
r ˜ 1 ( x ˜ ) + r ˜ 2 ( x ˜ ) K F x ˜ C x ˜ 2 ,
where C = C 1 + C 2 . Using (43) in (40) and recalling that x ˜ T Q ¯ x ˜ λ min Q ¯ x ˜ 2 one can find the following:
V ˙ λ min Q ¯ x ˜ 2 + 2 C x ˜ 2 P ¯ x ˜ = λ min Q ¯ x ˜ 2 1 2 C P ¯ λ min Q ¯ x ˜ .
Recalling (39), one can present (45) as follows:
V ˙ λ min Q ¯ λ max P ¯ V 1 2 C P ¯ λ min Q ¯ λ min P ¯ V .
In order to guarantee the asymptotic stability terms in the bracket in (46) have to satisfy
1 2 C P ¯ λ min Q ¯ λ min P ¯ V γ ,
where γ 0 , 1 is an assumed constant. It can be easily shown that inequality (47) holds if V V ¯ , while V ¯ satisfies
V ¯ = 1 γ 2 λ min Q ¯ 2 λ min P ¯ 4 P 2 C 2 .
Thus, if V ( 0 ) V ¯ , the following conservative bound of V can be considered:
V ˙ γ λ max Q ¯ λ max P ¯ V .
As a result, the closed-loop system (37) is locally exponentially stable. Furthermore, the convergence set can be conservatively estimated by the following set of initial conditions:
x ˜ ( 0 ) M ¯ = V ¯ λ max P ¯ .
 ☐
The LQR control strategy can also be employed for the feedback equivalent system Σ * . Recalling (7) and (8), the equilibrium point of Σ * can be defined by ξ = ξ 0 = ρ x 0 at ν = ν 0 = φ x 0 , u 0 . The linear approximation of Σ * can be represented as:
Σ ˜ * : ξ ˜ ˙ = A * ξ ˜ + B * ν ˜ ,
where ξ ˜ = ξ ξ 0 , ν ˜ = ν ν 0 , and A * = f * ξ | ν = ν 0 ξ = ξ 0 , and B * = f * ν | ν = ν 0 ξ = ξ 0 . Since Σ ¯ is controllable and ρ and φ are diffeomorphisms, system Σ ¯ * is also controllable. Hence, the following state feedback can be designed:
ν ˜ = K * ξ ˜ ,
where K * R m × n is a gain matrix selected to guarantee the Routh-Hurwitz stability of the closed-loop system. This feedback can also be applied to stabilize the nonlinear system Σ * in a neighborhood of ξ 0 . Referring to Formula (34), the following local stabilizer can be considered:
ν = K * ξ ˜ + ν 0 .
Gains K * can be evaluated using the LQR strategy applied with respect to system Σ ¯ * taking into account the redefined performance index:
J = 0 ξ ˜ T ( t ) ν ˜ T ( t ) W * ξ ˜ ( t ) ν ˜ ( t ) d t ,
where W * 0 is some weight matrix.
Here, the question can be raised about the methodology for designing control algorithms for equivalent systems using the LQR approach. In particular, it should be noted that feedback is designed for a specific choice of state and input signal. Thus, the properties of closed-loop systems for the same weight matrices defining the quality index will be different, which can be seen as an undesirable effect. Hence, one can ask how the feedback equivalence property can be used to support control design. To address this issue, one can define the following Taylor expansion of the maps (7) and (8) at x = x 0 and u = u 0 .
ξ ˜ = ρ x | x = x 0 x ˜ + r ξ x ˜ ,
ν ˜ = φ x | u = u 0 x = x 0 x ˜ + φ u | u = u 0 x = x 0 u ˜ + r ν x ˜ , u ˜ ,
while r ξ x ˜ R n and r ν x ˜ , u ˜ R m stand for higher order terms. Then, the following proposition can be considered.
Proposition 1 (Locally equivalent LQR-based feedback design).
Consider equivalent control systems Σ and Σ * and assume that the linear feedback (34) with gain K is designed based on the LQR approach, using the weight matrix described by (28). The control law given by:
u = φ 1 ( x , K * ρ x ρ x 0 + φ ( x 0 , u 0 ) ) ,
where
K * = H u 1 K H x P x
and
P x = ρ x | x = x 0 1 , H x = φ x | u = u 0 x = x 0 , H u = ( φ u | u = u 0 x = x 0 ) 1
are Jacobian matrices, locally equivalent to feedback (34) and provides a comparable performance according to (32) for all x ˜ ( 0 ) < ϵ , while ϵ is set small enough. The matrix parameterizing the performance index (54) is given by:
W * = Q * N * N * T R * ,
where
Q * = P x T Q + H x T H u T R H u H x 2 N H u H x P x , R * = H u T R H u , N * = P x T N H x T H u T R H u .
Proof. 
Taking into account the linear terms of (55), (56) and Jacobian matrices (59) in (52) one obtains:
H x x ˜ + H u 1 u ˜ = K * P x 1 x ˜ .
Computing u ˜ from (62) gives:
u ˜ = H u K * P x 1 + H x x ˜ .
Comparing (63) with (33) one concludes that: K = H u K * P x 1 + H x . Thus, the gain matrix K * satisfies (58).
Next, to find the optimal performance index in new states and inputs, the following inverse transformations based on (55), (56) with (59) can be written:
x ˜ u ˜ = P ξ ˜ ν ˜ , P = P x 0 H u H x P x H u .
Substituting (64) in (32) and computing the following product: P T W P yields the matrix (60) along with (61). ☐
Now it is worth comparing the control law (34), which is designed directly for the nominal system, and the control law (57), for which state feedback is determined in new states and new inputs. To facilitate the description, both structures are shown in Figure 2. It can be clearly seen that while the first stabilizer is fully linear, the second one, in general, is nonlinear.
According to the given proposition, it is possible to guarantee the same properties of the closed-loop system for both algorithms in a sufficiently small vicinity of the desired point x 0 . However, if the basing of attraction is larger for the transformed system Σ * stabilized by a linear feedback, an analogous property will be observed for the nominal system Σ , which is stabilized according to the control law (57). This property indicates that the nonlinear controller (57) can ensure a better performance of the closed-loop system.
Remark 3.
The stability of the closed loop system using the controller (57) can be proved in a similar way that is presented for the controller (34), however, Lyapunov function candidate (38) has to be defined in terms of auxiliary error ξ ˜ . Although the convergence set can be estimated with respect to ξ ˜ , it can be determined in the original space using the inverse of the state transformation (7).
Remark 4.
Analyzing the original Pendubot dynamics (1) and recalling definition of closed-loop dynamics (37) one can see that both nonlinear terms, r 1 and r 2 are present. However, in the case of transformed dynamics QNF and NF described by (14) and (23), respectively, term r 2 in new states ξ ˜ vanishes. This is due to constant input matrix obtained in new representations and it can improve the convergence set in the stabilization task.

4.2. Approximated Models

Stabilization of the Pendubot is considered at the upright position determined by:
q 0 = π 2 0 T , ω 0 = 0 R 2
and u 0 = 0 . Based on this assumption, the transformation matrices (59) were evaluated and collected in Table 1. These show that the transformation to the normal form is more complex, since new coordinate-like states are introduced and the H x term is nonzero, which is due to the presence of an additive term in (18), which is independent of the input u.
For each representation of the Pendubot dynamics, the corresponding linear approximated form can be computed; cf. Table 2. The equations obtained confirm the key role of the gravity component in ensuring the controllability of the linear forms. It is straightforward to show that stabilization under zero gravity by a smooth state feedback would not be possible. It is also easy to notice that linear forms trivially neglect the influence of the resulting centrifugal and Coriolis forces due to the presence of quadratic velocity components, cf. [41].

5. Results

In this section, a comparison of Pendubot stabilization controllers is considered. The research was conducted using both numerical simulation and experimental methods. In order to ensure that the results can be compared, the simulation model takes into account the properties of the laboratory system used in the experiments. Its parameters are summarised in the Table 3. Furthermore, the torque input u was saturated according to the DC motor model. The saturation level and resulting the maximal motor input voltage is equal to 10 V.
The laboratory system shown in Figure 1b is build based on Quanser’s—rotary double inverted pendulum, [42] and consists of the main unit (Rotary Servo Base Unit), which includes the motor, gear with the clearance erasing system and the encoder coupled with the motor, and the passive double pendulum module.

5.1. Simulations

Here three control design approaches are compared. The gains of the linear controller (34) are designed according to the LQR procedure applied to the approximated model of the Pendubot Σ ¯ (cf. Table 2) and taking advantage of the performance index (32) parameterized by the following weight matrices: Q = diag { 50 , 50 , 0.01 , 0.01 } , R = 100 and N = 0 R 4 . The gains of two non-linear controllers described by (57) are established using (58) along with (59) collected in Table 1.
At first attempt the convergence sets where established taking into account Lyapunov analysis considered in the proof of the local exponential stability in Section 4.1. Lyapunov function (38) was chosen for each closed-loop system taking into account the original and transformed errors defined by x ˜ and ξ ˜ , respectively. The quadratic bounds of the residual terms r 1 and r 2 were then numerically approximated in the assumed vicinity of the desired point. In this way, the constant C in (44) can be determined and the upper bound V ¯ can be computed. The set of feasible initial conditions can be found by searching for such states for which V < V ¯ . The sets obtained in the three cases are roughly illustrated in Figure 3. Since the state space is 4-dimensional, the set cannot be visualized on a 2D figure. Therefore, two velocity components q ˙ 1 and q ˙ 2 are replaced by q ˙ presented on the z-axis.
Although the obtained results confirm the local stability of the closed-loop systems, the main task of the research is to compare the attraction basin of each controller in more realistic conditions. To achieve such a comparison, the trajectories of closed-loop systems were evaluated. Such an analysis requires many simulation trails; thus, efficient implementation of simulation models is an important issue. As a result, simulations were carried out with the use of programming tools in the C++ language, including libraries for solving non-stiff differential equations.
For each controller, a discrete set of initial configurations is defined in the form of a two-dimensional grid. Each cell of the grid corresponds to an initial condition represented by ( q 1 ( 0 ) , q 2 ( 0 ) ) and zero velocity ω ( 0 ) . If, for the given condition, the state trajectory converges to the desired point, this trial is considered a positive result and the initial configuration considered is added to the set representing the basin of attraction. The results obtained are presented in Figure 4. Darker cells present an approximation of the basin of attraction, whereas white cells indicate that, for the corresponding initial condition, the control goal has not been accomplished properly. Table 4 presents a comparison of the algorithms considered. The Area [%] index specifies the percentage of positive results related to the entire searched grid.
As part of the extended analysis, the waveforms of the robot configuration and the control input obtained during the simulations for the particular choice of initial configuration are presented along with the experimental results in Section 5.2. To facilitate a comparison between the simulation and the experimental results, the control input is represented by the motor voltage signal instead of the torque u. To further quantify the performance of each algorithm for the chosen initial condition, the index (32) is computed and presented in Table 5.

5.2. Experiment

In the considered application for hardware implementation, a dedicated LabView environment and a driver with an amplifier provided by Quanser are used together with a PC computer, whose task is supervision, monitoring, and measurement registration.
Taking into account the determined basin of attraction for each simulated algorithm, presented in Figure 4, an experimental verification of these algorithms was carried out for a particular selection of initial conditions that belong to the obtained sets. Two different initial postures of the Pendubot were selected. In Scenario 1 the Pendubot initially is tilted to the right, while in Scenario 2 it is tilted to the left. The desired point has been chosen according to (65). For each of the cases considered, a table containing the simulation and experiment conditions, as well as a schematic visualization of the robot initial position. Additionally, the table contains both the controller parameters used during the simulation and the corresponding settings used in the experiment. To make the presentation more clear, the following cases are distinguished:
  • Case A: the linear controller designed for the nominal dynamics Σ is used, see Table 6 and Table 7;
  • Case B: the nonlinear controller designed based on transformed dynamics Σ NQV * is used, see Table 8 and Table 9;
  • Case C: the nonlinear controller designed based on transformed dynamics Σ NF * is used, see Table 10 and Table 11.

6. Discussion

The simulation results obtained show that the LQR design, taking advantage of equivalent systems, can improve the attraction basin in the task of stabilizing the Pendubot. Based on the conservative estimation of attraction sets by the Lyapunov method, the attraction basin for the classical version of LQR is the most limited. This is due to the dependence of the input matrix on the state, which introduces additional nonlinearity, shown in (40). In contrast, for the transformed systems considered Σ * , the input matrix is constant, leading to an increase in the attraction basin. It is worth noting that the largest volume of attraction set is obtained for system Σ NQV * , which is defined in terms of inertial quasi-velocities.
The results of the convergence analysis, based on extensive simulations, see Figure 4, indicate that the considered Lyapunov-based method is too restrictive. Comparing the areas occupied by the cells corresponding to feasible initial configurations presented in Table 4 one can conclude that algorithms based on Formula (57) make it possible to increase the area of acceptable initial configurations from 2.5 to 4.5 times compared to the nominal case. It is interesting that the basin of attraction is the largest for the system Σ NF * while the Lyapunov-based analysis suggests better characteristics with respect to system Σ NQV * .
The results of the simulations performed for the same initial conditions and presented in Figure 5a, Figure 6a, Figure 7a, Figure 8a, Figure 9a, Figure 10a, Figure 11a, Figure 12a, Figure 13a, Figure 14a, Figure 15a, Figure 16a confirm that the step response of closed-loop systems is similar. However, comparing Figure 5a and Figure 7a with Figure 9a, Figure 11a, Figure 13a and Figure 15a more thoroughly, one can state that the transient response obtained for the nominal feedback more clearly exhibits the characteristics inherent in non-linear systems, while for the nonlinear controllers the time plots are smoother and even more characteristic for linear systems. A similar conclusion can be drawn with respect to the analysis of control inputs presented in Figure 6a, Figure 8a, Figure 10a, Figure 12a, Figure 14a and Figure 16a. Similar values of the performance index shown in Table 5 also confirm that the dynamics of the closed-loop system is preserved.
Based on the outcomes presented in Section 5.2, it can be seen that the simulation results and their counterparts obtained in experiments, cf. Figure 5a,b, Figure 6a,b, Figure 7a,b, Figure 8a,b, Figure 9a,b, Figure 10a,b, Figure 11a,b, Figure 12a,b, Figure 13a,b, Figure 14a,b, Figure 15a,b, Figure 16a,b, are fairly similar. Thus, one can cautiously conclude that the mathematical model used to describe the real system is quite close to it; however, some uncertainties can affect the results of the experiment. The similarity manifests itself in terms of the signal amplitudes, but time parameters such as the regulation time is in most cases is longer during experiment. The differences can be explained by the occurrence of effects omitted in the object dynamics model, which include, e.g., static and dynamic friction (occurring in the drive system as well as during the influence of aerodynamic phenomena), backlash and spring effects, and measurement uncertainties. It was observed that the resolutions of the encoders were not enough to obtain high-quality velocity estimates. In particular, the velocity of the second joint cannot be accurately estimated. Taking into account these limitations, adjustments to the gains were required to ensure proper operation of the real system. Especially, in experiments it was necessary to decrease gains associated with velocity components due to the impact of insufficient quality of the velocity estimation. However, decreasing these gains introduces a higher oscillatory response of the closed-loop system, which can be clearly observed in Figure 5, Figure 7 and Figure 13.
In the simulation and experimental tests under consideration, attention should also be paid to the problem of input saturation, which has a significant impact on the attraction set. Despite this, it can be seen that even under input signal constraints, the alternative representation of the dynamics enlarges the convergence set.
There is another issue worth highlighting. Namely, the application of variable transformation may also have the negative effect of increasing the sensitivity of the control system to measurement noise. This may explain the appearance of higher noise in the input signal u for the cases including transformations (cf. Figure 10b, Figure 12b, Figure 14b and Figure 16b) than for the nominal case (cf. Figure 6b and Figure 8b).

7. Conclusions

This paper considers the issue of LQR design for nonlinear systems using a smooth state and input transformation. The proposed design methodology is considered in the Pendubot stabilization task. The properties of the controllers studied were investigated in a simulation environment using experimental tests. Despite some limitations and technical imperfections of the experimental stand, one can conclude that the considered methods to some extent are robust to unmodelled effects and make it possible to provide satisfactory results in real applications.
The results of the tests carried out allow for a hypothesis that the controller using quasi-velocities allows one to increase the range of stabilizer convergence while maintaining the same dynamics of the closed system at the desired point. This property results from the introduction of nonlinearities in the stabilizer equation, which have a positive effect on the properties of the closed-loop system.
To the best of the authors’ knowledge, this work compares for the first time the properties of LQR controllers using different representations of Pendubot dynamics. The detailed forms of transformations and linear approximations given can be regarded as ready-made procedures that can be applied to stabilize similar mechanical systems in robotics.
In the future, the control methodology discussed in this paper can be applied for trajectory tracking, making it possible to also consider swing-up control problems addressed for a class of underactuated systems.

Author Contributions

Conceptualization, D.P.; methodology, D.P. and P.H.; software, P.B. and P.P.; validation, P.B. and P.P.; formal analysis, D.P.; investigation, P.H.; resources, D.P.; data curation, P.P. and D.P.; writing—original draft preparation, D.P., P.P. and P.B.; writing—review and editing, P.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Poznan University of Technology under grant No. 0211/SBAD/0122.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kalman, R. Contributions to the theory of optimal control. Bol. Soc. Mat. Mex. 1960, 5, 102–119. [Google Scholar]
  2. Khalil, H. Nonlinear Systems; Prentice Hall: Upper Saddle River, NJ, USA, 1996. [Google Scholar]
  3. Wang, Y.; Zhang, Z.; Li, C.; Buss, M. Adaptive incremental sliding mode control for a robot manipulator. Mechatronics 2022, 82, 102717. [Google Scholar] [CrossRef]
  4. Cunha, A.; Shao, M.; Huang, Y.; Silberschmidt, V.V. Intelligent Manipulator with Flexible Link and Joint: Modeling and Vibration Control. Shock Vib. 2020, 2020, 4671358. [Google Scholar] [CrossRef] [Green Version]
  5. Batista, J.G.; Souza, D.A.; dos Reis, L.L.; Filgueiras, L.V.; Ramos, K.M.; Junior, A.B.; Correia, W.B. Performance Comparison Between the PID and LQR Controllers Applied to a Robotic Manipulator Joint. In Proceedings of the IECON 2019—45th Annual Conference of the IEEE Industrial Electronics Society, Lisbon, Portugal, 14–17 October 2019; Volume 1, pp. 479–484. [Google Scholar] [CrossRef]
  6. Mason, S.; Righetti, L.; Schaal, S. Full dynamics LQR control of a humanoid robot: An experimental study on balancing and squatting. In Proceedings of the 2014 IEEE-RAS International Conference on Humanoid Robots, Madrid, Spain, 18–20 November 2014; pp. 374–379. [Google Scholar] [CrossRef]
  7. Xin, X.; Tanaka, S.; She, J.; Yamasaki, T. New analytical results of energy-based swing-up control for the Pendubot. Int. J. Non-Linear Mech. 2013, 52, 110–118. [Google Scholar] [CrossRef]
  8. Toan, T.V.; Ha, T.T.; Do, T.V. Hybrid control for swing up and balancing pendubot system: An experimental result. In Proceedings of the 2017 International Conference on System Science and Engineering (ICSSE), Ho Chi Minh City, Vietnam, 21–23 July 2017; pp. 450–453. [Google Scholar] [CrossRef]
  9. Leines, M.T.; Yang, J.S. LQR control of an under actuated planar biped robot. In Proceedings of the 2011 6th IEEE Conference on Industrial Electronics and Applications, Beijing, China, 21–23 June 2011; pp. 1684–1689. [Google Scholar] [CrossRef]
  10. Alcala, E.; Puig, V.; Quevedo, J.; Escobet, T.; Comasolivas, R. Autonomous vehicle control using a kinematic Lyapunov-based technique with LQR-LMI tuning. Control Eng. Pract. 2018, 73, 1–12. [Google Scholar] [CrossRef]
  11. Çimen, T. State-Dependent Riccati Equation (SDRE) Control: A Survey. IFAC Proc. Vol. 2008, 41, 3761–3775. [Google Scholar]
  12. Bernat, J.; Kołota, J.; Stępień, S.; Superczyńska, P. Suboptimal control of nonlinear continuous-time locally positive systems using input-state linearization and SDRE approach. Bull. Pol. Acad. Sci. Tech. Sci. 2018, 66, 17–22. [Google Scholar] [CrossRef]
  13. Stępień, S.; Superczyńska, P. Modified Infinite-Time State-Dependent Riccati Equation Method for Nonlinear Affine Systems: Quadrotor Control. Appl. Sci. 2021, 11, 10714. [Google Scholar] [CrossRef]
  14. Giernacki, W.; Stępień, S.; Chodnicki, M.; Wróblewska, A. Hybrid Quasi-Optimal PID-SDRE Quadrotor Control. Energies 2022, 15, 4312. [Google Scholar] [CrossRef]
  15. Jacobsen, D.; Mayne, D. Differential Dynamic Programming; Elsevier: New York, NY, USA, 1970. [Google Scholar]
  16. Theodorou, E.; Tassa, Y.; Todorov, E. Stochastic Differential Dynamic Programming. In Proceedings of the 2010 American Control Conference, Baltimore, MD, USA, 30 June–2 July 2010; pp. 1125–1132. [Google Scholar] [CrossRef] [Green Version]
  17. Li, W.; Todorov, E. Iterative Linear Quadratic Regulator Design for Nonlinear Biological Movement Systems. In Proceedings of the the First International Conference on Informatics in Control, Automation and Robotics (ICINCO), Setúbal, Portugal, 25–28 August 2004; Araújo, H., Vieira, A., Braz, J., Encarnação, B., Carvalho, M., Eds.; pp. 222–229. [Google Scholar]
  18. Van den Berg, J. Iterated LQR smoothing for locally-optimal feedback control of systems with non-linear dynamics and non-quadratic cost. In Proceedings of the 2014 American Control Conference, Portland, OR, USA, 4–6 June 2014; pp. 1912–1918. [Google Scholar] [CrossRef]
  19. Respondek, W. Feedback classification of nonlinear control systems on R2 and R3. In Geometry of Feedback and Optimal Control; Jakubczyk, B., Respondek, W., Eds.; Dekker: New York, NY, USA, 1998; pp. 347–381. [Google Scholar]
  20. Jakubczyk, B.; Respondek, W. On linearization of control systems. Bull. Acad. Polon. Sci. Ser. Sci. Math. 1980, 28, 517–522. [Google Scholar]
  21. Hammami, M.; Rettab, N. On the region of attraction of dynamical systems: Application to Lorenz equations. Arch. Control Sci. 2020, 30, 389–409. [Google Scholar] [CrossRef]
  22. Zhao, K. Local exponential stability of four almost-periodic positive solutions for a classic Ayala-Gilpin competitive ecosystem provided with varying-lags and control terms. Int. J. Control 2022, 1, 1–13. [Google Scholar] [CrossRef]
  23. Valmorbida, G.; Anderson, J. Region of attraction estimation using invariant sets and rational Lyapunov functions. Automatica 2017, 75, 37–45. [Google Scholar] [CrossRef] [Green Version]
  24. El-Guindy, A.; Han, D.; Althoff, M. Estimating the region of attraction via forward reachable sets. In Proceedings of the 2017 American Control Conference (ACC), Seattle, WA, USA, 24–26 May 2017; pp. 1263–1270. [Google Scholar] [CrossRef] [Green Version]
  25. Ziętkiewicz, J. Linear quadratic control with feedback-linearized models. Stud. Autom. Inf. Technol. 2015, 40, 37–49. [Google Scholar]
  26. Owczarkowski, A.; Horla, D.; Zietkiewicz, J. Introduction of Feedback Linearization to Robust LQR and LQI Control—Analysis of Results from an Unmanned Bicycle Robot with Reaction Wheel. Asian J. Control 2019, 21, 1028–1040. [Google Scholar] [CrossRef]
  27. Spong, M.; Block, D. The Pendubot: A mechatronic system for control research and education. In Proceedings of the 1995 34th IEEE Conference on Decision and Control, New Orleans, LA, USA, 13–15 December 1995; Volume 1, pp. 555–556. [Google Scholar] [CrossRef]
  28. Prasad, L.B.; Tyagi, B.; Gupta, H.O. Optimal Control of Nonlinear Inverted Pendulum System Using PID Controller and LQR: Performance Analysis Without and With Disturbance Input. Int. J. Autom. Comput. 2014, 11, 661–670. [Google Scholar] [CrossRef] [Green Version]
  29. Cimborová, K.; Jadlovská, S. Modeling of benchmark underactuated systems via different approaches. IFAC-PapersOnLine 2020, 53, 8935–8940. [Google Scholar] [CrossRef]
  30. Li, S.; Moog, C.; Respondek, W. Maximal feedback linearization and its internal dynamics with applications to mechanical systems on R4. Int. J. Robust Nonlinear Control 2019, 29, 2639–2659. [Google Scholar] [CrossRef]
  31. Parulski, P.; Bartkowiak, P.; Pazderski, D. Evaluation of Linearization Methods for Control of the Pendubot. Appl. Sci. 2021, 11, 7615. [Google Scholar] [CrossRef]
  32. Meirovitch, L. Methods of Analytical Dynamics; McGraw–Hill: New York, NY, USA, 1970. [Google Scholar]
  33. Gutowski, R. Mechanika Analityczna; PWN: Warsaw, Poland, 1971. [Google Scholar]
  34. Jain, A.; Rodriguez, G. Diagonalized Lagrangian robot dynamics. IEEE Trans. Robot. Autom. 1995, 11, 571–584. [Google Scholar] [CrossRef]
  35. Rodriguez, G.; Jain, A.; Kreutz-Delgado, K. A Spatial Operator Algebra for Manipulator Modeling and Control. Int. J. Robot. Res. 1991, 10, 371–381. [Google Scholar] [CrossRef] [Green Version]
  36. Herman, P. A controller of the pendubot using quasi-velocities. In Proceedings of the 2008 16th MED Conference, Vancouver, BC, Canada, 27–31 October 2008; pp. 1066–1070. [Google Scholar] [CrossRef]
  37. Olfati-Saber, R. Nonlinear Control of Underactuated Mechanical Systems with Application to Robotics and Aerospace Vehicles. Ph.D. Thesis, Massachusetts Institute of Technology, Cambridge, MA, USA, 2001. AAI0803036. [Google Scholar]
  38. Olfati-Saber, R. Normal forms for underactuated mechanical systems with symmetry. IEEE Trans. Autom. Control 2002, 47, 305–308. [Google Scholar] [CrossRef] [Green Version]
  39. Spong, M.W. Partial feedback linearization of underactuated mechanical systems. In Proceedings of the IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS’94), Munich, Germany, 12–16 September 1994; Volume 1, pp. 314–321. [Google Scholar]
  40. Kwakernaak, H.; Sivan, R. Linear Optimal Control System; Wiley-Interscience: New York, NY, USA, 1972. [Google Scholar]
  41. Nowicki, M.; Respondek, W. A Mechanical Feedback Classification of Linear Mechanical Control Systems. Appl. Sci. 2021, 11, 10669. [Google Scholar] [CrossRef]
  42. Quanser. Rotary Double Inverted Pendulum. Available online: www.quanser.com/products/rotary-double-inverted-pendulum/ (accessed on 12 July 2022).
Figure 1. Scheme of the Pendubot: (a) basic model: m i —mass of i th link, I i —moment of inertia of i th link determined with respect to the axis normal to the plane and moving through the center of mass of the link, l i —length of i th link, l c i —distance between the revolute joint and the center of mass of i th link, g—gravity acceleration, (b) the real system used in experimental research.
Figure 1. Scheme of the Pendubot: (a) basic model: m i —mass of i th link, I i —moment of inertia of i th link determined with respect to the axis normal to the plane and moving through the center of mass of the link, l i —length of i th link, l c i —distance between the revolute joint and the center of mass of i th link, g—gravity acceleration, (b) the real system used in experimental research.
Energies 15 05146 g001
Figure 2. Locally equivalent control structures (linear parts are denoted by dotted lines): (a) Classical linear controller, (b) Nonlinear controller using linear feedback in new set of coordinates.
Figure 2. Locally equivalent control structures (linear parts are denoted by dotted lines): (a) Classical linear controller, (b) Nonlinear controller using linear feedback in new set of coordinates.
Energies 15 05146 g002
Figure 3. Area of convergence estimated in a conservative way based on the Lyapunov analysis: (a) Σ , (b) Σ NQV * , (c) Σ NF * .
Figure 3. Area of convergence estimated in a conservative way based on the Lyapunov analysis: (a) Σ , (b) Σ NQV * , (c) Σ NF * .
Energies 15 05146 g003
Figure 4. Area of convergence: (a) Σ , (b) Σ NQV * , (c) Σ NF * .
Figure 4. Area of convergence: (a) Σ , (b) Σ NQV * , (c) Σ NF * .
Energies 15 05146 g004
Figure 5. Case A1: Angular positions: (a) simulation, (b) experiment.
Figure 5. Case A1: Angular positions: (a) simulation, (b) experiment.
Energies 15 05146 g005
Figure 6. Case A1: Input signal: (a) simulation, (b) experiment.
Figure 6. Case A1: Input signal: (a) simulation, (b) experiment.
Energies 15 05146 g006
Figure 7. Case A2: Angular positions: (a) simulation, (b) experiment.
Figure 7. Case A2: Angular positions: (a) simulation, (b) experiment.
Energies 15 05146 g007
Figure 8. Case A2: Input signal: (a) simulation, (b) experiment.
Figure 8. Case A2: Input signal: (a) simulation, (b) experiment.
Energies 15 05146 g008
Figure 9. Case B1: Angular positions: (a) simulation, (b) experiment.
Figure 9. Case B1: Angular positions: (a) simulation, (b) experiment.
Energies 15 05146 g009
Figure 10. Case B1: Input signal: (a) simulation, (b) experiment.
Figure 10. Case B1: Input signal: (a) simulation, (b) experiment.
Energies 15 05146 g010
Figure 11. Case B2: Angular positions: (a) simulation, (b) experiment.
Figure 11. Case B2: Angular positions: (a) simulation, (b) experiment.
Energies 15 05146 g011
Figure 12. Case B2: Input signal: (a) simulation, (b) experiment.
Figure 12. Case B2: Input signal: (a) simulation, (b) experiment.
Energies 15 05146 g012
Figure 13. Case C1: Angular positions: (a) simulation, (b) experiment.
Figure 13. Case C1: Angular positions: (a) simulation, (b) experiment.
Energies 15 05146 g013
Figure 14. Case C1: Input signal: (a) simulation, (b) experiment.
Figure 14. Case C1: Input signal: (a) simulation, (b) experiment.
Energies 15 05146 g014
Figure 15. Case C2: Angular positions: (a) simulation, (b) experiment.
Figure 15. Case C2: Angular positions: (a) simulation, (b) experiment.
Energies 15 05146 g015
Figure 16. Case C2: Input signal: (a) simulation, (b) experiment.
Figure 16. Case C2: Input signal: (a) simulation, (b) experiment.
Energies 15 05146 g016
Table 1. Local transformation matrices: Σ NQV * − Equation (14), Σ NF * − Equation (23).
Table 1. Local transformation matrices: Σ NQV * − Equation (14), Σ NF * − Equation (23).
SystemTransformations
P x H x H u
Σ NQV *   I 2 × 2 0 2 × 2 0 2 × 2 L T q 0 1    0 1 × 4    a 1 a 2 a 3 2 a 2  
Σ NF *    0 0 1 0 1 + a 32 0 1 a 32 0 0 0 0 1 0 1 + a 32 0 1 a 32    a 4 a 2 a 5 a 3 a 1 a 2 a 3 2 a 5 a 3 a 1 a 2 a 3 2    a 1 a 2 a 3 2 a 2  
Table 2. Approximated models of nonlinear equivalent systems describing Pendubot dynamics: Σ − Equation (1), Σ NQV * − Equation (14), Σ NF * − Equation (23).
Table 2. Approximated models of nonlinear equivalent systems describing Pendubot dynamics: Σ − Equation (1), Σ NQV * − Equation (14), Σ NF * − Equation (23).
Nonlinear ModelLinear Model Σ ¯
DriftInput
Σ A = 0 2 × 2 I 2 × 2 D 1 q G q q 0 2 × 2 q = q 0 B = 0 2 × 1 D 1 q 0 b
Σ NQV * A * = 0 2 × 2 L T q 1 G σ q q 0 2 × 2 q = q 0 B * = 0 2 × 1 b
Σ NF * A * = 0 1 0 0 η θ 1 0 η θ 2 0 0 0 0 1 0 0 0 0 θ = θ 0 B * = 0 3 × 1 1
Table 3. Pendubot robot parameters.
Table 3. Pendubot robot parameters.
LinkMassLink LengthMass CenterInertia
i m i [ kg ] l i [ m ] l c i [ m ] I i [ kg m 2 ]
10.0970.200.16350.0069 *
20.1270.33650.17780.0048
* Inertia I1 takes into account the weight of the encoder menc = 0.141 kg.
Table 4. Comparison of areas occupied by the basin of attraction sets estimated based on numerical simulations.
Table 4. Comparison of areas occupied by the basin of attraction sets estimated based on numerical simulations.
     Algorithm          Area [%]     
Σ , (34)1.10
Σ NQV * , (57)2.91
Σ NF * , (57)4.98
Table 5. Comparison of values of performance index obtained for the same initial configuration.
Table 5. Comparison of values of performance index obtained for the same initial configuration.
     Algorithm          J     
Σ , (34)10.96
Σ NQV * , (57)10.45
Σ NF * , (57)11.21
Table 6. Case A1: Initial posture of the Pendubot and parameters of the linear controller.
Table 6. Case A1: Initial posture of the Pendubot and parameters of the linear controller.
Energies 15 05146 i001
Type Σ , (34)
( q 1 ( 0 ) , q 2 ( 0 ) ) ( 65 , 25 ) deg
SIM:K { 10.4 , 9.7 , 2.5 , 1.9 }
EXP:K { 9.8 , 9 , 0.6 , 0.3 }
Table 7. Case A2: Initial posture of the Pendubot and parameters of the linear controller.
Table 7. Case A2: Initial posture of the Pendubot and parameters of the linear controller.
Energies 15 05146 i002
Type Σ , (34)
( q 1 ( 0 ) , q 2 ( 0 ) ) ( 130 , 40 ) deg
SIM:K { 10.4 , 9.7 , 2.5 , 1.9 }
EXP:K { 10.4 , 9.7 , 2.5 , 1.9 }
Table 8. Case B1: Initial posture of the Pendubot and parameters of the nonlinear controller.
Table 8. Case B1: Initial posture of the Pendubot and parameters of the nonlinear controller.
Energies 15 05146 i003
Type Σ NQV * , (57)
( q 1 ( 0 ) , q 2 ( 0 ) ) ( 65 , 25 ) deg
SIM: K * { 94 , 88 , 26 , 184 }
EXP: K * { 150 , 150 , 11 , 100 }
Table 9. Case B2: Initial posture of the Pendubot and parameters of the nonlinear controller.
Table 9. Case B2: Initial posture of the Pendubot and parameters of the nonlinear controller.
Energies 15 05146 i004
Type Σ NQV * , (57)
( q 1 ( 0 ) , q 2 ( 0 ) ) ( 130 , 40 ) deg
SIM: K * { 94 , 88 , 26 , 184 }
EXP: K * { 150 , 150 , 11 , 100 }
Table 10. Case C1: Initial posture of the Pendubot and parameters of the nonlinear controller.
Table 10. Case C1: Initial posture of the Pendubot and parameters of the nonlinear controller.
Energies 15 05146 i005
Type Σ NF * , (57)
( q 1 ( 0 ) , q 2 ( 0 ) ) ( 65 , 25 ) deg
SIM: K * { 1187 , 236 , 314 , 26 }
EXP: K * { 1187 , 130 , 314 , 8 }
Table 11. Case C2: Initial posture of the Pendubot and parameters of the nonlinear controller.
Table 11. Case C2: Initial posture of the Pendubot and parameters of the nonlinear controller.
Energies 15 05146 i006
Type Σ NF * , (57)
( q 1 ( 0 ) , q 2 ( 0 ) ) ( 130 , 40 ) deg
SIM: K * { 1187 , 236 , 314 , 26 }
EXP: K * { 1187 , 130 , 314 , 8 }
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pazderski, D.; Parulski, P.; Bartkowiak, P.; Herman, P. Sub-Optimal Stabilizers of the Pendubot Using Various State Space Representations. Energies 2022, 15, 5146. https://doi.org/10.3390/en15145146

AMA Style

Pazderski D, Parulski P, Bartkowiak P, Herman P. Sub-Optimal Stabilizers of the Pendubot Using Various State Space Representations. Energies. 2022; 15(14):5146. https://doi.org/10.3390/en15145146

Chicago/Turabian Style

Pazderski, Dariusz, Paweł Parulski, Patryk Bartkowiak, and Przemysław Herman. 2022. "Sub-Optimal Stabilizers of the Pendubot Using Various State Space Representations" Energies 15, no. 14: 5146. https://doi.org/10.3390/en15145146

APA Style

Pazderski, D., Parulski, P., Bartkowiak, P., & Herman, P. (2022). Sub-Optimal Stabilizers of the Pendubot Using Various State Space Representations. Energies, 15(14), 5146. https://doi.org/10.3390/en15145146

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