abaqus提取結(jié)果
關(guān)注創(chuàng)建者：王靖雯 創(chuàng)建時(shí)間：2023-02-27
abaqus提取結(jié)果的視頻教程
sysweld軟件常用后處理提取結(jié)果操作
本教程是本人長期使用該軟件經(jīng)驗(yàn)積累，該教程基本覆蓋了常規(guī)結(jié)果的提取方法，以及怎樣對結(jié)果進(jìn)行處理，這是對寫論文出成果至關(guān)重要的。
￥20 43分鐘 80播放
查看
abaqus提取結(jié)果的實(shí)例教程
采用Python可以提取Abaqus計(jì)算結(jié)果的歷史輸出數(shù)據(jù)，提取出來后保存到本地并采用Matplitlib繪制折線圖，記錄一下代碼：以下三張圖片分別是讀取etotal歷史數(shù)據(jù)并保存到本地的代碼、繪制折線圖的代碼和etotal歷史變化折線圖。
上一篇技文《ABAQUS結(jié)果提取大于某值的區(qū)域體積-CAE方法》中帶大家使用ABAQUS CAE界面直接提取大于100e6 Mises等效應(yīng)力的區(qū)域體積，方法雖然比較好操作，但也存在明顯的缺點(diǎn)：
方法太繁瑣；統(tǒng)計(jì)歷程曲線時(shí)會(huì)顯得非常麻煩，因此我們找到了代步工具：Python。
精度較低；這是由于一個(gè)單元內(nèi)只需要一個(gè)積分點(diǎn)滿足數(shù)值要求，整個(gè)單元都會(huì)顯示出來，但一個(gè)單元可能存在多個(gè)積分點(diǎn)，比如演示案例中采用六面體二階減縮積分單元（C3D20R）存在8個(gè)積分點(diǎn)，單元內(nèi)部分積分點(diǎn)可能不滿足條件，也會(huì)被統(tǒng)計(jì)在內(nèi)，從而導(dǎo)致統(tǒng)計(jì)的體積偏大，我們對比下CAE方法和Python方法的結(jié)果差異如下圖，最少都差32%！
樣的誤差雖然可以通過提高網(wǎng)格密度來減緩，但并不能完全避免，因此也就有必要對這些單元內(nèi)所有積分點(diǎn)值進(jìn)行判斷，獲得單元內(nèi)滿足條件的積分點(diǎn)所占比例，再對其體積進(jìn)行加權(quán)相加。而這樣的操作過程顯然不適合手動(dòng)統(tǒng)計(jì)，再一次將目光投向了我們的代步工具：Python！
沒有Python基礎(chǔ)的小伙伴，建議先看曹金鳳姐姐的《Python語言在Abaqus中的應(yīng)用》或江丙云哥哥的《ABAQUS Python二次開發(fā)攻略》，不然會(huì)有些吃力。雖然代碼備注的已經(jīng)非常詳細(xì)，但還是需要一些背景支撐的。
【編程需求與分解】
為了與CAE過程進(jìn)行對比，我們還是將需求定義為：提取懸臂梁加載過程中Mises等效應(yīng)力大于100e6的區(qū)域體積。
展開 在有些場景下，需要將Abaqus計(jì)算結(jié)果提取出來保存為VTK格式，用于一些可視化和數(shù)據(jù)處理（VTK在可視化這方面還是強(qiáng)大）。最近寫了兩個(gè)腳本，一個(gè)是保存VTK非結(jié)構(gòu)網(wǎng)格格式的文件，一個(gè)是提取Abaqus計(jì)算結(jié)果中的網(wǎng)格信息、field信息，將兩個(gè)腳本結(jié)合起來就可以實(shí)現(xiàn)提取ODB中的信息，并保存成為VTK非結(jié)構(gòu)網(wǎng)格文件，當(dāng)然，還是需要有Abaqus環(huán)境。
注意：提取ODB信息的那個(gè)腳本目前僅僅是針對我自己提供的那個(gè)算例，有能力的朋友完全可以修改一下用于自己的情況。寫VTK非結(jié)構(gòu)網(wǎng)格文件的那個(gè)腳本目前應(yīng)該沒有太大的bug。ODB信息的提取還是需要有Abaqus環(huán)境，并不是解析二進(jìn)制文件。請看清楚再確定是否需要。
之所以做這個(gè)，還是希望慢慢完善，從而積累出一定數(shù)量的腳本，用于ODB數(shù)據(jù)的提取。
以下是abaqus中的網(wǎng)格和云圖，和用paraview查看的vtk文件的網(wǎng)格和云圖：
展開 ABAQUS 顯式explicit計(jì)算時(shí)所有結(jié)果都將寫入到.odb文件中，在GUI界面進(jìn)行數(shù)據(jù)提取時(shí)效率較低，且有些詳細(xì)結(jié)果不能直接從.odb中獲取。因此，需要利用python讀取.odb，然后進(jìn)行提取關(guān)心的數(shù)據(jù)，以此進(jìn)行下一步計(jì)算。
Python腳本操作ODB文件進(jìn)行計(jì)算結(jié)果數(shù)據(jù)的提取。
因?qū)?em>結(jié)果分析需求，需提取某單元集的每一分析步（包含每一子步）的平均應(yīng)力，目前網(wǎng)上雖有很多代碼是關(guān)于單元集的平均應(yīng)力的提取，但并未有針對每一子步都需要提取結(jié)果的代碼。故針對此需要編寫了python代碼。
abaqus提取結(jié)果的相關(guān)專題、標(biāo)簽、搜索
abaqus提取結(jié)果的最新內(nèi)容
一分鐘上手的ABAQUS結(jié)果二次處理工具6個(gè)月前
我的專業(yè)一直是結(jié)構(gòu)力學(xué)方向。陰差陽錯(cuò)之下，畢業(yè)后搞了三年的CFD方向。
剛接觸CFD那會(huì)我就很感慨，流體的人玩的真花。前處理，光網(wǎng)格搞了一堆花樣，結(jié)構(gòu)非結(jié)構(gòu)，多面體，附面層，還有什么泥瓦匠網(wǎng)格思路和雕刻家網(wǎng)格思路。為了收斂，前處理可以說是不擇手段。
后處理更是如此，各種流線、壁面、矢量、渦，穩(wěn)態(tài)的、動(dòng)態(tài)的。光專門用來做后處理的軟件工具都一堆，Tecplot、CFD-post、Insight
CFD結(jié)果映射到ABAQUS模型上6個(gè)月前
在流體到固的單向耦合仿真中，常常需要傳遞的物理量是壓力和溫度，它們可以作為結(jié)構(gòu)分析的載荷條件。
不同軟件之間數(shù)據(jù)格式不通，需要人為的去做數(shù)據(jù)轉(zhuǎn)換。更關(guān)鍵的是，在實(shí)際項(xiàng)目中，做CFD仿真的和做CAE仿真的壓根不是一個(gè)人。
這會(huì)帶來兩個(gè)問題：
(1) 流體壁面和結(jié)構(gòu)壁面模型坐標(biāo)系不一致，參數(shù)無法直接用；
(2) 流體壁面和結(jié)構(gòu)壁面玩?zhèn)€節(jié)點(diǎn)不一致，無法進(jìn)行點(diǎn)對點(diǎn)的參數(shù)傳遞。
所謂逢山開路遇水搭橋
Abaqus：激光輔助車削仿真結(jié)果對比解析6個(gè)月前
[圖片]
Abaqus平均應(yīng)力和應(yīng)變提取7個(gè)月前
利用python腳本對ODB文件中單元集里所有積分點(diǎn)的應(yīng)力及應(yīng)變進(jìn)行自動(dòng)提取并計(jì)算平均值
能夠得到每一幀的應(yīng)力和應(yīng)變平均值，并保存到CSV文件中
所得到的應(yīng)力包括S11,S22,S33,S12,S13,S23以及Mises七個(gè)應(yīng)力平均值，以及E11,E22,E33,E12,E13,E23六個(gè)應(yīng)變平均值
[圖片]
<p><span style="background-color: rgba(0, 0, 0, 0);">粘結(jié)滑移行為是結(jié)構(gòu)工程精細(xì)化有限元的重要部分。常見的應(yīng)用場景例如FRP加固結(jié)構(gòu)中的粘結(jié)界面、新舊混凝土疊合面、鋼筋和混凝土的粘結(jié)界面等。相比綁定/內(nèi)置，考慮粘結(jié)滑移的模型可能出現(xiàn)變形增大，峰值位移相對滯后，剛度相對較低，滯回曲線捏縮等宏觀現(xiàn)象，但對模型本身粘結(jié)力和滑移大小等微觀行為的查看使用者通常不知如何處理
如圖所示，采用軸對稱模型，Z軸為對稱軸，計(jì)算受力情況
1.模型說明
軸對稱模型，上面兩個(gè)線圈，圈數(shù)分別為100圈和150圈，下面為一個(gè)圓筒，材料為銅，那么在線圈中通過一個(gè)變化的電流，在下面的圓筒中就會(huì)產(chǎn)生感應(yīng)電流，感應(yīng)電流在在磁場中受到洛倫茲力的作用，查看隨時(shí)間變化的受力情況。
2.模型建立及加載
建立模型，給定材料，設(shè)置winding，添加電流為
模型描述：兩個(gè)件以10mm/ms的速度撞擊另外兩個(gè)用rigid約束起來的件，已設(shè)置重力加速度和剛性地面。各自配重10kg。自接觸使用*CONTACT_AUTOMATIC_SINGLE_SURFACE，包含所有shell單元，靜/動(dòng)摩擦系數(shù)0.2。
同時(shí)建立segment set使用*FORCE_TRANSDUCER_PENALTY提取自接觸力，一個(gè)件分面提取，一個(gè)件整體提取
也是在站內(nèi)找到別的大佬的評論看到的，分享給大家。直接使用查詢工具就可以，選擇“質(zhì)量屬性（mass properties）”
點(diǎn)擊后選擇幾何區(qū)域
在上邊選擇“面”
然后直接選擇需要得到面積的面點(diǎn)擊“完成”就可以啦！體積也是同樣的道理，在這就不演示啦！
附圖是站內(nèi)評論區(qū)看到的，再次感謝大佬
