ACCESS THE COMPLETE LIBRARY OF COURSES & TUTORIALS 
43. A strategy for avoiding a singular stiffness matrix
Drilling degrees of freedom and avoiding singularity
Please log in or enroll to access resources

Summary

In this lecture, we will cover the following:

  • How adding an arbitrary drilling stiffness to the local element stiffness matrix improves the stiffness matrix condition,
  • We examine why the drilling stiffness does not couple into the other degrees of freedom,
  • We reduce the theory to a compact implementation strategy.

We begin by introducing the central fix: adding a small arbitrary stiffness associated with the missing drilling rotation. In the local element stiffness matrix, this stiffness is placed only on the drilling diagonal terms. Because the off-diagonal coupling terms remain zero, the added value prevents the drilling mode from mapping to zero without contaminating the meaningful membrane, bending, or transverse shear behaviour.

We then work through how this local drilling stiffness appears once transformed into the global reference frame. Although the derivation involves expanding the stiffness and transformation matrices, the final implementation is relatively simple: calculate a drilling contribution based on the element normal, place it into the rotational blocks of the element stiffness matrix, and then assemble the structure as usual.

Next up

In the next lecture, we will implement this drilling stabilisation strategy directly in our custom shell solver and check whether it removes the rank deficiency.

Tags

drilling stabilisationarbitrary drilling stiffnessstiffness matrix transformationrank deficiency

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
44. Implementing drilling stabilisation