磨床电主轴动静压轴承静动特性的理论研究＊

王 苗 门 静 韩少燕 贾 谦 李 晓

（①西安交通大学城市学院机械工程系，陕西 西安 710018；②机器人与智能制造陕西省高校工程研究中心，陕西 西安 710018；③西安景辉信息科技有限公司，陕西 西安 710049）

制造业是国民经济的基础和基石，发展装备制造业已成为世界各国制造科技领域中抢占制高点的突破口[1]。研发先进制造设备对提升国家制造业水平有着非常重要的意义，随着社会的发展，数控机床作为制造业的基础，目前的发展趋势为高速化、智能化、高精密化和绿色化[2]。美国、德国和日本等国家在精密机床领域一直处于领先地位，且关键技术目前处于对外封锁状态，我国在机床技术的自主研究方面处于快速发展阶段[3]。电主轴在机床中是一整套组件，是由电主轴本身及其附件组成的系统工程[4]。高速电主轴可以获得特殊要求的加工精度和表面质量，在机械加工装备领域广泛应用，已成为金属切削的主要技术[5]。

针对电主轴的结构设计，刘洋等[6]利用优化设计理论，建立了机床主轴结构设计数学模型，获得了最优主轴结构设计方案，使得机床主轴在满足约束条件下，实现节省材料且提高刚度的目的。金翔等[7]构建了电主轴刚度计算方法，并设计了一种主轴动刚度测试方法，利用动刚度测试平台研究了高速电主轴支撑刚度的变化规律。贾谦等[8]针对精密铣磨床的电主轴的使用需求，分析了动静压轴承的基本结构和使用工况，建立了动静压轴承静动特性参数的计算模型，静特性参数包括温升、流量、承载力和功耗，动特性参数为刚度。熊万里等[9]阐述了气体悬浮电主轴的现状和问题，探讨了轴承类型、参数等对电主轴刚度的影响。范丽婷等[10]采用层次分析法确定电主轴运行状态各影响因素的权重，表征各因素的影响程度，并结合综合指数对电主轴运行状态进行定量评价。郭维祺等[11]提出以电主轴稳定性评价参数为优化目标，以电主轴尺寸为设计变量，来实现电主轴稳定性的优化设计。

本文针对精密机床电主轴，采用动静压滑动轴承作为径向支撑，通过研究从车身与底盘合装自动拧紧工艺效率的提升角度出发，引入超声检测技术进行车身底盘合装螺栓预紧力的超声检测技术研究静动特性分析。

1 电主轴动静压轴承的结构和使用工况

精密机床在航空航天领域中起着重要的作用，例如精密铣和磨床可高效、高精密地加工航空薄壁件、光学玻璃等结构件和精密观测仪器的配件。主轴系统的回转精度和刚度直接决定加工件的制造质量，因此电主轴支撑轴承的设计显得尤为重要。本文研究的精密机床加工精度在微米量级，图1a 所示为电主轴的内部布局示意图，可以看出电主轴主要由转子、定子、前轴承、后轴承及壳体等组成。电主轴支撑轴承的静动特性对其使用精度有重要的影响，较大温升或较小的支撑刚度都会降低加工精度。电主轴的支撑轴承种类繁多，主要有滚动轴承、液体滑动轴承及气体滑动轴承，本文研究的电主轴采用动静压轴承作为支撑，动静压轴承润滑系统的布局如图1b 所示。

图1 电主轴的内部布局

本文研究的高速高精密磨床电主轴基本参数指标见表1，可以看出因为是磨床使用的，所以该电主轴的功率不算很高，但对转速和回转精度要求较高。

表1 磨床电主轴基本参数指标

本文采用的动静压轴承在结构设计上结合了动压轴承和静压轴承的优点，从其结构上看是在动压滑动轴承的基体上添加了4 个静压油腔和4 个节流孔。动静压轴承的结构如图2 所示。动静压轴承所具有的显著特点包括支撑精度高、刚度大、吸振抗震性能好等。将动静压轴承应用于精密机床主轴有望增加机床在高速切削或磨削下的主轴回转精度。动静压轴承的主要结构参数件见表2。

表2 动静压轴承的结构参数

图2 动静压轴承的结构示意图

2 动静压轴承静动特性的理论计算方法

动静压轴承的二维模型如图3 所示，e表示偏心量，h表示油膜间隙，θ为以竖直方向为轴，偏心线所对应的角度，即偏位角，φ为自连心线算起的位置角。

图3 动静压轴承的二维模型图

本文在动静压轴承结构设计和性能的计算中，采用无限宽和无限窄轴承理论来对实际轴承形式进行简化，纳维-斯托克斯方程表示惯性力、压力和摩擦力黏滞力之间的平衡关系，是研究黏滞流体运动的基本方程，但无法求得通解，因此在工程中须进行简化假设：（1）计算时假设液体为层流，不考虑紊流等流态。与黏性力相比，质量力忽略不计；与剪切力相比流体的惯性力可省略。（2）与其他尺寸相比，假设流体膜厚度是小量，可不计流体膜的曲率，并把转动速度用平均速度替代。在流体膜的任何位置上，沿膜厚方向的黏度、压力和密度不变。（3）假设轴承表面无相对滑动，流体膜流动时，仅沿膜厚方向的压力是大量，其余均可忽略不计。（4）对于本文研究的动静压轴承，又做了不考虑转轴挠曲影响的假设。

由流动的连续性条件，得Reynolds 方程的一般形式为

式（1）右边的三项依次代表楔形间隙、切向速度的变化和表面的法向接近对油膜压力承载所起的作用。对于具有不可压缩流体和稳态载荷的动静压轴承，如果不考虑温度对粘度和密度的影响，表示压力分布的二维流动雷诺方程为

求轴承阻尼系数时使用的动态Reynolds 方程如下：

Reynolds 方程的边界条件为

式中：U为轴颈的线速度；η为润滑介质动力粘度；φin为轴承的进液边界；φout是出液边界；Pa为环境压力，这里取0 MPa；Ps为轴承供油压力，这里表示轴承的内外压力差；B为轴承的宽度。

对于圆柱轴承，油膜厚度h沿圆周方向变化，如图3 所示，其间隙函数可表示为

式中：h0为轴承半径间隙；φ为自连心线算起的位置角。

根据流量连续性原理，从轴承的一个油腔向轴向和周向流出的润滑油的流量等于从这个油腔处的节流器流进的润滑油的流量。以一个油腔为例，如图4 所示，由于该轴承结构和尺寸的对称性，只需取半个液腔的流量进行研究。根据流量连续方程，得单位时间内流经控制体的体积流量增量为

图4 轴承油腔流量示意图

沿方向φ和z单位宽度的体积流量分别为

本文的动静压轴承选用的是毛细管节流方式，采用毛细管节流器时，经节流器流入控制体的体积流量为

对于本文所研究的动静压轴承，取图3 所示的坐标系来研究。轴颈旋转时液膜在微弧Rdφ上作用的压力为p，其沿坐标轴x和y的分量分别为pRsinφdφ和pRcosφdφ，沿压力作用区积分，可得两方向承载力分量：

则总的承载量为

载荷角为

在轴承受到小扰动时，油膜力的增量与轴和轴承之间相对运动的关系可表示为

矩阵对角线上的动特性系数分别称为主刚度、主阻尼和正质量系数，非对角线上的系数则被称作交叉刚度、交叉阻尼和交叉质量系数。工程中，可以用4 个刚度系数和4 个阻尼系数表示轴承液膜力的扰动的线性模型，即：

本文动静压轴承静动特性的计算流程如图5 所示。

图5 轴承静动特性的计算流程图

3 动静压轴承静动态性能的计算结果分析

3.1 静态性能的计算结果分析

改变动静压轴承毛细管节流器的节流孔径d，供油压力ps为5 MPa、转速n为10 000 r/min、偏心率ε为0.4，润滑油为2#主轴油，利用图5 中的计算流程对动静压轴承的静、动特性进行了计算，计算结果如图6 所示。图6a 所示为轴承承载力F随d的变化情况，可以看出，随着d的增大，F先增大后减小，在d为0.5 mm 时达到最大值，可见d为0.5 mm 时静压作用明显，且与动压效应叠加获得最大承载力，当d继续增大时静压效应减弱，承载力下降。图6b 所示为轴承流量Q随d的变化曲线，可以看出Q随d的大而增大，且增大的趋势逐渐变缓。Q增大趋势变缓的原因在于动压效应大于静压效应，轴承的膜厚等参数受d的影响越来越小。

图6 轴承静态特性随节流孔径变化图

油腔的存在使得动静压轴承兼具了动压和静压效应。油腔的面积减为零，动静压轴承则蜕变为动压轴承。反之，当油腔面积参数等于整个表面积，轴承则变为了静压轴承。所以，油腔面积参数对动静压轴承的性能有直接影响。这里改变油腔面积A，观察轴承承载能力和动特性系数相应的变化，以对油腔面积参数的选取进行指导。轴承节流孔径d为0.5 mm、供油压力ps为5 MPa、转速n为10 000 r/min、偏心率ε为0.4，润滑油为2#主轴油，计算结果如图7 所示。图7a 和图7b 是轴承静特性随油腔面积A的变化曲线。从图7a 可以看出，承载力F随着A的增加呈线性增加，当A大于1 000 mm2之后F开始减小，因为过大的油腔面积会破坏动压油膜的完成性，所以降低了轴承的承载力。从图7b 可以看出，流量Q随着A增大而增大，但增大趋势有递增迹象。因为油腔面积越大进入封油面的润滑油越多，所以随着A增大会增加Q值。可知，在增大Q的前提之下，适当增加A有利于提高承载能力F。

图7 轴承静态特性随油腔面积变化图

3.2 动态特性的计算结果分析

图8a 和图8b 是轴承刚度、阻尼系数随着节流孔径d的变化曲线。从图8a 中可以看出，主刚度kxx先增大后减小，也在d为0.5 mm 时达到最大值。kyy随着直径增大而减小，但在d为0.5 mm 之后的减小趋势放缓，主刚度的kxx、kyy变化规律与轴承承载力变化规律一致。交叉刚度kxy为负值，但其绝对值大小也随着d增大而减小，且最终趋于一个最小值。交叉刚度kyx随d 着增大而减小，但减小趋势逐渐变缓。图8b 中可以看出，主阻尼cxx和cyy均逐渐减小，最终趋于平缓，变化规律与轴承承载力变化规律一致。而交叉阻尼cxy、cyx很小，变化幅度平缓，可以忽略。

图8 轴承动态特性随节流孔径变化图

图9a 和图9b 是刚度和阻尼系数随A的变化曲线。如图9a 所示，随着A的增加，主刚度kxx、kyy均近似线性增大，kxx在A大于1 000 mm2时达到最大值。交叉刚度kxy随A增加而递减，但在递减趋势逐渐平缓。kxy一直是负数，但其幅值保持递减趋势，几乎与kxy关于x轴对称。图9b 是阻尼系数的变化曲线，可以看出随着A增大，主阻尼和交叉阻尼都呈现递减趋势，油腔面积参数增大，轴承静压性质逐渐占主导地位，对轴承的稳定性有好的作用。

图9 轴承动态特性随油腔面积变化图

4 结语

（1）根据精密磨床的使用要求，采用动静压轴承作为支承轴承，分析了动静压轴承的结构和使用工况。

（2）建立了动静压轴承静动特性的理论分析模型和求解方法，计算了轴承的承载力F、流量Q、刚度系数k和阻尼系数c。

（3）随着d或A的增大，F先增大后减小，Q增大。kxx随着d的增大先增大后减小，kyy随着d增大而减小。随着A的增加，kxx、kyy均增大。