任意坐标正反算及隧道超欠挖程序最终版.doc

任意坐标正反算及隧道超欠挖程序最终版(2013-6-1)说明：一直以来程序正反算没有问题，只是高程不对，那超欠挖就扯谈了，经修改实践后，本程序决对能用。此程序以后不会再有什么改动。如果您的计算器不能运行或是计算错误，请认真核对，因为每一个语句命令我都反复核对计算器，运行顺畅无误。1-JIN YU(主程序)Lbl 4: “1ZS,2FS,3CQW”?N(选择计算模式，1为正算，2为反算，3为超欠挖)N=1=Goto 1:N=2=Goto 2:N=3=Goto 3Lbl 1: “K=”? S：“P=”? Z：Prog “SJ-PM”： Abs(S-O) W：Prog “SUB1-ZS”： “X=”：Locate4,4,X: “Y=”：Locate4,4,Y:F-90F:S K:Prog“SJ-GC”：“H=”：Locate4,4,H：“W=”: Locate4,4,FGoto 4(正算-输入待求点里程K=、输入待求点偏距P=、显示待求点里程X=、显示待求点里程Y=、显示待求点里程设计高程H=)Lbl 2: “K(L)=”? S: “X=”? X：“Y=”？ Y:Prog“SJ-PM”： X I： Y J：Prog SUB2-FS：O+WS: “K=：Locate4,4, S ：“P=：Locate4,4, Z：S K:Prog“SJ-GC” ：“H=”：Locate4,4, H:F-90F: “W=”: Locate4,4,FGoto 4(反算- 输入反算点大概里程K(L)=、输入实测点X=、 输入实测点Y=、显示实测点里程K=、显示实测点偏距P=、显示实测点的设计高程H=)Lbl 3: “K(L)=”? S: “X=”?X：“Y=”？Y:Prog“SJ-PM”： X I： Y J：Prog SUB2-FS：O+WS: “K=：Locate4,4, S ：“P=：Locate4,4, Z：S K:Prog“SJ-GC” ：“H=”：Locate4,4, H:F-90F: “W=”: Locate4,4,FProg“SJ-CQW” Goto4(超欠挖-输入实测点X=、 输入实测点Y=、显示实测点里程K=、显示实测点偏距P=、显示实测点的设计高程H=、输入实测高程H1=、显示超欠挖值及第几圆心上CQY？=)SJ-PM（子程序名-平面线形数据库）（JD1前直线段要素输入）If S 121000（线元起点里程）:Then 1064587.393U（线元起点X坐标）：523392.492V（线元起点Y坐标）：121000O（线元起点里程）：710 0607”G（线元起点方位角）：864.718 H（线元长度）：11045P（线元起点曲率半径）：11045R（线元终点曲率半径）：0 Q（线元左右偏标志：左-1右1直0）：IfEnd（JD1第一缓和曲线要素输入）If S 121864.718（线元起点里程）:Then 1064867.459U（线元起点X坐标）：524210.598V（线元起点Y坐标）：121864.718O（线元起点里程）：7100607”G（线元起点方位角）：90 H（线元长度）：11045P（线元起点曲率半径）：1200R（线元终点曲率半径）：1 Q（线元左右偏标志：左-1右1直0）：IfEnd（JD1圆曲线要素输入）If S 121954.718（线元起点里程）:Then 1064895.540U（线元起点X坐标）：524296.100V（线元起点Y坐标）：121954.718O（线元起点里程）：7301502”G（线元起点方位角）：425.529 H（线元长度）：1200P（线元起点曲率半径）：1200R（线元终点曲率半径）：1 Q（线元左右偏标志：左-1右1直0）：IfEnd（JD1第二缓和曲线要素输入）If S 122380.246（线元起点里程）:Then 1064944.123U（线元起点X坐标）：524716.604V（线元起点Y坐标）：122380.246O（线元起点里程）：9303405”G（线元起点方位角）：90 H（线元长度）：1200P（线元起点曲率半径）：11045R（线元终点曲率半径）：1 Q（线元左右偏标志：左-1右1直0）：IfEnd(两交点夹直线要素)If S 122470.246（线元起点里程）:Then 1064936.279U（线元起点X坐标）：524806.256V（线元起点Y坐标）：122470.246O（线元起点里程）：9504300”G（线元起点方位角）：1571.776 H（线元长度）：11045P（线元起点曲率半径）：11045R（线元终点曲率半径）：0 Q（线元左右偏标志：左-1右1直0）：IfEnd（JD2第一缓和曲线要素输入）If S 124042.022（线元起点里程）:Then 1064779.711U（线元起点X坐标）：526370.215V（线元起点Y坐标）：124042.022O（线元起点里程）：9504300”G（线元起点方位角）：90 H（线元长度）：11045P（线元起点曲率半径）：1200R（线元终点曲率半径）：-1 Q（线元左右偏标志：左-1右1直0）：IfEnd（JD2圆曲线要素输入）If S 124132.022（线元起点里程）:Then 1064771.867U（线元起点X坐标）：526459.867V（线元起点Y坐标）：124132.022O（线元起点里程）：9303405”G（线元起点方位角）：116.56 H（线元长度）：1200P（线元起点曲率半径）：1200R（线元终点曲率半径）：-1 Q（线元左右偏标志：左-1右1直0）：IfEnd（JD2第二缓和曲线要素输入）If S 124248.582（线元起点里程）:Then 1064770.269U（线元起点X坐标）：526576.370V（线元起点Y坐标）：124248.582O（线元起点里程）：8800010”G（线元起点方位角）：90 H（线元长度）：1200P（线元起点曲率半径）：11045R（线元终点曲率半径）：-1 Q（线元左右偏标志：左-1右1直0）：IfEnd（JD2后直线要互输入）If S 124338.582（线元起点里程）:Then 1064775.653U（线元起点X坐标）：526666.203V（线元起点Y坐标）：124338.582O（线元起点里程）：8505115”G（线元起点方位角）：5201.860 H（线元长度）：11045P（线元起点曲率半径）：11045R（线元终点曲率半径）：0 Q（线元左右偏标志：左-1右1直0）：IfEnd(注：有多个平面线形，依照上面的依次变更, 每多一个，就增加一个。每次只需要修改以上的数据或增加一个判断。)SJ-GC（子程序名-竖曲线数据库）(起点为直坡)If S 起点里程And S下一竖曲线起点:Then起点高程 A:起点里程B: 坡度（+或-）C:A+(S-B)CH：IfEndIf S 121000And S(121700-98/2):Then1108.766 A:121000B: 0.0168C:A+(S-B)CH：IfEnd(以下为竖曲线要素)If S本竖曲线起点And S后一竖曲线起点里程段:Then交点高程A: 交点里程B: 交点前坡（+或-）C: 交点后坡（+或-）D:交点半径R: Prog“SUB3-GC”IfEndIf S(121700-98/2) And S(122050-85/2):Then1120.53A:121700B:0.0168C:0.007D:10000R: Prog“SUB3-GC”IfEndIf S(122050-85/2) And S(122550-50/2):Then1122.98A: 122050B: 0.007C:-0.0015D:10000R: Prog“SUB3-GC”IfEndIf S(122550-50/2) And S(123250-45/2):Then1122.23A: 122550B: -0.0015C:-0.0065D:10000R: Prog“SUB3-GC”IfEndIf S(123250-45/2) And S(123650-120/2):Then1117.68A: 123250B: -0.0065C:-0.011D:10000R: Prog“SUB3-GC”IfEnd(注：1.有多个竖曲线，依照上面的依次变更，每多一个，就增加一个。SJ-CQW (子程序名-隧道数据库)“H1=”?FIf F-H5.845:Then (Abs (Z)2+(F-H-0.45)2 )-7.45W: “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圆心)(注：隧道断面有多个圆心，依照上面的依次变更, 每多一个，就增加一个。每次只需要修改以上的数据或增加一个判断。)详细数据见下图:SUB1-ZS(正算子程序)1P C： (P-R)(2HPR) D： 180 E： 0.1739274226 A： 0.3260725774 B： 0.0694318442 K： 0.3300094782 L： 1-L F：1-K M： U+W(Acos(G+QEKW(C+KWD)+Bcos(G+QELW(C+LWD)+Bcos(G+QEFW(C+FWD)+Acos(G+QEMW(C+MWD) 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) YSUB2-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.000001)：then Goto1：Else W+Z W：Goto 0：IfEndLbl 1： 0 Z：Prog SUB1-ZS： (J-Y)sin(F) ZSUB3-GC（高程子程序）C-DQAbs (RQ2)TR(Abs (Q)Q)R:If SB-T:Then 0K:Else If SB+T:Then 0K:DC:Else S-B+TK:IfEnd:IfEndA-(B-S)C-K2(2R)HReturn说明：(正算1秒，反算和超欠挖5秒完成)所有程序名。不同线路，只需改动SJ-PM，SJ-GC，SJ-CQW三个子程序里的内容，其它不变。1为正算，2为反算，3为超欠挖.K= 正算