[求助]这个程序应该怎么编？？

zzb0888

翅片管热分析，管内外均为水，输入进水的温度、压力、速度，求出水的情况
我的程序如下，不知道哪里错了，结果不对！！
/TITLE,ZZB
L=80E-3
L1=20E-3
LY1=16E-3
LY2=4E-3
LY3=10E-3
LY4=20E-3
LX2=1E-3
LX3=8E-3
/PREP7
toffst,273
SMART,OFF
ET,1,141,,,2      !设定单元类型
KEYOPT,1,1,2      !设定单元类型主要选项
     
RECTNG,0,L,0,LY1   !创建模型
RECTNG,0,L,LY1+LY2,LY1+LY2+LY4
RECTNG,0,L,LY1,LY1+LY2
RECTNG,(LX3-LX2)/2,(LX3+LX2)/2,LY1+LY2,LY1+LY2+LY3
/PNUM,AREA,1     !控制颜色显示  
APLOT            !显示选择的平面
AGEN,10,4,,,LX3  !复制面
NUMCMP,ALL       !压缩编号
/REPLOT         
AOVLAP,2,4       !合并面
NUMCMP,AREA      !压缩编号
*DO,I,4,12,1     !循环合并
AOVLOP,I,13
NUMCMP,AREA
/REPLOT
*ENDDO
AADD,2,3         !迭加面
NUMCMP,AREA
*DO,I,2,10,1     !循环迭加
AADD,2,14-I
NUMCMP,AREA
/REPLOT
*ENDDO
NUMCMP,AREA
RECTNG,-L1,0,0,LY1
RECTNG,-L1,0,LY1+LY2,LY1+LY2+LY4
RECTNG,-L1,0,LY1,LY1+LY2
RECTNG,L,L+L1,0,LY1
RECTNG,L,L+L1,LY1+LY2,LY1+LY2+LY4
RECTNG,L,L+L1,LY1,LY1+LY2
NUMMRG,ALL       !重新编号
NUMCMP,ALL
/REPLOT

ASEL,S,,,2,8,3    !网格
ASEL,A,,,1,7,3    !选面
ESIZE,LX2         !网格大小
TYPE,1
MAT,1
AMESH,ALL
ASEL,S,,,3,9,3
ESIZE,LX2
TYPE,1
MAT,2
AMESH,ALL
ALLSEL
NUMMRG,ALL
NUMCMP,ALL

LSEL,s,,,22     !选线
NSLL,S,1     !选择线上节点
D,ALL,TEMP,20
D,ALL,VX,2.07
D,ALL,VY,0
D,ALL,PRES,0.1E6
LSEL,S,,,59
LSEL,a,,,62
NSLL,S,1
D,ALL,VY,0
LSEL,S,,,37
NSLL,S,1
D,ALL,VX,2
D,ALL,TEMP,60
D,ALL,VY,0
D,ALL,PRES,0.1E6
LSEL,S,,,11
LSEL,A,,,1
LSEL,A,,,58
NSLL,S,1
D,ALL,VY,0
LSEL,S,,,7,33,26
LSEL,A,,,63
NSLL,S,1
D,ALL,VX,0
D,ALL,VY,0
ALLSEL
/REPLOT

MP,DENS,2,8O00
MP,KXX,2,65
MP,C,2,4178E-3

/SOLU
ALLSEL,ALL
FLDA,ITER,EXEC,250
nsel,s,loc,x,-L1,L+L1  
nsel,s,loc,y,0,LY1
ic,all,sp01,1.0   !sp01,02多种自由度
ic,all,sp02,0.0   
asel,s,,,2
asel,a,,,5
asel,a,,,8
nsla,s,1
ic,all,sp01,0.0     
ic,all,sp02,1.0

allsel,all
FLDA,PROT,DENS,CMIX      !定义流体类型为cmix     
FLDA,PROT,VISC,CMIX      
FLDA,PROT,COND,CMIX        
FLDA,PROT,SPHT,CMIX
         
FLDA,VARY,DENS,T       !允许流体变化
FLDA,VARY,VISC,T
FLDA,VARY,COND,T
FLDA,VARY,SPHT,T

FLDA,NOMI,DENS,1.0      !定义相同类型的NOMI系数
FLDA,NOMI,VISC,2.0
FLDA,NOMI,COND,3.0
FLDA,NOMI,SPHT,4.0

msprop,1,spht,constant,4.178
msprop,2,spht,constant,4.183
msprop,1,dens,constant,983.2
msprop,2,dens,constant,998.2
msprop,1,cond,constant,0.6
msprop,2,cond,constant,0.64
msprop,1,visc,constant,47e-5
msprop,2,visc,constant,100.5e-5
FLDA,CONV,temp,1.0E-4   !控制输出
FLDA,TERM,temp,1.0E-3
FLDA,CONV,PRES,1.0E-10   !控制输出
FLDA,TERM,PRES,1.0E-9    !设置控制
FLDA,OUTP,HFLM,T     !
FLDA,MIR,MOME,0.5    !设置松弛因子
FLDA,OUTP,SP01,T
FLDA,OUTP,SP02,T
solve