CN113312729A - 一种飞溅润滑螺旋锥齿轮润滑状态的识别方法及系统 - Google Patents
一种飞溅润滑螺旋锥齿轮润滑状态的识别方法及系统 Download PDFInfo
- Publication number
- CN113312729A CN113312729A CN202110704929.2A CN202110704929A CN113312729A CN 113312729 A CN113312729 A CN 113312729A CN 202110704929 A CN202110704929 A CN 202110704929A CN 113312729 A CN113312729 A CN 113312729A
- Authority
- CN
- China
- Prior art keywords
- oil film
- bevel gear
- film deposition
- meshing
- determining
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000005461 lubrication Methods 0.000 title claims abstract description 64
- 238000000034 method Methods 0.000 title claims abstract description 53
- 230000001050 lubricating effect Effects 0.000 title claims description 17
- 230000008021 deposition Effects 0.000 claims abstract description 88
- 239000012530 fluid Substances 0.000 claims abstract description 60
- 238000004088 simulation Methods 0.000 claims abstract description 25
- 230000005540 biological transmission Effects 0.000 claims abstract description 24
- 238000007781 pre-processing Methods 0.000 claims abstract description 16
- 238000004364 calculation method Methods 0.000 claims description 37
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 239000003921 oil Substances 0.000 description 160
- 239000010687 lubricating oil Substances 0.000 description 11
- 238000010586 diagram Methods 0.000 description 3
- 238000005381 potential energy Methods 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 239000002131 composite material Substances 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- ZLSWBLPERHFHIS-UHFFFAOYSA-N Fenoprop Chemical compound OC(=O)C(C)OC1=CC(Cl)=C(Cl)C=C1Cl ZLSWBLPERHFHIS-UHFFFAOYSA-N 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000012208 gear oil Substances 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000007654 immersion Methods 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Computational Mathematics (AREA)
- Gears, Cams (AREA)
Abstract
本发明公开一种飞溅润滑螺旋锥齿轮润滑状态的识别方法及系统,利用三维建模软件建立锥齿轮传动系统的三维模型;对锥齿轮传动系统的三维模型进行预处理,获得网格文件;将网格文件导入已设置仿真参数的流体力学计算软件,计算多份油膜沉积厚度;基于CFD方法,根据多份油膜沉积厚度计算螺旋锥齿轮入口处齿面的平均油膜沉积厚度;根据螺旋锥齿轮确定入口位置处的临界油膜厚度;根据平均油膜沉积厚度和临界油膜厚度判断是否为充分润滑。本专利通过理论与仿真相结合,进而准确确定滑螺旋锥齿轮润滑状态为充分润滑状态还是乏油状态。
Description
技术领域
本发明涉及润滑状态识别技术领域,特别是涉及一种飞溅润滑螺旋锥齿轮润滑状态的识别方法及系统。
背景技术
航空螺旋锥齿轮高速重载、飞行姿态多变工况条件下,齿面润滑效果难以判定,若进入乏油润滑状态则严重威胁直升机的飞行安全。飞溅润滑方式下,机匣内油-气两相流动规律不明,滑油在齿面沉积的同时受离心力作用产生飞溅,齿面上形成微米级厚度的润滑油膜。因为通过实验测量微米级别的油膜厚度较为困难,所以目前无法实现针对螺旋锥齿轮的润滑状态进行准确判断。
发明内容
本发明的目的是提供一种飞溅润滑螺旋锥齿轮润滑状态的识别方法及系统,以实现准确判断润滑状态。
为实现上述目的,本发明提供了一种飞溅润滑螺旋锥齿轮润滑状态的识别方法,所述方法包括:
利用三维建模软件建立锥齿轮传动系统的三维模型;所述三维模型包括箱体、多个出油口、通气孔、主动锥齿轮和从动锥齿轮;
对锥齿轮传动系统的三维模型进行预处理,获得网格文件;
将所述网格文件导入已设置仿真参数的流体力学计算软件,计算多份油膜沉积厚度;
基于CFD方法,根据多份油膜沉积厚度计算螺旋锥齿轮入口处齿面的平均油膜沉积厚度;
根据螺旋锥齿轮确定入口位置处的临界油膜厚度;
判断所述平均油膜沉积厚度是否大于或等于所述临界油膜厚度;如果所述平均油膜沉积厚度大于或等于所述临界油膜厚度,则确定润滑状态为充分润滑状态;如果所述平均油膜沉积厚度小于所述临界油膜厚度,则确定润滑状态为乏油状态。
可选地,所述对锥齿轮传动系统的三维模型进行预处理,获得网格文件,具体包括:
对所述三维模型进行简化,获得简化模型;
通过布尔运算根据所述简化模型构建流体域;
根据所述流体域确定锥齿轮传动系统对应的流体数值计算模型;
将所述流体数值计算模型导入网格划分软件进行网格划分,获得网格文件。
可选地,所述将所述流体数值计算模型导入网格划分软件进行网格划分,获得网格文件,具体包括:
将所述流体数值计算模型导入网格划分软件,采用非结构化四面体网格,应用全局网格划分方法进行全局网格划分;
对已进行全局网格划分的关键部位应用局部网格划分方法进行局部网格加密,获得网格文件;所述关键部位包括齿面间距、齿轮壁面、通气孔和导油管。
可选地,所述将所述网格文件导入已设置仿真参数的流体力学计算软件,计算多份油膜沉积厚度,具体包括:
设置流体力学计算软件的仿真参数;
将所述网格文件导入已设置仿真参数的流体力学计算软件中,截取对应离散齿面的法向截面;
计算所述法向截面对应的油液体积分数等值线坐标和齿面坐标;
将所述油液体积分数等值线坐标和所述齿面坐标代入距离公式,计算沿齿宽方向第i份油膜宽度、沿齿高方向第j份油膜沉积长度对应的油膜沉积厚度hcij,其中,i和j均为大于等于1的正整数。
可选地,所述根据螺旋锥齿轮确定入口位置处的临界油膜厚度,具体包括:
将螺旋锥齿轮等效为当量斜齿轮;
计算当量斜齿轮的接触区半宽;
根据所述接触区半宽确定入口位置处的临界油膜厚度。
本发明还提供一种飞溅润滑螺旋锥齿轮润滑状态的识别系统,所述系统包括:
三维模型构建模块,用于利用三维建模软件建立锥齿轮传动系统的三维模型;所述三维模型包括箱体、多个出油口、通气孔、主动锥齿轮和从动锥齿轮;
预处理模块,用于对锥齿轮传动系统的三维模型进行预处理,获得网格文件;
导入模块,用于将所述网格文件导入已设置仿真参数的流体力学计算软件,计算多份油膜沉积厚度;
平均油膜沉积厚度确定模块,用于基于CFD方法,根据多份油膜沉积厚度计算螺旋锥齿轮入口处齿面的平均油膜沉积厚度;
临界油膜厚度确定模块,用于根据螺旋锥齿轮确定入口位置处的临界油膜厚度;
判断模块,用于判断所述平均油膜沉积厚度是否大于或等于所述临界油膜厚度;如果所述平均油膜沉积厚度大于或等于所述临界油膜厚度,则确定润滑状态为充分润滑状态;如果所述平均油膜沉积厚度小于所述临界油膜厚度,则确定润滑状态为乏油状态。
可选地,所述预处理模块,具体包括:
简化单元,用于对所述三维模型进行简化,获得简化模型;
流体域构建单元,用于通过布尔运算根据所述简化模型构建流体域;
流体数值计算模型确定单元,用于根据所述流体域确定锥齿轮传动系统对应的流体数值计算模型;
网格划分单元,用于将所述流体数值计算模型导入网格划分软件进行网格划分,获得网格文件。
可选地,所述网格划分单元,具体包括:
全局网格划分子单元,用于将所述流体数值计算模型导入网格划分软件,采用非结构化四面体网格,应用全局网格划分方法进行全局网格划分;
局部网格划分子单元,用于对已进行全局网格划分的关键部位应用局部网格划分方法进行局部网格加密,获得网格文件;所述关键部位包括齿面间距、齿轮壁面、通气孔和导油管。
可选地,所述油膜沉积厚度确定模块,具体包括:
仿真参数设置单元,用于设置流体力学计算软件的仿真参数;
法向截面截取单元,用于将所述网格文件导入已设置仿真参数的流体力学计算软件中,截取对应离散齿面的法向截面;
坐标计算单元,用于计算所述法向截面对应的油液体积分数等值线坐标和齿面坐标;
油膜沉积厚度确定单元,用于将所述油液体积分数等值线坐标和所述齿面坐标代入距离公式,计算沿齿宽方向第i份油膜宽度、沿齿高方向第j份油膜沉积长度对应的油膜沉积厚度hcij,其中,i和j均为大于等于1的正整数。
可选地,所述临界油膜厚度确定模块,具体包括:
齿轮转换单元,用于将螺旋锥齿轮等效为当量斜齿轮;
接触区半宽计算单元,用于计算当量斜齿轮的接触区半宽;
临界油膜厚度确定单元,用于根据所述接触区半宽确定入口位置处的临界油膜厚度。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明基于“飞溅润滑-沉积特性-入口油膜厚度-润滑状态识别”的研究思路,通过CFD方法计算螺旋锥齿轮齿面的平均油膜沉积厚度,并与基于弹流润滑理论计算的临界油膜厚度进行对比,据此判断螺旋锥齿轮是否处于充分润滑状态,以提高确定的准确性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明飞溅润滑螺旋锥齿轮润滑状态的识别方法流程图;
图2为本发明锥齿轮传动系统结构简图;
图3为本发明流体域示意图;
图4为本发明分离后的齿轮啮合间隙示意图;
图5为本发明网格划分示意图;
图6为本发明螺旋锥齿轮箱体初始时刻的状态示意图;
图7为本发明齿轮油液体积分布示意图;
图8为本发明油膜沉积厚度计算原理图;
图9为本发明线接触区域的油膜分布示意图;
图10为本发明飞溅润滑螺旋锥齿轮润滑状态的识别系统结构图;
符号说明:1-箱体,2-主动锥齿轮,3-从动锥齿轮,4-通气孔,5-出油口。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种飞溅润滑螺旋锥齿轮润滑状态的识别方法及系统,以实现准确判断润滑状态。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
实施例1
如图1所示,本发明公开一种飞溅润滑螺旋锥齿轮润滑状态的识别方法,所述方法包括:
步骤S1:利用三维建模软件建立锥齿轮传动系统的三维模型。
步骤S2:对锥齿轮传动系统的三维模型进行预处理,获得网格文件。
步骤S3:将所述网格文件导入已设置仿真参数的流体力学计算软件,计算多份油膜沉积厚度。
步骤S4:基于CFD方法,根据多份油膜沉积厚度计算螺旋锥齿轮入口处齿面的平均油膜沉积厚度。
步骤S5:根据螺旋锥齿轮确定入口位置处的临界油膜厚度。
步骤S6:判断所述平均油膜沉积厚度是否大于或等于所述临界油膜厚度;如果所述平均油膜沉积厚度大于或等于所述临界油膜厚度,则确定润滑状态为充分润滑状态;如果所述平均油膜沉积厚度小于所述临界油膜厚度,则确定润滑状态为乏油状态。
下面对各个步骤进行详细论述:
步骤S1:利用三维建模软件建立锥齿轮传动系统的三维模型。如图2所示,所述三维模型包括箱体1、多个出油口5、通气孔4、主动锥齿轮2和从动锥齿轮3。
步骤S2:对锥齿轮传动系统的三维模型进行预处理,获得网格文件,具体包括:
步骤S21:对所述三维模型进行简化,获得简化模型,具体包括:
步骤S211:去除所述三维模型中齿轮过渡圆角和倒圆角,获得平滑模型。
步骤S212:合并所述平滑模型中箱体壁面狭长细缝结构,获得简化模型。
步骤S22:通过布尔运算根据所述简化模型构建流体域,具体为:减去所述简化模型中的主动锥齿轮和从动锥齿轮,获得流体域,如图3所示。
步骤S23:根据所述流体域确定锥齿轮传动系统对应的流体数值计算模型。
步骤S24:将所述流体数值计算模型导入网格划分软件进行网格划分,获得网格文件,具体包括:
步骤S241:将所述流体数值计算模型导入网格划分软件,采用非结构化四面体网格,应用全局网格划分方法进行全局网格划分。
步骤S242:对已进行全局网格划分的关键部位应用局部网格划分方法进行局部网格加密,获得网格文件;所述关键部位包括齿面间距、齿轮壁面、通气孔和导油管;所述齿面间距为所述简化模型中主动锥齿轮和从动锥齿轮拉开的间距,齿面间隙小于1mm;所述网格文件为.msh格式。本发明将主动锥齿轮和从动锥齿轮拉开小段间距,使得啮合区域有小于1mm的齿面间距,具体详见图4。
所述网格文件中啮合区域内的齿面间距至少有两层网格;除了啮合区域以外的区域为计算域,计算域划分较稀疏的网格,如图5所示,图5中(a)为网格截面图,图5中(b)为导油管处局部细化网格图。控制计算域的网格划分质量,网格偏斜度(skewness)作为网格质量的判断标准;对主从动锥齿轮的轮齿、轮辐及各个壁面命名英文名称,如big gear_face、big_gear_lunfu、small_gear_face、small_gear_lunfu等。
步骤S3:将所述网格文件导入已设置仿真参数的流体力学计算软件,计算多份油膜沉积厚度,具体包括:
步骤S31:设置流体力学计算软件的仿真参数,具体包括:设置瞬态模式、多相流模型、湍流模型、润滑油的物性参数和壁面条件,编写描述齿轮的旋转运动状态的自定义函数(UserDefinedFunction:简称UDF),设置合适的求解方法和松弛迭代因子,初始化并设置初始工况下锥齿轮的浸油深度,设置时间步长;所述润滑油的物性参数包括密度、动力粘度。
具体的,进入FLUENT软件,在General中设置瞬态模式、开启重力加速度:沿Z轴负方向9.81m/s。
在Models中开启VOF多相流模型并将空气设为主相,润滑油设为次相、开启k-εRNG湍流模型。
在Materials中添加润滑油的物性参数,其中:密度为875.15kg/m3、动力粘度0.025125Pa/s。
在Boundary Conditions中设置壁面条件,按照在网格文件中命名的边界:出油口、通气孔设为压力出口(pressure-outlet)、其余壁面设为无滑移壁面(wall)。
在Dynamic Mesh开启动网格模型,选择弹簧光顺模型(Smoothing)和重构模型(Remeshing)。
编写用户自定义函数(UserDefinedFunction,UDF)描述齿轮的旋转运动状态并通过UserDefinedFunctions-compile导入FLUENT。
在Dynamic Mesh中选择Create/Edit,选择齿轮壁面和对应的角速度函数名称,并输入齿轮的旋转中心。
在Methods中设置求解方法,选用PISO压力-速度耦合求解格式、其余离散格式采用默认即可。
在Controls中设置亚松弛因子。
通过Settingup domain-Adapt-Region-选择两个坐标点-Mark,则所选区域为两坐标点所在水平面之间的区域。
通过Initialization-Initialize初始化并通过Patch选择上一步设置的Region,将Phase选为润滑油,选择Volume Fraction-Value为1,点击Patch,如图6所示。
在Run Calculation中设置时间步长,时间步长通过库朗数(Courant Number)进行估算,如下:
VOF模型中要求最大流动速度在一个时间步长内通过最小网格尺寸的1/4(即:库朗数默认为0.25),根据所需齿轮转数设置迭代时间步并点击Calculate开始计算。
计算完成后,通过Results-Graphics-Contours-Contours of-Phase-选择润滑油相-选中齿轮壁面-Save/Display,可以观察到齿轮的油液体积分数,如图7所示。
步骤S32:将所述网格文件导入已设置仿真参数的流体力学计算软件FLUENT中,截取对应离散齿面的法向截面。
步骤S33:计算所述法向截面对应的油液体积分数等值线坐标和齿面坐标。
步骤S34:将所述油液体积分数等值线坐标和所述齿面坐标代入距离公式,计算沿齿宽方向第i份油膜宽度、沿齿高方向第j份油膜沉积长度对应的油膜沉积厚度hcij。
如图8所示,将油膜沿齿宽n等分,每一份的油膜宽度为Δb=bi/n;计算出沿齿宽方向每一份的油膜沉积长度lj;再将该油膜沿齿高方向m等分,每一份的油膜沉积长度为Δl=lj/m;计算出油膜宽度为Δb、油膜沉积长度为Δl的油膜沉积厚度hcij,即可得到每一份的油膜体积ΔVij=Δb×Δl×hcij;将即将进入啮合区的齿面上沉积的每一份油膜体积求和并除以该齿面的面积A,最终得到齿面的平均油膜沉积厚度hc。因此总结的具体步骤如下:
步骤S4:基于CFD方法,根据多份油膜沉积厚度计算螺旋锥齿轮入口处齿面的平均油膜沉积厚度,具体包括:
步骤S41:基于CFD方法,根据多份油膜沉积厚度计算各份油膜沉积厚度对应的油膜体积ΔVij,具体公式为:
ΔVij=Δb×Δl×hcij(2);
其中,Δl为每一份的油膜沉积长度,Δb为每一份的油膜宽度,hcij为沿齿宽方向第i份油膜宽度、沿齿高方向第j份油膜沉积长度对应的油膜沉积厚度,ΔVij为油膜沉积厚度hcij对应的油膜体积,又称沿齿宽方向第i份油膜宽度、沿齿高方向第j份油膜沉积长度对应的油膜体积。
步骤S42:将进入啮合区的齿面上沉积的每一份油膜体积求和并除以齿面的面积A,最终得到螺旋锥齿轮入口处齿面的平均油膜沉积厚度hc,具体公式为:
其中,n为油膜沿齿宽等分的份数,m为油膜沿齿高方向等分的份数,ΔVij为沿齿宽方向第i份油膜宽度、沿齿高方向第j份油膜沉积长度对应的油膜体积,bi为齿宽方向的油膜宽度,lj为齿高方向的油膜沉积长度。
步骤S5:根据螺旋锥齿轮确定入口位置处的临界油膜厚度。
为获得弹流润滑分析时所需的齿轮力学参数,需进行螺旋锥齿轮的几何参数分析及啮合时的运动学分析。由于螺旋锥齿轮齿形复杂,齿轮参数受加工方式影响大,不易直接进行几何及运动学计算,因此,本发明需要将螺旋锥齿轮等效为当量斜齿轮,为求解螺旋锥齿轮入口临界油膜厚度提供力学参数。因此本发明总结步骤如下:
步骤S51:将螺旋锥齿轮等效为当量斜齿轮;螺旋锥齿轮等效为当量斜齿轮的转换关系如表1所示。
表1螺旋锥齿轮等效为当量斜齿轮的转换关系
表中:R为螺旋锥齿轮的外锥距,m;δ1、δ2分别为主动锥齿轮和从动锥齿轮的分锥角,(°);β为螺旋角;z1、z2分别为螺旋锥齿轮中的主动锥齿轮和从动锥齿轮的齿数;ω1、ω2分别为螺旋锥齿轮中的主动锥齿轮和从动锥齿轮的角速度;ra、rb为螺旋锥齿轮中的主动锥齿轮和从动锥齿轮的节圆半径,ra=Rsinδ1,rb=Rsinδ2;z1'、z2'分别为当量斜齿轮中的主动锥齿轮和从动锥齿轮的齿数,z1'=z1/cosδ1,z2'=z1/cosδ2;ra'、rb'分别为当量斜齿轮中的主动锥齿轮和从动锥齿轮的节圆半径,ra'=Rtanδ1,rb'=Rtanδ2;ω1'、ω2'分别为当量斜齿轮中的主动锥齿轮和从动锥齿轮的角速度,ω1'=ω1cosδ1,ω2'=ω1sinδ1/tanδ2。
步骤S52:计算当量斜齿轮的接触区半宽,具体包括:
计算单齿面法向载荷分布主要有:接触线百分比法及最小势能法。接触线百分比法的基本原理为计算同时啮合各齿对对应的时变接触线长度ln占总时变接触线长度L的百分比;最小势能法是基于最小势能原理(minimum elastic potential energy:简称MEPE),得到考虑载荷在时间和空间上变化规律的载荷分布。斜齿轮的螺旋角逐渐增加后,两种单齿面法向载荷分布的计算方法差别较小。因此本发明采用接触线百分比法对单位接触线上载荷进行计算,具体步骤如下:
步骤S522:采用接触线百分比法,根据所述总时变接触线长度计算单位接触线上载荷,具体公式为:
其中,w为单位接触线上载荷,N,ln为同时啮合第n个齿对对应的时变接触线长度,L为当量斜齿轮的总时变接触线长度,Fn为当量斜齿轮的齿面法向载荷。
本发明当量斜齿轮的齿面法向载荷应与螺旋锥齿轮齿宽中点处的齿面法向载荷相等,因此当量斜齿轮的齿面法向载荷的具体计算公式为:
其中,Tin为输入扭矩,Nm;dm1为主动轮齿宽中点处的分度圆直径,m;α为压力角,β为螺旋角。
在单位接触线上载荷w作用下,弹性圆柱体接触时产生变形,接触线形成一个细长的面,其接触区半宽为bh;接触区应力呈现椭圆分布,最大接触应力为ph。由弹流润滑理论可知,接触点的综合曲率半径远大于油膜厚度,因此将接触区视为平面。根据Hertz弹性接触理论,接触区半宽bh及最大接触应力ph的计算方法为式(7)和式(9),且最大接触应力与接触区半宽bh为负相关。
步骤S523:根据所述单位接触线上载荷计算接触区半宽,具体公式为:
其中,bh为接触区半宽,w为单位接触线上载荷,Rc为综合曲率半径,m,E'为当量弹性模量,Pa。
计算当量弹性模量E'的具体公式为:
其中,E1、E2分别为主动锥齿轮及从动锥齿轮的弹性模量,Pa;μ1、μ2均为齿轮泊松比。
步骤S524:根据所述接触区半宽和所述单位接触线上载荷计算最大接触应力,具体公式为:
其中,ph为最大接触应力,bh为接触区半宽,w为单位接触线上载荷。
将螺旋锥齿轮等效为当量斜齿轮,对于斜齿轮来说,其啮合点处椭圆率通常大于8,因此可近似采用线接触计算方法分析齿轮的膜厚情况。根据Castele和Dowson理论,对于线接触问题,若无量纲入口位置坐标X=|xa|/bh≥4时,则为充分润滑,可以理解为:入口位置坐标X为接触区半宽的4倍时,即有充分的滑油可以到达距离啮合点的4倍接触区半宽,这是齿轮充分润滑与乏油润滑的临界位置xa,如图9所示。若入口位置坐标位于临界位置xa的左侧,说明齿轮啮合区处于充分润滑状态,否则,齿轮啮合区处于乏油状态,润滑不充分。因此,充分润滑临界点xa=4bh时,入口位置处的临界油膜厚度hin为表征齿轮啮合区充分润滑的临界油膜厚度。
步骤S53:根据所述接触区半宽确定入口位置处的临界油膜厚度,具体公式为:
其中,hin(t)为入口位置处的临界油膜厚度,hmin(t)为最小油膜厚度,Rc(t)为综合曲率半径,m;bh(t)为接触区半宽,xa(t)为齿轮即将进入啮合的入口区坐标值,Rc为Rc(t)的简写,bh为bh(t)的简写,hmin为hmin(t)的简写,hin为hin(t)的简写,xa为xa(t)的简写。
由于油膜出口区存在颈缩现象,所以油膜出口区的膜厚为最小油膜厚度,因此计算最小油膜厚度的具体公式为:
hmin(t)=2.65κ0.54[η0·ur(t)]0.7E′-0.03[Rc(t)]0.43[w(t)]-0.13(11);
其中,hmin(t)为最小油膜厚度,κ为润滑油的黏压系数,m2/N;η0为润滑油初始动力粘度,Pa/s;ur(t)为啮合点处的卷吸速度,m/s;E'为当量弹性模量;Rc(t)为综合曲率半径,m;w(t)为单位接触线载荷,N。
实施例2
如图10所示,本发明公开一种飞溅润滑螺旋锥齿轮润滑状态的识别系统,所述系统包括:
三维模型构建模块100,用于利用三维建模软件建立锥齿轮传动系统的三维模型;所述三维模型包括箱体、多个出油口、通气孔、主动锥齿轮和从动锥齿轮。
预处理模块200,用于对锥齿轮传动系统的三维模型进行预处理,获得网格文件。
导入模块300,用于将所述网格文件导入已设置仿真参数的流体力学计算软件,计算多份油膜沉积厚度。
平均油膜沉积厚度确定模块400,用于基于CFD方法,根据多份油膜沉积厚度计算螺旋锥齿轮入口处齿面的平均油膜沉积厚度。
临界油膜厚度确定模块500,用于根据螺旋锥齿轮确定入口位置处的临界油膜厚度。
判断模块600,用于判断所述平均油膜沉积厚度是否大于或等于所述临界油膜厚度;如果所述平均油膜沉积厚度大于或等于所述临界油膜厚度,则确定润滑状态为充分润滑状态;如果所述平均油膜沉积厚度小于所述临界油膜厚度,则确定润滑状态为乏油状态。
下面对各个模块进行详细论述:
作为一种可选的实施方式,本发明所述预处理模块200,具体包括:
简化单元,用于对所述三维模型进行简化,获得简化模型。
流体域构建单元,用于通过布尔运算根据所述简化模型构建流体域。
流体数值计算模型确定单元,用于根据所述流体域确定锥齿轮传动系统对应的流体数值计算模型。
网格划分单元,用于将所述流体数值计算模型导入网格划分软件进行网格划分,获得网格文件。
作为一种可选的实施方式,本发明所述网格划分单元,具体包括:
全局网格划分子单元,用于将所述流体数值计算模型导入网格划分软件,采用非结构化四面体网格,应用全局网格划分方法进行全局网格划分。
局部网格划分子单元,用于对已进行全局网格划分的关键部位应用局部网格划分方法进行局部网格加密,获得网格文件;所述关键部位包括齿面间距、齿轮壁面、通气孔和导油管。
作为一种可选的实施方式,本发明所述油膜沉积厚度确定模块400,具体包括:
仿真参数设置单元,用于设置流体力学计算软件的仿真参数。
法向截面截取单元,用于将所述网格文件导入已设置仿真参数的流体力学计算软件中,截取对应离散齿面的法向截面。
坐标计算单元,用于计算所述法向截面对应的油液体积分数等值线坐标和齿面坐标。
油膜沉积厚度确定单元,用于将所述油液体积分数等值线坐标和所述齿面坐标代入距离公式,计算沿齿宽方向第i份油膜宽度、沿齿高方向第j份油膜沉积长度对应的油膜沉积厚度hcij,其中,i和j均为大于等于1的正整数。
作为一种可选的实施方式,本发明所述临界油膜厚度确定模块500,具体包括:
齿轮转换单元,用于将螺旋锥齿轮等效为当量斜齿轮。
接触区半宽计算单元,用于计算当量斜齿轮的接触区半宽。
临界油膜厚度确定单元,用于根据所述接触区半宽确定入口位置处的临界油膜厚度。
与实施例1相同的部分不再一一论述,具体详见实施例1。
本发明公开的技术方案存在以下优点:
(1)、本发明基于“飞溅润滑-沉积特性-入口油膜厚度-润滑状态识别”的研究思路,通过CFD方法计算螺旋锥齿轮齿面油膜沉积厚度,并与基于弹流润滑理论计算的临界油膜厚度进行对比,据此判断螺旋锥齿轮是否处于充分润滑状态,以提高确定的准确性。
(2)本发明判断螺旋锥齿轮齿面润滑效果提供明确的技术途径,建立齿面宏观流场与微观油膜厚度间的映射关系,为完善螺旋锥齿轮润滑参数的设计提供优化方法。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (10)
1.一种飞溅润滑螺旋锥齿轮润滑状态的识别方法,其特征在于,所述方法包括:
利用三维建模软件建立锥齿轮传动系统的三维模型;所述三维模型包括箱体、多个出油口、通气孔、主动锥齿轮和从动锥齿轮;
对锥齿轮传动系统的三维模型进行预处理,获得网格文件;
将所述网格文件导入已设置仿真参数的流体力学计算软件,计算多份油膜沉积厚度;
基于CFD方法,根据多份油膜沉积厚度计算螺旋锥齿轮入口处齿面的平均油膜沉积厚度;
根据螺旋锥齿轮确定入口位置处的临界油膜厚度;
判断所述平均油膜沉积厚度是否大于或等于所述临界油膜厚度;如果所述平均油膜沉积厚度大于或等于所述临界油膜厚度,则确定润滑状态为充分润滑状态;如果所述平均油膜沉积厚度小于所述临界油膜厚度,则确定润滑状态为乏油状态。
2.根据权利要求1所述的飞溅润滑螺旋锥齿轮润滑状态的识别方法,其特征在于,所述对锥齿轮传动系统的三维模型进行预处理,获得网格文件,具体包括:
对所述三维模型进行简化,获得简化模型;
通过布尔运算根据所述简化模型构建流体域;
根据所述流体域确定锥齿轮传动系统对应的流体数值计算模型;
将所述流体数值计算模型导入网格划分软件进行网格划分,获得网格文件。
3.根据权利要求2所述的飞溅润滑螺旋锥齿轮润滑状态的识别方法,其特征在于,所述将所述流体数值计算模型导入网格划分软件进行网格划分,获得网格文件,具体包括:
将所述流体数值计算模型导入网格划分软件,采用非结构化四面体网格,应用全局网格划分方法进行全局网格划分;
对已进行全局网格划分的关键部位应用局部网格划分方法进行局部网格加密,获得网格文件;所述关键部位包括齿面间距、齿轮壁面、通气孔和导油管。
4.根据权利要求1所述的飞溅润滑螺旋锥齿轮润滑状态的识别方法,其特征在于,所述将所述网格文件导入已设置仿真参数的流体力学计算软件,计算多份油膜沉积厚度,具体包括:
设置流体力学计算软件的仿真参数;
将所述网格文件导入已设置仿真参数的流体力学计算软件中,截取对应离散齿面的法向截面;
计算所述法向截面对应的油液体积分数等值线坐标和齿面坐标;
将所述油液体积分数等值线坐标和所述齿面坐标代入距离公式,计算沿齿宽方向第i份油膜宽度、沿齿高方向第j份油膜沉积长度对应的油膜沉积厚度hcij,其中,i和j均为大于等于1的正整数。
5.根据权利要求1所述的飞溅润滑螺旋锥齿轮润滑状态的识别方法,其特征在于,所述根据螺旋锥齿轮确定入口位置处的临界油膜厚度,具体包括:
将螺旋锥齿轮等效为当量斜齿轮;
计算当量斜齿轮的接触区半宽;
根据所述接触区半宽确定入口位置处的临界油膜厚度。
6.一种飞溅润滑螺旋锥齿轮润滑状态的识别系统,其特征在于,所述系统包括:
三维模型构建模块,用于利用三维建模软件建立锥齿轮传动系统的三维模型;所述三维模型包括箱体、多个出油口、通气孔、主动锥齿轮和从动锥齿轮;
预处理模块,用于对锥齿轮传动系统的三维模型进行预处理,获得网格文件;
导入模块,用于将所述网格文件导入已设置仿真参数的流体力学计算软件,计算多份油膜沉积厚度;
平均油膜沉积厚度确定模块,用于基于CFD方法,根据多份油膜沉积厚度计算螺旋锥齿轮入口处齿面的平均油膜沉积厚度;
临界油膜厚度确定模块,用于根据螺旋锥齿轮确定入口位置处的临界油膜厚度;
判断模块,用于判断所述平均油膜沉积厚度是否大于或等于所述临界油膜厚度;如果所述平均油膜沉积厚度大于或等于所述临界油膜厚度,则确定润滑状态为充分润滑状态;如果所述平均油膜沉积厚度小于所述临界油膜厚度,则确定润滑状态为乏油状态。
7.根据权利要求6所述的飞溅润滑螺旋锥齿轮润滑状态的识别系统,其特征在于,所述预处理模块,具体包括:
简化单元,用于对所述三维模型进行简化,获得简化模型;
流体域构建单元,用于通过布尔运算根据所述简化模型构建流体域;
流体数值计算模型确定单元,用于根据所述流体域确定锥齿轮传动系统对应的流体数值计算模型;
网格划分单元,用于将所述流体数值计算模型导入网格划分软件进行网格划分,获得网格文件。
8.根据权利要求7所述的飞溅润滑螺旋锥齿轮润滑状态的识别系统,其特征在于,所述网格划分单元,具体包括:
全局网格划分子单元,用于将所述流体数值计算模型导入网格划分软件,采用非结构化四面体网格,应用全局网格划分方法进行全局网格划分;
局部网格划分子单元,用于对已进行全局网格划分的关键部位应用局部网格划分方法进行局部网格加密,获得网格文件;所述关键部位包括齿面间距、齿轮壁面、通气孔和导油管。
9.根据权利要求6所述的飞溅润滑螺旋锥齿轮润滑状态的识别系统,其特征在于,所述油膜沉积厚度确定模块,具体包括:
仿真参数设置单元,用于设置流体力学计算软件的仿真参数;
法向截面截取单元,用于将所述网格文件导入已设置仿真参数的流体力学计算软件中,截取对应离散齿面的法向截面;
坐标计算单元,用于计算所述法向截面对应的油液体积分数等值线坐标和齿面坐标;
油膜沉积厚度确定单元,用于将所述油液体积分数等值线坐标和所述齿面坐标代入距离公式,计算沿齿宽方向第i份油膜宽度、沿齿高方向第j份油膜沉积长度对应的油膜沉积厚度hcij,其中,i和j均为大于等于1的正整数。
10.根据权利要求6所述的飞溅润滑螺旋锥齿轮润滑状态的识别系统,其特征在于,所述临界油膜厚度确定模块,具体包括:
齿轮转换单元,用于将螺旋锥齿轮等效为当量斜齿轮;
接触区半宽计算单元,用于计算当量斜齿轮的接触区半宽;
临界油膜厚度确定单元,用于根据所述接触区半宽确定入口位置处的临界油膜厚度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110704929.2A CN113312729A (zh) | 2021-06-24 | 2021-06-24 | 一种飞溅润滑螺旋锥齿轮润滑状态的识别方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110704929.2A CN113312729A (zh) | 2021-06-24 | 2021-06-24 | 一种飞溅润滑螺旋锥齿轮润滑状态的识别方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113312729A true CN113312729A (zh) | 2021-08-27 |
Family
ID=77380347
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110704929.2A Pending CN113312729A (zh) | 2021-06-24 | 2021-06-24 | 一种飞溅润滑螺旋锥齿轮润滑状态的识别方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113312729A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114117673A (zh) * | 2021-11-25 | 2022-03-01 | 重庆大学 | 一种高线速度齿轮最优润滑喷嘴方位确定方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111985062A (zh) * | 2020-08-24 | 2020-11-24 | 哈尔滨工程大学 | 一种考虑三维表面加工粗糙度的柴油机正时齿轮润滑状态预测方法 |
CN112100938A (zh) * | 2020-09-11 | 2020-12-18 | 南京航空航天大学 | 一种飞溅润滑锥齿轮温度的数值计算方法及系统 |
-
2021
- 2021-06-24 CN CN202110704929.2A patent/CN113312729A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111985062A (zh) * | 2020-08-24 | 2020-11-24 | 哈尔滨工程大学 | 一种考虑三维表面加工粗糙度的柴油机正时齿轮润滑状态预测方法 |
CN112100938A (zh) * | 2020-09-11 | 2020-12-18 | 南京航空航天大学 | 一种飞溅润滑锥齿轮温度的数值计算方法及系统 |
Non-Patent Citations (6)
Title |
---|
叶友东等: ""弹流润滑条件下齿轮传动润滑油膜厚度影响参数研究"", 《机械设计》, vol. 29, no. 9, 30 September 2012 (2012-09-30), pages 73 - 76 * |
尹明虎: ""考虑润滑效应的人字齿轮-转子-滑动轴承系统动力学特性分析"", 《中国博士学位论文全文数据库 工程科技II辑》, 15 January 2019 (2019-01-15), pages 39 - 63 * |
曾韬等: "《螺旋齿轮设计与加工》", 31 May 1989, 哈尔滨工业大学出版社, pages: 48 - 49 * |
李玉哲等: ""基于CFD的某减速器飞溅润滑 仿真分析与验证"", 《航空动力学报》, vol. 35, no. 7, 31 July 2020 (2020-07-31), pages 1482 - 1488 * |
陆凤霞等: ""直升机中减飞溅润滑流场分析与优化方法"", 《航空学报》, vol. 41, no. 11, 13 April 2020 (2020-04-13), pages 123659 - 2 * |
陆凤霞等: ""直升机中减飞溅润滑流场分析与优化方法"", 《航空学报》, vol. 41, no. 11, 30 April 2020 (2020-04-30), pages 123659 - 1 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114117673A (zh) * | 2021-11-25 | 2022-03-01 | 重庆大学 | 一种高线速度齿轮最优润滑喷嘴方位确定方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ooi et al. | Reynolds averaged simulation of flow and heat transfer in ribbed ducts | |
Lositaño et al. | Steady wind performance of a 5 kW three-bladed H-rotor Darrieus Vertical Axis Wind Turbine (VAWT) with cambered tubercle leading edge (TLE) blades | |
CN111460699B (zh) | 平壁表面减阻功能微织构的设计方法 | |
Deng et al. | Lubrication mechanism in gearbox of high-speed railway trains | |
CN101567026B (zh) | 一种输液管道导流结构减振弯头优化设计方法 | |
CN113312729A (zh) | 一种飞溅润滑螺旋锥齿轮润滑状态的识别方法及系统 | |
Webb et al. | Parametric modelling of a spiral bevel gear using CFD | |
CN111079226B (zh) | 一种基于cfd的汽车变速器润滑油液位仿真方法 | |
Saha et al. | On scaling pipe flows with sinusoidal transversely corrugated walls: analysis of data from the laminar to the low-Reynolds-number turbulent regime | |
Larsen | Impact loads on circular cylinders | |
Jain et al. | Capturing the dynamic stall in H-type Darrieus wind turbines using different URANS turbulence models | |
Ameri | Implicit-les simulation of variable-speed power turbine cascade for low free-stream turbulence conditions | |
Dahlström et al. | Large eddy simulation applied to a high-Reynolds flow around an airfoil close to stall | |
Prabhakar et al. | Numerical investigations to assess the impact of shaft speed on the performance of scoop devices | |
Betancour et al. | Design and optimization of a runner for a gravitational vortex turbine using the response surface methodology and experimental tests | |
Lu et al. | Oil film deposition characteristics and judgment of lubrication effect of splash lubricated gears | |
Mubarok et al. | Numerical and analytical modeling of pressure drop through a geothermal two-phase orifice plate | |
CN114626266A (zh) | 轻度失油状态下直升机中减摩擦动力学确定方法及系统 | |
Yuan et al. | Numerical simulation of bubble motion in horizontal reducer pipelines | |
Sunny | Effect of turbulence in modeling the reduction of local drag forces in a computational automotive model | |
Yang et al. | Simulation and numerical calculation on pipeline leakage process | |
Li et al. | Effects of discharge-to-submergence ratio on evolution of air-core vortex | |
Yin et al. | Assessment of RANS to predict flows with large streamline curvature | |
Allan | A CFD investigation of wind tunnel interference on delta wing aerodynamics | |
Liu et al. | Influence of vibration on the lubrication effect of a splash-lubricated gearbox |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |