CN102799143A - 薄壁曲面零件加工质量缺陷区域计算方法 - Google Patents

薄壁曲面零件加工质量缺陷区域计算方法 Download PDF

Info

Publication number
CN102799143A
CN102799143A CN2012103038228A CN201210303822A CN102799143A CN 102799143 A CN102799143 A CN 102799143A CN 2012103038228 A CN2012103038228 A CN 2012103038228A CN 201210303822 A CN201210303822 A CN 201210303822A CN 102799143 A CN102799143 A CN 102799143A
Authority
CN
China
Prior art keywords
thin
surface part
wall curved
lathe
matrix
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
CN2012103038228A
Other languages
English (en)
Other versions
CN102799143B (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201210303822.8A priority Critical patent/CN102799143B/zh
Publication of CN102799143A publication Critical patent/CN102799143A/zh
Application granted granted Critical
Publication of CN102799143B publication Critical patent/CN102799143B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

本发明涉及薄壁曲面零件加工表面质量缺陷快速计算方法。包括如下步骤:1.产生机床铣削薄壁曲面零件的各坐标轴数控指令;2.计算产生机床加工薄壁曲面零件的刚度分布:根据步骤1中得到的机床各坐标轴数控指令,计算机床铣削薄壁曲面零件过程中各坐标轴合成的刚度矩阵,将刚度矩阵沿零件法线方向投影,得到加工过程中机床的刚度分布;3.计算产生薄壁曲面零件在加工过程中自身的刚度分布;4.计算薄壁曲面零件的加工质量缺陷区域。本发明的有益效果:通过机床和薄壁曲面零件的法向刚度分布计算,产生加工过程中的机床和零件自身的刚度分布综合对比图,从而得到薄壁曲面零件加工质量缺陷区域,针对加工质量缺陷区域设定不同的加工工艺。

Description

薄壁曲面零件加工质量缺陷区域计算方法
技术领域
本发明属于薄壁曲面零件的制造技术领域,尤其涉及薄壁曲面零件的五坐标数控机床铣削加工技术领域。
背景技术
五轴联动数控机床是指在一台机床上有五个坐标轴,包括三个平动坐标轴和两个旋转坐标轴,机床的各坐标轴在计算机数控系统的控制下按一定的速度同时到达某一个设定的点,协调运动进行加工,特别适用用于加工复杂曲面零件。通常薄壁曲面零件的加工质量控制最为困难,一方面数控机床坐标轴间的不同配合状态可能引起机床加工时刚度的实时变化,另一方面零件在切削力的作用下也可能引起零件的加工变形。机床和零件两方面因素叠加使得薄壁曲面零件的表面加工质量成为难点,极易引起曲面局部的表面质量缺陷,出现粗糙度、波纹度下降以及切削折痕,导致零件质量不合格,造成成本的浪费和效率的损失。现有方法中缺少对薄壁曲面零件加工前表面质量的缺陷预测,如果能在加工零件之前根据零件特性和机床运动状态预测出零件表面质量缺陷区域,就可以在工艺实施中对于缺陷区域设定不同的加工工艺方案,保障零件的加工质量和提高零件加工的成品率。
由于企业对薄壁曲面零件质量控制的需求,基于薄壁曲面零件表面质量缺陷区域的快速计算成为一技术要点。由此可在切削前预知零件加工表面质量缺陷区域,针对表面质量缺陷区域设定不同的加工工艺,从而为零件的工艺优化和精度保障提供依据。
发明内容
本发明的目的是为了快速计算五轴联动数控机床(以下简称机床)在切削零件时的加工质量缺陷区域从而提高机床加工薄壁曲面零件的表面质量,提出了一种薄壁曲面零件加工表面质量缺陷快速计算方法。
本发明的技术方案:薄壁曲面零件加工表面质量缺陷快速计算方法,其特征在于,包括如下步骤:
步骤1.产生机床铣削薄壁曲面零件的各坐标轴数控指令;
步骤2:计算产生机床加工薄壁曲面零件的刚度分布:根据步骤1中得到的机床各坐标轴数控指令,计算机床铣削薄壁曲面零件过程中各坐标轴合成的刚度矩阵,将刚度矩阵沿零件法线方向投影,得到加工过程中机床的刚度分布;
步骤3:计算产生薄壁曲面零件在加工过程中自身的刚度分布;
步骤4:计算薄壁曲面零件的加工质量缺陷区域:将步骤2中得到的机床铣削薄壁曲面零件过程中沿零件法线方向的刚度分布与步骤3中得到的薄壁曲面零件在不同位置处的法向刚度分布相叠加,得到加工过程中机床和薄壁曲面零件的刚度综合对比图,由此得出薄壁曲面零件的加工质量缺陷区域。
本发明的有益效果:本发明通过机床和薄壁曲面零件的法向刚度分布计算,产生加工过程中的机床和零件自身的刚度分布综合对比图,从而得到薄壁曲面零件加工质量缺陷区域,由此可在加工前计算出零件加工质量的缺陷区域,针对加工质量缺陷区域设定不同的加工工艺,从而为零件的工艺优化和精度保障提供依据,使薄壁曲面零件加工具备更好的精度控制效果,降低零件的加工成本,提供加工效率。
附图说明
图1是发明步骤1中S形曲面零件模型及切削路径示意图。
图2是本发明的机床V51030ABJ加工S形曲面刚度分布图。
图3是本发明的步骤3中S曲面有限元模型示意图。
图4是本发明的步骤3中S曲面在加工过程的法向变形示意图。
图5是本发明的步骤3中S形曲面不同各位置处的法向刚度示意图。
图6是本发明的五轴数控机床加工S形曲面刚度综合对比图。
具体实施方式
下面结合附图和具体实施例对本发明做进一步的说明:一种薄壁曲面零件加工表面质量缺陷快速计算方法,包括如下步骤:
步骤1.产生机床铣削薄壁曲面零件的各坐标轴数控指令:在计算机辅助制造CAM软件中创建加工零件模型,设定铣削路径,产生零件加工的前置指令(x、y、z、i、j、k),所述前置指令中x、y和z分别代表刀具加工点的空间位移坐标,i、j、k分别代表刀具加工点的空间向量的方向,再选定机床类型后生成机床各坐标轴的数控指令,所述数控指令可以表示为(X、Y、Z、A、B)或(X、Y、Z、A、C)或(X、Y、Z、B、C),X、Y、Z分别表示机床各坐标轴中平动轴
Figure BDA00002052614500021
Figure BDA00002052614500022
Figure BDA00002052614500023
数控指令,A、B和C分别表示机床各坐标轴中旋转轴
Figure BDA00002052614500024
Figure BDA00002052614500025
Figure BDA00002052614500026
数控指令。
本实施例中,以刀具两摆的五轴联动数控机床(型号为V51030ABJ)加工薄壁曲面零件(本实施例选择S形薄壁曲面零件为例进行说明),如图1所示在三维建模软件(如UG)中建立薄壁零件曲面模型,设定铣削路径,产生前置指令(x、y、z、i、j、k),生成S形薄壁曲面零件加工过程中机床各坐标轴的数控指令(X、Y、Z、A、B)。
对于本领域的普通技术人员来说,由前置指令(x、y、z、i、j、k)计算得到机床各坐标轴的数控指令(X、Y、Z、A、B)是一个公知过程,因此将本步骤视为现有技术而不再详细描述。
步骤2:计算产生机床加工薄壁曲面零件的刚度分布:根据步骤1中得到的机床各坐标轴数控指令,计算机床铣削薄壁曲面零件过程中各坐标轴合成的刚度矩阵,将刚度矩阵沿零件法线方向投影,得到加工过程中机床的刚度分布。
本步骤的具体过程包含如下步骤:
步骤21:计算产生机床各坐标轴刚度组成的刚度矩阵Kjoint,由公式1表示如下:
Figure BDA00002052614500031
公式1
式中k1,k2,k3,k4,k5分别机床各坐标轴的数控指令所对应的刚度,例如当数控指令的组合为(X、Y、Z、A、B),k1,k2,k3,k4,k5分别表示机床坐标轴
Figure BDA00002052614500034
Figure BDA00002052614500035
Figure BDA00002052614500036
所对应的刚度。
本实施例中,仍然以刀具两摆的五轴联动数控机床(型号为V51030ABJ)加工S形薄壁曲面零件为例,经测得机床五个坐标轴
Figure BDA00002052614500037
Figure BDA00002052614500038
Figure BDA00002052614500039
Figure BDA000020526145000310
Figure BDA000020526145000311
的刚度分别为k1=1.1100×109N/m、k2=1.10186×109N/m、k3=1.0990×109N/m、k4=7.85×107N/m、k5=7.85×107N/m,机床刚度矩阵Kjoint可由公式2表示如下:
K joint = 1.1100 0 0 0 0 0 1.10186 0 0 0 0 0 1.0990 0 0 0 0 0 0.0785 0 0 0 0 0 0.0785 × 10 9 公式2
步骤22:计算产生机床三个平动轴的变换矩阵和三个旋转轴中任意两个旋转轴的变换矩阵。
本步骤中,A1、A2、A3分别表示三个平动轴
Figure BDA000020526145000313
Figure BDA000020526145000314
Figure BDA000020526145000315
的变换矩阵,可由公式3表示。A4、A5、A6表示三个旋转轴
Figure BDA000020526145000316
的变化矩阵,可由公式4表示:
A 1 = 1 0 0 X 0 1 0 0 0 0 1 0 0 0 0 1 , A 2 = 1 0 0 0 0 1 0 Y 0 0 1 0 0 0 0 1 , A 3 = 1 0 0 0 0 1 0 0 0 0 1 Z 0 0 0 1 公式3
A 4 = 1 0 0 0 0 cos A - sin A 0 0 sin A cos A 0 0 0 0 1 , A 5 = cos B 0 sin B 0 0 1 0 0 - sin B 0 cos B 0 0 0 0 1 , A 6 = cos C - sin C 0 0 sin C cos C 0 0 0 0 1 0 0 0 0 1 公式4
本实施例中,以刀具两摆的机床(型号为V51030ABJ)加工S形薄壁曲面零件为例,那么,该机床的三个平动坐标轴的变换矩阵同公式3,两个旋转轴的变换矩阵为公式4中的A4和A5
步骤23:计算机床各坐标轴到刀具加工点的变换矩阵T1,T2,T3,T4,T5
本步骤中,T1表示由坐标轴至刀具加工点的变换矩阵,T2表示由坐标轴至刀具加工点的变换矩阵,T3表示由坐标轴
Figure BDA00002052614500046
至刀具加工点的变换矩阵,T4,T5表示由两个旋转坐标轴(
Figure BDA00002052614500047
Figure BDA00002052614500048
Figure BDA00002052614500049
中任两个)至刀具加工点的变换矩阵,可分别由公式5计算如下:
T1=A1A2A3A4A5,T2=A2A3A4A5,T3=A3A4A5,T4=A4A5,T5=A5    公式5
本实施例中,以刀具两摆的机床(型号为V51030ABJ)加工S形薄壁曲面零件为例,将数控指令(X、Y、Z、A、B)代入公式5后,T1,T2,T3,T4,T5分别由如下公式6表示:
T 1 = cos B 0 sin B X sin A sin B cos A - cos B sin A Y - cos A sin B sin A cos A cos B Z 0 0 0 1 T 2 = cos B 0 sin B 0 sin A sin B cos A - cos B sin A Y - cos A sin B sin A cos A cos B Z 0 0 0 1
T 3 = cos B 0 sin B 0 sin A sin B cos A - cos B sin A 0 - cos A sin B sin A cos A cos B Z 0 0 0 1 , T 4 = cos B 0 sin B 0 sin A sin B cos A - cos B sin A 0 - cos A sin B sin A cos A cos B 0 0 0 0 1
T 5 = cos B 0 sin B 0 0 1 0 0 - sin B 0 cos B 0 0 0 0 1 公式6
步骤24:根据步骤23中得到的变换矩阵T1,T2,T3,T4,T5,计算得到机床的Jacobi矩阵(雅可比矩阵)J的各元素值,其具体过程如下:
本步骤中,将步骤23中得到的变换矩阵T1,T2,T3,T4,T5中各元素值采用公式7中的统一的抽象符号表示,如T1,T2,T3,T4,T5统一抽象符号为变换矩阵T,nx表示T1,T2,T3,T4,T5的第一行第一列元素统一抽象符号,ox表示T1,T2,T3,T4,T5的第一行第二列元素统一抽象符号,……,以此类推得到全部的统一抽象符号。
T = n x o x a x p x n y o y a y p y n z o z a z p z 0 0 0 1 公式7
Jacobi矩阵J中每一列的元素由公式8产生。
J1i=(-nxpy+nypx),J2i=(-oxpy+oypx),J3i=(-axpy+aypx),公式8
J4i=nz,J5i=oz,J6i=az
公式8中i=1,2,3,4,5。当i=1时,代入变换矩阵T1中的元素值计算Jacobi矩阵J中第一列的元素J11,J21,J31,J41,J51,J61,当i=2时,代入变换矩阵T2中的元素值计算Jacobi矩阵J中第二列的元素J12,J22,J32,J42,J52,J62,依次类推,直至计算产生Jacobi矩阵J中的全部元素。
本实施例中,以刀具两摆的机床(型号为V51030ABJ)加工S形薄壁曲面零件为例,将数控指令(X、Y、Z、A、B)代入公式8后得到如公式9所示的具体的Jacobi矩阵:
J = X sin A sin B - Y cos B - Y cos B 0 0 0 X cos A 0 0 0 0 - Y sin B - X cos B sin A - Y sin B 0 0 0 - cos A sin B - cos A sin - cos A sin B - cos A sin B - sin B sin A sin A sin A sin A 0 cos A cos B cos A cos B cos A cos B cos A cos B cos B
公式9
步骤25:将步骤21和步骤24中得到的刚度矩阵Kjoint和Jacobi矩阵J代入如下公式10,计算得到机床铣削薄壁曲面零件过程中各坐标轴合成的刚度矩阵Kj
Kj=(J(Kjoint)-1(J)T)-1                                公式10
公式10中,(J)T中的代表Jacobi矩阵J的转置矩阵,符号-1表示对矩阵进行求逆运算。
步骤26:产生机床铣削薄壁曲面零件过程中沿零件法线方向的刚度分布:将步骤1中得到的前置指令(x、y、z、i、j、k)代入公式11中计算零件法线方向的单位矢量nj,将步骤25中得到的机床合成刚度矩阵Kj代入公式12中,得到机床铣削薄壁曲面零件过程中沿零件法线方向的刚度分布Kn(t)。
n j = ( ( x t 1 - x t 0 ) , ( y t 1 - y t 0 ) , ( z t 1 - z t 0 ) ) × ( i , j , k ) | ( ( x t 1 - x t 0 ) , ( y t 1 - y t 0 ) , ( z t 1 - z t 0 ) ) × ( i , j , k ) | 公式11
公式11中,xt0,yt0,zt0表示当前时刻的机床前置指令,xt1,yt1,zt1表示下一时刻的机床前置指令。
Kn(t)=Kjnj                                                    公式12
本实施例中,以刀具两摆的机床(型号为V51030ABJ)加工S形薄壁曲面零件为例,将机床铣削薄壁曲面零件过程中各时刻的合成刚度矩阵Kj与S形曲面法线方向单位矢量nj按公式12计算,即可得到机床加工曲面的刚度分布Kn(t),据此画出机床加工S形薄壁曲面零件过程中的刚度分布图,如图2所示,图中法向直线长短表示机床刚度在不同加工位置处的刚度大小,法线直线越长,刚度越大。
步骤3:计算产生薄壁曲面零件在加工过程中自身的刚度分布,本步骤的具体过程包含如下步骤:
步骤31:产生薄壁曲面零件的有限元模型:在有限元分析软件中建立切削零件的几何建模。该几何建模采用壳单元类型、映射方法划分网格,设定边界条件模拟加工中零件定位与装夹情况。对于本领域的普通技术人员来说,在有限元分析软件中进行切削零件的几何建模是一个公知过程,因此将本步骤视为现有技术而不再详细描述。
本实施例中,以刀具两摆的机床(型号为V51030ABJ)加工薄壁曲面零件为例,在有限元分析软件(如ANSYS)中建立如图3所示的几何模型,采用壳单元映射方法划分网格。
步骤32:计算得到薄壁曲面零件在不同位置处的法向刚度分布:在有限元分析软件中计算产生薄壁曲面零件切削过程中切削力作用下的法向变形,将不同位置处加载下的零件法向变形叠加,再将变形数值求倒数,即为零件在不同位置处的法向刚度。
本实施例中,以刀具两摆的机床(型号为V51030ABJ)加工S型薄壁曲面零件为例,将有限元软件(如ANSYS)计算得到的S形薄壁曲面零件不同位置的法向变形按照零件轮廓位置叠加,得到S型薄壁曲面零件在加工过程的法向变形分布,如图4所示,再将各位置处的法向变形数值求倒数,得到S型薄壁曲面零件不同位置处的法向刚度分布,如图5所示。
步骤4:计算薄壁曲面零件的加工质量缺陷区域:将步骤2中得到的机床铣削薄壁曲面零件过程中沿零件法线方向的刚度分布Kn(t)与步骤3中得到的薄壁曲面零件在不同位置处的法向刚度分布相叠加,得到加工过程中机床和薄壁曲面零件的刚度综合对比图,由此得出薄壁曲面零件的加工质量缺陷区域。
本实施例中,以刀具两摆的机床(型号为V51030ABJ)加工S形薄壁曲面零件为例,图6中S形薄壁曲面零件轮廓的外侧(图中上部)表示机床铣削零件过程中沿零件法线方向的刚度分布,S形薄壁曲面零件轮廓的内侧(图中下部)表示零件在加工过程中不同位置处的法向刚度分布,由图6中可以得出,在(1)、(2)、(3)、(4)位置是机床和S形薄壁曲面零件低刚度区域,也就是薄壁曲面零件的加工质量缺陷区域。由此可在切削前计算出零件加工表面质量缺陷区域,针对表面质量缺陷区域设定不同的加工工艺,从而为零件的工艺优化和精度保障提供依据。
尽管本发明的具体实施例中,仅仅选取型号为V51030ABJ的刀具两摆的机床(五轴联动数控机床)为例进行说明,但是,本领域的普通技术人员应该意识到,机床具体类型的选择不影响本发明的实施,所有其他类型的五轴联动数控机床均可以实施本发明的技术方案。
尽管本发明的具体实施例中,薄壁曲面零件采用S形薄壁曲面零件作为机床加工和零件的具体实施例,但是由于S形薄壁曲面零件具有开角加工区、闭角加工区和开闭角转换加工区,以及缘条型面与基座平面的夹角是变化的等加工特征,相比于只有开角加工区、型面与底平面成定角度的圆锥、圆柱、NAS件以及其他三维曲面零件,其对加工质量反映更为全面,另外在开闭角转换加工区、以及缘条型面与基座平面的变化夹角等加工特征对机床各坐标轴联动的刚度特性要求更高,因此S型薄壁曲面零件比NAS件和其他三维曲面零件加工要求更为严格,详细情况请见美国专利申请US2010004777A1的描述,因此本领域的普通技术人员应该意识到,对于其他更为简化的三维曲面零件也适用于本发明,因此不再对其他更为简化的三维曲面零件的实施例进行举例。
上述具体实施例中,对于属于现有技术的步骤的具体内容未详细介绍,但是这些并不影响本发明的实施,本领域的普通技术人员完全可以根据本发明介绍的背景技术或本领域的其他现有技术实施这些本申请中未详细公开的步骤。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。

Claims (9)

1.薄壁曲面零件加工质量缺陷区域计算方法,其特征在于,包括如下步骤:
步骤1.产生机床铣削薄壁曲面零件的各坐标轴数控指令;
步骤2:计算产生机床加工薄壁曲面零件的刚度分布:根据步骤1中得到的机床各坐标轴数控指令,计算机床铣削薄壁曲面零件过程中各坐标轴合成的刚度矩阵,将刚度矩阵沿零件法线方向投影,得到加工过程中机床的刚度分布;
步骤3:计算产生薄壁曲面零件在加工过程中自身的刚度分布;
步骤4:计算薄壁曲面零件的加工质量缺陷区域:将步骤2中得到的机床铣削薄壁曲面零件过程中沿零件法线方向的刚度分布Kn(t)与步骤3中得到的薄壁曲面零件在不同位置处的法向刚度分布相叠加,得到加工过程中机床和薄壁曲面零件的刚度综合对比图,由此得出薄壁曲面零件的加工质量缺陷区域。
2.根据权利要求1所述的薄壁曲面零件加工质量缺陷区域计算方法,其特征在于,上述步骤2包含如下具体步骤:
步骤21:计算产生机床各坐标轴刚度组成的刚度矩阵Kjoint
步骤22:计算产生机床三个平动轴的变换矩阵和三个旋转轴中任意两个旋转轴的变换矩阵;
步骤23:计算机床各坐标轴到刀具加工点的变换矩阵T1,T2,T3,T4,T5
步骤24:根据步骤23中得到的变换矩阵T1,T2,T3,T4,T5,计算得到机床的Jacobi矩阵J的各元素值;
步骤25:根据步骤21和步骤24中得到的刚度矩阵Kjoint和Jacobi矩阵计算得到机床铣削薄壁曲面零件过程中各坐标轴合成的刚度矩阵Kj
步骤26:产生机床铣削零件过程中沿零件法线方向的刚度分布。
3.根据权利要求2所述的薄壁曲面零件加工质量缺陷区域计算方法,其特征在于,上述步骤21的具体过程为:计算产生机床各坐标轴刚度组成的刚度矩阵Kjoint,由公式1表示如下:
Figure FDA00002052614400011
公式1。
4.根据权利要求2所述的薄壁曲面零件加工质量缺陷区域计算方法,其特征在于,上 述步骤22的具体过程为:计算产生机床三个平动轴的变换矩阵和三个旋转轴中任意两个旋转轴的变换矩阵;本步骤中,A1、A2、A3分别表示三个平动轴 
Figure FDA00002052614400021
Figure FDA00002052614400022
Figure FDA00002052614400023
的变换矩阵,可由公式3表示;A4、A5、A6表示三个旋转轴 
Figure FDA00002052614400024
的变化矩阵,可由公式4表示:
Figure FDA00002052614400025
Figure FDA00002052614400026
Figure FDA00002052614400027
公式3;
Figure FDA00002052614400028
Figure FDA000020526144000210
公式4。
5.根据权利要求2所述的薄壁曲面零件加工质量缺陷区域计算方法,其特征在于,上述步骤23的具体过程为:计算机床各坐标轴到刀具加工点的变换矩阵T1,T2,T3,T4,T5:本步骤中,T1表示由坐标轴 
Figure FDA000020526144000211
至刀具加工点的变换矩阵,T2表示由坐标轴 
Figure FDA000020526144000212
至刀具加工点的变换矩阵,T3表示由坐标轴 
Figure FDA000020526144000213
至刀具加工点的变换矩阵,T4,T5表示由两个旋转坐标轴( 
Figure FDA000020526144000214
Figure FDA000020526144000216
中任两个)至刀具加工点的变换矩阵,可分别由公式5计算如下:
T1=A1A2A3A4A5,T2=A2A3A4A5,T3=A3A4A5,T4=A4A5,T5=A5                公式5。
6.根据权利要求2所述的薄壁曲面零件加工质量缺陷区域计算方法,其特征在于,上述步骤24的具体过程为:本步骤中,将步骤23中得到的变换矩阵T1,T2,T3,T4,T5中各元素值采用公式7中的统一的抽象符号表示,如T1,T2,T3,T4,T5统一抽象符号为变换矩阵T,nx表示T1,T2,T3,T4,T5的第一行第一列元素统一抽象符号,ox表示T1,T2,T3,T4,T5的第一行第二列元素统一抽象符号,……,以此类推得到全部的统一抽象符号;
Figure FDA000020526144000217
公式7;
Jacobi矩阵J中每一列的元素由公式8产生;
J1i=(-nxpy+nypx),J2i=(-oxpy+oypx),J3i=(-axpy+aypx),公式8;
J4i=nz,J5i=oz,J6i=az
公式8中i=1,2,3,4,5;当i=1时,代入变换矩阵T1中的元素值计算Jacobi矩阵J中第一列的元素J11,J21,J31,J41,J51,J61,当i=2时,代入变换矩阵T2中的元素值计算Jacobi矩阵J中第二列的元素J12,J22,J32,J42,J52,J62,依次类推,直至计算产生Jacobi矩阵J中的全部元素。
7.根据权利要求2所述的薄壁曲面零件加工质量缺陷区域计算方法,其特征在于,上述步骤25的具体过程为:将步骤21和步骤24中得到的刚度矩阵Kjoint和Jacobi矩阵J代入如下公式10,计算得到机床铣削薄壁曲面零件过程中各坐标轴合成的刚度矩阵Kj
Kj=(J(Kjoint)-1(J)T)-1                                        公式10;
公式10中,(J)T中的代表Jacobi矩阵J的转置矩阵,符号-1表示对矩阵进行求逆运算。
8.根据权利要求2所述的薄壁曲面零件加工质量缺陷区域计算方法,其特征在于,上述步骤26的具体过程为:将步骤1中得到的前置指令(x、y、z、i、j、k)代入公式11中计算零件法线方向的单位矢量nj,将步骤25中得到的机床合成刚度矩阵Kj代入公式12中,得到机床铣削薄壁曲面零件过程中沿零件法线方向的刚度分布Kn(t);
Figure FDA00002052614400031
公式11;
公式11中,xt0,yt0,zt0表示当前时刻的机床前置指令,xt1,yt1,zt1表示下一时刻的机床前置指令;
Kn(t)=Kjnj                            公式12。
9.根据权利要求1所述的薄壁曲面零件加工质量缺陷区域计算方法,其特征在于,上述步骤3包括如下具体步骤:
步骤31:产生薄壁曲面零件的有限元模型:在有限元分析软件中建立切削零件的几何建模;
步骤32:计算得到薄壁曲面零件在不同位置处的法向刚度分布:在有限元分析软件中计算产生薄壁曲面零件切削过程中切削力作用下的法向变形,将不同位置处加载下的零件法向变形叠加,再将变形数值求倒数,即为零件在不同位置处的法向刚度。 
CN201210303822.8A 2012-08-24 2012-08-24 薄壁曲面零件加工质量缺陷区域计算方法 Active CN102799143B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210303822.8A CN102799143B (zh) 2012-08-24 2012-08-24 薄壁曲面零件加工质量缺陷区域计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210303822.8A CN102799143B (zh) 2012-08-24 2012-08-24 薄壁曲面零件加工质量缺陷区域计算方法

Publications (2)

Publication Number Publication Date
CN102799143A true CN102799143A (zh) 2012-11-28
CN102799143B CN102799143B (zh) 2014-09-03

Family

ID=47198272

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210303822.8A Active CN102799143B (zh) 2012-08-24 2012-08-24 薄壁曲面零件加工质量缺陷区域计算方法

Country Status (1)

Country Link
CN (1) CN102799143B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103116678A (zh) * 2013-02-22 2013-05-22 厦门大学 一种用于气囊抛光的气囊外形的设计方法
CN104597836A (zh) * 2015-01-23 2015-05-06 哈尔滨理工大学 外覆盖件模具四轴加工系统综合刚度性能进行切削的方法
CN105388842A (zh) * 2015-12-10 2016-03-09 华中科技大学 一种零件表面加工缺陷的定位方法
TWI760675B (zh) * 2020-01-06 2022-04-11 財團法人工業技術研究院 加工路徑缺陷檢測方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4990840A (en) * 1988-02-19 1991-02-05 The Cross Company Method and system for controlling a machine tool such as a turning machine
CN101488017A (zh) * 2009-02-24 2009-07-22 上海奈凯电子科技有限公司 基于机械视觉的数控机床加工刀具路径实时控制的方法
US20090307915A1 (en) * 2006-04-26 2009-12-17 Renishaw Plc Differential Calibration
US20100004777A1 (en) * 2007-01-16 2010-01-07 Chengdu Aircraft INdustrial (Group) Co., Ltd S-shape detection test piece and a detection method for detecting the precision of the numerical control milling machine
CN101870073A (zh) * 2010-06-11 2010-10-27 华中科技大学 基于工艺系统刚度特性的多轴数控加工刀具运动规划方法
JP4912824B2 (ja) * 2006-10-12 2012-04-11 株式会社シギヤ精機製作所 研削盤
CN102528554A (zh) * 2010-12-09 2012-07-04 中国科学院沈阳计算技术研究所有限公司 一种五轴加工奇异区域的轨迹优化方法
CN102592035A (zh) * 2012-03-20 2012-07-18 北京航空航天大学 一种车铣复合切削加工表面粗糙度及表面形貌仿真预测方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4990840A (en) * 1988-02-19 1991-02-05 The Cross Company Method and system for controlling a machine tool such as a turning machine
US20090307915A1 (en) * 2006-04-26 2009-12-17 Renishaw Plc Differential Calibration
JP4912824B2 (ja) * 2006-10-12 2012-04-11 株式会社シギヤ精機製作所 研削盤
US20100004777A1 (en) * 2007-01-16 2010-01-07 Chengdu Aircraft INdustrial (Group) Co., Ltd S-shape detection test piece and a detection method for detecting the precision of the numerical control milling machine
CN101488017A (zh) * 2009-02-24 2009-07-22 上海奈凯电子科技有限公司 基于机械视觉的数控机床加工刀具路径实时控制的方法
CN101870073A (zh) * 2010-06-11 2010-10-27 华中科技大学 基于工艺系统刚度特性的多轴数控加工刀具运动规划方法
CN102528554A (zh) * 2010-12-09 2012-07-04 中国科学院沈阳计算技术研究所有限公司 一种五轴加工奇异区域的轨迹优化方法
CN102592035A (zh) * 2012-03-20 2012-07-18 北京航空航天大学 一种车铣复合切削加工表面粗糙度及表面形貌仿真预测方法

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103116678A (zh) * 2013-02-22 2013-05-22 厦门大学 一种用于气囊抛光的气囊外形的设计方法
CN104597836A (zh) * 2015-01-23 2015-05-06 哈尔滨理工大学 外覆盖件模具四轴加工系统综合刚度性能进行切削的方法
CN105388842A (zh) * 2015-12-10 2016-03-09 华中科技大学 一种零件表面加工缺陷的定位方法
TWI760675B (zh) * 2020-01-06 2022-04-11 財團法人工業技術研究院 加工路徑缺陷檢測方法

Also Published As

Publication number Publication date
CN102799143B (zh) 2014-09-03

Similar Documents

Publication Publication Date Title
CN102298358B (zh) 五轴数控加工双nurbs刀具轨迹速度规划方法
CN105425727B (zh) 五轴侧铣加工刀具路径光顺方法
CN106502202B (zh) 一种球头铣刀与导向叶片接触区域的半解析建模方法
CN105242637A (zh) 一种航空薄壁叶片补偿加工方法
CN102566494B (zh) 一种基于刀具矢量平滑压缩的五轴数控插补方法
CN101870073B (zh) 基于工艺系统刚度特性的多轴数控加工刀具运动规划方法
CN105700466A (zh) 高速数控加工轨迹的曲率光顺方法
US20130144425A1 (en) Five-axis flank milling system for machining curved surface and a toolpath planning method thereof
CN102722137A (zh) 直纹面叶轮五坐标插铣加工方法
CN104460516A (zh) 一种基于后置处理五轴刀具半径补偿方法
CN102799143B (zh) 薄壁曲面零件加工质量缺陷区域计算方法
CN105243218A (zh) 机床热误差精度转化和模型建立方法
CN106873522A (zh) 一种非轴对称扫掠曲面的数控车削刀具轨迹规划方法
CN105458372A (zh) 基于非可展直纹面的侧铣误差补偿装置及其刀位规划方法
CN102608952A (zh) 对采用球头刀具的五轴联动机床平滑加工路径的方法
CN102096391B (zh) 汽轮机长叶片型面的双刀头式加工的刀具轨迹生成方法
CN105005265A (zh) 一种基于曲线光滑变形的数控加工进给率规划方法
CN104181865A (zh) 一种整体叶轮粗加工环形走刀路径规划方法
CN110032140A (zh) 一种五轴加工中球形刀刀轴矢量规划方法
CN102707664B (zh) 一种对五轴加工刀具平滑加工路径的方法
CN107239603A (zh) 五轴数控机床加工中基于精细积分的球头铣刀颤振稳定域叶瓣图建模方法
CN106735464B (zh) 一种工件的插铣加工方法
CN105549535A (zh) 五轴摇篮式数控机床未变形切屑的三维几何建模方法
CN105252232B (zh) 一种加工汽轮机动叶片叶身曲面的方法
CN103926873A (zh) 一种在路径规划阶段消除五轴奇异问题的方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant