帮忙看一下我做的水化热应力场命令流对不对？

作者 新星帝国～: 来源: 小木虫 400 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

返回小木虫查看更多

今日热帖