Hello, I am having somewhat difficulty understanding the concepts of rank deficiency and hourglassing in finite element methods. Essentially, I have been reading a book outlining this very briefly on half a page and I need a bit more information. As an example:(adsbygoogle = window.adsbygoogle || []).push({});

If we have a 2D elasticity problem with a bilinear element, the element stiffness matrix "k" will be an 8x8 matrix, composed of the integral of B(T) D B, where B is a 3x8 matrix and D is a 3x3 matrix. So far so good. In 2D, this matrix "k" is said to be of rank 5 because you subtract 3 rigid body motions from this 8x8 matrix. While I see that you can have these 3 motions, I don't quite understand why this reduces the rank of the matrix. Could anyone give me some insight into how this works? I know, this probably sounds ridiculously simple, but somehow I am having problems picturing it.

Secondly, when it comes to hourglassing: If you use one Gauss point only to integrate this "k"-matrix numerically, you have: k = 4J B(T)DB | xi, eta = 0, which makes sense. However, this k-matrix has rank 3, the same as the D-matrix, and I don't really know why it isn't 5, like before. This rank deficiency obviously then leads to 2 hourglassing modes. Does anyone have any advice of how to determine, not just here, but also in other cases, the rank of a stiffness matrix? That would be very helpful, thank you.

**Physics Forums - The Fusion of Science and Community**

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# FEM: Rank deficiency and hourglassing

Loading...

Similar Threads - Rank deficiency hourglassing | Date |
---|---|

FEM hourglass control 2D | Jun 16, 2016 |

Steam plant operates with ideal reheat -regenerative Ranking cycle. | May 24, 2011 |

**Physics Forums - The Fusion of Science and Community**