基于svr-pso的超声加工特种刀具优化设计方法
【专利摘要】本发明公开了一种基于SVR(支持向量回归)?PSO(粒子群算法)的超声加工特种刀具优化设计方法。本发明通过建立结构响应的近似模型，定量的建立刀具结构参数与结构响应的映射关系，为更加精确的对刀具结构进行优化设计提供了解决方法。
【专利说明】
基于SVR-PSO的超声加工特种刀具优化设计方法
技术领域
[0001] 本发明设及一种基于SVR(支持向量回归)-PSO(粒子群算法)的超声加工特种刀具 优化设计方法。
【背景技术】
[0002] 超声加工技术是特种加工技术之一，针对传统加工技术难W加工的复合材料，超 声加工技术往往能取得意想不到的效果。蜂窝复合材料因其结构特点和组分特点，是一种 典型的难加工材料。超声加工技术的应用有效地解决了蜂窝复合材料传统锐削加工中存在 的环境污染严重、加工表面质量差等问题。在超声加工组件中，超声加工特种刀具是极其重 要的一部分。刀具的结构特征和声学性能直接影响到蜂窝复合材料的表面加工质量和加工 效率，W及刀具本身的使用寿命。通过大量文献的阅读发现，目前学者们对超声加工刀具的 研究广泛地应用了有限元法和实验法，前者定性地总结了刀具结构参数对结构响应的变化 规律，为刀具参数优化提供了思路，后者通过实物加工应用，确定合理的结构参数。但由于 有限元法的时间成本和实验法的制作成本所限，难W获取刀具结构参数与结构响应的量化 关系，运往往使得最终的设计结果难W达到最优。
【发明内容】
[0003] 为了解决现有超声加工特种刀具设计方法中存在高时间成本和高制作成本，W及 难W量化刀具结构参数与结构响应关系的问题，本发明提供一种更加快捷、精确的超声加 工特种刀具结构的优化设计方法。
[0004] 为了解决上述技术问题，本发明采用如下技术方案:基于SVR-PSO的超声加工特种 刀具优化设计方法，包括W下步骤：
[0005] 步骤一:利用ABAQUS有限元二次开发技术建立刀具参数优化插件，利用此插件获 取一定数量的刀具结构参数变量和对应的结构响应变量的样本数据；
[0006] 步骤二：W上述样本数据为基础，运用支持向量回归方法得到刀具结构参数与结 构响应之间映射关系的近似模型；
[0007] 步骤=:改进粒子群算法，提升粒子群算法在优化过程中的收敛速度；
[0008] 步骤四：建立刀具结构设计的优化目标函数；
[0009] 步骤五:应用改进的粒子群算法对刀具结构进行优化，获取刀具设计的最优参数。
[0010] 所述步骤二是运用高斯径向基函数作为支持向量机回归模型的核函数，利用 巧thon平台下的SC化it-learn模块的svr函数对获取的样本数据进行训练得到近似模型。
[0011] 所述步骤=中改进后粒子群的实现公式：
[0012]
[0013]
[0014]
[0015] 式中代表第i个粒子第k次迭代时的位置参数，v代表此时的速度参数，1>^代表 第i个粒子第k次迭代时的速度参数，CO代表的是惯性权重，惯性权重数值上采用惯性权重 线性递减策略，pk'g代表历史全局最优解，乂代表粒子i迭代至第k次的历史最优解，ri及n 是范围[0，1 ]的随机数W使得每次迭代输入值该变量具有一定的随机性，Cl及C2分别代表经 验学习因子与社会学习因子。
[0016] 所述步骤四中优化目标函数为：
[0017] Y(l，d，m) = l-[0.30(l，d，m)+0.3(I-Uid，d，m))+0.3U2(l，d，m)
[001 引 +0.l(l-f(l，d，m))]
[0019] 式中I为刀具长度，m为材料代号，d为刀具厚度，〇(l，d，m)为最大应力，Ui(l，d，m) 为刀尖振幅，化(l，d，m)为刀尖横向位移，f(l，d，m)为固有频率。
[0020] 本发明通过建立结构响应的近似模型，定量的建立刀具结构参数与结构响应的映 射关系，为更加精确的对刀具结构进行优化设计提供了解决方法。
【附图说明】
[0021] 图1是训练样本的最大应力数据图；
[0022] 图2是训练样本的刀尖振幅数据图；
[0023] 图3是训练样本的刀尖横向位移数据图；
[0024] 图4是训练样本的固有频率数据图；
[0025] 图5是刀具优化设计实现流程图。
【具体实施方式】
[0026] 本发明基于SVR-PSO的超声加工特种刀具优化设计方法，包括W下步骤：
[0027] 步骤一:输入变量为刀具长度、刀具材料、刀具厚度，输出变量为最大应力、谐振频 率、刀尖振幅、刀尖横向位移。刀具长度变量选择25-55mm范围间隔2mm的15个，厚度为 1.5111111、2111111两种，材料选择￥￥2硬质合金与高速钢两种材料。利用有限元二次开发建立的超 声切割刀具组件仿真插件对刀具进行参数化遍历仿真获取了最大应力、刀尖振幅、刀尖横 向位移、固有频率等训练样本。对样本中不同维度的数据进行归一化处理，使得数据标准化 在[0，1]或[-1,1]范围内。
[0028] 步骤二：运用高斯径向基函数（RBF)作为支持向量机回归模型的核函数，利用 Python的SC化it-learn模块的svr函数对获取的样本数据进行训练得到近似模型。由对2mm 厚度材料为硬质合金的刀具最大应力的拟合值可知近似模型的样本点拟合误差不超过 0.1%。
[0029] 步骤针对本发明中的刀具而言，设计精度定为0.1mm已经能满足实际情况。因 此，在算法中粒子的随机生成值设定为0.1的整数倍，如粒子群算法中的飞行速度参数Vi取 0.1的整数倍，在其小于0.1时，将其飞行速度重设为最小值0.1。改进后粒子群的实现公式：
[0030] 獄
[0031] (6)
[0032] 式中，4代表第i个粒子第k次迭代时的位置参数，V代表此时的速度参数，皆代表 第i个粒子第k次迭代时的速度参数，CO代表的是惯性权重，惯性权重数值上采用惯性权重 线性递减策略，pk'g代表历史全局最优解，pf代表粒子i迭代至第k次的历史最优解，ri及n 是范围[0，1 ]的随机数W使得每次迭代输入值该变量具有一定的随机性，Cl及C2分别代表经 验学习因子与社会学习因子。经过改进后的粒子群算法较标准粒子群算法而言，收敛所需 的迭代次数要少于标准粒子群算法，提升了粒子群算法在优化过程中的收敛速度。
[0033] 步骤四：对有限元仿真获取的样本的输出量进行数据分析并确定目标函数如式 (2)
[0034] Y(l，d，m) = l-[0.30(l，d，m)+0.3(I-Uid，d，m))+0.3U2(l，d，m)
[0035] +〇.l(l-f(l,d,m))] (2)
[0036] 式中I为刀具长度，m为材料代号(0代表硬质合金，I代表高速钢），d为刀具厚度，曰 (l，d，m)为最大应力，Ui(l，d，m)为刀尖振幅，化（l，d，m)为刀尖横向位移，f(l，d，m)为固有频 率。输出函数都是[0，1]标准数据。所WY(l，d，m)的值为[0，1]内的无量纲常数。
[0037] 则优化目标函数为：
[00；3 引
[0039] (4)
[0040] 步骤五:在本发明实际问题中输入变量有=个参数，仅有刀具长度为连续变量，材 料与厚度均为离散变量，通过将=个结构参数映射到设定的一维连续变量空间（25,145)上 W解决离散变量优化问题，将此变量空间分为长度为30的四部分W对应4类刀具长度变化 范围（25,55)及变化长度，即（l，d，m)^x。运用支持向量机与粒子群算法组合优化设计方法 (SVR-POS)在一维连续变量空间对目标函数进行寻优，寻优后得出的结果与刀具参数相对 应的结果即为刀具结构参数最优值。
【主权项】
1. 基于SVR-PSO的超声加工特种刀具优化设计方法，其特征在于:包括以下步骤： 步骤一：利用ABAQUS有限元二次开发技术建立刀具参数优化插件，利用此插件获取一 定数量的刀具结构参数变量和对应的结构响应变量的样本数据； 步骤二：以上述样本数据为基础，运用支持向量回归方法得到刀具结构参数与结构响 应之间映射关系的近似模型； 步骤三:改进粒子群算法，提升粒子群算法在优化过程中的收敛速度； 步骤四：建立刀具结构设计的优化目标函数； 步骤五:应用改进的粒子群算法对刀具结构进行优化，获取刀具设计的最优参数。2. 根据权利要求1所述的基于SVR-PS0的超声加工特种刀具优化设计方法，其特征在 于:所述步骤二是运用高斯径向基函数作为支持向量机回归模型的核函数，利用Python的 scikit-learn模块的svr函数对获取的样本数据进行训练得到近似模型。3. 根据权利要求1所述的基于SVR-PS0的超声加工特种刀具优化设计方法，其特征在 于:所述步骤三中改进后粒子群的实现公式：式中，< 代表第i个粒子第k次迭代时的位置参数，v代表此时的速度参数，vf代表第i个 粒子第k次迭代时的速度参数，ω代表的是惯性权重，惯性权重数值上采用惯性权重线性递 减策略，Pk' g代表历史全局最优解'代表粒子i迭代至第k次的历史最优解，rLSr2是范围 [〇，1 ]的随机数以使得每次迭代输入值该变量具有一定的随机性，(^及(32分别代表经验学习 因子与社会学习因子。4. 根据权利要求1所述的基于SVR-PS0的超声加工特种刀具优化设计方法，其特征在 于:所述步骤四中目标函数为： Y(l，d，m) = l-[0.3〇(l，d，m)+0·3(l_Ui(l，d，m))+0·3U2( 1，d，m) +0.1(l-f(l,d,m))] 式中1为刀具长度，m为材料代号，d为刀具厚度，〇(l，d，m)为最大应力，?ΜΜ,πι)为刀尖 振幅，U2(l，d，m)为刀尖横向位移，f(l，d，m)为固有频率。
【文档编号】G06F17/50GK105956235SQ201610252149
【公开日】2016年9月21日
【申请日】2016年4月20日
【发明人】胡小平, 陈宁池, 董昕頔
【申请人】杭州电子科技大学