任意坐標(biāo)正反算及隧道超欠挖程序Microsoft Office Word 文檔

1、任意坐標(biāo)正反算及隧道超欠挖程序1-JIN YU(主程序) Lbl 4: “1ZS,2FS,3CQW”?N(選擇計(jì)算模式，1為正算，2為反算，3為超欠挖)N=1=>Goto 1:N=2=>Goto 2:N=3=>Goto 3：N1 Or N2 Or N3=> Goto 4Lbl 1: “K=”? S：“P=”? Z：If S>=起點(diǎn)里程And S<=終點(diǎn)里程：Then 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”

2、：“H=”：Locate4,4,H：“W=”: Locate4,4,F°Else Goto 1：IfEndGoto 4(正算-輸入待求點(diǎn)里程K=、輸入待求點(diǎn)偏距P=、顯示待求點(diǎn)里程X=、顯示待求點(diǎn)里程Y=、顯示待求點(diǎn)里程設(shè)計(jì)高程H=)Lbl 2:“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=

3、”: Locate4,4,F°Goto 4(反算-輸入實(shí)測(cè)點(diǎn)X=、 輸入實(shí)測(cè)點(diǎn)Y=、顯示實(shí)測(cè)點(diǎn)里程K=、顯示實(shí)測(cè)點(diǎn)偏距P=、顯示實(shí)測(cè)點(diǎn)的設(shè)計(jì)高程H=)Lbl 3: “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, HProg“SJ-CQW” Goto4(超欠挖-輸入實(shí)測(cè)點(diǎn)X=、 輸入實(shí)測(cè)點(diǎn)Y=、顯示實(shí)測(cè)點(diǎn)里程K=、顯示實(shí)測(cè)

4、點(diǎn)偏距P=、顯示實(shí)測(cè)點(diǎn)的設(shè)計(jì)高程H=、輸入實(shí)測(cè)高程H1=、顯示超欠挖值及第幾圓心上CQY？=)SJ-PM（子程序名-平面線形數(shù)據(jù)庫(kù)）If S 59227.681（線元起點(diǎn)里程）:Then 2599818.013U（線元起點(diǎn)X坐標(biāo)）：496887.918V（線元起點(diǎn)Y坐標(biāo)）：59227.681O（線元起點(diǎn)里程）：20802839”G（線元起點(diǎn)方位角）：100 H（線元長(zhǎng)度）：1×1045P（線元起點(diǎn)曲率半徑）：1×1045R（線元終點(diǎn)曲率半徑）：0 Q（線元左右偏標(biāo)志：左-1右1直0）：IfEndIf S 59327.681（線元起點(diǎn)里程）:Then 2599730.112U

5、（線元起點(diǎn)X坐標(biāo)）：496840.237V（線元起點(diǎn)Y坐標(biāo)）：59327.681O（線元起點(diǎn)里程）：20802839”G（線元起點(diǎn)方位角）：90 H（線元長(zhǎng)度）：1×1045P（線元起點(diǎn)曲率半徑）：750R（線元終點(diǎn)曲率半徑）：-1 Q（線元左右偏標(biāo)志：左-1右1直0）：IfEnd(注：有多個(gè)平面線形，依照上面的依次變更, 每多一個(gè)，就增加一個(gè)。每次只需要修改以上的數(shù)據(jù)或增加一個(gè)判斷。)SJ-GC（子程序名-豎曲線數(shù)據(jù)庫(kù)）If S變坡終點(diǎn)里程And S起點(diǎn)里程：Then大里程坡度A：小里程坡度B：變坡點(diǎn)里程O(píng)：變坡點(diǎn)高度G：半徑R：Prog“SUB3-GC”:IfEndIf S變坡終

6、點(diǎn)里程And S前一個(gè)變坡終點(diǎn)里程：Then大里程坡度A：小里程坡度B：變坡點(diǎn)里程O(píng)：變坡點(diǎn)高度G：半徑R：Prog“SUB3-GC”:IfEnd(注：1.有多個(gè)豎曲線，依照上面的依次變更，每多一個(gè)，就增加一個(gè)。每次只需要修改以上的數(shù)據(jù)或增加一個(gè)判斷。2.如整條線只有一個(gè)縱坡比喻為2，那么程序應(yīng)為If S終點(diǎn)里程And S起點(diǎn)里程：Then 0.02A：0.02B：起點(diǎn)里程O(píng)：起點(diǎn)高程G：1×1045R：Prog“SUB3-GC”:IfEnd)SJ-CQW (子程序名-隧道數(shù)據(jù)庫(kù))“H1=”?FIf F-H5.845:Then (Abs (Z)2+(F-H-0.45)2 )-7.45

7、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圓心)(注：隧道斷面有多個(gè)圓心，依照上面的依次變更, 每多一個(gè)，就增加一個(gè)。每次只需要修改以上的數(shù)據(jù)或增加一個(gè)判斷。)詳細(xì)數(shù)據(jù)見(jiàn)下圖:SUB1-ZS(正算子程序)1÷P C： (

8、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

9、 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) ZSU

10、B3-GC（高程子程序）Abs(O-S)L:A-BW:Abs(R*W÷2) T:O-TM:O+TPIf SM:ThenG-L*BH:Goto5: Else If S>=P:Then G+L*AH:Goto5:Else If SO And S>M:Then Goto3:Else If S>O And SP:Then Goto4:IfEnd: IfEnd: IfEnd: IfEndLbl3If W0：Then G+(M-S)2÷2÷R-L*BH:Goto5: Else If W0: Then G-(M-S)2÷2÷R-L*BH:G

11、oto5: IfEnd: IfEndLbl4If WO：Then G+( S-P)2÷2÷R+L*AH:Goto5: Else If WO: Then G-( S-P)2÷2÷R+L*AH:Goto5: IfEnd: IfEndLbl5HH: Return說(shuō)明：(正算1秒，反算和超欠挖5秒完成)所有程序名。不同線路，只需改動(dòng)SJ-PM，SJ-GC，SJ-CQW三個(gè)子程序里的內(nèi)容，其它不變。1為正算，2為反算，3為超欠挖.K= 正算時(shí)，輸入所求點(diǎn)里程。 反算時(shí)得出里程結(jié)果 X=，Y= 正算時(shí)得出結(jié)果。 反算時(shí)輸入實(shí)測(cè)坐標(biāo) P=（左負(fù)右正） 正算時(shí)輸入偏距。

12、 反算時(shí)得出偏距 H= 正算時(shí)得出高程。 反算時(shí)得出高程 H1= 超欠挖時(shí)輸入實(shí)測(cè)高程 CQY?= 得出超欠挖值及第幾圓心上casio fx-5800P坐標(biāo)正反算 Lbl 3:"1ZS,2FS"?NN=1=>Goto 1:N=2=>Goto 2Lbl 1:"X1="?M:"Y1="?F:"S= "?L: "W+V="?A:Rec(L,A):M+IC:F+JDCls"X2=":Locate 4,1,C:"Y2=":Locate 4,2,DGoto 3Lbl 2: "X1="?G:"Y1="?H:"X=2"?N:&