-
随着现代科学技术的飞速发展,机床制造精度的提高,热误差在众多的机床误差中所占比重也越来越大,尤其是在精密加工和大件加工中,热变形所引起的加工误差通常会占到总误差的40%~70%[1]。由于重型机床载荷量大、发热量大的特点,热误差占总误差的比例更大。因此,分析重型机床热特性对减小机床热误差,提高机床加工精度至关重要。
文献[2]采用结构优化设计结合温度控制法有效地降低了机床温升,减小了机床热变形。文献[3]分析了滚珠丝杠进给系统的温升和热变形,并使用有限元法分析了丝杠的变形特点。文献[4]首先通过有限元仿真主轴装配体,得到其温升及温升导致的结构变形。实验结果与有限元仿真结果基本一致。文献[5]采用有限元法分析了车削中心主轴箱温度场和变形场,理论分析了热误差组件相对刀具位置变化关系。文献[6]在热分析理论基础上,计算了液体静压轴承内外壁的温度和主轴箱体的温度,分析了主轴系统的稳态温度场和主轴的变形情况。文献[7]对机床主轴箱的主要热源的发热量进行了分析计算,运用有限元方法建立了主轴箱的温度场和变形场,并进行了实验验证。文献[8]通过实验和有限元数值仿真来分析主轴系统不同部件对系统温升和热变形特性的影响程度。文献[9]采用有限元分析方法计算了车削中心主轴箱的温度场和变形场,准确预测主轴箱热特性,有助于提高设计水平。文献[10]通过有限元分析获得主轴箱的热特性,然后采用激光三角测量仪和红外热成像仪测量主轴装配体的温度场和变形场。上述研究都是针对普通数控机床展开的,对于重型机床研究较少,重型机床由于其尺寸大、载荷重、转速慢等特点,存在的问题尚需研究解决。
本文以某型号重型卧式车床为研究对象,对主轴部件及其主轴箱进行了热和结构的耦合分析,并以主轴箱装配体稳态热分析的温度场为基础,计算出主轴的热变形,并讨论了轴承转速、液压油粘度和油膜间隙对主轴温升的影响,为主轴及主轴箱部件的热设计奠定基础,使设计者在设计阶段就可以预测机床主轴箱装配体的温升和热变形情况,并为结构优化设计提供必要的计算依据。
-
对于有限元计算,首先进行零件三维建摸,要求在尽量保持部件的原始结构的基础上,对部件进行适当简化;主轴系统三维模型的简化原则如下:忽略主轴及主轴箱上小的通孔及螺纹孔等;忽略主轴及主轴箱上的倒角和圆角等。
将主轴箱及主轴装配简化模型导入workbench中,采用四面体划分网格法中的Patch Conforming法进行分块网格划分。网格数量的多少将影响计算结果的精度和计算规模的大小。一般而言,网格数量增加,计算精度会有所提高,但同时计算规模也会增加,所以在确定网格数量时应权衡两个因素综合考虑。模型中单元数146 458个,节点数257 950个。建模及分网如图 1所示。
-
主轴系统内部的热源有电机发热、齿轮传动发热、轴承发热等[11]。本文的主轴箱为整体式结构,主轴系统和主传动系统分隔在箱体的两个腔内,便于控制主轴的温升和液压系统的回油,传动齿轮发热大部分热量被润滑油带走,传到主轴系统的热量很少。同时,电机通过皮带轮连接,其发热对机床热变形影响很小[5]。因此主轴静压轴承的发热成为主轴热变形的最主要热源。静压轴承模型如图 2所示。
1) 轴承发热量的计算
液体静压轴承在工作状态时,由于润滑油的粘度,轴和轴承表面的相对运动队油层起剪切作用,因而油液各层间的内摩擦力消耗功率并产生摩擦热,其计算公式为[12]:
$$ {P_f} = \sum {{F_f}} u = 4\eta \left( {\frac{{{A_s}}}{h} + \frac{{{A_r}}}{{h + {h_p}}}} \right){\left( {\frac{{\pi nD}}{{60}}} \right)^2} $$ (1) 式中, ${P_f}$ 为摩擦功耗; ${F_f}$ 为油膜的粘滞阻力;n为轴承转速; $\eta $ 为润滑油的动力粘度;As为油腔封油面面积;Ar为油腔面积;h为油膜平均间隙;hp为油腔深度;D为主轴直径。主轴静压轴承示意图如图 3所示。
由于液压油从静压轴承的油腔中向外流出后,压力降为零。因此,除了供油管路中功率损耗外,油泵的有效输入功率全部损耗在静压油腔中,这些功率全部转换成热能的形式释放[12]:
$$ {P_p} = {p_s}q $$ (2) 式中, ${P_p}$ 为油泵输入功率; ${p_s}$ 为供油压力;q为轴承流量。
静压轴承的热载荷主要是油膜施加在静压轴承上的热流密度。每个轴承与油膜接触面的热流密度为:
$$ q = \frac{N}{S} $$ (3) 式中,N为轴承的发热量(近似等于摩擦功耗与泵功耗之和);S为轴承中油层的总面积。
2) 散热条件分析
主轴箱内外表面与空气和润滑油接触发生对流换热。根据努谢尔特准则,换热系数为:
$$ h = \frac{{{\rm{Nu}}\lambda }}{L} $$ (4) 式中,λ为流体导热系数;Nu为努谢尔特数;L为放热壁定型尺寸。
所以,只要知道努谢尔特数Nu,便可计算换热系数。
① 自然对流散热
主轴箱外表面直接与空气接触,属于自然对流,其准则方程形式可表示为:
$$ {\rm{Nu}} = C{({\rm{Gr}}\Pr )^s} $$ (5) $$ {\rm{Gr}} = \frac{{g\beta {L^3}\Delta t}}{{{v^2}}} $$ (6) 式中,C、n为常数取值,如表 1所示[13];Gr为格拉晓夫准数;Pr为普朗特数;g为重力加速度;β为流体膨胀系数;L为放热壁定型尺寸;Δt为流体和壁面温差;v为流体运动粘度。
表 1 常数C、n的值
位置 适用范围 流态 C n 竖平壁 104≤(GrPr)≤109 层流 0.59 1/4 109≤(GrPr)≤1012 紊流 0.12 1/3 水平板朝上 105≤(GrPr)≤2×107 层流 0.54 1/4 2×107≤(GrPr)≤3×1010 紊流 0.14 1/3 水平板热面朝下 3×105≤(GrPr)≤3×1010 层流 0.27 1/4 ② 旋转轴对流换热
当机床主轴旋转时,会加速主轴表面和空气的对流,同时也加速轴承与润滑油的对流,其准则方程为:
$$ \begin{array}{c} {\rm{Nu}} = 0.133{{\mathop{\rm Re}\nolimits} ^{2/3}}{\Pr ^{1/3}}\\ {\mathop{\rm Re}\nolimits} < 4.3 \times {10^5},{\rm{ }}0.7 < \Pr < 670 \end{array} $$ (7) ③ 强制对流散热
由于机床在运行过程中轴承生热量较大,主轴静压轴承回油,产生热量交换。这属于强制对流换热,其努谢尔特准则方程可表示为[14]:
$$ {\rm{Nu}} = 0.15{\rm{R}}{{\rm{e}}^{0.33}}{\Pr ^{0.43}}{\rm{G}}{{\rm{r}}^{0.1}}{\left( {\frac{{\Pr }}{{{{\Pr }_w}}}} \right)^{0.25}}K(层流) $$ (8) 式中,Prw为壁面温度时的普朗特数;K为校正系数;Re为雷诺数,即:
$$ {\mathop{\rm Re}\nolimits} = \frac{{\omega L}}{v} $$ (9) 式中,ω为流速。
普朗特数为:
$$ \Pr = \frac{{{c_p}\rho \nu }}{\lambda } $$ (10) 式中, ${c_p}$ 为等压比热容; $\rho $ 为密度; $\nu $ 为流体运动粘度; $\lambda $ 为导热系数。
-
机床最高转速140 r/min,当空载状态下主轴转速为100 r/min,温度为22℃时计算主轴系统的稳态温度场与变形场,以及由温度引起的主轴在X、Y、Z方向的热变形和主轴的总变形。
选取油膜厚度为0.06 mm,供油压力为6 MPa,静压油选用46号液压油,密度为870 kg/m3,膨胀系数为0.000 9 k-1,等压比热容为1 890 J/kg∙k,导热系数为0.1,40℃时运动粘度范围41.4~50.6 cSt。
根据上述式(1)~式(2)计算各轴承油膜的热流密度,式(3)~式(10)计算对流换热系数。计算结果前轴承油膜热流密度为1 879.1 W/m2,后轴承油膜热流密度为1 677.8 W/m2;自然对流换热系数为2.7 W/m2∙k,主轴与空气接触部分旋转对流换热系数为10.2 W/m2∙k,静压油腔强制对流换热系数为80.5 W/m2∙k。
主轴箱材料为HT300,主轴材料为42CrMo,静压轴承材料为锡青铜,主轴系统材料物理性能参数详见表 2[12, 15]。
表 2 材料物理性能参数
材料名称 HT300 42CrMo 锡青铜 密度/kg·m3 7 350 7 850 8 800 比热/J×(kg·℃)-1 480 460 343 导热系数/W·(m·℃)-1 52.3 46.5 398 弹性模量/Gpa 135 210 110 热膨胀系数/℃-1×10-6 10.5 12.5 18.0 泊松比 0.25 0.3 0.34 -
达到平衡状态后主轴系统的温度分布如图 4所示,最大温升为11.8℃。剖面图如图 5所示,箱体上端温度较高,这是由于下端箱体与床身接触散热快的缘故。由于主轴箱的非对称性,其温度场也呈现出以主轴为中心发散的状态。
-
影响热变形大小由两个因素是结构的温度场的变化和结构的位移约束边界条件。主轴系统温度场变化由仿真计算得到,机床装配是主轴箱底部安装固定,其余部分均无约束,即位移边界条件可设为主轴系统的安装平面位移为零。
在对主轴箱体进行热-结构耦合分析时,先进行主轴箱体的热分析,然后将热分析的结果作为载荷条件,对主轴箱体的相关位置进行约束,进而通过仿真计算得到主轴箱体的结构变形分析结果。如图 6所示,主轴箱主要考虑影响主轴位置的热变形,主轴中心线在X向最大位移量约为12.3 μm,在Y向位移量约为33.6 μm,主轴在Z向热变形量为132.1 μm。主轴箱体的最大热变形为141.1 μm,发生在主轴端部位置。由此可知,主轴系统主要热变形是由主轴轴向热变形所致。
-
设置与仿真条件一致的实验条件,进行主轴系统温升和热误差测试实验,测量设备如图 7所示。由于温度测点无法布置到静压轴承表面,只是测量了轴承对应位置的主轴箱体表面的温升。测得的最大温升为4.2℃,对应位置仿真结果温升为6.2℃左右,仿真最大温升在轴承处,为11.8℃。实验测得热变形数据如图 8所示,X方向为5 μm,Y方向为30 μm,Z方向为107 μm,这与仿真结果规律一致,在数值上小于仿真结果,这是实验条件与仿真条件的差异造成的,实验环境受到外界影响因素较多,有些散热仿真中没有考虑,所以实验结果略小。可以验证仿真的正确性。
-
主轴轴承发热是机床中最主要的热源,其发热量的大小有许多因素影响。供油压力、轴承转速、油膜间隙、液压油粘度等因素对静压轴承发热量都有影响。因此,分析不同影响因素对主轴系统温度场的影响程度,对有效控制主轴系统温升及热变形有重要意义。
-
通过式(1)可以看到,轴承转速越高,其摩擦功耗就越大,摩擦发热情况就越严重,这将直接导致轴承的温度明显升高。其他条件不变,改变转速,由式(1)~式(10)可计算得到的不同转速下,主轴前后轴承的油膜热流密度及主轴表面与空气间的对流换热系数,如表 3所示。
表 3 不同转速下热流密度和对流系数
转速/r·min-1 热流密度/W·m-2 对流换热系数/W·(m2·k)-1 前轴承 后轴承 20 666.9 597.4 3.5 40 818.4 732.5 5.5 60 1 071.0 957.5 7.2 80 1 424.5 1 272.7 8.8 100 1 879.1 1 677.8 10.2 将不同转速时前后轴承的热流密度和主轴表面的对流系数分别加载到主轴系统的有限元分析模型中,通过软件仿真计算,可以得到在不同转速时主轴系统达到热平衡时的温度场及变形场分布情况。不同转速下主轴中心线Y向偏移量如图 9所示。
由图可知,随着主轴转速的增大,主轴中心线Y向偏移量逐渐从13.5 μm增加到33.6 μm。所以,工作时间允许的前提下,尽量避免主轴高速运转。
-
液压油的粘度通过影响静压轴承油垫粘滞阻力,从而影响轴承发热。同时粘度变化也会影响到油泵的发热功耗。在保证静压轴承足够承载能力和转速的情况下,选取几种不同粘度的液压油计算前后轴承油膜热流密度和液压油对主轴的对流换热系数。其他条件不变,改变液压油牌号,计算结果如表 4所示。
表 4 不同粘度下热流密度
液压油标号 热流密度/W·m -2 前轴承 后轴承 22 1 864.1 1 667.5 32 1 768.9 1 580.6 46 1 879.1 1 677.8 68 2 305.0 2 056.4 由图 10可得,随着液压油粘度的增加,主轴中心线Y向偏移量先是降低,然后逐渐增加。出现该现象的原因是液压油粘度变化导致油膜摩擦功耗和油泵压力功耗同时变化,且变化方向不一致。可见液压油粘度的选取相当重要,偏大偏小都不合适。因此先要找出偏移量最小时对应的液压油粘度,再选择适当的粘度,以降低温升和减小主轴系统热变形。
-
对于静压轴承,油膜间隙越大发热量越小,油膜散热效果越好。但是供油压力会增大,油泵功耗增加。其他条件不变,改变油膜间隙,计算不同油膜间隙时主轴前后轴承热流密度,结果如表 5所示。
表 5 不同间隙下热流密度
油膜间隙/mm 热流密度/W·m-2 前轴承 后轴承 0.04 2 072.0 1 847.6 0.05 1 870.1 1 668.5 0.06 1 879.1 1 677.8 0.07 2 062.5 1 843.0 由图 11可以看出,随着静压轴承油膜间隙的增大,主轴中心线Y向偏移量先是减小后增加,同样是由于间隙变化导致油膜摩擦功耗和油泵功耗变化不一致。所以在设计选取静压轴承时,满足承载要求情况下可以向降低热变形的方向适当调整油膜间隙。
-
在静压轴承转速、静压油粘度、油膜间隙一定的情况下,不同的载荷对应不同的静压轴承供油压力。改变静压轴承供油压力,油泵的功耗也随之改变,从而影响油膜发热。其他条件不变,改变供油压力,计算结果如表 6所示。
表 6 不同供油压力下热流密度
供油压力/Mpa 热流密度/W·m-2 前轴承 后轴承 4 1 536.7 1 370.9 5 1 690.8 1 509.0 6 1 879.1 1 677.8 7 2 101.7 1 877.3 由图 12可知,随着供油压力的增加,主轴中心线的偏移量近似线性增加。由此可知,随着机床载荷的增加,主轴偏移量也是近似线性增加。在载荷与轴承内径尺寸一定的情况下,偏心率越大,供油压力越小,所以在偏心率允许的范围内,尽量减小供油压力以降低温升。
Thermal Characteristics Analysis of the Spindle System of Heavy Duty Horizontal Lathe
-
摘要: 主轴系统的热特性对重型卧式车床的加工精度有重要的影响。以某型号重型卧式车床的主轴系统为主要研究对象,采用有限元热-固耦合方法,仿真计算了主轴系统达到热平衡状态后的温度场分布和热变形特性,分析了热变形导致主轴中心线偏移情况。结合主轴系统静压轴承的结构特点,分析了轴承转速、液压油粘度、油膜间隙和供油压力对主轴系统变形场的影响。结果表明,这些因素对主轴系统的变形场有不同程度的影响,为主轴系统优化设计和热误差补偿提供了理论依据。Abstract: The thermal characteristics of the spindle system has a significant influence on the machining accuracy of the heavy horizontal lathe. The major research object is the spindle system of heavy horizontal lathe. Based on the thermoset coupling analysis of the finite element method, the temperature field and thermal deformation field of the spindle system after reaching the thermal equilibrium state is simulated. And the shift of the spindle's center line that thermal deformation has led is analyzed. According to the structural characteristics of the hydrostatic bearing, the effect of the rotational speed of bearing, viscosity of hydraulic oil, the thickness of the hydrostatic oil film and the oil supply pressure to the thermal deformation of the spindle system are analyzed. The results show that all the factors have different-level influences on the deformation field of the spindle system, and therefore provide a reference for the optimal design and thermal error compensation of the spindle system.
-
表 1 常数C、n的值
位置 适用范围 流态 C n 竖平壁 104≤(GrPr)≤109 层流 0.59 1/4 109≤(GrPr)≤1012 紊流 0.12 1/3 水平板朝上 105≤(GrPr)≤2×107 层流 0.54 1/4 2×107≤(GrPr)≤3×1010 紊流 0.14 1/3 水平板热面朝下 3×105≤(GrPr)≤3×1010 层流 0.27 1/4 表 2 材料物理性能参数
材料名称 HT300 42CrMo 锡青铜 密度/kg·m3 7 350 7 850 8 800 比热/J×(kg·℃)-1 480 460 343 导热系数/W·(m·℃)-1 52.3 46.5 398 弹性模量/Gpa 135 210 110 热膨胀系数/℃-1×10-6 10.5 12.5 18.0 泊松比 0.25 0.3 0.34 表 3 不同转速下热流密度和对流系数
转速/r·min-1 热流密度/W·m-2 对流换热系数/W·(m2·k)-1 前轴承 后轴承 20 666.9 597.4 3.5 40 818.4 732.5 5.5 60 1 071.0 957.5 7.2 80 1 424.5 1 272.7 8.8 100 1 879.1 1 677.8 10.2 表 4 不同粘度下热流密度
液压油标号 热流密度/W·m -2 前轴承 后轴承 22 1 864.1 1 667.5 32 1 768.9 1 580.6 46 1 879.1 1 677.8 68 2 305.0 2 056.4 表 5 不同间隙下热流密度
油膜间隙/mm 热流密度/W·m-2 前轴承 后轴承 0.04 2 072.0 1 847.6 0.05 1 870.1 1 668.5 0.06 1 879.1 1 677.8 0.07 2 062.5 1 843.0 表 6 不同供油压力下热流密度
供油压力/Mpa 热流密度/W·m-2 前轴承 后轴承 4 1 536.7 1 370.9 5 1 690.8 1 509.0 6 1 879.1 1 677.8 7 2 101.7 1 877.3 -
[1] BRYAN J. International status of thermal error research[J]. CIRP Annals-Manufacturing Technology, 1990, 39(2):645-656. doi: 10.1016/S0007-8506(07)63001-7 [2] ZHANG J F, FENG P F, WU Z J, et al. Thermal structure design and analysis of a machine tool headstock[J]. Mechanika, 2013, 19(4):478-485. http://www.chemtech.ktu.lt/index.php/Mech/article/view/5044 [3] WU C, KUNG Y. Thermal analysis for the feed drive system of a CNC machine center[J]. International Journal of Machine Tools and Manufacture, 2003, 43(15):1521-1528. doi: 10.1016/j.ijmachtools.2003.08.008 [4] CREIGHTON E, HONEGGER A, TULSIAN A, et al. Analysis of thermal errors in a high-speed micro-milling spindle[J]. International Journal of Machine Tools and Manufacture, 2010, 50(4):386-393. doi: 10.1016/j.ijmachtools.2009.11.002 [5] 蔄靖宇, 赵海涛, 杨建国.车削中心主轴箱热误差仿真及特性分析[J].中国机械工程, 2009(18):2182-2186. http://www.cnki.com.cn/Article/CJFDTOTAL-ZGJX200918009.htm MAN Jing-yu, ZHAO Hai-tao, YANG Jian-guo. Simulation and analysis of thermal errors of a turning center[J]. China Mechanical Engineering, 2009(18):2182-2186. http://www.cnki.com.cn/Article/CJFDTOTAL-ZGJX200918009.htm [6] 李金华, 刘永贤, 韩家亮, 等.车削中心主轴系统热特性有限元分析的研究[J].组合机床与自动化加工技术, 2011(11):6-8. http://www.cnki.com.cn/Article/CJFDTOTAL-ZHJC201111003.htm LI Jin-hua, LIU Yong-xian, HAN Jia-liang, et al. Finite element analysis of thermal characteristics of the spindle system of a turning center[J]. Modular Machine Tool & Automatic Manufacturing Technique, 2011(11):6-8. http://www.cnki.com.cn/Article/CJFDTOTAL-ZHJC201111003.htm [7] 杨淇帆, 梁睿君, 叶文华.基于有限元的机床主轴箱部件的热特性分析[J].机械科学与技术, 2012(11):1860-1863. http://www.cnki.com.cn/Article/CJFDTOTAL-JXKX201211031.htm YANG Qi-fan, LIANG Rui-jun, YE Wen-hua. Thermal characteristics analysis of spindle box based on FEM method[J]. Mechanical Science and Technology for Aerospace Engineering, 2012(11):1860-1863. http://www.cnki.com.cn/Article/CJFDTOTAL-JXKX201211031.htm [8] BABU S R, RAJA V P, THYLA P R, et al. Prediction of transient thermo-mechanical behavior of the headstock assembly of a CNC lathe[J]. The International Journal of Advanced Manufacturing Technology, 2014, 74(1-4):17-24. doi: 10.1007/s00170-014-5916-2 [9] ZHANG Y M, GU R J, HAN J L. Thermal characteristic analysis on a high precision turning center's headstock[J]. Advanced Materials Research, 2013, 655-657:305-309. doi: 10.4028/www.scientific.net/AMR.655-657 [10] ZHANG Y M, CHEN Y, LI H. Temperature field analysis and thermal error testing for cnc machine tool's headstock[J]. Advanced Materials Research, 2014, 889-890:316-320. doi: 10.4028/www.scientific.net/AMR.889-890 [11] 陈兆年, 陈子辰.机床热态特性学基础[M].北京:机械工业出版社, 1998:25-33. CHEN Zhao-nian, CHEN Zi-chen. Basic of machine tool thermal state characteristics[M]. Beijing:China Machine Press, 1998:25-33. [12] 高士强.超高速磨削实验机床液体动静压混合轴承的热态特性研究[D].沈阳:东北大学, 2008. http://cdmd.cnki.com.cn/Article/CDMD-10145-2010258963.htm GAO Shi-qiang. Study on thermal properties of hybrid journal bearing for super-high speed grinding platform[D]. Shenyang:Northeastern University, 2008. http://cdmd.cnki.com.cn/Article/CDMD-10145-2010258963.htm [13] 刘殿崇.重型数控落地铣镗床滑枕组件热误差补偿技术研究[D].哈尔滨:哈尔滨工业大学, 2011. http://cdmd.cnki.com.cn/Article/CDMD-10213-1012001530.htm LIU Dian-chong. Research on thermal error compensation for ram components of NC heavy type milling-boring machine tool with adjustable stand[D]. Harbin:Harbin Institute of Technology, 2011. http://cdmd.cnki.com.cn/Article/CDMD-10213-1012001530.htm [14] 赵镇南.传热学[M].北京:高等教育出版社, 2008:172-286. ZHAO Zhen-nan. Heat transfer theory[M]. Beijing:Higher Education Press, 2008:172-286. [15] 曾正明.机械材料手册金属材料[M].北京:机械工业出版社, 2010:32-35. ZENG Zheng-ming. Mechanical materials manual-metallic material[M]. Beijing:China Machine Press, 2010:32-35.