abaqus溷凝土徐变计算.doc

abaqus混凝土徐变计算2.0此为1.0版的修正版，可以考虑混凝土弹性模量随时间，应力变化等情况下的徐变可以考虑徐变恢复注意，getvrm中的变量编号与坐标系有关否则，结果不正确        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(10),JARRAY(10),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(3)C        CC    Modifying Field variables to model development of Young’s C      modulus with age.Cc        IF(TIME(2).GT.1) THENc           FIELD(1)=TIME(2)c        ELSEc           FIELD(1)=0c        END IF        RETURN        ENDCCC        SUBROUTINE UEXPAN(EXPAN,DEXPANDT,TEMP,TIME,DTIME,PREDEF,     1 DPRED,STATEV,CMNAME,NSTATV,NOEL)C        INCLUDE 'ABA_PARAM.INC'C        CHARACTER*80 CMNAMEC        DIMENSION EXPAN(*),DEXPANDT(*),TEMP(2),TIME(2),PREDEF(*),     1 DPRED(*),STATEV(NSTATV),ARRAY(15)        REAL CINI(4,8),T0(8),CLAMDA(4),AN(4)        REAL CCLOAD(4),EMOD(8)CC    SWITCH, II=1, CREEP; II=2, SHRINKAGE, II=3, CREEP+SHRANKAGEC        KK=3 C            C===================== INTIAL VALUES =========================C ******* FOR CREEP (USING 1STOPT) *********C       C    INITIAL FUNCTION IS Y=1.9*T0^(-0.118)*(T-T0)/(61-0.51*FC+T-T0)CC   MODELING FUNCTION (SOFTWARE OF 1STOPT) ISC     Y=C1*(1-EXP(-LAMDA1*(X-T)))+C2*(1-EXP(-LAMDA2*(X-T)))+C       C3*(1-EXP(-LAMDA3*(X-T)))+C4*(1-EXP(-LAMDA4*(X-T)));   CC     WHERE T IS THE CALCULATING TIME, T0 IS THE LOADING TIMEC     THE FACTORS ARE LISTED AS FOLLOWINGC     LOADING TIME ARE 7,30,50,100,365,1000,3000(DAYS)C   ---------------------------------------------------------------C   | A |   7 | 30 | 50 | 100 | 365 | 500 | 1000 | 3000 |C   | C1 |C(1,1)|C(1,2)|C(1,3)|C(1,4)|C(1,5)|C(1,6)|C(1,7)|C(1,8)|C   | C2 |C(2,1)|C(2,2)|C(2,3)|C(2,4)|C(2,5)|C(2,6)|C(2,7)|C(2,8)|C   | C3 |C(3,1)|C(3,2)|C(3,3)|C(3,4)|C(3,5)|C(3,6)|C(3,7)|C(3,8)|C   | C4 |C(4,1)|C(4,2)|C(4,3)|C(4,4)|C(4,5)|C(4,6)|C(4,7)|C(4,8)|      C   ---------------------------------------------------------------        IF((CMNAME.EQ."MAT-GIRDER").OR.(CMNAME.EQ."MAT-SLAB"))THENC   ELASTIC MODULUSc          EMOD=1.0 C   INITIAL VALUES OF COEFFICIENT C1        CINI(1,1)=0.255636119E-4        CINI(1,2)=0.179033082E-4        CINI(1,3)=0.163880576E-4        CINI(1,4)=0.147710753E-4        CINI(1,5)=0.124416704E-4        CINI(1,6)=0.119344460E-4        CINI(1,7)=0.110286669E-4          CINI(1,8)=0.096165142E-4 C   INITIAL VALUES OF COEFFICIENT C2       CINI(2,1)=0.157591725E-4        CINI(2,2)=0.110425824E-4        CINI(2,3)=0.101134277E-4        CINI(2,4)=0.091232941E-4        CINI(2,5)=0.077200858E-4        CINI(2,6)=0.074402105E-4        CINI(2,7)=0.068149034E-4        CINI(2,8)=0.060733782E-4 C   INITIAL VALUES OF COEFFICIENT C3          CINI(3,1)=0.053080102E-4           CINI(3,2)=0.037166275E-4           CINI(3,3)=0.034018732E-4           CINI(3,4)=0.030634778E-4           CINI(3,5)=0.025744519E-4           CINI(3,6)=0.024744322E-4           CINI(3,7)=0.022776919E-4           CINI(3,8)=0.018527641E-4C   INITIAL VALUES OF COEFFICIENT C4          CINI(4,1)=0.120382947E-4           CINI(4,2)=0.084348462E-4           CINI(4,3)=0.077257019E-4           CINI(4,4)=0.069613683E-4          CINI(4,5)=0.058953358E-4           CINI(4,6)=0.056904171E-4           CINI(4,7)=0.051955726E-4           CINI(4,8)=0.045768253E-4     C              C   LOADING AGE(DAYS)          T0=(/7,30,50,100,365,500,1000,3000/)C   MODULUS OF ELASTIC AT T0'S AGE                    EMOD=(/2.546295,3.061405,3.147964,3.217927,     *     3.271749,3.277386,3.285052,3.290193/)C        C    ======== FOR SHRINKAGE =======C    TGSH-- TIME OF SHRINKAGE CALCULATION FOR "GIRDER" (DAYS)C    TDSH-- TIME OF SHRINKAGE CALCULATION FOR "DECK" (DAYS)C    480E-6-- ULTIMATE SHRINKAGE STRAIN IN AASHRO 2007(SI)C    FC --SPECIFIED COMPREI\SSIVE STRENGTH OF CONCRETE AT TIME OF C          PRESTRESSING FOR PRESTRESSIONED MEMBERS AND AT TIME OF C          INITIAL LOADING FOR NONPRESTRESSIONED MEMBERS, fci'IN AASHTO(2007)C          STANDARD VALUE IS 28MPaC。