隧道超欠挖程序

If CElse √((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]=”:-I cos(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 1 If O≠1:Then “L”? →Q:”PJ”? →J:IfEnd：如选择边线提示输入宽度及偏角，宽度始终输正值，偏角右法线输+90，左法线输-90 Prog “ZB”

X+Rec(Q,U+J) →X Y+J→Y “X=”:X◢ “Y=”:Y◢

If O=1:Then “QXFWJ=”:U▲DMS◢IfEnd：当计算中线时，会显示该点的切线方位角 Goto 1

11、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左转）

If Abs(C) >180:Then 360-Abs(C) →C :If END Abs(C) →C

L÷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]→T Z[7]-P+Z[17]→P

πRC÷180+(L+M) ÷2→Z[7]

注：其余公式请参照不等长缓和曲线计算公式

N=1右转，（左右转参数计算，

12、ZB（坐标计算子程序） Lb1 D 22→W:?K Lb1 2

If KZ[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] →S A→B

If S≤0:Then Z[14]+Scos(B) →X：（第一直线段坐标计算） Z[15]+Ssin(B) →Y B→U:Goto 7:IfEnd Z[14]→X:Z[15]→Y

If S≤L:Then 1→H:Goto 4:IfEnd：（第一缓和曲线段坐标计算） If S≤Z[7]-M:Then 1→H:Goto 5:IfEnd：（圆曲线段坐标计算） A+NC→B

If 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) →I S-S^(5)÷40÷R2÷L2→J

B+90NHS2÷π÷R÷L→U:Goto 6 Lb1 5（圆曲线计算公式） 90(2S-L) ÷π÷R→U R(1-cos(U))+Z[18]→I Rsin(U)+Z[16]→J

B+180N(S-L÷2)÷π÷R→U Lb1 6

Pol(J,I):NHJ+B+90(H-1)→J Rec(I,J)+X→X Y+J→Y Lb1 7

If U<0:Then U+360°→U:IfEnd

13、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 0

14、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 0

15、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)→X Y+J→Y

Pol(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 15

d：偏移两圆心（O2、O3）与隧道中线之水平距离 ？：线路中线与隧道中线之偏移距离

Z[8]：左右线判断变量（分离式双线隧道）及隧道左右偏判断变量（相对线路设计中心线）

平曲线及竖曲线数据库设置表（n个平交点）

1、平曲线数据库设置表

平曲线起第一个交第二个交„„„„ 第i个交点 „„„„ 第n-1个交第n个交点 平曲线终点 点 点 点 点 X坐标 Z[20] Z[21] Z[22] Z[20+i] Z[19+n] Z[20+n] Z[21+n] Y坐标 Z[22+n] Z[23+n] Z[24+n] Z[22+n+i] Z[21+2n] Z[22+2n] Z[23+2n] 平曲线半 Z[24+2n] Z[25+2n] Z[23+2n+i] Z[22+3n] Z[23+3n] 径 Ls1 Z[24+3n] Z[25+3n] Z[23+3n+i] Z[22+4n] Z[23+4n] Ls2 Z[24+4n] Z[25+4n] Z[23+4n+i] Z[22+5n] Z[23+5n] 曲线起点Z[24+5n] Z[25+5n] Z[26+5n] Z[24+5n+i] Z[23+6n] Z[24+6n] Z[25+6n] 里程 注：1、从Z[20]开始存储数据，是为了留出前面的扩展变量作为其他变量使用，以补充变量A-Z之不足，当然也可以从Z[30]、 Z[40]开始，留出更多的变量供使用，只需根据上面表格作相应调整。

2、曲线起点里程，对每一个交点来说，指的是对应的直缓点ZH（直园点ZY）里程。 3、Ls1—第一缓和曲线长度（没有时输0）；Ls2—第二缓和曲线长度（没有时输0）。

2、竖曲线数据库设置表（g个变坡点）

竖曲线起点 第一个变坡点 第二个变坡点 „„„„ 第i个变坡点 „„„„ 第m-1个变坡点 第m个变坡点 竖曲线终点

变坡点里程 Z[26+6n] Z[29+6n] Z[32+6n] Z[26+6n+3i] Z[23+6n+3g] Z[26+6n+3g] Z[29+6n+3g] 变坡点高程 Z[27+6n] Z[30+6n] Z[33+6n] Z[27+6n+3i] Z[24+6n+3g] Z[27+6n+3g] Z[30+6n+3g] 竖曲线半径 Z[28+6n]为空 Z[31+6n] Z[34+6n] Z[28+6n+3i] Z[25+6n+3g] Z[28+6n+3g] 注：在设置扩展变量时，扩展变量的数目不得小于30+6n+3g