分块瓦导轴承双向流固耦合动力特性数值模拟

王思淼　王文全　闫妍　王秀

摘要：基于双向流固耦合分析方法，对某混流式水轮发电机上导轴承进行动力学特性分析。讨论轴瓦数、轴承间隙和主轴偏心距离对轴瓦和主轴应力和形变量的影响。研究发现：适当地增加导轴承瓦块个数，可以使单个瓦片所受到的峰值应力降低，形变量减小，当轴瓦从4块增加到8块时，应力峰值可以降低31%，位移峰值降低33.5%;轴瓦应力对主轴与轴瓦之间的轴承油膜间隙厚度变化不敏感;主轴偏心距离增大，会使轴瓦的应力位移显著增加，增长趋势呈线性增长。主轴发生偏心时，轴瓦中心部位主要承载主轴偏心时的径向力。

关键词：导轴承;轴瓦数;轴承间隙;主轴偏心距;应力;流固耦合

中图分类号：TK730.3文献标志码：A

水力发电是可再生能源的重要组成部分，约占全球发电量的20%[1]。水力发电具有产电量大，产能稳定等优点，在能源供给中占据十分重要的位置。我国山川河流较多，水头落差大，水能资源丰富[2-3]，截至2020年底，我国常规水电装机达到3.38亿kW，年发电量1.35万亿kW·h，在建规模约4 800万kW;抽水蓄能装机3 149万kW，在建规模5 373万kW[4]。我国目前机电排灌保有量居世界第一，但仍有巨大的水力资源待开发利用[5-7]。

水力发电机组是水力发电站的核心动力机械设备，能够将水流冲击的能量转换为旋转机械能，通过主轴带动发电机转子转动产生电能。而水轮发电机组作为大型旋转机械，在运行中存在振动和摆动问题。水轮发电机组轴系的稳定性和振动特性是反映机组运行质量的重要指标，同时会影响到机组的安全运行和使用寿命[8]。过大的主轴振动会引起机组动静碰摩、加速零部件的磨损、造成构件疲劳损伤等故障，当机组的振动非常严重时，发电站轻则停机停产，重则毁坏机组，造成严重的破坏性事故，最终影响电网的安全运行[9]。所以，保证水轮发电机组安全稳定地运行是当前研究的重要课题。目前，水力发电电站大多是采用机组竖直安装的方式，即主轴系部分为铅垂方向安装，整个轴系包括主轴、安装在主轴上的转子、转轮等部件，工作状态下主轴所受到的外界载荷均由主轴上的导轴承承担[10]。水电站所使用的大部分机组均采用三导轴承布置，即只有上导、下导和水导3个导轴承，用来约束及承载整个轴系的工作载荷。导轴承则需要承担主轴上的全部径向载荷，导轴承的动力特性是关乎主轴轴系安全的关键因素之一[9，11-12]。所以，确保轴承拥有足够的承载能力和强度是导轴承研究的重要内容。

目前，国内外对导轴承的研究主要以轴瓦静力分析或单独分析油膜压力为主。杨晓明、马震岳等[13]采用有限元方法分析了立式可倾瓦导轴承的动力特性;张雷克等[14]讨论了密封间隙对轴系的影响;张金剑等[15]分析了质量偏心对轴系运动的影响。导轴承对主轴轴系的约束作用主要取决于导轴承的几何尺寸，而本文中采用双向流固耦合动力学计算方法，讨论轴瓦数、油膜间隙厚度和偏心距对轴瓦和轴颈应力分布的影响，得出较为贴近真实工况下的数值结果，为水轮发电机导轴承的优化设计提供有益的参考。

1控制方程及数值模型建立

1.1流固耦合控制方程

1.2导轴承模型的建立

本次计算所建立的模型参考立式水轮发电机组（水轮机组HLA551-LJ-43，发电机SF55-10/740），导轴承参数参考机组的上导轴承尺寸。为模拟主轴运转过程中发生偏心时轴瓦的受力情况，模型中轴系转动部分只取主轴部分长度，研究轴段全长取轴承高度的11倍长度，轴承位置设置于主轴所取轴段部分的中点处。

分块瓦油润滑导轴承三维模型结构和轴承瓦片布置结构见图1。

计算所用导轴承模型尺寸如下：轴瓦包角为40°，轴瓦瓦片厚度15 mm，支撑的楔板高度为4 mm，主轴与导轴承接触的轴颈部分直径100 mm，主轴直径90 mm，导轴承高度为80 mm。

导轴承润滑所使用的油膜材料参数设置：润滑油密度为960 kg/m³，动力粘度为0.048 Pa·s，不考虑温度变化。

為后续结果显示方便，这里对轴瓦进行编号，如图1（b）所示，以图示位置左边轴瓦编号为1，顺时针方向旋转编号，轴瓦编号依次为2、3、4、…，所有轴瓦沿圆周方向均匀分布安装，后文提到偏心方向均为沿半径指向1号瓦中心，图中表示瓦片均用“WA”字母表示。

1.3模型假设和边界条件

1）模型假设

根据实际模拟条件限制，对所分析的导轴承进行适当的模型简化假设。

（1）轴瓦楔板支撑位置假设为与楔板固结，即瓦片为固定瓦;

（2）油膜假设为牛顿液体，动力粘度沿油膜厚度方向稳定;

（3）油膜与主轴接触的壁面之间、油膜表面与主轴表面没有相对位移，不考虑由于摩擦而产生的热能。

2）边界条件

（1）支撑轴瓦的楔板固定端表面设置为混合固定边界，即该表面不存在变形与位移;

（2）油膜内表面设置为旋转边界，旋转速度与轴颈转速均为600 rpm;

（3）油膜除内表面外的其他表面设置为固定无滑移壁面;

（4）对偏心方式进行简化处理，偏心运动距离施加于研究轴段的上下两端面，模拟运转总时长设置为0.4 s，0～0.1 s主轴运转期间不发生偏心，0.1～0.2 s期间主轴轴线向1号轴瓦方向发生线性径向偏心，0.2～0.3 s期间主轴于设置的偏心距离位置运行，0.3～0.4 s期间主轴偏心线性回归与轴承同心位置。

3）模型网格数

网格划分采取分区域划分，油膜部分切分为圆环部分与轴瓦间油膜部分。圆环部分为主轴与轴瓦内面之间间隙区域。由于这部分几何尺寸很小，故单独划分网格，划分后网格单元数为378 840。固体部分网格单元数为635 026。模型总网格数为1 101 946。E5A545AB-5E52-4087-8057-AFF7EEBA172A

1.4计算模拟参数

改变导轴承3个主要参数进行模拟计算，分别是：瓦片数、油膜间隙厚度、不同主轴偏心距。对于瓦片数，选取5种瓦片数进行计算模拟（4块瓦片、5块瓦片、6块瓦片、7块瓦片、8块瓦片）;对于油膜间隙厚度，選取间隙分别为0.12、0.15、0.18 mm进行对比计算。根据瓦片数和油膜间隙的结果进行对比分析，最后选取油膜间隙0.15 mm，6块瓦片数的组合参数进行5种不同偏心距离的计算模拟。

2结果与讨论

2.1不同轴瓦数对导轴承动力特性的影响

模拟不同瓦块数时，油膜间隙取用0.12 mm。由于主轴发生偏心时，会对偏心方向的轴瓦产生最大压力，故只取最危险瓦片（即1号瓦片）结果为评断标准。图2为轴承瓦片数不同的情况下，主轴发生偏心时1号瓦片的应力大小。可以看出，主轴产生相同偏心距的情况下，轴承瓦片增加，可以使1号瓦片的最大应力值和平均应力值下降;当瓦片数从4增加到5、6、7、8时，相对于瓦片数为4的1号瓦片的应力峰值分别下降了4.71%、11.82%、20.96%、31.14%，应力均值分别下降了3.41%、8.78%、12.74%、19.16%。 应力峰值的显著下降可以使轴瓦在工作状态下更安全稳定。

图3是不同瓦片数情况下，1号瓦片受力发生位移的最大值和平均值。从图中可以看出，从5块瓦片增加到6块瓦片时，可以显著降低危险瓦片的最大位移值，但6块瓦片增加为7片时，形变量减小相对不明显;当瓦片数从4增加到5、6、7、8时，相对于瓦片数为4的1号瓦片的位移峰值分别下降了5.82%、19.17%、22.50%、33.50%，位移均值分别下降了12.33%、24.24%、29.26%、36.86%。

2.2油膜间隙厚度对导轴承动力特性的影响

关于油膜间隙这一参数，对不同的瓦块数都进行3种油膜间隙（0.12、0.15、0.18 mm）的模拟分析。与上文2.1一样，1号瓦片受最大压力，为最危险瓦片，所以以1号瓦片峰值应力和轴瓦最大变形量作为参考标准进行分析，结果如表1所示。由表1和表2中结果可以看出：油膜间隙的变化对轴瓦所受应力和变形影响较小;但轴瓦数目变多，会使整体偏心产生的力分解为多个分力，减小单块瓦片所受的应力峰值，有利于防止轴瓦因应力过大发生破坏。

2.3主轴发生不同偏心距离对导轴承动力特性的影响为探究不同偏心距离对轴瓦应力变形的影响，对油膜间隙0.15 mm、瓦块数为6的轴承模型，进行5个不同的偏心参数的模拟计算，偏心距离分别为0.01、0.02、0.03、0.04、0.05 mm。结果中提取6块瓦所受应力的最大值和平均值（图4），形变的最大值和平均值（图5）。

由图4可以看出，偏心距离越大，对轴瓦的压力会变大，相对于偏心距0.01 mm的应力峰值2.016 8 MPa，当偏心距增加至0.02、0.03、0.04、0.05 mm时，应力峰值分别增加了97.28%、195.25%、294.11%、397.82%，应力峰值增加量随偏心距离的增加大致呈线性增长;偏心距为0.01 mm时应力均值为0.355 051 MPa，偏心距增加后应力均值增长量为100.67%、192.92%、302.07%、405.15%，同样呈线性增长。如图5，随着主轴偏心距离的增大，轴瓦的位移也随之增加，偏心距为0.01 mm时最大位移为0.000 134 mm，偏心距增加至0.02、0.03、0.04、0.05 mm时，相对于0.01 mm位移最大值分别增加98.36%、178.60%、295.84%、397.26%。

图6为1号瓦片在主轴偏心0.01、0.03、0.05 mm时的应力云图对比图。从图中可以看出，轴瓦受到应力最大值主要集中分布在轴瓦中心的两侧位置处，靠近瓦面边缘和中心的位置受力很小。图7为轴颈表面在主轴偏心0.01、0.03、0.05 mm时的应力云图对比图。可以看出，轴瓦所对应的轴颈位置应力较大，轴瓦之间空隙位置对应的轴颈区域应力相对较小。

3结论

基于双向流固耦合分析方法，分析了轴瓦数、油膜间隙厚度、主轴偏心距离对轴瓦和轴颈应力和形变量的影响规律。主要结论如下：

1）当主轴发生偏心、对轴瓦产生过大的压力时，可以考虑采用适当增加瓦块数的方式来改善单个瓦块所承载的应力峰值。

2）轴瓦应力对油膜间隙厚度变化不敏感。

3）主轴偏心距离增大，会对轴瓦局部造成较大压力;偏心越大，会使轴瓦的形变峰值显著增加，应预防轴瓦变形过大而引起破坏。

4）主轴轴系发生偏心时，主要对轴瓦的中心产生较大的压力，对轴瓦边缘的压力很小，轴瓦中心主要承载主轴偏心时的径向力。参考文献：

[1]SAWIN J L， MARTINOT E. REN21 Global Status Report finds solid growth in renewables[J]. Power Engineering International， 2011， 19（8）： 26-30.

[2] 邢利钧. 清洁能源开发利用方式与储能技术展望[J]. 科技资讯， 2021， 19（1）： 56-58.

[3] 曾庆捷. 我国未来储能技术发展探讨[J]. 山东电力技术， 2017， 44（3）： 51-53.

[4] 《中国可再生能源发展报告》正式发布[J]. 建设机械技术与管理， 2021， 34（4）： 7.

[5] 孙洁， 蔡建国， 葛新峰， 等. 基于正交试验多目标的两叶片灯泡贯流式水轮机性能优化[J/OL]. 中国电机工程学报： 1-10[2021-12-15]. https：//kns.cnki.net/kcms/detail/detail.aspx？FileName=ZGDC2021081300K&DbName=DKFX2021.E5A545AB-5E52-4087-8057-AFF7EEBA172A

[6] ZHOU D Q，DENG Z Q.Ultra-low-head hydroelectric technology：a review[J].Renewable and Sustainable Energy Reviews，2017，78：23-30.

[7] CHEN H X， ZHOU D Q， KAN K， et al. Transient characteristics during the co-closing guide vanes and runner blades of a bulb turbine in load rejection process[J]. Renewable Energy， 2021， 165： 28-41.

[8] 彭悅蓉， 赖喜德， 张惟斌， 等. 考虑轴承及密封的水电机组轴系特性分析[J]. 中国农村水利水电， 2015， 6： 153-157.

[9] 彭悦蓉. 导轴承的液膜参数对立式水电机组转子动力特性影响[D].成都：西华大学， 2016.

[10]徐鹰健. 大型立式电机上导轴承数值模拟和工程实践[D]. 苏州：苏州大学， 2014.

[11]KANG Y， CHANG Y P， TSAI J W， et al. An investigation in stiffness effects on dynamics of rotor-bearing-foundation systems[J]. Journal of Sound and Vibration， 2000， 231（2）： 343-374.

[12]王正伟， 喻疆， 方源， 等. 大型水轮发电机组转子动力学特性分析[J]. 水力发电学报， 2005（4）： 62-66.

[13]杨晓明，马震岳，黄军义.导轴承与推力轴承耦合作用下水电机组的横向振动研究[J]. 水力发电学报， 2007（6）： 132-136.

[14]张雷克， 张金剑， 王雪妮， 等. 水电机组转子-转轮弯扭耦合振动特性分析[J]. 水力发电学报， 2021， 40（9）： 102-112.

[15]张金剑， 张雷克， 吴嵌嵌， 等. 动静偏心电磁激励下水电机组碰摩转子-轴承系统弯扭耦合振动特性分析[J]. 振动与冲击， 2021， 40（12）： 236-245.

（责任编辑：曾晶）

Numerical Simulation of Dynamic Characteristics of Bidirectional

Fluid-solid Coupling of Block-segment Guide Bearing

WANG Simiao WANG Wenquan YAN Yan WANG Xiu

（1.Faculty of Civil Engineering and Mechanics， Kunming University of Science and Technology， Kunming 650500， China;

2.School of Water Resources and Hydropower， Sichuan University， Chengdu 610065， China;

3.State Key Laboratory of Hydraulics and Mountain River Engineering， Sichuan University， Chengdu 610065， China）Abstract： Based on the bidirectional fluid-solid coupling analysis method， the dynamic characteristics of the guide bearing of a Francis hydro-generator set are analyzed. The effects of the bearing bush number， the bearing clearance and the eccentricity of the spindle on stress and deformation of the bearing bush and the spindle are discussed. The study found that appropriately increasing the number of guide bearing bushes can reduce the peak stress and deformation of a single bearing bush. When the bearing bushes are increased from 4 to 8， the peak stress can be reduced by 31% and the peak displacement can be reduced by 33.5%; the stress of the bearing bush is insensitive to the thickness change of bearing clearance between the spindle and the bearing bush; the increase in the eccentricity of the spindle will cause linear growth in local pressure on the bearing bush. When the spindle is eccentric， the center of the bearing bush mainly bears the radial force.

Key words： guide bearing; number of bearing bush; bearing clearance; spindle eccentricity; stress; fluid-solid couplingE5A545AB-5E52-4087-8057-AFF7EEBA172A