5800测量程序全线带高程分围岩级别计算20091025修改.doc

卡西欧5800测量程序一、 坐标计算源程序 1.主程序(QXjs-000) Lbl 0：1.SZ = XY：2.XY = SZ：?N：S=?S：Prog QXjs SUB0：D=?D: Prog B2：IfN=1：ThenGoto1：ElseGoto 2：Lbl 1：S =?S：Z =?Z：Abs(S-O) W：Prog SUB1：XS=：X YS=：YFWJ=：F-90DMSProg B1:H=：Z8Goto 0Lbl 2：X=?X：Y=?Y： XI：YJ：Prog SUB2：S：O+WZ：ZAbs (O+W)S :Prog B1:H=：Z8Abs (Z-2.4)Z: CG=?A: NX=?B: Prog RRRR:L=：L Goto 02.平曲线数据库子程序(QXjs-SUB0)：IfS*(本段线元终点里程): Then*G(线元起点方位角)：*O(线元起点里程)：*U(线元起点X)：*V(线元起点Y)：*P(线元起点曲率半径)：*R(线元终点曲率半径)：*H(线元起点至终点长度)：*Q(线元左为-1右+1)：Return:IfEndIfS*: Then*G：*O：*U：*V：*P：*R：*H：*H ：Return:IfEndIfS带长短链的线元，输入时从长短链处分开输线元要素，3. 正算子程序(SUB1) 1PC：(P-R)(2HPR) D：180E：0.1739274226A：0.3260725774B：0.0694318442K： 0.3300094785L：1-LF：1-KM：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)+90F： X+Zcos(F)X： Y+Zsin(F)Y（想计算任意偏角时，可以将A值直接写成数值，再最后一步里F后面加-90+A，就可以了，不过反算A值要输90，左偏-角右+角)4. 反算子程序(SUB2) G-90T： Abs(Y-V)cos(T)-(X-U)sin(T) W： 0ZLbl 0：Prog SUB1： T+QEW(C+WD) L： (J-Y)cos(L)-(I-X)sin(L)Z：IfAbsZ10(-6) ：ThenGoto1：Else W+ZW：IfEnd：Goto 0 Lbl 1：0Z：Prog SUB1： (J-Y)sin(F )Z5高程数据库（B2）IfD=*: Then *Z1(竖曲线变坡点里程) :*Z2(交点高程)：*Z3(前段坡率)：*Z4(后段坡率)：*Z5(切线长): *Z6(竖曲线半径)：*Z7(圆弧在上面为-1，圆弧在下面为+1): *Z9 (本段线元起点里程): *Z10(本段线元终点里程)：6 高程计算程序(B1)IfS Z9:Then Cls:Stop:IfEndIfSZ1-Z5: Then Z2-(Z1-S)Z3Z8 :Return:IfEndIfSZ1:Then Z2-(Z1-S)Z3Z7(Z5-Z1+S)2 (2Z6) Z8: Return:IfEndIfSZ1+Z5: Then Z2+(S-Z1)Z4- Z7(Z5+Z1-S)2 (2Z6) Z8: Return:IfEndIfSZ1+Z5: Then Z2+(S-Z1)Z4Z8 :IfEndIfS Z10:Then Cls:Stop:IfEnd7半径计算（RRRR）围岩类型IfB=3: Then IfAZ8+1.3：Then6.98-（Z-0.00）2+(A-Z8-2.27)2）L: Return:IfEndIfAZ8-0.78: Then(6.982-0.972)-0.28(Z8+1.3-A) 2.08-ZL: Return:IfEndIfZ4.827: Then2.91-(Z-3.711)2+(Z8+1.0715-A)2)L: Return:IfEndIfZ4.827: Then15.27-(Z8+12.87-A)2+Z2) L: Return:IfEnd: Return:IfEndIfB=4: Then IfAZ8+1.3：Then7.21-（Z-0.00）2+(A-Z8-2.27)2）L: Return:IfEndIfAZ8-0.78: Then(7.212-0.972)-0.28(Z8+1.3-A) 2.08-ZL: Return:IfEndIfZ4.875: Then3.06-(Z-3.711)2+(Z8+1.0715-A)2)L: Return:IfEndIfZ4.875: Then15.42-(Z8+12.87-A)2+Z2) L: Return:IfEnd: Return:IfEndIfB=5: Then IfAZ8+1.3：Then7.39-（Z-0.00）2+(A-Z8-2.27)2）L: Return:IfEndIfAZ8-0.78: Then(7.392-0.972)-0.28(Z8+1.3-A) 2.08-ZL: Return:IfEndIfZ5.042: Then3.49-(Z-3.711)2+(Z8+1.0715-A)2)L: Return:IfEndIfZ5.042: Then15.85-(Z8+12.87-A)2+Z2) L: Return:IfEnd: Return:IfEndIfB=2: Then IfAZ8+0.3: Then6.46-（Z-0.00）2+(A-Z8-2.27)2）L: Return:IfEndIfAZ8+0.3: Then2.3-(Z-3.958)2+(Z8+0.9916-A)2)L: Return:IfEnd: Return:IfEndIfB=1: Then IfAZ8+5.902:Then7.66-（Z+0.00）2+(A-Z8-2.07)2）L: Return:IfEndIfAZ8+1.3: Then7.26-(Z-0.35)2+( A-Z8-2.27)2）L: Return:IfEndIfAZ8-0.78:Then(7.262-0.972)-0.285(Z8+1.3-A) 2.08-Z+0.35L: Return:IfEnd二 说明 （1）以线路前进方向区分左右，当曲线在左时P、R取负值，当曲线半径在右时P、R取正值，当曲线半径无穷大（即直线）时，P、R以10的45次代替。 （2）当所求点位于中线时Z=0:左侧取负值右侧取正值 （3）当线元为圆曲线时无论起点终点和什么线元相接，半径均取圆曲线半径 （4）当线元为直线时，其起点、止点的曲率半径为无穷大，以 10的45次代替。 (5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半 径为无穷大，以10的45次代替；与圆曲线相接时，曲率半 径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大， 以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。 (6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值； 与圆曲线相 接时，曲率半径等 于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲 线相接时，曲率半径等于圆曲线的半径。 (7) 带长短链的线元，输入时从长短链处分开输线元要素，第一步先扩充变量输50DimZ，第二步再输