缓和曲线高程计算程序.doc

直竖曲线高程计算程序清单一、计算依据中桩计算时依据：变坡点里程桩号和高程，相邻坡段纵坡度，竖曲线半径；计算边桩时还需要中桩至边桩距离和路拱（横坡）.二、计算范围 计算竖曲线上各点高程时，只能在竖曲线范围内计算，竖曲线外则是直线或平曲线.三、计算竖曲线高程放样数据计算公及说明：1） 中桩高程计算公式及说明：Y=H+QxI+Q（T-x）^2/2/R 式中：H——变坡点高程（m）；x——变坡点到两侧坡段之纵坡，上坡为正，下坡为负；T——竖曲线切线长度（m）；R——竖曲线半径（m）；Q——控制竖曲线凸凹，凹取 1，凸取 -12） 边桩高程计算公式及说明：Z=y+BEB——半副路宽E——路拱3） 直竖结合计算线路高程计算公式及说明：公式：G=H-CP+ZF（T-AbsC）^2/2/R G——线路直线、曲线（平、竖、超高）上任意一点中桩的设计高程；H——变坡点高程；C=B-L ，B 为变坡点里程桩号；L——为线路上任意一点（所求点）里程桩号；P、F、Z——控制前纵坡度 I、后纵坡度 J 的条件；T——竖曲线切线长度；R——竖曲线半径四、直竖结合曲线高程计算程序清单：LbI O ：H ：B ： R：I ：J：{LNM}T=R Abs（J-I） /2C=B-LF=1I>J=>F=-1⊿L≤ 0 =>{HRBIJ}：Goto 0≠=>LZ=0：P=I≠=>LZ=1：P=I≠=>LZ=1：P=J≠=> Z=0：P=J⊿⊿⊿G=H-N-CP+ZF（T- AbsC）^2/2/R◢V=G+ME◢Goto 0程序中：N——路面层至各施工层的厚度（M） M——中桩至边桩的距离E——路拱 （横坡）V——左（右）边桩高程 其他字符含义请参照上面公式中的注释 使用说明1. 计算范围：前竖曲线终点至后竖曲线起点（直线--本竖曲线—直线），只要在这段范围内路线的直线、圆曲线、缓和曲线超高、竖曲线上任意一点的中桩设计高程都可以计算；边桩高程除超高段需另行计算外，其他直线、圆曲线、竖曲线都可以计算 。

2. 计算时以变坡点为起点，所需计算要素有该变坡点里程桩号相邻坡段的前纵坡度 I、后纵坡度 J 和变坡点所在竖曲线的半径 R3. L 为计算范围内任意一点里程桩号，计算过程中，只要输入 L 的里程桩就可以计算所需点的中桩高程此程序的三大功能特点：1. 当 L 输入 0 时，计算自动中止，需要重新输入起点要素（H、B、R、I、J）. 这一功能，可以帮助检查输入起算要素的正确性，或是进行下一个竖曲线计算时，不需再重新寻找文件名，方便操作.2. 在计算中桩高程的同时，可以计算与中桩同一横断面的左右边桩高程. 3. 在计算中桩高程的同时，可以计算出路面各结构层的中桩、边桩高程.特别注意：对于缓和曲线超高段，此程序只能计算出中桩高程，左右边桩高程需要另行计算（续见缓和曲线超高段中边桩高程放样计算程序）.实践证明，直竖结合曲线程序计算线路设计高程的方法在公路施工中非常适用.超高渐变率 P 按下式计算:P=B×(I2%-I1%)/Lc其中 B 为超高旋转轴至路缘带外缘的距离Lc 为超高渐变段长度I1、I2 为超高段起终点横坡根据超高方式图中给定的相应参数就可以计算各桩号横坡CASIO4500 坐标计算程序默认分类 2008-07-27 10:49:06 阅读 1388 评论 2 字号：大中小 订阅 CASIO4500 坐标计算程序竖曲线计算程序 程序名：SQX LBl 1：{DEIHRQPJ}：W=DEIHRQPJ：T=(I-H)R/2：N=T2/2R LBl 2：{A}：Z=A LBl 3：A>P=>Goto 4：△A>E=>Goto 6：△A>Q=>Goto 5：△X=D+(A-E)I▲ Goto 2 LBl 4：X= D+(A-E)H▲Goto 2 LBl 5：M=(Q-A)2/2RJ：X=D+(A-E)I+M▲ Goto 2 LBl 6：M=(P-A)2/2RJ：X=D+(A-E)H+M▲ Goto 2 注意：本程序在 CASIO4500 上调试运行通过。

如果是在4800 上运行请将源程序中的 V，W 分别用 I，J 替代即可 字母含义： D：变坡点高程 E：变坡点里程 J：曲线凸凹（凸为-1、凹为+1） L：开始坡度 H：结束坡度 A：待求点标高 R：竖曲线半径 Q：起点里程 P：终点里程==== 4 ====计算曲线上任一点放样极坐标 F1： QXJS L1 Abs W"FixmC=>"=7=>Fixm:Goto F△Fix 4:Z=E7A"A0"/Abs AN[9:A=Abs A:Q=R+S"Lsci"2/24R L2 T=S/2-Sxy3/240R2+Anstan .5A▲L=πRA/180 +S▲H"ZH"=J"JDT"-T▲V"HY"=H+S▲V"QZ"=H+.5L▲V"YH"=H+L-S▲L"HZ"=H+L▲D"xJ"K"yJ" B"Gra"O"XE"U"YE"G"Xm"M"Xm" L3 LbI F:{FN}:F≤H=>X=F-J:ZN:≠>F≤L-S=>I= L-S=>I= F-H:FProg Q:≠>Rec(R-ZN,90(2I-S)/πR: X=W-Qtan.5A:Q-V△≠>I=L-F:PoI(T-I,ZN:FProg Q:PoI(-X,Ans△X=Rec(Ans,A+W L4 W△△X=D+Rec(PoI(X,ZAns,W+B:Y=K+W: "EnD OK!":PoI(G-O,M-U:I=W:PoI(X-O,Y-U:W-I: AnsAns+360△I=Ans:60Frac Ans:I"≤Deg"=IntI+.01Int Ans+6mFrac Ans▲V"P"=V▲Goto F [主程序结束] F2： Q [子程序] L1 Y=2RS:W=180I2/πY:X=I-Ixy5/10Y2-T-ZNsin W:Ixy3/3Y-Ixy7/42Yxy3+ZNcosW [子程序结束] 输入说明： F1 为第一个程序，L1 为第一个程序的第一行，L2 为第二行，依次类推。

1、程序中 Z=E7A"A0"/Abs AN[9 处 E 用按键[EXP] 键输入；2、程序中 π 代表圆周率，也就是 3.1415926 的那个； 3、程序中 Ans 用按键[Ans]输入； 4、程序中 I"≦Deg"=Int I+.01Int Ans+6mFrac Ans▲处的m 用按键[SHIFT][5]输入； 操作说明： Fixm C=>? 首次启动时，需输入密码：“10000000”重启动功能控制符，输入“7”或“-7”可启动重启动功能 A0？ 曲线转角，要求以度分秒形式输入，左转输入“负” 值，右转输入“正”值 R? 曲线半径，直接输入即可 LSci? 缓和曲线长度，有缓和曲线时输入缓和曲线长度，圆曲线输入“0” T=S/2-Sxy3/2→ 程序输出曲线切线长度 L=ПRA/180+S 程序输出曲线长度 JDT? 曲线交点里程，直接输入平曲线交点里程，如K123+456.789 输入 123456.789 即可 ZH= 程序输出曲线起点[ZH]里程 HY= 程序输出 HY 里程 HY= 程序输出 HY 里程 QZ= 程序输出 QZ 里程 YH= 程序输出 YH 里程。

HZ= 程序输出曲线止点[HZ]里程 XJ? 交点大地坐标[X 坐标] YJ? 交点大地坐标[Y 坐标] Gra? 曲线起始方位角即由 ZH 指向实交点的方位角 XE? 输入放样置仪点大地坐标[X 坐标] YE? 输入放样置仪点大地坐标[Y 坐标] Xm? 输入放样后视点大地坐标[X 坐标] Ym? 输入放样后视点大地坐标[Y 坐标] F? 放样点里程桩号 N？ 放样中桩距，即边桩至中桩的距离无论曲线左转还是右转，均以左“负”右“正”输入，中桩输入“0” EnD OK! 程序提示：程序已完成放样点的坐标计算，正在进行放样元素计算，请稍候！ ≤Deg= 程序输出放样极角为 28 度 16 分 47 秒(度分秒化处理)，不要再人为转换，将仪器水平度盘拨至该角度即可放样 P= 程序输出放样极距即置仪点至放样点的水平距离按 EXE键循环到 F? 5 ====计算曲线上任一点大地坐标 F1： QXJS L1 Abs W"FixmC=>"=7=>Fixm:Goto F△Z=E7A"A0"/Abs AN[9:A=Abs A:Q=R+S"Lsci"2/24R L2 T=S/2-Sxy3/240R2+Anstan .5A▲L=πRA/180+S▲H"ZH"=J"JDT"-T▲V"HY"=H+S▲V"QZ"=H+.5L▲V"YH"=H+L-S▲L"HZ"=H+L▲D"xJ"K"yJ" B"Gra" L3 LbI F:{FN}:F≤H=>X=F-J:ZN:≠>F≤L-S=>I=F-H:FProg Q:≠>Rec(R-ZN,90(2I-S)/πR: X=W-Qtan.5A:Q-V△≠>I=L-F:PoI(T-I,ZN:FProg Q:PoI(-X,Ans△X=Rec(Ans,A+W L4 W△△X”X”=D+Rec(PoI(X,ZAns,W+B▲Y”Y”=K+W▲ Goto F [主程序结束] F2： Q [子程序] L1 Y=2RS:W=180I2/πY:X=I-Ixy5/10Y2-T-ZNsin W:Ixy3/3Y-Ixy7/42Yxy3+ZNcosW [子程序结束] 输入说明： F1 为第一个程序，L1 为第一个程序的第一行，L2 为第二行，依次类推。

1、程序中 Z=E7A"A0"/Abs AN[9 处 E 用按键[EXP] 键输入；2、程序中 π 代表圆周率，也就是 3.1415926 的那个； 3、程序中 Ans 用按键[Ans]输入； 4、程序中 I"≦Deg"=Int I+.01Int Ans+6mFrac Ans▲处的m 用按键[SHIFT][5]输入； 操作说明： Fixm C=>? 首次启动时，需输入密码：“10000000”重启动功能控制符，输入“7”或“-7”可启动重启动功能 A0？ 曲线转角，要求以度分秒形式输入，左转输入“负” 值，右转输入“正”值 R? 曲线半径，直接输入即可 LSci? 缓和曲线长度，有缓和曲线时输入缓和曲线长度，圆曲线输入“0”T=S/2-Sxy3/2→ 程序输出曲线切线长度 L=ПRA/180+S 程序输出曲线长度 JDT? 曲线交点里程，直接输入平曲线交点里程，如K123+456.789 输入 123456.789 即可 ZH= 程序输出曲线起点[ZH]里程 HY= 程序输出 HY 里程 HY= 程序输出 HY 里程 QZ= 程序输出 QZ 里程 YH= 程序输出 YH 里程 HZ= 程序输出曲线止点[HZ]里程。

XJ? 交点大地坐标[X 坐标] YJ? 交点大地坐标[Y 坐标] Gra? 曲线起始方位角即由 ZH 指向实交点的方位角 F? 放样点里程桩号 N？ 放样中桩距，即边桩至中桩的距离无论曲线左转还是右转，均以左“负”右“正”输入，中桩输入“0” X= 程序输出放样点坐标[X 坐标] Y= 程序输出放样点坐标[Y 坐标]按 EXE 键循环到 F?三、结束语 以上程序可以输入 Fx-4500PA 计算器内，只要程序容量允许以上程序均用于实际计算工作当中，结果可靠，操作简单这些程序在 CASIO4500 上调试运行通过如果是在 4800 上运行请将源程序中的 V，W 分别用 I，J 替代即可此外，该计算器还可进行一些其它的测量内外业计算，因篇幅关系仅将上述程序做以介。