公路中边桩计算程序
35→DIMZ
Fix 4
LBI 0
“K[JD]=”?B:“X[JD]=”?C:“Y[JD]=”?D:“LS1=”?G:“LS2=”?H: “R=”?R:“F1=”?W: “F2=”?K: “LY=”?P: “A=”?A:“O=”?O:
Prog“YS”
LBI 1
“K=”?M
IF M ≤Z[8]:Then Goto 2:ifend
IF M <Z[9] And M>Z[8]:Then Goto 3:ifend
IF M≥Z[9] :Then Goto 4:ifend
LBI 2
(M-Z[7])-(M-Z[7])^(5)÷(40×R2×G2)→Z[21]
(M-Z[7]) ^(3)÷(6×R×G)-(M-Z[7]) ^(7)÷(336×R^(3)G^(3)) →Z[22]
Z[12]+Z[21]COS(W)-O×Z[22]SIN(W)→N
Z[13]+Z[21]SIN(W)+O×Z[22]COS(W)→E
“B=”?V
“X=”:N+VCOS(W+O×90×(M-Z[7])2÷(πRG)+90) ◢
“Y=”:E+VSIN(W+O×90×(M-Z[7])2÷(πRG)+90) ◢
GOTO 1
LBI 3
90×(M-Z[8])÷(πR)→Z[21]
2×RSIN(Z[21]) →Z[22]
Z[20]+O×Z[21] →Z[23]
Z[16]+Z[22]COS(Z[23])→N
Z[17]+Z[22]SIN(Z[23])→E
“B=”?V
“X=”:N-VCOS(Z[20]+O×2×Z[21]-90)◢
“Y=”:E-VSIN(Z[20]+O×2×Z[21]-90)◢
GOTO 1
LBI 4
(Z[10]-M)-(Z[10] -M)^(5)÷(40R2H2)→Z[21]
(Z[10]-M)^(3)÷(6RH)-(Z[10]-M)^(7)÷(336R^(3)H^(3))→Z[22]
Z[18]+Z[21]COS(K)+O×Z[22]SIN(K)→N
Z[19]+Z[21]SIN(K)-O×Z[22]COS(K)→E
“B=”?V
“X=”:N+VCOS(K-O×90×(Z[10]-M)2÷(πRH)-90) ◢
“Y=”:E+VSIN(K-O×90×(Z[10]-M)2÷(πRH)-90) ◢
GOTO1
YS(子程序)
G÷2-G^(3)÷240÷R2+G^(5)÷34560÷R^(4) →Z[1]
H÷2-H^(3)÷240÷R2+H^(5)÷34560÷R^(4) →Z[2]
G2÷24÷R-G^(4)÷2688÷R^(3) →Z[3]
H2÷24÷R-H^(4)÷2688÷R^(3) →Z[4]
(Z[4]-Z[3])÷(2TAN(A÷2))+0.5×(Z[3]+Z[4]+2R)TAN(A÷2)+Z[1] →Z[5]
(Z[3]-Z[4])÷(2TAN(A÷2))+0.5×(Z[3]+Z[4]+2R)TAN(A÷2)+Z[2] →Z[6]
B-Z[5] →Z[7]
Z[7]+G→Z[8]
Z[8]+P→Z[9]
Z[9]+H→Z[10]
Z[10]-(G+H+P)÷2→Z[11]
C-Z[5]COS(W) →Z[12]
D-Z[5]SIN(W) →Z[13]
G-G^ (3)÷(40R2)→Z[14]
G2÷6÷R-G^ (4)÷336÷R^ (3) →Z[15]
Z[12]+Z[14]COS(W)-O×Z[15]SIN(W)→Z[16]
Z[13]+Z[14]SIN(W)+O×Z[15]COS(W)→Z[17]
C-Z[16]COS(K)→Z[18]
D-Z[16]SIN(K)→Z[19]
W+O×90G÷π÷R→Z[20]
Z[7] →List X[1]
Z[12] →List Y[1]
Z[13] →List freq[1]
Z[8] →List X[2]
Z[16] →List Y[2]
Z[17] →List freq[2]
Z[10] →List X[3]
Z[18] →List Y[3]
Z[19] →List freq[3]
K[JD] 交点桩号
X[JD] 交点X坐标
Y[JD] 交点Y坐标
Ls1 第一缓和曲线
Ls2 第二缓和曲线
LY 圆曲线
R 半径
F1 ZH→JD方位角
F2 HZ→JD方位角
A 转角
O 转向系数(左=-1右=1)
K 桩号
B 路半幅宽度(右边为正,左边为负。当B=0时计算中桩)