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
47. Max-displacement parameter sweep across models
Benchmarking against OpenSeesPy and Pynite
Please log in or enroll to access resources

Summary

In this lecture, we'll cover the following:

  • Systematic comparison of three plate models across a range of slab thicknesses.
  • Use of maximum mid-span displacement as the primary comparison metric.
  • Implementation of a parameter sweep over thickness.
  • Validation against OpenSeesPy and PyNite solutions as benchmarks.
  • Introduction of Navier’s analytical thin plate solution for additional comparison.
  • Identification of shear locking as the cause of inaccuracy in the custom model.

In this lecture, we explore how our three plate models behave across a wide range of slab thicknesses, moving beyond a single “spot” comparison to a full parameter sweep. We focus on maximum displacement at the plate centre as a consistent and fundamental measure of accuracy, since all other results stem from displacement. By restructuring our code into reusable utility functions, we efficiently run repeated simulations for varying thicknesses and compare outputs from our custom Reissner–Mindlin formulation against OpenSeesPy and PyNite.

We observe that the benchmark models agree closely across all thicknesses, giving us confidence in their validity. Our custom model performs well for thick plates but becomes overly stiff as the plate becomes thin, underpredicting displacement. By also plotting Navier’s analytical solution, we confirm that it aligns well in the thin plate regime but diverges for thicker plates due to neglecting transverse shear effects. This comparison highlights the gradual transition between thin and thick plate behaviour and reveals shear locking as the key issue affecting our custom formulation.

Next up

The parameter sweep has revealed shear locking as the source of the discrepancy. In the next lecture, we will examine what shear locking is and why it occurs.

Tags

parameter sweepReissner–Mindlin plateshear lockingNavier solutionfinite element validation

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
48. The role of shear-locking