CN105373664B - 一种特种传动曲面五轴铣削力建模方法 - Google Patents
一种特种传动曲面五轴铣削力建模方法 Download PDFInfo
- Publication number
- CN105373664B CN105373664B CN201510800855.7A CN201510800855A CN105373664B CN 105373664 B CN105373664 B CN 105373664B CN 201510800855 A CN201510800855 A CN 201510800855A CN 105373664 B CN105373664 B CN 105373664B
- Authority
- CN
- China
- Prior art keywords
- cutting
- axis
- tool
- milling
- curved surface
- 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
- 238000003801 milling Methods 0.000 title claims abstract description 119
- 238000000034 method Methods 0.000 title claims abstract description 67
- 230000005540 biological transmission Effects 0.000 title claims abstract description 52
- 238000005520 cutting process Methods 0.000 claims abstract description 135
- 238000010408 sweeping Methods 0.000 claims abstract description 57
- 238000003754 machining Methods 0.000 claims abstract description 30
- 238000004364 calculation method Methods 0.000 claims abstract description 9
- 241000220317 Rosa Species 0.000 claims description 8
- 238000005516 engineering process Methods 0.000 claims description 8
- 239000011159 matrix material Substances 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 8
- 230000002459 sustained effect Effects 0.000 claims description 7
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 claims description 2
- 238000013459 approach Methods 0.000 claims description 2
- 238000000205 computational method Methods 0.000 claims description 2
- 238000011438 discrete method Methods 0.000 claims description 2
- 238000013178 mathematical model Methods 0.000 claims description 2
- 238000003672 processing method Methods 0.000 claims description 2
- 230000009466 transformation Effects 0.000 claims description 2
- 241000452638 Parasaissetia nigra Species 0.000 claims 1
- 238000004519 manufacturing process Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 3
- 238000004886 process control Methods 0.000 description 2
- 238000007514 turning Methods 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
Landscapes
- Numerical Control (AREA)
Abstract
本发明一种特种传动曲面五轴铣削力建模方法属于特种传动曲面零件精密高效加工领域,特别涉及一种基于瞬时切削扫略面积及主轴转速的五轴加工铣削力建模方法。该方法首先基于刀具切削刃曲线及加工过程中瞬时未变形切削厚度模型,提出加工过程瞬时切削扫略面积计算方法;以此为基础,建立特种传动曲面五轴铣削加工过程瞬时切削扫略面积计算模型;最后基于瞬时切削扫略面积及主轴转速建立特种传动曲面五轴铣削力模型。本发明基于瞬时切削扫略面积及主轴转速,实现了不同主轴转速下特种传动曲面五轴铣削力预测,对主轴转速优选,应用范围广,提高特种传动曲面加工精度及效率。
Description
技术领域
本发明属于特种传动曲面零件精密高效加工领域,涉及一种基于瞬时切削扫略面积及主轴转速的五轴加工铣削力建模方法。
背景技术
研究特种传动曲面的通用五轴铣削加工技术,提高加工精度,对提高传动精度、进而提升高端装备性能具有重要意义。特种传动曲面五轴铣削过程中,铣削力作为铣削过程中的一个重要物理参数,研究铣削力建模方法对优化特种传动曲面五轴铣削加工工艺具有重要指导作用。
以弧齿锥齿轮球面渐开线齿形为例的特种传动曲面,鉴于曲面各点几何特征,如曲率、斜率、主法线等,往往具有非一致性。数控加工过程中曲面各点切削瞬时扫略面积具有时变性,导致铣削力剧烈波动,从而影响特种传动曲面加工质量及刀具寿命,甚至引起机床振动。另一方面,特种传动曲面五轴铣削过程中,主轴转速变化亦引起铣削力变化,主轴转速优选困难。目前,缺乏特种传动曲面五轴铣削力建模方法研究,难以对加工过程进行有效控制及采用最优主轴转速。因此,针对特种传动曲面五轴铣削,研究基于瞬时切削扫略面积及主轴转速的五轴加工铣削力建模方法,对实现特种传动曲面的精密高效加工具有重要意义。
刘长青等人专利公告号为CN103439917A的“基于特征的切削力预测方法”,基于CAD三维模型进行零件特征识别,根据特征识别结果得到加工区域,提取用来生成刀轨的驱动几何信息,并离散成一些点,通过实验测得剪切力作用系数以及刃口力系数,结合切削参数逐点计算出切削力。该方法为针对飞机结构件局部特征建立的切削力预测方法,难以实现曲面铣削力的预测。文献“Modeling and experimental validation of cuttingforces in five-axis ball-end milling based on true tooth trajectory”,He等,International Journal of Advanced Manufacturing Technology,2015,78(1-4),189-197,提出了一种基于刀齿实际运动轨迹的球头铣刀五轴铣削力预测模型,将预测得到的铣削力与实验所测铣削力对比,结果表明预测结果与实验结果吻合度较高。然而,该方法涉及的铣削力建模过程未涉及瞬时切削扫略面积和主轴转速对铣削力带来的影响,因而对特种传动曲面五轴加工过程有效控制及主轴转速优选仍存在较大限制。
发明内容
本发明旨在克服现有技术的缺陷,考虑瞬时切削扫略面积和主轴转速对铣削力带来的影响,发明一种特种传动曲面五轴铣削力建模方法。该方法是一种基于瞬时切削扫略面积及主轴转速的特种传动曲面五轴铣削力建模方法,实现了特种传动曲面五轴铣削力精确预测,为面向几何特征的特种传动曲面五轴铣削主轴转速优选及加工过程控制提供理论与技术支撑,提高特种传动曲面加工质量及效率。
本发明所采用的技术方案一种特种传动曲面五轴铣削力建模方法,其特征在于,首先基于刀具切削刃曲线及加工过程中瞬时未变形切削厚度模型,建立加工过程瞬时切削扫略面积计算方法;建立特种传动曲面五轴铣削加工过程瞬时切削扫略面积计算模型;最后,基于瞬时切削扫略面积及主轴转速建立特种传动曲面五轴铣削力模型。建模方法的具体步骤如下:
1)加工过程瞬时切削扫略面积计算
曲面铣削加工过程中,不同加工位置刀具切削刃接触工件的微元范围时刻变化,导致刀具切削刃参与切削的未变形切削厚度随之变化。为计算加工过程瞬时切削扫略面积,首先建立刀具切削刃模型。通常加工中,球头铣刀的切削刃曲线为恒导程球面螺旋线,以切削刃微元相对于刀尖位置的螺旋滞后角为参数的刀具切削刃曲线表达式为:
式中,ψj=ψ1-(j-1)2π/m,j=1,2,…m;R为刀具半径;ψ为切削刃的位置角;α为公称螺旋角,即在刀具球头和切削刃之间的螺旋角;k为微元切削刃一点处轴向位置角;ψ1为首个切削刃的位置角;ψj为第j个切削刃的位置角;m为刀具刃数。
过球头铣刀轴线对刀具进行剖面分析,依照离散法求取球头铣刀铣削加工过程中的未变形切削厚度近似数学模型,考虑每齿进给的水平分量,表达式为:
fc(θ)=R2(z)-R1(z)+FHsinθ (2)
式中,FH为每齿进给量的水平分量;z为轴向切削深度;R1(z)、R2(z)为相邻刀具铣削位置同一高度下的轴向微元铣削实际切削半径;θ为切削刃微元点的位置角,即切削刃微元点与刀尖的连线与刀具坐标系Xc轴正方向之间的夹角。
曲面五轴铣削加工与平面加工及斜坡加工的不同在于刀具切触区域与加工进给方向的不断变化。针对曲面五轴铣削加工,基于刀具切削刃模型及球头铣刀铣削加工过程未变形切削厚度模型,建立曲面加工过程瞬时切削扫略面积模型。
首先,建立坐标系。坐标系建立方法为:
坐标系1,原点Oc—刀尖点;Yc轴—铣削过程中进给方向水平投影为正方向;Zc轴—刀具轴线向上为正方向;Xc轴—与Yc轴、Zc轴成右手坐标系。
坐标系2,原点O'c—球头铣刀球心;Y′c轴—平行于斜面且沿斜面进给方向为正方向;Z'c轴—垂直于斜面,沿刀具轴线向上为正方向;X'c轴—与Y′c轴、Z'c轴成右手坐标系,且与Xc轴平行。
然后,建立曲面加工过程瞬时切削扫略面积模型。根据曲面各点几何特征的不同,曲面加工过程瞬时切削扫略面积计算分为上坡和下坡两种。
曲面点斜率为正时,为上坡,计算瞬时切削扫略面积;坐标系OcXcYcZc到坐标系O'cX'cY′cZ'c的齐次变换矩阵为:
式中,η为加工路径曲线一点切线与Yc轴夹角,与OcXcYcZc坐标系和O'cX'cY′cZ'c坐标系的Z轴夹角相等,为正值。
刀具切削刃曲线在OcXcYcZc坐标系中的坐标为:
经过齐次变换后,刀具切削刃曲线在O'cX'cY′cZ'c坐标系中的坐标矩阵为:
参与切削的刀具切削刃部分的最大轴向高度zmax及最小轴向高度zmin的约束条件为:
式中,aP为切深。参见附图2,P为刀具切削刃与曲面的交点,cost=(R-aP)/R,tanη=f′(x),其中,t为OcP与Z′c的夹角,为正值,f(x)为加工路径曲线,在OcO′cP平面内最大轴向高度zmax及最小轴向高度zmin可表示为:
则上坡时瞬时切削扫略面积可表示为:
为提高瞬时切削扫略面积计算精度,对模型进行补偿计算,进给量补偿值ΔfH为:
式中,f′(x)为加工路径曲线的一阶导数即曲线斜率;f″(x)为加工路径曲线的二阶导数即曲线斜率的变化速度;η′为下一刀具位置处加工路径曲线切线与Yc轴夹角,为正值。
补偿后的上坡瞬时切削扫略面积为:
曲面点斜率为负时,为下坡,计算瞬时切削扫略面积,坐标系OcXcYcZc到坐标系O'cX'cY′cZ'c的齐次变换矩阵为:
刀具切削刃曲线在OcXcYcZc坐标系中的坐标表示为式(4),刀具切削刃曲线在O'cX'cY′cZ'c坐标系中的坐标矩阵为:
P为刀具切削刃与曲面的交点,cost=(R-aP)/R,tanη=f′(x)。R1(z)、R2(z)为线OcO′c左侧连续两个刀具铣削位置同一高度下的轴向微元铣削实际切削半径,表示为:
R′1(z)、R′2(z)为线OcO′c右侧连续两个刀具铣削位置同一高度下的轴向微元铣削实际切削半径,表示为:
当t>η时,线OcO′c右侧参与切削的刀具切削刃部分最大轴向高度zmax的约束条件为:
式中,θ1为线OcO′c右端刀具切削刃微元点的位置角。OcPO′cQ平面内最大轴向高度z′max为:
z′max=R-Rcos(t-η) (16)
线OcO′c左侧最大轴向高度zmax约束条件为:
式中,θ2为线OcO′c右端刀具切削刃微元点的位置角。OcPO′cQ平面内最大轴向高度zmax为:
zmax=R-Rcosη (18)
则下坡时瞬时切削扫略面积可表示为:
补偿后的下坡瞬时切削扫略面积为:
另一种情况当t<η时,参与切削的刀具切削刃部分的最大轴向高度zmax及最小轴向高度zmin的约束条件同式(6)。
在OcPO′cQ平面内最大轴向高度zmax以及最小轴向高度zmin可表示为:
则下坡时瞬时切削扫略面积可表示为:
补偿后的下坡瞬时切削扫略面积为:
2)特种传动曲面五轴铣削瞬时切削扫掠面积计算
针对特种传动曲面五轴铣削瞬时切削扫掠面积计算,利用NUBRS曲面建模技术建立特种传动曲面模型。基于建立的曲面模型,计算切削点Cp处u方向的切矢和主法矢采用加工方法为侧铣,则切削点Cp处u方向的副法矢为:
切削点Cp处v方向的切矢、法矢求法与u向相同。
将切削点Cp曲面u向和v向曲线沿副法矢方向偏置切深ap距离,求解偏置曲线与刀具球端相交的点围成球面形状,将刀刃曲线代入与该区域求解交点,得出参与切削的刀具切削刃部分的最大轴向高度zmax、最小轴向高度zmin,则特种传动曲面五轴铣削瞬时切削扫掠面积如式(10)、式(20)、式(23)。
3)基于瞬时切削扫略面积及主轴转速的五轴铣削力建模
针对特种传动曲面,给定切深、每齿进给量、刀具螺旋角、刀具半径等工艺参数,利用特种传动曲面五轴铣削瞬时切削扫掠面积计算方法,计算出特种传动曲面五轴铣削瞬时切削扫掠面积,采用最小二乘拟合方法,建立具有5阶精度的铣削力模型,表示为:
F=K4S4+K3S3+K2S2+K1S+K0 (25)
式中,Ki(i=0,1,2,3,4)为切削系数。
根据实测切削力数据及瞬时切削扫略面积,得到不同主轴转速下的切削系数,表示为:
Fn=Kn4S4+Kn3S3+Kn2S2+Kn1S+Kn0 (26)
式中,Fn为不同主轴转速下的切削力,Kni(i=0,1,2,3,4)为不同主轴转速下的切削系数。基于Kni及相应的主轴转速n,利用最小二乘法建立切削系数与主轴转速间的函数模型,表示为:
将式(27)带入(25),最终获得基于瞬时切削扫略面积及主轴转速的特种传动曲面五轴铣削力模型,表示为:
本发明的显著效果和益处是基于瞬时切削扫略面积及主轴转速,发明了一种特种传动曲面五轴铣削力建模方法,可实现不同主轴转速下特种传动曲面五轴铣削力预测,对主轴转速优选,提高特种传动曲面加工精度及效率。
附图说明
图1——特种传动曲面五轴铣削力建模方法整体流程图
图2——上坡段微元铣削示意图;Oc、Xc、Yc、Zc-坐标系1,O'c、X'c、Y′c、Z'c-坐标系2,P-刀具切削刃与曲面的交点,P'-下一刀具铣削位置处切削刃与曲面的交点,A、B-刀具与曲面接触的边界极限点,Q-刀具切触切点,Q'-下一刀具铣削位置切触切点,η-加工路径曲线一点切线与Yc轴夹角,t-OcP与Z′c的夹角,aP-切深,fH、fV-每齿进给量的水平分量和竖直分量。
图3——下坡段微元铣削示意图;Oc、Xc、Yc、Zc-坐标系1,O'c、X'c、Y′c、Z'c-坐标系2,P-刀具切削刃与曲面的交点,P'-下一刀具铣削位置处切削刃与曲面的交点,A、B-刀具与曲面接触的边界极限点,Q-刀具切触切点,Q'-下一刀具铣削位置切触切点,η-加工路径曲线一点切线与Yc轴夹角,t-OcP与Z′c的夹角,aP-切深,fH、fV-每齿进给量的水平分量和竖直分量。
图4——主轴转速为4500rpm铣削力模型验证对比图;A-铣削力测量值,B-铣削力计算值
图5——主轴转速为8400rpm铣削力模型验证对比图;A-铣削力测量值,B-铣削力计算值
具体实施方案
结合附图和技术方案详细说明本发明的具体实施方式。
曲面五轴铣削加工与平面加工及斜坡加工的不同在于刀具切触区域与加工进给方向的不断变化,鉴于特种传动曲面各点几何特征,如曲率、斜率、主法线等,往往具有非一致性,数控加工过程中曲面各点切削瞬时扫略面积具有时变性,导致铣削力剧烈波动。另外,特种传动曲面五轴铣削过程中,主轴转速变化亦引起铣削力变化。据此,针对特种传动曲面五轴铣削力建模难题,发明了一种特种传动曲面五轴铣削力建模方法。图1是一种特种传动曲面五轴铣削力建模方法整体流程图。
以球头铣刀侧铣球面渐开线齿面为例,借助MATLAB软件计算并仿真,详细说明本发明实施过程。
首先,对球面渐开线齿面进行建模,给定加工工艺参数为切深aP=0.1mm,每齿进给量fH=0.02mm/z,刀具半径R=3mm,刀具公称螺旋角α=30°。
图2是上坡段微元铣削示意图,图3是下坡段微元铣削示意图。其中,OcXcYcZc为坐标系1,O'cX'cY′cZ'c为坐标系2,P为刀具切削刃与曲面的交点,P'为下一刀具铣削位置处切削刃与曲面的交点,A、B为刀具与曲面接触的边界极限点,Q为刀具切触切点,Q'为下一刀具铣削位置切触切点,η为加工路径曲线一点切线与Yc轴夹角,t为OcP与Z′c的夹角,aP为切深,fH、fV为每齿进给量的水平分量和竖直分量,
其次,利用式(10)、式(20)、式(23)计算得到切削区域瞬时切削扫略面积。
然后,利用式(25)建立具有5阶精度的铣削力预测模型,根据实测切削力数据及计算出的瞬时切削扫略面积,利用最小二乘法法得到主轴转速为4000rpm、6000rpm、7200rpm、8000rpm、8400rpm及9200rpm时铣削力预测模型为:
最后,利用最小二乘法建立切削系数与主轴转速之间的函数模型,如下:
将式(30)带入式(25),最终获得基于瞬时切削扫略面积及主轴转速的特种传动曲面五轴铣削力模型,如下:
F=(-775.23n4+2.16×107n3-2.20×1011n2+9.62×1014n-1.52×1018S4+
(24.72n4-6.90×105n3+7.02×109n2-3.07×1013n+4.83×1016)S3+
(-0.30n4+8244.57n3-8.38×107n2+3.66×1011n-5.77×1014)S2+ (31)
(1.57×10-3n4-43.75n3+4.45×105n2-1.94×109n+3.06×1012)S+
(-3.00×10-6n4+0.09n3-884.66n2+3.87×106n-6.09×109)
为验证模型的有效性,将主轴转速为4500rpm和8400rpm时测得的铣削力值与模型计算值对比,参见附图4和附图5。结果表明,铣削力预测值与实测值吻合,实现特种传动曲面五轴铣削力预测。
本发明针对特种传动曲面五轴铣削过程中缺乏对加工过程铣削力建模方法,建立了一种新型的基于瞬时切削扫略面积及主轴转速的特种传动曲面五轴铣削力建模方法,为面向几何特征的特种传动曲面五轴铣削主轴转速优选及加工过程控制提供理论与技术支撑,提高特种传动曲面加工质量及效率。
Claims (1)
1.一种特种传动曲面五轴铣削力建模方法,其特征在于,该方法首先基于刀具切削刃曲线及加工过程中瞬时未变形切削厚度模型,建立加工过程瞬时切削扫略面积计算方法;再建立特种传动曲面五轴铣削加工过程瞬时切削扫略面积计算模型;最后,依据基于瞬时切削扫略面积及主轴转速建立特种传动曲面五轴铣削力模型;建模方法的具体步骤如下:
1)加工过程瞬时切削扫略面积计算
曲面铣削加工过程中,不同加工位置刀具切削刃接触工件的微元范围时刻变化,导致刀具切削刃参与切削的未变形切削厚度随之变化;为计算加工过程瞬时切削扫略面积,首先建立刀具切削刃模型;通常加工中,球头铣刀的切削刃曲线为恒导程球面螺旋线,以切削刃微元相对于刀尖位置的螺旋滞后角为参数的刀具切削刃曲线表达式为:
式中,ψj=ψ1-(j-1)2π/m,j=1,2,…m;R为刀具半径;ψ为切削刃的位置角;α为公称螺旋角,即在刀具球头和切削刃之间的螺旋角;k为微元切削刃一点处轴向位置角;ψ1为首个切削刃的位置角;ψj为第j个切削刃的位置角;m为刀具刃数;
过球头铣刀轴线对刀具进行剖面分析,依照离散法求取球头铣刀铣削加工过程中的未变形切削厚度近似数学模型,考虑每齿进给的水平分量,表达式为:
fc(θ)=R2(z)-R1(z)+FHsinθ (2)
式中,FH为每齿进给量的水平分量;z为轴向切削深度;R1(z)、R2(z)为相邻刀具铣削位置同一高度下的轴向微元铣削实际切削半径;θ为切削刃微元点的位置角,即切削刃微元点与刀尖的连线与刀具坐标系Xc轴正方向之间的夹角;
曲面五轴铣削加工与平面加工及斜坡加工的不同在于刀具切触区域与加工进给方向的不断变化;针对曲面五轴铣削加工,基于刀具切削刃模型及球头铣刀铣削加工过程未变形切削厚度模型,建立曲面加工过程瞬时切削扫略面积模型;
首先,建立坐标系;坐标系建立方法为:
坐标系1,原点Oc——刀尖点;Yc轴——铣削过程中进给方向水平投影为正方向;Zc轴——刀具轴线向上为正方向;Xc轴——与Yc轴、Zc轴成右手坐标系;
坐标系2,原点Oc'——球头铣刀球心;Yc'轴——平行于斜面且沿斜面进给方向为正方向;Zc'轴——垂直于斜面,沿刀具轴线向上为正方向;Xc'轴——与Yc'轴、Zc'轴成右手坐标系,且与Xc轴平行;
然后,建立曲面加工过程瞬时切削扫略面积模型;根据曲面各点几何特征的不同,曲面加工过程瞬时切削扫略面积计算分为上坡和下坡两种;
曲面点斜率为正时,为上坡,计算瞬时切削扫略面积,坐标系OcXcYcZc到坐标系Oc'Xc'Yc'Zc'的齐次变换矩阵为:
式中,η为加工路径曲线一点切线与Yc轴夹角,与OcXcYcZc坐标系和Oc'Xc'Yc'Zc'坐标系的Z轴夹角相等,为正值;
刀具切削刃曲线在OcXcYcZc坐标系中的坐标为:
经过齐次变换后,刀具切削刃曲线在Oc'Xc'Yc'Zc'坐标系中的坐标矩阵为:
参与切削的刀具切削刃部分的最大轴向高度zmax及最小轴向高度zmin的约束条件为:
式中,aP为切深;P为刀具切削刃与曲面的交点,cost=(R-aP)/R,tanη=f′(x),其中t为OcP与Zc′的夹角,为正值,f(x)为加工路径曲线,在OcOc′P平面内最大轴向高度zmax及最小轴向高度zmin可表示为:
则上坡时瞬时切削扫略面积可表示为:
式中,θ为刀具切削刃微元点的位置角,且fH、fV为每齿进给量的水平分量和竖直分量,且fV=fHtanη;R1(z),R2(z)为相邻刀具铣削位置同一高度下的轴向微元铣削实际切削半径,且
为提高瞬时切削扫略面积计算精度,对模型进行补偿计算,进给量补偿值ΔfH为:
式中,f′(x)为加工路径曲线的一阶导数即曲线斜率;f″(x)为加工路径曲线的二阶导数即曲线斜率的变化速度;η′为下一刀具位置处加工路径曲线切线与Yc轴夹角,为正值;
补偿后的上坡瞬时切削扫略面积为:
曲面点斜率为负时,为下坡,计算瞬时切削扫略面积,坐标系OcXcYcZc到坐标系Oc'Xc'Yc'Zc'的齐次变换矩阵为:
刀具切削刃曲线在OcXcYcZc坐标系中的坐标表示为式(4),刀具切削刃曲线在Oc'Xc'Yc'Zc'坐标系中的坐标矩阵为:
P为刀具切削刃与曲面的交点,cost=(R-aP)/R,tanη=f′(x);R1(z)、R2(z)为线OcOc′左侧连续两个刀具铣削位置同一高度下的轴向微元铣削实际切削半径,表示为:
R1′(z)、R2′(z)为线OcOc′右侧连续两个刀具铣削位置同一高度下的轴向微元铣削实际切削半径,表示为:
当t>η时,线OcOc′右侧参与切削的刀具切削刃部分最大轴向高度zmax的约束条件为:
式中,θ1为线OcOc′右端刀具切削刃微元点的位置角;OcPOc′Q平面内最大轴向高度z′max为:
z′max=R-Rcos(t-η)(16)
线OcOc′左侧最大轴向高度zmax约束条件为:
式中,θ2为线OcOc′右端刀具切削刃微元点的位置角;OcPOc′Q平面内最大轴向高度zmax为:
zmax=R-Rcosη(18)
则下坡时瞬时切削扫略面积可表示为:
补偿后的下坡瞬时切削扫略面积为:
另一种情况当t<η时,参与切削的刀具切削刃部分的最大轴向高度zmax及最小轴向高度zmin的约束条件同式(6);
在OcPOc′Q平面内最大轴向高度zmax以及最小轴向高度zmin可表示为:
则下坡时瞬时切削扫略面积可表示为:
补偿后的下坡瞬时切削扫略面积为:
2)特种传动曲面五轴铣削瞬时切削扫掠面积计算
针对特种传动曲面五轴铣削瞬时切削扫掠面积计算,利用NUBRS曲面建模技术建立特种传动曲面模型;基于建立的曲面模型,计算切削点Cp处u方向的切矢和主法矢采用加工方法为侧铣,则切削点Cp处u方向的副法矢为:
切削点Cp处v方向的切矢、法矢求法与u向相同;
将切削点Cp曲面u向和v向曲线沿副法矢方向偏置切深ap距离,求解偏置曲线与刀具球端相交的点围成球面形状,将刀刃曲线代入与该区域求解交点,得出参与切削的刀具切削刃部分的最大轴向高度zmax、最小轴向高度zmin,则特种传动曲面五轴铣削瞬时切削扫掠面积如式(10)、式(20)、式(23);
3)基于瞬时切削扫略面积及主轴转速的五轴铣削力建模
针对特种传动曲面,给定切深、每齿进给量、刀具螺旋角、刀具半径工艺参数,利用特种传动曲面五轴铣削瞬时切削扫掠面积计算方法,计算出特种传动曲面五轴铣削瞬时切削扫掠面积,采用最小二乘拟合方法,建立具有5阶精度的铣削力模型,表示为:
F=K4S4+K3S3+K2S2+K1S+K0(25)
式中,Ki(i=0,1,2,3,4)为切削系数;
根据实测切削力数据及瞬时切削扫略面积,得到不同主轴转速下的切削系数,表示为:
Fn=Kn4S4+Kn3S3+Kn2S2+Kn1S+Kn0 (26)
式中,Fn为不同主轴转速下的切削力,Kni(i=0,1,2,3,4)为不同主轴转速下的切削系数;基于Kni及相应的主轴转速n,利用最小二乘法建立切削系数与主轴转速间的函数模型,表示为:
将式(27)带入(25),最终获得基于瞬时切削扫略面积及主轴转速的特种传动曲面五轴铣削力模型,表示为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510800855.7A CN105373664B (zh) | 2015-09-16 | 2015-11-18 | 一种特种传动曲面五轴铣削力建模方法 |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2015105861459 | 2015-09-16 | ||
CN201510586145 | 2015-09-16 | ||
CN201510800855.7A CN105373664B (zh) | 2015-09-16 | 2015-11-18 | 一种特种传动曲面五轴铣削力建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105373664A CN105373664A (zh) | 2016-03-02 |
CN105373664B true CN105373664B (zh) | 2018-07-13 |
Family
ID=55375859
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510800855.7A Active CN105373664B (zh) | 2015-09-16 | 2015-11-18 | 一种特种传动曲面五轴铣削力建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105373664B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106041183B (zh) * | 2016-07-05 | 2017-12-26 | 大连理工大学 | 薄壁曲面稳定铣削的极限轴向切深判定方法 |
CN107045556B (zh) * | 2016-10-14 | 2020-08-11 | 台州市德嘉机电科技有限公司 | 螺旋曲面磨削过程中的磨削力预测方法 |
CN106502202B (zh) * | 2017-01-06 | 2018-10-16 | 大连理工大学 | 一种球头铣刀与导向叶片接触区域的半解析建模方法 |
CN107168245B (zh) * | 2017-05-04 | 2019-08-23 | 武汉理工大学 | 一种考虑刃口效应的倒棱圆形车刀切削力精确预测方法 |
CN109719314B (zh) * | 2019-01-11 | 2020-06-05 | 深圳大学 | 一种车削加工中刀具的对准方法、车削装置 |
CN111759488B (zh) * | 2020-07-09 | 2021-08-24 | 山东大学 | 一种应用于根管预备的变截面镍钛根管锉的设计方法及***与制备 |
CN113255076B (zh) * | 2021-05-26 | 2023-06-02 | 西安理工大学 | 一种球头铣刀垂直加工时刀-工切触区域的识别方法 |
CN116309124B (zh) * | 2023-02-15 | 2023-10-20 | 霖鼎光学(江苏)有限公司 | 一种光学曲面模具的修正方法、电子设备及存储介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006053789A (ja) * | 2004-08-12 | 2006-02-23 | Hitachi Ltd | 多軸数値制御装置用のncポストプロセッサ装置 |
CN102622489A (zh) * | 2012-03-26 | 2012-08-01 | 上海交通大学 | 基于acis平台的五轴侧铣加工切削力预测方法 |
CN103390076A (zh) * | 2013-06-24 | 2013-11-13 | 河南理工大学 | 含跳动复杂路径精确未变形切削厚度的确定方法 |
CN103558808A (zh) * | 2013-09-28 | 2014-02-05 | 大连理工大学 | 复杂曲面五轴数控加工刀矢的运动学控制方法 |
CN103592891A (zh) * | 2013-09-28 | 2014-02-19 | 大连理工大学 | 运动学约束的复杂曲面五轴数控加工刀矢光顺方法 |
-
2015
- 2015-11-18 CN CN201510800855.7A patent/CN105373664B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006053789A (ja) * | 2004-08-12 | 2006-02-23 | Hitachi Ltd | 多軸数値制御装置用のncポストプロセッサ装置 |
CN102622489A (zh) * | 2012-03-26 | 2012-08-01 | 上海交通大学 | 基于acis平台的五轴侧铣加工切削力预测方法 |
CN103390076A (zh) * | 2013-06-24 | 2013-11-13 | 河南理工大学 | 含跳动复杂路径精确未变形切削厚度的确定方法 |
CN103558808A (zh) * | 2013-09-28 | 2014-02-05 | 大连理工大学 | 复杂曲面五轴数控加工刀矢的运动学控制方法 |
CN103592891A (zh) * | 2013-09-28 | 2014-02-19 | 大连理工大学 | 运动学约束的复杂曲面五轴数控加工刀矢光顺方法 |
Non-Patent Citations (3)
Title |
---|
多零件几何要素影响下的装配产品特性预测方法;贾振元 等;《机械工程学报》;20090715;第45卷(第7期);167-173 * |
特种传动曲面五轴铣削工艺方法及工艺参数优化研究;杨彦宇;《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》;20160315(第3期);B022-671 * |
高性能复杂曲面零件测量-再设计-数字加工一体化加工方法;贾振元 等;《机械工程学报》;20131031;第49卷(第19期);126-131 * |
Also Published As
Publication number | Publication date |
---|---|
CN105373664A (zh) | 2016-03-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105373664B (zh) | 一种特种传动曲面五轴铣削力建模方法 | |
CN102637216B (zh) | 一种复杂组合曲面的数控侧铣加工刀轨生成方法 | |
CN102681488B (zh) | 一种铣削加工工件表面形貌的建模方法 | |
CN100418027C (zh) | 一种螺旋线制导的曲面数控加工方法 | |
CN103955169B (zh) | 五轴数控侧铣加工铣削力预报方法 | |
Hao et al. | Analysis of milling surface roughness prediction for thin-walled parts with curved surface | |
CN110262397A (zh) | 车铣加工空间螺旋次摆线运动轨迹及瞬时切削力预测模型 | |
CN107369127A (zh) | 一种沿行切路径进给的三维类摆线抛光轨迹生成方法 | |
CN108109199A (zh) | 一种端铣加工平面成形表面三维建模方法 | |
CN103934569B (zh) | 一种基于选择性激光烧结的分层切片方法 | |
CN104318022A (zh) | 一种预测工件表面粗糙度并提高切削效率的方法 | |
CN103433487A (zh) | 一种提高激光快速成形金属零件表面平整度的方法 | |
CN108445834A (zh) | 一种大型复杂构件机器人加工离线轨迹规划方法 | |
CN112705794A (zh) | 一种用于摆线轮加工的剐齿刀具及其设计方法 | |
CN113704928B (zh) | 基于加工物理过程的倾斜铣削表面形貌创成方法 | |
CN109740269B (zh) | 一种渐开线蜗杆车削加工齿面三维建模方法 | |
CN109530768B (zh) | 一种叶片叶尖削边减薄区的加工方法 | |
CN108345733A (zh) | 螺旋铣刀平面加工成形表面三维建模方法 | |
Jiang et al. | Prediction and simulation of cutting force in hypoid gear machining using forming method | |
CN108229046A (zh) | 一种机械加工车端面工艺中已加工表面的三维建模方法 | |
Vyboishchik | Modelling topology of freeform surfaces with ball-end milling | |
CN105631139A (zh) | 一种劈型光栅刻划刀具参数化设计方法 | |
Fan | Cutting speed modelling in ball nose milling applications | |
Ma et al. | A new cutting force modeling method in high-speed milling of curved surface with difficult-to-machine material | |
CN114065427A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |