




版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
經(jīng)過進(jìn)一步的完善后,子程序的計(jì)算結(jié)果跟一般程序的計(jì)算結(jié)果吻合得比較好。目前在三維實(shí)體元中應(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ù)單獨(dú)準(zhǔn)備creep.inc第一行為文字說明,第二行為Kabir‘s公式中的lamda,其余行為kabir’s公式中的系數(shù)A,每一行對(duì)應(yīng)相
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 智能門禁系統(tǒng)安裝合同
- 教師職稱評(píng)定育人工作證明
- Tricyclohexylphosphine-Tricyclohexylphosphane-生命科學(xué)試劑-MCE
- 影視劇制作投資拍攝協(xié)議
- Actein-Standard-生命科學(xué)試劑-MCE
- 寓言故事愚公移山的教育意義深度解讀
- 美容美發(fā)產(chǎn)品使用效果免責(zé)承諾書
- 水電站落水孔清洗施工方案
- 蘇州玻璃鋼化糞池施工方案
- 2025年滁州c1貨運(yùn)上崗證模擬考試
- 艾梅乙的實(shí)驗(yàn)室診斷與溝通
- 《黑駿馬》讀書筆記思維導(dǎo)圖
- 煤礦建設(shè)安全規(guī)范解讀
- 電力安全工作規(guī)程線路部分
- 元宵節(jié)英文介紹ppt
- 盤扣式懸挑腳手架專項(xiàng)施工方案
- 勞動(dòng)防護(hù)用品知識(shí)考試試題(含答案)
- 高中教師業(yè)務(wù)知識(shí)考試 數(shù)學(xué)試題及答案
- 湘科版科學(xué)(2017)五年級(jí)下冊(cè)全冊(cè)教案
- GB/T 22751-2008臺(tái)球桌
- GB/T 13234-2009企業(yè)節(jié)能量計(jì)算方法
評(píng)論
0/150
提交評(píng)論