abaqus 數(shù)據(jù)輸出的案例
Abaqus的歷史數(shù)據(jù)輸出的問題
Abaqus的歷史數(shù)據(jù)輸出的問題
Abaqus結(jié)果輸出時(shí)有場(chǎng)輸出（Field Output）和歷史輸出(History Output)，場(chǎng)變量輸出用于描述某個(gè)量隨空間位置的變化，歷史變量用于描述某個(gè)量隨時(shí)間的變化，區(qū)別是場(chǎng)變量輸出大量的單元或節(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)的位移場(chǎng)肯定沒問題，因?yàn)楣?jié)點(diǎn)位移解是最先得到的節(jié)點(diǎn)解。
將Mises換成UT，如圖3所示。
同時(shí)，新建一個(gè)Element Set，因?yàn)镹ode Set無法輸出歷史數(shù)據(jù)，那么采用Element Set也肯定能夠輸出單元積分點(diǎn)的數(shù)據(jù)，如圖4所示。
展開 Abaqus Python 二次開發(fā)-歷程輸出數(shù)據(jù)提取的python實(shí)現(xiàn) ￥3.49
在Abaqus中，歷程變量的輸出，尤其是反力及位移，經(jīng)常是如下的格式。
此時(shí)，如果想提取數(shù)據(jù)，一些關(guān)于Abaqus 二次開發(fā)的書籍里提供的方式是：odb.step['Step-1'].historyRegions['Node PART-TRIMMED-MESH.288422']. historyOutputs['U1'].data 。這種方式，需要輸入 'Node PART-TRIMMED-MESH.288422' 這種很長(zhǎng)的字符串，繁瑣且容易出錯(cuò)。
在付費(fèi)內(nèi)容中，提供了一種便捷的方法，得到長(zhǎng)字符串，提取相關(guān)的數(shù)據(jù)，并將數(shù)據(jù)保存為csv文件。
展開 abaqus 獲取odb場(chǎng)輸出中某一集合的數(shù)據(jù)
獲取odb中某一個(gè)集合的數(shù)據(jù)，有代碼實(shí)現(xiàn)和cae實(shí)現(xiàn)
代碼實(shí)現(xiàn)：主要的是一個(gè)getSubset的函數(shù)，可以參考abaqus手冊(cè)
https://help.3ds.com/2022/english/dssimulia_established/SIMACAECMDRefMap/simacmd-c-odbintroreadregpyc.htm?
利用Python解析Abaqus的ODB數(shù)據(jù)并輸出為VTK格式 ￥199
在做結(jié)構(gòu)分析的時(shí)候，做出來的結(jié)果在一些情況下可能需要提取稱為其他格式的數(shù)據(jù)，比如vtk格式，能夠提供更靈活的數(shù)據(jù)分析能力。
于是開發(fā)了一些代碼用于這些數(shù)據(jù)的提取，部分Abaqus展示效果圖和在Paraview中的效果圖展示如下：
需要PC安裝了Abaqus軟件并且支持Python Script接口；目前本人使用的版本是Abaqus/CAE 2021版本。
目前支持C3D4、C3D10、C3D8、C3D8R單元類型。
提取ODB文件中所有frame的數(shù)據(jù)，每一幀對(duì)應(yīng)一個(gè)結(jié)果文件；提取方法為遍歷所有frame，獲取所有fieldOutput，獲取fieldOutput上的數(shù)據(jù)，SCALAR和VECTOR類型的數(shù)據(jù)是定義在節(jié)點(diǎn)上的，直接獲取POINT_DATA即可。
目前對(duì)于Tensor類型的數(shù)據(jù)，由于從場(chǎng)數(shù)據(jù)中提取得到的是積分點(diǎn)上的數(shù)據(jù)，目前的處理方式是計(jì)算所有積分點(diǎn)數(shù)據(jù)的平均值，作為CELL_DATA進(jìn)行處理；其他數(shù)據(jù)作為POINT_DATA進(jìn)行處理。
目前的腳本文件結(jié)果如下所示：
展開 
Abaqus_Python 后處理結(jié)點(diǎn)數(shù)據(jù)提取、場(chǎng)變量云圖輸出 ￥1.99
今天體會(huì)了Python腳本操作ODB文件進(jìn)行數(shù)據(jù)的輸出、處理等方面的快速高效。
首先介紹一下背景，同事的項(xiàng)目，計(jì)算case較多。每個(gè)case需要輸出幾個(gè)路徑上的部分節(jié)點(diǎn)值，而且需要針對(duì)不同的case、不同的幀、不同的視角截圖、保存云圖，如果人工出圖和提取節(jié)點(diǎn)數(shù)據(jù)，工作量會(huì)非常大，費(fèi)時(shí)會(huì)遠(yuǎn)遠(yuǎn)超出計(jì)算所用時(shí)間；并且，同事還要求我根據(jù)case名稱、組件名稱（模型有多個(gè)組件組成，且網(wǎng)格是在part上畫的，所以同一個(gè)節(jié)點(diǎn)編號(hào)可能對(duì)應(yīng)多個(gè)節(jié)點(diǎn)）等保存成不同的文件。
結(jié)點(diǎn)溫度數(shù)據(jù)沒有保存，因此沒有輸出。下圖是自動(dòng)輸出的png格式云圖。
比較倉(cāng)促，代碼也比價(jià)粗糙。回頭做一個(gè)詳細(xì)的介紹。
主要代碼見下。
展開 Python二次開發(fā)—?jiǎng)?chuàng)建輸出數(shù)據(jù)庫(kù)并添加數(shù)據(jù)
本期內(nèi)容轉(zhuǎn)載本人公眾號(hào)：
易木木響叮當(dāng)，源代碼可在公眾號(hào)內(nèi)回復(fù)：
創(chuàng)建數(shù)據(jù)庫(kù)，即可自動(dòng)獲取
。
Python二次開發(fā)更新啦~
最近由于木木學(xué)業(yè)壓力比較大，故更文速度有點(diǎn)小慢，請(qǐng)大家諒解諒解。廢話不多說，進(jìn)入今天的正題：
本期內(nèi)容將編寫腳本創(chuàng)建輸出數(shù)據(jù)庫(kù)，并向其中添加模型數(shù)據(jù)、場(chǎng)數(shù)據(jù)和歷史數(shù)據(jù)，在后處理模塊中顯示位移場(chǎng)變量計(jì)算結(jié)果，最后保存ODB文件。源腳本來自于曹老師的《Python語言在Abaqus中的應(yīng)用》，好書值得推薦！
編程思路：
導(dǎo)入相應(yīng)模塊。本例中將創(chuàng)建輸出數(shù)據(jù)庫(kù)，并向其中添加數(shù)據(jù)，于是應(yīng)添加：odbAccess、odbMaterial、odbSection、和abaqusConstants等模塊;
創(chuàng)建輸出數(shù)據(jù)庫(kù)時(shí)，數(shù)據(jù)庫(kù)不包含任何模型需要調(diào)用相應(yīng)函數(shù)構(gòu)造對(duì)象。一般情況下需要?jiǎng)?chuàng)建以下對(duì)象：材料屬性（Material模塊）、截面（Section模塊）、部件實(shí)例（Assembly模塊）、分析步和幀（Step模塊）、節(jié)點(diǎn)和單元（Mesh模塊）；
向輸出數(shù)據(jù)庫(kù)中添加場(chǎng)數(shù)據(jù)和歷程數(shù)據(jù)，本例中只向輸出數(shù)據(jù)庫(kù)中添加了節(jié)點(diǎn)位移，場(chǎng)輸出為默認(rèn)輸出；
創(chuàng)建完數(shù)據(jù)庫(kù)后，應(yīng)調(diào)用save( )方法保存數(shù)據(jù)庫(kù)文件，然后調(diào)用close( )方法關(guān)閉文件。
在Abaqus操作如下：點(diǎn)擊File，Run Script，
若CAE界面下方出現(xiàn)：
New output database has been created successfully!
You can view the odb_Create_ODB.odb file in the visualization module!
展開 關(guān)于LS-DYNA 的輸出數(shù)據(jù) ￥30
</p><p><br></p><p><br></p><p><strong style="color: rgb(26, 28, 30); background-color: rgb(252, 252, 252);">ASCII格式</strong><span style="color: rgb(26, 28, 30); background-color: rgb(252, 252, 252);">的結(jié)果文件，主要用于生成</span><strong style="color: rgb(26, 28, 30); background-color: rgb(252, 252, 252);"><em>二維XY曲線圖：</em></strong></p><p><strong> </strong> ASCII文件由關(guān)鍵字＊DATABASE_OPTION控制輸出，使用較多的是GLSTAT（總體統(tǒng)計(jì)和能量）、MATSUM（材料能量總結(jié)）、RBDOUT（主要記錄剛體的運(yùn)動(dòng)狀態(tài)和力學(xué)響應(yīng)，如剛體的位移、速度、加速度、角速度、角加速度等等）、RCFORC（接觸面反作用力），還有FSI（輸出dbfsi文件，耦合力輸出），NCFORC（輸出單個(gè)節(jié)點(diǎn)的接觸力）這六個(gè)文件。</p><p><br></p><p><br></p><p>二進(jìn)制的文件可以使用LS-PrePost軟件打開，輸出帶數(shù)據(jù)的動(dòng)畫等。 ASCII文件就是使用ASCII直接打開（各種數(shù)據(jù)曲線），使用MPP計(jì)算的文件，需要在Binout里打開（binout0000文件）。</p><p><br></p><p><br></p><p><br></p><p><br></p><p>收費(fèi)部分部分預(yù)覽，就是將常見的輸出關(guān)鍵字整合到一起了，方便使用。
展開 70基于matlab的BP神經(jīng)網(wǎng)絡(luò)多輸入單輸出數(shù)據(jù)結(jié)果預(yù)測(cè)，輸出結(jié)果包括均方根誤差，決定系數(shù)。 ￥25.9
基于matlab的BP神經(jīng)網(wǎng)絡(luò)多輸入單輸出數(shù)據(jù)結(jié)果預(yù)測(cè)，輸出結(jié)果包括均方根誤差，決定系數(shù)。數(shù)據(jù)可更換自己的，程序已調(diào)通，可直接運(yùn)行。
衍射光學(xué)元件的產(chǎn)生及其結(jié)構(gòu)數(shù)據(jù)的輸出
結(jié)構(gòu)數(shù)據(jù)的產(chǎn)生
2 點(diǎn)擊Ok按鈕
結(jié)構(gòu)數(shù)據(jù)的產(chǎn)生
2 出現(xiàn)一個(gè)保存對(duì)話框
2 選擇你要保存的文件夾輸入有效文件名
2 所有的數(shù)據(jù)文件可以有相同的名字而不同的擴(kuò)展名
2 文件包含的蝕刻掩模(etching mask)的數(shù)據(jù)將在文件名中有附加的掩模（層layer）數(shù).Mask 1表示最深的蝕刻步.
輸出數(shù)據(jù)結(jié)果
數(shù)據(jù)輸出結(jié)果
2 打開Windows Explorer到VL_Samples文件夾
2 雙擊DO.8_Generation_of_DOE_Fabrication_Data_02.xml file.
2 Internet Explorer將顯示文件
數(shù)據(jù)輸出結(jié)果
2 該XML文件包含了所有參數(shù)的匯總并連接所有產(chǎn)生的表面數(shù)據(jù)文件
結(jié)論
2 VirtualLab包含一個(gè)從傳輸器（transmission）中計(jì)算高度分布的結(jié)構(gòu)設(shè)計(jì)
2 高度分布可以被輸出。VirtualLab包含結(jié)構(gòu)數(shù)據(jù)的不同文件格式（GDSII，ASCII，Bitmap格式）。
展開 ANSYS數(shù)據(jù)輸出格式Fw.d的應(yīng)用
ANSYS數(shù)據(jù)輸出格式Fw.d的應(yīng)用
ANSYS計(jì)算后處理時(shí)經(jīng)常需要將得到的結(jié)果進(jìn)行輸出，輸出時(shí)需要采用一定的輸出格式指定輸出方式。數(shù)據(jù)輸出一般配合*CFOPEN和*VWRITE完成，*CFOPEN用于打開文件，*VWRITE用于寫數(shù)據(jù)。
APDL的輸出格式和Fortran一致其中F格式為：Fw.d。這個(gè)用的比較多，用于輸出浮點(diǎn)數(shù)據(jù)。
其中，w表示數(shù)據(jù)的總字符寬度，d表示小數(shù)部分所占的寬度，不夠的補(bǔ)零。例如F10.5表示輸出數(shù)據(jù)一共占10個(gè)寬度，其中小數(shù)部分占5個(gè)寬度，需要注意點(diǎn)號(hào)也占一個(gè)字符寬度，不夠10位的在數(shù)據(jù)的前面補(bǔ)空格。
展開 骨骼掃描數(shù)據(jù)成像與有限元網(wǎng)格輸出
Simpleware軟件是英國(guó)的Simpleware Ltd公司推出的一套集成逆向工程、材料工程、生物力學(xué)工程、有限元分析等多工業(yè)、多學(xué)科領(lǐng)域的統(tǒng)一解決方案的專業(yè)軟件，它通過對(duì)三維體圖或二維平面序列圖像（CT/MRI/Microscopy）進(jìn)行提取、分割、合并等操作，生成表面重建后的三維圖像，并將其輸出到第三方軟件進(jìn)行快速成型制造；同時(shí)提供多種體網(wǎng)格生成技術(shù)，直接生成多種有限元軟件或計(jì)算流體力學(xué)軟件的網(wǎng)格文件，并支持導(dǎo)出NURBS曲面描述的 CAD模型。
Simpleware軟件為骨科醫(yī)生進(jìn)行數(shù)字骨科實(shí)踐提供了理想的平臺(tái)
Simpleware股骨建模、植入物定位與有限元網(wǎng)格生成
使用Simpleware中的ScanIP模塊，將CT等掃描的股骨數(shù)據(jù)圖像文件，生成三維股骨模型。
中仿高級(jí)工程師Sean，與某大型醫(yī)院骨科王醫(yī)生密切配合，將生成的三維股骨，通過Simpleware 3D編輯工具修改模型，去除股骨頭，使用+CAD模塊導(dǎo)入和定位手術(shù)中的植入物CAD模型，從而模擬外科手術(shù)；使用+FE模塊，可以生成用戶需要的表面網(wǎng)格模型和體網(wǎng)格模型。
基于三重周期性隱函數(shù)，Simpleware為構(gòu)造立方體格子結(jié)構(gòu)，提供多種內(nèi)置單元，并能與醫(yī)學(xué)圖像完美融合，將修復(fù)區(qū)域替換為功能梯度材料。Simpleware強(qiáng)大易用的功能以及中仿的技術(shù)實(shí)施能力，都得到了王醫(yī)生的高度肯定。
文章來源：中仿科技
展開 
ANSYS后處理將數(shù)據(jù)以輸出到txt文本中應(yīng)用案例 ￥10
ANSYS后處理將數(shù)據(jù)以輸出到txt文本中，用到的主要命令為do循環(huán)、get命令；
循環(huán)命令*DO,Par,IVAL,FVAL,INC
Par循環(huán)變量的名稱，可以定義為i，j等
IVAL, Par循環(huán)變量的初始值
FVAL, Par循環(huán)變量的終止值
INC循環(huán)變量的增長(zhǎng)步長(zhǎng)，缺省值為1
例子1：*do,i,1,100
***********************************************************************************************************************
獲取命令*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM
Par變量的名稱，將獲取的值存到變量中去
Entity,獲取的對(duì)象類型可以為node、element、area等
ENTNUM為對(duì)象的ID號(hào)
Item1可以為單元表選項(xiàng)smic等
IT1NUM可以為單元表選項(xiàng)smic對(duì)應(yīng)的輸出序號(hào)，對(duì)于梁?jiǎn)卧梢杂脕?em>輸出載荷或應(yīng)力
例子1統(tǒng)計(jì)完單元數(shù)量賦值給變量number：*GET, number, ELEM,0,count
例子2提取單元最大的ID號(hào)賦值給變量number：*GET, number,ELEM,0,NUM,MAX
**********************************************************************************************************************
數(shù)組定義命令*DIM,Par,Type,IMAX,JMAX,KMAX,Var1,Var2,Var3,CSYSID
Par數(shù)組的名稱
Type
展開 umat等子程序用open和write輸出數(shù)據(jù),出現(xiàn)一直顯示運(yùn)行中，增量步不計(jì)算的現(xiàn)象。 ￥10
abaqus子程序中用open和write，close這些命令輸出自己感興趣的數(shù)據(jù)。 但有時(shí)候會(huì)一直顯示運(yùn)行中，abaqus沒有增量步計(jì)算。我發(fā)現(xiàn)大量論壇上很多人出現(xiàn)了這個(gè)問題，但一直未解決。如何解決這個(gè)問題，我嘗試了多種方法，調(diào)試了大量程序，最終發(fā)現(xiàn)可以通過以下方法解決。以下內(nèi)容包括
產(chǎn)生此現(xiàn)象”一直顯示運(yùn)行中，abaqus沒有增量步計(jì)算“的原因分析；
解決方案；
abaqus二次開發(fā)：后處理批量提取場(chǎng)輸出和歷程輸出結(jié)果（源碼帶注釋） ￥198
abaqus二次開發(fā)：后處理批量提取場(chǎng)輸出和歷程輸出結(jié)果（購(gòu)買后需要加V私聊獲取使用教程）
購(gòu)買后，請(qǐng)加V：wzd_1021_
更多插件內(nèi)容，請(qǐng)私聊
同時(shí)提供定制hypermesh/hyperview/abaqus前后處理插件開發(fā)。
6混淆矩陣，根據(jù)預(yù)測(cè)結(jié)果輸出混淆矩陣。將自己標(biāo)簽數(shù)據(jù)帶進(jìn)去即可得到結(jié)果，可更改顏色，更改預(yù)測(cè)的個(gè)數(shù) ￥8.9
混淆矩陣，根據(jù)預(yù)測(cè)結(jié)果輸出混淆矩陣。將自己標(biāo)簽數(shù)據(jù)帶進(jìn)去即可得到結(jié)果，可更改顏色，更改預(yù)測(cè)的個(gè)數(shù)，基于MATLAB平臺(tái)，程序已調(diào)通可直接運(yùn)行。
abaqus 數(shù)據(jù)輸出的相關(guān)專題、標(biāo)簽、搜索
abaqus數(shù)據(jù)輸出abaqus輸出數(shù)據(jù)abaqus 數(shù)據(jù)輸出abaqus中輸出數(shù)據(jù)abaqus輸出數(shù)據(jù)中文ansys 輸出數(shù)據(jù) Abaqus abaqus提取加速度數(shù)據(jù)用場(chǎng)輸出還是歷程輸出的數(shù)據(jù)abaqus場(chǎng)輸出與歷程輸出數(shù)據(jù)一樣abaqus輸出數(shù)據(jù)abaqus數(shù)據(jù)輸出abaqus能直接輸出歷史輸出數(shù)據(jù)嗎abaqus數(shù)據(jù)輸出excel
