I used linear_elasticity.py as a starting point for simulating stress in thin silicon layers. My manual meshing worked, and the displacements seem to be calculated properly. However, how do I get the stress in each element? I need it for determining where the layer will break.
The original node positions of the mesh are totally relaxated (no tension anywhere).