深入剖析UMAT线形粘弹性体.DOC

深入剖析UMAT线形粘弹性体 liuxiong simwe用户名：LXCAD [摘要]本文介绍了线形粘弹体的UMAT程序说明，与ZHANG chunyu一起学习UMAT为姐妹篇。 UMAT其实并不难学，你要把握几点即可： 必须提供准确的雅可比距阵，程序收敛速度快 必须用增量法更新应力 must update the stresses and solution-dependent state variables to their values at the end of the increment for which it is called; must provide the material Jacobian matrix 3，与本构方程相关的状态变量必须更新 把我分析的子程序过程多读几遍，相关的弹性力学，朔性力学概念弄懂，可能就理解更为清晰！我的目的就是让大家学的轻松！请大家鼓励！ 三单元体的固体模型(线形粘弹性) 下图为虎克体和开尔文体的串联 一维方向的应力与应变的行为 推广到实体 其中 PROPS(1) PROPS(2) PROPS(3) PROPS(4) PROPS(5) 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), 2 DDSDDT(NTENS),DRPLDE(NTENS), 3 STRAN(NTENS),DSTRAN(NTENS),TIME(2),PREDEF(1),DPRED(1), 4 PROPS(NPROPS),COORDS(3),DROT(3,3),DFGRD0(3,3),DFGRD1(3,3) DIMENSION DSTRES(6),D(3,3) C C EVALUATE NEW STRESS TENSOR C EV = 0. DEV = 0. DO K1=1,NDI EV = EV + STRAN(K1) !! 直接应力和 DEV = DEV + DSTRAN(K1) !!直接增量应力和 END DO C TERM1 = .5*DTIME + PROPS(5) !! TERM1I = 1./TERM1 !! TERM2 = (.5*DTIME*PROPS(1)+PROPS(3))*TERM1I*DEV !! TERM3 = (DTIME*PROPS(2)+2.*PROPS(4))*TERM1I !! C更新正应力 DO K1=1,NDI DSTRES(K1) = TERM2+TERM3*DSTRAN(K1) 1 +DTIME*TERM1I*(PROPS(1)*EV +2.*PROPS(2)*STRAN(K1)-STRESS(K1)) 本构方程 程序执行： ++ STRESS(K1) = STRESS(K1) + DSTRES(K1) END DO C 更新剪应力 TERM2= TERM2 = (.5*DTIME*PROPS(2) + PROPS(4))*TERM1I I1 = NDI DO K1=1,NSHR I1 = I1+1 本构方程 程序执行： *+** DSTRES(I1) = TERM2*DSTRAN(I1)+ 1 DTIME*TERM1I*(PROPS(2)*STRAN(I1)-STRESS(I1)) STRESS(I1) = STRESS(I1)+DSTRES(I1) END DO C雅可比距阵开始 C CREATE NEW JACO