王德明 第 1 页 2017-9-3匝道上回旋线任意点坐标方位角计算源程序：“ZDHQ ” C“A2”: H=1 M=0 ◣ H=2 M=C/R ◣ H=1 G=N“ZH（G）HY ”: ≠ G=N“ZH（ G）HY”+180D/Rл ◣ G <360 G“ZH（G ）YH”◢ ≠ G“ZH（G）YH”=G-360◢ ◣ Lbl 1：{ZLAB}：Z=0 Goto2 ◣ P=1 L=M+L-K ◣ P=-1 L=M-L+K ◣ A：B： T=G-90PM2/Cл： E=T：E=E+90PL 2/Cл：E <360 E◢ ≠ E=E-360 ◢ ◣ O=（L∧3-M∧3）/6C-（L ∧7-M∧7）/336C ∧3+（L∧11-M∧11）/42240C ∧5-（L∧15-M∧ 15）/9676800C∧7：F=L-M-（L∧5-M∧5）/40C 2+（L∧9-M∧ 9）/3456C∧ 4：U“XN”=W+PF cosT-OsinT ◢ Z“YN”=V+PFsinT+OcosT ◢ X=U： Y=Z： Prog“F”： Prog“BZ”：Goto 1：Lbl 2：Z=0 W=O：V=F：P=1 K=L-M+K ◣ P=-1 K=M-L+K ◣ F=G：“END”“XY” U“XP”=X+ Rec(S，Q ◢ Z“YP”=Y+J“ F ” Pol（X-A，Y-B：J＜0 J=J+360 ◢ ≠ J ◢ ◣ I ◢ “BZ ” X=U: Y=Z： O=U: F=Z: {S}: T=1 ◢ Q=E-90 : Prog“XY”◢ X=U: Y=Z:Prog“F”: {S}: T=T+1 ◢ X=O: Y=F: Q=E+90：Prog“XY”◢ X=U: Y=Z: Prog“F”参考〈土木工程施工测量手册〉中国建材工业出版社 2002 年 5 月版说明：A2 = L0R1R2 /（R 1-R2） 当 R1为无穷大时则 R1=100万 R 2A2 — 曲率半径变化率的平方H — H= 1 计算第一缓和曲线 H=2 计算第二缓和曲线 ZH（G）HY — H = 1 时过 ZH H = 2 时过 HY 起算方位角P — 起点 RK 半径＞终点 RL 半径时输入 1，反之输入- 1L — 回旋线待求点里程K — 回旋线起点里程E — 所求点切线方位角A、B — 测站点坐标W、 V — 回旋线起点坐标 S — 中线至边线距离R — 已算过的前圆曲线半径， （ 左负右正）D — 已算过的前圆曲线长度Z = 0 结束本曲线计算，显示 END求： ZH（G）HY — 过本曲线 ZH 或过 YH 点切线方位角E — 所求点切线方位角 XN、YN、J、I — 所求桩号点中桩及左、右边桩坐标 及 测站至所求点方位角、边长T=1 为所求桩号点左边桩坐标T=2 为所求桩号点右边桩坐标算例：计算第一缓和曲线：已知： A 2= 902 H= 1 ZH = 134.6108889P= -1 L= 50 K= 47.86W= 933.866 V= 312.967 S=20A= 900 B= 300成果： ZH（G）YH =134.6108889 过本曲线 ZH、YH 切线方位E = 134.5946919 （所求点切线方位角）中桩：XP= 932.363 =24°07 ′13.34 ″YP=314.491 S = 35.459左桩：XP= 946.605 =31°28 ′32.17 ″YP= 328.532 S = 54.645右桩：XP=918.121 = 1°25 ′07.93 ″YP=300.449 S = 18.127计算第二缓和曲线：A2= 1202 H= 2 R= -200HY（G）= 75.15295764 D= 58.78P= 1 L= 280 K= 272.985W= 905.365 V= 514.788 S= 20A= 900 B= 500成果：ZH（G）YH = 58.31372804 过本曲线 YH 切线方位角 E = 56.40197921 （所求点切线方位角）中桩：XP= 909.150 = 66°08 ′45.38 ″YP= 520.694 S = 22.626左桩：XP= 925.809 = 20°27 ′17.09 ″YP= 509.626 S = 27.546右桩：XP= 892.492 = 103°18 ′03.15 ″王德明 第 2 页 2017-9-3YP=531.761 S = 32.636。