The aim of this chapter is to provide an in‐depth survey of arbitrary Lagrangian–Eulerian (ALE) methods, including both conceptual aspects of the mixed kinematical description and numerical implementation details. Applications are discussed in fluid dynamics, nonlinear solid mechanics, and coupled problems describing fluid–structure interaction. The need for an adequate mesh‐update strategy is underlined, and various automatic mesh‐displacement prescription algorithms are reviewed. This includes mesh‐regularization methods essentially based on geometrical concepts, as well as mesh‐adaptation techniques aimed at optimizing the computational mesh according to some error indicator. Emphasis is then placed on particular issues related to the modeling of compressible and incompressible flow and nonlinear solid mechanics problems. This includes the treatment of convective terms in the conservation equations for mass, momentum, and energy, as well as a discussion of stress‐update procedures for materials with history‐dependent constitutive behavior.

The arbitrary Lagrangian-Eulerian (ALE) description in non-linear solid mechanics is nowadays standard for hypoelastic-plastic models. An extension to hyperelastic-plastic models is presented here. A fractional-step method - a common choice in ALE analysis - is employed for time-marching: every time-step is split into a Lagrangian phase, which accounts for material effects, and a convection phase, where the relative motion between the material and the finite element mesh is considered. In contrast to previous ALE formulations of hyperelasticity or hyperelastoplasticity, the deformed configuration at the beginning of the time-step, not the initial undeformed configuration, is chosen as the reference configuration. As a consequence, convecting variables are required in the description of the elastic response. This is not the case in previous formulations, where only the plastic response contains convection terms. In exchange for the extra convective terms, however, the proposed ALE approach has a major advantage: only the quality of the mesh in the spatial domain must be ensured by the ALE remeshing strategy; in previous formulations, it is also necessary to keep the distortion of the mesh in the material domain under control. Thus, the full potential of the ALE description as an adaptive technique can be exploited here. These aspects are illustrated in detail by means of three numerical examples: a necking test, a coining test and a powder compaction test.Peer ReviewedPostprint (author’s final draft

A key issue in Arbitrary Lagrangian-Eulerian (ALE) non-linear solid mechanics is the correct treatment of the convection terms in the constitutive equation. These convection terms, which re ect the relative motion between the finite element mesh and the material, are found for both transient and quasistatic ALE analyses. It is shown in this paper that the same explicit algorithms can be employed to handle the convection terms of the constitutive equation for both types of analyses. The most attractive consequence of this fact is that a quasistatic simulation can be upgraded from Updated Lagrangian (UL) to ALE without significant extra computational cost. These ideas are illustrated by means of two numerical examples.

We propose a new framework for fracture mechanics, based on the idea of an approximate fracture geometry representation combined with approximate interface conditions. Our approach evolves from the shifted interface method, and introduces the concept of an approximate fracture surface composed of the full edges/faces of an underlying grid that are geometrically close to the true fracture geometry. The original interface conditions are then modified on the surrogate fracture geometry, by way of Taylor expansions. The shifted fracture method does not require cut cell computations or complex data structures, since the behavior of the true fracture is mimicked with standard integrals on the approximate fracture surface. Furthermore, the energetics of the true fracture are represented within the accuracy of the underlying polynomial finite element approximation and independently of the grid topology. The computational framework is presented here in its generality and then applied in the specific context of cohesive zone models, with an extensive set of numerical experiments in two and three dimensions.

In this paper, numerical differentiation is applied to integrate plastic constitutive laws and to compute the corresponding consistent tangent operators. The derivatives of the constitutive equations are approximated by means of difference schemes. These derivatives are needed to achieve quadratic convergence in the integration at Gauss-point level and in the solution of the boundary value problem. Numerical differentiation is shown to be a simple, robust and competitive alternative to analytical derivatives. Quadratic convergence is maintained, provided that adequate schemes and stepsizes are chosen. This point is illustrated by means of some numerical examples.