modelContact

Information

<!DOCTYPE html><html>
<p>Depending on the shape, we use 1 (sphere), 2 (cylinder) or 4 (plane) points to describe the surfaces of the contact partners. These constitute potential contact points. For each of them the collision detection is performed. For this purpose, analytic solutions for simple geometries are provided in the library. As the contact region may alter with the moving bodies, the contact points will also move on the defined surface. </p>
<p>Then, the contact block calculates the appropriate force depending on the combination of surfaces. So, using it the respective combination of contact surface has to be chosen at first by setting the parameter contactDefinition. This will use the Modelica replaceable statement to define the appropiate components of the contact block. Then connect the contact interfaces of the two contact surfaces to the respective port of the contact block (first&nbsp;surface&nbsp;mentioned&nbsp;must&nbsp;be&nbsp;connected&nbsp;to&nbsp;port&nbsp;1).</p>
<p>In the case of a collision of the two connected surface (the contact condition holds for at least one contact point) a three-dimensional contact force is applied. It consists of both the normal force and the tangential friction. The respective directions can be obtained by means of the local coordinate systems in the contact points. As compared to more complex models, the continuous surface layer is replaced by a nonlinear spring/damper element. Consequently, the normal force Fn&nbsp;is determined by means of the penetration p&nbsp;and the penetration velocity. A continuous contact force model with hysteresis damping according to [1] is implemented. Nevertheless, choosing n1=1 and n=0&nbsp;one can get the linear Kelvin-Voigt model, where the coefficients are the spring and damping constant. Choosing n1=n2&nbsp;one will get the formulation according to [2].</p>
<p><img src="resources/IdealizedContact/Images/equations/Fn.jpg"/></p>
<p><br/>In order to calculate the friction forces without further discontinuous events, which would decrease the simulation speed and impede controller design, we use the continuously differentiable friction model of Makkar et al. [3]. They introduced the following function of the relative velocity&nbsp;to approximate the friction coefficient&nbsp;of the characteristic Stribeck curve.</p>
<p><img src="resources/IdealizedContact/Images/equations/mue.jpg"/></p>
<p>In doing so, no ideal static friction can be obtained because the actual force to be applied in the static state is independent from the relative velocity&nbsp;of the two bodies. Static friction is rather represented by sliding with very small relative velocities. To set the unknown constants gamma_i we use five parameters, which can be seen in the figure. The parameters mue_s&nbsp;and mue_k&nbsp;denote the coefficients of static and kinetic friction. The limit velocity v_e1&nbsp;and v_e2&nbsp;define the beginning of mixed and viscous friction. The latter is described by the proportionality factor k_v. The actual approximation can be monitored by calling the function <a href=" IdealizedContact.ContactBlock.plotFrictionCurve">plotFrictionCurve</a>.</p>
<p><br/><br/><img src="resources/IdealizedContact/Images/mue.jpg"/></p>
<p><br/>The complete vector of the contact force is then computed as follows.</p>
<p><img src="resources/IdealizedContact/Images/equations/Fcontact.jpg"/></p>
<p><b>Note:</b> The collision of two cylinders can lead to linear or punctiform contact regions. The calculation for these two cases is currently seperated in two blocks. Integration of the two blocks is in progress.</p>
<p><h4>References:</h4></p>
<p><br/>[1]&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; H. M. Lankarani, P. E. Nikravesh: Continuous Contact Force Models for Impact Analysis in Multibody Systems, Nonlinear Dynamics, 5, 1994 </p>
<p>[2]&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; K. H. Hunt, F. R. E. Crossley: Coefficient of restitution interpreted as damping in vibroimpact, ASME J. Appl. Mech, 1975 </p>
<p>[3]&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; C. Makkar, W. E. Dixon, W. G. Sawyer, G. Hu: A New Continuously Differentiable Friction Model for Control Systems Design, Proceedings of the 2005 IEEE/ASME International Conference on Advanced Intelligent Mechatronics, Monterey CA, July, 2005 </p>
</html>

Parameters

TypeNameDefaultDescription
Booleanexactfalse=true for exact contact point movement; =false for filtered movement
SI.Frequencyf10000filter frequency to filter contact point movement
Normal force
SI.TranslationalSpringConstantspringCoefficient1000000spring coefficient to calculate normal force
SI.TranslationalDampingConstantdampingCoefficient100000damper coefficient to calculate normal force
Realn11.5stiffness exponent
Realn2n1indentation exponent
SI.Distancep_max0.001maximum penetration depth
Stribeck curve
SI.CoefficientOfFrictionmue_k0.03coefficient of kinetic friction
SI.CoefficientOfFrictionmue_s0.04coefficient of static friction
SI.CoefficientOfFrictionmue_r0coefficient of rolling friction
Realk_v0gradient of viscous friction
SI.Velocityv_e10.01limit velocity of static friction
SI.Velocityv_e20.1limit velocity of kinetic friction
Friction model
Realgamma1(mue_s - tanh(2*v_e1/v_e2)*gamma4 - k_v*2*v_e1)/(tanh(2*v_e1/v_e1) - tanh(2*v_e1/v_e2)) - gamma4friction parameter 1
Realgamma22/v_e1friction parameter 2
Realgamma33/v_e2friction parameter 3
Realgamma4mue_k - k_v*v_e2friction parameter 4
Realgamma5gamma2friction parameter 5
Realgamma6k_vfriction parameter 6
Animation
Booleananimationtrue= true to animate contact points
Modelica.SIunits.RadiusradiusContactPoint0.0025radius of contact point animation
Modelica.Mechanics.MultiBody.Types.ColorcolorContactPoints1{0, 180, 0}color of contact points of body 1
Modelica.Mechanics.MultiBody.Types.ColorcolorContactPoints2{255, 0, 255}color of contact points of body 2

Connectors

TypeNameDefaultDescription
IdealizedContact.Interfaces.Contact_bPort1
IdealizedContact.Interfaces.Contact_bPort2

Components

TypeNameDefaultDescription
IdealizedContact.ContactBlock.PunctiformContact.SphereToRectanglecontactDefinition