等效塑性應變ansys的案例
EE LE PE NE PEEQ(等效塑性應變) PEMAG(塑性應變量)
一）名詞解釋：
EE 彈性應變
NE 名義應變
LE 對數應變【即真應變，對于單軸拉伸LE=ln（1+NE）】
PE 塑性應變
PEEQ （equivalent plastic strain）等效塑性應變
PEMAG （Plastic strain magnitude） 塑性應變量
PEEQ與PEMAG的區別是
PEMAG描述的是變形過程中某一時刻的塑性應變，與加載歷史無關，而PEEQ是整個變形過程中塑性應變的累積結果。
如果一個圓桿受單向拉伸至屈服，再通過單向壓縮使其恢復初始長度，則最終的PEMAG為0，而PEEQ是拉伸和壓縮過程中塑性應變的絕對值之和。
二）CAE模型
下圖示同樣尺寸的鋁板，2種拉伸工況
板的拉伸及其1/4模型
工況（1） 單調拉伸
依據property材料參數，理論計算算的板拉伸板0.06m時屈服。Load單調拉伸到位移為0.09m，輸出EE， LE ，PE， PEEQ ，PEMAG ，NE
工況（2） 拉伸到屈服再壓縮到屈服，再拉伸到屈服（反復加載），最大位移和工況（1）一樣是0.09m。
展開 LS-DYNA——等效塑性應變
等效塑性應變
等效塑性應變是一個單調增加的標量值，它是作為變形率張量的塑性分量(Dp)ij的函數遞增計算的。在張量表示法中，表示為：
epspl=integral over time of (depspl)=integral[sqrt(2/3(Dp)ij*(Dp)ij)]*dt
只要材料屈服，即只要應力狀態在屈服面上，等效塑性應變就會增長。
應變張量?
相反，當在*DATABASE_EXTENT_BINARY中設置STRFLG為 1 時，由LS-DYNA寫出的張量應變值不一定是單調增加的，因為它們反映的是當前的 總變形 狀態（彈性+塑性）。在LS-PrePost中顯示繪制 應變張量 ，請單擊Fcomp> Strain。
以張量表示的等效應變為sqrt(2/3(eps)ij*(eps)ij)；（見2006年LS-DYNA理論手冊第461頁）。這與等效塑性應變不是一回事。
其它應變也可以在LS-PrePost中進行繪制顯示：
FCOMP>Infin；（無窮小或工程應變）
FCOMP>Green
FCOMP>Almansi
等效應力，也稱為馮-米塞斯應力，定義如下：
sigvm=1/sqrt(2)*sqrt[(sigx-sigy)^2+(sigy-sigz)^2+(sigz-sigx)^2+6*sigxy^2+6*sigyz^2+6*sigzx^2]
展開 等效塑性應變為什么為0
為什么我輸出的等效塑性應變一直是0啊
UMAT (各項同性+J2流動+自定義屈服強度等效塑性應變關系+歐拉后推徑向返回) ￥10
Abaqus自帶有3維的各項同性+J2流動+自定義屈服強度等效塑性應變關系+歐拉后推徑向返回的UMat例子
在此基礎上我進行了一些修訂用于以下情況（附件中包含for和inp）
1. 2維平面應變+各項同性+J2流動+自定義屈服強度等效塑性應變關系+歐拉后推徑向返回
2. 2維平面應變+各項同性+J2流動+冪硬化+歐拉后推徑向返回
冪硬化本構更新在張純禹的power-law基礎上修改得到，涉及到牛頓迭代的方式進行屈服應力求解
其原始文件，一起上傳
附件如下：
適用于ansys的應變梯度塑性本構（CMSG）子程序（開源資源）
/blob/f4680eb4fe4febb1c8f3a270e2a958663b52a978/Source/usermatps.F
該程序以ansys為開發平臺，但里面的很多內容是相通的。
