, 胡森, 张义民东北大学 机械工程与自动化学院, 辽宁 沈阳 110819
收稿日期：2016-03-31
基金项目：国家自然科学基金资助项目 (51105062, 51305071)；教育部新世纪优秀人才支持计划项目 (NCET-13-0103)；中央高校基本科研业务费专项资金资助项目 (N140304003)。
作者简介：黄贤振 (1982-)，男，山东定陶人，东北大学副教授，博士;
张义民 (1958-)，男，吉林长春人，东北大学教授，教育部“****奖励计划”特聘教授。
摘要：铣削过程中的系统参数往往具有随机性, 严重影响了铣削加工的稳定性.利用BP神经网络综合分析了随机因素对铣削加工过程的影响, 提出了一种铣削加工再生型颤振稳定性可靠度计算方法.建立铣削加工再生型颤振动力学模型, 使用全离散法获取铣削稳定性叶瓣图.利用神经网络拟合极限轴向切深的函数表达式, 再分别使用蒙特卡罗法和一次二阶矩法进行可靠度计算.结果表明, 基于BP神经网络的方法兼具高效和精确的优点.
关键词：再生型颤振铣削稳定性叶瓣图神经网络可靠性
On the Calculation Method for the Reliability of Chatter Stability in CNC Milling
HUANG Xian-zhen
, HU Sen, ZHANG Yi-minSchool of Mechanical Engineering & Automation, Northeastern University, Shenyang 110819, China
Corresponding author: HUANG Xian-zhen, E-mail: xzhhuang@mail.neu.edu.cn
Abstract: The system parameters in the process of milling generally have randomness, which influences the stability of milling. The neural network was used to make a comprehensive analysis of the effects of random factors in milling, and a method for the reliability analysis of regenerative chatter stability in milling was proposed. A dynamic model for regenerative chatters of milling was established, and the stability lobe diagram was obtained by means of full discretization.The neural network was applied to approximate the functional relationship of the limit axial cutting depth, and the reliability was then calculated with the Monte Carlo method and first order second moment method, respectively. Finally, the results of an example were used to demonstrate the efficiency and accuracy of the proposed method based on the neural network.
Key Words: regenerative chattermillingstability lobe diagramneural networkreliability
铣削加工具有高精度、高表面质量以及高材料去除率等优点, 因而广泛应用于航天航空、汽车、船舶、能源、轨道交通等领域.然而近年来, 伴随着现代工业的飞速发展, 机械产品趋向于形状复杂、薄壁结构、难加工材料多、精度要求高等问题使得加工过程中更容易产生振动.轻者影响加工质量, 降低生产效率, 重者造成生产事故的发生, 危及人身安全, 因此有关铣削振动的研究变得极为重要.
颤振是一种不稳定的自激振动, 往往是由系统内部的“再生反馈”引起的.根据反馈方式的不同, 颤振主要分为再生颤振、模态耦合颤振、摩擦颤振等[1].就铣削而言, 再生颤振是引起系统失稳的主要因素[2].文献[3]提出了在频域中解析计算铣削稳定边界的单频率法.单频率法具有较高的计算效率, 但是不能预报小径向切深工况的稳定性.文献[4]随后提出了多频率法解决了这一问题.文献[5]建立了PTP (peak-to-peak) 时域仿真模型, 提出通过仿真计算得到的切削力峰-峰值作为颤振发生的依据.文献[6]改进了时域模型, 将预测的动态切削厚度和静态切削厚度的比值作为颤振发生的判据.文献[7]提出了时域有限元分析法判断铣削稳定性, 该方法可以同步预报铣削稳定性和加工误差.文献[8]在研究周期性变化系数的时滞微分方程稳定性的分析中首次提出了半离散法.文献[9]提出了一种基于直接积分方法的颤振稳定域求解方法, 称为全离散法.
上述方法中的切削动力学参数往往被认为是确定的, 然而实际加工中可以精确描述或控制的因素是不存在的, 大部分因素都是随机、不确定的[10].例如, 电压波动会使工件的进给量具有随机性; 测量、加工等误差会造成刀具几何尺寸的随机性.文献[11]考虑了影响参数的随机性, 采用单频率法判定铣削颤振稳定性, 应用蒙特卡罗法计算出铣削颤振可靠度.但是, 单频率法仅适合多齿刀具和大径向切深工况的稳定性判定, 本文采用全离散法判定铣削颤振稳定性, 并且利用人工神经网络综合分析随机因素对铣削颤振稳定性的影响, 相比蒙特卡罗法大大降低了抽样样本数量, 提高了计算效率.
1 铣削加工颤振稳定性分析常见的铣削加工动力学模型如图 1所示[12], 其中?j为第j齿的瞬时切削角 (rad), Ω为主轴转速 (r/min), ae为径向切深 (mm), vt表示工件相对刀具的进给方向, Ftj和Fnj分别是j齿受到的切向切削力 (N) 和径向切削力 (N).
图 1(Fig. 1)
 | 图 1 铣削加工动力学模型Fig.1 Dynamic model of milling |
图 1所示铣削加工动力学系统的控制方程的状态空间形式可以表示为[9, 13]
 | (1) |
 | (2) |
 | (3) |
对时间t进行离散化处理, 将时间周期T分为r个时间段, 即T= Δt, r∈Z.在任意时间区段上, 式 (1) 可以重新表示为
 | (4) |
 | (5) |
 | (6) |
 | (7) |
解方程 (4), 得到
 | (8) |
 | (9) |
 | (10) |
 | (11) |
 | (12) |
 | (13) |
 | (14) |
 | (15) |
 | (16) |
 | (17) |
2.1 BP神经网络人工神经网络具有良好的自学习能力、非线性映射能力和容错能力, 适合处理对大量数据进行分类、建立复杂的非线性映射等问题.而BP神经网络是人工神经网络最精华、最完美的部分.BP神经网络的设计主要包括输入输出层设计、隐含层设计以及转换函数设计.
2.1.1 输入输出层设计根据全离散法预报颤振稳定性的分析过程, 输入层采用9个神经元, 分别为径向切削力系数Kn, 切向切削力系数Kt, x方向上的模态阻尼cx，模态刚度kx和模态质量mx, y方向上的模态阻尼cy，模态刚度ky和模态质量my以及径向浸入比rc, 输出层为单神经元结构, 代表极限轴向切深aplim.
2.1.2 隐含层设计在三层网络中, 隐含层神经元个数n2和输入层神经元个数n1之间的近似关系为[15]
 | (18) |
2.1.3 数据归一化网络的转换函数为sigmoid函数, 该函数输出曲线两端平坦, 中间部分变化剧烈.因此获得输入和输出变量后, 要对其进行归一化处理, 将数据处理为区间 (0, 1) 之间的数据, 以加快收敛速度.
因此, 本文的BP神经网络模型的拓扑结构如图 2所示.
图 2(Fig. 2)
 | 图 2 BP神经网络拓扑结构Fig.2 Topological structure of BP neural network |
BP网络计算铣削颤振稳定性可靠度的基本思路为:由基本变量产生n组随机样本xi(i=1, 2, …, n), 通过仿真计算得到随机样本的极限轴向切深aplim(xi).利用数据{xi, aplim(xi)}(i=1, 2, …, n) 训练网格, 训练好的网络可以模拟计算出不同变量值的极限轴向切深.在此基础上, 使用蒙特卡罗、一次二阶矩等方法计算可靠度.具体步骤如图 3所示.
图 3(Fig. 3)
 | 图 3 基于BP神经网络的颤振稳定性可靠度计算流程Fig.3 Flow chart for calculating chatter stability reliability by the BP neural network |
2.2 蒙特卡罗法蒙特卡罗法具有通用性强的特点, 并且随着模拟次数的增加, 其计算精度也会提高.该方法的基本思想是:在输入变量的空间内生成随机样本, 然后代入系统模型求出结构响应值, 利用随机样本及其响应的概率特性和统计规律, 推断整体的概率特性和统计规律.基本步骤如下:
1) 根据基本变量的分布特性, 生成n个随机样本xi(i=1, 2, …, n).
2) 将n个随机样本代入可靠性分析模型中, 计算结构响应值g(xi).
3) 统计落入失效域的样本数, 记为nf.
4) 由失效频率近似求得失效概率Pf, 即
 | (19) |
2.3 改进一次二阶矩法实际工程中, 一次二阶矩法是最简单、最常用的方法.一次二阶矩法分为中心点法和设计验算点法.本文采用的是设计验算点法, 即改进一次二阶矩法.
铣削颤振可靠性分析的功能函数为
 | (20) |
设x*=(x1*, x2*, …, xn*)T为极限状态面上的一点, 在点x*处将式 (20) 按泰勒级数展开并取至一次项, 有
 | (21) |
 | (22) |
 | (23) |
 | (24) |
表 1(Table 1)
 表 1 基本变量的概率分布特性Table 1 Probability distribution characteristic of random variables
| 表 1 基本变量的概率分布特性 Table 1 Probability distribution characteristic of random variables |
当主轴转速Ω= 10 740 r/min时, 极限轴向切深最小, 此时发生颤振的可能性最大, 故选取此速度进行稳定性可靠度计算.利用BP神经网络拟合极限轴向切深函数表达式, 实际轴向切深取8.3×10-5 m, 再分别使用蒙特卡罗法和一次二阶矩法计算稳定性可靠度.表 2为不同方法计算可靠度结果的对比, 其中神经网络的训练样本数为500, 蒙特卡罗法的抽样样本数为106.由结果可以看出, 基于BP神经网络的可靠度计算方法在大量降低了抽样样本数量的情况下, 其精度是可以接受的, 其中蒙特卡洛法的精度要远高于改进一次二阶矩法, 这与BP神经网络拟合的极限轴向切深表达式非线性过高有一定的关系.在实际工程中, 可以通过降低实际轴向切深的取值提高铣削颤振稳定性可靠度.例如, 当实际轴向切深取8×10-5 m时, 可靠度会提高到0.993 5.
表 2(Table 2)
表 2 可靠度结果比较Table 2 Comparison of reliability results
| 表 2 可靠度结果比较 Table 2 Comparison of reliability results | ||||||||||
4 结语铣削加工中的振动是限制其优势发挥的主要障碍, 而颤振是其主要形式之一.本文分析了铣削加工再生型颤振稳定性, 基于BP神经网络, 对铣削加工颤振稳定性可靠度计算方法进行了研究, 综合分析了随机因素对铣削颤振稳定性的影响, 更加符合实际工程状况.采用上述方法获取了不同主轴转速下的铣削稳定性可靠度, 对充分发挥铣削优势、提高加工精度以及加工效率具有重要意义.
参考文献
| [1] | 刘强, 李忠群. 数控铣削加工过程仿真与优化:建模、算法与工程应用[M]. 北京: 航空工业出版社, 2010: 57-93. ( Liu Qiang, Li Zhong-qun. Simulation and optimization of CNC milling process[M]. Beijing: Aviation Industry Press, 2010: 57-93.) |
| [2] | Altintas Y, Stn G, Merdol D, et al. Chatter stability of milling in frequency and discrete time domain[J].CIRP Journal of Manufacturing Science & Technology, 2008, 1(1): 35–44. |
| [3] | Altintas Y, 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 |
| [4] | Merdol S D, Altintas Y. Multi frequency solution of chatter stability for low immersion milling[J].Journal of Manufacturing Science & Engineering, 2004, 126(3): 459–466. |
| [5] | Smith S, Tlusty J. Efficient simulation programs for chatter in milling[J].CIRP Annals:Manufacturing Technology, 1993, 42(1): 463–466.DOI:10.1016/S0007-8506(07)62486-X |
| [6] | Campomanes M L, Altintas Y. An improved time domain simulation for dynamic milling at small radial immersions[J].Journal of Manufacturing Science & Engineering, 2003, 125(3): 416–422. |
| [7] | Bayly P V, Halley J E, Mann B P, et al. Stability of interrupted cutting by temporal finite element analysis[J].Journal of Manufacturing Science & Engineering, 2003, 125(2): 220–225. |
| [8] | Insperger T, Stn G. Semi-discretization method for delayed systems[J].International Journal for Numerical Method in Engineering, 2002, 55(5): 503–518.DOI:10.1002/nme.v55:5 |
| [9] | Ding Y, Zhu L M, Zhang X J. A full-discretization method for prediction of milling stability[J].International Journal of Machine Tools & Manufacture, 2010, 50(5): 502–509. |
| [10] | Iglesias A, Munoa J, Ciurana J. Optimisation of face milling operations with structural chatter using a stability model based process planning methodology[J].International Journal of Advanced Manufacturing Technology, 2014, 70(1/2/3/4): 559–571. |
| [11] | 冯吉路, 孙志礼, 孙洪哲, 等. 铣削颤振及其可靠度分析[J].组合机床与自动化加工技术, 2015(12): 8–11. ( Feng Ji-lu, Sun Zhi-li, Sun Hong-zhe, et al. Investigation of the milling chatter and its probability[J].Modular Machine Tool & Automatic Manufacturing Technique, 2015(12): 8–11.) |
| [12] | Ding H, Ding Y, Zhu L M. On time-domain methods for milling stability analysis[J].Chinese Science Bulletin, 2012, 57(33): 4336–4345.DOI:10.1007/s11434-012-5499-y |
| [13] | Gradisek J, Kalveram M, Insperger T, et al. On stability prediction for milling[J].International Journal of Machine Tools & Manufacture, 2005, 45: 769–781. |
| [14] | Farkas M. Periodic motions[M]. New York: Springer-Verlag, 1994: 151-202. |
| [15] | 郁磊, 史峰, 王辉, 等. MATLAB智能算法[M]. 北京: 北京航空航天大学出版社, 2015: 1-37. ( Yu Lei, Shi Feng, Wang Hui, et al. MATLAB intelligent algorithm[M]. Beijing: Beihang University Press, 2015: 1-37.) |
