≥ 0.001: (循环语句。先计算测点对应线路坐标,然后进行坐标转换,给切垂距D假设一个小于 1mm 的精度,满足要求后精确计算缓和曲线长和偏距 S ZH +L→ S: (计算测点对应精确里程“ S=”:S ⊿“ C=”:C ⊿?(显示精确里程及偏距 HHQX2 (计算出缓和曲线里程及偏距) 程序清单:? X:?Y: ?W:S HZ-X→L:D O:S HZ -(L-L 5÷ (40 ×R 2×l 02 2)+L 9÷ (3456 ×R 4×l 02 4 )) → P:W× (L 3÷ (6×R×l 02)-L 7÷ (336 ×R 3×l 02 3 ))→ Q: 90×L 2÷( π× R×l 02) → K: (X-P) × cos(-W × K)+(Y-Q) × sin(-W × K)→ D: (Y-Q) × cos(-W × K)-(X-P) × sin(-W × K)→ C: L-D → L: LpWhile Abs(D) ≥ 0.001: S HZ -L→ S:“ S=”:S ⊿“ C=”:C ⊿?变量注释: X、Y ——以 ZH 点或 HZ 点里程为 X 坐标的切线右手坐标系中测点坐标;P、Q ——测点对应线路中线点的坐标;L ——缓和曲线长;l 01、l 02 ——入、出缓和曲线全长( 可自行选择变量); S ——测点对应里程; S HZ、S HZ ——设计 HZ 点、 ZH 点里程(可自行选择变量); C ——测点到线路中线的偏距( 左为-, 右为+);D ——切垂距;R ——圆曲线半径; K ——缓和曲线的切线角 W ——判别符号(左转曲线为-1 ,右转曲线为 1)。附注: 1. 程序中变量可以根据计算器或者本人习惯改变; 2. 使用前必须将设计交桩坐标系转换为以 ZH 点或 HZ 点里程为 X 坐标的切线右手坐标系。尹书云 2016-1-29 更新