python - 在 python 脚本 Abaqus 中，STRESS 数组为空

Question

我想在每个节点上提取模型顶部表面的应力，但无法正常完成。当我使用这个脚本时：

odb = visualization.openOdb('My.odb')
frame=odb.steps['AStep'].frames[-1]

dispNode = odb.rootAssembly.nodeSets['UPPER']

STRESS= frame.fieldOutputs['S'].getSubset(region=dispNode).values
COORD= frame.fieldOutputs['COORD'].getSubset(region=dispNode).values

print(STRESS)
print(COORD[1].data)

STRESS 返回一个空数组。

我怎样才能编辑我的脚本有压力和相应的坐标？

Answer 1

如果您只计算积分点上的压力值，您的代码将无法工作。节点上根本没有值，因此如果您在节点上请求值，您将得到一个空数组。

这是它应该如何工作的：

• 将您的集成点结果外推到节点
• 平均您的 ElementNodal 值。这就是它的工作原理：https ://stackoverflow.com/a/43175485/4045774
• 提取节点坐标（变形或未变形）
• 从您的点集中获取节点标签
• 使用您的点集中的节点标签找到相应的唯一节点值https://docs.scipy.org/doc/numpy/reference/generated/numpy.in1d.html

如果您需要一个小的示例代码，请随时询问。