自编CASIO 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)÷(2NPR) →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←┘” 中改为“If H>内轨顶标高+圆心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)等曲线要素赋予定值，可减少在工作中的输入量。

二，说明：1. 程序中内轨顶标高用包含S（里程）的式子表示。

          2. 程序中线中相对遂中距离采取左负右正

          3. 圆心相对遂中距离，若圆心与所对应的弧在遂中同一侧则取负号，若

在同一侧取正号。

三，算例：      例如本人现在所在隧道为兰渝铁路新城子隧道出口，目前正在施工右线大断面与左线小断面。

右线大断面所在的线元段为直线，其HZ点X坐标为3762283.97，Y坐标为483022.06，里程为2728.016，方位角为153-27-10.26，线元长度为1579.479，并且该段是以起点里程276701（内轨顶高程为1734.4526）终点里程277174（内轨顶高程为1731.6146）坡度为6‰的下坡。

该段标准断面的圆心及遂中，线中关系如下图：

那么该段的超欠挖程序主程序如下：

主程序（CQW-YOUDA）

3762283.97 →U：483022.06 →V：2728.016 →O：153°27°10.26°→G：1579.479 →N：×1045 →P：×1045 →R：0 →Q←┘

1÷P→C：(P-R)÷(2HPR) →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+0.006S-3396.909) ²+(Abs(Z-2.5)) ² )-7.36→M←┘

"CQW="：M◢

Goto 2

当输入

       X=3761963.219 Y=483171.325  H=1737.714   算得 S=276781.661  CQW=+0.118

 X=3761966.574 Y=483177.777  H=1743.5  算得 S=276781.543   CQW=+0.066

      X=3761969.783  Y=483184.431  H=1737.78  算得S=276781.6  CQW=+0.135

又如：

左线小断面所在的线元段为缓和曲线，其YH点X坐标为3762418.968，Y坐标为482984.424，里程为276293.087，方位角为151-32-35.26，线元长度为300，圆曲线半径为4500，线路为右偏。并且该段是以起点里程275390（内轨顶高程为1734.4526）终点里程276599（内轨顶高程为1735.4114）坡度为12.8‰的下坡。

该段标准断面的圆心，遂中，线中及圆心所对应弧的分段高度关系如下图：

那么该段的超欠挖程序主程序如下：

主程序（CQW-ZUOXIAO）

3762418.968→U：482984.424→V：276293.087→O：151°32°35.26°→G：300→N：4500→P：×1045 →R：1→Q←┘

1÷P→C：(P-R)÷(2HPR) →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◢

If H>5282.7-0.0128S：Then√( (H+0.0128S-5280.86) ²+( Z+0.15) ²)→M：

Else√( (H+0.0128S-5278.6768) ²+(Abs ( Z+0.15)+3.42) ²)-8.06→M：If End←┘

"CQW="：M◢

Goto 2

当输入

X=3762247.585 Y=483079.072 H=1740.795 算得 S=2788.755  CQW=+0.174

X=3762247.163 Y=483078.211 H=1743.471 算得S=2788.744  CQW=+0.097

X=3762245.235 Y=483074.534 H=1745.735 算得S=2788.808  CQW=+0.051

X=3762243.659 Y=483071.107 H=1742.710 算得S=2788.671  CQW=+0.116

X=3762243.347 Y=483070.514 H=1739.479 算得S=2788.683  CQW=+0.181