Which Q1P0 element? 2D/3D small or finite deformation?
Some elements do print all the results. For example, sld2d2.f (the small deformation mixed element does). All elements compute the stresses, however some then average to print a single point. Since the element you are using does not print all it will be necessary for you to modify the code yourself.
To do this:
Go to the directory where the element is located (e.g. /elements/solid3d contains the 3-d elements)
Edit the element routine (sld3d2.f for small fld3d2.f for finite - for 3d)
Locate the region where 'isw = 4' (may be a go to or an isw.eq.4 transfer)
Follow throught the computation of the stress (stresses are named sig... (maybe sigl).
Print values of the stress for eacy quadrature point (inside a do l = 1,lint loop).
You may need to add a format or you may be able to use the existing one.
