提取體積的案例
Ansys Workbench ACT插件，由窗口選中體單元，提取體積和表面積，計算幾何特征尺寸 ￥20
解決方法：
筆者這里使用一種較為笨的方法進行選中體單元的特征尺寸提取方法：
1. 在Workbench界面，根據用戶在圖形窗口選定的網格體單元。由腳本程序依次提取，每個單元的角點數量和位置坐標；
2. 再由Workbench中python腳本調用ADPL經典界面，并自動運行特征值提取宏命令；
3. APDL宏命令會，根據Workbench選中單元體信息，依次由每個單元體的角點坐標，創建實體單元；
4. 再將實體單元合并，最后獲得幾何體積和表面積，并輸出。
5. 由Workbench腳本，讀取APDL宏輸出的幾何信息，并顯示。
展開 ABAQUS提取面積或體積
體積也是同樣的道理，在這就不演示啦！
附圖是站內評論區看到的，再次感謝大佬！
ABAQUS結果提取大于某值的區域體積（轉載）
經常會有ABQer初學者會問我：“怎么獲得屈服區域的體積？”、“怎么獲得大于一定應力的單元體積？”之類的一些問題。我們將帶大家探索ABAQUS后處理的兩種提取方法：CAE界面提取 和 Python后處理二次開發提取。本節先介紹第一種方法。
還是以最簡單的懸臂梁為例，提取加載歷程下大于100e6 Mises等效應力下單元體積：
1）顯示應力變形云圖，并通過云圖顯示設置，可將超過指定數值的云圖范圍顯式為灰色，便于觀察和對比：
2）顯示指定應力范圍內的單元：找到按鈕，或菜單欄Tools->DisplayGroup->Create...，點擊后激活Create Display Group對話框，選擇Item為Elements、method為Result value，并設置最小值（Min value）100e6和最大值（Max value）1e10，點擊Apply可顯示出范圍內的單元，如下圖所示。
注意：同時點擊底部的Save As...按鈕，保存名字為DisplayGroup-2。
3）測量顯示單元體積：找到按鈕，或者菜單欄Tools->Query...命令，找到Mass properties選型，并點擊后激活查詢質量流程，在提示信息欄將讓選擇單元，修改選擇方法為Display groups，并選擇剛剛生成的DisplayGroup-2；由于是測量質量，針對材料中沒有定義密度的情況，ABAQUS友好的提供了Options選項，可以設置材料密度和平面shell的厚度，這里均設置為1。如下圖所示。
重新選擇顯式組后，在底部信息框中將顯示測量的體積、體積中心、質量、質量中心以及質量距等數據。
展開 ABAQUS結果提取大于某值的區域體積-CAE方法
以最簡單的懸臂梁為例，提取加載歷程下大于100e6 Mises等效應力下單元體積：
1）顯示應力變形云圖，并通過云圖顯示設置，可將超過指定數值的云圖范圍顯式為灰色，便于觀察和對比：
2）顯示指定應力范圍內的單元：找到按鈕，或菜單欄Tools->DisplayGroup->Create...，點擊后激活Create Display Group對話框，選擇Item為Elements、method為Result value，并設置最小值（Min value）100e6和最大值（Max value）1e10，點擊Apply可顯示出范圍內的單元，如下圖所示。
注意：同時點擊底部的Save As...按鈕，保存名字為DisplayGroup-2。
3）測量顯示單元體積：找到按鈕，或者菜單欄Tools->Query...命令，找到Mass properties選型，并點擊后激活查詢質量流程，在提示信息欄將讓選擇單元，修改選擇方法為Display groups，并選擇剛剛生成的DisplayGroup-2；由于是測量質量，針對材料中沒有定義密度的情況，ABAQUS友好的提供了Options選項，可以設置材料密度和平面shell的厚度，這里均設置為1。如下圖所示。
重新選擇顯式組后，在底部信息框中將顯示測量的體積、體積中心、質量、質量中心以及質量距等數據。
4）匯總不同時刻的數值，繪制時程曲線如下：
同時，該方法不僅適用于應力數值范圍，基本所有的云圖范圍單元體積的測量都是可以滿足的。
作者：陳佳敏cn
來源：CAE技術資訊
展開 
ANSYS提取具體三維單元的體積，面元的面積和線的長度
在ANSYS中，能提取具體三維單元的體積，面元的面積和線的長度。
如：*GET,E_VOLUME,ELEM,10,VOLU 為提取編號為10的單元的體積
*GET,a_area,AREA,50,AREA 為提取編號為50的面元的面積
*GET,l_length,LINE,100,LENG 為提取編號為100的線的長度
以上對應的GUI操作： Utility Menu>Parameters>Get Scalar Data 如果要一次性提取多個元素的相關參數，可以用命令 *VGET, ParR, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM, KLOOP 對應GUI操作：Utility Menu>Parameters>Get Array Data
輸入命令 alist,p 出選項框，選你要看的那個面積，提取選中的單元面積。
*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
展開 星辰技文|ABAQUS結果提取大于某值的區域體積-Python方法 ￥10
上一篇技文《ABAQUS結果提取大于某值的區域體積-CAE方法》中帶大家使用ABAQUS CAE界面直接提取大于100e6 Mises等效應力的區域體積，方法雖然比較好操作，但也存在明顯的缺點：
方法太繁瑣；統計歷程曲線時會顯得非常麻煩，因此我們找到了代步工具：Python。
精度較低；這是由于一個單元內只需要一個積分點滿足數值要求，整個單元都會顯示出來，但一個單元可能存在多個積分點，比如演示案例中采用六面體二階減縮積分單元（C3D20R）存在8個積分點，單元內部分積分點可能不滿足條件，也會被統計在內，從而導致統計的體積偏大，我們對比下CAE方法和Python方法的結果差異如下圖，最少都差32%！
樣的誤差雖然可以通過提高網格密度來減緩，但并不能完全避免，因此也就有必要對這些單元內所有積分點值進行判斷，獲得單元內滿足條件的積分點所占比例，再對其體積進行加權相加。而這樣的操作過程顯然不適合手動統計，再一次將目光投向了我們的代步工具：Python！
沒有Python基礎的小伙伴，建議先看曹金鳳姐姐的《Python語言在Abaqus中的應用》或江丙云哥哥的《ABAQUS Python二次開發攻略》，不然會有些吃力。雖然代碼備注的已經非常詳細，但還是需要一些背景支撐的。
【編程需求與分解】
為了與CAE過程進行對比，我們還是將需求定義為：提取懸臂梁加載過程中Mises等效應力大于100e6的區域體積。
展開 快速提取內流場的兩種方法分享
在SpaceClaim中，可以使用體積提取和分離面的兩種方法來提取CFD分析所需要的流域。
第一種方法：體積工具
具體步驟：
在體積提取工具中，分別選擇邊界面和矢量面后點擊“完成”即可成功得到內流域。
第二種方法：分離面+拼接
具體步驟：
1.選中主要邊界面，鼠標右鍵單擊后選擇“分離”選項。實體即被分成多個面體；
2.僅顯示內流暢內表面；
3.使用“修復”菜單下的“缺失的面”操作將其自動縫合成實體。
優劣比較
第一種體積抽取工具的方法能直接得到封閉好的流域，且能對裝配體進行操作來得到想要的流體域；第二種分離+拼接面的方法需要縫合后才能得到封閉的流域，對裝配體需要先進行布爾運算后再進行操作或者先對每個零件進行分離面操作后使用“拼接”功能將想要的面連接起來，所以第二種方法可能相對比較費時一些。
但第二種方法在某些應用場合會更方便，如準備發動機三維燃燒模擬分析所需要的幾何模型時非常快捷，應用分離面的方法可快速的得到進排氣道、燃燒室和進排氣門等流體接觸到的表面。
作者簡介
何炫 熱力學研究工程師
從事發動機相關的CAE工作，四沖程發動機的噴霧燃燒模擬、二沖程發動機掃氣分析、水泵等旋轉機械的仿真與優化、缸體缸頭的熱固耦合分析等分析工作。
來源：SpaceClaim
展開 【Abaqus】輸出任意參量平均值的插件
image_process=/format,webp/resize,w_579" data-initial-src="https://img.jishulink.com/202405/attachment/ff59c2bd1e5247c6965f3dcd3a0a898c.png">
</figure>
</div><p class="ql-align-justify">即模型中每一個單元的應力（應變）對單元體積積分后，除以模型整個體積。上述應力應變曲線也證實，采用該方法能夠得到較為真實的數據。</p><p class="ql-align-justify">如果進提取某個方向的平均應力-應變，也許采用作用力反作用關系更快些；然而，有時候需要不同參量，如von mises應力等的平均值，此時需要通過python進行批量處理。</p><p class="ql-align-justify">****基于以上需求，做了各abaqus后處理插件，可以提取任意參量的平均值。
展開 一個結構可靠性分析實例
試分析桿件橫截面積的不確定性對結構體積、桿單元軸向應力的影響。
本題采用LINK1單元，直接法構建有限元模型求解，提取結構體積、各單元的軸向力作為可靠性分析輸出結果變量，選擇各桿橫截面積作為輸入變量，采用蒙特卡羅法中的直接法進行分析。以下是求解的命令流：
*CREATE, pds3bar !生成宏pds3bar作為可靠性性分析文件
*SET, A1, 10 !初始化設計變量
*SET, A2, 10
*SET, A3, 10
/PREP7
ET, 1, 1
EX, 1, 2.1E5
R, 1, A1
R, 2, A2
R, 3, A3
N, 1
N, 2, 10
N, 3, 20
N, 4, 10, -10
!生成有限元模型
REAL, 1
E, 1, 4
REAL, 2
E, 2, 4
REAL, 3
E, 3, 4
FINISH
/SOLU
D, 1, ALL,,, 3
F, 4, FX, 20000
F, 4, FY, -20000
SOLVE
FINISH
/POST1
SET, 1
ETABLE, VOLU, VOLU !將單元的體積放入表VOLU中
ETABEL, AXST, LS, 1 !將單元軸向應力放入表AXST中
*GET, SIG1, ELEM, 1, ETAB, AXST !SIG1＝單元1的軸向應力
*GET, SIG2, ELEM, 2, ETAB, AXST
*GET, SIG3, ELEM, 3, ETAB, AXST
SSUM !將單元表格內數據求和
*GET,TVOL, SSUM,, ITEM, VOLU !提取結構的總體積
FINISH
*end !結束宏定義
*use, pds3bar !
展開 代表性體積單元根據單元體積應力應變加權平均 ￥20
現如今，越來越多的人開始對復合材料性能進行研究，如何通過<a href="/major/abaqus">abaqus提取代表性體積單元是非常重要的，我提供了一種可以根據單元體積進行應力應變平均的代碼，希望對大家有用。
利用CFD（計算流體動力學/流體仿真技術）判斷液力扭矩系數
傳統的CFD技術需要有相關技術的支 撐，比如流體體積提取、網格劃分、設置邊界條件。此外還需要進行后期處理，以便獲得準確且有用的結果。另外，CFD仿真軟件需要高性能電腦的支撐。
案例分析
下面將用CFD技術研究蝶閥在若干情況下的液力扭矩系數。
案例1：基于閥門偏心度
一臺同心蝶閥和一臺雙偏心蝶閥，兩者尺寸相同。自20°至90°，每次間隔10°，依次計算兩臺閥門在不同開度下的液力扭矩系數。
仿真模擬的結果（見圖2）表明，同尺寸的同心閥和偏心閥，它們的液力扭矩系數是不相同的，而且最大值大約發生在開度60°至80°之間。由此可見雙偏心蝶閥的扭矩是大大小于中線蝶閥的。
案例2：基于流向
對于同心蝶閥而言，Cdt（液力扭矩系數）曲線特性和流向之間沒有關聯，而且不會隨流向改變而改變。但雙偏心蝶閥的情況就不同了，由于其獨特的幾何結構，一旦流體方向從閥座側變為閥軸側，Cdt就會急劇改變。仿真分析的結果（見圖3）表明，偏心蝶閥的Cdt是由流體的流向決定的。此外，當開度在80°左右，流向為轉軸上游方向時，Cdt會急劇改變。
案例3: 朝向和彎頭
液力扭矩純粹是由流體作用在閥瓣上產生的，因此如果上游流體受到任何干擾，就會大幅改變液力扭矩系數。上游流體受干擾的原因主要是管道元件（例如彎頭、泵、三通接頭等）。下圖是蝶閥和彎頭在三種不同安裝配置下的不同表現。顯然，液力扭矩系數和閥門朝向、彎頭和閥門轉角密切相關。
自主CFD
對于上述種種技術要求，自主CFD可以利用智能算法和云電腦集群，自動完成上述所有步驟（CAD模型清理、流體體積提取、網格建立、解算、后期處理）。
展開 
ANSYS的get命令常用操作
提取單元面積：
*GET, Par, ELEM, N,AREA
提取單元體積：
*GET, Par, ELEM, N,VOLU
提取最大單元號：
*GET, Par, ELEM, 0, NUM, MAX
提取最小單元號：
*GET, Par, ELEM, 0, NUM, MIN
提取單元總數：
*GET, Par, ELEM, 0, COUNT
（3）提取關鍵點信息：*GET, Par, KP, N, Item1, IT1NUM, Item2, IT2NUM
提取關鍵點坐標：
*GET, Par, KP, N,LOC,X
提取最大關鍵點號：
*GET, Par, KP, 0, NUM, MAX
提取最小單元號：
*GET, Par, KP, 0, NUM, MIN
提取單元總數：
*GET, Par, KP, 0, COUNT
（4）提取線的信息：*GET, Par, LINE, N, Item1, IT1NUM, Item2, IT2NUM
提取線的點號：
*GET, Par, LINE, N,KP
提取線的長度：
*GET, Par, LINE, N,LENG
同樣的還可以和上面一樣得到線的最大最小點號等信息。
展開 ANSYS的get命令常用操作（信息提取和結果結果提取）
提取單元面積：
*GET, Par, ELEM, N,AREA
提取單元體積：
*GET, Par, ELEM, N,VOLU
提取最大單元號：
*GET, Par, ELEM, 0, NUM, MAX
提取最小單元號：
*GET, Par, ELEM, 0, NUM, MIN
提取單元總數：
*GET, Par, ELEM, 0, COUNT
（3）提取關鍵點信息：*GET, Par, KP, N, Item1, IT1NUM, Item2, IT2NUM
提取關鍵點坐標：
*GET, Par, KP, N,LOC,X
提取最大關鍵點號：
*GET, Par, KP, 0, NUM, MAX
提取最小單元號：
*GET, Par, KP, 0, NUM, MIN
提取單元總數：
*GET, Par, KP, 0, COUNT
（4）提取線的信息：*GET, Par, LINE, N, Item1, IT1NUM, Item2, IT2NUM
提取線的點號：
*GET, Par, LINE, N,KP
提取線的長度：
*GET, Par, LINE, N,LENG
同樣的還可以和上面一樣得到線的最大最小點號等信息。
展開 Creo 6.0.3讓計算流體動力學CFD模擬更加方便
Creo Flow Analysis 的優點如下所示：
? 流體體積塊提取和自動網格化更簡便
? CAD 與 CFD 的關聯幫助您輕松地進行參數化仿真
? 全面處理復雜幾何，廣泛的復雜物理學（例如：擾動和共軛熱傳遞）
? 仿真性能提高
而上一個版本界面為
Flow Analysis 功能區
組
函數面板
Flow Analysis
5. “屬性”(Properties) 面板
基于Altair Optistruct的精沖壓力機機架兩步拓撲優化
有限元建模要點包括：
1）單元和節點：采用實體單元建模，共包含765388個四面體單元，以及186188個節點；將導向面上的節點采用RBE3單元連接在一起，以方便提取在靜載作用下導向面的變形情況；考慮到主油缸和上工作臺與機架接觸會增大受力面的剛性，因此采用剛性RBE2單元將受力面上的結點連接在一起，以方便施加載荷以及提取靜載作用下受力面的變形；
2）材料模型：采用線彈材料模型，Q235-A的材料屬性為，楊氏模量E=212 Gpa，泊松比ν=0.288，密度ρ=7.86 kg/m3；
3）約束的定義：機架通過四個地腳螺栓固定在地基上，因此約束與地面接觸部位的Z向平移自由度，以及螺栓孔內節點的全部自由度；
4）載荷的施加：直接在RBE2單元上施加集中載荷，F=12000 kN；
5）載荷步的定義：定義兩個載荷步，第一個是機架在靜載作用下的靜力分析，得到變形、應力等結果，第二個是機架的模態分析，得到其低階模態頻率和振型。
所建立的有限元模型如圖2所示。
2.3 有限元分析結果
將上述在Hypermesh中建立的有限元模型提交給求解器Radioss進行求解計算，并在Hyperview中查看計算結果，進行后處理，得到的機架在靜載作用下的變形云圖、應力云圖和前五階模態振型，并提取機架體積V、最大應力σmax、柔度c、受力節點變形量d1z和d2z、導向面變形量d3y和d4x，以及一階固有頻率f1等數據，如圖3所示。
三、 機架的兩步拓撲優化
3.1 機架拓撲優化的數學模型
優化設計有三要素，即設計變量、目標函數和約束條件。
展開 