a54321 发表于 2004-3-16 17:08:17

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

请教各位大侠,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 发表于 2004-3-25 13:17:57

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

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

fossiler 发表于 2004-3-25 14:35:11

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

求解控制
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 发表于 2004-4-2 15:40:14

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

/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 发表于 2004-4-2 15:44:56

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

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

mack888 发表于 2004-8-20 15:10:57

mack888 发表于 2004-8-20 15:38:14

页: [1]
查看完整版本: 高压SF6断路器气流场求解的问题