在Abaqus ODB中处理张量数据

时间:2018-11-17 16:45:59

标签: python tensor abaqus

我想对odb的每个模式(模态分析的框架)的所有应变能密度(ELSE)值求和,并将其保存在新框架中。

我的代码在位移'U'和type = VECTOR上工作正常 但是,当我将其更改为“ ELSE”并键入= TENSOR_3D_FULL时,出现以下错误:

Traceback (most recent call last):    File "..\Python\modaldata.py", line 62, in <module> 
    newField.addData(field=strechtMode1)  OdbError: Fields are not compatible. They are either associated with different structural models or have incompatible attributes.

我是否必须使用另一种方法来处理张量数据?我找到了一些具有position = ELEMENT_NODAL的东西,但通过获取某些节点集的某些标量应力值而采用了完全不同的方法。我希望将每个元素的每个集成点的整个ELSE张量作为一个新框架。

这是我的代码:

#-ODB OEFFNEN------------------------------------------------------------------------#    
odb = openOdb(odbName) 
#-DIE MODEN EINZELN ABSPEICHERN UM SIE VERRECHNEN ZU KOENNEN-------------------------#
nmoden=len(odb.steps['Modalanalyse'].frames)

#------------------------------------------------------------------------------------#
#-FIELD VARIABLE ZIEHEN UND VERRECHNEN-----------------------------------------------#
#------------------------------------------------------------------------------------#

field='U' # ''ELSE' # U oder E oder ELSE oder whatever...
type=VECTOR #TENSOR_3D_FULL  # muss an das field angepasst werden. 

Modes_field=[]
for i in range(nmoden):
    Modes_field.insert(i,odb.steps['Modalanalyse'].frames[i].fieldOutputs[field])  
#-FIELD VARIABLE ALLER MODEN UEBERLAGERN---------------------------------------------#
sumAllModes_field = Modes_field[0]
for i in range(nmoden):
    sumAllModes_field = sumAllModes_field + Modes_field[i] 
#-NEUE FIELD VARIABLE MODE ABSPEICHERN-----------------------------------------------#
newFrame = odb.steps['Modalanalyse'].Frame(incrementNumber=0, frameValue=0.0)
newField = newFrame.FieldOutput(name=field, description='sum of ELSE of all modes', type=type)
newField.addData(field=sumAllModes_field)
#-ODB ABSPEICHERN--------------------------------------------------------------------#
odb.save()

1 个答案:

答案 0 :(得分:1)

ELSE是总能量。
因此,其字段输出类型为SCALAR。 这就是为什么出现错误的原因,您试图用SCALAR字段值填充TENSOR_3D_FULL字段。

我强烈建议您使用textRepr库中的prettyPrint来访问字段并更好地了解其属性。

from textRepr import prettyPrint as pp  

pp(odb.steps['Modalanalyse'].frames[-1].fieldOutputs['ELSE'])
pp(sumAllModes_field)
pp(newField)