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

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.


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 highperformance, 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.

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.
Consider The Storage of Energy Flow.Consider Energy Dissipation.
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).
Knowing that  =  2 , ẋ =  θ 2 , and ẍ =  θ 2 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.
From ( 1) until ( 5) is represented total system in the kinetic energy and the following equation is obtained: where  is the Lagrangian and is given by For the generalized coordinate , where  q +  q =  () − .(10) Considering that  =   θ 1 and θ 1 =  θ 2 as well as θ 2 = ẋ / and replacing them in (10), the following equation is obtained: Knowing that the following equation is obtained: For the generalized coordinate  2 , Mathematical Problems in Engineering Equation of Euler-Lagrange formalism is The mathematical model of the linear slide taking into account the dynamics of the DC motor is defined by

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 has Fixing the equation, one gets The state variables are ) .
In 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 0 to ∞, so that where (, ) 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 as In the present case,  is given by The state weighting matrix is proposed as  16) and ( 17) in state space, one gets  for the second case.

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.

Control Based on Differential Flatness
This third-order linear system (30), where its Kalman controllability matrix is calculated by the following expression:  = [, , Since 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 as Hence, we have chosen  1 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 is Plotting the locus of roots verifies that the system is stable when these are in the left half of the complex plane (Figure 4).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  ini = 0 m valid until  ini = 15 sec and the final desired value of the angular velocity is specified as  fin = 0.5 m to be reached at  fin = 45 sec; that is, we used where   (,  ini ,  fin ) is a polynomial function of time, exhibiting a sufficient number of zero derivatives at times  ini and  fin , while also satisfying   ( ini ,  ini ,  fin ) = 0 and   ( fin ,  ini ,  fin ) = 1.For instance, one such polynomial may be given by The 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 V aux : Therefore, the linear displacement-tracking controller is Under these circumstances, the closed loop system tracking error,  =  −   , satisfies the linear differential equation The appropriate choice of the constant coefficients { 2 ,  1 ,  0 }, 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 ( + )( 2 + 2   +  2  ) with  > 0, 0 <  < 1, and   > 0, is given by The virtual controller, V aux , 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.

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 tenthorder 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  = 0.01 and considering the matrix  = 1 − 10, 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 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.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.

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.

1 Figure 2 :
Figure 2: Schematic diagram of a linear slide system-DC motor.
definite.As for the scalar weighting for the control input, it is chosen as  = 0.01,  = 1 − 10.(29)3.2.Considering State Variable Motor Current CD.By representing mathematical model equations (

Figure 3 :
Figure 3: Virtual prototype of the linear slide.

Notation𝑛 : 1 : 2 :
Speed ratio : Force opposite to the movement of   2 : Coefficient of viscous friction : Pitch of the screw thread power : Mass to be displaced : Momentofinertia : Voltage   : Constant emf : Resistance of DC motor : Motorinductance  tao : Constant torque  1 : Angular position of DC motor θ Angular velocity of DC motor  2 : Angular position of power screw θ Angular velocity of power screw  1 : Torque delivered by the DC motor  2 : Torque delivered by the speed reducer : Displacement of the mass ẋ : Velocityofthemass : Force of displacement of the mass.