尹庭赟，裴 吉，袁寿其，王文杰



余热排出泵叶轮流固耦合特性分析

尹庭赟，裴 吉※，袁寿其，王文杰

（江苏大学国家水泵及系统工程技术研究中心，镇江 212013）

余热排出泵长期在高温高压环境下运行，其结构的可靠性对整个机组的安全运行有着重要的影响。该文采用单向流固耦合方法研究了几种常用材料及增加叶轮盖板厚度对转子结构动力学特性的影响，同时对比了常温和高温下叶轮的振动特性。结果表明，余热排出泵高温运行时，叶轮的第1阶固有频率为394.17 Hz，较常温下升高了2.28%，远高于叶轮叶片通过频率。在各阶模态振型下，1Cr13MoS叶轮固有频率最大，ZG225-450叶轮固有频率最小。设计工况和偏工况下叶轮的应力及变形变化趋势基本一致，叶轮的变形随半径的增加而不断增大，最大变形量出现在后盖板叶轮出口处。叶轮的应力分布不均，最大应力均出现在叶片尾缘与后盖板接合处。沿着前盖板和叶片接合线，应力的峰值随前盖板厚度增加而减小。沿着后盖板和叶片接合线，增加后盖板厚度明显减小了应力峰值。高温和常温应力变化趋势基本一致，但是高温下叶轮应力明显高于常温。尤其在叶轮进口处附近，高温应力相对于常温应力增加超过300%。该研究结果为下一步进行更加复杂的动力学分析、疲劳分析以及结构优化提供参考。

泵；叶轮；动力学；余热排出泵；盖板厚度

0 引 言

随着中国对核电的高度重视，国内高校和泵企业等正加紧研制开发核电用泵，实现核电用泵的国产化。在核电站内的泵类设备中，余热排出泵是重要的设备之一，可保护核电站的安全运行：当反应堆进行正常停堆时，排出堆芯余热和核主泵在一回路中产生的热量；当发生严重事故时，减小核电站受损程度[1-2]。由于余热排出泵高温高压运行环境的特殊性，其结构的可靠性对整个机组的安全运行有着重要的影响，因此，准确地对余热排出泵转子部件进行强度分析可为其结构设计及寿命延长提供依据。

近年来，流固耦合有限元分析已成为分析流体机械可靠性的重要方法。国内外不少学者尝试采用流固耦合的方法对转子系统进行振动特性分析，取得了许多有价值的成果。刘厚林等[3]分别运用单、双向流固耦合方法对导叶式离心泵叶轮进行静应力强度分析，指出单向耦合即可满足叶轮静力学分析。Kan等[4]采用单向流固耦合计算得到了水轮机转子在多工况下应力和变形变化规律。牟介刚等[5]采用单向流固耦合方法，分析了悬臂式离心泵在不同悬臂比、不同工况下叶轮应力和变形的变化规律，并对不同悬臂比离心泵的转子动力学特性进行研究。Schneider等[6]对多级离心泵进行了单向流固耦合计算，分析了叶轮结构设计参数和温度对叶轮应力和变形的影响。刘厚林等[7]对比分析了无预应力下、流固耦合作用力和旋转离心力作用下的余热排出泵转子的模态。刘建瑞等[8]考虑了水的附加质量对余热排出泵转子振动性能的影响，对转子进行水中模态分析。Zhou等[9-10]采用双向流固耦合计算方法，研究了考虑流固耦合作用的余热排出泵在不同工况下内部流场及转子应力和变形分布。考虑到余热排出泵实际高温高压的运行环境，目前很少有学者考虑对叶轮盖板进行加厚处理。同时，材料属性对转子动力学特性的影响研究较少，尤其是对于有特殊用途的余热排出泵而言，合理选择叶轮材料至关重要。

因此，本文应用单向流固耦合计算方法研究了余热排出泵实际高温运行环境下叶轮模态分布，详细分析了叶轮几种常用材料和不同盖板厚度下叶轮变形量和应力分布，以期为研究整个系统的动态特性提供参考。

1 流固耦合动力学方程

1.1 基于模态分析的静流体整体耦合法

对于具有个自由度的系统，其结构动力方程为

式中是系统的质量矩阵；是系统的阻尼矩阵；是系统的刚度矩阵；是系统承受的载荷，Pa；是节点的加速度矢量，m/s2；是节点的速度矢量，m/s；是节点的位移矢量，m[11]。

在静流体中，考虑到流固耦合作用对方程做了一些假设：流体是可压缩的，压力波动将会引起密度变化；流体是非黏性的（即无黏性耗散）以及非旋转的；流体没有运动。因此，流体将被视为声学流体，控制方程可用如下的声波方程表示

（3）

式中是声音在水中的传播速度，m/s；是由声波引起的流体压力以及作用于流体的力，Pa；Δ是拉普拉斯算子；是流体体积模量；是时间，s；0是流体平均密度，kg/m3[11]。

在流固耦合交界面，用如下的方程描述了流体与结构的相互作用

（5）

考虑到在交界面处压力对结构的影响，方程（1）改写为

其中，= -表示等效耦合刚度[11]。

因此，可以耦合流体和结构方程，基于单一求解器通过整体求解方法同时求解[11]，方程表达如下

在对叶轮进行模态分析时，自由振动方程系统视为非阻尼情况[11]，方程（7）简化如下

（8）

1.2 基于应力和变形分析的动流体分区耦合法

在分区耦合法中，有2种耦合求解流体动力学和结构响应策略：单向耦合和双向耦合[12]。双向耦合方法主要是针对具有强物理耦合效应的问题，即除了考虑流动对结构的影响，结构的变形或运动对流动的反作用也需要考虑[13]。在泵机械中，大多数情况下由于泵内不稳定流引起的叶轮的变形很小，而叶轮变形和振动对流场影响可以忽略不计[10,14]。因此，本文采用单向耦合方法，如图1所示，在每一个时间步长流体载荷通过交界面单向传递到结构。

2 模型基本参数及边界条件

2.1 计算模型

余热排出泵过流部件包括后掠扭曲叶轮、径向导叶和环形蜗壳3部分，其三维结构如图2所示。余热排出泵的主要性能和几何参数分别为设计流量d=910 m3/h，扬程d=77 m，转速=1 490 r/min，比转数s=104.5，叶轮进口直径j=270 mm，叶轮出口直径2=511 mm，叶轮叶片数i=5，导叶进口直径3=515 mm，导叶出口直径4=718 mm，导叶叶片数g=7，蜗壳基圆直径5=840 mm。为探讨增加盖板厚度对转子结构动力学特性的影响，设计了4种不同的叶轮结构：1）原始盖板；2）前盖板厚度增加50%；3）后盖板厚度增加50%；4）前后盖板厚度均增加50%，详细参数见表1。余热排出泵长期在高温高压环境恶劣环境下运行，叶轮材料的选取至关重要。因此本文选取5种常用材料[15]如表2所示，分析了不同材料叶轮的动力学特性。

2.2 网格划分

2.2.1 流体域网格

网格划分是数值模拟中的重要步骤之一，网格的好坏直接影响数值模拟的准确性。对余热排出泵计算域采用ICEM软件进行结构化网格划分，对壁面进行网格加密[16-17]，如图3所示。对数值模拟进行网格无关性检查，如图4所示，确定余热排出泵计算域的最终网格总数为5.20×106。

表1 叶轮结构参数

表2 叶轮材料参数

2.2.2 结构域网格

为研究湿模态下叶轮的动力学特性，将叶轮完全淹没在1.12 m×0.36 m×1.12 m的水体中[18]。如图5所示，运用ANSYS Workbench中自带的网格划分软件，对叶轮及水体部分进行网格划分[19]。分别基于Method和Hex Dominant网格类型通过Sizing命令对叶轮和水体域进行网格划分。经网格无关性检查，确定最终网格节点总数为129 220，网格单元总数为49 556。

2.3 边界条件

采用整体耦合法对叶轮进行模态分析时，为了获得更高的求解精度，用Fluid220和221高阶单元作为水体域单元，默认使用Solid186和187作为结构域单元。由于水体域单元不支持完全积分法，所以本文对水体域积分控制使用缩减积分法，而对流固耦合面的固体域积分控制使用完全积分法。考虑到湿模态分析时，构成结构动力学有限元方程的质量矩阵是非对称的，因此在求解中需选用非对称模态算法[20]。

假设流体在泵内部是不可压缩三维湍流流动，满足连续方程和动量方程，采用SST湍流模型求解雷诺时均方程[21-24]。采用分区耦合法对叶轮进行应力和变形量分析时，先对流体域进行了三维定常和非定常数值模拟：在定常计算过程中，边界条件设置为总压进口和质量流量出口，动静域之间交界面采用Frozen rotor，计算域壁面采用无滑移网格函数[25-26]，采用高阶求解精度，收敛残差RMS设置为10-5，计算迭代步数最大500；在非定常数值计算过程中，以定常计算结果作为初始值，边界条件设置不变，而动静域之间交界面改成Transient frozen rotor[27]。计算周期为6圈，叶轮旋转1度为1个时间步长，即5.747 13×10-5s。以非定常计算最后一步的计算结果作为初值对叶轮进行了6圈单向流固耦合计算：时间步长设置和非定常计算一致；流体场和固体场叶轮叶片、前后盖板设置为流固耦合交界面（fluid-structure interface）；有限元时间积分采用采用Hilber-Hughes- Taylor（HTT）离散方法，使用瑞利阻尼（Rayleigh damping）系数来定义结构瞬态阻尼效应[28]。如图1所示，设置叶轮后盖板的孔面为Fixed Support支撑面。本文分别模拟了设计工况下25 ℃常温水以及180 ℃高温水余热排出泵内部流动以对比实际运行环境下叶轮的动力学特性。

2.4 试验验证

为了验证数值模拟的准确性，对余热排出泵按比例缩小30%的模型泵进行了外特性试验验证。试验在江苏大学流体机械质量技术检验中心的开式试验台上进行，试验现场如图6所示。采用上海自仪九仪表有限公司生产的型号规格为LWGY-200A涡轮流量计进行测量流量，量程为0～1.6 MPa，精确度为0.5级；采用上海威尔泰仪器仪表有限公司生产的WT200智能压力变送器测量进出口压力，进出口处压力变送器的测量量程分别为−0.1～0.1 MPa和0～1.6 MPa，测量精度均为0.1级。如图7所示，根据相似换算原理[29-30]，得到了真实尺寸泵的试验和数值模拟性能曲线。由图7可知，试验扬程和效率略高于计算值，预测扬程的最大偏差为6.78%，预测效率的最大偏差为9%。但是试验结果趋势一致，计算模型较准确地预测了余热排出泵的外特性，数值计算结果较为可信。

3 结果及分析

3.1 模态分析

为了了解叶轮模态振型情况，首先对原始盖板的合金结构钢叶轮的计算结果提取分析。如图8所示，叶轮第1、2阶振型相似，都是在叶轮某2个对称位置的变形量最大，即叶轮是摆动变形。第3、4阶模态振型相似，都是叶轮沿着轴向的扭动变形。第5阶和第6阶模态振型相似，都是扭曲变形，可以看到2条明显的节径，且分布角度不同。第7、8、9、10阶振型为叶轮前后盖板沿着轴向的扭曲变形。通过对其他几种形式的叶轮前10阶模态振型分析发现，叶轮振型基本相同，但是固有频率有一定差别。

表3列出了不同材料、盖板厚度和温度下叶轮前10阶固有频率。由表3可知，改变叶轮盖板厚度对前6阶固有频率的影响可以忽略。从第7阶开始，叶轮后盖板的厚度对固有频率有很大影响，增加后盖板厚度可以提高叶轮的固有频率。5种材料的叶轮固有频率变化规律趋于一致。在常温各阶模态振型下，1Cr13MoS(B_3)叶轮固有频率最大，ZG225-450(B_2)叶轮固有频率最小。以合金结构钢(B_0)原始盖板叶轮为例，余热排出泵在高温运行时，叶轮的第1阶固有频率为394.17 Hz，较常温下升高了2.28%，远高于高于叶轮叶片通过频率，在安全运行范围内。

3.2 应力和变形分析

如图9、10所示分别为叶轮在3种工况下变形及应力分布，其中对变形量和应力进行了无量纲处理[7]，从图中可以看出在设计工况和偏工况下叶轮的应力及变形变化趋势基本一致。叶轮的变形随半径的增大而不断增大，最大总变形出现在后盖板叶轮出口处。叶轮的应力分布不均，最大应力均出现在叶片尾缘与后盖板接合处。

为了详细分析叶轮关键部位的应力变化，如图11定义了前后盖板和叶片接合处的4条路径，这4条路径在叶轮表面从进口前缘延伸到出口尾缘。

图12为不同工况下叶轮4条路径应力分布图。从图12中可以看出3种工况下4条路径应力变化趋势基本一致：在前后盖板与叶片压力面接合处，路径中间附近出现应力峰值，在前缘附近可以明显的看出大流量工况下应力大，小流量工况下应力小；在前后盖板与叶片吸力面接合处，最大应力均出现在尾缘处。其中，叶片吸力面与盖板接合处的应力峰值远高于与压力面接合处的应力峰值。

余热排出泵的实际运行环境对叶轮材料的选择提出了苛刻的要求，如图13所示为设计工况下5种叶轮材料的4条路径应力分布图。从图13中可以看出5种材料的4条路径应力变化趋势基本一致，但是在前缘附近合金结构钢(B_0)的应力明显高于其他4种材料。

如图14所示为设计工况下4种叶轮盖板厚度对应的4条路径应力分布图。从图14b和14d中可以看出，沿着前盖板和叶片接合线，应力的峰值随前盖板厚度增加而减小。从图14a和14c中可以看出，沿着后盖板和叶片接合线，增加后盖板厚度明显减小了应力峰值。

图15所示对比分析了设计工况下原始模型叶轮在常温和高温下4条路径应力。由图15可知，两种温度下四条路径应力变化趋势基本一致，但是高温下叶轮应力明显高于常温。尤其在叶轮进口处附近，高温应力相对于常温应力增加超过300%。在前盖板和叶片压力面接合处，常温时应力峰值在路径中间附近，而高温时峰值在尾缘。其余3种路径峰值位置没有改变。

4 结论与讨论

本文应用单向流固耦合计算方法分析了余热排出泵实际高温运行环境下叶轮动力学特性，并定量研究了叶轮不同结构设计参数和不同材料下的变形量和应力，结论如下：

1）余热排出泵高温运行时，叶轮的第1阶固有频率为394.17 Hz，较常温下升高了2.28%，远高于叶轮叶片通过频率。在各阶模态振型下，1Cr13MoS的叶轮固有频率最大，ZG225-450叶轮固有频率最小。改变叶轮盖板厚度对前6阶固有频率影响可以忽略。但是从第7阶开始，增加后盖板厚度可以提高叶轮的固有频率。

2）设计工况和偏工况下叶轮的应力及变形变化趋势基本一致。叶轮的变形随半径的增大而不断增大，最大总变形出现在后盖板叶轮出口处。叶轮的应力分布不均，最大应力均出现在叶片与后盖板接合的进口处。其中，沿着前盖板和叶片接合线，应力的峰值随前盖板厚度增加而减小。沿着后盖板和叶片接合线，增加后盖板厚度明显减小了应力峰值。

3）高温和常温应力变化趋势基本一致，但是高温下叶轮应力明显高于常温。尤其在叶轮进口处附近，高温应力相对于常温应力增加超过300%。在前盖板和叶片压力面接合处，常温时应力峰值在路径中间附近，而高温时峰值在尾缘。其余3种路径峰值位置没有改变。后文将在本文基础上深入分析高温高压条件下泵的转子动力学特性及叶轮使用寿命评估。

[1] Wang Wenjie, Pei Ji, Yuan Shouqi, et al. Application of different surrogate models on the optimization of centrifugal pump[J]. Journal of Mechanical Science and Technology, 2016, 30(2): 567－574.

[2] 王文杰，袁寿其，裴吉，等. 余热排出泵叶轮与导叶匹配的水力性能研究[J]. 华中科技大学学报：自然科学版，2015，43(7)：39－43.

Wang Wenjie, Yuan Shouqi, Pei Ji, et al. Hydraulic performance study on matching of impeller and diffuser of residual heat removal pump[J]. Journal of Huazhong University of Science & Technology: Natural Science Edition, 2015, 43(7): 39－43. (in Chinese with English abstract)

[3] 刘厚林，徐欢，吴贤芳，等. 基于流固耦合的导叶式离心泵强度分析[J]. 振动与冲击，2013，32(12)：27－30.

Liu Houlin, Xu Huan, Wu Xianfang, et al. Strength analysis of a diffuser pump based on fluid-structure interaction[J]. Journal of Vibration & Shock, 2013, 32(12): 27－30. (in Chinese with English abstract)

[4] Kan Kan, Zheng Yuan, Zhang Xin, et al. Numerical study on unidirectional fluid-solid coupling of Francis turbine runner[J]. Advances in Mechanical Engineering, 2015, 7(3): 1－14.

[5] 牟介刚，陈莹，谷云庆，等. 悬臂式离心泵流固耦合特性研究[J]. 哈尔滨工程大学学报，2016，37(8)：1111－1117.

Mou Jiegang, Chen Ying, Gu yunqing, et al. Research on fluid-structure interaction characteristics of cantilever centrifugal pump[J]. Journal of Harbin Engineering University, 2016, 37(8): 1111－1117. (in Chinese with English abstract)

[6] Schneider A, Will B C, Böhle M. Numerical evaluation of deformation and stress in impellers of multistage pumps by means of fluid structure interaction[C]//ASME 2013 Fluids Engineering Division Summer Meeting (FEDSM2013), 2013.

[7] 刘厚林，徐欢，王凯，等. 基于流固耦合的余热排出泵转

子模态分析[J]. 流体机械，2012(6)：28－32.

Liu Houlin, Xu Huan, Wang Kai, et al. Modal analysis for rotor of residual heat removal pump based on fluid-structure interaction[J]. Fluid Machinery, 2012(6): 28－32. (in Chinese with English abstract)

[8] 刘建瑞，陈斌，张金凤，等. 余热排出泵水中模态分析[J]. 排灌机械工程学报，2015，33(4)：290－295.

Liu Jianrui, Chen Bin, Zhang Jinfeng, et al. Mode analysis for rotor of residual heat removal pump in water[J]. Journal of Drainage and Irrigation Machinery Engineering (JDIME), 2015, 33(4): 290－295. (in Chinese with English abstract)

[9] Zhou Banglun, Yuan Jianping, Fu Yanxia, et al. Investigation of dynamic stress of rotor in residual heat removal pumps based on fluid-structure interaction[J]. Advances in Mechanical Engineering, 2016, 8(9): 1－12.

[10] Zhou Banglun, Yuan Jianping, Lu Jiaxing, et al. Numerical investigation of residual heat removal pumps based on fluid-structure interaction in 1000 MW nuclear power plants[C]//International Symposium on Fluid Machinery and Fluid Engineering. IET, 2014.

[11] Wu Yulin, Li Shengcai, Liu Shuhong, et al. Vibration of Hydraulic Machinery[M]. Germany: Springer Netherlands, 2015.

[12] Piperno S, Farhat C. Partitioned procedures for the transient solution of coupled aeroelastic problems-Part II: Energy transfer analysis and three-dimensional applications[J]. Computer Methods in Applied Mechanics and Engineering, 2001, 190(24): 3147－3170.

[13] Pei Ji, Benra F K, Dohmen H J. Application of different strategies of partitioned fluid-structure interaction simulation for a single-blade pump impeller[J]. Proceedings of the Institution of Mechanical Engineers, Part E: Journal of Process Mechanical Engineering, 2012, 226(4): 297－308.

[14] Yuan Shouqi, Pei Ji, Yuan Jianping. Numerical investigation on fluid structure interaction considering rotor deformation for a centrifugal pump[J]. Chinese Journal of Mechanical Engineering, 2011, 24(4): 539－545.

[15] Deng Xin, Yuan Shouqi, Yuan Jianping, et al. Antiseismic response research of horizontal residual heat removal pump in different seismic spectrum input directions[J]. Journal of Vibroengineering, 2014, 16(3): 1306－1317.

[16] 洪锋，袁建平，周帮伦. 余热排出泵叶轮内空化流动特性的数值分析[J]. 排灌机械工程学报，2016，34(3)：185－190.

Hong Feng, Yuan Jianping, Zhou Banglun. Numerical analysis of cavitating flow characteristics in impeller of residual heat removal pump[J]. Journal of Drainage and Irrigation Machinery Engineering (JDIME), 2016, 34(3): 185－190. (in Chinese with English abstract)

[17] 袁建平，朱钰雯，周帮伦，等. 余热排出泵内部压力脉动特性分析[J]. 排灌机械工程学报，2015，33(6)：475－480.

Yuan Jianping, Zhu Yuwen, Zhou Banglun, et al. Pressure fluctuation in residual heat removal pumps[J]. Journal of Drainage and Irrigation Machinery Engineering (JDIME), 2015, 33(6): 475－480. (in Chinese with English abstract)

[18] 施卫东，郭艳磊，张德胜，等. 大型潜水轴流泵转子部件湿模态数值模拟[J]. 农业工程学报，2013，29(24)：72－78.

Shi Weidong, Guo Yanlei, Zhang Desheng, et al. Numerical simulation on modal of large submersible axial-flow pump rotor[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(24): 72－78. (in Chinese with English abstract)

[19] 梁权伟，王正伟，方源. 考虑流固耦合的混流式水轮机转轮模态分析[J]. 水力发电学报，2004，23(3)：116－120.

Liang Quanwei, Wang Zhengwei, Fang Yuan. Modal analysis of Francis turbine with considering FSI[J]. Journal of Hydroelectrical Engineering, 2004, 23(3): 116－120. (in Chinese with English abstract)

[20] 李兵. ANSYS Workbench设计、仿真与优化[M]. 北京：清华大学出版社，2011.

[21] MENTER F R. Two-equation eddy-viscosity turbulence models for engineering applications[J]. AIAA journal, 1994, 32(8): 1598－1605.

[22] 任芸，刘厚林，舒敏骅，等. 湍流模型在离心泵偏工况性能预测中的适用性分析[J]. 流体机械，2012，40(10)：18－22.

Ren Yun, Liu Houlin, Shu Minhua, et al. Analysis of applicability of turbulence models in performance prediction for centrifugal pumps at the off-design conditions[J]. Fluid Machinery, 2012, 40(10): 18－22. (in Chinese with English abstract)

[23] 杨兴标，李红，陈超. 湍流模型对离心泵扬程预测准确性的影响[J]. 排灌机械工程学报，2015，33(8)：656－660.

Yang Xingbiao, Li Hong, Chen Chao. Influence of turbulence models on head prediction accuracy for centrifugal pumps[J]. Journal of Drainage and Irrigation Machinery Engineering (JDIME), 2015, 33(8): 656－660. (in Chinese with English abstract)

[24] 王福军. 流体机械旋转湍流计算模型研究进展[J]. 农业机械学报，2016，47(2)：1－14.

Wang Fujun. Research progress of computational model for rotating turbulent flow in fluid machinery[J]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(2): 1－14. (in Chinese with English abstract)

[25] 朱相源，江伟，李国君，等. 导叶式离心泵内部流动特性数值模拟[J]. 农业机械学报，2016，47(6)：34－41.

Zhu Xiangyuan, Jiang Wei, Li Guojun, et al. Numerical analysis of hydraulic performance in centrifugal pump with vane diffuser[J]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(6): 34－41. (in Chinese with English abstract)

[26] 孟根其其格，谭磊，曹树良，等. 离心泵蜗壳内非定常流动特性的数值模拟及分析[J]. 机械工程学报，2015，51(22)：183－190.

Meng Genqiqige, Tan Lei, Cao Shuliang, et al. Numerical simulation and analysis of unsteady flow characteristics in centrifugal pump volute[J]. Journal of Mechanical Engineering. 2015, 51(22): 183－190. (in Chinese with English abstract)

[27] 王秀礼，卢永刚，袁寿其，等. 基于流固耦合的核主泵汽蚀动力特性研究[J]. 哈尔滨工程大学学报，2015，36(2)：213－217.

Wang Xiuli, Lu Yonggang, Yuan Shouqi, et al. Dynamic characteristics analysis of the reactor coolant pump variation based on fluid-structure coupling[J]. Journal of Harbin Engineering University, 2015, 36(2): 213－217. (in Chinese with English abstract)

[28] 裴吉，袁寿其. 离心泵非定常流动特性及流固耦合机理[M]. 北京：机械工业出版社，2014.

[29] 关醒凡. 现代泵理论与设计[M]. 北京：中国宇航出版社，2011.

[30] 袁寿其，施卫东，刘厚林，等. 泵理论与技术[M]. 北京：中国宇航出版社，2014.

Analysis of fluid-structure interaction characteristics for impeller of residual heat removal pump

Yin Tingyun, Pei Ji※, Yuan Shouqi, Wang Wenjie

(,,212013,)

In recent years, nuclear power has drawn increased attention because of its high efficiency and low pollution. Thus, a rising number of nuclear power stations are being constructed. The safety of nuclear station operations is mainly guaranteed by the residual heat removal system. Residual heat removal pumps (RHRP) are operated when the nuclear main pump stops working and the nuclear station needs to be maintained. The RHRP works in a complex environment, and its work status directly affects the performance of the entire plant. To ensure the reliability of the RHRP, the vibration characteristics of the rotor were analyzed using fluid-structure interaction theory. Stress and deformation analysis by partitioned solution for an impeller in a moving fluid was performed, and modal analysis of the impeller by monolithic solution was conducted in still fluid. For the partitioned method, there are two strategies for coupled solutions of dynamic fluid and structure interaction, one-way coupling and two-way coupling. Two-way coupling is typically used for large structural deformations. One-way coupling is suitable for the small structural deformation cases. In pump machinery, the impeller vibration caused by unsteady flow results in small deformations. Additionally, the feedback of the impeller motion onto the flow is small and therefore, can be neglected for most cases. Consequently, one-way coupling has been chosen, in which dynamic forces are transferred to the structure through the interface in a single direction at every time step. To understand the influence of the impeller shroud thickness on the resulting vibration characteristics, three impeller modifications were investigated and compared to the initial geometry under different flow rates. Moreover, five commonly used materials for an impeller were also evaluated. The three-dimensional turbulent flow was modeled utilizing a SST-turbulence model, and the numerical results were verified by the experimental data. The results showed that due to local structural differences between the pumps used in the numerical calculation model and the test measurement, as well as other effects, such as mesh quality, it was inevitable that there would be differences between the numerical calculation and the test measurements. However, the overall external characteristics of the numerical simulation were generally consistent with the performance of the test measurements, indicating that the flow-field calculation model can accurately predict its performance. By comparing with impellers adapted from four other materials and different shroud thicknesses, the vibration modes of the impellers were basically same for each order; however, the natural frequencies differed to some extent. The first order frequency of original impeller rotor was 394.17 Hz at hot condition and increased by 2.28% compared with cold condition, which was higher than blade passing frequency. Natural frequency of 1Cr13MoS was the highest among employed materials for each order mode, while ZG225-450 was the lowest. At design and off-design flow rates, the stress and displacement fields were similar. The displacement grew from the hub to the outer diameter, and each blade passage had a local maximum on the rear shroud. Moreover, the higher equivalent stress values can be observed in the junction between blade and shroud. Under three operating points, the peak values of stresses occurred in the middle of the junction between shroud and blade pressure side. Decreasing the head caused a significant reduction at the beginning of the blade passage. The stresses along defined paths were almost independent of the front shroud thickness, but peak values could be significantly reduced with a thicker rear shroud. The trendy of stress distribution between hot and cold condition was basically same. However, stress of hot condition was higher than cold one. Especially nearby the leading edge of the impeller, stress of hot condition increased more than 300%, compared to cold one. The results provide a theoretical basis for improving system performance and further study for more complicated dynamic analysis and fatigue analysis.

pumps; impellers; kinetics; residual heat removal pump; shroud thickness

10.11975/j.issn.1002-6819.2017.09.010

TH38

A

1002-6819(2017)-09-0076-08

2016-10-27

2017-04-05

国家科技支撑计划项目（2014BAB08B01）；国家自然科学基金青年基金（51409123）；江苏省自然科学基金青年基金（BK20140554）；江苏大学“青年骨干教师培养工程”

尹庭赟，男，博士生，主要从事水力机械空化空蚀机理研究。镇江 江苏大学国家水泵及系统工程技术研究中心，212013。 Email：tingyun_YIN@ujs.edu.cn

裴 吉，男，博士，副教授，主要从事离心泵非定常流动特性及流固耦合机理。镇江 江苏大学国家水泵及系统工程技术研究中心，212013。Email：jpei@ujs.edu.cn

尹庭赟，裴 吉，袁寿其，王文杰. 余热排出泵叶轮流固耦合特性分析[J]. 农业工程学报，2017，33(9)：76－83. doi：10.11975/j.issn.1002-6819.2017.09.010 http://www.tcsae.org

Yin Tingyun, Pei Ji, Yuan Shouqi, Wang Wenjie. Analysis of fluid-structure interaction characteristics for impeller of residual heat removal pump[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(9): 76－83. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2017.09.010 http://www.tcsae.org