An augmented Lagrangian algorithm for contact mechanics based on linear regression

SUMMARY The penalty method for the solution of contact problems yields an approximate satisfaction of the contact constraints. Augmentation schemes, which can be adopted to improve the solution, either include the contact forces as additional unknowns or are strongly affected by the penalty parameter and display a poor convergence rate. In a previous investigation, an unconventional augmentation scheme was proposed, on the basis of estimating the ‘exact’ values of the contact forces through linear interpolation of a database extracted by the previous converged states. An enhanced version of this method is proposed herein. With respect to the original method, the enhanced one eliminates some numerical problems and improves the regularity of the convergence path by carrying out the estimate through linear regression methods. The resulting convergence rate is superlinear, and the method is quite insensitive to the penalty parameter. The main underlying concept is that, within the iterative solution of a non-linear problem, linear regression techniques may be used as a tool to ‘shoot’ faster to the final solution, on the basis of a set of intermediate data. The generality of this concept makes it potentially applicable to contact problems in more general settings, as well as to other categories of non-linear problems. Copyright © 2012 John Wiley & Sons, Ltd.

[1]  P. Wriggers.  Computational contact mechanics , 2012 .

[2]  P. Wriggers,et al.  A superlinear convergent augmented Lagrangian procedure for contact problems , 1999 .

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

[4]  A. Curnier,et al.  Large deformation frictional contact mechanics: continuum formulation and augmented Lagrangian treatment , 1999 .

[5]  J. C. Simo,et al.  Algorithms for static and dynamic multiplicative plasticity that preserve the classical return mapping schemes of the infinitesimal theory , 1992 .

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

[7]  Dimitri P. Bertsekas,et al.  Constrained Optimization and Lagrange Multiplier Methods , 1982 .

[8]  Mokhtar S. Bazaraa,et al.  Nonlinear Programming: Theory and Algorithms , 1993 .

[9]  Wolfgang A. Wall,et al.  A dual mortar approach for 3D finite deformation contact with consistent linearization , 2010 .

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

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

[12]  P. Wriggers,et al.  On augmented Lagrangian algorithms for thermomechanical contact problems with friction , 1995 .

[13]  Peter Wriggers,et al.  A note on tangent stiffness for fully nonlinear contact problems , 1985 .

[14]  Barbara I. Wohlmuth,et al.  A Primal-Dual Active Set Algorithm for Three-Dimensional Contact Problems with Coulomb Friction , 2008, SIAM J. Sci. Comput..

[15]  P. Wriggers,et al.  Application of augmented Lagrangian techniques for non‐linear constitutive laws in contact interfaces , 1993 .

[16]  J. Z. Zhu,et al.  The finite element method , 1977 .