基于区间理论的精密卧式加工中心几何精度优化分配方法与流程

本发明涉及数控机床几何精度优化分配领域，尤其是一种基于区间理论的精密卧式加工中心几何精度优化分配方法。
背景技术：
：精密数控机床是基础制造装备的工作母机，其设计和加工能力直接影响制造业的发展。高档数控机床质量的好坏，直接表现在机床的制造精度和装配精度，而机床精度设计和精度分配的优劣，对其质量好坏有着重大影响。对于机床精度分配问题，国内外学者进行了较为广泛和深入的探究，取得了一定的进展。主要有零部件公差加权欧式范数最大法、蒙特卡罗法、尺寸链法、基于公差—成本模型的精度分配方法等方法，这些方法大多针对并联机构，并且无法直接指导高档数控机床的制造装配过程。目前，国内的机床企业在进行精密数控机床精度设计和精度分配时，通常采用传统的类比、机床手册查询以及设计人员估算的方法，没有从根本上解决问题，设计的机床难免会过分依赖设计人员的经验，无法形成科学化、理论化、系统化的设计方法。精密数控机床待切削点到刀尖点的空间误差直接反映了机床的加工精度，而空间误差是通过各运动副几何误差耦合作用产生的。因此，如何在满足整机空间误差的前提下，合理、经济以及快速地优化分配几何误差是目前精密数控机床装配制造亟需解决的问题。技术实现要素：本发明的目的在于解决上述技术问题而提供一种基于区间理论的精密卧式加工中心几何精度优化分配方法，其能够使难于实现的误差源放松至最大的可行区间。为了解决上述技术问题，本发明采用如下技术方案：一种基于区间理论的精密卧式加工中心几何精度优化分配方法，包括以下步骤：步骤(1)：首先将数控机床看作由多个刚体串联而成的多体系统；其次使用多体系统运动学理论描述数控机床上各运动部件之间的拓扑关系；再次使用齐次坐标变换矩阵相乘表示数控机床上相邻部件连体坐标系之间的位姿转换关系；最后建立起数控机床各几何误差与机床空间误差之间的映射模型；步骤(2)：基于区间理论，根据区间扩张因子定义灵敏度，并进行灵敏度分析以反映各误差源对整机空间误差的影响程度；步骤(3)：以各误差源区间宽度为设计变量、整机空间体对角线误差为约束条件、各误差源灵敏度系数与区间宽度负二次方乘积的代数和为目标函数，建立区间优化分配模型；步骤(4)：基于遗传算法，运用matlab软件来求解从而实现几何精度优化分配。本发明的有益效果是：本发明针对精密卧式加工中心提出了基于区间理论的几何精度优化分配方法，基于多体系统运动学理论建立机床几何误差的整机空间误差映射模型，基于区间理论进行灵敏度分析并建立区间优化分配模型，最后基于遗传算法运用matlab软件来求解从而实现几何精度优化分配，使难于实现的误差源放松至最大的可行区间，大大降低了生产成本，提高了机床生产装配效率。附图说明图1为本发明基于区间理论的精密卧式加工中心几何精度优化分配方法的流程图；图2为u2000/800h精密卧式加工中心三维模型图；图3为遗传算法原理图；图4为各几何误差源精度优化分配结果。其中，1-床身2-x向滑台3-b向转台4-工件5-立柱6-主轴箱7-主轴8-刀具。具体实施方式下面结合附图和实施例对本发明作进一步详细的说明：如图1所示，本发明提供的基于区间理论的精密卧式加工中心几何精度优化分配方法，包括以下步骤：(1)建立机床各几何误差和机床空间误差之间的映射模型；数控机床的几何误差，根据是否与机床运动位置相关，可分为位置相关几何误差(pdges)和位置独立几何误差(piges)。本发明以u2000/800h精密卧式加工中心为例，模型如图2所示，对本发明的方法加以说明，上述加工中心在仅考虑平动轴几何误差的情况下共有21项几何误差：与三个轴线坐标位置相关的9项位置误差和9项转角误差，以及位置独立的三个轴线之间的3项垂直度误差，共计18项位置相关几何误差、3项位置独立几何误差，如表1所示。表1数控机床的21项几何误差坐标轴pdgespigesx轴δx(x),δy(x),δz(x),εx(x),εy(x),εz(x)--y轴δx(y),δy(y),δz(y),εx(y),εy(y),εz(y)εxyz轴δx(z),δy(z),δz(z),εx(z),εy(z),εz(z)εxz,εyz表中δx、δy、δz代表位置误差，下角标表示误差的方向；εx、εy、εz代表转角误差，下角标表示转角误差旋转轴线的方向；括号中的x、y、z表示x、y、z轴的运动方向。εxy、εyz、εzx分别表示坐标轴x与y、y与z、x与z之间的垂直度误差。“--”表示没有误差。数控机床是一种由多个刚体串联而成的多体系统，可使用多体系统理论来描述数控机床上各运动部件之间的拓扑关系，用齐次坐标变换矩阵相乘表示数控机床上相邻部件连体坐标系之间的位姿转换关系，从而建立起数控机床21项几何误差与机床末端位姿误差之间的映射模型(以下简称机床误差模型)。根据图2的u2000/800h卧式加工中心的拓扑结构得到相邻体的理想静止、运动特征矩阵以及理想静止、运动误差特征矩阵，如表2所示。表2相邻体的理想静止、运动矩阵以及理想静止、运动误差矩阵表中tijp、δtijp、tijs、δtijs分别表示相邻体的理想静止特征矩阵、理想运动特征矩阵、理想静止误差特征矩阵、理想运动误差特征矩阵；x、y、z分别表示机床运动坐标；xijp、yijp、zijp分别表示机床结构参数；i4×4表示四阶单位矩阵；i,j＝1,2,3,…,8；1表示床身，2表示x向滑台，3表示b向转台，4表示工件，5表示立柱，6表示主轴箱，7表示主轴，8表示刀具。设刀尖点在主轴坐标系下的齐次坐标值pt为：pt＝[00t1]t式中t表示刀具长度。设工件上待加工点在工件坐标系下的齐次坐标值pw为：pw＝[xwywzw1]txw,yw,zw表示工件上待加工点在工件坐标系下的坐标值。刀尖点在理想床身坐标系下的齐次坐标值pt-ideal为：pt-ideal＝t04pt04st45pt45st56pt67ppt工件上待加工点在理想床身坐标系下的齐次坐标值pw-ideal为：pw-ideal＝t01pt01st12pt12st23ppw在理想条件下，刀尖点和工件上待加工点是重合的，则有：pt-ideal＝pw-ideal刀尖点在实际床身坐标系下的齐次坐标值pt-actual为：pt-actual＝t04pδt04pt04sδt04st45pδt45pt45sδt45st56pδt56pt67pδt67ppt工件上待加工点在理想床身坐标系下的齐次坐标值pw-actual为：pw-actual＝t01pδt01pt01sδt01st12pδt12pt12sδt12st23pδt23ppw空间误差为：式中：a、分别为几何误差映射矩阵和几何误差向量。将上述机床几何误差模型按位置相关的几何误差和位置独立的几何误差写成如下形式：式中表示空间误差；ad、ai分别为位置相关的几何误差映射矩阵和位置独立的几何误差映射矩阵；εd、εi分别为位置相关的几何误差向量和位置独立的几何误差向量。(2)基于区间理论，进行对各误差源进行灵敏度分析；区间数定义为：式中上角标i表示区间数，后面出现的上角标i含义相同；分别表示区间的上下限。将几何误差空间模型区间化：式中上角标t表示转置。区间扩张因子ηi,j定义为：式中λ用于描述区间扩张因子的正负；分别表示第i个空间误差分量的区间上下限；表示第j个几何误差向量分量的区间上下限；ai,j表示几何误差映射矩阵的第i行，第j列的元素。灵敏度定义为：对于位置误差源，灵敏度定义：对于角度误差源，由于阿贝误差的影响，灵敏度定义：式中表示机床位于坐标且刀具长度为时的灵敏度；表示机床位于坐标且刀具长度为时的区间扩张因子；表示机床位于坐标原点且刀具长度为零时的区间扩张因子。(3)基于区间理论，建立区间优化分配模型：以各误差源区间宽度为设计变量、整机空间体对角线误差以及各误差源边界为约束条件、各误差源灵敏度系数与区间宽度负二次方乘积的代数和为目标函数，建立区间优化分配模型。式中min表示求最小值；f表示目标函数；xi表示区间变量的集合；n表示区间变量的总数；μi表示第i个区间变量的灵敏度；w表示区间宽度；表示第i个区间变量；δx,δy,δz分别表示x，y，z方向的空间误差；[xi]表示区间变量的约束区间。(4)基于遗传算法(ga)，运用matlab软件的遗传算法工具箱来求解从而实现几何精度优化分配；基于遗传算法，原理图如图3所示，运用r2014a版本的matlab软件的遗传算法工具箱，遗传参数采用默认值，进行仿真求解，运行结果如图4所示，各误差源几何精度优化分配结果如表3所示。表3各误差源几何精度优化分配结果误差源δx(x)δx(x)δx(x)δx(x)δx(x)δx(x)δx(x)δx(x)δx(x)经验值/um655655655优化值/urad8.97.47.48.97.47.48.97.47.4放宽/％48.3484848.3484848.34848误差源εx(x)εx(x)εx(x)εx(x)εx(x)εx(x)εx(x)εx(x)εx(x)经验值/um888888888优化值/urad11.51211.711.910.79.511.21111.9放宽/％43.85046.348.833.818.84037.548.8误差源εxyεxyεxy经验值/um101010优化值/urad14.911.112.8放宽/％491128尽管上面结合附图对本发明的优选实施例进行了描述，但是本发明并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，并不是限制性的，本领域的普通技术人员在本发明的启示下，在不脱离本发明宗旨和权利要求所保护的范围情况下，还可以做出很多形式，这些均属于本发明的保护范围之内。当前第1页12