如何用Python脚本将Abaqus场输出中的零值替换为极小非零值？

问题背景 #

需要将Abaqus场输出中U2（Y方向位移）的所有零值替换为极小非零值（1e-9），避免后续运算出现除以零问题。但运行脚本后，新生成的U2_no_zeros场与原U2场数值完全一致，零值未被替换。

用户原脚本如下：

from abaqusConstants import *
from odbAccess import *

# ***********************************************
odbPath="path_to_odb_file"
stepName="Step-1"
frameNumber=-1 #last frame in the stepName
sourceOutputFieldName='U' #displacement field
newOutputFieldName='U2_no_zeros'
# ************************************************
odb=session.openOdb(name=odbPath,readOnly=FALSE)
step=odb.steps[stepName]
frame=step.frames[frameNumber]
AllInstances=(odb.rootAssembly.instances.keys())
MyInstance=(AllInstances[-1])
instance1=odb.rootAssembly.instances[MyInstance]
sourceField=frame.fieldOutputs[sourceOutputFieldName]
subField=sourceField.getScalarField(componentLabel="U2")
Values=subField.bulkDataBlocks[0].data
NodeLabels=subField.bulkDataBlocks[0].nodeLabels
for value in Values:
  if value==0:
     value=1e-9
newField=frame.FieldOutput(name=newOutputFieldName, type=SCALAR, description="field")
newField.addData(position=NODAL, instance=instance1, labels=NodeLabels, data=Values)
odb.save()
odb.close()

错误原因 #

原脚本中for value in Values的遍历方式仅操作了循环内的临时变量value，并未修改Values数组本身的元素。因为Values是Abaqus返回的numpy数组，直接遍历元素只能读取，无法原地修改原数组内容。

解决方案 #

方式1：索引遍历修改（基础兼容） #

直接通过数组索引遍历，修改原数组的对应元素：

from abaqusConstants import *
from odbAccess import *

# ***********************************************
odbPath="path_to_odb_file"
stepName="Step-1"
frameNumber=-1 #last frame in the stepName
sourceOutputFieldName='U' #displacement field
newOutputFieldName='U2_no_zeros'
# ************************************************
odb=session.openOdb(name=odbPath,readOnly=FALSE)
step=odb.steps[stepName]
frame=step.frames[frameNumber]
AllInstances=(odb.rootAssembly.instances.keys())
MyInstance=(AllInstances[-1])
instance1=odb.rootAssembly.instances[MyInstance]
sourceField=frame.fieldOutputs[sourceOutputFieldName]
subField=sourceField.getScalarField(componentLabel="U2")
Values=subField.bulkDataBlocks[0].data
NodeLabels=subField.bulkDataBlocks[0].nodeLabels

# 修正：通过索引遍历修改原数组元素
for i in range(len(Values)):
    if Values[i] == 0:
        Values[i] = 1e-9

newField=frame.FieldOutput(name=newOutputFieldName, type=SCALAR, description="Y-direction displacement with zero values replaced by 1e-9")
newField.addData(position=NODAL, instance=instance1, labels=NodeLabels, data=Values)
odb.save()
odb.close()

方式2：numpy向量化操作（高效推荐） #

利用numpy的布尔索引批量替换零值，适合大数据量场景，效率更高：

from abaqusConstants import *
from odbAccess import *
import numpy as np

# ***********************************************
odbPath="path_to_odb_file"
stepName="Step-1"
frameNumber=-1 #last frame in the stepName
sourceOutputFieldName='U' #displacement field
newOutputFieldName='U2_no_zeros'
# ************************************************
odb=session.openOdb(name=odbPath,readOnly=FALSE)
step=odb.steps[stepName]
frame=step.frames[frameNumber]
AllInstances=(odb.rootAssembly.instances.keys())
MyInstance=(AllInstances[-1])
instance1=odb.rootAssembly.instances[MyInstance]
sourceField=frame.fieldOutputs[sourceOutputFieldName]
subField=sourceField.getScalarField(componentLabel="U2")
Values=subField.bulkDataBlocks[0].data
NodeLabels=subField.bulkDataBlocks[0].nodeLabels

# 修正：numpy批量替换零值
Values[Values == 0] = 1e-9

newField=frame.FieldOutput(name=newOutputFieldName, type=SCALAR, description="Y-direction displacement with zero values replaced by 1e-9")
newField.addData(position=NODAL, instance=instance1, labels=NodeLabels, data=Values)
odb.save()
odb.close()

额外注意事项 #

• 浮点数精度：若原U2值为接近0的极小浮点数（如1e-16），直接判断==0会遗漏，可改为abs(Values[i]) < 1e-12或np.isclose(Values, 0)来替换接近零的值。
• 多实例场景：原脚本仅处理最后一个实例，若模型包含多个实例，需添加循环遍历所有实例。

内容的提问来源于stack exchange，提问作者FEA-eng