CN116415462B - 基于浮动衬套的双层油膜润滑分析方法及系统 - Google Patents
基于浮动衬套的双层油膜润滑分析方法及系统 Download PDFInfo
- Publication number
- CN116415462B CN116415462B CN202310404676.6A CN202310404676A CN116415462B CN 116415462 B CN116415462 B CN 116415462B CN 202310404676 A CN202310404676 A CN 202310404676A CN 116415462 B CN116415462 B CN 116415462B
- Authority
- CN
- China
- Prior art keywords
- oil film
- roller
- floating
- determining
- displacement
- 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.)
- Active
Links
- 238000005461 lubrication Methods 0.000 title claims abstract description 68
- 238000004458 analytical method Methods 0.000 title claims abstract description 37
- 238000006073 displacement reaction Methods 0.000 claims abstract description 82
- 238000000034 method Methods 0.000 claims abstract description 27
- 230000002159 abnormal effect Effects 0.000 claims abstract description 14
- 230000008569 process Effects 0.000 claims abstract description 14
- 239000010410 layer Substances 0.000 claims description 149
- 230000005489 elastic deformation Effects 0.000 claims description 26
- 238000004590 computer program Methods 0.000 claims description 9
- 239000002355 dual-layer Substances 0.000 claims description 5
- 239000002131 composite material Substances 0.000 claims description 4
- 230000005856 abnormality Effects 0.000 abstract description 6
- 239000003921 oil Substances 0.000 description 193
- 238000004364 calculation method Methods 0.000 description 12
- 230000007246 mechanism Effects 0.000 description 8
- 230000001052 transient effect Effects 0.000 description 7
- 238000004891 communication Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 239000010687 lubricating oil Substances 0.000 description 4
- 230000001133 acceleration Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 125000004122 cyclic group Chemical group 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000005299 abrasion Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 239000010710 diesel engine oil Substances 0.000 description 1
- 230000020169 heat generation Effects 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 230000009191 jumping Effects 0.000 description 1
- 230000001050 lubricating effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002356 single layer Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- 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/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- 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/08—Thermal analysis or thermal optimisation
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Operations Research (AREA)
- Evolutionary Computation (AREA)
- Algebra (AREA)
- Geometry (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Sliding-Contact Bearings (AREA)
Abstract
本发明提供一种基于浮动衬套的双层油膜润滑分析方法及系统，属于柴油机技术领域，该方法包括：在凸轮转动过程中，确定每一时刻下浮动衬套内层油膜厚度、外层油膜厚度、滚轮位移以及浮动衬套位移；根据滚轮位移和浮动衬套位移确定滚轮的运动轨迹和浮动衬套的运动轨迹；当滚轮的运动轨迹和浮动衬套的运动轨迹均首尾闭合时，根据每一时刻下浮动衬套内层油膜厚度、外层油膜厚度确定浮动衬套内层油膜厚度最小值、外层油膜厚度最小值；当浮动衬套内层油膜厚度最小值、外层油膜厚度最小值中任一个小于预设油膜厚度阈值时，确定双层油膜润滑存在异常。本发明能够实现更准确的双层油膜润滑异常分析。
Description
技术领域
本发明属于柴油机技术领域，尤其涉及一种基于浮动衬套的双层油膜润滑分析方法及系统。
背景技术
滚轮-浮动衬套-滚轮销的组合结构作为柴油机供油机构中的关键部分，其润滑性能会直接影响柴油机的可靠性和寿命。在工作过程中，凸轮带动滚轮转动，滚轮销为固定件，而衬套是浮动的，浮动衬套内外可形成双层承载油膜，承载能力较强；同时浮动衬套具有一定的转速，可降低滚轮与滚轮销之间的相对运动速度，从而降低摩擦功耗与发热。因此，在传统柴油机供油机构较为典型的低速重载工况下，相较于单层油膜润滑的滚轮-滚轮销结构，滚轮-浮动衬套-滚轮销结构具有承载能力强、摩擦温升小等优势。不过，在实际工作中，浮动衬套常会出现异常抱死、磨损等润滑故障，影响供油机构可靠性；但目前针对供油机构滚轮-浮动衬套-滚轮销的润滑状态识别往往为稳定工况，因此需要建立更为完善的滚轮-浮动衬套-滚轮销瞬态润滑模型，将凸轮-滚轮运动学、动力学关系与滚轮-浮动衬套-滚轮销瞬态润滑相耦合，并基于此实现更为准确的滚轮-浮动衬套-滚轮销的双层油膜润滑分析。
发明内容
针对现有技术存在的问题，本发明提供一种基于浮动衬套的双层油膜润滑分析方法及系统，能够对浮动衬套的双层油膜润滑状态进行更准确的分析判断。
第一方面，本发明提供一种基于浮动衬套的双层油膜润滑分析方法，包括：
在凸轮转动过程中，确定每一时刻下浮动衬套内层油膜厚度和外层油膜厚度，以及确定每一时刻下滚轮位移和浮动衬套位移；
根据所述每一时刻下滚轮位移和浮动衬套位移分别确定滚轮的运动轨迹和浮动衬套的运动轨迹；
当所述滚轮的运动轨迹和浮动衬套的运动轨迹均首尾闭合时，根据所述每一时刻下浮动衬套内层油膜厚度和外层油膜厚度分别确定浮动衬套内层油膜厚度最小值和外层油膜厚度最小值；
当浮动衬套内层油膜厚度最小值和外层油膜厚度最小值中任一个小于预设油膜厚度阈值时，确定双层油膜润滑存在异常。
可选地，获取第一组合结构的结构参数，所述第一组合结构包括凸轮和滚轮；
根据所述第一组合结构的结构参数计算滚轮时变转速；
获取第二组合结构的结构参数，所述第二组合结构包括滚轮、浮动衬套和滚轮销；
根据所述第二组合结构的结构参数计算浮动衬套当前时刻的内层油膜厚度和外层油膜厚度；
根据所述滚轮时变转速、内层油膜厚度和外层油膜厚度求解内外双层油膜平均雷诺方程，当压力收敛时，得到当前时刻下浮动衬套内层油膜压力和外层油膜压力；
根据所述内层油膜压力和外层油膜压力分别确定浮动衬套内层弹性变形和外层弹性变形；
将所述内层弹性变形和外层弹性变形代入油膜厚度方程，当弹性变形收敛、浮动衬套内外层达到热平衡条件时，输出所述当前时刻下浮动衬套内层油膜压力和外层油膜压力；
确定当前时刻的滚轮位移与浮动衬套位移；
根据所述滚轮位移与浮动衬套位移确定当前时刻下滚轮运动轨迹和浮动衬套运动轨迹；
若滚轮运动轨迹和浮动衬套运动轨迹二者中任一个没有首尾闭合，则获取下一时刻的第二组合结构的结构参数；
根据所述下一时刻的第二组合结构的结构参数确定下一时刻的浮动衬套内层油膜厚度和外层油膜厚度。
可选地，根据所述第一组合结构的结构参数计算滚轮时变转速，具体包括：
根据凸轮型线参数计算凸轮压力角和凸轮曲率半径；
根据所述凸轮压力角和凸轮曲率半径计算凸轮表面速度；
根据所述凸轮表面速度计算滚轮表面速度；
根据所述滚轮表面速度计算滚轮时变转速。
可选地，确定每一时刻下滚轮位移和浮动衬套位移，具体包括：
获取第一组合结构的结构参数，所述第一组合结构包括凸轮和滚轮；
根据所述第一组合结构的结构参数计算滚轮时变载荷；
根据所述滚轮时变载荷、滚轮与浮动衬套之间的运动学方程以及浮动衬套转速方程确定每一时刻下滚轮位移和浮动衬套位移。
可选地，所述滚轮时变载荷包括弹簧对滚轮的作用力、惯性力和柱塞油腔液压力。
可选地，所述预设油膜厚度阈值为0.8微米到1.2微米范围内的任意值。
第二方面，本发明提供一种基于浮动衬套的双层油膜润滑分析系统，包括：
第一参数确定模块，用于在凸轮转动过程中，确定每一时刻下浮动衬套内层油膜厚度和外层油膜厚度，以及确定每一时刻下滚轮位移和浮动衬套位移；
运动轨迹确定模块，用于根据所述每一时刻下滚轮位移和浮动衬套位移分别确定滚轮的运动轨迹和浮动衬套的运动轨迹；
第二参数确定模块，用于当所述滚轮的运动轨迹和浮动衬套的运动轨迹均首尾闭合时，根据所述每一时刻下浮动衬套内层油膜厚度和外层油膜厚度分别确定浮动衬套内层油膜厚度最小值和外层油膜厚度最小值；
双层油膜润滑分析模块，用于当浮动衬套内层油膜厚度最小值和外层油膜厚度最小值中任一个小于预设油膜厚度阈值时，确定双层油膜润滑存在异常。
第三方面，本发明提供一种电子设备，包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述程序时实现所述的基于浮动衬套的双层油膜润滑分析方法。
第四方面，本发明提供一种非暂态计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现所述的基于浮动衬套的双层油膜润滑分析方法。
本发明在凸轮转动过程中，确定每一时刻下浮动衬套内层油膜厚度和外层油膜厚度，以及确定每一时刻下滚轮位移和浮动衬套位移；根据每一时刻下滚轮位移和浮动衬套位移分别确定滚轮的运动轨迹和浮动衬套的运动轨迹；当滚轮的运动轨迹和浮动衬套的运动轨迹均首尾闭合时，根据每一时刻下浮动衬套内层油膜厚度和外层油膜厚度分别确定浮动衬套内层油膜厚度最小值和外层油膜厚度最小值；当浮动衬套内层油膜厚度最小值和外层油膜厚度最小值中任一个小于预设油膜厚度阈值时，确定双层油膜润滑存在异常。即本发明利用了凸轮和滚轮之间的运动学和动力学关系，并结合滚轮、浮动衬套和滚轮销之间的瞬态(每一时刻)润滑参数(油膜厚度相关参数)进行浮动衬套的双层油膜润滑状态实时异常分析，能够实现更为准确的浮动衬套双层油膜润滑异常分析判断。
附图说明
为了更清楚地说明本发明或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明基于浮动衬套的双层油膜润滑分析方法的第一流程图；
图2是本发明基于浮动衬套的双层油膜润滑分析方法的第二流程图；
图3是本发明凸轮与滚轮间的动力学模型示意图；
图4是本发明凸轮与滚轮间的运动学模型示意图；
图5是本发明基于浮动衬套的双层油膜润滑分析系统模块图；
图6是本发明提供的电子设备的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合本发明中的附图，对本发明中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
下面结合图1至图4描述本发明的基于浮动衬套的双层油膜润滑分析方法，如图1所示，一种基于浮动衬套的双层油膜润滑分析方法，包括：
步骤101：在凸轮转动过程中，确定每一时刻下浮动衬套内层油膜厚度和外层油膜厚度，以及确定每一时刻下滚轮位移和浮动衬套位移。
其中，每一时刻中各个时刻间隔固定的时间步长，时间步长为预设值。根据每一时刻下第一组合结构(包括凸轮和滚轮)的结构参数，确定每一时刻下的滚轮位移、浮动衬套位移。利用当前时刻浮动衬套内层油膜厚度和外层油膜厚度求解内外双层油膜平均雷诺方程，当压力收敛时，得到当前时刻下浮动衬套内层油膜压力和外层油膜压力；根据内层油膜压力和外层油膜压力分别确定浮动衬套内层弹性变形和外层弹性变形；将内层弹性变形和外层弹性变形代入油膜厚度方程，当弹性变形收敛，且浮动衬套内外层达到热平衡条件时，得到当前时刻下浮动衬套内层油膜厚度和外层油膜厚度，然后计算当前时刻下滚轮位移与浮动衬套位移，根据当下时刻滚轮位移与浮动衬套位移确定当下时刻滚轮运动轨迹和浮动衬套运动轨迹，若滚轮运动轨迹和浮动衬套运动轨迹二者中有一个没有首尾闭合，则根据下一时刻的滚轮、浮动衬套和滚轮销之间的结构参数，确定下一时刻的浮动衬套内层油膜厚度和外层油膜厚度，以此类推，从而得到每一时刻下浮动衬套内层油膜厚度和外层油膜厚度。
步骤102：根据所述每一时刻下滚轮位移和浮动衬套位移分别确定滚轮的运动轨迹和浮动衬套的运动轨迹。
步骤103：当所述滚轮的运动轨迹和浮动衬套的运动轨迹均首尾闭合时，根据所述每一时刻下浮动衬套内层油膜厚度和外层油膜厚度分别确定浮动衬套内层油膜厚度最小值和外层油膜厚度最小值。
步骤104：当浮动衬套内层油膜厚度最小值和外层油膜厚度最小值中任一个小于预设油膜厚度阈值时，确定双层油膜润滑存在异常。
由于凸轮转动会带动滚轮转动，且滚轮、浮动衬套和滚轮销又组成一个联动机构，所以本发明能够确定在凸轮转动过程中，每一时刻下浮动衬套内层油膜厚度和外层油膜厚度，以及确定每一时刻下滚轮位移和浮动衬套位移；根据所述每一时刻下滚轮位移和浮动衬套位移分别确定滚轮的运动轨迹和浮动衬套的运动轨迹；当所述滚轮的运动轨迹和浮动衬套的运动轨迹均首尾闭合时，根据所述每一时刻下浮动衬套内层油膜厚度和外层油膜厚度分别确定浮动衬套内层油膜厚度最小值和外层油膜厚度最小值；当浮动衬套内层油膜厚度最小值和外层油膜厚度最小值中任一个小于预设油膜厚度阈值时，确定双层油膜润滑存在异常。即本发明利用了凸轮和滚轮之间的运动学和动力学关系，并结合滚轮、浮动衬套和滚轮销之间的瞬态(每一时刻)润滑参数(压力和温度相关参数)进行浮动衬套的双层油膜润滑状态实时异常分析，能够实现更为准确的浮动衬套双层油膜润滑异常分析判断。
在一个具体的实施例中，上述步骤101-104具体计算过程如图2所示，主要包括凸轮与滚轮间的动力学与运动学计算过程和基于滚轮、浮动衬套和滚轮销的双层油膜瞬态润滑分析过程两部分。
一方面，动力学与运动学计算过程包括计算滚轮时变载荷和滚轮时变转速两部分。
对于滚轮时变载荷，需要根据力的合成作用来计算。
在一个具体的实施例中，凸轮与滚轮间的动力学模型如图3所示，滚轮时变载荷的具体计算过程包括：
将机构当作刚性体，对于滚轮所受时变载荷F表达式为：
F＝FC+FN+FP
其中，FC为弹簧恢复力折算到滚轮的量，即弹簧对滚轮的作用力。
FC＝F0+c'·s(α)
式中，F0为弹簧预紧力(N)，c'为弹簧刚度(N/m)，s(α)为滚轮升程(m)，α为升程角。
惯性力FN由下式确定：
其中，ω1为凸轮转速(rad/s)，M为凸轮所驱动的质量(kg)，其计算公式为：
式中，Ma为弹簧质量(kg)，Mb为凸轮机构中除了柱塞弹簧、滚轮以外其他关联件质量(例如从动的机架等)(kg)，Mc为柱塞质量(kg)，Md为滚轮质量(kg)。
FP为柱塞油腔液压力，已知油液压力最小值与最大值，可通过线性插值的方法得到柱塞油腔压力变化情况：
其中，d为柱塞直径(m)，p1为最小油液压力(MPa)，p2为最大油液压力(MPa)，smax为滚轮升程最大值。
对于滚轮时变转速的计算，首先获取第一组合结构(包括凸轮和滚轮)的结构参数、运行工况等参数；然后根据凸轮型线参数计算凸轮压力角和曲率半径，进而计算凸轮表面速度，最后利用凸轮与滚轮之间的两表面速度关系计算得到滚轮表面速度，最后根据滚轮表面速度计算滚轮时变转速。
在一个具体的实施例中，凸轮与滚轮间的运动学模型如图4所示，图4表示出了凸轮机构中凸轮曲率半径、压力角等运动学参数，滚轮时变转速的具体计算过程包括：
对于供油凸轮曲率半径，其计算公式为：
其中，g(α)＝R0+R2+s(α)，α为升程角，R1为凸轮曲率半径(m)，R0为凸轮的基圆半径(m)，R2为滚轮外圈半径(m)。
对于凸轮与滚轮副综合曲率半径，其计算公式为：
第一组合结构(包括凸轮和滚轮)的压力角随着凸轮转角的变化而不同，可表示为：
式中，e为偏心距(m)，s为滚轮升程(m)。
凸轮与滚轮两表面之间的速度关系满足：
其中，u1、u2分别为凸轮表面线速度、滚轮表面线速度，ω1为凸轮时变转速(rad/s)，α0为推程角。因此，滚轮时变转速ω2可表示为：
ω2＝u2/R2
另一方面，基于滚轮、浮动衬套和滚轮销的双层油膜瞬态润滑分析过程包括：
根据第二组合结构(包括滚轮、浮动衬套和滚轮销)的结构参数、初始假定位置以及初始热变形量、弹性变形量计算初始时刻(t＝1)下浮动衬套内外层油膜厚度。
在一个具体的实施例中，浮动衬套内、外层油膜厚度方程分别为：
式中，hi、h0分别为浮动衬套内、外层油膜厚度；Ri、Ro分别为浮动衬套内、外圈半径；Rb外滚轮内圈半径；εi、εo分别为浮动衬套内、外层偏心率；分别为浮动衬套内、外层偏位角；Rj为滚轮销半径；δTi、δTo分别为内、外层热变形导致的间隙改变量，δPi、δPo分别为浮动衬套内、外层弹性变形量。
在得到浮动衬套内、外层油膜厚度的基础上，以得到的滚轮时变转速为基础，利用有限差分法求解内、外双层油膜平均雷诺方程，循环迭代直到满足压力收敛，得到内、外层油膜压力；其中，循环迭代过程中，压力边界条件采用雷诺边界条件。
在一个具体的实施例中，内、外双层油膜平均雷诺方程分别为：
式中，θ为周向坐标；y为轴向坐标pi、po分别为浮动衬套内、外层油膜压力；ωr为浮动衬套转速；ω2为滚轮转速；φx、φy为压力流量因子；φc为接触因子；φs为剪切流量因子；σi·σo分别为内、外层表面综合粗糙度；ρ为润滑油密度；μi、μo分别为内、外层考虑温黏效应润滑油黏度。
利用变形矩阵法，根据内、外层油膜压力对应计算浮动衬套内、外层各个节点的弹性变形，将内、外层弹性变形量分别代入浮动衬套内、外层油膜厚度方程，进行弹性变形收敛判断，若满足收敛条件，则跳出循环，若不满足收敛条件，则更新弹性变形量，重复之前的油膜压力计算。
在一个具体的实施例中，利用变形矩阵法，计算浮动衬套内、外层各个节点的弹性变形表示为：
其中，δpi(θ,y)、δpo(θ,y)分别为浮动衬套内、外层在节点(θ,y)处的弹性变形量；pi(θ,y)、po(θ,y)分别为浮动衬套内、外层在节点(θ,y)处的油膜压力；分别为浮动衬套内、外层弹性变形矩阵，即在表面(θ′,y′)节点作用单位压力，在(θ,y)节点产生的弹性变形。
利用内、外层油膜热平衡方程计算内、外层油膜平衡温度，即判断摩擦产热与端泄热是否满足相等条件，若不满足，则更新内、外层油膜温度，进而根据更新后的内、外层油膜温度计算内、外层热变形量，并将初始热变形量进行更新；若满足，则以滚轮时变载荷为基础，利用差分法求解滚轮、浮动衬套运动学方程和浮动衬套转速方程，得到下一时刻滚轮位移、浮动衬套位移和浮动衬套转速。
各个时刻进行循环迭代时，当前时刻的计算结果满足判断条件后，则将当前时刻加一个固定的时间步长，得到下一时刻；且各个时刻之间间隔的时间步长固定，时间步长为预设值。例如在t＝1时，若满足内外膜热平衡条件，则对t＝1加一个固定的时间步长，得到t＝2；当t＝2满足内外膜热平衡条件，则对t＝2加一个固定的时间步长，得到t＝3，以此类推。
在一个具体的实施例中，内、外层油膜热平衡方程分别为：
cpρQiΔTi＝Wi
cpρQoΔTo＝Wo
式中，cp为滑油比热；ρ为润滑油密度；ΔTi、ΔTo分别为内、外层温升；Qi、Qo分别为内、外层端泄流量，其计算公式为：
式中，y＝Li、y＝Lo分别指浮动衬套内层端侧位置、外层端侧位置。
Wi、Wo分别为内、外层摩擦损失，其计算公式为：
在一个具体的实施例中，内、外层热变形量计算公式分别为：
δTi＝αRRiΔTi-αJRjΔTi
δTo＝αBRbΔTo-αRRoΔTo
式中，δTi为内层热变形量，δTo为外层热变形量，αJ、αR、αB分别为滚轮销的热膨胀系数、浮动衬套的热膨胀系数和滚轮的热膨胀系数。
在一个具体的实施例中，滚轮运动学方程为：
浮动衬套运动学方程为：
浮动衬套转速方程为：
其中，m1是滚轮质量，m2是浮动衬套质量，F为滚轮所受时变载荷，Pxi、Pyi分别为浮动衬套内层水平、竖直方向油膜承载力，Pxo、Pyo分别为浮动衬套外层水平、竖直方向油膜承载力，Γi、Γo分别对应内、外层油膜摩擦力矩，I2为浮动衬套转动惯量，为浮动衬套角加速度，/>分别对应滚轮水平、竖直加速度，/>分别对应浮动衬套水平、竖直加速度。
根据各个时刻的滚轮位移、浮动衬套位移和浮动衬套转速可以得到滚轮的运动轨迹和浮动衬套的运动轨迹，当滚轮的运动轨迹和浮动衬套的运动轨迹均首尾闭合时，即轨迹稳定时，此时跳出循环，输出浮动衬套内层油膜厚度最小值、外层油膜厚度最小值；当浮动衬套内层油膜厚度最小值、外层油膜厚度最小值中任一个小于预设油膜厚度阈值时，确定双层油膜润滑存在异常；当滚轮的运动轨迹和浮动衬套的运动轨迹二者中有一个没有首尾闭合时，则更新滚轮初始位移、浮动衬套初始位移和浮动衬套初始转速，继续迭代。
在一个具体的实施例中，本发明还可以通过判断凸轮转角每转360°时对应得到的滚轮运动轨迹和浮动衬套运动轨迹是否稳定来判断是否需要跳出循环，当二者轨迹均首尾闭合(即稳定)时，跳出循环，输出浮动衬套内层油膜厚度最小值、外层油膜厚度最小值；当浮动衬套内层油膜厚度最小值、外层油膜厚度最小值中任一个小于预设油膜厚度阈值时，确定双层油膜润滑存在异常；当滚轮的运动轨迹和浮动衬套的运动轨迹中有一个没有首尾闭合时，则更新滚轮初始位移、浮动衬套初始位移和浮动衬套初始转速，继续迭代。
在一个具体的实施例中，预设油膜厚度阈值为0.8微米到1.2微米范围内的任意值。
根据本发明提供的实施例，本发明还公开了如下技术效果：
1、本发明考虑了浮动衬套内外层热弹变形、以及浮动衬套内外层油膜混合润滑效应，更符合实际润滑状态，所以基于此进行双层油膜润滑状态判断更为准确。
2、本发明将凸轮和滚轮之间的运动学、动力学关系与滚轮、浮动衬套、滚轮销之间的双层油膜的瞬态混合热弹流润滑相耦合，考虑了滚轮的时变转速与时变载荷，更能准确的反应各个时刻的滚轮、浮动衬套和滚轮销工作状态，基于此能够提高对浮动衬套双层油膜润滑状态分析的准确性。
下面对本发明提供的基于浮动衬套的双层油膜润滑分析系统进行描述，下文描述基于浮动衬套的双层油膜润滑分析系统与上文描述基于浮动衬套的双层油膜润滑分析方法可相互对应参照。
如图5所示，一种基于浮动衬套的双层油膜润滑分析系统，包括：
第一参数确定模块501，用于在凸轮转动过程中，确定每一时刻下浮动衬套内层油膜厚度和外层油膜厚度，以及确定每一时刻下滚轮位移和浮动衬套位移。
运动轨迹确定模块502，用于根据所述每一时刻下滚轮位移和浮动衬套位移分别确定滚轮的运动轨迹和浮动衬套的运动轨迹。
第二参数确定模块503，用于当所述滚轮的运动轨迹和浮动衬套的运动轨迹均首尾闭合时，根据所述每一时刻下浮动衬套内层油膜厚度和外层油膜厚度分别确定浮动衬套内层油膜厚度最小值和外层油膜厚度最小值。
双层油膜润滑分析模块504，用于当浮动衬套内层油膜厚度最小值和外层油膜厚度最小值中任一个小于预设油膜厚度阈值时，确定双层油膜润滑存在异常。
图6示例了一种电子设备的实体结构示意图，如图6所示，该电子设备可以包括：处理器(processor)610、通信接口(Communications Interface)620、存储器(memory)830和通信总线640，其中，处理器610，通信接口620，存储器630通过通信总线640完成相互间的通信。处理器610可以调用存储器630中的逻辑指令，以执行基于浮动衬套的双层油膜润滑分析方法，该方法包括：
在凸轮转动过程中，确定每一时刻下浮动衬套内层油膜厚度和外层油膜厚度，以及确定每一时刻下滚轮位移和浮动衬套位移。
根据所述每一时刻下滚轮位移和浮动衬套位移分别确定滚轮的运动轨迹和浮动衬套的运动轨迹。
当所述滚轮的运动轨迹和浮动衬套的运动轨迹均首尾闭合时，根据所述每一时刻下浮动衬套内层油膜厚度和外层油膜厚度分别确定浮动衬套内层油膜厚度最小值和外层油膜厚度最小值。
当浮动衬套内层油膜厚度最小值和外层油膜厚度最小值中任一个小于预设油膜厚度阈值时，确定双层油膜润滑存在异常。
此外，上述的存储器630中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
另一方面，本发明还提供一种计算机程序产品，所述计算机程序产品包括计算机程序，计算机程序可存储在非暂态计算机可读存储介质上，所述计算机程序被处理器执行时，计算机能够执行以执行基于浮动衬套的双层油膜润滑分析方法，该方法包括：
在凸轮转动过程中，确定每一时刻下浮动衬套内层油膜厚度和外层油膜厚度，以及确定每一时刻下滚轮位移和浮动衬套位移。
根据所述每一时刻下滚轮位移和浮动衬套位移分别确定滚轮的运动轨迹和浮动衬套的运动轨迹。
当所述滚轮的运动轨迹和浮动衬套的运动轨迹均首尾闭合时，根据所述每一时刻下浮动衬套内层油膜厚度和外层油膜厚度分别确定浮动衬套内层油膜厚度最小值和外层油膜厚度最小值。
当浮动衬套内层油膜厚度最小值和外层油膜厚度最小值中任一个小于预设油膜厚度阈值时，确定双层油膜润滑存在异常。
又一方面，本发明还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现以执行基于浮动衬套的双层油膜润滑分析方法，该方法包括：
在凸轮转动过程中，确定每一时刻下浮动衬套内层油膜厚度和外层油膜厚度，以及确定每一时刻下滚轮位移和浮动衬套位移。
根据所述每一时刻下滚轮位移和浮动衬套位移分别确定滚轮的运动轨迹和浮动衬套的运动轨迹。
当所述滚轮的运动轨迹和浮动衬套的运动轨迹均首尾闭合时，根据所述每一时刻下浮动衬套内层油膜厚度和外层油膜厚度分别确定浮动衬套内层油膜厚度最小值和外层油膜厚度最小值。
当浮动衬套内层油膜厚度最小值和外层油膜厚度最小值中任一个小于预设油膜厚度阈值时，确定双层油膜润滑存在异常。
以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下，即可以理解并实施。
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件。基于这样的理解，上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (9)
1.一种基于浮动衬套的双层油膜润滑分析方法，其特征在于，包括：
在凸轮转动过程中，确定每一时刻下浮动衬套内层油膜厚度和外层油膜厚度，以及确定每一时刻下滚轮位移和浮动衬套位移；
根据所述每一时刻下滚轮位移和浮动衬套位移分别确定滚轮的运动轨迹和浮动衬套的运动轨迹；
当所述滚轮的运动轨迹和浮动衬套的运动轨迹均首尾闭合时，根据所述每一时刻下浮动衬套内层油膜厚度和外层油膜厚度分别确定浮动衬套内层油膜厚度最小值和外层油膜厚度最小值；
当浮动衬套内层油膜厚度最小值和外层油膜厚度最小值中任一个小于预设油膜厚度阈值时，确定双层油膜润滑存在异常；
其中，确定每一时刻下浮动衬套内层油膜厚度和外层油膜厚度，包括：
根据第二组合结构的结构参数计算浮动衬套当前时刻的内层油膜厚度和外层油膜厚度，所述第二组合结构包括滚轮、浮动衬套和滚轮销；
根据所述滚轮位移与浮动衬套位移确定当前时刻下滚轮运动轨迹和浮动衬套运动轨迹，若滚轮运动轨迹和浮动衬套运动轨迹二者中任一个没有首尾闭合，则获取下一时刻的第二组合结构的结构参数；
根据所述下一时刻的第二组合结构的结构参数确定下一时刻的浮动衬套内层油膜厚度和外层油膜厚度。
2.根据权利要求1所述的基于浮动衬套的双层油膜润滑分析方法，其特征在于，确定每一时刻下浮动衬套内层油膜厚度和外层油膜厚度，具体包括：
获取第一组合结构的结构参数，所述第一组合结构包括凸轮和滚轮；
根据所述第一组合结构的结构参数计算滚轮时变转速；
获取第二组合结构的结构参数，所述第二组合结构包括滚轮、浮动衬套和滚轮销；
根据所述第二组合结构的结构参数计算浮动衬套当前时刻的内层油膜厚度和外层油膜厚度；
根据所述滚轮时变转速、内层油膜厚度和外层油膜厚度求解内外双层油膜平均雷诺方程，当压力收敛时，得到当前时刻下浮动衬套内层油膜压力和外层油膜压力；
根据所述内层油膜压力和外层油膜压力分别确定浮动衬套内层弹性变形和外层弹性变形；
将所述内层弹性变形和外层弹性变形代入油膜厚度方程，当弹性变形收敛、浮动衬套内外层达到热平衡条件时，输出所述当前时刻下浮动衬套内层油膜压力和外层油膜压力；
确定当前时刻的滚轮位移与浮动衬套位移；
根据所述滚轮位移与浮动衬套位移确定当前时刻下滚轮运动轨迹和浮动衬套运动轨迹；
若滚轮运动轨迹和浮动衬套运动轨迹二者中任一个没有首尾闭合，则获取下一时刻的第二组合结构的结构参数；
根据所述下一时刻的第二组合结构的结构参数确定下一时刻的浮动衬套内层油膜厚度和外层油膜厚度。
3.根据权利要求2所述的基于浮动衬套的双层油膜润滑分析方法，其特征在于，根据所述第一组合结构的结构参数计算滚轮时变转速，具体包括：
根据凸轮型线参数计算凸轮压力角和凸轮曲率半径；
根据所述凸轮压力角和凸轮曲率半径计算凸轮表面速度；
根据所述凸轮表面速度计算滚轮表面速度；
根据所述滚轮表面速度计算滚轮时变转速。
4.根据权利要求1所述的基于浮动衬套的双层油膜润滑分析方法，其特征在于，确定每一时刻下滚轮位移和浮动衬套位移，具体包括：
获取第一组合结构的结构参数，所述第一组合结构包括凸轮和滚轮；
根据所述第一组合结构的结构参数计算滚轮时变载荷；
根据所述滚轮时变载荷、滚轮与浮动衬套之间的运动学方程确定每一时刻下滚轮位移和浮动衬套位移。
5.根据权利要求4所述的基于浮动衬套的双层油膜润滑分析方法，其特征在于，所述滚轮时变载荷包括弹簧对滚轮的作用力、惯性力和柱塞油腔液压力。
6.根据权利要求1所述的基于浮动衬套的双层油膜润滑分析方法，其特征在于，所述预设油膜厚度阈值为0.8微米到1.2微米范围内的任意值。
7.一种基于浮动衬套的双层油膜润滑分析系统，其特征在于，包括：
第一参数确定模块，用于在凸轮转动过程中，确定每一时刻下浮动衬套内层油膜厚度和外层油膜厚度，以及确定每一时刻下滚轮位移和浮动衬套位移；
运动轨迹确定模块，用于根据所述每一时刻下滚轮位移和浮动衬套位移分别确定滚轮的运动轨迹和浮动衬套的运动轨迹；
第二参数确定模块，用于当所述滚轮的运动轨迹和浮动衬套的运动轨迹均首尾闭合时，根据所述每一时刻下浮动衬套内层油膜厚度和外层油膜厚度分别确定浮动衬套内层油膜厚度最小值和外层油膜厚度最小值；
双层油膜润滑分析模块，用于当浮动衬套内层油膜厚度最小值和外层油膜厚度最小值中任一个小于预设油膜厚度阈值时，确定双层油膜润滑存在异常；
其中，确定每一时刻下浮动衬套内层油膜厚度和外层油膜厚度，包括：
根据第二组合结构的结构参数计算浮动衬套当前时刻的内层油膜厚度和外层油膜厚度，所述第二组合结构包括滚轮、浮动衬套和滚轮销；
根据所述滚轮位移与浮动衬套位移确定当前时刻下滚轮运动轨迹和浮动衬套运动轨迹，若滚轮运动轨迹和浮动衬套运动轨迹二者中任一个没有首尾闭合，则获取下一时刻的第二组合结构的结构参数；
根据所述下一时刻的第二组合结构的结构参数确定下一时刻的浮动衬套内层油膜厚度和外层油膜厚度。
8.一种电子设备，包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现如权利要求1至6任一项所述的基于浮动衬套的双层油膜润滑分析方法。
9.一种非暂态计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至6任一项所述的基于浮动衬套的双层油膜润滑分析方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310404676.6A CN116415462B (zh) | 2023-04-14 | 2023-04-14 | 基于浮动衬套的双层油膜润滑分析方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310404676.6A CN116415462B (zh) | 2023-04-14 | 2023-04-14 | 基于浮动衬套的双层油膜润滑分析方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116415462A CN116415462A (zh) | 2023-07-11 |
| CN116415462B true CN116415462B (zh) | 2023-11-17 |
Family
ID=87054417
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310404676.6A Active CN116415462B (zh) | 2023-04-14 | 2023-04-14 | 基于浮动衬套的双层油膜润滑分析方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116415462B (zh) |
Citations (25)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS62282715A (ja) * | 1986-05-30 | 1987-12-08 | Kawasaki Steel Corp | 圧延機用ロ−ルの油膜軸受系統の異常診断方法 |
| JPH03160325A (ja) * | 1989-11-17 | 1991-07-10 | Nippon Seiko Kk | 転がり軸受の潤滑状態監視方法 |
| JPH11294364A (ja) * | 1998-04-09 | 1999-10-26 | Hitachi Ltd | 軸受け油膜厚さ制御方法 |
| KR20050008890A (ko) * | 2003-07-14 | 2005-01-24 | 두산중공업 주식회사 | 패드형 가이드 베어링의 유막 특성 계산 방법 |
| JP2009150434A (ja) * | 2007-12-19 | 2009-07-09 | Ntn Corp | 転がり軸受の保持器およびその設計方法 |
| JP2011058535A (ja) * | 2009-09-08 | 2011-03-24 | Nippon Steel Corp | 圧延ロール軸受部に供給する潤滑油の冷却方法 |
| CN102589489A (zh) * | 2012-02-09 | 2012-07-18 | 西安交通大学 | 一种圆柱滚子轴承润滑油膜厚分布的检测方法 |
| CN103615994A (zh) * | 2013-11-06 | 2014-03-05 | 西安交通大学 | 滚动轴承润滑油膜厚度测试装置及其测试方法 |
| EP2734755A1 (en) * | 2011-07-20 | 2014-05-28 | Aktiebolaget SKF | Method of determining when a critical film thickness will be reached in a grease-lubricated seal |
| CN104517029A (zh) * | 2014-10-28 | 2015-04-15 | 北京理工大学 | 一种润滑油粘度选取方法 |
| CN104776998A (zh) * | 2015-03-26 | 2015-07-15 | 北京工业大学 | 一种基于动态刚度系数和阻尼系数的转子轴心轨迹求解方法 |
| CN106842202A (zh) * | 2017-02-28 | 2017-06-13 | 黄晓霞 | 一种基于遥感影像确定溢油事故源点的方法及装置 |
| CN109002569A (zh) * | 2018-01-12 | 2018-12-14 | 哈尔滨理工大学 | 一种建立静压油膜运动层边界条件的方法 |
| CN109186473A (zh) * | 2018-09-06 | 2019-01-11 | 常州大学 | 滑滚工况下滚子摩擦副的油膜厚度以及油膜拖动系数测试装置 |
| CN110206608A (zh) * | 2019-07-08 | 2019-09-06 | 杭州新坐标科技股份有限公司 | 一种带储油功能的滚轮摇臂 |
| JP2019203561A (ja) * | 2018-05-24 | 2019-11-28 | 株式会社ジェイテクト | 油膜厚さ算出装置 |
| KR20200045153A (ko) * | 2018-10-22 | 2020-05-04 | 현대자동차주식회사 | 엔진의 유로차압방식 캠 샤프트 윤활시스템 |
| CN112710212A (zh) * | 2020-12-21 | 2021-04-27 | 重庆江增船舶重工有限公司 | 一种提高浮动套内孔油楔检测效率的加工和检测方法 |
| CN113503197A (zh) * | 2021-07-12 | 2021-10-15 | 哈尔滨工程大学 | 一种考虑结构振动的船用凸轮-挺柱副弹流润滑分析方法 |
| CN113686231A (zh) * | 2021-08-23 | 2021-11-23 | 华能威海发电有限责任公司 | 一种汽轮机发电机油膜厚度的测量方法 |
| WO2022071164A1 (ja) * | 2020-09-29 | 2022-04-07 | 日本精工株式会社 | 油膜の状態検出方法、状態検出装置、およびプログラム |
| CN114543716A (zh) * | 2022-03-16 | 2022-05-27 | 哈尔滨工业大学 | 基于超声波反射系数实部的油膜厚度测量方法 |
| CN114611433A (zh) * | 2022-03-22 | 2022-06-10 | 郑州大学 | 一种耦合流态与粗糙度的动静压浮环轴承模型计算方法 |
| JP2022094088A (ja) * | 2020-12-14 | 2022-06-24 | Ntn株式会社 | 油膜パラメータの算出方法 |
| CN115013296A (zh) * | 2022-05-24 | 2022-09-06 | 北京航空航天大学 | 一种轴向柱塞泵滑靴副油膜厚度确定方法及系统 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| BE394136A (zh) * | ||||
| US20170102292A1 (en) * | 2015-10-08 | 2017-04-13 | United Technologies Corporation | Systems and methods for monitoring lubricant film thickness of a journal bearing in an epicyclic gear system of an operating gas turbine engine |
| EP3187768B1 (en) * | 2015-12-17 | 2023-03-15 | Trane International Inc. | System and method for dynamically determining refrigerant film thickness and dynamically controlling refrigerant film thickness at rolling-element bearing of an oil free chiller |
| US10760977B2 (en) * | 2017-05-17 | 2020-09-01 | General Electric Company | System and method for detecting lubricated bearing condition |
| DE102017223418B4 (de) * | 2017-12-20 | 2023-05-25 | Zf Friedrichshafen Ag | Modellbasiertes Verfahren und System zur Zustandsüberwachung eines Gleitlagers, insbesondere für Windkraftanlagen |
-
2023
- 2023-04-14 CN CN202310404676.6A patent/CN116415462B/zh active Active
Patent Citations (25)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS62282715A (ja) * | 1986-05-30 | 1987-12-08 | Kawasaki Steel Corp | 圧延機用ロ−ルの油膜軸受系統の異常診断方法 |
| JPH03160325A (ja) * | 1989-11-17 | 1991-07-10 | Nippon Seiko Kk | 転がり軸受の潤滑状態監視方法 |
| JPH11294364A (ja) * | 1998-04-09 | 1999-10-26 | Hitachi Ltd | 軸受け油膜厚さ制御方法 |
| KR20050008890A (ko) * | 2003-07-14 | 2005-01-24 | 두산중공업 주식회사 | 패드형 가이드 베어링의 유막 특성 계산 방법 |
| JP2009150434A (ja) * | 2007-12-19 | 2009-07-09 | Ntn Corp | 転がり軸受の保持器およびその設計方法 |
| JP2011058535A (ja) * | 2009-09-08 | 2011-03-24 | Nippon Steel Corp | 圧延ロール軸受部に供給する潤滑油の冷却方法 |
| EP2734755A1 (en) * | 2011-07-20 | 2014-05-28 | Aktiebolaget SKF | Method of determining when a critical film thickness will be reached in a grease-lubricated seal |
| CN102589489A (zh) * | 2012-02-09 | 2012-07-18 | 西安交通大学 | 一种圆柱滚子轴承润滑油膜厚分布的检测方法 |
| CN103615994A (zh) * | 2013-11-06 | 2014-03-05 | 西安交通大学 | 滚动轴承润滑油膜厚度测试装置及其测试方法 |
| CN104517029A (zh) * | 2014-10-28 | 2015-04-15 | 北京理工大学 | 一种润滑油粘度选取方法 |
| CN104776998A (zh) * | 2015-03-26 | 2015-07-15 | 北京工业大学 | 一种基于动态刚度系数和阻尼系数的转子轴心轨迹求解方法 |
| CN106842202A (zh) * | 2017-02-28 | 2017-06-13 | 黄晓霞 | 一种基于遥感影像确定溢油事故源点的方法及装置 |
| CN109002569A (zh) * | 2018-01-12 | 2018-12-14 | 哈尔滨理工大学 | 一种建立静压油膜运动层边界条件的方法 |
| JP2019203561A (ja) * | 2018-05-24 | 2019-11-28 | 株式会社ジェイテクト | 油膜厚さ算出装置 |
| CN109186473A (zh) * | 2018-09-06 | 2019-01-11 | 常州大学 | 滑滚工况下滚子摩擦副的油膜厚度以及油膜拖动系数测试装置 |
| KR20200045153A (ko) * | 2018-10-22 | 2020-05-04 | 현대자동차주식회사 | 엔진의 유로차압방식 캠 샤프트 윤활시스템 |
| CN110206608A (zh) * | 2019-07-08 | 2019-09-06 | 杭州新坐标科技股份有限公司 | 一种带储油功能的滚轮摇臂 |
| WO2022071164A1 (ja) * | 2020-09-29 | 2022-04-07 | 日本精工株式会社 | 油膜の状態検出方法、状態検出装置、およびプログラム |
| JP2022094088A (ja) * | 2020-12-14 | 2022-06-24 | Ntn株式会社 | 油膜パラメータの算出方法 |
| CN112710212A (zh) * | 2020-12-21 | 2021-04-27 | 重庆江增船舶重工有限公司 | 一种提高浮动套内孔油楔检测效率的加工和检测方法 |
| CN113503197A (zh) * | 2021-07-12 | 2021-10-15 | 哈尔滨工程大学 | 一种考虑结构振动的船用凸轮-挺柱副弹流润滑分析方法 |
| CN113686231A (zh) * | 2021-08-23 | 2021-11-23 | 华能威海发电有限责任公司 | 一种汽轮机发电机油膜厚度的测量方法 |
| CN114543716A (zh) * | 2022-03-16 | 2022-05-27 | 哈尔滨工业大学 | 基于超声波反射系数实部的油膜厚度测量方法 |
| CN114611433A (zh) * | 2022-03-22 | 2022-06-10 | 郑州大学 | 一种耦合流态与粗糙度的动静压浮环轴承模型计算方法 |
| CN115013296A (zh) * | 2022-05-24 | 2022-09-06 | 北京航空航天大学 | 一种轴向柱塞泵滑靴副油膜厚度确定方法及系统 |
Non-Patent Citations (2)
| Title |
|---|
| 华德良.船用柴油机供油凸轮-滚轮副弹流润滑分析.全文. * |
| 郭怀谦.船用柴油机供油机构滚轮衬套混合热弹流润滑特性研究.全文. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116415462A (zh) | 2023-07-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109753723B (zh) | 一种向心滚动轴承疲劳寿命计算方法 | |
| Lu et al. | The analysis of secondary motion and lubrication performance of piston considering the piston skirt profile | |
| CN116415462B (zh) | 基于浮动衬套的双层油膜润滑分析方法及系统 | |
| Kushwaha et al. | Valve-train dynamics: a simplified tribo-elasto-multi-body analysis | |
| Liu et al. | Effects of bushing profiles on the elastohydrodynamic lubrication performance of the journal bearing under steady operating conditions | |
| Xu et al. | Modeling and simulation of aero-hydraulic pump wear failure | |
| CN115470584B (zh) | 应用于行星齿轮与滚动轴承耦合系统的动力学建模方法 | |
| Nuruzzaman et al. | Study on pressure distribution and load capacity of a journal bearing using finite element method and analytical method | |
| CN112069706B (zh) | 一种坦克炮塔座圈摩擦学和动力学耦合分析方法 | |
| CN111177898A (zh) | 基于bp神经网络的滚动轴承-转子系统耦合性能求解方法 | |
| McClure | Numerical modeling of piston secondary motion and skirt lubrication in internal combustion engines | |
| CN112836445B (zh) | 一种用于求解滚动轴承的动态接触混合润滑参数的方法 | |
| CN112861283B (zh) | 一种计算曲轴与轴承耦合特性的方法及系统 | |
| CN115688297A (zh) | 一种应用于往复机械的动载滑动轴承设计方法 | |
| Raju et al. | Bearing life evaluation of wheel hub ball bearing based on finite element analysis | |
| Mancò et al. | Miniature gerotor pump prototype for automotive applications | |
| Savchuk et al. | Providing of Sliding Bearings Reliability of Transmissions Gear Wheels of Transport Cars by Optimization of Assembly Tolerances | |
| Meng et al. | A numerical model for piston pin lubrication in internal combustion engines | |
| Meng | Modeling of piston pin lubrication in internal combustion engines | |
| Henry | Bell helicopter advanced rotorcraft transmission (ART) program | |
| Wang et al. | Nonlinear tribo-dynamic performance and transient stability for marine dynamically loaded offset-halves journal bearings | |
| CN116186933A (zh) | 一种考虑活塞与缸套热变形的裙部润滑特性仿真方法 | |
| CN117057177A (zh) | 一种轴颈动态偏转的径向滑动轴承摩擦动力学特性计算方法 | |
| CN114021281B (zh) | 基于影响系数法和润滑影响的轮毂轴承复合变形计算方法 | |
| Jang et al. | Dynamic EHL film thickness in cam and follower contacts of various valve lifts |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |