ABAQUS子程序VUMAT的版本问题(已解决)
本帖最后由 zhanghaidong 于 2012-7-4 20:47 编辑版本是: Abaqus6.11-1
采用vumat(Johnson-Cook)模拟长方体受波冲击加载,出现错误,提示是版本的问题.
Error in job cell1: Problem during linking - Single Precision Abaqus/Explicit User Subroutines. This error may be due to a mismatch in the Abaqus user subroutine arguments. These arguments sometimes change from release to release, so user subroutines used with a previous release of Abaqus may need to be adjusted.
Job cell1 aborted due to errors.
改用例题的各向同性硬化VUMAT,依然出现上述错误.
恳请各位大侠帮忙分析下.
问题已解决,不会是版本问题。
子程序中有一个硬化函数vuhard(……)与abaqus自带的冲突,变成vuhardk(……)就好了。另附vumat和inp 你好,程序里的硬化程序vuhard为什么会冲突啊,abaqus里不是有vuhard这个子程序吗? by1986188 发表于 2012-7-18 13:40 static/image/common/back.gif
你好,程序里的硬化程序vuhard为什么会冲突啊,abaqus里不是有vuhard这个子程序吗? ...
你可以自己试试,自定义的是要避开原来的 楼主我想问一下,子程序算出来所有的结果都为0是怎么回事? guidingfa 发表于 2014-1-2 19:29
楼主我想问一下,子程序算出来所有的结果都为0是怎么回事?
解决了吧 本帖最后由 xukai439 于 2014-3-18 19:11 编辑
我也出现了同样的问题,我是一个初学者,现在只是算一个简单的模型的Mises把整个流程走通,网上找了个简单的程序,调用了出现了同样的问题,怎么回事
SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,RPL,DDSDDT,
1 DRPLDE,DRPLDT,STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,
2 CMNAME,NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,
3 PNEWDT,CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,KSTEP,KINC)
INCLUDE 'ABA_PARAM.INC'
CHARACTER*8 CMNAME
DIMENSION STRESS(NTENS),STATEV(NSTATV),DDSDDE(NTENS,NTENS),
1 DDSDDT(NTENS),DRPLDE(NTENS),STRAN(NTENS),DSTRAN(NTENS),
2 TIME(2),PREDEF(1),DPRED(1),PROPS(NPROPS),COORDS(3),DROT(3,3),
3 DFGRD0(3,3),DFGRD1(3,3)
CUMAT FOR ISOTROPIC ELASTICITY
CCANNOT BE USED FOR PLANE STRESS
C ----------------------------------------------------------------
CPROPS(1) - E
CPROPS(2) - NU
C ----------------------------------------------------------------
C
IF (NDI.NE.3) THEN
WRITE(6,*)'THIS UMAT MAY ONLY BE USED FOR ELEMENTS
1 WITH THREE DIRECT STRESS COMPONENTS'
CALL XIT
ENDIF
C
C ELASTIC PROPERTIES
EMOD=PROPS(1)
ENU=PROPS(2)
EBULK3=EMOD/(1-2*ENU)
EG2=EMOD/(1+ENU)
EG=EG2/2
EG3=3*EG
ELAM=(EBULK3-EG2)/3
C
C ELASTIC STIFFNESS
C
DO K1=1,NDI
DO K2=1,NDI
DDSDDE(K2, K1)=ELAM
END DO
DDSDDE(K1,K1)=EG2+ELAM
END DO
DO K1=NDI+1,NTENS
DDSDDE(K1 ,K1)=EG
END DO
C
C CALCULATE STRESS
C
DO K1=1,NTENS
DO K2=1,NTENS
STRESS(K2)=STRESS(K2)+DDSDDE(K2, K1)*DSTRAN(K1)
END DO
END DO
C
RETURN
END 请问楼主,子程序算出来的结果为零是怎么回事? 感谢楼主分享!初学者。 zhanghaidong 发表于 2012-7-18 20:46
你可以自己试试,自定义的是要避开原来的
你好,我也在编写vumat子程序,子程序里面还调用了vuhard子程序,想用原来的vuhard(系统默认的)子程序,该怎么调用呢?谢谢! zhanghaidong 发表于 2012-7-18 20:46
你可以自己试试,自定义的是要避开原来的
如果想用原来的vuhard子程序,该怎么办呢?不用自己定义的 新手一个,想问一下楼主,用VUMAT编写JC本构模型需要从哪入手? 请问怎么知道是哪个函数很Abaqus函数重了呢? 感谢楼主分享!初学者
感谢楼主分享! 感谢楼主分享!
页:
[1]