elastic.for--编法1

2022-07-09 00:01:55 阅读：179 来源： 互联网

标签：DO DDSDDE elastic -- TWO NTENS MU 编法 LAMBDA

      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,KSTEP,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)

C
      PARAMETER(ZERO=0.D0,ONE=1.D0,TWO=2.D0)
      DOUBLE PRECISION E,NU,LAMBDA,MU
C
C ELASTIC PROPERTIES
C
      E  = PROPS(1)      ! 弹性模量
      NU = PROPS(2)      ! 泊松比

      LAMBDA = E*NU/((ONE + NU)*(ONE - TWO*NU))
      MU = E/(TWO*(ONE + NU))

C     Jacobi矩阵初始化
      DO K1=1,NTENS
        DO K2=1,NTENS
           DDSDDE(K1,K2)=ZERO
        END DO
      END DO            

      DDSDDE(1,1) = LAMBDA + TWO*MU
      DDSDDE(2,2) = LAMBDA + TWO*MU
      DDSDDE(3,3) = LAMBDA + TWO*MU
      DDSDDE(4,4) = MU
      DDSDDE(5,5) = MU
      DDSDDE(6,6) = MU
      DDSDDE(1,2) = LAMBDA
      DDSDDE(1,3) = LAMBDA
      DDSDDE(2,3) = LAMBDA
      DDSDDE(2,1) = LAMBDA
      DDSDDE(3,1) = LAMBDA
      DDSDDE(3,2) = LAMBDA

C     CALCULATE STRESS
      DO K1=1,NTENS
         DO K2=1,NTENS
            STRESS(K1) = STRESS(K1) + DDSDDE(K1,K2)*DSTRAN(K2)
         END DO
      END DO

      RETURN
      END

标签：DO,DDSDDE,elastic,--,TWO,NTENS,MU,编法,LAMBDA
来源： https://www.cnblogs.com/Nazorine/p/16459946.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

ICode9

elastic.for--编法1