Abaqus 混凝土徐变收缩子程序 2.1.docx

经过进一步的完善后，子程序的计算结果跟一般程序的计算结果吻合得比较好。目前在三维实体元中应用效果良好。注意：外部数据文件的数据由其它拟合软件获得。 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,MATLAYO, 3 LACCFLA)C INCLUDE ABA_PARAM.INCC 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(*), 1 COORD(*)CC Reading instantaneous thermal strain in direction 11(x axial)C??? Storing the thermal strain in state variableC CALL GETVRM(THE,ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP, 1 MATLAYO,LACCFLA) STATEV(1)= ARRAY(1)CC Reading instantaneous elastic x axial strain C??? Storing the thermal strain in state variableC??? NOTE: ARRAY(1)--X AXIALC????????? ARRAY(2)--Y AXIALC????????? ARRAY(3)--Z AXIALC CALL GETVRM(EE,ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP, 1 MATLAYO,LACCFLA) STATEV(2)= ARRAY(1)C CALL GETVRM(SDV,ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP, 1 MATLAYO,LACCFLA) STATEV(3)= ARRAY(3) STATEV(4)= ARRAY(4)CC Modifying Field variables to model development of Young’s C????? modulus with age.C IF(TIME(2).GT.1) THEN FIELD(1)=TIME(2) ELSE FIELD(1)=0 END IF RETURN ENDCC SUBROUTINE UEXTERNALDB(LOP,LRESTART,TIME,DTIME,KSTEP,KINC)C INCLUDE ABA_PARAM.INCC DIMENSION TIME(2) CHARACTER*80 TITLEC COMMON/MYGIRD/GIRD_INI(4,10),GIRD_LAMDA(4) COMMON/MYSLAB/SLAB_INI(4,10),SLAB_LAMDA(4)CC Abaqus uses the FORTRAN unit numbers outlined in the table below. Unless noted otherwise, C????? you should not try to write to these FORTRAN units from user subroutines.C For Abaqus/Standard, you should specify unit numbers 15–18 or unit numbers greater than 100CC Read the creep coefficients from external dababase at the beginning of computing.C IF(LOP.EQ.0) THEN OPEN(106,STATUS=OLD,FILE=D:\ABAQUSWORK\creep.in) READ(106, *)TITLE READ(106, *)(GIRD_LAMDA(I),I=1,4) DO J=1,10 READ(106, *)(GIRD_INI(I,J),I=1,4) END DOC READ(106, *)TITLE READ(106, *)(SLAB_LAMDA(I),I=1,4) DO J=1,