汽轮机可倾瓦轴承油膜特性研究综述

本文格式为Word版，下载可任意编辑——汽轮机可倾瓦轴承油膜特性研究综述

本文介绍了国内外研究者对于滑动轴承所做的研究工作，以及现在对于汽轮机可倾瓦滑动轴承数值分析的一些空白，给出了对于用数值模拟研究汽轮机可倾瓦滑动轴承油膜的建议。

汽轮机；可倾瓦；滑动轴承；数值模拟

Thispaperintroducesthestudiesconductedbyresearchersfromthedomesticandforeignfortheslidingbearingandthegapsofnumericalanalysisforthesteamturbinetiltingslidingbearing，whichgivesadvicesofresearchingsteamturbinetiltingtileslidingbearingoilfilmusingthenumericalsimulation.

Steamturbine；Tiltingtile；Slidingbearing；Thenumericalsimulation

0前言

滑动轴承有轴颈和轴瓦组成，相比滚动轴承来说，其生产本金低、工艺简朴、可以反复使用。液体在楔形空间内滚动时会产生很强的压力，滑动轴承就是利用这个原理工作的。轴颈的旋转带动油液一起旋转，随着旋转方向，油液流入由轴颈和轴瓦组成的楔形空间内，进而使轴颈在一个相对稳定的位置旋转，同时轴颈与轴瓦之间形成了一层薄薄的油膜[1]。滑动轴承的布局有大量种，但是其外形大致体验了圆形轴瓦、椭圆轴瓦、可倾瓦轴瓦，不同轴瓦油的膜特性不同。汽轮机运行时经常发生油膜失稳，引发了人们对轴瓦的研究，油膜在圆形轴瓦中运行时会形成一个高压区和一个空穴区，且轴系振动较大。而椭圆轴瓦由于其特有的椭圆度，轴系振动有所改善。可倾瓦滑动轴承的研发对于轴系润滑来说具有划时代的意义，其瓦块随着轴颈位置不断调整，不仅保证了油膜的稳定性，同时也使轴承具有更大的承载力[2]。

对于可倾瓦滑动轴承的研究，大多数研究者使用Reynolds方程计算单个瓦块，然后再举行叠加，这样计算的结果与实际油膜特性有较大偏差。由于他们疏忽好多因素，例如轴瓦间隙的影响。计算流体力学被近期的大量研究者使用，但是对于整个可倾瓦滑动轴承油楔的研究还很少见，还没有合理的得出可倾瓦滑动轴承的油膜特性以及在外界因素影响下油膜特性的变化。基于数值模拟方法可以透露这些特性。本文总结了可倾瓦滑动轴承的研究概况与研究方法，并对汽轮机可倾瓦滑动轴承的数值分析给出了建议。

1可倾瓦滑动轴承的研究概况

1.1国外研究概况

对可倾瓦滑动轴承的研究，国外开头的对比早，更加近些年进展很快。Yang[3]通过试验研究了汽轮机可倾瓦滑动轴承轴瓦疲乏破坏对轴瓦振动的影响，提出了槽式轴瓦对轴瓦举行了提升，并验证了提升后的可倾瓦轴承轴瓦性能显著提高。但是仅仅分析了轴瓦性能，没有对油楔中油膜举行分析；MatthewCha[4]等利用数值模拟方法对雷诺方法举行有限元求解，分析了支点对可倾瓦滑动轴承非线性油膜力的影响，得出了油膜厚度与油膜压力的关系。但是对于可倾瓦轴承其他因素的影响，没有举行分析。Kyung-BoBang[5]等通过测验研究了传统六瓦可倾瓦轴承油膜温度与功率损耗之间的关系，通过变更可倾瓦轴承轴颈转速、轴承载荷和润滑油进口速度的大小，得出了这三者分别对轴承功率损耗的影响大小。但是，没能给出其他可倾瓦功率损耗与轴瓦温度的关系。White和Chan[6]研究了轴颈振动频率与动力特性的关系，得出动力特性的频率与支点偏置度和负荷成反比，但是计算仅限于Sommerfeld数范围内，导致结论使用的局限性；S.Wen等[7]用含有Newton-Raphson方法的数值计算方法对静态可倾瓦滑动轴承TEHD性能举行研究，察觉这种方法比单个瓦块叠加收敛速度快。但是，没能很好的考虑瓦块与瓦块之间的间隙，还有润滑油温升等问题；A.M.EI-Butch等[8]用TEHD模型研究了轴线偏斜对油膜厚度的影响，察觉轴线偏斜对油膜产生的影响可以用瓦块的塑性变形举行中和。但是可倾瓦滑动轴承瓦块有点支撑和线支撑，而对于点支撑的瓦块未举行分析，并且求解时疏忽较多其他因素的影响；D.c.zachariadisl[9]通过求解Reynolds方程，研究轴向楔形与角偏斜的关系，得出油膜刚度与阻尼系数在很大程度上受角偏斜的影响。但是的研究仅限于处于静态平衡时的油楔，而旋转机械的滑动轴承油楔是在轴颈旋转时才能形成的，故没能考虑动态油楔的角偏斜现象，这与实际处境误差较大。

1.2国内研究概况

随着国内工业化进程的快速进展，国内动压滑动轴承在好多方面也取得了不小成就。

朱均等[10-11]利用一种类似于有限元法的计算方法，研究了轴承动力学系数，虽然这种方法计算精度可以达成要求精度，且所用计算时间比有限元法缩短好多倍，但是对于轴承动力学的计算结果和实际值相差较大。因此，这种方法没能在动压滑动轴承上有效利用，需要进一步提升；张勇斌、温诗铸等[12]在使用有限元法对油膜压力举行求解时，为了提高求解精度和速度，首次使用了多重网格算法。通过这种方法，研究了热效应与轴承性能之间的关系，进而得出了油膜特性与热效应之间的关系。但是，有限元法使用时需要输入的初始参数较大；姜歌东、谢友柏等[13]在针对转子动力系统的研究中，使用了自己提出的一种新的利用时域多工况识别方法来对滑动轴承油膜动特性举行研究，并且论证了这种方法的可行性与稳当性。但是，在实际处境中，汽轮机轴系的运行受大量因素影响，而文中未举行分析；赵文芳等[14]针对可倾瓦滑动轴承中存在的金属材料问题举行研究，系统总结了国内外轴瓦的多方面研究现状，得出了对可倾瓦轴承的研究中不成疏忽轴瓦金属材料与轴瓦粗糙度的影响。重点综述了可倾瓦轴承的金属材料，然而对于轴承油楔中油膜运动处境没有提及；吕延军等[15]对Reynolds方程举行修正时使用了变分约束原理，通过计算修正后Reynolds方程来研究瓦块的非线性油膜力。使用了八节点等参有限元法研究了单瓦块油膜力，通过叠加得到了整周油膜压力。通过计算Poincar映射和Runge-Kutta方法，研究了轴瓦支点对刚性转子系统的不平衡响应的影响。但是，对于轴承性能的研究仅限于单个瓦块的计算，叠加结果与实际值有确定误差；纪峰等[16]研究了实际生产中汽轮机可倾瓦轴承，在举行动态特性的模型计算时疏忽了瓦块变形，进而得出在研究可倾瓦轴承动特性时理应考虑瓦块变形，而不理应考虑瓦块的摇摆频率是否与转子转动频率一致。

1.3数值模拟研究概况

对于可倾瓦滑动轴承油膜的研究，轴瓦间隙、油膜间的剪切应力等通常都不被考虑。传统的只计算Reynolds方程或者计算分析轴承理论和非线性油膜力的方法已经不能得志实际生产需要，由于这两种方法产生的误差较大。

20世纪初，我国航天航空工业飞速进展，同时也带动了偏微分方程理论、数值计算方法、网格生成、计算机硬件工业等的空前进展，以及计算流体动力学CFD的逐步成熟。计算流体力学相比于测验流体力学和理论流体力学有大量优点：研究问题时，假设条件较少，可以模拟繁杂流场，应用广阔；时间短，效率高，效益好；可以考虑大量其他影响因素举行研究，应用范围很广。因此，在对于油膜特性的研究时，大多研究者采用计算流体力学[17]。

高庆水等[18]在分析圆柱轴承压力特性时首次使用Fluent数值模拟，论证了各种方法计算结果与实际的误差大小，但是没能给出一个研究其他滑动轴承的最优数值模拟方法；涂林等[19]在使用Fluent中的RNGk-ε模型模拟油膜特性时修正了湍动黏度，使模拟结果更接近实际值，但是对于小间隙的油膜流场计算误差较大。张楚等[20]基于气液两相流原理模拟计算油膜特性，得出了滑动轴承油膜分布规律，但是该方法不能很好的计算流体本身的滚动特性；而文献[21]使用单个模型对滑动轴承油膜举行数值计算，得出非线性油膜力与轴系的关系。但是对于轴承的间隙比没有很好考虑，因而研究的结果没有很精确的计算出实际应用当中的油膜特性。

2CFD中数学模型的选取

在流体运动中，当流体惯性力和粘性力在一个数量级时，湍流就产生了。汽轮机轴瓦进油口进油为压力进油，流体是三维的随机扰动流体，必然会发生湍流现象。对模型举行一系列假设，将理论与阅历结合，建立包括雷诺平均方程和脉动方程并且用来描写湍流平均量的封闭方程组，将该方程组称为湍流模型。

CFD中有大量湍流模型，包括Spalart-Allmaras模型、k-ε模型、k-ω模型、雷诺应力模型（RSM）、大涡模拟模型（LES）。其中，Spalart-Allmaras模型是相对简朴的单方程模型，不能用于繁杂的流场；雷诺应力模型（RSM）适用于雷诺应力明显具有各向异性的流场中，譬如龙卷风、燃烧室等带有猛烈旋转的流场；大涡模拟模型（LES）适用于大涡布局受流场影响较大的流场中。

k-ε模型假定湍流为各向同性的平匀湍流，不适用于汽轮机滑动轴承润滑油的运动。由于汽轮机滑动轴承进油为压力进油，对轴颈有确定冲击作用。并且当轴颈高速旋转，润滑油在轴颈和轴瓦组成的狭小间隙内做粘性剪切湍流运动。在这种存在冲击和不完全湍流的处境下，使用k-ε模型的结果会和实际结果悬殊很大。对k-ε模型举行一些提升可以得到RNAk-ε模型，它在计算功能上大大强于k-ε模型，同样由于油膜的冲击和涡旋的存在，RNAk-ε模型也不能用于汽轮机可倾瓦滑动轴承油膜特性的研究。

而k-ω模型包含了流体的剪切流分散，剪切应力输运k-ω模型简称SST（ShearStressTransport）k-ω湍流模型，它是将k-ω模型和k-ε模型都乘以一个混合函数，然后再相加的结果。这种模型具有k-ω模型近壁区计算和k-ε模型远场计算两大优点。SSTk-ω模型增加了横向耗散导数项，考虑了湍流剪切应力，使用的湍流常数也和其他模型不同。这使得SSTk-ω模型应用范围更加广泛，它可以用于带逆压梯度的滚动计算，并且能够很合理的考虑流体的涡流黏度。ShearStressTransport（SST）湍流模型对于研究汽轮机滑动轴承油膜特性的模拟计算更适合。

3结论与建议

1）在对油膜模型举行计算时，根据可倾瓦油楔中油膜的运动环境，可用ShearStressTransport（SST）湍流模型来计算，由于该种模型可以考虑到不成疏忽的狭小间隙内的油膜剪切应力，该模型对于油膜特性的模拟特别接近实际油膜特性。

2）对于可倾瓦油膜特性的研究，务必考虑到油膜温度和压力的关系，将两者耦合计算会使结果更加接近实际油膜流场。

3）可倾瓦滑动轴承的油膜压力、油膜温度、油膜厚度等都关系到汽轮机转子旋转的稳定性，虽然有关研究者做过圆形轴瓦、椭圆轴瓦、三油楔轴瓦轴承油膜对转子稳定性的研究，但是用N-S方程数值模拟计算对于可倾瓦的研究甚少。

[1]毕世华，李乃宏，郑钢铁.油膜轴承动态特性参数及转子不平衡的现场识别[J].航空学报，1998，9（3）：378-383.

[2]王致清.粘性流体力学[M].哈尔滨：哈尔滨工业大学出版社，1990.

[3]SeongheonYang，ChaesilKim，WonchangLee.PreventionofFlutteringFatigueDamageinaTiltingPadJournalBearing[J].TribologyInternational，2022，42：816-822.

[4]MatthewCha，PatrikIsaksson，SergeiGlavatskih.InfluenceofPadComplianceonNonlinearDynamicCharacteristicsofTiltingDadJournalBearings[J].TribologyInternational，2022，57：46-53.

[5]Kyung-BoBang，Jeong-HunKima，Yong-JooCho.ComparisonofPowerLossandPadTemperatureforLeadingEdgeGrooveTiltingPadJournalBearingsandConventionalTiltingPadJournalBearings[J].TribologyInternational，2022，43：1287-1293.

[6]White，M.F.andChanS.H.TheSubSynchronousDynamicBehaviorofTilting-PadJournalBearings[J].ASMEJ.ofTribology，1992，11（4）：167-173.

[7]ChangQYangP，MengY，WenS.ThermoElasticHydroDynamicanalysisoftheStaticPerfomanceofTilting-PadJournalBearingswiththeNewton-Raphsonmethod[J].TribologyInternational，2022，35（4）：225-234.

[8]EI-ButchAM，AshotirNM.TransientanalysisofMisalignedElasticTilting-PadJournalBearing[J].TribologyInternational，2022，38（1）：41-48.

[9]ZaehariadisDC.AXialWedgeEffectinTiltedHydrodynamicJournalBearings[J].TribologyTransactions，2022，51（3）：173-181.

[10]李小江，朱均.300MW汽轮发电机组径向可倾瓦轴承性能研究[J].润滑与密封，1996（1）：29-32.

[11]王凤才，徐华，朱均，等.大型动静力径向可倾瓦轴承热弹流体润滑研究[J].中国机械工程，2022，12（2）：140-143.

[12