大功率调速型液力偶合器温度场分析
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
大功率调速型液力偶合器温度场分析
刘建;高欣欣;范丽丹
【摘要】利用计算流体力学CFD(Computational Fluid Dynamics)方法对大功率调速型液力偶合器工作腔温度场进行数值计算,为了证明YOCQZ465液力偶合器具有较好的换热能力,建立两个几何模型:一个是YOCQZ465型液力偶合器流道模型,另一个是假设泵轮和涡轮上没有进出油孔的流道模型,通过对两者的计算结果分析可知,在泵轮外环内缘和涡轮外环外缘上开进出油孔能够加快工作腔内液体与冷却油的热交换,使液力传动油的温度被控制在合理工作温度范围内.将YOCQZ465型液力偶合器工作腔温度场计算结果与试验数据进行对比,结果表明,液力偶合器工作腔温度场的计算结果与试验数据基本吻合,证明所使用的数值计算方法具有一定的可靠性.
【期刊名称】《长春工程学院学报(自然科学版)》
【年(卷),期】2017(018)002
【总页数】5页(P54-58)
【关键词】液力偶合器;温度场;CFD;数值模拟
【作者】刘建;高欣欣;范丽丹
【作者单位】吉林工程技术师范学院,长春 130052;吉林工程技术师范学院,长春130052;吉林工程技术师范学院,长春 130052
【正文语种】中文
【中图分类】TH137.332
近年来,国内外学者对调速型液力偶合器内部流场做了大量研究,通过分析流场特性准确预测液力偶合器性能[1-6]。
但有关液力偶合器工作腔温度场的研究和取得
的成果非常少[7]。
在液力偶合器传递扭矩过程中,液力传动油不断克服黏性力做
功而转换成热能,液力偶合器传递功率越大,产生的热能越多,流体温升越高。
温度升高将导致液力传动油黏度降低,液力偶合器传递扭矩的能力也随之降低。
如果不能合理地控制油温,不仅会导致传动能力降低,而且会加速密封元件老化,导致漏油和渗油现象发生。
因此,研究液力偶合器工作腔温度场及热交换规律,尤其对于大功率液力偶合器具有十分重要的工程应用价值。
为了确保液力偶合器工作腔温度场符合液力传动油工作温度要求,必须保证液力偶合器工作液体外循环流量足够大,以使偶合器工作腔内的热油及时排出,冷却油及时充入。
YOCQZ465大功率调速型液力偶合器泵轮外环内缘周向均布23个直径
为10 mm的进油孔,涡轮外环外缘周向均布16个直径为18 mm的出油孔,如
图1所示。
这种结构能够保证液力传动油的循环流量足够大,从而实现工作腔内
部的对流换热。
为了更好地说明此种结构具有较好的降温能力,本文以
YOCQZ465调速型液力偶合器为原型,假设一种泵轮和涡轮上没有进出油孔的工作腔,对它和YOCQZ465调速型液力偶合器工作腔温度场进行数值计算。
1.1 控制方程
液力偶合器流道可以看作旋转的变截面曲线管道,液力传动油在泵轮与涡轮间周而复始不停地做螺旋环流运动传递能量。
调速型液力偶合器通常在部分充液状态下工作,因而其内部为气液两相流动,并且伴有涡流、回流、二次流等不规则流动现象,使流动恶化,增大流道内的液流损失,这部分能量全部转化为热能[8-11]。
本文利用计算流体动力学CFD进行数值模拟计算温度场,为了简化分析,做下列假设:
1)忽略各叶轮间的液流泄漏;
2)由于调速型液力偶合器的工作油温被控制在60~100 ℃范围内,通常采用的8
号液力传动油的黏温性较好,其物理性质在此温度范围内变化较小,因此,假设液力传动油的物理性质不随温度变化;
3)不考虑其他热源对偶合器内部的热辐射影响。
在以上假设基础上,利用两相流的连续性方程、动量方程、能量方程、湍流k-ε模型来描述液力偶合器内部液力传动油的对流换热过程。
1)连续方程:
▽
式中为质量平均速度;ρm为混合密度。
2)动量方程,通过对所有相各自的动量方程求和来获得。
▽▽P+[μm(▽vm+▽▽
式中:n为相数;F为体积力,N;▽为哈密尔顿算子;μm为混合黏度,pa·s;vdr,k为第k相的飘移速度,m/s。
第p相的体积分数方程为
▽▽
式中:αp为第p相的体积分数;ρp为第p相的密度,kg/m3;为第p相的飘移速度
3)能量方程
▽▽·(keff▽T)+Sh。
式中:E为总能量,J;Cp为流体比热容,J/(kg·℃);T为温度,℃;P为压力,Pa;keff为有效热传导。
4)湍流动能k方程:
湍流动能耗散率ε方程:
式中的常数一般取作:C1ε=1.44,C2ε=1.92,σk=1.0,σε=1.3。
1.2 模型建立
根据YOCQZ465调速型液力偶合器模型抽取其内部流道模型如图2所示。
图2(a)为假设工作腔内没有进出油孔的流道模型,图2(b)为YOCQZ465调速型流道模型。
由于液力偶合器泵轮和涡轮中的流动是非稳态的,工作液体处于湍流流动状态,存在强烈的脉动和漩涡,流体间进行着强迫对流传热过程[12-15],为了减少计算误差,对液力偶合器的全流道模型进行三维非稳态数值计算,泵轮流道和涡轮流道的交界面采用滑动网格技术。
将流道模型导入Gambit中划分网格后,设置边界条件,对于图2(a)设置泵轮和涡轮之间的交界面为网格分界面,其余表面设置为壁面;对于图2(b),除了把泵轮和涡轮之间的交界面设置为网格分界面,其余表面
设置为壁面外,还需要设置进油孔为压力入口,根据试验条件设置进口油温为
60 ℃,出油孔为自由流条件。
采用8号液力传动油为工作介质,计算泵轮转速为6 000 r/min,牵引工况下(i=0.8),充液率为qc=80%时液力偶合器工作腔的温度场。
1.3 计算结果分析
通过数值模拟计算得到液力偶合器工作腔温度场如图3~4所示,图3为假设泵轮和涡轮上没有进出油孔的液力偶合器,图4为泵轮和涡轮上开有进出油孔的YOCQZ465型液力偶合器。
图3中液力偶合器工作腔温度场分布规律,泵轮的温度梯度随着径向尺寸的增大
而递减,涡轮的温度梯度随径向尺寸的增大而递增,最高温度集中在涡轮叶片的压力面处。
8号液力传动油的最佳工作油温为60~95 ℃(温度场云图中的温度为开
氏温度,开氏温度等于摄氏温度加273),而图3中液力偶合器工作腔内最低温度
为77 ℃,最高温度为160 ℃,已经接近液力传动油的闪点。
图4中液力偶合器工作腔的温度场被控制在60~99 ℃范围内,符合液力传动油的最佳工作温度。
图4温度分布情况与图3有较大差别,泵轮的温度梯度总体上仍
然随着径向尺寸的增大逐渐递减,涡轮的温度梯度随径向尺寸的增大而递增,但在开有进油孔的流道内工作油温度较低,而没有进油孔的流道内工作油温度较高。
为了研究YOCQZ465调速型液力偶合器工作腔的热交换过程,本文截取具有进出油孔的流道截面温度场云图,如图5所示。
冷却油从进油孔涌入泵轮与内部流体进行热交换,使流道内工作油温度迅速降低。
由于涡轮外环外缘处出油孔径较大,流动阻力较小,部分高温流体能够从此处流出液力偶合器工作腔,而其余未能流出工作腔的高温流体将继续沿涡轮流道流动,进入泵轮后与低温流体混合,再一次进入循环流动。
另外,液力偶合器工作腔内存在大量的涡流、二次流等不规则流动,尤其在涡轮中,这可以使流体间的热交换进一步加快,所以涡轮流道内温度梯度没有泵轮明显。
从图5中可见,此出油孔油温分布在85~90 ℃之间,经CFD计算YOCQZ465调速型液力偶合器16个出油孔平均油温为85.8 ℃。
本试验对象为YOCQZ465调速型液力偶合器,用于热电厂给水泵的调速节能。
液力偶合器的输入轴与电机相连,电机输入转速为1 500 r/min,经前增速,泵轮输入转速为6 000 r/min,输出轴与给水泵相连。
如图6所示。
分别在液力偶合器的工作油进口和工作油出口处设置温度传感器,如图7所示。
数据采集时间选定为电网负荷相对较稳定的中午11:30—12:30,牵引工况(i=0.8)下,与数值模拟的条件相同。
试验所采集的液力偶合器进出口油温如图8所示,液力偶合器的工作油进口温度在60 ℃左右,出口温度在89 ℃左右波动,这种波动主要是由电网负荷的变化导致的。
CFD计算中设置的进口油温为试验数据的平均值,CFD计算出油孔油温平均值与试验数据基本吻合,但比试验数据偏低,这是由于在数值计算中省略了液力传动油在油路中产生的摩擦热和来自液力偶合器内部其他热源的辐射热。
1)本文利用CFD数值计算得到液力偶合器工作腔温度场分布的详细情况,通过对比分析可知,YOCQZ465型大功率液力偶合器具有较好的冷却能力,在泵轮和涡轮外环上开进出油孔能够加快冷却油与工作腔内高温流体的对流换热,使液力传动油温度控制在合理温度范围内。
2)液力偶合器的温度场计算值与试验值基本吻合,表明本文所采用的计算方法能够较准确地表述液力偶合器工作腔内换热过程,具有一定的通用性。
【相关文献】
[1] 何延东.基于CFD的大功率调速型液力偶合器设计[D].长春:吉林大学,2009.
[2] 何延东,马文星,刘春宝.液力偶合器部分充液流场数值模拟与特性计算[J].农业机械学报,2009,40(5):24-28.
[3] 何延东,马文星,邓洪超,等.基于CFD的调速型液力偶合器设计方法[J].农业机械学
报,2010,41(6):31-36.
[4] Lu xiu-Quan,Ma wen-Xing,Fan li-Dan,et al.Visualization experiment with PIV and analysis of flow field in hydrodynamic coupling[J].Applied Mechanics and
Materials,2010,29-32:1327-1333.
[5] Fan Lidan, Ma Wenxing,Lu Xiuquan,et al.PIV measurement of incompletely liquid-filled hydrodynamic coupling[C]//2010 IEEE International Conference on Information and Automation, ICIA.Harbin,China:IEEE,2010:1448-1451.
[6] 柴博森,马文星,卢秀泉,等.基于粒子跟踪测速技术的液力偶合器内部流速测定方法[J].农业工程学报,2011,27(7):140-145.
[7] 麻剑锋.旋转曲线管道内湍流流动结构和传热特性研究[D].杭州:浙江大学,2007.
[8] 流体力学.林建忠,阮晓东.[M].北京:清华大学出版社,2005.
[9] 金文,张鸿雁.灌水器内流道流场Micro-PIV试验分析[J].农业工程学报,2010,26(2):12-17.
[10] 陶文铨.数值传热学[M].西安:西安交通大学出版社,2004.
[11] 廖伟丽,姬晋廷,肖微,等.轴流转桨式水轮机转轮进口前流场的数值分析[J].农业工程学
报,2007,23(10):106-111.
[12] 范丽丹,马文星,柴博森,等.液力偶合器气液两相流动的数值模拟与粒子图像测速[J].农业工程学报,2011,27(11):66-70.
[13] 杨培志,顾小松,傅俊萍.YW25G型空调硬卧列车车厢内气流数值计算[J].中南大学学报:自然科学版,2005,36(5):888-891.
[14] 谢东,王汉青.不同气流组织下夏季空调室内热舒适环境模拟[J].建筑热能通风空
调,2008,27(3):66-69.
[15] 袁文麒,刘遂庆.管道充水工况下气液两相流瞬态数值模拟研究[J].同济大学学报:自然科学版,2010,38(5):709-715.。