Nonlinear finite elements/Solution procedure

Solution procedure

edit

The finite element system of equations is

 

Let's assume that the mass matrix is diagonal. Let us also assume that we would like to solve this problem using an explicit method that uses central differencing.

Let the time step size be  . Let us also assume that the time step is constant. Then, at time step  , the time is  .

Let

 

Let us take a half step to compute the velocity at the middle of the timestep. Then,

 

We will use this half-step velocity to compute the acceleration

 

Now,

 

Therefore,

 

Algorithm

edit
  1. Initialize:
    1. Set   and  .
    2. Set the initial velocities ( ) at the nodes.
    3. Set the initial accelerations ( ) at the nodes.
    4. Set the initial stresses ( ) at the element Gauss points.
    5. Compute the lumped mass   at the nodes.
  2. Set the displacements and velocities at the first time step:
    1. Displacement:  
    2. Velocity:  
  3. Apply essential BCs:  
  4. Update the nodal displacements:  
  5. Set   and  .
  6. Loop through the following steps until  .
  7. For each element (at the Gauss points):
    1. Compute the strain measure:  
    2. Compute the stress:  
  8. For each node:
    1. Compute the internal force  .
    2. Compute the external force  .
    3. Compute the total force  .
    4. Compute the acceleration  :  
    5. Update the velocity  .  
    6. Apply the essential boundary conditions.
    7. Update the displacement  .  
  9. Update the counters:  ,  .

Stability of the explicit algorithm

edit

If the time step   is too large, the algorithm may not be stable and may give unreasonable results. To provide a check on the time step size, we use the CFL (Courant-Friedrichs-Lewy) condition to determine  . This condition (in 1-D) states that

 

where   is the initial length of the element, and   is speed of sound in the material (wavespeed) given by