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
40. Visualising OpenSeesPy moments
Benchmarking against OpenSeesPy and Pynite
Please log in or enroll to access resources

Summary

In this lecture, we'll cover the following:

  • Visualising heat maps on irregular grids using triangulation.
  • Using matplotlib.tri and tricontourf to generate filled contour plots.
  • Handling different bending moment components (Mx,My,Mxy)(M_x, M_y, M_{xy}).
  • Adding interactive widgets to control plotted data and display options.
  • Overlaying triangulation, Gauss points, and element grids for interpretation.

In this lecture, we focus on how to visualise data defined on an irregular grid, which commonly arises when working with finite element meshes. We introduce the use of triangulation via matplotlib.tri, allowing us to convert irregularly spaced data into a form suitable for contour plotting. By using tricontourf, we can generate smooth, filled contour plots that effectively represent quantities such as bending moments across the structure.

We implement an interactive plotting function that allows us to switch between different moment components (Mx,My,Mxy)(M_x, M_y, M_{xy}), while also controlling visual features such as triangulation, Gauss points, and the element grid. We explore how the data is prepared differently depending on the moment type, particularly noting the treatment of absolute values for bending moments versus torsional moments. Through this, we gain insight into how results from analysis (e.g. Gauss point data) map onto the physical structure.

Finally, we enhance the visualisation by overlaying triangulation lines and sampling points to better understand how the irregular grid is constructed and how results are distributed. This provides a clear and flexible framework for interpreting finite element results, setting us up to extend the same approach to shear forces in the next lecture.

Next up

Next, we will extend the same visualisation approach to transverse shear forces, incorporating masking techniques to focus on relevant regions.

Tags

triangulationtricontourfirregular grid visualisationbending momentsmatplotlib widgets

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
41. Visualising OpenSeesPy shears