搜索
经过进一步的完善后，子程序的计算结果跟一般程序的计算结果吻合得比较好。目前在三维实体元中应用效果良好。注意：外部数据文件的数据由其它拟合软件获得。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(15),JARRAY(15),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(1)CCALLGETVRM('SDV',ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP,1MATLAYO,LACCFLA)STATEV(3)=ARRAY(3)STATEV(4)=ARRAY(4)CCModifyingFieldvariablestomodeldevelopmentofYoung’sCmoduluswithage.CIF(TIME(2).GT.1)THENFIELD(1)=TIME(2)ELSEFIELD(1)=0ENDIFRETURNENDCCSUBROUTINEUEXTERNALDB(LOP,LRESTART,TIME,DTIME,KSTEP,KINC)CINCLUDE'ABA_PARAM.INC'CDIMENSIONTIME(2)CHARACTER*80TITLECCOMMON/MYGIRD/GIRD_INI(4,10),GIRD_LAMDA(4)COMMON/MYSLAB/SLAB_INI(4,10),SLAB_LAMDA(4)CCAbaqususestheFORTRANunitnumbersoutlinedinthetablebelow.Unlessnotedotherwise,CyoushouldnottrytowritetotheseFORTRANunitsfromusersubroutines.CForAbaqus/Standard,youshouldspecifyunitnumbers15–18orunitnumbersgreaterthan100CCReadthecreepcoefficientsfromexternaldababaseatthebeginningofcomputing.CIF(LOP.EQ.0)THENOPEN(106,STATUS='OLD',FILE='D:\ABAQUSWORK\creep.in')READ(106,*)TITLEREAD(106,*)(GIRD_LAMDA(I),I=1,4)DOJ=1,10READ(106,*)(GIRD_INI(I,J),I=1,4)ENDDOCREAD(106,*)TITLEREAD(106,*)(SLAB_LAMDA(I),I=1,4)DOJ=1,10READ(106,*)(SLAB_INI(I,J),I=1,4)ENDDOCLOSE(106)CcOPEN(107,FILE='D:\ABAQUSWORK\CREEP.OUP')cWrite(107,100)(SLAB_LAMDA(I),I=1,4)cDOJ=1,10cwrite(107,100)(SLAB_INI(I,J),I=1,4)cENDDOcCLOSE(107)c100FORMAT(1X,4E16.7)ENDIFRETURNENDCSUBROUTINEUEXPAN(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,10),T0(10),CLAMDA(4),AN(4)REALCLOAD(4),EMOD(10)cCOMMON/MYGIRD/GIRD_INI(4,10),GIRD_LAMDA(4)COMMON/MYSLAB/SLAB_INI(4,10),SLAB_LAMDA(4)CCSWITCH,II=1,CREEP;II=2,SHRINKAGE,II=3,CREEP+SHRANKAGECKK=3TGSH=7.0TDSH=41.0CC=====================INTIALVALUES=========================C*******FORCREEP(USING1STOPT)*********CCINITIALFUNCTIONISY=1.9*T0^(-0.118)*(T-T0)/(61-0.51*FC+T-T0)CCEQUIVALENTFUNCTION(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,T0ISTHELOADINGTIMECTHEFACTORSARELISTEDASFOLLOWINGCNOTE:FORTHECREEPCOMPLIANCEISRELATIVETOMODULUSOFELASTIC,ALLTHECCOEFFICIENTOFKABIRWILLBEDIFFERENTACCORDINGTOVARYINGMODULUS.CLOADINGTIMEARE7,14,28,60,90,180,365,600,1000,3000(DAYS)C------------------------------------------------------------------------------C|A|7|14|28|60|90|180|365|600|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(1,9)|C(1,10)|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(2,0)|C(2,10)|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(3,9)|C(3,10)|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(4,9)|C(4,10)|C-----------------------------------------------------------------------------CCREEPONLYFORCONCRETEMATERIALCCLOADINGAGE(DAYS)T0=(/7,14,28,60,90,180,365,600,1000,3000/)CIF(CMNAME.EQ.MAT-GIRDER)THENDOI=1,4CLAMDA(I)=GIRD_LAMDA(I)ENDDOCDOJ=1,10DOI=1,4CINI(I,J)=GIRD_INI(I,J)ENDDOENDDOCMODULUSOFELASTICATT0'SAGEEMOD=(/25463,28486,30467,31708,32099,*32506,32717,32799,32851,32902/)ENDIFCIF(CMNAME.EQ.MAT-SLAB)THENDOI=1,4CLAMDA(I)=SLAB_LAMDA(I)ENDDOCDOJ=1,10DOI=1,4CINI(I,J)=SLAB_INI(I,J)ENDDOENDDOCMODULUSOFELASTICATT0'SAGEEMOD=(/23818,26647,28499,29660,30026,*30406,30604,30681,30729,30777/)ENDIFCC========FORSHRINKAGE=======CTGSH--TIMEOFSHRINKAGECALCULATIONFORGIRDER(DAYS)CTDSH--TIMEOFSHRINKAGECALCULATIONFORDECK(DAYS)C480E-6--ULTIMATESHRINKAGESTRAININAASHRO2007(SI)CFC--SPECIFIEDCOMPREI\SSIVESTRENGTHOFCONCRETEATTIMEOFCPRESTRESSINGFORPRESTRESSIONEDMEMBERSANDATTIMEOFCINITIALLOADINGFORNONPRESTRESSIONEDMEMBERS,fci'INAASHTO(2007)CSTANDARDVALUEIS28MPaCCKF--FACTORFORTHEEFFECTOFCONCRETESTRENGTHCCKTD--TIMEDEPENDENT(DEVELOPMENT)FACTOR(THEFIRSTTWOTURMS)CRHSH--AMBIENTRELERTIVEHUMIDITYCORRECIONFACTORFORSHRINKAGESHU=480E-6RH=70.0FC_GIRD=40.0T_PRESTRESS=7.0FC=FC_GIRD*T_PRESTRESS/(4+0.85*T_PRESTRESS)IF(CMNAME.EQ.MAT-SLAB)THENRH=90FC=0.8*FCENDIFCKF=35.0/(7.0+FC)CKTD=61.0-0.58*FCRHSH=2.0-0.014*RHCC=================INTERPOSITION=====================C****COMPUTINGTHECOEFFICIENTSOFKABIRSERIES****CCSTATEV(1)--THERMALSTRAINCSTATEV(2)--ELASTICSTRAINOFCURRENTINCREMENTCSTATEV(3)--ELASTICSTRAINOFPREVIOUSINCREAMENTCSTATEV(4)--DTIMEOFPREVIOUSINCREAMENTCDELTEE--INCREMENTO