| 查看: 2218 | 回复: 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 |
回复此楼» 猜你喜欢
为什么网上很多人说本周 12号出结果
已经有4人回复
-
基金中了
已经有15人回复
-
帮忙看看fileCode
已经有8人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有3人回复
-
FileCode能看出啥？
已经有15人回复
-
filecode
已经有8人回复
-
确定了，国自然21号放榜
已经有7人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有11人回复
-
关于代码变化问题，想知道的进来
已经有24人回复
-
静等基金结果
已经有16人回复
100