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