5800正反算、边坡放样、隧道测量程序.doc

正算主程序(ZS) FIX4：?S：?Z：Prog “PM-SJ”： Abs(S-O) W：ProgSUB1： XS=：XYS=：YF-90F(需要时可以让他显示)：Prog“SQX-SJ”: Prog“SQX”：“H=”：H反算主程序(FS)FIX4: ?S：?X：？Y：Prog“PM-SJ”： X I： Y J：ProgSUB2： S=：O+WSZ=：ZProg“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：ifPB：thenGoto1：ElseGoto2Lb1 1：ifPB：thenL+C（P-B）U：P-B F：IfEndifPB+G：thenL+CG+K+D（P-B-G）U：P-B-G F：IfEndifPB+G+M：thenL+GC+2k+MD+E（P-B-G-M）U：P-B G-MF：IfEnd：Goto3Lb1 2：ifPB：thenL+I（B -P）U： B -P F：IfEndifPB-Q：thenL+IQ+K+J（B-Q-P）U：B-Q-PF：IfEndifPB-Q-R：thenL+IQ+2k+JR+O（B-Q-R-P）U：B-Q-R-PF：IfEnd：Goto3Lb1 3：U-AbsZV：ifZ0：thenZ-VZ：Else z+vz：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.5A：1.65B：0.89G：2.27C：3.3D：2.41E：0.51J：“H1”：F：ifFH+A：then E-AbsZW：IfEnd：ifFH+A And FH+A+B：then（D-（F-H-A）-G-AbsZW：IfEnd：ifFH+A+B：then（C-（F-H-A-J）-AbsZW：IfEnd：“W=”：W （水平方向的超欠挖：正为欠,负为超）Goto0正算子程序(SUB1)1P 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-90T： Abs(Y-V)cosT-(X-U)sin（T）) W： 0 Z：Lbl0：ProgSUB1： T+QEW(C+WD) L： (J-Y)cosL-(I-X)sinL Z：ifAbsZ1E-6：thenGoto1：ElssW+Z W：Goto0：IfEndLbl1： 0 Z：ProgSUB1： (J-Y)sinF Z子程序（平面数据库）PM-SJifS 500（线元起点里程）AndS769.256（线元止点里程）：Then19942.837 U（线元起点X坐标）：28343.561 V（线元起点Y坐标）：500 O（线元起点里程）：1251631G（线元起点方位角）：269.265 H（线元长度）：10（45）P（线元起点曲率半径）：10（45）R（线元终点曲率半径）：0 Q（线元左右偏标志：左负右正）：IfEndifS 769.256（线元起点里程）AndS806.748（线元止点里程）：Then19787.34 U（线元起点X坐标）：28563.378V（线元起点Y坐标）：769.256O（线元起点里程）：1251631G（线元起点方位角）：37.492H（线元长度）：10（45）P（线元起点曲率半径）：221.75R（线元终点曲率半径）：-1Q（线元左右偏标志：左负右正）：IfEndifS 806.748（线元起点里程）AndS919.527（线元止点里程）：Then19766.566U（线元起点X坐标）：28584.574V（线元起点Y坐标）：806.748O（线元起点里程）：1202554.07G（线元起点方位角）：112.779H（线元长度）：221.75P（线元起点曲率半径）：221.75R（线元终点曲率半径）：-1Q（线元左右偏标志：左负右正）：IfEndifS 919.527（线元起点里程）AndS999.812（线元止点里程）：Then19736.072 U（线元起点X坐标）：28701.893 V（线元起点Y坐标）：919.527O（线元起点里程）：911730.63G（线元起点方位角）：80.285H（线元长度）：221.75P（线元起点曲率半径）：9579.228R（线元终点曲率半径）：-1Q（线元左右偏标志：左负右正）：IfEndifS 999.812（线元起点里程）：Then19744.038 U（线元起点X坐标）：28781.659V（线元起点Y坐标）：999.812O（线元起点里程）：804050G（线元起点方位角）：100H（线元长度）：10（45）P（线元起点曲率半径）：10（45）R（线元终点曲率半径）：0 Q（线元左右偏标志：左负右正）：IfEnd子程序（竖曲线数据库）SQX-SJifS999.812（竖曲线终点里程）：then0.0357（前坡度：上坡正下坡负）E：0.0227（后坡度：上坡正下坡负）：600R（竖曲线半径）：780G（变坡点里程）：110C（变坡点高程）：IfEnd子程序（竖曲线计算公式）SQXE-FJ：ifJ0：then-RR：ElssRR：Abs（R*J2）T：ifSG-T：thenG-SL：C-LEH：IfEndifSG-T And SG： thenG-SL：（S+T-G）2RN：C-LE+N H：IfEndifSG And SG+T： then S-GL：（G +T- S）2RN：C-LF+N H：IfEndifSG+T ：thenS-GL： 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、起点坐标X0Y0、起点切线方位角F0、线元长度LS、起点曲率半径R0、止点曲率半径RN、线元左右偏标志Q）如下：S0X0Y0F0LSR0RNQ500.00019942.83728343.5611251631.00269.2561E451E450769.25619787.34028563.3781251631.0037.4921E45221.75-1806.74819766.56628594.5741202554.07112.779221.75221.75-1919.52719736.07228701.893911730.6380.285221.759579.228-1999.81219744.03828781.659804050.00100.0001E451E4501、正算（注意：略去计算方式及线元要素输入，请自行根据所求点所在的线元输入线元要素）S=700Z=-5计算得XS=19831.41785YS=28509.72590 H=107.144S=700Z=0计算得XS=19827.33592YS=28506.83837 H=107.144S=700Z=5计算得XS=19823.25398YS=28503.95084 H=107.144 S=780Z=-5计算得XS=19785.25749YS=28575.02270 H=112.759S=780Z=0计算得XS=19781.15561YS=28572.16358 H=112.759S=780Z=5计算得XS=19777.05373YS=28569.30446 H=112.759S=870Z=-5计算得XS=19747.53609YS=28654.13091 H=112.759S=870Z=0计算得XS=19742.68648YS=28652.91379 H=112.759S=870Z=5计算得XS=19737.83688YS=28651.69668 H=112.759S=940Z=-5.123计算得XS=19741.59118YS=28722.05802 H=130.475S=940Z=0计算得XS=19736.47687YS=28722.35642 H=130.475S=940Z=3.009计算得XS=19733.47298YS=28722.53168 H=130.4752、反算X=19831.418Y=28509.726计算得S=699.9999974Z=-5.00018164X=19827.336Y=28506.838计算得S=699.9996493Z=0.000145136 H=107.144X=19823.25398Y=28503.95084计算得S=699.9999985Z=5.000003137X=19785.25749Y=28575.02270计算得S=780.0000035Z=-5.000001663X=19781.15561Y=28572.16358计算得S=780.0000025Z=-0.000002979H=112.759X=19777.05373Y=28569.30446计算得S=780.0000016Z=4.99999578X=19747.536Y=28654.131计算得S=870.0001137Z=-4.99941049X=19