空心转子热流分析命令流运行不出来

漩涡迷人

空心转子的热流分析，空腔内流体为空气，隔油套与转子间隙内流体也为空气。
计算区域为：转子，流体，隔油套。
求整个温度场分布，以及流体的流速。
命令流不完全正确，但一步步复制到命令窗口也可建模，网格划分，但最后结果运行不出来。
请各位高手帮我看看是不是边界加载不正确，非常感谢！
/clear,all
/TITLE,Thermal analysis of rotor after stopped
keyw,flotran CFD,1
/units,SI
/PREP7
!build model(2D)
k,1,,10e-3
k,,51e-3,10e-3
k,,55e-3,10e-3
k,,71e-3,10e-3
k,,126e-3,10e-3
k,,141e-3,10e-3
k,,396e-3,10e-3
k,,411e-3,10e-3
k,,432e-3,10e-3
k,,480e-3,10e-3
k,,496e-3,10e-3
k,,512e-3,10e-3
k,13,512e-3,26e-3
k,14,496e-3,26e-3
k,15,496e-3,21e-3
k,16,480e-3,21e-3
k,17,480e-3,26e-3
k,18,432e-3,26e-3
k,19,432e-3,32e-3
k,20,411e-3,32e-3
k,21,411e-3,34e-3
k,22,396e-3,34e-3
k,23,396e-3,27e-3
k,24,141e-3,27e-3
k,25,141e-3,34e-3
k,26,126e-3,34e-3
k,27,126e-3,27e-3
k,28,71e-3,27e-3
k,29,71e-3,23e-3
k,30,55e-3,23e-3
k,31,55e-3,27e-3
k,32,51e-3,27e-3
k,33,51e-3,19e-3
k,34,0,19e-3
k,35,0,0
k,36,512e-3,0
*do,i,1,33
l,i,i+1
*enddo
l,34,1
lplot
l,2,33
l,3,30
l,4,29
l,5,27
l,6,24
l,7,23
l,8,20
l,9,18
l,10,16
l,11,15
l,28,31
l,25,22
l,17,14
/pnum,line,1
al,34,1,35,33
al,35,32,31,30,36,2
al,30,45,28,29
al,36,29,37,3
al,37,28,27,38,4
al,38,26,25,24,39,5
al,24,46,22,23
al,39,23,40,6
al,40,22,21,20,41,7
al,41,19,18,42,8
al,42,17,16,43,9
al,16,47,14,15
al,43,15,44,10
al,14,13,12,11,44
k,37,126e-3,40e-3
k,38,411e-3,40e-3
k,39,411e-3,41e-3
k,40,411e-3,42e-3
k,41,126e-3,42e-3
k,42,126e-3,41e-3
l,26,37
l,37,38
l,38,21
l,39,40
l,40,41
l,41,42
l,39,42
al,48,49,50,21,46,25
al,51,52,53,54
l,38,39
l,37,42
al,49,55,54,56
l,12,36
l,35,36
l,1,35
al,1,2,3,4,5,6,7,8,9,10,11,57,58,59
! 2D mesh
allsel,all
! input et mp
et,1,fluid141
!et,2,surf151
!keyopt,2,8,2
mp,kxx,1,0.0296
mp,kxx,2,49.8
mp,DENS,2,7840
mp,C,2,465
mp,kxx,3,10        !kxx_NdFeB
mp,kxx,4,11.4      !护套 kxx_gh
mp,kxx,5,0.25      !油套
type,1
lsel,s,,,30,36,6
lccat,all
lsel,s,,,32,35,3
lccat,all
lsel,s,,,26,38,12
lccat,all
lsel,s,,,24,39,15
lccat,all
lsel,s,,,22,40,18
lccat,all
lsel,s,,,28,37,9
lccat,all
lsel,s,,,20,41,21
lccat,all
lsel,s,,,18,42,24
lccat,all
lsel,s,,,16,43,27
lccat,all
lsel,s,,,14,44,30
lccat,all
lsel,s,,,21,25,4
lsel,a,,,46
lccat,all
lsel,s,,,1,11
lccat,all
aglue,all
allsel,all
asel,s,,,15
aatt,4,,1
allsel,all
asel,s,,,7
aatt,3,,1
allsel,all
asel,s,,,17,18
aatt,1,,1
allsel,all
asel,s,,,16
aatt,5,,1
allsel,all
asel,s,,,7,15,8
asel,a,,,16,18
asel,inve
aatt,2,,1
size1=1e-3
size2=5e-3
lsel,s,,,48,50,2
lsel,a,,,57,59,2
lesize,all,size1
lsel,inve
lesize,all,size2
allsel,all
mshkey,2
mshape,0,2d
amesh,all
allsel,all
ldele,51,61
finish
/solo
!施加合金护套,永磁体，轴芯热源
allsel,all
asel,s,,,15
bfa,all,hgen,5.5222e+006
allsel,all
asel,s,,,7
bfa,all,hgen,1.1475e+006
allsel,all
asel,s,,,7,15,8
asel,a,,,16,18
asel,inve
bfa,all,hgen,2.2614e+005
!施加间隙内气体温度，压强边界条件
allsel,all
LSEL,S,,,56
DL,ALL,,TEMP,30
dl,all,,pres,0.22e+6
lsel,s,,,55
dl,all,,pres,1.03e+5

!施加空腔内气体温度，压强边界条件
allsel,all
LSEL,S,,,59
DL,ALL,,TEMP,25
dl,all,,pres,0.22e+6
lsel,s,,,57
dl,all,,pres,1.03e+5
!施加隔油套外壁的温度边界条件
LSEL,S,,,52
DL,ALL,,TEMP,50
!施加模型外各个端部的对流换热边界条件
lsel,s,,,34
sfl,all,conv,183.3,,30,
lsel,s,,,32
sfl,all,conv,223.5,,30,
lsel,s,,,26,48,22
sfl,all,conv,270.9,,30,
lsel,s,,,20,50,30
sfl,all,conv,272.6,,30,
lsel,s,,,18
sfl,all,conv,247.6,,30,
lsel,s,,,12
sfl,all,conv,201,,30,
lsel,s,,,51,53,2
sfl,all,conv,10,,30,
!施加模型外各个外表面的对流换热边界条件
lsel,s,,,33
sfl,all,conv,295.7,,30,
lsel,s,,,27,31,4
lsel,a,,,45
sfl,all,conv,332.5,,30,
lsel,s,,,19
sfl,all,conv,351.9,,30,
lsel,s,,,13,17,4
lsel,a,,,47
sfl,all,conv,328.3,,30,

allsel,all
/solo
fldata1,solu,flow,1
fldata1,solu,temp,1
fldata1,solu,turb,1
fldata2,ITER,EXEC,100
fldata7,PROT,DENS,CONSTANT
fldata8,NOMI,DENS,1.029
fldata7,PROT,VISC,CONSTANT
fldata8,NOMI,VISC,20.02E-6
fldata7,PROT,COND,CONSTANT
fldata8,NOMI,COND,0.0296
fldata7,PROT,SPHT,CONSTANT
fldata8,NOMI,SPHT,1.009
fldata24,TURB,MODL,2
fldata24,TURB,RATT,2
fldata18,METH,TEMP,6
SOLVE
FINISH
/POST1
SET,LAST
PLNSOL,TEMP,,0,
PLVECT,V,,,,VECT,ELEM,ON,0
FINISH