Nonlinear finite elements/Lagrangian and Eulerian descriptions
In our study of the axially loaded bar, we have assumed that the area of cross section remains constant during the deformation. Therefore, the physical stress is not affected by the area of cross section of the rod. However, in most materials, the area of cross section changes with deformation. This change is most prominent during large deformations. In this handout, we will discuss a Lagrangian finite element formulation for large deformations.
There are two main ways of approaching problems that involve the motion of deformable materials - the Lagrangian way and the Eulerian way. These approaches are distinguished by three important aspects:
- The mesh description.
- The stress tensor and momentum equation (kinetics).
- The strain measure (kinematics).
Some essential definitions are given below:
- Spatial or Eulerian coordinates ( ):
- These coordinates are used to locate a point in space with respect to a fixed basis. You can think of these coordinates as the ones you are familiar with.
- Material or Lagrangian coordinates ( ):
- These coordinates are used to label material points. If we sit on a material point, the label does not change with time. We do start with a reference label which we usually choose as the initial spatial coordinates of a material point.
- Motion or Deformation :
- A motion or deformation is defined as a mapping between the initial and the current configuration. We usually write this relationship as
- The initial configuraton is
- Displacement :
- The displacement is defined as the difference between the reference and the current configuration. We write
- Velocity :
- The velocity of a material point is the derivative of the motion with fixed.
- Acceleration :
- The acceleration of a material point is the derivative of the velocity with fixed.
An example edit
Let us consider a one-dimensional problem. Assume that the motion is
The inverse of the map gives us in terms of , i.e.,
Then the displacement of the material point is
The velocity of the material point is
We could alternatively wish to express the velocity in terms of the spatial coordinate . In that case, we have
Note that and give us different ways of expressing the same quantity. To keep the number of symbols reasonably small, we usually represent both functions with the same symbol ( ) and determine which form we are talking about on the basis of the arguments.
We can find the Lagrangian and Eulerian versions of the acceleration in a similar manner. Thus,
Lagrangian and Eulerian Meshes edit
The Lagrangian and Eulerian descriptions can be visualized in terms of the corresponding meshes (see Figure 1).
We can think of the Lagrangian mesh as being drawn on the body. The mesh deforms with the body. Both the nodes and the material points change position as the body deforms. However, the position of the material points relative to the nodes remains fixed.
On the other hand, the Eulerian mesh is a background mesh. The body flows through the mesh as it deforms. The nodes remain fixed and the materials points move through the mesh. The position of a material point relative to the nodes varies with the motion.
Some features, advantages, and disadvantages of the two descriptions are given below.
- This description is based on Nonlinear Finite Elements of Continua and Structures by T. Belytschko, W.K. Liu, and B. Moran.
The Arbitrary Lagrangian-Eulerian (ALE) approach has features of both resolves some of the problems discussed above.