隧道超欠挖计算原理与程序设计

隧道超欠挖5800计算程序ZBFS(坐标反算通用版《老版本》?N:?E:?A:?R:?C:?T:?G:?M:?FLbl 0:?X:?Y:If G → P:Else -1 → P:IfEndIf P=1:Then F→U:Else P=-1: Then F+A→U:IfEndIf A>0:Then 1→Z:Else -1→Z:IfEndPol(X-N,Y-E:If G>M:Then Goto 3:IfEndG+PT+Rec(I,J-U→S:J→ＤIf S G:Then ＂K=＂:S→K◢＂Q=＂:ZD=Q◢Goto 9Else Goto 4:IfEndLbl 3:If S G:Then ＂K=＂:S→K◢＂Q=＂:ZD→Q◢Goto 9Lbl 4:C2÷24÷Ｒ-C∧(4÷2688÷R∧(3→W:C÷２-C∧(3÷240÷R2→V:90C÷π÷R→B:Pol(P(S-G-V,ZD-R-W:If J B-90:Then Goto 5:IfEndP(S-G+P(S-G∧(5÷40÷R2÷C2→Z:Z∧(3÷(6RC→H:If A>0:Then D-H→W:Else –D-H→W:IfEndIf RC2WC:Then RC÷W→Z:Else ((RC-(R2C2-2RCWZ÷W→Z:Ifend＂K=＂:G+PZ→K◢K-PZ∧(5÷(40R2C2→V:Z∧(3÷(6RC→H:If A<0:Then -1→Z:Else1→Z:IfEnd((S-V2+(ZD-H2→W:If ZD H:Then ＂Q=＂:W→Q◢Else ＂Q=＂:-W→Q◢IfEndGoto 9Lbl 5(90+Jπ÷180-C÷(2R→J:If RJ P(M-G-2C: Then Goto 6:Ifend＂K=＂:G+P(C+JR→K◢＂Q=＂:R-I→Q◢Goto 9Lbl 6:Pol(P(S-G-T,ZD:M-P(T-Rec(I,J-Abs(A→O:If G>M:Then Goto 7:IfEndIf O>M:Then ＂k=＂:O→K◢＂Q=＂J→Q◢Goto 9Else Goto 8Lbl 7:If O M:Then ＂K=＂:O→K◢＂Q=＂:J→Q◢IfendGoto 9Lbl 8:P(M-O+P(M-O∧(5÷(40R2C2→Z:Z∧(3÷(6RC→H:J-H→W:If RC<2WZ:Then RC÷W→Z:Else((RC-((R2C2-2RCWZ÷W→Z:IfEnd＂K=＂:M-PZ→K◢K-PZ∧(5÷(40R2C2→V:Z∧(3÷(6RC→H:((O-V2+(J-H2→W:If J H:Then＂Q=＂:W→Q◢Else ＂Q=＂:-W→Q◢IfEndGoto 9Lbl 9:?Z:?S:＂Q=＂357.957（变坡点高程）+(K-26966（里程）×0.02098（坡度）+S→O◢O-Z-S→B:＂P=＂: ((5.3+Q2+B2-S→P◢Goto 0ZBZS(坐标正算通用版?D:?A:?R:?C:C2÷24÷Ｒ-C∧(4÷2688÷Ｒ∧(3→P:C÷２-C∧(3÷240÷R2→Q:90C÷π÷R→B(R+Ptan(Abs(A÷2+Q→T:(R+P÷cos(A÷２-R→J:(Abs(-A-(2BπR÷180+2C→L: ＂G＂:D-T→G◢＂H＂:G+C→H◢G+L÷2→W: ＂Z＂:G+L-C→Z◢＂M＂:G+L→M◢?N:?Z:?I:?F:If A<0:Then-1→S:Else 1→S:IfEndF+A÷2+90S→U:J+R→V: ＂B＂:N+Vcos(U→B: ＂O＂:E+Vsin(U→ＯLbl 1:?K:If K G:Then G+T-K→L:F+180→V:F→U:Goto 2:IfEndIf K H:Then K-G→L:L-L∧(5÷90÷(RC2→V:30SL2÷（πRC→L:F+180→P:F+L→Q:F+3L→U:Goto4:IfEndIf K Z:Then F+A÷2+180+90S+(K-W×180S÷(Rπ→L:L+90S→U:Goto 5:IfEndIf K M:Then M-K→L:L-L∧(5÷90÷(RC2→V:30SL2÷（πRC→L:F+A→P:F+A+180-L→Q:F+A-3L→U:Goto 4:IfEndIf K>M:Then K-M+T→L:F+A→V:V→U:Goto 2:IfEndLbl 2:＂X=＂:N+Lcos(V→X◢＂Y=＂:E+Lsin(V→Y◢Goto 6Lbl 3:?Q:＂X=＂:X+Qcos(U+90→X◢＂Y=＂:Y+Qsin(U+90→Y◢Goto 1Lbl 4:＂X=＂:N+Tcos(P+Vcos(Q→X◢＂Y=＂:E:+Tsin(P+Vsin(Q→Y◢Goto 6Lbl 5:＂X=＂:B+Rcos(L→X◢＂Y=＂:O+Rsin(L→Y◢Goto 6Lbl 6:I=1Goto 3◢Goto 1使用说明：反算说明：1、N：交点X坐标， E：交点Y坐标。

自编卡西欧fx-5800P隧道超欠挖程序作者：何威该程序由“卡西欧fx-4800P隧道超欠挖程序”改编而成，其运行原理与操作方法与其相同。

1.主程序（WSK-CQW）"XO"?U："YO"?V："SO"?O："FO"?G："LS"?N："RO"?P："RN"?R：?Q←┘1÷P→C：(P-R)÷(2N PR) →D：180÷π→E←┘Lbl2：?X：?Y：?H：X→I：Y→J←┘G-90→T：(Y-V)cos(T)-(X-U)sin(T)→W：Abs(W)→W：0→Z←┘LbI 4：Prog"SUB1"：T+QEW(C+WD)→L：(J-Y)cos(L)-(I-X)sin(L)→Z←┘If Abs(Z)＜10∧(-6）：Then 0→Z ：Prog"SUB1"：(J-Y)÷sin(F)→Z：Else W+Z→W：Goto 4：IfEnd：O+W→S←┘"Z="：Z◢"S="：S◢√( (H-内轨顶标高-圆心距离内轨顶高差)²+(Abs(Z+线中相对遂中距离)+圆心相对遂中距离)²)-R→M←┘"CQW="：M◢Goto 2←┘2.正算子程序（SUB1）0.1739274226→Z[1]：0.3260725774→Z[2]：0.0694318442→K：0.3300094782→L←┘1-L→F:1-K→M←┘U+W(Z[1]cos(G+QEKW(C+KWD))+Z[2]cos(G+QELW(C+LWD))+Z[2]cos(G+QEFW( C+FWD))+Z[1]cos(G+QEMW(C+MWD))→X←┘V+W(Z[1]sin(G+QEKW(C+KWD))+Z[2]sin(G+QELW(C+LWD))+Z[2]sin(G+QEFW(C+ FWD))+ Z[1]sin(G+QEMW(C+MWD))→Y←┘G+QEW(C+WD)+90→F：X+Zcos(F)→X：Y+Zsin(F)→Y一．补充1.输入显示说明X0 ？线元起点的X坐标Y0 ？线元起点的Y坐标S0 ？线元起点里程F0 ？线元起点切线方位角LS ？线元长度R0 ？线元起点曲率半径RN ？线元止点曲率半径Q ？线元左右偏标志(左偏Q=-1，右偏Q=1，直线段Q=0）S ？所求点的里程Z ？所求点距中线的边距(左侧取负，值右侧取正值，在中线上取零)CQW? 所求点的超欠挖2. 若隧道所测部位有两种圆心时，利用圆心所对应的弧的分段高度划分可设置两个选择条件，于是只需将主程序中“√( (H-内轨顶标高-圆心距离内轨顶高差)²+(Abs(Z+线中相对遂中距离)+圆心相对遂中距离)²)-R→M←┘”中改为“IfH>内轨顶标高+圆心O1分段高度：Then√( (H-内轨顶标高-圆心O1距离内轨顶高差)²+(Abs(Z+线中相对遂中距离)+圆心O1相对遂中距离)²)-R1→M ：Else√( (H-内轨顶标高-圆心O2距离内轨顶高差)²+ (Abs(Z+线中相对遂中距离)+圆心O2相对遂中距离)²)-R2→M：If End←┘”3. 若所测隧道断面无需在经常在线元之间转换时，也可以将主程序中U(线元起点X坐标)，V(线元起点Y坐标)，O(线元起点里程)，G(线元起点切线方位角)，N(线元长度)，P(线元起点半径)，R(线元终点半径)，Q(曲线左偏为-1右偏为1直线为0)等曲线要素赋予定值，可减少在工作中的输入量。

隧道测量超欠挖计算程序(5800)曲线：QXLK“ZF”？U:“RP”?O:“ZY”? K:D“X 0”?D:E“YO”?E :?V:?H:?B:?R:?M:?N:“OZ-ZXX”?A:“R1”?P:“R2”？Q←┘Lb1 8←┘?X:?Y←┘Pol((X-D)，(Y-E) )) ←┘“S ”:= -V(I-O)=→S◢ (字母O)“L ”:=K+πO sin-1(sin(90+V(J-U )))÷180→L◢RAbs(N-M)÷2→T←┘B-L→C←┘1→F←┘IF M＞N:Then-1→F:THEND←┘IF L＜B-T:Then Goto 1:ELSE IFL＜B:THEN Goto 2:ELSE IFL＜B+T:THENGoto3:ELSE IF L＞B+T:THEN Goto4:IFEND: :IFEND: :IFEND: :IFEND←┘Lb1 1←┘“G ”:=H-CM→G◢Goto5←┘Lb1 2←┘“G =”:H-C M +F(T-Abs（C)2÷（2R）→G◢Goto5←┘Lb1 3←┘“G= ”:=H-C N +F(T-Abs（C)2÷（2R）→G◢Goto5←┘Lb1 4←┘“G=” :H-CN→G◢Lb1 5←┘?Z←┘IF Z＞G+5.92:THEN Goto6:ELSE Goto7:IFEND←┘Lb1 6←┘“W=” :√(S2+(Z-G-3.1) 2)-P→W◢Goto8←┘Lb1 7←┘“W=” :√((A+Abs（S)2+(Z-G-1.7) 2)-Q→W◢Goto8←┘说明：Z F？方位角RP？半径ZY？圆曲线起点桩号X0？圆心X坐标Y0？圆心Y坐标V？转向H ? 变坡点高程B? 变坡点桩号R? 竖曲线半径M? 前纵坡N? 后纵坡OZ-ZXX? 圆心距中心线距离R 1 ? 拱顶半径R 2 ？内圆半径XO ? 实测坐标YO ? 实测坐标J ?L ? 桩号S ? 平距G ? 设计高程Z ? 实测高程W ?超欠挖(“+”为超挖“-”为欠挖)算例索子沟二号隧道圆曲线段参数设置：ZF? 324 35 53.4RP？500 ZY？ 29117.543 X O？60807.872 Y O？12833.4328 V？1 H ? 1717.4 B? 29750 R? 12000 M? 0.65％N? -0.1028％ OZ-ZXX? 2.2 R1? 4.286 R2?6.893X=60869.64 Y=12332.5154 H=1720.211计算L=29487.827 S=-4.711 W=0.56966直线：ZXLK“ZF”？U:“QSDZH”? K: “X0”?D: “YO”?E :?H:?B:?R:?M:?N:“OZ-ZXX”?A:“R1”?P:“R2”？Q←┘Lb1 8←┘?X:?Y←┘Pol((X-D)，(Y-E) )) ←┘“L ”:K+COS (J-U)×I→L◢“S ”:SIN (J-U) ×I→S◢ (RAbs(N-M)÷2→T←┘B-L→C←┘1→F←┘IF M＞N:Then-1→F:THEND←┘IF L＜B-T:Then Goto 1:ELSE IFL＜B:THEN Goto 2:ELSE IFL＜B+T:THENGoto3:ELSE IF L＞B+T:THEN Goto4:IFEND: :IFEND: :IFEND: :IFEND←┘Lb1 1←┘“G ”:=H-CM→G◢Goto5←┘Lb1 2←┘“G =”:H-C M +F(T-Abs（C)2÷（2R）→G◢Goto5←┘Lb1 3←┘“G= ”:=H-C N +F(T-Abs（C)2÷（2R）→G◢Goto5←┘Lb1 4←┘“G=” :H-CN→G◢Lb1 5←┘?Z←┘IF Z＞G+5.92:THEN Goto6:ELSE Goto7:IFEND←┘Lb1 6←┘“W=” :√(S2+(Z-G-3.1) 2)-P→W◢Goto8←┘Lb1 7←┘“W=” :√((A+Abs（S)2+(Z-G-1.7) 2)-Q→W◢Goto8←┘说明：Z F？直线方位角 QSDZH? 直援点圆心桩号X O？圆心X坐标 Y O？圆心Y坐标H ? 变坡点高程 B? 变坡点桩号R? 竖曲线半径 M? 前纵坡N? 后纵坡 OZ-ZXX? 圆心距中心线距离R 1 ? 拱顶半径 R 2 ？内圆半径XO ? 实测坐标 YO ? 实测坐标J ? L ? 桩号 S ? 平距G ? 设计高程 Z ? 实测高程 W ?超欠挖(“+”为超挖“-”为欠挖)算例索子沟二号隧道直线段参数设置：ZF? 4 12 28.7 QSDZH? 29463.203 X O？ 60844.56 Y O？12334.7807 H ? 1717.4 B? 29750 R? 12000 M? 0.65％N? -0.1028％OZ-ZXX? 2.2 R1? 4.286 R2? 6.893X=61056.762 Y=12348.355 H=1723.69 计算L=29675.829 S=-2.033 G=1716.888. W=-0.063隧道超欠挖计算程序正算主程序(ZS)：Lb1 0：?S：?Z：Prog“PM-SJ”：Abs(S-O)→ W：Prog "SUB1"："XS="：X ◢"YS="：Y◢F-90→F:S→ K:Prog“SQX”：“H=”：H◢Goto0反算主程序(FS)Lb1 0: ?S：?X：？Y：Prog“PM-SJ”：X→ I：Y→ J：Prog "SUB2"："S="：O+W→S◢"Z="：Z◢S→ K:Prog“SQX”：“H=”：H◢Goto0隧道3心圆放样主程序（CQW）Lb1 1：Fix3:7.315→R：6.19→P:“H 1”?F:?Z:F-H→F:Abs(5.72-Z)→ZIf F≥6.319:Then √(Z2+(F-0.715)2 )-R→W:IfEnd:If F≥1.577 AND F＜6.319 Then √((Z-0.723)2+(F-1.577) 2)-P→W:IfEnd:If≤1.577: Then Z-(P+0.723)→W:IfEnd:“W=”: W◢Goto1R----第一个圆圆心P----第二个圆圆心F----实测高程H----路面纵断设计高程Z----由反算主程序反算得到边距（不需修改）程序中右线输入Abs(5.72-Z)→Z，左线输入Abs(5.72+Z)→ZCQW----计算结果（+超，-欠）隧道二衬断面检测主程序（CQJC）Lb1 1：Fix3:6.625→R：5.5→P: “H1”?F:?Z:F-H→F:Abs(5.72-Z)→ZIf F≥5.79:Then √(Z2+(F-0.715) 2)-R→W:IfEnd:If F＜5.79 Then √(Z-0.723)2+(F-1.577)2)-P→W:IfEnd:“W=”: W◢Goto1R----第一个圆圆心P----第二个圆圆心F----实测高程H----路面纵断设计高程Z----由反算主程序反算得到边距（不需修改）程序中右线输入Abs(5.72-Z)→Z，左线输入Abs(5.72+Z)→ZCQW----计算结果（+超，-欠）正算子程序(SUB1)1÷P→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+ZcosF→X：Y+ZsinF→Y反算子程序(SUB2)G-90→T：Abs((Y-V)cosT-(X-U)sin（T）)→W：0→Z：Lbl 0：Prog "SUB1"：T+QEW(C+WD)→L：(J-Y)cosL-(I-X)sinL→Z：ifAbsZ<1E-6：thenGoto1：ElssW+Z →W：Goto0：IfEndLbl 1：0→Z：Prog "SUB1"：(J-Y)÷sinF→Z子程序（平面线形数据库）PM-SJifS≥45978.226（线元起点里程）Then2214.419→U（线元起点X坐标）：4802.542→V（线元起点Y坐标）：45798.226→O（线元起点里程）：280049’54”→G（线元起点方位角）：200→H（线元长度）：1300→P（线元起点曲率半径）：1×1045→R（线元终点曲率半径）：1→Q（线元左右偏标志：左负右正）：IfEndifS≥45998.226（线元起点里程）Then2262.012→U（线元起点X坐标）：4608.341→V（线元起点Y坐标）：45998.226→O（线元起点里程）：285014’20”→G（线元起点方位角）：238.741→H（线元长度）：1×1045→P（线元起点曲率半径）：1×1045→R（线元终点曲率半径）：0→Q（线元左右偏标志：左负右正）：IfEnd子程序（竖曲线计算公式）SQXLbI 0:578.318→Z[1]:46080→B:32000→R:160→T:0.025→I:0.035→J:？K:B-K→C: 1→F:I>J＝＞-1→FIf K<B-T then 0→A: I→P:Goto 1: IfEnd: If K<B then 1→A: I→P: Goto 1: IfEnd: If K<B+T then 1→A: J→P: Goto 1 :IfEnd: If K>B then 0→A: J →P: Goto 1: IfEndLbI 1: Z[1]-CP+AF(T-Abs (C))2÷2÷R→H:“H”:H◢Goto 0Z——变坡点高程 B——变坡点桩号R——半径 T——切线长I——前纵坡度 J——后纵坡K——待求点桩号 H——待求点高程说明：仪器架至测站点上定向后，观测掌子面任意点，测得数据进入反算主程序FS计算得出：对应里程桩号和边距及对应里程路面纵断设计高程。

曲线（含直线）任意里程中边桩坐标正反算（CASIO fx-4800p&fx-4850）J-SQX(竖曲线数据输入)“J-SQX”：{NHUVMQP}：Z[9]=N“SJD”：Z[10]=H“JDZ”：Z[11]=U“I1”：Z[12]=V“I2”：Z[13]=M“R”：Z[14]=Q“QD”：Z[15]=P“ZD”：“TO J-JS MS”J-PQX(平曲线数据输入,自动切换到J-JSMS)A“JD”B“JDX”C“JDY”F“FJ”O“ZJ：Z-，Y+”RE“LS1”K“LS2”：E＜1=＞E=1E-9⊿K＜1=＞K=1E-9⊿Z[1]=EE÷(24R)-E^4÷(2688RRR)：Z[2]=E÷2-EEE÷(240RR)：X=(EE-KK)÷(24R)÷sinAbsO：Z[3]“T1”=(R+Z[1])tan(AbsO÷2)+Z[2]-X▲Z[4]“T2”=(R+KK÷(24R)-K^4÷(2688RRR))tan(AbsO÷2)+K÷2-KKK÷(240RR)+X▲L=AbsOπR÷180+(E+K)÷2▲J=tan-1((R+Z[1])÷(Z[3]-Z[2]))：X“E”=(R+Z[1])÷sinJ-R▲X=A-Z[3]：Y=X+E：E＜1=＞X“ZY”▲≠=＞X“ZH”▲Y“HY”▲⊿Y“QZ”=X+(L-K-E)÷2+E▲Y=X+L-K：X=X+L：K＜1=＞X“YZ”▲≠=＞Y“YH”▲X“HZ”▲⊿Prog“J-JSMS”J-JSMS(放样模式主程序)“1-ZS,2-FS,4-DMFY” Lb1 0：{Z}：Z“MS”≤1=＞Goto 1：≠=＞Z“MS”≤4=＞Goto 2⊿⊿Lb1 1：{PDW}：PD“BZ”W“BJ”：Prog“JP”：X“X=”▲Y“Y=”▲Goto 0⊿Lb1 2：{XYW}：XYW“BJ”： Prog“JF”：P“P=”▲D“BZ=”▲Z=3=＞Prog“JS”：Prog“DMFY”⊿(运行竖曲线高程计算程序)Z=4=＞Z[26] “H”=12.417+(P-75360)*5.2/1000: Prog“DMFY”⊿（运行单面坡比高程计算,语句中12.417为起点桩号DK75+360的高程,可以根据实际情况进行调整；5.2/1000为单面上坡率,上坡输正值，下坡输负值）Goto 0JS(竖曲线计算主程序)P＜Z[14] =＞Prog“J-SQX”⊿P＞Z[15] =＞ Prog“J-SQX”⊿N=Z[9]：U=Z[11]：V=Z[12]：H=Z[13]：G=H Abs(U-V)÷200：P＜N=＞I=U：M=N-G：M＞P=＞M=P⊿≠=＞I=V：M=N+G：M＜P=＞M=P⊿⊿J=(P-M)2÷(2H)：U-V＞0=＞J=-J⊿H=Z[10]+(P-N)×I÷100+JZ[26] “H”=HJP(平曲线正算子程序)FixmLb1 1：J=F：X=B-Z[3]cosF：Y=C-Z[3]sinF：G=EP≤A-Z[3] =＞I=A-Z[3]-P：M=-I：N=0：H=F+W：Goto5≠=＞P≤A-Z[3]+E=＞I=P-A+Z[3]：H=90II÷(REπ)：O＜0=＞H=-H⊿H=H+W+F： Goto 3≠=＞P≤A-Z[3]+L-K=＞I=P-A+Z[3]-E： Goto 4：≠=＞Goto 2⊿⊿Lb1 2：X=B+Z[4]cos（F+O）：Y=C+Z[4]sin（F+O）：J=F+O+180：G=KP≤A-Z[3]+L=＞I= A-Z[3]+L-P：H=90II÷(REπ)：O＞0=＞H=-H⊿H=H+J+W+180： Goto 3≠=＞I=P-A+Z[3]-L：M=-I：N=0：H=J+W+180： Goto 5Lb1 3：M=I-I^5÷(40RRGG)：N=III÷(6GR)-I^7÷(336RRRGGG)： Goto 5Lb1 4：H=(E+2I)×90÷(πR)：M=RsinH+Z[2]：N=R(1-cosH)+Z[1]：O＜0=＞H=-H⊿H=J+W+H： Goto 5Lb1 5：P≤A-Z[3]+L-K=＞O＜0=＞N=-N⊿≠＞O＞0=＞N=-N⊿⊿Goto 6Lb1 6：X=X+Mcos J-Nsin J+Dcos H：Y=Y+Msin J+N cosJ+Dsin HJF(平曲线反算子程序)FixmU=X：V=Y：D=0：J=F-W：P=A+(Y-C)cos J-(X-B)sin J-4=＞Goto 2：≠=＞P=P+I： Goto 1 Lb1 1：Prog“JP”：J=H-180：I=(V-Y)cos J-(U-X)sin J：Abs I＜1E⊿Lb1 2：D=(V-Y)÷sin HDMFY隧道开挖断面(超欠挖情况)计算程序{HIJMNG}：Z[20]= H“C”：Z[21]=I“A”：Z[22]= J“H1”：Z[23]= G“GC”：Z[24]= M“R1”：Z[25]= N“R2”D≤2.3=＞D“ZD”=Abs(D-2.3)▲≠=＞D＞2.3=＞D“YD”= D-2.3▲⊿⊿(C为圆心O1与O2的高差值,A为圆心O1与O2的宽度差值,H1为圆心O1至设计高程的高度,GC为实测高程,R1为上部第1个半径,R2为下部第2个半径)(显示值ZD为至隧道中心线左边的宽度,YD为至隧道中心线右边的宽度,2.3为设计线路距隧道中心线的宽度,可以实际情况进行改变调整)Z[26]“H”▲(显示所求桩号点设计高程值)Z[27]“O1” =Z[26]+ Z[22](第1个圆心高程)Z[28]“O2”=Z[26]+ Z[22]+ Z[20] (第2个圆心高程)Z[29]= Z[24]×Sin 60▲(此数据为上半弧60度范围的宽度值，可根据实际情况调整；数值可不显示) D>Z[29] =＞Goto 1⊿Z[30]=√（（Z[24]）2-D2）Z[31]= √（Abs(（Z[24]）2-(Z[23]-Z[27])2)）Z[32] “HGD”=Z[27]+ Z[30]- Z[23] ▲Z[33] “SKD”=Z[31]-D ▲Prog“J-JSMS”Lb1 1Z[23]<Z[27] =＞Goto 2⊿Z[31]= √(（Z[25]）2-(Z[23]-Z[28])2)Z[33] “SKD”=Z[21]+ Z[31]-D ▲Prog“J-JSMS”Lb1 2Z[31]=Z[25]-（Z[27]-Z[23]）*0.35/(0.78+Z[20]+Z[22])(0.35为下边墙往边墙底内缩的距离，0.78为设计高程至边墙底部的高度)Z[33] “SKD”= Z[31]-D ▲Prog“J-JSMS”一、程序简介1、本套程序共有2个主程序，5个子程序。

隧道隧隧道隧道的施工放样程序及 CAD 计算超欠挖量隧道测量的程序及运用:在测量隧道中由于时代的变化、科学的进步，我们运用的计算工具也在不断的变化。

在如今我们测量 工作中一般运用的是 CASI04500 、4800、4850等型号的科学计算器还是一种有编程功能的计算器。

在隧洞测量时测量人员要根据现场的要求来进行编程，边角程序如:边角后方交会BJHFJH1L 1 ABCD : Lbl5 : {KSP}||L 2 pol(C-A,D-B)||L 3 Q=9O(1-K)+K SIN -1 (S SIN P/V) l 4 T=W+180-P-QE 6 ||L 7 Goto5说明:2、 K=-1 角度P 是以测边方向为起始方向，顺时针观测另一个已知点方向的右角。

注：理想图形要求实测的 S 边相对于已知边 P i P 2越短越好，角P 越接近180 °越好。

L6 Rec (S,T) : X=A+V 丄 Y=B+W 丄1、测边的已知点作为 P i 测边对角为锐角时 K=1 (A,B )，未测边的已知点作为 P 2( C,D )。

，测边对角为钝角时P.(氐怕 f倉F /ZZBFS|L 1 AB : Fixm : {CD}L2 pol(C-A,D-B) JL3 W < 0 W=W+360L4 IntW +0.01lnt(60 FracW )+0.006 Frac(60 FracW)说明:1、本程序用于计算直角坐标值已知的两点间的边长和坐标方位角。

V ”和“ W ”中。

“W ”的单位为：度“。

”。

亠——Pl (毎坐标反算2、起算点和目标点的坐标分别为（ A , B )、( C , D )。

3、起算点改变时应重新调用程序以改变 A 、B 的值。

4、边长值和方位角值分别自动存放在“ P, CD 〉。

隧道程序卡西欧计算器5800隧道正反算、隧道超欠挖计算程序在隧道中的应用，速度有点慢概述：本图采用蓝星2.3版本程序隧道断面程序由三个子程序构成1.1-SDZS 隧道正算，具有三维坐标功能。

2.1-SDFS 隧道反算。

3.1-CQW 隧道超欠挖检测程序。

其余均为子程序，本程序还有许多需要改进的地方，欢迎大家提供宝贵的意见。

有兴趣的可以把三个程序汇总。

本程序为三维坐标全能程序，正在调试阶段，期待和大家一起提高，一起进步。

隧道超欠挖计算程序1.正算主程序1-SDZSLb1 0：“KP”?S（待求桩号）：“BIANZHU”?Z（边距）：Prog “PM-SJ”：Abs(S-O) →W：Prog "SUB1"：X：Y: F-90→F:S→K:Prog“SQX”：“XS=”（X正算坐标）:Locate6,4,X: “YS=”(Y正算坐标):Locate6,4, Y: “BC=”(方位角):Locate6,4, F°:“H=”（高程）:Locate6,4, H◢Goto 02.反算主程序1-SDFSLb1 0: “KP”?S：“CZX”?X：CZY”？Y（测站坐标）：Prog “PM-SJ”：X→I：Y→J：Prog "SUB2"：O+W→S："S="：Locate6,4,S: "Z="：Locate6,4,Z: S→K:Prog“SQX”：F-90→F: “BC=”:Locate6,4, F°: “H=”：H◢3.隧道3心圆检测程序1-CQW程序说明：为了精简，不含数据库，直接把断面的相关数据插入程序内，一些相关的数据，请仔细去推敲它的来历1. 0.729为设计高程到仰拱临界面的高差2. 10.6为设计标高到O3的高差3. 12为O3的半径4. 0.131为设计高程到O2临界点的高差5. 3.498为O2临界宽度6．0.59为设计高程到O2的高差7. 1.8为设计高程到O1的高差8. 1.4为O2半径9. 5.1为O1半径Fix3:Lb1 1：Prog“SDFS”:“HC”?F：“Z”? Z:If F＜H-0.729:Then √(Z2+(10.6+H-F)2 )-12→W:Eles If F≥H-0.729 AND F ≤H+0.131:Then Abs(Z)-3.498→Z:√((Z2+(H+0.59-F)2-1.4→W :Eles √((F-H-1.8)2+Z2)-5.1→W: IfEnd: IfEnd“CQHD”?（衬砌厚度）A:W=”: W-A→W◢Goto14.正算子程序(SUB1)1÷P→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))+B sin(G+QEFW(C+FWD))+Asin(G+QEMW(C+MWD))) →Y：G+QEW(C+WD)+90→F：X+ZcosF→X：Y+ZsinF→Y5.反算子程序(SUB2)G-90→T：Abs((Y-V)cosT-(X-U)sin（T）) →W：0→Z：Lbl 0：Prog "SUB1"：T+QEW(C+WD) →L：(J-Y)cosL-(I-X)sinL →Z：ifAbsZ<1E-6：thenGoto1：ElssW+Z→W：Goto 0：IfEndLbl 1：0→Z：Prog "SUB1"：(J-Y)÷sinF→Z6.子程序（平面线形数据库）PM-SJIf S ≥98509.97（线元起点里程）Then 2963.804→U（线元起点X坐标）：2770.297→V（线元起点Y坐标）：98509.97→O（线元起点里程）：200°13′51.45″→G（线元起点方位角）：35 →H （线元长度）：10^(45)→P（线元起点曲率半径）：164.4→R（线元终点曲率半径）：1 →Q（线元左右偏标志：左负右正）：IfEndIf S ≥98544.97（线元起点里程）Then2931.426→U（线元起点X坐标）：2757.04→V（线元起点Y坐标）：98544.97→O（线元起点里程）：206°19′57.88″→G（线元起点方位角）：27.891 →H（线元长度）：164.4→P（线元起点曲率半径）：164.4→R（线元终点曲率半径）：1 →Q（线元左右偏标志：左负右正）：IfEndIf S ≥98572.86（线元起点里程）Then2907.597→U（线元起点X坐标）：2742.613→V（线元起点Y坐标）：98572.86→O（线元起点里程）：216°3′10.13″→G（线元起点方位角）：35 →H（线元长度）：164.4→P（线元起点曲率半径）：10^(45)→R（线元终点曲率半径）：1 →Q（线元左右偏标志：左负右正）：IfEnd子程序（竖曲线计算公式）SQX101.43→Z[1](变坡点高程):98792→B(变坡点桩号):720→R(半径):43.092→T(切线长):-.0512→J（后纵坡）:0.0676→I（前纵坡）: K:B-K→C : 1→U:I>J＝＞-1→UIf K<B-T then 0→A: I →P:Goto 1: IfEnd: If K<B then 1→A: I→P: Goto 1: IfEnd: If K<B+T then 1→A: J→P: Goto 1 :IfEnd: If K>B then 0→A: J→P: Goto 1: IfEndLbI 1: Z[1]-CP+AU(T-Abs (C))2÷2÷R→H: Return平曲线要素表角度采用十进制小数模式，请采用普通模式版图观看才能正常显示。

5800隧道超欠挖程序，带有垂直方向，水平方向，圆心方向超欠计算1-CQW4→DimZLb1 1Z[4] →S：“K(L)=”? S:Z[1] →A：“X=”? A：A→Z[1]：Z[2] →B ：“Y=”?B：B→Z[2]:Lb1 2：Cls“1→RD,2→,CD,3→ZD,4→YD”? ZZ=1=> Prog “SJ—PMRD”：Z=2=> Prog “SJ—PMCD”：Z=3=> Prog “SJ—PMZD”：Z=4=> Prog “SJ—PMYD”Prog "SUB2-FS"If Abs（N）>0.001:Then N+S→S:Goto2:IfEndN+S→S : S→Z[4]：“K="：Locate4,4, S ：“P="：Locate4,4, W:S→ K: Z→Z[3]Z[3]=1=> Prog “SJ—GCRD”：Z[3]=2=> Prog “SJ—GCCD”：Z[3]=3=> Prog “SJ—GCZD”：Z[3]=4=> Prog “SJ—GCYD”:H→A: “H(CG) =”? A: H- A →H:“H=”：Locate4,4, H◢Cls“1→KW,2→CZ,3→EC”? ZZ=1=> Prog “SJ-CQW-KW”Z=2=> Prog “SJ-CQW-CZ”Z=3=> Prog “SJ-CQW-EC”Prog "1-CQW"Return(超欠挖-输入实测点X=、输入实测点Y=、显示实测点里程K=、显示实测点偏距P=、输入层高H（CG）=、显示实测点的设计高程H=、输入实测高程H（SC）=、显示超欠挖值及第几圆心上CQ→Y？=)SUB1-ZS(正算子程序)-1÷P→ 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：Abs(S-O) →WU+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QEMW(C+MW D))+Acos(G+QEFW(C+F WD)))→ XV+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEMW(C+MW D))+Asin(G+QEFW(C+F WD)))→ YG+QEW(C+WD)→ FReturnSUB2-FS(反算子程序)-Prog "SUB1-ZS"Pol((Z[1]-X),(Z[2]-Y))I cos(J-F) →NI sin (I-F) →WReturnSUB3-GC（高程子程序）C-D→QAbs (RQ÷2)→TR(Abs (Q)÷Q)→R:If S≤B-T:Then 0→K:Else If S≥B+T:Then 0→K:D→C:Else S-B+T→K:IfEnd:IfEndA-(B-S)C-K2÷(2R)→HReturnSJ--PMRD（子程序名-平面线形数据库）SJ--PMCDSJ--PMZDSJ--PMYD（JD1前直线段要素输入）If S ≥0（线元起点里程）:Then 24163.6449→U（线元起点X坐标）：21069.7772→V （线元起点Y坐标）：0→O（线元起点里程）：2°00′58.5″→G（线元起点方位角）：70.2756 →H（线元长度）：10^(45)→P（线元起点曲率半径）：10^(45)→R（线元终点曲率半径）：0 →Q（线元左右偏标志：左-1右1直0）：IfEnd（JD1第一缓和曲线要素输入）If S ≥70.2756（线元起点里程）:Then 24233.8770→U（线元起点X坐标）：21072.2497→V（线元起点Y坐标）：70.2756→O（线元起点里程）：2°00′58.5″→G （线元起点方位角）：60 →H（线元长度）：10^(45)→P（线元起点曲率半径）：500→R （线元终点曲率半径）：-1→Q（线元左右偏标志：左-1右1直0）：IfEnd（JD1圆曲线要素输入）If S ≥130.2756（线元起点里程）:Then 24293.8605→U（线元起点X坐标）：21073.1610→V（线元起点Y坐标）：130.2756→O（线元起点里程）：358°34′42.6″→G （线元起点方位角）：41.525 →H（线元长度）：500→P（线元起点曲率半径）：500→R （线元终点曲率半径）：-1 →Q（线元左右偏标志：左-1右1直0）：IfEnd（JD1第二缓和曲线要素输入）If S ≥171.8005 （线元起点里程）:Then 24335.2821→U（线元起点X坐标）：21070.4093→V（线元起点Y坐标）：171.8005→O（线元起点里程）：353°49′12.4″→G （线元起点方位角）：60 →H（线元长度）：500→P（线元起点曲率半径）：10^(45)→R （线元终点曲率半径）：-1 →Q（线元左右偏标志：左-1右1直0）：IfEnd(两交点夹直线要素)If S ≥231.8005 （线元起点里程）:Then 24394.6180→U（线元起点X坐标）：21061.5714→V（线元起点Y坐标）：231.8005→O（线元起点里程）：350°22′56.5″→G （线元起点方位角）：334.132 →H（线元长度）：500→P（线元起点曲率半径）：10^(45)→R（线元终点曲率半径）：0 →Q（线元左右偏标志：左-1右1直0）：IfEnd（JD2第一缓和曲线要素输入）If S ≥565.9321 （线元起点里程）:Then 24724.0533→U（线元起点X坐标）：21005.7473→V（线元起点Y坐标）：565.9321→O（线元起点里程）：350°22′56.5″→G（线元起点方位角）：45→H（线元长度）：10^(45)→P（线元起点曲率半径）：1000→R （线元终点曲率半径）：1 →Q（线元左右偏标志：左-1右1直0）：IfEnd（JD2圆曲线要素输入）If S ≥610.9321 （线元起点里程）:Then 24768.4749→U（线元起点X坐标）：20998.5621→V（线元起点Y坐标）：610.9321→O（线元起点里程）：351°40′17.5″→G （线元起点方位角）：61.505 →H（线元长度）：1000→P（线元起点曲率半径）：1000→R （线元终点曲率半径）：1 →Q（线元左右偏标志：左-1右1直0）：IfEnd（JD2第二缓和曲线要素输入）If S ≥672.4372 （线元起点里程）:Then 24829.5668→U（线元起点X坐标）：20991.5298→V（线元起点Y坐标）：672.4372→O（线元起点里程）：355°11′43.8″→G （线元起点方位角）：45 →H（线元长度）：1000→P（线元起点曲率半径）：10^(45)→R （线元终点曲率半径）：1 →Q（线元左右偏标志：左-1右1直0）：IfEnd（两交点夹直线要素输入）If S ≥717.4372 （线元起点里程）:Then 24874.4592→U（线元起点X坐标）：20988.4338→V（线元起点Y坐标）：717.4372→O（线元起点里程）：356°29′04.8″→G （线元起点方位角）：62.767 →H（线元长度）：10^(45)→P（线元起点曲率半径）：10^(45)→R（线元终点曲率半径）：0 →Q（线元左右偏标志：左-1右1直0）：IfEnd（JD3第一缓和曲线要素输入）If S ≥780.2044 （线元起点里程）:Then 24937.1083→U（线元起点X坐标）：20984.5852→V（线元起点Y坐标）：780.2044→O（线元起点里程）：356°29′04.8″→G （线元起点方位角）：60 →H（线元长度）：10^(45)→P（线元起点曲率半径）：200→R （线元终点曲率半径）：1 →Q（线元左右偏标志：左-1右1直0）：IfEnd（JD3圆曲线要素输入）If S ≥840.2044 （线元起点里程）:Then 24997.0444→U（线元起点X坐标）：20983.9041→V（线元起点Y坐标）：840.2044→O（线元起点里程）：5°04′44.5″→G （线元起点方位角）：258.124→H（线元长度）：200→P（线元起点曲率半径）：200→R （线元终点曲率半径）：1 →Q（线元左右偏标志：左-1右1直0）：IfEnd（JD3第二缓和曲线要素输入）If S ≥1098.3286（线元起点里程）:Then 25175.6814→U（线元起点X坐标）：21145.0467→V（线元起点Y坐标）：1098.3286→O（线元起点里程）：79°01′34.2″→G （线元起点方位角）：60 →H（线元长度）：200→P（线元起点曲率半径）：10^(45)→R （线元终点曲率半径）：1 →Q（线元左右偏标志：左-1右1直0）：IfEnd（JD3后直线要素输入）If S ≥1158.3286（线元起点里程）:Then 25181.1594→U（线元起点X坐标）：21204.7358→V（线元起点Y坐标）：1158.3286→O（线元起点里程）：87°37′13.9″→G （线元起点方位角）：0 →H（线元长度）：10^(45)→P（线元起点曲率半径）：10^(45)→R （线元终点曲率半径）：0 →Q（线元左右偏标志：左-1右1直0）：IfEnd(注：有多个平面线形，依照上面的依次变更, 每多一个，就增加一个。

隧道超欠挖计算公式及原理
隧道超欠挖是指在隧道开挖过程中，挖掘的截面积超过了设计要求的截面积，这可能会导致隧道结构的不稳定和安全隐患。

为了计算隧道超欠挖，可以使用以下公式：
超欠挖量 = (实际挖掘截面积设计截面积) / 设计截面积。

这个公式可以帮助工程师计算出隧道的超欠挖量，从而评估隧道结构的稳定性和安全性。

隧道超欠挖的原理涉及到土力学和结构力学的知识。

在隧道开挖过程中，地下土体会受到不同方向的应力，当挖掘的截面积超过设计要求时，会导致土体的变形和应力分布的改变。

这可能会对隧道结构产生不利影响，如引起地表沉降、隧道结构变形等。

因此，及时准确地计算隧道的超欠挖量对于确保隧道结构的安全和稳定至关重要。

除了计算超欠挖量外，工程师还需要根据具体情况采取相应的补救措施，如加固隧道结构、调整开挖方法等，以确保隧道的安全施工和使用。

在实际工程中，隧道超欠挖的计算和处理需要综合考
虑地质条件、工程技术和安全要求，以确保隧道工程的顺利进行和安全运营。