在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

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