一种考虑几何误差的装配位姿高精度线性化求解方法技术

本发明专利技术公开一种考虑几何误差的装配位姿高精度线性化求解方法，步骤一、测量零件上待装配面的形貌并进行参数化处理，得到每个待装配面面形的数学描述；步骤二、根据待装配面的空间平面关系、形貌误差和受力情况得到待装配面装配位姿的最优化表达式；步骤三、求解所述最优化表达式，得到待装配面的最优装配位姿，能够在装配和测量过程中能够快速地、准确地、定量化地求解零部件之间触位姿关系，提高零部件位姿预测与控制的精度与效率。

A High Precision Linearization Method for Assembly Posture Considering Geometric Errors

【技术实现步骤摘要】
一种考虑几何误差的装配位姿高精度线性化求解方法
本专利技术属于数字化装配领域，具体涉及一种考虑几何误差的装配位姿高精度线性化求解方法。
技术介绍
在机械制造过程中，经常会因装配误差问题导致废品率高、生产周期长、生产成本提高。因此在装配前开展装配工艺模拟与工艺参数优化的研究具有重要的意义。因为零部件之间接触位姿关系的计算与调整优化是装配工艺模拟与工艺参数优化的最为基础的环节，也是影响装配质量的决定性因素，其中零部件之间接触位姿关系是调整优化的基础，所以能够准确快速求解零部件之间接触位姿关系的方法便是关键所在。但是现有的接触位姿求解方法未考虑接触面实际形貌误差影响，求解方法存在缺陷。
技术实现思路
有鉴于此，本专利技术提供了一种考虑几何误差的装配位姿高精度线性化求解方法，能够在装配和测量过程中能够快速地、准确地、定量化地求解零部件之间触位姿关系，提高零部件位姿预测与控制的精度与效率。实现本专利技术的技术方案如下：一种考虑几何误差的装配位姿高精度线性化求解方法，包括以下步骤：步骤一、测量零件上待装配面的形貌并进行参数化处理，得到每个待装配面面形的数学描述；步骤二、根据待装配面的空间平面关系、形貌误差和受力情况得到待装配面装配位姿的最优化表达式；步骤三、求解所述最优化表达式，得到待装配面的最优装配位姿。进一步地，记待装配的两个面为f1和f2，其中，f1为运动面，f2不动面，所述最优化表达式为：其中，ρi为第i点的分布力，θx和θy分别表示x轴和y轴方向的转动量；dz表示z轴方向的平移量；(xi,1,yi,1,zi,1)是待装配面f1上点集,(xi,2,yi,2,zi,2)是待装配面f2上的点集。进一步地，步骤三的求解过程中，将实时更新的xi,1，yi,1和zi,1带入最优化表达式中，求解该优化问题。进一步地，迭代计算的过程如下：1、将待装配的两个面f1和f2参数化；2、建立空间运动矩阵H(θx,θy,dz)；设定装配误差δ，装配误差为装配空间运动后转动角度θx,θy绝对值之和，如果δ无限小则装配稳定，装配完成；设定Sθx＝0,Sθy＝0,Sdz＝0，Sθx,Sθy分别为计算过程中x和y轴方向的旋转累积量，Sdz为计算过程中z轴方向的移动累积量；3、利用式(11)求解θx,θy,dz；4、令Sθx:＝Sθx+θx,Sθy:＝Sθy+θy,Sdz:＝Sdz+dz，更新保存旋转角度和移动量，“：＝”符号表示左边的数值为右边对应数值的更新值；5、判断|θy|+|θx|是否小于δ，如果是，则装配转动角度非常小，装配稳定，转7，否则转6；6、利用式(12)(13)(14)更新f1上点坐标为P1(xi,1,yi,1,zi,1,1)，完成一次迭代，并转3进行下一次迭代计算；xi,1:＝xi,1+zi,1·θy(12)yi,1:＝yi,1-zi,1·θx(13)zi,1:＝zi,1-xi,1·θy+yi,1θx+dz(14)7、停止计算，实际装配位姿为H(Sθx,Sθy,Sdz)。有益效果：(1)本专利技术通过分析待装配零件空间平面关系、装配面的形貌特征和外力作用，得到了对装配参数的最优化表达式，该最优化表达式为线性化模型，从而这就将零件装配这样一个非线性的问题转换为线性化问题来进行求解，提高了装配问题的计算速度；且计算所得的装配误差低。(2)本专利技术方法在求解装配的最优化问题中，通过对装配面上点数据空间位置的迭代更新，找到满足设定误差的装配位姿。本专利技术适用于所有精度的装配问题，可将装配误差控制在任意小的范围内。附图说明图1本专利技术方案流程图。具体实施方式下面结合附图并举实施例，对本专利技术进行详细描述。本专利技术提供了一种考虑几何误差的装配位姿高精度线性化求解方法，适用于刚性零件精密装配中的位姿计算，面向装配面为平面或者任意曲面的装配关系。不限于精密装配，装配角度可以变化更大甚至任意。在装配过程中，每一个装配关系包含了两个装配面，即两个零件的外表面。零件表面非理想平面，有不规则的凹凸、缺陷等形貌特征，这导致了在微观尺度上，两平面的实际的接触形式是点接触与平面接触的混合接触形式。正是这些实际接触点与接触面的高低、位置等导致同一批次的每个装配件有不同的装配位姿，构成了装配误差的一部分。在装配精度要求高的场合下，这些形貌特征造成的接触位姿关系将会对装配精度产生很大影响。且在装配体包含的零件较多时，所有装配关系中接触位姿关系引起的装配误差将会耦合，或累加增大，或者抵消减小，由此接触位姿关系对装配精度的不确定度影响进一步扩大，所以接触位姿关系是在精密装配中必须考虑的因素。但计算得到接触位姿关系同时也要兼顾计算速度，因为如果计算不够快速，机械系统的装配优化无法实现。所以要做到装配工艺模拟与工艺参数优化，就需要一种方法可以准确且快速地求解零部件之间接触位姿关系。因此可以通过测量装配平面，参数化描述形貌特征，迭代求解接触位姿关系，优化参数制定装配方案，从而控制装配误差、保证装配体精度指标。如图1所示，本专利技术方法具体包括以下步骤：步骤一、测量零件上的待装配面，对待装配面进行参数化处理，得到每个待装配面面形的数学描述；本步骤的目的是获得单个装配面的真实形貌和位姿信息。实际操作方法如下：先选择零件上的装配面进行测量，例如，使用三坐标测量仪测得装配面上呈网格分布的点的坐标，坐标系汇总为坐标集Rn＝{(xnl,ynl,znl)l}，其中n是装配面编号n＝1,2，所述装配面记为装配中接触的两个面，l是装配面上被测点的编号。再对坐标集Rn进行曲线插值，根据曲面插值得到装配前零件上装配面的面型Pn＝(xi,n,yi,n,zi,n,1)，xl,n,yl,n,zl,n为待装配面n上i位置的点数据，多出的1是为了后续进行坐标的空间变换方便而添加的。如此便对每个零件上的每个装配面进行参数化处理，得到每个装配面的数学描述。步骤二、根据待装配面的空间平面、形貌误差和受力情况得到待装配面装配位姿的最优化表达式；本步骤目的是考虑装配面形貌特征造成的接触位姿关系，求解两个接触面之间的位姿关系。首先不考虑接触关系和受力情况，仅考虑空间平面关系。有公式：H·P1＝P2(1)装配面1用f1表示，装配面2用f2表示，其中，f1为运动面，f2不动面，P1(xi,1,yi,1,zi,1,1)是f1上点集,P2(xi,2,yi,2,zi,2,1)是面f2上的点集。f2与f1的位姿关系通过空间运动变换矩阵进行描述，f1上的点集P1经过空间变换H变成点集P2。H是空间运动矩阵，为空间平面相对于装配坐标系的三维空间变换矩阵，H＝MxMyMz；Mx、My、Mz分别为测量坐标系相对于装配坐标系的x、y和z轴方向的转动和平移，其中，θx，θy和θz分别表示x、y和z轴方向的转动量；dx，dy和dz分别表示x、y和z轴方向的平移量；实际装配时需要限制自由度，这里限制x方向的移动，y方向的移动，z方向的转动即：dx＝0；dy＝0；θz＝0.则空间运动矩阵H简化为：将式(5)带入式(1)，展开有：面向精密装配过程，因装配面的表面质量较好，满足平面连续假设，并且无突变，所以装配角θx,θy非常小，由将式(6)简化得：忽略二阶及其以上无穷小，式(7)进一步简化得：装配时坐标系应该置于f1的y轴的中心处，此时zi,1值最小，对于本文档来自技高网...

【技术保护点】
1.一种考虑几何误差的装配位姿高精度线性化求解方法，其特征在于，包括以下步骤：步骤一、测量零件上待装配面的形貌并进行参数化处理，得到每个待装配面面形的数学描述；步骤二、根据待装配面的空间平面关系、形貌误差和受力情况得到待装配面装配位姿的最优化表达式；步骤三、求解所述最优化表达式，得到待装配面的最优装配位姿。

【技术特征摘要】
1.一种考虑几何误差的装配位姿高精度线性化求解方法，其特征在于，包括以下步骤：步骤一、测量零件上待装配面的形貌并进行参数化处理，得到每个待装配面面形的数学描述；步骤二、根据待装配面的空间平面关系、形貌误差和受力情况得到待装配面装配位姿的最优化表达式；步骤三、求解所述最优化表达式，得到待装配面的最优装配位姿。2.如权利要求1所述的一种考虑几何误差的装配位姿高精度线性化求解方法，其特征在于，记待装配的两个面为f1和f2，其中，f1为运动面，f2不动面，所述最优化表达式为：其中，ρi为第i点的分布力，θx和θy分别表示x轴和y轴方向的转动量；dz表示z轴方向的平移量；(xi,1,yi,1,zi,1)是待装配面f1上点集,(xi,2,yi,2,zi,2)是待装配面f2上的点集。3.如权利要求2所述的一种考虑几何误差的装配位姿高精度线性化求解方法，其特征在于，步骤三的求解过程中，将实时更新的xi,1，yi,1和zi,1带入最优化表达式中，求解该优化问题。4.如权利要求3所述的一种考虑几何误差的装配位姿高精度线性化求解方法，其特征在于，迭代计算的过程如下...

【专利技术属性】
技术研发人员：金鑫，刘志华，张之敬，尚可，王子夫，
申请(专利权)人：北京理工大学，
类型：发明
国别省市：北京,11