User:Runsen.zhang/sandbox
In contact mechanics, a unilateral contact, also called unilateral constraint, denotes a mechanical constraint which prevents penetration between two rigid/flexible bodies. This kind of constraints is omnipresent in non-smooth multibody dynamics applications, such as granular flows[1], Legged robot, vehicle dynamics, particle damping, imperfect joints[2], landing rocket. In these applications, the unilateral constraints result in impacts happening, so it is important to give a proper method to deal with the unilateral constraints.
Modelling of the unilateral constraints
[edit]There are mainly two kinds of method to model the unilateral constraints. The first kind of method is inherited from continuum mechanics, and the second kind of method is based on the non-smooth theory.
Continuum mechanics based method
[edit]![](http://upload.wikimedia.org/wikipedia/commons/thumb/4/4e/Hertz_contact_animated.gif/170px-Hertz_contact_animated.gif)
In this method, normal forces generated from the unilateral constraints are obtained according to the local material properties of bodies. In detail, contact force models are derived from the continuum mechanics, and expressed as functions of gap and the impact velocity between bodies.The animation of classic Hertz contact law is shown the left figure, where the contact is explained by the local deformation of bodies. And more contact models can refer to the introduction in contact mechanics, and some reviews[3][4][5].
Non-smooth theory based method
[edit]In the non-smooth method, unilateral interactions between bodies are fundamentally modelled by the Signorini condition[6] for the non-penetration condition, and impact laws to define the impact process[7]. The Signorini condition can be concluded as the complementarity problem:
,
where means the distance between two bodies and denotes the contact force generated by the unilateral constraints, as shown in the below figure. Also, using the proximal point of convex theory, the Signorini condition can be equivalently expressed[6][8] as:
,
where denotes an auxiliary parameter, and represents the proximal point in the set to the variable [9], expressed as:
.
Both two expressions above represent the dynamic behaviours of unilateral constraints: when the normal distance is above zero, the contact is open, which means that there is no contact force between bodies, ; while when the normal distance is equal to zero, the contact is closed, resulting in .
![](http://upload.wikimedia.org/wikipedia/en/7/72/Contact_dynamics_unilateral.jpg)
During the implementation of non-smooth theory based method, velocity Signorini condition or the acceleration Signorini condition are actually employed in most cases. The velocity Signorini condition is expressed as[6][10]:
,
where means the relative normal velocity after impact. This velocity Signorini condition should be understood together with the Signorini . The acceleration Signorini condition is considered under the close contact (), as[8]:
,
where the over dot means the derivative with respect to time.
When using this method for unilateral constraints between two rigid bodies, the Signorini condition alone is not enough for the impact process, so impact laws, which give the information for states before the impact and after the impact, are needed to handle the impact process[6]. For example, when the Newton restitution law is employed, a coefficient of the restitution will be defined as: , where denotes the relative normal velocity before impact, to give information for the impact process.
Frictional unilateral constraints
[edit]For frictional unilateral constraints, the normal contact forces are described by the methods above, while the friction forces are commonly described by the Coulomb's friction law. The Coulomb's friction law can be concluded as: when the tangential velocity is not equal to zero, namely the sliding state, the friction force is proportional to the normal contact force ; when the tangential velocity is equal to zero, namely the stick state, the friction force is nor more than the maximum of the static friction force. This relationship can be summarised using the maximum dissipation principle[6], as
where
represents the friction cone, and denotes the kinematic friction coefficient. Similar as the normal contact force, using the proximal point of convex theory, the formulation above can be equivalently expressed as[6]:
,
where denotes an auxiliary parameter.
Solution techniques
[edit]If the unilateral constraints are modelled by the contact models, the contact forces generated from the unilateral constraints can be calculated directly from the explicit mathematical formula, which can refer to the reviews listed in Section Continuum mechanics based method. If the non-smooth theory based method is employed, there are mainly two kinds of formulations to solve the Signorini conditions: the nonlinear/linear complementarity problem (N/LCP) formulation and the augmented Lagrangian formulation. Compared to the solution preprocessing the contact models, the non-smooth method is more tedious, but requires less computation cost. More comparisons between solutions using contact models and non-smooth theory can refer to the paper[11].
N/LCP formulations
[edit]In this way, solution of dynamics equations with unilateral constraints is transformed into the solution of N/LCPs. In detail, for the frictionless unilateral constraints or the unilateral constraints with planar friction, the problem is transformed into LCPs, while for frictional unilateral constraints, the problem is transformed into NCPs. To solve LCPs, the pivoting algorithm, originating from the algorithm of Lemek and Dantzig, is the most popular method[8]. But unfortunately, numerical experiences show that the pivoting algorithm may fail when handling systems with a large number of unilateral contacts, even using the deepest optimizations[12]. For the NCPs, using a polyhedral approximation can transform the NCPs into a set of LCPs, and then be solved by the LCP solver [13] . Besides, other methods, like NCP-functions[14][15], cone complementarity problems (CCP) based method[16][17], are also employed to solve NCPs.
Augmented Lagrangian formulation
[edit]Different from the N/LCP formulations, the augmented Lagrangian formulation uses the proximal functions listed above, . Together with dynamics equations, this formulation is solved by root finding algorithms. A comparative study between LCP formulations and the augmented Lagrangian formulation can be found in the reference [18].
See also
[edit]- Multibody dynamics – Tool to study dynamic behavior of interconnected rigid or flexible bodies
- contact dynamics – Motion of multibody systems
- contact mechanics – Study of the deformation of solids that touch each other
- discrete element method – Numerical methods for computing the motion and effect of a large number of small particles
- Non-smooth mechanics – Modeling approach in mechanics
- Collision response
References
[edit]- ^ Anitescu, Mihai; Tasora, Alessandro (26 November 2008). "An iterative approach for cone complementarity problems for nonsmooth dynamics". Computational Optimization and Applications. 47 (2): 207–235. doi:10.1007/s10589-008-9223-4.
- ^ Flores, Paulo (7 March 2010). "A parametric study on the dynamic response of planar multibody systems with multiple clearance joints". Nonlinear Dynamics. 61 (4): 633–653. doi:10.1007/s11071-010-9676-8.
- ^ Machado, Margarida; Moreira, Pedro; Flores, Paulo; Lankarani, Hamid M. (July 2012). "Compliant contact force models in multibody dynamics: Evolution of the Hertz contact theory". Mechanism and Machine Theory. 53: 99–121. doi:10.1016/j.mechmachtheory.2012.02.010.
- ^ Gilardi, G.; Sharf, I. (October 2002). "Literature survey of contact dynamics modelling". Mechanism and Machine Theory. 37 (10): 1213–1239. doi:10.1016/S0094-114X(02)00045-9.
- ^ Alves, Janete; Peixinho, Nuno; da Silva, Miguel Tavares; Flores, Paulo; Lankarani, Hamid M. (March 2015). "A comparative study of the viscoelastic constitutive models for frictionless contact interfaces in solids". Mechanism and Machine Theory. 85: 172–188. doi:10.1016/j.mechmachtheory.2014.11.020.
- ^ a b c d e f Jean, M. (July 1999). "The non-smooth contact dynamics method". Computer Methods in Applied Mechanics and Engineering. 177 (3–4): 235–257. doi:10.1016/S0045-7825(98)00383-1.
- ^ Pfeiffer, Friedrich (14 March 2012). "On non-smooth multibody dynamics". Proceedings of the Institution of Mechanical Engineers, Part K: Journal of Multi-body Dynamics. 226 (2): 147–177. doi:10.1177/1464419312438487.
- ^ a b c Pfeiffer, Friedrich; Foerg, Martin; Ulbrich, Heinz (October 2006). "Numerical aspects of non-smooth multibody dynamics". Computer Methods in Applied Mechanics and Engineering. 195 (50–51): 6891–6908. doi:10.1016/j.cma.2005.08.012.
- ^ Jalali Mashayekhi, Mohammad; Kövecses, József (August 2017). "A comparative study between the augmented Lagrangian method and the complementarity approach for modeling the contact problem". Multibody System Dynamics. 40 (4): 327–345. doi:10.1007/s11044-016-9510-2. ISSN 1384-5640.
- ^ Tasora, A.; Anitescu, M. (January 2011). "A matrix-free cone complementarity approach for solving large-scale, nonsmooth, rigid body dynamics". Computer Methods in Applied Mechanics and Engineering. 200 (5–8): 439–453. doi:10.1016/j.cma.2010.06.030.
- ^ Pazouki, Arman; Kwarta, Michał; Williams, Kyle; Likos, William; Serban, Radu; Jayakumar, Paramsothy; Negrut, Dan (2017-10-13). "Compliant contact versus rigid contact: A comparison in the context of granular dynamics". Physical Review E. 96 (4): 042905. doi:10.1103/PhysRevE.96.042905. ISSN 2470-0045.
- ^ Anitescu, Mihai; Tasora, Alessandro (26 November 2008). "An iterative approach for cone complementarity problems for nonsmooth dynamics". Computational Optimization and Applications. 47 (2): 207–235. doi:10.1007/s10589-008-9223-4.
- ^ Xu, Ziyao; Wang, Qi; Wang, Qingyun (December 2017). "Numerical method for dynamics of multi-body systems with two-dimensional Coulomb dry friction and nonholonomic constraints". Applied Mathematics and Mechanics. 38 (12): 1733–1752. doi:10.1007/s10483-017-2285-8. ISSN 0253-4827.
- ^ Mangasarian, O. L. (July 1976). "Equivalence of the Complementarity Problem to a System of Nonlinear Equations". SIAM Journal on Applied Mathematics. 31 (1): 89–92. doi:10.1137/0131009. ISSN 0036-1399.
- ^ Fischer, A. (January 1992). "A special newton-type optimization method". Optimization. 24 (3–4): 269–284. doi:10.1080/02331939208843795. ISSN 0233-1934.
- ^ Melanz, Daniel; Fang, Luning; Jayakumar, Paramsothy; Negrut, Dan (June 2017). "A comparison of numerical methods for solving multibody dynamics problems with frictional contact modeled via differential variational inequalities". Computer Methods in Applied Mechanics and Engineering. 320: 668–693. doi:10.1016/j.cma.2017.03.010.
- ^ Negrut, Dan; Serban, Radu; Tasora, Alessandro (2018-01-01). "Posing Multibody Dynamics With Friction and Contact as a Differential Complementarity Problem". Journal of Computational and Nonlinear Dynamics. 13 (1): 014503. doi:10.1115/1.4037415. ISSN 1555-1415.
- ^ Jalali Mashayekhi, Mohammad; Kövecses, József (August 2017). "A comparative study between the augmented Lagrangian method and the complementarity approach for modeling the contact problem". Multibody System Dynamics. 40 (4): 327–345. doi:10.1007/s11044-016-9510-2. ISSN 1384-5640.
Further reading
[edit]Open-source software
[edit]Open-source codes and non-commercial packages using the non-smooth based method:
- Siconos – Open source scientific software for modeling non-smooth dynamical systems
- Chrono, an open source multi-physics simulation engine, see also project website
Books and articles
[edit]- Acary V., Brogliato B. Numerical Methods for Nonsmooth Dynamical Systems. Applications in Mechanics and Electronics. Springer Verlag, LNACM 35, Heidelberg, 2008.
- Brogliato B. Nonsmooth Mechanics. Communications and Control Engineering Series Springer-Verlag, London, 1999 (2dn Ed.)
- 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