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

隧道超欠挖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坐标。

加强施工管理
培训技术人员
对施工技术人员进行培训， 提高其对隧道施工工艺、 技术和安全的认识。
规范施工操作
制定详细的施工操作规程， 规范施工人员的作业行为， 避免因操作不当导致的超 欠挖。
强化现场监控
采用实时监控系统，对隧 道施工过程进行实时监测 和记录，及时发现和纠正 施工中的问题。
定期检测和维护
混凝土衬砌法
总结词
在隧道内表面浇筑一层混凝土衬砌，以增加隧道结构的稳定性和耐久性。
详细描述
对于隧道内表面存在较大超挖或破损的情况，可采用混凝土衬砌法进行治理。浇筑前需 对隧道内表面进行处理，如清洗、修整等，以确保混凝土与原有结构的良好结合。该方 法可有效提高隧道结构的整体性和耐久性，但施工周期较长，且对隧道通风和排水有一
05
隧道超欠挖的预防措施
提高设计精度
01
02
03
精确地质勘察
在隧道设计前进行详细的 地质勘察，了解隧道围岩 的分布、性质和稳定性， 为设计提供准确依据。
优化设计方案
根低 超欠挖的风险。
引入先进技术
采用BIM技术、地质雷达 等先进技术手段，提高设 计精度和可靠性。
施工方法不当
例如，爆破参数选择不当、钻孔精度不足等。
施工管理不到位
例如，施工队伍技术水平不足、施工监测不及时 等。
地质条件复杂
例如，岩层节理裂隙发育、软弱围岩等。
02
隧道超欠挖的检测方法
直接测量法
总结词
直接测量法是一种传统的隧道超欠挖检测方法，通过实地测 量隧道断面的轮廓尺寸，并与设计断面进行比较，计算超欠 挖的数值。
加强施工监控
总结词
施工监控是及时发现和解决隧道超欠挖问题的重要手段。

自编卡西欧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(注：有多个平面线形，依照上面的依次变更, 每多一个，就增加一个。

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

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

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

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

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

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

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

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

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

加强数据安全保障
采用加密技术保障数据安全，确保用户数据 不被泄露和损坏。
谢谢您的聆听
THANKS
增加多隧道管理功能
支持同时管理多个隧道项目，满足大型工程项 目的需求。
强化数据导入导出功能
支持多种格式的数据导入导出，方便与其他软件进行数据交换和共享。
使用体验提升
简化操作流程
优化用户界面，减少操作步骤，降低使用难 度。
提供智能提示和帮助文档
自动提示操作步骤和注意事项，提供详细的 帮助文档和在线支持。
05
隧道超欠挖程序优化建议
算法优化
提高计算精度
采用更精确的数学模型和算法，减少计算误差，提高 超欠挖量的计算准确性。
优化算法效率
改进算法结构，减少计算时间，提高程序的运行速度 。
引入智能算法
结合人工智能和机器学习技术，自动学习和优化算法 参数，提高计算精度和效率。
功能完善
增加数据可视化功能
提供更直观的图表和数据展示，方便用户理解 和分析超欠挖情况。
《隧道超欠挖程序》ppt课件
汇报人：文小库
2024-01-03
CONTENTS
• 隧道超欠挖概述 • 隧道超欠挖程序介绍 • 隧道超欠挖程序应用案例 • 隧道超欠挖程序效果评估 • 隧道超欠挖程序优化建议
01
隧道超欠挖概述
隧道超欠挖定义
隧道超欠挖是指在隧道施工过程中， 实际开挖断面超过了设计断面或者未 达到设计断面的现象。
统计规律分析
通过统计分析法，发现隧道超欠 挖程序处理后的数据分布更加规 律，符合正态分布。
专家评估反馈
专家评估结果显示，隧道超欠挖 程序效果良好，具有较高的实用 价值和应用前景。
评估结论

（换行即为回车符号）坐标转换程序ＳDe g：Fix 3: ＂ZJD=X＂?A:＂ZJD=Y＂?BLbl 0:＂1.ZS＂:＂2.FS＂:＂N=＂?N2874433.193→Z[1] U坐标X570174.274→Z[2] V坐标Y232692.401→Z[3] O里程100.644833→Z[4] G角度470→Z[5] H10^(45) →P曲率5500→R半径1→Q1÷P→C(P-R)÷(2HPR) →D180÷Π→EIF N=1:Then Goto 1:Else Goto 2:IfEndLbl 1:＂S=＂?S:＂Z=＂?ZAbs(S-O) →W:Prog ＂S1＂＂X=＂:X◢＂Y=＂:Y◢Pol(X-A,Y-B):If J<0:Then ＂FWJ=＂:J+360→J DMS◢＂D=＂:I◢Goto 0Lbl 2:(＂R=＂?R:＂D=＂?D＂X=＂A+D*Cos(R) →I◢＂Y=＂B+D*Sin(R) →J◢)＂X=＂?X:＂Y=＂?Y:X→I:Y→J:Prog ＂S2＂＂S=＂:O+W→S◢＂Z=＂：Z◢Prog ＂GC＂:＂LGD=＂:A◢＂QSGC=＂: ?HH-A→BIF B≤-1.54:Then“Error”:IFEndIF B≥-1.54 And B≤2.37:Then 5.41+(1.54+B)*(5.62-5.41)÷3.91→D:Abs(Z)→Z:Z-D→M:IFEnd IF B≥2.37 And B≤4.88:Then A+2.37→E:Abs(Z)+0.49→Z:√(6.11²-(H-E)²)→D:Z-D→M: IFEndIF B≥4.88 And B≤7.72:Then A+2.32→E:Abs(Z)+0.6→Z:√(6.23²-(H-E)²)→D:Z-D→M:IFEndIF B≥7.72:Then A+3.36→E:√(5.03²-Z²)→D:H-E-D→M:“Shang”:IFEnd“-W,+N=”:M ◢Goto 0坐标正算程序S1（计算测点X,Y坐标）0.1739274226→Z[1]0.3260725774→Z[2]0.0694318442→Z[3]0.3300094782→Z[4]0.669991→F:0.930568→MU+W(Z[1]cos(G+QEZ[3]W(C+Z[3]WD))+Z[2]cos(G+QEZ[4]W(C+Z[4]WD))+Z[2]cos(G+QEFW(C+FWD))+ Z[1]cos(G+QEMW(C+MWD)))→XV+W(Z[1]sin(G+QEZ[3]W(C+Z[3]WD))+Z[2]sin(G+QEZ[4]W(C+Z[4]WD))+Z[2]sin(G+QEFW(C+FWD))+ Z[1]sin(G+QEMW(C+MWD)))→YG+QEW(C+WD)+90→FX+Zcos(F)→XY+Zsin(F)→Y坐标反算程序S2(计算里程，偏距)G-90→TAbs((Y-V)cos(T)-(X-U)sin(T))→W:0→ZLbl 3Prog ＂S1＂:T+QEW(C+WD)→L(J-Y)cos(L)-(I-X)sin(L)→ZIf Abs(Z)<10^(-6):Then Goto 4:Else W+Z→W:Goto 3:IfEndLbl 4:0→Z:Prog＂S1＂:(J-Y)÷sin(F)→Z（曲线数据库）If S≥25900 And S≤26615.555:Then25900→O:11587.421→U:1847.983→V:101。

三心圆隧道超欠挖,不足之处望高手指点二、CASIOfx-5800P计算器编程方法：SDCQWLbI 0←┘Fix9:31→DimZ:“R1”?→Z[1]:“R2”?→Z[2]:“DO2”?→Z[3]:“HO1”?→Z[4]:“HO2”?→Z[5]:“B”?→Z[6]:“H”? →Z[7]:“XA”?→Z[8]:“YA”?→Z[9]:“T”? →Z[10] ←┘LbI 1←┘“X”? →Z[11]:“Y”?→Z[12]:“Z”?→Z[13]:( Z[1]- Z[4]+ Z[7]-( Z[1]- Z[1]хcos(Z[6]))) →Z[14] ←┘LbI 2←┘If Z[13]≥Z[14]:Then√((Z[11]-Z[8])2+(Z[12]-Z[9])2)→Z[15]:Z[13]-(Z [7]-Z[4]) →Z[16]:√(Z[15]2-Z[16]2)→Z[17]:Z[17]-Z[1]→Z[18]:“CQ=”: Z[18]▲Goto 1:Else:Goto 3:IfEnd ←┘LbI 3←┘If Z[13]< Z[14]:Then Z[8]+Z[3]хcos(Z[10]-90)→Z[19]:Z[9]+Z[3]хsin(Z [10]-90)→Z[20]: Z[8]+Z[3]хcos(Z[10]+90)→Z[21]:Z[9]+Z[3]хsin(Z[10] +90)→Z[22]:“L”? →Z[23]:If Z[23]= 1:Then √((Z[11]-Z[19])2+(Z[12]-Z[20])2)→Z[24]:Z[13]-(Z[7]+Z[5])→Z[25]:√(Z[24]2+Z[25]2)→Z[26]:Z[26]-Z[2]→Z[27]: “CQ=”:Z[27]▲Else √((Z[11]-Z[21])2+(Z[12]-Z[22])2)→Z[28]:Z[13]-(Z[7]+Z[5])→Z[2 9]:√(Z[28]2+Z[29]2)→Z[30]:Z[30]-Z[2]→Z[31]: “CQ=”:Z[31]▲IfEnd←┘Goto 1←┘三、程序说明本程序按本标段实际情况编制,适用于三心圆隧道,由于第三半径包含范围小及所处水沟范围内,未编入程序内,该程序配合全站仪使用,隧道断面测量时,将全站仪置在该断面里程点中线上,测出断面上任一点三维坐标,数据输入计算R1→半径1R2→半径2DO2→隧道中线至圆心O2的横向距离HO1→隧道设计标高点至圆心O1的竖向距离HO2→隧道设计标高点至圆心O2的竖向距离B→半径R1与隧道中线的夹角H→里程点隧道设计高程XA→里程点隧道中线X坐标YA→里程点隧道中线Y坐标T→里程点切线方位角X→里程点隧道断面上任一点X坐标Y→里程点隧道断面上任一点Y坐标Z→里程点隧道断面上任一点Z坐标四、使用方法1、执行程序SDCQW，按FILE键调出“SDCQW”文件，按EXE键屏幕显示“R1？”，依次输入基本数据：R2：（若测开挖断面半径需加上二衬厚度、预留沉降量及初支厚度）DO2: 隧道中线至圆心O2的横向距离HO1: 隧道设计标高点至圆心O1的竖向距离HO2: 隧道设计标高点至圆心O2的竖向距离B: 半径R1与隧道中线的夹角H: 待测断面里程点隧道设计高程XA: 待测断面里程点隧道中线纵坐标YA: 待测断面里程点隧道中线横坐标T：待测断面里程点线路中线切线方位角2 输入全站仪采集数据X：待测断面上任一点纵坐标Y：待测断面上任一点横坐标Z：待测断面上任一点高程3 计算超欠数据程序运行，依据输入高程Z自动判断该点在哪一半径范围内，在第一半径范围内，直接算出超欠数据“CQ＝”(“—“欠，“+”超)，若测点在第二半径范围内，会出现“L？”（测点在隧道中线左侧？右侧？），在隧道中线左侧直接输入1,在右侧输入其它任何数字，计算超欠挖数据“CQ＝”五、算例R1→9.1R2→5.55半径2DO2→1.818HO1→1.495HO2→1.554B→30°48′36″H→366.002XA→4477933.954YA→504250.721T→101°06′39.18″全站仪采集数据1(第一半径范围内)X→4477935.751Y→504251.073Z→373.372全站仪采集数据2（隧道中线左测）X→4477941.020Y→504252.108Z→368.920全站仪采集数据2（隧道中线右测）X→4477926.800Y→504249.316Z→368.385键操作屏幕显示操作说明[FILE] Prog List 调出SDCQW文件[EXE] R1？输入半径19.1 9.1[EXE] R2？输入半径25.55 5.55[EXE] DO2? 输入隧道中线至圆心O2的横向距离1.818 1.818[EXE] HO12? 输入隧道设计标高点至圆心O1的竖向距离1.495 1.495[EXE] HO2？输入隧道设计标高点至圆心O2的竖向距离1.554 1.554[EXE] B? 输入半径R1与隧道中线的夹角30°48′36″30°48′36″[EXE] H？输入里程点隧道设计高程366.002 366.002[EXE] XA0？输入里程点隧道中线X坐标4477933.954 4477933.954[EXE] YA？输入里程点隧道中线Y坐标504250.721 504250.721[EXE] T？输入里程点切线方位角101°6′39.18″101°6′39.18[EXE] X? 输入测得里程点隧道断面X坐标4477935.751 4477935.751[EXE] Y? 输入测得里程点隧道断面Y坐标504251.073 504251.073[EXE] Z? 输入测得里程点隧道断面Z坐标373.372 373.372[EXE] CQ=-0.047854508 输出超欠数据结果[EXE] X? 输入测得里程点隧道断面X坐标4477941.020 *******.020[EXE] Y? 输入测得里程点隧道断面Y坐标504252.108 504252.108[EXE] Z? 输入测得里程点隧道断面Z坐标368.920 368.920[EXE] L?程序自动判断测点在第二半径范围内，测点在隧道中线左侧？右侧？），在隧道中线左侧直接输入1,在右侧输入其它任何数字1 1[EXE] CQ=0.002970720 输出超欠数据结果[EXE] X? 输入测得里程点隧道断面X坐标4477926.800 4477926.800[EXE] Y? 输入测得里程点隧道断面Y坐标504249.316 504249.316[EXE] Z? 输入测得里程点隧道断面Z坐标368.385 368.385[EXE] L? 判断点位，隧道中线右侧0 0[EXE] CQ=-0.014906344 输出超欠数据结果。

5800正反算、隧道超欠挖程序.doc-正算主程序(BBC)-Lb1 0-?S：?Z：Prog“SJK”：Abs(S-O) → W：Prog "SUB1"："XS="：X◢-"YS="：Y◢F-90→F(需要时可以让他显示，不要时不必输入)：Goto 0-反算主程序(AAFS)-Lb1 0-?S：? X：？Y：“SCH”?→Z[3]：Prog“SJK”（数据库）：X→ I：Y→J：Prog "SUB2"（反算）："S="：O+W→S◢"Z="：Z-2.3→Z[2]◢IfS<5782:Then286.252+(S-5558)*5.6/1000-(S-5558)2/40000→Z[1]:IfEndIfS≧5782:Then286.252-(S-5782)*5.6/1000→Z[1]:IfEnd “SJH”:Z[1]+0→Z[1] ◢IfZ[3]≧Z[1]+0.3:Then:6.86-√((Z[2]2+(Z[3]-Z[1]-2.27)2)→Z[4] :IfEndIfZ[3]∠Z[1]+0.3:Then:6.86-√(((AbsZ[2]-2.85)2+(Z[3]-Z[1]-1.12)2)→Z[4] :IfEndZ[4]→”CQW”◢“GJ”:Z[1]+2.27-Z[3]→Z[9] ◢“TC”:Z[3]-(Z[1]-0.83+Abs(Z[2]*2/100))→Z[5] ◢“ER”:Z[3]-Z[1]-8.73→Z[6] ◢Goto0正算子程序(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)cos(T)-(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：then0→ Z：Prog "SUB1"：(J-Y)÷sinF→ ZElesW+Z → W：Goto0：IfEnd←┘-子程序SJKifS≥5600（线元起点里程）AndS＜5786.5807（线元止点里程）：Then2723.6697→U（线元起点X坐标）：1586.0304→V（线元起点Y坐标）：5600→O（线元起点里程）：191°10°08.3°→G（线元起点方位角）：186.5807→H（线元长度）：10^（45）→P（线元起点曲率半径）：10^（45）→R（线元终点曲率半径）：0→Q（线元左右偏标志：左负右正）：IfEnd←┘-ifS≥5786.5807（线元起点里程）AndS＜6066.5807（线元止点里程）：Then2540.622 →U（线元起点X坐标）：1549.889→V（线元起点Y坐标）：5786.5807→O（线元起点里程）：191°10°08.3°→G（线元起点方位角）：280→H（线元长度）：10^（45）→P（线元起点曲率半径）：5500→R（线元终点曲率半径）：1→Q（线元左右偏标志：左负右正）：IfEnd三、使用说明-1、规定-(1)以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。

• Lbl 3: “X=”?X：“Y=”？Y:Prog“SJ-Βιβλιοθήκη M”： X→ I： Y→ J：Prog
"SUB2-FS"：O+W→S: “K="：Locate4,4, S ：“P="：Locate4,4, Z：S→ K:Prog“SJ-GC” ：“H=”：Locate4,4, H◢ • Prog“SJ-CQW” ◢ • Goto4
• Lbl 2:“X=”? X：“Y=”？ Y:Prog“SJ-PM”： X→ I： Y→ J：Prog
"SUB2-FS"：O+W→S: “K="：Locate4,4, S ：“P="：Locate4,4, Z：S→ K:Prog“SJ-GC” ：“H=”：Locate4,4, H:F-90→F: “W=”: Locate4,4,F°◢ • Goto 4
隧道超欠挖程序
• 1-JIN YU(主程序)
•
Lbl 4: “1→ZS,2→FS,3→CQW”?N(选择计算模式，1为正算，2为 反算，3为超欠挖) N=1=>Goto 1:N=2=>Goto 2:N=3=>Goto 3 • Lbl 1: “K=”? S：“P=”? Z：Prog “SJ-PM”： Abs(S-O) → W： Prog “SUB1-ZS”： “X=”：Locate4,4,X: “Y=”：Locate4,4,Y:F90→F:S→ K:Prog“SJ-GC”：“H=”：Locate4,4,H：“W=”: Locate4,4,F°◢ • Goto 4
• (注：隧道断面有多个圆心，依照上面的依次变更, 每多一个，就
增加一个。每次只需要修改以上的数据或增加一个判断。)详细 数据见下图:

坐标正反算及超欠挖程序主程序（ZBZFS） "1.SZ => XY"："2.XY => SZ"：?N： If N=3:Then Goto 3：IfEnd If N=4 :Then Prog "0-ZZZB": IfEnd "K(L)= "？S：Prog "SJ-PM" 1÷P→C:(P-R) ÷（2HPR）→D:180÷π →E：If N=1:Then Goto 1：Else Goto 2:IfEnd Lbl 1："s"?s:"z"?z：Abs(S-O)→W：Prog "SUB1":"XS=":X→X◢ "YS=":Y→Y◢ "FS=":F-90→F:F►DMS◢ Goto 1←┘ Lbl 2: "X"?X: "Y"?Y:X→I:Y→J：Prog "SUB2"："S=":O+W→S◢ "Z=":Z→Z◢ "MS:M=1，2"?B: If B=1:Then Prog "SDCQW": IfEnd:Goto 2 LbI 3 "X0"?U： "Y0"?V： "S0"?O： "F0"?G： "LS"?H： "R0"?P： "RN"?R： "Q"?Q： Goto 2: IfEnd: Goto 3 使用说明 数据库主要输入整条线路的曲线要素，再输入完成的情况下，选择正反算，按 1 正算，输入 要算的里程判断你要算的里程在那段曲线要输里面， 如 412348， 然后输入计算里程和偏距， 就可以算出改里程的坐标。

隧道的施工放样程序及CAD计算超欠挖量隧道测量的程序及运用：在测量隧道中由于时代的变化、科学的进步，我们运用的计算工具也在不断的变化.在如今我们测量工作中一般运用的是CASIO4500、4800、4850等型号的科学计算器还是一种有编程功能的计算器。

在隧洞测量时测量人员要根据现场的要求来进行编程，边角程序如：边角后方交会BJHFJH说明：1、测边的已知点作为P1(A,B),未测边的已知点作为P2（C,D)。

测边对角为锐角时K=1,测边对角为钝角时2、K=-1。

3、角度P是以测边方向为起始方向,顺时针观测另一个已知点方向的右角。

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

坐标反算ZBFS说明：1、本程序用于计算直角坐标值已知的两点间的边长和坐标方位角.2、起算点和目标点的坐标分别为(A，B)、(C，D）。

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

4、边长值和方位角值分别自动存放在“V”和“W"中。

“W”的单位为：度“°”。

隧洞断面图如上的程序如下：直线断面放样程序(2）ZXFY2说明:1。

本程序用于计算直线段的如图断面样式的隧洞程系放样程系。

2。

坐标A，B，H,是测算出来的坐标数据。

3。

已知的坐标X,Y是从图纸上的起算点坐标。

4。

J是方位角，是隧洞的轴线方向。

5。

M是偏中，V是实际高程，W是实际测量出来的顶拱位置。

后方交会3HFJHCX1、本程序用于利用3个合适的已知点进行方向后方交会法计算测站坐标。

2、观测、计算时将3个已知点按顺时针方向对应排列，已知点的直角坐标分别为（A，B）、（C，D）和（E，F).对应3个已知点的方向值分别为O、P、Q。

3、L3至L9行的作用是当两相邻方向间的夹角出现直角或平角时将导致不能计算时进行自动处理.4、为提高解算精度和防止错误，宜尽可能使测站点与3个已知点组成较理想的图形，如采取测站点靠近3个已知点组成的三角形的中心区域、避免出现“危险园”图形和增加已知点组成多组后交图形比较计算等措施。

隧道超欠挖测量程序SDCX---子程序6.625→R: 5.5→P: 0.2→A: 0.715→U:1.577→V:0.723→I↙If M≥H+5.79: Then √(Z2+(M-H-U)2)-R-A→W: IfEnd↙If M＜H+5.79: Then √((Z-I)2+(M-H-V)2)-P-A→W: IfEndSDCQW---主程序Lb1 0: ?X：？Y：Prog“PMSJ”： X→ I： Y→ J：Prog "XY2"： O+W→K:"Z="：Z◢Prog“SQX”：“H=”：H◢?M:Prog"SDCX"↙“W=”：W◢Goto 0↙M----实测高程PMSJ---子程序Lbl 1 ： IF S<*** （线元终点里程）:Then*** → G （线元起点方位角） :*** → O （线元起点里程） :*** → U （线元起点 X ）:*** → V （线元起点 Y ） :*** → P （线元起点曲率半径）:*** → R （线元终点曲率半径）: *** → H （线元起点至终点长度）： 0 或 1 、 -1 → Q ：Return:IfEnd ↙Lbl 1 ： IF S<*** （线元终点里程）:Then*** → G（线元起点方位角）:*** → O （线元起点里程） :*** → U （线元起点 X ）:*** → V （线元起点 Y ） :*** → P （线元起点曲率半径）:*** → R （线元终点曲率半径）: *** → H （线元起点至终点长度）： 0 或 1 、 -1 → Q ： Return:IfEndXY1———正算子程序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))+Bc os(G+QEFW(C+FWD))+Acos(G+QEMW(C+MWD))) → X：V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bs in(G+QEFW(C+FWD))+Asin(G+QEMW(C+MWD))) → Y：G+QEW(C+WD)+90→ F： X+ZcosF→ X： Y+ZsinF→ YXY2——反算子程序G-90→T： Abs((Y-V)cosT-(X-U)sin（T）) → W： 0→ Z↙Lbl 0：Prog "XY1"： T+QEW(C+WD)→L：(J-Y)cosL-(I-X)sinL → Z：if AbsZ<10-6：then Goto1：Else W+Z → W：Goto 0：IfEnd↙Lbl 1： 0→ Z：Prog "XY1"： (J-Y)÷sinF→ Z↙SQX——子程序竖曲线Lb1 1↙If k<67549.755 AND K≥66894.3 :Then -0 .00052→A : 0→B : 67394.3→S : 67.37→G : 600000→R : Goto 2 : IfEnd↙If k<68708.391 AND K≥67549.755 :Then-0.0048→A : -0.00052→B : 68494.3→S :66.8→G : 100000→R : Goto 2 : IfEnd↙If k<本段竖曲线终点里程 AND K≥前一竖曲线终点里程 :Then -0.0048→前坡（大里程向）A : -0.00052→后坡（小里程向）B : 68494.3→竖曲线交点里程S :66.8→交点高程G : 100000→曲率半径R :Goto 2 : IfEnd↙………依次类推，计算原始数据完成输入，坡度换算成小数。