限矩型液力偶合器内部两相流动数值模拟_唐欣怡.pdf

年第 卷 月第 期机 械 科 学 与 技 术 ：收稿日期：作者简介：唐欣怡（），硕士研究生，研究方向为流体机械及工程，通信作者：李志鹏，教授，硕士生导师，唐欣怡，李志鹏，王阳，等限矩型液力偶合器内部两相流动数值模拟机械科学与技术，（）：限矩型液力偶合器内部两相流动数值模拟唐欣怡，李志鹏，王阳，李铁梅，孙军红（长沙理工大学 能源与动力工程学院，长沙；湖南中特液力传动机械有限公司，湖南益阳）摘要：限矩型液力偶合器内部存在着复杂的气液两相流动，为掌握限矩型液力偶合器内部流场分布及转矩特性变化，分别采用（）与 两相流模型对偶合器内部流场进行模拟计算，并对其转矩特性进行监测。通过 模型获得的偶合器内部气液两相分布情况和通过 模型得到的其内部压力速度分布图能够较好地反映偶合器内部流场变化情况。仿真结果表明：模型能够较好地模拟出偶合器因流态转变而造成的转矩跌落情况，模型不能模拟出该效果，但在高、低转速比工况下模拟的转矩值仍具有一定参考意义。关键词：限矩型液力偶合器；气液两相流动；转矩特性；两相流模型中图分类号：文献标志码：文章编号：（），（，；，）：，（），；，：；动压泄液式限矩型液力偶合器主要通过调节其内部充液率实现轻载启动、过载保护、协调多动力机运行等功能，内部流场特性决定了其转矩输出特性，因此研究限矩型液力偶合器内部流场对偶合器内部流道结构优化具有重要指导意义。目前针对限矩型液力偶合器内部流场的研究方法主要有实验与（）数值模拟。国内外学者主要通过实验方式对液力偶合器内部流场进行可视化研究，等通过 射线断层扫描成像技术对 充液率不同转速比的偶合器内部流场进行实验，得到其工作腔不同断机 械 科 学 与 技 术第 卷：面的气液两相分布情况。等通过在泵轮叶片两面安装平面阵列传感器的方式，得出了 充液率下转速比 、时泵轮叶片上压力面和吸力面的气液两相分布情况。李兴忠等利用（）技术对偶合器内部流场进行测量，分析了在不同工况条件下不同流道区域的速度场和涡量场。卢秀泉等结合 技术对液力偶合器进行了可视化与外部特性同步测试实验，对调速工况下偶合器涡轮内部瞬态流动特性进行了研究。由于实验周期较长、平台搭建成本较高且操作较为复杂，所以 技术在偶合器内部流场模拟研究方面得到了广泛的应用。等采用标准 湍流模型，利用 平均法对矩形流道的液力偶合器内部流场求解，发现仿真结果与实验测试结果吻合度较高，且解释了液力偶合器扭矩传递的物理过程。等利用 研究方法分析了影响限矩形偶合器性能的参数，奠定了将 技术运用到液力偶合器优化的基础。等利用 技术模拟液力偶合器两相瞬态流场，研究了不同充液条件、不同的叶片倾斜角度和速度比下液力偶合器的流动特性，得出叶片角度会影响流场速度分布情况、液相体积分数会影响传动性能的结论。张德生等将 技术应用到液力偶合器的流道腔型优化中，设计提出优化后的流道结构参数。赵继云等采用 模型对阀控充液型液力偶合器在不同工况下的气液两相分布情况进行了模拟，揭示了阀控充液型液力偶合器内部环流变化规律，并分析了挡圈对环流形态的影响。王阳等利用 模型对液力偶合器内部流场进行瞬态分析，总结了其流场 结 构 的 变 化 与 规 律。卢 秀 权 等利 用 模型实现了对限矩形液力偶合器内部流场分布情况的预测，且通过试验获得的外特性对比，对数值模拟方法的适用情况进行评价。目前，针对限矩型液力偶合器内部两相流的研究主要基于单个两相流模型，不同工况下的偶合器内部气液两相流场差异较大，单一的两相流模型并不能完全适用于全工况限矩型液力偶合器内部流场进行的模拟，因此本文分别利用 与 两种两相流模型对 型限矩型液力偶合器工作腔进行瞬态分析，通过 模型得到其气液两相分布，通过 模型得到其压力速度场分布，从而揭示其内部流场特性，并得出两种两相流模型的适用性规律，为偶合器内部流场模拟提供方法。物理模型的建立本文选取 型限矩型液力偶合器为研究对象，对泵轮、涡轮及前辅腔之间工作腔中两相运动情况进行研究。此型号液力偶合器循环圆直径为 ，泵轮 枚叶片，涡轮 枚叶片，将其整体流道结构提取如图）所示，网格模型如图）所示。工作腔内油液为 号液力传动油，设置其密度为 ，黏度为 ，油液与空气之间的张力系数设置为 。图 几何模型为实现对 型液力偶合器工作腔内部的数值求解，结合偶合器的工作流场特点，首先假设：）偶合器工作腔为密闭结构，无工作油液的流入流出；）偶合器叶轮为刚体，工作中无变形，流道空间不改变；）工作油液的黏性系数不变；因为泵轮、涡轮及前辅腔之间存在复杂的耦合运动，为保证数值计算的精度及收敛稳定性，本文使用高质量的六面体网格对偶合器流道进行划分，如图）所示，网格数量为，网格质量均大于。两相流模型分别选用 与 模型，湍流模型采用（）方第 期 唐欣怡，等：限矩型液力偶合器内部两相流动数值模拟：程，压力、速度耦合方式采用（）算法对偶合器瞬态流场进行数值计算。计算分析模型的建立 两相流模型限矩型液力偶合器内部的流动是三维的、不可压缩的非定常两相流动，选取的两相流模型是否合适对模拟结果的准确性有着较大影响。本文分别选用 模型与 模型对限矩型液力偶合器内部流场进行瞬态分析。模型是一种简化的多相流模型，可以通过求解混合物的动量、连续性和能量方程、次相的体积分数方程和相对速度的代数表达式来模拟相。其连续性方程为（）（）式中：为质量平均速度，；为混合物密度，；为第 相体积分数。通过将各相的动量方程求和来获得混合相的动量方程，其可以表示为（）（）（），()（）式中：为相的个数；为体积力，；为混合物黏度，；，为次相的漂移速度，。模型是在固定欧拉网格下的一种表面跟踪方法。前提是两种或多种流体不会相互混合。在 模型中，在连续性方程中引入了相体积分数这一变量，从而实现了对计算单元相界面的追踪。在每个控制体中，将第 相流体的体积分数记为，其连续性方程为（）（）（）（）式中：为第 相的密度，；为第 相的速度，；为从 相向 相的传质，（）；为从 相向 相的传质，（）；为源项，其默认值为；为算子符号；为时间，。在 模型中，求解动量方程所得到的速度场被各相共用，动量方程取决于通过属性 和 的所有相的体积分数。（）（）（）（）式中：为相密度，；为动力黏度系数；为速度，；为压强，；为重力加速度，。滑移网格模型偶合器在运行过程中始终存在滑差，为了更好地模拟泵轮与涡轮之间的耦合作用，采用滑移网格模型对泵轮与涡轮区域进行数据交换。运用滑移网格时，计算域至少包含两个或多个存在相对运动的子域，每个子域存在一个或多个与其他相邻子域的交界面，称为，本文 设置如图 所示。图 泵轮和涡轮流道的 示意图 初始化在使用 模型求解时，因为偶合器内部始终处于部分充液，因此需要对其进行初始化设置。恰当的初始化方法能够加快收敛速度，提高模拟准确度，本文试用了两种初始化方式，并对其收敛情况进行比较。方法一是以真实初始状态确定初始化气液两相分布，如图）所示。方法二是根据对结果的预估进行初始化设置，定义流道外围为液相，流道内部为气相，如图）所示。以转速比 充液率 的情况为例，分别使用两种初始化方法对液力偶合器内部流动进行模拟，并对其收敛速度及结果进行比较。图）为运用方法 初始化得到的泵轮转矩图，图）为运用方法 初始化得到的泵轮转矩图，可以看出运用方法二模拟所得的泵轮转矩在约在 时保持稳定，不再发生过大波动，从收敛时间和收敛情况来看都优于方法一所得到的结果。所以为保证模拟准确度和加速收敛过程，本文选用方法二进行偶合器流场初始化。机 械 科 学 与 技 术第 卷：图 充液率下初始化图图 转矩监测曲线 流场分析与特性预测限矩型液力偶合器内部气液两相流动较为复杂，在不同转速比和充液率下，产生不同形态的环流，从而改变偶合器外特性输出。本文模拟了、充液率的偶合器在 种典型工况下的工作腔中流场运动情况，通过分析使用 模型模拟得出的气液两相体积分布云图与使用 模型模拟得出的压力分布云图和速度矢量图来揭示限矩型偶合器的流场运动特性。启动工况图 为转速比 时的启动工况下，种典型充液率（，）的气液两相体积分布云图。图 种充液率下轴向截面气液两相体积分布云图（）图 中蓝色区域为气相，红色区域为液相，其它颜色为不同混合程度的气液混合物。在启动工况下，大部分油液停留在前辅腔中，参与流动的小部分油液由泵轮进入涡轮后，受到涡轮离心压力阻抗的作用很快折回泵轮，且高速流动的油液会受到涡轮叶片的阻碍，在涡轮入口处形成油液混合空气的漩涡，因此涡轮流道内出现大片由液相包裹着的气相区域。随着充液率的增加，液相中夹杂气相的比例降低，气相区域面积随之减小。启动工况下，泵轮高速旋转，涡轮静止，由图 可以看出，泵轮出入口与涡轮出入口有着明显的压力集中，且此时腔内压力分布不均衡，出现大面积负压区域。此工况下工作腔内的油液做大环流运动，沿壁面处流道有着较高的压力分布，气体被包裹在泵涡轮中部，因此这些区域出现明显的负压，且负压区偏向涡轮侧分布较大。第 期 唐欣怡，等：限矩型液力偶合器内部两相流动数值模拟：图 种充液率下轴向截面压力分布云图（）图 为启动工况下的速度矢量图。由于此时泵轮出口、涡轮入口处与涡轮出口、泵轮入口处两侧压力差达到最大，因此这两处速度值也达到最大。偶合器内介质在经过泵轮牵引及离心力的作用下以较高的速度进入涡轮，在经由涡轮外壁做环流运动后大部分进入前辅腔，少部分进入泵轮开始下一轮循环流动。由于高速旋转的泵轮对回流介质产生的阻碍作用，因此在涡轮内形成了一个涡流区域，整个工作腔内呈现出一种大环流的流动状态。图 种充液率下轴向截面速度矢量图（）牵引工况牵引工况时，涡轮转速升高，泵轮与涡轮间的滑差减小，涡轮对油液的离心阻抗减小，如图 所示，种充液率下工作腔内都出现明显的气液两相分界面。在低充液率时，油液在泵轮腔上部及涡轮腔内大部分区域循环，气体主要分布在前辅腔及泵轮中下部。在充液率为 时，泵涡轮腔内充满大量油液，随着环流区的扩大，部分油液由涡轮工作腔进入前辅腔形成环流，部分气体与油液混合，分布在流道外沿参与循环。图 种充液率下轴向截面气液两相体积分布云图（）机 械 科 学 与 技 术第 卷：图 为牵引工况下的压力分布云图，图 为速度矢量图。随着转速比的增加，工作腔内的压力逐渐趋向于带状分布，泵轮与涡轮之间的两侧压力差减小，因此最大速度值会相对减小，但最大速度依然会出现在两个压力差较大的地方，即泵轮出口、涡轮入口处与涡轮出口、泵轮进口处。同时转速比的增加使得泵轮与涡轮间的阻碍作用降低，泵轮与涡轮内的涡流减弱，因此，油液在进入涡轮上部大圆弧段速度达到最大值之后受到壁面阻碍，速度逐渐减小，且涡流中心逐渐向工作腔中心移动，环流形态由大环流慢慢向小环流转变，被环流包裹的油液滞流区域逐渐缩小，从速度矢量图可以看出，滞流区随充液率的升高逐渐缩小。图 种充液率下轴向截面压力分布云图（）图 种充液率下轴向截面速度矢量图（）额定工况当偶合器达到额定工况时，进入涡轮的油液受到涡轮离心压力阻抗很快折回泵轮，此时腔内为小环流形态。如图 所示，种充液率下偶合器气液分界面都十分明显，且几乎呈现水平，这是由于泵轮与涡轮间滑差很小，液体在泵轮和涡轮中所受离心压力的差值也很小，所以气液分界面没有发生明显倾斜，且可以看出，气体全部位于工作腔内侧，几乎无气液混合物参与循环，此时偶合器已经达到平稳运行状态。图 种充液率下轴向截面气液两相体积分布云图（）第 期 唐欣怡，等：限矩型液力偶合器内部两相流动数值模拟：图 为额定工况下的压力分布云图。此时腔内油液环流作用减弱，环流形态稳定，压力呈现出基本水平的层状分布，且与工作腔半径呈明显线性关系，并随充液率的升高流道整体压力增大。图 种充液率下轴向截面压力分布云图（）图 为额定工况下的速度矢量图。在额定工况时，泵轮与涡轮趋于同步运转，油液进入涡轮后能量损失较少，所以泵轮出口速度和涡轮出口速度基本相等。此时泵轮与涡轮间压力基本相同，因此工作腔内介质流速较小，且 种充液率下的速度分布和速度值无明显区别。由图 可看出，到达额定工况后，偶合器内部环流区域减小，这是因为内部工作介质从泵轮流出后，由于涡轮此时转速较高，