import csv
from odbAccess import openOdb
from abaqusConstants import *
#Extract Mises
def getMISES(output_path_mises, odb_path):
odb = openOdb(odb_path)
mises_data = odb.steps['Belastung'].frames[-1].fieldOutputs['S']
#max mises is always in this instance
instance = odb.rootAssembly.instances['PRUEFKOERPER_4-1']
node_coords = {node.label: node.coordinates for node in instance.nodes}
with open(output_path_mises, 'w') as csv_file:
csv_writer = csv.writer(csv_file)
csv_writer.writerow(['NodeLabel', 'x', 'y', 'z', 'MISES'])
for block in mises_data.bulkDataBlocks:
for node_label, mises_value in zip(block.nodeLabels, block.data):
if mises_value > 300.0: #Filter
coords = node_coords.get(node_label, [None, None, None])
csv_writer.writerow([node_label, coords[0], coords[1], coords[2], '{:.4f}'.format(float(mises_value))])
odb.close()
您的方法很少。
mises_data = odb.steps['Belastung'].frames[-1].fieldOutputs['S']
stress_data = odb.steps['Belastung'].frames[-1].fieldOutputs['S']
mises_data = stress_data.getScalarField(invariant=MISES)
access nodal信息以应对压力:由于压力是在集成点写的,因此使用
mises_data.bulkDataBlocks[0].nodeLabels