Unilateral non‐linear dynamic contact of thin‐walled structures using a primal‐dual active set strategy

The efficient modelling of three-dimensional contact problems is still a challenge in non-linear implicit structural analysis. We use a primal-dual active set strategy (SIAM J. Optim. 2003; 13:865-888), based on dual Lagrange multipliers (SIAM J. Numer. Anal. 2000; 38:989-1012) to handle the non-linearity of the contact conditions. This allows us to enforce the contact constraints in a weak, integral sense without any additional parameter. Due to the biorthogonality condition of the basis functions, the Lagrange multipliers can be locally eliminated. We perform a static condensation to achieve a reduced system for the displacements. The Lagrange multipliers, representing the contact pressure, can be easily recovered from the displacements in a variationally consistent way. For the application to thin-walled structures we adapt a three-dimensional non-linear shell formulation including the thickness stretch of the shell to contact problems. A reparametrization of the geometric description of the shell body gives us a surface-oriented shell element, which allows the application of contact conditions directly to nodes lying on the contact surface. Shell typical locking phenomena are treated with the enhanced-assumed-strain-method and the assumed-natural-strain-method. The discretization in time is done with the implicit Generalized-α method (J. Appl. Mech. 1993; 60:371-375) and the Generalized Energy-Momentum Method (Comp. Methods Appl. Mech. Eng. 1999; 178:343-366) to compare the development of energies within a frictionless contact description. In order to conserve the total energy within the discretized frictionless contact framework, we follow an approach from Laursen and Love (Int. J. Numer. Methods Eng. 2002; 53:245-274), who introduced a discrete contact velocity to update the velocity field in a post-processing step. Various examples show the good performance of the primal-dual active set strategy applied to the implicit dynamic analysis of thin-walled structures.

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

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

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

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

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

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

[7]  E. Ramm,et al.  On the physical significance of higher order kinematic and static variables in a three-dimensional shell formulation , 2000 .

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

[9]  T. Laursen,et al.  DESIGN OF ENERGY CONSERVING ALGORITHMS FOR FRICTIONLESS DYNAMIC CONTACT PROBLEMS , 1997 .

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

[11]  J. Oden,et al.  Contact Problems in Elasticity: A Study of Variational Inequalities and Finite Element Methods , 1987 .

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

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

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

[15]  E. Ramm,et al.  Three‐dimensional extension of non‐linear shell formulation based on the enhanced assumed strain concept , 1994 .

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

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

[18]  F. Armero,et al.  Formulation and analysis of conserving algorithms for frictionless dynamic contact/impact problems , 1998 .

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

[20]  Wolfgang A. Wall,et al.  Parallel multilevel solution of nonlinear shell structures , 2005 .

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

[22]  Barbara Wohlmuth,et al.  Stable Lagrange multipliers for quadrilateral meshes of curved interfaces in 3D , 2007 .

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

[24]  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.

[25]  P. Tallec,et al.  Energy-controlling time integration methods for nonlinear elastodynamics and low-velocity impact , 2006 .

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

[27]  J. C. Simo,et al.  An augmented lagrangian treatment of contact problems involving friction , 1992 .

[28]  Jintai Chung,et al.  A Time Integration Algorithm for Structural Dynamics With Improved Numerical Dissipation: The Generalized-α Method , 1993 .

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

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