CN106156477B - 薄壁件动态铣削稳定性叶瓣图高精度预测方法 - Google Patents
薄壁件动态铣削稳定性叶瓣图高精度预测方法 Download PDFInfo
- Publication number
- CN106156477B CN106156477B CN201510205459.XA CN201510205459A CN106156477B CN 106156477 B CN106156477 B CN 106156477B CN 201510205459 A CN201510205459 A CN 201510205459A CN 106156477 B CN106156477 B CN 106156477B
- Authority
- CN
- China
- Prior art keywords
- thin
- milling
- time
- differential equation
- cutting
- 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.)
- Expired - Fee Related
Links
- 238000003801 milling Methods 0.000 title claims abstract description 102
- 238000010586 diagram Methods 0.000 title claims abstract description 12
- 238000013277 forecasting method Methods 0.000 title abstract 2
- 238000000034 method Methods 0.000 claims abstract description 60
- 239000011159 matrix material Substances 0.000 claims abstract description 26
- 238000006073 displacement reaction Methods 0.000 claims abstract description 18
- 238000012546 transfer Methods 0.000 claims abstract description 16
- 230000010354 integration Effects 0.000 claims abstract description 7
- 238000005520 cutting process Methods 0.000 claims description 82
- 238000013016 damping Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000003754 machining Methods 0.000 claims description 8
- 238000013178 mathematical model Methods 0.000 claims description 7
- 125000004432 carbon atom Chemical group C* 0.000 claims description 6
- 230000003068 static effect Effects 0.000 claims description 6
- 101100001676 Emericella variicolor andK gene Proteins 0.000 claims description 3
- 101150089047 cutA gene Proteins 0.000 claims description 3
- 230000000087 stabilizing effect Effects 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims 1
- 238000012545 processing Methods 0.000 description 6
- 238000006243 chemical reaction Methods 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000007796 conventional method Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005316 response function Methods 0.000 description 1
Landscapes
- Numerical Control (AREA)
Abstract
薄壁件动态铣削稳定性叶瓣图高精度预测方法属于动态铣削领域。该方法首先将薄壁件动态铣削过程的动力学模型简化为单自由度二阶时滞微分方程。然后，采用降阶原理将该时滞微分方程降阶为一阶时滞微分方程组。再次，在机床主轴旋转的一个周期内对该微分方程组进行离散。对于每个离散的小区间，将该时滞微分方程组分解为两个部分，分别采用插值法与积分法对其进行简化，得到离散区间首尾时间节点处位移值的递推公式。由此，建立一个周期内最后时间节点的位移与第一个时间节点的位移之间的联系传递矩阵。最后，利用该矩阵与Floquet原理，计算动态铣削系统稳定性叶瓣图。利用本发明提供的方法选择加工工艺参数能够稳定、高效、高精度地加工薄壁件。
Description
技术领域
薄壁件动态铣削稳定性叶瓣图高精度预测方法属于动态铣削领域。
背景技术
数控加工技术是现代制造技术的基础，它的广泛应用使普通机械被数控机械所代替，这使制造业发生了根本的变化，在数控加工中，以加工系统稳定、被加工表面的高精度与高材料去除率为主要特征的高性能加工广泛应用于飞机的大梁、隔板、火箭的整流罩、武器（导弹）的战斗部壳体等大材料去除量的关键薄壁结构和零件的制造中。在动态切削过程中存在刀具切削刃与工件之间的相对振动，这种相对振动的主要表现形式之一为自激颤振，若这种自激颤振控制不当，极易诱发切削系统失稳，进而损坏加工零件，甚至破坏加工设备。 所以对自激颤振进行研究以及对极限切削条件（稳定性叶瓣图）进行预测进而合理的选择加工工艺参数来指导实际加工显得十分必要。张俊、赵万华、卢秉恒等人发明的专利一种高速铣削稳定性快速判定方法，即： CN101905340A；利用耦合得到的机床整机频响函数特性，绘制出主轴转数与轴向切深的稳定性极限图，最后为切削进行稳定性判定。该发明大大提高了判断速度，适合于在线的切削参数优化。然而，该发明也存在缺陷，即：该发明通过频域法计算切削稳定性，造成时域/频域转化过程中的精度丢失，使得该发明损失计算精度而追求计算速度。实际上，如果不进行时域/频域转化，则不存在时/频域转化精度丢失误差，进而可以提高预测精度。
发明内容
本发明要解决的技术难题是克服现有技术的缺陷，采用数值计算方法预测圆周螺旋立铣刀在薄壁件铣削过程中的稳定性叶瓣图，解决传统发明中存在的利用频域描述预测稳定性叶瓣图时/频域转化过程中精度丢失的缺陷，提高动态铣削系统稳定性叶瓣图预测精度。
本发明采用的技术方案是薄壁件动态铣削系统稳定性叶瓣图高精度的数值预测方法。首先，考虑薄壁件动态铣削过程，将该动态过程的动力学模型简化为一单自由度的二阶时滞微分方程。然后，采用降阶原理，将该单自由度的二阶时滞微分方程降阶为一一阶时滞微分方程组。再次，在机床主轴旋转的一个周期内对该一阶时滞微分方程组进行离散，对于每个离散的小区间，将该一阶时滞微分方程组分解为两个部分（Part I 与 Part II）。分别采用插值法与积分法对一阶时滞微分方程组中的Part I 与 Part II进行简化，得到该方程解的数学表达式，进而得到一个离散区间首尾时间节点处位移值的递推公式。由此，由主轴旋转的一个周期内最后的时间节点的位移与第一个时间节点的位移建立联系，得到传递矩阵。最后，利用该传递矩阵与Floquet原理，即可计算出动态铣削系统稳定性叶瓣图。与传统的发明相比，利用本文发明的薄壁件动态铣削系统稳定性叶瓣图的数值预测方法精度更高，根据该预测的稳定性叶瓣图选择合适的工艺参数对薄壁件进行加工，从而能够稳定、高效、高精度的加工薄壁件。具体步骤如下：
1. 利用现有的商业设备获取薄壁件在加工过程中沿着刚度最弱方向的模态参数，包含自然频率ω、模态质量M、模态刚度K和阻尼比ξ以及沿着铣刀切向与径向静态切削力系数K t 与K r ，x向指的是铣刀的进给方向，y向指的是薄壁件刚度最小的方向，即薄壁方向；
2. 构建动态铣削薄壁件过程的数学模型——可以简化为一单自由度的二阶时滞微分方程，该方程可以用下式表示：
其中，M 为薄壁件铣削区域沿着y方向的模态质量，C为薄壁件铣削区域沿着y方向的阻尼，K为薄壁件铣削区域沿着y方向的模态刚度，a为薄壁件铣削区域铣刀轴向切削深度，K c (t)为薄壁件铣削过程中沿着y方向的动态切削力系数，a为轴向切削深度，可以有下式表示：
N为刀具所含切削刃个数，φ j 为第j个切削刃t时刻的位置角，K t 和K r 分别为铣削系统沿着铣刀切向与径向静态切削力系数，g(φ j )函数为一单位跃阶函数，用于确定切削刃是否参与切削，为1时表示此刻第j个切削刃正在参与切削，否则表示第j个切削刃此刻没有参与切削，可以由下式表示：
上式中，φ st 与φ ex 分别表示切削刃的切入切出角，y(t)为薄壁件铣削系统的振动位移；
3. 根据建立的薄壁件动态铣削过程的动力学模型，采用降阶原理，将该动力学模型的二阶时滞微分方程降阶为一一阶时滞微分方程组，该一阶时滞微分方程组可以用下式表示：
具体采用的降阶公式如下：
这时，薄壁件动态铣削数学模型可以由一阶时滞微分方程组表示；
4. 对铣刀旋转一个时间周期进行均匀离散，由于铣刀安装在主轴上，所以铣刀旋转一个周期实际上为主轴旋转一个周期，可以用T表示，具体由下式计算得到：
同时，设离散区间数为m，则每个区间长度可以用h=T/m表示，每个时间节点分别为t 0、t 1、t 2、…、t m ；
5. 在每个小的时间区间内对一阶时滞微分方程组进行求解，例如在第i个时间区间内，该一阶时滞微分方程组的解可用下式表示；
为了书写方便，可以把x(t i+1)简化为x i+1，ε∈（t i t i+1]为积分参数。于是，上式可以简化为：
；
6. 对一阶时滞微分方程组的解进行分部分（第一部分： Part I 与 第二部分：Part II）处理，具体将上式的积分部分分解为两个部分，分别为：
；
7. 对第6步得到的第一部分与第二部分分别采用插值法与积分法进行离散，这里对第一部分采用三阶插值算法，第二部分采用积分算法。对于第一部分：
上式中，ζ∈（0 h], a、b、c和d分别是系数，可以由下式表示：
所以，将上述两式代入第一部分的公式可以得到：
对于第二部分：
由于T是循环周期，所以可以得出：
所以第二部分可以进一步简化：
；
8. 对一阶时滞微分方程组的解计算公式进行整合，对第一部分以及第二部分的公式进行整合并带回到一阶时滞微分方程组的解，得到：
合并同类项得到：
；
9. 对第8步方程进行整合，得到整合方程实际上为一递推公式，对该公式进行变换，得到矩阵递推公式：
；
10. 利用第9步得到的递推矩阵公式，从而可以推得一个周期内Y0与Ym之间的联系，可以由下式表示：
定义D为传递矩阵，计算薄壁件在铣削过程中的稳定性极限图时求解传递矩阵D的特征值，根据Floquet原理，若全部特征值绝对值的最大值小于1，则该系统稳定。同时，在改特征矩阵中，未知数为主轴转速n与轴向切削深度a。给定主轴转速n，并给定传递矩阵D特征值绝对值最大值为1，由此可得对应的最大轴向切削深度a。 进而获取了薄壁件在铣削过程中的切削稳定性叶瓣图。
本发明具有以下明显效果：能够有效的提高薄壁件铣削加工叶瓣图的预测精度。
附图说明
图1，某薄壁零件的局部视图；
图2，薄壁零件动态铣削示意图，其中k y 表示薄壁件沿着y轴方向的模态刚度，c y 为薄壁件沿着y轴方向的模态阻尼，φ j 为第j个切削刃t时刻的位置角，j、j-1表示铣刀的第j个与第j-1个切削刃，“第j个切削刃路径”与“第j-1个切削刃路径”表示第j个与第j-1个切削刃切削过程中运动轨迹；
图3，薄壁件铣削过程中的坐标系定义，其中O表达坐标系的原点，x、y、z分别表示该坐标系的x、y与z轴；
图4，主轴旋转一个周期的离散示意图，T表示主轴旋转的周期，m代表对主轴旋转一个周期进行离散的区间数，t 0、t 1、t 2、... ... t m-2、t m-1、 t m 分别表示m个离散区间的时间节点，x(t 0)、x(t 1 )、x(t 2 )、... ... 、x(t m-2)、x(t m-1)、x(t m )表示对应于时间节点处当前切削刃(设为第j个)的位移值，h表示每个离散小区间的时间长度；
图5，对时滞微分方程组的时间滞后项简化示意图，其中m为对主轴旋转一个周期进行离散的区间数，在该示意图中取m=7，x(t 5-T)表示第5个离散区间右端点对应的时间节点(第6个时间节点)处相对于当前切削刃的前一切削刃(设为第j-1个)的时间滞后位移值，x(t -2)表示相对于当前周期的前一周期倒数第2个时间节点处的当前切削刃的前一切削刃(设为第j-1个)的时间滞后位移值；
图6，利用传统方法预测的薄壁件铣削稳定性叶瓣图，其中由圆圈组成的非连续曲线为实际的薄壁件铣削过程中稳定性叶瓣图，连续曲线为利用传统方法计算的薄壁件铣削过程中的稳定性叶瓣图；
图7，利用本发明提出的方法预测薄壁件铣削过程稳定性极限图, 其中由圆圈组成的非连续曲线为实际的薄壁件铣削过程中稳定性叶瓣图，连续曲线为利用本发明提出的方法计算的薄壁件铣削过程中的稳定性叶瓣图。
具体实施方式
下面结合附图和技术方案，详细说明本发明的具体实施。薄壁件（如图１所示）在数控加工中，常需要优化进给率、定制加工参数等，从而达到加工系统稳定，保证高精度的同时提高铣削效率，这一切的必要条件为高精度的预测薄壁件在加工过程中的稳定性叶瓣图。由于薄壁件铣削动态铣削系统（如图２所示）的动力学模型可以简化为一二阶时间滞后微分方程，不能用解析方法进行预测，为了提高薄壁件铣削系统稳定性叶瓣图的预测精度，本发明提出一种薄壁件动态铣削稳定性叶瓣图高精度预测方法。首先，考虑薄壁件动态铣削过程，将该动态过程的动力学模型简化为一单自由度的二阶时滞微分方程。然后，采用降阶原理，将该单自由度的二阶时滞微分方程降阶为一一阶时滞微分方程组。再次，在机床主轴旋转的一个周期内对该一阶时滞微分方程组进行离散，对于每个离散的小区间，将该一阶时滞微分方程组分解为两个部分（Part I 与 Part II）。分别采用插值法与积分法对一阶时滞微分方程组中的Part I 与 Part II进行简化，得到该方程解的数学表达式，进而得到一个离散区间首尾时间节点处位移值的递推公式。由此，由主轴旋转的一个周期内最后的时间节点的解与第一个时间节点的解建立联系，得到传递矩阵。最后，利用该传递矩阵与Floquet原理，即可计算出动态铣削系统稳定性叶瓣图。该模型预测薄壁件在铣削过程中的稳定性叶瓣图的精度更高；方法的具体步骤如下：
1. 利用现有的商业设备获取薄壁件在加工过程中沿着刚度最弱方向的模态参数，包含自然频率ω、模态质量M、模态刚度K和阻尼比ξ以及沿着铣刀切向与径向静态切削力系数K t 与K r ，x向与y向如图3所示，图3所示的为薄壁件铣削过程中的坐标系定义，其中O表达坐标系的原点，x、y、z分别表示该坐标系的x、y与z轴；
2. 构建动态铣削薄壁件过程的数学模型，如图2所示，图2所示的为薄壁零件动态铣削示意图，其中k y 表示薄壁件沿着y轴方向的模态刚度，c y 为薄壁件沿着y轴方向的模态阻尼，φ j 为第j个切削刃t时刻的位置角，j、j-1表示铣刀的第j个与第j-1个切削刃，“第j个切削刃路径”与“第j-1个切削刃路径”表示第j个与第j-1个切削刃切削过程中运动轨迹；此时，动态铣削薄壁件过程的数学模型可以简化为一单自由度的二阶时滞微分方程，该方程可以用下式表示：
其中，M 为薄壁件铣削区域沿着y方向的模态质量，C为薄壁件铣削区域沿着y方向的阻尼，K为薄壁件铣削区域沿着y方向的模态刚度，a为薄壁件铣削区域铣刀轴向切削深度，K c (t)为薄壁件铣削过程中沿着y方向的动态切削力系数，a为轴向切削深度，可以有下式表示：
N为刀具所含切削刃个数，φ j 为第j个切削刃t时刻的位置角（如图3所示），K t 和K r 分别为铣削系统沿着铣刀切向与径向静态切削力系数，g(φ j )函数为一单位跃阶函数，用于确定切削刃是否参与切削，为1时表示此刻第j个切削刃正在参与切削，否则表示第j个切削刃此刻没有参与切削，可以由下式表示：
上式中，φ st 与φ ex 分别表示切削刃的切入切出角，y(t)为薄壁件铣削系统的振动位移；
3. 根据建立的薄壁件动态铣削过程的动力学模型，采用降阶原理，将该动力学模型的二阶时滞微分方程降阶为一一阶时滞微分方程组，该一阶时滞微分方程组可以用下式表示：
具体采用的降阶公式如下：
这时，薄壁件动态铣削数学模型可以由一阶时滞微分方程组表示；
4. 对铣刀旋转一个时间周期进行均匀离散，如图4所示，图4所示的为主轴旋转一个周期的离散示意图，T表示主轴旋转的周期，m代表对主轴旋转一个周期进行离散的区间数，t 0、t 1、t 2、... ... t m-2、t m-1、 t m 分别表示m个离散区间的时间节点，x(t 0)、x(t 1 )、x(t 2 )、... ... 、x(t m-2)、x(t m-1)、x(t m )表示对应于时间节点处当前切削刃(设为第j个)的位移值，h表示每个离散小区间的时间长度；由于铣刀安装在主轴上，所以铣刀旋转一个周期实际上为主轴旋转一个周期，可以用T表示，具体由下式计算得到：
同时，设离散区间数为m，则每个区间长度可以用h=T/m表示，每个时间节点分别为t 0、t 1、t 2、…、t m ；
5. 在每个小的时间区间内对一阶时滞微分方程组进行求解，例如在第i个时间区间内，该一阶时滞微分方程组的解可用下式表示；
为了书写方便，可以把x(t i+1)简化为x i+1，ε∈（t i t i+1]为积分参数。于是，上式可以简化为：
6. 对一阶时滞微分方程组的解进行分部分（第一部分： Part I 与 第二部分：Part II）处理，具体将上式的积分部分分解为两个部分，分别为：
；
7. 对第6步得到的第一部分与第二部分分别采用插值法与积分法进行离散，这里对第一部分采用三阶插值算法，第二部分采用积分算法。 对于第一部分：
上式中，ζ∈（0 h], a、b、c和d分别是系数，可以由下式表示：
所以，将上述两式代入第一部分的公式可以得到：
对于第二部分：
由于T是循环周期，如图5所示，图5所示的为对时滞微分方程组的时间滞后项简化示意图，其中m为对主轴旋转一个周期进行离散的区间数，在该示意图中取m=7，x(t 5-T)表示第5个离散区间右端点对应的时间节点(第6个时间节点)处相对于当前切削刃的前一切削刃(设为第j-1个)的时间滞后位移值，x(t -2)表示相对于当前周期的前一周期倒数第2个时间节点处（当前切削刃的前一切削刃为第j-1个）的时间滞后位移值；所以可以得出：
所以第二部分可以进一步简化：
；
8. 对一阶时滞微分方程组的解计算公式进行整合，对第一部分以及第二部分的公式进行整合并带回到一阶时滞微分方程组的解，得到：
合并同类项得到：
；
9. 对第8步方程进行整合，得到整合方程实际上为一递推公式，对该公式进行变换，得到矩阵递推公式：
；
10. 利用第9步得到的递推矩阵公式，从而可以推得一个周期内Y0与Ym之间的联系，可以由下式表示：
定义D为传递矩阵，计算薄壁件在铣削过程中的稳定性极限图时求解传递矩阵D的特征值，根据Floquet原理，若全部特征值绝对值的最大值小于1，则该系统稳定。同时，在改特征矩阵中，未知数为主轴转速n与轴向切削深度a。给定主轴转速n，并给定传递矩阵D特征值绝对值的最大值为1，由此可得对应的最大轴向切削深度a。 进而获取了薄壁件在铣削过程中的切削稳定性叶瓣图。
实施例：选用图6，图7所示的实施例来说明本发明的特点：
在图6与图7所示的实施例中，选择的加工参数与刀具参数完全相同，铣刀所含切削刃的个数为N=2，铣刀直径为12mm，径向切削深度为12mm，设离散区间数为m=200的情况下的薄壁件稳定性叶瓣图为真实值。
薄壁件的模态参数与铣削过程中的铣削力系数如分别为:
由此可以算出模态刚度的值:
阻尼为:
分别利用传统方法与本发明提出的方法(m=40)进行计算薄壁件在切削过程中的稳定性叶瓣图，在计算过程中主轴转速范围与铣刀轴向切削深度范围分别为:
经过计算得出两种方法预测的稳定性叶瓣图，图6为传统方法预测的稳定性叶瓣图，其中不连续曲线为真实的叶瓣图曲线，连续曲线为利用传统方法预测的叶瓣图曲线，发现二者不能很好的贴合。
图7为利用本发明提出的预测方法预测的在铣削薄壁件过程中的稳定性叶瓣图与真实稳定性叶瓣图的比较，其中不连续曲线为真实的叶瓣图曲线，连续曲线为利用本发明提出的方法预测的叶瓣图曲线，可以发现二者贴合非常紧密。
本实施例说明本发明提出的方法预测的稳定性叶瓣图精度更高，说明了本发明专利的必要性。
Claims (1)
1.一种薄壁件动态铣削稳定性叶瓣图高精度预测方法，其特征是，首先，考虑薄壁件动态铣削过程，将该动态过程的动力学模型简化为一单自由度的二阶时滞微分方程；然后，采用降阶原理，将该单自由度的二阶时滞微分方程降阶为一一阶时滞微分方程组；再次，在机床主轴旋转的一个周期内对该一阶时滞微分方程组进行离散，对于每个离散的小区间，将该一阶时滞微分方程组分解为两个部分Part I 与 Part II；分别采用插值法与积分法对一阶时滞微分方程组中的Part I 与 Part II进行简化，得到该方程解的数学表达式，进而得到一个离散区间首尾时间节点处位移值的递推公式；由此，由主轴旋转的一个周期内最后的时间节点的位移与第一个时间节点的位移建立联系，得到传递矩阵；最后，利用该传递矩阵与Floquet原理，即可计算出动态铣削系统稳定性叶瓣图；具体步骤如下：
（1）利用现有的商业设备获取薄壁件在加工过程中沿着刚度最弱方向的模态参数，包含自然频率ω、模态质量M、模态刚度K和阻尼比ξ以及沿着铣刀切向与径向静态切削力系数K t 与K r ，x向指的是铣刀的进给方向，y向指的是薄壁件刚度最小的方向，即薄壁方向；
（2） 构建动态铣削薄壁件过程的数学模型——可以简化为一单自由度的二阶时滞微分方程，该方程可以用下式表示：
其中，M 为薄壁件铣削区域沿着y方向的模态质量，C为薄壁件铣削区域沿着y方向的阻尼，K为薄壁件铣削区域沿着y方向的模态刚度，a为薄壁件铣削区域铣刀轴向切削深度，K c (t)为薄壁件铣削过程中沿着y方向的动态切削力系数，a为轴向切削深度，可以有下式表示：
N为刀具所含切削刃个数，φ j 为第j个切削刃t时刻的位置角，K t 与K r 分别为铣削系统沿着铣刀切向与径向静态切削力系数，g(φ j )函数为一单位跃阶函数，用于确定切削刃是否参与切削，为1时表示此刻第j个切削刃正在参与切削，否则表示第j个切削刃此刻没有参与切削，可以由下式表示：
上式中，φ st 与φ ex 分别表示切削刃的切入切出角，y(t)为薄壁件铣削系统的振动位移；
（3） 根据建立的薄壁件动态铣削过程的动力学模型，采用降阶原理，将该动力学模型的二阶时滞微分方程降阶为一一阶时滞微分方程组，该一阶时滞微分方程组可以用下式表示：
具体采用的降阶公式如下：
这时，薄壁件动态铣削数学模型可以由一阶时滞微分方程组表示；
（4） 对铣刀旋转一个时间周期进行均匀离散，由于铣刀安装在主轴上，所以铣刀旋转一个周期实际上为主轴旋转一个周期，可以用T表示，具体由下式计算得到：
T＝60/n
同时，设离散区间数为m，则每个区间长度可以用h=T/m表示，每个时间节点分别为t 0、t 1、t 2、…、t m ；
（5） 在每个小的时间区间内对一阶时滞微分方程组进行求解，在第i个时间区间内，该一阶时滞微分方程组的解可用下式表示；
为了书写方便，可以把x(t i+1)简化为x i+1，ε∈（t i t i+1]为积分参数，于是，上式可以简化为：
;
（6） 对一阶时滞微分方程组的解进行分部分处理，具体将上式的积分部分分解为两个部分，分别为：
；
（7） 对第(6)步得到的第一部分与第二部分分别采用插值法与积分法进行离散，这里对第一部分采用三阶插值算法，第二部分采用积分算法，对于第一部分：
上式中，ζ∈（0 h], a、b、c和d分别是系数，可以由下式表示：
所以，将上述两式代入第一部分的公式可以得到：
对于第二部分：
由于T是循环周期，所以可以得出：
所以第二部分可以进一步简化：
；
（8） 对一阶时滞微分方程组的解计算公式进行整合，对第一部分以及第二部分的公式进行整合并带回到一阶时滞微分方程组的解，得到：
合并同类项得到：
；
（9） 对第(8)步方程进行整合，得到整合方程实际上为一递推公式，对该公式进行变换，得到矩阵递推公式：
Yi+1＝DiYi
Yi+1＝[Xi+1 Xi Xi-1 … Xi-m+1 Xi-m]T
；
（10） 利用第9步得到的递推矩阵公式，从而可以推得一个周期内Y0与Ym之间的联系，可以由下式表示：
Ym＝DY0＝Dm-1…D2D1D0Y0
定义D为传递矩阵，计算薄壁件在铣削过程中的稳定性极限图时求解传递矩阵D的特征值，根据Floquet原理，若全部特征值绝对值的最大值小于1，则该系统稳定；同时，在该传递矩阵中，未知数为主轴转速n与轴向切削深度a；给定主轴转速n，并给定传递矩阵D特征值绝对值的最大值为1，由此可得对应的最大轴向切削深度a；进而获取了薄壁件在铣削过程中的切削稳定性叶瓣图。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510205459.XA CN106156477B (zh) | 2015-04-28 | 2015-04-28 | 薄壁件动态铣削稳定性叶瓣图高精度预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510205459.XA CN106156477B (zh) | 2015-04-28 | 2015-04-28 | 薄壁件动态铣削稳定性叶瓣图高精度预测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106156477A CN106156477A (zh) | 2016-11-23 |
| CN106156477B true CN106156477B (zh) | 2019-11-05 |
Family
ID=57347322
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510205459.XA Expired - Fee Related CN106156477B (zh) | 2015-04-28 | 2015-04-28 | 薄壁件动态铣削稳定性叶瓣图高精度预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106156477B (zh) |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106843147B (zh) * | 2016-12-15 | 2020-03-17 | 太原科技大学 | 一种基于Hamming公式预测铣削稳定性的方法 |
| CN106647625B (zh) * | 2016-12-15 | 2018-12-21 | 太原科技大学 | 一种基于Gear公式预测铣削稳定性的方法 |
| CN106802630B (zh) * | 2017-03-14 | 2019-01-04 | 东北大学 | 基于薄壁件的铣削颤振稳定性预测的并行频域方法 |
| CN106940746A (zh) * | 2017-03-14 | 2017-07-11 | 东北大学 | 基于薄壁件的铣削颤振稳定性预测的并行时域方法 |
| CN106934170B (zh) * | 2017-03-22 | 2019-08-09 | 大连理工大学 | 基于球头铣刀与工件接触区域的颤振稳定域叶瓣图建模方法 |
| CN107346356B (zh) * | 2017-06-27 | 2020-07-03 | 西北工业大学 | 盒型薄壁件铣削稳定性预测方法 |
| CN108519758B (zh) * | 2018-03-21 | 2020-12-22 | 南京工业大学 | 一种避免干涉和颤振的多轴数控加工球头刀轴矢量规划方法 |
| CN108520117B (zh) * | 2018-03-26 | 2020-07-10 | 华中科技大学 | 一种利用全离散法获取稳定性叶瓣图的方法 |
| CN109740264B (zh) * | 2019-01-07 | 2022-08-05 | 南京航空航天大学 | 一种利用牛顿和埃尔米特插值法的铣削稳定域预测方法 |
| CN109909806B (zh) * | 2019-03-22 | 2020-06-19 | 南京理工大学 | 一种提升机器人铣削稳定域的方法 |
| CN112417616A (zh) * | 2020-11-20 | 2021-02-26 | 北京信息科技大学 | 一种铣削加工稳定性预测方法、系统及存储介质 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102873381A (zh) * | 2012-09-29 | 2013-01-16 | 西安交通大学 | 一种基于动力学模型的高速铣削工艺参数优化方法 |
| CN103559550A (zh) * | 2013-09-09 | 2014-02-05 | 西北工业大学 | 多模态耦合下的铣削稳定域预测方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8011864B2 (en) * | 2005-05-26 | 2011-09-06 | University Of Connecticut | Method for facilitating chatter stability mapping in a simultaneous machining application |
| US8296919B2 (en) * | 2007-06-04 | 2012-10-30 | The Boeing Company | Increased process damping via mass reduction for high performance milling |
-
2015
- 2015-04-28 CN CN201510205459.XA patent/CN106156477B/zh not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102873381A (zh) * | 2012-09-29 | 2013-01-16 | 西安交通大学 | 一种基于动力学模型的高速铣削工艺参数优化方法 |
| CN103559550A (zh) * | 2013-09-09 | 2014-02-05 | 西北工业大学 | 多模态耦合下的铣削稳定域预测方法 |
Non-Patent Citations (6)
| Title |
|---|
| An improved prediction of stability lobes using nonlinear thin wall dynamics;O.B. Adetoro et al;《Journal of Materials Processing Technology》;20100401;第210卷(第6-7期);969-979页 * |
| Spectral Method for Prediction of Chatter Stability in Low Radial Immersion Milling;Ye Ding et al;《2011 IEEE International Conference on Robotics and Automation》;20110513;4359-4363页 * |
| Stability Lobes Prediction in Thin Wall Machining;O.B.Adetoro et al;《Proceedings of the World Congress on Engineering 2009》;20090703;第I卷;1-6页 * |
| Toolpath dependent stability lobes for the milling of thin-walled parts;Sébastien Seguy et al;《 International Journal of Machining and Machinability of Materials》;20081231;第4卷(第4期);377-392页 * |
| 铣削加工过程稳定性分析;杨建中等;《中国机械工程》;20130228;第24卷(第3期);360-365页 * |
| 铣削过程颤振稳定域建模研究现状与展望;李忠群等;《制造技术与机床》;20131015(第10期);45-48,53页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106156477A (zh) | 2016-11-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106156477B (zh) | 薄壁件动态铣削稳定性叶瓣图高精度预测方法 | |
| CN107457609B (zh) | 基于刚度变化的铣削颤振抑制方法及铣削颤振优化系统 | |
| CN107480352B (zh) | 一种铣削加工工艺参数的可靠性优化方法 | |
| CN113168491B (zh) | 一种无颤振铣削加工表面形貌仿真方法 | |
| CN109909806B (zh) | 一种提升机器人铣削稳定域的方法 | |
| Jin et al. | 3D stability lobe considering the helix angle effect in thin-wall milling | |
| Comak et al. | Stability of milling operations with asymmetric cutter dynamics in rotating coordinates | |
| CN111177860A (zh) | 一种提升钛合金薄壁件铣削稳定域的方法 | |
| CN108804749B (zh) | 薄壁件铣削过程阻尼建模方法 | |
| CN105414616A (zh) | 螺旋铣孔过程中切削力预报和稳定性判别方法 | |
| CN110188311A (zh) | 基于刀齿切削时程精细积分的高速加工稳定域预测方法 | |
| Liu et al. | Chatter stability prediction in milling using time-varying uncertainties | |
| CN110750891B (zh) | 一种平行同步正交车铣颤振稳定性叶瓣图预测方法 | |
| CN107045578B (zh) | 一种基于npso算法的容屑槽加工砂轮位姿求解方法 | |
| CN110162733B (zh) | 基于整体离散策略的铣削稳定性分析方法 | |
| Shi et al. | A chatter mitigation technique in milling based on H∞-ADDPMS and piezoelectric stack actuators | |
| Zhao et al. | A framework for accuracy enhancement in milling thin-walled narrow-vane turbine impeller of NiAl-based superalloy | |
| Esfandi et al. | Robot assisted machining of thin-walled structures | |
| CN111176209A (zh) | 型腔螺旋铣削加工进给率与转速离线规划方法 | |
| CN109740264B (zh) | 一种利用牛顿和埃尔米特插值法的铣削稳定域预测方法 | |
| Ma et al. | Effect of cutter runout on chatter stability of milling process | |
| CN114309744A (zh) | 一种改善机器人铣边稳定性的方法 | |
| CN109702554A (zh) | 一种自适应振动塑形铣削颤振抑制方法 | |
| Chen et al. | FEM-based prediction and control of milling deformation for a thin-wall web of Ti-6Al-4V alloy | |
| CN111299668B (zh) | 一种不等齿距铣刀的齿间角确定方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20191105 Termination date: 20210428 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |