隧道开挖程序

1、线路坐标计算里程I=0：J=0：Lb1 0:xyz:LMXYZ:Pol(9162.8184-X,73723.7612-Y):J=J+180:N=7564.628+7060X(J-308-47-47.9)/180:E=7060-I+6.07Goto 1:Lbi 1:Pol(13000.466-X,67750.3949-Y):J=J+180:Rec(I,J-38-47-47.9):N7564.628N=6812.934+I:E=J+6.07Goto 2LBI 2:”N=”:N“E=”:EGoto 3LBI 3: N8104.721H=99913.00052-(1000002-(N-10968.68

2、72)2): N6540H=-1.224-(N-6540)*0.0286H=H-G:“H=”:HPol(H+1.531-Z,1.597-Abs E):J=J+180:J180A=(H+0.95-Z)2+E2)-(7.4+M):J70A=I-5.7-M: A=(H+0.95-Z)2+E2)-(7.4+M)L=1“R=”：AGoto 4LBI 4：E4.87V=H-0.48+E*0.02-Z: L=2“V=”：VGoto 0程序中字母表示的意思：L 输入2 程序运行一段路面，输入其他任意值计算开挖M 程序中是根据二衬净空编制的程序M输入0计算二衬，比如衬砌厚度是30公分，M输入0.3计算开挖XYZ

3、 是大地坐标N 计算出来的里程E 是到隧道中线的距离H是某里程上的设计标高“R=”是最终结果正超、负欠程序中第一段为圆曲线坐标反算， 第二段为直线的坐标反算 第三段是 线路高程计算和隧道开挖二衬超前挖计算，第四段是隧道路基计算，原本还有第五段仰拱开挖，被小子删掉了 程序中的数据是本工地上的数据，第一段中9162.8184,73723.7612，为圆曲线中的圆心坐标，7564.628为圆曲线起点里程，308-47-47.9为圆心至ZY点的方位角， 第二段中13000.466 67750.3949为直线起点的坐标，6812.934直线起点里程，38-47-47.9线路前进方位角第三段的中的高程计算

4、中第一段是竖曲线计算，大家很可能不能理解，给大家解释一下我的计算方式，第一步在CAD中根据里程为X坐标，高程为Y坐标，把竖曲线绘制出来，标注出圆心的里程和高程 ，我的程序中99913.00052就是这段竖曲线的圆心高程，1W是半径，10968.6872是竖曲线圆心的里程，程序中的N表示计算点的里程，更具勾股定理计算出半径对应的一条直角边长减去圆心高程就得出计算点的设计高程，第四段是隧道超欠挖计算，计算隧道超欠有很多种方式计算，以上程序是用角度来判定以上衬砌断面是3心圆组成，选其中任意的一个圆心创建一个直角坐标系， 上诉程序中用5.7半径圆创建的直角坐标系，POL(设计标高+第二圆心的高度-实测

5、高程，隧道中线到第二圆心的距离-绝对值实测的偏距（E）)而得到一个直角坐标系，这个坐标系主要的目的是计算衬砌中半径角度（象限角） ，隧道中线偏距为0，正北方向角为0度。现在的隧道均为对称形，左右完全一样，所以判定隧道超欠挖只需要判断一边即可判定的角度计算方式，如（以上5.7半径的圆正北方向角度为70度，而不是70°加28°）不管隧道中有几个半径圆，都是一样的道理，卡西欧计算器中角度判断必须由于最大角度开始判定，所以程序中第一步必须大于等于180°推出7.4半径的圆（J180A=(H+0.95-Z)2+E2)-(7.4+M)），其次大于等于70°（也可以用

6、勾股定理，而以上直角坐标系中，POL计算中默认变量I的结果就是实测半径J70A=I-5.7-M，最后大于等于0°（则返回到 7.4半径圆中A=(H+0.95-Z)2+E2)-(7.4+M)）。从而循环。上诉的角度判定不管开挖实测坐标点在任何位置，计算器都能计算出象限角，从而判定出在那一个半径圆心上，用勾股定理，计算实测半径减去设计半径，得到的结果正为超挖，负为欠挖。不管隧道中是几个圆心，同样的道理，例如5.7半径下面还有一个小的圆，在角度判定中最好用下面的小圆创建直角坐标系（你选择第二个圆心创建直角坐标系也可以），同样的道理，大于等于180°开始推出7.4半径的圆，大于等于那个小圆的角度(例如105°）推出小圆，其次才是大于等于70°推出5.7半径的圆，最后大于等于0°返回了7.4半径圆以上程序中没有表示大于等于0&#