资源描述:
2014 年 12 月 第 42 卷 第 23 期 机床与液压 MACHINE TOOL & HYDRAULICS Dec 2014 Vol 42 No 23 DOI10.3969/ j issn 1001-3881 2014 23 001 收稿日期 2013-10-15 基金项目 国家自然科学基金资助项目 (51275305); 高等学校博士学科点专项科研基金资助项目 (20110073110041); 国 家科技重大专项项目 (2011ZX04015-031) 作者简介 张景然 (1988), 女, 硕士研究生, 研究方向为数控机床精密加工与测试。 E - mail wenzhijingran@ sj⁃ tu edu cn。 基于模拟退火遗传算法优化 BP 网络的数控机床温度 布点优化及热误差建模 张景然, 沈牧文, 杨建国 (上海交通大学机械与动力工程学院, 上海 200240) 摘要 在热误差建模中, 温度测点的优化选择至关重要。 提出了运用相关性方法, 分析测点温度与主轴热漂移之间的 关系, 找到相关性较高的测点位置, 实现温度布点的优化选择。 在此基础上采用模拟退火遗传算法 (GSA) 优化 BP 神经 网络的方法建立热误差模型, 并通过实验验证。 结果表明 优化的热误差模型能够跳出局部最优而达到全局最优解, 得到 的误差模型拟合值更加接近实测误差值; 基于 GSA 优化的 BP 网络模型较传统的神经网络模型有较高的精度及更强鲁棒性。 关键词 数控机床; 温度布点优化; 主轴热漂移; 热误差建模 中图分类号 TH161 文献标识码 A 文章编号 1001-3881 (2014) 23-001-4 CNC Machine Tool Temperature Measuring Point Optimization and Thermal Error Modeling Based on Stimulated Annealing and Optimized BP Network by Genetic Algorithm ZHANG Jingran, SHEN Muwen, YANG Jianguo (School of Mechanical Engineering, Shanghai Jiaotong University, Shanghai 200240, China) Abstract Optimization selection of the temperature measuring point is crucial during thermal error modeling. A method using correlation analysis was present to analyze the relationship between the spindle thermal drift and point temperature of measurement. The temperature distribution points of the optimal choice were achieved by finding a higher correlation measuring point of location. On basis of this, by using simulated annealing and genetic algorithm (GSA) optimized BP neural network method, the thermal error model was established, and was verified by experiments. The results show optimized thermal error model can escape from local optimal and achieve global optimal solution. The resulting error model can fit values more closer to the actual measured error values. Based on simulated an⁃ nealing genetic algorithm (GSA) optimization, BP neural network model has higher accuracy and greater robustness than that of the traditional neural network model. Keywords CNC machine tools; Temperature measuring point optimization; Thermal drift of spindle; Thermal error modeling 0 前言 随着现代工业的不断发展, 制造业对于机械产品 精度的要求越来越高, 机械产品的精度很大一部分取 决于数控机床的加工精度, 在众多引起的数控机床加 工误差中, 热误差可达总误差的 40% ~80%[1]。 数控 机床的误差补偿可能获得比母机更高的精度[2], 有效 的补偿取决于误差数学模型。 热误差建模方法有许多 种, 如粗集方法、 最小二乘法、 遗传算法、 神经网络 法以及多种算法协同优化的建模方法。 本文作者利用相关性分析温度值和误差值之间的 相互联系程度, 找出最优的测点, 减少温度测点的数 目, 提高了建模效率。 利用模拟退火遗传算法优化 BP 神经网络[3]对数控机床热误差进行建模, 消除单 一建模方式的弊端, 跳出了神经网络建模容易陷入局 部最优解的情况, 达到全局最优, 同时提高热误差建 模的精确性和鲁棒性。 1 基于模拟退火遗传算法 (GSA) 改进的 BP 网 络 1 1 模拟退火算法 (SA) 的原理 模拟退火算法是基于对固体退火降温过程而建立 起的一种全局优化的方法。 加热固体至充分高, 然后让其慢慢冷却, 升温 时, 固体内能增加, 内部粒子变为无序状, 慢慢冷却 时, 粒子渐渐变得有序, 最终在每一个温度下都能达 到平衡的状态, 在常温时达到基态。 根据 Metropolis 准则, 固体 粒 子 在 温 度 T 时 达 到 平 衡 的 概 率 为 e -ΔE/ KT, E 表示温度 T 时的内能, ΔE 为内能的增减 量, K 为 Boltzmann 常数。 在实际应用中, 将内能 E 近似为目标函数 f, 温度 T 则近似为控制参数, 从给 定解开始, 在其邻域中产生新解, 算法允许目标函数 在一定范围接受使目标函数恶化解, 由较高温度初始 值开始, 持续进行 “产生新解计算目标函数差 判断是否接受新解接受或者舍弃” 的迭代, 经过 大量计算, 结合突跳性概率求得给定 T 的目标函数 的全局最优解。 1 2 遗传算法 (GA) 的原理 遗传算法源于生物学的新达尔文主义, 模仿自然 选择和基因进化论。 GA 是基于自然选择以及自然遗 传的搜索算法, 把最优问题的搜索空间映射为遗传空 间, 首先随机产生初始群体, 根据目标函数的要求选 择适应度函数计算个体对环境的适应能力, 然后按照 个体适应度对染色体进行选择, 交叉、 变异得到新一 代群体。 因此遗传算法对算子的基本操作分为选择、 交叉和变异三步。 如此迭代、 进化, 最终得到满足目 标函数要求和收敛条件的适应度最高的群体, 从而得 到最优解[4]。 1 3 BP 神经网络的原理 BP 神经网络是按照误差传播算法训练的多层前 馈网络, 利用输出后的误差来估计输出层直接前导层 的误差, 再用这个误差估计更前一层的误差, 如此一 层一层反传, 就获得了各层的误差估计。 神经网络的 学习过程是在外界输入样本的刺激下不断改变网络连 接权值, 从而使网络模型的实际输出不断接近期望的 输出。 BP 神经网络算法其实质是以误差的平方和为 目标函数, 使用梯度法通过反向传播不断调整网络权 值和阈值, 最终求其误差平方和的最小值。 BP 神经 网络结构包括输入样本输入层各隐层输出层, 最后使网络输出的误差值减少到可接受的程度。 1 4 基于模拟退火遗传算法优化的 BP 神经网络 BP 神经网络采用误差反传算法有一定的局限性, 用梯度法求其最小值实质是一个无约束的非线性最优 化过程, 在网络结构较大时, 不仅计算时间较长, 除 非误差平方和的函数正定, 否则很容易陷入局部最小 值而得不到全局最优解。 模拟退火遗传算法 (GSA) 是在基本遗传算法的编码和优点基础上, 引入模拟退 火算法寻优搜索过程, SA 则对 GA 的种群进行进一 步的优化, 当温度 T 充分大时, 相应的接受概率趋近 于 1, SA 赋予优化过程较高的概率突跳性, 此时算 法在全局进行搜索, 当温度 T 较小时, 相应的概率趋 近于 0, 则进行局部搜索。 因而 GSA 既能解决 BP 网 络易陷入局部最优而得不到全局最优这一问题, 还能 使相应的收敛速度得到改善[5]。 基于 GSA 优化的 BP 神经网络模型, 主要利用 GSA 算法取代神经网络的梯度法, 作为神经网络训 练方法。 分为创建 BP 网络和模拟退火遗传算法的权 值和阈值的选定。 2 数控机床温度布点优化及热误差测量 2 1 温度及热误差测量 研究对象是在国产 ETC3650 数控车削中心, 布置 8 个温度传感器测量温度, 分别是 T1(主轴表面), T2 (主轴前端面), T3(主轴前轴承), T4(主轴后轴承), T5(主轴箱内侧), T6(主轴箱外壳), T7(导轨), T8 (丝杠)。 实验过程中, 主轴转速 n 为 2 000 r/ min, 采 用 API 动态测试仪测试数控车削中心主轴径向 (两个 方向在图 3 中分别记录为 x1、 y1)、 主轴轴向 (z) 三 个方向的热漂移误差值。 记录下 0~120 min 的数据, 每分钟采集一次, 一共 120 组数据。 图 1 为 API 主轴 热漂移测量图, 图 2 为主轴箱温度场的测量图, 图 3 为测量 120 min 主轴径向 (两个方向) 及主轴轴向, 共 3 个方向的热漂移曲线。 图 1 主轴热漂移测量图 图 2 主轴箱温度场测试 2机床与液压第 42 卷 图 3 API 测试过程主轴热漂移误差曲线 图中, x1为 y 轴径向, y1为 x 轴径向, z 为主轴 轴向。 2 2 温度测点优化 最优温度测点满足两个条件(1) 在机床主要 的发热源位置;(2) 温度的改变最能反映出机床热 误差的变化情况[6]。 如果最优温度测点通过经验和试 凑的方法来获取, 效率不高。 此处用相关性分析法来 分析所得数据。 相关性分析是对两个或者多个具备相 关性的变量进行分析, 衡量两个变量间的关联程度。 例如温度 T 与主轴 x 方向误差, 用相关性系数 ρx,T来 表示 T 与 x 相互的密切程度[7]。 计算公式为 ρx, T= cov(x, T) D(x)D(T) (1) 通过 MATLAB 编程求得各温度 T 与主轴热漂移 之间的相关系数。 T1~ T8与主轴热漂移 x ( y1), y (x1), z 3 个方向之间的相关性系数, 见表 1。 表 1 温度和主轴热误差之间的相关性系数 方向T1T2T3T4T5T6T7T8 x0.988 7160.988 2030.968 2900.999 1900.991 6780.988 9920.981 9470.988 109 y0.952 3360.998 5170.996 6560.983 0150.956 5240.950 2690.936 2190.947 210 z0.833 2180.945 3980.960 6550.904 6120.840 3170.831 3430.807 3600.826 681 从上表可以看出, T2、 T3、 T4、 T5四个测点温度 与主轴热漂移误差的相关性最大。 则这四个点为最优 敏感热源点, 也就是最优温度测点。 选出 T2、 T3、 T4、 T5的温度值和 x、 y、 z 3 个方向的热漂移误差 值。 3 基于模拟退火遗传算法优化 BP 神经网络的热 误差建模 3 1 基于 GSA 优化的 BP 神经网络的具体实现 基于 GSA 优化的 BP 网络包括结构确定和 GSA 算法的权值和阈值选定。 算法流程图如图 4 所示。 图 4 算法流程图 3 2 模拟退火遗传 (GSA) 算法设计 网络中有 4 组输入参数 T1, T2, T3, T4, 每组参 数有 120 组数据, 这些参数包括所有的网络权值。 遗 传算法总群规模 480, 进化次数为 1 000, 交叉概率 为 0 5, 变异概率为 0 1, 模拟退火算法初始温度为 1 000 ℃, 退温速率为 0 99。 BP 网络训练参数 最 大步数 1 000, 学习速率 0 01, 目标 0 001[8]。 根据改进后的 Paul L Stoffa 适应度拉伸方法, 个 体适应度 F 为 F = ef i/ t ∑ n i = 1 ef i/ t fi= 1 ∑ g k = 1 (Ls - L r) 2 t = t0(0 99h-1)(2) 式中 fi为改进前第 i 个个体的适应度, g 为输出神 经元个数, Ls为期望输出值, Lr为实际输出值, t0为 初始温度, h 为遗传进化代数[9]。 从而通过选择、 交 叉、 变异来找出最具适应度的个体。 选择概率计算方法作为排序法进行排序。 排序后 个体被选中的概率为 Pi Pi= r(1 - q) (a-1) r = q 1 - (1 - q) n (3) 式中 q 为个体选择概率, n 为种群个数, a 为排序 后个体在种群中的位置。 交叉操作为整数交叉算子法, 在参与交叉的个体 X (x1, x2, , xn) 及 Y (y1, y2, , yn) 中随机 选择一个整数 0<m<n, 即 X (x1, x2, , xm, ym+1, , yn) 和 Y (y1, y2, , ym, xm+1, , xn)。 变异操作方法可表示为 初始个体为 X=x1, x2, 3第 23 期张景然 等 基于模拟退火遗传算法优化 BP 网络的数控机床温度布点优化及热误差建模 , xm, , xn变异后个体是 x2, , ym, , xn, 变异点 Xm取值范围为 [Imin, Imax], 变异后 ym= xm+ f(t, Imax - x m), a = 0 xm- f(t, xm - I min), a = 1 { (4) 其中 f (t, y) = y [c (1-d) ] 3, 式中 t 为当前进化 代数, a 和 c 为 random [0, 1], d 为当前代数与最 大进化代数之比, y 取值为[Imin,xm]和[xm,Imax]。 而后进行退火操作, 其中接收函数为, min{1, e -Δ/ tk} > random[0, 1] tk = μt k-1 { (5) 3 3 BP 网络模型的设计 根据 GSA 算法找到最优群体, 即得到对应的最 优权值和阈值。 文中采用 3 层的 BP 网络进行建模。 4 个关键的温度测点作为 BP 神经网络的输入, 输入 层 4 个节点, 通过定理可计算中间隐层节点数为 9 个。 输出层 1 个节点。 该模型记为 4-9-1。 权值个数 利用公式得出 49+19= 45, 阈值个数 9+1= 10。 编 码后的串长度为 45+10=55。 权值和阈值在[-1,1]范 围[10]。 输入与输出的关系, y = f(∑ m t = 1 ωixi - θ) (6) 其中 xi为输入, y 为输出, ωi为权值, θ 为阈值。 3 4 基于模拟退火遗传算法优化的 BP 神经网络 训练及结果 通过 Matlab 的 BP 神经网络库, 建立 BP 网络, 并根据以上步骤进行优化, 并对 BP 网络进行训练, 达到满足精度要求误差 110 -5 mm 的拟合曲线。 图 58 的实际模型是由最优测点 T2、 T3、 T4、 T5 的温度值为横坐标, 图 9 的实际模型是以主轴温升测 试时间为横坐标, 图 59 都是以 x、 y、 z 3 个方向的 热漂移误差为纵坐标, 拟合模型的纵坐标是经过优化 的 BP 神经网络建立的 x、 y、 z 3 个方向的拟合误差, 图中, 下标 r 代表实测的误差, 下标 s 代表拟合的误 差。 平滑一些的曲线是拟合曲线。 两者比较表明了文 中建模方法的精确性。 相比单一的 BP 神经网络建模, 拟合值更加接近于实际测量值, 残差更小。 图 5 误差随 T2测点温度的变化曲线 图 6 误差随 T3测点温度的变化曲线 图 7 误差随 T4测点温度的变化曲线 图 8 误差随 T5测点温度的变化曲线 图 9 误差随时间的变化曲线 (下转第 50 页) 4机床与液压第 42 卷 图 8 AGC 四缸同步跟随曲线 4 结束语 提出的耦合误差计算方法, 与基于模糊神经网络 和偏差耦合控制结构相结合的控制方案, 研究了对影 响系统协同性能的耦合误差的补偿, 可以有效地补偿 在中厚板矫直机动作中存在的位移协同误差, 提高中 厚板的矫直精度。 经过两种控制方案实验曲线的分析 比较, 得知相对于无偏差耦合控制的四路阀控缸系 统, 采用基于模糊神经网络算法的偏差耦合控制结构 时, 系统的控制特性表现在稳态性能高, 响应快, 抗 干扰性能增强, 实现了很好的控制效果, 验证了基于 模糊神经网络算法的偏差耦合控制结构, 在实际的四 路阀控缸协同控制系统中实现了良好的控制效果。 参考文献 [1] HOLMES D G,LIPO T A.Implementation of a Controlled Rectifier UsingAC⁃AC Matrix Converter Theory[J].IEEE Trans.Power Electron,1992,7(1)240-250. [2] SHIH Y T,CHEN C S,LEE A C.A Novel Cross⁃coupling Control Design for Bi⁃axis Motion[J].Int J of Machine Tool and Manufacture,2002,42(1)1539-1548. [3] 肖本贤,郭福权,王群京,等.基于模糊神经网络的轮廓 误差附加补偿控制研究[J].系统仿真学报,2003,15 (12)1733-1736. [4] 张承慧,石庆升,程金,等.一种多电机同步传动模糊神 经网络控制器的设计[J].控制与决策,2007,22(1)30- 34. [5] 李宏杰,黄庆学,琚长江,等.全液压十一辊矫直机的多 嵌入式协同控制系统[J].塑性工程学报,2013,20(2) 73-77. [6] 邢宗义,张媛,侯远龙,等.电液伺服系统的建模方法研 究与应用[J].系统仿真学报,2009,21(6)1719-1725. [7] 曹春平,王波,胥小勇,等.基于神经网络多电机偏差耦 合同步控制研究[J].控制工程,2013,20(3)415-418. [8] PEREZ⁃PINAL,F J,CALDERON G,ARAUJO⁃VARGAS I. Relative Coupling Strategy[C].IEMDC’03.IEEE Interna⁃ tional,Electric Machines and Drives Conference,2003,2 1162-1166. [9] 黄镇海,王桂荣,刘芳璇.多通道液压伺服系统的同步设 计[J].机床与液压,2012,40(11) 29-31. [10] 顾德英,吴成赛,侯娇.基于补偿模糊神经网络的 BLD⁃ CM 伺服控制[J].东北大学学报自然科学版,2013,34 (1)13-16. (上接第 4 页) 4 结束语 利用相关性分析, 找出与主轴热漂移关联性最大 的最优敏感测点, 极大地减少了温度测点的数量, 避 免了繁琐的测量过程。 通过基于模拟退火遗传算法优 化 BP 神经网络进行热误差建模, 根据测得的温度 值, 建立预测模型。 与传统单一的 BP 神经网络模型 进行比较, 文中建模的残差较小, 表明基于 GSA 优 化的 BP 网络热误差模型的精度更高鲁棒性更强。 参考文献 [1] KONO Daisuke,MATSUBARA Atsushi,YAMAJI Iwao,et al.High⁃precision Machining by Measurement and Competi⁃ tion of Motion Error[J].International Journal of Machine Tools and Manufacturer,2008,48(10)1103-1110. [2] ZHANG Yi,YANG Jianguo. Modeling for Machine Tool Thermal Erro Based on Grey Model Preprocessing Neural Network[J].Journal of Mechanical Engineering,2011,47 (7)134-139. [3] 吕琼帅,王世卿.基于遗传模拟退火算法优化的 BP 神经 网络[J].计算机与现代化,2011,26(6)91-94. [4] 黄洪钟,赵正佳,姚新胜,等.遗传算法原理、实现及其在 机械工程中的应用研究与展望[J].机械设计(设计理论 与方法),2000(3)1-5. [5] 王慧琳,胡树根,王耘.基于模拟退火遗传算法优化的 BP 网络在质量预测中的应用[J].轻工机械,2011,29 (4)26-35. [6] 王秀山,杨建国,余永昌,等.双转台五轴数控机床热误 差建模、检测及补偿实验研究[J].中国机械工程,2009, 20(4)405-408. [7] 张志勇,精通 MATLAB6.5 版[M].北京北京航空航天 大学出版社,2004. [8] 任小洪,徐卫东,刘立新,等.基于遗传算法优化 BP 神经 网络的数控机床热误差补偿[J].制造业自动化,2011, 33(5)41-43. [9] 刘月娥,何东健,李峥嵘.一种用于 BP 网络优化的并行 模拟退火遗传算法[J]. 计算机应用,2006,26(1) 204-206. [10] 侯媛彬,杜京义,汪梅.神经网络[M].西安西安电子 科技大学出版社,2007. 05机床与液压第 42 卷
展开阅读全文