CN107169681A - 误差分区表征方法 - Google Patents
误差分区表征方法 Download PDFInfo
- Publication number
- CN107169681A CN107169681A CN201710515169.4A CN201710515169A CN107169681A CN 107169681 A CN107169681 A CN 107169681A CN 201710515169 A CN201710515169 A CN 201710515169A CN 107169681 A CN107169681 A CN 107169681A
- Authority
- CN
- China
- Prior art keywords
- error
- area
- msub
- equipotential
- tolerance
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract 11
- 238000012512 characterization method Methods 0.000 claims abstract 3
- 238000010586 diagram Methods 0.000 claims 6
- 230000001133 acceleration Effects 0.000 claims 2
- 230000009194 climbing Effects 0.000 claims 2
- 238000013459 approach Methods 0.000 claims 1
- 238000009826 distribution Methods 0.000 claims 1
- 238000005259 measurement Methods 0.000 claims 1
- 238000009827 uniform distribution Methods 0.000 claims 1
- 238000013316 zoning Methods 0.000 claims 1
- 238000004519 manufacturing process Methods 0.000 abstract 1
- 238000005192 partition Methods 0.000 abstract 1
- 238000003908 quality control method Methods 0.000 abstract 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/04—Manufacturing
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- Data Mining & Analysis (AREA)
- Economics (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Physics (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Pure & Applied Mathematics (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Entrepreneurship & Innovation (AREA)
- Mathematical Optimization (AREA)
- Operations Research (AREA)
- Tourism & Hospitality (AREA)
- Probability & Statistics with Applications (AREA)
- General Health & Medical Sciences (AREA)
- Bioinformatics & Computational Biology (AREA)
- Quality & Reliability (AREA)
- Health & Medical Sciences (AREA)
- Evolutionary Biology (AREA)
- Manufacturing & Machinery (AREA)
- Game Theory and Decision Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Algebra (AREA)
- Life Sciences & Earth Sciences (AREA)
- Primary Health Care (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- General Factory Administration (AREA)
Abstract
本发明公开了一种针对误差的分区表征方法,该方法首先通过测量获得当前加工对象的实际面并与理论面进行比对,再将其离散为基本单元并获得所有基本单元的误差,然后根据所有基本单元的误差规律进行区域划分,最后对每个误差区域进行表征。和原有的误差表征方法相比,本发明能将误差表征得更加清楚,为制造过程中的质量控制以及误差溯源提供了更加详细的依据。
Description
技术领域
本发明涉及一种误差表征方法,具体地说是在产品制造中根据误差分布规律划分区域,再使用相应的指标去衡量,从而清楚地表征误差区域性特征的方法。
背景技术
在产品制造中,工艺水平、设备性能、工件材料等因素均会导致工件产生误差。目前,描述工件相关形状精度可以使用如圆度、圆柱度、轮廓度等形位公差,但是它们是基于整体工件而言的指标,只能用于检验产品是否合格,无法表达局部的误差情况。对于复杂曲面工件而言,如飞机蒙皮,实际面可能包含局部误差、梯度误差、均匀误差在内的多种误差区域,形位公差无法表达误差的区域性特征。所以,对于制造过程中的质量控制或误差溯源而言,形位公差所提供的信息量还不够。
针对获得复杂曲面工件误差的问题,论文《基于自适应采样的曲面加工误差在机测量方法》(机器仪表学报2016年第1期)获取加工曲面上少量测点的坐标数据,拟合加工曲面再计算理论面和实际面的法向距离,获得了自由曲面的详细误差并绘制出三维网格图。这种方法虽然提供了非常详细的误差信息,但是这些误差信息没有经过整理和分类,难以使用。
因此,需要一种经过整理分类,容易使用并且能对产品质量控制,误差溯源提供足够信息量的误差表征方法,可以清楚的表征误差区域性特征。
发明内容
本发明的目的是针对已有的误差表征方法无法清楚地表征区域性分布特征的问题,发明了一种根据误差分布规律划分区域,再使用相应的指标去衡量,从而清楚地表征误差区域性特征的误差分区表征方法。
本发明的技术方案是:
一种误差分区表征方法,其特征在于先通过测量获得当前加工对象的实际面并与理论面进行比对,再将其离散为基本单元并获得所有基本单元的误差,然后在理论面上根据所有基本单元的误差将理论面划分为合格区域和超差区域,再根据误差规律将超差区域进一步细分为多个误差区域,最后对每个误差区域进行表征。所述的基本单元可以是点、三角网格、四边形网格等能将误差离散化的形式中的单元。合格区域为:区域内所有基本单元的误差均在给定公差范围内。超差区域为:区域内所有基本单元的误差均超出给定公差范围。
所述的根据误差规律将超差区域进一步细分为多个误差区域是指使用误差等势图法将理论面根据误差规律划分为多个子面。误差等势图为:以工件理论面为基准,将理论面向外法向和内法向两个方向每次偏移距离ds去截取实际面得到多组空间曲线,再将所有曲线投影到理论面上得到误差等势图。所述的ΔE是偏移距离,其值根据误差分布范围而定。同一条曲线上的所有基本单元的误差是相等的,由等势图的定义可知,在区域误差最大值的局部范围内,任意两条相邻曲线所代表误差的差值均为ΔE。并且,以一条曲线所代表的误差值来标识这条曲线,例如误差值为ΔE的等势线标识为ΔE等势线。
所述的多个误差区域的误差特征包括均匀误差、梯度误差和局部误差。均匀误差为:工件实际面在某一区域内所有基本单元的误差均接近于某一值,呈现出均匀分布的特点,表现为实际面的整体偏移。所述的梯度误差为:工件实际面在某一区域内基本单元的误差朝某一方向不断增大,在接近该区域误差最大值的局部范围内呈现出加速增大的趋势,最终在实际面边界处达到该区域的梯度误差最大值δTMax,表现为实际面在这一边界处的翘曲。所述的局部误差为:工件实际面在某一区域内基本单元的误差向区域内部一点不断增大,并在接近该区域误差最大值的局部范围内呈现出减速增大的趋势,最终在区域内某一基本单元处达到该区域的局部误差最大值δLMax,表现为实际面上这一区域隆起形成凸包或塌陷形成凹坑。
所述的误差等势图法,其实施方法为:先根据误差等势图中的Tu,Td等势线划分出合格区域,再在剩下的区域中,根据误差近乎均匀不变的特点划分出均匀误差区域,根据误差在该区域误差最大值的局部范围内加速上升的特点划分出梯度误差区域,根据误差在该区域误差最大值的局部范围内减速上升的特点划分出局部误差区域。所述的Tu为公差上限,Td为公差下限,公差范围为[Td,Tu]。具体步骤为:
步骤一、在误差等势图中获得误差为Tu,Td的所有等势线。
步骤二、以这些等势线为边界将曲面划分为k个区域A1,A2,A3…Ai…Ak。
步骤三、对于其中的任意一个区域Ai而言,若其中不包含误差超出[Td,Tu]范围的等势线,则Ai是合格区域。若包含误差超出[Td,Tu]范围的等势线,则Ai是超差区域。用这种方法得到所有的超差区域。设有m个超差区域,m≤k,则所有的超差区域记为B1,B2,B3…Bm。
步骤四、对于任意一个超差区域Bi而言,从误差最大的等势线开始,连续获得相邻的三个最小法向间距记为dmin1,dmin2和dmin3。其中,最小法向间距为相邻的两条等势线之间的最小法向间距。若dmin1>∈则认为等势线近似均匀不变,区域Bi为均匀误差区域。∈根据实际工件大小而定,当两条等势线的最小法向间距大于∈的时候就认为,等势线足够稀疏,误差值近似均匀不变。若dmin1<dmin2<dmin3且dmin3<∈,则表明等势线在区域误差最大值的局部范围内越来越密集,误差值加速上升,区域Bi为梯度误差区域。若dmin1>dmin2>dmin3且dmin1<∈则表明,等势线在区域误差最大值的局部范围内越来越稀疏,误差值减速上升,区域Bi为局部误差区域。
步骤五、迭代步骤四的过程直至完成B1,B2,B3…Bi…Bm的误差区域判定,完成根据误差类型的曲面划分。
步骤五、迭代步骤四的过程直至完成B1,B2,B3…Bi…Bm的误差区域判定,完成根据误差类型的曲面划分。
所述的对每个误差区域进行表征是指对每个误差区域使用相应的指标衡量其特征。其中,均匀误差区域以平均误差δjμ,区域面积Sj为指标。梯度误差区域以梯度误差最大值δTMax,梯度误差爬升加速度aT和区域面积ST为指标,局部误差区域以局部误差最大值δLMax,局部误差爬升加速度al和区域面积Sl为指标。其中,δjμ、δTMax和δLMax可通过比较误差基本单元得到,aT和al可通过等势图中的误差密度变化情况来表示,误差密度为相邻两等势线间所代表的误差的差值去除以它们的最小法向间距。
其中,dmin1,dmin2和dmin3为任意一个超差区域中,从误差最大的等势线开始,连续获得三个相邻的最小法向间距。ΔE为相邻两个等势线所代表误差的差值。
所述的误差可以是机械加工、零件成形等一切产品制造方法中产生的工件实际面和理论面的偏移,针对的工件可以是复杂曲面工件、简单曲面工件、平面工件等工件。
本发明的有益效果是:
1、本发明使用根据误差分布规律划分区域的方法将复杂的误差比对结果信息进行了整理和分类,使之容易使用的同时又保证了其足够的信息量。
2、本发明使用的误差等势图法可以高效地实现根据误差规律划分区域的功能。
3、本发明使用的误差区域表征指标准确地提供了表征特定误差区域的误差分布情况所需要的主要信息。
附图说明
图1为本发明实现功能的简单示意图。图中:1.1代表实际曲面,1.2代表理论曲面,1.3代表梯度误差区域,1.4代表合格区域,1.5代表局部误差区域。
图2为具体验证的飞机蒙皮零件示意图。图中:2.1代表理论面,2.2代表实际面,L和D是两个方向上的尺寸。
图3为某一蒙皮零件的误差等势图。图中:Td代表误差值为公差下限的等势线,Tu代表误差值为公差上限的等势线,δmax1代表第一个误差区域中误差值最大的等势线,δmax1-1是与δmax1相邻的等势线。δmax2代表第二个误差区域中误差值最大的等势线,从δmax2开始向外统计等势线,依次为δmax2-1,δmax2-2,δmax2-3。δmax3代表第三个误差区域中误差值最大的等势线,从δmax3开始向外统计等势线,依次为δmax3-1,δmax3-2,δmax3-3和δmax3-4。
图4为初步划分区域的结果,其中Td代表误差值为公差下限的等势线,Tu代表误差值为公差上限的等势线,A1~A7为划分的七个区域。
图5为划分出超差区域的结果,其中Td代表误差值为公差下限的等势线,Tu代表误差值为公差上限的等势线,B1~B3为划分的三个超差区域,阴影部分为合格区域。
图6为将B1~B3三个超差区域放大的示意图。其中,等势线的符号和之前是一样的,其中dmin1,dmin2和dmin3为三个相邻的最小法向间距,在图7~9中会详细解释。∈是一个用于判定均匀误差的指标。
图7为细分的第一个误差区域的详细等势图,其中,δmax1代表这个区域中误差值最大的等势线,δmax1-1是与δmax1相邻的等势线。dmin1为δmax1和δmax1-1的最小法向间距。∈是一个用于判定均匀误差的指标。
图8为细分的第二个误差区域的详细等势图,其中,δmax2代表这个区域中误差值最大的等势线,从δmax2开始向外统计等势线,依次为δmax2-1,δmax2-2,δmax2-3。dmin1为δmax2和δmax2-1的最小法向间距,dmin2为δmax2-1和δmax2-2的最小法向间距,dmin3为δmax2-2和δmax2-3的最小法向间距,∈是一个用于判定均匀误差的指标。
图9是细分的第三个误差区域等势图,其中,δmax3代表第这个区域中误差值最大的等势线,从δmax3开始向外统计等势线,依次为δmax3-1,δmax3-2,δmax3-3和δmax3-4。dmin1为δmax3和δmax3-1的最小法向间距,dmin2为δmax3-1和δmax3-2的最小法向间距,dmin3为δmax3-2和δmax3-3的最小法向间距,∈是一个用于判定均匀误差的指标。
具体实施方式
下面结合附图和示例零件对本发明作进一步说明。
一种误差分区表征方法,具体步骤是:
首先通过测量获得当前加工对象的实际面并与理论面进行比对,再将其离散为基本单元并获得所有基本单元的误差;
其次,在理论面上根据所有基本单元的误差将理论面划分为合格区域和超差区域;
第三,根据误差规律将超差区域进一步细分为多个误差区域,最后对每个误差区域进行表征;所述的基本单元为能将误差离散化的形式中的点单元、三角网格单元或四边形网格单元;合格区域为:区域内所有基本单元的误差均在给定公差范围内;超差区域为:区域内所有基本单元的误差均超出给定公差范围。其中:根据误差规律将超差区域进一步细分为多个误差区域是指使用误差等势图法将理论面根据误差规律划分为多个子面;所述的误差等势图为:以工件理论面为基准,将理论面向外法向和内法向两个方向每次偏移距离ΔE去截取实际面得到多组空间曲线,再将所有曲线投影到理论面上得到误差等势图;所述的ΔE是偏移距离,其值根据误差分布范围而定;同一条曲线上的所有基本单元的误差是相等的,由等势图的定义可知,在任一区域误差最大值的局部范围内,任意两条相邻曲线所代表误差的差值均为ΔE;并且,以一条曲线所代表的误差值来标识这条曲线,误差值为ΔE的等势线标识为ΔE等势线。将超差区域进一步细分为多个误差区域是指将超差区域的误差细分为均匀误差、梯度误差和局部误差;所述的均匀误差为:工件实际面在某一区域内所有基本单元的误差均接近于某一值,呈现出均匀分布的特点,表现为实际面的整体偏移;所述的梯度误差为:工件实际面在某一区域内基本单元的误差朝某一方向不断增大,在接近该区域误差最大值的局部范围内呈现出加速增大的趋势,最终在实际面边界处达到该区域的梯度误差最大值δTMax,表现为实际面在这一边界处的翘曲;所述的局部误差为:工件实际面在某一区域内基本单元的误差向区域内部一点不断增大,并在接近该区域误差最大值的局部范围内呈现出减速增大的趋势,最终在区域内某一基本单元处达到该区域的局部误差最大值δLMax,表现为实际面上这一区域隆起形成凸包或塌陷形成凹坑。所述的误差等势图法是根据误差等势图中的Tu,Td等势线划分出合格区域,再在剩下的区域中,根据误差近乎均匀不变的特点划分出均匀误差区域,根据误差在该区域最大值的局部范围内加速上升的特点划分出梯度误差区域,根据误差在该区域最大值的局部范围内减速上升的特点划分出局部误差区域;所述的Tu为公差上限,Td为公差下限,公差范围为[Td,Tu],具体步骤为:
步骤一、在误差等势图中获得误差为Tu,Td的所有等势线;
步骤二、以这些等势线为边界将曲面划分为k个区域A1,A2,A3…Ai…Ak;
步骤三、对于其中的任意一个区域Ai而言,若其中不包含误差超出[Td,Tu]范围的等势线,则Ai是合格区域;若包含误差超出[Td,Tu]范围的等势线,则Ai是超差区域;用这种方法得到所有的超差区域;设有m个超差区域,m≤k,则所有的超差区域记为B1,B2,B3…Bm;
步骤四、对于任意一个超差区域Bi而言,从误差最大的等势线开始,连续获得相邻的三个最小法向间距记为dmin1,dmin2和dmin3;其中,最小法向间距为相邻的两条等势线之间的最小法向间距;若dmin1>∈则认为等势线近似均匀不变,区域Bi为均匀误差区域,∈根据实际工件大小而定,当两条等势线的最小法向间距大于∈的时候就认为等势线足够稀疏,误差值近似均匀不变;若dmin1<dmin2<dmin3且dmin3<∈,则表明等势线在区域误差最大值的局部范围内越来越密集,误差值加速上升,区域Bi为梯度误差区域;若dmin1>dmin2>dmin3且dmin1<∈则表明,等势线在区域误差最大值的局部范围内越来越稀疏,误差值减速上升,区域Bi为局部误差区域;
步骤五、迭代步骤四的过程直至完成B1,B2,B3…Bi…Bm的误差区域判定,完成根据误差类型的曲面划分。
所述的对每个误差区域进行表征是对每个误差区域使用相应的指标衡量其特征;其中,均匀误差区域以平均误差δjμ,区域面积Sj为指标;梯度误差区域以梯度误差最大值δTMax,梯度误差爬升加速度aT和区域面积ST为指标,局部误差区域以局部误差最大值δLMax,局部误差爬升加速度al和区域面积Sl为指标;其中,δjμ、δTMax和δLMax可通过比较误差基本单元得到,Sj、ST和Sl在分区域完成之后易求得,而aT和al可通过等势图中的误差密度变化情况来表示,误差密度为相邻两等势线间所代表的误差的差值去除以它们的最小法向间距;
其中,dmin1,dmin2和dmin3为任意一个超差区域中,从误差最大的等势线开始,连续获得三个相邻的最小法向间距,ΔE为相邻两个等势线所代表误差的差值。
详述如下:
如图2,这是一张飞机蒙皮零件,其尺寸为4000mm×2000mm。其中,实际面(2.2)和理论面(2.1)有多种类型的误差,使用误差等势图法将理论面根据误差规律划分区域,具体步骤如下:
步骤一:绘制误差等势图。
先测量获得这张蒙皮的实际面2.2并与理论面2.1进行比对,再将其离散为大量的点,并获得所有点的误差。对于这张蒙皮零件而言,如果一个点的误差在[-0.2,+0.2]范围内,则这个点是合格的。
以理论面2.1为基准,将理论面向外法向方向分别偏移-0.2mm,-1.2mm,-2.2mm,-3.2mm……-10.2mm;+0.2mm,+1.2mm,+2.2mm,+3.2mm,+4.2mm,+5.2mm……+10.2mm。再使用偏移出来的面分别去截取实际面2.2,得到空间中多组曲线图,再将其全部投影到理论面上,得到如附图3所示等势图,图中所有数据如下表:
步骤二、根据误差大小将理论面划分为合格区域和超差区域。
如图3所示,获得误差为Tu,Td的所有等势线,此时的Tu为+0.2mm,Td为-0.2mm,以这些等势线为边界将曲面划分为了7个区域,记为A1,A2,A3…A7,得到如附图4所示的图。
对照附图3,将A1,A2,A3…A7一一对应,只有A3,A6和A7内部存在误差值超出[-0.2,+0.2]范围的等势线。将这三个区域重新记为B1,B2,B3,得到如附图5所示的图。其中B1,B2和B3就是超差区域,阴影部分为合格区域。
步骤三、判定所有误差区域的误差类型。
对照附图3和附图4,将B1,B2和B3所在区域的详细等势线图放大分析,过程如附图6。
如附图7所示,δmax1是这个误差区域中误差值最大的等势线,其误差值为-2.2mm,而它和与之相邻的等势线δmax1-1之间的最小法向间距dmin1为352.8mm。根据这张蒙皮的尺寸4000mm×2000mm以及其以往的误差范围,将∈定为45mm,dmin1>∈,所以B1是均匀误差区域。
如图8所示,δmax2是这个区域中误差最大的等势线,其误差值为+4.2mm,从δmax2开始向外统计等势线,依次为δmax2-1,δmax2-2和δmax2-3,如图所示,相邻之间的最小法向间距为,dmin1=84.6,δmax2=67.2mm,dmin3=50.4mm。∈>dmin1>δmax2>dmin3,表明在靠近最大值的局部范围内,等势线越来越稀疏,变化率越来越小,减速上升形成凸包,满足局部误差区域判定准则,B2为局部误差区域。
如图9所示,δmax3是这个区域中误差最大的等势线,其误差值为+5.2mm,从δmax3开始向统计等势线,依次为δmax3-1,δmax3-2和δmax3-3,如图所示,相邻之间的最小法向间距为:dmin1=50.4mm,δmax2=92.4mm,dmin3=117.6mm。dmin1<δmax2<dmin3<∈,表明在靠近最大值的局部范围内,等势线越来越密集,变化率越来越大,误差加速增大,在边界形成翘曲,满足梯度误差区域判定准则,B3为梯度误差区域。
步骤四、表征所有误差区域。
B1为均匀误差区域,比较其中所有的误差点,得到该区域的平均误差δjμ为+1.8mm,该区域面积Sj为9525.56cm2;B2为局部误差区域,比较其中所有的误差点,得到该区域的局部误差最大值δLMax为+4.3mm,该区域面积Sl为4386.86cm2,局部误差爬升加速度al采用下式计算。
计算出al为-1.19×10-4mm/mm2,其中加速度单位mm/mm2的含义为:每毫米变化的误差密度值。;B3为梯度误差区域,比较其中所有的误差点,得到该区域的梯度误差最大值δTMax为+5.8mm,该区域面积ST为5715.36cm2,梯度误差爬升加速度aT采用下式计算
计算出aT为1.23×10-4
最后,统计误差分区表征结果如下表
误差分区表征结果统计表
本发明未涉及部分与现有技术相同采用现有技术加以实现。
Claims (5)
1.误差分区表征方法,其特征在于先通过测量获得当前加工对象的实际面并与理论面进行比对,再将其离散为基本单元并获得所有基本单元的误差,然后在理论面上根据所有基本单元的误差将理论面划分为合格区域和超差区域,再根据误差规律将超差区域进一步细分为多个误差区域,最后对每个误差区域进行表征;所述的基本单元为能将误差离散化的形式中的点单元、三角网格单元或四边形网格单元;合格区域为:区域内所有基本单元的误差均在给定公差范围内;超差区域为:区域内所有基本单元的误差均超出给定公差范围。
2.根据权利要求1的方法,其特征是所述的根据误差规律将超差区域进一步细分为多个误差区域是指使用误差等势图法将理论面根据误差规律划分为多个子面;所述的误差等势图为:以工件理论面为基准,将理论面向外法向和内法向两个方向每次偏移距离ΔE去截取实际面得到多组空间曲线,再将所有曲线投影到理论面上得到误差等势图;所述的ΔE是偏移距离,其值根据误差分布范围而定;同一条曲线上的所有基本单元的误差是相等的,由等势图的定义可知,在任一区域误差最大值的局部范围内,任意两条相邻曲线所代表误差的差值均为ΔE;并且,以一条曲线所代表的误差值来标识这条曲线,误差值为ΔE的等势线标识为ΔE等势线。
3.根据权利要求1所述的方法,其特征是所述的将超差区域进一步细分为多个误差区域是指将超差区域的误差细分为均匀误差、梯度误差和局部误差;所述的均匀误差为:工件实际面在某一区域内所有基本单元的误差均接近于某一值,呈现出均匀分布的特点,表现为实际面的整体偏移;所述的梯度误差为:工件实际面在某一区域内基本单元的误差朝某一方向不断增大,在接近该区域误差最大值的局部范围内呈现出加速增大的趋势,最终在实际面边界处达到该区域的梯度误差最大值δTMax,表现为实际面在这一边界处的翘曲;所述的局部误差为:工件实际面在某一区域内基本单元的误差向区域内部一点不断增大,并在接近该区域误差最大值的局部范围内呈现出减速增大的趋势,最终在区域内某一基本单元处达到该区域的局部误差最大值δLMax,表现为实际面上这一区域隆起形成凸包或塌陷形成凹坑。
4.根据权利要求2所述的方法,其特征是所述的误差等势图法是根据误差等势图中的Tu,Td等势线划分出合格区域,再在剩下的区域中,根据误差近乎均匀不变的特点划分出均匀误差区域,根据误差在该区域最大值的局部范围内加速上升的特点划分出梯度误差区域,根据误差在该区域最大值的局部范围内减速上升的特点划分出局部误差区域;所述的Tu为公差上限,Td为公差下限,公差范围为[Td,Tu],具体步骤为:
步骤一、在误差等势图中获得误差为Tu,Td的所有等势线;
步骤二、以这些等势线为边界将曲面划分为k个区域A1,A2,A3…Ai…Ak;
步骤三、对于其中的任意一个区域Ai而言,若其中不包含误差超出[Td,Tu]范围的等势线,则Ai是合格区域;若包含误差超出[Td,Tu]范围的等势线,则Ai是超差区域;用这种方法得到所有的超差区域;设有m个超差区域,m≤k,则所有的超差区域记为B1,B2,B3…Bm;
步骤四、对于任意一个超差区域Bi而言,从误差最大的等势线开始,连续获得相邻的三个最小法向间距记为dmin1,dmin2和dmin3;其中,最小法向间距为相邻的两条等势线之间的最小法向间距;若dmin1>∈则认为等势线近似均匀不变,区域Bi为均匀误差区域,∈是一个数,根据实际工件大小而定,当两条等势线的最小法向间距大于∈的时候就认为等势线足够稀疏,误差值近似均匀不变;若dmin1<dmin2<dmin3且dmin3<∈,则表明等势线在区域误差最大值的局部范围内越来越密集,误差值加速上升,区域Bi为梯度误差区域;若dmin1>dmin2>dmin3且dmin1<∈则表明,等势线在区域误差最大值的局部范围内越来越稀疏,误差值减速上升,区域Bi为局部误差区域;
步骤五、迭代步骤四的过程直至完成B1,B2,B3…Bi…Bm的误差区域判定,完成根据误差类型的曲面划分。
5.根据权利要求1所述的方法,其特征是所述的对每个误差区域进行表征是对每个误差区域使用相应的指标衡量其特征;其中,均匀误差区域以平均误差δjμ,区域面积Sj为指标;梯度误差区域以梯度误差最大值δTMax,梯度误差爬升加速度aT和区域面积ST为指标,局部误差区域以局部误差最大值δLMax,局部误差爬升加速度al和区域面积Sl为指标;其中,δjμ、δTMax和δLMax可通过比较误差基本单元得到,Sj、ST和Sl在分区域完成之后易求得,而aT和al可通过等势图中的误差密度变化情况来表示,误差密度为相邻两等势线间所代表的误差的差值去除以它们的最小法向间距;
<mrow>
<msub>
<mi>a</mi>
<mi>T</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>&Delta;</mi>
<mi>E</mi>
</msub>
<msub>
<mi>d</mi>
<mrow>
<mi>min</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>-</mo>
<mfrac>
<msub>
<mi>&Delta;</mi>
<mi>E</mi>
</msub>
<msub>
<mi>d</mi>
<mrow>
<mi>min</mi>
<mn>3</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>/</mo>
<msub>
<mi>d</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<mn>2</mn>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>a</mi>
<mi>l</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>&Delta;</mi>
<mi>E</mi>
</msub>
<msub>
<mi>d</mi>
<mrow>
<mi>min</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>-</mo>
<mfrac>
<msub>
<mi>&Delta;</mi>
<mi>E</mi>
</msub>
<msub>
<mi>d</mi>
<mrow>
<mi>min</mi>
<mn>3</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>/</mo>
<msub>
<mi>d</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<mn>2</mn>
</mrow>
</msub>
</mrow>
其中,dmin1,dmin2和dmin3为任意一个超差区域中,从误差最大的等势线开始,连续获得三个相邻的最小法向间距,ΔE为相邻两个等势线所代表误差的差值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710515169.4A CN107169681B (zh) | 2017-06-29 | 2017-06-29 | 误差分区表征方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710515169.4A CN107169681B (zh) | 2017-06-29 | 2017-06-29 | 误差分区表征方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107169681A true CN107169681A (zh) | 2017-09-15 |
CN107169681B CN107169681B (zh) | 2018-01-16 |
Family
ID=59827152
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710515169.4A Active CN107169681B (zh) | 2017-06-29 | 2017-06-29 | 误差分区表征方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107169681B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110110289A (zh) * | 2019-05-24 | 2019-08-09 | 大连市市政设计研究院有限责任公司 | 一种降低土样分析超差率的方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6256546B1 (en) * | 1998-09-28 | 2001-07-03 | General Electric Company | System and method for numerical control processing of an in-processing part |
CN102371504A (zh) * | 2010-08-18 | 2012-03-14 | 中国科学院沈阳计算技术研究所有限公司 | 一种五轴数控侧铣加工用刀具半径补偿方法 |
EP1476797B1 (en) * | 2001-10-11 | 2012-11-21 | Laser Projection Technologies, Inc. | Method and system for visualizing surface errors |
CN102880114A (zh) * | 2012-10-18 | 2013-01-16 | 南京航空航天大学 | 飞机结构件数控加工在线自适应刀轨补偿方法 |
CN103256916A (zh) * | 2013-06-10 | 2013-08-21 | 陈磊磊 | 一种基于最小区域的零件平面度误差评定方法 |
CN103292773A (zh) * | 2013-06-18 | 2013-09-11 | 陈磊磊 | 一种基于最小区域的对称度误差评定方法 |
CN103292769A (zh) * | 2013-06-19 | 2013-09-11 | 陈磊磊 | 一种基于最小区域的平面倾斜度误差评定方法 |
CN104021242A (zh) * | 2014-05-15 | 2014-09-03 | 南京航空航天大学 | 基于零件特征的数控机床加工能力评定方法 |
CN104550380A (zh) * | 2014-12-11 | 2015-04-29 | 广东工业大学 | 一种基于水火加工平台的船体外板成型检测方法 |
CN105446264A (zh) * | 2015-12-24 | 2016-03-30 | 南京航空航天大学 | 基于特征的机床精度优化设计方法 |
-
2017
- 2017-06-29 CN CN201710515169.4A patent/CN107169681B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6256546B1 (en) * | 1998-09-28 | 2001-07-03 | General Electric Company | System and method for numerical control processing of an in-processing part |
EP1476797B1 (en) * | 2001-10-11 | 2012-11-21 | Laser Projection Technologies, Inc. | Method and system for visualizing surface errors |
CN102371504A (zh) * | 2010-08-18 | 2012-03-14 | 中国科学院沈阳计算技术研究所有限公司 | 一种五轴数控侧铣加工用刀具半径补偿方法 |
CN102880114A (zh) * | 2012-10-18 | 2013-01-16 | 南京航空航天大学 | 飞机结构件数控加工在线自适应刀轨补偿方法 |
CN103256916A (zh) * | 2013-06-10 | 2013-08-21 | 陈磊磊 | 一种基于最小区域的零件平面度误差评定方法 |
CN103292773A (zh) * | 2013-06-18 | 2013-09-11 | 陈磊磊 | 一种基于最小区域的对称度误差评定方法 |
CN103292769A (zh) * | 2013-06-19 | 2013-09-11 | 陈磊磊 | 一种基于最小区域的平面倾斜度误差评定方法 |
CN104021242A (zh) * | 2014-05-15 | 2014-09-03 | 南京航空航天大学 | 基于零件特征的数控机床加工能力评定方法 |
CN104550380A (zh) * | 2014-12-11 | 2015-04-29 | 广东工业大学 | 一种基于水火加工平台的船体外板成型检测方法 |
CN105446264A (zh) * | 2015-12-24 | 2016-03-30 | 南京航空航天大学 | 基于特征的机床精度优化设计方法 |
Non-Patent Citations (2)
Title |
---|
丁杰雄 等: ""一种五轴机床检验试件轮廓误差的处理与显示技术研究"", 《组合机床与自动化加工技术》 * |
周静 等: ""数控加工误差主动补偿方法"", 《计算机集成制造系统》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110110289A (zh) * | 2019-05-24 | 2019-08-09 | 大连市市政设计研究院有限责任公司 | 一种降低土样分析超差率的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107169681B (zh) | 2018-01-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101702200B (zh) | 一种机载激光雷达点云数据的自动分类方法 | |
GB2588013A (en) | Operation of a vehicle using multiple motion constraints | |
US20200271502A1 (en) | Method and device for acquiring volume of structure, non-transitory computer-readable storage medium and printer | |
CN110895333A (zh) | 一种基于多普勒频率的77g车载雷达数据快速聚类方法 | |
CN109000656B (zh) | 基于空间聚类的水下地形匹配导航适配区选择方法 | |
Hassan et al. | Benchmarking of pattern recognition techniques for online tool wear detection | |
CN107169681A (zh) | 误差分区表征方法 | |
CN109597354A (zh) | 一种三角网格模型的多约束数控加工轨迹生成方法 | |
Wang et al. | Artificial Intelligence for advanced non-conventional machining processes | |
CN109901512A (zh) | 一种基于加工参数标准化的车削工时定额方法 | |
CN106815447B (zh) | 基于历史数据的复杂结构件加工特征智能定义与分类方法 | |
CN114116829A (zh) | 异常数据分析方法、异常数据分析系统和存储介质 | |
CN103218837B (zh) | 一种基于经验分布函数的不等组距直方图的绘制方法 | |
CN107122475A (zh) | 大数据异常点检测方法及其系统 | |
CN108537249B (zh) | 一种密度峰值聚类的工业过程数据聚类方法 | |
CN112116590B (zh) | 一种复杂曲面零件超差区域加工边界识别方法 | |
CN112750119A (zh) | 一种面向白玻璃盖板表面微弱缺陷的检测与测量方法 | |
CN105930602A (zh) | 一种基于最优加权小波包熵的颤振检测方法 | |
CN111061246B (zh) | 一种机械产品装配过程故障模式、影响及危害性分析方法 | |
CN110866317B (zh) | 基于产品历史信息及设计行为分析的设计工具包重构方法 | |
CN110102507A (zh) | 珍珠分级方法 | |
CN104715137A (zh) | 一种考虑犹豫信息权重获取的决策方法 | |
EP4068026A1 (en) | Data novelty detection method and apparatus | |
Latiff et al. | Dimensional accuracy evaluation of rapid prototyping fused deposition modeling process of FDM200mc machine on basic engineering profiles | |
CN113902711A (zh) | 待测物检出缺陷的聚类方法、装置、设备及存储介质 |
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 |