Skip to main content

Foundation of FEM; choosing a numerical method

Short answer

FEM is founded on potential energy stationarity, using piecewise shape functions to form [K]{U}={F}. It is less suitable for large deformation because Lagrangian mesh distortion causes element tangling and accuracy loss. Method choice depends on geometry and material continuity; e.g., FEM for structural stress, DEM for granular flow.

Explain the foundation of the finite element method. Why is this method less appropriate for large deformation problems? How do you choose a numerical method for different problems? Illustrate with examples.

Foundation of the finite element method

FEM rests on the idea that a complex continuum can be approximated by a collection of simple pieces (elements) whose behaviour is known. Its mathematical basis is:

  1. Variational / energy principle: the correct displacement field makes the total potential energy Π=U−W\Pi=U-W stationary (minimum), or equivalently satisfies the virtual work equation. Weighted residual (Galerkin) methods can also be used.
  2. Piecewise approximation: displacement in each element is u=[N]{d}u=[N]\{d\} with shape functions that are continuous across elements.
  3. Element equations: from {ε}=[B]{d}\{\varepsilon\}=[B]\{d\} and {σ}=[D]{ε}\{\sigma\}=[D]\{\varepsilon\}, [k]=∫[B]T[D][B] dV[k]=\int [B]^T[D][B]\,dV.
  4. Assembly and solution: compatibility at nodes and equilibrium give [K]{U}={F}[K]\{U\}=\{F\}.

Why FEM is less suitable for large deformation

  • Elements are attached to the material (Lagrangian); under very large strain the mesh distorts, elements become tangled or negative in volume, and accuracy falls.
  • Small-strain formulation [B][B] and constant [K][K] assume small displacement; for large deformation the stiffness changes with geometry and needs non-linear incremental solution and re-meshing.
  • Material separation, fracture and flow are not represented by continuous elements. Methods such as DEM, SPH or ALE formulations are better.

Choosing a numerical method

ProblemSuitable method
Stress in frame, plate, damFEM
Seepage, unsteady channel/pipe flow on a gridFDM
Infinite/semi-infinite soil, tunnel, crackBEM
Granular soil, jointed rock, rockfallDEM
Dam-break, debris flow, large deformationSPH / mesh-free

Selection depends on geometry, material behaviour (linear or non-linear, continuum or discontinuum), boundary conditions, required accuracy and computing resources.

✓ Answer checked by Claude Opus 5.5 (high effort)

Related questions

All Computational Techniques in Civil Engineering chapter 1 questions →