A formulation of conserving impact system based on localized Lagrange multipliers

This paper proposes a new finite element method of frictional impact of elastic bodies. The formulation introduces a contact frame that is placed in between contacting bodies and represents the contact surface. The nonpenetration condition and the slip–stick condition are defined between the contacting body and the contact frame with the aid of the independent localized Lagrange multipliers representing the contact forces. The position of the contact frame and the local coordinate of the contacting node along the contact frame are also treated as the independent variable, which enables the exact satisfaction of the constraint conditions without the deficiency or redundant constraint. The energy and momentum conservation algorithm is applied to the proposed impact system. For the case of frictional impact, the linear momentum is exactly conserved and the angular momentum is approximately conserved with negligible error. Copyright © 2006 John Wiley & Sons, Ltd.     

A finite deformation mortar contact formulation using a primal–dual active set strategy

In recent years, nonconforming domain decomposition techniques and, in particular, the mortar method have become popular in developing new contact algorithms. Here, we present an approach for 2D frictionless multibody contact based on a mortar formulation and using a primal–dual active set strategy for contact constraint enforcement. We consider linear and higher‐order (quadratic) interpolations throughout this work. So‐called dual Lagrange multipliers are introduced for the contact pressure but can be eliminated from the global system of equations by static condensation, thus avoiding an increase in system size. For this type of contact formulation, we provide a full linearization of both contact forces and normal (non‐penetration) and tangential (frictionless sliding) contact constraints in the finite deformation frame. The necessity of such a linearization in order to obtain a consistent Newton scheme is demonstrated. By further interpreting the active set search as a semi‐smooth Newton method, contact nonlinearity and geometrical and material nonlinearity can be resolved within one single iterative scheme. This yields a robust and highly efficient algorithm for frictionless finite deformation contact problems. Numerical examples illustrate the efficiency of our method and the high quality of results. Copyright © 2009 John Wiley & Sons, Ltd.     

An implicit finite wear contact formulation based on dual mortar methods

Finite deformation contact problems with frictional effects and finite shape changes due to wear are investigated. To capture the finite shape changes, a third configuration besides the well‐known reference and spatial configurations is introduced, which represents the time‐dependent worn state. Consistent interconnections between these states are realized by an arbitrary Lagrangean–Eulerian formulation. The newly developed partitioned and fully implicit algorithm is based on a Lagrangean step and a shape evolution step. Within the Lagrangean step, contact constraints as well as the wear equations are weakly enforced following the well‐established mortar framework. Additional unknowns due to the employed Lagrange multiplier method for contact constraint enforcement and due to wear itself are eliminated by condensation procedures based on the concept of biorthogonality and the so‐called dual shape functions. Several numerical examples in both 2D and 3D are provided to demonstrate the performance and accuracy of the proposed numerical algorithm. Copyright © 2016 John Wiley & Sons, Ltd.     

Development of a finite element contact analysis algorithm for charged‐hydrated soft tissues with large sliding

This paper focuses on a finite element analysis of contact phenomena with large sliding between charged‐hydrated biological soft tissues, such as articular cartilages, based on the triphasic theory. The impenetrability constraint between the contacting bodies and the continuity of the interstitial fluid and ion phases at the contact surfaces are imposed by applying a Lagrange multiplier approach with the contact pressure, chemical potential of the fluid and electrochemical potentials of ions as Lagrange multipliers. A node‐to‐segment one‐pass approach is adopted to cope with large deformations and sliding between the contact surfaces. To pass the contact patch test, contact boundary integrations are performed on both the master and slave contact surfaces. On the other hand, the degrees of freedom of the multipliers at the master nodes are eliminated by projecting the master nodes onto the slave surface to avoid overconstraint. The effectiveness of the proposed algorithm is verified by a couple of numerical examples, in which continuous distributions of displacement, fluid flow, ionic molar flow and Lagrange multipliers on or across the contact surface are confirmed. Copyright © 2008 John Wiley & Sons, Ltd.     

A large deformation frictional contact formulation using NURBS‐based isogeometric analysis

This paper focuses on the application of NURBS‐based isogeometric analysis to Coulomb frictional contact problems between deformable bodies, in the context of large deformations. A mortar‐based approach is presented to treat the contact constraints, whereby the discretization of the continuum is performed with arbitrary order NURBS, as well as C⁰‐continuous Lagrange polynomial elements for comparison purposes. The numerical examples show that the proposed contact formulation in conjunction with the NURBS discretization delivers accurate and robust predictions. Results of lower quality are obtained from the Lagrange discretization, as well as from a different contact formulation based on the enforcement of the contact constraints at every integration point on the contact surface. Copyright © 2011 John Wiley & Sons, Ltd.     

Finite deformation frictional mortar contact using a semi‐smooth Newton method with consistent linearization

A two‐dimensional, finite deformation frictional contact formulation with Coulomb's law is presented. The approach considers multibody contact and is based on a mortar formulation. The enforcement of contact constraints is realized with dual Lagrange multipliers. These alternative multiplier spaces are constructed in a way that the multipliers can easily be eliminated from the global system of equations by static condensation such that the system size does not increase. Friction kinematic variables are formulated in an objective way and enter non‐smooth complementarity functions for expressing the contact constraints. An active set strategy is derived by applying a semi‐smooth Newton method, which treats contact nonlinearities, material and geometrical nonlinearities in one single iterative scheme. By further carrying out a consistent linearization for both normal and frictional contact forces and constraints, a robust and highly efficient algorithm for linear and higher‐order (quadratic) interpolation is achieved. Efficiency of the proposed method and quality of results are demonstrated in several examples. Copyright © 2010 John Wiley & Sons, Ltd.     

Improved implicit integrators for transient impact problems—geometric admissibility within the conserving framework

The value of energy and momentum conserving algorithms has been well established for the analysis of highly non‐linear systems, including those characterized by the nonsmooth non‐linearities of an impact event. This work proposes an improved integration scheme for frictionless dynamic contact, seeking to preserve the stability properties of exact energy and momentum conservation without the heretofore unavoidable compromise of violating geometric admissibility as established by the contact constraints. The physically motivated introduction of a discrete contact velocity provides an algorithmic framework that ensures exact conservation locally while remaining independent of the choice of constraint treatment, thus making full conservation equally possible in conjunction with a penalty regularization as with an exact Lagrange multiplier enforcement. The discrete velocity effects are incorporated as a post‐convergence update to the system velocities, and thus have no direct effect on the non‐linear solution of the displacement equilibrium equation. The result is a robust implicit algorithmic treatment of dynamic frictionless impact, appropriate for large deformations and fully conservative for a range of geometric constraints. Copyright © 2001 John Wiley & Sons, Ltd.     

A finite element method for contact using a third medium

The numerical simulation of contact problems is nowadays a standard procedure in many engineering applications. The contact constraints are usually formulated using either the Lagrange multiplier, the penalty approach or variants of both methodologies. The aim of this paper is to introduce a new scheme that is based on a space filling mesh in which the contacting bodies can move and interact. To be able to account for the contact constraints, the property of the medium, that imbeds the bodies coming into contact, has to change with respect to the movements of the bodies. Within this approach the medium will be formulated as an isotropic/anisotropic material with changing characteristics and directions. In this paper we will derive a new finite element formulation that is based on the above mentioned ideas. The formulation is presented for large deformation analysis and frictionless contact.     

Development of contact algorithm for three‐dimensional numerical manifold method

This paper customizes a contact detection and enforcing scheme to fit the three‐dimensional (3‐D) numerical manifold method (NMM). A hierarchical contact system is established for efficient contact detection. The mathematical mesh, a unique component in the NMM, is utilized for global searching of possible contact blocks and elements, followed by the local searching to identify primitive hierarchies. All the potential contact pairs are then transformed into one of the two essential entrance modes: point‐to‐plane and crossing‐lines modes, among which real contact pairs are detected through a unified formula. The penalty method is selected to enforce the contact constraints, and a general contact solution procedure in the 3‐D NMM is established. Because of the implicit framework, an open‐close iteration is performed within each time step to determine the correct number of contact pairs among multi‐bodies and to achieve complete convergence of imposed contact force at corresponding position. The proposed contact algorithm extensively utilizes most of the original components of the NMM, namely, the mathematical mesh/cells and the manifold elements, as well as the external components associated with contacts, such as the contact body, the contact facet and the contact vertex. In particular, the utilization of two mutually approaching mathematical cells is efficient in detecting contacting territory, which makes this method particularly effective for both convex and non‐convex bodies. The validity and accuracy of the proposed contact algorithm are verified and demonstrated through three benchmark problems. Copyright © 2013 John Wiley & Sons, Ltd.     

A general conjugate‐gradient‐based predictor–corrector solver for explicit finite‐element contact

A Lagrange‐multiplier based approach is presented for the general solution of multi‐body contact within an explicit finite element framework. The technique employs an explicit predictor step to permit the detection of interpenetration and then utilizes a corrector step, whose solution is obtained with a pre‐conditioned matrix‐free conjugate gradient projection method, to determine the Lagrange multipliers necessary to eliminate the predicted penetration. The predictor–corrector algorithm is developed for deformable bodies based upon the central difference method, and for rigid bodies from momentum and energy conserving approaches. Both frictionless and Coulomb‐based frictional contact idealizations are addressed. The technique imposes no time‐step constraints and quickly mitigates velocity discontinuities across closed interfaces. Special attention is directed toward contact between rigid bodies. Algorithmic moment arms conserve the translational and angular momentums of the system in the absence of external loads. Elastic collisions are captured with a two‐phase predictor–corrector approach and a geometrically approximate velocity jump criterion. The first step solves the inelastic contact problem and identifies inactive constraints between rigid bodies, while the second step generates the necessary velocity jump condition on the active constraints. The velocity criterion is shown to algorithmically preserve the system kinetic energy for two unconstrained rigid bodies. Copyright © 1999 John Wiley & Sons, Ltd. This paper was produced under the auspices of the U.S. Government and it is therefore not subject to copyright in the U.S.     

A yield‐limited Lagrange multiplier formulation for frictional contact

An analogy with rigid plasticity is used to develop a constitutive framework for quasi‐static frictional contact between finitely deforming solids. Within this setting, a Lagrange multiplier method is used to impose a sharp distinction between stick and slip. The scope of the multipliers is limited by a constitutively defined ‘yield’ function and a finite element‐based predictor–corrector scheme is employed to efficiently determine the regions of stick and slip and the associated tractions. Selected simulations of planar quasi‐static problems are presented to validate the method and illustrate its capabilities. Copyright © 2000 John Wiley & Sons, Ltd.     

A dual mortar approach for 3D finite deformation contact with consistent linearization

In this paper, an approach for three‐dimensional frictionless contact based on a dual mortar formulation and using a primal–dual active set strategy for direct constraint enforcement is presented. We focus on linear shape functions, but briefly address higher order interpolation as well. The study builds on previous work by the authors for two‐dimensional problems. First and foremost, the ideas of a consistently linearized dual mortar scheme and of an interpretation of the active set search as a semi‐smooth Newton method are extended to the 3D case. This allows for solving all types of nonlinearities (i.e. geometrical, material and contact) within one single Newton scheme. Owing to the dual Lagrange multiplier approach employed, this advantage is not accompanied by an undesirable increase in system size as the Lagrange multipliers can be condensed from the global system of equations. Moreover, it is pointed out that the presented method does not make use of any regularization of contact constraints. Numerical examples illustrate the efficiency of our method and the high quality of results in 3D finite deformation contact analysis. Copyright © 2010 John Wiley & Sons, Ltd.     

Confidence structural robust optimization by non‐linear semidefinite programming‐based single‐level formulation

Structural robust optimization problems are often solved via the so‐called Bi‐level approach. This solution procedure often involves large computational efforts and sometimes its convergence properties are not so good because of the non‐smooth nature of the Bi‐level formulation. Another problem associated with the traditional Bi‐level approach is that the confidence of the robustness of the obtained solutions cannot be fully assured at least theoretically. In the present paper, confidence single‐level non‐linear semidefinite programming (NLSDP) formulations for structural robust optimization problems under stiffness uncertainties are proposed. This is achieved by using some tools such as S‐procedure and quadratic embedding for convex analysis. The resulted NLSDP problems are solved using the modified augmented Lagrange multiplier method which has sound mathematical properties. Numerical examples show that confidence robust optimal solutions can be obtained with the proposed approach effectively. Copyright © 2010 John Wiley & Sons, Ltd.     

Algorithms for the analysis of 3D finite strain contact problems

From the constraint imposition aspects in 3D to friction regularization, various ideas are exposed in this paper. A variation of the Rockafellar Lagrangian is proposed which results in continuous second‐order derivatives if Lagrange multiplier estimates are greater or equal than one. This fact allows the adoption of a full second‐order (i.e. Lagrange–Newton) method avoiding sequential unconstrained minimization techniques. An algorithm for global and local contact detection is presented which is developed for dealing with large step sizes typical of implicit methods. A modified constraint definition to deal with non‐smooth situations is presented. Aspects of friction implementation, including a regularization scheme which ensures stepwise objectivity, are detailed. Finally, several illustrative examples are carried out with success. Copyright © 2004 John Wiley & Sons, Ltd.     

Decomposition contact response (DCR) for explicit finite element dynamics

We propose a new explicit contact algorithm for finite element discretized solids and shells with smooth and non‐smooth geometries. The equations of motion are integrated in time with a predictor‐corrector‐type algorithm. After each predictor step, the impenetrability constraints and the exchange of momenta between the impacting bodies are considered and enforced independently. The geometrically inadmissible penetrations are removed using closest point projections or similar updates. Penetration is measured using the signed volume of intersection described by the contacting surface elements, which is well‐defined for both smooth and non‐smooth geometries. For computing the instantaneous velocity changes that occur during the impact event, we introduce the decomposition contact response method. This enables the closed‐form solution of the jump equations at impact, and applies to non‐frictional as well as frictional contact, as exemplified by the Coulomb frictional model. The overall algorithm has excellent momentum and energy conservation characteristics, as several numerical examples demonstrate. Copyright © 2005 John Wiley & Sons, Ltd.     

DESIGN OF ENERGY CONSERVING ALGORITHMS FOR FRICTIONLESS DYNAMIC CONTACT PROBLEMS

This paper proposes a formulation of dynamic contact problems which enables exact algorithmic conservation of linear momentum, angular momentum, and energy in finite element simulations. It is seen that a Lagrange multiplier enforcement of an appropriate contact rate constraint produces these conservation properties. A related method is presented in which a penalty regularization of the aforementioned rate constraint is utilized. This penalty method sacrifices the energy conservation property, but is dissipative under all conditions of changing contact so that the global algorithm remains stable. Notably, it is also shown that augmented Lagrangian iteration utilizing this penalty kernel reproduces the energy conserving (i.e. Lagrange multiplier) solution to any desired degree of accuracy. The result is a robust, stable method even in the context of large deformations, as is shown by some representative numerical examples. In particular, the ability of the formulation to produce accurate results where more traditional integration schemes fail is emphasized by the numerical simulations. © 1997 by John Wiley & Sons, Ltd.     

Interface problems with quadratic X‐FEM: design of a stable multiplier space and error analysis

The aim of this paper is to propose a procedure to accurately compute curved interfaces problems within the extended finite element method and with quadratic elements. It is dedicated to gradient discontinuous problems, which cover the case of bimaterials as the main application. We focus on the use of Lagrange multipliers to enforce adherence at the interface, which makes this strategy applicable to cohesive laws or unilateral contact. Convergence then occurs under the condition that a discrete inf‐sup condition is passed. A dedicated P1 multiplier space intended for use with P2 displacements is introduced. Analytical proof that it passes the inf‐sup condition is presented in the two‐dimensional case. Under the assumption that this inf‐sup condition holds, a priori error estimates are derived for linear or quadratic elements as functions of the curved interface resolution and of the interpolation properties of the discrete Lagrange multipliers space. The estimates are successfully checked against several numerical experiments: disparities, when they occur, are explained in the literature. Besides, the new multiplier space is able to produce quadratic convergence from P2 displacements and quadratic geometry resolution. Copyright © 2014 John Wiley & Sons, Ltd.     

A Nitsche‐type formulation and comparison of the most common domain decomposition methods in isogeometric analysis

This paper provides a detailed elaboration and assessment of the most common domain decomposition methods for their application in isogeometric analysis. The methods comprise a penalty approach, Lagrange multiplier methods, and a Nitsche‐type method. For the Nitsche method, a new stabilized formulation is developed in the context of isogeometric analysis to guarantee coercivity. All these methods are investigated on problems of linear elasticity and eigenfrequency analysis in 2D. In particular, focus is put on non‐uniform rational B‐spline patches which join nonconformingly along their common interface. Thus, the application of isogeometric analysis is extended to multi‐patches, which can have an arbitrary parametrization on the adjacent edges. Moreover, it has been shown that the unique properties provided by isogeometric analysis, that is, high‐order functions and smoothness across the element boundaries, carry over for the analysis of multiple domains. Copyright © 2013 John Wiley & Sons, Ltd.     

A stable energy‐conserving approach for frictional contact problems based on quadrature formulas

A common approach for the numerical simulation of non‐linear multi‐body contact problems is the use of Lagrange multipliers to model the contact conditions. The stability of standard algorithms is improved by introducing a modified mass matrix which assigns no mass to the potential contact nodes. By this, the spurious algorithmic oscillations in the multiplier do not occur any more, which facilitates the application of the primal–dual active set strategy to dynamical contact problems. The new mass matrix is calculated via a modified quadrature formula that needs no extra computational cost. In addition the conservation properties of the underlying algorithm are transferred to the modified mass version. Different numerical examples for frictional two‐body contact problems illustrate the improvement in the results for the contact stresses. Copyright © 2007 John Wiley & Sons, Ltd.