1 Introduction

In the closest form of human-robot interaction, robots and humans complete common tasks next to each other at the same time [3]. At such workplaces, accidental contacts such as clamping or unconstrained collisions constitute a risk of human injury that technical measures cannot eliminate entirely, since there is always the possibility of a technical failure or foreseeable misuse. As a metric and reference for acceptable risks, ISO/TS 15066 defines biomechanical force and pressure limits for multiple body parts. The technical risk mitigation measures of the robot must reliably prevent these limits from being exceeded [7].

During the validation measurement, the velocity of the robot must be reduced until the measured contact forces and pressures do not exceed the applicable limit values [12]. The proper use of the PFMD requires mounting it on a rigid structure that avoids any undesired displacements of the PFMD. This requirement, however, does not reflect the conditions of an unconstrained contact in which the human body part hit by the robot can move freely.

Falco et al. and Oberer-Treitz presented moveable PFMDs that are capable of simulating the dynamics of free collisions [5, 10]. They attached a PFMD to a linear guide rail and adhered additional weight, so that the total mass of the movable parts reflect exactly the mass of the human body part under test. The US American de-facto standard RIA T R15.806-2018 incorporated the concept of a free-moving PFMD and specifies that it should be used to evaluate free collisions instead of a fixed PFMD. However, a movable PFMD has the major limitation that it can only measure collisions within the horizontal plane.

This article presents a method that allows to convert the results recorded with a fixed PFMD into values that reflect the results taken with a free-moving PFMD. Validation experiments with three collaborative robots (UR3e, UR10e, Doosan M0607) and two different PFMDs confirm that our conversion method provides correct and sufficiently accurate results.

2 Technical Approach

The development of the proposed conversion method can be demonstrated by the analysis of an accurate and simplified impact model. The accurate model in Fig. 1 (top) reproduces the impact behavior of a compliant robot that collides with a biofidelic PFMD. The collaborative robot is modeled as a two mass oscillator [4, 6]. The first mass represents the contribution of the robot’s drive inertia \(m_D\) and the second the mass of its links \(m_L\) , while \(c_T\) describes the joints elasticity. The other part of the model covers the PFMD. The parameters \(c_D\) and \(d_D\) represent the damping behavior and elasticity of the damping material of the PFMD. The damper is attached to a bumper of mass \(m_B\). A mechanical spring of stiffness \(c_S\) absorbs the impact energy from the bumper and transmits the generated force to a load cell that records the impact force \(F_I(t)\). To simulate a properly fixed PFMD the housing is considered to be immovable. To simulate an unconstrained collision the load cell is attached to a body of mass \(m_M\) with static friction \(F_R\).

Fig. 1.
figure 1

Simulation and conversion model of a fixed and free measurement device

To convert the maximum contact force \(\hat{F}_C\) from a fixed PFMD into the maximum contact force \(\hat{F}_F\) of a movable PFMD, the complexity of the simulation model must be reduced. This ultimately leads to the model shown in Fig. 1 (bottom). This model assumes that \(m_D\), \(m_L\), and \(c_T\) are included in an apparent robot mass \(m_R^*\). Due to its simplified structure, the comparison of the energies stored in elasticities and masses of both devices (fixed and free) lead to the following conversion

$$\begin{aligned} \hat{F}_{F} = \hat{F}_{C} \sqrt{\frac{m_H}{m_H + m_R^*}} \ . \end{aligned}$$
(1)

The expression reflects the relative difference between the maximum force of a free \(\hat{F}_F\) and a clamped collision \(\hat{F}_C\). In contrast to \(m_H\), the apparent robot mass \(m_R^*\) is unknown, but can be estimated with the momentum of the impact

$$\begin{aligned} p = m_R^* v_I = \int _{t_0}^{t_{max}} F_C(t) dt\ , \end{aligned}$$
(2)

where \(v_I\) is the impact velocity, \(t_0\) the time of initial contact and \(t_{max}\) the time at which the clamped measured impact force \(F_C(t)\) reaches \(\hat{F}_C\).

3 Model-Based Results

A simulation with the accurate model (Fig. 1, on the top) was performed to study the conversion method. Therefore, the applied model parameters were precisely adjusted to the robot and PFMD, that we later used in the experimental tests for model validation. The robot masses \(m_L\) and \(m_D\) have been calculated with the method presented by Khatib [8, 9]. This method projects the mass matrix of the robot manipulator \(\mathbf {M} ^{+}\) (incl. the contribution of drives and links) to a directional point mass \(m_{R}\) that equals the sum of \(m_L\) and \(m_D\)

$$\begin{aligned} m_{R} = m_L + m_D \ . \end{aligned}$$
(3)

In this context, \(m_L\) denotes the directional point mass that can be calculated with the mass matrix \(\mathbf {M} \), which does not contain contributions from the joint inertias. Then, \(m_D\) can easily be calculated by subtracting \(m_L\) from \(m_R\). The robot’s effective stiffness \(c_T\) was obtained in a similar way, but based on the stiffness matrix \(\mathbf {K} _q\) instead of the \(\mathbf {M} \) [2, 11]. The compliant configuration of the PFMD used for the simulation corresponded to the Back of the Hand (ID 25 in ISO/TS 15066). A spring of \({75\,\mathrm{\text {N}\text {m}^{-1}\text {m}}}\) combined with a damping material of shore hardness 70 (SH70) simulates the biomechanical characteristics of this specific body part [12]. Both stiffness parameters, \(c_D\) for the damping material and \(c_S\) for the spring, were previously measured with a testing system and the same impactors that were used in the experiments. The damping parameter \(d_D\) was estimated in impact experiments at various collision speeds. Also, the friction \(F_R\) of the free PFMD was determined by measurements. Figure 2 shows the impact force over impact velocity obtained from the simulation with a PFMD for fixed and free conditions. The solid line in yellow indicates the results for the fixed PFMD. The results for the free PFMD (solid lines with red diamonds and blue boxes for \({6\,\mathrm{\text {kg}}}\) and \({22\,\mathrm{\text {kg}}}\)) clearly show that the impact force decreases in case \(m_H\) can move freely. The dashed lines in red and blue are the results obtained with the conversion method according to Eq. (1). For both masses (\(m_H = {6\,\mathrm{\text {kg}}}\) and \({22\,\mathrm{\text {kg}}}\)) the curves confirm a good match with the simulation. The intersection points of all lines with the black horizontal line indicate the maximum allowable robot velocity when assuming a force limit of \({280\,\mathrm{\text {N}}}\) (here: Back of the Hand). For the given scenario, the maximum allowable robot velocity for the fixed PFMD is \({271\,\mathrm{\text {m}\text {m}\text {s}^{-1}}}\). In contrast, the simulation of a free collision with a body part of \({22\,\mathrm{\text {kg}}}\) (which is far beyond the expected weight of an outreached hand) enables an increase of the safe velocity by \({62\,\mathrm{\%}}\). For the same configuration, the conversion method estimates an increase of \({60\,\mathrm{\%}}\). The considerable low difference between both results demonstrates that the proposed conversion method obtains sufficiently accurate results.

Fig. 2.
figure 2

Result of the simulated impact force compared to the predicted force of the conversion method for an UR10e robot with a payload of \({66\,\mathrm{\%}}\)

4 Experiments

As the analysis of the simulation data has shown, the mass of the free moving Pressure and Force Measurement Device (PFMD) significantly affects the maximum impact force. Therefore, it was necessary to gather results from collision experiments with a freely moving PFMD of different total masses. Figure 3 shows the test setup that was developed for the experiment. It is based on the concept of a freely movable PFMD as illustrated in RIA T R15.806-2018.

The PFMD was positioned on a movable carrier that was placed on low-friction linear guides with a static friction force of approximately \({2\,\mathrm{\text {N}}}\). Tests with the simulation model revealed that a friction force within this range has a negligible influence on the measurement result. A large industrial robot was used to hold the assembly consisting of guides and the PFMD. Additional weights were mounted on the moveable PFMD to represent different body masses. For the measurement of the impact force, two different PFMDs with a similar design were used as listed in Table 1. The collisions were executed by three collaborative robots (UR3e, UR10e, Doosan M0607) that were equipped with a cylindrical impactor of \({50\,\mathrm{\text {m}\text {m}}}\) in diameter as a tool simulant. The robots were tested in the joint configuration listed in Table 2. In the test with free collisions, the total mass of the freely moving PFMD ranged from 4 to \({22\,\mathrm{\text {kg}}}\) (see Table 3.). For the tests with a fixed PFMD the movable carrier and the system’s ground plate were connected. For a precise validation of the conversion method, the experiments aimed to obtain accurate results for any given combination of damping material and spring. Therefore, every test was performed with combinations that correspond to three different body parts, namely the Abdomen (10), the Sternum (8) and the Back of the Hand (25) (IDs according to ISO/TS 15066). Furthermore, three different payloads were tested with the UR10e in order to exclude possible negative influences of the robot’s payload on the accuracy of the conversion method.

Table 1. Characteristics of the biofidelic PFMDs used for the experimental validation
Fig. 3.
figure 3

Experimental test-setup to measure the impact force of a collaborative robot with a movable biofidelic measurement device

Table 2. Joint configurations of the robots at the impact position
Table 3. Test configurations of the different robots

During the experiments the collaborative robot executes a linear movement at different velocities (starting from 50 and increasing to \({650\,\mathrm{\text {m}\text {m}\text {s}^{-1}}}\)) perpendicular towards the PFMD. The robot’s starting position was chosen so that the desired velocity was ensured to be reached before colliding with the PFMD. The safety functions of the robot monitored the external forces acting on the robot. They triggered an emergency stop, if the magnitude of these forces exceeded the configured threshold (see Table 3).

Fig. 4.
figure 4

Comparison of maximum forces from tests with a fixed and moveable PFMD with those obtained with the conversion method (UR10e with Payload \({6.9\,\mathrm{\text {kg}}}\), Measurement Device: KOLROBOT)

5 Experimental Results

The analysis of the experimental results compares the maximum impact forces of the moveable PFMD with those that were obtained with the conversion method. For the conversion method, the impact force recorded by the fixed PFMD was used. Figure 4 shows one sample of the obtained experimental results. Each of the three diagrams illustrates the results measured with different configurations of the examined PFMD. The markers indicate mean values of the measurement results from tests with the fixed and moveable PFMD. The corresponding mean standard deviation over all measurements shown in Fig. 4 was \(\pm {3.4\,\mathrm{\text {N}}}\) for the fixed and \(\pm {2.8\,\mathrm{\text {N}}}\) for the moveable measurement device. The solid lines represent the results obtained with the conversion method. The results confirm a linear relationship between the impact force and the impact velocity. Moreover, the results demonstrate that the impact forces of a moveable PFMD are significantly lower than those of a fixed PFMD. This effect is enhanced by decreasing the total mass of the human body \(m_H\). Furthermore, the forces obtained with the conversion method are slightly higher than those measured with the movable PFMD. In conclusion, the conversion method achieves almost the same accuracy in the results from experimental tests as from simulations.

6 Discussion

A sensitivity analysis regarding the input parameters of the conversion method is performed to study the influences of uncertainties on the overall accuracy. Of particular interest is the following part of Eq. (1)

$$\begin{aligned} c = \sqrt{\frac{m_H}{m_H + m_R^*}} \ . \end{aligned}$$
(4)

that calculates a factor \(\le 1\) to convert \(\hat{F}_C\) into \(\hat{F}_F\). To determine the sensitivity with respect to the input parameters the partial derivatives for \(m_H\) and \(m_R^*\) must be calculated

$$\begin{aligned} \begin{aligned} \frac{\partial c}{\partial m_H}&= \frac{m_R^*}{2(m_H+m_R^*)^2 \sqrt{\frac{m_H}{m_H + m_R^*}}} \\ \frac{\partial c}{\partial m_R^*}&= -\frac{m_H}{2(m_H+m_R^*)^2 \sqrt{\frac{m_H}{m_H + m_R^*}}} \ . \end{aligned} \end{aligned}$$
(5)

The summation of the partial derivatives gives the total error

$$\begin{aligned} \Delta c = \frac{\partial c}{\partial m_H} \Delta m_H + \frac{\partial c}{\partial m_R^*} \Delta m_R^* \ . \end{aligned}$$
(6)

Figure 5 shows the courses of sensitivity for each parameter. It indicates that the collision masses \(m_H\) and \(m_R^*\) have a non-linear effect on the conversion result. The absolute magnitude of the sensitivity for \(m_H\) and \(m_R^*\) decreases with increasing mass. However, the sensitivity of \(m_R^*\) with respect to c is negative denoting that an underestimation of \(m_R^*\) will result in a higher scaling factor or higher impact force (\(\Delta m_R^* < 0\) leads to \(\Delta c > 0\) and then to \(\Delta \hat{F}_F > 0\)). Since \(m_H\) is assumed to be given, the estimation of \(m_R^*\) must be analyzed in more depth. The partial derivatives of Eq. (2) with respect to the input parameters momentum p and impact velocity v are

$$\begin{aligned} \begin{aligned} \frac{\partial m_R^*}{\partial v}&= -\frac{p}{v^2} \\ \frac{\partial m_R^*}{\partial p}&= \frac{1}{v} \ , \end{aligned} \end{aligned}$$
(7)
Fig. 5.
figure 5

Sensitivity analysis of the conversion with respect to the input parameters human mass \(m_H\) and apparent mass of the robot \(m_R^*\)

that further give the following expression for the total error

$$\begin{aligned} \Delta m_R^* = \frac{\partial m_R^*}{\partial v} \Delta v + \frac{\partial m_R^*}{\partial p} \Delta p \ . \end{aligned}$$
(8)

Equation (7) indicates a proportional influence of p and a negative non-linear one of v. Since \(\frac{\partial m_R^*}{\partial v}\) describes a hyperbolic function, the absolute magnitude of the sensitivity of \(m_R^*\) decreases with increasing v. This behavior indicates that the higher the impact velocity, the better the estimation accuracy. To study whether the sensitivity has a significant influence on the conversion method, the following relative error

$$\begin{aligned} \Delta x_{rel} = \frac{\hat{F}_{conv}}{\hat{F}_{meas}}-1 \end{aligned}$$
(9)

and absolute error

$$\begin{aligned} \Delta x_{abs} = \hat{F}_{conv} - \hat{F}_{meas} \end{aligned}$$
(10)

is introduced. Both can be calculated from the converted force \(\hat{F}_{conv}\), obtained with the conversion method, and the impact force \(\hat{F}_{meas}\) measured with a movable PFMD. The green bars in Fig. 6 display the distribution of the errors in a histogram. It includes all data from the experimental tests, which incorporate more than 1000 individual results. The error distribution shows that \({91\,\mathrm{\%}}\) of the impact forces obtained with the conversion method exceed the forces measured in tests with the movable PFMD. The mean errors are \(\Delta x_{rel}={13.4\,\mathrm{\%}}\) and \(\Delta x_{abs}={12.2\,\mathrm{\text {N}}}\). Considering this tendency, it can be concluded that the conversion method yields to slightly increased impact forces. Consequently, the conversion method can be considered as conservative as higher forces ultimately result in lower but safer robot velocities. The observations confirm that the parameters with positive influence on \(\Delta c\) dominate. Since \(\hat{F}_C\) and \(m_H\) can be considered as accurately known or given, the conservative tendency of the converted values must be caused by a systematic underestimation of \(m_R^*\). Consequently, the technique to estimate \(m_R^*\) results in too small values, which derive from \(\Delta v > 0\) or \(\Delta p < 0\) originating from Eq. (8). We assume that the calculation of p tends to yield smaller values, since the exact time of initial contact \(t_0\) is difficult to detect. Considering this result, it must be highlighted that the conversion method can underestimate the impact forces even though the probability of this effect is very low with \({10\,\mathrm{\%}}\). An evaluation of the dataset shows that a safety factor of \(S=1.1\) multiplied with \(\hat{F}_{conv}\) is sufficient to reduce the probability of underestimations to less than \({1\,\mathrm{\%}}\). This safety factor additionally leads to lower velocities of about \({9\,\mathrm{\%}}\) for the conversion method.

Fig. 6.
figure 6

Distribution of the relative and absolute error of the conversion method for different approaches to calculate the apparent mass \(m_R^*\) (the green bars use Eq. (2), the red bars Eq. (11) to calculate \(m_R^*\); the yellow line indicates \({100\,\mathrm{\%}}\) accuracy; each error on the right hand-side of the yellow line indicates that the conversion method yields to conservative and, thus, safer results)

7 Simplification of the Conversion Method

The error analysis has shown that underestimations of p lead to conservative estimates and lower robot velocities. This finding can be used to further simplify the estimation of the apparent robot mass \(m_R^*\). Especially for end-users who do not have direct access to the measurement data or the knowledge to apply Eq. (2), it might be difficult to obtain p. An approximation of the integral can be obtained when considering \(p_{app}\) as the area of a triangle

$$\begin{aligned} p_{app} = \frac{1}{2}(t_{max}-t_0)\hat{F}_C \ , \end{aligned}$$
(11)

where \(\hat{F}_C\) is the maximum contact force at time \(t_{max}\). Figure 7 compares the triangle approximation to the exact integral. The area of the triangle always corresponds to a slightly smaller momentum \(\Delta p\). According to Eq. (7), a negative momentum error \(\Delta p < 0\) leads to an underestimation of \(m_R^*\), consequently to higher collision forces \(\hat{F}_F\) (see Eq. (5)) and ultimatley to lower, but safe robot velocities. The simplified conversion method that considers \(p_{app}\) was applied to all experimental results to calculate \(m_R^*\). The errors compared to the measurements with the movable PFMDs are shown as red bars in Fig. 6. As expected, the simplified conversion method provides a higher relative \(\Delta x_{rel}={15.9\,\mathrm{\%}}\) and absolute mean error \(\Delta x_{abs}={15.1\,\mathrm{\text {N}}}\), whereas the shape of the error distribution remains almost identical.

Fig. 7.
figure 7

Approximation of the momentum using the area of a triangle (Force signal of a UR10e robot with a payload of \({66\,\mathrm{\%}}\) at an impact velocity of \({150\,\mathrm{\text {m}\text {m}\text {s}^{-1}}}\) colliding with a fixed measurement device with Spring \({10\,\mathrm{\text {N}\text {m}^{-1}\text {m}}}\) and Damper SH10)

8 Conclusion and Outlook

This article presents a new method to convert impact forces recorded with a fixed Pressure and Force Measurement Device (PFMD) into such that would occur in a free collision. The technique has been validated with three different collaborative robots and two different PFMDs. The results obtained from the simulation and experiments show a significant increase of the allowable robot velocities. The findings of this study have a decisive contribution to increase the overall productivity of collaborative robots without compromising their ability to reduce the risk of human injuries. Another contribution of the conversion method to the robot community is the opportunity to evaluate unconstrained human-robot collisions of any direction even in confined spaces. Both advantages outweigh the slightly lower, but more conservative allowable robot velocity that the conversion method provides.

The moveable PFMD intends to represent a free collision between a robot and a part of the human body. One limitation is that it cannot reproduce the exact human behavior since it only has one Degree of Freedom (DoF), while the human body is a multi-body system with several DoFs. However, the short impact duration of free collisions (\({30\dots 200\,\mathrm{\text {m}\text {s}}}\)) and the elasticity of the human joints suggest that the affected human body part initially behaves like a one DoF system of specific apparent mass.

Further work will focus on the estimation of the robots velocity, which is necessary to calculate the apparent mass \(m_R^*\). In fact, the translational velocity of the robot must be assumed as unknown, especially if the robot executes Point-to-Point movements or if the collision point is not at the end-effector. However, it is always possible to use the velocity threshold as it is configured in the safety controller, since a safety-rated function will always ensure that the robot is never exceeding this value.