实例2 一圆筒形的罐有一接管,罐外径为3英尺，壁厚为0.2英尺，接管外径为0.5英尺，壁厚为0.1英尺，罐与接管的轴线垂直且接管远离罐的端部。如图所示：     罐内流体温度为华氏450度,与罐壁的对流换热系数年为250BUT/hr-ft2-oF,接管内流体的温度为华氏100度,与管壁的对流换热系数随管壁温度而变。接管与罐为同一种材料，它的热物理性能如下表所示： 温度 70 200 300 400 500 oF 密度 0.285 0.285 0.285 0.285 0.285 lbm/in3 导热系数 8.35 8.90 9.35 9.8 10.23 Btu/hr-ft-oF 比热 0.113 0.117 0.119 0.122 0.125 Btu/lbm-oF 对流系数* 426 405 352 275 221 Btu/hr-ft2-oF *接管内壁对流系数 求罐与接管的温度分布。 以下分别列出LOG文件及菜单操作 /prep7/title,Steady-state thermal analysis of pipe junction/units,bin       ！使用英制单位et,1,90        ！定义热单元mp,dens,1,.285      ！密度mptemp,,70,200,300,400,500     ！建立温度表mpdata,kxx,1,,8.35/12,8.90/12,9.35/12,9.80/12,10.23/12 ！导热系数mpdata,c,1,,0.133,0.177,0.119,0.122,0.125   ！比热 mpdata,hf,2,,426/144,405/144,352/144,275/144,221/144 ！接管对流系数！定义几何模型参数ri1=1.3        ！罐内半径ro1=1.5       ！罐外半径z1=2        ！罐长ri2=0.4        ！接管内半径ro2=0.5       ！接管外半径z2=2        ！接管长！建立几何模型cylind,ri1,ro1,,z1,,90      ！1/4罐体wprota,0,-90     ！将工作平面旋转到垂直于接管轴线cylind,ri2,ro2,,z2,-90    ！1/4接管wpstyl,defa     ！将工作平面恢复到默认状态vovlap,1,2     ！进行OVERLAP布尔操作/pnum,volu,1     ！打开实体编号/view,,-3,-1,1     ！定义显示角度/type,,4/title, Volumes used in building pipe/tank junctionvplot      ！显示实体vdele,3,4,,1     ！删除多余实体！划分网格asel,,loc,z,z1     ！选择罐上Z=Z1的面 asel,a,loc,y,0     ！添加选择罐上Y=0的面cm,aremote,area    ！创建名为AREMOTE的面组/pnum,area,1/pnum,line,1/title,lines showing the portion being modeledaplot/noeraselplot/eraseaccat,all     ！组合罐远端的面及线，为映射划分网！格作准备lccat,12,7lccat,10,5lesize,20,,,4     ！在接管壁厚方向分4等分lesize,40,,,6     ！在接管长度方向分6等分lesize,6,,,4     ！在罐壁厚方向分4等分allsel      ！选择EVERYTHINGesize,0.4     ！设定默认的单元大小mshape,0,3d     ！选择3D映射网格mshkey,1     save      ！保存数据文件vmesh,all     ！划分网格，产生节点与单元/pnum,defa/title, elements in portion being modeledeplot      ！显示单元finish！加载求解/solu      antype,static     ！定义为稳态分析nropt,auto     ！设置求解选项为Program-chosen ！Newton-Raphsontunif,450     ！设定初始所有节点温度csys,1      ！变为柱坐标nsel,s,loc,x,ri1     ！选择罐内表面的节点sf,all,conv,250/144,450   ！定义对流边界条件cmsel,,aremote    ！选择AREMOTE面组nsla,,1      ！选择属于AREMOTE面组的节点d,all,temp,450     ！定义节点温度wprota,0,-90     ！将工作平面旋转到垂直于接管轴线 cswpla,11,1     ！创建局部柱坐标nsel,s,loc,x,ri2     ！选择接管内壁的节点sf,all,conv,-2,100    ！定义对流边界条件 allsel      ！选择EVERYTHING/pbc,temp,,1     ！显示所有温度约束/psf,conv,,2     ！显示所有对流边界/title,Boundary conditionsnplot      ！显示节点wpstyle,defa     ！工作平面恢复默认状态csys,0      ！变为直角坐标autots,on     ！打开自动步厂长nsubst,50     ！设定子步数量 kbc,0      ！设定为阶越outpr,nsol,last     ！设置输出solve      ！进行求解finish！进入后处理/post1/title,Temperature contrours at pipe/tank junctionplnsol,temp     ！显示温度彩色云图finish/exit,all

评论