5800线元法程序

1. 主程序

Lbl 4："1.SZ => XY"："2.XY => SZ"：？N：？S：Prog“SUB0”↙

1÷P→C: (P-R)÷(2HPR) →D:180÷∏→E:N=1 => Goto1：Goto2↙

Lbl 1:? Z：Abs(S-O) →W：Prog "SUB1"："XS="：X◢

"YS="：Y◢F-90→F：“FS=”：F▲DMS◢：Goto4↙

Lbl 2：？X：？Y：X→I：Y→J：Prog“SUB2”：O+W→S：“S=”：S◢“Z=”：Z◢Goto4↙

2. 正算子程序(SUB1)

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+Zcos（F）→X：Y+Zsin（F）→Y

2. 反算子程序(SUB2)

G-90→T：Abs((Y-V)cosT-(X-U)sinT) →W：0→Z：Lbl 0：Prog "SUB1"：

T+QEW(C+WD) →L：(J-Y)cos（L）-(I-X)sin（L）→Z：Abs（Z）<1E-6=>Goto1：W=W+Z：Goto 0←┘Lbl 1：Z=0：Prog "SUB1"：(J-Y)÷sin（F）→Z

SUB0 数据库子程序

Goto1↙同时保存多个曲线时的指针

Lbl 1：IF S<***（线元终点里程）:Then***→G（线元起点方位角）:***→O（线元起点里程）:***→U（线元起点X）:***→V（线元起点Y）:***→P（线元起点曲率半径）:***→R（线元终点曲率半径）: ***→H（线元起点至终点长度）：0或1、-1→Q：Return:IfEnd↙

Lbl 1：IF S<***（线元终点里程）:Then***→G（线元起点方位角）:***→O（线元起点里程）:***→U（线元起点X）:***→V（线元起点Y）:***→P（线元起点曲率半径）:***→R（线元终点曲率半径）: ***→H（线元起点至终点长度）：0或1、-1→Q：Return:IfEnd

………………………..

为了便于解读，每增加一个线元增加一行语句，每增加一条曲线增加一个Lbl，每增加一个工程增加一个文件。

三、使用说明

1、规定

(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，

Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。

(2) 当所求点位于中线时，Z=0；当位于中线左铡时，Z取负值；当位于中线中线右

侧时，Z取正值。

(3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。

(4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆

弧的半径。

(5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45

次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半

径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。

(6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的

值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等

于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。

2、输入与显示说明

输入部分：

1. SZ => XY

2. XY = > SZ

N ? 选择计算方式，输入1表示进行由里程、边距计算坐标；输入2表示由坐标反算里程和边距。

S ？正算时所求点的里程，反算时为所求点的近似里程，反算时输入的近似里程不能夸线元，若夸线元的话，需进行第二次反算，并以第一次的结果做为近似里程，否则反算的里程及边距不对，Z ？正算时所求点距中线的边距(左侧取负，值右侧取正值，在中线上取零)

X ？反算时所求点的X坐标

Y ？反算时所求点的Y坐标

显示部分：

XS=×××正算时，计算得出的所求点的X坐标

YS=×××正算时，计算得出的所求点的Y坐标

FS=×××正算时，所求点对应的中线点的切线方位角

S=×××反算时，计算得出的所求点的里程

Z=×××反算时，计算得出的所求点的边距

分两步计算：1、计算曲线要素 P=S²÷(24R);Q=S÷2-S³÷(240R²);T=(P+R)tan|0.5B|+Q;L=πRB÷

180+S ZH(C)=A-T;HY(D)=C+S;YH(H)=C+L-S;HZ(G)=C+L

式中：F----方位角；A----交点桩号；B-----路线偏角（左偏为负，右偏为正），R----半径；S---缓和曲线长

T----切线长；L-----曲线长

2、计算方位角：ZH=F;HY=F+S÷2R;HZ=F+B;YH=HZ-S÷2R;(当F＜0时S÷2R为负，大于0时为正)