By coading u can do it

algorithm which was the base of my proper algorithm for python abaqus features extraction (I attach part of it). You should take into acccount that stresses are determined at integration points, not at the nodes. For viewing them at nodes, abaqus uses an extrapolation algorithm.

Therefore, my advice is to extract it at integration points or at the element centroid (this is the way used on my algorithm).

I hope you'll find it useful, it took some hard work. You'll only need to add the file writing on it.

# Get ABAQUS interface

from abaqus import *

from abaqusConstants import *

# Open the odb

myOdb = session.openOdb(name=odbPath)

session.viewports['Viewport: 1'].setValues(displayedObject=myOdb)

# Get the frame repository for the step, find number of frames (starts at frame 0)

frames = myOdb.steps[stepName].frames

numFrames = len(frames)

# Isolate the instance, get the number of nodes and elements

myInstance = myOdb.rootAssembly.instances[instanceName]

numNodes = len(myInstance.nodes)

numElements = len(myInstance.elements)

for el in range(0,numElements):

   # Isolate current and previous element's stress field

Stress=myOdb.steps[stepName].frames[fr].fieldOutputs['S'].getSubset(region=myInstance.elements[el],position=CENTROID, elementType='C3D8H').values

   sz = len(Stress)

   for ip in range(0,sz):

      Sxx = Stress[ip].data[0]

      Syy = Stress[ip].data[1]

      Szz = Stress[ip].data[2]

      Sxy = Stress[ip].data[3]

      Sxz = Stress[ip].data[4]

      Syz = Stress[ip].data[5]

myOdb.close()

By coading u can do it

algorithm which was the base of my proper algorithm for python abaqus features extraction (I attach part of it). You should take into acccount that stresses are determined at integration points, not at the nodes. For viewing them at nodes, abaqus uses an extrapolation algorithm.

Therefore, my advice is to extract it at integration points or at the element centroid (this is the way used on my algorithm).

I hope you'll find it useful, it took some hard work. You'll only need to add the file writing on it.

# Get ABAQUS interface

from abaqus import *

from abaqusConstants import *

# Open the odb

myOdb = session.openOdb(name=odbPath)

session.viewports['Viewport: 1'].setValues(displayedObject=myOdb)

# Get the frame repository for the step, find number of frames (starts at frame 0)

frames = myOdb.steps[stepName].frames

numFrames = len(frames)

# Isolate the instance, get the number of nodes and elements

myInstance = myOdb.rootAssembly.instances[instanceName]

numNodes = len(myInstance.nodes)

numElements = len(myInstance.elements)

for el in range(0,numElements):

   # Isolate current and previous element's stress field

Stress=myOdb.steps[stepName].frames[fr].fieldOutputs['S'].getSubset(region=myInstance.elements[el],position=CENTROID, elementType='C3D8H').values

   sz = len(Stress)

   for ip in range(0,sz):

      Sxx = Stress[ip].data[0]

      Syy = Stress[ip].data[1]

      Szz = Stress[ip].data[2]

      Sxy = Stress[ip].data[3]

      Sxz = Stress[ip].data[4]

      Syz = Stress[ip].data[5]

myOdb.close()

Thank you so much Manigandan Sekar.

I agree with Manigandan answer

Thank you Dr. Chiad