一种基于动态加工特征的切削参数分段和变切深优化方法
一种基于动态加工特征的切削参数分段和变切深优化方法
【技术领域】
[0001 ]本发明属于数控加工技术领域，设及对加工工件的切削参数优化方法，为一种基 于动态加工特征的切削参数分段和变切深优化方法。
【背景技术】
[0002] 飞机结构件具有尺寸大、富含薄壁结构、材料去除高、加工精度高等特点，由于W 上特点的存在，飞机结构件数控加工一直是个难题。在飞机结构件的数控加工过程中，常常 存在由于切削参数选择不当而引起工件变形、颤振的问题，严重时造成工件报废。在实际生 产中，为了保证加工质量，一般选用保守的切削参数加工，却降低了生产效率。进一步的，对 于复杂结构特征，加工变形与颤振问题仍然难W避免。
[0003] 由于切削参数在加工中的重要性，近几十年，切削参数优化在传统加工中一直是 一个比较热口的研究。最初，大多数切削参数选择的方法是基于经验;现在，实验方法被广 泛应用于一些特殊的加工状态W获得切削参数，另外有基于人工智能算法进行切削参数优 化，包括基于神经网络的优化方法和基于遗传算法的优化方法，对切削速度、切深、转速多 目标进行优化，并考虑切削力、切削功率、刀具寿命的约束。
[0004] 加工过程中存在切削力和工艺稳定等影响参数优化的因素，W.X.Tang等人的论文 "Prediction of chatter stability in high-speed finishing end milling considering multi-mode dynamics. journal of materials processing technology, 2009,209(5) :2585-2591."中，建立了切削力和工艺稳定的分析模型，并被应用于切削参数 优化过程;Alan等人的论文''Analytical Prediction of Part Dynamics for Machining Stability Analysis . International Journal of Automation Technology(Special Issue on Mod&Sim Cutting Process,2010,4(3):259-267/'中，考虑了加工工艺的改变， 在不同切削深度处，预测工件动态特性，构建稳定区域，并在加工稳定的约束下，选定一个 主轴转速，使切深在一定范围内变化。
[000引但上述基于经验或实验获取参数，一般只适用于简单的工件，并且得到的参数比 较保守，成本高、花费时间长;人工智能算法，有些算法容易形成局部优化或得不到最优解， 而且一般都是针对工件的最终状态，很少考虑工件的中间状态;基于稳定域的切削参数优 化，一般直接在稳定域区间内取点，并没有考虑切削参数的AI算法优化。
[0006] 在数控加工过程中，恒定的进给速度会影响机床运动的平稳性，同时在加工时，不 能充分利用机床的性能，严重影响了加工效率，尤其是在转角加工时，进给速度方向急剧变 化会导致机床的平稳性下降。为了提高加工效率，增强机床的平稳性，需根据机床的运动特 性对转角加工速度进行优化。
【发明内容】
[0007] 本发明要解决的问题是:在飞机结构件的数控加工过程中，对切削参数选择不当 会引起工件变形、颤振的问题，现有技术的解决方案不能满足需求，加工效率低，机床的平 稳性不足，需根据机床的运动特性进行优化。
[0008] 本发明的技术方案为：一种基于动态加工特征的切削参数分段和变切深优化方 法，首先有限元分析工件模型，获取工件的初始频率响应函数矩阵，并基于工件的工艺数据 得到工件中间加工状态，然后基于加工后工件的修改矩阵及获取的工件初始频率响应函数 矩阵，预测工件的实时频率响应函数，实验获取刀具刀尖点频率响应函数，进而得到中间加 工状态的频率响应函数，再计算颤振稳定性叶瓣图，在叶瓣图内选择切深范围，同时考虑加 工的约束条件，建立遗传算法优化模型，得到最优切深、切宽、进给和转速切削参数，切深一 定时，对每一层刀轨，根据机床的运动特性优化转角处进给速度，实现切削参数分段优化。
[0009] 本发明包括W下步骤：
[0010] 步骤1、导入工件模型，并用有限元分析，得到工件的初始频率响应函数；
[0011] 步骤2、提取工件的工艺数据，包括工件的加工工艺参数信息、刀具的几何信息、机 床信息、加工操作顺序W及上一步工序加工余量，所述加工工艺参数信息包括切深、切宽、 转速和进给，由工艺数据得到工件中间加工状态；
[0012] 步骤3、加工后工件结构改变，由几何参数改变量得到修改矩阵，基于修改矩阵和 工件初始频率响应函数，预测加工过程中工件的实时频率响应函数，并由实验获取刀具刀 尖点频率响应函数，进而得到中间加工状态的频率响应函数；
[0013] 步骤4、根据预测的中间加工状态的频率响应函数计算颤振稳定性叶瓣图，根据稳 定性叶瓣图选择稳定切削区域内切削深度范围；
[0014] 步骤5、建立优化目标函数，计算约束条件，约束条件包括切削力约束、机床速度和 进给率约束、切深和切宽约束、刀具刚度约束、机床功率约束W及转角进给速率约束；
[001引步骤6、通过遗传算法，W加工效率为目标，结合约束条件，选择最优适合切深、切 宽、转速和进给；
[0016] 步骤7、在切深一定时，对每一层刀轨，根据机床的运动特性优化转角处进给速度， 在转角处根据机床的运动特性优化进给速度，在直线进给处设置机床允许的最大进给速 度，实现切削参数分段优化。
[0017] 本发明基于动态加工特征的切削参数分段和变切深优化方法，同时考虑动态加工 特征的中间加工状态和切削参数的人工智能算法优化，基于动态加工特征可W对加工过程 中不同状态下的工件进行参数优化，根据工件的不同状态下的动态特性，设置不同的切深， 保证加工质量，提高加工效率;切深一定时，对每一层刀轨，根据刀轨的几何变化特点及机 床的运动特性，分段优化，保证加工效率和机床运动的平稳性。
【附图说明】
[0018] 图1为本发明的基于动态特征的切削参数分段和变切深优化方法流程图。
[0019] 图2为本发明的在颤振稳定性叶瓣图内分区域选择切削深度图。
【具体实施方式】
[0020] 本发明提供了一种基于动态加工特征的切削参数分段和变切深优化方法，首先用 有限元软件分析工件模型，获取工件的初始频率响应函数，基于工件的工艺信息得到工件 中间加工状态，然后基于工件结构改变后的修改矩阵及获取的工件初始频率响应函数，预 测工件的实时频率响应函数，并实验获取刀具刀尖点频率响应函数，计算中间加工状态的 频率响应函数，进而建立稳定性叶瓣图，在叶瓣图内选择切深范围和转速范围，同时考虑其 余约束条件，基于遗传算法建立优化模型，然后根据机床的运动特性优化转角处进给速度。
[0021] 下面结合附图对本发明的技术方案进行详细说明。
[0022] 图1是本发明的基于动态特征的切削参数分段和变切深优化方法流程图，如图所 示，包括W下各步骤：
[0023] 1、导入工件模型，并用有限元软件进行分析，获取工件的初始频率响应函数；
[0024] 2、提取所选工件的工艺信息，刀具的几何信息，机床信息，加工工艺参数信息包括 切深、切宽、转速和进给切削参数，上一步工序加工余量，加工操作顺序。用G_D表示中间状 态的工件几何尺寸，D_F表示最终状态的工件几何尺寸，f (ms)表示由加工顺序决定的工件 尺寸，P_P表示加工工艺参数信息，f(fs)表示一个动态加工特征的所有加工顺序所决定的 工件尺寸，得到中间状态的表达式：
[002引 G_D = D_F+f(ms)+P_P+f(fs) (2.1)
[0026] 3、预测中间加工状态的频率响应函数。基于工件修改后的矩阵及步骤1中获取的 工件结构初始频率响应函数，预测加工过程中工件的动态特性，计算工件实 时频率响应函 数。其计算公式如下：
[0027] 工件的初始频率响应函数矩阵为：
[002引[a] = [[K]-c〇2[M]+i 阳]]-1 (3.1)
[0029] 式中，[M]，[H]，[K]分别是系统的质量矩阵，结构阻尼矩阵及刚度矩阵，ω为激励 频率，i为复数的虚部表示。工件结构的改变可W用修改矩阵[ΔΜ]，[Δ扣，[ΔΚ]来表示，工 件结构改变后的实时频率响应函数矩阵[GwDrkpiwe( ω )]表示为：
[0030] [Gworkpiece(W)] = [[K] + [AK]]-C〇2[[M] + [AM]] + i[W] + [AH]]]-l (3.2)
[0031 ]其中，[ΔΜ]，[ΔΗ]，[ΔΚ]可由实际加工中去除材料的体积及尺寸确定。刀尖点频 率响应函数矩阵可由实验测得，为[GtDDi( ω )]。则中间状态频率响应函数[G( ω )]可W表示 为：
[0032] [G( ω ) ] = [GtD0l( W )+GwDrkpiece( W ) ] (3.3)
[0033] 4、计算颤振稳定性叶瓣图，选择叶瓣图稳定切削区域内转速范围和切削深度范 围：
[0034] 无颤振条件下临界轴向切深：
[0035]
(4.1)
[0036] 式中，N为刀具刀齿数目，Kt为切向力系数，Ar为动态锐削系统特征值Λ的实数部 分，由下式确定：
[0037] (4.2)
[003引 Λ = Λκ+? 八 I [0039] Λ I为特征值的虚部。
[0040]其中：
[0041 ]运里Gxx和Gyy分别表示X和y方向上中间加工状态的频率响应函数，为中间状态频 率响应函数矩阵[G( ω )]中的元素，axxQyyaxyayx为刀具定向动态锐削力系数，ω。为颤振频 率，i为复数的虚部表示。
[0042] K表示特征值虚部和实部的比值：
[0043]
(4,4)
[0044] Ar为特征值Λ的实部，Λι为特征值Λ的虚部，ω。为颤振频率，T为刀齿切削周期。
[0045] 式(4.4)可W用来建立颤振频率和主轴转速之间的关系：
[004引ε为内外调制间的相位差，k为叶瓣数，k = 0，1，2···，η为主轴转速，Ν为刀齿数，Τ为 刀齿切削周期。
[0049] 如图2,叶瓣图下方区域为稳定切削区域，从叶瓣图稳定区域右侧开始依次选取转 速范围及对应的切削深度范围，选择叶瓣图临界点切深处的转速值m与最高点切深处的转 速值m，则选择转速范围的最小值为：
[0050]
(4.巧^
[0051 ]其在叶瓣图线上对应切深为切深范围的最大值apmax，apmax在同一稳定区域叶瓣图 线上对应的转速值为转速范围的最大值nmax，切深范围的最小值apmin取为0。
[0052] 5、计算约束条件，用M(X)表示生产效率，却表示切深，iie表示切宽，fz表示每齿进给 量，N表示齿数，得到的遗传算法优化目标函数公式为：
[0053] M(X)=apXaeXfzXN (5.1)
[0054] 1)用Fn表示法向力，Fmax表示最大预测切削力，根据刀具信息和切削参数得到切削 进给率和切削厚度，再计算出切入、切出角，根据切削力的分析，测量剪切力在切向、径向与 轴向的作用系数分别为KtT、KrT和KaT，与刃口力系数分别为Kw、Kr。和Ka。，Φ为刀具切削角，切 入角为恥，切出角为4e，c为切削进给率，切深为ap，转速η，切削力分解成X向、y向、Ζ向的切削 力，得到切削力的预测公式为：
[0055]
[0056]
[0058] 根据力的合成得到Fmax，得到切削力的约束条件公式为：
[0059] Fn < Fmax 巧.3)
[0060] 2)用Vc表示切削速度，ft表示进给率，VCmax、VCmin分别表示最大最小切削速度， f tmax、f tmin分别表示最大最小进给率，η表示主轴转速，f Z表示每齿进给量，D表示刀具直径， 得到机床的速度和进给率的约束条件为：
[0061] VCmin < Vc < VCmax
[0062] ftmin < ft < ftmax (5.4)
[0063] ft、fz、Vc、N和n的关系式为：
[0064] ft = fzXnXN
[0065] Vc = 3iXDXn/1000 巧.5)
[0066] 3)用ap表示切深，ae表示切宽，apmax、apmin分别表示最大和最小的切深，aemax、aemax 分别表示切宽的最大值和最小值，切深和切宽的约束表达式为：
[0067] apmin < ap < apmax
[0068] aemin < ae < aemax (5.6)
[0069] 4)用Fr表示径向切削力，1表示外伸长度，E表示弹性模量，I表示转动惯量，δ表示 晓度，Smax晓度的最大值，得到刀具曲率的约束公式为：
[0070]
[0071]
[0072] 5)用Pq表示机床实际运转功率，Pe表示机床的额定功率，η表示机床的效率，Ft表示 力，d表示直径，得到机床功率的约束公式为：
[0073]
[0074]
[007引6)ACC表示给定机床的最大曲线加速度，Vcorner表示转角速率，Rc表示转角面的曲 率，得到的刀具进给曲率的约束条件为：
[0076]
[0077] 6、利用遗传算法进行优化，根据得到的切削力、机床的速度和进给率、切深和切 宽、机床功率、等约束条件，采用遗传算法，进行全局优化。其中初始种群设置为1〇〇,4个变 量分别是切深、切宽、转速、进给，交叉概率是0.6，变异概率是0.08，W最大效率为目标，通 过迭代500次得到最终优化结果，再进行工艺决策，得到优化的切削参数，获得稳定切削深 度。
[0078] 7、在切深一定时，对每一层刀轨，根据刀轨的几何变化特点，对切削参数进行分段 优化。在加工中，转角处进给速度的最大值可用公式来表示：
[0079] (7.1)
[0080] 式中，Vimax表示第i个刀位点处速度的最大值;Vmax表示机床允许的最大进给;aipmax 表示在第i个刀位点处法向最大加速度值，可W通过刀轨在第i个刀位点处法向同Χ、Υ轴的 夹角W及机床Χ、γ方向的最大加速度值来确定;化表示在第i个刀位点处刀轨的曲率半径， 根据刀轨来确定。
[0081] 本发明基于动态加工特征构建特征的中间加工状态，考虑特征的几何形状和中间 加工状态性质随切削加工过程的变化，根据特征中间加工状态的刚度，采用变切深的切削 参数优化策略;在变切深切削参数优化的基础上，对每一层刀轨考虑机床的动力学特性，根 据特征的几何形状和刀轨的几何变化特点，对切削参数进行分段优化，实现提高加工效率 的目的。
【主权项】
1. 一种基于动态加工特征的切削参数分段和变切深优化方法，其特征是首先有限元分 析工件模型，获取工件的初始频率响应函数矩阵，并基于工件的工艺数据得到工件中间加 工状态，然后基于加工后工件的修改矩阵及获取的工件初始频率响应函数矩阵，预测工件 的实时频率响应函数，实验获取刀具刀尖点频率响应函数，进而得到中间加工状态的频率 响应函数，再计算颤振稳定性叶瓣图，在叶瓣图内选择切深范围，同时考虑加工的约束条 件，建立遗传算法优化模型，得到最优切深、切宽、进给和转速切削参数，切深一定时，对每 一层刀轨，根据机床的运动特性优化转角处进给速度，实现切削参数分段优化。2. 根据权利要求1所述的一种基于动态加工特征的切削参数分段和变切深优化方法， 其特征是包括以下步骤： 步骤1、导入工件模型，并用有限元 分析，得到工件的初始频率响应函数； 步骤2、提取工件的工艺数据，包括工件的加工工艺参数信息、刀具的几何信息、机床信 息、加工操作顺序以及上一步工序加工余量，所述加工工艺参数信息包括切深、切宽、转速 和进给，由工艺数据得到工件中间加工状态； 步骤3、加工后工件结构改变，由几何参数改变量得到修改矩阵，基于修改矩阵和工件 初始频率响应函数，预测加工过程中工件的实时频率响应函数，并由实验获取刀具刀尖点 频率响应函数，进而得到中间加工状态的频率响应函数； 步骤4、根据预测的中间加工状态的频率响应函数计算颤振稳定性叶瓣图，根据稳定性 叶瓣图选择稳定切削区域内切削深度范围； 步骤5、建立优化目标函数，计算约束条件，约束条件包括切削力约束、机床速度和进给 率约束、切深和切宽约束、刀具刚度约束、机床功率约束以及转角进给速率约束； 步骤6、通过遗传算法，以加工效率为目标，结合约束条件，选择最优适合切深、切宽、主 轴转速和进给； 步骤7、在切深一定时，对每一层刀轨，根据机床的运动特性优化转角处进给速度，在转 角处根据机床的运动特性优化进给速度，在直线进给处设置机床允许的最大进给速度，实 现切削参数分段优化。3. 根据权利要求1或2所述的一种基于动态加工特征的切削参数分段和变切深优化方 法，其特征是中间加工状态的频率响应函数的获取具体为： 工件的初始频率响应函数矩阵为： [α] = [[Κ]-ω2[Μ]+?[Η]]-1 (3.1) 式中，[Μ]，[Η]，[Κ]分别是工件的质量矩阵，结构阻尼矩阵及刚度矩阵，ω为激励频率， i为复数的虚部表示，工件结构的改变用修改矩阵[ΛΜ]，[ ΛΗ]，[ ΛΚ]来表示，由此得到工 件结构改变，即加工过程中的实时频率响应函数[Gwcirkpi_( ω )]，表示为： [Gw〇rkpiece(〇 )] = [[[Κ]+[ΔΚ]]-ω2[[Μ]+[ΔΜ]]+?[[Η]+[ΔΗ]]]_1 (3.2) 其中，[ΛΜ]，[ΛΗ]，[ΛΚ]由加工中去除材料的体积及尺寸确定； 由实验测得刀尖点频率响应函数[Gt?i( ω )]，则中间加工状态的频率响应函数[G( ω )] 表不为： [G(W)] - [Gtool ( 〇 ) +Gworkpiece (〇)] (3.3)。4. 根据权利要求1或2所述的一种基于动态加工特征的切削参数分段和变切深优化方 法，其特征是计算颤振稳定性叶瓣图，选择叶瓣图稳定切削区域内切削深度范围具体为： 无颤振条件下临界轴向切深aiim为：式中，N为刀具刀齿数目，Kt为切向力系数，Λ R为动态铣削系统特征值Λ的实数部分，由 下式确定：Λ I为特征值的虚部； ?0 - Gxx ( I W c)Gyy ( I W c) (axxayy_a XyOyx) (4.3) ?1 一 ClxxGxx ( i W c ) +OyyGyy ( i W c ) Gxx和Gyy分别表示x和y方向上中间加工状态的频率响应函数，为中间状态频率响应函数 矩阵[G( ω )]中的元素，€^、€^、€^和€^为刀具定向动态铣削力系数，(0。为颤振频率，：[为复 数的虚部表示； K表示特征值虚部和实部的比值：T为刀齿切削周期； 由式(4.4)建立颤振频率和主轴转速之间的关系：ε为内外调制间的相位差，k为叶瓣数，k = 0,1,2···，η为主轴转速； 叶瓣图下方区域为稳定切削区域，从叶瓣图稳定区域右侧开始依次选取转速范围及对 应的切削深度范围，选择叶瓣图临界点切深处的转速值m与最高点切深处的转速值η2，则选 择转速范围的最小值为：其在叶瓣图线上对应切深为切深范围的最大值apmax，apmax在同一稳定区域叶瓣图线上 对应的转速值为转速范围的最大值nmax，切深范围的最小值apmin取为0。5.根据权利要求1或2所述的一种基于动态加工特征的切削参数分段和变切深优化方 法，其特征是根据优化变量得到优化目标函数： 用M(X)表示生产效率，aP表示切深，ae表示切宽，fz表示每齿进给量，N表示刀具刀齿数 目，得到的遗传算法优化目标函数公式为： M(X)=apXaeXfzXN (5.1) 1) 用Fn表示法向力，Fmax表示最大预测切削力，根据刀具信息和切削参数得到切削进给 率和切削厚度，再计算出切入、切出角，根据切削力的分析，测量剪切力在切向、径向与轴向 的作用系数分别为Ku、UPK aT，与刃口力系数分别为Kt。、Kr。和Κμ,Φ为刀具切削角，切入角 为也，切出角为也，c为切削进给率，切深Sa p，主轴转速η，切削力分解成X向、y向、z向的切削 力，得到切削力的预测公式为：I 根据力的合成得到Fmax，切削力的约束条件公式为： Fn < Fmax (5.3) 2) 用Vc表示切削速度，ft表示进给率，Vcmax、Vcmin分别表示最大最小切削速度，ft max、 ftmin分别表示最大最小进给率，η表示主轴转速，fZ表示每齿进给量，D表示刀具直径，得到 机床的速度和进给率的约束条件为： VCmin ^ Vc ^ VCmax ftmin < ft < ftmax (5.4) ft、fz、Vc、N和n的关系式为： ft = fz XnXN Vc = JiXDXn/1000 (5.5) 3) 用ap表不切深，ae表不切宽，apmax、apmin分别表不最大和最小的切深，ae max、aemax分别 表示切宽的最大值和最小值，切深和切宽的约束表达式为： ciPmin ^ cip ^ ciPmax ￡1βιη?η - ￡1β ^ clGmax (5.6) 4) 用Fr表示径向切削力，1表示外伸长度，E表示弹性模量，I表示转动惯量，δ表示烧度， Smax挠度的最大值，得到刀具曲率的约束公式为：5) 用Pq表示机床实际运转功率，Pe3表示机床的额定功率，η表示机床的效率，Ft表示力，d 表示直径，得到机床功率的约束公式为：6 )Acc表示给定机床的最大曲线加速度，Vc_^表示转角速率，R。表示转角面的曲率，得 到的刀具进给曲率的约束条件为：6. 根据权利要求1或2所述的一种基于动态加工特征的切削参数分段和变切深优化方 法，其特征是利用遗传算法进行优化时，根据约束条件，采用遗传算法进行全局优化，其中 初始种群设置为100，4个变量分别是切深、切宽、转速和进给，交叉概率为0.6，变异概率为 0.08，以最大效率为目标，迭代得到最终优化结果，再进行工艺决策，得到优化的切削参数， 获得稳定切削深度。7. 根据权利要求1或2所述的一种基于动态加工特征的切削参数分段和变切深优化方 法，其特征是切削参数分段优化具体为： 在加工中，转角处进给速度的最大值用公式(7.1)来表示：式中，Vimax表示第i个刀位点处速度的最大值;Vmax表示机床允许的最大进给;a iPmax表示 在第i个刀位点处法向最大加速度值，aiPmax通过刀轨在第i个刀位点处法向同X、Y轴的夹角 以及机床Χ、Υ方向的最大加速度值来确定;Pi表示在第i个刀位点处刀轨的曲率半径，根据 刀轨来确定。
【专利摘要】一种基于动态加工特征的切削参数分段和变切深优化方法，首先用有限元软件分析工件模型，获取工件的初始频率响应函数，基于工件的工艺信息得到工件中间加工状态，然后基于工件结构改变后的修改矩阵及获取的工件初始频率响应函数，预测工件的实时频率响应函数，实验获取刀具刀尖点频率响应函数，计算中间加工状态的频率响应函数，进而建立稳定性叶瓣图，在叶瓣图内选择切深范围和转速范围，同时考虑其余约束条件，基于遗传算法建立优化模型，然后根据机床的运动特性优化转角处进给速度。本发明基于动态加工特征优化切削参数，保证加工质量并提高了加工效率，并且由于分段优化过程考虑了机床动态特性，保证了机床运动的平稳性。
【IPC分类】G06F17/50
【公开号】CN105488282
【申请号】CN201510864151
【发明人】李迎光, 刘长青
【申请人】南京航空航天大学
【公开日】2016年4月13日
【申请日】2015年11月30日
【技术领域】
[0001 ]本发明属于数控加工技术领域，设及对加工工件的切削参数优化方法，为一种基 于动态加工特征的切削参数分段和变切深优化方法。
【背景技术】
[0002] 飞机结构件具有尺寸大、富含薄壁结构、材料去除高、加工精度高等特点，由于W 上特点的存在，飞机结构件数控加工一直是个难题。在飞机结构件的数控加工过程中，常常 存在由于切削参数选择不当而引起工件变形、颤振的问题，严重时造成工件报废。在实际生 产中，为了保证加工质量，一般选用保守的切削参数加工，却降低了生产效率。进一步的，对 于复杂结构特征，加工变形与颤振问题仍然难W避免。
[0003] 由于切削参数在加工中的重要性，近几十年，切削参数优化在传统加工中一直是 一个比较热口的研究。最初，大多数切削参数选择的方法是基于经验;现在，实验方法被广 泛应用于一些特殊的加工状态W获得切削参数，另外有基于人工智能算法进行切削参数优 化，包括基于神经网络的优化方法和基于遗传算法的优化方法，对切削速度、切深、转速多 目标进行优化，并考虑切削力、切削功率、刀具寿命的约束。
[0004] 加工过程中存在切削力和工艺稳定等影响参数优化的因素，W.X.Tang等人的论文 "Prediction of chatter stability in high-speed finishing end milling considering multi-mode dynamics. journal of materials processing technology, 2009,209(5) :2585-2591."中，建立了切削力和工艺稳定的分析模型，并被应用于切削参数 优化过程;Alan等人的论文''Analytical Prediction of Part Dynamics for Machining Stability Analysis . International Journal of Automation Technology(Special Issue on Mod&Sim Cutting Process,2010,4(3):259-267/'中，考虑了加工工艺的改变， 在不同切削深度处，预测工件动态特性，构建稳定区域，并在加工稳定的约束下，选定一个 主轴转速，使切深在一定范围内变化。
[000引但上述基于经验或实验获取参数，一般只适用于简单的工件，并且得到的参数比 较保守，成本高、花费时间长;人工智能算法，有些算法容易形成局部优化或得不到最优解， 而且一般都是针对工件的最终状态，很少考虑工件的中间状态;基于稳定域的切削参数优 化，一般直接在稳定域区间内取点，并没有考虑切削参数的AI算法优化。
[0006] 在数控加工过程中，恒定的进给速度会影响机床运动的平稳性，同时在加工时，不 能充分利用机床的性能，严重影响了加工效率，尤其是在转角加工时，进给速度方向急剧变 化会导致机床的平稳性下降。为了提高加工效率，增强机床的平稳性，需根据机床的运动特 性对转角加工速度进行优化。
【发明内容】
[0007] 本发明要解决的问题是:在飞机结构件的数控加工过程中，对切削参数选择不当 会引起工件变形、颤振的问题，现有技术的解决方案不能满足需求，加工效率低，机床的平 稳性不足，需根据机床的运动特性进行优化。
[0008] 本发明的技术方案为：一种基于动态加工特征的切削参数分段和变切深优化方 法，首先有限元分析工件模型，获取工件的初始频率响应函数矩阵，并基于工件的工艺数据 得到工件中间加工状态，然后基于加工后工件的修改矩阵及获取的工件初始频率响应函数 矩阵，预测工件的实时频率响应函数，实验获取刀具刀尖点频率响应函数，进而得到中间加 工状态的频率响应函数，再计算颤振稳定性叶瓣图，在叶瓣图内选择切深范围，同时考虑加 工的约束条件，建立遗传算法优化模型，得到最优切深、切宽、进给和转速切削参数，切深一 定时，对每一层刀轨，根据机床的运动特性优化转角处进给速度，实现切削参数分段优化。
[0009] 本发明包括W下步骤：
[0010] 步骤1、导入工件模型，并用有限元分析，得到工件的初始频率响应函数；
[0011] 步骤2、提取工件的工艺数据，包括工件的加工工艺参数信息、刀具的几何信息、机 床信息、加工操作顺序W及上一步工序加工余量，所述加工工艺参数信息包括切深、切宽、 转速和进给，由工艺数据得到工件中间加工状态；
[0012] 步骤3、加工后工件结构改变，由几何参数改变量得到修改矩阵，基于修改矩阵和 工件初始频率响应函数，预测加工过程中工件的实时频率响应函数，并由实验获取刀具刀 尖点频率响应函数，进而得到中间加工状态的频率响应函数；
[0013] 步骤4、根据预测的中间加工状态的频率响应函数计算颤振稳定性叶瓣图，根据稳 定性叶瓣图选择稳定切削区域内切削深度范围；
[0014] 步骤5、建立优化目标函数，计算约束条件，约束条件包括切削力约束、机床速度和 进给率约束、切深和切宽约束、刀具刚度约束、机床功率约束W及转角进给速率约束；
[001引步骤6、通过遗传算法，W加工效率为目标，结合约束条件，选择最优适合切深、切 宽、转速和进给；
[0016] 步骤7、在切深一定时，对每一层刀轨，根据机床的运动特性优化转角处进给速度， 在转角处根据机床的运动特性优化进给速度，在直线进给处设置机床允许的最大进给速 度，实现切削参数分段优化。
[0017] 本发明基于动态加工特征的切削参数分段和变切深优化方法，同时考虑动态加工 特征的中间加工状态和切削参数的人工智能算法优化，基于动态加工特征可W对加工过程 中不同状态下的工件进行参数优化，根据工件的不同状态下的动态特性，设置不同的切深， 保证加工质量，提高加工效率;切深一定时，对每一层刀轨，根据刀轨的几何变化特点及机 床的运动特性，分段优化，保证加工效率和机床运动的平稳性。
【附图说明】
[0018] 图1为本发明的基于动态特征的切削参数分段和变切深优化方法流程图。
[0019] 图2为本发明的在颤振稳定性叶瓣图内分区域选择切削深度图。
【具体实施方式】
[0020] 本发明提供了一种基于动态加工特征的切削参数分段和变切深优化方法，首先用 有限元软件分析工件模型，获取工件的初始频率响应函数，基于工件的工艺信息得到工件 中间加工状态，然后基于工件结构改变后的修改矩阵及获取的工件初始频率响应函数，预 测工件的实时频率响应函数，并实验获取刀具刀尖点频率响应函数，计算中间加工状态的 频率响应函数，进而建立稳定性叶瓣图，在叶瓣图内选择切深范围和转速范围，同时考虑其 余约束条件，基于遗传算法建立优化模型，然后根据机床的运动特性优化转角处进给速度。
[0021] 下面结合附图对本发明的技术方案进行详细说明。
[0022] 图1是本发明的基于动态特征的切削参数分段和变切深优化方法流程图，如图所 示，包括W下各步骤：
[0023] 1、导入工件模型，并用有限元软件进行分析，获取工件的初始频率响应函数；
[0024] 2、提取所选工件的工艺信息，刀具的几何信息，机床信息，加工工艺参数信息包括 切深、切宽、转速和进给切削参数，上一步工序加工余量，加工操作顺序。用G_D表示中间状 态的工件几何尺寸，D_F表示最终状态的工件几何尺寸，f (ms)表示由加工顺序决定的工件 尺寸，P_P表示加工工艺参数信息，f(fs)表示一个动态加工特征的所有加工顺序所决定的 工件尺寸，得到中间状态的表达式：
[002引 G_D = D_F+f(ms)+P_P+f(fs) (2.1)
[0026] 3、预测中间加工状态的频率响应函数。基于工件修改后的矩阵及步骤1中获取的 工件结构初始频率响应函数，预测加工过程中工件的动态特性，计算工件实 时频率响应函 数。其计算公式如下：
[0027] 工件的初始频率响应函数矩阵为：
[002引[a] = [[K]-c〇2[M]+i 阳]]-1 (3.1)
[0029] 式中，[M]，[H]，[K]分别是系统的质量矩阵，结构阻尼矩阵及刚度矩阵，ω为激励 频率，i为复数的虚部表示。工件结构的改变可W用修改矩阵[ΔΜ]，[Δ扣，[ΔΚ]来表示，工 件结构改变后的实时频率响应函数矩阵[GwDrkpiwe( ω )]表示为：
[0030] [Gworkpiece(W)] = [[K] + [AK]]-C〇2[[M] + [AM]] + i[W] + [AH]]]-l (3.2)
[0031 ]其中，[ΔΜ]，[ΔΗ]，[ΔΚ]可由实际加工中去除材料的体积及尺寸确定。刀尖点频 率响应函数矩阵可由实验测得，为[GtDDi( ω )]。则中间状态频率响应函数[G( ω )]可W表示 为：
[0032] [G( ω ) ] = [GtD0l( W )+GwDrkpiece( W ) ] (3.3)
[0033] 4、计算颤振稳定性叶瓣图，选择叶瓣图稳定切削区域内转速范围和切削深度范 围：
[0034] 无颤振条件下临界轴向切深：
[0035]
(4.1)
[0036] 式中，N为刀具刀齿数目，Kt为切向力系数，Ar为动态锐削系统特征值Λ的实数部 分，由下式确定：
[0037] (4.2)
[003引 Λ = Λκ+? 八 I [0039] Λ I为特征值的虚部。
[0040]其中：
[0041 ]运里Gxx和Gyy分别表示X和y方向上中间加工状态的频率响应函数，为中间状态频 率响应函数矩阵[G( ω )]中的元素，axxQyyaxyayx为刀具定向动态锐削力系数，ω。为颤振频 率，i为复数的虚部表示。
[0042] K表示特征值虚部和实部的比值：
[0043]
(4,4)
[0044] Ar为特征值Λ的实部，Λι为特征值Λ的虚部，ω。为颤振频率，T为刀齿切削周期。
[0045] 式(4.4)可W用来建立颤振频率和主轴转速之间的关系：
[004引ε为内外调制间的相位差，k为叶瓣数，k = 0，1，2···，η为主轴转速，Ν为刀齿数，Τ为 刀齿切削周期。
[0049] 如图2,叶瓣图下方区域为稳定切削区域，从叶瓣图稳定区域右侧开始依次选取转 速范围及对应的切削深度范围，选择叶瓣图临界点切深处的转速值m与最高点切深处的转 速值m，则选择转速范围的最小值为：
[0050]
(4.巧^
[0051 ]其在叶瓣图线上对应切深为切深范围的最大值apmax，apmax在同一稳定区域叶瓣图 线上对应的转速值为转速范围的最大值nmax，切深范围的最小值apmin取为0。
[0052] 5、计算约束条件，用M(X)表示生产效率，却表示切深，iie表示切宽，fz表示每齿进给 量，N表示齿数，得到的遗传算法优化目标函数公式为：
[0053] M(X)=apXaeXfzXN (5.1)
[0054] 1)用Fn表示法向力，Fmax表示最大预测切削力，根据刀具信息和切削参数得到切削 进给率和切削厚度，再计算出切入、切出角，根据切削力的分析，测量剪切力在切向、径向与 轴向的作用系数分别为KtT、KrT和KaT，与刃口力系数分别为Kw、Kr。和Ka。，Φ为刀具切削角，切 入角为恥，切出角为4e，c为切削进给率，切深为ap，转速η，切削力分解成X向、y向、Ζ向的切削 力，得到切削力的预测公式为：
[0055]
[0056]
[0058] 根据力的合成得到Fmax，得到切削力的约束条件公式为：
[0059] Fn < Fmax 巧.3)
[0060] 2)用Vc表示切削速度，ft表示进给率，VCmax、VCmin分别表示最大最小切削速度， f tmax、f tmin分别表示最大最小进给率，η表示主轴转速，f Z表示每齿进给量，D表示刀具直径， 得到机床的速度和进给率的约束条件为：
[0061] VCmin < Vc < VCmax
[0062] ftmin < ft < ftmax (5.4)
[0063] ft、fz、Vc、N和n的关系式为：
[0064] ft = fzXnXN
[0065] Vc = 3iXDXn/1000 巧.5)
[0066] 3)用ap表示切深，ae表示切宽，apmax、apmin分别表示最大和最小的切深，aemax、aemax 分别表示切宽的最大值和最小值，切深和切宽的约束表达式为：
[0067] apmin < ap < apmax
[0068] aemin < ae < aemax (5.6)
[0069] 4)用Fr表示径向切削力，1表示外伸长度，E表示弹性模量，I表示转动惯量，δ表示 晓度，Smax晓度的最大值，得到刀具曲率的约束公式为：
[0070]
[0071]
[0072] 5)用Pq表示机床实际运转功率，Pe表示机床的额定功率，η表示机床的效率，Ft表示 力，d表示直径，得到机床功率的约束公式为：
[0073]
[0074]
[007引6)ACC表示给定机床的最大曲线加速度，Vcorner表示转角速率，Rc表示转角面的曲 率，得到的刀具进给曲率的约束条件为：
[0076]
[0077] 6、利用遗传算法进行优化，根据得到的切削力、机床的速度和进给率、切深和切 宽、机床功率、等约束条件，采用遗传算法，进行全局优化。其中初始种群设置为1〇〇,4个变 量分别是切深、切宽、转速、进给，交叉概率是0.6，变异概率是0.08，W最大效率为目标，通 过迭代500次得到最终优化结果，再进行工艺决策，得到优化的切削参数，获得稳定切削深 度。
[0078] 7、在切深一定时，对每一层刀轨，根据刀轨的几何变化特点，对切削参数进行分段 优化。在加工中，转角处进给速度的最大值可用公式来表示：
[0079] (7.1)
[0080] 式中，Vimax表示第i个刀位点处速度的最大值;Vmax表示机床允许的最大进给;aipmax 表示在第i个刀位点处法向最大加速度值，可W通过刀轨在第i个刀位点处法向同Χ、Υ轴的 夹角W及机床Χ、γ方向的最大加速度值来确定;化表示在第i个刀位点处刀轨的曲率半径， 根据刀轨来确定。
[0081] 本发明基于动态加工特征构建特征的中间加工状态，考虑特征的几何形状和中间 加工状态性质随切削加工过程的变化，根据特征中间加工状态的刚度，采用变切深的切削 参数优化策略;在变切深切削参数优化的基础上，对每一层刀轨考虑机床的动力学特性，根 据特征的几何形状和刀轨的几何变化特点，对切削参数进行分段优化，实现提高加工效率 的目的。
【主权项】
1. 一种基于动态加工特征的切削参数分段和变切深优化方法，其特征是首先有限元分 析工件模型，获取工件的初始频率响应函数矩阵，并基于工件的工艺数据得到工件中间加 工状态，然后基于加工后工件的修改矩阵及获取的工件初始频率响应函数矩阵，预测工件 的实时频率响应函数，实验获取刀具刀尖点频率响应函数，进而得到中间加工状态的频率 响应函数，再计算颤振稳定性叶瓣图，在叶瓣图内选择切深范围，同时考虑加工的约束条 件，建立遗传算法优化模型，得到最优切深、切宽、进给和转速切削参数，切深一定时，对每 一层刀轨，根据机床的运动特性优化转角处进给速度，实现切削参数分段优化。2. 根据权利要求1所述的一种基于动态加工特征的切削参数分段和变切深优化方法， 其特征是包括以下步骤： 步骤1、导入工件模型，并用有限元 分析，得到工件的初始频率响应函数； 步骤2、提取工件的工艺数据，包括工件的加工工艺参数信息、刀具的几何信息、机床信 息、加工操作顺序以及上一步工序加工余量，所述加工工艺参数信息包括切深、切宽、转速 和进给，由工艺数据得到工件中间加工状态； 步骤3、加工后工件结构改变，由几何参数改变量得到修改矩阵，基于修改矩阵和工件 初始频率响应函数，预测加工过程中工件的实时频率响应函数，并由实验获取刀具刀尖点 频率响应函数，进而得到中间加工状态的频率响应函数； 步骤4、根据预测的中间加工状态的频率响应函数计算颤振稳定性叶瓣图，根据稳定性 叶瓣图选择稳定切削区域内切削深度范围； 步骤5、建立优化目标函数，计算约束条件，约束条件包括切削力约束、机床速度和进给 率约束、切深和切宽约束、刀具刚度约束、机床功率约束以及转角进给速率约束； 步骤6、通过遗传算法，以加工效率为目标，结合约束条件，选择最优适合切深、切宽、主 轴转速和进给； 步骤7、在切深一定时，对每一层刀轨，根据机床的运动特性优化转角处进给速度，在转 角处根据机床的运动特性优化进给速度，在直线进给处设置机床允许的最大进给速度，实 现切削参数分段优化。3. 根据权利要求1或2所述的一种基于动态加工特征的切削参数分段和变切深优化方 法，其特征是中间加工状态的频率响应函数的获取具体为： 工件的初始频率响应函数矩阵为： [α] = [[Κ]-ω2[Μ]+?[Η]]-1 (3.1) 式中，[Μ]，[Η]，[Κ]分别是工件的质量矩阵，结构阻尼矩阵及刚度矩阵，ω为激励频率， i为复数的虚部表示，工件结构的改变用修改矩阵[ΛΜ]，[ ΛΗ]，[ ΛΚ]来表示，由此得到工 件结构改变，即加工过程中的实时频率响应函数[Gwcirkpi_( ω )]，表示为： [Gw〇rkpiece(〇 )] = [[[Κ]+[ΔΚ]]-ω2[[Μ]+[ΔΜ]]+?[[Η]+[ΔΗ]]]_1 (3.2) 其中，[ΛΜ]，[ΛΗ]，[ΛΚ]由加工中去除材料的体积及尺寸确定； 由实验测得刀尖点频率响应函数[Gt?i( ω )]，则中间加工状态的频率响应函数[G( ω )] 表不为： [G(W)] - [Gtool ( 〇 ) +Gworkpiece (〇)] (3.3)。4. 根据权利要求1或2所述的一种基于动态加工特征的切削参数分段和变切深优化方 法，其特征是计算颤振稳定性叶瓣图，选择叶瓣图稳定切削区域内切削深度范围具体为： 无颤振条件下临界轴向切深aiim为：式中，N为刀具刀齿数目，Kt为切向力系数，Λ R为动态铣削系统特征值Λ的实数部分，由 下式确定：Λ I为特征值的虚部； ?0 - Gxx ( I W c)Gyy ( I W c) (axxayy_a XyOyx) (4.3) ?1 一 ClxxGxx ( i W c ) +OyyGyy ( i W c ) Gxx和Gyy分别表示x和y方向上中间加工状态的频率响应函数，为中间状态频率响应函数 矩阵[G( ω )]中的元素，€^、€^、€^和€^为刀具定向动态铣削力系数，(0。为颤振频率，：[为复 数的虚部表示； K表示特征值虚部和实部的比值：T为刀齿切削周期； 由式(4.4)建立颤振频率和主轴转速之间的关系：ε为内外调制间的相位差，k为叶瓣数，k = 0,1,2···，η为主轴转速； 叶瓣图下方区域为稳定切削区域，从叶瓣图稳定区域右侧开始依次选取转速范围及对 应的切削深度范围，选择叶瓣图临界点切深处的转速值m与最高点切深处的转速值η2，则选 择转速范围的最小值为：其在叶瓣图线上对应切深为切深范围的最大值apmax，apmax在同一稳定区域叶瓣图线上 对应的转速值为转速范围的最大值nmax，切深范围的最小值apmin取为0。5.根据权利要求1或2所述的一种基于动态加工特征的切削参数分段和变切深优化方 法，其特征是根据优化变量得到优化目标函数： 用M(X)表示生产效率，aP表示切深，ae表示切宽，fz表示每齿进给量，N表示刀具刀齿数 目，得到的遗传算法优化目标函数公式为： M(X)=apXaeXfzXN (5.1) 1) 用Fn表示法向力，Fmax表示最大预测切削力，根据刀具信息和切削参数得到切削进给 率和切削厚度，再计算出切入、切出角，根据切削力的分析，测量剪切力在切向、径向与轴向 的作用系数分别为Ku、UPK aT，与刃口力系数分别为Kt。、Kr。和Κμ,Φ为刀具切削角，切入角 为也，切出角为也，c为切削进给率，切深Sa p，主轴转速η，切削力分解成X向、y向、z向的切削 力，得到切削力的预测公式为：I 根据力的合成得到Fmax，切削力的约束条件公式为： Fn < Fmax (5.3) 2) 用Vc表示切削速度，ft表示进给率，Vcmax、Vcmin分别表示最大最小切削速度，ft max、 ftmin分别表示最大最小进给率，η表示主轴转速，fZ表示每齿进给量，D表示刀具直径，得到 机床的速度和进给率的约束条件为： VCmin ^ Vc ^ VCmax ftmin < ft < ftmax (5.4) ft、fz、Vc、N和n的关系式为： ft = fz XnXN Vc = JiXDXn/1000 (5.5) 3) 用ap表不切深，ae表不切宽，apmax、apmin分别表不最大和最小的切深，ae max、aemax分别 表示切宽的最大值和最小值，切深和切宽的约束表达式为： ciPmin ^ cip ^ ciPmax ￡1βιη?η - ￡1β ^ clGmax (5.6) 4) 用Fr表示径向切削力，1表示外伸长度，E表示弹性模量，I表示转动惯量，δ表示烧度， Smax挠度的最大值，得到刀具曲率的约束公式为：5) 用Pq表示机床实际运转功率，Pe3表示机床的额定功率，η表示机床的效率，Ft表示力，d 表示直径，得到机床功率的约束公式为：6 )Acc表示给定机床的最大曲线加速度，Vc_^表示转角速率，R。表示转角面的曲率，得 到的刀具进给曲率的约束条件为：6. 根据权利要求1或2所述的一种基于动态加工特征的切削参数分段和变切深优化方 法，其特征是利用遗传算法进行优化时，根据约束条件，采用遗传算法进行全局优化，其中 初始种群设置为100，4个变量分别是切深、切宽、转速和进给，交叉概率为0.6，变异概率为 0.08，以最大效率为目标，迭代得到最终优化结果，再进行工艺决策，得到优化的切削参数， 获得稳定切削深度。7. 根据权利要求1或2所述的一种基于动态加工特征的切削参数分段和变切深优化方 法，其特征是切削参数分段优化具体为： 在加工中，转角处进给速度的最大值用公式(7.1)来表示：式中，Vimax表示第i个刀位点处速度的最大值;Vmax表示机床允许的最大进给;a iPmax表示 在第i个刀位点处法向最大加速度值，aiPmax通过刀轨在第i个刀位点处法向同X、Y轴的夹角 以及机床Χ、Υ方向的最大加速度值来确定;Pi表示在第i个刀位点处刀轨的曲率半径，根据 刀轨来确定。
【专利摘要】一种基于动态加工特征的切削参数分段和变切深优化方法，首先用有限元软件分析工件模型，获取工件的初始频率响应函数，基于工件的工艺信息得到工件中间加工状态，然后基于工件结构改变后的修改矩阵及获取的工件初始频率响应函数，预测工件的实时频率响应函数，实验获取刀具刀尖点频率响应函数，计算中间加工状态的频率响应函数，进而建立稳定性叶瓣图，在叶瓣图内选择切深范围和转速范围，同时考虑其余约束条件，基于遗传算法建立优化模型，然后根据机床的运动特性优化转角处进给速度。本发明基于动态加工特征优化切削参数，保证加工质量并提高了加工效率，并且由于分段优化过程考虑了机床动态特性，保证了机床运动的平稳性。
【IPC分类】G06F17/50
【公开号】CN105488282
【申请号】CN201510864151
【发明人】李迎光, 刘长青
【申请人】南京航空航天大学
【公开日】2016年4月13日
【申请日】2015年11月30日
最新回复(0)