标签: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. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。