ABAQUS蠕變裂紋的案例
abaqus蠕變基本設置及雙曲正弦函數損傷蠕變本構CREEP子程序 ￥59.9
該部分為abaqus蠕變計算基本流程
ABAQUS蠕變問題計算流程.pdf
付費部分為使用CREEP子程序建立雙曲正弦函數蠕變損傷子程序，含到達預設損傷值（假設為1.0）后終止計算，和USDFLD子程序控制材料參數（該子程序可用于損傷后的材料退化，如蠕變第三階段或者蠕變疲勞分析，若不需要場變量控制可對該部分代碼進行刪除），相關理論請參考附件sci文獻。可提供關于CREEP子程序的幫助文件學習的相關指導
abaqus蠕變分析例子 ￥2
對粘彈性材料試件進行拉伸蠕變分析，得到試件的蠕變曲線
包括材料屬性，分析步設置及最后得到的蠕變曲線。
ABAQUS土體蠕變模擬
最近做了粉黏土的蠕變實驗，接下來要進行數模的模擬與實驗數據對應，做一個簡單的三軸蠕變試驗模擬，但是材料模型選取上有沒有大佬提出一些建議，是用自定義粘彈性子程序umat，還是DP creep 再或者是自帶蠕變經驗模型，希望有大佬可以探討一下，umat代碼有償
ABAQUS直裂紋、斜裂紋圍道積分計算裂紋尖端J積分
之前算過一個關于裂紋擴展的問題，當時創建裂紋選擇的是contour intergral，后來又有人咨詢我裂紋尖端J積分的計算問題。我才恍然大悟，其實圍道積分方法還是適用于計算裂紋尖端在某時刻的J積分，至于動態擴展問題，還是交給XFEM吧（雖然也不太好）。
計算了幾種情況下的裂紋尖端J積分，包括直裂紋、斜裂紋以及裂紋尖端傾斜等三種情況。
部分試件的應力分布及J積分結果如圖所示：
ABAQUS蠕變UMAT本構
creep.for
免費分享不帶損傷的UMAT蠕變本構，需要Sinh損傷本構的請移步高溫合金蠕變損傷本構/UMAT子程序/Sinh 蠕變損傷本構/論文復現_本構模型 損傷模型-技術鄰
Abaqus蠕變分析（step by step）
事實上蠕變是非常復雜的，這里僅給出了abaqus中的簡單流程，足以解決常規工程問題。
為了簡化塑料結構蠕變問題的計算(如降低蠕變應變與其他非彈性應變的耦合程度)，可以將該分析問題分成一個靜態加載的過程，然后再進行蠕變過程的分析。
1.靜態加載過程的計算
靜態加載過程就是一與時間無關的加載過程，使用ABAQUS/Standard時主要是在中設置，如圖1所示。
2.蠕變過程的計算
在通過步驟1的靜態分析后，結構中將產生一個應力場，接下來可以進行蠕變過程的計算。蠕變過程的計算主要分為兩個過程：獲得該結構材料的蠕變模型參數和建立蠕變分析步。
1) 獲得材料的蠕變模型參數
目前ABAQUS蠕變模型有三種，分別是Power-law model和Hyperbolic-sine law model。其中Power-law model有兩種形式為Time hardening form和Strain hardening form。其中Time hardening form形式最為簡單，對于簡單的蠕變過程（如蠕變過程應力變化范圍不太大）是比較適用的，式(1)為其微分形式：
由于圖2中表征的是蠕變應變與時間和等效應力的關系，故必須對公式(1)積分，積分結果見公式(2)：
表征材料蠕變特性的三個參數確定后，通過ABAQUS/CAE的添加材料的蠕變特性，如圖3所示：
2) 蠕變計算
由于蠕變是一個時間相關的過程，因此必須計入時間。同時蠕變又是一個慣性效應不明顯的過程，即結構的加速度效應不用考慮。針對這些ABAQUS提供了專門針對這一類型的分析步。
蠕變計算分析步設置在中完成，見圖4。
其中蠕變應變的容差設置將影響增量步的大小，容差設的很小，增量步也將降低。
展開 Abaqus蠕變K-R子程序代碼 ￥20
Abaqus蠕變K-R子程序-forture源代碼及學習文檔
ABAQUS——CREEP蠕變子程序 ￥10
首先給出蠕變定義時，材料對應的輸入參數為：
引自：《Creep and Swelling》
其中，前兩種的公式內涵如下
時效硬化：比較好編寫，求導前的函數較為簡單
子程序為： time hardening下面片段
應變硬化公式為：
對上式進行整理可以得到：
程序中DECRA(1)為下式：
子程序為： C strain hardening 下面片段
比較上述四個模型：相同蠕變模式下的結果相同。
子程序+案例模型（cae、inp）
abaqus雙曲正弦金屬蠕變子程序 ￥20
金屬蠕變子程序，內含相應的碩士學位論文，共計7個參數，使用的是雙曲正弦蠕變本構方程，可以用作子程序學習以及金屬蠕變仿真的參考。
ABAQUS裂紋尖端應變、裂紋擴展模擬及問題
前幾天有人問我ABAQUS做焊點分析，我一看他給我的一片文獻，其實是用ABAQUS做裂紋擴展分析。之前也沒接觸過裂紋分析，于是照貓畫虎做了個算例，但是裂紋沒有擴展。
ABAQUS做裂紋有三種方法：contour integral，擴展有限元及VCCT法，這里用了contour integral法。
如圖所示，V形楔形處有一個預制裂紋，是采用Interaction模塊的assign seam設定的，裂紋的擴展面及方向是通過crack來設定的，類型為contour integral。材料模型定義了塑性應力-應變關系，彈性參數、GTN參數、脆性失效參數等。模型上的兩個孔，一個固支、一個勻速拉。預期當裂紋尖端的單元變形達到某一個值時將刪除單元。
您看見了就給個意見唄。
步驟：
建立模型，進行適當的partition
定義材料：分別定義了elastic彈性參數、plastic真實應力-應變關系、GTN模型參數、脆性失效參數（包括一個叫演化參數）。
定義預制裂紋、定義裂紋擴展面、方向，定義失效單元的generation。
邊界條件，提交job，查看結果。
結果：預期模型在塑性變形不是很大時就會產生裂紋擴展，但是模型產生了很大塑性變形后仍然沒有發生失效。
Mises應力場：
x方向正應力場
x方向真實應變場
x方向塑性應變場
裂紋尖端應變的結果還是挺漂亮的，雖然正確性有待考證，如果裂紋出來了就完美了，可惜裂紋沒出來。
展開 Abaqus蠕變本構定義插件——Anand、Darveaux、Double Power ￥200
Anand模型、Darveaux和雙冪律模型是電子封裝領域常用于描述焊料蠕變行為的本構模型，Abaqus內置了這三類本構模型，但大多版本的CAE界面中卻沒有為其提供定義的方式，需要自行更改關鍵字。本文介紹了一種可以直接在用戶界面定義三類蠕變本構參數的插件，可避免更改關鍵字過程的輸入錯誤。
Anand模型
Darveaux模型
雙冪律（Double Power）模型
插件需要購買激活，歡迎咨詢！
QQ: 654870152
利用Abaqus的UMAT子程序仿真木材蠕變現象
蠕變現象
蠕變現象簡單講，就是載荷不變，材料或者結構變形隨著時間的推移而逐漸增加的現象。引起蠕變的原因有很多，溫度、材料本構、載荷水平等等。從微觀機理上可以追溯到晶體結構。當然這不是我們做工程的該考慮的。
我們要考慮的是，如何用一個本構來描述這種變形特點。下面這個圖具備相當的代表性。
通常這類問題一個顯著的特點就是和時間相關，類似于一個生長現象。搞自然科學的，看到這基本都能猜到，這個本構一定要有自然常數。下面我們將以木材蠕變為例，介紹下在ABAQUS UMAT中如何實現蠕變仿真。
2. 本構理論
文獻[1]給出了木材蠕變過程中本構：
3. 算例
3.1 模型
考慮懸臂梁模型，如下圖。
3.2 邊界條件
根據蠕變的定義，模型必須現有一個穩定的載荷，因此可以分成兩個step。第一個step，完成力加載，第二個step保持載荷，實現蠕變變形的生長。
需要指出的是，蠕變通常需要在較長的時間尺度上才能有明顯的效果。比如我們要觀察100天的變形情況，那么這個時候，ABAQUS設定的總時間還是1，在UMAT里面要乘以相應的系數，給出物理時間，才能有效的實現蠕變效果。
3.3 結果
最終得到懸臂梁端部位移如下圖所示。從圖中可以看出，在最初的幾天，蠕變變形較大，隨著時間的推移，變形增加的幅度放緩，符合蠕變的特點。
參考文獻
[1] 《濕度變化和荷載作用下膠合木曲梁的工作性能研究》
[2] 《旋切板膠合木的蠕變及其對結構穩定性的影響》
最后，有abaqus二次開發相關需求歡迎大家通過我們的微信公眾號聯系我們。
微信公眾號：320科技工作室
展開 Abaqus粘彈性材料蠕變測試仿真案例講解
Abaqus粘彈性材料蠕變測試仿真案例講解
abaqus考慮混凝土蠕變流變的三點彎曲梁跨中撓度模擬
abaqus考慮混凝土蠕變流變的三點彎曲梁跨中撓度模擬。 模型考慮了混凝土的蠕變效應，蠕變規律依據文獻取值，并與文獻進行了對比。
1.模型 2.蠕變子程序
土木工程博士畢業，具有abaqus的10使用經驗，精通各種模型及二次開發，可以幫助解決各種模型問
基于ABAQUS蠕變儲層稠油蒸汽吞吐開發過程數值模擬
因為這些層位遠離頻繁加載，所以可不考慮蠕變。
式（1）
表1 地表土層與泥巖層Drucker-Prager模型參數
蠕變特性顯著的中間層
D1到D7使用改進的Drucker-Prager Cap塑性模型模擬，材料屬性數據列于表2中。因其考慮蠕變，因此須k=1.0，且剪切破壞面與蓋帽間不可有過渡區域（即0）；硬化/軟化行為通過屈服應力與體積應變的關系曲線定義，數據列于表2中；初始蓋帽屈服面的位置設定為0.02，如果應力位于蓋表面之外，Abaqus會自動調整蓋屈服面的位置。
表2 D1到D7層Drucker-Prager Cap模型參數
使用式（2）所示Singh-Mitchell蠕變模型對固結蠕變進行模擬。 蠕變模型參數表3所示，其為溫度的函數，在cap creep model中指定，如圖3所示。
式（2）
式中:
表3 蠕變模型參數
圖3 蠕變模型參數定義
儲層滲透系數定義為溫度的函數，另外為了實現熱力耦合，還需要定義熱膨脹系數，本例熱膨脹系數定義為：5.76E–6 1/°C (3.2E–6 1/°F)，對于熱膨脹系數的含義，以線性膨脹為例如圖4所示。
展開 