1 回答
TA贡献1810条经验 获得超4个赞
我终于找到了可以帮助很多人的问题的解决方案(至少它帮助了我的同事)。如果您想通过 Abaqus 脚本了解表面的应力,max9111提出的外部节点平均应力值的解决方案是完美且唯一的解决方案。但是,代码应修改如下才能运行,至少在 Abaqus 6.16 上:
第一部分没有变化
Field = session.odbs['ModelName.odb'].steps['StepName'].frames[0].fieldOutputs['S']
Field = Field.getSubset(position = ELEMENT_NODAL)
Values=Field.bulkDataBlocks[0].data
NodeLabels=Field.bulkDataBlocks[0].nodeLabels
NodeLabels_unique, unq_idx = np.unique(NodeLabels, return_inverse=True)
Values_Averaged=np.zeros((NodeLabels_unique.size,Values.shape[1]))
unq_counts = np.bincount(unq_idx)
由于 Abaqus 使用 bulkDataBlocks 存储数据的方式,Values 变量的创建方式发生了变化:
for i in xrange(0,Values.shape[1]):
ValuesTemp = [item[i] for item in Values]
unq_sum = np.bincount(unq_idx, weights=ValuesTemp)
Values_Averaged[:,i] = unq_sum / unq_counts
max_ind=np.unravel_index(np.argmax(Values_Averaged),Values_Averaged.shape)
print("The max stress is at NodeLabel "+str(NodeLabels_unique[max_ind[0]])+ " its value is "+ str(Values_Averaged[max_ind]) +" MPa.")
添加回答
举报