CN114237154A - 基于ont-gcn时空模型的热误差预测模型及建模方法和霾-边-雾-云误差补偿系统 - Google Patents

基于ont-gcn时空模型的热误差预测模型及建模方法和霾-边-雾-云误差补偿系统 Download PDF

Info

Publication number
CN114237154A
CN114237154A CN202111421117.3A CN202111421117A CN114237154A CN 114237154 A CN114237154 A CN 114237154A CN 202111421117 A CN202111421117 A CN 202111421117A CN 114237154 A CN114237154 A CN 114237154A
Authority
CN
China
Prior art keywords
error
representing
thermal
axis
coordinate system
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202111421117.3A
Other languages
English (en)
Other versions
CN114237154B (zh
Inventor
马驰
刘佳兰
桂洪泉
王时龙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing University
Original Assignee
Chongqing University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Chongqing University filed Critical Chongqing University
Priority to CN202111421117.3A priority Critical patent/CN114237154B/zh
Publication of CN114237154A publication Critical patent/CN114237154A/zh
Application granted granted Critical
Publication of CN114237154B publication Critical patent/CN114237154B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/18Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
    • G05B19/404Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by control arrangements for compensation, e.g. for backlash, overshoot, tool offset, tool wear, temperature, machine construction errors, load, inertia
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/33Director till display
    • G05B2219/33133For each action define function for compensation, enter parameters

Landscapes

  • Engineering & Computer Science (AREA)
  • Human Computer Interaction (AREA)
  • Manufacturing & Machinery (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Numerical Control (AREA)

Abstract

本发明公开了一种基于ONT‑GCN时空模型的热误差预测模型及其建模方法和霾‑边‑雾‑云误差补偿系统。本发明基于ONT‑GCN时空模型的热误差预测模型,利用LSTMN神经网络捕捉传感器收集的数据的时间依赖性,利用GCN神经网络捕捉拓扑结构的空间特征,可以将捕获的空间特征和时序特征保留在ONT‑GCN单元中;ONT‑GCN单元的独特排序特性允许保留重要的热误差信息,因此,ONT‑GCN单元可以有序地传递热误差的时间和空间特征,从而提高热误差的预测精度。本发明的霾‑边‑雾‑云误差补偿系统,云计算具有强大的计算能力,用于解决耗时问题;通过设置雾计算层和边缘层等分布式计算层,可缓解工业互联网的带宽压力;通过霾计算层以加快信息获取过程,从而以保证系统执行效率和机床的加工精度。

Description

基于ONT-GCN时空模型的热误差预测模型及建模方法和霾- 边-雾-云误差补偿系统
技术领域
本发明属于机械误差分析技术领域,具体的为一种基于ONT-GCN时空模型的热误差预测模型及霾-边-雾-云误差补偿系统。
背景技术
机械传动是装备制造业的基础,蜗轮是机械传动中最常见和最基本的零件之一。蜗轮的几何精度受加工系统多源误差的影响。几何误差和热误差是两个最重要的误差。对于蜗轮加工机床来说,影响热误差的热源很多,包括电机、轴承、滚动导轨、滚珠丝杠等。主轴是对热误差影响显着的关键部件。减少热误差的方法主要有以下三种:(1)温度控制方法;(2)热误差预防方法;(3)热误差补偿方法。温度控制方法通过控制热误差来减少热误差,但实施成本高,技术实施难度也高;温度传递的滞后导致温度控制的滞后,所以它不是减少热误差的最佳选择。热误差预防方法在设计和制造阶段选择合理的材料和机床结构设计,从而降低对热流的误差敏感性,提高热稳定性,然而,该方法的成本随着控制精度的提高呈指数增长,因此,它不是常用的措施。热误差补偿方法是探索热误差与其他变量之间的关系,因此,热误差补偿方法是实际应用中最实用的措施,因其简单、高效、低成本的独特特点而被广泛使用。
为了表征了热误差的变化过程,提出了多种数学模型,主要有数值模拟模型和经验模型,传统的经验模型基于统计方法,无法挖掘潜在关系,缺乏数据特征的提取能力。近来的一些经验模型使用了机器学习技术,与传统的经验模型不同,考虑了误差产生机理,说明了热误差是运行时间的函数。采用具有强记忆能力的长短期记忆网络(LSTMN)进行建模。结果表明,LSTMN是一个出色的模型,这一结论也反映在机床热误差建模中。但误差数据的时间序列分析导致传感器的空间信息被忽略,从而热误差数据不受传感器网络的限制。
经验误差模型的训练非常耗时。云计算具有强大的计算能力,被广泛用于解决耗时问题,但由于工业互联网高延迟和有限带宽等问题,云计算不能直接应用于误差控制,必须有雾计算或边缘计算等分布式计算来缓解工业互联网的带宽压力,然而,边缘-雾-云架构中忽略了终端设备的设计和配置。
发明内容
有鉴于此,本发明的目的在于提供一种基于ONT-GCN时空模型的热误差预测模型及霾-边-雾-云误差补偿系统,能够有序地传递热误差的时间和空间特征,从而提高热误差的预测精度。
为达到上述目的,本发明提供如下技术方案:
本发明首先提出了一种基于ONT-GCN时空模型的热误差预测模型,所述ONT-GCN神经网络的表达式为:
ft=δ(Wf[ht-1,f(A,Xt)]+bf)
it=δ(Wi[ht-1,f(A,Xt)]+bi)
Figure BDA0003376641850000011
ot=δ(Wo[ht-1,f(A,Xt)]+bo)
Figure BDA0003376641850000021
Figure BDA0003376641850000022
cumax(·)=cumsum(softmax(·))
Figure BDA0003376641850000023
Figure BDA0003376641850000024
Figure BDA0003376641850000025
Figure BDA0003376641850000026
Figure BDA0003376641850000027
其中,ft表示遗忘门;it表示输入门;ot表示输出门;
Figure BDA0003376641850000028
表示主遗忘门;
Figure BDA0003376641850000029
表示主输入门;ht-1表示前一个单元的输出值,与f(A,Xt)结合增加三个门的权重,Xt表示特征序列,A表示邻接矩阵,f(A,Xt)表示Xt经过图卷积之后的序列;
Figure BDA00033766418500000210
表示t时刻临时单元的内存状态;ct-1表示t-1时刻临时单元的内存状态;ct表示内存在t时刻的单元状态;Wf、Wi、Wc和Wo表示权重矩阵;bf、bi、bc和bo表示偏置向量;σ表示激活函数;ht表示t时刻的隐藏状态;
Figure BDA00033766418500000223
Figure BDA00033766418500000224
分别表示
Figure BDA00033766418500000211
Figure BDA00033766418500000212
的权重矩阵;
Figure BDA00033766418500000213
Figure BDA00033766418500000214
分别表示
Figure BDA00033766418500000215
Figure BDA00033766418500000216
的偏置向量;comsum表示求和的函数;wt1、wt2和wt3分别表示高等级、中等级和低等级,其中,wt1*ct-1存储低等级信息,受当前输入信息的影响;
Figure BDA00033766418500000217
存储高等级信息,受历史信息影响;
Figure BDA00033766418500000218
存储中等级信息,受当前输入信息和历史信息影响;
f(A,Xt)表示Xt经过图卷积之后的序列,并作为ON-LSTMN模型的输入向量,且:
Figure BDA00033766418500000219
其中,
Figure BDA00033766418500000220
I表示单位矩阵;A表示邻接矩阵;
Figure BDA00033766418500000221
X表示特征序列;W0表示第一层的权重参数;σ表示激活函数;
热误差预测模型需要时空特性,因此,热误差的时空预测可以表示为传感器网络G与特征序列Xt之间映射关系f的学习,下一时刻T的热误差序列可以表示为:
[Xt+1,…,Xt+T]=f(G;(Xt-n,…,Xt-1,Xt))
其中,n表示历史时间序列的长度。
本发明还提出了一种基于ONT-GCN时空模型的热误差预测模型的建模方法,包括如下步骤:
11)构建传感器网络:在机床上布置温度传感器和位移传感器,并将各个传感器分别标记在机床的三维模型上;
12)生成时空图数据集:需要两个条件来预测当前节点的热误差信号,分别为:①位移测量点的历史热误差信号;②相邻传感器的位置和历史信号;
传感器网络的空间特征由邻接矩阵表征,邻接矩阵为:
Figure BDA00033766418500000222
其中,wi,j表示边权重,与节点i和节点j的位置有关;δ表示控制wi,j分布的阈值;di ,j为节点i和节点j之间的距离,是机床三维模型上测量得到的;
13)构建如上所述基于ONT-GCN时空模型的热误差预测模型;
14)以邻接矩阵和传感器的时间数据作为热误差预测模型的输入,预测热误差。
本发明还提出了一种霾-边-雾-云误差补偿系统,其特征在于:包括霾计算层、边缘计算层、雾计算层和云计算层;
所述霾计算层包括:
温度传感器:用于测量温度;
位移传感器:用于测量位移;
采集卡:用于采集所述温度传感器和位移传感器测量得到的温度信号和位移信号;
树莓派:所述树莓派上连接有用于将温度信号和位移信号转换为数字信息的A/D转换器;
网关:用于经A/D转换器转换得到的数据传输至所述边缘计算层;
所述边缘计算层包括用于存储经所述霾计算层传入的数据的微型数据中心、用于过滤数据的交换机和用于将过滤后的数据传输至雾计算层和云计算层的网关和路由器,其中,历史数据传输至云计算层,实时误差数据传输至雾计算层;
所述雾计算层内设有分析模块,所述分析模块内设有如上所述基于ONT-GCN时空模型的热误差预测模型,将以温度传感器、位移传感器的位置信息生成的邻接矩阵和以各个传感器采集的数据作为特征序列作为热误差预测模型的输入预测热误差,再结合误差映射模型得到误差补偿值,并将误差补偿值传输至云计算层;
所述云计算层包括误差补偿执行模块和热误差模型更新模块,所述误差补偿执行模块根据所述雾计算层输入的误差补偿值确定需要补偿的变量,并将补偿指令发送至霾计算层进行误差补偿;所述热误差模型更新模块根据边缘计算层输入的历史数据更新热误差预测模型,并将更新后的热误差预测模型传输至雾计算层。
本发明的有益效果在于:
本发明基于ONT-GCN时空模型的热误差预测模型,利用LSTMN神经网络捕捉传感器收集的数据的时间依赖性,利用GCN神经网络捕捉拓扑结构的空间特征,从而可以将捕获的空间特征和时序特征保留在ONT-GCN单元中;ONT-GCN单元的独特排序特性允许保留重要的热误差信息,因此,ONT-GCN单元可以有序地传递热误差的时间和空间特征,从而提高热误差的预测精度。
本发明的霾-边-雾-云误差补偿系统,云计算具有强大的计算能力,用于解决耗时问题,但由于工业互联网高延迟和有限带宽等问题,云计算不能直接应用于误差控制;通过设置雾计算层和边缘层等分布式计算层,可缓解工业互联网的带宽压力,但雾计算层和边缘层忽略了终端设备的设计和配置,因此设置了霾计算层以加快信息获取过程,霾计算层可以被认为是一个极端的边缘层,从而以保证系统执行效率和机床的加工精度。
附图说明
为了使本发明的目的、技术方案和有益效果更加清楚,本发明提供如下附图进行说明:
图1为本发明霾-边-雾-云误差补偿系统实施例的框架图;
图2为简化的一维主轴的结构示意图;
图3为LSTM神经网络的结构图;
图4为ON-LSTMN的结构图;
图5为ONT-GCN的时空预测过程;
图6为ONT-GCN单元的示意图;
图7为机床坐标系传动链的示意图;(a)为机床坐标系;(b)为坐标系传动链;
图8为涡轮蜗杆的加工示意图;
图9为X轴热误差示意图;
图10为Z轴热误差示意图;
图11为圆柱滚刀的渐开线侧面示意图;
图12为敏感性分析的流程图;
图13为齿面误差分量敏感度指标的曲线图;
图14为热误差时空预测原理图;
图15为温度传感器现场布置;
图16为蜗轮加工机床传感器网络;(a)为蜗轮蜗杆机床三维模型图;(b)为反映距离示意图;
图17为传感器邻接矩阵;
图18为800mm/min以下的温度场和热伸长率;(a)测量温度;(b)热误差;
图19为1500mm/min下的温度场和热伸长率;(a)测量温度;(b)热误差;
图20为预测性能对比的曲线图;
图21为蜗轮蜗杆加工实验图;
图22为齿面误差的示意图;
图23为齿面误差曲线图;(a)不执行霾-边-雾-云系统;(b)使用霾-边-雾-云系系统进行几何误差控制;(c)使用霾-边-雾-云系系统进行几何-热误差控制。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明,以使本领域的技术人员可以更好的理解本发明并能予以实施,但所举实施例不作为对本发明的限定。
如图1所示,为本发明霾-边-雾-云误差补偿系统实施例的框架图。本实施例的霾-边-雾-云误差补偿系统,包括霾计算层、边缘计算层、雾计算层和云计算层。本实施例以蜗轮蜗杆机床为例,对霾-边-雾-云误差补偿系统的具体实施方式进行详细说明。
所述霾计算层包括:温度传感器:用于测量温度;位移传感器:用于测量位移;采集卡:用于采集所述温度传感器和位移传感器测量得到的温度信号和位移信号;树莓派:所述树莓派上连接有用于将温度信号和位移信号转换为数字信息的A/D转换器;网关:用于经A/D转换器转换得到的数据传输至所述边缘计算层。具体的,霾计算层是终端,主要负责硬件放置、数据采集和数据预处理。温度Humirel Ni1000SOT和位移传感器Micro-EpsiloneddyDT3010-A布置在霾计算层。温度和位移传感器的放置非常重要,因为合理的温度传感器网络可以反映蜗轮加工机床的热状态,正确放置的位移传感器可以准确测量热误差。被测信号由树莓派4采集,再传输到A/D转换器PCF8591模块进行数据转换,具体的,A/D转换器安装在树莓派4上,是典型的终端设备。数据通过树莓派4到达网关,然后发送到边缘层,所有的计算层都由ESP8266Wi-Fi编程的Linux传输。
所述边缘计算层包括用于存储经所述霾计算层传入的数据的微型数据中心、用于过滤数据的交换机和用于将过滤后的数据传输至雾计算层和云计算层的网关和路由器,其中,历史数据传输至云计算层,实时误差数据传输至雾计算层。从雾层传输到边缘层的数据存储在微数据中心,数据通过交换机,其型号为西门子SCALANCE X-200,检测出异常值,如果数据很大,它将被滤波并到达网关和路由器。其中,历史数据直接发送到云计算层数据库进行存储,实时错误数据发送到雾计算层进行错误建模和分析。
所述雾计算层内设有分析模块,所述分析模块内设有如上所述基于ONT-GCN时空模型的热误差预测模型,将以温度传感器、位移传感器的位置信息生成的邻接矩阵和以各个传感器采集的数据作为特征序列作为热误差预测模型的输入预测热误差,再结合误差映射模型得到误差补偿值,并将误差补偿值传输至云计算层。ONT-GCN的建模过程在雾计算层完成。此外,ONT-GCN模型的鲁棒性可能会在霾-边-雾-云系统运行一段时间后下降。因此需要云计算层将重新训练的模型传输到雾计算层。ONT-GCN的输入有两个关键因素,即:(1)温度和位移传感器的位置信息用于生成邻接矩阵;(2)传感器采集的数据作为特征序列。数据从边缘计算层传输到雾计算层,根据蜗轮蜗杆机床的传递参数建立齿面误差映射模型,再结合齿面误差映射模型进行的灵敏度分析,确定关键误差项。如此就可以得到补偿值,传送到云计算层生成补偿指令。上述任务的执行平台为第11代Intel Core i5-1135G7,由CUDA编程。
所述云计算层包括误差补偿执行模块和热误差模型更新模块,所述误差补偿执行模块根据所述雾计算层输入的误差补偿值确定需要补偿的变量,并将补偿指令发送至霾计算层进行误差补偿;所述热误差模型更新模块根据边缘计算层输入的历史数据更新热误差预测模型,并将更新后的热误差预测模型传输至雾计算层。云计算层根据雾计算层提供的补偿值,确定需要补偿的变量,通过FB3输入补偿指令给西门子840D的霾计算层终端进行补偿。霾计算层接收来自云计算层的控制指令,通过伺服控制来控制各轴的实际位置,最后通过位置反馈调整实现重新控制。历史误差数据传输到配置有Hadoop的云服务器,再由云服务器用于重新训练ONT-GCN模型,ONT-GCN模型更新后的参数传递到雾层。另外,云计算层还提供监控任务,通过云服务器识别异常情况,实现各层的工作状态的监控,。
几何-热误差建模控制系统架构的合理布局,有利于提高补偿效率,缓解工业互联网的带宽压力。补偿建模架构如图1所示,用于构建霾-边-雾-云系统的功能划分见表1。
表1构建雾边雾云系统的功能划分
Figure BDA0003376641850000051
1、基于ONT-GCN时空模型的热误差预测模型及其建模方法的具体实施过程如下:
1.1、时空行为
将机床主轴简化为一维细长轴,如图2所示。主轴左端是自由的,由恒定的热流加热,右端固定。一维主轴的热微分方程为:
Figure BDA0003376641850000061
其中
Figure BDA0003376641850000062
其中c表示比热容,ρ表示密度,k表示热导率;
Figure BDA0003376641850000063
Figure BDA0003376641850000064
是拉普拉斯算子;l表示轴的长度;t表示时间;x表示距离;
热微分方程的初始条件如下:
T(x,0)=Tt
其中,T表示环境温度。为方便起见,坐标原点设置在自由端,因此边界条件为:
Figure BDA0003376641850000065
其中,h是一个正常数;温度T被离散化为:
T(x,t)=X(x)Time(t)
得到:
Figure BDA0003376641850000066
只有当上式等于一个常数时,才是合理的,即:
Figure BDA0003376641850000067
Figure BDA0003376641850000068
Figure BDA0003376641850000069
根据热微分方程的初始条件,Time=Tte-wt;根据X的特征方程,可以选择解为
Figure BDA00033766418500000610
根据边界条件A=0,并得到
Figure BDA00033766418500000611
为了避免平凡的解决方案,有一个方程
Figure BDA00033766418500000612
如果
Figure BDA00033766418500000613
Figure BDA00033766418500000614
该方程有无数个解,满足(k-1/2)π<vk<kπ。那么就可以得到
Figure BDA00033766418500000615
得到一个特殊的解为:
Figure BDA00033766418500000616
其中,Bk是常数。因此,热变形为:
Figure BDA0003376641850000071
表明热误差数据的时空行为对于指导经验模型的选择具有重要意义。为了进一步了解热误差的时间和空间特性,将热传导模型离散为n个节点,热传导的表达式如下:
Figure BDA0003376641850000072
其中T=[T1,T2,T3,…Tn]T;L=D-A,其中A表示邻接矩阵,D表示度矩阵;L表示拉普拉斯矩阵。GCN使用拉普拉斯矩阵来提取空间特征,GCN在提取特征过程中使用的Laplacian矩阵与热传导过程中使用的Laplacian矩阵是一致的,热扩散定律与GCN中的扩散定律是一致的。GCN可以描述空间热传导的特征,用于提取热误差的空间特征,再加上ON-LSTMN可以挖掘数据的时间特征,因此,所提出的ONT-GCN可以捕获热误差的时间和空间特征。
1.2、ON-LSTMN神经网络
LSTMN模型的结构如图3所示,遗忘门ft、输入门it、输出门ot如下图所示:
ft=σ(Wf[ht-1,xt]+bf)
Figure BDA0003376641850000073
Figure BDA0003376641850000074
其中,xt表示输入向量;ht-1表示前一个单元的输出值,与xt结合增加三个门的权重;
Figure BDA0003376641850000075
表示t时刻临时单元的内存状态;ct-1表示t-1时刻临时小区的小区状态;ct表示内存在t时刻的单元状态;Wf,Wi,Wc和Wo表示权重矩阵;bf,bi,bc和bo表示偏置向量;σ表示激活函数,是sigmoid;ht表示t时刻的隐藏状态。ON-LSTMN在LSTMN的基础上设计了主遗忘门
Figure BDA0003376641850000076
和主输入门
Figure BDA0003376641850000077
ft=δ(Wf[ht-1,xt]+bf)
it=δ(Wi[ht-1,xt]+bi)
Figure BDA0003376641850000078
ot=δ(Wo[ht-1,xt]+bo)
Figure BDA0003376641850000079
Figure BDA00033766418500000710
cumax(·)=cumsum(softmax(·))
Figure BDA00033766418500000711
Figure BDA00033766418500000712
Figure BDA00033766418500000713
Figure BDA00033766418500000714
ht=ot*tanh(ct)
其中,
Figure BDA0003376641850000081
Figure BDA0003376641850000082
分别是
Figure BDA0003376641850000083
Figure BDA0003376641850000084
的权重矩阵,
Figure BDA0003376641850000085
Figure BDA0003376641850000086
表示它们的偏置向量;wt1、wt2和wt3分别表示高等级、中等级和低等级,其中,wt1*ct-1存储低等级信息,受当前输入信息的影响;
Figure BDA0003376641850000087
存储高等级信息,受历史信息影响;
Figure BDA0003376641850000088
存储中等级信息,受当前输入信息和历史信息影响;comsum是一个求和的函数。ON-LSTMN的结构如图4所示。
ON-LSTMN对神经元进行排序后,信息水平由位置的前后表示。wt1*ct-1存储低等级信息,受当前输入信息的影响。
Figure BDA0003376641850000089
存储高等级信息,受历史信息影响。
Figure BDA00033766418500000810
存储的是中等级信息,受当前输入信息和历史信息影响。信息分级非常有利于热误差的预测。热误差数据包含长期趋势信息和短期瞬态信息。趋势信息是由于机器运行过程中温度持续升高引起的,瞬态信息是由环境干扰引起的,与当前时刻的输入密切相关。因此,需要合理保留热误差的长期和短期信息。ON-LSTMN可以将高等级信息保留相当远的距离,将热误差的趋势信息作为高等级信息。低等级信息可以在输入信息的每一步更新,热误差的瞬态信息作为低等级信息。重叠部分由
Figure BDA00033766418500000811
处理。ON-LSTMN捕获了热误差的趋势特征和电流输入特征,显着提高了基于ON-LSTMN的热误差预测精度。
1.3、ONT-GCN神经网络
为表征热误差的时间和空间行为,热误差的空间特征可以通过传感器网络的图G=(V,E)来反映,其中V是传感器的集合,E是传感器通道的集合。GCN的方程为:
Figure BDA00033766418500000812
其中,
Figure BDA00033766418500000813
I表示单位矩阵;A表示邻接矩阵;
Figure BDA00033766418500000814
X表示特征序列;W0表示第一层的权重参数;σ表示激活函数;
热误差时空预测模型需要时空特性,即(1)空间特性由传感器网络G反映;(2)位移传感器信号序列反映时间特性。因此,热误差的时空预测问题可以表示为传感器网络G与特征序列Xt之间映射关系f的学习,下一时刻T的热误差序列可以表示为:
[Xt+1,…,Xt+T]=f(G;(Xt-n,…,Xt-1,Xt))
其中,n表示历史时间序列的长度。
为了捕捉传感器网络的时间和空间依赖性,提出了基于ON-LSTMN和GCN的ONT-GCN模型。ONT-GCN的时空预测过程如图5所示。邻接矩阵和传感器的时间数据作为ONT-GCN的输入,GCN用于提取传感器的空间特征,ON-LSTMN模型用于提取序列的时间特征,如此即可实现数据预测了。
传输过程如图6所示,图6左边是热误差的时空预测,右边是ONT-GCN信元。GC是图卷积过程。ht-1和ht分别是t-1时刻和t时刻的输出。W和b分别是权重和偏差。将ON-LSTMN中的xt更新为f(X,A),则ft,it,
Figure BDA00033766418500000815
ot,
Figure BDA00033766418500000816
Figure BDA00033766418500000817
的六个门被改写为:
ft=δ(Wf[ht-1,f(A,Xt)]+bf)
it=δ(Wi[ht-1,f(A,Xt)]+bi)
Figure BDA0003376641850000091
ot=δ(Wo[ht-1,f(A,Xt)]+bo)
Figure BDA0003376641850000092
Figure BDA0003376641850000093
如此,所提出的ONT-GCN单元可以通过f(A,Xt)捕获传感器网络的空间依赖性,然后使用ft,it,
Figure BDA0003376641850000094
ot,
Figure BDA0003376641850000095
Figure BDA0003376641850000096
的六个独特门来有序地挖掘热误差的时间和空间依赖性,从而有效地完成了热误差的时空预测任务。
1.4、基于ONT-GCN时空模型的热误差预测模型及其建模方法
热误差时空预测框架如图14所示,本实施例基于ONT-GCN时空模型的热误差预测模型的建模方法,包括如下步骤:
11)构建传感器网络:在机床上布置温度传感器和位移传感器,并将各个传感器分别标记在机床的三维模型上。在蜗轮蜗杆机床上合理布置温度和位移传感器,采用9个测温点和1个位移测点,如图15所示。以上9个测温点可以完整反映整个机床的温度场分布并挖掘它们对位移传感器的依赖性。时空预测可以胜任挖掘相关性的任务,将传感器节点标记在蜗轮机床三维模型上,传感器网络如图16(a)所示,S代表位移传感器。由于X和Y方向的位移传感器相对较近,因此将它们视为同一节点。可以看出,S接近于T7、T8和T9,所以热误差主要受它们影响。
12)生成时空图数据集:根据时空图的假设,需要两个条件来预测当前节点的热误差信号,分别为:①位移测量点的历史热误差信号;②相邻传感器的位置和历史信号。图16(b)是反映距离的示意图,传感器网络的空间特征由邻接矩阵表征,导热通道在材料方向,不同部位的导热系数不同,因此,邻接矩阵应根据热导率进行加权。为简化模型,本实施例仅考虑材料方向的直线距离,计算邻接矩阵为:
Figure BDA0003376641850000097
其中,wi,j表示边权重,与节点i和节点j的位置有关;δ表示控制wi,j分布的阈值;di ,j为节点i和节点j之间的距离,是机床三维模型上测量得到的,如表5所示。
表5每两个节点之间的距离
Figure BDA0003376641850000098
从而得到邻接矩阵,如图17所示。可以看出,T8和T9比其他传感器对热误差的影响更大,因为它们更靠近位移传感器。以上9个测温点T8、T9温升明显,说明T8、T9是引起热误差的主要热源。
图18为转速800r/min,运行时间100分钟左右得到的图;图19是在1500r/min的速度和大约45分钟的运行时间下得到的。这两个实验的结果用于获得热特性的时序特性。图18用作训练集,图19用作测试集。速度和运行时间的不同导致温度场和热变形的不同趋势。但是,蜗轮蜗杆机床的热误差产生机制没有改变,即热误差的时空依赖性不会改变。时空预测模型用于学习这种依赖性,变化的条件可以测试所提出的ONT-GCN模型的鲁棒性。
13)构建如上所述基于ONT-GCN时空模型的热误差预测模型。收敛时间与硬件和软件平台有关,详情如下:使用Windows 10操作系统,CPU为第11代Intel Core i5,RAM为16.0GB。该程序由Python编程语言开发,开发环境为PyCharm Community Edition2020.2.5,管理使用Anaconda。为了便于机器学习建模,使用Pytorch机器学习库进行编程。
为了测试所提出的ONT-GCN模型的有效性,使用MLR、标准LSTMN和ON-LSTMN与ONT-GCN进行比较。超参数是影响ONT-GCN模型预测性能和鲁棒性的重要因素,通过贝叶斯优化算法进行优化。由于传感器节点的信号电平不同,在训练模型之前需要对特征序列进行归一化,然后对预测结果进行反向归一化。为了比较,LSTMN、ON-LSTMN和ONT-GCN被设置为具有表6中列出的相同超参数。
表6模型超参数
Figure BDA0003376641850000101
此外,MLR作为传统的热误差建模方法也用于与ONT-GCN进行比较。通过模糊聚类选择了三个温度变量T5、T8和T9,调用MATLAB中的regress toolbox,确定X方向热误差的回归系数为b0=0.9974,b1=84.677,b2=-36.593and b3=28.687;b0=-0.3514,b1=133.667,b2=-35.048,and b3=-64.987为Y方向。那么,建立的MLR模型为:
Figure BDA0003376641850000102
结果表明,MLR只将温度作为输入,并没有考虑两个重要的关键点,即(1)热误差序列的时间特性,以及(2)传感器节点的空间特性。因此,MLR无法表征错误的时间和空间特征。
14)以邻接矩阵和传感器的时间数据作为热误差预测模型的输入,预测热误差。
X-和Y-方向的预测结果如图20所示。MLR模型只能捕捉到热误差的趋势,而不能捕捉到热误差的瞬态特征。MLR受温度影响很大,未经过滤的温度波会在热误差的预测结果中产生噪声。时间序列预测模型可以更好地捕捉热误差数据的非线性特征,而且,具有空间特征提取能力的时间序列模型比其他误差预测模型具有更优异的预测性能。
然后计算评价指标,如表7至8所示。ONT-GCN的RMSE小于ON-LSTMN模型的RMSE,因为ONT-GCN模型基于热误差产生机制考虑了温度传感器的空间特性。ONT-GCN模型比ON-LSTMN模型消耗更多的时间,因为卷积操作增加了计算量。ON-LSTMN模型的RMSE小于LSTMN模型,因为ON-LSTMN模型保留了基于有序神经元的热误差的长期趋势和短期瞬态特性。ON-LSTMN模型比LSTMN模型花费更多的时间,因为额外的
Figure BDA0003376641850000103
Figure BDA0003376641850000104
增加了计算量。LSTMN模型的RMSE小于MLR模型,因为LSTMN模型可以反映热误差的时间特征。LSTMN模型比MLR需要更多的时间,因为LSTMN需要根据输入不断改变门的状态,而MLR模型是固定参数的模型。因此,热误差的时间特性是首要考虑的问题。引入温度传感器的空间特性可以进一步探索温度场之间的关系,提高热误差的预测精度。
表7X方向热误差预测性能评价
Figure BDA0003376641850000111
表8Y方向热误差预测性能评价
Figure BDA0003376641850000112
2、误差映射模型
本实施例的误差映射模型为根据蜗轮蜗杆机床的传递参数建立的齿面误差与几何-热误差的映射模型。
蜗轮蜗杆机床的结构和传动顺序分别如图7(a)和图7(b)所示,传动链0-3-4-5-6-7-8从床身传送到滚刀端,传动链0-1-2从床身传送到工件端。蜗轮蜗杆机床有六个主要的运动轴,功能见表2。X、Y、Z、A、B、C轴的误差是影响机床加工精度的主要误差源,因此,有必要重点分析这些误差。
表2机床运动轴功能特性
Figure BDA0003376641850000113
齿轮加工时,滚刀与工件的位置关系如图8所示,根据齿轮啮合原理,滚刀(B轴)的旋转角度B和工件(C轴)的旋转角度
Figure BDA0003376641850000114
有固定的比例关系。蜗轮和滚刀的模数为mn,螺旋角为β和λ,齿数分别为z1和z2,工件和滚刀在P点具有相同的切向速度V0,即:
Figure BDA0003376641850000115
其中,ωg和ωh分别表示工件和滚刀的角速度。那么工件转角
Figure BDA0003376641850000116
与滚刀转角B的比例关系表示为:
Figure BDA0003376641850000117
具体的,齿面误差与几何-热误差的映射模型的创建方法为:
21)根据蜗轮蜗杆机床的坐标系传动链,得到理想状态下滚刀坐标系H相对于工件坐标系G的齐次变换矩阵
Figure BDA0003376641850000118
在理想状态下,蜗轮蜗杆机床各运动副均无运动误差,每两个相邻运动对之间的齐次坐标变换矩阵(齐次坐标变换矩阵)表示为:
Figure BDA0003376641850000121
Figure BDA0003376641850000122
理想状态下滚刀坐标系H相对于工件G坐标系的齐次坐标变换矩阵
Figure BDA0003376641850000123
为:
Figure BDA0003376641850000124
其中,
Figure BDA0003376641850000125
Figure BDA0003376641850000126
分别表示理想状态下C轴相对于工件坐标系的齐次坐标变换矩阵和滚刀坐标系相对于B轴的齐次坐标变换矩阵,且
Figure BDA0003376641850000127
Figure BDA0003376641850000128
表示表示理想状态下床身R相对于C轴坐标系的齐次坐标变换矩阵,
Figure BDA0003376641850000129
Figure BDA00033766418500001210
表示表示理想状态下C轴相对于床身R坐标系的矩阵;
Figure BDA00033766418500001211
表示表示理想状态下X轴相对于床身R坐标系的距离为x的矩阵;
Figure BDA00033766418500001212
表示表示理想状态下Z轴相对于X轴的距离为z的矩阵;
Figure BDA00033766418500001213
表示表示理想状态下A轴相对于Z轴的矩阵;
Figure BDA00033766418500001214
表示表示理想状态下Y轴相对于A轴的矩阵;
Figure BDA00033766418500001215
表示表示理想状态下B轴相对于Y轴的矩阵;
22)在实际状态下,结合理想状态下滚刀坐标系相对于工件坐标系的齐次变换矩阵
Figure BDA00033766418500001216
考虑各运动副的误差,得到滚刀坐标系H相对于工件坐标系G的齐次变换矩阵
Figure BDA00033766418500001217
在实际加工过程中,制造、装配、调试、运动等过程中都会产生各种误差,实际状态下的齐次坐标变换矩阵与理想状态下的齐次坐标变换矩阵存在差异。实际状态下
Figure BDA00033766418500001218
的齐次坐标变换矩阵为:
Figure BDA00033766418500001219
其中,
Figure BDA00033766418500001220
Figure BDA00033766418500001221
分别表示P坐标系到Q坐标系的运动误差矩阵、几何误差矩阵和热误差矩阵。
对于X轴相对于床身的齐次坐标变换矩阵,几何误差矩阵为
Figure BDA00033766418500001222
运动误差矩阵
Figure BDA00033766418500001223
表示为:
Figure BDA00033766418500001224
其中,εxx、εyx和εzx分别表示x轴沿x方向的定位误差,x轴沿y方向的直线度误差,x轴沿z方向的直线度误差;εxx、ωyx和εzx分别表示x轴绕x方向的角度误差,x轴绕y方向的角度误差和x轴沿z方向的角度误差。
而且,X轴是机床的最大部件,需要考虑其热误差。床身与工作台刚性连接,因此X轴和C轴一起分析。床身的加热主要来自两个方面:一是加工过程中冷却液与金属屑形成的高温混合物所产生的热量,二是立柱在X轴方向移动时因摩擦产生的热量。床身向上弯曲变形引起X轴上的两个位置误差,分别是X方向的位置误差
Figure BDA00033766418500001225
和Z方向的位置误差
Figure BDA00033766418500001226
如图9所示。
Figure BDA00033766418500001227
Figure BDA0003376641850000131
其中,L表示床身上表面长度;H表示立柱内侧面长度;β表示床身热变形引起X轴在Y方向的角度误差;ΔL表示床身上表面在X轴方向总的受热膨胀量。
那么热误差矩阵为:
Figure BDA0003376641850000132
则X轴相对于床身坐标系的齐次坐标变换矩阵
Figure BDA0003376641850000133
为:
Figure BDA0003376641850000134
对于Z轴相对于X轴的齐次坐标变换矩阵,几何误差矩阵
Figure BDA0003376641850000135
和运动误差矩阵
Figure BDA0003376641850000136
表示为:
Figure BDA0003376641850000137
其中,Sxz表示z轴沿x方向的垂直度误差;δxz表示z轴沿x方向的直线度误差;δyz表示z轴沿y方向的直线度误差;δzz表示z轴沿z方向的定位误差;εxz表示z轴绕x方向的角度误差;εyz表示z轴绕y方向的角度误差;εzz表示z轴绕z方向的角度误差。
此外,Z轴受热误差影响很大,因为立柱是蜗轮机床的第二大部件。立柱向外弯曲变形,导致Z轴产生两个位置误差,分别是X方向的位置误差
Figure BDA0003376641850000138
和Z方向的位置误差
Figure BDA0003376641850000139
如图10所示。
Figure BDA00033766418500001310
Figure BDA00033766418500001311
其中,H'表示立柱上表面长度;W表示立柱外侧面长度;α表示立柱热变形引起Z轴在Y方向的角度误差;ΔH'表示表示立柱在Z轴方向总的受热膨胀量。
热误差矩阵为:
Figure BDA00033766418500001312
可以得到Z轴相对于X轴的齐次坐标变换矩阵
Figure BDA00033766418500001313
Figure BDA0003376641850000141
对于A轴相对于Z轴的齐次坐标变换矩阵,几何误差矩阵
Figure BDA0003376641850000142
和运动误差矩阵
Figure BDA0003376641850000143
表示为:
Figure BDA0003376641850000144
其中,δ'ya表示a轴沿y方向的位置误差;δ'za表示a轴沿z方向的位置误差;ε'ya表示a轴在xz平面内的平行度误差;ε'za表示a轴在xy平面内的平行度误差;δxa表示a轴沿x方向的直线度误差;δya表示a轴沿y方向的直线度误差;δza表示a轴沿z方向的直线度误差;εxa表示a轴绕x方向的角度误差;εya表示a轴绕y方向的角度误差;εza表示a轴绕z方向的角度误差;γ表示a轴绕x方向旋转的角度。
则A轴相对于Z轴的齐次坐标变换矩阵
Figure BDA0003376641850000145
为:
Figure BDA0003376641850000146
对于Y轴相对于A轴的齐次坐标变换矩阵,几何误差矩阵
Figure BDA0003376641850000147
和运动误差矩阵
Figure BDA0003376641850000148
表示为:
Figure BDA0003376641850000149
其中,Sxy表示y轴沿x方向的垂直度误差;Szy表示y轴沿z方向的垂直度误差;δxy表示y轴沿x方向的直线度误差;δyy表示y轴沿y方向的定位误差;δzy表示y轴沿z方向的直线度误差;εxy表示y轴绕x方向的角度误差;εyy表示y轴绕y方向的角度误差;εzy表示y轴绕y方向的角度误差。
则Y轴相对于A轴的齐次坐标变换矩阵
Figure BDA00033766418500001410
为:
Figure BDA00033766418500001411
对于B轴相对于Y轴的齐次坐标变换矩阵,几何误差矩阵
Figure BDA00033766418500001412
和运动误差矩阵
Figure BDA00033766418500001413
表示为:
Figure BDA00033766418500001414
其中,δxb表示b轴沿x方向的直线度误差;δyb表示b轴沿y方向的直线度误差;δzb表示b轴沿z方向的直线度误差;εxb表示b轴绕x方向的角度误差;εyb表示b轴绕y方向的角度误差;εzb表示b轴绕z方向的角度误差;δ'xb表示b轴沿x方向的位置误差;δ'zb表示b轴沿z方向的位置误差;ε'xb表示b轴在yz平面内的平行度误差;ε'zb表示b轴在xy平面内的平行度误差。
则B轴到Y轴的齐次坐标变换矩阵
Figure BDA0003376641850000151
为:
Figure BDA0003376641850000152
对于床身C轴相对于R坐标系的齐次坐标变换矩阵,几何误差矩阵
Figure BDA0003376641850000153
和运动误差矩阵
Figure BDA0003376641850000154
表示为:
Figure BDA0003376641850000155
其中,δ'xc表示c轴沿x方向的位置误差;δ'yc表示表示c轴沿y方向的位置误差;ε'xc表示c轴在yz平面内的平行度误差;ε'yc表示c轴在xz平面内的平行度误差;δxc表示c轴沿x方向的直线度误差;δyc表示c轴沿y方向的直线度误差;δzc表示c轴沿z方向的直线度误差;εxc表示c轴绕x方向的角度误差;εyc表示c轴绕y方向的角度误差;εzc表示c轴绕z方向的角度误差。
则C轴到床身的齐次坐标变换矩阵
Figure BDA0003376641850000156
为:
Figure BDA0003376641850000157
在实际状态下,考虑各运动副的误差,滚刀坐标系H相对于工件坐标系G的齐次坐标变换矩阵
Figure BDA0003376641850000158
为:
Figure BDA0003376641850000159
其中,
Figure BDA00033766418500001510
Figure BDA00033766418500001511
分别表示实际状态下C轴相对于工件坐标系的齐次坐标变换矩阵和滚刀坐标系相对于B轴的齐次坐标变换矩阵,且
Figure BDA00033766418500001512
Figure BDA00033766418500001513
表示实际状态下X轴相对于床身坐标系的齐次坐标变换矩阵;
Figure BDA00033766418500001514
表示实际状态下Z轴相对于X轴的齐次坐标变换矩阵;
Figure BDA00033766418500001515
表示实际状态下A轴相对于Z轴的齐次坐标变换矩阵;
Figure BDA00033766418500001516
表示实际状态下Y轴相对于A轴的齐次坐标变换矩阵;
Figure BDA00033766418500001517
表示实际状态下B轴相对于Y轴的齐次坐标变换矩阵;
Figure BDA00033766418500001518
表示实际状态下C轴相对于床身R坐标系的齐次坐标变换矩阵,且
Figure BDA00033766418500001519
将热误差和几何误差映射到工具的空间误差为:
Figure BDA00033766418500001520
其中,[a,b,c,1]T和[i,j,k,0]T分别表示滚刀中心位置和轴姿态的齐次坐标;δ和ε分别表示位置误差和姿态误差,δx、δy和δz分别表示x、y、z方向的位置误差分量;εx、εy和εz分别表示x、y、z方向的姿态误差分量。
23)创建齿面误差与几何-热误差的映射模型
渐开线齿轮滚刀端面齿廓为标准渐开线,如图11所示。齿廓基圆半径rbh对称,σ0为基齿槽半角,u为渐开线展开角。那么坐标系OPXPYP中渐开线a0b0的参数方程为:
Figure BDA0003376641850000161
其中,
Figure BDA00033766418500001618
表示滚刀的转角。
顺时针旋转时,得到角度θ,得到渐开线齿面a1b1,参数方程为:
Figure BDA0003376641850000162
Figure BDA0003376641850000163
其中,Ph表示滚刀的螺旋参数,Ph=mnz/[2·sin(π/2-λ)];mn和z分别表示齿轮滚刀的法向模数和头数。齿面的法向量为:
Figure BDA0003376641850000164
在理想状态下,理想齿面在工件坐标系中的坐标向量和单位法向量表示为:
Figure BDA0003376641850000165
Figure BDA0003376641850000166
其中,
Figure BDA0003376641850000167
表示理想状态下工件坐标系G相对于滚刀坐标系H的齐次变换矩阵,且
Figure BDA0003376641850000168
Figure BDA0003376641850000169
表示理想齿面在工件坐标系中的坐标向量;rh表示滚刀齿面的渐开线参数方程;
Figure BDA00033766418500001610
表示理想齿面在工件坐标系中的单位法向量;Nh表示滚刀齿面的单位法向量;
则在考虑了热误差和几何误差的实际状态下,工件坐标系中的坐标向量和单位法向量表示为:
Figure BDA00033766418500001611
Figure BDA00033766418500001612
其中,
Figure BDA00033766418500001613
表示实际状态下工件坐标系G相对于滚刀坐标系H的齐次变换矩阵,且
Figure BDA00033766418500001614
Figure BDA00033766418500001615
表示实际齿面在工件坐标系中的坐标向量;
Figure BDA00033766418500001616
表示实际齿面在工件坐标系中的单位法向量;
从而得到齿面误差与几何-热误差的映射模型为:
Figure BDA00033766418500001617
其中,n表示离散值;齿面误差定义为E=[δ;ε]=[δxyzxyz]T,其中δxy和δz分别表示x、y、z方向的齿面位置误差分量;εxyz分别表示齿面单位法向量在x、y、z方向的分量。
进一步,基于所述齿面误差与几何-热误差的映射模型,采用Morris方法识别出齿面误差的关键误差项,所述雾计算层中,结合所述齿面误差与几何-热误差的映射模型得到所述关键误差项的误差补偿值,并将关键误差项的误差补偿值传输至云计算层。
具体的,对蜗轮蜗杆加工误差敏感性分析的方法如下。
敏感性分析方法常用于统计学或经济学领域,该方法的目标是识别和分析对目标结果有重要影响的敏感因素。基于齿面误差与几何-热误差的映射模型,采用Morris方法识别出齿面误差的关键误差项和敏感部位,为后续研究提供理论支持。敏感性分析流程如图12所示。
循环次数N极大地影响Morris方法的运行效率。经过多次测试,选择N=50来平衡实验效率和求解精度。选取直线轴和旋转轴的几何误差和热误差G=[x1,x2,...,xm]T作为模型输入参数,以6个齿面误差E=[δxyzxyz]T作为分析目标。输出是每个几何误差和热误差的均值μi和标准差σi。具体的,以δx为分析对象,将Morris方法应用于机床的几何误差。齿面误差映射模型的具体过程如下:
(1)求解齐次变换矩阵、渐开线滚刀模型、几何和热误差G=[x1,x2,...,x55]T
(2)利用机床的齐次坐标变换矩阵将渐开线滚刀方程转换为工件坐标系,计算渐开线扩展角u的取值范围,并根据扩展角取值范围将转换后的方程离散化,获得一定数量的离散点坐标方程并计算法向量方程Ng
(3)测量每个运动轴的热和几何误差分布。测量结果表明,热误差、几何位移误差和几何角度误差分别在[0,100μm]、[0,20μm]和[0,29mdeg]范围内均匀分布。热误差、几何位移误差、几何角度误差的误差采样区间分别生成为{0,Δ1,2Δ1,...,100},{0,Δ1,2Δ1,...,20}和{0,Δ3,2Δ3,...,29}。其中Δ1=1/(q-1)×100,Δ2=1/(q-1)×20和Δ3=1/(q-1)×29,生成几何误差基向量
Figure BDA0003376641850000171
向量元素的值在采样间隔内随机采样。
(4)误差个数为m=55,构造并生成随机化矩阵
Figure BDA0003376641850000172
和随机化采样矩阵
Figure BDA0003376641850000173
然后将
Figure BDA0003376641850000174
的每一行元素作为齿面误差与几何-热误差,并且任何相邻的两行都只有某个误差元素与Δ不同。
(5)将步骤(4)中得到的随机抽样矩阵
Figure BDA0003376641850000175
的第ith行代入接触条件作为输入参数,求解二元非线性方程组。得到滚刀加工面与待加工蜗轮齿面的接触点组。在输入误差参数的条件下,这些接触点形成的轨迹被用作加工的齿面方程。
(6)计算每个误差项的基本效果,重复步骤(4)至(5)。独立计算每个误差的N个基本效应,得到其均值μi和标准差σi
本实施例的滚刀和蜗轮蜗杆加工参数设置见表3。
表3滚刀、蜗轮加工参数设置
Figure BDA0003376641850000176
Figure BDA0003376641850000181
误差敏感性分析结果如图13所示。横坐标为由51个几何误差和4个热误差项组成的误差序列,纵坐标为均值μi和标准偏差σi。比较误差序列中每个误差项的灵敏度指标μi和σi的大小,根据具体的齿面成分识别出关键误差项。也可以通过对运动轴的整体评价来判断某一运动轴误差项的敏感性,以齿面误差δx为例进行分析。根据图13,关键误差项通过比较均值μi确定,分别为δxxyxxyyyxzyzxayaya',εybxcycyc',δxc',δt xx,andδt zx,,通过比较标准差σi确定耦合误差,分别为εyxzxzy,Sxyyzza,,εya',εzbzb',εyczcyc',δt zzt xzt xx,andδt zx
综合评价表明,敏感元件的顺序为X>C>A>Z。因此,应综合考虑上述误差项,调整C轴和X轴的精度,以减小齿面误差。可以得到每个齿面误差的关键误差和敏感部位,如表4所示。因此,找出关键误差并根据计算出的补偿值进行补偿,可以提高蜗轮蜗杆加工机床的加工精度。
表4关键错误和敏感部分
Figure BDA0003376641850000182
3、系统构架验证
为了测试所提出架构的有效性,表9列出了以ONT-GCN模型作为错误预测模型的不同系统架构的执行时间。霾-云、霾-边缘-云、霾-雾-云和本实施例提出的霾-边-雾-云架构系统的总执行时间分别为184s、179s、164s和153s。霾-边-雾-云系统架构为四层架构。综合发挥霾计算层、雾计算层和边缘计算层服务器的优势。对于霾-边-雾-云系统,在边缘层进行数据异常检查。将训练有素的ONT-GCN模型嵌入到雾层中以实现热误差预测。雾层中嵌入了齿面误差与几何热误差的映射模型和灵敏度分析模型。云层用于实现ONT-GCN误差模型的再训练。ONT-GCN误差模型更新后的参数传递到雾层,实现训练好的ONT-GCN模型的参数更新。霾-云系统架构只有两层。为保证比较的公平性,将数据异常检查和模型预测再训练放在云层,导致云层计算量增加;此外,由于缺少边缘层和雾层,云层的数据传输时间和时间消耗显着增加,导致霾-云架构最耗时。霾-边-云架构分为三层,边缘层分担霾-边-雾-云系统架构中雾层的工作,承担数据异常检查任务,导致边缘层花费的时间增加,云层花费的时间比霾-云架构的云层花费的时间减少;边缘层的边缘服务器具有强大的计算能力和数据传输能力,导致霾-边-云架构的执行时间比霾-云架构短。霾-雾-云架构也有三层,雾层分担边缘层的任务,承担数据异常检查任务;雾服务器的计算能力和数据传输能力强于边缘服务器,因此,雾层可以承担数据异常检查和错误预测等任务,所以霾-雾-云架构的效率要高于雾-边-云架构。对于所提出的霾-边-雾-云系统架构,与霾-雾-云架构相比,数据传输时间减少,并且所提出的系统架构在所有架构中效率最高。
表9ONT-GCN不同框架的时间消耗
Figure BDA0003376641850000183
Figure BDA0003376641850000191
3.1、加工实验验证
随着提出的霾-边-雾-云系统的实施,进行了蜗轮的加工实验,如图21所示。求解了网格点坐标和单位法向量。蜗轮机参数如下:加工的蜗轮为右手蜗轮,蜗杆头数z1为1,蜗轮齿数z2为90,齿廓角a为20°,模数m为5,分度圆的直径d1为320mm,修正系数x2为0.042,齿顶高度系数h为1,头部间隙系数为0.2。
如图22所示,可以清楚地获得整个齿面的误差分布,得到蜗轮齿面误差测量的最终结果。根据网格平面图,齿面有5列9行网格点,每条网格有45个网格点,如图22所示。红色实线代表实际齿面点数据,蓝色交错虚线表示理论齿面点数据,上下两侧为齿面齿根,中间部分为齿面齿顶。两个网格用于表示整个牙齿的左侧和右侧的误差。每个测量点的误差值用该点垂直黑色实线的长度表示。箭头符号表示左右两侧齿面误差的方向。
给定相同的实验条件,多次测量同一被测齿面,并通过多次测量计算平均误差值。
Figure BDA0003376641850000192
其中,xi表示第i次测量得到的误差;n为误差测量的次数。
每个测量网格点的单个数据的标准不确定度由贝塞尔方程计算。
Figure BDA0003376641850000193
标准不确定度为:
Figure BDA0003376641850000194
由测量重复性引起的标准不确定度ux的计算属于A类不确定度,即
Figure BDA0003376641850000195
自由度为v=n-1。
计算正态误差,如表10中所列。
表10齿面法向误差
Figure BDA0003376641850000196
Figure BDA0003376641850000201
如图23所示,得到齿面误差。不实施霾-边-雾-云系统时,齿面第一、二、三、四、五行几何误差在范围内分别为[5.75μm,23μm]、[5.75μm,21.85μm]、[0,23μm]、[-0.23μm,17.25μm]和[5.75μm,19.55μm]。采用霾-边-雾-云系统进行几何误差控制时,齿面第一、二、三、四、五行几何误差在[2.95μm,13.47μm]范围内,[5.26μm,13.31μm],[0.12,14.13μm],[-0.28μm,11.27μm],和[4.60μm,12.81μm]。采用霾-边-雾-云系统进行几何-热误差控制时,齿面第一、二、三、四、五行几何误差在[3.27μm,12.87μm],[2.35μm,11.70μm],[0.23,13.28μm],[-0.20μm,8.90μm]和[2.23μm,10.35μm]范围内。结论是,霾-边-雾-云系统的实施减少了齿面几何误差,因此,所提出的系统架构的效率得到验证。此外,嵌入模型,包括ONT-GCN模型、渐开线滚刀模型以及齿面误差与几何热误差之间的映射模型,都得到了验证。更重要的是,采用雾边-雾云系统进行几何误差控制的齿面几何误差大于采用雾边-雾云系统进行几何-热误差控制的齿面几何误差。控制。因此,证明了几何热误差控制的必要性。
4、结论
霾-边-雾-云计算系统基于ONT-GCN时空模型进行几何和热误差补偿与预测。基于ONT-GCN时空模型提出深度学习热误差模型。此外,提出了齿面误差与几何-热误差的映射模型和灵敏度分析模型,并将其嵌入到霾-边-雾-云计算系统的雾层中。证明了热补偿和几何误差补偿的必要性。最后,通过执行所提出的霾-边-雾-云计算系统,实现了同时的几何-热误差控制。主要结论如下:
(1)在误差控制方面,设计了霾-边-雾-云计算系统来补偿几何误差和热误差。霾-云系统、雾-边-云系统、霾-雾-云系统和霾-边-雾-云系统的总执行时间分别为184s、179s、164s和153s。霾-边-雾-云计算系统实现了合理的任务分配,所提出的霾-边-雾-云系统的执行效率高于霾-云系统、霾-边-云系统和霾-雾-云系统。
(2)求解一维热传导的偏微分方程,揭示误差机制。结果表明,热变形是运行时间和空间的函数。应建立一个能够捕捉时空特征的模型。错误机制为ONT-GCN模型提供了理论依据。为提高热误差的预测精度,构建了蜗轮蜗杆机床传感器网络,反映热误差与温度的空间位置关系,并提出了用于时空预测的ONT-GCN模型。一个传感器网络可以充分反映蜗轮蜗杆机床的热状态,网络拓扑作为ONT-GCN的输入之一。提议的ONT-GCN模型可以保留和捕获ONT-GCN单元中的空间和时序特征。ONT-GCN单元的独特排序特性允许保留重要的热误差信息。因此,ONT-GCN单元可以有序地传递热误差的时间和空间特征,从而比其他机器学习模型更高的预测精度。MLR、LSTMN、ON-LSTMN和ONT-GCNN模型的预测准确率分别为93.18%、94.56%、95.36%和97.35%。MLR、LSTMN、ON-LSTMN和ONT-GCN的RMSE分别为0.7642、0.6931、0.4822和0.4167。
(3)考虑热误差产生的机理,考虑X轴和Z轴的剧烈热变形,根据几何关系反映热误差的变化。然后获得几何误差和热误差的转换矩阵。齿面误差与几何热误差的映射模型及灵敏度分析模型。使用Morris方法识别关键误差项,然后通过映射模型对关键误差项进行补偿。当采用霾-边-雾-云系统进行几何误差控制时,齿面第一、二、三、四、五排的几何误差减小到[2.95μm,13.47μm],[5.26μm,13.31μm],[0.12,14.13μm],[-0.28μm,11.27μm]和[4.60μm,12.81μm]范围内。当采用霾-边-雾-云系统进行几何-热误差控制时,齿面第一、二、三、四、五排的几何误差减小到[3.27μm,12.87μm]、[2.35μm,11.70μm]、[0.23,13.28μm]、[-0.20μm,8.90μm]和[2.23μm,10.35μm]范围内。
以上所述实施例仅是为充分说明本发明而所举的较佳的实施例,本发明的保护范围不限于此。本技术领域的技术人员在本发明基础上所作的等同替代或变换,均在本发明的保护范围之内。本发明的保护范围以权利要求书为准。

Claims (6)

1.一种基于ONT-GCN时空模型的热误差预测模型,其特征在于:所述ONT-GCN神经网络的表达式为:
ft=δ(Wf[ht-1,f(A,Xt)]+bf)
it=δ(Wi[ht-1,f(A,Xt)]+bi)
Figure FDA0003376641840000011
ot=δ(Wo[ht-1,f(A,Xt)]+bo)
Figure FDA0003376641840000012
Figure FDA0003376641840000013
cumax(·)=cumsum(softmax(·))
Figure FDA0003376641840000014
Figure FDA0003376641840000015
Figure FDA0003376641840000016
Figure FDA0003376641840000017
ht=ot*tanh(ct)
其中,ft表示遗忘门;it表示输入门;ot表示输出门;
Figure FDA0003376641840000018
表示主遗忘门;
Figure FDA0003376641840000019
表示主输入门;ht-1表示前一个单元的输出值,与f(A,Xt)结合增加三个门的权重,Xt表示特征序列,A表示邻接矩阵,f(A,Xt)表示Xt经过图卷积之后的序列;
Figure FDA00033766418400000110
表示t时刻临时单元的内存状态;ct-1表示t-1时刻临时单元的内存状态;ct表示内存在t时刻的单元状态;Wf、Wi、Wc和Wo表示权重矩阵;bf、bi、bc和bo表示偏置向量;σ表示激活函数;ht表示t时刻的隐藏状态;
Figure FDA00033766418400000111
Figure FDA00033766418400000112
分别表示
Figure FDA00033766418400000113
Figure FDA00033766418400000114
的权重矩阵;
Figure FDA00033766418400000115
Figure FDA00033766418400000116
分别表示
Figure FDA00033766418400000117
Figure FDA00033766418400000118
的偏置向量;comsum表示求和的函数;wt1、wt2和wt3分别表示高等级、中等级和低等级,其中,wt1*ct-1存储低等级信息,受当前输入信息的影响;
Figure FDA00033766418400000119
存储高等级信息,受历史信息影响;
Figure FDA00033766418400000120
存储中等级信息,受当前输入信息和历史信息影响;
f(A,Xt)表示Xt经过图卷积之后的序列,并作为ON-LSTMN模型的输入向量,且:
Figure FDA00033766418400000121
其中,
Figure DEST_PATH_FDA00035244144100000120
I表示单位矩阵;A表示邻接矩阵;
Figure DEST_PATH_FDA00035244144100000121
X表示特征序列;W0表示第一层的权重参数;σ表示激活函数;
热误差预测模型需要时空特性,因此,热误差的时空预测可以表示为传感器网络G与特征序列Xt之间映射关系f的学习,下一时刻T的热误差序列可以表示为:
[Xt+1,,Xt+T]=f(G;(Xt-n,,Xt-1,Xt))
其中,n表示历史时间序列的长度。
2.一种基于ONT-GCN时空模型的热误差预测模型的建模方法,其特征在于:包括如下步骤:
11)构建传感器网络:在机床上布置温度传感器和位移传感器,并将各个传感器分别标记在机床的三维模型上;
12)生成时空图数据集:需要两个条件来预测当前节点的热误差信号,分别为:①位移测量点的历史热误差信号;②相邻传感器的位置和历史信号;
传感器网络的空间特征由邻接矩阵表征,邻接矩阵为:
Figure FDA0003376641840000021
其中,wi,j表示边权重,与节点i和节点j的位置有关;δ表示控制wi,j分布的阈值;di,j为节点i和节点j之间的距离,是机床三维模型上测量得到的;
13)构建如权利要求1所述基于ONT-GCN时空模型的热误差预测模型;
14)以邻接矩阵和传感器的时间数据作为热误差预测模型的输入,预测热误差。
3.一种霾-边-雾-云误差补偿系统,其特征在于:包括霾计算层、边缘计算层、雾计算层和云计算层;
所述霾计算层包括:
温度传感器:用于测量温度;
位移传感器:用于测量位移;
采集卡:用于采集所述温度传感器和位移传感器测量得到的温度信号和位移信号;
树莓派:所述树莓派上连接有用于将温度信号和位移信号转换为数字信息的A/D转换器;
网关:用于经A/D转换器转换得到的数据传输至所述边缘计算层;
所述边缘计算层包括用于存储经所述霾计算层传入的数据的微型数据中心、用于滤波数据的交换机和用于将滤波后的数据传输至雾计算层和云计算层的网关和路由器,其中,历史数据传输至云计算层,实时误差数据传输至雾计算层;
所述雾计算层内设有分析模块,所述分析模块内设有如权利要求1所述基于ONT-GCN时空模型的热误差预测模型,将以温度传感器、位移传感器的位置信息生成的邻接矩阵和以各个传感器采集的数据作为特征序列作为热误差预测模型的输入预测热误差,再结合误差映射模型得到误差补偿值,并将误差补偿值传输至云计算层;
所述云计算层包括误差补偿执行模块和热误差模型更新模块,所述误差补偿执行模块根据所述雾计算层输入的误差补偿值确定需要补偿的变量,并将补偿指令发送至霾计算层进行误差补偿;所述热误差模型更新模块根据边缘计算层输入的历史数据更新热误差预测模型,并将更新后的热误差预测模型传输至雾计算层。
4.根据权利要求3所述的霾-边-雾-云误差补偿系统,其特征在于:所述误差映射模型为根据蜗轮蜗杆机床的传递参数建立的齿面误差与几何-热误差的映射模型。
5.根据权利要求4所述的霾-边-雾-云误差补偿系统,其特征在于:所述齿面误差与几何-热误差的映射模型的创建方法为:
21)根据蜗轮蜗杆机床的坐标系传动链,得到理想状态下滚刀坐标系H相对于工件坐标系G的齐次变换矩阵
Figure FDA0003376641840000031
为:
Figure FDA0003376641840000032
其中,
Figure FDA0003376641840000033
Figure FDA0003376641840000034
分别表示理想状态下C轴相对于工件坐标系的齐次坐标变换矩阵和滚刀坐标系相对于B轴的齐次坐标变换矩阵,且
Figure FDA0003376641840000035
Figure FDA0003376641840000036
表示表示理想状态下床身R相对于C轴坐标系的齐次坐标变换矩阵,
Figure FDA0003376641840000037
Figure FDA0003376641840000038
表示表示理想状态下C轴相对于床身R坐标系的矩阵;
Figure FDA0003376641840000039
表示表示理想状态下X轴相对于床身R坐标系的距离为x的矩阵;
Figure FDA00033766418400000310
表示表示理想状态下Z轴相对于X轴的距离为z的矩阵;
Figure FDA00033766418400000311
表示表示理想状态下A轴相对于Z轴的矩阵;
Figure FDA00033766418400000312
表示表示理想状态下Y轴相对于A轴的矩阵;
Figure FDA00033766418400000313
表示表示理想状态下B轴相对于Y轴的矩阵;
22)在实际状态下,结合理想状态下滚刀坐标系相对于工件坐标系的齐次变换矩阵
Figure FDA00033766418400000314
考虑各运动副的误差,得到滚刀坐标系H相对于工件坐标系G的齐次变换矩阵
Figure FDA00033766418400000315
为:
Figure FDA00033766418400000316
其中,
Figure FDA00033766418400000317
Figure FDA00033766418400000318
分别表示实际状态下C轴相对于工件坐标系的齐次坐标变换矩阵和滚刀坐标系相对于B轴的齐次坐标变换矩阵,且
Figure FDA00033766418400000319
Figure FDA00033766418400000320
表示实际状态下X轴相对于床身坐标系的齐次坐标变换矩阵;
Figure FDA00033766418400000321
表示实际状态下Z轴相对于X轴的齐次坐标变换矩阵;
Figure FDA00033766418400000322
表示实际状态下A轴相对于Z轴的齐次坐标变换矩阵;
Figure FDA00033766418400000323
表示实际状态下Y轴相对于A轴的齐次坐标变换矩阵;
Figure FDA00033766418400000324
表示实际状态下B轴相对于Y轴的齐次坐标变换矩阵;
Figure FDA00033766418400000325
表示实际状态下C轴相对于床身R坐标系的齐次坐标变换矩阵,且
Figure FDA00033766418400000326
将热误差和几何误差映射到工具的空间误差为:
Figure FDA0003376641840000041
其中,[a,b,c,1]T和[i,j,k,0]T分别表示滚刀中心位置和轴姿态的齐次坐标;δ和ε分别表示位置误差和姿态误差,δx、δy和δz分别表示x、y、z方向的位置误差;εx、εy和εz分别表示x、y、z方向的姿态误差;
23)创建齿面误差与几何-热误差的映射模型
在理想状态下,理想齿面在工件坐标系中的坐标向量和单位法向量表示为:
Figure FDA0003376641840000042
Figure FDA0003376641840000043
其中,
Figure FDA0003376641840000044
表示理想状态下工件坐标系G相对于滚刀坐标系H的齐次变换矩阵,且
Figure FDA0003376641840000045
Figure FDA0003376641840000046
表示理想齿面在工件坐标系中的坐标向量;rh表示滚刀齿面的渐开线参数方程;
Figure FDA0003376641840000047
表示理想齿面在工件坐标系中的单位法向量;Nh表示滚刀齿面的单位法向量;
则在考虑了热误差和几何误差的实际状态下,工件坐标系中的坐标向量和单位法向量表示为:
Figure FDA0003376641840000048
Figure FDA0003376641840000049
其中,
Figure FDA00033766418400000410
表示实际状态下工件坐标系G相对于滚刀坐标系H的齐次变换矩阵,且
Figure FDA00033766418400000411
Figure FDA00033766418400000412
表示实际齿面在工件坐标系中的坐标向量;
Figure FDA00033766418400000413
表示实际齿面在工件坐标系中的单位法向量;
从而得到齿面误差与几何-热误差的映射模型为:
Figure FDA00033766418400000414
其中,n表示离散值。
6.根据权利要求5所述的霾-边-雾-云误差补偿系统,其特征在于:基于所述齿面误差与几何-热误差的映射模型,采用Morris方法识别出齿面误差的关键误差项,所述雾计算层中,结合所述齿面误差与几何-热误差的映射模型得到所述关键误差项的误差补偿值,并将关键误差项的误差补偿值传输至云计算层。
CN202111421117.3A 2021-11-26 2021-11-26 基于ont-gcn时空模型的热误差预测模型及建模方法和霾-边-雾-云误差补偿系统 Active CN114237154B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111421117.3A CN114237154B (zh) 2021-11-26 2021-11-26 基于ont-gcn时空模型的热误差预测模型及建模方法和霾-边-雾-云误差补偿系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111421117.3A CN114237154B (zh) 2021-11-26 2021-11-26 基于ont-gcn时空模型的热误差预测模型及建模方法和霾-边-雾-云误差补偿系统

Publications (2)

Publication Number Publication Date
CN114237154A true CN114237154A (zh) 2022-03-25
CN114237154B CN114237154B (zh) 2023-06-30

Family

ID=80751423

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111421117.3A Active CN114237154B (zh) 2021-11-26 2021-11-26 基于ont-gcn时空模型的热误差预测模型及建模方法和霾-边-雾-云误差补偿系统

Country Status (1)

Country Link
CN (1) CN114237154B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115475900A (zh) * 2022-09-19 2022-12-16 江兴(淮安)汽车部件有限公司 一种传动轴自动化精密锻造的方法和系统
CN115616978A (zh) * 2022-10-20 2023-01-17 重庆大学 基于sq-lstma的热误差预测模型、预测方法和控制系统
CN115793568A (zh) * 2023-02-06 2023-03-14 常州克迈特数控科技有限公司 一种基于三维建模的机床零件误差补偿控制系统

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105241382A (zh) * 2015-10-16 2016-01-13 合肥工业大学 时栅转台热误差测量系统及其测量计算方法
CN106896360A (zh) * 2017-04-21 2017-06-27 南京航空航天大学 一种sar信号处理算法的fpga实现方法
GB201907165D0 (en) * 2019-05-21 2019-07-03 Oxford Nanopore Tech Ltd Method for determining a consensus sequence of a target polymer
WO2020025984A1 (en) * 2018-08-01 2020-02-06 Pantazis Alexandros Method of use of a lidar device and operatively associated lidar data processing unit for providing real-time monitoring of meteorological parameters
CN111077869A (zh) * 2019-12-31 2020-04-28 国电九江发电有限公司 一种大数据智能控制布袋除尘器优化控制方法及系统
CN111090238A (zh) * 2018-10-24 2020-05-01 宝马股份公司 智能家居设备的控制方法、装置、设备、系统及存储介质
CN111240268A (zh) * 2020-01-14 2020-06-05 重庆大学 基于slstm神经网络的轴系统热误差建模方法和热误差补偿系统
CN112098715A (zh) * 2020-09-18 2020-12-18 宁波智诚祥科技发展有限公司 基于5g和修正gcn图神经网络的电能监控与预警系统
CN113051831A (zh) * 2021-04-01 2021-06-29 重庆大学 机床热误差自学习预测模型建模方法及基于数字孪生的机床热误差控制方法
CN113156822A (zh) * 2021-04-22 2021-07-23 重庆大学 基于Mist-边缘-雾-云计算的热误差预测系统及热误差补偿系统
CN113210911A (zh) * 2021-06-03 2021-08-06 重庆大学 基于图卷积网络的白车身点焊变形预测模型构建方法
CN113591395A (zh) * 2021-08-11 2021-11-02 重庆大学 热误差预测模型建模方法及基于霾-边缘-雾-云计算的智能热误差控制系统框架
US20220092391A1 (en) * 2021-12-07 2022-03-24 Santiago Miret System and method of using neuroevolution-enhanced multi-objective optimization for mixed-precision quantization of deep neural networks
CN115903516A (zh) * 2022-12-19 2023-04-04 重庆大学 基于stfgcn的时空热误差预测模型与时空热误差控制系统

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105241382A (zh) * 2015-10-16 2016-01-13 合肥工业大学 时栅转台热误差测量系统及其测量计算方法
CN106896360A (zh) * 2017-04-21 2017-06-27 南京航空航天大学 一种sar信号处理算法的fpga实现方法
WO2020025984A1 (en) * 2018-08-01 2020-02-06 Pantazis Alexandros Method of use of a lidar device and operatively associated lidar data processing unit for providing real-time monitoring of meteorological parameters
CN111090238A (zh) * 2018-10-24 2020-05-01 宝马股份公司 智能家居设备的控制方法、装置、设备、系统及存储介质
GB201907165D0 (en) * 2019-05-21 2019-07-03 Oxford Nanopore Tech Ltd Method for determining a consensus sequence of a target polymer
CN111077869A (zh) * 2019-12-31 2020-04-28 国电九江发电有限公司 一种大数据智能控制布袋除尘器优化控制方法及系统
CN111240268A (zh) * 2020-01-14 2020-06-05 重庆大学 基于slstm神经网络的轴系统热误差建模方法和热误差补偿系统
CN112098715A (zh) * 2020-09-18 2020-12-18 宁波智诚祥科技发展有限公司 基于5g和修正gcn图神经网络的电能监控与预警系统
CN113051831A (zh) * 2021-04-01 2021-06-29 重庆大学 机床热误差自学习预测模型建模方法及基于数字孪生的机床热误差控制方法
CN113156822A (zh) * 2021-04-22 2021-07-23 重庆大学 基于Mist-边缘-雾-云计算的热误差预测系统及热误差补偿系统
CN113210911A (zh) * 2021-06-03 2021-08-06 重庆大学 基于图卷积网络的白车身点焊变形预测模型构建方法
CN113591395A (zh) * 2021-08-11 2021-11-02 重庆大学 热误差预测模型建模方法及基于霾-边缘-雾-云计算的智能热误差控制系统框架
US20220092391A1 (en) * 2021-12-07 2022-03-24 Santiago Miret System and method of using neuroevolution-enhanced multi-objective optimization for mixed-precision quantization of deep neural networks
CN115903516A (zh) * 2022-12-19 2023-04-04 重庆大学 基于stfgcn的时空热误差预测模型与时空热误差控制系统

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
GUI HQ 等: "Mist-edge-fog-cloud computing system for geometric and thermal error prediction and compensation of worm gear machine tools based on ONT-GCN spatial -temporal model", 《MECHANICAL SYSTEMS AND SIGNAL PROCESSING》, vol. 184, pages 109682 *
YE DAI 等: "Modeling of thermal error elecric spindle based on KELM ameliorated by snake optimization", 《CASE STUDIES IN THERMAL ENGINEERING》, vol. 40, pages 102504 *
杨军 等: "时序分析在电主轴热误差建模中的应用", 《计算机集成制造系统》, vol. 21, no. 5, pages 1359 - 1367 *
申建广;陶涛;梅雪松;刘星;马星星;杜志国;: "一种数控滚齿机工作台动态特性建模方法及实验分析", 西安交通大学学报, no. 12, pages 6 - 12 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115475900A (zh) * 2022-09-19 2022-12-16 江兴(淮安)汽车部件有限公司 一种传动轴自动化精密锻造的方法和系统
CN115475900B (zh) * 2022-09-19 2023-10-27 江兴(淮安)汽车部件有限公司 一种传动轴自动化精密锻造的方法和系统
CN115616978A (zh) * 2022-10-20 2023-01-17 重庆大学 基于sq-lstma的热误差预测模型、预测方法和控制系统
CN115793568A (zh) * 2023-02-06 2023-03-14 常州克迈特数控科技有限公司 一种基于三维建模的机床零件误差补偿控制系统

Also Published As

Publication number Publication date
CN114237154B (zh) 2023-06-30

Similar Documents

Publication Publication Date Title
Aivaliotis et al. Degradation curves integration in physics-based models: Towards the predictive maintenance of industrial robots
Liu et al. Thermally-induced error compensation of spindle system based on long short term memory neural networks
Ding et al. State of AI-based monitoring in smart manufacturing and introduction to focused section
Davim Statistical and computational techniques in manufacturing
CN114237154A (zh) 基于ont-gcn时空模型的热误差预测模型及建模方法和霾-边-雾-云误差补偿系统
Ren et al. MCTAN: A novel multichannel temporal attention-based network for industrial health indicator prediction
Papananias et al. A Bayesian framework to estimate part quality and associated uncertainties in multistage manufacturing
Gui et al. Mist-edge-fog-cloud computing system for geometric and thermal error prediction and compensation of worm gear machine tools based on ONT-GCN spatial–temporal model
Majeske Approval criteria for multivariate measurement systems
CN113591395B (zh) 热误差预测模型建模方法及基于霾-边缘-雾-云计算的智能热误差控制系统框架
CN114002998B (zh) 齿形磨床几何误差和热误差同步控制方法
Guo et al. Study of static thermal deformation modeling based on a hybrid CNN-LSTM model with spatiotemporal correlation
CN110007645A (zh) 一种基于动力学及深度神经网络的进给系统混合建模方法
Yuan et al. Correlation analysis-based thermal error control with ITSA-GRU-A model and cloud-edge-physical collaboration framework
Wegener et al. Dawn of new machining concepts:: Compensated, intelligent, bioinspired
CN113589757A (zh) 成形磨削综合误差模型创建方法及磨齿机几何-热误差控制方法
Chengyang et al. Thermal error modeling of rotary axis based on convolutional neural network
Li et al. A wiener-based remaining useful life prediction method with multiple degradation patterns
CN111985149A (zh) 一种基于卷积网络的五轴机床旋转轴热误差建模方法
Lee et al. In-line predictive monitoring framework
Khezri et al. Evolutionary cost-tolerance optimization for complex assembly mechanisms via simulation and surrogate modeling approaches: application on micro gears
Yao et al. Synthetic error modeling for NC machine tools based on intelligent technology
Zhao et al. Modeling and prediction of full-term thermal error in linear axis of machine tools based on MSTGCN-A
CN112990601B (zh) 基于数据挖掘的蜗轮加工精度自愈系统及方法
CN115903516A (zh) 基于stfgcn的时空热误差预测模型与时空热误差控制系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant