ANSYS線性本構(gòu)的案例
巖土非線性本構(gòu)UMAT開發(fā) ￥666
巖土非線性本構(gòu)UMAT開發(fā)
等效線性黏彈性本構(gòu)的應用
有個困惑就是，做試驗可以獲得不同圍壓下土體的最大剪切模量，以及剪切模量比與剪應變，阻尼比和剪應變的曲線，然后根據(jù)不同圍壓與最大剪切模量的關(guān)系就知道了公式中的k和n，做模擬的時候，在材料屬性輸入k,n,v,w，關(guān)鍵字中輸入各土體單元的震前圍壓，剪切模量比，阻尼比，最大剪應變。如果假定震前圍壓為100，那在迭代過程用所用到的剪切模量比與剪應變及阻尼比的曲線就是100kpa所對應的曲線。但是如果考慮震前圍壓，就是先做靜力分析求出各單元的有效應力作為關(guān)鍵字輸入中的第一列，那這樣的話基本一層土是一個應力，也就是一層土一種圍壓，一種圍壓對應一個最大剪切模量和關(guān)系曲線，迭代的時候不可能取每層土對應圍壓的下土體的關(guān)系曲線，那么要用哪個圍壓下的關(guān)系曲線？對于正常固結(jié)土，最大的圍壓是密度*g*h，做土力學實驗獲得以上關(guān)系曲線是根據(jù)土體深度來加的圍壓，那如果現(xiàn)在土體密度是2，模擬土體厚度為60m，最大圍壓就是1200kpa，迭代的時候要用圍壓1200kpa對應的關(guān)系曲線嗎？可是考慮靜應力之后，每層土一個圍壓，只有最下層土體圍壓才是1200
展開 Chaboche各向同性非線性隨動硬化行為的材料本構(gòu)模型計算matlab程序 ￥475
Chanboche模型是一種用于描述材料各向同性非線性隨動硬化行為的材料本構(gòu)模型。該模型由Chanboche在1981年提出，其基本形式包括各向同性部分和隨動硬化本構(gòu)部分。
具體而言，Chanboche模型各向同性本構(gòu)部分可以用以下方程表示：
dR(p)=b(Q-R)dp
非線性隨動硬化模型可以用以下方程表示：
dx=(2/3)cdεp-rxdp
本程序已經(jīng)在上一個帖子基礎(chǔ)上進一步完善，實現(xiàn)可直接輸入試驗拉伸循環(huán)曲線，計算本構(gòu)參數(shù)，黑色線為計算結(jié)果，紅色為試驗循環(huán)拉伸應力應變曲線。
展開 flac3d的interface界面單元非線性本構(gòu)模型開發(fā)代碼 ￥12
<p>基于fish語言的flac3d的interface界面單元非線性本構(gòu)模型開發(fā)實例</p>
非線性本構(gòu)關(guān)系在ABAQUS中的實現(xiàn)
非線性本構(gòu)關(guān)系在ABAQUS中的實現(xiàn)-闞前華，康國政，徐祥
歡迎大家引用和閱讀！
ABAQUS umat 非線性等向硬化本構(gòu)模型（Voce 硬化模型） ￥129
<p class="ql-align-justify">本資源包含一份 PDF 文檔和可直接編譯運行的 Fortran UMAT 代碼，具體內(nèi)容為：</p><p class="ql-align-justify">非線性等向硬化本構(gòu)模型（Voce硬化模型） + 隱式積分 + 徑向返回</p><p class="ql-align-justify">完整公式推導 + Fortran 源碼直接編譯</p><p class="ql-align-justify">完整的算法一致切線模量推導與實現(xiàn)</p><p class="ql-align-justify">PDF 包含規(guī)范化的本構(gòu)方程、隱式積分、徑向返回與一致切線模量推導，可供初學者學習。配套 UMAT 代碼可直接在 ABAQUS 編譯運行，采用全隱式積分搭配一致切線模量，收斂速度極快、計算精度極高，適合初學者快速入門。</p><p class="ql-align-justify">下圖展示了部分PDF內(nèi)容，及umat計算結(jié)果與abaqus內(nèi)置模型對比，可以發(fā)現(xiàn)umat收斂速度極快，與abaqus內(nèi)置模型幾乎一致。
展開 基于回映算法的Chanboche各向同性非線性隨動硬化本構(gòu)matlab程序 ￥369
<p>Chanboche模型是一種用于描述材料各向同性非線性隨動硬化行為的材料本構(gòu)模型。該模型由Chanboche在1981年提出，其基本形式包括各向同性部分和隨動硬化本構(gòu)部分。</p><p>具體而言，Chanboche模型各向同性本構(gòu)部分可以用以下方程表示：</p><p>dR(p)=b(Q-R)dp</p><p>非線性隨動硬化模型可以用以下方程表示：</p><p>dx=(2/3)cdεp-rxdp</p><p>程序基于3個背應力分量編寫，效果參見鏈接<a href="https://www.bilibili.com/video/BV1B54y1F7gS/?vd_source=9f1dda2358e63ace0b661e56fe417806" rel="noopener noreferrer" target="_blank">https://www.bilibili.com/video/BV1B54y1F7gS/?vd_source=9f1dda2358e63ace0b661e56fe417806</a>，程序為回映算法核心算法，可以修改此程序?qū)崿F(xiàn)基于試驗數(shù)據(jù)的本構(gòu)參數(shù)計算，不太會編程的可移步我的另外一個帖子，具體的<a href="https://www.yqgqt.org.cn/major/matlab" rel="noopener noreferrer" target="_blank">matlab程序</a>如下：</p>
展開 ABAQUS umat 非線性混合硬化本構(gòu)模型（Chaboche 硬化模型 ） ￥239
<p>本資源包含一份 PDF 文檔和可直接編譯運行的 Fortran UMAT 代碼，具體內(nèi)容為：</p><p>Chaboche硬化本構(gòu)模型 + 隱式積分 + 徑向返回</p><p>完整公式推導 + Fortran 源碼直接編譯</p><p>任意個數(shù)背應力分量 + 解析一致切線模量</p><p>PDF 包含規(guī)范化的本構(gòu)方程、隱式積分、徑向返回與一致切線模量推導，可供初學者學習。配套 UMAT 代碼可直接在 ABAQUS 編譯運行，采用全隱式積分搭配一致切線模量，收斂速度極快、計算精度極高，適合初學者快速入門。</p><p>下圖展示了部分PDF內(nèi)容，及umat計算結(jié)果與abaqus內(nèi)置模型對比，可以發(fā)現(xiàn)umat收斂速度極快，與abaqus內(nèi)置模型幾乎一致。
展開 平紋復合材料VUMAT子程序本構(gòu)介紹（hashin準則及線性損傷演化方法）
接下來主要介紹層內(nèi)的損傷本構(gòu)關(guān)系。
1. 平紋復合材料損傷剛度矩陣
（1）
（2）
（3）
其中，d代表損傷系數(shù)，L、T以及Z代表三個垂直的方向，t、c代表拉伸，壓縮損傷，例如dlt代表縱向拉伸損傷。
2. 損傷起始準則（hashin準則）
(4)
其中，f1t,f1c代表縱向纖維拉伸和壓縮損傷，f2t,f2c代表橫向纖維拉伸和壓縮損傷，f3代表厚度方向上的失效，其中應變大于0時為拉伸失效，小于0時為壓縮失效。
3. 損傷后損傷演化模型（線性退化模型）
（5）
當滿足損傷起始準則后，損傷演化開始起作用。其中1t、1c、2t、2c、3t、3c的失效模式下對應的損傷系數(shù)分別為d1t、d1c、d2t、d2c、d3t、d3c。其中，εii為當前應變，εiimax為初始失效應變，εif為最終失效應變。當i為1時，本構(gòu)關(guān)系如下圖1所示。
圖1. 線性退化示意圖
至此，平紋復合材料的本構(gòu)介紹完畢，隨后就是將這些公式利用fortran語言寫成VUMAT代碼。詳細VUMAT后臺私信。
展開 康國政老師的書《非線性本構(gòu)關(guān)系在ABAQUS中的實現(xiàn)》第7章的幾處疑似印刷錯誤
《非線性本構(gòu)關(guān)系在ABAQUS中的實現(xiàn)》是由闞前華、康國政、徐祥三位老師著的重要的關(guān)于采用abaqus編寫本構(gòu)關(guān)系的著作。該書出版后獲得了廣大abaqus本構(gòu)開發(fā)研究者的廣泛好評和推薦。
圖片來源（實體書拍照）
該書第一章描述了全量和增量兩種形式的本構(gòu)關(guān)系，本構(gòu)關(guān)系的張量表示，非線性方程組的求解策略和本構(gòu)關(guān)系的有限元實現(xiàn)過程及abaqus的用戶材料子程序接口，后續(xù)幾章內(nèi)容分別描述了非線性彈性，彈塑性，黏塑性，超彈性，循環(huán)彈塑性和循環(huán)黏塑性，耦合損傷循環(huán)塑性，大變形彈塑性循環(huán)本構(gòu)，晶體塑性循環(huán)和應變梯度塑性等內(nèi)容。涉及的本構(gòu)關(guān)系種類多樣，推導詳細，對筆者在開發(fā)彈塑性材料本構(gòu)方面有十分大的啟發(fā)。
本文主要針對該書第7章“循環(huán)彈塑性本構(gòu)關(guān)系”中的幾處疑似印刷錯誤進行討論。本文作者在彈塑性本構(gòu)方面尚未達到入門水平，無論是在本構(gòu)關(guān)系的知識廣度還是深度上，都與本書作者三位老師相距甚遠，本文指出的幾處疑似印刷錯誤僅僅從書中內(nèi)容印刷出發(fā)，不涉及本構(gòu)關(guān)系的具體理論糾正，指出的也僅僅是本人的個人看法，非常可能本人理解有誤，歡迎三位作者和其他讀者批評指正。
第一處：
第101頁式（7-26）為：
第102頁式（7-27）為：
從式（7-26）推導至式（7-27），式（7-27）中的分母有誤，實際上式7-27應當是：
第二處：
書中（7-29）式：
實際應為：
第三處：
書中第（7-30）式：
上式應為：
第4處：
上式應為：
第5處：
上式應為：
以上，即是《非線性本構(gòu)關(guān)系在ABAQUS中的實現(xiàn)》第7章循環(huán)彈塑性本構(gòu)關(guān)系的部分疑似印刷錯誤。
展開 基于回映算法的Chanboche各向同性非線性隨動硬化本構(gòu)編程重要參考資料 ￥599
<p><strong style="color: rgb(27, 27, 27); background-color: rgb(255, 255, 255);">Chaboche各向同性非線性隨動硬化行為的材料本構(gòu)模型計算matlab程序+</strong>基于回映算法的Chanboche各向同性非線性隨動硬化本構(gòu)<strong style="color: rgb(27, 27, 27); background-color: rgb(255, 255, 255);">本人編程所用的資料，沒有程序，只有資料，都是干貨</strong></p>
展開 
基于單個單元的有限元模型對Chaboche各向同性非線性隨動硬化本構(gòu)模型進行了仿真驗證 ￥149
<p>可以使用單個單元對計算出來的本構(gòu)進行驗證，這是對chaboche各向同性非線性隨動硬化本構(gòu)進行驗證，格式不被允許，下載后后綴改成<a href="https://www.yqgqt.org.cn/major/cae" rel="noopener noreferrer" target="_blank">cae</a>即可，abaqus2020版本以上打開，詳情可查看視頻https://www.bilibili.com/video/BV1Qc411p7E3/?vd_source=9f1dda2358e63ace0b661e56fe417806</p><div contenteditable="false" width="100%"><div><img src="https://img.jishulink.com/upload/202305/d126c60f514f41e499e1de172b8e5049.jpg" title="單個單元滯回環(huán)曲線.jpg" alt="單個單元滯回環(huán)曲線.jpg" style="max-width:760px;" data-mobile-src="https://img.jishulink.com/upload/202305/d126c60f514f41e499e1de172b8e5049.jpg?image_process=/format,webp/resize,w_400" data-pc-src="https://img.jishulink.com/upload/202305/d126c60f514f41e499e1de172b8e5049.jpg?
展開 如何在ANSYS中擬合橡膠材料曲線？ 附Ansys橡膠材料的粘彈性本構(gòu)模型下載
STEP 3：從hyperelastic模型本構(gòu)中拖動需要擬合的材料本構(gòu)模型到材料中，此時可以在材料橡膠本構(gòu)模型中發(fā)現(xiàn)curve fitting選項。
STEP 4：右鍵curve fitting，選擇solve curve fit，擬合好后，然后選擇copy calculated values to property，擬合參數(shù)便復制到定義的橡膠本構(gòu)模型中了。另外，擬合的曲線和實驗曲線均會在圖片中顯示出來，可以對比其重合度，測試哪種本構(gòu)更適合。
下載地址：Ansys橡膠材料的粘彈性本構(gòu)模型
展開 如何在ANSYS中擬合橡膠材料曲線？ 附Ansys橡膠材料的粘彈性本構(gòu)模型下載
STEP 3：從hyperelastic模型本構(gòu)中拖動需要擬合的材料本構(gòu)模型到材料中，此時可以在材料橡膠本構(gòu)模型中發(fā)現(xiàn)curve fitting選項。
STEP 4：右鍵curve fitting，選擇solve curve fit，擬合好后，然后選擇copy calculated values to property，擬合參數(shù)便復制到定義的橡膠本構(gòu)模型中了。另外，擬合的曲線和實驗曲線均會在圖片中顯示出來，可以對比其重合度，測試哪種本構(gòu)更適合。
下載地址：Ansys橡膠材料的粘彈性本構(gòu)模型
展開 【ANSYS】橡膠材料本構(gòu)擬合與拉扭試驗驗證
01 引子
橡膠材料是典型的超彈性材料，要獲取超彈性材料本構(gòu)模型（常見有Mooney-Rivlin、Ogden、Yeoh等），一般需要做一系列標準橡膠試驗并進行數(shù)據(jù)擬合。
本例演示了ANSYS對超彈性材料的曲線擬合能力，并通過有限元分析與拉扭試驗的對比，驗證所建立的本構(gòu)模型的有效性。
常見的橡膠標準拉伸試驗
02 案例介紹
現(xiàn)需要一個本構(gòu)模型來匹配硫化天然橡膠材料在各種變形模式下的100%工程應變的行為。
本例中，已通過試驗（單軸、雙軸和平面拉伸試驗）獲取了橡膠的實驗數(shù)據(jù)。使用這些數(shù)據(jù)，通過超彈性擬合能力確定本構(gòu)模型的參數(shù)，可以擬合3參、5參和9參的Mooney-Rivlin超彈性模型。
試驗數(shù)據(jù)
同時對橡膠進行了拉扭實驗（將條形試件的兩端夾入測試儀器中，然后將試樣拉伸到原尺寸長度的50%，并將試樣的一端扭四圈）。試樣與ASTM D1043中規(guī)定的試樣相似，如下圖所示：
拉扭試驗條形試件
使用擬合得出的Mooney-Rivlin超彈性模型（5參為例）對拉扭試驗就行有限元分析，并與試驗結(jié)果相對比，據(jù)此判斷前面擬合得出的本構(gòu)模型能否反映橡膠材料的真實行為。
模型采用SOLID186單元，兩端夾鉗區(qū)域采用MPC算法綁定到定位點。
有限元模型示意圖
按照拉扭試驗的加載順序：
step1：對兩端夾持區(qū)域施加試件厚度25%的壓縮位移，模擬夾具對試件的夾持作用。
step2：通過移動一側(cè)的夾持區(qū)域（剛性接觸面），同時固定另一側(cè)夾持區(qū)域，模擬拉伸到50%的拉伸狀況。
展開 