在下写了一个ANSTYS CFD程序，怎么也调不出来，恳请高手指点

wangjite

finish
/clear
rad1=5.89*2
rad2=3.815*2
h1=9.8
h2=h1+3.09
h3=h2+9.65
h4=h3+3.81

/PREP7
ET,1,FLUID141,,,2                             !定义单元类型


K,1 ,rad1,,
K,2 ,rad1,h1,
K,3 ,rad1+2,h1,
K,4 ,rad1+2,h2,
K,5 ,rad1,h2,
K,6 ,rad1,h3,
K,7 ,rad1+1.5,h3,
K,8 ,rad1+1.5,h4,
K,9 ,2.075,h4,
K,10 ,,h3,
K,11 ,,h2,
K,12 ,-,h1,
K,13 ,,,


l,1,2
*repeat,12,1,1
l,13,1
ASKIN,1,12
ASKIN,3,11
ASKIN,5,10
ASKIN,9,7

MSHKEY,1
ESIZE,0.4,0,  
MSHAPE,0,2D
AMESH,1
allsel

MSHKEY,1
ESIZE,0.4,0,  
MSHAPE,0,2D
AMESH,2
allsel

MSHKEY,1
ESIZE,0.4,0,  
MSHAPE,0,2D
AMESH,3
allsel


MSHKEY,1
ESIZE,0.4,0,  
MSHAPE,0,2D
AMESH,4
allsel
ALLSEL,ALL  
save


LSEL,S,,,1
LSEL,A,,,2
LSEL,A,,,4
LSEL,A,,,5
LSEL,A,,,6
LSEL,A,,,8
LSEL,A,,,9
LSEL,A,,,10
LSEL,A,,,11
LSEL,A,,,12
LSEL,A,,,13

DL,ALL,,VX,
DL,ALL,,VY,

LSEL,S,,,3
DL,ALL,,VY,
dl,all,,vx,-5,1

LSEL,S,,,7

wangjite

！接上
DL,ALL,,PRES,,
ALLSEL
FINISH
/SOLU
FLDATA12,PROT,DENS,AIR-IN            ! Set fluid properties
FLDATA12,PROT,VISC,AIR-IN
FLDATA2,ITER,EXEC,40          ! Set FLOTRAN execution controls
FLDATA15,PRES,REFE,14.7
FLDATA14,TEMP,NOMI,70
FLDATA14,TEMP,TTOT,70
FLDATA14,TEMP,BULK,70
TOFF,460

/SOLU
SOLVE                            ! Solve model
FINISH
/POST1
PLNSOL,VX