Dear Abdelhak Khechai,

Please check it out through the use of manual in Abaqus. 

Dear Walled

Thank you for your reply

Yes, I tried to find it in the manual, but I just found how I can extract the stress, displacements for each node.... no material orientation for each element 

Hi Abdelhak,

You can access the "localCoordSystem" variable to extract the material orientation for each element (refer to the bottom of this thread: http://www.eng-tips.com/viewthread.cfm?qid=207919)

So, something like:

fileName='pathToYourOutputDatabase/yourODBfile.odb'

outputDatabase = session.openOdb(name= fileName )

myFrame=outputDatabase.steps['step_name'].frame[frame_number]

field=myFrame.fieldOutputs['U']

values=field.values

for iValue in values:

    print "Element=",iValue.elementLabel

    print "IP=",iValue.integrationPoint

    # This is the part you are looking for:

    print "Local system=", iValue.localCoordSystem

You may also have "LOCALDIR#" fields that you can access. You can check in CAE, or in ABAQUS Python:

fieldNames = myFrame.fieldOutputs.keys()

print fieldNames

Ken

P.S., I recommend that you do a simple test with a single element first to make sure that you are getting the correct values :)

Dear Kenneth I Aycock

I will try to apply your script. If there is positive results I will tell you.

Thank you very much

Hi Kenneth I Aycock and Abdelhak Khechai ,

I'm wondering have you solved this problem to export the deformed material orientation of each element.

I've tried the method above but it doesn't work even though I assign the material orientation for different element...I just get nones for all the integrationPoint , localCoordSystem and elementLabel.

I provide my python command as below and a figure and an .inp as attached.

"""Write Material Orientation to the file "output_fieldvalues.txt" """

fileName='C:/temp/PythonScripting/Job-1.odb'

outputDatabase = session.openOdb(name= fileName )

myFrame=outputDatabase.steps['Step-1'].frames[17]

field=myFrame.fieldOutputs['U']

values=field.values

OutputFile = open('output_fieldvalues.txt','w')

for iValue in values:

OutputFile.write("Element={}, IP={}, Local coordinate system={}.\n".format(iValue.elementLabel, iValue.integrationPoint, iValue.localCoordSystem))

OutputFile.write('Done')

OutputFile.close()

Hope someone can help.

The output I get from typing "print values[1]" is the following

({'baseElementType': '', 'conjugateData': None, 'conjugateDataDouble': 'unknown', 'data': array([0.109024919569492, -0.0102628562599421, -6.86978660332312e-16], 'f'), 'dataDouble': 'unknown', 'elementLabel': None, 'face': None, 'instance': 'OdbInstance object', 'integrationPoint': None, 'inv3': None, 'localCoordSystem': None, 'localCoordSystemDouble': 'unknown', 'magnitude': 0.109506890177727, 'maxInPlanePrincipal': None, 'maxPrincipal': None, 'midPrincipal': None, 'minInPlanePrincipal': None, 'minPrincipal': None, 'mises': None, 'nodeLabel': 2, 'outOfPlanePrincipal': None, 'position': NODAL, 'precision': SINGLE_PRECISION, 'press': None, 'sectionPoint': None, 'tresca': None, 'type': VECTOR})