# Impact Behavior of a Rotating Rigid Body with Impact and Viscous Friction

Authors: Dorian Cojocaru; Dan B. Marghitu
Journal: Mathematical Problems in Engineering
 (2020)

---

## Abstract

The impact between a rotating link and a solid flat surface is considered. For the impact, we consider three distinct periods: elastic period, elastoplastic period, and restitution period. A Hertzian contact force is considered for the elastic period. Nonlinear contact forces developed from finite element analysis are used for the remaining two phases. The tangential effect is taken into account considering a friction force that combines the Coulomb dry friction model and a viscous friction function of velocity. Simulations results are obtained for different friction parameters. An experimental setup was designed to measure the contact time during impact. The experimental and simulation results are compared for different lengths of the link.


---

## Body

## 1. Introduction

The study of the mechanical impact has a very long history starting from Maclaurin, Newton, and Poisson. They have produced impressive number of studies for different applications, and a diversity of practical benches were used to verify the theoretical approaches.An important improvement for this technical domain was produced through the development of new modern fields of engineering as computer aided design, robotics, and biomechanics. Modern sensors, high speed image acquisition, and processing devices and increased computer resources better support the ability for practical verification of the theoretical models. The diversity of applications must take into consideration the shape of impacting bodies, the materials, and the dynamics conditions including angles, speeds, accelerations and frictions.Different surfaces for the impact were considered: rod with a flat surface[1], sphere with a flat [2, 3], rotor system into a film damper [4], elastic wedges [5], sphere on a beam [6], anisotropic bisinusoidal surface and a rigid base [7], cylindrical body with a rigid plane [8], superball [9], tennis ball with a racket [10], and 3D-printed polymers under collision with a rigid rod [11]. Different types of behaviors depending of the nature of the constituent materials of the impacting bodies were studied: elastoplastic impact [2, 3, 6, 8, 12–14], nonlinear elastic [1, 15], elastic [16], elastic-perfectly plastic [17], and nearly complete elastic [14, 18]. The material properties were studied from the point of view of the influence of the impact behavior: incompressible isotropic [5], nanoparticles [15], 3D-printed polymers [11], transmembrane domains [19], stable CuO nanoparticle enhanced lubricants [16], Al-6061 instrumented spherical microindentation and microstructurally graded samples [20], and various finishes [21]. The direction of the impact is another factor to be considered in modeling the impact: normal [6], oblique [10, 13, 22], normal, and tangential [7].Friction could have an important effect in many practical applications and by consequence the theoretical model must be adapted. Studies have considered the effect of friction in combination with other influences such as strain hardening [12], wear [15], lubrication [16], and deformation [22]. Better solutions for applications involving impacts with friction were proposed as the treatment of impact for the case of planar multibody mechanical systems [23] and the use of an efficient contact detection algorithm [24].A friction model was proposed in order to make possible a correct detection of slipping, stop sliding, and reverse slipping [25]. Both cases of open and closed loop are analyzed. In order to avoid the use of Newton’s hypothesis, Poisson’s kinetic coefficient of restitution is proposed. Cartesian coordinates and Cartesian momenta are utilized. The final velocities are calculated momentum equations for impact. The theory is verified on some practical applications: the impact of a two-link pendulum and the impact of automotive back wheel.The dynamics of a class of multibody mechanical systems, involving contact and friction, was also analyzed [26]. The authors assume as main contribution the development of a new method for detecting contact among the components of a system of solid bodies. The geometry of impacting system is very important from the method used to solve the problems. Analytical and numerical calculations are used for simple and complex geometries. The needed information is related to the force arising between the contacting bodies. The new method was tested by practical applications, and a corresponding methodology was proposed. A comparison between the sliding and rolling friction was performed by Alaci et al. [27]. Following the Amonton–Coulomb law, the proportionality between the friction force and the normal force was assumed for the case of the sliding friction. The dependency between the rolling friction and the normal force has been demonstrated to be linear only for specific domains. A technique has been determined in order to obtain the critical value of the coefficient of rolling friction.A specialized software has been used for the friction behavior in Ghaednia and Jackson [15]. Specific tribology application has been targeted, but the results were applicable for friction modeling under the conditions of industrial sheet metal forming processes. Considering the influence of the temperature, Nevmerzhitskiy et al. [28] developed a friction model of industrial robot joint.For the stance phase of the human gait, estimation for the friction force, considering both dry and viscous components, was studied by Oje and Mayo [29]. The authors stated that it is very important to have a correct estimation for the contact parameters, both normal and tangential. The researchers still look for the best estimation of tangential forces. They did not use the components of the tangential force but the modulus of this force.Zobova [30] considered the dry friction for an element in motion on a plane with friction. In this application, the plane suffers a deformation during the motion. In order to obtain the resulting Coulomb friction force, an integration process was developed for the infinitesimal components. For this case, the estimation of normal force reactions must take into consideration the position, the orientation, the mass center velocity, and the rotation of the element.The problem of the elastic impact is considered by Thornton et al. [31]. The application deals with spheres obliquely impacting a plane, a wall. The authors propose four different contact force models and compare the obtained results. The goal is to calibrate the normal spring stiffness by using a tangential to normal spring stiffness. The conclusion was that a linear spring can support very good results, and a nonlinear spring asks for incremental computation. A similar approach has been developed by Palli et al. [32].A model for a single-link flexible robotic arm is studied by Chapnik et al. [33]. The model took into consideration the effects of beam damping and hub inertia. The Coulomb and viscous hub friction are computed. The model should be used also in applications with elastic impacts.An application where face seals are designed to be in contact at standstill is presented in Green [34]. The viscous and dry frictions heating are the cause of the deformation of the seal faces. In order to obtain both, a stable operation and a favorable dynamic behavior, it is needed that the upper deformations form a converging gap for radial flow. The proposed solution is based on an elastoplastic contact model. A viscous heating model was used to compute the time-dependent deformation. A plastic deformation behavior was studied by Fu [35].Verifying through simulation of the abovementioned conclusions, we are proposing a trustful mathematic support for the design of mechatronics and robotic systems subjected to collisions during service. Taking into consideration the effect of friction should contribute to increasing the precision and safety for the corresponding applications. The mechanical impact study is important for robotic applications, including robotic systems for health applications. In this case one of the impacting bodies is a robotic link. The authors of this paper acquired an important experience studying different such cases as kinematic of rebound [36, 37], the impact between a rotating link and a massive surface [38], the impact between a planar link and a granular material [39], simultaneous impact of two links chain [19], and the variation of the coefficients of restitutions for different type of impact [26, 40, 41].Using the computer controlled wheelchair people with special needs can be supported also by a robotic arm picking up objects, turning on or off switches, opening or closing doors, making tea, and pulling drawers [42]. Taking into account the particularities of such applications, it is necessary to take into account the mechanical impact when designing robotic structures. In this way an adapted mechanical structure will be obtained, but it will also be possible to design and implement a more efficient control system. In order for the evaluation of the real situation to be as efficient as possible, the consideration of friction is an important factor.The results presented in this paper were not intended for a particular robotic application but are useful to support the design process in several types of robotic applications. Some of these applications will be presented below.The obtained results are useful for supporting the design of robotic arms that involve coupling with a pendulum [43].The study of balance behavior is interesting for several types of mobile robots. The intrinsic ability of a robotic structure, including those that use couplings such as those studied in this paper, to function in a controllable manner, both in the absence of collisions and in the presence directly influences the overall performance of the system. A significant application for this aspect is the simplest case of a self-balancing mobile robot. In this situation, 2D balancing is desired for a flat robot with a double pendulum that makes a single point contact with a support surface [44].By copying the abilities of some living creatures (e.g., a dog [45]), robotic structures have been designed that can make jumps to overcome obstacles. To make a jump, the robot must be designed in such a way that it can go through three successive actions: take-off, evolution in the air, and landing. The corresponding robotic structures are based on mechanisms that include structures such as the one studied in the paper for the implementation of hip, knee, and ankle joints.Robots with navigational skills on terrain with complex geometry can be designed modularly. The performance of the modular design can be improved by increasing the number of modules, including modules such as the one presented in the paper [46]. Sustained research efforts are being made for the use of robotic structures in medical applications where knee prostheses, hand prostheses, and exoskeletons are needed [47]. These robotic structures include couplings of the type studied in this paper. A representative example is when people with special needs who use wheelchairs also benefit from exoskeletons for walking assistance. In this way, these people can sit, walk, and climb stairs.Most of the applications presented above assume by their nature the existence of collisions and the influence of friction. The robotic structures can be controlled by people with special needs by manipulating a joystick but also by eye gaze [48], by gesture [49], by head motion [50], or by signals made with the fingers. In these applications it is absolutely necessary to seriously treat collisions between controlled robotic structures and obstacles.The main objective of this paper is to propose a model for the mechanical impact, meant to support a more accurate design procedure for mechatronics and robotics systems. We present a mathematical model taking into consideration the viscous friction during impact. The relations between the impact and the friction are presented. Simulation and experimental results are developed.


## 2. Mathematical Model

For this research, in order to study the impact between a compound pendulum and a flat surface, we use a Coulomb model for the case of viscous friction. To the best of our knowledge, this frictional model has not been employed for the impact problems.In Figure1, a compound pendulum during the impact is represented, where ϕ is the angle of impact. It is needed to define a three-axis inertial reference frame. The compound link is moving in a plane, and we designated it as the xy-plane. The upper mentioned plane has two axes, the vertical one designated as y-axis (having the vertically direction downward as the positive sense), and the horizontal one designed as x-axis (having the sense from left to right as a positive sense). The unit vector of x-axis is ı and the unit vector of y-axis is j. The plane of motion has a normal axis, z, with k as unit vector.Figure 1
Rotating robotic link with impact.Having a uniform link, the mass centerC has the same coordinates as the geometric center. There is a point fixed to the ground denominated the pivot point O. There is another point T being the end of the link. This point gets in contact with a fixed plane surface. Because the length of the rotating link is L, the distance between the contact point T and the pivot point O is L, and the distance between the mass center C and the pivot point O is L/2. In consequence, the vector from O to C is rOC=r=xC1+yCj, where xC is the horizontal coordinate of the mass center and yC is the vertical coordinate of the mass center xC=L/2cosϕ,yC=L/2sinϕ. In consequence, the position vector of T isjrT=xT1+yTj, where xT and yT are the coordinates of TxT=Lcosϕ,yT=Lsinϕ. Hence,(1)rT=Lcosϕı+Lsinϕj.Taking into account that, as it was already stated, the pendulum is moving in a plane, making a pure rotation around the pivot point, it follows that the angular velocity and angular acceleration vectors follow thez-direction being perpendicular to this plane. We have only one moving rigid body and for the angular velocity referring to this body, the rod, we will use the following expression ω=dϕ/dt=ϕ˙, where ω is the angular velocity of the rotating link. Further, making the step from the angular velocity to angular acceleration of the kinematic element, we will use the following expression: α=d2ϕ/dt2=ϕ¨, where the angular acceleration of the rotating link is α. It is needed to compute the mass moment of inertia of the rod rotating around the fixed pivot point O. Considering the second moment of inertia around the center of mass, one way to solve this problem is to use the transfer theorem. Thus, IO=IC+mL/22. The mass of the rotating link is m. There is a moment supporting the rotation of the rod around the pivot point. This moment is determined by the weight of the link acting through its mass center and could be computed by the vectorial product of the position vector r crossed into G=mgj, where g is the gravitational acceleration.At collision instant an impact normal force,Fn with the point of application at T is acting upward on y−axis. The elastic deformation due to the normal impact force is given by(2)δ=Lsinϕ−Lsinϕ0,where ϕ0 is the initial impact angle.The sum of the moment with respect toO is sum of the vector product of the position vector r crossed into G and the vector product of the position vector rT crossed into the impact force:(3)IOαk=∑MO=r×G+rT×Fn+Ff,where the friction force during impact is Ff. The impact interval is divided into compression phase and restitution phase. The compression phase is initiated at the contact of the impact point with the solid surface and continues until the normal velocity of T is zero. The compression phase is divided into an elastic interval and next an elastoplastic interval. The restitution interval starts at the maximum compression and ends until the impact force is zero.


## 3. Impact Force

For the normal impact force, a formulation developed by Jackson and Green [51, 52] is employed.The Hertz theory is considered for the elastic phase during the impact and the restitution phase. The formula for the reduced modulus of elasticity is(4)E−1=1−ν12E1+1−ν22E2,and the reduced radius is defined as(5)1R=1R1+1R2,where E1 and E2 are the modulus of elasticity of the rod and the flat surface, R1 and R2=∞ are the radii of curvature of the rod and the flat, and ν1 and ν2 are the Poisson ratios of the rod and the flat surface. The impact force for the elastic phase is given by(6)Fn=Fe=43ER0.5δ1.5.The elastic phase is ending when(7)δ∗=δδy>1.9,where(8)δy=RπCjSy2E2andCj=1.295e0.736ν.The elastic displacement at which the yield starts isδy and the yield strength of the weaker material is Sy.For the elastoplastic phase, the force is calculated with a modified Jackson and Green expression [13, 14].(9)Fn=Fep=Pce−0.17δ/δy5/12δδy1.5+4HCjSy1−e−1/78δ/δy5/9δδy1.1,where(10)B=0.14e23Sy/E,a=Rδyδ1.9δyB,HSy=2.84−0.921−cosπaR,Pc=43RE2πCjSy23.The average normal pressure isH.For the restitution, an elastic impact force is considered.(11)Fn=Fr=43ER0.5δ−δr1.5,where the permanent deformation is calculated from a modified curve-fitting from the finite element models [13, 22].(12)δr=δm0.81−δm+5.56.5−2.The maximum deformation at the end of the compression phase isδm. The radius of curvature for the restitution phase is(13)Rr=1δm−δr33Pm4E2,where Pm=Fepδ=δm is the maximum elastoplastic force.The friction is given by a combination of the dry Coulomb force and a viscous force:(14)Ff=−μFnvT⋅11vT⋅1−kvvT⋅1,where μ is the kinematic coefficient of friction and kv is the viscous coefficient of friction. The velocity of the impact point is defined as(15)vT=ωk×rT.


## 4. Simulation Results

As shown above, the purpose of the paper was to provide a solution that can support the process of designing robotic structures. These structures have applicability in several fields, which have also been presented above. Under these conditions, no specific customization was made for a robotic structure. One could apply the method for designing a small stepping robot that impacts a metal surface, and the impact also involves friction.During the simulation, many experiments were performed with different conditions, giving multiple values to the quantities involved in the mathematical model. The paper included, by way of example, a particular experiment which was considered representative, but no particular solution was sought for a specific application.To find the equations of motion, symbolical MATLAB was used. The numerical solution for the equations of motion was obtained using an ODE45 function.The air resistance has a small value and was neglected during impact. For this application with impacts with rebounds at low velocity, the air resistance can be neglected even before and after impact.We considered only the rigid body motion of the link. The elastic deformations were taken into account only at contact point. For future studies, we intend to include the longitudinal and transverse vibrations of the link.Taking into account the assumed objectives and the type of robotic coupling chosen, the parameters to be calculated by simulation for the different experiments were identified.Input variables that differentiate the experiments were chosen: geometric dimensions of the bars (length and diameter), properties of the material from which the bars are made (density, Young’s modulus, Poisson ratio, and Yield strength), and nature of friction (coefficient of friction and coefficient of viscous friction).Based on the mathematical model presented above and taking into account the behavior of the robotic couple, the kinematic and dynamic quantities for a design process were calculated for the period and point of the impact: the impact angle, angular velocity, elastic displacement, velocity (normal and tangential), kinetic energy variation, normal force, velocities of the impact point function of the coefficient of friction (tangential, normal, and angular), and rebound velocity of the impact point as a function of the viscous coefficient of friction (tangential, normal, and angular).Considering the above, the following describes a particular experiment. The values obtained by simulation are commented in terms of variation and dependencies.The motion of a homogeneous link, as shown in Figure1, with the length L=0.25m and the diameter d=2R=0.01m is simulated. The link and the flat are made of steel with the same characteristics. The density of the link is ρ=7800kg/m3 and Young’s modulus is E=210109Pa. The Poisson ratio is ν=0.29 and the Yield strength is SY=1.12109Pa. For the first simulation, the coefficient of friction is μ=0.3 and the coefficient of viscous friction is kv=80. In Figure 2, the impact angle ϕ during the impact is shown. The elastic compression is very small. The elastoplastic phase is dominant and the impact angle is increasing until the maximum compression. The last phase is the restitution phase. The classical stereomechanics considered constant position during impact. The rotation of the link for the contact period is depicted in Figure 3. The angular velocity changes the sign at the maximum compression. The initial impact angular velocity is ωti=7.67rad/s and the final rebound angular velocity is ωtf=−3.15rad/s. At the maximum compression, the angular velocity is zero.Figure 2
Impact angle,ϕ, during impact.Figure 3
Angular velocity,dϕ/dt, during impact.The elastic displacementδ is depicted in Figure 4. The maximum elastic displacement is at maximum compression and the permanent deformation is δr=4.1710−5m. The y-component of the velocity of T for the impact period is represented in Figure 5. The initial perpendicular velocity is vTyti=1.66m/s. The final normal rebound velocity is vTytf=−0.68m/s. The x-component of velocity of the impact point T during the collision is depicted in Figure 6. The rotating element contacts the solid surface and starts sliding. The initial tangential velocity is vTxti=−0.96m/s. The final tangential rebound velocity is vTxtf=0.39m/s. The tangential velocity changes sign at maximum compression.Figure 4
Elastic displacement,δ, during impact.Figure 5
Normal velocity of the pointT during impact.Figure 6
Tangential velocity of the pointT during impact.The variation of the kinetic energy during impact is depicted in Figure7. The normal force during impact is shown in Figure 8. For the next simulations, the viscous coefficient of friction is neglected. Next, the tangential velocities of the impact point as a function of the coefficient of friction μ are represented in Figure 9. The tangential rebound velocity is decreasing with the increasing of the coefficient of friction.Figure 7
Kinetic energy variation during impact.Figure 8
Normal force during impact.Figure 9
Tangential velocities of the impact point function of the coefficient of frictionμ.The normal velocities of the impact point as a function of the coefficient of frictionμ are represented in Figure 10. The magnitude of the normal rebound velocity is decreasing with the increase of the coefficient of friction. If the coefficient of friction is increasing, the energy loss is increasing and the final magnitude of the rebound angular velocity is decreasing. The decreasing of the magnitude of the angular velocity during impact with the increasing of the coefficient of friction μ is shown in Figure 11.Figure 10
Normal velocities of the impact point function of the coefficient of frictionμ.Figure 11
Angular velocity of the link as a function of the coefficient of frictionμ.For the next simulations, the coefficient of friction wasμ=0.1. In Figure 12 the rebound tangential velocity function of the coefficient of viscous friction kv is represented. The final tangential velocity is decreasing with respect to the viscous coefficient of friction. The magnitude of the normal rebound velocity is decreasing with respect to the viscous coefficient of friction as shown in Figure 13. The magnitude of the rebound angular velocity of the link is decreasing with respect to the viscous coefficient of friction as depicted in Figure 14.Figure 12
Tangential rebound velocity of the impact point as a function of the viscous coefficient of friction.Figure 13
Normal rebound velocity of the impact point as a function of the viscous coefficient of friction.Figure 14
Angular rebound velocity of the link as a function of the viscous coefficient of friction.


## 5. Experimental Setup

The mathematical model was obtained by combining in a manner adapted to the experiment in question. To validate this model, two stages were completed: numerical simulation and experimental validation.Using the numerical simulation, the mathematical model was refined, and, using the latest version of this model, the validation was done by experimental testing. In the experimental stage, bars made of different materials and having different geometric dimensions were used.Experimental data are acquired and stored in digital format. This data is then processed offline. As a result, processing time is not a problematic condition. Numerical simulation also does not pose problems in terms of processing cost, neither in time, nor in memory or in computing capacity.Figure15(a) shows the main scheme of the practical experiment. The whole structure is metallic ensuring electrical conductivity. When the coupling makes mechanical contact with the horizontal base surface, an electrical contact is also made. The resulting electrical signal is variable for a period of time due to the nature of the mechanical impact with friction. A digital oscilloscope records the signal (Figure 15(b)) which is then after that analyzed offline. To increase the accuracy of the procedure, a robotic arm with a gripper will release the coupling which, by falling, will impact the horizontal surface. A digital video system will take pictures during the fall of the couple. In this way, it is possible to determine with precision the positions in the 2D plane of evolution of the couple.Figure 15
(a) Experimental setup; (b) example of an image recorded by the digital oscilloscope.
(a)(b)The following is the information that defines one of the sets of experiments that have been performed. For the impacting link, four dimensions are used:L1=0.10m, L2=0.15m, L3=0.18m, and L4=0.24m. The diameter of the bars is d=0.008m. The initial position of the rotating link is defined by the angle ϕ=70°. The bar is released from rest in a horizontal position. For the simulation, the following coefficients are used: μ=0.3 and kv=100.In Figure16, the mean experimental data and the simulation results are shown. The experimental and the simulation results are in good agreement. The impact time is increasing with the length of the link. For longer links, the permanent deformation is larger and the impact time will be longer.Figure 16
Time during impact for different lengths,L, of the rotating link.


## 6. Conclusions

The collision with viscous friction of a robotic link with a solid surface was developed. The three phases during collision are elastic compression, elastic plastic compression, and restitution. The rebound velocities were calculated function of the coefficient of friction and the viscous coefficient of friction. The rebound angular velocity of the link is increasing with respect to the viscous coefficient of friction. The rebound tangential velocity is decreasing with respect to the viscous coefficient of friction. To validate the theoretical model, experimental results are obtained. The duration of the impact for experimental and simulation values are similar. Verifying through simulation and experiments of the abovementioned conclusions, we are proposing a trustful mathematic support for the design of mechatronics and robotic systems subjected to collisions during service, taking into consideration that the effect of friction should contribute to increasing the precision and safety for the corresponding applications.



---

*Source: 5471629-2020-11-03.xml*