Journal of Robotics

Volume 2018 (2018), Article ID 9324623, 15 pages

https://doi.org/10.1155/2018/9324623

## A Family of Hyperbolic-Type Explicit Force Regulators with Active Velocity Damping for Robot Manipulators

^{1}Facultad de Ciencias de la Electrónica, Benemérita Universidad Autónoma de Puebla, 72570 Puebla, PUE, Mexico^{2}Departamento de Ingeniería Robótica, Universidad Autónoma de Aguascalientes, Campus Sur, Av. Prolongación Mahatma Gandhi No. 6601, 20340 Aguascalientes, AGS, Mexico^{3}Facultad de Ingeniería, Universidad Autónoma de San Luis Potosí, Av. Manuel Nava 8, Zona Universitaria, 78290 San Luis Potosí, SLP, Mexico

Correspondence should be addressed to Fernando Reyes-Cortés; xm.paub.oerroc@seyer.odnanref

Received 28 December 2017; Accepted 1 March 2018; Published 16 April 2018

Academic Editor: L. Fortuna

Copyright © 2018 Fernando Reyes-Cortés 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.

#### Abstract

This paper addresses the explicit force regulation problem for robot manipulators in interaction tasks. A new family of explicit force-control schemes is presented, which includes a term driven by a large class of saturated-type hyperbolic functions to handle the force error. Also, an active velocity damping term with the purpose of obtaining energy dissipation on the contact surface is incorporated plus compensation for gravity. In order to ensure asymptotic stability of the closed-loop system equilibrium point in Cartesian space, we propose a strict Lyapunov function. A force sensor placed at the end-effector of the robot manipulator is used in order to feed back the measure of the force error in the closed-loop, and an experimental comparison of the performance -norm between 5 explicit force control schemes, which are the classical proportional-derivative (PD), arctangent, and square-root controls and two members of the proposed control family, on a two-degree-of-freedom, direct-drive robot manipulator, is presented.

#### 1. Introduction

The use of robotic systems has increased during the last four decades. Nowadays, several industrial applications continue to benefit from robotic systems in order to automate repetitive or hazardous maneuvers. It is well known that, in applications involving an automatic process in which the robot is moving freely within its workspace or executes tasks with no interaction with the environment, a control system based on position or motion control represents a viable alternative. However, in applications where the manipulator performs interaction tasks, the manipulator encounters environmental constraints and the interaction forces are not negligible [1–3].

In this situation, the use of motion-control strategies will produce an unsuccessful execution of the task; the adequate implementation of an interaction task requires, in addition to motion/position control strategies, force control and a detailed and accurate modeling of the system and the environment [4–9]. A particular type of interaction control is the impedance control scheme which can regulate environmental forces according to the specified mechanical impedance such as inertial, stiffness, and damping effects. Impedance control belongs to the category of indirect force control, which may not require measurements of contact forces; in addition, the user cannot specify a desired contact force [2, 4, 9–12].

In contrast, force sensing and force control are essential for the successful execution of tasks involving interaction between the manipulator and the environment. The explicit force control strategies are the most suitable option for controlling interaction forces between the robot manipulator and the environment. Such control schemes have a large number of applications in medical and industrial environments, as they are capable of increasing the speed and quality in applications such as deburring, grinding of surfaces, assembly processes, and other operations [3, 10, 13–16].

Some types of force control can be cited, such as [17] where two auxiliary, automatically controlled wings of a monohull ship used to attenuate the roll due to the waves are presented. These wings are variably dipped into the water in order to introduce a force that reduces the oscillations on the ship. Two different compensators for the roll effect have been designed, the first one using classical frequency domain techniques and the second being an adaptive, linear quadratic compensator with a multilayer perceptron neural network. Simulations and experimental tests performed on the ship are presented. In [18], an impedance function is proposed to specify a desired force directly in a robust robot force tracking impedance scheme, which employs a neural network as compensator to cancel out all uncertainties; simulation studies are presented. In [19], the robot dynamics are decomposed into force, position, and redundant joint subspaces; then, using that decomposition in a neural network-based adaptive control scheme for hybrid force/position in Lyapunov sense is proposed. Numerical simulation studies were performed in the model of a two-link, rigid robot manipulator. In [20], an approach based on hybrid force/position with fuzzy PID structure with unspecified robot dynamics in the presence of external disturbances, including experimental results on a PUMA-560 robot, is presented.

Unlike the impedance control schemes, the concept of explicit force control is to consider a desired force rather than the position or velocity references. Then the force error is defined as the difference between a desired force and the actual interaction force, measured with a force sensor mounted on the end-effector, enabling the possibility of ensuring force regulation. The explicit force control selects the directions in which the end-effector position should be controlled and those directions in which the forces controlled by the end-effector are applied on some environment, for a given task. This type of control strategies has two feedback loops for separated, position, and force control and is referred to as hybrid force/position control [1, 21–34].

The original hybrid force/position control concept was proposed in [1]; the basic idea is to divide the robot manipulator space into two orthogonal spaces using of a selection matrix. Position control and force control are applied on two spaces to achieve the capability of compliance control. Explicit force control involves the direct command and measurement of force values in order to control the contact force to a desired value.

Considering that explicit force control is an important topic as research area, it currently represents an open problem of interest for the scientific community and several research efforts have focused on this topic during the last four decades. As a result, different control schemes have been proposed in the literature; for example, we can cite the fundamental works [1, 5, 24, 26, 32]. Progressively this theme has been growing; for instance, [35] describes the design and implementation of a hybrid force/motion scheme to control the interaction force of a robot manipulator, using a PD controller with gain-scheduled linear parameter-varying controller; however, although the reference for this control scheme does not include a stability proof, it presents experimental results on a six-degree-of-freedom manipulator. The paper [36] proposes the design of a robust adaptive neural network-based hybrid position/force control with PD-type structure for manipulators. It includes a stability analysis in the sense of Lyapunov, but only simulation results with a model of a two-link robot manipulator are presented. In [37], the robustness and stability of disturbance observer based on explicit linear force control systems are analyzed. The validity of the proposals is verified with experimental results on a linear dc motor.

In [38], explicit force regulators for robot manipulators are presented. The control structure includes a nonlinear term formed by square-root-type functions to drive the force error, with an active velocity damping plus compensation of gravity forces. The stability analysis is performed in Lyapunov sense; in order to obtain asymptotic stability, LaSalle’s theorem is used and experimental results on a two-degree-of-freedom robot are presented. In [39], an optimization-based approach for the regulation of excessive or insufficient forces at the end-effector level is introduced. Their approach minimizes the interaction force error at the robot end-effector and considers the linear impedance model and penalty functions in order to handle the constraints on the interaction forces; stability and convergence properties are analyzed taking into account force limits. The proposed scheme is validated via experimental results for a robotic grasping task. The work [40] presents the state of the art on compliant control algorithms applied to both stiff and soft joints based on electric motors for robot manipulators; there is a widely used control structure for all the works presented in this survey, which is the simple PD control.

Recently, in [16] a double active observer control is proposed to control the desired interaction and motion compensation under contact with the heart, relying on robot force control techniques. In [34], a position-based force explicit control strategy based on online trajectory prediction, with proportional-integral-type control structure, is presented. In this work, the experimental results show that the proposed method has a good control effect on the contact force. In [41], the regulation of the clamping force of a robot manipulator gripper is analyzed and discussed, when applied to rescue activities in the ruins after disasters; for this problem, a hybrid force/position control on a cascade PID-scheme and fuzzy sliding mode is presented; the control algorithms were validated with experimental results, whose force errors were acceptable. Force trajectories during isometric force-matching task involving contractions are described in [42]; a PID control with different gains and frequencies is introduced to examine the nature of the contractions in many individuals, which vary substantially. This variability can be explained by discrete PID with varying model parameters to adjust the muscle commands. The position and force control task presented in [43] has been defined considering the manipulator and environment models; asymptotic stability of the control systems is demonstrated, and experimental study of the issue is presented.

Most of the previous works use explicit force control strategies with PD or PID schemes; however, in general, there is a lack of proof on asymptotic stability analysis with different control structures; only a limited number of works have proposed strict Lyapunov functions. Therefore, our main motivation is in the theoretical and practical aspects of the problem. The objective of this paper is to present the design and analysis of a new explicit force control family with a term driven by a large class of saturated-type hyperbolic functions to handle the force error for the problem of force regulation in interaction tasks. This control structure includes also an active velocity damping term with the purpose of obtaining energy dissipation on the contact surface. To ensure global asymptotic stability of the closed-loop system equilibrium point in Cartesian space, we propose a strict Lyapunov function. Additionally, we present in this work another contribution, consisting of the mathematical development of a Cartesian dynamics model property (on the vector of centripetal and Coriolis forces in Cartesian space, see Appendix), which is useful for stability analysis.

The paper includes also an experimental comparison of the performance obtained between 5 explicit force control schemes; they are the classical proportional-derivative (PD), arctangent, square-root controls, and two members of the proposed control family on a two-degree-of-freedom, direct-drive robot manipulator. A force sensor placed at the end-effector of the robot manipulator is used in order to feed back the measure of the force error in the closed-loop, which is given as the difference between a desired and actual force. The Cartesian explicit force control is converted to joint torque signal using the transpose Jacobian introduced by Takegaki and Arimoto [44].

The paper is organized as follows. In Section 2, the joint-space dynamic model of robot manipulators and its relationship with the Cartesian dynamic model are presented; it also describes the most important properties used in the stability analysis of the control schemes. In Section 3, the new family of explicit force regulators is presented, as well as the design of the strict function required to carry out the asymptotic stability analysis in Lyapunov sense. In Section 4, the experimental results on a two-degree-of-freedom, direct-drive robot manipulator including a performance comparison of the control schemes under evaluation are presented. Finally, our conclusions are in Section 5.

#### 2. Dynamic Model of the Robot Manipulator

It is a well-known fact that the dynamic model of a robot manipulator with rigid links is described by a second-order system formed by nonlinear differential equations. This model describes a relationship between the generalized joint torques/forces applied by the actuators and the motion of the robot, which incorporates the position, velocity, and acceleration joints. In joint space, the dynamic model of a robot manipulator with degrees of freedom is defined by [32]where denote the joint displacements, velocities, and acceleration values, respectively, stands for the vector of generalized forces/torques applied, is the symmetric positive-definite inertia matrix, is the matrix of centripetal and Coriolis torques, is the vector of gravitational torques, and, finally, is the friction torque vector. For analysis purposes, only the viscous-friction model is considered in this paper; that is, , where is a diagonal matrix formed by the viscous-friction coefficients of each joint.

Moreover, the above model (1) can be expressed in terms of coordinates associated with a coordinate system attached to the tool by using a coordinate transformation [6, 25, 45]. For the formulation of the model in Cartesian space, it is assumed that the manipulator performs in a nonsingular configuration.

Let us consider the relationship between the joint velocities and the end-effector Cartesian velocities , given by , where is referred to as the analytical Jacobian matrix; that is, , where is the direct kinematic. A well-known relationship exists between the applied forces on the end-effector and the applied torques in the joints, which is referred to as the Jacobian Transposed Controller [44]:

The Cartesian dynamic model of a robot manipulator is [45]where the term represents the external forces applied to the end-effector, such as the forces appearing in an interaction task. The terms , and are described asWe use the notation to indicate the time derivative of the Jacobian matrix, , and the inverse of transpose Jacobian matrix is denoted by .

Next, some well-known properties that are useful in the analysis and design of our controller scheme are presented.

*Property 1. *The inertia matrix is symmetric and definite positive (on a workspace without singularity).

*Property 2. *The matrix is upper and lower bounded bywhere and stand for the maximum and minimum eigenvalues of matrix , respectively.

*Property 3. *; the matrix is skew-symmetric and it satisfies

*Property 4. *The norm of vector is bounded, as indicated below:for some . As part of the contributions of this article, a brief demonstration of this property is included in the Appendix.

The following decoupled model is considered for modeling the compliant environment [3]:where is a diagonal matrix referred to as the environment’s stiffness matrix and the term denotes the contact position with the environment, while represents the actual position of the end-effector.

#### 3. A New Family of Explicit Force Regulator Schemes

In this section, a new family of force control schemes with hyperbolic-type structure for robot manipulators is presented. Consider the following family with large class of explicit force control schemes plus gravity compensation, given bywhere is a positive integer, are the proportional and derivative gains, which are positive-definite diagonal matrices, and is the force error vector, which is defined as the difference between a desired force and the actual force on the end-effector, measured from the force sensor; . The applied torques on the robot joints are derived from (2); that is, .

The first term in the proposed control scheme (10) represents a large class of hyperbolic-type regulators to drive the force error; this one is derived from energy shaping of the artificial potential energy [44]. The second term acts as an active velocity damping; the next component includes the Cartesian gravity forces ; for implementation purposes, partial knowledge of these gravity forces is required, and the last term represents the interaction forces.

The closed-loop equation expressed in state variables is obtained by combining the dynamics model (3) and control law (10):In order to demonstrate that the state origin is an equilibrium point, it is necessary to consider the Jacobian matrix with full rank; for example, ; in other words, the Cartesian dynamic model requires a singularity-free workspace. The following assumptions are taken into account:(a)The Jacobian matrix is full rank and its spectral norm is bounded for all , where is a positive constant, .(b)According to the environment’s stiffness and friction models, the stiffness and viscous-friction matrices are diagonal and positive-definite matrices, respectively; then .

The first component of the closed-loop equation (11) satisfies , due to the fact that the matrix is a diagonal definite positive matrix. For the second component in (11), we use Property 1 of the inertia matrix; then the first term inside the brackets composed of the proportional gain and the vector of hyperbolic functions results in the vector , since is a positive-definite matrix and each component , and . We have also if and in the singularity-free workspace. Therefore, the closed-loop system (11) is autonomous differential equation and the state space origin exists and it is a local equilibrium point.

##### 3.1. Control Problem Statement

The problem of explicit force control consists of finding a Cartesian control scheme such that the force error vector and Cartesian velocity converge asymptotically stably to equilibrium point of the closed-loop equation (11) at least for sufficiently small initial conditions ; then the following objective is satisfied:

##### 3.2. Proposition

Consider the closed-loop system (11), under conditions sufficiently small; then the equilibrium point is locally asymptotically stable and the control objective (12) is achieved.

*Proof. *To carry out the local stability analysis for the dynamic system (11), we use Lyapunov’s direct method through the following strict Lyapunov function (a strict Lyapunov function is a positive-definite function, whose time derivative along the trajectories of closed-loop system yields a negative-definite function), which is composed by the sum of the robot manipulator kinetic and artificial potential energy, and a cross-term between force error and velocity :where is a positive constant number defined in the following interval:where is a positive constant, is the minimum eigenvalue of the product , the minimum eigenvalue of the environment’s inverse stiffness matrix is , and are the maximum eigenvalues of the proportional and derivative gains and viscous-friction matrix, respectively. It is important to note that the constant is only required for analysis purposes and its numerical value is not required for the control scheme.

The artificial potential energy function is a continuously differentiable and positive-definite function, given byThe candidate Lyapunov function (13) corresponds to the total energy of the closed-loop (11) and can be rewritten asNow, we will give sufficient conditions to make a positive-definite function. The first term of the right side of (16) is a definite positive function with respect to and , because inertia matrix is definite positive. Next, we provide upper and lower bounds on the second and third terms of the Lyapunov function (16).

The artificial energy function (15) has a unique and global minimum at , and it is not a decreasing function; then there exists positive constant , such that . Therefore, the second and third terms of function (16) satisfyOn the left-hand side of the inequality (17), to ensure that is a positive expression, must be in the interval , which is satisfied by the first term of (14). Hence, the Lyapunov function given by (13) is a definite positive function in local form.

The total derivative of the candidate Lyapunov function (13) with respect to time isSubstituting the closed-loop equation (11), we also use Property 3; then it yields,

where, after cancelling some terms, we obtainNote that the second term of the right side in (20) corresponds to the gradient of the artificial potential energy ; it is a vector composed of saturated functions with hyperbolic-type structure, for each th component = ; note also that ; then ; there exist positive constants , such that ; therefore it satisfiesIn other words, the second term of function (20) satisfies the following condition:where and stand for the maximum and minimum eigenvalues of matrix , respectively. Using Property 4, the fourth term from (23) satisfies . Next, by taking (22) and Cartesian dynamic model property (2), the expression for (20) takes the formThe time derivative of the Lyapunov function becomesIn order to ensure that matrix is positive definite, the number must satisfy the second condition of the left side of (14); then the element ; also the third condition of the left side in the expression (14) must be satisfied; then the determinant of matrix is positive. Considering that satisfies those considerations, we can conclude that is a negative-definite function; note that, due to the fact that both conditions and are satisfied, the proposed Lyapunov function is a strict Lyapunov function. Therefore, according to Lyapunov’s direct method, we conclude local asymptotic stability of the origin of the closed-loop equation (11), which means that both state variables and asymptotically converge to zero, as .

#### 4. Experimental Results

The experimental results were obtained by using an experimental platform known as “Rotradi” [46], shown in Figure 1. This platform, designed and built at the Robotics Laboratory of Benemérita Universidad Autónoma de Puebla (BUAP) in Mexico, is a two-degree-of-freedom, direct-drive robot. Rotradi consists of two 6061 aluminum links, actuated by brushless, direct-drive servo actuators, in order to drive the joints without a gear reduction. The motors used in Rotradi are DM-1150A and DM-1015B models from Parker Compumotor, for the shoulder and elbow joints, respectively. Since the servomotors are operated in torque mode, they act as a torque source; in this mode, an analog voltage is used as a reference of the torque signal, in proportional scale. The features of the servomotors are shown in Table 1. A motion-control board of Precision MicroDynamics Inc. is included in the robot system; this device is designed for reading the encoders and generating reference voltages. The control algorithms, written in C code, run in real time with a 2.5 ms sample period, on a Pentium computer.