Nonlinear finite elements/Solution procedure

From Wikiversity
Jump to navigation Jump to search

Solution procedure[edit | edit source]

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 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



Algorithm[edit | edit source]

  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 | edit source]

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