ansys提取桿單元應(yīng)力的案例
ANSYS分析VS理論解 | 簡單托架應(yīng)力和變形分析（桿單元實(shí)例）
(4) 查看各單元應(yīng)力：
①定義軸向應(yīng)力單元表：Main
Menu >General Postproc >Element Table>Define
Table，→Lab：輸入Stress_I →Item：選擇By sequence num →Comb：選擇LS，在LS后面輸入“1”→OK
→Apply →Lab：輸入Stress_J →Item：選擇By sequence num →Comb：選擇LS，在LS后面輸入“2”→OK
→Close。
③軸力列表顯示：Main Menu >General Postproc >Element Table>List Element Table→選擇FN→OK→記錄各個(gè)單元的軸力→File →Close。
④畫軸力圖：Main Menu>General Postproc>Plot Results>Contour Plot>Line Elem Res(見圖1.5)→LabI選擇Stress_I，LabJ選擇Stress_J→OK。
5．退出ANSYS軟件
Utility Menu >File >Exit →Quit-No Save →OK
來源：ANSYS學(xué)習(xí)與應(yīng)用公眾號(hào)，版權(quán)歸作者所有。
展開 ANSYS中桿單元和殼單元的單元耦合問題
在比較復(fù)雜的結(jié)構(gòu)的有限元分析中，不同的結(jié)構(gòu)部件通常使用不同類型的單元來模擬。
通常情況下，不同類型的單元的各個(gè)節(jié)點(diǎn)的自由度數(shù)目是不同的，不同類型單元的連接節(jié)點(diǎn)處的自由度的耦合問題，是一個(gè)比較令人頭疼的問題。
在ANSYS中通常可以用耦合命令CP來耦合不同類型單元在連接節(jié)點(diǎn)處的自由度(DOF)。
也可以用CE命令來認(rèn)為添加自由度之間的約束方程來達(dá)到耦合的目的。
下面是一個(gè)簡單的算例，使用了CE命令來耦合連接節(jié)點(diǎn)處的自由度。
模型是航天器的機(jī)翼的一個(gè)Section的某一個(gè)隔框。上下表皮是薄殼結(jié)構(gòu)，用Shell63單元來模擬，在上下表皮之間有起支撐作用的桿件，用link8單元來模擬。
建模的時(shí)候，link8單元和shell63單元在連接有各自獨(dú)立的節(jié)點(diǎn)。即：link8單元和shell63單元的節(jié)點(diǎn)在連接處是重合的，但是，節(jié)點(diǎn)編號(hào)是各自獨(dú)立的。
link8單元在每個(gè)節(jié)點(diǎn)有 ux,uy,uz3個(gè)平動(dòng)自由度；
shell63在每個(gè)節(jié)點(diǎn)有ux,uy,uz這3個(gè)平動(dòng)自由度和rotx,roty,rotz這3個(gè)轉(zhuǎn)個(gè)自由，共6個(gè)自由度。
在耦合節(jié)點(diǎn)處，兩個(gè)耦合節(jié)點(diǎn)的ux,uy,uz自由度應(yīng)該是相等的。
這個(gè)等式可以用CE命令來描述。
完整的命令流如下：
finish
/clear,start
/prep7
!定義第一種材料屬性；
mp,ex,1,30e6
mp,prxy,1,0.3
!定義shell63單元和實(shí)常數(shù)；
et,1,shell63
r,1,1e-3
!建立幾何模型；
rectng,31.8,33.2,0,0.3556
agen,2,1,1,1,0,0,1
a,1,4,8,5
a,6,7,3,2
KL,7,0.5, ,
KL,3,0.5, ,
在關(guān)鍵點(diǎn)處生成節(jié)點(diǎn)；
nkpt,100,4 !與編號(hào)為117的節(jié)點(diǎn)耦合
nkpt,101,9 !
展開 ANSYS各類型單元連接專題講解（二）之桿與梁殼體單元的連接
前一篇文章主要介紹了單元之間連接的主要原則，今天開始主要從具體方面講解連接方法。
按照桿、梁、殼、實(shí)體的順序，先說說桿單元與各單元的連接方法。
那么什么時(shí)候需要用到桿單元與各種單元的連接呢？水哥稍微列舉下實(shí)際工程中需要考慮此類連接的例子。
案例一：工業(yè)廠房
此類結(jié)構(gòu)一般橫向跨度較大，屋頂采用鋼結(jié)構(gòu)形式，在具體模擬屋架時(shí)，此時(shí)各個(gè)桿件可看成鉸接，采用桿單元模擬。而下方框架柱則采用梁單元進(jìn)行模擬，在相交部位則需要用到桿單元與梁單元的連接。
案例二：門廳鋼結(jié)構(gòu)雨棚
在具體模擬該結(jié)構(gòu)時(shí)，雨棚上方拉桿采用桿單元模擬，而下方的鋼梁采用梁單元模擬，混凝土框架柱可采用實(shí)體單元模擬。
一直以來，桿單元一般用于模擬桁架結(jié)構(gòu)的時(shí)候比較多，其特點(diǎn)是桿件兩端不考慮承受彎矩作用，節(jié)點(diǎn)只有平動(dòng)自由度，是所有單元中最為簡單的一種。
桿單元分為2D桿單元和3D桿單元，2D桿單元節(jié)點(diǎn)只有Ux和Uy兩個(gè)平動(dòng)自由度，而3D桿單元除了這兩個(gè)，還有Uz。其他單元，梁單元、殼單元、體單元都包含了這三個(gè)自由度，且具有相同的物理意義，按照前面一篇文章所介紹的連接總則，桿單元與其他單元連接時(shí)只需要共用節(jié)點(diǎn)即可，無需建立約束方程。
下面是一個(gè)簡單的類似雨棚案例，注意本案例各構(gòu)件尺寸僅為演示操作需要所擬，未經(jīng)仔細(xì)推敲，各工程大佬可忽略。
某屋外雨棚平面簡化模型如上，長度為4m，折算荷載為10 KN/m，雨棚梁采用工字型鋼I40，系桿截面面積為238.64mm^2，材料均為Q235，采用ANSYS模擬該結(jié)構(gòu)。
下面為建模過程
!
展開 ABAQUS提取單元平均應(yīng)力/應(yīng)變 ￥10
利用python讀取odb文件（可一次讀取多個(gè)odb）生成csv（excel）文件。提供源文件，注釋詳細(xì)，可根據(jù)需要進(jìn)行修改。
ANSYS Beam188提取彎矩為例介紹ANSYS定義單元表提取數(shù)據(jù) (解決彎矩圖鋸齒狀） ￥20
在ANSYS中有些數(shù)據(jù)無法直接訪問，需要通過定義單元表完成單元的結(jié)果的訪問。下面就以Beam188單元提取彎矩為例介紹ANSYS定義單元表提取數(shù)據(jù)的詳細(xì)過程。
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)單元），彎矩圖也就成了鋸齒形，于是去問了度娘一波，各路盆友給出解決方法，然而并沒有起作用的，于是乎我又想起來了“幫助文檔大法”，于是認(rèn)認(rèn)真真將Beam188的幫助文檔閱讀了一遍，功夫不負(fù)有心人，最終。。。
展開 基于ANSYS APDL的兩端固定桿的單元生死仿真【轉(zhuǎn)載】
希望有所收獲
【問題描述】
一根兩端固定的桿如下圖所示。
材料數(shù)據(jù)如下
為了闡述如何使用ANSYS的單元生死技術(shù)，決定把該桿等分為3個(gè)單元，然后通過控制中間單元的生死，進(jìn)行如下的熱應(yīng)力仿真
（1）設(shè)置所有單元的材料參考溫度是0度，給所有節(jié)點(diǎn)施加100度，并保持所有單元都存活，做1次仿真
（2）設(shè)置所有單元的材料參考溫度是0度，給所有節(jié)點(diǎn)施加100度，殺死中間單元，做1次仿真
（3）設(shè)置其它單元的材料參考溫度是0度，給所有節(jié)點(diǎn)施加100度，激活中間單元，并設(shè)置該單元的材料參考溫度是100度，做1次仿真
（4）設(shè)置其它單元的材料參考溫度是0度，給所有節(jié)點(diǎn)施加0度，保持中間單元存活，并設(shè)置該單元的材料參考溫度是100度，做1次仿真
通過上述四次仿真，以說明
（1）如何使用單元的生死技術(shù)
（2）當(dāng)單元激活時(shí)，會(huì)根據(jù)節(jié)點(diǎn)溫度和該單元的材料參考溫度之差來確定它的初始熱應(yīng)變。
【問題分析】
1.該例子來自于ANSYS15 APDL的認(rèn)證算例《VM194 Element Birth/Death in a Fixed Bar》為了更清晰的闡明思路，本文對(duì)其進(jìn)行了較大幅度的調(diào)整。
2.單元生死技術(shù)的使用，關(guān)鍵是首先要?jiǎng)?chuàng)建出所有的單元，然后在需要?dú)⑺栏?em>單元時(shí)使用EKILL命令，而在需要激活時(shí)使用ELIVE命令。
3.使用LINK180來建模桿。
4.創(chuàng)建2種材料。這兩種材料的彈性模量和泊松比一樣，但是參考溫度不一樣。一個(gè)參考溫度是0度，一個(gè)是100度。
5.先創(chuàng)建4個(gè)節(jié)點(diǎn)，然后創(chuàng)建3個(gè)單元。
6.固定兩個(gè)端節(jié)點(diǎn)，并給所有節(jié)點(diǎn)固定Z方向自由度，借此模擬二維桿件。7.按照題目要求進(jìn)行先后四次的計(jì)算和后處理，以考察生死單元的使用。
8.本文采用APDL命令進(jìn)行講解。
【求解過程】
1.
展開 基于ANSYS Workbench2024R2 桿單元不同載荷下的瞬態(tài)分析 ￥50
<p>基于ANSYS Workbench2024R2 桿單元不同載荷下的瞬態(tài)分析</p><p>預(yù)應(yīng)力分析</p><p><br></p><figure style="text-align: center;" class="ql-align-center">
<figure class="figure-image" contenteditable="false" data-img="https://img.jishulink.com/202503/attachment/4968e839a1834fbf9d63a7f4a426758e.png" style="display: inline-block;" data-regular="true">
<img src="https://img.jishulink.com/202503/attachment/4968e839a1834fbf9d63a7f4a426758e.png" style="" width="622" data-mobile-src="https://img.jishulink.com/202503/attachment/4968e839a1834fbf9d63a7f4a426758e.png?image_process=/format,webp" data-pc-src="https://img.jishulink.com/202503/attachment/4968e839a1834fbf9d63a7f4a426758e.png?
展開 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靜力分析（桿單元） ￥1.25
作者介紹： 力學(xué)碩士，有七年的結(jié)構(gòu)有限元分析經(jīng)驗(yàn)
～～～～～～～～～～～～～～～～～～～～～～～～
在ANSYS中，桁架結(jié)構(gòu)（只承受拉壓，不承受彎矩）要使用桿單元（link單元）進(jìn)行分析。在新版的ANSYS中，一般都推薦使用link180單元，該單元有兩個(gè)節(jié)點(diǎn)，每個(gè)節(jié)點(diǎn)有三個(gè)平移自由度。對(duì)于本文的平面三角桁架分析，有如下注意事項(xiàng)：
1 link180是三維桿，分析平面問題，需要約束一個(gè)自由度，一般為Z向。
2 桁架結(jié)構(gòu)的建模，可以直接從節(jié)點(diǎn)單元開始，因?yàn)殍旒艿拿扛?em>桿都只劃分為一個(gè)單元。
3 link180單元的截面雖然可以用sectype和secdata來定義，但計(jì)算本質(zhì)還是轉(zhuǎn)化為實(shí)常數(shù)。
4 對(duì)于桿結(jié)構(gòu)，荷載都施加在節(jié)點(diǎn)上，桿單元不能施加線荷載。
對(duì)于線模型（桿結(jié)構(gòu)，梁結(jié)構(gòu)，管結(jié)構(gòu)），SECTYPE和SECDATA是很重要的命令：
當(dāng)命令sectype的type是link的時(shí)候，secdata定義桿截面面積。
如果讀者想詳細(xì)了解SECTYPE和SECDATA，可以輸入help, sectype或者h(yuǎn)elp, secdata。如下圖：
然后按一下鍵盤的enter，軟件會(huì)跳出help文件，詳細(xì)解釋sectype。
后文目錄：
一：建模
二：求解
三：后處理
四：源文件
展開 ansys workbench rst 文件應(yīng)力、位移和坐標(biāo)結(jié)果提取
采用python語言提取rst 文件結(jié)果提取
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>
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)。簡要流程為：首先，讀取每一個(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中整體、單元剛度和質(zhì)量矩陣的提取
指定輸出單元矩陣
/SOLU
SOLVE
finish
/OUTPUT, TERM ! 將輸出信息送到output windows中
! 這時(shí)用編輯器打開cp.out文件，可以看到按單元寫出的質(zhì)量、剛度等矩
陣
ANSYS中整體、單元剛度和質(zhì)量矩陣的提取.rar
ANSYS提取具體三維單元的體積，面元的面積和線的長度
在ANSYS中，能提取具體三維單元的體積，面元的面積和線的長度。
如：*GET,E_VOLUME,ELEM,10,VOLU 為提取編號(hào)為10的單元的體積
*GET,a_area,AREA,50,AREA 為提取編號(hào)為50的面元的面積
*GET,l_length,LINE,100,LENG 為提取編號(hào)為100的線的長度
以上對(duì)應(yīng)的GUI操作： Utility Menu>Parameters>Get Scalar Data 如果要一次性提取多個(gè)元素的相關(guān)參數(shù)，可以用命令 *VGET, ParR, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM, KLOOP 對(duì)應(yīng)GUI操作：Utility Menu>Parameters>Get Array Data
輸入命令 alist,p 出選項(xiàng)框，選你要看的那個(gè)面積，提取選中的單元面積。
*cfopen,'area','txt',
*GET,MaxEleNum,ELEM,,NUM,MAX
*GET,MinEleNum,ELEM,,NUM,MIN
*do,i,MinEleNum,MaxEleNum,1
*if,esel(i),eq,1,then
*get,volu,elem,i,volu
*vwrite,i,volu
(f5.0,f15.12)
*end if
*enddo
*cfcols
展開 『分享』ANSYS中整體、單元剛度和質(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提供的庫文件，將其
引入項(xiàng)目中（也可直接扔進(jìn)debug目錄）
.full文件由子空間迭代模態(tài)分析獲得。
運(yùn)行編譯后的可執(zhí)行文件.exe
2、單元剛度和質(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,, !
展開 