Abstract
With the increasing importance of collaborative robots in industrial manufacturing, their economic efficiency is becoming more and more important. Today, one approach to prevent collaborative robots from injuring humans is to assure that the robot cannot exceed biomechanical limits in the event of an accidental collision or clamping. The ability of the robot to avoid collision forces beyond the limits must be validated with a biofidelic measurement device that mimics the biomechanical behavior of the human. For reliable use, the measurement devices must be attached to a rigid frame. Consequently, the test setup is solely able to simulate the contact dynamics and biomechanical consequences of a clamping contact in which the human body part cannot move. Free collisions that allow the human to move freely reduce the collision forces, but can only be evaluated with such a measurement device. This technical limitation leads to slower robots and thus a loss of productivity. The study presents a method that increases the efficiency of safe collaborative robots by adding a new validation procedure. The presented method incorporates a model-based conversion of measurements that enables safety experts to validate robots in free collisions. Data from experimental tests with a collaborative robot and a biofidelic measurement device show a good fit of the model-based prediction and thus confirm our approach. This new approach has great potential to increase the productivity of collaborative robots, since our method will allow them to move at faster but still safe speeds.
You have full access to this open access chapter, Download conference paper PDF
Similar content being viewed by others
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\).
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
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
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\)
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.
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.
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).
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)
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
The summation of the partial derivatives gives the total error
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
that further give the following expression for the total error
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
and absolute error
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.
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
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.
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.
Change history
17 June 2021
The original version of the book was inadvertently published without abstracts in chapters 19, 20. The erratum chapter and the book have been updated with the changes.
28 March 2021
The original version of the chapter “A New Conversion Method to Evaluate the Hazard Potential of Collaborative Robots in Free Collisions” was previously published as non-open access. This has now been changed to the copyright holder “The Author(s)” and Open Access under the terms of the Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0/). The chapter and the book have been updated with the change.
References
7th International Conference Safety of Industrial Automated Systems, Montreal, Canada, October 2012
Albu-Schäffer, A., Fischer, M., Schreiber, G., Schoeppe, F., Hirzinger, G.: Soft robotics: what Cartesian stiffness can obtain with passively compliant, uncoupled joints? In: Proceedings of the 2004 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS), Sendai, Japan, February 2004, vol. 4, pp. 3295–3301 (2004). https://doi.org/10.1109/IROS.2004.1389925
Behrens, R., Saenz, J., Vogel, C., Elkmann, N.: Upcoming technologies and fundamentals for safeguarding all forms of human-robot collaboration. In: Proceedings of the 8th International Conference on the Safety of Industrial Automated Systems, pp. 18–23, Königswinter, Germany, November 2015
Bicchi, A., Tonietti, G.: Fast and soft-arm tactics. IEEE Rob. Autom. Mag. 11(2), 22–33 (2004). ISSN 1070-9932. https://doi.org/10.1109/MRA.2004.1310939
Falco, J., Marvel, J., Norcross, R.: Collaborative robotics: measuring blunt force impacts on humans. In: Proceedings of the 7th International Conference on the Safety of Industrial Automated Systems (2012)
Haddadin, S., Krieger, K., Mansfeld, N., Albu-Schäffer, A..: On impact decoupling properties of elastic robots and time optimal velocity maximization on joint level. In: Proceedings, pp. 5089–5096, Vilamoura, Portugal, October 2012. https://doi.org/10.1109/IROS.2012.6385913
Huelke, H.-J., Ottersbach, M.: How to approve collaborating robots - the IFA force pressure measurement system. In: Proceedings of the 7th International Conference on the Safety of Industrial Automated Systems SIA, pp. 5089–5096 (2012)
Khatib, O.: A unified approach for motion and force control of robot manipulators: the operational space formulation. IEEE J. Rob. Autom. 3(1), 43–53 (1987). https://doi.org/10.1109/JRA.1987.1087068
Khatib, O.: Inertial properties in robotic manipulation: an object-level framework. Int. J. Rob. Res. 14(1), 19–36 (1995). ISSN: 0278-3649. https://doi.org/10.1177/027836499501400103
Oberer-Treitz, S.: Abschaetzung der Kollisionsfolgen von Robotern zur Bewertung des sicheren Einsatzes in der Mensch-Roboter-Kooperation. Dissertation, Universitaet Stuttgart (2017)
Salisbury, J.-K.: Active stiffness control of a manipulator in Cartesian coordinates (1980)
Fachbereich Holz und Metall der DGUV. Dguv-information - collaborative robot systems: design of systems with power and force limiting function (2017)
Acknowledgement
Parts of this work have been funded by the ROSSINI-Project, which has received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement No 818087. The biofidelic measurement device PRMS (PILZ Robot Measurement System) used for the experiments was provided by courtesy of PILZ.
Author information
Authors and Affiliations
Corresponding author
Editor information
Editors and Affiliations
Rights and permissions
Open Access This chapter is licensed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license and indicate if changes were made.
The images or other third party material in this chapter are included in the chapter's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the chapter's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.
Copyright information
© 2021 The Author(s)
About this paper
Cite this paper
Herbster, S., Behrens, R., Elkmann, N. (2021). A New Conversion Method to Evaluate the Hazard Potential of Collaborative Robots in Free Collisions. In: Siciliano, B., Laschi, C., Khatib, O. (eds) Experimental Robotics. ISER 2020. Springer Proceedings in Advanced Robotics, vol 19. Springer, Cham. https://doi.org/10.1007/978-3-030-71151-1_20
Download citation
DOI: https://doi.org/10.1007/978-3-030-71151-1_20
Published:
Publisher Name: Springer, Cham
Print ISBN: 978-3-030-71150-4
Online ISBN: 978-3-030-71151-1
eBook Packages: Intelligent Technologies and RoboticsIntelligent Technologies and Robotics (R0)