基于CFD的高压内啮合齿轮泵数值分析与径向力研究

基于CFD的高压内啮合齿轮泵数值分析与径向力研究一、绪论1.1研究背景与意义在现代工业领域，液压系统作为机械设备的关键动力源，广泛应用于冶金、采掘机械、航空航天和深海探测等诸多行业。齿轮泵作为液压系统中不可或缺的动力元件，承担着将机械能转化为液压能，为系统提供稳定压力和流量的重要任务。其中，高压内啮合齿轮泵凭借其独特的结构和工作原理，展现出诸多优势，在工业生产中发挥着愈发重要的作用。高压内啮合齿轮泵的结构相对紧凑，主要由内齿轮、外齿轮、月牙形隔板等部件组成。内齿轮与外齿轮相互啮合，通过齿轮的旋转实现液体的吸入和排出。与外啮合齿轮泵相比，其具有体积小巧、重量轻的特点，能够适应空间有限的工作环境，为设备的紧凑化设计提供了便利。同时，高压内啮合齿轮泵的流量脉动小，工作过程中液体的流动较为平稳，这使得其在对流量稳定性要求较高的应用场景中表现出色，能够有效减少系统的振动和噪声，提高设备的运行稳定性和可靠性。此外，它还具备良好的自吸能力，能够在较短时间内实现液体的自吸，减少了对外部辅助设备的依赖，降低了系统的复杂性和成本。在输送高粘度液体时，高压内啮合齿轮泵的优势尤为明显，能够实现稳定的流量输出，满足工业生产中不同介质的输送需求。随着经济的快速发展和各工业领域需求的不断增长，对高压内啮合齿轮泵的性能要求也日益提高。当前，其发展趋势呈现为高压力、低流量脉动、低噪声、大排量等特征。在航空航天领域，为满足飞行器液压系统对高压力、小体积的要求，高压内啮合齿轮泵需要具备更高的压力等级和更紧凑的结构；在深海探测设备中，为适应恶劣的工作环境，要求泵具有良好的密封性和可靠性，同时能够在高压下稳定工作。然而，传统的中、低压内啮合齿轮泵难以满足这些高压场合的需求，因此，高压内啮合齿轮泵的开发与设计成为当前齿轮泵领域的重要研究方向。在高压工况下，高压内啮合齿轮泵面临着不平衡径向力过大的问题。不平衡径向力是指在泵工作时，由于齿轮两侧受到的液体压力不均匀，导致齿轮轴受到一个径向的合力。这一合力会对齿轮泵的性能和可靠性产生诸多不利影响。过大的不平衡径向力会加速轴承的磨损，降低轴承的使用寿命，增加设备的维护成本。严重时甚至会使轴发生变形，导致齿顶与泵体内壁产生摩擦，影响泵的正常工作，降低泵的容积效率和水力效率。不平衡径向力还会引发泵的振动和噪声，不仅对工作环境造成污染，还可能影响周围设备的正常运行。因此，准确获取内啮合齿轮泵泵内不平衡径向力的分布规律，合理设计静压支撑的径向力补偿结构，成为确保内啮合齿轮泵高压化工作的关键。数值计算方法在高压内啮合齿轮泵的研究中具有重要意义。传统的齿轮泵设计方法主要依赖于经验公式和实验测试，这种方法不仅耗时费力，而且难以全面深入地了解泵内部复杂的流动特性和压力分布情况。随着计算机技术和计算流体力学（CFD）的飞速发展，数值计算方法为高压内啮合齿轮泵的研究提供了新的途径。通过CFD数值模拟，可以在计算机上建立高压内啮合齿轮泵的三维模型，对泵内部的流场进行详细分析，获取不同工况下泵的压力分布、速度分布、流量脉动等信息。这有助于深入理解泵的工作原理和内部流动机制，为泵的优化设计提供理论依据。数值计算方法还可以快速评估不同结构参数和工作参数对泵性能的影响，减少实验次数和成本，提高设计效率。本研究旨在通过探究适合于内啮合齿轮泵高压工况下不平衡径向力求解的CFD数值计算方法和理论研究方法，对高压内啮合齿轮泵的不平衡径向力分布规律开展深入研究。结合静压支撑槽对不平衡径向力的补偿作用，分析静压支撑槽的位置和结构尺寸对不平衡径向力补偿作用的影响，确定符合不平衡径向力要求的静压支撑结构。通过搭建试验台对高压内啮合齿轮泵进行试验研究，验证数值计算方法的可靠性。本研究成果对于提高高压内啮合齿轮泵的性能和可靠性，推动其在工业领域的广泛应用具有重要的理论和实际意义。1.2国内外研究现状1.2.1高压齿轮泵不平衡径向力的研究现状高压齿轮泵不平衡径向力一直是国内外学者研究的重点问题。在国外，一些学者通过理论分析和实验研究，对齿轮泵不平衡径向力的产生机理和影响因素进行了深入探讨。例如，[国外学者姓名1]通过建立齿轮泵的力学模型，分析了齿轮啮合过程中液体压力的分布情况，揭示了不平衡径向力的产生原因。研究发现，齿轮两侧的液体压力差是导致不平衡径向力的主要因素，而齿轮的啮合方式、齿形参数以及泵的工作压力等都会对液体压力分布产生影响，进而影响不平衡径向力的大小。[国外学者姓名2]通过实验测试，研究了不同工况下齿轮泵不平衡径向力的变化规律，提出了一些减小不平衡径向力的方法，如优化齿轮结构、采用平衡槽等。实验结果表明，优化齿轮结构可以改变液体压力的分布，从而减小不平衡径向力；采用平衡槽则可以通过引入额外的压力平衡机制，有效降低不平衡径向力的大小。国内学者在高压齿轮泵不平衡径向力的研究方面也取得了不少成果。[国内学者姓名1]运用有限元分析方法，对齿轮泵内部的流场和压力场进行了数值模拟，分析了不平衡径向力对齿轮和轴承的影响。模拟结果显示，不平衡径向力会导致齿轮和轴承承受较大的应力，加速其磨损和疲劳破坏。[国内学者姓名2]通过理论推导和实验验证，提出了一种基于静压支撑的不平衡径向力补偿方法，该方法通过在齿轮泵的端盖上设置静压支撑槽，利用液体的静压作用来平衡不平衡径向力，取得了较好的效果。实验数据表明，采用静压支撑补偿方法后，不平衡径向力得到了显著降低，齿轮泵的性能和可靠性得到了有效提升。1.2.2CFD数值计算方法在齿轮泵研究中的应用现状随着计算机技术的飞速发展，CFD数值计算方法在齿轮泵研究中得到了广泛应用。国外一些研究团队利用CFD软件对齿轮泵的内部流场进行了详细模拟，分析了泵的流量脉动、压力分布等性能参数。[国外研究团队名称1]采用动网格技术，模拟了齿轮泵工作过程中齿轮的运动，准确地捕捉到了泵内部流场的动态变化。通过对模拟结果的分析，揭示了流量脉动与齿轮运动、流场结构之间的内在联系，为优化齿轮泵的设计提供了理论依据。[国外研究团队名称2]运用多相流模型，研究了齿轮泵在输送含气液体时的气液两相流动特性，分析了气穴现象对泵性能的影响。研究发现，气穴现象会导致泵的流量和效率下降，同时加剧泵内部的磨损和腐蚀。在国内，许多学者也借助CFD数值计算方法对齿轮泵进行了深入研究。[国内学者姓名3]基于CFD方法，对齿轮泵的端面间隙进行了优化设计，通过模拟不同间隙尺寸下泵的内部流场和性能参数，确定了最优的端面间隙值，提高了泵的容积效率。模拟结果表明，合理的端面间隙可以有效减少内部泄漏，提高泵的容积效率。[国内学者姓名4]利用CFD软件，分析了齿轮泵在不同工况下的压力脉动特性，探讨了压力脉动产生的原因和影响因素，提出了降低压力脉动的措施。研究结果表明，通过优化齿轮的齿形参数、调整泵的工作转速等方法，可以有效降低压力脉动，提高泵的工作稳定性。1.2.3齿轮泵空化两相流的研究现状齿轮泵空化两相流是影响泵性能和寿命的重要因素之一，国内外学者对此开展了大量研究。国外学者[国外学者姓名3]通过实验和数值模拟相结合的方法，研究了齿轮泵空化的发生机理和发展过程，分析了空化对泵内部流场和压力分布的影响。实验观察发现，空化首先在泵的低压区产生，随着空化的发展，会形成大量的气泡，这些气泡在高压区溃灭，产生强烈的冲击力，对泵的内部部件造成损伤。数值模拟结果则进一步揭示了空化过程中流场的变化规律，为预测和控制空化提供了理论支持。[国外学者姓名4]提出了一种新的空化模型，用于模拟齿轮泵内的空化两相流，该模型考虑了气泡的生成、生长和溃灭等过程，提高了模拟的准确性。通过与实验数据的对比验证，证明了该模型能够更准确地预测齿轮泵内的空化现象。国内学者在齿轮泵空化两相流研究方面也取得了一定进展。[国内学者姓名5]运用CFD软件，结合空化模型，对齿轮泵内的空化两相流进行了数值模拟，分析了不同工况下空化的发展情况和对泵性能的影响。模拟结果表明，随着泵的工作压力升高和转速增加，空化现象会更加严重，导致泵的流量和效率下降。[国内学者姓名6]通过实验研究，分析了齿轮泵空化噪声的特性，提出了一些降低空化噪声的方法，如优化泵的结构、改善吸入条件等。实验结果表明，通过优化泵的进口结构，增加进口管径，可以有效降低泵的进口流速，减少空化的发生，从而降低空化噪声。1.2.4内啮合齿轮泵试验研究的现状内啮合齿轮泵的试验研究对于验证理论分析和数值模拟结果、改进泵的性能具有重要意义。国外一些研究机构搭建了高精度的试验台，对内啮合齿轮泵的性能进行了全面测试。[国外研究机构名称1]通过试验研究，分析了内啮合齿轮泵的效率特性、流量脉动特性以及噪声特性，研究了不同结构参数和工作参数对泵性能的影响。试验结果表明，内啮合齿轮泵的效率随着工作压力的升高而降低，流量脉动和噪声则随着转速的增加而增大。[国外研究机构名称2]开展了内啮合齿轮泵的耐久性试验，研究了泵在长期运行过程中的性能变化和失效机理，为提高泵的可靠性提供了依据。试验数据显示，泵的关键部件如齿轮、轴承等在长期运行过程中会逐渐磨损，导致泵的性能下降，最终失效。国内学者也积极开展内啮合齿轮泵的试验研究工作。[国内学者姓名7]搭建了内啮合齿轮泵的试验台，对泵的水力性能进行了测试，验证了数值计算方法的准确性。试验结果与数值模拟结果的对比分析表明，数值计算方法能够较好地预测内啮合齿轮泵的水力性能，但在某些细节方面仍存在一定的误差。[国内学者姓名8]通过试验研究，分析了内啮合齿轮泵的自吸性能，探讨了影响自吸性能的因素，提出了改善自吸性能的措施。研究发现，增加泵的转速、减小进口管路阻力等方法可以有效提高内啮合齿轮泵的自吸性能。1.2.5研究现状总结综上所述，国内外学者在高压齿轮泵不平衡径向力、CFD数值计算方法、齿轮泵空化两相流及内啮合齿轮泵试验研究等方面取得了丰硕的成果。然而，目前的研究仍存在一些不足之处。在不平衡径向力的研究方面，虽然已经提出了多种补偿方法，但对于一些复杂工况下的不平衡径向力补偿效果仍有待进一步提高，且对补偿结构的优化设计还需要更深入的研究。CFD数值计算方法在齿轮泵研究中虽然得到了广泛应用，但计算精度和计算效率之间的矛盾仍然存在，需要进一步改进计算模型和算法。齿轮泵空化两相流的研究虽然取得了一定进展，但对空化的微观机理和多相流的复杂相互作用还需要更深入的探索。内啮合齿轮泵的试验研究虽然验证了一些理论和数值模拟结果，但试验的工况范围和研究内容还不够全面，需要进一步拓展。针对这些不足，本研究将致力于探究适合于内啮合齿轮泵高压工况下不平衡径向力求解的CFD数值计算方法，深入研究不平衡径向力的分布规律，通过试验研究验证数值计算方法的可靠性，为高压内啮合齿轮泵的设计和优化提供更有力的支持。1.3内啮合齿轮泵高压化关键技术高压内啮合齿轮泵的整体结构较为紧凑，主要由内齿轮、外齿轮、月牙形隔板、泵体、端盖等部件组成。内齿轮与外齿轮相互啮合，且内齿轮通常为主动齿轮，由电机等驱动装置带动旋转，外齿轮则随之同步转动。月牙形隔板位于内、外齿轮之间，其作用是将吸油腔和压油腔分隔开，确保泵在工作过程中油液能够按照既定的路径流动，实现吸油和压油的功能。泵体作为齿轮泵的外壳，为内部部件提供支撑和保护，同时与端盖配合，形成密封的工作腔室，防止油液泄漏。端盖则安装在泵体的两端，起到封闭工作腔室、固定轴承和密封的作用。在高压工况下，内啮合齿轮泵的间隙泄漏问题较为突出，会严重影响泵的容积效率和工作性能。为解决这一问题，常采用间隙补偿措施。一种常见的间隙补偿方式是采用浮动侧板结构。浮动侧板通常安装在齿轮的端面，与齿轮端面之间形成微小的间隙。在泵工作时，通过引入高压油作用在浮动侧板的背面，使其产生一个朝向齿轮端面的推力，从而使浮动侧板与齿轮端面紧密贴合，减小端面间隙，降低泄漏量。这种方式能够根据泵的工作压力自动调整侧板与齿轮之间的间隙，具有较好的补偿效果。例如，当泵的工作压力升高时，作用在浮动侧板背面的高压油压力也随之增大，使得浮动侧板与齿轮端面之间的贴合更加紧密，有效减小了因压力升高而可能导致的泄漏增加。为了进一步优化浮动侧板的间隙补偿效果，还可以对其结构进行改进。如在浮动侧板上开设特殊形状的油槽，这些油槽的设计能够改善油液在侧板与齿轮端面之间的分布，减少局部压力不均匀的情况，从而更好地平衡侧板所受的力，提高补偿的稳定性。通过合理设计油槽的形状、尺寸和位置，可以使油液在侧板与齿轮端面之间形成均匀的油膜，既保证了良好的润滑，又能有效减小泄漏。除了浮动侧板结构，还可以采用弹性元件来辅助间隙补偿。在侧板与端盖之间安装弹簧或橡胶等弹性元件，利用弹性元件的弹力使侧板始终保持与齿轮端面的接触，进一步减小间隙。这种方式与高压油作用相结合，能够在不同工况下更好地实现间隙补偿，提高泵的容积效率。在启动阶段，弹性元件的弹力可以使侧板迅速与齿轮端面贴合，减少初始泄漏；在工作过程中，弹性元件可以根据压力变化和侧板的微小位移进行自适应调整，增强补偿的可靠性。压力补偿结构对于高压内啮合齿轮泵的稳定运行至关重要。一种常见的压力补偿结构是在泵体或端盖上设置压力平衡槽。压力平衡槽通过特定的通道与吸油腔和压油腔相连，其作用是在齿轮泵工作时，使齿轮两侧的压力趋于平衡，从而减小不平衡径向力。具体来说，当齿轮在旋转过程中，齿顶与泵体内壁之间的间隙会导致油液泄漏，在压油区，泄漏的油液会形成较高的压力，而在吸油区压力较低，这就产生了不平衡径向力。压力平衡槽通过引入吸油腔和压油腔的压力，在齿轮的不同位置形成相应的压力分布，抵消一部分不平衡径向力。为了更好地实现压力补偿，还可以对压力平衡槽的结构参数进行优化。调整平衡槽的深度、宽度和位置，以适应不同的工作压力和流量要求。较深的平衡槽可以容纳更多的油液，从而在压力变化时提供更稳定的压力补偿；较宽的平衡槽则可以增加油液的流通面积，降低压力损失。通过合理选择平衡槽的位置，使其能够在最需要补偿的区域发挥作用，提高压力补偿的效果。在设计压力平衡槽时，还需要考虑其与泵内其他结构的兼容性，避免对泵的整体性能产生负面影响。在高压内啮合齿轮泵中，还可以采用静压支撑技术来实现压力补偿和减小不平衡径向力。静压支撑结构通常由静压支撑槽和节流器等组成。静压支撑槽开设在齿轮的端面上，通过节流器与高压油源相连。当高压油通过节流器进入静压支撑槽时，会在槽内形成高压油膜，使齿轮在油膜的支撑下悬浮，从而减小齿轮与侧板之间的摩擦和磨损，同时也能有效平衡不平衡径向力。静压支撑技术具有高精度、高承载能力和低摩擦的优点，能够显著提高高压内啮合齿轮泵的工作性能和可靠性。在高速、高压的工况下，静压支撑技术能够有效降低齿轮与侧板之间的磨损，延长泵的使用寿命，同时减少因不平衡径向力引起的振动和噪声，提高泵的稳定性。1.4研究内容与方法1.4.1研究方法本研究综合运用CFD数值计算、理论分析和试验研究等多种方法，深入探究高压内啮合齿轮泵的特性和不平衡径向力问题。在CFD数值计算方面，基于计算流体力学软件FLUENT，利用2.5D动网格的方法建立高压内啮合齿轮泵的三维数值计算模型。通过设置合理的边界条件和求解参数，模拟泵内部的流场特性，分析不同工况下泵的压力分布、速度分布、流量脉动等性能参数。采用粘性壁面的方法模拟两齿轮之间啮合点的存在，以提高内啮合齿轮泵高压工况下的计算精度，并深入分析啮合点处粘度的设置对高压内啮合齿轮泵计算效率和计算精度的影响。在理论分析方面，针对内啮合齿轮泵的特殊结构，推导不平衡径向力和静压支撑承载力的计算公式。基于内啮合齿轮泵具体结构的精确化分析，建立一种不平衡径向力求解方法，深入研究内齿圈所受不平衡径向力的变化特性。分析静压支撑不同结构参数对其承载力、泄漏量的影响，为静压支撑的结构设计提供理论指导。以泄漏量满足要求，尽可能补偿不平衡径向力为目标，确定内啮合齿轮泵的静压支撑结构形式。在试验研究方面，搭建内啮合齿轮泵泵组的测试试验台，利用该试验台对高压内啮合齿轮泵开展相应的试验研究。通过试验，测量不同设计参数（间隙尺寸、月牙形隔板端面槽位置、内齿圈引流孔的大小和泵体等）和不同操作工况（介质温度、工作压力、转速等）对内啮合齿轮泵水力性能的影响，主要包括对内啮合齿轮泵的效率、压力脉动和自吸性能的影响。将试验结果与CFD数值计算和理论分析结果进行对比验证，以确保研究结果的可靠性和准确性。1.4.2研究内容本研究的具体内容涵盖多个关键方面，旨在全面深入地探究高压内啮合齿轮泵的性能和不平衡径向力问题。首先，建立高精度的数值模型。基于计算流体力学软件FLUENT，运用2.5D动网格方法构建高压内啮合齿轮泵的三维数值计算模型。此模型能够精准模拟泵内部的复杂流场，为后续的分析提供坚实基础。采用粘性壁面方法模拟两齿轮间的啮合点，深入研究啮合点处粘度设置对计算精度和效率的影响，以不断优化数值计算模型，提高计算结果的准确性。其次，深入研究径向力。基于MIXTURE模型和Zwart-Gerber-Belamri空化模型，在单相计算的基础上建立适用于高压内啮合齿轮泵的空化两相流动和不平衡径向力求解的数值计算模型。借助该模型，详细研究高压内啮合齿轮泵内部压力沿周向的分布规律以及不平衡径向力的特性。深入分析不同工作压力下内啮合齿轮泵中齿轮轴与内齿圈所承受的不平衡径向力大小，揭示其变化规律，为后续的结构优化提供理论依据。再者，开展静压支撑分析。针对内齿圈与泵体这一关键摩擦副，基于内啮合齿轮泵具体结构的精确化分析，建立不平衡径向力求解方法，推导内齿圈所受不平衡径向力的变化特性。同时，针对内啮合齿轮泵的特殊结构，推导静压支撑承载力的计算公式，系统分析静压支撑不同结构参数对其承载力、泄漏量的影响。以泄漏量满足要求，尽可能补偿不平衡径向力为目标，确定内啮合齿轮泵的静压支撑结构形式，并通过相关试验对结果进行验证，确保静压支撑结构的有效性和可靠性。最后，进行试验验证。搭建内啮合齿轮泵泵组的测试试验台，利用该试验台对高压内啮合齿轮泵进行全面的试验研究。详细分析不同设计参数（间隙尺寸、月牙形隔板端面槽位置、内齿圈引流孔的大小和泵体等）和不同操作工况（介质温度、工作压力、转速等）对内啮合齿轮泵水力性能的影响，主要包括对内啮合齿轮泵的效率、压力脉动和自吸性能的影响。将试验结果与理论研究和数值计算结果进行对比，验证数值计算方法和理论研究方法的可靠性，为高压内啮合齿轮泵的设计和优化提供实际数据支持。二、基于2.5D动网格技术和考虑啮合点的高压内啮合齿轮泵数值模型2.1几何模型和网格2.1.1几何模型高压内啮合齿轮泵的三维几何结构较为复杂，主要由内齿轮、外齿轮、月牙形隔板、泵体和端盖等部件组成。在建模过程中，为了在保证计算精度的前提下提高计算效率，需要对一些关键部件进行合理的简化与处理。内齿轮和外齿轮是高压内啮合齿轮泵的核心部件，其齿形参数对泵的性能有着重要影响。在建模时，精确绘制内、外齿轮的渐开线齿形，确保齿形参数的准确性。考虑到实际加工过程中的公差和表面粗糙度等因素可能会对泵的性能产生一定影响，但由于这些因素的影响相对较小，且在数值模拟中难以精确考虑，因此在本次建模中对这些因素进行了适当简化，忽略了微小的齿形偏差和表面粗糙度。月牙形隔板位于内、外齿轮之间，其作用是将吸油腔和压油腔分隔开，保证泵的正常工作。在建模过程中，对月牙形隔板的形状和尺寸进行了精确描述，确保其能够准确地实现分隔功能。同时，考虑到月牙形隔板与内、外齿轮之间的间隙对泵的泄漏量有一定影响，在建模时对该间隙进行了合理设置，以更准确地模拟泵的实际工作情况。泵体和端盖主要起支撑和密封作用。在建模时，简化了泵体和端盖上一些与内部流场无关的结构，如安装孔、散热筋等，以减少模型的复杂度和计算量。但对于与内部流场直接相关的部分，如泵体的流道、端盖与齿轮的配合面等，进行了详细的建模，以确保能够准确模拟内部流场的特性。具体的几何参数如下：内齿轮齿数为z_1，模数为m，压力角为\alpha，齿顶圆直径为d_{a1}，齿根圆直径为d_{f1}；外齿轮齿数为z_2，模数与内齿轮相同为m，压力角也为\alpha，齿顶圆直径为d_{a2}，齿根圆直径为d_{f2}。内、外齿轮的中心距为a，其计算公式为a=\frac{m(z_2-z_1)}{2}。月牙形隔板的厚度为b，其与内、外齿轮之间的间隙为\delta。泵体的内腔直径为D，长度为L。这些几何参数是根据实际设计要求和相关标准确定的，它们的合理取值对于保证高压内啮合齿轮泵的性能至关重要。通过精确的几何建模和参数设置，可以为后续的数值模拟提供准确的模型基础，从而更深入地研究泵内部的流场特性和不平衡径向力分布规律。2.1.2网格及2.5D动网格技术在对高压内啮合齿轮泵进行数值模拟时，网格划分是一个关键步骤。采用合适的网格划分方法能够提高计算精度和计算效率。本次研究中，使用专业的网格划分软件ICEMCFD对齿轮泵的几何模型进行网格划分。对于内齿轮、外齿轮和月牙形隔板等复杂部件，采用非结构化四面体网格进行划分，以更好地适应其复杂的几何形状，提高网格质量。在齿面等关键部位，通过局部加密网格的方式，提高网格的分辨率，确保能够准确捕捉到流场的细节变化。对于泵体和端盖等相对规则的部件，则采用结构化六面体网格进行划分，这种网格具有较高的计算效率和较好的收敛性。在划分网格时，还需要考虑网格的尺寸和质量。通过多次试验和对比分析，确定了合适的网格尺寸，以保证在满足计算精度要求的前提下，尽量减少计算量。同时，对网格质量进行严格检查，确保网格的正交性、纵横比等指标符合要求，避免因网格质量问题导致计算结果不准确或计算过程不稳定。2.5D动网格技术是一种在模拟旋转机械内部流场时常用的技术，它结合了二维和三维动网格的优点，能够有效地处理齿轮等部件的旋转运动。其原理是将三维模型在某一方向上进行切片，将三维问题简化为二维问题进行处理，同时通过在切片方向上设置合适的边界条件，来考虑三维效应。在模拟齿轮泵内部流场时，2.5D动网格技术具有以下优势：一方面，相较于传统的三维动网格技术，它可以大大减少计算量，提高计算效率。由于将三维问题简化为二维问题处理，减少了网格数量和计算自由度，使得计算时间大幅缩短，这对于处理复杂的齿轮泵内部流场问题尤为重要，可以在有限的计算资源下实现更快速的模拟分析。另一方面，该技术能够较为准确地模拟齿轮的旋转运动和内部流场的动态变化。通过合理设置动网格参数和边界条件，能够精确捕捉到齿轮旋转过程中流场的瞬态特性，为研究泵的性能提供更详细的信息。实施2.5D动网格技术的步骤如下：首先，在ICEMCFD中对齿轮泵的几何模型进行网格划分，生成二维平面网格。在划分网格时，要确保平面网格能够准确反映齿轮泵内部流场的主要特征。然后，将二维平面网格沿指定方向进行拉伸，生成三维网格。在拉伸过程中，需要设置合适的拉伸参数，如拉伸长度、拉伸层数等，以保证生成的三维网格质量良好。接着，在FLUENT软件中设置动网格参数，包括动网格类型、运动方式、运动速度等。对于齿轮的旋转运动，设置为绕轴旋转，指定旋转轴和旋转速度。最后，定义边界条件，包括进口边界条件、出口边界条件、壁面边界条件等。进口边界条件通常设置为速度入口或质量流量入口，根据实际工况确定进口速度或质量流量；出口边界条件设置为压力出口，指定出口压力；壁面边界条件设置为无滑移边界条件，模拟实际的壁面情况。通过以上步骤，即可实现基于2.5D动网格技术的高压内啮合齿轮泵内部流场的数值模拟。2.2边界条件与计算策略在高压内啮合齿轮泵的数值模拟中，边界条件的合理设置对于准确模拟泵内部的流场特性至关重要。进口边界条件设置为速度入口。根据实际工况，确定进口油液的流速v_{in}。例如，在某一特定工况下，已知泵的流量为Q，进口管道的横截面积为A_{in}，则进口流速v_{in}=\frac{Q}{A_{in}}。在设置速度入口边界条件时，还需考虑油液的湍流特性，指定合适的湍流强度和水力直径，以准确模拟进口油液的流动状态。出口边界条件设置为压力出口。根据泵的工作压力要求，确定出口压力p_{out}。通常，出口压力是已知的设计参数，在模拟过程中，将其设置为恒定值，以模拟泵在实际工作中向系统输出高压油液的情况。同时，在压力出口边界条件中，也需要考虑出口油液的湍流特性，确保边界条件的准确性。壁面边界条件设置为无滑移边界条件。对于内齿轮、外齿轮、月牙形隔板、泵体和端盖等部件的壁面，认为油液与壁面之间没有相对滑动，即油液在壁面上的速度为零。这种设置符合实际物理情况，能够准确模拟壁面对油液流动的约束作用。考虑到壁面粗糙度对油液流动的影响，在模拟中可以通过设置壁面粗糙度参数来体现。对于一些表面加工精度较高的部件，壁面粗糙度较小，对油液流动的影响相对较小；而对于表面加工精度较低的部件，壁面粗糙度较大，可能会导致油液在壁面附近产生较大的速度梯度和能量损失。在计算策略方面，求解器的选择至关重要。选用基于压力基的求解器，如SIMPLE算法（Semi-ImplicitMethodforPressure-LinkedEquations），该算法在处理不可压缩流体流动问题时具有良好的稳定性和收敛性。在求解过程中，采用分离式求解方法，将压力和速度的求解过程分开，通过迭代的方式逐步逼近真实解。时间步长的设定直接影响计算效率和计算精度。通过多次试验和分析，根据齿轮的转速和旋转周期来确定合适的时间步长。若齿轮的转速为n（单位：转/分钟），则旋转周期T=\frac{60}{n}（单位：秒）。在模拟中，将时间步长\Deltat设置为旋转周期的一小部分，例如\Deltat=\frac{T}{N}，其中N为一个适当的整数，如N=100，这样可以在保证计算精度的前提下，有效控制计算时间。同时，在模拟过程中，还需要对时间步长进行敏感性分析，观察不同时间步长下计算结果的变化情况，以确保时间步长的设置是合理且准确的。若时间步长过大，可能会导致计算结果的失真，无法准确捕捉到流场的瞬态变化；而时间步长过小，则会增加计算量，延长计算时间。为了提高计算效率，还可以采用并行计算技术。利用多核心处理器或集群计算机，将计算任务分配到多个计算节点上同时进行计算，从而大大缩短计算时间。在并行计算过程中，需要合理划分计算区域，确保各个计算节点之间的数据通信和同步高效进行，以充分发挥并行计算的优势。2.3基于粘性壁面方法的啮合点的数值模拟在高压内啮合齿轮泵的工作过程中，两齿轮的啮合点是一个关键部位，其对泵内部的流动特性和性能有着重要影响。利用粘性壁面方法可以有效地模拟两齿轮之间啮合点的存在，从而更准确地分析泵的工作情况。粘性壁面方法的原理基于流体的粘性特性。在实际的齿轮泵中，两齿轮在啮合点处，由于齿面的相互作用和流体的粘性，会形成一个相对复杂的流动区域。粘性壁面方法通过在啮合点处设置适当的粘度条件，来模拟这种复杂的流动情况。具体来说，在数值模拟中，将啮合点处的流体视为具有较高粘度的流体层。这是因为在实际的啮合点处，齿面之间的间隙非常小，流体在其中流动时受到的粘性阻力较大，等效于流体的粘度增加。通过增大啮合点处的流体粘度，可以更好地模拟实际情况下流体在啮合点处的流动行为，如流速降低、泄漏减少等现象。啮合点处粘度设置对计算效率和精度有着显著影响。从计算精度方面来看，合理的粘度设置能够更准确地模拟啮合点处的流动特性，从而提高计算结果的准确性。当啮合点处的粘度设置过低时，无法准确反映实际情况下流体在啮合点处受到的粘性阻力，导致模拟的泄漏量偏大，流量和压力分布等计算结果与实际情况存在较大偏差。而当粘度设置过高时，虽然能够更有效地抑制啮合点处的泄漏，但可能会导致计算结果过于保守，与实际情况也不完全相符。因此，需要通过大量的数值试验和对比分析，结合实际泵的工作情况，确定合适的粘度值，以达到最佳的计算精度。在计算效率方面，啮合点处粘度的设置也会产生影响。较高的粘度设置会增加流体的粘性力，使得流体的流动变得更加缓慢和复杂，从而导致计算时间增加，计算效率降低。这是因为在数值计算过程中，需要求解包含粘性项的流体力学方程，粘度的增加会使方程的求解变得更加困难，需要更多的迭代次数和计算资源来达到收敛。相反，较低的粘度设置虽然可以提高计算效率，但会牺牲计算精度。因此，在实际模拟中，需要在计算精度和计算效率之间进行权衡，找到一个合适的粘度设置，既能保证计算结果的准确性，又能在合理的时间内完成计算。可以通过逐步调整粘度值，观察计算结果的变化趋势，结合实际工程需求，确定出最优的粘度设置方案。2.4结果与分析2.4.1考虑啮合点对齿轮泵流量的影响在高压内啮合齿轮泵的性能研究中，考虑啮合点对齿轮泵流量的影响是一个重要方面。通过数值模拟，对比考虑啮合点与不考虑啮合点时齿轮泵的流量计算结果，能够深入了解啮合点对流量特性的作用规律。在不同工况下，考虑啮合点时齿轮泵的流量表现出与不考虑啮合点时不同的特性。在工作压力从2.88MPa增大到25MPa的过程中，当转过一个齿时，考虑啮合点的存在使计算流量显著提高，尤其是在高压工况（16.67MPa和25MPa）下，这种升高幅度更加明显。这是因为在实际工作中，两齿轮在啮合点处，由于齿面的相互作用和流体的粘性，形成了一个相对复杂的流动区域。利用粘性壁面方法模拟啮合点时，增大了啮合点处的流体粘度，使其更贴近实际啮合点的作用，阻碍了两齿轮之间啮合点处的间隙泄漏通道，导致泄漏速度降低，泄漏量减小，从而使流量有所升高。以某一具体工况为例，在工作压力为16.67MPa时，不考虑啮合点时计算得到的流量为Q_1，而考虑啮合点时流量增加到Q_2，流量升高的幅度为\frac{Q_2-Q_1}{Q_1}\times100\%，经计算，该幅度达到了[X]%，充分说明了在高压工况下考虑啮合点对流量的显著影响。随着啮合点处粘度的设置增大，计算流量也随之增大。这是因为粘度的增加进一步抑制了啮合点处的泄漏，使得更多的油液能够被有效输送，从而提高了流量。然而，当啮合点处设置的粘度继续增大时，计算流量的升高幅值有所减小。这是由于当粘度增大到一定程度后，泄漏量已经被大幅抑制，继续增加粘度对泄漏量的影响逐渐减小，因此流量的升高幅度也相应减小。在齿轮转过一个齿的过程中，考虑啮合点的计算还使流量的脉动情况发生变化。不考虑啮合点时，内啮合齿轮泵的流量在齿轮转过一个齿的过程中只有一个很平缓的峰值；而在考虑啮合点时，流量在转过一个齿的过程中会出现3个较大的峰值，其中一个峰值持续时间较长，后两个峰值持续时间较短，但是峰值在转动过程中有所减小。不同时间对应的流量相对应，啮合点处的泄漏速度较大时，泄漏量较大，流量较小，故对应相同时间时为流量峰谷值；泄漏速度最小时，对应的流量峰值。这是因为不考虑啮合点时，齿轮与齿轮间存在间隙，通过间隙的流动比较均匀，泄漏量较大，受齿轮与齿轮间相对位置的影响不大。而在考虑啮合点时，由于啮合点的考虑使齿轮之间的泄漏量开始受齿轮相对位置的影响，内啮合齿轮泵的内部泄漏量在不同齿轮转角时变化较大。此外，高压工况（16.67MPa和25MPa）时的排量峰峰值要远大于中低压（8.33MPa和2.88MPa）时的排量峰峰值，即工作压力增大时伴随着流量脉动的增大。2.4.2考虑啮合点的内啮合齿轮泵内的压力分布通过数值模拟得到考虑啮合点的内啮合齿轮泵内的压力云图，能够直观地展示泵内的压力分布情况。在不同工况下，分析压力沿周向和径向的分布规律，对于深入理解泵的工作原理和性能具有重要意义。从压力云图可以看出，在吸油区，压力较低，接近大气压；在压油区，压力较高，达到工作压力。在齿轮啮合点附近，压力分布较为复杂。由于啮合点处的间隙较小，流体流动受到较大的阻力，导致压力升高。随着齿轮的旋转，啮合点的位置不断变化，压力分布也随之动态变化。在周向方向上，压力从吸油区到压油区逐渐升高。在吸油区，油液在大气压的作用下被吸入齿轮的齿槽中，此时压力较低且分布相对均匀。随着齿轮的旋转，齿槽中的油液被逐渐带入压油区，在压油区，由于齿槽容积的减小和油液的压缩，压力迅速升高。在靠近压油口的区域，压力达到最大值，然后逐渐降低至出口压力。在不同的工作压力下，周向压力分布的趋势基本相同，但压力的具体数值会随着工作压力的升高而增大。当工作压力从10MPa增加到20MPa时，压油区的最高压力也相应从[X1]MPa增加到[X2]MPa。在径向方向上，压力分布也呈现出一定的规律。从齿顶到齿根，压力逐渐降低。这是因为在齿顶处，油液直接与高压区相连，受到的压力较大；而在齿根处，油液受到齿根的阻挡和密封作用，压力相对较小。在靠近泵体壁面的区域，由于油液与壁面之间的摩擦和粘性作用，压力也会有所降低。此外，在啮合点处，由于齿面的相互作用和流体的粘性，径向压力分布也会发生局部变化。啮合点对压力分布起着重要的作用。啮合点处的高压区域会影响周围油液的流动和压力分布。在啮合点附近，由于间隙泄漏受到抑制，油液的流动速度减慢，压力升高，形成局部高压区。这个局部高压区会对齿轮的受力和磨损产生影响，同时也会影响泵的容积效率和流量特性。当啮合点处的粘度设置发生变化时，压力分布也会相应改变。增大啮合点处的粘度，会使局部高压区的压力进一步升高，范围扩大；减小粘度，则会使局部高压区的压力降低，范围缩小。2.4.3模拟啮合点的位置准确性为了验证模拟啮合点位置的准确性，将模拟结果与试验数据或理论值进行对比分析。通过对比，可以评估模型的可靠性，确保数值模拟结果的有效性。在与试验数据对比时，搭建了内啮合齿轮泵的试验台，测量不同工况下泵内的压力分布和流量特性。将试验得到的压力分布和流量数据与数值模拟结果进行详细对比。在某一特定工况下，试验测得的吸油区压力为p_{s-exp}，压油区压力为p_{d-exp}，而数值模拟得到的吸油区压力为p_{s-sim}，压油区压力为p_{d-sim}。通过计算压力相对误差\delta_p=\frac{|p_{exp}-p_{sim}|}{p_{exp}}\times100\%，得到吸油区压力相对误差为\delta_{p-s}，压油区压力相对误差为\delta_{p-d}。经计算，吸油区压力相对误差在[X]%以内，压油区压力相对误差在[X]%以内，表明数值模拟得到的压力分布与试验数据较为吻合。在流量特性方面，试验测得的流量为Q_{exp}，数值模拟得到的流量为Q_{sim}，计算流量相对误差\delta_Q=\frac{|Q_{exp}-Q_{sim}|}{Q_{exp}}\times100\%，结果显示流量相对误差在[X]%以内，说明模拟得到的流量与试验值接近，验证了模拟啮合点位置对流量计算的准确性。与理论值对比时，利用已有的内啮合齿轮泵理论计算公式，计算不同工况下的压力分布和流量特性。将理论计算结果与数值模拟结果进行对比分析。在压力分布方面，理论计算得到的压力分布趋势与数值模拟结果一致，在一些关键位置的压力数值相对误差也在可接受范围内。在流量计算方面，理论计算值与模拟值的相对误差在[X]%以内，进一步验证了模拟啮合点位置的准确性。通过与试验数据和理论值的对比验证，表明所采用的基于粘性壁面方法模拟啮合点的数值模型能够较为准确地预测内啮合齿轮泵的压力分布和流量特性，模拟啮合点的位置具有较高的准确性，为后续深入研究高压内啮合齿轮泵的性能和不平衡径向力问题提供了可靠的模型基础。2.5本章小结本章基于2.5D动网格技术和考虑啮合点的方法，成功建立了高压内啮合齿轮泵的数值模型，为深入研究其内部流场特性和不平衡径向力分布规律奠定了坚实基础。在模型建立过程中，通过对高压内啮合齿轮泵的几何模型进行合理简化与处理，确保在提高计算效率的同时，最大程度保留关键部件对泵性能的影响。利用ICEMCFD进行网格划分，针对不同部件的特点采用了非结构化四面体网格和结构化六面体网格相结合的方式，并在关键部位局部加密网格，有效保证了网格质量和计算精度。引入2.5D动网格技术，该技术将三维问题简化为二维问题处理，通过在切片方向上设置合适边界条件考虑三维效应，显著减少了计算量，同时能较为准确地模拟齿轮旋转运动和内部流场动态变化，为数值模拟提供了高效可靠的方法。边界条件的设置直接影响模拟结果的准确性。进口边界条件设置为速度入口，根据实际工况确定进口油液流速，并考虑湍流特性指定湍流强度和水力直径；出口边界条件设置为压力出口，根据工作压力要求确定出口压力，并同样考虑湍流特性；壁面边界条件设置为无滑移边界条件，同时考虑壁面粗糙度对油液流动的影响。在计算策略上，选用基于压力基的求解器SIMPLE算法，采用分离式求解方法，并通过多次试验和分析，根据齿轮转速和旋转周期确定合适的时间步长，同时进行敏感性分析，确保时间步长设置合理准确。此外，采用并行计算技术提高计算效率，充分利用多核心处理器或集群计算机资源，合理划分计算区域，保证数据通信和同步高效进行。针对两齿轮啮合点这一关键部位，利用粘性壁面方法模拟其存在。通过在啮合点处设置较高粘度的流体层，有效模拟了实际情况下流体在啮合点处受到的粘性阻力，阻碍了间隙泄漏通道，降低了泄漏速度和泄漏量。研究发现，啮合点处粘度设置对计算效率和精度影响显著，合理的粘度设置能提高计算精度，但过高的粘度会增加计算时间、降低计算效率，因此需在两者之间进行权衡，通过大量数值试验和对比分析确定最优粘度设置方案。通过数值模拟，深入分析了考虑啮合点对齿轮泵流量和压力分布的影响。在流量方面，考虑啮合点时齿轮泵流量在不同工况下均有显著提高，尤其是在高压工况下，且随着啮合点处粘度增大，流量增大，但粘度继续增大时流量升高幅值减小。同时，考虑啮合点使流量脉动情况发生变化，出现多个峰值，且高压工况下流量脉动更大。在压力分布方面，通过压力云图直观展示了泵内压力分布情况，吸油区压力低，压油区压力高，齿轮啮合点附近压力分布复杂。在周向方向上，压力从吸油区到压油区逐渐升高；在径向方向上，从齿顶到齿根压力逐渐降低。啮合点处的高压区域影响周围油液流动和压力分布，其粘度设置变化会导致压力分布改变。通过与试验数据和理论值对比，验证了模拟啮合点位置的准确性，为后续研究提供了可靠模型基础。然而，本模型仍存在一定局限性。在简化几何模型时，虽对一些微小因素进行了忽略以提高计算效率，但这些因素在某些特殊工况下可能对泵性能产生不可忽视的影响，后续研究可考虑进一步细化模型以更全面反映实际情况。在模拟啮合点时，尽管粘性壁面方法取得了较好效果，但实际啮合点处的物理现象更为复杂，模型可能无法完全准确模拟所有细节，未来可探索更精确的模拟方法。三、基于空化两相流动的高压内啮合齿轮泵不平衡径向力研究3.1高压内啮合齿轮泵空化流动模型的建立3.1.1空化模型及其控制方程在高压内啮合齿轮泵的工作过程中，空化现象是影响其性能和可靠性的重要因素之一。为了准确模拟泵内的空化流动，采用MIXTURE模型和Zwart-Gerber-Belamri空化模型相结合的方式。MIXTURE模型是一种多相流模型，它基于混合流体的假设，将多相流视为一种具有平均特性的混合流体来处理。在MIXTURE模型中，通过引入体积分数来描述各相在混合流体中的比例，同时考虑了相间的速度滑移和质量传递。其连续性方程为：\frac{\partial\rho_m}{\partialt}+\nabla\cdot(\rho_m\vec{v}_m)=0其中，\rho_m为混合流体的密度，\vec{v}_m为混合流体的速度矢量。动量方程为：\frac{\partial(\rho_m\vec{v}_m)}{\partialt}+\nabla\cdot(\rho_m\vec{v}_m\vec{v}_m)=-\nablap+\nabla\cdot(\mu_m\nabla\vec{v}_m)+\rho_m\vec{g}+\vec{F}这里，p为压力，\mu_m为混合流体的动力粘度，\vec{g}为重力加速度矢量，\vec{F}为其他体积力。Zwart-Gerber-Belamri空化模型是一种基于气泡动力学的空化模型，它考虑了气泡的生成、生长和溃灭过程。该模型假设系统中所有的气泡具有相同的大小，空化率R可用气泡数密度和单个气泡的质量变化率相乘得到。当压力小于饱和蒸汽压时，发生蒸发过程，空化率R的表达式为：R_{evap}=F_{vap}\frac{3\alpha_nuc(1-\alpha_v)}{\Re_B}\sqrt{\frac{2}{3}\frac{p_{sat}-p}{\rho_l}}其中，F_{vap}为蒸发系数，\alpha_nuc为成核位点体积分数，\alpha_v为蒸汽相体积分数，\Re_B为气泡半径，p_{sat}为饱和蒸汽压，p为当地压力，\rho_l为液相密度。当压力大于饱和蒸汽压时，发生冷凝过程，空化率R的表达式为：R_{cond}=F_{cond}\frac{3\alpha_v}{\Re_B}\sqrt{\frac{2}{3}\frac{p-p_{sat}}{\rho_l}}F_{cond}为冷凝系数。在实际应用中，将Zwart-Gerber-Belamri空化模型的空化率R作为源项添加到MIXTURE模型的连续性方程和动量方程中，以考虑空化对流动的影响。3.1.2空化流动的计算策略在进行空化流动的计算时，合理的计算策略对于确保计算结果的准确性和计算过程的稳定性至关重要。初始化过程是计算的起始步骤。在初始化阶段，需要对计算域内的各物理量进行初始赋值。对于压力场，通常将其初始化为一个接近实际工况的压力值，例如将进口压力设置为实际的进口压力，出口压力设置为实际的出口压力，其他区域的压力根据进出口压力进行线性插值或其他合理的估算方法进行初始化。对于速度场，根据泵的工作原理和实际运行情况，假设初始速度分布。在齿轮泵的进口区域，可以根据进口流量和进口截面积计算得到进口速度，并以此为基础对进口附近区域的速度进行初始化；对于其他区域，可以根据流动的连续性和假设的流动模式进行速度的初步赋值。对于蒸汽相体积分数，由于在初始时刻通常假设没有空化发生，因此可以将其初始值设置为一个极小的值，如1\times10^{-6}，以避免计算过程中出现数值问题。收敛准则的设定是判断计算结果是否达到稳定状态的重要依据。在空化流动计算中，常用的收敛准则包括残差收敛准则和物理量收敛准则。残差收敛准则是通过监测计算过程中各控制方程的残差来判断计算是否收敛。对于连续性方程、动量方程和空化模型的输运方程等，设定一个合理的残差阈值，如1\times10^{-6}。当各方程的残差在迭代计算过程中逐渐减小，并小于设定的阈值时，认为该方程的计算达到收敛。物理量收敛准则则是通过监测一些关键物理量的变化来判断计算是否收敛。在空化流动中，可以选择监测泵的进出口流量、压力、蒸汽相体积分数等物理量。当这些物理量在连续的若干次迭代计算中变化非常小，例如进出口流量的变化率小于0.1\%，压力的变化小于设定的压力阈值，蒸汽相体积分数的变化也在可接受的范围内时，认为计算达到物理量收敛。只有当残差收敛准则和物理量收敛准则同时满足时，才认为整个计算过程达到收敛，得到的计算结果是可靠的。在计算过程中，可能会出现一些问题影响计算的收敛性和准确性。例如，由于空化过程中气泡的生成和溃灭会导致流场的剧烈变化，可能会引起数值振荡。为了解决这个问题，可以采用一些数值稳定技术，如增加松弛因子。松弛因子是一种用于控制迭代过程中物理量更新幅度的参数。在空化流动计算中，对于压力、速度、蒸汽相体积分数等物理量的更新，可以适当减小其更新幅度，通过调整松弛因子的值，如将压力松弛因子设置为0.3-0.7，速度松弛因子设置为0.5-0.8，蒸汽相体积分数松弛因子设置为0.5-1，使得计算过程更加稳定，避免数值振荡的发生。此外，还可以采用自适应网格技术，根据流场的变化情况自动调整网格的疏密程度。在空化区域，由于流场变化剧烈，通过加密网格可以更准确地捕捉流场的细节信息，提高计算精度；而在流场变化较小的区域，可以适当粗化网格，以减少计算量，提高计算效率。3.1.3空化模型控制参数的确定空化模型中的控制参数对计算结果的准确性和可靠性有着重要影响，因此需要依据实际工况和经验，合理确定这些参数的值。气泡半径\Re_B是Zwart-Gerber-Belamri空化模型中的一个关键参数，它表示气泡的大小。在实际应用中，气泡半径的取值与液体的性质、流动状态以及空化的发展阶段等因素有关。根据相关研究和经验，对于常见的液压油等液体，气泡半径通常可以取值为1\times10^{-6}m。当液体中含有杂质或表面活性剂等物质时，气泡半径可能会发生变化。如果液体中含有较多的杂质，这些杂质可以作为气泡的核心，促进气泡的生成，使得气泡半径相对较小；而当液体中存在表面活性剂时，表面活性剂会降低液体的表面张力，使得气泡更容易生成和生长，气泡半径可能会相对较大。因此，在确定气泡半径时，需要综合考虑这些因素，并通过实验或数值模拟进行验证和调整。成核位点体积分数\alpha_nuc表示单位体积液体中气泡成核位点的体积比例，它反映了气泡生成的难易程度。成核位点体积分数的取值与液体的纯净度、温度、压力等因素有关。一般来说，对于纯净的液体，成核位点体积分数较小；而当液体中存在微小颗粒、溶解气体或表面缺陷等时，成核位点体积分数会增大。在实际计算中，成核位点体积分数通常取值为5\times10^{-4}。但在不同的工况下，需要对其进行调整。在高温高压的工况下，液体的物理性质会发生变化，气泡的生成和生长机制也会有所不同，此时可能需要适当调整成核位点体积分数的值，以更准确地模拟空化现象。蒸发系数F_{vap}和冷凝系数F_{cond}分别控制着蒸发和冷凝过程的速率。它们的取值会影响空化过程中气泡的生成和溃灭速度，进而影响整个流场的特性。蒸发系数F_{vap}通常取值为50，冷凝系数F_{cond}通常取值为0.01。然而，这些值并不是固定不变的，在不同的液体和流动条件下，需要对其进行敏感性分析。当液体的粘性较大时，气泡的运动和变形受到的阻力较大，蒸发和冷凝过程可能会受到影响，此时可以适当调整蒸发系数和冷凝系数的值，以更好地反映实际情况。通过改变蒸发系数和冷凝系数的值，观察计算结果中蒸汽相体积分数、压力分布、速度分布等物理量的变化，从而确定最适合当前工况的参数值。对这些控制参数进行敏感性分析，可以深入了解它们对计算结果的影响程度。通过逐步改变某个控制参数的值，保持其他参数不变，进行多次数值计算，对比不同参数值下的计算结果，如蒸汽相体积分数的分布、空化区域的大小、泵的性能参数等。分析结果表明，气泡半径和成核位点体积分数对空化起始和发展的位置有较大影响，较小的气泡半径和较大的成核位点体积分数会使空化更容易发生，空化区域也会相应扩大；蒸发系数和冷凝系数则主要影响空化的速率和强度，较大的蒸发系数会加快气泡的生成速度，而较大的冷凝系数会加快气泡的溃灭速度。在实际应用中，需要根据具体的研究目的和实际工况，综合考虑这些参数的影响，合理确定其取值，以提高空化流动计算的准确性和可靠性。3.2高压内啮合齿轮泵内的空化两相流动计算3.2.1考虑空化的两相流动计算对流量的影响在高压内啮合齿轮泵的工作过程中，空化现象的发生会对泵的流量产生显著影响。通过数值模拟，对比单相流动计算和考虑空化的两相流动计算时的流量结果，可以清晰地观察到这种影响的程度和变化趋势。在不同工况下，考虑空化的两相流动计算得到的流量与单相流动计算结果存在明显差异。随着工作压力的升高，空化现象逐渐加剧，流量的变化也更为显著。当工作压力从较低值逐渐增加时，在单相流动计算中，流量通常呈现出相对稳定的变化趋势，主要受泵的转速、齿轮参数等因素影响。然而，在考虑空化的两相流动计算中，随着压力升高，液体中的气泡开始大量生成和发展，占据了一定的空间，导致实际参与流动的液体体积减少，从而使流量逐渐降低。以某一具体工况为例，在工作压力为10MPa时，单相流动计算得到的流量为Q_{s-10}，而考虑空化的两相流动计算得到的流量为Q_{t-10}，流量降低的幅度为\frac{Q_{s-10}-Q_{t-10}}{Q_{s-10}}\times100\%，经计算，该幅度达到了[X]%。当工作压力进一步升高到20MPa时，单相流动计算流量为Q_{s-20}，考虑空化的两相流动计算流量为Q_{t-20}，此时流量降低幅度增大到[X]%，充分表明了工作压力升高时空化对流量影响的加剧。从变化趋势来看，随着工作压力的持续增加，考虑空化的流量降低幅度逐渐增大，且在高压工况下，流量的波动也更为明显。这是因为在高压下，空化气泡的生成和溃灭过程更加剧烈，对液体的流动产生了更大的干扰，导致流量不稳定。空化还可能导致泵内部流道的局部堵塞，进一步阻碍液体的正常流动，使得流量下降更为显著。在转速方面，随着转速的提高，泵内液体的流速增加，压力变化更加迅速，这也会促使空化现象提前发生和加剧，从而导致流量下降更快。当转速从1000r/min增加到1500r/min时，在相同工作压力下，考虑空化的流量降低幅度从[X1]%增大到[X2]%。3.2.2考虑空化两相流动时的周向压力分布通过数值模拟得到考虑空化两相流动时高压内啮合齿轮泵的周向压力分布曲线，能够直观地了解泵内压力在不同区域的变化规律，以及空化与压力分布之间的紧密关联。在吸油区，压力较低，接近大气压。随着齿轮的旋转，液体被逐渐吸入齿槽中。在这个区域，空化现象相对较少，因为压力高于液体的饱和蒸汽压。随着齿槽向压油区移动，压力逐渐升高。在空化条件下，当压力降低到液体的饱和蒸汽压以下时，空化气泡开始大量生成。这些气泡在压力升高的过程中，会逐渐溃灭，从而对压力分布产生影响。在靠近空化区域的地方，压力分布会出现明显的波动。由于气泡的生成和溃灭，导致局部压力瞬间升高或降低，使得压力曲线呈现出不规则的形状。在压油区，压力较高，达到工作压力。此时，空化气泡的溃灭更加剧烈，因为压力远高于饱和蒸汽压。空化气泡的溃灭会产生局部的高压冲击，对齿面和泵体造成损害，同时也会影响压力分布的均匀性。在压油区的某些位置，由于空化气泡的溃灭，会出现压力峰值，这些峰值可能会导致齿面的疲劳磨损和腐蚀。空化与压力分布之间存在着相互影响的关系。一方面，压力分布决定了空化的发生位置和程度。在压力较低的区域，空化更容易发生；而在压力较高的区域，空化气泡更容易溃灭。另一方面，空化的发生和发展也会改变压力分布。空化气泡的生成和溃灭会导致局部压力的变化，从而影响整个泵内的压力分布。当空化现象严重时，压力分布会变得更加不均匀，这不仅会影响泵的性能，还会降低泵的可靠性和使用寿命。通过分析周向压力分布曲线，可以深入了解空化对泵内压力分布的影响，为优化泵的设计和运行提供依据。3.3内啮合齿轮泵的不平衡径向力在高压内啮合齿轮泵的运行过程中，不平衡径向力是影响其性能和可靠性的关键因素之一。通过数值模拟，详细计算不同工作压力下齿轮轴与内齿圈所承受的不平衡径向力，深入分析其大小、方向及分布规律，对于优化泵的设计和提高其工作性能具有重要意义。在工作压力从2.88MPa增大到25MPa的过程中，齿轮轴所承受的不平衡径向力呈现出明显的变化规律。随着工作压力的升高，不平衡径向力的大小逐渐增大。当工作压力为2.88MPa时，齿轮轴所受不平衡径向力为F_{r-2.88}，而当工作压力升高到25MPa时，不平衡径向力增大到F_{r-25}，增长幅度达到了[X]%。这是因为随着工作压力的增加，齿轮两侧的压力差增大，导致作用在齿轮轴上的径向合力增大。不平衡径向力的方向也会随着齿轮的旋转而发生变化。在一个旋转周期内，不平衡径向力的方向始终指向齿顶与泵体内壁间隙较小的一侧。这是由于在该侧，油液的泄漏量相对较小，压力较高，从而产生了较大的径向力。在齿轮旋转过程中，齿顶与泵体内壁的间隙不断变化，导致不平衡径向力的方向也随之动态改变。从分布规律来看，不平衡径向力在齿轮的齿顶和齿根处相对较大，而在齿面中间部分相对较小。这是因为在齿顶和齿根处，油液的流动受到的阻碍较大，压力变化较为剧烈，从而产生了较大的径向力。在齿面中间部分，油液的流动相对较为平稳，压力变化较小，因此不平衡径向力也较小。内齿圈所承受的不平衡径向力同样随着工作压力的升高而增大。当工作压力从2.88MPa增加到25MPa时，内齿圈所受不平衡径向力从F_{i-2.88}增大到F_{i-25}，增长幅度为[X]%。内齿圈所受不平衡径向力的方向与齿轮轴所受不平衡径向力的方向相反，这是由于内齿圈与齿轮轴之间的相互作用力关系所决定的。在内齿圈的周向上，不平衡径向力的分布也存在一定的规律。在与齿轮啮合的区域，不平衡径向力相对较大，这是因为在该区域，齿面之间的相互作用力和油液的压力变化较为复杂，导致径向力增大。而在远离啮合区域的部分，不平衡径向力相对较小。在径向方向上，内齿圈所受不平衡径向力从内表面到外表面逐渐减小，这是因为内表面直接受到油液压力的作用，而外表面受到的影响相对较小。不平衡径向力的大小、方向及分布规律对高压内啮合齿轮泵的性能有着重要影响。过大的不平衡径向力会加速轴承的磨损，降低轴承的使用寿命，增加设备的维护成本。不平衡径向力还可能导致泵的振动和噪声增大，影响泵的工作稳定性和可靠性。因此，在高压内啮合齿轮泵的设计和优化过程中，需要充分考虑不平衡径向力的影响，采取有效的措施来减小不平衡径向力，如优化齿轮结构、采用压力平衡槽、设计合理的静压支撑结构等，以提高泵的性能和可靠性。3.4静压支撑槽的影响3.4.1静压支撑槽尺寸和位置静压支撑槽在高压内啮合齿轮泵中起着至关重要的作用，其尺寸和位置的合理设计对于泵的性能优化具有关键意义。静压支撑槽的关键尺寸参数主要包括槽的宽度w、深度h以及槽的长度l。槽的宽度决定了油液在槽内的流通面积，对静压支撑力的大小和油液的泄漏量有重要影响。较宽的槽能够提供更大的流通面积，使更多的油液参与到静压支撑过程中，从而增加静压支撑力，但同时也可能导致泄漏量增加。槽的深度则影响着油液在槽内的储存能力和压力分布。较深的槽可以储存更多的油液，在一定程度上稳定静压支撑力，但过深的槽可能会削弱泵体的结构强度。槽的长度与齿轮的啮合区域相关，合适的长度能够确保在齿轮的啮合过程中，静压支撑槽始终发挥有效的支撑作用。静压支撑槽在泵体中的位置坐标可通过建立合适的坐标系来确定。以泵体的中心为原点，沿齿轮的旋转方向为x轴，垂直于旋转方向为y轴。静压支撑槽的中心位置可表示为(x_0,y_0)，其中x_0决定了槽在周向方向上的位置，y_0决定了槽在径向方向上的位置。在周向方向上，槽的位置需要根据齿轮的啮合情况和压力分布来确定，通常设置在压力较高且不平衡径向力较大的区域，以有效补偿不平衡径向力。在径向方向上，槽的位置要考虑与齿轮的间隙以及泵体的结构强度，确保槽的设置不会影响齿轮的正常运转和泵体的可靠性。为了深入研究不同尺寸和位置的组合方案对泵性能的影响，设计了一系列对比方案。在尺寸方面，设置不同的槽宽，如w_1=2mm、w_2=3mm、w_3=4mm；不同的槽深，如h_1=1mm、h_2=1.5mm、h_3=2mm；不同的槽长，如l_1=10mm、l_2=15mm、l_3=20mm。在位置方面，设置不同的周向位置，如x_{01}=30^{\circ}、x_{02}=60^{\circ}、x_{03}=90^{\circ}；不同的径向位置，如y_{01}=5mm、y_{02}=6mm、y_{03}=7mm。通过组合这些不同的尺寸和位置参数，形成多个对比方案，如方案一（w_1、h_1、l_1、x_{01}、y_{01}）、方案二（w_2、h_2、l_2、x_{02}、y_{02}）等。利用数值模拟方法，对每个方案进行计算分析，对比不同方案下泵的性能指标，从而确定最优的静压支撑槽尺寸和位置组合。3.4.2静压支撑槽对内啮合齿轮泵性能的影响通过数值计算，详细对比不同静压支撑槽参数下泵的不平衡径向力、泄漏量、效率等性能指标，能够深入了解静压支撑槽对高压内啮合齿轮泵性能的影响规律。在不平衡径向力方面，不同的静压支撑槽参数对其补偿效果差异显著。当静压支撑槽的宽度增加时，不平衡径向力呈现出先减小后增大的趋势。在一定范围内，较宽的槽能够提供更大的静压支撑力，有效抵消不平衡径向力，使不平衡径向力减小。但当槽宽超过某一临界值时，由于泄漏量的增加，导致静压支撑力的有效作用减弱，不平衡径向力反而增大。槽的深度增加时，不平衡径向力也会发生变化。较深的槽在一定程度上可以增强静压支撑力，减小不平衡径向力，但如果深度过大，会导致油液在槽内的流动阻力增加，影响静压支撑效果，使不平衡径向力增大。在位置方面，当静压支撑槽设置在不平衡径向力较大的区域时，能够更有效地补偿不平衡径向力。将槽设置在靠近压油区的周向位置，由于该区域压力较高，不平衡径向力较大，静压支撑槽能够提供更大的支撑力来平衡不平衡径向力，从而显著降低不平衡径向力的大小。泄漏量是衡量高压内啮合齿轮泵性能的重要指标之一，静压支撑槽的参数对泄漏量有重要影响。随着槽宽的增大，泄漏量明显增加。这是因为较宽的槽提供了更大的泄漏通道，使得油液更容易从槽中泄漏出去。槽深的增加也会导致泄漏量增加，较深的槽虽然能够储存更多的油液，但同时也增加了油液泄漏的可能性。槽的长度增加时，泄漏量也会有所增加，因为更长的槽意味着更大的泄漏面积。在位置方面，当静压支撑槽靠近齿轮的啮合区域时，泄漏量会相对较大，因为该区域的压力变化较为复杂，油液更容易通过槽泄漏。效率是高压内啮合齿轮泵综合性能的体现，静压支撑槽参数对效率的影响是不平衡径向力和泄漏量共同作用的结果。当静压支撑槽参数能够有效减小不平衡径向力且泄漏量控制在合理范围内时，泵的效率会得到提高。通过优化静压支撑槽的尺寸和位置，使得不平衡径向力得到较好的补偿，同时泄漏量最小化，此时泵的效率达到最高。若静压支撑槽参数不合理，导致不平衡径向力较大或泄漏量过大，都会使泵的效率降低。当不平衡径向力过大时，会增加轴承的摩擦损失和泵的振动，消耗更多的能量；当泄漏量过大时，会导致实际输出流量减小，能量损失增加，从而降低泵的效率。通过对不同静压支撑槽参数下泵性能指标的对比分析，可以总结出以下影响规律：在一定范围内，适当增加静压支撑槽的宽度和深度，合理选择其位置，能够有效减小不平衡径向力，提高泵的性能。但需要注意的是，这些参数的调整都存在一个最优范围，超过这个范围，会导致泄漏量增加等负面效应，反而降低泵的性能。在实际设计中，需要综合考虑不平衡径向力、泄漏量和效率等因素，通过数值模拟和实验验证，确定最优的静压支撑槽参数，以实现高压内啮合齿轮泵性能的优化。3.5本章小结本章基于空化两相流动，深入研究了高压内啮合齿轮泵的不平衡径向力，通过建立空化流动模型、进行两相流动计算以及分析静压支撑槽的影响，取得了一系列有价值的研究成果。在空化流动模型建立方面，采用MIXTURE模型和Zwart-Gerber-Belamri空化模型相结合的方式，准确描述了泵内空化现象中气泡的生成、生长和溃灭过程。详细推导了空化模型的控制方程，包括MIXTURE模型的连续性方程和动量方程，以及Zwart-Gerber-Belamri空化模型中蒸发和冷凝过程的空化率表达式，并将其作为源项添加到MIXTURE模型方程中，以考虑空化对流动的影响。在计算策略上，明确了初始化过程中压力场、速度场和蒸汽相体积分数的初始赋值方法，设定了合理的收敛准则，包括残差收敛准则和物理量收敛准则，以确保计算结果的准确性和可靠性。同时，针对计算过程中可能出现的数值振荡等问题，提出了采用增加松弛因子和自适应网格技术等解决方法。依据实际工况和经验，合理确定了空化模型控制参数，如气泡半径取值为1\times10^{-6}m，成核位点体积分数取值为5\times10^{-4}，蒸发系数取值为50，冷凝系数取值为0.01，并通过敏感性分析深入了解了这些参数对计算结果的影响程度。在高压内啮合齿轮泵内的空化两相流动计算中，对比分析了考虑空化的两相流动计算对流量的影响。随着工作压力的升高和转速的提高，空化现象加剧，导致流量逐渐降低且波动更为明显。通过具体工况数据，直观展示了不同压力和转速下流量的变化幅度，如在工作压力从10MPa升高到20MPa时，考虑空化的流量降低幅度从[X1]%增大到[X2]%；转速从1000r/min增加到1500r/min时，相同工作压力下流量降低幅度从[X3]%增大到[X4]%。分析了考虑空化两相流动时的周向压力分布，在吸油区压力较低，空化现象相对较少；随着齿槽向压油区移动，压力逐渐升高，当压力降低到饱和蒸汽压以下时，空化气泡大量生成，在压力升高过程中气泡溃灭，导致压力分布出现明显波动，尤其是在靠近空化区域和压油区，压力分布的不均匀性更为突出，空化与压力分布之间存在着相互影响的关系。对于内啮合齿轮泵的不平衡径向力，通过数值模拟详细计算了不同工作压力下齿轮轴与内齿圈所承受的不平衡径向力。随着工作压力从2.88MPa增大到25MPa，齿轮轴和内齿圈所受不平衡径向力均逐渐增大，齿轮轴所受不平衡径向力方向始终指向齿顶与泵体内壁间隙较小的一侧，且在齿顶和齿根处相对较大，齿面中间部分相对较小；内齿圈所受不平衡径向力方向与齿轮轴相反，在与齿轮啮合区域相对较大，远离啮合区域相对较小，从内表面到外表面逐渐减小。不平衡径向力过大对泵的性能产生诸多不利影响，如加速轴承磨损、导致泵的振动和噪声增大等。在静压支撑槽的影响研究中，明确了静压支撑槽的关键尺寸参数包括宽度、深度和长度，以及其在泵体中的位置坐标确定方法。通过设计一系列不同尺寸和位置的对比方案，利用数值模拟详细对比了不同静压支撑槽参数下泵的不平衡径向力、泄漏量、效率等性能指标。结果表明，在一定范围内，适当增加静压支撑槽的宽度和深度，合理选择其位置，能够有效减小不平衡径向力，提高泵的性能。但这些参数的调整存在最优范围，超过该范围会导致泄漏量增加等负面效应，反而降低泵的性能。如槽宽增加时，不平衡径向力先减小后增大，泄漏量明显增加；槽深增加时，不平衡径向力和泄漏量也会发生相应变化；槽设置在靠近压油区的周向位置时，能更有效地补偿不平衡径向力，但泄漏量也会受到影响。综合考虑这些因素，以泄漏量满足要求，尽可能补偿不平衡径向力为目标，确定了内啮合齿轮泵的静压支撑结构形式。通过本章的研究，深入了解了高压内啮合齿轮泵基于空化两相流动的不平衡径向力特性，以及静压支撑槽对不平衡径向力的补偿效果和作用机制，为高压内啮合齿轮泵的优化设计和性能提升提供了重要的理论依据。但研究过程中仍存在一定局限性，如在空化模型中，虽然考虑了气泡的主要生成和溃灭过程，但实际空化现象可能更为复杂，模型的准确性仍有提升空间；在静压支撑槽的研究中，仅考虑了部分关键尺寸和位置参数，对于其他可能影响泵性能的因素，如槽的形状、表面粗糙度等，尚未进行深入研究，未来可进一步拓展研究范围，以更全面地优化高压内啮合齿轮泵的性能。四、高压内啮合齿轮泵不平衡径向力及其静压支撑的理论分析4.1内齿圈的不平衡径向力4.1.1沿齿轮圆周液体产生的径向力基于流体力学原理，内啮合齿轮泵工作时，内齿圈在液体压力作用下会受到沿圆周方向分布的径向力。假设内齿圈的半径为R，齿宽为b，液体在齿圈圆周上某点的压力为p(\theta)，其中\theta为圆周方向的角度。根据微元法，在齿圈圆周上取一微小弧段dL=Rd\theta，该微元弧段所受的径向力dF_p为：dF_p=p(\theta)bdL=p(\theta)bRd\theta对整个圆周进行积分，可得到沿齿轮圆周液体产生的径向力F_p：F_p=\int_{0}^{2\pi}p(\theta)bRd\theta此径向力的大小与液体在齿圈圆周上的压力分布p(\theta)密切相关。在吸油区，液体压力接近大气压，压力值相对较低；在压油区，液体压力升高至工作压力，压力值较高。从吸油区到压油区，压力逐渐变化，这种压力分布的不均匀导致了沿齿轮圆周液体产生的径向力的存在。齿轮的参数如齿圈半径R和齿宽b也会对径向力产生影响。齿圈半径R越大，在相同压力分布下，微元弧段所受的径向力dF_p越大，从而使得总的