ansys提取單元的應變的案例
ANSYS如何提取能量結果（應變能，應變能密度，應變能時程）？ ￥100
<div contenteditable="false" width="100%">
<p>對于靜力分析，常提取結構的變形、<a href="https://www.yqgqt.org.cn/qa/4700" class="jsk-anchor">應力</a>、應變和約束反力等結果，相關方法可查看，而對于動力分析，常提取結構的位移、速度、加速度、反應譜等計算結果。而能量是表征物理系統做功的量度，是<a href="https://www.yqgqt.org.cn/major/Ansys" class="jsk-anchor">ANSYS</a>重要的計算結果之一。應變能(Strain Energy)是應力和應變結果計算出來的，由于變形而儲存在結構內的能量，包括由于材料塑性而產生的塑性應變能。</p>
<p>在<a href="https://www.yqgqt.org.cn/major/Ansys" class="jsk-anchor">ANSYS</a>中，/POST1中觀察整個模型在指定時刻的結果，而在/POST26中，可以觀察到指定節點在整個持時范圍的響應。本文分別從這兩個方面對ANSYS中能量的提取方法進行介紹。
展開 ABAQUS提取單元平均應力/應變 ￥10
利用python讀取odb文件（可一次讀取多個odb）生成csv（excel）文件。提供源文件，注釋詳細，可根據需要進行修改。
ANSYS Beam188提取彎矩為例介紹ANSYS定義單元表提取數據 (解決彎矩圖鋸齒狀） ￥20
在ANSYS中有些數據無法直接訪問，需要通過定義單元表完成單元的結果的訪問。下面就以Beam188單元提取彎矩為例介紹ANSYS定義單元表提取數據的詳細過程。
1. 首先需要知道在哪里定義單元表：Main Menu>General Postproc>Element Table>Define Table>add
2. 定義你想要的數據，這里以Beam188的彎矩為例
2.1 啟動ANSYS幫助菜單, 在索引框輸入Beam188然后搜索, 在單元輸出介紹找到彎矩的名稱（代號）。
2.2 回到ANSYS界面，比如要輸出Mz, 則需要在添加SMISC,3 和SMISC,16 ,如圖
3. 輸出數據：Main Menu>General Postproc>Element Table> List E T, 選擇前面定義的SMISC,3 和SMISC,16 輸出單元I和J節點的Mz數值，如圖
4. 顯示彎矩云圖：Main Menu>General Postproc>Plot Results>Contour Plot>Line Elem Res, 這里要注意要在LabI 選SMISC,3 LabJ 選SMSCI,16。
輸出彎矩到這就結束了，小編突然發現，輸出的彎矩值在每個單元的I和J處是一樣的（Beam188為2節點單元），彎矩圖也就成了鋸齒形，于是去問了度娘一波，各路盆友給出解決方法，然而并沒有起作用的，于是乎我又想起來了“幫助文檔大法”，于是認認真真將Beam188的幫助文檔閱讀了一遍，功夫不負有心人，最終。。。
展開 ANSYS瞬態分析全時程結構響應最大值的提取方法（變形、應力、應變、能量） ￥100
</p><p>同樣的方法，可以提取全時程最大的位移、應力、應變、能量等結果。</p><p><img src="https://img.jishulink.com/upload/202302/623025b5c0d646b9973cd2adc6c6037f.png" alt="1.png"></p><p>收費內容為相關命令流。</p>
ANSYS中整體、單元剛度和質量矩陣的提取
指定輸出單元矩陣
/SOLU
SOLVE
finish
/OUTPUT, TERM ! 將輸出信息送到output windows中
! 這時用編輯器打開cp.out文件，可以看到按單元寫出的質量、剛度等矩
陣
ANSYS中整體、單元剛度和質量矩陣的提取.rar
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
展開 『分享』ANSYS中整體、單元剛度和質量矩陣的提取
、整體剛度和質量矩陣的提取。
該功能需要進行二次開發，由ansys形成
的二進制文件.full提取整體剛度和質量
矩陣。
基于ansys的一個用戶開發程序例子編
了一個程序（附件中）。
開發環境：compaq fortran 6.5
運行環境：win2000。
一個主文件：self.for，
另一個文件matrixout.f90用于矩陣輸出
binlib.lib為ansys提供的庫文件，將其
引入項目中（也可直接扔進debug目錄）
.full文件由子空間迭代模態分析獲得。
運行編譯后的可執行文件.exe
2、單元剛度和質量矩陣的提取。
/DEBUG命令。詳細說明可由以下轉載文章
finish
/clear
PI=3.1415926
w1=3
w2=10
w3=6
w4=1.2
r=.8
t=0.08
/PREP7
!*
ET,1,SHELL63
R,1,t
ET,2,MASS21
R,2,500,500,500,2000,2000,2000,
!*
UIMP,1,EX, , ,2e11
UIMP,1,NUXY, , ,0.3,
UIMP,1,DAMP, , ,0.2,
UIMP,1,DENS, , ,7800,
BLC4,0,0,w2,w1
ESIZE,1.5,0,
AMESH,all
NSEL,S,LOC,X,0.0
D,all, , , , , ,ALL, , , , ,
allsel,all
SFA,all,1,PRES,12
FINISH
/OUTPUT,cp,out,, !
展開 Ansys Workbench ACT插件，由窗口選中體單元，提取體積和表面積，計算幾何特征尺寸 ￥20
Ansys Workbench ACT插件，由窗口選中體單元，提取體積和表面積，計算幾何特征尺寸
問題：
在FKM關于結構疲勞評估計算方法中指出：零部件特征尺寸，影響疲勞結果評估。原因是材料的應力壽命曲線是由標準試樣進行試驗測試獲得的。當零部件的特征尺寸與測試樣件不一致時，需要考慮零部件的特征尺寸這一因素。（一般而言，當零部件的尺寸大于材料標準測試樣件時，零部件的表面或內部缺陷發生的概率會增加，從而導致零部件尺寸越大，疲勞壽命越低）
對與規則幾何形狀的零部件，有相應的經典公式提供特征尺寸的計算；例如圓形細長桿的特征尺寸是直徑；薄板零部件的特征尺寸是板厚等；但是實際工作中的零部件幾何形狀千差萬別，沒有統一的經典公式可以提供特征尺寸的計算；在FKM手冊中給出了一個通用公式，用于估計零部件疲勞危險區域的局部特征尺寸；
FKM關于循環載荷的疲勞評估中，提及可以使用循環載荷下的有限元應力結果進行疲勞損傷估計。此時，除了需要由應力結果估計危險疲勞區域，提取危險點的應力結果外，還需要給出危險疲勞區域的特征尺寸。在Ansys Workbench中，用戶可以方便的查看應力結果云圖，從而大體評估出危險疲勞區域。并且用戶可以通過選取高應力區域的單元體，再通過特征尺寸一般計算公式，來估計高應力區域的特征尺寸，進行進行合理的FKM疲勞評估。
但是，Ansys Workbench中，當用戶選中了某個/某些體單元后，在選擇信息欄中并不能直接給出單元體積和表面的有效信息輸出。并且通過查詢資料，即使在APDL經典界面中對與體單元也是僅僅只能輸出體積（沒有體單元表面的輸出）；并且對與FKM特征尺寸的一般計算公式中，關于表面積A，也并不是指每個體單元所有面的表面積的總和。
展開 【Ls-dyna】Hypermesh&Ls-dyna聯合仿真時如何設置輸出單元應變和查看單元應變？
通常，使用Ls-dyna進行瞬態分析時，默認的輸出控制中是不包含單元的應變值的。但是，如果關心單元上的應變，或者需要對單元應變有更全面的了解，那么就需要在建模和分析中進行輸出控制，使計算結果中包含單元的應變。那么，問題來了，如何在Hypermesh中如何設置可以輸出單元的應變？如何在LS-POST中顯示單元的應變云圖？ 接下來，依次解答。
如何在Hypermesh中如何設置可以輸出單元的應變？
這里需要首先說明的是，本人使用Ls-dyna以來，一直使用的前處理不是Ls-prepost前處理，而是前處理軟件Hypermesh，后處理則使用Ls-prepost。所以，這里只說明Hypermesh&Ls-dyna聯合仿真時如何在Hypermesh中設置控制應變輸出。
前處理軟件工作界面
后處理軟件工作界面
首先啟動Hypermesh，點擊“user Profiles”，選擇Ls-dyna，其他默認，進入Ls-dyna分析模塊。
然后，在軟件的面板區域選擇“analysis”，并點擊“control card”。
點擊面板區域的Next，一直到出現“database-extent-binary”，點擊該按鈕
在關鍵字*Database-Extent-Binary的定義中，將第一行第四個參數【STRFLAG】的值設置為1，表示在二進制結果文件d3plot中輸出單元的應變。
在Hypermesh中按照上述方式就可以在d3plot中輸出單元的應變，由于d3plot是二進制文件，所以只能借助后處理LS-Prepost查看應變結果。
那么，如何在LS-POST查看應變，顯示應變云圖呢？
打開LS-PrePost，點擊【Fcomp】，選擇【strain】，就可以觀察計算后的應變云圖。
展開 Abaqus平均應力和應變提取 ￥80
利用python腳本對ODB文件中單元集里所有積分點的應力及應變進行自動提取并計算平均值
能夠得到每一幀的應力和應變平均值，并保存到CSV文件中
所得到的應力包括S11,S22,S33,S12,S13,S23以及Mises七個應力平均值，以及E11,E22,E33,E12,E13,E23六個應變平均值
Abaqus插件——平均應力應變提取 ￥60
通過該插件可實現:
1）提取所有幀的任意單元集合的平均應力（事先定義單元集合，如圖中的SET-1）
2）提取所有幀的任意區域的x、y、z方向的平均應變（事先定義節點集合，如圖中的SET-2）
3）將以上數據保存至excel文件（excel文件名為odb文件名稱+_Stress_Strain.csv）
*************************注意事項******************************
1、插件使用過程中，如有任何問題請發郵件至shenz1hao@126.com
2、插件僅做學習交流使用，尊重原創者，切勿以營利目的傳播
*****************************************************************
********************插件安裝及使用*******************************
1、電腦路徑下輸入 %homepath%\abaqus_plugins并回車
2、將Stress-Strain文件夾解壓至當前目錄下
3、打開abaqus，菜單欄中點擊plug-ins，里面找出Stress-Strain
4、輸入相應參數（hx、hy、hz表示x、y、z方向模型長度，當以上三參數取1時輸出的為該方向位移）
*****************************************************************
展開 
Python提取場輸出結果計算平均應力應變
###########################################################################
寫在前面，最近回看了一些自己之前編寫的和這個py代碼，發現自己寫的可能和公式表達的有偏差，在我的測試過程中，提取的是單元積分點處的應力應變值，由于只有一個積分點（這個積分點在質心，與centroid提取得到的結果一樣），因此自然的把這個值當成了單元的平均應力或平均應變，然后進行計算，但是我現在覺得這并不是獲取單元平均應力的方式，也就是代碼并沒有實現所謂的提取平均應力應變的功能，希望有大神可以指點迷津。
如果是有多個積分點的話，是不是應該對每個積分點權重進行積分，加權平均這樣得到單元的平均應力，然后乘單元體積，將所有單元的值求和再除模型的總體積，就得到整個RVE模型的平均應力。
###########################################################################
最近寫了一個簡單的python讀取abaqus結果中的場輸出數據，想通過均勻化計算方法來計算所定義集合的平均應力應變曲線，之前是手動提取了各個數據導出，然后用excel、matlab處理，但是很慢，而且很費勁，于是就想著用Python來處理結果。
有需要的同學可以下載附件文件，打開abaqus，file→run script，選擇腳本文件即可運行。
average.zip
均勻化計算方法：
參考文獻：馬思鳴. 精沖用碳鋼微觀組織對宏觀力學性能及精沖性能影響研究[D]. 上海：上海交通大學，2017.
我在網上找了挺久挺多的關于Python提取場輸出結果的，
主要對以下幾篇帖子進行了參考：
http://forum.simwe.com/forum.php?
展開 有限元2D單元妙用 平面應力與平面應變 廣義平面應變 硬干涉 ￥10
平面應力單元還可以跟軸對稱單元結合，模擬出變厚度模型。比如對葉盤的分析。需要注意的是，在ANSYS里面，當我們將平面應力和軸對稱單元結合的時候，平面應力單元的厚度應該設置為所有圓周分布葉片厚度的總和。如下圖。
平面應變單元：
代表性體積單元根據單元體積應力應變加權平均 ￥20
現如今，越來越多的人開始對復合材料性能進行研究，如何通過<a href="/major/abaqus">abaqus提取代表性體積單元是非常重要的，我提供了一種可以根據單元體積進行應力應變平均的代碼，希望對大家有用。
基于Python腳本提取復合材料應力應變均勻化有效彈性模量 ￥75
通過均勻化方法使用Python腳本從細觀尺度得到材料宏觀的等效模量
https://www.bilibili.com/video/BV1sb4y1273j
具體腳本使用方法及理論介紹詳見視頻教程
