A mortar based contact formulation for non-linear dynamics using dual Lagrange multipliers

Many existing algorithms for the analysis of large deformation contact problems use a so-called node-to-segment approach to discretize the contact interface between dissimilar meshes. It is well known, that this discretization strategy may lead to problems like loss of convergence or jumps in the contact forces. Additionally it is popular to use penalty methods to satisfy the contact constraints. This necessitates a user defined penalty parameter the choice of which is somehow arbitrary, problem dependent and might influence the accuracy of the analysis. In this work, a frictionless segment-to-segment contact formulation is presented that does not require any user defined parameter to handle the non-linearity of the contact conditions. The approach is based on the mortar method enforcing the compatibility condition along the contact interface in a weak integral sense. The application of dual spaces for the interpolation of the Lagrange multiplier allows for a nodal decoupling of the contact constraints. A local basis transformation in combination with a primal-dual active set strategy enables the exact enforcement of the contact constraints via prescribed incremental boundary conditions. Due to the biorthogonality condition of the basis functions the Lagrange multipliers can be locally eliminated. A static condensation leads to a reduced system of equations to be solved solely for the unknown nodal displacements. Thus the size of the system of equations remains constant during the whole calculation. The discrete Lagrange multipliers, representing the contact pressure, can be easily recovered from the displacements in a variationally consistent way. For the analysis of dynamic contact problems the proposed contact description is combined with the implicit Generalized Energy-Momentum Method. Several numerical examples illustrate the performance of the suggested contact formulation.

[1]  P. Alart,et al.  A mixed formulation for frictional contact problems prone to Newton like solution methods , 1991 .

[2]  Gerhard A. Holzapfel,et al.  Subdivision schemes for smooth contact surfaces of arbitrary mesh topology in 3D , 2004 .

[3]  J. C. Simo,et al.  A perturbed Lagrangian formulation for the finite element solution of contact problems , 1985 .

[4]  C. Bernardi,et al.  A New Nonconforming Approach to Domain Decomposition : The Mortar Element Method , 1994 .

[5]  David J. Benson,et al.  Sliding interfaces with contact-impact in large-scale Lagrangian computations , 1985 .

[6]  Nathan M. Newmark,et al.  A Method of Computation for Structural Dynamics , 1959 .

[7]  David J. Benson,et al.  A single surface contact algorithm for the post-buckling analysis of shell structures , 1990 .

[8]  R. Taylor,et al.  A mixed formulation for the finite element solution of contact problems , 1992 .

[9]  T. Laursen,et al.  A framework for development of surface smoothing procedures in large deformation frictional contact analysis , 2001 .

[10]  Marc Garbey,et al.  Asymptotic and numerical methods for partial differential equations with critical parameters , 1993 .

[11]  Peter Wriggers,et al.  Mortar based frictional contact formulation for higher order interpolations using the moving friction cone , 2006 .

[12]  T. Laursen,et al.  A mortar‐finite element formulation for frictional contact problems , 2000 .

[13]  T. Laursen,et al.  Improved implicit integrators for transient impact problems—geometric admissibility within the conserving framework , 2002, International Journal for Numerical Methods in Engineering.

[14]  M. Puso,et al.  A mortar segment-to-segment contact method for large deformation solid mechanics , 2004 .

[15]  P. W. Christensen,et al.  Formulation and comparison of algorithms for frictional contact problems , 1998 .

[16]  T. Laursen Computational Contact and Impact Mechanics , 2003 .

[17]  Tod A. Laursen,et al.  Formulation and treatment of frictional contact problems using finite elements , 1992 .

[18]  Jacques-Louis Lions,et al.  Nonlinear partial differential equations and their applications : Collège de France seminar , 1982 .

[19]  Kazufumi Ito,et al.  The Primal-Dual Active Set Strategy as a Semismooth Newton Method , 2002, SIAM J. Optim..

[20]  P. Wriggers,et al.  Smooth C1‐interpolations for two‐dimensional frictional contact problems , 2001 .

[21]  M. Crisfield,et al.  Energy‐conserving and decaying Algorithms in non‐linear structural dynamics , 1999 .

[22]  Tod A. Laursen,et al.  Two dimensional mortar contact methods for large deformation frictional sliding , 2005 .

[23]  Anthony T. Patera,et al.  Domain Decomposition by the Mortar Element Method , 1993 .

[24]  M. Puso,et al.  A 3D contact smoothing method using Gregory patches , 2002 .

[25]  Ekkehard Ramm,et al.  Unilateral non‐linear dynamic contact of thin‐walled structures using a primal‐dual active set strategy , 2007 .

[26]  Faker Ben Belgacem,et al.  The mortar finite element method for contact problems , 1998 .

[27]  Ekkehard Ramm,et al.  Generalized Energy–Momentum Method for non-linear adaptive shell dynamics , 1999 .

[28]  Rolf Krause,et al.  Monotone Multigrid Methods on Nonmatching Grids for Nonlinear Multibody Contact Problems , 2003, SIAM J. Sci. Comput..

[29]  Peter Wriggers,et al.  Frictionless 2D Contact formulations for finite deformations based on the mortar method , 2005 .

[30]  M. Schäfer,et al.  A fast and robust iterative solver for nonlinear contact problems using a primal‐dual active set strategy and algebraic multigrid , 2007 .

[31]  Tod A. Laursen,et al.  A mortar segment-to-segment frictional contact method for large deformations , 2003 .

[32]  Peter Wriggers,et al.  Computational Contact Mechanics , 2002 .

[33]  Barbara I. Wohlmuth,et al.  A Mortar Finite Element Method Using Dual Spaces for the Lagrange Multiplier , 2000, SIAM J. Numer. Anal..

[34]  P. Wriggers,et al.  FINITE ELEMENT FORMULATION OF LARGE DEFORMATION IMPACT-CONTACT PROBLEMS WITH FRICTION , 1990 .

[35]  Patrick Hild,et al.  Numerical Implementation of Two Nonconforming Finite Element Methods for Unilateral Contact , 2000 .

[36]  Barbara Wohlmuth,et al.  A primal–dual active set strategy for non-linear multibody contact problems , 2005 .

[37]  Yvon Maday,et al.  Coupling finite element and spectral methods: first results , 1990 .

[38]  P. Wriggers,et al.  A segment-to-segment contact strategy , 1998 .