Abaqus 混凝土徐变收缩子程序.docx
- 文档编号:7262227
- 上传时间:2023-01-22
- 格式:DOCX
- 页数:12
- 大小:18.14KB
Abaqus 混凝土徐变收缩子程序.docx
《Abaqus 混凝土徐变收缩子程序.docx》由会员分享,可在线阅读,更多相关《Abaqus 混凝土徐变收缩子程序.docx(12页珍藏版)》请在冰豆网上搜索。
Abaqus混凝土徐变收缩子程序
Abaqus混凝土徐变收缩子程序
经过进一步的完善后,子程序的计算结果跟一般程序的计算结果吻合得比较好。
目前在三维实体元中应用效果良好。
注意:
外部数据文件的数据由其它拟合软件获得。
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
CReadinginstantaneousthermalstrainindirection11(xaxial)
CStoringthethermalstraininstatevariable
C
CALLGETVRM('THE',ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP,1MATLAYO,LACCFLA)
STATEV
(1)=ARRAY
(1)
C
CReadinginstantaneouselasticxaxialstrain
CStoringthethermalstraininstatevariable
CNOTE:
ARRAY
(1)--XAXIAL
CARRAY
(2)--YAXIAL
CARRAY(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
CModifyingFieldvariablestomodeldevelopmentofYoung’s
Cmoduluswithage.
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
CAbaqususestheFORTRANunitnumbersoutlinedinthetablebelow.Unlessnotedotherwise,
CyoushouldnottrytowritetotheseFORTRANunitsfromusersubroutines.
CForAbaqus/Standard,youshouldspecifyunitnumbers15–18orunitnumbersgreaterthan100
C
CReadthecreepcoefficientsfromexternaldababaseatthebeginningofcomputing.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
cOPEN(107,FILE='D:
\ABAQUSWORK\CREEP.OUP')
cWrite(107,100)(SLAB_LAMDA(I),I=1,4)
cDOJ=1,10
cwrite(107,100)(SLAB_INI(I,J),I=1,4)
cENDDO
cCLOSE(107)
c100FORMAT(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
CSWITCH,II=1,CREEP;II=2,SHRINKAGE,II=3,CREEP+SHRANKAGEC
KK=3
TGSH=7.0
TDSH=41.0
C
C=====================INTIALVALUES=========================
C*******FORCREEP(USING1STOPT)*********
C
CINITIALFUNCTIONISY=1.9*T0^(-0.118)*(T-T0)/(61-0.51*FC+T-T0)C
CEQUIVALENTFUNCTION(SOFTWAREOF1STOPT)IS
CY=C1*(1-EXP(-LAMDA1*(X-T)))+C2*(1-EXP(-LAMDA2*(X-T)))+
CC3*(1-EXP(-LAMDA3*(X-T)))+C4*(1-EXP(-LAMDA4*(X-T)));
C
CWHERETISTHECALCULATINGTIME,T0ISTHELOADINGTIME
CTHEFACTORSARELISTEDASFOLLOWING
CNOTE:
FORTHECREEPCOMPLIANCEISRELATIVETOMODULUSOFELASTIC,ALLTHE
CCOEFFICIENTOFKABIRWILLBEDIFFERENTACCORDINGTOVARYINGMODULUS.
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-----------------------------------------------------------------------------
CCREEPONLYFORCONCRETEMATERIAL
C
CLOADINGAGE(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
CMODULUSOFELASTICATT0'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
CMODULUSOFELASTICATT0'SAGE
EMOD=(/23818,26647,28499,29660,30026,
*30406,30604,30681,30729,30777/)
ENDIF
C
C========FORSHRINKAGE=======
CTGSH--TIMEOFSHRINKAGECALCULATIONFOR"GIRDER"(DAYS)
CTDSH--TIMEOFSHRINKAGECALCULATIONFOR"DECK"(DAYS)
C480E-6--ULTIMATESHRINKAGESTRAININAASHRO2007(SI)
CFC--SPECIFIEDCOMPREI\SSIVESTRENGTHOFCONCRETEATTIMEOF
CPRESTRESSINGFORPRESTRESSIONEDMEMBERSANDATTIMEOF
CINITIALLOADINGFORNONPRESTRESSIONEDMEMBERS,fci'IN
AASHTO(2007)
CSTANDARDVALUEIS28MPa
CCKF--FACTORFORTHEEFFECTOFCONCRETESTRENGTH
CCKTD--TIMEDEPENDENT(DEVELOPMENT)FACTOR(THEFIRSTTWOTURMS)
CRHSH--AMBIENTRELERTIVEHUMIDITYCORRECIONFACTORFOR
SHRINKAGE
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
CSTATEV
(1)--THERMALSTRAIN
CSTATEV
(2)--ELASTICSTRAINOFCURRENTINCREMENT
CSTATEV(3)--ELASTICSTRAINOFPREVIOUSINCREAMENT
CSTATEV(4)--DTIMEOFPREVIOUSINCREAMENT
CDELTEE--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
4CONTINUE
TA=T0(JJ)
TB=T0(JJ+1)
EMA=EMOD(JJ)
EMB=EMOD(JJ+1)
C
CCOMPUTNGREALCOEFFICIENTOFC,STORINGINCLOAD(I)C
IF(TCUR.LT.T0
(1))THEN
DO5I=1,4
CLOAD(I)=CINI(I,1)
5CONTINUE
ELSEIF((TCUR.GE.T0
(1)).AND.(TCUR.LT.T0(8)))THENDELTAT=TB-TA
CVARINGMODULUSOFELASTIC
EMODL=(TCUR-TA)/DELTAT*(EMB-EMA)+EMA
CEMODL=EMA
CCOMPUTINGCOEFFICIENTOFC
DO10I=1,4
CA=CINI(I,JJ)
CB=CINI(I,JJ+1)
CINTERPOSITIONFORANYLOADINGTIME
CLOAD(I)=(TCUR-TA)/DELTAT*(CB-CA)+CA
10CONTINUE
ELSE
DO20I=1,4
CLOAD(I)=CINI(I,8)
20CONTINUE
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
30CONTINUE
ELSE
DO40I=1,4
AN(I)=STATEV(I+4)*EXP(-CLAMDA(I)*DTPRE)+*DELTSTR*CLOAD(I)
STATEV(I+4)=AN(I)
40CONTINUE
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
CCOMPUTINGSHRINKAGE(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
CFOROUTPUT,KK=1FORCREEPONLY
CKK=2FORSHRINKAGEONLY
CKK=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
CSTORINGCURRENTSTRAININTOSTATEV(3)
CSTORINGCURRENTDTIMEINTOSTATEV(4)
C
STATEV(3)=EECUR
STATEV(4)=DTIME
RETURN
END
C
c外部数据文件,数据单独准备creep.in
c第一行为文字说明,第二行为Kabir‘s公式中的lamda,其余行为kabir’s公式中的系数A,每一行对应相应的加载龄期。
*2009.10*FORVARYINGMODULUS*CREEPCOEFFICIENTSFOR"MAT-GIRDER"**0.0200975390.0011267560.0617200550.005655398
3.198793E-057.108837E-061.425333E-052.013975E-05
2.237047E-054.972322E-069.971452E-061.408240E-05
1.854555E-054.122967E-068.285373E-061.169368E-05
1.603102E-053.556803E-067.149851E-061.011927E-05
1.503312E-053.332869E-066.705964E-069.485050E-06
1.363170E-053.018741E-066.073185E-068.589683E-06
1.246732E-052.749243E-065.519788E-067.821465E-06
1.169726E-052.579682E-065.200560E-067.362794E-06
1.096448E-052.408366E-064.896739E-066.944929E-06
9.820035E-062.193812E-064.148105E-065.977530E-06
*2009.10**CREEPCOEFFICIENTSFOR"MAT-SLAB"****
0.0011048110.0594278710.0055105690.01946376
7.23804E-061.39894E-052.02925E-053.17882E-05
5.06217E-069.78774E-061.41908E-052.22287E-05
4.19748E-068.13038E-061.17833E-051.84308E-05
3.62171E-067.01451E-061.01943E-051.59352E-05
3.39417E-066.58011E-069.55456E-061.49428E-05
3.07387E-065.96081E-068.65557E-061.35456E-05
2.80141E-065.41592E-067.87781E-061.23924E-05
2.62613E-065.10482E-067.42096E-061.16222E-05
2.45367E-064.80627E-066.99733E-061.08954E-05
2.22476E-064.11174E-066.04019E-069.70505E-06
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- Abaqus 混凝土徐变收缩子程序 混凝土 收缩 子程序
![提示](https://static.bdocx.com/images/bang_tan.gif)