5800隧道超欠挖程序

平曲线坐标、反算桩号计算程序prog "PQX"{mn}:M"X0":N"Y0":Defm 7:Z[1]=0: Z[2] "H1"=1.940: Z[3]”RS”=6.430: Z[4]”JK”=0.6:Z[5]”RZ”=6.08Lbl 1:{ABGHVWK}:K"JL":G"JX":H"JY":V"I0":W"J0":A"R":B"LS":FixmP=B^2/24/A:Q=B/2-B^3/240/A^2:T=(A+P)tan(Abs W/2)+Q:C=K-T:D=C+B:F=D+π*A*Abs W/180:E=F-Bw<0=>s=-1:≠=>s=1△Lbl 2:{L}:L:L=0=>GOTO 1△Lbl 6:L≤C=>O=K-L:R=G-OcosV:U=H-OsinV:Z=V:GOTO 3△L≥F=>O=L-F+T:Z=V+W:R=G+OcosZ:U=H+OsinZ:GOTO 3△L<D=>O=L-C:prog "XY":Z=V+SR:R=G+(I-T)cosV-SjsinV:U=H+(I-T)sinV+jscosV:GOTO 3△L>E=>O=F-L:prog "XY":O=V+W:Z=O-SR:R=G+(T-I)cosO-SjsinO:U=H+(T-I)sinO+SJcosO:GOTO 3△O=180(L-D+B/2)/π/A:I=AsinO+Q:J=P+A(1-cosO):Z=V+SO:R=G-(T-I)cosV-SjsinV:U=H-(T-I)sinV+SJcosVLbl 3:Z[1]=0.002=>GOTO 5:≠=>O=0△R"XL"◢U"YL"◢Z"ZL"◢Lbl 4:{O}:O"YC":O=0=>GOTO 2△O=-1=>GOTO 5△X=R-OsinZ:Y=U+OcosZ:pol(X-M,Y-N:J<0=>J=J+360△X"XZ"◢Y"YZ"◢J"A0"◢I"DD"◢GOTO 4Lbl 5:{XY}:X"XF":X=-1=>Z[1]=0:GOTO 4△Y"YF":pol(X-R,Y-U:O=Icos(J-Z):AbsO≤0.001=>O"DL"◢L+O◢O=Isin(J-Z):O"YC"◢GOTO 7:≠=>Z[1]=0.002:L=L+O:GOTO 6△prog "XY"U=AB:I=O-O^5/40/U^2:J=O^3/6/U-O^7/336/U^3:R=90O^2/π/U隧道断面检查程序Lbl 7: {ZJ}::J”XJJ”:Z”H2”:Z[6]”NGD”=780.78+(622640-L)0.003: Z[6]”NGD”◢I= Z[6]+ Z[2]+ Z[3]cos60°:I≤Z=> Z[7]”CQ”=√((Z- Z[6]- Z[2])^2+(O-J/2)^2)- Z[3]:Z[7]”CQ”◢GOTO 5△I>Z=> Z[7]”CQ”=√((Z- Z[6]- Z[2]- Z[4]/2/tan60°)^2+(Abs(O-J/2) - Z[4]/2)^2)- Z[5]: Z[7]”CQ”◢GOTO 5△以上为隧道上导及中导计算程序。

隧道超欠挖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-5800正反算、隧道超欠挖计算程序正算主程序(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◢Goto 0反算主程序(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◢Goto 0隧道3心圆放样主程序（SD-3XY）正算子程序(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+F WD))+Acos(G+QEMW(C+MWD))) → X：V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEFW(C+F WD))+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：Goto 0：IfEndLbl 1：0→ Z：Prog "SUB1"：(J-Y)÷sinF→ Z子程序（平面线形数据库）PM-SJifS ≥？（线元起点里程）AndS＜？（线元止点里程）：Then？→U（线元起点X坐标）：？→V（线元起点Y坐标）：？→O（线元起点里程）：？→G（线元起点方位角）：？→H（线元长度）：？→P（线元起点曲率半径）：？→R （线元终点曲率半径）：？→Q（线元左右偏标志：左负右正）：IfEnd子程序（竖曲线计算公式）SQX。

全线坐标正反算及隧道超欠挖程序(线元法)DFL(主程序)Lbl 5: “1→ZS,2→FS,3→初支CQW，4→二衬CQW”?N(选择计算模式，1为正算，2为反算，3为初支超欠挖，4为二衬超欠挖)N=1=>Goto 1:N=2=>Goto 2:N=3=>Goto 3:N=4=>Goto 4Lbl 1: “K=”?S：“P=”?Z：Prog “SJ-D”： Abs(S-O) → W：Prog “ZS”：“X=”：Locate4,4,X:“Y=”：Locate4,4,Y:F-90→F:S→ K:Prog“SJ-GC”：“H=”：Locate4,4,H：“W=”: Locate4,4,F°◢Goto 5(正算-输入待求点里程K=、输入待求点偏距P=、显示待求点里程X=、显示待求点里程Y=、显示待求点里程设计高程H=)Lbl 2:“X=”?X：“Y=”？Y:Prog“SJ-D”： X→ I： Y→J：Prog "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 5(反算-输入实测点X=、输入实测点Y=、显示实测点里程K=、显示实测点偏距P=、显示实测点的设计高程H=)Lbl 3:“X=”?X：“Y=”？Y:Prog“SJ-D”： X→ I： Y→J：Prog "FS"：O+W→S: “K="：Locate4,4, S：“P="：Locate4,4, Z：S→ K:Prog“SJ-GC”：“H=”：Locate4,4, H◢If K≥19910（XS-VA支护类型起点桩号） And K＜19950（XS-VA支护类型终点桩号）:Then Prog “XS-VA”◢IfEnd:If K≥19950（XS-VC支护类型起点桩号） And K＜20000（XS-VC支护类型终点桩号）:Then Prog “XS-VC”◢IfEnd:(隧道断面匹配、依照上面的依次变更, 每多一个支护类型，就增加一个。

正算主程序(ZS) 检验正确的朋友请加我QQ：624787716FIX4：?S：?Z：Prog “PM-SJ”：Abs(S-O)→W：Prog "SUB1"："XS="：X◢"YS="：Y◢F-90→F(需要时可以让他显示)：Pro g“SQX-SJ”:Prog“SQX”：“H=”：H◢反算主程序(FS)FIX4: ?S：?X：？Y：Prog“PM-SJ”：X→I：Y→J：Prog "SUB2"："S="：O+W→S ◢"Z="：Z◢Prog“SQX-SJ”：Prog“SQX”：“H=”：H◢边坡放样主程序（BP-FY）Lb1 0：Prog“ZS”：“H-BG”（中桩与坡脚起算点高差值，比中桩高正，反之负）?A：H +A→B：?P（实测点高程）：?L（坡脚起算点到中桩的距离）：0.75（挖方时一级坡度）→C：1（挖方时二级坡度）→D：1.5（挖方时三级坡度）→E：8（挖方时一级坡高）→G：10（挖方时二级坡高）→M：15（挖方时三级坡高）→N：2（平台宽度）→K：1（填方时一级坡度）→I：1.5（填方时二级坡度）→J：2（填方时三级坡度）→O：2（填方时一级坡高）→Q：8（填方时二级坡高）→R：10（填方时三级坡高）→T：ifP＞B：thenGoto1：Else Goto2Lb1 1：ifP＞B：thenL+C（P-B）→U：P-B→F：IfEnd←┘ifP＞B+G：thenL+CG+K+D（P-B-G）→U：P-B-G→F：IfEnd←┘ifP＞B+G+M：thenL+GC+2k+MD+E（P-B-G-M）→U：P-B–G-M→F：IfEn d：Goto3←┘Lb1 2：ifP≤B：thenL+I（B -P）→U： B -P→F：I fEnd←┘ifP≤B-Q：thenL+IQ+K+J（B -Q-P）→U：B-Q-P→F：IfEnd←┘ifP≤B-Q-R：thenL+IQ+2k+JR+O（B-Q-R-P）→U：B-Q-R-P→F：IfEn d：Goto3←┘Lb1 3：U-AbsZ→V：ifZ＜0：thenZ-V→Z：Else z+v→z：IfEnd←┘“Z=”：Z◢计算得出正确的宽度，路线左为负，右为正。

1-1－SUIDAO－ZCX (隧道主程序)2-Lbl 5:" ZS=1,FS=2, SD=3, QH=4 "?N3-N=1=>Goto 1:N=2=>Goto 2:N=3=>Goto 3:N=4=>Goto 4:4-Lbl 1: “(JSZH)K=”?S（计算桩号）：“（＋Z,－Y）Z=”?Z（计算偏距，左正右负）：Prog “SJ－PM”（平曲线数据库）：Abs(S－O) → W：Prog "SUB 1－ZS"（正算子程序）：“X="：Locate4,4,X:"Y="：Locate4,4,Y：F－90→F （该处F－90为计算正交90°边桩，要算斜交边桩即可改为？N：F－N→F）：S→ K：Prog“SJ－GC”（竖曲.线数据库）：“H=”：Locate4,4,H◢5-Goto 56-Lbl 2：“X=”?X：“Y=”？Y：Prog“SJ－PM”（平曲线数据库）：X→ I：Y→ J：Prog "SUB2－FS"（反算子程序）：O+W→S ：“K="：Locate4,4, S ：“（＋Z,－Y）P="：Locate4,4, Z：S→ K：Pro g“SJ－GC”（竖曲线数据库）：“H=”：Locate4,4, H◢7-Goto 58-Lbl 3 ：“X=”?X：“Y=”？Y：Prog“SJ－PM”（平曲线数据库）：X→ I：Y→ J：Prog "SUB2－FS"（反算子程序）：O+W→S：“K="：Locate4,4, S ：“（＋Z,－Y）P="：Locate4,4, Z：S→ K：Prog“SJ－GC” ：“H=”：Locate4, 4, H◢9-Prog“SJ－CQW”（隧道超欠挖数据库）: “（SJ-GD）H="：Locate4,4,F:“R:(+ CW,-QW)=”：Locate4,4,W◢10-Goto511-Lbl 4 ：“K=”?S（计算桩号）：“（＋Z,－Y）P=”?Z（计算偏距，左正右负）：Prog“SJ－PM”（平曲线数据库）：Abs(S－O) → W：Prog "SUB1－ZS"（正算子程序）：F－90→F（该处F－90为计算正交90°边桩，要算斜交边桩即可改为？N：F－N→F）：“JJ=”?J （涵洞轴线于路线主线右夹角）：0→A：0→B：Prog "QHZBJS"（桥涵坐标计算程序）：12-Goto5以上为第一个主程序13-SJ-PM（子程序名-平曲线形数据库）14-If S ≥38990.919（线元起点里程）:Then 9502.991→U（线元起点X坐标）：9616.568→V（线元起点Y坐标）：38990.919→O（线元起点里程）：298°11’04.8”→G（线元起点方位角）：130 →H（线元长度）：1×1045→P（线元起点曲率半径）：850→R（线元终点曲率半径）：－1 →Q（线元左右偏标志：左-1右1）：IfEnd15-If S ≥39120.919:Then9561.437→U：9500.485→V：39120.919→O：293°48’11.62”→G：295.75 →H：850→P：850→R：－1 →Q：IfEnd16-If S ≥39416.670:Then9631.804→U：9214.762→V：39416.670→O：273°52’03.4”→G：130→H：850→P：1×1045→R：－1 →Q：IfEnd17-If S ≥39546.67:Then9633.951→U：9084.813→V：39546.67→O：269°29’10.2”→G：486.016 →H：1×1045→P：1×1045→R：0→Q：IfEnd18-If S ≥40032.687:Then9629.593→U：8598.816→V：40032.687→O：369°29’10.2”→G：524.34 →H：2500→P：2500→R：－1 →Q：IfEnd19-If S ≥40557.027:Then9570.142→U：8078.824→V：40557.027→O：257°28’09.1”→G：295.341 →H：1×1045→P：1×1045→R：0→Q：IfEnd20-If S ≥40852.366:Then9506.063→U：7790.519→V：40852.366→O：257°28’09.1”→G：130 →H：1×1045→P：800→R：1 →Q：IfEnd21-If S ≥340982.366:Then9481.312→U：7662.936→V：40982.366→O：262°07’28.12”→G：169.903 →H：800→P：800→R：1 →Q：IfEnd22-If S ≥41152.27:Then9467.011→U：7493.435→V：41152.27→O：274°17’35.205”→G：130 →H：800→P：1×1045→R：1 →Q：IfEnd以上为第二个数据子程序23-SJ－GC（子程序名-竖曲线数据库）24-If S＜39537.6（变坡终点里程）And S≥39062.4（变坡起点里程）：Then0.02 98（大里程坡度）→A：0.01（小里程坡度）→B：39300（变坡点里程）→O：379.858（变坡点高度）→G：－24000（竖曲线半径，凸为正号，凹为负号，直线过度纵坡段半径为1×1045）→R：Prog“SUB3－GC”（高程子程序）:If End25-If S＜40048.9And S≥39537.6：Then0.0298→A：0.0298→B：39537.6→O：386.938→G：1×1045→R：Prog“SUB3－GC”:IfE nd26-If S＜40551.1And S≥40048.9：Then－0.026→A：0.0298→B：40300→O：409.658→G：9000→R：Prog“SUB3－GC”:IfEnd27-If S＜40924.5And S≥40551.1：Then－0.026→A：－0.026→B：40551.1→O：403.129→G：1×1045→R：Prog“SUB3－GC”:IfEnd28-If S＜41275.5And S≥40924.5：Then－0.013→A：－0.026→B：41100→O：388.858→G：－27000→R：Prog“SUB3－GC”:IfEnd29-If S＜41670And S≥41275.5：Then－0.013→A：－0.013→B：41275.5→O：386.576→G：1×1045→R：Prog“SUB3－GC”:IfEnd以上为第三个数据子程序30-SJ－CQW (子程序名-隧道超欠挖数据库)31-Fix 3 ：“H1=”?F（输入实测的高程）：“CQHD=”?C（输入超欠厚度）：5.6+C（圆心1的半径，拱顶圆弧）→R：8.1+C（圆心2的半径，左右两个）→P：1+C（圆心2的半径，左右两个）→M：15+C（圆心2的半径，左右两个）→N：F－H→F：Abs（6.8－Z）→Z（该处的0为数值，即道路中心线到圆心1的竖轴线的水平距离）32-If F＞1.4（圆心1的弧线到设计路面高程的距离）:Then √(Z2+(F－1.4)2)－R（该处的0.980为圆心1点到设计路面高程的距离，圆心点高于设计标高该处为－号，反则为+号）→W（即计算出来的超欠挖值，正为超挖，负为欠挖）: IfEnd:33-If F≥－0.161（即为2.3两个左右圆弧的下限底点到设计高程的距离）AndF ＜1.4（圆心1的弧线到设计路面高程的距离）: Then √((Z+2.1)2+( F－1.4)2)－P→W :IfEnd:（该处的+1.4722为圆心点2.3到隧道中心线的水平距离，注意该处的符号判别，左边的圆弧圆心重到中心线右边的距离为+号，反则为－号；0.13为圆心点2.3到设计路面高程的距离，圆心点高于设计标高该处为－号，反则为+号）34-若If F≥－0.916（即为2.3两个左右圆弧的下限底点到设计高程的距离）And F＜-0.161（圆心1的弧线到设计路面高程的距离）: Then √((Z-4.4669)2+( F +1.3683)2)－M→W :IfEnd:（该处的+1.4722为圆心点2.3到隧道中心线的水平距离，注意该处的符号判别，左边的圆弧圆心重到中心线右边的距离为+号，反则为－号；0.13为圆心点2.3到设计路面高程的距离，圆心点高于设计标高该处为－号，反则为+号）35-If F<－0.916（即为2.3两个左右圆弧的下限底点到设计高程的距离）: Then √((Z-0)2+(- F+13.3)2)－N→W :IfEnd:（该处的+1.4722为圆心点2.3到隧道中心线的水平距离，注意该处的符号判别，左边的圆弧圆心重到中心线右边的距离为+号，反则为－号；0.13为圆心点2.3到设计路面高程的距离，圆心点高于设计标高该处为－号，反则为+号以上为第四个数据子程序36-SUB1－ZS (正算子程序)1÷P→ C：(P-R)÷(2HPR) → D：180÷π→ E：0.1739274226→ A：0.326 0725774→ 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+QE FW(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+Zcos(F)→ X：Y+Zsin(F)→ Y以上为第五个正算子程序37-SUB2－FS (反算子程序)38-G-90→T：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：Got o 0：IfEndLbl 1：0→ Z：Prog "SUB1-ZS"：(J-Y)÷sin(F)→ Z以上为第六个反算子程序39-SUB3－GC （高程子程序）O-S→L:A-B→W:Abs(R*W÷2) →T:O-T→M:O+T→PIf S≤M:ThenG-L*B→H:Goto5:Else If S≤O:Then Goto3:Else If S≤P:Then Go to4:IfEnd: IfEnd: IfEndLbi3If W＞O：Then G+(M-S)2÷2÷R-L*B→H:Goto5: Else If W＜O: Then G-(M-S)2÷2÷R-L*B→H:Goto5: IfEnd: IfEndLbi4If W＞O：Then G+( S-P)2÷2÷R-L*A→H:Goto5: Else If W＜O: Then G-( S -P)2÷2÷R-L*A→H:Goto5: IfEnd: IfEndLbi5H→H: Return以上为第七个高程子程序40-QHZBJS（桥涵坐标计算子程序）41-Lbl 0：42-“ZXJL(+D,－X)=”?A（中线上的距离即涵洞宽度尺寸）：“BXJL（＋Z,－Y）=”?B（中线到边线上的距离即涵洞长度尺寸）：43-180°－J→I：44-180°－I－90→M：45-A÷cos（M）→T：46-X+cos（F）×T→C（该处是计算的中线上的坐标点，改为“XZ=”：X+cos（F）×T→C◢即可显示出来）：47-Y+sin（F）×T→D（该处是计算的中线上的坐标点，改为“YZ=”：Y+sin（F）×T→D◢即可显示出来）：48-“XB=”（显示计算的边桩）：C+cos（F－I）×B→H◢49-“YB=”（显示计算的边桩）：D+sin（F－I）×B→V◢50-0→D："JIXU"?D（是否继续运算输入0继续运算，输入1结束桥涵坐标的计算）：D=0＝﹥Goto 0以上为第八个桥涵子程序51-说明：(正算1秒，反算和超欠挖5秒完成)所有程序名。

自编卡西欧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)等曲线要素赋予定值，可减少在工作中的输入量。

FX5800计算器公路测量常用程序集2.2 版一、程序功能主要功能：采用线元法与交点法相结合计算多条线路坐标正反算，可算任意复杂线型及立交匝道，包括C型，S型、卵型、回头曲线等；极坐标放样，全线路基边坡开挖口及坡脚放样计算、路基任意点设计高程、横坡、设计半幅宽度.隧道欠超挖放样计算。

新版本优化：1、调整一次显示结果；2、交点法中考虑了不对称缓和曲线；3、修改原版本部分地方笔误.(红色为修改处)二、源程序1.总主程序（1、坐标计算放样程序2、坐标反算程序；3、高程计算查阅程序；4、路基半幅标准宽度查阅程序；5、路基边坡及开挖口放样程序；6、路基标准距离放样；7、桥梁锥坡计算放样程序；8、极坐标计算程序；9、隧道超欠挖计算程序）运行后输入1~9数子则选择1至9的程序，返回时，在桩号输入-1,返回选择选择计算类型。

输入-2,返回选择线路。

坐标计算中输入-3,则显示本段曲线要素。

程序名：0ZCXLbI 0:17→DimZ：Norm 2:1→A：＂A:XY=1，ZD=2 ，GC=3，GD=4，BP=5，FM=6，ZP=7，JS=8，SD=9＂?A:A=1=>Goto 1:A=2=>Goto 2:A=3=>Goto 3:A=4=>Goto 4:A=5=>Goto 5:A=6=>Goto 6:A=7=>Goto 7:A=8=>Goto 8:A=9=>Goto 9：LbI 1:Prog ＂DX＂:LbI A:Prog＂QX＂:90→B: ＂PJ1＂?B:B→C: ＂PJ2＂?C:B→Z[1]:C→Z[8]:LbI B:1→F: ＂KM＂?Z:Z= -1=>Goto 0:Z=-2=>Goto A:Z=-3=>Goto X:Prog＂KM＂:?D:Porg＂THB＂:O→L: ＂L0＂?L:Z[2]+Z[1]-Z[8] →E:X+L cos(E) →X:Y+Lsin(E) →Y:Prog＂XY＂:Prog＂JS＂:Goto B：LbI 2:2→F:90→Z[1]:Prog＂QX＂:LbI C: ＂KM＂?Z:Z=-1=>Goto 0:Z=-2=>Goto 2:Z=-3=>Goto X:Prog＂KM＂: ＂X O＂?X: ＂Y0＂?Y:Prog＂THB＂:Porg＂ZD＂:Goto C:LbI 3:Prog＂QX＂:0→B: ＂H-B＂?B:B→Z[9]:LbI D: ＂KM＂? Z:Z=-1=>Goto 0:Z=-2=>Goto 3:?D:Porg＂H＂:Fix 3: ＂ H=＂: Lcoate 6,4,H-Z[9]：＂ I=＂:Locate 6,4,I： Goto D：LbI 4:Prog＂QX＂:LbI E: ＂KM＂?Z:Z=-1=>Goto 0:Z=-2=>Got o 4:?D:Prog＂GD＂:Fix 3: ＂ SJGD=＂:Locate7,4,L： Goto E：LbI 5:Prog＂QX＂:0.5→B:＂TH-GD＂?B:B→Z[17]:LbI F:2→F: 90→Z[1]:＂KM＂?Z:Z=-1 =>Goto 0:Z=-2 =>Goto 5:Prog＂KM＂: ＂X0＂?X: ＂Y0＂?Y:0→M:”M0”?M: M→Z[4]:Prog＂3FB X＂：Fix 2: P→D ：Abs（D）-S→O: ＂ L0=＂:Locate 6,4,O：Prog＂ZD＂: Z[4]→T:＂ TW=＂:Lcoate 6,4,T◢ Goto F:LbI 6:Prog＂DX＂:LbI G:Prog＂QX＂:LbI H:1→F:90→Z[1]:＂KM＂?Z:Z=-1=>Goto 0:Z=-2=>Goto G:Prog＂KM＂:?D:Prog＂THB＂:Prog＂XY＂:Prog＂JS＂:Prog＂H＂:0→M: ＂M0＂?M:Fix 2:H-M→T: ＂ TW=＂:Locate 6,4,T◢ Goto H:LbI 7:Prog＂DX＂:LbI I:Prog＂QX＂:0→D: ＂LD:Z-,Y+＂? D:Abs(D) →R: ＂LR＂?R: ＂Z0＂?Z: D→Z[6]:R→Z[7]:Z→Z[8]:LbI J:0→Y: ＂L0:SZ+,DZ-＂?Y:Y=-1=>Goto 0:Y=-2=>Goto I:Prog＂ZP＂:Goto J:LbI 8:Prog＂DX＂:LbI K:?X:X=-1=>Goto 0:X=-2=>Goto 8:?Y:Prog＂JS＂:Goto K:LbI 9:Porg＂DX＂:LbI L:Prog”QX”:LbI M:2→F:90→Z[1]: ＂KM＂?Z:Z=-1=>Goto 0:Z=-2=>Goto L:Z=-3=>Goto X: Prog＂SD1＂: 0.5→A:“CQHD＂?A:“X0＂?X:＂Y0＂?Y:0→M=＂M0＂?M:M→Z[8]:A→Z[15]:Prog＂5SD＂： Goto M:LbI X:Mat F◢ Goto 02、次程序:路基开挖边线及填方坡脚线放样程序（输入填方放宽值、大概桩号及测量坐标、地面标高计算出偏移距离、桩号、距中距离、填挖高度）程序名：3FBXProg”THB”:D→P:Prog”GD”:If D<0:Then 0.75-L→D:Goto H:Else L-0.75→D:Goto H:IfEndLbI H:Prog“H”:H-0.03-Z[4] →Z[4]: Prog “W1”:If Z[4] <0:Then –Z[4] →G:Goto W:Else Z[4] →G:Goto T:LbI W:Z[8]+Z[9]→A: If G>A:Then Goto 1:Else If G>Z[8]:T hen Goto 2:Else Goto 3:IfEnd:LbI 1:L+Z[10]+Z[11]+Z[12]+(G-A)×Z[7]+Z[9]×Z[6]+Z[8]×Z[5] →S:Goto Z:LbI 2:L+Z[10]+Z[11]+(G-Z[8])×Z[6]+Z[8]×Z[5] →S:Goto Z:LbI 3:L+Z[10]+G×Z[5] →S:Goto Z:LbI T:L+Z[17]→T:If G>Z[15]:Then (T+Z[16]+(G-Z[15])×Z [14]+Z[15]×Z[13])→S:Goto Z:Else (T+G×Z[13])→S:Goto Z:LbI Z3.次程序3.隧道超欠挖值计算放样程序（输入隧道线路，大概桩号、输入衬砌厚度、测量三维坐标，计算准确桩号及位置、计算欠超挖值）程序名：5SDProg”KM”:Prog“THB”: D→Z[13]:0→D:Prog“H”:Cos(Z[1 0]) ×Z[4] →E: H+Z[5] →Z[5]: H+Z[7] →Z[7]: E+Z[5] →Z[10]:If Z[8]>Z[10]:Then Goto R:Else Goto S:IfEnd:LbI R:√((Z[13]-Z[9]) 2+(Z[8]-Z[5]) 2 )- Z[4]-Z[15] →L:G oto L:LbI S:Z[4]-Z[6] →S:√(S 2-(Z[7]-Z[5]) 2 ) ×S÷A bs(S)→T:Abs(Z[13]-Z[9])-T→T:√(T 2+(Z[8]-Z[7]) 2 )- Z[6] -Z[15] →L:Goto L:LbI L: Fix 2:” L0=”: Locate 6,4,L：Prog”ZD”:“ H0=”: Locate 6,4,Z[8]-H◢4：极坐放样计算程序(计算放样点至置仪点方位角及距离)程序名：JSX：Y：Z[11]→K：Z[12]→L：Pol(X-K, Y-L):IF J<0:Then J+3 60→J:IFEnd:Int(J)+0.01Int(60Frac(J))+0.006Frac(60Frac (J)) →J: Fix 4:” FWJ=”: Locate 6,4,J：F ix 3:” S=”:Locate 6,4,I◢5．路基宽度子程序6：路基标准半幅宽度计算程序(对于设计有加宽渐变的有用，如路基宽度无变化，则把此程序直接输入半幅宽度值至L)程序名： GD1→S: Prog “G1”:Z-C→E:(B-A)*E/S+A→L:6．坐标计算次程序（ＴＨＢ）程序名：THBIf Q>1:Then Goto J: Else 1÷P→Ｃ：(P-R)÷(2HPR)→Ｓ：180÷Π→E：If F=1:Then Abs(Ｚ－0) →W：Prog "Ａ"：Ｇoto 2:Else Ｘ→I：Ｙ→J：Prog ＂Ｂ＂：Ｏ＋Ｗ→Z：Ｄ→Ｄ：Ｇoto 2: IfEnd:LbI J: If F=1:Then Prog "Z"：Ｇoto 1:Else Prog ＂ZX＂：Ｇoto 2: IfEnd:LbI 1:I+D×cos(Z[2]+Z[1]) →X: J+D×sin(Z[2]+Z[1]) →Y: LbI 27. 线元法正算子程序(Ａ)程序名：A0.1184634425→A： 0.2393143352→B：0.2844444444→N：0.0 46910077→K： 0.2307653449→L：0.5→M： U+W(Acos(G+QEKW (C+KWＳ))+Bcos(G+QELW(C+LWＳ))+Ncos(G+QEMW(C+MWＳ))+Bco s(G+QE(1-L)W(C+(1-L)WS))+Acos(G+QE(1-K)W(C+(1-K)WS))) →X：V+W(Asin(G+QEKW(C+KWＳ))+Bsin(G+QELW(C+LWＳ))+Nsin(G+QE MW(C+MWS))+Bsin(G+QE(1-L)W(C+(1-L)WS))+Asin(G+QE(1-K)W (C+(1-K)WＳ))) →Y：G+QEW(C+WＳ)→Ｚ［２］：X+Ｄcos（Ｚ［２］+ Z[1]）→Ｘ：Y+Ｄsin（Ｚ［２］+ Z[1]）→Ｙ8. 线元法反算子程序(Ｂ)程序名：BG-90→Ｔ： Abs((Y-V)cos(T)-(X-U)sin(T)) →Ｗ： 0→Ｄ：L bl 0：Prog "Ａ"： T+QEW(C+WＳ) →Ｌ： (J-Y)cos(L)-(I-X) sin(L)→Ｄ：IF Abs(Ｄ)<0.01:Then Goto 1：Else W+Ｄ→W:G oto 0: IfEnd:←┘Lbl 1:0→D：Prog "Ａ"：(J-Y)÷sin(Ｚ［２］+90) →D9. 交点法正算子程序(Z)程序名：ZH2÷R÷24-H∧(4)÷2688÷R∧(3)→A:H÷2-H∧(3)÷240÷R2→B:（（H2-N2）÷24÷R）÷Sin(Abs(P))-((H∧(4)-N∧(4))÷2688÷R ∧(3))÷sin(Abs(P))→E:(R+A)tan(Abs(P)÷2)+B-E→T:P÷Abs(P) →W:0→M:H→C:If Z≤O-T:Then Z-O→S:G→Z[2]:Goto 2: IfEnd:If Z≤O-T+H:Then Z-O+T→S:Prog “HX”:G+WK→Z[2]:Goto 4:IfEnd:If Z≤O-T+ΠR×Abs(P)÷180+H÷2-N÷2: Then 180(Z-O+T-0. 5H)÷R÷Π→S : A+R(1-Cos(S))→B：H÷2-H∧(3)÷240÷R2+Rsin(S)→A:R→M:G+WS→Z[2]:Goto 4: IfEnd:O-T+ΠR×Abs(P)÷180+H÷2+N÷2-Z→S:(R+N2÷R÷24-N∧(4)÷2688÷R∧(3))tan(Abs(P)÷2)+N÷2-N∧(3)÷240÷R2+E→T:N →H:Prog “HX”:G+P→S:S-WK→Z[2]:U+(T-A)Cos(S)-WBSin(S)→I:V+(T-A)Sin(S)+WBcos(S)→J:Goto 3:LbI 4:U+(A-T)cos(G)-WBsin(G)→I:V+(A-T)Sin(G)+WBcos(G)→J: Goto 3:LbI 2:U+Scos(Z[2])→I:V+Ssin(Z[2]) →J: LbI 3:C→H10. 交点法缓和段转化子程序(HX)程序名：HXS-S∧(5)÷40÷R2÷H2+S∧(9)÷3456÷R∧(4)÷H∧(4)→A:S∧(3)÷6÷R÷H-S∧(7)÷336÷R∧(3)÷H∧(3)+S∧(11)÷42240÷R∧(5)÷H∧(5)→B:90S2÷Π÷R÷H→K:RH÷S→M11. 交点法反算子程序(ZX)程序名：ZXZ:0→D:LbI 0:Prog “Z”:Pol(X-I,Y-J):J-Z[2] →J:Isin(J) →S:Icos(J) →I:If Abs(I)<0.1:Then Z+I→Z:S→D:Goto 2:E lse Goto 1: LbI 1:If M=0:Then Z+I→Z:Goto 0:Eles Pol(M-WS,I):(JMΠ)/180→I:Z+I→Z:Goto 0:IfEndLbI 212．高程计算子程序（H）程序名：HP rog “S1”:C-T→F:Z-F→S:C+T→E:G-TI→Q:If T=O:Then Q+ SI→H:Goto 0:Else If Z<F:Then Q+SI→H:Goto 0:Else If Z ≤E:Then Q+SI+S2÷2÷R→H:Goto 0:LbI 0:H:If D=0:Then Goto I:Else Prog “I”:H+V→H:Goto I:LbI I13.高程超高计算程序（I）程序名：IIf Z[3]=1:Then Prog “I1”:Goto 1:Else Prog “I2”:Goto 1: IfEnd:LbI 1: If W=1:Then Goto Z:Else Goto X: IfEndLbI Z:If S=0:Then Abs(D)×M→V:Goto 2:Else Abs(D)×((N-M)×(Z-C)÷S+M)→V:Goto 2:IfEnd:LbI X:If S=0:Then Abs(D)×M→V:Goto 2:Else Abs(D)×(((3 ((Z-C)÷S)2-2((Z-C)÷S)∧(3))×(N-M))+M)→V:Goto 2:IfEn d:LbI 2:Abs(D)→E:V÷E→I:I(E-K)→V14．导线点子程序（DX）程序名：DXZ[11] →K:Z[12] →L:“X Z”?K:”YZ”?L:K→Z[11]:L→Z[12] 15．线路选择子程序(线路选择输0时。

5800隧道三心圆超欠挖程序点击：176 次添加日期：2012-07-22 [ 打印] [ 收藏] [ 关闭]【举报老板技巧】充分掌握黑心老板拒发薪水证据,绝对保护举报者的身份信息，请加QQ：55136901、ECDM(二衬断面检查子程序，结合后面附图理解)“SCGC”?→CIf C<H+0.3:Then √((H+1.1-C)2+(Abs(Z[4]-2.5)-3.932)2)-2.5→A:Abs(Z[4]-2.5)-(3.932+√(2.52-(H+1.1-C)2)→B“R=”:A◢“B=”:B◢Else √((C-H-2.43)2+(Z[4]-2.5)2)-6.7→A: “R=”:A◢IfEnd（所用字母：C H A B）2、FSFYCX（反算放样程序：输入实测点X、Y坐标及实测高程C，求隧道超欠挖值或边坡开挖值）Lb1 B“SCX”?E: “SCY”?Q:E→Z[1]:Q→Z[2]:提示输入实测X、Y坐标Prog “ZB”Lb1 C(X-Z[1]) sin(U-90)-(Y-Z[2])cos(U-90)→Z[10]：试算点与实测点的连线在线路中线上的投影长度，当计算精度小于或等于0.001即1mm时输出计算结果，否则继续计算直至满足精度为止。

If Abs(Z[10])≤0.001:Then K+Z[10]→K:Pol(Z[1]-X,Z[2]-Y):”Z[4]=”:-Icos(J-U+90)→Z[4]◢“K=”:K◢Prog “GCJS”: Prog “ECDM”:Goto B:Else K+Z[10]→K: Prog“ZB”:Goto C:IfEnd(所用字母：E Q X Y U K)3、GCJS（高程计算子程序）26+6D→Z:Lb1 1:Prog “SQXYS”:K-O→L:If K≤O:Then F+AL→H:Else If K≤O+2T:Then F+AL+WLL÷R÷2÷Abs(W)→H:Else Z+3→Z:If Z≤23+6D+3G:Then Goto 1:Else Z-3→Z:Prog “SQXYS”:Z[Z+4]+B(K-O-T)→H:IfEnd(所用字母：Z O K F A L H T R W D)4、GCJSZCX（高程计算主程序）Lb1 C“K=”? →K:Prog “GCJS”:”H=”:H◢Goto C(所用字母：K H)5、SQXYS（竖曲线要素计算子程序）(Z[Z+4]-Z[Z+1])÷(Z[Z+3]-Z[Z])→A:(Z[Z+7]-Z[Z+4])÷(Z[Z+6]-Z[Z+3])→B:B-A→W:Z[Z+5]→R:Abs(0.5RW)→T:Z[Z+3]-T→O:Z[Z+4]-AT→F:A、B：竖曲线前后坡度计算，其余请参看竖曲线计算公式(所用字母：A B W T F Z R)6、KWDM（开挖断面计算子程序，结合后面附图理解）“KWBJ-R”? →R:”SCGC”?C:If C≤H+2.43-Rsin(13°10°9°):Then Abs(Z[4]-2.5)-(R-0.45(H+2.43-RSin(13°10°9°)-C) ÷2.1)→B:”B=”:B◢Else If C≤H+2.43+R:Then√((Z[4]-2.5)2)+(C-H-2.43)2)-R→A:Abs(Z[4]-2.5)-√(R2-(C-H-2.43)2)→B:C-H-2.43-√(R2-(Z[4]-2.5)2)→V:“R=”:A◢“B=”:B◢“V=”:V◢Else C-H-2.43-√(R2-(Z[4]-2.5) 2)→V: “V=”:V◢IfEnd（所用字母：R C H A B V）7、SJK 30（平曲线及竖曲线数据库程序，投影面高程为30米，红色部分为竖曲线）7→D: 10→G :30+6D+3G→DimZ:10→Z[1]:3433109.5170→Z[20]:3429407.741→Z[21]:3426833.089→Z[22]:3424404.600→Z[23]:3421023.289→Z[24]:3417221.367→Z[25]:3414448.180→Z[26]:3403462.4596→Z[27]:3405098.748→Z[28]:502494.2176→Z[29]:502639.171→Z[30]:503107.984→Z[31]:503310.589→Z[32]:503948.482→Z[33]:504343.928→Z[34]:505035.224→Z[35]:509994.4824→Z[36]:509997.8856→Z[37]:9000→Z[38]:9000→Z[39]:9000→Z[40]:10000→Z[41]:9000→Z[42]:9000→Z[43]:7000→Z[44]:490→Z[45]:490→Z[46]:490→Z[47]:430→Z[48]:490→Z[49]:490→Z[50]:670→Z[51]:490→Z[52]:490→Z[53]:490→Z[54]:430→Z[55]:490→Z[56]:490→Z[57]:670→Z[58]:1139000→Z[59]:1141824.092→Z[60]:1144638.005→Z[61]:1147045.518→Z[62]:1150566.149→Z[63]:1154138.336→Z[64]:1156817.127→Z[65]:1170598.039→Z[66]:1173840.340→Z[67]:1138090→Z[68]:77.040→Z[69]:1140000→Z[71]:71.310→Z[72]:30000→Z[73]:1141600→Z[74]:61.710→Z[75]:30000→Z[76]:1142900→Z[77]:45.460→Z[78]:30000→Z[79]:1143950→Z[80]:42.100→Z[81]:30000→Z[82]:1145400→Z[83]:47.755→Z[84]:30000→Z[85]:1147800→Z[86]:51.355→Z[87]:30000→Z[88]:1149300→Z[89]:44.005→Z[90]:30000→Z[91]:1153500→Z[92]:53.245→Z[93]:30000→Z[94]:1156300→Z[95]:72.125→Z[96]:30000→Z[97]:1156300→Z[98]:79.325→Z[99]:30000→Z[100]:1157980→Z[101]:49.085→Z[102]（所用字母：D、G）8、SJK90（平曲线及竖曲线数据库程序，投影面高程为90米，红色部分为竖曲线）2→D: 2→G :30+6D+3G→DimZ:2→Z[1]:44381.669→Z[20]:39644.955→Z[21]:36290.557→Z[22]:33141.9901→Z[23]:1072.7059→Z[24]:2635.9640→Z[25]:2370.9500→Z[26]:2494.2410→Z[27]:7000→Z[28]:9000→Z[29]:670→Z[30]:490→Z[31]:670→Z[32]:490→Z[33]:1127535.872→Z[34]:1130778.084→Z[35]:1135073.865→Z[36]:1139000→Z[37]:1132900→Z[38]:102.26→Z[39]:1134650→Z[41]:67.260→Z[42]:30000→Z[43]:1138000→Z[44]:77.310→Z[45]:30000→Z[46]:1139955→Z[47]:71.445→Z[48]:（所用字母：D、G）9、SWKZSJK5（石武客专数据库，红色部分为竖曲线，）7→D: 8→G: 30+6D+3G→DimZ:71955.5972→Z[20]:71428.348→Z[21]:64928.597→Z[22]:60648.194→Z[23]:58890.534→Z[24]:56131.938→Z[25]:50547.524→Z[26]:47965.040→Z[27]:47303.7078→Z[28]:652.0165→Z[29]:709.292→Z[30]:1017.189→Z[31]:938.295→Z[32]:812.152→Z[33]:813.346→Z[34]:462.228→Z[35]:456.788→Z[36]:517.0987→Z[37]:12000→Z[38]:11000→Z[39]:12000→Z[40]:12000→Z[41]:12000→Z[42]:12000→Z[43]:9000→Z[44]:330→Z[45]:370→Z[46]:330→Z[47]:330→Z[48]:330→Z[49]:330→Z[50]:490→Z[51]:330→Z[52]:370→Z[53]:330→Z[54]:330→Z[55]:330→Z[56]:330→Z[57]:490→Z[58]:1099872.681→Z[59]:1099872.681→Z[60]:1106362.974→Z[61]:1110706.287→Z[62]:1112355.010→Z[63]:1115166.386→Z[64]:1120776.803→Z[65]:1123224.207→Z[66]:1124551.665→Z[67]:1104639→Z[68]:113.4094→Z[69]:1104900→Z[71]:109.390→Z[72]:30000→Z[73]:1105900→Z[74]:111.390→Z[75]:30000→Z[76]:1106900→Z[77]:93.390→Z[78]:25000→Z[79]:1109400→Z[80]:73.390→Z[81]:30000→Z[82]:1112050→Z[83]:98.035→Z[84]:25000→Z[85]:1113900→Z[86]:84.190→Z[87]:30000→Z[88]:1116500→Z[89]:113.310→Z[90]:30000→Z[91]:1117700→Z[92]:110.910→Z[93]:30000→Z[94]:1122700→Z[95]:165.910→Z[96]:（所用字母：D、G）10、XLZB（线路坐标计算主程序）“1.ZX 2.BX”?→O: 提示选择：计算线路中线还是边线0→Q：Lb1 1If O≠1:Then “L”? →Q:”PJ”? →J:IfEnd：如选择边线提示输入宽度及偏角，宽度始终输正值，偏角右法线输+90，左法线输-90Prog “ZB”X+Rec(Q,U+J) →XY+J→Y“X=”:X◢“Y=”:Y◢If O=1:Then “QXFWJ=”:U▲DMS◢IfEnd：当计算中线时，会显示该点的切线方位角Goto 111、PQXYS（平曲线要素计算子程序）Z[V+2D+3] →R (调入半径)Z[V+3D+3] →L（调入第一缓和曲线长度）Z[V+4D+3] →M（调入第二缓和曲线长度）Pol(Z[V+1]-Z[V],Z[V+D+3]-Z[V+D+2])（交点前方位角，前进方向）J→A:Pol(Z[V]-Z[V-1],Z[V+D+2]-Z[V+D+1])（交点后方位角，前进方向）A -J→C:J→A:If C<-180:Then 1→N:ElseI if C<0: Then -1→N :Else if C<180:Then 1→N:Else -1→N:IfEnd （左右转参数计算，N=1右转，N=-1左转）If Abs(C) >180:Then 360-Abs(C) →C :If ENDAbs(C)→CL÷2-L^(3)÷240÷R2→Z[16]M÷2-M^(3)÷240÷R2→Z[17]L2÷24÷R→Z[18]M2÷24÷R→Z[19](Z[19]-Z[18])÷2÷tan(C÷2)→P(Z[18]+Z[19]+2R)tan(C÷2)÷2→Z[7]P+Z[7]+Z[16]→TZ[7]-P+Z[17]→PπRC÷180+(L+M) ÷2→Z[7]注：其余公式请参照不等长缓和曲线计算公式12、ZB（坐标计算子程序）Lb1 D22→W:?KLb1 2If K<Z[24+5D] Or K>Z[25+6D]:Then “OUT OF RANGE”◢Goto D: IfEnd:(里程计算超限提示)If K≤Z[W+5D+4]:Then W-1→V:Else W+1→W:Goto 2:IfEnd:（判断所输入里程对应的交点编号）Prog “PQXYS”Z[V]-Rec(T,A)→Z[14]Z[V+D+2]-J→Z[15]If K>1138800.074 And K<1139000:Then K+0.074→K:IfEnd（断链处理）If K=1139000 And D=2:Then K+0.074→K:IfEnd（断链处理）K-Z[V+5D+4] →SA→BIf S≤0:Then Z[14]+Scos(B)→X：（第一直线段坐标计算）Z[15]+Ssin(B) →YB→U:Goto 7:IfEndZ[14]→X:Z[15]→YIf S≤L:Then 1→H:Goto 4:IfEnd：（第一缓和曲线段坐标计算）If S≤Z[7]-M:Then 1→H:Goto 5:IfEnd：（圆曲线段坐标计算）A+NC→BIf S≤Z[7]:Then Z[7]-S→S:Z[V]+Rec(P,B)→X:Z[V+D+2]+J→Y:-1→H:M→L:Goto 4:Z[V+3D+3] →L:IfEnd：（第二缓和曲线段坐标计算）Z[V]+Rec(P+S-Z[7],B) →X:Z[V+D+2]+J→Y:B→U:Goto 7：（第二直线段坐标计算）LB1 4（缓和曲线计算公式）S^(3)÷6÷R÷L-S^(7)÷336÷R^(3)÷L^(3)→IS-S^(5)÷40÷R2÷L2→JB+90NHS2÷π÷R÷L→U:G oto 6Lb1 5（圆曲线计算公式）90(2S-L) ÷π÷R→UR(1-cos(U))+Z[18]→IRsin(U)+Z[16]→JB+180N(S-L÷2)÷π÷R→ULb1 6Pol(J,I):NHJ+B+90(H-1)→JRec(I,J)+X→XY+J→YLb1 7If U<0:Then U+360°→U:IfEnd13、ZBFS（坐标反算：已知两点坐标求距离及方位角）Lb1 0“X1”?X: “Y1”?Y:”X0”?D:”Y0”?G:Pol(X-D,Y-G):If J<0:Then J+360→J:IfEnd“JL=”:I◢“FWJ=”:J▲DMS◢Goto 014、ZBZS（坐标正算：已知起点坐标、距离及方位角求终点坐标）Lb1 0“X0=”?C:”Y0=”?D:”JL=”?E:”FWJ=”?F:Rec(E,F):C+I→X:D+J→Y:“X=”:X◢“Y=”:Y◢Goto 015、ZSFYCX（正算放样程序：已知里程、宽度求距离及方位角）“ZJX=”?→Z[1]：提示输入置镜点X坐标“ZJY=”?→Z[2]：提示输入置镜点Y坐标Lb1 A“FY:1.ZX 2.BX”?→O：提示选择是放样中线还是边线，1为中线、2为边线If O≠1:Then “KD=”?→Q:”PJ=”?→J:IfEnd：放样边线提示输入宽度及偏角Prog “ZB”X+Rec(Q,U+J)→XY+J→YPol(X-Z[1],Y-Z[2])If J<0:Then J+360→J:IfEnd“S=”:I◢“A=”:J▲DMS◢Goto A三心圆开挖放样程序“ZDKW”（也可用于多心圆、同心圆）：请结合后面所附示意图理解！｛CR｝:C“△H”:Z[9]+C≤H+h1+RSinα＝＞A=√((Abs(Z[4]±Z[8]×?)+d)2+(Z[9]+C－H－h2)2) －R－r1▲B=Abs(Z[4] ±Z[8]×?)+d－√(（R+ r1）2－(Z[9]+C－H－h2)2) ▲V=Z[9]+C-H-h2-√((R+r1)2－(Z[4] ±Z[8]×?)2) ▲≠＞Z[9]+C≤H+h1+R ＝＞A=√((Abs(Z[4] ±Z[8]×?) 2+(Z[9]+C－H－h1)2) －R▲B=Abs(Z[4] ±Z[8]×?)－√(R2－(Z[9]+C－H－h1)2) ▲V=Z[9]+C-H-h1-√(R2－(Z[4] ±Z[8]×?)2) ▲≠＞V=Z[9]+C－H－h1－√(R2－Abs(Z[4] ±Z[8]×?) 2) ▲△△C (显示为△H)：实测点至仪器中心的高差（如直接输入高程，则将C“△H”改为C“SCGC”，且程序中的Z[9]+C就直接用C代替）R：开挖半径（O1）Z[9]：仪高（直接输实测高程时不用）A：实测点半径与设计半径之差值（径向），+超挖，—欠挖B：实测点与设计点之差值（水平方向），+超挖，—欠挖V：实测点与设计点之竖向距离r1：O1与O2(O3)之差h1：圆心O1至设计高程点的距离,h2：圆心O2（O3）至设计高程点的距离H：实测点设计高程，由高程子程序带入Z：须预先置入扩展变量中，Defm 15d：偏移两圆心（O2、O3）与隧道中线之水平距离？：线路中线与隧道中线之偏移距离Z[8]：左右线判断变量（分离式双线隧道）及隧道左右偏判断变量（相对线路设计中心线）平曲线及竖曲线数据库设置表（n个平交点）1、平曲线数据库设置表注：1、从Z[20]开始存储数据，是为了留出前面的扩展变量作为其他变量使用，以补充变量A-Z之不足，当然也可以从Z[30]、 Z[40]开始，留出更多的变量供使用，只需根据上面表格作相应调整。

隧道超欠挖计算原理与程序设计摘要：随着铁路、公路建设的快速发展，我国隧道建设的施工技术也大幅度提高，由于工期紧迫促成测量放样的过程也是一个重要环节，如果再采用普通的尺距法不仅仅降低了放样效率还造成了轮廓线的放样精度，导致开挖掘进造成隧道超欠挖，所以现在通过卡西欧编程计算器结合带红外线的全站仪进行配套操作，这样不仅提高了施测效率还保证了放样精度，现通过本文给大家讲解隧道超欠挖的计算原理和程序设计。

关键词：隧道超欠挖计算原理程序设计何为超欠挖？：隧道超欠挖分为（超挖和欠挖），超挖即为隧道开挖轮廓线大于隧道设计轮廓线，欠挖即为隧道开挖轮廓线小于隧道设计轮廓线。

超欠挖的影响：隧道超欠挖不止直接影响到了施工进度、安全质量，还会让开挖费用增加，更重要的是由此造成了过量超填混凝土的费用。

超挖在实际施工中由于重视不够或方法不当，以至于在施工过程中会不知不觉地提高工程成本，从而也减少了应得的利润。

超欠挖是如何产生的？在目前的隧道施工中，掘进技术有两种方法，一种是传统的“钻爆法（开挖台阶法）”；一种是“全断面掘进法（盾构掘进法）”。

受各种条件的制约，“钻爆法”仍是隧道施工的主要掘进方法。

所以隧道超欠挖的形成也是不可避免的，下面讲述一下形成超欠挖的三种情况。

1、岩层变化：由于隧道开挖过程中随着岩层的变化，地质条件和围岩裂隙的发生会出现不可避免的超欠现象，所以岩体是超欠挖的主要因素之一。

2、爆破方式：由于工作面（掌子面）是一个不平整的岩体面，导致钻孔间距控制不当或间距过大、过小，容易影响其他孔位的爆破效果，或者由于装药结构控制不当和掏槽不合理也会造成隧道超欠现象。

3、测量放线：由于隧道测量放线过程中能见度低，操作有限，测量人员进入隧道测量时导致前后视照准误差，同时因为掌子面的凹凸不平画轮廓线时也会产生偏离现象。

如何正确实施隧道测量工作？隧道测量工作由公司专业测量人员组织成立测量小组，根据设计院给定的坐标控制点和高程控制点进行建立导线控制网。

隧道开挖断面测量步骤（适用于Casio fx-5800P）步骤如下：1.架站，后视,记下测站坐标（X,Y,H）;2.求出断面桩号，计算断面中心点O0坐标（X0,Y0）、切线方位角V;3.运行程序：程序名：“ZBZH”?A:?B:?X:?Y:W→WPol(A-X,B-Y):Rec(I,J-W)Cls:”N”:Locate 3,1,I:“E”:Locate 3,2,J◢A 测站坐标X值B 测站坐标Y值V 第2步中记下的切线方位角V值X,Y 第2步中记下的洞轴线坐标（X0,Y0）记下所求得的N值与E值;4.转动全站仪，将水平角值调到与切线方位角V相等，然后置零，再将测站坐标（X,Y）值改成第3、4步所求得的（N,E）所对应的值，并将目标高改为0。

（注意！切勿改动仪器高！）5.开始放开挖轮廓点，将所测得的N,E,Z代入以下程序运行：程序名：“DMCL”1.“1.ZD,2.YD”?U:If U=1:Then *→A:* →B:* →C:* →D:* →R:IfEnd2.If U=2:Then *→A:* →B:* →C:* →D:* →R:IfEnd3.R÷cos((tan-1 A- tan-1 B)÷2)-R→G4.√((R+G)2-R2)→T5.Abs(Tcostan-1A)→Q6.Abs(Rsintan-1A)→O7.C-Q→S8.C+ Abs (Tcostan-1B)→K9.“DMLX”?T:”R”?P:”O-H”?X:1.5→L:0.25→M T－断面类型，P－断面半径，X－设计高与圆心高高差，L、M如图所示。

10.Lbl 0:?N:?E:?Z:Abs(E)→E11.If N≤S:Then D-(C-N)A→H:IfEnd12.If N≥K:Then D+( N-C)B→H:IfEnd13.If N<K:Then If A＞B:Then 1→F:If A＜0:Then 1→W:Else -1→ W:IfEnd:IfEnd14.If A＜B:Then -1→F:If A＞0:Then 1→W:Else -1→W:IfEnd:IfEnd15.D-F(√((R+G)2-(O+WQ)2) -√(R2-(O+W(N-S))2))→H:IfEnd16.Z-(H+X)→V :T=1=>Goto 1: T=2=>Goto 217.Lbl 1:If V>0:Then Pol(E,V):"SHANG": P-I→Y◢18.Goto 0:IfEnd19.If V<0:Then Pol(E+L,V):"XIA":P+L-I→Y◢20.Goto 0:IfEnd 正常段21.Lbl 2:If V<0:Then Pol(E+M,V):"XIA":P+L-I→Y◢22.Goto 0:IfEnd23.If V>0:Then Pol(E-L+M,V)24.If J<50°:Then "SHANG":P-I→Y◢25.Goto 0:IfEnd26.If J≥50°:Then Pol(E,V+L-M): "SHANG":P+(L-M)÷cos50°-I→Y◢27.Goto 0:IfEnd 加宽段图一：图二：。

3.GL-LUNKUO 超欠挖主程序（最快版……里程精度要求高）Lbl 0?P:P→L:Prog“GLZUOBIAO-1”↙?X:?Y:?C↙Pol ((X-S),(Y-T)):If V>180:Then V→V:Goto 1:Else V+180→V:Goto 1↙Lbl 1↙“L=”:I×cos(J-V)+L→L◢Prog“GLZUOBIAO-1”↙Pol ((X-S),(Y-T)):If V>180:Then V→V:Goto 2: Else V+180→V:Goto 2↙Lbl 2↙“FAJU=”:I×s in(J-V)→T◢“LICHENG=”:I×cos(J-V)+L→L◢Goto 3↙Lbl 3↙If L<1156000:Then Abs(7.725-T)→T: Goto V:Else Abs(7.725-(-T))→T: Goto V:IfEnd↙Lbl V↙If L≤156300: Then Prog“GC1”: Goto A: Else Goto Z: IfEnd↙Lbl Z↙If L≤160000: Then Prog“GC3”: Goto A: Else Pros“GC2”:Goto A: IfEnd↙Lbl A↙If C≥H+1.55: Then Goto B: Else Goto C: IfEnd↙Lbl B↙If tan- 1((C-H+0.43)÷T)-45≤0: Then Goto D: Else Goto E: IfEnd↙Lbl D↙“GBR=”:√((C-H-1.55)2+(T-1.98)2→R◢ ( √为开方符号)Goto 0↙Lbl E↙“GZR=”:√((C-H+0.43)2+T2 )→R◢Goto 0↙Lbl C↙If tan- 1(Abs(H+1.155-C) ÷(T-1.98))≤-19.80138≤0: Then Goto F: Else Goto G: IfEnd↙Lbl F↙“BQR=” : √((H+1.55-C) 2+(T-1.98)2 ) →R◢Goto 0↙Lbl G↙If tan- 1(T÷(H+23-C))-15.13444≤0: Then Goto H: Else Goto I: IfEnd↙Lbl H↙“YGR=”: √(T2 +(H+23-C) 2→R◢Goto 0↙Lbl I↙T-6.2139→M◢“R=”:√(M2 +(H-C+0.0255) 2)→R◢Goto 0↙CQW超欠挖主程序（精简版……里程精度相对放宽）Lbl 0:Cls:Fix 3:Deg:“L”?L:“X=”?T:“Y=”?S:“H=”?A:Prog“GL-FS”↙“LICHENG”:L◢“PIANJU”:J◢Goto 1↙Lbl 1:If L<1156800:Then Abs(8.975-J) →J:Goto 2: Else Abs(8.975-(-J)) →J:Goto 2: IfEnd↙Lbl 2:L→K:Prog“GCSJK”:Goto 3↙Lbl 3:If A≥H+1.55:Then Goto 4:Else Goto7: IfEnd↙Lbl 4:If tan- 1((A-H+1.68) ÷J)-45≦0: Then Goto 5:Else Goto 6: IfEnd↙Lbl 5:“R=”:√((A-1.55)2+(J-3.23)2→R◢ ( √为开方符号)Goto 0↙Lbl 6:“R=”:√((A-H+1.68)2+J2 )→R◢Goto 0↙Lbl 7: If tan- 1(Abs(H+1.55-A)÷(J-3.23))-20.0944≦0: Then Goto 8:Else Goto 9: IfEnd↙Lbl 8:“R=”:√((H+1.55-A)2+(J-3.23)2→R◢Goto 0↙Lbl 9: If tan- 1(J÷(H+27.622-A))-15.1344≦0: Then Goto 10:Else Goto 11: IfEnd↙Lbl 10:“R=”:√(J2+(H+27.622-A)2 )→R◢Goto 0↙Lbl 11:J-6.2139→M◢“R=”:√(M2+(H-A+0.0255)2 )→R◢Goto 0↙GC1 （高程1？）IF L≦147221:Then 52.403→A:-0.005→I:24000→R:147221→O (非零，为大写o):146997.56→Q:Goto 1:Else Goto A:IfEnd↙LbI 1↙If L<0:Then 0→K: Goto 0: Else 1→K: Goto 0:IfEnd↙LbI A↙IF L≦147444.44:Then 52.403→A:-0.02362→I:24000→R:147221→O (非零，为大写:147444.44→Q:1→K:Goto 0:Else Goto B:IfEnd↙LbI B↙IF L≦148000:Then 70.803→A:-0.02362→I:160000→R:148000→O (非零，为大写o):147574.4→Q:Goto 2:Else Goto C:IfEnd↙LbI 2↙If L<0:Then 0→K: Goto 0: Else -1→K: Goto 0:IfEnd↙LbI C↙IF L≦148425.6:Then 70.803→A:-0.0183→I:160000→R:148000→O (非零，为大写:148425.6→Q:-1→K:Goto 0:Else Goto D:IfEnd↙LbI D↙IF L≦149000:Then 89.103→A:-0.0183→I:50000→R:149000→O (非零，为大写o):148707.5→Q:Goto 3:Else Goto E:IfEnd↙LbI 3↙If L<0:Then 0→K: Goto 0: Else 1→K: Goto 0:IfEnd↙LbI E↙IF L≦149292.5:Then 89.103→A:0.03→I:50000→R:149000→O(非零，为大写:149295.5→Q:1→K:Goto 0:Else Goto F:IfEnd↙LbI F↙IF L≦149900:Then 116.103→A:-0.03→I:22086.957→R:149900→O (非零，为大写o):149646→Q:Goto 4:Else Goto G:IfEnd↙LbI 4↙If L<0:Then 0→K: Goto 0: Else -1→K: Goto 0:IfEnd↙LbI G↙IF L≦150154:Then 116.103→A:0.007→I:22086.957→R:149000→O (非零，为大写:150154→Q:-1→K:Goto 0:Else Goto H:IfEnd↙LbI H↙IF L≦150450:Then 119.953→A:-0.007→I:16000→R:150450→O (非零，为大写o):150154→Q: -1→K:Goto 0:Else Goto I:IfEnd↙LbI I↙IF L≦150746:Then 119.953→A:-0.03→I:16000→R:150450→O (非零，为大写o):150746→Q: -1→K:Goto 0:Else Goto J:IfEnd↙LbI J↙IF L≦151200:Then 97.453→A:-0.03→I:27000→R:151200→O(非零，为大写o):150961.05→Q: Goto 5:Else Goto K:IfEnd↙LbI 5↙If L<0:Then 0→K: Goto 0: Else 1→K: Goto 0:IfEnd↙LbI K↙IF L≦151438.95:Then 97.453→A:-0.0123→I:27000→R:151200→O (非零，为大写o):151438.95→Q: 1→K: Goto 0:Else Goto L:IfEnd↙LbI L↙IF L≦151800:Then 90.074→A:0.0123→I:17070.922→R:151800→O (非零，为大写o):151438.95→Q: 1→K: Goto 0:Else Goto M:IfEnd↙LbI M↙IF L≦152161.05:Then 90.074→A:0.03→I:17070.922→R:151800→O (非零，为大写o):152161.05→Q: 1→K: Goto 0:Else Goto N:IfEnd↙LbI N↙IF L≦152650:Then 115.574→A:-0.03→I:12000→R:152650→O (非零，为大写o):152338→Q: Goto 6:Else Goto O:IfEnd↙LbI 6↙If L<0:Then 0→K: Goto 0: Else 1→K: Goto 0:IfEnd↙LbI O↙IF L≦152962:Then 115.574→A:-0.022→I:12000→R:152650→O (非零，为大写o):152962→Q: -1→K:Goto 0:Else Goto P:IfEnd↙LbI P↙IF L≦153350:Then 100.174→A: 0.022→I:18000→R:153350→O (非零，为大写o):153101.15→Q: Goto 7:Else Goto Q:IfEnd↙LbI 7↙If L<0:Then 0→K: Goto 0: Else 1→K: Goto 0:IfEnd↙LbI Q↙IF L≦153598.85:Then 100.174→A: 0.00565→I:18000→R:153350→O (非零，为大写o):153598.85→Q: 1→K:Goto 0:Else Goto R:IfEnd↙LbI R↙IF L≦155500:Then 112.329→A: -0.00565→I:19000→R:155500→O (非零，为大写o):15587.675→Q: Goto 8:Else Goto S:IfEnd↙LbI 8↙If L<0:Then 0→K: Goto 0: Else 1→K: Goto 0:IfEnd↙LbI S↙IF L≦155712.325:Then 112.329→A: 0.028→I:19000→R:155500→O (非零，为大写o):155712.325→Q: 1→K:Goto 0:IfEnd↙LbI 0↙A+I×Abs(L- O)+K×(L-Q) 2÷2÷R→H↙GC2 （高程2？）IF L≦1156300:Then 134.729→A:-0.028→I:140000→R:1156300→O (非零，为大写o):1155843.6→Q:Goto 9: Else Goto A:IfEnd↙LbI 9↙If L<0:Then 0→K: Goto 0: Else -1→K: Goto 0:IfEnd↙LbI A↙IF L≦1156756.4:Then 134.729→A:-0.02148→I:140000→R:1156300→O (非零，为大写o):1156756.4→Q: -1→K :Goto 0: Else Goto B:IfEnd↙LbI B↙IF L≦1157270:Then 155.569→A:-0.02148→I:25000→R:157270→O (非零，为大写o):1156958.625→Q: Goto 2: Else Goto C:IfEnd↙LbI 2↙If L<0:Then 0→K: Goto 0: Else -1→K: Goto 0:IfEnd↙LbI C↙IF L≦1157581.375:Then 155.569→A:-0.00343→I:25000→R:1157270→O (非零，为大写o):1157581.375→Q: -1→K :Goto 0: Else Goto D:IfEnd↙LbI D↙IF L≦1159800:Then 146.895→A: 0.00343→I:60000→R:115900→O (非零，为大写o):1159122.9→Q: Goto 3: IfEnd↙LbI 3↙If L<0:Then 0→K: Goto 0: Else -1→K: Goto 0:IfEnd↙LbI 0↙A+I×Abs(L- O)+K×(L-Q) 2÷2÷R→H↙GC3 （高程3？）IF L≦156300:Then 134.729→A:-0.028→I:110000→R:156300→O (非零，为大写o):1155830.85→Q:Goto 1: Else Goto A:IfEnd↙LbI 1↙If L<0:Then 0→K: Goto 0: Else -1→K: Goto 0:IfEnd↙LbI A↙IF L≦156769.15:Then 134.729→A:0.01947→I:110000→R:156300→O (非零，为大写o):156769.15→Q: -1→K :Goto 0: Else Goto B:IfEnd↙LbI B↙IF L≦157300:Then 154.199→A:-0.01947→I:25000→R:157300→O (非零，为大写o):157010.25→Q: Goto 2: Else Goto C:IfEnd↙LbI 2↙If L<0:Then 0→K: Goto 0: Else -1→K: Goto 0:IfEnd↙LbI C↙IF L≦157589.75:Then 154.199→A:-0.00371→I:25000→R:157300→O (非零，为大写o):157589.75→Q: -1→K :Goto 0: Else Goto D:IfEnd↙LbI D↙IF L≦159720:Then 145.215→A:0.00371→I:62000→R:159720→O (非零，为大写o):159091.01→Q: Goto 4: Else Goto E:IfEnd↙LbI 4↙If L<0:Then 0→K: Goto 0: Else -1→K: Goto 0:IfEnd↙LbI 0↙A+I×Abs(L- O)+K×(L-Q) 2÷2÷R→H↙GCSJK （高程数据库？）IF K≦159091.01:Then154.199→G:157300→B:25000→D:289.750→U:0.01947→I: -0.00371→W:B-K→C:1→F:I>W =>-1→F:Goto 3:IfEnd↙IF K≦161364.5:Then145.215→G:159720→B:62000→D:628.990→U:-0.00371→I: -0.024→W:B-K→C: I>W =>-1→F:Goto 3:IfEnd↙LbI 3↙IF K<B-U:Then 0→Z:I→P:Goto 4: IfEnd↙IF K<B+U:Then 1→Z:I→P:Goto 4: IfEnd↙IF K<B+U:Then 1→Z:W→P:Goto 4: IfEnd↙IF K>B:Then 0→Z:W→P:Goto 4: IfEnd↙LbI 4:G-CP+Z×F(U-Abs(C)) 2÷2÷D→H↙GL-FS （什么反算？）20→V：Do:Prog“PMSJK”:(P-R)÷(2HPR) →D:L-O→O:OD→D:1÷P→C:If V>10: Then Prog“GL-JY-ZS”:Else Prog“GL-ZS”:IfEnd:T-X→I:S-Y→J:Pol(I,J):J-Q →J:Rec(I,J):I÷(1-CJ-2DJ)→V:L+V→L:Abs(V)→V:LpWhile V>0.001GL-JY-ZS （什么正算？）0.7886751346→E:1-E→F:Q+E r O(C+ED)→E:Q+E r O(C+FD)→F:X+O(cos(E)+cos(F))÷2→X:Y+O(sin(E)+sin(F))÷2→Y:Q+O r(C+D)→QGL -ZS （什么正算？）0.1739274226→I:0.3260725774→J:0.0694318442→K:0.3300094782→E:1-E→F:1-E→F:1-K→M:Q+K r O(C+KD)→K:Q+E r O(C+ED)→E:Q+F r O(C+FD)→F:Q+M r O(C+MD)→M:X+O(Icos(K)+Jcos(E)+ Icos(M))→X:Y+O(Isin(K)+Jsin(E)+ Isin(M))→Y:Q+O r(C+D)→Q:GLZUOBIAO-1 （表示什么？）If L≤154246.087:Then Goto A:Else Goto 1:IfEnd↙LbI 1↙If L≤160350:Then Goto B:Else Goto C:IfEnd↙LbI A↙Prog“YZB”↙Goto 5↙LbI B↙Prog“ZB2”↙Goto 5↙LbI C↙Prog“ZZB”↙Goto 5↙LbI 5↙If W<1:Then V-180→V: Goto 3: Else Goto 3:IfEnd↙LbI 3↙If V<90:Then 360+V→V: Else V→V: IfEnd↙PMSJK (平面数据库?)IFL≦156981.234:Then590888.498→X:533021.975→Y:10^(46)→P:1545.627→R:17 5→H:156806.234→O:289°14′14.5″→Q:Return:IfEnd↙L≦157432.396:Then590949.257→X:532857.886→Y:1545.627→P:1545.627→R:4 51.162→H:156981.234→O:292°28′51.4″→Q:Return:IfEnd↙IFL≦157607.396:Then591179.742→X:532471.904→Y:1545.627→P:10^(46)→R:17 5→H:157432.396→O:309°12′19.3″→Q:Return:IfEnd↙IFL≦158977.376:Then591295.382→X:532340.588→Y:10^(46)→P:10^(46)→R:136 9.98→H:157607.396→O:312°26′56.2″→Q:Return:IfEnd↙YQX （右曲线？）Abs(L-O)→H:S 2÷(24×R)-S^(4)÷(2688×R^(3))→P:S÷2-S^(3)÷(240R2)→Q: (H-S)×180÷R÷∏+S×180÷（2×∏×R→U：K×U+Z→V:Rsin(U)+Q→E:R(1-cos (U))+P→N:“S=”:A+Ecos(Z)-KNsin(Z)→S◢“T=”:B+Esin(Z)+ KNcos(Z) →T◢↙↙Return上面∏为3.1415926……YXSQX （右线竖曲线？）LbI 0↙Cls:“K=”?K:IF K<159091.01:Then 154.199→G：157300→B：25000→R：289.750→T:0.01947→I:-0.00371→J:B-K→C:1-F:I>J=>-1→F:Goto 3:IfEng↙IF K<161600:Then 145.215→G：159720→B：62000→R：628.990→T:-0.00371→I: -0.024→J:B-K→C:I>J=>-1→F:Goto 3: IfEng↙LbI 3: IF K<B-T:Then 0→A:I→ P: Goto 4: IfEng:IF K<B:Then 1→A:I→ P: Goto 4: IfEng:IF K<B+T:Then 1→A:J→ P: Goto 4: IfEng:IF K>B:Then 0→A:J→ P: Goto 4: IfEng↙LbI 4:G-CP+AF(T-Abs(C)) 2÷2÷R→H:“H”:H◢Goto 0↙ZXSQX （左线竖曲线？）LbI 6↙Cls:“K=”?K:IF K<159122.9:Then 155.569→G：157270→B：25000→R：311.375→T:0.02148→I:-0.00343→J:B-K→C:1→F:I>J=>-1→F:Goto 7:IfEng↙IF K<160477:Then 146.895→G：159800→B：60000→R：677.1→T:-0.00343→I: -0.026→J:B-K→C: -1→F:Goto 7: IfEng↙LbI 7: IF K<B-T:Then 0→A:I→ P: Goto 8: IfEng:IF K<B-T:Then 1→A:I→ P: Goto 8: IfEng:IF K<B+T:Then 1→A:J→ P: Goto 8: IfEng:IF K>B:Then 0→A:J→ P: Goto 8: ↙LbI 8:G-CP+AF(T-Abs(C)) 2÷2÷R→H:“H”:H◢Goto 6↙ZB2 （表示什么？）If L≦154426.087:Then Goto A:Else Goto 1:IfEng↙If L≦155115.597:Then Goto B:Else Goto 2:IfEng↙LbI 2↙If L≦155265.597:Then Goto C:Else Goto 3:IfEng↙LbI 3↙If L≦155620.15:Then Goto D:Else Goto 4:IfEng↙LbI 4↙If L≦155790.15:Then Goto E:Else Goto 5:IfEng↙LbI 5↙If L≦156656.233:Then Goto F:Else Goto 6:IfEng↙LbI 6↙If L≦156806.233:Then Goto G:Else Goto 7:IfEng↙LbI 7↙If L≦157432.399:Then Goto H:Else Goto 8:IfEng↙LbI 8↙If L≦157607.399:Then Goto I:Else Goto 9:IfEng↙LbI 9↙If L≦158977.376:Then Goto J:Else Goto 0:IfEng↙LbI 0↙If L≦160298.935:Then Goto K↙LbI A↙58874.9732→A:534215.4827→B:1608.999→R:180→S:154426.087→O:183°37′56.5″→Z:-1→K:-1→W↙LbI B↙588974.9732→A:534215.4827→B:1000→R:150→S:154426.087→O:3°37′56.5″→Z:-1→K:-1→W↙Goto T↙LbI C↙589760.5335→A:533988.4355→B:1000→R:150→S:155265.597→O:144°07′34. 8″→Z: 1→K:-1→W↙Goto T↙LbI D↙589760.5335→A:533988.4355→B:1522.495→R:170→S:155265.597→O:324°07′34.8″→Z: 1→K:1→W↙Goto T↙LbI E↙590216.856→A:533733.3811→B:1522.495→R:170→S:155790.15→O:157°28′0 8.8″→Z: -1→K:-1→W↙Goto T↙LbI F↙590216.856→A:533733.3811→B:1028.84→R:150→S:155790.15→O:337°28′08 .8″→Z: -1→K:1→W↙Goto T↙LbI G↙590888.4984→A:533021.9777→B:1028.84→R:150→S:156806.233→O:109°14′14.1″→Z: 1→K:-1→W↙Goto T↙LbI H↙590888.4984→A:533021.9777→B:1545.627→R:175→S:156806.233→O:289°14′14.1″→Z: 1→K:1→W↙Goto T↙LbI I↙591295.3826→A:532340.5881→B:1545.627→R:175→S:157607.399→O:132°26′56.2″→Z: -1→K:-1→W↙Goto T↙LbI J↙59129.3826→A:532340.5881→B:0→R:0→S:157607.399→O:312°26′56.2″→V: 1→K:1→W↙Goto R↙LbI K↙592220.0262→A:531329.7104→B:2097.935→R:350→S:158977.376→O:312°26′56.2″→Z: -1→K:1→W↙Goto T↙LbI R↙“S=”:A+(L-O)×cos(V)→S◢“T=”:B+(L-O)×sin(V)→T◢Goto W↙LbI T↙Abs(L-O)→H↙If H-S≧0:Then Goto U:Else Goto V↙LbI U↙Prog“YQX”↙S=S:T=T:V=V↙Goto W↙LbI V↙H-H^(5)÷40÷R2÷S2→M: H^(3)÷6÷R÷S-H^(7)÷336÷R^(3)÷S^(3)→N↙K×(H2×180÷2÷∏÷R÷S)+Z→V:“S=”:A+M×cos(Z)-N×K×sin(Z)→S◢“T=”:B+M×sin(Z)+N×K×cos(Z)→T◢Goto W↙LbI W↙8↙ZZB （表示什么？）If L≦1156653.357:Then Goto A:Else Goto 1:Ifend↙LbI 1↙If L≦1156803.357:Then Goto B:Else Goto 2:Ifend↙LbI 2↙If L≦1157450.193:Then Goto C:Else Goto 3:Ifend↙LbI 3↙If L≦1157660.193:Then Goto D:Else Goto 4:Ifend↙LbI 4↙If L≦1158956.336:Then Goto E:Else Goto 5:Ifend↙LbI 5↙If L≦1160306.154:Then Goto F:Ifend↙LbI A↙590216.856→A:533733.3811→B:1000.32→R:150→S:1155790.15→O:337°28′0 8.8″→Z: -1→K:1→W↙Goto T↙LbI B↙590877.7692→A:533018.36→B:1000.32→R:150→S:1156803.357→O:108°01′3 6.5″→Z: 1→K:-1→W↙Goto T↙LbI C↙590877.7692→A:533018.36→B:1517.525→R:210→S:1156803.357→O:288°01′36.5″→Z: 1→K:1→W↙Goto T↙LbI D↙591304.6519→A:532286.014→B:1517.525→R:210→S:115760.193→O:132°26′55.7″→Z: -1→K:-1→W↙Goto TLbI E↙591304.6519→A:532286.014→B:0→R:0→S:115760.193→O:312°26′55.7″→V: 1→K: 1→W↙Goto TLbI F↙592179.4593→A:531329.6152→B:2240.219→R:350→S:1158956.336→O:312°26′55.7″→Z: -1→K: 1→W↙Goto TLbI R↙“S=”:A+(L-O)×cos(V)→S◢“T=”:B+(L-O)×sin(V)→T◢Goto W↙LbI T↙Abs(L-O)→H↙If H-S≧0:Then Goto U:Else Goto V↙LbI U↙Prog“YQX”↙S=S:T=T:V=V↙Goto W↙LbI V↙H-H^(5)÷40÷R2÷S2→M: H^(3)÷6÷R÷S-H^(7)÷336÷R^(3)÷S^(3)→N↙K×(H2×180÷2÷∏÷R÷S)+Z→V:“S=”:A+M×cos(Z)- N×K×sin(Z)→S◢“T=”:B+M×sin(Z)+N×K×cos(Z)→T◢Goto W↙LbI W↙8↙（完）……。

公路路线座标正反算（积分公式）通用程序公路路线座标正反算由于现在计算机普及，计算机功能日益强大，宜采用较简单的积分公式，便于计算机处理。

单线元通用积分公式如下M = (1.0/Re-1.0/Rs)/Ls;x=∫｛cos(Ta + L/Rs + 0.5*M *L*L)，0，L｝; y=∫｛sin(Ta +L/Rs + 0.5*M *L*L)，0，L｝; a(i)= Ta +L/Rs + 0.5*M *L*LRs:缓和曲线起点半径Re:缓和曲线止点半径Rs，Re （NE坐标系下，右偏为正，左偏为负）Ta:缓和曲线起点的真北方位角Ls:不完整缓和曲线长度。

此公式为缓和曲线在坐标系下任意位置的通用积分公式，能完全适应缓和曲线左偏、右偏、Rs >Re 、Rs <Re 等各种情况，不必先凑成完整缓和曲线，降低算法的复杂程度。

虽然此公式是由缓和曲线推导出来，也可和于直线与圆曲线，可降低计算机编程的复杂程度。

Fx-5800计算机程序QXJS-000主程序Lbl 4:“1.SZ=>NE”:“2.NE=>SZ”:?Q:Z[12]→S：?S：Prog “QXJS-SUB0”↙Lbl 0:Q=1 => Goto1:Q=2 => Goto2:↙Lbl 1:?Z:Prog “QXJS-SUB1”:“N=”:N◢:“E=”:E ◢:“F=”:F DMS◢: Goto3↙Lbl 2: “N=”:?B: “E=”:?C:B→N: C→E:Prog“QXJS-SUB2”:“S=”:S◢: “Z=”:Z ◢: Goto4↙Lbl 3：J＜0＝＞J＝J+360：“FW＝”：J DMS◢Goto4↙S→Z[12]:Z→Z[13] Prong “SD”QXJS-SUB0数据库子程序IF S<***（线元终点里程）:Then***→A（线元起点方位角）:***→O（线元起点里程）:***→U（线元起点X）:***→V（线元起点Y）:***→P（线元起点曲率半径）:***→R（线元终点曲率半径）: ***→L （线元起点至终点长度）: Return:IfEnd↙IF S<35693.064:Then257°13′1″→A:35643.064→O: 3140794.382→U: 514646.536→V:10^(45)→P:965→R:50→L: Return:IfEnd↙QXJS-SUB1正算子程序0.5（1÷R-1÷P）÷L →D:S-O→X↙U+∫(cos(A+(X÷P+DX2)×180÷π,0,X)→N↙V+∫(sin(A+(X÷P+DX2)×180÷π,0,X)→E↙A+(X÷P+DX2)×180÷π→F↙N+Zcos(F+90) →N:E+Zsin(F+90) →E QXJS-SUB2反算子程序Lbl 1:0→Z：1→Q：Prog “QXJS-SUB0”: Prog “QXJS-SUB1”↙Pol(N-B+10^(-46),E-C+10^(-46)):Isin(F-9 0-J) →W:S+W→S↙Abs(W)>0.0001 => Goto1↙Lbl 2: 0→Z：Prog “QXJS-SUB1”:(C-E) ÷sin(F+90) →Z隧道超欠挖程序程序名：SDProg“SD1”: 0.84→A:”CQHD”?A:A→Z[15]: 0→M:”Hc1=”?M:M→Z [8]:Prog“SQX”:Cos(Z [10]) ×Z[4] →E: H +Z[5] →Z[5]: H+Z[7]→Z[7]: E+Z[5] →Z [10]:If Z[8]>Z[10]: Then Goto R:Else Goto S:IfEnd:LbI R:√((Z[13]-Z [9]) 2+(Z[8]-Z[5]) 2 ) - Z[4]-Z[15] →L:Go to L:LbI S:Z[4]-Z[6] →S:√(S 2-(Z[7]-Z[5]) 2 ) ×S÷Abs(S)→T:A bs(Z[13]-Z[9])-T→T:√(T 2+(Z[8]-Z[7]) 2 )- Z[6]-Z[15] →L:Go to L:LbI L: Fix 3:＂CQW=＂: Locate 6, 4,L：＂H=＂: Locate 6,4, Z[8]-H＂Ｋ=＂：Locate 6,4,Ｚ[12]：＂Z=＂：Locate 6,4,Z [13]：◢竖曲线计算程序程序名：SQXZ[12]→KProng “S1”0.5R Abs(I-J)→T If I-J﹤0: Then -1→V:Else 1→V：IfEnd K-Z→B:If K≤Z-T:Then BI+A →H:Return:IfEnd:If K≤Z:Then A+BI-V (T+K-Z)2÷（2R）→H: Return:IfEndIf K≤Z+T:Then A+BJ -V(T-K+Z)2÷（2R）→H: Return:IfEndIf K﹥Z+T:Then A+BJ →H:Return:IfEnd程序名：S1(线路1高程竖曲线要素子程序)LbI 1:If Z≤26157.8: Then 200000→R:0.00 3→I:-0.0017→J:705 18.232→Z:380.342→A:Goto 0:IfEnd:If Z≤27421.915:The n 300000→R: 0.003→I:-0.0017→J:7051 8.232→Z:380.342→A: Goto 0:IfEnd:If Z≤27889.029:The n 1000000→R: 0.003→I:-0.0017→J:7051 8.232→Z:380.342→A: Goto 0:IfEnd:LbI 0程序字母说明：R-竖曲线半径；I－竖曲线前纵坡; J－竖曲线后纵坡；Z-竖曲线交点桩号；A-交点桩号高程（未竖曲线调整的）程序名：SD1(线路1隧道参数子程序)LbI 1:5.53→Z[4]:1. 62→Z[5]:8→Z[6]:1. 62→Z[7]:-5.375→Z [9]: 90→Z[10]:LbI 0子程序中字母表示说明：Z[4]－隧道上拱半径；Z[5]－隧道上拱圆心至设计面高度；Z[6]－隧道中腰部半径；Z[7]－中腰圆心至设计面高度；Z[9]－设计隧道中线与设计路基中线偏移值，左为负数，右为正数；Z[10]－隧道上部拱跨半幅角度三、使用说明1、规定(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当曲线半径在左时，P、R取负值，当曲线半径在右时，P、R 取正值，当曲线半径为无穷大（即直线）时，P、R以10的45次代替。

隧道超欠挖5800版1. 主程序(NTSSD--CQW)Lbl 1：5→DimZ：（增加Z型5个变量）600→S：（输入近似里程）Prog“SUB0”↙(定义所在线元数据库)1÷P→C: (P-R)÷(2HPR) →D:180÷∏→E:“X=”？X：“Y=”？Y：X→I：Y→J：(输入实测XY坐标,)“H=”？→Z[1]: (输入实测高程)“R=”？→Z[2]: (输入断面理论半径)Prog“SUB2”：(调用反算子程序,由反算两次调用正算子程序)O+W→S：“K=”：S◢“BIANJULI=”：Z◢(得到桩号及离中心桩距离) Abs(Z)→B： (将变量值另存,利于后面计算)617.42→C： (定义此段竖曲线变坡点桩号，也可？C：输入)71.1→D： (定义此段竖曲线此段竖曲线变坡点高程)1→I： (定义此段竖曲线后坡点坡度)-0.6→J： (定义此段竖曲线前坡点坡度)1500→R： (定义此段竖曲线变坡半径值)-1→Q： (判断是凸凹曲线类型, -1→Q：I<J=>1→Q )240→L： (定义此段竖曲线曲线长)120→T： (定义此段竖曲线切线长)0.48→E： (定义此段竖曲线外距)497.42→M： (定义此段竖曲线变坡曲线起点桩号)737.42→N： (定义此段竖曲线变坡曲线终点桩号)S→K： (将变量值另存,利于后面计算)Prog“SUB3”： (调用高程计算子函数)Fix 3：（输出保留3位小数）1.5+Z→Z： (将路面设计值加减设计值得到开挖圆心高程)开平方根((Z[1]-Z)2+B2)→X：(得到实测点到圆心的距离)“CQW=”:X-Z[2]→A：◢(将实测点半径减理论半径即为超欠挖)ClrMemory:Goto 12. SUB0 数据库子程序Goto1↙ 同时保存多个曲线时的指针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:IfEnd……………..为了便于解读，每增加一个线元增加一行语句，每增加一条曲线增加一个Lbl，每增加一个工程增加一个文件3. 正算子程序(SUB1)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+Zcos（F）→X： Y+Zsin（F）→Y4. 反算子程序(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)cos（L）-(I-X)sin（L）→Z：IF Abs（Z）<0.000001：Then0→Z:Prog "SUB1" :(J-Y)÷sin（F）→Z↙Else W+Z→W:Goto4:IfEnd5高程计算子程序(SUB3)IF K≤ M: Then D+(K-C)I/100→Z: Return:IfEnd↙IF K> M And K≤C :Then (K-M)2/(2R) →P:D+(K-C)I/100+QP→Z:Return:IfEnd↙IF K> C And K≤N :Then (K-M)2/(2R) →P:D+(K-C)J/100+QP→Z:Return:IfEnd↙IF K> N: Then D+(K-C)J/100→Z: Return:IfEnd↙经反复核对的，绝对能用，5800P上运行要4秒出结果。

三心圆隧道超欠挖,不足之处望高手指点二、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 输出超欠数据结果。

1. 主程序(TYQXJS)Lbl 4："1SZ => XY，2XY => SZ"？N：？S：Prog“SUB0”↙1÷P→C: (P-R)÷(2HPR) →D:180÷∏→E:N=1 => Goto1：Goto2：↙Lbl 1:”DZ”? Z：Abs(S-O) →W：Prog "SUB1"："XS"：X◢"YS"：Y◢Prog “SQX”（竖曲线的程序名字）F-90→F：“FS=”：F▲DMS◢Prog “CQW”（隧道超欠挖的程序名字）Goto4↙Lbl 2：”X”？X：“Y”？Y：X→I：Y→J：Prog“SUB2”：O+W→S：“S”：S◢“Z”：Z◢Prog“CQW”Goto4↙2. 正算子程序(SUB1)0.1739274226→A：0.3260725774→B：0.0694318442→K：0.3300094782→L：1-L→F：1-K→MG+QEKW(C+KWD)→Z[1]G+QELW(C+LWD)→Z[2]G+QEFW(C+FWD)→Z[3]G+QEMW(C+MWD)→Z[4]A×cos(Z[1])→XX+Bcos(Z[2])→XX+Bcos(Z[3])→XX+Acos(Z[4])→XU+WX→XAsin(Z[1])→YY+Bsin(Z[2])→YY+Bsin(Z[3])→YY+Asin(Z[4])→YV+WY→YG+QEW(C+WD)+90→FX+Zcos(F)→XY+Zsin(F)→Y2. 反算子程序(SUB2)G-90→T：(Y-V)cos(T)-(X-U)sin(T) →W：Abs(W) →W:0→ZLbl 0：Prog "SUB1"T+QEW(C+W D) →L：(J-Y)cos（L）-(I-X)sin（L）→Z：IF Abs（Z）<0.000001：Then0→Z:Prog "SUB1"(J-Y)÷sin（F）→ZElse W+Z→WGoto0:IfEndSUB0 数据库子程序Goto1↙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:IfEnd……………..为了便于解读，每增加一个线元增加一行语句，每增加一条曲线增加一个Lbl，每增加一个工程增加一个文件。

正算主程序(ZS)FIX4：?S：?Z：Prog “PM-SJ”： Abs(S-O)→ W：Prog "SUB1"： "XS ="：X◢"YS="：Y◢F-90→F(需要时可以让他显示)：Prog“SQX-SJ”: Prog“SQX”：“H=”：H◢反算主程序(FS)FIX4: ?S：?X：？Y：Prog“PM-SJ”： X→ I： Y→ J：Prog "SUB2"： "S="：O+W→S◢"Z="：Z◢Prog“SQX-SJ”：Prog“SQX”：“H=”：H◢边坡放样主程序（BP-FY）Lb1 0：Prog“ZS”：“H-BG”（中桩与坡脚起算点高差值，比中桩高正，反之负）?A：H+A →B：?P（实测点高程）：?L（坡脚起算点到中桩的距离）：0.75（挖方时一级坡度）→C：1（挖方时二级坡度）→D：1.5（挖方时三级坡度）→E：8（挖方时一级坡高）→G：10（挖方时二级坡高）→M：15（挖方时三级坡高）→N：2（平台宽度）→K：1（填方时一级坡度）→I：1.5（填方时二级坡度）→J：2（填方时三级坡度）→O：2（填方时一级坡高）→Q：8（填方时二级坡高）→R：10（填方时三级坡高）→T：ifP＞B：thenGoto1：ElseGoto2Lb1 1：ifP＞B：thenL+C（P-B）→U：P-B→F：IfEnd←┘ifP＞B+G：th enL+CG+K+D（P-B-G）→U：P-B-G→F：IfEnd←┘ifP＞B+G+M：thenL+GC +2k+MD+E（P-B-G-M）→U：P-B–G-M→F：IfEnd：Goto3←┘Lb1 2：ifP≤B：thenL+I（B -P）→U： B -P→F：IfEnd←┘ifP≤B-Q：thenL+IQ+K+J（B-Q-P）→U：B-Q-P→F：IfEnd←┘ifP≤B-Q-R：thenL+I Q+2k+JR+O（B-Q-R-P）→U：B-Q-R-P→F：IfEnd：Goto3←┘Lb1 3：U-AbsZ→V：ifZ＜0：thenZ-V→Z：Else z+v→z：IfEnd←┘“Z=”：Z◢计算得出正确的宽度，路线左为负，右为正。