CN107169681B - 误差分区表征方法 - Google Patents

误差分区表征方法 Download PDF

Info

Publication number
CN107169681B
CN107169681B CN201710515169.4A CN201710515169A CN107169681B CN 107169681 B CN107169681 B CN 107169681B CN 201710515169 A CN201710515169 A CN 201710515169A CN 107169681 B CN107169681 B CN 107169681B
Authority
CN
China
Prior art keywords
error
region
msub
equipotential lines
mrow
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
CN201710515169.4A
Other languages
English (en)
Other versions
CN107169681A (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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201710515169.4A priority Critical patent/CN107169681B/zh
Publication of CN107169681A publication Critical patent/CN107169681A/zh
Application granted granted Critical
Publication of CN107169681B publication Critical patent/CN107169681B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30Computing 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)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Quality & Reliability (AREA)
  • General Health & Medical Sciences (AREA)
  • Evolutionary Biology (AREA)
  • Manufacturing & Machinery (AREA)
  • Algebra (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Game Theory and Decision Science (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的误差区域判定,完成根据误差类型的曲面划分。
所述的对每个误差区域进行表征是指对每个误差区域使用相应的指标衡量其特征。其中,均匀误差区域以平均误差δ,区域面积Sj为指标。梯度误差区域以梯度误差最大值δTMax,梯度误差爬升加速度aT和区域面积ST为指标,局部误差区域以局部误差最大值δLMax,局部误差爬升加速度al和区域面积Sl为指标。其中,δ、δ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的误差区域判定,完成根据误差类型的曲面划分。
所述的对每个误差区域进行表征是对每个误差区域使用相应的指标衡量其特征;其中,均匀误差区域以平均误差δ,区域面积Sj为指标;梯度误差区域以梯度误差最大值δTMax,梯度误差爬升加速度aT和区域面积ST为指标,局部误差区域以局部误差最大值δLMax,局部误差爬升加速度al和区域面积Sl为指标;其中,δ、δ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为均匀误差区域,比较其中所有的误差点,得到该区域的平均误差δ为+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 (3)

1.一种误差分区表征方法,先通过测量获得当前加工对象的实际面并与理论面进行比对,再将其离散为基本单元并获得所有基本单元的误差,然后在理论面上根据所有基本单元的误差将理论面划分为合格区域和超差区域,再根据误差规律将超差区域进一步细分为多个误差区域,最后对每个误差区域进行表征;所述的基本单元为能将误差离散化的形式中的点单元、三角网格单元或四边形网格单元;合格区域为:区域内所有基本单元的误差均在给定公差范围内;超差区域为:区域内所有基本单元的误差均超出给定公差范围;所述的根据误差规律将超差区域进一步细分为多个误差区域是指使用误差等势图法将理论面根据误差规律划分为多个子面;所述的误差等势图为:以工件理论面为基准,将理论面向外法向和内法向两个方向每次偏移距离ΔE去截取实际面得到多组空间曲线,再将所有曲线投影到理论面上得到误差等势图;所述的ΔE是偏移距离,其值根据误差分布范围而定;同一条曲线上的所有基本单元的误差是相等的,由等势图的定义可知,在任一区域误差最大值的局部范围内,任意两条相邻曲线所代表误差的差值均为ΔE;并且,以一条曲线所代表的误差值来标识这条曲线,误差值为ΔE的等势线标识为ΔE等势线;其特征是所述的误差等势图法是根据误差等势图中的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的误差区域判定,完成根据误差类型的曲面划分。
2.根据权利要求1所述的方法,其特征是所述的将超差区域进一步细分为多个误差区域是指将超差区域的误差细分为均匀误差、梯度误差和局部误差;所述的均匀误差为:工件实际面在某一区域内所有基本单元的误差均接近某一值,呈现出均匀分布的特点,表现为实际面的整体偏移;所述的梯度误差为:工件实际面在某一区域内基本单元的误差朝某一方向不断增大,在接近该区域误差最大值的局部范围内呈现出加速增大的趋势,最终在实际面边界处达到该区域的梯度误差最大值δTMax,表现为实际面在这一边界处的翘曲;所述的局部误差为:工件实际面在某一区域内基本单元的误差向区域内部一点不断增大,并在接近该区域误差最大值的局部范围内呈现出减速增大的趋势,最终在区域内某一基本单元处达到该区域的局部误差最大值δLMax,表现为实际面上这一区域隆起形成凸包或塌陷形成凹坑。
3.根据权利要求1所述的方法,其特征是所述的对每个误差区域进行表征是对每个误差区域使用相应的指标衡量其特征;其中,均匀误差区域以平均误差δ,区域面积Sj为指标;梯度误差区域以梯度误差最大值δTMax,梯度误差爬升加速度aT和区域面积ST为指标,局部误差区域以局部误差最大值δLMax,局部误差爬升加速度al和区域面积Sl为指标;其中,δ、δTMax和δLMax可通过比较误差基本单元得到,δj、ST和Sl在分区域完成之后易求得,而aT和al可通过等势图中的误差密度变化情况来表示,误差密度为相邻两等势线间所代表的误差的差值去除以它们的最小法向间距;
<mrow> <msub> <mi>a</mi> <mi>T</mi> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mfrac> <msub> <mi>&amp;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>&amp;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> <mn>1</mn> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mfrac> <msub> <mi>&amp;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>&amp;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为相邻两个等势线所代表误差的差值。
CN201710515169.4A 2017-06-29 2017-06-29 误差分区表征方法 Active CN107169681B (zh)

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 CN107169681A (zh) 2017-09-15
CN107169681B true 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)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110110289B (zh) * 2019-05-24 2023-03-17 大连市市政设计研究院有限责任公司 一种降低土样分析超差率的方法

Citations (7)

* Cited by examiner, † Cited by third party
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
CN103256916A (zh) * 2013-06-10 2013-08-21 陈磊磊 一种基于最小区域的零件平面度误差评定方法
CN103292773A (zh) * 2013-06-18 2013-09-11 陈磊磊 一种基于最小区域的对称度误差评定方法
CN103292769A (zh) * 2013-06-19 2013-09-11 陈磊磊 一种基于最小区域的平面倾斜度误差评定方法
CN104550380A (zh) * 2014-12-11 2015-04-29 广东工业大学 一种基于水火加工平台的船体外板成型检测方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102880114B (zh) * 2012-10-18 2014-12-10 南京航空航天大学 飞机结构件数控加工在线自适应刀轨补偿方法
CN104021242B (zh) * 2014-05-15 2017-02-15 南京航空航天大学 基于零件特征的数控机床加工能力评定方法
CN105446264B (zh) * 2015-12-24 2018-05-22 南京航空航天大学 基于特征的机床精度优化设计方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
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 中国科学院沈阳计算技术研究所有限公司 一种五轴数控侧铣加工用刀具半径补偿方法
CN103256916A (zh) * 2013-06-10 2013-08-21 陈磊磊 一种基于最小区域的零件平面度误差评定方法
CN103292773A (zh) * 2013-06-18 2013-09-11 陈磊磊 一种基于最小区域的对称度误差评定方法
CN103292769A (zh) * 2013-06-19 2013-09-11 陈磊磊 一种基于最小区域的平面倾斜度误差评定方法
CN104550380A (zh) * 2014-12-11 2015-04-29 广东工业大学 一种基于水火加工平台的船体外板成型检测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"一种五轴机床检验试件轮廓误差的处理与显示技术研究";丁杰雄 等;《组合机床与自动化加工技术》;20121030(第10期);摘要、第2章 *
"数控加工误差主动补偿方法";周静 等;《计算机集成制造系统》;20100930;第16卷(第9期);第1903页左栏第10-21行、第1.2节 *

Also Published As

Publication number Publication date
CN107169681A (zh) 2017-09-15

Similar Documents

Publication Publication Date Title
CN106273446B (zh) 一种用于3d打印的切片路径生成方法及系统
CN107169681B (zh) 误差分区表征方法
CN111016179B (zh) 一种基于增材制造的变层厚剖分计算方法
CN110738092B (zh) 一种发票文本检测方法
Nayak et al. A modified triangle box-counting with precision in error fit
CN103823928A (zh) 基于扫描偏置的三维打印零件支撑自动生成方法
CN105772721A (zh) 基于增材制造的光斑大小可变的扫描路径生成方法
CN111462119B (zh) 一种基于机器视觉的宽厚板剪切排样方法
CN102073867A (zh) 一种遥感图像分类方法及装置
CN111097906B (zh) 基于多激光器的扫描分配方法、装置以及三维物体制造设备
CN102601227A (zh) 适用于笔记本电脑键盘框架的冲压模具及冲压方法
CN103927765B (zh) 用于定位图像中条码区域的方法和装置
JP2012019066A (ja) 荷電粒子ビーム描画方法及び装置
CN112200870B (zh) 基于孪生网络的分类和位置损失相结合的单目标跟踪方法
CN107256556A (zh) 一种基于灰度跳变思想的太阳能电池组件单元分割方法
US2559915A (en) Half-tone screen determiner
CN106790812A (zh) 移动终端中框料板
CN105373672A (zh) 一种四边形单元实时加密-减密方法
US11759891B2 (en) Method for depositing horizontal rib on unequal-height section by laser cladding
CN113128574A (zh) 场景缩减方法、装置及终端设备
CN113553690A (zh) 电路板的加工方法及加工设备、计算机可读存储介质
US9519960B2 (en) Method for estimating a number of objects
CN106547758A (zh) 一种数据分箱的方法和设备
CN106228126B (zh) 一种基于视频分析的影院上座人数统计方法、装置及电子设备
CN103310447A (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