0

我是编程新手，我需要编写一个代码来创建一个 odb，以可视化单个步骤的六个字段输出的最大值。对于每个单个节点，我想代表每个字段输出的最大值。现在我创建了一个导致此错误的代码：

myFieldOutput.addData(position=INTEGRATION_POINT, instance=instance1, data=maxStress)

来自序列的 omu_PrimArray 失败 - 维度。

这是代码：

from abaqus import *
from abaqusConstants import *
import visualization
myViewport = session.Viewport(name='Max_Stess',
origin=(10, 10), width=150, height=100)
# Open the tutorial output database.
myOdb = visualization.openOdb(path='PROVA_04_10.odb', readOnly=False)
# Associate the output database with the viewport.
myViewport.setValues(displayedObject=myOdb)
instance1 = myOdb.rootAssembly.instances['MODIFICA_FUNZIONANTE']
# Create variables that refer to the first steps.
firstStep = myOdb.steps['sigma_equivalenti']

frame = firstStep.frames[-1]
#creo nuovo step

# lettura degli stress dallo step sigma

sigma_eq_1 = frame.fieldOutputs['sigma_eq_1']
sigma_eq_2 = frame.fieldOutputs['sigma_eq_2']
sigma_eq_3 = frame.fieldOutputs['sigma_eq_3']
sigma_eq_4 = frame.fieldOutputs['sigma_eq_4']
sigma_eq_5 = frame.fieldOutputs['sigma_eq_5']
sigma_eq_6 = frame.fieldOutputs['sigma_eq_6']

#calcolo del vettore degli sforzi massimi

maxStress=[]
for i in range(len(sigma_eq_1.values)):
   v1=sigma_eq_1.values[i]
   v2=sigma_eq_2.values[i]
   v3=sigma_eq_3.values[i]
   v4=sigma_eq_4.values[i]
   v5=sigma_eq_5.values[i]
   v6=sigma_eq_6.values[i]
   maxStress.append((max(v1,v2,v3,v4,v5,v6),)) #max(valori)


myFieldOutput = frame.FieldOutput(name='tensioni_max',description='calcolo delle tensioni massime', type=SCALAR)
myFieldOutput.addData(position=INTEGRATION_POINT, instance=instance1, data=maxStress)

#visualizzazione della variabile delle sigma massima sul programma
myViewport.odbDisplay.setPrimaryVariable(field=myFieldOutput,outputPosition=INTEGRATION_POINT)
myViewport.odbDisplay.display.setValues(plotState=(CONTOURS_ON_DEF,))
4

1 回答 1

1

您可以使用field参数将字段输出数据添加为FieldOutput object. 但是，您应该从现有的字段输出创建字段数据。
在您的情况下，您可以执行以下操作：

  1. 访问字段输出（不是字段输出数据），就像您做的那样。
sigma_eq_1 = frame.fieldOutputs['sigma_eq_1']
sigma_eq_2 = frame.fieldOutputs['sigma_eq_2']
sigma_eq_3 = frame.fieldOutputs['sigma_eq_3']
sigma_eq_4 = frame.fieldOutputs['sigma_eq_4']
sigma_eq_5 = frame.fieldOutputs['sigma_eq_5']
sigma_eq_6 = frame.fieldOutputs['sigma_eq_6']
  1. 然后，操作这些现场输出数据。您可以对这些数据进行任何可行的数学运算，例如。加减等 在这种情况下，为了找到最大值，我们有一个名为maxEnvelope.
maxS = maxEnvelope([sigma_eq_1, sigma_eq_2, sigma_eq_3, sigma_eq_4, sigma_eq_5, sigma_eq_6])

请注意，此命令将列表作为参数。并且该命令在原始数据可用的相同位置计算新的字段数据。在这种情况下，它在集成点可用，因此将在集成点写入。
3. 现在，您可以使用fieldarugment 添加字段输出数据。

myFieldOutput = frame.FieldOutput(name='tensioni_max',description='calcolo delle tensioni massime', type=SCALAR)
myFieldOutput.addData(field=maxS[0])

PS：由于您已将可视化模块导入为import visualization，因此您应该使用maxEnvelope方法 asvisualization.maxEnvelope(...)

于 2021-10-16T05:29:07.550 回答