蒸发循环泵泵轴性能的仿真分析.pdf

第48 卷第2 期2024年2 月doi:10.11832/j.issn.1000-4858.2024.02.015（1.四川大学机械工程学院，四川成都6 1 0 0 6 5；2.四川省自贡工业泵有限责任公司，四川自贡6 43 0 0 0）摘要：为了探究极限工况下泵轴受力变形对轴流泵性能的影响，采用ANSYS软件中流固耦合技术与CFD模块相结合的方法,模拟ZWX-1600型蒸发循环泵在不同工况下工作时，泵轴的受力情况,并将泵轴转速以及受到的扭矩、径向力、轴向力等拟合成关系曲线，基于以上数据使用nCode软件进行寿命预测,最终形成高可靠设计方案，为轴流泵的设计和制造提供参考。关键词：轴流泵；泵轴；流固耦合；流体仿真中图分类号：TH137;TH312文献标志码：B文章编号：1 0 0 0-48 58(2 0 2 4)0 2-0 1 1 6-0 8Simulation Analysis of Pump Shaft Performance ofEvaporation Circulating PumpMA Rong-ke,LIU Xiao-zhen,LI Bo-wei,FU Bo,LI Yu?(1.School of Mechanical Engineering,Sichuan University,Chengdu,Sichuan 610065;2.Zigong,Sichuan Industrial Pump Co.,Ltd.,Zigong,Sichuan 643000)Abstract:In order to investigate the influence of pump shaft deformation under extreme operating conditions on theperformance of axial flow pumps,this paper employs a method that combines the fluid-solid coupling technique andthe CFD module in ANSYS software to simulate the stress situation of the pump shaft of the ZWX-1600 evaporatingcirculation pump under different working conditions.The relationship curves of pump shaft speed,torque,radialforce,and axial force are fitted based on the data obtained.Using nCode software,the life prediction is performedbased on the above data,and a highly reliable design scheme is ultimately developed.This study provides valuablereferences for the design and manufacture of axial flow pumps.Key words:axial flow pump,pump shaft,fluid-structure interaction,fluid simulation引言强制循环泵（也称卧式轴流泵）是利用叶轮旋转所产生的沿着泵轴方向的水平推力进行工作。主要应用在真空制盐、隔膜法烧碱、熔盐、造纸、废水处理等工业活动中。为了进一步提高设备的生产能力,将强制循环泵换热器的传热系数提高，即为轴流式蒸发循环泵。蒸发循环泵作为强制循环泵系统的重要组成部分，有着流量大、耐高温、扬程低等特点。随着工况条件要求的不断提高，目前蒸发循环泵的设计研发逐渐向大口径,高可靠性等方向发展 。图1 为ZWX-1600型蒸发循环泵结构示意图。泵轴是轴流泵传递扭矩的关键零件之一，需要承受叶轮液压与气动Chinese Hydraulics&Pneumatics蒸发循环泵泵轴性能的仿真分析马荣柯，刘小振，李伯巍，傅波，李宇所产生的轴向力和径向力，同时在各种载荷的作用力下,其结构会产生一定量的变形 2。在泵轴与驱动电机的联接处,其直径为泵轴结构上最小,而且通常会设计键槽,这进一步影响了泵轴的刚度和稳定性。当泵轴发生弯曲变形时，可能会导致叶片与壁面碰撞，引发泵的故障和损坏,甚至可能对生产过程带来重大的安全隐患 3。因此,针对轴流泵泵轴的变形问题，必须收稿日期:2 0 2 3-0 6-1 9基金项目：四川省科技计划（2 0 2 3 YFG0226）;川大-自贡校地科技合作专项资金项目（2 0 2 0 CDZG-1）作者简介：马荣柯（1 9 9 8 一），男，四川自贡人,硕士研究生，主要研究方向为数字化设计与制造。Vol.48 No.2February.2024修回日期：2 0 2 3-0 7-2 72024年第2 期121.叶轮2.泵体3.机械密封4.主轴5.轴承箱6.轴承组件图1 ZWX-1600型蒸发循环泵结构示意图Fig.1 ZwX-1600 evaporation circulating pumpstructure diagram进行深人的分析和研究，以提高轴流泵的工作效率和安全性。轴流泵内部流场和叶轮泵轴之间存在着流固耦合(FSI)作用 4:大型轴流泵在工作时,水泵转轴旋转带动叶轮旋转，流体从吸水管进人高速旋转的叶轮，流经叶片后得到能量并往后输送。流体流经泵轴时,泵轴将受到各种力的联合作用，包括流体压力、离心力、重力、黏性力以及叶轮的反作用力等的作用,叶轮和泵轴的变形以及动力学响应反过来又会影响流场的分布，从而使流体载荷的分布和大小发生改变。因此，为了使最后的计算结果能给轴流泵的优化设计提供更可靠的数据参考，可将泵轴叶轮与流场进行双向流固耦合联合求解。目前，流固耦合技术已广泛运用在流体机械静力学以及动力学仿真研究中,并取得了一定成果。杨丹丹等 5以轴流泵为例,采用沉浸边界与流固耦合及大涡模拟相结合的方法，对泵内流场及鱼体通过泵流道的运动过程进行了数值模拟;杨庆俊等 6 1 针对泵源与机体共同激励下的管路振动，建立了流固耦合的总体动力学方程，发现液压系统在泵源谐波激励下的振动响应表现为宽频域的强迫振动;许正萱等 7 运用单向流固耦合技术，求解试验转速和泵压条件下螺杆泵的漏失量，进一步得到其容积效率；冀宏等 8 运用流固耦合技术分析了液压电机泵在工作时全域温度场分布特征；显文雄等 9 基于流固耦合的算法理论，研究了不同工况下预应力对轴流泵叶轮模态的影响。本研究采用流固耦合技术与CFD模块相结合,基于ANSYS计算了ZWX-1600型蒸发循环泵在不同工况下泵轴的受力情况,通过数据分析,得到泵轴转速与所受的扭矩、径向力、轴向力之间的关系曲线,并使用液压与气动nCode 软件进行了寿命预测,为泵轴的设计提供参考。1物理模型与数学模型1.1基本参数ZWX系列轴流式蒸发循环泵一般应用在大流量、低扬程场合 1 0)。ZWX-1600型蒸发循环泵基本参数及工况条件如表1 所示。基于SolidWorks三维建模软件3456117对叶轮及流体域部分进行三维建模，其主要计算域如图2 所示。表1 基本参数及工况条件Tab.1 Basic parameters and conditions of work conditions泵型号介质名称设计流量/m3 h-1扬程/m黏度/MPas工作温度/泵进出口直径/mm泵转速/r min-1泵进出口流速/ms=1泵轴材质Fig.2Calculation domain model1.2网格模型三维模型的网格划分对仿真结果有决定性影响，网格质量会直接作用于仿真计算的敛散性与精确性。基于SolidWorks软件对叶轮及泵轴部分进行三维建模，并将建立的模型导人ANSYSMesh模块中进行网格划分。采用六面体网格，对叶轮、泵轴、进水管和出水管分别进行离散。为了模拟流体介质在进出叶轮流域时的真实情况，对叶轮前端的流域及其管道进行了加长处理，使流体介质能均匀地流入叶轮流域。流体介质通过叶轮后，经过9 0 的圆弧管道流出。叶轮在工作时,其叶片的正面受力面和背面的厚度较薄，ZWXIIENZ-1600-A1-220盐浆(氯化钠）30000320002.832noml13(max)13016002204.14 4.422205整轴图2 计算域模型118网格划分具有一定难度。因此在对其网格划分时,先对叶片部分进行加密处理，从而保证网格质量满足要求。确定网格数量时，在保证结果精度的情况下,还需要考虑计算机性能,因此需要进行网格无关性验证。本研究选取了5组网格，在设计工况下，进行仿真求解各组的扬程,得到的结果如图3 a所示。从图中可以看出,随着网格数量的增加,扬程变化幅度较小,趋向于定值。从50 8 万到8 8 9 万网格，结果数据相差约3.4%；从8 8 9 万到1 2 53 万网格，结果数据相差约1.6%；故可认为该仿真结果已经收敛，网格无关性验证完毕。综合考虑计算机性能与仿真精度要求，最终确定泵轴网格和计算区域网格共8 8 9 万个，如图3 a、图3 b所示。其中泵轴部分1 8 8 万个，流体域部分7 0 1 万个。3.02.8上2.6上u/$2.42.22.01.801.3数学模型轴流泵工作时,叶片表面所受到的力将直接作用到泵轴的键槽上，而叶轮作为整个泵装置的一部分,其液压与气动受力状态与进出水流道以及其他部分流动状态息息相关。因此要准确求解泵轴在工作时的受力状态,必须对泵装置进行流体仿真。在计算模型中，整个轴流泵装置内的流域被划分为叶轮附近的旋转流域和除此之外的流域。选用参考系来模拟其内部流动,对于非定常不可压缩流体,在以恒定角速度绕X轴旋转的相对直角坐标系下，连续性方程为：3U:=0动量方程：+0 x;1卫atax;x式中，T,为雷诺应力。选用以雷诺平均方程为基础的RNG k-8湍流模型：at10520406080网格数量a）网格无关性验证a)Mesh independence validationb)泵轴网格b)Pump shaft meshc)计算域网格c)Computational domain mesh图3 网格模型Fig.3Mesh model第48 卷第2 期(1)au;(2)akG+G一p8atax;Ci。%(Ce +Ca.C.)-Ca p100120P81401式中，G由层流速度梯度而产生的流动能G一一由浮力产生的湍流动能Cle,C2e,Cse,u常系数k。k 方程和方程的流Prandtl数,其取值分别为 C1。=1.44,C2e=1.92,C3e=1.52,C=0.09,0=1.0,g。=1.3一紊流黏性系数上述各方程构成了封闭的非线性偏微分方程组，将上述方程组变换至贴体坐标系中,通过ANSYS软件中CFD分析与SIMPLEC算法结合计算，即可求出目标收敛解 1 2 2数值模拟2.1边界条件在CFD模块中，载荷主要由外部作用力和边界条件构成。在实际仿真过程中，载荷可以直接施加在实(3)(4)2024年第2 期体模型上或者有限单元模型上。本研究根据实际情况出发，将载荷施加在泵轴有限单元模型上。边界条件的设置会直接影响求解的准确性 1 3,主要边界条件如表2 所示，表中步长时间根据正常工况下叶轮转过1 所需要的时间确定，经过计算设置为0.0 1 0 1 8 5s。由于叶轮和泵轴都处于转动状态，因此需要运用MRF技术来模拟叶轮和泵轴的转动,对转动域开启MRF功能，设置转轴为X轴。添加新材料盐液，密度1450kg/m,黏度2 MPas,将泵轴转速和人口处流速设置为变量。为保证仿真方案的可行性和正确性，在进行仿真之前，还需要对该模型进行残差曲线计算。一般情况下，对于RNGk-s流模型，残差曲线默认设置为1 0-3。当收敛后的残差值达到设定值，即可认为计算结果收敛。本研究计算完的残差曲线如图4所示，当迭代1 0 0 次左右残差曲线开始逐渐平缓，最终计算结果达到了预先设置的数量级,故计算收敛。表2 边界条件设置Tab.2Boundary condition settings流体流模型进口面出口面壁面动静耦合面步长/s收敛精度RenidualcontioutyXieloctyyhniocryzweloctyKepsilon液压与气动固体形变后对流场的影响，两门力学交叉，便有了流固耦合力学这一门力学分支 1 4。流体与固体之间的相互作用是流固耦合力学的重要特征。在流体域中，固体受流体载荷的影响，表面会发生微小形变和位移,这些形变和位移又会反馈给流场，从而改变流体的运动轨迹和载荷大小。因此,在流体力学研究中,需要考虑到流场和固体场之间的相互作用，及动态变化。采用流固耦合技术，可以有效地分析复杂的流体域问题,并为工程实践提供有力支持。此外，流固耦合力学还具有广泛的应用前景，在航空、汽车、轮机以及建筑结构等领域都有大量应用。利用ANSYS软件中CFD模块和Workbench功能互相耦合的方法，对流体域和结构域进行了有限元分析。在计算过程中,设置相同的计算时间和时间步长，首先计算流体域，并将计算后得到的结果（主要为压强）传递给固体域进行分析计算。固体域计算后得到的结果(包括变形、速度、位移和加速度等)再反馈给流体域进行分析。这样不断进行反复计算和迭代，直到整个耦合场的求解结束。非定常，不可压缩图5是泵轴流固耦合分析的工作流程图，首先计RNG k-8算出Fluent模块里轴流泵内流场的仿真结果,接着将给定初始速度流体域内各个节点的速度、应力等物理量导人Static自由出流Structure模块，进行耦合求解，最终得到收敛解，就能无滑移求解出泵轴等的应力和应变。最后将泵轴的极限受力情况导入ANSYS nCode模块进行疲劳分析,预测使用Stage寿命。0.0101851 10-31e+0210+0110+00119中P10-011e-0216-0316-04图4残差曲线Fig.4Residual curve2.2流固耦合分析在流体力学研究中，流体与固体之间总存在着相互作用。在研究流场对变形固体的作用的同时，考虑50100150200250300350400450Iterations图5流固耦合计算流程Fig.5 Fluid-structure interaction calculation process3仿真分析结果3.1流场仿真结果设定不同的流体域人口速度以及不同的转速，以此模拟不同工况。图6 是在人口速度为4.42 m/s，转120PresureContour32.084e+05.054e+824e+2.605e+.303e+040.0000+00IPa)a)Pressure distribution contour map on back of impellerPressurgContour2.000e+00Pab)Pressure distribution contour map on front of impeller8920+04D040390+043270-040477203.2160+04+04(Pa液压与气动速为2 2 0 r/min情况下的仿真结果。图6 a为叶轮反面所受到的应力云图，图6 b为叶轮正面所受到的应力云图。可以看出叶片正反两面压差很大，压力分布沿着叶片从进口处到出口处呈逐渐递增趋势，说明叶片对流体作用明显。图6 c为流体域中泵轴所受到的应力分布图,由图可知,流体域对泵轴的最大应力发生在泵轴与管壁交接处，流体对泵轴底部有较大冲击压力。图6 d为流场流线图。从流线可以看出，仿真结果中流体流动较为规律，且在设计工况下，出口处流体速度大a）叶轮背面受力云图致为4m/s，与设计出口速度接近，表明该流场模拟仿真结果较为准确。3.2泵轴受力分析结果泵轴主要受到电机的扭矩、叶轮的重力、空气对叶轮和泵轴的作用力,受到的约束主要为两处轴承支撑，另外双列圆锥滚子轴承能够限制轴双向轴向位移,因此在此轴肩位置处对轴向位移也有约束。此次分析主要关注旋转区域与叶轮和泵轴接触位置的应力情况，通过对叶轮与泵轴接触部位的应力分析,能够将作用m于叶轮的应力等效到泵轴轴端处,将作用于泵轴的力b）叶轮正面受力云图直接施加到泵轴上。另外泵轴在转动过程中有恒定的转速，稳态中将恒定的转速简化为惯性力，不会出现转动，但是转速会产生离心力，因此还需要给整个泵轴一个转动速度，速度设置为2 2 0 r/min对应流体分析时的速度。载荷及约束如图7 所示。国CyindricalSupport-0.mmJMoment2.0016e+007NmmMomem2:2.0836+007NmmDDisplecementUStandard EarthGrailyi9B06mmyssFoed SuppontRotationalVelecityeImgortedPressure:2第48 卷第2 期c）流体域内泵轴受力云图c)Force distribution contour map of pump shaft within fluid domain1.8180+011.3620+019.000+00图7 泵轴载荷及约束4.540e+00Fig.7Pump shaft load and constraints0.000e+00通过 CFD-POST进行后处理可以得到液体对叶轮的XY两个方向的作用力以及液体对轴的XY两个方向作用力。再通过直接优化模块，将转速设置为变量，范围为0 2 2 0 r/min,将对叶轮的力以及对泵轴的力d)流体矢量图分别设置为目标,通过系统自带的 Screening算法自动d)Fluid vector diagram生成8 组样本点。将样本点分别做成相对应的曲线图6 Fluent仿真结果图,得到在流体域入口速度为4.42 m/s和0 m/s条件Fig.6Fluent simulation result下，泵轴转速与泵轴受力情况的关系，如图8 所示。2024年第2 期656055504540352000001000000Z-100000-200000-300000-400000L8000060000U.N/400001200000-200004000350030002500Z20015001000500液压与气动分析图8 可知，泵轴的最大受力情况出现在流体-0 m/s-.4.42 m/s0100n/r-minla)最大应力a)Maximum stress0m/s-4.42 m/s0100n/r:min-1b）最大轴向力bMaximum axial force0 m/s-4.42 m/s0100n/r-min-l)最大扭矩c)Maximum torque121域人口处流体流速为0,泵轴转速为2 2 0 r/min时。该情况下得到泵轴的受力及变形情况如图9 所示。62.732Max55.76248.79341.82334.85427.88420.91513.9456.97572000.0061772Min6.6193Max5.88385.148344129Lutomat2.84192.200414710.735-48oMinMax200Fig.9Force and deformation of pump shaft一般情况下，泵轴最容易被破坏或疲劳失效的位置位于轴颈过渡圆角处，该处用于安装联轴器，需要承受较大的力矩和弯矩，并且该处直径尺寸最小 1 5。图9 中泵轴最大应力点与最大应变点与实际情况相符，最大等效应力为6 2.7 3 2 MPa,最大变形数值为6.6 1 9 3 mm。将应力分布导人ANSYSnCode模块进行疲劳分析,寿命预测结果如图1 0 所示，泵轴最危险处受力循环次数大于3 1 0 次，满足实际生产需要。2004结论本研究采用ANSYS软件中流固耦合技术与CFD0 m/s4.42 m/sMAya）最大受力处云图a)Cloud map of maximum stress pointb）最大变形处云图b)Cloud map of maximum radial force图9 泵轴受力及变形情况BeredoutedfL.1400141.300-0131-4828-0121-004-0111.0270-0102.1976-0092.805e-002.856+-0072.256-0062.7120-005AENode145632700d)Maximum radial force图8 泵轴受力情况曲线图Fig.8Curve chart of pump shaft stress100n/r-min-1d)最大径向力200图1 0 nCode 寿命预测结果Fig.10nCode life prediction results122模块相结合的方法，对ZWX-1600型轴流泵内部流场进行数值模拟，分析了轴流泵流场内泵轴的受力情况，得到如下结论：（1）针对轴流泵主要流体域部分，对其进行了流体仿真分析。在给定的设计工况下，泵轴压力云图和流场线速度分布较为合理，仿真结果与预期性能基本符合,表明仿真结果具有较高的准确性；（2）运用CFD模块与流固耦合技术,对各个工况下泵轴受力情况进行模拟，将求解出的多组样本点拟合成曲线，更加直观的反映出泵轴受力与转速等的关系。在极限工况下,泵轴与电机链接键槽处出现应力集中,最大应力不到6 3 MPa,远小于泵轴材料许用应力 1 6。泵轴最大变形出现在轴端处,与泵轴整体5m左右长度相比,变形幅度6.6 mm处在合理范围内。仿真结果说明泵轴结构设计合理，强度可靠；（3）将泵轴极限工况下的受力情况导入nCode模块进行疲劳仿真，结果表明,泵轴最危险处受力循环次数大于3 1 0 次，满足工厂的实际生产活动中，对泵轴的寿命需求,为轴流泵的设计和制造提供了有价值的参考。参考文献：1梁武科,侯聪,董玮,等.卧式轴流泵空化特性的数值分析J.排灌机械工程学报,2 0 2 0,3 8（8）：7 6 4-7 6 9.LIANG Wuke,HOU Cong,DONG Wei,et al.NumericalAnalysis of Cavitation Characteristics of Horizontal Axial-flowPump J.Journal of Drainage and Irrigation MachineryEngineering,2020,38(8):764-769.2邢晋,张建军.闭式循环水泵泵轴断裂分析及处理 J.电力学报,2 0 2 1,3 6(3);2 54-2 6 0.XING Jin,ZHANG Jianjun.Analysis and Treatment ofPump Shaft Fracture of Closed Circulating Water Pump J.Journal of Electric Power,2021,36(3):254-260.3胡兆祥,彭思桐，李清泉，等.某核电厂常规岛分离水泵中段螺栓和泵轴断裂失效分析及对策研究 J.水泵技术，2021(4):32 37.HU Zhaoxiang,PENG Sitong,LI Qingquan,et al.FractureFailure Analysis and Countermeasures of the Middle SectionBolt and Pump Shaft in the Separation Water Pump forConventional Island of a Nuclear Power Plant J.PumpTechnology,2021(4):32-37.4刘松，吴先梅,彭修乾,等.基于双向流固耦合的输流圆管液压与气动应力应变响应分析 J.振动与冲击，2 0 2 1,40（2 0）：7 3-79.LIU Song,WU Xianmei,PENG Xiuqian,et al.Analysis ofStress-strain Responses of a Liquid-filled Pipe Based onTwo-way Fluid-structure Interaction J.Jo u n a l o fVibration and Shock,2021,40(20):73-79.5 林杨丹丹,赵明翔,叶维祥，等.基于沉浸边界法的轴流泵内鱼体运动数值模拟 J.水力发电学报，2 0 2 1,40（7）：9 5-104.YANG Dandan,ZHAO Mingxiang,YE Weixiang,et al.Numerical Simulations of Fish Movement in Axial PumpUsing Immersed Boundary Method J.Jo u r n a l o fHydroelectric Engineering,2021,40(7):95-104.6杨庆俊，董日治，罗小梅，等.泵源与机体共同激励下液压管路振动特性 J.液压与气动,2 0 2 1,45（6）：1 6 3-1 7 0.YANG Qingjun,DONG Rizhi,LUO Xiaomei,et al.Vibration Characteristics of Hydraulic Pipe Network UnderMulti-source Excitation J.C h in e s e H y d r a u lic s&Pneumatics,2021,45(6):163-170.7许正萱，张帆，陈，等.流固耦合作用下斜流泵转子动力学特性研究 J.农业机械学报，2 0 2 2，53（S02）：1 7 9-187.XUZhengxuan,ZHANG Fan,CHEN Ke,et al.RotorDynamic Characteristics of Oblique Flow Pump under Actionof Fluid-structure Interaction J.T r a n s a c t i o n s o f t h eChinese Society for Agricultural Machinery,2022,53(S02):179-187.8 冀宏,邢晖晖，孙飞，等.基于流固耦合传热分析的液压电机泵温度场特征 J.液压与气动,2 0 2 3,47（3）：1-8.JI Hong,XING Huihui,SUN Fei,et al.Temperature FieldCharacteristics of Hydraulic Motor Pump Based on Fluid-structure Coupling Heat Transfer Analysis J.ChineseHydraulics&Pneumatics,2023,47(3):1-8.9晃文雄，王均儒，何彬，等.半开式复合叶轮离心泵不同叶顶间隙的水力特性分析 J.液压与气动,2 0 1 8，（1 1）：8 9-93.CHAO Wenxiong,WANG Junru,HE Bin,et al.HydraulicCharacteristics Analysis of Centrifugal Pump with UnshroudedComposite Impeller at Different Tip Clearances J.ChineseHydraulics&Pneumatics,2018,(11):89-93.10曹永，肖国华，李宇.ZWX500型蒸发循环泵进水短管空蚀空化分析 J.机械制造,2 0 1 6,54（3）：51-52.第48 卷第2 期2024 年第2 期CAO Yong,XIAO Guohua,LI Yu.Analysis of CavitationCavitation of Influent Short Pipe of ZWX500 EvaporativeCirculation Pump J.Machinery,2016,54(3):51-52.11何金辉，李明广，陈锦剑，等.考虑动态流体网格的颗粒-流体耦合算法J.上海交通大学学报，2 0 2 1,55（6)：645-651.HE Jinhui,LI Mingguang,CHEN Jinjian,et al.Particle-Fluid Coupling Algorithm Considering Dynamic Fluid MeshJ.Journal of Shanghai Jiaotong University,2021,55(6):645-651.12 宋文武,符杰,李庆刚,等.ZWX系列轴流泵水力模型的研究 J.中国农村水利水电,2 0 0 7，（1 0）：9 9-1 0 1,1 0 4.SONG Wenwu,FU Jie,LI Qinggang,et al.Research onHydraulic Model of the ZWX Series Axial Pump J.China Rural Water and Hydropower,2007,(10):99-101,104.13黄龙龙,王艳华,刘勇,等,基于 CFD 和FSI的摆线泵仿真分析 J.液压与气动,2 0 2 0，（7）：6 9-7 5.HUANG Longlong,WANG Yanhua,LIU Yong,et al.Simulation Analysis of Cycloid Pump Based on CFD andFSI J.Chinese Hydraulics&Pneumatics,2020,(7):引用本文：马荣柯，刘小振，李伯巍，等.蒸发循环泵泵轴性能的仿真分析 J.液压与气动,2 0 2 4,48(2）：1 1 6-1 2 3.MA Rongke,LIU Xiaozhen,LI Bowei,et al.Simulation Analysis of Pump Shaft Performance of Evaporation Circulating Pump J.Chinese Hydraulics&Pneumatics,2024,48(2):116-123.液压与气动69-75.14 施卫东,王传龙,谭林伟,等.基于数值模拟的单叶片离心泵能量损失分析J.排灌机械工程学报，2 0 2 3,41(4):325-330.SHI Weidong,WANG Chuanlong,TAN Linwei,et al.Energy Loss Analysis of Single-blade Centrifugal PumpBased on Numerical Simulation J.Journal of Drainageand Irrigation Machinery Engineering,2023,41(4):325-330.1 5王维，尤保健，陆荣，等.基于有限元法的多级离心泵结构强度及流固耦合分析 J.水泵技术,2 0 2 0，（3）：2 9-32.WANG Wei,YOU Baojian,LU Rong,et al.StructuralStrength and Fluid-structure Interaction Analysis of MultistageCentrifugal Pump Based on Finite Element Method J.Pump Technology,2020,(3):29-32.16刘立伟，姚益轩，杨建炜，等.铀矿用不锈钢泵轴断裂原因分析 J.表面技术,2 0 1 9,48（1 2）：2 8 9-2 9 5.LIU Liwei,YAO Yixuan,YANG Jianwei,et al.Causes ofFractures of Stainless Steel Pump Shaft Used in UraniumMine J.Surface Technology,2019,48(12):289-295.123