python脚本Abaqus中的应力数组为空

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)

应力返回空数组。在

如何编辑我的脚本，使其具有应力及其对应的坐标？？在

Answer 1

如果只计算积分点上的应力值，代码就不能工作。节点上根本没有值，因此如果在节点上请求值，则会得到一个空数组。在

这就是它的工作原理：

• 将结果外推到积分节点
• 求元素节点值的平均值。工作原理如下：https://stackoverflow.com/a/43175485/4045774
• 提取节点坐标（变形或未变形）
• 从点集中获取节点标签
• 使用点集中的节点标签，找到相应的唯一节点值https://docs.scipy.org/doc/numpy/reference/generated/numpy.in1d.html

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