ansys命令流_ANSYS管道稳态热仿真（材料热物理性能随温度可变）_附命令流

一圆筒形的罐有一接管,罐外径为3英尺，壁厚为0.2英尺，接管外径为0.5英尺，壁厚为0.1英尺，罐与接管的轴线垂直且接管远离罐的端部。如图1所示：

图1 模型示意图

罐内流体温度为华氏450度,与罐壁的对流换热系数为250BUT/hr-ft2-oF,接管内流体的温度为华氏100度,与管壁的对流换热系数随管壁温度而变。接管与罐为同一种材料，它的热物理性能如表1所示，求罐与接管的温度分布。

温度	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

表表1 材料热物理性能

*接管内壁对流系数

以下只给出命令流，菜单操作从略。

/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

3D映射网格图：

温度结果云图：