1、平面度误差的测量及数据处理学校代码： 保密学 号： 保密课程设计说明书题 目: 平面度误差的测量及数据处理学生姓名：学 院：班 级：指导教师： 摘要平面度误差是将被测实际表面与理想平面进行比较，两者之间的线值距离即为平面度误差值；或通过测量实际表面上若干点的相对高度差，再换算以线值表示的平面度误差值。本文就平面度误差的数学模型与按最小二乘法建立理想平面（评定基准）的数学模型展开分析讨论；并结合案列分析，得出比较客观地评定平面度误差或者测量较大平面的平面度误差，最小二乘法是最佳方法1关键词：最小二乘法；平面度误差；最佳方法AbstractFlatness error is measured by
2、 actual surface with an ideal plane are compared, 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 t
3、he mathenatical 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；B
4、est method第一章 平面度的测量方法1.1 引言平面度误差是指被测实际表面对其理想平面的变动量。理想平面是评定平面度误差的评定基准，而评定基准的方位不同求得的平面度误差值也不同。若用水平仪、自准直仪按节距法测量实际表面上各点相对于测量基准的平面度误差时，确定评定基准的方法可用：简便法、最小二乘法、和最小包容区域法2。本文着重分析介绍最小二乘法来确定评定基准，从而求得平面度误差值。最小二乘法能准确而充分地利用全部的原始观测数据提供的信息，比较客观地评定出不需要经过多次试算的平面度误差，而且可直接运用于电子计算机运算，列入MATLAB软件，是的平面度误差的计算达到迅速、准确、可靠7。1.2
5、 平面度误差的测量平面度误差的测量方法很多，常用的有如下所列的方法：光隙法：将被测直线和测量基线间形成的光隙，与标准光隙比较，测量不同方向的若干个截面上的直线度误差，取其中最大值，作为平面度误差近似值的方法。该方法适用于磨削或研磨加工的小平面的平面度误差测量。指示器法：将被测零件支撑在平板上，平板工作面为测量基准，按一定的方式布点，如图3所示，用指示器对被测面上各点进行测量并记录所测数据，然后，按一定的方法评定其误差值。该方法适用于中、小平面的平面度误差测量。光轴法：以几何光轴建立测量基面，测出被测面相对测量基面的偏离量，进而评定平面度误差值的方法，该方法适用于一般精度大平面的平面度误差测量。
6、干涉法：利用光波干涉原理，根据干涉条纹形状、条数，来确定平面度误差值的方法，该方法适用于精研表面的平面度误差测量。三坐标测量机：三坐标测量机是综合利用精密机械、微电子、光栅和激光干涉仪等先进技术的测量仪器，其原理是在三个相互垂直的方向上，有导向机构、测长元件、数显装置，有一个能够放置工件的工作台，测头可以以手动或机动方式，轻快地移动到被测点上，由读数设备和数显装置，把被测点的坐标值显示出来的一种测量设备4。除上面介绍的几种方法外，还有液面法、自准直仪法等，娄底工贸职业中专学校实验室在教学中，常采用光学合像水平仪测定平板的平面度误差。第二章 平面度的评定平面度误差是被测实际表面相对理想平面的变动
7、量。根据GBT1 1337-2004平面度误差检测对形状误差的定义，理想平面的位置应符合最小条件，平面度误差大小，等于包容实际表面且距离为两平行平面之间的宽度。平面度误差的评定方法有：最小包容区域法，对角线法，三远点法和最小二乘法5。2.1 最小区域法它是两平行包容平面与实际被测要素的接触状态符合平面度误差判别法中某一准则时，此时两平行平面之问的距离，为平面度误差。用最小条件评定平面度误差有3种判断准则：三角形准则，交叉准则，直线准则。平面度误差的最小区域判别方法是：由两个平行平面包容被测实际表面时，至少有3点或4点相接触，接触点的高低分布，有下列3种形式之一者，即符合最小区域。三角形准则：如
8、果被测实际表面上有3个最高(低)点及一个最低(高)点分别与两个包容平面相接触，并且最高(低)点能投影到3个最低(高)点之间，形成最小包容，称为三角形准则。交叉准则：两个高点与两个低点某些实际表面具有鞍形的形状特征，其与上下包容面各有两个接触点，若两高点的投影位于两低点连线之两侧，形成最小包容，称为交叉准则。直线准则：如果被测表面上的同一截面内有两个最高(低)点与一个低(高)点分别和两个平行的包容面相接触，并且有一个最低(高)点的投影要落在两高(或两低)点连线之上，此时也形成最小包容，称为直线准则。对角线法：评定基准平面是通过实际平面的一条对角线和平行于另一条对角线的平面，各测点对此平面的偏差中
9、最大正值与最大负值的绝对值之和为被测实际平面的平面度误差值。三点法：评定基准平面是通过被测实际平面上相距较远且不在一条直线上的三点建立的平面。偏离此平面的最大值和最小值的绝对值之和为平面度误差这三种方法虽然都是针对同一被测平面，但由于评定基准平面选取不同，得到的误差值也不同6。2.2 最小二乘法测量平面度误差的原理最小二乘平面是个理想的平面，它使从实际被测轮廓上各点到该平面的距离的平方和为最小。因此，最小二乘法的目标函数，即各测点到最小二乘基准平面的距离的平方和F（A,B,C）= （2-1）其满足最小化时，F（A,B,C）的（A,B,C）即最小二乘基准平面的法向量， (2-2)即为最小二乘法平
10、面度误差值3。 2.2.1 建立被测实际表面的数学模型平面度误差是指被测实际表面不平的程度，而平面在空间直角坐标系中，它的一般方程为;AX+BY+CZ+D=0; 假设测得的点的坐标是（）（=1,2，,）,其中A,B,C为待定参数。根据最小二乘法原理得到最小二乘平面的待定系数A,B,C为A= （2-3） B= （2-4) C= (2-5)由于现在知道A,B,C 的系数，即确定了最小二乘平面。接下来就是就能知道改平面的上方的点和下方的点到该平面的距离，由于在该平面上下方都有点，所以点的本身就带有矢量性。点到平面的距离公式为 （2-6）根据最小二乘法的原理有： （2-7）由于观测到的数据位于最小二乘
11、平面（基准平面）的上下两侧，所以有点到平面的距离有最大值和最小值的两个点！即平面度为 （2-8）第三章 用MATLA实现的过程3.1 软件编程本文是通过MATLAB软件求得最小二乘平面（基准平面）的系数A,B,C，将最小二乘平面（基准平面）拟合在三维坐标中和求出基准平面的平面度误差。根据观测到的数据如下：x（mm）y（mm）z（mm）110.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 7
12、4.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.464 -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
13、.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.007 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;注解：将数据放入到excel表格中
14、如下图所示clear all;close all;% 读入数据data = xlsread(数据.xls);% 3列对应x,y,zxdata = data(:,2);ydata = data(:,3);注解：将保存好的数据文件与MATLAB的代码文件放到同一个文件夹下。然后就可以启动MATLAB软件打开代码文件运行一下，就可以出来结果图啦。祝你好运！zdata = data(:,4);% 画出原始数据点，红色figure;plot3(xdata,ydata,zdata,r*);hold on;% 初始值 a b c 都为0x0 = 0 0 0;% 调用最小二乘法非线性拟合函数x =lsqnon
15、lin(x) myfun(x,xdata,ydata,zdata),x0)% 求平面度误差% 平面方程的系数a = x(1);b = x(2);c = x(3);% 根据公式求 didi = (zdata - a*xdata - b*ydata - c)/sqrt(a2+b2+1);dmax = max(di);dmin = min(di);f = dmax - dmin % 就是要求的平面度误差% 画拟合的平面xx = linspace(min(xdata),max(xdata),200);yy = linspace(min(ydata),max(ydata),200);X,Y = mesh
16、grid(xx,yy);Z = a*X+b*Y+c;surf(X,Y,Z);colormap;holdon;grid on;% 按照这个函数去拟合，使得 式(2)成立function F = myfun(x,xdata,ydata,zdata) F = (zdata - x(1)*xdata - x(2)*ydata - x(3)/sqrt(1 + x(1)2+x(2)2);得到结果为：图1 最小二乘平面系数与平面度误差图2 MATLAB最小二乘法平面方程三维图3.2 平面度误差的最小二乘法评定及其评定结果的不确定度测量不确定度是指测量结果变化的不肯定，是表征被测量的真值在某个量值范围的一个估
17、计，是测量结果含有的一个参数，用以表示被测量值的分散性。这种测量不确定度的定义表明，一个完整的测量结果应该包含被测量值的估计与分散性参数两部分9。根据测量不确定度定义，在测量实践中如何对测量不确定度进行合理的评定，这是必须解决的基本问题。对于一个实际测量过程，影响测量结果的精度有多方面的因素，因此测量不确定一般包含若干个分量，各确定度分量不论性质如何，皆可用两类方法进行评定，即A类评定与B类评定。其中一些分量由一系列观测数据的统计分析来评定，称为A类评定；另一些分量不是用一系列观测数据的统计分析法，而是基于经验或其他信息所认定的概率分布来评定，称为B类评定。所有的不确定度分量均用标准差表征，它
18、们或是由随机误差而引起的，或是由系统误差而引起的，都对测量结果的分散性产生相应的影响10。现在将平面方程的系数A,B,C的表达式分别展开 令V= （3-1） P=，q=，整理后得到 （3-2） 设按最小二乘法评定公式得到最大偏离值点和最小偏离值的坐标分别为,带入公式得到 令 , , 再令 得到展开 转化为矩阵形式，令 得到因被测点的X,Y坐标值在平面度的非敏感方向上，其不确定度可以忽略，因此上式为坐标值Z的线性函数，若被测点的测量结果间互相独立，且具有相同的标准差，即有相同的不确定度则可以用不确定度的公式，得到用最小二乘法评定的平面度的不确定度为若各点被测点直接测量结果间有关，且各自的确定度不
19、同，则要计入协方差的影响，为将表达式简化起见，将矩阵增广，即当时，令,当时，这时候最小二乘法评定的平面度误差的不确定度为 由于测量了30组数据，所以n=30。由excel表格就可求出，所以P=571.190、q=965.116 u=13762.216 w=35543.224 t=18360.108 v=3.86x。由于最小二乘平面的方程系数已经算出来，即A=0.0003，B=-0.0001，C=-0.0058,在利用点到excel计算出最小二乘平面的距离 （3-3）所以能够得到偏离最小二乘平面的最大值和最小值所对应的坐标（17.690,14.664,0.003）此为最大值，（29.966,47
20、.174，-0.007）为最小值，将数据带入到公式中得到的结果不确定度为。3.3 小结：最小二乘法评定平面度误差从数理统计分析观点来看，能够准确地和充分地处理原始观测数据所提供的信息，比较客观地评定平面度误差，有利于采用电子计算机进行运算，尤其适用于测量较大平面的平面度误差，测量较大平面时，因其布点较多，若不用最小二乘法，则往往需要经过许多次试算，使得计算变得繁琐。而采用最小二乘法只需简单的四则运算，不仅较易掌握，且快迅速准确地求得其平面度误差。参考文献1陈永鹏,基于MATLAB的平面度评定方法D.中国工程物理研究院.20042李高峰,陈晓怀，王宏涛,李红莉.基于MATLAB的平面度评定D.合
21、肥工业大学. 20123刘志平，石林英.最小二乘法及其MATLAB实现D.淮海工学院 理学院.20084张少南，曹树文.矩阵最小二乘法评定平面度误差J.大庆石油学院学报.1989年6月第十三卷第2期5张善镬，车仁，生胡远.利用计算机处理平面度误差的一种新方法迭代逼近法J.仪器仪表学报.1986年8月第7卷第3期6余勃强,崔亚辉.刘西安中大平面的平面度检测D.西安理工大学.20077王伟.一种利用Matlab软件计算平面度的方法D.广东省计量科学研究院东莞分院.20108邹春龙,黄浩，王生怀,邓小雯.平面度误差三维测量系统的研究J.仪表技术与传感器. 2014年第9期9陈飞,闵淑辉,徐学林.旋转
22、法处理平面度误差测量数据D.中南林学院.2004310陆晓珩.平面度误差测量不确定度评定D.甘肃省计量研究院.2008.1011Radlovaki,VladanHadistevi, Miodragtrbac, BrankoDeli, MilanKamberovi, new methodBundle of plains through one pointJ. University of Novi Sad, Faculty of Technical Sciences, Department of Industrial Engineering and Management, Trg Dositeja Obradovia 6.Jan2016, Vol. 43.
