ABAQUS用户子程序`UMAT`在线弹性问题模拟中的应用

阅读数 849

需要将abaqus与vs和fortran进行关联，使用的是abaqus2022版本，vs2019和intel ONEAPI 2022版本。参考关联链接保姆级Abaqus2022+VS2019+Fortran关联配置 - 知乎 (zhihu.com)

创建线弹性模型简单的几何模型

问题: 悬臂梁（长10，宽1，高1）一端四个点固定约束，另外一个端两个点载荷-1000。各项同性材料，E=1e9，v=0.3。

模型

使用UMAT进行分析

与abaqus常规模型分析步骤一样，下面不一样的分析。

定义材料 General->User Material 使用默认的参数添加弹性模量1e9和泊松比0.3.

定义材料

将单元默认的减缩积分C3D8R改成C3D8单元（如果使用C3D8R进行计算的时候会出现报错。）

分配单元类型

创建job，指定UMAT的fortran子程序

创建job

查看输出和结果

test_umat

对比常规输入的材料计算结果，与`UMAT`定义的材料一样。

常规计算结果

`UMAT` 子程序代码

*USER SUBROUTINES
      SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,
     1 RPL,DDSDDT,DRPLDE,DRPLDT,
     2 STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME,
     3 NDI,NSHR,NTENS,NSTATEV,PROPS,NPROPS,COORDS,DROT,PNEWDT,
     4 CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,KSTEP,KINC)
C
    INCLUDE 'ABA_PARAM.INC'
C
      CHARACTER*8 CMNAME
      DIMENSION STRESS(NTENS),STATEV(NSTATEV),
     1 DDSDDE(NTENS,NTENS),DDSDDT(NTENS),DRPLDE(NTENS),
     2 STRAN(NTENS),DSTRAN(NTENS),TIME(2),PREDEF(1),DPRED(1),
     3 PROPS(NPROPS),COORDS(3),DROT(3,3),DFGRD0(3,3),DFGRD1(3,3)
C
      PARAMETER (ONE=1.0D0, TWO=2.0D0)
      write(6,*) PROPS(1) , PROPS(2)
      E=PROPS(1)
      ANU=PROPS(2)
      ALAMDA=ANU*E/ (ONE+ANU)/(ONE-TWO*ANU)
      AMU=E/TWO/(ONE+ANU)
      DO I=1,NTENS
       DO J=1,NTENS
        DDSDDE(I,J)=0.0D0
       ENDDO
      ENDDO
      DDSDDE(1,1)=ALAMDA+TWO*AMU
      DDSDDE(2,2)=DDSDDE(1,1)
      DDSDDE(3,3)=DDSDDE(1,1)
      DDSDDE(4,4)=AMU
      DDSDDE(5,5)=AMU
      DDSDDE(6,6)=AMU
      DDSDDE(1,2)=ALAMDA
      DDSDDE(1,3)=ALAMDA
      DDSDDE(2,3)=ALAMDA      
      DDSDDE(2,1)=DDSDDE(1,2)
      DDSDDE(3,1)=DDSDDE(1,3)
      DDSDDE(3,2)=DDSDDE(2,3)
C      
      DO I=1,NTENS
        DO J=1,NTENS
        STRESS(I)=STRESS(I)+DDSDDE(I,J)*DSTRAN(J)
        ENDDO
      ENDDO
      RETURN
      END