ABAQUS子程序VUMAT的版本问题(已解决)

zhanghaidong

版本是: 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,依然出现上述错误.

恳请各位大侠帮忙分析下.

zhanghaidong

问题已解决，不会是版本问题。
子程序中有一个硬化函数vuhard（……）与abaqus自带的冲突，变成vuhardk（……）就好了。另附vumat和inp

by1986188

你好，程序里的硬化程序vuhard为什么会冲突啊，abaqus里不是有vuhard这个子程序吗？

zhanghaidong

by1986188 发表于 2012-7-18 13:40
你好，程序里的硬化程序vuhard为什么会冲突啊，abaqus里不是有vuhard这个子程序吗？ ...

你可以自己试试，自定义的是要避开原来的

guidingfa

楼主我想问一下，子程序算出来所有的结果都为0是怎么回事？

zhanghaidong

guidingfa 发表于 2014-1-2 19:29
楼主我想问一下，子程序算出来所有的结果都为0是怎么回事？

解决了吧

xukai439

我也出现了同样的问题，我是一个初学者，现在只是算一个简单的模型的Mises把整个流程走通，网上找了个简单的程序，调用了出现了同样的问题，怎么回事[code]      
    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)
C  UMAT FOR ISOTROPIC ELASTICITY
C  CANNOT BE USED FOR PLANE STRESS
C ----------------------------------------------------------------
C  PROPS(1) - E
C  PROPS(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  

chengziuu

请问楼主，子程序算出来的结果为零是怎么回事？

ffubei

感谢楼主分享！初学者。

helj

zhanghaidong 发表于 2012-7-18 20:46
你可以自己试试，自定义的是要避开原来的

你好，我也在编写vumat子程序，子程序里面还调用了vuhard子程序，想用原来的vuhard（系统默认的）子程序，该怎么调用呢？谢谢！

helj

zhanghaidong 发表于 2012-7-18 20:46
你可以自己试试，自定义的是要避开原来的

如果想用原来的vuhard子程序，该怎么办呢？不用自己定义的

20110503

新手一个，想问一下楼主，用VUMAT编写JC本构模型需要从哪入手？

liyushu

请问怎么知道是哪个函数很Abaqus函数重了呢？

pcstudying

感谢楼主分享！初学者

zhuhao3802

感谢楼主分享！

llzyq

感谢楼主分享！