CN109002569A - 一种建立静压油膜运动层边界条件的方法 - Google Patents
一种建立静压油膜运动层边界条件的方法 Download PDFInfo
- Publication number
- CN109002569A CN109002569A CN201810028374.2A CN201810028374A CN109002569A CN 109002569 A CN109002569 A CN 109002569A CN 201810028374 A CN201810028374 A CN 201810028374A CN 109002569 A CN109002569 A CN 109002569A
- Authority
- CN
- China
- Prior art keywords
- oil film
- static pressure
- pressure oil
- film thickness
- boundary condition
- 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
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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Bearings And Hydrostatic Bearings (AREA)
Abstract
静压轴承间隙油膜的膜厚在轴承运行过程中会随着工作转速的不同而发生变化,传统上采取定油膜厚度的研究方法不能够满足实际工况对数值计算结果正确性的要求。本发明通过分析静压轴承工作台运行过程中静压油膜厚度的变化情况,结合润滑介质的粘温方程,建立出静压油膜运动层边界条件,目的是为了缩减理论计算周期,提高静压油膜润滑特性计算的精度。步骤A、定义初始静止时油膜厚度值及运行中油膜厚度变化规律。步骤B、推导静压润滑油粘温表达式。步骤C、定义静压油膜运动区域边界条件。步骤D、编写静压油膜模型运动层程序及变粘度程序。本发明应用于一种建立静压油膜运动层边界条件的方法。
Description
技术领域
本发明涉及了一种建立静压油膜运动层边界条件的方法。
背景技术
静压轴承间隙油膜的膜厚在轴承运行过程中会随着工作转速的不同而发生变化,传统上采取定油膜厚度的研究方法不能够满足实际工况对数值计算结果正确性的要求,若由人工一次次地建立不同油膜厚度的模型并加以仿真模拟,仿真工作量将会变得非常巨大,也会影响课题研究进度和研究周期,对研究人员的技术要求也比较高。基于此,本研究通过分析静压轴承工作台运行过程中静压油膜厚度的变化情况,结合润滑介质的粘温方程,建立出静压油膜运动层边界条件,目的是为了缩减理论计算周期,提高静压油膜润滑特性计算的精度。
发明内容
一种建立静压油膜运动层边界条件的方法
一种建立静压油膜运动层边界条件的方法按以下步骤实现:
步骤A、定义初始静止时油膜厚度值及运行中油膜厚度变化规律:
依据静压油膜载荷条件及结构尺寸确定初始设计油膜厚度值,设定与工作台相接触的边界层油膜为运动边界并以一定的速度向下运动,根据油膜厚度间隙和工况要求给出运动时间,进而得到油膜厚度变化规律。
步骤B、推导静压润滑油粘温表达式:
考虑到静压润滑油膜随着温度变化而变薄,温度对润滑油粘度有影响,在设定流体粘度属性时,采用变粘度即推导出粘度随温度变化的表达式如下:
式中为动力粘度,单位:10-2Pa·s;为润滑油平均温度,单位:K。
步骤C、定义静压油膜运动区域边界条件:
在进行静压油膜动态计算时,需要定义运动区域边界条件。针对静压油膜区域,定义与旋转工作台接触的边界层网格向下运动即形状不变,只是质心向下移动,与边界层相连接的四周边界网格会发生挤压变形,将四周边界定义为变形区域。
步骤D、编写静压油膜模型运动层程序及变粘度程序。
发明效果:
本发明主要是提出一种建立静压油膜运动层边界条件的方法,通过分析静压轴承工作台运行过程中静压油膜的动态变化情况,结合润滑介质的粘温方程,建立出静压油膜运动层边界条件,编写出静压油膜运动层程序及变粘度程序用于后期的求解计算。本发明的研究方法代替了静压油膜动态计算前处理的整个过程,应用此研究方法进行静压油膜求解计算,计算效率可提高近30%,且此方法可以推广到不同结构尺寸静压油膜研究中,大大地缩减了理论计算周期,并提高了静压油膜性能计算精度。
附图说明
图1是静压油膜三维模型图。
具体实施方式
重型静压轴承热油携带判定方法按以下步骤实现:
步骤A、定义初始静止时油膜厚度值及运行中油膜厚度变化规律:
依据静压油膜载荷条件及结构尺寸确定初始设计油膜厚度值,设定与工作台相接触的边界层油膜运动并以一定的速度向下运动,根据油膜厚度间隙和工况要求给出运动时间,进而得到油膜厚度变化规律。
步骤B、推导静压润滑油粘温表达式:
考虑到静压润滑油膜随着温度变化而变薄,温度对润滑油粘度有影响,在设定流体粘度属性时,采用变粘度即推导出粘度随温度变化的表达式如下:
式中为动力黏度,单位:10-2Pa·s,为润滑油平均温度,单位:K。
步骤C、定义静压油膜运动区域边界条件:
在进行静压油膜动态计算时,需要定义运动区域边界条件。针对静压油膜区域,定义与旋转工作台接触的边界层网格向下运动即形状不变,只是质心向下移动,与边界层相连接的四周边界网格会发生挤压变形,将四周边界定义为变形区域。
步骤D、编写静压油膜运动层程序及变粘度程序用于后期的求解计算:
采用VS2012编译的用于控制边界层网格运动的程序如下:
DEFINE_CG_MOTION
{
if(time<=16.0)
vel[2]=-0.000004;
else
vel[2]=0;
}
采用VS2012编译的用于反映静压油膜变粘度边界条件程序如下:
DEFINE_PROPERTY
{
real mu_lam;//液压油粘度值
real a_x=1.877563e-12;//粘温系数
real b_x=11849.26282;
real c_x=192.8886355;
real temp=C_T(c,t); //提取流场温度
if (temp>330)
mu_lam=0.01;
else if (temp>293)
mu_lam=a_x*exp(b_x/(temp+c_x));//计算流体粘度值
else
mu_lam=0.075;
return mu_lam;
}
本实施方式效果:
本发明主要是提出一种建立静压油膜运动层边界条件的方法,通过分析静压轴承工作台运行过程中静压油膜的动态变化情况,结合润滑介质的粘温方程,建立出静压油膜运动层边界条件,编写出静压油膜运动层程序及变粘度程序用于后期的求解计算。本发明的研究方法代替了静压油膜动态计算前处理建模及边界设置的整个过程,应用此研究方法进行静压油膜求解计算,计算效率可提高近30%,且此方法可以推广到不同结构尺寸静压油膜研究中,大大地缩减了理论计算周期,并提高了静压油膜性能计算精度。
Claims (2)
1.一种建立静压油膜运动层边界条件的方法,其特征在于通过分析静压轴承工作台运行过程中静压油膜厚度的变化情况,结合润滑介质的粘温方程,建立静压油膜运动层边界条件,目的是为了缩减求解计算周期,提高静压油膜润滑特性计算精度;
一种建立静压油膜运动层边界条件的方法按以下步骤实现:
步骤A、定义初始静止时油膜厚度值及运行中油膜厚度变化规律:
依据静压油膜载荷条件及结构尺寸确定初始设计油膜厚度值,设定与工作台相接触的边界层油膜为运动边界并以一定的速度向下运动,根据油膜厚度间隙和工况要求给出运动时间,进而得到油膜厚度变化规律;
步骤B、推导静压润滑油粘温表达式:
考虑到静压润滑油膜随着温度变化而变薄,温度对润滑油粘度有影响,在设定流体粘度属性时,采用变粘度即推导出粘度随温度变化的表达式如下:
式中为动力粘度,单位为10-2Pa·s;为润滑油平均温度,单位为K;
步骤C、定义静压油膜运动区域边界条件:
在进行静压油膜动态计算时,需要定义运动区域边界条件,针对静压油膜区域,定义与旋转工作台接触的边界层网格向下运动即形状不变,只是质心向下移动,与边界层相连接的四周边界网格会发生挤压变形,将四周边界定义为变形区域;
步骤D、编写静压油膜模型运动层程序及变粘度程序。
2.根据权利要求1所述的一种建立静压油膜运动层边界条件的方法,其特征在于将编写出的程序直接用于静压油膜性能求解计算,代替了静压油膜动态计算前处理建模及动态边界设置等过程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810028374.2A CN109002569A (zh) | 2018-01-12 | 2018-01-12 | 一种建立静压油膜运动层边界条件的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810028374.2A CN109002569A (zh) | 2018-01-12 | 2018-01-12 | 一种建立静压油膜运动层边界条件的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109002569A true CN109002569A (zh) | 2018-12-14 |
Family
ID=64574141
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810028374.2A Pending CN109002569A (zh) | 2018-01-12 | 2018-01-12 | 一种建立静压油膜运动层边界条件的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109002569A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113158369A (zh) * | 2021-04-19 | 2021-07-23 | 哈尔滨理工大学 | 静压推力轴承油垫封油边的油膜流量模拟监测方法 |
CN116415462A (zh) * | 2023-04-14 | 2023-07-11 | 哈尔滨工程大学 | 基于浮动衬套的双层油膜润滑分析方法及*** |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1005516A1 (en) * | 1997-12-02 | 2000-06-07 | Infineum USA L.P. | Use of molybdenum complexes in lubricating oil compositions for diesel engines |
CN102141084A (zh) * | 2011-03-30 | 2011-08-03 | 哈尔滨理工大学 | 基于膜厚可变的静压推力轴承间隙油膜温度与厚度关系数值模拟方法 |
CN103246786A (zh) * | 2013-05-24 | 2013-08-14 | 哈尔滨理工大学 | 一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法 |
CN104615801A (zh) * | 2014-12-03 | 2015-05-13 | 哈尔滨理工大学 | 重型静压轴承临界润滑状态转速值确定方法 |
CN105257701A (zh) * | 2015-11-16 | 2016-01-20 | 哈尔滨理工大学 | 多油垫式重型静压轴承热油携带因子及热油携带方程 |
EP3042945A1 (en) * | 2013-09-02 | 2016-07-13 | Nippon Steel & Sumitomo Metal Corporation | Composition for forming lubricating coating film, and threaded joint for steel pipe |
CN107506562A (zh) * | 2017-09-29 | 2017-12-22 | 西安科技大学 | 一种水润滑橡胶轴承双向热流固耦合计算方法 |
-
2018
- 2018-01-12 CN CN201810028374.2A patent/CN109002569A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1005516A1 (en) * | 1997-12-02 | 2000-06-07 | Infineum USA L.P. | Use of molybdenum complexes in lubricating oil compositions for diesel engines |
CN102141084A (zh) * | 2011-03-30 | 2011-08-03 | 哈尔滨理工大学 | 基于膜厚可变的静压推力轴承间隙油膜温度与厚度关系数值模拟方法 |
CN103246786A (zh) * | 2013-05-24 | 2013-08-14 | 哈尔滨理工大学 | 一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法 |
EP3042945A1 (en) * | 2013-09-02 | 2016-07-13 | Nippon Steel & Sumitomo Metal Corporation | Composition for forming lubricating coating film, and threaded joint for steel pipe |
CN104615801A (zh) * | 2014-12-03 | 2015-05-13 | 哈尔滨理工大学 | 重型静压轴承临界润滑状态转速值确定方法 |
CN105257701A (zh) * | 2015-11-16 | 2016-01-20 | 哈尔滨理工大学 | 多油垫式重型静压轴承热油携带因子及热油携带方程 |
CN107506562A (zh) * | 2017-09-29 | 2017-12-22 | 西安科技大学 | 一种水润滑橡胶轴承双向热流固耦合计算方法 |
Non-Patent Citations (3)
Title |
---|
ZHANG YAN-QIN 等: "Multi-pad hydrostatic bearing oil film temperature characteristics under different working mediums", 《2016 11TH INTERNATIONAL CONFERENCE ON COMPUTER SCIENCE & EDUCATION (ICCSE)》 * |
于晓东 等: "高速重载静压推力轴承温度场速度特性", 《哈尔滨理工大学学报》 * |
张艳芹 等: "基于动网格方法的不同油膜厚度下静压轴承承载特性分析", 《哈尔滨理工大学学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113158369A (zh) * | 2021-04-19 | 2021-07-23 | 哈尔滨理工大学 | 静压推力轴承油垫封油边的油膜流量模拟监测方法 |
CN113158369B (zh) * | 2021-04-19 | 2023-11-28 | 哈尔滨理工大学 | 静压推力轴承油垫封油边的油膜流量模拟监测方法 |
CN116415462A (zh) * | 2023-04-14 | 2023-07-11 | 哈尔滨工程大学 | 基于浮动衬套的双层油膜润滑分析方法及*** |
CN116415462B (zh) * | 2023-04-14 | 2023-11-17 | 哈尔滨工程大学 | 基于浮动衬套的双层油膜润滑分析方法及*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | Effect of surface roughness parameters on mixed lubrication characteristics | |
Zhang et al. | Analysis of the flow dynamics characteristics of an axial piston pump based on the computational fluid dynamics method | |
Castilla et al. | Three-dimensional numerical simulation of an external gear pump with decompression slot and meshing contact point | |
CN103226633B (zh) | 一种静压油垫流场特性仿真优化方法 | |
CN109002569A (zh) | 一种建立静压油膜运动层边界条件的方法 | |
Abbassi et al. | Exact solution of three-dimensional unsteady stagnation flow on a heated plate | |
CN104091022A (zh) | 一种速度滑移工况下静压主轴***流固耦合仿真分析方法 | |
Markopoulos et al. | Cutting mechanics and analytical modeling | |
De Boer et al. | Two-scale EHL: Three-dimensional topography in tilted-pad bearings | |
Manring et al. | The control torque on the swash plate of an axial-piston pump utilizing piston-bore springs | |
Balci | The effect of punch speed on frictional contact mechanics of finite-thickness graded layer resting on the rigid foundation | |
Kim et al. | Two Dimensional Cfd Analysis Of A Hydraulic Gear Pump | |
Quan et al. | Numerical analysis of the water film characteristics in the eccentric state of a radial piston pump | |
Wu et al. | A mixed EHL analysis method for grease and formulas for film thickness and asperity load | |
CN113849933B (zh) | 一种大功率船用柴油机润滑***一维仿真方法 | |
CN104517029B (zh) | 一种润滑油粘度选取方法 | |
CN106503391B (zh) | 一种考虑流固耦合的矩形静压油垫快捷计算方法 | |
Han et al. | Improved Implicit SPH Method for simulating free surface flows of power law fluids | |
Wang et al. | A study of the effect of model geometry and lubricant rheology upon the elastohydrodynamic lubrication performance of metal-on-metal hip joints | |
CN101738166B (zh) | 一种高压润滑油界面滑移长度测量方法 | |
Grün et al. | Computational fluid dynamics of the lubricant flow in the sealing gap of rotary shaft seals | |
Xu et al. | Boundary Lubrication in Transient Elliptical Contact: Part 1—Theoretical Formulation and Results | |
Jiang et al. | An integrated model of hydrodynamic lubricating for piston/cylinder interface | |
Chung et al. | The mechanics of clearance in a non-Newtonian lubrication layer | |
CN110119576A (zh) | 一种基于adina二次开发的点接触弹流润滑仿真分析方法 |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20181214 |
|
WD01 | Invention patent application deemed withdrawn after publication |