Investigating frictional contact behavior for soft material robot simulations

The ability to interact safely with the environment is known as one of the major advantages of soft robots (SRs). Due to their low material stiffness, these continuously deformable robots offer inherent flexibility. These advantages make them suitable for application that involve human-robot collaboration in industrial settings as well as medical application such as minimally invasive surgery. To date only few research groups have analyzed the contact and frictional behavior of soft robots. In fact, the contact behavior is often oversimplified or neglected. Motivated by the idea to bridge this gap, this work presents measurements and the resulting coefficient of friction (COF) for silicone rubbers that are widely used in the field of SRs and different contact partners which depend on contact pressure and ambient temperature. From these measurements, a more representative contact model is established and used to more accurately simulate soft material robots’ frictional contact behavior. Moreover the influence of friction and therefore the need to implement frictional behavior is demonstrated for a typical application of a SR.


Introduction
Inspired by nature, soft robots (SRs) are an aspiring field of research.Compared to their rigid counterparts, SRs are defined by their low material and structural stiffness which makes them suitable to interact safely in complex or unknown environments.Particularly, SRs are inherently safe to use in human-robot interaction due to their softness and the associated integrated compliance.The idea behind this bioinspired technology is to mimic specific functionalities such as the flexibility of an octopus or the movement of a worm [1].These advantages make SRs well suited to use in medical applications [2] or in highly constrained environments as for example shown in [3].Soft materials have enabled the development of numerous types of actuators, for example, fluid-driven actuators [4].In the past years, the field of SRs developed rapidly, presenting new approaches of manufacturing, modeling, and controlling these robots.Yet the field has not been able to significantly contribute to industrial applications according to [5].One approach to enable new and meaningful developments which will allow the SRs field to go beyond fundamental research is to establish toolboxes and open source resources which future research can build upon [6][7][8].As the main motivation behind SRs is strongly tied to their favourable environmental contact behavior, understanding the phenomena related to contact and in particular friction is essential.[9] mention contact modeling as a challenge that has to be addressed in future research in the field of soft material robotics.In Fig. 1 three SRs that represent typical fields in soft robotics are shown.While there exist various designs of SRs, they all experience similar contact and frictional behavior when interacting with their respective environments.

On the importance of modeling friction
While the soft robotics community acknowledges the importance of friction in modeling SRs [9,13], the frictional behavior of commonly used materials, e.g.[8,[14][15][16][17], has not been investigated in a systematic way.
In order to clearly demonstrate the influence of modeling friction, the contact behavior is investigated using a toy example.For this purpose, a cube with an edge length of 30 mm and placed at a distance of 1 mm from a rigid plate is modeled in the commercially available software ABAQUS, as shown in Fig. 2 (I-III).The material for the cube is a soft silicon rubber (Sect.4.1.1)that represents the low material stiffness which is typical for SRs.The cube's clamped side is then moved by x=1.7 mm so that it comes into contact with the plate.Afterwards the cube is moved along the plate's longitudinal axis by y=2 cm in a static simulation.For the contact setup, the penalty contact is used for both the normal and tangential contact behavior.Here the COF based on the Fig. 1 The figure gives an overview of typical bioinspired SRs (a) [10] whose functionalities rely on frictional behavior such as crawling (b) [11] or gripping (c) [12].During contact, the COF can vary depending on e.g.temperature, sliding distance or contact pressure.The schematic SR (d) showcases a contact scenario where a friction map, that captures these dependencies is used Fig. 2 Simulation steps for contact simulation of a toy example (I-III).In the first step (A), the soft cube is placed at a distance from a wall.In the second step (B), the cube is brought into contact with the wall by applying a displacement boundary condition as indicated by the arrow in the x direction.In the last step (C), relative motion between the cube and the wall is implemented.The results show the soft cube sliding on a rigid wall with non-frictional contact a-c and frictional contact d-f Coulomb friction model can be defined as described in [18].For this example a COF of 0.5 is chosen to represent the frictional behavior.
In Fig. 2a-f, the results of this simulation are shown for non-frictional contact (a-c) and frictional contact (d-f).The results demonstrate, that the deformation of the cube strongly depends on the frictional behavior.While no tangential forces are transmitted without friction, the tangential forces for the frictional contact lead to a significant deformation of the cube.For frictional contact, the displacement differs from that of the non-frictional case by about 2 mm for a 20 mm displacement, which corresponds to a 10% difference.Also, the total area in contact varies by nearly 20% between the non-frictional and frictional simulations according to the history output variable CNAREA.With regard to applications such as gripping, this result underlines the need to include friction, as this not only enables the gripping mechanism itself, but also significantly affects the contact surface.
These simulations demonstrates the effect of friction on a toy example to highlight the importance of implementing frictional behavior for tasks such as gripping.In the field of SRs, contact modeling is only implemented in some models so far.Some research groups use non-frictional contact for a simple approximation [19,20].According to [21], non-frictional contact can be assumed for medical applications as insufficient data is available and the COF is sufficiently small to be neglected.While this assumption holds for very low COFs, the authors acknowledge that suitable data is not available so far.More recently [22] presented a dynamic model based on discrete differential geometry that accounts for surface roughness as well as frictional contact for SRs.In this paper, the authors measure and implement a single COF.In [13], an overview is given among other things on the modeling and characterization of materials including friction.The overview shows that friction is only investigated and used in few works so far.For different SRs, the assumption of non-frictional contact is no longer possible as their operating principle rely on frictional behavior as demonstrated with the toy example.For example the locomotion of snakelike robots or crawlers as well as various grippers rely on friction [23][24][25][26].
In [27], a framework for dynamic modeling of robots that mimic the function of biological cilia is discussed.The framework is based on a Kirchhoff rod model where the frictional force is incrementally applied if contact occurs to ensure smooth increase.At the same time, depending on the robot and application, friction needs to be minimized in some cases in order to minimize energy loss while obtaining the locomotion [26,28].These conflicting goals demonstrate that each robot's frictional characteristics can be optimized to fulfill a designated task.
Therefore, there is a need for simulations which employ contact models that use reasonable experimentally-informed frictional behavior.Especially when aiming for application that rely on frictional behavior, a good estimation of the coefficient of friction (COF) is required to, e.g., forecast the normal force needed to pick up an object or the actuation needed to realize a specified locomotion.On a broader scale, in [29] contact modeling for flexible systems that undergo large deformation is presented.In contrast to approaches that introduce artificial factors to ensure a smooth increase of contact forces such as the penalty method, this work uses a non-smooth approach.In this case the non-smooth behavior of a system is modeled without the need to smooth the increase of contact forces for numerical reasons.In the case of contact, it can be assured that there is no penetration between contacting bodies.Overall, there are works that consider frictional behavior with an estimated COF, but to the best of the author's knowledge, there exists no work that implements a COF which depends on influencing parameters, such as ambient temperature or contact pressure, in the field of soft robotics.

Contributions
The main contribution of this paper is to provide a friction map containing COFs for different influencing parameter with the aim of encouraging researchers to implement frictional behavior in their SRs models.We present frictional measurements and highlight the importance of friction modeling.
Following the example of the Soft Robotics Materials Database [8], which is a database for unified material constitutive models, and to carry forward the idea of open sourcing knowledge we provide an open source database of COFs, making them accessible to the soft robotics community. 1The COF is sensitive to different factors such as the ambient temperature, the sliding speed between contact partners, the contact pressure as well as the contact partners' materials, and surfaces.Therefore, we present a friction map that lists the COF of different contact partners depending on contact pressure and ambient temperature.The overall aim is to publish a growing friction map that researchers can contribute to in order to provide reasonable COFs for various scenarios.We demonstrate the effectiveness of implementing frictional contact behavior in a finite element analysis (FEA).As a starting point, measurements for the silicone rubbers Ecoflex 00-50 and Dragonskin 30 are presented.According to the manufacturer SmoothOn, Ecoflex 00-50 and Dragonskin 30 have a shore hardness of 50 Shore 00 and 30 Shore A respectively [30,31].

Contact modeling
The COF can depend on the ambient temperature, the sliding speed between contact partners, the contact pressure as well as the contact partners' materials, and surfaces.Therefore, a friction map that lists the COFs between different contact partners depending on contact pressure and ambient temperature is required.The influence of the contact pressure is particularly important for gripping tasks since these tasks require friction.In this case the contact pressure or normal force needs to be controlled in order to lift objects.Therefore this work will focus on the influence of the contact pressure.The influence of the sliding speed is not yet considered in this work.The friction map is then implemented in a finite element simulation to demonstrate the necessity of implementing frictional behavior, using a soft pneumatic actuator (SPA) used for a gripping task as an example.
The frictional behavior of elastomers is very complex including effects such as hysteresis and adhesion and is strongly dependent on factors such as the surface condition and roughness [32,33].For dry contact, a rather simple contact model can be used to describe the frictional behavior between solid bodies [34]: the coulomb friction model.It states that the frictional force F R is proportional to the normal Force F N times a coefficient -the coefficient of friction ( ) which leads to: Frictional behavior can be subdivided into two different cases [34]: static friction: and kinetic friction: The idea of using a coefficient to describe frictional behavior can be enhanced to use a COF that depends on multiple environmental influences such as ambient temperature or sliding distance [35,36].In Fig. 3, a qualitative example measurement of the COF over the sliding distance is shown.The COF rises steadily over measurement time and the related glide path until a maximum is reached.This point marks the maximum frictional forces that can be transmitted under the given condition.Up to this point, we assume gliding in the contact zone.In this case, only static friction occurs which is expressed by the static coefficient of friction S .Thereafter, the COF settles at a constant value that can be lower or equal to the static COF as described for example in [37,38].In this case, the contact partners move relative to one another and the frictional behavior is described by the s t a t i c f r i c t i o n kinetic friction F Fig. 3 Example friction measurement divided into static and kinetic friction Vol.: (0123456789) kinetic coefficient of friction K .The static frictional force F R,stat increases until the maximal force F max is reached which corresponds to the maximal transmittable force.At this point, relative motion can be detected and the contact partners begin to slide relative to each other.In the following measurements, the average COF refers to the average kinetic COF K .
The COF depends on the two specific contact partners and can furthermore depend on: To address these dependencies, a friction map is created which provides a COF for different contact partners depending on the contact pressure and ambient temperature.

Friction measurement
The following section describes the conducted measurements and provides details on the measurement setup.

Experimental setup
All measurements were conducted with a pin-ondisc measurement setup on a MCR 302e rheometer (Anton Paar, Austria) with the T-PID/44 measurement cell (Anton Paar, Austria) as shown in Fig. 4. The frictional behavior of two different contact partners is determined for a translational motion on a circular path at a given ambient temperature.All measurements are conducted in a temperature-controlled chamber to ensure a constant ambient temperature, where the ambient temperature refers to the temperature that is set for a specific measurement.Three sample pins (radius 3 mm, height 6 mm) of the same material are mounted in the measurement cell.A disk (radius 30 mm, height 2 mm) mounted in the rheometer serves as a counterpart to represent the contact partners for the pins.During the measurement, the contact pressure as well as the sliding speed can be adjusted.In preparation for the measurement, the samples are mounted and placed in the temperature controlled chamber at least 10 min before the measurement to ensure the whole sample adopts to the ambient temperature of the following experiment.Note that the time is increased according to the difference between room and ambient temperature.
In a first measurement step, the measuring cell is moved so that contact is established between the sample pins and the disk.The pins are moved until the resulting contact pressure reaches a preset target pressure.After that, the friction measurement starts during which the preset contact pressure from the previous step is maintained at a constant value.The measuring cell with the sample pins starts to rotate at a constant rotational speed of 0.005 rpm.The resulting sliding speed between the measurement cell and the disk is 0.007 mm/s.Preliminary measurements have shown that good and reproducible measurement results can be achieved at the selected rotational speed.The sampling rate of the measurement is 10 Hz.
The sample pins are made of the softer of the two contact partners' materials, which are the silicone rubbers Ecoflex 00-50 and Dragonskin 30, to prevent from measuring ploughing.In cases where a high COF occurs between the pin and disk, it must be ensured that the COF corresponds to the friction between the contact partners, and not between the pin and clamping.For this purpose, the material is glued on a 3D printed adapter.A commercially available silicone adhesive (Loctite 495) is used for this purpose.Care must be taken to ensure that the adhesive remains only in the adhesive zone and does not run into the contact surface, thus causing a change in the measurement.Measurements are conducted on clean samples.For the sample disks the following material samples were used: These materials are chosen to represent a broad range of different contact partners for SRs.In order to account for self-contact or contact between two SRs measurements are conducted where both samples, i.e. the pins and the disk, are made out of the silicones.The friction is very high in these cases which is a challenge for the measurement.In this case, special care must be taken to ensure that the sample pins remain mounted correctly through out the measurements.First, the measurements were performed at a reference temperature of 20 • C. Note that an overview of all measurements can be seen in Fig. 5.All measurements show the transition from static to kinetic COF over the sliding distance, for different influencing parameter.

Influence of contact pressure on coefficient of friction
First, the influence of the contact pressure on the COF was investigated at an ambient temperature of 20 • C. Therefore measurements of both pin materials with the above described contact partners were performed with different normal forces.The contact pressure can be calculated using the normal forces as follows: where F N,Pin = 1 3 F N is the normal force for each pin and A Pin = r 2 = 28.27mm 2 is the contact area.
The contact pressure is set based on Eq. 4 via the normal force of the measurement cell.The resulting contact pressures are 11.79kPa, 35.37kPa, and 58.95 kPa.In Fig. 5a, a measurement for Ecoflex and acrylic glass is shown as an example.The measurement was conducted at an ambient temperature of 20 • C. It can be seen, that increasing normal forces lead to decreasing COFs.Each measurement is repeated two times in order to ensure repeatability.The resulting COF can be seen in table 1.

Influence of the surface topology on the coefficient of friction
The surface topology can influence the frictional behavior.In order to showcase the dependency of the COF on the surface topology, two samples of the same material with different surfaces were used.For the first measurement, an aluminium disk cut out of a sheet is used.In Fig. 5e, measurements at the two different ambient temperatures are shown.For the second measurement the aluminium disk was lathed, showing the ring pattern from the manufacturing process.The surface roughness was not further investigated.The example is made to show the influence of the surface texture to create awareness that the surface topology itself can influence the frictional behavior.
The measurement results can be seen in Fig. 5d.The results show that the kinetic COF with a mean value of 0.56 is slightly smaller for the lathed aluminium sample disk which is 0.57 at an ambient temperature  C, the lathed aluminum's mean COF of 0.39 is significantly smaller than that of the smooth aluminum sheet, which is 0.5.This result demonstrates that surface topology has an influence on the frictional behavior and should be considered when using non-smooth surfaces.

Coefficient of friction for self-contact
To determine the COF for the self-contact case, the COF of Ecoflex and Dragonskin was measured.In Fig. 5f, the COF at an ambient temperature of 20 • C and for different normal forces are shown.The measurements show that the COF of the silicone rubbers is significantly higher than the COF with any other contact partner.For a normal force of 1 N, the average COF is 0.96.The measurements show a comparably high COF for silicone contact, which should be taken into account for simulations involving self-contact or contact of multiple SRs.

Friction map
In the previous sections, individual COF measurements were presented and discussed.We summarize our measurement results in a friction map in Table 2.For readability, we limit the table to a subset of entries and refer the reader to the open access friction map for the full set of measurements.The friction map may be used to find the appropriate COF for a specific simulation.Furthermore, the friction map can be used to interpolate a COF if the needed COF was not listed, providing the basis for a more reliable contact simulation.

Contact simulation
In order to demonstrate the difference between frictional and non-frictional contact and moreover the influence of varying COFs, we conducted a series of measurements in finite element simulation (ABAQUS Dassault Systemes).

Soft robot example
For this paper, we simulate the behavior of a soft pneumatic actuator, a common testbed for soft robotics research [39].The SPA has three pneumatic chambers radially distributed around the actuator's longitudinal axis as shown in Fig. 6.Each chamber is fiber-reinforced to prevent the actuator from ballooning.In this example, the main body is made of a soft silicone rubber named Ecoflex 00-50.At either end, the body is supported with two stiff caps to prevent the cross section from significantly deforming.Both caps are made of Dragonskin 30.For the following simulation, a second actuator was added to the model to account for a common application and to form a gripper.

Finite element simulation
For the simulation, both silicones are modeled with a hyperelastic Ogden material model.The supporting fibers are modeled using beam elements and elastic material properties.More details on the model and specifically the material models can be found in [39].
A cube (mass= 54 g) is modeled as the object to pick up as shown in Fig. 7. Contact properties are defined between the SPAs and the cube.For normal and tangential contact, a penalty formulation was used.For the simulation, both pneumatic actuators are fixed at the top end while the bottom end can move freely.The first static simulation step consists of fixing the cube on its top surface and placing it 4 mm away from each actuator as shown in Fig. 7.In an additional static step, the actuators' chambers are inflated using a pressure load so that the SPA bends towards the cube and gravity is applied in the z-direction as shown in Fig. 8.In a last dynamic implicit step, the boundary condition of the cube is changed to allow for displacement in the z-direction induced by gravity.In this step the cube is only held by the gripping actuators.The simulation time for this step is 1 sec.As demonstrated in Sect.3.2, the contact pressure and the resulting normal force influence the COF.To demonstrate the impact of this effect, the described simulations are performed with different chamber pressures which result in varying contact pressure between the actuator and the cube.In ABAQUS, COFs can be defined depending on the slip-rate, contact pressure, and ambient temperature.The friction map for the Dragonskin-steel contact pair (given in table 2) is implemented in ABAQUS for the tangential interaction property.A finite sliding formulation with a surface-to-surface contact discretization was chosen.The simulation is performed at an ambient temperature of 20 • C so the COF is implemented for a constant temperature.Values for the COF between the measured discrete values are linearly interpolated [18].In ABAQUS, a critical shear stress crit = p , which is defined by the COF and the contact pressure p, represents the maximal transmittable tangential force.It marks the transition before the state of sticking changes to gliding.If contact occurs, a constraint is applied to each node that is detected to be in contact.

Scenario: single chamber actuation (a) frictional and non-frictional contact, 80 kPa
In the first step of the analysis, the simulation is run with frictional and non-friction contact with an input pressure of 80 kPa in one chamber as shown in Fig. 7 on the left side (1).It is shown that the cube can only be held due to friction.In the case of non-frictional contact, the cube drops as expected (Fig. 8A).For frictional contact the cube is held by the SPAs as shown in Fig. 8B.
In the following simulation, input pressures of 50 kPa, 80 kPa, and 100 kPa are used to vary the contact pressure.The contact pressure as well as the displacement of the cube in the z-direction are evaluated.As expected, the contact pressure increases with increasing chamber pressure.
(b) frictional contact, 50 kPa For the lowest input pressure of 50 kPa, the cube drops because the COF is too low to reach a sufficient tangential contact force at this contact pressure.In contrast to non-frictional contact the cube drops more slowly after 0.15 sec compared to 0.08 sec.
(c) frictional contact, 80 kPa For an input pressure of 80 kPa, the cube moved 2.8 mm in the z-direction.This displacement indicates that the friction is still not sufficient to hold the cube but slows down the fall.
(d) frictional contact, 100 kPa For an input pressure of 100 kPa, a displacement of 1.5 mm is detected after a simulation time of 1 second.After 5 seconds, the displacement settles at a value of 1.7 mm when one chamber is actuated.In this case, the displacement is caused by a change in the cube's equilibrium point due to its weight.

Scenario: multi chamber actuation
In the next simulation, the actuators were rotated by 180 • and two chambers were each pressurized to 100 kPa for each actuator as shown in Fig. 7 on the right side (2).The actuation of two chambers leads to an additional increase in the contact pressure.
(a) frictional contact, 100 kPa For an input pressure of 100 kPa, a displacement of 1.2 mm is detected.It is shown that the maximum transmissible tangential force increases with increasing contact pressure, resulting in a lower displacement in the z-direction.

Conclusions
In this paper, the frictional behavior of commonly employed material in soft robotics is assessed experimentally.The coefficient of friction is provided in a friction map for common material combinations.Moreover, the resulting COFs were implemented in a finite element simulation to demonstrate how the choice of a COF influences the performance of a SRs with respect to a specified task.The simulation experiments demonstrate the importance of friction for SRs which come in contact with their environment.In this paper, we presented the influence of actuation pressure with respect to the frictional behavior of the SPA.For future work the input pressure can be optimally adjusted with respect to the object's weight, which can inform optimal control and actuation parameters which in turn can reduce actuation effort as well as material wear on the SPA.The presented friction map enables researchers to implement a pressure and temperature dependant COF for various materials.The COF is a function of velocity and also strongly dependent on the surface properties.Establishing and conducting experimental protocols to include these dependencies into the friction map and adding more measurements are subject to future work.The open source nature will enable sharing of results among the soft robotics research community.

3 . 3
Influence of the ambient temperature on the coefficient of friction Next, the influence of the ambient temperature on the COF was investigated.The measurement was conducted at two different ambient temperatures of 10 • C and 20 • C with the same normal force.In Fig.5b, a measurement for Dragonskin and Steel with a normal force F N = 3 N at ambient temperatures of 10 • C and 20 • C is shown as an example.The average COF for the lower ambient temperature of 10 • C is 0.41 and 0.45 at the higher ambient temperature of 20 • C.This difference was for example also observed in the case where a normal force of F N = 5 N was applied, as shown in Fig.5c.In this measurement, the average COF at the lower ambient temperature of 10 • C is 0.31 and 0.36 at the higher ambient temperature of 20 • C.

Fig. 5
Fig. 5 Overview of the conducted measurements: a Ecoflex sample pin on an acrylic glass disk to demonstrate the influence of the contact pressure at a reference temperature of 20 • C. b Dragonskin sample pin on a steel disk to demonstrate the influence of the ambient temperature with a normal force of 3 N. c Dragonskin sample pin on a steel disk to demonstrate the influence of the ambient temperature with a normal force of 5 N. d Dragonskin sample pin on a lathed disk subjected to a normal force of 3N and ambient temperatures of 10 • C and 20 • C. e Dragonskin sample pin on a smooth aluminum disk made from a sheet subjected to a normal force of 3 N and ambient temperatures of 10 • C and 20 • C. f Ecoflex sample pin on a Dragonskin disk to demonstrate the frictional behavior of self-contact with two silicone rubbers at an ambient temperature 20 • C for various normal forces

Fig. 6 Fig. 7
Fig.6 Schematic figure of the soft pneumatic actuator with three fiber-reinforced chambers according to[39]

Fig. 8
Fig.8 Example results for a simulation with non-frictional contact A and frictional contact B

Table 1
Dependency of the coefficient of friction on the contact pressure for Ecoflex and acrylic glass example at an ambient temperature of 20 • C • C. At a lower ambient temperature of 10

Table 2
Abstract of the friction map