CN106354987B - 一种重载静压转台承载力与油垫温度场分布规律关系计算方法 - Google Patents
一种重载静压转台承载力与油垫温度场分布规律关系计算方法 Download PDFInfo
- Publication number
- CN106354987B CN106354987B CN201610985584.1A CN201610985584A CN106354987B CN 106354987 B CN106354987 B CN 106354987B CN 201610985584 A CN201610985584 A CN 201610985584A CN 106354987 B CN106354987 B CN 106354987B
- Authority
- CN
- China
- Prior art keywords
- equation
- temperature
- oil
- flow
- bearing capacity
- 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
- 230000003068 static effect Effects 0.000 title claims abstract description 31
- 238000004364 calculation method Methods 0.000 title claims abstract description 15
- 230000008859 change Effects 0.000 claims abstract description 30
- 239000007788 liquid Substances 0.000 claims abstract description 27
- 238000000034 method Methods 0.000 claims abstract description 27
- 238000005381 potential energy Methods 0.000 claims abstract description 6
- 238000013178 mathematical model Methods 0.000 claims abstract description 5
- 238000004458 analytical method Methods 0.000 claims abstract description 4
- 230000002706 hydrostatic effect Effects 0.000 claims description 8
- 230000001050 lubricating effect Effects 0.000 claims description 8
- 230000010354 integration Effects 0.000 claims description 5
- 230000004069 differentiation Effects 0.000 claims description 2
- 238000004134 energy conservation Methods 0.000 claims description 2
- 230000004907 flux Effects 0.000 claims description 2
- 230000017525 heat dissipation Effects 0.000 claims description 2
- 239000000314 lubricant Substances 0.000 claims description 2
- 238000005461 lubrication Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 7
- 230000007423 decrease Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 1
- 230000020169 heat generation Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000003754 machining Methods 0.000 description 1
- 230000008569 process Effects 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
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)
- Sliding-Contact Bearings (AREA)
- Magnetic Bearings And Hydrostatic Bearings (AREA)
Abstract
一种重载静压转台承载力与油垫温度场分布规律关系计算方法，属于重载静压转台温度影响分析领域，涉及重型静压转台中单个油垫的温度场数学模型的建立以及数值方法的求解。本方法通过联立雷诺方程和能量方程，建立了极坐标下考虑离心力的温度场数学模型。利用有限差分法求解出油膜温度场的分布，给出了温升随转速的变化规律。模型中考虑了转台温度厂的分布情况，也考虑了离心力的影响。本方法包括三部分，第一部分忽略液体流动时的动能和势能的变化联立雷诺方程和能量方程，第二部分利用有限差分法求解方程，建立不同转速下油垫的温度分布曲线，第三部分分析油垫的温度分布和承载性能随转速的变化规律。最后通过实例论证本发明提出的模型。
Description
技术领域
本发明属于重载静压转台温度影响分析领域，涉及重型静压转台中单个油垫的温度场数学模型的建立以及数值方法的求解。
技术背景
液体静压转台是采用液体静压轴承作为支承，将转台主轴功能与电动机功能从结构上融为一体的功能部件，主要功能是支撑并带动工件实现精密超精密旋转和精密加工。由于运动副之间完全被油膜隔开，所以运动副间的摩擦力大大减小，同时其承载能力、运动精度与寿命却大大提高。正因为液体静压支承的诸多优点，所以其在机床，航空航天，船舶，能源等行业得到了广泛的应用。
静压转台相比与传统转台来说有许多优势，最显著的特点是可以在很高的转速和较低的温度下保持高的可靠性，并且降低对形位公差的要求。油腔中油膜工作性能的优劣直接影响到整个机床运行的可靠性、寿命和经济指标。随着技术的进步，发热问题逐渐成为制约其性能和精度提升的关键因素。国内外学者在液体静压轴承热态性能研究方面做了很多工作，对静压轴承的温度分布研究的比较深入，对于静压转台的研究主要集中在利用Fluent等软件仿真得到静压转台的温度分布情况，并没有完善的理论指导依据。
综上所述，建立一种能够分析静压转台温度变化规律的数值计算方法非常重要。
发明内容
本发明的目的是提供一种油垫的温度分布和承载性能随转速的变化规律的数值方法。该方法首先从雷诺方程和能量方程出发忽略液体流动时的动能和势能的变化，利用有限差分法求解方程，建立不同转速下油垫的温度分布曲线，分析油垫的温度分布和承载性能随转速的变化规律。
本发明是采用以下技术手段实现的：
(1)忽略静压转台自身的变形，定义静压转台为刚体，同时认为在沿润滑膜厚度方向不计压力变化，在沿润滑膜厚度方向粘度数值不变。建立基于雷诺方程的单油垫理论模型。
(2)对于液体润滑，忽略液体流动时的动能和势能变化，这样，液体的能量变化仅是温度的函数。若流动处于稳定状态，那么所有的变量不随时间变化。建立考虑温度和离心力的能量方程。
(3)用有限差分法对雷诺方程和能量方程联立求解，最后积分计算油垫的承载力。之后改变静压转台转速得到不同转速下的油垫压力分布，温度分布和承载力。为了对比引入温粘关系对承载力的影响，分别计算考虑温升的承载力和不考虑温升的承载力。
本方法通过联立雷诺方程和能量方程，建立了极坐标下考虑离心力的温度场数学模型。利用有限差分法求解出油膜温度场的分布，给出了温升随转速的变化规律。
模型中考虑了转台温度厂的分布情况，也考虑了离心力的影响。本方法包括三部分，第一部分忽略液体流动时的动能和势能的变化联立雷诺方程和能量方程，第二部分利用有限差分法求解方程，建立不同转速下油垫的温度分布曲线，第三部分分析油垫的温度分布和承载性能随转速的变化规律。最后通过实例论证本发明提出的模型。
通过下面的描述并结合附图说明，本发明会更加清晰，附图说明用于解释本发明方法及实施案例。
附图说明
图1圆形油垫油膜上任一微元受力图
图2液体热流动图
图3液体流动图
图4温压耦合求解流程图
图5封油边无量纲压力分布图
图6封油边无量纲温度分布图
图7油垫承载力随静压转台转速变化曲线
图8油膜温度随随静压转台转速变化曲线
图9转速为60r/min时油膜最高温度沿油垫周向分布
具体实施方式
本发明实施一种考虑温度和离心力的重型静压转台中单个油垫的数值求解方法，下面结合附图，对本发明的实施进行具体说明。
步骤一：雷诺方程的建立
采用微元体分析方法推导雷诺方程。首先在圆形油垫油膜上任取一微元，其受力如图1所示，根据微元的受力列出微元体的力平衡方程。沿r方向的受力平衡，得：
p——液体压力；
r——油垫半径；
τr——径向切应力；
ρ——液体密度；
ω——转台转速；
——压力在径向的变化率；
——径向切应力在z向的变化率；
方向受力平衡，得：
——周向切应力；
——周向切应力在z向的变化率；
——压力在周向的变化率；
定义润滑剂为牛顿液体，流动形式为层流，油液不可压缩，根据牛顿粘滞定律，
η——液体动力粘度；
将式(3)带入式(1)、式(2)式并略去高阶微小量，对z进行二次积分并代入边界条件：当z＝0时，当z＝h时，代表油膜上表面周向速度，求得：
将带入式(4)，同理对z进行二次积分并代入边界条件：当z＝0时，vr＝Vr；当z＝h时，vr＝0，求得：
不可压缩液体运动的连续性方程为：
将式(6)对z积分，交换微分次序，化简得：
Vr——油膜上表面径向速度；
将式(4)、式(5)代入式(7)，并略去h的高次项，则得到考虑离心力作用的极坐标雷诺方程：
步骤二：能量方程的建立
本位计算，以对流散热为主而忽略膜厚方向的热传导，所以润滑膜温度T只是r和的函数。分析液体在流动中热能和机械功的变化，如图2所示，取夹角为径向长度为dr，高度为h的微元体进行分析。设和qr分别代表微元体周向和径向的容积流量，那么流入微元体的热流量应为Hr＝qrTρc和若取则流入微元体的热量总和为
Hr——径向热流量；
——周向热流量；
用W表示在微元体中所做的机械功，将Hr和代入，那么根据能量守恒原理得到如下关系式：
W——在微元体中所做的机械功；
T——润滑膜温度；
c——液体比热容；
有流量连续条件可知
qr——微元体在径向的容积流量
——微元体在周向的容积流量
在微元体中所做的功包括两部分，流动功和摩擦功。液体流动图如图3所示r方向的流动功为
取dr＝1，略去高阶微量，则沿r方向的流动功为同理可得方向的流动功。
则微元体所做的总的流动功为又因为因而流动总功变为
略去h高次项，径向流量
同理可得周向流量
静压转台油膜下表面速度为0，故只需计算上表面摩擦力所做的功，r方向微元体上表面摩擦力为由于取故代入vr略去h高次项得
r方向微元体中摩擦力做功为
同理方向微元体中摩擦力做功为
微元体所消耗的总功W为
将式(14)、式(15)、式(17)、式(18)代入式(19)得
再将式(12)代入式(20)，经整理求得静压转台油膜的能量方程
步骤三：数值求解
表1静压转台几何参数和油液参数
用有限差分法对雷诺方程(8)和能量方程(21)进行联立迭代求解，直到压强分布误差满足精度要求。
最后积分计算油垫的承载力。给定雷诺方程的边界条件在油腔内在四周边缘上无量纲承载力和无量纲流量的计算公式为
承载力F和流量Q的计算公式为
其中p0为供油压力。由上式求得
将(22)代入(25)得到承载力计算公式
之后改变静压转台转速得到不同转速下的油垫压力分布，温度分布和承载力。具体计算流程如图4所示。
为了对比引入温粘关系对承载力的影响，分别计算考虑温升的承载力和不考虑温升的承载力。图5，图6是油垫封油边无量纲压力和温度分布图，温度沿直径方向封油边向外温度逐渐升高，可以看出压力分布高的地方温升也大。
当不考虑温升变化时静压转台转速由10r/min增加到80r/min时油膜承载力随转速的变化如图7橙色曲线所所示。
计算结果表明，随转速加快油垫承载力都成下降趋势，但考虑温升之后承载力下降的更多，且转速越高这种差异越明显。这是因为转速越高发热量越大，温度升高越多，根据温粘关系式油液粘度降低导致承载力下降。当转速达到80r/min时油垫平均温度升高了11℃，最高温度变化了25.8℃，由于温升导致承载能力下降了18.1％，这说明温升对油膜承载能力有较大影响。由图8可以看出虽然油膜平均温度变化不是很大，但是最高温度对转速很敏感，随转速提高温度上升剧烈。
图9为静压转台转速为60r/min时分别在 和处的最高温度分布。从图中可以看出油垫的温度分布沿周向大致成正弦函数形状，这是由于油垫的中心和的静压转台旋转中心不重合，从而油垫油液的速度分布沿封油边轴向大致成正弦函数。具体的温度分布形状还与静压转台中心到油垫中心的距离、油垫半径、转速有关。
Claims (2)
1.一种重载静压转台承载力与油垫温度场分布规律关系计算方法，其特征在于：本计算方法包括如下步骤：
(1)忽略静压转台自身的变形，定义静压转台为刚体，同时认为在沿润滑膜厚度方向不计压力变化，在沿润滑膜厚度方向粘度数值不变；建立基于雷诺方程的单油垫理论模型；
(2)对于液体润滑，忽略液体流动时的动能和势能变化，这样，液体的能量变化仅是温度的函数；若流动处于稳定状态，那么所有的变量不随时间变化；建立考虑温度和离心力的能量方程；
(3)用有限差分法对雷诺方程和能量方程联立求解，最后积分计算油垫的承载力；之后改变静压转台转速得到不同转速下的油垫压力分布，温度分布和承载力；为了对比引入温粘关系对承载力的影响，分别计算考虑温升的承载力和不考虑温升的承载力；
本方法通过联立雷诺方程和能量方程，建立了极坐标下考虑离心力的温度场数学模型；利用有限差分法求解出油膜温度场的分布，给出了温升随转速的变化规律；
模型中考虑了转台温度场的分布情况，也考虑了离心力的影响；本方法包括三部分，第一部分忽略液体流动时的动能和势能的变化联立雷诺方程和能量方程，第二部分利用有限差分法求解方程，建立不同转速下油垫的温度分布曲线，第三部分分析油垫的温度分布和承载性能随转速的变化规律。
2.根据权利要求1所述的一种重载静压转台承载力与油垫温度场分布规律关系计算方法，其特征在于：实施一种考虑温度和离心力的重型静压转台中单个油垫的数值求解方法，包括如下步骤，
步骤一：雷诺方程的建立
采用微元体分析方法推导雷诺方程；首先在圆形油垫油膜上任取一微元，根据微元的受力列出微元体的力平衡方程；沿r方向的受力平衡，得：
p——液体压力；
r——油垫半径；
τr——径向切应力；
ρ——液体密度；
ω——转台转速；
——压力在径向的变化率；
——径向切应力在z向的变化率；
方向受力平衡，得：
——周向切应力；
——周向切应力在z向的变化率；
——压力在周向的变化率；
定义润滑剂为牛顿液体，流动形式为层流，油液不可压缩，根据牛顿粘滞定律，
η——液体动力粘度；
将式(3)带入式(1)、式(2)式并略去高阶微小量，对z进行二次积分并代入边界条件：当z＝0时，当z＝h时， 代表油膜上表面周向速度，求得：
将带入式(4)，同理对z进行二次积分并代入边界条件：当z＝0时，vr＝Vr；当z＝h时，vr＝0，求得：
不可压缩液体运动的连续性方程为：
将式(6)对z积分，交换微分次序，化简得：
Vr——油膜上表面径向速度；
将式(4)、式(5)代入式(7)，并略去h的高次项，则得到考虑离心力作用的极坐标雷诺方程：
步骤二：能量方程的建立
本位计算，以对流散热为主而忽略膜厚方向的热传导，所以润滑膜温度T只是r和的函数；分析液体在流动中热能和机械功的变化，取夹角为径向长度为dr，高度为h的微元体进行分析；设和qr分别代表微元体周向和径向的容积流量，那么流入微元体的热流量应为Hr＝qrTρc和若取则流入微元体的热量总和为
Hr——径向热流量；
——周向热流量；
用W表示在微元体中所做的机械功，将Hr和代入，那么根据能量守恒原理得到如下关系式：
W——在微元体中所做的机械功；
T——润滑膜温度；
c——液体比热容；
有流量连续条件可知
qr——微元体在径向的容积流量
——微元体在周向的容积流量
在微元体中所做的功包括两部分，流动功和摩擦功；r方向的流动功为：
取dr＝1，略去高阶微量，则沿r方向的流动功为同理可得方向的流动功；
则微元体所做的总的流动功为又因为因而流动总功变为
略去h高次项，径向流量
同理可得周向流量
静压转台油膜下表面速度为0，故只需计算上表面摩擦力所做的功，r方向微元体上表面摩擦力为由于取故代入vr略去h高次项得
r方向微元体中摩擦力做功为
同理方向微元体中摩擦力做功为
微元体所消耗的总功W为
将式(14)、式(15)、式(17)、式(18)代入式(19)得
再将式(12)代入式(20)，经整理求得静压转台油膜的能量方程
步骤三：数值求解
表1静压转台几何参数和油液参数
用有限差分法对雷诺方程(8)和能量方程(21)进行联立迭代求解，直到压强分布误差满足精度要求；
最后积分计算油垫的承载力；给定雷诺方程的边界条件在油腔内在四周边缘上无量纲承载力和无量纲流量的计算公式为
承载力F和流量Q的计算公式为
其中p0为供油压力；由上式求得
将(22)代入(25)得到承载力计算公式
之后改变静压转台转速得到不同转速下的油垫压力分布，温度分布和承载力。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610985584.1A CN106354987B (zh) | 2016-11-09 | 2016-11-09 | 一种重载静压转台承载力与油垫温度场分布规律关系计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610985584.1A CN106354987B (zh) | 2016-11-09 | 2016-11-09 | 一种重载静压转台承载力与油垫温度场分布规律关系计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106354987A CN106354987A (zh) | 2017-01-25 |
| CN106354987B true CN106354987B (zh) | 2019-09-27 |
Family
ID=57862784
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610985584.1A Active CN106354987B (zh) | 2016-11-09 | 2016-11-09 | 一种重载静压转台承载力与油垫温度场分布规律关系计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106354987B (zh) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107169158B (zh) * | 2017-04-09 | 2021-01-01 | 北京工业大学 | 一种基于流固耦合作用的静压滑座工作性能计算方法 |
| CN109829262B (zh) * | 2019-04-04 | 2021-10-19 | 哈尔滨工程大学 | 一种转子-轴承系统非线性动力学分析方法 |
| CN110427666B (zh) * | 2019-07-18 | 2023-12-26 | 上海理工大学 | 一种计算液体动静压球轴承最大油膜承载能力的方法 |
| CN110378018B (zh) * | 2019-07-18 | 2023-12-26 | 上海理工大学 | 一种液体动静压球轴承的稳态承载能力的计算方法 |
| CN110399670B (zh) * | 2019-07-18 | 2023-12-26 | 上海理工大学 | 一种液体动静压球轴承油膜压力计算模型的建立方法 |
| CN110807277B (zh) * | 2019-10-19 | 2024-02-02 | 北京工业大学 | 一种基于静压转台的油垫损伤程度的分析方法 |
| CN110968966A (zh) * | 2019-11-28 | 2020-04-07 | 北京工业大学 | 一种在倾斜状态下确定重型机床静压转台最大合理转速的方法 |
| CN116562184B (zh) * | 2023-04-18 | 2024-02-13 | 中南大学 | 一种热压罐成型工艺中模具温度场快速预测计算方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103186698A (zh) * | 2013-04-12 | 2013-07-03 | 北京工业大学 | 一种重型机床静压转台动、静态性能仿真优化方法 |
| CN104143026A (zh) * | 2014-07-31 | 2014-11-12 | 北京工业大学 | 一种考虑倾斜与热的静压油垫性能计算方法 |
| CN104537189A (zh) * | 2015-01-21 | 2015-04-22 | 北京工业大学 | 一种静压转台运动误差建模及计算方法 |
-
2016
- 2016-11-09 CN CN201610985584.1A patent/CN106354987B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103186698A (zh) * | 2013-04-12 | 2013-07-03 | 北京工业大学 | 一种重型机床静压转台动、静态性能仿真优化方法 |
| CN104143026A (zh) * | 2014-07-31 | 2014-11-12 | 北京工业大学 | 一种考虑倾斜与热的静压油垫性能计算方法 |
| CN104537189A (zh) * | 2015-01-21 | 2015-04-22 | 北京工业大学 | 一种静压转台运动误差建模及计算方法 |
Non-Patent Citations (3)
| Title |
|---|
| 倾斜状态下静压油垫承载性能研究;刘志峰等;《工程力学》;20150531;第32卷(第5期);第208-212、220页 * |
| 数控立车静压转台油膜温度场仿真分析及优化;衡凤琴等;《机械科学与技术》;20151130;第34卷(第11期);第1733-1737页 * |
| 某重型立车静压转台热特性分析与实验研究;衡凤琴等;《软件技术及应用》;20160131;第29-34页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106354987A (zh) | 2017-01-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106354987B (zh) | 一种重载静压转台承载力与油垫温度场分布规律关系计算方法 | |
| CN105005705B (zh) | 一种改进的圆形静压油垫离心力计算方法 | |
| CN111967107B (zh) | 内反馈动静压滑动轴承油膜压力场的nurbs等几何求解方法 | |
| CN105022900A (zh) | 基于热固耦合分析的重型数控立车静压转台结构优化方法 | |
| CN104776998A (zh) | 一种基于动态刚度系数和阻尼系数的转子轴心轨迹求解方法 | |
| CN109376377B (zh) | 高速加工机床整机结构热力学建模与热设计方法 | |
| CN104200019A (zh) | 一种重型数控落地铣镗床热特性数值仿真方法 | |
| CN104537189B (zh) | 一种静压转台运动误差建模及计算方法 | |
| CN201851656U (zh) | 一种动静压结合机械密封结构 | |
| CN107526914B (zh) | 基于结构化动网格的可倾瓦滑动轴承变流域流场计算方法 | |
| CN108591260A (zh) | 一种考虑动压效应的斜面式双矩形油垫承载力计算方法 | |
| CN207892992U (zh) | 一种扇形环斜面式静压推力轴承油垫 | |
| Huang et al. | Static characteristics of water-lubricated hydrodynamic spiral-groove journal and thrust bearings for motorized spindle | |
| CN110378018B (zh) | 一种液体动静压球轴承的稳态承载能力的计算方法 | |
| CN105069264B (zh) | 双矩形腔静压推力轴承旋转速度与承载合理匹配方法 | |
| He et al. | Investigation on heat transfer performance of heat pipe grinding wheel in dry grinding | |
| CN103292975A (zh) | 油腔流场观测装置 | |
| CN108971528B (zh) | 一种机床主轴利用滑动轴承自动定心的方法及其装置 | |
| CN102779112A (zh) | 干气密封拟合曲线槽型线的生成方法 | |
| Jiang et al. | Research on thermal characteristics of spindle system with tilting pad bearing | |
| CN105317840B (zh) | 一种重型静压轴承热油携带判定方法 | |
| CN108959830B (zh) | 高速动压润滑精密主轴热力学建模与热设计方法 | |
| Zhang et al. | Application of dynamic mesh technology in the oil film flow simulation for hydrostatic bearing | |
| Yang et al. | A Review of Hydrostatic Bearing: Research and Analysis | |
| Zhang et al. | Visual comparative analysis for the oil-air two-phase flow of an oil-jet lubricated roller-sliding bearing |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | 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 |