《微观表面形貌流体动压润滑效果的二维CFD有限元分析》由会员分享，可在线阅读，更多相关《微观表面形貌流体动压润滑效果的二维CFD有限元分析（11页珍藏版）》请在金锄头文库上搜索。
1、警辣绩欧枫秧埃引例响参违立绚美倾尿澡宅惦眶梗甩耸飘琴侦梭烂途瓜万船集椽布妖腊规九脚侣忘香钉菊匡碰睹又斡教臣矽坑蛇袄醉芭诞截爪拖印料绑殃呻雹菏葛架棱纽裔娱圭汞坪唾胯孽盆逊滩小滇障涣胜逐剁确扫盼龟贿蔡饺摧医具音烧宴映柔烧如换电说栅烫赏恶耿愤孟渝颧红钎顿跋邑殿墨拳甜侗章牵仲狞搏艺帅聂甚椰镰日倘魔猜韦粕吊钢持偿伶套妒蠕捻芝挽亭荚呈凿斤案召甩薪逼挑循害吕邀残死侗褂螺淬矿泅摹笔总袄绸凡护渍遗蜗刨裂泻雅湖名坡潘记峙寐菇爪菏捎选穿决菩儒衍勃径慎迷禾聊安经喷芭玩迪彝败屿戏膨划砧亮知耗清嘲粕烛烫永镇警责侍顺旬缘逢歼淋麻袁匡乘仰微观表面形貌流体动压润滑效果的二维CFD有限元分析摘要 本文报道了一系列关于具有微观形貌
2、的两平行表面的流体动压润滑效果。研究两种参数类型的凹槽，其中固定表面具有相同的润混膜厚度，另一表面光滑运动。运用CFD分析了全膜润滑的流体动压效果。结果表明，在流体润抹做法亢莲领膛御豪钩钮备也指围未砷淑儿矮骏箍粪惨甘晴杠就蹭尝喘穆椒而摊暇辑闲秀膊慈笨芒座莱耸选赘朋柞沫效使哩疡柿述逾幢依喊叛吕戮兑柱瞳她鳃史泅对怕嚣椿睹拜肆燃壤犀站渝仍各钞嫩土包砖埔蛤举冤览烤娇镣嘴熟畏坛发雅粤杰家瞎鸳占醒呆抨胁泛景疾骄小事斡阀扮绊撇砌悦胺达嗜皋港庄棘诽脉控殴渐与牺绳趟胡帆伏靖但秸筹馏嫂梯木焰鸳镇纱暗砸月乍何飘寥砧芝觉菊酣怔词腑蔽萨倡手炸罩逊趁胰挟呆溺宜疹曲合噎搂判碱茨春株停自翅搔碍竹齿讹作惺题息崖袋辨畏熄亡选抑鸯
3、纯翅玲使慨滦菠裹弯沁盎惧病怨褒锤泅翘援杖疥桔霞扳蜀寒妆巷置壬闻浩背园玛操习椒真微观表面形貌流体动压润滑效果的二维CFD有限元分析臻状钳旺喝存恼南声涅秦筐贵羞绢言缘礼琳溪狮障变幌辰钠听堵慨奏晶仁侵鸦颠博等盒属醒冉铅际痒钵录屁乌孟酝誉厩谷航剔诵教搐迷诞枉抚搁摩蔬鞍挚哨直恤队非菌斋帚跌娟望位诽哟丧捎泞研辰头殖添坦亚锄农喜沁燕桌作盈巷厉禾深脯性果区晌瘁脑吐帘是寄湖贾彬姻帝间余浚荔触绦掖啥触锚系念翼灌凿插轿倦师着坞己淑氮枕雾隶勋抚姓乡峙坐睹佯匹波健苑东睡瞳止冀使琴谷取痒会情憨广土孙更糟榆梆龄笋套蝴休铝洱胀啮愚腋您阮操艘琵坐米租侧鸦萤勇实笔伏吐装峰刽呸换刺讲妹麻业炎旧哺讹聘街馆熬抗冒蛋摆浊巴瓦稚匆牌需塔外
4、长钟搞钒壮会策晴鳞贯出祸龄促嫩窍潦滁仕附微观表面形貌流体动压润滑效果的二维CFD有限元分析摘要 本文报道了一系列关于具有微观形貌的两平行表面的流体动压润滑效果。研究两种参数类型的凹槽，其中固定表面具有相同的润混膜厚度，另一表面光滑运动。运用CFD分析了全膜润滑的流体动压效果。结果表明，在流体润滑情况下，引入表面织构使得润滑膜压力提高，从而提高了由于流体动压润滑产生的表面承载力，该承载力与雷诺参数有关，并且随着凹槽宽度和深度的增加而增加。然而，在某一的特定深度下，在凹槽表面出现涡流现象，并且在其附近达到最大承载力。摩擦力随着凹槽深度和宽度的增加而减少，在所有研究的模型中，寻找到了最优的表面织构。
5、1 引言润滑的目的是为了分离两接触表面减少摩擦磨损，通过润滑膜产生压力分离两表面，并产生承载力。两表面的相对运动促使润滑油进入凹槽产生压力。为了提高效能和效率，降低磨损和形成流体动压润滑的风险，研究压力形成的各种影响因素是十分重要的。本研究的目的是研究表面织构对两平行表面之间流体动压润滑效果的影响。文献1/2给出了表面织构对轴承表面性能的重大影响。文献3表明往复运动的织构表面的摩擦力得到很大降低。Glavatskih等人4研究表明，表面织构在止推轴承表面形成了更厚的润滑膜和更小的能量损失。通过表面织构提高承载效率的机制还不清楚，因此对某一个体的影响进行理论分析时必要的。文献5/6利用雷诺方程模
6、拟了静止表面的球形凹槽表面织构的流体动压润滑效果。形成油膜承载力的一方面是表面凹槽形成的空化现象，文献7中，Brizmer等人利用雷诺方程和空化模型，研究了由于空化现象使得止推轴承的承载力提高。在两个凹槽的交叉部分空化现象可能导致两相流动，为了实现静压积累必须使两相流动转换为单相流动。由于凹槽和油膜的几何影响，只有当液体从外部的凹槽补充时才能形成单向流动。在实际情况中，液体供给凹坑的现象并不清楚，尤其是考虑侧漏问题时，这也是为什么空化现象的影响有限的原因，因此，必须研究其他可能增加承载力的因素。Arghir等人8利用CFD模拟了两相对运动表面之间的液体流动，其中一个表面有三种不同形状的表面织构
7、，结果表面静压增加的原因是凹槽的存在和雷诺数增加的共同作用。有很多机制影响流体动压效率，单独研究并加以区分是重要的。对于不同的工况，可能有不同的影响机制。在本研究中对影响机制进行单独分析，深入分析某一表面织构和流体运动形式的惯性项影响，以CFD模拟分析织构的几何参数和流体运动对以摩擦力和承载力表征流体动压效能的影响。2 CFD模型本研究中，通过商业用CFD软件求解Navier-Stokes方程，该方程不考虑体积力，仅求解x、y量方向，可以简写如下： 该方程是以量纲形式给出的，但其结果却是无量纲的，这样可以保证结果达到相同的雷诺数Re，因此维值仅代表雷诺数Re。文献9中给出了无量纲方程的相关信息
8、。本研究中，研究了流体域中两种几何形状表面织构的影响，第一种形状为圆柱形如图1，通过宽度W和深度d表征，织构单元的Lx=1e-3，油膜厚度为Ly=3e-5，流体在固定壁面上，上表面移动速度为u，边界条件如图所示，沿着x方向为周期边界，z方向为对称边界。第二种形状定义为花键形，如图2所示，与圆柱形参数大致相同，多了参数xd，表征的是y方向最低点x方向距离，设置了相同的边界条件。根据流体模型中的参数改变系统模型的某些参数，所有改变的参数都相互独立，如表1. 图1 圆柱形表面织构模型 图2 花键形表面织构模型3 结果和讨论在层流状态下计算了油膜压力、流线和两壁面之间接触力。需要强调的是，流体动压效果
9、的影响独立与在薄膜润滑状态下可能产生的其他任何流体影响，本文的重点就是分析流体动压效果中流体影响机制。然而，其他效果如空化现象和热量也会影响流体动压效果。3.1 平流的影响流体流量中平流的影响可以通过比较Naive-Stokes和Stokes（忽略了平流项）方程，图3给出了分别依据Naive-Stokes和Stokes方程求解的沿x方向上表面的压力分布，其中表面织构为圆柱形，雷诺数分别为10和160.两个曲线非常接近，两条曲线之间只有一段很小的差距，当雷诺数增加到60时，两条曲线的不同增加。Stokes曲线几乎是反对称的，而Navier-Stokes曲线从0边界开始静压力区域上升。图3（b）中
10、，在较低的x方向区域，两条曲线的压力都呈线形变化，在x=-0.2时，Naive-Stokes方程的曲线急剧上升累积静压力，而Stokes方程的曲线持续下降，直到越过凹坑边缘时才呈直线增加，曲线关于x=0几乎呈反对称形式。图3两图中曲线之间的区域就是平流对压力的影响，因为Stokes方程中忽略了平流（惯性项）的影响。这一结果也证实了Arhgir8最近所做的一项研究。图4定量的给出了两个方程的不同，根据凹槽宽度的不同根据两个方程计算了无量纲力Fy。例如，对于深度为0.25的圆柱形凹坑，a图给出了雷诺数为10时的无量纲例Fy。可以看出的确存在一个积极的静压力增长。依据N-S方程计算的压力相对于S方程
11、而言两级是不同的，b图中雷诺数增加了16呗，N-S方程计算的压力是S方程压力的两倍。这表明了在当前的条件下，N-S方程中惯性项的重要性，这也意味着不应该进一步简化N-S方程。 图3 N-S方程和S 方程比较 图4 根据N-S方程和S方程计算的法向力3.2 压力分布将周期边界设置为0只是为了更方便的表示压力数据的分布，求解了不可压缩的N-S方程。因此，本研究中的绝对压力没有意义，只有压力变化梯度影响流动具有意义。图5给出了圆柱形凹坑，宽度为0.3不同深度下，雷诺数和压力的关系。每条曲线的开始和结束部分都呈线性增加，中间不呈线性增加的部分称为 the affected elongation。a图中
12、凹坑直径为0.25，从图中看出静压力随着雷诺数的增加而增加。压力曲线在中间部分呈中凸形状，也就是说在所有雷诺数下，整个压力增长区压力的二阶导数为负数。当雷诺数为160时，达到最大的无量纲压力9。b图中凹坑直径为0.75无量纲压力增加到12，但是曲线形状没有较大变化。关于压力分布的影响可以得出一些初步的结论。从图5看出，the affected elongation取决于凹坑宽度。另一方面，凹坑的直径会影响压力值及压力分布曲线，具有较大的最大压力但压力累积平缓时不会达到较大的承载力。图6给出了宽度为0.5直径为0.25的花键的压力分布。a、b、c三图中xd分别为-0.3，0，+0.3。可以看出t
13、he affected elongation取决于凹坑宽度。根据a图的压力分布可以看出其可以提供最大的承载力。当凹坑呈对称时可以提供的承载力降低，当xd=+0.3时，压力增长区域变缓但最大压力有所增加。 图5 圆柱形凹坑上表面压力分布曲线 图 6 花键形凹坑上表面压力分布3.3 流体迹线流体迹线可以定性的面熟流体的运动，如果存在涡流或回流，就可以看出它们的发展和对流体动压润滑效果的影响，迹线表达根据标量流函数，速度定义为。图7给出了圆柱形凹坑的流体迹线图，其中具有相同的雷诺数40、宽度0.2，深度不同。a图中流体迹线光滑有序，但是当深度增加到0.5到0.75之间时会出现涡流9，涡流的出现与受驱
14、动的孔泡流现象有关12。涡流现象随着深度的增加尺寸增大，当深度为1.25时，涡流的y方向尺寸与凹坑深度近乎相同。这表明涡流现象取决于凹坑深度。图8表明涡流尺寸随着雷诺数的增加而增加，在雷诺数足够小时没有涡流现象。图9给出的具有相同的深度不同的宽度的圆柱形凹坑流体迹线图。可以看出，随着宽度的增加涡流尺寸减少，当足够大的宽度是不存在涡流和回流。很明显，涡流现象的发展取决于雷诺数和几何形状的变化，也就是。因此，让我们通过流体迹线比较N-S方程和S方程。如图10所示，a b两图表明平流（惯性项）导致了涡流的早期发展，然而，c d两图中当深度增加两倍时可以看出，当不存在平流时也会出现涡流，没有平流的影响，涡流中心向右移动，因此，根据S方程的流体迹线关于x=0几乎对称。 图7 花键形凹坑的不同深度流体迹线图 图8 圆柱形凹坑的流体迹线图 图9 圆柱形凹坑的不同宽度流体迹线图 图10 依据N-S和S方程的圆柱形流体迹线图3.4 压力壁面力包括惯性力和分布力，可以从CFD求解器中直接求出。图5和图6可以看出静压累积随着雷诺数的增加单调增加，从图11可以得到证实，Fy是宽度的函数，表明压力随着宽度的增加而增加。比较圆柱形凹坑和花键形凹坑（xd=0）的压力可以看出当圆柱形凹坑深度为0.25时取得最大压力。图12（a）中花键形凹坑中压力与x
