Output stress at Gauss points in umacr

[K.Li]

Dear FEAP users:

I am trying to modify the FEAP umacr9 routine for ParaView interface to output stress at Gauss points, instead of on the nodes. So could anyone please tell me where are the stresses at Gauss points located? Is it mr() or hr()?

Thanks a lot!

[Prof. R.L. Taylor]

Gauss point stresses are not stored.  They are computed and output in each element routine under isw=4.  Thus you need to recompute them in order to output to your file.  Hopefully, you are only using one element type!  Otherwise it is a tedious task.

[K.Li]

Thank you for the information, Prof Taylor!