3

If a member is very slender then its length will be large compared to its cross-sectional area. I've heard and read that slender members can potentially cause the stiffness matrix of a beam/frame element to become ill-conditioned in FEM.

Question: How/why can a slender member cause the stiffness matrix to become ill-conditioned?

Here are my own thoughts regarding this question, which may or may not be correct - I'd like some confirmation. A long slender member will have larger axial stiffness than flexural stiffness. In the stiffness matrix, the axial stiffness term has the length of the member ($L$) in the denominator. However, the flexural stiffness terms typically have $L^2$ and $L^3$ in their denominators. Therefore, given constant cross-sectional properties ($A$ and $I$), as $L$ grows the flexural stiffness terms approach zero much quicker than axial stiffness. Is this the reason for the potential for ill-conditioning? Or is there something I'm missing?

pauloz1890
  • 1,217
  • 9
  • 11

2 Answers2

2

You have the right idea, but it doesn't really make sense to say "the ratio of bending to axial stiffness depends on $L^2$" because $L$ isn't dimensionless. If you measured the length of same beam in mm and km you would reach two very different conclusions about the ratio.

A better version of the argument is to say that the axial stiffness is of order to $A/L$ and the bending stiffness $I/L^3$. For a typical beam cross section, $I$ is of order $At^2$ where $t$ is the "thickness" or "diameter" of the beam. So the ratio of axial to bending stiffness is of the order of $t^2/L^2$ (or $A/L^2$ if you prefer that) which is a dimensionless quantity.

In some special situations this doesn't matter at all. For example, if the axis of the beam is aligned with one of the coordinate axes of the model, the axial and bending stiffness terms are not cross-coupled in any way, and there is nothing to degrade the answers.

In practice, even with FE software that uses a simple-minded equation solver working in 64-bit arithmetic (double precision, 16 significant figures), there shouldn't be any serious problems with a ratio of the order of $10^6$, but if you exceed $10^8$ you are starting to live dangerously and at $10^{10}$, all bets are off unless you have some solid evidence that the results are meaningful.

If you have structural analysis FE code that works entirely in 32-bit arithmetic, my advice would be "don't use it for anything except toy problems" - but that advice may be over-pessimistic for some non-structural applications of FE.

alephzero
  • 12,578
  • 1
  • 16
  • 28
0

Yes, you're right. As you imply, in a long, slender beam the bending stiffnesses will be orders of magnitude less than the axial stiffness and mixing numbers that are orders of magnitude different in any numerical calculation can cause problems.

rdt2
  • 314
  • 2
  • 4