5800正反算及隧道超欠挖资料

1、5800正反算及隧道超欠挖B任意坐标正反算及隧道超欠挖程序1-JIN YU(主程序)Lbl 4:"1 f ZS,2 f FS,3 f CQW"?NN=1=>Goto 1:N=2=>Goto 2:N=3=>Goto 3Lbl 1:“ K=” ？S“P=” ?Zrog “ SJ-PMT : Abs(S-O) f W: Prog "SUB1-ZS":“ X=" Locate4,4,X:"Y=" : Locate4,4,Y:F-90f F:S f K:Prog “ SJGC :“ H=” Locate4,4,H 丄

2、Goto 4Lbl 2:“ X=” ？X“ Y=”Y:Prog“SJPM” :XfI:YfJ: Prog "SUB2-FS" : O +Wf S:“ K=" S 丄“ P="ZJSf K:Prog “ SJGC':“ H=”HJGoto 4Lbl 3:“ X=” ?X “ Y=”Y:Prog“ SJPM :XfI:YfJ: Prog "SUB2-FS" : O+Wf S:“ K=" S J“ P="ZJSf K:Prog “ SJGC':“ H=”HJProg “ SCQW 丄Goto4SJ?-PM

3、 （子程序名-平面线形数据库）If S > 59227.681（线元起点里程）:Then 2599818.013 fU（线元起点 X 坐标）：496887.918 fV （线元起点 Y 坐标）：59227.681 -O （线元起点里程）：208028 ' 39 fG （线元起点方位角）：100 fH （线元长度）：1 X 1045f（线元起点曲率半径）：1 X 1045 fR （线元终点曲率半径）：0 fQ （线元左右偏标志：左-1右1）： IfEndIf S > 59327.681（线元起点里程）:Then 2599730.112 fU （线元起点 X 坐标）：49684

4、0.237 fV （线元起点 Y 坐标）：59327.681 -O （线元起点里程）：208028 ' 39 fG （线元起点方位角）：90 fH （线元长度）：1 X 1045fp线元起点曲率半径）：750fR （线兀终点曲率半径）：-1 fQ （线兀左右偏标志：左-1右1 ）： IfEndSJ-GC?(子程序名-竖数据库)If Sv变坡终点里程 And S浚坡起点里程：Then大里程坡度 F :小里程坡度B : 变坡点里程 0:变坡点高度 G :半径R : Prog “SUB3GC :lfEnd(注：1.有多个竖曲线，依照上面的依次变更，每多一个，就增加一个。每次只需要 修改以上的

5、数据或增加一个判断，子程序不用变动。2.如整条线只有一个纵坡比喻为2%,那么程序应为If Sv终点里程 And S刁起点里程：Then 0.02 A : 0.02 B : 起点里程0:起点高程G : 1X 1045R: Prog “SUB3GC :IfEnd)SJ-CQW (子程序名 -隧道数据库 )“H1=”?F: ”P=”?ZIf F-H > 5.845:Then V(Abs (Z)2+(F-H-0.45)2 )-7.45 W: IfEnd:If F-H< 5.845 And F-H >-1.169: Then V(Abs (Z)-1)2+( Abs (F-H-1.5)

6、2)-6W : IfEnd:If F-HW-1.169: Then V(Abs (Z) 2+(H-F+14.05) 2)-16.5 W: IfEnd:“+CQ-,CQ=”: WSUB1-ZS(正算子程序)?1- P C : (P-R) - (2HPR) D : 180 -nE : 0.1739274226 A : 0.3260725774 B: 0.0694318442 K: 0.3300094782 L: 1-L F: 1-K M : U+W(A cos(G+QEKW(C+KWD)+Bcos(G+QELW(C+LWD)+Bcos(G+QEFW(C+FWD)+Aco s(G+QEMW(C+MW

7、D) :X：V+W(Asin(G+QEKW(C+KWD)+Bsin(G+QELW(C+ LWD)+Bsin(G+QEFW(C+FWD)+Asin(G+QEMW(C+MWD):Y： G+QEW(C+WD)+90:F： X+Zcos(F) :X： Y+Zsin(F) Y:SUB2-FS(反算子程序)?G-90T： Abs(Y-V)cos(T)-(X-U)sin(T) :W： 0:Z：Lbl 0：Prog "SUB1-ZS": T+QEW(C+WD) L :(J-Y)cos(L)-(I- X)sin(L) Z: If Abs(Z)<1E-6(注：此处 1E-6 可输入 0.

8、000001): then Goto1: Else W+Z W: Goto 0: IfEndLbl 1: Of Z: Prog "SUB1-ZS" :(J-Y) - sin(F) 吃SUB3-GC (高程子程序)O-Sf L:A -Bf W:Abs(R*W 2) f T:O-Tf M: O +TfPIf S< M:ThenGL*Bf H:Goto5:Else If SO:ThenGoto3:Else If S< P:ThenGoto4:lfEnd: IfEnd: IfEndLbi3If W > O : Then G+(M-S)2 -2R-L*Bf H:Go

9、to5:Else If W v O: Then G-(M-S)2 吃申-L*B f H:Goto5: IfEn d: IfE ndLbi4If W > O : Then G+( S-P)2 乞申-L*Af H:Goto5:Else If W v O: Then G-( S-P)2 -2曲-L*Af H:Goto5: IfEnd: IfEndLbi5HfH: Return说明: （正算 1 秒，反算和超欠挖 5秒完成 ）所有程序名。不同线路，只需改动SJ-PM，SJ-GC，SJ-CQW 三个子程序里的内容，其它不变。1 为正算， 2为反算， 3 为超欠挖 .K= 正算时，输入所求点里程。反

10、算时得出里程结果X=，Y= 正算时得出结果。反算时输入实测坐标P= 正算时输入偏距。反算时得出偏距H= 正算时得出高程。反算时得出高程H1=超欠挖时输入实测高程得出超欠挖值+CQ,-CQ=注，超欠挖输入时显示里程，编距，设计高程，不需输入。只输入实测量高程H1谢谢发贴： +302 分 【 测量空间 】0复化辛卜生公式及公路中线坐标计算中 的应用时间:2010-06-05 09:12:33要收藏 投稿指南来源:本站 作者:未知我要投稿我在公路中线坐标计算中，我们通常采用切线支距公式来计算曲线上各点的坐标。 但当在不同的曲线上计算时就需用不同的计算公式，这为计算也带来不便。在设有缓和曲线的圆曲线半

11、径较小或是卵形曲线上的坐标计算时，如公式选用不当就会出现较大计算误差，即便是能对切线支距公式进行多项展开， 也会增加计算的 难度。而用复化辛卜生公式不仅能解决不同曲线线型或直线上的坐标计算问题， 而且用复化辛卜生公式计算完全是可逆的(即：可顺前进方向也可逆向计算)， 尤其在计算第二缓和曲线和卵形曲线时显得尤为方便。用辛卜生公式计算坐标的精度可由人为或程序自行判断， 其计算结果完全能保证 坐标计算的精度要求。因此，可以说复化辛卜生公式是一个计算公路中线坐标的 万能公式。下面本人就该公式在公路中线坐标计算中的具体应用进行实例解析。Ij-Ya+ I、复化辛卜生公式rrlnrl(sim+4Ssinak

12、- + 2Ssinak+sinao)k=OE公式(1)Xa曲线旋点x坐标Ya曲线危点湮标 為曲线元起点切线方位角-曲2品点处徊线方妙 曲线元kn等分点处的切线方位角哉曲线元t待求点处切线方位角式中:H= (Zi - Za) /na 产 H + *+pJ(Zi-Z型(公式2)P.=Pa +Pb'Pa(z ,-za)Zb-Za（公式3）Zi 待求点桩号Za 曲线兀起点桩号Zb 曲线兀终点桩号p A 曲线兀起点曲率P B曲线元终点曲率a i曲线上任意一点处切线方位角的计算方法有以下三种方法：1 利用公式（3）求得曲率代入公式（2）计算2利用曲线元上已知起点和终点曲率用内插法求得曲率代入公式（

13、2）计算3 利用切线角公式计算二、算例例：已知雅（安）攀（枝花）高速公路西昌西宁立交 A匝道一卵形曲线（卵形曲 线相关参数见图一，其计算略。），相关设计数据见下表。现用辛卜生公式来计 算卵形曲线中桩坐标。图一已知相关设计数据见下表:主点桩号坐标(m切线方位角(0 )XY0 ”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.5YH2A

14、K0+384.0329922.31610007.909337 04 54.2HZAK0+444.0329981.36310000.0000 00 00(一) 由+271.881推算Zi=+223.715的坐标，n取2等分用公式、公式(2)计算+247.798处曲线及方位角：P +247.798=1 宁 75+(1 - 50-1 - 75)(247.798 -271.881)- (223.715 -271.881)=0.01666666666666667a+247.798 =71° 24' 18.5 ”+(0.016666667+1 - 75)(247.798 -271.881

15、) X 180- n 十2=50° 42' 26.37”其它各点依次代入公式计算，结果见下表： 曲率及切线方位角计算表桩号n等分点处曲线曲率n等分点 处切线方 位角0 , ”2n等分点处曲线曲率2n等分点处切线方 位角0 , ”公式计算内插法计算公式计算内插法计算+271.8871 2410.013333333=1/7518.5+259.840.0150.01561 37052.2+247.790.016666660.0166666650 4287726.4+235.750.018333330.01833333338 387330.96+223.7150.02=1/5025

16、24 36切线方位角图示1将计算出的数据代入公式(1)求得+223.715中桩坐标如下：X=9880.438+(271.881- 223.715) - 2-6X (cos71 ° 24' 18.5 ” +4(cos61 ° 37'52.22 ” +cos38° 38' 0.96 ”)+2cos50° 42' 26.37” + cos25 ° 24' 35.99 ”)=9910.5975 (设计值：9910.603)Y=10100.904+(223.715- 271.881) - 2-6X (sin71 &

17、#176; 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等分用公式计算+247.798处曲线及方位角：P +247.798=1 宁 50+(1 - 75-1 - 50)(247.798 -223.715) - (271.881 -223.715)=.0

18、1666666666666667&+247.798 =205° 24' 33.6 ” +(0.016666667+1 - 50)(247.798 -223.715) X 180- n 十2=230° 42' 23.98 ”其它各点依次代入公式计算，结果见下表：曲率及切线方位角计算表桩号n等分点处曲线曲率n等分点 处切线方 位角0 , ”2n等分点处曲线曲率2n等分点处切线方 位角0 , ”公式计算内插法计算公式计算内插法计算+223.7150.02=1/50205 2433.6+235.750.018333330.018333333218 37733

19、58.6+247.790.016666660.01666666230 42866 二24.0+259.840.0150.015241 3749.8+271.8810.013333333=1/75251 2418.5切线方位角图示2X=9910.603+(271.881- 223.715) - 2-6X (cos205 ° 24' 33.6 ” +4(cos218 ° 37 '58.87” +cos241° 37' 49.83 ”)+2cos230° 42' 23.98 ” + cos251 ° 24' 1

20、6.11 ”)=9880.4431 (设计值：9880.438)Y=10136.791+(271.881- 223.715) - 2-6X (sin205 ° 24' 33.6 ” +4(sin218 ° 37 '58.87” +sin241 ° 37' 49.83 ”)+2sin230 ° 42' 23.98 ” + sin251 ° 24' 16.11 ”)=10100.9008 (设计值：10100.904)由上可知，利用复化辛卜生公式计算路线坐标时可顺向或逆向计算。(三) 前面正、反向的计算均是从

21、卵形曲线的部分推算的，现在我们从卵形曲线所在完整缓和曲线的起点HZ =+368.213点来推算终点+223.715和中间点 +271.881.2ai=aA-90L -( n Rb)公式(4)1.计算Zi=+223.715中桩坐标曲线元上等分点曲率及切线方位角计算表(n取4等分)计算计算+368.2130108 12 02.1+350.1510.00250.0025106 54106 5425.1625.22+332.0890.0050.005103 01103 0134.3334.58+314.0260.00750.007596 3396 3328.8328.6487 3087 30+295.

22、9640.010.0109.979.97+277.90210.01250.012575 5137.2275 5137.54+259.8400.0150.01561 3761 3750.5851.35II44 4844 48+241.7770.01750.0175148.25 |47.8+223.71525 2425 240.02=1/5033.5833.58切线方位角图示3X=9890.293+(368.213- 223.715) -4-6X (cos108 ° 12' 02.1” +4(cos106 ° 54 '25.22 ” +cos96° 3

23、3' 28.64 ”+ cos75 ° 51' 37.54 ” + cos44 ° 48' 47.8 ”)+2(cos103 ° 01' 34.58 ” + cos87 ° 30' 9.97” +cos61 ° 37' 51.35 ”)+ cos25 ° 24' 33.58 ”)=9910.5963 (设计值：9910.603)Y=10006.838+(368.213- 223.715) -4-6X (sin108 ° 12' 02.1” +4(sin 106

24、° 54' 25.22 ” +si n96° 33' 28.64 ”+ sin75 ° 51' 37.54 ” + sin44 ° 48' 47.8 ”)+2(sin 103 ° 01' 34.58 ” +sin87 ° 30' 9.97” +sin61 ° 37' 51.35 ”)+ sin25 ° 24' 33.58 ”)=10136.7925 (设计值：10136.791)2.计算zi=+271.881中桩坐标曲线元上等分点曲率及切线方位角计算表

25、(n取4等分)桩号n等分点处曲线曲率n等分点处切线方位角2n等分点处曲 线曲率2n等分点处切线方位角0 , ”内插法计算公式 算)计公式(4)计算内插法计算公式 算)计公式(4)计算+368.2130108 12 02.1+356.1720.001666667107 3732.49+344.1300.003333333105 5402.97+332.0890.005103 0134.58+320.0470.00666666799 0005.6+308.0060.00833333393 4938.42+295.9640.01187 3009.97+283.9230.01166666780 014

26、4.02+271.8810.013333333=1/7571 2416.09+223.7150.02=1/5025 2425 24切线方位角图示4X=9890.293+(368.213- 271.881) -4-6X (cos108 ° 12' 2.1” +4(cos107° 37' 32.49 ” +cos103° 01' 34.58 ”+ cos93 ° 49' 38.42 ” + cos80 ° 01' 44.02 ”)+2(cos105 ° 54' 02.97 ” + cos99

27、 ° 00' 05.6 ” +cos87 ° 30' 09.97 ”)+ cos71 ° 24' 16.09 ”)=9880.4416 (设计值：9880.438)Y=10006.838+(368.213- 271.881) -4-6X (sin 108 ° 12' 2.1” +4(sin 107 ° 37' 32.49 ” +sin 103° 01' 34.58 ”+ sin93 ° 49' 38.42 ” + sin80 ° 01' 44.02 ”

28、)+2(sin105 ° 54' 02.97” + sin99 ° 00' 05.6 ” +sin87 ° 30' 09.97 ”)+ sin71 ° 24' 16.09 ”)=10100.9008 (设计值：10100.904)由此我们不难看出，利用复化辛卜生公式计算公路中桩坐标不仅顺逆向很方便, 而且可从中任取一段进行计算。而计算中的关键是计算各等分点处的切线方位 角，同时应分清 n 与 2n 等分点，注意不要混淆。对于 n 的取值，从上我们可以 看出：对一般缓和曲线取 2-3 即可，对卵形曲线视推算基准点而言不少于

29、4。其它如圆曲线、直线等中桩坐标的计算方法是一样的，在此就不一一简述。5800计算器全线坐标计算放样程序（修 改版）时间:2010-06-19 10:23:30来源：本站 作者：简单de回忆 我要投稿我要收藏投稿指南“XLZBJSC”丄（第一个程序）Lb1 0 /Cls : Fix 4 : 30f Dimz /“XHS="?G （后视点 X） :"YHS="?L （后视点 Y） :"XZJ="?M （ 置镜点X） :"YZJ="?N （置镜点 Y） :Pol（G-M,L-N）:"DH=":l （ 后视距）

30、丄J<0=>J+36（f J:"FH=":J ?DMS丄（后视方位角）Lbl 1 / （ If（如果的意思） And （和字的意思）The n （然后的意思）“K二”？K丄（计算里程）（下面有色的是数据库）If （如果）K< 51760.052本曲线 缓直点 桩号And （和）K> 51048.785上 一个曲线缓直点桩号：Then （然后）本曲线缓直点桩号51760.052f Z1: 上一个曲线缓直点桩号51048.785 f Z2 : 1 - 1 O （注： 左 偏曲线输 入-1f O右偏曲线输入1f O ）:偏角12 ' 23 '

31、; 19.5 ' f A : 半径3289.486 f R :第一缓和曲线 长度0 fZ6: 第二缓和曲线 长度0 f Z7:交点X坐标 f B :交点丫坐标 f C :小里程向交点 方位角 f E :交点向大里程 方位角 fF : Goto 2 : IfEnd / （重兰字是输入的数字） （曲线段分段输入）补充直线段输入如下If （如果）K<本段直线终点里程And （和）K>本段直线起点里 程:Then （然后）1f O本段直线终点里程f Z3:终点坐标Xf Z16: 终点 坐标 YfZ17: 方位角 f E:Goto 4:IfEnd /Lb1 2 / （曲线要素计算）Z

32、6 - 2- Z6A 3 - （240RA2）+ Z6A 5- （34560*RA4） f Z8 / （ M1不输）Z7 - 2- Z7a 3- （240RA2）+ Z7a 5- （34560 RA4） f Z9 / （M2不输）Z6a 2- （24R） - Z6a 4-（2688RA3） f Z10 /（ P1 不输）Z7a 2- （24R） - Z7a 4-（2688RA3） f Z11 /（ P2 不输）（n AR ） - 180+0. 5 X （ Z6+ Z7） f S / （曲线总长）90 X Z6 - （ R X n ） f Z14 / （第一缓和曲线总偏角）（X十以后自己改）9

33、0 X Z7 - （ R X n ） f Z15 / （第二缓和曲线总偏角,可以省略）Z8+（ R+Z10）TAN（A/2）-（Z10-Z11 ）/SIN Af Z12/ （切线 T1）Z9+（ R+Z11）TAN（A/2）+（Z10-Z11 ）/SIN Af Z13/ （切线 T2）B+ Z12*COS （E+180） f Z16 /（ ZH 点 X ）C+ Z12*SIN（E+180） f Z17/（ ZH 点 Y ）Z1- Sf Z3 / （ZH 点里程）Z3+ Z6 -Z4 / (HY 点里程) Z1-Z7 -Z5 / (YH 点里程)GOTO 3/LB1 3 / (判断里程点与曲线关

34、系)If K < Z3 And K> Z2 : Then Goto 4 : IfEnd/If K < Z4 A ndK> Z3 : T henGoto 5:IfEnd/If K < Z5 A ndK> Z4 : T henGoto 6:IfEnd/If K < Z1 A ndK> Z5 : T henGoto 7:IfEnd/LB1 4 / (里程小于直缓点直线独立坐标)K- Z3 X : 0 Y : E T : Prog“TYZBC” ： Goto 1 /Lb1 5 / (第一缓和曲线独立坐标)K- Z3 H /H-H A5/(40*RA2* Z6 A2)+H A9/(3456*RA4* Z6A4) X /屮3/(6*R* Z6)-屮7/(336*只八3* Z6F3) Y /90*屮2/( R* n * Z6) T /IF O >0 : Then T +E T : Else E -T T : T<0=>360+T T : IfEnd /PRrog “TYZ