线路测量通用公式的推导及编程（很经典的论文）..docVIP

线路测量通用公式的推导及编程 【摘 要】用定积分的方法推导出求线路中线点坐标的积分通式，采用复化辛卜生（Simpson）公式给出其数值积分式。在此基础上推导出了求地面点与线路中线相对位置关系的公式。并进行工程实例计算，给出基于CASIO fx-4500P的计算器程序，大大提高了工作效率。 【关键词】测量；通用公式；编程；曲线元在线路中线及边线点坐标计算中，我们通常采用切线支距等公式计算各点坐标，但在不同线型上计算时就需要使用不同的公式，计算公式不具备通用性，给计算带来了不便。在缓和曲线计算时，螺旋线方程式的取项数对计算精度的影响随缓和曲线长度及圆曲线半径的不同而不同，取项不当会出现较大计算误差，取项太多又会增加计算难度。在道路的设计、施工、竣工阶段，经常需要测定地面点与线路中线的位置关系，以便为设计、施工及竣工提供准确数据。用常规方法及公式测定计算地面点与线路中线的位置关系难度相当大。针对以上问题，笔者用定积分的方法给出了求线路中线点的积分通式，并采用复化辛卜生公式给出其数值积分式，进而推导出了求地面点与线路中线相对位置关系的公式。并进行了工程实例计算及编程。 1 求线路中线坐标的通用公式推导 线路形式多种多样，但都是由直线、圆曲线及缓和曲线组成。而这3种线型都具有曲率随弧长作线性变化这一共性，故可用曲线元代表3种线型中的任意一种。由此易推出已知弧长l求任意点i的坐标的积分通式如下： （1） 式中：XA、YA为曲线元起点的坐标，αi为曲线元上任意点i的切线方位角，其计算通式为： （2） 式中，αA为曲线元起点切线方位角，PA为曲线元起点曲率，PB为曲线元终点曲率，L为曲线元长度，l为任意点距曲线元起点的弧长。±表示曲线元的偏向，左偏取“-”，右偏取“+”。αA、αi单位均为度。 由于用式（2）代入式（1）所表示的积分是不可积的，故采用复化辛卜生公式给出其数值积分式为： （3） 式中，m为积分区间等分数n的一半，α2K-1为曲线元上n等分点（偶数点）处的切线方位角，α2K为曲线元上n/2等分点（奇数点）处的切线方位角。 式（3）即为求线路中线坐标的复化辛卜生公式，也就是通用公式，它具有如下特点：1、适用于直线、圆曲线及各种类型的缓和曲线（完整或不完整）；2、计算精度可人为控制，即随积分区间等分数n的取值不同，其计算精度也不同，可以根据实际情况选用合适的n值以保证计算精度和提高计算效率，3、计算方向是可逆的，即可以沿里程递增方向计算也可以沿里程递减方向计算。 2　线路边线点坐标的计算 用式（3）求得线路中线点坐标后，根据式（2）所求得的切线方位角，易得出求边线点坐标的计算公式： （4） 式中，S为所求边线点与中线点之间的距离，Q为边线点方向与线路切线方向的夹角。±表示边线点相对于中线的方向，左边线点取“-”，右边线点取“+”。 3　地面点与线路中线相对位置关系的公式推导 3.1　计算原理 从线路中线坐标的通用公式推导可知，求地面点与线路中线相对位置关系即为求地面点与线路中线曲线元的相对位置关系。 如图1，AB为一曲线元，点P为一地面点，点P的坐标为已知或可用仪器测得。过点P作垂线PP′，lP为垂足点P′到起点A之间的弧长，DP为点P与线路中线的垂距，由图可知，点P与线路中线相对位置可由2个唯一的量lP和DP确定。由于点P可能位于中线的左边或右边，所以为了区分，DP也分正负，当P点位于线路左边时取负值，反之取正值。 图1 地面点与曲线元相对位置关系 由此，根据DP的符号，可以确定点P相对于线路的边向，而DP的值则确定了点P与中线的垂距，点P相对于线路的里程LP可由曲线元起点的里程LA和弧长lP计算得： LP=LA+lP 因此，确定地面点与线路中线的相对位置关系就转化为求DP和lP了。 直接求解DP和lP十分困难，在此采用趋近原理来达到计算目的。 如图2，已知曲线元起点A的坐标（XA，YA），过点A的切线方位角αA，起点A和终点B的曲率KA和KB，曲线元的长度L，地面点P的坐标（XP，YP），现求lP 和DP。 图2 DP与LP的趋近计算原理 在lP 和DP两个量中，如果求得lP ，即可用式（3）求得P′的坐标，由此可方便求出DP。因此问题的关键在于求lP ，在此趋近也就是指lP 趋近。由图可知，曲线元的法线有无穷条，但通过点P的法线却只有一条（点P为圆心时除外），即PP′，因此，只要求得P点到曲线元的某一法线的距离为0，则问题得解。 取点P到起点A法线的垂距l1的绝对值作为lP的近似值代入式（3）可求得曲线元上一点P1，再以P点到通过P1点法线的垂距l2加上l1作为新的lP的近似值代入式（3）。为了确保循环收敛，l2的符号必须满足：当P1点位于P′点后（l1＜