1#
用COSLO fx-4800p计算非对称曲线中(边)桩上任意点坐标程序
绥德县农村公路管理站 李刚
1.主点曲线要素计算程序(主程序)文件名 FDCQX ZB FY
L b1 0：{R,O,M,N,J}:
R〝R=〞：O〝PJ=°〞：M〝Ls(1) =〞：N〝Ls(2) =〞：
P=M2÷24R－M^4÷2688R^3：
Q=N2÷24R－N^4÷2688R^3：
S=M÷2－M^3÷240R2：
T=N÷2－N^3÷240R2：
Z[1] 〝T1〞=S+(R+Q－(R+P)cosAbsO)÷sinAbsO◢
Z[2] 〝T2〞=T+(R+P－(R+Q)cosAbsO)÷sinAbsO◢
L〝LY〞：=RπAbsO÷180+M÷2+N÷2：
H〝EH〞=(R+(P+Q) ÷2) ÷cos(AbsO÷2)-R◢
D= Z[1]+ Z[2] －L◢{CABZ}:
V〝ZH〞=C〝JD〞－Z[1] ◢
Z[3] 〝HY〞=V+M◢
Z[4] 〝QZ〞= Z[3]+(L－M－N) ÷2◢
Z[5] 〝YH〞=V+L－N◢
Z[6] 〝HZ〞= Z[5] +N◢
A〝JD(X)=〞:B〝JD(Y)=〞:Z〝JD(FWJ)°=〞◢
L b1 1:{D}：D〝ZZ：1,BZ：2〞：D=1＝＞Prog〝ZZ ZB FY〞◢≠＞D=2＝＞Prog〝BZ ZB FY〞◢≠＞Goto0◣◣Goto1↙
2.中桩坐标计算及放样程序 文件名 ZZ ZB FY
1)直线（上一曲线终点HZ点至计算曲线起点ZH点之间的直线上任意点的坐标）
{K}：K〝K=〞：KX〝ZZ(X)〞=A+ Z[1] cos(Z+180)+I◢
Y〝ZZ(Y)〞=B+ Z[1] sin(Z+180)+J◢
G=Z：
2)本计算曲线内ZH~HY段曲线内任意点的坐标
≠＞K≤Z[3] ＝＞Pol((K－V －(K－V)^5÷(40R2M2)，(K－V)^3÷(6RM) －(K－V)^7÷336R^3M^3) ：Rec(I，Z+ (O÷AbsO)J ) ：
X〝ZZ(X)〞=A+ Z[1] cos(Z+180)+I◢
Y〝ZZ(Y)〞= B+ Z[1] sin(Z+180)+J◢
G=Z+(O÷AbsO )(90(K－V)2÷(πRM)) ：
3) 本计算曲线内HY~QZ段曲线内任意点的坐标
≠＞K≤Z[4] ＝＞Pol((Rsin(90(2(K－Z[3])+M) ÷πR)+S) ，R(1－cos(90(2(K－Z[3])+M) ÷πR))+P) ：Rec(I，Z+ (O÷AbsO )J) ：
X〝ZZ(X)〞=A+ Z[1] cos(Z+180)+I◢
Y〝ZZ(Y)〞= B+ Z[1] sin(Z+180)+J◢
G=Z+(O÷AbsO )(90M÷(πR)+180(K－Z[3] )÷πR) ：
4) 本计算曲线内QZ~YH段曲线内任意点的坐标
≠＞ K≤Z[5] ＝＞Pol((Rsin(90(2(Z[5] －K)+N) ÷πR)+T)，R(1－cos(90(2(Z[5] －K)+N) ÷πR))+Q) ：Rec(I，Z+O－180－(O÷AbsO )J) ：
X〝ZZ(X)〞=A+ Z[2] cos(Z+O)+I◢
Y〝ZZ(Y)〞= B+ Z[2] sin(Z+O)+J◢
G= Z+O－(O÷AbsO )(90N÷(πR)+180(Z[5] －K)÷πR) ：
5) 本计算曲线内YH~HZ段曲线内任意点的坐标
≠＞ K≤Z[6] ＝＞Pol((Z[6]－K－(Z[6]－K)^5÷40R2N2)) ， (Z[6]－K)^3÷6RN) ：Rec(I，Z+O－180－(O÷AbsO )J) ：
X〝ZZ(X)〞=A+ Z[2] cos(Z+O)+I◢
Y〝ZZ(Y)〞= B+ Z[2] sin(Z+O)+J◢
G=Z+O－(O÷AbsO )(90(Z[6] －K)2÷πRN) ：
6)直线（计算曲线终点HZ点至下一曲线起点ZH点之间的直线上任意点的坐标）
≠＞K> Z[6] ＝＞Rec(K－Z[6]，Z+O ) ：
X〝ZZ(X)〞=A+ Z[2]cos(Z+O) +I◢
Y〝ZZ(Y)〞=B+Z[2]sin(Z+O)+J◢
G=Z+O◣◣◣◣◣◣↙
Pol(X－E〝CZ(X)=〞，Y－F〝CZ(Y)=〞) ：I〝ZZFY(JL)=〞◢J<0＝＞J = J +360: ≠＞J = J ◣J〝ZZFY(F**)°=〞◢
3.边桩坐标计算及放样程序 文件名 *Z ZB FY
{L，U}：Rec(L〝ZZ----BZ(JL=)〞，U〝JIA JIAO〞+G ) ：H〝BZ(X)〞=X+I◢W〝BZ(Y)〞=Y+J◢
Pol(H－E，W－F) ：I〝BZFY(JL)=〞◢J<0＝＞J = J +360: ≠＞J = J ◣
J〝BZFY(FWJ)°=〞◢↙
附注：QZ~YH段曲线内任意点的坐标用Z[6]〝HZ〞计算，将90代入括号
≠＞ K≤Z[5] ＝＞Pol((Rsin(180(Z[6] －K－N) ÷πR+90N÷πR)+T)，R(1－cos(180(Z[6] －K－N) ÷πR+90N÷πR))+Q) ：Rec(I，Z+O－180－(O÷AbsO )J) ：X〝ZZ(X)〞=A+ Z[2] cos(Z+O)+I◢ Y〝ZZ(Y)〞= B+ Z[2] sin(Z+O)+J◢G= Z+O－(O÷AbsO )(90N÷(πR)+180(Z[5] －K)÷πR) ：
程序符号说明
文件名FDCQX ZB FY的符号
R=? 输入曲线半径
PJ°=? 输入偏角（左负，右正）
Ls(1)=？ 输入第一曲线的缓和曲线长
Ls(2)=？ 输入第二曲线的缓和曲线长
T1= 输出第一曲线的切线长
T2= 输出第二曲线的切线长
LY= 输出圆曲线长
EH= 输出外距
D= 输出校正值
JD=? 输入交点桩号（例如：K20+200.123,输20200.123）
ZH= 输出第一回旋线起点
HY= 输出第一回旋线终点或圆曲线起点
QZ= 输出圆曲线中点
YH= 输出圆曲线终点或第二回旋线起点
HZ= 输出第二回旋线终点
JD(X)=? 输入交点的X 坐标
JD(Y)=? 输入交点的Y 坐标
JD(FWJ)°=? 输入本站交点与上一交点的方位角
输入1进入中桩坐标计算及放样程序（ZZ ZB FY），输入 2进入边桩坐标计算及放样程序（BZ ZB FY），输入其他数返回非对称曲线主点计算程序（FDCQX ZB FY）。
文件名 ZZ ZB FY的符号（G代表切线方位角）
K=? 输入需要加桩桩号（例如：K20+200.123,输20200.123）
I? 内部运算不考虑
ZZ(X)= 输出加桩点的X坐标
J? 内部运算不考虑
ZZ(Y)= 输出加桩点的 Y坐标
CZ(X)=? 输入测站点的X坐标
CZ(X)=? 输入测站点的Y坐标
ZZFY(JL)= 输出测站点距放样点的距离（中桩放样数值）
ZZFY(F**)°= 输出测站点距放样点的方位角（中桩放样数值）
文件名 *Z ZB FY的符号
ZZ----BZ(JL)=? 输入边桩距中桩的距离
JIA JIAO°=? 输入放样点与路线的夹角
（左负，右正，例如垂直则输入+90，－90）
BZ(X)= 输出边桩的X坐标
BZ(Y)= 输出边桩的Y坐标
BZFY(JL)= 输出测站点距放样点的距离（边桩放样数值）
BZFY(FWJ)°= 输出测站点距放样点的方位角（边桩放样数值）
修改时间:2005年8月26日星期五晚20：30于吴子高速公路N-18标驻地办