大型水轮发电机通风发热综合计算

  • 格式:pdf
  • 大小:190.33 KB
  • 文档页数:4

第20卷第11期2000年11月中 国 电 机 工 程 学 报Proceedings of the CSEEVol.20No.11Nov.2000文章编号:025828013(2000)1120006204大型水轮发电机通风发热综合计算温嘉斌1,孟大伟1,鲁长滨2(哈尔滨理工大学,黑龙江省哈尔滨市150040;2.哈尔滨大电机研究所,黑龙江省哈尔滨市150040)SY NTHETIC CALCU LATION FOR THE VENTI LATION AN D HEATINGOF LARGE WATER WHEE L GENERAT ORWEN Jia2bin1,M EN G Da2wei1,L U Chang2bin2(1.Harbin University of Science and Technology,Harbin150040,China;2.Harbin Research Institute of Large Electric Machinery,Harbin150040,China)ABSTRACT:This paper takes the thermal deformation research and cooling technology of a large water wheel generator as back2 ground,and makes a synthetical calculation of the three2dimen2 sional temperature field and the ventilation system of the rotor of an electric machine.In the past,when the heating problem of a large water wheel generator was concerned,the tempera2 ture field was individually solved regardless of the interactions a2 mong thermal source,ventilation structure,temperature and its distribution.In this paper,a synthetical algorithm related to the ventilation and heat of an electric machine is adopted.The research work in this paper is very important to raise material u2 tilization during the machine design and manufacture and relia2 bility during the operation of the machine.KE Y WOR DS:large water wheel generator;calculation of venti2 lation;temperature field摘要:以大型水轮发电机热变形及冷却技术研究为背景,进行了大型水轮发电机转子三维温度场及电机通风系统的综合计算研究。

以往在进行大型水轮发电机发热分析时,都是孤立地求解温度场,而很少考虑热源及通风结构与温度及其分布之间的相互影响。

文中采用一种综合算法,将电机通风发热联系起来。

该研究工作对于提高电机设计制造中材料利用率以及电机运行中的可靠性,具有重要意义。

关键词:大型水轮发电机;通风计算;温度场中图分类号:TM301.4 文献标识码:A基金项目:国家自然科学基金(59493705)和机械工业技术发展基金(94J Z005)联合资助项目。

1 引言电机内的各种物理场是相互影响、相互制约的,是有一定耦合关系的综合场。

就电机部件内温度及其分布而言,除了其自身的材料性质之外,主要依赖于热源和电机通风情况。

一般孤立地求解电机内温度场而不考虑它们之间的相互影响,这将由于初始条件的变化而使计算结果失去准确性。

因而要想准确地计算出电机的温度分布,就必须将这些因素的相互影响在计算时加以综合考虑。

2 综合计算的总体思路及整体流程图电机内部的通风与温度场是通过对流换热的方式联系在一起的。

文中转子三维温度场计算采用有限元计算方法、电机通风计算采用通风网络法。

温度对热源的影响反映在电损耗上,由于材料的电阻率是温度的函数,因而电损耗也是温度的函数。

这个影响很容易计及,只需对温度场几次迭代计算就可以实现。

温度对材料导热系数的影响同样处理。

具体的实施过程如图1所示。

3 大型水轮发电机通风计算电机通风冷却系统的通风计算,理论上可以通过求解N2S方程及流体连续性方程,求解系统内冷却空气的三元流动问题,但实际上由于系统内冷却空气的过流情况十分复杂,边界条件不易确定,所以一般都化为带有集中参数的等值风路求解。

本文以网络矩阵法对漫湾电站SF250248/12200水轮发电机做了通风计算。

该电机的通风系统采用转子供风、密闭自循环、双路径向通风方式。

将其通风冷却系统内各个流道的基本等值风路组合起来,建立通风计算网络示于图2。

该通风计算网络由76个节点、145个等效支路及70个回路构成。

图1 综合计算整体流程图Fig.1 I nteger flow chart of the synthesis calculation 应用回路流量法,可列出方程组:AX =F(1)式中 X =[x 1,x 2,…,x m ]T 是回路流量矢量;A 是函数风阻矩阵;F 为风压升矢量,等于回路i 中各个风扇风压升的代数和。

表1列出部分支路风速的计算结果,并给出了实测值(来源于漫湾电站SF250248/12200型水轮发电机现场测试报告)及误差。

4 大型水轮发电机转子温度场计算对于大型水轮发电机转子温度场的计算,目前国内所采用的剖分单元都为正三棱柱单元,剖分时无法准确描述励磁绕组端部的圆弧部分,会对计算结果带来误差;同时,所选定的计算区域都未包括转子磁轭,而是将磁轭与磁极的交界面作为绝热面处理,忽略了两者间的热传递。

本文转子温度场求解区域包括了转子磁轭和励磁绕组端部的圆弧部分,并用六面体等参元及其退化的五面体对求解区域剖分,整体剖分图见图3。

4.1 温度场计算区域表面散热系数计算大型水轮发电机通风发热综合计算是通过温度场计算区域表面散热系数作为媒介连接起来的,因而正确地描述出表面散热系数与各部分风速的关系是十分重要的。

本文借用了文[1][2]给出的部分表面散热系数的计算方法;对于没有现成计算方法图2 通风计算网络Fig.2 C alculation netw ork of ventilation7 第11期 温嘉斌等: 大型水轮发电机通风发热综合计算表1 定子通风沟风速计算值与实测值对比T ab.1 Contrast of calculation and measure wind speed values in ventilation channels of stator支路号计算值/(m/s )实测值(平均)/(m/s )相对误差/%风沟编号风沟风速实测值/(m/s )13219.80417.95010.331216.319.613418.13917.280 4.973456714.619.617.715.818.713620.15918.7337.61891011121316.318.219.018.220.919.813819.14320.200 5.2315171919.819.321.514018.86118.270 3.2321232520.021.213.614117.86615.23317.2829313314.714.816.213918.14818.733 3.1235373916.719.220.313719.16619.375 1.074143454618.320.719.219.313519.59818.920 3.58474849505116.719.821.320.016.813317.15116.950 1.18525318.215.7图3 三维温度场整体剖分图Fig.3 G enerating meshes of the whole temperature f ield的表面散热系数,本文将依据流体力学和传热学理论推导出相应的计算公式。

(1)转子轭部表面和磁轭通风道内表面散热系数以往在大型水轮发电机转子温度场计算中,所选定的计算区域都未包括转子磁轭,因此转子轭部区域表面散热系数没有现成计算方法。

本文根据外层壁换热理论推导出了相应的计算公式[4]。

转子轭部表面散热系数的计算式为αl =ρcu ∞S t (2)式中 斯坦顿数的表达式为S t =c f2/[c f2(13.2P r -10.16)+0.9](3)其中c f /2可由如下经验式求出c f2=0.0287R -0.2ex(R ex ≤7×106)0.185(lg R ex )-2.584 (R ex >7×106)转子磁轭通风道内表面散热系数的计算式为α2=λD hN u (4)式中 D h 为当量直径;N u 为努谢尔特数,用下式计算N u =0.152R 0.9e P r0.833[2.25ln (0.114R 0.9e )+13.2P r -5.8](5)(2)转子极靴表面散热系数[2]α3=28(1+u 0.5δ)(6)式中 u δ为气隙平均风速。

(3)励磁绕组表面散热系数文[1]通过整体通风散热模型的测试得到大量实测数据,对这些实测数据做数学拟合处理后,得到不同入风角、不同位置的励磁绕组表面散热系数计算公式。

本文通过对这些计算公式进行插值求得表面单元散热系数。

4.2 材料特性及热传导方程在求解区域内共有8种媒质:转子磁极、阻尼条、励磁绕组、磁极压板、上下托板、极身绝缘、阻尼环和转子磁轭。

其中:由于匝间绝缘的影响归算到了其径向导热系数中,励磁绕组径向(y 轴)导热系数比周向(x 轴)和轴向(z 轴)的导热系数要小得8 中 国 电 机 工 程 学 报 第20卷多,为热的各向异性媒质;转子磁极和转子磁轭的轴向导热系数比叠片方向的导热系数亦要小,也为热的各向异性媒质;其他5种导热媒质中,由于各方向导热系数差别很小,均为热的各向同性媒质。

在这些媒质内,满足热传导方程:5x(k ix 5Tx)+5y(k iy5Ty)+55z(k iz 5T5z)=-qν(i=1,…,8)(7)式中 k1x、k1y和k1z表示励磁绕组在x、y、z方向上的导热系数;k2x、k2y、k2z和k3x、k3y、k3z分别为转子磁极和转子磁轭在x、y、z方向上的导热系数。