金屬蠕變子程序的案例
abaqus雙曲正弦金屬蠕變子程序 ￥20
金屬蠕變子程序，內(nèi)含相應(yīng)的碩士學(xué)位論文，共計7個參數(shù)，使用的是雙曲正弦蠕變本構(gòu)方程，可以用作子程序學(xué)習(xí)以及金屬蠕變仿真的參考。
abaqus蠕變基本設(shè)置及雙曲正弦函數(shù)損傷蠕變本構(gòu)CREEP子程序 ￥59.9
該部分為abaqus蠕變計算基本流程
ABAQUS蠕變問題計算流程.pdf
付費部分為使用CREEP子程序建立雙曲正弦函數(shù)蠕變損傷子程序，含到達(dá)預(yù)設(shè)損傷值（假設(shè)為1.0）后終止計算，和USDFLD子程序控制材料參數(shù)（該子程序可用于損傷后的材料退化，如蠕變第三階段或者蠕變疲勞分析，若不需要場變量控制可對該部分代碼進(jìn)行刪除），相關(guān)理論請參考附件sci文獻(xiàn)。可提供關(guān)于CREEP子程序的幫助文件學(xué)習(xí)的相關(guān)指導(dǎo)
ABAQUS——CREEP蠕變子程序 ￥10
首先給出蠕變定義時，材料對應(yīng)的輸入?yún)?shù)為：
引自：《Creep and Swelling》
其中，前兩種的公式內(nèi)涵如下
時效硬化：比較好編寫，求導(dǎo)前的函數(shù)較為簡單
子程序為： time hardening下面片段
應(yīng)變硬化公式為：
對上式進(jìn)行整理可以得到：
程序中DECRA(1)為下式：
子程序為： C strain hardening 下面片段
比較上述四個模型：相同蠕變模式下的結(jié)果相同。
子程序+案例模型（cae、inp）
Abaqus蠕變K-R子程序代碼 ￥20
Abaqus蠕變K-R子程序-forture源代碼及學(xué)習(xí)文檔
蠕變損失子程序（k-r,l-m,sinh）
提供蠕變損傷子程序
利用Abaqus的UMAT子程序仿真木材蠕變現(xiàn)象
蠕變現(xiàn)象
蠕變現(xiàn)象簡單講，就是載荷不變，材料或者結(jié)構(gòu)變形隨著時間的推移而逐漸增加的現(xiàn)象。引起蠕變的原因有很多，溫度、材料本構(gòu)、載荷水平等等。從微觀機理上可以追溯到晶體結(jié)構(gòu)。當(dāng)然這不是我們做工程的該考慮的。
我們要考慮的是，如何用一個本構(gòu)來描述這種變形特點。下面這個圖具備相當(dāng)?shù)拇硇浴?通常這類問題一個顯著的特點就是和時間相關(guān)，類似于一個生長現(xiàn)象。搞自然科學(xué)的，看到這基本都能猜到，這個本構(gòu)一定要有自然常數(shù)。下面我們將以木材蠕變為例，介紹下在ABAQUS UMAT中如何實現(xiàn)蠕變仿真。
2. 本構(gòu)理論
文獻(xiàn)[1]給出了木材蠕變過程中本構(gòu)：
3. 算例
3.1 模型
考慮懸臂梁模型，如下圖。
3.2 邊界條件
根據(jù)蠕變的定義，模型必須現(xiàn)有一個穩(wěn)定的載荷，因此可以分成兩個step。第一個step，完成力加載，第二個step保持載荷，實現(xiàn)蠕變變形的生長。
需要指出的是，蠕變通常需要在較長的時間尺度上才能有明顯的效果。比如我們要觀察100天的變形情況，那么這個時候，ABAQUS設(shè)定的總時間還是1，在UMAT里面要乘以相應(yīng)的系數(shù)，給出物理時間，才能有效的實現(xiàn)蠕變效果。
3.3 結(jié)果
最終得到懸臂梁端部位移如下圖所示。從圖中可以看出，在最初的幾天，蠕變變形較大，隨著時間的推移，變形增加的幅度放緩，符合蠕變的特點。
參考文獻(xiàn)
[1] 《濕度變化和荷載作用下膠合木曲梁的工作性能研究》
[2] 《旋切板膠合木的蠕變及其對結(jié)構(gòu)穩(wěn)定性的影響》
最后，有abaqus二次開發(fā)相關(guān)需求歡迎大家通過我們的微信公眾號聯(lián)系我們。
微信公眾號：320科技工作室
展開 ABAQUS/Standard 用戶材料子程序實例 － Johnson-Cook 金屬本構(gòu)模型
ABAQUS_Standard用戶材料子程序實例.pdf
