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

Abaqus 混凝土徐变收缩子程序 经过进一步的完善后，子程序的计算结果跟一般程序的计算结果吻合 得比较好。目前在三维实体元中应用效果良好。注意：外部数据文件的数 据由其它拟合软件获得。 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.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(*), 1 COORD(*) C C Reading instantaneous thermal strain in direction 11(x axial) C Storing the thermal strain in state variable C CALL GETVRM('THE',ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP, 1 MATLAYO,LACCFLA) STATEV(1)= ARRAY(1) C C Reading instantaneous elastic x axial strain C Storing the thermal strain in state variable C NOTE: ARRAY(1)--X AXIAL C ARRAY(2)--Y AXIAL C ARRAY(3)--Z AXIAL C 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) C C 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 END C C SUBROUTINE UEXTERNALDB(LOP,LRESTART,TIME,DTIME,KSTEP,KINC) C INCLUDE 'ABA_PARAM.INC' C DIMENSION TIME(2) CHARACTER*80 TITLE C COMMON/MYGIRD/GIRD_INI(4,10),GIRD_LAMDA(4) COMMON/MYSLAB/SLAB_INI(4,10),SLAB_LAMDA(4) C C 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 ForAba