CASIO5800計算器公路豎曲線全線高程計算程序Word版

1、CASIO5800計算器公路豎曲線全線高程計算程序 “SQXJSCX”Lb1   1CLs:Fix   3“K=”?k     (計算點里程輸入)If   k<67549.755   AND   K66894.3 :Then   -0.00052A : 0B : 67394.3S : 67.37G : 600000R : Goto   2 : IfEndIf   k<68708.391

2、   AND   K67549.755 :Then   -0.0048A : -0.00052B : 68494.3S :66.8G : 100000R : Goto   2 : IfEndIf   k<本段豎曲線終點里程   AND   K前一豎曲線終點里程 :Then   -0.0048前坡（大里程向）A : -0.00052后坡（小里程向）B : 68494.3豎曲線交點里程S :66.8交點高程G : 100000曲率半徑

3、R : Goto   2 : IfEnd依次類推，計算原始數(shù)據(jù)完成輸入，坡度換算成小數(shù)。Lb1   2S-KL ： A-BW : Abs（R*W/2）T : S-TC : S+TD If   KC : Then   G-L*B H : Goto 5 : Else   If   KS : Then   Goto3 : Else   If   KD : Then   Goto4  

4、 : IfEnd : IfEnd : IfEndLb1   3If   W>0 : Then   G+(C-K) 2/2/R-L*BH   : Goto 5 : Else   If W<0 : Then   G-(C-K) 2/2/R-L*BH   : Goto 5 : IfEnd : IfEndLb1   4If   W>0 : Then   G+(K-D) 2/2/R-L*A

5、H   : Goto 5 : Else   If W<0 : Then   G-(K-D) 2/2/R-L*AH   : Goto 5 : IfEnd : IfEnd1 / 181   5H-0.304H     (0.304為中心路肩高差)“H=”：HGoto 1Norm 1:50C:12345J"RESET PW"?I:I=J=>500DimZ   為數(shù)據(jù)庫增加額外變量500個，在SET、SETPFD

6、YS程序根據(jù)實際再增減變量50ZC+22"PASSWORDS"?I:IZC+39      重設要素保護密碼Cls:StopFileName:DATLOCK      要素保護密碼確定認程序Cls:Norm 1:50C"PASSWORDS"?J:Cls:JZC+39=>Stop    2008-08-07修改FileName:SHELL(外殼程序)50C   &

7、#160;         ;在擴充變量預留前50個給別的程序用。如不夠就適量加大。RESET，SHELL，SET，SETPFDYS，這幾個程序中C值必需一致ZC+35I:"STATION-N"?I:IZC+35      設置測站N坐標ZC+36I:"STATION-E"?I:IZC+36      設置測站E坐標ZC+44I:"STATION-

8、Z"?I:IZC+44設置測站Z坐標ZC+45I:"STATION-HI"?I:IZC+45     設置儀高ZC+41I:"GC-DH"?I:IZC+41     放樣點高差常數(shù)Norm 1ZC+23N當前分段要素N坐標ZC+24E當前分段要素E坐標ZC+25M當前分段要素起點樁號ZC+26H當前分段要素起點方位角（正北）單位：弧度ZC+27A當前分段要素起點曲率有左偏負右偏正（注意不半徑）ZC+28R當前分段要素終點曲率有左偏負右偏正（注意不半徑）ZC+29L當

9、前分段要素長度ZC+3WZC+32I:"DAT1 2 3"?I:IZC+32  平曲線要素數(shù)庫類型選擇1為內(nèi)置式，2文件式，3實時輸入(查看當前要素值)0IZC+4I:"1 2 3"?I:IZC+4   選擇本程序模式默認為坐標正算，1坐標反算，2橫向邊仰坡放樣，3，縱向邊仰坡放樣(隧道進出口用到)Lbi 1ZC+4II=1=>Prog"TURNZH"I=2=>Prog"HXBYP"I=3=>Prog"ZXBYP"Lbi ADeg:Nor

10、m 1:ClsZC+1G:"ZH"?G:GZC+1 ;G樁求樁號ZC+2B:"JL"?B:BZC+2 ;B橫向距離，左正右負Lbi BProg "ZBJS"     ;坐標正算Prog "GCJS"     ;高程計算Fix 3:Cls 設置三位小數(shù)"ZH":Locate 4,1,G 第一行顯示樁號"X":Locate 3,2,X  第二行顯示X(N)坐標&

11、quot;Y":Locate 3,3,Y；第三行顯示Y(E)坐標Prog "GCJS"     ;高程計算"Z":Locate 2,4,Z+ZC+41  ；第四行顯示Z坐標Locate 10,4,ZC+2         ；第四行顯示橫向距離0IDo:I+1I:I=1000=>Goto 1:LpWhile Not(Getkey=57 Or Getkey=27) ;鎖定鍵盤，并在幾十秒后自動返回Getk

12、ey=57=> Goto 1Z+ZC+41+ZC+49-ZC+44-ZC+45ZC+48Pol(ZC+35-X,ZC+36-Y)Cls:"DL":Locate 3,1,I:Locate 12,1,ZC+45"DH":Locate 3,2,ZC+48:Locate 12,2,ZC+49"FWJ"J<0=>J+360J:JDMS   2008-08-07 修改JZC+46IZC+47Goto 1FileName：ZBJS  坐標計算程序Prog"READDAT"

13、RadG-MQIF AR=0 :Then If A=R :Then 1J:Else 3J :IfEnd:Else If A=R :Then 2J:Else 3J:IfEnd:IfEndIf J=1 :Then HF:N+QCos(H)X:E+QSin(H)Y:IfEnd   直線段直接計算If J=2 :Then H+QRF:Rec(Abs(2Sin(Abs(0.5QR)÷R),H+QR÷2):N+IX:E+JY:IfEnd  圓弧段直接計算If J=3 :Then Goto 5 :IfEndX+BCos(F+W)XY+BSin(F+

14、W)YDeg:ReturnLbi 5     用五點通用坐標計算計算緩和段0.5(R-A)÷LKAQIKQ2J0.0469100770P:H+IP+JP2U0.2307653449P:H+IP+JP2V0.5P         :H+IP+JP2D0.7692346551P:H+IP+JP2F0.9530899230P:H+IP+JP2T0.1184634425I0.2393143352J0.2844444444ON+Q(ICos(U)+JCos(V)+OCos(D

15、)+JCos(F)+ICos(T)XE+Q(ISin(U)+JSin(V)+OSin(D)+JSin(F)+ISin(T)YH+AQ+KQ2FX+BCos(F+W)XY+BSin(F+W)YDeg:ReturnFileName:TURNZH(坐標反算)Norm 1ZC+37I:"DQD-N"?I:IZC+37     輸入待求點N坐標I=-1=>Prog"INFWJDLDH"ZC+38I:"DQD-E"?I:IZC+38     輸入待求點E坐標0B

16、：M+L÷2G：Prog"ZBJS"Lbi S:RadZC+37-XI:ZC+38-YJIf I=0 And J=0 :Then Goto A:IfEndPol(I,J)Lbi A:Rec(I,J-F):JZC+2G+IG:IF Abs(I)>0.0001 :Then Prog"ZBJS":Goto S    GZC+1FileName：INFWJDLDH實測坐標（用方位角，距離，高差）輸入程序Deg:Norm 1:ClsZC+46J:"DQD-FWJ"?J:JZC+46ZC+4

17、7I:"DQD-DL"?I:IZC+47ZC+48K:"DQD-DZ"?K:KZC+48ZC+49P:"RHT"?P:PZC+49Rec(I,J)ZC+35+IZC+37ZC+36+JZC+38ZC+44+ZC+45+K-PZC+43FileName:HXBYP   橫向邊仰坡放樣程序Cls:Norm 1ZC+40I:"QPD-DL"?I:IZC+40    輸入起坡點與中樁距離常數(shù)ZC+41I:"QPD-DZ"?I:IZC+41

18、60;      輸入起坡點與中樁高差常數(shù)ZC+42I:"i"?I:IZC+42            輸入邊仰坡坡度，左仰坡(路塹)為正，右仰坡(路塹)為負，左邊坡為負，右邊坡為正Cls:ZC+43I:"DQD-Z"?I:IZC+43     輸入實測高程Prog"TURNZH"Prog "GCJS"(ZC

19、+43-Z-ZC+41)ZC+42+ZC+40BFix 3:B-ZC+2 顯示與上一次測點的偏差BZC+2FileName:ZXBYP  縱向仰坡放樣程序(隧道進口使用)Cls:Norm 1:ZC+33I:"QPD-ZH"?I:IZC+33    輸入起坡點樁號ZC+34I:"QPD-Z"?I:IZC+34       輸入起坡點高程ZC+42I:"i"?I:IZC+42    

20、;        仰坡時(隧道進出口仰坡)，進口為正，出口為負)Cls:ZC+43I:"DQD-Z"?I:IZC+43     輸入實測高程Prog"TURNZH"ZC+42(ZC+43-ZC+34)+ZC+33GFix 3:G-ZC+1 顯示與上一次測點的偏差GZC+1FileName:SET 設置程序1Lbi S2Norm 1350C450ZC+2250K6ZC+5I:"1 Or 2"?I:IZC+5 

21、 選擇1線或者2線7180ZC+3÷J:"XZJJ"?J:J÷180ZC+3;斜樁夾角8If I=1:Then ZC+8ZC+6:ZC+9ZC+79ZC+12ZC+16:ZC+13ZC+17:ZC+18ZC+2010Else ZC+10ZC+6:ZC+11ZC+711ZC+14ZC+16:ZC+15ZC+17:ZC+19ZC+2012IfEnd  13ZC+7+3ZC+31140I:"1PQX,2SQX,3YSDS"?I   ;1輸入平曲線要素，2輸入豎曲線要素，3輸入1、2線夾平豎曲

22、線交點(變坡點)數(shù)量以讓程序規(guī)劃內(nèi)存建立要素數(shù)據(jù)為，其它值退出本程序15I0=>Prog"DATLOCK"             2008-07-28 18：23修改16I=1 => Goto 1:I=2 => Goto 2:I=3 => Goto 3  2008-07-28 18：23修改17Prog"SETPFDYS"18Stop19Lbi 1 ;平曲線要素輸入部分20ZC+6D21-1

23、ZD+322For 1J To ZC+1623Cls:"PQX":Locate 9,1,J24ZD+1I:"N"?I:IZD+1      請輸入N坐標25ZD+2I:"E"?I:IZD+2   ;請輸入E坐標26ZD+3I:"R"?I:IZD+3      請輸入半徑R27ZD+4I:"LS1"?I:IZD+4   ;請輸入緩和曲線1長

24、度28ZD+5I:"LS2"?I:IZD+5        請輸入緩和曲線2長度29D+5D30Next311ZD-232Goto S33Lbi 2                             &#

25、160; ;豎曲線要素輸入部分34ZC+7D35-1ZD+336For 1J To ZC+1737Cls:"SQX":Locate 9,1,J38ZD+1I:"BPDZH"?I:IZD+1     ;請輸入變坡點樁號    2008-08-07修改39ZD+2I:"BPDGC"?I:IZD+2   ;請輸入變坡點高程     2008-08-07修改40ZD+3I:"R"?I:

26、IZD+3         ;請輸入半徑R     2008-08-07修改41D+3D42Next431ZD44Goto S45Lbi 3                             ;

27、要素點數(shù)輸入部分46ZC+12I:"1PQXDS"?I:IZC+12  輸入1線平曲線點數(shù),最小值為347ZC+13I:"1SQXDS"?I:IZC+13  輸入1線豎曲線點數(shù),最小值為348ZC+18I:"1JDZH"?I:IZC+18     ;輸入1線平曲線第一個交點的交點樁號49ZC+14I:"2PQXDS"?I:IZC+14  輸入2線平曲線點數(shù),最小值為350ZC+15I:"2SQXDS&quo

28、t;?I:IZC+15  輸入2線豎曲線點數(shù),最小值為351ZC+19I:"2JDZH"?I:IZC+19  輸入2線平曲線第一個交點的交點樁號52C+ZC+22ZC+8               1線路平曲線要素指針偏移基數(shù)53ZC+8+5ZC+12ZC+9           

29、60;  1線路豎曲線要素指針偏移基數(shù)54ZC+9+3ZC+13ZC+10         2線路平曲線要素指針偏移基數(shù)        55ZC+10+5ZC+14ZC+11      2線路豎曲線要素指針偏移基數(shù)56ZC+11+3ZC+15ZC+21     要素數(shù)據(jù)庫結(jié)束指針57ZC+21+1ZC+3058ZC+2

30、1DimZ                 增加額外變量59-1K60Goto SFileName:SETPFDYS   交點要素=>分段要素1Norm 1250C3Rad4ZC+6+5D5ZC+21+1F6F+6DimZ7-1ZF          設置分段要素開始標志8ZC+20V  &#

31、160;      ; 從數(shù)據(jù)庫讀入當前線第一個平曲線交點樁號90G10Pol(ZD+1-ZD-4,ZD+2-ZD-3)  ；計算第一條直線的距離I，計算方位角J11JH 12ZD-4ZF+1          直線段13ZD-3ZF+214V-IZF+315For 1G To ZC+16-2  16Pol(ZD+1-ZD-4,ZD+2-ZD-3)    ；計算第本交

32、點與下一交點的距離I，計算方位角J17JH       2008-08-07 修改18Pol(ZD+6-ZD+1,ZD+7-ZD+2)19IW:J-HK                         ；計算轉(zhuǎn)角K，負值是左轉(zhuǎn)角，正值是右轉(zhuǎn)角   2008-08-07 修改20I

33、F Abs(K)>:Then If K>0:Then K-2K:Else K+2K:IfEnd:IfEnd21ZD+3R  從數(shù)據(jù)庫讀入半徑R22ZD+4L  從數(shù)據(jù)庫讀入緩和曲線1長度23ZD+5M  從數(shù)據(jù)庫讀入緩和曲線2長度24L÷2÷RB  計算緩和曲線1的B0角25M÷2÷RO  計算緩和曲線2的B0角26(Abs(K)-B-O)RP  計算圓曲線長度27L2÷24÷R-L(4) ÷2

34、384÷R(3) Q      內(nèi)移值計算方法228L÷2-L(3) ÷240÷R2S     ;切線增長值計算方法229M2÷24÷R-M(4) ÷2384÷R(3) J      內(nèi)移值計算方法230M÷2-M(3) ÷240÷R2I     ;切線增長值計算方法231Abs(

35、7;2- Abs(K) A32Rtan(Abs(K÷2)+S+J÷Cos(A)-QTan(A)T 計算切線133Rtan(Abs(K÷2)+I+Q÷Cos(A)-JTan(A)U 計算切線234H<0=>H+2H351A36K<0 =>-1A37HZF+4380ZF+5390ZF+640F+6F41F+6DimZ42Rec(T,H)緩和段143ZD+1-IN44ZD+2-JE45NZF+146EZF+247V-TZF+348L=0=>Goto A49HZF+4500ZF+551A÷RZF+652

36、F+6F53F+6DimZ54Lbi A圓弧段550I:0J56L=0=>Goto B57Pol(L-L(3)÷40÷R2,L2÷6÷R-L(4)÷336÷R(4)58Lbi B59Rec(I,H+AJ)60N+IN61E+JE62NZF+163EZF+264V-T+LZF+365H+ABZF+466A÷RZF+567A÷RZF+668F+6F69F+6DimZ70M=0=>Goto C        緩和段271Rec(2

37、RSin(0.5P÷R),H+AB+AP÷R÷2)72N+IN73E+JE74NZF+175EZF+276V-T+L+PZF+377H+AB+AP÷RZF+478A÷RZF+5790ZF+6           此處刪除2008-08-06 修改80F+6F81F+6DimZ82Lbi C          此處插入2008-08-06 修改83R

38、ec(U,H+K)84ZD+1+IZF+1        此處2008-08-07 修改85ZD+2+JZF+2        此處2008-08-07 修改86V-T+L+P+MZF+387V-T+L+P+M-U+WV        此處2008-08-07 修改88D+5D89Next90H+KZF+4910ZF+5920ZF+693F+9DimZ941

39、ZF+795VZF+9   要素終點樁號96DegFileName：READDATZC+32II=3=>Prog"READDAT3"If G<M Or G>M+L  :Then ;驗證要素是否可用，否則重設要素I=1=>Prog"READDAT1"I=2=>Prog"READDAT2"Prog"REALDAT3"  IfEndFileName:READDAT1 平曲線分段要素讀取程序1ZC+30F0JLbi 1 &#

40、160;         搜索要素If G<ZF+3:Then-1=ZF=>Goto E:F-6F:-1J:Goto 1Else 1=ZF+1=>Goto E:-1=J=>Goto 2:1J:F+6F:Goto 1IfEndLbi 2FZC+30ZF+1ZC+23當前分段要素N坐標ZF+2ZC+24當前分段要素E坐標ZF+3ZC+25當前分段要素起點樁號ZF+4C+26當前分段要素起點方位角（正北）ZF+5ZC+27當前分段要素起點曲率有左偏負右偏正（注意不半徑）ZF+6ZC+28 

41、; 當前分段要素終點曲率有左偏負右偏正（注意不半徑）ZF+9-ZF+3ZC+29  當前分段要素長度ReturnLbi ECls"ERROR"StopREADDAT2文件請輸入內(nèi)容為以下內(nèi)容以便測試結(jié)果與本文一致。÷180JGoto 1Lbi 1If G<30008.396:Then Cls:Stop:IfEnd ;判斷待求樁號是否小于本要素起點，如小于就停止運行程序If G<30661.038(終點里程):Then 2851422.090(起點X)N:513896.440(起點Y)E:30008.396M:213

42、6;3252J(起點方位角)H:0(起點曲率)A:0(終點曲率)R:652.642起點至終點長度L:Goto E:IfEndIf G<30921.038:Then 2850878.162N:513535.768E:30661.038M:213°3252JH:0A:1÷1300R:260L:Goto E:IfEndIf G<31426.367:Then 2850666.478N:513385.004E:30921.038M:219°1639JH:1÷1300A:1÷1300R:505.328L:Goto E:IfEndIf G<31686.367:Then 2850346.483N:512998.015E:31426.367M