CN101673318B - 透平膨胀机径向静压气体轴承最佳静态参数设计方法 - Google Patents
透平膨胀机径向静压气体轴承最佳静态参数设计方法 Download PDFInfo
- Publication number
- CN101673318B CN101673318B CN2009101851363A CN200910185136A CN101673318B CN 101673318 B CN101673318 B CN 101673318B CN 2009101851363 A CN2009101851363 A CN 2009101851363A CN 200910185136 A CN200910185136 A CN 200910185136A CN 101673318 B CN101673318 B CN 101673318B
- Authority
- CN
- China
- Prior art keywords
- bearing
- coefficient
- equation
- static parameter
- gas
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Magnetic Bearings And Hydrostatic Bearings (AREA)
Abstract
Description
技术领域
本发明涉及轴承设计领域,可以用于确定透平膨胀机径向节流小孔静压气体轴承的最佳静态参数,适用于各种工质气体透平膨胀机径向节流小孔静压气体轴承的静态设计,亦可以用于径向节流小孔静压气体轴承的性能计算。
背景技术
气体轴承具有转速高、摩擦小和污染少的特点,是适用于透平膨胀机中的最理想形式,也是其核心部件。静压气体轴承相比动压气体轴承具有承载能力高、运行稳定的特点。受到国内加工水平的限制,国内目前还没有动压气体轴承应用的事例,节流小孔静压气体轴承是当前国内用于透平膨胀机中最常用的形式。但对于高转速、有特殊用途的透平膨胀机的设计国内还没有完全掌握,气体轴承设计方法的核心技术被国外少数几个跨国公司所垄断,例如国内的低温氦透平膨胀机主要还得依靠进口,其主要原因就是国内有关透平膨胀机静压气体轴承的设计方法都很零散,特别是低温氦透平膨胀机气体轴承的设计还是基于空气轴承的设计理念和方法,并且本身的设计精度就不高,往往只考虑几个性能参数就确定最终设计方案,设计出来的参数跟实际情况差别很大,例如在低温氦透平膨胀机设计中,忽略的气膜摩擦功率损耗产生的热量对工作轮端的效率影响很大。
发明内容
本发明的目的是提供一种透平膨胀机径向静压气体轴承最佳静态参数设计方法,针对目前透平膨胀机径向节流小孔静压气体轴承的设计方法精度不高,设计方法在推导过程中忽略的因素很多。用一个程序计算得到四个反映径向节流小孔静压气体轴承的性能参数,在保证承载能力和刚度系数最大,耗气量和摩擦损耗最小的前提下,最终确定径向节流小孔静压气体轴承的最佳设计方案,从而获得很高的计算精度。
为了达到上述目的,本发明所采用的技术方案为:
透平膨胀机径向静压气体轴承最佳静态参数设计方法,其特征在于:将轴颈轴承分解成许多“等效狭缝”,假设气膜厚度方向为y,气体无环向流动,只是沿x方向即轴向流向端面,节流小孔的扩散效应和环向流动影响的修正系数fw,节流小孔节流后的压力为pdi,端面的环境压力为pa。其它静态参数如下:供气压力ps,轴承长度L,轴颈轴径D,节流小孔孔径d,工质运动粘性系数μ,偏心率e,平均气膜厚度h0,每排小孔个数n,小孔到端面距离l1、l2,轴承承载能力W,轴承刚度系数KW,气体消耗量摩擦功率损耗FT,转子的转速ω,下标i为每个节流小孔的参数,下标1、2表示各排节流小孔的参数。
列出气体运动的微分方程,从微分方程入手,得到承载能力系数CW、轴承刚度系数KW、总流量和总摩擦功率损耗FT的方程式,并通过承载能力系数CW、轴承刚度系数KW、总流量和总摩擦功率损耗FT的方程式得到在不同偏心率条件下透平膨胀机径向节流小孔静压气体轴承的最佳静态参数,具体步骤如下:
(1)选定节流小孔类型,所述节流小孔形式有简单小孔和环面小孔两种,在计算节流小孔面积时两种节流小孔有所区别,简单小孔节流面积计算公式为: 环面小孔节流面积计算公式为:A=πdh,另外确定供气工质,根据供气工质可确定运动气体的运动粘性系数μ、密度ρ、绝热指数k;
(2)假设一组给定的静态参数:偏心率ε,小孔直径d,供气压力ps,供气温度Ts,排气压力pa,排气温度Ta,轴承长度L,轴承直径D,小孔到断面的距离l1、l2,以及一组间隙hj值,通过每个hj可以计算得到一个对应的压降比βi,进而确定每个间隙hj值对应的轴承刚度系数KW·i值,并做出KW·i~hj曲线图,求得拐点,此值即为初步确定的最佳间隙hop;
(3)保持间隙hop和步骤(2)所述静态参数中的某些静态参数不变,改变其中一个静态参数,代入方程
可以确定承载能力系数CW,将所求压降比βi代入方程
可以确定轴承刚度系数KW,将所求压降比βi代入方程
(4)重复步骤(3),通过改变步骤(2)所述的其中一个静态参数,做出各个变化的静态参数与承载能力系数CW、轴承刚度系数KW、总流量和总摩擦功率损耗FT之间的关系曲线,综合考虑CW、KW、FT,遵循CW和KW要求尽可能大,和FT要求可能小的原则,确定最终的最佳静态参数;
(5)根据步骤(4)所得到的最佳静态参数,代入修正系数方程
根据选取的设计余量选择方程fw1或fw2,可计算得到修正系数值,再根据步骤(3)所得到的承载能力系数CW、轴承刚度系数KW、总流量和总摩擦功率损耗FT的值,承载能力系数乘以修正系数fw,求得最终的轴承的CW、KW、FT。
本发明假设气膜厚度方向为y,气体无环向流动,只是沿x方向即轴向流向端面,从轴承间隙内气体的雷诺方程入手,在将轴颈轴承分解成许多“等效狭缝”的假设前提下,运动方程可以表示如下:
(1)式两次对y积分,并且利用边界条件y=0, y=hi,u=0;
(2)
质量连续方程:
气膜内满足等温条件:
(2)式代入(3)式,由边界条件x=l1,x=l2时,p=pa,利用(4)式消去ρ再对x从0到x积分:
(5)
(6)
如果忽略转速影响,(5)、(6)式变为:
(7)
(8)
设在任意的等分点中,x为固定值的全宽度上压力不变,单位轴向长度的合力作用在过该等分点半径上,则第i等分的总合力为:
(9)
其中:
其中:pdi1、pdi2分别表示各排节流小孔后的压力
假设节流小孔均匀分布,且有一对小孔分布在垂直方向,设(9)式的合力与基准线的夹角为θi,把(7)式和(8)式代入并积分,得总合力为:
(10)
其中:
如果是单排供气孔,L-l1-l2=0,下面各公式中亦相同。
用CW表示无因次承载能力系数:
(12)
各节流小孔处的气膜厚度为:
hi=hop(1-εcosθi) (13)
其中:hop为确定的轴承间隙;ε为偏心率;
每等分的轴承刚度系数为:
(14)
总的刚度系数为:
(15)
小孔节流过程是一个近似气体绝热过程,从伯努力方程入手,求得节流小孔中的流速为:
(16)
从单个小孔出中流出的气体质量流量为:
(17)
总的流量为:
(18)
其中N为总的小孔数目
小孔节流的堵塞条件为:
(19)
其中βk为滞止压比,综合(17)式和(19)式,经过小孔的质量流量为:
(20)
对于简单小孔形式: 环面小孔形式:A=πdhi,φ为速度系数,一般取0.85,并且在设计过程中会尽量避免出现小孔堵塞。
(21)
(21)式中只有βi是未知数,其他参数均为已知量,l分别取l1和l2,通过一个解非线性方程的小程序即可求得βi1和βi2。
以上径向节流小孔静压气体轴承的承载能力、耗气量和轴承刚度系数已经求解完毕,承载能力乘以修正系数fw,选fw1或是fw2由选取的设计余量而定,就可以得到最终轴承的承载能力。
(22)
下面将计算径向节流小孔静压气体轴承的摩擦功率损耗,由于气体黏度系数相比油的黏度系数小3个数量级,x方向上压力变化引起的摩擦功率损耗很小,可忽略,则(2)式两边对y求导得:
由牛顿应力公式,得:
摩擦功率力矩为:
摩擦功率损耗为:
本发明与现有技术相比的优点在于:本发明能实现一个程序计算得到所有反映径向节流小孔静压气体轴承的四个性能参数,这样确定的最佳静态参数在工程应用中精度很高,可以适用于高转速、特殊用途的透平膨胀机径向节流小孔静压气体轴承设计;并且本发明还可以用于在确定静态参数的前提下,计算不同偏心率下的径向节流小孔静压气体轴承承载能力、轴承刚度系数、气体消耗量和摩擦功率损耗。
附图说明
图1为一种透平膨胀机节流小孔静压气体轴承的结构图,其中红色箭头为气流流动方向。
图2为利用本发明计算的一台氦透平膨胀机径向节流小孔静压气体轴承的结构图。
图3为利用本发明改进后的此台氦透平膨胀机径向节流小孔静压气体轴承的结构图。
图4为本发明方法的设计流程图。
具体实施方式
如图2为现有的一台进口俄罗斯油气混合氦透平膨胀机径向节流小孔静压气体轴承,具体参数如下:单排供气,每排供气小孔n=8,小孔直径d=0.5mm,供气压力ps=19bar,直径气膜厚度h0=0.045mm,轴承长度L=28.5mm,轴颈直径D=12mm,额定转速180krpm。根据以上参数,按照本发明设计方法的设计流程图4,气体轴承的性能参数计算结果如下:总流量 刚度系数KW=2.33,承载能力W=1.06×102N,摩擦功率损耗FT=8.13W,以上结果跟俄罗斯透平公司提供的性能参数很接近,计算得到的流量跟实际测量值也很吻合,但从KW·i~hj图判断,直径气膜厚度h0偏离最佳间隙hop,其他参数也偏离最佳设计参数,是导致目前此氦气透平膨胀机运行不稳定的主要原因。
图3为利用本发明改进后的俄罗斯氦透平膨胀机径向节流小孔静压气体轴承,按照本发明设计流程图4改进的静压气体轴承参数如下:直径气膜厚度h0=0.04mm,双排供气,每排供气小孔n=8,小孔直径d=0.45mm,小孔到断面的距离l1=9.5mm,l2=15mm,改进后的性能计算结果如下:总流量 增加34.8%;刚度系数KW=2.92,提高25.3%;承载能力W=1.29×102N,提高21.7%;摩擦功率损耗FT=9.15W,增加11.15%。
Claims (1)
1.透平膨胀机径向静压气体轴承最佳静态参数设计方法,其特征在于:将轴颈轴承分解成许多“等效狭缝”,假设气膜厚度方向为y,气体无环向流动,只是沿x方向即轴向流向端面,节流小孔的扩散效应和环向流动影响的修正系数fw,节流小孔节流后的压力为pdi,端面的环境压力为pa。其它静态参数如下:供气压力ps,轴承长度L,轴颈轴径D,节流小孔孔径d,工质运动粘性系数μ,偏心率e,平均气膜厚度h0,每排小孔个数n,小孔到端面距离l1、l2,轴承承载能力W,轴承刚度系数KW,气体消耗量 摩擦功率损耗FT,转子的转速ω,下标i为每个节流小孔的参数,下标1、2表示各排节流小孔的参数;
列出气体运动的微分方程,从微分方程入手,得到承载能力系数CW、轴承刚度系数KW、总流量 和总摩擦功率损耗FT的方程式,并通过承载能力系数CW、轴承刚度系数KW、总流量 和总摩擦功率损耗FT的方程式得到在不同偏心率条件下透平膨胀机径向节流小孔静压气体轴承的最佳静态参数,具体步骤如下:
(1)选定节流小孔类型,所述节流小孔形式有简单小孔和环面小孔两种,在计算节流小孔面积时两种节流小孔有所区别,简单小孔节流面积计算公式为: 环面小孔节流面积计算公式为:A=πdh,另外确定供气工质,根据供气工质可确定运动气体的运动粘性系数μ、密度ρ、绝热指数k;
(2)假设一组给定的静态参数:偏心率ε,小孔直径d,供气压力ps,供气温度Ts,排气压力pa,排气温度Ta,轴承长度L,轴承直径D,小孔到断面的距离l1、l2,以及一组间隙hj值,通过每个hj可以计算得到一个对应的压降比βi,进而确定每个间隙hj值对应的轴承刚度系数KW·i值,并做出KW·i~hj曲线图,求得拐点,此值即为初步确定的最佳间隙hop;
(3)保持间隙hop和步骤(2)所述静态参数中的某些静态参数不变,改变其中一个静态参数,代入方程
可以确定承载能力系数CW,将所求压降比βi代入方程
可以确定轴承刚度系数KW,将所求压降比βi代入方程
(4)重复步骤(3),通过改变步骤(2)所述的其中一个静态参数,做出各个变化的静态参数与承载能力系数CW、轴承刚度系数KW、总流量 和总摩擦功率损耗FT之间的关系曲线,综合考虑CW、KW、 FT,遵循CW和KW要求尽可能大, 和FT要求可能小的原则,确定最终的最佳静态参数;
(5)根据步骤(4)所得到的最佳静态参数,代入修正系数方程
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009101851363A CN101673318B (zh) | 2009-09-30 | 2009-09-30 | 透平膨胀机径向静压气体轴承最佳静态参数设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009101851363A CN101673318B (zh) | 2009-09-30 | 2009-09-30 | 透平膨胀机径向静压气体轴承最佳静态参数设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101673318A CN101673318A (zh) | 2010-03-17 |
CN101673318B true CN101673318B (zh) | 2012-01-04 |
Family
ID=42020540
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009101851363A Expired - Fee Related CN101673318B (zh) | 2009-09-30 | 2009-09-30 | 透平膨胀机径向静压气体轴承最佳静态参数设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101673318B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102779112A (zh) * | 2012-05-31 | 2012-11-14 | 兰州石化职业技术学院 | 干气密封拟合曲线槽型线的生成方法 |
CN104331553B (zh) * | 2014-10-29 | 2017-07-11 | 浙江大学 | 考虑缺陷的大型透平膨胀机叶轮叶片结构优化设计方法 |
TWI641768B (zh) * | 2017-06-09 | 2018-11-21 | 國立清華大學 | 薄膜節流器及液靜壓軸承模組 |
CN108804842B (zh) * | 2018-06-15 | 2022-06-03 | 上海卫星工程研究所 | 一种基于系统工程的气体静压轴承工程设计方法 |
CN110848259B (zh) * | 2019-12-17 | 2024-07-02 | 中国工程物理研究院机械制造工艺研究所 | 一种节流效应可调控的静压气体止推轴承 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101280801A (zh) * | 2008-04-25 | 2008-10-08 | 西安交通大学 | 一种具有径向支撑箔片的动压气体轴承结构 |
CN101532536A (zh) * | 2009-04-22 | 2009-09-16 | 上海启元空分技术发展有限公司 | 一种静压定向辐射槽止推气体轴承 |
-
2009
- 2009-09-30 CN CN2009101851363A patent/CN101673318B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101280801A (zh) * | 2008-04-25 | 2008-10-08 | 西安交通大学 | 一种具有径向支撑箔片的动压气体轴承结构 |
CN101532536A (zh) * | 2009-04-22 | 2009-09-16 | 上海启元空分技术发展有限公司 | 一种静压定向辐射槽止推气体轴承 |
Non-Patent Citations (1)
Title |
---|
刘立强,熊联友,侯予 等.透平膨胀机新型切向小孔气体轴承稳定性的实验研究.《低温工程》.中国航天科技集团公司第101研究所,1999,(第110期),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN101673318A (zh) | 2010-03-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Xie et al. | Theoretical and experimental exploration into the fluid structure coupling dynamic behaviors towards water-lubricated bearing with axial asymmetric grooves | |
CN101673318B (zh) | 透平膨胀机径向静压气体轴承最佳静态参数设计方法 | |
Picardo et al. | Rotordynamic coefficients for a tooth-on-stator labyrinth seal at 70 bar supply pressures: measurements versus theory and comparisons to a hole-pattern stator seal | |
Amiri et al. | Experimental investigation of the hydraulic loads on the runner of a Kaplan turbine model and the corresponding prototype | |
Zhang et al. | Experimental study of the leakage and rotordynamic coefficients of a long-smooth seal with two-phase, mainly oil mixtures | |
CN110096784B (zh) | 一种具有轴向压差的径向滑动轴承的快速计算与设计方法 | |
Smith | Journal bearings in turbomachinery | |
San Andrés et al. | On the Leakage and Rotordynamic Force Coefficients of Pump Annular Seals Operating With Air/Oil Mixtures: Measurements and Predictions | |
Zuo et al. | Comparative performance analysis of conical hydrostatic bearings compensated by variable slot and fixed slot | |
Spakovszky et al. | Scaling laws for ultra-short hydrostatic gas journal bearings | |
Zhang et al. | Application of a novel rotordynamic identification method for annular seals with arbitrary elliptical orbits and eccentricities | |
CN101498214B (zh) | 一种偏心圆管环空幂律流体流态的测算方法 | |
CN110378018B (zh) | 一种液体动静压球轴承的稳态承载能力的计算方法 | |
Kaneko et al. | Experimental study on static and dynamic characteristics of liquid annular convergent-tapered damper seals with honeycomb roughness pattern | |
CN103823788B (zh) | 一种线性变形下静压转台承载能力分析方法 | |
Sari et al. | Performance of undershot waterwheel in pico scale with difference in the blades number | |
CN108520123B (zh) | 基于全流场计算的大功率离心泵转子动力特性分析方法 | |
Duffin et al. | Paper 4: Some Experimental and Theoretical Studies of Journal Bearings for Large Turbine-Generator Sets | |
CN104933219A (zh) | 一种水润滑推力轴承的全流态模拟参数获取方法 | |
Ligrani et al. | Miniature viscous disk pump: Performance variations from non-Newtonian elastic turbulence | |
Wu et al. | Hydrodynamic characteristics of cylindrical spiral grooves for pump annular seals | |
CN110427666B (zh) | 一种计算液体动静压球轴承最大油膜承载能力的方法 | |
Martin | Extended expansion of the Reynolds equation | |
Li et al. | On the lubrication performance of journal bearing coupled with the axial movement of journal | |
Benti et al. | Cost-Effective Design Modification of a Sleeve Bearing with Large Bearing Clearance |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120104 Termination date: 20130930 |