最小二乘拟合法在工程测量上的应用论文_杨丽

摘 要 利用全站仪对地铁盾构横断面进行观测，根据测量数据通过拟合分析其变形情况，本文利用最小二乘法对采集的数据进行处理，根据所得的结果分析其断面的变形情况。

关键词 平面拟合；圆拟合；平整度；圆度变形；最小二乘

1.引言

对于空间圆形物体，在实际和设计时会存在差异或者是在长期使用后会发生形变而影响

其使用情况，因此必须对空间圆形物体进行检测，判断其是否满足工程要求[1]。

在地铁、隧道等的施工过程中，管片以及钢模拼接时经常要进行设计断面与实际断面的比较检测，由于其断面都是一个空间圆形，所以都可以看作为是空间圆形物体的检测问题。本文以盾构横断面为例，阐述检测空间圆形物体变形的方法及步骤，传统的空间圆的拟合方法，把三维空间的圆形看成是二维的平面。本文利用最小二乘拟合法对地铁盾构横断面的实测数据进行处理并分析断面的变形情况。

2.观测及算法原理

2.1 数据采集

在断面四周适当位置安置棱镜，棱镜的位置可根据工程的要求分布，在面前任意找一点安置全站仪，确定测站坐标和零方向，这样就形成了一个自定义的坐标系 o-xyh即可进行测量，测量各点的平面坐标和高程，若要提高精度，可以增加测回数[3]。

2.2 最小二乘拟合法

2.2.1 平面拟合

平面方程为

按最小二乘求出平面方程的四个参数a、b、c、d来确定此空间平面，其中权阵p可以取为单位阵或者根据工程要求确定。

2.2.2 坐标转换

为了方便计算，将测量坐标系中的点位坐标转换为平面坐标系中坐标。在求得所有测定

利用上述转换方法，可以将投影点的坐标由测量坐标系转换到平面坐标系中，用于平面圆的拟合。

2.2.3 圆拟合

经过以上转换，所有测定点都由o-xyh坐标系中的点归算到了平面上，以表示各投影点在oo-xxyyhh中的平面坐标。

设圆的方程为：

按便可求得圆心坐标和圆半径，权阵可选为单位阵或者也可根据实际情况确定，式中的就是各点至圆周的距离也就是空间圆度的变形量。

3.算例分析

表1是在盾构断面测定的8个坐标点，为了拟合的准确性，需要在圆周上均匀的布点。

　表1 原始数据

表2 两种拟合方法结果比较

根据平面拟合算法编程实现空间平面拟合，得出此空间圆所在的平面方程为：

所得圆的参数为圆心（10.986949308，7.554189997，1.995575827）、半径5.012236628m。此方法求出来的半径大于断面的设计半径5m，说明断面的半径增大了12.2mm。与设计圆心（10.9893，7.535，2.000）相比，在X、Y、Z三个方向的坐标差分别为-2.351mm、+19.19mm、-4.423mm,所以实际圆心位置偏离设计圆心位置的距离为19.83296mm。

表2中的点面距反映了该空间平面的平整度，点圆距反映了圆度变形，所以表中的点面距和点圆距就反应了盾构断面的变形情况。

4.结论

本文介绍了空间圆形物体的检测方法，利用全站仪观测盾构断面的空间坐标，拟合出空间平面方程和圆的方程，各观测坐标至平面的距离即可反映物体平面的平整度，各点至圆周的距离就是该圆的变形程度，实例中断面半径比设计时增大了12.2mm且圆心偏离设计位置19.83296mm，各点的变形情况如图2所示。

参考文献：

[1]王解先,季凯敏.工业测量拟合[M].北京:测绘出版社,2008.

[2]潘国荣，陈晓龙.空间圆形物体拟合新方法[J].大地测量与地球动力学.2008年4月第28卷第2期：92-94

[3]王解先，赵向阳.圆形轨道变形测量[J].工程勘察.2003年第4期：60-62.

[4]李清玉.三维坐标测量数据处理及可视化技术研究[D].哈尔滨工业大学，2004.

[5]陈基伟.工业测量数据拟合研究[D].同济大学，2005.

[6]王清泉，闫保旭，王解先.长方形立面检测方法研究[J].测绘通报.2008年第8期：48-50.

[7]田洪祯，田林.空间长方形设备监测[J].科技信息.2007年第11期：46-47.

作者简介：杨丽(1985~)女，江苏淮安人，上海市城市建设工程学校（上海市园林学校），测量专业课教师

论文作者:杨丽

论文发表刊物:《城镇建设》2019年17期

论文发表时间:2020/4/13

标签：断面论文;  平面论文;  坐标论文;  空间论文;  测量论文;  盾构论文;  圆形论文;  《城镇建设》2019年17期论文;