限矩型液力偶合器全充液工况的流固耦合研究*

2020-05-25任豪宗王时英张嘉华

机电工程 2020年5期

任豪宗，王时英*，吕明，张嘉华

（1.太原理工大学机械与运载工程学院，山西太原 030024；2.精密加工山西省重点实验室，山西太原 030024）

0 引言

液力偶合器是一种以液体为工作介质来传递能量的动力装置，主要分为普通型、限矩型和调速型。限矩型偶合器又被称作安全型偶合器，主要用于解决启动困难、过载保护和减缓冲击等问题。如果设备在在高速运转过程中，出现突然卡住的现象，会产生瞬时的过载，此时动力机和工作机所有的动能都被突然释放出来，将造成原动机和工作机的破坏。因此，就需要限矩型偶合器来进行过载保护，避免原动机和工作机受破坏。另外，原动机经常有载甚至重载启动，采用偶合器可大大改善电动机和工作机的启动性能[1]。

国内外许多学者已对偶合器进行了广泛的研究。何延东等[2]利用UG和CFD，针对调速型液力偶合器的3种结构优化方案进行了数值模拟，将结果与基型进行对比得到了合理的优化方案；刘春宝等[3]采用液力变矩器瞬态模型以及多流动区域的耦合算法，深入分析了液力变矩器的流动特性并分析了原因；邹波等[4]为了分析液力减速器在高转速下的叶片强度，采用单向耦合技术和全局守恒差值法，对减速器全充液工况进行了流固耦合分析，得到了减速器应力分布云图和变形云图，为液力元件的有限元分析提供了有效的方法；魏巍等[5]研究了变矩器叶片强度，采用流固耦合的方法分析了工作液体对叶片产生的应力与应变情况，并根据材料特性对某型液力变矩器失效情况进行了预测。

本文将在D483限矩型液力偶合器的基础上对偶合器全充液工况下3种典型转速比工况各自的流场情况和叶片受力特性进行分析，并将叶片工作时的最大受力值与叶片材料的极限强度进行对比，以排除偶合器因为工作液的正应力而破损的可能，并为偶合器的改进提供一种有效途径。

1 流固耦合的分析流程

流固耦合流程图如图1所示。

图1 流固耦合流程图

图1中的具体流程：

（1）首先利用SolidWorks建立偶合器几何模型，运用布尔运算提取几何模型的互补模型作为流道模型；

（2）然后对流道模型进行网格划分，利用ANSYS Fluent进行流场数值模拟得出速度、压力场及矢量图；

（3）再将流场分析数据导入静力学分析软件中，对叶轮几何模型进行网格划分、边界条件的设置以及离心载荷的加载；

（4）最后通过耦合求解，得出各工作轮叶片所受到的的应力和应变情况，并根据耦合结果对各工作轮叶片的强度进行对比分析。

2 流道模型与几何模型建立

根据已有的偶合器模型，笔者建立限矩型液力偶合器几何模型，并提取流道。考虑到偶合器的循环结构，选取周期性模型作为计算区域，并对偶合器模型作一定的简化。

流道模型和几何模型的网格划分如图2所示。

图2 流道模型和几何模型的网格划分

为保证计算的准确度，笔者采用高质量六面体网格作为网格单元。流道模型的网格数目为116 226，几何模型的网格数目为85 103。

3 流场数值模拟

3.1 基本假设和边界条件

基本假设如下：

（1）液力偶合器工作温度为等温工况，液体的内能不变；

（2）工作腔完全封闭，工作液不发生泄露；

（3）相同工况，叶轮的每个流道特性相同[6-9]。

边界条件如下：

（1）根据实际工作情况，泵轮转速恒为1 900 r/min；

（2）由于偶合器循环结构，流道的循环切割面采用周期性边界条件，其他表面均采用壁面边界条件；

（3）采用SIMPLE算法和一阶上游迎风格式，计算收敛准则设置为迭代的残差小于10-3。

3.2 仿真及结果分析

笔者对制动工况i=0、牵引工况i=0.5、额定工况i=0.97这3种典型工况的瞬态特性进行了仿真分析，得到了速度和压力分布。

压力云图如图3所示。

图3 压力云图

由图3可知：

（1）偶合器在转速比i=0时，其工作特点是：涡轮不转动，不输出动力，η=0。从压力云图可知：此时工作腔内压力分布没有规律，出现大面积负压区域，这是造成能量耗散以及效率低的主要原因，也是偶合器结构优化的重点方向；

（2）随着转速比的提高，泵轮与涡轮的转速差减小，涡轮对来自泵轮的高速工作液的阻滞作用减弱，使得流动更加趋于规律性，因此，在牵引工况i=0.5时，流道内的压力出现较均匀的压力梯度分布；

（3）当转速比i=0.97时，泵轮和涡轮都具有较高转速，产生较大的离心力，压力场逐渐增大并出现层状分布，且越靠近流场外环产生的压力越大。

速度及矢量云图如图4所示。

图4 速度及矢量云图

由图4可知：

（1）i=0时，工作液体在离心力的作用下沿泵轮壁面向外环作高速运动，之后进入处于静止状态的涡轮，在涡轮处受到运动阻滞，只有一部分工作液沿着涡轮壁面运动，速度有所降低，而涡轮内远离壁面的工作液速度急剧降低，并且出现回流、旋涡以及大面积的低速区，无法参与到液力传动的大循环中，导致能量损失，效率低。如图4（a）中的A处出现了回流、旋涡等；

（2）i=0.5时，涡轮与泵轮的转速差减小，涡轮对来自泵轮的高速工作液的阻滞减弱，同时涡轮产生一定的离心力，大部分的液体能够参与到液力传动的大循环中，仅在涡轮中间位置出现小面积的低速区，能量损耗小，传动效率提高；

（3）i=0.97时，涡轮与泵轮的转速基本一致，涡轮入口的速度与泵轮出口相同，因此工作腔内的速度沿径向层状分布，且越靠近外环压力越大[10-11];

（4）在制动工况时涡轮工作腔内的矢量图如图4（d）所示。由图4（d）可知：B处和C处分别出现了旋涡以及回流等损耗能量的现象，这些现象导致了偶合器工作效率的降低，可以通过改变叶轮的数目、辅腔的体积以及过流孔的数目和直径等，来优化偶合器的结构，提高偶合器的限矩性能。

4 叶轮强度分析

4.1 网格划分及边界条件的设定

本文模型为限矩型液力偶合器，具体参数为：

泵轮叶片数为46，涡轮叶片数为46，叶轮材料为铝合金（ZL101A），工作液为水，密度ρ=998.2 kg/m3,泵轮恒转速为1 900 r/min。

液力偶合器为循环结构，只需取周期结构边界条件作分析。为了计算结果的准确性，笔者同样选取高质量的六面体网格。

在全充液工况时，高转速液力偶合器的工作载荷主要来自流体作用下的壁面正压力、在高速旋转下叶轮产生的离心载荷以及工作液在能量转化过程中消耗能量产生的热量。由于偶合器在稳定的工况下的温度基本不变，产生的热量可以忽略不计。

对限矩型液力偶合器的叶轮结构进行强度分析，常采用流固耦合（FSI）的分析方法。流固耦合法一般分为两种：单向流固耦合和双向流固耦合。单向流固耦合是指在液力传动过程中，只考虑工作液运动对叶片的作用，而不考虑工作液体与结构相互之间的作用，即在耦合界处的数据传递是单向的；同理，双向流固耦合要考虑流固耦合的相互作用，即耦合界处的数据传递是双向的。

本研究主要研究工作液对叶片的影响，因此，选用单向流固耦合分析。

设置边界条件：

在偶合器几何模型的循环切割面施加循环约束；分别在制动工况i=0、牵引工况i=0.5以及额定工况i=0.97下，施加泵轮和涡轮对应的离心载荷；

载荷施加：

在限矩型液力偶合器全充液工况下，对转速比i=0、i=0.5、i=0.97这3种工况进行CFD数值模拟，得到液力偶合器内部流场的压力和速度分布，接下来将计算数据导入，并施加于液力偶合器，对其进行强度计算。

4.2 仿真及结果分析

笔者在Workbench的Static Structural中进行流固耦合的仿真分析，得到了3种典型工况下叶轮的总体变形量和等效应力。

总体变形量和等效应力如图5所示。

图5 总体变形量和等效压力

根据图5，从3种工况下的总体变形量和等效应力的对比可以分析出：

（1）制动工况时，最大等效应力为33.87 MPa；（2）牵引工况时，最大等效应力为19.55 MPa；（3）制动工况时，最大等效应力为12.65 MPa。

随着涡轮转速的提高，泵轮和涡轮的转速差变小，转速比变大，叶轮对工作液的阻滞作用越来越弱，工作液对叶轮壁面产生的正应力逐渐减弱，叶轮的最大变形量和等效应力逐渐减小。同时，偶合器内部的压力趋于均匀，并且沿径向出现层状分布特性，越靠近外环，压力越大，D处泵轮出口和涡轮进口的交界，此处为较大应力区，是最容易发生损坏的部位；涡轮E处一直都是受工作液应力较大位置，并且随转速比的提高，应力越来越大。

从流固耦合的结果可知：在制动工况i=0时，整个叶珊系统的等效应力最大，最大等效应力为33.87 MPa。

液力偶合器的正常工作温度为90 ℃～100 ℃，经查阅相关材料手册，可得叶轮材料的强度极限σb为270 MPa，循环交变载荷下叶轮的疲劳强度为119.5 MPa。将其与最大等效应力对比可得：叶轮材料的强度能够满足正常工作要求。

泵轮局局部破损图如图6所示。