请教带有表面效应单元的２D温度场分析

zhaoxutao

各位高手好：
我做了个随时间变化的２D温度场的分析，要加载热流密度和热对流，使用表面效应单元，将热流密度加在节点上，热对流及系数加在表面单元上，但是为什么只存在一种，为什么前者被后者给覆盖了呀，那表面效应单元好象没有用，好多书提到要用表面效应单元的，不知道错在那里，小弟不会命令流，所以写不下来，敬请各位高手指点　

fmz

热流密度和对流系数是不可以加载在同一单元上的.不知道你做什么题目,一般想在加热流密度的单元上对流忽略行不行.

zhaoxutao

谢谢fmz的热心帮助，我试试看

zhaoxutao

这是一个路基温度场的模拟，需要施加对流边界和热流密度，
对流及系数、热流是随时间变化的，我是个新手，使用了笨笨
的方法才把命令流写出来，有些错误，尤其是在使用表面
效应单元加载时候，还是出现，覆盖的提示，不知道为什么，
请高手指点或者给个解题的新思路，小弟在此谢谢了
/filname,roadbed
/title,temperature of roadbed
/prep7

et,1,plane55
et,2,surf151 表面效应单元
et,3,plane35

keyopt,2,4,1
keyopt,2,8,1

toffst,273
这里有我有7个温度值,但是只能接受6个,不知道为什么??
mptemp,1,-10,-5,-2,-1,-0.5,15

mpdata,kxx,1,1,2.2,2.2,2.2,2.2,2.2,1.68
mpdata,dens,1,1,1870,1870,1870,1870,1870,1870
mpdata,c,1,1,856.52,856.52,856.52,856.52,856.52,1102.51
mpdata,enth,1,1,0,8.88,15.59,19.53,23.91,102.06

mpdata,kxx,2,1,1.19,1.19,1.19,1.19,1.19,0.95
mpdata,dens,2,1,1708,1708,1708,1708,1708,1708
mpdata,c,2,1,1169.79,1169.79,1169.79,1169.79,1169.79,1439.69
mpdata,enth,2,1,0,16.91,38.24,59.60,70.88,185.46

mpdata,kxx,3,1,0.56,0.56,0.56,0.56,0.56,0.39
mpdata,dens,3,1,1050,1050,1050,1050,1050,1050
mpdata,c,3,1,1465.71,1465.71,1465.71,1465.71,1465.71,2063.14
mpdata,enth,3,1,0,11.20,22.37,31.85,45.48,157.77
关键点
k,,0,0
k,,10,0
k,,13,0
k,,23,0
k,,26,0
k,,36,0
k,,0,7
k,,10,7
k,,13,7
k,,23,7
k,,26,7
k,,36,7
k,,0,10
k,,10,10
k,,13,10
k,,23,10
k,,26,10
k,,36,10
k,,13,12
k,,23,12
创建面
a,1,2,8,7
a,2,3,9,8
a,3,4,10,9
a,4,5,11,10
a,5,6,12,11
a,7,8,14,13
a,8,9,15,14
a,9,10,16,15
a,10,11,17,16
a,11,12,18,17
a,14,15,19
a,15,16,20,19
a,16,17,20
asel,s,,,all
aglue,all
allsel

划分网格
lsel,s,,,4
lsel,a,,,2
lsel,a,,,6
lsel,a,,,9
lsel,a,,,12
lsel,a,,,15
lesize,all,,,14
allsel

lsel,s,,,19
lsel,a,,,17
lsel,a,,,20
lsel,a,,,22
lsel,a,,,24
lsel,a,,,26
lesize,all,,,6
allsel

lsel,s,,,1
lsel,a,,,14
lsel,a,,,3
lsel,a,,,16
lsel,a,,,18
lsel,a,,,27
lesize,all,,,20
allsel

lsel,s,,,5
lsel,a,,,7
lsel,a,,,21
lsel,a,,,11
lsel,a,,,13
lsel,a,,,25
lesize,all,,,10
allsel

lsel,s,,,8
lsel,a,,,10
lsel,a,,,23
lsel,a,,,31
lesize,all,,,30
allsel

lsel,s,,,28
lsel,a,,,30
lesize,all,,,7
allsel

lsel,s,,,29
lsel,a,,,32
lesize,all,,,12
allsel

asel,s,,,1,5
type,1
mat,3
amesh,all
allsel

asel,s,,,6,10
type,1
mat,2
amesh,all
allsel

asel,s,,,12
type,1
mat,1
amesh,all
allsel

asel,s,,,11
asel,a,,,13
type,3
mat,1
amesh,all
allsel

lsel,s,,,18
lsel,a,,,29
lsel,a,,,31
lsel,a,,,32
lsel,a,,,27
nsll,s,1
type,2
esurf,all
allsel

创建数组一部分
*dim,midu,table,5,6,1
*dim,xishu,table,5,5,1
xishu(0,0)=0,1,2,3,4,5
xishu(0,1)=1,864000,1728000,2592000,3456000,4320000
xishu(0,2)=2,10.16958093,9.344549721,9.344549721,9.344549721,8.760563148
xishu(0,3)=3,9.55,8.74,8.74,8.74,8.35
xishu(0,4)=4,9.6,8.69,8.69,8.69,7.99
xishu(0,5)=5,10.21334882,9.371578132,9.371578132,9.371578132,8.856408719
midu(0,0)=0,1,2,3,4,5
midu(0,1)=1,864000,1728000,2592000,3456000,4320000
midu(0,2)=2,5.419129355,4.587292542,3.458110426,2.065910948,0.453013699
midu(0,3)=3,61,43,43,43,18.03
midu(0,4)=4,60,51,51,51,40.42
midu(0,5)=5,55,12,12,12,-8
midu(0,6)=6,129.6,90.43,90.43,90.43,51.18

/solu
antype,trans
trnopt,full
eqslv,jcg 这个选择jacobi共轭梯度迭代的命令也不好使??
outers,all,all
nsubst,1 子步的设定在这里合理还是在*do之后合理呢??
autots,on
kbc,0

初始温度分层施加
nsel,s,loc,y,12,10
ic,all,temp,4
allsel

nsel,s,loc,y,9.5
ic,all,temp,6.11
allsel

nsel,s,loc,y,9
ic,all,temp,2.62
allsel

nsel,s,loc,y,8.5
ic,all,temp,0.55
allsel

nsel,s,loc,y,8
ic,all,temp,-0.4
allsel

nsel,s,loc,y,7.5
ic,all,temp,-1.35
allsel

nsel,s,loc,y,7
ic,all,temp,-1.57
allsel

nsel,s,loc,y,6.5
ic,all,temp,-1.86
allsel

nsel,s,loc,y,6
nsel,a,loc,y,5.5
ic,all,temp,-2.03
allsel

nsel,s,loc,y,5
ic,all,temp,-2.09
allsel

nsel,s,loc,y,4.5
ic,all,temp,-2.12
allsel

nsel,s,loc,y,4
ic,all,temp,-2.1
allsel

nsel,s,loc,y,3.5
ic,all,temp,-2.15
allsel

nsel,s,loc,y,3
nsel,a,loc,y,0
ic,all,temp,-2
allsel

nsel,s,loc,y,2.5
nsel,a,loc,y,2
nsel,a,loc,y,1.5
nsel,a,loc,y,1
nsel,a,loc,y,0.5
ic,all,temp,-2.25
allsel

*do,t,1,5,1

time,t*864000
施加对流条件
lsel,s,,,18
lsel,a,,,27
nsll,s,1
sf,all,conv,%xishu(t,2)%,%midu(t,2)%
allsel

lsel,s,,,29
nsll,s,1
sf,all,conv,%xishu(t,3)%,%midu(t,2)%
allsel

lsel,s,,,32
nsll,s,1
sf,all,conv,%xishu(t,4)%,%midu(t,2)%
allsel

lsel,s,,,31
nsll,s,1
sf,all,conv,%xishu(t,5)%,%midu(t,2)%
allsel

施加热流密度
lsel,s,,,18
lsel,a,,,27
nsll,s,1
ESLN
ESEL,R,TYPE,,2
sfe,all,,hflux,,%midu(t,3)%
allsel

lsel,s,,,29
nsll,s,1
ESLN
ESEL,R,TYPE,,2
sfe,all,,hflux,,%midu(t,4)%
allsel

lsel,s,,,32
nsll,s,1
ESLN
ESEL,R,TYPE,,2
sfe,all,,hflux,,%midu(t,5)%
allsel

lsel,s,,,31
nsll,s,1
ESLN
ESEL,R,TYPE,,2
sfe,all,,hflux,,%midu(t,6)%
allsel

solve
*enddo

aiyouguai

高手在哪里 没发现问题