CN108427837B - 一种装配粗差判定方法及系统 - Google Patents

一种装配粗差判定方法及系统 Download PDF

Info

Publication number
CN108427837B
CN108427837B CN201810182694.3A CN201810182694A CN108427837B CN 108427837 B CN108427837 B CN 108427837B CN 201810182694 A CN201810182694 A CN 201810182694A CN 108427837 B CN108427837 B CN 108427837B
Authority
CN
China
Prior art keywords
assembly deviation
assembly
simulation data
data
uncertainty
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810182694.3A
Other languages
English (en)
Other versions
CN108427837A (zh
Inventor
朱永国
霍正书
周结华
徐雪峰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanchang Hangkong University
Original Assignee
Nanchang Hangkong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanchang Hangkong University filed Critical Nanchang Hangkong University
Priority to CN201810182694.3A priority Critical patent/CN108427837B/zh
Publication of CN108427837A publication Critical patent/CN108427837A/zh
Application granted granted Critical
Publication of CN108427837B publication Critical patent/CN108427837B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/18Manufacturability analysis or optimisation for manufacturability

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Complex Calculations (AREA)

Abstract

本发明公开一种装配粗差判定方法及系统。该方法包括:获取装配配合面的实测装配偏差量;分析装配偏差的变动范围;对变动范围进行有限元仿真,得到装配偏差仿真数据;去除装配偏差仿真数据中的极限值,得到置信装配偏差仿真数据;将置信装配偏差仿真数据进行分组,得到多个装配偏差仿真数据组;计算每个装配偏差仿真数据组的不确定度区间,得到装配偏差合格区间;将每个装配偏差仿真数据组和实测装配偏差量组成的实测数据集合融合,得到多个融合数据集合;计算包含同一实测装配偏差量的所有融合数据的平均不确定度;确定位于装配偏差合格区间之外的平均不确定度所对应的实测装配偏差量为粗差。本发明的技术方案能够提高粗差识别和筛选的准确度。

Description

一种装配粗差判定方法及系统
技术领域
本发明涉及飞机制造技术领域,特别是涉及一种装配粗差判定方法及系统。
背景技术
薄壁件变形回弹、铆接干涉等不利因素,导致飞机装配偏差存在较大的不确定性。近年来,随着激光跟踪仪、局域GPS、CCD工业相机等数字化测量设备的快速发展,在飞机研制过程中引入数字化测量设备已经成为国内外飞机制造公司的共识,利用数字化测量设备可测得装配尺寸数据。充分利用和挖掘装配尺寸数据中所包含的信息,对于筛选出处于关键地位的装配偏差影响因素、抓住飞机装配精度控制的关键要素、揭示装配偏差的传播机理、减少装配过程异常波动、控制和改进产品装配质量具有重要的意义。然而飞机生产批量小等因素导致无法观测到装配过程中各类偏差输入、传递与输出的大量完备检测数据,因此实际的测量数据呈现出小样本、高维、不完备等特点,导致难以用传统的统计方法对实际的测量数据进行建模和分析。经典的测量数据粗差识别准则是基于统计判别法等数理统计方法,明显不适用于飞机小样本装配粗差识别。随着信息论研究的深入和发展,信息论与测量理论之间不断相互渗透发展,形成了测量信息论。虽然测量信息论能挖掘小样本测量数据信息,但由于复杂结构产品装配偏差不确定性大,如果仅基于测量信息论进行粗差识别,粗差判定的置信区间难以划定。
发明内容
本发明的目的是提供一种装配粗差判定方法及系统,提高粗差判定的置信区间的准确度,提高粗差识别和筛选的准确度。
为实现上述目的,本发明提供了如下方案:
一种装配粗差判定方法,包括:
获取装配配合面的实测装配偏差量;
利用先验经验和知识,分析装配偏差的变动范围;
利用蒙特卡洛法对所述装配变差的变动范围进行装配偏差有限元仿真,得到装配偏差仿真数据;
利用顺序统计、秩估计和信息熵相结合的方法去除所述装配偏差仿真数据中的极限值,得到置信装配偏差仿真数据;
以所述实测装配偏差量的数量为一个数据组内的数据的数量,将所述置信装配偏差仿真数据进行分组,得到多个装配偏差仿真数据组;
利用顺序统计、秩估计和信息熵结合的方法计算每个所述装配偏差仿真数据组的不确定度区间,得到装配偏差合格区间;
将每个所述装配偏差仿真数据组和所述实测装配偏差量组成的实测数据集合融合,得到多个融合数据集合;每个所述融合数据集合均由一个装配偏差仿真数据组和实测数据集合相加形成;
计算包含同一所述实测装配偏差量的所有融合数据的平均不确定度;
将所述平均不确定度与所述装配偏差合格区间对比,确定包含在所述装配偏差合格区间内的所述平均不确定度所对应的实测装配偏差量为合格偏差,确定位于所述装配偏差合格区间之外的所述平均不确定度所对应的实测装配偏差量为粗差。
可选的,所述获取装配配合面的实测装配偏差量,具体包括:
获取数字化测量设备采集的装配配合面图像;
对所述装配配合面图像进行处理,得到装配特征的实测数据;
将所述装配特征的实测数据与装配特征的理论数据进行对比,得到实测装配偏差量。
可选的,所述利用顺序统计、秩估计和信息熵相结合的方法去除所述装配偏差仿真数据中的极限值,得到置信装配偏差仿真数据,具体包括:
将所述装配偏差仿真数据按从小到大的顺序排列,得到有序装配偏差仿真数据集合;
依次计算所述有序装配偏差仿真数据集合中每个数据的平均秩估计;
根据所述有序装配偏差仿真数据集合中每个数据的平均秩估计计算所述有序装配偏差仿真数据集合的信息熵;
利用所述有序装配偏差仿真数据集合的信息熵计算所述有序装配偏差仿真数据集合的不确定度;
根据所述有序装配偏差仿真数据集合的中所有数据的平均值和所述有序装配偏差仿真数据集合的不确定度确定所述有序装配偏差仿真数据集合的置信区间;
判断所述装配偏差仿真数据中是否存在位于所述置信装配偏差仿真数据区间之外的数据,得到第一判断结果;
若所述第一判断结果表示不存在,则确定所述装配偏差仿真数据为置信装配偏差仿真数据;
若所述第一判断结果表示存在,则将所述装配偏差仿真数据中位于所述置信装配偏差仿真数据区间之外的数据去除之后返回计算所述有序装配偏差仿真数据集合的步骤。
可选的,所述利用顺序统计、秩估计和信息熵结合的方法计算每个所述装配偏差仿真数据组的不确定度区间,得到装配偏差合格区间,具体包括:
将每个所述装配偏差仿真数据组中的数据按从小到大的顺序排列,得到有序装配偏差仿真数据组;
依次计算每个所述有序装配偏差仿真数据组中每个数据的平均秩估计;
根据每个所述有序装配偏差仿真数据组中每个数据的平均秩估计计算每个所述有序装配偏差仿真数据组的信息熵;
利用每个所述有序装配偏差仿真数据组的信息熵计算每个所述有序装配偏差仿真数据组的不确定度;
在每个所述有序装配偏差仿真数据组的不确定度中,筛选出最小值作为所述装配偏差合格区间的最小值,筛选出最大值作为所述装配偏差合格区间的最大值,得到所述装配偏差合格区间。
可选的,所述计算包含同一所述实测装配偏差量的所有融合数据的平均不确定度,具体包括:
计算每个所述融合数据集合的不确定度;
将包含同一所述实测装配偏差量的所有融合数据集合划分到同一个组中,得到多个融合数据组;
求取包含在同一个融合数据组中的所有融合数据集合的不确定度的平均值,得到各个融合数据组的平均不确定度。
本发明还公开一种装配粗差判定系统,包括:
获取模块,用于获取装配配合面的实测装配偏差量;
变动范围确定模块,用于利用先验经验和知识,分析装配偏差的变动范围;
仿真模块,用于利用蒙特卡洛法对所述装配变差的变动范围进行装配偏差有限元仿真,得到装配偏差仿真数据;
置信数据筛选模块,用于利用顺序统计、秩估计和信息熵相结合的方法去除所述装配偏差仿真数据中的极限值,得到置信装配偏差仿真数据;
分组模块,用于以所述实测装配偏差量的数量为一个数据组内的数据的数量,将所述置信装配偏差仿真数据进行分组,得到多个装配偏差仿真数据组;
合格区间确定模块,用于利用顺序统计、秩估计和信息熵结合的方法计算每个所述装配偏差仿真数据组的不确定度区间,得到装配偏差合格区间;
融合模块,用于将每个所述装配偏差仿真数据组和所述实测装配偏差量组成的实测数据集合融合,得到多个融合数据集合;每个所述融合数据集合均由一个装配偏差仿真数据组和实测数据集合相加形成;
平均不确定度计算模块,用于计算包含同一所述实测装配偏差量的所有融合数据的平均不确定度;
粗差判定模块,用于将所述平均不确定度与所述装配偏差合格区间对比,确定包含在所述装配偏差合格区间内的所述平均不确定度所对应的实测装配偏差量为合格偏差,确定位于所述装配偏差合格区间之外的所述平均不确定度所对应的实测装配偏差量为粗差。
可选的,所述获取模块,具体包括:
图像获取单元,用于获取数字化测量设备采集的装配配合面图像;
图像处理单元,用于对所述装配配合面图像进行处理,得到装配特征的实测数据;
偏差计算单元,用于将所述装配特征的实测数据与装配特征的理论数据进行对比,得到实测装配偏差量。
可选的,所述置信数据筛选模块,具体包括:
第一排列单元,用于将所述装配偏差仿真数据按从小到大的顺序排列,得到有序装配偏差仿真数据集合;
第一平均秩估计计算单元,用于依次计算所述有序装配偏差仿真数据集合中每个数据的平均秩估计;
第一信息熵计算单元,用于根据所述有序装配偏差仿真数据集合中每个数据的平均秩估计计算所述有序装配偏差仿真数据集合的信息熵;
第一不确定计算单元,用于利用所述有序装配偏差仿真数据集合的信息熵计算所述有序装配偏差仿真数据集合的不确定度;
置信区间确定单元,用于根据所述有序装配偏差仿真数据集合的中所有数据的平均值和所述有序装配偏差仿真数据集合的不确定度确定所述有序装配偏差仿真数据集合的置信区间;
第一判断单元,用于判断所述装配偏差仿真数据中是否存在位于所述置信装配偏差仿真数据区间之外的数据,得到第一判断结果;
置信数据确定单元,用于若所述第一判断结果表示不存在,则确定所述装配偏差仿真数据为置信装配偏差仿真数据;
返回单元,用于若所述第一判断结果表示存在,则将所述装配偏差仿真数据中位于所述置信装配偏差仿真数据区间之外的数据去除之后返回计算所述有序装配偏差仿真数据集合的步骤。
可选的,所述合格区间确定模块,具体包括:
第二排列单元,用于将每个所述装配偏差仿真数据组中的数据按从小到大的顺序排列,得到有序装配偏差仿真数据组;
第二平均秩估计计算单元,用于依次计算每个所述有序装配偏差仿真数据组中每个数据的平均秩估计;
第二信息熵计算单元,用于根据每个所述有序装配偏差仿真数据组中每个数据的平均秩估计计算每个所述有序装配偏差仿真数据组的信息熵;
第二不确定计算单元,用于利用每个所述有序装配偏差仿真数据组的信息熵计算每个所述有序装配偏差仿真数据组的不确定度;
合格区间确定单元,用于在每个所述有序装配偏差仿真数据组的不确定度中,筛选出最小值作为所述装配偏差合格区间的最小值,筛选出最大值作为所述装配偏差合格区间的最大值,得到所述装配偏差合格区间。
可选的,所述平均不确定度计算模块,具体包括:
第三不确定度计算单元,用于计算每个所述融合数据集合的不确定度;
分组单元,用于将包含同一所述实测装配偏差量的所有融合数据集合划分到同一个组中,得到多个融合数据组;
平均不确定度计算单元,用于求取包含在同一个融合数据组中的所有融合数据集合的不确定度的平均值,得到各个融合数据组的平均不确定度。
根据本发明提供的具体实施例,本发明公开了以下技术效果:本发明利用先验经验和知识分析装配偏差的变动范围,然后对所述装配变差的变动范围进行装配偏差有限元仿真,有效的实现了装配偏差先验经验和知识的挖掘,为置信装配偏差仿真数据的计算提供了有力支撑,提高了置信装配偏差仿真数据计算的准确度,进一步提高了粗差判定的准确度;同时将实测装配偏差量和装配偏差仿真数据融合到一起,避免了仅基于一种数据进行粗差判定的局限性,有效解决了飞机等复杂结构产品装配过程仅依靠实测数据进行粗差判定时装配偏差合格区间难以确定的难题,可在装配累积偏差不确定性大的情况下实现装配测量粗差判定,提高了在装配累积偏差不确定性大的情况下的装配测量粗差判定的准确度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明装配粗差判定方法实施例的方法流程图;
图2为本发明装配粗差判定系统实施例的系统结构图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种装配粗差判定方法及系统,提高粗差判定的置信区间的准确度,提高粗差识别和筛选的准确度。本发明针对装配偏差测量数据少、累积偏差值不确定性大的特点,提出测量信息论和有限元法相融合的装配粗差判定方法。首先,基于先验经验与知识,利用有限元分析和蒙特卡罗法获取装配偏差先验信息。然后,挖掘偏差检测数据,形成测量数据信息。再将偏差先验信息与测量数据信息相融合,制定测量信息论和有限元分析法相融合的装配粗差判定法则。最终实现多源信息利用和数据融合,形成基于知识的小样本装配偏差数据挖掘方法。本发明特点在于:1)基于先验经验与知识,利用有限元分析和蒙特卡罗法获取装配偏差先验信息;2)将偏差先验信息与测量数据信息相融合,制定测量信息论和有限元分析法相融合的装配粗差判定法则。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明装配粗差判定方法实施例的方法流程图。
该装配粗差判定方法,包括:
步骤101:获取装配配合面的实测装配偏差量。
该步骤101具体包括:
获取数字化测量设备采集的装配配合面图像;对所述装配配合面图像进行处理,得到装配特征的实测数据;将所述装配特征的实测数据与装配特征的理论数据进行对比,得到实测装配偏差量。数字化测量设备包括CCD工业相机组等。图像处理过程包括预处理、图像分割、图像二值化、边缘检测和图像目标提取等操作。该步骤101中装配偏差试验次数记为n。
步骤102:利用先验经验和知识,分析装配偏差的变动范围;装配偏差包括零件制造误差、工装定位误差、薄壁件变形回弹、铆接变形、装配基准偏差等。
步骤103:利用蒙特卡洛法对所述装配变差的变动范围进行m次装配偏差有限元仿真,m=kn(k≥10),得到装配偏差仿真数据ym={y1a,y2a,…,yma},其中a表示装配偏差仿真数据的维度,表示装配偏差的不同描述方面。
步骤104:为了使仿真得到的数据更具稳健性,避免极端值对计算不确定度的影响,针对各装配偏差的不同描述方面,利用顺序统计、秩估计和信息熵相结合的方法去除所述装配偏差仿真数据中的极限值,得到置信装配偏差仿真数据。
该步骤104具体包括:
设描述装配的某一方面偏差仿真数据为
Figure BDA0001589349280000081
将所述装配偏差仿真数据
Figure BDA0001589349280000082
按从小到大的顺序排列,得到有序装配偏差仿真数据集合
Figure BDA0001589349280000083
依次计算所述有序装配偏差仿真数据集合
Figure BDA0001589349280000084
中每个数据的平均秩估计
Figure BDA0001589349280000085
根据所述有序装配偏差仿真数据集合
Figure BDA0001589349280000086
中每个数据的平均秩估计
Figure BDA0001589349280000087
计算所述有序装配偏差仿真数据集合的信息熵
Figure BDA0001589349280000088
利用所述有序装配偏差仿真数据集合
Figure BDA0001589349280000089
的信息熵
Figure BDA00015893492800000810
计算所述有序装配偏差仿真数据集合的不确定度
Figure BDA00015893492800000811
根据所述有序装配偏差仿真数据集合的中所有数据的平均值和所述有序装配偏差仿真数据集合的不确定度确定所述有序装配偏差仿真数据集合的置信区间;将置信度确定为0.95,置信度为0.95以上的置信区间为
Figure BDA00015893492800000812
Figure BDA00015893492800000813
表示数据组
Figure BDA00015893492800000814
的平均值;
判断所述装配偏差仿真数据中是否存在位于所述置信装配偏差仿真数据区间之外的数据,得到第一判断结果;
若所述第一判断结果表示不存在,则确定所述装配偏差仿真数据为置信装配偏差仿真数据;
若所述第一判断结果表示存在,则将所述装配偏差仿真数据中位于所述置信装配偏差仿真数据区间之外的数据去除之后返回计算所述有序装配偏差仿真数据集合的步骤。
步骤105:以所述实测装配偏差量的数量为一个数据组内的数据的数量,将所述置信装配偏差仿真数据进行分组,得到多个装配偏差仿真数据组。
首先针对每一种装配偏差描述,将步骤104中得到的置信装配偏差仿真数据再次进行顺序统计,记为β12,…,βl(其中l<m),然后将β12,…,βl的第1至
Figure BDA0001589349280000091
项分成
Figure BDA0001589349280000092
组,“[]”表示取整运算,每组n个数据,并记为N1,N2,…,Nt,其中
Figure BDA0001589349280000093
并将剩下的
Figure BDA0001589349280000094
项数据记为Np,其中p=t+1,从而得到多个装配偏差仿真数据组N1,N2,…,Np
步骤106:利用顺序统计、秩估计和信息熵结合的方法计算每个所述装配偏差仿真数据组的不确定度区间,得到装配偏差合格区间。
该步骤106具体包括:计算装配偏差仿真数据组的不确定度和根据装配偏差仿真数据组的不确定度计算不确定度区间。
其中计算装配偏差仿真数据组的不确定度的过程与步骤步骤104中计算不确定度的原理相同,具体过程为:
将每个所述装配偏差仿真数据组中的数据按从小到大的顺序排列,得到有序装配偏差仿真数据组;
依次计算每个所述有序装配偏差仿真数据组中每个数据的平均秩估计;
根据每个所述有序装配偏差仿真数据组中每个数据的平均秩估计计算每个所述有序装配偏差仿真数据组的信息熵;
利用每个所述有序装配偏差仿真数据组的信息熵计算每个所述有序装配偏差仿真数据组的不确定度。多个装配偏差仿真数据组N1,N2,…,Np分别对应的不确定度为Ui(i=1,2,…,p)。
根据装配偏差仿真数据组的不确定度计算不确定度区间的具体过程为:
在每个所述有序装配偏差仿真数据组的不确定度中,筛选出最小值作为所述装配偏差合格区间的最小值,筛选出最大值作为所述装配偏差合格区间的最大值,得到所述装配偏差合格区间[Umin,Umax]。
步骤107:将每个所述装配偏差仿真数据组和所述实测装配偏差量组成的实测数据集合融合,得到多个融合数据集合Mbc=Nb+{xc},b=1,2,L,p;c=1,2,…,n。其中x表示实验数据,b表示有限元仿真数据演变出来的小样本的次序,即装配偏差仿真数据组的次序,c表示实验数据的次序;每个所述融合数据集合均由一个装配偏差仿真数据组和实测数据集合相加形成。
步骤108:计算包含同一所述实测装配偏差量的所有融合数据的平均不确定度。
该步骤108具体包括:具体包括:
计算每个所述融合数据集合Mbc的不确定度Ubc;融合数据集合的不确定度的计算过程与步骤步骤104中计算不确定度的原理相同;
将包含同一所述实测装配偏差量xi的所有融合数据集合划分到同一个组中,得到多个融合数据组;
求取包含在同一个融合数据组中的所有融合数据集合的不确定度的平均值,得到各个融合数据组的平均不确定度
Figure BDA0001589349280000101
并以
Figure BDA0001589349280000102
作为判断xi是否是粗差的依据。
步骤109:将所述平均不确定度
Figure BDA0001589349280000103
与所述装配偏差合格区间[Umin,Umax]对比,确定包含在所述装配偏差合格区间内的所述平均不确定度所对应的实测装配偏差量为合格偏差,确定位于所述装配偏差合格区间之外的所述平均不确定度所对应的实测装配偏差量为粗差,并将粗差予以剔除,实现装配粗差源的量化识别。
图2为本发明装配粗差判定系统实施例的系统结构图。
参见图2该装配粗差判定系统,包括:
获取模块201,用于获取装配配合面的实测装配偏差量;该获取模块201,具体包括:图像获取单元、图像处理单元和偏差计算单元,其中:图像获取单元,用于获取数字化测量设备采集的装配配合面图像;图像处理单元,用于对所述装配配合面图像进行处理,得到装配特征的实测数据;偏差计算单元,用于将所述装配特征的实测数据与装配特征的理论数据进行对比,得到实测装配偏差量。
变动范围确定模块202,用于利用先验经验和知识,分析装配偏差的变动范围。
仿真模块203,用于利用蒙特卡洛法对所述装配变差的变动范围进行装配偏差有限元仿真,得到装配偏差仿真数据。
置信数据筛选模块204,用于利用顺序统计、秩估计和信息熵相结合的方法去除所述装配偏差仿真数据中的极限值,得到置信装配偏差仿真数据;该置信数据筛选模块204,具体包括:第一排列单元、第一平均秩估计计算单元、第一信息熵计算单元、第一不确定计算单元、置信区间确定单元、第一判断单元、置信数据确定单元和返回单元。其中:
第一排列单元,用于将所述装配偏差仿真数据按从小到大的顺序排列,得到有序装配偏差仿真数据集合;
第一平均秩估计计算单元,用于依次计算所述有序装配偏差仿真数据集合中每个数据的平均秩估计;
第一信息熵计算单元,用于根据所述有序装配偏差仿真数据集合中每个数据的平均秩估计计算所述有序装配偏差仿真数据集合的信息熵;
第一不确定计算单元,用于利用所述有序装配偏差仿真数据集合的信息熵计算所述有序装配偏差仿真数据集合的不确定度;
置信区间确定单元,用于根据所述有序装配偏差仿真数据集合的中所有数据的平均值和所述有序装配偏差仿真数据集合的不确定度确定所述有序装配偏差仿真数据集合的置信区间;
第一判断单元,用于判断所述装配偏差仿真数据中是否存在位于所述置信装配偏差仿真数据区间之外的数据,得到第一判断结果;
置信数据确定单元,用于若所述第一判断结果表示不存在,则确定所述装配偏差仿真数据为置信装配偏差仿真数据;
返回单元,用于若所述第一判断结果表示存在,则将所述装配偏差仿真数据中位于所述置信装配偏差仿真数据区间之外的数据去除之后返回计算所述有序装配偏差仿真数据集合的步骤。
分组模块205,用于以所述实测装配偏差量的数量为一个数据组内的数据的数量,将所述置信装配偏差仿真数据进行分组,得到多个装配偏差仿真数据组。
合格区间确定模块206,用于利用顺序统计、秩估计和信息熵结合的方法计算每个所述装配偏差仿真数据组的不确定度区间,得到装配偏差合格区间;该合格区间确定模块206,具体包括:第二排列单元、第二平均秩估计计算单元、第二信息熵计算单元、第二不确定计算单元和合格区间确定单元。其中:
第二排列单元,用于将每个所述装配偏差仿真数据组中的数据按从小到大的顺序排列,得到有序装配偏差仿真数据组;
第二平均秩估计计算单元,用于依次计算每个所述有序装配偏差仿真数据组中每个数据的平均秩估计;
第二信息熵计算单元,用于根据每个所述有序装配偏差仿真数据组中每个数据的平均秩估计计算每个所述有序装配偏差仿真数据组的信息熵;
第二不确定计算单元,用于利用每个所述有序装配偏差仿真数据组的信息熵计算每个所述有序装配偏差仿真数据组的不确定度;
合格区间确定单元,用于在每个所述有序装配偏差仿真数据组的不确定度中,筛选出最小值作为所述装配偏差合格区间的最小值,筛选出最大值作为所述装配偏差合格区间的最大值,得到所述装配偏差合格区间。
融合模块207,用于将每个所述装配偏差仿真数据组和所述实测装配偏差量组成的实测数据集合融合,得到多个融合数据集合;每个所述融合数据集合均由一个装配偏差仿真数据组和实测数据集合相加形成。
平均不确定度计算模块208,用于计算包含同一所述实测装配偏差量的所有融合数据的平均不确定度;所述平均不确定度计算模块208,具体包括:第三不确定度计算单元、分组单元和平均不确定度计算单元。其中:
第三不确定度计算单元,用于计算每个所述融合数据集合的不确定度;
分组单元,用于将包含同一所述实测装配偏差量的所有融合数据集合划分到同一个组中,得到多个融合数据组;
平均不确定度计算单元,用于求取包含在同一个融合数据组中的所有融合数据集合的不确定度的平均值,得到各个融合数据组的平均不确定度。
粗差判定模块209,用于将所述平均不确定度与所述装配偏差合格区间对比,确定包含在所述装配偏差合格区间内的所述平均不确定度所对应的实测装配偏差量为合格偏差,确定位于所述装配偏差合格区间之外的所述平均不确定度所对应的实测装配偏差量为粗差。
根据本发明提供的具体实施例,本发明公开了以下技术效果:本发明利用先验经验和知识分析装配偏差的变动范围,然后对所述装配变差的变动范围进行装配偏差有限元仿真,有效的实现了装配偏差先验经验和知识的挖掘,为置信装配偏差仿真数据的计算提供了有力支撑,提高了置信装配偏差仿真数据计算的准确度,进一步提高了粗差判定的准确度;同时将实测装配偏差量和装配偏差仿真数据融合到一起,避免了仅基于一种数据进行粗差判定的局限性,有效解决了飞机等复杂结构产品装配过程仅依靠实测数据进行粗差判定时装配偏差合格区间难以确定的难题,可在装配累积偏差不确定性大的情况下实现装配测量粗差判定,提高了在装配累积偏差不确定性大的情况下的装配测量粗差判定的准确度。
对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。

Claims (10)

1.一种装配粗差判定方法,其特征在于,包括:
获取装配配合面的实测装配偏差量;
利用先验经验和知识,分析装配偏差的变动范围;
利用蒙特卡洛法对所述装配偏 差的变动范围进行装配偏差有限元仿真,得到装配偏差仿真数据;
利用顺序统计、秩估计和信息熵相结合的方法去除所述装配偏差仿真数据中的极限值,得到置信装配偏差仿真数据;
以所述实测装配偏差量的数量为一个数据组内的数据的数量,将所述置信装配偏差仿真数据进行分组,得到多个装配偏差仿真数据组;
利用顺序统计、秩估计和信息熵结合的方法计算每个所述装配偏差仿真数据组的不确定度区间,得到装配偏差合格区间;
将每个所述装配偏差仿真数据组和所述实测装配偏差量组成的实测数据集合融合,得到多个融合数据集合;每个所述融合数据集合均由一个装配偏差仿真数据组和实测数据集合相加形成;
计算包含同一所述实测装配偏差量的所有融合数据的平均不确定度;
将所述平均不确定度与所述装配偏差合格区间对比,确定包含在所述装配偏差合格区间内的所述平均不确定度所对应的实测装配偏差量为合格偏差,确定位于所述装配偏差合格区间之外的所述平均不确定度所对应的实测装配偏差量为粗差。
2.根据权利要求1所述的一种装配粗差判定方法,其特征在于,所述获取装配配合面的实测装配偏差量,具体包括:
获取数字化测量设备采集的装配配合面图像;
对所述装配配合面图像进行处理,得到装配特征的实测数据;
将所述装配特征的实测数据与装配特征的理论数据进行对比,得到实测装配偏差量。
3.根据权利要求1所述的一种装配粗差判定方法,其特征在于,所述利用顺序统计、秩估计和信息熵相结合的方法去除所述装配偏差仿真数据中的极限值,得到置信装配偏差仿真数据,具体包括:
将所述装配偏差仿真数据按从小到大的顺序排列,得到有序装配偏差仿真数据集合;
依次计算所述有序装配偏差仿真数据集合中每个数据的平均秩估计;
根据所述有序装配偏差仿真数据集合中每个数据的平均秩估计计算所述有序装配偏差仿真数据集合的信息熵;
利用所述有序装配偏差仿真数据集合的信息熵计算所述有序装配偏差仿真数据集合的不确定度;
根据所述有序装配偏差仿真数据集合的中所有数据的平均值和所述有序装配偏差仿真数据集合的不确定度确定所述有序装配偏差仿真数据集合的置信区间;
判断所述装配偏差仿真数据中是否存在位于所述置信装配偏差仿真数据区间之外的数据,得到第一判断结果;
若所述第一判断结果表示不存在,则确定所述装配偏差仿真数据为置信装配偏差仿真数据;
若所述第一判断结果表示存在,则将所述装配偏差仿真数据中位于所述置信装配偏差仿真数据区间之外的数据去除之后返回计算所述有序装配偏差仿真数据集合的步骤。
4.根据权利要求1所述的一种装配粗差判定方法,其特征在于,所述利用顺序统计、秩估计和信息熵结合的方法计算每个所述装配偏差仿真数据组的不确定度区间,得到装配偏差合格区间,具体包括:
将每个所述装配偏差仿真数据组中的数据按从小到大的顺序排列,得到有序装配偏差仿真数据组;
依次计算每个所述有序装配偏差仿真数据组中每个数据的平均秩估计;
根据每个所述有序装配偏差仿真数据组中每个数据的平均秩估计计算每个所述有序装配偏差仿真数据组的信息熵;
利用每个所述有序装配偏差仿真数据组的信息熵计算每个所述有序装配偏差仿真数据组的不确定度;
在每个所述有序装配偏差仿真数据组的不确定度中,筛选出最小值作为所述装配偏差合格区间的最小值,筛选出最大值作为所述装配偏差合格区间的最大值,得到所述装配偏差合格区间。
5.根据权利要求1所述的一种装配粗差判定方法,其特征在于,所述计算包含同一所述实测装配偏差量的所有融合数据的平均不确定度,具体包括:
计算每个所述融合数据集合的不确定度;
将包含同一所述实测装配偏差量的所有融合数据集合划分到同一个组中,得到多个融合数据组;
求取包含在同一个融合数据组中的所有融合数据集合的不确定度的平均值,得到各个融合数据组的平均不确定度。
6.一种装配粗差判定系统,其特征在于,包括:
获取模块,用于获取装配配合面的实测装配偏差量;
变动范围确定模块,用于利用先验经验和知识,分析装配偏差的变动范围;
仿真模块,用于利用蒙特卡洛法对所述装配偏 差的变动范围进行装配偏差有限元仿真,得到装配偏差仿真数据;
置信数据筛选模块,用于利用顺序统计、秩估计和信息熵相结合的方法去除所述装配偏差仿真数据中的极限值,得到置信装配偏差仿真数据;
分组模块,用于以所述实测装配偏差量的数量为一个数据组内的数据的数量,将所述置信装配偏差仿真数据进行分组,得到多个装配偏差仿真数据组;
合格区间确定模块,用于利用顺序统计、秩估计和信息熵结合的方法计算每个所述装配偏差仿真数据组的不确定度区间,得到装配偏差合格区间;
融合模块,用于将每个所述装配偏差仿真数据组和所述实测装配偏差量组成的实测数据集合融合,得到多个融合数据集合;每个所述融合数据集合均由一个装配偏差仿真数据组和实测数据集合相加形成;
平均不确定度计算模块,用于计算包含同一所述实测装配偏差量的所有融合数据的平均不确定度;
粗差判定模块,用于将所述平均不确定度与所述装配偏差合格区间对比,确定包含在所述装配偏差合格区间内的所述平均不确定度所对应的实测装配偏差量为合格偏差,确定位于所述装配偏差合格区间之外的所述平均不确定度所对应的实测装配偏差量为粗差。
7.根据权利要求6所述的一种装配粗差判定系统,其特征在于,所述获取模块,具体包括:
图像获取单元,用于获取数字化测量设备采集的装配配合面图像;
图像处理单元,用于对所述装配配合面图像进行处理,得到装配特征的实测数据;
偏差计算单元,用于将所述装配特征的实测数据与装配特征的理论数据进行对比,得到实测装配偏差量。
8.根据权利要求6所述的一种装配粗差判定系统,其特征在于,所述置信数据筛选模块,具体包括:
第一排列单元,用于将所述装配偏差仿真数据按从小到大的顺序排列,得到有序装配偏差仿真数据集合;
第一平均秩估计计算单元,用于依次计算所述有序装配偏差仿真数据集合中每个数据的平均秩估计;
第一信息熵计算单元,用于根据所述有序装配偏差仿真数据集合中每个数据的平均秩估计计算所述有序装配偏差仿真数据集合的信息熵;
第一不确定计算单元,用于利用所述有序装配偏差仿真数据集合的信息熵计算所述有序装配偏差仿真数据集合的不确定度;
置信区间确定单元,用于根据所述有序装配偏差仿真数据集合的中所有数据的平均值和所述有序装配偏差仿真数据集合的不确定度确定所述有序装配偏差仿真数据集合的置信区间;
第一判断单元,用于判断所述装配偏差仿真数据中是否存在位于所述置信装配偏差仿真数据区间之外的数据,得到第一判断结果;
置信数据确定单元,用于若所述第一判断结果表示不存在,则确定所述装配偏差仿真数据为置信装配偏差仿真数据;
返回单元,用于若所述第一判断结果表示存在,则将所述装配偏差仿真数据中位于所述置信装配偏差仿真数据区间之外的数据去除之后返回计算所述有序装配偏差仿真数据集合的步骤。
9.根据权利要求6所述的一种装配粗差判定系统,其特征在于,所述合格区间确定模块,具体包括:
第二排列单元,用于将每个所述装配偏差仿真数据组中的数据按从小到大的顺序排列,得到有序装配偏差仿真数据组;
第二平均秩估计计算单元,用于依次计算每个所述有序装配偏差仿真数据组中每个数据的平均秩估计;
第二信息熵计算单元,用于根据每个所述有序装配偏差仿真数据组中每个数据的平均秩估计计算每个所述有序装配偏差仿真数据组的信息熵;
第二不确定计算单元,用于利用每个所述有序装配偏差仿真数据组的信息熵计算每个所述有序装配偏差仿真数据组的不确定度;
合格区间确定单元,用于在每个所述有序装配偏差仿真数据组的不确定度中,筛选出最小值作为所述装配偏差合格区间的最小值,筛选出最大值作为所述装配偏差合格区间的最大值,得到所述装配偏差合格区间。
10.根据权利要求6所述的一种装配粗差判定系统,其特征在于,所述平均不确定度计算模块,具体包括:
第三不确定度计算单元,用于计算每个所述融合数据集合的不确定度;
分组单元,用于将包含同一所述实测装配偏差量的所有融合数据集合划分到同一个组中,得到多个融合数据组;
平均不确定度计算单元,用于求取包含在同一个融合数据组中的所有融合数据集合的不确定度的平均值,得到各个融合数据组的平均不确定度。
CN201810182694.3A 2018-03-06 2018-03-06 一种装配粗差判定方法及系统 Active CN108427837B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810182694.3A CN108427837B (zh) 2018-03-06 2018-03-06 一种装配粗差判定方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810182694.3A CN108427837B (zh) 2018-03-06 2018-03-06 一种装配粗差判定方法及系统

Publications (2)

Publication Number Publication Date
CN108427837A CN108427837A (zh) 2018-08-21
CN108427837B true CN108427837B (zh) 2021-09-10

Family

ID=63157951

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810182694.3A Active CN108427837B (zh) 2018-03-06 2018-03-06 一种装配粗差判定方法及系统

Country Status (1)

Country Link
CN (1) CN108427837B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109460599B (zh) * 2018-10-30 2022-06-28 南昌航空大学 一种装配特征偏差的传递量化分析方法及系统
CN114547920B (zh) * 2022-04-28 2022-08-12 南昌航空大学 一种飞机结构件装配偏差仿真模型可信度评价方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4734702A (en) * 1986-02-25 1988-03-29 Litton Systems, Inc. Passive ranging method and apparatus
CN102411735A (zh) * 2011-09-09 2012-04-11 河海大学常州校区 可重构装配系统重构规划方案评价方法
CN104376211A (zh) * 2014-11-14 2015-02-25 合肥工业大学 一种坐标测量机测量不确定度评定方法
CN106066901B (zh) * 2016-04-22 2019-05-07 中南大学 一种gnss自动化变形监测的基准点稳定性分析方法
CN106570252B (zh) * 2016-10-26 2019-05-24 中国运载火箭技术研究院 一种基于概率技术的热防护系统设计方法
CN107138586B (zh) * 2017-07-11 2018-10-19 南昌航空大学 一种用于密封多层管的推头结构
CN107609227B (zh) * 2017-08-23 2020-10-30 北京理工大学 一种基于最大熵理论的装配工艺优化方法

Also Published As

Publication number Publication date
CN108427837A (zh) 2018-08-21

Similar Documents

Publication Publication Date Title
CN107885928B (zh) 考虑测量误差的步进应力加速性能退化可靠性分析方法
TWI660277B (zh) 資訊處理裝置及資訊處理方法
CN103954940B (zh) 雷达网基于交叉定位点聚类的集中式与分布式压制干扰鉴别方法
CN108427837B (zh) 一种装配粗差判定方法及系统
CN110544047A (zh) 一种不良数据辨识方法
CN113807004A (zh) 基于数据挖掘的刀具寿命预测方法、装置及系统
CN116224930A (zh) 一种数控磨床产品的加工控制方法及系统
CN115422617A (zh) 一种基于cad的框架图像尺寸测量方法、装置及介质
CN115794916A (zh) 多源数据融合的数据处理方法、装置、设备和存储介质
CN109684713B (zh) 基于贝叶斯的复杂系统可靠性分析方法
García-Alfonso et al. Comparison of uncertainty analysis of the Montecarlo and Latin Hypercube algorithms in a camera calibration model
CN103902798A (zh) 数据预处理方法
JP2019528523A (ja) 計測−収率の相関性の分析方法及びシステム
CN109460599B (zh) 一种装配特征偏差的传递量化分析方法及系统
CN106779354B (zh) 用于飞行器性能评估的Bayes数据融合评估方法
CN115859211A (zh) 一种基于三温不确定度估算模型的地表温度产品融合方法
CN113722917B (zh) 基于垂直剖面上裂变径迹长度分布的热史模拟方法及系统
CN112328951B (zh) 一种分析样品的实验数据的处理方法
CN115752243A (zh) 一种测量数据融合方法
CN113065604B (zh) 一种基于dtw-dbscan算法的空中目标分群方法
RU2586025C2 (ru) Способ автоматической кластеризации объектов
Villa et al. IUGS–IUPAC recommendations and status reports on the half-lives of 87Rb, 146Sm, 147Sm, 234U, 235U, and 238U (IUPAC Technical Report)
CN113312587A (zh) 基于arima预测和回归预测的传感器采集数据缺失值处理方法
CN109146840B (zh) 基于高斯概率统计的数据几何定位精度评估方法及装置
US6766265B2 (en) Processing tester information by trellising in integrated circuit technology development

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant