数控车床综合热误差建模及工程应用-数控机床市场网

数控车床综合热误差建模及工程应用

2018-1-26 来源：机械制造系统国家实验室河北工程大学作者：孙志超侯瑞生陶涛杨军梅雪松王新访问量：

摘　要：针对车床实际工程应用中主轴与进给轴综合误差对工件加工精度产生影响问题，建立包含工件膨胀效应的主轴与进给轴综合热误差模型，并进行实际切削验证．以精密车床为研究对象，综合分析车床主轴、进给轴和工件在实际加工中的相互影响关系，并建立三者之间的综合热误差多元线性回归模型（ＭＬＲＡ）．实验结果表明：含有工件膨胀效应系数的综合热误差模型符合实际工况，有效提高了车床的加工精度．主轴热误差模型的预测精度达８５％以上，进给轴预测精度达７０％以上，实际加工中工件误差由１５ μｍ降低到５ μｍ左右．综合热误差模型显著提高了高精密数控车床的加工精度．

关键词：数控车床；主轴；进给轴；工件膨胀效应；热误差建模；热误差补偿

车床在轴类及盘类零件加工中占有显要位置，我国数控车床主要存在精度低、精度保持性差等问题，影响机床精度的关键因素之一热误差占据机床总体误差的４０％～７０％［１］，而对于高精密数控车床来说所占比重更大．近年来，国内外针对机床热特性的研究不胜枚举，也取得了一些良好的效果．杨军等［２－５］利用模糊聚类选择温度变量，建立了机床主轴热误差的ＢＰ神经网络模型、多元线性回归模型、最小二乘支持向量机模型及时间序列模型，并在多种工况下验证模型的准确性及鲁棒性；还有学者通过实验反求热流密度和热辐射等边界条件，提高机床热变形的仿真精度［６－８］；Ｂｏｓｓｍａｎｎｓ等［９－１０］利用有限差分模型分析并预测了电主轴热源的分布机理；Ａｇｕａｄｏ［１１］提出机床空间误差的测量方法；Ｈｅｉｓｅｌ等［１２］研究丝杠温度场分布，并建立了进给轴热误差模型；Ｇｕｏ和Ｓｈｅｎ等［１３－１４］利用不同的算法对ＢＰ神经网络模型进行了优化，提高了模型精度；径向基函数ＲＢＦ（ＲａｄｉｕｓＢａｓｉｓＦｕｎｃｔｉｏｎ）的神经网络模型被广泛用于机床热误差建模中［１５－１６］；Ｃｈｅｎ等［１７－１９］建立了主轴系统的多元线性回归模型．

现有文献都是针对主轴或进给轴单个系统分别建立模型，而且都是在理论上验证模型的准确性和鲁棒性，没有进行过实际的加工验证，实际加工中主轴与进给轴相互依赖共同影响工件的加工精度．本文针对ＨＴＣ５５０／５００车床建立主轴与进给轴的综合热误差模型，并进行补偿应用，并用实际加工来验证模型的准确性．

１、　Ｓｉｅｍｅｎｓ系统热误差补偿方式及热特性实验

１．１　Ｓｉｅｍｅｎｓ系统热误差补偿方式

Ｓｉｅｍｅｎｓ开放了热误差补偿接口，一定温度下所开放的补偿模型为线性模型，主轴的热误差只与温度相关，与坐标位置无关；进给轴热误差不仅与温度相关，且与坐标位置相关．热误差模型原理图如图１所示．

图１　热误差补偿原理：温度θ 下热误差的近似拟合线

Ｓｉｅｍｅｎｓ内部模型：

１．２　热特性实验

１．２．１　实验原理及方法

以精密数控车床为研究对象，测试设备包括：ＲＥＮＩＳＨＡＷ激光干涉仪测量进给轴误差；自主设计的温度与位移同步采集系统测得温度及变形数据；传感器选用高精密温度传感器ＰＴ１００和高精密电涡流传感器．采用五点法测量主轴空间变形［２０］，原理如图２所示．

图２　主轴热误差测量原理示意图

　　Ｓ１、Ｓ３为测量主轴Ｘ方向的热误差，Ｓ２、Ｓ４为测量主轴Ｙ向的热误差，Ｓ５测量主轴Ｚ向的热误差；文献４中详细介绍了利用激光干涉仪测量进给轴热误差的测量方法及注意事项，冷态下第一次测量进给轴误差为机床进给系统几何误差，进给系统连续往复运行２０ｍｉｎ后测量误差值，此误差值减去几何误差作为此时进给系统热误差［４］．

１．２．２　热特性分析

进给轴电机、轴承、丝杠螺母副等摩擦生热导致丝杠温度升高，丝杠向自由端方向发生热膨胀引起丝杠导程变化．丝杠导程变化导致半闭环控制系统产生误差，进给轴的热误差变化如图３所示．冷态下第１次测量值为机床的几何误差，故冷态下机床的热误差为０ μｍ．由图３可以看出，进给轴热误差不仅与温度相关，而且与坐标位置相关，随坐标值的增大而增大；负向热误差变化相对较小，正向热误差变化相对较大，由此判断正向为进给轴自由端，即丝杠热膨胀的方向．

图３　进给轴热误差曲线

主轴内部结构如图４（ａ）所示，轴承摩擦生热，引起主轴及外壳温度升高，导致主轴发生热伸长Δｚ和热升高Δｈ．主轴热变形曲线如图４（ｂ）所示，主轴Ｚ向热误差Δｚ最高达４４ μｍ，Ｘ向热误差Δｘ最高达１４ μｍ，随着温度的升高，主轴的热变形随之增大，停机后随着温度的降低主轴的热变形随之减小．

图４　主轴结构示意图及主轴Ｘ／Ｚ向热误差

机床热特性实验主轴及进给轴温度变化见图５．

图５　主轴与进给轴温度场变化

主轴传感器ＰＴ１００配置前端３个、中部２个、后端３个，前端最高温度达３５．９ ℃、后端３３．８ ℃、中部３２．８ ℃，其中前部最高温差１３．２ ℃、后端１２．１ ℃、中部１０．８ ℃．主轴内部结构如图４（ａ）所示，循环空气冷却中空式结构，前端３个轴承后端１个，导致前端发热量最大，后端次之，中部最小，温度场变化与结构相符合．由图５（ａ）和图４（ｂ）对比可看出，曲线变化规律一致，变形与温度之间具有一定的线性关系．

　　进给轴传感器ＰＴ１００配置主要在电机、轴承及丝杠螺母座上，通过螺母座温度间接反映丝杠温度变化．其中，电机温度变化最大，床身温度变化最小；前轴承温度大于后轴承；Ｚ轴螺母座温度大于Ｘ轴螺母座．

２、热误差建模及补偿实现

２．１　综合热误差建模

在车床加工过程中，主轴与进给轴热误差相互耦合共同影响工件的精度，因此需要建立主轴与进给轴的综合热误差模型．

２．１．１　主轴热误差模型

主轴热特性实验中，电涡流传感器的安装位置影响测量结果，以Ｘ向热误差测量为例说明．图６为主轴热特性实验传感器安装主轴轴向视图

图６　传感器安装主轴轴向示意图

　　传感器的安装支架安装在刀塔上．实验过程非恒温，环境温度的升高导致丝杠温度升高，并伴随着热伸长，Ｘ轴丝杠热伸长导致刀塔位置发生变化，从而引起传感器相对于测量芯棒的位置变化，导致传感器测量主轴Ｘ向热误差就包含了Ｘ轴丝杠热变形误差．因此，主轴Ｘ向热误差建模中要消除Ｘ轴丝杠的热变形误差，处理方法：

２．１．２　进给轴热误差模型

由于轴承及丝杠螺母副摩擦发热，丝杠温度升高导致进给系统产生热误差．然而，在实际加工中工件也会发生热胀冷缩现象，同样会影响工件的加工精度．图７为工件与丝杠变形示意图．相同温度下所有工件的长度均为Ｌ，在相同温升条件下工件１、２、３的膨胀量分别为Δ１、Δ２、Δ３，丝杠在相同长度Ｌ上的膨胀量为ΔＳ．当工件的膨胀系数＜丝杠的膨胀系数，即当ΔＳ＞ Δ１时，丝杠的膨胀量大于工件的膨胀量，此时的补偿量为丝杠与工件膨胀量的差值，方向为丝杠膨胀反方向；当丝杠的膨胀系数等于工件的膨胀系数，即当ΔＳ＝ Δ２时，丝杠的膨胀量与工件膨胀量相同，此时丝杠的膨胀量刚好补偿了工件的膨胀量，不需要对丝杠的膨胀量进行补偿；当丝杠的膨胀系数小于工件的膨胀系数，即ΔＳ＜ Δ３时，丝杠的膨胀量小于工件的膨胀量，此时的补偿量亦为丝杠与工件膨胀量的差值，方向为丝杠膨胀方向．因此，进给轴的热误差补偿要考虑工件的膨胀效应，补偿方法：

图７　工件与丝杠变形示意图

２．１．３　综合热误差模型

选取主轴及床身温度为温度变量，结合ＭＬＲＡ方法得到如下主轴的热误差模型：

式中： θ１、θ２、θ３、θ４分别为床身、主轴、Ｘ轴螺母和Ｚ轴螺母温度；ｔａｎ β （θ）ｓｐｅｅｄ－Ｘ、ｔａｎ β （θ）ｓｐｅｅｄ－Ｙ分别为

Ｘ、Ｙ轴丝杠膨胀系数；ＰＸ、ＰＹ为进给轴坐标；Ｐ０Ｘ、Ｐ０Ｙ为进给轴参考点坐标值．模型中将２０ ℃作为参考温度是因为ＧＢ中将２０ ℃ 作为检测时标准环境温度，机床定位精度检测标准环境温度亦为２０ ℃．

　　图８、９为主轴及Ｘ／Ｚ进给轴热误差模型预测值与实验值的对比图．建立模型预测精度评价标准均方根误差值ＲＭＳＥ及预测精度η［５］．其中Ｒ为均方根误差值，ｙｉ为实验测量值，ｙ～ｉ为模型预测值．主轴Ｘ／Ｚ方向热误差模型的Ｒ和η 分别为２．５、５．２ μｍ和８９．４％、８８．７％；Ｘ／Ｚ进给轴热误差模型的Ｒ和η 分别为２．４、５．１ μｍ和８４．５％、８２．７％．说明热误差模型有一定准确性，应用效果还需进一步验证．

图８　主轴热误差模型预测值与测量值的比较

图９　进给轴热误差模型预测值与测量值的比较

机床主轴系统与进给轴系统为相对独立的个体，实际加工中二者缺一不可，轴与进给轴的热误差相互关联共同影响工件的加工精度．得到主轴与进给轴热误差的相互关系对于模型的建立尤为重要．主轴系统由于轴承及加工摩擦生热造成主轴系统温度升高，随之产生热变形Δｌｚ、Δｈ，丝杠受热发生膨胀导致进给系统产生热误差Δｘ、Δｚ，如图１０所示．由于Ｘ轴是倾斜式安装，主轴热变形Δｈ在机床Ｘ方向产生分量Δｌｘ，方向与Ｘ进给轴相同；主轴热伸长Δｌｚ方向与Ｚ进给轴方向相同．因此机床Ｘ／Ｚ方向的热误差模型应该是主轴与进给轴的综合热误差模型．结合式（１）～（４）及文章２．２．１、２．２．２节分析得到机床在Ｘ／Ｚ方向的综合热误差模型．

图１０　主轴与进给轴热误差耦合示意图

Ｘ方向热误差模型：

２．２　补偿实现及实际加工分析

２．２．１　补偿实现

Ｓｉｅｍｅｎｓ８２８Ｄ热误差补偿总体方案如图１１所示，由ＰＬＣ直接获取机床热源温度值，在ＰＬＣ内计算补偿参数，最后ＰＬＣ通过数据接口ＤＢ１２００将补偿参数写入ＮＣ系统内，系统根据补偿参数及插补指令计算正确的电机指令从而达到补偿效果，提高机床的加工精度．

　　冷态下测量机床的定位精度，随后同时运转主轴及进给轴系统，模拟实际加工主轴与进给轴热误差耦合现象，验证综合热误差模型的准确性，直到机床达到热平衡状态．测量机床热误差补偿前后的定位精度，结果如图１２所示．热补偿前Ｘ／Ｚ轴定位精度分别为１９．８ μｍ、２７．２ μｍ；热补偿后Ｘ／Ｚ轴定位精度分别为６．９、９．１ μｍ，热补偿后Ｘ／Ｚ轴定位精度分别提高了６５．２％、６８．４％，表明热误差综合模型有一定的补偿效果．

图１１　Ｓｉｅｍｅｎｓ８２８Ｄ热误差补偿总体方案

图１２　Ｘ／Ｚ轴热平衡下热误差补偿前后对比

２．２．２　实际加工分析

加工工件如图１３所示，严格按照工程实际在有无热误差补偿状态下按图纸要求进行加工，两种状态下各加工一天，对工件按照加工顺序做编号．将加工好的工件置于２０ ℃的恒温环境中８ｈ以上，按编号使用三坐标测量仪测量工件Ｒ１、Ｒ２的直径，记录于表格，比较有无热误差补偿状态下的工件误差。

实际结果如图１４所示．由图１４可知，在有无热误差补偿状况下工件误差首先负向变大而后向正向变化，这是由于Ｘ轴丝杠的安装在Ｘ负向有预拉伸，丝杠温升初始时首先要消耗预拉伸量，因此导致工件误差负向变化．图１４（ａ）所示预拉伸消耗之后工件误差正向有明显变化，跨度１５ μｍ，这便是热误差造成的影响；图１４（ｂ）所示预拉伸消耗之后工件误差有了明显改善，跨度５ μｍ左右，由此证明热误差补偿的准确性．