abaqus 損傷演化的案例
關于內聚單元的損傷演化參數設定
<p>關于abaqus里cohesive單元損傷演化階段，我們可設定關于位移的線性、指數、自定義表格演化以及關于能量的線性、指數演化，其中線性與指數演化關系明確，abaqus幫助文檔中也有說明，但關于損傷D-位移E表格演化并沒有相應的幫助文檔，有時我們需要根據實驗數據確定這段演化關系，下面我將和大家分享如何確定損傷D與位移E的關系。</p><p><br></p><p><br></p><p><br></p><figure style="text-align: center;" class="ql-align-center">
<figure class="figure-image" contenteditable="false" data-img="https://img.jishulink.com/202503/attachment/fe893c991a46473683e10503eeeae1ce.png" style="display: inline-block;">
<img src="https://img.jishulink.com/202503/attachment/fe893c991a46473683e10503eeeae1ce.png" data-mobile-src="https://img.jishulink.com/202503/attachment/fe893c991a46473683e10503eeeae1ce.png?image_process=/format,webp" data-pc-src="https://img.jishulink.com/202503/attachment/fe893c991a46473683e10503eeeae1ce.png?
展開 abaqus vumat子程序損傷演化階段（應力更新）問題
我想編寫一個關于Johnson-Cook損傷的vumat子程序，然后做了一個單軸拉伸的案例，但是進入損傷階段后力和位移曲線震蕩嚴重
我想知道損傷演化階段應力如何更新，只需要改變彈塑性階段的剛度矩陣嗎
ABAQUS-纖維強化復合材料的損傷演化和網格移動
ABAQUS-纖維強化復合材料的損傷演化和網格移動.doc
ABAQUS VUMAT子程序 PUCK損傷起始準則+指數演化方法 ￥58
ABAQUS PUCK損傷起始準則+指數演化方法
復合材料本構
損傷矩陣
損傷應力應變關系
PUCK準則
指數退化方法
VUMAT代碼如下
ABAQUS損傷參數
無論斷裂還是切削等問題，考慮材料的損傷行為，基本上就是考慮損傷起始和損傷演化，這兩個參數基本上可以通過實驗數據計算而來。
對于單軸拉伸試驗，其典型的真應力-應變曲線包括了初始彈性階段，塑性階段，剛度下降階段和最終斷裂階段。最后兩個階段為損傷耗散過程，而對于材料損傷的定義也是基于這兩個階段的特征值進行的。
在ABAQUS中，材料損傷通常以一個損傷起始判據來定義材料的失效初始化。這個判據可以是材料在失效時的應力與應力狀態（基于不同的斷裂準則），也可以是損傷本構（如JC失效模型）。
在材料，或某一個單元發生了損傷后，此處的剛度會下降，而承擔的載荷（如應力）會隨之減少，并發生重新分配。則需要對材料在發生了失效時的演變過程進行定義，即損傷演變準則。對單軸拉伸而言，此階段與拉伸曲線在達到了最高值（抗拉極限或起裂應變處）之后的下降段所對應。單元的剛度沿損傷演化規律下降，最終完全失效，在分析中可以將其刪除（單元刪除法）或允許分離（黏聚力單元）或允許裂紋完全擴展（XFEM)。
實際分析中更關注失效參數的獲取，對簡單問題可以采用單軸拉伸的真應力應變曲線來計算，一個很好的例子如下：
Ductile_Damage_004.pdf
來自于木蟲上的一個問題
Abaqus損傷演化中的失效位移如何得到？ - 仿真模擬 - 小木蟲 - 學術 科研 互動社區 (muchong.com)
原文（感謝Ronald Heinz Norbert Wagner大佬）
How to perform Element deletion in ABAQUS using ductile damage criteria ? (researchgate.net)
展開 輕骨料混凝土細觀損傷演化分析
輕骨料混凝土細觀損傷演化分析
1研究現狀
混凝土材料的準脆性斷裂破壞是裂紋萌生、穩定擴展和失穩擴展的過程。為描述和分析混凝土的斷裂行為，1961年，Kaplan［1］開始將線彈性斷裂力學方法引入到了混凝土結構的分析中。由于混凝土存在應變軟化現象，混凝土出現裂縫后，骨料與砂漿間仍存在著齒合的粘結效應，且混凝土的微觀裂縫和亞臨界裂縫尺寸相較于一般的金屬材料大幾個量級，故彈塑性斷裂力學的COD和J積分理論亦不適用于描述混凝土在細觀尺度上的斷裂損傷［2］。
近年來在細觀尺度上通過數值模擬來分析和研究混凝土的宏觀斷裂逐漸成為熱點。目前國內外學者提出了許多用于模擬混凝土斷裂損傷過程的細觀力學模型，主要分為連續和非連續裂縫模型兩類［3-4］。經典連續裂縫模型包括彌散裂紋模型［5］、鈍裂縫帶模型［6］、旋轉裂縫模型［7］等。經典的非連續裂縫模型包括格構模型［8］、剛體彈簧模型［9］、隨機力學特征模型［10］等。這些裂縫模型均難以對混凝土內部裂縫的萌生、擴展以及貫通的全過程進行直觀展示，且模型網格敏感性較大，收斂性較差。內聚力模型因能克服上述缺陷，是目前模擬材料斷裂使用較多的模型。
輕骨料混凝土具有輕質高強、保溫隔熱性能好、抗火性好和抗震性能良好等諸多優點，目前已成為僅次于普通混凝土用量最大的一種新型混凝土，得到了國內外的廣泛應用，國內外對輕骨料混凝土的力學性能、耐久性能等開展了較為廣泛的試驗研究，而有關輕骨料混凝土細觀層次上斷裂損傷演化分析和數值模擬研究仍然很少。基于連續損傷力學的內聚力模型是采用在實體單元之間嵌入內聚力單元的方法來模擬損傷以及斷裂行為。目前采用內聚力模型來模擬輕骨料混凝土單軸拉壓下損傷斷裂行為的還未見相關報道。
展開 蘭大《Nature》子刊：超導帶材損傷演化檢測技術突破性進展！
依照磁光法的磁通圖像演化過程特征，進一步提出了對磁通運動與損傷機制認知研究的新途徑，取得的主要成果有：
1、發現了如圖3所示的一種點狀磁通運動的新模式及其運動時間尺度的變化特征，即隨著環境溫度從40K升高到60K和77K，運動時間尺度從 30ms增加至60ms和100ms；
2、隨著應變的增大，磁通運動由點狀結構變化為紡錘狀的穿透模式，其速度分布具有多級特征（如圖4所示），位于（6 μm/s，1059.3 μm/s）區間。
3、提出了磁光穿透圖像與內部損傷關聯的研究途徑。亦即，將實驗樣品采用逐層化學腐蝕處理，由磁光圖像定位出穿透位置處的SEM掃描結果進行對比（見圖5），由此得到了帶材內部損傷與裂紋擴展的特征規律。
基于這一新的實驗研究途徑及其實驗結果，發現了YBCO二代超導材料在使役環境下材料損傷的磁場敏感性，即穿透深度遠大于貫穿裂紋，進而揭示出這種損傷模式完全不同于傳統推測的貫穿裂紋破壞模式。其次，通過逐層腐蝕和SEM掃描，還發現了這種YBCO二代帶材從基底向上傳播的新破壞模式，且損傷尖端存在一種非晶化現象，如圖6所示。
圖5. 損傷與磁光圖像的對應。
圖6.
展開 平紋復合材料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為最終失效應變。
展開 濕熱環境下碳纖維復合材料宏-細觀損傷演化Vumat子程序。
濕熱環境下碳纖維復合材料宏-細觀損傷演化Vumat子程序。 感興趣的話和我私聊吧！
ABAQUS后處理之提取分層損傷面積/分層面積/基體損傷面積（ABAQUS+Photoshop） ￥28
ABAQUS后處理之提取損傷面積（ABAQUS+Photoshop聯合使用）
為了定量描述損傷程度，提取載荷造成的損傷面積變得尤為重要，下面介紹損傷面積的提取方法。
1. 去除單元網格，以及邊緣
2. 突出顯示損傷區域，建立損傷與未損區域色差
3. 導出圖片
ABAQUS混凝土損傷塑性模型損傷因子對本構關系影響 附c40~c45混凝土損傷因子ABAQUS輸入
但是ABAQUS塑性損傷模型除了能模擬單調加載的混凝土行為外，更重要的功能就是模擬循環、動態荷載下的混凝土反應，在結構的抗震性能分析能起到很好的作用。
在動荷載作用下，混凝土在受力過程中拉伸和壓縮都會產生損傷造成的裂縫開展，從而導致材料剛度退化。CDP 模型就假定混凝土材料主要因為拉伸開裂和壓縮破碎而破壞，拉伸和壓縮采用不同的損傷因子來描述這種剛度退化，詳見圖 1、圖 2。
圖中E0是材料初始未受損的彈性剛度。損傷變量dc和dt分別為壓縮和拉伸條件下的損傷因子，表示彈性剛度的退化。損傷后的彈性模量為（1-dc）E0，或（1-dt）E0。損傷因子dc或dt=0時表示沒有損傷，dc或dt=1時表示材料失去強度。
那么混凝土的塑性損傷本構模型中的損傷因子到底對混凝土的應力-應變曲線有什么影響呢？讓我們采用100mm*100mm*300mm的混凝土棱柱體模型來做個測試看一下。
依然采用C110級混凝土的本構關系，混凝土的屈服應力和非彈性應變表格如下。子選項中損傷參數和非彈性應變關系的表格也在圖中給出。
但是注意上圖中紅色框部分默認是不填的，即下圖中的混凝土壓縮損傷——拉伸恢復因子wt，混凝土拉伸損傷——壓縮復原因子wc，默認是不填的。
因為CDP模型假定混凝土從拉伸到壓縮時裂縫會閉合，剛度會恢復；從壓縮到拉伸時裂縫仍然存在，剛度不會恢復。因此在ABAQUS中不填的話默認wt(拉伸剛度恢復因子)=0，wc(壓縮剛性恢復因子)=1.
下圖為損傷因子和剛度恢復因子在混凝土載荷循環中對混凝土本構模型的影響。
展開 
abaqus調用damask實現FCC，BCC，HCP多晶織構演化和應力應變場分布模擬
FCC------以鋁為代表，參數使用原始abaqus提供的參數
織構演化模擬模型使用包含1000個單元的1*1*1mm立方體，其中每個單元表示一個特定取向的單晶，初始織構使用軟件生成1000組隨機取向，并分配給不同的單元，模型和初始織構如下圖所示，
利用周期性邊界條件分別模擬多晶沿著ND方向拉伸，壓縮，以及沿著ND方向進行平面應變壓縮時的織構
RD拉伸織構：
RD壓縮織構：
ND平面應變壓縮織構：
BCC------以鐵素體為代表，參數使用原始abaqus提供的參數
織構演化模擬模型使用包含1000個單元的1*1*1mm立方體，其中每個單元表示一個特定取向的單晶，初始織構使用軟件生成1000組隨機取向，并分配給不同的單元，模型和初始織構如下圖所示，
、
利用周期性邊界條件分別模擬多晶沿著ND方向拉伸，壓縮，以及沿著ND方向進行平面應變壓縮時的織構
拉伸織構：
壓縮織構：
平面應變壓縮織構：
HCP------以鎂為代表，參數使用原始abaqus提供的參數
織構演化模擬模型使用包含1000個單元的1*1*1mm立方體，其中每個單元表示一個特定取向的單晶，初始織構使用軟件生成1000組隨機取向，并分配給不同的單元，模型和初始織構如下圖所示，
利用周期性邊界條件分別模擬多晶沿著ND方向拉伸，壓縮，以及沿著ND方向進行平面應變壓縮時的織構
拉伸織構：
壓縮織構：
平面應變壓縮織構：
多晶局部應力應變場分布模擬與宏觀應力應變響應。以FCC－鋁為例子。BCC與HCP同理。
展開 Abaqus調用damask實現軋制變形中FCC，BCC織構演化分析------案例六
Abaqus調用damask實現軋制變形中FCC，BCC織構演化分析
案例實操一
1，使用abaqus建立20*20*20（mm）的立方塊
2，對立方塊進行單元劃分共包含1000個單元
3，假設每個單元代表一個單獨的晶粒，通過腳本隨機賦予每個單元材料屬性
4，施加對應的邊界提交（60%的下壓量）
5，提交與后處理材料數據
包含1000個晶粒的有限元模型
材料的初始取向分布
FCC軋制后的取向分布情況
BCC軋制后的取向分布情況
abaqus調用damask實現FCC，BCC，HCP多晶織構演化和應力應變場分布模擬
FCC------以鋁為代表，參數使用原始abaqus提供的參數
織構演化模擬模型使用包含1000個單元的1*1*1mm立方體，其中每個單元表示一個特定取向的單晶，初始織構使用軟件生成1000組隨機取向，并分配給不同的單元，模型和初始織構如下圖所示，
利用周期性邊界條件分別模擬多晶沿著ND方向拉伸，壓縮，以及沿著ND方向進行平面應變壓縮時的織構
RD拉伸織構：
RD壓縮織構：
ND平面應變壓縮織構：
BCC------以鐵素體為代表，參數使用原始abaqus提供的參數
織構演化模擬模型使用包含1000個單元的1*1*1mm立方體，其中每個單元表示一個特定取向的單晶，初始織構使用軟件生成1000組隨機取向，并分配給不同的單元，模型和初始織構如下圖所示，
、
利用周期性邊界條件分別模擬多晶沿著ND方向拉伸，壓縮，以及沿著ND方向進行平面應變壓縮時的織構
拉伸織構：
壓縮織構：
平面應變壓縮織構：
HCP------以鎂為代表，參數使用原始abaqus提供的參數
織構演化模擬模型使用包含1000個單元的1*1*1mm立方體，其中每個單元表示一個特定取向的單晶，初始織構使用軟件生成1000組隨機取向，并分配給不同的單元，模型和初始織構如下圖所示，
利用周期性邊界條件分別模擬多晶沿著ND方向拉伸，壓縮，以及沿著ND方向進行平面應變壓縮時的織構
拉伸織構：
壓縮織構：
平面應變壓縮織構：
多晶局部應力應變場分布模擬與宏觀應力應變響應。以FCC－鋁為例子。BCC與HCP同理。
展開 ABAQUS UMAT - 混凝土塑性損傷模型的實現 ￥1500
混凝土塑性損傷模型在工程上應用較為廣泛，同類型的本構模型多內置于各類仿真軟件中，供用戶模擬混凝土結構的破壞和受力情況。本文根據Peter Grassl 和 Milan Jira′sek 2006年的文章《Damage-plastic model for concrete failure》進行本構模型代碼復現，并對文中的模型進行了一些簡化。
UMAT代碼和INPUT文件見付費內容
