自编5800隧道超欠挖综合程序(可用)(非常好).doc

任意坐标正反算及隧道超欠挖程序(2010-10-22)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:“X=”? X：“Y=”？ Y:Prog“SJ-PM”： X I： YJ：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(反算-输入实测点X=、 输入实测点Y=、显示实测点里程K=、显示实测点偏距P=、显示实测点的设计高程H=)Lbl 3: “X=”?X：“Y=”？Y:Prog“SJ-PM”： X I： YJ：Prog SUB2-FS：O+WS: “K=：Locate4,4, S ：“P=：Locate4,4, Z：S K:Prog“SJ-GC” ：“H=”：Locate4,4, HProg“SJ-CQW” Goto4(超欠挖-输入实测点X=、 输入实测点Y=、显示实测点里程K=、显示实测点偏距P=、显示实测点的设计高程H=、输入实测高程H1=、显示超欠挖值及第几圆心上CQY？=)SJ-PM（子程序名-平面线形数据库）If S 59227.681（线元起点里程）:Then 2599818.013U（线元起点X坐标）：496887.918V（线元起点Y坐标）：59227.681O（线元起点里程）：20802839”G（线元起点方位角）：100 H（线元长度）：11045P（线元起点曲率半径）：11045R（线元终点曲率半径）：0 Q（线元左右偏标志：左-1右1直0）：IfEndIf S 59327.681（线元起点里程）:Then 2599730.112U（线元起点X坐标）：496840.237V（线元起点Y坐标）：59327.681O（线元起点里程）：20802839”G（线元起点方位角）：90 H（线元长度）：11045P（线元起点曲率半径）：750R（线元终点曲率半径）：-1 Q（线元左右偏标志：左-1右1直0）：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)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.1739274226A： 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)Goto 1:N=2=Goto 2Lbl 1:X1=?M:Y1=?F:S= ?L: W+V=?A:Rec(L,A):M+IC:F+JDClsX2=:Locate 4,1,C:Y2=:Locate 4,2,DGoto 3Lbl 2: X1=?G:Y1=?H:X=2?N:Y2=?EPol(N-G,E-H)If J0:Then