CN104598666A - 大尺度静压支承临界载荷参数获取方法 - Google Patents
大尺度静压支承临界载荷参数获取方法 Download PDFInfo
- Publication number
- CN104598666A CN104598666A CN201410710729.8A CN201410710729A CN104598666A CN 104598666 A CN104598666 A CN 104598666A CN 201410710729 A CN201410710729 A CN 201410710729A CN 104598666 A CN104598666 A CN 104598666A
- Authority
- CN
- China
- Prior art keywords
- hydrostatic support
- critical
- oil film
- large scale
- load
- 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.)
- Granted
Links
Landscapes
- Magnetic Bearings And Hydrostatic Bearings (AREA)
Abstract
工作台逆时针旋转时,油膜受到剪切流与压差流的影响,从而产生油膜润滑为零,导致润滑失效的情况。本发明通过理论推导,计算出大尺度静压支承油膜承载载荷的理论值。并利用Fluent软件对理论计算的临界载荷值进行验证。步骤A、理论推导大尺度静压支承油膜承载能力表达式。步骤B、计算出不同转速对应的油腔压力及理论临界外载荷值。步骤C、利用前处理软件,对静压支承油膜建立三维仿真模型步骤。D、应用流体计算软件Fluent,对静压支承油膜进行仿真分析,得到大尺度静压支承的仿真结果。步骤E、通过此数值方法验证大尺度静压支承油膜载荷临界状态数学模型的正确性。本发明应用于大尺度静压支承临界载荷参数获取方法研究。
Description
技术领域
本发明涉及了大尺度静压支承在不同转速下油膜临界载荷参数获取方法研究。
背景技术
大尺度多油垫静压支承作为重型数控装备的核心部件,其润滑参数是决定支承乃至整个设备可靠运行的重要指标。由于我国大型零部件加工和装配精度较低,使静压轴承在重载时易产生临界润滑烧瓦的情况。通过此研究方法可以得到不同转速下油膜临界载荷参数,为提高大尺度静压支承运行的可靠性,防止发生零润滑状态即临界润滑导致支承失效提供了理论依据。
发明内容
本发明是要解决在工作台逆时针旋转时,油膜受到剪切流与压差流的影响,从而产生油膜润滑为零,导致润滑失效的情况。本发明通过理论推导,计算出大尺度静压支承油膜承载载荷的理论值。通过前处理软件,得到静压支承油膜模型,利用Fluent软件对理论计算的临界载荷值进行验证,大尺度静压支承临界载荷参数获取按以下步骤实现:
步骤A、理论推导大尺度静压支承油膜承载能力表达式如下:
(1)
式中 w为工作台承载载荷,为工作台旋转角速度,为流体动力粘度,为油膜厚度,为静压支承结构参数,A e 为有效承载面积;
步骤B、基于推导出的理论公式,得出工作台转速与油腔平均压力的临界关系,表1列出了当转速分别为0.5r/min、1r/min、1.5r/min、2r/min、2.5r/min、3r/min、3.5r/min、4r/min、4.5r/min、5r/min时对应的油腔压力及理论临界外载荷值;
表1 静压支承临界润滑时载荷理论值
转速(r/min) | 0.5 | 1 | 1.5 | 2 | 2.5 | 3 | 3.5 | 4 | 4.5 | 5 |
油腔压力(MPa) | 0.05 | 0.11 | 0.16 | 0.22 | 0.28 | 0.33 | 0.39 | 0.44 | 0.50 | 0.55 |
理论临界载荷W(t) | —— | —— | 14.41 | 42.95 | 65.48 | 89.02 | 115.55 | 139.09 | 162.62 | —— |
步骤C、利用前处理软件,对静压支承油膜建立三维仿真模型,并对模型进行网格划分,最后将前处理结果保存chengzai. mesh文件;
步骤D、应用流体计算软件Fluent,对静压支承油膜进行仿真分析,最终得到大尺度静压支承的仿真结果;
步骤E、根据大尺度静压支承的仿真结果,通过此数值方法方法验证大尺度静压支承油膜载荷临界状态数学模型的正确性。
发明效果
本发明方法主要是基于推导出的大尺度静压支承油膜承载能力的理论公式,计算出静压支承临界润滑时载荷的理论值。根据计算结果得出当工作台转速等于4.5r/min时,理论计算得出的静压支承临界载荷值已高于支承本身的最大承载值。当工作台转速等于0.5r/min和1r/min时,通过理论推导出的公式计算出对应临界润滑状态下的油腔压力小于空载下油腔压力值,故此时静压支承不会发生临界润滑。而当工作台转速处于大于等于4.5r/min时,在此转速范围内静压支承将不会发生临界润滑现象;当工作台转速小于等于1r/min时,在此转速范围的静压支承将不会发生临界润滑现象。最后再通过Fluent软件对计算结果进行验证,得出各工况下对应油膜流体域流动情况,即对应的截面face1的流量情况。通过仿真结果得出流量值均从负方向趋近零,而剪切流量稍大于压差流量,理论计算得到的临界载荷值稍小一些,与模拟对应的差值主要来源于理论上对模型结构尺寸简化带来的误差。说明此发明可以得到大尺度静压支承间隙油膜处于临界润滑状态时,进行对应载荷参数,这对有效预防大尺度静压支承润滑失效,避免发生临界润滑产生烧瓦情况有很大的实用意义。
附图说明
图1是本发明流程图,图2是本发明单个油垫封油边处face1位置图,图3是本发明步骤D的具体流程图,图4是本发明在不同转速与对应载荷下face1的临界流量仿真结果。
具体实施方式
大尺度静压支承临界载荷参数获取按以下步骤实现:
步骤A、理论推导大尺度静压支承油膜承载能力表达式如下:
(1)
式中 w为工作台承载载荷,为工作台旋转角速度,为流体动力粘度,为油膜厚度,为静压支承结构参数,A e 为有效承载面积;
步骤B、基于推导出的理论公式,得出工作台转速与油腔平均压力的临界关系,表A列出了当转速分别为0.5r/min、1r/min、1.5r/min、2r/min、2.5r/min、3r/min、3.5r/min、4r/min、4.5r/min、5r/min时对应的油腔压力及理论临界外载荷值;
表A 静压支承临界润滑时载荷理论值
转速(r/min) | 0.5 | 1 | 1.5 | 2 | 2.5 | 3 | 3.5 | 4 | 4.5 | 5 |
油腔压力(MPa) | 0.05 | 0.11 | 0.16 | 0.22 | 0.28 | 0.33 | 0.39 | 0.44 | 0.50 | 0.55 |
理论临界载荷W(t) | —— | —— | 14.41 | 42.95 | 65.48 | 89.02 | 115.55 | 139.09 | 162.62 | —— |
步骤C、利用前处理软件,对静压支承油膜建立三维仿真模型,并对模型进行网格划分,最后将前处理结果保存chengzai. mesh文件;
步骤D、应用流体计算软件Fluent,对静压支承油膜进行仿真分析,最终得到大尺度静压支承的仿真结果;
步骤D1、运行仿真软件Fluent,读入网格文件chengzai. mesh文件,检查并显示网格;
步骤D2、定义求解设置,Define-Models-Solver,单击Steady采用稳态计算,其他选项采用默认值;
步骤D3、流体物理属性设置,单击Define-Materials设置流体密度和粘度等油膜参数的;
步骤D4、边界条件的设置,点击Define-Boundary Conditions选项,设置质量流入口,压力出口,周期边界和旋转壁面的转速;
步骤D5、初始化流场,点击Solve-Initialize-Init并保存case文件;
步骤D6、监视平面设置,点击Surface-Plane选取三个坐标点确定监视平面同时命名为face1;
步骤D7、监视器设置,点击Solve-Monitors-Residual参数为默认参数并采用质量流量进行监视;
步骤D8、进行800步迭代计算,点击Iterate进行迭代;
步骤E、根据大尺度静压支承的仿真结果,通过此数值方法方法验证大尺度静压支承油膜载荷临界状态数学模型的正确性。
本实施方式效果
本发明方法主要是基于推导出的大尺度静压支承油膜承载能力的理论公式,计算出静压支承临界润滑时载荷的理论值。根据计算结果得出当工作台转速等于4.5r/min时,理论计算得出的静压支承临界载荷值已高于支承本身的最大承载值。当工作台转速等于0.5r/min和1r/min时,通过理论推导出的公式计算出对应临界润滑状态下的油腔压力小于空载下油腔压力值,故此时静压支承不会发生临界润滑。而当工作台转速处于大于等于4.5r/min时,在此转速范围内静压支承将不会发生临界润滑现象;当工作台转速小于等于1r/min时,在此转速范围的静压支承将不会发生临界润滑现象。最后再通过Fluent软件对计算结果进行验证,得出各工况下对应油膜流体域流动情况,即对应的截面face1的流量情况。通过仿真结果得出流量值均从负方向趋近零,而剪切流量稍大于压差流量,理论计算得到的临界载荷值稍小一些,与模拟对应的差值主要来源于理论上对模型结构尺寸简化带来的误差。说明此发明可以得到大尺度静压支承间隙油膜处于临界润滑状态时,进行对应载荷参数,这对有效预防大尺度静压支承润滑失效,避免发生临界润滑产生烧瓦情况有很大的实用意义。
Claims (6)
1.大尺度静压支承临界载荷参数研究,其特征在于首先推导出的大尺度静压支承油膜承载能力的理论公式,计算出静压支承临界润滑时载荷的理论值;再利用流体计算软件Fluent得出仿真结果,进而验证大尺度静压支承油膜载荷临界状态数学模型的正确性,大尺度静压支承临界载荷参数获取按以下步骤实现:
步骤A、理论推导大尺度静压支承油膜承载能力表达式如下: (1)
式中 w为工作台承载载荷,为工作台旋转角速度,为流体动力粘度,为油膜厚度,为静压支承结构参数,Ae为有效承载面积;
步骤B、基于推导出的理论公式,得出工作台转速与油腔平均压力的临界关系,表1列出了当转速分别为0.5r/min、1r/min、1.5r/min、2r/min、2.5r/min、3r/min、3.5r/min、4r/min、4.5r/min、5r/min时对应的油腔压力及理论临界外载荷值;
表1静压支承临界润滑时载荷理论值
步骤C、利用前处理软件,对静压支承油膜建立三维仿真模型,并对模型进行网格划分,最后将前处理结果保存chengzai. mesh文件;
步骤D、应用流体计算软件FLUENT,对静压支承油膜进行仿真分析,最终得到大尺度静压支承的仿真结果;
步骤E、根据大尺度静压支承的仿真结果,通过此数值方法方法验证大尺度静压支承油膜载荷临界状态数学模型的正确性。
2.根据权利要求1所述的大尺度静压支承临界载荷参数研究,其特征在于步骤A中通过理论推导,得出大尺度静压支承油膜承载载荷的表达式。
3.根据权利要求1所述的大尺度静压支承临界载荷参数研究,其特征在于所述步骤B中根据推导出的理论公式,计算出相应的油腔压力和理论临界外载荷值。
4.根据权利要求1所述的大尺度静压支承临界载荷参数研究,其特征在于所述步骤C中利用前处理软件,对静压支承油膜建立三维仿真模型,并对模型进行网格划分,最后将前处理结果保存chengzai. mesh文件。
5.根据权利要求1所述的大尺度静压支承临界载荷参数研究,其特征在于所述步骤D中运用Fluent软件,对静压支承油膜进行仿真分析,最终得到大尺度静压支承的仿真结果,其操作过程具体为:
步骤D1、运行仿真软件Fluent,读入网格文件chengzai. mesh文件;
步骤D2、定义求解设置;
步骤D3、流体物理属性设置;
步骤D4、边界条件的设置;
步骤D5、初始化流场;
步骤D6、监视平面设置;
步骤D7、监视器设置;
步骤D8、进行800步迭代计算。
6.根据权利要求1所述的大尺度静压支承临界载荷参数研究,其特征在于所述步骤E中根据大尺度静压支承的仿真结果,通过此数值方法方法验证大尺度静压支承油膜载荷临界状态数学模型的正确性。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410710729.8A CN104598666B (zh) | 2014-12-01 | 2014-12-01 | 大尺度静压支承临界载荷参数获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410710729.8A CN104598666B (zh) | 2014-12-01 | 2014-12-01 | 大尺度静压支承临界载荷参数获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104598666A true CN104598666A (zh) | 2015-05-06 |
CN104598666B CN104598666B (zh) | 2017-12-08 |
Family
ID=53124446
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410710729.8A Expired - Fee Related CN104598666B (zh) | 2014-12-01 | 2014-12-01 | 大尺度静压支承临界载荷参数获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104598666B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104615801A (zh) * | 2014-12-03 | 2015-05-13 | 哈尔滨理工大学 | 重型静压轴承临界润滑状态转速值确定方法 |
CN106870911A (zh) * | 2016-08-31 | 2017-06-20 | 沈阳鼓风机集团压力容器有限公司 | 用于鼓风机及压缩机辅机的润滑油系统的设计方法及装置 |
CN110017988A (zh) * | 2019-04-23 | 2019-07-16 | 哈尔滨理工大学 | 径向倾斜对静动压混合油垫综合润滑性能影响的研究方法 |
CN113158369A (zh) * | 2021-04-19 | 2021-07-23 | 哈尔滨理工大学 | 静压推力轴承油垫封油边的油膜流量模拟监测方法 |
-
2014
- 2014-12-01 CN CN201410710729.8A patent/CN104598666B/zh not_active Expired - Fee Related
Non-Patent Citations (2)
Title |
---|
于晓东: "重型静压推力轴承力学性能及油膜态数值模拟研究", 《 中国博士学位论文全文数据库》 * |
张艳芹: "重型静压推力轴承综合物理场研究", 《 中国博士学位论文全文数据库》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104615801A (zh) * | 2014-12-03 | 2015-05-13 | 哈尔滨理工大学 | 重型静压轴承临界润滑状态转速值确定方法 |
CN104615801B (zh) * | 2014-12-03 | 2018-01-19 | 哈尔滨理工大学 | 重型静压轴承临界润滑状态转速值确定方法 |
CN106870911A (zh) * | 2016-08-31 | 2017-06-20 | 沈阳鼓风机集团压力容器有限公司 | 用于鼓风机及压缩机辅机的润滑油系统的设计方法及装置 |
CN106870911B (zh) * | 2016-08-31 | 2019-09-17 | 沈阳鼓风机集团压力容器有限公司 | 用于鼓风机及压缩机辅机的润滑油系统的设计方法及装置 |
CN110017988A (zh) * | 2019-04-23 | 2019-07-16 | 哈尔滨理工大学 | 径向倾斜对静动压混合油垫综合润滑性能影响的研究方法 |
CN113158369A (zh) * | 2021-04-19 | 2021-07-23 | 哈尔滨理工大学 | 静压推力轴承油垫封油边的油膜流量模拟监测方法 |
CN113158369B (zh) * | 2021-04-19 | 2023-11-28 | 哈尔滨理工大学 | 静压推力轴承油垫封油边的油膜流量模拟监测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104598666B (zh) | 2017-12-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104598666A (zh) | 大尺度静压支承临界载荷参数获取方法 | |
Guan et al. | Theoretical study of flow ripple for an aviation axial-piston pump with damping holes in the valve plate | |
CN107341288B (zh) | 一种通过调整轴承标高控制联合循环机组振动的优化方法 | |
Zhou et al. | A rotor unbalance response based approach to the identification of the closed-loop stiffness and damping coefficients of active magnetic bearings | |
CN102141084B (zh) | 基于膜厚可变的静压推力轴承间隙油膜温度与厚度关系数值模拟方法 | |
CN103632006B (zh) | 一种多级离心泵湿转子临界转速的确定方法 | |
CN105608297A (zh) | 一种角接触球轴承动刚度的数值模拟方法 | |
CN104268342A (zh) | 一种基于有限元的轴承保持架的振动特性分析方法 | |
CN103246786B (zh) | 一种建立静压推力轴承热力耦合变形后润滑油膜模型的方法 | |
CN105527955B (zh) | 一种飞机质量特性建模方法 | |
CN105893695A (zh) | 一种合缸状态下小型汽轮机通流间隙误差补偿方法 | |
CN104615801B (zh) | 重型静压轴承临界润滑状态转速值确定方法 | |
CN111753380A (zh) | 一种火箭发动机液氧涡轮泵的间隙密封模型的建模方法 | |
CN104091021B (zh) | 一种计算倾斜时油垫承载力的方法 | |
CN103226635A (zh) | 基于三维动网格的旋转叶轮机械非定常流场的计算方法 | |
Larsen et al. | On the numerical simulation of nonlinear transient behavior of compliant air foil bearings | |
CN107526914B (zh) | 基于结构化动网格的可倾瓦滑动轴承变流域流场计算方法 | |
Li et al. | Application of computational fluid dynamics and fluid structure interaction techniques for calculating the 3D transient flow of journal bearings coupled with rotor systems | |
CN105095583B (zh) | 一种微尺度下的静压主轴模态分析方法 | |
CN104537189A (zh) | 一种静压转台运动误差建模及计算方法 | |
CN104331539B (zh) | 一种核电站管道热分层效应疲劳评价方法及系统 | |
Temis et al. | Finger seal design based on fluid-solid interaction model | |
Nowald et al. | Influence of circumferential grooves on the non-linear oscillations of turbocharger rotors in floating ring bearings | |
CN103150434A (zh) | 减振器环形阀片复合应力的计算方法 | |
Guo et al. | The FEA contact analysis of high pressure packers |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171208 Termination date: 20181201 |
|
CF01 | Termination of patent right due to non-payment of annual fee |