复化辛普森公式应用.doc

在公路中线坐标计算中，我们通常采用切线支距公式来计算曲线上各点的坐标。但当在不同的曲线上计算时就需用不同的计算公式，这为计算也带来不便。在设有缓和曲线的圆曲线半径较小或是卵形曲线上的坐标计算时，如公式选用不当就会出现较大计算误差，即便是能对切线支距公式进行多项展开，也会增加计算的难度。而用复化辛卜生公式不仅能解决不同曲线线型或直线上的坐标计算问题，而且用复化辛卜生公式计算完全是可逆的（即：可顺前进方向也可逆向计算），尤其在计算第二缓和曲线和卵形曲线时显得尤为方便。用辛卜生公式计算坐标的精度可由人为或程序自行判断，其计算结果完全能保证坐标计算的精度要求。因此，可以说复化辛卜生公式是一个计算公路中线坐标的万能公式。下面本人就该公式在公路中线坐标计算中的具体应用进行实例解析。一、复化辛卜生公式式中：h=（ziza）/n(公式2)(公式3)zi 待求点桩号za曲线元起点桩号zb曲线元终点桩号a曲线元起点曲率b曲线元终点曲率ai曲线上任意一点处切线方位角的计算方法有以下三种方法：1利用公式（3）求得曲率代入公式（2）计算2利用曲线元上已知起点和终点曲率用内插法求得曲率代入公式（2）计算3利用切线角公式计算二、算例例：已知雅（安）攀（枝花）高速公路西昌西宁立交a匝道一卵形曲线（卵形曲线相关参数见图一，其计算略。），相关设计数据见下表。现用辛卜生公式来计算卵形曲线中桩坐标。图一已知相关设计数据见下表：主点桩号坐 标（m）切线方位角（）xy ”zhak0+0909987.40310059.37892 17 26.2hy1ak0+1609968.98110125.341132 23 51.6yh1ak0+223.7159910.60310136.791205 24 33.6hy2ak0+271.8819880.43810100.904251 24 18.5yh2ak0+384.0329922.31610007.909337 04 54.2hzak0+444.0329981.36310000.0000 00 00（一）由+271.881推算zi=+223.715的坐标，n取2等分用公式(3)、公式（2）计算+247.798处曲线及方位角:+247.798=175+(150-175)(247.798-271.881) (223.715-271.881)=0.01666666666666667a+247.798=712418.5” +(0.016666667+175)(247.798-271.881)1802=504226.37”其它各点依次代入公式计算，结果见下表:曲率及切线方位角计算表桩号n等分点处曲线曲率n等分点处切线方位角 ”2n等分点处曲线曲率2n等分点处切线方位角 ”公式(3)计 算内插法计 算公式(3)计 算内插法计 算+271.8810.013333333=1/7571 24 18.5+259.8400.0150.01561 37 52.2+247.7980.0166666670.01666666750 42 26.4+235.7570.0183333330.018333333338 38 0.96+223.7150.02=1/5025 24 36切线方位角图示1将计算出的数据代入公式（1）求得+223.715中桩坐标如下：x=9880.438+(271.881-223.715)26(cos712418.5”+4(cos613752.22”+cos38380.96”)+2cos504226.37”+ cos252435.99”)=9910.5975(设计值：9910.603)y=10100.904+(223.715-271.881)26(sin712418.5”+4(sin613752.22”+sin38380.96”)+2sin504226.37”+ sin252435.99”)=10136.7945(设计值：10136.791)（二）由+223.715推算zi=+271.881的坐标，n取2等分用公式(3)计算+247.798处曲线及方位角:+247.798=150+(175-150)(247.798-223.715)(271.881-223.715)=.01666666666666667a+247.798=2052433.6”+ (0.016666667+150)(247.798-223.715)1802=2304223.98”其它各点依次代入公式计算，结果见下表:曲率及切线方位角计算表桩号n等分点处曲线曲率n等分点处切线方位角 ”2n等分点处曲线曲率2n等分点处切线方位角 ”公式(3)计 算内插法计 算公式(3)计 算内插法计 算+223.7150.02=1/50205 24 33.6+235.7570.0183333330.0183333333218 37 58.6+247.7980.0166666660.016666666230 42 24.0+259.8400.0150.015241 37 49.8+271.8810.013333333=1/75251 24 18.5切线方位角图示2x=9910.603+(271.881-223.715)26(cos2052433.6”+4(cos2183758.87”+cos2413749.83”)+2cos2304223.98”+ cos2512416.11”)=9880.4431(设计值：9880.438)y=10136.791+(271.881-223.715)26(sin2052433.6”+4(sin2183758.87”+sin2413749.83”)+2sin2304223.98”+ sin2512416.11”)=10100.9008(设计值：10100.904)由上可知，利用复化辛卜生公式计算路线坐标时可顺向或逆向计算。（三）前面正、反向的计算均是从卵形曲线的部分推算的，现在我们从卵形曲线所在完整缓和曲线的起点hz=+368.213点来推算终点+223.715和中间点+271.881.ai=aa-90l2(rls)公式(4)1.计算zi=+223.715中桩坐标曲线元上等分点曲率及切线方位角计算表(n取4等分)桩 号n等分点处曲线曲率n等分点处切线方位角 ”2n等分点处曲线曲率2n等分点处切线方位角 ”公式(3)计 算内插法计 算公式(2)计算公式(4)计算公式(3)计 算内插法计 算公式(2)计算公式(4)计算+368.2130108 12 02.1+350.1510.00250.0025106 54 25.16106 54 25.22+332.0890.0050.005103 01 34.33103 01 34.58+314.0260.00750.007596 33 28.8396 33 28.64+295.9640.010.0187 30 09.9787 30 9.97+277.9020.01250.012575 51 37.2275 51 37.54+259.8400.0150.01561 37 50.5861 37 51.35+241.7770.01750.017544 48 48.2544 48 47.8+223.7150.02=1/5025 24 33.5825 24 33.58切线方位角图示3x=9890.293+(368.213-223.715)46(cos1081202.1”+4(cos1065425.22”+cos963328.64”+ cos755137.54”+ cos444847.8”)+2(cos1030134.58”+ cos87309.97” +cos613751.35”)+ cos252433.58”)=9910.5963(设计值：9910.603)y=10006.838+ (368.213-223.715)46(sin1081202.1”+4(sin1065425.22”+sin963328.64”+ sin755137.54”+ sin444847.8”)+2(sin1030134.58”+ sin87309.97” +sin613751.35”)+ sin252433.58”)=10136.7925(设计值：10136.791)2.计算zi=+271.881中桩坐标曲线元上等分点曲率及切线方位角计算表(n取4等分)桩号n等分点处曲线曲率n等分点处切线方位角 ”2n等分点处曲线曲率2n等分点处切线方位角 ”内插法计 算公式(2)计算公式(4)计算内插法计 算公式(2)计算公式(4)计算+368.2130108 12 02.1+356.1720.001666667107 37 32.49+344.1300.003333333105 54 02.97+332.0890.005103 01 34.58+320.0470.00666666799 00 05.6+308.0060.00833333393 49 38.42+295.9640.0187 30 09.97+283.9230.01166666780 01 44.02+271.8810.013333333=1/7571 24 16.09+223.7150.02=1/5025 24 33.5825 24 33.58切线方位角图示4x=9890.293+(368.213-271.881)46(cos108122.1”+4(cos1073732.49”+cos1030134.58”+ cos934938.42”+ cos800144.02”)+2(cos1055402.97”+ cos990005.6” +cos873009.97”)+ cos712416.09”)=9880.4416(设计值：9880.438)y=10006.838+ (368.213-271.881)46(sin108122.1”+4(sin1073732.49”+sin1030134.58”+ sin934938.42”+ sin800144.02”)+2(sin1055402.97”+ sin990005.6” +sin8730