| 查看: 2231 | 回复: 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人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有13人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有15人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有14人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有17人回复
-
基金系统什么内容也没有
已经有10人回复
-
学科评审组评审是指会评吗？
已经有5人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
20