天堂亚洲精品_国语精品视频_精品中文字幕在线_韩国在线视频一区_亚洲校园激情春色_国产精品久久久_91麻豆精品久久久久蜜臀_国产一区二区不卡_久久久久久久久久久久电影_国产成人一区二区

首次發(fā)布 | 模擬表演應(yīng)用程序

簡(jiǎn)介：大家好，我是模擬秀專欄作者九千。 我畢業(yè)于上海交通大學(xué)，獲得博士學(xué)位。 我目前從事機(jī)械仿真研究。 我在有限元模擬方面有 9 年的經(jīng)驗(yàn)。 擅長(zhǎng)材料本構(gòu)開發(fā)、ABAQUS子程序開發(fā)、ABAQUS-Python 2.0二次開發(fā)、Hypermesh-ABAQUS聯(lián)合仿真。 今后我將在Simulation Show為您提供Abaqus原創(chuàng)視頻課程、培訓(xùn)和技術(shù)咨詢服務(wù)。

眾所周知，Abaqus提供了豐富的子程序接口，其中UMAT、UHARD、UHYPER等子程序接口用于實(shí)現(xiàn)用戶自定義的材料本構(gòu)結(jié)構(gòu)。 但UMAT等子程序的編程開發(fā)需要用戶具有連續(xù)介質(zhì)力學(xué)的專業(yè)知識(shí)儲(chǔ)備。 以及一定的 Fortran 編程技能。 市面上有很多介紹UMAT等素材用戶子程序開發(fā)的文檔和書籍，但往往不夠系統(tǒng)，尤其是理論公式的推導(dǎo)過于跳躍，很少解釋理論公式與代碼的對(duì)應(yīng)關(guān)系。 另一方面，自媒體時(shí)代，仿真領(lǐng)域涌現(xiàn)了大量的教學(xué)視頻，但UMAT相關(guān)的講解視頻卻很少。

本文基于 Abaqus 各向同性硬化彈塑性 UMAT 示例。 本例收錄于仿真秀官網(wǎng)-九千CAE制作的UMAT系列視頻《Abaqus UMAT材質(zhì)子程序開發(fā)59講：經(jīng)典材質(zhì)本構(gòu)UMAT編程實(shí)例詳解》。 ，歡迎大家點(diǎn)擊文末閱讀原文進(jìn)行預(yù)覽。 以下為正文：

1. 各向同性彈塑性本構(gòu)

1. 彈性本構(gòu)

各向同性線彈性本構(gòu)是廣義胡克定律：

2. 塑料構(gòu)成

各向同性硬化塑性本構(gòu)采用經(jīng)典的米塞斯屈服準(zhǔn)則廣義胡克定律，即材料進(jìn)入屈服狀態(tài)后，米塞斯等效應(yīng)力等于隨后的屈服應(yīng)力：

利用屈服函數(shù)作為塑性勢(shì)函數(shù)（關(guān)聯(lián)流動(dòng)定律），可以建立應(yīng)力應(yīng)變關(guān)系，即等效應(yīng)變率與應(yīng)力偏差成正比：

等效應(yīng)變率在時(shí)域的積分即為等效應(yīng)變率：

2. 離散形式的應(yīng)力求解

在彈塑性UMAT的實(shí)現(xiàn)中，ABAQUS求解器將傳入增量步的總應(yīng)變?cè)隽繌埩俊?在UMAT中，我們需要根據(jù)彈塑性本構(gòu)量求解彈性應(yīng)變?cè)隽繌埩亢退苄詰?yīng)變?cè)隽繌埩浚缓笄蠼鈶?yīng)力增量來完成應(yīng)力更新。 對(duì)于彈塑性本構(gòu)，通常采用試驗(yàn)應(yīng)力-徑向回歸法來實(shí)現(xiàn)上述求解過程。 為此，首先假設(shè)當(dāng)前應(yīng)變?cè)隽烤鶠閺椥詰?yīng)變?cè)隽?strong>廣義胡克定律，由此計(jì)算出的應(yīng)力稱為試驗(yàn)應(yīng)力。 通過將試驗(yàn)應(yīng)力與當(dāng)前材料狀態(tài)對(duì)應(yīng)的后續(xù)屈服強(qiáng)度進(jìn)行比較，可以判斷材料是否發(fā)生屈服。 如果試驗(yàn)應(yīng)力小于隨后的屈服強(qiáng)度，則試驗(yàn)應(yīng)力即為所需應(yīng)力。 否則，需要采用徑向回歸法求解。 經(jīng)過推導(dǎo)，我們可以得到米塞斯應(yīng)力與試驗(yàn)應(yīng)力的等效應(yīng)力之間的關(guān)系：

將上式代入收益率函數(shù)可得：

進(jìn)一步，通過牛頓迭代法求解上式得到等效塑性應(yīng)變?cè)隽浚渌兞靠赏ㄟ^下式更新：

3.雅可比矩陣

對(duì)于彈性增量步長(zhǎng)，雅可比矩陣為

對(duì)于彈塑性增量步，需要對(duì)應(yīng)力應(yīng)變關(guān)系進(jìn)行微分，得到雅可比矩陣：

在

4. UMAT實(shí)現(xiàn)彈塑性本構(gòu)

基于上述應(yīng)力求解理論方法和導(dǎo)出的雅可比矩陣，通過Fortran實(shí)現(xiàn)UMAT子程序。 對(duì)應(yīng)的流程圖如下：

為了調(diào)用UMAT計(jì)算，需要在ABAQUS中建立與子程序相匹配的相應(yīng)材料。 相應(yīng)設(shè)置如下（通過User Material設(shè)置本構(gòu)參數(shù)，通過Depvar指定狀態(tài)變量個(gè)數(shù)）：

由于ABAQUS內(nèi)置了彈塑性本構(gòu)模型，因此選擇相同的本構(gòu)參數(shù)和有限元模型，利用內(nèi)置模型（圖左）和編寫的子程序（圖右）進(jìn)行求解問題。 結(jié)果如下圖所示。 兩者計(jì)算得到的米塞斯應(yīng)力與等效塑性應(yīng)變是一致的，說明了子程序的準(zhǔn)確性。

4、ABAQUS UMAT材質(zhì)子程序開發(fā)59講

UMAT不僅可以實(shí)現(xiàn)上述各向同性硬化彈塑性選區(qū)，還可以實(shí)現(xiàn)線彈性、超彈性、后續(xù)硬化彈塑性等選區(qū)。 如果有人想了解更多本文理論公式和UMAT代碼示例的詳細(xì)推導(dǎo)過程，請(qǐng)點(diǎn)贊或在文末觀看。 截圖發(fā)送至本公眾號(hào)并回復(fù)【UMAT】即可獲贈(zèng)作者原創(chuàng)視頻課程《ABAQUS UMAT Materials》《子程序開發(fā)59講》第一章全部資料以及umat-isotropic-hardening-plas2第 4 章中的 .html 信息。

以下是我的課程《ABAQUS UMAT材質(zhì)子程序開發(fā)講座59：經(jīng)典材質(zhì)本構(gòu)UMAT編程詳解實(shí)例》的大綱，強(qiáng)烈推薦大家學(xué)習(xí)。

限時(shí)特價(jià)（僅限10人）

優(yōu)惠券后價(jià)格：399元（價(jià)值：599元）

課程可隨時(shí)回放，可開具發(fā)票

講師提供VIP群知識(shí)圈問答及模型下載

掃描二維碼立即嘗試

作者以經(jīng)典的線彈性本構(gòu)、Neo-Hookean超彈性本構(gòu)、Mooney-Rivlin超彈性本構(gòu)、彈塑性本構(gòu)、線性運(yùn)動(dòng)硬化彈塑性本構(gòu)為例，詳細(xì)講解了與上述選區(qū)相關(guān)的編程。 實(shí)現(xiàn)了理論公式的推導(dǎo)和具體的UMAT編程實(shí)例，旨在幫助機(jī)械仿真領(lǐng)域的初學(xué)者達(dá)到ABAQUS材料子程序開發(fā)的入門水平。 還為學(xué)習(xí)者提供了答疑解惑的知識(shí)圈和VIP交流群。

1.你會(huì)得到

（1）掌握各種齊次線彈性UMAT的開發(fā)能力

（2）掌握Neo-Hookean超彈性UMAT開發(fā)能力

（3）掌握UHYPER實(shí)現(xiàn)Mooney-Rivlin超彈本構(gòu)的開發(fā)能力

（4）掌握各向同性硬化彈塑性UMAT的開發(fā)能力

（5）掌握UHARD實(shí)現(xiàn)各種各向同性硬化彈塑性本構(gòu)的開發(fā)能力

（6）掌握線性運(yùn)動(dòng)硬化彈塑性UMAT開發(fā)能力

（7）提供所有課程相關(guān)資料及VIP用戶交流群等服務(wù)

2. 哪些人適合學(xué)習(xí)？

（1）力學(xué)相關(guān)專業(yè)的學(xué)生和工程師

(2)學(xué)習(xí)結(jié)構(gòu)設(shè)計(jì)仿真工程師

(3) ABAQUS軟件用戶和學(xué)習(xí)者

（4）所有對(duì)UMAT子程序二次開發(fā)感興趣的工程師

（超過）