CN115310218B 点接触混合弹流润滑油膜厚度的确定方法 （武汉大学）.docxVIP

(19)国家知识产权局

(12)发明专利

(10)授权公告号CN115310218B(45)授权公告日2025.07.11

(21)申请号202210814003.3

(22)申请日2022.07.11

(65)同一申请的已公布的文献号申请公布号CN115310218A

(43)申请公布日2022.11.08

(73)专利权人武汉大学

地址430072湖北省武汉市武昌区珞珈山

街道八一路299号

(72)发明人孙韵韵余欣李鸿光

(74)专利代理机构武汉科皓知识产权代理事务所(特殊普通合伙)42222

专利代理师罗敏清

(51)Int.CI.

GO6F30/17(2020.01)

GO6F30/23(2020.01)

G06F119/14(2020.01)

(56)对比文件

CN108984933A,2018.12.11

CN110119576A,2019.08.13

审查员张敏

权利要求书3页说明书8页附图2页

(54)发明名称

点接触混合弹流润滑油膜厚度的确定方法

(57)摘要

CN115310218B本发明提供一种点接触混合弹流润滑油膜厚度的确定方法，包括：测量形成点接触的第一物体与第二物体的三维轮廓高度，依据轮廓高度计算油膜初始厚度h?,并根据外加载荷F设定油膜的初始压力分布p。(x,y);根据h?和p。(x,y)推导润滑油膜厚度方程，再引入表征油膜厚度与压力之间关系的Reynolds方程，对膜厚方程和Reynolds方程进行无量纲化在进行离散；根据步骤2计算接触界面每个节点的油膜厚度h和压力p,从而计算固体部分的真实接触面积；根据得到的真实接触面积，计算固体接触力F和混合润滑界面总接触力F;判定计算载荷F和外加载荷F之间是否满足收敛条件，如不满足，则进行更新

CN115310218B

开始

开始

输入最荷尺、速度以突微连形能D、G

设定初始条件：油膜厚度肠，压力分布p(J)

粘压和留压方糕

油膜厚度分布Mxy)润滑油粘度》、密度pReynolda方

油膜压力分布px)

遍历计算节点，计算微凸体接触比例2分利理诊

P(xJ)PK)

计算相总峰最大接狱面和

计算液体鞍力F计昇周体接力F,

计算总裁荷F=F?+F.

收数性判断

输出油膜厚度分布(xy)

结来

更新初始油膜厚度：

ho=fF)

油膜压力分布：

CN115310218B权利要求书1/3页

2

1.一种点接触混合弹流润滑油膜厚度的确定方法，其特征在于，包括步骤如下：

S1:测量形成点接触的第一物体与第二物体的三维轮廓高度，其中，第一物体为具有曲面的物体，根据得到的三维轮廓高度建立粗糙表面的轮廓高度函数，再依据轮廓高度函数计算油膜初始厚度h?,并根据外加载荷F设定油膜的初始压力分布p。(x,y);

S2:计算初始压力下的润滑油粘度η和密度p,并根据油膜初始厚度h?和初始压力分布p?(x,y)推导润滑油膜厚度方程，再根据计算的润滑油粘度η、密度p以及油膜厚度方程引入表征油膜厚度与压力之间关系的Reynolds方程，对膜厚方程和Reynolds方程进行无量纲化，并通过有限差分法离散化无量纲化后的方程；

S3:根据步骤S2的膜厚方程计算接触界面每个节点的油膜厚度h,并根据步骤S2离散化后的无量纲化方程计算接触界面每个节点的压力p,得到混合弹流润滑油膜厚度分布h(x,y)和压力分布p(x,y),并根据压力分布p(x,y)计算固体接触节点占总节点数的比例λ,从而得到固体部分的真实接触面积；

S4:根据步骤S3得到的真实接触面积，计算单个粗糙峰的最大接触面积a,再结合最大接触面积a计算固体接触力F,根据固体接触力和液体接触力得到混合润滑界面总接触力Fn;

S5:判定计算载荷F和外加载荷F之间是否满足收敛条件，如不满足，则根据计算载荷F修正油膜初始厚度h?和油膜压力分布p。(x,y),迭代计算直到载荷收敛，最终得到点接触混合润滑界面的油膜厚度分布h(x,y);

其中，步骤S1具体包括：

采用测量仪器获得第一物体和第二物体的三维轮廓高度，建立接触界面的自相关函数，对自相关函数进行傅氏变换得到真实表面轮廓的功率谱，建立功率谱函数与空间频率