复化辛普森公式应用

《复化辛普森公式应用》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、在公路中线坐标计算中，我们通常采用切线支距公式来计算曲线上各点的坐标。但当在不同的曲线上计算时就需用不同的计算公式，这为计算也带来不便。在设有缓和曲线的圆曲线半径较小或是卵形曲线上的坐标计算时，如公式选用不当就会出现较大计算误差，即便是能对切线支距公式进行多项展开，也会增加计算的难度。而用复化辛卜生公式不仅能解决不同曲线线型或直线上的坐标计算问题，而且用复化辛卜生公式计算完全是可逆的（即：可顺前进方向也可逆向计算），尤其在计算第二缓和曲线和卵形曲线时显得尤为方便。用辛卜生公式计算坐标的精度可由人为或程

2、序自行判断，其计算结果完全能保证坐标计算的精度要求。因此，可以说复化辛卜生公式是一个计算公路中线坐标的万能公式。下面本人就该公式在公路中线坐标计算中的具体应用进行实例解析。一、复化辛卜生公式式中：H=（Zi－ZA）/n(公式2)(公式3)Zi—待求点桩号ZA—曲线元起点桩号ZB—曲线元终点桩号ρA—曲线元起点曲率ρB—曲线元终点曲率a i曲线上任意一点处切线方位角的计算方法有以下三种方法：1．利用公式（3）求得曲率代入公式（2）计算2．利用曲线元上已知起点和终点曲率用内插法求得曲率代入公式（2）计算3

3、．利用切线角公式计算二、算例例：已知雅（安）攀（枝花）高速公路西昌西宁立交A匝道一卵形曲线（卵形曲线相关参数见图一，其计算略。），相关设计数据见下表。现用辛卜生公式来计算卵形曲线中桩坐标。图一已知相关设计数据见下表：主点桩号坐标（m）切线方位角（θ）XY°’”ZHAK0+0909987.40310059.378921726.2HY1AK0+1609968.98110125.3411322351.6YH1AK0+223.7159910.60310136.7912052433.6HY2AK0+271.88

4、19880.43810100.9042512418.5YH2AK0+384.0329922.31610007.9093370454.2HZAK0+444.0329981.36310000.00000000（一）由+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.01666666666666667a+247.

5、798=71°24’18.5”+(0.016666667+1÷75)(247.798-271.881)×180÷π÷2=50°42’26.37”其它各点依次代入公式计算，结果见下表:曲率及切线方位角计算表桩号n等分点处曲线曲率n等分点处切线方位角°’”2n等分点处曲线曲率2n等分点处切线方位角°’”公式(3)计算内插法计算公式(3)计算内插法计算+271.8810.013333333=1/75712418.5　　+259.840　　0.0150.015613752.2+247.7980.0166666

6、670.016666667504226.4　　+235.757　　0.0183333330.018333333338380.96+223.7150.02=1/50252436　　切线方位角图示1将计算出的数据代入公式（1）求得+223.715中桩坐标如下：X=9880.438+(271.881-223.715)÷2÷6×(cos71°24’18.5”+4(cos61°37’52.22”+cos38°38’0.96”)+2cos50°42’26.37”+cos25°24’35.99”)=9910.597

7、5 (设计值：9910.603)Y=10100.904+(223.715-271.881)÷2÷6×(sin71°24’18.5”+4(sin61°37’52.22”+sin38°38’0.96”)+2sin50°42’26.37”+sin25°24’35.99”)=10136.7945 (设计值：10136.791)（二）由+223.715推算Zi=+271.881的坐标，n取2等分用公式(3)计算+247.798处曲线及方位角:ρ+247.798=1÷50+(1÷75-1÷50)(247.798-

8、223.715)÷(271.881-223.715)=.01666666666666667a+247.798=205°24’33.6”+(0.016666667+1÷50)(247.798-223.715)×180÷π÷2=230°42’23.98”其它各点依次代入公式计算，结果见下表:曲率及切线方位角计算表桩号n等分点处曲线曲率n等分点处切线方位角°’”2n等分点处曲线曲率2n等分点处切线方位角°’”公式(3)计算内插法计算公式(3)计算内插法计算+2