斷裂模擬ansys單元的案例
粘結單元cohesive 模擬斷裂
通過建立零厚度粘聚單面cohesive Element 模擬混凝土或復合材料層合板開裂：
實例可詳見：
http://forum.simwe.com/thread-1134530-1-1.html
交流郵箱：1551647104@qq.com；
ABAQUS模擬梁單元斷裂的本構方程（VUMAT） ￥3
* TEMPNEW, STRETCHNEW, DEFGRADNEW, FIELDNEW,
C WRITE ONLY -
* STRESSNEW, STATENEW, ENERINTERNNEW, ENERINELASNEW )
C
INCLUDE 'VABA_PARAM.INC'
C
DIMENSION COORDMP(NBLOCK,*), CHARLENGTH(NBLOCK), PROPS(NPROPS),
1 DENSITY(NBLOCK), STRAININC(NBLOCK,NDIR+NSHR),
2 RELSPININC(NBLOCK,NSHR), TEMPOLD(NBLOCK),
3 STRETCHOLD(NBLOCK,NDIR+NSHR),
4 DEFGRADOLD(NBLOCK,NDIR+NSHR+NSHR),
5 FIELDOLD(NBLOCK,NFIELDV), STRESSOLD(NBLOCK,NDIR+NSHR),
6 STATEOLD(NBLOCK,NSTATEV), ENERINTERNOLD(NBLOCK),
7 ENERINELASOLD(NBLOCK), TEMPNEW(NBLOCK),
8 STRETCHNEW(NBLOCK,NDIR+NSHR),
9 DEFGRADNEW(NBLOCK,NDIR+NSHR+NSHR),
1 FIELDNEW(NBLOCK,NFIELDV),
2 STRESSNEW(NBLOCK,NDIR+NSHR), STATENEW(NBLOCK,NSTATEV),
3 ENERINTERNNEW
展開 單元刪除模擬斷裂失效實例-----小球沖擊鋁板
很多人對單元刪除實現斷裂的模擬有認識誤區，以為定義了塑性就能算斷裂。 不否認從初學者的角度來說，參數太多，很容易讓人手足無措，無從下手。
但是斷裂問題本來就是一個非常復雜的問題，控制的材料參數非常多。所以，如果你真的要做斷裂分析，一定要好好看看manual里面progressive damge部分。有些軟件可能經過簡單的定義就能算出結果，但是算出結果是否能用？是否得到科研界和業界人士承認？ 這些年abaqus在科研界的廣泛使用，其實能說明一些問題。
這里提供一個鋁板被沖擊的實例，以前與實驗對比過，效果還可以，現在提供inp
小球沖擊模型
幾個常見需要注意的問題：
1）本例的小球為了得到等效質量（實際的沖擊體只有接觸的沖頭是圓形的，后面是柱體），沒有采用本身密度。 并其后設置為剛體。
2）property
材料參數是難點。很多人為材料參數哪兒來？ 唯一靠譜的途徑是自己做實驗，但是大家都知道實驗費用不菲而且還不一定有實驗條件。所以國內很多simulator其實是找老外做的實驗結果。
對于絕大部分常見材料，沖擊模擬所需的材料參數，baidu or google 材料編號+Johnson-Cook 都可以搜到（不過本版還是很多人問某某鋼的材料參數，搜索能力啊）
（這個inp里面材料參數來自abaqus manual）另外，實在找不到的材料參數去ls-dyna版逛逛，會有收獲的。
3 step--time：動態計算，時間是物理時間。很多人問time設為多少比較合適。現在流行的各種沖擊工況，幾乎都有現成方法或者經驗公式來估算沖擊歷時。
如果你實在是沒有任何經驗和參考數據，可以看看我的經驗。
就梁和板的沖擊而言，我個人經驗：
低速沖擊無擊穿，時間可以用pi*[sqrt（M1+M2/4）/kb] 近似估算.
展開 單元刪除模擬斷裂失效實例-----小球沖擊鋁板
附件：
1 inp文件
b.rar
2 載荷-位移曲線
3 動畫
這個中面斷裂圖（注意：是后處理顯示剖面，并非在前處理取半結構分析）
----------------------------------------------------分割線-----------------------------------------------------
aba版用單元刪除做斷裂分析的帖子小結：
2.1.16 Progressive failure analysis of thin-wall aluminum extrusion under quasi-static and dynamic loads
這個是manual里面經典實例，我附上manual里面模擬三點完全斷裂的inp：
three point bending_ alextrusion.inp
threepointbending_alextrusion.rar
---------------------------------------------------------------------------------------------------------------------------------------------
技術鄰推薦：
ANSYS SCDM基礎培訓視頻與教學模型
基于ANSYS的3D打印熱分析
cae論壇精華帖：CAE有限元分析介紹
展開 
ABAQUS模擬梁單元斷裂的本構程序（VUMAT） ￥5
完整的代碼并且有說明，尤其適合模擬納米材料斷裂（調整判據，不斷裂也行）的行為。根據自身需要，修改相應的參數即可。
混凝土細觀斷裂模擬-XFEM-擴展有限單元法
利用 隨機生成混凝土骨料，建立2D或者3D模型，采用XFEM模擬裂縫擴展
3D彎曲斷裂例子可見：
http://forum.simwe.com/thread-1134530-1-1.html
ANSYS Mechanical 2022 新功能：單元、接觸、斷裂力學、并行計算
本期是ANSYS Mechanical 2022 功能更新之單元、接觸、斷裂力學、并行計算。
文末領取學習資料
下面我們看看具體的更新內容：
一、單元部分
增強單元性能加強
面增強單元的彎曲剛度
使用單軸剛度單元進行反向求解
耦合單元的增強
運動副單元增強
二、接觸部分
基于Dual Shape函數的接觸算法
新的自適應小滑移選項
殼-實體組裝件的準確性改進
螺栓預緊支持通用軸對稱單元
網格獨立點焊增強功能
瞬態動力學精度改進：HHT算法
力矩收斂參考值計算穩健性改進
三、斷裂力學
基于應力比率的疲勞裂紋閉合
Paris定律與裂縫閉合效應相結合
應力比率（R）相關的疲勞裂紋擴展規律
靜態裂紋擴展的溫度/時間相關斷裂準則
自適應裂紋初始化/插入
3D界面單元
動態裂紋擴展尺寸控制
四、求解器效率提升
資源預測增強
分布式求解增強
文章篇幅有限
下圖微信掃碼領取完整版學習資料
展開 cohesive單元分層開膠斷裂模擬-雙懸臂梁剝離DCB（三維模型） ￥2
cohesive單元分層開膠斷裂模擬-雙懸臂梁剝離DCB（三維模型）
ABAQUS損傷斷裂(例2) 盾構機掘進巖石破碎模擬(單元刪除技術) ￥66.67
盾構機掘進時的巖石破碎模擬(含單元刪除技術)
采用顯示動力學:
盾構機刀盤模擬為剛體，在轉動掘進的過程中破碎巖石。在巖石達到破碎應力后，采用單元刪除技術刪除掉已破碎的巖石單元。
模型概況：
模型的建立及邊界條件的設置
模型的網格劃分
模擬的結果：
掘進時的動態效果
模擬后的結果
隱藏掉盾構機刀盤后巖石破碎時的應力分布
隱藏掉盾構機刀盤后巖石破碎時的等效塑性應變分布
Abaqus利用梁單元模擬螺栓連接 附基于ABAQUS對螺栓斷裂問題仿真分析下載
螺栓的模擬在Abaqus也有幾種不同的處理方式。
（1）建立三維實體的螺栓模型，包括螺紋結構；
（2）建立三維實體的螺栓模型，忽略螺紋結構；
（3）建立三維實體的螺栓模型，由Abaqus自帶的螺紋接觸定義方式設置螺紋接觸；
（4）利用梁單元或者桿單元模擬螺栓。
本次以梁單元模擬螺栓為例，簡單闡述其應用。利用梁單元模擬螺栓與實體螺栓相比優勢比較明顯，模型簡單、接觸定義簡單、收斂容易，同時梁單元也能有效反應螺栓的受力情況，在很多情況下比較適用。
螺栓的模擬通常需要考慮預緊力的作用，利用CAE方法模擬螺栓預緊力的過程主要由三個載荷步完成，下面的例子會涉及。
建立如下所示的模型，三個部件，兩塊板和一根梁，其中梁是一個3D wire，建立一條線即可。
圖1
材料屬性定義的時候，梁單元需要指定梁截面，如下圖所示。
圖2
梁的截面形狀可以根據需要指定，本次為圓形截面，半徑為10，如下圖所示。
圖3
同時，梁單元還需要指定方向，通過菜單欄Assign-Beam Section Orientation，給出其中的n1向量，這里注意，梁的軸向是由向量t表示的，n1和n2兩個向量決定梁截面，其中t向量和n1、n2兩個向量決定的平面垂直。
本次定義n1向量為0,0,-1，最終梁的方向定義完成如下所示。
圖4
之后利用Interaction模塊下面的Constraint將梁與相關位置建立MPC連接，如下所示。
圖5
梁單元的兩端節點分別與螺栓螺帽位置處的節點進行MPC連接，連接形式可以由多種，這里選擇Beam連接。
定義多個載荷步，其中前三個載荷步用于施加螺栓預緊力。
展開 基于ANSYS ls-dyna拉伸斷裂實驗模擬
基于ANSYS ls-dyna拉伸斷裂實驗模擬
作者：大龍貓 微信公眾號：CAE_ANSYS
拉伸斷裂實驗是測試材料的經典實驗，可以測量材料的應力應變曲線，測量材料的抗拉強度，作為經典的實驗如何獲取其模擬過程呢？仿真分析軟件AYSYS在默認的情況下，無論受力多大都不會被拉斷，其主要原因是算法的問題。
Ansys Workbench使用非線性彈簧單元模擬配合間隙 ￥10
問題：
工程中兩個零部件之間經常會有配合間隙，Ansys Workbench中可以使用combin39號非線性單元，通過控制不同行程的彈簧剛度來模擬間隙配合。
模型示例：
設定支座與軸有1mm的配合間隙，在一端施加X向100N作用力，查看運動位移。
計算步驟：
1. 在間隙配合位置，建立jiont連接，放開X向平動自由度。
2. 在間隙配合位置，建立spring連接，同時插入Commands 命令。
ET,_sid,39,0,0,0,1
R,_sid,0.95,1,1.05,10000
3. 查看計算結果，當運動至0.95mm后spring彈簧剛度值陡增限制了X向運動。
建議：
? 同一個連接區域不建議使用兩個重復的連接關系，即jiont連接和spring連接不要使用同一個區域。
? 本文對配合區域進行分段處理，中間為spring連接，兩側為jiont連接
? 使用Remote Point點創建連接，需要打開Beta選項。
? 這種等效方式并不能良好的反應間隙配合位置的應力狀態，需要校核配合區域的應力狀態還是需要使用接觸連接。
展開 ANSYS的生死單元模擬焊接過程
ANSYS的生死單元模擬焊接過程
1 概述
焊接模擬計算在CAE仿真是比較大的一塊內容，也是比較復雜的一個過程，幾個比較關鍵的問題是熱源函數的描述、單元的融覆、熱源的移動等等，通過單純的GUI操作，無論使ANSYS還是Abaqus都不大可能完成這個過程，通常需要借助軟件的內置語言。
本次主要介紹單元生死的應用，單元生死主要用于單元缺失的場合，比如凝固溶解過程，斷裂過程，焊接過程等等，這些過程都是非線性或者時間歷程過程，計算需要很多子步和迭代，為了在此過程中避免一遍一遍修改單元，便引入生死單元的概念，通俗的講就是通過一些方法讓單元失效，具體的改變是單元的彈性模量的改變，當單元死時，修改其彈性模量為非常小的值，讓其在求解過程中不起作用。
詳細地說，激活單元死這個狀態時，ANSYS程序將單元剛度矩陣乘以很小的因子，程序默認值為1E-6，死單元的單元載荷為0，從而不對載荷向量生效，同樣的，死單元的質量、阻尼、比熱等等參數也設置為0，單元的應力應變也因此為0。
2 前處理
前處理包括單元定義、材料定義和建模，單元定義是需要注意單元屬性，此次定義13號二維耦合單元，具有溫度和位移自由度。
材料屬性包括結構參數和熱參數，具體包含彈性模量，泊松比，屈服強度，塑性屬性，材料密度，熱膨脹系數，熱傳導系數，比熱容。焊接時溫度較高，定義材料通常需要定義多個溫度下的值。
展開 ANSYS生死單元之焊接過程模擬
在ansys計算過程中，如果需要向模型中加入（或刪除）實體，模型中對應實體部位的單元就“存在”（或消亡）。單元生死選項就用于在這種情況下殺死或重新激活選擇的單元。
例如，在焊接分析過程中，隨著高溫焊料的加入，坡口處的單元需要不斷地被激活；在材料斷料分析中，隨著裂紋的延伸，斷裂處的單元需要不斷的被殺死；在隧道挖掘和橋梁建立分析中，材料也需要不斷的被殺死或激活。因此，單元的生死應用技術廣泛的存在于ansys仿真分析中，是一項應用非常廣泛的技術。
單元的生死并不是ansys程序將殺死單元對應的實體從模型中刪除，或者激活重新生成材料，而是通過將其剛度矩陣，或者傳導矩陣（對應于不同的分析），乘以很小的因子（ESTIF），默認值為1E-6。死單元的單元載荷將為0，從而不對載荷向量生效，等效于將單元殺死；
同樣，當一個單元被重新激活時，其剛度，單元載荷等恢復其原始的數值，重新激活的單元也沒有應變記錄，在熱分析里面沒有熱量存儲。需要注意的是，生死單元對大部分單元可以應用，然而對某些單元卻是不可用的。
在一些情況下，單元生死狀態可以根據ansys的計算結果決定。如在斷裂分析中，我們需要將應力值大于材料屈服強度的單元殺死，可以利用Etable選擇相應的單元進行殺死，繼而返回到求解器進行求解，如果如此循環，則可觀察到裂紋的生長過程。
可以在大多數靜態和非線性瞬態分析中使用單元生死，其基本分析與相應的分析過程是一致的，主要包括三個步驟：建模，施加載荷并求解，查看結果。
現通過ansys焊接過程，講解生死單元的應用。
兩個平板進行對接，采用V型坡口。在焊接的過程中，焊料不斷加入坡口，進行焊接。平板溫度采用20℃，焊料溫度采用1500℃。
展開 ANSYS單元生死功能模擬門式剛架施工例子
ANSYS單元生死功能模擬門式剛架施工
! Simulate the construction of a frame with the element active/kill
! function of ANSYS
! 施工分為三步
! The construction is divided into 3 steps
! 1: 建立立柱和臨時支撐
! 1: Install the column and temporary support
! 2: 安裝橫梁
! 2: Install the beams
! 3: 去掉臨時支撐
! 3: Remove the temporary support
! 作者：陸新征，清華大學土木系
! Author: Lu Xinzheng Dept. Civil Engrg. of Tsinghua University
[Money=20]
FINISH
/CLEAR
/UNITS,SI
/PREP7
!*
SECTWIDTH=300 !構件截面寬度300MM
SECTHEIGHT=600 !構件截面高度600MM
SECTAREA=SECTWIDTH*SECTHEIGHT
SECTIYY=SECTWIDTH**3*SECTHEIGHT/12.
SECTIZZ=SECTWIDTH*SECTHEIGHT**3/12.
SPAN=24E3 !跨度24M
COLUMNHEIGHT=8E3 !柱子高度8M
SLOP=3E3 !
展開 