- 积分
- 0
- 注册时间
- 2016-11-10
- 仿真币
-
- 最后登录
- 1970-1-1
|
SUBROUTINE USDFLD(FIELD,STATEV,PNEWDT,DIRECT,T,CELENT,
1 TIME,DTIME,CMNAME,ORNAME,NFIELD,NSTATV,NOEL,NPT,LAYER,
2 KSPT,KSTEP,KINC,NDI,NSHR,COORD,JMAC,JMATYP,
3 MATLAYO,LACCFLA)
C
INCLUDE 'ABA_PARAM.INC'
C
CHARACTER*80 CMNAME,ORNAME
CHARACTER*3 FLGRAY(15)
DIMENSION FIELD(NFIELD),STATEV(NSTATV),DIRECT(3,3),
1 T(3,3),TIME(2)
DIMENSION ARRAY(15),JARRAY(15),JMAC(*),JMATYP(*),COORD(*)
C
REAL*8 A1,A3,Ms,Mf,Q1
Q1=-0.011
A1=723.
A3=898.
Ms=558.
Mf=200.
C
CALL GETVRM('TEMP',ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP,
1 MATLAYO,LACCFLA)
C
TEM=ARRAY(1)
C
IF(KSTEP.EQ.1)THEN
IF(KINC.EQ.0.OR.KINC.EQ.1)THEN
STATEV(1)=0.
STATEV(2)=0.
STATEV(3)=0.
FIELD(1)=0.
END IF
IF(TEM.LT.A1)THEN
STATEV(2)=0.
ELSE IF(TEM.LT.A3)THEN
STATEV(2)=(TEM-A1)/(A3-A1)
ELSE
STATEV(2)=1.
END IF
END IF
C
IF(KSTEP.EQ.2)THEN
IF(TEM.LT.Ms.AND.STATEV(2).NE.0.)THEN
STATEV(2)=1.-EXP(Q1*(Ms-TEM))
END IF
END IF
STATEV(3)=(STATEV(2)-STATEV(1))/DTIME
FIELD(1)=STATEV(1)
STATEV(1)=STATEV(2)
RETURN
END
|
|