【讨论】表面单元加载的问题？

hangch

在学习利用表面单元加载过程中的一个例子,出现如下警告
*** WARNING ***   SUPPRESSED MESSAGE    CP=      7.721   TIME= 16:22:59
  Element 6750 face 1 has convection applied.                              
  This element heat flux will override the convection.      
但是表面单元的功能不就是能够在一个表面上同时施加两种不同的载荷吗,拖过我对整个表面施加热流载荷,那么上表面的对流不就都被覆盖,就不起作用了嘛,不知道问题在哪里?如何解决?请高手指教!   
  
附算例和结果            
/CLEAR,NOSTART                                                !清零
/FILNAM,PLASMASPRAY                                           !制定工作文件
/TITLE,RESIDUAL STRESS ANALYSIS OF THE ORIGINAL PROTOTYPE     !制定工作标题
/UNITS,SI  
/PREP7  
ET,1,SOLID90  
ET,2,Surf152
DOF,TEMP  
KEYOPT,2,3,0  
KEYOPT,2,4,1  
KEYOPT,2,5,0  
KEYOPT,2,6,0  
KEYOPT,2,7,0  
KEYOPT,2,8,4 !根据固体表面温度计算对流系数  
KEYOPT,2,9,0
MP,EX,1,1.2E11  
MP,ALPX,1,7.962E-6  
MP,REFT,1,20  
MP,PRXY,1,0.3  
MP,DENS,1,4000  
MP,C,1,1.5   
MP,KXX,1,0.4  
MP,HF,1.6
BLOCK,0,0.041,0,0.036,0,0.004        
SAVE  
ESIZE,0.001  
VMESH,ALL  
allsel,all
  
asel,s,loc,z,0.004
type,2 !定义属性  
amesh,all  
esel,s,type,,2 !选择表面效应单元  
sf,all,conv,10,25 !确定对流边界条件
  
/PNUM,DEFA  
EPLOT  
FINISH  
  
/SOLU
!开始计算   
  ANTYPE,TRANS  
  TUNIF,150                         !初始温度150度  
  KBC,1  
   
  *SET,DELT,1  
   TIMINT,ON  
   AUTOTS,OFF  
   
  
  TIME,50  
  AUTOTS,on  
  DELTIM,20,10,200
  KBC,0  
  
nsel,s,loc,x,0.015,0.025  
nsel,r,loc,y,0.015,0.021  
nsel,r,loc,z,0.004  
sf,all,hflu,1000
allsel,all
  SOLVE  
  /REPLOT  
      
  
OUTRES,ALL,ALL,  
  
FINISH

hangch

我把热流扩展到整个平面，果然对流被全部覆盖，结果如下,真的很苦恼,这到底是为什么???

xbzhang

加HFLUX时，试试这个：
  
esln
esel,r,type,2
而后再
SF
仅供参考。

hangch

还是有冲突，警告有重叠部分，是不是都必须加载到表面单元上呢

xbzhang

你的KEYOPT有问题，KEYOPT（8）=4,  则include convection (ignore heat flux)，看152帮助.当然加不上热流率了。
我给你改了代码，不知道能否满足你的计算要求。不满足的话提出你的物理模型，要做一个什么计算。一定要根据壁面温度估计CONVECTIVE，很有意思的哦，大家再努力。
我算了10TIME，还有一个图形。
  
/CLEAR,NOSTART !清零  
/FILNAM,PLASMASPRAY !制定工作文件  
/TITLE,RESIDUAL STRESS ANALYSIS OF THE ORIGINAL PROTOTYPE !制定工作标题  
/UNITS,SI  
/PREP7  
ET,1,SOLID90  
ET,2,Surf152  
DOF,TEMP  
KEYOPT,2,3,0  
KEYOPT,2,4,1  
KEYOPT,2,5,0  
KEYOPT,2,6,0  
KEYOPT,2,7,0  
  
KEYOPT,2,8,1
  
KEYOPT,2,9,0  
MP,EX,1,1.2E11  
MP,ALPX,1,7.962E-6  
MP,REFT,1,20  
MP,PRXY,1,0.3  
MP,DENS,1,4000  
MP,C,1,1.5  
MP,KXX,1,0.4  
MP,HF,1.6  
BLOCK,0,0.041,0,0.036,0,0.004  
SAVE  
ESIZE,0.001  
VMESH,ALL  
allsel,all  
  
asel,s,loc,z,0.004  
type,2 !定义属性  
amesh,all  
esel,s,type,,1 !选择90单元。  
sf,all,conv,10,25 !确定对流边界条件  
  
/PNUM,DEFA  
EPLOT  
FINISH  
  
/SOLU  
!开始计算  
  ANTYPE,TRANS  
  TUNIF,150 !初始温度150度  
  KBC,1  
    
  *SET,DELT,1  
   TIMINT,ON  
   AUTOTS,OFF  
    
  
  TIME,10  
  AUTOTS,on  
  DELTIM,20,10,200  
  KBC,0  
  
nsel,s,loc,x,0.015,0.025  
nsel,r,loc,y,0.015,0.021  
nsel,r,loc,z,0.004  
ESLN
ESEL,R,TYPE,,2
sfe,all,,hflux,,1000  
  
allsel,all  
  SOLVE  
  /REPLOT  
  
OUTRES,ALL,ALL,  
  
FINISH

xbzhang

从这个意思上讲，152单元使模型既实现了加载对流，又实现了加载热流率。不知我的理解对不对。

zhaoxutao

太好了，我也遇到了同样的问题，好不让人烦恼啊，我也去试试．
在次谢谢大侠了．．．．

zhaoxutao

esel,s,type,,1 !选择90单元。
sf,all,conv,10,25 !确定对流边界条件

/PNUM,DEFA
EPLOT
．
．

/SOLU
!开始计算
  ANTYPE,TRANS
．
．
nsel,s,loc,x,0.015,0.025
nsel,r,loc,y,0.015,0.021
nsel,r,loc,z,0.004
ESLN
ESEL,R,TYPE,,2
sfe,all,,hflux,,1000　　施加 热流密度

我想知道，确定对流边界条件 在/solu 之前,而施加热流密度在/solu之后,为什么呢,
  这与二者都在/solu之后有什么区别呢?
请大侠们给予指点,谢谢了.　

xbzhang

我自己的想法:
对于静力分析而言，应该是没什么区别，加载可以在前处理中进行，也可以在求解模块中进行。
对于与时间有关的物理过程的有限元分析，载荷常常与时间有关，那么就要在求解模块中了，如果边界条件也与时间有关，那么边界条件也要在求解模块中设定。

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

xbzhang

1.你的命令流里的中文没有加注释符号，这在ANSYS中会出现错误，而且后面的命令会不执行，从而导致错误。以后要养成注释文字在前面加注释符号的习惯。
2.我运行了你的程序，没有发现有错误。
3.你提到的：!这里有我有7个温度值,但是只能接受6个,不知道为什么??
建议马上看看MPTEMP和MPDATA命令的帮助，你会知道为什么，怎么作？下面的命令流供你很快知道怎么作？替你偷懒！

MPTEMP,1,0+TCAP,100+TCAP,300+TCAP,450+TCAP,550+TCAP,600+TCAP
MPTEMP,7,720+TCAP,800+TCAP,1200+TCAP

MPDATA,ALPX,1,,1.0E-5,1.1E-05,1.20E-05,1.3E-05,1.40E-05,1.40E-05      !Coefficient of thermal expansion
MPDATA,ALPX,1,,1.40E-05,1.40E-05,1.50E-05

以下是我改过的命令流：
/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

xbzhang

这是温度云图。供你参考。

zhaoxutao

谢谢你的热心帮助，以后有问题还请多多指点，，

zhaoxutao

xbzhang你好,我运行了你给我修改的"路基温度场"程序,把所有荷载步加全后,运行一段时间后,出现不收敛的现象,我有点找不着头脑了,能不能给我讲讲怎么选择"求解器"中的选项,,小弟在此谢谢了.

xbzhang

我这里没有出现不收敛啊，我用的是5.7。而且错误文件里一个警告都没有。
一定要选择求解器吗，默认就可以了。

lion2004

不错不错， 研究一下

xbzhang

把你现在的程序贴出来吧。

zhaoxutao

(1) 首先象各位道歉，上次我给的命令流有些错误,就是焓值错了,这次改过来了.
不收敛的情况得到纠正了,我是把非线性求解中的"线性搜索"打开了,还有把
收敛公差修改了。结果出来了，比上次的好多了，这样看着才象路基温度场，上次的
太不象话了

(2) 但是我不知道，修改公差对结果有多大的影响??应该修改到什么程度才是最好的，
有没有计算公式或者参考数值???

(3)在施加边界条件是时候，xbzhang  运行说没有错误提示，但是在我的机子上还是出现如下错误提示：

*** WARNING ***                         CP =      12.875   TIME= 12:49:30
Element 6264 face 1 has convection applied.                             
  This element heat flux will override the convection.   
               
但是查看加载的数据的时候，对流和热流密度
都存在，是程序版本的原因,又不太可能啊,很想知道为什么???????

(4) 还有个问题就是，我的命令很不熟练，在输入的时候，打开菜单验证，但是有的
命令好象没有起作用。比如：定义荷载时间，自动时间步和子步数
time,3   autots,on    nsubst,50
但是在打开菜单的时候，只有前两项，子步数没有显示??????

很想知道为什么，是不是在运行的时间会有呢??

(5)以下是命令流，请各位高手给予指点，先谢谢了。

/filname,roadbed
/title,temperature of roadbed 2m
/prep7

et,1,plane55
et,2,surf151
et,3,plane35

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

toffst,273
mptemp,1,-10,-5,-2,-1,-0.5,0
mptemp,7,15

mpdata,kxx,1,1,2.2,2.2,2.2,2.2,2.2,1.68
mpdata,kxx,1,7,1.68
mpdata,dens,1,1,1870,1870,1870,1870,1870,1870
mpdata,dens,1,7,1870
mpdata,c,1,1,856.52,856.52,856.52,856.52,856.52,1102.51
mpdata,c,1,7,1102.51
mpdata,enth,1,1,0,8880000,15590000,19530000,23910000,71130000
mpdata,enth,1,7,102060000

mpdata,kxx,2,1,1.19,1.19,1.19,1.19,1.19,0.95
mpdata,kxx,2,7,0.95
mpdata,dens,2,1,1708,1708,1708,1708,1708,1708
mpdata,dens,2,7,1708
mpdata,c,2,1,1169.79,1169.79,1169.79,1169.79,1169.79,1439.69
mpdata,c,2,7,1439.69
mpdata,enth,2,1,0,16910000,38240000,59600000,70880000,148570000
mpdata,enth,2,7,185460000

mpdata,kxx,3,1,0.56,0.56,0.56,0.56,0.56,0.39
mpdata,kxx,3,7,0.39
mpdata,dens,3,1,1050,1050,1050,1050,1050,1050
mpdata,dens,3,7,1050
mpdata,c,3,1,1465.71,1465.71,1465.71,1465.71,1465.71,2063.14
mpdata,c,3,7,2063.14
mpdata,enth,3,1,0,11200000,22370000,31850000,45480000,125280000
mpdata,enth,3,7,157770000

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,,,28
allsel

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

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

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

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

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

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

meshkey,2
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,20,6,1
*dim,xishu,table,20,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

xishu(6,0)=6,7,8,9,10
xishu(6,1)=5184000,6048000,6912000,7776000,8640000
xishu(6,2)=8.760563148, 8.760563148,8.522678851,8.522678851,8.522678851
xishu(6,3)=8.35,8.35,8.12,8.12,8.12
xishu(6,4)=7.99,7.99,7.35,7.35,7.35
xishu(6,5)=8.856408719,8.856408719,8.546581297,8.546581297,8.546581297

xishu(11,0)=11,12,13,14,15
xishu(11,1)=9504000,10368000,11232000,12096000,12960000
xishu(11,2)= 9.005836281,9.005836281,9.005836281,10.26664366,10.26664366
xishu(11,3)=9.63,9.63,9.63,9.79,9.79
xishu(11,4)=8.27,8.27,8.27,8.47,8.47
xishu(11,5)=10.20317293,10.20317293,10.20317293,10.34177396,10.34177396

xishu(16,0)=16,17,18,19,20
xishu(16,1)=13824000,14688000,15552000,16416000,17280000
xishu(16,2)= 10.26664366,11.38593646,11.38593646,11.38593646,12.39516568
xishu(16,3)=9.79,10.85,10.85,10.85,12.92
xishu(16,4)=8.47,9.56,9.56,9.56,11.29
xishu(16,5)=10.34177396,11.47887248,11.47887248,11.47887248,13.61611889

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

midu(6,0)=6,7,8,9,10
midu(6,1)=5184000,6048000,6912000,7776000, 8640000
midu(6,2)=-1.331555942,-3.233556414,-5.195178082,-7.15679975,-9.058800222
midu(6,3)=18.03,18.03,-3.78,-3.78,-3.78
midu(6,4)=40.42,40.42,21.99,21.99,21.99  
midu(6,5)=-8,-8,-26,-26,-26
midu(6,6)=51.18,51.18,19.85,19.85,19.85

midu(11,0)=11,12,13,14,15
midu(11,1)=9504000,10368000,11232000,12096000,12960000
midu(11,2)=-10.84336986,-12.45626711,-13.84846659,-14.97764871,-15.80948552
midu(11,3)=-12,-12,-12,1.814,1.814
midu(11,4)=13.87,13.87,13.87,22.84,22.84
midu(11,5)=-30,-30,-30,-18,-18
midu(11,6)=8.58,8.58,8.58, 30.68,30.68

midu(16,0)=16,17,18,19,20
midu(16,1)=13824000,14688000,15552000,16416000,17280000
midu(16,2)=-16.31868377,-16.48975342,-16.31747829,-15.80707456,-14.97403227
midu(16,3)=1.814,33.91,33.91,33.91,34
midu(16,4)=22.84,51.43,51.43,51.43,41.29
midu(16,5)=-18,0,0,0,3
midu(16,6)=30.68,63.73,63.73,63.73,62.85

!施加荷载
/solu
antype,trans
trnopt,full
eqslv,jcg       !菜单没有显示?????
outres,all,all
nsubst,2        !菜单没有显示?????
autots,on       !菜单没有显示?????
kbc,0
lnsrch,on          !线性搜索
cnvtol,heat,,0.1  !修改公差

!初始温度
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,20,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

zhaoxutao

怎么没有人回答我的问题呢,是不是有点幼稚 啊