辛普生法则在公路匝道中资料.ppt

辛普生法则在公路匝道中及边桩坐标计算中的应用 目前在公路匝道中、边桩计算中采用的计算公式较多，但大部分采用泰勒级数展开公式进行计算。当在不同曲线元上计算就需要用不同的计算公式，这给使用普通计算器计算带来很大不便。当遇有卵形曲线或缓和曲线的圆曲线半径较小时，往往因公式选用不当带来计算误差，虽然泰勒级数展开公式可以进行多项次展开，但公式过于繁杂，同时也会增加计算难度。而采用辛普生法则进行计算就非常简便，该公式可以适用于任意的直线段、圆曲线段、缓和曲线段上中、边桩，并且其计算是可逆的。下面就该法则在公路工程测量计算中的具体应用进行解析。 一、曲线元上任意点的切线方位角推导 曲线平面线形组合大致分基本型、S型、卵形、凸形、复合型和C型等，无论哪种组合型式，其平面线形计算的复杂性均主要体现在缓和曲线的计算上。1、卵形曲线中的缓和曲线 如图（1），在半径分别为R1、R2（R1 R2）的两段圆曲线之间插入缓和曲线段AQB，·设曲线MAQB为一段缓和曲线。M、A、B的半径分别为∞、R1、R2，切线分别为ME、AJ、BF；CAQ1是半径为R1的圆弧段并在A点与缓和曲线公切，Q2BD是半径为R2的圆弧段并在B点与缓和曲线公切。令MA=L1，MB=L2， LS＝L1－L2。缓和曲线参数为：A2=R1R2 LS /（R1-R2） 公式（1） 则有：L1= A2/ R1= R2LS /（R1－R2）； L2= A2/ R2= R1LS /（R1－R2）； βA=L12/2A2＝R2LS/（2R1（R1－R2））；βB＝L22/2A2＝R1LS/（2R2（R1－R2））；可以推出：α=βB -βA ＝(1/R1+1/R2)LS/2 公式（2） （α：缓和曲线AB段转向角。 对于曲线元上任意点P，P点处曲率半径为1/Rp，该点处转向角为αP ，曲线元起点A到P点的弧长为L，根据式（2）有： αP＝（1/R1+1/RP) L /2 (3)由式（1）得：（R1L＋A2）RP＝R1A2，可以求得1/RP＝1/R1＋L/A2；将1/RP代入式（3）得： αP＝L/R1＋L2（1/R2－1/R1）/（2LS） 公式 (4)故缓和曲线上任意一点P的切线方位角：TP＝TA±αP＝TA±（L/R1＋L2（1/R2－1/R1）/（2LS）） 公式(5) 2、其它曲线元 对于基本型缓和曲线，R1＝∞（或R2＝∞）代入式（4）得：αP＝L2/（2R2LS）；与我们常用公式一样。对于圆曲线段，R1＝R2，代入式（4）得：αP＝L/R1；与我们常用公式一样。对于直线段，R1＝ R2＝∞，代入式（4）得：αP＝0。 因此可以得出这样的结论：任何曲线元上任意一点的切线方位角均可以用式（5）进行计算。 二、中、边桩坐标计算 1、使用辛普生积分法则进行计算 方位角TP＝TA±αP＝TA＋K(L/R1＋L2（1/R2－1/R1）/（2LS）)； (方位角TP计算结果为弧度；K为判断符，左偏取－1，右偏取＋1)。中桩坐标为：XP＝XA＋∫0L cos（T）dx 公式（6）YP＝YA＋∫0L sin（T）dx 公式（7）casio fx-4800p等计算器具有积分功能，其计算采用了辛普生法则进行积分计算，式（6）和式（7）转换成计算器计算模型为：XP＝XA＋∫（cos(TA＋K(L/R1＋L2(1/R2-1/R1)/(2LS))，0，L,n)；YP＝YA＋∫（sin(TA＋K(L/R1＋L2(1/R2-1/R1)/(2LS))，0，L,n)。式中n为1～9的整数，积分分割数N=2n，其他符号意义同上述。 值得注意的是，按此公式在计算器中编计算程序时，各种计算器将使用某些变量储存数据，例如casio fx-4800p使用(K,L,M,N)分别储存积分区间、分割数、积分值，这几个变量在编程中不能使用；另外计算时要将计算器角度变换到弧度状态。本公式可以适用于直线、圆曲线、缓和曲线等一切曲线元计算，是一项综合性公式。2、casio fx-4800p程序程序名：L-R （中、边桩坐标及切线方位角计算） Prog “Y-S”Q=Z[10(T－1)+1]E=Z[10(T－1)+2] F=Z[10(T－1)+3]A=Z[10(T－1)+4] C=Z[10(T－1)+5] J=Z[10(T－1)+6] O=Z[10(T－1)+7]  S=Z[10(T－1)+8]P=AbsH－QP0 == P=P≠== Goto 2: ◣Lbl 2PS == P=P◢≠