Python script for Tsai-Wu usage factor

Hi,

I have a Python script used to iterate through all elements to find the maximum von Mises stress. I would like to have a similar script to loop through all elements to find the maximum Tsai-Wu usage factor, checking top mid and bottom of each ply in the laminate.

The script should report:

* max Tsai-Wu usage factor

* element number

* Ply number

* Ply location (top, mid or bottom)

I am not able to figure out how to do this from the existing Abaqus documentation. Anyone who have a similar script or can give some suggestions?

Thanks!

Best regards,

Geir Olav