Abaqus的歷史數(shù)據(jù)輸出的問題

Abaqus結(jié)果輸出時(shí)有場輸出（Field Output）和歷史輸出(History Output)，場變量輸出用于描述某個(gè)量隨空間位置的變化，歷史變量用于描述某個(gè)量隨時(shí)間的變化，區(qū)別是場變量輸出大量的單元或節(jié)點(diǎn)上的計(jì)算結(jié)果，寫入odb文件的頻率低，用與生成后出的各個(gè)圖。歷史變量輸出少量單元或節(jié)點(diǎn)上的計(jì)算結(jié)果，寫入Odb的頻率高，用于生成X-Y圖。歷史變量允許單獨(dú)輸出某個(gè)獨(dú)立分量，經(jīng)常用到的就是這個(gè)獨(dú)立分量的輸出，比如輸出結(jié)構(gòu)最大點(diǎn)的位移變化曲線，或者盈利最大點(diǎn)的應(yīng)力變化過程等。

在輸出歷史數(shù)據(jù)時(shí)，大多數(shù)的輸出結(jié)果比如位移和應(yīng)力等，不能實(shí)現(xiàn)整個(gè)模型全部輸出，也就是需要指定set（參考集），需要設(shè)置相關(guān)的set后，結(jié)果輸出才會(huì)有這些數(shù)據(jù)。

定義Set的對(duì)話框如圖1所示，定義Set的類型時(shí)有三個(gè)，按照以往的理解，一般經(jīng)常用到的應(yīng)該是Node，也即節(jié)點(diǎn)集，想要輸出這個(gè)節(jié)點(diǎn)在求解過程中的Mises應(yīng)力變化過程，選擇某個(gè)節(jié)點(diǎn)，通過Step模塊的History Output Requests，如圖2所示。

圖1

圖2

選擇剛才設(shè)置的Set，選擇需要輸出的結(jié)果項(xiàng)，此處為Mises。之后計(jì)算，計(jì)算后提取結(jié)果時(shí)，發(fā)現(xiàn)并沒有輸出該節(jié)點(diǎn)的Mises歷史數(shù)據(jù)。

筆者很快想到了，有可能是不能輸出節(jié)點(diǎn)的應(yīng)力，因?yàn)樯婕暗焦?jié)點(diǎn)解和單元積分點(diǎn)解的問題，這個(gè)之前已經(jīng)討論過。基于此考慮，那么輸出該節(jié)點(diǎn)的位移場肯定沒問題，因?yàn)楣?jié)點(diǎn)位移解是最先得到的節(jié)點(diǎn)解。

將Mises換成UT，如圖3所示。

同時(shí)，新建一個(gè)Element Set，因?yàn)?/span>Node Set無法輸出歷史數(shù)據(jù)，那么采用Element Set也肯定能夠輸出單元積分點(diǎn)的數(shù)據(jù)，如圖4所示。

圖3

圖4

之后進(jìn)行計(jì)算，計(jì)算完成后，進(jìn)入后處理模塊，點(diǎn)擊Create XY data，如圖5所示，選擇創(chuàng)建方式為ODB history output，彈出對(duì)話框如圖6所示。

圖5

圖6中第一行是Element Set的單元積分點(diǎn)的歷史數(shù)據(jù)，此次因?yàn)槭且粋€(gè)平面模型而且采用的縮減積分單元，每個(gè)單元只有一個(gè)積分點(diǎn)。

后面兩行是Node Set的節(jié)點(diǎn)位移U1和U2。

圖6

以上是歷史輸出的操作，但是也明顯存在了一個(gè)問題，如果想要輸出某節(jié)點(diǎn)的應(yīng)力怎么辦？

雖然采用單元積分點(diǎn)的值其實(shí)已經(jīng)能夠滿足通常的使用要求，但是單元積分點(diǎn)畢竟沒有節(jié)點(diǎn)直觀。

這里有兩種方法，第一種是通過Python語言，寫后處理程序，輸出節(jié)點(diǎn)的所有相關(guān)數(shù)據(jù)。

另外一種是通過修改inp文件，這種方法筆者沒有試過，是網(wǎng)上提到的一種方法，具體操作如下：

1.記住要輸出的單元編號(hào)，例如為27，在INP文件的*Assembly部分定義一個(gè)element型的set：

*Elset, elset=Set-Name, instance=Part-Name

 27,

2.在*End Step前加上：

*Output, history, frequency=10

*Element Output, elset=Set-Elem-forMeshPart

MISES,