ansys 彈塑性本構(gòu)的案例
并不簡單的彈塑性本構(gòu)子程序
在寫彈塑性本構(gòu)之前，我對(duì)塑性流動(dòng)是干嘛使的沒有直觀概念。寫的時(shí)候我才明白，由于只能先算出來等效塑性應(yīng)變，沒有流動(dòng)方向的話，就無法把它轉(zhuǎn)換到各個(gè)應(yīng)變分量，不知道應(yīng)變分量就無法計(jì)算應(yīng)力。這玩意從數(shù)學(xué)上講，是一個(gè)轉(zhuǎn)換公式。
我們目前重工業(yè)上大部分的結(jié)構(gòu)材料還是金屬，盡管ABAQUS中有自帶的JC模型，但是如果要模擬更復(fù)雜的情況，學(xué)會(huì)寫彈塑性本構(gòu)就十分必要。
本期就給一個(gè)彈塑性VUMAT拉伸失效的案例，結(jié)合單元?jiǎng)h除技術(shù)，模擬結(jié)構(gòu)破壞過程。
本構(gòu)模型
采用經(jīng)典老演員JC模型描述本案例的彈塑性本構(gòu)：
為了模擬結(jié)構(gòu)破壞，采用如下準(zhǔn)則判斷單元完全失效，滿足其一即可：
（1）材料Mises應(yīng)力達(dá)到極限值；
（2）材料極限應(yīng)變達(dá)到極限值。
子程序結(jié)構(gòu)
子程序的基本結(jié)構(gòu)如下：
1.初始化準(zhǔn)備工作
程序首先進(jìn)行初始化準(zhǔn)備工作，讀入材料的彈性參數(shù)、強(qiáng)度參數(shù)、硬化參數(shù)以及應(yīng)變率相關(guān)參數(shù)，然后構(gòu)建彈性剛度矩陣，為后續(xù)計(jì)算奠定基礎(chǔ)。
2.進(jìn)入材料點(diǎn)循環(huán)
接下來進(jìn)入材料點(diǎn)循環(huán)，對(duì)每個(gè)積分點(diǎn)逐一進(jìn)行計(jì)算。對(duì)于每個(gè)材料點(diǎn)，程序首先讀取上一步的狀態(tài)變量，包括累積的等效塑性應(yīng)變、應(yīng)力狀態(tài)以及背應(yīng)力等內(nèi)部變量。
3.失效判斷
程序隨后進(jìn)行失效判斷，檢查材料是否滿足失效準(zhǔn)則。判斷依據(jù)包括兩個(gè)方面：一是等效塑性應(yīng)變是否超過極限應(yīng)變閾值，二是等效應(yīng)力是否達(dá)到破壞強(qiáng)度。一旦滿足任一失效條件，程序?qū)⒉牧蠘?biāo)記為失效狀態(tài)，并大幅降低其剛度以模擬材料的承載能力喪失。
4.本構(gòu)響應(yīng)計(jì)算階段
在本構(gòu)響應(yīng)計(jì)算階段，程序考慮了應(yīng)變率效應(yīng)和材料硬化特性，更新當(dāng)前的屈服應(yīng)力。同時(shí)計(jì)算應(yīng)力偏量，得到米塞斯等效應(yīng)力和塑性流動(dòng)方向，這些是判斷材料是否屈服的關(guān)鍵參數(shù)。
5.彈塑性判別
然后進(jìn)行彈塑性判別。
展開 ABAQUS umat 理想彈塑性本構(gòu)模型 ￥99
<p class="ql-align-justify"><span style="color: rgb(15, 17, 21);">本資源包含一份 PDF 文檔和可直接編譯運(yùn)行的 Fortran UMAT 代碼，具體內(nèi)容為：</span></p><p class="ql-align-justify">理想彈塑性本構(gòu) + 隱式積分 + 徑向返回</p><p class="ql-align-justify">完整公式推導(dǎo) + Fortran 源碼直接編譯</p><p class="ql-align-justify">von Mises 屈服+ 一致切線模量全實(shí)現(xiàn)</p><p class="ql-align-justify">PDF 包含規(guī)范化的本構(gòu)方程、隱式積分、徑向返回與一致切線模量推導(dǎo)，可供初學(xué)者學(xué)習(xí)。配套 UMAT 代碼可直接在 ABAQUS 編譯運(yùn)行，采用全隱式積分搭配一致切線模量，收斂速度極快、計(jì)算精度極高，<span style="background-color: rgba(0, 0, 0, 0);">適合初學(xué)者快速入門。</span></p><p class="ql-align-justify"><span style="background-color: rgba(0, 0, 0, 0);">下圖展示了</span><span style="color: rgb(25, 27, 31);">部分</span><span style="background-color: rgba(0, 0, 0, 0);">PDF內(nèi)容，及umat計(jì)算結(jié)果與abaqus內(nèi)置模型對(duì)比，可以發(fā)現(xiàn)umat收斂速度極快，與abaqus內(nèi)置模型幾乎一致。
展開 材料本構(gòu)彈塑性力學(xué)知識(shí)三
在連續(xù)介質(zhì)力學(xué)中，所有問題(包括運(yùn)動(dòng)、應(yīng)力、應(yīng)變以及守恒定律等)既可用物體變形前的初始構(gòu)形B為參照構(gòu)形(取x1為自變量)來描述，又可用物體變形后的新構(gòu)形，B'為參照構(gòu)形(取x1*為自變量)來描述，前者稱為拉格朗日(LagrangeJ L)描述，后者稱為歐拉(Euler L)描述。
在固體力學(xué)中，我們常采用拉格朗日描述;在流體力學(xué)中采用歐拉描述更為方便;而對(duì)大變形問題及一般的物理定律，采用拉格朗日坐標(biāo)來建立它的數(shù)學(xué)表達(dá)式更為方便，但在求解具體問題時(shí)，又常以歐拉描述更方便，所以兩種描述都要采用。
—End—
CAE仿真與數(shù)值模擬微信公眾號(hào)，主要介紹CAE仿真與數(shù)值模擬的知識(shí)與應(yīng)用公眾號(hào)主要介紹CAE仿真與數(shù)值模擬的知識(shí)與應(yīng)用。通過論壇，博客，論文，案例等為大家?guī)碇R(shí)食糧。仿真軟件:abaqus、ansys、flunet、comsol、hypermesh、moldflow等，涉及領(lǐng)域有機(jī)械材料土木物理等。
展開 材料本構(gòu)彈塑性力學(xué)知識(shí)二
彈塑性材料：固體材料在受力后產(chǎn)生變形，從變形開始到破壞一般要經(jīng)歷彈性變形和塑性變形這兩個(gè)階段。根據(jù)材料力學(xué)性質(zhì)的不同，有的彈性階段較明顯，而塑性階段很不明顯，像鑄鐵等脆性材料，往往經(jīng)歷彈性階段后就破壞。有的則彈性階段很不明顯，從開始變形就伴隨著塑性變形，彈塑性變形總是耦連產(chǎn)生，像混凝土材料就是這洋。而大部分固體材料都呈現(xiàn)出明顯的彈性變形階段和塑性變形階段。今后我們主要是討論這種有彈性與塑性變形階段的固體材料，并統(tǒng)稱為彈塑性材料。
鮑辛格效應(yīng)：由于預(yù)加塑性拉伸荷載而使壓縮屈服應(yīng)力降低的現(xiàn)象稱為Bauschinger效應(yīng)。正是由于這種效應(yīng)，塑性變形時(shí)一種各向異性的過程，Bauschinger效應(yīng)是一種由塑性應(yīng)變引起的特殊的方向各向異性的形式，因?yàn)樵诤罄^逆向荷載作用下，一個(gè)方向的初始塑性變形會(huì)減小其反方向的屈服一個(gè)應(yīng)力。在多軸應(yīng)力情況下，與這種現(xiàn)象對(duì)應(yīng)的是具有不同方向屈服應(yīng)力之間的相互影響和橫向效應(yīng)，某一方向的預(yù)加應(yīng)變達(dá)到塑性范圍將會(huì)改變其所有方向的屈服應(yīng)力值。因此Bauschinger效應(yīng)對(duì)于多維問題更重要，包括荷載方向有明顯改變的復(fù)雜應(yīng)力歷史，比如應(yīng)力改變符號(hào)和循環(huán)荷載的情況。
彈性變形與塑性變形的區(qū)別：卸除載荷后。變形可以完全恢復(fù)，是彈性變形的基本特征，而變形的不可恢復(fù)性是塑性變形的基本特征。彈性與塑性的基本區(qū)別不在于它們的應(yīng)力一應(yīng)變關(guān)系是否線性。
例如，在比例極限與彈性極限之間的AB曲線段，應(yīng)力與應(yīng)變不再成比例，進(jìn)入了非線性階段，但在B點(diǎn)以前卸除載荷，變形仍將完成恢復(fù)，屬于彈性變形階段。因此，彈性和塑性的基本區(qū)別在于卸載后，是否保留一個(gè)永久變形(塑性應(yīng)變〕。
在彈性變形階段，應(yīng)力與應(yīng)變之間呈一一對(duì)應(yīng)的關(guān)系。
展開 
運(yùn)用ABAQUS軟件對(duì)冰材料彈塑性本構(gòu)模型改進(jìn)及驗(yàn)證（附源文件） ￥1300
<p class="ql-align-justify"><strong>內(nèi)容：</strong></p><p class="ql-align-justify">基于參考文獻(xiàn)通過ABAQUS建立了冰材料彈塑性本構(gòu)模型；對(duì)比已有試驗(yàn)，對(duì)比裂紋演化現(xiàn)象和沖擊載荷曲線，驗(yàn)證了冰材料本構(gòu)模型的有效性。</p><p class="ql-align-justify"><img src="https://img.jishulink.com/202507/attachment/7b0d26ab81f645dc98e8b15335447247.png" width="1027"></p><p class="ql-align-justify"><br></p><p class="ql-align-justify"><br></p><p class="ql-align-justify"><br></p><p class="ql-align-justify"><br></p><figure style="text-align: center;" class="ql-align-center"><img src="https://img.jishulink.com/202510/attachment/7cbe0c886d1d4de59fdee40d233200d8.png" style="" width="616" data-mobile-src="https://img.jishulink.com/202510/attachment/7cbe0c886d1d4de59fdee40d233200d8.png?
展開 材料本構(gòu)彈塑性力學(xué)知識(shí)一
彈塑性材料
大多數(shù)材料往往都同時(shí)具有彈性和塑性性質(zhì),特別是在塑性變形階段,變形中既有可恢復(fù)的彈性變形,又有不可恢復(fù)的塑性變形;因此有時(shí)又稱彈塑性材料
彈性設(shè)計(jì)方法:
是以彈性分析為基礎(chǔ)的結(jié)構(gòu)設(shè)計(jì),假定材料為理想彈性地,相應(yīng)地這種設(shè)計(jì)觀點(diǎn)便以分析結(jié)果的實(shí)際使用范圍作為設(shè)計(jì)的失效準(zhǔn)則,即認(rèn)為應(yīng)力[嚴(yán)格地說是應(yīng)力的某一函數(shù)值]達(dá)到一定限值[彈性界限],將進(jìn)入塑性變形階段時(shí),材料將破壞.
塑性設(shè)計(jì)方法:
結(jié)構(gòu)中如果有一處或一部分材料"破壞",則認(rèn)為結(jié)構(gòu)失效(喪失所規(guī)定的效用).由于一般的結(jié)構(gòu)都處于非均勻受力狀態(tài)。當(dāng)高應(yīng)力點(diǎn)或高應(yīng)力區(qū)的材料到達(dá)彈性界限時(shí)、結(jié)構(gòu)的大部分材料仍處于彈性界限之內(nèi);而實(shí)際材料在應(yīng)力超過彈性界限以后并不實(shí)際發(fā)生破壞，仍具有一定的繼續(xù)承受應(yīng)力(載荷)的能力，只不過剛度相對(duì)地降低。因此彈性設(shè)計(jì)方法不能充分發(fā)揮材料的潛力，導(dǎo)致材料的某種浪費(fèi)。實(shí)際上，當(dāng)結(jié)構(gòu)內(nèi)的局部材料進(jìn)入塑性變形階段，在繼續(xù)增加外載時(shí)，結(jié)構(gòu)的外力(應(yīng)力)分布規(guī)律與彈性階段不同，即所謂內(nèi)力(應(yīng)力)重分布；這種重分布總的是使內(nèi)力(應(yīng)力)分布更趨均勻，使原來處于低應(yīng)力區(qū)的材料承受更大的應(yīng)力，從而更好地發(fā)揮材料的潛力，提高結(jié)構(gòu)的承載能力。顯然，以塑性分析為基礎(chǔ)的設(shè)計(jì)比彈性設(shè)計(jì)更為優(yōu)越。但是，塑性設(shè)計(jì)允許結(jié)構(gòu)有更大的變形，以及完全卸載后結(jié)構(gòu)將存在殘余變形。因此，對(duì)于剛度要求較高及不允許出現(xiàn)殘余變形的場合、這種設(shè)計(jì)方法不適用。
展開 彈塑性本構(gòu)關(guān)系的部分推導(dǎo)【樣版】
彈塑性本構(gòu)關(guān)系的部分推導(dǎo)【樣版】
木材三維彈塑性本構(gòu)子程序開發(fā)
問題介紹
木材的本構(gòu)模型是采用連續(xù)體單元建模模擬木材彈塑性響應(yīng)的基礎(chǔ)，然而木材復(fù)雜的力學(xué)性質(zhì)常常為其本構(gòu)模型的建立帶來困難。木材力學(xué)性質(zhì)的復(fù)雜性主要表現(xiàn)在：
不同方向的強(qiáng)度值和剛度值各不相同；
2. 同一方向的抗拉強(qiáng)度和抗壓強(qiáng)度之間存在差異；
3. 不同形式荷載作用下材料的響應(yīng)不同，壓力作用下材料的表現(xiàn)以延性為主， 而拉力和剪力作用下材料的破壞呈脆性。
木材在復(fù)雜應(yīng)力狀態(tài)下的彈塑性本構(gòu)模型。以經(jīng)典彈塑性力學(xué)為框架，該本構(gòu)模型建立在如下四個(gè)基本假設(shè)的基礎(chǔ)之上：
木材在彈性階段是理想的橫觀各向同性材料；
2. 材料的屈服符合簡化的 Hashin 屈服準(zhǔn)則；
3. 材料在受拉和受剪屈服之前是理想線彈性的，屈服之后進(jìn)入塑性流動(dòng)階段；
4. 材料受壓初始屈服之前是理想線彈性的，屈服之后進(jìn)入應(yīng)變硬化階段，隨 著屈服面的轉(zhuǎn)移到達(dá)最終屈服面后進(jìn)入完全的塑性流動(dòng)。
二。子程序編寫流程
本工作室在三維hashing模型的基礎(chǔ)上，利用Abaqus軟件平臺(tái)，開發(fā)了完整的木材的彈塑性本構(gòu)umat子程序，包含木材完整的彈性、塑性、強(qiáng)化以及軟化階段。編寫子程序的流程如下：
三。結(jié)果驗(yàn)證
通過如下圖的木材模型進(jìn)行驗(yàn)證：
該模型在受壓、受剪及受拉的工況下，應(yīng)力應(yīng)變曲線如下所示：
該子程序還有以下特征：
能計(jì)算靜力非線性
2. 收斂性好
3. 能計(jì)算復(fù)雜應(yīng)力狀態(tài)
附件為本子程序參考的文獻(xiàn)，供大家學(xué)習(xí)探討~
2. 木材的力學(xué)性質(zhì)試驗(yàn)研究及數(shù)值模擬方法.pdf
最后，大家有關(guān)于Abaqus二次開發(fā)的相關(guān)需求可以添加管理員扣扣：3045552826，微信：CAE320，同時(shí)也歡迎大家關(guān)注“320科技工作室”的微信公眾號(hào)，掃一掃二維碼即可關(guān)注~~
展開 各向同性硬化von Mises率無關(guān)彈塑性本構(gòu)理論以及umat源代碼 ￥99
各向同性硬化von Mises率無關(guān)彈塑性本構(gòu)理論以及umat源代碼
1 本構(gòu)理論
1.1 率形式
對(duì)于各向同性線彈性材料，其本構(gòu)方程為：
式中假設(shè)了應(yīng)變張量可以分解為彈性應(yīng)變和塑性應(yīng)變兩部分：
因此塑性本構(gòu)的關(guān)鍵在于計(jì)算塑性應(yīng)變的演化。對(duì)于率無關(guān)彈塑性的本構(gòu)理論，需要確定以下三個(gè)部分：
(1)：屈服條件
(2)：流動(dòng)法則
(3)：硬化法則
在此采用的是 von Mises 屈服條件：
式中后繼屈服應(yīng)力是等效塑性應(yīng)變的函數(shù)：
流動(dòng)法則為：
式中流動(dòng)方向的表達(dá)式為：
硬化法則為：
1.2 Return-mapping算法
上述的本構(gòu)方程均為率形式。在增量步中，給定增量應(yīng)變：
首先假設(shè)該增量應(yīng)變?nèi)珵閺椥詰?yīng)變，計(jì)算試驗(yàn)狀態(tài)下的一些物理量：
試驗(yàn)狀態(tài)下的應(yīng)力
試驗(yàn)狀態(tài)下的屈服函數(shù)值：
利用該試驗(yàn)屈服函數(shù)值來判斷在該增量步下是否發(fā)生了塑性屈服。如果：
則說明試驗(yàn)狀態(tài)即為真實(shí)狀態(tài)，即可進(jìn)行更新：
反之則需要進(jìn)行塑性更正，即需要計(jì)算塑性乘子的增量，利用以下非線性方程組進(jìn)行計(jì)算：
可以將該非線性方程組簡化至一個(gè)非線性方程，過程如下，將該方程組中的第一式分解為球量和偏量兩部分：
因此可以計(jì)算應(yīng)力為：
將上式中的第二式整理得到：
可以得到兩個(gè)張量的方向相同：
因此偏應(yīng)力可以用試驗(yàn)狀態(tài)的信息表示出來：
代入到最后一個(gè)一致性方程中可得：
即可利用牛頓迭代法對(duì)上述非線性方程進(jìn)行求解，得到塑性乘子增量。
展開 如何將三維彈塑性本構(gòu)應(yīng)用于平面應(yīng)力問題中
1 本構(gòu)理論
本文講解如何將三維的率無關(guān)彈塑性理論應(yīng)用到平面應(yīng)力問題中。對(duì)于平面應(yīng)變和軸對(duì)稱問題，由于是相應(yīng)的應(yīng)變分量為0，因?yàn)榭梢灾苯邮褂萌S的本構(gòu)，只需將相應(yīng)的應(yīng)變分量設(shè)為0作為本構(gòu)的輸入即可。然后，對(duì)于平面應(yīng)力問題，是相應(yīng)的應(yīng)力分量為0，由于本構(gòu)是由應(yīng)變驅(qū)動(dòng)求得對(duì)應(yīng)的應(yīng)力，相應(yīng)應(yīng)力分量為0相當(dāng)于對(duì)系統(tǒng)施加了相應(yīng)的約束，因此三維的本構(gòu)理論不可直接應(yīng)用于平面應(yīng)力問題中，需要將相應(yīng)的約束考慮其中進(jìn)行求解。
1.1 平面應(yīng)力理論
對(duì)于線彈性情況，由三維本構(gòu)方程推導(dǎo)平面應(yīng)力方程如下：
1.2 應(yīng)力更新算法
采用一種嵌套迭代的方法進(jìn)行應(yīng)力更新。我們將平面外應(yīng)變?nèi)匀蛔鳛?em>本構(gòu)的輸入，此時(shí)可調(diào)用三維的本構(gòu)方程，得到對(duì)應(yīng)的應(yīng)力。如果得到的平面外應(yīng)力不為0，則使用牛頓迭代法對(duì)平面外應(yīng)變進(jìn)行更新，持續(xù)此過程，直至滿足平面應(yīng)力假設(shè)。
展開 有關(guān)隧道襯砌用cdp彈塑性損失本構(gòu)的問題
我在隧道二次襯砌中采用了cdp模型，并進(jìn)行了開挖模擬，但是開挖完成后損傷參數(shù)為0，沒有變化，這是為什么呢，很它這里的警告有關(guān)嗎，我看別的文獻(xiàn)里面參數(shù)值雖小，但是是有值的呀
各向同性彈塑性本構(gòu)的vumat源代碼：通過修改umat ￥99
在檢查時(shí)，傳入vumat的totalTime和stepTime都為0，根據(jù)用戶給定的本構(gòu)關(guān)系，程序進(jìn)行計(jì)算并得到初始的穩(wěn)定時(shí)間增量。如果這個(gè)穩(wěn)定時(shí)間增量太大，就會(huì)導(dǎo)致計(jì)算不穩(wěn)定（不收斂），所以需要給出彈性的計(jì)算過程，以保證得到一個(gè)比較合適的初始穩(wěn)定時(shí)間增量。
vumat與umat的更多對(duì)比見下圖
2 代碼修改
從umat的源代碼出發(fā)，作出相應(yīng)修改可以得到vumat。
首先增加應(yīng)力和應(yīng)變分量的轉(zhuǎn)換函數(shù)
!*******************************************************************************
! transfer_strain_vumat2umat：將vumat接口中的應(yīng)變變量轉(zhuǎn)化為umat接口的應(yīng)變變量
! vumat的應(yīng)變變量為張量應(yīng)變，且存儲(chǔ)順序不同
! 變量說明
! 輸入:
! strain_vumat : vumat接口中應(yīng)變張量的存儲(chǔ)方式(e11,e22,e33,e12,e23,e31)
!
! 輸出：
! strain_umat ：umat接口中應(yīng)變張量的存儲(chǔ)方式(e11,e22,e33,2*e12,2*e13,2*e23)
!
展開 隨動(dòng)硬化von Mises率無關(guān)彈塑性本構(gòu)理論以及umat源代碼 ￥99
</p><p>4 測試</p><p>4.1 一個(gè)單元加卸載測試</p><p>設(shè)置Abaqus自帶線性隨動(dòng)硬化的本構(gòu)為：</p><p><img src="https://img.jishulink.com/msimage/202402/500ecbc35544cd0f21d7b44893591563.png"></p><p>使用umat設(shè)置的材料參數(shù)為：</p><p><img src="https://img.jishulink.com/msimage/202402/75334bcebf74f8bcc98e8eebfe627740.png"></p><p>分別代表?xiàng)钍夏Ａ俊⒉此杀龋跏记?yīng)力，以及等效塑性應(yīng)變與隨動(dòng)屈服應(yīng)力的數(shù)據(jù)點(diǎn)。對(duì)于線性隨動(dòng)硬化模型，可以選取三個(gè)數(shù)據(jù)點(diǎn)，保證三點(diǎn)處于同一直線上，對(duì)最后一組數(shù)據(jù)點(diǎn)進(jìn)行一個(gè)特殊處理，可以選取一個(gè)很大的塑性應(yīng)變值，以保證計(jì)算過程中的等效塑性應(yīng)變都落在這三個(gè)數(shù)據(jù)點(diǎn)點(diǎn)，由此插值得到便滿足線性關(guān)系。
展開 Abaqus Umat (子程序4) 彈塑性本構(gòu),米塞斯模型(Mises Model) ￥10
1）米塞斯模型為經(jīng)典的彈塑性本構(gòu)，主要用來模擬金屬材料在外荷載作用下的彈塑性行為
2）具體為金屬在各向均勻受壓狀態(tài)下不會(huì)產(chǎn)生塑性變形，只有在剪切作用下會(huì)發(fā)生塑性變形
3該Fortran代碼為Abaqus的外接子程序(user subroutine)，可用于學(xué)習(xí)最簡單的彈塑性本構(gòu)的編寫過程
米塞斯模型的適用范圍及屈服面形狀
所編寫米塞斯模型UMAT子程序
abaqus三維復(fù)合材料彈塑性+漸進(jìn)損傷本構(gòu)模型-3D VUMAT ￥145
對(duì)于纖維增強(qiáng)復(fù)合材料的模擬，在<a href="/major/ABAQUS中，集成了二維Hashin失效準(zhǔn)則與多種損傷演化準(zhǔn)則，但缺少三維的復(fù)合材料本構(gòu)模型。
參考一篇已發(fā)表的SCI文章，使用Fortran語言建立三維平紋織物復(fù)合材料彈塑性、漸進(jìn)損傷本構(gòu)模型-Vumat子程序。平紋織物復(fù)合材料在1方向和2方向絲束性能近似相同。
該程序是博士期間學(xué)習(xí)復(fù)材子程序的小部分總結(jié)，編程結(jié)構(gòu)并不是非常漂亮及完美，但確保能順利運(yùn)行，且單元驗(yàn)證結(jié)果與理論公式一致，介意請勿拍。
程序中塑性迭代部分并非主流的牛頓-拉夫遜和梯度下降方法，但經(jīng)過驗(yàn)證能夠適用于該模型，介意請勿拍。
附件內(nèi)容：1. inp算例模型(低速?zèng)_擊工況，1/4模型，層間使用cohesive element) 2. 子程序 3 .使用方法 4.參考論文名稱
首先介紹該子程序的使用方法與效果
1. 在ABAQUS中建立三維復(fù)合材料模型，這里建立一個(gè)簡單的方塊。賦給材料方向，1,2方向分別表示絲束的方向，3方向表示垂直于1，2的方向，也就是面外方向。
2. 建立材料屬性
3. 建立顯示Explicit計(jì)算時(shí)間步，時(shí)間0.005，在場輸出中勾選輸出 SDV和 STATUS.
4. 劃分網(wǎng)格，賦給Explicit 3D stress單元類型，邊界條件根據(jù)需要設(shè)定即可，此處不再贅述。此處劃分為一個(gè)單元，使用12方向往復(fù)加載卸載。建立Job，提交模型前在Job中選擇該子程序，設(shè)置雙精度計(jì)算。
5. 查看結(jié)果，等效塑性應(yīng)變在卸載時(shí)沒有變化，再次加載時(shí)剪切應(yīng)力按照原來的路徑返回，剪切損傷在卸載時(shí)也保持不變。
6. 將該子程序應(yīng)用在低速?zèng)_擊模型中，可以順利運(yùn)行。
接下來簡要介紹該子程序的相關(guān)理論，子程序、參考的論文名稱以及輸入材料參數(shù)的對(duì)應(yīng)含義打包在附件中。
展開 