# Extract max displacement from .odb odb = openOdb(path=jobName + '.odb') lastFrame = odb.steps['ApplyLoad'].frames[-1] displacement = lastFrame.fieldOutputs['U'] max_disp = max([value.dataDouble for value in displacement.values]) output_file.write(f'E, max_disp\n') odb.close()

# batch_postprocess.py import os from odbAccess import * odb_folder = './simulation_results/' output_lines = []

# Get last step, last frame step = odb.steps.values()[-1] frame = step.frames[-1] # Reaction force (RF) rf_field = frame.fieldOutputs['RF'] total_RF = sum([value.data[1] for value in rf_field.values]) # Y-direction sum output_lines.append(f"filename, Total RF = total_RF:.2f") odb.close() with open('reaction_forces.csv', 'w') as f: f.write("ODB_File, Total_RF_Y\n") f.write("\n".join(output_lines))

for E in modulus_values: modelName = f'Model_E_int(E)' mdb.Model(name=modelName) # ... (beam creation similar to Example 1, but with E variable) # Omitted for brevity – reuse Example 1 with dynamic model naming