基于实测数据的舱门边框轮廓铣切余量测算方法技术

本发明专利技术公开了一种基于实测数据的舱门边框轮廓铣切余量测算方法，包括：通过三维扫描设备获取实测舱门边框表面三维点云数据；对实测舱门边框表面三维点云数据预处理；提取实测数据交点和定位孔圆心；基于舱门交点和定位孔特征将实测三维数据与CAD数模进行由粗到精配准；计算实测舱门边框轮廓点；计算配准后实测边框轮廓点至CAD数模的欧式距离，此距离即为舱门边框理论铣切余量。所提出的方法，可以准确计算舱门边框轮廓点的铣切余量，指导零部件精确加工，从而实现飞机舱门的精确装配。

Calculation method of milling allowance of hatch frame based on measured data

【技术实现步骤摘要】
基于实测数据的舱门边框轮廓铣切余量测算方法
本专利技术涉及航天测量
，具体而言涉及一种基于实测数据的舱门边框轮廓铣切余量测算方法。
技术介绍
飞机是航天领域的运输和作业的重要交通工具，舱门边框是属于飞机最大部件机体的一部分，制作舱门边框时，要对初加工成的边框计算铣切余量，输入数控机床铣切掉多余的金属来实现精准装配，然而铣切余量的计算大多由人工完成，测量每个基准点周围的铣切余量，测量工作量大，定位不准，测量精度低，容易导致错误测量，需要多次和铣切，而且边框的大小和材质导致舱门的铣切难度高，浪费人力和物力资源。使用三维扫描仪采集舱门边框的三维点云数据，将舱门边框点云数据与CAD模型进行精准匹配，就能准确计算整个舱门边框轮廓点的铣切余量，指导零部件精确加工，从而实现飞机舱门的精确装配。
技术实现思路
本专利技术目的在于提供一种基于实测数据的舱门边框轮廓铣切余量测算方法，可以准确计算舱门边框轮廓点的铣切余量，指导零部件精确加工，从而实现飞机舱门的精确装配，解决了舱门边框轮廓点铣切余量计算需求。为达成上述目的，结合图1，本专利技术提出一种基于实测数据的舱门边框轮廓铣切余量测算方法，所述测算方法包括：S1：采用三维扫描设备获取实测舱门边框表面三维点云数据；S2：对获取的实测舱门边框表面三维点云数据进行预处理，去除离群点；S3：选取舱门交点和定位圆孔处的点云，用最小二乘法拟合舱门交点圆孔特征；S4：基于提取的舱门交点圆孔和定位孔特征将实测三维数据与CAD数模粗配准后，再进行加权精配准；S5：计算得到实测舱门边框轮廓点；S6：计算得到配准后实测边框轮廓点至CAD数模的欧式距离，将此距离作为舱门边框理论铣切余量。进一步的实施例中，步骤S2中，所述对获取的实测舱门边框表面三维点云数据进行预处理，去除离群点的过程包括以下步骤：S21：对于扫描获取的舱门边框表面三维点云数据，依次针对每个点选取k个邻近点，计算该点到其对应的k个邻近点的平均距离，计算得到的所有点到其对应的k个邻近点的平均距离分布符合高斯分布，计算得到该分布的均值和方差；S22：依次判断每个点到其对应的k个邻近点的平均距离是否在均值和方差确定的范围之外，如果在均值和方差确定的范围之外，将该点判定为离群点，从点云数据集中删除该点。进一步的实施例中，步骤S3中，所述选取舱门交点和定位圆孔处的点云，用最小二乘法拟合舱门交点圆孔特征是指，在舱门边框轮廓三维数据中选取交点圆孔和定位圆孔附近的点，采用最小二乘法将选取点拟合成圆，得到圆心和半径。进一步的实施例中，步骤S3中，所述在舱门边框轮廓三维数据中选取交点圆孔和定位圆孔附近的点，采用最小二乘法将选取点拟合成圆，得到圆心和半径的过程包括以下步骤：S31：在其中一个交点圆孔或定位圆孔附近选取n个点Ni，假设圆心为O，半径为R；S32：设αi为点Ni到圆心O的距离的平方，令S33：求F()对圆心坐标和半径的偏导数，令偏导数为0；S34：对前述方程进行求解得到圆心O的坐标和半径R的取值；S35：重复步骤S31-S34，直至计算得到所有交点圆孔和所有定位圆孔的特征信息。进一步的实施例中，步骤S4中，所述基于提取的舱门交点圆孔和定位孔特征将实测三维数据与CAD数模粗配准的过程包括以下步骤：S401：根据舱门边框轮廓三维数据中的圆孔圆心O1，在CAD模型上提取对应位置圆孔的圆孔圆心O2；S402：求出O1到O2的刚性变换矩阵；S403：采用变换矩阵对整个三维数据进行刚性变换。进一步的实施例中，步骤S402中，采用SVD分解法得到刚性变换矩阵。进一步的实施例中，步骤S4中，所述进行加权精配准的过程包括以下步骤：S411：令实测三维点云数据为P，CAD数模为Q；S412：随机选取实测三维点云数据中的n个点生成目标点云集合P1；设pi∈P1，寻找Q上最接近pi的点qi使得||pi-qi||最小，设qi对应的集合为Q1，集合P1和Q1中所有对应点的平均距离为dk：S413：计算P1到Q1的变换矩阵R，对P1集进行坐标变换得到新的点集P′，求P′到Q1的平均距离dk+1：式中，p′i∈P′，为p变换后的对应点，q′i为Q1中距p′i最近的点。令Δd＝dk+1-dk作为精配缩减误差，D为粗配误差：式中，o1i为舱门边框轮廓三维数据中的圆孔圆心，o2i为CAD模型上提取对应位置圆孔的圆孔圆心，各有4个点。S414：判断αD+βΔd否满足αD+βΔd<γ，γ>0，α，β为误差权重，避免铣切余量分布不均匀时导致4个圆心定位点大幅度偏移：如果满足不等式，对P变换，结束流程，否则，转入S412，重新执行精配准流程。进一步的实施例中，步骤S5中，所述计算得到实测舱门边框轮廓点的过程包括以下步骤：计算所有点的法线矢量，依次判断每个点的邻近点中是否存在两个点的法矢方向夹角为60°-120°，如果存在，判定该点为轮廓点；其中，优先访问已确定轮廓点所对应的未被访问的邻近点。进一步的实施例中，步骤S5中，所述计算得到实测舱门边框轮廓点的过程还包括以下步骤：S51：删去距离过近的点，对优先访问队列进行滤波；S52：搜索优先访问队列：若优先访问队列不为空，取队首点为选取点b，若优先访问队列为空，随机获取其中一个未访问点为选取点b，若无未访问点，转入S55；S53：对选取点b进行标记访问，在半径r范围内寻找邻近点并记录，选取点b与对应的最近点b0组成基准向量bb0，其余邻近点bi与b0组成向量b0bi，当bbi与基准向量成30°-150°夹角时，计算点b所有bbi与bb0的外积向量，记为xbbi；S54：若在xbbi中存在两个向量的夹角在60°-120°之间，将选取点b记为轮廓点，选取点b的全部未标记访问的邻近点加载至优先访问队列，转入S52，直至轮廓点数量达到预估数量且累计访问点数达到设定比例或者连续多次访问非轮廓点；S55：结束流程。进一步的实施例中，所述计算得到配准后实测边框轮廓点至CAD数模的欧式距离，将此距离作为舱门边框理论铣切余量的过程包括以下步骤：计算所有舱门边框轮廓点和其邻近轮廓点的欧式距离的平均值，将其作为该点位置的铣切余量。以上本专利技术的技术方案，与现有相比，其显著的有益效果在于：(1)可以准确计算舱门边框轮廓点的铣切余量，指导零部件精确加工，从而实现飞机舱门的精确装配，解决了舱门边框轮廓点铣切余量计算需求(2)采用滤波器去除离群点声，提高了轮廓点估计的准确性，使用定位点粗配准，使用ICP精配准，挺高了铣切余量计算的精准度，快速得到整个舱门边框轮廓每个点的铣切余量。应当理解，前述构思以及在下面更加详细地描述的额外构思的所有组合只要在这样的构思不相互矛盾的情况下都本文档来自技高网...

【技术保护点】
1.一种基于实测数据的舱门边框轮廓铣切余量测算方法，其特征在于，所述测算方法包括：/nS1：采用三维扫描设备获取实测舱门边框表面三维点云数据；/nS2：对获取的实测舱门边框表面三维点云数据进行预处理，去除离群点；/nS3：选取舱门交点和定位圆孔处的点云，用最小二乘法拟合舱门交点圆孔特征；/nS4：基于提取的舱门交点圆孔和定位孔特征将实测三维数据与CAD数模粗配准后，再进行加权精配准；/nS5：计算得到实测舱门边框轮廓点；/nS6：计算得到配准后实测边框轮廓点至CAD数模的欧式距离，将此距离作为舱门边框理论铣切余量。/n

【技术特征摘要】
1.一种基于实测数据的舱门边框轮廓铣切余量测算方法，其特征在于，所述测算方法包括：
S1：采用三维扫描设备获取实测舱门边框表面三维点云数据；
S2：对获取的实测舱门边框表面三维点云数据进行预处理，去除离群点；
S3：选取舱门交点和定位圆孔处的点云，用最小二乘法拟合舱门交点圆孔特征；
S4：基于提取的舱门交点圆孔和定位孔特征将实测三维数据与CAD数模粗配准后，再进行加权精配准；
S5：计算得到实测舱门边框轮廓点；
S6：计算得到配准后实测边框轮廓点至CAD数模的欧式距离，将此距离作为舱门边框理论铣切余量。


2.根据权利要求1所述的基于实测数据的舱门边框轮廓铣切余量测算方法，其特征在于，步骤S2中，所述对获取的实测舱门边框表面三维点云数据进行预处理，去除离群点的过程包括以下步骤：
S21：对于扫描获取的舱门边框表面三维点云数据，依次针对每个点选取k个邻近点，计算该点到其对应的k个邻近点的平均距离，计算得到的所有点到其对应的k个邻近点的平均距离分布符合高斯分布，计算得到该分布的均值和方差；
S22：依次判断每个点到其对应的k个邻近点的平均距离是否在均值和方差确定的范围之外，如果在均值和方差确定的范围之外，将该点判定为离群点，从点云数据集中删除该点。


3.根据权利要求1所述的基于实测数据的舱门边框轮廓铣切余量测算方法，其特征在于，步骤S3中，所述选取舱门交点和定位圆孔处的点云，用最小二乘法拟合舱门交点圆孔特征是指，
在舱门边框轮廓三维数据中选取交点圆孔和定位圆孔附近的点，采用最小二乘法将选取点拟合成圆，得到圆心和半径。


4.根据权利要求3所述的基于实测数据的舱门边框轮廓铣切余量测算方法，其特征在于，步骤S3中，所述在舱门边框轮廓三维数据中选取交点圆孔和定位圆孔附近的点，采用最小二乘法将选取点拟合成圆，得到圆心和半径的过程包括以下步骤：
S31：在其中一个交点圆孔或定位圆孔附近选取n个点Ni，假设圆心为O，半径为R；
S32：设αi为点Ni到圆心O的距离的平方，令
S33：求F()对圆心坐标和半径的偏导数，令偏导数为0；
S34：对前述方程进行求解得到圆心O的坐标和半径R的取值；
S35：重复步骤S31-S34，直至计算得到所有交点圆孔和所有定位圆孔的特征信息。


5.根据权利要求1所述的基于实测数据的舱门边框轮廓铣切余量测算方法，其特征在于，步骤S4中，所述基于提取的舱门交点圆孔和定位孔特征将实测三维数据与CAD数模粗配准的过程包括以下步骤：
S401：根据舱门边框轮廓三维数据中的圆孔圆心O1，在CAD模型上提取对应位置圆孔的圆孔圆心O2；
S402：求出O1到O2的刚性变换矩阵；
S403：采用变换矩阵对整个三维数据进行刚性变换。


6.根据权利要求5所述的基于实测数据的舱门边框轮廓铣切余量测算方法，其特征在于，步骤S402中，采用SVD分解法得到刚性变换矩阵。


7.根据权利要...

【专利技术属性】
技术研发人员：汪俊，陈红华，魏泽勇，张沅，
申请(专利权)人：南京航空航天大学，
类型：发明
国别省市：江苏;32