自由曲线外形零件圆周铣削过程中瞬时铣削力建模方法.pdf

摘要
申请专利号：	CN201010535532.7	申请日：	2010.11.04
公开号：	CN102004818A	公开日：	2011.04.06
当前法律状态：	授权	有效性：	有权
法律详情：	授权\|\|\|实质审查的生效IPC(主分类):G06F 17/50申请日:20101104\|\|\|公开
IPC分类号：	G06F17/50	主分类号：	G06F17/50
申请人：	西北工业大学
发明人：	张卫红; 杨昀; 党建卫; 王怡婷; 万敏
地址：	710072 陕西省西安市友谊西路127号
优先权：
专利代理机构：	西北工业大学专利中心 61204	代理人：	黄毅新
PDF下载：	PDF下载

内容摘要

本发明公开了一种自由曲线外形零件圆周铣削过程中瞬时铣削力建模方法，用于解决现有的自由曲线外形零件圆周铣削过程中瞬时铣削力建模方法计算工作量大的技术问题。技术方案是以理论刀轨上的等效刀具位置点对应的等效进给方向、等效法向量和等效曲率代替实际刀轨上的刀具位置点对应的进给方向、法向量和曲率，避免用非二阶连续的实际刀轨近似理论刀轨时所导致的预测铣削力的突变的现象；在考虑偏心的同时，推导出瞬时未变形切屑厚度的解析模型，无需用数值计算的方法计算瞬时未变形切屑厚度。利用本发明的方法，用MATLAB在个人计算机(Intel Core(TM)2Duo Processor，2.4GHz，2GB)计算10个周期的瞬时铣削力的时间为1.9s，比采用文献2的方法计算时间233.3s，计算效率提高122.79倍。

权利要求书

1.一种自由曲线外形零件圆周铣削过程中瞬时铣削力建模方法，其特征在于包括下述步骤：(1)选定立铣刀和工件几何参数，包括立铣刀的半径RAD、螺旋角β₀、刀齿数N_f、刀具偏心参数ρ和λ，得到工件表面矢量方程p(u)＝[X(u)，Y(u)，Z(u)]^T和工件毛坯边界的矢量方程P_W(v)＝[X_W(v)，Y_W(v)，Z_W(v)]^T；设定切削参数，包括进给率V_f、轴向切削深度R_z、径向切削深度R_r、刀具主轴转速n；输入由CAM系统产生的加工工件的实际刀轨信息和采样时间间隔T_s；(2)根据步骤(1)设定的采样时间间隔T_s和实际刀轨信息，按照下式在实际刀轨上计算出各采样时刻对应的刀具位置点；式中，t表示采样时刻；p_en表示直线刀轨段的起点；p_st表示直线刀轨段的终点；o_CTP表示圆弧刀轨段的圆心；R_CTP表示圆弧刀轨段的半径；(3)根据步骤(1)给定的工件表面矢量方程，通过偏置，计算出理论刀轨的矢量方程p_e(u)＝[X_t(u)，Y_t(u)，Z_t(u)]^T，式中，Xt(u)=X(u)+dY′(u)[X′(u)]2+[Y′(u)]2;]]>Yt(u)=Y(u)+dX′(u)[X′(u)]2+[Y′(u)]2]]>Z_t(u)等于实际刀轨信息中刀心点的Z坐标值；(4)按照下式，将步骤(2)中得到的刀具位置点p_a(t)投影到理论刀轨上，得到等效刀具位置点p_e(u(t))，(p_e(u(t))-p_a(t))×n_a(t)＝0解非线性方程组，得到参数u(t)，将参数u(t)带入理论刀轨矢量方程p_e(u)中，即得到等效刀具位置点p_e(u(t))；(5)将步骤(4)中解得的参数u(t)代入下式，即得到等效进给方向、等效法向量以及等效曲率，以等效进给方向、等效法向量以及等效曲率作为刀具位置点处的实际进给方向、实际法向量以及实际曲率；f_e(ut))＝[X′_t(t) Y′_t(t) 0]^Tn_e(u(t))＝[0 0 1]^T×f_e(u(t))Ke(u(t))=Xt′(u(t))Yt′′(u(t))-Xt′′(u(t))Yt′(u(t))((Xt′(u(t)))2+(Yt′(u(t)))2)32]]>同时计算出等效进给方向在整体坐标系中的角位置；θ(t)=arccos(fe(u(t))·i|fe(u(t))|)]]>式中，i＝[1 0 0]^T；(6)将刀具参与切削的区域沿轴向划分为N个等高梁段，根据步骤(1)给定的工件毛坯边界矢量方程、立铣刀几何参数和步骤(2)中得到的刀具位置点，通过下式计算刀刃片{i，j}的切入切出角和切出角；(a)切入角：式中，p_W(v_i，j，en(t))为满足方程|p_a(t)-p_W(v)|＝RAD_i，j²的点，RADi,j=RAD+ρcos[λ-tan(β0)RADz-2(i-1)πNf],i=1,2,···Nf,]]>N_f是刀刃数，j＝1，2，…，N，z是刀刃片{i，j}中点的Z向高度；(b)切出角：刀具切入工件阶段；刀具在其他切削阶段；式中，p_D(t，m)为满足方程组的点，m表示当前刀齿之前的第m个刀齿，m＝1，...，N_fp_W(v_{i，j，B，ex}(t))为满足方程|p_a(t)-p_W(v)|＝RAD_i，j²的点；(7)通过下式计算作用在刀刃片{i，j}的铣削力：式中，K_i，j，T(t)，K_i，j，R(t)分别是与h_i，j(t)相关的切向合径向瞬时铣削力系数，Δa是等高梁段的高度，h_i，j(t，m)是利用步骤(5)中得到的等效曲率K_e(u(t))，通过下式计算凸型曲面；凹型曲面；式中，RTP(t)=|1Ke(u(t))|,]]>ft=VfnNf;]]>(8)将各个刀刃上的力转化到X_S，Y_S和Z_S方向：式中，是刀具在t时刻处与刀刃片{i，j}对应的切削角度，被定义为从Y_S向顺时针到刀刃片{i，j}的中点所转过的角度；(9)对于每个侧刃，将作用在所有刀刃片上的微元力求和，求得t时刻局部坐标系下作用于各个侧刃的铣削合力：FXS(t)=Σi,jFi,j,XS(t)]]>FYS(t)=Σi,jFi,j,YS(t)]]>(10)将局部坐标系下的铣削合力转化到X，Y和Z方向：FX(t)=FXS(t)cosθ(t)-FYS(t)sinθ(t)]]>FY(t)=FXS(t)sinθ(t)+FYS(t)cosθ(t)]]>

说明书

自由曲线外形零件圆周铣削过程中瞬时铣削力建模方法

技术领域

本发明涉及一种圆周铣削过程中瞬时铣削力建模方法，特别是一种自由曲线外形零件圆周铣削过程中瞬时铣削力建模方法。

背景技术

圆周铣削是通过立铣刀在材料表面切除多余的材料层来获得理想的工件形状、尺寸以及表面光洁度的机械加工方法，它是机械制造中加工各种模具、汽车零件以及航空航天零件的最常用的方式之一。在圆周铣削过程中，经常出现刀具磨损、刀具破损、加工颤振以及加工变形等诸多问题，严重影响零件的加工精度和加工效率。为了克服以上问题，需要建立圆周铣削加工中铣削力预测模型。诸多学者已经在具有代表性的平面铣削方面做了广泛而深入的研究。由于自由曲线外形零件的圆周铣削在实际加工中也经常遇到，一些学者在这方面做了相关的研究。

文献1“Z.C.Wei，M.J.Wang，R.G.Ma，L.Wang，Modeling of process geometry in peripheralmilling of curved surfaces，Journal of Materials Processing Technology 210(2010)799-806.”公开了一种自由曲线外形零件圆周铣削过程中铣削力的模型，并系统地给出了铣削力预测的整个流程。然而该模型忽略了刀具偏心对铣削力的影响，并且无法避免利用非二阶连续的实际刀轨预测铣削力时所导致的预测铣削力的突变。

文献2“K.A.Desai，P.K.Agarwal，P.V.M.Rao，Process geometry modeling with cutter runoutfor milling of curved surfaces，International Journal of Machine Tools and Manufacture 49(2009)1015-1028.”公开了一种自由曲线外形零件圆周铣削过程中铣削力的模型，并且考虑了刀具偏心对铣削力的影响。该模型中计算与铣削力相关的瞬时未变形切屑厚度的基本方法如下：

1.计算当前切削刀齿的实际半径值和前m个刀齿的实际半径值。

2.计算直线与圆弧的交点。其中，该直线由当前的切触点和与当前切触点相对应的刀齿位置角确定，圆弧是以对应于当前切触点的前m个刀齿轨迹上的切触点所对应的刀具位置点为圆心，前m个刀齿的实际半径值为半径。

3.计算以上交点到当前刀具位置点的距离。

4.计算可能的瞬时未变形切屑厚度，即以当前刀齿的实际半径值减去交点到当前刀具位置点的距离。

5.重复步骤1-4，直到m等于刀齿数N_f结束。

6.选出N_f个可能的瞬时未变形切屑厚度的最小值，并将该最小值与0比较，取大者作为瞬时未变形切屑厚度。

从以上步骤可以看出，该方法中计算可能的瞬时未变形切削厚度是通过数值计算的方法完成的。这种方法求一次瞬时未变形切屑厚度需要求解N_f个非线性方程，因此计算工作量大。如，用MATLAB在个人计算机(Intel Core(TM)2Duo Processor，2.4GHz，2GB)计算10个周期的瞬时铣削力的时间为233.3s。

发明内容

为了克服现有的自由曲线外形零件圆周铣削过程中瞬时铣削力建模方法计算工作量大的不足，本发明提供一种自由曲线外形零件圆周铣削过程中瞬时铣削力建模方法，该方法首先通过采样时间间隔计算出各时刻所对应的刀具位置点，同时计算出每个刀具位置点对应的几何参数，然后通过数学推导建立局部坐标系下切入/切出角的数值计算模型和瞬时未变形切屑厚度的解析模型，最后建立局部坐标系下瞬时铣削力与瞬时未变形切屑厚度的关系，并将局部坐标系下的瞬时铣削力投影到整体坐标系下，可以减少自由曲线外形零件圆周铣削过程中瞬时铣削力建模方法的计算工作量。

本发明解决其技术问题所采用的技术方案是：一种自由曲线外形零件圆周铣削过程中瞬时铣削力建模方法，其特点是包括下述步骤：

(1)选定立铣刀和工件几何参数，包括立铣刀的半径RAD、螺旋角β₀、刀齿数N_f、刀具偏心参数ρ和λ，得到工件表面矢量方程p(u)＝[X(u)，Y(u)，Z(u)]^T和工件毛坯边界的矢量方程p_W(v)＝[X_W(v)，Y_W(v)，Z_W(v)]^T；设定切削参数，包括进给率V_f、轴向切削深度R_z、径向切削深度R_r、刀具主轴转速n；输入由CAM系统产生的加工工件的实际刀轨信息和采样时间间隔T_s；

(2)根据步骤(1)设定的采样时间间隔T_s和实际刀轨信息，按照下式

在实际刀轨上计算出各采样时刻对应的刀具位置点；

式中，t表示采样时刻；p_en表示直线刀轨段的起点；p_st表示直线刀轨段的终点；o_CTP表示圆弧刀轨段的圆心；R_CTP表示圆弧刀轨段的半径；

(3)根据步骤(1)给定的工件表面矢量方程，通过偏置，计算出理论刀轨的矢量方程p_e(u)＝[X_t(u)，Y_t(u)，Z_t(u)]^T，

式中，

Xt(u)=X(u)+dY′(u)[X′(u)]2+[Y′(u)]2;]]>

Yt(u)=Y(u)+dX′(u)[X′(u)]2+[Y′(u)]2]]>

Z_t(u)等于实际刀轨信息中刀心点的Z坐标值；

(4)按照下式，将步骤(2)中得到的刀具位置点p_a(t)投影到理论刀轨上，得到等效刀具位置点p_e(u(t))，

(p_e(u(t))-p_a(t))×n_a(t)＝0

解非线性方程组，得到参数u(t)，将参数u(t)带入理论刀轨矢量方程p_e(u)中，即得到等效刀具位置点p_e(u(t))；

(5)将步骤(4)中解得的参数u(t)代入下式，即得到等效进给方向、等效法向量以及等效曲率，以等效进给方向、等效法向量以及等效曲率作为刀具位置点处的实际进给方向、实际法向量以及实际曲率；

f_e(u(t))＝[X′_t(t) Y′_t(t) 0]^T

n_e(u(t))＝[0 0 1]^T×f_e(u(t))

Ke(u(t))=Xt′(u(t))Yt′′(u(t))-Xt′′(u(t))Yt′(u(t))((Xt′(u(t)))2+(Yt′(u(t)))2)32]]>

同时计算出等效进给方向在整体坐标系中的角位置；

θ(t)=arccos(fe(u(t))·i|fe(u(t))|)]]>

式中，i＝[1 0 0]^T；

(6)将刀具参与切削的区域沿轴向划分为N个等高梁段，根据步骤(1)给定的工件毛坯边界矢量方程、立铣刀几何参数和步骤(2)中得到的刀具位置点，通过下式计算刀刃片{i，j}的切入切出角和切出角；

(a)切入角：

式中，p_W(v_i，j，en(t))为满足方程|p_a(t)-p_W(v)|＝RAD_i，j²的点，

RADi,j=RAD+ρcos[λ-tan(β0)RADz-2(i-1)πNf],i=1,2,···Nf,]]>N_f是刀刃数，j＝1，2，…，N，z是刀刃片{i，j}中点的Z向高度；

(b)切出角：

刀具切入工件阶段；

刀具在其他切削阶段；

式中，

p_D(t，m)为满足方程组的点，

m表示当前刀齿之前的第m个刀齿，m＝1，…，N_f

p_W(v_{i，j，B，ex}(t))为满足方程|p_a(t)-p_W(v)|＝RAD_i，j²的点；

(7)通过下式计算作用在刀刃片{i，j}的铣削力：

式中，K_i，j，T(t)，K_i，j，R(t)分别是与h_i，j(t)相关的切向合径向瞬时铣削力系数，

Δa是等高梁段的高度，

h_i，j(t，m)是利用步骤(5)中得到的等效曲率K_e(u(t))，通过下式计算

凸型曲面；

凹型曲面；

式中，RTP(t)=|1Ke(u(t))|,]]>ft=VfnNf;]]>

(8)将各个刀刃上的力转化到X_S，Y_S和Z_S方向：

式中，是刀具在t时刻处与刀刃片{i，j}对应的切削角度，被定义为从Y_S向顺时针到刀刃片{i，j}的中点所转过的角度；

(9)对于每个侧刃，将作用在所有刀刃片上的微元力求和，求得t时刻局部坐标系下作用于各个侧刃的铣削合力：

FXS(t)=Σi,jFi,j,XS(t)]]>

FYS(t)=Σi,jFi,j,YS(t)]]>

(10)将局部坐标系下的铣削合力转化到X，Y和Z方向：

FX(t)=FXS(t)cosθ(t)-FYS(t)sinθ(t)]]>

FY(t)=FXS(t)sinθ(t)+FYS(t)cosθ(t)]]>

本发明的有益效果是：以理论刀轨上的等效刀具位置点对应的等效进给方向、等效法向量和等效曲率代替实际刀轨上的刀具位置点对应的进给方向、法向量和曲率，避免用非二阶连续的实际刀轨近似理论刀轨时所导致的预测铣削力的突变的现象；在考虑偏心的同时，推导出瞬时未变形切屑厚度的解析模型，无需用数值计算的方法计算瞬时未变形切屑厚度。利用本发明的方法，用MATLAB在个人计算机(Intel Core(TM)2Duo Processor，2.4GHz，2GB)计算10个周期的瞬时铣削力的时间为1.9s，比采用文献2的方法计算时间233.3s，计算效率提高122.79倍。

下面结合附图和具体实施方式对本发明进一步说明。

附图说明

图1是实施例1的铣削力随时间变化图。

图2是图1中3.7-3.75s时的放大图。

图3是图1中8.6-8.65s时的放大图。

图4是实施例1的刀轨曲率随时间变化图和进给方向的角位置θ(t)随时间变化图。

图5是实施例2的铣削力随时间变化图。

图6是图5中3.2-3.25s时的放大图。

图7是图5中5.7-5.75s时的放大图。

图中，实线代表测量力，点画线代表预测力，点线代表切入阶段，虚线代表稳定切削阶段。

具体实施方式

实施例1：(1)选定半径RAD＝6mm、螺旋角β₀＝30度的三齿硬质合金立铣刀在三坐标立铣床上对铝合金Al7050进行顺铣切削，刀具偏心参数ρ＝0.0026和λ＝31.8°，主轴转速n＝2000RPM，进给率V_f＝300mm/min，轴向切削深度R_z＝10mm，径向切削深度R_r＝3mm，采样时间间隔T_s＝0.0002s，输入工件欲得到表面的矢量方程p(u)＝[X(u)，Y(u)，Z(u)]^T，式中

X(u)＝20+105u(1-u)²+15u²(1-u)+40u³

Y(u)＝5+90u²(1-u)+30u³ u∈[0，1]

0≤Z(u)≤10

毛坯边界的矢量方程p_W(v)＝[X_W(v)，Y_W(v)，Z_W(v)]^T，式中

XW(v)=20+105v(1-v)2+15v2(1-v)+40v3-36v(1-v)49-364v+1184v2-1640v3+820v4]]>

YW(v)=5+90v2(1-v)+30v3+21(1-v)2-36v(1-v)+21v249-364v+1184v2-1640v3+820v4]]>v∈[0，1]

0≤Z_W(v)≤10

X_W(v)＝20

Y_W(v)＝8(1+v)-5v v∈[-1，0]

0≤Z_W(v)≤10

X_W(v)＝60

Y_W(v)＝38(v-1)+35(2-v)v∈[1，2]

0≤Z_W(v)≤10

实际刀轨信息

…

GOT0/13.9002，16.9496，0.0000

CIRCLE/19.9000，16.9996，0.0000，0.0000000，0.0000000，-1.0000000，6.0000，0.0100，0.5000，12.0000，0.0000

GOTO/19.9500，10.9998，0.0000

CIRCLE/20.2133，52.2325，0.0000，0.0000000，0.0000000，-1.0000000，41.2336，0.0100，0.5000，12.0000，0.0000

GOTO/24.7369，11.2478，0.0000

CIRCLE/22.4035，32.7858，0.0000，0.0000000，0.0000000，-1.0000000，21.6640，0.0100，0.5000，12.0000，0.0000

GOTO/28.0904，11.8815，0.0000

CIRCLE/25.3711，22.0905，0.0000，0.0000000，0.0000000，-1.0000000，10.5649，0.0100，0.5000，12.0000，0.0000

GOTO/30.4537，12.8285，0.0000

CIRCLE/28.0568，17.2549，0.0000，0.0000000，0.0000000，-1.0000000，5.0337，0.0100，0.5000，12.0000，0.0000

GOTO/32.2469，14.4654，0.0000

CIRCLE/25.8212，18.5149，0.0000，0.0000000，0.0000000，-1.0000000，7.5952，0.0100，0.5000，12.0000，0.0000

GOTO/33.0908，16.3149，0.0000

GOTO/33.2840，16.9667，0.0000

GOTO/33.4579，17.6565，0.0000

GOTO/33.6215，18.4079，0.0000

GOTO/33.7774，19.2173，0.0000

GOTO/34.4212，22.9586，0.0000

GOTO/34.6224，24.0036，0.0000

GOTO/34.8583，25.0871，0.0000

CIRCLE/58.6472，19.7687，0.0000，0.0000000，0.0000000，1.0000000，24.3761，0.0100，0.5000，12.0000，0.0000

GOTO/36.3976，29.7262，0.0000

CIRCLE/52.0438，22.6689，0.0000，0.0000000，0.0000000，1.0000000，17.1642，0.0100，0.5000，12.0000，0.0000

GOTO/40.1300，35.0248，0.0000

CIRCLE/52.4252，22.0901，0.0000，0.0000000，0.0000000，1.0000000，17.8460，0.0100，0.5000，12.0000，0.0000

GOTO/44.3945，38.0271，0.0000

CIRCLE/55.0895，16.6900，0.0000，0.0000000，0.0000000，1.0000000，23.8675，0.0100，0.5000，12.0000，0.0000

GOTO/49.4660，39.8855，0.0000

CIRCLE/58.0277，4.7902，0.0000，0.0000000，0.0000000，1.0000000，36.1246，0.0100，0.5000，12.0000，0.0000

GOTO/54.8228，40.7723，0.0000

CIRCLE/59.7855，-13.1822，0.0000，0.0000000，0.0000000，1.0000000，54.1822，0.0100，0.5000，12.0000，0.0000

GOTO/59.9500，40.9998，0.0000

GOTO/60.0000，41.0000，0.0000

(2)根据步骤(1)设定的采样时间间隔T_s＝0.0002s，按照下式，在实际刀轨上计算出刀具位置点

式中

t表示采样时刻

fa(t)=pen-pst|pen-pst|]]>

[R]=cosαs-sinαs0sinαscosαs0001]]>

αs=VfTsRCTP]]>

p_en和p_st分别表示直线刀轨段的起点和终点，o_CTP和R_CTP分别表示圆弧刀轨段的圆心和半径。

(3)根据步骤(1)给定的欲得到表面的矢量方程，通过偏置，计算出理论刀轨的矢量方程p_e(u)＝[X_t(u)，Y_t(u)，Z_t(u)]^T，式中

Xt(u)=20+105u(1-u)2+15u2(1-u)+40u3-72u(1-u)49-364u+1184u2-1640u3+820u4]]>

Yt(u)=5+90u2(1-u)+30u3+42(1-u)2-72u(1-u)+42u249-364u+1184u2-1640u3+820u4]]>u∈[0，1]。

Z_t(u)＝0

(4)按照下式，将步骤(2)中得到的刀具位置点p_a(t)投影到理论刀轨上，得到等效刀具位置点，

(p_e(u(t))-p_a(t))×n_a(t)＝0

解非线性方程组，得到参数u(t)，将参数u(t)带入理论导轨矢量方程p_e(u)中，即得到等效刀具位置点p_e(u(t))。

(5)将步骤(4)中解得的参数u(t)代入下式，即得到等效进给方向、等效法向量以及等效曲率，以等效进给方向、等效法向量以及等效曲率作为刀具位置点处的实际进给方向、实际法向量以及实际曲率。

f_e(u(t))＝[X′_t(t) Y′_t(t) 0]^T

n_e(u(t))＝[0 0 1]^T×f_e(u(t))

Ke(u(t))=Xt′(u(t))Yt′′(u(t))-Xt′′(u(t))Yt′(u(t))((Xt′(u(t)))2+(Yt′(u(t)))2)32]]>

同时计算出等效进给方向在整体坐标系中的角位置

θ(t)=arccos(fe(u(t))·i|fe(u(t))|)]]>

式中，i＝[1 0 0]^T。

(6)将刀具参与切削的区域沿轴向划分为N个等高梁段，根据步骤(1)给定的毛坯的边界的矢量方程、立铣刀几何参数和步骤(2)中得到的刀具位置点，通过下式计算刀刃片{i，j}的切入切出角。

(a)切入角的计算

式中p_W(v_i，j，en(t))为满足方程|p_a(t)-p_W(v)|＝RAD_i，j²的点，

RADi,j=RAD+ρcos[λ-tan(β0)RADz-2(i-1)πNf],i=1,2,···Nf,]]>N_f是刀刃数，j＝1，2，…，N，z是刀刃片{i，j}中点的Z向高度。

(b)切出角的计算

刀具切入工件阶段

刀具在其他切削阶段

式中，

p_D(t，m)为满足方程组的点，

m表示当前刀齿之前的第m个刀齿，m＝1，…，N_f

p_W(v_{i，j，B，ex}(t))为满足方程|p_a(t)-p_W(v)|＝RAD_i，j²的点。

(7)通过下式计算作用在刀刃片{i，j}的铣削力：

式中

K_i，j，T(t)，K_i，j，R(t)分别是与h_i，j(t)相关的切向合径向瞬时铣削力系数，

Δa是等高梁段的高度，

h_i，j(t，m)利用步骤5)中得到的等效曲率K_e(u(t))，通过下式计算

凸型曲面

凹型曲面

式中，RTP(t)=|1Ke(u(t))|,]]>ft=VfnNf.]]>

(8)将各个刀刃上的力转化到X_S，Y_S和Z_S方向：

式中，是刀具在t时刻处与刀刃片{i，j}对应的切削角度，被定义为从Y_S向顺时针到刀刃片{i，j}的中点所转过的角度。

(9)对于每个侧刃，将作用在所有刀刃片上的微元力求和，求得t时刻局部坐标系下作用于各个侧刃的铣削合力：

FXS(t)=Σi,jFi,j,XS(t)]]>

FYS(t)=Σi,jFi,j,YS(t)]]>

(10)将局部坐标系下的铣削合力转化到X，Y和Z方向：

FX(t)=FXS(t)cosθ(t)-FYS(t)sinθ(t)]]>

FY(t)=FXS(t)sinθ(t)+FYS(t)cosθ(t)]]>

通过以上方法，即得到自由曲线外形零件圆周铣削过程中瞬时铣削力的模型。

从图1、2、3可以看出，本发明的方法有效地考虑刀具偏心对铣削力的影响，其结果与实际测量结果吻合，验证了本发明的有效性。

从图1、2、3、4可以看出，在稳定切削阶段，切削力、刀轨曲率和进给方向的角位置θ(t)随时间变化连续，没有文献1中用非二阶连续的实际刀轨近似理论刀轨时所导致的切削过程中的几何量和预测铣削力的突变的现象。

实施例2：(1)选定半径RAD＝6mm、螺旋角β₀＝30度的三齿硬质合金立铣刀在三坐标立铣床上对铝合金Al7050进行顺铣切削，刀具偏心参数ρ＝0.0026和λ＝31.8°，主轴转速n＝2000RPM，进给率V_f＝450mm/min，轴向切削深度R_z＝10mm，径向切削深度R_r＝3mm，采样时间间隔T_s＝0.0002s，输入工件欲得到表面的矢量方程p(u)＝[X(u)，Y(u)，Z(u)]^T，式中

X(u)＝20+105u(1-u)²+15u²(1-u)+40u³

Y(u)＝5+90u²(1-u)+30u³ u∈[0，1]

0≤Z(u)≤10

毛坯边界的矢量方程p_W(v)＝[X_W(v)，Y_W(v)，Z_W(v)]^T，式中

XW(v)=20+105v(1-v)2+15v2(1-v)+40v3-36v(1-v)49-364v+1184v2-1640v3+820v4]]>

YW(v)=5+90v2(1-v)+30v3+21(1-v)2-36v(1-v)+21v249-364v+1184v2-1640v3+820v4]]>v∈[0，1]

0≤Z_W(v)≤10

X_W(v)＝20

Y_W(v)＝8(1+v)-5v v∈[-1，0]

0≤Z_W(v)≤10

X_W(v)＝60

Y_W(v)＝38(v-1)+35(2-v)v∈[1，2]

0≤Z_W(v)≤10

实际刀轨信息

…

GOTO/13.9002，16.9496，0.0000

CIRCLE/19.9000，16.9996，0.0000，0.0000000，0.0000000，-1.0000000，6.0000，0.0100，0.5000，12.0000，0.0000

GOTO/19.9500，10.9998，0.0000

CIRCLE/20.2133，52.2325，0.0000，0.0000000，0.0000000，-1.0000000，41.2336，0.0100，0.5000，12.0000，0.0000

GOTO/24.7369，11.2478，0.0000

CIRCLE/22.4035，32.7858，0.0000，0.0000000，0.0000000，-1.0000000，21.6640，0.0100，0.5000，12.0000，0.0000

GOTO/28.0904，11.8815，0.0000

CIRCLE/25.3711，22.0905，0.0000，0.0000000，0.0000000，-1.0000000，10.5649，0.0100，0.5000，12.0000，0.0000

GOTO/30.4537，12.8285，0.0000

CIRCLE/28.0568，17.2549，0.0000，0.0000000，0.0000000，-1.0000000，5.0337，0.0100，0.5000，12.0000，0.0000

GOT0/32.2469，14.4654，0.0000

CIRCLE/25.8212，18.5149，0.0000，0.0000000，0.0000000，-1.0000000，7.5952，0.0100，0.5000，12.0000，0.0000

GOTO/33.0908，16.3149，0.0000

GOTO/33.2840，16.9667，0.0000

GOTO/33.4579，17.6565，0.0000

GOTO/33.6215，18.4079，0.0000

GOTO/33.7774，19.2173，0.0000

GOTO/34.4212，22.9586，0.0000

GOTO/34.6224，24.0036，0.0000

GOTO/34.8583，25.0871，0.0000

CIRCLE/58.6472，19.7687，0.0000，0.0000000，0.0000000，1.0000000，24.3761，0.0100，0.5000，12.0000，0.0000

GOTO/36.3976，29.7262，0.0000

CIRCLE/52.0438，22.6689，0.0000，0.0000000，0.0000000，1.0000000，17.1642，0.0100，0.5000，12.0000，0.0000

GOTO/40.1300，35.0248，0.0000

CIRCLE/52.4252，22.0901，0.0000，0.0000000，0.0000000，1.0000000，17.8460，0.0100，0.5000，12.0000，0.0000

GOTO/44.3945，38.0271，0.0000

CIRCLE/55.0895，16.6900，0.0000，0.0000000，0.0000000，1.0000000，23.8675，0.0100，0.5000，12.0000，0.0000

GOTO/49.4660，39.8855，0.0000

CIRCLE/58.0277，4.7902，0.0000，0.0000000，0.0000000，1.0000000，36.1246，0.0100，0.5000，12.0000，0.0000

GOTO/54.8228，40.7723，0.0000

CIRCLE/59.7855，-13.1822，0.0000，0.0000000，0.0000000，1.0000000，54.1822，0.0100，0.5000，12.0000，0.0000

GOTO/59.9500，40.9998，0.0000

GOTO/60.0000，41.0000，0.0000

...

(2)根据步骤(1)设定的采样时间间隔T_s＝0.0002s，按照下式，在实际刀轨上计算出刀具位置点

式中

t表示采样时刻

fa(t)=pen-pst|pen-pst|]]>

[R]=cosαs-sinαs0sinαscosαs0001]]>

αs=VfTsRCTP]]>

p_en和p_st分别表示直线刀轨段的起点和终点，o_CTP和R_CTP分别表示圆弧刀轨段的圆心和半径。

(3)根据步骤(1)给定的欲得到表面的矢量方程，通过偏置，计算出理论刀轨的矢量方程p_e(u)＝[X_t(u)，Y_t(u)，Z_t(u)]^T，式中

Xt(u)=20+105u(1-u)2+15u2(1-u)+40u3-72u(1-u)49-364u+1184u2-1640u3+820u4]]>

Yt(u)=5+90u2(1-u)+30u3+42(1-u)2-72u(1-u)+42u249-364u+1184u2-1640u3+820u4]]>u∈[0，1]。

Z_t(u)＝0

(4)按照下式，将步骤(2)中得到的刀具位置点p_a(t)投影到理论刀轨上，得到等效刀具位置点，

(p_e(u(t))-p_a(t))×n_a(t)＝0

解非线性方程组，得到参数u(t)，将参数u(t)带入理论导轨矢量方程p_e(u)中，即得到等效刀具位置点p_e(u(t))。

f_e(u(t))＝[X′_t(t) Y′_t(t) 0]^T

n_e(u(t))＝[0 0 1]^T×f_e(u(t))

Ke(u(t))=Xt′(u(t))Yt′′(u(t))-Xt′′(u(t))Yt′(u(t))((Xt′(u(t)))2+(Yt′(u(t)))2)32]]>

同时计算出等效进给方向在整体坐标系中的角位置

θ(t)=arccos(fe(u(t))·i|fe(u(t))|)]]>

式中，i＝[1 0 0]^T。

(6)将刀具参与切削的区域沿轴向划分为N个等高梁段，根据步骤(1)给定的毛坯的边界的矢量方程、立铣刀几何参数和步骤2)中得到的刀具位置点，通过下式计算刀刃片{i，j}的切入切出角。

(a)切入角的计算：

式中p_W(v_i，j，en(t))为满足方程|p_a(t)-p_W(v)|＝RAD_i，j²的点，

RADi,j=RAD+ρcos[λ-tan(β0)RADz-2(i-1)πNf],i=1,2,···Nf,]]>N_f是刀刃数，j＝1，2，…，N，z是刀刃片{i，j}中点的Z向高度。

(b)切出角的计算：

刀具切入工件阶段

刀具在其他切削阶段

式中，

p_D(t，m)为满足方程组的点，

m表示当前刀齿之前的第m个刀齿，m＝1，…，N_f

p_W(v_{i，j，B，ex}(t))为满足方程|p_a(t)-p_W(v)|＝RAD_i，j²的点。

(7)通过下式计算作用在刀刃片{i，j}的铣削力：

式中，K_i，j，T(t)，K_i，j，R(t)分别是与h_i，j(t)相关的切向合径向瞬时铣削力系数，

Δa是等高梁段的高度，

h_i，j(t，m)利用步骤5)中得到的等效曲率K_e(u(t))，通过下式计算

凸型曲面

凹型曲面

式中，RTP(t)=|1Ke(u(t))|,]]>ft=VfnNf.]]>

(8)将各个刀刃上的力转化到X_S，Y_S和Z_S方向：

式中，是刀具在t时刻处与刀刃片{i，j}对应的切削角度，被定义为从Y_S向顺时针到刀刃片{i，j}的中点所转过的角度。

(9)对于每个侧刃，将作用在所有刀刃片上的微元力求和，求得t时刻局部坐标系下作用于各个侧刃的铣削合力：

FXS(t)=Σi,jFi,j,XS(t)]]>

FYS(t)=Σi,jFi,j,YS(t)]]>

(10)将局部坐标系下的铣削合力转化到X，Y和Z方向：

FX(t)=FXS(t)cosθ(t)-FYS(t)sinθ(t)]]>

FY(t)=FXS(t)sinθ(t)+FYS(t)cosθ(t)]]>

通过以上方法，即得到自由曲线外形零件圆周铣削过程中瞬时铣削力的模型。

从图5、6、7可以看出，本发明的方法有效地考虑刀具偏心对铣削力的影响，其结果与实际测量结果吻合，验证了本发明的有效性。

利用本发明的方法，用MATLAB在个人计算机(Intel Core(TM)2Duo Processor，2.4GHz，2GB)计算10个周期的瞬时铣削力的时间为1.9s，比采用文献2的方法计算时间233.3s，计算效率提高122.79倍。

资源描述

《自由曲线外形零件圆周铣削过程中瞬时铣削力建模方法.pdf》由会员分享，可在线阅读，更多相关《自由曲线外形零件圆周铣削过程中瞬时铣削力建模方法.pdf（28页珍藏版）》请在专利查询网上搜索。

1、10申请公布号CN102004818A43申请公布日20110406CN102004818ACN102004818A21申请号201010535532722申请日20101104G06F17/5020060171申请人西北工业大学地址710072陕西省西安市友谊西路127号72发明人张卫红杨昀党建卫王怡婷万敏74专利代理机构西北工业大学专利中心61204代理人黄毅新54发明名称自由曲线外形零件圆周铣削过程中瞬时铣削力建模方法57摘要本发明公开了一种自由曲线外形零件圆周铣削过程中瞬时铣削力建模方法，用于解决现有的自由曲线外形零件圆周铣削过程中瞬时铣削力建模方法计算工作量大的技术问题。技术方案是以。

2、理论刀轨上的等效刀具位置点对应的等效进给方向、等效法向量和等效曲率代替实际刀轨上的刀具位置点对应的进给方向、法向量和曲率，避免用非二阶连续的实际刀轨近似理论刀轨时所导致的预测铣削力的突变的现象；在考虑偏心的同时，推导出瞬时未变形切屑厚度的解析模型，无需用数值计算的方法计算瞬时未变形切屑厚度。利用本发明的方法，用MATLAB在个人计算机INTELCORETM2DUOPROCESSOR，24GHZ，2GB计算10个周期的瞬时铣削力的时间为19S，比采用文献2的方法计算时间2333S，计算效率提高12279倍。51INTCL19中华人民共和国国家知识产权局12发明专利申请权利要求书4页说明书16页附。

3、图7页CN102004831A1/4页21一种自由曲线外形零件圆周铣削过程中瞬时铣削力建模方法，其特征在于包括下述步骤1选定立铣刀和工件几何参数，包括立铣刀的半径RAD、螺旋角0、刀齿数NF、刀具偏心参数和，得到工件表面矢量方程PUXU，YU，ZUT和工件毛坯边界的矢量方程PWVXWV，YWV，ZWVT；设定切削参数，包括进给率VF、轴向切削深度RZ、径向切削深度RR、刀具主轴转速N；输入由CAM系统产生的加工工件的实际刀轨信息和采样时间间隔TS；2根据步骤1设定的采样时间间隔TS和实际刀轨信息，按照下式在实际刀轨上计算出各采样时刻对应的刀具位置点；式中，T表示采样时刻；PEN表示直线刀轨段的。

4、起点；PST表示直线刀轨段的终点；OCTP表示圆弧刀轨段的圆心；RCTP表示圆弧刀轨段的半径；3根据步骤1给定的工件表面矢量方程，通过偏置，计算出理论刀轨的矢量方程PEUXTU，YTU，ZTUT，式中，ZTU等于实际刀轨信息中刀心点的Z坐标值；4按照下式，将步骤2中得到的刀具位置点PAT投影到理论刀轨上，得到等效刀具位置点PEUT，PEUTPATNAT0解非线性方程组，得到参数UT，将参数UT带入理论刀轨矢量方程PEU中，即得到等效刀具位置点PEUT；5将步骤4中解得的参数UT代入下式，即得到等效进给方向、等效法向量以及等效曲率，以等效进给方向、等效法向量以及等效曲率作为刀具位置点处的实际进给。

5、方向、实际法向量以及实际曲率；FEUTXTTYTT0TNEUT001TFEUT权利要求书CN102004818ACN102004831A2/4页3同时计算出等效进给方向在整体坐标系中的角位置；式中，I100T；6将刀具参与切削的区域沿轴向划分为N个等高梁段，根据步骤1给定的工件毛坯边界矢量方程、立铣刀几何参数和步骤2中得到的刀具位置点，通过下式计算刀刃片I，J的切入切出角和切出角；A切入角式中，PWVI，J，ENT为满足方程|PATPWV|RADI，J2的点，NF是刀刃数，J1，2，N，Z是刀刃片I，J中点的Z向高度；B切出角刀具切入工件阶段；刀具在其他切削阶段；式中，PDT，M为满足方程组的。

6、点，M表示当前刀齿之前的第M个刀齿，M1，NFPWVI，J，B，EXT为满足方程|PATPWV|RADI，J2的点；7通过下式计算作用在刀刃片I，J的铣削力权利要求书CN102004818ACN102004831A3/4页4式中，KI，J，TT，KI，J，RT分别是与HI，JT相关的切向合径向瞬时铣削力系数，A是等高梁段的高度，HI，JT，M是利用步骤5中得到的等效曲率KEUT，通过下式计算凸型曲面；凹型曲面；式中，8将各个刀刃上的力转化到XS，YS和ZS方向式中，是刀具在T时刻处与刀刃片I，J对应的切削角度，被定义为从YS向顺时针到刀刃片I，J的中点所转过的角度；9对于每个侧刃，将作用在所有。

7、刀刃片上的微元力求和，求得T时刻局部坐标系下作用于各个侧刃的铣削合力10将局部坐标系下的铣削合力转化到X，Y和Z方向权利要求书CN102004818ACN102004831A4/4页5权利要求书CN102004818ACN102004831A1/16页6自由曲线外形零件圆周铣削过程中瞬时铣削力建模方法技术领域0001本发明涉及一种圆周铣削过程中瞬时铣削力建模方法，特别是一种自由曲线外形零件圆周铣削过程中瞬时铣削力建模方法。背景技术0002圆周铣削是通过立铣刀在材料表面切除多余的材料层来获得理想的工件形状、尺寸以及表面光洁度的机械加工方法，它是机械制造中加工各种模具、汽车零件以及航空航天零件的最。

8、常用的方式之一。在圆周铣削过程中，经常出现刀具磨损、刀具破损、加工颤振以及加工变形等诸多问题，严重影响零件的加工精度和加工效率。为了克服以上问题，需要建立圆周铣削加工中铣削力预测模型。诸多学者已经在具有代表性的平面铣削方面做了广泛而深入的研究。由于自由曲线外形零件的圆周铣削在实际加工中也经常遇到，一些学者在这方面做了相关的研究。0003文献1“ZCWEI，MJWANG，RGMA，LWANG，MODELINGOFPROCESSGEOMETRYINPERIPHERALMILLINGOFCURVEDSURFACES，JOURNALOFMATERIALSPROCESSINGTECHNOLOGY2102。

9、010799806”公开了一种自由曲线外形零件圆周铣削过程中铣削力的模型，并系统地给出了铣削力预测的整个流程。然而该模型忽略了刀具偏心对铣削力的影响，并且无法避免利用非二阶连续的实际刀轨预测铣削力时所导致的预测铣削力的突变。0004文献2“KADESAI，PKAGARWAL，PVMRAO，PROCESSGEOMETRYMODELINGWITHCUTTERRUNOUTFORMILLINGOFCURVEDSURFACES，INTERNATIONALJOURNALOFMACHINETOOLSANDMANUFACTURE49200910151028”公开了一种自由曲线外形零件圆周铣削过程中铣削力的模型。

10、，并且考虑了刀具偏心对铣削力的影响。该模型中计算与铣削力相关的瞬时未变形切屑厚度的基本方法如下00051计算当前切削刀齿的实际半径值和前M个刀齿的实际半径值。00062计算直线与圆弧的交点。其中，该直线由当前的切触点和与当前切触点相对应的刀齿位置角确定，圆弧是以对应于当前切触点的前M个刀齿轨迹上的切触点所对应的刀具位置点为圆心，前M个刀齿的实际半径值为半径。00073计算以上交点到当前刀具位置点的距离。00084计算可能的瞬时未变形切屑厚度，即以当前刀齿的实际半径值减去交点到当前刀具位置点的距离。00095重复步骤14，直到M等于刀齿数NF结束。00106选出NF个可能的瞬时未变形切屑厚度的最。

11、小值，并将该最小值与0比较，取大者作为瞬时未变形切屑厚度。0011从以上步骤可以看出，该方法中计算可能的瞬时未变形切削厚度是通过数值计算的方法完成的。这种方法求一次瞬时未变形切屑厚度需要求解NF个非线性方程，因此计算工作量大。如，用MATLAB在个人计算机INTELCORETM2DUOPROCESSOR，24GHZ，2GB计算10个周期的瞬时铣削力的时间为2333S。说明书CN102004818ACN102004831A2/16页7发明内容0012为了克服现有的自由曲线外形零件圆周铣削过程中瞬时铣削力建模方法计算工作量大的不足，本发明提供一种自由曲线外形零件圆周铣削过程中瞬时铣削力建模方法，该。

12、方法首先通过采样时间间隔计算出各时刻所对应的刀具位置点，同时计算出每个刀具位置点对应的几何参数，然后通过数学推导建立局部坐标系下切入/切出角的数值计算模型和瞬时未变形切屑厚度的解析模型，最后建立局部坐标系下瞬时铣削力与瞬时未变形切屑厚度的关系，并将局部坐标系下的瞬时铣削力投影到整体坐标系下，可以减少自由曲线外形零件圆周铣削过程中瞬时铣削力建模方法的计算工作量。0013本发明解决其技术问题所采用的技术方案是一种自由曲线外形零件圆周铣削过程中瞬时铣削力建模方法，其特点是包括下述步骤00141选定立铣刀和工件几何参数，包括立铣刀的半径RAD、螺旋角0、刀齿数NF、刀具偏心参数和，得到工件表面矢量方程。

13、PUXU，YU，ZUT和工件毛坯边界的矢量方程PWVXWV，YWV，ZWVT；设定切削参数，包括进给率VF、轴向切削深度RZ、径向切削深度RR、刀具主轴转速N；输入由CAM系统产生的加工工件的实际刀轨信息和采样时间间隔TS；00152根据步骤1设定的采样时间间隔TS和实际刀轨信息，按照下式00160017在实际刀轨上计算出各采样时刻对应的刀具位置点；0018式中，T表示采样时刻；PEN表示直线刀轨段的起点；PST表示直线刀轨段的终点；OCTP表示圆弧刀轨段的圆心；RCTP表示圆弧刀轨段的半径；00193根据步骤1给定的工件表面矢量方程，通过偏置，计算出理论刀轨的矢量方程PEUXTU，YTU，Z。

14、TUT，0020式中，002100220023ZTU等于实际刀轨信息中刀心点的Z坐标值；00244按照下式，将步骤2中得到的刀具位置点PAT投影到理论刀轨上，得到等效刀具位置点PEUT，0025PEUTPATNAT0说明书CN102004818ACN102004831A3/16页80026解非线性方程组，得到参数UT，将参数UT带入理论刀轨矢量方程PEU中，即得到等效刀具位置点PEUT；00275将步骤4中解得的参数UT代入下式，即得到等效进给方向、等效法向量以及等效曲率，以等效进给方向、等效法向量以及等效曲率作为刀具位置点处的实际进给方向、实际法向量以及实际曲率；0028FEUTXTTYTT。

15、0T0029NEUT001TFEUT00300031同时计算出等效进给方向在整体坐标系中的角位置；00320033式中，I100T；00346将刀具参与切削的区域沿轴向划分为N个等高梁段，根据步骤1给定的工件毛坯边界矢量方程、立铣刀几何参数和步骤2中得到的刀具位置点，通过下式计算刀刃片I，J的切入切出角和切出角；0035A切入角00360037式中，PWVI，J，ENT为满足方程|PATPWV|RADI，J2的点，0038NF是刀刃数，J1，2，N，Z是刀刃片I，J中点的Z向高度；0039B切出角0040刀具切入工件阶段；0041刀具在其他切削阶段；0042式中，0043说明书CN102004。

16、818ACN102004831A4/16页90044PDT，M为满足方程组的点，0045M表示当前刀齿之前的第M个刀齿，M1，NF00460047PWVI，J，B，EXT为满足方程|PATPWV|RADI，J2的点；00487通过下式计算作用在刀刃片I，J的铣削力004900500051式中，KI，J，TT，KI，J，RT分别是与HI，JT相关的切向合径向瞬时铣削力系数，0052A是等高梁段的高度，005300540055HI，JT，M是利用步骤5中得到的等效曲率KEUT，通过下式计算0056凸型曲面；00570058凹型曲面；0059说明书CN102004818ACN102004831A5/。

17、16页100060式中，00618将各个刀刃上的力转化到XS，YS和ZS方向006200630064式中，是刀具在T时刻处与刀刃片I，J对应的切削角度，被定义为从YS向顺时针到刀刃片I，J的中点所转过的角度；00659对于每个侧刃，将作用在所有刀刃片上的微元力求和，求得T时刻局部坐标系下作用于各个侧刃的铣削合力00660067006810将局部坐标系下的铣削合力转化到X，Y和Z方向006900700071本发明的有益效果是以理论刀轨上的等效刀具位置点对应的等效进给方向、等效法向量和等效曲率代替实际刀轨上的刀具位置点对应的进给方向、法向量和曲率，避免用非二阶连续的实际刀轨近似理论刀轨时所导致的预。

18、测铣削力的突变的现象；在考虑偏心的同时，推导出瞬时未变形切屑厚度的解析模型，无需用数值计算的方法计算瞬时未变形切屑厚度。利用本发明的方法，用MATLAB在个人计算机INTELCORETM2DUOPROCESSOR，24GHZ，2GB计算10个周期的瞬时铣削力的时间为19S，比采用文献2的方法计算时间2333S，计算效率提高12279倍。0072下面结合附图和具体实施方式对本发明进一步说明。附图说明0073图1是实施例1的铣削力随时间变化图。0074图2是图1中37375S时的放大图。0075图3是图1中86865S时的放大图。0076图4是实施例1的刀轨曲率随时间变化图和进给方向的角位置T随时。

19、间变化图。0077图5是实施例2的铣削力随时间变化图。0078图6是图5中32325S时的放大图。0079图7是图5中57575S时的放大图。说明书CN102004818ACN102004831A6/16页110080图中，实线代表测量力，点画线代表预测力，点线代表切入阶段，虚线代表稳定切削阶段。具体实施方式0081实施例11选定半径RAD6MM、螺旋角030度的三齿硬质合金立铣刀在三坐标立铣床上对铝合金AL7050进行顺铣切削，刀具偏心参数00026和318，主轴转速N2000RPM，进给率VF300MM/MIN，轴向切削深度RZ10MM，径向切削深度RR3MM，采样时间间隔TS00002S。

20、，输入工件欲得到表面的矢量方程PUXU，YU，ZUT，式中0082XU20105U1U215U21U40U30083YU590U21U30U3U0，100840ZU100085毛坯边界的矢量方程PWVXWV，YWV，ZWVT，式中00860087V0，100880ZWV100089XWV200090YWV81V5VV1，000910ZWV100092XWV600093YWV38V1352VV1，200940ZWV100095实际刀轨信息00960097GOT0/139002，169496，000000098CIRCLE/199000，169996，00000，00000000，00000000。

21、，10000000，60000，00100，05000，120000，000000099GOTO/199500，109998，000000100CIRCLE/202133，522325，00000，00000000，00000000，10000000，412336，00100，05000，120000，000000101GOTO/247369，112478，000000102CIRCLE/224035，327858，00000，00000000，00000000，10000000，216640，00100，05000，120000，000000103GOTO/280904，118815，000。

22、000104CIRCLE/253711，220905，00000，00000000，00000000，10000000，105649，00100，05000，120000，00000说明书CN102004818ACN102004831A7/16页120105GOTO/304537，128285，000000106CIRCLE/280568，172549，00000，00000000，00000000，10000000，50337，00100，05000，120000，000000107GOTO/322469，144654，000000108CIRCLE/258212，185149，00000，。

23、00000000，00000000，10000000，75952，00100，05000，120000，000000109GOTO/330908，163149，000000110GOTO/332840，169667，000000111GOTO/334579，176565，000000112GOTO/336215，184079，000000113GOTO/337774，192173，000000114GOTO/344212，229586，000000115GOTO/346224，240036，000000116GOTO/348583，250871，000000117CIRCLE/586472，1。

24、97687，00000，00000000，00000000，10000000，243761，00100，05000，120000，000000118GOTO/363976，297262，000000119CIRCLE/520438，226689，00000，00000000，00000000，10000000，171642，00100，05000，120000，000000120GOTO/401300，350248，000000121CIRCLE/524252，220901，00000，00000000，00000000，10000000，178460，00100，05000，120000，0。

25、00000122GOTO/443945，380271，000000123CIRCLE/550895，166900，00000，00000000，00000000，10000000，238675，00100，05000，120000，000000124GOTO/494660，398855，000000125CIRCLE/580277，47902，00000，00000000，00000000，10000000，361246，00100，05000，120000，000000126GOTO/548228，407723，000000127CIRCLE/597855，131822，00000，0000。

26、0000，00000000，10000000，541822，00100，05000，120000，000000128GOTO/599500，409998，000000129GOTO/600000，410000，0000001302根据步骤1设定的采样时间间隔TS00002S，按照下式，在实际刀轨上计算出刀具位置点01310132式中说明书CN102004818ACN102004831A8/16页130133T表示采样时刻0134013501360137PEN和PST分别表示直线刀轨段的起点和终点，OCTP和RCTP分别表示圆弧刀轨段的圆心和半径。01383根据步骤1给定的欲得到表面的矢量方程，。

27、通过偏置，计算出理论刀轨的矢量方程PEUXTU，YTU，ZTUT，式中01390140U0，1。0141ZTU001424按照下式，将步骤2中得到的刀具位置点PAT投影到理论刀轨上，得到等效刀具位置点，0143PEUTPATNAT00144解非线性方程组，得到参数UT，将参数UT带入理论导轨矢量方程PEU中，即得到等效刀具位置点PEUT。01455将步骤4中解得的参数UT代入下式，即得到等效进给方向、等效法向量以及等效曲率，以等效进给方向、等效法向量以及等效曲率作为刀具位置点处的实际进给方向、实际法向量以及实际曲率。0146FEUTXTTYTT0T0147NEUT001TFEUT0148014。

28、9同时计算出等效进给方向在整体坐标系中的角位置01500151式中，I100T。01526将刀具参与切削的区域沿轴向划分为N个等高梁段，根据步骤1给定的毛坯的边界的矢量方程、立铣刀几何参数和步骤2中得到的刀具位置点，通过下式计算刀刃片I，J的切入切出角。说明书CN102004818ACN102004831A9/16页140153A切入角的计算01540155式中PWVI，J，ENT为满足方程|PATPWV|RADI，J2的点，0156NF是刀刃数，J1，2，N，Z是刀刃片I，J中点的Z向高度。0157B切出角的计算0158刀具切入工件阶段0159刀具在其他切削阶段0160式中，01610162。

29、PDT，M为满足方程组的点，0163M表示当前刀齿之前的第M个刀齿，M1，NF01640165PWVI，J，B，EXT为满足方程|PATPWV|RADI，J2的点。01667通过下式计算作用在刀刃片I，J的铣削力016701680169式中0170KI，J，TT，KI，J，RT分别是与HI，JT相关的切向合径向瞬时铣削力系数，0171A是等高梁段的高度，0172说明书CN102004818ACN102004831A10/16页1501730174HI，JT，M利用步骤5中得到的等效曲率KEUT，通过下式计算0175凸型曲面01760177凹型曲面01780179式中，01808将各个刀刃上的力。

30、转化到XS，YS和ZS方向018101820183式中，是刀具在T时刻处与刀刃片I，J对应的切削角度，被定义为从YS向顺时针到刀刃片I，J的中点所转过的角度。01849对于每个侧刃，将作用在所有刀刃片上的微元力求和，求得T时刻局部坐标系下作用于各个侧刃的铣削合力01850186说明书CN102004818ACN102004831A11/16页16018710将局部坐标系下的铣削合力转化到X，Y和Z方向018801890190通过以上方法，即得到自由曲线外形零件圆周铣削过程中瞬时铣削力的模型。0191从图1、2、3可以看出，本发明的方法有效地考虑刀具偏心对铣削力的影响，其结果与实际测量结果吻合，。

31、验证了本发明的有效性。0192从图1、2、3、4可以看出，在稳定切削阶段，切削力、刀轨曲率和进给方向的角位置T随时间变化连续，没有文献1中用非二阶连续的实际刀轨近似理论刀轨时所导致的切削过程中的几何量和预测铣削力的突变的现象。0193实施例21选定半径RAD6MM、螺旋角030度的三齿硬质合金立铣刀在三坐标立铣床上对铝合金AL7050进行顺铣切削，刀具偏心参数00026和318，主轴转速N2000RPM，进给率VF450MM/MIN，轴向切削深度RZ10MM，径向切削深度RR3MM，采样时间间隔TS00002S，输入工件欲得到表面的矢量方程PUXU，YU，ZUT，式中0194XU20105U1。

32、U215U21U40U30195YU590U21U30U3U0，101960ZU100197毛坯边界的矢量方程PWVXWV，YWV，ZWVT，式中01980199V0，102000ZWV100201XWV200202YWV81V5VV1，002030ZWV100204XWV600205YWV38V1352VV1，202060ZWV100207实际刀轨信息02080209GOTO/139002，169496，000000210CIRCLE/199000，169996，00000，00000000，00000000，10000000，60000，00100，05000，120000，0000002。

33、11GOTO/199500，109998，000000212CIRCLE/202133，522325，00000，00000000，00000000，10000000，412336，00100，05000，120000，000000213GOTO/247369，112478，00000说明书CN102004818ACN102004831A12/16页170214CIRCLE/224035，327858，00000，00000000，00000000，10000000，216640，00100，05000，120000，000000215GOTO/280904，118815，000000216C。

34、IRCLE/253711，220905，00000，00000000，00000000，10000000，105649，00100，05000，120000，000000217GOTO/304537，128285，000000218CIRCLE/280568，172549，00000，00000000，00000000，10000000，50337，00100，05000，120000，000000219GOT0/322469，144654，000000220CIRCLE/258212，185149，00000，00000000，00000000，10000000，75952，00100，05。

35、000，120000，000000221GOTO/330908，163149，000000222GOTO/332840，169667，000000223GOTO/334579，176565，000000224GOTO/336215，184079，000000225GOTO/337774，192173，000000226GOTO/344212，229586，000000227GOTO/346224，240036，000000228GOTO/348583，250871，000000229CIRCLE/586472，197687，00000，00000000，00000000，10000000，24。

36、3761，00100，05000，120000，000000230GOTO/363976，297262，000000231CIRCLE/520438，226689，00000，00000000，00000000，10000000，171642，00100，05000，120000，000000232GOTO/401300，350248，000000233CIRCLE/524252，220901，00000，00000000，00000000，10000000，178460，00100，05000，120000，000000234GOTO/443945，380271，000000235CIRCL。

37、E/550895，166900，00000，00000000，00000000，10000000，238675，00100，05000，120000，000000236GOTO/494660，398855，000000237CIRCLE/580277，47902，00000，00000000，00000000，10000000，361246，00100，05000，120000，000000238GOTO/548228，407723，000000239CIRCLE/597855，131822，00000，00000000，00000000，10000000，541822，00100，05000。

38、，120000，000000240GOTO/599500，409998，000000241GOTO/600000，410000，000000242说明书CN102004818ACN102004831A13/16页1802432根据步骤1设定的采样时间间隔TS00002S，按照下式，在实际刀轨上计算出刀具位置点02440245式中0246T表示采样时刻0247024802490250PEN和PST分别表示直线刀轨段的起点和终点，OCTP和RCTP分别表示圆弧刀轨段的圆心和半径。02513根据步骤1给定的欲得到表面的矢量方程，通过偏置，计算出理论刀轨的矢量方程PEUXTU，YTU，ZTUT，式中0。

39、2520253U0，1。0254ZTU002554按照下式，将步骤2中得到的刀具位置点PAT投影到理论刀轨上，得到等效刀具位置点，0256PEUTPATNAT00257解非线性方程组，得到参数UT，将参数UT带入理论导轨矢量方程PEU中，即得到等效刀具位置点PEUT。02585将步骤4中解得的参数UT代入下式，即得到等效进给方向、等效法向量以及等效曲率，以等效进给方向、等效法向量以及等效曲率作为刀具位置点处的实际进给方向、实际法向量以及实际曲率。0259FEUTXTTYTT0T0260NEUT001TFEUT02610262同时计算出等效进给方向在整体坐标系中的角位置说明书CN10200481。

40、8ACN102004831A14/16页1902630264式中，I100T。02656将刀具参与切削的区域沿轴向划分为N个等高梁段，根据步骤1给定的毛坯的边界的矢量方程、立铣刀几何参数和步骤2中得到的刀具位置点，通过下式计算刀刃片I，J的切入切出角。0266A切入角的计算02670268式中PWVI，J，ENT为满足方程|PATPWV|RADI，J2的点，0269NF是刀刃数，J1，2，N，Z是刀刃片I，J中点的Z向高度。0270B切出角的计算0271刀具切入工件阶段0272刀具在其他切削阶段0273式中，02740275PDT，M为满足方程组的点，0276M表示当前刀齿之前的第M个刀齿，M。

41、1，NF02770278PWVI，J，B，EXT为满足方程|PATPWV|RADI，J2的点。02797通过下式计算作用在刀刃片I，J的铣削力02800281说明书CN102004818ACN102004831A15/16页200282式中，KI，J，TT，KI，J，RT分别是与HI，JT相关的切向合径向瞬时铣削力系数，0283A是等高梁段的高度，028402850286HI，JT，M利用步骤5中得到的等效曲率KEUT，通过下式计算0287凸型曲面02880289凹型曲面02900291式中，02928将各个刀刃上的力转化到XS，YS和ZS方向029302940295式中，是刀具在T时刻处与刀。

42、刃片I，J对应的切削角度，被定义为从YS向顺时针到刀刃片I，J的中点所转过的角度。02969对于每个侧刃，将作用在所有刀刃片上的微元力求和，求得T时刻局部坐标系说明书CN102004818ACN102004831A16/16页21下作用于各个侧刃的铣削合力02970298029910将局部坐标系下的铣削合力转化到X，Y和Z方向030003010302通过以上方法，即得到自由曲线外形零件圆周铣削过程中瞬时铣削力的模型。0303从图5、6、7可以看出，本发明的方法有效地考虑刀具偏心对铣削力的影响，其结果与实际测量结果吻合，验证了本发明的有效性。0304利用本发明的方法，用MATLAB在个人计算机I。

43、NTELCORETM2DUOPROCESSOR，24GHZ，2GB计算10个周期的瞬时铣削力的时间为19S，比采用文献2的方法计算时间2333S，计算效率提高12279倍。说明书CN102004818ACN102004831A1/7页22图1说明书附图CN102004818ACN102004831A2/7页23图2说明书附图CN102004818ACN102004831A3/7页24图3说明书附图CN102004818ACN102004831A4/7页25图4说明书附图CN102004818ACN102004831A5/7页26图5说明书附图CN102004818ACN102004831A6/7页27图6说明书附图CN102004818ACN102004831A7/7页28图7说明书附图CN102004818A。

展开阅读全文