常规岛汽机旁路调节阀流激共振研究.pdf

流 体 机 械第 51 卷第 9 期2023 年 9 月 27 收稿日期：2022-06-02 修稿日期：2023-04-13基金项目：深圳市科技计划资助项目（JSGG20210629144536003）doi：10.3969/j.issn.1005-0329.2023.09.005常规岛汽机旁路调节阀流激共振研究颜炳良1，金 辉1，黄 冲2，沈恒云3，胡迎港3，李树勋3（1.重庆川仪调节阀有限公司，重庆 400707；2.中广核工程有限公司，广东深圳 518124；3.兰州理工大学 石油化工学院，兰州 730050）摘 要：针对核电常规岛汽机旁路阀小开度下易发生振动的问题，设计了一种汽机旁路阀，基于流激共振理论及流固耦合方法，对其阀芯组件进行流激振动分析。选取 30%开度工况，采用CFD 软件进行瞬态流场仿真分析，以获得阀芯组件流固耦合处的流体压力脉动时域信息，将时域信息进行快速傅里叶变换（FFT）得到流体作用于阀芯组件的激振频率。采用ANSYSWorkbench 软件对典型开度进行流固耦合模态分析，得到阀芯组件的模态频率。通过流体激振频率与模态频率对比进行进一步分析研究。结果表明：典型工况下该型调节阀流体激振的频率区为 8481272Hz，阀芯组件结构 1 阶固有模态频率为 3170.1Hz，随着阶数增大依次递增，6 阶模态频率与流体激振频率均不接近，因此该型汽机旁路阀在典型工况下不会发生流激共振，验证了该设计的有效性，可为同类常规岛汽机旁路调节阀防共振结构设计提供参考。关键词：旁路调节阀；流固耦合；流激共振；数值模拟中图分类号：TH134 文献标志码：A Research on flow-induced resonance of conventional island steam turbine bypass control valveYANBingliang1，JINHui1，HUANGChong2，SHENHengyun3，HUYinggang3，LIShuxun3（1.ChongqingChuanyiControlValveCo.，Ltd.，Chongqing 400707，China；2.ChinaGeneralNuclearPowerCorporation，Shenzhen 518124，China；3.SchoolofPetrochemicalTechnology，LanzhouUniversityofTechnology，Lanzhou 730050，China）Abstract：Fortheproblemthatthebypassvalveofthesteamturbinebypassvalveinthenuclearpowerconventionalislandispronetovibrationatsmallopenings，asteamturbinebypassvalvewasdesigned，andthentheflow-inducedvibrationanalysisofthevalvecoreassemblywascarriedoutbasedontheflow-inducedresonancetheoryandthefluid-solidcouplingmethod.Theworkingconditionof30%openingwasselected，andCFDsoftwarewasusetoconducttransientflowfieldsimulationanalysistoobtainthetime-domaininformationoffluidpressurefluctuationatthefluid-solidcouplingofthevalvecoreassembly，andperformfastFouriertransform（FFT）onthetime-domaininformationtoobtaintheexcitationfrequencyatwhichthefluidactsonthespoolassembly.ANSYSWorkbenchsoftwarewasusedtoconductfluid-solidcouplingmodalanalysisontypicalopeningstoobtainthemodalfrequencyofthevalvecoreassembly.Furtheranalysisandresearchwerecarriedoutbycomparingthefluidexcitationfrequencywiththemodalfrequency.Theresultsshowthatundertypicalworkingconditions，thefluidexcitationfrequencyrangeofthistypeofcontrolvalveis8481272Hz，thefirst-ordernaturalmodefrequencyofthevalvecoreassemblystructureis3170.1Hz，andasitsuccessivelyincreaseswiththeincreaseofthenumberoforders，andthesixth-ordermodalfrequenciesarenotclosetotheexcitationfrequencyofthefluid，sothistypeofturbinebypassvalvewillnotexperienceflow-inducedresonanceundertypicaloperatingconditions，whichverifiesthefeasibilityofthisdesignmethodandcanbeusedfortheanti-resonancestructuredesignofconventionalislandturbinebypasscontrolvalvesforreference.Key words：bypassregulatingvalve；fluid-structurecoupling；flowinducedresonance；numericalsimulation0 引言汽机旁路阀主要作用是保护常规岛高温加热器，同时也能作为汽轮机的保护装置。若发生事故，也可作为一回路的冷源，以保证反应堆安全、经济运行1。汽机旁路调节阀是一种典型的高压差调节阀，由于其结构的特殊性，小开度下易发生振动，发生振动会对阀门内部结构造成破坏，影28FLUID MACHINERYVol.51，No.9，2023响正常使用，严重时可能导致阀杆断裂2。因此，对汽机旁路阀进行流激共振分析十分必要。国内外学者对调节阀进行了结构优化和仿真研究。朱禄3针对先导式调节阀启闭及调节过程中阀芯组件振动等问题，采用压力脉动频域分析法对先导式调节阀进行流激共振模态分析。王伟波等4针对核电常规岛热力系统关键阀门振动问题，采用FLUENT 软件对各工况模型进行瞬态流场仿真分析，并对调节阀各工况进行流激共振分析。TECZA 等5对汽轮机节流阀进行FEA与CFD 流固耦合分析，通过结构改进使其正常运行。ENGEDA6对汽轮机蒸汽控制阀进行实验研究，结果表明采用对称结构可以提高阀门的稳定性。刘丽等7通过对气体管道阀门小开度工况下的振动管路进行模态分析发现，管道振动是由于阀后压力脉动引起。以上学者虽然对不同类型的调节阀进行了流激共振分析研究，但研究对象节流件结构相对单一，而本文设计的汽机旁路阀采用节流套筒和迷宫盘片组合的方式进行节流，并在阀后设有消声孔板，结构复杂，对采用该设计的汽机旁路阀进行研究十分必要，可为调节阀防共振结构设计提供参考。因此采用流固耦合分析方法，对DN125 汽机旁路阀进行流激共振数值模拟研究，首先对其实际工况下稳态流场进行分析得到内部介质对阀芯组件表面的作用力，再进行瞬态模拟得到其固有频率，最后进行流固耦合模态分析判断典型工况下该设计结构是否发生共振。1 内部流场分析汽机旁路阀在小开度下极易发生振动，30%开度为典型工况中开度最小位置，因此选取汽机旁路阀 30%开度下的流场进行分析。进行流固耦合分析需要先进行内部稳态流场模拟和瞬态流场模拟，以获得固体域所受压力场及流体激振频率。1.1 结构设计及流道模型建立DN125 汽机旁路阀设计参数见表 1。汽机旁路阀通过阀杆带动套筒上下移动，以改变套筒节流小孔及迷宫盘片的流通面积来实现对流体的控制，阀杆采用先导结构，通过先导阀芯的开启平衡介质力，减小驱动力矩。DN125 汽机旁路阀的三维模型，如图 1 所示。表 1 DN125 汽机旁路阀设计参数Tab.1 DesignparametersofDN125steamturbinebypassvalve名称数值公称通径/mm125公称压力/MPa15设计温度/-29350注：介质为水蒸气。图 1 DN125 汽机旁路阀三维模型Fig.1 3DmodelofDN125steamturbinebypassvalve为防止使用过程中可能会出现的流激共振问题，该汽机旁路阀设计时采用多孔式与盘片迷宫式组合套筒进行节流，并在套筒及阀后增加多孔消声组件及消声孔板。多孔式套筒在流量调节范围内采用特殊的多孔斜排列设计，1%开度变化引起的流量系数Cv变化波动较小，如图 2 所示。小开度下结合多层迷宫槽设计，满足了 1.5%行程开度变化 3kg/s 的设计要求。图 2 多孔式套筒示意Fig.2 Schematicdiagramofporoussleeve采用ANSYSDM 软件对三维模型进行反向建模，为确保内部流动为湍流，在阀门前后分别建立长度为 2 倍和 6 倍阀门口径的直管段。汽机旁路阀 30%开度流道模型如图 3 所示。图 3 汽机旁路阀 30%开度流道模型Fig.3 Flowpathmodelofsteamturbinebypassvalvewith30%opening29颜炳良，等：常规岛汽机旁路调节阀流激共振研究1.2 网格划分与网格无关性检验流道模型网格由ANSYSMesh软件划分生成，由于汽机旁路阀内腔形状和流动状态复杂，采用自动网格划分方法生成阀体流域网格，并在套筒处、节流孔进行局部加密，两端直管段采用六面体网格划分，以节约计算资源。选取流量值为目标进行网格无关性检验，划分流体网格单元数从 4296782 到 5451959，在压差为 100kPa 的边界条件下，计算 3 种不同网格的流体流量值见表 2。表 2 流道网格无关性检验Tab.2 Runnergridindependencetest网格类型节点数单元数计算流量/（kg s-1）网格 188229642967823.55网格 296953047231963.63网格 3111949854519593.65根据表 2 网格无关性检验结果，选取网格数为 472 万的网格模型进行模拟计算。DN150 汽机旁路阀 30%开度网格模型如图 4 所示，套筒、阀芯处网格结构剖视如图 5 所示。图 4 流道模型网格结构Fig.4 Gridstructurediagramofrunnermodel图 5 套筒、阀芯处网格结构Fig.5 Cross-sectionalviewofthegridstructureatthesleeveandvalvecore1.3 流体力学控制方程计算汽机旁路阀稳态流场时采用湍流模型为标准k-方程，湍流动能k 与湍流动能耗散率方程分别为式（1）及式（2）：()+()=+ktxU kjj ()+()=+xkxPPjjjtkjk+k kb（1）()+()=+txU kjj()+()=+xxkC Pjtkjk+1+()CC Pb21（2）式中，Pkb和Pb代表浮力的影响；C1，C2，k，的值分别为 1.44，1.92，1.0 和 1.3；Pk为黏性力和浮力的产生项。1.4 稳态流场模拟根据 30%开度实际工况，设定边界条件为进口压力为 6.86MPa，出口压力为 0.15MPa，温度为284.4，介质为水蒸气。实际工况下饱和水蒸气物理参数性能见表 3。表 3 实际工况下饱和水蒸汽的物理参数性能Tab.3 Physicalparametersandpropertiesofsaturatedsteamunderactualworkingconditions介质名称水摩尔质量/（g mol-1）18温度T/284.4压力/MPa6.86密度/（kg m-3）35.7蒸汽黏度/（Pa S）18.9考虑介质可压缩性，采用实际气体状态方程，计算流域内介质在不同位置的压力、密度及黏度等其他热力学参数。选择稳态压力基分离求解器（SIMPLE 算法）作为模型求解器，采用标准k-双方程湍流模型，离散格式采用highresolution 格式。30%开度流域内压力、流速、流线云图如图 6所示。由图 6（a）可知，在 30%开度时，阀门入口及套筒部位压力相对较高，低压区主要集中在阀芯下部的流道区域，出口段压力分布相对均匀；由图 6（b）可知，阀门入口速度分布相对稳定，在速度为 66.9m/s 以内时，套筒底部速度分布不均匀，中心区域流速较低，套筒两侧形成高速区，流道出30FLUID MACHINERYVol.51，No.9，2023口阀体部分区域流速较高，结合图 6（c）流线可知，阀芯底部及流道出口区域流线分布较混乱，有涡流产生。（a）压力云图（b）流速云图（c）流线图 6 汽车旁路阀 30%开度流场云图Fig.6 Clouddiagramoftheflowfieldofthebypassvalveoftheautomobileat30%opening文献8基于计算流体力学（CFD）与流固耦合模态分析方法，对DN250 多级降压调节阀的涡激振动特性进行了分析，并用实际试验验证了仿真模型。本文的旁路调节阀结构与该文献相似，采用的流场计算方法及后续的研究方法与该文献原理一致，因此，可认为本文的仿真模型具有一定的可信度。1.5 瞬态流场模拟对 30%开度下的内部流场进行瞬态流场模拟，以监测阀芯组件的受力及激振频率初始数据。将稳态模拟结果作为瞬态模拟的初始值，以时均N-S 方程为流体流动基本控制方程，湍流模型采用标准k-湍流模型。采用具有二阶精度的迎风格式，以基于离散单元中心有限体积法流域离散的方法，通过Couple 耦合求解方法进行流场的速度压力求解。边界条件设置为压力进口和压力出口，壁面采用无滑移壁面。迭代收敛控制方程采用RMS 方法，残差判定收敛边界设置为 110-4，参考压力为大气压 0.1MPa，设置时间步长为0.0001s，计算总时长为 0.1s。采用CEL 语言监测汽机旁路阀阀芯受力随时间变化规律，每一个时间步保存一个数据，对数据进一步处理，可得阀芯面受力时域曲线如图 7所示。图 7 阀芯面受力时域曲线Fig.7 Forcetimedomaincurveofvalvecoresurface由图 7 可知，在时间为 00.02s 内，阀内介质逐渐充满，汽机旁路阀阀芯力呈上升趋势，在时间为 0.02s 后，阀芯受力在一定范围内波动，并随着时间变化呈现出随机性。为获得流体湍流脉动引起的激振频率，通过快速傅里叶变换（FFT），计算阀芯受力频域幅值分布如图 8 所示。图 8 阀芯面受力频域曲线Fig.8 Forcefrequencydomaincurveofvalvecoresurface由图 8 阀芯受力频域幅值可知，在 0200Hz内，由于流体未完全通过节流件，阀芯受力不稳定，振幅大。阀芯受力幅值在 1200Hz 内有峰值点，其中最大频率在 1060Hz 左右。其他频率下压力脉动幅值较小。依据结构共振原理，工程中常把 0.81.2 倍的激振频率范围称为共振区9，由此得到流体激振的共振频率区为 8481272Hz。其中 800Hz 在共振区间附近，若结构的固有频率也为 800Hz，则其振幅为 0.3，阀门出现大幅振动，即认为发生流激共振。312 流固耦合前处理2.1 阀芯组件材料性能根据文献 10查询阀芯组件的材料参数，见表 4。表 4 先导式汽机旁路阀零部件材料参数Tab.4 Materialparametersofpilot-operatedsteamturbinebypassvalveparts零件材料密度/（kg m-3）弹性模量E/（GPa）泊松比先导阀芯316L80271950.31阀杆66079202010.31阀芯/套筒F6a80272010.312.2 流固耦合简化模型建立汽机旁路阀在实际工作过程中的阀芯面轴向及横向受力具有一定的相关性，主阀芯和先导阀芯之间通过弹簧连接。在ANSYSWorkbenchDM中，将实体域与流道模型全部保留，如图 9 所示。图 9 汽机旁路阀流固耦合实体模型Fig.9 Fluid-structureinteractionentitymodelofsteamturbinebypassvalve为简化计算，使用质量集中法，将弹簧质量叠加在先导阀芯上，简化为先导阀芯、阀杆、活塞及套筒组成的阀芯组件，阀芯组件多自由度力学模型如图 10（a）所示，实体模型如图 10（b）所示。（a）简化力学模型 （b）阀芯组件剖视图 10 阀芯组件简化示意Fig.10 Simplifiedschematicdiagramofvalvecoreassembly2.3 固体域边界载荷施加流体对固体域的影响由稳态结果导入，再按照实际的配合情况对阀芯组件固体域施加边界载荷。根据阀芯组件实际工况下的工作位置，对简化模型施加静力场约束，如图 11（a）所示，简化模型固体域与稳态流场压力场的耦合如图11（b）（c）所示。（a）施加约束 （b）压力耦合 （c）耦合剖面图 11 固体域载荷示意Fig.11 Schematicdiagramofsoliddomainload3 流固耦合模态分析模态分析是将流场信息导入静力场的一种有预应力的分析方法11-14。3.1 流激共振分析理论基于流激共振产生的原理，通过改变系统的质量、阻尼以及刚度等参数可以达到改变结构模态频率和模态振型的目的；错开结构固有频率与流体湍流产生的激励载荷频率，可以提高汽机旁路阀的稳定性。结构固有频率可以通过模态分析得到，模态分析的核心内容是确定结构系统动态特性的模态参数，对于振动系统，其自由振动方程为：Mx tCx tKx t()+()+()=0 （3）式中，M为质量矩阵；C为阻尼矩阵；K为结构刚度矩阵。忽略结构阻尼，得到振动方程：Mx tKx t()+()=0（4）模态分析中假设了结构的线性特性，振动为简谐振动，求解振动方程得到：KMii =20（5）颜炳良，等：常规岛汽机旁路调节阀流激共振研究32FLUID MACHINERYVol.51，No.9，2023每个特征值相应的特征向量是i，其表示为模态振型。方程的根为i2，即特征值，其平方根为i，是结构自然圆周频率，通过圆周频率得到结构自然频率为：fii=/2（6）由式（3）（6）可知，结构模态分析即求解振动方程的特征值与特征向量。特征值与模态频率相关，特征向量决定结构的模态振型。3.2 模态分析对简化得到的阀芯组件模型分别在流体域和结构域定义流固耦合面，联合ANSYSWorkbench中流场分析（Fluent）、静力学分析模块进行阀芯组件流固耦合仿真分析，在此基础上基于ANSYSWorkbenchModel 模块，采用Lanczons 方法对阀芯组件的固有频率进行分析。通过模态分析得到汽机旁路阀 30%开度阀芯组件前 6 阶模态频率及模态振型分别见表 5。表 5 30%开度阀芯组件前 6 阶模态频率Tab.5 Thefirstsixordersofmodalfrequenciesat30%openingofthevalvecoreassembly阶数模态频率/Hz13170.123617.033622.444150.454662.665315.8由表 5 可知，第一阶模态频率为 3243.8Hz，远小于其他阶频率。图 12 示出阀芯组件前 6 阶模态振型。由图 12 可知，主振型表现为阀芯的弯曲、扭转、拉伸等特征，其中 1 阶振型变形主要集中在接近先导阀芯的阀杆部分，表现为阀芯沿Y轴扭转，其他阶振型表现为主阀芯的横向摆动及阀杆的弯曲，其中 6 阶振型变形主要集中在接近先导阀芯的阀杆部分，表现出阀芯沿Z 轴的拉伸。当阀芯组件的固有频率接近阀内流体激振主频率时，极易诱发汽机旁路阀组件振动，若 2 个频率恰好相等，就会发生共振，进而引起流量和压力失调。基于汽机旁路阀流场瞬态计算和流固耦合模态分析，将流体激振主频率和模态频率进行对比分析，如图 13 所示。从图中可知，阀门前、后压差为6.71MPa 时，该型汽机旁路阀在 30%开度时流体激振载荷共振区频率范围为 8481272Hz，随着阶数增大依次递增，对比表 5，其前 6 阶模态频率均不在流体激振共振区内。说明该工况下汽机旁路阀采用多孔与迷宫盘片组合套筒时不会发生流激共振，可为防共振结构设计提供参考。（a）1 阶模态 （b）2 阶模态 （c）3 阶模态 （d）4 阶模态 （e）5 阶模态 （f）6 阶模态图 12 阀芯组件前 6 阶模态振型Fig.12 Thefirstsixordersofmodalvibrationshapesofthevalvecoreassembly图 13 流体激振频率与结构固有频率对比分析Fig.13 Comparativeanalysisoffluidexcitationfrequencyandstructuralnaturalfrequency4 结论（1）设计的汽机旁路阀在 30%开度，压差为 6.71MPa 时，流体激振载荷共振区频率范围为8481272Hz。33（2）阀芯组件前 6 阶模态频率分别为 3170.1，3617.0，3622.4，4150.4，4662.6，5315.8Hz，随着阶数增大依次递增。其中 1 阶振型变形主要集中在接近先导阀芯的阀杆部分，表现为阀芯沿Y 轴扭转。（3）通过激振频率与模态频率对比分析，前六阶模态频率均不在流体激振共振区内。因此汽机旁路阀结构在 30%开度，压差为 6.71MPa 时，采用多孔与迷宫盘片组合套筒时，阀芯组件不发生流激共振，模拟结果可为高压差调节阀防共振结构提供设计依据。参考文献：1 赵福强，祝洪青，易朝晖.海阳核电一期 21253MW机组主要热力系统设计特点 J.节能技术，2011，29（4）：336-339.ZHAOFQ，ZHUHQ，YICH.HaiyangnuclearpowerphaseI2Mainthermalsystemdesignfeaturesof1253MWunit J.EnergySavingTechnology，2011，29（4）：336-339.2 李树勋，王天龙，徐晓刚，等.压降套筒式蒸汽疏水阀振动特性研究 J.振动与冲击，2018，37（4）：147-152.LISX，WANGTL，XUXG，etal.StudyonvibrationcharacteristicsofhighpressuredropsleevesteamtrapJ.VibrationandShock，2011，29（4）：336-339.3 朱禄.先导式调节阀流动及流激振动数值模拟研究D.兰州：兰州理工大学，2019.ZHUL.NumericalSimulationofflowandflow-inducedvibrationinpilotcontrolvalveJ.Lanzhou：LanzhouUniversityofTechnology，2019.4 王伟波，蒋永兵，祝欣慰，等.多孔式套筒调节阀开孔设计及流激共振特性研究 J.仪器仪表学报，2020，41（9）：169-179.WANGWB，JIANGYW，ZHUXW，etal.OrificedesignandStudyonflowresonancecharacteristicsofporoussleeveregulatingvalveJ.JournalofInstrumentation，2020，41（9）：169-179.5 TECZAJ，CHOCHUAG，MOLLR.Analysisoffluid-structureinteractioninasteamturbinethrottlevalveC/AsmeTurboExpo：PowerforLand，Sea&Air，2010.6 ENGEDAA.PerformancestudyandinstabilityanalysisofsteamturbinecontrolvalveC/AsmeTurboExpo：TurbineTechnicalConference&Exposition，2012.7 刘丽，张小斌，邱利民.大流量气体管道中阀门诱发振动机理研究 J.低温工程，2016（4）：50-55.LIUL，ZHANGXB，QIULM.StudyonmechanismofvalveinducedvibrationinlargeflowgaspipelineJ.CryogenicEngineering，2016（4）：50-55.8 徐东涛，葛长榕，孟祥瑞，等.基于FLUENT 软件的多级降压调节阀涡振特性分析 J.流体机械，2022，50（8）：52-58.XUDT，GECR，MENGXR，etal.Analysisofvortexvibrationcharacteristicsofmultistagestep-downregulatorbasedonFLUENTsoftwareJ.FluidMachinery，2022，50（8）：52-58.9 孙丹凤，马咏梅，冯成德，等.气力输送管道固有频率及共振分析J.机械设计与制造工程，2015，44（10）：7-10.SUNDF，MAYM，FENGCD，etal.NaturalfrequencyandresonanceanalysisofpneumaticconveyingpipelineJ.MechanicalDesignandManufacturingEngineering，2015，44（10）：7-10.10 Americansocietyofmechanicalengineersboilerandpressurevesselcode，sectionII，partDS.2019.11 李存程，殷鸣，王玲，等.压气机叶片流固耦合共振判断和谐响应分析J.计算机辅助工程，2020，29（1）：14-18.LICC，YINM，WANGL，etal.ResonancejudgmentandharmonicresponseanalysisofcompressorbladefluidstructurecouplingJ.ComputerAidedEngineering，2020，29（1）：14-18.12 陈文飞，陈虎，王杜，等.生物质锅炉炉管开裂失效分析J.压力容器，2020，37（3）：41-49.CHENWF，CHENH，WANGD，etal.FailureanalysisofboilertubecrackinginbiomassboilerJ.PressureVesselTechnology，2020，37（3）：41-49.13 王晓桥，张永红，徐拥军，等.碳纤维复合材料修复带缺陷储气井的应力分析及试验研究J.压力容器，2020，37（3）：65-71.WANGXQ，ZHANGYH，XUYJ.etalStressanalysisandexperimentalstudyofundergroundnaturegasstoragewellwithdefectsrepairedwithcarbonfibercompositeJ.PressureVesselTechnology，2020，37（3）：65-71.14 田晓洁，谢大帅，刘贵杰，等.基于ANSYS 的气液两相流海洋立管流固耦合特性分析J.振动与冲击，2021，40（7）：260-267.TIANXJ，XIEDS，LIUG，etal.Fluidsolidcouplingcharacteristicsofgas-liquidtwo-phaseflowinmarineriserbasedonANSYS J.VibrationandShock，2021，40（7）：260-267.作者简介：颜炳良（1980），男，硕士，高级工程师，主要从事控制阀的研究开发工作，通信地址：400707 重庆市北培区人民村 1 号重庆川仪调节阀有限公司，E-mail：。本文引用格式：颜炳良，金辉，黄冲，等.常规岛汽机旁路调节阀流激共振研究 J.流体机械，2023，51（9）：27-33.YANBl，JINH，HUANGC，etal.Researchonflow-inducedresonanceofconventionalislandsteamturbinebypasscontrolvalveJ.FluidMachinery，2023，51（9）：27-33.颜炳良，等：常规岛汽机旁路调节阀流激共振研究