The Finite Element Method: From Integration to Implementation

This self-study workbook provides a first-principles presentation of the Finite Element Method (FEM). Beginning with piecewise numerical integration as a 1-D finite element model, the text establishes the Method of Weighted Residuals, the weak form, and variational principles without unstated mathematical leaps. We derive element formulations for spring networks, 2-D trusses, Euler--Bernoulli flexural beams requiring C¹ Hermite continuity, 2-D Constant Strain Triangles (CST), and isoparametric bilinear quadrilaterals (Q4) with 2× 2 Gauss quadrature. Physical phenomena such as shear locking and advective instability are motivated mechanically before introducing their numerical remedies. The entire computational pipeline---from CAD boundary representation (STEP) and Gmsh unstructured meshing to boundary condition enforcement and post-processing---is demonstrated end-to-end on a stress concentration benchmark. High-performance sparse assembly (Coordinate format to Compressed Sparse Row) and direct linear solvers are implemented from scratch in pure Python. Every chapter features hand-worked numerical examples verified by an independent audit suite, complete exercises with detailed analytical solutions, and extensions to coupled multiphysics, Streamline-Upwind Petrov--Galerkin (SUPG) fluid stabilization, the Extended Finite Element Method (XFEM), and Isogeometric Analysis (IGA).