调压阀流体空化特性仿真及影响因素分析.docxVIP

? ? 调压阀流体空化特性仿真及影响因素分析 ? ? 林腾蛟，陈少勋，赵俊渝 （1.重庆大学 机械传动国家重点实验室，重庆 400044；2.重庆齿轮箱有限责任公司，重庆 402263） 0 引言 调压阀是液压系统中的核心控制元件，用于控制系统压力、流量和流动方向，对液压系统性能好坏起到关键作用。在调压阀内部复杂流道和较高流速作用下，内流场易发生空化现象，导致流体产生空化气泡。空化气泡形态和空间位置的变化将改变流场特性，引起调压阀异常振动和噪声，降低阀门的稳定性和可靠性，从而影响液压系统性能，因此有必要对调压阀流体空化现象开展研究。 为了充分了解流场空化机理和空化流动特性，国内外学者针对空化现象建模开展了一系列研究。SINGHAL等［1］采用均相流假设建立了一种基于输运方程的流场空化模型。黄彪等［2］分别采用Kubota空化模型和Singhal空化模型进行空化流动的数值模拟，结合试验得出后者能更准确模拟空化流动。宋玉等［3-4］引入混合密度修正函数，并考虑气体蒸发和凝结影响对空化模型进行修正，所得计算结果与试验数据相符。刘德民等［5］考虑不凝结气体和表面张力等因素，对水轮机流场空化过程进行数值模拟，试验验证了模型的合理性。 在调压阀狭窄流道中，流体以较高速度流动，导致流场呈现高雷诺数湍流流动［6］，因此对湍流模型的选取将直接影响空化计算结果的准确性。CHATTOPADHYAY等［7］基于标准k-ε和Realizable k-ε湍流模型，仿真分析了滑阀内流场的空化现象。COUTIER等［8］针对是否考虑流体可压缩性，采用不同湍流模型对空化现象进行流场模拟，指出了流体可压缩性对湍流模型的影响。王松林等［9］考虑流体可压缩性影响，修正湍流模型，计算结果与试验数据吻合良好。 利用上述空化模型和湍流模型，诸多学者对阀内流体空化特性开展了研究。ANTONI等［10-20］分别研究了开度、背压、阀门结构参数对阀内流场压力分布、流速分布、湍流漩涡和空化特性的影响。 目前，关于调压阀流场空化现象研究已取得了一定进展，但多为二维计算模型且对空化气泡形态演变和空化流动特性研究还有待深入。本文以船用二级调压阀为研究对象，建立三维内流道多相流瞬态仿真模型，分析流体空化流动形态及流场特性，研究开度、流量和背压对调压阀流体空化现象的影响规律。 1 数值理论模型 1.1 控制方程 在调压阀内流场中，忽略两相间的相对滑移速度，将液相与空泡相作为混合流体相进行研究。流体的空化流动控制方程［1，16］如下： （1）连续性方程。 空泡相： 液相： 混合流体相： 式中 αv，αl——空泡相与液相体积分数； ρv，ρl——空泡相与液相密度，kg/m3； Re，Rc——蒸汽产生速率和蒸汽凝结速率； ρm——混合流体密度，kg/m3； （2）动量守恒方程。 式中 P ——压力，Pa； μm——混合流体动力黏度，Pa·s。 1.2 空化模型 基于 Singhal空化模型［1］，考虑流体相变、不凝结气体、气动阻力和空泡表面张力等因素影响，蒸汽产生速率Re和蒸汽凝结速率Rc可表示为： 式中 Ce，Cc——蒸汽产生系数和蒸汽凝结系数； k ——湍流动能，m2/s2； σ ——表面张力，N/m； Pv——修正的饱和蒸汽压，Pa； fv——蒸汽质量分数； fg——不凝结气体质量分数。 1.3 湍流模型 调压阀节流口处雷诺数大于临界值，流动形式为湍流。因此，本文采用标准k-ε湍流模型［20］进行计算，湍流输运方程为： 式中 ρ ——流体密度，kg/m3； xi，xj——位移分量，m； ui——流体速度，m/s； μ，μt——层流黏度和湍流黏度，Pa·s； σk——湍流动能的湍流普朗特数； Gk—— 由平均速度梯度引起的湍流动能k的生成项； Gb——浮力引起的湍流动能k的生成项； ε ——湍流耗散率，m2/s3； σε——湍流耗散率的湍流普朗特数； YM—— 可压缩湍流中脉动扩张引起总耗散率的波动； Sk，Sε——根据计算工况定义的源项； C1ε，C2ε，C3ε——经验系数。 2 计算模型及边界条件 2.1 计算模型 针对图1所示的调压阀结构建立三维实体模型，并抽取流体区域，而后对流域模型进行网格划分。为提高计算精度并降低计算时长，采用六面体网格进行模型离散化处理。考虑到节流口速度和压力变化较大，对该区域进行网格加密。流体计算域网格共计1 078 648个节点，722 848个单元。图2示出了调压阀流场计算域模型。 图1 调压阀结构示意Fig.1 Structural diagram of pressure regulating valve 图2 调压阀三维流场计算域模型Fig.2 3D flow field calculation domain model of pressure