4800公路测量放样计算程序
CASIO fx-4800P计算器程序
DA DI——N E (公路弯道坐标计算)
K“JD”:R:S“LS”:F“PJ——L-R+”:A“ZH——N”:B“ZH——E”:C“JD——N”:D“JD——E”
M“ZH”=K-(R+S2÷(24R))tan(Abs F÷2)-S÷2+S^3÷(240R2)▲
N“HY”=M+S▲
O“QZ”=M+(πRAbs F÷180+S)÷2▲
Z[1]“YH”=M+πRAbs F÷180▲
Q“HZ”= Z[1]+S▲
T=K-M▲
L=Q-M▲
E=(R+S2÷(24R))÷cos(F÷2)-R▲
Lb1 0:{G}:G“ZHUANG HAO”:G<O==>Z=G-M≠=>Z=Q-G:△
{U}:U“I——B OUT+IN-”
Z<S+0.005==>H=tan-1(24Z2S2R2-Z^6)÷(48Z^3R^3-6Z^4SR)):V=Z-Z^5÷(40S2R2)+Usin H:W=Z^3÷(6SR)-Z ^7÷(336Z^3R^3)-Ucos H:≠=>H=180(Z-S÷2)÷(πR):V=Rsin H+S÷2-S^3÷(240R2)+Usin H:W=R-Rcos H+ S2÷(24R)-Ucos H:△
F<0==>G<O==>X=V:Y=W:≠=>X=T+Tcos F-Vcos F-Wsin Abs F:Y=Tsin Abs F-Vsin Abs F+Wcos F:△≠=>G<O==>X=V:Y=-W:≠=>X=T+Tcos F-Vcos F-Wsin Abs F:Y=Vsin Abs F-Tsin Abs F-Wcos F:△△
C=A==>D≥B==>P=90:≠=>P=270:△≠=>C<A==>P= tan-1((D-B)÷(C-A))+180:≠=>P= tan-1((D-B)÷(C-A)):△△
I“COM——N”=A+√(X2+Y2)cos(P-tan-1(Y÷X))▲
J“COM——E”=B+√(X2+Y2)sin(P-tan-1(Y÷X))▲
G=G+20
Goto 0
程序运行
输入:JD? ——交点桩号
R? ——圆曲线半径
LS? ——
缓和曲线长度
PJ——L-R+? ——偏角 右偏为正
ZH——N? ——ZH点N坐标
ZH——E? ——ZH点E坐标
JD——N? ——JD点N坐标
JD——E? ——JD点E坐标
输出:ZH= ——ZH点桩号
HY= ——HY点桩号
QZ= ——QZ点桩号
YH= ——YH点桩号
HZ= ——HZ点桩号
T= ——切线长
L= ——曲线长
E= ——外矢距
再输入:ZHUANG HAO? ——输入任一点桩号
I——B OUT+IN-? ——距中线距离(外侧为正)
再输出:COM——N= ——计算出N坐标
COM——E= ——计算出E坐标
以此循环——再输入——再输出
程序中各字母所代表的含义:
A——ZH点N坐标
B——ZH点E坐标
C——JD点N坐标
D——JD点E坐标
E——外矢距
F——偏角(右偏为正)
G——弯道内任一点桩号
H——弯道内任一点切线角
I——计算出的N坐标
J——计算出的E坐标
K——交点桩号
L——曲线长
M——ZH点桩号
N——HY点桩号
O——QZ点桩号
P——ZH--JD
方位角 Q——HZ点桩号
R——圆曲线半径
S——缓和曲线长度
T——切线长
U——距中线距离(外侧为正)
V——支距
W——支距
X——弯道内坐标
Y——弯道内坐标
Z——任一点曲线长度
Z[1]——YH点桩号
另个程序
公路逐桩坐标计算程序
(可以计算对称、不对称缓和曲线)
Lb1 0
Z=?
V=?
W=V+2:Fixm
{K}
Lb1 1
K>Z[W+5Z+4]=>W=W+1:Goto 1⊿ (判断桩号在哪个交点范围,就是该交点曲线起点至下一交点曲线起点)
S=K-Z[W+5Z+3] (计算该桩号与曲线起点的距离)
R=Z[W+2Z+2]:L=Z[W+3Z+2]:E=Z[W+4Z+2] (读取该交点曲线要素R、Ls1 、Ls2)
Pol(Z[W]-Z[W-1],Z[W+Z+2]-Z[W+Z+1]) (计算该交点与下一交点直线方位角)
J<0=>J=J+360⊿
A=J
Pol(Z[W-1]-Z[W-2],Z[W+Z+1]-Z[W+Z]) (计算该交点与上一交点直线方位角)
J<0=>J=J+360⊿
C=A-J:A=J (计算偏角)
W=V+2=>Goto2⊿ (如果桩号在起点与第一交点曲线起点之间,则转 Lb1 2 )
I=Abs(tan(c÷2))
M=L÷2-L^3÷240R^2:N=E÷2-E^3÷240R^2
P=L^2÷6R-L^4÷336R^3-R(1-cos(90L÷πR))
Q=E^2÷6R-E^4÷336R^3-R(1-cos(90E÷πR))
D=(P-Q)I÷2 : F=(P+Q+2R)I÷2
M=F+M-D:Q=F+N+D
N=πRAbsC÷180+(L+E)÷2
X=Z[W-1]-McosA
Y=Z[W+Z+1]-MsinA
M=Z[W-1]+Qcos(A+C)
V=Z[W+Z+1]+Qsin(A+C)
Q=AbsC÷C
S≤L=>P=0:Goto3⊿ (如果桩号在第一缓和曲线内,则转 Lb1 3)
S≤N-E=>S=S-L:Goto4⊿(如果桩号在圆曲线内,则转 Lb1 4)
S≤N=>S=N-S
Q=-Q:A=A+C-180:X=M:Y=V:L=E:P=180:Goto3 ⊿ (如果桩号在第二缓和曲线内,则转 Lb1 3)
P=A+C:S=S-N:D=M+ScosP:F=V+SsinP
Goto6 (如果桩号在直线内,则转 Lb1 6)
Lb1 2
P=A+C
D=Z[W-1]+ScosP
F=Z[W+Z+1]+SsinP:Goto6
Lb1 3
I=S-S^5÷40R^2÷L^2+S^9÷3456R^4÷L^4
J=Q(S^3÷6RL-S^7÷336R^3÷L^3)
P=P+A+90QS^2÷πRL:Goto5
Lb1 4
M=90(2S+L)÷πR
I=RsinM+L÷2-L^3÷240R^2
J=Q(L^2÷24R+R(1-cosM))
P=A+QM
Lb1 5
D=X+IcosA-jsinA:F=Y+JcosA+IsinA
Lb1 6
D〃X=〃◢ (结果显示X坐标)
F〃Y=〃◢ (结果显示Y坐标)
P〃AT=〃◢ (结果显示该桩号方位角)
{BO}:B〃S〃O〃⊿〃 (输入边桩距离,交角)
P=P+O
L〃XB〃=D+BcosP◢ (结果显示边桩X坐标)
M〃YB〃=F+BsinP◢ (结果显示边桩Y坐标)
以上是坐标计算程序,括号内是程序计算的大致原理及说明,中间部分为直线、圆曲线、缓和曲线计算的各种公式,大家也知道,书上也有。
该程序是最基本的,如再加几条语句就可以处理断链。
该程序还需变量输入,把整条线路曲线要素一次输入到4800中,就可以计算线路任意中、边桩坐标。