ansys單元開裂的案例
COH2D4T單元傳熱+開裂模擬（但是單元不自動刪除?。?/span>
之前一個帖子說了如何實現粘聚力單元的傳熱問題（詳細的大家可以翻閱一下我的那個帖子），本來以為aba終于帶來了個福利，coh可以傳熱了，按常規思維，coh傳熱是后加入的功能，那么它模擬開裂的功能應該一直都在啊，但是...，就是那么氣人，它可以傳熱了，但是不會開裂了，那我要你有何用呢，下面附上兩個簡單例子實錘證明，翻遍幫助文檔 ，也沒發現有一個例子coh既傳熱又能開裂的。歡迎大家一起嘗試，探討，也更希望大家能找到coh既能傳熱又能開裂的解決辦法。
1 簡單的拉伸模型
使用coh2d4單元時可以開裂：
使用coh2d4t單元時bu開裂：
2 復雜的熱力模型
使用coh2d4單元時可以開裂：
使用coh2d4t單元時不開裂：
關于此問題的其他討論帖子
https://www.yqgqt.org.cn/content/post/1788990
ABAQUS斷裂模擬收徒 ，快速學會各種ABAQUS斷裂模擬方法 **/人（將有機會享有各種插件以及程序，價值**、專門定制視頻、全程親自教學、各種模型調試及解答問題等等，傾囊相教）
展開 批量插入0厚度cohesive單元模擬開裂
</p><p>批量插入cohesive單元的原理（下圖是在三角形單元和四面體單元間嵌入0厚度的cohesive單元，這里為了方便觀察，把cohesive畫為有限厚度，實際插入的是0厚度的）：</p><div contenteditable="false" width="100%"><div><img title="1.jpg" style="max-width:760px;" alt="1.jpg" src="https://img.jishulink.com/upload/201808/f49d76c4d26f4beeb3b5f73fc3ba5789.jpg" data-mobile-src="https://img.jishulink.com/upload/201808/f49d76c4d26f4beeb3b5f73fc3ba5789.jpg?image_process=/format,webp/resize,w_400" data-pc-src="https://img.jishulink.com/upload/201808/f49d76c4d26f4beeb3b5f73fc3ba5789.jpg?
展開 基于huang.for結合cohesive單元模擬晶粒之間的晶界開裂
<p> 已有大佬做了基于cohesive+泰森多邊形(Voronoi)插件模擬晶粒之間的晶界開裂的視頻，這個視頻中晶粒模型主要是通過homtools插件建立的，cohesive單元的建立是通過Cohesive_generator_2D3D插件實現的。附上視頻鏈接：</p><div contenteditable="false" width="100%">
<figure class="figure-link" data-title="基于cohesive+泰森多邊形(Voronoi)插件模擬晶粒之間的晶界開裂--Abaqus平臺" data-link="https://www.bilibili.com/video/BV1LV411y7CL/?share_source=copy_web&vd_source=ee5e911cda47c9e62824b381dae143c2" data-regular="true">
<a href="https://www.bilibili.com/video/BV1LV411y7CL/?share_source=copy_web&vd_source=ee5e911cda47c9e62824b381dae143c2" target="_blank" class="figure-link-a" rel="nofollow">基于cohesive+泰森多邊形(Voronoi)插件模擬晶粒之間的晶界開裂--Abaqus平臺</a>
</figure>
</div><p><br></p><p> 在Neper中建立晶粒模型劃分網格時也可以批量插入0厚度cohesive單元，我之前進行了嘗試，這里給出一個示例。
展開 一維單元模擬混凝土構件開裂的解決方案
一、算例背景及分析技術
鋼筋混凝土結構開裂是廣泛存在的現象，準確的說從鋼筋混凝土結構服役開始就進入帶裂工作狀態，只是這種裂縫肉眼難以辨別并且對結構安全沒有影響。但裂縫作為鋼筋混凝土構件工作狀態的重要表征指標，是結構損傷的表現，也許是破壞的先兆，也許是耐久
性不足的預警，更是災后調查和受力機理揭示的重要線索，見圖1。
圖1 鋼筋混凝土結構裂縫
本案例使用ABAQUS對一根鋼筋混凝土受拉構件進行裂縫估算分析，根據鋼筋應力狀態計算等效裂縫寬度，并采用Python腳本在ODB結果文件中創建裂縫場變量，實現在ABAQUS中使用一維單元進行快速分析，并在后處理模塊 顯示等效裂縫寬度的目標。
案例涉及的相關技術：
①ABAQUS梁單元Rebar積分點插入；
②利用Python腳本提取ABAQUS場變量數據；
③利用Python腳本創建ABAQUS場變量數據。
計算報告編寫采用操作引導式，希望能為讀者使用ABAUQS場變量創建提供有益參考。操作分析要點為：
①ABAQUSABAQUS梁單元Rebar積分點插入；
②ABAQUS場變量輸出Python腳本getSubset()函數應用；
③ABAQUS場變量編輯Python腳本addData()()函數應用。
二、計算任務
1.模型裝配及接觸連接
計算模型取自《混凝土結構：混凝土結構設計原理》（第六版）習題8-3。
計算模型為鋼筋混凝土屋架下弦按軸心受拉構件，見圖2。模型只包含1個part。截面寬200mm，截面高160mm。因為案例模型較為簡單，混凝土梁采用B21單元模擬，鋼筋通過在與混凝土單元共節點建立鋼筋箱型截面單元實現。
展開 
ANSYS鋼筋混凝土結構開裂計算介紹 附ANSYS土木工程應用實例下載
懸臂梁由于其特殊性，是屬于開裂計算中比較難以處理的一種結構，這里得到了比較合理的結果，這說明對于其它類型的結構，ANSYS技術同樣是可以處理的。
下載地址：ANSYS土木工程應用實例
ansys之——混凝土開裂
finish
/clear
/title, fixed - fixed concrete beam example
/prep7
et,1,65
mp,ex,1,3e7 ! steel rebar (units are pounds, inches)
mp,ex,2,1e6, ! concrete
mp,dens,2,.00025
tb,concr,2
tbdata,1,.3,.5,200,4000 ! shear coeffs, tensile and compress strength
r,1,1,.03,0,0 ! mat 1 (steel), 3 percent reinforcement in x dir
r,2,1,.01,0,0 ! mat 1 (steel), 1 percent
r,3,1,.04,0,0 ! mat 1 (steel), 4 percent
block,,100,,5,,5
block,,100,5,10,,5
block,,100,10,15,,5
!vovlap,all
NUMMRG,KP, , , ,LOW
numcmp,volu
esize,5
mat,2 ! concrete material
real,1 ! rebar
vmesh,1
real,2
vmesh,2
real,3
vmesh,3
nsel,s,loc,x
d,all,all
nsel,s,loc,x,100
d,all,all
nsel,all
fini
/solu
nsel,s,loc,y,15
sf,all,pres,100
nsel,all
OUTRES,ALL,ALL,
nsub,10
solve
展開 ANSYS鋼筋混凝土結構開裂計算介紹
1．概述
開裂計算是工程中比較關心的問題，但一直是有限元分析的一個難點，涉及到材料本構、計算收斂性等諸多問題。ANSYS＋CivilFEM提供了鋼筋混凝土結構開裂計算功能，其中土木專用模塊CivilFEM提供的非線性混凝土計算適用于混凝土梁結構的非線性計算（包括開裂），可以直接通過截面定義鋼筋，從而模擬鋼筋混凝土梁。
但對于更一般的結構，用梁單元來模擬不一定合適，需要采用更一般的單元，ANSYS提供了專用的鋼筋混凝土實體單元SOLID65來模擬鋼筋混凝土結構，該單元材料采用混凝土材料模型，可定義混凝土的開裂、壓碎準則。
另外可以定義鋼筋方向和體積率，可用來模擬鋼筋混凝土的破壞。本文將通過算例對ANSYS＋CivilFEM開裂計算的效果進行探討，并針對一些計算難點提出初步的解決方案。
2．CivilFEM開裂計算
CivilFEM適合于梁結構開裂分析，另外為了與后面SOLID65單元開裂計算結果進行比較，先探討了CivilFEM的開裂計算。
CivilFEM開裂計算需要考慮的要點：
1、激活CivilFEM非線性模塊（~CFACTIV,NLC,Y），這是CivilFEM非線性計算的前提。
2、即使事實上為小變形，也必須打開幾何非線性效應（NLGEOM,ON），否則無法激活非線性迭代。
3、通常應該關閉求解控制（SOLCONTROL,OFF），由于CivilFEM非線性計算通過修改實常數的等效方法，自動求解控制反而可能導致發散。
4、在收斂不好的情況下，可以增加子步數、打開自動步長（AUTOTS,ON）或可以給定一個比較大的迭代數（NEQIT,NUM）,以改善收斂，線性搜索有時也可以改善收斂（LNSRCH,ON）。
展開 基于ansys的梁單元、實體單元徐變精細化分析（含各參數解釋） ￥25
2、改網格模型，改成自己對應的網格模型，網格用ansys，hypermesh，ansa等前處理軟件都沒問題。
3、改材料參數，改成你想要的徐變模型，對著規范或者是你做出來的試驗擬合曲線。
以上即可實際應用。
ANSYS單元類型選擇方法 附ansys結構單元與材料應用手冊下載
六、單元類型選擇方法
7.進行完前面的選擇工作，單元類型就基本上已經定位在2-3種單元類型上了，接下來打開這幾種單元的幫助手冊，進行以下工作：
仔細閱讀其單元描述，檢查是否與分析問題的背景吻合、
了解單元所需輸入的參數、單元關鍵項和載荷考慮；
了解單元的輸出數據；
下載地址：ansys結構單元與材料應用手冊
ANSYS中桿單元和殼單元的單元耦合問題
在比較復雜的結構的有限元分析中，不同的結構部件通常使用不同類型的單元來模擬。
通常情況下，不同類型的單元的各個節點的自由度數目是不同的，不同類型單元的連接節點處的自由度的耦合問題，是一個比較令人頭疼的問題。
在ANSYS中通?？梢杂民詈厦頒P來耦合不同類型單元在連接節點處的自由度(DOF)。
也可以用CE命令來認為添加自由度之間的約束方程來達到耦合的目的。
下面是一個簡單的算例，使用了CE命令來耦合連接節點處的自由度。
模型是航天器的機翼的一個Section的某一個隔框。上下表皮是薄殼結構，用Shell63單元來模擬，在上下表皮之間有起支撐作用的桿件，用link8單元來模擬。
建模的時候，link8單元和shell63單元在連接有各自獨立的節點。即：link8單元和shell63單元的節點在連接處是重合的，但是，節點編號是各自獨立的。
link8單元在每個節點有 ux,uy,uz3個平動自由度；
shell63在每個節點有ux,uy,uz這3個平動自由度和rotx,roty,rotz這3個轉個自由，共6個自由度。
在耦合節點處，兩個耦合節點的ux,uy,uz自由度應該是相等的。
這個等式可以用CE命令來描述。
完整的命令流如下：
finish
/clear,start
/prep7
!定義第一種材料屬性；
mp,ex,1,30e6
mp,prxy,1,0.3
!定義shell63單元和實常數；
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, ,
在關鍵點處生成節點；
nkpt,100,4 !與編號為117的節點耦合
nkpt,101,9 !
展開 ANSYS各類型單元連接專題講解(五)之3D梁單元與殼單元剛接
例如采用ANSYS模擬一個多層混凝土框架結構，一般除計算整體指標外，我們在計算具體荷載作用時（如風荷載、地震作用、恒載、活載等），樓板一般采用彈性版，此時可用殼單元模擬，主梁、次梁采用梁單元模擬，此時變為梁單元包含在殼面內的情況，當然此類情況是否需要考慮截面偏置，可根據具體工程而定。
對這中梁單元包含在殼單元面內的情況，只需要將梁單元與殼單元共用節點即可，而無須格外建立約束方程。
三、梁單元在殼單元內但不包含
此種情況為梁與殼位于同一面內，但其中面不包含梁線，適用于多尺度建模分析（如下圖）。梁單元與殼單元的連接在端部可以通過剛性梁和剛性區域兩種方式連接。剛性梁采用MPC184單元，剛性區域采用Cerig命令，具體使用方法下期文章討論。
展開 
ANSYS APDL實體單元和殼單元（不共節點）之間的連接 ￥100
實體單元和殼單元之間的連接是ANSYS中常見的問題。即使兩種單元之間共節點，但單元之間不連續（實體單元每個節點有3個平動自由度，而殼單元每個節點有3個平動自由度和3個轉動自由度），對于兩種單元之間面面接觸，可直接定義剛域，本文主要采用MPC法對實體-殼單元的連接方法進行說明。
1 單元類型
算例模型中，實體單元采用SOLID45，殼單元采用SHELL63，接觸位置不共節點。對于兩種單元之間的連接，通過目標單元TARGE170和接觸單元CONTA175實現，定義約束為實體-殼約束，接觸單元為MPC算法，接觸類型為綁定接觸。
2 有限元模型和綁定接觸
圖1 底部固定約束，殼單元施加均布荷載
圖2 目標單元和接觸單元
3 計算結果
圖3 von Mises stress
圖4 X-Component of displacement
付費內容為相關命令流。
展開 ANSYS梁單元與實體單元的耦合與約束方程
ANSYS梁單元與實體單元的耦合與約束方程
By長安CAE
1 概述
在ANSYS計算過程中，有時候會遇到不同單元之間進行連接，由于不同的單元自由度不同，連接時通常需要通過耦合和約束方程建立節點自由度的關系，保證結果的準確性。
耦合可以理解成是將耦合的對象某個自由度作相等處理，而約束方程則不局限于相等這個關系，其可以描述具有某種關系的自由度。如圖1所示，為梁單元與平面單元的連接。如果不采用約束方程，力矩的傳遞無法完成，因為平面單元沒有轉動自由度。
圖1 梁單元與平面單元連接
為使節點2具有力矩傳遞的能力，要求1、2、3節點之間的自由度滿足以下關系：
ROTZ2 = (UY3 - UY1)/10
再通過CE命令，即可將此關系通過約束方程的形式施加給1、2、3節點。
2 命令
查看ANSYS的幫助文檔，查詢CE命令的解釋，如圖2所示。
圖2 ANSYS的CE命令解釋
CE, NEQN, CONST, NODE1, Lab1, C1, NODE2, Lab2, C2, NODE3, Lab3, C3
其中，NEQT表示常數，用于區別約束方程，一般可以用數字1、2、3表示即可，表示第幾個約束方程；
CONST表示方程的常數項，一般為0；
NODE1，表示第一個節點；
Lab1，表示自由度標簽，對于結構而言，就是三個平移和三個轉動自由度；
C1，表示該自由度的系數；
同理，后面的也一樣。
展開 在ANSYS 中3維坐標下的 shell structure 使用2D 平面單元劃分，應該使用哪個單元型號的單元
在ANSYS 中3維坐標下的 shell structure 使用2D 平面單元（僅考慮平面內的位移）劃分，應該使用哪個單元型號的單元？
ANSYS特殊單元——Follw201（隨動荷載）單元
ANSYS的Follw201單元是ANSYS的幾個特殊單元（比如mesh200）之一，稱為隨動荷載單元。都知道在ansys里面施加壓力載荷pressure時，其實載荷是可以隨動的，也就是能夠一直保持著面的法線方向，而施加集中了或者力矩時則不能保證。Follw201單元便能解決這個問題。
Follw201單元是一個單節點的3D單元，具有六個自由度，只能夠覆蓋在既有單元節點上，而且節點必須具有3個平移自由度和3個轉動自由度，也即是只能用在梁單元和殼單元上，實體單元僅有三個自由度。
Follw201單元主要用于幾何非線性分析問題中，在這類問題的分析過程中幾何會發生比較大的變形，面或線的法線方向可能發生比較大的變化，施加的載荷的方向是否隨動對結果的影響非常大。也就是用到此單元時會配對使用Nlgeon,on命令以打開大變形開關。如下圖所示為單元示意圖。
圖1
每個單元有兩個面，面1用于設定集中力的大小，面2用于設置力矩的大小，面的方向在應用時是通過單元的實常數進行定義的。
另外還需要注意，有限元求解的時候大部分是求解對稱矩陣，但隨動荷載單元的應用則包括了隨動荷載剛化效應，使剛度矩陣為非對稱的，因此需要采用非對稱求解器進行計算。
下面是具體應用，建一根梁單元，在梁的端部施加隨動集中力。
/prep7
!定義參數
EE=207E3
B=10
LCD=300
AA=B*B
IZ=B**4/12
PHZ=EE*IZ/LCD/LCD
!定義單元和材料
!201單元不需要定義材料
et,1,beam4
et,2,follw201
mp,ex,1,ee
mp,prxy,1,0.3
!定義實常數，實常數1設置梁單元的參數
r,1,aa,iz,iz,b,b
r,2,,1.0
!
展開 