# International Journal of Robotic Engineering

(ISSN: 2631-5106)

Volume 5, Issue 1

Research Article

#### DOI: 10.35840/2631-5106/4124

# Fractional Order Control for the Grasping by Coiling

Mircea Ivanescu^{1*}, Mihaela Florescu^{2} and Cristian Vladu^{3}

### Table of Content

### Figures

**Figure 10:** Circle criterion for continuum arm with.....

Circle criterion for continuum arm with high gravitational component.

### References

- Deghami M, Moosavin S (2014) Dynamics modeling of a continuum robotic arm with a contact point in planar grasp. Journal of Robotics 308283.
- Robinson G, Davies JBC (1999) Continuum robots-A state of the art. IEEE Int Conf on Robotics and Automation, Detroit, USA, 2849-2854.
- Gravagne Ian A, Walker ID (2000) On the kinematics of remotely-actuated continuum robots. Proc 2000 IEEE Int Conf on Rob and Aut, San Francisco, 2544-2550.
- Gravagne IA, Walker ID (2000) Kinematic transformations for remotely-actuated planar continuum robots. Proc 2000 IEEE Int Conf on Rob and Aut, San Francisco, 19-26.
- Gravagne IA, Walker ID (2002) Uniform regulation of a multi-section continuum manipulators. Proc 2002 IEEE Int Conf on Rob and Aut, Washington, USA, 1519-1525.
- Gravagne I, Walker ID (2001) Manipulability and force ellipsoids for continuum robot manipulators. Proc 2001 IEEE/RSJ International Conference on Intelligent Robors and Systems, Maui, Hawaii, USA, 304-310.
- Chirikjian GS, Burdick JW (1990) An obstacle avoidance algorithm for hyper-redundant manipulators. Proc IEEE Int Conf on Rob and Aut, Cincinnati, Ohio, USA, 625-631.
- Mochiyama H, Kobayashi H (1999) The shape Jacobian of a manipulator with hyper degrees of freedom. Proc 1999 IEEE Int Conf on Robotics and Automation, Detroit, USA, 2837-2842.
- Jinglin L, Jing X (2011) Determining grasping configurations for a spatial continuum manipulator. 2011 IEEE/RSJ Int Conf on Intelligent Robots and Systems, San Francisco, USA, 4207-4213.
- Walker IM, Hannan M (1999) A novel elephant's trunk robot. AIM 99: 410-415.
- Jones B, Walker ID (2006) Practical kinematics for real-time implementation of continuum robots. IEEE Trans Robotics 22: 1087-1099.
- Kapadia A, Walker I, Dawson DM (2009) A model-based sliding mode controller for extensible continuum robots. Recent Advances in Signal Processing, Robotics and Automation, 103-120.
- Rucker DC, Webster RJIII, Chirikjian GS, Cowan NJ (2010) Equilibrium conformations of concentric-tube continuum robots. Int Journal of Robotic Research 29: 1263-1280.
- Popescu N, Popescu D, Ivanescu MA (2013) Spatial weight error control for a class of hyper- redundant robots. IEEE Trans on Robotics 29: 1043-1050.
- Ivanescu M, Popescu N, Popescu D (2015) The shape control of a tentacle arm. Robotica Cambridge Journal 33: 684-703.
- Ivanescu M, Popescu D, Popescu N (2015) A decoupled sliding mode control for a continuum arm. Adv Robotics, Special Issue: Continuum Robots and Manipulation 29: 831-845.
- Hammond IIIF, Kramer R, Wan Q, Howe D, Wood R (2014) Soft tactile sensor arrays for force feedback. IEEE Sensors Journal 14: 1443-1452.
- Diethelm K (2004) The analysis of fractional differential equations. Springer- Verlag, London.
- Petras I (2011) Fractional-order nonlinear systems, modeling, analysis and simulation. Higher Education Press, Springer-Verlag, London.
- Monje C, Chen YQ, Vinagre B, Hue D, Feliu V (2010) Fractional-order systems and controls. Springer Verlag London Limited.
- Norelys AC, Manuel DM, Javier AG (2014) Lyapunov functions for fractional order systems. Commun Nonlinear Sci Numer Simulat 19: 2951-2957.
- Khalil H (2003) Nonlinear Systems, Prentice Hall, New Jersey, USA.
- Li Y, Chen YQ, Podlubny I (2010) Stability of fractional order nonlinear systems: Lyapunov direct method and generalized Mittag-Leffler stability. Comput Math Appl 59: 1810-1821.
- Argarwal R, Hristova S, Regan DO (2015) Lyapunov functions and strict stability of Caputo fractional differential equations. Advances in Difference Equations 346: 2-20.
- Zhao Y, Wang Y, Liu Z (2015) Lyapunov function method for linear fractional order systems. Proc of the 34th Chinese Control Conference, 1457-1463.
- Dadras S, Soodeh D, Malek H, Chen Y (2017) A note on the lyapunov stability of fractional order nonlinear systems. In: Proc of the ASME Cleveland, Orlando, 123-129.
- Slotine JE, Weiping Li (1991) Applied nonlinear control. Prentice Hall, USA.
- Al-Saggaf UM, Mehedi IM, Rachid M, Maamar B (2017) Rotary flexible joint control by fractional order controllers. International Journal of Control, Automation and Systems 15: 2561-2569.
- Rhong L, X Peng, Zhang B (2018) A reduced-order fault detection filter design for polytopic uncertain continuous-time markovian jump systems with time-varying delays. International Journal of Control, Automation and Systems 16: 2021-2032.
- Khimani D, Patil M (2018) High performance super-twisting control for state delay systems. International Journal of Control, Automation and Systems 16: 2063-2073.
- Wu B, Wang CL, Hu YJ (2018) Stability analysis for time-delay systems with nonlinear disturbances via new generalized integral inequalities. International Journal of Control, Automation and Systems 16: 2772-2780.
- Heymann N, Podlubni I (2014) Physical interpretation of initial conditions for fractional differential with Riemann-Liouville fractional derivatives. Rheologica Acta 23: 45-63.

**Author Details**

Mircea Ivanescu^{1*}, Mihaela Florescu^{2} and Cristian Vladu^{3}

^{1}Department of Mechatronics, University of Craiova, Romania

^{2}Department of Mechanical Engineering, University of Craiova, Romania

^{3}Department of Electrical Engineering, University of Craiova, Romania

**Corresponding author**

Mircea Ivanescu, Department of Mechatronics, University of Craiova, Romania.

Accepted: May 14, 2020 | Published Online: May 16, 2020

Citation: Ivanescu M, Florescu M, Vladu C (2020) Fractional Order Control for the Grasping by Coiling. Int J Robot Eng 5:024.

Copyright: © 2020 Ivanescu M, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

## Abstract

The paper studies the grasping control by coiling for a class of grippers with continuum arms. First, the fractional order model of continuum arm is proposed. The viscoelastic components are assimilated by the fractional Kevin-Voigt model and the fractional-order dynamics are inferred. A hybrid control technique with two control loops for position and force control, respectively, is proposed. The Lyapunov method for position control is applied. Sector-type constraints for input torque are implemented and frequential conditions, Popov Circle Criterion, that ensure asymptotic stability, are inferred. A conventional PD controller is proposed for the force control. Numerical simulations illustrate the performances of the control system.

## Keywords

Continuum arm, Fractional-order model, Control

## Introduction

In this paper, a class of arms with continuum elements that performs the grasping function by coiling is analysed. This function is often met in the animal world as octopus tentacle, the elephant's trunk or constrictor snakes. In robotic applications, this function is used to manipulate objects with geometric shapes that contain borders with concave areas (Figure 1), to objects with certain fragility and resistance to manipulate, etc. In this area of application, the architectures of continuum arm type play an important factor in determining the control laws. The specialized literature contains a large number of papers addressing the study of these models known as systems with distributed parameters. A summary of these results can be found in [1]. Among the pioneering works we can mention [2-4] which introduce the basic concepts of these models, the state equations. The dynamic models and the control laws for constrained motions are inferred in [5,6]. In [7] are discussed kinematic algorithms to avoid the obstacles of hyper-redundant manipulators. An extension of Jacobian based methods for infinite dimensional systems is presented in [8]. In [9] are studied the grasping configurations of continuum arms. Similarities of these systems with certain classes of biological models are analysed in [10]. Kinematic control for real time implementation is discussed in [11]. Sliding mode control techniques for hyper-redundant robots are studied in [12]. Architectures based on concentric tubes for the continuum arm are presented in [13]. Spatial weight error control for a class of hyper-redundant robots is treated in [14]. Shape control of the tentacle robots is presented in [15,16].

This paper studies the grasping control by coiling for a class of grippers with continuum arms. First, the fractional order model of continuum arm is proposed. The viscoelastic components are assimilated by the fractional Kevin-Voigt model and the fractional-order dynamics are inferred. A hybrid control technique with two control loops for position and force control, respectively, is proposed. The Lyapunov method for position control is applied. Sector-type constraints for input torque are implemented and frequential conditions that ensure asymptotic stability, similar to the Popov Circle Criterion for classical systems, are inferred. A conventional PD controller is proposed for the force control. Numerical simulations illustrate the performances of the control system.

The paper is structured as follows: Section 2 presents technological principles, section 3 treats the dynamic model, section 4 describes the control system, section 5 verifies the control techniques by numerical simulations, section 6 presents discussions and section 7 resumes the conclusions.

## Technological Principles

The main component of the grasping system is the continuum arm that manipulates the object by exerting a sufficiently high pressure on the contact surface so that it can execute its displacement on a predetermined path. We denote this system as "Continuum Arm Grasping System (CAGS)". CAGS can achieve any position and orientation in 3-D space to surround the manipulated load on the surface of different geometric shapes. CAGS constitutes the terminal component of a conventional robot intended for current part handling functions. The present work is focused only on the control of the grasping operations.

The main component of the CAGS mechanical structure is a long flexible backbone. It is surrounded by an elastic envelope and a viscoelastic solution fill the whole system ensuring a large flexibility (Figure 2). The actuation system consists of antagonistic Bowden cables with DC motor driving. The mechanical architecture is composed of several segments with independent Bowden cable driving system that allows a decoupled control strategy for each arm segment. The measuring system consists of two components: A measuring system of the geometric parameters of the orientation of the arm and one for measuring the compression forces during the grasping operation. The first system is achieved by measuring the lengths of the Bowden cables [16]. The second system is consisted by an array of Force Sensitive Resistors (FSR) (Figure 2). They are constructed of several thin flexible layers that varies in resistance as pressure is applied and released [17].

The backbone of the arm is defined by a 3-dimensional curve *C* and a length parameter *s*, $s\in \text{}\Omega $, $\Omega \text{=}\left[0.\text{}l\right]$, where *l* is the length of the arm. The arm orientation is given by two angles ${u}_{1}$, ${u}_{2}\text{}\in \text{}\Gamma \left(\Omega \right)\subset {L}_{2}\left(\Omega \right)$, defined by a right-handed orthonormal basis vector $\{{i}_{x},{i}_{y},{i}_{z}\}$, where $\{{i}_{x},{i}_{y},{i}_{z}\}$ is parallel to a base coordinate frame, Figure 3.

## Dynamic Model

The dynamic behaviour of the CAGS is determined by its structure with distributed mass and by the viscoelastic component characterized by the Kevin-Voigt fractional order model defined by the stress-strain $\left(f,\text{}\varrho \right)$ relation

$$f\left(t\right){\text{=c}}_{e}\varrho \left(t\right){\text{+c}}_{v}\left({}_{t}{}^{C}D{}_{0}^{\beta}\varrho \right)\left(t\right)\text{(1)}$$

Where $\left({}_{t}{}^{C}D{}_{0}^{\beta}\varrho \right)\left(t\right)$ is the Caputo fractional derivative of order β. The dynamic model is derived by using Lagrange equations [5,14] for this class of hyper-redundant arms. We consider a light weight arm with following parameters: ${I}_{\rho}$ is inertial rotational density, E is equivalent elastic modulus and *I* is the moment of inertia. It is assumed that during the grasping operation, the arm is subjected to a disturbance, the moment *m*, distributed along the contact surface. In terms of the orientation variables ${u}_{1},{u}_{2}$ this model can be rewritten as

$${I}_{\rho}\ddot{u}=EI{u}_{ss}-{C}_{v}{D}^{\beta}u-{C}_{e}u+h\text{(2)}$$

Where

$u\text{=}u\left(t,s\right)$, ${u}_{ss}\text{=}\frac{\partial {}^{2}u\left(t,s\right)}{\partial s{}^{2}}$, $u\text{=}{\left[{u}_{1}{u}_{2}\right]}^{T}$, ${u}_{1},{u}_{2}\in {L}_{2}\left(\Omega \right)$, $D{}^{\beta}u{\text{=}}_{t}^{C}D{}_{0}{}^{\beta}u\left(t,s\right)$, ${I}_{\rho}\text{=}diag\left({I}_{\rho 1},{I}_{\rho 2}\right)$, $EI\text{=}diag\left(E{I}_{1},E{I}_{2}\right)$, ${C}_{v}\text{=}diag\left({c}_{v1},{c}_{v2}\right)$, ${C}_{e}\text{=}diag\left({c}_{e1},{c}_{e2}\right)$, ${C}_{m}\text{=}diag\left({c}_{m1},{c}_{m2}\right)$, ${D}^{\beta}u$ is the shortened notation for the Caputo fractional derivative. The initial conditions are

$$u\left(0,s\right)\text{=}{u}_{0}\left(s\right)$$

The boundary conditions are determined by the driving moments,

$$u\left(t,0\right)\text{=0;}{u}_{s}\left(t,l\right)\text{=}\tau \left(t,l\right)\text{(3)}$$

Where, $\tau \text{=}{\left[{\tau}_{1},{\tau}_{2}\right]}^{T}$. Technical constraints impose the sector constraints on the driving torques,

$${\alpha}^{\ast}u\text{}\le \text{}\tau \text{}\le \text{}{\beta}^{\ast}u\text{(4)}$$

The nonlinear term $h\left(u\right)$ is determined by the gravitational component and satisfies the constraints [14,15].

$$\Vert uh\left(u\right)\Vert \text{}\le \text{}\chi {u}^{T}u\text{(5)}$$

Where, $\chi $ is a positive constant. The system output is evaluated as a weighted average

$$y\left(t\right)\text{=}{\displaystyle \underset{0}{\overset{1}{\int}}\psi \left(s\right)}u\left(t,s\right)ds\text{(6)}$$

Where $y\text{=}\left[{y}_{1},{y}_{2}\right]{}^{T}$ and $\psi \left(s\right)$ is an eigenfunction of the operator $\frac{\partial {}^{2}}{\partial s{}^{2}}$. The aim of the control problem is to ensure the movement to the target position represented by the curve ${C}_{t}$,

$${u}_{t}\left(s\right)\text{={}u\left(s\right):s\in {C}_{t}\}\text{(7)}$$

The evolution error will be

$$\epsilon \left(t\right)\text{=}{\displaystyle \underset{0}{\overset{1}{\int}}\psi \left(s\right)}\left({u}_{t}\left(s\right)\text{-}u\left(t,s\right)\right)ds\text{(8)}$$

The dynamic model (2)-(3) becomes

$${I}_{\rho}{\ddot{\epsilon}}_{i}\left(t\right)\text{=-}{C}_{v}D{}^{\beta}\epsilon {}_{i}\left(t\right)\text{-}\alpha \text{}{\epsilon}_{i}\left(t\right)\text{+}EI{\tau}_{i}\left(t\right)\text{-}h{}^{\ast},i\text{=1,2(9)}$$

Where

$$a\text{=}\left(\frac{{\pi}^{2}}{4l{}^{2}}EI\text{+}{C}_{e}\right)\text{(10)}$$

$${\epsilon}^{i}\left(0\right)\text{=}{\epsilon}_{0}^{i}\text{(11)}$$

The fractional order model will be [18-21] $\left(for\text{}\beta \text{=0}\text{.5}\right)$.

$${\epsilon}_{1}^{i}\left(t\right)\text{=}{e}_{1}^{i},D{}^{\frac{1}{2}}e{}_{1}^{i}\left(t\right)\text{=}{e}_{2}^{i}\left(t\right),D{}^{\frac{1}{2}}e{}_{2}^{i}\left(t\right)\text{=}{e}_{3}^{i}\left(t\right),\mathrm{....},i\text{=1,2(12)}$$

The following fractional order model can be inferred from (9)-(12) (for simplification of notations the index i is eliminated starting from here),

$$D{}^{\frac{1}{2}}e\left(t\right)\text{=}Ae\left(t\right)\text{+}b\tau \left(t\right)\text{+}d{m}^{\ast}\left(t\right)\text{+}{h}^{\ast}\text{(13)}$$

$$z\left(t\right)\text{=}{c}^{T}e\left(t\right)\text{(14)}$$

Where

$$e\left(t\right)\text{=}{\left[{e}_{1}\left(t\right){e}_{2}\left(t\right){e}_{3}\left(t\right){e}_{4}\left(t\right)\right]}^{T}\text{(15)}$$

$$A\text{=}\left[\begin{array}{cccc}0& 1& 0& 0\\ 0& 0& 1& 0\\ 0& 0& 0& 1\\ -\text{}\frac{1}{{I}_{\rho}}& \left(\frac{{\pi}^{2}}{4l{}^{2}}EI\text{+}{C}_{e}\right)& -\text{}\frac{1}{{I}_{\rho}}{C}_{v}0& 0\end{array}\right]\text{(16)}$$

$$c\text{=}{\left[\begin{array}{cccc}1& 0& 0& 0\end{array}\right]}^{T};b\text{=}{\left[\begin{array}{cccc}0& 0& 0& \frac{1}{{I}_{\rho}}\end{array}\right]}^{T}\text{(17)}$$

$${h}^{*}\text{=}{\left[\begin{array}{cccc}0& 0& 0& \frac{1}{{I}_{\rho}}{h}^{\ast}\end{array}\right]}^{T}\text{(18)}$$

## Grasping Control System

The control of the grasping operation is a hybrid control system, position-force control that comprises two phases. In the first phase, the free movement of the arm, without contact with the object, is controlled to the target position ${u}_{t}$ (Figure 4), in the second phase; the grasping operation is performed, the arm exerting on the object the needed compression force. The control system is shown in Figure 5.

Two cases will be analyzed: The control system for a free movement to a desired position and the control system for the control of the grasping compression forces.

### Grasping position control

Dynamic model is rewritten in terms of error (13)-(18) as

$${D}^{\frac{1}{2}}e\left(t\right)\text{=}Ae\left(t\right)\text{+}b\tau \left(t\right)\text{+}{h}^{\ast}\left(e\right)\text{(19)}$$

Consider the control law

$$\tau \left(t\right)\text{=-}ke\left(t\right)\text{(20)}$$

The force constraints (4) determine sector type restrictions on the gain k,

$$k\in \left[{k}_{\alpha},{k}_{\beta}\right]\text{(21)}$$

Where ${k}_{\alpha},{k}_{\beta}$ define the new sector parameters.

*Algorithm 1*. The dynamic system (19), with the control variable (20), (21) is asymptotically stable if:

a) ${A}^{\ast}\text{=}\left(A\text{-}R\right)$ is Hurwitz stable, where $$R\text{=}{R}^{T}\text{0}$$ (22)

b) $Re\left(G\left(j\omega \right)\right)\text{}\ge \text{-}\sqrt{\sigma},\forall \omega \in \left(\text{-}\infty \text{,+}\infty \right)\text{(23)}$

c) ${\alpha}^{\ast}\text{=-2}{\lambda}_{\mathrm{max}\left(PR\right)}\text{0(24)}$

Where the following notations were used: $G\left(j\omega \right)\text{=}C{}^{T}\left(j\omega I\text{-}\left({A}^{\ast}\text{+}\chi I\right)\right){}^{-1}b$, $Q\text{=}qq{}^{T}$, *P* is a solution of Lyapunov equation, *I* is the unit matrix, $\varrho \text{}=\text{}\Vert \left(q\text{+}k\sqrt{\sigma}c\right){\left(q\text{+}k\sqrt{\sigma}c\right)}^{T}\Vert $ and $\sigma $ is a constant that satisfies the condition

$$k\sigma \text{}\le \text{1(25)}$$

Where, from (21) yields

$$\sigma \in \left[{\sigma}_{\beta},{\sigma}_{\alpha}\right]\text{(26)}$$

The frequency inequality (23) can be interpreted in terms of the constraints (21), (26) as [22],

$$Re\left[\frac{{\sigma}_{\beta}\text{+}G\left(j\omega \right)}{{\sigma}_{\alpha}\text{+}G\left(j\omega \right)}\right]\text{0,}\forall \omega \in \left(-\infty ,\text{+}\infty \right)\text{(27)}$$

Consider that $C\left({\sigma}_{\alpha},{\sigma}_{\beta}\right)$ to be the circle in the complex plane (Figure 6) defined by the points $-\text{}{\sigma}_{\alpha}\text{+}j\text{0}$ and $-\text{}{\sigma}_{\beta}\text{+}j\text{0}$. In this case, the condition (26) can be defined as: all points of plot of $G\left(j\omega \right)$ must be outside the circle $C\left({\sigma}_{\alpha},{\sigma}_{\beta}\right)$, or as the plot of $G\left(j\omega \right)$ does not intersect $C\left({\sigma}_{\alpha},{\sigma}_{\beta}\right)$ (Figure 6).

Because $G\left(j\omega \right)$ defines the frequency characteristics of the arm and $C\left({\sigma}_{\alpha},{\sigma}_{\beta}\right)$ represents the input constraints, the stability condition can be rewritten as:" the both domains must not intersect".

$$G\cap C\left({\sigma}_{\alpha},{\sigma}_{\beta}\right)\text{=}\varphi \text{(28)}$$

The relation (28) can be considered as an extension of the Popov's Circle Criterion for this class of fractional order models [22].

*Convergence Algorithm Proof.* The following Lyapunov function is used

$$V\left(e\right)\text{=}{e}^{T}Pe\text{(29)}$$

Where, $P\text{=}{P}^{T}\text{0}$. This Lyapunov function verifies asymptotic stability conditions [19,23,24]. From [21- 26], the fractional derivative of (29) becomes

$${D}^{\beta}V\left(e\right)\text{}\le \text{}\left({D}^{\beta}{e}^{T}\right)Pe\text{+}{e}^{T}P\left({D}^{\beta}e\right)\text{(30)}$$

By evaluating this inequality along the solutions of (19), results

$${D}^{\beta}V\left(e\right)\text{}\le \text{}{e}^{T}\left({A}^{T}P\text{+}PA\right)e\text{+}2{e}^{T}Pb\tau \text{(31)}$$

Considering (22), the relation (31) becomes

$${D}^{\beta}V\left(e\right)\text{}\le \text{}{e}^{T}\left(\left(A\text{-}R\right){}^{T}P\text{+}P\left(A\text{-}R\right)\right)e\text{+}2{e}^{T}Pb\tau \text{+}2{e}^{T}PRe\text{(32)}$$

From (19), (21), (25), this inequality is rewritten as

$${D}^{\beta}V\left(e\right)\text{}\le \text{}{e}^{T}\left(\left(A\text{-}R\right){}^{T}P\text{+}P\left(A\text{-}R\right)\right)e\text{+}2e\left(Pb\text{-}\frac{1}{2}c\right)\tau \text{-}\sigma \tau {}^{2}\text{}\text{+}2{e}^{T}PRe\text{(33)}$$

Applying Yakubovici-Kalman-Popov (YKP) Lemma [22], yields

$${e}^{T}\left(\left(A\text{-}R\right){}^{T}P\text{+}P\left(A\text{-}R\right)\right)e\text{=-}qq{}^{T}\text{}\text{(34)}$$

$$Pb\text{-}\frac{1}{2}c\text{=}\sqrt{\sigma}q\text{(35)}$$

Substituting this result into (33), after simple calculations, considering the control law (19), it follows that

$${D}^{\beta}V\left(e\right)\text{}\le \text{-}{e}^{T}\left(q\text{+}k\sqrt{\sigma}c\right)\left(q\text{+}k\sqrt{\sigma}c\right){}^{T}e\text{+}{\rho}_{PR}{e}^{T}e\text{(36)}$$

$${D}^{\beta}V\left(e\right)\text{}\le \text{-}\alpha {}^{\ast}\Vert e\Vert {}^{2}\text{}\text{}QED\text{(37)}$$

### Grasping force control

Consider that the stress-strain $\left(f,\varrho \right)$ relation of the Kevin-Voigt model (1) can be approximated by

$\left(f,\varrho \right)\text{}\approx \text{}\left(f,\Delta u\right)\text{(38)}$

Assume that $f\left(t,s\right)$ is the measured distributed force, $f\text{=}{\left[{f}_{1,}{f}_{2}\right]}^{T}$ , ${f}_{d}\left(s\right)$ is the desired, needed grasping force. The weighted average of the error is evaluated from (6) as

$${e}_{f}\left(t\right)\text{=}{\displaystyle \underset{0}{\overset{l}{\int}}\psi \left(s\right)}\left({f}_{d}\left(s\right)\text{-}f\left(t,s\right)\right)\text{}ds\text{(39)}$$

The force dynamics can be approximated from (1), (38), (40) using the same procedure as developed in (8),

$${\ddot{e}}_{fi}\text{=}{a}_{f}{e}_{fi}\text{+}{b}_{f}{\tau}_{i}\text{+}{h}^{\ast},i\text{=1,2(40)}$$

$${e}_{f}\text{=}{\left[{e}_{f1},{e}_{f2}\right]}^{T}\text{(41)}$$

Where ${a}_{f},{b}_{f}$ are the equivalent coefficients obtained from (8).

*Algorithm 2.* The force control system is asymptotically stable with the following PD control law

$${\tau}_{i}\left(t\right)\text{=-}{k}_{f1i}{e}_{fi}\left(t\right)\text{-}{k}_{f2i}{\dot{e}}_{fi},\text{}i=1,2\text{(42)}$$

Where ${k}_{f1i},{k}_{f2i}$ are positive gains of the controller that satisfy the conditions:

$${a}_{f}{p}_{3}\text{+}{k}_{f1i}{b}_{f}{p}_{2}\text{-}\chi \text{0(43)}$$

$${k}_{f2i}b{p}_{2}\text{-}{p}_{1}\text{}0\text{,}i=1,2\text{(44)}$$

Where ${p}_{1},{p}_{2},{p}_{3}$ are positive constants ${p}_{1}\text{}{p}_{3},{p}_{2}\text{}{p}_{3}$. that verify

$${p}_{1}\text{-}{a}_{f}{p}_{2}\text{-}{k}_{f1i}{b}_{f}{p}_{2}\text{-}{k}_{f2i}{b}_{f}{p}_{3}\text{+}{p}_{2}\chi =0\text{(45)}$$

Proof: The proof is simple and is obtained by Lyapunov conventional techniques [22,27].

## Numerical Simulations

Consider a grasping control problem for a 2-D motion in a horizontal plane (Figure 7). The following parameters characterize the mechanical system: the length of the arm segment $l=0.15m$, the rotational inertial density ${I}_{\rho}=0.032\text{}kg.m{}^{2}$ , the bending stiffness $EI=1.2\text{}N.m{}^{-1}$, equivalent Kelvin-Voigt model coefficients: ${c}_{v}\text{=0}\text{.26}Nm/rad$, ${c}_{e}\text{=2}\text{.1}Nm\text{}ra{d}^{-1}$. The fractional order model (12)-(15) will be

$$A\text{=}\left[\begin{array}{cccc}0& 1& 0& 0\\ 0& 0& 1& 0\\ 0& 0& 0& 1\\ -\text{249}& -\text{56}\text{.8}& 0& 0\end{array}\right]\text{(46)}$$

$$b\text{=}{\left[\begin{array}{cccc}0& 0& 0& 130\end{array}\right]}^{T};\text{}c\text{=}{\left[\begin{array}{cccc}1& 0& 0& 0\end{array}\right]}^{T}\text{(47)}$$

From (6), the constraints of the compression forces $m\left(s\right)$ are assumed as

$${\alpha}_{m}\text{=}100\text{}Nm,\text{}{\beta}_{m}\text{}=\text{}5000\text{}Nm\text{(48)}$$

The system is fractional order stable but is not Hurwitz stable [28-30]. The matrix $R\text{}=\text{}diag\left(4.5\text{}4.5\text{}4.5\text{}4.5\right)$ is selected to satisfy the condition (22) and the behaviour of the nonlinear model is characterized by $\chi \text{=2}$. A control law (20) is applied with ${k}^{i}\text{}=\text{}4$, ${k}_{\alpha}\text{}=\text{}2$, ${k}_{\beta}\text{}=\text{5}$, ${\sigma}_{\alpha}\text{}=\text{}0.5$, ${\sigma}_{\beta}\text{}=\text{}0.2$. The condition (24) is verified for ${\lambda}_{\mathrm{max}}{}_{\left(Q\right)\text{}}\text{}=\text{}14.5,{\lambda}_{\mathrm{max}}{}_{\left(P\right)\text{}}\text{}=\text{}0.725,\Vert PR\Vert \text{}=\text{}2.47,{\rho}_{R}\text{}=\text{}2.7$, MATLAB/SIMULINK and techniques based on the Mittag-Leffler functions are used for the simulation [31,32]. A weighting eigenfunction $\psi \left(s\right)\text{}=\text{}sin\left(\frac{\pi x}{2l}\right)$ is used. Initial position is ${u}_{0}\left(s\right)\text{}=\text{}\frac{\pi}{2},s\in \left[0,0.15\right]$ and the target position is defined as

$${u}_{t}\left(s\right)\text{}=\text{}s\frac{\pi}{6}\text{(49)}$$

In Figure 8 are shown the error trajectories for position ${e}_{1}\left(t\right)\text{}=\text{}e\left(t\right)$ and velocity ${e}_{3}\left(t\right)\text{}=\text{}\dot{e}\left(t\right)$. The good quality of motion can be remarked from this imagine.

The grasping force control is obtained for ${f}_{d}\left(s\right)\text{}=\text{}150\text{}Nm,s\in \left[0,0.15\right]$ with a spatial weighting average ${f}_{d}\text{}\sim \text{24}Nm$. A control law (42) with ${k}_{f1}\text{}=\text{}18,{k}_{f2}\text{}=\text{}0.7$ is used. The trajectory of the force error ${e}_{f}\left(t\right)$ (39) is shown in Figure 9.

## Discussion

The dynamic model of the continuum arm is associated with the light weight model, for which the gravitational component has a small influence. If this component increases, the stability conditions are getting worse.

In Figure 10 are shown two plots of the linear and nonlinear continuum arm models, respectively. We remark that the nonlinear component moves the plot of $\overline{G}\left(j\omega \right)$ to the left, toward the circle *C* which leads to a decrease in stability. For high value of nonlinearity coefficient $\chi $, when $\overline{G}\left(j\omega \right)$ intersects the circle *C*, the system becomes unstable.

## Conclusions

The paper studies the grasping control by coiling for a class of grippers with continuum arms. First, the fractional order model of continuum arm is proposed. The viscoelastic components are assimilated by the fractional Kevin-Voigt model and the fractional-order dynamics are inferred. A hybrid control technique with two control loops for position and force control, respectively, is proposed. The Lyapunov method for position control is applied. Sector-type constraints for input torque are implemented and frequential conditions that ensure asymptotic stability, similar to the Popov Circle Criterion for classical systems, are inferred. A conventional PD controller is proposed for the force control. Numerical simulations illustrate the performances of the control system.