本人编写了一个ABAQUS各向同性随动强化子程序，跑了以下发生报错，求指点
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
CHARACTER80 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),STRESS_A(NTENS),STRESS_D(NTENS),
5 STRESS_TR(NTENS),STRESS_TR_D(NTENS),STRESS_TR_A(NTENS),
6 DSTRESS_B(NTENS),DSTRAN_P(NTENS),DSTRAN_E(NTENS),
7 S_S(NTENS,NTENS)
DOUBLE PRECISION EE,NEU,SGM_Y,alfa,beta,RR,lamda,GG
DOUBLE PRECISION STRESS_V,AA,CC,DP
DOUBLE PRECISION STRESS_TR_V,STRESS_TR_MISES
Integer i,j
C***********************************************************************
CINPUT_ELASTIC_PROPERTIES
C****************************************************************
EE=PROPS(1) !young's modulus
NEU=PROPS(2) !Possion_radio
SGM_Y=PROPS(3) !yield_stress
alfa=PROPS(4) !material_parameter
beta=PROPS(5) !material_parameter
RR=PROPS(6) !material_parameter
lamda=EENEU/((1.0+NEU)(1.0-2.0NEU)) !Lame_constant
GG=EE/(2.0(1.0+NEU)) !shear_modulus
C************************************************************************
CCalculate_the_elastic_matrix_of_stifness
C********************************************************************
DDSDDE=0.0
DO i=1,NDI
DO j=1,NDI
DDSDDE(i,j)=lamda
END DO
DDSDDE(i,i)=lamda+2.0*GG
DDSDDE(i+NDI,i+NDI)=GG
END DO
C************************************************************************
CCalculate_the_trying_stress
C*************************************************************
DO i=1,NTENS
STRESS_TR(i)=STRESS(i)
END DO
DO i=1,NTENS
DO j=1,NTENS
STRESS_TR(i)=STRESS_TR(i)+DDSDDE(i,j)*DSTRAN(j)
END DO
END DO
C************************************************************************
Ccalculate_the_stress_try_mises
C********************************************************
CSTATEV(1)~STATEV(NTENS)_为rij背应力*************
STRESS_TR_V=0.0
DO i=1,NDI
STRESS_TR_V=STRESS_TR_V+STRESS_TR(i)
END DO
STRESS_TR_V=STRESS_TR_V/3.0
DO i=1,NDI
STRESS_TR_D(i)=STRESS_TR(i)-STRESS_TR_V !STRESS_D偏应力张量
END DO
DO i=1+NDI,NTENS
STRESS_TR_D(i)=STRESS_TR(i)
END DO
DO i=1,NTENS
STRESS_TR_A(i)=STRESS_TR_D(i)-STATEV(i) !STRESS_A_=STRESS_D-背应力
END DO
STRESS_TR_MISES=0.0
DO i=1,NDI
STRESS_TR_MISES=STRESS_TR_MISES+STRESS_TR_A(i)2.0
END DO
DO i=1+NDI,NTENS
STRESS_TR_MISES=STRESS_TR_MISES+2.0STRESS_TR_A(i)*2.0
END DO
STRESS_TR_MISES=sqrt(1.5STRESS_TR_MISES)
P=STATEV(NTENS+1) !hardning_parameter_Odkvista 有效塑性应变
IF (STRESS_TR_MISES. LE. SGM_Y) THEN
DO i=1,NTENS
STRESS(i)=STRESS_TR(i)
END DO
DP=0.0
ELSE
AA=1.0+RR/(3.0*GG)
DO i=1,NTENS
AA=AA-alfaSTATEV(i)STRESS_TR_A(i)/(2.0GGSTRESS_TR_MISES
1 (1.0+betaP))
END DO
DP=0.0
DO i=1,NTENS
DP=DP+STRESS_TR_A(i)DSTRAN(i)/(AASTRESS_TR_MISES) !计算有效塑性应变增量
END DO
DO i=1,NTENS
DSTRAN_P(i)=1.5DPSTRESS_TR_A(i)/STRESS_TR_MISES !计算塑性应变增量
END DO
DO i=1,NTENS
DSTRAN_E(i)=DSTRAN(i)-DSTRAN_P(i)
END DO
DO i=1,NTENS
DO j=1, NTENS
STRESS(i)=STRESS(i)+DDSDDE(i,j)DSTRAN_E(j) !更新应力
END DO
END DO
DO i=1,NTENS
DSTRESS_B(i)=2.0RRDSTRAN_P(i)/3.0-alfaSTATEV(i)DP !计算背应力增量
1 /(1+betaP)
END DO
DO i=1,NTENS
STATEV(i)=STATEV(i)+DSTRESS_B(i) !更新背应力
END DO
C*********************************************************************
C更新应力**********
C************************************************************************
CC=2.0AASTRESS_TR_MISES2.0/3.0
C**********************************************************************
C更新雅克比矩阵****
C************************************************************************
DO i=1,NTENS
DO j=1,NDI
S_S(i,j)=2.0GGSTRESS_TR_A(i)STRESS_A(j)/CC
END DO
DO j=1+NDI,NTENS
S_S(i,j)=GGSTRESS_TR_A(i)*STRESS_A(j)/CC
END DO
END DO
DO i=1,NTENS
DO j=1,NTENS
DDSDDE(i,j)=DDSDDE(i,j)-S_S(i,j)
END DO
END DO
DO i=1,NTENS
DO j=1,NTENS
STRESS(i)=STRESS(i)+DDSDDE(i,j)*DSTRAN(j)
END DO
END DO
END IF
STATEV(NTENS+1)=P+DP !statev7保存有效塑性应变
RETURN
END
报错内容为libirc.lib(fast_mem_ops.obj) : warning LNK4210: 存在 .CRT 节；可能有未处理的静态初始值设定项或结束符
End Linking Abaqus/Standard User Subroutines
你好，我是有问必答小助手，非常抱歉，本次您提出的有问必答问题，技术专家团超时未为您做出解答
本次提问扣除的有问必答次数，将会以问答VIP体验卡（1次有问必答机会、商城购买实体图书享受95折优惠）的形式为您补发到账户。
因为有问必答VIP体验卡有效期仅有1天，您在需要使用的时候【私信】联系我，我会为您补发。