# Collision-Free Compliance Control for Redundant Manipulators: An Optimization Case

^{1}Guangdong Key Laboratory of Modern Control Technology, Guangdong Institute of Intelligence Manufacturing, Guangzhou, China^{2}School of Engineering, Swansea University, Swansea, United Kingdom

Force control of manipulators could enhance compliance and execution capabilities, and has become a key issue in the field of robotic control. However, it is challenging for redundant manipulators, especially when there exist risks of collisions. In this paper, we propose a collision-free compliance control strategy based on recurrent neural networks. Inspired by impedance control, the position-force control task is rebuilt as a reference command of task-space velocities, by combing kinematic properties, the compliance controller is then described as an equality constraint in joint velocity level. As to collision avoidance strategy, both robot and obstacles are approximately described as two sets of key points, and the distances between those points are used to scale the feasible workspace. In order to save unnecessary energy consumption while reducing impact of possible collisions, the secondary task is chosen to minimize joint velocities. Then a RNN with provable convergence is established to solve the constraint-optimization problem in realtime. Numerical results validate the effectiveness of the proposed controller.

## 1. Introduction

Industry 4.0 is becoming a label of modern industry combining traditional manufacturing and increasingly technological world. As an important executor, robot manipulator must be more flexible and intelligent, to satisfy production requirements which is more personalized and customized (Gonzalez et al., 2018). Among various kinds of robot manipulators, redundant manipulators have become an important branch of robotics due to its flexibility (Zhang, 2015). This enables robots to fulfill more complicated tasks and has been a hot topic in the field of robotic control.

With the increasing popularity of robot manipulators, traditional position control based applications (such as welding, painting and so on) can hardly meet complex production tasks (He et al., 2015), for instance, in pure position control based structures, the interaction between robot and workpieces is usually ignored, which could probably lead to high security risk, since excessive system stiffness would lead to the unpredictable responses (Cai and Xiang, 2018). Therefore, aiming at enhancing the execution ability of the system, precise control of contact force is required to ensure compliance to external environment. Accordingly a series of control methods are proposed, depending on different robotic structure and control signals. By imitating flexible joints and muscles of animals, compliance units are introduced into the robots, such as series elastic actuators (SEA), variable stiffness actuators, etc. In Pan et al. (2018b), a compliance controller is designed for SEA based systems, and a modified command-filtered back-stepping control strategy (CFBC) based on adaptive mechanism is then proposed to overcome the discontinuous friction and complexity problem of traditional back-stepping based methods. By adjusting the compliance of joint angles, precise control of torque output is realized. As to the interaction between the robot and workpieces, Hogan proposes a basic idea of impedance control, in which the robot and environment usually bear as an impedance and admittance, respectively (Hogan, 1985). Generally speaking, the contact force and relative movement of the robot and workpieces can be described as a combination of mass-spring-damper systems. Therefore, the contact force can be controlled by designing motion commands indirectly. Another representative approach is hybrid position-force control, the controller is usually designed in the torque loop of the joint space, in which both contact forces and movement of the robot are modeled based on dynamic analysis. Then the controller can be described as a combination of control efforts which achieve position and force control, respectively (Raibert and Craig, 1981). Similar research can be found in literature such as (Khatib, 1987; Pan et al., 2018a, 2019; Zhao et al., 2018a,b).

During the operating process, since the manipulators are usually required to keep in touch with the workpieces, it is possible that the robot would collide with the environment. Besides, the workspace of a robot as also limited (Khatib, 1986). For example, in a production line with multiple manipulators, each robot is located at a fixed position, in order to avoid interference, the robot's workspace is limited by hardware (fences, barriers, etc.) or software constraints(pre-planned space). In situations such as human-machine collaboration, the robot must not collide with human. Therefore, it is crucial to avoid obstacles during the operating process. In present reports, the desired trajectory is generally obtained by off-line programming, which is limited by programming efficiency. To realize obstacle avoidance control in realtime, artificial potential field based methods are widely used. The basic idea of is that the target bears as an attractive pole while the obstacle creates repulsion on the robot, then the robot will be controlled to converge to the target without colliding with obstacles (Wang et al., 2018). In Csiszar et al. (2011), a modified method is proposed, which describes the obstacles by different geometrical forms, both theoretical conduction and experimental tests validate the proposed method. Considering the local minimum problem that may caused by multi-link structures, in Badawy (2016), a two minima is introduced to construct potential field, such that a dual attraction between links enables faster maneuvers comparing with traditional methods. Other improvements to artificial potential field method can be found in Tsai et al. (2001), Tsuji et al. (2002). A series of pseudo-inverse methods are constructed for redundant manipulators in Sciavicco and Siciliano (1988), in which the control efforts consists of a minimum-norm particular solution and homogeneous solutions, and the collision can be avoided by calculating a escape velocity as homogeneous solutions. By understanding the limited workspace, the obstacle avoidance can be described in forms of inequalities, which opens a new way in realtime collision avoidance. In Zhang and Wang (2004), the robot is regarded as the sum of several links, and the distances between the robot and obstacle is obtained by calculating distances between points and links. Then Guo and Zhang (2012) improves the method by modifying obstacle avoidance MVN scheme, and simulation results show that the modified control strategy can suppress the discontinuity of angular velocities effectively.

In terms with compliance control problem of a robot, the controller efforts should be designed according to the desired commands and system characteristics. That is so say, the robot must follow a constraint that achieves compliance control, and at the same time, the inequality constraints are ensured to avoid obstacles. It is obvious that the control problem involves several constraints, including equality constraints and inequality ones. Using the thought of constraint-optimization, the control problem with multiple constraints can be well handled. Recently, the applications of recurrent neural networks for robotic control have been studied extensively, and have shown great efficient for real-time processing (Wang et al., 2015; Jin et al., 2017; Xu et al., 2019a). In those literatures, analysis in dual space and a convex projection are introduced to handle inequality constraints.

Recently, taking advantage of parallel computing, neural networks are used to solve the constraint-optimization, and have shown great efficiency in real-time processing. In Zhang et al. (2004), Li et al. (2017), Yang et al. (2018b), controllers are established in joint velocity/acceleration level, to fulfill kinematic tracking problem for robot manipulators. In Xu et al. (2019b), tracking problem with model uncertainties is considered, and an adaptive RNN based controller is proposed for a 6DOF robot Jaco_{2}. Discussions on multiple robot systems, parallel manipulators, time-delay systems using RNN can be found in Zhang et al. (2018), Li et al. (2019), Xu et al. (2019b).

From the previous observations, we propose a RNN based collision-free compliance control strategy for redundant manipulators. The remainder of this paper is organized as follows. In section 2, the control objective including the position-force control as well as collision avoidance is pointed out, and then rewritten as a QP problem. In section 3, the RNN based controller is proposed, and the stability of the system is also analyzed. A number of numerical experiments on a 4-DOF redundant manipulator including model uncertainties and narrow workspace are carried out to further verify the proposed control strategy. section 5 concludes the paper. The contributions of this paper are summarized as below

• To the best of the author's knowledge, there is few research on compliance control using recurrent neural networks, the study in this paper is of great significance in enriching the theoretical frame of RNN.

• The proposed controller is capable of handling compliance control, as well as avoiding obstacles in realtime, which does make sense in industrial applications, besides, physical constraints are also guaranteed.

• Comparing to traditional neural-network-based controllers used in robotics, not only control errors but model information is considered, therefore, the proposed RNN has a simple structure, and the global convergence can be ensured.

## 2. Problem Formulation

### 2.1. Robot Kinematics and Impedance Control

Without loss of generality, we consider series robot manipulators with redundant DOFs, and the joints are assumed as rotational joints. Let θ ∈ ℝ^{n} be the vector of joint angles, the description of the end-effector in the cartesian space is:

where *x* ∈ ℝ^{m} is the coordination of the end-effector. In the velocity level, the forward kinematic model can be formulated as:

in which *J*(θ) = ∂*x*/∂θ is Jacobian matrix. As to redundant manipulators, *J* ∈ ℝ^{m × n}, *rank*(*J*) < *n*.

In industrial applications, position control based operation mode has many limitations: due to the lack of compliance, pure kinematic control methods may cause unexpected consequences, especially when the robot is in contact with the environment. To enhance the compliance and achieve precise control of contact force, according to impedance control technology, the interaction between robot and environment can be described as a damper-spring system, as shown in Figure 1 (Senoo et al., 2017).

where, *K*_{p} and *K*_{d} are interaction coefficients, and Δ*x* = *x* − *x*_{d} is the difference between the actual response *x* and desired trajectory *x*_{d}. The basic idea of impedance control methods is shown in Equation (2.1). By referring to Equations (2) and (3), we have:

When the real values of *K*_{p} and *K*_{d} are known, *F* can be obtained by adjusting the velocity ẋ of the end-effector according to Equation (4).

### 2.2. Obstacle Avoidance Scheme

In the process of robot force control, there is a risk of collision as the robot may contact with workpieces. Besides, robot manipulators usually work in a limited workspace restricted by fences, which are used to isolated robots from humans or other robots. This problem could be even more acute in tasks which requires collaboration of multiple robots. Therefore, obstacle avoidance problem must be taken into consideration. When collision does not happens, the distance between robot and obstacles keep positive. By describing the robot and obstacles as two separated sets, namely *S*_{A} = {*A*_{1}, ⋯ , *A*_{a}}, *S*_{B} = {*B*_{1}, ⋯ , *B*_{b}}, where *A*_{i}, *i* = 1, ⋯ , *a* and *B*_{j}, *j* = 1, ⋯ , *b* are points on the robot and obstacles, respectively. Then the sufficient and necessary conditions of obstacle avoidance problem is that the intersection of A and B is an empty set. That is to say, for any point pair *A*_{i} on the robot and *B*_{j} on the obstacle, the distance between *A*_{i} and *B*_{j} is always positive, i.e., $||{A}_{i}{B}_{j}{||}_{2}^{2}>0$, for all *i* = 1, ⋯ , *a*, *j* = 1, ⋯ , *b*, where $||\u2022{||}_{2}^{2}$ is the Euclidean norm of vector *A*_{i}*B*_{j}. For sake of safety, let *d* > 0 be a proper value describing the minimum distance between robot and obstacles, the collision can be avoided b ensuring $||{A}_{i}{B}_{j}{||}_{2}^{2}\ge d$.

*Remark 1*. In fact, both *S*_{A} and *S*_{B} consist of infinite points. However, by evenly selecting representative points from the robot link and obstacles, *S*_{A} and *S*_{B} can be simplified significantly. Besides, the safety distance *d* can be appropriately increased. Despite that this treatment will sacrifice some workspace of the robot (the inequality $||{A}_{i}{B}_{j}{||}_{2}^{2}\ge d$ would into account some areas that collisions do not happen, due to a bigger *d* is considered), this sacrifice is meaningful: the number of inequality constraints can be reduced greatly, which is helpful for constraint description and solution.

In real applications, the key points of the robot manipulator is easy to select. Cylindrical envelopes are usually used to describe the robotic links, then the key points can be selected on the axes of the cylinders uniformly, and the distance between those points can be defined the same as the radius of the cylinder. As to the obstacles with irregular shapes, the key points can be selected based on image processing techniques, such as edge detection, corrosion, etc.

### 2.3. Problem Reformulation in QP Type

From the above description, the purpose of this paper is to build a collision-free force controller for redundant manipulators, to achieve precise force control along a predefined trajectory, in the sense that *F* → *F*_{d}, *x* → *x*_{d}, and $||{A}_{i}{B}_{j}{||}_{2}^{2}\ge d$ for all *i* = 1, ⋯ , *a*, *j* = 1, ⋯ , *b*.

As to a redundant manipulator, there exist redundant DOFs, which can be used to enhance the flexibility of the robot. When the robot gets close to the obstacles, the robot must avoid the obstacle without affecting the contact force and tracking errors. On the other hand, when there is no risk of collision, the robot may work in an economic way, by minimizing the joint velocities, energy consumption can be reduced effectively. Therefore, by defining an objective function as $||\stackrel{.}{\theta}{||}_{2}^{2}$, the control objective can be summarized as:

where $||\stackrel{.}{\theta}{||}_{2}^{2}$ is the Euclidean norm of $\stackrel{.}{\theta}$. It is noteworthy that in actual industrial applications, the robot is also limited by its own physical structures. For instance, the joint angles are limited in a fixed range, and the upper/lower bounds of joint velocities are also constrained due to actuator saturation. By combing (Equation 4), the control objective rewrites to:

with θ^{−}, θ^{+}, ${\stackrel{.}{\theta}}^{-}$, ${\stackrel{.}{\theta}}^{+}$ being the upper/lower bounds of joint angles and velocities, respectively. However, the optimization problem is described in different levels, i.e., joint speed level or joint angle level, which remains challenging to solve (Equation 6) directly. Therefore, we will rewrite this formula in velocity level. As to the key points *A*_{i} on the robot, let *x*_{Ai} be the coordination of *A*_{i} in the cartesian space, both *x*_{Ai} and ẋ_{Ai} are available:

where *f*_{Ai}(•) is the forward kinematics of point *A*_{i}, and *J*_{Ai} is the corresponding Jacobian matrix from *A*_{i} to joint space. Let us consider the following equality:

in which *k* is a positive constant. It is obviously that the equilibrium point of Equation (8) is $||{A}_{i}{B}_{j}{||}_{2}^{2}=d$. By letting $\frac{\text{d}}{\text{d}t}(||{A}_{i}{B}_{j}{||}_{2}^{2})\ge 0$, the inequality (5d) can be readily guaranteed. Taking the time-derivative of $||{A}_{i}{B}_{j}{||}_{2}^{2}$ yields:

where, $\stackrel{\u20d7}{|{B}_{j}{A}_{i}|}={({A}_{i}-{B}_{j})}^{\text{T}}/||\stackrel{.}{\theta}{||}_{2}^{2}$ is a unit vector from *B*_{j} to *A*_{i}, and Ḃ_{j} is the velocity of key point *B*_{j} on the obstacles. By Equations (9) and (6c), the inequality description of obstacle avoidance strategy is

*Remark 2*. In this part, we have shown the basic idea of obstacle avoidance scheme in velocity level, whose equilibrium point is described in Equation (8). It is notable that the right-hand side of Equation (8) is only a common form to realize obstacle avoidance. Generally speaking, the right-hand side of Equation (8) may have different forms, such as $k(||{A}_{i}{B}_{j}{||}_{2}^{2}-d)$, $k{(||{A}_{i}{B}_{j}{||}_{2}^{2}-d)}^{3}$, etc. From Equation (10), the value of the response velocity to avoid obstacles is related to the two parts, the first part is the difference between the actual and safety distance, the other part depends on the movement of the obstacles.

In terms of the physical constraints of joint angles, according to escape velocity method, inequalities (6d) and (6e) can be uniformly described as $\mathrm{\text{max}}(\alpha ({\theta}^{-}-\theta ),{\stackrel{.}{\theta}}^{-})\le \stackrel{.}{\theta}\le \mathrm{\text{min}}({\stackrel{.}{\theta}}^{+},\alpha ({\theta}^{+}-\theta ))$. So far, the position-force control problem together with obstacle avoidance strategy in velocity level is as below

where (11c) is a rewritten inequality considering (6d) and (6e) based on escape velocity scheme (Zhang et al., 2004), ${J}_{o}=[\underset{b}{\underbrace{{\stackrel{\u20d7}{|{B}_{1}{A}_{1}|}}^{\text{T}}{J}_{A1};\cdots \phantom{\rule{0.3em}{0ex}};{\stackrel{\u20d7}{|{B}_{b}{A}_{1}|}}^{\text{T}}{J}_{Ab}}},\cdots \phantom{\rule{0.3em}{0ex}},\underset{b}{\underbrace{{\stackrel{\u20d7}{|{B}_{1}{A}_{a}|}}^{\text{T}}{J}_{Aa}^{\text{T}};\cdots \phantom{\rule{0.3em}{0ex}};{\stackrel{\u20d7}{|{B}_{b}{A}_{a}|}}^{\text{T}}{J}_{Ab}}}]$ ∈ ℝ^{ab × n} is the concatenated form of *J*_{Ai} considering all pairs between *A*_{i} and *B*_{j}, $B={\left[{B}_{11},\cdots \phantom{\rule{0.3em}{0ex}},{B}_{1b},\cdots \phantom{\rule{0.3em}{0ex}},{B}_{a1},\cdots \phantom{\rule{0.3em}{0ex}},{B}_{ab}\right]}^{\text{T}}\in {\mathbb{R}}^{ab}$ is the vector of upper-bounds, in which $-k(||{A}_{i}{B}_{j}{||}_{2}^{2}-d)-{\stackrel{\u20d7}{|{B}_{j}{A}_{i}|}}^{\text{T}}{\u1e02}_{j}$. From the definition of *J*_{o}, *B*, inequality (11d) in equivalent to ${\stackrel{\u20d7}{|{B}_{1}{A}_{1}|}}^{\text{T}}{J}_{A1}(\theta )\stackrel{.}{\theta}\ge k(||{A}_{1}{B}_{1}{||}_{2}^{2}-d)+{\stackrel{\u20d7}{|{B}_{1}{A}_{1}|}}^{\text{T}}{\u1e02}_{1}$,…${\stackrel{\u20d7}{|{B}_{b}{A}_{a}|}}^{\text{T}}{J}_{Aa}(\theta )\stackrel{.}{\theta}\ge k(||{A}_{a}{B}_{b}{||}_{2}^{2}-d)+{\stackrel{\u20d7}{|{B}_{b}{A}_{a}|}}^{\text{T}}{\u1e02}_{b}$, which is the cascading form of the inequality description (10) for all points pairs *A*_{i}*B*_{j}, i.e., if (11d) hold, the obstacle avoidance can be achieved. It is notable that a larger number of key points do help to describe the information of the obstacle more clearly, but it would lead to a computational burden, since the number of inequality constraints also increases. Therefore, the distance of the key points on the obstacle can be selected similar to those of the manipulator.

## 3. RNN Based Controller Design

In section II, we have transform the compliance control as well as obstacle avoidance problem into a constraint-optimization one. However, because that the QP problem described in Equation (11) contains equality and inequality constraints, moreover, both (Equations 11b,d) are nonlinear, it is difficult to solve directly, especially in industrial applications in realtime. Based on the parallel computation ability, a RNN is established to solve (Equation 11) online, and the stability of the closed-loop system is also discussed.

### 3.1. RNN Design

In terms with the QP problem (Equation 11), although the analytical solution can be hardly obtained, by defining a Lagrange function as:

where λ_{1} and λ_{2} are state variables, respectively. According to Karush-Kuhn-Tucker (KKT) conditions, the inherent solution of Equation (11) satisfies:

where, *P*_{Ω}(*x*) = argmin_{y ∈ Ω}||*y*−*x*|| is a projection operator of $\stackrel{.}{\theta}$ to convex Ω, and $\Omega =\left\{\stackrel{.}{\theta}\in {\mathbb{R}}^{n}|\mathrm{\text{max}}(\alpha ({\theta}^{-}-\theta ),{\stackrel{.}{\theta}}^{-})\le \stackrel{.}{\theta}\le \mathrm{\text{min}}({\stackrel{.}{\theta}}^{+},\alpha ({\theta}^{+}-\theta ))\right\}$. In Equation (13c), the operation function (•)^{+} is defined as a mapping to the non-negative space. Equation (13c) can be rewritten as:

When ${J}_{o}\stackrel{.}{\theta}\le B$, the inequality (Equation 11d) holds, then λ_{2} stays zero. Instead, if the inequality reaches a critical state, λ_{2} becomes positive to ensure ${J}_{o}\stackrel{.}{\theta}=B$. In order to obtain the inherent solution in real time, a recurrent neural network is built as follows:

with ϵ being a positive constant scaling the convergence of Equation (15).

The proposed RNN based algorithm is shown in Algorithm 1. Based on escape velocity method, the convex set of joint speed can be obtained based on the positive constant α and physical constraints θ^{−}, θ^{+}, ${\stackrel{.}{\theta}}^{-}$, ${\stackrel{.}{\theta}}^{-}$. After initializing state variables λ_{1} and λ_{2}, the reference velocity can be obtained based on the desired command and actual responses according to Equation (4). then the output of RNN (which is also the control command) can be calculated based on Equation (15a), at the same time, both λ_{1} and λ_{2} can be updated according to Equations (15b) and (15c).

In real applications, the nonlinear system can be hardly approximated completely. Therefore, the approximate error is inevitable, which would influence the performance of the proposed controller. However, the approximate error is a small value of higher order, and the influence can be suppressed based on the negative feedback scheme in the outer-loop, as shown in Equation (4).

*Remark 3*. The output dynamics of the proposed RNN is given in Equation (15a), in which the projection operator *P*_{Ω}(•) plays an important rule in handling physical constraints (Equation 11c), the updating of $\stackrel{.}{\theta}$ depends on three parts: the first part $-\stackrel{.}{\theta}/||\stackrel{.}{\theta}{||}_{2}^{2}$ in used to optimize the objective function $||\stackrel{.}{\theta}{||}_{2}^{2}$, and the second item ${J}^{\text{T}}{\lambda}_{1}$ guarantees the equality constraint (Equation 11b) by adjusting the dual state variable λ_{1} according to Equation (15b), and the last item $-{J}_{o}^{\text{T}}{\lambda}_{2}$ ensures the inequality constraint (Equation 11d). The RNN consists of three kinds of nodes, namely, $\ddot{\theta}$, λ_{1} and λ_{2}, with the number of neurons being *n* + *ab* + *m*.

It is remarkable that the proposed controller is based on the information of system models such as *J*, *J*_{o}, *K*_{p}, etc., which is helpful to reduce computational cost. As to the constraint-optimization problem (Equation 11), the main challenge is to solve it in real-time, since the parameters in constraints (Equations 11b, 11d) are time varying. From Equation (15), the control effort is obtained by calculating its updating law, which is based on the historical data and model information, i.e., it is no longer necessary to solve the solution of Equation (11) as every step, and the computational cost is thus reduced. In the following section, we will also show the convergence of the RNN based controller.

In this paper, we mainly concern the obstacle avoidance problem in force control tasks. It is notable that force control is mainly based on the idea of impedance control theory, which is similar to existing methods in Huang et al. (2019), Zhang and Xia (2019). The main challenge of the proposed control scheme lies in the limitation of sampling ability of cameras, which are used to capture the obstacles. To handle the measurement noise or disturbances, a larger safety distance *d* can be introduced to ensure the performance of obstacle avoidance.

### 3.2. Stability Analysis

*Lemma 1:* (Convergence for a class of neural networks) (Gao, 2003) A dynamic neural network is said to converge to its equilibrium point if it satisfies:

where κ > 0 and ϱ > 0 are constant parameters, and *P*_{S} = argmin_{y∈S}||*y* − *x*|| is a projection operator to closed set *S*.

*Definition 1:* For a given function *F*(•) which is continuously differentiable, with its gradient defined as ∇*F*, if ∇*F* + ∇*F*^{T} is positive semi-definite, *F*(•) is called a monotone function.

About the stability of the closed-loop system, we offer the following theorem.

*Theorem 1:* Given the collision-free position-force controller based on a recurrent neural network, the RNN will converge to the inherent solution (optimal solution) of Equation (11), and the stability of the closed-loop system is also ensured.

*Proof:* Define a vector ξ as $\xi =\left[\stackrel{.}{\theta};{\lambda}_{1};{\lambda}_{2}\right]\in {\mathbb{R}}^{n+m+ab}$, according to Equation (15), the time derivative of ξ satisfies:

in which ϵ > 0, and $F(\xi )={\left[{F}_{1}(\xi ),{F}_{2}(\xi ),{F}_{3}(\xi )\right]}^{\text{T}}$, where ${F}_{1}=\stackrel{.}{\theta}/||\stackrel{.}{\theta}{||}_{2}^{2}-{J}^{\text{T}}{\lambda}_{1}+{J}_{o}^{\text{T}}{\lambda}_{2}$, ${F}_{2}=J\stackrel{.}{\theta}-{K}_{d}^{-1}F+{K}_{p}{K}_{d}^{-1}\Delta x-{\u1e8b}_{\text{d}}$, ${F}_{3}=-{J}_{o}\stackrel{.}{\theta}+B$. By calculating the gradient of *F*(ξ), we have:

It is obviously that ∇*F*(ξ) is positive definite. According to Definition 1, *F*(ξ) is a monotone function. From the description of (17), the projection operator *P*_{S} can be formulated as *P*_{S} = [*P*_{Ω}; *P*_{R}; *P*_{Λ}], in which *P*_{Ω} is defined in (13a), *P*_{R} can be regarded as a projection operator of λ_{1} to *R*, with the upper and lower bounds being ±∞, and ${P}_{\Lambda}={(\u2022)}^{+}$ is a special projection operator to closed set ${\mathbb{R}}_{+}^{ab}$. Therefore, *P*_{S} is a projection operator to closed set $\left[\Omega ;{\mathbb{R}}^{m};{\mathbb{R}}_{+}^{ab}\right]$. Based on Lemma 1, the proposed neural network (15) is stable and will globally converge to the optimal solution of (11).

Notable that the equality constraint 11(b) describes the impedance controller, and the convergence can be found in Na et al. (2015). Similarly, the establishment of inequality constraint enables obstacle avoidance during the whole process. The proof is completed.

*Remark 4*. It is remarkable that the original impedance controller described in 11(b) bears similar with traditional methods in Yang et al. (2018a) the main contribution of the proposed controller is that the controller can not only realize the force control, but also realize the obstacle avoidance, besides, the control strategy is capable of handling inequality constraints, including joint angles, and velocities.

## 4. Numerical Results

In this part, we carry out a series of numerical simulations on a planar 4-DOF robot, aiming at verifying the validity of the proposed control scheme. Firstly, a pure force control experiment is done to show the effectiveness of the force controller, and then the control scheme is further verified by examining the system response after introduction of obstacles. Then we check the control performance in more general situations, including model uncertainties and multiple obstacles.

### 4.1. Simulation Settings

First of all, the planar robot used in the simulation is show in Figure 2. The D-H parameters are also listed in Figure 2B. It is remarkable that in force control tasks, the end-effector is required to keep in touch with workpieces, which makes it necessary to distinguish between the necessary contact and the unnecessary collisions. In this paper, the proposed controller is capable of handling this problem by selecting the key points properly. Therefore, the end-effector is not considered as a key point, to make it possible to contact with the obstacles (or external environment). In order to avoid obstacles, the set of key points of the robot is defined as *A*_{1}, ⋯ , *A*_{7}, in which *A*_{1}, *A*_{3}, *A*_{5}, and *A*_{7} locate at the center of the links, and *A*_{2}, *A*_{4}, and *A*_{6} are defined to be at *J*2, *J*3, and *J*4, as shown in Figure 2A. The lower and upper bounds of joint angles and joint velocities are defined as ${\theta}_{i}^{-}=-3$rad, ${\theta}_{i}^{+}=3$rad, ${\stackrel{.}{\theta}}_{i}^{-}=-1$rad/s, ${\stackrel{.}{\theta}}_{i}^{+}=1$rad/s for *i* = 1…4, respectively. The safety margin is selected as 0.01 m. The coefficients describing the contact force are selected as *K*_{d} = 50, *K*_{p} = 5000. For simplicity, let ${b}_{0}={K}_{d}^{-1}F-{K}_{p}{K}_{d}^{-1}\Delta x+{\u1e8b}_{\text{d}}$.

**Figure 2**. The robot to be simulated in this paper. **(A)** Is the physical structure and the location of key-points. **(B)** Is D-H parameters.

### 4.2. Force Control Without Obstacles

First of all, an ideal case where there is no obstacles in the workspace is considered, and the parameters *K*_{d} and *K*_{p} are assumed to be known. The robot is wished to offer a constant contact force on a given plane. The contact force is set to be 20N, while the direction of contact force is aligned with the y-axis of the tool coordination system. In this example, the y-axis of is [1, −1]^{T} in the base coordination. The pre-defined path on the contact plane is *x*_{d} = [0.4+0.1*cos*(0.5*t*), 0.5+0.1*cos*(0.5*t*)]. The initial state of the robot system is set as ${\theta}_{0}={\left[1.57,-0.628,-0.524,-0.524\right]}^{\text{T}}$rad, ${\stackrel{.}{\theta}}_{0}={\left[0,0,0,0\right]}^{\text{T}}$rad/s. The control gains of the proposed RNN controller are α = 8,ϵ = 0.02, respectively. Numerical results are shown in Figure 3. The tracking error along the contact plane is given in Figure 3B, the transient is about 1s. At the beginning stage, since the end-effector is not in contact with the surface, the contact force stays zero before 0.5s. As the end-effector approaches the surface, the contact force converges to 20N, showing the convergence of both positional and force errors. The Euclidean norm of joint velocities (which is also output of the established RNN) is shown in Figure 3D, $||\stackrel{.}{\theta}||$ changes periodically, with the same cycle as the expected trajectory. The time history of the end-effector's motion trajectory and the corresponding joint configurations are shown in Figure 3A, in which the red arrow indicates the direction of the contact force, and the blue arrow shows the direction of the end-effector's free-motion. All in all, the proposed controller can achieve the position-force control precisely.

**Figure 3**. Numerical results of compliance control without obstacles. **(A)** Is the robot's tracking path and the corresponding joint configurations. **(B)** Is the profile of position error along the free-motion direction. **(C)** Is the profile of contact force. **(D)** Is the profile of $||\stackrel{.}{\theta}{||}_{2}^{2}$.

### 4.3. Force Control With Single Obstacles

In this section, a stick obstacle is introduced into the workspace, which is defined as *x* = −0.05 m. The initial states and expected values of *x*_{d}, *F*_{d} are the same as section 4.2.

*Remark 5*. In Equation (10), we have shown the basic idea of calculating the distance between the robot and obstacles, i.e., by abstracting key points form the robot and obstacles, the distances can be the robot and obstacle can be described approximately at a set of point-to-point distances. In this example, the distance can be obtained in a simpler way. However, the obstacle avoidance strategy is essentially consistent with (Equation 10).

Simulation results are given in Figures 4, 5. The output of RNN is shown in Figure 4E, when simulation begins, $\stackrel{.}{\theta}$ reaches its maximum value, driving the end-effector to move toward the desired path. And then the robot slows down quickly (after *t*≈0.5s), the robot move smoothly, as a result, the position error successfully converges to 0, and simultaneously, the contact force converges to 20N. It is notable that at *t* = 1.2 s, the key point *A*_{2} of the robot gets close to the obstacle, as shown in Figure 4F. Based on the obstacle avoidance strategy (Equation 15c), the state variable λ_{2}(2) becomes positive, and then the output of the RNN varies with λ_{2} (Figure 5B). Correspondingly, an error (about 1 × 10^{−3} m) occurs in the positional tracking, and so as the contact force (force error is about 2N). However, the RNN converges to the new equilibrium point(since the equilibrium point would change when the inequality constraint works), and both *e*_{x} and *e*_{f} converges to 0. By comparing Figures 3A, 4A, after introducing the obstacle, the robot is capable of adjusting its joint configuration to avoid the obstacle. The distances between the key points *A*_{1}−*A*_{7} to the obstacle are shown in Figure 4D, a minimum value of about 0.01 m is ensured during the whole process. Using impedance model, the force control problem is transferred into a kinematic control one by modifying the reference speed (Equation 4). Consequently, the resulting trajectory *x*_{r} together with *x*_{d} are as shown in Figures 5D,E. As an important index in the proposed control scheme, the norm of joint speed $||\stackrel{.}{\theta}{||}_{2}^{2}$ is wished as small as possible. Therefore, we introduce a comparative simulation, in which the solution is obtained based on pseudo-inverse of Jacobian matrix and physical limitations are not considered. Comparative curves of the objective functions are as shown in Figure 5F. The RNN based controller can optimize the objective function, it is remarkable that a difference appears at about *t* = 1.2−5 s, which is mainly caused by obstacle avoidance (which is not considered in JMPI based method). Since the output of RNN $\stackrel{.}{\theta}$ is used to approximate the reference speed *b*_{0}, the approximate error $||J\stackrel{.}{\theta}-{b}_{0}{||}_{2}^{2}$ is shown in 55(C), demonstrating the effectiveness of the established RNN.

**Figure 4**. Control performance of the proposed controller while avoiding a wall obstacle. **(A)** Is the robot's tracking path and the corresponding joint configurations. **(B)** Is the profile of position error along the free-motion direction. **(C)** Is the profile of contact force. **(D)** Is the profile of joint angles. **(E)** Is the profile of joint velocities. **(F)** Is the profile of the closest distance to the obstacle of each key points *A*_{i}, *i* = 1, ⋯ , 7.

**Figure 5**. Simulation results of the established RNN while avoiding a wall obstacle. **(A)** Is the profile of λ_{1}. **(B)** Is the profile of λ_{2}. **(C)** Is the profile of $||J\stackrel{.}{\theta}-{b}_{0}{||}_{2}^{2}$. **(D)** Is the profiles of the desired and reference trajectory along x-axis. **(E)** Is the profiles of the desired and reference trajectory along y-axis. **(F)** Is the profiles of the objective function of the proposed controller and JPMI based method.

### 4.4. Force Control With Uncertain Parameters

In this example, we check the control performance of the proposed control scheme in presence of model uncertainties. Similar with previous simulations, the initial states of the robot are also ${\theta}_{0}={\left[1.57,-0.628,-0.524,-0.524\right]}^{\text{T}}$rad, ${\stackrel{.}{\theta}}_{0}={\left[0,0,0,0\right]}^{\text{T}}$rad/s. In real implementations, the interaction model is usually unknown, and the nominal values of *K*_{d} and *K*_{p} are not accurate. Without loss of generality, we select the nominal values of *K*_{d} and *K*_{p} as ${\widehat{K}}_{d}=80$, ${\widehat{K}}_{p}=4000$, respectively.In order to handle model uncertainties in the interaction coefficients, an extra node is introduced into (15). Then the modified RNN can be formulated as:

in which *W* = [*K*_{p}; *K*_{d}], η = [*x* − *x*_{d}; ẋ − ẋ_{d}], and the positive coefficient *K*_{in} scaling the updating rate is defined as *K*_{in} = *diag*(500, 20). Simulation results are shown in Figures 6, 7. Although the exact values of *K*_{d} and *K*_{p} are unknown, the closed-loop system is still stable, which can be shown from the convergence of tracking error *e*_{x} and contact force *F* in Figures 6A,B. The change curves of joint angles and joint velocities with respect to time are shown in Figures 6C,D, in which the bounded-ness of joint angles and velocities are guaranteed. The observed interaction coefficients ${\widehat{K}}_{d}$ and ${\widehat{K}}_{p}$ are shown in Figure 6E, indicating that both ${\widehat{K}}_{d}$ and ${\widehat{K}}_{p}$ converge to their real values. Figure 7A shows the distances between the key points and the obstacle, it is obvious that all key points keep at a safe distance from the obstacle (the closest key point is *A*_{2}). Euclidean norm of ${b}_{0}-J\stackrel{.}{\theta}$ is illustrated in Figure 7C, despite fluctuation occurs at about *t* = 1.5 s, the proposed controller could handle model uncertainties. The impedance model based reference trajectory and the original desired trajectory are shown in Figures 7D,E. Although *x*_{r} and *x*_{d} are different, the tracking error *e*_{x} along the direction of free motion and force error *e*_{F} converges to zero, as shown in Figures 6A,B. The objective function $||\stackrel{.}{\theta}{||}_{2}^{2}$ to be optimized is given in Figure 7F. the convergence of the established RNN is shown in Figure 7C, despite the uncertain parameters, using the adaptive updating law, the established RNN is capable of learning the optimal solution. The spikes is mainly because of the change of λ_{2} when obstacle avoidance scheme is activated.

**Figure 6**. Control performance of the proposed controller while avoiding a wall obstacle with uncertain *K*_{p} and *K*_{d}. **(A)** Is the robot's tracking path and the corresponding joint configurations. **(B)** Is the profile of position error along the free-motion direction. **(C)** Is the profile of contact force. **(D)** Is the profile of joint angles. **(E)** Is the profile of joint velocities. **(F)** Is the profile of the closest distance to the obstacle of each key points *A*_{i}, *i* = 1, ⋯ , 7.

**Figure 7**. Simulation results of the established RNN while avoiding a wall obstacle with uncertain *K*_{p} and *K*_{d}. **(A)** Is the profile of λ_{1}. **(B)** Is the profile of λ_{2}. **(C)** Is the profile of $||J\stackrel{.}{\theta}-{b}_{0}{||}_{2}^{2}$. **(D)** Is the profiles of the desired and reference trajectory along x-axis. **(E)** Is the profiles of the desired and reference trajectory along y-axis. **(F)** Is the profiles of the objective function of the proposed controller and JPMI based method.

### 4.5. Manipulation in Narrow Space

In this part, we discuss a more general case of motion-force control task, in which the workspace is defined in a limited narrow space. The robot is limited by two parallel lines, namely, *y*_{1} = 0.15 and *y*_{2} = −0.15 m. Considering the safety distance, all key points except *A*_{8} must satisfy the workspace description −0.14 ≤ *y* ≤ 0.14 m. The initial joint angles are set to be ${\theta}_{0}={\left[0.393,-1.05,1.57,-0.785\right]}^{\text{T}}$rad, and ${\stackrel{.}{\theta}}_{0}={\left[0,0,0,0\right]}^{\text{T}}$rad/s. The desired path is selected as ${x}_{\text{d}}={\left[0.8+0.1cos(0.5t),-0.15\right]}^{\text{T}}$ m, and the expected contact force is *F*_{d} = 20N, with the direction vector being [0, −1]^{T}. Simulation results are given in Figures 8, 9. When simulation begins, the initial position error is about 0.1 m, and the converges to zero, with the transient being about 0.5s. Simultaneously, the contact force also converges to 20N. In Figure 9A, minimum distances between the key points to *y*_{1} and *y*_{2} are represented by blue and red curves, respectively. The tracking trajectory and the corresponding joint configurations are shown in Figure 8A. During *t* = 1 − 1.5 s and *t* = 6 − 13 s, point *A*_{2} gets close to *y*_{1}, during *t* = 4 − 7 s, *A*_{4} is close to *y*_{2}. Remarkable that there exist fluctuations in positional and force errors at *t* = 1 s and *t* = 4 s (i.e., when *A*_{2} and *A*_{4} get close to the bounds), respectively. Similar to the previous simulations, the reference trajectories are given in Figures 7C,D, and the objective functions are shown in Figure 7E. Using the proposed RNN controller, the robot can realize both position and force control in limited narrow space.

**Figure 8**. Control performance of the proposed controller in a narrow workspace. **(A)** Is the robot's tracking path and the corresponding joint configurations. **(B)** Is the profile of position error along the free-motion direction. **(C)** Is the profile of contact force. **(D)** Is the profile of joint angles. **(E)** Is the profile of joint velocities. **(F)** Is the profile of the closest distance to the obstacle of each key points *A*_{i}, *i* = 1, ⋯ , 7.

**Figure 9**. Simulation results of the established RNN in a narrow workspace. **(A)** Is the profile of λ_{1}. **(B)** Is the profile of λ_{2}. **(C)** Is the profiles of the desired and reference trajectory along x-axis. **(D)** Is the profiles of the desired and reference trajectory along y-axis. **(E)** Is the profiles of the objective function of the proposed controller and JPMI based method.

### 4.6. Comparisons

In this part, comparisons among the proposed control scheme and existing methods are given to show the superiority of the RNN based strategy. The comparisons are shown in Table 1. In Guo and Zhang (2012), a RNN based controller is designed for redundant manipulators, both obstacle avoidance and physical constraints are considered. However, the controller only focus on kinematic control problem. In Nanayakkara et al. (2001) and Csiszar et al. (2011), force control together with obstacle avoidance are taken into account, but the physical constraints are ignored. Xu et al. (2019a) develop an adaptive admittance control strategy, which is capable of dealing with force control under model uncertainties, physical constraints and real-time optimization. It is remarkable that the proposed strategy focus on real-time obstacle avoidance in force control tasks, and the controller is capable of ensuring the boundedness of joint angles and velocities. At the same time, simulations have shown the potential of optimization ability of norm of joint speed.

## 5. Conclusions

In this paper, a novel collision-free compliance controller is constructed based on the idea of QP programming and neural networks. Different with existing methods, in this paper, the control problem is described from an optimization perspective, and the compliance control and collision avoidance are formulated as equality or inequality constraints. The physical constraints such as limitations of joint angles and velocities are also taken into consideration. Before ending this paper, it is worth pointing out that it is the first RNN based compliance control method, which considers collision avoidance problem in realtime, and also shows great potential in handling physical limitations. In this paper, simple numerical simulations in MATLAB are carried out to verify the efficiency of the proposed controller. In the future, we will check the control framework with different impedance models in physically realistic simulation environments, and then consider machine vision technology and system delay problem on physical experimental platforms.

## Data Availability

All datasets analyzed for this study are included in the manuscript and the supplementary files.

## Author Contributions

XZ designed the control strategy and writes the paper. ZX did the numerical experiments. SL revised the paper and proposed amendments.

## Funding

This work was supported by the GDAS' Foundation of National and Provincial Science and Technology Talent (Grant No. 2018GDASCX-0603), Guangdong Province Applied Science and Technology Research Project (Grant No. 2015B090922010), Guangdong Province Public Welfare Research and Capacity Building Project (Grant No. 2017A010102017), Guangdong Province Science and Technology Major Projects (Grant No. 2016B090910003), Guangdong Province Science and Technology Major Projects (Grant No. 2017B010116005), Guangdong Province Key Areas R&D Program (Grant No. 2019B090919002).

## Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## Acknowledgments

The authors would like to thank the reviewers for their valuable comments and suggestions.

## References

Badawy, A. (2016). Dual-well potential field function for articulated manipulator trajectory planning. *Alexand. Eng. J.* 55, 1235–1241. doi: 10.1016/j.aej.2016.03.042

Cai, M., and Xiang, Z. (2018). Adaptive finite-time control of a class of non-triangular nonlinear systems with input saturation. *Neural Comput. Appl.* 29, 565–576. doi: 10.1007/s00521-016-2540-0

Csiszar, A., Drust, M., Dietz, T., Verl, A., and Brisan, C. (2011). Dynamic and interactive path planning and collision avoidance for an industrial robot using artificial potential field based method. *Mechatronics* 1, 413–421. doi: 10.1007/978-3-642-23244-2-50

Gao, X. (2003). Exponential stability of globally projected dynamic systems. *IEEE Trans. Neural Netw.* 14, 426–431. doi: 10.1109/TNN.2003.809409

Gonzalez, A., Alves, M., Viana, G., Carvalho, L., and Basilio, J. (2018). Supervisory control-based navigation architecture: a new framework for autonomous robots in industry 4.0 environments. *IEEE Trans. Indust. Inform.* 14, 1732–1743. doi: 10.1109/TII.2017.2788079

Guo, D., and Zhang, Y. (2012). A new inequality-based obstacle-avoidance mvn scheme and its application to redundant robot manipulators. *IEEE Trans. Syst. Man Cybernet. Part C* 42, 1326–1340. doi: 10.1109/TSMCC.2012.2183868

He, W., Ge, S. S., Li, Y. N., and Chew, E. (2015). Neural network control of a rehabilitation robot by state and output feedback. *J. Intell. Robot. Syst.* 80, 15–31. doi: 10.1007/s10846-014-0150-6

Hogan, N. (1985). Impedance control - an approach to manipulation. i - theory. ii - implementation. iii - applications. *Asme Trans. J. Dynam. Syst. Measure. Control B* 107, 304–313. doi: 10.1115/1.3140713

Huang, B., Li, Z., Wu, X., Ajoudani, A., Bicchi, A., and Liu, J. (2019). Coordination control of a dual-arm exoskeleton robot using human impedance transfer skills. *IEEE Trans. Syst. Man Cybernet. Syst.* 49, 954–963. doi: 10.1109/TSMC.2017.2706694

Jin, L., Li, S., and La, H. (2017). Manipulability optimization of redundant manipulators using dynamic neural networks. *IEEE Trans. Indust. Electr.* 64, 4710–4720. doi: 10.1109/TIE.2017.2674624

Khatib, O. (1986). Real-time obstacle avoidance system for manipulators and mobile robots. *Int. J. Robot. Res.* 5, 90–98.

Khatib, O. (1987). A unified approach for motion and force control of robot manipulators: the operational space formulation. *IEEE J. Robot. Automat.* 3, 43–53.

Li, S., Zhang, Y., and Jin, L. (2017). Kinematic control of redundant manipulators using neural networks. *IEEE Trans. Neural Netw. Learn. Syst.* 28, 2243–2254. doi: 10.1109/TNNLS.2016.2574363

Li, Y., Li, S., and Hannaford, B. (2019). A model based recurrent neural network with randomness for efficient control with applications. *IEEE Trans. Indust. Inform.* 15, 2054–2063. doi: 10.1109/TII.2018.2869588

Na, J., Nasiruddin, M., Guido, H., Ren, X., and Phil, B. (2015). Robust adaptive finite-time parameter estimation and control for robotic systems. *Int. J. Robust Nonlinear Control* 25, 345–3071. doi: 10.1002/rnc.3247

Nanayakkara, D., Kiguchi, K., Murakami, T., Watanabe, K., and Izumi, K. (2001). “Skillful adaptation of a 7-dof manipulator to avoid moving obstacles in a teleoperated force control task,” in *ISIE 2001. 2001 IEEE International Symposium on Industrial Electronics Proceedings (Cat. No. 01TH8570)* (Busan: IEEE), 1982–1988.

Pan, Y., Li, X., Wang, H., and Yu, H. (2018a). Continuous sliding mode control of compliant robot arms: a singularly perturbed approach. *Mechatronics* 52, 127–134. doi: 10.1016/j.mechatronics.2018.04.005

Pan, Y., Li, X., and Yu, H. (2019). Efficient pid tracking control of robotic manipulators driven by compliant actuators. *IEEE Trans. Control Syst. Techn.* 27, 915–955. doi: 10.1109/TCST.2017.2783339

Pan, Y., Wang, H., Li, X., and Yu, H. (2018b). Adaptive command-filtered backstepping control of robot arms with compliant actuators. *IEEE Trans. Control Syst. Techn.* 26, 1149–1156. doi: 10.1109/TCST.2017.2695600

Raibert, M. H., and Craig, J. J. (1981). Hybrid position/force control of manipulators. *Asme J. Dynam. Syst. Measure Control* 102, 126–133.

Sciavicco, L., and Siciliano, B. (1988). A solution algorithm to the inverse kinematic problem for redundant manipulators. *IEEE J. Robot. Automat.* 1, 403–410.

Senoo, T., Koike, M., Murakami, K., and Ishikawa, M. (2017). Impedance control design based on plastic deformation for a robotic arm. *IEEE Robot. Automat. Lett.* 2, 209–2061. doi: 10.1109/LRA.2016.2587806

Tsai, C., Lee, J., and Chuang, J. (2001). Path planning of 3-d objects using a new workspace model. *IEEE Trans. Systems Man Cybernet. Part C* 31, 405–410. doi: 10.1109/5326.971669

Tsuji, T., Tanaka, Y., Morasso, P., Sanguineti, V., and Kaneko, M. (2002). Bio-mimetic trajectory generation of robots via artificial potential field with time base generator. *IEEE Trans. Syst. Man. Cybernet. Part C* 32, 426–439. doi: 10.1109/TSMCC.2002.807273

Wang, H., Liu, X., and Liu, K. (2015). Robust adaptive neural tracking control for a class of stochastic nonlinear interconnected systems. *IEEE Trans. Neural Netw. Learn. Syst.* 27, 510–523. doi: 10.1109/TNNLS.2015.2412035

Wang, S., Zhang, J., and Zhang, J. (2018). Artificial potential field algorithm for path control of unmanned ground vehicles formation in highway. *Electr. Lett.* 54, 1166–1168. doi: 10.1049/el.2018.5018

Xu, Z., Li, S., Zhou, X., and Cheng, T. (2019a). Dynamic neural networks based adaptive admittance control for redundant manipulators with model uncertainties. *Neurocomputing* 357, 271–281. doi: 10.1016/j.neucom.2019.04.069

Xu, Z., Li, S., Zhou, X., Yan, W., Cheng, T., and Huang, D. (2019b). Dynamic neural networks based kinematic control for redundant manipulators with model uncertainties. *Neurocomputing* 329, 255–266. doi: 10.1016/j.neucom.2018.11.001

Yang, C., Jiang, Y., He, W., Na, J., Li, Z., and Xu, B. (2018a). Adaptive parameter estimation and control design for robot manipulators with finite-time convergence. *IEEE Trans. Indust. Electr.* 65, 8112–8123. doi: 10.1109/TIE.2018.2803773

Yang, C., Peng, G., Li, Y., Cui, R., Cheng, L., and Li, Z. (2018b). Neural networks enhanced adaptive admittance control of optimized robot-environment interaction. *IEEE Trans. Cybernet.* 49, 2568–2579. doi: 10.1109/TCYB.2018.2828654

Zhang, T., and Xia, J. (2019). Interconnection and damping assignment passivity-based impedance control of a compliant assistive robot for physical human-robot interactions. *IEEE Robot. Automat. Lett.* 4, 538–545. doi: 10.1109/LRA.2019.2891434

Zhang, Y. (2015). Singularity-conquering tracking control of a class of chaotic systems using zhang-gradient dynamics. *IET Control Theory Appl.* 9, 871–881. doi: 10.1049/iet-cta.2014.0931

Zhang, Y., Chen, S., Li, S., and Zhang, Z. (2018). Adaptive projection neural network for kinematic control of redundant manipulators with unknown physical parameters. *IEEE Trans. Indust. Electr.* 65, 4909–4920. doi: 10.1109/TIE.2017.2774720

Zhang, Y., Ge, S., and Lee, T. (2004). A unified quadratic-programming-based dynamical system approach to joint torque optimization of physically constrained redundant manipulators. *IEEE Trans. Systems Man Cybernet Part B* 34, 2126–2132. doi: 10.1109/TSMCB.2004.830347

Zhang, Y., and Wang, J. (2004). Obstacle avoidance for kinematically redundant manipulators using a dual neural network. *IEEE Trans. Syst. Man Cybernet. Part B* 34, 752–759. doi: 10.1109/TSMCB.2003.811519

Zhao, Z., He, X., Ren, Z., and Wen, G. (2018a). Boundary adaptive robust control of a flexible riser system with input nonlinearities. *IEEE Trans. Syst. Man Cybernet. Syst.* 1–10. doi: 10.1109/TSMC.2018.2882734

Keywords: recurrent neural network, compliance control, redundant manipulator, obstacle avoidance, zeroing neural network

Citation: Zhou X, Xu Z and Li S (2019) Collision-Free Compliance Control for Redundant Manipulators: An Optimization Case. *Front. Neurorobot.* 13:50. doi: 10.3389/fnbot.2019.00050

Received: 30 April 2019; Accepted: 24 June 2019;

Published: 11 July 2019.

Edited by:

Yongping Pan, National University of Singapore, SingaporeReviewed by:

Predrag S. Stanimirovic, University of Niš, SerbiaMuhammad Fahad, Stevens Institute of Technology, United States

Jing Na, Kunming University of Science and Technology, China

Copyright © 2019 Zhou, Xu and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Zhihao Xu, zh.xu@giim.ac.cn