Abstract

The design and tuning of a simple feedback strategy with delay to stabilize a class of underactuated mechanical systems with dead time are presented. A linear time-invariant (LTI) model with time delay of fourth order and a Proportional Retarded (PR) controller are considered. The PR controller is shown as an appealing alternative to the application of observer-based controllers. This paper gives a step forward to obtain a better understanding of the effect of output delays and related phenomena in mechatronic systems, making it possible to design resilient control laws under the presence of uncertain time delays in measurements and obtain an acceptable performance without using a derivative action. The Furuta pendulum is a standard two-degrees-of-freedom benchmark example from the class of underactuated mechanical systems. The configuration under study includes an inherent output delay due to wireless communication used to transmit measurements of the pendulum’s angular position. Our approach offers a constructive design and a procedure based on a combination of root loci and Mikhailov methods for the analysis of stability. Experiments over a laboratory platform are reported and a comparison with a standard linear state feedback control law shows the advantages of the proposed scheme.

1. Introduction

Underactuated mechatronic systems have been a subject of intensive research during the past three decades, where partial feedback linearization, normal forms, and energy-based methods have offered successful results as well as breakthrough paradigms [19]. Experiments with pendulums having one or more degrees of underactuation play an important role in control theory. One important challenge concerning the construction of such equipment is how to send the measurements from the encoders located in the underactuated links without adding restrictions to the movement; one possible solution is wireless communication [10, 11]. Nowadays, wireless sensor networks are a feasible option when mobility and flexibility are requirements. However, when the sensor signals are transmitted through a wireless network, the time-delay phenomena are unavoidable, representing an important as well as an open challenge [1113]. Time delays in measurements are often neglected. Moreover, when only positions are instrumented, the corresponding velocities are often required to be reconstructed, either with an observer or with the use of on-line differentiators or filters (see, e.g., [14, 15]). It is well known that noise, time delays, sampling, and unmodeled dynamics deteriorate the performance of such estimators. Hence, it is of relevance to understand the effect of time delays in measurements as well as the potential of such delays as a second degree of freedom for control design (see [16, 17]). For this aim, a PR controller design will be presented together with the Mikhailov stability approach. Furthermore, a Furuta pendulum case study, where the underactuated variable is subject to an inherent output delay due to wireless signal transmission, is presented. The analysis of time delays in control systems is not a new topic; particularly, the addition of time delays in the control design has been studied for many years (see, e.g., [1827]). In particular, the introduction of time delays in the control law avoids the use of techniques for estimating the derivative action, which is not an easy task under the presence of noise and delays in measurements and increases the cost and complexity of the overall controller design and tuning of the gains [16]. On the other hand, the implementation of a retarded control law requires only memory registers and is relatively simple to program, since it is based on a sum of a proportional signal and the corresponding delayed one. However, the stability analysis of a closed-loop system with a retarded feedback control law requires dealing with infinite-dimensional systems even if the initial plant model is not. This analysis is studied in the framework of two approaches: time domain and frequency domain [28, 29]. In the first approach, the analysis is based on the well-known Lyapunov’s criteria using linear matrix inequalities (LMIs) via convex optimization. However, this approach provides only sufficient conditions of stability or stabilization, which are usually conservative for some systems, while for other systems the LMIs are unfeasible. Besides, in our experience for similar problems, the parameters tuning by a frequency method using LMIs sometimes is not feasible. In the second approach, the analysis is based on the root locus of its corresponding characteristic equation in the complex plane. Unlike the temporal approach, the frequency approach can provide necessary and sufficient conditions that are not conservative. In [3032], some numerical methods to approximate the roots of these quasi-polynomial or analytic functions are displayed. Furuta pendulum is a well-known academic example that is used as a benchmark for the study of underactuated mechanical systems [33, 34]. Recently, in [17], the control design of a PR controller for the stabilization of Furuta pendulum was introduced; however, the effect of the inherent delay, present in the unactuated variable, was not studied. In contrast to the previous work [17], this paper presents a more detailed analysis for the tuning of control gains; also new experiments and graphs are presented, as well as experimental tests of robustness under external disturbances. In addition, the organization of the document is better and more natural; also a more comprehensive bibliographic review is performed and the introduction presents a better justification of the topic of study. Besides, to the best of our knowledge, the stability analysis including the inherent time delay as well as the enhancement of the PR control parameters has not been presented for Furuta pendulum. In addition, the Mikhailov procedure [35] is described and formalized for the present case study; concerning application of the method to systems with delay, see, for example, [36, 37]. The proposed control design is composed of two elements: a PD controller for the actuated variable and PR controller for the stabilization of the unactuated link. The main difference with respect to other works is that the stability analysis incudes two delays, one is inherent of the wireless communication and the other is introduced in the control law for stabilization purposes. Moreover, the time derivative of the unactuated variable or its on-line estimation is not required in the proposed control law. Thus, the characteristic equation analyzed is an analytical function or quasi-polynomial with two transcendental terms of the formwhere is the inherent delay of the angular position measurements due to the wireless communication, (), , are polynomials with real coefficients, , , τ are the parameters of PR controller, and , are the parameter of PD controller. The first stability analysis is using the D-partition method [38] and the root continuity property [39, 40]. A growing body of literature has analyzed the stability of quasi polynomials (see [23, 26, 4144]). Particularly in [26], the delay phenomena have been studied for a fully actuated cart-pendulum system, introducing a full-state feedback controller with feedforward term, assuming that all velocities and positions are measurable.

The present paper focuses on the delay phenomena in an underactuated system called Furuta pendulum. The main result introduces a methodology for the optimal selection of the controller gains and a procedure based on Mikhailov method for the study of complete stability under the presence of measurement delays. Here, it is assumed that only angular positions are available for measurement. For comparison purposes, a design of a state feedback controller, which includes the on-line estimation of velocities, is presented. All the proposed algorithms are verified with simulations and experiments over a laboratory test bench, obtaining the desired performance.

The outline of this paper is distributed as follows. In Section 2, the mathematical model of the Furuta pendulum is introduced. The proposed control design is presented in Section 3. Experimental results are found in Section 4. The conclusions of this study are drawn in Section 5.

2. Mathematical Model

The Furuta pendulum test bench consists of a passive pendulum attached to the end of an actuated rotating arm as shown in Figure 1. The coordinates of the centers of mass of both links are given by = and = , where is the angular position of the pendulum and is the angle of the rotational arm; and both variables represent the generalized coordinates that describe the motion. The following notation is used: is the mass of the arm, is the pendulum mass, and are the moments of inertia of the arm and of the pendulum, respectively, and are the lengths of the arm and the pendulum, respectively, and is the acceleration due to gravity. Note that with the chosen coordinated system the unstable equilibrium point is located at = . Note that the equilibrium does not depend on ϕ and every real value of corresponds to an unstable equilibrium. A model of the Furuta pendulum can be obtained applying the Euler-Lagrange formalism (see, e.g., [4547]). The corresponding Lagrangian is given by , with + + and leading to the following expression:where the related terms are , , , , and The elements of the inertia matrix are , , and ; the Coriolis matrix is formed by , , and ; the parameters of viscous friction are included in matrix . For the particular laboratory test bench in Figure 1, static friction does not play an important role, because it is small and bounded, and it is neglected in the plant model.

The external torque that allows the control of the arm is produced by a DC motor. The motor induces a torque , where is the normalized input voltage of the PWM amplifier, is constant, and .

There are two control problems commonly associated with the Furuta pendulum: stabilization around one of the two equilibrium points and the swing-up [48, 49]. In order to study the inherent delay due to the wireless communication together with the design of a retarded control law, here we are going to focus on the problem of stabilization.

In this context, the equation of motion (2) can be linearized around the unstable equilibrium point = , obtaining the following controllable linear system: where is the available output subject to communication delay The total delay involved in the sensor communication can be divided into three parts: the communication delay from sensor to controller, the computation delay at the controller, and the communication delay from controller to actuator. Here, any delay from the sensor to the controller is assumed to be small relative to the sampling period, since short delays are typical for optical encoders. Also, the computation delay is relatively small compared with the communication delay that normally is affected by the physical distance between the communication nodes. As a result, is modeled as a static delay within a bounded set; that is, . The other system parameters are , , , , , , , and . Here, the physical parameters of the experimental platform given in Table 1 are employed.

3. Controller Design and Tuning

Using classical techniques, stabilization of the Furuta pendulum requires measurements of variables θ and ϕ and the estimation of the corresponding time derivatives and (see, e.g., [33, 48, 49]). However, when wireless communication is considered, delays in the measurement deteriorate the accuracy of estimation of time derivatives. The system under study includes two encoders for measurements of and , where the communication for is via a wireless link. Since the wireless communication introduces a variable time delay in the measurement, the estimation of time derivative for variable is not a trivial task and we aim to avoid this by using a PR control action. The control methodology includes the design of a PR controller to stabilize the passive degree of freedom associated with together with a low-gain PD controller for variable . The D-partition method, root locus, and method of Mikhailov are implemented for gain design and stability analysis. Some preliminaries of Mikhailov method and second-order PR controllers are given in Appendices A and B, respectively.

3.1. Control Scheme

Consider the following control law: where ; represents the inherent delay of the angular position measurements due to the wireless communication, and , and τ are parameters of design. For application of this control law, a method for derivative estimation is necessary for variable ; in this case, either a high-gain observer or a super twisting differentiator can been used in order to obtain (see [15, 50]). Thus, the closed-loop system, (3)–(5), is and its representation in state space is given as follows:whereThus, the characteristic quasi polynomial of system (7) isNote that, for and , it reduces to Clearly, the first factor of the above quasi polynomial is identical to the one for a decoupled one-degree-of-freedom mechanical system under a PR feedback (see (A.3) from Appendix A). Then, as a first step, let us assume small values for , and , proceeding with the design of , and with the help of Lemmas A.1 and A.2 given in Appendix B, placing some of the roots into the optimal locations for the reduced system. Lemmas A.1 and A.2 bring -stability to the dynamics of the variable . Then, choosing the right values for and , one should not expect -stability for the full fourth-order system but may expect asymptotic stability under sufficiently small and appropriately chosen and , which appears only in the terms associated with in (9). This is expected due to continuous dependence of the roots of a quasi polynomial on the coefficients and will be quantified in the next section with the analysis of root locus.

As a second step, the Mikhailov analysis will provide the stability conditions under the presence of an inherent sufficiently small delay

3.2. Tuning and Stability Analysis

In this section, the tuning of the controller gains proposed in the previous section and an analysis of stability of closed-loop system (7) with are presented.

Consider the second-order quasi polynomial, roots of which cover all the roots of (10) except for and : where , and are tuned for the reduced system (see Lemma A.2). This choice provides three dominant roots in the same point at , as can be seen in Figure 2. Suppose that α and β are small and ; then quasi polynomial (9) has five dominant roots, three close to and two close to zero, as can be seen in Figure 3. Therefore, system (7) is close to be marginally stable.

To determine stability of quasi polynomial (9), in Figure 4, a root locus is presented. Here, α and β are varied. Observe that if = 0 and is varied or and is varied, then quasi polynomial (9) remains unstable. The case where ≠ 0 and ≠ 0 can be handled with the same method; however, for some combination of values, it can be exhaustive to find the stability condition or quasi polynomial (9) can turn unstable. Meanwhile, if the value of = is increased simultaneously, then quasi polynomial (9) can be stabilized relatively easier.

Next, the root locus of quasi polynomial (9) is presented when α = β > 0 are simultaneously varied as can be seen in Figure 5. Here, two of the three dominant roots in are moved to the right and one root is moved to the left, while the two dominant roots at zero and at are moved to the left. Clearly, there is a breaking point such that quasi polynomial (9) becomes unstable; this occurs when Thus, the stability of system (3) in closed loop with control law (5) is ensured, when . The above analysis enables a proper selection of the gains α and β such that quasi polynomial (9) is stable. Without loss of generality, the gains have been set at α = β = 0.5. Here, the five dominant roots are , , and (see roots marked with blue asterisk in Figure 5). Note that if , then quasi polynomial (9) is unstable.

3.3. Robustness Analysis under the Presence of

Now, Mikhailov’s criterion [35] to corroborate the stability under the presence of and is implemented (see Appendix B). For this purpose, consider the next closed curve in the complex plane: where is the interval on the imaginary axis and is a semicircle of radius in the right-half complex plane, where is sufficiently large (we take ). Mikhailov’s criterion is used to analyze robustness when . In Figures 6(a) and 6(b), one clockwise turn and one counterclockwise turn around zero are depicted. In Figures 6(c), 6(d), and 6(e), two clockwise turns and two counterclockwise turns around zero are observed. It should be mentioned that such behavior is verified for . Nevertheless, the figures have been omitted to save the space. Thus, due to continuous dependence of the roots on the parameter and due to the above observations, it can be concluded that the change of argument of quasi polynomial (9) is when . However, in Figure 6(f), two clockwise turns and one counterclockwise turn around zero are observed when . Here, . Therefore, quasi polynomial (9) remains stable if .

On the other hand, in Figure 7, the root locus of quasi polynomial (9) when and is shown. Observe that (9) remains stable if and ; and it is unstable if and .

Therefore, we have the following result.

Proposition 1. Consider a quasi polynomial of the form (9) and a positive constant . Then, the control law (5) stabilizes a fourth-order system of the form (6) ifand , and are such that and . Here, follows from (A.5) and is given in Section 3.2; .

4. Experiments

The experimental results were obtained in the Robotics and Control Lab at the Department of Applied Physics and Electronics, Umeå University, Umeå, Sweden. The system used for the experiments is from the PendCon Advanced series [10] with the configuration for rotary pendulum or Furuta pendulum (see Figure 1). In this hardware setup, signals from the encoders with resolution of 4096 pulses per revolution are transmitted wirelessly, and for the actuated link the motor MAXON RE40 (40 mm) is used, which is a graphite brushless DC motor (150 W) with nominal voltage of 24 V. The real-time controller is on a dedicated target PC that communicates with the host PC; it has one PCI slot for the encoder card and four interfaces for RS-232 signals; this PC has Intel Core i5 processor with real-time workshop; the host-target environment is a PC with MATLAB real-time workshop compiler installed on Intel Core i5 processor with Window 7. Both communicate with the Conversion Unit via a Matlab-compatible Quatech RS-232 card. The control system is developed in MATLAB R2015a/Simulink and downloaded to the target PC.

The control goal is to stabilize the pendulum at its unstable equilibrium point , that is, vertically pointing up. For comparison purposes, a state feedback controller is designed: with . In order to have a reasonable comparison, the gain is designed such that the poles in closed loop of the linearized system are located at and . Note that these values correspond with the dominant complex conjugate poles obtained with the PR controller. In Figure 8, the performance of regulation of the angular positions θ() and () with the PR and state feedback controller are presented. Here, the initial position is close to the unstable equilibrium point . Figure 9 presents the control signals.

The design and tuning of the controller are done as follows:(1)Consider a system of the form (3)-(4) in closed loop with a control law of the form (5) and its corresponding quasi polynomial of 4th order , specified in (9).(2)Factor the quasi polynomial of 4th order in one quasi polynomial of 2nd order and one polynomial of 2nd order, when and (see (10)):(3)Use Lemma A.2 to -stabilize the quasi polynomial of 2nd order and obtain the parameters of the PR controller (see (11)):(4)Obtain α, β, and such that quasi polynomial (9) remains stable, as depicted in Figures 6 and 7.

Remark 2. It is important to remark that there is an inherent unknown delay , in the measurement of pendulum position; that is, at time , we obtain . However, if remains inside a small range, stability is still preserved as we have shown.

In order to compute the estimation of the first derivative for both position variables, θ and , two methods were explored: high gain and super twisting differentiators (see [15, 50]).

These methods are designed to be used on-line. For comparison purposes, it is obvious that better velocity estimation can be achieved with an off-line method when both previous and future values of the position are used. Then, we propose to postprocess the measured position to obtain an off-line estimation. For this purpose, the measured signal is fitted with a smoothing spline and next the off-line estimation of velocity is obtained as an analytical differentiation of the spline (see, e.g., [51], page 194). Using this off-line estimation as the true value of the velocity, we compare the performances of high gain (with different gains: HG1 and HG2) and super twisting (ST) differentiators, computing the estimation errors as is shown in Table 2. From Table 2, it is clear that, due to the wireless communication, the estimation of is deteriorated, resulting in a bigger error and decreased closed-loop performance when we use this estimation with a state feedback controller.

Remark 3. With the proposed PR controller, the use of the estimation of is not required. This is desirable in applications where the estimation of velocity deteriorates the overall performance.

An additional experiment was performed in order to test the robustness against external perturbations. For this aim, several impacts were applied to the pendulum at the time instants , as can be seen in Figures 10 and 11.

From the experimental results, the recovery of performance with the proposed PR controller is verified. Besides, it is important to highlight that a better transient, a smoother response, and chattering attenuation in the control signal are observed with the proposed approach under the presence of delays in measurements. A video of the recorded experiments is available in the following link: https://www.youtube.com/watch?v=F9QbPVkagwE.

5. Conclusion

Presence of time delays, quantization errors of encoders, and errors of the wireless communication make the design of a robust control law difficult. The design and stability analysis of a PR action-based controller for stabilization of the Furuta pendulum under the presence of time delays in measurements is introduced in this work. With the proposed approach, the design of observers for the passive variable is avoided simplifying the control algorithm and recovering an acceptable performance without relying on accurate estimates of the derivative for feedback. Since a retarded action is included in the control law, the closed-loop dynamics are represented by a characteristic quasi polynomial and the stability analysis was successfully done with Mikhailov’s criterion and a root locus of the corresponding quasi polynomial. Furthermore, the stability analysis includes an inherent time delay due to the wireless communication. Besides, a two-step constructive procedure for tuning the four parameters of a controller is presented. Experimental results and a comparison with a state feedback controller show the advantages of the proposed methodology. We are aware that a possible extension of the proposed controller is the addition of a PR controller in the actuated variable as well; however, this is considered as a future work, since the order of the quasi polynomial will increase significantly, losing the scope of the actual contribution. Additionally, extension for a wider class of underactuated mechanical systems and the inclusion of nonlinearities is considered for future work.

Appendix

A. PR Controller for a Second-Order Plant

Consider the next second-order system, and the PR control law, where the design parameters are as follows: is the proportional gain, is the retarded gain, and τ is the time delay. Thus, the characteristic quasi polynomial of closed-loop system (A.1)-(A.2) is We say that the triplet (σ-stabilizes quasi polynomial (A.3) ifwith , where is the real part of In [16], the parameter space of control law [13] to σ-stabilize system [12] is determined.

Lemma A.1 (see [16]). Let be such that and , whereThen, the parametric equations delimiting the regions of σ-stabilizability of closed-loop system (A.1)-(A.2) are given by the following:
Lower Boundaryfor defined as for the upper boundary.
Upper Boundary

In Figure 12, the upper boundary (solid line) and lower bounds (dashed line) using and the parameters given in Table 1 are shown. These bounds delimit the -stability regions for quasi polynomials of the form (A.3).

Also, note that these regions collapse into a single point (red asterisk) which correspond to the unique maximum achievable decay . The following result characterized analytically this point.

Lemma A.2 (see [16]). Let be such that . Then, the maximum exponential decay of the closed-loop system (A.1)-(A.2) is given in (A.5). Moreover, the values that σ-stabilize system (A.1)-(A.2) with the best possible exponential decay areFor illustrative purposes, for Figure 12, we take .

B. Summary of Mikhailov’s Criterion

Stability conditions stated in terms of hodographs (frequency plots) have a long history. In their most general form, these conditions were given by Cauchy’s principle of the argument around 1829. The graphical use of the argument principle was introduced to the engineering community by Nyquist in 1932 and Mikhailov in 1938 for linear time-independent systems; see [36, 37, 5254] for extensions to the case with time delays. Mikhailov’s criterion gives necessary and sufficient conditions for the stability of a quasi polynomial for a time-delay system. Let us formulate one of some possible versions of this criterion.

Theorem B.1. Consider the following quasi polynomial for a time-delay system: where is a closed curve in the complex plane , , , and , for . If has no roots on the curve , then the number of its roots inside of iswhere denotes the changes of the argument of along C.

Immediately, the above theorem is used to determine stability or instability of quasi polynomials of Section 3, exploiting the well-known fact that the number of unstable roots is always finite (see [55], page 9) and therefore can be trapped by a large enough semicircle in the right-half plane.

Let denote one closed curve in the complex plane, as depicted in Figure 13. Here, is the interval on the imaginary axis and is the semicircle of the radius in the right-half (left-half) complex plane. It follows from the above theorem thatwhere and denote the changes of argument of along and , respectively. Assuming that the radius is sufficiently large, to claim stability, one must just verify that , that is, verify that the number of the turns around the origin in the clockwise direction is equal to the number of turns in the counterclockwise direction.

It should be noted, however, that application of Mikhailov’s method in the case of quasi polynomials must be done with great care, since there is typically no monotonicity of rotation of the curve along , as in Figure 13, which can be proven only for the case of polynomials with a finite number of roots.

Nomenclature

Angular position of the pendulum
Angle of the rotational rod
Mass of the rotational rod
Pendulum mass
Moment of inertia of the rod
Moment of inertia of the pendulum
Length of the horizontal rod
Pendulum length
Distance to rod center of gravity
Distance to pendulum center of gravity
Gravity acceleration
Induced torque.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.