ansys內聚力單元劃分的案例
任意單元間插入零厚度內聚力單元的新腳本
腳本使用方法
（1）首先創建一個要插入內聚力的零件，劃分網格
（2）將該零件添加至裝配中（此步驟不能遺漏）
（3）寫出inp文件：創建計算任務名，例如Job-1，則在本地寫出Job-1.inp
（4）為了防止再被各種盜版，該腳本添加了license和時間、單元數量限制，運行之前需要將名為license.key的許可證復制到abaqus當前工作中。
（5）在Abaqus中運行腳本，會提示輸入inp文件名，如下圖所示。待處理的inp文件須放入工作目錄下。
程序執行完成以后，會在本地生成名為“cohesive.inp”的文件，里面是單元和節點信息。并且將原始網格單元和新生成的內聚力單元分別創建了單元集合，方便賦予屬性。
同時為了方便賦予屬性，程序自動將inp導入到了Abaqus CAE中，并創建了新的model、part和單元集。
零厚度內聚力單元生成后，就可以建立完整的有限元模型去做隨機裂紋擴展分析了，下面是一個脆性材料開裂的案例。
展開 【CAE案例】熱裂紋仿真：引入內聚力單元以及應用前景
作為一個例子，當溫度為1400°C時的內聚力法則的參數k0=10-5，σc=22Mpa，δc=0.04mm
內聚力法（通過對試驗的調整來確定）允許在有裂紋和非無裂紋試樣之間有一個明顯的分界。
05 結論
本研究通過對熱裂紋產生原因的分析，在SPAR項目框架內展開宏觀建模，并使用通用結構仿真軟件中的熱力學計算進行測試，最后引入內聚力單元以更深一步對熱裂紋進行分析和研究。
更多資訊可登錄格物CAE官方網站
https://cae.yuansuan.cn/
遠算科技在bilibili、頭條、知乎、技術鄰定期發布課程視頻等內容
或關注微信公眾號遠算云仿真
敬請關注
展開 如何提高運用內聚力單元時的收斂性
最早是看高雁飛老師的一篇文章 cohesive，2004 他引入了一個一個粘性系數，來提高收斂性。這也是abaqus 內置的 viscous regularization的方法。也就是第一種方法（1) viscous regularization
(2) automatic stabilization
也就是施加自動材料或者模型阻尼
（3） non-default solution controls.
詳細的介紹大家可以看看 advance in numerical model in adhesive joints 和高老師文章.
Abaqus批量插入內聚力單元插件源碼分享
之前在Simwe論壇上分享過一個Abaqus中批量插入內聚力單元的插件http://forum.simwe.com/thread-1287781-1-1.html，這里把當時的源碼分享給大家，功能比較單一，有需要的可以下載下來自行修改及使用。
插入內聚力單元的效果
注意：目前只支持實體的六面體單元，生成的Cohesive單元為零厚度。下載的插件解壓后放到abaqus工作文件夾的abaqus_plugins文件夾中，重啟abaqus就能在工具欄plugins下面看到Cohesive按鈕。
這個插件的運行速度還算比較快，像包含256000個單元的球體，生成COHESIVE單元只需要60幾秒。幾萬個單元的話，幾秒中就能生成完畢。
V3.2.rar
Insert_Coh.zip
有關于Abaqus、LS-DYNA等子程序開發，復合材料仿真的問題可以聯系QQ1653004885或者關注CAE320公眾號
展開 
abaqus 2020 熱力耦合內聚力單元初探
最近嘗試了一下abaqus 2020 熱力耦合內聚力單元，給的是零厚度。但是，好像和我自己編的子程序存在差異，用零厚度熱力耦合內聚力單元（coh2d4T）結果也不是很合理，不知道是什么原因。沒有看到相應幫助文檔。好像不對稱，邊界是對稱的，挺奇怪的。模型中間是零厚度的內聚力。
我的子程序結果
具體的內聚力傳熱可以看看我們的文章：
Analysis of delamination and heat conductivity of epoxy
impregnated pancake coils using a cohesive zone model
展開 內聚力單元二次開發子程序
支持二維，三維(四面體和六面體)。
自研批量生成0厚度內聚力單元插件 ￥400
該插件基于python語言批量生成內聚力單元，使用該插件需要配置python語言環境。使用者可隨時反饋使用過程中遇到的任何問題，發送至郵箱shenz1hao@126.com，承諾修正一切問題。
插件功能具體包含：
（1）三種插入方式：兩種材料內部+邊界、僅邊界、邊界+一種材料內部
（2）任意單元形式：三角形單元、四邊形單元、六面形單元、楔形單元、四面體單元以及這些單元之間的任意組合
（3）單元數量及INP生成方式無限制
abaqus 內聚力單元插件 ，不限制網格數量 ￥5
[圖片]
考慮三軸度影響的內聚力單元本構
考慮三軸度影響的內聚力單元本構
內聚力單元被廣泛應用于模擬各種斷裂行為。本文采用三軸度相關的內聚力本構模型模擬了材料的斷裂失效行為。有效三軸度參數
其中
αc為雙軸度比的臨界值
其中
三軸度相關的內聚力模型本構由三部分組成：線彈性階段，彈塑性階段及軟化階段，如下式所示。
其中
根據上述本構，編寫三軸度相關內聚力模型UMAT子程序，通過單胞計算得到的結果如圖所示。
在進行真實試件計算時，需要將相鄰實體單元的三軸度參數傳遞給內聚力單元，如圖所示。
Abaqus中具體實施時，通過UVARM子程序獲得當前實體單元的應力狀態，并計算得到對應的三軸度，然后通過公共變量實現UVARM和UMAT之間的數據交換，將實體單元的相關參數傳遞給對應的Cohesive單元。計算得到的結果如下圖所示。
展開 基于hypermesh零厚度內聚力單元創建3D！ ￥50
背景: 內聚力單元的使用越來越頻繁,內聚力單元也廣泛的應用于脆性材料甚至塑性材料的損傷破壞問題的有限元描述!那么如果用hm來構建內聚力單元,特別是零厚度的內聚力單元呢? 如果是用過ABAQUS的朋友可能對這個問題不屑一顧，的確ABAQUS是我目前用過在內聚力單元建立方面最好的通用前處理器軟件，特別是對于零厚度單元的構建操作極其簡易。但是畢竟并不是所有的人都用abaqus，其中很多人用慣了hm，希望能直接在hm中構建cohesive zone element， 下面就講下在hm中零厚度內聚力單元的構建方法。關鍵字： hypermesh， 零厚度內聚力單元
目標：
1.實現平板零厚度內聚力單元的構建,板子上下厚0.1m，中間插入一層零厚度單元。
2.曲面零厚度內聚力單元的構建
在如下兩者構型中中插入零厚度內聚力單元（ 1種平板結構，1種曲面體結構）
需要分如下兩步操作：
Step 1先生成含有厚度內聚力單元模型
step 2 偏置有厚度內聚力單元生成零厚度內聚力單元模型
如下（隱去另一側單元便于顯示）
下面是step by step的教程和hm的練習模型
hm的版本是2017
展開 【ABAQUS建模】內聚力單元模擬復合材料分層（附cae文件）
創建分層網格：根據復合材料的分層結構，使用ABAQUS提供的網格劃分工具創建相應的分層網格。確保每個層級都被適當地劃分，并且層間接觸良好。
定義單元類型：根據復合材料的性質，選擇適當的單元類型。對于復合材料，常用的單元類型包括二維殼單元（例如S4R、S8R）和三維實體單元（例如C3D8）等。確保所選的單元類型適合您的分析目的和模型幾何。
定義內聚力模型：對于復合材料的分層界面，可使用ABAQUS中的內聚力模型來模擬分層的粘合特性。選擇適當的內聚力模型（例如表面內聚力模型或體積內聚力模型），并設置相關的參數，如強度、剛度和失效準則等。
施加邊界條件和加載：根據您的分析需求，在模型中定義適當的邊界條件和加載。這包括約束邊界條件、施加的載荷或位移等。確保邊界條件和加載方式與實際情況相符。
設置分析步驟和求解器選項：在ABAQUS中設置適當的分析步驟和求解器選項，以便執行所需的分析。這包括選擇合適的加載步驟、求解器類型和收斂準則等。
注意事項：
確保幾何模型的準確性，包括分層結構的幾何形狀和尺寸。
選擇適當的材料模型，確保準確地描述復合材料的各向異性性質。
注意分層網格的質量，避免扭曲或退化的單元。
確保內聚力模型的選擇合理，并根據實際情況設置參數。
在施加邊界條件和加載時，仔細考慮所需的約束和實際載荷情況。
檢查模型的收斂性和穩定性，調整求解器選項以獲得準確和可靠的結果。
展開 
一直特別好奇為什么有人寫UMAT的內聚力單元
個人覺得內聚力單元是一類特殊的單元，其形函數和求的位移量并不匹配。為什么有人寫UMAT的內聚力單元還能成功？是我理解的局限還是有人亂來？算了，以后再說，目前也不感興趣了
ABAQUS讀懂用戶手冊系列—修煉Cohesive內功：內聚力單元/接觸基礎知識點
</p><p class="ql-align-center"><strong>02</strong>內聚力單元的易錯點——網格行為</p><p>根據本人項目經歷與同學們反饋的問題，喵星人總結了使用內聚力單元時可能出現的易錯點：</p><p>1、由于內聚力單元必須劃分或指派粘結層網格，為了識別粘結層方向，必須掃掠指派厚度方向，若未正常掃略網格，則按照網格節點序號進行識別厚度方向，這將使得內聚力粘結方向可能存在問題。</p><p>2、被粘結的部件接觸面需與粘結層綁定or來自同一部件。部分同學會在設置內聚力單元后，又把粘結層和被粘結物體間設置了一層內聚力接觸，這樣相當于串聯了兩次粘結行為，屬于概念型錯誤。</p><p class="ql-align-center"><strong>03</strong>內聚力接觸的易錯點——初始化過程</p><p>這個過程ABAQUS用戶手冊中特別指出容易弄錯：</p><p>The initial contact status as a function of position along a cohesive contact interface can fundamentally affect simulation results. The intent for this example is that the block is initially touching the wall with the cohesive status initialized to bonded.
展開 LS-DYNA有厚度/零厚度內聚力單元建模案例（共節點法） ￥9.9
本案例為有厚度/零厚度內聚力單元建模示范，操作軟件為LS-DYNA自帶的前處理器Ls-Prepost。
效果如下圖所示：
無厚度
有厚度
付費內容為：
操作過程PDF文件
步驟詳細，演示形式為靜態圖片與動態圖片；給出了詳細的內聚力單元層的材料參數及出處。
2. 對應的K文件及計算結果
展開 結構膠內聚力單元（cohesive element）本構在ABAQUS和LS-DYNA中的應用 ￥5.99
本文以一簡單的膠粘兩塊板的模型說明內聚力單元在ABAQUS和LS-DYNA軟件中的應用。
ABAQUS采用的隱式計算，在顯式計算中的設置與隱式計算的設置相同，LS-DYNA采用的顯式計算，由于給的時間較短，可以看出板有明顯的抖動，僅供交流學習，感興趣的可以下載后邊的.inp和.k文件
