我做了一个有关流固偶合的例子，中间有流场和电热场，贴出来跟大家讨论一下

我在做一个有关流固偶合的问题。现将问题和程序付上。
有一个房间，在地板上铺导电的混凝土，混凝土通电后就会发热，通过对流辐射向空气传热，求通电10000秒后室内空气的温度分布。房间是用保温材料做成的。
程序做完后有一个问题就是结果不收敛。望有人可以看下我的问题出在哪？

/prep7
et,1,141
keyopt,1,4,1
et,2,PLANE67 ！定义单元类型电热单元
ET,3,PLANE55
RECTNG,0.3,1.5,0.32,2.3
type,1
LESIZE,4,,,50
LESIZE,2,,,50
LESIZE,1,,,30
LESIZE,3,,,30
ALLSEL
AMESH,ALL
ALLSEL,BELOW,AREA
NSEL,R,LOC,y,0.32
SF,ALL,FSIN,1 ！设置耦合标签
ALLSEL,all
lsel,s,,,2
NSLL,S,1
SF,ALL,FSIN,2
ALLSEL,all
lsel,s,,,3
NSLL,S,1
SF,ALL,FSIN,3
ALLSEL,all
lsel,s,,,4
NSLL,S,1
SF,ALL,FSIN,4

ALLSEL,all
lsel,s,,,1,4,1
NSLL,S,1
SF,ALL,CONV,12,40 !加对流

ALLSEL,all !设置无滑移边界条件
lsel,s,,,2,4,1
Dl,ALL,1,vx,0
Dl,ALL,1,vy,0
Dl,ALL,1,ux,0
Dl,ALL,1,uy,0

ALLSEL
RECTNG,0,0.3,0,2.6
RECTNG,1.5,1.8,0,2.6
RECTNG,0.3,1.5,2.3,2.6
RECTNG,0.3,1.5,0,0.3
AADD,2,3,4,5

/PREP7
FLST,2,3,4,ORDE,3
FITEM,2,5
FITEM,2,9
FITEM,2,21
LCOMB,P51X, ,0

/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,1.08222638492,1
/REP,FAST
FLST,2,3,4,ORDE,3
FITEM,2,7
FITEM,2,11
FITEM,2,22
LCOMB,P51X, ,0

RECTNG,0.3,1.5,0.3,0.32 !正确
MPTEMP,1,0
MPDATA,KXX,2,,1.28
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,2,,2147
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,C,2,,1210
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,RSVX,2,,0.61
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,3,,0.028
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,3,,38
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,C,3,,1413

CM,_Y,AREA
ASEL, , , , 2
CM,_Y1,AREA
CMSEL,S,_Y
!*
CMSEL,S,_Y1
AATT, 2, , 2, 0,
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
CM,_Y,AREA
ASEL, , , , 6
CM,_Y1,AREA
CMSEL,S,_Y
!*
CMSEL,S,_Y1
AATT, 3, , 3, 0,
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1

LESIZE,5,,,34 ！划分网格
LESIZE,7,,,34
LESIZE,13,,,30
LESIZE,11,,,30
LESIZE,8,,,54
LESIZE,10,,,54
LESIZE,23,,,50
LESIZE,24,,,50
LESIZE,9,,,5
LESIZE,19,,,30

ALLSEL
AMESH,ALL
ALLSEL
LSEL,S, , , 11 ！设置耦合标签
NSLL,S,1
SF,ALL,FSIN,1
ALLSEL
LSEL,S, , , 24 ！设置耦合标签
NSLL,S,1
SF,ALL,FSIN,2
ALLSEL
LSEL,S, , , 13 ！设置耦合标签
NSLL,S,1
SF,ALL,FSIN,3
ALLSEL
LSEL,S, , , 23 ！设置耦合标签
NSLL,S,1
SF,ALL,FSIN,4

ALLSEL,ALL
FLST,2,1,4,ORDE,1
FITEM,2,12 ！选择线12施加电压为0
/GO
!*
DL,P51X, ,VOLT,
FLST,2,1,4,ORDE,1
FITEM,2,9
/GO
!*
DL,P51X, ,VOLT,50 ！施加电压

ALLSEL,ALL
lsel,s,,,5
lsel,A,,,10
lsel,A,,,7
lsel,A,,,8
lsel,A,,,19
lsel,A,,,23
lsel,A,,,13
lsel,A,,,24
NSLL,S,1
dl,all,6,VX,0 !定义无滑移边界条件
dl,all,6,VY,0
dl,all,6,uX,0 !定义无滑移边界条件
dl,all,6,uY,0
ALLSEL,ALL
lsel,s,,,6
lsel,A,,,9
lsel,A,,,11
lsel,A,,,12
NSLL,S,1
dl,all,2,VX,0 !定义无滑移边界条件
dl,all,2,VY,0
dl,all,2,uX,0 !定义无滑移边界条件
dl,all,2,uY,0

ALLSEL,all
lsel,s,,,11
NSLL,S,1
SF,ALL,CONV,12,40 ！加对流换热系数
FINISH
ALLSEL,ALL
/AUX12 ！加辐射
!*
STEF,0.567E-07
TOFFST,273
RADOPT,0.10000000149,0.1E-02,0,1000,0.10000000149,0.10000000149
SPCTEMP,1,40
!*
FINISH
/SOL
FLST,5,1,4,ORDE,1
FITEM,5,11
CM,_Y,LINE
LSEL,R, , ,P51X
CM,_Y1,LINE
CMSEL,S,_Y
CMDELE,_Y
!*
STEF,5.67e-008,
!*
SFL,_Y1,RAD,0.8, ,15
!*
CMDELE,_Y1
!*
FLST,5,3,4,ORDE,3
FITEM,5,3
FITEM,5,-4
FITEM,5,24
CM,_Y,LINE
LSEL,R, , ,P51X
CM,_Y1,LINE
CMSEL,S,_Y
CMDELE,_Y
!*
STEF,5.67e-008,
!*
SFL,_Y1,RAD,0.6, ,15
!*
CMDELE,_Y1

/SOLU ！求解控制22
FLDATA1,SOLU,TRAN,1
FLDATA1,SOLU,FLOW,1
FLDATA1,SOLU,TEMP,1
FLDATA1,SOLU,TURB,0
FLDATA1,SOLU,COMP,0
FLDATA1,SOLU,VOF,0
FLDATA1,SOLU,SFTS,0
FLDATA1,SOLU,IVSH,0
FLDATA1,SOLU,SWRL,0
FLDATA1,SOLU,SPEC,0
FLDATA1,SOLU,ALE,1
FLDATA1,SOLU,RDSF,1

!*
FLDATA,ALGR,SEGR,SIMPLEF
FLDATA,ALGR,HFLM,TEMP
!*
*SET,_z4,-1
*SET,_z14,1
!*
FLDATA12,PROP,DENS,4
FLDATA13,VARY,DENS,1
FLDATA12,PROP,VISC,4
FLDATA13,VARY,VISC,1
FLDATA12,PROP,COND,4
FLDATA13,VARY,COND,0
FLDATA12,PROP,SPHT,4
FLDATA13,VARY,SPHT,0
!*
FLDATA7,PROT,DENS,AIR-SI!空气密度
FLDATA8,NOMI,DENS,-1
FLDATA9,COF1,DENS,0
FLDATA10,COF2,DENS,0
FLDATA11,COF3,DENS,0
FLDATA7,PROT,VISC,AIR-SI!空气速度
FLDATA8,NOMI,VISC,-1
FLDATA9,COF1,VISC,0
FLDATA10,COF2,VISC,0
FLDATA11,COF3,VISC,0
FLDATA12,PROP,IVIS
FLDATA7,PROT,COND,AIR-SI
FLDATA8,NOMI,COND,-1
FLDATA9,COF1,COND,0
FLDATA10,COF2,COND,0
FLDATA11,COF3,COND,0
FLDATA7,PROT,SPHT,AIR-SI
FLDATA8,NOMI,SPHT,-1
FLDATA9,COF1,SPHT,0
FLDATA10,COF2,SPHT,0
FLDATA11,COF3,SPHT,0
!*
FLDATA15,PRES,REFE,101350,
FLDATA16,BULK,BETA,1e+015,
FLDATA17,GAMM,COMP,1.4,
FLDATA14,TEMP,NOMI,40,
FLDATA14,TEMP,TTOT,40,
FLDATA14,TEMP,BULK,40, !体积温度
TOFFSET,273, !温度偏移
!*
ACEL,0,9.8,0,
FLDATA26,STAB,MOME,1e-6, !修改过
FLDATA26,STAB,PRES,1e-6,
FLDATA26,STAB,TEMP,1e-6,
FLDATA26,STAB,TURB,1e-6,
FLDATA26,STAB,VISC,1e-6,
!*
FLDATA40,WADV,TEMP,1
FLDATA25,RELX,VX,0.5,
FLDATA25,RELX,VY,0.5,
FLDATA25,RELX,VZ,0.5,
FLDATA25,RELX,PRES,0.5,
FLDATA25,RELX,TEMP,0.8,
FLDATA25,RELX,ENKE,0.5,
FLDATA25,RELX,ENDS,0.5,
!*
!*
FSAN,1
FSOR,SOLID
FSTR,FLUID,TRAN
FSTR,SOLID,TRAN

!*
FSTI,10000
FSDT,100
fsre,all,1 !松弛系数
!*
FSIT,100 !交错迭代最大次数
!*
FSCO,ALL,1.0E-5, !指定收敛值
FSOU,1
fsin,cons !保守传递

DELTIM,100 ! Time increment for SOLID analysis
FLDA, TIME, STEP, 100
FLDATA,ALGR,SEGR,SIMPLEN
FLDATA,ALGR,HFLM,TEMP
SOLVE
SAVE
FINISH
/SOLU
FSRS,TIME,10000 ! Restart time for FSI analysis
FSTI,12500 ! End time for FSI analysis
FSDT,100 ! Time increment for FSI analysis
SOLVE
SAVE
FINISH