大家好

調用UEL(一個模擬粘結滑移的用戶單元）的abaqus模型，如何用Python后處理UEL。
比如如何查看用戶單元（4節點平面單元）的積分點的mises應力？

我嘗試訪問了輸出數據庫（odb）,編寫如下script:

#!/usr/bin/python

#-*-coding: UTF-8-*-

from abaqus import*

from odbAccess import*

odb=openOdb(path='verifiedmodel.odb')

#單元集合P1就是用戶單元集合

topCenter=odb.rootAssembly.instances['NS-1'].elementSets['P1']

stressField = odb.steps['Step-1'].frames[-1].fieldOutputs['S']

field = stressField.getSubset(region=topCenter, 

    position=INTEGRATION_POINT, elementType = 'U1')

fieldValues = field.values

for v in fieldValues: 

    print 'Element label = ', v.elementLabel,     

    if v.integrationPoint:

        print 'Integration Point = ', v.integrationPoint

    else:

        print 

    for component in v.data:

        print '%-10.5f' % component

    print

運行沒有結果。

但是把單元集合換成C2（混凝土單元）或者是S1（鋼筋單元），運行是有結果的。

另外，我試著僅運行下列腳本程序：

stressField = odb.steps['Step-1'].frames[-1].fieldOutputs['S']

print stressField 

結果只有T3D2(鋼筋的類型)，C3D8R（混凝土的類型），就是沒有U1（用戶單元類型）。這是為什么呢？

整個模型是通過運行INP文件，并且調用用戶子程序運行出來的，verifiedmodel.odb是該模型的輸出數據庫，是不是說這個ODB不包含用戶單元的信息。

眾所周知，abaqus的CAE界面的可視化窗口是不會顯示用戶單元模型的。

在書中找了資料，莊茁的有本書上說，STA文件會有顯示用戶單元信息，我查看了，并沒有。

在網上找了一個網名為（窗臺上的叔本華），講了要通過編寫python腳本來查看用戶單元運行的信息，上面所列的程序就是我的嘗試，并沒有得到結果。

所以想請教大家，如何得到運行后用戶單元的信息呢？用python編寫腳本可行嗎，是程序出了問題嗎，還是ODB根本就不包含用戶單元的信息呢？