Abaqus-混凝土徐变收缩子程序-2.1.docx
- 文档编号:30852721
- 上传时间:2024-03-02
- 格式:DOCX
- 页数:9
- 大小:20.44KB
Abaqus-混凝土徐变收缩子程序-2.1.docx
《Abaqus-混凝土徐变收缩子程序-2.1.docx》由会员分享,可在线阅读,更多相关《Abaqus-混凝土徐变收缩子程序-2.1.docx(9页珍藏版)》请在冰豆网上搜索。
经过进一步的完善后,子程序的计算结果跟一般程序的计算结果吻合得比较好。
目前在三维实体元中应用效果良好。
注意:
外部数据文件的数据由其它拟合软件获得。
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
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- Abaqus 混凝土 收缩 子程序 2.1
![提示](https://static.bdocx.com/images/bang_tan.gif)