| 查看: 1435 | 回复: 8 | ||
[求助]
帮忙看一下我做的水化热应力场命令流对不对？
|
|
finish/clear/filname,xianghan/title,hydration heat of concretekeyw,pr_therm,1/prep7!单位，长度：m，力:ncha=0.3$d=0.4$b=6$h=3.6$l=2hs=0.4$bs=0.15$ts=0.01$r1=0.05$r2=0.025$r3=0.02k,1k,2,b+d+dk,3,b+d+d,h+d+dk,4,0,h+d+dk,5,d+cha,dk,6,b+d-cha,dk,7,b+d,d+chak,8,b+d,d+h-chak,9,b+d-cha,d+hk,10,d+cha,d+hk,11,d,d+h-chak,12,d,d+cha*do,i,1,3l,i,i+1*enddol,4,1*do,i,5,11l,i,i+1*enddol,12,5al,allvext,1,,,,,let,1,solid70keyopt,1,2,0 Ec=3.25e10 *dim,exx,array,103 *do,i,3,103,1t=(i-2)/24exx(i)=Ec*exp(0.06*(1-(28/t)**0.5)) *enddo *do,i,3,103,1mp,dens,i,2400mp,kxx,i,9500 mp,kyy,i,9500mp,kzz,i,9500mp,c,i,930mp,prxy,i,0.167mp,alpx,i,0.8e-5mp,ex,i,exx(i)*enddovsel,allvatt,3,,1allselWPROTA,,,90$WPOFF,,,0.2$VSBW,ALLWPOFF,,,0.2$VSBW,ALLWPOFF,,,0.3$VSBW,ALLWPOFF,,,2.7$VSBW,ALLWPOFF,,,2.7$VSBW,ALLWPOFF,,,0.3$VSBW,ALLWPOFF,,,0.2$VSBW,ALLWPOFF,,0.2$WPROTA,,90$VSBW,ALLWPOFF,,,-0.2$VSBW,ALLWPOFF,,,-0.3$VSBW,ALLWPOFF,,,-1.5$VSBW,ALLWPOFF,,,-1.5$VSBW,ALLWPOFF,,,-0.3$VSBW,ALLWPOFF,,,-0.2$VSBW,ALLWPCSYS,-1,0!网格划分ESIZE,0.2$VATT,3,1,1$MSHKEY,1VMESH,ALL*get,NUM,elem,,count NUMMRG,NODEnumcmp,node !=========================（二）热分析求解获取温度场=========================!!可以根据材料编号和实常数号选择实体，方便!初始温度，即浇筑温度allsel,allinitial=20/soluantype,transNROPT,FULL !牛顿拉夫逊法TIMINT,ON !打开时间积分效应AUTOTS,ON !打开自动时间步长allsel,ALL,NODEIC,ALL,temp,initial !设置节点初始温度*dim,hef1,ARRAY,168 hef1(1)=0.511,0.274,0.204,0.175,0.181,0.181,0.183,0.19,0.204,0.23hef1(11)=0.273,0.349,0.48,0.697,1.026,1.462,1.957,2.415,2.744,2.902hef1(21)=2.915,2.862,2.833,2.77,2.637,2.46,2.26,2.055,1.859,1.684hef1(31)=1.53,1.399,1.289,1.195,1.117,1.051,0.995,0.954,0.902,0.863hef1(41)=0.828,0.795,0.766,0.74,0.714,0.691,0.669,0.648,0.628,0.609hef1(51)=0.59,0.572,0.554,0.537,0.522,0.503,0.487,0.47,0.454,0.439hef1(61)=0.424,0.409,0.394,0.38,0.366,0.352,0.339,0.327,0.315,0.303hef1(71)=0.293,0.282,0.12240,0.12140,0.11789,0.11353,0.10878,0.10423,0.09837,0.09268hef1(81)=0.08902,0.08561,0.08263,0.07945,0.07643,0.07327,0.06991,0.06688,0.06303,0.05802hef1(91)=0.02948,0.01151,0.01599,0.04059,0.05279,0.05642,0.05686,0.05619,0.05424,0.05563hef1(101)=0,05502,0.02433,0.018,0.01278,0.01285,0.00962,0.00600,0.00303,0.00002,0.00206hef1(111)=0.00375,0.00513,0.00428,0.00460,0.00907,0.01136,0.01209,0.01130,0.01093,0.01103hef1(121)=0.00964,0.00879,0.00880,0.00910,0.00922,0.00915,0.00267,0.00179,0.00234,0.00206hef1(131)=0.00147,0.00111,0.00100,0.00129,0.00189,0.00256,0.00374,0.00045,0.00167,0.00336hef1(141)=0.00418,0.00294,0.01842,0.03276,0.03856,0.04180,0.03701,0.03135,0.02704,0.02379hef1(151)=0.01992,0.01594,0.01358,0.01168,0.01000,0.00816,0.00654,0.00527,0.00367,0.00200hef1(161)=0.00011,0.00167,0.00155,0.00328,0.00455,0.00442,0.00215,0.00008!1h至100h，钢模*do,h,1,100,1time,hheat=h+2 !设置瞬态分析步时间shuihr=hef1(heat)*3600*1350 !水化热生热率，不是总热量，总热量对t求导。W/m3esel,s,mat,,heat !混凝土材料号从3开始，选择混凝土。BFE,all,HGEN,,shuihr !在单元上施加水化热放热 !BFE命令为在单元上施加体荷载的命令,HGEN为热分析时体荷载的标识符allselasel,s,,,18asel,a,,,5asel,a,,,33asel,a,,,42asel,a,,,51asel,a,,,59asel,a,,,66asel,a,,,55 !选择预制顶面；空气nsla,s,1 sf,all,conv,35330,20allselasel,u,,,18asel,u,,,5asel,u,,,33asel,u,,,42asel,u,,,51asel,u,,,59asel,u,,,66asel,u,,,55 !选择现浇侧面；钢模板nsla,s,1 sf,all,conv,35174,20allsel,allSOLVEesel,s,mat,,heat !在之前材料号（混凝土号从3开始）的基础上选择混凝土，再改变其属性进入下一轮;mpchg,(h+3),all *enddo!=========================（三）结构分析求解获取应力场=========================!/prep7 allsel,allsfdele,all,all !删除节点上的面荷载bfdele,all,all !删除面上的体载荷 et,2,solid65 !温度场分析用SOLID70，应力分析改为SOLID185，单元性质改变，定义SOLID185esel,s,type,,1emodif,all,type,2 !将混凝土单元改为SOLID185*do,i,3,103,1mp,dens,i,2400mp,prxy,i,0.167mp,alpx,i,0.8e-5mp,ex,i,exx(i)*enddoesel,s,mat,,103 mpchg,3,all allselasel,s,,,54asel,a,,,65asel,a,,,48asel,a,,,35asel,a,,,29asel,a,,,22asel,a,,,3asel,a,,,17asel,a,,,26asel,a,,,47asel,a,,,39asel,a,,,46nsla,s,1 d,all,uy,0allselasel,s,,,231asel,a,,,248asel,a,,,208asel,a,,,176asel,a,,,163asel,a,,,141asel,a,,,129asel,a,,,107asel,a,,,86asel,a,,,247asel,a,,,207asel,a,,,145asel,a,,,158asel,a,,,115asel,a,,,126asel,a,,,80asel,a,,,164asel,a,,,138asel,a,,,167asel,a,,,153nsla,s,1d,all,ux,0allselasel,s,,,232asel,a,,,239asel,a,,,236asel,a,,,249asel,a,,,252asel,a,,,214asel,a,,,209asel,a,,,212asel,a,,,28asel,a,,,177asel,a,,,182asel,a,,,161asel,a,,,154asel,a,,,121asel,a,,,76asel,a,,,14asel,a,,,130asel,a,,,134asel,a,,,123asel,a,,,108asel,a,,,116asel,a,,,104asel,a,,,102asel,a,,,99asel,a,,,91asel,a,,,95asel,a,,,1asel,a,,,70asel,a,,,119asel,a,,,67asel,a,,,53asel,a,,,127asel,a,,,88asel,a,,,62asel,a,,,81asel,a,,,139asel,a,,,159asel,a,,,180asel,a,,,146asel,a,,,170asel,a,,,30asel,a,,,57asel,a,,,142asel,a,,,243asel,a,,,245asel,a,,,20asel,a,,,185asel,a,,,205asel,a,,,229asel,a,,,222asel,a,,,225asel,a,,,218asel,a,,,4asel,a,,,58asel,a,,,31asel,a,,,171asel,a,,,246asel,a,,,244asel,a,,,143asel,a,,,206asel,a,,,186asel,a,,,21asel,a,,,147asel,a,,,181asel,a,,,160asel,a,,,150asel,a,,,120asel,a,,,71asel,a,,,128asel,a,,,64asel,a,,,68asel,a,,,82asel,a,,,63asel,a,,,89asel,a,,,92asel,a,,,96asel,a,,,103asel,a,,,100asel,a,,,15asel,a,,,84asel,a,,,140asel,a,,,124asel,a,,,135asel,a,,,131asel,a,,,105asel,a,,,117asel,a,,,109asel,a,,,155asel,a,,,162asel,a,,,183asel,a,,,178asel,a,,,237asel,a,,,240asel,a,,,233asel,a,,,215asel,a,,,253asel,a,,,250asel,a,,,34asel,a,,,213asel,a,,,210asel,a,,,230asel,a,,,226asel,a,,,223asel,a,,,219 nsla,s,1d,all,uz,0allselasel,s,,,7asel,a,,,9asel,a,,,11asel,a,,,13nsla,s,1 d,all,ux,0,,,uyallsel/soluallsel,allantype,staticacel,0,10,0 *do,i,1,100,1 ldread,temp,i,,,,'xianghan','rth' allsel solveallselesel,s,mat,,i+2mpchg,i+3,all*enddofinishsave 发自小木虫Android客户端 |
回复此楼
youyuwanglin
木虫 (著名写手)
- 应助: 22 (小学生)
- 金币: 17575.7
- 散金: 50
- 红花: 11
- 帖子: 2344
- 在线: 338.1小时
- 虫号: 1371413
- 注册: 2011-08-18
- 专业: 防灾工程
2楼2017-10-01 12:55:56
3楼2017-10-01 18:37:10
zuiaixuexi
至尊木虫 (文坛精英)
- 应助: 41 (小学生)
- 金币: 27261.7
- 散金: 400
- 红花: 51
- 沙发: 63
- 帖子: 31364
- 在线: 360.5小时
- 虫号: 4350002
- 注册: 2016-01-11
- 性别: GG
- 专业: 会计学
4楼2017-10-01 18:51:31
5楼2017-10-11 15:23:22
6楼2017-10-11 15:52:02
|
是这样的： finish /clear /filname,xianghan /title,hydration heat of concrete keyw,pr_therm,1 /prep7 !单位，长度：m，力:n cha=0.3$d=0.4$b=6$h=3.6$l=2 hs=0.4$bs=0.15$ts=0.01$r1=0.05$r2=0.025$r3=0.02 k,1 k,2,b+d+d k,3,b+d+d,h+d+d k,4,0,h+d+d k,5,d+cha,d k,6,b+d-cha,d k,7,b+d,d+cha k,8,b+d,d+h-cha k,9,b+d-cha,d+h k,10,d+cha,d+h k,11,d,d+h-cha k,12,d,d+cha *do,i,1,3 l,i,i+1 *enddo l,4,1 *do,i,5,11 l,i,i+1 *enddo l,12,5 al,all vext,1,,,,,l et,1,solid70 *do,i,3,103,1 mp,dens,i,2400 mp,kxx,i,9500 mp,kyy,i,9500 mp,kzz,i,9500 mp,c,i,930 *enddo vsel,all vatt,3,,1 allsel WPROTA,,,90$WPOFF,,,0.2$VSBW,ALL WPOFF,,,0.2$VSBW,ALL WPOFF,,,0.3$VSBW,ALL WPOFF,,,2.7$VSBW,ALL WPOFF,,,2.7$VSBW,ALL WPOFF,,,0.3$VSBW,ALL WPOFF,,,0.2$VSBW,ALL WPOFF,,0.2$WPROTA,,90$VSBW,ALL WPOFF,,,-0.2$VSBW,ALL WPOFF,,,-0.3$VSBW,ALL WPOFF,,,-1.5$VSBW,ALL WPOFF,,,-1.5$VSBW,ALL WPOFF,,,-0.3$VSBW,ALL WPOFF,,,-0.2$VSBW,ALL WPCSYS,-1,0 !网格划分 ESIZE,0.2$VATT,3,1,1$MSHKEY,1 VMESH,ALL *get,NUM,elem,,count NUMMRG,NODE numcmp,node !=========================（二）热分析求解获取温度场=========================! !可以根据材料编号和实常数号选择实体，方便 !初始温度，即浇筑温度 allsel,all initial=20 /solu antype,trans NROPT,FULL !牛顿拉夫逊法 TIMINT,ON !打开时间积分效应 AUTOTS,ON !打开自动时间步长 allsel,ALL,NODE IC,ALL,temp,initial !设置节点初始温度 *dim,hef1,ARRAY,168 hef1(1)=0.511,0.274,0.204,0.175,0.181,0.181,0.183,0.19,0.204,0.23 hef1(11)=0.273,0.349,0.48,0.697,1.026,1.462,1.957,2.415,2.744,2.902 hef1(21)=2.915,2.862,2.833,2.77,2.637,2.46,2.26,2.055,1.859,1.684 hef1(31)=1.53,1.399,1.289,1.195,1.117,1.051,0.995,0.954,0.902,0.863 hef1(41)=0.828,0.795,0.766,0.74,0.714,0.691,0.669,0.648,0.628,0.609 hef1(51)=0.59,0.572,0.554,0.537,0.522,0.503,0.487,0.47,0.454,0.439 hef1(61)=0.424,0.409,0.394,0.38,0.366,0.352,0.339,0.327,0.315,0.303 hef1(71)=0.293,0.282,0.12240,0.12140,0.11789,0.11353,0.10878,0.10423,0.09837,0.09268 hef1(81)=0.08902,0.08561,0.08263,0.07945,0.07643,0.07327,0.06991,0.06688,0.06303,0.05802 hef1(91)=0.02948,0.01151,0.01599,0.04059,0.05279,0.05642,0.05686,0.05619,0.05424,0.05563 hef1(101)=0,05502,0.02433,0.018,0.01278,0.01285,0.00962,0.00600,0.00303,0.00002,0.00206 hef1(111)=0.00375,0.00513,0.00428,0.00460,0.00907,0.01136,0.01209,0.01130,0.01093,0.01103 hef1(121)=0.00964,0.00879,0.00880,0.00910,0.00922,0.00915,0.00267,0.00179,0.00234,0.00206 hef1(131)=0.00147,0.00111,0.00100,0.00129,0.00189,0.00256,0.00374,0.00045,0.00167,0.00336 hef1(141)=0.00418,0.00294,0.01842,0.03276,0.03856,0.04180,0.03701,0.03135,0.02704,0.02379 hef1(151)=0.01992,0.01594,0.01358,0.01168,0.01000,0.00816,0.00654,0.00527,0.00367,0.00200 hef1(161)=0.00011,0.00167,0.00155,0.00328,0.00455,0.00442,0.00215,0.00008 !1h至100h，钢模 *do,h,1,100,1 time,h heat=h+2 !设置瞬态分析步时间 shuihr=hef1(heat)*3600*1350 !水化热生热率，不是总热量，总热量对t求导。W/m3 esel,s,mat,,heat !混凝土材料号从3开始，选择混凝土。 BFE,all,HGEN,,shuihr !在单元上施加水化热放热 !BFE命令为在单元上施加体荷载的命令,HGEN为热分析时体荷载的标识符 allsel asel,s,,,18 asel,a,,,5 asel,a,,,33 asel,a,,,42 asel,a,,,51 asel,a,,,59 asel,a,,,66 asel,a,,,55 !选择预制顶面；空气 nsla,s,1 sf,all,conv,35330,20 allsel asel,u,,,18 asel,u,,,5 asel,u,,,33 asel,u,,,42 asel,u,,,51 asel,u,,,59 asel,u,,,66 asel,u,,,55 !选择现浇侧面；钢模板 nsla,s,1 sf,all,conv,35174,20 allsel,all SOLVE esel,s,mat,,heat !在之前材料号（混凝土号从3开始）的基础上选择混凝土，再改变其属性进入下一轮; mpchg,(h+3),all *enddo !=========================（三）结构分析求解获取应力场=========================! /prep7 allsel,all sfdele,all,all !删除节点上的面荷载 bfdele,all,all !删除面上的体载荷 et,2,solid65 !温度场分析用SOLID70，应力分析改为SOLID185，单元性质改变，定义SOLID185 esel,s,type,,1 emodif,all,type,2 !将混凝土单元改为SOLID185 Ec=3.25e10 *dim,exx,array,103 *do,i,3,103,1 t=(i-2)/24 exx(i)=Ec*exp(0.06*(1-(28/t)**0.5)) *enddo *do,i,3,103,1 mp,dens,i,2400 mp,prxy,i,0.167 mp,alpx,i,0.8e-5 mp,ex,i,exx(i) mp,reft,i,20 *enddo esel,s,mat,,103 mpchg,3,all allsel asel,s,,,54 asel,a,,,65 asel,a,,,48 asel,a,,,35 asel,a,,,29 asel,a,,,22 asel,a,,,3 asel,a,,,17 asel,a,,,26 asel,a,,,47 asel,a,,,39 asel,a,,,46 nsla,s,1 d,all,uy,0 allsel asel,s,,,231 asel,a,,,248 asel,a,,,208 asel,a,,,176 asel,a,,,163 asel,a,,,141 asel,a,,,129 asel,a,,,107 asel,a,,,86 asel,a,,,247 asel,a,,,207 asel,a,,,145 asel,a,,,158 asel,a,,,115 asel,a,,,126 asel,a,,,80 asel,a,,,164 asel,a,,,138 asel,a,,,167 asel,a,,,153 nsla,s,1 d,all,ux,0 allsel asel,s,,,232 asel,a,,,239 asel,a,,,236 asel,a,,,249 asel,a,,,252 asel,a,,,214 asel,a,,,209 asel,a,,,212 asel,a,,,28 asel,a,,,177 asel,a,,,182 asel,a,,,161 asel,a,,,154 asel,a,,,121 asel,a,,,76 asel,a,,,14 asel,a,,,130 asel,a,,,134 asel,a,,,123 asel,a,,,108 asel,a,,,116 asel,a,,,104 asel,a,,,102 asel,a,,,99 asel,a,,,91 asel,a,,,95 asel,a,,,1 asel,a,,,70 asel,a,,,119 asel,a,,,67 asel,a,,,53 asel,a,,,127 asel,a,,,88 asel,a,,,62 asel,a,,,81 asel,a,,,139 asel,a,,,159 asel,a,,,180 asel,a,,,146 asel,a,,,170 asel,a,,,30 asel,a,,,57 asel,a,,,142 asel,a,,,243 asel,a,,,245 asel,a,,,20 asel,a,,,185 asel,a,,,205 asel,a,,,229 asel,a,,,222 asel,a,,,225 asel,a,,,218 asel,a,,,4 asel,a,,,58 asel,a,,,31 asel,a,,,171 asel,a,,,246 asel,a,,,244 asel,a,,,143 asel,a,,,206 asel,a,,,186 asel,a,,,21 asel,a,,,147 asel,a,,,181 asel,a,,,160 asel,a,,,150 asel,a,,,120 asel,a,,,71 asel,a,,,128 asel,a,,,64 asel,a,,,68 asel,a,,,82 asel,a,,,63 asel,a,,,89 asel,a,,,92 asel,a,,,96 asel,a,,,103 asel,a,,,100 asel,a,,,15 asel,a,,,84 asel,a,,,140 asel,a,,,124 asel,a,,,135 asel,a,,,131 asel,a,,,105 asel,a,,,117 asel,a,,,109 asel,a,,,155 asel,a,,,162 asel,a,,,183 asel,a,,,178 asel,a,,,237 asel,a,,,240 asel,a,,,233 asel,a,,,215 asel,a,,,253 asel,a,,,250 asel,a,,,34 asel,a,,,213 asel,a,,,210 asel,a,,,230 asel,a,,,226 asel,a,,,223 asel,a,,,219 nsla,s,1 d,all,uz,0 allsel asel,s,,,7 asel,a,,,9 asel,a,,,11 asel,a,,,13 nsla,s,1 d,all,ux,0,,,uy allsel /solu allsel,all antype,static acel,0,10,0 *do,i,1,100,1 time,i ldread,temp,i,,,,'xianghan','rth' allsel solve allsel esel,s,mat,,i+2 mpchg,i+3,all *enddo finish save |
7楼2017-10-11 16:00:35
8楼2017-10-11 16:02:25
|
finish /clear /filname,xianghan /title,hydration heat of concrete keyw,pr_therm,1 /prep7 !单位，长度：m，力:n cha=0.3$d=0.4$b=6$h=3.6$l=2 hs=0.4$bs=0.15$ts=0.01$r1=0.05$r2=0.025$r3=0.02 k,1 k,2,b+d+d k,3,b+d+d,h+d+d k,4,0,h+d+d k,5,d+cha,d k,6,b+d-cha,d k,7,b+d,d+cha k,8,b+d,d+h-cha k,9,b+d-cha,d+h k,10,d+cha,d+h k,11,d,d+h-cha k,12,d,d+cha *do,i,1,3 l,i,i+1 *enddo l,4,1 *do,i,5,11 l,i,i+1 *enddo l,12,5 al,all vext,1,,,,,l et,1,solid70 *do,i,3,103,1 mp,dens,i,2400 mp,kxx,i,9500 mp,kyy,i,9500 mp,kzz,i,9500 mp,c,i,930 *enddo vsel,all vatt,3,,1 allsel WPROTA,,,90$WPOFF,,,0.2$VSBW,ALL WPOFF,,,0.2$VSBW,ALL WPOFF,,,0.3$VSBW,ALL WPOFF,,,2.7$VSBW,ALL WPOFF,,,2.7$VSBW,ALL WPOFF,,,0.3$VSBW,ALL WPOFF,,,0.2$VSBW,ALL WPOFF,,0.2$WPROTA,,90$VSBW,ALL WPOFF,,,-0.2$VSBW,ALL WPOFF,,,-0.3$VSBW,ALL WPOFF,,,-1.5$VSBW,ALL WPOFF,,,-1.5$VSBW,ALL WPOFF,,,-0.3$VSBW,ALL WPOFF,,,-0.2$VSBW,ALL WPCSYS,-1,0 !网格划分 ESIZE,0.2$VATT,3,1,1$MSHKEY,1 VMESH,ALL *get,NUM,elem,,count NUMMRG,NODE numcmp,node !=========================（二）热分析求解获取温度场=========================! !可以根据材料编号和实常数号选择实体，方便 !初始温度，即浇筑温度 allsel,all initial=20 /solu antype,trans NROPT,FULL !牛顿拉夫逊法 TIMINT,ON !打开时间积分效应 AUTOTS,ON !打开自动时间步长 allsel,ALL,NODE IC,ALL,temp,initial !设置节点初始温度 *dim,hef1,ARRAY,168 hef1(1)=0.511,0.274,0.204,0.175,0.181,0.181,0.183,0.19,0.204,0.23 hef1(11)=0.273,0.349,0.48,0.697,1.026,1.462,1.957,2.415,2.744,2.902 hef1(21)=2.915,2.862,2.833,2.77,2.637,2.46,2.26,2.055,1.859,1.684 hef1(31)=1.53,1.399,1.289,1.195,1.117,1.051,0.995,0.954,0.902,0.863 hef1(41)=0.828,0.795,0.766,0.74,0.714,0.691,0.669,0.648,0.628,0.609 hef1(51)=0.59,0.572,0.554,0.537,0.522,0.503,0.487,0.47,0.454,0.439 hef1(61)=0.424,0.409,0.394,0.38,0.366,0.352,0.339,0.327,0.315,0.303 hef1(71)=0.293,0.282,0.12240,0.12140,0.11789,0.11353,0.10878,0.10423,0.09837,0.09268 hef1(81)=0.08902,0.08561,0.08263,0.07945,0.07643,0.07327,0.06991,0.06688,0.06303,0.05802 hef1(91)=0.02948,0.01151,0.01599,0.04059,0.05279,0.05642,0.05686,0.05619,0.05424,0.05563 hef1(101)=0,05502,0.02433,0.018,0.01278,0.01285,0.00962,0.00600,0.00303,0.00002,0.00206 hef1(111)=0.00375,0.00513,0.00428,0.00460,0.00907,0.01136,0.01209,0.01130,0.01093,0.01103 hef1(121)=0.00964,0.00879,0.00880,0.00910,0.00922,0.00915,0.00267,0.00179,0.00234,0.00206 hef1(131)=0.00147,0.00111,0.00100,0.00129,0.00189,0.00256,0.00374,0.00045,0.00167,0.00336 hef1(141)=0.00418,0.00294,0.01842,0.03276,0.03856,0.04180,0.03701,0.03135,0.02704,0.02379 hef1(151)=0.01992,0.01594,0.01358,0.01168,0.01000,0.00816,0.00654,0.00527,0.00367,0.00200 hef1(161)=0.00011,0.00167,0.00155,0.00328,0.00455,0.00442,0.00215,0.00008 !1h至100h，钢模 *do,h,1,100,1 time,h heat=h+2 !设置瞬态分析步时间 shuihr=hef1(heat)*3600*1350 !水化热生热率，不是总热量，总热量对t求导。W/m3 esel,s,mat,,heat !混凝土材料号从3开始，选择混凝土。 BFE,all,HGEN,,shuihr !在单元上施加水化热放热 !BFE命令为在单元上施加体荷载的命令,HGEN为热分析时体荷载的标识符 allsel asel,s,,,18 asel,a,,,5 asel,a,,,33 asel,a,,,42 asel,a,,,51 asel,a,,,59 asel,a,,,66 asel,a,,,55 !选择预制顶面；空气 nsla,s,1 sf,all,conv,35330,20 allsel asel,u,,,18 asel,u,,,5 asel,u,,,33 asel,u,,,42 asel,u,,,51 asel,u,,,59 asel,u,,,66 asel,u,,,55 !选择现浇侧面；钢模板 nsla,s,1 sf,all,conv,35174,20 allsel,all SOLVE esel,s,mat,,heat !在之前材料号（混凝土号从3开始）的基础上选择混凝土，再改变其属性进入下一轮; mpchg,(h+3),all *enddo !=========================（三）结构分析求解获取应力场=========================! /prep7 allsel,all sfdele,all,all !删除节点上的面荷载 bfdele,all,all !删除面上的体载荷 et,2,solid65 !温度场分析用SOLID70，应力分析改为SOLID185，单元性质改变，定义SOLID185 esel,s,type,,1 emodif,all,type,2 !将混凝土单元改为SOLID185 Ec=3.25e10 *dim,exx,array,103 *do,i,3,103,1 t=(i-2)/24 exx(i)=Ec*exp(0.06*(1-(28/t)**0.5)) *enddo *do,i,3,103,1 mp,dens,i,2400 mp,prxy,i,0.167 mp,alpx,i,0.8e-5 mp,ex,i,exx(i) mp,reft,i,20 *enddo esel,s,mat,,103 mpchg,3,all allsel asel,s,,,54 asel,a,,,65 asel,a,,,48 asel,a,,,35 asel,a,,,29 asel,a,,,22 asel,a,,,3 asel,a,,,17 asel,a,,,26 asel,a,,,47 asel,a,,,39 asel,a,,,46 nsla,s,1 d,all,uy,0 allsel asel,s,,,231 asel,a,,,248 asel,a,,,208 asel,a,,,176 asel,a,,,163 asel,a,,,141 asel,a,,,129 asel,a,,,107 asel,a,,,86 asel,a,,,247 asel,a,,,207 asel,a,,,145 asel,a,,,158 asel,a,,,115 asel,a,,,126 asel,a,,,80 asel,a,,,164 asel,a,,,138 asel,a,,,167 asel,a,,,153 nsla,s,1 d,all,ux,0 allsel asel,s,,,232 asel,a,,,239 asel,a,,,236 asel,a,,,249 asel,a,,,252 asel,a,,,214 asel,a,,,209 asel,a,,,212 asel,a,,,28 asel,a,,,177 asel,a,,,182 asel,a,,,161 asel,a,,,154 asel,a,,,121 asel,a,,,76 asel,a,,,14 asel,a,,,130 asel,a,,,134 asel,a,,,123 asel,a,,,108 asel,a,,,116 asel,a,,,104 asel,a,,,102 asel,a,,,99 asel,a,,,91 asel,a,,,95 asel,a,,,1 asel,a,,,70 asel,a,,,119 asel,a,,,67 asel,a,,,53 asel,a,,,127 asel,a,,,88 asel,a,,,62 asel,a,,,81 asel,a,,,139 asel,a,,,159 asel,a,,,180 asel,a,,,146 asel,a,,,170 asel,a,,,30 asel,a,,,57 asel,a,,,142 asel,a,,,243 asel,a,,,245 asel,a,,,20 asel,a,,,185 asel,a,,,205 asel,a,,,229 asel,a,,,222 asel,a,,,225 asel,a,,,218 asel,a,,,4 asel,a,,,58 asel,a,,,31 asel,a,,,171 asel,a,,,246 asel,a,,,244 asel,a,,,143 asel,a,,,206 asel,a,,,186 asel,a,,,21 asel,a,,,147 asel,a,,,181 asel,a,,,160 asel,a,,,150 asel,a,,,120 asel,a,,,71 asel,a,,,128 asel,a,,,64 asel,a,,,68 asel,a,,,82 asel,a,,,63 asel,a,,,89 asel,a,,,92 asel,a,,,96 asel,a,,,103 asel,a,,,100 asel,a,,,15 asel,a,,,84 asel,a,,,140 asel,a,,,124 asel,a,,,135 asel,a,,,131 asel,a,,,105 asel,a,,,117 asel,a,,,109 asel,a,,,155 asel,a,,,162 asel,a,,,183 asel,a,,,178 asel,a,,,237 asel,a,,,240 asel,a,,,233 asel,a,,,215 asel,a,,,253 asel,a,,,250 asel,a,,,34 asel,a,,,213 asel,a,,,210 asel,a,,,230 asel,a,,,226 asel,a,,,223 asel,a,,,219 nsla,s,1 d,all,uz,0 allsel asel,s,,,7 asel,a,,,9 asel,a,,,11 asel,a,,,13 nsla,s,1 d,all,ux,0,,,uy allsel /solu allsel,all antype,static acel,0,10,0 *do,i,1,100,1 time,i ldread,temp,i,,,,'xianghan','rth' allsel solve allsel esel,s,mat,,i+2 mpchg,i+3,all *enddo finish save |
9楼2017-10-11 16:03:03
