Stiffness of tri and tet elements
Stiffness of tri and tet elements
(OP)
I took an FEA methods course one time that went over element types. The course said that quad and hex element meshes were the best provided certain aspect ratios were not exceeded, and that tri and tet element meshes were "stiff" and tended to overestimate stresses. This was in the 80's. Have the software programs progressed to overcome this one-time weakness of tri and tet elements? The reason I ask is that I recently used an FEA code that meshes only in tet elements. The code, however, was a CFD code. When solving CFD problems is there an inherent weakness with tet elements like there once was with stress problems? TIA for your thoughts.
Tunalover





RE: Stiffness of tri and tet elements
Tunalover
RE: Stiffness of tri and tet elements
http://www.colorado.edu/engineering/CAS/courses.d/...
Cheers
Greg Locock
New here? Try reading these, they might help FAQ731-376: Eng-Tips.com Forum Policies http://eng-tips.com/market.cfm?
RE: Stiffness of tri and tet elements
RE: Stiffness of tri and tet elements
ANSYS CFD consists of Fluent and Icepak, both of which are Finite Difference Method programs where element types and aspect ratios are not as important as with FEA software. Am I just refining your statement?
Tunalover
RE: Stiffness of tri and tet elements
Thanks for the link. That's a pretty thorough and lengthy paper so I just read the conclusion. The conclusion said basically "brick elements are better than tetrahedron elements" but didn't really explain why. But in the body of the paper it said that tetrahedral elements are stiffer (as I was taught) but in the brief scan of the paper I didn't see anything in the way of validation of that statement.
Bruce Jackson
aka Tunalover
Tunalover
RE: Stiffness of tri and tet elements
Only looking at the case of linear elements, we see that the tetrahedron element has only one integration point in the entire element. I have seen some people even call these "constant strain" elements because there cannot be a strain/stress gradient throughout the element, if it only has one integration point.
Comparing this to the linear hexahedral (8 node element), there are two integration points in each orthogonal direction for the master element. Clearly, there can now be a strain/stress gradient in the element.
The way stress is computed in FEA requires that integration point value is mapped to the nodes on the elements, and then we acquire our nodal stresses.
The answer I think you are looking for is that tet elements are usually not considered sufficient for most analyses that require a reasonable level of accuracy for calculating stress. That being said, the computational cost is different.
I try to use tet elements in places where A) nice hex mapped meshes are not easy to create, and B) we don't care about high strains, etc.
That being said, what I just said applies only to structural FEA. I have very little experience with CFD, and I realize that the fundamental differential equations are different, but I would be willing to extend this analogy if you want a simple answer.
RE: Stiffness of tri and tet elements
One thing in common in finite element is data is computed at the integration point and result between each element is computed by taking average from integration point of two elements.
Since hex element has more faces than tri element, a single hex element has more connections with other elements. Therefore hex element gives smoother result than tri. Also since tri provides a more discrete result, the structural result will look "stiffer". This is the stiffing effect when less elements are used in structural analysis
P.S. This conclusion is based on my education knowledge and experience. I have no mathematical prove, so please feel free to correct me
RE: Stiffness of tri and tet elements
For my employer years ago I purchased Fluent and Icepak and attended the training. I also studied CFD in grad school and wrote my own programs. At both places it was clear that the codes were Finite Difference Method, not FEA. Unless they rewrote the codes when ANSYS purchased Fluent and Icepak (which is highly unlikely) these two programs are based on the Finite Difference Method. As for your reasoning as to why tet elements are stiffer than hex elements, I'm going to stick with csk62's explanation as it does agree with my understanding of FEA with regard to Gauss integration points and how they are used to compute stresses.
Tunalover
RE: Stiffness of tri and tet elements
TET4s suffer a similar limitation.
TRIA6 and TET10 don't.
another day in paradise, or is paradise one day closer ?