ACCESS THE COMPLETE LIBRARY OF COURSES & TUTORIALS 
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
48. The role of shear-locking
Benchmarking against OpenSeesPy and Pynite
Please log in or enroll to access resources

Summary

In this lecture, we'll cover the following:

  • What shear locking is and how it manifests in finite element results.
  • The relationship between bending and shear stiffness and their dependence on plate thickness.
  • How shear stiffness dominates as thickness approaches zero, causing locking.
  • The concept of selective integration as a remedy and how it is implemented.
  • Practical rules for applying selective integration in plate elements.
  • Potential issues such as spurious zero-energy modes (hourglassing).

In this lecture, we examine shear locking, a phenomenon where finite elements become artificially too stiff, leading to underestimated deflections. We clarify that this is not a numerical instability, as solutions still converge, but rather a modelling issue that produces incorrect results. By analysing the stiffness formulation, we see that bending stiffness scales with thickness cubed, while shear stiffness scales linearly with thickness. As the plate becomes very thin, shear stiffness dominates, effectively suppressing transverse shear deformation and causing the locking effect.

We then explore selective integration as a practical and effective solution. By reducing the integration order for the shear component (while retaining full integration for bending), we intentionally weaken the shear contribution, preventing it from overwhelming the response as thickness decreases. This approach works well across both thin and thick plates, forming a reliable default strategy. However, we also note a potential drawback: the introduction of non-physical, zero-energy deformation modes (hourglassing), which must be recognised and managed in more advanced analyses.

Next up

In the next lecture, we will implement the selective integration fix for shear locking and re-run the comparison to confirm that the issue has been resolved.

Tags

shear lockingReissner-Mindlin platesselective integrationGauss quadraturehourglassing

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
49. Adapting for shear-locking and comparing again