湖南新光工程有限公司三板溪项目部 李甲生

程序目的：依平曲线要素计算直线、圆曲线、缓和曲线的任意中桩、左、右桩坐标。
程序说明：K0：起始桩号 X0：起始X坐标 Y0：起始Y坐标 ALF：起始方位角 R：半径 LS：缓和曲线长 N：曲线左转N=1，右转N=2 K：待求桩号 LL、LR：左、右桩距离 Q：左、右桩与中线斜交角 求得XZ、YZ、XL、YL、XR、YR分别为中桩、左、右桩坐标。


一、 直线段
文件名：ZX (COMP)
程式： L”K0”:O”X0”:P”Y0”:W”ALF”:
Lbl 0:{K}:X”XZ”=O+(K-L)cosW◢
Y”YZ”=P+(K-L)sinW◢
{B}:S”XL”=X-B”LL”cos(W+Q) ◢
T”YL”=Y-Bsin(W+Q) ◢
{C}:U”XR”=X+C”LR”cos(W+Q) ◢
V”YR”=Y+Csin(W+Q) ◢
Goto 0


注：在程序执行过程中，赋给的要素变数的值被固定不变，可对变数（K、LL、LR）赋予不同值，迅速求得所需坐标。

二、 圆曲线段
文件名：YQX (COMP)
程式： L”K0”:O”X0”:P”Y0”:W”ALF”: Lbl 1:{K}:
J=(-1)^N (K-L)÷R×180÷π:D=2Rsin((-1)^N J÷2):
X”XZ”=O+Dcos(W+J÷2) ◢
Y”YZ”=P+Dsin(W+J÷2) ◢
{B}:S”XL”=X-B”LL”cos(W+J+Q) ◢
T”YL”=Y-Bsin(W+J+Q) ◢
{C}:U”XR”=X+C”LR”cos(W+J+Q) ◢
V”YR”=Y+Csin(W+J+Q) ◢
Goto 1


注：若没有直接HY点方位角，则ALFHY=ALFZH±Ls/2/R×180/π,(左转-，右转+)。



三、 缓和曲线
文件名：HHQX (COMP)
程式： L”K0”:O”X0”:P”Y0”:W”ALF”: M”LS”:Lbl 3:{K}:
I=(-1)^N×(K-L)^2÷M÷R÷6×180÷π:
D=(K-L)-(K-L)^5÷90÷(RM)^2:
X”XZ”=O+Dcos(W+I)◢
Y”YZ”=P+Dsin(W+I) ◢
{B}:S”XL”=X-B”LL”cos(W+3I+Q) ◢
T”YL”=Y-Bsin(W+3I+Q) ◢
{C}:U”XR”=X+C”LR”cos(W+3I+Q) ◢
V”YR”=Y+Csin(W+3I+Q) ◢
Goto 3


注：1、坐标计算方法是根据偏角法原理；
2、缓和曲线（ZH~HY或YH~HZ）以ZH（或HZ）为起始点；
3、平曲线左转（ZH~HY段N=1，YH~HZ段N=2），曲线右转（ZH~HY段N=2，YH~HZ段N=1）。



Casio fx-4800P计算机具有强大的储存记忆、程序设计功能，能为众多专业技术工作提供大量复杂程序计算。笔者根据工程测量中的需要汇编以上程序，本着严密准确、计算程式简洁，避免冗长复杂化的计算过程。以上三部分平曲线计算程序有待继续研究成更简便的三合一程序。笔者另编写Casio fx-4800P方位角计算程序、竖曲线设计高程计算程序，及依坐标求算里程桩号及至中桩距离的程序，Excel坐标计算程序、Excel平曲线要素表计算程序、Excel竖曲线计算程序均具大批量计算并打印功能， 程序准确实用、方法简捷，下次见面。望大家不吝赐教，谢谢！