Hashin損傷建模的案例
平紋復合材料VUMAT子程序本構介紹（hashin準則及線性損傷演化方法）
因此在建立均質化模型時，平紋復合材料的剛度矩陣，損傷起始準則，損傷演化方法以及退化的剛度矩陣與單向復合材料具有明顯的差異。主要體現為平紋復合材料在面內的兩個方向均有纖維，因為對于平紋復合材料的失效模式主要有：經向拉伸/壓縮損傷，緯向拉伸/壓縮損傷以及厚度方向上的拉伸/壓縮損傷,此外還可以通過在層間插入cohesive單元考慮層間分層失效。接下來主要介紹層內的損傷本構關系。
1. 平紋復合材料損傷剛度矩陣
（1）
（2）
（3）
其中，d代表損傷系數，L、T以及Z代表三個垂直的方向，t、c代表拉伸，壓縮損傷，例如dlt代表縱向拉伸損傷。
2. 損傷起始準則（hashin準則）
(4)
其中，f1t,f1c代表縱向纖維拉伸和壓縮損傷，f2t,f2c代表橫向纖維拉伸和壓縮損傷，f3代表厚度方向上的失效，其中應變大于0時為拉伸失效，小于0時為壓縮失效。
3. 損傷后損傷演化模型（線性退化模型）
（5）
當滿足損傷起始準則后，損傷演化開始起作用。其中1t、1c、2t、2c、3t、3c的失效模式下對應的損傷系數分別為d1t、d1c、d2t、d2c、d3t、d3c。其中，εii為當前應變，εiimax為初始失效應變，εif為最終失效應變。
展開 隧道光面爆破局部建模損傷分析 ￥50
1.CAD中對隧道爆破炮孔布置圖進行輔助線切割，然后通過REG命令生成封閉面域
2.導出iges格式，將文件導入ANSYS/APDL軟件中進行巖石區域建模及網格劃分。鉆孔區域采用映射網格劃分，鉆孔外巖石區域采用掃掠劃分方式，單元類型為solid164單元，模型厚度方向擴展200cm，采用三維建模方法進行分析。模型網格劃分好后導出k文件，后續操作通過k文件導入ls-prepost中進行炸藥，堵塞及空氣的分區及材料參數、邊界條件、求解等設置。
3.巖石采用RHT模型，炸藥模型中采用*MAT_HIGH_EXPLOSIVE_BURN 材料模型，空氣材料采用*MAT-NULL 材料模型描述，空氣的狀態方程采用*EOS- LINEAR-POLYNOMAIAL 描述。
4.計算結果如下：
損傷破壞圖
展開 復合材料沖擊損傷自動建模程序
composite impact.rar
初學python做了一個復合材料沖擊損傷的自動建模程序，
程序使用說明如下：
1.程序還有待完善。目前只支持矩形板，矩形網格。
2.可以選擇3D實體單元及2D連續殼單元，但是如果用3D實體單元需要自己編寫VUMAT，這里就不上傳VUMAT了，大家自己努力編吧。
3.自動生成的有限元模型是沒有邊界條件的。自己需要手動添加邊界條件
4.輸出變量需要自己去設置。
5.需要先創建材料或者從自己的材料庫中導入材料。
6.鋪層參數可以從文本文件直接導入，在鋪層表格出點擊右鍵會出現導入對話框。
7.可以根據需要選擇是否增加cohesive elements。
補充說明：第一次上傳的附件有一個問題，就是用CompositeLayup直接生成實體單元是不能用在abaqus/explicit中的，提交計算會出錯，現在已經將程序完善，附件已經更新。
直接將壓縮包解壓縮后放到plugin文件夾下即可。
composite_imoact_en.rar
展開 轉載，復合材料分析中Cohesive單元建模及損傷簡介
為了更好的理解K，我們把寫成：
這里我們用來代替1，其中L可以理解為建模厚度，即建模時cohesive interface的幾何厚度；為實際厚度，即cohesive interface的真實厚度，這個厚度在cohesive section中定義。可以理解為幾何剛度，即模型中cohesive interface所具有的剛度；為cohesive interface的真實剛度。當為1時，計算界面剛度就采用幾何剛度，當為0.001時，計算時界面剛度變為1000。舉個小例子，如果界面的實際厚度為0.01，而在建模時就是按照這個厚度建立的，在定義material-section時又specify這層的厚度為0.01，實際上就等于把界面剛度提高了2個數量級，模擬結果當然是不對的，這時定義section時應采用默認厚度1。
ABAQUS在cohesive建模中使用了很“人性化”的設計，實際問題中界面可能很薄，有的只有0.001mm，甚至更小。有些問題cohesive單元的interface還可能是0厚度（比如crack問題），而相對來說整體模型也許很大，如果不引入這兩個厚度，我們就要在很大的模型中去創建這個很小的界面這是一個很麻煩的事情。引入這兩個厚度，在建模時我們就可以用有限的厚度來代替這個很小的界面厚度，只要在section中定義這個就好了。
進入property界面，點擊Material→Creat，在彈出的Edit Material對話框中，可以編輯新創建的cohesive材料的名稱，然后點擊Mechanical→Elasticity→Elastic→Traction，在空格中輸入相應的剛度。
2.1.損傷準則
初始損傷準則
初始損傷對應于材料開始退化，當應力或應變滿足于定義的初始臨界損傷準則，則此時退化開始。
展開 
論文帖一個：samcef層壓復合材料損傷建模
DAMAGE MODELING OF LAMINATED COMPOSITES: VALIDATION OF THE INTRA_LAMINAR LAW OF SAMCEF AT THE COUPON LEVEL FOR UD PLIES
文章中，使用samcef對層壓復合材料進行建模，并通過與實驗測試結果的對比驗證了模型，對不同方向的單層復合材料進行了深入研究，對參數辨識方法進行了一定討論。結果認為僅僅通過少量的測試獲取參數，就能通過建模仿真預測出復合材料的特性。
DAMAGE MODELING OF LAMINATED COMPOSITES VALIDATION.pdf
展開 西北工業大學《Composite Structures》:評價固化殘余應力對復合材料損傷行為影響的建模框架
2.2 損傷失效模型
由于單向復合材料的橫向損傷行為不受纖維失效的影響，因此主要研究基體和界面的失效。樹脂基體在載荷作用下依次表現出彈性、塑性和損傷，使用擴展線性Drucker-Prager定義其屈服行為，采用Ductile韌性準則描述其損傷行為。界面采用零厚度Cohesive單元進行建模，使用雙線性本構模型進行表征。
圖3 基體本構方程
圖4 界面本構方程
2.3 固化過程-損傷失效一體化分析框架
通過復合材料多尺度熱-化-力耦合建模分析，得到了固化過程中的細觀殘余應力，為進一步研究細觀殘余應力對復合材料力學行為的影響，將殘余應力作為初始預定義場引入復合材料的橫向加載過程中進行損傷預測，建立了如下所示的固化過程-損傷失效一體化分析框架。
展開 