Abaqus ODB中张量数据的处理

Question

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

我的代码可以很好地处理位移'U‘和type=VECTOR，但是当我将它改为'ELSE’和type=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()

Answer 1

ELSE是一种总能量。

因此，它的字段输出类型是标量。这就是你得到错误的原因，你试图用标量字段值填充TENSOR_3D_FULL字段。

我强烈推荐使用textRepr库中的prettyPrint，以便访问字段并更好地了解其属性。

from textRepr import prettyPrint as pp  

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