数控机床热误差补偿的建模方法
数控机床热误差补偿的建模方法
【技术领域】
[0001] 本发明属于精密机床加工技术领域，尤其涉及一种数控机床热误差补偿的建模方 法。
【背景技术】
[0002] 工件的加工精度取决于刀具和工件之间相对移动的精度。如果机床连接部件因温 度变化而产生变形，则会影响到零件的加工精度，通常我们把由热变形引起的机床误差成 为热误差。研宄表明，热误差占到整个机床定位误差的65%-75%，因此对机床热误差进行 补偿是提高机床精度的一种可行办法。
[0003]目前，较为常用的热误差补偿方法有限元建模补偿法、多元线性回归方法等，但是 有限元建模法由于边界条件确定和获得精确的热传递特性困难的因素应用受到一定的限 制；线性回归的方法建模容易但其鲁棒性差。近些年，虽然神经网络、模糊推理系统、自适应 神经网络模糊推理系统（ANFIS)及以上几种方法的复合已经成功应用于热误差补偿系统 的建模，但其训练网络时的数据过多造成对训练时间增加，且假如温度测点的位置和数量 没有合理选取，还会照成系统的鲁棒性差，ANFIS随着输入点数目的增加，模糊规则数目会 指数级的递增，对计算机的要求会大幅度地提高。
【发明内容】
[0004] 本发明所要解决的技术问题是：提供一种可使训练时间短、鲁棒性好的数控机床 热误差补偿的建模方法。
[0005] 为解决上述技术问题本发明所采用的技术方案是：数控机床热误差补偿的建模方 法，步骤为：
[0006]A、在机床关键位置布置n个温度传感器，形成n个数据点{Xi, x2,…，xn}，在机 床定位误差测量处设置位移传感器；
[0007] B、采集温度信号和位移信号，得到n个温度传感器测得的温度T随时间t的变化 量{ (t)，T2 (t)，…，Tn (t) }，以及机床定位误差量Y(t);
[0008]C、将所测数据{ ^⑴心⑴，…，Tn(t) }以及Y(t)进行减法聚类处理：
[0009]对于n个数据点{Xpx2,…，xn }，数据点xj处的密度指标定义为：
[0011] 其中Ya是一个正数，常数ya定义了数据点Xj的一个邻域；
[0012] 在计算每个数据点的密度指标后，选择具有最高密度指标的数据点作为第一个聚 类中心，令x。^选中的第一个聚类中心，Mel为其密度指标，那么每个数据点Xi的密度指标 可以用公式：
[0014] 进行修正，其中Mik为修正后的密度指标；Yb是一个正数，常数Yb定义了一个密 度指标显著减小的邻域，且Yb大于Ya;
[0015] 修正了每个数据点的密度指标后，选定下一个聚类中心xe2，再次修正数据点的所 有密度指标，该过程不断反复，直到产生需要数量的聚类中心。
[0016] D、经减法聚类处理后得到的聚类中心数量为c个，通过FCM聚类计算新的聚类中 心点，确定新的聚类中心点为建模使用的热关键点，所采用的FCM聚类算法为：
[0017] 用公式3计算新的隶属矩阵U，
[0019] 其中，Uu的值介于0和1之间，且满足公式4要求：
[0021] 根据公式5计算价值函数，
[0023]其中(^为模糊组i的聚类中心，du= | |ci_Xj | |为第i个聚类中心与第j个数据 点间的欧几里德距离，me[1，<-)为加权指数，然后根据公式6构造新的目标函数，可求 得使公式5达到最小值的必要条件：
[0025]其中A』是公式4中的n个约束式的拉格朗日乘子，然后对所有输入参量求导，使 公式5达到最小值的必要条件为：
[0027] 根据公式5计算价值函数J，若J小于某个确定的阈值或J相对于上次价值函数值 的改变量小于某个阈值，则算法停止，否则用公式3计算新的隶属矩阵U，循环往复可得到 优化后的热关键点；
[0028] E、利用优化后的热关键点对热误差补偿系统进行建模。
[0029] 进一步的是：步骤E中对热误差补偿系统进行建模时，采用灰色系统进行建模，根 据灰模型得到的热误差预测值进行补偿；
[0030] 灰色系统建模方法为：
[0031] 设置一阶的灰模型GM(1，N)，其具有N个变量Xi(i= 1，2，…，N)，每个变量有n 个初始序列值
[0033] 进行一阶累加运算后序列为：
[0036] 灰模型可表达为：
[0038] 其中zf'幻定义为：
[0040] 分别称系数a和bj为系统发展参数和驱动参数；
[0041] 参数列丨义根据最小二乘法可得到：
[0043] 其中
[0045] 灰模型GM(1，N)的近似时间响应式为：
[0048] 通过公式10得出4(1\序列后，通过公式11进行一次累减还原，可得到GM(1，N)模 型对于热误差的预测值，
[0050] 本发明的有益效果是：经过减法聚类处理后，可以获得初步分组，可以减少FCM聚 类处理的运算时间，提高系统的响应性能，同时经过FCM聚类处理后，具有更好的鲁棒性； 利用减法聚类处理后，不会出现在单独使用FCM聚类处理时出现的局部最小点，使系统具 有更高的精度；本发明的方法简便易行，可以很方便的移植到其他的数控机床上，可以广泛 应用于数控机床的热误差补偿建模。
【附图说明】
[0051] 图1是本发明的工作流程图；
[0052] 图2是FCM聚类处理的流程图；
[0053] 图3是本发明实施后的补偿效果图。
【具体实施方式】
[0054] 下面结合附图对本发明作进一步说明，文中提到的"FCM聚类"全称为"模糊C平 均聚类"。
[0055]如图1、图2所示，本发明步骤为：
[0056] A、在机床关键位置布置n个温度传感器，形成n个数据点{Xp x2,…，xn}，在机 床定位误差测量处设置位移传感器；布置温度传感器时，应该在机床关键位置及其附近布 置尽量多的温度传感器，以保证样本空间足够大；
[0057] B、采集温度信号和位移信号，得到n个温度传感器测得的温度T随时间t的变化 量{ (t)，T2 (t)，…，Tn (t) }，以及机床定位误差量Y(t);
[0058]C、将所测数据{TiahTjt)，…，Tn(t) }以及Y(t)进行减法聚类处理：
[0059]对于n个数据点{Xp x2,…，xn}，数据点x j处的密度指标定义为：
[0061] 其中Ya是一个正数，常数Ya定义了数据点X」的一个邻域；显然如果一个数据点 若有多个临近的数据点，则该数据点具有高密度值，半径73以外的数据点对该点的密度指 标贡献甚微。
[0062] 在计算每个数据点的密度指标后，选择具有最高密度指标的数据点作为第一个聚 类中心，令x。^选中的第一个聚类中心，Mel为其密度指标，那么每个数据点Xi的密度指标 可以用公式：
[0064] 进行修正，其中Mik为修正后的密度指标；yb是一个正数，显然靠近第一个聚类中 心Xc;1的数据点的密度指标将显著减小，这样使得这些点不太可能成为下一个聚类中心。常 数yb定义了一个密度指标显著减小的邻域，yb应大于ya，以避免出现相聚很近的聚类中 心。
[0065] 修正了每个数据点的密度指标后，选定下一个聚 类中心xe2，再次修正数据点的所 有密度指标，该过程不断反复，直到产生需要数量的聚类中心。
[0066] D、经减法聚类处理后得到的聚类中心数量为c个，通过FCM聚类计算新的聚类中 心点，确定新的聚类中心点为建模使用的热关键点，所采用的FCM聚类算法为：
[0067] 用公式3计算新的隶属矩阵U，
[0069] 其中，Uu的值介于0和1之间，且满足公式4要求：
[0071] 根据公式5计算价值函数，
[0073]其中(^为模糊组i的聚类中心，du= | |ci_Xj | |为第i个聚类中心与第j个数据 点间的欧几里德距离，me[1，<-)为加权指数，然后根据公式6构造新的目标函数，可求 得使公式5达到最小值的必要条件：
[0075]其中A』是公式4中的n个约束式的拉格朗日乘子，然后对所有输入参量求导，使 公式5达到最小值的必要条件为：
[0077] 根据公式5计算价值函数J，若J小于某个确定的阈值或J相对于上次价值函数值 的改变量小于某个阈值，则算法停止，否则用公式3计算新的隶属矩阵U，循环往复可得到 优化后的热关键点；FCM聚类算法具体流程如图2所示；
[0078] E、利用优化后的热关键点对热误差补偿系统进行建模。
[0079] 本发明可应用现有技术中常见的神经网络、ANFIS系统等进行建模，作为对本发明 的改进，步骤E中对热误差补偿系统进行建模时，采用灰色系统进行建模，根据灰模型得到 的热误差预测值进行补偿，可使热误差补偿系统具有更高的精度和鲁棒性。
[0080] 灰色系统建模方法为：
[0081] 设置一阶的灰模型GM(1，N)，其为一个多变量（即多个输入量）的灰模型，GM(1， N)意味着这个灰模型包含一个因变量和N-1个自变量，即具有N个变量Xi(i= 1，2，…， N)，每个变量有n个初始序列值
[0083] 为了使序列有规律性，进行累加生成运算，为了简化，进行一阶累加运算后序列 为：
[0086] 灰模型可表达为：
[0088] 其中zf1)(幻定义为：
[0090] 分别称系数a和bj为系统发展参数和驱动参数；
[0091] 参数列丨表根据最小二乘法可得到：
[0093] 其中
[0095] 灰模型GM(1，N)的近似时间响应式为：
[0098] 通过公式10得出if1}序列后，通过公式11进行一次累减还原，可得到GM(1，N)模 型对于热误差的预测值，
[0100] 本发明采用经减法聚类处理和FCM聚类处理得到的热关键点，以及灰色系统建模 技术实施，补偿效果如图3所示。
[0101] 本发明经减法聚类处理后再进行FCM聚类，避免了单独使用FCM聚类时，随机输入 初始中心点，耽误了大量的处理时间，通过减法聚类，可以在最短时间内将类分好，确定好 聚类中心，然后把预处理好的分类和聚类中心交给FCM聚类处理，从而减少处理时间；另外 一方面使用减法聚类极好的解决了FCM聚类中的收敛性能严重依赖于聚类的初始点问题， FCM聚类的收敛速度随初始点选择的不同有很大变化，且可能陷入局部极小点，从而影响精 密数控机床制造精度与加工效率。本发明具有更高的精度和鲁棒性，方法简便易行，可以很 方便的移植到其他的数控机床上，可以广泛应用于数控机床的热误差补偿建模。
【主权项】
1.数控机床热误差补偿的建模方法，其特征在于，步骤为： A、在机床关键位置布置η个温度传感器，形成η个数据点{ X1, X2,…，Xn }，在机床定 位误差测量处设置位移传感器； Β、采集温度信号和位移信号，得到η个温度传感器测得的温度T随时间t的变化量{ T1 (t)，T2 (t)，…，Tn (t) }，以及机床定位误差量Y (t); C、 将所测数据（}以及Y(t)进行减法聚类处理： 对于η个数据点{ X1, X2,…，xn }，数据点X j处的密度指标定义为：其中Ya是一个正数，常数γ a定义了数据点h的一个邻域； 在计算每个数据点的密度指标后，选择具有最高密度指标的数据点作为第一个聚类中 心，令Xc;1为选中的第一个聚类中心，Mcl为其密度指标，那么每个数据点^的密度指标可以 用公式：进行修正，其中Mik为修正后的密度指标；γ b是一个正数，常数γ b定义了一个密度指 标显著减小的邻域，且Yb大于γ a; 修正了每个数据点的密度指标后，选定下一个聚类中心xe2，再次修正数据点的所有密 度指标，该过程不断反复，直到产生需要数量的聚类中心。 D、 经减法聚类处理后得到的聚类中心数量为c个，通过FCM聚类计算新的聚类中心点， 确定新的聚类中心点为建模使用的热关键点，所采用的FCM聚类算法为： 用公式3计算新的隶属矩阵U，其中，Uu的值介于O和1之间，且满足公式4要求：根据公式5计算价值函数，其中CiS模糊组i的聚类中心，d ij= I I c i-Xj I I为第i个聚类中心与第j个数据点间 的欧几里德距离，me [1，m)为加权指数，然后根据公式6构造新的目标函数，可求得使 公式5达到最小值的必要条件：其中\是公式4中的η个约束式的拉格朗日乘子，然后对所有输入参量求导，使公式 5达到最小值的必要条件为：根据公式5计算价值函数J，若J小于某个确定的阈值或J相对于上次价值函数值的改 变量小于某个阈值，则算法停止，否则用公式3计算新的隶属矩阵U，循环往复可得到优化 后的热关键点； Ε、利用优化后的热关键点对热误差补偿系统进行建模。2.如权利要求1所述的数控机床热误差补偿的建模方法，其特征在于：步骤E中对热 误差补偿系统进行建模时，采用灰色系统进行建模，根据灰模型得到的热误差预测值进行 补偿； 灰色系统建模方法为： 设置一阶的灰模型GM(1，Ν)，其具有N个变量Xi (i = 1，2，…，Ν)，每个变量有η个初 始序列值进行一阶累加运算后序列为：灰模型可表达为：分别称系数a和h为系统发展参数和驱动参数； 参数列根据最小二乘法可得到： Θ^{Β'Β) 'B1Y (9), 其中灰模型GM(1，Ν)的近似时间响应式为：通过公式10得出序列后，通过公式11进行一次累减还原，可得到GM(1，N)模型对 于热误差的预测值，
【专利摘要】本发明涉及一种数控机床热误差补偿的建模方法，步骤为：A、在机床关键位置布置n个温度传感器，形成n个数据点﹛x1,x2,…,xn﹜,在机床定位误差测量处设置位移传感器；B、采集温度信号和位移信号，得到n个温度传感器测得的温度T随时间t的变化量﹛T1(t),T2(t),…,Tn(t)﹜，以及机床定位误差量Y(t)；C、将所测数据﹛T1(t),T2(t),…,Tn(t)﹜以及Y(t)进行减法聚类处理；D、经减法聚类处理后得到的聚类中心数量为c个，通过FCM聚类计算新的聚类中心点，确定新的聚类中心点为建模使用的热关键点；E、利用优化后的热关键点对热误差补偿系统进行建模。本发明具有更高的精度和鲁棒性，方法简便易行，可以广泛应用于数控机床的热误差补偿建模。
【IPC分类】G05B19/404
【公开号】CN104898562
【申请号】CN201510204608
【发明人】魏弦
【申请人】攀枝花学院
【公开日】2015年9月9日
【申请日】2015年4月27日
【技术领域】
[0001] 本发明属于精密机床加工技术领域，尤其涉及一种数控机床热误差补偿的建模方 法。
【背景技术】
[0002] 工件的加工精度取决于刀具和工件之间相对移动的精度。如果机床连接部件因温 度变化而产生变形，则会影响到零件的加工精度，通常我们把由热变形引起的机床误差成 为热误差。研宄表明，热误差占到整个机床定位误差的65%-75%，因此对机床热误差进行 补偿是提高机床精度的一种可行办法。
[0003]目前，较为常用的热误差补偿方法有限元建模补偿法、多元线性回归方法等，但是 有限元建模法由于边界条件确定和获得精确的热传递特性困难的因素应用受到一定的限 制；线性回归的方法建模容易但其鲁棒性差。近些年，虽然神经网络、模糊推理系统、自适应 神经网络模糊推理系统（ANFIS)及以上几种方法的复合已经成功应用于热误差补偿系统 的建模，但其训练网络时的数据过多造成对训练时间增加，且假如温度测点的位置和数量 没有合理选取，还会照成系统的鲁棒性差，ANFIS随着输入点数目的增加，模糊规则数目会 指数级的递增，对计算机的要求会大幅度地提高。
【发明内容】
[0004] 本发明所要解决的技术问题是：提供一种可使训练时间短、鲁棒性好的数控机床 热误差补偿的建模方法。
[0005] 为解决上述技术问题本发明所采用的技术方案是：数控机床热误差补偿的建模方 法，步骤为：
[0006]A、在机床关键位置布置n个温度传感器，形成n个数据点{Xi, x2,…，xn}，在机 床定位误差测量处设置位移传感器；
[0007] B、采集温度信号和位移信号，得到n个温度传感器测得的温度T随时间t的变化 量{ (t)，T2 (t)，…，Tn (t) }，以及机床定位误差量Y(t);
[0008]C、将所测数据{ ^⑴心⑴，…，Tn(t) }以及Y(t)进行减法聚类处理：
[0009]对于n个数据点{Xpx2,…，xn }，数据点xj处的密度指标定义为：
[0011] 其中Ya是一个正数，常数ya定义了数据点Xj的一个邻域；
[0012] 在计算每个数据点的密度指标后，选择具有最高密度指标的数据点作为第一个聚 类中心，令x。^选中的第一个聚类中心，Mel为其密度指标，那么每个数据点Xi的密度指标 可以用公式：
[0014] 进行修正，其中Mik为修正后的密度指标；Yb是一个正数，常数Yb定义了一个密 度指标显著减小的邻域，且Yb大于Ya;
[0015] 修正了每个数据点的密度指标后，选定下一个聚类中心xe2，再次修正数据点的所 有密度指标，该过程不断反复，直到产生需要数量的聚类中心。
[0016] D、经减法聚类处理后得到的聚类中心数量为c个，通过FCM聚类计算新的聚类中 心点，确定新的聚类中心点为建模使用的热关键点，所采用的FCM聚类算法为：
[0017] 用公式3计算新的隶属矩阵U，
[0019] 其中，Uu的值介于0和1之间，且满足公式4要求：
[0021] 根据公式5计算价值函数，
[0023]其中(^为模糊组i的聚类中心，du= | |ci_Xj | |为第i个聚类中心与第j个数据 点间的欧几里德距离，me[1，<-)为加权指数，然后根据公式6构造新的目标函数，可求 得使公式5达到最小值的必要条件：
[0025]其中A』是公式4中的n个约束式的拉格朗日乘子，然后对所有输入参量求导，使 公式5达到最小值的必要条件为：
[0027] 根据公式5计算价值函数J，若J小于某个确定的阈值或J相对于上次价值函数值 的改变量小于某个阈值，则算法停止，否则用公式3计算新的隶属矩阵U，循环往复可得到 优化后的热关键点；
[0028] E、利用优化后的热关键点对热误差补偿系统进行建模。
[0029] 进一步的是：步骤E中对热误差补偿系统进行建模时，采用灰色系统进行建模，根 据灰模型得到的热误差预测值进行补偿；
[0030] 灰色系统建模方法为：
[0031] 设置一阶的灰模型GM(1，N)，其具有N个变量Xi(i= 1，2，…，N)，每个变量有n 个初始序列值
[0033] 进行一阶累加运算后序列为：
[0036] 灰模型可表达为：
[0038] 其中zf'幻定义为：
[0040] 分别称系数a和bj为系统发展参数和驱动参数；
[0041] 参数列丨义根据最小二乘法可得到：
[0043] 其中
[0045] 灰模型GM(1，N)的近似时间响应式为：
[0048] 通过公式10得出4(1\序列后，通过公式11进行一次累减还原，可得到GM(1，N)模 型对于热误差的预测值，
[0050] 本发明的有益效果是：经过减法聚类处理后，可以获得初步分组，可以减少FCM聚 类处理的运算时间，提高系统的响应性能，同时经过FCM聚类处理后，具有更好的鲁棒性； 利用减法聚类处理后，不会出现在单独使用FCM聚类处理时出现的局部最小点，使系统具 有更高的精度；本发明的方法简便易行，可以很方便的移植到其他的数控机床上，可以广泛 应用于数控机床的热误差补偿建模。
【附图说明】
[0051] 图1是本发明的工作流程图；
[0052] 图2是FCM聚类处理的流程图；
[0053] 图3是本发明实施后的补偿效果图。
【具体实施方式】
[0054] 下面结合附图对本发明作进一步说明，文中提到的"FCM聚类"全称为"模糊C平 均聚类"。
[0055]如图1、图2所示，本发明步骤为：
[0056] A、在机床关键位置布置n个温度传感器，形成n个数据点{Xp x2,…，xn}，在机 床定位误差测量处设置位移传感器；布置温度传感器时，应该在机床关键位置及其附近布 置尽量多的温度传感器，以保证样本空间足够大；
[0057] B、采集温度信号和位移信号，得到n个温度传感器测得的温度T随时间t的变化 量{ (t)，T2 (t)，…，Tn (t) }，以及机床定位误差量Y(t);
[0058]C、将所测数据{TiahTjt)，…，Tn(t) }以及Y(t)进行减法聚类处理：
[0059]对于n个数据点{Xp x2,…，xn}，数据点x j处的密度指标定义为：
[0061] 其中Ya是一个正数，常数Ya定义了数据点X」的一个邻域；显然如果一个数据点 若有多个临近的数据点，则该数据点具有高密度值，半径73以外的数据点对该点的密度指 标贡献甚微。
[0062] 在计算每个数据点的密度指标后，选择具有最高密度指标的数据点作为第一个聚 类中心，令x。^选中的第一个聚类中心，Mel为其密度指标，那么每个数据点Xi的密度指标 可以用公式：
[0064] 进行修正，其中Mik为修正后的密度指标；yb是一个正数，显然靠近第一个聚类中 心Xc;1的数据点的密度指标将显著减小，这样使得这些点不太可能成为下一个聚类中心。常 数yb定义了一个密度指标显著减小的邻域，yb应大于ya，以避免出现相聚很近的聚类中 心。
[0065] 修正了每个数据点的密度指标后，选定下一个聚 类中心xe2，再次修正数据点的所 有密度指标，该过程不断反复，直到产生需要数量的聚类中心。
[0066] D、经减法聚类处理后得到的聚类中心数量为c个，通过FCM聚类计算新的聚类中 心点，确定新的聚类中心点为建模使用的热关键点，所采用的FCM聚类算法为：
[0067] 用公式3计算新的隶属矩阵U，
[0069] 其中，Uu的值介于0和1之间，且满足公式4要求：
[0071] 根据公式5计算价值函数，
[0073]其中(^为模糊组i的聚类中心，du= | |ci_Xj | |为第i个聚类中心与第j个数据 点间的欧几里德距离，me[1，<-)为加权指数，然后根据公式6构造新的目标函数，可求 得使公式5达到最小值的必要条件：
[0075]其中A』是公式4中的n个约束式的拉格朗日乘子，然后对所有输入参量求导，使 公式5达到最小值的必要条件为：
[0077] 根据公式5计算价值函数J，若J小于某个确定的阈值或J相对于上次价值函数值 的改变量小于某个阈值，则算法停止，否则用公式3计算新的隶属矩阵U，循环往复可得到 优化后的热关键点；FCM聚类算法具体流程如图2所示；
[0078] E、利用优化后的热关键点对热误差补偿系统进行建模。
[0079] 本发明可应用现有技术中常见的神经网络、ANFIS系统等进行建模，作为对本发明 的改进，步骤E中对热误差补偿系统进行建模时，采用灰色系统进行建模，根据灰模型得到 的热误差预测值进行补偿，可使热误差补偿系统具有更高的精度和鲁棒性。
[0080] 灰色系统建模方法为：
[0081] 设置一阶的灰模型GM(1，N)，其为一个多变量（即多个输入量）的灰模型，GM(1， N)意味着这个灰模型包含一个因变量和N-1个自变量，即具有N个变量Xi(i= 1，2，…， N)，每个变量有n个初始序列值
[0083] 为了使序列有规律性，进行累加生成运算，为了简化，进行一阶累加运算后序列 为：
[0086] 灰模型可表达为：
[0088] 其中zf1)(幻定义为：
[0090] 分别称系数a和bj为系统发展参数和驱动参数；
[0091] 参数列丨表根据最小二乘法可得到：
[0093] 其中
[0095] 灰模型GM(1，N)的近似时间响应式为：
[0098] 通过公式10得出if1}序列后，通过公式11进行一次累减还原，可得到GM(1，N)模 型对于热误差的预测值，
[0100] 本发明采用经减法聚类处理和FCM聚类处理得到的热关键点，以及灰色系统建模 技术实施，补偿效果如图3所示。
[0101] 本发明经减法聚类处理后再进行FCM聚类，避免了单独使用FCM聚类时，随机输入 初始中心点，耽误了大量的处理时间，通过减法聚类，可以在最短时间内将类分好，确定好 聚类中心，然后把预处理好的分类和聚类中心交给FCM聚类处理，从而减少处理时间；另外 一方面使用减法聚类极好的解决了FCM聚类中的收敛性能严重依赖于聚类的初始点问题， FCM聚类的收敛速度随初始点选择的不同有很大变化，且可能陷入局部极小点，从而影响精 密数控机床制造精度与加工效率。本发明具有更高的精度和鲁棒性，方法简便易行，可以很 方便的移植到其他的数控机床上，可以广泛应用于数控机床的热误差补偿建模。
【主权项】
1.数控机床热误差补偿的建模方法，其特征在于，步骤为： A、在机床关键位置布置η个温度传感器，形成η个数据点{ X1, X2,…，Xn }，在机床定 位误差测量处设置位移传感器； Β、采集温度信号和位移信号，得到η个温度传感器测得的温度T随时间t的变化量{ T1 (t)，T2 (t)，…，Tn (t) }，以及机床定位误差量Y (t); C、 将所测数据（}以及Y(t)进行减法聚类处理： 对于η个数据点{ X1, X2,…，xn }，数据点X j处的密度指标定义为：其中Ya是一个正数，常数γ a定义了数据点h的一个邻域； 在计算每个数据点的密度指标后，选择具有最高密度指标的数据点作为第一个聚类中 心，令Xc;1为选中的第一个聚类中心，Mcl为其密度指标，那么每个数据点^的密度指标可以 用公式：进行修正，其中Mik为修正后的密度指标；γ b是一个正数，常数γ b定义了一个密度指 标显著减小的邻域，且Yb大于γ a; 修正了每个数据点的密度指标后，选定下一个聚类中心xe2，再次修正数据点的所有密 度指标，该过程不断反复，直到产生需要数量的聚类中心。 D、 经减法聚类处理后得到的聚类中心数量为c个，通过FCM聚类计算新的聚类中心点， 确定新的聚类中心点为建模使用的热关键点，所采用的FCM聚类算法为： 用公式3计算新的隶属矩阵U，其中，Uu的值介于O和1之间，且满足公式4要求：根据公式5计算价值函数，其中CiS模糊组i的聚类中心，d ij= I I c i-Xj I I为第i个聚类中心与第j个数据点间 的欧几里德距离，me [1，m)为加权指数，然后根据公式6构造新的目标函数，可求得使 公式5达到最小值的必要条件：其中\是公式4中的η个约束式的拉格朗日乘子，然后对所有输入参量求导，使公式 5达到最小值的必要条件为：根据公式5计算价值函数J，若J小于某个确定的阈值或J相对于上次价值函数值的改 变量小于某个阈值，则算法停止，否则用公式3计算新的隶属矩阵U，循环往复可得到优化 后的热关键点； Ε、利用优化后的热关键点对热误差补偿系统进行建模。2.如权利要求1所述的数控机床热误差补偿的建模方法，其特征在于：步骤E中对热 误差补偿系统进行建模时，采用灰色系统进行建模，根据灰模型得到的热误差预测值进行 补偿； 灰色系统建模方法为： 设置一阶的灰模型GM(1，Ν)，其具有N个变量Xi (i = 1，2，…，Ν)，每个变量有η个初 始序列值进行一阶累加运算后序列为：灰模型可表达为：分别称系数a和h为系统发展参数和驱动参数； 参数列根据最小二乘法可得到： Θ^{Β'Β) 'B1Y (9), 其中灰模型GM(1，Ν)的近似时间响应式为：通过公式10得出序列后，通过公式11进行一次累减还原，可得到GM(1，N)模型对 于热误差的预测值，
【专利摘要】本发明涉及一种数控机床热误差补偿的建模方法，步骤为：A、在机床关键位置布置n个温度传感器，形成n个数据点﹛x1,x2,…,xn﹜,在机床定位误差测量处设置位移传感器；B、采集温度信号和位移信号，得到n个温度传感器测得的温度T随时间t的变化量﹛T1(t),T2(t),…,Tn(t)﹜，以及机床定位误差量Y(t)；C、将所测数据﹛T1(t),T2(t),…,Tn(t)﹜以及Y(t)进行减法聚类处理；D、经减法聚类处理后得到的聚类中心数量为c个，通过FCM聚类计算新的聚类中心点，确定新的聚类中心点为建模使用的热关键点；E、利用优化后的热关键点对热误差补偿系统进行建模。本发明具有更高的精度和鲁棒性，方法简便易行，可以广泛应用于数控机床的热误差补偿建模。
【IPC分类】G05B19/404
【公开号】CN104898562
【申请号】CN201510204608
【发明人】魏弦
【申请人】攀枝花学院
【公开日】2015年9月9日
【申请日】2015年4月27日
最新回复(0)