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
60. Custom mesh Reissner-Mindlin plate analysis
Custom mesh finite element analysis
Please log in or enroll to access resources

Summary

In this lecture, we'll cover the following:

  • Importing a pre-generated custom mesh into a finite element code.
  • Defining support conditions using restrained nodes and degree-of-freedom masks.
  • Constructing the global force vector from element contributions.
  • Assembling the global stiffness matrix from element stiffness matrices.
  • Applying boundary conditions and solving the system for nodal displacements.

In this lecture, we walk through the process of integrating a previously generated mesh into our own custom finite element code and preparing it for analysis. We begin by recreating the model setup, material properties, geometry and mesh, before defining support conditions through a set of restrained nodes. We then translate these supports into restrained degrees of freedom using a masking approach, ensuring that only the relevant vertical displacement constraints are enforced in line with the simplified degrees of freedom per node.

We then move on to assembling the global force vector by looping through each element, computing its equivalent nodal forces, and mapping these contributions into the correct global locations. A similar assembly process is followed for the global stiffness matrix, where each element’s stiffness contribution is computed and inserted using appropriate indexing. After constructing the full system, we apply boundary conditions by reducing both the stiffness matrix and force vector, and finally solve for nodal displacements. This completes the core finite element solution process, setting the stage for post-processing.

Next up

Next, we will compute and visualise the reaction forces from the custom analysis to verify equilibrium and assess the overall response.

Tags

custom mesh integrationboundary condition applicationfinite element solution

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
61. Custom mesh analysis results - reactions