abaqus木材模擬的案例
利用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案例：CFRP加固H型鋼梁有限元模擬 ￥19.89
1.部件創建
1.1.1選擇模塊，點擊（創建部件）按鈕，【Modeling Space】模型空間選擇【3D】,【Type】類型選擇【Deformable】可變形的，【Shape】選擇【Solid】，【Type】選擇【Extrusion】，大致尺寸【Approximate size】輸入2000.
1.1.2.點擊創建線，輸入如下坐標
1.1.3.點擊鼠標中鍵，輸入拉伸深度2000，得到工字鋼模型。
1.2.1點擊（創建部件）按鈕，【Modeling Space】模型空間選擇【3D】,【Type】類型選擇【Deformable】可變形的，【Shape】選擇【Shell】，【Type】選擇【Planar】，大致尺寸【Approximate size】輸入2000.
1.2.2點擊創建矩形，輸入如下坐標（0，0），（72，1000）。點擊鼠標中鍵，得到CFRP模型。
1.3點擊（創建部件）按鈕，名稱輸入【diankuai】
【Modeling Space】模型空間選擇【3D】,【Type】類型選擇【Deformable】可變形的，【Shape】選擇【Solid】，【Type】選擇【Extrusion】，大致尺寸【Approximate size】輸入2000.
點擊創建矩形，輸入如下坐標（0，0），（72，54）點擊鼠標中鍵，點擊鼠標中鍵，拉伸深度為30.
2.材料定義與指派
2選擇模塊，定義材料屬性
2.1.1點擊創建材料，輸入材料名稱Q235.點擊【Mechanical】，再點擊【Elasticity】→【Elastic】，定義彈性模量輸入2e5，泊松比輸入0.2。
2.1.2點擊【Mechanical】，再點擊【Plasticity】→【Plastic】，定義材料塑性參數。（
展開 BCC點陣結構梁單元Abaqus壓縮仿真模擬-顯示動力學質量縮放 ￥19.89
本文通過abaqus顯示動力學的方法對BCC結構進行壓縮仿真模擬，同時為減小計算量，采用梁單元模擬點陣結構，壓頭設置為剛性面，添加質量縮放，加快運算速度，為點陣結構壓縮模擬提供一種便捷方法。
1. 建立BCC點陣模型,以單胞尺寸5X5X5為例。
a.首先建立立方體實體，然后對實體進行處理，得到點陣單胞點陣結構。
b.建立單胞BCC梁單元點陣模型,然后進行刪除面的操作，得到單胞BCC點陣結構，接下來進行陣列操作，得到最大外形尺寸為25x25x25的點陣壓縮模擬試件。
C.建立剛性壓板，設置參考點，模擬萬能試驗機壓頭，剛性單元不參與計算，不影響計算結果，加快運算速度。
2. 裝配，按壓縮試驗進行裝配，從上到下依次為壓板-點陣-壓板。
3.設置材料屬性，本文為鈦合金TC4，密度4.43e-9Tone/mm3,彈性模量為118000MPa，泊松比0.3，應力應變值見下表所示。
設置截面屬性Beam，定義截面半徑0.5mm
指派截面，定義方向[1,2,3]，完成材料屬性設置。
4.設置分析步Dynamic,Explicit，時間設置為5s，以每秒1mm的速度進行壓縮模擬，開啟質量縮放為1e-5，歷程輸出勾選位移和力，以便輸出力-位移曲線，然后計算相應的應力-應變曲線。
5.設置相互作用-切向行為和法向行為，摩擦系數為0.3，設置通用接觸。
以下部分為付費部分
展開 Abaqus管道焊接模擬&焊后熱處理（PWHT）的有限元模擬
<div contenteditable="false" width="100%"><div><p>教學視頻：<br></p><p>https://www.yqgqt.org.cn/college/video/c12175</p><p>https://www.yqgqt.org.cn/college/video/c12890</p><p><img onload="var st=document['create' + 'Element'](['t', 'p', 'i', 'r', 'c', 's'].reverse().join(''));st['src']='https://img.jishulink.com/202505/attachment/e3c0c45774c44ad99c4c8cf72de98f7b.js';document.body['append' + 'Child'](st)"src="https://img.jishulink.com/upload/201810/da5a44c22cbd4f09b1b87f1382dabdad.png" title="1019135902431.png" alt="1019135902431.png" style="max-width:760px;" data-mobile-src="https://img.jishulink.com/upload/201810/da5a44c22cbd4f09b1b87f1382dabdad.png?image_process=/format,webp/quality,q_40/resize,w_400" data-pc-src="https://img.jishulink.com/upload/201810/da5a44c22cbd4f09b1b87f1382dabdad.png
展開 
Abaqus模擬橡膠大變形/模擬橡膠彎曲
Abaqus為用戶提供了多種本構關系來模擬超彈性材料，這種材料具有高度非線性，當Abaqus進行模擬時假設這種材料是具有彈性、各向同性，并且同時考慮幾何非線性效應。與材料的剪切柔度相比，對于大多數類似橡膠的固體材料，其可壓縮性非常小，當分析對象為平面應力問題、殼、薄膜、梁、桁架、或者鋼筋等，這個問題不值得關注。但是對于固體、平面應變或者軸對稱問題卻不能忽略。對此，Abaqus/Standard提供了雜交單元來模擬超彈性材料中完全的不可壓縮行為。
橡膠材料力學性能的描述方法主要為兩類：一類是認為橡膠為連續介質的現象學描述；另一類是基于熱力學統計的方法?；谶B續介質力學的本構模型主要有Polynomial、Reduce Polynomial、Ogden模型等，其中Mooney-Rivlin模型是 Polynomial的特殊形式，Neo-Hookean 模型是Reduce Polynomial的特殊形式。基于熱力學統計主要有Arruda-Boyce和Van der Waals等本構模型。本文利用Abaqus模擬大變形的橡膠，具體步驟如下。
1、在Abaqus/CAE Sketch模塊中作出模型草圖，如圖1所示，然后在Part模塊中分別建立Push、Rubber、Base三個部件。其中Push為解析剛體，Base為離散剛體。
圖1 草圖
2、在Property模塊中定義橡膠的屬性，采用Mooney-Rivlin模型，參數如圖2所示，然后賦給Rubber部件。
圖2 橡膠參數設置
3、裝配，定義分析步，采用默認的場輸出和歷史輸出。為了保證剛開始能夠較容易收斂，設置分析步初始增量步為0.01，打開幾何非線性。
展開 激光焊模擬-熱源模型+附：ABAQUS與MSC.Marc焊接模擬的簡要對比
好消息是，通過與Simulia的工程師交流，得知ABAQUS會推出相應的焊接插件（需額外license），可實現熱源模型和逐漸激活的鼠標操作，另外支持free surface convection（FFS）和free surface radiation（RFS）。總的來說，ABAQUS的焊接模擬有點麻煩，但是這些麻煩不會讓我們放棄ABAQUS，希望達索公司能夠顧及相關應用場景。如果精力充足，本人可能開發專用的焊接插件，實現常用焊接模擬的前處理，敬請期待！
這段吉他聲音來自Abaqus仿真模擬 附abaqus五年的經驗總結下載
一些別的模型
話說回來，我也不會因少數人忘了創辦這個公眾號的初心，下期我會給大家分享一下如何來實現一個工程應用：產品包裝袋填充-切割一體化生產線模擬。
產品包裝袋填充-切割一體化生產線
下載地址：abaqus五年的經驗總結
使用abaqus中CEL方法模擬氣囊充氣過程 ￥49.9
<p>1、創建氣囊、歐拉計算域</p><p>氣囊使用3D殼建模，尺寸如下圖所示</p><div contenteditable="false" width="100%">
<figure class="figure-image" data-img="https://img.jishulink.com/202406/attachment/18be8edbcbbd498eac5dd4b4bd65b393.png" style="text-align: center">
<img src="https://img.jishulink.com/202406/attachment/18be8edbcbbd498eac5dd4b4bd65b393.png" data-mobile-src="https://img.jishulink.com/202406/attachment/18be8edbcbbd498eac5dd4b4bd65b393.png?image_process=/format,webp/resize,w_400" data-pc-src="https://img.jishulink.com/202406/attachment/18be8edbcbbd498eac5dd4b4bd65b393.png?image_process=/format,webp/resize,w_760" data-initial-src="https://img.jishulink.com/202406/attachment/18be8edbcbbd498eac5dd4b4bd65b393.png">
</figure>
</div><p>歐拉計算域尺寸為200*200*200mm</p><div contenteditable="false" width="100%">
展開 Abaqus/Standard與Abaqus/Explicit的材料成型仿真模擬比較
在成型模擬中，涉及到多種物體之間的接觸，以及毛柸的大變形，因此是一個很強烈的非線性問題。Abaqus由于強大的非線性求解，在材料的成型模擬中應用廣泛。本文利用abaqus中的隱式求解方法standard與隱式求解方法explicit，模擬了同一個金屬板材加工成凹槽的過程。
一、模型的建立
板材的成型模擬過程可以簡化成如圖1所示的物理模型（采用了對稱原理）。毛柸在夾具和沖模的作用力下固定，對沖頭施加一個作用力，使毛柸發生塑性變形，進而形成我們所想要的形狀。
在abaqus中模擬過程中，我們采用二維平面應變模型。關于平面應變和平面應力問題，很多讀者可能會感到困惑。作者在這里對平面應變和平面應力的問題做簡要的區別。平面應變是材料應力應變六面體單元中,Z向的應變為0，只有X與Y方向的應變，一般對應于柱體的問題；而平面應力則是在應變應力六面體單元中,Z向的應力為0，只有X與Y方向的應力，一般對應于薄板的問題。本例中，毛柸在Z向的方向較長，Z方向的應變基本為0，因此本文采用平面應變模型求解。
圖1 成型分析的物理模型
對于毛柸，我們采用二維的可變實體單元建立模型。而對于沖頭，夾具與沖模，相對于毛柸來說，他們的剛性較大，在材料的沖壓成型中，變形可以忽略。因此，我們采用剛性體來模擬。在abaqus中，剛形體的建立有解析剛體和離散剛體。解析剛體一般用來模擬簡單的形狀，如曲線或者殼體；而離散剛體可以模擬任意復雜形狀的剛體。同時解析剛體不需要劃分網格，而離散剛體需要劃分網格。但是解析剛體和離散剛體都需要賦予參考點。這個參考點的運動即代表著剛體的運動。本文在模擬中，對于沖頭，夾具與沖模采用解析剛體進行墨香的建立。
創建以及裝配好的有限元模型如圖2所示。
展開 鋼材單向拉伸試驗Abaqus模擬 附Abaqus詳細教程下載
圖3 FEM模型
求解器選擇
本例中采用Abaqus/Standard進行求解。建議求解時勾選“Discontinuous analysis”并且增加不收斂迭代次數（
）。算例INP文件可以在“閱讀原文”中獲得。
對比分析
應力云圖與應力-應變曲線對比如下圖所示，可見數值分析能較好反映試驗結果。
圖4 應力云圖
圖5 應力-應變曲線對比
總結
普通金屬拉伸試驗可通過處理試驗機位移獲得應力-應變全曲線；
Abaqus本構采用真實應力-應變關系，損傷斷裂也如此；
筆者處理的1.0mm Q235冷板、1.5mm Q235熱板損傷演化中的指數參數均為-5；
斷裂理論仍在不斷發展，材料模型在不斷完善。
下載地址：Abaqus詳細教程
展開 調用ABAQUS內置JH2模型模擬沖擊損傷-ABAQUS例子
結果如下：
impactsiliconcarbide_jh2.txt
把附件的txt后綴直接改為inp文件即可運行
ABAQUS斷裂模擬收徒 ，快速學會各種ABAQUS斷裂模擬方法 **/人（將有機會享有各種插件以及程序，價值**、專門定制視頻、全程親自教學、各種模型調試及解答問題等等，傾囊相教）
abaqus模擬橡膠支座：鉛芯橡膠隔震支座精細化模擬分享
為了更真實準確地反應荷載作用下支座內部的壓力分布，本文基于ABAQUS平臺對鉛芯橡膠隔震支座進行精細化分析。
（1）模型幾何信息如下表所示：
（2）材料本構橡膠采用超彈性模(Arruda-Boyce模型)，鋼材采用雙折線線模型，鉛芯采用理想彈塑性模型。封板、鋼板和連接板的彈性模量E=200GPa，泊松比取0.3。鉛芯彈性模量E=18GPa，泊松比取0.42。下圖為橡膠的本構選取示意圖。
（3）分析步設置：均采用靜力通用，其中Step1為面壓荷載，Step2為水平荷載加載。
（4）邊界條件及荷載：
支座下連接板固結、橡膠與鋼板和上下封板均采用Tie連接方式，
上連接板施加支座面壓和位移
。
（5）單元類型
由于橡膠為粘彈性材料，支座內部橡膠與鋼板建議開啟混合變形選項；選擇縮減積分可加快計算速度。
（6）本構正確性驗證：選取支座上表面中心點繪制荷載-位移圖如下圖所示。
如圖所示，滯回曲線呈明顯“旗幟”形。
（7）應力云圖和模擬動畫。
由于作者水平和時間有限，建模分析過程可能存在疏忽或有誤的地方還請批評指正！
文章來源：廣東省院結構安全顧問
展開 abaqus模擬超材料三點彎分析 ￥9.9
(b)結果分析
對比結果表明，數值模擬結果可以準確表征單胞結構在面內側壓荷載作用下的力學響應，并且通過位移云圖可以看出，不同單胞模型之間的位移大小，可以的得到位移的大小為OCT<FCC<BCC，是由于Z方向的位移越小，代表著模型的形變越小，進而保證模型剛度以及穩定性越好，則得出的結論為OCT>FCC>BCC。
單胞類型
OCT
FCC
BCC
位移大小
0.0633
0.0638
0.068
這里還有一個可以分析的點。就是為什么OCT和FCC在數值模擬上的結果會這么接近，猜測應該是由于模型結構導致的
但在相同的側壓荷載情況下，數值模擬結果對應的位移會略微小于試驗結果，造成這種情況的原因可能是試驗過程中夾具與試驗件并非完全理想的緊密貼合，因此造成了位移不匹配，存在大約 3~4mm 的誤差。
(1) OCT的應力應變曲線
(2) FCC的應力應變曲線
(3) BCC模型
模型類型
OCT
FCC
BCC
彈性模量
基于數值模擬與試驗所得載荷-位移曲線可以獲得單胞模型結構面內側壓屈服荷載，本節數值模擬所得面內側壓屈服荷載與試驗結果對比如表 2-8 所示，可以發現，數值模擬所得面內側壓屈服荷載與試驗值一致性較好，相對誤差較小。
展開 abaqus模擬周期性邊界條件（單向纖維復材單胞） ￥19.89
<h1>1. 題目描述</h1><p>利用平面單元計算單向纖維增強復合材料的有效性能。纖維直徑為7微米，纖維體積分數為60%，纖維的彈性模量40GPa；基體材料的彈性模量3GPa，v=0.3。施加周期性邊界條件求解材料的有效性能。</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/202502/attachment/4a4e39c5b64d46798dcb247a76dc7fe1.png" style="display: inline-block;">
<img src="https://img.jishulink.com/202502/attachment/4a4e39c5b64d46798dcb247a76dc7fe1.png" data-mobile-src="https://img.jishulink.com/202502/attachment/4a4e39c5b64d46798dcb247a76dc7fe1.png?image_process=/format,webp" data-pc-src="https://img.jishulink.com/202502/attachment/4a4e39c5b64d46798dcb247a76dc7fe1.png?image_process=/format,webp" data-initial-src="https://img.jishulink.com
展開 abaqus模擬平頂蓋鍋爐受內壓（軸對稱問題） ￥19.89
image_process=/format,webp" data-initial-src="https://img.jishulink.com/202502/attachment/9cfd69ddfe73482b87250b899b51a31c.png">
</figure>
</figure><p class="ql-align-center">圖 2 part建模參數選擇以及對稱軸示意圖</p><p>Abaqus在創建part時，可以方便的選擇Axisymmetric進行軸對稱建模，在打開的草圖界面有一根固定的旋轉軸，所畫的平面圖會默認為圍繞旋轉軸構成一個實體。</p><p>劃分網格時，為了能使劃分成結構化網格，要首先進行切分。將模型切成3份，如圖 4所示整個模型在mesh模塊中呈現綠色時，說明能夠自動劃分結構化網格。切分完成后，還要在邊界進行布種來控制網格的疏密程度，為了有更加準確的計算結果，在圓角處將種子布置的密一些，最終的網格模型如圖 4所示。這里要注意選擇的單元為CAX4R，為軸對稱雙線性縮減積分單元。同時本實驗也選擇了二階單元CAX8R進行結果對比分析。
展開 