1、CASIO 计算器竖曲线计算程式R “ B J “ T “ Q X “ Q “ Q Z L “H “ Q Z H “ C “ C P “ E “ H P “: A = Q - T : B = Q + T Lbl 0 : I : I “ S Q “ I B G = H + E ( I - Q ) = S = I - Q : I D = B - I : P = S E : Goto 1 Lbl 1 : C 0 Goto 2 = Goto 3 Lbl 2 : E 0 Goto 4 = Goto 6 Lbl 3 : E 0 Goto 7 = Goto 5 Lbl 4 : C E Goto 6 = Go
2、to 7 Lbl 5 : C E Goto 7 = Goto 6 Lbl 6 : G “ G H “ = H + P - D2 2 R Goto 0 Lbl 7 : G “ G H “ = H + P + D2 2 R Goto 0 BJ 半径 QX 切线长 QZL 交点里程QZH 交点高程 CP 前坡 HP 后坡SQ 所求里程公路三维坐标控制系统高程控制部分(for 4850)公路三维坐标控制系统高程控制部分(for CASIO4850)主程序 DHSZ:SZSS1N=1:Prog“HDAT”:GOTO 1SS2N=2:Prog“HDAT”:GOTO 1SS3N=3:Prog“HDAT”:G
3、OTO 1Lbl 1:S I Q= F+G(S- E S I Q= F+G(S- E)+(-1)J(S-I)2/(2 MR= Q+11.75 P-0.03:”DH=”: Q= Q Z=0Q= Q Z0Q= R H:H”H=”: T= Q-H“ H=”: T= T ABKL: T=Abs T: U= T / 8 : T 8W= TA+14+LT 8 W=8A+B(T-8)+14+2Int U+LT 16 W=8A+8B+K(T-16)+14+2 Int U+L“DW”: W= W “ W”: D= W-Abs Z竖曲线数据库 HDATN=1E= xxxx.xxx: F= xxxx.xxx: G=
4、 xxxx.xxx: J= xxxx.xxx: M= xxxx.xxx:O= xxxx.xxxN=2E= xxxx.xxx: F= xxxx.xxx: G= xxxx.xxx: J= xxxx.xxx: M= xxxx.xxx: O= xxxx.xxx N=3E= xxxx.xxx: F= xxxx.xxx: G= xxxx.xxx: J= xxxx.xxx: M= xxxx.xxx: O= xxxx.xxxI= E- O说明：(1)变量赋值说明E 变坡点桩号F 变坡点高程G 变坡点前段坡比I 该段竖曲线起点桩号即高差改正值为零处桩号J 竖曲线系数 1 为凸 2 为凹M 竖曲线半径O 竖曲线切
5、线长P 路面横坡（已在坐标反算程序中赋值）Q 路线设计高程，当 Z=0 时为中桩设计高程（或中分带边缘设计高程），当 Z0时为路基边缘设计高程R 路基边缘设计高程T：“H”高差V 边坡分级，第 n 级坡W 开口、坡脚设计宽度（根据高差及坡比求得）可计算任意高度处的设计宽度D 宽差（当其为负則向中桩方向移相应宽度，当其为正则向外移相应宽度）S、Z 所测点对应桩号及距中桩处的距离A、B、K 横断面一、二、三级坡比L 附加宽度（以便根据横断面设计实际情况调整，象坠落台、护坡道、路堑处边沟设置宽度的不固定值设置及填方段为使路基满足设计的压实要求超填的宽度）H 实测高程(2)编制思路该程序编者根据工地的
6、实际情况编制在此仅供一个编制思路，用户需根据实际情况进行修改。注意：程序中有下划线处数据说明（单位：cm）11.75 为中央分隔带边缘距硬路肩宽度，该段横坡为一定值（超高段需加入超高计算子程序，另外在设计宽度计算时加入加宽子程序即可）003 为土路肩内外侧高差，因土路肩与行车道及硬路肩横坡不同，故分开计算。（用户需根据实际情况修改）2 为平台设置宽度，第 n 级坡的平台设置总宽度即为 2n. （用户需根据实际情况修改）8 每级边坡设置高度.14 路面设计宽度（路线中线处距土路肩边缘设计宽度）该程序为公路三维控制系统中的一部分，即高程控制部分。该程序可与 hanniu_1973老师的曲线任意里程
7、中边桩坐标正反算(CASIO fx-4850P)程序（第三次修改版）及大部分坐标正反算程序联合使用（修改变量名后）。曲线任意里程中边桩坐标正反算(CASIO fx-4850P 计算器)程序（修改版）曲线任意里程中边桩坐标正反算(CASIO fx-4850P 计算器)程序（修改版）说明：1. 在此特别感谢佚名对本人原程序的扩展，并在“工程测量员” 网站中发布。2.程序中存有下主要问题：（1）Pol()数得到角度不会大于 360 度，“:F360=F=F-36 ”是多余的；（2）“FDMS”在程序中是不可能实现的；（3）“F=1-L:M=1-K”是多余的，造成 M0.5；（4）正算时，当 Z XY
8、“：“2.XY = SZ“：N：Prog”DAT1”:U“X0“：V“Y0“：O“S0“：G“F0“：H“LS“：P“R0“：R“RN“ ：Q：C=1P：D=(P-R)(2HPR)：E=180：N=1=Goto 1：Goto 2Lbl 1：SZT：SZ:T”ANG” ： W=Abs(S-O)：Prog “SUB1“：“XS=“:X=X“YS=“:Y=YGoto 3Lbl 2：XY ：XY：I=X：J=Y ：Prog “SUB2“：“S=“:S=O+W“Z=“:Z=ZGoto 2Lbl 3:J=0:I=Pol(X-Z1,Y-Z2):F=J:FF=F+360:F360=F=F-36 ”DIST=”
9、: I”FW=”:FDMS Goto 12. 正算子程序 (SUB1)A=0.1184634425:B=0.2393143352:N=0.2844444444:K=0.0469100770:L=0.2307653449:M=0.5:F=1-L:M=1-K:X=U+W(Acos(G+QEKW(C+KWD)+Bcos(G+QELW(C+LWD)+Ncos(G+QEMW(C+MWD)+Bcos(G+QE(1-L)W(C+(1-L)WD)+Acos(G+QE(1-K)W(C+(1-K)WD): Y=V+W(Asin(G+QEKW(C+KWD)+Bsin(G+QELW(C+LWD)+Nsin(G+QEM
10、W(C+MWD)+Bsin (G+QE(1-L)W(C+(1-L)WD)+Asin (G+QE(1-K)W(C+(1-K)WD)：F=G+QEW(C+WD)+T：X=X+ZcosF：Y=Y+ZsinF3. 反算子程序 (SUB2)T=G-90：W=Abs(Y-V)cosT-(X-U)sinT)：Z=0：Lbl 0：Prog “SUB1“：L=T+QEW(C+WD)：Z=(J-Y)cosL-(I-X)sinL：AbsZGoto1：W=W+Z：Goto 0Lbl 1：Z=0：Prog “SUB1“：Z=(J-Y)sinF4.曲线元要素数据库：DAT1J=1=U=xxxx.xxx:V=xxxx.xx
11、x:O=xxxx.xxx:G=xxxx.xxx:H=xxxx.xxx:P=xxxx.xxx:R=xxxx.xxx:Q=xxJ=2=U=xxxx.xxx:V=xxxx.xxx:O=xxxx.xxx:G=xxxx.xxx:H=xxxx.xxx:P=xxxx.xxx:R=xxxx.xxx:Q=xx.J=n=U. (注：如有多个曲线元要素继续添加入数据库 DAT1 中)三、使用说明1、规定(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时， Q=1；当线元为直线时，Q=0。(2) 当所求点位于中线时， Z=0；当位于中线左铡时，Z 取负值；当位于中线中线
12、右侧时，Z 取正值。(3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以 10 的 45 次代替。(4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。(5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以 10 的 45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以 10 的 45 次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。(6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆
13、曲线相接时，曲率半径等于圆曲线的半径。(7)曲线元要素数据库（ DAT1）可根据线型不同分为各个线元段输入到 DAT1 中，即分为直线段、缓和曲线、圆曲线等。2、输入与显示说明输入部分：1. SZ = XY2. XY = SZ N ? 选择计算方式，输入 1 表示进行由里程、边距计算坐标 ；输入 2 表示由坐标反算里程和边距。X0 ？线元起点的 X 坐标 Y0 ？线元起点的 Y 坐标 S0 ？线元起点里程F0 ？线元起点切线方位角LS ？线元长度R0 ？线元起点曲率半径RN ？线元止点曲率半径Q ？ 线 元左右偏标志 (左偏 Q=-1，右偏 Q=1，直线段 Q=0)S ？ 正算时所求点的里程Z
14、 ？正算时所求点距中线的边距(左侧取负，值右侧取正值，在中线上取零)ANG？正算边桩时左右边桩连线与线路中线的右交角J？曲线元数据库曲线段判断系数（J=1、2n）X ？反算时所求点的 X 坐标Y ？反算时所求点的 Y 坐标显示部分：XS= 正算时，计算得出的所求点的 X 坐标YS= 正算时，计算得出的所求点的 Y 坐标S= 反算时，计算得出的所求点的里程Z= 反算时，计算得出的所求点的边距四、算例某匝道的由五段线元（直线+完整缓和曲线+圆曲线+ 非完整缓和曲线 +直线）组成，各段线元的要素（起点里程 S0、起点坐标 X0 Y0、起点切线方位角 F0、线元长度LS、起点曲率半径 R0、止点曲率半
15、径 RN、线 元左右偏标志 Q）如下：S0 X0 Y0 F0 LS R0 RN Q500.000 19942.837 28343.561 125 16 31.00 269.256 1E45 1E45 0769.256 19787.340 28563.378 125 16 31.00 37.492 1E45 221.75 -1806.748 19766.566 28594.574 120 25 54.07 112.779 221.75 221.75 -1919.527 19736.072 28701.893 91 17 30.63 80.285 221.75 9579.228 -1999.812
16、 19744.038 28781.659 80 40 50.00 100.000 1E45 1E45 0 1、正算（注意：略去计算方式及线元要素输入，请自行根据所求点所在的线元输入线元要素）S=700 Z=-5 计算得 XS=19831.41785 YS=28509.72590S=700 Z=0 计算得 XS=19827.33592 YS=28506.83837S=700 Z= 5 计算得 XS=19823.25398 YS=28503.95084S=780 Z=-5 计算得 XS=19785.25749 YS=28575.02270S=780 Z=0 计算得 XS=19781.15561 Y
17、S=28572.16358S=780 Z= 5 计算得 XS=19777.05373 YS=28569.30446S=870 Z=-5 计算得 XS=19747.53609 YS=28654.13091S=870 Z=0 计算得 XS=19742.68648 YS=28652.91379S=870 Z= 5 计算得 XS=19737.83688 YS=28651.69668S=940 Z=-5.123 计算得 XS=19741. 59118 YS=28722.05802S=940 Z=0 计算得 XS=19736.47687 YS=28722.35642S=940 Z= 3.009 计算得 X
18、S=19733.47298 YS=28722.531682、 反算X=19831.418 Y=28509.726 计算得 S=699.9999974 Z= -5 .00018164X=19827.336 Y=28506.838 计算得 S=699.9996493 Z= 0.000145136X=19823.25398 Y=28503.95084 计算得 S=699.9999985 Z= 5.000003137X=19785.25749 Y=28575.02270 计算得 S=780.0000035 Z= -5 .000001663X=19781.15561 Y=28572.16358 计算得
19、S=780.0000025 Z=- 0.000002979X=19777.05373 Y=28569.30446 计算得 S=780.0000016 Z= 4.99999578X=19747.536 Y=28654.131 计算得 S=870.0001137 Z= -4.99941049X=19742.686 Y=28652.914 计算得 S=870.0003175 Z=- 0.00041814X=19737.837 Y=28651.697 计算得 S=870.0002748 Z= 4.999808656X=19741.5912 Y=28722.0580 计算得 S=939.9999786 Z= -5.123024937X=19736.4769 Y=28722.3564 计算得 S=939.9999862 Z=- 0.000027710X=19733.4730 Y=28722.5317 计算得 S=940.0000238 Z= 3.00898694
