CN105005705B - 一种改进的圆形静压油垫离心力计算方法 - Google Patents
一种改进的圆形静压油垫离心力计算方法 Download PDFInfo
- Publication number
- CN105005705B CN105005705B CN201510463344.0A CN201510463344A CN105005705B CN 105005705 B CN105005705 B CN 105005705B CN 201510463344 A CN201510463344 A CN 201510463344A CN 105005705 B CN105005705 B CN 105005705B
- Authority
- CN
- China
- Prior art keywords
- mrow
- lubricating pad
- formula
- oil
- msubsup
- 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
Landscapes
- Lubricants (AREA)
- Magnetic Bearings And Hydrostatic Bearings (AREA)
Abstract
一种改进的圆形静压油垫离心力计算方法，在该方法中先推导改进的雷诺方程，再计算油垫油膜的压力分布，最后积分计算油垫的承载力。该方法的主要特点是推导改进的雷诺方程时考虑了油垫中心与转台回转中心不重合影响，在计算油垫油膜压力分布时不仅计算了封油边的压力分布而且计算了油垫油腔内部的压力分布，最终得到的油垫的承载力公式能更准确的计算油垫中心与转台回转中心不重合时离心力对油垫承载力的影响。该方法得到的油垫压力分布及承载力公式对油垫油膜温度分布的计算及转台的提速设计有积极作用。
Description
技术领域
本发明属于静压转台设计分析领域，涉及一种圆形静压油垫离心力计算方法，更具体是一种既考虑了转台转动对于油腔内部压力分布的影响也考虑了油垫中心与转台回转中心不重合情况的转台圆形支承油垫离心作用力的计算方法。
背景技术
静压转台(Hydrostatic Rotary Table)是重型高档数控机床的关键配件，其承载原理是用有压力的流体使有相对运动的两个表面分开并借助流体静压来承载。由于运动副之间完全被油膜隔开，所以运动副间的摩擦力大大减小，同时其承载能力、运动精度与寿命却大大提高。静压转台的支承系统由支承油垫，预紧油垫及径向轴承组成。预紧油垫一般为环形油垫而支承油垫一般为扇形和圆形。对于环形预紧油垫及扇形支承油垫而言它们均属于油垫中心与转台回转中心重合的情况，而对于圆形支承油垫则属于油垫中心与转台回转中心不重合的情况。当前各类设计手册中只有计算油垫中心与转台回转中心重合时的离心力，并且在计算过程中认为油腔内部压力处处相等并作为一个边界条件来求解封油边的压力分布，而事实上转台转动后支承油垫不仅其封油边的压力分布会发生改变而且油腔内部的压力分布也会发生改变，进而影响油垫的承载能力。在工程中通常认为离心力对于油垫承载能力影响很小所以在转台的设计过程中大多没有专门计算离心力的影响，而是在设计过程中取一定的设计余量来保证转台的承载性能。然而不论离心力对于油垫承载能力影响有多小，总需要一个定量的指标来精确衡量它的大小，而且转台转动过程中还会产生大量的摩擦热，这一方面影响油垫的承载另一方面会使转台产生热变形进而影响转台的加工精度，而要准确的计算油垫的摩擦热首先就要计算油垫的压力分布。这也是本发明的原因及意义所在。
发明内容
本发明旨在提供一种圆形油垫离心作用力的计算方法。该方法的主要特点是在计算油垫离心力的过程中充分考虑了油垫偏置安装的影响，并计算了油腔内部压力分布的变化。
本发明是采用以下技术手段实现的：
1、将油垫上表面的转速分解为周向和径向速度并将其作为速度边界条件代入简化的N-S方程中得到改进的雷诺方程。
2、对于油腔内部，认为油腔中心部位的压力等于油垫进油口压力并且此处油膜压力分布函数连续可导，以此为边界条件求解改进的雷诺方程得到油腔内部的压力分布函数。同样认为封油边油膜压力与油腔油液压力连续且封油边外侧压力为零，再次求解修正的雷诺方程得到封油边的压力分布函数。
3、对于油腔内部及封油边内部的压力分布函数积分得到油垫的承载力函数。函数中有转台转速的就是离心力项。
本发明的特点在于在推导修正雷诺方程时考虑了转台回转中心与油垫中心不重合的情况，求解油垫的压力分布时考虑了油腔内部压力分布的变化。下面将结合附图更清晰的阐明本发明的方法及实施例。
附图说明
图1静压转台结构简图
图2支承油垫位置简图
图3支承油垫结构简图
具体实施方式
步骤(1)，改进雷诺方程的建立
定义油垫的润滑符合薄膜润滑理论，油膜厚度大小远小于油膜其他方向的尺寸；油液的流态为层流，油液属于牛顿流体不可压缩且其粘度是一个常量；油液本身质量及惯性力的影响可忽略不计；忽略油膜曲率的影响；忽略温度的影响；油膜方向油液的流速远小于其它方向的流速；油膜厚度方向速度的梯度占主导地位；那么极坐标系下的连续性方程与N-S方程可以简化为：
公式(1)为连续性方程，公式(2)是N-S方程，式中r,z为极坐标系的坐标参量，ur,为油液在r,方向的流速，η为油液的粘度，ρ为油液的密度，p为油液压力。对于应用圆形支承油垫的转台而言油垫的中心位置与转台中心位置并不重合而是有一定的偏距，这个偏距为RL(如图1所示)。
图2为静压转台支承油垫与回转中心的示意图，支承油垫的结构简图如图3所示，R1位油垫的内径，R2为油垫的内径，h为油垫的油膜厚度，hc为油腔深度，Q0为油垫的供油量。当转台以Ω的转速转动后支承油垫上方导轨面的线速度为将其向r方向和方向分解得到
所以对方程(3)两边同时对z进行二次积分并代入边界条件：z＝h,z＝0,ur＝0，以得到
将(5)式代入(2)式并代入边界条件：z＝h,z＝0,ur＝0对z进行二次积分得到油膜的径向速度ur：
油膜的径向速度(6式)对油垫区域进行积分就可以得到流出油垫油液的流量。
将(5)式和(6)式代入(1)式并对z在油膜厚度范围内积分就得到了雷诺方程。
当RL＝0时即没有偏心时，公式(8)与当前静压支承书籍中得到的雷诺方程相同，表明该公式的正确性。
步骤(2)油垫压力分布的计算
由于支承油垫与转台的回转中心距离较大，使得油垫处导轨的线速度较大。所以应当考虑转台转动对于油腔压力的影响。对于油腔内部令h＝hc(hc为油腔深度),求解(8)式并带入边界条件存在；可得到油腔内部的压力分布：
上式中p0为油垫进油口压力。对于油腔外部也就是封油边部分，同样求解(8)式并代入边界条件r＝R1,p＝p1(r＝R1)；r＝R2,p＝0油膜的压力分布就可以得到：
认为油液不可压缩则流入油垫的油液的流量可以分为两部分一部分是油泵供给的为Q0，另一部分是上导轨面挤压油膜产生的为流入油垫的油液量应当与流出油垫的油液量相等于是有：
将(7)式与(9)式代入(10)式并化简方程可以得到油垫的进油口压力为：
步骤(3)油垫承载力的计算
在再在油垫表面对式(9)积分就可以得到油垫的承载力为：
上式中支承油垫的承载力可分为三部分，一部分是由泵产生(式中为：)，一部分由导轨对于油垫的挤压效应产生(式中为：)，另一部分由转台的转动产生既离心力项(式中为：)。
Claims (1)
1.一种改进的圆形静压油垫离心力计算方法，在该方法中先推导改进的雷诺方程，再计算油垫油膜的压力分布，最后积分计算油垫的承载力；
其特征在于：其步骤如下，
步骤(1)改进雷诺方程的建立
若油垫的润滑符合薄膜润滑理论，油膜厚度大小远小于油膜其他方向的尺寸；油液的流态为层流，油液属于牛顿流体不可压缩且其粘度是一个常量；油液本身质量及惯性力的影响可忽略不计；忽略油膜曲率的影响；忽略温度的影响；油膜方向油液的流速远小于其它方向的流速；油膜厚度方向速度的梯度占主导地位；那么极坐标系下的连续性方程与N-S方程可以简化为：
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>p</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>z</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
公式(1)为连续性方程，公式(2)是N-S方程，式中r,z为极坐标系的坐标参量，ur,为油液在r,方向的流速，η为油液的粘度，ρ为油液的密度，p为油液压力；对于应用圆形支承油垫的转台而言油垫的中心位置与转台中心位置并不重合而是有一定的偏距，这个偏距为RL；
支承油垫的结构为圆形阶梯结构，R1位油垫的内径，R2为油垫的内径，h为油垫的油膜厚度，hc为油腔深度，Q0为油垫的供油量；当转台以Ω的转速转动后支承油垫上方导轨面的线速度为将其向r方向和方向分解得到
所以对方程(3)两边同时对z进行二次积分并代入边界条件：
z＝h,z＝0,ur＝0，以得到
将(5)式代入(2)式并代入边界条件：z＝h,z＝0,ur＝0对z进行二次积分得到油膜的径向速度ur：
油膜的径向速度对油垫区域进行积分就得到流出油垫油液的流量；
将(5)式和(6)式代入(1)式并对z在油膜厚度范围内积分就得到了雷诺方程；
步骤(2)油垫压力分布的计算
由于支承油垫与转台的回转中心距离较大，使得油垫处导轨的线速度较大；所以应当考虑转台转动对于油腔压力的影响；对于油腔内部令h＝hc，hc为油腔深度，求解(8)式并带入边界条件存在；
得到油腔内部的压力分布：
上式中p0为油垫进油口压力；对于油腔外部也就是封油边部分，同样求解(8)式并代入边界条件r＝R1,p＝p1，r＝R1；r＝R2,p＝0油膜的压力分布就可以得到：
认为油液不可压缩则流入油垫的油液的流量可以分为两部分一部分是油泵供给的为Q0，另一部分是上导轨面挤压油膜产生的为流入油垫的油液量应当与流出油垫的油液量相等于是有：
<mrow>
<mi>Q</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>Q</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<msubsup>
<mi>&pi;R</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
将(7)式与(9)式代入(10)式并化简方程可以得到油垫的进油口压力为：
<mrow>
<msub>
<mi>p</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>6</mn>
<msub>
<mi>&eta;Q</mi>
<mn>0</mn>
</msub>
</mrow>
<mrow>
<msup>
<mi>&pi;h</mi>
<mn>3</mn>
</msup>
</mrow>
</mfrac>
<mi>ln</mi>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mfrac>
<mrow>
<mn>3</mn>
<msup>
<mi>&rho;&Omega;</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>R</mi>
<mi>L</mi>
<mn>2</mn>
</msubsup>
<mi>ln</mi>
<mo>(</mo>
<mfrac>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>)</mo>
<mo>+</mo>
<msubsup>
<mi>R</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
<mn>20</mn>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mn>3</mn>
<mi>&eta;</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>R</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>(</mo>
<mrow>
<msubsup>
<mi>h</mi>
<mi>c</mi>
<mn>3</mn>
</msubsup>
<mo>-</mo>
<msup>
<mi>h</mi>
<mn>3</mn>
</msup>
</mrow>
<mo>)</mo>
<mo>-</mo>
<msubsup>
<mi>R</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>h</mi>
<mi>c</mi>
<mn>3</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mo>&part;</mo>
<mi>h</mi>
</mrow>
<mrow>
<msup>
<mi>h</mi>
<mn>3</mn>
</msup>
<msubsup>
<mi>h</mi>
<mi>c</mi>
<mn>3</mn>
</msubsup>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤(3)油垫承载力的计算
在油垫表面对式(9)积分就得到油垫的承载力为：
上式中支承油垫的承载力可分为三部分，一部分是由泵产生，式中为：一部分由导轨对于油垫的挤压效应产生，式中为：另一部分由转台的转动产生既离心力项，式中为：
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510463344.0A CN105005705B (zh) | 2015-07-31 | 2015-07-31 | 一种改进的圆形静压油垫离心力计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510463344.0A CN105005705B (zh) | 2015-07-31 | 2015-07-31 | 一种改进的圆形静压油垫离心力计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105005705A CN105005705A (zh) | 2015-10-28 |
| CN105005705B true CN105005705B (zh) | 2017-12-01 |
Family
ID=54378375
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510463344.0A Active CN105005705B (zh) | 2015-07-31 | 2015-07-31 | 一种改进的圆形静压油垫离心力计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105005705B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106503391B (zh) * | 2016-11-09 | 2019-10-29 | 北京工业大学 | 一种考虑流固耦合的矩形静压油垫快捷计算方法 |
| CN107061494A (zh) * | 2016-11-09 | 2017-08-18 | 北京工业大学 | 一种用于非对称工件装夹的静压转台支撑垫铁布局方法 |
| CN106439005B (zh) * | 2016-12-07 | 2018-07-31 | 北京工业大学 | 一种带圆弧槽封油边静压转台圆形油垫 |
| CN109977573A (zh) * | 2019-04-01 | 2019-07-05 | 哈尔滨理工大学 | 一种双矩形腔油垫油膜压力场重构方法 |
| CN110069871B (zh) * | 2019-04-26 | 2022-11-04 | 哈尔滨理工大学 | 静压回转工作台偏载承载力计算方法 |
| CN110378018B (zh) * | 2019-07-18 | 2023-12-26 | 上海理工大学 | 一种液体动静压球轴承的稳态承载能力的计算方法 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100980600B1 (ko) * | 2008-06-24 | 2010-09-07 | 이더블유에스 코리아 주식회사 | 유정압 회전테이블 장치 |
| CN104036079B (zh) * | 2014-06-09 | 2017-06-16 | 北京工业大学 | 一种基于灵敏度分析的定量式静压转台优化设计方法 |
| CN104091021B (zh) * | 2014-07-11 | 2017-05-17 | 北京工业大学 | 一种计算倾斜时油垫承载力的方法 |
| CN104143026B (zh) * | 2014-07-31 | 2017-07-14 | 北京工业大学 | 一种考虑倾斜与热的静压油垫性能计算方法 |
-
2015
- 2015-07-31 CN CN201510463344.0A patent/CN105005705B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105005705A (zh) | 2015-10-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105005705B (zh) | 一种改进的圆形静压油垫离心力计算方法 | |
| CN104143026B (zh) | 一种考虑倾斜与热的静压油垫性能计算方法 | |
| Srinivasan | Analysis of static and dynamic load on hydrostatic bearing with variable viscosity and pressure | |
| CN106354987B (zh) | 一种重载静压转台承载力与油垫温度场分布规律关系计算方法 | |
| CN105022900B (zh) | 基于热固耦合分析的重型数控立车静压转台结构优化方法 | |
| CN103823940B (zh) | 一种基于整体动力学模型的静压转台动态响应计算方法 | |
| CN104537189B (zh) | 一种静压转台运动误差建模及计算方法 | |
| CN103909469A (zh) | 一种轧辊磨床砂轮主轴薄膜节流动静压轴承 | |
| CN104091021B (zh) | 一种计算倾斜时油垫承载力的方法 | |
| CN104036079A (zh) | 一种基于灵敏度分析的定量式静压转台优化设计方法 | |
| CN101673318B (zh) | 透平膨胀机径向静压气体轴承最佳静态参数设计方法 | |
| CN113434983A (zh) | 一种滑动轴承转子系统非线性动力学特性快速计算方法 | |
| Yu et al. | Analysis of lubrication characteristics of dynamic-static pressure hybrid thrust bearing considering key factors under eccentric loads | |
| CN110378018B (zh) | 一种液体动静压球轴承的稳态承载能力的计算方法 | |
| CN103292975A (zh) | 油腔流场观测装置 | |
| CN105317840B (zh) | 一种重型静压轴承热油携带判定方法 | |
| CN104933219A (zh) | 一种水润滑推力轴承的全流态模拟参数获取方法 | |
| Zhang et al. | Application of dynamic mesh technology in the oil film flow simulation for hydrostatic bearing | |
| CN206495889U (zh) | 一种动静压推力轴承矩形弧腔扇形油垫 | |
| CN107061494A (zh) | 一种用于非对称工件装夹的静压转台支撑垫铁布局方法 | |
| CN105257701A (zh) | 多油垫式重型静压轴承热油携带因子及热油携带方程 | |
| CN106777637A (zh) | 一种考虑流固耦合影响的静压转台变形的计算方法 | |
| Ni et al. | Study on oil starvation at the oil seal edge of the oil film of hydrostatic bearing with double rectangular cavity | |
| Tóth | CFD Analyses of the Pressure Distribution in Hydrostatic Journal Bearings with Different Recess Shapes | |
| Zhang et al. | Temperature Field of Hydrostatic Supporting disk in Different Viscosity and Rotational Speed |
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 |