Open Access
Issue
JNWPU
Volume 39, Number 1, February 2021
Page(s) 111 - 118
DOI https://doi.org/10.1051/jnwpu/20213910111
Published online 09 April 2021

© 2021 Journal of Northwestern Polytechnical University. All rights reserved.

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

整体叶盘是新一代航空发动机实现结构创新与性能跨越的核心零部件[1-2]。整体叶盘将叶片与轮盘集成一体, 省去了中间榫接结构, 具有减轻质量、简化结构、减小气流损失、提高推重比、改善安全性和可靠性等优点。整体叶盘的叶片型面复杂弯扭大, 叶盘通道窄而深, 且毛坯材料通常为钛合金、高温合金等难加工材料。使其加工难度极大, 被誉为制造业皇冠上的明珠[3-4]。

多轴数控铣削广泛应用于整体叶盘加工, 侧铣常用于整体叶盘型面精加工[5]。由于整体叶盘通道较深, 刀具的悬长较大, 径向刚度较小, 加工过程中容易发生颤振。颤振导致工件表面质量下降、刀具磨损加剧、生产效率受限[6-10]。因此, 抑制颤振对加工过程至关重要。使用稳定性叶瓣图是避免颤振最有效的方法。稳定性叶瓣图为提高加工效率和参数优化提供了理论依据[11-17]。

传统的稳定性叶瓣图不考虑随机变量对切削稳定性的影响[18-21], Duncan等[22]首次研究了随机变量对稳定性叶瓣图的影响; 根据标准偏差, 利用结构参数的平均值、上界和下界确定3条稳定性叶瓣图的曲线; 但没有给出表示不确定度影响的量化指标。为了解决此问题, Liu等[23]研究了随机变量对车削系统颤振稳定性预测的影响, 并首次提出了可靠性叶瓣图。可靠性叶瓣图可用于确定颤振和无颤振区域, 代替了传统的稳定性叶瓣图。颤振可靠性定义为加工系统的稳定性概率, 可用可靠性叶瓣图来确定。

本文应用可靠性叶瓣图来判断侧铣加工系统的稳定性概率。将结构参数和主轴转速视为随机变量, 以单自由度和二自由度铣削的轴向极限切削深度构造临界稳定状态下的性能函数; 采用一次二阶矩法求解铣削加工系统的可靠性模型, 确定系统的稳定性概率。

1 侧铣颤振的动态模型

1.1 侧铣加工动态模型

侧铣加工系统有2个自由度, 如图 1所示。假设侧铣刀有Nt个刀齿, 则动态切屑厚度可以表示为[24]

式中:

Δx(t)=x(t)-x(t-T);Δy(t)=y(t)-y(t-T);T是两齿之间的切削间隔时间; ϕstϕex分别是刀具的切入角与切出角。

X, Y向的切削力可以表示为[24]

式中:

式中, KT, KR分别为切向和径向切削力系数, b为轴向切深。

当刀具旋转时, 方向系数随着时间而变化。[A(t)]是以刀齿切削频率(ω(t)=2πnNt/60)为周期的, 因此可以展开为傅里叶级数。利用傅里叶级数的平均分量, [A(t)]可以近似表示为[A0][24]

其中积分函数可以表示为[24]

动态铣削公式(3)可以简化为[24]

thumbnail 图1

侧铣加工过程的动态模型

1.2 颤振稳定性叶瓣图

刀具-工件接触区域的传递函数矩阵可以表示为[24]

对于单自由度铣削系统, 无颤振时极限轴向切深可以表示为[24]

式中: ωc为给定颤振频率, mXkXcX分别为X向的模态质量、模态刚度和模态阻尼。

对于二自由度铣削系统, 无颤振极限轴向切深可以表示为[24]

式中:

2 铣削颤振可靠性模型

将铣削加工系统的极限状态作为临界稳定状态。"失效"表示不稳定的加工状态, 刀具的振动随着切削力无限地增加。"可靠"表示稳定的加工状态, 刀具的振动会在阻尼作用下减弱。

铣削过程中刀具和工件的相互作用产生颤振。颤振频率与结构参数有关。铣刀的结构参数m, c, k和主轴转速S是铣削系统基本的随机变量。

考虑颤振的铣削系统极限状态函数可以表示为[25-26]

动态铣削过程的可靠性是在给定的时间和参数下不发生颤振的概率, 该可靠性可以定义为[25]

式中: X为随机变量向量, X=(m, c, k, ωc); 结构参数m, c, k是独立的变量且服从正态分布; 颤振频率ωcm, c, k相关, 它的分布规律与随机变量m, c, kS有关。

3 一次二阶矩法

一次二阶矩法将失效面泰勒展开点的性能函数线性化。下面将介绍基于一次二阶矩法的铣削系统可靠性评估方法。

铣削系统的极限状态函数为[25]

对于一个二自由度铣削系统, 基本的随机变量X=(x1, x2xi)T是不相关的, 且服从正态分布, 其中xi分别为mX, mY, cX, cY, kX, kYωc。颤振频率ωcmX, mY, cX, cY, kX, kY相关。铣削系统的相关矩阵可以表示为

式中: ρxixj是变量xixj的相关系数。

标准偏差ζmX, ζmY, ζcX, ζcY, ζkX, ζkY, ζSc分别用ζ1, ζ2, ζ3, ζ4, ζ5, ζ6, ζ7表示。铣削系统的协方差矩阵为

矩阵D是7阶正定对称矩阵。该矩阵有7个实特征值和7个线性无关的正交特征向量。

矩阵B的列向量由矩阵D的正则化特征向量构成。对铣削系统的随机向量X做正交变换[25]

极限状态函数可以表示为独立随机变量Y的函数[25]

基于设计点法, 变量Y的导数为[25]

将平均值X*=(m, c, k, Sc)作为初始设计点, Y*的初始值为[25]

在空间Y中, Z=0是通过点Y*的极限状态切平面。根据独立随机变量线性组合的性质, 可以得到Z的均值和标准差[25]

铣削系统的可靠性指标为[25]

定义变量Yi的敏感度系数为[25]

Y*可以表示为[25]

X*可以表示为[25]

设置误差为ε=10-6, 一直迭代直到2次迭代的偏差小于ε。将β值代入公式(29)可计算得到Y*的值。

铣削系统的可靠性概率定义为[25]

式中,Φ(·)是标准正态分布的累积函数。

4 应用

4.1 整体叶盘侧铣加工实验

实验在五轴加工中心上进行。利用模态冲击锤(PCB 086C03)、加速度传感器(Kistler 8778A500)、振动信号采集卡(NI 9233)和分析软件(Cutpro V9.0)完成侧铣加工系统模态参数的识别, 如图 2所示。

模态实验沿着侧铣刀的XY向进行。X, Y向的频响函数如图 34所示。表 12列出了X, Y向模态参数的多次实验测量结果。侧铣加工系统结构参数平均值和标准差计算结果如表 3所示[27-28]。

thumbnail 图2

振动信号采集过程

thumbnail 图3

侧铣刀X向频响函数

thumbnail 图4

侧铣刀Y向频响函数

表1

侧铣刀X向模态参数测量值

表2

侧铣刀Y向模态参数测量值

表3

侧铣刀模态参数平均值及标准差

4.2 侧铣颤振频率的分布

对于二自由度的侧铣加工系统, 结构参数mX, mY, cX, cY, kX, kY取为实验测量结果的平均值。图 5给出了侧铣加工系统的稳定性叶瓣图和颤振频率图, 稳定性叶瓣图如黑色实线所示, 颤振频率随主轴转速变化图如图蓝色点线所示。

thumbnail 图5

侧铣加工系统稳定性叶瓣图与颤振频率图

4.3 侧铣加工系统的可靠度计算

侧铣加工系统稳定运行条件必须满足fX(X)<0。当已知切削深度时, 侧铣加工系统的可靠度可由一次二阶矩法[24, 26]和蒙特卡罗法得到。初始点设为侧铣加工系统随机变量平均值。一次二阶矩法和蒙特卡罗法得到的可靠度比较结果如表 4所示。结果表明, 一次二阶矩法的计算结果与蒙特卡罗法的计算结果吻合较好。相对误差最大的是7.43%。因此, 一次二阶矩法可以应用于侧铣加工系统颤振可靠度的计算。

表4

给定切削深度主轴转速6 000 r/min的可靠度

4.4 可靠性叶瓣图

可靠度值设置为0.99。可靠性叶瓣图与稳定性叶瓣图对比如图 6所示, 黑色实线为稳定性叶瓣图, 蓝色虚线为可靠性叶瓣图。蓝色虚线下方是可靠区域, 可靠度大于0.99, 蓝色虚线上方是不可靠区域, 可靠度小于0.99。黑色实线下方与蓝色虚线上方的区域是稳定但不可靠区域, 可靠度小于0.99。蓝色虚线下方的区域是即稳定又可靠的区域, 可靠性大于0.99。黑色实线上方的区域是即不稳定也不可靠区域, 可靠性小于0.99。

thumbnail 图6

侧铣加工系统稳定性与可靠性叶瓣图对比图

4.5 实验验证

侧铣验证实验在五轴加工中心上进行, 主轴转速设置为6 000, 8 000和14 000 r/min。切削深度为10~60 mm。当主轴转速为6 000 r/min, 切削深度为20 mm时, 加工表面光滑, 可靠度值为0.999 86。当主轴转速为8 000 r/min, 切削深度为20 mm时, 加工表面呈波纹状, 如图 7所示, 可靠度值趋近于0。

在可靠性叶瓣图上得到了关键可靠性值和极限可靠性值, 如图 8中的黑色圆点所示。当主轴转速为6 000 r/min, 切削深度为26.6 mm时, 可靠度为0.99, 切削深度继续增加时, 可靠度值下降, 直到发生颤振, 可靠度为零; 当主轴转速为15 000 r/min, 切削深度为36.8 mm时, 可靠度为0.99, 随着切削深度的增加, 可靠度值会一直下降到零。验证实验的结果在图 8中标出, 五角形代表稳定点即无颤振发生, 三角形代表不稳定点即发生颤振。实验结果与预测结果保持一致。实验结果表明可靠性叶瓣图可用于预测侧铣加工系统的稳定性概率。

thumbnail 图7

不同转速加工表面示意图

thumbnail 图8

验证实验结果

5 结论

本文建立了考虑随机变量影响的侧铣加工系统颤振可靠性预测模型。通过模态实验获得了侧铣加工系统X向和Y向的结构参数。利用结构参数的平均值绘制了侧铣加工系统的稳定性叶瓣图和颤振频率图。利用一次二阶矩法计算了侧铣加工系统的可靠性指标和可靠性概率, 并与模特卡罗法的计算结果进行了对比。利用可靠性叶瓣图对侧铣加工系统的颤振区域与非颤振区域进行识别。可以得出以下结论:

1) 颤振频率是一个随机变量, 它与结构参数m, c, k和主轴转速S有关;

2) 用一次二阶矩法和蒙特卡罗法对侧铣加工系统的可靠性概率计算结果具有较高的一致性。利用一次二阶矩法可以计算侧铣系统的可靠性概率。

3) 可靠性叶瓣图可较好地预测侧铣加工系统的颤振可靠性; 可靠性叶瓣图可用于侧铣加工系统的稳定性概率估计。

References

  1. Huang Chunfeng. Modern aeroengine integral blisk and its manufacturing technology[J]. Aeronautics and Astronautics, 2004(2): 1–8 [Article] (in Chinese) [Google Scholar]
  2. Ren Junxue, Zhang Dinghua, Wang Zengqiang, et al. Research on the NC machining technique of blisk[J]. Acta Aeronautica et Astronautica Sinica, 2004, 25(2): 205–208 [Article] (in Chinese) [Google Scholar]
  3. Ren Junxue, Jiang Zhennan, Yao Changfeng, et al. Process for 4-axis high efficiency slot plunge milling of open blisk[J]. Acta Aeronautica et Astronautica Sinica, 2008, 29(6): 287–293 [Article] (in Chinese) [Google Scholar]
  4. Ren Junxue, Tian Weijun, Tian Rongxin, et al. A study of the rough milling technique of blisk-tunnel[J]. Mechanical Science and Technology for Aerospace Engineering, 2008, 27(10): 106–110 [Article] (in Chinese) [Google Scholar]
  5. Shi Yaoyao, Duan Jihao, Zhang Junfeng, et al. Blisk disc manufacturing process technology[J]. Aeronautical Manufacturing Technology, 2012, 399(3): 26–31 [Article] (in Chinese) [Google Scholar]
  6. Ren Junxue, Yang Dawang, Yao Changfeng, et al. A cutter-axis vector control method for open blisk 4-axis NC machining based on the control curve[J]. Acta Aeronautica et Astronautica Sinica, 2012, 33(8): 1515–1523 [Article] (in Chinese) [Google Scholar]
  7. Ren Junxue, He Qingong, Yao Chanfeng, et al. Tool axis orientation planning method of fixed axis in each cutting liner for closed blisk tunnel five-axis machining[J]. Acta Aeronautica et Astronautica Sinica, 2012, 33(10): 1923–1930 [Article] (in Chinese) [Google Scholar]
  8. Wang Zengqiang. Machining technology of aeroengine blisk[J]. Aeronautical Manufacturing Technology, 2013, 429(9): 38–43 [Article] (in Chinese) [Google Scholar]
  9. Wang Zengqiang. Key manufacturing technology of aeroengine for large aircraft[J]. Aeronautical Manufacturing Technology, 2013(9): 38–43 [Article] (in Chinese) [Google Scholar]
  10. Zhao P B, Shi Y Y. Experimental research on compound and efficient powerful milling process of blisk[J]. Advanced Materials Research, 2011, 314/315/316: 1150–1153 [Article] [Google Scholar]
  11. Zhang Nan, Shi Yaoyao, Yang Cheng, et al. An instantaneous cutting force model for disc mill cutter based on the machining blisk-tunnel of aero-engine[J]. International Journal of Advanced Manufacturing Technology, 2018, 99(1/2/3/4/5/6/7/8/9/10/11/12/13/14): 233–246 [Google Scholar]
  12. Xin Hongmin, Shi Yaoyao, Zhao Tao. Compound efficient and powerful milling machine tool of blisk[J]. International Journal of Advanced Manufacturing Technology, 2018, 98(9/10/11/12): 1–9 [Article] [Google Scholar]
  13. Zhang Nan, Shi Yaoyao. A 3-D instantaneous cutting force prediction model of indexable disc milling cutter for manufacturing blisk-tunnels considering run-out[J]. The International Journal of Advanced Manufacturing Technology, 2019, 103(9/10/11/12): 4029–4039 [Article] [Google Scholar]
  14. Zhang Nan, Shi Yaoyao. Improvement of cutting force and material removal rate for disc milling TC17 blisk tunnels using GRA-RBF-PSO method[J]. Journal of Mechanical Engineering Science, 2019, 233(16): 5556–5567 [Article] [Google Scholar]
  15. Zhang Nan, Jiang Guangjun, Zhou Jie. Reliability-based interval optimization for the disc-mill cutter machining TC17 blisk-tunnel[J]. International Journal of Performability Engineering, 2019, 15(1): 35–44 [Article] [Google Scholar]
  16. Quintana Guillem, Ciurana Joaquim. Chatter in machining processes: a review[J]. International Journal of Machine Tools & Manufacture, 2011, 51(5): 363–376 [Article] [Google Scholar]
  17. Siddhpura M, Paurobally R. A review of chatter vibration research in turning[J]. International Journal of Machine Tools and Manufacture, 2012, 61: 27–47 [Article] [Google Scholar]
  18. Altinta Yusuf, Budak E. Analytical prediction of stability lobes in milling[J]. CIRP Annals-Manufacturing Technology, 1995, 44(1): 357–362 [Article] [Google Scholar]
  19. Altinta Y, Budak E. Analytical prediction of chatter stability in milling-part I: general formulation[J]. Journal of Dynamic Systems, Measurement, and Control, 1998, 120: 22–30 [Article] [Google Scholar]
  20. Insperger Tama, Stepan Gabor. Semi-discretization method for delayed systems[J]. International Journal for numerical methods in engineering, 2002, 55(5): 503–518 [Article] [Google Scholar]
  21. Ding Ye, Zhu Limin, Zhang Xiaojian, et al. A full-discretization method for prediction of milling stability[J]. International Journal of Machine Tools and Manufacture, 2010, 50(5): 502–509 [Article] [Google Scholar]
  22. Duncan G S, Kurdi M H, Schmitz T L, et al. Uncertainty propagation for selected analytical milling stability limit analyses[C]//34th North American Manufacturing Research Conference, 2006 [Google Scholar]
  23. Liu Yu, Li Tianxiang, Liu Kuo, et al. Chatter reliability prediction of turning process system with uncertainties[J]. Mechanical Systems and Signal Processing, 2016, 66/67: 232–247 [Article] [Google Scholar]
  24. Altintas Yusuf. Manufacturing automation: metal cutting mechanics, machine tool vibrations, and CNC design[M]. Cambridge: Cambridge University Press, 2000: 149–158 [Google Scholar]
  25. Liu Yu, Meng Linlin, Liu Kuo, et al. Chatter reliability of milling system based on first-order second-moment method[J]. International Journal of Advanced Manufacturing Technology, 2016, 87(1/2/3/4): 1–9 [Article] [Google Scholar]
  26. Liu Yu, Wang Zhenyu, Liu Kuo, et al. Chatter stability prediction in milling using time-varying uncertainties[J]. International Journal of Advanced Manufacturing Technology, 2016, 89(9/10/11/12): 1–10 [Article] [Google Scholar]
  27. Mohammadi Yaser, Azvar Milad, Budak Erhan. Suppressing vibration modes of spindle-holder-tool assembly through FRF modification for enhanced chatter stability[J]. Cirp Annals, 2018, 67(1): 397–400 10.1016/j.cirp.2018.03.003 [Google Scholar]
  28. Hu Sen, Huang Xianzhen, Zhang Yimin, et al. Reliability analysis of the chatter stability during milling using a neural network[J]. International Journal of Aerospace Engineering, 2016(1): 1–10 [Article] [Google Scholar]

All Tables

表1

侧铣刀X向模态参数测量值

表2

侧铣刀Y向模态参数测量值

表3

侧铣刀模态参数平均值及标准差

表4

给定切削深度主轴转速6 000 r/min的可靠度

All Figures

thumbnail 图1

侧铣加工过程的动态模型

In the text
thumbnail 图2

振动信号采集过程

In the text
thumbnail 图3

侧铣刀X向频响函数

In the text
thumbnail 图4

侧铣刀Y向频响函数

In the text
thumbnail 图5

侧铣加工系统稳定性叶瓣图与颤振频率图

In the text
thumbnail 图6

侧铣加工系统稳定性与可靠性叶瓣图对比图

In the text
thumbnail 图7

不同转速加工表面示意图

In the text
thumbnail 图8

验证实验结果

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.