平面度误差的测量及数据处理

1、-. z. -. -可修编- .学校代码： *学 号： *课程设计说明书题目: 平面度误差的测量及数据处理学生：学院：班级：指导教师：摘要平面度误差是将被测实际外表与理想平面进展比拟，两者之间的线值距离即为平面度误差值；或通过测量实际外表上假设干点的相对高度差，再换算以线值表示的平面度误差值。本文就平面度误差的数学模型与按最小二乘法建立理想平面评定基准的数学模型展开分析讨论；并结合案列分析，得出比拟客观地评定平面度误差或者测量较大平面的平面度误差，最小二乘法是最正确方法1关键词：最小二乘法；平面度误差；最正确方法AbstractFlatness error is measured by act

2、ual surface with an ideal plane are pared, the line between the two values of distance, which is the flatness error values, or by measuring the actual surface on several points of relative height difference, conversion to line value representation of flatness error value .This paper studies the math

3、enatical model of flatness error and ideal plane made by least square method .With illust ration of practical cases, the author reaches the conclusion that least sauare method is the best one in judging and measuring larger planes flatness error.Keywords: Leastaquare method ; Flatness error；Best met

4、hod目录 TOC o 1-3 h z u HYPERLINK l _Toc440573341第一章平面度的测量方法 PAGEREF _Toc440573341 h 3HYPERLINK l _Toc4405733421.1引言 PAGEREF _Toc440573342 h 3HYPERLINK l _Toc4405733431.2平面度误差的测量 PAGEREF _Toc440573343 h 3HYPERLINK l _Toc440573349第二章平面度的评定 PAGEREF _Toc440573349 h 3HYPERLINK l _Toc44057335021最小区域法 PAGER

5、EF _Toc440573350 h 3HYPERLINK l _Toc4405733562.2最小二乘法测量平面度误差的原理 PAGEREF _Toc440573356 h 3HYPERLINK l _Toc4405733572.21建立被测实际外表的数学模型 PAGEREF _Toc440573357 h 3HYPERLINK l _Toc440573358第三章用MATLA实现的过程 PAGEREF _Toc440573358 h 3HYPERLINK l _Toc4405733593.1软件编程 PAGEREF _Toc440573359 h 3HYPERLINK l _Toc4405

6、733603.2平面度误差的最小二乘法评定及其评定结果的不确定度 PAGEREF _Toc440573360 h 3HYPERLINK l _Toc4405733613.3小结 PAGEREF _Toc440573361 h 3-. z.第一章平面度的测量方法1.1 引言平面度误差是指被测实际外表对其理想平面的变动量。理想平面是评定平面度误差的评定基准，而评定基准的方位不同求得的平面度误差值也不同。假设用水平仪、自准直仪按节距法测量实际外表上各点相对于测量基准的平面度误差时，确定评定基准的方法可用：简便法、最小二乘法、和最小包容区域法2。本文着重分析介绍最小二乘法来确定评定基准，从而求得平面度

7、误差值。最小二乘法能准确而充分地利用全部的原始观测数据提供的信息，比拟客观地评定出不需要经过屡次试算的平面度误差，而且可直接运用于电子计算机运算，列入MATLAB软件，是的平面度误差的计算到达迅速、准确、可靠7。1.2 平面度误差的测量平面度误差的测量方法很多，常用的有如下所列的方法： = 1 * GB3 光隙法：将被测直线和测量基线间形成的光隙，与标准光隙比拟，测量不同方向的假设干个截面上的直线度误差，取其中最大值，作为平面度误差近似值的方法。该方法适用于磨削或研磨加工的小平面的平面度误差测量。 = 2 * GB3 指示器法：将被测零件支撑在平板上，平板工作面为测量基准，按一定的方式布点，如

8、图3所示，用指示器对被测面上各点进展测量并记录所测数据，然后，按一定的方法评定其误差值。该方法适用于中、小平面的平面度误差测量。 = 3 * GB3 光轴法：以几何光轴建立测量基面，测出被测面相对测量基面的偏离量，进而评定平面度误差值的方法，该方法适用于一般精度大平面的平面度误差测量。 = 4 * GB3 干预法：利用光波干预原理，根据干预条纹形状、条数，来确定平面度误差值的方法，该方法适用于精研外表的平面度误差测量。 = 5 * GB3 三坐标测量机：三坐标测量机是综合利用精细机械、微电子、光栅和激光干-. z.涉仪等先进技术的测量仪器，其原理是在三个相互垂直的方向上，有导向机构、测长元件、

9、数显装置，有一个能够放置工件的工作台，测头可以以手动或机动方式，轻快地移动到被测点上，由读数设备和数显装置，把被测点的坐标值显示出来的一种测量设备4。除上面介绍的几种方法外，还有液面法、自准直仪法等，工贸职业中专学校实验室在教学中，常采用光学合像水平仪测定平板的平面度误差。第二章平面度的评定平面度误差是被测实际外表相对理想平面的变动量。根据GBT1 1337-2004平面度误差检测对形状误差的定义，理想平面的位置应符合最小条件，平面度误差大小，等于包容实际外表且距离为两平行平面之间的宽度。平面度误差的评定方法有：最小包容区域法，对角线法，三远点法和最小二乘法5。2.1 最小区域法它是两平行包容

10、平面与实际被测要素的接触状态符合平面度误差判别法中*一准则时，此时两平行平面之问的距离，为平面度误差。用最小条件评定平面度误差有3种判断准则：三角形准则，穿插准则，直线准则。平面度误差的最小区域判别方法是：由两个平行平面包容被测实际外表时，至少有3点或4点相接触，接触点的上下分布，有以下3种形式之一者，即符合最小区域。 = 1 * GB3 三角形准则：如果被测实际外表上有3个最高(低)点及一个最低(高)点分别与两个包容平面相接触，并且最高(低)点能投影到3个最低(高)点之间，形成最小包容，称为三角形准则。 = 2 * GB3 穿插准则：两个高点与两个低点*些实际外表具有鞍形的形状特征，其与上下

11、包容面各有两个接触点，假设两高点的投影位于两低点连线之两侧，形成最小包容，称为穿插准则。 = 3 * GB3 直线准则：如果被测外表上的同一截面有两个最高(低)点与一个低(高)点分别和两个平行的包容面相接触，并且有一个最低(高)点的投影要落在两高(或两低)点连线之上，此时也形成最小包容，称为直线准则。 = 4 * GB3 对角线法：评定基准平面是通过实际平面的一条对角线和平行于另一条对角线的平面，各测点对此平面的偏差中最大正值与最大负值的绝对值之和为被测实际平面的平面度误差值。 = 5 * GB3 三点法：评定基准平面是通过被测实际平面上相距较远且不在一条直线上的三点建立的平面。偏离此平面的最

12、大值和最小值的绝对值之和为平面度误差这三种方法虽然都是针对同一被测平面，但由于评定基准平面选取不同，得到的误差值也不同6。2.2 最小二乘法测量平面度误差的原理最小二乘平面是个理想的平面，它使从实际被测轮廓上各点到该平面的距离的平方和为最小。因此，最小二乘法的目标函数，即各测点到最小二乘基准平面的距离的平方和FA,B,C=2-1其满足最小化时，FA,B,C的A,B,C即最小二乘基准平面的法向量， (2-2)即为最小二乘法平面度误差值3。2.2.1 建立被测实际外表的数学模型平面度误差是指被测实际外表不平的程度，而平面在空间直角坐标系中，它的一般方程为;A*+BY+CZ+D=0; 假设测得的点的

13、坐标是=1,2，,其中A,B,C为待定参数。根据最小二乘法原理得到最小二乘平面的待定系数A,B,C为A=2-3 B=2-4) C=(2-5)由于现在知道A,B,C 的系数，即确定了最小二乘平面。接下来就是就能知道改平面的上方的点和下方的点到该平面的距离，由于在该平面上下方都有点，所以点的本身就带有矢量性。点到平面的距离公式为2-6根据最小二乘法的原理有：2-7由于观测到的数据位于最小二乘平面基准平面的上下两侧，所以有点到平面的距离有最大值和最小值的两个点！即平面度为2-8第三章用MATLA实现的过程3.1 软件编程本文是通过MATLAB软件求得最小二乘平面基准平面的系数A,B,C，将最小二乘平

14、面基准平面拟合在三维坐标中和求出基准平面的平面度误差。根据观测到的数据如下：*mmymmzmm110.132 43.982 -0.003 219.641 44.092 -0.001 330.444 44.217 -0.002 430.618 29.359 -0.002 517.741 29.210 -0.001 65.450 29.068 -0.007 74.912 46.453 -0.012 815.887 46.580 -0.004 929.745 46.740 -0.002 1030.104 15.746 0.003 1118.721 15.614 -0.004 125.738 15.4

15、64 -0.010 1317.630 47.032 -0.010 1429.966 47.174 -0.007 155.989 44.945 -0.011 1616.808 45.071 -0.005 1728.761 45.209 -0.002 1828.901 33.125 0.000 1917.108 32.988 -0.007 203.975 32.836 -0.008 214.188 14.508 -0.008 2217.690 14.664 0.003 2330.909 14.817 0.004 2430.140 32.118 -0.003 2518.378 31.981 -0.0

16、07 264.371 31.820 -0.004 2718.421 31.982 -0.008 2830.889 32.127 -0.004 2931.111 13.176 0.001 3017.422 13.018 -0.005 现在将MATLAB代码输入软件中：function testlsqnonlin111clc;注解：将数据放入到e*cel表格中如以下图所示clear all;close all;% 读入数据data = *lsread(数据.*ls);% 3列对应*,y,z*data = data(:,2);ydata = data(:,3);注解：将保存好的数据文件与MATLAB

17、的代码文件放到同一个文件夹下。然后就可以启动MATLAB软件翻开代码文件运行一下，就可以出来结果图啦。祝你好运！zdata = data(:,4);% 画出原始数据点，红色figure;plot3(*data,ydata,zdata,r*);hold on;% 初始值 a b c 都为0*0 = 0 0 0;% 调用最小二乘法非线性拟合函数* =lsqnonlin(*) myfun(*,*data,ydata,zdata),*0)% 求平面度误差% 平面方程的系数a = *(1);b = *(2);c = *(3);% 根据公式求 didi = (zdata - a*data - b*ydat

18、a - c)/sqrt(a2+b2+1);dma* = ma*(di);dmin = min(di);f = dma* - dmin % 就是要求的平面度误差% 画拟合的平面* = linspace(min(*data),ma*(*data),200);yy = linspace(min(ydata),ma*(ydata),200);*,Y = meshgrid(*,yy);Z = a*+b*Y+c;surf(*,Y,Z);colormap;holdon;grid on;% 按照这个函数去拟合，使得式(2)成立function F = myfun(*,*data,ydata,zdata) F

19、= (zdata - *(1)*data - *(2)*ydata - *(3)/sqrt(1 + *(1)2+*(2)2);得到结果为：图1 最小二乘平面系数与平面度误差图2 MATLAB最小二乘法平面方程三维图3.2 平面度误差的最小二乘法评定及其评定结果的不确定度测量不确定度是指测量结果变化的不肯定，是表征被测量的真值在*个量值围的一个估计，是测量结果含有的一个参数，用以表示被测量值的分散性。这种测量不确定度的定义说明，一个完整的测量结果应该包含被测量值的估计与分散性参数两局部9。根据测量不确定度定义，在测量实践中如何对测量不确定度进展合理的评定，这是必须解决的根本问题。对于一个实际测量

20、过程，影响测量结果的精度有多方面的因素，因此测量不确定一般包含假设干个分量，各确定度分量不管性质如何，皆可用两类方法进展评定，即A类评定与B类评定。其中一些分量由一系列观测数据的统计分析来评定，称为A类评定；另一些分量不是用一系列观测数据的统计分析法，而是基于经历或其他信息所认定的概率分布来评定，称为B类评定。所有的不确定度分量均用标准差表征，它们或是由随机误差而引起的，或是由系统误差而引起的，都对测量结果的分散性产生相应的影响10。现在将平面方程的系数A,B,C的表达式分别展开令V=3-1P= QUOTE ，q=， QUOTE ， QUOTE 整理后得到 QUOTE 3-2 QUOTE QU

21、OTE QUOTE 设按最小二乘法评定公式得到最大偏离值点和最小偏离值的坐标分别为,带入公式得到 QUOTE 令 QUOTE , QUOTE , 再令 QUOTE QUOTE QUOTE 得到 QUOTE 展开 QUOTE 转化为矩阵形式，令 QUOTE QUOTE 得到 QUOTE 因被测点的*,Y坐标值在平面度的非敏感方向上，其不确定度可以忽略，因此上式为坐标值Z的线性函数，假设被测点的测量结果间互相独立，且具有一样的标准差，即有一样的不确定度则可以用不确定度的公式，得到用最小二乘法评定的平面度的不确定度为假设各点被测点直接测量结果间有关，且各自确实定度不同，则要计入协方差的影响，为将表达

22、式简化起见，将矩阵增广，即当时，令,当时，这时候最小二乘法评定的平面度误差的不确定度为 QUOTE 由于测量了30组数据，所以n=30。由e*cel表格就可求出，所以P=571.190、q=965.116 u=13762.216 w=35543.224 t=18360.108 v=3.86*。由于最小二乘平面的方程系数已经算出来，即A=0.0003，B=-0.0001，C=-0.0058,在利用点到e*cel计算出最小二乘平面的距离3-3所以能够得到偏离最小二乘平面的最大值和最小值所对应的坐标17.690,14.664,0.003此为最大值，29.966,47.174，-0.007为最小值，将数据带入到公式中得到的结果不确定度为 QUOTE 。3.3 小结：最小二乘法评定平面度误差从数理统计分析观点来看，能够准确地和充分地处理原始观测数据所提供的信息，比拟客观地评定平面度误差，有利于采用电子计算机进展运算，尤其适用于测量较大平面的平面度误差，测量较大平面时，因其布点较多，假