2. 内蒙古工业大学 机械工程学院, 内蒙古 呼和浩特 010062;
3. 内蒙古工业大学 航空学院, 内蒙古 呼和浩特 010062
整体叶盘是新一代航空发动机实现结构创新与性能跨越的核心零部件[1-2]。整体叶盘将叶片与轮盘集成一体, 省去了中间榫接结构, 具有减轻质量、简化结构、减小气流损失、提高推重比、改善安全性和可靠性等优点。整体叶盘的叶片型面复杂弯扭大, 叶盘通道窄而深, 且毛坯材料通常为钛合金、高温合金等难加工材料。使其加工难度极大, 被誉为制造业皇冠上的明珠[3-4]。
多轴数控铣削广泛应用于整体叶盘加工, 侧铣常用于整体叶盘型面精加工[5]。由于整体叶盘通道较深, 刀具的悬长较大, 径向刚度较小, 加工过程中容易发生颤振。颤振导致工件表面质量下降、刀具磨损加剧、生产效率受限[6-10]。因此, 抑制颤振对加工过程至关重要。使用稳定性叶瓣图是避免颤振最有效的方法。稳定性叶瓣图为提高加工效率和参数优化提供了理论依据[11-17]。
传统的稳定性叶瓣图不考虑随机变量对切削稳定性的影响[18-21], Duncan等[22]首次研究了随机变量对稳定性叶瓣图的影响; 根据标准偏差, 利用结构参数的平均值、上界和下界确定3条稳定性叶瓣图的曲线; 但没有给出表示不确定度影响的量化指标。为了解决此问题, Liu等[23]研究了随机变量对车削系统颤振稳定性预测的影响, 并首次提出了可靠性叶瓣图。可靠性叶瓣图可用于确定颤振和无颤振区域, 代替了传统的稳定性叶瓣图。颤振可靠性定义为加工系统的稳定性概率, 可用可靠性叶瓣图来确定。
本文应用可靠性叶瓣图来判断侧铣加工系统的稳定性概率。将结构参数和主轴转速视为随机变量, 以单自由度和二自由度铣削的轴向极限切削深度构造临界稳定状态下的性能函数; 采用一次二阶矩法求解铣削加工系统的可靠性模型, 确定系统的稳定性概率。
1 侧铣颤振的动态模型 1.1 侧铣加工动态模型侧铣加工系统有2个自由度, 如图 1所示。假设侧铣刀有Nt个刀齿, 则动态切屑厚度可以表示为[24]
(1) |
式中:
(2) |
Δx(t)=x(t)-x(t-T);Δy(t)=y(t)-y(t-T);T是两齿之间的切削间隔时间; ϕst和ϕex分别是刀具的切入角与切出角。
X, Y向的切削力可以表示为[24]
(3) |
式中:
(4) |
(5) |
式中, KT, KR分别为切向和径向切削力系数, b为轴向切深。
当刀具旋转时, 方向系数随着时间而变化。[A(t)]是以刀齿切削频率(ω(t)=2πnNt/60)为周期的, 因此可以展开为傅里叶级数。利用傅里叶级数的平均分量, [A(t)]可以近似表示为[A0][24]
(6) |
其中积分函数可以表示为[24]
(7) |
动态铣削公式(3)可以简化为[24]
(8) |
刀具-工件接触区域的传递函数矩阵可以表示为[24]
(9) |
对于单自由度铣削系统, 无颤振时极限轴向切深可以表示为[24]
(10) |
式中: ωc为给定颤振频率, mX、kX、cX分别为X向的模态质量、模态刚度和模态阻尼。
对于二自由度铣削系统, 无颤振极限轴向切深可以表示为[24]
(11) |
式中:
(12) |
(13) |
将铣削加工系统的极限状态作为临界稳定状态。"失效"表示不稳定的加工状态, 刀具的振动随着切削力无限地增加。"可靠"表示稳定的加工状态, 刀具的振动会在阻尼作用下减弱。
铣削过程中刀具和工件的相互作用产生颤振。颤振频率与结构参数有关。铣刀的结构参数m, c, k和主轴转速S是铣削系统基本的随机变量。
(14) |
动态铣削过程的可靠性是在给定的时间和参数下不发生颤振的概率, 该可靠性可以定义为[25]
(15) |
式中: X为随机变量向量, X=(m, c, k, ωc); 结构参数m, c, k是独立的变量且服从正态分布; 颤振频率ωc与m, c, k相关, 它的分布规律与随机变量m, c, k和S有关。
3 一次二阶矩法一次二阶矩法将失效面泰勒展开点的性能函数线性化。下面将介绍基于一次二阶矩法的铣削系统可靠性评估方法。
铣削系统的极限状态函数为[25]
(16) |
对于一个二自由度铣削系统, 基本的随机变量X=(x1, x2…xi)T是不相关的, 且服从正态分布, 其中xi分别为mX, mY, cX, cY, kX, kY和ωc。颤振频率ωc与mX, mY, cX, cY, kX, kY相关。铣削系统的相关矩阵可以表示为
(17) |
式中: ρxixj是变量xi和xj的相关系数。
标准偏差ζmX, ζmY, ζcX, ζcY, ζkX, ζkY, ζSc分别用ζ1, ζ2, ζ3, ζ4, ζ5, ζ6, ζ7表示。铣削系统的协方差矩阵为
(18) |
矩阵D是7阶正定对称矩阵。该矩阵有7个实特征值和7个线性无关的正交特征向量。
矩阵B的列向量由矩阵D的正则化特征向量构成。对铣削系统的随机向量X做正交变换[25]
(19) |
(20) |
(21) |
极限状态函数可以表示为独立随机变量Y的函数[25]
(22) |
基于设计点法, 变量Y的导数为[25]
(23) |
将平均值X*=(m, c, k, Sc)作为初始设计点, Y*的初始值为[25]
(24) |
在空间Y中, Z=0是通过点Y*的极限状态切平面。根据独立随机变量线性组合的性质, 可以得到Z的均值和标准差[25]
(25) |
(26) |
铣削系统的可靠性指标为[25]
(27) |
定义变量Yi的敏感度系数为[25]
(28) |
新Y*可以表示为[25]
(29) |
新X*可以表示为[25]
(30) |
设置误差为ε=10-6, 一直迭代直到2次迭代的偏差小于ε。将β值代入公式(29)可计算得到Y*的值。
铣削系统的可靠性概率定义为[25]
(31) |
式中,Φ(·)是标准正态分布的累积函数。
4 应用 4.1 整体叶盘侧铣加工实验实验在五轴加工中心上进行。利用模态冲击锤(PCB 086C03)、加速度传感器(Kistler 8778A500)、振动信号采集卡(NI 9233)和分析软件(Cutpro V9.0)完成侧铣加工系统模态参数的识别, 如图 2所示。
模态实验沿着侧铣刀的X和Y向进行。X, Y向的频响函数如图 3和4所示。表 1和2列出了X, Y向模态参数的多次实验测量结果。侧铣加工系统结构参数平均值和标准差计算结果如表 3所示[27-28]。
序号 | fnX/Hz | kX/(N·m-1) | mX/kg | cX/(Ns·m-1) |
1 | 1 022 | 4.33×106 | 0.204 38 | 58.43 |
2 | 1 025 | 4.08×106 | 0.198 75 | 58.75 |
3 | 1 029 | 4.27×106 | 0.210 94 | 58.60 |
4 | 1 020 | 4.41×106 | 0.209 35 | 58.14 |
5 | 1 025 | 4.27×106 | 0.190 53 | 58.95 |
6 | 1 022 | 4.39×106 | 0.200 60 | 59.17 |
7 | 1 023 | 4.51×106 | 0.211 83 | 58.49 |
8 | 1 026 | 4.40×106 | 0.192 06 | 59.07 |
9 | 1 025 | 4.37×106 | 0.215 39 | 58.42 |
10 | 1 028 | 4.50×106 | 0.206 57 | 58.68 |
序号 | fnY/Hz | kY/(N·m-1) | mY/kg | cY/(Ns·m-1) |
1 | 1 000 | 8.33×106 | 0.384 09 | 106.43 |
2 | 1 000 | 8.26×106 | 0.392 15 | 107.75 |
3 | 1 001 | 8.21×106 | 0.410 33 | 105.60 |
4 | 1 002 | 8.40×106 | 0.405 14 | 107.14 |
5 | 1 000 | 8.37×106 | 0.395 90 | 106.95 |
6 | 1 001 | 8.30×106 | 0.400 80 | 105.17 |
7 | 1 001 | 8.41×106 | 0.385 72 | 106.49 |
8 | 1 002 | 8.48×106 | 0.391 09 | 107.07 |
9 | 1 002 | 8.36×106 | 0.427 60 | 107.42 |
10 | 1 001 | 8.45×106 | 0.405 29 | 106.68 |
参数 | fn/Hz | k/(N·m-1) | m/kg | c/(Ns·m-1) | |
X | μX | 1 024.5 | 4.35×106 | 0.204 04 | 58.67 |
σX | 2.66 | 1.19×105 | 0.084 20 | 0.31 | |
Y | μY | 1 001.00 | 8.36×106 | 0.399 81 | 106.67 |
σY | 0.78 | 0.80×105 | 0.013 10 | 0.79 |
对于二自由度的侧铣加工系统, 结构参数mX, mY, cX, cY, kX, kY取为实验测量结果的平均值。图 5给出了侧铣加工系统的稳定性叶瓣图和颤振频率图, 稳定性叶瓣图如黑色实线所示, 颤振频率随主轴转速变化图如图蓝色点线所示。
4.3 侧铣加工系统的可靠度计算侧铣加工系统稳定运行条件必须满足fX(X)<0。当已知切削深度时, 侧铣加工系统的可靠度可由一次二阶矩法[24, 26]和蒙特卡罗法得到。初始点设为侧铣加工系统随机变量平均值。一次二阶矩法和蒙特卡罗法得到的可靠度比较结果如表 4所示。结果表明, 一次二阶矩法的计算结果与蒙特卡罗法的计算结果吻合较好。相对误差最大的是7.43%。因此, 一次二阶矩法可以应用于侧铣加工系统颤振可靠度的计算。
序号 | Pr(FOSM) | Pr(Monte Carlo) | 相对误差/% |
1 | 1.000 000 | 1.000 000 | 0 |
2 | 1.000 000 | 1.000 000 | 0 |
3 | 1.000 000 | 1.000 000 | 0 |
4 | 0.999 990 | 0.999 988 | 0.000 2 |
5 | 0.998 969 | 0.998 931 | 0.003 8 |
6 | 0.987 965 | 0.987 852 | 0.011 4 |
7 | 0.786 742 | 0.786 976 | -0.029 7 |
8 | 0.365 226 | 0.365 798 | -0.157 0 |
9 | 0.087 965 | 0.087 674 | 0.330 8 |
10 | 0.003 578 | 0.003 312 | 7.434 3 |
可靠度值设置为0.99。可靠性叶瓣图与稳定性叶瓣图对比如图 6所示, 黑色实线为稳定性叶瓣图, 蓝色虚线为可靠性叶瓣图。蓝色虚线下方是可靠区域, 可靠度大于0.99, 蓝色虚线上方是不可靠区域, 可靠度小于0.99。黑色实线下方与蓝色虚线上方的区域是稳定但不可靠区域, 可靠度小于0.99。蓝色虚线下方的区域是即稳定又可靠的区域, 可靠性大于0.99。黑色实线上方的区域是即不稳定也不可靠区域, 可靠性小于0.99。
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中标出, 五角形代表稳定点即无颤振发生, 三角形代表不稳定点即发生颤振。实验结果与预测结果保持一致。实验结果表明可靠性叶瓣图可用于预测侧铣加工系统的稳定性概率。
5 结论本文建立了考虑随机变量影响的侧铣加工系统颤振可靠性预测模型。通过模态实验获得了侧铣加工系统X向和Y向的结构参数。利用结构参数的平均值绘制了侧铣加工系统的稳定性叶瓣图和颤振频率图。利用一次二阶矩法计算了侧铣加工系统的可靠性指标和可靠性概率, 并与模特卡罗法的计算结果进行了对比。利用可靠性叶瓣图对侧铣加工系统的颤振区域与非颤振区域进行识别。可以得出以下结论:
1) 颤振频率是一个随机变量, 它与结构参数m, c, k和主轴转速S有关;
2) 用一次二阶矩法和蒙特卡罗法对侧铣加工系统的可靠性概率计算结果具有较高的一致性。利用一次二阶矩法可以计算侧铣系统的可靠性概率。
3) 可靠性叶瓣图可较好地预测侧铣加工系统的颤振可靠性; 可靠性叶瓣图可用于侧铣加工系统的稳定性概率估计。
[1] |
黄春峰. 现代航空发动机整体叶盘及其制造技术[J]. 航空与航天, 2004(2): 1-8.
HUANG Chunfeng. Modern aeroengine integral blisk and its manufacturing technology[J]. Aeronautics and Astronautics, 2004(2): 1-8. (in Chinese) |
[2] |
任军学, 张定华, 王增强, 等. 整体叶盘数控加工技术研究[J]. 航空学报, 2004, 25(2): 205-208.
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. (in Chinese) |
[3] |
任军学, 姜振南, 姚倡锋, 等. 开式整体叶盘四坐标高效开槽插铣工艺方法[J]. 航空学报, 2008, 29(6): 287-293.
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. (in Chinese) |
[4] |
任军学, 田卫军, 田荣鑫, 等. 开式整体叶盘通道侧铣粗加工技术的研究[J]. 机械科学与技术, 2008, 27(10): 106-110.
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. (in Chinese) |
[5] |
史耀耀, 段继豪, 张军锋, 等. 整体叶盘制造工艺技术综述[J]. 航空制造技术, 2012, 399(3): 26-31.
SHI Yaoyao, DUAN Jihao, ZHANG Junfeng, et al. Blisk disc manufacturing process technology[J]. Aeronautical Manufacturing Technology, 2012, 399(3): 26-31. (in Chinese) |
[6] |
任军学, 杨大望, 姚倡锋, 等. 基于控制线的开式整体叶盘叶片四轴数控加工刀轴控制方法[J]. 航空学报, 2012, 33(8): 1515-1523.
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. (in Chinese) |
[7] |
任军学, 何卿功, 姚倡锋, 等. 闭式整体叶盘通道五坐标分行定轴加工刀轴矢量规划方法[J]. 航空学报, 2012, 33(10): 1923-1930.
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. (in Chinese) |
[8] |
王增强. 航空发动机整体叶盘加工技术[J]. 航空制造技术, 2013, 429(9): 38-43.
WANG Zengqiang. Machining technology of aeroengine blisk[J]. Aeronautical Manufacturing Technology, 2013, 429(9): 38-43. (in Chinese) |
[9] |
王增强. 大飞机发动机关键制造技术: 航空发动机整体叶盘加工技术[J]. 航空制造技术, 2013(9): 38-43.
WANG Zengqiang. Key manufacturing technology of aeroengine for large aircraft[J]. Aeronautical Manufacturing Technology, 2013(9): 38-43. (in Chinese) |
[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. |
[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. |
[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. DOI:10.1007%2Fs00170-018-2225-1 |
[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. DOI:10.1007/s00170-019-03780-5 |
[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. DOI:10.1177/0954406219848473 |
[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. |
[16] | QUINTANA Guillem, CIURANA Joaquim. Chatter in machining processes: a review[J]. International Journal of Machine Tools & Manufacture, 2011, 51(5): 363-376. |
[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. DOI:10.1016/j.ijmachtools.2012.05.007 |
[18] | ALTINTA Yusuf, BUDAK E. Analytical prediction of stability lobes in milling[J]. CIRP Annals-Manufacturing Technology, 1995, 44(1): 357-362. DOI:10.1016/S0007-8506(07)62342-7 |
[19] | ALTINTA Y, BUDAK E. Analytical prediction of chatter stability in milling-part Ⅰ: general formulation[J]. Journal of Dynamic Systems, Measurement, and Control, 1998, 120: 22-30. DOI:10.1115/1.2801317 |
[20] | INSPERGER Tama, STEPAN Gabor. Semi-discretization method for delayed systems[J]. International Journal for numerical methods in engineering, 2002, 55(5): 503-518. DOI:10.1002/nme.505 |
[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. DOI:10.1016/j.ijmachtools.2010.01.003 |
[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 |
[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. DOI:10.1016/j.ymssp.2015.06.030 |
[24] | ALTINTAS Yusuf. Manufacturing automation: metal cutting mechanics, machine tool vibrations, and CNC design[M]. Cambridge: Cambridge University Press, 2000: 149-158. |
[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. DOI:10.1007/s00170-016-8523-6 |
[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. |
[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. DOI:10.1016/j.cirp.2018.03.003 |
[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. |
2. School of Mechanical Engineering, Inner Mongolia University of Technology, Hohhot 010051, China;
3. School of Aeronautics, Inner Mongolia University of Technology, Hohhot 010051, China