abaqus壓痕損傷的案例
Abaqus用戶子程序USDFLD實例詳解-復合材料層壓板漸進損傷強度
Abaqus用戶子程序USDFLD實例詳解--復合材料層壓板漸進損傷強度分析
1、用USDFLD子程序漸進損傷分析方法分析層壓板的強度。
2、層壓板由25層單層板組成，單層板厚度為0.15mm。層壓板的鋪層順序、單層板的材料屬性見下表。
3、層壓板幾何尺寸：
厚度 T = 3.75mm;
寬度 W = 12mm；
長度 L >1.5W。
L可以任取，主要研究材料性能，為研究大部件性能做參考。
4、加載方式：沿長度方向施加壓縮載荷，在本例中施加位移載荷，讀取相應的力。
詳細步驟：
code.zip
Abaqus用戶子程序USDFLD實例詳解.rar
展開 Abaqus考慮拉壓不對稱的樹脂彈塑性損傷本構vumat子程序開發
本文通過在屈服準則中引入拉壓非對稱參量，研究了樹脂的拉壓不對稱彈塑性損傷行為。
由于樹脂的屈服行為與靜水壓力相關，這里采用下式所示的拋物面屈服準則。
式中J2為偏應力的第二不變量，I1為應力第一不變量，σt和σc為拉壓屈服應力
采用非關聯塑性流動準則，如下所示。
式中σvm為mises等效應力，P為靜水壓力，α為材料參數
損傷萌生準則如下所示
式中J2和I1為無損應力下的不變量。
為了降低模型的網格依賴性，損傷演化采用特征長度相關的指數模型
式中，rm為損傷內變量，am為特征長度相關的材料參數。
Melro的文章中給出了通過Simpson積分和弦截法計算Am的方法，實際計算發現通過該方法計算的am效果不是太理想，因此本文未對am進行迭代，直接采用其初值進行仿真計算，如下所示。
計算流程如下
計算流程圖
根據上文的彈塑性損傷模型編寫了vumat子程序，并通過單胞模型進行了驗證，計算結果如下圖所示。
abaqus單胞模型
拉伸載荷下的應力應變曲線
壓縮載荷下的應力應變曲線
展開 Abaqus通過VUMAT子程序實現混凝土拉壓不對稱彈塑性損傷本構模型
在考慮混凝士等準脆性材料的非彈性力學行為方面，連續損傷力學模型可以通過不同的方式來描述材料剛度和強度的退化以及單邊效應。真正意義上的彈塑性損傷本構模型：不僅考慮卸載時不可恢復塑性變形的影響，而且還應該考慮損傷和塑性的雙向耦合效應。
彈性階段應力應變滿足如下關系
通過對應力進行譜分解，可得
式中，σ 為名義應力，d 為損傷，d=1-exp(-εp/ρ0)，公式右端σ為有效應力的正負分解。
拉壓屈服函數如下所示
屈服后，塑性流動由下式定義
按照彈性預測-塑性修正-損傷修正的流程，通過在主應力空間進行譜分解，結合徑向返回算法，本文編寫了混凝土彈塑性損傷的VUMAT子程序。
通過對單胞的單向拉壓模擬可以計算得到混凝土的應力應變響應如下圖所示。
拉伸損傷演化過程
壓縮損傷演化過程
不同圍壓下的應力應變曲線
可以發現，隨著圍壓增大，混凝土壓縮強度提高
展開 Abaqus通過VUMAT子程序實現混凝土拉壓不對稱彈塑性損傷本構模型
在考慮混凝士等準脆性材料的非彈性力學行為方面，連續損傷力學模型可以通過不同的方式來描述材料剛度和強度的退化以及單邊效應。真正意義上的彈塑性損傷本構模型：不僅考慮卸載時不可恢復塑性變形的影響，而且還應該考慮損傷和塑性的雙向耦合效應。
彈性階段應力應變滿足如下關系
通過對應力進行譜分解，可得
式中，σ 為名義應力，d 為損傷，d=1-exp(-εp/ρ0)，公式右端σ為有效應力的正負分解。
拉壓屈服函數如下所示
屈服后，塑性流動由下式定義
按照彈性預測-塑性修正-損傷修正的流程，通過在主應力空間進行譜分解，結合徑向返回算法，本文編寫了混凝土彈塑性損傷的VUMAT子程序。
通過對單胞的單向拉壓模擬可以計算得到混凝土的應力應變響應如下圖所示。
拉伸損傷演化過程
壓縮損傷演化過程
不同圍壓下的應力應變曲線
可以發現，隨著圍壓增大，混凝土壓縮強度提高
最后歡迎通過公眾號"320科技工作室"聯系我們.
展開 
