1#
用CASIO fx—4500计算器
计算道路平曲线中桩及其边桩的大地坐标
一：已知条件：
缓和曲线长：L1 ；平曲线半径：R ；导线边方位角：A ；交点里程桩号：J ；边桩横距：S ；直缓点、缓直点的大地坐标：X0、Y0。
二：计算程序：
ZBJS (程序名)
L1 Lbi 0
L2 E“L1”：R：A：B：J：S
L3 T=Rtan（Abs B/2）+E2/（24R）×tan
（Abs B/2）+E/2-E xY3/（240R2）▲
L4 L=Abs（RBπ/180）+Ｅ▲
L5 W=28.6479Ｅ/Ｒ：B<0 => W= -WΔ
L6 Z“ZH”=J-T▲H“HY”=Z+Ｅ▲
G“QZ”=Z+L/2▲D“YH”=Z+L-Ｅ▲
I“HZ”=Z+L▲
L7 Lbi 1:{N}:NGoto 0△
N>I=>Goto 0△
L8 N≤H =>F=N-Z:GotO 2△
L9 N≤D =>F=N-H:Goto 3△
L10 N >D =>F=I-N:Goto 2△
L11 Lbi 2:K=W×(F/E)2/3:C=F：Goto 4
L12 Lbi 3:K=90F/(πR)):
C＝2Rsin K:B<0 => K= -K△Goto 5
L13 Lbi 4:N>D = Goto 6△ U“X0”：V“Y0”
L14 X=U+C cos（A+K）▲Y=V+C sin（A+K）▲
L15 X[3]“X1”=X+S Cos（A+3K+90）▲
Y[3]“Y1”=Y+S Sin（A+3K+90）▲
L16 X[5]“X2”=X+S Cos（A+3K-90）▲
Y[5]“Y2”=Y+S Sin（A+3K-90）▲
L17 X[7]=X：Y[7]=Y：Goto 1
L18 Lbi 5：U=X[7]：V=Y[7]
L19 X = U+C Cos（A+W+K）▲
Y = V+C Sin（A+W+K）▲
L20 X[3]“X1”=X+S Cos（A+W+2K+90）▲
Y[3]“Y1”=Y+S Sin（A+W+2K+90）▲
L21 X[5]“X2”=X+S Cos（A+W+2K-90）▲
Y[5]“Y2”=Y+S Sin（A+W+2K-90）▲
Goto 1
L22 Lbi 6：P“X0”：Q“Y0”
L23 X= P+C Cos（A+B+180-K）▲
Y= Q+C Sin（A+B+180-K）▲
L24 X[3]“X1”=X+S Cos（A+B+90-3K）▲
Y[3]“Y1”=Y+S Sin（A+B+90-3K）▲
L25 X[5]“X2”=X+S Cos（A+B-90-3K）▲
Y[5]“Y2”=X+S Cos（A+B-90-3K）▲
Goto 1
三、符号标识：
L1：缓和曲线；R：平曲线半径；A：导线边方位角；B：转角；J：交点里程；S：边桩横距；T：切线长；L：曲线长；ZH：直缓点里程；HY：缓圆点里程；QZ：曲中点里程；YH：圆缓点里程；HZ：缓直点里程；N：曲线上任一点里程桩号；X0、Y0：直缓点、缓
直点大地坐标；X、Y曲线上任一点大地坐标；
X1、Y1、X2、Y2：对应点边桩大地坐标。
四、算例：
某弯道平曲线缓和曲线长200m，半径1800m，起始导线边方位角87о09′06″，左转角15ο26′20″，交点里程102K+044.126，直缓点大地坐标：（38310.860,40822.807）,缓直点大地坐标(38435.931,41493.211)。边桩横距30m。
运行本程序，根据计算器提示逐一输入数据，可显示：
1:曲线元素：切线长T=344.107m;
曲线长L=685.027m
2:曲线主点里程:ZH=101+700.019;
HY=101+900.019;QZ=102+042.533;
YH=102+185.046;HZ=102+385.046
3、中桩及边桩的大地坐标：
里程桩号 中桩点的大地坐标 边桩点的大地坐标 备注
X Y X1 Y1 X2 Y2
101+700.019 38310.860 40822.807 38280.897 40824.298 38340.823 40821.316 ZH点
+740 38312.876 40862.737 38282.917 40864.294 38342.836 40861.180
+780 38315.071 40902.677 38285.122 40904.434 38345.019 40900.920
+800 38316.290 40922.640 38286.351 40924.547 38346.230 40920.734
+840 38319.084 40962.546 38289.173 40964.852 38348.995 40960.240
+880 38322.498 41002.411 38292.633 41005.248 38352.364 40999.574
101+900.019 38324.496 41022.342 38294.662 41025.494 38354.330 41019.189 HY点
+940 38329.138 41062.051 38299.381 41065.865 38358.894 41058.237
+980 38334.664 41101.667 38304.999 41106.141 38364.328 41097.193
102+000 38337.756 41121.426 38308.143 41126.230 38367.369 41116.623
102+042.533 38345.062 41163.326 38315.571 41168.828 38374.553 41157.824 QZ点
+080 38352.316 41200.083 38322.946 41206.198 38381.687 41193.969
102+100 38356.501 41219.641 38327.201 41226.081 38385.802 41213.200
+140 38365.552 41258.609 38336.372 41265.699 38394.672 41251.519
+160 38370.357 41278.016 38341.287 41285.429 38399.426 41270.603
102+185.046 38376.715 41302.242 38347.751 41310.059 38405.678 41294.424 YH点
+200 38380.660 41316.612 38351.759 41324.661 38409.560 41308.564
+260 38397.554 41374.199 38368.872 41382.992 38426.237 41365.407
+300 38409.516 41412.371 38380.938 41421.498 38438.094 41403.244
+340 38421.837 41450.427 38393.325 41459.760 38450.348 41441.095
+360 38428.079 41469.427 38399.586 41478.816 38456.572 41460.039
102+385.046 38435.931 41493.211 38407.446 41502.624 38464.416 41483.798 HZ点
五、说明：
1：运行程序时，先扩充储存区：MODE Ans 4 EXE
2：当平曲线为左转弯时，转角B用其负值输入。
3：计算主曲线的点的坐标时，请务必先在前一桩号计算出“HY点”的坐标。
原作者为：hfsh