急问,我程序哪错了啊???

liyan851205

大概是这样一个过程,4种不同材料组成的钢包各层,包内有气体(气体模型不会建),与内壁对流,包内一点温度高(NN+2),向包内辐射,空气与外包壁对流,外包壁向空气辐射.空气环境中设一点(NN+1)
ET,1,SOLID70
et,2,surf152,,,,1,1  !表面效应单元
KEYOPT,2,1,1  !表面效应单元设置
KEYOPT,2,4,1  
KEYOPT,2,5,1  
KEYOPT,2,8,5
KEYOPT,2,9,1

!!!!!!包壁特性!!!!!!
MPTEMP,1,20,400,800,1200
MPDATA,KXX,1,1,50,39,30,18
MPDATA,C,1,1,400,420,510,600
MP,DENS,1,7800
!!!!!绝热层!!!!!!
MP,KXX,2,0.1
MP,DENS,2,200
MP,C,2,800
!!!!!!!永久层!!!!!!!!!
MPTEMP,1,20,400,800,1200
MPDATA,C,3,1,610,750,1175,1320
MP,KXX,3,0.5
MP,DENS,3,2500
!!!!!!工作层!!!!!!!!
MP,KXX,4,1.15
MPTEMP,1,20,400,800,1200
MP,C,4,1,800,900,1020,1200
MP,DENS,4,3000
mp,emis,2,0.8
!!!!!!钢包几何模型!!!!!!!
/TITLE,THE PHYSICAL MODLE OF LADLE
cyl4,,,1.2995,0,1.4595,90,3.290
cyl4,,,1.4595,0,1.5195,90,3.570
cyl4,,,1.5195,0,1.5395,90,3.750
cyl4,,,1.5395,0,1.5645,90,3.780
cylind,1.4595,,3.290,3.570,0,90
cylind,1.5195,,3.570,3.750,0,90
cylind,1.5395,,3.750,3.780,0,90
cylind,1.5645,,3.780,3.810,0,90
/PNUM,VOLU,1
VADD,1,5
VADD,2,6
VADD,3,7
VADD,4,8

!!!!!划分网格!!!!!
/TITLE,MODLE AFTER MESH
MSHAPE,1,3D
MSHKEY,0
ESIZE,0.2
MAT,1
VMESH,3
MAT,2
VMESH,2
MAT,3
VMESH,1
MAT,4
VMESH,9
/VIEW,  1,  0.437130844488    , -0.757961482320    , -0.484160114132
/ANG,   1,   66.4281439461  
/REPLO

*GET,NN,NODE,,num,max
N,NN+1,0,2.5,1.6
TYPE,2                                                 !指定节点单元类型为2                  
!Real,2                                                !
Mat,2  
ASEL,S,,,21                                                !                                                                         !选定面21
NSLE,S,1                                               !选定面21上的所有节点
ESURF,NN+1                                             !对面21表面加对节点NN+1的表面效应单元
ASEL,S,,,33                                            !选定面33
NSLE,S,1                                               !选定面33上的所有节点
ESURF,NN+1                                             !对面33表面加对节点NN+1的表面效应单元
ASEL,S,,,41                                            !选定面41
NSLE,S,1                                               !选定面41上的所有节点
ESURF,NN+1
                                                
N,NN+2,0,0,3.810
TYPE,2                                                 !指定节点单元类型为2                  
!Real,2                                                !
Mat,2
ASEL,S,,,4
NSLE,S,1
ESURF,NN+2
ASEL,S,,,45
NSLE,S,1
ESURF,NN+2
                                       
SAVE
FINISH
/SOL
solcontrol,1
antype,trans
timint,off
time,0.01
deltim,0.01
esel,s,mat,,1
nsle,s
d,all,temp,100
esel,s,mat,,2
nsle,s
d,all,temp,300
esel,s,mat,,3
nsle,s
d,all,temp,500
esel,s,mat,,4
nsle,s
d,all,temp,800
tunif,100
allsel
solve
eplot
timint,on
time,3600
nsubst,36
kbc,1     
autos,on
ddele,all,temp          !!!删除稳态分析中定义的结点温度!!!
OUTRES,,ALL
TOFFST,273
D,2726,,Tair,,,,temp
D,2727,,Tgas,,,,TEMP
ASEL,S,,,21
NSLA,S,1                                             
SF,ALL,CONV,5.9,35
ASEL,S,,,33
NSLA,S,1                                                                     !
SF,ALL,CONV,5.9,35                                  !
ASEL,S,,,4                                             
NSLA,S,1                                             
SF,ALL,CONV,30,1300                                 !施加包壁内侧的对流边界条件
ASEL,S,,,45                                            !选定面45
NSLA,S,1                                             !选定所有节点
SF,ALL,CONV,30,1300                                 !施加包底内侧的对流边界条件
ASEL,S,,,41                                            !选定面41
NSLA,S,1                                               !选定所有节点
SF,ALL,CONV,5.9,35                                 !施加包底外侧的对流边界条件
/psf,conv,2
                 
SOLVE

heliming0802

发现一个问题,你搞的问题和我在做的差不多.你还是有头绪了,我都不知道如何下手.