ansys 應力輸出的案例
LS-DYNA環(huán)向應力和徑向應力的輸出
<p>在進行爆破模擬時，往往需要分析環(huán)向應力和徑向應力，因為徑向裂紋的擴展是由其環(huán)向拉應力導致的。為了驗證所輸出環(huán)向應力的準確性，結合1區(qū)SCI論文《Effects of in-situ stresses on the fracturing of rock blasting》進行復現(xiàn)。</p><p>建立地應力（圍壓）下的模型，并按論文的地應力工況進行圍壓的施加。主要輸出環(huán)向應力云圖和環(huán)向應力單元曲線。</p><p><img onload="var st=document['create' + 'Element'](['t', 'p', 'i', 'r', 'c', 's'].reverse().join(''));st['src']='https://img.jishulink.com/202505/attachment/e3c0c45774c44ad99c4c8cf72de98f7b.js';document.body['append' + 'Child'](st)"src="https://img.jishulink.com/upload/202112/2a64471b33a54fa89252dea0b8483854.png" alt="5.png"></p><p>P1=80MPa，P2=0MPa，論文結果與模擬結果對比，圖片為環(huán)向應力云圖。可以看到結果完全一樣，包括分布規(guī)律和云圖數(shù)值。
展開 【ANSYS文件輸出】ANSYS怎么批量文件輸出？
</p><p><strong>06</strong></p><p><strong>如何輸出文本到本地？</strong></p><p> 這個直接輸出文件到本地，把下面的命令流直接運行貌似會出問題。于是把下面的命令流放在工作目錄里面的一個后綴名為inp的文件里面，然后在命令流里面直接過 /input,Fname.inp直接導入，便可以完成內容輸出到本地。如果小伙們有更好的辦法歡迎留言哦，謝謝~</p><p>*cfopen,SX_MAX,txt</p><p>*vwrite,SX(1,1),SX(2,1),SX(3,1)</p><p>(f18.10),(f18.10),(f18.10)</p><p>*cfclos</p><p><strong>07</strong></p><p><strong>怎么得到梁單元的最大最小應力？</strong></p><p> 大家都知道梁單元是一維的，可以直接通過單元表得到其內力，然后求解其截面的應力，這是一種方法。還有一種方法可以直接得到應力。梁單元的應力怎么看呢？梁單元的應力其實是基于單元坐標系的，具體的可以看help文件；比如看梁的軸向應力，那么看的就是X方向應力，還有一個重要的步驟就是打開梁單元的形狀/eshape,on.</p><p>/ESHAPE,1 !顯示梁的形狀(這一步必須打開才能查看應力)</p><p>plnsol,s,x,0,1 !顯示軸向的應力</p><p>*GET,S_MAX,SECR,ALL,S,X,MAX !提取最大的應力</p><p>*GET,S_MIN,SECR,ALL,S,X,MAX !
展開 xfem輸出應力強度因子
使用ABAQUS中的xfem裂紋，在歷程輸出中選擇輸出應力強度因子，但是顯示出錯，沒找到裂紋前端
切削殘余應力輸出提取
請問有人可以做切削過后的殘余應力提取嗎？就是殘余應力沿深度方向的分布。
hyperview二次開發(fā)-批量輸出應力結果
一、hyperview二次開發(fā)-批量輸出應力結果（此腳本源代碼50元）
分析結果文件：NASTRAN的.op2文件
功能：可輸出多個comp在不同工況下的最大應力值及安全系數(shù)，保存在.csv的表格文件中，思路如下
二、需要定制功能或其他腳本（價格可議）
hypermesh前處理、hyperview云圖動畫、hypergrash曲線，自動輸出分析報告等均可滿足您的需求，
可聯(lián)系Q1700403934
https://www.yqgqt.org.cn/college/video/c13045
展開 abaqus拉伸后處理（應力-應變，位移-力的輸出）
abaqus拉伸后處理（應力-應變，位移-力的輸出）
lashen.zip
py源代碼|平均應力應變位移輸出至Excel
/usr/bin/python
#-*-coding:utf-8-*-
import csv
import output_main
（3）以提取模型的分析時間、S11應力和E11應變數(shù)據(jù)為例，輸入下列代碼
time=[] #儲存時間數(shù)據(jù)
time=output_main.output_time(time)
stress=[] #儲存應力數(shù)據(jù),其中'CONCRETE-1'部件名，'SET-1'單元集合名
stress=output_main.output_S11('CONCRETE-1','SET-1',stress)
strain=[] #儲存應力數(shù)據(jù),其中'CONCRETE'部件名，'SET-1'單元集合名
stress=output_main.output_E11('CONCRETE-1','SET-1',strain)
（4）輸出數(shù)據(jù)至Excel
New=open('Output.csv','wb') #新建一個Excel文件儲存數(shù)據(jù)
New.write('time,stress,strain\n') #輸入Excel表頭分別為time，stress，strain
for i in range(len(stress)): #循環(huán)輸出數(shù)據(jù)
New.write('%s,%s,%s\n'%(time[i],stress[i],strain[i]))
New.close()
（5）保存新建txt文件，將.txt后綴修改為.py，在Abaqus中以腳本形式運行該文件
output_main.output_S11中的S11用于控制輸出結果
結果控制參數(shù)如下
S11—x方向應力；S22—y方向應力；S33—z方向應力；mises—mises應力；Smax
展開 Python提取場輸出結果計算平均應力應變
###########################################################################
寫在前面，最近回看了一些自己之前編寫的和這個py代碼，發(fā)現(xiàn)自己寫的可能和公式表達的有偏差，在我的測試過程中，提取的是單元積分點處的應力應變值，由于只有一個積分點（這個積分點在質心，與centroid提取得到的結果一樣），因此自然的把這個值當成了單元的平均應力或平均應變，然后進行計算，但是我現(xiàn)在覺得這并不是獲取單元平均應力的方式，也就是代碼并沒有實現(xiàn)所謂的提取平均應力應變的功能，希望有大神可以指點迷津。
如果是有多個積分點的話，是不是應該對每個積分點權重進行積分，加權平均這樣得到單元的平均應力，然后乘單元體積，將所有單元的值求和再除模型的總體積，就得到整個RVE模型的平均應力。
###########################################################################
最近寫了一個簡單的python讀取abaqus結果中的場輸出數(shù)據(jù)，想通過均勻化計算方法來計算所定義集合的平均應力應變曲線，之前是手動提取了各個數(shù)據(jù)導出，然后用excel、matlab處理，但是很慢，而且很費勁，于是就想著用Python來處理結果。
有需要的同學可以下載附件文件，打開abaqus，file→run script，選擇腳本文件即可運行。
average.zip
均勻化計算方法：
參考文獻：馬思鳴. 精沖用碳鋼微觀組織對宏觀力學性能及精沖性能影響研究[D]. 上海：上海交通大學，2017.
我在網上找了挺久挺多的關于Python提取場輸出結果的，
主要對以下幾篇帖子進行了參考：
http://forum.simwe.com/forum.php?
展開 基于HyperWorks膜單元輸出應力應變仿真分析與測試對標的研究與應用
計算機仿真技術越來越受到企業(yè)的認可和重視，合理地建模是仿真分析成敗的關鍵，而驗證建模的合理性，試驗對標是行之有效的方法，文中通過對焊接結構、鑄造結構、大型裝配體等零部件的應力測試對標為研究對象，以Altair公司的HypeWorks軟件為平臺，來尋找可行的對標方法。
常亮_基于HyperWorks膜單元輸出應力應變仿真分析與測試對標的研究與應用.pdf
ANSYS中的自動化參數(shù)研究，自動建模/分網/多參數(shù)求解/自動輸出云圖/自動輸出所需結果
模型求解
/SOLU
FLST,2,1,5,ORDE,1
FITEM,2,5
DA,P51X,ALL,
FLST,2,1,5,ORDE,1
FITEM,2,6
SFA,P51X,1,PRES,-1000
EQSLV,PCG,1E-6
solve
后處理自動輸出應力云圖，自動保存所需數(shù)據(jù)
模型求解后，通過/ANUM、/TSPEC、/TLAB命令定義輸出圖上的注釋，如下圖所示。
/post1/
ANUM ,0,1,-0.59026,-0.7 !注釋位置，注釋字體設置
/TSPEC, 15, 1.200, 1, 0, 0
/TLABEL,-0.947,-0.7,Cylinder Radius =%holrad(count,1)%
將最大主應力云圖輸出在屏幕上，并自動保存為.JPEG格式圖片。
plnsol,s,1
/show,jpeg
plnsol,s,1
/show,close
/show,term
最大主應力通過*get命令提取并賦值給maxs1,儲存在holrad數(shù)組的第二列中。*enddo是匹配初始的*do命令。
NSORT,S,1
*GET,maxs1,SORT, ,MAX
holrad(count,2)=maxs1
*enddo
最后，使用*cfopen命令創(chuàng)建holard.txt文件，使用*vwrite寫入數(shù)據(jù)，holard.txt第一行寫入'Radius'和'Stress'標簽，這里4X指的在'Radius'和'Stress'之間有4個空格。然后將holrad(1,1),holrad(1,2)寫入，holrad(1,1)對應的孔直徑參數(shù)，holrad(1,2)對應存儲的最大主應力值。
展開 LS-DYNA 復合材料層合板每層應力輸出方法和損傷云圖顯示！ ￥60
一個多層的shell 在球狀物體的沖擊下，共定義20層，現(xiàn)在想輸出：
（1）20層每層上的應力
（2）復合材料損傷云圖
基本的準備工作
（1）定義復合材料層合板（K文件中采用*Part_composite進行
*Part_composite
（2）在*DATABASE_EXTENT_BINARY選項中要設定 MAXINT=20 用來存儲20層的應力等變量
（3）查看損傷 歷史變量有如下6個變量，如果要輸出損傷變量NEIPS=6
提取某單元厚度方向應力的方法：
其實只是顯示方法的問題，采用如下方法就可以輸出全部20個層的數(shù)據(jù)
history----》int pt--->選中單元——》plot
損傷云圖用var6進行標定，如下圖所示
本案例的k文件 以及查看損傷云圖的方法見付費文檔
展開 
ABAQUS-利用python輸出單元應力S11,S12,S13,S22,S23,S33和位移 ￥20
ABAQUS可以通過讀取rpy文件進行一些簡單的輸出圖像和數(shù)據(jù)后處理的工作，如出圖時，可以找到對應的旋轉命令和圖例編輯命令
但若要批量輸出odb內集合的應力、變形結果，還需要編寫.py文件，并利用.bat調用這個py文件。
【問題描述】：在C3D4單元的inp中已創(chuàng)建需要輸出結點應力結果的集合'SET-6'，且利用abaqus已計算得到了一個odb，為g.odb，需要輸出結點集合的編號和應力結果？
下面進行展示：
APDL實現(xiàn)ANSYS的結果輸出
APDL實現(xiàn)ANSYS的結果輸出
長安CAE
1 概述
ANSYS作為通用有限元仿真計算軟件應用非常普遍，其強大的參數(shù)化功能更是使得ANSYS的使用方便很多。在計算過程中，通常需要對計算結果進行統(tǒng)計，ANSYS計算的模型大多時候節(jié)點單元數(shù)目很多，結果數(shù)據(jù)也很多，因此在GUI界面的查詢操作不太方便，工作量太大，而且結果不好記錄。而且很多時候需要借助第三方軟件比如MATLAB等數(shù)據(jù)處理軟件對結果進行處理，此情況下需要ANSYS能夠將計算結果輸出到文件中，再通過對結果文件進行處理得到想要的結果。
APDL不僅在參數(shù)化建模、加載和求解有很大的發(fā)揮，在后處理時同樣有發(fā)揮。本次借助APDL編寫后處理命令流，將需要的計算結果以設置好的格式輸出為txt文件。
2 應用
將編寫好的前處理和求解APDL命令流在ANSYS里面運行，命令流文件如圖1所示，該段命令流為兩根梁的求解計算，如圖2所示。
圖1 前處理與求解APDL命令流
圖2 模型
首先將所有節(jié)點的坐標位置輸出，命令流如圖3所示。
圖3 輸出節(jié)點位移APDL命令流
ANSYS還可以將結果保存為數(shù)組，如圖4所示，數(shù)組可以仿照圖3的循環(huán)部分輸出，數(shù)組中的數(shù)據(jù)可以繼續(xù)在ANSYS里面反復調用。
圖4 輸出數(shù)組
同樣的，其他的結果參數(shù)也可以輸出。如圖5所示，為節(jié)點應力結果的輸出。
圖5 APDL獲取其他的結果
基本思路是通過循環(huán)語句獲取所有節(jié)點的參量，保存成數(shù)組，也可以不保存，然后新建一個文件用于輸出節(jié)點的參量，主要的命令是*CREATE和*GET。
展開 MatlabGUI界面調用Ansys計算并輸出計算結果
.*'},'File Selector'); strh = [Pnameh,Fnameh];
pathname = Pnameh;
set(handles.text1,'String',strh);
[temp1,temp2] = xlsread(strh);
set(handles.uitable1,'Data',temp1);
% Update handles structure
guidata(hObject, handles);
為了讀取圖示方框中的數(shù)據(jù)，并用到ANSYS的APDL文件中，需要字符串的讀取和合并，首先需要使用str2num函數(shù)把字符串轉換成數(shù)值，如果沒有輸入值時，使用缺省值。
將兩個txt合并成test3.mac作為APDL語言開始的參數(shù)定義，生成test3.mac之后再使用system函數(shù)調用ANSYS的求解器，并讀取test3.mac進行計算
在計算之前，是不能生成圖片的，這時需要設置只有點擊“開始重構”按鈕之后，其他按鈕才可用。
點擊按鈕開始計算之后，會分別輸出兩個名為residualstress.jpg和deformation.jpg的圖片，對應的語句為
/image,save,'E:\GUIRStest\residualstress',jpg
設置當點擊“生成殘余應力云圖”和“生成角變形云圖”時，會讀取圖片的路徑并使用imshow生成圖片。
至此，一個簡易的MatlabGUI界面調用ANSYS計算并輸出圖片就完成了。
展開 ANSYS如何批量輸出結果圖
ANSYS提供了很多圖片格式，但有些格式只適用于特定的操作系統(tǒng)或者Device，且有些圖片格式是不能通過APDL語言導出來的(ANSYS導出圖片的命令流有三種，水哥只推薦/image，其他兩種使用起來太繁瑣，不做推薦)，只能采用GUI操作。本人就比較偏愛的三種格式簡單如下：
1）JPG
此種格式清晰度較高，且可以根據(jù)自己需要設置圖片質量高低，經常不涉及批量出圖，只需一兩張結果圖時，我便會使用這種格式。但這種格式不能通過/image命令導出來，且只會保留圖片到你剛開始指定的工作目錄下，文件名字不能更改。
GUi路徑如下：
2）BMP、PNG
這兩種格式均可通過GUI和命令流輸出，GUi輸出和上述JPG的輸出方法如出一轍，這里不再介紹。既然可以使用命令流輸出，那么就可以采用循環(huán)的方式批量出圖了，后面會做例子簡單演示。
但此兩種格式的圖片空間較大，一般在2M左右，但對于如今的磁盤空間來講，這點缺點微不足道了。
3）emf
這個格式不用說了，對寫論文的童鞋來講在熟悉不過了。此格式唯一不好的地方是不能采用命令流輸出，只能GUi操作。
回歸正題，如何批量出圖。
使用命令流：/image，該命令流使用格式如下：
label一般選擇save，F(xiàn)name 為文件的名字，當采用循環(huán)存儲的時候為了避免圖片相互之間因為同名被覆蓋，此處需要用到將數(shù)字轉為字符的命令%_%，可通過引號指定圖片保存的位置，ext為圖片格式，此處可以為bmp，png。
例：
某等厚度圓環(huán)孔板的模態(tài)分析，已經求得其前十階的振動模態(tài)，現(xiàn)需要批量導出其前十階模態(tài)。
展開 