Mathematical Problems in Engineering

Mathematical Problems in Engineering / 2015 / Article
Special Issue

Control Problem of Nonlinear Systems with Applications

View this Special Issue

Research Article | Open Access

Volume 2015 |Article ID 589184 |

Fabio Abel Gómez Becerra, Víctor Hugo Olivares Peregrino, Andrés Blanco Ortega, Jesús Linares Flores, "Optimal Controller and Controller Based on Differential Flatness in a Linear Guide System: A Performance Comparison of Indexes", Mathematical Problems in Engineering, vol. 2015, Article ID 589184, 10 pages, 2015.

Optimal Controller and Controller Based on Differential Flatness in a Linear Guide System: A Performance Comparison of Indexes

Academic Editor: Wenguang Yu
Received27 Aug 2015
Revised20 Nov 2015
Accepted23 Nov 2015
Published27 Dec 2015


The use of linear slide system has been augmented in recent times due to features granted to supplement electromechanical systems; new technologies have allowed the manufacture of these systems with low coefficients of friction and offer a variety of types of sliding. In this paper, we present a comparison between the performance indexes of two techniques of control applying optimal control LQR (Linear Quadratic Regulator) acronym for STIs in English and the technique of differential flatness controller. The use of linear slide bolt of potency takes into account the dynamics of the DC motor; the Euler-Lagrange formalism was used to establish the mathematical model of the slide. Cosimulation via the MATLAB/Simulink-ADAMS virtual prototype package, including realistic measurement disturbances, is used to compare the performance indexes between the LQR controller versus differential flatness controller for the position tracking of linear guide system.

1. Introduction

The linear slides have been used in different electromechanical systems as actuators; these have great advantages, especially those consisting of screws-and-nuts; some advantages are as follows: the effect of gravity at the beginning of the movement introduces no disturbance of change of position by not overcoming the power screw, the degree of accuracy in positioning is very high, it is possible to track both soft paths laws of classic and modern control, and the force developed by using this type of device is very high usually requiring actuators (DC motor) of lower power than those requiring other types of systems. Blanco Ortega et al. [1] proposed a rehabilitation ankle device based on an table which uses two linear slidings, one on each axis, through a Generalized Proportional Integral (GPI) controller. Valdivia et al. [2] proposed a rehabilitation ankle TobiBot, which covers only the movements of dorsiflexion/plantarflexion, and a degree of freedom; it is controlled by an outline of PID control, to perform the movements, and uses a linear slide based on a power screw. Blanco Ortega et al. [3, 4] in turn presented a rehabilitation of the design and construction of an ankle rehabilitation based on a parallel robot of 3 degrees of freedom, which provides the movements of dorsiflexion/plantarflexion and inversion/eversion made by the ankle and uses a PID control technique to perform rehabilitation ankle movements by using three linear power sliding screws. Blanco Ortega et al. [3, 4] have presented an ankle rehabilitation machine, using a linear slide to effect movements dorsiflexion/ankle plantarflexion; the control technique using a PID was smooth trajectory tracking.

As can be seen, linear slides have been widely used as actuators in various robotic devices from a simple table. Various control techniques have been implemented for accurate position with smooth tracking paths. The present work proposes control technique for differential flatness, compared with optimal control, showing both its advantages and disadvantages over the other, considering the control technique for differential flatness as an option with highly acceptable results in tracking soft paths, and highlighting also the lack of optimization for optimum control trajectories mechanical tracking systems of this type; the randomness of the matrix appears when weighing matrixes and . Control technique by feedback employed in this work is based on the concept of differential flatness, which come from differential plan systems (Fliess et al. [5]). This was made known fifteen years ago in France by professor Fliess and his collaborators. Differential flatness has had important uses within the areas of robotics, control processes, aerospace systems, optimization systems, trajectory planning in linear and nonlinear aspects, and systems of infinite dimensions described in partially controlled differential equations with border conditions (Fliess et al. [5, 6]; Linares-Flores and Sira-Ramírez [7]; and Sira-Ramírez and Agrawal [8]). Thounthong et al. [9] suggest the use of a differential based on flatness of a fuel cell system and a hybrid supercapacitors source achieving robustness, stability, and efficiency of the controlled system controller. Jörgl and Gattringer [10] proposed the control of a conveyor belt using a control law based on differential flatness reducing the trajectory tracking error compared with traditional drivers, flatness theory has been used in a variety of nonlinear systems in various engineering disciplines, Thounthong and Pierfederici [11], such as the inverted pendulum control and aircraft vertical rise and fall, Fliess et al. [12]. Danzer et al. [13] proposed driver in such a control system of the pressure of a cathode and oxygen excess of a chemical system. Gensior et al. [14] used a control of tracking of a DC voltage boost converter. Song et al. [15] have shown that based on flatness control is robust and provides improved performance monitoring transience compared to a traditional method of linear control (PI). A nonlinear system is flat if there is a set of independent variables (differentially equal in number to the number of entries) so that all the state variables and input variables can be expressed in terms of those output by Syed et al. [16], Rabbani et al. [17], and Agrawal et al. [18].

In order to make a comparison between the performance index of different controllers, this paper also presents an optimal control law applied to the same linear slide (Figure 1); optimal control theory focuses on the design of controllers to perform their objective and concurrently satisfy physical constraints to optimize predetermined performance criteria (Hassani and Lee [19]). With the new trend of seeking a high-performance, sustainable manufacturing, pollution awareness and finding ways for greater energy efficiency, greater emphasis on the optimal design of control systems is made. The optimality design criteria may include minimum fuel, low energy, minimum time (Lewis [20]). Because of this the focus in recent years has been directed to the use of various techniques of optimal control; Yu and Hwang [21] have presented an LQR (Linear Quadratic Regulator) approach for determining a control law PID optimal in order to control the speed of a DC motor; this contribution proposes a systematic approach to design a speed control of a DC motor based on an identification model and LQR design with a nonlinear increase with feedforward compensator. Ruderman et al. [22] propose a methodical approach LQR state feedback control of a DC motor. Likewise LQR optimal control has been used in conjunction with other control techniques such as that of Prasad et al. [23] who proposed a control system high nonlinearity such as the inverted pendulum. It is presented with a linearized dynamics using a PID controller LQR bringing out results in a robust control scheme for optimal system control.

The contributions of this paper are as follows: (1) the design of an LQR position tracking controller; (2) the design of a differential flatness position tracking controller; these tracking controllers are compared on the performance index of position tracking error. Thus, we conclude that the two controllers have a high capacity for the position tracking of the linear guide system. This paper is organized as follows. Section 2 presents the mathematical model of linear slide system. In this section, we obtained the mathematical model via Euler-Lagrange formalism and incorporate the dynamic of DC motor to model. The position control based on linear quadratic regulator is presented in Section 3. In Section 4, we present the design of tracking controller based on differential flatness. In Section 5, the simulations results are obtained through the MATLAB/Simulink-ADAMS virtual prototype package, and we compare the performance indexes of the position tracking error of both controllers. In Section 6, the results are shown in the experiment using a physical slider and data acquisition card, using the LabView software with a graphical interface. Finally, in Section 7, we give the conclusions of all of the work.

2. Mathematical Model of Linear Slide

The linear slide control in this work is formed with a motor coupled CD to a power screw through a gearbox speed, screw, and rotating, linearly moving mass , as shown in Figures 1 and 2.

Motor mathematical model of linear guide system was obtained by applying the Euler-Lagrange formalism. It considers the dynamics of the DC motor. The generalized coordinates are and .

Consider Figure 2 and the notation shown in Notation.

2.1. Electric Motor

Considering Figure 2, the coenergy by storage of effort is given by (1), the coenergy storage of effort by (2), and the energy dissipation by (3).

Coenergy by Storage of Effort. ConsiderThe Storage of Energy Flow. Consider

Energy Dissipation. Consider

2.2. Linear Slide

The kinetic energy of the system of linear slide is given by (4) and the dissipated energy is represented by (7).

Kinetic Energy. ConsiderKnowing that , , and and substituting in (4), the following equation is obtained, with kinetic energy remaining in function of the angular velocity of the power screw (see Figure 2):Dissipated Energy. ConsiderFrom (1) until (5) is represented total system in the kinetic energy and the following equation is obtained:where is the Lagrangian and is given byFor the generalized coordinate ,where

Considering that and as well as and replacing them in (10), the following equation is obtained:Knowing thatthe following equation is obtained:For the generalized coordinate ,Equation of Euler-Lagrange formalism isThe mathematical model of the linear slide taking into account the dynamics of the DC motor is defined by

3. Optimal Control

3.1. Considering Acceleration State Variable

Based on the mathematical model (16) and (17), is derived and from (17) it is clear that Substituting in (16), one hasFixing the equation, one getsThe state variables arePlacing the state variables in the matrix form, one getsIn control systems, often we want to select the control vector such that a given performance index is minimized. A quadratic performance index, where the integration limits are to , so thatwhere is a quadratic function or a Hermitian function of and , produces linear control laws; that is to say,For weight matrixes for semipositive definite and , the optimal control system is based on minimizing the performance index. This requires numerically solving the Riccati algebraic equation: for a symmetric positive definite matrix .

Finally, gains are calculated asIn the present case, is given by The state weighting matrix is proposed aswhich meets nonnegative definite. As for the scalar weighting for the control input, it is chosen as

3.2. Considering State Variable Motor Current CD

By representing mathematical model equations (16) and (17) in state space, one getswithfor the first case and for the second case.

3.3. LQR Optimal Control in Virtual Prototype

A test was performed using a virtual prototype of the linear slide, Figure 3, in room ADAMS MSC together with Matlab-Simulink in order to test the effectiveness of optimal control LQR; the results of this experiment can be seen in Section 5.

4. Control Based on Differential Flatness

This third-order linear system (30), where its Kalman controllability matrix is calculated by the following expression: , is given asSince the determinant is nonzero, then the system is controllable and, therefore, is differentially flat (Sira-Ramírez and Agrawal [8]). The flat output of a linear system input output (I/O) is obtained by multiplying the inverse matrix of controllability by the state vector , associated with the system. Column vector obtained by multiplying the last line is chosen to obtain the flat output (Linares-Flores and Sira-Ramírez [7]). In particular for reducing motor-drive CD, flat output system is calculated asHence, we have chosen as the flat output. This flat output provides the following differential parametrization of the system variables: and the control input:

Considering (36) and in order to verify the stability of the system, applying the Laplace transform, the following holds:

Taking into account the values of Table 1, the system transfer function isPlotting the locus of roots verifies that the system is stable when these are in the left half of the complex plane (Figure 4).

Parameter Value

= speed ratio 0.19
= force opposite to the movement of 0.01 N
= coefficient of viscous friction2
= pitch of the screw thread power0.00196 m
= mass to be displaced10 kg
= moment of inertia0.0000014 kg m2
= voltage12 volts
= constant emf0.022 (Vs)/rad
= resistance of DC motor5.3 ohm
= motor inductance0.00058 henries
= constant torque90 (N-m)/A

Hence, we have the fact that all state variables and the control input are in terms of and its successive derivatives, where it denotes the position of the mass as . The differential parametrization above allows obtaining the equilibrium for the system in terms of the equilibrium values of the flat output and the disturbance inputs. Thus, From (36), we design the average controller based on differential flatness property. Thus, we replace the higher-order derivative of the flat output by a virtual controller (see Slotine and Li [24]) resulting in the following:For the tracking controller design, we use a nominal desired linear displacement profile that exhibits a rather smooth start for the motor linear slide system. This is specified using an interpolating Bézier polynomial of 10th order where the initial linear displacement is set to be  m valid until  sec and the final desired value of the angular velocity is specified as  m to be reached at  sec; that is, we usedwhere is a polynomial function of time, exhibiting a sufficient number of zero derivatives at times and , while also satisfying and . For instance, one such polynomial may be given byThe differential parametrization of the control input, , in terms of shows that the proposed flat output trajectory tracking task is that of controlling the third derivative of by means of :Therefore, the linear displacement-tracking controller isUnder these circumstances, the closed loop system tracking error, , satisfies the linear differential equationThe appropriate choice of the constant coefficients , as coefficients of the third-order Hurwitz polynomial, guarantees the asymptotic exponential stability to zero of the tracking error, . One such choice, yielding a characteristic polynomial of the form with , , and , is given byThe virtual controller, , achieves a smooth start for the linear displacement of the system. If we apply an unknown constant load torque on the system, we have to include a term of integral action into the virtual control (43). Thus, we minimize the tracking error near to zero.

5. Simulation Results

These results were obtained using the parameter values of Table 1.

In Figure 5, one can see the movement of the mass 0.5 meters in 60 seconds following a path originated by a tenth-order Bézier polynomial; the tracking error can be seen as zero in the whole path; in this case, the acceleration is taken as a state variable.

Figure 6 shows the graphs of results taking into account the current as state variable; considering the matrix and considering the matrix , Figure 7, in both cases the tracking error can be seen as zero in the whole path.

Figure 8 shows the results of simulation based on differential flatness controller observed; it also shows that it is capable of moving the mass to a position 0.5 meters along a desired path of a Bézier polynomial of tenth order.

The results of the performance index are as follows: the two controllers are shown in Figures 9, 10, and 11; the indexes of performance of the optimal controller for two values of the matrix are shown.

Figure 12 shows the results using a virtual prototype environment ADAMS MSC and MATLAB-Simulink; we can see that the offset for this test was set at 0.7 meters corresponding to the desired position; a tenth-order Bézier polynomial was used as desired path.

6. Results with Linear Slide

Experimentation was held using a physical linear slide (Figure 13); the laws of LQR control are implemented based on differential flatness with the same parameter values and gains controllers; the results of the experiment can be seen in Figures 14, 15, and 16, where it can be seen that in all three cases the controller is able to zero the position error; it is clear that the experiment was carried out with an acquisition card myRIO data of National Instruments; the measurement units were centimeters so that they should be taken into account in the interpretation of Figures 14, 15, and 16.

It is verified by the results the correct operation for achieving desired trajectories and references.

On the -axis is the number of samples; to determine the time it must be multiplied by the sampling period, that is, 0.1 s.

7. Conclusions

By applying control techniques and optimal LQR and based on differential flatness it can be seen that the results for the trajectory tracking are highly achievable in both cases, with regard to the performance indices of each controller for the optimal controller; a great disadvantage exists since it is necessary that the values of the weighting matrix be too low to achieve a speed of response adapted to the needs of the plant to controlling; this means that the gains are higher further input values, such as the voltage that must be adjusted by trial and error of said matrix with the purpose of respecting the voltages of the actuators.

Figures 12 and 13 show that the values of the performance index are minimizing in the measure that the value of the weighting matrix, , is minor compared to the index controller performance differential flatness, if it is minor, but this is not guaranteed to be optimal because the values of the control force are lower than those required for the good functioning of DC motors.

In conclusion it follows that the use of control laws differential flatness is by far better than the use of laws of optimal control; the concept of optimality, in this particular case, is lost due to the randomness representing search for the value of the matrix to obtain appropriate response speed without sacrificing system actuators.


Speed ratio
:Force opposite to the movement of
:Coefficient of viscous friction
:Pitch of the screw thread power
:Mass to be displaced
:Moment of inertia
:Constant emf
:Resistance of DC motor
:Motor inductance
:Constant torque
:Angular position of DC motor
:Angular velocity of DC motor
:Angular position of power screw
:Angular velocity of power screw
:Torque delivered by the DC motor
:Torque delivered by the speed reducer
:Displacement of the mass
:Velocity of the mass
:Force of displacement of the mass.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


  1. A. Blanco Ortega, F. A. Gómez-Becerra, L. G. Vela Valdés, and R. O. Delgado Arcega, “A generalized proportional integral controller for an ankle rehabilitation machine based on an XY table,” in Proceedings of the IEEE International Conference on Mechatronics, Electronics and Automotive Engineering (ICMEAE '13), pp. 152–157, IEEE, Morelos, Mexico, November 2013. View at: Publisher Site | Google Scholar
  2. C. H. G. Valdivia, J. L. C. Escobedo, A. B. Ortega, M. A. O. Salazar, and F. A. G. Becerra, “Diseño y control de un sistema interactivo para la rehabilitación de tobillo: TobiBot,” Ingeniería Mecánica. Tecnología y Desarrollo, vol. 5, no. 1, pp. 255–264, 2014. View at: Google Scholar
  3. A. Blanco Ortega, J. Isidro Godoy, E. Quintero Mármol, and L. G. Vela Valdés, “Robot paralelo para rehabilitación asistida de tobillo,” in X Congreso Internacional sobre Innovación y Desarrollo Tecnológico (CIINDET '13), Cuernavaca, Mexico, March 2013. View at: Google Scholar
  4. A. Blanco Ortega, H. R. Azcaray Rivera, R. F. Vázquez Bautista, and L. J. Morales Mendoza, “Máquina de Rehabilitación de Tobillo: prototipo virtual y físico,” in Proceedings of the 10th Congreso Internacional sobre Innovación y Desarrollo Tecnológico (CIINDET '13), vol. 1, Cuernavaca Morelos, México, March 2013. View at: Google Scholar
  5. M. Fliess, J. Lévine, P. Martin, and P. Rouchon, “Flatness and defect of non-linear systems: introductory theory and examples,” International Journal of Control, vol. 61, no. 6, pp. 1327–1361, 1995. View at: Publisher Site | Google Scholar | MathSciNet
  6. M. Fliess and R. Marquez, “Continuous-time linear predictive control and flatness: a module-theoretic setting with examples,” International Journal of Control, vol. 73, no. 7, pp. 606–623, 2000. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  7. J. Linares-Flores and H. Sira-Ramírez, “DC motor velocity control through a DC-to-DC power converter,” in Proceedings of the 43rd IEEE Conference on Decision & Control, vol. 5, pp. 5297–5302, Nassau, Bahamas, December 2004. View at: Publisher Site | Google Scholar
  8. H. Sira-Ramírez and S. K. Agrawal, Diěerentially Flat Systems, Marcel Dekker, Inc, 2004.
  9. P. Thounthong, S. Pierfederici, J.-P. Martin, M. Hinaje, and B. Davat, “Modeling and control of fuel cell/supercapacitor hybrid source based on differential flatness control,” IEEE Transactions on Vehicular Technology, vol. 59, no. 6, pp. 2700–2710, 2010. View at: Publisher Site | Google Scholar
  10. M. Jörgl and H. Gattringer, “Dynamical modeling and flatness based control of a belt drive system,” Proceedings in Applied Mathematics and Mechanics, vol. 14, no. 1, pp. 887–888, 2014. View at: Publisher Site | Google Scholar
  11. P. Thounthong and S. Pierfederici, “A new control law based on the differential flatness principle for multiphase interleaved DC–DC converter,” IEEE Transactions on Circuits and Systems II: Express Briefs, vol. 57, no. 11, pp. 903–907, 2010. View at: Publisher Site | Google Scholar
  12. M. Fliess, J. Lévine, P. Martin, and P. Rouchon, “A Lie-Bäcklund approach to equivalence and flatness of nonlinear systems,” IEEE Transactions on Automatic Control, vol. 44, no. 5, pp. 922–937, 1999. View at: Publisher Site | Google Scholar | MathSciNet
  13. M. A. Danzer, J. Wilhelm, H. Aschemann, and E. P. Hofer, “Model-based control of cathode pressure and oxygen excess ratio of a PEM fuel cell system,” Journal of Power Sources, vol. 176, no. 2, pp. 515–522, 2008. View at: Publisher Site | Google Scholar
  14. A. Gensior, H. Sira-Ramirez, J. Rudolph, and H. Guldner, “On some nonlinear current controllers for three-phase boost rectifiers,” IEEE Transactions on Industrial Electronics, vol. 56, no. 2, pp. 360–370, 2009. View at: Publisher Site | Google Scholar
  15. E. Song, A. F. Lynch, and V. Dinavahi, “Experimental validation of nonlinear control for a voltage source converter,” IEEE Transactions on Control Systems Technology, vol. 17, no. 5, pp. 1135–1144, 2009. View at: Publisher Site | Google Scholar
  16. F. U. Syed, M. L. Kuang, M. Smith, S. Okubo, and H. Ying, “Fuzzy gain-scheduling proportional–integral control for improving engine power and speed behavior in a hybrid electric vehicle,” IEEE Transactions on Vehicular Technology, vol. 58, no. 1, pp. 69–84, 2009. View at: Publisher Site | Google Scholar
  17. T. Rabbani, S. Munier, D. Dorchies, P.-O. Malaterre, A. Bayen, and X. Litrico, “Flatness-based control of open-channel flow in an irrigation canal using SCADA,” IEEE Control Systems Magazine, vol. 29, no. 5, pp. 22–30, 2009. View at: Publisher Site | Google Scholar | MathSciNet
  18. S. K. Agrawal, K. Pathak, J. Franch, R. Lampariello, and G. Hirzinger, “A differentially flat open-chain space robot with arbitrarily oriented joint axes and two momentum wheels at the base,” IEEE Transactions on Automatic Control, vol. 54, no. 9, pp. 2185–2191, 2009. View at: Publisher Site | Google Scholar
  19. K. Hassani and W.-S. Lee, “Optimal tuning of linear quadratic regulators using quantum particle swarm optimization,” in Proceedings of the International Conference on Control, Dynamic Systems, and Robotics (CDSR '14), Paper no. 59, Ottawa, Canada, May 2014. View at: Google Scholar
  20. F. L. Lewis, Optimal Feedback Control: Practical Performance and Design Algorithms for Industrial and Aerospace Systems, UTA Research Institute The University of Texas at Arlington, USA, 2012.
  21. G.-R. Yu and R.-C. Hwang, “Optimal PID speed control of brush less DC motors using LQR approach,” in Proceedings of the IEEE International Conference on Systems, Man and Cybernetics (SMC '04), pp. 473–478, Hague, The Netherlands, October 2004. View at: Publisher Site | Google Scholar
  22. M. Ruderman, J. Krettek, F. Hoffmann, and T. Bertram, “Optimal state space control of DC motor,” in Proceedings of the 17th IFAC World Congress, International Federation of Automatic Control, Seoul, Republic of Korea, July 2008. View at: Google Scholar
  23. L. B. Prasad, B. Tyagi, and H. O. Gupta, “Optimal control of nonlinear inverted pendulum system using PID controller and LQR: performance analysis without and with disturbance input,” International Journal of Automation and Computing, vol. 11, no. 6, pp. 661–670, 2014. View at: Publisher Site | Google Scholar
  24. J.-J. E. Slotine and W. Li, Applied Nonlinear Control, Prentice Hall, 1991.

Copyright © 2015 Fabio Abel Gómez Becerra et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.