船用三通调节阀流-固耦合噪声数值模拟.pdf

第 43 卷 第 1 期Vol.43,No.12024 年 1 月Journal of Applied AcousticsJanuary,2024 研究报告 船用三通调节阀流-固耦合噪声数值模拟施红旗1周 旭1周爱民1张博浩2李树勋2(1 武汉第二船舶设计研究所武汉430205)(2 兰州理工大学石油化工学院兰州730050)摘要：针对船用PN10DN32三通调节阀噪声声压频谱、声指向性等声学特性规律不明确，噪声声压级是否满足使用要求的问题，基于流-固耦合理论，同时考虑流-固耦合面及流体域内的脉动声学激励源，开展阀门噪声数值模拟研究。分别对三通调节阀在80%及60%开度阀外1 m处的噪声进行数值模拟，分析研究噪声声压频谱特性及声指向性规律。结果表明：80%及60%开度下的噪声声压级分别为49.14 dB(A)、50.79 dB(A)，均小于60 dB(A)的噪声限制，满足使用要求。该文为船用三通调节阀噪声数值模拟提供了理论及方法参考。关键词：三通调节阀；流-固耦合；声指向性；声压频谱特性；噪声数值模拟中图法分类号:TH134文献标识码:A文章编号:1000-310X(2024)01-0142-09DOI:10.11684/j.issn.1000-310X.2024.01.017Numerical simulation of fluid-solid coupling noise in marine three-waycontrol valveSHI Hongqi1ZHOU Xu1ZHOU Aimin1ZHANG Bohao2LI Shuxun2(1 Wuhan Second Ship Design and Research Institute,Wuhan 430205,China)(2 School of Petrochemical Engineering,Lanzhou University of Technology,Lanzhou 730050,China)Abstract:Aiming at the problem that the acoustic characteristics of Marine PN10DN32 three-way regulatingvalve noise spectrum and acoustic directivity are not clear,and whether the sound pressure level of noise meetsthe requirements of use,based on the theory of fluid-structure interaction,considering the fluid-structureinteraction surface and the pulsating acoustic excitation source in the fluid domain,the numerical simulationof valve noise is carried out.The noise at 1 m outside the 80%and 60%opening of the three-way regulatingvalve was numerically simulated,and the characteristics of noise sound pressure spectrum and the law of sounddirectivity were analyzed.The results show that the sound pressure levels of noise at 80%and 60%opening are49.14 dB(A)and 50.79 dB(A)respectively,which are both lower than the noise limit of 60 dB(A)and meet therequirements of operation.It provides theoretical and methodological reference for the numerical simulation ofMarine three-way regulating valve noise.Keywords:Three-way control valve;Fluid-solid coupling;Acoustic directivity;Sound pressure spectrumcharacteristics;Numerical simulation of noise2022-09-08收稿;2023-01-18定稿作者简介:施红旗(1973),男,湖北武汉人,本科,研究员,研究方向:船舶辅机系统声学。通信作者 E-mail:第43卷 第1期施红旗等：船用三通调节阀流-固耦合噪声数值模拟1430 引言随着我国对船舶总体性能要求的提升，船用阀门作为船舶使用装备中重要的压力元件，学者们对其技术要求也不断提高。三通调节阀作为船用阀门中重要的调节元件，是保障系统正常安全运行的关键所在，因此准确分析实际工作工况下的噪声，对三通调节阀乃至整个系统安全运行具有重要意义。近年来，国内外学者对阀门噪声进行了较多研究。文献12做了超临界多级套筒调节阀空化抑制模拟研究，另外针对蒸汽疏水阀在高温高压工况的过热蒸汽导致的振动噪声问题，设计了节流降压消声器；文献3以一种小流量调节阀为研究对象，计算了在内部湍流作用下阀门壁面结构振动导致的外部噪声；文献4用数值模拟与实验验证的办法研究了阀门应力与振动激励的关系；文献5提出了一种分析流动引起的噪声问题的仿真技术，该技术结合了计算流体动力学(Computational fluiddynamics,CFD)、边界元法和气动声学时间反转源定位方法确定了涡旋脱落频率下气动声源的性质；文献6 进行了实验和模拟，发现产生的噪声并不取决于流入速度条件，而是取决于管道的声学和振动频率特性。以上研究分析了相关声学理论及两通阀噪声数值模拟，但是尚无对于合流型三通阀的噪声数值模拟方法研究。为此，基于流-固耦合理论，同时考虑流-固耦合面及流体域脉动声源，开展船用三通调节阀实际工作工况下80%及60%开度下的噪声数值模拟研究。1 三通调节阀噪声理论分析机械振动噪声、流体动力噪声、流体空化噪声是阀门振动噪声的主要来源，由于阀门内部零件加工、装配精度均满足设计要求，故该阀不会产生机械振动噪声。三通调节阀介质为淡水，主通与旁通压力1 MPa左右，远高于淡水的饱和蒸气压力，故不会产生空化噪声。综上可得三通调节阀研究过程中主要考虑流体动力噪声。流体动力噪声指由流体中物体的阻碍作用导致流体流场剧烈变化从而产生辐射的噪声称之为流体动力噪声。噪声声源体为流-固耦合作用的偶极子。偶极子声源是指两个距离很近、相位相反，并以相同幅值振动的单极子声源构成。模型图如图1所示。+-图1偶极子声源模型图Fig.1 Dipole sound source model diagram偶极子声源产生声压的空间部分可表示为p(r,)=jkcAejkr4r2(1+jkr)cos,(1)其中：为介质密度；c为声波传播速度；为相位角；r为声源半径。kr=1tan,(2)式(2)中：kr为方便计算引入的相位角。三通调节阀在偶极子声学激励下，噪声产生机理如图2所示。图中AWPF(Acoustic wall pressurefluctuation)指声学压力波动；TWPF(Turbulentwall pressure fluctuation)指湍流压力波动；阀门在工作时，流体介质处于湍流状态。流体介质向四周传播，辐射出流动噪声与阀及阀控管系内壁面声-固耦合产生振动噪声。此外，流体介质除了内部湍流产生的噪声外，还包括流体介质与阀及阀控管系内壁面的振动辐射噪声。综合考虑以上两种声源，采用流-固耦合理论开展三通调节阀噪声数值模拟研究。流-固耦合技术路线图如图3所示。?AWPFAWPF?图2三通调节阀噪声产生机理Fig.2 Noise generation mechanism of three-wayregulating valve2分析模型2.1三通调节阀结构模型建立三通阀的模型结构，三通调节阀的技术要求如表1所示。模型三维视图如图4所示。1442024 年 1 月?CFD?PN10DN32?CFD?CFD?80%?60%?PN10DN32?图3流-固耦合方法路线图Fig.3 Fluid-solid coupling method roadmap表1三通调节阀技术要求Table 1 Technical requirements of three-way regulating valve名称技术要求工作介质淡水工作压力PN10流量及控制特征三通合流流量特性主控制回路等百分比旁控制回路线性噪声要求6 60 dB(A)1.?;2.?;3.?;4.?;5.?45312图4三通调节球阀三维模型Fig.4 Three-dimensional model of three-way reg-ulating ball valve2.2三通调节阀分析工况及材料参数三通调节阀实际工作工况参数和主要材料性能参数分别见表2、表3所示。表2三通调节阀实际工作工况参数Table 2 Actual working parameters ofthree-way regulating valve公称通径公称压力适用介质DN32PN10淡水表3三通调节阀阀体材料性能参数Table 3 Material performance parametersof three-way control valve body零件阀体/端盖法兰球体/阀杆材料CF800Cr22Ni5Mo3N密度/(kgm3)80307850弹性模量/GPa195210泊松比0.310.333三通调节阀实际工作工况下80开度噪声分析3.1基于k-模型三通调节阀80开度定常流动分析为实现高精度CFD流场模拟，在保证计算精度的前提下，对三通调节阀的三维模型进行合理简化。在80%开度下基础上，阀前与阀后分别增加2倍与6倍公称直径长度的管道，保证流体域充分发展。根据三通阀模型进行反向建模得到流道模型，并根据四面体和六面体各自的网格优势对流道模型进行网格划分，同时对阀芯与流体接触区域的网格局部第43卷 第1期施红旗等：船用三通调节阀流-固耦合噪声数值模拟145加密处理。阀芯流道模型网格如图5所示。以出口压力为目标进行网格无关性检验。划分网格单元数从1137556到1835247，计算3种不同网格的质量流量值。从网格1到网格2压力的变化值为3.7%，从网格2到网格3压力的变化值为0.93%，压力变化值很小，最终流道全开网格节点数为257536，单元数为1538455。图5阀芯流道模型网格剖视图Fig.5Grid profile of valve core flow channelmodel使用有限元仿真模拟软件进行流场仿真计算，具体边界条件数值如表4所示。采用标准 k-双方程湍流模型，对湍流方程的离散格式采用highresolution格式。表4三通调节阀实际工作工况80%度工况参数Table4Three-waycontrolvalveactualworking condition 80%opening conditionparameter table工况实际工作工况开度80%主通压力P1/MPa1.08旁通压力P2/MPa1.07质量流量Q(t/h)3三通调节阀80%开度下的内部流动，其压力、速度云图分别选取ZX 截面云图，三维流线云图选取斜二测云图，具体如图6、图7、图8所示。由图6可知，三通调节阀在实际工作工况80%开度下，由于阀芯节流作用压力波动较大，局部压降达到0.02 MPa，阀后压力分布均匀。由图7可知，阀门在阀芯节流处流速最大，约为5.46 m/s。结合图8速度流线图可知，阀芯内部节流及流道出口区域速度流线分布混乱，流体介质处于湍流，并在阀芯内部流域有涡流产生。因此，三通调节阀阀芯内部及流道出口区域的湍流涡旋，是三通调节阀产生流致噪声的主要来源。?/MPa图6ZX 截面压力分布图Fig.6 Pressure distribution map of ZX section?/(mSs-1)图7ZX 截面流速分布图Fig.7 Velocity distribution of ZX section?/(mSs-1)图8三维流线斜二侧视图Fig.8 Three-dimensional streamline oblique two-side view3.2三通调节阀80开度非定常流动计算在定常流动收敛的基础上，以三通调节阀实际工作工况下80%开度CFD定常流动计算结果为初始条件，数值模拟该工况下的CFD非定常流动。选择时均N-S方程为流体流动基本控制方程，采用可实现k-双方程为湍流模型，选择标准壁面函1462024 年 1 月数。边界条件设置主通压力进口为0.58 MPa和流量出口为0.3961 kg/s，旁通压力进口为0.565 MPa和流量出口为0.6393 kg/s，壁面采用无滑移壁面。迭代收敛控制方程采用均方根(Root mean square,RMS)方法，精度设置为10 106，参考压力为大气压0.1 MPa。非定常流动模拟时，时间间隔取1 103s。选择Acoustic声学求解模型，并导出*.CGNS格式的声场信息(时域压力脉动)，作为流-固耦合模拟三通调节阀流致噪声的声学激励。3.3三通调节阀流-固耦合噪声数值模拟方法验证综合流-固耦合三通调节阀噪声声学数值模拟方法，及文献7中通过50A47H-16C型微启式弹簧载荷安全阀噪声试验验证了流-固耦合声学数值模拟方法的正确性与精确度，噪声试验现场图如图9所示，试验结果表明数值模拟噪声声压级与试验实测噪声声压级的偏差在2%以内，这也间接验证了提出的大口径管线球阀数值模拟内漏噪声仿真方法的可行性与精确性。由表5可知，DN250三通调节阀经实验测得实测噪声声压级与数值模拟噪声声压级仅相差1.85 dB(A)，偏差为3.5%，符合标准规定的5 dB(A)精确度要求。此外，数值模拟噪声声压级略小于实测噪声声压级，这是由于数值模拟时并未虑外部环境噪声。故DN32三通调节阀数值模拟结果可信。图9噪声试验现场图Fig.9 Field diagram of noise experiment表5DN250三通调节阀试验噪声声压级Table 5Sound pressure level of DN250three-waycontrolvalveexperimentalnoise实验开度试验声压级/dB(A)模拟声压级/dB(A)60%52.9051.053.4三通调节阀实际工作工况下80开度振动分析开展三通调节阀的振动分析，将其三维模型导入LMS Virtual.Lab软件中进行网格划分。80%开度下三通调节阀的结构网格模型如图10所示。以网格映射方式将流场外壁面的压力脉动信息转移到壳体内壁面，并将时域压力脉动信息进行傅里叶变换，转换后的采样频率为20 5000 Hz。xyz图1080%开度下阀门结构网格模型Fig.10Grid model of valve structure at 80%opening由图11可知，80%开度下三通调节阀X、Y、Z三向振动峰值频率均出现在4560 Hz，其中X、Y、Z三向振动峰值频率重合。振动幅值大会产生较大的振动噪声。将三向振动峰值频率下的压力幅值(dB)云图列出。80%开度下三通调节阀及阀控管系内壁面峰值频率处的压力脉动幅值(dB)云图如图12示。峰值频率下的阀-阀控管系内壁面压力幅值分布呈阶梯状，主通及旁通进口管道压力脉动幅值较小，在阀芯节流以及阀后管道部位由于流体介质的湍流作用压力脉动幅值较大，这也与CFD定常流动分析结果相吻合。100020003000400050000?/(mSs-2),?0.020.040.060.08?/Hz.X.Y.Z图11三通调节阀振动监测点频域加速度曲线图Fig.11Frequency-domain acceleration curve ofvibration monitoring point of three-way controlvalve第43卷 第1期施红旗等：船用三通调节阀流-固耦合噪声数值模拟147Pressure amplitude dB(RMS).1 occurrence 57238232227221215210204198192187181On boundaryxzy图1280%开度下阀-阀控管道内壁面压力幅值云图Fig.12 Pressure amplitude nephogram of inner wallof valve-valve control pipeline at 80%opening分析结果取振动加速度(m/s2)。取如图13所示三通调节阀及阀控管系振动模型中黄色圆点为振动监测点，再使用LMS软件对流致振动进行定量分析，流致振动信息的耦合面位置是流体和固体接触的固体壁面内部。求解得到振动监测点的X、Y、Z 三向频域振动加速度(m/s2)。图14为三通调节阀80%开度下三向振动峰值频率的振动加速度(m/s2)云图。80%开度下三通调节阀的振动加速度云图取峰值频率4560 Hz。由图14可知，80%开度下三通调节阀在阀后管道及阀前法兰处的部分振动加速度较大，在阀体底部和阀前管道处振动较弱。xyz图13三通调节阀及阀控管系三维振动模型Fig.13Three-dimensional vibration model ofthree-way regulating valve and valve control pip-ing systemTranslational acceleration component.1 occurrence 574.752.530.3011.924.156.388.610.813.115.317.5On boundary/ms2xyz图1480%开度下峰值频率振动加速度云图Fig.14Nephogram of peak frequency vibrationacceleration at 80%opening表明阀门在工作工况时，流体介质流经三通调节阀-阀芯节流部位时，产生的流致振动对节流部位附近的零件及阀后管道影响较大，产生了较大的振动加速度。阀门振动监测点的编号为22779，监测点的频域振动加速度幅值如图11所示。由图14可知，X、Y、Z 三向振动加速度均在4560 Hz 处达到峰值，三向振动加速度幅值分别为0.079 m/s2、0.009 m/s2、0.057 m/s2。3.5三通调节阀实际工作工况下80开度噪声分析建立阀-阀控管系的声场模型，并对其进行网格划分。80%开度下三通调节阀噪声声学分析模型如图15所示，80%开度下三通调节阀噪声分析模型主要包括平面场网格1、平面场网格2、80%开度三通调节阀声学网格、声指向性曲线Line1、声指向性曲线Line2以及6个声压检测点组成。取距离管壁外1 m处点作为声压监测点。采用6个声压监测点求算术平均值的方法获取三通调节阀流致噪声。xyz?图15三通调节阀80%开度下噪声声学分析模型图Fig.15 Acoustic analysis model of noise at 80%opening of three-way control valve将振动数据转移到三通调节阀声学网格上作为声学边界条件，基于Lighthill声类比理论，以Helmholtz声学波动方程计算三通调节阀的流致噪声。对三通调节阀80%开度下流致噪声进行流-固耦合数值模拟。图16为三通调节阀80%开度下阀-阀控管系在不同频率下的声压(dB)分布图。80%开度下三通调节阀声压监测点在1700 Hz、3060 Hz和3780 Hz频率的声压幅值较高，是噪声来源的主要成分。由图16可知，三通调节阀产生的流致噪声经过阀-阀1482024 年 1 月控管系辐射到声场网格上，声压(dB)在声场网格水平面的分布呈现出前后基本对称的趋势，而在竖直面呈现出轴对称的趋势；辐射到声场网格的声压主要集中在阀体节流处与阀门下游管道出口处。Pressure amplitude dB(RMS).1 occurrence 2294.287.881.374.868.361.955.448.942.536.029.5Pressure amplitude dB(RMS).1 occurrence 3984.276.468.660.752.945.137.229.421.613.75.92Pressure amplitude dB(RMS).1 occurrence 4882.975.968.961.854.847.840.733.726.719.612.6On boundaryOn boundaryOn boundaryxyzzxyzxy(a)1700 Hz(b)3060 Hz(c)3780 Hz图16三通调节阀-阀阀控管系不同频率下声压分布图Fig.16 Sound pressure distribution of valve-valve con-trol system at different frequencies of three-way con-trol valve声学分析模型中6个声压监测点的编号分别为937、41034、28489、942、27305、40603。分别对噪声监测点不同频率下的声压级进行A计权，得到20 5000 Hz下声压频谱图，80%开度下三通调节阀不同监测点的声压频谱图如图17所示。由图17可知，各个监测点的声压频谱曲线图不完全一致，但总体声压级分布趋势一致。80%开度下三通阀调节阀的噪声成分在1000 Hz以下声压级较小，在1700 Hz、3060 Hz、3780 Hz频率下噪声声压级较大。根据式(3)编写声压级叠加程序，分别对80%开度下三通调节阀不同声压监测点的声压值进行叠加，得到各个监测点在20 5000 Hz频率范围内的总噪声声压级，详见表6。?/Hz10002000300040005000-140-120-100-80-60-40-2002040?/(NSm-2),dB 40603 27350 942 28489 937 41034图1780%开度下三通调节阀不同监测点的声压频谱图Fig.17 The sound pressure spectrum of different mon-itoring points of three-way control valve at 80%open-ing degree表680%开度下三通调节阀各个监测点处声压级Table 6 Sound pressure level at each moni-toring point of three-way control valve under80%opening工况声压监测点声压级/dB(A)实际工作工况80%开度93749.854103449.612848948.2794248.622735048.804060349.73平均声压级49.14声压级叠加计算公式为LpAe=10lg(10n110+10ni10+10nN10),(3)式(3)中：ni为第i个采样频率点的声压；N 为频率采样数。由表6可知，三通调节阀80%开度下不同声压监测点的声压级差别较小，声压级相差不超过2 dB(A)，取6个声压监测点声压级的算术平均值为三通调节阀80%开度下噪声声压级。三通调节阀80%开度下噪声声压级为49.10 dB(A)，满足要求的噪声低于60 dB(A)的规定，故PN10DN32 三通调节阀在实际工作工况下80%开度的噪声满足要求。第43卷 第1期施红旗等：船用三通调节阀流-固耦合噪声数值模拟149图18为80%开度下三通调节阀及阀控管系在1700 Hz、3060 Hz、3780 Hz频率下的声指向性曲线图。如图18所示，在三通调节阀外部设置两条圆周型声指向性曲线，分别命名Line1和Line2。声指向性曲线按圆周方向360分布，监测点间隔10，共由36个监测点组成。将声指向性曲线的声压值转化极坐标。图18(a)Line1指向性曲线反映了三通调节阀在XOZ 平面的噪声指向性。由图18(a)可以看出，1700 Hz与3780 Hz下三通调节阀噪声声压级在XOZ 平面圆周方向分布均匀；在3060 Hz下，阀门主通管道方向与阀门管道方向，即0与170方向噪声声压值较小。图18(b)Line2指向性曲线反映了三通调节阀XOY 平面的噪声指向性。由图18(b)可以看出，1700 Hz与3060 Hz下声指向性分布均匀；在3780 Hz下，三通调节阀在阀后管道，即0方向噪声声压值大于其他方位声压值。(a)Line1?(b)Line2?0O30O60O90O120O150O180O210O240O270O300O330O01020304050600102030405060?/dB(A)?/dB(A)1700 Hz 3060 Hz 3780 Hz0O30O60O90O120O150O180O210O240O270O300O330O010203040506001020304050601700 Hz2420 Hz3060 Hz图18三通调节阀80%开度下声指向性曲线图Fig.18 Acoustic directivity curve of three-way con-trol valve with 80%opening4三通调节阀实际工作工况下60开度噪声分析采用与80%开度相同的数值模拟方法，分析三通阀60%开度声学特性。得到各个监测点在20 5000 Hz频率范围内的总噪声声压级，详见表7。由表7可知，三通调节阀在实际工作工况下60%开度的噪声同样满足要求。表760%开度下三通调节阀各个监测点处声压级Table 7 Sound pressure level at each moni-toring point of three-way control valve under60%opening工况声压监测点声压级/dB(A)实际工作工况60%开度2868851.212955952.082913950.642911949.872953951.712954949.2平均声压级50.79声压频谱图如图19所示。由图19可知，各个监测点的声压频谱曲线图不完全一致，但总体声压级分布趋势一致。60%开度下三通阀调节阀的噪声成分在500 Hz以下声压级较小，在1700 Hz、3540 Hz、4100 Hz频率下声压级较大。?/Hz10002000300040005000-40-200204060 28688 29559 29139 29119 29539 29549?/(NSm-2),dB图19三通调节阀不同监测点的声压频谱图Fig.19 The sound pressure spectrum of differentmonitoring points of three-way control valve图20为60%开度下三通调节阀及阀控管系在1700 Hz、3540 Hz和4100 Hz频率下的声指向性曲线图。在60%开度三通调节阀外部设置两条圆周型声指向性曲线，分别命名Line3和Line4。声指向性曲线按圆周方向360分布，监测点间隔10，共1502024 年 1 月由36个监测点组成。图20(a)Line3指向性曲线反映了XOZ 平面的噪声指向性。由图20(a)可以看出，1700 Hz、3540 Hz 以及4100 Hz下三通调节阀噪声在XOZ 平面120方向附近噪声声压级较小，其余方向声压级分布较为均匀。图20(b)Line4指向性曲线反映了XOY 平面的噪声指向性。由图20(b)可以看出，1700 Hz下三通调节阀在XOY 平面345方向噪声声压级较大，3540 Hz下三通调节阀在150方向噪声声压值大于其他方向。(a)Line3?(b)Line4?/dB(A)?/dB(A)0O30O60O90O120O150O180O210O240O270O300O330O01020304050010203040501700 Hz3540 Hz4100 Hz0O30O60O90O120O150O180O210O240O270O300O330O0102030405001020304050 1700 Hz 3540 Hz 4100 Hz图20三通调节阀60%开度下声指向性曲线图Fig.20Acoustic directivity curve of three-waycontrol valve with 60%opening5 结论根据以上对PN10DN32三通调节阀在实际工作工况下80%及60%开度的流-固耦合噪声数值模拟分析，可得出如下结论：(1)采用流-固耦合数值模拟方法，求得三通调节阀实际工作工况下80%与60%开度的噪声声压级分别为49.14 dB(A)、50.79 dB(A)。噪声声压级均小于60 dB(A)的噪声限制。(2)由两个方向上的声指向性曲线分析可得：80%开度下，在XOZ 平面，1700 Hz与3780 Hz下三通调节阀噪声声压级在XOZ 平面圆周方向分布均匀，且0与170方向噪声声压值较小；在XOY平面，1700 Hz与3060 Hz下声指向性分布均匀，且0方向噪声声压值大于其他方位声压值。60%开度下，1700 Hz、3540 Hz以及4100 Hz下三通调节阀噪声在XOZ 平面120方向附近噪声声压级较小，其余方向声压级分布较为均匀；1700 Hz下三通调节阀在XOY 平面345方向噪声声压级较大，3540 Hz下三通调节阀在150方向噪声声压值大于其他方向。(3)三通调节阀噪声实验验证证明数值模拟分析方法结果可行性与精确度较高，对今后三通调节阀声学特性研究具有指导意义。参考文献1 娄燕鹏.高压降多级降压疏水阀及阀控管道振动噪声特性研究D.兰州:兰州理工大学,2016.2 李树勋,康云星,孟令旗,等.多级降压疏水调节阀流致噪声数值模拟研究J.振动与冲击,2020,39(14):116121.Li Shuxun,Kang Yunxing,Meng Lingqi,et al.Flow-induced noise of a high pressure drop control valveJ.Journal of Vibration and Shock,2020,39(14):116121.3 孙长周.内部湍流作用下调节阀外部噪声的耦合计算研究D.上海:华东理工大学 2017.4 Yonezawa K,Ogi K,Takino T,et al.Experimental andnumerical investigation of flow induced vibration of steamcontrol valveC/Asme Joint USEuropean Fluids Engi-neering Summer Meeting Collocated with InternationalConference on Nanochannels,2010.5 Croaker P,Mimani A,Doolan C,et al.A compu-tational flow-induced noise and time-reversal techniquefor analysing aeroacoustic sourcesJ.The Journal of theAcoustical Society of America,2018,143(4):23012312.6 Mori M,Masumoto T,Ishihara K.Study on acoustic,vi-bration and flow induced noise characteristics of T-shapedpipe with a square cross-sectionJ.Applied Acoustics,2017,120:137147.7 潘伟亮.微启式弹簧载荷安全阀振动特性与气动噪声预测研究D.兰州:兰州理工大学,2019.