CN103440416A - 基于扩展误差流的叶片加工过程误差预测方法 - Google Patents
基于扩展误差流的叶片加工过程误差预测方法 Download PDFInfo
- Publication number
- CN103440416A CN103440416A CN2013103767650A CN201310376765A CN103440416A CN 103440416 A CN103440416 A CN 103440416A CN 2013103767650 A CN2013103767650 A CN 2013103767650A CN 201310376765 A CN201310376765 A CN 201310376765A CN 103440416 A CN103440416 A CN 103440416A
- Authority
- CN
- China
- Prior art keywords
- delta
- centerdot
- error
- coordinate system
- formula
- 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.)
- Pending
Links
Images
Abstract
本发明公开了一种基于扩展误差流的叶片加工过程误差预测方法,用于解决现有误差流方法以刚体零件装配多工序过程或者加工多工序过程为研究对象仅考虑几何误差源的技术问题。技术方案是采用模型驱动进行过程监控,对几何定位误差源与受力变形误差源产生的误差及时进行预测,对误差耦合和传递机理进行了研究,对加工过程误差耦合建模方案进行了设计;通过网格划分的思想对叶片零件进行了微元划分,将变形分析转换为微元坐标系的坐标变换;建立了基于扩展误差流的多源多工序误差耦合模型,在此基础上给出了基于误差耦合模型的综合误差预测模型。帮助工作人员及时发现叶片的制造缺陷,进行校正和调整,减少返工次数,提高了加工质量和效率。
Description
技术领域
本发明涉及一种叶片加工过程误差预测方法,特别是涉及一种基于扩展误差流的叶片加工过程误差预测方法。
背景技术
由于叶片零件的特殊性,其成型加工与传统高刚性零件的成型加工有着本质不同,对各种误差源的波动变化更加敏感,工序间误差的传递、累积等交互作用现象更加复杂。因此研究叶片零件加工过程中误差耦合和传播情况,分析误差的耦合和传递形式、建立误差源与误差的对应关系是进行加工过程质量监控、溯源与调整的关键与前提,对叶片零件的质量改进有重要的意义。
文献Liu J.,Jin J.and Shi J.State space modeling for3-dimensional variationpropagation in rigid-body multistage assembly processes[J],IEEE Trans.Autom.Sci.Eng.,2009,274-290针对刚体零件构建了多工序装配过程的状态空间模型,给出了模型的系数矩阵,解决了装配制造系统中多个偏差流影响产品质量的问题,但是装配过程误差相对简单,是多个误差的简单的叠加,由于加工过程对材料进行切削,使得误差受到多个因素的耦合影响,因此这种方法不适用于多工序的加工过程。
文献Abellan-Nebot J.,Liu J.,Romero F.Limitations of the current state spacemodeling approach in multistage machining processes due to operation variations[C].3rdManufacturing Engineering Society International Conference,Alcoy,Spain,2009利用状态空间方程建立描述多工序加工过程的空间模型,针对基准、夹具等误差源建立了刚体零件的通用的误差流模型,但是没有考虑到物理误差源的影响,并且由于刚体零件不具有局部的变形,因此其不适用与具有局部变形的弱刚性零件的多工序误差建模。
发明内容
为了克服现有误差流技术以刚体零件装配多工序过程或者加工多工序过程为研究对象仅考虑几何误差源的局限性,本发明提供一种基于扩展误差流的叶片加工过程误差预测方法。该方法采用模型驱动进行过程监控,对几何定位误差源与受力变形误差源产生的误差及时进行预测,对误差耦合和传递机理进行了研究,对加工过程误差耦合建模方案进行了设计;通过网格划分的思想对叶片零件进行了微元划分,将变形分析转换为微元坐标系的坐标变换;对加工过程误差与误差源因素的映射关系进行了深入的研究分析,从加工工序流的角度采用扩展误差流的理论,建立了基于扩展误差流(SoV)的多源多工序误差耦合模型,在此基础上给出了基于误差耦合模型的综合误差预测模型。帮助工作人员及时发现叶片的制造缺陷,进行校正和调整,减少返工次数,以提高加工质量和效率。
本发明解决其技术问题所采用的技术方案是:一种基于扩展误差流的叶片加工过程误差预测方法,其特点是包括以下步骤:
步骤一、确定叶片加工过程误差源分别为定位基准精度产生的误差;安装产生的误差;刀具位姿变化产生的误差;切削力产生的变形误差;噪声误差。对误差源进行建模,并对加工过程中的因素约定如下:多工序加工过程中将一道工序的加工过程分为加工前和加工中两个阶段。基准误差、夹具安装几何误差以及刀具倾斜误差对工件产生定位误差,定义为几何定位误差。加工中的弹性变形是由切削力引起的。
步骤二、基于微元刚体坐标变换的几何定位误差计算。
定义五个坐标系,分别是全局坐标系{O}、工件坐标系{ws}、特征坐标系{fs}、微元坐标系{ps}以及刀具坐标系{T}。
根据各个对象在空间的位姿关系,由全局坐标系{O}、工件坐标系{ws}、特征坐标系{Fs}、微元坐标系{ps}以及刀具坐标系共同构成空间运动链:
式中,表示刀具坐标系{T}相对于全局坐标系{O}的坐标变换矩阵,表示工件坐标系{ws}相对于全局坐标系{O}的坐标变换矩阵,表示特征坐标系{fs}相对于工件坐标系{ws}的坐标变换矩阵,表示微元坐标系{ps}相对于特征坐标系{fs}的坐标变换矩阵表示刀具坐标系{T}相对于微元坐标系{ps}的坐标变换矩阵。
运用微分运动原理得到:
式中,Δ表示微分。
式(2)的左边变换为
则得到:
将式(3)改写为:
忽略高阶项后得到以下表达式:
从而得到刀具在微元坐标系下的偏差:
由此得到刀具在特征坐标系下和在工件坐标系下的偏差分别为:
式中,表示微元坐标系相对于全局坐标系的旋转矩阵的名义值,表示特征坐标系相对于全局坐标系的旋转矩阵的名义值,表示工件坐标系相对于全局坐标系的旋转矩阵的名义值, 和表示反对称矩阵,为刀具相对于全局坐标系的偏差量,为工件相对于全局坐标系的偏差量,为特征相对于全局坐标系的偏差量,表示微元相对于全局坐标系的偏差量。
在装夹好工件时,接触点在工件坐标系{ws}、特征坐标系{fs}、微元坐标系{ps}与定位元件坐标系是重合的则工件与元件l的第i个接触点有四种方法描述:
则该式中的前三个等式与最后一个相等,即Fi=fi
采用一阶泰勒级数展开,省去高阶项得到:
将式(10)写成:
则得到:
由于夹具固定在机床上,假定忽略掉夹具的位置和姿态误差,则得到:
从而得到:
假设工件上第i个接触点的几何位置Δrpsi足够小,则忽略工件上基准面接触点的位置误差:
假定在i个接触点在法向上有误差Δrni,并且坐标系{li}与法向ni重合,那么
对于有m个定位元件的定位系统,采用m个等式的矩阵进行下述描述:
从而将式(17)写成:
式中:
Wpsl=NGl;
N=diag(n1...nm)∈R3m*m;
Δrpsl=(Δrn1...Δrnm)T∈Rm*1。
假设工件上第i个接触元件的几何位置Δrpsi足够小,则工件上夹具元件的位置误差可以忽略:
若基准偏差为小偏差,则式(14)的第二个等式表示为:
对于有m个定位点的工件,采用m个等式的矩阵进行下述描述:
从而将式(19)写成:
式中:
Wpsd=-NGl;
N=diag(n1...nm)∈R3m*m;
Δrpsd=(Δrn1...Δrnm)T∈Rm*1。
取广义逆后,则基准d偏差、定位l偏差与刀具c偏差在微元坐标系产生的误差分别为:
Δrpsc为等效夹具误差。
对于特征坐标系采用坐标变化得到:
则基准偏差、定位偏差与刀具偏差在特征坐标系产生的误差分别为:
则基准偏差、定位偏差与刀具偏差在工件坐标系产生的误差分别为:
式中,ΔΨli表示第i个定位销在全局坐标系下的位移量,表示微元坐标系相对于全局坐标系的旋转矩阵的名义值,Δrli表示工件与第i个定位销相接触定位销的偏差,特征坐标系相对于全局坐标系的旋转矩阵的名义值,Δrpsi表示工件与第i个定位销相接触的微元的偏差。
将式(13)、式(25)和式(29)带入到式(6)中,得:
式中:
Em=Eps=Ups;
Ews=Uws;
Efs=Ufs。
由式(13)得:
根据式(6)、式(7)、式(8)、式(33)和式(34)得到几何因素产生的误差。
工件坐标系下微元的几何偏差为:
特征坐标系下微元的几何偏差:
式中,
微元坐标系下的几何偏差:
式中:
步骤三、基于微元刚体坐标变换的受力变形误差计算。
式中,ki表示i点的刚度。
根据切削力实验求得切削力系数与切深之间的回归关系,得到不同切深下的切削力系数。在固定接触角和轴向切削深度下,改变进给速度和径向切深进行一组铣削实验,测量每个刀齿周期的平均力。一个周期内一个齿的三个方向的平均切削力如式(38)所示:
式中,N为刀齿数,fz为每齿进给量,da为轴向切深,KT,KA,KR为切削力系数θstart,θexit分别为切入切出角。
对不同径向切削深度下的切削力系数进行三次曲线拟合,得到切削力系数关于径向切削深度de的表达式为:
根据铣削力系数的表达式,得到不同径向切深下的铣削力系数,代入式(38),则得到不同径向切深与轴向切深下的工件的平均切削力。
工件表面误差Δ由两部分形成,工件沿曲面法向的变形Δw和刀具沿工件曲面法向的变形Δt,铣刀在轴向有较高的刚性,忽略刀具变形:
Δ=Δw (40)
当工件和刀具均为刚体时,轴向切削深度为da,当工件和刀具发生变形时,轴向切削深度变为d′a,从而得到:
da′=da+Δ轴向几何+Δ轴向让刀变形 (41)
径向切深变为:
d′e=de+Δ径向几何+Δ径向让刀变形 (42)
若忽略刀具变形对切入切出角的影响,将式(39)和式(41)代入式(38),得到工件发生变形时的平均切削力。
设刀轴在进给方向与曲面法向组成的XZ平面上的投影与Z轴正方向夹角为0,在切削刃与工件表面某微元的切触点P(u,v)处,由于有切削力的存在,工件将发生变形,工件的偏移量形成该微元处的尺寸误差。随着刀具的旋转,所产生的表面误差将沿着刀具进给的方向进行分布,且切触点处的误差一方面与切削力的大小有关,另一方面与切触点处工件的刚度K(u,v)相关。刀具在垂直于刀轴方向刚性较弱,在和作用下会发生弹性变形,但弹性变形产生的误差在随后的加工中被切除,不影响Z方向的加工表面误差。考虑对于工件变形的影响,工件在平面法向的受力为:
工件在平面法向的变形为:
则微元在工件坐标系下的偏差写成:
工件坐标系下切削力转换到微元坐标系下,则微元坐标系下切削力产生的变形为:
步骤四、基于扩展误差流的加工过程误差预测。
依据工序的输出值P和影响工序质量的误差源u建立扩展SoV多工序加工过程。具体的字母表示的含义如下:
假定误差为小误差,则依据状态空间得到扩展SoV流模型:
式中,Pt表示工序k加工质量特性波动的测量值或者预测值;μt表示工序k质量特性波动的真实值,为不可观测的状态值;wt表示工序k制造系统的状态噪声项,Wt表示其方差;vt表示工序k测量噪声项,Vt表示其方差;D0表示t=0时刻关于工序质量的初始信息集合;m0表示在D0条件下对工序质量的均值的一个估计值;C0表示关于均值m0的方差,是对m0一种不确定性的度量。并且认为,对所有的t和s,当t≠s时,vt和vs、wt和ws以及vt和ws都相互独立。
则多种误差源影响下的耦合误差表示为:
式中:
则得到扩展SoV误差耦合模型式(47)中的系数:
为了描述KPCs和KCCs之间的关系,将式(47)的状态方程带入到测量等式中,得到如下显式表达式,即为误差的预测模型:
本发明的有益效果是:该方法采用模型驱动进行过程监控,对几何定位误差源与受力变形误差源产生的误差及时进行预测,对误差耦合和传递机理进行了研究,对加工过程误差耦合建模方案进行了设计;通过网格划分的思想对叶片零件进行了微元划分,将变形分析转换为微元坐标系的坐标变换;对加工过程误差与误差源因素的映射关系进行了深入的研究分析,从加工工序流的角度采用扩展误差流的理论,建立了基于扩展误差流(SoV)的多源多工序误差耦合模型,在此基础上给出了基于误差耦合模型的综合误差预测模型。帮助工作人员及时发现叶片的制造缺陷,进行校正和调整,减少返工次数,提高了加工质量和效率。
下面结合附图和实施例对本发明作详细说明。
附图说明
图1是本发明方法的流程图。
图2是本发明方法中的多种误差源驱动下的误差建模方案示意图。
图3是本发明方法中工艺系统示意图。
图4是本发明方法中切削过程切深变化的示意图。
图5是本发明方法中工件表面误差形成的示意图。
图6是本发明方法中叶片零件多工序加工过程误差流表示方法的示意图。
图7是本发明方法中工件尺寸和定位销布局图。
具体实施方式
参照图1-7。以某航空发动机制造厂加工某种叶片(材料钛合金TC4)零件为例,依据步骤一对叶片加工过程进行误差分析预测。对本发明进行详细描述,验证本发明对某类型叶片多工序加工过程误差的预测效果。
对叶片零件的截面线进行等精度离散,可以得到微元坐标点,并将叶片零件简化为薄板零件,将对应的微元坐标点也映射到叶片简化模型,可以得到薄板上与叶片零件等进度离散处理后点相对应,用薄板上的这些对应点近似分析叶片上的对应点。
在简化模型中,ps1到ps4为微元。考虑了几种可能产生误差源的情况,假定过程由两道工序构成,第一道工序加工部位1,第二道工序加工部位2,第一道工序的加工特征是第二道工序的加工基准。
在简化模型上构建全局坐标系、工件坐标系、刀具坐标系、加工特征坐标系以及微元坐标系,为简化计算,让全局坐标系和工件坐标重合,加工特征面的中心为特征坐标系。
1、确定叶片加工过程误差源。
本实施例考虑影响叶片零件加工的关键几何定位误差与受力变形误差,分别为:(1)定位基准精度产生的误差;(2)安装产生的误差(夹具)误差;(3)刀具位姿变化产生的误差;(4)切削力产生的变形误差;(5)噪声误差。
本实施例对以上引起误差的主要误差源进行建模,并对加工过程中的因素做如下的约定:
(1)多工序加工过程中将一道工序的加工过程分为安装好(加工前)和加工中两个阶段。
(2)基准误差、夹具安装几何误差以及刀具倾斜误差对工件产生定位误差,定义为几何定位误差,它是加工前引起的叶片零件刚体误差的主要原因。
(3)加工中的弹性变形主要是由切削力引起的。
2、基于微元刚体坐标变换的几何定位误差计算。
为了描述一个叶片零件的加工过程,定义了5个坐标系,分别为全局坐标系{O}、工件坐标系{ws},特征坐标系{fs},和微元坐标系{ps}以及刀具坐标系{T}。
根据各个对象在空间的位姿关系,由全局坐标系{O}、工件坐标系{ws},特征坐标系{Fs},和微元坐标系{ps}以及刀具坐标系共同构成空间运动链可知:
表示刀具坐标系{T}相对于全局坐标系{O}的坐标变换矩阵,表示工件坐标系{ws}相对于全局坐标系{O}的坐标变换矩阵,表示特征坐标系{fs}相对于工件坐标系{ws}的坐标变换矩阵,表示微元坐标系{ps}相对于特征坐标系{fs}的坐标变换矩阵表示刀具坐标系{T}相对于微元坐标系{ps}的坐标变换矩阵。
运用微分运动原理得到:
式中,Δ表示微分。
式(2)的左边变换为
则得到:
将式(3)改写为:
因此得到以下表达式(忽略高阶项):
从而得到刀具在微元坐标系下的偏差:
由此得到刀具在特征坐标系下和在工件坐标系下的偏差分别为:
式中,表示微元坐标系相对于全局坐标系的旋转矩阵的名义值,表示特征坐标系相对于全局坐标系的旋转矩阵的名义值,表示工件坐标系相对于全局坐标系的旋转矩阵的名义值, 和表示反对称矩阵,为刀具相对于全局坐标系的偏差量,为工件相对于全局坐标系的偏差量,为特征相对于全局坐标系的偏差量,表示微元相对于全局坐标系的偏差量。
在装夹好工件时,理想情况下接触点在工件坐标系、特征坐标系以及微元坐标系与定位元件坐标系是重合的则工件与元件l的第i个接触点有四种方法描述:
则该式中的前三个等式与最后一个相等,即Fi=fi
采用一阶泰勒级数展开,省去高阶项得到:
将式(10)写成:
则得到:
由于夹具固定在机床上,假定忽略掉夹具的位置和姿态误差,则得到:
从而得到:
假设工件上第i个接触点的几何位置Δrpsi足够小,则忽略工件上基准面接触点的位置误差:
假定在i个接触点在法向上有误差Δrni,并且坐标系{li}与法向ni重合,那么
对于有m个定位元件的定位系统,采用m个等式的矩阵进行下述描述:
从而将式(17)写成:
式中:
Wpsl=NGl;
N=diag(n1...nm)∈R3m*m;
Δrpsl=(Δrn1...Δrnm)T∈Rm*1。
假设工件上第i个接触元件的几何位置△rpsi足够小,则工件上夹具元件的位置误差可以忽略:
若基准偏差为小偏差,则式(14)的第二个等式表示为:
对于有m个定位点的工件,采用m个等式的矩阵进行下述描述:
从而将式(19)写成:
式中:
Wpsd=-NGl;
N=diag(n1...nm)∈R3m*m;
△rpsd=(△rn1...△rnm)T∈Rm*1。
取广义逆后,则基准d偏差、定位l偏差与刀具c偏差在微元坐标系产生的误差分别为:
△rpsc为等效夹具误差。
对于特征坐标系采用坐标变化得到:
则基准偏差、定位偏差与刀具偏差在特征坐标系产生的误差分别为:
则基准偏差、定位偏差与刀具偏差在工件坐标系产生的误差分别为:
式中,△Ψli表示第i个定位销在全局坐标系下的位移量,表示微元坐标系相对于全局坐标系的旋转矩阵的名义值,△rli表示工件与第i个定位销相接触定位销的偏差,特征坐标系相对于全局坐标系的旋转矩阵的名义值,△rpsi表示工件与第i个定位销相接触的微元的偏差。
将式(13)、式(25)和式(29)带入到式(6)中,得:
式中:
Em=Eps=Ups;
Ews=Uws;
Efs=Ufs。
由式(13)得:
据以上分析,根据式(6)、式(7)、式(8)、式(33)和式(34),从而得到几何因素产生的误差。
工件坐标系下微元的几何偏差为:
特征坐标系下微元的几何偏差:
式中,
微元坐标系下的几何偏差:
式中:
假定简化模型的两道工序中工件的微元及其对应的特征的关键质量特性如表1所示,并分别给出了其名义方向和位置。
表1微元的名义方向和位置
两道工序各自的定位元件位置和方向向量如表2所示。
表2两道工序各自的定位元件的位置和方向向量
为了方便计算,两道工序的坐标系和定位元件信息相同,则根据各坐标系之间的相对位置关系得到转换矩阵相同,为:
根据表2得到以下定位相关的矩阵:
定位销的姿态矩阵相同,为:
式中,ni=[0 0 1](i=1~3),ni=[0 1 0](i=4~5),n6=[-1 0 0]。
定位矩阵相同,为:
式中,
第一道工序几何定位误差预测。
a.基准引起的误差。
基准误差是上一道工序引起的误差,假定第一道工序要加工的特征的前几道工序没有偏差,则将其的基准误差设置为0。
b.夹具几何偏差造成的工件加工误差。
对于定位点1,夹具元件1产生的误差为Rf=[0,0.25,0,0,0,0]T,则由式(35) 计算得到第一道工序的第一个微元ps1的在工件坐标系的误差为:
[0.0024,0.0269,0.00383,-0.0008,0.0002,-0.0000]T
由式(37) 计算在微元坐标系的误差为[-0.0026,-0.0452,0.00124,0.0008,0.0000,0.0002]T。
同理得第二个微元ps2在工件坐标系误差为[0.0024,0.00109,0.00383,-0.0008,0.0002,-0.0000]T。
在特征坐标系的误差为[-0.0027,-0.00303,0.00112,0.0008,0.0000,0.0002]T。
在微元坐标系的误差为[-0.0028,-0.00276,0.00113,0.0008,0.0000,0.0002]T。
同理得第三个微元ps3在工件坐标系误差为[0.0024,0.00109,0.0383,-0.0008,0.0002,-0.0000]T。
在特征坐标系的误差为[-0.0027,-0.00301,0.00112,0.0008,0.0000,0.0002]T。
在微元坐标系的误差为[-0.0028,-0.00466,0.00134,0.0008,0.0000,0.0002]T。
可以看出微元ps1受到的影响最大,而微元ps2和ps3的误差经过计算近似为0。
c.切削力引起的工件的误差。
因为是刚体,切削力引起的误差假定为0。
d.刀具位姿引起的误差。
刀具位姿误差设置为0。
设置噪声为0,则通过扩展SoV预测模型公式 (50) 可以得到微元ps1的预测误差值为:
[-0.0026,-0.0452,0.00124,0.0008,0.0000,0.0002]T
第二道工序几何定位误差预测
由第一道工序可以知道第二道工序基准ps1产生的微元偏差在微元坐标系下为:
[-0.0026,-0.0452,0.00124,0.0008,0.0000,0.0002]T,基准ps2和ps3的偏差近似为[0,0,0,0,0,0]。
第二道工序定位基准1在工序二的工件坐标系下产生的偏差为:
式中,
式中,
3、基于微元刚体坐标变换的受力变形误差计算。
受力变形是耦合误差,是叶片零件加工质量的重要影响因素,加工中切削参数的变化会导致它发生变化。本实施例只考虑静态切削力引起的静态变形,一般采用未变形切削力平均值的方法进行分析。前面所描述的几何定位误差源基准、夹具元件以及刀具位姿等产生的几何定位误差影响切削深度使得切削力发生变化,当进行切削时,会产生让刀变形,影响零件最终的加工质量。本实施例通过辨识切削力系数,建立切削力系数与切削参数的回归关系,求解平均铣削力,然后采用UG高级仿真求解叶片微元的刚度,最后根据胡克定律求解受力变形。
(1)加工过程中切深的变化。
在加工前,几何定位误差影响了工件的位置,使切深发生了变化;在加工过程中,切深是不断变化的,切深的变化影响着切削力,切削力又会影响切深。因此有必要分析加工过程切深的变化以及它对最终质量的影响。
式中,ki表示i点的刚度。
由于切深是不断变化的,通过切深可以求得切削力,从而得到由切削力产生的让刀变形。为此,需要进行切削力的求解。本实施例采用平底刀端铣的方式,进行铣削力系数的标定,进而求取平均铣削力,得到平底刀的微元划分、参数以及切削力的辨识,得到随径向切深变化的切削力系数,进而得到不同径向切深下的平均铣削力,从而计算平均铣削力,进而求得对应的铣削力变形。
(2)铣削加工的平均切削力计算。
根据切削力实验求得切削力系数与切深之间的回归关系,得到不同切深下的切削力系数。在固定接触角和轴向切削深度下,改变进给速度和径向切深进行一组铣削实验,测量每个刀齿周期的平均力。实验获得的平均切削力与从切削力表达式中求得的平均切削力相等,用于辨识切削力系数。由于一个齿周期内每个刀齿切除的材料总量是一个常数,与螺旋角无关,因此平均切削力与螺旋角无关。一个周期内一个齿的三个方向的平均切削力如式(38)所示:
式中,N为刀齿数,fz为每齿进给量,da为轴向切深,KT,KA,KR为切削力系数θstart,θexit分别为切入切出角。
铣槽铣削实验是最为方便快捷的切削力系数标定方法,此时,刀具的切入角θstart和切出角θexit分别为0和π。主轴转速为2000rpm,进给速度为320mm/min,轴向切深为2mm,刀具螺旋角为35°,刀具半径为5mm,刀齿数为4时,设置径向切削深度分别为1mm、2mm、3mm、4mm、5mm、6mm、7mm、8mm、9mm、10mm时,通过切削力实验计算来标定得到平底刀径向切削力系数。在本组试验中,选用的铣刀为四齿硬质合金平底铣刀,直径是10mm,采用平底刀端铣方式,因此径向切深取值从1mm增至10mm,增量为1mm。试验所用机床为JOHNFORD VMC-850四轴立式数控加工中心,工件材料选用钛合金TC4。
对不同径向切削深度下的切削力系数进行三次曲线拟合,得到切削力系数关于径向切削深度de的表达式为:
根据铣削力系数的表达式,得到不同径向切深下的铣削力系数,代入式(38),则得到不同径向切深与轴向切深下的工件的平均切削力。
(3)基于平均切削力的微元加工变形计算。
叶片零件的加工尺寸和形状误差影响因素比较多,如刀具工件接触区的热变形,工件残余应力变形等。本实施例只考虑由切削力引起的弹性变形,具体指在刀具沿着预先规划的路径运动时,由于切削力的变化以及刀具和工件之间的相对刚度变化,使得工件尺寸偏离了期望值的这种相对位移引起的尺寸和形状误差。在加工过程中,切削力使得刀具和叶片零件产生弹性变形,走刀过后弹性变形恢复,致使部分材料未被切除,造成零件表面的加工误差。加工面法向的分力是决定工件表面误差的主要因素,所以本实施例在计算叶片零件变形时将以法向分力为主。
工件表面误差Δ由两部分形成,工件沿曲面法向的变形Δw和刀具沿工件曲面法向的变形Δt,通常铣刀在轴向有较高的刚性,忽略刀具变形:
Δ=Δw (40)
当工件和刀具均为刚体时,轴向切削深度为da,当工件和刀具发生变形时,轴向切削深度变为da′,从而得到:
da′=da+Δ轴向几何+Δ轴向让刀变形 (41)
径向切深变为:
de′=de+Δ径向几何+Δ径向让刀变形 (42)
若忽略刀具变形对切入切出角的影响,将式(39)和式(41)代入式(38),得到工件发生变形时的平均切削力。
设刀轴在进给方向与曲面法向组成的平面(XZ平面)上的投影与曲面法向(Z轴正方向)夹角为0,在切削刃与工件表面某微元的切触点P(u,v)处,由于有切削力的存在,工件将发生变形,工件的偏移量形成该微元处的尺寸误差。随着刀具的旋转,所产生的表面误差将沿着刀具进给的方向进行分布,且切触点处的误差一方面与切削力的大小有关,另一方面与切触点处工件的刚度K(u,v)相关。刀具在垂直于刀轴方向刚性较弱,在和作用下会发生弹性变形,但弹性变形产生的误差在随后的加工中被切除,不影响Z方向的加工表面误差,因此本实施例忽略刀具在这两个方向的变形。因此,只考虑对于工件变形的影响,工件在平面法向的受力为:
工件在平面法向的变形为:
则微元在工件坐标系下的偏差写成:
工件坐标系下切削力转换到微元坐标系下,则微元坐标系下切削力产生的变形为:
设定第一道工序的切削力引起的误差为0,则只需计算第二道工序的切削力产生的误差。刀具半径为5mm,径向切深为2mm,轴向切深与标定系数时保持一致为2mm,采用逆铣方式,得到切入切出角分别为0和Arccos0.6。由于定位产生偏差,使得轴向切深发生变化:
径向切深变为1.933mm,则由式(39) 计算出切削力系数,由式(38)
计算微元ps4的雅克比矩阵:
得到
又由于
从而得到微元ps4总偏差为:
4、基于扩展误差流的加工过程误差预测。
叶片零件多源多工序加工过程误差耦合是指依据工件的多工序过程,由于误差源的耦合影响作用使得误差在本道工序积累和向下道工序传递的过程。依据工序的输出值P和影响工序质量的误差源u建立扩展SoV多工序加工过程。具体的字母表示的含义如下:
假定误差为小误差,则依据状态空间得到扩展SoV流模型:
式中,Pt表示工序k加工质量特性波动的测量值或者预测值;μt表示工序k质量特性波动的真实值,为不可观测的状态值;wt表示工序k制造系统的状态噪声项(或状态误差项),Wt表示其方差;vt表示工序k测量噪声项(或测量误差项),Vt表示其方差;D0表示t=0时刻关于工序质量的初始信息集合;m0表示在D0条件下对工序质量的均值的一个估计值;C0表示关于均值m0的方差,是对m0一种不确定性的度量。并且认为,对所有的t和s,当t≠s时,vt和vs、wt和ws以及vt和ws都相互独立。假设测量噪声和状态噪声服从高斯核拉普拉斯分布。在该波动模型中,观测方程反映了工序k对制造过程质量特征的观测状况,状态方程反映了工序k制造过程的质量波动状况。其中系数矩阵Ak-1,Bk和是需要求解和计算的,根据上道工序、本道输入误差源和测量系统进行确定。
则多种误差源影响下的耦合误差表示为:
式中:
则得到扩展SoV误差耦合模型式(47)中的系数:
当工序被加工后,则该道工序的输出值可以测量得到,可以估计出状态值来。但是当工序还没有被加工,测量值无法知道,为了预测该工序的加工质量,需要将扩展SoV模型变换成预测模型,对该道零件的加工质量进行预测。
为了描述KPCs和KCCs之间的关系,将式(47)的状态方程带入到测量等式中,得到如下显式表达式,即为误差的预测模型:
式中,是追踪基准误差、夹具几何误差、刀具位姿误差以及切削力误差的状态转移矩阵,i=1,...k-1,并且初始状态向量μ0表示的是一个零件的KQCs在进入第一道工序加工前的初始偏差。这些初始偏差信息来自于过去加工过程的历史实验。
设置夹具元件的误差为0,刀具位姿误差为0,噪声为0,则通过扩展SoV预测得到最终微元ps4在工件坐标系下的误差值为:
[-0.0005,-0.0516,0.1274,0.0002,0,0,]T
可以看出由于几何定位误差与切深引起的切削力的变化,从而造成了y方向和z方向产生了偏差。通过本实施例得到误差值为[-0.0003,-0.0520,0.01563,0.0001,0.0001,0.0001]T,与预测值很接近,符合所设定的误差源的误差产生规律,因此预测是合理的。
Claims (1)
1.一种基于扩展误差流的叶片加工过程误差预测方法,其特征在于包括以下步骤:
步骤一、确定叶片加工过程误差源分别为定位基准精度产生的误差;安装产生的误差;刀具位姿变化产生的误差;切削力产生的变形误差;噪声误差;对误差源进行建模,并对加工过程中的因素约定如下:多工序加工过程中将一道工序的加工过程分为加工前和加工中两个阶段;基准误差、夹具安装几何误差以及刀具倾斜误差对工件产生定位误差,定义为几何定位误差;加工中的弹性变形是由切削力引起的;
步骤二、基于微元刚体坐标变换的几何定位误差计算;
定义五个坐标系,分别是全局坐标系{O}、工件坐标系{ws}、特征坐标系{fs}、微元坐标系{ps}以及刀具坐标系{T};
根据各个对象在空间的位姿关系,由全局坐标系{O}、工件坐标系{ws}、特征坐标系{Fs}、微元坐标系{ps}以及刀具坐标系共同构成空间运动链:
式中,表示刀具坐标系{T}相对于全局坐标系{O}的坐标变换矩阵,表示工件坐标系{ws}相对于全局坐标系{O}的坐标变换矩阵,表示特征坐标系{fs}相对于工件坐标系{ws}的坐标变换矩阵,表示微元坐标系{ps}相对于特征坐标系{fs}的坐标变换矩阵表示刀具坐标系{T}相对于微元坐标系{ps}的坐标变换矩阵;
运用微分运动原理得到:
式中,Δ表示微分;
式(2)的左边变换为
则得到:
将式(3)改写为:
忽略高阶项后得到以下表达式:
(5)
从而得到刀具在微元坐标系下的偏差:
由此得到刀具在特征坐标系下和在工件坐标系下的偏差分别为:
式中,表示微元坐标系相对于全局坐标系的旋转矩阵的名义值,表示特征坐标系相对于全局坐标系的旋转矩阵的名义值,表示工件坐标系相对于全局坐标系的旋转矩阵的名义值, 和表示反对称矩阵,为刀具相对于全局坐标系的偏差量,为工件相对于全局坐标系的偏差量,为特征相对于全局坐标系的偏差量,表示微元相对于全局坐标系的偏差量;
在装夹好工件时,接触点在工件坐标系{ws}、特征坐标系{fs}、微元坐标系{ps}与定位元件坐标系是重合的则工件与元件l的第i个接触点有四种方法描述:
则该式中的前三个等式与最后一个相等,即Fi=fi
采用一阶泰勒级数展开,省去高阶项得到:
将式(10)写成:
则得到:
由于夹具固定在机床上,假定忽略掉夹具的位置和姿态误差,则得到:
从而得到:
假设工件上第i个接触点的几何位置Δrpsi足够小,则忽略工件上基准面接触点的位置误差:
假定在i个接触点在法向上有误差Δrni,并且坐标系{li}与法向ni重合,那么
对于有m个定位元件的定位系统,采用m个等式的矩阵进行下述描述:
从而将式(17)写成:
式中:
Wpsl=NGl;
N=diag(n1...nm)∈R3m*m;
Δrpsl=(Δrn1...Δrnm)T∈Rm*1;
假设工件上第i个接触元件的几何位置Δrpsi足够小,则工件上夹具元件的位置误差可以忽略:
若基准偏差为小偏差,则式(14)的第二个等式表示为:
对于有m个定位点的工件,采用m个等式的矩阵进行下述描述:
从而将式(19)写成:
式中:
Wpsd=-NGl;
N=diag(n1...nm)∈R3m*m;
Δrpsd=(Δrn1...Δrnm)T∈Rm*1;
取广义逆后,则基准d偏差、定位l偏差与刀具c偏差在微元坐标系产生的误差分别为:
Δrpsc为等效夹具误差;
对于特征坐标系采用坐标变化得到:
则基准偏差、定位偏差与刀具偏差在特征坐标系产生的误差分别为:
则基准偏差、定位偏差与刀具偏差在工件坐标系产生的误差分别为:
式中,ΔΨli表示第i个定位销在全局坐标系下的位移量,表示微元坐标系相对于全局坐标系的旋转矩阵的名义值,Δrli表示工件与第i个定位销相接触定位销的偏差,特征坐标系相对于全局坐标系的旋转矩阵的名义值,Δrpsi表示工件与第i个定位销相接触的微元的偏差;
将式(13)、式(25)和式(29)带入到式(6)中,得:
式中:
Em=Eps=Ups;
Ews=Uws;
Efs=Ufs;
由式(13)得:
根据式(6)、式(7)、式(8)、式(33)和式(34)得到几何因素产生的误差;
工件坐标系下微元的几何偏差为:
特征坐标系下微元的几何偏差:
式中,
微元坐标系下的几何偏差:
式中:
步骤三、基于微元刚体坐标变换的受力变形误差计算;
式中,ki表示i点的刚度;
根据切削力实验求得切削力系数与切深之间的回归关系,得到不同切深下的切削力系数;在固定接触角和轴向切削深度下,改变进给速度和径向切深进行一组铣削实验,测量每个刀齿周期的平均力;一个周期内一个齿的三个方向的平均切削力如式(38)所示:
式中,N为刀齿数,fz为每齿进给量,da为轴向切深,KT,KA,KR为切削力系数θstart,θexit分别为切入切出角;
对不同径向切削深度下的切削力系数进行三次曲线拟合,得到切削力系数关于径向切削深度de的表达式为:
根据铣削力系数的表达式,得到不同径向切深下的铣削力系数,代入式(38),则得到不同径向切深与轴向切深下的工件的平均切削力;
工件表面误差Δ由两部分形成,工件沿曲面法向的变形Δw和刀具沿工件曲面法向的变形Δt,铣刀在轴向有较高的刚性,忽略刀具变形:
Δ=Δw (40)
当工件和刀具均为刚体时,轴向切削深度为da,当工件和刀具发生变形时,轴向切削深度变为d′a,从而得到:
da′=da+Δ轴向几何+Δ轴向让刀变形 (41)
径向切深变为:
d′e=de+Δ径向几何+Δ径向让刀变形 (42)
若忽略刀具变形对切入切出角的影响,将式(39)和式(41)代入式(38),得到工件发生变形时的平均切削力;
设刀轴在进给方向与曲面法向组成的XZ平面上的投影与Z轴正方向夹角为0,在切削刃与工件表面某微元的切触点P(u,v)处,由于有切削力的存在,工件将发生变形,工件的偏移量形成该微元处的尺寸误差;随着刀具的旋转,所产生的表面误差将沿着刀具进给的方向进行分布,且切触点处的误差一方面与切削力的大小有关,另一方面与切触点处工件的刚度K(u,v)相关;刀具在垂直于刀轴方向刚性较弱,在和作用下会发生弹性变形,但弹性变形产生的误差在随后的加工中被切除,不影响Z方向的加工表面误差;考虑对于工件变形的影响,工件在平面法向的受力为:
工件在平面法向的变形为:
则微元在工件坐标系下的偏差写成:
工件坐标系下切削力转换到微元坐标系下,则微元坐标系下切削力产生的变形为:
步骤四、基于扩展误差流的加工过程误差预测;
依据工序的输出值P和影响工序质量的误差源u建立扩展SoV多工序加工过程;具体的字母表示的含义如下:
假定误差为小误差,则依据状态空间得到扩展SoV流模型:
式中,Pt表示工序k加工质量特性波动的测量值或者预测值;μt表示工序k质量特性波动的真实值,为不可观测的状态值;wt表示工序k制造系统的状态噪声项,Wt表示其方差;vt表示工序k测量噪声项,Vt表示其方差;D0表示t=0时刻关于工序质量的初始信息集合;m0表示在D0条件下对工序质量的均值的一个估计值;C0表示关于均值m0的方差,是对m0一种不确定性的度量;并且认为,对所有的t和s,当t≠s时,vt和vs、wt和ws以及vt和ws都相互独立;
则多种误差源影响下的耦合误差表示为:
ΔP=ΔPI+ΔPII
式中:
则得到扩展SoV误差耦合模型式(47)中的系数:
为了描述KPCs和KCCs之间的关系,将式(47)的状态方程带入到测量等式中,得到如下显式表达式,即为误差的预测模型:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2013103767650A CN103440416A (zh) | 2013-08-27 | 2013-08-27 | 基于扩展误差流的叶片加工过程误差预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2013103767650A CN103440416A (zh) | 2013-08-27 | 2013-08-27 | 基于扩展误差流的叶片加工过程误差预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103440416A true CN103440416A (zh) | 2013-12-11 |
Family
ID=49694109
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2013103767650A Pending CN103440416A (zh) | 2013-08-27 | 2013-08-27 | 基于扩展误差流的叶片加工过程误差预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103440416A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104268350A (zh) * | 2014-09-30 | 2015-01-07 | 中国西电电气股份有限公司 | 一种仿真预测和实际生产相集成的闭环质量控制仿真方法 |
CN104298171A (zh) * | 2014-09-24 | 2015-01-21 | 北京理工大学 | 一种介观尺度零件加工误差预测以及控制方法 |
CN106709193A (zh) * | 2016-12-30 | 2017-05-24 | 西北工业大学 | 基于学习算法的航空发动机薄壁叶片加工误差补偿方法 |
CN107111306A (zh) * | 2014-12-05 | 2017-08-29 | 赛峰飞机发动机公司 | 基于定心系数的分析的零件制造方法 |
CN107729287A (zh) * | 2017-09-14 | 2018-02-23 | 广东工业大学 | 一种用于柔性材料卷对卷加工过程张力预测的sov方法 |
CN112272803A (zh) * | 2018-06-08 | 2021-01-26 | 赫克斯冈技术中心 | 利用刀具路径修正制造多个零件的方法和设备 |
CN113076977A (zh) * | 2021-03-22 | 2021-07-06 | 华中科技大学 | 一种多工序加工中的误差溯源方法及其应用 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5822208A (en) * | 1994-04-12 | 1998-10-13 | Bay Instrumentation & Technology Co. | Method and apparatus for predicting errors in a manufacturing process |
CN1364675A (zh) * | 2002-02-05 | 2002-08-21 | 华中科技大学 | 镗削尺寸误差预测补偿控制方法及其装置 |
US20050193284A1 (en) * | 2004-02-06 | 2005-09-01 | Fujitsu Limited | Electronic device, failure prediction method, and computer product |
CN101436057A (zh) * | 2008-12-18 | 2009-05-20 | 浙江大学 | 数控机床热误差贝叶斯网络补偿方法 |
CN103197598A (zh) * | 2013-03-13 | 2013-07-10 | 西北工业大学 | 基于非参数多元统计过程控制的叶身型面加工误差控制方法 |
-
2013
- 2013-08-27 CN CN2013103767650A patent/CN103440416A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5822208A (en) * | 1994-04-12 | 1998-10-13 | Bay Instrumentation & Technology Co. | Method and apparatus for predicting errors in a manufacturing process |
CN1364675A (zh) * | 2002-02-05 | 2002-08-21 | 华中科技大学 | 镗削尺寸误差预测补偿控制方法及其装置 |
US20050193284A1 (en) * | 2004-02-06 | 2005-09-01 | Fujitsu Limited | Electronic device, failure prediction method, and computer product |
CN101436057A (zh) * | 2008-12-18 | 2009-05-20 | 浙江大学 | 数控机床热误差贝叶斯网络补偿方法 |
CN103197598A (zh) * | 2013-03-13 | 2013-07-10 | 西北工业大学 | 基于非参数多元统计过程控制的叶身型面加工误差控制方法 |
Non-Patent Citations (3)
Title |
---|
张智海等: ""基于铣削力/力矩模型的铣削表面几何误差模型"", 《机械工程学报》 * |
李昊: ""航空发动机叶片铣削过程变形控制研究"", 《中国优秀硕士学位论文全文数据库 工程科技I辑》 * |
王佩等: ""基于动态Bayesian网络的叶片加工质量监控与溯源"", 《航空学报》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104298171A (zh) * | 2014-09-24 | 2015-01-21 | 北京理工大学 | 一种介观尺度零件加工误差预测以及控制方法 |
CN104298171B (zh) * | 2014-09-24 | 2017-02-15 | 北京理工大学 | 一种介观尺度零件加工误差预测以及控制方法 |
CN104268350A (zh) * | 2014-09-30 | 2015-01-07 | 中国西电电气股份有限公司 | 一种仿真预测和实际生产相集成的闭环质量控制仿真方法 |
CN107111306A (zh) * | 2014-12-05 | 2017-08-29 | 赛峰飞机发动机公司 | 基于定心系数的分析的零件制造方法 |
CN107111306B (zh) * | 2014-12-05 | 2019-05-14 | 赛峰飞机发动机公司 | 基于定心系数的分析的零件制造方法 |
CN106709193A (zh) * | 2016-12-30 | 2017-05-24 | 西北工业大学 | 基于学习算法的航空发动机薄壁叶片加工误差补偿方法 |
CN106709193B (zh) * | 2016-12-30 | 2019-07-19 | 西北工业大学 | 基于学习算法的航空发动机薄壁叶片加工误差补偿方法 |
CN107729287A (zh) * | 2017-09-14 | 2018-02-23 | 广东工业大学 | 一种用于柔性材料卷对卷加工过程张力预测的sov方法 |
CN107729287B (zh) * | 2017-09-14 | 2021-05-18 | 广东工业大学 | 一种用于柔性材料卷对卷加工过程张力预测的sov方法 |
CN112272803A (zh) * | 2018-06-08 | 2021-01-26 | 赫克斯冈技术中心 | 利用刀具路径修正制造多个零件的方法和设备 |
CN113076977A (zh) * | 2021-03-22 | 2021-07-06 | 华中科技大学 | 一种多工序加工中的误差溯源方法及其应用 |
CN113076977B (zh) * | 2021-03-22 | 2022-04-19 | 华中科技大学 | 一种多工序加工中的误差溯源方法、设备及存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103440416A (zh) | 基于扩展误差流的叶片加工过程误差预测方法 | |
CN102591257B (zh) | 面向参数曲线刀具轨迹的数控系统轮廓误差控制方法 | |
CN104289748B (zh) | 一种大型薄壁蒙皮自适应等壁厚铣削系统及其加工方法 | |
CN105069257A (zh) | 一种自由曲面微细铣削切削力建模方法 | |
CN103218475B (zh) | 一种基于在机测评的复杂空间型面误差反馈补偿方法 | |
CN103777570B (zh) | 基于nurbs曲面的加工误差快速检测补偿方法 | |
CN104759942A (zh) | 一种薄壁零件的铣削变形在线测量与补偿加工方法 | |
CN106709193B (zh) | 基于学习算法的航空发动机薄壁叶片加工误差补偿方法 | |
CN105234743A (zh) | 一种五轴加工中心刀具变形误差补偿方法 | |
CN102880114B (zh) | 飞机结构件数控加工在线自适应刀轨补偿方法 | |
CN105222712A (zh) | 一种改进的大曲率半径曲面零件法向测量与调整方法 | |
CN108544181B (zh) | 一种整体叶盘损伤叶片的修复方法 | |
CN101870073B (zh) | 基于工艺系统刚度特性的多轴数控加工刀具运动规划方法 | |
CN105205221A (zh) | 一种重型数控机床精度可靠性分析方法 | |
CN105108215A (zh) | 一种自由曲面微细铣削让刀误差预测及补偿方法 | |
CN105426697A (zh) | 一种螺旋铣刀五轴加工铣削力精确预测方法 | |
CN102855362B (zh) | 一种基于计算机辅助设计的检验管路装配的方法 | |
Wang et al. | FEM-based simulation for workpiece deformation in thin-wall milling | |
Guo et al. | Cutting force modeling for non-uniform helix tools based on compensated chip thickness in five-axis flank milling process | |
CN104537172A (zh) | 一种通过优化夹具布局控制孔组复合位置度误差的方法 | |
Hou et al. | The variable radial depth of cut in finishing machining of thin-walled blade based on the stable-state deformation field | |
Mickaël et al. | 3D ISO manufacturing specifications with vectorial representation of tolerance zones | |
CN103438844A (zh) | 基于粒子群算法的复杂曲面工件定位方法 | |
Zhao et al. | Measured data-driven shape-adaptive machining via spatial deformation of tool cutter positions | |
Chatelain | A level-based optimization algorithm for complex part localization |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20131211 |