5800缓和曲线计算程序

正算主程序(ZS)
Lb1 0
FIX4：?S：?Z：Prog “PM-SJ”： Abs(S-O) → W：Prog "SUB1"： "XS="：X◢
"YS="：Y◢F-90→F(需要时可以让他显示，不要时不必输入)：Prog“SQX-SJ”: Prog“SQX”：“H=”：H◢
Goto 0
S代表正算桩号
Z代表偏距左负右正
XS所求点的X坐标
YS所求点的Y坐标
H所求桩号设计高程

反算主程序(FS)
Lb1 0
FIX4：?X：?Y：Prog “PM-SJ”：X→I:Y→J: Prog "SUB2"： "S="：O+W→S◢
"Z="：Z◢
Prog“SQX-SJ”: Prog“SQX”：“H=”：H◢
Goto 0

X测点的X坐标
Y测点的Y坐标
S测点的桩号
Z测点的偏距左负右正
H测点所在桩号的设计高程

正算子程序(SUB1)：
1÷P→ 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-90→T
Abs((Y-V)cosT-(X-U)sin（T）) → W： 0→ Z：Lbl 0：Prog "SUB1"： T+QEW(C+WD) →L： (J-Y)cosL-(I-X)sinL → Z：ifAbsZ<1E-6：thenGoto1：ElseW+Z → W：Goto 0：IfEnd←┘：
Lbl 1:0→ Z：Prog "SUB1"： (J-Y)÷sinF→ Z

子程序（平面线形数据库）PM-SJ
直线
ifS ≥500（线元起点里程）AndS＜769.256（线元止点里程）：Then19942.837 →U（线元起点X坐标）：28343.561 →V（线元起点Y坐标）：500 →O（线元起点里程）：125°16°31°→G（线元起点方位角）：269.265 →H（线元长度）：10^（45）→P（线元起点曲率半径）：10^（45）→R（线元终点曲率半径）：0 →Q（线元左右偏标志：左负右正）：IfEnd←┘

缓和曲线
ifS ≥769.256（线元起点里程）AndS＜806.748（线元止点里程）：Then19787.34 →U（线元起点X坐标）：28563.378→V（线元起点Y坐标）：769.256→O（线元起点里程）：125°16°31°→G（线元起点方位角）：37.492→H（线元长度）：10^（45）→P（线元起点曲率半径）：221.75→R（线元终点曲率半径）：-1→Q（线元左右偏标志：左负右正）：IfEnd←┘

圆曲线
ifS ≥806.748（线元起点里程）AndS＜919.527（线元止点里程）：Then19766.566→U（线元起点X坐标）：28594.574→V（线元起点Y坐标）：806.748→O（线元起点里程）：120°25°54.07°→G（线元起点方位角）：112.779→H（线元长度）：221.75→P（线元起点曲率半径）：221.75→R（线元终点曲率半径）：-1→Q（线元左右偏标志：左负右正）：IfEnd←┘

缓和曲线
ifS ≥919.527（线元起点里程）AndS＜999.812（线元止点里程）：Then19736.072 →U（线元起点X坐标）：28701.893 →V（线元起点Y坐标）：919.527→O（线元起点里程）：91°17°30.

63°→G（线元起点方位角）：80.285→H（线元长度）：221.75→P（线元起点曲率半径）：10^（45）→R（线元终点曲率半径）：-1→Q（线元左右偏标志：左负右正）：IfEnd←┘

直线
ifS ≥999.812（线元起点里程）：Then19744.038 →U（线元起点X坐标）：28781.659→V（线元起点Y坐标）：999.812→O（线元起点里程）：80°40°50°→G（线元起点方位角）：100→H（线元长度）：10^（45）→P（线元起点曲率半径）：10^（45）→R（线元终点曲率半径）：0 →Q（线元左右偏标志：左负右正）：IfEnd←┘

子程序（竖曲线数据库）SQX-SJ
ifS≤999.812（竖曲线终点里程）：then0.0357（前坡度：上坡正下坡负）→E：0.0227（后坡度：上坡正下坡负）→F：600→R（竖曲线半径）：780→G（变坡点里程）：110→C（变坡点高程）：IfEnd←┘

子程序（竖曲线计算公式）SQX
E-F→J：ifJ＞0：then-R→R：ElssR→R：Abs（R*J÷2）→T：IfEnd←┘
ifS≤G-T：thenG-S→L：C-LE→H：IfEnd←┘
ifS＞G-T And S≤G： thenG-S→L：（S+T-G）2÷2÷R→N：C-LE+N → H：IfEnd←┘
ifS＞G And S≤G+T： then S-G→L：（G +T- S）2÷2÷R→N：C+LF+N → H：IfEnd←┘
ifS＞G+T ：thenS-G→L： 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) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。 以上程序为个人改编而得，如有好的建议及不妥处请留言。如有好的程序，也发一份我共享啊！