一种数控机床滚珠丝杠进给系统热误差预测方法技术方案

本发明专利技术属于热误差预测方法技术领域，提出了一种数控机床滚珠丝杠进给系统热误差预测方法。利用自适应实时模型(ARTM)预测滚珠丝杠进给系统的热误差。本方法首先建立实验测量工作台表面温度随时间的变化，然后建立一种自适应实时模型用来预测滚珠丝杠的瞬时温度分布和热误差分布。利用有限元结合蒙特卡洛方法(MC)确定滚珠丝杠进给系统两个轴承、移动丝母、滑轨的发热率。然后依据有限元计算数据，建立进给速度与时间的指数方程，用来描述测量表面点和运动对中心之间的温度差随时间的变化，最后建立一种数值预测算法用来预测滚珠丝杠进给系统的热误差。通过监测两个轴承座和运动螺母侧面的表面温度，使用数值预测算法快速高精度地预测相应的热误差。

A thermal error prediction method for ball screw feed system of CNC machine tools

The invention belongs to the technical field of thermal error prediction method, and puts forward a thermal error prediction method for ball screw feed system of numerical control machine. The adaptive real-time model (ARTM) is used to predict the thermal error of the ball screw feed system. This method first establishes the change of the surface temperature of the experimental worktable with time, and then establishes an adaptive real-time model to predict the instantaneous temperature distribution and the thermal error distribution of the ball screw. Using finite element method and Monte Carlo method (MC), the heating rate of two bearings, moving wire and slide rail of ball screw feed system is determined. Then according to the finite element calculation data, the exponential equation of feed speed and time is established to describe the change of temperature difference between the measuring surface point and the motion between the center and the time. Finally, a numerical prediction algorithm is established to predict the heat error of the ball screw feed system. By monitoring the surface temperature of the two bearing seats and the nut of the movement nut, the numerical prediction algorithm is used to predict the corresponding thermal error quickly and accurately.

【技术实现步骤摘要】
一种数控机床滚珠丝杠进给系统热误差预测方法
本专利技术属于热误差预测方法
，涉及数控机床滚珠丝杠进给系统热误差建模方法。
技术介绍
滚珠丝杠作为数控机床进给系统的关键部件，工作时轴承、丝杠螺母等部件摩擦生热引起丝杠温度的上升，其温度的变化导致机械结构的变形,从而产生热误差降低了进给系统的定位精度。由于热误差的存在，单件的加工精度易不合格；批量加工产品的一致性较差。因此，机床热误差问题一直是近期国内外研究的重点，已经成为了提高数控机床定位精度的重要研究方向。近年来，国内外研究学者建立了几种滚珠丝杠热误差模型的有效方法，但是都存在其局限性。相关研究主要集中两个方面：第一种，对滚珠丝杠进给系统的结构进行简化，应用传热学理论为基础近似计算系统的热边界条件，建立滚珠丝杠进给系统热误差有限元模型。第二种，实测大量数据进行处理，建立滚珠丝杠进给系统热误差经验模型。通过对目前研究的分析。
技术实现思路
为了克服现有技术中存在的问题，本专利技术提出了一种数控机床滚珠丝杠进给系统热误差预测方法。该方法利用自适应实时模型(ARTM)预测滚珠丝杠进给系统的热误差。本方法首先建立实验测量工作台表面温度随时间的变化，然后建立一种自适应实时模型用来预测滚珠丝杠的瞬时温度分布和热误差分布。利用有限元结合蒙特卡洛方法(MC)确定滚珠丝杠进给系统两个轴承、移动丝母、滑轨的发热率。然后依据有限元计算数据，建立进给速度与时间的指数方程，用来描述测量表面点和运动对中心之间的温度差随时间的变化，最后建立一种数值预测算法用来预测滚珠丝杠进给系统的热误差。本专利技术的具体技术方案为：步骤1，滚珠丝杠进给系统温度测量确定进给系统表面温度检测点，其中，温度检测点包括两个丝杠支撑轴承的表面点、丝母法兰的表面点和丝杠轴一侧导轨表面均匀分布点；T型贴片热电偶通过磁性底座分别固定在两个丝杠支撑轴承的表面点、丝母法兰的表面点上，按周期进行连续采样；通过红外线测量仪对测量点进行定时检测；步骤2，建立其有限元模型；步骤2-1，根据目标数控机床进给系统的几何尺寸，建立其有限元模型，模型包括伺服电机、两个支撑轴承、丝杠、螺母和两个滑轨，各部分均选择固体单元，而丝杠与两个支撑轴承、丝母之间的连接选用接触单元；步骤2-2，确定系统的热源；系统包括固定热源和移动热源；固定热源为伺服电机、轴承一和轴承二；依据系统结构，将伺服电机发热率合并到轴承一，记为QB1，轴承二的发热率记为QB2；移动热源为两个滑轨和一个丝母，在此只考虑各热源传入模型部分的发热率，两个滑轨发热率为Qg1＝Qg2＝Qg，丝母的发热率Qn；步骤2-3，对流换热系数的确定；在丝杠旋转过程，丝杠表面与空气产生对流传热，则对流换热系数为：其中，Nu为努塞尔数；λfluid为空气的导热系数；d为丝杠的直径。步骤3，用蒙特卡洛方法实现发热率的有限元计算；利用x1，x2，x3和x4分别表示发热率QB1，QB2，Qn和Qg，则蒙特卡洛模拟计算为优化如下的目标函数：其中，Tij是在第j个采样时间步长的给定点i的温度，i＝1,2,…,15，j＝1,2,…,N，上标EM和MC分别代表实验测量值和有限元结合蒙特卡洛的模拟值；蒙特卡洛方法具体为：首先确定采样次数为k，并给定每个参数xl的取值[xLl,xUl]区间；在每一次采样过程中，确定各变量的数值：xl＝xLl+random(0,δl),δl＝xUl-xLl,l＝1,2,3,4.(5)每次采样后进行有限元的计算，提取各检测点的温度，依据式(4)求出目标函数F，并将k次采样计算的最小目标函数记为F*＝min(F0,F1,…,Fk-1)，而对应参数记为和完成k次蒙特卡洛模拟后，将F*与给定的计算精度ε进行比较，如果达到给定计算精度，则输出发热率和结束模拟计算；否则，判断取值区间大小：min(δ1,δ2,δ3,δ4)<δ0(6)其中，δ0为给定常数。如果式(6)满足，则计算失败，计算过程结束；否则，修改各变量的随机抽样区间如式(7)，并进行重新模拟计算。步骤4，运动副中心与表面温度关系的确定描述测量表面点和运动副对中心温度差：ΔT(t,v)＝av(1-e-t/bv)(8)其中，v是工作台的进给速度，t是工作时间，a和b是通过曲线拟合得到的常数；在机床工作过程中，检测支撑轴承或丝母表面的温度，计算得到其中心温度：Tc＝Tm+ΔT(t,v)(9)其中，Tc运动对的中心温度，Tm是其表面上给定点的测量温度；步骤5，滚珠丝杠实时热误差预测；步骤5-1建立滚珠丝杠热传导模型将丝杠简化为一维杆，其热传导方程：其中，T(x,t)是时间t和位置x的函数，表示t时刻丝杠上距离热源为x的位置点的温度变化；κ是导热系数；ρ是滚珠丝杠的材料密度；c是比热容；h是对流系数；Tf(t)是工作环境空气的温度；对滚珠丝杠T(x,t)的函数划分空间步长s和时间步长τ，xk＝k·s，k＝1,2,…,M，tj＝j·τ，j＝1,2,…,N；用表示一阶导数和二阶导数的有限差分方程如下：把差分方程式(11)和式(12)代入热传导方程式(10)，得到：其中A＝ρc/κ；B＝4h/(κd)，当(As2-Bs2+2τ)/(As2)≥0且(τ/s2)≤(A-B)/2时，等式成立；依据式(9)，两个轴承中心的温度如下式表示：Tcbl＝Tmbl+ΔTbl(t,v)l＝1,2.(14)移动螺母中心的温度表示为：Tcn＝Tmn+ΔTn(t,v)(15)步骤5-2，滚珠丝杠热误差计算；在加工过程中，由于摩擦热引起的内部温度上升导致滚珠丝杠轴的伸长；热伸长值如下式进行计算:ΔL＝L·α·ΔT(16)其中，L是滚珠丝杠的长度，因此滚珠丝杠在X轴的热伸长量：其中ΔL(x)表示t时刻的热延伸量，ΔT(x,t)表示滚珠丝杠的温度变化，α是线性热膨胀系数。本专利技术的有益效果：本文提出利用蒙特卡洛模拟所得发热率有限元计算，可用于预测滚珠丝杠进给驱动系统的发热率，温度分布和热误差。提出的自适应实时模型(ARTM)可用于预测滚珠丝杠进给驱动系统的热误差。通过监测两个轴承座和运动螺母侧面的表面温度，使用数值预测算法快速高精度地预测相应的热误差。附图说明图1是滚珠丝杠进给系统温度测量点位置图；图2是结合蒙特卡洛法的有限元模拟流程图；图3是试验检测温度与利用蒙特卡洛法进行有限元计算所得温度结果变化图；EM代表实验测量数据，FEM代表有限元计算数据；图4是轴承1在三种进给速度下温度差计算有限元结果和拟合曲线的对比。图5是滚珠丝杠进给系统热源简化图；图6是滚珠丝杠温度测量点位置图；图7是进给速度为10m/min时，有限元方法；自适应实时模型预测计算；试验结果三者之间的温升对比图；图8是进给速度为10m/min时，有限元方法，自适应实时模型预测计算，试验结果三者之间热误差对比图。具体实施方式以下结合附图和具体实施例对技术方案进行详细阐述。一种数控机床滚珠丝杠进给系统热误差预测方法，包括如下步骤：步骤1，滚珠丝杠进给系统温度测量确定进给系统表面温度检测点，如图1所示。其中，T1和T8为两个丝杠支撑轴承表面点，T9为丝母法兰表面点；其余在丝杠轴一侧导轨表面均匀分布点，如T2-T7和T10-T15。用磁性底座分别将三个T型贴片热电偶固定在T1，T8和T9点，有一台笔记本文档来自技高网...

【技术保护点】
1.一种数控机床滚珠丝杠进给系统热误差预测方法，其特征在于，包括如下步骤：步骤1，滚珠丝杠进给系统温度测量确定进给系统表面温度检测点，其中，温度检测点包括两个丝杠支撑轴承的表面点、丝母法兰的表面点和丝杠轴一侧导轨表面均匀分布点；T型贴片热电偶通过磁性底座分别固定在两个丝杠支撑轴承的表面点、丝母法兰的表面点上，按周期进行连续采样；通过红外线测量仪对测量点进行定时检测；步骤2，建立其有限元模型；步骤2‑1，根据目标数控机床进给系统的几何尺寸，建立其有限元模型，模型包括伺服电机、两个支撑轴承、丝杠、螺母和两个滑轨，各部分均选择固体单元，而丝杠与两个支撑轴承、丝母之间的连接选用接触单元；步骤2‑2，确定系统的热源；系统包括固定热源和移动热源；固定热源为伺服电机、轴承一和轴承二；依据系统结构，将伺服电机发热率合并到轴承一，记为QB1，轴承二的发热率记为QB2；移动热源为两个滑轨和一个丝母，两个滑轨发热率为Qg1＝Qg2＝Qg，丝母的发热率Qn；步骤2‑3，对流换热系数的确定；在丝杠旋转过程，丝杠表面与空气产生对流传热，则对流换热系数为：

【技术特征摘要】
1.一种数控机床滚珠丝杠进给系统热误差预测方法，其特征在于，包括如下步骤：步骤1，滚珠丝杠进给系统温度测量确定进给系统表面温度检测点，其中，温度检测点包括两个丝杠支撑轴承的表面点、丝母法兰的表面点和丝杠轴一侧导轨表面均匀分布点；T型贴片热电偶通过磁性底座分别固定在两个丝杠支撑轴承的表面点、丝母法兰的表面点上，按周期进行连续采样；通过红外线测量仪对测量点进行定时检测；步骤2，建立其有限元模型；步骤2-1，根据目标数控机床进给系统的几何尺寸，建立其有限元模型，模型包括伺服电机、两个支撑轴承、丝杠、螺母和两个滑轨，各部分均选择固体单元，而丝杠与两个支撑轴承、丝母之间的连接选用接触单元；步骤2-2，确定系统的热源；系统包括固定热源和移动热源；固定热源为伺服电机、轴承一和轴承二；依据系统结构，将伺服电机发热率合并到轴承一，记为QB1，轴承二的发热率记为QB2；移动热源为两个滑轨和一个丝母，两个滑轨发热率为Qg1＝Qg2＝Qg，丝母的发热率Qn；步骤2-3，对流换热系数的确定；在丝杠旋转过程，丝杠表面与空气产生对流传热，则对流换热系数为：其中，Nu为努塞尔数；λfluid为空气的导热系数；d为丝杠的直径；步骤3，用蒙特卡洛方法实现发热率的有限元计算；步骤4，运动副中心与表面温度关系的确定描述测量表面点和运动副对中心温度差：ΔT(t,v)＝av(1-e-t/bv)(8)其中，v是工作台的进给速度，t是工作时间，a和b是通过曲线拟合得到的常数；在机床工作过程中，检测支撑轴承或丝母表面的温度，计算得到其中心温度：Tc＝Tm+ΔT(t,v)(9)其中，Tc运动对的中心温度，Tm是其表面上给定点的测量温度；步骤5，滚珠丝杠实时热误差预测；步骤5-1建立滚珠丝杠热传导模型将丝杠简化为一维杆，其热传导方程：其中，T(x,t)是时间t和位置x的函数，表示t时刻丝杠上距离热源为x的位置点的温度变化；κ是导热系数；ρ是滚珠丝杠的材料密度；c是比热容；h是对流系数；Tf(t)是工作环境空气的温度；对滚珠丝杠T(x,t)的函数划分空间步长s和时间步长τ，xk＝k·s，k＝1,2,…,M，tj...

【专利技术属性】
技术研发人员：赵春雨，李朕均，李铁军，张义民，闻邦椿，
申请(专利权)人：东北大学，
类型：发明
国别省市：辽宁,21