ANSYS提取整體應(yīng)力場(chǎng)
關(guān)注創(chuàng)建者：王靖雯 創(chuàng)建時(shí)間：2023-03-08
ANSYS提取整體應(yīng)力場(chǎng)的視頻教程
增材仿真+生死單元+ansys apdl+熱力耦合+溫度場(chǎng)+應(yīng)力場(chǎng)
介紹：運(yùn)用ANSYS二次開(kāi)發(fā) APDL語(yǔ)言編輯出參數(shù)化程序來(lái)建立模型、控制和劃分網(wǎng)格、 定義材料參數(shù)、施加載荷與邊界條件、分析控制以及求解等完成有限元分析全部過(guò)程。在模擬成型過(guò)程中，通過(guò)改變溫度載荷的位置來(lái)模擬噴嘴的掃描移動(dòng)，利用生死單元循環(huán)算法技術(shù)控制單元“生死”的激活來(lái)模擬材料的堆積增加，通過(guò)控制單元激活的時(shí)間間隔控制成型速度
免費(fèi) 20.0833秒 100播放
查看
ansys fluent電路板強(qiáng)制對(duì)流換熱、熱應(yīng)力、模態(tài)、ncode隨機(jī)振動(dòng)及正弦振動(dòng)疲勞-多場(chǎng)耦合
熱應(yīng)力計(jì)算、熱應(yīng)力對(duì)模態(tài)的影響與不考慮熱應(yīng)力進(jìn)行對(duì)比分析； ncode進(jìn)行隨機(jī)振動(dòng)疲勞以及正弦振動(dòng)疲勞分析注意事項(xiàng)，S-N曲線的估計(jì)方法，以及后處理等操作
￥39.9 2小時(shí)24分鐘 237播放
查看
ANSYS提取整體應(yīng)力場(chǎng)的實(shí)例教程
這時(shí)用編輯器打開(kāi)cp.out文件，可以看到按單元寫(xiě)出的質(zhì)量、剛度等矩
陣
ANSYS中整體、單元?jiǎng)偠群唾|(zhì)量矩陣的提取.rar
、整體剛度和質(zhì)量矩陣的提取。
該功能需要進(jìn)行二次開(kāi)發(fā)，由ansys形成
的二進(jìn)制文件.full提取整體剛度和質(zhì)量
矩陣。
基于ansys的一個(gè)用戶(hù)開(kāi)發(fā)程序例子編
了一個(gè)程序（附件中）。
開(kāi)發(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ì)說(shuō)明可由以下轉(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,, !
展開(kāi) ###########################################################################
寫(xiě)在前面，最近回看了一些自己之前編寫(xiě)的和這個(gè)py代碼，發(fā)現(xiàn)自己寫(xiě)的可能和公式表達(dá)的有偏差，在我的測(cè)試過(guò)程中，提取的是單元積分點(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)在覺(jué)得這并不是獲取單元平均應(yīng)力的方式，也就是代碼并沒(méi)有實(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)力。
###########################################################################
最近寫(xiě)了一個(gè)簡(jiǎn)單的python讀取abaqus結(jié)果中的場(chǎng)輸出數(shù)據(jù)，想通過(guò)均勻化計(jì)算方法來(lái)計(jì)算所定義集合的平均應(yīng)力應(yīng)變曲線，之前是手動(dòng)提取了各個(gè)數(shù)據(jù)導(dǎo)出，然后用excel、matlab處理，但是很慢，而且很費(fèi)勁，于是就想著用Python來(lái)處理結(jié)果。
有需要的同學(xué)可以下載附件文件，打開(kāi)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?
展開(kāi) 在對(duì)結(jié)構(gòu)進(jìn)行時(shí)程分析后，我們經(jīng)常提取的是全時(shí)程最大von Mises stress。
那么如何提取某一個(gè)節(jié)點(diǎn)的von Mises stress呢？
首先明確ANSYS的節(jié)點(diǎn)附加在單元上，可以通過(guò)選擇單元上節(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)命令流。
</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>
ANSYS提取整體應(yīng)力場(chǎng)的相關(guān)專(zhuān)題、標(biāo)簽、搜索
ANSYS提取整體應(yīng)力場(chǎng)的最新內(nèi)容
問(wèn)題：
Ansys workbench進(jìn)行諧響應(yīng)仿真計(jì)算的后處理結(jié)果中，提供了單一頻率下的Von Mises應(yīng)力查看功能和應(yīng)力頻響曲線功能，但是應(yīng)力頻響曲線的應(yīng)力列表中沒(méi)有Von Mises應(yīng)力查看項(xiàng)。因?yàn)閂on Mises應(yīng)力太常用，所以這就給我們?cè)谡麄€(gè)掃頻范圍內(nèi)，定位Von Mises應(yīng)力的最大頻率和應(yīng)力值帶來(lái)一定的困難。如下所示。
需求：
希望后處理結(jié)果中可以在應(yīng)力響應(yīng)曲線中
采用python語(yǔ)言提取rst 文件結(jié)果提取
在對(duì)結(jié)構(gòu)進(jìn)行時(shí)程分析后，我們經(jīng)常提取的是全時(shí)程最大von Mises stress。
那么如何提取某一個(gè)節(jié)點(diǎn)的von Mises stress呢？
首先明確ANSYS的節(jié)點(diǎn)附加在單元上，可以通過(guò)選擇單元上節(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, 定義變量。
<p>在<a href="https://www.yqgqt.org.cn/qa/3655" rel="noopener noreferrer" target="_blank">ANSYS結(jié)構(gòu)</a>動(dòng)力分析時(shí)，時(shí)程分析（瞬態(tài)分析）的后處理經(jīng)常想要提取全時(shí)程結(jié)構(gòu)響應(yīng)的最大值及對(duì)應(yīng)的時(shí)間步。在<a href="https://www.yqgqt.org.cn/major/Ansys" rel="noopener
###########################################################################
寫(xiě)在前面，最近回看了一些自己之前編寫(xiě)的和這個(gè)py代碼，發(fā)現(xiàn)自己寫(xiě)的可能和公式表達(dá)的有偏差，在我的測(cè)試過(guò)程中，提取的是單元積分點(diǎn)處的應(yīng)力應(yīng)變值，由于只有一個(gè)積分點(diǎn)（這個(gè)積分點(diǎn)在質(zhì)心，與centroid提取得到的結(jié)果一樣），因此自然的把這個(gè)值當(dāng)成了單元的平均應(yīng)力或平均應(yīng)變
PCB熱-應(yīng)力可靠性和多場(chǎng)耦合分析培訓(xùn)班
培訓(xùn)背景
電路的集成規(guī)模越來(lái)越大，I/O數(shù)越來(lái)越多，PCB互連密度不斷加大，隨之帶來(lái)許多PCB及集成電路封裝可靠性問(wèn)題。ANSYS專(zhuān)門(mén)針對(duì)PCB設(shè)計(jì)分析解決方案，可以快速?gòu)腅CAD中直接導(dǎo)入PCB熱物參數(shù)，從而能在Mechanical中進(jìn)行準(zhǔn)確的PCB板熱力、疲勞、隨機(jī)振動(dòng)、跌落等可靠性問(wèn)題的仿真。ANSYS針對(duì)集成電路封裝也提供強(qiáng)大解決方案
、整體剛度和質(zhì)量矩陣的提取。
該功能需要進(jìn)行二次開(kāi)發(fā)，由ansys形成
的二進(jìn)制文件.full提取整體剛度和質(zhì)量
矩陣。
基于ansys的一個(gè)用戶(hù)開(kāi)發(fā)程序例子編
了一個(gè)程序（附件中）。
開(kāi)發(fā)環(huán)境：compaq fortran 6.5
運(yùn)行環(huán)境：win2000。
一個(gè)主文件：self.for，
另一個(gè)文件matrixout.f90用于矩陣輸出
binlib.lib
1、整體剛度和質(zhì)量矩陣的提取。
該功能需要進(jìn)行二次開(kāi)發(fā)，由ansys形成
的二進(jìn)制文件.full提取整體剛度和質(zhì)量
矩陣。
基于ansys的一個(gè)用戶(hù)開(kāi)發(fā)程序例子編
了一個(gè)程序（附件中）。
開(kāi)發(fā)環(huán)境：compaq fortran 6.5
運(yùn)行環(huán)境：win2000。
一個(gè)主文件：self.for，
另一個(gè)文件matrixout.f90用于矩陣輸出
binlib.lib為ansys提供的庫(kù)文件
