CN115422670B - 一种基于切削力与夹紧力时空学习的薄壁件加工变形预测方法 - Google Patents
一种基于切削力与夹紧力时空学习的薄壁件加工变形预测方法 Download PDFInfo
- Publication number
- CN115422670B CN115422670B CN202210950004.0A CN202210950004A CN115422670B CN 115422670 B CN115422670 B CN 115422670B CN 202210950004 A CN202210950004 A CN 202210950004A CN 115422670 B CN115422670 B CN 115422670B
- Authority
- CN
- China
- Prior art keywords
- clamping force
- time
- deformation
- cutting force
- model
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 68
- 238000003754 machining Methods 0.000 title claims abstract description 19
- 238000012545 processing Methods 0.000 claims abstract description 40
- 239000011159 matrix material Substances 0.000 claims abstract description 36
- 230000008569 process Effects 0.000 claims abstract description 24
- 238000013528 artificial neural network Methods 0.000 claims abstract description 13
- 238000012549 training Methods 0.000 claims abstract description 13
- 125000004122 cyclic group Chemical group 0.000 claims abstract description 12
- 238000013507 mapping Methods 0.000 claims description 26
- 230000006870 function Effects 0.000 claims description 25
- 230000008859 change Effects 0.000 claims description 15
- 230000009466 transformation Effects 0.000 claims description 14
- 230000004913 activation Effects 0.000 claims description 9
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 230000002596 correlated effect Effects 0.000 claims description 5
- 230000000875 corresponding effect Effects 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims description 4
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 230000003213 activating effect Effects 0.000 claims 2
- 238000005452 bending Methods 0.000 abstract description 8
- 238000013527 convolutional neural network Methods 0.000 description 5
- 239000007787 solid Substances 0.000 description 4
- 238000004088 simulation Methods 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000014759 maintenance of location Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/24—Sheet material
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Abstract
本发明提出一种基于切削力与夹紧力时空学习的薄壁件加工变形预测方法,该方法将实际加工过程中采集到的切削力与夹紧力作为主要外载荷,首先建立切削力与夹紧力的时变关联模型;接着将被加工零件进行几何模型参数化,获得可以表示零件几何结构的参数化矩阵,用来反映零件的弯曲刚度;之后针对变形形成过程中上述外载荷、几何与变形之间的时序与空间关联特征,建立卷积网络与循环神经网络结合的时空学习模型,将切削力、夹紧力、工件几何与变形信息作为数据样本训练该网络,实现加工变形预测。本发明可以用在对切削力与夹紧力敏感的加工变形中,训练完成后变形预测模型可以达到较高的效率,同时适用性较好。
Description
技术领域
本发明涉及一种薄壁件加工变形预测方法,特别是一种基于切削力与夹紧力时空学习的薄壁件加工变形预测方法。
背景技术
薄壁结构零件在航空工业中被广泛使用,以减轻重量、提高结构强度和性能,如航空整体壁板、整体框、梁、壁板等。由于具有材料去除率高,薄壁结构和低刚度特性,这些零件在加工过程中很容易发生变形。在实际加工中,影响变形的因素很多,包括切削力、夹具系统、工件的几何形状和刚度,以及工件的残余应力等,这些因素相互作用使得加工变形更加复杂难以预测,容易造成几何误差甚至零件报废,严重影响了加工质量、加工效率。
对加工变形的预测不仅是后续进行加工变形控制的基础,也是保证加工质量、提升加工效率的前提。针对加工过程中由于材料去除工件内部应力释放而产生的整体加工变形相关研究很多,主要是以数值模拟和分析模型为主的加工过程仿真模型,且都需要获得精确的初始残余应力大小和分布。但是,由于残余应力分布不均匀、现有残余应力测量方法的限制,以及理论假设与实际加工场景较大的差异,导致现有变形预测方法很难准确预测变形,或者模型适用性受到限制。
由于加工过程中,由刀具、工件、机床和夹具组成的是一个动态平衡系统,切削力与夹紧力是工件加工过程中受到的主要的外载荷,而且外载荷与工件内部的变形抗力处于动态平衡。也就是说,对于这种平衡系统,可以用外界载荷代替难以测量的工件内部的变形抗力。实际上,利用夹具来减小或抑制变形已有很多研究,但往往都是特定加工场景下的专用夹具,或者需要复杂高效的加工过程仿真或解析模型,方法的通用性和预测效果仍受限制。
发明内容
鉴于此,为解决现有技术存在的问题,本发明提供一种基于切削力与夹紧力时空学习的薄壁件加工变形预测方法。该方法将实际加工过程中采集到的切削力与夹紧力作为主要外载荷,首先建立切削力与夹紧力的时变关联模型;接着将被加工零件进行几何模型参数化,获得可以表示零件几何结构的参数化矩阵,用来反映零件的弯曲刚度;之后针对变形形成过程中上述外载荷、几何与变形之间的时序与空间关联特征,建立卷积网络与循环神经网络结合的时空学习模型,将切削力、夹紧力、工件几何与变形信息作为数据样本训练该网络,实现加工变形预测。
本发明的技术方案为:
所述一种基于切削力与夹紧力时空学习的薄壁件加工变形预测方法,包括以下步骤:
步骤1:根据薄壁件加工过程中实际采集的切削力与夹紧力变化数据,建立切削力与夹紧力的时变模型,进而建立表示切削力与夹紧力关联变化的时变关联模型;
步骤2:采用曲面映射结合保刚性变化将被加工的薄壁件进行几何模型参数化;
步骤3:建立基于时空学习的变形预测模型:
首先,将关联后的切削力与夹紧力,以及工件几何参数化模型输入卷积网络进行空间关系学习:
HΔt=fc(Wc(L(Δt),D)+bc)
其中,HΔt是经过卷积网络学习后的输出;fc是激活函数,Wc,bc是需要学习的网络参数;
接着,将若干个按照时间序列排列,且经过卷积网络学习后的输出,经过扩展处理后,输入循环神经网络中,对与变形相关的时序关系进行学习:
将网络输出的变形和实际加工变形的均方差作为误差函数:
其中,Et是在时间t的损失值,E是整个训练数据的损失值;Defn为经过n个时间间隔测得的真实工件变形值,N为以真实变形值为标签的训练样本数量;
根据误差计算GRU网络中每个权重的梯度和偏置梯度,同时采用Aadm算法更新网络参数,直到损失函数收敛,完成变形预测模型的建立;
步骤4:利用训练好的变形预测模型实现变形预测。
进一步的,步骤1中,切削力时变模型为
F(t)=(u(t),v(t),w(t),f(t))
其中,u(t),v(t),w(t)表示随着刀具移动的刀具坐标系,其原点位于刀尖,u轴、v轴和w轴分别代表瞬时进给方向、工件表面法向量和刀具轴向量,f(t)表示随时间变化的切削力;
夹紧力时变模型为
cF(t)=(xi,yi,zi,cf(t))
其中,xi,yi,zi表示第i个夹紧力测量位置在工件坐标系下对应的坐标,cf(t)表示随时间变化的夹紧力;
加工过程中切削力与夹紧力空间转化关系为:
[x,y,z,1]T=Trans[u,v,w,1]T
其中,Trans表示刀具坐标系向工件坐标系转换的齐次变换矩阵,上标T表示转置;在时间上,采用相邻刀位点之间的时间间隔,分别对切削力与夹紧力取平均值,实现时间上的一致性。
进一步的,齐次变换矩阵Trans表示为:
式中,3×3的矩阵RT表示刀具坐标系相对于工件坐标系的坐标旋转矩阵,由刀具坐标轴矢量决定,包括瞬时进给方向、工件切削表面法向量和刀轴矢量;3×1的矩阵QT表示刀具原点相对于机床的位置向量;γ和η分别表示透视变换和比例系数,这里η取1。
进一步的,切削力与夹紧力关联变化的时变关联模型以关联矩阵L(t)表示:
其中,L(t)是切削力与夹紧力关联后的外载荷模型矩阵表达形式,(x(t),y(t),z(t))是经过齐次变换后的切削力坐标。
进一步的,步骤2中,采用曲面映射结合保刚性变化实现薄壁件几何模型参数化的过程为:
首先采用曲面映射中的网格生成方法将薄壁件几何模型映射到参数域:
然后根据保刚性参数化方法,将参数化后的几何模型上每个三角形曲面块不变形地投影到一个单独的局部坐标系中,然后将分散在各自局部坐标系中的三角形块粘合在一起,形成一个完整的二维展开平面,得到薄壁件几何模型的二维数组表达式D(ui,vi),(ui,vi)表示二维平面上点的坐标。
进一步的,映射函数采用Hermite超限插值方法来实现工件几何模型到参数域的映射。
有益效果
本发明将加工过程中切削力与夹紧力的变化作为施加在工件上的主要外载荷,来表示难以测量的且随着加工变化的工件残余应力分布状态,用工件几何参数化模型来反映工件的弯曲刚度。即将影响加工变形但难测量或者难以获得的加工过程物理量,转换为用其他几何-物理量来共同描述。接着充分考虑到工件变形形成过程中时序空间相关性,采用卷积网络与循环神经网络结合的时空学习来对切削力夹紧力、工件几何与变形进行建模,实现变形预测。该方法可以用在对切削力与夹紧力敏感的加工变形中,训练完成后变形预测模型可以达到较高的效率,同时适用性较好。
本发明的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
具体实施方式
本发明提供一种基于切削力与夹紧力时空学习的薄壁件加工变形预测方法,该方法将实际加工过程中采集到的切削力与夹紧力作为主要外载荷,首先建立切削力与夹紧力的时变关联模型;接着将被加工零件进行几何模型参数化,获得可以表示零件几何结构的参数化矩阵,用来反映零件的弯曲刚度;之后针对变形形成过程中上述外载荷、几何与变形之间的时序与空间关联特征,建立卷积网络与循环神经网络结合的时空学习模型,将切削力、夹紧力、工件几何与变形信息作为数据样本训练该网络,实现加工变形预测。
本发明具体包括下述步骤:
步骤1:根据薄壁件加工过程中实际采集的切削力与夹紧力变化数据,建立切削力与夹紧力的时变模型,进而建立表示二者关联变化的时变关联模型。
切削力时变模型可以表示为:
F(t)=(u(t),v(t),w(t),f(t))
其中u(t),v(t),w(t)表示随着刀具移动的刀具坐标系,其原点位于刀尖,u轴、v轴和w轴分别代表瞬时进给方向、工件表面法向量和刀具轴向量。f(t)表示随时间变化的切削力。
对于施加在工件上的夹紧力来说,其作用位置固定,夹紧力随时间变化,夹紧力时变模型可以表示为:
cF(t)=(xi,yi,zi,cf(t))
其中xi,yi,zi表示第i个夹紧力测量位置在工件坐标系下对应的坐标,cf(t)表示随时间变化的夹紧力。
显然,切削力与夹紧力同时作用于工件,但作用过程中切削力的作用位置随着刀具位置变化而变化,夹紧力作用位置则固定在工件坐标系下某几个确定的位置,即二者施加的坐标系不同,为了准确描述加工过程中工件所受到的复杂外力状态,给出二者转化关系:
[x,y,z,1]T=Trans[u,v,w,1]T
其中,Trans表示刀具坐标系向工件坐标系转换的齐次变换矩阵,上标T为向量的转置。这样,在空间上将切削力与夹紧力关联起来。在时间上,采用相邻刀位点之间的时间间隔,分别对切削力以及夹紧力取平均值,来保证时间上的一致性。
步骤2:将工件几何模型参数化,获得表示工件几何结构的参数矩阵,反映零件的弯曲刚度,作为后续变形预测模型的输入。因此,工件几何模型参数化需要能够在不丢失工件几何拓扑关系的同时进行变形表示。这里采用曲面映射结合保刚性变化来实现几何模型参数化。
对于工件几何模型,首先采用曲面映射中的网格生成方法将工件几何模型映射到参数域,可以表示为:
接着,利用保刚性参数化方法,将参数化后的几何模型上每个三角形曲面块不变形地投影到一个单独的局部坐标系中,然后将分散在各自局部坐标系中的三角形块粘合在一起,形成一个完整的二维展开平面,即可得到几何模型的二维数组表达式D(ui,vi),(ui,vi)表示二维平面上点的坐标。
步骤3:建立基于时空学习的变形预测模型:
切削力与夹紧力,工件变形状态以及内部变形抗力在工件本身的几何空间上处于动态平衡。这样将工件几何状态作为空间基础,切削力与夹紧力以及工件变形状态在该空间上不断叠加且动态平衡。因此,可以建立一种同时包含对加工过程空间动态平衡,以及时间上不断叠加变化且平衡的时空学习网络。
首先,步骤1中将切削力与夹紧力关联后,考虑二者的空间坐标,表示为矩阵形式L(t);接着,步骤2中工件几何模型参数化后表示为矩阵D;然后,建立深度卷积网络,包含两个独立的浅层卷积过程,之后用多层卷积将二者结果融合。可以表示为:
HΔt=CNN(L(Δt),D)
其中,HΔt表示以时间间隔t为切削力、夹紧力的取值间隔,结合工件几何参数化后的矩阵D,经过卷积神经网络后得到的输出;这里Δt取两个相邻刀位点的时间间隔;CNN表示卷积操作,L(Δt)表示时间间隔Δt下的切削力与夹紧力关联矩阵。
再将若干个经过时间间隔Δt卷积后的输出输入循环神经网络,这里采用GRU门控循环单元网络,可以表示为:
通过卷积网络与循环神经网络结合,得到时空学习模型,将切削力、夹紧力、工件几何与变形信息作为数据样本训练该时空学习模型,得到具有时序空间关联学习的变形预测模型。之后利用训练好的变形预测模型实现变形预测。
下面详细描述本发明的实施例,所述实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
1、切削力与夹紧力的时变关联。
首先,建立切削力时变模型,可以表示为:
F(t)=(u(t),v(t),w(t),f(t))
其中,u(t),v(t),w(t)表示随着刀具移动的刀具坐标系,其原点位于刀尖,u轴、v轴和w轴分别代表瞬时进给方向、工件表面法向量和刀具轴向量。f(t)表示随时间变化的切削力。
建立夹紧力时变模型,可以表示为:
cF(t)=(xi,yi,zi,cf(t))
其中,xi,yi,zi表示第i个夹紧力测量位置在工件坐标系下对应的坐标,cf(t)表示随时间变化的夹紧力。
接着,将刀具坐标系向工件坐标系转换。这里采用的齐次变换矩阵表示从切削力坐标到机床坐标的映射。设FCS为刀具坐标系,WCS为工件坐标系(或者机床坐标系),两者之间的转换关系可以用下式矩阵Trans表示:
式中,3×3的矩阵RT表示FCS相对于MCS的坐标旋转矩阵,主要由刀具坐标轴矢量决定,包括瞬时进给方向、工件切削表面法向量和刀轴矢量;3×1的矩阵QT表示刀具原点相对于机床的位置向量;γ和η分别表示透视变换和比例系数,这里η取1。
之后,将二者在统一的工件坐标系下进行关联,获得切削力与夹紧力时变关联的外载荷模型,以关联矩阵L(t)表示。
其中,L(t)是切削力与夹紧力关联后的外载荷模型矩阵表达形式,(x(t),y(t),z(t))是经过齐次变换后的切削力坐标。
2、工件几何模型的参数化。
首先采用曲面映射中的网格生成方法将几何模型映射到参数域,可以表示为:
其中,S(x,y,z)表示在笛卡尔坐标系下的三维几何实体模型,D(u,v,w)表示参数域坐标系下的几何模型,F(x,y,z,u,v,w)表示映射函数,符号表示一对一映射关系。具体来说,映射函数采用Hermite超限插值方法来实现工件几何模型到参数域的映射。
接着,根据保刚性参数化方法,将参数化后的几何模型上每个三角形曲面块不变形地投影到一个单独的局部坐标系中,然后将分散在各自局部坐标系中的三角形块粘合在一起,形成一个完整的二维展开平面,即可得到几何模型的二维数组表达。
由于二维平面上的点不一定均匀,因此进一步对在二维平面上的几何模型进行均匀采样,使得几何模型在二维平面上可以用矩阵表示,经过采样后的工件几何参数化模型具体可以表示为:
其中,(ui,vi)表示二维平面上点的坐标,i对应参数化后几何模型上三角切片顶点的索引号。
3、基于时空学习的变形预测。
首先,将关联后的切削力与夹紧力,以及工件几何参数化模型输入卷积网络进行空间关系学习。
HΔt=fc(Wc(L(Δt),D)+bc)
其中,HΔt是经过卷积网络学习后的输出;fc是激活函数,这里采用ReLU激活函数,具体为f(x):=max(0,x);Wc,bc是需要学习的网络参数。
接着,将若干个按照时间序列排列,且经过卷积网络学习后的输出,经过扩展处理后,输入循环神经网络中,对与变形相关的时序关系进行学习。
其中,为GRU网络在第n个时间间隔Δt后,预测得到的变形输出;fg是激活函数,这里采用Sigmoid函数作为激活函数,具体为sigmoid(x)=1/(1+e-x);是网络前一刻的输出,是HΔt经过扩展处理后符合GRU网络的输入;Wg,U是需要学习的网络参数。
将网络输出的变形和实际加工变形的均方差作为误差函数,具体如下:
其中,Et是在时间t的损失值,E是整个训练数据的损失值;Defn为经过n个时间间隔测得的真实工件变形值,N为以真实变形值为标签的训练样本数量。
根据误差计算GRU网络中每个权重的梯度和偏置梯度,同时采用Aadm算法更新网络参数,直到损失函数收敛,完成变形预测模型的建立。
本发明首先根据加工过程中内应力与外载荷在工件上处于动平衡状态,将切削力与夹紧力的变化用来表征难测量的工件内应力状态变化。虽然切削力与夹紧力可以通过传感器采集,但对于工件受外载荷来说,工件受到的夹紧力位置相对不变,而受到的切削力随着刀具轨迹不断变化,即切削力与夹紧力作用过程并不一致。
因此,采用齐次变换矩阵将刀具坐标系下的切削力转化到工作坐标系下,并将坐标转化后的切削力与夹紧力进行动态关联,形成带有空间坐标信息且随着时间变化,具有时空一致性的外载荷模型。
其次,由于工件几何是三维的复杂结构,很难直接进行变形相关的数据建模与计算。考虑到工件是实体结构,工件表面信息足以表达零件不同位置的弯曲刚度,因此本发明通过超限插值与保刚度结合的方法,将工件三维几何模型通过网格生成,保刚度映射到二维平面并均匀采样,从而实现用参数化后的工件几何模型作为工件弯曲刚度的表征。
这样参数化后的几何模型既满足后续变形预测模型输入,同时又减少计算量。该方法可以灵活的控制半精加工及精加工的余量,同时也是利用数据驱动模型中对复杂几何模型进行参数化的一次尝试。
之后,在加工过程中工件作为一个连续整体,其加工变形过程在工件几何结构上是连续的,在时序上具有自相关性。这里针对加工变形过程中的时空相关性特点,提出设计全卷积神经网络对工件几何、外载荷及变形信息进行空间相关性学习。与卷积神经网络的区别在于,将卷积神经网络中的全连接层改为卷积层,可以使卷积网络从全局更好地学习整个工件几何-载荷-变形的空间关联特征(关系),更符合实际加工过程中工件整体变形的连续性。
同时,将不同时间段工件几何、外载荷以及变形信息,经过全卷积网络学习后的空间关联特征,看作工件变形演变过程的一个个时间切片,用循环神经网络中的GRU网络来建立这些切片之间的时序关系,模型参数较少且效率高,短时序的网络也符合实际易变形工件工序间变形预测与控制要求。
综上所述,本发明通过建立包含实际加工过程切削力与夹紧力的外载荷模型,代替难以测量的工件残余应力分布的变化,进一步将几何模型参数化来表示工件弯曲刚度。之后针对加工变形过程的时空相关性特点,将全卷积神经网络与循环神经网络结合,建立具有时序空间关联学习的变形预测模型。利用实际采集的切削力与夹紧力,而不是解析或者仿真的数据,将二者动态关联后更能反映工件真实变形状态,也能提高模型的通用性,且训练完成后的变形预测模型计算效率也更高,相比直接用端到端网络模型预测精度提升20%左右。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (6)
1.一种基于切削力与夹紧力时空学习的薄壁件加工变形预测方法,其特征在于:包括以下步骤:
步骤1:根据薄壁件加工过程中实际采集的切削力与夹紧力变化数据,建立切削力与夹紧力的时变模型,进而建立表示切削力与夹紧力关联变化的时变关联模型;
其中切削力时变模型为
F(t)=(u(t),v(t),w(t),f(t))
其中,u(t),v(t),w(t)表示随着刀具移动的刀具坐标系,其原点位于刀尖,u轴、v轴和w轴分别代表瞬时进给方向、工件表面法向量和刀具轴向量,f(t)表示随时间变化的切削力;
夹紧力时变模型为
cF(t)=(xi,yi,zi,cf(t))
其中,xi,yi,zi表示第i个夹紧力测量位置在工件坐标系下对应的坐标,cf(t)表示随时间变化的夹紧力;
加工过程中切削力与夹紧力空间转化关系为:
[x,y,z,1]T=Trans[u,v,w,1]T
其中,Trans表示刀具坐标系向工件坐标系转换的齐次变换矩阵,上标T表示转置;在时间上,采用相邻刀位点之间的时间间隔,分别对切削力与夹紧力取平均值,实现时间上的一致性;
步骤2:采用曲面映射结合保刚性变化将被加工的薄壁件进行几何模型参数化;
步骤3:建立基于时空学习的变形预测模型:
首先,将关联后的切削力与夹紧力,以及工件几何参数化模型输入卷积网络进行空间关系学习:
HΔt=fc(Wc(L(Δt),D)+bc)
其中,HΔt是经过卷积网络学习后的输出;fc是激活函数,Wc,bc是需要学习的网络参数;L(Δt)表示时间间隔Δt下的切削力与夹紧力关联矩阵,D表示工件几何参数化后的矩阵;
接着,将若干个按照时间序列排列,且经过卷积网络学习后的输出,经过扩展处理后,输入循环神经网络中,对与变形相关的时序关系进行学习:
将网络输出的变形和实际加工变形的均方差作为误差函数:
其中,Et是在时间t的损失值,E是整个训练数据的损失值;Defn为经过n个时间间隔测得的真实工件变形值,N为以真实变形值为标签的训练样本数量;
根据误差计算GRU网络中每个权重的梯度和偏置梯度,同时采用Aadm算法更新网络参数,直到损失函数收敛,完成变形预测模型的建立;
步骤4:利用训练好的变形预测模型实现变形预测。
4.根据权利要求1所述一种基于切削力与夹紧力时空学习的薄壁件加工变形预测方法,其特征在于:步骤2中,采用曲面映射结合保刚性变化实现薄壁件几何模型参数化的过程为:
首先采用曲面映射中的网格生成方法将薄壁件几何模型映射到参数域:
然后根据保刚性参数化方法,将参数化后的几何模型上每个三角形曲面块不变形地投影到一个单独的局部坐标系中,然后将分散在各自局部坐标系中的三角形块粘合在一起,形成一个完整的二维展开平面,得到薄壁件几何模型的二维数组表达式D(ui,vi),(ui,vi)表示二维平面上点的坐标。
5.根据权利要求4所述一种基于切削力与夹紧力时空学习的薄壁件加工变形预测方法,其特征在于:映射函数采用Hermite超限插值方法来实现工件几何模型到参数域的映射。
6.根据权利要求1所述一种基于切削力与夹紧力时空学习的薄壁件加工变形预测方法,其特征在于:激活函数fc采用ReLU激活函数,激活函数fg采用Sigmoid函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210950004.0A CN115422670B (zh) | 2022-08-09 | 2022-08-09 | 一种基于切削力与夹紧力时空学习的薄壁件加工变形预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210950004.0A CN115422670B (zh) | 2022-08-09 | 2022-08-09 | 一种基于切削力与夹紧力时空学习的薄壁件加工变形预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115422670A CN115422670A (zh) | 2022-12-02 |
CN115422670B true CN115422670B (zh) | 2023-05-16 |
Family
ID=84195701
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210950004.0A Active CN115422670B (zh) | 2022-08-09 | 2022-08-09 | 一种基于切削力与夹紧力时空学习的薄壁件加工变形预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115422670B (zh) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114444347A (zh) * | 2021-12-31 | 2022-05-06 | 天津大学 | 基于有限元分析的大型曲面薄壁件镜像铣削变形预测方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111079690B (zh) * | 2019-12-27 | 2022-05-20 | 华中科技大学 | 基于堆栈稀疏自动编码网络的主轴和工件振动预测方法 |
CN111723440B (zh) * | 2020-05-12 | 2022-04-12 | 华中科技大学 | 一种薄壁件加工精度预测混合建模方法 |
CN111783241B (zh) * | 2020-06-16 | 2022-09-30 | 大连理工大学 | 一种薄壁微铣削变形预测方法 |
CN112784451A (zh) * | 2020-11-27 | 2021-05-11 | 北京工业大学 | 一种基于有限元和支持向量机的薄壁件加工变形预测方法 |
CN112668227B (zh) * | 2020-12-31 | 2024-02-02 | 华中科技大学 | 一种薄壁件让刀变形误差预测模型建立方法及其应用 |
CN113414638B (zh) * | 2021-06-04 | 2023-02-10 | 西北工业大学 | 一种基于铣削力时序图深度学习的变工况铣刀磨损状态预测方法 |
CN113703395B (zh) * | 2021-07-07 | 2022-10-21 | 西北工业大学 | 一种面向加工变形控制的变工况铣削加工夹紧力预测方法 |
-
2022
- 2022-08-09 CN CN202210950004.0A patent/CN115422670B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114444347A (zh) * | 2021-12-31 | 2022-05-06 | 天津大学 | 基于有限元分析的大型曲面薄壁件镜像铣削变形预测方法 |
Non-Patent Citations (1)
Title |
---|
薄壁件加工变形控制快速仿真平台开发;陈蔚芳;陈华;楼佩煌;郑会龙;;计算机集成制造系统(第02期);第321页-第327页 * |
Also Published As
Publication number | Publication date |
---|---|
CN115422670A (zh) | 2022-12-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111964575B (zh) | 一种面向移动机器人铣削加工的数字孪生建模方法 | |
CN110262394B (zh) | 一种数控加工中轮廓误差的补偿方法 | |
CN111596614B (zh) | 基于云边协同的运动控制误差补偿系统及方法 | |
US20070106421A1 (en) | Method for optimising the performance of a robot | |
CN111723440B (zh) | 一种薄壁件加工精度预测混合建模方法 | |
CN109822576B (zh) | 一种机器人加工虚拟夹具的生成方法 | |
Gohari et al. | A digital twin for integrated inspection system in digital manufacturing | |
CN112433507B (zh) | 基于lso-lssvm的五轴数控机床热误差综合建模方法 | |
CN110103071B (zh) | 一种变形复杂零件的数字化寻位加工方法 | |
Gorkavyy et al. | Intelligent system for prognostication and optimization of power expenses of technological processes at robotized productions | |
CN108304687A (zh) | 一种预测薄壁复杂曲面回转件车削加工变形的方法 | |
CN114925462B (zh) | 一种基于切削力与刚度关联演变的薄壁件加工变形预测方法 | |
CN115422670B (zh) | 一种基于切削力与夹紧力时空学习的薄壁件加工变形预测方法 | |
CN113799137A (zh) | 一种基于神经网络的移动加工机器人快速视觉定位方法 | |
Geng et al. | Analysis of nonlinear error caused by motions of rotation axes for five-axis machine tools with orthogonal configuration | |
Suneel et al. | Intelligent tool path correction for improving profile accuracy in CNC turning | |
CN109446721B (zh) | 基于标识符软件线程执行顺序排列的机床工艺交互算法 | |
CN107728578A (zh) | 一种基于加工变形监测数据的加工顺序自适应调整方法 | |
CN112157654A (zh) | 一种机器人加工系统定位误差的优化方法 | |
CN116756851A (zh) | 基于nffd背景网格的参数化网格变形方法及系统 | |
CN108197809B (zh) | 基于维度最优转换的权值共享深度网络的实时调度方法 | |
CN113642129B (zh) | 一种工件校正载荷快速施加及变形有限元获取方法 | |
CN115407680A (zh) | 一种并行化的高精度时域气动弹性控制耦合仿真方法 | |
Mohring et al. | Simulation aided design of intelligent machine tool components | |
CN109725595A (zh) | 工件的加工路径的补偿方法、加工方法及工件 |
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 |