基于神经网络的薄壁件加工变形预测方法

本资料包含pdf文件1个，下载需要1积分

减杏工误差、提高零件的表面质量是制造企业追求的 目标 。在航空、航天、兵器等国防工业领域的高性能零件广泛采用复杂薄壁结构，其精度要求越来越高，如运用在雷达上的军品铝质鼻件和薄筒件 、飞机上的翼肋和发动机上的涡轮 叶片等薄壁件，都具有尺寸大、刚度弱、表面精度高等特点，加工中采用数控铣削或数控车削成形。在这些结构件的加工中，由于其弱刚度、残余应力、过大切削力、切削热、不合理的刀具结构与切削参数等因素，引起工件产生加工变形 ，极端情况甚至报废 。为了减小工件的变形，需对薄壁件切削加工进行深入研究。

有限元模拟能够揭示实验方法难 以获得的切削力、应力、变形、材料失效应力等加工数据，目前许多学者利用有限元分析方法对薄壁件加工变形分析进行了大量研究。张臣等 视参与切削的切削单元受力之和为实际切削力，建立球头铣削过程中刀具的受力变形模 型，研究 了球头铣刀刀具 变形 引起 的加工误差仿真分析方法。梅中义等 研究了残余应力对加工变形的影响，采用单元生死方法实现切削过程的仿真。然而，这些研究工作仅限于工况条件已知的工件变形有限元分析 ，不能实 现工况未知条件下工件变形的预测。

在假定铣削速度、铣削深度、进给量和铣削宽度引起的工件变形量已知条件下，唐东红等 建立了铣削速度、铣削深度、进给量与工件变形之间的神经网络模型，通过均匀实验设计方法获得训练样本，最终利用 Levenberg-Marquardt算法实现对工件变形的预测。李目 在建立 BP神经网络非线性映射模型时，将铣削的主轴转速和每齿进给量作为 BP神经网络输入特征预测工件加工变形。而刘新玲等根据切削刃对工件作用时间的先后对已加工表面上相应节点进行加载、卸载，利用载荷步方法对动态切削力引起的工件变形进行有限元分析，通过 BP神经网络算法建立了摆线轮铣削力与加工变形的预测模型。上述工件变形的预测方法仅考虑了铣削工艺参数的影响，而忽略了刀具的几何参数。

薄壁件切削加工过程 中引起加工变形的影响因素多种多样，多个影响因素的组合则更是成千上万。

显然，仅仅根据加工变形的影响因素利用有限元方法分析计算工件的加工变形，是不切实际的。为此，本文针对薄壁件铣削过程建立了加工变形的有限元仿真模型，分析了刀具几何参数对加工变形的影响规律 ，并利用实验数据验证 了有 限元分析模 型的有效性¤助 BP神经网络的非线性映射能力，建立了铣刀前角、螺旋角与工件加工变形的非线性映射模型 ，从而实现了利用有限的变形分析结果预测任意前角 、螺旋角引起的加工变形 。该方法不仅为实际加工实现智能化提供条件，还为薄壁件切削参数优化和路径补偿提供有效数据 。

1 薄壁件加工变形的有限元仿真技术加工过程仿真技术是-项基于工程力学与有限元方法的虚拟制造技术。本文通过切削过程的合理建模、有限元分析计算及数值结果的静动态可视化后置处理，无须对工件进行试切即可获得切削温度 、切削力和切削应力 的分布形式以及工件和刀具 的变形状况的相关数据 。

1.1 有限元分析模型如图 1所示为薄壁件的铣削加工示意图，薄壁部分尺寸分别为：长 z50 mm，高 h25 mm，壁厚t：4 mm，工件材料为铝合金 7050-T7451.加工过程采用螺旋角 卢30。、前 角 。20。的双齿 圆柱立铣刀，刀具直径 d10 mm，刀具材料为硬质合金YG8.加工参数为径 向铣削宽度 o 2 mm，轴 向铣削深度。 25 mm，每齿进给量 0.125 mm，主轴转速n1 200 r/min，摩擦系数 0.3。

图 1 薄壁件铣削仿真示意图Fig.1 The simulation diagram of the thin-walledworkpiece milling1.2 工件材料的本构模型1.2.1 材料的本构关系Johnson.Cook准则利用变量乘积关系分别描述应变 、应变率和温度的影响 ，常用于模拟金属材料从低应变率到高应变率下的动态行为。Johnson-Cook842 兵 工 学 报 第34卷准则具体表述为 ( n -( ，(1)式中： 为累积塑性应变；s 为累积塑性应变率；s。为基准塑性应变率；0为温度；0o为基准温度；0 为融化温度 ；A、B、C、m和 n为材料属性。

1.2.2 材料的失效准则随着刀具的切入，当工件材料的等效塑性应变达到 1.0时，材料失效，实现被切单元与工件脱离。其失效准则是把每个工件单元定义为等效塑性应变，即D：∑ ， (2)式中： 为 Johnson-Cook通过大量的实验总结的材料失效准则 ，考虑了应变、应变率、温度和压力参数。Johnson-Cook失效准则提供了材料到达失效点时等效塑性应变的计算方法，具体形式为： [D。D2exp(D。 )]f 1D41 芒1。

0 D5( ) ， ㈩式中：屈服应力的变化率 or： ，O'y为屈服应力，为非应变率下的屈服应力；D (i1，2，3，4，5)为材料参数。

1.3 材料模型参数铝合金 7050-T7451的密度为2 820 kg/mm ，弹性模量为 71.7 GPa，泊松 比为 0.33.另外，7050-T7451的热导率、线热膨胀系数以及比热容参数如表 1～表3所示 。

表 1 热导率Tab.1 Thermal conductivity表 2 线膨胀系数Tab.2 Thermal expansion由于Johnson.Cook材料本构模型引入了应变率表 3 比热容Tab.3 Specific heat强化、应变硬化以及热软化参数，综合反映了金属材料大应变、高应变率和高温加载下的本构关系，对各种材料参数均适用且形式简单，所以本文选用该模型。7050-T7451的 Johnson.Cook材料本构模 型参数如表4所示 71。Johnson-Cook材料失效准则参数如表 5所示 。

表4 材料模型参数Tab.4 Constants of the Johnson.Cook constitutive modelDl D2 D3 D4 D5O所用刀具为硬质合金 YG8，其化学成分为 92%的 WC和 8% 的 Co，密 度 为 14.4～14.6 kg/m ，YG8的耐磨性 良好，使用强度和冲击韧性非常优秀。适用于铸铁、有色金属及其合金、非金属材料的粗加工，在铝合金材料的铣削过程中变形量很校YG8的主要物理与热力学性能如表6所示 。

表 6 刀具 的物理 与热力学性能Tab.6 Material parameters of the cutting tool 。

杨氏模 泊松 线膨胀系数/ 比热容/ 热导率/材料量 比 10-6C/MPa - (J.kg-.℃- )(W.m～.K-，) · 。· - ·1.4 仿真结果分析根据有限元仿真计算结果，获得工件的变形云图，如图 2所示 。

实验是在 Micron UCP710五轴加工中心上进行的，铣削方式为单向走刀、顺铣，铣削参数与仿真过程相同 。由于工件在 方向上变形最大，故于铣削完成后，在已加工表面上取 Y10 mm，25 mm，40 mm 3条垂直线，分别沿 轴方向任意取 0.5 mm，2.5 mm，5 mm，7.5 mm，10 mm，12.5 mm，15 mm，18 mm，22 mil共 9个 点。利用 Global Classic 9128型三坐标测量机分别对 9个点加工前与加工后的第 7期 基于神经网络的薄壁件加工变形预测方法图 2 薄壁件的变形云图Fig.2 Deformation nephogram of the thin-waled workpiece轴坐标值进行测量，其差值即为加工变形的大小，其变形 曲线如 图 3所示 。另外 ，在加工表面 。8 mm处，沿 Y方 向分别取 Y1 mm，5 mm，10 mm，15 mm，20 mm，25 mm，30 mm，35 mm，40 mm ，45 mm，49 mm共 11处工件变形数据，绘制出工件的加工变形状态，如图4所示▲行对比分析后可知 ，无论是仿真结果 的大小 ，还是变形趋势，都与实验值具有较好的吻合度 ，从而验证 了图 1中薄壁件有限元分析模型的正确性。

2 薄壁件加工变形的神经网络预测方法由第 1节内容可知，利用有限元方法可以分析给定刀具前角与螺旋角情况下薄壁件的变形状态。

然而 ，刀具前角 、螺旋角 的组合情况多种多样 ，不可能利用有限元方法--计算与分析，否则费时费力也找不出刀具结构与加工变形之间的复杂关系。因此，本节将建立基于 BP神经网络的薄壁件加工变形预测方法，实现对刀具角度的有效性评估与修改，最终保证零件的加工质量，缩短制造周期，降低制造成本 。

2.1 BP网络结构的确定BP神经网络是-种具有 3层或 3层 以上 的前馈神经网络，包括输入层、隐藏层和输出层，可任意的精度逼近任意的连续函数 。Hecht-Nielsen证明具有 1个隐藏层的3层前馈型网络可以逼近任何多变量函数，故本文采用1个3层 BP神经网络进行铣削变形的预测。BP网络结构如图 5所示 ，输入层有刀具前角和螺旋角2个神经元，输出层有 1个神经元 ，即需预测的由刀具结构产生 的最大变形量 ，再根据 Kolmogorov定理，并综合考虑网络的学习速度和吕吕血i制图3 垂直线的加工变形模拟值与实验结果的比较Fig.3 Comparison of the simulated values of the vertical linedeformations with the experimental results图4 水平线的加工变形模拟值与实验结果的比较Fig.4 Comparison of the simulated values of the horizontalline deformations with the experimental results泛化能力，最终确定隐含层包含4个神经元。

2.2 训练样本的选择根据训练样本的相容性、遍历性和致密性 3要素要求，螺旋角的采样间隔 15。，分别取 30。，45。，兵 工 学 报 第34卷输入层 隐藏层 输出层量图 5 BP网络结构 图Fig.5 Configuration of BP neural network60。，而前角的采样间隔5。，分别取 -10。，-5。，0。，5。，1O。，15。，20。，25。，30。，最后 以 l4种工况下铣刀引起的工件最 大变形作为 网络 的训练样本 。利用第 1节已经验证过的有限元分析方法，容易获得工件的最大变形，如表 7所示。

表 7 不同铣刀角度下的切削变形Tab.7 Workpiece deformations under diferent angleof miling cuter2.3 神经网络的训练由于基于高斯-牛顿法改进得到的 Levenberg-Marquardt算法在实际应用中具有很快的收敛速度，故选择 Levenberg-Marquardt算法作为前馈型神经网络的学习算法对权系数进行更新和计算，实现输入与输出间的非线性映射。由于系统是非线性的，初始权值对于学习能否达到局部最型是否能够收敛的影响很大，这就要求初始权值在输入累加时使每个神经元 的状态值接近于 0，输入样本要进行归-化处理，使那些比较大的样本数据输入仍落在传递函数梯度大的地方，可以采用(4)式进行归-化处理，使各个样本数据都落在[0，1]。

互 ： ， (4) - - 斗，㈨- min式中： (i1，2，，14)为第 个输入样本； 与"m ax分别为样本最小值和最大值。

样本归-化 后，就 可进行 网络训 练。根据 BP神经网络的-般设计原则及大量实验，训练中隐藏层神经元的传递函数为双曲正切函数，实现不同样本的空间分割，而输出层神经元的传递函数为线性函数，输出网络的识别结果。初始化网络权系数在[0，1]之间，学习误差为 0.000 001.将表 7中的14个样本输入网络进行训 练，网络在 63步收敛到要求精度，神经网络训练过程误差曲线如图6所示。

图6 网络训练结果Fig.6 Results of network training利用仿真函数 sim对上述网络进行预测 ，并将预测值与仿真值对比，如表 8所示。由此可见，测试样本的网络输出值与有限元计算结果之间的误差不超过 3%。

3 结论1)本文以薄壁件铣削过程有限元仿真模型为基础，利用 BP神经网络卓越的学习能力，有效地建立了铣刀角度与加工变形之间的非线性逻辑关系。通过表8数据可以看出，该网络结构预测结果误差合理。

2)BP神经网络具有较高的精度和良好的泛化能力，只要选择合适的网络结构、传递函数和训练参数，合理地设计训练样本，神经网络模型就可以有效地用来预测薄壁件的变形量。这些利用 BP神经网络预测的数据还可以作为样本数据，采用遗传算法第 7期 基于神经网络的薄壁件加工变形预测方法 845表 8 样本测试数据输出值与有限元计算结果比较Tab.8 Comparison of the predicted results with thesimulated values对铣刀的前角与螺旋角进行优化，有效地实现薄壁件铣削变形的控制。