请教熟悉命令的高手：把一个管内传热的例子转换为ＧＵＩ形式
劳驾高手把下面的例子分析一下，并用ＧＵＩ形式展出供大家共同学习！谢谢，一起努力求解圆管中流体以及管壁的温度分布，考虑流体与管内壁的对流，管外壁与
周围环境对流。已知条件如下：
管内流体的物理属性以及入口条件
导热系数 0.6W/(m*c)
密度 1000kg/m3
比热 4183J/(kg*c)
入口温度 20c
入口流速 3m/s
壁管的物理属性和几何属性
导热系数 100W/(m*c)
管外径 5/16inch
壁厚 0.083inch
长度 12inch
管外对流条件
环境温度 100c
对流系数 10000w/(m2*c)
管内对流条件
对流系数 15000w/(m2*c)
/prep7
!转换数据为SI
pi=acos(-1)
ri=((5/16)-0.083)*0.0254 =0.0058293 !管内半径
r0=(5/16)*0.0254=0.0079375 !管外半径
l=12*0.0254 =0.3048 !管长度
nthick=1 !径向单元个数
size=(r0-ri)/nthick =0.0021082 !径向单元长度
nlength=nint(l/size) =145 !长度方向单元个数
ktube=100 !管导热系数
kwater=0.6 !流体导热系数
denswatr=1000 !流体密度
cwater=4183 !流体比热
vinput=3 !入口流速
carea=pi*(ri**2) !流通面积
massflow=vinput*carea*denswatr !质量流率
tinlet=20 !流体入口温度
tbulk=100 !管外环境温度
hi=15000 !管内壁与流体的对流系数
h0=10000 !管外壁对流系数
!定义三类单元
et,1,116,1,1 !KEYOPT(1)＝1自由度为温度
!KEYOPT(2)＝1使用SURF151模拟管内对流
et,2,151,,,1,1,1 !KEYOPT(3)＝1，轴对称
!KEYOPT(4)＝1，无中间节点
!KEYOPT(5)＝1，有额外节点，在FLUID116上
keyopt,2,8,2 !KEYOPT(8)＝2，考虑对流
et,3,55,,,1 !KEYOPT(3)＝1，轴对称
r,1,2*ri,carea,1 !FLUID116的实常数
mp,kxx,1,kwater
mp,c,1,cwater
mp,dens,1,denswatr
mp,kxx,3,ktube
!建模
k,1,0,0,0
k,2,0,-l,0 !注意不是-1，而是-L
k,3,0,-(l+size),0 !注意不是-(1+size)，而是-(L+size)
l,1,2
l,2,3 !线1和2模拟管内流体流动
cm,waterlin,lines !创建组元，用于生成表面单元
rectng,ri,r0,0,-l !管壁,-L
type,1
real,1
mat,1
esize,,nlength
lmesh,1
esize,,1
lmesh,2
lsel,s,loc,x,(ri+r0)/2 !选择管两端
lesize,all,,,nthick !厚度方向单元个数
esize,,nlength !单元大小
type,3
mat,3
allsel,all
amesh,1
lsel,s,loc,x,ri !选择管内壁
cm,tubelin,lines
type,2
lfsurf,'tubelin','waterlin' !使用宏LFSURF创建表面效应单元
/solu
dk,1,temp,tinlet !流体入口温度
lsel,s,loc,x,r0
sfl,all,conv,h0,,tbulk !管外壁对流边界
esel,s,ename,,surf151
sfe,all,,conv,0,hi !管内壁对流系数
esel,s,ename,,fluid116
sfe,all,,hflux,,massflow !管内流体质量流率，以热流密度方式加载
allsel,all
solve
finish
/post1
nsel,s,loc,x,0 !选择FLUID116
prnsol,temp !列出流体温度分布
nsel,s,loc,x,r0 !选择管外壁节点
prnsol,temp !列出管外壁节点温度 终于已经把该命令搞懂了．不过有出了个问题．我模拟出来的图片竟然整个管壁的温度都相同，没有温度剃度，为什么？请高手指正 这个命令流可以用不？和我 想做的很相视啊，不过不怎么懂！还请楼主多指点，另外您说的问题解决了吗？
[ 本帖最后由 fengkai 于 2006-10-31 10:13 编辑 ]
页:
[1]