變形梯度的案例
結構有限元計算中的Locking和F-Bar單元
比如說在做熱-變形耦合分析時，溫度（在三角形單元中線性分布）和應變（在三角形單元中均勻分布）計算精度帶來的誤差。
對于該問題的進一步解釋建議參看文獻[2]，在哪里給出了通俗和準確的解釋。
1.2 如何減輕鎖死
可以通過單元技術來減輕鎖死的影響。如Selective reduce integration, Incompatible，Assumed strain，Mixed formulation，MITC等等。對于近似不可壓縮材料的體積鎖死，最常見的對應方法是采用B-bar單元，這種單元的新式變種是F-bar單元，這是本文下面要講的內(nèi)容。
2. B-bar單元[3,4]
一般的，單元內(nèi)任意一點的應變ε表為單元節(jié)點i的位移u的函數(shù)。在B-bar單元中，將應變分解，其中為單元的平均體積應變或單元中心的體積應變值為偏差應變值。相應地，單元內(nèi)任意一點的應變用下式計算
這里將體積應變做平均計算的方法緩和了體積鎖死。這種方法實現(xiàn)起來非常簡單，計算費用不高，效果也不錯，因此得到非常廣泛的應用。但是這種方法視乎不能滿足LBB條件[5]。
3. F-bar單元[6,7]
相應于B-bar單元對應變的加法分解，單元中的任意一點的變形梯度（deformation gradient）可以分解為體積變形梯度和偏差變形梯度兩個部分 F=Fs Fv。 理論上來說這種分解方式是嚴密的，而B-bar單元中的應變加法分解是一種近似，只是在應變很小時與變形梯度的乘法分解相近。因此把B-bar單元推廣到 F-bar單元是非常自然的，在變形較大時采用F-bar單元精度更好。
在F-bar單元中，單元中的任意一點的變形梯度表述為，J=detF. 其中
可以是單元中Jocobian J的平均值 ，也可以指單元中心的Jocobian值。
展開 一種新穎的多尺度晶體塑性實現(xiàn)方案-------Direct FE2
通過耦合積分點和多晶RVE模型實現(xiàn)尺度的模擬效果，宏觀模型的積分點提供變形梯度用于微觀RVE模型的邊界條件，微觀模型通過邊界條件計算應力，狀態(tài)變量，并返回一致性雅可比矩陣，模擬效果如下：
Direct FE2 對應的二維模型和三維模型如下圖所示
施加X方向的單軸拉伸，二維和三維的變形結果如下圖所示：
二維模擬效果：
三維模擬效果：
兩類多晶平均場均勻化方案
（一）taylor均勻化方案
Taylor均勻化模型在abaqus實現(xiàn)時，通常一個積分點包含多個晶粒，其中積分點響應由晶粒響應的體積平均響應給出，積分內(nèi)的各個晶粒的變形梯度均等于宏觀變形梯度。通過變形梯度計算各個晶粒的應力響應，晶粒彼此之間不存在應力平衡關系，但應變協(xié)調(diào)，計算得到各個晶粒的得到柯西應力，之后得到體平均柯西應力，更新積分點響應。以及整體的織構，使用該模型通常用于反應多晶聚集體整體的響應，以及織構演化分析。并且由于等應變假設，Taylor模型通常相較于有限元方法，整體的應力響應會更高，織構演化預測更加鋒銳。但計算成本顯著低于有限元離散。
以FCC軋制軋制模擬為例子（并行過程織構演化以及應力應變響應情況）
（二）粘塑性自洽方案（VPSC）
VPSC 全稱 Visco Plastic Self Consistent，指的是特定的機械狀態(tài) (VP) 和使用的方法 (SC)。VPSC 是為應用于低對稱材料（六邊形、三角形、正交、三角形）而開發(fā)的，理論的基本出發(fā)點與前述的晶體塑性本構一致，不同的是在于處理晶粒之間相互作用時，不適用taylor等應變假設，而是將晶粒視為無限大基體的夾雜，基體的應力與應變響應等于各個晶粒的平均響應，滿足應變協(xié)調(diào)，并通過制定晶粒之間相互作用方案（切線，仿射，割線……）滿足晶粒與基體的應力平衡，特別適用于大規(guī)模成型過程中的取向演化分析，尤其是hcp對稱性較差的結構，但無法考慮晶粒與晶粒之間以及晶粒內(nèi)部的響應情況分析，往往用于定性分析滑移系統(tǒng)開動情況，孿晶演化，宏觀應變硬化等。
展開 ABAQUS中的單元選擇
特點：在一階單元中引入一個增強單元變形梯度的附加自由度。這種對變形梯度的增強允許一階單元在單元域上對于變形梯度有一個線性變化。
不足：對單元的扭曲很敏感，在使用時必須小心以確保單元扭曲是非常小的。
（4）雜交單元：
應用：當材料行為是不可壓縮（泊松比=0.5）或非常接近于不可壓縮（泊松比>0.475）時，如橡膠材料，采用雜交單元。
特點：對于具有不可壓縮材料性質(zhì)的任何單元，一個純位移的數(shù)學公式是不適宜的，壓應力不能由節(jié)點位移計算。雜交單元包含一個可以直接確定單元壓應力的附加自由度，節(jié)點的位移場則主要用來計算偏應變和偏應力。
在選取單元類型時需要綜合考慮以下幾方面的問題：
1）如果需要得到的是節(jié)點應力，盡量不要選用線性減縮積分單元；
2）如果使用線性減縮積分單元，應注意避免出現(xiàn)沙漏模式，常采用網(wǎng)格細化來解決；
3）在定義了接觸和彈塑性材料的區(qū)域后，不要使用C3D20、C3D8R、C3D10等二次單元；
4）完全積分單元容易出現(xiàn)剪切閉鎖和體積閉鎖問題，一般情況下盡量不要使用；
5）對于ABAQUS/Standard分析，如果能夠劃分四邊形(Quad)或六面體(Hex)網(wǎng)格，建議盡量使用非協(xié)調(diào)單元(如C3D8I)，同時注意保證關鍵部位的單元形狀是規(guī)則的。
6）如果無法劃分六面體(Hex)網(wǎng)格，則應使用修正的二次四面體單元(C3D10M)，它適用于接觸和彈塑性問題，只是計算代價較大；
7）有些適用于ABAQUS/Standard分析的單元類型不能用于ABAQUS/Explicit分析中(例如非協(xié)調(diào)單元)。
展開 
ABAQUS中實體單元的應用
特點：在一階單元中引入一個增強單元變形梯度的附加自由度。這種對變形梯度的增強允許一階單元在單元域上對于變形梯度有一個線性變化。
不足：對單元的扭曲很敏感，在使用時必須小心以確保單元扭曲是非常小的。
（4）雜交單元：
應用：當材料行為是不可壓縮（泊松比=0.5）或非常接近于不可壓縮（泊松比>0.475）時，如橡膠材料，采用雜交單元。
特點：對于具有不可壓縮材料性質(zhì)的任何單元，一個純位移的數(shù)學公式是不適宜的，壓應力不能由節(jié)點位移計算。雜交單元包含一個可以直接確定單元壓應力的附加自由度，節(jié)點的位移場則主要用來計算偏應變和偏應力。
基于ABAQUS中如此豐富詳細的實體單元劃分，在使用時應尤其注意。
對于三維問題應盡量地采用六面體單元(磚型)。它們會以最低的成本給出最好的結果。當幾何形狀復雜時，可采用四面體單元和楔形單元。這些單元C3D4和C3D6的一階模式是較差的單元(需要細化網(wǎng)格以取得較好的精度)。
某些前處理包含了自由劃分網(wǎng)格算法，用四面體單元劃分任意幾何體的網(wǎng)格。對于小位移無接觸的問題，在ABAQUS/Standard中的二次四面體單元(C3D10)能夠給出合理的結果。這個單元的另一種模式是修正的二次四面體單元(C3D10M)，它適用于ABAQUS/Standard和ABAQUS/Explicit，對于大變形和接觸問題，這種單元是強健的，展示了很小的剪切和體積自鎖。但是，無論采用何種四面體單元，所用的分析時間都長于采用等效網(wǎng)格的六面體單元。
展開 晶體塑性有限元仿真入門(2)--BCC、FCC、HCP晶格材料以及多相材料的有限元模擬
其他模塊的設置與第1節(jié)"FCC晶格材料的變形模擬-單晶體"的創(chuàng)建過程類似，在此不再贅述。
后處理界面
應力應變分布如圖3.4所示。
圖3.4 后處理應力應變分布
4. BCC晶格材料的變形模擬-多晶體
所有模塊的設置與第2節(jié)"FCC晶格材料的變形模擬-多晶體"的創(chuàng)建過程類似，只是需要按照第3節(jié)的操作修改BCC晶格材料的材料參數(shù)，構建好的inp文件見附件，下面看看模型的變形結果。
后處理界面
應力應變分布：
圖4.1 后處理應力應變分布
圖4.2 后處理一些SDV結果的分布
5. HCP晶格材料的變形模擬-單晶體
塑性部分考慮滑移和孿晶的變形梯度如圖5.1所示，總的塑性變形梯度來自等式右側的兩項：分別代表由滑移和孿晶對塑性變形梯度的貢獻。與考慮滑移和孿晶的變形相比，滑移為主的塑性變形梯度忽略了孿晶對塑性變形梯度的貢獻，如圖5.2所示。
圖5.1 包含滑移、孿晶的變形梯度分解示意圖
圖5.2 滑移變形為主的梯度分解示意圖
Huang的晶體塑性子程序是基于滑移原理編寫的，因此使用該子程序進行HCP晶格材料變形的模擬時需保證變形是以滑移為主的前提條件，如論文《純鈦塑性變形行為的晶體塑性有限元模擬》、《單晶純鈦的細觀力學性能模擬》和《鈦合金雙態(tài)組織高溫拉伸行為的晶體塑性有限元研究》里明確指出忽略了HCP晶格材料的孿生變形。
如考慮孿生和滑移兩者的貢獻，需修改Huang基于滑移原理的子程序，或者換DAMASK平臺進行模擬， 鏈接里是DAMASK平臺模擬鎂合金的變形。
圖5.3 DAMASK平臺模擬鎂合金變形的參數(shù)設置
下面，將采用Huang的基于滑移原理晶體塑性子程序，開展HCP晶格材料以滑移為主的塑性變形模擬。
展開 如何才能選出適合于分析的單元類型?
特點：在一階單元中引入一個增強單元變形梯度的附加自由度。這種對變形梯度的增強允許一階單元在單元域上對于變形梯度有一個線性變化。
不足：對單元的扭曲很敏感，在使用時必須小心以確保單元扭曲是非常小的。
（4）雜交單元：
應用：當材料行為是不可壓縮（泊松比=0.5）或非常接近于不可壓縮（泊松比>0.475）時，如橡膠材料，采用雜交單元。
特點：對于具有不可壓縮材料性質(zhì)的任何單元，一個純位移的數(shù)學公式是不適宜的，壓應力不能由節(jié)點位移計算。雜交單元包含一個可以直接確定單元壓應力的附加自由度，節(jié)點的位移場則主要用來計算偏應變和偏應力。
在選取單元類型時需要綜合考慮以下幾方面的問題：
1）如果需要得到的是節(jié)點應力，盡量不要選用線性減縮積分單元；
2）如果使用線性減縮積分單元，應注意避免出現(xiàn)沙漏模式，常采用網(wǎng)格細化來解決；
3）在定義了接觸和彈塑性材料的區(qū)域后，不要使用C3D20、C3D8R、C3D10等二次單元；
4）完全積分單元容易出現(xiàn)剪切閉鎖和體積閉鎖問題，一般情況下盡量不要使用；
5）對于ABAQUS/Standard分析，如果能夠劃分四邊形(Quad)或六面體(Hex)網(wǎng)格，建議盡量使用非協(xié)調(diào)單元(如C3D8I)，同時注意保證關鍵部位的單元形狀是規(guī)則的。
6）如果無法劃分六面體(Hex)網(wǎng)格，則應使用修正的二次四面體單元(C3D10M)，它適用于接觸和彈塑性問題，只是計算代價較大；
7）有些適用于ABAQUS/Standard分析的單元類型不能用于ABAQUS/Explicit分析中(例如非協(xié)調(diào)單元)。
文章轉(zhuǎn)自有限元在線博客，分享給大家學習交流
展開 基于黃umat梯度結構晶粒變形模擬------案例十九 ￥99
? 基于黃umat梯度結構晶粒變形模擬
案例實操
1，建立包含896個晶粒的梯度多晶模型
2，對多晶模型賦予對應的材料屬性
3，X0方向固定，施加X1方向50%工程應變的拉伸載荷
4，保留晶界形狀，使用CPE3單元
5，提交與后處理材料數(shù)據(jù)
梯度晶粒幾何模型
模型載荷示意圖
不同時刻材料的對數(shù)應變分布
不同時刻材料的應力分布
材料的等效塑性應變的分布
根據(jù)應力應變分布情況可以清晰的看出，梯度晶粒結構應力應變分布更加均勻，不容易集中于某些區(qū)域，從而避免更早的發(fā)生頸縮失效，提高材料的延性。從而提高材料的服役壽命。
展開 Acta Mater：梯度納米結構高強高韌金屬動態(tài)變形機理
均勻結構由于其結構均勻，動態(tài)變形和剪切帶演化機理相對成熟，但缺乏加工硬化能力，在動態(tài)條件下缺乏均勻塑性。近年來研究表明，多級結構是高強度金屬突破其塑性瓶頸的一個重要途徑。
多級結構在準靜態(tài)條件下具有優(yōu)于均勻結構的力學性能，那么其在動態(tài)條件下是否同樣具有優(yōu)越的力學性能？各級結構及其協(xié)調(diào)變形如何影響動態(tài)力學性能？多級結構動態(tài)變形行為的微結構機理是什么？近期，中國科學院力學研究所、北卡州立大學、約翰霍普金斯大學的科研人員合作，在以上科學問題的研究中取得進展。相關研究成果發(fā)表在Materials Research Letters和Acta Materialia上。
論文鏈接：
https://www.sciencedirect.com/science/article/pii/S1359645418301174
https://www.tandfonline.com/doi/full/10.1080/21663831.2017.1334715
針對梯度結構，科研人員設計了一套新的動態(tài)剪切試驗手段，首次揭示了梯度納米結構的動態(tài)剪切變形機理：由于各層之間在動態(tài)變形過程中發(fā)生應變分配，產(chǎn)生了額外的加工硬化，能夠延遲剪切帶在納米晶表層的萌生，以及限制剪切帶從表面到芯部的擴展（其傳播速度相比均質(zhì)結構低一個數(shù)量級），梯度納米結構金屬能夠獲得比均質(zhì)結構優(yōu)越的動態(tài)剪切性能，同時發(fā)現(xiàn)廣為人知的剪切帶萌生的最大應力準則在梯度結構中不再適用。
圖1.梯度納米結構的動態(tài)剪切變形機理
通過冷軋和低溫短時退火，科研人員在低層錯能金屬中熵合金中獲得多尺度晶粒結構，研究發(fā)現(xiàn)多尺度晶粒之間的變形協(xié)調(diào)和應變分配能夠促進加工硬化，動態(tài)變形過程中發(fā)生了晶粒細化，能夠延緩剪切帶的萌生，促進動態(tài)剪切塑性，獲得了迄今為止報道的最優(yōu)越的動態(tài)剪切性能。
展開 . : 動態(tài)化學鍵助推水凝膠梯度結構、三維變形重復設計
自變形水凝膠在軟機器人、柔性電子等領域具有重大應用價值，如何控制凝膠內(nèi)部組成與應力分布是實現(xiàn)水凝膠三維變形的關鍵。通過光刻法可以構筑復雜的梯度結構，從而實現(xiàn)可控三維變形。但是，水凝膠網(wǎng)絡的梯度結構通常隨著化學反應的完成而永久固定。要實現(xiàn)不同的形狀變化并獲得新的三維構型，需要重新制備含有不同梯度結構的水凝膠。
浙江大學高分子系鄭強教授、吳子良研究員團隊通過引入動態(tài)化學鍵，實現(xiàn)了水凝膠梯度結構以及三維構型的重復調(diào)控。該團隊利用膠束聚合將疏水性香豆素單元引入親水性聚電解質(zhì)網(wǎng)絡，得到的水凝膠材料能夠通過光照射實現(xiàn)網(wǎng)絡結構的可逆調(diào)控。其中，香豆素單元可以在365 nm、254 nm紫外光下發(fā)生二聚、解聚反應（圖1）。
在水凝膠中，基于光致可逆交聯(lián)實現(xiàn)梯度結構與三維變形的重復設計需要解決以下問題：
（1）疏水性單元在水凝膠中的濃度較低，難以保證較高的光化學反應效率；
（2）解交聯(lián)導致凝膠劇烈溶脹，難以回復至初始狀態(tài)，不利于梯度結構的完全擦除。
針對以上問題，該團隊選取了十六烷基三甲基氯化銨（CTAC）作為表面活性劑，聚合后帶相反電荷的膠束與聚丙烯酸鏈段形成聚電解質(zhì)/表面活性劑復合物（PESC），提高香豆素單元的局部濃度及光化學反應效率。同時，PESC作為物理交聯(lián)，有效提高了水凝膠的力學性能。PESC的強度具有pH響應性，在中性條件下較弱，有利于體現(xiàn)香豆素不同狀態(tài)（交聯(lián)、解交聯(lián)）對凝膠溶脹度的影響；PESC在酸性條件下較強，可以消除含有不同狀態(tài)香豆素水凝膠溶脹度的差異。
圖1.
展開 有限變形晶體塑性快速傅里葉變化CPFFT實現(xiàn)
有限應變運動學 (Finite Strain Kinematics)
在有限變形框架下，總變形梯度被分解為彈性和塑性兩部分。文章強調(diào)了在參考構型下求解第一類 Piola-Kirchhoff 應力平衡的重要性，這確保了在大旋轉(zhuǎn)、大應變工況下計算的物理準確性。
應力共軛與本構更新
為了保證能量守恒，文章在晶體本地坐標系下采用 Mandel 應力作為滑移驅(qū)動力，并配合隱式時間積分更新塑性變形梯度。
文章的模擬效果如下：
需要注意的是當前的這代積分方案和damask的快速傅里葉變化方案計算效果基本保持一致，整體也是使用fortran語言編寫，并使用vtk格式用于輸出，使用paraview可視化。
使用類似的思想，我們可以根據(jù)文章的公式實現(xiàn)對應的CPFFT的計算方案。這里展示使用matlab實現(xiàn)對應的CPFFT方案，matlab的顯著優(yōu)勢可以很容易和相場和再結晶去結合，因此后續(xù)非常容易擴展。
使用FFT作為邊值問題的求解器，使用固定點迭代完成內(nèi)部的晶體塑性迭代。使用經(jīng)典的位錯密度模型計算硬化和熱激活流動方程計算滑移系的剪切變形。
初始RVE模型使用neper建模，建立一個包含100個晶粒的多晶模型：
matlab導入幾何模型網(wǎng)格：
并沿著X方向進行1.0%的拉伸變形，所有量綱使用m-s-pa。
拉伸變形結束后的累計剪切滑移結果：
拉伸變形結束后的統(tǒng)計儲存位錯密度分布結果：
拉伸變形結束后的幾何必須位錯密度分布結果：
展開 
Damask 2.0.3聯(lián)合abaqus進行晶體塑性有限元模擬的簡單介紹
對于本例中所采用的現(xiàn)象學晶體塑性本構關系，共有114個狀態(tài)變量，分別如下：
SDV1：HomogenizationCount
SDV2：GrainCount
SDV3：CrystalliteCount
SDV4~7：quaternion（取向四元數(shù)）
SDV8~11：deviation from initial orientation as axis (1-3) and angle in degree (4) in crystal reference coordinates（晶體參考坐標系中與初始取向的偏差（以軸（1-3）和角度（4）為單位））
SDV12~20：f，deformation gradient tensor（變形梯度張量）
SDV21~29：fe，elastic deformation gradient tensor（彈性變形梯度張量）
SDV30~38：fp，plastic deformation gradient tensor（塑性變形梯度張量）
SDV39~47：p，first Piola-Kichhoff stress tensor（第一Piola-Kichhoff應力張量）
SDV48~56：s，second Piola-Kichhoff stress tensor（第二Piola-Kichhoff應力張量）
SDV57-65：lp，plastic velocity gradient tensor（塑性速度梯度張量）
SDV66：ConstitutiveCount
SDV67~78：resistance_slip（滑移系阻力）
SDV79~90：shearrate_slip（滑移系剪切應變率）
SDV91~102：resolvedstress_slip（滑移系臨界分切應力）
SDV103~114：accumulatedshear_slip
展開 umat子程序編寫常用的fortran函數(shù)分享（一）
計算2*2矩陣的跡：
subroutine trace2x2(a,aii)
implicit none
real(8), intent(in) :: a(2,2)
real(8), intent(out) :: aii
aii = a(1,1)+a(2,2)
return
end subroutine trace2x2
計算3*3矩陣的跡：
subroutine trace3x3(a,aii)
implicit none
real(8), intent(in) :: a(3,3)
real(8), intent(out) :: aii
aii = a(1,1)+a(2,2)+a(3,3)
return
end subroutine trace3x3
將向量（9*1：變形梯度）轉(zhuǎn)化為3*3矩陣存儲
subroutine vecmat9(dvin,dmout)
implicit none
real(8), intent(in) :: dvin(9)
real(8), intent(out) :: dmout(3,3)
integer :: i
dmout(1,1) = dvin(1)
dmout(1,2) = dvin(2)
dmout(1,3) = dvin(3)
dmout(2,1) = dvin(4)
dmout(2,2) = dvin(5)
dmout(2,3
展開 ABAQUS中實體單元的應用
由于單元變形沒有產(chǎn)生應變能，所以這種彎曲的變形模式是一個零能量模式。由于單元在此模式下沒有剛度，所以不能抵抗此種形式的位移。在粗網(wǎng)格中，這種零能量模式會通過網(wǎng)格擴展出去，從而產(chǎn)生無意義的結果，這就是所謂的沙漏問題。
可在ABAQUS中對減縮積分單元引入少量的人工“沙漏剛度”以限制沙漏模式的擴展。當模型中有更多的單元時，這種剛度在限制沙漏模式方面是更有效的，這意味著只要采用合理的細網(wǎng)格，線性減縮積分單元會給出可接受的結果。對許多應用而言，采用細網(wǎng)格的線性減縮積分單元所產(chǎn)生的誤差是在一個可接受的范圍內(nèi)的。這個結果說明當用這類單元來模擬承受彎曲載荷的結構時，在厚度方向上至少應采用四個單元。當在梁的厚度方向只有一個線性減縮積分單元時，所有的積分點都位于中性軸上，從而該模型將不能抵抗彎曲載荷。(這種情況在表4－2中用*標出)。
因為線性減縮積分單元對變形的魯棒性，因此可在變形很大的模擬中采用剖分較細的此類單元。
二次減縮積分單元也有沙漏模式。然而在正常網(wǎng)格中這種模式幾乎不可能擴展出去，并且在網(wǎng)格足夠細時基本上不會造成什么問題。由于沙漏問題,C3D20R單元的1′6網(wǎng)格計算發(fā)散；若在寬度方向上變?yōu)閮蓚€單元，即2×6網(wǎng)格，就不會發(fā)散，但對于更細的網(wǎng)格，即便在寬度方向上只有一個單元也不會發(fā)散。即使在復雜應狀態(tài)下，二次減縮積分單元對鎖閉并不敏感。因此一般來說，除了大應變的大位移問題和一些接觸分析問題外，這些單元是應力/位移模擬最佳選擇。
4.1.3 非協(xié)調(diào)單元
非協(xié)調(diào)單元是克服完全積分的一階單元的剪力鎖閉問題的一種嘗試。既然剪力鎖閉是由于單元的位移場不能模擬與彎曲相關的運動學而引起的，那么可以考慮把增強單元變形梯度的附加自由度引入到一階單元中去。對變形梯度的加強使一階單元在單元中的變形梯度呈線性變化，如圖4－9(a)所示。
展開 2019.07.29 關于VUMAT的一些參數(shù)
tempOld（nblock）：增量開始時每個材料點的溫度
stretchOld（nblock，ndir+nshr）:增量開始時每個材料點由極坐標分解F=R·U定義的的拉伸張量U
defgradOld（nblock，ndir+2*nshr）：增量開始時每個材料點的變形梯度張量。Stored in 3D as (
F
11
, F22, F33, F12, F23, F31, F21, F32, F13) and in 2D as (F11, F22, F33, F12, F21)
fieldOld（nblock，nfieldv）：增量開始時每個材料點用戶指定的場變量值
stressOld（nblock，ndir+nshr）：增量開始時每個材料點的應力張量
stateOld（nblock，nstatev）：增量開始時每個材料點的狀態(tài)變量值
enerInternOld（nblock）：增量開始時每個材料點的單位質(zhì)量內(nèi)能
enerInelasOld（nblock）：增量開始時每個材料點單位質(zhì)量耗散的非彈性能
tempNew（nblock）：增量結束時每個材料點的溫度
stretchNew（nblock，ndir+nshr）:增量結束時每個材料點由極坐標分解F=R·U定義的的拉伸張量U
defgradNew（nblock，ndir+2*nshr）：增量結束時每個材料點的變形梯度張量
fieldNew（nblock，nfieldv）：增量結束時每個材料點用戶指定的場變量值
展開 