ANSYS提取整體應(yīng)力場(chǎng)的案例
ANSYS中整體、單元?jiǎng)偠群唾|(zhì)量矩陣的提取
這時(shí)用編輯器打開cp.out文件，可以看到按單元寫出的質(zhì)量、剛度等矩
陣
ANSYS中整體、單元?jiǎng)偠群唾|(zhì)量矩陣的提取.rar
『分享』ANSYS中整體、單元?jiǎng)偠群唾|(zhì)量矩陣的提取
、整體剛度和質(zhì)量矩陣的提取。
該功能需要進(jìn)行二次開發(fā)，由ansys形成
的二進(jìn)制文件.full提取整體剛度和質(zhì)量
矩陣。
基于ansys的一個(gè)用戶開發(fā)程序例子編
了一個(gè)程序（附件中）。
開發(fā)環(huán)境：compaq fortran 6.5
運(yùn)行環(huán)境：win2000。
一個(gè)主文件：self.for，
另一個(gè)文件matrixout.f90用于矩陣輸出
binlib.lib為ansys提供的庫(kù)文件，將其
引入項(xiàng)目中（也可直接扔進(jìn)debug目錄）
.full文件由子空間迭代模態(tài)分析獲得。
運(yùn)行編譯后的可執(zhí)行文件.exe
2、單元?jiǎng)偠群唾|(zhì)量矩陣的提取。
/DEBUG命令。詳細(xì)說明可由以下轉(zhuǎn)載文章
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,, !
展開 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è)簡(jiǎn)單的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?
展開 ANSYS如何提取某一節(jié)點(diǎn)的應(yīng)力時(shí)程 ￥100
在對(duì)結(jié)構(gòu)進(jìn)行時(shí)程分析后，我們經(jīng)常提取的是全時(shí)程最大von Mises stress。
那么如何提取某一個(gè)節(jié)點(diǎn)的von Mises stress呢？
首先明確ANSYS的節(jié)點(diǎn)附加在單元上，可以通過選擇單元上節(jié)點(diǎn)的方法提取節(jié)點(diǎn)應(yīng)力。
1 確定節(jié)點(diǎn)所在單元，顯示節(jié)點(diǎn)編號(hào)。
例單元號(hào)8560，節(jié)點(diǎn)號(hào)8678。
2 進(jìn)入TimeHist Postpro, 定義變量。
3變量顯示。
付費(fèi)內(nèi)容為相關(guān)命令流。
ansys workbench rst 文件應(yīng)力、位移和坐標(biāo)結(jié)果提取
采用python語言提取rst 文件結(jié)果提取
Ansys Workbench諧響應(yīng)掃頻結(jié)果后處理，提取Von Mises掃頻曲線和應(yīng)力幅值 ￥10
問題：
Ansys workbench進(jìn)行諧響應(yīng)仿真計(jì)算的后處理結(jié)果中，提供了單一頻率下的Von Mises應(yīng)力查看功能和應(yīng)力頻響曲線功能，但是應(yīng)力頻響曲線的應(yīng)力列表中沒有Von Mises應(yīng)力查看項(xiàng)。因?yàn)閂on Mises應(yīng)力太常用，所以這就給我們?cè)谡麄€(gè)掃頻范圍內(nèi)，定位Von Mises應(yīng)力的最大頻率和應(yīng)力值帶來一定的困難。如下所示。
需求：
希望后處理結(jié)果中可以在應(yīng)力響應(yīng)曲線中，有一項(xiàng)Von Mises應(yīng)力選項(xiàng)。實(shí)現(xiàn)每個(gè)掃頻點(diǎn)的最大Von Mises應(yīng)力和掃頻頻率的曲線圖顯示，從而一眼就看出產(chǎn)品在整個(gè)掃頻范圍內(nèi)，哪個(gè)頻率下結(jié)構(gòu)的等效應(yīng)力最大。而后再通過應(yīng)力云圖查看這個(gè)頻率下的Von Mises應(yīng)力。
解決方法：
利用APDL命令實(shí)現(xiàn)。簡(jiǎn)要流程為：首先，讀取每一個(gè)掃頻點(diǎn)的最大Von Mises應(yīng)力值。記下應(yīng)力值、頻率值和最大節(jié)點(diǎn)號(hào)。再統(tǒng)計(jì)記錄的所有掃頻點(diǎn)的Von Mises應(yīng)力值，提取整個(gè)掃頻過程中最大應(yīng)力值及其頻率。并將結(jié)果寫出到txt文件。進(jìn)一步提取這個(gè)最大Von Mises應(yīng)力點(diǎn)對(duì)應(yīng)的整個(gè)掃頻范圍內(nèi)的Von Mises應(yīng)力曲線。
這個(gè)樣就可以在txt文檔中直接看到所有掃頻點(diǎn)下，結(jié)構(gòu)的等效應(yīng)力幅值；以及全頻段中最大Von Mises應(yīng)力所在節(jié)點(diǎn)的等效應(yīng)力掃頻曲線圖。
效果展示如下：
在結(jié)果文件夾中，會(huì)生成一個(gè)txt結(jié)果文件和一張Von Mises應(yīng)力曲線圖。如此我們可以直觀注意到，在當(dāng)前掃頻范圍內(nèi)，結(jié)構(gòu)在78.95Hz時(shí)應(yīng)力最大約為17.552Mpa。
結(jié)果后處理問題示例：
Ansys workbench進(jìn)可以查看某個(gè)頻率下的 Von Mises應(yīng)力幅值
Ansys workbench進(jìn)掃頻應(yīng)力響應(yīng)曲線中，應(yīng)力選項(xiàng)卻沒有Von Mises應(yīng)力選型，只能按三個(gè)方向來分別查看。
展開 ANSYS瞬態(tài)分析全時(shí)程結(jié)構(gòu)響應(yīng)最大值的提取方法（變形、應(yīng)力、應(yīng)變、能量） ￥100
</p><p>同樣的方法，可以提取全時(shí)程最大的位移、應(yīng)力、應(yīng)變、能量等結(jié)果。</p><p><img src="https://img.jishulink.com/upload/202302/623025b5c0d646b9973cd2adc6c6037f.png" alt="1.png"></p><p>收費(fèi)內(nèi)容為相關(guān)命令流。</p>
【11月15-16日 上海】ANSYS官方培訓(xùn)—PCB熱-應(yīng)力可靠性和多場(chǎng)耦合分析培訓(xùn)班
PCB熱-應(yīng)力可靠性和多場(chǎng)耦合分析培訓(xùn)班
培訓(xùn)背景
電路的集成規(guī)模越來越大，I/O數(shù)越來越多，PCB互連密度不斷加大，隨之帶來許多PCB及集成電路封裝可靠性問題。ANSYS專門針對(duì)PCB設(shè)計(jì)分析解決方案，可以快速?gòu)腅CAD中直接導(dǎo)入PCB熱物參數(shù)，從而能在Mechanical中進(jìn)行準(zhǔn)確的PCB板熱力、疲勞、隨機(jī)振動(dòng)、跌落等可靠性問題的仿真。ANSYS針對(duì)集成電路封裝也提供強(qiáng)大解決方案，可以快速準(zhǔn)確進(jìn)行集成電路熱應(yīng)力問題、封裝翹曲、焊球疲勞問題、裂紋預(yù)測(cè)及擴(kuò)展等可靠性分析。
本次培訓(xùn)從解決PCB及集成電路封裝結(jié)構(gòu)可靠性基礎(chǔ)功能入手，逐步深入到ANSYS解決PCB及集成電路封裝結(jié)構(gòu)可靠性高級(jí)解決方案，并將演示國(guó)外專家解決集成電路封裝可靠性問題的多層次模型方案。
為了解決集成電路封裝結(jié)構(gòu)可靠性仿真需求，提升相關(guān)科技工作者的技術(shù)水平，普及ANSYS軟件高級(jí)功能。因此，ANSYS公司特開辦“PCB熱-應(yīng)力可靠性和多場(chǎng)耦合分析培訓(xùn)班”。
培訓(xùn)合格者發(fā)放ANSYS技術(shù)培訓(xùn)認(rèn)證證書。
展開 