复化辛卜生式在公路中线坐标计算中的应用1.doc

复化辛卜生（辛甫生）公式在公路中线坐标计算中的应用 【摘 要】在现有公路中线坐标计算的公式中，对于直线、圆曲线、缓和曲线(或卵形曲线)等中桩坐标的计算需用切线支距公式等不同的公式。而用复化辛卜生公式不仅能解决不同线形的坐标计算所需多个公式问题，而且计算的方向完全是可逆的，其计算精度也可根据需要自行判断决定。 【关键词】辛卜生公式 中线坐标计算 应用 在公路中线坐标计算中，我们通常采用切线支距公式等不同公式来计算直线或曲线上各点的坐标。但当计算不同的线形时就需用不同的计算公式，这就为计算带来不便。尤其是在设有较小圆曲线半径（即相对于半径越小而缓和曲线越长时）的缓和曲线（或卵形曲线）上的坐标计算时，如公式选用不当或切线支距公式所带项数较少时就会出现较大的计算误差，即便是能对切线支距公式进行多项展开，也会增加计算的难度。而用复化辛卜生公式不仅能解决不同曲线线型和直线上的坐标计算问题，而且用复化辛卜生公式计算完全是可逆的（即：可顺前进方向也可逆向计算），尤其在计算第二缓和曲线和卵形曲线时显得尤为方便。 用辛卜生公式计算坐标的精度可由人为或程序自行判断，其计算结果也完全能保证施工放样时坐标计算的精度要求。因此，可以说复化辛卜生公式是一个计算公路中线各种线形坐标的通用公式。下面笔者就该公式在公路中线坐标计算中的具体应用进行实例解析（本文以卵形曲线坐标的计算为例）。 一、复化辛卜生公式 X A — 曲线元起点X坐标 Y A — 曲线元起点Y坐标 a A — 曲线元起点切线方位角 a K+ — 曲线元上2n等分点处的切线方位角 a K — 曲线元上n等分点处的切线方位角 a i—曲线元上待求点处切线方位角，可用如下3种方法计算： 1．利用公式（3）求得曲率代入公式（2）计算 2．利用曲线元上已知推算起点和终点曲率用内插法求得曲率代入公式（2）计算 3．利用切线角公式计算 注1：aA为推算起点的切线方位角，以实际推算方向为前进方向，而不是指路线设计的前进方向。 式中： H=│Zi-ZA│÷n ai=aA±90(ρi +ρA)(Zi-ZA)÷π (公式2) 注2： 当桩号顺推算曲线元前进方向递增且为左转时取- 当桩号顺推算曲线元前进方向递增且为右转时取+ 当桩号顺推算曲线元前进方向递减且为左转时取+ 当桩号顺推算曲线元前进方向递减且为右转时取- ρi=ρA+(ρB-ρA) (Zi-ZA)÷(ZB-ZA) (公式3) Z i — 待求点桩号 Z A — 曲线元推算起点桩号 Z B — 曲线元终点桩号 Ρ i — 曲线元上待求点曲率 ρ A — 曲线元推算起点曲率 ρ B — 曲线元终点曲率 二、算例 例：已知雅（安）攀（枝花）高速公路西昌西宁立交A匝道一卵形曲线（相关参数见图一、图二、图三，其计算略。），相关设计数据见表一。现用辛卜生公式来计算卵形曲线中桩坐标：  图一 表一：卵形曲线相关设计数据 主点 桩号 坐 标 （m） 切线方位角 （θ） X Y ° ’ ” ZH AK0+090 9987.403 10059.378 92 17 26.2 HY1 AK0+160 9968.981 10125.341 132 23 51.6 YH1 AK0+223.715 9910.603 10136.791 205 24 33.6 HY2 AK0+271.881 9880.438 10100.904 251 24 18.5 YH2 AK0+384.032 9922.316 10007.909 337 04 54.2 HZ AK0+444.032 9981.363 10000.000 0 00 00 图二 图三  （一）由+271.881推算Zi=+223.715的坐标（逆向），n取2等分 用公式(3)、公式（2）计算+247.798处曲线曲率及切线方位角: ρ+247.798=1÷75+(1÷50-1÷75) ×(247.798-271.881) ÷(223.715-271.881) =0.016666667 a+247.798=71°24’18.5” +90×(0.016666667+1÷75) ×(247.798-271.881)÷π =50°42’26.37” 其它各点依次代入公式计算，结果见表二： 表二：曲率及切线方位角计算表 桩号 n等分点处曲线曲率 n等分点处 切线方位角 ° ’ ” 2n等分点处曲线曲率 2n等分点处 切线方位角 ° ’ ” 公式(3) 计 算 内插法 计 算 公式(3) 计 算 内插法 计 算 +271.881 0.013333333=1/75 71 24 18.5 +259.840 0.015 0.015 61