2022年abaqus混凝土徐变计算子程序. .pdf

此为 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.INC C CHARACTER*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(*)C C Reading instantaneous thermal strain in direction 11(x axial)C Storing the thermal strain in state variable C CALL GETVRM(THE,ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP,1 MATLAYO,LACCFLA)STATEV(1)=ARRAY(1)C C Reading instantaneous elastic x axial strain C Storing the thermal strain in state variable C NOTE:ARRAY(1)-X AXIAL C ARRAY(2)-Y AXIAL C ARRAY(3)-Z AXIAL C CALL GETVRM(EE,ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP,1 MATLAYO,LACCFLA)STATEV(2)=ARRAY(3)C C C Modifying Field variables to model development of Youngs C modulus with age.C c IF(TIME(2).GT.1)THEN c FIELD(1)=TIME(2)c ELSE c FIELD(1)=0 c END IF 名师资料总结-精品资料欢迎下载-名师精心整理-第 1 页，共 8 页 -RETURN END C C C SUBROUTINE UEXPAN(EXPAN,DEXPANDT,TEMP,TIME,DTIME,PREDEF,1 DPRED,STATEV,CMNAME,NSTATV,NOEL)C INCLUDE ABA_PARAM.INC C CHARACTER*80 CMNAME C 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)C C SWITCH,II=1,CREEP;II=2,SHRINKAGE,II=3,CREEP+SHRANKAGE C KK=3 C C=INTIAL VALUES=C*FOR CREEP(USING 1STOPT)*C C INITIAL FUNCTION IS Y=1.9*T0(-0.118)*(T-T0)/(61-0.51*FC+T-T0)C C MODELING FUNCTION(SOFTWARE OF 1STOPT)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 WHERE T IS THE CALCULATING TIME,T0 IS THE LOADING TIME C THE FACTORS ARE LISTED AS FOLLOWING C 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)|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)THEN C ELASTIC MODULUS c EMOD=1.0 名师资料总结-精品资料欢迎下载-名师精心整理-第 2 页，共 8 页 -C VALUES OF LAMDA CLAMDA(1)=0.0193907053463775 CLAMDA(2)=0.00537298967078406 CLAMDA(3)=0.00105487934454054 CLAMDA(4)=0.0597818282775493 C INITIAL VALUES OF COEFFICIENT C1 CINI(1,1)=0.255636119E-4 CINI(1,2)=0.179033082E-4 CINI(1,3)=0.163880576E-4 CINI(1,4)=0.147710753E-4 CINI(1,5)=0.124416704E-4 CINI(1,6)=0.119344460E-4 CINI(1,7)=0.110286669E-4 CINI(1,8)=0.096165142E-4 C INITIAL VALUES OF COEFFICIENT C2 CINI(2,1)=0.157591725E-4 CINI(2,2)=0.110425824E-4 CINI(2,3)=0.101134277E-4 CINI(2,4)=0.091232941E-4 CINI(2,5)=0.077200858E-4 CINI(2,6)=0.074402105E-4 CINI(2,7)=0.068149034E-4 CINI(2,8)=0.060733782E-4 C INITIAL VALUES OF COEFFICIENT C3 CINI(3,1)=0.053080102E-4 CINI(3,2)=0.037166275E-4 CINI(3,3)=0.034018732E-4 CINI(3,4)=0.030634778E-4 CINI(3,5)=0.025744519E-4 CINI(3,6)=0.024744322E-4 CINI(3,7)=0.022776919E-4 CINI(3,8)=0.018527641E-4 C INITIAL VALUES OF COEFFICIENT C4 CINI(4,1)=0.120382947E-4 CINI(4,2)=0.084348462E-4 CINI(4,3)=0.077257019E-4 CINI(4,4)=0.069613683E-4 CINI(4,5)=0.058953358E-4 CINI(4,6)=0.056904171E-4 CINI(4,7)=0.051955726E-4 CINI(4,8)=0.045768253E-4 C C LOADING AGE(DAYS)T0=(/7,30,50,100,365,500,1000,3000/)名师资料总结-精品资料欢迎下载-名师精心整理-第 3 页，共 8 页 -C MODULUS OF ELASTIC AT T0S AGE EMOD=(/2.546295,3.061405,3.147964,3.217927,*3.271749,3.277386,3.285052,3.290193/)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-SPECIFIED 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 28MPa C CKF-FACTOR FOR THE EFFECT OF CONCRETE STRENGTH C CKTD-TIME DEPENDENT(DEVELOPMENT)FACTOR(THE FIRST 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 COEFFICIENTS OF KABIR SERIES*C C STATEV(1)-THERMAL STRAIN C STATEV(2)-ELASTIC STRAIN OF CURRENT INCREMENT C STATEV(3)-ELASTIC STRAIN OF PREVIOUS INCREAMENT C STATEV(4)-DTIME OF PREVIOUS INCREAMENT C DELTEE -INCREMENT OF ELASIC STRAIN TCUR=TIME(2)THE=STATEV(1)EECUR=STATEV(2)EEPRE=STATEV(3)DTPRE=STATEV(4)C 名师资料总结-精品资料欢迎下载-名师精心整理-第 4 页，共 8 页 -IF(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(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(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 IF 名师资料总结-精品资料欢迎下载-名师精心整理-第 5 页，共 8 页 -C C 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-TA C VARING MODULUS OF ELASTIC 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 IF C DELTEE=EECUR-EEPRE DELTSTR=DELTEE*EMODL C =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.0 30 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 C C WRITE(2,(10X,5F10.6)TIME(2),(STATEV(I+4),I=1,4)名师资料总结-精品资料欢迎下载-名师精心整理-第 6 页，共 8 页 -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 DELTSTR C WRITE(2,(6X,4F10.6)TIME(2),DTIME,DTPRE,DELTSTR c 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(-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 END IF C END IF!CORESPONG TO LINE 78S IF C C =SHRINKAGE=C C COMPUTING SHRINKAGE(FORMULA 5.4.2.3.3-1 IN AASHTO(2007)C 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 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 EXPANSH=SHU*RHSH*CKF*(-TSH/(CKTD+TSH)+TPSH/(CKTD+TPSH)END IF 名师资料总结-精品资料欢迎下载-名师精心整理-第 7 页，共 8 页 -END IF C WRITE(2,(10X,4F10.6)RHSH,CKF,TSH,TPSH,EXPANSH C C FOR OUTPUT,KK=1 FOR CREEP ONLY C KK=2 FOR SHRINKAGE ONLY C 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 IF C C STORING CURRENT STRAIN INTO STATEV(3)C STORING CURRENT DTIME INTO STATEV(4)C STATEV(3)=EECUR STATEV(4)=DTIME C CLOSE(2)RETURN END C名师资料总结-精品资料欢迎下载-名师精心整理-第 8 页，共 8 页 -