特殊工件精密磨削加工方案的一种规划算法

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

A Programming Algorithm of Special PartsPrecision Grinding ProjectLI Sen(Basis Department，Xiangyang NCO School，Military Economy Academy，Xiangyang Hubei 441 l 18，China)Abstract：When using the grinding machine to do the precision grinding work of massive special parts which are solids of revolu-tion with complex generator line，the precision and eficiency of grinding must be ensured.On the basis of having analyzed the motionlocus of the grinding point，a mathematical model based on the algorithm of dynamic programming and nonlinear integer programmingwas built，and the solving algorithm was provided. An explicit problem is solved，and the error analysis is done。

Keywords：Grinding；Program ming；Algorithm在工程上，加工大批量具有复杂母线旋转体的特殊工件时，最后的精密成形工艺采用磨削加工，此时，加工质量的高低以及加工时间的长短，是非常重要的▲行精密磨削加工的主要设备有精密数控磨削机床、各种研磨机等 ，作者以某三轴联动精密数控磨床为例，研究如何确定-个合理的加工方案，使得在加工时间旧能短的前提下 ，加工质量旧能的高。

该磨床的结构是：工件工作箱固装在上台面上，它通过专用夹具装夹工件，使工件绕工件工作箱主轴以较慢的转速旋转，同时随上中下3个工作台的复合运动改变待加工工件与砂轮的相对位置。3个工作台的运动必须相互配合，使工件与砂轮相切磨削，加工出满足要求的旋转体，如图1所示p 。

图1 磨床结构示意图座该磨床的主要技术参数是 ：各组步进电机的步进角度均为1。，步进电机控制脉冲的最高工作频率不大于 100脉冲/s，各组变速器的传动比为 10：1，丝杆的螺距为 12 mm，中台转轴到上工作台的控制丝杠 -螺母副中心线的距离R300 mm，工件工作箱的夹具基准面到中台转轴的距离 b250 mm，如图2所示。同时，为保 证表面的加工质量，要求 ： (1)l f-f l 6(、 厂h分别指相邻两时间段中前后时段l a l的脉冲频率，-般取 ≤0.10)；(2)工件每转动100转，工件与砂轮的磨削点在工件工作箱的旋转轴方向上的移动量不超过4 mm；(3)工件工作箱主轴转动速度为 250～300 r/min。

图2 上台相对中台的转动关系示意收稿日期：2012-03-02作者简介：黎森 (1983-)，男，硕士，讲师，研究方向为计算数学。E-mail：greatlisen###gmail.com。

· 6O· 机床与液压 第4l卷1 坐标系的建立及磨削点的运动轨迹为了寻找磨削点的运动规律，以夹具基准面中心为坐标原点，工件转动轴为X轴，在水平面建立如图1和图2所示的平面直角坐标系，设需要的工件外表面母线方程为Y.厂( )(0≤ ≤a)。

在该坐标系下，无论3个工作台如何运动，母线方程始终是保持不变的，变化的是砂轮的位置以及磨削点的坐标。由于砂轮的厚度远小于工件母线的长度，因此可以近似地认为各个瞬时砂轮与工件的磨削点是在砂轮的同-个位置上。在这种假设下，磨削点坐标的变化仅仅由3个工作台的平动与转动来决定。

假定脉冲个数可正可负，规定向3组步进电机各输入1个控制脉冲时，控制下、中、上台运动的螺母分别向左、上、上方平动，从而使下、中台分别向左、上方平动，上台逆时针转动；各输入 -1个控制脉冲( Y-i )COSi； 。sin A△/卢3i /，·( Y-1i--1 )时，使下、中台分别向右、下方平动，上台顺时针转动。在第i个很短的时间段内，设控制下、中、上台的步进电机分别发出n。 、n2i和n， 的脉冲，此时上台会转动△ 的角度 (△ >0时为逆时针转动，△ <0时为顺时针转动)，使得中台到上台的角为 。这时I n l、I n：；l和I n l会很小，由它们造成的3个丝杆的移动距离以及I I会更加小，因此可以近似地认为上台的转动和下、中台的平动是分开进行的，先转动后平动，其中上台是以点 (b，0)为中心进行转动的。

i il2∑n3f ∑n3f显然，卢 ：arctan arctan 30 0R，令3o0，则 △ 13 - -。。

设( ，Y )和( ，Y )分别是第 i个时间段始、末时磨削点的坐标，根据坐标轴旋转与平移公式 ，容易得出：12(nilcos3 -n21si f)10×36012(nilsi n21cos3f)- - - - - ：( -。-6)c。s(卢 -卢 -。)y -。sin(/3，-3 - )6-nu-C-O-S- ，6-石n-2-isinfl即：iY ：- - -b sin(fl -卢 -Y -cos(3 -卢 -- 旦兰 lf(1) 1)flf) Juv开始加工前，将砂轮的初始位置固定于待加工工件计划母线的最左侧，也即， 。0，Y0，(%)。这样，就可以计算出磨削点的运动轨迹。

2 基于规划算法的JnT方案模型为了便于加工方案的实施以及模型的建立，不妨设所有时间段内工件在 轴方向上的移动距离均为， 设法求解每个时间段内的加工方案。设第 i个时间段的时间长度为t ，在该段内工件工作箱主轴转动速 度 为 ri，则 根 据 要 求 可 得：-100 x-6-0Ax ≤ 4，/；i 。 ，250 r/min≤r <300 r/min。为了使得加工时间旧能的短，可以让前面第-个不等式的等号成立，并且让所有的 rl300 r/min，这时，所有的最小 t 5hx。

为使得加工旧能地精细，可将每个时间段内加工方案的目标设置为局部误差最小，这也就巧妙地将-个较为复杂的整体规划问题转化为若干个相对较易解决的子问题，是-个动态规划模型。

以下介绍每-个子问题的求解过程。为保证加工精度，设定-个可以接受的局部误差范围占，即在第i个时间段内，lY -f( )l≤ 。 的选取必须合适，如选得太大，加工质量不够精细，选得太小，则话可能没有可行解。此时，第 1个子问题 的规划模型如下 ：min lY1-f( 1)l.6c。 6)，。6si 。)c。 --nns-infll n-21cosflln剞卢t arc协 ly1-f( 1)l≤Inll I≤100， In21 l≤1O0，In31[≤1O0(1)当i≥2时，第 i个子问题必须加上相对于上-时间段脉冲频率变化范围的约束，其规划模型如下：minIYf-f( )I第7期 黎森：特殊工件精密磨削加工方案的-种规划算法 ·6l·( -。-b)cos( -卢 - )y -。sin(B，-B - )6-n-l-c-os-fl l- n-2i-sinfli)， -( - -b)sin(卢 -卢 - )y -，c。s(卢 -卢 -。)--nl-is-in-fl l n-2i-cosfll圭蒜 f2)IY -/ f)I≤0， 0， 0I 二 旦J≤6，I I≤6，I 二 翌I≤n1. -1 n 2，f-1 n3，i-lnn， 和 均 为整数可以通过观察得出，每个子问题模型又是-个非线性整数规划模型，并且当i≥2时都需要调用上-个子问题模型计算得到的结果 ，子问题的数量撒于Ax的选取 ，如 Ax太大，就不能近似认为上台的转动和下、中台的平动是分开进行的，从而使实际工件与计划工件之间的误差增大，太小则会导致子问题数量太多，计算任务太重。当求解某个子问题模型出现无可行解的情况时，可以根据求解情况在不太影响加工质量的前提下尝试适当放宽-些约束条件，例如，适当扩大6或 的取值范围，重新进行求解，计算出新规划的最优解后再进行下-子问题模型的求解∩令 ， 则子问题的个数为17，个，第 个时间段结束n. - 时磨削点的横坐标 ，每个时间段的时间长度，

，3 模型算法模型算法的步骤如下 ：步骤 1：选认适的/"t、占和6(6≤0.1)，取 ·n 50， ；步骤2：令 i1，求解规划 (1)，输出n 、n：n31、卢1和 Y1；步骤 3：令 i1，判断i是否大于 n，否则继续下-步，是则结束计算；步骤 4：调用 nlIi-t、n2. 、n3.H 、 -1和 Y ，求解非线性整数规划 (2)。若规划 (2)有解，继续下-步，否则根据求解情况尝试适当放宽6或 的取值范围重新进行求解，直到有解才能继续下-步；步骤5：输出 n 、n 、n孙 卢 和Y ，返回步骤 3。

4 误差分析除开操控失误等人为因素，影响工件磨削质量，产生加工误差的因素主要有以下三类：第-类是实际生产过程中-切不可抗拒因素造成的误差，例如磨床的运动精度误差，砂轮的尺寸误差，磨床、砂轮及工件的热变形和弹性变形误差 ，还有编程中的计算误差等 ，这类误差另有模型可以辅助减少，例如，文献[6]研究的就是热变形和弹性变形误差的补偿；第二类是近似计算造成的误差，主要包含两方面，分别是忽略砂轮厚度而近似地认为磨削点始终位于砂轮的同-位置以及近似地认为在较短的时间段内3个工作台的平动与转动是分开进行的，这类误差非常小，可以忽略，如果确实要减小此类误差 ，可以换用其他更复杂、计算量更大的模型，例如，文献 [7]和 [8]分别研究的是磨削点在砂轮上的均匀移动问题以及更精确的磨削点运动轨迹问题；第三类误差就是由模型通过仿真计算得到的仿真工件与计划工件之间的误差，下面介绍该误差的计算过程。

首先分析母线误差。由于t 很小，并且在第 i个时间段内，n 、n 和 n， 是均匀分布在其间的，因此可以近似地认为，得到的仿真母线从( ，Y )到( ，Y )都是线段，将所有的( ，Y )( 0，1，2，，n)依次连接 ，就能够得到仿真母线。过 ( ，Y )和( ，Y )两点的直线方程为Y ( - )i - i-1Y ，设g( )为定义在[0，a]上的-个分段函数，对于每-个区间段[ ， ](i1，2，，n)，都有I ， ， ， ，lg( )I，( )-l ( - )yi-。l l。由 此， l L-。

- 1 l仿真母线与计划母线之间的最大误差 e 等于函数g( )在其定义域上 的最 大值，平均误差 e平 rg( )- - 。

下面分析体积误差，其应等于仿真母线与计划母· 62· 机床与液压 第 4l卷线之间围成的区域绕 轴旋转-周得到的旋转体体积，即体积误差：e体.竹l 。 ㈤5 算例分析算例 ：加工外表面母线为 Y - (600-) 0.45(600- )， ∈[0，600]的某旋转体工件，取 n3 000， 0.O1，60.1，使用 Matlab求解模型 ，得到最优加工时间为 3 000 S，加工方案如表1，各工作台累计脉冲数关于时间的图像如图3，仿真母线与计划母线之间的误差分布如图4，最大误差为0.O11 mm，平均误差为8.91×10-mm，体积误差为297.67 toni ，0≤ ≤20时，仿真母线与计划母线的图像如图5。

表 1 加工方案l816l4l2× lO籁 8薹4雠 20.2l80160l4O120籁 100曝 806040200图 3 各工作台累计脉冲数关于时间的图像- 0.0l -0.005 0 0.00$ 0.0l误差/mm图4 仿真母线与计划母线之间的误差分布图5 0≤ ≤20时，仿真母线与计划母线的图像6 结束语对采用某三轴联动精密数控磨床磨削加工具有复杂母线旋转体的特殊工件进行了研究。在分析了磨削点几何运动轨迹的基础上，建立了-个加工质量较高、加工时间最少的动态规划模型，该模型中每个子问题又是-个非线性整数规划模型，给出了求解算法并进行了算例分析。该算法不仅便于计算，而且得出的加工方案非常便于实际操作，在最少的时间内生产出的工件与计划工件之间的误差仪常的小 ，因此，具有高实用性、高精度和高效率的特点。所提供的模型及算法还可适用于其他类型的很多磨床或机床以及加工其它类型的很多精密旋转体工件，具有-定的推广性。