资源描述:
《abaqus溷凝土徐变计算》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库。
1、abaqus混凝土徐变计算2.0此为1.0版的修正版,可以考虑混凝土弹性模量随时间,应力变化等情况下的徐变。可以考虑徐变恢复。注意,getvrm中的变量编号与坐标系有关。否则,结果不正确。 SUBROUTINEUSDFLD(FIELD,STATEV,PNEWDT,DIRECT,T,CELENT, 1TIME,DTIME,CMNAME,ORNAME,NFIELD,NSTATV,NOEL,NPT,LAYER, 2KSPT,KSTEP,KINC,NDI,NSHR,COORD,JMAC,
2、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(10),JARRAY(10),JMAC(*),JMATYP(*),
3、1COORD(*)CC Readinginstantaneousthermalstrainindirection11(xaxial)C StoringthethermalstraininstatevariableC CALLGETVRM('THE',ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP, 1MATLAYO,LACCFLA) STATEV(1)=ARRAY(1) CC Readinginstantaneouselas
4、ticxaxialstrainC StoringthethermalstraininstatevariableC NOTE:ARRAY(1)--XAXIALC ARRAY(2)--YAXIALC ARRAY(3)--ZAXIALC CALLGETVRM('EE',ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP, 1MATLAYO,LACCFLA) STATEV(2)=ARRAY(3)C CC
5、ModifyingFieldvariablestomodeldevelopmentofYoung’sC moduluswithage.Cc IF(TIME(2).GT.1)THENc FIELD(1)=TIME(2)c ELSEc FIELD(1)=0c ENDIF RETURN ENDCCC SUBROUTINEUEXPAN(EXPAN,DEXPANDT,TEMP,TIME,DTIME,
6、PREDEF, 1DPRED,STATEV,CMNAME,NSTATV,NOEL)C INCLUDE'ABA_PARAM.INC'C CHARACTER*80CMNAMEC DIMENSIONEXPAN(*),DEXPANDT(*),TEMP(2),TIME(2),PREDEF(*), 1DPRED(*),STATEV(NSTATV),ARRAY(15) REALCINI(4,8),T0(8),CLAMDA(4),AN(4) REA
7、LCCLOAD(4),EMOD(8)CC SWITCH,II=1,CREEP;II=2,SHRINKAGE,II=3,CREEP+SHRANKAGEC KK=3C C=====================INTIALVALUES=========================C*******FORCREEP(USING1STOPT)*********C C INITIALFUNCTIONISY=1.9*T0^(-0.118)*(T-T0)/(61-
8、0.51*FC+T-T0)CC MODELINGFUNCTION(SOFTWAREOF1STOPT)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 WHERETISTHECALCULATINGTIME,T0ISTHELOADINGTI