I am using Fenics Mechanics to simulate the stress in a solid body. The example code provided by Fenics Mechanics only outputs displacement. Anyone could suggest how I could output stress? The two lines used for output are as follows:
problem = fm.SolidMechanicsProblem(config)
solver = fm.SolidMechanicsSolver(problem, fname_disp=“results/displacement.pvd”)
Many thanks!