24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2167  |  回复: 0

德群

新虫 (小有名气)

[求助] 混凝土冻融循环下温度应力的模拟

各位大神好,本人想要模拟混凝土受冻融情况下的温度应力和温度应变,下面的是我的程序。我看资料自己的理解是将温度作为荷载加到混凝土上,然后由于混凝土线膨胀系数的不同产生约束,导致最终产生应力和应变。我选择time history variable里面有stress 和elastic strain 这个是温度产生的温度应力和由此产生的应变么
FINISH
/CLEAR
/FILNAME,SHIYAN quater
/PREP7
ET,1,SOLID70
MP,KXX,1,1.8
MP,KYY,1,1.8
MP,KZZ,1,1.8
MP,DENS,1,2400
MP,C,1,950

MP,KXX,2,0.61
MP,KYY,2,0.61
MP,KZZ,2,0.61
MP,DENS,2,996
MP,C,2,4185

BLOCK,0,-57.5E-3,0,-57.5E-3,0,450E-3
WPAVE,0,0,400E-3
CSYS,4
VSBW,ALL
WPCSYS,-1,0

WPROTA,0,-90,0
WPAVE,0,0,-50E-3
CSYS,4
VSBW,ALL
WPCSYS,-1,0

WPROTA,0,0,90
WPAVE,0,0,-50E-3
CSYS,4
VSBW,ALL
WPCSYS,-1,0

ALLSEL,ALL
VSEL,S,LOC,Z,0,400E-3
VSEL,R,LOC,X,-50E-3,0
VSEL,R,LOC,Y,-50E-3,0
VATT,1,,1
ALLSEL,ALL
VSEL,U,MAT,,1,1,1
VATT,2,,1
ALLSEL,ALL

ESIZE,75E-4
MSHKEY,1
MSHAPE,0
VMESH,ALL

/SOLU
ALLSEL
VSEL,S,MAT,,1                  
NSLV,S,1
D,ALL,TEMP,5.08              
ALLSEL,ALL
VSEL,S,MAT,,2
NSLV,S,1
D,ALL,TEMP,2.88                  
ALLSEL,ALL
ANTYPE, TRANS
TIMINT,OFF
TIME,0.01
!  DELTIM,0.01,    0.02,0.005     !***  设置为一个子步即可,不用 DELTIM 命令,改为下一句
nsubst, 1,1,1
SOLVE

*DIM,SHIJ,,1000
/INPUT,SJ,TXT
*DIM,WENDU,,1000
/INPUT,WD,TXT

/SOLU
*DO,I,1,1000,1
TIME,SHIJ(I,1,1)*60
TIMINT,ON
AUTOTS,ON
KBC,1
NSUBST,1,1
DDELE,ALL,TEMP
ALLSEL,ALL
ASEL,S,AREA,,26
ASEL,A,AREA,,22
ASEL,A,AREA,,4
ASEL,A,AREA,,18
ASEL,A,AREA,,31
ASEL,A,AREA,,34
ASEL,A,AREA,,38
ASEL,A,AREA,,41
NSLA,S,1
D,ALL,TEMP,WENDU(I,1,1)
OUTRES,ALL,1
ALLSEL,ALL
SOLVE
*ENDDO
FINISH

/PREP7
ETCHG,TTS
MP,EX,1,3E10
MP,ALPX,1,10E-6
MP,ALPY,1,10E-6
MP,ALPZ,1,10E-6
MP,NUXY,1,0.2  

MP,EX,2,2.18E9
MP,ALPX,2,0.0002
MP,ALPY,2,0.0002
MP,ALPZ,2,0.0002
MP,NUXY,2,0.2


ALLSEL,ALL
ASEL,S,AREA,,13
NSLA,S,1
D,ALL,UZ,0  
ALLSEL,ALL
ASEL,S,AREA,,5
NSLA,S,1
DSYM,SYMM,X,
ALLSEL,ALL
ASEL,S,AREA,,15
NSLA,S,1
DSYM,SYMM,Y
ALLSEL,ALL         

*DIM,SHIJIAN,,1000
/INPUT,JIANSHI,TXT
/SOLU
ANTYPE,TRANS
*DO,I,1,1000,1
TIME,SHIJIAN(I,1,1)*60
TIMINT,ON
AUTOTS,ON
NSUBST,1,1
NLGEOM,ON
NROPT,FULL

ALLSEL,ALL
ESEL,S,MAT,,2
EKILL,ALL
ALLSEL,ALL
ESEL,S,MAT,,1
NSLE,S
LDREAD,TEMP,,,SHIJIAN(I,1,1)*60,,SHIYAN quater,Rth              
ALLSEL,ALL
OUTRES,ALL
SOLVE
*ENDDO
FINISH

混凝土冻融循环下温度应力的模拟
time history variable.png@wuming524
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 德群 的主题更新
信息提示
请填处理意见