5800正反算、边坡放样、隧道超欠挖

正算主程序(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：ElseGoto2 Lb1 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：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+IQ+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◢ 计算得出正确的宽度，路线左为负，右为正。 “V=”V◢ 宽度调整值，正为远离中线，负为向中线靠近。 “PT=”：F◢ 距离即将到达的填挖方平台的高度。 注：由于路基的宽度和超高情况较为复杂，因此该程序的缺点是不能自动计算路基宽度和路基边缘坡脚起算点的设计高程，需要手算坡脚起算点与中桩的宽度及与中桩的高差。 算例；S? 700（放样点桩号） Z？-15（试算宽度） X=19839.5817 Y=28515.5010 H=107.144 H-BG？ -0.15（边缘起算点的高程比中桩低0.15） p?150（实测高程）L? 12（起算边缘点距中桩的距离） Z=-69.509（计算得出正确的宽度）V=54.509（远离中线方向距测点54.509米），PT=25.006（距离平台高度）。 然后再次进入下一轮试算，直到V趋近于零。 隧道3心圆放样主程序（SD-3XY） Lb1 0：Prog“FS”：1.5→A：1.65→B：0.89→G：2.27→C：3.3→D：2.41→E：0.51→J：“H1”：F：ifF≤H+A：then E-AbsZ→W：IfEnd：ifF＞H+A And F≤H+A+B：then√（D²-（F-H-A）²）-G-AbsZ→W：IfEnd：ifF＞H+A+B：then√（C²-（F-H-A-J）²）-AbsZ→W：IfEnd：“W=”：W◢ （水平方向的超欠挖：正为欠,负为超）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)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：IfEnd←┘ Lbl 1： 0→ Z：Prog "SUB1"： (J-Y)÷sinF→ Z 子程序（平面数据库）PM-SJ ifS ≥500（线元起点里程）AndS＜769.256（线元止点里程）：Then19942.837 →U（线元起点X坐标）：28343.561 →V（线元起点Y坐标）：500 →O（线元起点里程）：125°16°31°→G（线元起点方位角）：269.265 →H（线元长度）：10^（45）→P（线元起点曲率半径）：10^（45）→R（线元终点曲率半径）：0 →Q（线元左右偏标志：左负右正）：IfEnd←┘ ifS ≥769.256（线元起点里程）AndS＜806.748（线元止点里程）：Then19787.34 →U（线元起点X坐标）：28563.378→V（线元起点Y坐标）：769.256→O（线元起点里程）：125°16°31°→G（线元起点方位角）：37.492→H（线元长度）：10^（45）→P（线元起点曲率半径）：221.75→R（线元终点曲率半径）：-1→Q（线元左右偏标志：左负右正）：IfEnd←┘ ifS ≥806.748（线元起点里程）AndS＜919.527（线元止点里程）：Then19766.566→U（线元起点X坐标）：28584.574→V（线元起点Y坐标）：806.748→O（线元起点里程）：120°25°54.07°→G（线元起点方位角）：112.779→H（线元长度）：221.75→P（线元起点曲率半径）：221.75→R（线元终点曲率半径）：-1→Q（线元左右偏标志：左负右正）：IfEnd←┘ ifS ≥919.527（线元起点里程）AndS＜999.812（线元止点里程）：Then19736.072 →U（线元起点X坐标）：28701.893 →V（线元起点Y坐标）：919.527→O（线元起点里程）：91°17°30.63°→G（线元起点方位角）：80.285→H（线元长度）：221.75→P（线元起点曲率半径）：9579.228→R（线元终点曲率半径）：-1→Q（线元左右偏标志：左负右正）：IfEnd←┘ ifS ≥999.812（线元起点里程）：Then19744.038 →U（线元起点X坐标）：28781.659→V（线元起点Y坐标）：999.812→O（线元起点里程）：80°40°50°→G（线元起点方位角）：100→H（线元长度）：10^（45）→P（线元起点曲率半径）：10^（45）→R（线元终点曲率半径）：0 →Q（线元左右偏标志：左负右正）：IfEnd←┘ 子程序（竖曲线数据库）SQX-SJ ifS≤999.812（竖曲线终点里程）：then0.0357（前坡度：上坡正下坡负）→E：0.0227（后坡度：上坡正下坡负）：600→R（竖曲线半径）：780→G（变坡点里程）：110→C（变坡点高程）：IfEnd←┘ 子程序（竖曲线计算公式）SQX E-F→J：ifJ＞0：then-R→R：ElssR→R：Abs（R*J÷2）→T：ifS≤G-T：thenG-S→L：C-LE→H：IfEnd←┘ ifS＞G-T And S≤G： thenG-S→L：（S+T-G）²÷2÷R→N：C-LE+N → H：IfEnd←┘ ifS＞G And S≤G+T： then S-G→L：（G +T- S）²÷2÷R→N：C-LF+N → H：IfEnd←┘ ifS＞G+T ：thenS-G→L： C+LF → H：IfEnd←┘ 三、使用说明 1、规定         (1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时， Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。         (2) 当所求点位于中线时，Z=0；当位于中线左铡时，Z取负值；当位于中线中线右 侧时，Z取正值。         (3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。         (4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆 弧的半径。          (5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45 次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半 径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。          (6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的 值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等 于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。 2、输入与显示说明       S ？ 正算时所求点的里程，反算时估算点里程，反算时尽量不要超出真实测点所在的线元里程终点，即尽量向小里程估计，如超出线元终点计算器将死机。       Z ？正算时所求点距中线的边距(左侧取负，值右侧取正值，在中线上取零)       X ？反算时所求点的X坐标       Y ？反算时所求点的Y坐标        显示部分：        XS=×××    正算时，计算得出的所求点的X坐标        YS=×××    正算时，计算得出的所求点的Y坐标        S=×××    反算时，计算得出的所求点的里程        Z=×××    反算时，计算得出的所求点的边距 H=××× 正反算时，计算得出所求点的中桩设计高程（如弯道沿边线超高时要另外考虑超高段中桩高程） 四、算例         某匝道的由五段线元（直线+完整缓和曲线+圆曲线+非完整缓和曲线+直线）组 成，各段线元的要素（起点里程S0、起点坐标X0 Y0、起点切线方位角F0、线元长度 LS、起点曲率半径R0、止点曲率半径RN、线 元左右偏标志Q）如下： S0   X0       Y0                   F0            LS R0 RN     Q 500.000    19942.837    28343.561     125 16 31.00    269.256     1E45 1E45      0 769.256    19787.340    28563.378     125 16 31.00      37.492     1E45 221.75     -1 806.748    19766.566    28594.574     120 25 54.07     112.779    221.75 221.75     -1 919.527    19736.072    28701.893       91 17 30.63       80.285    221.75 9579.228   -1 999.812    19744.038    28781.659       80 40 50.00     100.000    1E45 1E45       0          1、正算           （注意：略去计算方式及线元要素输入，请自行根据所求点所在的线元输入线元 要素）             S=700  Z=-5    计算得  XS=19831.41785   YS=28509.72590 H=107.144             S=700  Z=0     计算得  XS=19827.33592   YS=28506.83837 H=107.144             S=700  Z= 5    计算得  XS=19823.25398   YS=28503.95084 H=107.144             S=780  Z=-5    计算得  XS=19785.25749   YS=28575.02270 H=112.759             S=780  Z=0     计算得  XS=19781.15561   YS=28572.16358 H=112.759             S=780  Z= 5    计算得  XS=19777.05373   YS=28569.30446 H=112.759                         S=870  Z=-5    计算得  XS=19747.53609   YS=28654.13091 H=112.759             S=870  Z=0     计算得  XS=19742.68648   YS=28652.91379 H=112.759             S=870  Z= 5    计算得  XS=19737.83688   YS=28651.69668 H=112.759                     S=940  Z=-5.123    计算得  XS=19741. 59118  YS=28722.05802 H=130.475             S=940  Z=0            计算得  XS=19736.47687   YS=28722.35642 H=130.475             S=940  Z= 3.009    计算得  XS=19733.47298   YS=28722.53168 H=130.475                     2、 反算              X=19831.418   Y=28509.726   计算得  S=699.9999974  Z= -5 .00018164              X=19827.336   Y=28506.838   计算得  S=699.9996493  Z= 0.000145136 H=107.144              X=19823.25398   Y=28503.95084    计算得  S=699.9999985  Z= 5.000003137            X=19785.25749   Y=28575.02270  计算得  S=780.0000035  Z= -5 .000001663            X=19781.15561   Y=28572.16358   计算得  S=780.0000025  Z=- 0.000002979H=112.759            X=19777.05373   Y=28569.30446   计算得  S=780.0000016  Z= 4.99999578                       X=19747.536   Y=28654.131   计算得  S=870.0001137  Z= -4.99941049            X=19742.686   Y=28652.914   计算得  S=870.0003175  Z=- 0.00041814 H=130.475            X=19737.837   Y=28651.697   计算得  S=870.0002748  Z= 4.999808656           X=19741.5912   Y=28722.0580   计算得  S=939.9999786  Z= -5.123024937           X=19736.4769   Y=28722.3564   计算得  S=939.9999862  Z=- 0.000027710 H=130.475           X=19733.4730   Y=28722.5317   计算得  S=940.0000238  Z= 3.00898694