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
63. What's going on with our deflection?
Custom mesh finite element analysis
Please log in or enroll to access resources

Summary

In this lecture, we'll cover the following:

  • Why the displacement patterns from the previous lecture occur.
  • How selective (reduced) integration for shear introduces zero-energy modes.
  • The role of the strain–displacement matrix and single-point sampling.
  • What it means for a stiffness matrix to be rank deficient.
  • Physical interpretation of zero-energy (hourglass) displacement patterns.

In this lecture, we explore the origin of the unusual zigzag displacement patterns observed previously and identify them as zero-energy modes caused by selective integration. We revisit how the shear stiffness matrix is constructed and show that using a one-point Gauss integration scheme effectively samples the strain–displacement matrix only at the element centre. This assumption forces the transverse shear strain to appear constant across the element, even though in reality it can vary, leading to a loss of critical information.

We then see that this under-integration causes the element stiffness matrix to become rank deficient, meaning it admits non-zero displacement patterns without any applied force. These non-physical modes, also known as hourglass modes, correspond to deformations with no associated strain energy. Through simple element and multi-element examples, we build intuition for how these patterns arise and why they manifest as characteristic zigzag shapes. We conclude that although selective integration alleviates shear locking, it introduces these spurious modes, motivating the need for an improved integration strategy in the next lecture.

Next up

In the next lecture, we will develop a detailed plan for fixing the Reissner–Mindlin element by addressing the root cause of shear locking without reintroducing zero-energy modes.

Tags

zero-energy modesselective integrationshear lockingrank deficiencyhourglass modes

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
64. A plan for fixing our Reissner-Mindlin element