EMA=EMOD(JJ)
EMB=EMOD(JJ+1) C
C COMPUTNG REAL COEFFICIENT OF C, STORING IN CLOAD(I) C
IF(TCUR.LT.T0(1))THEN DO 5 I=1,4
CLOAD(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*(EMB-EMA)+EMA C EMODL=EMA
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
CLOAD(I)=(TCUR-TA)/DELTAT*(CB-CA)+CA 10 CONTINUE ELSE
DO 20 I=1,4
CLOAD(I)=CINI(I,8) 20 CONTINUE END IF C
DELTEE=EECUR-EEPRE DELTSTR=DELTEE*EMODL C
IF((TCUR-T0(1)).LT.1E-5)THEN DO 30 I=1,4 STATEV(I+4)=0.0 30 CONTINUE ELSE DO 40 I=1,4
AN(I)=STATEV(I+4)*EXP(-CLAMDA(I)*DTPRE)+ * DELTSTR*CLOAD(I) STATEV(I+4)=AN(I) 40 CONTINUE END IF
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
STATEV(10)=EXPANCR C
C ======================SHRINKAGE========================= C
C COMPUTING SHRINKAGE (FORMULA 5.4.2.3.3-1 IN AASHTO(2007)) C
IF (CMNAME.EQ.\ 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)) ELSE
EXPANSH=0.0
END IF
STATEV(11)=EXPANSH END IF C
IF (CMNAME.EQ.\ 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)) ELSE
EXPANSH=0.0 END IF
STATEV(11)=EXPANSH END IF 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 RETURN END C
c 外部数据文件，数据单独准备 creep.in
c 第一行为文字说明，第二行为Kabir‘s公式中的lamda，其余行为kabir’s公式中的系数A，每一行对应相应的加载龄期。
*2009.10*FOR VARYING MODULUS*CREEP COEFFICIENTS FOR \ 0.020097539 0.001126756 0.061720055 0.005655398 3.198793E-05 7.108837E-06 1.425333E-05 2.013975E-05 2.237047E-05 4.972322E-06 9.971452E-06 1.408240E-05 1.854555E-05 4.122967E-06 8.285373E-06 1.169368E-05 1.603102E-05 3.556803E-06 7.149851E-06 1.011927E-05 1.503312E-05 3.332869E-06 6.705964E-06 9.485050E-06 1.363170E-05 3.018741E-06 6.073185E-06 8.589683E-06 1.246732E-05 2.749243E-06 5.519788E-06 7.821465E-06 1.169726E-05 2.579682E-06 5.200560E-06 7.362794E-06 1.096448E-05 2.408366E-06 4.896739E-06 6.944929E-06 9.820035E-06 2.193812E-06 4.148105E-06 5.977530E-06 *2009.10**CREEP COEFFICIENTS FOR \0.001104811 0.059427871 0.005510569 0.01946376 7.23804E-06 1.39894E-05 2.02925E-05 3.17882E-05 5.06217E-06 9.78774E-06 1.41908E-05 2.22287E-05 4.19748E-06 8.13038E-06 1.17833E-05 1.84308E-05 3.62171E-06 7.01451E-06 1.01943E-05 1.59352E-05 3.39417E-06 6.58011E-06 9.55456E-06 1.49428E-05 3.07387E-06 5.96081E-06 8.65557E-06 1.35456E-05 2.80141E-06 5.41592E-06 7.87781E-06 1.23924E-05 2.62613E-06 5.10482E-06 7.42096E-06 1.16222E-05 2.45367E-06 4.80627E-06 6.99733E-06 1.08954E-05 2.22476E-06 4.11174E-06 6.04019E-06 9.70505E-06