Extracting stresses from Python

RE: Extracting stresses from Python

Stefanos Papavasileiou — Tue, 09 Mar 2021 16:10:31 GMT

Dear Douglas,

The difference highly depends on the size of the generated elements.

The smaller they are the more accurate the node results will be. I would recommend checking directly on stress points if you want high accuracy on results that are purely computed on integration points, e.g. stresses.

RE: Extracting stresses from Python

Douglas Tun — Mon, 01 Mar 2021 23:02:13 GMT

Hi Stefanos,

Thanks for your comments. I tried the stress point comment, but there is a minimal stress point compared to node points.

Is there a huge difference in results?

Thanks,
Doug

RE: Extracting stresses from Python

Stefanos Papavasileiou — Thu, 25 Feb 2021 18:23:27 GMT

Dear Douglas,

Firstly, one difference arises from the fact that you are checking results for stresses that are primarily computed on stresses but via nodes (extrapolated results).

You can run the command:
soilS1 = g_o.getresults(g_o.Phases[-1], g_o.ResultTypes.Soil.SigmaTotal1, 'stress point')

The getsingleresult command works in a different way. You specify a position and the value will get extrapolated from the closest stress points available.

Therefore checking stress results on stress points is the most accurate, but when looking at nodes and on a random location requires an extrapolation process which depending on the mesh density may affect the results.