Why do the shape of finite elements matter?

Question

I have used FEA for a couple of years now, but using it and using it correctly are two different things, safety factor is not the solution to everything. I have the feeling I won't be using it right unless I have a clear answer to that question:

I am aware elements must be close to their ideal shape (based on the Jacobian only?) in order to get accurate results.. But why? Since I understand it comes from a coordinate transform, unless two vectors of the element become colinear shouldn't the results be accurate no matter its shape?

A step-by-step answer based on an illustrated example (arbitrary stress distribution) would be ideal, especially given that it is a relatively common question (but never well answered from what I have seen).

P.S: This is a "repost" of the Physics stackexchange, I had initially posted it in the wrong one. — Mister Mystère, May 04 '14 at 14:27

score 11 · Accepted Answer · answered May 04 '14 at 15:19

11

A good introduction to how issues of element shape influence quality and ease of solution, with pictures, is Jon Shewchuk's "What Is a Good Linear Finite Element? Interpolation, Conditioning, Anisotropy, and Quality Measures"

http://www.cs.berkeley.edu/~jrs/papers/elemj.pdf

answered May 04 '14 at 15:19

Peter Brune

1,675
12
20

Thank you for the answer, this one looks extremely interesting. I will have a look at it shortly; in the meantime I am still open to "short versions" as a summary and a complement. – Mister Mystère May 05 '14 at 12:22

score 3 · Answer 2 · edited Jul 10 '14 at 15:55

You're right, the Jacobian would be invertible unless you have colinear/coplanar elements, i.e., degenerate elements. But this is the case when you have exact arithmetic. In real FEM programs you use floating point arithmetic and then, you can have Jacobians that are nearly zero. Probably, you can still invert this $2\times2$ or $3\times3$ matrix, but the error will accumulate in your model as a whole. Finally... getting an ill-conditioned matrix.

Why do the shape of finite elements matter?

2 Answers2