线元法

线元法

1：ZHY主程序

“1.SZ→XZ”：“2.XY→SZ”： If N=1：Then goto 1 ：Else goto 2：lfEnd Lbl 1：‘S=’?S：“Z=”?Z：Prog〝QXSJK〞: 1÷P→C :（P－R）÷（2HPR）→D : 180÷π→E: Abs（S－O）→W:Prog〝SUB 1〞:〝XS〞: X→X▲

〝YS=〞: Y→Y▲ 〝FS=〞: F－90→F▲ Prog〝GAO〞: goto 1

LBI 2 :〝X=〞?X :〝Y=〞?Y :X→I:Y→J：Prog〝QXSJK〞: 1÷P→C :（P－R）÷（2HPR）→D : 180÷π→E: :Prog〝SUB 2〞: 〝S〞: O+W→S▲ 〝Z=〞: Z▲

Prog〝GAO〞: 〝SC〞?C 〝R=〞:5.-√((C-G)2+Z2) ▲ goto 2

“L=”：?L :199.965-(L-49040)*0.006+0.00 →H ◢ 设计高程 R-√（P²+G²）

2：SUB1（缓和曲线正算子程序，不能运行）

0.1739274226→A : 0.3260725774→B : 0.0694318442→K : 0.3300094782→L : 1－L→F : 1－K→M : U+W（A Cos（G+QEKW（C+KWD））+B Cos（G+QELW（C+LWD））+ B Cos（G+QEFW（C+FWD））+A Cos（G+QEMW（C+MWD）））→X :V+W（A Sin（G+QEKW（C+KWD））+B Sin（G+QELW（C+LWD））+ B Sin（G+QEFW（C+FWD））+A Sin（G+QEMW（C+MWD）））→Y : G+QEW（C+WD）+90→F : X+Z cos（F）→X : Y+Z Sin（F）→Y

3 : SUN 2（缓和曲线反算子程序，不能运行）

G－90→T : Abs（（Y－V）Cos（T）－（X－U）Sin（T））→W : 0→Z : LBI 0 : Prog

〝SUB 1〞: T+QEW（C+WD）→L :（J－Y）Cos（L）－（I－X）Sin（L）→Z : If Abs（Z）< 1×10－6 : Then Goto 1 : ELse W+Z→W : Goto 0 : Ifend

LBI 1 : 0→Z : Prog〝SUB 1〞:（J－Y）÷Sin（F）→Z

4 : QX SJK （缓和曲线数据库子程序，不能运行）

If S≥7000 And S<8552.052 : Then 39351.594→U : 72418.7097→V : 7000→O :257.8746719→G : 2000→H : 1×1045→P : 1×1045→R : 0→Q : ELse If S≥8552.052 And S<8752.052 : Then 39025.584→U : 70901.283→V : 8552.052→O: 257.8746719→G : 200→H : 1×1045 →P : 1800→R : 1→Q : ELse If S≥8752.052 And S<9900.413 : Then 387.2071→U : 70705.0275→V : 8752.052→O : 261.0577708→G : 1148.361→H : 1800→P : 1800→R : 1→Q : ELse If S≥9900.413 : Then 39170.3263→U : 69590.9925→V : 9900.413→O : 297.6112367→G : 200→H : 1800→P : 1×1045→R : 1→Q : Ifend : Ifend : Ifend : Ifend : Return

注意：0 : 表示零。O ：表示大写字母O。子程序结束后一定要加上Return

5：SQXBG主程序名

Lbi1：?H：Prog＂SJK＂：C-D→F：Abs(RF÷2)→T：Rabs（F÷F）→R:

IF H≤B-T:Then0→K：Eise if H≥B+T：Then 0→K：D→C：Eise if H≤B：Then H-B+T→ K：Eise T-(H-B) →K：D→C：IFEnd：IFEnd：IFEnd A-（B-H）C-K2÷（2R）→G▲ Goto1 子程序：SJK

If H≤第二竖曲线起点桩号（5569.108）: Then A=第一竖曲线交点高程（820.968）：B=第一竖曲线交点桩号（4980）：C=第一竖曲线前坡度（.03）：D=第一竖曲线后坡度（-.02902）：R=第一竖曲线半径（13000）：

ELse If H≤第三竖曲线起点桩号（7174.4）: Then A=第二竖曲线交点高程（798.331）：B=第二竖曲线交点桩号（5760）：C=第二竖曲线前坡度（-.02902）：D=第二竖曲线后坡度（.0187）：R=第二竖曲线半径（8000）：

ELse If H≤第四竖曲线起点桩号（78.584）: Then A=第三竖曲线交点高程（828.438）：B=第三竖曲线交点桩号（7370）：C=第三竖曲线前坡度（.0187）：D=第三竖曲线后坡度（.035）：R=第三竖曲线半径（24000）：: Ifend : Ifend : Ifend……………………继续添加要素

Lbi 1

1、A：交点高程 2、B：交点桩号

3、C：前坡度（上坡为正，下坡为负） 4、D：坡度（上坡为正，下坡为负） 5、R：竖曲线半径 6、H：待求点桩号 7、E:下一竖曲起点桩号

说明：如果只需要输入一段竖曲线，输入到半径的时候，后面直接加入 ：Goto1⊿ 如下：

H≤第二竖曲线起点桩号=＞A=第一竖曲线交点高程（820.968）：B=第一竖曲线交点桩号（4980）：C=第一竖曲线前坡度（.03）：D=第一竖曲线后坡度（-.02902）：R=第一竖曲线半径（13000）：Goto1⊿ H≤7174.4=＞A=第二竖曲线交点高程（798.331）：B=第二竖曲线交点桩号（5760）：C=D：D＝第二竖曲线后坡度（？？？？？）：R=第二竖曲线半径（8000）：Goto1⊿

三、使用说明 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表示由坐标反算

里程和边距。 X0 ？线元起点的X坐标 Y0 ？线元起点的Y坐标 S0 ？线元起点里程 F0 ？线元起点切线方位角 LS ？线元长度

R0 ？线元起点曲率半径 RN ？线元止点曲率半径

Q ？ 线 元左右偏标志(左偏Q=-1，右偏Q=1，直线段Q=0） S ？ 正算时所求点的里程 Z ？正算时所求点距中线的边距(左侧取负，值右侧取正值，在中线上取零) X ？反算时所求点的X坐标

Y ？反算时所求点的Y坐标 显示部分：

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

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

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

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

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