# Exoscarne: Assistive Strategies for an Industrial Meat Cutting System Based on Physical Human-Robot Interaction

Harsh Maithani, Juan Antonio Corrales Ramon, Laurent Lequievre, Youcef Mezouar, Matthieu Alric

<figure><img src="image_1.png" /></figure>

To cite this version:

Harsh Maithani, Juan Antonio Corrales Ramon, Laurent Lequievre, Youcef Mezouar, Matthieu Alric.

Exoscarne: Assistive Strategies for an Industrial Meat Cutting System Based on Physical Human-Robot Interaction. Applied Sciences, 2021, 11 (9), pp.3907. [hal-03209657](https://uca.hal.science/hal-03209657)

HAL Id: hal-03209657

https://uca.hal.science/hal-03209657

Submitted on 27 Apr 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Article

# Exoscarne: Assistive Strategies for an Industrial Meat Cutting System Based on Physical Human-Robot Interaction

Harsh Maithani <fnref n="1" />*, Juan Antonio Corrales Ramon <fnref n="2" />, Laurent Lequievre <fnref n="1" />, Youcef Mezouar <fnref n="1" />
and Matthieu Alric <fnref n="3" />

**Abstract:** Musculoskeletal disorders of the wrist are common in the meat industry. A proof of concept of a physical human-robot interaction (pHRI)-based assistive strategy for an industrial meat cutting system is demonstrated which can be transferred to an exoskeleton later. We discuss how a robot can assist a human in pHRI, specifically in the context of an industrial project i.e for the meat cutting industry. We developed an impedance control-based system that enables a KUKA LWR robot to provide assistive forces to a professional butcher while simultaneously allowing motion of the knife (tool) in all degrees of freedom. We developed two assistive strategies — a force amplification strategy and an intent prediction strategy - and integrated them into an impedance controller.

**Keywords:** physical human-robot interaction; assistive robotics; collaborative robots

<figure><img src="image_4.png" /></figure>

check for
updates

**Citation:** Maithani, H.; Corrales, J.A.; Lequievre, L.; Mezouar, Y.; Alric, M. Exoscarne: Assistive Strategies for an Industrial Meat Cutting System Based on Physical Human-Robot Interaction. Appl. Sci. **2021**, *1*, 0. https://doi.org/

Academic Editor: Emanuele Carpanzano

Received: 28 March 2021

Accepted: 20 April 2021

Published:

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_5.png" /></figure>

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

## 1. Introduction and Motivation

The agri-food industry, and particularly the meat industry, is one of the most dangerous industries when it comes to employee safety. Among the various occupations, slaughtering, cutting and meat processing operations require specific dexterity to handle sharp tools or dangerous machines. However, they also require physical strength to carry heavy loads such as pallets or carcasses, quarters or muscles of meat, or to perform deboning work tasks and cutting quarters into pieces of meat. Similarly, they require performing repetitive movements or working in a cold refrigerated room and humid environment. In fact, accidents at work are common and can occur at any time. Thus, the rate of these accidents and their frequency are among the highest, among all professions combined. In addition to work-related accidents, musculoskeletal disorders (MSDs) accounted for almost 91% of occupational illness cases in 2019 [1], with 842490 days of temporary interruption of work for all sectors. A study carried out in the Brittany region of France showed that the agri-food industries are the most risky sectors in terms of MSDs [1]. For instance, at the French level, 30% of the declared MSDs are recorded in the meat sector.

Not only is this a major problem for employees, it is also a big problem for companies and society at large. In 2019, compensation for MSDs generated two billion euros in fees (social security estimate) in France, with an average of more than 21,000 euros for each MSD stoppage, not counting the daily allowances [2]. In some companies with 10 to 20% absenteeism, MSDs disrupt production and generate additional costs. In the meat sector, MSDs involve personnel at all stages. Thus, operators involved in meat cutting, represented by boners, parers and slicers, as well as those located at manufacturing and packaging stations are affected [3]. MSDs affecting the wrist, hand and fingers represent approximately 50% of all MSDs in the meat industry.

<footnote n="1">¹ Université Clermont Auvergne, CNRS, Clermont Auvergne INP, Institut Pascal, F-63000 Clermont-Ferrand, France; harshmaithani09@gmail.com (H.M.); laurent.lequievre@uca.fr (L.L.); youcef.mezouar@sigma-clermont.fr (Y.M.)</footnote>

<footnote n="2">² Centro Singular de Investigación en Tecnoloxías Intelixentes (CiTIUS), Universidade de Santiago de Compostela, 15782 Santiago de Compostela, Spain; juanantonio.corrales@usc.es (J.A.C.R)</footnote>

<footnote n="3">³ ADIV, ZAC des Gravanches, 10 Rue Jacqueline Auriol, 63100 Clermont-Ferrand, France; matthieu.alric@adiv.fr (M.A.)</footnote>

<footnote>* Correspondence: harshmaithani09@gmail.com</footnote>
The arduousness of tasks related to physical effort, the repeatability of movements and the agri-food environment (cold, humidity and hygiene) encountered in the meat sector deters the recruitment of young people, and ultimately few people, trained initially for this sector, remain there. This difficulty in recruiting or retaining young people leads to a significant shortage of manpower in these sectors. Technological evolution of certain workstations which would allow on the one hand a reduction of the arduousness, opening to the women certain activities which were until now reserved for the men, and on the other hand a revaluation of the trades, could change the outlook of the individuals in this sector of activity, give a positive image and promote its attractiveness.

Faced with these findings, the meat sector could soon integrate assistive and cobotic equipments to encourage companies to improve the quality of life at work of their employees. Cobotics, also known as collaborative robotics, is a technology that uses robotics, mechanics, electronics, and cognitive science to assist humans in their tasks. For the meat sector, the advantages of having cobotics are numerous:

1. Reduction of occupational risks by reducing the arduousness of the operators' tasks and improving their working conditions.

2. Overcome a major shortage of skilled labour by reducing the lack of image and attractiveness, particularly among young people (hardship in low-value jobs, etc.), which will help maintain or even develop jobs in the sector.

3. Improve the competitiveness of companies by reducing the direct and indirect costs of work stoppages and by increasing productivity.

4. Improve the safety of products by reducing the direct handling of products by operators or by integrating cleaning systems (e.g., sterilization of tools online between each operation).

In this paper, we propose the development of a proof-of-concept assistive strategy implemented through a collaborative robot in meat cutting tasks in order to reduce the musculoskeletal disorders on the wrist of human operators working in the meat industry. The developed impedance control strategy enables a KUKA LWR robot to provide assistive forces to a professional butcher while simultaneously allowing motion of the knife (tool) in all degrees of freedom. Previous robotic systems for autonomous meat handling [4] required one or several robots for performing very specific meat cutting operations. For instance, the ARMS system [5-7] was based on the separation of beef shoulder muscles, the GRIBBOT system [8] was applied to chicken breast fillet harvesting while the DEXDEB system [9,10] was useful for ham deboning. Therefore these robotic systems could not be reused for other meat handling tasks since the quality of their cut was not enough for other types of meat pieces. The new proposed robotic system (called Exoscarne) works on the principle of pHRI (physical Human-Robot Interaction) to solve this lack of generalization. From one side, the expertise of the skilled butcher is kept since the cutting trajectory is defined by the human operator, who is holding the tool at the same time as the robot. From the other side, the robot carries the load of the tool and increases the cutting force when touching the meat so that the effort applied by the human is smaller. This system results in a greater flexibility for different meat cutting tasks and can adapt itself to on-the-fly decisions made by the user.

Therefore, our new Exoscarne system, which will be described in the next sections (see Section 3 for its software components and Section 4 for its hardware components), is able to:

1. Compensate the weight of the cutting tool.

2. Permit the user to move the tool in all 6 degrees of freedom.

3. Permit impedance shaping according to the operation being performed by the user.

4. Provide assistive forces during the meat cutting operations as per the user's convenience.

5. Allow the user to perform the operations autonomously.
## 2. Background on Robot Assistance

One of the primary motivations of using robots for pHRI is their ability to share physical loads with their human partners. When the physical load of an object is shared or when the object is manipulated for a task, a natural division of effort between the human and the robot occurs. For example, load sharing can be for human-robot cooperative manipulation [11] or during rehabilitation [12]. Mortl et al. [13] proposed effort sharing policies for load sharing of an object by a human and multiple robots.

Some researchers focused on the larger question of selection of an 'assistance strategy' for a pHRI task. Dumora et al. [14] considered large object manipulation tasks in pHRI and proposed a library of robot assistances. Medina et al. [15] proposed a dynamic strategy selection between model-based and model-free strategies. The strategy selection is based on the concept of disagreement between the human and the robot, which in turn depends on the interaction force.

In the literature, the only example of load sharing of tool for a cooperative pHRI task, similar to meat cutting, was in [16] in which the author used a robot for assistive welding by supporting the weight of the welding equipment. In fact, most existing IADs, "Intelligent Assist devices" (i.e., active cobotics systems for human assistance) [17], are used in the automotive industry for the quasi-static collaborative transportation of heavy loads (e.g., motors, doors ...). However, these solutions are not suitable for the meat industry [18] since they can only assist through specific directions in the work-space [19] (while meat cutting requires complex 6D trajectories), they do not integrate safety solutions for handling dangerous tools (such as the knife for meat cutting) and they do not take into account the important dynamic non-linear effects of the meat cutting operations [20]. To the best of our knowledge, there is no prior pHRI related work which uses a robot for assisting a human for cutting meat by taking into account not only a classical force amplification strategy (such as in [20]) but also an intent prediction module in order to reduce the final forces to be applied by the human operator.

## 3. Methodology

### 3.1. pHRI Assistive Strategies

The main interaction controllers for pHRI are impedance and admittance control (see Section 3.2 for their mathematical definition). Controller stability issues are common with admittance control [21,22]. When a human holds the tool at the end-effector, it results in a coupled system that can lose stability if the human operator stiffens his arm muscles, leading to robot vibrations. Hence for this task we chose impedance control. Investigating the stability issues of admittance control is a field in itself and there are several heuristic methods in the literature [20]. Impedance control was possible as the robot had a torque sensor in each joint, enabling torque control. Admittance control is preferred with robots that have only position control, by using an external FT sensor. As the task involved motion in the cartesian space hence we used the cartesian impedance controller that is explained later. We devised two assistive strategies:

1. Force amplification strategy

In this strategy we amplify the forces applied by the user on the knife's handle (see Section 4 for the experimental setup), detected by the FT sensor and input to our control scheme. The control diagram is shown in Figure 1. The appropriate parameter $\eta$ has to be determined, as explained in Section 5.3.

2. Intent prediction strategy

In this strategy we predict the forces to be applied by the user using RNN-LSTM networks explained in Section 3.3. The control diagram is shown in Figure 2.
<figure><img src="image_6.png" /><figcaption>Figure 1. Impedance controller diagram with amplification module for the force amplification strategy.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>Figure 2. Impedance controller diagram with intent prediction module for intent prediction strategy.</figcaption></figure>

Both Figures 1 and 2 are identical except for the module of their respective assistive strategies. In the force amplification strategy, the amplification module amplifies the human user's force input $F_h$ at each time step through robot assistance $\eta F_h$. In the intent prediction strategy, a trained RNN-LSTM network takes the human user's force input $F_h$ at each time step to anticipate the user's input for the next time step $F_{pred}$ and provides this force via robot assistance. In both cases the user can haptically sense this assistance being provided by the robot.

## 3.2. Impedance Control

The forward kinematics of a robotic manipulator is written as [23]:

$$x(t) = f(q) \qquad (1)$$

where $x(t) \in \mathbb{R}^n$ and $q \in \mathbb{R}^n$ are the pose (i.e., position/orientation) of the end-effector in the Cartesian space and the joint angle coordinates in the joint space, respectively. Differential kinematics is obtained by deriving (1) with respect to time:

$$\dot{x}(t) = J(q)\dot{q} \qquad (2)$$

where $J(q) \in \mathbb{R}^{n \times n}$ is the Jacobian matrix. Differentiating (2) again results in the acceleration of the end-effector:

$$\ddot{x}(t) = \dot{J}(q)\dot{q} + J(q)\ddot{q} \qquad (3)$$

The robot arm dynamics in the joint space is given by:

$$M(q)\ddot{q} + C(q, \dot{q})\dot{q} + G(q) = \tau - J^T(q)F \qquad (4)$$
where $M(q) \in \mathbb{R}^{n \times n}$ is the symmetric positive-definite inertia matrix; $C(q, \dot{q}) \dot{q} \in \mathbb{R}^n$ is the Coriolis and Centrifugal forces; $G(q) \in \mathbb{R}^n$ is the gravitational force; $\tau \in \mathbb{R}^n$ is the vector of control input; $F \in \mathbb{R}^n$ denotes the net force exerted by the robot on the environment at the end-effector ($F = F_r - F_{ext}$). Otherwise, $F_r \in \mathbb{R}^n$ is the force exerted by the robot on the environment at the end-effector while $F_{ext} \in \mathbb{R}^n$ is the external force exerted by the environment on the robot at the end-effector. In our pHRI task, the environment is the human, specifically the human hand that is in contact with the robot. The robot dynamics can be written in the Cartesian space as:

$$M_x(q)\ddot{x} + C_x(q, \dot{q})\dot{x} + G_x(q) = u - F \quad (5)$$

where

$$\begin{aligned}
M_x(q) &= J^{-T}(q)M(q)J^{-1}(q), \\
C_x(q, \dot{q}) &= J^{-T}(q)(C(q, \dot{q}) - M(q)J^{-1}(q)\dot{J}(q))J^{-1}(q), \\
G_x(q) &= J^{-T}(q)G(q), \quad u = J^{-T}(q)\tau, \quad F = F_r - F_{ext}
\end{aligned}$$

As explained in Section 3.1, the most common interaction controllers for pHRI are admittance and impedance control. Both controllers are based on a target impedance model for the robot and they only differ in terms of input and output (see Figures 3 and 4). Therefore, the target robot impedance model can be represented as a mass-damper-spring system:

$$M_d(\dot{x}_d - \ddot{x}) + B_d(\dot{x}_d - \dot{x}) + K_d(x_d - x) = F_d = F_r \quad (6)$$

where $M_d, B_d, K_d$ are the virtual inertia, damping and stiffness of the robot, respectively. $F_d$ is the desired force and $x_d$ can be interpreted as the rest position of this virtual mass-damper-spring system. For pHRI tasks in which the human touches the robot, the human limb (arm+hand) can also be modeled as a mass-damper-spring system:

$$M_h \ddot{x}_h + B_h \dot{x}_h + K_h(x_h - x_{hd}) = F_h \quad (7)$$

where $M_h, B_h, K_h$ are the limb inertia, damping and stiffness respectively and $F_h$ is the force applied by the human to the robot. The limb impedance values are not fixed and depend from person to person, as well as on the task being carried out. $x_h$ is the position of the human wrist in the robot frame and $x_{hd}$ can be interpreted as the desired target position. Discussion on the limb impedance parameters or their calculation are not in the purview of this paper. In fact, the investigation of the limb impedance parameters would have been necessitated if an admittance controller was used, to tackle the potential controller stability issues. Since we have used a torque-based impedance controller for our experiments, it does not face such issues. While real time measurements of limb impedance values could be used for interpreting the intent of the human operator during the operation, this would have required attaching EMG sensors on the arms of the human operator, such as in [24], thereby reducing the practicality of this research for an industrial meat cutting scenario. In our experiments the human intent is interpreted via the force measurements read by the FT sensor, while staying comfortable and intuitive for the human operator.

When no assistance strategy is required and we want the robot to freely follow the motion of the human, we can set $K_d = 0$ or $x = x_d$. This “direct teaching mode” is essentially based on the elimination of the spring component of the robot impedance model (“minimum impedance” strategy in Section 5.4) and avoids any restoring forces. If the human holds the tool rigidly then we can assume that the forces are transmitted completely to the robot. If the hand is close enough to the end-effector, we can assume that their positions, velocities and forces are equivalent: $x_d = x_{hd}, \dot{x}_h = \dot{x}$ and $\ddot{x}_h = \ddot{x}$.

In admittance control (popularly called position-based impedance control) the input is the force applied by the environment and the output is displacement/velocity as shown in Figure 3. In pHRI tasks, the force is applied by the human on the tool at the end-effector ($F_{ext} = F_h$) and it is measured by the FT sensor. The difference between the desired force
for the robot ($F_d = F_r$) and this external force is applied as input to the impedance model in (6) in order to obtain the desired position for the end-effector $x_d$ as output:

$$M_d(\ddot{x}_d - \ddot{x}) + B_d(\dot{x}_d - \dot{x}) + K_d(x_d - x) = F_d - F_h \quad (8)$$

On the contrary, the impedance control (torque-based impedance control) has as input the end-effector displacement/velocity and the desired force as output, as shown in Figure 4. In this case, the desired robot impedance model is the same as in (6) but the difference between the desired robot force and the external/human force is transformed into joint torques for the robot control.

<figure><img src="image_8.png" /><figcaption>Figure 3. Admittance control (position-based impedance control).</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>Figure 4. Torque-based impedance control.</figcaption></figure>

### 3.3. Long Short-Term Memory Model

For the intent prediction module, we use RNN-LSTM units [25]. RNNs (Recurrent Neural Networks) are based on processing sequential data, especially temporal data as they have an internal memory. In fact, they are useful for making predictions using time-series data [23]. RNNs can be improved by using what are called LSTM units in order to solve the vanishing gradient problem (i.e., gradients tend to disappear in RNNs when backpropagating errors in too long sequences). In fact, each LSTM unit has a special structure composed by 3 gates to control what information to keep and what to forget so that more stable errors are backpropagated (see Figure 5):

1. An input gate (i).

2. An output gate (o).

3. A forget gate (f).

As a result, RNNs with LSTM units are able to learn long-term dependencies within data sequences that were not possible only with RNNs. Given an input sequence $\bar{x}_1, \bar{x}_2, ..., \bar{x}_t$ the LSTM unit maps the input sequence to a sequence of hidden states $h_1, h_2, ..., h_t$ (which are also the outputs) by passing information through a combination of gates (see Figure 6): The Input gate (for updating the cell) is:

$$i_t = \sigma_g(W_i \bar{x}_t + R_i h_{t-1} + b_i) \quad (9)$$

The Forget gate (for resetting the cell/forgetting) is:

$$f_t = \sigma_g(W_f \bar{x}_t + R_f h_{t-1} + b_f) \quad (10)$$
The Cell candidate (for adding information to the cell) is:

$$g_t = \sigma_c(W_g \bar{x}_t + R_g h_{t-1} + b_g) \qquad (11)$$

where $\sigma_c$ is the state activation function (here it is the hyperbolic tangent function: $\sigma_c = tanh(\bar{x})$).

The Output gate is:

$$o_t = \sigma_g(W_o \bar{x}_t + R_o h_{t-1} + b_o) \qquad (12)$$

where $\sigma_g$ is the gate activation (here it is the sigmoid function: $\sigma(x) = (1 + e^{-\bar{x}})^{-1}$)

The Memory Cell state at timestep t is:

$$c_t = f_t \odot c_{t-1} + i_t \odot g_t \qquad (13)$$

Here $\odot$ is the Hadamard product (element-wise multiplication of vectors). The memory cell selectively retains information from previous timesteps by controlling what to remember via the forget gate $f_t$.

<figure><img src="image_10.png" /><figcaption>Figure 5. A single LSTM unit.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>Figure 6. An unrolled Recurrent Neural Network with LSTM units.</figcaption></figure>

The Hidden state (also called Output state) at time step t is:

$$h_t = o_t \odot \sigma_c(c_t) \qquad (14)$$

The hidden state is passed as input to the next timestep and thus, it is possible to stack numerous LSTM units (see Figure 6). $W_i, W_f, W_g$ are the learnable input weights; $R_i, R_f, R_g$ are the learnable recurrent weights and $b_i, b_f, b_g$ are the learnable bias. By using memory cells and hidden states, LSTM units are able to retain information. The sigmoid function is a good activation function for the 3 gates (In, Out and Forget) since it outputs a value between 0 and 1. However for the memory cell, the values should be able to increase or decrease (which is not possible with the sigmoid activation function whose output is
always non-negative). Therefore, the hyperbolic tangent function (tanh) is used as the activation function for the memory cell.

# 4. Experimental Setup

## 4.1. KUKA LWR Robot

The KUKA LWR is a Torque-controlled Flexible Robot with 7 degrees of freedom. The design and control concepts of the robot have been discussed in [26,27]. There are two control modes available. The joint position control is implemented at a frequency of 3 kHz (decentralized control) [28]. Inverse kinematics and the cartesian impedance control mode run at a frequency of 1 kHz. The KUKA LWR has a torque sensor at each joint that enables torque control and impedance control. It also has motor side position sensors, as well as link side position sensors. Due to friction it is difficult for robots to implement torque control only with motor current commands [28].

The default coordinate system of the KUKA LWR is Right Handed System for which $det(R) = 1$. The dimensions of the robot can be retrieved from the official manufacturer documentation. Using these dimensions we can obtain the DH parameters (refer Table 1). To simulate the robot, an alternate form of representation called URDF is shown in Table 2 and its corresponding 3D model is shown in Figure 7. The Unified Robotic Description Format (URDF) is an XML file format used in ROS to describe all elements of a robot. The robot interacts with the external PC through the Fast Research Interface [29] through three modes—(a) Joint position control (b) Joint impedance control mode and (c) Cartesian impedance control mode.

<figure><img src="image_12.png" /><figcaption>Figure 7. Exoscarne 3D simulated model with the KUKA LWR arm in Home configuration.</figcaption></figure>

<table><caption>Table 1. DH parameters of KUKA LWR 4+ robot.</caption><thead><tr><th></th><th>Joints</th><th>d<sub>i</sub>(m)</th><th>q<sub>i</sub>(rad)</th><th>a<sub>i</sub></th><th>α<sub>i</sub>(rad)</th><th>q<sub>min</sub></th><th>q<sub>max</sub></th><th>τ<sub>max</sub>(Nm)</th></tr></thead><tbody><tr><td>J1</td><td>A1</td><td>0.3105</td><td>q<sub>1</sub></td><td>0</td><td>π/2</td><td>-170</td><td>170</td><td>176</td></tr><tr><td>J2</td><td>A2</td><td>0</td><td>q<sub>2</sub></td><td>0</td><td>-π/2</td><td>-120</td><td>120</td><td>176</td></tr><tr><td>J3</td><td>E1</td><td>0.4</td><td>q<sub>3</sub></td><td>0</td><td>-π/2</td><td>-170</td><td>170</td><td>100</td></tr><tr><td>J4</td><td>A3</td><td>0</td><td>q<sub>4</sub></td><td>0</td><td>π/2</td><td>-120</td><td>120</td><td>100</td></tr><tr><td>J5</td><td>A4</td><td>0.39</td><td>q<sub>5</sub></td><td>0</td><td>π/2</td><td>-170</td><td>170</td><td>100</td></tr><tr><td>J6</td><td>A5</td><td>0</td><td>q<sub>6</sub></td><td>0</td><td>-π/2</td><td>-120</td><td>120</td><td>38</td></tr><tr><td>J7</td><td>A6</td><td>0.078</td><td>q<sub>7</sub></td><td>0</td><td>0</td><td>-170</td><td>170</td><td>38</td></tr></tbody></table>
<table><caption>**Table 2.** URDF description of KUKA LWR 4+ robot.</caption><thead><tr><th>Joints</th><th>x(m)</th><th>y(m)</th><th>z(m)</th><th>r</th><th>p</th><th>y</th><th>Axis</th></tr></thead><tbody><tr><td>A1</td><td>0</td><td>0</td><td>0.11</td><td>0</td><td>0</td><td>0</td><td>[0 0 1]</td><td>+z axis</td></tr><tr><td>A2</td><td>0</td><td>0</td><td>0.2005</td><td>0</td><td>0</td><td>0</td><td>[0 −1 0]</td><td>−y axis</td></tr><tr><td>E1</td><td>0</td><td>0</td><td>0.2</td><td>0</td><td>0</td><td>0</td><td>[0 0 1]</td><td>+z axis</td></tr><tr><td>A3</td><td>0</td><td>0</td><td>0.2</td><td>0</td><td>0</td><td>0</td><td>[0 1 0]</td><td>+y axis</td></tr><tr><td>A4</td><td>0</td><td>0</td><td>0.2</td><td>0</td><td>0</td><td>0</td><td>[0 0 1]</td><td>+z axis</td></tr><tr><td>A5</td><td>0</td><td>0</td><td>0.19</td><td>0</td><td>0</td><td>0</td><td>[0 −1 0]</td><td>−y axis</td></tr><tr><td>A6</td><td>0</td><td>0</td><td>0.078</td><td>0</td><td>0</td><td>0</td><td>[0 0 1]</td><td>+z axis</td></tr></tbody></table>

To fulfill our objectives we used the cartesian impedance control mode available with the KUKA LWR robot. The KUKA manual [30] states that the control law for the cartesian impedance controller is

$$ \tau_{cmd} = J^T(K_c(x_{desired} - x_{current}) + F_{cmd}) + D(d_c) + f_{dynamics}(q, \dot{q}, \ddot{q}) \quad (15) $$

where $q \in \mathbb{R}^n$ is the joint position vector, $K_d$ is the stiffness matrix in the end-effector frame, $D_d$ is the normalized damping parameter in the end-effector frame, $x$ and $x_d$ are the current and the desired pose of the end-effector respectively in the global frame. The translational stiffness $K_x$, $K_y$, $K_z \in [0.01, 5000]$ N/m and rotational stiffness $K_{Az}$, $K_{By}$, $K_{Cx} \in [0.01, 300]$ N/m-rad.

The KUKA LWR has an inbuilt external tool calibration functionality. Using this feature the robot can account for external tool dynamics as well, thereby enabling gravity compensation for the tool.

In this work we used the KUKA Fast Research Interface (FRI [29]), ROS [31], Kinematics and Dynamics library (KDL [32]), the MATLAB toolbox by Peter Corke (RCV [33]), and the MATLAB Robotics Toolbox (RTB [34]). For the intent prediction strategy, the high level program was written in MATLAB and connected to the network via ROS (Robot Operating System), by using the MATLAB Robotics Toolbox.

## 4.2. ATI FT Sensors

In order for the robot to provide assistive forces as per user comfort we had to measure how much forces are being applied by the user. For this, we used two 6-axis ATI Gamma force-torque sensors (see Figure 8): one mounted below the joystick (sensor B) and another mounted on the end-effector of the arm (sensor A, see Figure 9). Both FT (Force-Torque) sensors provide a 6-dimensional wrench in the sensor frame at 1000 Hz.

<figure><img src="image_13.png" /><figcaption>**Figure 8.** Photos of the force-torque (FT) sensor.</figcaption></figure>
<figure><img src="image_14.png" /><figcaption>Figure 9. Attaching the FT sensor A to the KUKA LWR 4+ robot end-effector.</figcaption></figure>

The relation between the different frames associated to these sensors (A and B), to the tool (i.e., tool center point—TCP—or tip of the knife) and the robot world frame O are shown in Table 3 and visualized in Figure 10. To use the impedance controller all the forces must be in the same frame. This requires transformation of the sensed forces in the sensor frame to the end-effector frame of the robot. The equation for transformation of forces from one frame to another is:

$$ \begin{bmatrix} A F_A \\ A M_A \end{bmatrix} = \begin{bmatrix} B R & 0 \\ [A t_x]_A^B R & B R \end{bmatrix} \begin{bmatrix} B F_B \\ B M_B \end{bmatrix} \quad (16) $$

where $[t_x] = \begin{bmatrix} 0 & -t_z & t_y \\ t_z & 0 & -t_x \\ -t_y & t_x & 0 \end{bmatrix}$ or

$$ A F_A = B_A T_f B F_B \quad (17) $$

<figure><img src="image_15.png" /><figcaption>Figure 10. Visualization of all frames relating the tool -knife- with the FT sensors A/B.</figcaption></figure>
<table><caption>**Table 3.** Joystick frames when the robot is in Home configuration. Sensor A is always aligned with the robot end-effector frame.</caption><thead><tr><th>Robot World Frame</th><th>Sensor A (FT13855)</th><th>Sensor B (FT13953)</th><th>Relation</th></tr></thead><tbody><tr><td>X<sub>o</sub></td><td>X<sub>A</sub></td><td>-Z<sub>B</sub></td><td>X<sub>o</sub>|| X<sub>A</sub>|| -Z<sub>B</sub></td></tr><tr><td>Y<sub>o</sub></td><td>Y<sub>A</sub></td><td>X<sub>B</sub></td><td>Y<sub>o</sub>|| X<sub>A</sub>|| X<sub>B</sub></td></tr><tr><td>Z<sub>o</sub></td><td>Z<sub>A</sub></td><td>Y<sub>B</sub></td><td>Z<sub>o</sub>|| Z<sub>A</sub>|| Y<sub>B</sub></td></tr></tbody></table>

### 4.3. Allen-Bradley Joystick

The joystick (i.e., the knife handle) is an Allen-Bradley 440J-N enabling switch. It has 8 electrical connections and when the joystick switch is pushed or released, combinations of these 8 connections are activated as shown in Figure 11.

<figure><img src="image_16.png" /><figcaption>**Figure 11.** Allen-Bradley joystick used as the handle of the knife.</figcaption></figure>

The 8 electrical wires from the joystick are connected to an Arduino Uno circuit board which was then connected to a laptop (see Figure 12 for the complete connection diagram). The *rosserial* ROS package and the *ros_lib* library is used to integrate the Arduino Uno with the robot network via ROS.

<figure><img src="image_17.png" /><figcaption>**Figure 12.** Arduino circuit diagram for connecting the joystick to a laptop.</figcaption></figure>
By using this joystick, the user will be able to communicate his intention clearly with regards to the meat cutting operation. See Section 5.2 for a more detailed description of how the user will use this joystick for a safe and unobtrusive physical human-robot interaction.

### 4.4. Meat Cutting Equipment

To fully understand the practical issues in implementing the technical aspects of the project we first performed some meat cutting trials using the robot in the gravity compensation mode. A special holding tool was developed to accommodate a joystick which would be held by the user and the cutting tool. The project involved meat operations on chicken, pork and beef for which special cutting tools were made. The joystick was an industrial joystick which we customized for the project and it allowed a very natural and intuitive user interface. The blades are professional butcher blades customized for the project by machining in the laboratory. A customized pneumatic hanger was developed for holding the pork vertically. Specific cutting tools and holding apparatus for each meat product or carcass to cut were devised as shown in Figure 13 below.

<figure><img src="image_18.png" /><figcaption>(a) Horizontal cutting</figcaption></figure>

<figure><img src="image_19.png" /><figcaption>(b) Vertical cutting</figcaption></figure>

<figure><img src="image_20.png" /><figcaption>(c) Knife</figcaption></figure>

<figure><img src="image_21.png" /><figcaption>(d) Peeler</figcaption></figure>

<figure><img src="image_22.png" /><figcaption>(e) Chicken holder</figcaption></figure>

<figure><img src="image_23.png" /><figcaption>(f) Pork</figcaption></figure>

<caption>Figure 13. Exoscarne equipment.</caption>

## 5. Experiments

The complete network diagram is shown in Figure 14. All the PCs of this network communicate via ROS. The two assistive strategies are implemented in the right laptop, which recovers FT sensor information from the two middle PCs and sends robot commands to the left PC.
<figure><img src="image_24.png" /><figcaption>**Figure 14.** Network diagram of the Exoscarne system. From left to right: control box of the KUKA LWR 4+ robot connected to one PC with FRI/KDL libraries; data acquisition boxes for FT sensors A and B connected to two PCs and an Arduino board for joystick acquisition connected to a laptop.</figcaption></figure>

We performed experiments to address the following questions:

1. Which FT sensor should be used as a source of input for the control scheme?

2. What should the impedance shaping strategy performed through the joystick buttons be?

3. What should the amplification factor be for the force amplification strategy?

4. Comparison of the force amplification strategy and the intent prediction strategy.

### 5.1. *Comparison of FT Sensors*

As shown in Figure 10, there were two FT sensors on the cutting tool. We wanted to confirm if 2 FT sensors provide an advantage over a single FT sensor, as well as determine which FT sensor is better as a source of input for the control scheme.

A piece of thick foam was used as the material to be cut. This foam reproduces the same type of shearing forces as meat cutting. Two experiments were performed. In the first one the robot was commanded to apply forces gradually from 0 to 50 N in the Y-direction of the world frame, which was parallel to the cutting direction of the knife. In the second experiment, a human cut the foam multiple times in the Y-direction, as shown in Figure 15.

Based on the sensor readings from both the sensors, it was observed that for free motion in space (without any cutting), both the sensors record the same readings. However, for cutting motion, sensor A (refer to Figure 10) is blind to the cutting forces. Sensor B recorded forces when the robot cut the foam autonomously, but the recorded forces were feeble and almost negligible to the forces commanded to the robot. When a human cut the foam using impedance controller the sensor B registered adequate forces.
<figure><img src="image_25.png" /><figcaption>Figure 15. Comparison of sensor A and sensor B.</figcaption></figure>

Hence it was determined that it is reasonable to consider that the forces sensed by sensor B (below the joystick) capture the human applied forces on the knife and thus, the user's intention. During actual meat cutting experiments, both the sensors were alternated as sensor inputs for the controller in numerous trials and the user responded that he preferred the system behaviour when sensor B was the sensor input.

## 5.2. Impedance Shaping for Cutting

While cutting meat, the user performs an active meat cutting operation followed by an inactive repositioning of the knife for the next stroke. Furthermore, the user may want to use both of his hands for some other task and as such need the robot to stay still in the last position. The joystick has two buttons (see Figure 10): the first lateral grey button with three possible positions (i.e., a dead-man's switch) and the second upper black button having two states (i.e., a pushbutton). This gives us a total of 6 possible combinations shown in Table 4. The corresponding electrical connections of both buttons are shown in Figure 11.

<table><caption>Table 4. States of joystick buttons.</caption><thead><tr><th rowspan="3">Button 1</th><th colspan="3">Button 2</th></tr><tr><th colspan="2">State 0</th><th>State 1</th></tr><tr><th>Position 0</th><th>Robot is stiff + no amplification</th><th>Robot is stiff + amplification</th></tr></thead><tbody><tr><th>Position 1</th><td colspan="2">Robot is free, but no amplification<br>(used for positioning of knife)</td><td>Robot is free + amplification<br>(used for cutting)</td></tr><tr><th>Position 2</th><td colspan="2">Robot is free, but no amplification<br>(used for positioning of knife)</td><td>Robot is free + amplification<br>(used for cutting)</td></tr></tbody></table>

The button configurations from Table 4 are interpreted in the impedance shaping algorithm as shown in c. We can see that the last two rows are identical. Originally we tried to use position 2 with an impedance relation (i.e $K, D \leftarrow F_h$); however, the user preferred to keep the interface simple and make the position 1 and position 2 identical for operation. This shows that developing user-friendly interface is essential for the adoption of pHRI over human-only or robot-only equipment. The constant values of Table 5 are:

$$K_{min} = [0.1 \ 0.1 \ 0.1 \ 0.1 \ 0.1 \ 0.1],$$

$$K_{max} = [5000 \ 5000 \ 5000 \ 300 \ 300 \ 300],$$

$$D_{min} = [0.01 \ 0.01 \ 0.01 \ 0.01 \ 0.01 \ 0.01],$$

$$D_{max} = [1.0 \ 1.0 \ 1.0 \ 1.0 \ 1.0 \ 1.0]$$
<table><caption>**Table 5.** Impedance shaping using joystick.</caption><thead><tr><th></th><th colspan="2">Button 2</th></tr><tr><th></th><th>State 0</th><th>State 1</th></tr></thead><tbody><tr><th rowspan="3" scope="rowgroup">Button 1</th><th scope="row">Position 0</th><td>K<sub>max</sub>, D<sub>max</sub>, η = 0</td><td>K<sub>max</sub>, D<sub>max</sub>, η</td></tr><tr><th scope="row">Position 1</th><td>K<sub>min</sub>, D<sub>min</sub>, η = 0<br>(used for positioning of knife)</td><td>K<sub>min</sub>, D<sub>min</sub>, η<br>(used for cutting)</td></tr><tr><th scope="row">Position 2</th><td>K<sub>min</sub>, D<sub>min</sub>, η = 0<br>(used for positioning of knife)</td><td>K<sub>min</sub>, D<sub>min</sub>, η<br>(used for cutting)</td></tr></tbody></table>

### 5.3. Tuning the Amplification Factor

In the force amplification strategy we amplify the forces applied by the user on the joystick, detected by the FT sensor and input to our control scheme (see Figure 1). For the robot to provide assistive forces we had to determine the amplification factor $\eta$ for each degree of freedom as shown in Equation (18):

$$F_{cmd} = \begin{bmatrix} e^{F_x} \\ e^{F_y} \\ e^{F_z} \\ e^{F_{Az}} \\ e^{F_{By}} \\ e^{F_{Cx}} \end{bmatrix} = \begin{bmatrix} \eta_x \\ \eta_y \\ \eta_z \\ \eta_{Az} \\ \eta_{By} \\ \eta_{Cx} \end{bmatrix} \begin{bmatrix} e^{F_x} \\ e^{F_y} \\ e^{F_z} \\ e^{\tau_z} \\ e^{\tau_y} \\ e^{\tau_x} \end{bmatrix} \quad (18)$$

While a high amplification factor would ease the load on the user, it could also give the user the perception that he is no longer in control of the operation and reduce his comfort with the system (as was realized during the experiments). This is the reason we cannot have the robot simply apply the highest forces possible.

The tuning of the comfortable amplification factors was a continuous process where several iterations of meat cutting were performed by a professional butcher. It was decided collectively that there would be only two amplification factors: one common $\eta_f$ for forces and one $\eta_\tau$ for torques (i.e., $\eta_f = \eta_x = \eta_y = \eta_z$ and $\eta_\tau = \eta_{Az} = \eta_{By} = \eta_{Cx}$). This iterative tuning process finished when the butcher found that the assistance behavior was comfortable, as explained in Figure 16.

The meat cutting operation involves both the application of forces and torques, as such the force amplification factor $\eta_f$ cannot be determined independently of the torque amplification factor $\eta_\tau$. In earlier experiments the user felt that $\eta_\tau = 3$ was ideal for him and hence to determine $\eta_f$ a series of consecutive experiments were done as shown in Table 6 enabling the user to make a subjective comparison.

At the end of the experiment it was concluded that the user preferred $\eta_f \in [10, 20]$ and $\eta_\tau = 3$. With $\eta_f = 20$ and $\eta_\tau = 4$, the user found the system to be too reactive and he felt he was no longer in control. These 6 experiments (see Figures 17, 19, 21, 23, 25 and 27, respectively) and the temporal evolution of the cutting forces applied by the user after applying these optimal force amplification factors (see Figures 18, 20, 22, 24, 26 and 28, respectively) are shown in the next section.
<figure><img src="image_26.png" /><figcaption>**Figure 16.** Iterative procedure for tuning the two amplification factors for forces and torques.</figcaption></figure>

<table><caption>**Table 6.** Determining the force amplification factor $\eta_f$.</caption><thead><tr><th>S.No.</th><th>Meat</th><th>Part</th><th>Side</th><th>Operation</th><th>FT Sensor Input</th><th>Force Amplification η<sub>f</sub></th><th>Torque Amplification η<sub>τ</sub></th></tr></thead><tbody><tr><td>1</td><td>Pork</td><td>Square</td><td>Right</td><td>Cutting</td><td>B</td><td>20</td><td>3</td></tr><tr><td>2</td><td>Pork</td><td>Square</td><td>Right</td><td>Cutting</td><td>B</td><td>20</td><td>3</td></tr><tr><td>3</td><td>Pork</td><td>Square</td><td>Right</td><td>Cutting</td><td>B</td><td>15</td><td>3</td></tr><tr><td>4</td><td>Pork</td><td>Square</td><td>Right</td><td>Cutting</td><td>B</td><td>15</td><td>3</td></tr><tr><td>5</td><td>Pork</td><td>Square</td><td>Right</td><td>Cutting</td><td>B</td><td>10</td><td>3</td></tr><tr><td>6</td><td>Pork</td><td>Square</td><td>Right</td><td>Cutting</td><td>B</td><td>20</td><td>3</td></tr></tbody></table>
5.4. *Meat Cutting with Force Amplification Strategy*

**Experiment 1—Cobot assisted pork cutting**

A knife mounted on the end-effector of a robot is inherently dangerous and hence the knife was always covered with a sheath when not in use to avoid accidental injury. When experiments were being performed, even if the user could stop the robot at any time by releasing the button, a second person always had his/her hand on the emergency robot switch to disable the robot if something goes wrong. The user also wore a protective gear around his arms. In addition, as a precaution, the robot was wrapped with a plastic sheet to prevent minute meat pieces from entering the internal structure of the robot.

<figure><img src="image_27.png" /><figcaption>**Figure 17.** Experiment 1—cobot assisted pork cutting.</figcaption></figure>

<figure><img src="image_28.png" /><figcaption>**Figure 18.** Forces applied by the user in Experiment 1.</figcaption></figure>

The terminology A cut, B cut etc. are author defined to refer to the iteration of the cutting and not the cutting methods. The images shown in the figures do not represent one single cut in progress, but instead one image of each cut. They can be interpreted as cut 1, cut 2, etc.
Experiment 2—Cobot assisted pork cutting

<figure><img src="image_29.png" /><figcaption>**Figure 19.** Experiment 2—cobot assisted pork cutting.</figcaption></figure>

<figure><img src="image_30.png" /><figcaption>**Figure 20.** Forces applied by the user in Experiment 2.</figcaption></figure>

In experiment 1, it took the butcher (the user) 4 cuts to make 1 slice. Nevertheless, in experiment 2, with the same section but from a different meat, it took him 9 cuts to make 1 slice. This is due to the natural variation in the body composition from one animal meat to another and also how much forces the user wanted to apply for a single cut.
Experiment 3—Cobot assisted pork cutting

<figure><img src="image_31.png" /><figcaption>**Figure 21.** Experiment 3—cobot assisted pork cutting.</figcaption></figure>

<figure><img src="image_32.png" /><figcaption>**Figure 22.** Forces applied by the user in Experiment 3.</figcaption></figure>
Experiment 4—Manual pork cutting

<figure><img src="image_33.png" /><figcaption>Figure 23. Experiment 4—manual pork cutting.</figcaption></figure>

<figure><img src="image_34.png" /><figcaption>Figure 24. Forces applied by the user in Experiment 4.</figcaption></figure>
Experiment 5—Manual pork cutting

<figure><img src="image_35.png" /><figcaption>Figure 25. Experiment 5—manual pork cutting.</figcaption></figure>

<figure><img src="image_36.png" /><figcaption>Figure 26. Forces applied by the user in Experiment 5.</figcaption></figure>
Experiment 6—Manual pork cutting

<figure><img src="image_37.png" /><figcaption>Figure 27. Experiment 6—manual pork cutting.</figcaption></figure>

<figure><img src="image_38.png" /><figcaption>Figure 28. Forces applied by the user in Experiment 6.</figcaption></figure>
Experiment 7—Foam cutting with intent prediction module

For meat cutting experiments with force amplification strategy we had a professional butcher as the user. However for verifying the intent prediction strategy we used a foam block and multiple users in order to perform the training of the LSTM network (see Figure 29 for the corresponding experimental setup). This foam reproduces the same type of shearing forces as meat cutting.

<figure><img src="image_39.png" /><figcaption>Figure 29. Foam cutting along the global y-direction.</figcaption></figure>

As explained earlier, our intent prediction module uses RNN-LSTM units. For each user, we performed sample trials with foam cutting to collect the training dataset. The LSTM network was trained on this dataset to predict force values, similar to the prediction of force values of the Natural Motion dataset (NM-F) in [23].

For each user we took 90 percent of the sample dataset as the training dataset. The entire architecture consisted of 4 layers- an input layer, an LSTM layer, a fully connected layer and a regression layer. We tested the prediction accuracy with combinations of different hyperparameters such as the number of epochs, sequence length and learning rate. However the results were almost the same i.e., a prediction accuracy of 0.3 (root mean square error).

The number of features was 1 as we had only a single variable—force applied. Similarly as only 1 output was expected, the number of responses was 1. The number of hidden units was taken as 200 and the maximum number of epochs was set to 250. The sequence length for the input layer was 25 time steps (0.2 s). For the output layer, we used Stochastic Gradient Descent algorithm with a learning rate of 0.01. The sigmoid function was used as the activation function for the 3 gates—In, Out and Forget in the LSTM units as it outputs a value between 0 and 1. However for the memory cell, the values should be able to increase or decrease which is not possible with the sigmoid function as the output is always non-negative, hence we used the hyperbolic tangent function (tanh) as the activation function for the memory cell.

Figure 30 shows the plot of cutting forces applied by a user with and without the intent prediction module (for 30 cm cutting of the foam). When the intent prediction module was turned off, we had $K_{min}$ and $D_{min}$ only (minimum impedance). With the intent prediction module the user applied only 20 percent of the forces as compared to when the module was turned off.

We also compared the force amplification strategy with the intent prediction strategy with 5 users as shown in Figure 31. For each user sample trials were conducted with force amplification strategy for them to decide which amplification factor they are comfortable
with. All the users stated that the intent prediction module made the cutting more intuitive than the force amplification strategy.

<figure><img src="image_40.png" /><figcaption>**Figure 30.** Cutting forces applied by the user using the intent prediction module for 30 cm cutting of the foam. The blue line is the force applied with *K*<sub>min</sub> and *D*<sub>min</sub> only (i.e. minimum impedance), while the red line is the force applied with *K*<sub>min</sub> and *D*<sub>min</sub> and the intent prediction module.</figcaption></figure>

<figure><img src="image_41.png" /><figcaption>**Figure 31.** Comparison of different strategies with 5 users.</figcaption></figure>

## 6. Conclusions and Future Work

In this paper, we demonstrated a proof of concept of two pHRI-based assistive strategies for an industrial meat cutting system. It was determined that sensor B below the joystick gives better reactivity and hence this sensor should be used as the sensor input. From the previous experiments, it is evident that the forces applied by the user are approximately 30% small with the cobot and the force amplification strategy than with the manual meat cutting operation. Blades with lengths 20 cm and 10 cm were tested, and the one with length 10 cm was adjudged to give better cutting performance, because it was stiffer, with less bending and mechanical compliance.

The cartesian impedance controller runs at a frequency of 1000 Hz, as well as the FT sensor. At this frequency the user found it intuitive and useful to operate the tool for the meat cutting task without any issue. However, it is known that the human central nervous
system operates at a lower frequency than a robot controller, and if he is coupled to the system via a tool it is possible for the user to 'perceive' a loss of control if the system is too reactive, an observation that capped the upper limit of the amplification in the force amplification strategy.

In the foam cutting experiment, it was shown that the intent prediction strategy was better than the force amplification strategy, especially with regards to intuitiveness. In a pHRI system such as this one, not only should the robot provide assistance as a machine, but also the interface should be intuitive, natural and easy to use.

Our contributions are:

1. We followed a systematic methodology to develop a user friendly pHRI controller system for meat cutting.

2. We developed two assistive strategies: a force amplification strategy and an intent prediction strategy.

3. The developed system allowed the user to move the knife in all 6 degrees of freedom.

4. Using impedance shaping, the impedance values were altered between the cutting and non-cutting re-positioning movement of the knife, which was found as comfortable by the user.

For future work, we would like to compare the force amplification strategy and the intent prediction strategy on meat cutting tasks. Furthermore, in the current experiments we had only one professional butcher and it would be interesting to see what are the experimental results with more professional users. The next step of project Exoscarne would involve the design and development of a specific exoskeleton for meat cutting and transferring the controller that was developed in this work.

**Author Contributions:** Conceptualization: H.M., J.A.C.R., Y.M.; methodology: H.M., J.A.C.R., L.L., Y.M., M.A.; software, H.M., L.L.; validation and analysis: H.M., M.A.; writing—original draft preparation: H.M.; writing—review and editing: J.A.C.R., L.L., Y.M., M.A.; supervision: J.A.C.R., Y.M.; funding acquisition: J.A.C.R., Y.M., M.A. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received funding from the French government research program Investissements d'Avenir through the UMTs ACTIA Mécarnéo-AgRobErgo and the project Exoscarne (Call P3A-ICF2A-2I2A, FranceAgriMer) and from the European Union's Horizon 2020 research and innovation programme under grant agreement n° 869855 (SoftManBot project).

**Informed Consent Statement:** Informed consent was obtained from all subjects involved in the study.

**Conflicts of Interest:** The authors declare no conflict of interest.

**Abbreviations**

The following abbreviations are used in this manuscript:

pHRI Physical human-robot interaction

MSD Musculoskeletal disorders

ARMS A multi arms Robotic system for Muscle Separation

RNN Recurrent neural networks

LSTM Long Short Term Memory

ROS Robot Operating System

FT Force-Torque

URDF Unified Robot Description Format

FRI Fast Research Interface

DH Denavit-Hartenberg

LWR Light Weight Robot

KDL Kinematics and Dynamics Library

IAD Intelligent Assist Devices
References

1. *Statistiques Accidents du Travail et Maladies Professionnelles*; Carsat Bretagne: Rennes, France, 2020.

2. *Rapport Annuel 2019, L'Assurance Maladie-Risques Professionnels, Éléments Statistiques et Financiers*; Caisse nationale de l'Assurance Maladie des travailleurs salariés: Paris, France, 2020.

3. *État de Santé des Salariés de la Filière Viande du Régime Agricole en Bretagne*; Institut de Veille Sanitaire: Saint-Maurice, France, 2019.

4. De Medeiros Esper, I.; From, P.J.; Mason, A. Robotisation and intelligent systems in abattoirs. *Trends Food Sci. Technol.* **2021**, *108*, 214–222, doi:10.1016/j.tifs.2020.11.005.

5. Long, P.; Khalil, W.; Martinet, P. Modeling control of a meat-cutting robotic cell. In Proceedings of the 2013 16th International Conference on Advanced Robotics (ICAR), Montevideo, Uruguay, 25–29 November 2013; pp. 1–6, doi:10.1109/ICAR.2013.6766471.

6. Long, P.; Khalil, W.; Martinet, P. Robotic cutting of soft materials using force control image moments. In Proceedings of the 2014 13th International Conference on Control Automation Robotics Vision (ICARCV), Singapore, 10–12 December 2014; pp. 474–479, doi:10.1109/ICARCV.2014.7064351.

7. Long, P.; Khalil, W.; Martinet, P. Force/vision control for robotic cutting of soft materials. In Proceedings of the 2014 IEEE/RSJ International Conference on Intelligent Robots and Systems, Chicago, IL, USA, 14–18 September 2014; pp. 4716–4721, doi:10.1109/IROS.2014.6943233.

8. Misimi, E.; Øye, E.R.; Eilertsen, A.; Mathiassen, J.R.; Åsebø, O.B.; Gjerstad, T.; Buljo, J.; Øystein Skotheim. GRIBBOT—Robotic 3D vision-guided harvesting of chicken fillets. *Comput. Electron. Agric.* **2016**, *121*, 84–100, doi:10.1016/j.compag.2015.11.021.

9. Wei, G.; Stephan, F.; Aminzadeh, V.; Dai, J.S.; Gogu, G. DEXDEB—Application of DEXtrous Robotic Hands for DEBoning Operation. In *Gearing up and Accelerating Cross-Fertilization Between Academic and Industrial Robotics Research in Europe*; Röhrbein, F., Veiga, G., Natale, C., Eds.; Springer International Publishing: Cham, Switzerland, 2014; pp. 217–235.

10. Alric, M.; Stephan, F.; Sabourin, L.; Subrin, K.; Gogu, G.; Mezouar, Y. Robotic solutions for meat cutting and handling. In *European Workshop on Deformable Object Manipulation*; Innorobo: Lyon, France, 2014.

11. Lawitzky, M.; Mörtl, A.; Hirche, S. Load sharing in human-robot cooperative manipulation. In Proceedings of the 19th International Symposium in Robot and Human Interactive Communication, Viareggio, Italy, 13–15 September 2010; pp. 185–191, doi:10.1109/ROMAN.2010.5598627.

12. Reinkensmeyer, D.J.; Wolbrecht, E.; Bobrow, J. A computational model of human-robot load sharing during robot-assisted arm movement training after stroke. In Proceedings of the 2007 29th Annual International Conference of the IEEE Engineering in Medicine and Biology Society, Lyon, France, 22–26 August 2007; pp. 4019–4023.

13. Mörtl, A.; Lawitzky, M.; Kucukyilmaz, A.; Sezgin, M.; Basdogan, C.; Hirche, S. The role of roles: Physical cooperation between humans and robots. *Int. J. Robot. Res.* **2012**, *31*, 1656–1674, doi:10.1177/0278364912455366.

14. Dumora, J.; Geffard, F.; Bidard, C.; Aspragathos, N.A.; Fraisse, P. Robot assistance selection for large object manipulation with a human. In Proceedings of the 2013 IEEE International Conference on Systems, Man, and Cybernetics, Manchester, UK, 13–16 October 2013; pp. 1828–1833.

15. Medina, J.R.; Lawitzky, M.; Molin, A.; Hirche, S. Dynamic strategy selection for physical robotic assistance in partially known tasks. In Proceedings of the 2013 IEEE International Conference on Robotics and Automation, Karlsruhe, Germany, 6–10 May 2013; pp. 1180–1186.

16. Erden, M.S.; Marić, B. Assisting manual welding with robot. *Robot. Comput.-Integr. Manuf.* **2011**, *27*, 818–828. doi: 10.1016/j.rcim.2011.01.003.

17. Colgate, J.E.; Peshkin, M.; Klostermeyer, S.H. Intelligent assist devices in industrial applications: A review. In Proceedings 2003 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS 2003), Las Vegas, NV, USA, 27–31 October 2003; Volume 3, pp. 2516–2521, doi:10.1109/IROS.2003.1249248.

18. Paxman, J.; Liu, D.; Wu, P.; Dissanayake, G. *Cobotics for Meat Processing: An Investigation into Technologies Enabling Robotic Assistance for Workers in the Meat Processing Industry*; Technical Report; University of Technology: Sydney, Australia, 2006.

19. Campeau-Lecours, A.; Otis, M.J.D.; Gosselin, C. Modeling of physical human-robot interaction: Admittance controllers applied to intelligent assist devices with large payload. *Int. J. Adv. Robot. Syst.* **2016**, *13*, doi:10.1177/1729881416658167.

20. Lamy, X.; Collédani, F.; Geffard, F.; Measson, Y.; Morel, G. Overcoming human force amplification limitations in comanipulation tasks with industrial robot. In Proceedings of the 2010 8th World Congress on Intelligent Control and Automation, Jinan, China, 7–9 July 2010; pp. 592–598, doi:10.1109/WCICA.2010.5553839.

21. Duchaine, V.; Gosselin, C.M. Investigation of human-robot interaction stability using Lyapunov theory. In Proceedings of the 2008 IEEE International Conference on Robotics and Automation, Pasadena, CA, USA, 19–23 May 2008; pp. 2189–2194.

22. Dimeas, F.; Aspragathos, N. Online stability in human-robot cooperation with admittance control. *IEEE Trans. Haptics* **2016**, *9*, 267–278.

23. Maithani, H.; Ramon, J.A.C.; Mezouar, Y. Predicting Human Intent for Cooperative Physical Human-Robot Interaction Tasks. In Proceedings of the 2019 IEEE 15th International Conference on Control and Automation (ICCA), Edinburgh, UK, 6–9 July 2019; pp. 1523–1528, doi:10.1109/ICCA.2019.8899490.

24. Grafakos, S.; Dimeas, F.; Aspragathos, N. Variable admittance control in pHRI using EMG-based arm muscles co-activation. In Proceedings of the 2016 IEEE International Conference on Systems, Man, and Cybernetics (SMC), Budapest, Hungary, 9–12 October 2016; pp. 001900–001905.

25. Hochreiter, S.; Schmidhuber, J. Long short-term memory. *Neural Comput.* **1997**, *9*, 1735–1780.
26. Albu-Schäffer, A.; Haddadin, S.; Ott, C.; Stemmer, A.; Wimböck, T.; Hirzinger, G. The DLR lightweight robot: Design and control concepts for robots in human environments. *Ind. Robot. Int. J.* **2007**, *34*, 376–385.

27. Bischoff, R.; Kurth, J.; Schreiber, G.; Koeppe, R.; Albu-Schäffer, A.; Beyer, A.; Eiberger, O.; Haddadin, S.; Stemmer, A.; Grunwald, G.; et al. The KUKA-DLR Lightweight Robot arm-a new reference platform for robotics research and manufacturing. In *Proceedings of the ISR 2010 (41st International Symposium on Robotics) and ROBOTIK 2010 (6th German Conference on Robotics)*, VDE, Munich, Germany, 7–9 June 2010; pp. 1–8.

28. Albu-Schaffer, A.; Hirzinger, G. Cartesian impedance control techniques for torque controlled light-weight robots. In *Proceedings 2002 IEEE International Conference on Robotics and Automation (Cat. No. 02CH37292)*, Washington, DC, USA, 11–15 May 2002; Volume 1, pp. 657–663.

29. Schreiber, G.; Stemmer, A.; Bischoff, R. The fast research interface for the kuka lightweight robot. In *Proceedings of the IEEE Workshop on Innovative Robot Control Architectures for Demanding (Research) Applications How to Modify and Enhance Commercial Controllers (ICRA 2010)*, Anchorage, Alaska, USA, 3-7 May 2010; pp. 15–21.

30. *KUKA LWR Document Manual*; KUKA Roboter GmbH, Augsburg, Germany 2010.

31. Quigley, M.; Conley, K.; Gerkey, B.; Faust, J.; Foote, T.; Leibs, J.; Wheeler, R.; Ng, A.Y. ROS: An open-source Robot Operating System. In *ICRA Workshop on Open Source Software*; IEEE ICRA 2009: Kobe, Japan, 12-17 May 2009; Volume 3, p. 5.

32. Available online: http://docs.ros.org/indigo/api/orocos_kdl/html/ (accessed on 21 April 2021).

33. Corke, P. Robotics, Vision and Control: Fundamental Algorithms in MATLAB® Second Edition, Completely Revised; Springer: Berlin/Heidelberg, Germany, 2017; Volume 118.

34. Available online: https://in.mathworks.com/help/robotics/ (accessed on 21 April 2021).