卡西歐5800公路、鐵路、隧道、橋梁測量計算程序(全)

坐標正算 ZBZS已知一個起算的坐標點，根據(jù)測得的方位角和平距推算另一點的坐標值。程序如下:“X0=”?B: 儀器點X坐標“Y0=”?C: 儀器點Y坐標Lbl0:“FWJ=”?F: 觀測方位角“PJ=”?D: 觀測平距Rec(D,F):B+I→X:C+J→Y:Cls:Fix3:“Xp=”: 待求點X坐標“Yp=”: 待求點Y坐標Locate7,1,X: 顯示在第一行Locate7,2,Y◢ 顯示在第二行Goto0注解：（B，C）是儀器點的坐標，F(xiàn)為方位角，D為平距。在本程序中，采用了循環(huán)（多次）運算多點與儀器點的關系計算法。如果不必要，只算一個，那么把Lbl0和Goto0刪除就行了。坐標反算 ZBFS根據(jù)兩已知坐標點求算其方位角及距離。程序如下:“X0=”?B:“Y0=”?C:“Xp=”?X:“Yp=”?Y:

儀器點X坐標儀器點Y坐標待求點X坐標待求點Y坐標ol(X-B,Y-C):J<0=>J+360→J:I→D:J→F:Cls:Fix3:“FWJ=”: 計算的方位角F?DMS◢

按度分秒顯示“PJ=”: 計算的平距Locate10,3,D◢ 平距顯示在第三行注解:執(zhí)行程序計算時,首先詢問四變量兩個點的坐標（B，C）、（X，Y）。為了明了顯示結果具體是何意思，特地使用字符串“PJ=”表示平距，字符串“FWJ=”表示方位角，但要注意是字符串而非變量名。直線段坐標計算 ZXZBJS單一直線段任意里程位各點的坐標計算。程序如下:直線起點里程→Q:直線起點X坐標→E:直線起點Y坐標→F:直線起點前進方位角→T:Lbl0:2→DimZ:“KM=”?L:“L-1,R+1”?K:“PJ=”?S:

待求點里程待求點位于前進方向左側為 -1,右側為+1待求點到線路設計中心的平距Rec(Abs(Q-L),T):E+I→Z[1]:F+J→Z[2]:IfK=1:ThenRec(S,T+90):ElseRec(S,T+270):IfEnd:Z[1]+I→X:Z[2]+J→Y:Cls:Fix3:“Xp=”: 計算的X坐標“Yp=”: 計算的Y坐標Locate7,1,X: X坐標顯示在第一行Locate7,2,Y◢ Y坐標顯示在第二行Goto0注解:如果本程序不需要反復計算其它里程位的坐標值

,可以將

Lbl0

和Goto0刪除掉。直線段里程偏中距反算 ZXFS程序如下:?Q: 直線起點里程?M: 直線起點X坐標?N: 直線起點Y坐標?T: 直線起點方位角“Xp=”?X:

觀測點

X坐標“Yp=”?Y:

觀測點

Y坐標Pol(X-M,Y-N):J<0=>J+360→J:Rec(I,J-T):Q+I→L:J→S:Cls;Fix3:“KM=”：“DISP=”:Locate7,1,L:Locate10,2,S◢注解:最好將?M:?N:?Q:?T:放入數(shù)據(jù)程序中。如果是單一的直線，也可作用賦值語句。圓曲線段坐標計算 YQXZBJS單一圓曲線段任意里程位各點的坐標計算。程序如下:圓曲線起點里程→Q:圓曲線起點X坐標→E:圓曲線起點Y坐標→F:圓曲線起點前進切線方位角→ T:圓曲線半徑→R:圓曲線轉向方向(左轉取-1,右轉取+1)→P:Lbl0:6→DimZ:“KM=”?L:“L-1,R+1”?K:“PJ=”?S:

待求點里程待求點位于前進方向左側為 -1,右側為+1待求點到線路設計中心的平距Abs(L-Q)→Z[1]:180Z[1]÷R÷π→Z[2]:Pol(Rsin(Z[2]),R(1-cos(Z[2]))):J→Z[3]:Rec(I,T+Z[3]P):E+I→Z[4]:F+J→Z[5]:T+(Z[3]+90Z[1]÷R÷π)P+90→Z[6]:IfK=1:ThenRec(S,Z[6]):ElseRec(S,Z[6]+180):IfEnd:Z[4]+I→X:Z[5]+J→Y:Cls:Fix3:“Xp=”: 計算的X坐標“Yp=”: 計算的Y坐標Locate7,1,X: X坐標顯示在第一行Locate7,2,Y◢ Y坐標顯示在第二行Goto0注解:如果本程序不需要反復計算其它里程位的坐標值

,可以將

Lbl0

和Goto0

刪除掉。圓曲線段里程偏中距反算 YQXFS程序如下

:?V:?W:?M:?N:?Q:?T:?R:?P:“Xp=”?X:“Yp=”?Y:

圓心X坐標圓心Y坐標圓曲線起點X坐標圓曲線起點Y坐標圓曲線起點里程圓曲線起點切線方位角圓曲線半徑圓曲線轉向方向(-1左,+1右)觀測點X坐標觀測點Y坐標Pol(X-V,Y-W):J<0=>J+360→J:(R-I)P→S:Rec(R,J):Pol(V+I-M,W+J-N):J<0=>J+360→J:Q+(J-T)PRπ÷90→LCls;Fix3:“KM=”：“DISP=”:Locate7,1,L:Locate10,2,S◢注解:最好將?V:?W:?M:?N:?Q:?T:?R:?P:放入數(shù)據(jù)程序中。如果是單一的直線，也可作用賦值語句。在(R-I)P±a→S中有個常數(shù)±a。常根據(jù)隧道特點將設計中心線位置移至到隧道中心線位置：路線左轉時左線用-a，右線用+a；路線右轉時左線用+a，右線用-a。應用時注意判斷經(jīng)驗算后才能正常使用。當Q不是ZY（或HY）樁號而是YZ（或YH）時，Q+(J-T)PRπ÷90→L修改為Q-(J-T)PRπ÷90→L；同時P值符號也要反向，T值也要注意方向（應指向園弧段）。緩和曲線段坐標計算 HHQXZBJS第一緩和曲線段任意里程位各點的坐標計算。程序如下:曲線起點里程→Q:曲線起點X坐標→E:曲線起點Y坐標→F:曲線起點前進切線方位角→ T:圓曲線半徑→R:緩和曲線長度→U:曲線轉向方向(左轉取-1,右轉取+1)→P:Lbl0:5→DimZ:“KM=”?L:“L-1,R+1”?K:“PJ=”?S:

待求點里程待求點位于前進方向左側為 -1,右側為+1待求點到線路設計中心的平距Abs(L-Q)→Z[1]:Pol(Z[1]-Z[1]^(5)÷(40R2U2)+Z[1]^(9)÷(3456(RU)^(4)),Z[1]^(3)÷(6RU)-Z[1]^(7)÷(336(RU)^(3))+Z[1]^(11)(42240(RU)^(5))):÷J→Z[2]:Rec(I,T+Z[2]P):E+I→Z[3]:F+J→Z[4]:T+Z[2]P+60PZ[1]2÷(RUπ)+90→Z[5]:IfK=1:ThenRec(S,Z[5]):ElseRec(S,Z[5]+180):IfEnd:Z[3]+I→X:Z[4]+J→Y:Cls:Fix3:“Xp=”: 計算的X坐標“Yp=”: 計算的Y坐標Locate7,1,X: X坐標顯示在第一行Locate7,2,Y◢ Y坐標顯示在第二行Goto0注解:如果本程序不需要反復計算其它里程位的坐標值 ,可以將Lbl0和Goto0刪除掉。這里是按第一緩和曲線編制的 ,在應用于第二緩和曲線時 ,起算點和轉向方向都要改變,同時原路線方向上的左側變成了右側 ,原路線方向上的右側變成了左側,等等,這點要注意。交點法中邊樁坐標正反算程序JDZFS 主程序1→DimZ:Cls:Fix3:“1.ZS=>XY,2.XY=>FS”?N:“1.Z,2.H1,3.Y,4.H2”?I:→Z[1]:Prog“SJ”:1/O→K:(O-R)/(2UOR)→T:180/π→H:IfN=1:ThenGoto1:ElseGoto2:IfEnd:Lbl1:“KM=”?L: 里程“DISP=”?S: 偏移量Abs(L-Q)→G:Prog“ZS”:Cls:“X0”:“Y0”:Locate6,1,X:Locate7,2,Y◢Goto1:Lbl2:“X0=”?X:X→I:“Y0=”?Y:Y→J:Prog“FS”:Q+G→L:Cls:“KM=”:“DISP=”:Locate8,1,L:Locate10,2,S◢Goto2SJ數(shù)據(jù)子程序IfZ[1]=1:Then43263.193→Q:線元直線起點樁號(HZ)3313243.205→E:線元直線起點X坐標485714.386→F:線元直線起點Y坐標186-30-41.37→W:線元直線起點方位角10^(45→O:線元直線起點半徑10^(45→R:線元直線終點半徑3362.610→U:線元直線起點到終點的長度0→P:線元直線轉向方向（直線為0，左為-1，右為1）IfEndIfZ[1]=2:Then46625.803→Q:線元第一緩和曲線起點樁號(ZH)3309902.287→E:線元第一緩和曲線起點X坐標485333.057→F:線元第一緩和曲線起點Y坐標186-30-41.37→W:線元第一緩和曲線起點方位角10^(45→O:線元第一緩和曲線起點半徑10000→R:線元第一緩和曲線終點半徑470→U:線元第一緩和曲線起點到終點的長度-1→P:線元第一緩和曲線轉向方向（直線為0，左為-1，右為1）IfEndIfZ[1]=3:Then47095.803→Q:線元圓曲線起點樁號(HY)3309434.927→E:線元圓曲線起點X坐標485283.419→F:線元圓曲線起點Y坐標185-09-54.14→W:線元圓曲線起點方位角10000→O:線元圓曲線起點半徑10000→R:線元圓曲線終點半徑591.076→U:線元圓曲線起點到終點的長度-1→P:線元圓曲線轉向方向（直線為0，左為-1，右為1）IfEndIfZ[1]=4:Then47686.879→Q:線元第二緩和曲線起點樁號(YH)3308845.022→E:線元第二緩和曲線起點X坐標485247.631→F:線元第二緩和曲線起點Y坐標181-46-42.33→W:線元第二緩和曲線起點方位角10000→O:線元第二緩和曲線起點半徑10^(45→R:線元第二緩和曲線終點半徑470→U:線元第二緩和曲線起點到終點的長度-1→P:線元第二緩和曲線轉向方向（直線為0，左為-1，右為1）IfEndZS正算子程序0.1739274226→A:0.3260725774→B:0.0694318442→C:0.3300094782→D:1-D→V:1-C→M:E+G(Acos(W+PHCG(K+CGT))+Bcos(W+PHDG(K+DGT))+Bcos(W+PHVG(K+VGT))+Acos(W+PHMG(K+MGT))) →X:F+G(Asin(W+PHCG(K+CGT))+Bsin(W+PHDG(K+DGT))+Bsin(W+PHVG(K+VGT))+Asin(W+PHMG(K+MGT))) →Y:W+PHG(K+GT)+90→V:X+Scos(V→X:Y+Ssin(V→YFS反算子程序W-90→Z:Abs((Y-F)cos(Z)-(X-E)sin(Z)) →G:0→S:Lbl0:Prog“ZS”:Z+PHG(K+GT)→D:(J-Y)cos(D)-(I-X)sin(D) →S:IfAbs(S)<10^(-6):ThenGoto1:ElseG+S→G:IfEnd:Goto0:Lbl1:0→S:Prog“ZS”:(J-Y)/sin(V)→S交點法坐標計算 JDFZBJS第一種：單一交點法計算任意里程位的坐標。程序如下:待求段的交點X坐標→V:待求段的交點Y坐標→W:前一段的交點X坐標→M:前一段的交點Y坐標→N:直緩點ZH的里程→Q:緩直點HZ的里程→Z:待求段圓曲線半徑→R:待求段緩和曲線長度→U:待求段的切線長度→T:待求段轉向方向（左轉取 -1，右轉取+1）→P:待求段轉向偏角→A:10→DimZ:Cls:“KM=”?L:“L-1,R+1”?K:“PJ=”?S:

待求點里程待求點位于前進方向左側為 -1,右側為+1待求點到線路設計中心的平距Lbl0:Pol(V-M,W-N):J→Z[1]:Rec(T,Z[1]+180):V+I→Z[2]:W+J→Z[3]:Lbl1:IfL≤Q:ThenRec(Q-L,Z[1]+180):ElseGoto2:IfEnd:Z[2]+I→Z[4]:Z[3]+J→Z[5]:IfK=1:ThenRec(S,Z[1]+90):ElseRec(S,Z[1]+270):IfEnd:Goto6:Lbl2:Abs(Q-L)→Z[6]:IfU≥Z[6]:ThenPol(Z[6]-Z[6]?(5)÷(40R2U2)+Z[6]^(9)÷(3456(RU)^(4)),Z[6]^(3)÷(6RU)-Z[6]^(7)÷(336(RU)^(3))+Z[6]^(11)÷(42240(RU)^(5))):ElseGoto3:IfEnd:J→Z[7]:Rec(I,Z[1]+Z[7]P):Z[2]+I→Z[4]:Z[3]+J→Z[5]:Z[1]+Z[7]P+60PZ[6]2÷(RUπ)+90→Z[10]:IfK=1:ThenRec(S,Z[10]):ElseRec(S,Z[10]+180):IfEnd:Goto6:Lbl3:Abs(Q-L)→Z[6]:IfL≤Z-U:Then Pol(Rsin(90(2Z[6]-U)÷(Rπ))+U-Rsin(90U÷(Rπ))-U^(3)÷2 2(40R),R-Rcos(90(2Z[6]-U)÷(Rπ))+U÷(24R)):IfEnd:J→Z[7]:Rec(I,Z[1]+Z[7]P):Z[2]+I→Z[4]:Z[3]+J→Z[5]:Z[1]+90P(2Z[6]-U)÷(Rπ)+90→Z[10]:IfK=1:ThenRec(S,Z[10]):ElseRec(S,Z[10]+180):IfEnd:Goto6:Lbl4:IfL≥Z-UAndL≤Z:ThenZ[1]+AP→Z[1]:ElseGoto5:IfEnd:Rec(T,Z[1]):V+I→Z[2]:W+J→Z[3]:Z→Q:-P→P:-K→K:Z[1]+180→Z[1]:Goto2:Lbl5:IfL≥Z:ThenZ[1]+AP→Z[1]:IfEnd:Rec(T,Z[1]):V+I→Z[2]:W+J→Z[3]:Rec(L-Z,Z[1]):Z[2]+I→Z[4]:Z[3]+J→Z[5]:IfK=1:ThenRec(S,Z[1]+90):ElseRec(S,Z[1]+270):IfEnd:Goto6:Lbi6:Z[4]+I→Z[8]:Z[5]+J→Z[9]:Cls:Fix3:“Xp=”: 計算的X坐標“Yp=”: 計算的Y坐標Locate7,1,Z[8]: X坐標顯示在第一行Locate7,2,Z[9]◢ Y坐標顯示在第二行Stop注解:在本程序中考慮了ZH點前和HZ點后的直線情況。程序最后用Stop結束一個里程位的計算,如果要計算的點多,可以在程序最前面(即程序開始處)加個語句行LblA:,在程序最后面將Stop改寫成GotoA。但這樣在停止計算時轉出來較麻煩。第二種：單一交點法計算任意里程位的坐標。主程序JD24→DimZ:Cls:交點里程樁號→K:交點X坐標→X:交點Y坐標→Y:第一緩和曲線長度→B:第二緩和曲線長度→C:圓曲線半徑→R:ZH點方位角→M:偏角(帶符號)→O:M+O→N: HZ點方位角Prog“JDA”Cls:“T1=”: 第一切線長度“T2=”: 第二切線長度“L=”: 曲線總長度“LY=”: 圓曲線長度Locate9,1,S:Locate9,2,T:Locate9,3,L:Locate9,4,Q◢Cls:“E=”: 外距“K(ZH)=”: 直緩點里程Locate11,1,E:Locate8,2,Z[1]◢Cls:“K(HY)=”: 緩圓點里程“K(QZ)=”: 曲中點里程“K(YH)=”: 圓緩點里程“K(HZ)=”: 緩直點里程Locate8,1,Z[2]:Locate8,2,Z[3]:Locate8,3,Z[4]:Locate8,4,Z[5]◢Lbl0:“JSLC”?P: 待求點里程“(-L,Z=0,+R)”?D: 邊樁偏移值(左用負,右用正,中用0)Prog“JDB”IfD<0:ThenCls:“X(L)=”: 左邊樁X坐標“Y(L)=”: 左邊樁Y坐標Locate6,1,F:Locate7,2,G◢Goto0:IfEndIfD=0:ThenCls:“X(Z)=”:“Y(Z)=”:

中樁中樁

X坐標Y坐標Locate6,1,F:Locate7,2,G◢“QXFWJ(Z)=”:Z?DMS◢Goto0:IfEndIfD>0:ThenCls:“X(R)=”: 右邊樁X坐標“Y(R)=”: 右邊樁Y坐標Locate6,1,F:Locate7,2,G◢Goto0:IfEnd子程序——JDAIfO<0: 負為曲線左轉Then-1→W:Else1→W:IfEndWO→A:B2/24/R-B^4/2688/R^3→Z[6]:C2/24/R-C^4/2688/R^3→Z[7]:B/2-B^3/240/R2→Z[8]:C/2-C^3/240/R2→Z[9]:Z[8]+((R+Z[7]-(R+Z[6])cos(A))/sin(A)) →S:Z[9]+((R+Z[6]-(R+Z[7])cos(A))/sin(A)) →T:RAπ/180+(B+C)/2→L:RAπ/180-(B+C)/2→Q:(R+(Z[6]+Z[7])/2/cos(A/2)-R→E:K-S→Z[1]:Z[1]+B→Z[2]:Z[2]+Q/2→Z[3]:Z[1]+L-C→Z[4]:Z[4]+C→Z[5]子程序——JDBX-Scos(M)→Z[19]:Y-Ssin(M)→Z[20]:X+Tcos(N)→Z[21]:Y+Tsin(N)→Z[22]:IfP>Z[1]:ThenGoto1:IfEndZ[1]-P→L:X-(S+L)cos(M)+Dcos(Z+90)→F:Y-(S+L)sin(M)+Dsin(Z+90)→G:M→Z:Goto5Lbl1:IfP>Z[2]:ThenGoto2:IfEndP-Z[1]→L:L→Z[12]:B→Z[13]:Prog“JDC”:Z[19]+Z[14]cos(M)-WZ[15]sin(M)+Dcos(Z+90) →F:Z[20]+Z[14]sin(M)+WZ[15]cos(M)+Dsin(Z+90) →G:M+90WL2/(BRπ)→Z:Goto5Lbl2:IfP>Z[4]:ThenGoto3:IfEndP-Z[1]→L:90(2L-B)/R/π→Z[11]:Rsin(Z[11])+Z[8]→Z[14]:R(1-cos(Z[11]))+Z[6]→Z[15]:Z[19]+Z[14]cos(M)-WZ[15]sin(M)+Dcos(Z+90) →F:Z[20]+Z[14]sin(M)+WZ[15]cos(M)+Dsin(Z+90) →G:M+WZ[11]→Z:Goto5Lbl3:IfP>Z[5]:ThenGoto4:IfEndZ[5]-P→L:L→Z[12]:C→Z[13]:Prog“JDC”:Z[21]-Z[14]cos(N)-WZ[15]sin(N)+Dcos(Z+90)→F:Z[22]-Z[14]sin(N)+WZ[15]cos(N)+Dsin(Z+90) →G:N-90WL2/(CRπ)→Z:Goto5Lbl4:P-Z[5]→L:X+(T+L)cos(N)+Dcos(Z+90)→F:Y+(T+L)sin(N)+Dsin(Z+90)→G:N→Z:Goto5Lbl5:360Frac((Z+360)/360)→Z子程序——JDCIfZ[12]=0:Then0→Z[14]:0→Z[15]:Else2Z[12]-Z[12]^5/40/(RZ[13]) +Z[12]^9/3456/(RZ[13])^4→Z[14]:+Z[12]^11/42240/(RZ[13])^5→Z[15]:IfEnd注意：本程序在計算邊樁前，必須先算中樁后才能算邊樁，不能直接先算邊樁，否則結果會出錯誤值！交會法坐標計算 JHZBJS由4800P改編程序如下:?Q:?V:?X:?Y:I=0:Pol(X-Q,Y-V):IfJ<0:ThenJ+360→J:IfEnd：?D:?M:?A:(I+D+M) ÷2→W:“F0=”:J+2Sin-1(√((W-D)(W-I)÷(ID)))→F◢X→G:Y→H:Rec(D,F):Q+I→X:V+J→Y:X→B:Y→C:G→X:H→Y:Pol(X-B,Y-C):“F=”:IfJ<0:ThenJ+360→J:IfEnd◢A-2sin-1(√((W-D)(W-M)÷(MD)))◢注解:在5800P中只對4800P中第三程序進行應用修改,即采用三角形三條邊為已知條件求算儀器點位的坐標值及提取零方向(變量F±180)或終方向方位角(變量J)；在本程序中重新定義了零方向的坐標為（Q，V），終方向的坐標仍是（X，Y）。沒有象4800P那樣多次調(diào)用“O-S”和“S-O”是不希望使用時把握不住那些是要輸入的變量值，那些是想要的結果值。如果之后要知道儀器點的坐標 ,可以提取變量 B和C,程序設計已轉入。零方向與終方向形成的夾角 A意義不大，在這里未利用?？梢詫崪y A值與計算出的兩方位角值之差進行比較，在使用棱鏡桿觀測時完全不必復核 A值，只有在使用人手扶棱鏡時才需注意，防止偏差過大造成大的錯誤。如果不需要完全掌握A的偏差,如心算,則?A:和:A-2sin-1(√((W-D)(W-M)÷（MD))）◢可以不必輸入計算器中。另外新編一個是通過觀測顯示的方位、平距和高差綜合計算測點坐標的程序如下：Cls:Fix3:28→DimZ:“X1=”?I: 第一方向點 X坐標I→Z[1]:“Y1=”?I: 第一方向點 Y坐標I→Z[2]:“H1=”?I: 第一方向點高程I→Z[3]:“X2=”?I: 第二方向點 X坐標I→Z[4]:“Y2=”?I: 第二方向點 Y坐標I→Z[5]:“H2=”?I: 第二方向點高程I→Z[6]:Pol(Z[4]-Z[1],Z[5]-Z[2]):I→Z[21]: 第一方向點到第二方向點的理論距離J→Z[22]: 第一方向點到第二方向點的方位角Z[22]<0=>360+Z[22] →Z[22]:“D1=”?I: 在第一方向點歸零后觀測的平距I→Z[7]:“h1=”?I: 在第一方向點歸零后觀測的高差I→Z[8]:“v1=”?I: 第一方向點的鏡高I→Z[11]:“D2=”?I: 在第二方向點觀測的平距I→Z[9]:“h2=”?I: 在第二方向點觀測的高差I→Z[10]:“v2=”?I: 第二方向點的鏡高I→Z[12]:“a2=”?I: 在第二方向點觀測的水平角度I→Z[13]:-1222→Z[23]:Z[22]+cos((Z[21]+Z[7]-Z[9])/(2Z[7]Z[21]))Z[22]+180-cos-1((Z[21]2+Z[9]2-Z[7]2)/(2Z[9]Z[21])) →Z[24]:Rec(Z[7],Z[23]):Z[1]+I→Z[14]:Z[2]+J→Z[15]:Rec(Z[9],Z[24]):Z[4]+I→Z[16]:Z[5]+J→Z[17]:(Z[14]+Z[16])/2→Z[18]:計算出的儀器點X坐標(Z[15]+Z[17])/2→Z[19]:計算出的儀器點Y坐標(Z[3]+Z[11]-Z[8]+Z[6]+Z[12]-Z[10])/2→Z[20]:儀器點的高程(未減儀高)Pol(Z[1]-Z[18],Z[2]-Z[19]):I→Z[25]:計算出的儀器點到第一點的平距J→Z[26]:計算出的儀器點到第一點的方位角Z[26]<0=>360+Z[26]→Z[26]:Pol(Z[4]-Z[18],Z[5]-Z[19]):I→Z[27]:計算出的儀器點到第二點的平距J→Z[28]: 計算出的儀器點到第二點的方位角Z[28]<0=>360+Z[28] →Z[28]:Cls:“X0=”:“Y0=”:“Z0=”:Locate7,1,Z[18]:Locate8,2,Z[19]:Locate11,3,Z[20]◢Cls:“FWJ1=”:Z[26]?DMS◢“JL1=”:Locate11,4,Z[25]◢“FWJ2=”:Z[28]?DMS◢“JL2=”:Locate11,4,Z[27]◢Stop這里強調(diào)第一方向點要歸零，第二方向點 Z[13]記錄了水平角但未在程序中參與計算比較，可以用 Z[28]-Z[26]再與Z[13]比較看誤差多大。這是在不會使用全站儀交會法的情況下使用。高程計算 GCJS求算任意里程位的設計高程。程序如下:變坡點里程→O:變坡點高程→G:豎曲線半徑→R:路線后縱坡→E:路線前縱坡→F:切線長度→T:4→DimZ:Cls:“KM=”?L: 待求點里程IfL≤O:Then-0.01E→Z[4]:Else-0.01F→Z[4]:IfEnd:Abs(O-L)→Z[1]:G+Z[1]Z[4]→Z[2]:IfL≤O-TOrL>O+T:ThenZ[2]→Z[2]:IfEnd:IfL>O-TAndL≤O:Then(L-O+T)2÷(2R)→Z[3]:IfE≤F:ThenZ[2]+Z[3]→Z[2]:ElseZ[2]-Z[3]→Z[2]:IfEnd:IfEnd:IfL>OAndL≤O+T:Then(O+T-L)2÷(2R)→Z[3]:IfE≤F:ThenZ[2]+Z[3]→Z[2]:ElseZ[2]-Z[3]→Z[2]:IfEnd:IfEnd:Z[2]→HFix3:“SJGC=”: 計算的設計高程Locate10,4,H: 設計高程顯示在第四行注解:本程序未考慮橫坡值。橫坡計算 HPJS對于有橫坡的計算樣式有多種情況程序如下:超高起點里程→Q:超高終止里程→Z:超高長度→U:超高后的橫坡→M:超高前的橫坡→N:曲線轉向方向→P:“KM=”?L:“PJ=”?S:

,以下是雙線高速公路超高的樣式。IfL≤QOrL≥Z:Then-N→I: 未超高段IfEnd:IfL≥Q+UAndL≤Z-U:Then-PSM÷Abs(PS)→I:IfEnd:IfL>QAndL<Q+U:Then(Q+U-L)(M-N) ÷U-M→I:

全超高段第一超高段IfEnd:IfL>Z-UAndL<Z:Then(L-Z+U)(M-N)

÷U-M→I:

第二超高段IfEndFix3:“HP=”:

計算的設計橫坡Locate10,4,I◢ 設計橫坡顯示在第四行注解:本程序只適用于分離式隧道中，它只考慮了左轉左側和右轉右側的橫坡情況，即要注意S要與P取同一正負符號。5800P隧道施工放樣程序第一種放樣程序 FY-01隧道情況：我隊四個中小隧道全部位于直線段內(nèi)，只有一個隧道處于豎曲線變坡點位置。各個隧道圍巖襯砌類型，洞門明洞設計不一，這里只放置了暗洞內(nèi)的內(nèi)容。由于實際施工中里程位置及尺寸大小都經(jīng)常自行變位，所以現(xiàn)場放樣中要隨時注意變動。FY-5800P 主程序＂X0=＂？X：實測坐標X值＂Y0=＂？Y：實測坐標Y值＂Z0=＂？C：實測高程Z值Prog＂ZLS＂:轉入直線里程反算子程序Prog＂SQX＂:轉入豎曲線子程序→B:xx→Q:B是拱部初襯變數(shù),Q是仰拱初襯變數(shù)；使用前根據(jù)放樣或檢查尺寸來定IfL>45090AndLIfL>45125AndLIfL>45135AndLIfL>45330AndLIfL>45345AndLIfL>45510AndLIfL>45540AndLIfL>45550AndL

≤45125:ThenProg＂5b＂:IfEnd:≤45135:ThenProg＂4a＂:IfEnd:≤45330:ThenProg＂3a＂:Stop:IfEnd:≤45345:ThenProg＂4a＂:IfEnd:≤45383:ThenProg＂5b＂:IfEnd:≤45540:ThenProg＂5b＂:IfEnd:≤45550:ThenProg＂4a＂:IfEnd:≤45650:ThenProg＂3a＂:Stop:IfEnd:IfL>45650AndL≤45660:ThenProg＂4a＂:IfEnd:IfL>45660AndL≤45680:ThenProg＂5b＂:IfEnd:IfL>46085AndL≤46100:ThenProg＂5c＂:IfEnd:IfL>46100AndL≤46110:ThenProg＂4c＂:IfEnd:IfL>46110AndL≤46275:ThenProg＂3a＂:Stop:IfEnd:IfL>46275AndL≤46285:ThenProg＂4b＂:IfEnd:IfL>46285AndL≤46310:ThenProg＂5b＂:IfEnd:Prog＂45DMJS＂: 轉入斷面放樣計算子程序Stop。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。ZLS-5800P 子程序(直線段里程反算 )313242.908→M:直線起點X坐標485714.351→N:直線起點Y坐標43263.493→Q:直線起點里程186°30＇41.37＂→T:直線起點方位角Pol(X-M,Y-N):這里的X、Y自動在主程序中提取J<0=>360+J→J:Rec(I,J-T):Q+I→L:這里反映是由小里程向大里程方向計算的結果存放在變量L里J-2.5→S:設計中線位于隧道中心線的左側2.5米處,計算的偏中距存放在變量S里Cls:Fix3:＂KM=＂:里程標示＂DISP=＂:偏中距標示Locate7,1,L:顯示計算的里程結果到第一行Locate10,2,S◢顯示偏中距結果到第二行,負在左,正在右。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。SQX-5800P 子程序(豎曲線設計高程計算 )46150→A:豎曲線變坡點里程68.787→G:豎曲線變坡點高程30000→R:豎曲線半徑13.86→E:豎曲線變坡點小里程方向的縱坡-9.7→F:豎曲線變坡點大里程方向的縱坡,設計資料上是沒有負號的,這里要加上353.385→T豎曲線切線長5→DimZ:IfL≤A:自動由子程序ZLS中提取LThen-0.001E→Z[1]:高鐵是千分比坡度,公路是百分比坡度,這里用的是高鐵Else0.001F→Z[1]:高鐵是千分比坡度,公路是百分比坡度,這里用的是高鐵IfEnd:G+Abs(A-L)Z[1] →Z[2]:IfL≤A-T:ThenZ[2]→Z[4]:IfEnd:IfL>A-TAndL ≤A:Then(L-A+T)2÷(2R)→Z[3]:IfE≤F:ThenZ[2]+Z[3]→Z[4]:ElseZ[2]-Z[3]→Z[4]:IfEnd:IfEnd:IfL>AAndL ≤A+T:Then(A+T-L)2÷(2R)→Z[3]:IfE≤F:ThenZ[2]+Z[3]→Z[4]:ElseZ[2]-Z[3]→Z[4]:IfEnd:IfEnd:IfL>A+T:ThenZ[2]→Z[4]:IfEnd:Z[4]+2.43→Z[5]: 將計算出來的設計高程（無橫坡的）轉換算到隧道圓心高度C-Z[5]→D: C值直接由主程序中提取 ,求算出實測點與圓心的高差關系Fix3:＂h+-=＂: 表示實測點與圓心的高差關系Locate10,3,D◢ 設置在第三行顯示 ,負在下,正在上。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。斷面計算過程中定義的 14個變量含義:進入斷面放樣工作,開挖或支護厚度(變數(shù))斷面凈空半徑與二襯厚度之和圓心到檢修道的垂直高度二襯邊緣邊墻起點高到檢修道的垂直高度檢修道到邊墻腳的垂直高度仰拱圓心到拱部圓心的垂直高度仰拱半徑與襯砌厚度之和這里是設計拱墻噴砼厚度這里是設計仰拱噴砼厚度襯砌仰拱小半徑邊墻基礎設計寬度拱部部分自檢修道到拱部圓心位的夾角襯砌仰拱小半徑所夾角度仰拱開挖或支護厚度(變數(shù))以下表格的內(nèi)容就是建立數(shù)據(jù)子程序的 (3A、4A、4B、4C、5B、5C)變量ⅢaⅣaⅣbⅣcⅤbⅤcB0~0.150~0.250~0.250~0.250~0.280~0.28T0.100.250.250.280.28R7.057.057.107.107.157.20H2.132.132.132.132.132.13E0.500.500.500.500.500.50F1.271.501.501.501.501.50A12.98112.98112.98112.98112.98112.981W17.6117.7117.7617.7617.8117.86V0.150.250.250.250.280.28M2.832.882.882.932.98N0.80.550.550.520.50U18°40′18″18°40′18″18°40′18″18°40′18″18°40′18″P55°00′47″55°00′47″55°00′47″55°00′47″55°00′47″Q0~0.100~0.250~0.250~0.280~0.28。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。45DMJS-5800P15→DimZ:Pol(Abs(S),D):J→Z[1]:I→Z[2]:(R2-(H-E)2)→Z[3]:-1sin((H-E)/R)→Z[4]:Rec(R+V,Z[4]):I→Z[5]:J→Z[6]:Rec(R-M,U):I→Z[15]:表示小圓心到大圓心的水平距離J→Z[7]:表示小圓心到大圓心的垂直高度H+F→Z[8]:I+√(Abs((M+T)2-(Z[8]-Z[7])2))→Z[9]:(Z[5]-Z[9]-N)/(H-Z[6]+F)→Z[10]:Z[10]表示邊墻的斜率Rec(R+B,Z[4]):I→Z[11]:J→Z[12]:Rec(M+T,U+P):I+Z[15]→Z[13]:表示大圓水平邊緣J+Z[7]→Z[14]:表示大圓垂直邊緣IfZ[1]≥-Z[4]： 表示上半部圓部分Then＂R+-=＂： 表示上半部圓部分實測與設計差Locate10，4,Z[2]-R-B◢ 差值顯示到第四行上

,負為欠

,正為超IfEnd:IfD>-H-FAndD≤-Z[6]AndAbs(S)>Z[9]: 表示下半部邊墻部分Then＂BQ+-=＂： 表示邊墻實測與設計差Locate10，4，Abs(S)-Z[11]+(Abs(D)-Z[12])Z[10] ◢差值顯示到第四行上

,負為欠

,正為超IfEnd：IfAbs(S)>Z[9]AndAbs(S)<Z[9]+NAndD<-H:Then＂QJ+-=＂:Locate10，4,-(D+H+F)◢IfEnd:IfAbs(S)>Z[13]AndAbs(S)<Z[9]AndD<-Z[6]:Then

＂YG+-=＂:

表示仰拱實測與設計差Locate10,4,Abs(D)-Z[7]-

√(Abs((M+Q)

2-(Abs(S)-Z[15])

2))◢差值顯示到第四行上

,負為欠,正為超IfEnd:IfAbs(S) ≤Z[13]AndD<-H:Then＂YG+-=＂:

表示仰拱實測與設計差Locate10,4,Abs(D)+A-

2 2√((W+Q)-S)◢差值顯示到第四行上

,負為欠,正為超IfEnd:3a-5800P(它是個特殊情況 ,無小圓心)→B:進入斷面放樣工作,開挖或支護厚度(變數(shù))7.05→R:斷面凈空半徑與二襯厚度之和2.13→H: 圓心到檢修道的垂直高度0.5→E: 二襯邊緣邊墻起點高到檢修道的垂直高度1.27→F: 檢修道到邊墻腳的垂直高度12.981→A: 仰拱圓心到拱部圓心的垂直高度17.61→W: 仰拱半徑與襯砌厚度之和0.15→V: 這里的0.15是設計拱墻噴砼厚度12→DimZ:Pol(Abs(S),D):J→Z[1]:I→Z[2]:(R2-(H-E)2)→Z[3]:sin-1((H-E)/R)→Z[4]:Rec(R+V,Z[4]):I→Z[5]:J→Z[6]:(W2-(H+F+A)2)→Z[7]:Z[5]-Z[7]→Z[8]:H+F-Z[6]→Z[9]:Z[8]/Z[9]→Z[10]:Z[10]Rec(R+B,Z[4]):I→Z[11]:J→Z[12]:IfZ[1]≥-Z[4]ThenR+-=Locate104,Z[2]-R-B◢,,IfEnd:IfD>-H-FAndD≤-Z[6]:ThenBQ+-=Locate104Abs(S)-Z[11]+(Abs(D)-Z[12])Z[10]◢,,IfEndIfD<-HAndAbs(S)<Z[7]:,,BQThenYG+-=:22Locate10,4,A-√(W-S)-D◢,,IfEnd:ReturnaI bⅤc襯砌斷面圖。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。第二種放樣程序 FY-02這個程序主要應用在電纜溝及道床頂面的測量工作。高程抄平程序 CPJS變坡點里程→O:變坡點高程→G:豎曲線半徑→R:路線后縱坡→E:路線前縱坡→F:切線長度→T:“DG=”?A: 對于全站儀是儀器點的高程值 ,對于水準儀是后視點的高程值“YG=”?B: 對于全站儀是儀器高度 ,對于水準儀是后視讀數(shù)“FYG+-=”?P: 放樣點到圓心的垂直高度 (圓心以下為負 ,以上為正)Lbl0:4→DimZ:Cls:“KM=”?L: 待求點里程IfL≤O:Then-0.001E→Z[4]:Else-0.001F→Z[4]:IfEnd:Abs(O-L)→Z[1]:G+Z[1]Z[4]→Z[2]:IfL≤O-TOrL>O+T:ThenZ[2]→Z[2]:IfEnd:IfL>O-TAndL ≤O:Then(L-O+T)2÷(2R)→Z[3]:IfE≤F:ThenZ[2]+Z[3]→Z[2]:ElseZ[2]-Z[3]→Z[2]:IfEnd:IfEnd:IfL>OAndL≤O+T:Then(O+T-L)2÷(2R)→Z[3]:IfE≤F:ThenZ[2]+Z[3]→Z[2]:ElseZ[2]-Z[3]→Z[2]:IfEnd:IfEnd:Z[2]+2.43→H: 2.43表示內(nèi)軌到圓心的垂直高度“READ=”?K: 前視讀數(shù)(全站儀觀測的高差之相反數(shù) ,水準儀的前視的讀數(shù)值 )Fix3:“GCCZ+-=”: 計算出要量取的高程差值Locate10,4,A+B-K-H-P◢ 高程差值顯示在第四行 ,負低于設計,正高于設計IfK=0: 以K為0值作為本站結束標志 ,否則繼續(xù)下一點觀測計算ThenStop:ElseGoto0:IfEnd這個程序主要要求必須先知道里程才能計算。 K輸入0后仍會先算出結果才能真正結束。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。第三種放樣程序 FY-03QXJS-000 主程序Fix3:Deg:Lbl4:“1.SZ=>NE”:“2.NE=>SZ”:“1,2=”?→Q:“LICHENG”?→S:Prog“QXJSSUB0”Lbl0:IfQ=1:ThenGoto1:IfEnd:IfQ=2:ThenGoto2:IfEndLbl1:“OUTJL=”?→Z:“OUTJD=”?→G:Prog“QXJSSUB1”:“XZB=”:N→N◢“YZB=”:E→E◢“QXFWJ=”:F→F:F?DMS◢Goto4Lbl2:“XZB=”?→B:“YZB=”?→C:B→N:C→E:Prog“QXJSSUB2”:“LICHENG=”:S→S◢“OUTJL=”:Z±a→Z◢“CDG=”?→G:“GC=”:G-(78.7156-0.0058(7900-S)+2.27)→D◢?R:“R+-=”:(Z2+D2)-RGoto4QXJSSUB0 白石坪2#出口數(shù)據(jù)子程序Goto1:Lbl1:IfS<3270.353:Then7913.106→U:2120.363→V:2850.353→O:90-48-35.26→A:420→L:10^(45→P:5500→R:Return:IfEndIfS<4976.210Then7901.826→U2540.185→V:3270.353→O:92-59-50.86→A:1705.857→L:5500→P:5500→R:Return:IfEndIfS<5396.210Then7551.980→U4202.804→V:4976.210→O:110-46-05.09→A:420→L:5500→P:10^(45→R:Return:IfEndIf