数控机床热误差补偿方法
技术领域
本发明涉及一种数控机床热误差补偿方法。
背景技术
随着现代机械制造技术的飞速发展,高柔性、高效率、高质量、高精度、高集成和高智能加工技术已成为现代机械制造中最重要的组成部分和发展方向,并成为提升一个国家工业发展水平的重要手段,得到了世界各国的广泛重视,进而对数控机床的制造、加工精度提出了越来越高的要求。然而随着数控机床精度的进一步提高,热误差逐渐成为了影响数控机床加工精度的主要因素。研究表明,在导致机床变形的众多误差源中,机床的几何误差和机床各零部件由于受热发生变形而引起的热误差是影响机床加工精度最重要的两个因素,分别为机床全部误差的20%和45%,且越是精密的数控机床,由机床内部与外部热源共同作用导致的热误差在机床总误差中所占的比例越大。随着现代制造技术的发展,机床的几何误差已经得到了有效的控制和补偿,其在机床总误差中所占百分比已经大大减小。因此在现代精密数控机床的生产加工过程中,热误差已经成为了阻碍数控机床加工精度进一步提高的最重要因素。如何减小机床的热误差、改善机床的加工精度,已经受到世界各国学者的广泛关注,总体来说,减小机床热误差的技术主要包括误差防止技术和误差补偿技术。前者主要依靠优化机床结构,减少热源生热;后者主要依靠误差检测以及误差反馈补偿。在机床热误差补偿中,建立高精度的热误差补偿模型是难点,而选择温度测点的数量和位置又是热误差建模的基础。所以合理又可行的温度测点选取方案和误差补偿模型建立显得尤为重要。
通过国内外研究学者对机床热误差补偿技术的大量研究,机床热误差补偿技术已经取得了很大的进展,但目前还未见高水平的热误差补偿技术在数控机床实际工业生产制造中大批量应用。传统的方法繁琐且可执行性差,较难获得合适的温度补偿点。
发明内容
本发明所要解决的技术问题是提供一种数控机床热误差补偿方法,本方法在有限元仿真分析的基础上基于粒子群优化算法,得到热误差补偿的关键温度点,以关键温度点处的温度测量值作为数控机床热误差补偿依据,简化了热误差温度测量点的选取,有效减小机床的热误差,提高机床的加工精度。
为解决上述技术问题,本发明数控机床热误差补偿方法包括如下步骤:
步骤一、建立数控机床主轴结构的有限元仿真分析几何模型,几何模型建立过程忽略对主轴结构影响较小的螺孔,油管,注油孔,细微倒角、圆角特征,采用简单圆环代替主轴结构的轴承、轴承套特征;
步骤二、对建立的几何模型进行网格划分和扫掠划分,主轴结构的主轴、主轴箱复杂实体采用网格划分,主轴结构的轴承、轴承套简单实体采用扫掠划分,各部件之间的接触关系选用粘接,各部件的材料根据机床主轴结构实际材料选定;
步骤三、综合测量和计算主轴结构工作时的环境温度,冷却液流量,进出油温度数据,以此对建立的几何模型施加边界条件;
步骤四、对主轴结构布置热关键点并进行温度测量,一个温度传感器测量环境温度,三个温度传感器测量主轴后端轴承处温度,十三个温度传感器测量主轴后端轴承与主轴前端轴承之间的主轴温度并间隔布置,三个温度传感器测量主轴前端轴承处温度,各热关键点依次编号;
步骤五、对主轴结构的热关键点进行有限元仿真分析,有限元仿真分析结果作为粒子群优化算法的参数,利用粒子群优化算法筛选出四个热关键点参数,给定目标函数:
(1)
式(1)中为各热关键点测得的温度值,为有限元仿真分析中各热关键点测得温度值的最小二乘支持向量机的预测值,为样本数目,分别为筛选的四个热关键点采集的温度值,分别为有限元仿真分析中最小二乘支持向量机的正则化参数和校核参数;
步骤六、将式(1)的目标函数最小作为粒子群优化算法目标,采用粒子群优化算法对参数进行优化,选择各热关键点的一部分温度值和热变形数据作为粒子群优化算法的训练集,选择各热关键点再次测得的温度值和热变形数据作为粒子群优化算法的检测集,给定最大迭代次数、学习因子、最大权重因子、最小权重因子,群体规模为每个粒子的位置六维向量,六维向量分别代表参数;
步骤七、根据各热关键点的当前温度参数,通过式(1)计算每个粒子的适应度值,将与粒子群优化算法的当前个体最优解的适应度值比较,如,则用该粒子的适应度值取代当前个体最优解,即;
步骤八、将每个粒子的最好适应度值与全局适应度值进行比较,如,则取,同时记忆个体与群体所对应最佳适应度值的位置和;
步骤九、根据式(2)和式(3)更新粒子的位置和速度,
(2)
(3)
其中:i为粒子数,取1、2、3、4…,为第i个粒子的位置,为第i个粒子的飞行速度,为第i个粒子在第d维的飞行速度,为第i个粒子在第d维的个体最优解,为全部粒子在第d维的最优解,k为迭代次数,c1和c2为学习因子,取值在[0,2]区间,rand()为0到1的随机数,为权重因子;
步骤十、重复步骤六至步骤九,直至每个粒子的适应度值随着迭代次数的增加而收敛至设定值,其浮动范围小于0.1%,选择收敛速度最快的一组粒子并输出结果,得到最优的四个热关键点的序号、温度值和各个方向的热形变量;
步骤十一、选取最优的四个热关键点的参数作为热误差补偿的参数,数控机床控制系统建立误差补偿模型,数控机床控制系统根据热误差补偿参数通过误差补偿模型对主轴结构的热误差作出补偿。
进一步,上述各热关键点的温度传感器是磁力吸附的Pt100温度传感器。
由于本发明数控机床热误差补偿方法采用了上述技术方案,即本方法对数控机床主轴结构建立有限元仿真分析几何模型并作简化处理,根据主轴结构的工作状态及工作环境,加载合理的边界条件,将简化后的几何模型进行热力学和静力学的有限元仿真分析,在分析结果中取出一定数量点的温度以及热变形量,将这些数据用粒子群优化算法进行筛选,最终选取最优的四个点的温度值作为数控机床热误差补偿参数并建立误差补偿模型,实现数控机床热误差补偿。本方法在有限元仿真分析的基础上基于粒子群优化算法,得到热误差补偿的关键温度点,以关键温度点处的温度测量值作为数控机床热误差补偿依据,简化了热误差温度测量点的选取,有效减小机床的热误差,提高机床的加工精度。
附图说明
下面结合附图和实施方式对本发明作进一步的详细说明:
图1为本发明数控机床热误差补偿方法的流程图;
图2为本方法中粒子群优化算法的流程图;
图3为本方法中主轴结构热关键点布置示意图。
具体实施方式
如图1、图2和图3所示,本发明数控机床热误差补偿方法包括如下步骤:
步骤一、建立数控机床主轴结构的有限元仿真分析几何模型,几何模型建立过程忽略对主轴结构影响较小的螺孔,油管,注油孔,细微倒角、圆角特征,采用简单圆环代替主轴结构的轴承、轴承套特征;
步骤二、对建立的几何模型进行网格划分和扫掠划分,主轴结构的主轴、主轴箱复杂实体采用网格划分,主轴结构的轴承、轴承套简单实体采用扫掠划分,各部件之间的接触关系选用粘接,各部件的材料根据机床主轴结构实际材料选定;其中粘接(Bonded)是分析有限元中的一种接触方式,即将两个接触部件绑定在一起;
步骤三、综合测量和计算主轴结构工作时的环境温度,冷却液流量,进出油温度数据,以此对建立的几何模型施加边界条件;
步骤四、对主轴结构21布置热关键点并进行温度测量,一个温度传感器1测量环境温度,三个温度传感器2、3、4测量主轴后端轴承22处温度,十三个温度传感器5、6、7、8、9、10、11、12、13、14、15、16、17测量主轴后端轴承22与主轴前端轴承23之间的主轴温度并间隔布置,三个温度传感器18、19、20测量主轴前端轴承23处温度,各热关键点依次编号;
步骤五、对主轴结构的热关键点进行有限元仿真分析,有限元仿真分析结果作为粒子群优化算法的参数,利用粒子群优化算法筛选出四个热关键点参数,给定目标函数:
(1)
式(1)中为各热关键点测得的温度值,为有限元仿真分析中各热关键点测得温度值的最小二乘支持向量机的预测值,为样本数目,分别为筛选的四个热关键点采集的温度值,分别为有限元仿真分析中最小二乘支持向量机的正则化参数和校核参数;
步骤六、将式(1)的目标函数最小作为粒子群优化算法目标,采用粒子群优化算法对参数进行优化,选择各热关键点的一部分温度值和热变形数据作为粒子群优化算法的训练集,该一部分温度值和热变形数据为在有限元仿真分析软件中,分别随机选取主轴前端轴承上的一个测试点、主轴后端上的一个测试点、以及主轴轴身上的两个测试点得到的测试数据,热变形数据可以从有限元仿真分析结果中取出热形变量;选择各热关键点再次测得的温度值和热变形数据作为粒子群优化算法的检测集,给定最大迭代次数、学习因子、最大权重因子、最小权重因子,群体规模为每个粒子的位置六维向量,六维向量分别代表参数;
步骤七、根据各热关键点的当前温度参数,通过式(1)计算每个粒子的适应度值,将与粒子群优化算法的当前个体最优解的适应度值比较,如,则用该粒子的适应度值取代当前个体最优解,即;
步骤八、将每个粒子的最好适应度值与全局适应度值进行比较,如,则取,同时记忆个体与群体所对应最佳适应度值的位置和;
步骤九、根据式(2)和式(3)更新粒子的位置和速度,
(2)
(3)
其中:i为粒子数,取1、2、3、4…,为第i个粒子的位置,为第i个粒子的飞行速度,为第i个粒子在第d维的飞行速度,为第i个粒子在第d维的个体最优解,为全部粒子在第d维的最优解,k为迭代次数,c1和c2为学习因子,取值在[0,2]区间,rand()为0到1的随机数,为权重因子;权重因子可以通过下式计算得到:
(4)
式(4)中,为实验要求的最大和最小权重因子,为迭代次数,为最大迭代次数;
步骤十、重复步骤六至步骤九,直至每个粒子的适应度值随着迭代次数的增加而收敛至设定值,其浮动范围小于0.1%,选择收敛速度最快的一组粒子并输出结果,得到最优的四个热关键点的序号、温度值和各个方向的热形变量;设定值为全局的最优适应度值;
步骤十一、选取最优的四个热关键点的参数作为热误差补偿的参数,数控机床控制系统建立误差补偿模型,数控机床控制系统根据热误差补偿参数通过误差补偿模型对主轴结构的热误差作出补偿。通常误差补偿模型可采用最小二乘法建立,其表达式为,其中:为热误差补偿值,为常数,t为热关键点与环境的温度差值,三个常数根据四个热关键点的参数,即各个方向的热变形量以及温度差值计算得到。
进一步,上述各热关键点的温度传感器是磁力吸附的Pt100温度传感器。
本方法中各热关键点选取主要是根据主轴结构工作时的热源以及温度场情况,在主轴上的各热关键点序号分别为2#至20#。主要包括主轴的前后轴承以及主轴上的关键点;有限元仿真分析主要是对于在不同转速下的主轴结构进行热力学和静力学的有限元仿真,以确定在不同转速下主轴结构的温度场和热变形量。本方法的粒子群优化算法中,将各热关键点的温度传感器采集到的数据一部分用于训练数控机床的最小二乘支持向量机热误差补偿模型,可设置粒子群优化算法的粒子群规模为40、粒子向量维数为6、学习因子、最大权重因子、最小权重因子、最大迭代次数;粒子群优化算法仅需约25步就可找到全局最优解,收敛速度非常快,最终适应度值收敛到0.5873,经粒子群优化算法搜索到的热关键点和最小二乘支持向量机的参数分别为3、9、15、19、95.5774、0.55,即选择位于主轴后端轴承处的3#热关键点、位于主轴上的9#热关键点和15#热关键点、位于主轴前端轴承处的19#热关键点的温度值用于建立数控机床热误差补偿模型,作为热误差补偿的依据。
下表为数控机床的主轴结构采用本方法作热误差补偿前后的参数对照表:
主轴结构坐标方向 | 补偿前最大误差(μm) | 补偿后最大误差(μm) | 补偿前平均误差(μm) | 补偿后平均误差(μm) |
X向 | 8.23 | 0.96 | 4.70 | 0.22 |
Y向 | 57.96 | 2.43 | 25.22 | 0.56 |
Z向 | 29.68 | 1.50 | 21.73 | 0.28 |
从上表可知,本方法可以有效减少数控机床的热误差,从而改善数控机床的加工精度,提高了数控机床制造及加工精度的水平。