2

I'm trying to implementing a Finite Element Analysis algorithm. I solve K u = f to get the displacement u, and then calculate strain with u, then calculate the stress. Finally, I use the stress to calculate the Von Mises Stress, and visualize this. From the result I find the strain is not continuous between tetrahedrons.

I use 10 nodes tetrahedron as the element, so the displacement is a second-order polynomial in every element. The displacement should be enforced to be continuous between tetrahedrons. And the strain, which is the first order derivatives of the displacements should be continuous inside every tetrahedron. But I'm not sure: is this true across the interface between tetrahedrons?

LorenzoDonati4Ukraine-OnStrike
  • 6,875
  • 4
  • 37
  • 56
Yue Xie
  • 23
  • 3
  • one bit of advice, for output visualization you might treat your tets as 4-node. Extrapolpate stresses to the corners and do a linear interpolation. Even if you think that an unsatisfactory approach it will still be a worthwhile exercise to implement. – agentp Apr 25 '13 at 20:25

2 Answers2

2

Only the components of strain tangent to the adjoining face are guaranteed continuous. This follows from the displacement continuity, when you take derivatives in the direction of the interface they are the same.

Commercial FEM programs typically do some post process averaging to make the other components look continuous. Note the strain components normal to an element boundary are only expected to be continuous if the underlying constitutive model is continuous, so such averaging is not always appropriate.

LorenzoDonati4Ukraine-OnStrike
  • 6,875
  • 4
  • 37
  • 56
agentp
  • 6,849
  • 2
  • 19
  • 37
1

You should not compute the stress and strain at the nodes but inside the elements. You can choose for example 4 Gauss points and compute the values there. You then have to think about a scheme on how to get the values computed at the Gauss points onto the tet nodes.

There is a Mathematica application example which illustrates this. Unfortunately the web page is no longer available, but the notebooks are here. You'll find the example in the application example section under Finite Element Method, Structural Mechanics 3D (in the old HelpBrowser). If you have difficulties I could convert it to PDF and send it you.

  • Thank you very much for reply. But I can't find the example in the url you give. This problem is driving me crazy, would you please to send me the pdf? My email is xie.full@gmail.com, do you mind if I send you some picture to illustrate my problem. – Yue Xie Apr 23 '13 at 14:13