彈塑性變形行為

Abaqus默認的塑性材料特性應(yīng)用金屬材料的經(jīng)典塑性理論，采用Mises屈服面來定義各向同性屈服。

金屬材料的彈塑性變形行為可以簡述如下（見圖1）：在小應(yīng)變時，材料性質(zhì)基本為線彈性，彈性模量E為常數(shù)；應(yīng)力超過屈服應(yīng)力（yield stress）后，剛度會顯著下降，此時材料的應(yīng)變包括塑性應(yīng)變（plastic strain）和彈性應(yīng)變（elastic strain）兩部分；在卸載后，彈性應(yīng)變消失，而塑性應(yīng)變是不可恢復(fù)的；如果再次加載，材料的屈服應(yīng)力會提高，即所謂的加工硬化（work hardening）。

在單向拉伸/壓縮試驗中得到的數(shù)據(jù)通常是以名義應(yīng)變ε_nom和名義應(yīng)力σ_nom表示的，其計算公式為：

其中，Δι 是試樣的長度變化量，ι₀ 是試樣的初始長度，F 是載荷，A₀ 是試樣的初始截面面積。為了準確地描述大變形過程中截面面積的改變，需要使用真實應(yīng)變ε_true（又稱對數(shù)應(yīng)變）和真實應(yīng)力σ_true，它們與名義應(yīng)變ε_nom 和名義應(yīng)力σ_nom 之間的換算公式為：

其中：ι 是試樣的當前長度， A 是試樣當前的截面面積。對于拉伸試驗，上式中的ε_nom 是正值；對于壓縮試驗，ε_nom 是負值。

真實應(yīng)變ε_true 是由塑性應(yīng)變ε_pι 和彈性應(yīng)變ε_eι 兩部分構(gòu)成的。在Abaqus中定義塑性材料參數(shù)時，需要使用塑性應(yīng)變ε_pι ，其表達式如下： 

上述各個量在Abaqus分析結(jié)果中所對應(yīng)的變量為：

○ 真實應(yīng)力σ_true ：Mises應(yīng)力S, Mises；

○ 真實應(yīng)變 ε_true ：對于幾何非線性問題（*STEP, NLGEOM=Yes），Abaqus在ODB文件中默認的輸出變量是對數(shù)應(yīng)變LE；對于幾何線性問題（*STEP, NLGEOM=No），Abaqus默認的輸出變量是總應(yīng)變E；

○ 塑性應(yīng)變 ε_pι ：等效塑性應(yīng)變PEEQ，塑性應(yīng)變量PEMAG，塑性應(yīng)變分量PE；

○ 彈性應(yīng)變 ε_eι ：彈性應(yīng)變EE；

○ 名義應(yīng)變ε_nom ：名義應(yīng)變NE。

需要說明的是：在比例加載時（即加載過程中主應(yīng)力方向和比值不變），大多數(shù)材料的PEMAG和PEEQ相等。這兩個量的區(qū)別在于，PEMAG描述的是變形過程中某一時刻的塑性應(yīng)變，與加載歷史無關(guān)，而PEEQ是整個變形過程中塑性應(yīng)變的累積結(jié)果。例如，單向拉伸一個圓柱體，使其發(fā)生塑性應(yīng)變，再通過單向壓縮使其恢復(fù)初始長度，則最終的PEMAG為0，而PEEQ是拉伸和壓縮過程中塑性應(yīng)變的累積。

等效塑性應(yīng)變PEEQ大于0表明材料發(fā)生了屈服。在工程結(jié)構(gòu)中，等效塑性應(yīng)變一般不應(yīng)超過材料的破壞應(yīng)變（failure strain）。對于金屬成形等大變形問題，應(yīng)根據(jù)生產(chǎn)工藝要求來確定許可的等效塑性應(yīng)變量。

Abaqus/Standard無法準確模擬構(gòu)件因塑性變形過大而破壞的過程，此問題應(yīng)使用Abaqus/Explicit來分析。

彈塑性分析方法

在Abaqus中進行彈塑性分析時，最主要的操作是按照下面介紹的方法來定義塑性材料數(shù)據(jù)（即應(yīng)力應(yīng)變關(guān)系）。如果模型的位移較大，則應(yīng)將幾何非線性參數(shù)Nlgeom設(shè)置為ON，對應(yīng)的關(guān)鍵詞為*STEP, NLGEOM=Yes。

需要注意的是：彈塑性分析中并不一定總要考慮幾何非線性?！皫缀畏蔷€性”的含義是位移的大小對結(jié)構(gòu)的響應(yīng)發(fā)生影響，例如大位移、大轉(zhuǎn)動、初始應(yīng)力、幾何剛性化和突然翻轉(zhuǎn)等。

在Abaqus中定義塑性材料數(shù)據(jù)的方法是：按照公式2~4，將單向拉伸或壓縮試驗得到的名義應(yīng)力σ_nom 和名義應(yīng)變ε_nom 數(shù)據(jù)換算為真實應(yīng)力σ_true 和塑性應(yīng)變ε_pι ，然后在Abaqus/CAE或INP文件中給出一系列由真實應(yīng)力和塑性應(yīng)變所構(gòu)成的數(shù)據(jù)點，Abaqus將自動在各數(shù)據(jù)點之間進行線性插值。在Abaqus中定義塑性材料的關(guān)鍵詞為：

*MATERIAL,  NAME = <材料名稱 >

……

*PLASTIC

<屈服點處的真實應(yīng)力 >,  0

<真實應(yīng)力 >,  <塑性應(yīng)變 >

……

需要注意的是：關(guān)鍵詞 *PLASTIC下面第一行中的第二項數(shù)據(jù)必須為0，其含義為：在屈服點處的塑性應(yīng)變?yōu)?。如果此處的值不為0，在運行時會出現(xiàn)以下錯誤信息:“***ERROR: THE PLASTIC STRAIN AT FIRST YIELD MUST BE ZERO”。

例如，下面的語句定義了圖2所示的塑性材料數(shù)據(jù)（鋼材的單軸壓縮試驗）：

*Material, name=Material-1

*Elastic

210000., 0.3

*Plastic

418, 0

500, 0.01581

605, 0.02983

695, 0.056

780,    0.095

829,     0.15

882,     0.25

908,     0.35

921,     0.45

932,     0.55

955,     0.65

988,    0.75

1040,    0.85

本例中，最大塑性應(yīng)變?yōu)?.85，相應(yīng)的真實應(yīng)力為1040，其含義為：材料的Mises應(yīng)力達到1040后，材料變?yōu)槔硐胨苄裕▓D6-2中的虛線部分），即材料會持續(xù)變形，直到應(yīng)力降至小于或等于1040。換言之，在理想塑性狀態(tài)下，應(yīng)力和應(yīng)變值不是一一對應(yīng)的，而這有可能會造成收斂問題，因此，在設(shè)定關(guān)鍵詞 *PLASTIC的塑性數(shù)據(jù)時，應(yīng)盡可能讓其中最大的真實應(yīng)力和塑性應(yīng)變大于模型中可能出現(xiàn)的應(yīng)力應(yīng)變值。

需要注意的是：關(guān)鍵詞 *PLASTIC下面各個數(shù)據(jù)行中的塑性應(yīng)變必須按照遞增的順序排列，否則在運行時會出現(xiàn)以下錯誤信息:“***ERROR: THE INDEPENDENT VARIABLES MUST BE ARRANGED IN ASCENDING ORDER”。

比較分析發(fā)現(xiàn)：當塑性應(yīng)變較小時，由單向拉伸試驗和單向壓縮試驗得到的真實應(yīng)力-真實應(yīng)變關(guān)系曲線基本一致。當塑性應(yīng)變很大時，單向拉伸試驗中的試樣會出現(xiàn)頸縮，而單向壓縮試驗中摩擦力的影響變大，試樣變成圓鼓形，此時兩種試驗結(jié)果在塑性應(yīng)變很大時都不準確，用戶應(yīng)仔細考察大變形分析結(jié)果的準確性。

另外，對于同一個模型用戶可以混合使用彈塑性材料和線彈性材料。為縮短分析時間，可以只將所關(guān)心的重要部位設(shè)置為彈塑性材料，而將不重要的部位設(shè)置為線彈性材料（前提是這樣的設(shè)置不會影響對重要部位的分析精度）。

SIMULIA提供了一套先進的仿真產(chǎn)品組合，其中包括Abaqus、Isight、fe-safe、Tosca、Simpoe-Mold、SIMPACK、CST Studio Suite、XFlow、PowerFLOW等。SIMULIA社區(qū)是尋找SIMULIA軟件最新資源以及與其他用戶展開協(xié)作的最佳場所。作為打開創(chuàng)新思維和知識積累之門的鑰匙，SIMULIA社區(qū)為您提供了您需要的各種工具，隨時隨地幫助您擴大知識面。