一种异形钢塔节段姿态快速定位测量方法与流程

本发明涉及土建工程技术领域，特别涉及一种异形钢塔节段姿态快速定位测量方法。

背景技术：

新世纪以来，伴随我国城市化快速发展的进程，城市桥梁建设规模与日俱增，桥梁已被纳入城市建筑的构成元素，除满足基本的交通功能外，被赋予更多的象征功能和文化内涵，桥梁美学需求更加突出，具有复杂空间造型美学感染力的异形钢塔斜拉桥受到广泛关注，桥梁自身的安全性、美学对几何外观的敏感性均对复杂造型的钢塔建造精度提出了严格的要求，如直立焊接式钢塔架设偏位控制标准为h/4000(h:塔高)，栓接式钢塔架设偏位控制标准为h/3000。

逐节段架设是钢塔建造的常用施工方法，在现场节段架设过程中，首先对节段下端口与前续已安装节段上端口进行匹配使得各壁板错边量满足焊接条件，此时需对该节段上端口姿态进行复核，确保该节段架设几何形态处于误差允许范围之内，其中轴线偏离度是衡量节段架设质量的一项重要指标，而上端口轴心点是一个人为设定的空间虚拟点，不具备直接测量的条件，一种可行的方法是对表征上端口架设质量的各特征点坐标进行复核，利用加权评定的方式或6d最优空间姿态推算上端口轴心的偏离度，按照仪器架设位置的不同其又可分为内控法与外控法两种方式。

内控法是将仪器直接架设在当前节段之上，通过照准地面已知坐标点的控制点利用后方交会法获取仪器架设点的坐标然后对节段上各特征点坐标进行采集，因所有待测点均位于由同一套控制点组成的测量控制网之内、且测量仪器距待测点距离较近，该方法具有各待测点可由一台仪器同时测量且相对位置测量精度高的优点，但该方法需要在节段之上布设测站，当上端口倾斜时则难以直接实施，另外较高时易受如大风等环境的影响；

外控法是将仪器架设在地面控制点之上，在节段上端口待测点之上布设靶标，对其坐标进行采集，由于受结构壁板或现场其他施工因素的遮挡干扰，需要结合场地条件变换测站方可完成所有特征点的坐标采集，因受待测点位于不同测量控制网、一台仪器无法完成同时测量等因素的制约，测量精度相对内控法较低。大型空间异形结构的定位是工程测量的难题，如董伟东针对国家体育场“鸟巢”倾斜柱研究了解析法与旋转矩阵法确定特征点坐标的测设方法，卓新等针对浙江美术馆空间多杆件相贯节点提出了异形构件的虚拟四面体顶点测量定位方法。

针对桥梁空间造型复杂，采用内控法对特征点坐标进行采集实施难度大，又因其所呈现的空间扭曲造型导致需要变换不同测站方能完成上端口坐标的采集进而影响其测量精度及效率。

因此，如何针对桥梁空间造型复杂的节段姿态进行准确快速定位，目前还没有很好的方法。

技术实现要素：

本发明的目的在于可针对桥梁空间造型复杂的节段姿态进行准确快速定位，提出的一种异形钢塔节段姿态快速定位测量方法，可提高钢塔架设测控效率，可解决现有技术条件下难以准确定位测量的问题。

为了解决上述技术问题，本发明实施例提供一种异形钢塔节段姿态快速定位测量方法，包括以下步骤：

s1、在工厂内根据结构几何形态进行节段几何特征点及测控点布局；构建临时坐标系并采集几何形态特征点的坐标数据；所述特征点包括：所述节段上端口控制点、下端口控制点、节段上端口轴心点、节段下端口轴心点；所述测控点为临近节段上端口预设距离处布设的三个测点；所述特征点及测控点均布设在节段外壁板相交处；

s2、构建第一局部坐标系；根据三点获取临时坐标系与第一局部坐标系传递关系；

s3、节段架设完毕后构建大地坐标系以及第二局部坐标系；根据三点获取大地坐标系与第二局部坐标系传递关系；

s4、根据三点获取临时坐标系与大地坐标系传递关系；

s5、根据上端口坐标在临时坐标系及大地坐标系之间的传递关系，获取节段几何形态特征点在大地坐标系下的坐标；将上端口轴心坐标与设计坐标对比判定节段定位状态。

在一个实施例中，所述步骤s1基于测量误差椭圆曲线及随机分布特性，依据内点罚函数法寻找最优三点布局方案;

所述步骤s1包括：在工厂内构建临时坐标系，根据布局结果对所述节段几何特征点坐标进行采集。

在一个实施例中，所述步骤s2，基于临近节段上端口处布设的三个测点，构建第一局部坐标系；由三点获得临时坐标系与第一局部坐标系传递关系。

在一个实施例中，所述步骤s3包括：节段架设完毕后，构建大地坐标系，确定所述三点在大地坐标系下的坐标值；并根据三点构建第二局部坐标系；由三点获得大地坐标系与第二局部坐标系传递关系。

在一个实施例中，所述步骤s4，包括：根据三点获取临时坐标系与大地坐标系传递关系，对于上端口任一点其，据此获取临时坐标系与大地坐标系之间传递关系。

在一个实施例中，所述步骤s5包括：基于加权评定或6d最优空间姿态推算的方式建立三点与上端口轴心点的关系，将上端口轴心坐标与设计坐标对比判定节段定位状态。

在一个实施例中，所述的第一局部坐标系和第二局部坐标系，原点位于点，其x轴方向由点指向点，其z轴方向为三个测点组成的平面的法方向，y轴按右手法则坐标系设置。

在一个实施例中，所述的由三点获得临时坐标系与第一局部坐标系传递关系，包括：

根据公式(1)将所述节段上任一实测特征点坐标由临时坐标系变换至第一局部坐标系；

(1)

公式(1)中及分别为由临时坐标系为变换至第一局部坐标系的旋转矩阵及平移向量，为第一局部坐标系中的坐标值，为临时坐标系中的坐标值；计算方法见公式(2)；

(2)

公式(2)中为坐标系的x轴、y轴、z轴在坐标系中的单位坐标向量，计算方法见公式(3)；

(3)

公式(3)中及均由点在坐标系中的坐标计算得来。

在一个实施例中，所述的由三点获得大地坐标系与第二局部坐标系传递关系，包括：

根据公式(4)将所述节段上任一实测特征点坐标由大地坐标系为变换至第二局部坐标系；

(4)

公式(4)中及分别为由大地坐标系为变换至第二局部坐标系的旋转矩阵及平移向量，为第二局部坐标系中的坐标值；为大地坐标系中的坐标值；

其计算方法分别见公式(5)、(6)，其中及均由点在坐标系中的坐标计算得来；

(5)

(6)。

在一个实施例中，对于上端口任一点其，据此获取临时坐标系与大地坐标系之间传递关系，包括：

根据三点获取临时坐标系与大地坐标系传递关系，获取节段几何形态特征点在大地坐标系下的坐标，见公式(7)；

(7)。

所述为大地坐标系中的坐标值；为临时坐标系中的坐标值；

及分别为由大地坐标系为变换至第二局部坐标系的旋转矩阵的逆矩阵及平移向量；

及分别为由临时坐标系为变换至第一局部坐标系的旋转矩阵及平移向量。

本发明的优点在于，提出了一种异形钢塔节段姿态快速定位测量方法，针对空间造型复杂的桥塔节段，采用三点法进行桥塔节段姿态定位快速测量，可实现节段厂内与架设现场坐标传递与姿态衔接，便于现场架设过程中的快速定位与姿态测量，提高了钢塔架设测控效率。

本发明的其它特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点可通过在所写的说明书、权利要求书、以及附图中所特别指出的结构来实现和获得。

下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。

附图说明

附图用来提供对本发明的进一步理解，并且构成说明书的一部分，与本发明的实施例一起用于解释本发明，并不构成对本发明的限制。在附图中：

图1为本发明实施例提供的异形钢塔节段姿态快速定位测量方法流程图；

图2为永定河大桥的高塔横桥向及纵桥向视图；

图3为本发明实施例提供的节段几何形态控制特征点布局示意图；

图4为本发明实施例提供的厂内节段几何形态坐标系布局示意图；

图5为本发明实施例提供的大地坐标系下测量状态示意图；

图6为本发明实施例提供的单点误差椭圆示意图；

图7为本发明实施例提供的多点误差分布示意图；

图8为本发明实施例提供的多布局误差循环求解流程图；

图9a为本发明实施例提供的控制网形及测量平差示意图；

图9b为本发明实施例提供的gtn15上端口形态示意图；

图10a-f为本发明实施例提供的布局示意图；

图11a为本发明实施例提供的误差比值与缩放因子α曲线图；

图11b为本发明实施例提供的误差比值与夹角比值曲线图；

图12a为本发明实施例提供的误差比值与面积比值曲线图；

图12b为本发明实施例提供的近优状态下误差比值δ变化规律曲线图；

图13a为本发明实施例提供的又一误差比值与面积关系曲线图；

图13b为本发明实施例提供点位偏离横隔板对误差的影响曲线图；

图14a为本发明实施例提供的高塔合龙后轴心线顺桥向偏差分布示意图；

图14b为本发明实施例提供的高塔合龙后轴心线横桥向偏差分布示意图。

具体实施方式

下面将参照附图更详细地描述本公开的示例性实施例。虽然附图中显示了本公开的示例性实施例，然而应当理解，可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了能够更透彻地理解本公开，并且能够将本公开的范围完整的传达给本领域的技术人员。

参照图1所示，本发明实施例提供一种异形钢塔节段姿态快速定位测量方法，包括以下步骤s1~s5；

s1、在工厂内根据结构几何形态进行节段几何特征点及测控点布局；构建临时坐标系并采集几何形态特征点的坐标数据；所述特征点包括：所述节段上端口控制点、下端口控制点、节段上端口轴心点、节段下端口轴心点；所述测控点为临近节段上端口预设距离处布设的三个测点；所述特征点及测控点均布设在节段外壁板相交处；

s2、构建第一局部坐标系；根据三点获取临时坐标系与第一局部坐标系传递关系；

s3、节段架设完毕后构建大地坐标系以及第二局部坐标系；根据三点获取大地坐标系与第二局部坐标系传递关系；

s4、根据三点获取临时坐标系与大地坐标系传递关系；

s5、根据上端口坐标在临时坐标系及大地坐标系之间的传递关系，获取节段几何形态特征点在大地坐标系下的坐标；将上端口轴心坐标与设计坐标对比判定节段定位状态。

本实施例中，针对桥梁空间造型复杂的节段姿态，提出了基于构建坐标系提出了三点法钢塔节段姿态定位技术，可实现节段厂内与架设现场坐标传递与姿态衔接，便于现场架设过程中的快速定位与姿态测量，提高了钢塔架设测控效率。

下面本发明实施例以永定河大桥为例，对本发明的技术方案进行详细说明。

永定河大桥工程概况：

长安街西延永定河特大桥位于“神州第一街”长安街与北京“母亲河”永定河交汇处，主桥采用五跨高低双塔斜拉刚构连续梁组合桥体系，全长639m，主跨280m，为世界已建的最高拱形钢塔斜拉桥，塔柱为倾斜拱形造型，桥梁中线与河道中线斜交角为57.40°为适应河道水流方向，塔柱双肢非一致倾斜，形成“迈步”效果，呈空间扭曲造型。

该桥钢塔为双肢非一致倾斜、空间不对称扭转的变截面全焊接钢箱拱，其中高塔高123.780m，重约9850t，南北塔肢倾斜角分别为71.8°和62.0°，塔根顺桥向距离为25.100m，南北肢各划分15个节段，含合龙段整塔共计31个节段，钢塔外形及节段划分见图2所示，各节段均在桥址处钢结构总拼厂内组拼焊接，钢塔基于支架辅助安装法进行架设，架设过程中节段自重主要由自身承担，支架提供拼装过程的临时支撑力，同时为索塔线形的可能调整提供反力支撑。

钢塔设计轴线偏位控制标准为h/4000,塔顶处允许偏位约为31mm，高程允许偏差为±2n（n为节段号）且不超过20mm。该桥梁单肢钢塔顺桥向投影为一变截面线性倾斜的轮廓，横桥向投影为一变截面椭圆轮廓，节段轴线为由顺桥向倾斜直线与横桥向的椭圆曲线组成的复合空间曲线，各节段端口面为与此处空间曲轴线切线垂直的平面，在大地坐标系下端口为一空间倾斜平面。

将本发明实例提供的方法具体应用在永定河大桥节段定位上：

几何形态特征点的选择应反映索塔节段的三维特征，且基于特征点能够对节段上下端口的制造加工状态进行评定，另外轴线的偏离度是衡量索塔制造架设的一个综合指标，即几何形态控制点还需包括节段上下端口的轴心线点，单个节段几何形态控制点布局见图3。单个节段几何形态控制点共计13个，为节段上端口控制点，为节段下端口控制点，控制点均布设在可以直接获取坐标的节段外壁板相交处，为节段上端口轴心点，为节段下端口轴心点，上下端口轴心点为虚拟空间位置，在节段几何形态采集时并不直接涉及，为节段三点定位设置的测控点，其布设位置相对灵活，一般布设在上端口第一道横隔板与四边壁板相交的外侧附近，主要因为该处因横隔板的刚性约束各测点之间的相对位置比较稳定，且测点距上端口一般为0.5m~1.5m左右，上端口相对横隔板的变形可以忽略不计。

基于三点法节段姿态定位技术，首先在厂内对节段制造几何姿态进行采集时，需要在节段上端口临近位置布设三个测点，设厂内坐标采集时构建的临时坐标系为，基于此三个测点构建一个局部坐标系，可以将节段上端口所有特征点封装至由此三点构建的局部坐标系之内，局部坐标系原点位于点，其x轴方向由点指向点，其轴方向为三个测点组成的平面的法方向，y轴按右手法则坐标系设置，厂内节段坐标系布局见图4所示。

厂内任一实测特征点坐标由临时坐标系为变换至局部坐标系的计算方法见公式(1)。

(1)

公式(1)中及分别为由临时坐标系为变换至局部坐标系的旋转矩阵及平移向量，计算方法见公式(2)。

(2)

公式(2)中为坐标系的x轴、y轴、z轴在坐标系中的单位坐标向量，计算方法见公式(3)。

(3)

公式(3)中及均由点在坐标系中的坐标计算得来。

当节段架设后，仍可以按照上述思路进行转换，设大地坐标系为，基于此三个测点构建一个局部坐标系，与在厂内构建局部坐标系的法则相同，任一实测特征点坐标由临时坐标系为变换至局部坐标系的计算方法见公式(4)。

(4)

公式(4)中及分别为由大地坐标系为变换至局部坐标系的旋转矩阵及平移向量，其计算方法分别见公式(5)、(6)，其中及均由点在坐标系中的坐标计算得来。

(5)

(6)

由于至上端口较近，可以忽略结构局部变形所引起的上端口临近区域不同点之间相对位置的改变，因此对于上端口任一点其，进而可以推导出任一点坐标在及两个坐标系间的传递关系，见公式(7)。

(7)

图5描述了现场架设过程中坐标系的关系，当在节段制造几何形态采集过程中建立了上端口特征点与点的关系后，基于加权评定或6d最优空间姿态推算的方式可进一步建立点与上端口轴心点的关系，在现场架设过程中仅需对点进行坐标采集即可按照公式(7)获得大地坐标系下上端口的轴线坐标，将其坐标与目标坐标进行对比，即可判定上端口轴心点的定位状态，进而为节段姿态评定或调整提供依据。事实上通过合理的布设测点可以做到在同一测站完成三个点坐标的采集，避免了对上端口所有特征点采集时需要面临变换测站的问题，有效提高了测控效率，这些测点不但可以用在节段姿态定位环节中，还可以为节段的变形观测提供依托。

进一步地，为了验证上述过程中对节段定位的准确度，作为一种测量方法，仅给出其结果的计算公式是不够的，尚需对其误差进行研究，以确定其适用的具体场合。

误差分析计算过程：

将公式(1)带入公式(7)可得坐标在厂内临时坐标系为与大地坐标系之间的转换计算方法，见公式(8)。

(8)

由公式(8)可知，基于三点定位法所推算的坐标与临时坐标系为向局部坐标系坐标系转换的旋转矩阵、平移向量、坐标测量值及由大地坐标系为向局部坐标系坐标系转换的旋转矩阵、平移向量均存在关联，坐标的误差同样来源上述坐标直接测量结果或基于坐标的计算的旋转矩阵r及平移向量t，由于厂内为近距测量，不受测量环境、基准网变形等因素的影响，测量精度较高，可以认为由厂内临时坐标系为向局部坐标系坐标系转换误差可以忽略，进一步的可假设临时坐标系为与大地坐标系为同一坐标系，只是临时坐标系为所测结果为无误差状态，则可求得大地坐标系基于点坐标推测上端口某点的坐标偏差，具体见公式(9)。

(9)

点坐标基于全站仪测量，全站仪架设在控制网观测站上，受大气折光、仪器精度等因素的影响，单点坐标存在测量误差，误差分布是以待定点为极点、ψ为极角、为长度的极坐标点的轨迹，见图6，该轨迹可以描述各方向位差，误差曲线关于两个极轴(e轴和f轴)对称，由于误差曲线作图的不便，一般用形状与误差曲线相似椭圆代替，误差椭圆与误差曲线在两个极值方向完全重合，、e、f称为点位误差椭圆的参数，忽略任意ψ极角上误差曲线与误差椭圆的偏差（如od在ψ方向上误差椭圆的向径与误差曲线向径od相差很小，忽略的影响），任一点的点位方差总是等于两个相互垂直方向上的方差分量之和，与坐标系的选择无关，基于随机的分布特点可以认为测量点随机分布在误差椭圆上，即点随机分布在其误差椭圆上，如图7所示。

由公式(9)可知，基于点对某个如端口轴心点进行坐标推算的偏差与旋转矩阵、平移向量存在关联，由公式(5)、(6)可知，与点直接相关，因此不同点布局方案会导致不同的推算点的偏差。当点在节段端口上确定后，在一定的测量条件下对其坐标进行测量并进行平差获取各点坐标的误差椭圆曲线，基于当前误差椭圆曲线，可以寻找由此三点推算的某个如端口轴心点的最大误差，此误差可以评定基于当前点推算点的适用性，一种特点测点布局方案下寻找点最大误差的数学计算模型见公式(10)。

(10)

公式(10)中目标函数中为点椭圆误差轨迹上的点坐标，其平面上需满足位于椭圆上的约束条件，高程处于点高程误差范围内，即公式(10)目标函数需满足的约束条件见公式(11)。

(11)

公式(11)中第一部分为平面上中心点为点()、长轴2a、短轴2b、焦点间距2c、长轴与x轴夹角的椭圆方程。

由公式(10)、公式(11)构造的数学优化问题约束函数包含等式及不等式约束，目标函数为复杂的非线性函数，目标问题决策变量数目为3，难以采用数学解析方式获取目标函数的极值,对于含有约束条件的非线性目标函数最优化问题，惩罚函数法是一种广泛有效的间接求解方法，其核心思想是将约束非线性规划转化为一系列无约束问题，求解这一系列无约束问题,可得到约束非线性规划的解，具体可借助数学计算软件matlab的fmincon函数进行求解，通过变换点的布局可以分析点误差变化规律及与相关参数的敏感性，进而为确定点最优误差及优化点的布局提供指导，基于多点布局方案的点误差循环求解流程见图8，在matlab中编制循环程序可实现该计算流程。

本实例中，比如选取高塔北肢第15节段（gtn15）上端口轴心坐标测量进行误差分析，全站仪架设测站编号为d2，后视点分别为a2、b1，观测站距点最远距离约为302m，采用leicats60全站仪对坐标进行测量，测量结果平差由leica多测回数据分析系统完成，实测坐标及平差结果见表1，控制网形及误差椭圆分布见图9a，上端口、点实际布局见图9b，由测量结果可知误差椭圆参数均非常相近，这是由于三点在控制网内较为集中有关，三点之间最大距离仅为6.38m，三点最大长半轴为2.1mm，最大短半轴为1.6mm，误差椭圆长轴与大地坐标系夹角=147.5°。

表1gtn15测量及平差结果

基于上述所测误差椭圆参数，本节研究不同布局方案下所推测最大偏差及其影响因素，以便通过优化布局，控制最大偏差。

点设置在横隔板与四周壁板交接的外缘，结合平面控制的布设及通视条件设置点在角点间的行进路线，共考察了6种工况，参照图10a-f所示对应case1-case6，各工况设置要点如下：

①case1点设置在与中间，、自点开始分别沿相同的距离缩放因子α行进至、点；

②case2点设置在与中间，、分别自点开始沿相同的距离缩放因子α行进至、点；

③case3点设置在点，、分别自点开始沿相同的距离行进至，直至点行进至点，整个过程中三点所构三角形夹角维持不变；

④case4点设置在点，自点开始往点行进，点往点行进，整个过程中三点所构三角形面积维持不变；

⑤case5点设置在点，、分别自点开始沿相同的距离缩放因子α行进至、点，在α=1时，三点所构三角形面积最大；

⑥case6在中α=1、三点所构三角形面积最大时，沿竖向移动点的位置，分析点局部偏离横隔板对误差的影响。

case1及case2计算结果见表2，表中为当前步骤三点所构三角形面积，为当前三角构形最大最小角度比值，为所有步骤最大构形面积与当前构形面积比值，为当前步骤推算坐标最大偏差与三点最大点位偏差的比值，其中，误差比值δ随缩放因子α变化规律见图11a，可见伴随α的增大δ比值呈指数形式衰减，究其根本在于三点构形面积在不断增大、三点构形夹角不断趋于均匀的缘故，δ比值与夹角比值及面积比值的关系分别见图11b、图12a，可见随着三点构形夹角趋于均匀或面积逐渐增大，其误差比值δ呈线性减小，最优状态下case2状况下δ=1.32，已非常接近单点测量位差，位差绝对值为4.6mm，能够满足钢塔高空定位精度需求，说明基于三点法进行异形索塔定位是可行的。

表2case1及case2计算结果

需要说明的是虽然结果表明case2工况最优状态下位差绝对值为4.6mm，已能够满足钢塔高空定位精度需求，但由于端口异形及壁板遮挡的缘故需要变换测站方能完成三点坐标测量，这一方面会对测量工作带来较大不便性，同时测量周期的加长、测量控制网的不一致还会导致其他的误差，故该工况仅是一个计算工况，不具备实施条件，为此需要研究无需变换测站的布局方案，case3~case6工况均满足无需变换测站即可完成三点坐标测量的条件。

表3case3及case5计算结果

case3描述了构形夹角比值维持不变状态下误差比值δ随面积比值的变化规律，case4描述了构形面积维持不变状态下误差比值δ随夹角比值的变化规律，图12a给出了在近优状态下误差比值δ与夹角比值或面积比值的关系，近优状态下误差比值δ与夹角比值或面积比值基本呈线性关系，在1.0~4.0的夹角比值或面积比值取值范围内，相对于夹角比值，面积比值更为敏感，即变化速率更快，最优误差比值δ=1.85，可见若需进一步降低推算坐标误差比值δ的大小，在一定夹角比值的区间内寻找面积更大的构形。

该桥梁gtn14上端口中一种相对于case3更能达到面积最大的方案即为case5方案，图13a给出了在近优状态下误差比值δ与面积绝对值的关系，可以发现误差比值δ与面积呈非线性变化关系，虽然case5的构形夹角不如case3均匀，但由于其最优状态面积明显大于case3工况，其最优误差比值较case3所计算的结果更小，δ=1.50，位差绝对值为5.3mm，可以满足钢塔高空定位精度需求。

表4case4及case6计算结果

事实上，由于局部遮挡的影响，中个别点难以布设在横隔板与壁板的交接处，需要上下移动，case6分析了三点中某点沿竖向移动时误差比值δ的变化规律，具体见图13b，可以发现在局部上下调整点的位置其误差比值δ未发生明显改变，算例中上下各移动60mm误差比值δ变化仅为0.001，可见误差比值δ与某点的上下位置的局部调整关联性不强，现场点布设过程中在局部高度范围内可以灵活调整，此为该方法的应用带来了较大的便利性。

结合上述计算结果及误差比值δ影响因素分析，可以发现为使得采用三点法实施异形索塔定位点位推算误差最小，应设法使得构形面积最大、构形夹角均匀，按照该思路即可确定永定河大桥桥梁gtn15上端口最优布局方案，对应case5工况中α=1.0的状态，此时误差比值δ=1.50，位差绝对值为5.3mm，可以满足钢塔高空定位精度需求。

高塔合龙后顺桥向及横桥向分布见图14a-14b，图中n标识北肢，s标识南肢，0标识无偏差状态，u标识偏差上限，l标识偏差下限，e标识实际偏差，为提高分辨效果，图中对各偏差分项均扩大了200倍，实测结果表明钢塔轴心线偏差满足小于h/4000的要求。

本实施例中，北京永定河大桥桥梁钢塔为双肢非一致倾斜、空间不对称扭转的变截面全焊接钢箱拱，钢塔具有异形、超高、节段数目多、端口倾斜、架设精度要求高等特点，本文基于创建局部坐标系的思路提出了适用于异形钢塔节段姿态快速定位的三点法测量技术，并对其误差影响因素进行了分析，优化了点位布局方案，具有如下技术效果：

⑴能够实现现场架设过程中的节段快速定位与姿态测量，提高了钢塔架设测控效率

⑵基于测量误差椭圆曲线及随机分布特性，在matlab中编制程序循环调用带有约束条件的优化函数fmincon可实现对不同点位布局方案的最大误差分析。

⑶三点构形面积大小及三点构形夹角均匀程度是轴心坐标推算误差主要影响因素，优化点位布局后，推算点位差与三点最大位差比值δ=1.50，比如在永定河大桥桥梁测控网形下，推算点位差绝对值为5.3mm，可以满足钢塔高空定位精度需求。

⑷可实现钢塔成形后轴线偏差小于h/4000几何姿态控制要求。

最后所应说明的是，以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，对本发明的技术方案进行修改或者等同替换，都不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。