abaqus提取體積的案例
ABAQUS提取面積或體積
體積也是同樣的道理，在這就不演示啦！
附圖是站內(nèi)評(píng)論區(qū)看到的，再次感謝大佬！
ABAQUS結(jié)果提取大于某值的區(qū)域體積（轉(zhuǎn)載）
經(jīng)常會(huì)有ABQer初學(xué)者會(huì)問我：“怎么獲得屈服區(qū)域的體積？”、“怎么獲得大于一定應(yīng)力的單元體積？”之類的一些問題。我們將帶大家探索ABAQUS后處理的兩種提取方法：CAE界面提取 和 Python后處理二次開發(fā)提取。本節(jié)先介紹第一種方法。
還是以最簡單的懸臂梁為例，提取加載歷程下大于100e6 Mises等效應(yīng)力下單元體積：
1）顯示應(yīng)力變形云圖，并通過云圖顯示設(shè)置，可將超過指定數(shù)值的云圖范圍顯式為灰色，便于觀察和對(duì)比：
2）顯示指定應(yīng)力范圍內(nèi)的單元：找到按鈕，或菜單欄Tools->DisplayGroup->Create...，點(diǎn)擊后激活Create Display Group對(duì)話框，選擇Item為Elements、method為Result value，并設(shè)置最小值（Min value）100e6和最大值（Max value）1e10，點(diǎn)擊Apply可顯示出范圍內(nèi)的單元，如下圖所示。
注意：同時(shí)點(diǎn)擊底部的Save As...按鈕，保存名字為DisplayGroup-2。
3）測量顯示單元體積：找到按鈕，或者菜單欄Tools->Query...命令，找到Mass properties選型，并點(diǎn)擊后激活查詢質(zhì)量流程，在提示信息欄將讓選擇單元，修改選擇方法為Display groups，并選擇剛剛生成的DisplayGroup-2；由于是測量質(zhì)量，針對(duì)材料中沒有定義密度的情況，ABAQUS友好的提供了Options選項(xiàng)，可以設(shè)置材料密度和平面shell的厚度，這里均設(shè)置為1。如下圖所示。
重新選擇顯式組后，在底部信息框中將顯示測量的體積、體積中心、質(zhì)量、質(zhì)量中心以及質(zhì)量距等數(shù)據(jù)。
展開 ABAQUS結(jié)果提取大于某值的區(qū)域體積-CAE方法
以最簡單的懸臂梁為例，提取加載歷程下大于100e6 Mises等效應(yīng)力下單元體積：
1）顯示應(yīng)力變形云圖，并通過云圖顯示設(shè)置，可將超過指定數(shù)值的云圖范圍顯式為灰色，便于觀察和對(duì)比：
2）顯示指定應(yīng)力范圍內(nèi)的單元：找到按鈕，或菜單欄Tools->DisplayGroup->Create...，點(diǎn)擊后激活Create Display Group對(duì)話框，選擇Item為Elements、method為Result value，并設(shè)置最小值（Min value）100e6和最大值（Max value）1e10，點(diǎn)擊Apply可顯示出范圍內(nèi)的單元，如下圖所示。
注意：同時(shí)點(diǎn)擊底部的Save As...按鈕，保存名字為DisplayGroup-2。
3）測量顯示單元體積：找到按鈕，或者菜單欄Tools->Query...命令，找到Mass properties選型，并點(diǎn)擊后激活查詢質(zhì)量流程，在提示信息欄將讓選擇單元，修改選擇方法為Display groups，并選擇剛剛生成的DisplayGroup-2；由于是測量質(zhì)量，針對(duì)材料中沒有定義密度的情況，ABAQUS友好的提供了Options選項(xiàng)，可以設(shè)置材料密度和平面shell的厚度，這里均設(shè)置為1。如下圖所示。
重新選擇顯式組后，在底部信息框中將顯示測量的體積、體積中心、質(zhì)量、質(zhì)量中心以及質(zhì)量距等數(shù)據(jù)。
4）匯總不同時(shí)刻的數(shù)值，繪制時(shí)程曲線如下：
同時(shí)，該方法不僅適用于應(yīng)力數(shù)值范圍，基本所有的云圖范圍單元體積的測量都是可以滿足的。
作者：陳佳敏cn
來源：CAE技術(shù)資訊
展開 星辰技文|ABAQUS結(jié)果提取大于某值的區(qū)域體積-Python方法 ￥10
上一篇技文《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ì)的體積偏大，我們對(duì)比下CAE方法和Python方法的結(jié)果差異如下圖，最少都差32%！
樣的誤差雖然可以通過提高網(wǎng)格密度來減緩，但并不能完全避免，因此也就有必要對(duì)這些單元內(nèi)所有積分點(diǎn)值進(jìn)行判斷，獲得單元內(nèi)滿足條件的積分點(diǎn)所占比例，再對(duì)其體積進(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)行對(duì)比，我們還是將需求定義為：提取懸臂梁加載過程中Mises等效應(yīng)力大于100e6的區(qū)域體積。
展開 
ANSYS提取具體三維單元的體積，面元的面積和線的長度
在ANSYS中，能提取具體三維單元的體積，面元的面積和線的長度。
如：*GET,E_VOLUME,ELEM,10,VOLU 為提取編號(hào)為10的單元的體積
*GET,a_area,AREA,50,AREA 為提取編號(hào)為50的面元的面積
*GET,l_length,LINE,100,LENG 為提取編號(hào)為100的線的長度
以上對(duì)應(yīng)的GUI操作： Utility Menu>Parameters>Get Scalar Data 如果要一次性提取多個(gè)元素的相關(guān)參數(shù)，可以用命令 *VGET, ParR, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM, KLOOP 對(duì)應(yīng)GUI操作：Utility Menu>Parameters>Get Array Data
輸入命令 alist,p 出選項(xiàng)框，選你要看的那個(gè)面積，提取選中的單元面積。
*cfopen,'area','txt',
*GET,MaxEleNum,ELEM,,NUM,MAX
*GET,MinEleNum,ELEM,,NUM,MIN
*do,i,MinEleNum,MaxEleNum,1
*if,esel(i),eq,1,then
*get,volu,elem,i,volu
*vwrite,i,volu
(f5.0,f15.12)
*end if
*enddo
*cfcols
展開 Ansys Workbench ACT插件，由窗口選中體單元，提取體積和表面積，計(jì)算幾何特征尺寸 ￥20
解決方法：
筆者這里使用一種較為笨的方法進(jìn)行選中體單元的特征尺寸提取方法：
1. 在Workbench界面，根據(jù)用戶在圖形窗口選定的網(wǎng)格體單元。由腳本程序依次提取，每個(gè)單元的角點(diǎn)數(shù)量和位置坐標(biāo)；
2. 再由Workbench中python腳本調(diào)用ADPL經(jīng)典界面，并自動(dòng)運(yùn)行特征值提取宏命令；
3. APDL宏命令會(huì)，根據(jù)Workbench選中單元體信息，依次由每個(gè)單元體的角點(diǎn)坐標(biāo)，創(chuàng)建實(shí)體單元；
4. 再將實(shí)體單元合并，最后獲得幾何體積和表面積，并輸出。
5. 由Workbench腳本，讀取APDL宏輸出的幾何信息，并顯示。
展開 ABAQUS螺栓滑移量批量提取工具，可以同時(shí)對(duì)多個(gè)工況多個(gè)面進(jìn)行提取 ￥100
適用所有螺栓滑移量提取，腳本與使用方法付費(fèi)
螺栓滑移量提取1.pptx
ABAQUS后處理之提取分層損傷面積/分層面積/基體損傷面積（ABAQUS+Photoshop） ￥28
ABAQUS后處理之提取損傷面積（ABAQUS+Photoshop聯(lián)合使用）
為了定量描述損傷程度，提取載荷造成的損傷面積變得尤為重要，下面介紹損傷面積的提取方法。
1. 去除單元網(wǎng)格，以及邊緣
2. 突出顯示損傷區(qū)域，建立損傷與未損區(qū)域色差
3. 導(dǎo)出圖片
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' 這種很長的字符串，繁瑣且容易出錯(cuò)。
在付費(fèi)內(nèi)容中，提供了一種便捷的方法，得到長字符串，提取相關(guān)的數(shù)據(jù)，并將數(shù)據(jù)保存為csv文件。
展開 ABAQUS水力壓裂模擬|XFEM和Cohesive方法-多縫、交叉縫、體積縫、轉(zhuǎn)向縫、縫間干擾
ABAQUS水力壓裂模擬|XFEM和Cohesive方法
by 星辰北極星
關(guān)鍵字：單縫、多縫、交叉縫、體積縫、轉(zhuǎn)向縫、縫間干擾、儲(chǔ)隔層
水力壓裂，對(duì)于石油工程的朋友并不陌生，它是石油開采和增產(chǎn)的重要手段；也廣泛應(yīng)用于地?zé)衢_采、地基處理等領(lǐng)域。由于畢業(yè)于石油大學(xué)，所以有很多機(jī)會(huì)接觸這方面的問題，也關(guān)注著ABAQUS在壓裂領(lǐng)域的應(yīng)用。這個(gè)專題將分享自己在水力壓裂仿真中的一些積累，希望大家喜歡。
【主要內(nèi)容】
一、課程概述
二、仿真要點(diǎn)介紹
2.1 ABAQUS水力壓裂模擬常用仿真方法
2.2 地應(yīng)力平衡分析（Geostatic）
2.3 滲流-位移耦合分析（Soils）
2.4 材料與單位制講解
2.5 特殊的輸出需求與定義
2.6 交叉裂縫處理
三、實(shí)例講解
3.1 基于Cohesive單元的二維水力壓裂模擬
3.2 基于Cohesive單元的三維水力壓裂模擬
3.3 水力裂縫與天然裂縫相交模擬-Cohesive單元法
3.4 裂縫發(fā)育地層的水力壓裂模擬-Cohesive單元法
3.5 基于XFEM的水力裂縫轉(zhuǎn)向模擬
3.6 基于XFEM的水平井多段壓裂裂縫的縫間干擾問題研究
視頻地址：https://i.xue.taobao.com/detail.htm?spm=a2174.7765247.0.0.OHNzvF&courseId=89321
【二維水力壓裂模擬（Cohesive）】
通過這個(gè)簡單的案例講述采用Cohesive單元模擬水力壓裂的基本技巧，讓大家掌握注液、停泵憋壓等基本設(shè)置，以及前后處理的一些技巧。
展開 ABAQUS實(shí)用子程序SPRINC提取主應(yīng)力
實(shí)用子程序SPRINC
在ABAQUS中使用UMAT子程序時(shí)有時(shí)會(huì)使用到最大主應(yīng)力進(jìn)行計(jì)算。通過查閱幫助文檔，ABAQUS實(shí)用子程序SPRINC可以在UMAT中計(jì)算最大主應(yīng)力和最大主應(yīng)變，SPRIND可以計(jì)算最大主應(yīng)力和最大主應(yīng)變的方向。
下面是ABAQUS幫助文檔關(guān)于實(shí)用子程序SPRINC的介紹：
SPRINC (calculate principal values)
Interface
CALL SPRINC(S,PS,LSTR,NDI,NSHR)
Variables to be provided to the utility routine
S
Stress or strain tensor.
LSTR
An identifier. LSTR=1 indicates that S contains stresses; LSTR=2 indicates that S contains strains.
NDI
Number of direct components.
NSHR
Number of shear components.
Variables returned from the utility routine
PS(I), I=1,2,3
The three principal values.
展開 
利用Python提取ABAQUS的.odb結(jié)果中的輪軌接觸信息 ￥500
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ì)算。
批量提取Abaqus指定Step-x下 Set節(jié)點(diǎn)集的反力RF（Reaction force） ￥40
<h2>摘要</h2><p>本文介紹如何使用Python腳本二次開發(fā)來批量提取ABAQUS輸出數(shù)據(jù)庫（ODB）文件中指定Step-x下的Set節(jié)點(diǎn)集的反力RF（Reaction force）。通過詳細(xì)的步驟說明、代碼示例和圖片展示，您將學(xué)會(huì)如何使用該腳本，自動(dòng)化輸出CSV文件包含（Node Label；Step Name、Increment、Step Time，RF1（X），<span style="color: rgb(25, 27, 31);">RF2（Y），RF3（Z）</span>）。</p><p><br></p><p>如果還需要提取<strong><em><u>節(jié)點(diǎn)變形量 </u></em></strong>或者<strong><em><u>節(jié)點(diǎn)坐標(biāo) </u></em></strong>可以關(guān)注：<a href="https://www.yqgqt.org.cn/post/1939915" rel="noopener noreferrer" target="_blank">利用Python腳本 批量提取Abaqus的odb文件中Set節(jié)點(diǎn)集 初始節(jié)點(diǎn)坐標(biāo)、指定Step下的變形量、變形后節(jié)點(diǎn)坐標(biāo)_CAE Abaqus提取odb-技術(shù)鄰 (jishulink.com)</a></p><h2>1. 問題描述</h2><p>在工程仿真和分析領(lǐng)域，提取ABAQUS輸出數(shù)據(jù)庫（ODB）文件中的節(jié)點(diǎn)集的<span style="color: rgb(25, 27, 31);">反力RF（Reaction force）</span>是一項(xiàng)常見任務(wù)。然而，手動(dòng)提取這些數(shù)據(jù)是一項(xiàng)繁瑣且容易出錯(cuò)的工作。
展開 Abaqus平均應(yīng)力和應(yīng)變提取 ￥80
利用python腳本對(duì)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)變平均值
Abaqus中的特征值提取
AMS特征值求解器是一種高效的，針對(duì)大規(guī)模問題的能提取大量特征值的方法，主要適用于1百萬自由度以上的模型及500階模態(tài)以上。
它包含3個(gè)求解步驟：
（1）生成子結(jié)構(gòu)；
（2）獲得特征值；
（3）從縮減的向量中獲得全部特征向量。在Abaqus中采用AMS特征值求解器的句法（提取100階頻率）
*STEP
*FREQUENCY,ELGENSOLVER=AMS
,,100
對(duì)于傳統(tǒng)的結(jié)構(gòu)，Lanczos是默認(rèn)的特征值提取方法。然而對(duì)于大規(guī)模問題，相比AMS方法，Lanczos方法效率較低。使用Lanczos方法，可以指定最大頻率的提取或提取的頻率數(shù)量，也可以指定最小頻率的提取。Lanczos特征值求解器允許計(jì)算到特征值真正的誤差限制時(shí)才終止，可以滿足正常的終止原則。對(duì)于多數(shù)問題，相對(duì)誤差為1.E-12數(shù)量級(jí)，因此Lanczos求解器的計(jì)算結(jié)果精度一般要比子空間迭代法高。而子空間迭代法的終止條件是通過判斷從這一次迭代到下一次迭代過程中特征值的相對(duì)變化來實(shí)現(xiàn)的，如果相對(duì)變化小于1.E-5則認(rèn)為已經(jīng)收斂，結(jié)束計(jì)算。
abaqus中的特征值提取.pdf
展開 