Abaqus-混凝土徐變收縮子程序

經(jīng)過進一步的完善后，子程序的計算結(jié)果跟一般程序的計算結(jié)果吻合得比較好。目前在三維實體元中應(yīng)用效果良好。注意：外部數(shù)據(jù)文件的數(shù)據(jù)由其它擬合軟件獲得。

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)

C

INCLUDE'ABA_PARAM.INC'

C

CHARACTER*80CMNAME,ORNAME

CHARACTER*3FLGRAY(15)

DIMENSIONFIELD(NFIELD),STATEV(NSTATV),DIRECT(3,3),

1T(3,3),TIME(2)

DIMENSIONARRAY(15),JARRAY(15),JMAC(*),JMATYP(*),

1COORD(*)

C

C

Readinginstantaneousthermalstrainindirection11(xaxial)

C

Storingthethermalstraininstatevariable

C

CALLGETVRM('THE',ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP,

1MATLAYO,LACCFLA)

STATEV(1)=ARRAY(1)

C

C

Readinginstantaneouselasticxaxialstrain

C

Storingthethermalstraininstatevariable

C

NOTE:ARRAY(1)--XAXIAL

C

ARRAY(2)--YAXIAL

C

ARRAY(3)--ZAXIAL

C

CALLGETVRM('EE',ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP,

1MATLAYO,LACCFLA)

STATEV(2)=ARRAY(1)

C

CALLGETVRM('SDV',ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP,

1MATLAYO,LACCFLA)

STATEV(3)=ARRAY(3)

STATEV(4)=ARRAY(4)

C

C

ModifyingFieldvariablestomodeldevelopmentofYoung’s

C

moduluswithage.

C

IF(TIME(2).GT.1)THEN

FIELD(1)=TIME(2)

ELSE

FIELD(1)=0

ENDIF

RETURN

END

C

C

SUBROUTINEUEXTERNALDB(LOP,LRESTART,TIME,DTIME,KSTEP,KINC)

C

INCLUDE'ABA_PARAM.INC'

C

DIMENSIONTIME(2)

CHARACTER*80TITLE

C

COMMON/MYGIRD/GIRD_INI(4,10),GIRD_LAMDA(4)

COMMON/MYSLAB/SLAB_INI(4,10),SLAB_LAMDA(4)

C

C

AbaqususestheFORTRANunitnumbersoutlinedinthetablebelow.Unlessnotedotherwise,

C

youshouldnottrytowritetotheseFORTRANunitsfromusersubroutines.

C

ForAbaqus/Standard,youshouldspecifyunitnumbers15–18orunitnumbersgreaterthan100

C

C

Readthecreepcoefficientsfromexternaldababaseatthebeginningofcomputing.

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)

DOJ=1,10

READ(106,*)(GIRD_INI(I,J),I=1,4)

ENDDO

C

READ(106,*)TITLE

READ(106,*)(SLAB_LAMDA(I),I=1,4)

DOJ=1,10

READ(106,*)(SLAB_INI(I,J),I=1,4)

ENDDO

CLOSE(106)

C

c

OPEN(107,FILE='D:\ABAQUSWORK\CREEP.OUP')

c

Write(107,100)(SLAB_LAMDA(I),I=1,4)

c

DOJ=1,10

c

write(107,100)(SLAB_INI(I,J),I=1,4)

c

ENDDO

c

CLOSE(107)

c100

FORMAT(1X,4E16.7)

ENDIF

RETURN

END

C

SUBROUTINEUEXPAN(EXPAN,DEXPANDT,TEMP,TIME,DTIME,PREDEF,

1DPRED,STATEV,CMNAME,NSTATV,NOEL)

C

INCLUDE'ABA_PARAM.INC'

C

CHARACTER*80CMNAME

C

DIMENSIONEXPAN(*),DEXPANDT(*),TEMP(2),TIME(2),PREDEF(*),

1DPRED(*),STATEV(NSTATV),ARRAY(15)

REALCINI(4,10),T0(10),CLAMDA(4),AN(4)

REALCLOAD(4),EMOD(10)

c

COMMON/MYGIRD/GIRD_INI(4,10),GIRD_LAMDA(4)

COMMON/MYSLAB/SLAB_INI(4,10),SLAB_LAMDA(4)

C

C

SWITCH,II=1,CREEP;II=2,SHRINKAGE,II=3,CREEP+SHRANKAGE

C

KK=3

TGSH=7.0

TDSH=41.0

C

C=====================INTIALVALUES=========================

C*******FORCREEP(USING1STOPT)*********

C

C

INITIALFUNCTIONISY=1.9*T0^(-0.118)*(T-T0)/(61-0.51*FC+T-T0)

C

C

EQUIVALENTFUNCTION(SOFTWAREOF1STOPT)IS

C

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)));

C

C

WHERETISTHECALCULATINGTIME,T0ISTHELOADINGTIME

C

THEFACTORSARELISTEDASFOLLOWING

C

NOTE:FORTHECREEPCOMPLIANCEISRELATIVETOMODULUSOFELASTIC,ALLTHE

C

COEFFICIENTOFKABIRWILLBEDIFFERENTACCORDINGTOVARYINGMODULUS.

C

LOADINGTIMEARE7,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

-----------------------------------------------------------------------------

C

CREEPONLYFORCONCRETEMATERIAL

C

C

LOADINGAGE(DAYS)

T0=(/7,14,28,60,90,180,365,600,1000,3000/)

C

IF(CMNAME.EQ."MAT-GIRDER")THEN

DOI=1,4

CLAMDA(I)=GIRD_LAMDA(I)

ENDDO

C

DOJ=1,10

DOI=1,4

CINI(I,J)=GIRD_INI(I,J)

ENDDO

ENDDO

C

MODULUSOFELASTICATT0'SAGE

EMOD=(/25463,28486,30467,31708,32099,

*

32506,32717,32799,32851,32902/)

ENDIF

C

IF(CMNAME.EQ."MAT-SLAB")THEN

DOI=1,4

CLAMDA(I)=SLAB_LAMDA(I)

ENDDO

C

DOJ=1,10

DOI=1,4

CINI(I,J)=SLAB_INI(I,J)

ENDDO

ENDDO

C

MODULUSOFELASTICATT0'SAGE

EMOD=(/23818,26647,28499,29660,30026,

*

30406,30604,30681,30729,30777/)

ENDIF

C

C

========FORSHRINKAGE=======

C

TGSH--TIMEOFSHRINKAGECALCULATIONFOR"GIRDER"(DAYS)

C

TDSH--TIMEOFSHRINKAGECALCULATIONFOR"DECK"(DAYS)

C

480E-6--ULTIMATESHRINKAGESTRAININAASHRO2007(SI)

C

FC--SPECIFIEDCOMPREI\SSIVESTRENGTHOFCONCRETEATTIMEOF

C

PRESTRESSINGFORPRESTRESSIONEDMEMBERSANDATTIMEOF

C

INITIALLOADINGFORNONPRESTRESSIONEDMEMBERS,fci'INAASHTO(2007)

C

STANDARDVALUEIS28MPa

C

CKF--FACTORFORTHEEFFECTOFCONCRETESTRENGTH

C

CKTD--TIMEDEPENDENT(DEVELOPMENT)FACTOR(THEFIRSTTWOTURMS)

C

RHSH--AMBIENTRELERTIVEHUMIDITYCORRECIONFACTORFORSHRINKAGE

SHU=480E-6

RH=70.0

FC_GIRD=40.0

T_PRESTRESS=7.0

FC=FC_GIRD*T_PRESTRESS/(4+0.85*T_PRESTRESS)

IF(CMNAME.EQ."MAT-SLAB")THEN

RH=90

FC=0.8*FC

ENDIF

CKF=35.0/(7.0+FC)

CKTD=61.0-0.58*FC

RHSH=2.0-0.014*RH

C

C=================INTERPOSITION=====================

C****COMPUTINGTHECOEFFICIENTSOFKABIRSERIES****

C

C

STATEV(1)--THERMALSTRAIN

C

STATEV(2)--ELASTICSTRAINOFCURRENTINCREMENT

C

STATEV(3)--ELASTICSTRAINOFPREVIOUSINCREAMENT

C

STATEV(4)--DTIMEOFPREVIOUSINCREAMENT

C

DELTEE

--INCREMENTOFELASICSTRAIN

TCUR=TIME(2)

THE=STATEV(1)

EECUR=STATEV(2)

EEPRE=STATEV(3)

DTPRE=STATEV(4)

C

DO4IT=1,9

IF((TCUR.GE.T0(IT)).AND.(TCUR.LT.T0(IT+1)))THEN

JJ=IT

ENDIF

4

CONTINUE

TA=T0(JJ)

TB=T0(JJ+1)

EMA=EMOD(JJ)

EMB=EMOD(JJ+1)

C

C

COMPUTNGREALCOEFFICIENTOFC,STORINGINCLOAD(I)

C

IF(TCUR.LT.T0(1))THEN

DO5I=1,4

CLOAD(I)=CINI(I,1)

5

CONTINUE

ELSEIF((TCUR.GE.T0(1)).AND.(TCUR.LT.T0(8)))THEN

DELTAT=TB-TA

C

VARINGMODULUSOFELASTIC

EMODL=(TCUR-TA)/DELTAT*(EMB-EMA)+EMA

C

EMODL=EMA

C

COMPUTINGCOEFFICIENTOFC

DO10I=1,4

CA=CINI(I,JJ)

CB=CINI(I,JJ+1)

CINTERPOSITIONFORANYLOADINGTIME

CLOAD(I)=(TCUR-TA)/DELTAT*(CB-CA)+CA

10

CONTINUE

ELSE

DO20I=1,4

CLOAD(I)=CINI(I,8)

20

CONTINUE

ENDIF

C

DELTEE=EECUR-EEPRE

DELTSTR=DELTEE*EMODL

C

IF((TCUR-T0(1)).LT.1E-5)THEN

DO30I=1,4

STATEV(I+4)=0.0

30

CONTINUE

ELSE

DO40I=1,4

AN(I)=STATEV(I+4)*EXP(-CLAMDA(I)*DTPRE)+

*

DELTSTR*CLOAD(I)

STATEV(I+4)=AN(I)

40

CONTINUE

ENDIF

C

EXPANCR=STATEV(5)*(1-EXP(-CLAMDA(1)*DTIME))+STATEV(6)*

*

(1-EXP(-CLAMDA(2)*DTIME))+STATEV(7)*(1-EXP(-CLAMDA(3)*

*

DTIME))+STATEV(8)*(1-EXP(-CLAMDA(4)*DTIME))

C

IF(TCUR.LT.T0(1))THEN

EXPANCR=0.0

ENDIF

STATEV(10)=EXPANCR

C

C

======================SHRINKAGE=========================

C

C

COMPUTINGSHRINKAGE(FORMULA5.4.2.3.3-1INAASHTO(2007))

C

IF(CMNAME.EQ."MAT-GIRDER")THEN

IF(TIME(2).GT.TGSH)THEN

TSH=TIME(2)-TGSH

TPSH=TSH-DTIME

IF(TPSH.LT.0)THEN

TPSH=0.0

ENDIF

EXPANSH=SHU*RHSH*CKF*(-TSH/(CKTD+TSH)+TPSH/(CKTD+TPSH))

ELSE

EXPANSH=0.0

ENDIF

STATEV(11)=EXPANSH

ENDIF

C

IF(CMNAME.EQ."MAT-SLAB")THEN

IF(TIME(2).GT.TDSH)THEN

TSH=TIME(2)-TDSH

TPSH=TSH-DTIME

IF(TPSH.LT.0)THEN

TPSH=0.0

ENDIF

EXPANSH=SHU*RHSH*CKF*(-TSH/(CKTD+TSH)+TPSH/(CKTD+TPSH))

ELSE

EXPANSH=0.0

ENDIF

STATEV(11)=EXPANSH

ENDIF

C

C

FOROUTPUT,KK=1FORCREEPONLY

C

KK=2FORSHRINKAGEONLY

C

KK=3FORSUMOFCREEPANDSHRINKAGE

IF((KK-1).LT.1E-6)THEN

EXPAN(1)=EXPANCR

ELSEIF((KK-2).LT.1E-6)THEN

EXPAN(1)=EXPANSH

ELSE

EXPAN(1)=EXPANCR+EXPANSH

ENDIF

C

C

STORINGCURRENTSTRAININTOSTATEV(3)

C

STORINGCURRENTDTIMEINTOSTATEV(4)

C

STATEV(3)=EECUR

STATEV(4)=DTIME

RETURN

END

Cc

外部數(shù)據(jù)文件，數(shù)據(jù)單獨準備creep.inc第一行為文字說明，第二行為Kabir‘s公式中的lamda，其余行為kabir’s公式中的系數(shù)A，每一行對應(yīng)相