搜索
此为1.0版的修正版，可以考虑混凝土弹性模量随时间，应力变化等情况下的徐变。可以考虑徐变恢复。注意，getvrm中的变量编号与坐标系有关。否则，结果不正确。SUBROUTINEUSDFLD(FIELD,STATEV,PNEWDT,DIRECT,T,CELENT,1TIME,DTIME,CMNAME,ORNAME,NFIELD,NSTATV,NOEL,NPT,LAYER,2KSPT,KSTEP,KINC,NDI,NSHR,COORD,JMAC,JMATYP,MATLAYO,3LACCFLA)CINCLUDE'ABA_PARAM.INC'CCHARACTER*80CMNAME,ORNAMECHARACTER*3FLGRAY(15)DIMENSIONFIELD(NFIELD),STATEV(NSTATV),DIRECT(3,3),1T(3,3),TIME(2)DIMENSIONARRAY(10),JARRAY(10),JMAC(*),JMATYP(*),1COORD(*)CCReadinginstantaneousthermalstrainindirection11(xaxial)CStoringthethermalstraininstatevariableCCALLGETVRM('THE',ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP,1MATLAYO,LACCFLA)STATEV(1)=ARRAY(1)CCReadinginstantaneouselasticxaxialstrainCStoringthethermalstraininstatevariableCNOTE:ARRAY(1)--XAXIALCARRAY(2)--YAXIALCARRAY(3)--ZAXIALCCALLGETVRM('EE',ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP,1MATLAYO,LACCFLA)STATEV(2)=ARRAY(3)CCCModifyingFieldvariablestomodeldevelopmentofYoung’sCmoduluswithage.CcIF(TIME(2).GT.1)THENcFIELD(1)=TIME(2)cELSEcFIELD(1)=0cENDIFRETURNENDCCCSUBROUTINEUEXPAN(EXPAN,DEXPANDT,TEMP,TIME,DTIME,PREDEF,1DPRED,STATEV,CMNAME,NSTATV,NOEL)CINCLUDE'ABA_PARAM.INC'CCHARACTER*80CMNAMECDIMENSIONEXPAN(*),DEXPANDT(*),TEMP(2),TIME(2),PREDEF(*),1DPRED(*),STATEV(NSTATV),ARRAY(15)REALCINI(4,8),T0(8),CLAMDA(4),AN(4)REALCCLOAD(4),EMOD(8)CCSWITCH,II=1,CREEP;II=2,SHRINKAGE,II=3,CREEP+SHRANKAGECKK=3CC=====================INTIALVALUES=========================C*******FORCREEP(USING1STOPT)*********CCINITIALFUNCTIONISY=1.9*T0^(-0.118)*(T-T0)/(61-0.51*FC+T-T0)CCMODELINGFUNCTION(SOFTWAREOF1STOPT)ISCY=C1*(1-EXP(-LAMDA1*(X-T)))+C2*(1-EXP(-LAMDA2*(X-T)))+CC3*(1-EXP(-LAMDA3*(X-T)))+C4*(1-EXP(-LAMDA4*(X-T)));CCWHERETISTHECALCULATINGTIME,T0ISTHELOADINGTIMECTHEFACTORSARELISTEDASFOLLOWINGCLOADINGTIMEARE7,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))THENCELASTICMODULUScEMOD=1.0CVALUESOFLAMDACLAMDA(1)=0.0193907053463775CLAMDA(2)=0.00537298967078406CLAMDA(3)=0.00105487934454054CLAMDA(4)=0.0597818282775493CINITIALVALUESOFCOEFFICIENTC1CINI(1,1)=0.255636119E-4CINI(1,2)=0.179033082E-4CINI(1,3)=0.163880576E-4CINI(1,4)=0.147710753E-4CINI(1,5)=0.124416704E-4CINI(1,6)=0.119344460E-4CINI(1,7)=0.110286669E-4CINI(1,8)=0.096165142E-4CINITIALVALUESOFCOEFFICIENTC2CINI(2,1)=0.157591725E-4CINI(2,2)=0.110425824E-4CINI(2,3)=0.101134277E-4CINI(2,4)=0.091232941E-4CINI(2,5)=0.077200858E-4CINI(2,6)=0.074402105E-4CINI(2,7)=0.068149034E-4CINI(2,8)=0.060733782E-4CINITIALVALUESOFCOEFFICIENTC3CINI(3,1)=0.053080102E-4CINI(3,2)=0.037166275E-4CINI(3,3)=0.034018732E-4CINI(3,4)=0.030634778E-4CINI(3,5)=0.025744519E-4CINI(3,6)=0.024744322E-4CINI(3,7)=0.022776919E-4CINI(3,8)=0.018527641E-4CINITIALVALUESOFCOEFFICIENTC4CINI(4,1)=0.120382947E-4CINI(4,2)=0.084348462E-4CINI(4,3)=0.077257019E-4CINI(4,4)=0.069613683E-4CINI(4,5)=0.058953358E-4CINI(4,6)=0.056904171E-4CINI(4,7)=0.051955726E-4CINI(4,8)=0.045768253E-4CCLOADINGAGE(DAYS)T0=(/7,30,50,100,365,500,1000,3000/)CMODULUSOFELASTICATT0'SAGEEMOD=(/2.546295,3.061405,3.147964,3.217927,*3.271749,3.277386,3.285052,3.290193/)CC========FORSHRINKAGE=======CTGSH--TIMEOFSHRINKAGECALCULATIONFORGIRDER(DAYS)CTDSH--TIMEOFSHRINKAGECALCULATIONFORDECK(DAYS)C480E-6--ULTIMATESHRINKAGESTRAININAASHRO2007(SI)CFC--SPECIFIEDCOMPREI\SSIVESTRENGTHOFCONCRETEATTIMEOFCPRESTRESSINGFORPRESTRESSIONEDMEMBERSANDATTIMEOFCINITIALLOADINGFORNONPRESTRESSIONEDMEMBERS,fci'INAASHTO(2007)CSTANDARDVALUEIS28MPaCCKF--FACTORFORTHEEFFECTOFCONCRETESTRENGTHCCKTD--TIMEDEPENDENT(DEVELOPMENT)FACTOR(THEFIRSTTWOTURMS)CRHSH--AMBIENTRELERTIVEHUMIDITYCORRECIONFACTORFORSHRINKAGETGSH=1.0TDSH=50.0SHU=480E-6RH=70.0FC=28.0CKF=35.0/(7.0+FC)CKTD=61.0-0.58*FCRHSH=2.0-0.014*RHIF(CMNAME.EQ.MAT-SLAB)THENCKF=35.0/(7.0+0.8*FC)CKTD=61.0-0.58*0.8*FCENDIFCC=================INTERPOSITION=====================C****COMPUTINGTHECOEFFICIENTSOFKABIRSERIES****CCSTATEV(1)--THERMALSTRAINCSTATEV(2)--ELASTICSTRAINOFCURRENTINCREMENTCSTATEV(3)--ELASTICSTRAINOFPREVIOUSINCREAMENTCSTATEV(4)--DTIMEOFPREVIOUSINCREAMENTCDELTEE--INCREMENTOFELASICSTRAINTCUR=TIME(2)THE=STATEV(1)EECUR=STATEV(2)EEPRE=STATEV(3)DTPRE=STATEV(4)CIF((TCUR.GE.T0(1)).AND.(TCUR.LT.T0(2)))THENJJ=1TA=T0(JJ)TB=T0(JJ+1)EMA=EMOD(JJ)EMB=EMOD(JJ+1)ELSEIF((TCUR.GE.T0(2)).AND.(TCUR.LT.T0(3)))THENJJ=2TA=T0(JJ)TB=T0(JJ+1)EMA=EMOD(JJ)EMB=EMOD(JJ+1)ELSEIF((TCUR.GE.T0(3)).AND.(TCUR.LT.T0(4)))THENJJ=3TA=T0(JJ)TB=T0(JJ+1)EMA=EMOD(JJ)EMB=EMOD(JJ+1)ELSEIF((TCUR.GE.T0(4)).AND.(TCUR.LT.T0(5)))THENJJ=4TA=T0(JJ)TB=T0(JJ+1)EMA=EMOD(JJ)EMB=EMOD(JJ+1)ELSEIF((TCUR.GE.T0(5)).AND.(TCUR.LT.T0(6)))THENJJ=5TA=T0(JJ)TB=T0(JJ+1)EMA=EMOD(JJ)EMB=EMOD(JJ+1)ELSEIF((TCUR.GE.T0(6)).AND.(TCUR.LT.T0(7)))THENJJ=6TA=T0(JJ)TB=T0(JJ+1)EMA=E