ansys提取應(yīng)力的代碼的案例
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)附加在單元上，可以通過(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)命令流。
ansys workbench rst 文件應(yīng)力、位移和坐標(biāo)結(jié)果提取
采用python語(yǔ)言提取rst 文件結(jié)果提取
Ansys Workbench諧響應(yīng)掃頻結(jié)果后處理，提取Von Mises掃頻曲線和應(yīng)力幅值 ￥10
問(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)曲線中，有一項(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)力最大。而后再通過(guò)應(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è)掃頻過(guò)程中最大應(yīng)力值及其頻率。并將結(jié)果寫(xiě)出到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é)果后處理問(wèn)題示例：
Ansys workbench進(jìn)可以查看某個(gè)頻率下的 Von Mises應(yīng)力幅值
Ansys workbench進(jìn)掃頻應(yīng)力響應(yīng)曲線中，應(yīng)力選項(xiàng)卻沒(méi)有Von Mises應(yīng)力選型，只能按三個(gè)方向來(lái)分別查看。
展開(kāi) 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>
py源代碼|平均應(yīng)力應(yīng)變位移輸出至Excel
當(dāng)提取模型文件計(jì)算結(jié)果數(shù)據(jù)量較大時(shí)，依次設(shè)置歷程輸出任務(wù)繁重。實(shí)際上結(jié)合公眾號(hào)(有限元與力學(xué))提供的腳本文件(output_constants.py和output_main.py)，只需幾行代碼即可高效提取有用數(shù)據(jù)。
話不多說(shuō)實(shí)現(xiàn)方式如下：（不需要下載任何軟件）
（1）以記事本形式打開(kāi)output_constants.py文件，修改odbName和stepName
odbName='RS1' #odb文件名
stepName='Step-1' #step名
（2）在odb文件同目錄下新建一個(gè)txt文件，將以下內(nèi)容作為表頭固定輸入
#!
展開(kāi) ANSYS Beam188提取彎矩為例介紹ANSYS定義單元表提取數(shù)據(jù) (解決彎矩圖鋸齒狀） ￥20
在ANSYS中有些數(shù)據(jù)無(wú)法直接訪問(wèn)，需要通過(guò)定義單元表完成單元的結(jié)果的訪問(wèn)。下面就以Beam188單元提取彎矩為例介紹ANSYS定義單元表提取數(shù)據(jù)的詳細(xì)過(guò)程。
1. 首先需要知道在哪里定義單元表：Main Menu>General Postproc>Element Table>Define Table>add
2. 定義你想要的數(shù)據(jù)，這里以Beam188的彎矩為例
2.1 啟動(dòng)ANSYS幫助菜單, 在索引框輸入Beam188然后搜索, 在單元輸出介紹找到彎矩的名稱（代號(hào)）。
2.2 回到ANSYS界面，比如要輸出Mz, 則需要在添加SMISC,3 和SMISC,16 ,如圖
3. 輸出數(shù)據(jù)：Main Menu>General Postproc>Element Table> List E T, 選擇前面定義的SMISC,3 和SMISC,16 輸出單元I和J節(jié)點(diǎn)的Mz數(shù)值，如圖
4. 顯示彎矩云圖：Main Menu>General Postproc>Plot Results>Contour Plot>Line Elem Res, 這里要注意要在LabI 選SMISC,3 LabJ 選SMSCI,16。
輸出彎矩到這就結(jié)束了，小編突然發(fā)現(xiàn)，輸出的彎矩值在每個(gè)單元的I和J處是一樣的（Beam188為2節(jié)點(diǎn)單元），彎矩圖也就成了鋸齒形，于是去問(wèn)了度娘一波，各路盆友給出解決方法，然而并沒(méi)有起作用的，于是乎我又想起來(lái)了“幫助文檔大法”，于是認(rèn)認(rèn)真真將Beam188的幫助文檔閱讀了一遍，功夫不負(fù)有心人，最終。。。
展開(kāi) 拓?fù)鋬?yōu)化結(jié)構(gòu)MISES應(yīng)力分布3D顯示MATLAB代碼 ￥100
拓?fù)鋬?yōu)化結(jié)構(gòu)MISES應(yīng)力分布3D顯示MATLAB代碼
將應(yīng)力拖到屈服面上程序代碼
應(yīng)力曲線沿路徑提取 ￥8
通過(guò)鉆柱長(zhǎng)度和角度繪制鉆柱簡(jiǎn)化
序號(hào)
角度（°）
長(zhǎng)度（m）
1
98.22
9.49
2
98.54
9.47
3
99.58
9.47
4
100.3
9.40
5
100.33
3.00
根據(jù)鉆柱簡(jiǎn)化模型基本參數(shù)建立模型
對(duì)鉆柱添加邊界條件和載荷約束如圖2所示
載荷及邊界條件
鉆桿頭部
中間鉆桿
鉆桿尾部
約束
X軸自由度
Y軸自由度
Z軸自由度
Y軸自由度
X軸自由度
Y軸自由度
Z軸自由度
繞X軸轉(zhuǎn)動(dòng)自由度
載荷
鉆壓5t
扭矩1000Nm
重力9.8m/s2（整個(gè)模型）
無(wú)
如果有需要文檔的同學(xué)，可以給我留言，備注信息。
ABAQUS實(shí)用子程序SPRINC提取主應(yīng)力
調(diào)用的時(shí)候直接在UMAT子程序中輸入：CALL SPRINC(S,PS,LSTR,NDI,NSHR)
需要提供給該實(shí)用子程序的變量包括：
實(shí)用子程序返回UMAT的變量
實(shí)例：
CALL SPRINC(S,PS,LSTR,NDI,NSHR) 可以直接寫(xiě)成
CALL SPRINC(STRESS,PS,1,3,3)或CALL SPRINC(STRESS,PS,1, NDI,NSHR)
這行代碼的意思是將UMAT中的應(yīng)力分量數(shù)組STRESS提供給實(shí)用子程序SPRINC計(jì)算；1表示LSTR=1，表明提供給實(shí)用子程序SPRINC的為應(yīng)力分量數(shù)組；后面的兩個(gè)“3”分別表示正應(yīng)力分量個(gè)數(shù)為3和切應(yīng)力分量個(gè)數(shù)為3，也可以不用改變直接使用NDI和NSHR，因?yàn)樵谌S模型的UMAT中NDI和NSHR已經(jīng)被賦予了值為3。PS數(shù)組儲(chǔ)存返回的計(jì)算得到的主應(yīng)力，這個(gè)需要自己在umat程序中定義數(shù)組PS.
下面是關(guān)于返回主應(yīng)力順序的驗(yàn)證：
CALL SPRINC(STRESS,PS,1,3,3)
STATEV(4)=PS(1)
STATEV(5)=PS(2)
STATEV(6)=PS(3)
需要注意以下點(diǎn)：
1.PS數(shù)組一般需要在程序中自定義數(shù)組變量；
2.S數(shù)組一般也需要自定義數(shù)組變量，且需要為S數(shù)組賦予應(yīng)力分量或應(yīng)變分量的值；
3.如果SPRINC或SPRIND返回的主應(yīng)力為零，檢查輸入的變量或數(shù)組是否未定義或者未賦值，其次檢查輸出數(shù)組是否未定義。
展開(kāi) Abaqus平均應(yīng)力和應(yīng)變提取 ￥80
利用python腳本對(duì)ODB文件中單元集里所有積分點(diǎn)的應(yīng)力及應(yīng)變進(jìn)行自動(dòng)提取并計(jì)算平均值
能夠得到每一幀的應(yīng)力和應(yīng)變平均值，并保存到CSV文件中
所得到的應(yīng)力包括S11,S22,S33,S12,S13,S23以及Mises七個(gè)應(yīng)力平均值，以及E11,E22,E33,E12,E13,E23六個(gè)應(yīng)變平均值
切削殘余應(yīng)力輸出提取
請(qǐng)問(wèn)有人可以做切削過(guò)后的殘余應(yīng)力提取嗎？就是殘余應(yīng)力沿深度方向的分布。
ABAQUS提取單元平均應(yīng)力/應(yīng)變 ￥10
利用python讀取odb文件（可一次讀取多個(gè)odb）生成csv（excel）文件。提供源文件，注釋詳細(xì)，可根據(jù)需要進(jìn)行修改。
Abaqus插件——平均應(yīng)力應(yīng)變提取 ￥60
通過(guò)該插件可實(shí)現(xiàn):
1）提取所有幀的任意單元集合的平均應(yīng)力（事先定義單元集合，如圖中的SET-1）
2）提取所有幀的任意區(qū)域的x、y、z方向的平均應(yīng)變（事先定義節(jié)點(diǎn)集合，如圖中的SET-2）
3）將以上數(shù)據(jù)保存至excel文件（excel文件名為odb文件名稱+_Stress_Strain.csv）
*************************注意事項(xiàng)******************************
1、插件使用過(guò)程中，如有任何問(wèn)題請(qǐng)發(fā)郵件至shenz1hao@126.com
2、插件僅做學(xué)習(xí)交流使用，尊重原創(chuàng)者，切勿以營(yíng)利目的傳播
*****************************************************************
********************插件安裝及使用*******************************
1、電腦路徑下輸入 %homepath%\abaqus_plugins并回車
2、將Stress-Strain文件夾解壓至當(dāng)前目錄下
3、打開(kāi)abaqus，菜單欄中點(diǎn)擊plug-ins，里面找出Stress-Strain
4、輸入相應(yīng)參數(shù)（hx、hy、hz表示x、y、z方向模型長(zhǎng)度，當(dāng)以上三參數(shù)取1時(shí)輸出的為該方向位移）
*****************************************************************
展開(kāi) Matlab“稍作修改”Abaqus-odb結(jié)果 | 提取三維裝配體Mises應(yīng)力
數(shù)據(jù)及函數(shù)文件獲取方式：
Matlab“稍作修改”Abaqus-odb結(jié)果 | 提取三維裝配體Mises應(yīng)力
https://mp.weixin.qq.com/s/l45KSgNxkJ-oCCrDMrUgHA
ansys提取應(yīng)力的代碼的相關(guān)專題、標(biāo)簽、搜索
ansys提取應(yīng)力的代碼ansys提取面上的應(yīng)力ansys 提取點(diǎn)的應(yīng)力ansys提取代碼ansys代碼提取節(jié)點(diǎn)號(hào)ansys應(yīng)力代碼 Ansys ansyswb代碼提取應(yīng)力強(qiáng)度因子ansysworkbench提取應(yīng)力應(yīng)變結(jié)果代碼基于pythonansysworkbench提取應(yīng)力應(yīng)變結(jié)果代碼基于python語(yǔ)言ansysworkbench提取應(yīng)力應(yīng)變結(jié)果代碼abaqus 宏錄制完怎么加個(gè)提取應(yīng)力的代碼ansyswb代碼提取數(shù)據(jù)
