zhanghaidong 发表于 2012-7-2 16:21:13

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

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







zhanghaidong 发表于 2012-7-4 15:01:21

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

by1986188 发表于 2012-7-18 13:40:57

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

zhanghaidong 发表于 2012-7-18 20:46:33

by1986188 发表于 2012-7-18 13:40 static/image/common/back.gif
你好,程序里的硬化程序vuhard为什么会冲突啊,abaqus里不是有vuhard这个子程序吗? ...

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

guidingfa 发表于 2014-1-2 19:29:58

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

zhanghaidong 发表于 2014-3-8 21:42:16

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

解决了吧

xukai439 发表于 2014-3-18 19:09:43

本帖最后由 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

chengziuu 发表于 2014-8-6 17:16:12

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

ffubei 发表于 2014-11-3 14:08:25

感谢楼主分享!初学者。

helj 发表于 2014-12-25 10:58:44

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

你好,我也在编写vumat子程序,子程序里面还调用了vuhard子程序,想用原来的vuhard(系统默认的)子程序,该怎么调用呢?谢谢!

helj 发表于 2014-12-25 10:59:50

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

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

20110503 发表于 2017-4-17 10:41:45

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

liyushu 发表于 2020-8-4 13:14:19

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

pcstudying 发表于 2021-3-21 11:49:58

感谢楼主分享!初学者

zhuhao3802 发表于 2021-9-14 22:49:39


感谢楼主分享!

llzyq 发表于 2021-11-6 13:39:18

感谢楼主分享!
页: [1]
查看完整版本: ABAQUS子程序VUMAT的版本问题(已解决)