高压SF6断路器气流场求解的问题

a54321

请教各位大侠，SF6断路器中的气流场是变粘性，变密度，跨音速，具有湍流，激波的流体，我现在想求解触头在不同开距下的速度，密度等的等值线图，请大家帮忙看看我设置流体特性的命令流哪里错了fldata12,prop,ufrq,10            !设置特性更新频率的标志
fldata12,prop,ivis,10
fldata13,vary,dens,t             !设置气体特性变化标志
fldata13,vary,visc,t
fldata13,vary,cond,t
fldata13,vary,spht,t
  
fldata7,prot,dens,gas            !密度
fldata8,nomi,dens,6.14
fldata9,cof1,dens,298
fldata10,cof2,dens,1.2367E5                          
!fldata11,cof3,dens,1
!fldata12,prop,ufro
fldata13,vary,dens,t
fldata14,temp,nomi,313
  
fldata7,prot,visc,gas            !粘性
fldata8,nomi,visc,0.02
fldata9,cof1,visc,298
fldata10,cof2,visc,1450
!fldata11,cof3,visc,0.5
!fldata12,prop,ufro
fldata13,vary,visc,t
fldata14,temp,nomi,313
  
fldata7,prot,cond,gas             !热传导率
fldata8,nomi,cond,0.013
fldata9,cof1,cond,298
fldata10,cof2,cond,1309
!fldata11,cof3,cond,0.5
!fldata12,prop,ufro
fldata13,vary,cond,t
fldata14,temp,nomi,313
  
fldata7,prot,spht,gas             !比热对于气体和液体为常数
fldata8,nomi,spht,96.6
fldata9,cof1,spht,293
fldata10,cof2,dens,298
!fldata11,cof3,dens,1
!fldata12,prop,ufro
fldata13,vary,spht,t
fldata14,temp,nomi,313
  
!以上确定了流体的状态如可压缩、变粘性、传热、变密度
  
!设置FLOTRAN分析参数（如迭代的次数、时间步长等）
!fldata1,solu,flow,t              !层流
fldata1,solu,turb,t               !湍流
fldata1,solu,temp,t               !求解温度方程
fldata1,solu,comp,t               !求解可压缩方程
fldata1,solu,swrl,t               !求解涡流方程
  
fldata2,iter,exec,100             !全局迭代次数
fldata2,iter,appe,100             !附加在结果文件之上的全局求解次数
  
!求解控制
flda,meth,pres,1                   !三对角矩阵求解压力
flda,meth,temp,1  
flda,meth,enke,1  
flda,meth,ends,1  
  
!提高收敛性和稳定性的工具，惯性松弛的设置
fldata26,stab,turb,1
fldata26,stab,pres,0.2            !压力松弛因子
fldata26,stab,mome,0.2            !动量松弛因子
  
!激活湍流模型
flda,turb,modl,1                  !k-e方程求解
fldata24,turb,inin,0.01
fldata24,turb,insf,0.01           !k-e二方程的求解系数
fldata24,turb,rati,10000          !分析开始时的湍流比
  
!限制速度峰值
fldata31,capp,temp,t
现在一求解就出现温度负值，郁闷啊！

fossiler

你的程序流不全呀
可以把前边的也贴出来
我看看,怎么办呀

fossiler

求解控制  
flda,meth,pres,1 !三对角矩阵求解压力  
flda,meth,temp,1  
flda,meth,enke,1  
flda,meth,ends,1  
  
是不是应该是fldata18呀,而且除了pres的求借器外其他都是默认的
  
!提高收敛性和稳定性的工具，惯性松弛的设置  
fldata26,stab,turb,1  
fldata26,stab,pres,0.2 !压力松弛因子  
fldata26,stab,mome,0.2 !动量松弛因子  
  
松弛因子好象比默认只大好几个数量机呀,
这个数值怎么确定呀,很大的
  
  t是不是在前边定义的参数呀
  
可以多多交流呀

a54321

/prep7   
rectng,,.675,x,1.068+x   
k,5,,.938+x  
k,6,.025 ,.938+x
k,7,,.908+x  
k,8,.0215,.908+x
k,9,,.819+x  
k,10,.0215 ,.819+x   
k,11,,.816+x
k,12,.0205,.816+x   
k,13,,.774+x
k,14,.0205,.774+x   
k,15,,.766+x
k,16,.0175 ,.766+x   
k,17,,.7+x   
k,18,.0175 ,.71452175972+x   
k,19,.007583333329  ,.701130479385+x
k,33,.025  ,1.068+x  
l,5,6   
l,6,33   
l,6,8   
l,4,5   
l,5,7   
l,7,8   
l,7,9   
l,9,10   
l,8,10   
l,9,11   
l,11,12  
l,10,12  
l,11,13  
l,13,14  
l,12,14  
l,13,15  
l,15,16  
l,14,16  
l,15,17  
l,16,18  
!k,19,.007583333329  ,.701130479385+x   
larc,18,19,15,.014   
larc,19,17,15,.026   
k,20,.09,1.068+x
k,21,.11,1.068+x
k,22,.09 ,.6955051026+x  
k,23,.11 ,.735+x
k,24,.125,.735+x
k,25,.14  ,.72+x
k,26,.098 ,.6857071436+x
k,27,.105,.685+x
l,20,21  
l,20,22  
l,21,23  
l,23,24  
larc,24,25,22,.015   
larc,22,26,23,.01   
larc,25,26,27   
k,28,.17 ,.685+x
k,29,.17 ,.710+x
!k,29,.145,.71+x
k,30,.27 ,.785+x
k,31,.27,1.068+x
!larc,28,29,30,.025  
circle,29,.025   
larc,28,30,31,.1
l,30,31  
k,40,.02 ,x  
k,41,.02 ,.402+x
k,42,.018,.408+x
k,43,.013,.415+x
k,44,.013,.425+x
k,45,.0175,.431+xx   
k,46,.0175,.4871239034+x
k,47,.02240747729 ,.4956239034+x
k,48,.0235,.4970816275+x
k,49,.029  ,.5+x
k,50,.043,.486+x
k,51,.043,.441+x
k,52,.058,.426+x
l,40,41  
l,41,42  
l,42,43  
l,43,44  
l,44,45  
l,45,46  
l,46,47  
larc,47,48,50,.005   
larc,48,50,49   
l,50,51  
larc,51,52,50,.015   
k,54,.1315,x
k,55,.1315,.387+x   
k,56,.1165 ,.402+x   
k,57,.1165,.412+x   
k,58,.1125,.416+x   
k,59,.094,.416+x
k,60,.09,.42+x   
l,54,55  
larc,55,56,54,.015   
l,56,57  
larc,57,58,56,.004   
l ,58,59
larc,59,60,58,.004   
!   
!   
k,61,.03826764281,.645+x
k,62,.046,.645+x
k,63,.0339 ,.633+x   
k,64,.025,.6176+x   
k,65,.034,.602+x
k,66,.025,.588+x
k,67,.032,.576+x
k,68,.023,.562+x
k,69,.02875,.552+x   
k,70,.021 ,.54+x
k,71,.021 ,.53+x
k,72,.03 ,.53+x  
k,73,.046,.537+x
l,73,62  
l,62,61  
l,61,63  
l,63,64  
l,64,65  
l,65,66  
l,66,67  
l,67,68  
l,68,69  
l,69,70  
l,70,71  
l,71,72  
k,90,.03,.522+x  
l,72,90  
lfillt,57,58,.005   
lfillt,58,59,.003   
lfillt,59,60,.003   
lfillt,61,60,.003   
lfillt,61,62,.003   
lfillt,63,62,.003   
lfillt,63,64,.003   
lfillt,65,64,.003   
lfillt,65,66,.003   
lfillt,67,68,.002   
lfillt,68,69,.003   
k,93,.073,.51+x  
k,94,.073 ,.4655+x   
k,95,.071,.462+x
k,96,.071,.437+x
k,97,.078,.437+x
k,98,.078,.415+x
k,99,.056,.5001834096+x  
k,100,.056,.436+x   
k,101,.071,.415+x   
l,73,93  
l,93,94  
l,95,94  
l,95,96  
l,96,97  
l,97,98  
l,90,99  
l,99,100
l,100,101   
l,98,101
lfillt,57,81,.01
lfillt,82,81,.015   
lfillt,84,85,.005   
lfillt,87,88,.015   
lfillt,89,88,.01
k,112,.09,.455+x
k,113,.082,.463+x   
k,114,.08,.463+x
k,115,.074,.4602915026+x
k,116,.073,.4576457513+x
l,60,112
larc,112,113,60,.008
l,113,114   
larc,114,115,112,.008   
larc,115,116,112,.004   
lang,89,52,125   
k,118,.17,.5+x   
k,119,.145,.475+x   
k,120,.27,.4+x   
k,121,.27,x  
k,122,.145,x
k,123,.145,1.068+x   
l,121,120   
larc,120,118,121,.1  
larc,118,119,122,.025   
l,119,122   
l,4,33   
l,20,21  
l,123,31
l,40,54  
l,122,121   
l,123,35
lglue,38,36,111  
lang,90,116,90   
linl,113,93  
l,116,124   
larc,106,124,116,.005   
larc,107,124,116,.005   
FLST,2,42,4  
FITEM,2,86   
FITEM,2,90   
FITEM,2,38   
FITEM,2,101  
FITEM,2,89   
FITEM,2,95   
FITEM,2,88   
FITEM,2,85   
FITEM,2,114  
FITEM,2,113  
FITEM,2,84   
FITEM,2,83   
FITEM,2,82   
FITEM,2,92   
FITEM,2,94   
FITEM,2,81   
FITEM,2,87   
FITEM,2,91   
FITEM,2,69   
FITEM,2,80   
FITEM,2,68   
FITEM,2,79   
FITEM,2,67   
FITEM,2,66   
FITEM,2,57   
FITEM,2,78   
FITEM,2,65   
FITEM,2,77   
FITEM,2,64   
FITEM,2,76   
FITEM,2,63   
FITEM,2,75   
FITEM,2,62   
FITEM,2,74   
FITEM,2,61   
FITEM,2,73   
FITEM,2,60   
FITEM,2,59   
FITEM,2,72   
FITEM,2,71   
FITEM,2,58   
FITEM,2,70   
AL,P51X  
FLST,2,31,4  
FITEM,2,109  
FITEM,2,40   
FITEM,2,41   
FITEM,2,42   
FITEM,2,43   
FITEM,2,44   
FITEM,2,45   
FITEM,2,46   
FITEM,2,47   
FITEM,2,48   
FITEM,2,49   
FITEM,2,50   
FITEM,2,102  
FITEM,2,101  
FITEM,2,90   
FITEM,2,38   
FITEM,2,86   
FITEM,2,85   
FITEM,2,114  
FITEM,2,93   
FITEM,2,99   
FITEM,2,100  
FITEM,2,98   
FITEM,2,97   
FITEM,2,96   
FITEM,2,56   
FITEM,2,55   
FITEM,2,54   
FITEM,2,53   
FITEM,2,52   
FITEM,2,51   
AL,P51X  
FLST,2,5,4   
FITEM,2,105  
FITEM,2,106  
FITEM,2,104  
FITEM,2,103  
FITEM,2,110  
AL,P51X  
FLST,2,5,4   
FITEM,2,39   
FITEM,2,112  
FITEM,2,36   
FITEM,2,111  
FITEM,2,108  
AL,P51X  
FLST,2,7,4   
FITEM,2,28   
FITEM,2,29   
FITEM,2,30   
FITEM,2,31   
FITEM,2,33   
FITEM,2,32   
FITEM,2,27   
AL,P51X  
FLST,2,17,4  
FITEM,2,107  
FITEM,2,8   
FITEM,2,6   
FITEM,2,7   
FITEM,2,9   
FITEM,2,13   
FITEM,2,11   
FITEM,2,16   
FITEM,2,14   
FITEM,2,17   
FITEM,2,19   
FITEM,2,22   
FITEM,2,20   
FITEM,2,23   
FITEM,2,24   
FITEM,2,25   
FITEM,2,26   
AL,P51X  
FLST,2,7,5,ORDE,2   
FITEM,2,1   
FITEM,2,-7   
AOVLAP,P51X  
k,1000,0.056,0.441   
l,111,1000   
l,51,1000   
al,50,102,89,95,1,3   
asba,8,1,sepo
l,111,1000
al,50,102,89,95,1,4
aadd,1,3,2  
numcmp,areaet,1,141,,,1
lesize,4,,,100
lesize,28,,,100
lesize,115,,,200
et,1,141,,,1
lesize,4,,,100
lesize,28,,,100
lesize,115,,,200
asel,s,,,6
aatt,1,,1
mshape,1,2D
mshkey,0
!smrtize,1
amesh,6
aplot
fldata12,prop,ufrq,10            !设置特性更新频率的标志
fldata12,prop,ivis,10
fldata13,vary,dens,t             !设置气体特性变化标志
fldata13,vary,visc,t
fldata13,vary,cond,t
fldata13,vary,spht,t
  
fldata7,prot,dens,gas            !密度
fldata8,nomi,dens,6.14
fldata9,cof1,dens,298
fldata10,cof2,dens,1.2367E5                          
!fldata11,cof3,dens,1
!fldata12,prop,ufro
fldata13,vary,dens,t
fldata14,temp,nomi,313
  
fldata7,prot,visc,gas            !粘性
fldata8,nomi,visc,0.02
fldata9,cof1,visc,298
fldata10,cof2,visc,1450
!fldata11,cof3,visc,0.5
!fldata12,prop,ufro
fldata13,vary,visc,t
fldata14,temp,nomi,313
  
fldata7,prot,cond,gas             !热传导率
fldata8,nomi,cond,0.013
fldata9,cof1,cond,298
fldata10,cof2,cond,1309
!fldata11,cof3,cond,0.5
!fldata12,prop,ufro
fldata13,vary,cond,t
fldata14,temp,nomi,313
  
fldata7,prot,spht,gas             !比热对于气体和液体为常数
fldata8,nomi,spht,96.6
fldata9,cof1,spht,293
fldata10,cof2,dens,298
!fldata11,cof3,dens,1
!fldata12,prop,ufro
fldata13,vary,spht,t
fldata14,temp,nomi,313
  
!以上确定了流体的状态如可压缩、变粘性、传热、变密度
fldata1,solu,turb,t               
fldata1,solu,temp,t               
fldata1,solu,comp,t               
fldata1,solu,swrl,t               
  
fldata2,iter,exec,100              
fldata2,iter,appe,100              
  
!求解控制
flda,meth,temp,1  
flda,meth,enke,1  
flda,meth,ends,1  
  
!激活湍流模型
flda,turb,modl,1                  
fldata24,turb,inin,0.01
fldata24,turb,insf,0.01            
fldata24,turb,rati,10000           
  
!限制速度峰值
fldata31,capp,temp,t
  
!施加边界条件
lsel,s,,,115
nsll,s,1
d,all,vx,0                          
lsel,s,,,4
nsll,s,1
d,all,vx,0
d,all,vy,2
!d,all,pres,3e5                     
lsel,s,,,116
lsel,a,,,119
nsll,s,1
d,all,pres,0                        
  
asel,s,,,6
lsla,s
lsel,u,,,4
lsel,u,,,115
lsel,u,,,116
lsel,u,,,119
nsll,s,1
d,all,vx,0
d,all,vy,0                                                   
d,all,pres,6e5 fini
/solu
alls
solve

a54321

希望大家给点建议，我的程序流文件现在结果不正确，我不知道错在那里，温升过高，发散！不过我认为模型的建立和剖分很不错哦！我希望我们能共同进步！同时谢谢斑竹的加分鼓励。