关于利用线元法计算匝道坐标及excel辅助计算的可行性探讨.doc

关于利用线元法计算匝道坐标及 excel辅助计算的可行性探讨 范利民 （江苏恒基路桥有限公司） 随着我国公路建设的发展，公路工程道路线型呈现出多样化的特点，尤其是互通立交匝道平面线型的各种线元使用灵活且复杂，主要使用直线、缓和曲线、圆曲线线型变化组合而成，它的平曲线要素成果不像道路主线那样标准的“直曲平面曲线要素表”，而且多采用小半径曲线，因此它的中边桩坐标计算也不宜采用标准的交点转角法来计算。本文介绍采用线元法计算则会使计算工作变得方便简洁。 一、线元的概念 匝道平曲线的线型多是利用直线、缓和曲线、圆曲线这几种基本线型经过所需变化相互连接组合，就像堆积木一样一个接一个拼接累加而成，我们把每个组成模块称为线元。线元法计算坐标的基本原理就是从曲线的起点开始，利用起终点的已知线元参数，按照各线元连接顺序逐步计算出整条线路的中边桩坐标。此处所需的已知线元参数有，起终点的中桩坐标、起点的走向方位角、起终点半径、起终点桩号以及曲线的转角方向（左转或右转）。 二、线路坐标计算 1、首先我们要弄清楚回旋线上任意点的走向方位角。如图1 设：回旋线起点A的曲率为ρA，其里程为DKA;回旋曲线终点B的曲率为ρB，其里程为DKB；AX’Y’为以A点为原点，以A点的切线为X’轴的局部坐标系，AXY为线路坐标系。因为回旋线上任意点的曲率半径Ri和该点距离起点的距离li成反比。 故此，任意点曲率ρi=1/Ri=li/C(其中C为常数)，为线性变化。 由此可知：里程为DKi处的曲率-------(1)，（其中，A、B的曲率路线右转时为正值，左转时为负值） 由图1可知： ---------(2) 将（1）代入（2）可得--------(3) 若A点处在线路坐标中的走向方位角为，则DKi处切线方位角 -------（4） 将（3）代入（4）可得：-----(5) 2、计算局部坐标系下DKi处的坐标 由图1可知-------（6），在区间DKA和DKi之间对式（6）进行定积分即 利用辛普森公式展开得，为公式简洁作如下定义： 令:； 式（7）中：；；；（当E360时E+360。当E360时E-360） ；； 则有： ---------(7) 至此,若A点在线路坐标系中的坐标为XA、YA。则可得里程DKi处的线路坐标 ； 边桩坐标式中：L为边桩距中桩的距离，Z为符号常数，当所求为左边桩时Z=-1；所求为右边桩时Z=1。 由此，可依次解算出整条线路所需点位的坐标值，由于手工计算工作量甚大，可考虑采用软件辅助计算。 三、线元法在EXCEL工作表中的应用 借助EXCEL工作表的超强函数功能，我们可以将线元法坐标的计算在实践中加以推广应用，帮助我们解决工作中的实际问题。以吉莲高速AP1合同段敖城互通A匝道坐标计算为例，其部分线元要素如下表： 要素桩 桩号 X Y 方位角 长度 线元参数 ZH 　 　 　 　 　 A=120 HY AK1+657.954 2984058.147 514401.51 187-52-57 　 66.325 R=左400 YH AK1+724.279 2983991.997 514397.89 178-22-55.7 36 A=120 HZ AK1+760.279 2983956.062 514399.985 175-48-13.8 　 　 将线元要素输入自编线元法坐标计算简易程序EXCEL表格如下： 输入起点X（A) 输入起点Y（B) 起点小数度方位角( C) 起点曲率半径（D) 终点曲率半径（E) 2984058.147 514401.51 187.8825 -0.0025 -0.0025 输入起点桩号(F) 输入终点桩号（G) 输入右边桩距离（W1) 输入左边桩距离（W2) 边桩夹角(小数度） 1657.954 1724.279 3 -3 90 输入边桩夹角（度分秒） 输入起点半径 输入终点半径 转向 起点方位角度分秒 90 400 400 左 187.5257 说明：如果左边桩输入负值，右边桩输入正值。 程序内部函数运算所需点位坐标结果如下： 里程(H) X（中桩） Y（中桩） X(右边桩） Y(右边桩） 1657.954 2984058.147 514401.510 2984058.558 514398.538 1660 2984056.120 514401.235 2984056.516 514398.261 1670 2984046.192 514400.038 2984046.514 514397.055 1680 2984036.237 514399.090 2984036.484 514396.100 1690