一种基于铣削力预测的曲面加工轨迹生成方法
【专利摘要】一种基于铣削力预测的曲面加工轨迹生成方法，涉及金属切削加工领域。它能限制加工过程中铣削力的波动，改善工件表面质量，提高加工效率。将工件曲面进行离散并划分；建立离散点力的矢量图并计算铣削力；任选一个离散点作为初始点，计算该初始点最近的八个离散点的位置；计算所述的初始点与最近离散点同方向铣削力的差值，并对差值进行比较，得到该初始切削点优化连接方向；将该优化连接方向上的最近离散点作为下一个切削点，计算距离所述的下一个切削点最近的八个离散点的位置；使用前述方法，得到该切削点优化连接方向，重复上一步，直至将所有离散点均连接过为止，得到实际加工中的刀路轨迹。本发明用于自由曲面切削加工。
【专利说明】
-种基于铁削力预测的曲面加工轨迹生成方法
技术领域
[0001] 本发明设及金属切削加工领域，尤其设及一种基于锐削力预测的曲面加工轨迹生 成方法。
【背景技术】
[0002] 在航空、造船、汽车制造等领域中，会经常设及一些如轮机叶片、舰船螺旋奖、汽车 外形等具有特殊过渡曲线的曲面，运种曲面称作自由曲面。自由曲面可满足现代产品的精 美外形及性能要求，但是，由于其对于精度较高的要求，使得自由曲面成为数控加工领域研 究的热点和难点。自由曲面数控编程是其数字化设计与制造集成技术的重要环节，而自由 曲面加工刀具路径轨迹规划技术作为自由曲面数控编程的核屯、技术，具有举足轻重的地 位。针对自由曲面进行刀具路径轨迹规划旨在提高自由曲面的加工质量和加工效率，传统 刀具路径轨迹规划方法都是在只考虑刀具与工件之间几何关系的情况下进行路径规划的， 而在实际切削加工过程中，各种物理因素对加工过程的影响也不容小勵，例如切削力、切削 溫度、切削振动、刀具磨损、不同工艺参数对加工质量的影响及危害等。运些因素对切削加 工的经济性影响重大，在加工过程中考虑物理因素的影响势在必行。
【发明内容】
[0003] 本发明的目的是提供一种可W限制加工过程中锐削力的波动，改善工件的表面质 量，并有效降低加工刀具路径长度，提高曲面加工效率的基于锐削力预测的曲面加工轨迹 生成方法。
[0004] 本发明旨在分析传统刀具路径轨迹规划方法的基础上，弥补传统刀具轨迹规划方 法只考虑刀具与工件相对几何关系，而未曾关注切削加工过程中各种物理因素对加工过程 影响的不足，将锐削力的变化引入到刀具路径轨迹规划中，并将刀具路径轨迹规划问题转 化为TSP(旅行商）问题，提出一种基于锐削力预测的曲面加工轨迹生成新方法，目的是限制 加工过程中锐削力波动，提高曲面加工效率，并且改善工件表面质量。
[0005] 实现上述目的，本发明的技术方案是：
[0006] -种基于锐削力预测的曲面加工轨迹生成方法，所述的方法包括W下步骤：
[0007] 步骤1:将工件的自由曲面进行离散，并根据每个离散点曲率变化划分该自由曲 面；
[000引步骤2:对每个离散点W水平向左为基准按顺时针方向W45°为间隔分别建立八个 进给方向力的矢量图并计算该八个方向的锐削力大小；
[0009] 步骤3:在所有的离散点中任选一个离散点作为初始切削点，计算距离该初始切削 点最近的八个离散点的位置；
[0010] 步骤4:利用遗传算法计算所述的初始切削点与其最近的八个离散点同方向锐削 力的差值，并对锐削力的差值进行比较，得到该初始切削点优化连接方向，并将该初始切削 点与优化连接方向上的最近离散点连接；
[0011] 步骤5:将该优化连接方向上的最近离散点作为下一个切削点，计算距离所述的下 一个切削点最近的八个离散点的位置；
[0012] 步骤6:利用遗传算法计算所述的下一个切削点与其最近的八个离散点同方向锐 削力的差值，并对锐削力的差值进行比较，得到该切削点优化连接方向后，执行步骤5,直至 将所有离散点均连接过为止，得到实际加工过程中应用的刀路轨迹。
[0013] 本发明相对于现有技术的有益效果是：
[0014] 1.本发明将锐削力的变化引入到刀具路径轨迹规划中，解决了传统刀具轨迹规划 方法只考虑刀具与工件相对几何关系，而未曾关注切削加工过程中物理因素对加工过程影 响的问题，限制了加工过程中锐削力的波动，保证了切削过程的平稳。
[0015] 2.本发明在考虑锐削力的变化基础上采用TSP(旅行商）问题来优化刀路轨迹，能 够快速有效地找到最优刀具路径，提高了曲面加工效率，加工的工件表面光顺性更好，改善 了工件的表面质量。
【附图说明】
[0016] 图1为本发明基于锐削力预测的曲面加工轨迹生成方法流程框图；
[0017] 图2为二维曲面离散仿真结果图；
[0018] 图3为特征曲面划分效果图；
[0019] 图4为离散点主要进给方向示意图；
[0020] 图5为离散点整体受力图；
[0021] 图6为图5的局部放大图；
[0022] 图7为刀路轨迹仿真结果图。
【具体实施方式】
[0023] 为了有助于本领域技术人员理解，下面结合附图和相关实施例对本发明方案做进 一步的阐述，W下实施例并不构成对本发明的不当限定。
【具体实施方式】 [0024] 一:如图1所示，本实施方式披露了一种基于锐削力预测的曲面加工 轨迹生成方法，所述的方法包括W下步骤：
[0025] 步骤1:将工件的自由曲面进行离散，并根据每个离散点曲率变化划分该自由曲 面；
[0026] 步骤2:对每个离散点W水平向左为基准按顺时针方向W45°为间隔分别建立八个 进给方向力的矢量图并计算该八个方向的锐削力大小；
[0027] 步骤3:在所有的离散点中任选一个离散点作为初始切削点，计算距离该初始切削 点最近的八个离散点的位置；
[0028] 步骤4:利用遗传算法计算所述的初始切削点与其最近的八个离散点同方向锐削 力的差值，并对锐削力的差值进行比较，得到该初始切削点优化连接方向，并将该初始切削 点与优化连接方向上的最近离散点连接；
[0029] 步骤5:将该优化连接方向上的最近离散点作为下一个切削点，计算距离所述的下 一个切削点最近的八个离散点的位置；
[0030] 步骤6:利用遗传算法计算所述的下一个切削点与其最近的八个离散点同方向锐 削力的差值，并对锐削力的差值进行比较，得到该切削点优化连接方向后，执行步骤5,直至 将所有离散点均连接过为止，得到实际加工过程中应用的刀路轨迹。如图7所示，图中的短 折线代表刀具走过的路径，长线段代表经过曲面划分后的边界线，圆形代表抬刀，=角形代 表进刀。
【具体实施方式】 [0031] 二:一所述的一种基于锐削力预测的曲面加工轨迹生 成方法，步骤1中，将工件的自由曲面进行离散，并根据每个离散点曲率变化划分该自由曲 面的具体方法步骤是：
[0032] 步骤1-1:将所述的工件自由曲面看做是由多个简单曲面拼凑而成的，每个所述的 简单曲面的形式描述为P(u，v)，其中P(u，v)是曲面的参数方程，u、v是满足0《u、v《l的两 个参数；
[0033] 本实施方式中，所选曲面参数方程如下：
[0034]
[0035] 所选曲面实例的MATLAB二维仿真离散结果如图2所示；
[0036] 步骤1-2:根据每个离散点的曲率不同将工件的自由曲面整体划分为凸曲面、凹曲 面和马鞍面。
[0037] 工件的自由曲面划分完成后，利用MATLAB软件对所选的曲面（即凸曲面、凹曲面或 马鞍面)进行曲面细分的仿真。图3为实现的曲面分片效果，其中圆圈代表凹面、五角星代表 马鞍面、立角形代表凸曲面，图内的直线代表经过曲面划分后的边界线。
【具体实施方式】 [0038] 一所述的一种基于锐削力预测的曲面加工轨迹生 成方法，步骤2中，对每个离散点W水平向左为基准按顺时针方向W45°为间隔分别建立八 个进给方向力的矢量图并计算八个方向的锐削力大小的具体方法步骤是：
[0039] 步骤2-1:将每个简单曲面近似看做为任意进给方向的斜平面微元加工，对运个简 单曲面上的离散点建立W水平向左为基准按顺时针方向W 45°为间隔八个进给方向力的矢 量图，如图4所示；
[0040] 步骤2-2:由于微元锐削力计算公式是在水平面上建立的，而斜平面和平面之间的 转换只需进行相应的坐标变换就能实现，计算八个进给方向的锐削力大小所用锐削力模型 为下式：
[0041]
123 其中，j-表示第j个刀齿； 2 化-表示锐刀齿数； 3 O-表示锐刀旋转角度；
[0045] 0 -表示球头刃线处任意一点的位置角。
[0046] 仿真结果如图5和图6所示，图中圆内线段所指的方向代表该离散点受力方向，线 段的长短代表该离散点在线段所指方向上受力的大小，线段越长代表该离散点所受的该方 向的力越大，反之，线段越短代表该离散点所受的该方向的力越小。
【具体实施方式】 [0047] 四：一所述的一种基于锐削力预测的曲面加工轨迹生 成方法，步骤3中，在所有的离散点中任选一个离散点作为初始切削点，计算距离该初始切 削点最近的八个离散点的位置的具体方法步骤是：
[0048] 步骤3-1:将刀具(锐刀）路径优化问题可归结为TSP(旅行商）问题，其中刀具相当 于旅行商，每个刀位点相当于旅行商要经过的城市，而最优走刀路径的目标函数即为所有 路径长度最短；
[0049] 步骤3-2:对路径长度最短的计算通过W下方式实现：
[0050] 若对于城市W= (vi，V2,…乂。)的一个访问顺序为(ti，t2,…心），其中i G V= U， 2,…，n}，且TSP(旅行商）问题的数学模型为下式：
[0化1 ]
[0052] 其中，W-表示城市列表；
[0053] t广表示第i个访问城市；
[0054] d(titw)-表示城市ti和城市心込间的距离；
[0化日]k表示通过所有城市且每个城市仅通过一次时的距离。
【具体实施方式】 [0056] 五:四所述的一种基于锐削力预测的曲面加工轨迹生 成方法，步骤4中，利用遗传算法计算所述的初始切削点与其最近的八个离散点同方向锐削 力的差值，并对锐削力的差值进行比较，得到该初始切削点优化连接方向的具体方法为：
[0057] 利用遗传算法计算所述的初始切削点八个方向的每个同方向力差值的大小，计算 初始切削点与其最近八个离散点每个同方向锐削力的差值，并对初始切削点与其最近八个 离散点同方向锐削力的差值进行比较，选择八个切削力差值中最小的切削力的差值方向作 为优化连接方向。
[0058] 遗传算法的计算原理可参考文献（1.卫蔵，李建勇，王恒.基于遗传算法的PCB数控 钻孔路径优化[J].计算机工程与应用，2008，25:229-232.2.李明海，邢桂华.用MATLAB实现 中国旅行商问题的求解[J].微计算机应用，2004,02:218-222.);
[0059] 虽然本发明已W较佳实施例披露如上，然而并非用W限定本发明的，本领域技术 人员还可W在本发明精神内做其他变化，W及应用到本发明未提及的领域中，当然，运些依 据本发明精神所做的变化都应包含在本发明所要求保护的范围内。
【主权项】
1. 一种基于铣削力预测的曲面加工轨迹生成方法，其特征在于:所述的方法包括以下 步骤： 步骤1:将工件的自由曲面进行离散，并根据每个离散点曲率变化划分该自由曲面； 步骤2:对每个离散点以水平向左为基准按顺时针方向以45°为间隔分别建立八个进给 方向力的矢量图并计算该+八个方向的铣削力大小； 步骤3:在所有的离散点中任选一个离散点作为初始切削点，计算距离该初始切削点最 近的八个离散点的位置； 步骤4:利用遗传算法计算所述的初始切削点与其最近的八个离散点同方向铣削力的 差值，并对铣削力的差值进行比较，得到该初始切削点优化连接方向，并将该初始切削点与 优化连接方向上的最近离散点连接； 步骤5:将该优化连接方向上的最近离散点作为下一个切削点，计算距离所述的下一个 切削点最近的八个离散点的位置； 步骤6:利用遗传算法计算所述的下一个切削点与其最近的八个离散点同方向铣削力 的差值，并对铣削力的差值进行比较，得到该切削点优化连接方向后，执行步骤5,直至将所 有离散点均连接过为止，得到实际加工过程中应用的刀路轨迹。2. 根据权利要求1所述的一种基于铣削力预测的曲面加工轨迹生成方法，其特征在于： 步骤1中，将工件的自由曲面进行离散，并根据每个离散点曲率变化划分该自由曲面的具体 方法步骤是： 步骤1-1:将所述的工件自由曲面看做是由多个简单曲面拼凑而成的，每个所述的简单 曲面的形式描述为P(u，v)，其中P(u，v)是曲面的参数方程，u、v是满足0<U、V<1的两个参 数； 步骤1-2:根据每个离散点的曲率不同将工件的自由曲面整体划分为凸曲面、凹曲面和 马鞍面。3. 根据权利要求1所述的一种基于铣削力预测的曲面加工轨迹生成方法，其特征在于： 步骤2中，对每个离散点以水平向左为基准按顺时针方向以45°为间隔分别建立八个进 给方向力的矢量图并计算八个方向的铣削力大小的具体方法步骤是： 步骤2-1:将每个简单曲面近似看做为任意进给方向的斜平面微元加工，对这个简单曲 面上的离散点建立以水平向左为基准按顺时针方向以45°为间隔八个进给方向力的矢量 图； 步骤2-2:由于微元铣削力计算公式是在水平面上建立的，而斜平面和平面之间的转换 只需进行相应的坐标变换就能实现，计算八个进给方向的铣削力大小所用铣削力模型为下 式： 其中，j_表示第j个刀齿；Nf-表示铣刀齿数； Φ-表示铣刀旋转角度； ? -表示球头刃线处任意一点的位置角。4. 根据权利要求1所述的一种基于铣削力预测的曲面加工轨迹生成方法，其特征在于： 步骤3中，在所有的离散点中任选一个离散点作为初始切削点，计算距离该初始切削点最近 的八个离散点的位置的具体方法步骤是： 步骤3-1:将刀具(铣刀）路径优化问题可归结为TSP(旅行商）问题，其中刀具相当于旅 行商，每个刀位点相当于旅行商要经过的城市，而最优走刀路径的目标函数即为所有路径 长度最短； 步骤3-2:对路径长度最短的计算通过以下方式实现： 若对于城市w= ( VI，V2，…vn)的一个访问顺序为T彡（ti，t2，"_tn)A*ieV={l，2，···， η}，且TSP(旅行商）问题的数学模型为下式： 其中，W-表不城市列表；ti-表示第i个访问城市； d(titi+i)_表示城市ti和城市ti+i之间的距离； L-表示通过所有城市且每个城市仅通过一次时的距离。5. 根据权利要求4所述的一种基于铣削力预测的曲面加工轨迹生成方法，其特征在于： 步骤4中，利用遗传算法计算所述的初始切削点与其最近的八个离散点同方向铣削力的差 值，并对铣削力的差值进行比较，得到该初始切削点优化连接方向的具体方法为： 利用遗传算法计算所述的初始切削点八个方向的每个同方向力差值的大小，计算初始 切削点与其最近八个离散点每个同方向铣削力的差值，并对初始切削点与其最近八个离散 点同方向铣削力的差值进行比较，选择八个切削力差值中最小的切削力的差值方向作为优 化连接方向。
【文档编号】G05B19/19GK105955195SQ201610321865
【公开日】2016年9月21日
【申请日】2016年5月16日
【发明人】吴明阳, 陈勇, 程耀楠, 张芮, 赵旭
【申请人】哈尔滨理工大学