CN107045578B - 一种基于npso算法的容屑槽加工砂轮位姿求解方法 - Google Patents

一种基于npso算法的容屑槽加工砂轮位姿求解方法 Download PDF

Info

Publication number
CN107045578B
CN107045578B CN201710276576.4A CN201710276576A CN107045578B CN 107045578 B CN107045578 B CN 107045578B CN 201710276576 A CN201710276576 A CN 201710276576A CN 107045578 B CN107045578 B CN 107045578B
Authority
CN
China
Prior art keywords
grinding wheel
chip pocket
pose
algorithm
npso
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
Application number
CN201710276576.4A
Other languages
English (en)
Other versions
CN107045578A (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.)
KUNSHAN YOUYUANSHENG MACHINERY Co.,Ltd.
Original Assignee
Jiangsu University of Science and Technology
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 Jiangsu University of Science and Technology filed Critical Jiangsu University of Science and Technology
Priority to CN201710276576.4A priority Critical patent/CN107045578B/zh
Publication of CN107045578A publication Critical patent/CN107045578A/zh
Application granted granted Critical
Publication of CN107045578B publication Critical patent/CN107045578B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design

Landscapes

  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Polishing Bodies And Polishing Tools (AREA)

Abstract

本发明提出了一种基于NPSO算法的容屑槽加工砂轮位姿求解方法，属于数控刀具制造领域，用于整体式刀具容屑槽刃磨工艺的制定。首先，基于容屑槽不同刃磨精度需求，建立了描述砂轮位姿是否合适的目标函数；其次，根据实际经验设置了砂轮初始位姿参数；然后，根据砂轮位姿对容屑槽结构的影响程度，设置了砂轮位姿更新速度；最后通过NPSO算法求解获得砂轮位姿参数。本发明方法在已知砂轮形状和容屑槽形状的条件下，实现砂轮位姿参数的快速、精确求解，为整体刀具容屑槽刃磨工艺的制定提供技术和方法支撑，经实际验证，该方法具有较高的鲁棒性和可靠性。

Description

一种基于NPSO算法的容屑槽加工砂轮位姿求解方法
技术领域
本发明涉及整体刀具制造领域，具体为一种基于NPSO算法的容屑槽加工砂轮位姿求解方法。
背景技术
容屑槽是整体立铣刀、钻头等刀具关键结构之一，不同的被加工材料及结构特征需设计与之相匹配的刀具容屑槽形状，形状尺寸各异的容屑槽不断涌现，高效高精度的容屑槽制造工艺成为整体数控刀具制造的核心技术之一。容屑槽通过砂轮刃磨成形，主要刃磨策略可分为两类：其一，设计与容屑槽相匹配的成形砂轮，采用成形砂轮生产制造所需容屑槽，其二，采用标准砂轮，通过五轴联动数控工件磨床调整成形砂轮位姿，制造所需容屑槽。成形砂轮制造精度高，但同时具有高成本、高周期的特点。根据砂轮形状合理设置砂轮位姿，进而获得容屑槽设计形状是降低整体刀具生产成本和周期的重要方法。现有砂轮位姿求解方法通常采用循环迭代思想，较难保证求解精度，且对复杂形状容屑槽适应性较差。为此，本发明提出了一种基于NPSO的砂轮位姿求解智能算法。
发明内容
发明目的：为解决基于标准砂轮的容屑槽刃磨砂轮位姿求解难题，本发明提出了一种基于NPSO算法的容屑槽加工砂轮位姿求解方法。
技术方案：一种基于NPSO算法的容屑槽加工砂轮位姿求解方法，包括以下步骤：
(1)建立描述砂轮位姿参数是否合适的目标函数；
(2)设置砂轮初始位姿参数及位姿更新速度；
(3)通过NPSO算法对砂轮位姿参数进行求解。
步骤(1)中的目标函数为：
其中，γ0、φ、rc分别为目标容屑槽的前角、槽宽、芯径，xi为砂轮位姿参数矩阵，fγ0(xi)、fφ(xi)、分别为砂轮在位姿xi时加工获得容屑槽的前角、槽宽、芯径，fdist(j)为描述容屑槽的第j个计算点与目标点之间的距离，ncore为芯径点序号，n为容屑槽端截面线点的个数，u()为罚函数，α1,α2,α3,α4,α5为协调系数。
进一步地，砂轮位姿参数矩阵xi＝[Δαx,Δx,Δy,Δz]，其中，Δαx、Δx、Δy、Δz分别为砂轮从刀具坐标系原点处绕刀具坐标系XT轴旋转角度以及沿刀具坐标轴XT、YT、ZT移动距离。
进一步地，步骤(2)中砂轮初始位姿参数设置为：
其中，Δαx取值为容屑槽螺旋角减5°，Δx取值为砂轮半径加rc/2，Δy取值为-rc/2，Δz取值为容屑槽导程值的一半，rand()为[0,1]范围内的随机数，k为迭代次数，np为种群中粒子个数。
步骤(2)中砂轮初始位姿更新速度为：
其中，[vΔαx,vΔx,vΔy,vΔz]T＝[1,1,1,1]T，vΔαx、vΔx、vΔy、vΔz分别为Δαx、Δx、Δy、Δz对应的更新速度。
进一步地，步骤(3)的NPSO算法中，每轮迭代计算并评估每个粒子的适应度值更新当前个体最优解和邻域最优解i＝1时，i-1取np，i＝np时，i+1取1。
进一步地，步骤(3)的NPSO算法中，根据如下公式更新每个粒子速度和位置：
式中，ω为惯性权重，λ1和λ2为加速常数。
进一步地，采用线性减小惯性权重策略，ω(k)＝ωmax-(ωmax-ωmin)·k/nk，其中ωmax为最大惯性权重，ωmin为最小惯性权重，nk为最大迭代次数。
有益效果：本发明提出的方法具有较好的通用性和求解精度，适合立铣刀、钻头等整体式刀具容屑槽，容屑槽的形状不受限制，且可以根据不同的刃磨精度要求求解相应的砂轮位姿策略，采用环形小生境粒子群智能算法，根据砂轮位姿对容屑槽形状的影响规律设置算法参数，增加了算法的收敛速度及全局最优解的搜索能力。使用该方法在已知砂轮形状和容屑槽形状的条件下，实现砂轮位姿参数的快速、精确求解，为整体刀具容屑槽刃磨工艺的制定提供技术和方法支撑，经实际验证，该方法具有较高的鲁棒性和可靠性。
附图说明
图1为本发明实施例的方法流程图；
图2为本发明实施例中的砂轮形状示意图；
图3为本发明实施例中的容屑槽形状示意图；
图4为不同迭代次数获得容屑槽形状示意图，图(a)～(e)分别对应迭代次数k为1、50、70、100、150时的容屑槽形状。
具体实施方式
本发明实施例公开的一种基于NPSO算法的容屑槽加工砂轮位姿求解方法，在建立砂轮位姿参数与容屑槽结构参数之间关系的基础上，基于环形拓扑小生境粒子群优化算法(NPSO)求得砂轮位姿参数。首先，基于容屑槽不同刃磨精度需求，建立了描述砂轮位姿是否合适的目标函数；其次，根据实际经验设置了砂轮初始位姿参数；然后，根据砂轮位姿对容屑槽结构的影响程度，设置了砂轮位姿更新速度；最后通过NPSO算法求解获得砂轮位姿参数。如图1所示，砂轮位姿求解过程具体包括以下步骤：
步骤1：建立目标函数
其中，γ0、φ、rc分别为目标容屑槽的前角、槽宽、芯径，xi为代表砂轮位姿参数的4×1阶矩阵，其值为xi＝[Δαx,Δx,Δy,Δz]，Δαx、Δx、Δy、Δz分别为砂轮从刀具坐标系原点处绕刀具坐标系XT轴旋转角度以及沿刀具坐标轴XT、YT、ZT移动距离，fγ0(xi)、fφ(xi)、分别为砂轮在位姿xi时加工获得容屑槽参数值，fdist(j)为描述容屑槽的第j个计算点与目标点之间的距离，ncore为芯径点序号，n为描述容屑槽端截面线点的个数，取值范围为[10,50]，u()为罚函数，当frc(xi)大于刀具半径时u()取值为10000，α1,α2,α3,α4,α5为满足前角、齿间距、芯径、前刀面轮廓、后刀面轮廓不同设计需求而设定的协调系数，满足α1+α2+α3+α4+α5＝1；
步骤2：初始化粒子位置和速度，其中初始种群位置为
其中，Δαx取值为容屑槽螺旋角减5°，Δx取值为砂轮半径加rc/2，Δy取值为-rc/2，Δz取值为容屑槽导程值的一半，rand()为[0,1]范围内的随机数，k为迭代次数，np为种群中粒子个数，取值范围为[30,200]。
初始种群速度为：
其中，[vΔαx,vΔx,vΔy,vΔz]T＝[1,1,1,1]T，vΔαx、vΔx、vΔy、vΔz分别为Δαx、Δx、Δy、Δz对应的速度；
步骤3：计算并评估每个粒子的适应度值更新当前个体最优解和邻域最优解pBestk 3,i，
其中，为第i个粒子在第k次迭代后对应的适应度函数；为粒子i在1～k次迭代中最小的适应度值，即当前个体最优解；pBestk 3,i为在第k次迭代后粒子i的邻域最优解，根据环形拓扑结构，即为pBesti-1 k，和pBesti+1 k三者最优；
步骤4：根据自身历史最优位置和邻域最优位置，更新每个粒子速度和位置：
式中，ω为惯性权重，拟采用线性减小惯性权重策略，即ω(k)＝ωmax-(ωmax-ωmin)·k/nk，式中ωmax为最大惯性权重，取值0.9，ωmin为最小惯性权重，取值0.4，nk为最大迭代次数，取值大于200，k为已迭代次数；λ1和λ2为加速常数，根据砂轮位姿各参数对容屑槽形状的不同影响程度，取值为λ1＝λ2＝1.494·[1.5，0.5，1，1]；
步骤5：获取全局最优解，比较粒子对应的邻域最优解，其中最小值为全局最优解，即：
其中，gBest＝[Δαx′,Δx′,Δy′,Δz′]即为加工目标容屑槽的砂轮位姿参数。
下面结合具体实例进一步描述本发明：
本实施实例选择双斜面型砂轮，具体尺寸为gR＝75mm，gb＝20mm，gb1＝5mm，ga1＝10°，ga2＝70°，gr1＝1mm，gr2＝1mm，gr3＝1mm，如图2所示。需要加工的容屑槽端截面参数为γ0＝6.41°，φ＝88.64°，rc＝5.01mm，直径为20mm，螺旋角为45°，导程为62.8mm，形状如图3所示，描述容屑槽端截面线各点的坐标为：
步骤1：建立目标函数
其中，u()为罚函数，当大于10时u()取值为10000。
步骤2：初始化粒子位置和速度，其中初始种群位置为
其中，rand()为[0,1]范围内的随机数，k为迭代次数，np为种群中粒子个数，取值为50。
初始种群速度为：
步骤3：计算并评估每个粒子的适应度值更新当前个体最优解和邻域最优解pBestk 3,i，
其中，为第i个粒子在第k次迭代后对应的适应度函数；为粒子i在1～k次迭代中最小的适应度值，即当前个体最优解；pBestk 3,i为在第k次迭代后粒子i的邻域最优解，根据环形拓扑结构，即为pBesti-1 k，和pBesti+1 k三者最优。
步骤4：根据自身历史最优位置和邻域最优位置，更新每个粒子速度和位置。
式中，ω为惯性权重，拟采用线性减小惯性权重策略，即ω(k)＝ωmax-(ωmax-ωmin)·k/nk，式中ωmax为最大惯性权重，取值0.9，ωmin为最小惯性权重，取值0.4，nk为最大迭代次数，取值250，k为已迭代次数；λ1和λ2为加速常数，根据砂轮位姿各参数对容屑槽形状的不同影响程度，取值为λ1＝λ2＝1.494·[1.5，0.5，1，1]。
步骤5：获取全局最优解，比较粒子对应的邻域最优解，其中最小值为全局最优解，即：
求得gBest＝[0.628,77,-15,20]即为加工目标容屑槽的砂轮位姿参数，随着迭代次数增加，获得容屑槽形状变化过程如图4所示，可知循环150次获得砂轮位姿能够满足要求。

Claims (1)

1.一种基于NPSO算法的容屑槽加工砂轮位姿求解方法，其特征在于，包括以下步骤：
(1)建立描述砂轮位姿参数是否合适的目标函数；
(2)设置砂轮初始位姿参数及位姿更新速度；
(3)通过NPSO算法对砂轮位姿参数进行求解；
所述步骤(1)中的目标函数为：
其中，γ0、φ、rc分别为目标容屑槽的前角、槽宽、芯径，xi为砂轮位姿参数矩阵，fγ0(xi)、fφ(xi)、分别为砂轮在位姿xi时加工获得容屑槽的前角、槽宽、芯径，fdist(j)为描述容屑槽的第j个计算点与目标点之间的距离，ncore为芯径点序号，n为容屑槽端截面线点的个数，u()为罚函数，α1,α2,α3,α4,α5为协调系数；
砂轮位姿参数矩阵xi＝[Δαx,Δx,Δy,Δz]，其中，Δαx、Δx、Δy、Δz分别为砂轮从刀具坐标系原点处绕刀具坐标系XT轴旋转角度以及沿刀具坐标轴XT、YT、ZT移动距离；
步骤(2)中砂轮初始位姿参数设置为：
其中，Δαx取值为容屑槽螺旋角减5°，Δx取值为砂轮半径加rc/2，Δy取值为-rc/2，Δz取值为容屑槽导程值的一半，rand()为[0,1]范围内的随机数，k为迭代次数，np为种群中粒子个数；
步骤(2)中砂轮初始位姿更新速度为：
其中，[vΔαx,vΔx,vΔy,vΔz]T＝[1,1,1,1]T，vΔαx、vΔx、vΔy、vΔz分别为Δαx、Δx、Δy、Δz对应的更新速度；
步骤(3)的NPSO算法中，每轮迭代计算并评估每个粒子的适应度值fobj k,i(xi)，更新当前个体最优解和邻域最优解其中，k为迭代次数，i＝1,…,np，np为种群中粒子个数，i＝1时，i-1取np，i＝np时，i+1取1；
步骤(3)的NPSO算法中，根据如下公式更新每个粒子速度和位置：
式中，ω为惯性权重，k为迭代次数，λ1和λ2为加速常数；
采用线性减小惯性权重策略，ω(k)＝ωmax-(ωmax-ωmin)·k/nk，其中ωmax为最大惯性权重，ωmin为最小惯性权重，nk为最大迭代次数。
CN201710276576.4A 2017-04-25 2017-04-25 一种基于npso算法的容屑槽加工砂轮位姿求解方法 Active CN107045578B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710276576.4A CN107045578B (zh) 2017-04-25 2017-04-25 一种基于npso算法的容屑槽加工砂轮位姿求解方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710276576.4A CN107045578B (zh) 2017-04-25 2017-04-25 一种基于npso算法的容屑槽加工砂轮位姿求解方法

Publications (2)

Publication Number Publication Date
CN107045578A CN107045578A (zh) 2017-08-15
CN107045578B true CN107045578B (zh) 2018-12-21

Family

ID=59545780

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710276576.4A Active CN107045578B (zh) 2017-04-25 2017-04-25 一种基于npso算法的容屑槽加工砂轮位姿求解方法

Country Status (1)

Country Link
CN (1) CN107045578B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108098515B (zh) * 2017-12-12 2020-07-28 科德数控股份有限公司 一种使用多种成型砂轮加工钻头槽型的方法
CN110990966B (zh) * 2019-11-08 2023-08-22 江苏科技大学 一种刀具复杂容屑槽磨制砂轮轨迹确定方法
CN113642220B (zh) * 2021-08-26 2023-09-22 江苏科技大学 一种基于rbf与mopso的船舶焊接工艺优化方法
CN114048565A (zh) * 2021-11-11 2022-02-15 江苏科技大学 一种基于机器学习及大数据的整体刀具刃磨工艺求解方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102183205A (zh) * 2011-01-19 2011-09-14 北京航空航天大学 一种大型零部件最佳装配位姿匹配的方法
CN104741994A (zh) * 2015-03-25 2015-07-01 华南理工大学 一种任意曲面砂轮用于曲面精密磨削的方法
CN105676776A (zh) * 2016-01-15 2016-06-15 广东工业大学 基于粒子群算法的轮廓曲线数控代码生成方法及数控机床

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI414376B (zh) * 2011-04-14 2013-11-11 Nat Univ Tsing Hua 一種五軸曲面側銑加工系統及其刀具路徑規劃方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102183205A (zh) * 2011-01-19 2011-09-14 北京航空航天大学 一种大型零部件最佳装配位姿匹配的方法
CN104741994A (zh) * 2015-03-25 2015-07-01 华南理工大学 一种任意曲面砂轮用于曲面精密磨削的方法
CN105676776A (zh) * 2016-01-15 2016-06-15 广东工业大学 基于粒子群算法的轮廓曲线数控代码生成方法及数控机床

Also Published As

Publication number Publication date
CN107045578A (zh) 2017-08-15

Similar Documents

Publication Publication Date Title
CN107045578B (zh) 一种基于npso算法的容屑槽加工砂轮位姿求解方法
CN110052713B (zh) 零件增减材复合制造工艺
CN103128312B (zh) 一种薄壁盘零件控制变形加工工艺
CN104015016A (zh) 高精度薄壁深腔零件的加工方法
CN103454973B (zh) 参数化数控立铣刀螺旋槽数控加工编程方法
CN106125666A (zh) 以切削力波动为约束的曲面加工刀具轨迹规划方法
CN102806380A (zh) 开式整体叶盘通道复合粗加工方法
CN102019572A (zh) 采用组合螺旋式抛光路径的抛光工艺
CN105269052A (zh) 精密凸轮轴加工方法
CN103645674A (zh) 一种整体叶轮叶片的粗-半精-精铣混合路径生成方法
CN103586519A (zh) 梯型槽分层铣削粗加工方法
CN107505913B (zh) 基于整体叶盘通道四轴数控加工的最大适用刀具半径计算方法
CN202411558U (zh) 粗精一体加工的金刚石铰刀
CN103163837B (zh) 一种面向五轴数控加工的刀具轨迹自适应管理系统及方法
CN103962807B (zh) 一种滑枕柔性加工的方法
CN103111529B (zh) U型钣金零件回弹模胎的加工方法
CN104439468A (zh) 适用于整体叶盘结构型面分层铣削成型工艺的铣刀
CN102689164A (zh) 大长径比钛合金环槽铆钉的加工工艺
CN103831600A (zh) 一种多角度孔钻模加工的方法
CN110039406A (zh) 一种单晶硅光学复杂表面的超精密加工工具及加工方法
CN103419540A (zh) 基于路径单元的变形曲面自适应投影加工方法
CN103084809A (zh) 一种三棱孔加工方法
CN101587348A (zh) 三维刃口轮廓的螺旋线插补加工方法
CN101905489B (zh) 金刚石薄膜涂层刀具对精密细小石墨电极的加工方法
CN106862628A (zh) 数控机床主轴定向自制插刀对型腔内直角加工方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20170815

Assignee: Zhenjiang Kaituo Machinery Co., Ltd

Assignor: JIANGSU University OF SCIENCE AND TECHNOLOGY

Contract record no.: X2020980007284

Denomination of invention: A method of solving the grinding wheel's pose in chip groove machining based on NPSO algorithm

Granted publication date: 20181221

License type: Common License

Record date: 20201029

EE01 Entry into force of recordation of patent licensing contract
EC01 Cancellation of recordation of patent licensing contract

Assignee: Zhenjiang Kaituo Machinery Co., Ltd

Assignor: JIANGSU University OF SCIENCE AND TECHNOLOGY

Contract record no.: X2020980007284

Date of cancellation: 20201223

EC01 Cancellation of recordation of patent licensing contract
TR01 Transfer of patent right

Effective date of registration: 20211220

Address after: 215300 room 7, 1188 Hongqiao Road, Yushan Town, Kunshan City, Suzhou City, Jiangsu Province

Patentee after: KUNSHAN YOUYUANSHENG MACHINERY Co.,Ltd.

Address before: 212003, No. 2, Mengxi Road, Zhenjiang, Jiangsu

Patentee before: JIANGSU University OF SCIENCE AND TECHNOLOGY

TR01 Transfer of patent right