如何把这个例题转变成瞬态的？

cuizhikun66

该题已经做过，但是它是稳态的．请高手指点怎么把他变为瞬态的？
求解圆管中流体以及管壁的温度分布，考虑流体与管内壁的对流，管外壁与
周围环境对流。已知条件如下：
管内流体的物理属性以及入口条件
导热系数          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                !流1
  
vinput=3                   !入口流速
carea=pi*(ri**2)=0.0001068           !流通面积
massflow=vinput*carea*denswatr  =0.3204 !质量流率
  
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                !列出管外壁节点温度



谢谢各位指点！

[ 本帖最后由 xinyuxf 于 2007-7-12 14:48 编辑 ]

cuizhikun66

我试做一次，但是结果和稳态没区别．不知道错在哪？