CASIO5800P计算器正反算及对称隧道超欠挖.doc

任意线路坐标正反算及对称隧道超欠挖程序由于现在计算机普及，计算机功能日益强大，宜采用较简单的积分公式，便于计算机处理。单线元通用积分公式如下M = (1.0/Re-1.0/Rs)/Ls;x=cos(Ta + L/Rs + 0.5*M *L*L)，0，L;y=sin(Ta +L/Rs + 0.5*M *L*L)，0，L;a(i)= Ta +L/Rs + 0.5*M *L*LRs:缓和曲线起点半径Re:缓和曲线止点半径Rs，Re （NE坐标系下，右偏为正，左偏为负）Ta:缓和曲线起点的真北方位角Ls:不完整缓和曲线长度。此公式为缓和曲线在坐标系下任意位置的通用积分公式，能完全适应缓和曲线左偏、右偏、Rs Re 、Rs NE”:“2.NE=SZ”:“1,2=”？Q：Lbl4:“LICHENG=”？S:Prog“QXJSSUB0”Lbl0:If Q=1:Then Goto1:IfEnd：If Q=2：Then Goto2：IfEnd（1为正算，2为反算及超欠挖）Lbl 1: “OUT JL=”? Z：Prog “QXJSSUB1”： “XZB=”：NN“YZB=”:EE“QXFWJ=”:FF:FDMSProg “SJ-GC”： “H=”：HGoto 4(正算-输入待求点里程LICHENG=、输入待求点偏距OUT JL=、显示待求点里程XZB=、显示待求点里程YZB=、显示切线方位角QXFWJ=、显示待求点里程设计高程H=)Lbl 2:“XZB=”?B：“YZB=” ?C: B N： CE：Prog QXJSSUB2 ： “LICHENG=： S S“OUT JL=:ZZProg“SJ-GC”：“H=”:HProg“QXJS-CQW”Goto 4(反算-输入近似里程LICHENG=、输入实测点XZB=、 输入实测点YZB=、显示实测点里程LICHENG=、显示实测点偏距OUT JL=、显示实测点的设计高程H=、输入实测高程HZB=、显示超欠挖值及第几圆心上CQY？=)QXJSSUB0（子程序名-平面线形数据库）If S 49056.209（线元终点里程）:Then 2863794.951U（线元起点X坐标）：484174.994V（线元起点Y坐标）：48908.857O（线元起点里程）：21902004”A（线元起点方位角）：147.352 L（线元长度）：650P（线元起点曲率半径）：1000R（线元终点曲率半径）：Return：IfEnd(注：有多个平面线形，依照上面的依次变更, 每多一个，就增加一个。每次只需要修改以上的数据或增加一个判断。SJ-GC（子程序名-竖曲线数据库）If S变坡终点里程（下一变坡点的起点里程）And S（这一变坡点起点里程）：Then大里程坡度A：小里程坡度B：变坡点里程O：变坡点高度G：半径R：Prog“SUB3-GC”:IfEnd(注：1.有多个竖曲线，依照上面的依次变更，每多一个，就增加一个。每次只需要修改以上的数据或增加一个判断。2.如整条线只有一个纵坡比喻为2，那么程序应为If S终点里程And S起点里程：Then 0.02A：0.02B：起点里程O：起点高程G：11045R：Prog“SUB3-GC”:IfEnd)QXJS-CQW (子程序名-隧道数据库)“HZB=”?FIf F-H5.845（为圆心1分段高度）:Then (Abs (Z)2+(F-H-0.45（圆心1高度）)2 )-7.45（圆心1半径）W: “CQY1=”:W: IfEnd:(显示超欠和第1圆心)If F-H5.845 And F-H-1.169: Then (Abs (Z)-1)2+( Abs (F-H-1.5) 2)-6W : “CQY2=”:W:IfEnd: (显示超欠和第2圆心)If F-H-1.169: Then (Abs (Z) 2+(H-F+14.05) 2)-16.5W: “CQY3=”:W: IfEnd(显示超欠和第3圆心)(注：隧道断面有多个圆心，依照上面的依次变更, 每多一个，就增加一个。图示中隧道中线与测设基线重合所以为ABS(Z)，如果不重合就是右线：ABS(4.375-Z),左线：ABS(4.375+Z)，见下图:QXJSSUB1(正算子程序)0.5（1R-1P）LD:S- OXU+(cos(A+（XP+DX2）180,0，X)NV+(sin(A+(XP+DX2)180,0，X)EA+(XP+DX2) 180FN+Zcos（F+90）N:E+Zsin(F+90) EReturnQXJSSUB2(反算子程序)Lbl 1:0Z：1Q：Prog“QXJSSUB0”: 0.5（1R-1P）LD:S-OXU+(cos(A+(XP+DX2)180,0,X)NV+(sin(A+(XP+DX2)180,0,X)EA+(XP+DX2)180FN+Zcos(F+90) N:E+Zsin(F+90) E:Pol(N-B+10(-46), E-C+10(-46):Isin(F-90-J) W:S+WSIfAbs(W)0.0001 :Then Goto1:IfEndLbl 2: 0Z：Prog“QXJSSUB1”:(C-E) sin(F+90) ZReturnSUB3-GC（高程子程序）O-SL:A-BW:Abs(RW2) T:O-TM:O+TPIf SM:Then G-LBH:Goto5:Else If SO:Then Goto3:Else If SP:Then Goto4:IfEnd: IfEnd: IfEndLbl3：If W0：Then G+(M-S)22R-LBH:Goto5: Else If W0: Then G-(M-S)22R-LBH:Goto5: IfEnd: IfEndLbl4：If W0：Then G+( S-P)22R-LAH:Goto5: Else If W0: Then G-( S-P)22R-LAH:Goto5: IfEnd: IfEndLbl5HH: Return说明：(正算1秒，反算和超欠挖5秒完成)所有程序名。不同线路，只需改动QXJSSUB0，SJ-GC，QXJS-CQW三个子程序里的内容，其它不变。1为正算，2为反算和超欠