壳单元与实体单元 UMAT 有哪些差异？

浏览: 2248

我在 Abaqus 中使用壳单元（如 S4R）调用 UMAT 时发现结果与实体单元不同，想确认具体差异。

下面是一个简单的线弹性实体单元umat，如果我想改为用壳单元使用，应该如何修改？

非常感谢大佬解答，解决了可以支付宝发您红包，

我的微信：13635374241


      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     线弹性例子
C     计算杨氏模量E泊松比XNU剪切模量G拉梅常数Lamda
      E = PROPS(1)
      Xmu = PROPS(2)
      G = E/2/(1+Xmu)
      Lamda = Xmu*E/(1+Xmu)/(1-2*Xmu)
C
C     雅克比矩阵初值为0
      DDSDDE(i,j) = 0
C
C     填入雅克比矩阵左上角为lamda
      DO i=1,NDI
          DO j=1,NDI
              DDSDDE(i,j)=Lamda
          END DO
      END DO
C
C     左上角主对角元素改为Lamda+2G
      DO i=1,NDI
          DDSDDE(i,i)=Lamda+2*G
      END DO
C
C     填入雅克比矩阵右下角
      DO i=NDI+1,NTENS
          DDSDDE(i,i)=G
      END DO
C
C     由应变增量计算应力
      DO i=1,NTENS
          DO j=1,NTENS
              STRESS(i)=STRESS(i)+DDSDDE(i,j)*DSTRAN(j)
          END DO
      END DO

      RETURN
      END



全部回答 (0)

默认 最新