虹藻 发表于 2017-10-8 16:04:33

大神们帮忙写一个弹性本构的umat

杨氏模量为28000泊松比为0.2。我的邮箱为1014580620@qq.com，万分感谢。

虹藻 发表于 2017-10-9 10:42:22

已经解决了！

1482524380 发表于 2017-10-25 14:17:51

umat和fortran写的程序是一样的吗？请指教，谢谢

黑明哥丶 发表于 2017-10-28 19:32:16

楼主能把解决方法发我下吗 谢谢啦 807519796@qq.com

划船不用桨 发表于 2018-1-19 20:01:59

楼主能把你的弹性本构的umat发我一下吗？ 1298690928@qq.com
多谢多谢！

hipeilei 发表于 2018-1-19 20:19:12

划船不用桨 发表于 2018-1-19 20:01
楼主能把你的弹性本构的umat发我一下吗？
多谢多谢！

附上代码。
      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,NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT,
   4 CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,JSTEP,KINC)
C
      INCLUDE 'ABA_PARAM.INC'
C
      CHARACTER*80 CMNAME
      DIMENSION STRESS(NTENS),STATEV(NSTATV),
   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),
   4 JSTEP(4)
C   ELASTIC USER SUBROUTINE
      PARAMETER (ONE=1.0D0, TWO=2.0D0)
       E=PROPS(1)
       ANU=PROPS(2)
       ALAMBDA=E/(ONE+ANU)/(ONE-TWO*ANU)
       BLAMBDA=(ONE-ANU)
         CLAMBDA=(ONE-TWO*ANU)
             DO I=1,NTENS
              DO J=1,NTENS
              DDSDDE(I,J)=0.0D0
              ENDDO
             ENDDO
                DDSDDE(1,1)=(ALAMBDA*BLAMBDA)
                DDSDDE(2,2)=(ALAMBDA*BLAMBDA)
                DDSDDE(3,3)=(ALAMBDA*BLAMBDA)
                DDSDDE(4,4)=(ALAMBDA*CLAMBDA)
                DDSDDE(5,5)=(ALAMBDA*CLAMBDA)
                DDSDDE(6,6)=(ALAMBDA*CLAMBDA)
                DDSDDE(1,2)=(ALAMBDA*ANU)
                DDSDDE(1,3)=(ALAMBDA*ANU)
                DDSDDE(2,3)=(ALAMBDA*ANU)
                DDSDDE(2,1)=(ALAMBDA*ANU)
                DDSDDE(3,1)=(ALAMBDA*ANU)
                DDSDDE(3,2)=(ALAMBDA*ANU)
         DO I=1,NTENS
              DO J=1,NTENS
              STRESS(I)=STRESS(I)+DDSDDE(I,J)*DSTRAN(J)
              ENDDO
             ENDDO
      RETURN
      END
材料里面设置成User material，输入280000和0.2两个常数即可。

划船不用桨 发表于 2018-1-22 10:53:59

hipeilei 发表于 2018-1-19 20:19
附上代码。
      SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,
   1 RPL,DDSDDT,DRPLDE,DRPLDT ...

谢谢楼主了

simon1106 发表于 2018-4-13 18:09:24

谢谢大神:)
页: [1]
查看完整版本: 大神们帮忙写一个弹性本构的umat