深海液压阀变形与动态特性的分段计算方法研究：基于双非线性流固耦合的全海深阀芯-阀套运动副分析

《Science Progress》：A segmented computational method for deformation and dynamic characteristics analysis of deep-sea hydraulic valve

　　

摘要：液压系统是深海潜水器的主要驱动机制。随着潜水器下潜，环境压力从0上升至110 MPa，诱导阀芯（spool）与阀套（sleeve）产生弹性变形，从而改变阻尼力并损害流动稳定性。为应对这一挑战，研究人员提出了一种分段计算方法（Segmentation Calculation Method, SCM），该方法在预测径向变形时能够高效捕捉压力分布与截面几何形状之间的双重非线性耦合。SCM通过将解析流固耦合（Fluid–Structure Interaction, FSI）集成到单次建模过程中，消除了对迭代多物理场仿真的需求，在保持高保真度的同时降低了计算成本。与三维有限元分析（Finite Element Analysis, FEA）的验证对比显示高度一致，在441个评估点上的决定系数（R² = 0.87），对应径向变形中89.6%的方差被解释。尽管尚缺乏实验验证，但与FEA及计算流体动力学（Computational Fluid Dynamics, CFD）结果的一致性支持了该模型的可靠性。利用SCM，研究人员进一步比较了由钢材、氧化铝（Al₂O₃）陶瓷和氮化硅（Si₃N₄）陶瓷制成的阀芯-阀套运动副的阻尼特性，为全海深液压阀的设计提供了实用见解。

论文解读

一、研究背景与问题提出

液压系统作为现代工程中至关重要的动力传输与控制技术之一，被广泛应用于深海大功率、大规模作业设备中，如水下遥控潜水器（Remotely Operated Vehicle, ROV）和载人潜水器。这些深海装备通过压力补偿器来适应极高的海水环境压力，使液压系统内外压力保持平衡。在全海深环境下，最大环境压力可达110 MPa，对应约11000米的挑战者深渊深度。

在此极端压力条件下，液压阀中的运动副——特别是滑阀中的阀芯（spool）与阀套（sleeve）——会发生显著的弹性变形。阀芯与阀套之间设计有适当的配合间隙以保障相对运动，但在超高环境压力和工作压力的共同作用下，间隙高度的变化直接影响阻尼特性，进而导致深海装备运行平稳性、密封性能以及速度/力输出的下降。此外，阀芯上均压槽的存在使得运动副间隙配合中的压力分布呈现非线性状态，进一步增加了问题的复杂性。

已有研究虽然取得了若干进展，但存在明显局限。部分研究基于特定假设推导了变形公式，如假设滑阀截面直径为常数，却忽略了密封时配合间隙对阻尼特性的影响；另一部分研究关注了油液属性变化引起的阻尼特性变化，却忽略了超高压力引起的间隙配合变形对阻尼特性的影响。因此，针对深海液压系统性能的全面分析，迫切需要开展配合间隙变形与阻尼特性之间关系的深入研究。

二、研究内容与重要意义

本研究以一台具有固定节流孔和可调压力补偿器的双向流量控制阀为对象，该阀在全海深条件下表现出流动不稳定性。研究人员从深海变形和液压流体属性变化两个角度出发，全面探索了该阀的阻尼特性。研究提出了一种双非线性变形的分段计算方法（Segmentation Calculation Method, SCM），获取了不同材料阀芯、阀套及配合间隙的变形规律。同时，基于描述液压油密度和粘度随压力变化的方程，建立了全海深及全工况条件下与变形、液压流体粘度和密度相关的阻尼系数方程。

研究表明，在深海条件下，高性能陶瓷（Al₂O₃和Si₃N₄）的性能退化显著小于合金钢，为全海深双向流量控制阀的动态分析奠定了基础。该研究为深海液压系统的设计和性能优化，特别是在复杂工况下的流动稳定性与材料选择方面，提供了重要参考。

三、关键技术方法

本研究采用的主要技术方法包括：（1）利用ANSYS/Fluent进行流场仿真，建立配合间隙中的压力分布模型，提出压力衰减率（k_p）概念来描述台肩处的压力下降规律；（2）提出分段计算方法（SCM），将阀芯和阀套沿轴向分别划分为19段和30段微元段，在每段内假定等载荷和等截面，通过剪力与弯矩的平衡方程组合各段，基于胡克定律和铁摩辛柯梁理论（Timoshenko Beam Theory）计算径向变形，并引入有效剪切相互作用系数（k_i）处理几何不连续处的剪切流重分布；（3）利用ANSYS软件中SOLID186单元进行三维有限元分析（FEA），在径向间隙区域将网格密度细化至约2 μm，以验证SCM的计算精度；（4）基于Vacca和Franzoni建立的粘度-压力关系，结合变形计算结果，建立全海深阻尼系数方程。研究使用的液压油为ISO VG 10液压油，参考温度为40°C。

四、主要研究结果

（一）配合间隙压力分布模型

通过CFD仿真发现，压力仅在均压槽之间的台肩位置下降，且各处台肩的压力衰减率一致。基于此，研究人员提出了压力衰减率（k_p）的概念，建立了无需依赖复杂仿真即可获取间隙中各位置压力值的解析表达式。在不同环境压力Pa（0、55、110 MPa）和工作压力Pw（31.5、20、10 MPa）条件下，该压力分布模型的计算结果与CFD仿真结果几乎重合，验证了模型的有效性。

（二）分段计算方法（SCM）的构建与验证

针对阀芯与阀套内外表面非线性压力分布与截面非线性变化的双重非线性问题，SCM将阀芯分为19段、阀套分为30段进行受力分析，建立了各段的力与力矩平衡方程。通过与三维FEA在441个评估点上的对比验证，SCM的决定系数R²达到0.87，解释了89.6%的径向变形方差。敏感性分析表明，对杨氏模量和载荷进行10.4%平均相对偏差的扰动后，所有预测结果的变异系数小于1.8%，证实了SCM的高重复性和低输出离散度。FEA与SCM的对比还揭示：随着Pa增大（海深增加），阀芯和阀套的变形更加显著；变形形态随Pa变化呈现不同模式；SCM与FEA之间的绝对差异极小（变形量在1 μm以内）；杨氏模量越小的材料变形越大，且随着杨氏模量增加，SCM与FEA的差异进一步缩小。

（三）全海深变形系数分析

为描述不同海深和工况下变形的整体程度，研究人员定义了变形系数（e），即各分析单元变形量与段长乘积之和。结果表明：随着Pw减小和Pa增大，阀芯的总体膨胀趋势逐渐减小并过渡到压缩；阀套主要表现为压缩趋势，Pw对阀套的影响相对较小；配合间隙的总体变形主要为压缩；材料的杨氏模量越小，变形越显著，但模量持续增加时对变形的抑制效果逐渐减弱。研究人员提出将最小后变形间隙3 μm作为避免接触卡死的标准。

（四）全海深阻尼特性分析

基于多台肩阀芯的阻尼力计算方法和随海深变化的液压油动力粘度方程，研究人员获得了全海深范围内粘性阻尼系数在Pa（0~110 MPa）和Pw（3.5~31.5 MPa）变化下的变化规律。结果表明：随着海深增加（Pa升高），阀芯间隙的粘性阻尼系数呈上升趋势，且随工作压力增加而相应增加；变形较小的材料（如Si₃N₄陶瓷），其粘性阻尼系数对Pa和Pw变化的敏感度相对较低。

五、讨论与结论

本研究通过FEA揭示，局部均压槽中的压力衰减仅发生在台肩处，衰减率与总压降和台肩总宽度之比相关，从而验证了所提出的压力分布模型。为有效处理压力分布与截面形状的双重非线性变化，研究人员提出了SCM方法，该方法通过单次建模显著节省计算时间，同时实现高计算精度，填补了载荷与截面双重非线性变化的液压运动副分析方法的空白。通过SCM计算的阀芯、阀套及配合间隙在全海深各工况下的变形规律，并经过FEA进一步验证。此外，基于多台肩阀芯-阀套间隙配合的阻尼系数计算方法，获取了运动副间隙阻尼系数的变化规律，为建立全海深阀芯运动模型提供了有力支撑。尽管该方法尚缺乏有效的实验验证手段，但考虑到仿真中综合假设以及变形相对于实际尺寸较小的事实，仿真结果能够稳健支持所提出的假设。研究深入探索了深海阀芯-阀套运动副的变形和阻尼变化，为建立全海深全工况条件下的完整运动模型奠定了坚实基础，对深海液压系统性能的全面理解具有积极意义。