Nonlinear optimal control for robotic exoskeletons with electropneumatic actuators
Authors: Rigatos Gerasimos G. 
Published: May 2024

Purpose
To provide high torques needed to move a robot’s links, electric actuators are followed by a transmission system with a high transmission rate. For instance, gear ratios of 100:1 are often used in the joints of a robotic manipulator. This results into an actuator with large mechanical impedance (also known as nonback-drivable actuator). This in turn generates high contact forces when collision of the robotic mechanism occur and can cause humans’ injury. Another disadvantage of electric actuators is that they can exhibit overheating when constant torques have to be provided. Comparing to electric actuators, pneumatic actuators have promising properties for robotic applications, due to their low weight, simple mechanical design, low cost and good power-to-weight ratio. Electropneumatically actuated robots usually have better friction properties. Moreover, because of low mechanical impedance, pneumatic robots can provide moderate interaction forces which is important for robotic surgery and rehabilitation tasks. Pneumatic actuators are also well suited for exoskeleton robots. Actuation in exoskeletons should have a fast and accurate response. While electric motors come against high mechanical impedance and the risk of causing injuries, pneumatic actuators exhibit forces and torques which stay within moderate variation ranges. Besides, unlike direct current electric motors, pneumatic actuators have an improved weight-to-power ratio and avoid overheating problems.

Design/methodology/approach
The aim of this paper is to analyze a nonlinear optimal control method for electropneumatically actuated robots. A two-link robotic exoskeleton with electropneumatic actuators is considered as a case study. The associated nonlinear and multivariable state-space model is formulated and its differential flatness properties are proven. The dynamic model of the electropneumatic robot is linearized at each sampling instance with the use of first-order Taylor series expansion and through the computation of the associated Jacobian matrices. Within each sampling period, the time-varying linearization point is defined by the present value of the robot’s state vector and by the last sampled value of the control inputs vector. An H-infinity controller is designed for the linearized model of the robot aiming at solving the related optimal control problem under model uncertainties and external perturbations. An algebraic Riccati equation is solved at each time-step of the control method to obtain the stabilizing feedback gains of the H-infinity controller. Through Lyapunov stability analysis, it is proven that the robot’s control scheme satisfies the H-infinity tracking performance conditions which indicate the robustness properties of the control method. Moreover, global asymptotic stability is proven for the control loop. The method achieves fast convergence of the robot’s state variables to the associated reference trajectories, and despite strong nonlinearities in the robot’s dynamics, it keeps moderate the variations of the control inputs.

Findings
In this paper, a novel solution has been proposed for the nonlinear optimal control problem of robotic exoskeletons with electropneumatic actuators. As a case study, the dynamic model of a two-link lower-limb robotic exoskeleton with electropneumatic actuators has been considered. The dynamic model of this robotic system undergoes first approximate linearization at each iteration of the control algorithm around a temporary operating point. Within each sampling period, this linearization point is defined by the present value of the robot’s state vector and by the last sampled value of the control inputs vector. The linearization process relies on first-order Taylor series expansion and on the computation of the associated Jacobian matrices. The modeling error which is due to the truncation of higher-order terms from the Taylor series is considered to be a perturbation which is asymptotically compensated by the robustness of the control algorithm. To stabilize the dynamics of the electropneumatically actuated robot and to achieve precise tracking of reference setpoints, an H-infinity (optimal) feedback controller is designed. Actually, the proposed H-infinity controller for the model of the two-link electropneumatically actuated exoskeleton achieves the solution of the associated optimal control problem under model uncertainty and external disturbances. This controller implements a min-max differential game taking place between: (i) the control inputs which try to minimize a cost function which comprises a quadratic term of the state vector’s tracking error and (ii) the model uncertainty and perturbation inputs which try to maximize this cost function. To select the stabilizing feedback gains of this H-infinity controller, an algebraic Riccati equation is being repetitively solved at each time-step of the control method. The global stability properties of the H-infinity control scheme are proven through Lyapunov analysis.

Research limitations/implications
Pneumatic actuators are characterized by high nonlinearities which are due to air compressibility, thermodynamics and valves behavior and thus pneumatic robots require elaborated nonlinear control schemes to ensure their fast and precise positioning. Among the control methods which have been applied to pneumatic robots, one can distinguish differential geometric approaches (Lie algebra-based control, differential flatness theory-based control, nonlinear model predictive control [NMPC], sliding-mode control, backstepping control and multiple models-based fuzzy control). Treating nonlinearities and fault tolerance issues in the control problem of robotic manipulators with electropneumatic actuators has been a nontrivial task.

Practical implications
The novelty of the proposed control method is outlined as follows: preceding results on the use of H-infinity control to nonlinear dynamical systems were limited to the case of affine-in-the-input systems with drift-only dynamics. These results considered that the control inputs gain matrix is not dependent on the values of the system’s state vector. Moreover, in these approaches the linearization was performed around points of the desirable trajectory, whereas in the present paper’s control method the linearization points are related with the value of the state vector at each sampling instance as well as with the last sampled value of the control inputs vector. The Riccati equation which has been proposed for computing the feedback gains of the controller is novel, so is the presented global stability proof through Lyapunov analysis. This paper’s scientific contribution is summarized as follows: (i) the presented nonlinear optimal control method has improved or equally satisfactory performance when compared against other nonlinear control schemes that one can consider for the dynamic model of robots with electropneumatic actuators (such as Lie algebra-based control, differential flatness theory-based control, nonlinear model-based predictive control, sliding-mode control and backstepping control), (ii) it achieves fast and accurate tracking of all reference setpoints, (iii) despite strong nonlinearities in the dynamic model of the robot, it keeps moderate the variations of the control inputs and (iv) unlike the aforementioned alternative control approaches, this paper’s method is the only one that achieves solution of the optimal control problem for electropneumatic robots.

Social implications
The use of electropneumatic actuation in robots exhibits certain advantages. These can be the improved weight-to-power ratio, the lower mechanical impedance and the avoidance of overheating. At the same time, precise positioning and accurate execution of tasks by electropneumatic robots requires the application of elaborated nonlinear control methods. In this paper, a new nonlinear optimal control method has been developed for electropneumatically actuated robots and has been specifically applied to the dynamic model of a two-link robotic exoskeleton. The benefit from using this paper’s results in industrial and biomedical applications is apparent.

Originality/value
A comparison of the proposed nonlinear optimal (H-infinity) control method against other linear and nonlinear control schemes for electropneumatically actuated robots shows the following: (1) Unlike global linearization-based control approaches, such as Lie algebra-based control and differential flatness theory-based control, the optimal control approach does not rely on complicated transformations (diffeomorphisms) of the system’s state variables. Besides, the computed control inputs are applied directly on the initial nonlinear model of the electropneumatic robot and not on its linearized equivalent. The inverse transformations which are met in global linearization-based control are avoided and consequently one does not come against the related singularity problems. (2) Unlike model predictive control (MPC) and NMPC, the proposed control method is of proven global stability. It is known that MPC is a linear control approach that if applied to the nonlinear dynamics of the electropneumatic robot, the stability of the control loop will be lost. Besides, in NMPC the convergence of its iterative search for an optimum depends on initialization and parameter values selection and consequently the global stability of this control method cannot be always assured. (3) Unlike sliding-mode control and backstepping control, the proposed optimal control method does not require the state-space description of the system to be found in a specific form. About sliding-mode control, it is known that when the controlled system is not found in the input-output linearized form the definition of the sliding surface can be an intuitive procedure. About backstepping control, it is known that it cannot be directly applied to a dynamical system if the related state-space model is not found in the triangular (backstepping integral) form. (4) Unlike PID control, the proposed nonlinear optimal control method is of proven global stability, the selection of the controller’s parameters does not rely on a heuristic tuning procedure, and the stability of the control loop is assured in the case of changes of operating points. (5) Unlike multiple local models-based control, the nonlinear optimal control method uses only one linearization point and needs the solution of only one Riccati equation so as to compute the stabilizing feedback gains of the controller. Consequently, in terms of computation load the proposed control method for the electropneumatic actuator’s dynamics is much more efficient.