不過值得指出的是文中引入 cohesive 單元主要用于裂紋路徑的可視化表達，而其插入區域和參數設置并未像 GTN 參數那樣得到充分展開，因此這一部分更適合作為輔助性的裂紋表征手段，而非全文最核心的機理貢獻。 使用作者提出的完整積分框架，并基于顯式vumat實現，同時使用基于損傷變量的單元刪除方案同時引入ALE自適應網格方案可以更好的預測梯度效應。

若輸入角度數目少于設定的總層數，系統將自動以 0° 鋪層補齊所缺層信息，避免因輸入遺漏導致模型鋪層角度不完整。

而上、下壓板幾何特征為2D-Shell，采用S4R單元（4節點曲殼單元）對上、下壓板進行網格劃分，考慮到上、下壓板網格大小不影響計算的準確度，網格基本尺寸定義為15 mm×15 mm。 為了模擬CFRP圓管中各個鋪層之間的連接關系，通過定義通用接觸Cohesive Behavior和Damage參數來模擬CFRP圓管中鋪層之間的粘膠層。

由于這些限制，近20年來LS-DYNA研發團隊開發了一些新的方法，例如自適應法、XFEM、cohesive elements，以及Node splitting方法進行材料分析，這些方法適用于裂紋擴展問題分析，卻不適合于涉及到材料大面積損傷的加工過程分析，或者說不適用于材料加工過程中部分材料去除問題，基于有限元的方法都不太適用這類問題。

大家好，我在使用COHESIVE單元模擬晶體壓縮斷裂失效的時候，總是出現奇異值錯誤，COHESIVE單元不能按照設定的單元參數進行單元刪除。我想知道有什么方法可以解決這種錯誤？

更新內容： 修復：面集合嵌入時，部分位置錯誤的嵌入cohesive單元； 新增：原實體單元集合的邊界Cohesive單元創建單獨的OutCohe集合； 新增：平面二階四邊形和二階三角形單元邊界嵌入Cohesive單元。

如圖所示，工況一：兩個單元之間用一定厚度的cohesive單元相連，底部的單元固定不動，頂部的單元向右水平移動，兩單元發生純剪切破壞并導致cohesive單元被破壞刪除。 工況二：在工況一的基礎上在頂部單元的上表面施加壓強荷載使兩個單元相互擠壓，兩個單元之間發生壓剪破壞，cohesive單元即使達到失效標準也無法被刪除。

二、 Cohesive element原理 ABAQUSd提供一種粘結單元（cohesive element），用以模擬兩個部分之間的粘性連接，一般來說，它要求粘結材料尺寸和強度都小于粘結部分（比如多層復合材料的膠粘層），進而可以利用cohesive element模擬材料的斷裂。從本質上講，利用cohesive element模擬材料的斷裂其實是單元刪除方法的一種。

三、粘聚力模型的實踐與互動 以上介紹了Cohesive Element和Cohesive Surface的具體設置方式。作為學習型仿真工程師不能只看不練，否則只是過眼云煙。 接下來對幾個設置關鍵點進行闡述，讀者可以打開INP文件進行一一對照。

有些模版源碼可能沒有用上，手工去掉或保留都可以，不影響真正開發部分的源碼的執行。用戶可以在這個變量里加入自己的源程序文件。 MY_INC包含LS-DYNA用戶開發包的所有的頭文件，用戶可以添加自己的頭文件，但不建議刪除已有的頭文件。