I use Python to generate my models and read the .odb files after the simulation. I am particularly interested in the coordinates of a set of nodes that I call T_N. It was working fine in Abaqus 6.12, but when I use it 6.13, it does not get all the nodes of a given subset (using getSubset). I believe it is not a problem during the simulation since the results in Abaqus are the same. It also works fine if I try to get all the nodes at the same time instead of one subset.

I have rewritten the function the read the database following the documentation and I am using this:

def myNodesReader(odb):

step = odb.steps['Hit']

timeFrameInitial =step.frames[0]

try:

displacementInitial = timeFrameInitial.fieldOutputs['COORD']

except:

err ="The analysis must include a field output request \\ for variable COORD."

print err

sys.exit(1)

try:

targetNodesInitial = odb.rootAssembly.nodeSets['T_N']

targetDisplacementInitial = displacementInitial.getSubset(region=targetNodesInitial)

except:

err ="T_N must be defined"

print err

sys.exit(1)

outFile = open('InitialNodes_TN.txt', 'w')

format = '%6i, %14.7e, %14.7e, %14.7e\\n'

for value in targetDisplacementInitial.values: outFile.write(format % ((value.nodeLabel,) + tuple(value.data)))

outFile.close()

These are the results I get from Abaqus 6.13:

This is from Abaqus 6.12:

Please let me know if my problem is not clear.

Thank you very much in advance!