This paper presents a complete derivation and implementation of the Arbitrary Lagrangian Eulerian (ALE) formulation for the simulation of large deformation quasi-static and dynamic problems. While most of the previous work done on ALE for dynamic applications was mainly based on operator split and explicit calculations, this work derives the quasi-static and dynamic ALE equations using a fully coupled implicit approach. Full expression for the ALE virtual work equations and finite element matrices are given. Time integration relations for the dynamic equations are also derived. Several quasistatic and dynamic large deformation applications are solved and presented.
IntroductionThe Arbitrary Lagrangian Eulerian (ALE) formulation has emerged in recent years as a technique that can alleviate many of the drawbacks of the traditional Lagrangian and Eulerian formulations [2,5,7,11,16]. Using ALE, the computational grid need not adhere to the material nor be fixed in space but can be moved arbitrarily. The grid is continuously moved to optimize element shapes independently from material deformation. A proper ALE formulation should reduce to both the Lagrangian and Eulerian formulations at any degree of freedom as desired. Combining the merits of both the Lagrangian and Eulerian formulations, ALE can easily describe different types of boundary conditions and prevent mesh distortion.The ALE equations are derived by substituting the relationship between the material time derivative and grid time derivative into the continuum mechanics governing equations. This substitution gives rise to convective terms in the ALE equations which account for the transport of material through the grid. ALE is usually termed a coupled formulation since material deformation and convective effects are coupled in the same equations. A survey of the ALE literature [17], however, shows that the majority of ALE analyses, whether quasi-static or dynamic, are based on the computationally convenient operator split technique. In this approach, material deformation and convective effects are treated separately. Thus each time step may be divided into two steps: a regular Lagrangian step followed by an Eulerian step. The main advantage of this technique over the fully coupled approach is the reduction in the cost of implementation of ALE to current Lagrangian codes as the Lagrangian step is unchanged and only the Eulerian step algorithm needs to be added. Moreover, the decoupling of the Lagrangian and Eulerian steps results in simpler equations to be solved. However, from the theoretical point of view, the fully coupled ALE approach represents a true kinematic description in which material deformation is described relative to a moving reference configuration.In this work, a complete treatment for the fully coupled implicit ALE formulation is presented. Virtual work equations are first derived from the basic principles of continuum mechanics. Next, finite element descritization of the virtual work equations is performed. Full expression for the resulting...