道路维座标计算.docVIP

09-10-2　楼主：yaochong　作者：yaochong 本文已被浏览　3554　次 5800道路三维坐标计算 Fx-5800道路三维坐标计算器程序作者：姚 崇QXJS-000 主程序（文件名）“1.SZ=NE”:“2.NE=SZ”:?Q↙Lbl 4:?S:“PJ=”?K:Prog“QXJS-SUB0”↙Lbl 0:Q=1 = Goto1:Q=2 = Goto2:↙Lbl 1:?Z:Prog“QXJS-SUB1”: “F=”:F◢DMS◢ “N=”:N◢ “E=”:E◢Prog“GCJS”: Goto4↙Lbl 2: “N=”:?B: “E=”:?C:B→N: C→E:Prog“QXJS-SUB2”: “S=”:S◢ “Z=”:Z◢??Goto4↙ QXJS-SUB0??数据库子程序（文件名）Goto1↙?????????????? 同时保存多个曲线时的指针Lbl 1IF S***（线元终点里程）:Then***→A（线元起点方位角）:***→O（线元起点里程）:***→U（线元起点X）:***→V（线元起点Y）:***→P（线元起点曲率半径）:***→R（线元终点曲率半径）: ***→L（线元起点至终点长度）: Return:IfEnd↙IF S***:Then***→A:***→O:***→U:***→V:***→P:***→R: ***→L: Return:IfEnd↙………………………..为了便于解读，每增加一个线元增加一行语句，每增加一条曲线增加一个Lbl，每增加一个工程增加一个文件。QXJS-SUB1?? 正算子程序（文件名）0.5（1÷R-1÷P）÷L→D:S-O→X↙U+∫(cos(A+(X÷P+DX2)×180÷π,0,X)→N↙V+∫(sin(A+(X÷P+DX2)×180÷π,0,X)→E↙A+(X÷P+DX2)×180÷π→F↙N+Zcos(F+K) →N:E+Zsin(F+K) →EQXJS-SUB2??反算子程序（文件名）Lbl 1:0→Z：1→Q：Prog“QXJS-SUB0”: Prog“QXJS-SUB1”↙Pol(N-B+10^(-46), E-C+10^(-46)):Isin(F-K-J) →W:S+W→S↙Abs(W)0.0001 = Goto1↙Lbl 2: 0→Z：Prog“QXJS-SUB1”:(C-E) ÷sin(F+K) →ZGCJS(高程计算—文件名)Lbl 0： Prog“SJK”： 0.005R Abs(I -L) →T：T2÷2÷R→E:Z[1]-T→A：Z[1]+T→B：S-A→C：B-S→D：S-Z[1]→E←┘If L I：Then 1→J：Else -1→J：IfEnd←┘If SA：Then 0→C：Else SB =0→D：IfEnd←┘If S≦Z[1]：Then I→Y：C→G：Else S Z[1]=L→Y：D→G： IfEnd←┘Lbl 1：Fix 3：“ZP=”：H+EY÷100+JG2÷2÷R◢ Goto 0←┘SJK(高程计算子程序—数据库)If S≦下一竖曲线起点And S 上一竖曲线终点：Then??……→Z[1](变坡桩号): ……→H（变坡高程）: ……→R（竖曲线半径）:……→I（第一纵坡坡度，带正负号）: ……→L(第二纵坡坡度，带正负号): IfEnd←┘坐标正反算输入与使用使用说明1、规定(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当曲线半径在左时，P、R取负值，当曲线半径在右时，P、R取正值，当曲线半径为无穷大（即直线）时，P、R以10的45次代替。(2) 当所求点位于中线时，Z=0；当位于中线左铡时，Z取负值；当位于中线中线右侧时，Z取正值。当PJ=90时，表示与线路正交(3) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。(5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次方代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。(6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。2、输入与显示说明输入部分：??????1. SZ = XY??????2. XY = SZ ??????Q ???选择计算方式，输入1表示进行由里程、边距计算坐标 ；输入2表示由坐标反算里程和边距。高