FINITE ELEMENT ANALYSIS OF SHELLS - EARLY ACCESS 
9. A review of the element stiffness matrix
Expanding the mechanics for flat shell elements
Please log in or enroll to access resources

Summary

In this lecture, we'll cover:

  • A quick review of the local element stiffness matrix and the force–displacement relationship,
  • How to distinguish between nodal actions and equivalent nodal actions from distributed loading,
  • How the element stiffness matrix is formed using the BTDBB^T D B integral over the element area,
  • How equivalent nodal actions are obtained from distributed forces and moments using shape functions,
  • How the stiffness matrix can be split into membrane, bending, and shear contributions,

In this lecture, we revisit the local element stiffness matrix equation, mainly to refresh your memory before moving on to coding. We start with the familiar force equals stiffness times displacement relationship, where the local stiffness matrix acts on the local nodal displacement vector and the right-hand side includes both directly applied nodal actions and equivalent nodal actions arising from distributed loads such as self-weight. We then restate the stiffness matrix as the area integral of BTDBB^T D B, noting that the thickness integration has already been handled when deriving the constitutive matrix.

We then break the formulation into its membrane, bending, and shear parts. Because the constitutive matrix is effectively block diagonal for the homogeneous material considered here, the element stiffness matrix can be written as a simple superposition of independent membrane, bending, and shear stiffness matrices, with no coupling at the element level. We also identify the equivalent nodal force vector as the area integral of the shape function matrix multiplied by the vector of distributed actions.

Next up

In the next lecture, we will look at how the stiffness matrix integral is actually computed in practice, using the Jacobian and Gauss numerical integration to evaluate the expression at a set of sampling points.

Tags

membrane stiffness matrixequivalent nodal actionsshell element formulationshape function matrixnumerical integration setup

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 2 - Shells

Expanding from plate to shell elements - build a workflow that unlocks the behaviour of 3D shell structures

After completing this course...

  • You will understand how we make the leap from Reissner-Mindlin plate elements to shell elements and what extra modelling fidelity that provides.
  • You will be comfortable using a combination of GMSH and the open-source 3D modelling software, Blender, to generate custom finite element meshes.
  • You will be able to use OpenSeesPy to model shell structures, as an alternative to your own custom finite element solver.
  • You will have a much greater understanding of what commercial finite element packages are doing, behind the UI, allowing you to authoritatively interrogate their results.
Next Lesson
10. Quick recap of numerical integration and the Jacobian