Python脚本在Abaqus odb文件中批量提取节点坐标与变形量

摘要

本文介绍如何使用Python脚本二次开发来批量提取ABAQUS输出数据库（ODB）文件中指定Step下的Set节点集变形量。通过详细的步骤说明、代码示例和图片展示，您将学会如何使用该脚本，自动化输出CSV文件包含（Node Label；Step Name、Increment、Step Time，U1，U2）。

如果还需要按Increment提取每个增量下的变形后的节点坐标的话，在提取变形量的基础上，与初始坐标进行简单的计算就可以求得坐标。（备注：该代码只提取了x，y方向的变形量）

1. 问题描述

在工程仿真和分析领域，提取ABAQUS输出数据库（ODB）文件中的节点集变形量是一项常见任务。然而，手动提取这些数据是一项繁琐且容易出错的工作。因此，需要一种自动化的方法来批量提取指定步骤下按节点集组织的变形量数据。

2. 实例展示

假设我们有一个名为`example.odb`的ODB文件，其中包含名为`Step-x`的步骤和名为`Set-x`的节点集。运行以上代码后，脚本会自动将该步骤下节点集的变形量提取出来，并保存为`NodalDisplacement.csv`文件。（图片展示的是名为`Step-4`的步骤和名为`SET-NODE`的节点集）

利用Python脚本批量提取Abaqus的odb文件中Set节点集初始节点坐标、指定Step下的变形量、变形后节点坐标的图1

3. 使用方法

3.1运行前注意事项

1. 在脚本中，设置待处理的ODB文件放置目录。

2. 在脚本中，设置要提取的步骤名称`step_name`和节点集名称`set_name`。

3. 在脚本中，设置一下保存的文件名称。（提取的数据将以CSV格式保存在待处理的ODB文件放置目录下）

3.2运行测试代码

打开Abaqus，输入测试代码，按部就班完成测试，方便后期查询问题。

1.识别路径下odb文件是否存在。

from odbAccess import openOdb
odb_path = r"E:\Abaqus work directory\rubber1.odb" #输入读取路径下的odb文件#
odb = openOdb(path=odb_path)
if odb is not None:
   print("odb found:", odb_path)
else:
   print("odb not found:", odb_path)

2.识别odb文件中关注的Step-4分析步否存在。

from odbAccess import openOdb
odb_path = r"E:\Abaqus work directory\rubber1.odb"
odb = openOdb(path=odb_path)


if odb is not None:
    print("Found:", odb_path)
    step_name = 'Step-4'
    try:
        step = odb.steps[step_name]
        print("Step found:", step_name)
    except KeyError:
        print("Step not found:", step_name)
else:
    print("Not found:", odb_path)

3.识别odb文件中关注的Set组分析步否存在。

from odbAccess import *
odb_path = r"E:\Abaqus work directory\rubber1.odb"
my_odb = openOdb(path=odb_path)


for node_set_name in my_odb.rootAssembly.nodeSets.keys():
    node_set = my_odb.rootAssembly.nodeSets[node_set_name]
    print("Node Set: {}".format(node_set_name))