Admin: 1 revision imported: Mitul Imported from Wikipedia

2021-04-11T19:43:22Z

1 revision imported: Mitul Imported from Wikipedia

setup>SimoneD89: changed category to the correct wikiproject

2021-04-05T09:31:55Z

changed category to the correct wikiproject

New page

{{short description|Motion of multibody systems}}
{{expert needed|1=Systems|date=June 2018}}
'''Contact dynamics''' deals with the motion of [[multibody system]]s subjected to [[unilateral contact]]s and [[friction]]. Such systems are omnipresent in many multibody dynamics applications. Consider for example
* Contacts between wheels and ground in [[vehicle dynamics]]
* Squealing of brakes due to friction induced oscillations
* Motion of many particles, spheres which fall in a funnel, mixing processes (granular media)
* [[Clockwork]]s
* Walking machines
* Arbitrary machines with limit stops, friction.
*Anatomic tissues (skin, iris/lens, eyelids/anterior ocular surface, joint cartilages, vascular endothelium/blood cells, muscles/tendons, et cetera)

In the following it is discussed how such mechanical systems with unilateral contacts and friction can be modeled and how the time evolution of such systems can be obtained by [[Numerical ordinary differential equations|numerical integration]]. In addition, some examples are given.

==Modeling==
The two main approaches for modeling mechanical systems with unilateral contacts and friction are the regularized and the non-smooth approach. In the following, the two approaches are introduced using a simple example. Consider a block which can slide or stick on a table (see figure 1a). The motion of the block is described by the equation of motion, whereas the friction force is unknown (see figure 1b). In order to obtain the friction force, a separate force law must be specified which links the friction force to the associated velocity of the block.

[[Image:contact dynamics block.jpg|frame|center|Figure 1: Block which can slide or stick on a table. Figure a) depicts the model, figure b) the equation of motion with unknown friction force]]

===Non-smooth approach===

A more sophisticated approach is the non-smooth approach, which uses set-valued force laws to model mechanical systems with unilateral contacts and friction. Consider again the block which slides or sticks on the table. The associated set-valued friction law of type Sgn is depicted in figure 3. Regarding the sliding case, the friction force is given. Regarding the sticking case, the friction force is set-valued and determined according to an additional algebraic [[Constraint (mathematics)|constraint]].

[[Image:contact dynamics setvalued.jpg|frame|center|Figure 3: Set-valued force law for friction]]

To conclude, the non-smooth approach changes the underlying mathematical structure if required and leads to a proper description of mechanical systems with unilateral contacts and friction. As a consequence of the changing mathematical structure, [[Impact force|impacts]] can occur, and the time evolutions of the positions and the velocities can not be assumed to be [[Smooth function|smooth]] anymore. As a consequence, additional impact equations and impact laws have to be defined. In order to handle the changing mathematical structure, the set-valued force laws are commonly written as [[inequality (mathematics)|inequality]] or [[Inclusion (set theory)|inclusion]] problems. The evaluation of these inequalities/inclusions is commonly done by solving linear (or nonlinear) [[Linear complementarity problem|complementarity problems]], by [[quadratic programming]] or by transforming the inequality/inclusion problems into projective equations which can be solved iteratively by [[Jacobi method|Jacobi]] or [[Gauss–Seidel method|Gauss–Seidel]] techniques.
The non-smooth approach provides a new modeling approach for mechanical systems with unilateral contacts and friction, which incorporates also the whole classical mechanics subjected to bilateral constraints. The approach is associated to the classical [[Differential algebraic equation|DAE]] theory and leads to robust integration schemes.

==Numerical integration==

The integration of regularized models can be done by standard stiff solvers for ordinary differential equations. However, oscillations induced by the regularization can occur. Considering non-smooth models of mechanical systems with unilateral contacts and friction, two main classes of integrators exist, the event-driven and the so-called time-stepping integrators.

===Event-driven integrators===

Event-driven integrators distinguish between smooth parts of the motion in which the underlying structure of the differential equations does not change, and in events or so-called switching points at which this structure changes, i.e. time instants at which a unilateral contact closes or a stick slip transition occurs. At these switching points, the set-valued force (and additional impact) laws are evaluated in order to obtain a new underlying mathematical structure on which the integration can be continued. Event-driven integrators are very accurate but are not suitable for systems with many contacts.

===Time-stepping integrators===

So-called time-stepping integrators are dedicated numerical schemes for mechanical systems with many contacts. The first time-stepping integrator was introduced by J.J. Moreau. The integrators do not aim at resolving switching points and are therefore very robust in application. As the integrators do work with the integral of the contact forces and not with the forces itself, the methods can handle both non-impulsive motion and impulsive events like impacts. As a drawback, the accuracy of time-stepping integrators is low. This lack can be fixed by using a step-size refinement at switching points. Smooth parts of the motion are processed by larger step sizes, and higher order integration methods can be used to increase the integration order.

==Examples==

This section gives some examples of mechanical systems with unilateral contacts and friction. The results have been obtained by a non-smooth approach using time-stepping integrators.

===Granular materials===
Time-stepping methods are especially well suited for the simulation of granular materials. Figure 4 depicts the simulation of mixing 1000 disks.

[[Image:contact dynamics mixer.jpg|frame|center|Figure 4: Mixing a thousand disks]]

===Billiard===
Consider two colliding spheres in a billiard play. Figure 5a shows some snapshots of two colliding spheres, figure 5b depicts the associated trajectories.

[[Image:contact dynamics billiard.jpg|frame|center|Figure 5: a) Snapshot. b) Trajectories of the two spheres]]

===Wheely of a motorbike===
If a motorbike is accelerated too fast, it does a wheelie. Figure 6 shows some snapshots of a simulation.

[[Image:contact dynamics motorbike.jpg|frame|center|Figure 6: Wheely of a motorbike]]

===Motion of the woodpecker toy===

The woodpecker toy is a well known benchmark problem in contact dynamics. The toy consists of a pole, a sleeve with a hole that is slightly larger than the diameter of the pole, a spring and the woodpecker body. In operation, the woodpecker moves down the pole performing some kind of pitching motion, which is controlled by the sleeve. Figure 7 shows some snapshots of a simulation.

[[Image:contact dynamics woodpecker.jpg|frame|center|Figure 7: Simulation of the woodpecker toy]]

A simulation and visualization can be found at https://github.com/gabyx/Woodpecker.

== See also ==
*[[Multibody system|Multibody dynamics]]
*[[Contact mechanics]]: Applications with unilateral contacts and friction. Static applications (contact between deformable bodies) and dynamic applications (Contact dynamics).
*[[Lubachevsky-Stillinger algorithm]] of simulating compression of large assemblies of hard particles

==References==
{{Unreferenced|date=May 2008}}

== Further reading ==
*Acary V. and Brogliato, B. Numerical Methods for Nonsmooth Dynamical Systems. Applications in Mechanics and Electronics. ''Springer Verlag, LNACM 35, Heidelberg, 2008.''
*Brogliato B. Nonsmooth Mechanics. Models, Dynamics and Control ''Communications and Control Engineering Series'' Springer-Verlag, London, 2016 (third Ed.)
* Drumwright, E. and Shell, D. Modeling Contact Friction and Joint Friction in Dynamic Robotic Simulation Using the Principle of Maximum Dissipation. ''Springer Tracks in Advanced Robotics: Algorithmic Foundations of Robotics IX'', 2010
*Glocker, Ch. ''Dynamik von Starrkoerpersystemen mit Reibung und Stoessen'', volume 18/182 of ''VDI Fortschrittsberichte Mechanik/Bruchmechanik.'' VDI Verlag, Düsseldorf, 1995
*Glocker Ch. and Studer C. Formulation and preparation for Numerical Evaluation of Linear Complementarity Systems. ''Multibody System Dynamics'' 13(4):447-463, 2005
*Jean M. The non-smooth contact dynamics method. ''Computer Methods in Applied mechanics and Engineering'' 177(3-4):235-257, 1999
*Moreau J.J. '' Unilateral Contact and Dry Friction in Finite Freedom Dynamics,'' volume 302 of '' Non-smooth Mechanics and Applications, CISM Courses and Lectures''. Springer, Wien, 1988
*Pfeiffer F., Foerg M. and Ulbrich H. Numerical aspects of non-smooth multibody dynamics. ''Comput. Methods Appl. Mech. Engrg'' 195(50-51):6891-6908, 2006
*Potra F.A., Anitescu M., Gavrea B. and Trinkle J. A linearly implicit trapezoidal method for integrating stiff multibody dynamics with contacts, joints and friction. ''Int. J. Numer. Meth. Engng'' 66(7):1079-1124, 2006
*Stewart D.E. and Trinkle J.C. An Implicit Time-Stepping Scheme for Rigid Body Dynamics with Inelastic Collisions and Coulomb Friction. ''Int. J. Numer. Methods Engineering'' 39(15):2673-2691, 1996
*Studer C. ''Augmented time-stepping integration of non-smooth dynamical systems'', PhD Thesis ETH Zurich, ETH E-Collection, to appear 2008
*Studer C. ''Numerics of Unilateral Contacts and Friction—Modeling and Numerical Time Integration in Non-Smooth Dynamics'', Lecture Notes in Applied and Computational Mechanics, Volume 47, Springer, Berlin, Heidelberg, 2009

== External links==
* [http://www.zfm.ethz.ch/e/res/mul/index.htm Multibody research group], Center of Mechanics, ETH Zurich.
* [http://www.amm.mw.tum.de/index.php?id=34&L=1 Lehrstuhl für angewandte Mechanik] TU Munich.
* [http://www.inrialpes.fr/bipop/index.html BiPoP Team], INRIA Rhone-Alpes, France,
* [[Siconos]] software. An open-source software dedicated to the modeling and the simulation or nonsmooth dynamical systems, especially mechanical systems with contact and Coulomb's friction
* [http://www.cs.rpi.edu/~trink/rigid_body_dynamics.html Multibody dynamics], Rensselaer Polytechnic Institute.
* [http://www.zfm.ethz.ch/dynamY/index.html dynamY software]
* [http://www.lmgc.univ-montp2.fr/~dubois/LMGC90 LMGC90 software]
* [https://www.migflow.be/ MigFlow software]
* [http://code.google.com/p/solfec/ Solfec software]
* [https://gabyx.github.io/GRSFramework/ GRSFramework] Granular Rigid Body Simulation Framework developed at [http://www.imes.ethz.ch/ IMES] in Ch. Glocker's group (High-Performance Computing with MPI), 2016
* [https://github.com/projectchrono/chrono/ Chrono], an open source multi-physics simulation engine, see also project [http://projectchrono.org website] 2017

[[Category:Mechanics]]
[[Category:Dynamical systems]]

← Older revision	Revision as of 19:43, 11 April 2021
(No difference)

Contact dynamics - Revision history

Admin: 1 revision imported: Mitul Imported from Wikipedia

setup>SimoneD89: changed category to the correct wikiproject