**Moderator:** bernd

Wait that’s not right. This is a pit of 8x8x4m deep.

Now 4x4x4m:

Now 38% stronger, but still not enough. This pit will collapse when reaching 3m depth or so.

Now 38% stronger, but still not enough. This pit will collapse when reaching 3m depth or so.

yeah this is what I like ... real 3D excavations !

most engineering tools only allow 2D excavation and soil analysis.

there is one ")" to much above at the end ...HarryvL wrote: ↑Tue Mar 19, 2019 9:52 amWith as an obvious simplification:

Code: Select all

`normal=stress[:3] shear=stress[3:] pressure=np.average(normal) sig_mises=np.sqrt(1.5*np.linalg.norm(normal-pressure)**2+3.0*np.linalg.norm(shear)**2))`

Code: Select all

```
normal=stress[:3]
shear=stress[3:]
pressure=np.average(normal)
sig_mises=np.sqrt(1.5*np.linalg.norm(normal-pressure)**2+3.0*np.linalg.norm(shear)**2)
```

like this: https://github.com/berndhahnebach/FreeC ... diff=split ?

I added a unit test lately. It still is ok The von mises test only coúld be run on linux with:

Code: Select all

`./bin/FreeCAD --run-test "femtest.testresult.TestResult.test_stress_von_mises"`

Yes, like that . It looks more Pythonic and should be slightly faster.

BTW: what is best to name the tuple for the 6 stress components. Ist stresstuple fine, or just stress (I don't like this), or stresstensor or whatever ? We gone make some use of it in femconcrete branch too. https://github.com/berndhahnebach/FreeC ... f6e4e5L473

stress_tensor sounds fine. I don't like to name variables after type. If down the line you change the type for some reason you'll have to change variable name as well.