第1期 李凤芹等:螺旋槽于气密封特性有限元分析 67 J, ,4 E一 、、一c= 图1螺旋槽于气密封的数学模型 图2螺旋槽干气密封部件 2.3基本方程 流体动力学基本方程纳维埃~斯托克斯方程 (Navie-Stokes Equation,简称N-S方程)的基本形式为: x一吉 +等未( 擎+ + a Du Du Du Du ...z y一古警+ ( (争+ + 争 =警+“ +“ Oy+“。警 z一丢 + ( + + a “ 、Du a“ Du Du 十“z 十“ 十 根据假设理论和所建的数学模型,又有气体的 连续性方程 + +一D puz=0和理想气体的 oZ o oZ 状态方程P = T,得到简化柱坐标方程为: 刍( )+r ( )=6 r2亘 (1) 已知第一类边界条件C1:当r=r 时,P=P ; 当r=ro时,P=P0。第二类周期性边界条件C2:P (r, )= (r,0+2 ̄r/N,) 2.4控制方程的离散和求解 为了使有限元分析具有一般性和避免计算中出 现天文数字的运算,减小数值计算误差,首先对方程 引入无量纲变量:P= / ;R=r/r H=h/ho代 人式(1),得无量纲化N—S方程: Ra 、( R )+ ’RaR、( “~aR )= P h a(P.O0 H) (2) 在利用加权余量法,则可将求解方程(2)的问题 变为一个求解积分方程的问题。 设未知压力场P的函数A(P): A(P)= ( )+ (RPH )一 =0 (3) 由于密封面上槽形不规则,很难得到上述方程 的解析解,为了得到比较精确的解析解,我们采用了 目前比较流行的有限元分析法,在密封环面区域上 对方程求解,为了获得有限元解,首先要获得有限元 方程,下面简要的介绍有限元方程的推导过程。将 计算区域分为N个单元,那么每个单元的压力P 可用单元节点的压力值P 和形函数N 来近似表 示,即 P= N尸 (i=1,2,…, ) (4) 把式(4)代人式(3)后,得到余量E,又已知压 力边界条件,采用前面的加权余量法推导有限元方 程,取权函数叫,,得到近似积分形式为: 1 wjEda=0( =1,2,…,”) (5) 最后得到有限元方程的单元有限元方程式: A +B =0( ,j.,k=1,2,…, )(6) 其中: =I『『H ( + 警)RdRdO B 一I『 删 N为计算区域节点数,△e为单元区域 单元有限元方程是有限元计算的基础,再按有 限元中的总体合成方法,可建立总体有限元方程: A +B! :0 7 根据文献[1]给出的实验用螺旋槽几何尺寸及 运行条件为fo=77.78mm;r =58.42mm;rg= 69mm;rb 61.3mm;口:15。;N =12;r=1;h = 5gm;h =3.05gm;N=10386r/min;Po= 4.5852MPa;P =0.1013MPa。又有上面的有限元 维普资讯 http://www.cqvip.com
西华大学学报・自然科学版 2008记 分析流程。便可得到端面间的某一点的气膜压力和 沿密封面的气膜压力分布,如图3所示。从图可见, 密封端面间的压力从外径岛内径逐渐上升,大约在 求出密封压力分布后,在此基础上便可求得密 封的承载能力、气膜刚度、功耗以及密封的泄漏量等 Bd芝 趟 性能参数,因而既可进行槽形参数的优化,得到良好 7 6 5 4 3 2 l 0 螺旋槽底部(即槽与密封坝的交界处)膜压力最高, 的气膜特性。 这与前面所示的密封机理描述的一致,图4即为所 3 结束语 计算的密封动环端面位置处沿半径方向的压力值与 文献的实验值的比较,可以看出二者比较接近。 从理论上讲,干气密封能够做到零泄露和零磨 损l2],但在实际应用中,也存在一定的磨损和泄露, 干气密封设计的根本点就在于如何选择最佳几何参 数,使密封达到或接近理想状态。 本文建立了合理的数学模型,推导了螺旋槽干 气密封气体稳态流动的有限元离散方程,并求解,最 后得到密封端面间的气膜压力分布图,并和文献中 的实验值相比较,从图中可以看到,螺旋槽干气密封 能产生较好的流体动压效应;该文的研究为以后槽 形参数的优化奠定了基础。 图3气膜压力分布云图 参考文献 [1]胡丹梅.低速干气端面密封参数优化及性能研究[D].山东: 石油大学,1996. [2]Gabriel R P.Fundamentals of sr)irla Groove non-contactin4g Face s ̄l[J].Lubrication Engineerign,1994,50(2):152—224. [3]王汝美.实用机械密封技术问答[M].2版.北京:中国石化 出版社,2004. [4]顾永泉.流体动密封[M].东营:石油大学出版社,1990. 58 60 62 64 66 68 70 72 74 76 78 [5]李人宪.有限元法基础[M].北京:国防工业出版社,2002. 半径/mm [6]王和顺.干气密封运行状态稳定性的研究[D].成都:西南交 图4气膜压力沿半径方向的分布图 通大学。2006. (编校:夏书林) ●● m●●● ●¨●¨◆●|¨●_I● (上接56页) [7]国家环境保护总局,水和废水监测分析方法编委会.水和废 参考文献 水监测分析方法[M].4版,北京:中国环境科学出版社,2002. [1]彭近新,陈慧君.水质富营养化与防治[M].北京:中国环境 [8]吴薇,刘雷,黄魁,等.国内外利用藻类去除污水中氮 科学出版社,1988. 磷的研究现状及发展趋势[J].江西科学,2007,25(1):47—50. [2]何少林,黄翔峰,李旭东,等.高效藻类塘藻菌共生系统的培 [9]李裕奇.概率论与数理统计(下)[M].北京:国防工业出版 养[J],中国给水排水,2006,22(5):74—76. 社。2004. [3]蔡艳荣,李孝民.稳定塘多池结构对藻类生长的影响[』].渤 f 10]Amah,Ginachi ljooma.Constructed Treatment Wedands for 海大学学报:自然科学版,2004,25(1):19 21. Reducign Nutrient Loadign in The Lower Malibu Creek Watershed. [4]Craggs R.J,Wastewater Treatment by AIgal Turf Scrubbing Dissertation Abstracts International【J].University of Califomia,Los [J].Water Science and Technology,.2001,44(11—12):427—433. Angeles,2004,65—08(B):4213. [5]MeDougal,Rhonda Louise.Algal Primary Production in [11]凌晓欢,况琪军,邱昌恩,等.两种藻类对水体氮、磷去除效 Prairie Wetlands:The Effects of Nutrients.Irradinace.Temperature 果[J].武汉大学学报:理学版,2006,52(4):487—491. nad Aquatic.Maerophytes.Dissertation Abstracts International[J].Uni— [12 3刘霞,杜贵森.藻类植物与水体富营养化控制[J].首都师 versity of Manitoba(Canada),2002,64.06(B):2501. 范大学学报,2002,23(4):56 60. [6]胡家骏,周群荚.环境工程微生物学[M].北京:高等教育出 [13]叶麟.三峡水库香溪河库湾富营养化及春季水华研究 版社,1988. [D].中国科学院研究生院,2006. (编校:李蓝)
因篇幅问题不能全部显示,请点此查看更多更全内容