abaqus材料子程序.docxVIP

各向同性材料损伤本构模型 SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD, + RPL,DDSDDT,DRPLDE,DRPLDT, + STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME, + NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT, + CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,KSTEP,KINC) INCLUDE ABA_PARAM.INC CHARACTER*80 CMNAME DIMENSION STRESS(NTENS),STATEV(NSTATV), + DDSDDE(NTENS,NTENS),DDSDDT(NTENS), + DRPLDE(NTENS),STRAN(NTENS),DSTRAN(NTENS), + TIME(2),PREDEF(1),DPRED(1),PROPS(NPROPS), + COORDS(3),DROT(3,3),DFGRD0(3,3),DFGRD1(3,3) DIMENSION STRANT(6),TSTRANT(4),PT(1) DIMENSION OLD_STRESS(6) DIMENSION DOLD_STRESS(6),D_STRESS(6) DIMENSION C(6,6),CD(6,6),DSTRESS(6),BSTRESS(6),ROOT(3), + DFMNDE(6),DDMDE(6),DCDDM(6,6),ATEMP1(6), ATEMP2(6) PARAMETER (ZERO=0.D0,ONE=1.D0,TWO=2.D0,FOUR=4.D0,HALF = 0.5D0)C startC IF (NPROPS.LT.2) THENC WRITE(7,*) ** ERROR: UMAT REQUIRES *NPROPS=2C STOPC END IF E11 =PROPS(1) V12 =PROPS(2) G12 =PROPS(1)/TWO/(ONE+PROPS(2))C Critical values of stresses XT=PROPS(3) XC=PROPS(4) XS=PROPS(5) GX=PROPS(6) !Fracture energy in matrix ETA=0.001C Current strain DO I = 1, NTENSSTRANT(I) = STRAN(I) + DSTRAN(I) END DOC Stiffness DO I = 1, 6 DO J = 1, 6C(I,J)=ZERO END DO END DO ATEMP = (1+V12)*(1-TWO*V12)C(1,1) = E11*(1-V12)/ATEMPC(2,2) = E11*(1-V12)/ATEMPC(3,3) = E11*(1-V12)/ATEMPC(1,2) = E11*V12/ATEMPC(1,3) = E11*V12/ATEMPC(2,3) = E11*V12/ATEMPC(4,4) = G12C(5,5) = G12C(6,6) = G12 DO I = 2, 6 DO J = 1, I-1C(I,J) = C(J,I) END DO END DOC Critical values of strains XET=XT/(C(1,1)-2*V12*C(1,2)) XEC=XC/(C(1,1)-2*V12*C(1,2)) XES=XS/C(4,4) DMOLD = STATEV(1) C Strain initiation criterion A11 = STRANT(1)**TWO+STRANT(2)**TWO+STRANT(3)**TWO A12 = A11 / XET / XEC A21 = STRANT(1)+STRANT(2)+STRANT(3) A22 = (XEC - XET) / XEC / XE