CN112966344B - 一种三通弹簧阀阀芯幅频特性计算方法 - Google Patents
一种三通弹簧阀阀芯幅频特性计算方法 Download PDFInfo
- Publication number
- CN112966344B CN112966344B CN202110225223.8A CN202110225223A CN112966344B CN 112966344 B CN112966344 B CN 112966344B CN 202110225223 A CN202110225223 A CN 202110225223A CN 112966344 B CN112966344 B CN 112966344B
- Authority
- CN
- China
- Prior art keywords
- valve
- frequency
- amplitude
- fluctuation
- flow
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Abstract
本发明的目的在于提供一种三通弹簧阀阀芯幅频特性计算方法,包括如下步骤:(1)对三通弹簧阀进行三维建模、网格离散、求解参数的设置;(2)选取波动流量和计算频段,从小到大设置不同弹簧阀入口波动频率,进行单频激励的阀芯振动响应计算;(3)获取阀芯的幅频特性曲线,分析得到该流量工况下的阀芯的幅频特征参数;(4)重复步骤(2)和步骤(3),更改不同的直流量和波动量,进行不同流量工况的幅频特性计算,得到弹簧阀全工况的动态幅频特性参数。本发明能够明晰单频激励工况下阀内的流动特性及动态参数,对优化阀芯“质量‑弹簧”系统的质量、刚度和阻尼提供技术支持,并且可以大幅缩减弹簧阀动态特性的实验成本和研制周期。
Description
技术领域
本发明涉及的是一种管路系统单元动态特性分析方法,具体地说是非对称弹簧阀阀芯的动态特性分析方法。
背景技术
管路系统包括泵、阀、管路等单元,管路系统的运行工作点是由泵的扬程性能曲线和管路阻抗特性曲线(包括阀门和管线)共同决定的。在系统进行流量调节时,泵和管路的工作点发生改变,泵的性能曲线和管路的阻抗特性曲线匹配不合理,导致管路系统出现流动失稳现象,产生流量脉动,该脉动流量经过系统的某些单元被放大后,导致管路系统振动过大,严重影响了管路系统运行的稳定性和可靠性。
管路系统中的三通弹簧阀安装在离心泵后,起到调节系统流量、防止回流的作用,该弹簧阀由弹簧、节流套、阀杆、阀瓣等组成,阀杆为空心筒结构,在其底端和中间位置的周向均布4个矩形通流口,节流套为圆柱套筒结构,固定在空排通道上,其周向也均布4个矩形通流口连通阀杆相应矩形通流口,这些矩形口连通阀内流域和空排通道。在阀瓣开启或关闭的过程中,阀杆和节流套产生相对运动,二者矩形通流口的重合面积发生变化,从而控制由进口流入空排口的流量。在小流量工况上游来流仅由空排口回流到泵入口,主出口没有流体通过,大流量工况阀瓣被顶开,空排口关闭,流体全部由主出口流出,当流体从主出口回流时,在流体力和弹簧力的共同作用下阀瓣快速关闭,防止流体回流。弹簧阀作为一个具有“质量-弹簧”系统的自动调节单元,当来流的脉动频率接近阀芯的“质量-弹簧”系统的固有频率时,导致阀瓣共振,阀瓣产生大幅的往复运动,主出口和空排口的流量波动加剧。主出口的波动流量直接影响下游设备的稳定性,空排口连接到泵入口,造成泵入口的流动失稳,在失稳的来流工况下,泵的振动特性恶化。弹簧阀空排口的波动流量和主出口的波动流量在上下游均导致管路系统流动失稳,影响系统运行的稳定性和可靠性,所以在设计阶段需要弄清弹簧阀在不同来流情况下的幅频特性,避免设备在某些工况出现共振问题。
由于三通弹簧阀的“质量-弹簧”系统起到联调主出口和空排口流量的作用,其动态幅频特性很难通过理论计算得到,现阶段对于这种复杂结构的动态特性的研究主要采用实验方法。实验法需要对阀体和阀芯进行打孔,并布置位移传感器,测量阀芯在不同激励下的幅频特性,然后利用实验数据对阀芯的质量和刚度进行优化,最终得到满足全工况指标的弹簧阀设计。这个过程进行了反复的加工实验,耗费了大量人力物力,而且在一些特殊设计情况下无法对设备进行打孔测试。
发明内容
本发明的目的在于提供突破传统的基于实验的幅频特性测试方法,能够在弹簧阀设计阶段获得波动边界下的动态参数,通过优化“质量-弹簧”系统的结构参数,保证弹簧阀在全工况下满足性能指标,达到弹簧阀定量声学设计目的的一种三通弹簧阀阀芯幅频特性计算方法。
本发明的目的是这样实现的:
本发明一种三通弹簧阀阀芯幅频特性计算方法,其特征是:
(1)对三通弹簧阀进行三维建模、网格离散、求解参数的设置;
(2)选取波动流量和计算频段,从小到大设置不同弹簧阀入口波动频率,进行单频激励的阀芯振动响应计算;
(3)获取阀芯的幅频特性曲线,分析得到该流量工况下的阀芯的幅频特征参数;
(4)重复步骤(2)和步骤(3),更改不同的直流量和波动量,进行不同流量工况的幅频特性计算,得到弹簧阀全工况的动态幅频特性参数。
本发明还可以包括:
1、所述对三通弹簧阀进行三维建模、网格离散、求解参数的设置包括:
对壁面网格进行加密处理以满足湍流模型和近壁面函数对第一层网格高度的要求,并且利用进出口压力和流量的实验数据修正计算模型的参数,通过UDF动网格技术时时输出阀芯的动力学参数,为幅频特性的分析做准备。
2、所述单频激励的阀芯振动响应计算包括:
阀芯简化为一个质量-弹簧系统,其在单频激励下的运动微分方程如下:
其解为:
式中,m表示阀芯质量,c表示质量-弹簧系统的阻尼,k表示刚度,w表示激励力的频率,A表示激励力的在系统刚度下的等效位移,H(w)为质量-弹簧系统对激励力的放大系数,该参数表征系统的幅频特性。
3、所述步骤(3)具体为:将步骤(2)的不同入口流量波动频率的计算结果进行汇总处理,得到弹簧阀阀芯的动态幅频特性曲线,并分析放大系数最大的工况弹簧阀主出口和空排口的流量波动、压力波动和内流场的流动状态。
本发明的优势在于:本发明将动网格技术和幅频特性处理方法相结合运用到CFD数值模拟中,实现了弹簧阀在全工况动态幅频特性的计算,利用CFD的方法能够明晰单频激励工况下阀内的流动特性及动态参数,对优化阀芯“质量-弹簧”系统的质量、刚度和阻尼提供技术支持,并且可以大幅缩减弹簧阀动态特性的实验成本和研制周期。
附图说明
图1a为弹簧阀剖视图,图1b为弹簧阀局部放大图;
图2为阀芯“质量-弹簧”系统示意图;
图3为弹簧阀单频激励计算示意图;
图4为阀芯“质量-弹簧”系统幅频特性。
具体实施方式
下面结合附图举例对本发明做更详细地描述:
结合图1a-4,本发明包括如下步骤:
步骤一、对三通弹簧阀进行三维建模、网格离散、求解参数的设置。
参考图1弹簧阀结构示意图,进行模型建立、网格离散和求解参数的设置等前处理操作。由于前处理对后续阀芯振动响应的计算的准确性有较大影响,本发明对三维模型进行精细化建模,对壁面网格进行加密处理以满足湍流模型和近壁面函数对第一层网格高度的要求,并且利用进出口压力和流量的实验数据修正计算模型的参数,保证计算模型的准确性的前提下,通过UDF动网格技术时时输出阀芯的动力学参数,为幅频特性的分析做准备。
步骤二、根据设计要求选取波动流量和计算频段,从小到大设置不同弹簧阀入口波动频率,进行单频激励的阀芯振动响应计算。
参考图2,阀芯可以简化为一个“质量-弹簧”系统,其在单频激励下的运动微分方程如下:
其解为:
式中,m表示阀芯质量,c表示“质量-弹簧”系统的阻尼,k表示刚度,w表示激励力的频率,A表示激励力的在系统刚度下的等效位移,H(w)为“质量-弹簧”系统对激励力的放大系数,该参数表征系统的幅频特性。
参考图3弹簧阀单频激励计算示意图,根据设计要求选取入口流量波动幅值和波动频率,确定入口流量的直流量和波动量后,波动频率按照从小到大的顺序施加到弹簧阀入口,由于弹簧阀入口施加的正弦流量波动,该波动量引起阀内压力波动,导致阀瓣受到与入口同频的波动力,从而进行该流量工况下的弹簧阀阀芯的动态频响特性的计算。
步骤三、得到阀芯的幅频特性曲线,分析得到该流量工况下的阀芯的幅频特征参数。
将步骤二的不同入口流量波动频率的计算结果进行汇总处理,参考图4阀芯“质量-弹簧”系统的幅频特性,得到弹簧阀阀芯的动态幅频特性,并分析放大系数最大的工况弹簧阀主出口和空排口的流量波动、压力波动和内流场的流动状态。
步骤四、重复步骤二和步骤三,更改不同的直流量和波动量,进行不同流量工况的幅频特性计算,得到弹簧阀全工况的动态幅频特性参数。
Claims (1)
1.一种三通弹簧阀阀芯幅频特性计算方法,其特征是:
(1)对三通弹簧阀进行三维建模、网格离散、求解参数的设置:
对壁面网格进行加密处理以满足湍流模型和近壁面函数对第一层网格高度的要求,并且利用进出口压力和流量的实验数据修正计算模型的参数,通过UDF动网格技术实时输出阀芯的动力学参数,为幅频特性的分析做准备;
(2)选取波动流量和计算频段,从小到大设置不同弹簧阀入口波动频率,进行单频激励的阀芯振动响应计算:
阀芯简化为一个质量-弹簧系统,其在单频激励下的运动微分方程如下:
其解为:
式中,m表示阀芯质量,c表示质量-弹簧系统的阻尼,k表示刚度,w表示激励力的频率,A表示激励力的在系统刚度下的等效位移,H(w)为质量-弹簧系统对激励力的放大系数,该参数表征系统的幅频特性;
根据设计要求选取入口流量波动幅值和波动频率,确定入口流量的直流量和波动量后,波动频率按照从小到大的顺序施加到弹簧阀入口,由于弹簧阀入口施加的正弦流量波动,该波动量引起阀内压力波动,导致阀瓣受到与入口同频的波动力,从而进行该流量工况下的弹簧阀阀芯的动态频响特性的计算;
(3)获取阀芯的幅频特性曲线,分析得到该流量工况下的阀芯的幅频特征参数:
将步骤(2)的不同入口波动频率的计算结果进行汇总处理,得到弹簧阀阀芯的动态幅频特性曲线,并分析放大系数最大的工况弹簧阀主出口和空排口的流量波动、压力波动和内流场的流动状态;
(4)重复步骤(2)和步骤(3),更改不同的直流量和波动量,进行不同流量工况的幅频特性计算,得到弹簧阀全工况的动态幅频特性参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110225223.8A CN112966344B (zh) | 2021-03-01 | 2021-03-01 | 一种三通弹簧阀阀芯幅频特性计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110225223.8A CN112966344B (zh) | 2021-03-01 | 2021-03-01 | 一种三通弹簧阀阀芯幅频特性计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112966344A CN112966344A (zh) | 2021-06-15 |
CN112966344B true CN112966344B (zh) | 2022-08-02 |
Family
ID=76276067
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110225223.8A Active CN112966344B (zh) | 2021-03-01 | 2021-03-01 | 一种三通弹簧阀阀芯幅频特性计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112966344B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102954064A (zh) * | 2011-08-23 | 2013-03-06 | 上海宝钢设备检修有限公司 | 力马达直动式三通伺服阀的检测系统及方法 |
CN112287498A (zh) * | 2020-10-27 | 2021-01-29 | 哈尔滨工程大学 | 一种三通弹簧阀阀芯开度数值模拟研究方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7775086B2 (en) * | 2006-09-01 | 2010-08-17 | Ut-Battelle, Llc | Band excitation method applicable to scanning probe microscopy |
-
2021
- 2021-03-01 CN CN202110225223.8A patent/CN112966344B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102954064A (zh) * | 2011-08-23 | 2013-03-06 | 上海宝钢设备检修有限公司 | 力马达直动式三通伺服阀的检测系统及方法 |
CN112287498A (zh) * | 2020-10-27 | 2021-01-29 | 哈尔滨工程大学 | 一种三通弹簧阀阀芯开度数值模拟研究方法 |
Non-Patent Citations (2)
Title |
---|
A study of a beam-like electromagnetic vibration absorber;Liu Zihao 等;《Journal of Vibration and Control》;20150817;第22卷(第11期);第1-10页 * |
油气分离式弹簧系统的建模与特性分析;熊新 等;《中国公路学报》;20161115;第29卷(第11期);第153-160页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112966344A (zh) | 2021-06-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105909536B (zh) | 一种离心泵气液两相流性能测试系统及其测试方法 | |
CN108087155B (zh) | 一种大流量液体输送系统频率特性试验系统的试验方法 | |
CN109753716B (zh) | 基于流场仿真的核/火电汽轮机组流体激励数值计算方法及系统 | |
CN203837917U (zh) | 一种锥壁孔式不等温进气掺混装置 | |
Wu et al. | Research on cavitation phenomena in pilot stage of jet pipe servo-valve with a rectangular nozzle based on large-eddy simulations | |
CN112966344B (zh) | 一种三通弹簧阀阀芯幅频特性计算方法 | |
EP2927465B1 (en) | Method for determining the instantaneous mass flow rate of a gas, corresponding device and computer program | |
Wang et al. | Structure optimization and flow field simulation of plate type high speed on-off valve | |
CN110909510B (zh) | 一种压力脉动衰减器三维仿真方法 | |
Młynarczyk et al. | The application of nozzles for the attenuation of volumetric compressor pressure pulsation | |
CN110162861A (zh) | 一种基于有限元的涡轮泵充浸液湿模态分析方法 | |
Morita et al. | Flow induced vibration of a steam control valve in middle-opening condition | |
Engeda | Performance study and instability analysis of steam turbine control valve | |
Ruprecht et al. | Simulation of the water hammer in a hydro power plant caused by draft tube surge | |
Cyklis et al. | The CFD based estimation of pressure pulsation damping parameters for the manifold element | |
Ozawa et al. | Estimation of fluid transients in a pipe using kalman filter based on optimized finite element model | |
CN116305608A (zh) | 一种考虑入口管路的三通弹簧阀阀芯动力学数值模拟研究方法 | |
Ashraf Abdelfattah et al. | Analysis of optimum performance of air vessels used in damping water hammer pressure wave | |
Simola et al. | Emulating pump system static head using PID-controlled flow-regulating valve | |
CN108729911A (zh) | 用于资源生产系统的优化装置、系统和方法 | |
Carta et al. | Transfer matrices of single volute centrifugal pumps | |
Xu et al. | Comparative Analysis of Monitoring Methods for Vortex-induced Vibration of Multistage Pressure Reducing Valves | |
JP4881751B2 (ja) | 圧力脈動解析装置 | |
Rammohan et al. | Effect of multi jets on cavitation performance of globe valves | |
Sunden et al. | Smart Flow Control Processes in Micro Scale |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |