《abaqus混凝土徐变计算子程序(共8页).docx》由会员分享，可在线阅读，更多相关《abaqus混凝土徐变计算子程序(共8页).docx（8页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、精选优质文档-倾情为你奉上此为1.0版的修正版，可以考虑混凝土弹性模量随时间，应力变化等情况下的徐变。可以考虑徐变恢复。注意，getvrm中的变量编号与坐标系有关。否则，结果不正确。 SUBROUTINE USDFLD(FIELD,STATEV,PNEWDT,DIRECT,T,CELENT, 1 TIME,DTIME,CMNAME,ORNAME,NFIELD,NSTATV,NOEL,NPT,LAYER, 2 KSPT,KSTEP,KINC,NDI,NSHR,COORD,JMAC,JMATYP,MATLAYO, 3 LACCFLA)C INCLUDE ABA_PARAM.INCC CHARACT
2、ER*80 CMNAME,ORNAME CHARACTER*3 FLGRAY(15) DIMENSION FIELD(NFIELD),STATEV(NSTATV),DIRECT(3,3), 1 T(3,3),TIME(2) DIMENSION ARRAY(10),JARRAY(10),JMAC(*),JMATYP(*), 1 COORD(*)CC Reading instantaneous thermal strain in direction 11(x axial)C Storing the thermal strain in state variableC CALL GETVRM(THE,
3、ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP, 1 MATLAYO,LACCFLA) STATEV(1)= ARRAY(1) CC Reading instantaneous elastic x axial strain C Storing the thermal strain in state variableC NOTE: ARRAY(1)-X AXIALC ARRAY(2)-Y AXIALC ARRAY(3)-Z AXIALC CALL GETVRM(EE,ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP, 1 MATLAYO,LACC
4、FLA) STATEV(2)= ARRAY(3)C CC Modifying Field variables to model development of Youngs C modulus with age.Cc IF(TIME(2).GT.1) THENc FIELD(1)=TIME(2)c ELSEc FIELD(1)=0c END IF RETURN ENDCCC SUBROUTINE UEXPAN(EXPAN,DEXPANDT,TEMP,TIME,DTIME,PREDEF, 1 DPRED,STATEV,CMNAME,NSTATV,NOEL)C INCLUDE ABA_PARAM.I
5、NCC CHARACTER*80 CMNAMEC DIMENSION EXPAN(*),DEXPANDT(*),TEMP(2),TIME(2),PREDEF(*), 1 DPRED(*),STATEV(NSTATV),ARRAY(15) REAL CINI(4,8),T0(8),CLAMDA(4),AN(4) REAL CCLOAD(4),EMOD(8)CC SWITCH, II=1, CREEP; II=2, SHRINKAGE, II=3, CREEP+SHRANKAGEC KK=3 C C= INTIAL VALUES =C * FOR CREEP (USING 1STOPT) *C C
6、 INITIAL FUNCTION IS Y=1.9*T0(-0.118)*(T-T0)/(61-0.51*FC+T-T0)CC MODELING FUNCTION (SOFTWARE OF 1STOPT) ISC 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); CC WHERE T IS THE CALCULATING TIME, T0 IS THE LOADING TIMEC THE FACTORS ARE LISTED AS FOLLO
7、WINGC LOADING TIME ARE 7,30,50,100,365,1000,3000(DAYS)C -C | A | 7 | 30 | 50 | 100 | 365 | 500 | 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 | 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 | 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
8、)|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 - IF(CMNAME.EQ.MAT-GIRDER).OR.(CMNAME.EQ.MAT-SLAB)THENC ELASTIC MODULUSc EMOD=1.0 C VALUES OF LAMDA CLAMDA(1)=0.63775 CLAMDA(2)=0. CLAMDA(3)=0. CLAMDA(4)=0.75493C INITIAL VALUES OF COEFFICIENT C1 CINI(1,1)=0.E-4 CINI(1,2)=0.E-4 CIN
9、I(1,3)=0.E-4 CINI(1,4)=0.E-4 CINI(1,5)=0.E-4 CINI(1,6)=0.E-4 CINI(1,7)=0.E-4 CINI(1,8)=0.E-4 C INITIAL VALUES OF COEFFICIENT C2 CINI(2,1)=0.E-4 CINI(2,2)=0.E-4 CINI(2,3)=0.E-4 CINI(2,4)=0.E-4 CINI(2,5)=0.E-4 CINI(2,6)=0.E-4 CINI(2,7)=0.E-4 CINI(2,8)=0.E-4 C INITIAL VALUES OF COEFFICIENT C3 CINI(3,1)
10、=0.E-4 CINI(3,2)=0.E-4 CINI(3,3)=0.E-4 CINI(3,4)=0.E-4 CINI(3,5)=0.E-4 CINI(3,6)=0.E-4 CINI(3,7)=0.E-4 CINI(3,8)=0.E-4C INITIAL VALUES OF COEFFICIENT C4 CINI(4,1)=0.E-4 CINI(4,2)=0.E-4 CINI(4,3)=0.E-4 CINI(4,4)=0.E-4 CINI(4,5)=0.E-4 CINI(4,6)=0.E-4 CINI(4,7)=0.E-4 CINI(4,8)=0.E-4 C C LOADING AGE(DAY
11、S) T0=(/7,30,50,100,365,500,1000,3000/)C MODULUS OF ELASTIC AT T0S AGE EMOD=(/2.,3.,3.,3., * 3.,3.,3.,3./)C C = FOR SHRINKAGE =C TGSH- TIME OF SHRINKAGE CALCULATION FOR GIRDER (DAYS)C TDSH- TIME OF SHRINKAGE CALCULATION FOR DECK (DAYS)C 480E-6- ULTIMATE SHRINKAGE STRAIN IN AASHRO 2007(SI)C FC -SPECI
12、FIED COMPREISSIVE STRENGTH OF CONCRETE AT TIME OF C PRESTRESSING FOR PRESTRESSIONED MEMBERS AND AT TIME OF C INITIAL LOADING FOR NONPRESTRESSIONED MEMBERS, fciIN AASHTO(2007)C STANDARD VALUE IS 28MPaC CKF -FACTOR FOR THE EFFECT OF CONCRETE STRENGTHC CKTD-TIME DEPENDENT(DEVELOPMENT) FACTOR (THE FIRST
13、 TWO TURMS) C RHSH-AMBIENT RELERTIVE HUMIDITY CORRECION FACTOR FOR SHRINKAGE TGSH=1.0 TDSH=50.0 SHU=480E-6 RH=70.0 FC=28.0 CKF=35.0/(7.0+FC) CKTD=61.0-0.58*FC RHSH=2.0-0.014*RH IF (CMNAME.EQ.MAT-SLAB) THEN CKF=35.0/(7.0+0.8*FC) CKTD=61.0-0.58*0.8*FC END IF C C=INTERPOSITION=C *COMPUTING THE COEFFICI
14、ENTS OF KABIR SERIES*CC STATEV(1)-THERMAL STRAINC STATEV(2)-ELASTIC STRAIN OF CURRENT INCREMENTC STATEV(3)-ELASTIC STRAIN OF PREVIOUS INCREAMENTC STATEV(4)-DTIME OF PREVIOUS INCREAMENTC DELTEE -INCREMENT OF ELASIC STRAIN TCUR=TIME(2) THE=STATEV(1) EECUR=STATEV(2) EEPRE=STATEV(3) DTPRE=STATEV(4)C IF(
15、TCUR.GE.T0(1).AND.(TCUR.LT.T0(2)THEN JJ=1 TA=T0(JJ) TB=T0(JJ+1) EMA=EMOD(JJ) EMB=EMOD(JJ+1) ELSE IF(TCUR.GE.T0(2).AND.(TCUR.LT.T0(3)THEN JJ=2 TA=T0(JJ) TB=T0(JJ+1) EMA=EMOD(JJ) EMB=EMOD(JJ+1) ELSE IF(TCUR.GE.T0(3).AND.(TCUR.LT.T0(4)THEN JJ=3 TA=T0(JJ) TB=T0(JJ+1) EMA=EMOD(JJ) EMB=EMOD(JJ+1) ELSE IF(
16、TCUR.GE.T0(4).AND.(TCUR.LT.T0(5)THEN JJ=4 TA=T0(JJ) TB=T0(JJ+1) EMA=EMOD(JJ) EMB=EMOD(JJ+1) ELSE IF(TCUR.GE.T0(5).AND.(TCUR.LT.T0(6)THEN JJ=5 TA=T0(JJ) TB=T0(JJ+1) EMA=EMOD(JJ) EMB=EMOD(JJ+1) ELSE IF(TCUR.GE.T0(6).AND.(TCUR.LT.T0(7)THEN JJ=6 TA=T0(JJ) TB=T0(JJ+1) EMA=EMOD(JJ) EMB=EMOD(JJ+1) ELSE IF(
17、TCUR.GE.T0(7).AND.(TCUR.LT.T0(8)THEN JJ=7 TA=T0(JJ) TB=T0(JJ+1) EMA=EMOD(JJ) EMB=EMOD(JJ+1) END IFCC COMPUTNG REAL COEFFICIENT OF C, STORING IN CCLOAD(I)C IF(TCUR.LT.T0(1)THEN DO 5 I=1,4 CCLOAD(I)=CINI(I,1)5 CONTINUE ELSE IF(TCUR.GE.T0(1).AND.(TCUR.LT.T0(8)THEN DELTAT=TB-TAC VARING MODULUS OF ELASTI
18、C EMODL=(TCUR-TA)/DELTAT*(EMODA-EMODB)+EMODA C COMPUTING COEFFICIENT OF C DO 10 I=1,4 CA=CINI(I,JJ) CB=CINI(I,JJ+1)C INTERPOSITION FOR ANY LOADING TIME CCLOAD(I)=(TCUR-TA)/DELTAT*(CB-CA) +CA 10 CONTINUE ELSE DO 20 I=1,4 CCLOAD(I)=CINI(I,8)20 CONTINUE END IFC DELTEE=EECUR-EEPRE DELTSTR=DELTEE*EMODLC
19、=C COMPUTING THE COEFFICIENT INCLUDING STRESS AN(I)C OPEN(2,ACCESS=APPEND,FILE=C:CREEP.TXT)C IF(TCUR-1.0).LT.1E-5)THEN DO 30 I=1,4 STATEV(I+4)=0.030 CONTINUE END IF C DO 40 I=1,4 AN(I)=STATEV(I+4)*EXP(-CLAMDA(I)*DTPRE)+ * DELTSTR*CCLOAD(I) STATEV(I+4)=AN(I)40 CONTINUE CC WRITE(2,(10X,5F10.6)TIME(2),
20、(STATEV(I+4),I=1,4) c WRITE(2,*) TIME(2) CCLOAD(I)C WRITE(2,(2X,4F10.6)(CINI(I,J),J=1,8),I=1,4) c WRITE(2,*) TIME(2) DTIME DTPRE DELTSTRC WRITE(2,(6X,4F10.6)TIME(2),DTIME, DTPRE, DELTSTRc WRITE(2,*) STATEV(I)C WRITE(2,(10X,5F10.6)TIME(2),(STATEV(I+4),I=1,4)c WRITE(2,*) C EXPANCR=STATEV(5)*(1-EXP(-CL
21、AMDA(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 END IFC END IF ! CORESPONG TO LINE 78S IFC C =SHRINKAGE=CC COMPUTING SHRINKAGE (FORMULA 5.4.2.3.3-1 IN AASHTO(2007) CC IF (CMNAME.EQ.MAT-GI
22、RDER) THEN IF (TIME(2).GT.TGSH)THEN TSH=TIME(2)-TGSH TPSH=TSH-DTIME IF(TPSH.LT.0) THEN TPSH=0.0 END IF EXPANSH=SHU*RHSH*CKF*(-TSH/(CKTD+TSH)+TPSH/(CKTD+TPSH) END IF END IF 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 END IF EXPAN
23、SH=SHU*RHSH*CKF*(-TSH/(CKTD+TSH)+TPSH/(CKTD+TPSH) END IF END IF C WRITE(2,(10X,4F10.6)RHSH,CKF,TSH,TPSH,EXPANSH CC FOR OUTPUT, KK=1 FOR CREEP ONLYC KK=2 FOR SHRINKAGE ONLYC KK=3 FOR SUM OF CREEP AND SHRINKAGE IF(KK-1).LT.1E-6) THEN EXPAN(1)=EXPANCR ELSE IF(KK-2).LT.1E-6)THEN EXPAN(1)=EXPANSH ELSE EXPAN(1)=EXPANCR+EXPANSH END IFC C STORING CURRENT STRAIN INTO STATEV(3)C STORING CURRENT DTIME INTO STATEV(4)C STATEV(3)=EECUR STATEV(4)=DTIMEC CLOSE(2) RETURN ENDC 专心-专注-专业