Abstract

The trajectory tracking task of a magnetic levitation system connected to a beam mechanism is solved by means of a nested saturation control with a feedforward term. The flatness property of the system allows to use the nested saturation control technique and the feedforward control to stabilize the output tracking error around the equilibrium. The closed-loop error dynamics is proven to be locally exponentially stable. Numerical simulations prove the effectiveness of the proposal.

1. Introduction

Underactuated systems are currently an active research area in automatic control design. On the one hand, their amount of advantages such as lighter structures, economic designs, the possibility of compensating a failure on a fully actuated system, among others, has allowed their use in many applications such as mobile robotics, marine systems, aerospace robotic systems, cranes, flexible robots, and fault compensation systems [14]. On the other hand, the control problems of regulation and trajectory tracking of underactuated systems lead to complex solutions which cannot be solved by classical schemes [5]; some of the complexities are due to the difficulties to find general properties concerning their capacity of being linearizable, relations between degrees of freedom and independent control actuators [6], passivity [7, 8], ill-defined relative degrees [9], and so on. Most of the solutions for this class of systems deal with the regulation problem, where energy-based schemes are the most popular. The trajectory planning tasks and their tracking demand more complex strategies. This is due to the fact that not all joint trajectories are attainable for this class of systems [10]. To overcome this problem, the capacity of finding a set of variables that can parameterize the system trajectories (differentially flatness [1113]) allows to establish a trajectory planning and a subsequent controller synthesis.

In particular, magnetic levitation systems have been extensively analysed and used in the design of frictionless bearings [14], vibration isolation [15], manipulation and micromanipulation systems [16, 17], drug and people transporting [18, 19], energy harvesting [2022], etc., where noncontact motion control or avoiding the use of lubrication systems is desirable [2325]. The control of this class of systems has been tackled from a wide variety of nonlinear approaches [26, 27]. The magnetic levitation in a beam balance configuration has been approached by feedback linearization [28]. In [29], the regulation problem of a beam balance system is solved though a Lyapunov-based control considering restrictions in states and input. Hu et al. [30] developed a saturated linear feedback to solve the same problem achieving a larger attraction region with respect to linear techniques. Passivity-based control is also implemented though the interconnection and damping assignment [31]. In [32], a sliding mode controller is used for the regulation problem using two magnetic actuators in an agonist-antagonist configuration. The problem of trajectory tracking for a ball levitation system controlled by magnetic bearings has been solved by a flatness-based control [3335], where path planning is used to avoid the input constraints in the control design although the natural saturations are not part of the control design in contrast with other schemes. The idea of using the saturations as a part of the control design in combination with the advantages of a flatness-based control seems to be a good alternative to solve the problem of trajectory tracking enhancing the control response in case the control saturation may arise. The nested saturation function approach was introduced by A. Teel in [36] to solve the stabilization of a chain of integrator of arbitrary order, and it can also be used for trajectory tracking for a class of trajectories restricted. Other works that solve this problem are presented in [3739]. Besides, this technique has been used for stabilization of a class of underactuated systems [4043].

In this sense, this article deals with the problem of trajectory tracking of a beam balance levitation system by means of a nested saturation function approach in combination with the flatness property of the system. The proposal consists in using the flatness of the system to transform the system in an integrator chain, and then, a controller based on a nested saturated function and a feedforward term is used to force the output trajectories to converge towards the reference trajectory. The closed-loop solution is proven to be locally exponentially stable by means of the second method of Lyapunov.

The remainder of the article is given as follows: the dynamical model of the magnetic levitation system, the flatness property, and the problem formulation are given in Section 2. The control proposal and its stability proof are provided in Section 3. Some numerical simulations which show the behaviour of the proposal are presented in Section 4, and finally, some concluding remarks are stated.

2. System Model

Consider the magnetic levitation system, shown in Figure 1, where a beam can rotate freely, and its movement is affected by the gravity effects and the control input derived from a magnetic force generated on a controlled electromagnet coil. The dynamical model is given as follows:where is the angular position of the pendulum, is its corresponding angular velocity, and is the current of the inductor. represents the electromagnetic force proportionality constant. is the inertia moment of the beam, is the distance between the pivot and the center of the permanent magnet, denotes the mass of the beam, is the gravity constant, and is the radial distance of the center of mass of the beam to the rotational center. and are the inductance and resistance parameters of the electromagnetic subsystem. Finally, is the control input, applied on the electromagnetic system. It is assumed that the output variable is the position of the beam .

Last system admits a canonical controllable form by means of the following change of coordinates [44]:where and .

We have the following inverse transform:

The dynamics of (1) in terms of coordinate transformation (2) is given as follows:where

This system can be controllable except when or in the set .

2.1. Problem Formulation

The following problem formulation is established: given a smooth admissible reference trajectory of the beam angle, denoted as , devise an output feedback control law such that the trajectory tracking error remains after a settling time on a vicinity of the origin of the phase variable plane, as small as allowed by the system nonlinearities and the control gains to be designed.

3. Nested Function Control Design

From (4), it can be shown that the system is differentially flat [11, 12] with flat output . That is, system (4) has the formwith

The following linearizing control inputyields the following linear controllable in the Brunovsky canonical form

Let us define the flat output error as , where . In order to express the dynamics in terms of the tracking error coordinates, the following feedforward input term is introduced in the last expression:

Thus, using (8) and (10) in (4), the following expression is obtained:

Last system admits the following canonical representation:

Let us introduce the following linear transformation [36, 45]:

Last transformation leads system (12) to

Therefore, as system (12) is expressed with a chain of integrator, we suggested a controller based on nested saturation function. This technique, proposed in [36], has been used for controlling a wide class of the underactuated system [4042, 46]. Thus, our stability problem will be solved as follows. First, a linear transformation is used to directly propose a trajectory tracking controller; then, it is shown that the proposed controller guarantees the boundedness of all states, and after a finite time, it is possible to ensure that all states converge to zero [43].

Let us define the following saturation function:

Definition 1. The linear saturation function is defined asFinally, the trajectory tracking controller can be proposed aswhere is a positive constant, and are fixed parameters, and is a linear saturation function.

3.1. Boundedness of All States

We show in three simple steps that the closed-loop solution of the proposed closed-loop systems (14) and (16) ensures that all the states are bounded. Moreover, the bound of each state directly depends on the designed parameters of the controller.

Step 1. We define a positive definite functionThen, differentiating (17) and using the third differential equation of (14), we have the time derivative of given byIt is clear that when ; therefore, there is a finite time such thatwhere .

Step 2. Now, we proceed to analyze the behaviour of the state . Hence, we introduce a positive definite function . Differentiating , we obtain after substituting (16) into second differential equation of (14)where and are selected such that . Evidently, if , then , and there is finite after whichwhen this condition is satisfied, and the control turns out to be

Step 3. Substituting (22) into first differential equation of (14), we obtainNow, we define an auxiliary positive definite function . By differentiating along the trajectories of (23), we obtainwhere must be chosen such that . If , then , and hence, there is finite after whichConsequently, is also bounded. So, all previous constraints on parameters and can be summarized asManipulating the last inequalities, we can select the control parameters as follows:where .

3.2. Convergence of All States to Zero

We shall prove that the closed-loop system given by (14) and (16) is asymptotically stable and locally exponentially stable, provided that the controller parameters satisfy (26).

We must note that, after , the control law is no longer saturated, that is,and the closed-loop system can be expressed as

Let us define the following Lyapunov function:with . Now, differentiating along the trajectories of (29), we obtainwhere

Note that , and therefore, M is positive definite Therefore, is negative definite, and the vector state locally exponentially converges to zero after .

From the above discussion, we have the following.

Proposition 1. Consider the magnetic levitation system as described (1) in a closed loop with controller (8). Then, the closed-loop system is asymptotically stable and locally exponentially stable provided that the control parameters and satisfy the inequalities.

4. Numerical Simulations

To test the performance of the controller, we carried out some numerical simulations using MATLAB program, and the results were obtained based on the numerical method of Runge–Kutta of fourth order with the fixed step of 0.001 s. The physical parameters of the system are , (kg ), (m), (kg), (m/s2), (m), (H), and (), and the controller parameter values were set as , , and .

The first experiment shows the behaviour of the system, for a rest-to-rest trajectory, and in this case, we use a reference trajectory-type Bezier polynomial . The initial conditions were set as , , and . In Figure 2, the output tracking response is shown, and we can see that the controller effectively brings the system from initial position () to the final rest position after 7 seconds. The behaviour for the states and is shown in Figure 3. Finally, Figure 4 depicts the control input.

The second experiment was designed under the same setup, but in this case, we use a reference sinusoidal trajectory . The result of the experiment is presented in Figure 5. From this figure, we can see that the controller was capable of performing the task of tracking a sinusoidal trajectory. In Figure 6, the behaviour for the states and is shown. Finally, Figure 7 presents the control input.

In Figure 8, the output tracking response for a sinusoidal trajectory with different initial conditions is shown. In this case, , and we can see that the controller was capable of performing the task of tracking a sinusoidal trajectory adequately. Some tests were carried out in order to verify the sensitivity of the control scheme with respect to the initial conditions, in which a set of different initial conditions in an admissible operation range was used. The results indicate that the response is not affected by the initial conditions, and there is a lack of overshooting effects which are typical in high-gain schemes.

4.1. A Comparison Test

In order to assess the behaviour of the control proposal, a comparison test was made against a robust control strategy, consisting in a sliding mode controller (see [47]). The sliding mode controller was set to be

The control parameters used in the test were , , and . Figure 9 shows the tracking results, where the sliding mode results (denoted by the subindex ) converged faster without any overshooting effects. Figure 10 depicts the states and , in which both schemes had similar values (reachable for experimental implementations). Figure 11 shows the advantages of the proposal since the voltage input of the sliding mode showed the high-frequency behaviour, which is a classic problem in sliding mode control implementations which may demand high control efforts. Notice that the proposal has a better energy management, leading to similar results with a smoother control input. This can be noticed with the obtained performance index which is shown in Figure 12.

5. Conclusion

A nested saturation-based controller for the trajectory tracking task on a beam levitation system was introduced and proven to be effective. The flatness property of the magnetic levitation system permits to express the system as a third-order integrator chain, allowing the use of nested saturation functions to design a trajectory tracking controller.

The proposed scheme forces the system output to track the reference trajectory even in the presence of the input constraints and ensures that all states converge exponentially to the desired trajectory. The stability analysis is based on the second method of Lyapunov using a simple candidate function. Some computer simulations showed the effectiveness of the proposal in the tracking of a rest-to-rest and sinusoidal trajectories.

On the contrary, as a future research, some other nonlinearities such as hysteresis in the magnetic actuator [48] can be addressed. This effect is important for the design of maglev systems in superconductor-based transportation system applications.

Data Availability

No data were used to support this study.

Conflicts of Interest

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

Acknowledgments

This article was supported by the Secretaría de Investigación y Posgrado of the Instituto Politécnico Nacional (SIP-IPN) under the research grants 20195901, 20196059, 20196079, 20196855 ,20201768, 20201829 and 20201830. J. Alejandro Aguirre-Anaya is a scholarship holder of the CONACYT.