一种基于Bull公式的铣削稳定性预测方法技术

本发明专利技术属于数控技术领域，尤其涉及一种基于Bull公式的铣削稳定性预测方法，主要用于多轴加工中铣削稳定性的预测。该方法和其它预测铣削稳定性的一阶、二阶半离散或全离散方法相比，不仅能更准确的预测铣削稳定性，还具有更高的预测效率。

A prediction method of milling stability based on Bull formula

The invention belongs to the field of numerical control technology, in particular to a method for predicting the milling stability based on the Bull formula, which is mainly used for the prediction of the milling stability in multi axis machining. Compared with other first-order or two order semi discrete or fully discrete methods for predicting milling stability, this method not only can predict milling stability more accurately, but also has higher prediction efficiency.

【技术实现步骤摘要】
一种基于Bull公式的铣削稳定性预测方法
：本专利技术属于数控
，尤其涉及一种基于Bull公式的铣削稳定性预测方法，主要用于多轴加工中铣削稳定性的预测。
技术介绍
：铣削是一种加工方法，能够从金属试样中切除大量的金属材料。该方法能够保证作业过程的精确性、高效性。铣削颤振常常发生在铣削加工过程中,是一种自激振动现象。铣削颤振的出现将会加剧刀具磨损，降低加工中心的使用寿命，恶化工件表面加工质量，严重制约加工中心的使用效率。因此，应及早预测铣削过程中铣刀的稳定性，选取合理的加工参数，避免发生铣削颤振现象。按照颤振机理的不同，加工过程中的颤振分为再生颤振、模态耦合颤振、摩擦颤振和力-热颤振。就铣削加工过程而言，再生颤振是引起加工过程失稳的主要因素。考虑再生效应的动态铣削过程一般可用包含周期系数矩阵的时滞微分方程组来描述。描述铣削过程再生效应动力学的时滞微分方程组具有无穷维状态空间，分析其运动稳定性具有相当的难度。到目前为止，国内外学者针对这一类时滞微分方程组，已提出许多近似方法以预报由临界加工参数构成的铣削过程稳定性边界，但目前铣削稳定性预测方法仍然存在预测准确度不高、效率低和通用性差的问题。
技术实现思路
：为了解决铣削稳定性预测方法存在的问题，本专利技术提出了一种基于Bull公式的铣削稳定性预测方法。和其它预测铣削稳定性的一阶、二阶半离散或全离散方法相比，新方法不仅能更准确的预测铣削稳定性，还具有更高的预测效率。一种基于Bull公式的铣削稳定性预测方法，包括以下步骤：1)将包含再生效应的铣削过程的动力学方程组通过变换转化为空间状态形式在结构动力学框架下，包括再生效应的铣削过程的动力学模型可由以下时滞微分方程组描述其中，M、C和K分别是刀具的模态质量、模态阻尼和模态刚度矩阵；q(t)为刀具模态坐标，且振型系数在刀尖点处归一，Kc(t)为周期系数矩阵，T为时滞量且等于刀齿切削周期，T＝60/(NΩ),且N为刀具齿数，Ω为刀具主轴转速，单位为rpm。令和x(t)＝[q(t)p(t)]T。通过变换，式(1)可以转换为如下的空间状态形式：其中，A表示常数矩阵；H(t)表示周期为T的考虑再生效应的系数矩阵，且H(t)＝H(t+T)。其中令G(t)＝A-H(t),F(t)＝G(t)x(t)+H(t)x(t-T),则(2)变为2)将刀齿切削周期T离散为等距的m个小区间，则任意一个时间小区间为[ti,ti+1]，其中ti＝t1+(i-1)τ(i＝1,2,…m+1)，由(3)式和牛顿-莱布尼兹公式得：利用数值积分中的Bull公式可得:将F(ti),F(ti+1),F(ti+2),F(ti+3),F(ti+4)代入(5)式得：此外，在t＝t3,t4时，由Simpson公式得：将F(ti),F(ti+1),F(ti+2)代入(7)得在t＝t2时，利用梯形公式可得将F(t1),F(t2)代入(9)得最后考虑到状态t＝t1和t＝tm+1之间的联系，可得x(t1)＝x(tm+1-T)(11)结合(6)、(8)、(10)、(11)式得其中3)由(12)式得系统的Floquet转移矩阵为Ψ＝L-1E4)计算转移矩阵Ψ的特征值，根据Floquet理论通过特征值的模判定系统的稳定性。一般判定方法为：附图说明图1为径向切深与刀具直径之比为1时的稳定性Lobe图，图2为径向切深与刀具直径之比为0.5时的稳定性Lobe图。具体实施方式为了使本专利技术的目的、技术方案及优点更加清楚明白，以下以结合附图及实施例子，对本专利技术进行进一步的详细描述。此处所描述的具体实施例子仅用于解释本专利技术，并不用于限定本专利技术。1.柔性刀具-刚性工件下的x方向上单自由度铣削动力学方程为：式中，mt、ζ、ωn为系统的模态参数，分别表示为刀具的模态质量、阻尼比、刀具的自然圆频率；ap为轴向切削深度；h(t)为切削力系数，即式中，Kt、Kn分别为切向和法向切削力系数，φj(t)为第j个刀齿的角位移，且N为刀具齿数，Ω为刀具主轴转速。g(φj(t))函数定义为：上式中，φst和φex分别表示刀具的切入角和切出角。对于顺铣，φst＝arccos(2a/D-1),φex＝π,逆铣时，φst＝0，φex＝arccos(1-2a/D)，其中a/D表示径向切深与刀具直径之比。令通过变换，则式(1)可改写为：其中令G(t)＝A-H(t),F(t)＝G(t)x(t)+H(t)x(t-T),则(2)变为2.将刀齿切削周期T离散为等距的m个小区间，则任意一个时间小区间为[ti,ti+1]，其中ti＝t1+(i-1)τ(i＝1,2,…m+1)，由(3)式和牛顿-莱布尼兹公式得：利用数值积分中的Bull公式可得:将F(ti),F(ti+1),F(ti+2),F(ti+3),F(ti+4)代入(5)式得：此外，在t＝t3,t4时，由Simpson公式得：将F(ti),F(ti+1),F(ti+2)代入(7)得在t＝t2时，利用梯形公式可得将F(t1),F(t2)代入(9)得最后考虑到状态t＝t1和t＝tm+1之间的联系，可得x(t1)＝x(tm+1-T)(11)结合(6)、(8)、(10)、(11)式得其中3)由(12)式得系统的Floquet转移矩阵为Ψ＝L-1E4)计算转移矩阵Ψ的特征值，根据Floquet理论通过特征值的模判定系统的稳定性。一般判定方法为：给定工艺参数：逆铣加工，刀齿数N＝2，模态质量mt＝0.03993kg，阻尼比ζ＝0.011，自然圆频率ωn＝922×2πrad/s，切削力系数Kt＝6×108N/m2，Kn＝2×108N/m2，离散数m＝50，刀具主轴为转速Ω从5000rpm到10000rpm，轴向切削深度ap从0到4mm。将上述步骤和参数通过Matlab得到稳定性Lobe图来预测铣削过程中的稳定性，选取径向切深与刀具直径之比a/D＝1,0.5，得到稳定性Lobe图如图1，2所示。本文档来自技高网...

【技术保护点】
一种基于Bull公式的铣削稳定性预测方法，其特征在于，包括以下步骤：1)将包含再生效应的铣削过程的动力学方程组通过变换转化为空间状态形式

【技术特征摘要】
1.一种基于Bull公式的铣削稳定性预测方法，其特征在于，包括以下步骤：1)将包含再生效应的铣削过程的动力学方程组通过变换转化为空间状态形式在结构动力学框架下，包括再生效应的铣削过程的动力学模型可由以下时滞微分方程组描述其中，M、C和K分别是刀具的模态质量、模态阻尼和模态刚度矩阵；q(t)为刀具模态坐标，且振型系数在刀尖点处归一，Kc(t)为周期系数矩阵，T为时滞量且等于刀齿切削周期，T＝60/(NΩ),且N为刀具齿数，Ω为刀具主轴转速，单位为rpm；令和x(t)＝[q(t)p(t)]T，通过变换，式(1)可以转换为如下的空间状态形式：其中，A表示常数矩阵；H(t)表示周期为T的考虑再生效应的系数矩阵，且H(t)＝H(t+T)；其中令G(t)＝A-H(t),F(t)＝G(t)x(t)+H(t)x(t-T),则(2)变为2)将刀齿切削周期T离散为等距的m个小区间，则任意一个时间小区间为[ti,ti+1]，其中ti＝t1+(i-1)τ(i＝1,2,…m+1)，由(3)式和牛顿-莱布尼兹公式得：利用数值积分中的Bull公式可得:将F(ti),F(ti+1),F(ti+2),F(ti+3),F(ti+4)代入(5)式得：

【专利技术属性】
技术研发人员：智红英，杜娟，闫献国，张唐圣，陈峙，郭宏，陈高华，曹启超，
申请(专利权)人：太原科技大学，
类型：发明
国别省市：山西,14