气穴现象

基于Fluent的油液减振器气穴现象研究

陈齐平;吴明明;康盛;魏佳成;刘禹;田玥

【摘要】为了抑制油液减振器气穴现象的产生,分析了减振器复原和压缩过程中气

穴的生成过程.基于CFD数值方法,建立了较高精度的减振器三维流体模型和流体网

格模型,在Fluent流体软件中进行了仿真分析,获得了油液减振器气穴产生的具体位

置和分布情况,分析研究了不同温度对气穴产生的影响,并进行了试验验证.结果表

明:气穴现象主要分布于减振器节流阀片周围,随着活塞速度的增大,气穴现象更加明

显;其他因素不变,油液的温度越低,气穴现象更加严重.此方法为降低油液减振器异常

振动噪声提供了一定的参考价值.

【期刊名称】《液压与气动》

【年(卷),期】2019(000)006

【总页数】7页(P35-41)

【关键词】油液减振器;气穴现象;Fluent;分布情况;温度

【作者】陈齐平;吴明明;康盛;魏佳成;刘禹;田玥

【作者单位】华东交通大学载运工具与装备教育部重点实验室,江西南昌330013;

华东交通大学载运工具与装备教育部重点实验室,江西南昌330013;华东交通大学

载运工具与装备教育部重点实验室,江西南昌330013;华东交通大学载运工具与装

备教育部重点实验室,江西南昌330013;华东交通大学载运工具与装备教育部重点

实验室,江西南昌330013;华东交通大学载运工具与装备教育部重点实验室,江西南

昌330013

【正文语种】中文

【中图分类】TH137;U463.213;O351.2

引言

油液减振器是汽车上的重要部件，广泛运用于汽车的悬架系统和转向系统[1]，减

振器通过其内部复原阀产生的油液阻尼力来衰减车轮与车身之间的振动，可以缓和

由不平路面所引起的车身振动和冲击，从而提高汽车乘坐的舒适性。然而当油液减

振器长时间工作在振动冲击强，温度变化大的环境当中时，减振器内部会循环产生

高压和低压区域[2]，加之油液中溶解了一定量的空气，当减振器工作到一定时间，

气穴现象将会随之产生。气穴现象会严重影响油液减振器的工作性能，除产生振动、

噪声外、还会使系统的工作效率降低，缩短减振器的工作寿命[3]，因此对油液减

振器气穴现象进行研究具有重要的意义。

空穴现象的研究是油液减振器噪声防治中不可忽略的一个流程，国内外已经开展了

大量的研究，周长城[4]等建立了油液减振器气穴产生机理的参数模型，提出了减

小气穴产生的阀系结构设计措施，然而没有获取气穴产生的具体位置和分布情况；

周长城[5]等还对油液减振器产生气穴的影响因素和油液使用前、后的物理和化学

特性进行了深入研究，但此研究没有详细分析油液减振器复原和压缩过程中气穴的

生成过程；安成光[6]等通过数值模拟和阻尼参数影响分析，确定了节流气穴现象

在压缩和复原行程出现的状态，但忽视了不同温度对油液减振器气穴现象的影响情

况；LUOTH[7]等采用Simulink仿真分析了油液减振器的压缩回收过程，确定了

缓冲区空化噪声、节流孔长度与油液运动黏度的关系，但此方法没有对油液减振器

的复原过程进行分析和仿真；ECY[8]等基于动力学分析，提出一种油液减振器稳

健设计的方法，并推导出了防止气蚀的设计约束条件，美中不足的是此研究只分析

了减振器节流阀片未开阀时的气穴现象，忽略了开阀后气穴现象的分布情况；

ECRRF[9]等在模拟中采用了单相空化模型，利用CFD数值方法研究了油液减振

器的流动和空化现象，并提出了一种油液减振器结构改进的设计方案，不足的是此

方案没有对CFD的数值方法原理进行深入研究。

针对上述油液减振器气穴现象研究存在的问题，本研究基于CFD的数值方法对油

液减振器的气穴现象进行研究，分析减振器复原和压缩过程中气穴的生成过程，建

立油液减振器三维流体模型和流体网格模型，在Fluent流体软件中进行仿真分析，

并进行试验验证。

1油液减振器的气穴现象

1.1气穴现象产生机理

通常在油液减振器的油液中溶解了一定量的空气[10]，空气的溶解和分解受到油液

温度和减振器内部压力的影响，具体为：温度越高，压力越大，空气的溶解量就会

增加，反之溶解量就会下降。具体的空化现象产生过程如图1和图2所示，由图

1和图2可知，在油液减振器伸张过程中，活塞杆向上推进导致活塞上方A腔的

油液压力变大，空气的溶解度增大，空气大量溶解在油液中；当压力上升到一定数

值时，达到油液减振器复原阀节流阀片的开阀压力，节流阀片逐渐打开阀口，油液

通过复原阀和补偿阀向B腔流动，由于节流阀片处存在节流损失从而导致压力骤

减，当压力下降到油液的饱和蒸气压时，节流阀片周围油液中将会有气体析出，并

逐渐转变成大气泡。由于油液中混有大气泡，导致油液呈现不联系状态，这种现象

称为油液减振器的气穴现象。

1.复原阀(伸张阀)2.压缩阀3.补偿阀4.流通阀5.工作筒6.储油筒7.防尘垫8.导向

座9.油封10.活塞杆图1油液减振器结构

图2气穴现象形成原理

2CFD数值方法

CFD无论具有什么形式，都需要以经典流体力学的3个基本控制方程——连续性

方程、动量方程、能量方程为基础[11]。在实际的计算时，还需要考虑不同的流态，

如层流和湍流等。

2.1连续性方程

在流场中，当流体从一控制面流入控制体时，同时也会从另一控制面流出控制体，

在这一过程中控制体内的流体质量会发生改变[12]。根据物理学上的质量守恒定律，

流入控制体的质量与流出控制体的质量之差会等于控制体内部增加的流体质量，因

此可以得出流体连续性方程的积分形式，流体具体的数学积分形式为：

(1)

式中，ρ为液体密度；V为控制体积；S为控制面积。等式左边第一项表示控制

体内部流体质量的增加率；第二项表示通过控制面流出整个控制体的质量净流量。

在直角坐标下可以将此积分形式化为微分形式：

式中，u、v、w分别表示流体速度的3个分量。

2.2动量方程

在一个给定的流体系统当中，流体动量随时间的变化率等于作用在其表面的外力之

和[13]，这三维是物理学上流体的动量守恒规律，其微分表达式为：

(3)

式中，pxx，pxy，pxz，pyy，pyz和pzz分别表示流体内应力张量的分量；Fbx，

Fby和Fbz表示单位流体质量的质量力在x，y，z3个方向的分量。

2.3能量方程

根据热力学第一定律可得到:流体微团内能量变化率=流入微团内的净热流量+体积

力和表面力对微团做功的功率，所以可以用有关的流体物理量将流体的能量方程表

示出来，具体为:

(4)

式中，为流体流动的速度；keff是有效热传导系数，keff=k+kt；其中kt为湍流

热传导系数；Jj′是组分J的扩散流量；Sh包括了化学反应热以及其他用户自定义

的体积热源项；方程右边3项分别描述了热传导、组分扩散和黏性耗散带来的能

量输运。

3油液减振器建模分析

3.1模型假设与简化

考虑到减振器三维流体模型计算量庞大，在不偏离实际运用的基础上，对减振器流

体模型进行了假设和简化，具体的简化措施为：

(1)忽略减振器内部油液的可压缩性[14]。忽略一些对减振器阻尼特性影响不大的

结构，如：垫片、密封环和螺母等。

(2)忽略减振器叠加阀片之间的摩擦，忽略活塞杆与导向座之间的摩擦[15]。

(3)减振器储油筒和工作缸做刚体处理，采用减振器的1/4模型进行仿真[16]。

3.2减振器复原阀流体模型

建立如图3所示的复原阀流体模型。由图3可知，根据复原阀流场特点，该模型

由上油腔、核心腔和下油腔构成[17]，将模型的对称工具面设置为了流体的对称面，

复原阀流体模型的上油腔为流体的入口，下油腔为流体的出口，核心腔为油液冲击

复原阀的开阀区域，周围壁面设置为无滑移的壁面。

图3复原阀流体模型

3.3减振器复原阀网格模型

减振器复原阀网格的模型如图4所示，此流体网格模型在ICEM中进行划分。由

图4可知，由于油液减振器核心腔流场变化剧烈，因此减振器的核心腔划分为致

密的四面体网格；上油腔和下油腔流场变化平缓，为了减少计算量，上油腔和下油

腔被划分为六面体网格，网格数量合计为484931个。为了仿真结果更加精确，

在活塞和复原阀之间添加预留液体间隙层[18]。减振器油液的所有的变量建立在拐

点节点处，计算数值准确，计算变量插值时具有较好的二阶精度。

图4复原阀流体网格

3.4材料和边界条件设置

在Fluent流体软件中进行减振器气穴现象分析，设定参考压强p为0[19]，所用

油液的运动黏度γ为13.28mm2/s，密度ρ为880kg/m3，黏度指数vi为196；

求解方法采用标准的k-w模型，由于气泡可以溶解于油液，同时为了确保计算结

果精确，因此使用开启多相流Eulerian模型；设置油液为第一相，空气为第二相，

并设置两相的空化作用，选择的空化模型为Zwart-Gerber-Belamri[20],空气压力

设置为0.0342MPa,其他选项保持默认值；采用速度型流体入口和压力型流体出

口，且入口处水体积分数设为100%；使用UDF(用户自定义函数)的动网格宏对叠

加阀片进行动网格自定义编程，同时使用UDF来实现入口速度大小随时间作抛物

线变化[21]。

3.5CFD求解设置

为了缩短仿真时间，同时得到更加准确的流体仿真结果，有必要对在Fluent中进

行仿真计算的求解设置。

(1)为了仿真可以更好的收敛，同时兼顾仿真的准确性，取松弛因子为0.6；

(2)采用适用于瞬态流体分析的PISO算法对仿真进行求解，提高仿真结果的精确

性[22]；

(3)仿真总时间设置为1s，时间步长0.001s，总仿真步数为1000步。

4CFD仿真结果分析

4.1空穴现象分布位置分析

当t=0.21s,v=0.37m/s时，油液减振器具体的气相分布云图和示功图如5所示，

由图5的气相分布云图可知，由于油液冲击节流阀片，节流阀片开始有轻微的开

阀，由于节流损失的存在性，节流阀片和活塞孔附近的压力下降[23]，导致油液减

振器的气穴现象首先在节流阀片和活塞孔附近生成，但气穴现象不明显，由图5

的示功图可知，示功图比较饱满，油液阻尼力F较小，大约为60N，此时油液减

振器工作正常；当t=0.48s,v=0.62m/s时，具体的气相分布云图和示功图如图6

所示，由图6的气相分布云图可知，此时节流阀片已部分开阀，由于油液速度和

节流损失的增大，导致油液压力下降变快，因此气穴现象较明显，同时气穴沿着节

流阀片上边缘分布，由图6的示功图可知，示功图开始产生畸变，油液阻尼力F

较大，大约为400N，此时油液减振器工作稍微异常；当t=1s,v=1.04m/s时，

油液减振器具体的气相分布云图和示功图如图7所示，由图7的气相分布云图可

知，节流阀片已完全打开，油液节流损失达到最大，节流阀片处的压力降到最小，

同时油液速度也达到最大，此时气穴现象最明显，气穴布满节流阀片周围，且有向

下腔蔓延的趋势，由图7的示功图可知，示功图已畸变，油液阻尼力F达到最大，

大约为800N，此时油液减振器工作异常，工作环境恶劣，油液减振器极易损坏。

图8为实验中测试的气穴分布状况，由图8可知，油液减振器气穴在节流阀片周

围分布，且随着活塞杆速度增大气穴现象更加严重，验证了仿真结果。

图5油液减振器气相分布云图和示功图(t=0.21s,v=0.37m/s)

图6油液减振器气相分布云图和示功图(t=0.48s,v=0.62m/s)

图7油液减振器气相分布云图和示功图(t=1s,v=1.04m/s)

图8实验中测试的气穴分布状况

4.2温度对气穴现象产生的影响

减振器气穴现象的产生受到多因素的影响，其中影响最大的因素为温度[24]。油液

减振器油液中通常溶解了一定量的空气，当油液的压力低于空气的饱和蒸气压时，

油液中的气体将会逐渐散发出来，此时气穴现象随之产生。其中油液的压力变化受

到油液温度的影响，具体为：在其他因素恒定的情况下，温度越高，油液压力越大，

则不容易产生气穴现象；温度越低，油液压力越小，则气穴现象容易产生[25]。当

车辆行驶在凹凸不平的路面上时，油液减振器的活塞杆速度会越来越大，此时减振

器油液温度也会越来越高，导致油液的压力越来越大，气穴现象也将会逐渐消失，

但减振器长期工作在高速情况下很容易损坏。

为了验证上诉理论，在Fluent流体软件中进行减振器气穴现象分析，为更好的分

析温度对油液减振器气穴现象产生的影响，在速度v=1.04m/s，温度分别为

T=300K、T=330K和T=350K的情况下进行了多相流仿真，此时气穴现象分布

明显，更容易进行分析和研究。得到的仿真结果如图9所示，图9为油液减振器

在不同温度下的气相分布云图，由图9可知，当T=300K时，油液减振器气穴现

象分布最明显，气穴布满节流阀片周围，且有向下油腔蔓延的趋势。由此证明：在

其他因素恒定的情况下，温度越低，油液减振器的气穴现象越明显，验证了上述理

论。

图9油液减振器不同温度下的气相分布云图

5油液减振器试验分析

为了验证油液温度对减振器气穴现象产生的影响，同时为了验证Fluent多相流仿

真的准确性，进行了减振器示功试验。

本研究的试验设备为上海交通大学仪器工程系生产的型号为QJ-4A-10的减振器

伺服示功机，为了减少试验误差，试验结果采用多个减振器试验的平均值，本试验

在DELLPrecisionT3620工作站上进行计算，具体的实验设备如图10所示。

图10减振器伺服示功机

油液减振器示功图的畸变程度反应了气穴的生成情况，具体为：畸变程度越大，则

油液气穴现象越明显。图11为其他因素一定，油液减振器在油液温度T分别为

290，310，330，350K下的示功图，由图11可知，其他因素一定的情况下，当

油液温度T为290K时，示功图面积最大且产生了严重的畸变，说明此温度下油

液减振器气穴现象最明显；其他因素一定的情况下，当油液温度T为350K时，

示功图面积最小且产生了轻微的畸变，说明此温度下油液减振器气穴现象不明显。

因此验证了其他因素一定的情况下，油液温度越低，油液减振器的气穴现象越明显。

图11油液减振器在不同温度下的示功图

6结论

(1)本研究搭建了较高精度的减振器三维流体模型,获得了油液减振器气穴产生的

具体位置和分布情况，分析研究了不同温度对气穴产生的影响，并进行了试验验证；

(2)气穴现象主要分布于减振器节流阀片周围，随着活塞速度的增大，气穴现象更

加明显；其他因素不变，油液的温度越低，气穴现象更加严重；

(3)采用CFD数值方法可以获得油液减振器气穴产生的分布情况，同时分析了温

度对油液减振器气穴产生的影响，对抑制油液减振器气穴现象的产生和解决油液减

振器的噪声问题具有重要的实际参考价值。

参考文献：

【相关文献】

[1]吕振华,李世民.筒式液阻减振器动态特性模拟分析技术的发展[J].清华大学学

报,2002,42(11):1532-1536.

LVZhenhua,pmentofSimulationandAnalysisTechnologyforDynamic

CharacteristicsofCylin-dricalLiquidDamper[J].JournalofTsinghuaUniversity,

2002,42(11):1532-1536.

[2]董小闵,于建强,杨茂举，等.考虑温度因素的磁流变减振器的优化设计与实验[J].振动与冲

击,2016,35(8):54-59.

DONGXiaomin,YUJianqiang,YANGMaoju,lDesignandExperimentof

MagnetorheologicalDamperConsi-deringTemperature[J].VibrationandImpact,

2016,35(8):54-59.

[3]SACRAMENTOG,tionToolforShockAbsorberNoiPredictionin

TimeandFrequencyDomains[J].InternationalJournalofVehicleNoiandVibration,

2007,3(3):217.

[4]周长城,王国丽,徐伟.减振器液压气穴与气体反弹现象的研究[J].液压与气动,2007,(8):3-7.

ZHOUChangcheng,WANGGuoli,chonHydraulicCavitationandGas

ReboundofShockAbsorber[J].ChineHydraulics&Pneumatics,2007,(8):3-7.

[5]周长城,李迪.减振器液压气穴机理分析与研究[J].起重运输机械,2008,(1):53-56.

ZHOUChangcheng,isandRearchonHyd-raulicCavitationMechanismof

ShockAbsorber[J].LiftingandTransportationMachinery,2008,(1):53-56.

[6]安成光,曹阳,张建武.双筒式液压减振器节流孔气穴现象和噪声分析[J].上海交通大学学

报,2018,52(3):297-304.

ANChengguang,CAOYang,tionPhe-nomenonandNoiAnalysis

ofThrottlingHoleofDouble-cylinderHydraulicShockAbsorber[J].JournalofShanghai

JiaotongUniversity,2018,52(3):297-304.

[7]LUOTH,JINRC,JIANGJ,licModelandCavitationsofHydraulicBuffer[J].

JournalofChongqingJiaotongUniversity,2014.

[8]ECY,SHL,TWY,cAnalysisofaDouble-tubeShockAbsorberforRobust

Design[J].JsmeIntern-ationalJournalSerCMechanicalSystemsMachineElements&

Manufacturing,1997,40(2):335-345.

[9]HERRF,ngShockAbsorberCavitationwithOptimizedDesignbyUsing

NumericalModelingofOilFlow[J].ArchiveofAppliedMechanics,2017.

[10]SYRAKOSA,DIMAKOPOULOSY,ticalStudyoftheFlowina

FluidDamperContainingHighViscositySiliconeOil:EffectsofShear-thinningand

Viscoelasticity[J].PhysicsofFluids,2018,30(3):030708.

[11]BHUYAND,KUMARK.3DCADModellingandCom-putationalFluidAnalysisof

PistonValveofTwinTubeShockAbsorbers[J].MaterialsTodayProceedings,

2017,4(8):7420-7425.

[12]nearContactAnalysisofanAutomotiveShockAbsorber

ShimsUsingFluidStructureInteractionTechnique[J].IndiaAltairTechnologyCon-ference,

2013.

[13]石运序，张建旭，刘源，等.基于Fluent的微型排气回收高效节能涡轮设计[J].液压与气

动,2017(9):38-41.

SHIYunxu,ZHANGJianxu,LIUYuan,ofaMiniatureExhaustRecovery

EfficientandEnergy-savingTurbineBadonFluent[J].ChineHydraulics&Pneu-matics,

2017(9):38-41.

[14]FUKUMORIY,HAYASHIR,OKANOH,nCoupledShockAbsorberSystem

UsingFourElectromagn-eticDampers[C]//l

ofPhysicsConferenceSeries,2016.

[15]税永波，戎红俊，丁渭平.汽车液压减振器抗异响稳健性设计[J].液压与气动,2018(9):20-25.

SHUIYongbo,RONGHongjun,nessDesignofAutomobile

HydraulicShockAbsorberAgainstAbnormalNoi[J].ChineHydraulics&Pneumatics,

2018(9):20-25.

[16]SHAMSM,EBRAHIMIR,RAOUFIA,-FEAAnalysisofHudraulicShock

AbsorberValveBehavior[J].InternationalJournalofAutomotiveTechnology,

2007,8(5):615-622.

[17]WANGR,GUF,CATTLEYR,ing,TestingandAnalysisofaRegenerative

HydraulicShockAbsorberSystem[J].Energies,2016,9(5):386.

[18]贺李平,顾亮,龙凯,等.基于流—固耦合的汽车减振器动态特性仿真分析[J].机械工程学

报,2012,48(13):96-101.

HELiping,GULiang,LONGKai,tionAn-alysisofDynamicCharacteristicsof

AutomotiveShockAbsorbersBadonFluidSolidCoupling[J].ChineJournalof

MechanicalEngineering,2012,48(13):96-101.

[19]YANGJ,HUH,LIUK,lAnalysisofFluid-solidCouplingofaHigh-speed

MarineEngineExhaustPipeBadonTransition[J].ShipEngineering,2017.

[20]earModelingwithCavitationPheno-menonofaNovelShock

AbsorberforAbove-kneeProst-hesis[J].InternationalJournalofInformationAcquisi-tion,

2010,7(3):243-257.

[21]刘华坪,陈浮,马波.基于动网格与UDF技术的阀门流场数值模拟[J].汽轮机技

术,2008,50(2):106-108.

LIUHuaping,CHENFu,calSimulationofValveFlowFieldBadon

DynamicGridandUDFTec-hnology[J].SteamTurbineTechnology,2008,50(2):106-108.

[22]LIANGL,TIANL,ZHANGYQ,beShockAbsorberThermo-mechanical

CouplingSimulation[J].AdvancedMaterialsRearch,2012,(566):669-675.

[23]NOCHNICHENKOI,teristicsofThrottlesinHydraulicShock

AbsorberConsideringTemperatureCh-angesofFluid[J].MechanicsandAdvanced

Technolo-gies，2017，80(2):39-44.

[24]CHENK,YUXP,ZHENGHM,sticThermodynamicModelandOil

TemperatureChangeofHydraulicShockAbsorber[J].ChineJournalofComputational

Mechanics,2017,34(4):523-528.

[25]NIASARM,CLEMENCER,WANGX,ofTemperatureonSurfaceDischarge

inOil[C]//,2012:96-99.