FINITE ELEMENT ANALYSIS OF SHELLS - EARLY ACCESS 
Section 4
Expanding to a full plate element solver
21. Section overview - Expanding to a full plate element solver
01:28 (Preview)
22. Procedurally generating a rectangular mesh
24:30
23. Defining plate constraints
11:08
24. Defining the self-weight force vector
10:35
25. Building the structure stiffness matrix
10:05
26. Solving the system and extracting reaction forces
28:13
27. Plotting the plate displacements
18:10
28. Building an evaluation grid for stress resultants
10:31
29. Calculating the moments and shears
22:00
30. Visualising the plate bending moments
14:13
31. Extracting shear forces
29:04
32. Visualising the plate shear forces
12:21
33. Adding strip and edge masking to the shear plot
26:04
34. Adding magnitude clipping to the shear plot
10:40
35. Building an interpolation utility function
09:53
15. A primer on numerical integration
Virtual Work and Calculating the Element Stiffness Matrix
Please log in or enroll to access resources

Summary

In this lecture, we'll cover the following:

  • Introduction to numerical integration using Gauss–Legendre (Gauss quadrature)
  • Evaluating integrals in finite element formulations that are difficult to solve analytically
  • Core idea of approximating integrals via weighted sums of function values at selected sampling points
  • Minimising integration error by choosing optimal weights and sample locations
  • General formulation of Gauss quadrature and its role in stiffness matrix computation

In this lecture, we introduce Gauss quadrature as a practical method for numerical integration, particularly in the context of finite element analysis where analytical integration is often infeasible. We explore the fundamental idea of approximating an integral by evaluating a function at selected sampling points, applying appropriate weightings, and summing the results. This is framed through the need to integrate expressions such as those arising in stiffness matrix formulation.

We then focus on how Gauss quadrature achieves high accuracy by approximating the integrand as a polynomial and carefully selecting sampling points and weights to minimise error. Through a worked example with two sampling points, we see how these values are determined and how they relate to the degree of the approximating polynomial. We also outline the general form of the quadrature rule and briefly discuss how the number of sampling points affects accuracy. This method is foundational knowledge that will be applied to compute element stiffness matrices in subsequent lectures.

Next up

In the next lecture, we will apply this numerical integration technique directly to the plate element, transforming the stiffness matrix integral into a form that can be computed in code.

Tags

Gauss quadraturenumerical integrationsampling points and weightsfinite element stiffness matrixpolynomial approximation

Please log in or enroll to continue

If you've already enrolled, please log in to continue.

Finite Element Analysis of Plate and Shell Structures: Part 1 - Plates

An analysis pipeline for thick and thin plate structures, a roadmap from theory to toolbox

After completing this course...

  • You will understand how Reissner-Mindlin theory enables us to accurately capture both thin and thick plate behaviour.
  • You will understand how to turn the fundamental mechanics of plate behaviour into a custom finite element solver written in Python.
  • You will have developed meshing workflows that utilise the powerful open-source meshing engine, GMSH.
  • In addition to using your own custom finite element code, you will be comfortable validating your results using OpenSeesPy and Pynite.
Next Lesson
16. Numerical integration applied to our element