Hi all
I am working on composite material using Abaqus and I want to know if there is a possibility to extract the material orientation for each element using Python script ??
if you have some ideas please tell me
Thank you
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})
- Badges
- Science topic