Research Article | Open Access
Integral Sliding Modes with Nonlinear -Control for Time-Varying Minimum-Phase Underactuated Systems with Unmatched Disturbances
This paper presents a methodology for controlling nonlinear time-varying minimum-phase underactuated systems affected by matched and unmatched perturbations. The proposed control structure consists of an integral sliding mode control coupled together with a global nonlinear -control for rejecting vanishing and nonvanishing matched perturbations and for attenuating the unmatched ones, respectively. It is theoretically proven that, using the proposed controller, the origin of the free-disturbance nonlinear system is asymptotically stabilized, while the matched disturbances are rejected whereas the -gain of the corresponding nonlinear system with unmatched perturbation is less than a given disturbance attenuation level with respect to a given performance output. The capability of the designed controller is verified through a flexible joint robot manipulator typically affected by both classes of external perturbations. In order to assess the performance of the proposed controller, an existing sliding modes controller based on a nonlinear integral-type sliding surface is also implemented. Both controllers are then compared for trajectory tracking tasks. Numerical simulations show that the proposed approach exhibits better performance.
Much research in recent years has focused on the stabilization and control of mechanical systems operating under uncertain conditions such as external disturbances, uncertain parameters of the plant, and parasitic dynamics. These problems are present in real-world applications revealing, for example, instability, limit cycles, steady-state error, poor repeatability, or imprecisions.
In spite of the rich and diverse literature on the matter (see, e.g., [1–3]), the unmatched disturbances are still a challenging problem, faced by control engineers, that adversely affect the performance of any system to be controlled. This kind of disturbances cannot be trivially neglected since they can be aroused by unavoidable noise in the measurements or perturbing the output as well. Moreover, disturbances acting in the nonactuated part of an underactuated mechanical system (e.g., pendulums, car-like robots, biped robots, and unmanned aerial vehicles) are a typical example where the unmatched disturbances must be counteracted. Indeed, the problem becomes more complicated for the motion control of this kind of systems since unmodeled dynamics can emerge (see, e.g., ).
Sliding modes are long recognized as a powerful control method to reject vanishing and nonvanishing uniformly bounded matched disturbances and plant uncertainties. However, unmatched disturbances are not counteracted. On the other hand, nonlinear control has the capability of attenuating both matched and unmatched disturbances [5, 6]. There have been many results dealing with unmatched disturbances; however, integral sliding modes have begun to receive a growing interest. For example, Kumar et al.  solve the regulation problem for a Stewart robot using a smooth integral sliding mode (ISM) controller, which drives the position error to the origin in finite time, while the closed-loop system is demonstrated to be robust against matched disturbances only. Mahieddine et al.  propose a sliding mode controller that allows attenuating matched and unmatched uncertainties in nonlinear systems and also reducing the chattering in the control signal. Han et al.  developed output feedback-based sliding mode control schemes for linear time-delayed systems considering both matched and unmatched uncertainties. A detailed revision of ISM addressing the unmatched disturbances is presented in .
In relevant works, Osuna et al.  make a -gain analysis for hybrid mechanical systems operating under unilateral constraints and admitting sliding modes and collision phenomena. Rubagotti et al.  prove that the definition of a suitable sliding manifold and the generation of sliding modes upon it guarantee that matched disturbances are completely rejected while the unmatched ones are not amplified. Besides, a linear control with nonlinear compensation is proposed against matched disturbances. Cao and Xu  present a nonlinear ISM controller where the unmatched disturbances are not amplified, but a linear controller is also used. Castaños and Fridman  show the robustness properties of an integral sliding mode controller ensuring rejection of the matched disturbances, while unmatched perturbations are not amplified using also a linear control. Galvan-Guerra and Fridman  proposed an output ISM for linear time-variant systems, where the main goal is to eliminate the matched perturbation. In , an adaptive ISM control for a class of nonlinear uncertain and invariant systems is proposed to eliminate the quantization sensitivity parameters and matched perturbations, which is accomplished by using an integral sliding function from the local dynamics of the plant. Chen et al.  propose a nonlinear ISM fault tolerant control where an optimal control is used against matched disturbances.
The aforementioned literature includes a linear control to attenuate unmatched disturbances to easily find a suboptimal solution through an algebraic Riccati equation (ARE). On the other hand, a physical phenomenon is better described by its nonlinear dynamic equations. However, the local solution of these dynamic equations is required for the ARE to be solved. Besides, hard computational work is also entailed for verifying the Hamilton-Jacobi-Isaacs inequality and obtaining a global solution.
In this paper, an ISM control combined together with a nonlinear control is presented. The proposed controller allows rejecting matched bounded disturbances and attenuating the effect of the unmatched ones. The synthesized controller, that admits a time-varying input matrix, was applied for solving the tracking control problem for 1 degree of freedom (DOF) flexible joint robot (FJR) manipulator, which consists of a single link interconnected by an elastic revolute joint. The formulation of the nonlinear -control problem is confined to nonautonomous affine systems, and it requires a controller design that guarantees both the internal asymptotic stability of the closed-loop system and its dissipativity with respect to admissible external disturbances. In contrast to previous works, a strict Lyapunov function was proposed in this paper to ensure a global solution of the control problem, by means of the verification of the Hamilton-Jacobi-Isaacs inequality, thus avoiding a hard numerical computation of a partial differential equation or a solution of the corresponding differential Riccati equation [18, 19] where the linearization around the equilibrium point of the plant is required. In order to assess the performance of the proposed controller, we implemented a controller, introduced by Cao and Xu in , based on a nonlinear integral-type sliding surface. The results of both, the proposed controller and the Cao-Xu controller, were compared in a trajectory tracking task.
This paper contributes to the following:(i)Presenting the design of a new controller by combining ISM control and a nonlinear control for time-varying minimum-phase underactuated systems affected by both matched and unmatched disturbances(ii)Developing a rigorous stability analysis, with a global solution to control and verifying the Hamilton-Jacobi-Isaacs inequality(iii)Detailing a procedure to implement the Cao-Xu controller in a 1-DOF FJR manipulator(iv)Presenting a comparative analysis of both controllers by means of numerical simulations with a trajectory tracking task
This paper is organized as follows. Section 2 presents the problem statement and synthesis of ISM control and control for a class of time-varying systems. In Section 3, the combined ISM and nonlinear tracking control is developed for a single pendulum with elastic joint affected by matched and unmatched perturbation. Here, the details regarding the implementation of the Cao-Xu controller are also presented. In Section 4, the performance of both controllers is evaluated in a simulation study. Finally, conclusions are provided in Section 5.
2. Controller Design
This work aims to design a controller for nonlinear time-varying minimum-phase underactuated systems being affected by matched and unmatched perturbations. Let us denote as the desired reference signal and as the output of the nonlinear system to be controlled. The control problem can now be defined as follows.
Control Problem. Given a smooth reference signal () , find a control law such that holds for the free-disturbance case (i.e., ), and the -gain of the perturbed nonlinear system is less than a positive constant level with respect to a given performance output to be controlled.
The block diagram of the closed-loop system with the proposed controller is depicted in Figure 1. The performance output is nothing other than a vector of variables, including the nonmeasurable ones, where the disturbances must be attenuated. The effect of the matched perturbation is canceled out by means of an ISM controller . However, the resulting closed-loop system is still affected by the unmatched perturbation . Then, using an -control, the -gain of the closed-loop system is made less than a given attenuation level .
The next three definitions are provided in order to clarify the -gain concept.
Definition 1 (space , ). The space consists in the set of all piecewise-continuous function such that
Definition 2 ( norm, ). norm of all piecewise-continuous inputs is given by
Definition 3 (-gain, ). Let and be piecewise-continuous functions, denoting the output and the input of a system, respectively. Viewing and as finite-energy signals, -gain is then defined as the ratio between norm of the output and norm of the input bounded by ; that is,and equivalently
Functions in represent signals having finite energy over the infinite time interval and therefore the number in inequalities (4) and (5) can be interpreted as (an upper bound of) the ratio between the energies of output and input.
The steps required to design the controller are detailed in the following.
Consider a nonautonomous nonlinear system of the form where is the time variable, is the state vector, is the control input, , and ; the elements and represent the matched and unmatched perturbations, respectively, which are uniformly bounded in ; that is, where , are positive constants, known a priori.
The functions , , and are assumed to be piecewise continuous in for all and continuously differentiable in for all . The matrix is assumed to be full rank, with image orthogonal to the image of , that is, for all and for all . From now on, it is assumed that fulfills the following assumption.
Assumption 1. for all .
Let us assume the next structure for the controller in (6) is as follows:where will be designed as an ISM controller, with the aim of rejecting the matched perturbation . The nominal control will be designed to drive the system trajectories to the origin, while attenuating the effect due to the unmatched perturbation .
2.1. Design of the Integral Sliding Modes Control
Let us define the sliding surface as with being a constant matrix.
Equation (9) can be seen as a penalizing factor of the difference between the actual trajectories and the trajectories of the system (6) in the absence of matched perturbations and in the presence of the control , projected along . The sliding mode begins from the initial time ; that is, for all . Now we assume the following.
Assumption 2. The product is nonsingular for all .
In order to drive the trajectories of (6) to the sliding surface (9), the following feedback controller is proposed:with the matrix satisfying the condition
Notice that, in (14), the unmatched disturbances are still present. Then, a nonlinear -controller is considered to attenuate the effects due to .
2.2. Nonlinear -Control Design
Let us consider the nonautonomous nonlinear systems of the form where is the performance output to be controlled and is the output of the system (16).
In the forthcoming analysis, the following assumptions are considered .
Assumption 3. The functions , , , , , and are assumed to be piecewise continuous in for all and continuously differentiable in for all .
Assumption 4. , , and for all .
Assumption 5. and .
Assumption 3 guarantees the well-posedness of the system (16), while being enforced by integrable exogenous inputs, and allows including nonsmooth nonlinearities. Assumption 4 ensures that the origin is an equilibrium point of the nondriven () disturbance-free () dynamic system (16). Assumption 5 is a simplified assumption inherited from the standard -control problem .
By considering the full-information case, the static state-feedback controller is said to be an admissible controller if equilibrium point of the closed-loop system (16)-(17) is asymptotically stable when . Besides, given the disturbance attenuation level , the system (16)-(17) is said to have -gain less than if the response , resulting from for initial state , satisfies for all and all piecewise-continuous functions . Thus, an admissible controller of the form (17) constitutes a solution of the -control problem if there exists a neighborhood of the origin such that inequality (18) is satisfied for all and all piecewise-continuous functions for which the state trajectory of the corresponding closed-loop system, starting from the initial point , remains in for all .
Let us consider the following hypothesis .
Hypothesis 1. For some positive and some positive definite function , there exists a locally Lipschitz continuous positive definite decrescent, radially unbounded solution of the Hamilton-Jacobi-Isaacs (HJI) inequality specified with
Theorem 6. Consider the system (16), with Assumptions 1–3. Let Hypothesis 1 be satisfied. Then, the static state-feedback controller is globally admisible and, with respect to a given output , the -gain of the closed-loop system (16), (21) is less than .
Summarizing, the following result is obtained.
Theorem 7. Consider system (6) with Assumptions 1–5 and Hypothesis 1. Let us assume that (7) holds. Then, the controller specified by the ISM controller together with the static state-feedback -control makes the equilibrium point of the unmatched disturbance-free closed-loop system (6), (10), and (21) asymptotically stable. Besides, the -gain of the perturbed nominal system (16) is less than with respect to output .
Proof. The conditions for the trajectories of (6) to converge to the manifold in (9) and the sliding mode to exist on this manifold may be derived based on the Lyapunov function whose time derivative along the solution of the closed-loop system (6) and (10), is given by Substituting (10) into the latter equation yieldsand, then, proves to be negative definite for any satisfying (11). The proof is completed by following the same line of reasoning of Theorem 22 from , which allows concluding asymptotic stability using -control for the full-information case employing system (16).
3. Application to a Flexible Joint Robot Manipulator
In order to support the applicability and performance of the proposed controller (22) given in Theorem 7, let us consider a 1-DOF FJR manipulator described by the following differential equations : where are the link position and velocity, respectively; , are the joint position and velocity, respectively; is the link mass; is the gravity constant; is the spring stiffness constant; is the rotor inertia; is the torque applied at the joint; and , are the matched and unmatched perturbations, respectively.
The tracking control problem for this system can be established as follows.
Control Problem. Given the smooth desired reference signal , find a control law such that, for any initial condition , the limit holds for the disturbance-free system (28)-(29), and the -gain of the perturbed nominal system (28)-(29) is less than with respect to a given performance output .
In order to apply the proposed controller stated in Theorem 7 to a FJR manipulator, let us consider the block diagram depicted in Figure 2. First, the dynamic equation (29) is decoupled from (28) using the controller . This decoupling controller combines the ISM control (10) and the -control (24) together with a properly defined reference signal and its first and second time derivatives. The reference signal is generated using the FJR manipulator output measurements and the outputs generated by and ISM controllers. In the following, the steps required to obtain the controller are described.
3.1. Stability Analysis
The first step to follow consists of designing a control law that decouples the dynamic equation (29) from (28). To this end, let us consider the control law where are positive constants, is the joint position error, is the joint velocity error, is the second time derivative of the joint reference signal , and is an auxiliary control input. Both and are to be defined later.
Let us consider the following change of coordinates: with .
Now, let us define the reference signal as follows:
Let us select the matrix in (9) as
From (36)-(37), it is clear that Assumptions 1 and 2 are satisfied; that is, and is nonsingular. Then, by assuming that satisfy inequalities (7), controller (10) drives the trajectories of system (35) towards the sliding surface given in (9). Then, the closed-loop system (35) and (10) can be expressed in the form (16), with given by
Since the full-information case is considered, the performance output to be controlled and the output are selected as with a positive constant. Then, the elements are given by
Now, let us define the smooth functions and as with , , being positive constants. Besides, the constants , are defined as Functions and are positive definite if the inequalities hold.
Using (38) and (42), the functions in (20) are given by and the left hand side of the HJI inequality (19), defined as , is given by with After completing squares, the expression (46) can be upper bounded as follows: and, then, will be negative semidefinite if fulfill the following inequalities: Then, the following theorem can be stated.
Theorem 8. Let us consider the system (28)-(29), with the reference signal defined by (34). Assume that there exist positive constants , , with , being defined by (43). If the inequalities (44) and (49) hold, controller (31), with , where is given by (10) and is given by accomplishes the control problem (30) for the disturbance-free system (28)-(29) and makes the -gain of the perturbed nominal system (28)-(29) less than with respect to the output given in (40).
3.2. Comparison with a Nonlinear Controller
The performance of the proposed controller was compared with that corresponding to the nonlinear controller developed by Cao and Xu , modified to deal also with matched and unmatched uncertainties. The next development details the procedure to implement the Cao-Xu controller in a 1-DOF FJR manipulator.
Under the same line of reasoning presented in Section 3.1, the dynamics of (52) are decoupled from (51) using the following control law: Substituting (53) into (52) yields which can be written in state-space form as
By a proper selection of the gains , , the origin of the system (55) is exponentially stable, that is, converges exponentially to .
Now, we will prove that if , this implies the exponential convergence of to the desired link reference signal . With this aim, let us define the desired rotor reference signal as follows:
A proper selection of the gains , allows concluding the exponential convergence of the origin . Thus, the controller (53), with defined as in (56), exponentially stabilizes the unperturbed nonlinear system (51)-(52), as was specified for the controller presented in .
Thus, the controller presented in  to be applied to the latter 1-DOF FJR manipulator is where is the nonlinear integral-type sliding surface defined as and is a constant row vector.
The switching surface must also fulfill the following inequality: provided , and .
4. Simulation Results
This section presents the numerical simulation results obtained using the proposed controller and the Cao-Xu controller. In the following the proposed method and the Cao-Xu controller will be referred to as the iSMH controller and the Cao-Xu controller, respectively. Simulations were performed using the software MATLAB-SIMULINK®, using the solver ode45 Dormand-Prince. The parameters of the model (28)-(29) were set to [kg], [Nmrad], , and , and all the initial conditions were set to zero.
The numerical simulations considered two trajectories to follow:(1)A sinusoidal desired link reference signal described by (2)A combination of sinusoids for the desired link reference signal; that is,
In order to provide a set of realistic gains for tuning the controllers, the parameters corresponding to an existing platform of Quanser  were considered. Specifically, the Maxon 273759 DC motor was taken as reference model, which has a torque constant [Nm/A]. Besides, it was assumed that the motor drives an harmonic gearbox with gear ratio of . The previous motor’s specifications allow obtaining the required voltage for each controller, in other words, the voltage demanded by each controller is obtained by using the following equation:
The controller gains for both, the iSMH and Cao-Xu controllers, was tuned to obtain a similar behavior for the link and rotor position errors. These gains fulfill the constraints inherent of each controller. Furthermore, it is assumed that the magnitude of the maximum allowable voltage was 10 [volts]. In all simulations, the matched perturbation was selected as a square signal with amplitude of [rad] and frequency of [rad/s]. The unmatched perturbation was chosen as a sinusoidal signal with unitary amplitude and frequency of [rad/s].
The link position error and joint position error for both controllers are depicted in Figure 3. Notice that the position errors are similar as was prespecified. However, it can be observed that a more oscillatory behavior is exhibited by the Cao-Xu controller.
In Figure 4, the behavior of the sliding surfaces , and the voltage demanded by each controller are shown. It is clear that the system trajectories are closer to the ideal sliding mode using the iSMH controller. Besides, the voltage demanded by the Cao-Xu controller is bigger than that required by the proposed controller.
The link position error and joint position error for the iSMH and Cao-Xu controllers are depicted in Figure 5. Again, the behavior of the position errors is similar for both controllers. However, it is observed again that a more oscillatory behavior is exhibited by the Cao-Xu controller.
In Figure 6, the behavior of the sliding surfaces , and the voltage demanded by each controller are shown. It is clear that the system is closer to the ideal sliding mode using the iSMH controller, and the voltage demanded by te Cao-Xu controller is slightly greater than that required by the proposed controller.
The previous simulation results show that, although similar responses may be obtained using any of the controllers, the iSMH requires less energy to achieve the control objective, while the oscillatory behavior is significantly less than that of the Cao-Xu controller. The reason for such behavior is that the proposed controller attenuates the effect of the unmatched uncertainty by means of the -controller. To illustrate this, Figure 7 shows the effect of the -control in the performance of the closed-loop system. The first seconds show the behavior of the joint position error, link position error, and the sliding surface when the -control is absent. After this time interval, the behavior of these signals is shown but now turning on the -control. It is observed that the vibration at the joint is considerably reduced by introducing the -control. The unmatched perturbation affects the performance of the ISM control. However, when the -control is applied, the effect of on the performance of the closed-loop system is reduced, which shows the advantage of using the proposed controller.
The previous numerical simulations validate the proposed methodology for controlling underactuated nonlinear time-varying systems affected by matched and unmatched perturbations. Besides, it was observed that it was convenient to use the combination of the ISM control with a -control because the later reduces vibrations and allows the ISM controller to have a better performance, even in presence of unmatched disturbances.
It is important to remark that in practice, a saturation function can also be used instead of the signum function given in (10). From sliding mode control theory analysis, signum function allows satisfying the robustness property and the reachability condition, that is, finite time convergence. However, the chattering phenomenon emerges, caused by fast switching in real applications. Instead, saturation functions are normally brought into play in order to avoid such undesirable phenomenon (see , e.g.).
In this paper, a controller for time-varying minimum-phase underactuated systems affected by matched and unmatched perturbations was presented. The proposed controller uses a combination of an ISM controller with -control. A Lyapunov approach was used to prove theoretically asymptotic stability of the equilibrium of the disturbance-free closed-loop system when applying the proposed controller and the boundedness of the solutions of the closed-loop system for the perturbed case. It also included an example using a flexible joint robot manipulator to show how the proposed controller can be applied to an underactuated system. In order to assess the performance of the proposed controller, an existing nonlinear controller based on a nonlinear integral-type sliding surface was implemented. Numerical simulations have shown that the proposed controller has the best performance despite the presence of matched and unmatched perturbations.
The authors declare that there is no conflict of interests regarding the publication of this paper.
This work was supported in part by CONACYT Project Cátedras 1537. C. Chavez would like to thank Programa para el Desarrollo Profesional Docente, Universidad Autónoma de Baja California, and the program of high quality scholarship for doctor studies.
- P. Ioannou and J. Sun, Robust Adaptive Control, Dover, Mineola, NJ, USA, 2012.
- R. Yedavalli, Robust Control of Uncertain Dynamic Systems: A Linear State Space Approach, Springer, London, UK, 2012.
- K. Zhou and J. Doyle, Essentials of Robust Control, Prentice Hall, New Jersey, NJ, USA, 1998.
- L. T. Aguilar, I. Boiko, L. Fridman, and R. Iriarte, Self-oscillations in dynamic systems, Systems & Control: Foundations & Applications, Springer International Publishing, Switzerland, Europe, 2015.
- V. I. Utkin, Sliding modes in control and optimization, Communications and Control Engineering Series, Springer, Berlin, Germany, 1992.
- Y. Orlov, Discontinuous Systems: Lyapunov Analysis and Robust Synthesis under Uncertainy Conditions, Springer, London, UK, 2014.
- P. R. Kumar, A. Chalanga, and B. Bandyopadhyay, “Smooth integral sliding mode controller for the position control of Stewart platform,” ISA Transactions, vol. 58, pp. 543–551, 2015.
- S. Mahieddine-Mahmoud, M. Ouriagli, L. Chrifi-Alaoui, and P. Bussy, “Robust control of nonlinear systems with both matched and unmatched disturbances,” in Proceedings of the 18th Mediterranean Conference on Control and Automation (MED '10), pp. 1043–1048, Marrakech, Morocco, June 2010.
- X. Han, E. Fridman, and S. K. Spurgeon, “Sliding-mode control of uncertain systems in the presence of unmatched disturbances with applications,” International Journal of Control, vol. 83, no. 12, pp. 2413–2426, 2010.
- L. Fridman, J. Barbot, and F. Plestan, Recent Trends in Sliding Mode Control, IET, London, UK, 2016.
- T. Osuna, O. E. Montano, and Y. Orlov, “Nonlinear -gain analysis of hybrid systems in the presence of sliding modes and impacts,” Mathematical Problems in Engineering, vol. 2016, Article ID 9074096, 10 pages, 2016.
- M. Rubagotti, A. Estrada, F. Castaños, A. Ferrara, and L. Fridman, “Integral sliding mode control for nonlinear systems with matched and unmatched perturbations,” IEEE Transactions on Automatic Control, vol. 56, no. 11, pp. 2699–2704, 2011.
- W.-J. Cao and J.-X. Xu, “Nonlinear integral-type sliding surface for both matched and unmatched uncertain systems,” IEEE Transactions on Automatic Control, vol. 49, no. 8, pp. 1355–1360, 2004.
- F. Castaños and L. Fridman, “Analysis and design of integral sliding manifolds for systems with unmatched perturbations,” IEEE Transactions on Automatic Control, vol. 51, no. 5, pp. 853–858, 2006.
- R. Galvan-Guerra and L. Fridman, “Robustification of time varying linear quadratic optimal control based on output integral sliding modes,” IET Control Theory and Applications, vol. 9, no. 4, pp. 563–572, 2015.
- B. C. Zheng and J. H. Park, “Adaptive integral sliding mode control with bounded gain performance of uncertain quantized control systems,” IET Control Theory and Applications, vol. 9, no. 15, pp. 2273–2282, 2015.
- C.-C. Chen, S. S.-D. Xu, and Y. Liang, “Study of nonlinear integral sliding mode fault-tolerant control,” IEEE/ASME Transactions on Mechatronics, vol. 21, no. 2, pp. 1160–1168, 2016.
- C. Chávez, L. Aguilar, and J. Mérida, “Analysis and synthesis of global nonlinear controller for robot manipulators,” Mathematical Problems in Engineering, vol. 2015, Article ID 410873, 9 pages, 2015.
- Y. V. Orlov and L. T. Aguilar, Advanced H∞ Control, Systems & Control: Foundations & Applications, Springer, New York, NY, USA, 2014.
- H. Khalil, Nonlinear Control, Pearson Education, Upper Saddle River, NJ, USA, 2015.
- L. Fridman, A. Poznyak, and F. J. Bejarano, Robust output LQ optimal control via integral sliding modes, Systems & Control: Foundations & Applications, Springer, New York, NY, USA, 2014.
- V. Utkin, J. Guldner, and J. Shi, Sliding Mode Control in Electromechanical Sytems, Automation and Control Engineering Series, CRC Press, 2nd edition, 2009.
- R. Kelly, B. Santibañez, and A. Loria, Control of Robot Manipulators in Joint Space, Springer, London, UK, 2005.
- S. He, D. Lin, and J. Wang, “Chattering-free adaptive fast convergent terminal sliding mode controllers for position tracking of robotic manipulators,” Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, vol. 230, no. 4, pp. 514–526, 2016.
Copyright © 2017 Roger Miranda-Colorado 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.