GETVRM
關注創建者：匿名 創建時間：2021-08-13
GETVRM的視頻教程
abaqus中利用UVARM子程序創建自定義場變量
CALL GETVRM('S',ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP, 1 MATLAYO,LACCFLA) JERROR = JERROR + JRCD UVAR(1) = (ARRAY(1)+ARRAY(2)+ARRAY(3))/3.0 UVAR(2) = COORD(1) 完成計算，后處理中會發現輸出變量中會出現
￥1.9
查看
跟兵哥零基礎學ABAQUS的子程序二次開發
內容安排（實時更新）： 1Fortran編譯環境的安裝和關聯；（贈送VS和Fortran安裝包） 2Fortran語法介紹和VS調試入門； 3Fortran語法； 4位移邊界條件子程序DISP； 5載荷邊界條件子程序DLOAD； 6在VDLOAD中實現高斯分布載荷 7工具子程序GETVRM； 8重定義場變量子程序USDFLD； 9VGETVRM介紹顯式隱式區別； 10-11VUDFLD
￥99 11小時32分鐘 107894播放
查看
GETVRM的實例教程
SDEG云圖是有數值的但是調用一直為零麻煩各位大佬幫忙看看 以下是源代碼 SUBROUTINE USDFLD(FIELD,STATEV,PNEWDT,DIRECT,T,CELENT, 1 TIME,DTIME,CMNAME,ORNAME,NFIELD,NSTATV,NOEL,NPT,LAYER, 2 KSPT,KSTEP,KINC,NDI,NSHR,COORD,JMAC,JMATYP,MATLAYO, 3 LACCFLA) C INCLUDE 'ABA_PARAM.INC' C CHARACTER*80 CMNAME,ORNAME CHARACTER*3 FLGRAY(15) DIMENSION FIELD(NFIELD),STATEV(NSTATV),DIRECT(3,3), 1 T(3,3),TIME(2) DIMENSION ARRAY(15),JARRAY(15),JMAC(*),JMATYP(*), 1 COORD(*) C real alfa,K0,D,DAMAGED_K,K PARAMETER(TOLER=1.0E-6) alfa=4.0 K0=1.0E-6 C Absolute value of current strain: CALL GETVRM('SDEG',ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP, 1 MATLAYO,LACCFLA) D=ARRAY(1) IF (D.LE.TOLER) THEN D=0 END IF DAMAGED_K=EXP(alfa*D) K=K0*DAMAGED_K C field variable FIELD(1)=D C state variable STATEV(1)=D STATEV(2)=K C If error, write comment to .DAT file: IF(JRCD.NE
展開 ('TEMP',array,jarray,flgray,jrcd, 1 jmac, jmtyp, matlayo, laccflg)C TEM = ARRAY(1)C C Store temperature as a solution dependent state variableC STATEV(2) = TEM CC Value of current strain :C call getvrm('NEP',array,jarray,flgray,jrcd, 1 jmac, jmtyp, matlayo, laccflg)C DSTRAIN = abs(ARRAY(1))CC Maximum value of strain up to this point in time:C call getvrm('SDV',array,jarray,flgray,jrcd, 1 jmac, jmtyp, matlayo, laccflg)C EPSMAX = ARRAY(1)CC Store maximum value of strain as a solution dependent state variableC STATEV(1) = MAX(DSTRAIN,EPSMAX)C C MINrage value of strainC if(time(2).eq.1) thenCC Store MINrage value of strain as a solution dependent state variableC STATEV(3) = STATEV(1)CC Write StrainMIN to msg fileCC
展開 ('TEMP',array,jarray,flgray,jrcd, 1 jmac, jmtyp, matlayo, laccflg)C TEM = ARRAY(1)C C Store temperature as a solution dependent state variableC STATEV(2) = TEM CC Value of current strain :C call getvrm('NEP',array,jarray,flgray,jrcd, 1 jmac, jmtyp, matlayo, laccflg)C DSTRAIN = abs(ARRAY(1))CC Maximum value of strain up to this point in time:C call getvrm('SDV',array,jarray,flgray,jrcd, 1 jmac, jmtyp, matlayo, laccflg)C EPSMAX = ARRAY(1)CC Store maximum value of strain as a solution dependent state variableC STATEV(1) = MAX(DSTRAIN,EPSMAX)C C MINrage value of strainC if(time(2).eq.1) thenCC Store MINrage value of strain as a solution dependent state variableC STATEV(3) = STATEV(1)CC Write StrainMIN to msg fileCC
展開 比如，在USDFLD場變量子程序或者UMAT自定義材料子程序開發過程中，當需要獲取積分點上的應力、應變等數據時，就可以通過“Utility Routines”中的GETVRM程序實現，如下：
CALL GETVRM('VAR',ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP,MATLAYO,
LACCFLA)
02
案例參考指南
案例參考指南中包含了大量Abaqus分析 案例，其中更多的案例是跟子程序無關的，要想找到自己想要學習的案例，通過關鍵字進行搜索即可，如上圖所示，想要學習UMAT子程序開發，只需要搜索"UMAT"關鍵字即可。可以發現在靜力分析部分有57條與UMAT這個關鍵字相關的信息，可以逐個打開，也可以打開對應關鍵字信息最多的章節查看。
比如打開損傷與斷裂一章，其中有一節是纖維-金屬復合材料（GLARE層板）開孔板失效分析的案例，記住一句話，案例參考指南中必定有模型文件和源代碼！運氣好的時候連參數化python建模腳本都有。
那么，打開這一節內容，看一看案例指南的套路：
(1)首先是交代當前案例問題的背景，并做內容概述
(2)模型信息介紹，包括幾何、材料、分析類別、邊界載荷等
(3)理論詳解：因為這個案例是要自定義一種復合材料的損傷本構模型，所以要把用到的本構關系講解清楚、包含了損傷起始判據、損傷演化、剛度折減等，想學復合材料UMAT子程序開發的同學看到這里眼前是不是一亮？
剛度折減和雅可比矩陣都給你準備好了。
(4)模型介紹完了以后會展示一些計算結果。
展開 材料實際應力可表示為
因此，需要在代碼中通過GETVRM調用相應的應力分量S11、S22、S33、S12、S13和S23。
除此之外，還需要在前處理中設置：
（1） NUVARM為材料屬性中定義的Number of user-defined variables at each material point數量，設置界面如下圖1所示。如NUVARM=3，則在子程序中需要定義UVAR(1)、UVAR(2)和UVAR(3)。
圖1
（2） 為能夠在后處理中顯示該自定義變量，還需要在Step模塊中場變量輸出請求中勾選UVARM，如下圖2所示。
圖2
下面在實例中應用UVARM子程序以在后處理中顯示安全系數云圖。下圖3是由軸承座、軸承和軸組成的裝配體，所用的材料有QT600（屈服強度850MPa）和45鋼（355MPa）。在Properties中設置UVARM的數量為3個；在Step場變量輸出請求中勾選UVARM。在軸兩端分別施加向下的10000N載荷。
圖3
如下圖4所示的是UARM1云圖，其對應的是45鋼的安全系數云圖，即屈服強度（355）與材料實際應力之比。軸承與軸區域的安全系數顯示是正確的。然而，材質為QT600的軸承座安全系數并不正確。
圖4
如下圖5所示的是UARM2云圖，其對應的是QT600的安全系數云圖，即屈服強度（850）與材料實際應力之比。軸承座區域的安全系數顯示的是正確的。然而，材質為45鋼的軸承和軸安全系數并不正確。
圖5
從上面UARM1和UVARM2云圖中可以看出，兩種材質在不同的安全系數云圖中顯示數值不一樣。因此，需要單獨顯示相應材質的結構件才能獲得正確的安全系數值。
展開 
GETVRM的最新內容
('NEP',array,jarray,flgray,jrcd, 1 jmac, jmtyp, matlayo, laccflg)C DSTRAIN = abs(ARRAY(1))CC Maximum value of strain up to this point in time:C call getvrm('SDV',array,jarray,flgray,jrcd
('NEP',array,jarray,flgray,jrcd, 1 jmac, jmtyp, matlayo, laccflg)C DSTRAIN = abs(ARRAY(1))CC Maximum value of strain up to this point in time:C call getvrm('SDV',array,jarray,flgray,jrcd
DIRECT(3,3), 1 T(3,3),TIME(2) DIMENSION ARRAY(15),JARRAY(15),JMAC(*),JMATYP(*), 1 COORD(*) C real alfa,K0,D,DAMAGED_K,K PARAMETER(TOLER=1.0E-6) alfa=4.0 K0=1.0E-6 C Absolute value of current strain: CALL GETVRM
③USDFLD可以訪問計算結果數據（利用GETVRM子程序），即材料屬性可以定義為結果數據的函數
2.可以采用表格輸入方法或其他用戶子程序,進行相關函數fi的定義
3.USDFLD可以與CREEP、HETVAL、UHYPEL、UEXPANUHARD、UMAT、UMATHT、UTRS、UINTER
材料實際應力可表示為
因此，需要在代碼中通過GETVRM調用相應的應力分量S11、S22、S33、S12、S13和S23。
除此之外，還需要在前處理中設置：
（1） NUVARM為材料屬性中定義的Number of user-defined variables at each material point數量，設置界面如下圖1所示。
比如，在USDFLD場變量子程序或者UMAT自定義材料子程序開發過程中，當需要獲取積分點上的應力、應變等數據時，就可以通過“Utility Routines”中的GETVRM程序實現，如下：
CALL GETVRM('VAR',ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP,MATLAYO,
LACCFLA)
02
Q: 在abaqus中，如果采用umat,利用自己的本構，如何讓abaqus明白這種材料的彈塑性應變，也就是說，如何讓程序返回彈性應變與塑性應變，好在output中輸出，我曾想用最笨地方法，在uvarm中定義輸出，利用getvrm獲取材料點的值，但無法獲取增量應力，材料常數等，研究了幫助中的例子，umatmst3.inp,umatmst3.for,他采用mises J2 流動理論，我在output
Q: 在abaqus中，如果采用umat,利用自己的本構，如何讓abaqus明白這種材料的彈塑性應
變，也就是說，如何讓程序返回彈性應變與塑性應變，好在output中輸出，我曾想用最笨
地方法，在uvarm中定義輸出，利用getvrm獲取材料點的值，但無法獲取增量應力，材料常
數等，研究了幫助中的例子，umatmst3.inp,umatmst3.for,他采用mises J2 流動理論，我
