abaqus數(shù)據(jù)探針的案例
abaqus cae數(shù)據(jù)文件損壞如何修復(fù)? abaqus data
將jnl文件復(fù)制一份，重命名后綴為.py的文件，打開abaqus，file-->run script, 運行一遍該py文件。但要注意你原始導(dǎo)入的幾何或孤立網(wǎng)格文件還在原來的文件夾位置下。
abaqus數(shù)據(jù)傳遞 ￥2
最近在學(xué)習(xí)abaqus數(shù)據(jù)傳遞功能。abaqus中數(shù)據(jù)傳遞的方法大體有三種，一種是我們常用的重啟動，一種是數(shù)據(jù)傳遞，還有一種是提取初始應(yīng)力場再導(dǎo)入。
1.重啟動的步驟如下：
1)在原模型中設(shè)置restart輸出請求；
2）在重啟動模型中設(shè)置重啟動請求：
單擊菜單Model/edit attributes,選擇重啟動模型名稱，設(shè)置重啟動分析步名稱以及重啟動迭代步；
基于Python向Abaqus導(dǎo)入txt、dat數(shù)據(jù)（附abaqus中python二次開發(fā)課程）
這次推送聚焦于解決采用Python向Abaqus里導(dǎo)入txt、dat數(shù)據(jù)的問題（dat文件只需要將txt文件的后綴名改為dat就可以生成dat文件），Abaqus基于Python讀入txt、dat數(shù)據(jù)主要有read()、readlines()、readlines()、numpy.loadtxt()函數(shù)，導(dǎo)入的txt文件放在工作目錄temp中。
1. read() 函數(shù)
可以一次性讀取txt文件中的內(nèi)容，并以字符串的結(jié)果返回來，如下列代碼所示：
with open("xyplot.dat", "r") as f: #打開文件
data = f.read() #讀取文件
print(data)
其中，
r是讀取人工書寫的數(shù)據(jù)，書寫的時候是什么樣子，讀出來就是什么樣。
rb是讀取二進(jìn)制文件，非人工書寫的數(shù)據(jù)，如.jpeg等這些。
將該命令在Abaqus/CAE的命令行接口（Command Line Interface, CLI）運行后，結(jié)果如圖1所示（圖中#注釋后的亂碼為作者安裝abaqus采用英文版，沒有漢化所致）：
圖1 read()讀取txt、dat數(shù)據(jù)到Abaqus
需要注意的是，以上data數(shù)據(jù)類型都是為字符，如下圖所示。調(diào)用data[0]為字符，因此data[0][0]也為字符，以及后續(xù)的data[20][0]，data[60][0]等也都為字符。
展開 Abaqus Python 二次開發(fā)-歷程輸出數(shù)據(jù)提取的python實現(xiàn) ￥3.49
在Abaqus中，歷程變量的輸出，尤其是反力及位移，經(jīng)常是如下的格式。
此時，如果想提取數(shù)據(jù)，一些關(guān)于Abaqus 二次開發(fā)的書籍里提供的方式是：odb.step['Step-1'].historyRegions['Node PART-TRIMMED-MESH.288422']. historyOutputs['U1'].data 。這種方式，需要輸入 'Node PART-TRIMMED-MESH.288422' 這種很長的字符串，繁瑣且容易出錯。
在付費內(nèi)容中，提供了一種便捷的方法，得到長字符串，提取相關(guān)的數(shù)據(jù)，并將數(shù)據(jù)保存為csv文件。
展開 
abaqus 運行完的ODB 數(shù)據(jù)提取
為什么之前運行完的ODB 文件可以正常提取數(shù)據(jù)，現(xiàn)在再對這個文件，在相同的位置卻提出不出來數(shù)據(jù)，就是一條橫線，這是什么問題呢
Abaqus的歷史數(shù)據(jù)輸出的問題
Abaqus的歷史數(shù)據(jù)輸出的問題
Abaqus結(jié)果輸出時有場輸出（Field Output）和歷史輸出(History Output)，場變量輸出用于描述某個量隨空間位置的變化，歷史變量用于描述某個量隨時間的變化，區(qū)別是場變量輸出大量的單元或節(jié)點上的計算結(jié)果，寫入odb文件的頻率低，用與生成后出的各個圖。歷史變量輸出少量單元或節(jié)點上的計算結(jié)果，寫入Odb的頻率高，用于生成X-Y圖。歷史變量允許單獨輸出某個獨立分量，經(jīng)常用到的就是這個獨立分量的輸出，比如輸出結(jié)構(gòu)最大點的位移變化曲線，或者盈利最大點的應(yīng)力變化過程等。
在輸出歷史數(shù)據(jù)時，大多數(shù)的輸出結(jié)果比如位移和應(yīng)力等，不能實現(xiàn)整個模型全部輸出，也就是需要指定set（參考集），需要設(shè)置相關(guān)的set后，結(jié)果輸出才會有這些數(shù)據(jù)。
定義Set的對話框如圖1所示，定義Set的類型時有三個，按照以往的理解，一般經(jīng)常用到的應(yīng)該是Node，也即節(jié)點集，想要輸出這個節(jié)點在求解過程中的Mises應(yīng)力變化過程，選擇某個節(jié)點，通過Step模塊的History Output Requests，如圖2所示。
圖1
圖2
選擇剛才設(shè)置的Set，選擇需要輸出的結(jié)果項，此處為Mises。之后計算，計算后提取結(jié)果時，發(fā)現(xiàn)并沒有輸出該節(jié)點的Mises歷史數(shù)據(jù)。
筆者很快想到了，有可能是不能輸出節(jié)點的應(yīng)力，因為涉及到節(jié)點解和單元積分點解的問題，這個之前已經(jīng)討論過。基于此考慮，那么輸出該節(jié)點的位移場肯定沒問題，因為節(jié)點位移解是最先得到的節(jié)點解。
將Mises換成UT，如圖3所示。
同時，新建一個Element Set，因為Node Set無法輸出歷史數(shù)據(jù)，那么采用Element Set也肯定能夠輸出單元積分點的數(shù)據(jù)，如圖4所示。
展開 ABAQUS批量提交Job與Python讀取ODB結(jié)果應(yīng)力應(yīng)變數(shù)據(jù)
問題2
如何從ODB結(jié)果文件中讀取計算結(jié)果數(shù)據(jù)
第二個問題是大家經(jīng)常要面對的，常規(guī)的數(shù)據(jù)提取采用ABAQUS自帶的后處理工具已經(jīng)能滿足需求了，但是如果我們需要進(jìn)行一些寫出文本，數(shù)據(jù)篩選，數(shù)據(jù)的二次計算處理。。。，我們就需要用程序語言的幫助了。
針對今天的問題——
提取結(jié)果并寫出到文本
，如果用自帶的工具，可以導(dǎo)出abaqus.rpt文件，不過這種方式并不能滿足我們較為苛刻的要求，用python程序會的自由度會更大一些。
另外需要說明的是，如果我們僅僅只是看看某單元或者節(jié)點的各種計算結(jié)果，或者簡要的數(shù)據(jù)處理，臨時用一用，可以采用如下方式，不必費心寫復(fù)雜的程序，不值當(dāng)浪費時間，
青春易逝，頭發(fā)寶貴
。當(dāng)然了，通過阿信這里給出的基本程序改一改，能省幾根頭發(fā)。
限于時間關(guān)系，這里只列出基本的應(yīng)力提取程序，寫入文本比較簡單，隨意百度即可實現(xiàn)，不再給出。
展開 星辰技文|Abaqus中提取裂縫數(shù)據(jù)并用matplotlib庫繪圖
目前在Abaqus中，基于全局或局部嵌入Cohesive單元，以模擬非均質(zhì)材料的裂縫擴展的方法已經(jīng)相當(dāng)普遍。我想POLARIS_InsertCohElem插件起到不小的作用。
后處理方面，也推出的POLARIS_CrackGeo插件提取Cohesive單元和XFEM模擬獲得的裂縫數(shù)據(jù)。但如何出圖來展示裂縫形態(tài)，成為插件用戶的一大痛點。
在Abaqus中，Cohesive單元模擬的裂縫路徑可以通過顯示特征邊的方式進(jìn)行展示，但這種方法有幾個方面的缺點：
1. 雖然可以疊加顯示出裂縫周邊實體單元的應(yīng)力、孔壓等場量輸出，但在表現(xiàn)裂縫自身場量結(jié)果時沒有線圖直觀，如下例中顯示縫寬的效果對比；
2. 很難像線圖那樣凸顯天然裂縫或顆粒邊界以及其它特殊位置的Cohesive單元；
3. 由于全局嵌入Cohesive單元，在單元共節(jié)點位置存在孔洞，部件變形后，孔洞呈現(xiàn)出小黑點的形式，無法去除，線圖就完全不存在這些問題；
常用的編程繪圖工具，目前以Matlab和Python matplotlib為主，Abaqus2021版本之后就已經(jīng)內(nèi)置了matplotlib庫，因此本文以matplotlib庫為基礎(chǔ)，帶大家繪制POLARIS_CrackGeo插件提取的裂縫線圖。
展開 ABAQUS odb數(shù)據(jù)解析系列文章---總綱
“使用python進(jìn)行ABAQUS后處理的二次開發(fā)是非常便捷的, 可以通過自動解析實現(xiàn)在一定程度上提高工作效率. 此系列文章將著手于odb文件的結(jié)構(gòu)特征, 進(jìn)而使得在后處理插件開發(fā)過程中少走彎路, 少做重復(fù)性工作.”
python是一種面向?qū)ο蟮母呒壵Z言, 那么在使用python處理odb的時候就要著手于對象, 將odb文件的對象譜系梳理清晰. 在完成上述工作之后, 再進(jìn)行后處理二次開發(fā)的時候就可以直接按圖索驥, 快速找到自己需要使用的方法了.
01
—
解析過程中的重要工具
ABAQUS為其定義的所有類提供了兩個特別有用的方法:
object.__methods__: 可以查看該類內(nèi)部的所有方法
object.__members__: 可以查看該類內(nèi)部的所有屬性
這兩個方法為梳理ABAQUS對象譜系起到了很大的作用, 避免了使用python自帶的dir方法(因其將對象內(nèi)部所有的方法與屬性同時返回, 導(dǎo)致解析困難).
ABAQUS還提供了優(yōu)化了的print方法, 該方法位于ABAQUS的textRepr中, 名為prettyPrint.
展開 Abaqus利用region讀取結(jié)果數(shù)據(jù)
Abaqus利用region讀取結(jié)果數(shù)據(jù)
長安CAE
Abaqus計算完成后，為輸出某部分模型的結(jié)果數(shù)據(jù)，需要用到region，并且配套使用的是getSubset函數(shù)。
示例一創(chuàng)建一個變量center，將一個叫做PUNCH的節(jié)點set賦值給這個變量，通過整個模型的位移場displacement得到該節(jié)點set的輸出數(shù)據(jù)，代碼如下：
center = odb.rootAssembly.instances['PART-1-1'].nodeSets['PUNCH']
centerDisplacement = displacement.getSubset(region=center)
centerValues = centerDisplacement.values
for v in centerValues:
print v.nodeLabel, v.data
其中的displacement是已經(jīng)創(chuàng)建好的位移場數(shù)據(jù)，centerDisplacement變量得到該節(jié)點集的所有數(shù)據(jù)，centerValues變量得到該節(jié)點集的數(shù)值數(shù)據(jù)，通過循環(huán)輸出節(jié)點號和當(dāng)前數(shù)據(jù)。
展開 Abaqus利用region讀取結(jié)果數(shù)據(jù)
單元質(zhì)心點
轉(zhuǎn)自公眾號——ABAQUS大世界
旨在分享，若侵即刪.
用Python從Abaqus導(dǎo)出txt、dat數(shù)據(jù)
__biz=MzIzNzcyMDE3OA==&mid=2247484080&idx=1&sn=1b34a15a8bf79e8f9cdf901f43aac9b1&chksm=e8c505addfb28cbb54225ef35bf7c6750ee21673d696b63fe13df9067320f1215cae6861c6a5&token=667787181&lang=zh_CN#rd</strong></a></p><p><strong>以下內(nèi)容均來自于本人微信公眾號（圖片版可能缺失內(nèi)容，完整版請關(guān)注查看公眾號內(nèi)容），有興趣請關(guān)注一波，共同學(xué)習(xí)ANSYS和ABAQUS有限元技術(shù)，我理解為兩把神器，屠龍刀和倚天劍，用好可行走有限元江湖，做個大俠。</strong></p><p><strong>同濟大學(xué)優(yōu)秀畢業(yè)生創(chuàng)立，研究Abaqus、ANSYS在土木領(lǐng)域的科研編程和數(shù)值模擬技術(shù)，涉及到多尺度、多物理場耦合等計算力學(xué)問題，以及基于Python、子程序、INP的Abaqus建模，和ANSYS APDL、Workbench二次開發(fā)。
展開 abaqus 獲取odb場輸出中某一集合的數(shù)據(jù)
獲取odb中某一個集合的數(shù)據(jù)，有代碼實現(xiàn)和cae實現(xiàn)
代碼實現(xiàn)：主要的是一個getSubset的函數(shù)，可以參考abaqus手冊
https://help.3ds.com/2022/english/dssimulia_established/SIMACAECMDRefMap/simacmd-c-odbintroreadregpyc.htm?
如何在abaqus中得到荷載-位移的數(shù)據(jù)
在variables選項卡中的position下拉框里選擇unique nodal 在下面的變量里勾選RF或RT（反力）、U（位移）一般只選某個方向的（如2方向）；在elements/nodes選項卡中的method選擇Node sets，右邊選擇你定義的參考點 點擊Save
5.這時在操作框里XY Data欄下會有兩個數(shù)據(jù)，他們是參考點處的反力和位移隨時間的變化，在右邊的operators里有一個函數(shù)combine(x,x),點一下這個函數(shù)會出現(xiàn)在expression欄里，將兩個數(shù)據(jù)位移和反力用add to expression添加到combine函數(shù)的括號里，注意位移在前，反力在后，中間的逗號是英文的“,”
6.將expression另存為（save as按鈕）一個新的名字，可以用plot expression查看曲線，也可以在主窗口的XY Data manager用plot查看，用edit讀取數(shù)值
如果覺得位移和反力的符號是相反的，可以在第5步combine之前將兩個數(shù)據(jù)反號另存為新的數(shù)據(jù)之后combine
展開 晶體塑性模擬，EBSD數(shù)據(jù)導(dǎo)入abaqus
[圖片]
