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