5800公路测量程序(修改).doc

CASIO-5800P计算器以其轻巧灵活、便于携带、功能强大、经济实用的特点，深受广大测量爱好者的喜爱。尤其在公路施工测量方面，全站仪与CASIO-5800P计算器的组合是目前公路施工测量应用最广泛的工具。为使其功能得到更好的应用，特编辑整理此程序，供广大公路测量爱好者参考。不妥之处敬请指教！ 编者：管东利 2009-3-15简 要 说 明本程序具有以下特点：1、 简洁实用，结构合理，功能齐全。由项目程序、主程序和子程序构成本程序系统。2、 有功能强大的路线坐标计算程序，同时配备了其它常用的应用程序。如：后方交会、面积计算、点到直线的距离及两直线交点坐标、距离方位反算等。3、 路线坐标计算：可以将整个项目的多条路线的多条曲线要素一次编入子程序(P)；竖曲线要素编入子程序(G),横坡度编入子程序（I）。做到一次编辑，项目程序(如：YL2)多次运行，直到本项目结束，避免每次计算都必须输入曲线要素的重复劳动，也减少了计算出错机会，极大地提高工作效率。这一点也是本程序的精华所在。4、 可对每条曲线进行坐标计算：即可使用曲线要素法计算，也可使用主点桩坐标法进行计算。其中主点桩坐标法可以用来计算任何曲线，尤其适用于对付互通立交。5、 与设计文件紧密结合。几乎所有已知曲线要素都可以在设计文件中直接找到，避免了重新计算的麻烦。6、 变量名力求与设计、习惯或数学公式相匹配，易于联想记忆。如：半径-R、直缓点-ZH、水平角-HR、里程桩号-LX等。7、 功能齐全的子程序，使编制新的应用程序简便易行。CASIO-5800P公路施工测量应用程序(GDL-2009)编制：管东利YL2(程序名：永蓝2段)Lb11:Prog“1”:Prog“P”:Prog“2”: Goto1输入：LX : 计算点桩号（对应变量V）; S: 计算点距该中桩距离（对应变量S）; W: 计算点与中桩连线与路线前进方向的夹角（对应变量W）计算并显示：计算点坐标X ,Y 说明:本程序为项目程序，先将本项目的平曲线要素编辑到子程序P，一次输入，长期使用。YLGC(永蓝高程)Lb11:Prog“1”:Prog“G”: Prog“I”: “H=”:H+I(S-1H Goto1输入：LX : 计算点桩号（对应变量V）; S: 计算点距该中桩距离（对应变量S）; W: 计算点与中桩连线与路线前进方向的夹角（对应变量W）计算并显示： 计算点对应的路面设计高程H说明:本程序为项目程序，先将本项目的竖曲线要素分别编辑到子程序G和子程序I中，一次输入，长期使用。(注：S要大于或等于零；W为900；1为高程设计线到路中线的距离，一般为中央分隔带宽度的一半，要与设计图纸相符。)ZX(直线段)Prog“8”:Lb1 1:Prog“1”:Prog“Z”:Prog“2”:Goto1已知：LA: 第一点桩号（对应变量K）; X1: 第一点X坐标（对应变量A）; Y1: 第一点Y坐标（对应变量B）;LN: 第二点桩号（对应变量N）; X2: 第二点X坐标（对应变量C）; Y2: 第二点Y坐标（对应变量D）输入：LX: 计算点桩号（对应变量V）; S: 计算点距该中桩距离（对应变量S）; W: 计算点与中桩连线与路线前进方向的夹角（对应变量W）计算并显示：计算点坐标X ,Y说明： 对应变量，在编辑子程序“P”和“G”时用的到，要对照输入相应的值。 YQX(圆曲线)Prog“8”:?R:Prog“7”:Lb11:Prog“1”:Prog“Y”:Prog“2”:Goto 1已知：LA: 第一点桩号（对应变量K）; X1: 第一点X坐标（对应变量A）; Y1: 第一点Y坐标（对应变量B）;LN: 第二点桩号（对应变量N）; X2: 第二点X坐标（对应变量C）; Y2: 第二点Y坐标（对应变量D）； R: 圆曲线半径（对应变量R）; Q: 曲线转向(Y=1右转 Z=-1左转) （对应变量Q） 输入：LX: 计算点桩号; S: 计算点距中桩距离; W: 中桩与计算点连线与路线前进方向夹角计算并显示：计算点坐标X ,Y示意图如下：HHQX(缓和曲线)“ZH”？K：“HZ”？N：“XD”？ C：“YD”? D： “XH”? A：“YH”?B:“R”?R:“T”?T: “LS”? U:Prog“7” :Lb1 1:Prog “1” :Prog “H”:Prog“2” :Goto1已知：ZH:直缓点桩号（对应变量K）; HZ:缓直点桩号（对应变量N ） ; XD:交点X坐标（对应变量C）; YD: 交点Y坐标（对应变量D）; XH:后方交点X坐标（对应变量A）; YH: 后方交点Y坐标（对应变量B）; R: 圆曲线半径（对应变量R）;T: 切线长（对应变量T）; Ls: 缓和曲线长 (对应变量U）; Q: Y=1右转 Z=-1左转（对应变量Q）输入：LX: 计算点桩号（对应变量V）; S: 计算点距该中桩距离（对应变量S）; W: 计算点与中桩连线与路线前进方向的夹角（对应变量W） 计算结果：计算点坐标X ,Y说明: 1.适用于标准的缓和曲线,既第一,第二缓和曲线长度相等。2.计算范围为：含两端直线段在内的整条曲线。 FHQX(复合曲线)Prog“8”: “R1”? H :“R2”? R: Prog “7”:Lb1 1: Prog “1” :Prog “F”: Prog “2”:Goto1已知：LA: 第一点桩号（对应变量K）; X1: 第一点X坐标（对应变量A）; Y1: 第一点Y坐标（对应变量B）;LN: 第二点桩号（对应变量N）; X2: 第二点X坐标（对应变量C）; Y2: 第二点Y坐标（对应变量D）； Q: 曲线转向(Y=1右转 Z=-1左转) （对应变量Q）; R1: 第一点曲线半径（桩号LA点) （对应变量H）; R2: 第二点曲线半径(桩号LN点) （对应变量R）（如：缓圆点或圆缓点的曲线半径等于圆曲线半径R； 直缓点或缓直点的曲线半径为+,可输入1G(千兆)）；其他项同HHQX说明：1、LA和LN在缓和曲线段上即可，不必是起点或终点；对应的曲率半径计算要正确。2、缓和曲线曲率半径计算公式：R=A2/L （A为缓和曲线参数；L为计算点至直缓点之间的缓和曲线长度）；3、LA可以大于LNSQX(竖曲线)“K1”?A:“K2”?B:“LD”?C:“HD”?D“R”?R:Lb11:“LX”?V:Pros“S”:“H=”HHGoto1已知：K1: 前一段路线纵坡（对应变量A） K2: 后一段路线纵坡（对应变量B） LD: 竖曲线交点桩号（对应变量C） HD： 交点高程（对应变量D） R: 竖曲线半径（对应变量R） 输入：LX 计算点桩号计算:LX桩号的路线设计高程 HLX S W(桩号反算)Lb10:“LC”？V:“X”？X：“Y”？Y:XZ1:YZ2:DO:Prog“P”: Z1E:Z2G: 0= (X-A) 2+(Y-B)2= E+.0001E: Pol(E-X,G-Y:Icos(J-M)K: V+KV:Lpwhile .005Abs(K:“LX=”:V“S=”:ISJ-MW: 180Abs(W =W-360Int (W360+.5WAbs(WW:“W=”:W DMSZ1 X: Z2 Y:Goto 0 已知：猜测桩号LC, 已知点P坐标（X,Y）。计算：已知点P（X,Y）相对于路线的位置。即对应的桩号LX,距中桩距离S，位于路线的哪一侧（W=90在右侧；W=-90在左侧）。说明：1、本程序在路线曲线要素已编入子程序“P”，即项目程序正常运行时方可运行。2、输入程序时注意区别：字母O和数字0、G和1G(千兆)、K和1k(千)、M和1M(兆)。HDZB(涵洞坐标)Lb11:“LX”？V:“W”？W：“LZ”？ E： “LY”？G: “B/2”？ B：Sin(W)=0=Goto1:VZ1:WZ2: EZ3: GZ4: BZ5:Prog“P”: XA: YB: Rec(Z4,M+W:X+IC: Y+JD: 0K: XA: Z4N:-Z3V: -Z2W: Abs(Z5 Sin(WS: Prog“Z”: Rec(S,M+W: IE: JG:“X1=”:X+E“Y1=”:Y+G“X4=”:X-E“Y4=”:Y-G0V: Prog“Z”: “X2=”:X+E“Y2=”:Y+G“X5=”:X-E“Y5=”:Y-GZ4V: Prog“Z”: “X3=”:X+E“Y3=”:Y+G“X6=”:X-E“Y6=”:Y-GZ1V: Z2W: Z3E: Z4G: Z5B: Goto 1 输入：LX ： 涵中心桩号（对应变量V） W ： 涵交角（对应变量W）LZ ： 左涵长（对应变量E） LY ： 右涵长（对应变量G） B/2 ： 涵正宽的一半（对应变量B）计算:涵放样点坐标（X1,Y1）-(X6,Y6)说明：本程序在路线曲线要素已编入子程序“P”，即项目程序正常运行时方可运行。HCJS (涵长计算)Lb11:“LX”？U:“W”？W： “B/2”？ B：“Z”？Z：“H+D”？P：BZ1:BS:Sin(W)=0=Goto1:Abs(SSin(WE:U+Ecos(WV:Prog“G”:Prog“I”: H+I(S-1-.75)-.03H: If8H-Z-P:Then(S+1.5(H-Z-P)Sin(WS:Else(S+1.5(H-Z-P)+.25(H-Z-P-8)+2) Sin(WS:IfEnd: “L=”:SZ1B:Goto 1 输入：LX ： 涵中心桩号（对应变量V） W ： 涵交角（对应变量W）B/2 ： 计算一侧的路面宽度（对应变量B） Z：涵底标高（对应变量Z） H+D ： 涵底至盖板顶高度（对应变量P）计算:对应W值一侧的涵长L（对应变量S）(计算另一侧涵长，将W输入对应的值，一般为W+180)说明：本程序在竖曲线要素已编入子程序“G”，横坡度要素编入子程序“I”时方可运行,内容也要根据边坡设计图纸做适当修改。TTBX (填土边线)“m1”?E: “m2”?F: “m3”?G: “B1”?O: “B2”?P: EZ3: FZ4: GZ5: OZ6: PZ7:Lb11:“LX”？V:“W”？W： “B/2”？ B：“Z”？Z：BZ1:ZZ2: BS: 1Abs(Sin(W =Goto1:Prog“G”: Prog“I”: H+I(S-1-.75)-.03H:HZ=Goto1:If HZ+8:ThenS+ Z3 (H-Z)S: Goto2: IfEnd: If HZ+20:ThenS+ 8Z3 + Z4 (H-Z -8)+OS: Goto2: IfEnd: If HZ+32:Then S+ 8Z3 + 12Z4 +Z5 (H-Z -20)+O+ Z7S: Goto2: IfEnd: HZ+32=Goto1: Lb12: “TH=”:H-Z“S=”:SProg“P”: Prog“2”: Z1B: Z2Z: Z3E: Z4F: Z5G: Z6O: Z7P: Goto 1 输入：一级、二级、三级坡比分别为m1、m2、m3；一级、二级平台宽度分别为B1、B2 LX ： 路中心桩号（对应变量V） W ：在路的哪一侧，交角90（对应变量W） B/2 ： 半幅路面宽度（对应变量B） Z：实测地面标高（对应变量Z）计算:计算显示填土高度“TH=”;填土边线到中桩距离“S=”;对应W值一侧填土边线坐标(X,Y)说明：本程序在在路线曲线要素已编入子程序“P”，竖曲线要素已编入子程序“G”，横坡度要素编入子程序“I”时方可运行,内容也要根据边坡设计图纸做适当修改。(本次边坡设计见下图，距中1m为高程设计线，土路肩宽度0.75m.暂时编辑三级坡以下的。)HJ(后方交汇)Prog“3”:“H”?H:“IH”?K:“RH”?R:“HR1”?O:“HD1”?L:“HR2”?P:“HD2”?S:“VD”?T:Prog“6”:IQ:P-OW:Rec(S,W:Pol(I-L,J:QIQ:Rec(QS,W:0U:0V:QLE:0G:IM:JN:Prog“B”:“Z=”:H-T-K+RZ已知：两已知点坐标（X1,Y1）;（X2,Y2）, H： 任意一点高程 ；实测：IH： 仪器高；RH： 棱镜高； HR1： 第一点水平角读数；HD1： 第一点水平距离读数；HR2： 第二点水平角读数； HD2： 第二点水平距离读数； VD： 垂直高差读数(对应H) 计算并显示 ： 所求测站点坐标(X,Y,Z) FWFS(方位反算) “X0”? A: “Y0” ? B : Lb1 1: “X”?X:“Y”?Y:Pol(X-A,Y-B: “L=”:IL“HR=”:J+180(1-AbsJJDMSGoto 1已知：两点坐标P0(X0,Y0)；P(X,Y)计算:P0点到P点的距离L 和 方向角HR（显示度、分、秒）。MJ(面积)Prog“3”:Prog“6”: 0S: IL: 0E:Lb11: “X”？X：“Y”？Y:.5(AD+BX+CY-DX-AY-CBK: “S=”: S+KS(X-C)2+(Y-D) 2G:(X-A)2+(Y-B) 2H: “L=”: L-E+G+HLXC: YD: HE:Goto 1已知:多边形各角点坐标，按顺时针或按逆时针顺序输入各角点坐标, 计算闭合多边形面积S（为累计值）及周长L。示意图如下:JD(交点)Prog“3”:Prog“6”:Lb11:“X3”？E：“Y3”？G:Pol(E-A,G-B: “S=”：Isin(J-FS“X4”？M：“Y4”？N：Pol(M-E,N-G: J-FW:IfInt(W180)W180:Then“XD=”: (CEB-CMB+AMD-AED+CGM-AGM+AEN-CEN)(EB-MB+MD-ED+CG-AG+AN-CNX“YD=”: (DGA-DNA+BNC-BGC+DNE-BNE+BMG-DMG) (GA-NA+NC-GC+DE-BE+BM-DMY IfEnd:Goto 1已知：第一条直线上两点坐标分别为P1（X1,Y1）、P2（X2,Y2）; 第二条直线上两点坐标分别为P3（X3,Y3）、P4（X4,Y4）。计算并显示： P3（X3,Y3）到直线（P1P2）的垂直距离S; 以及两直线交点坐标(XD,YD)。示意图如下：子程序 1 “Lx”？V：“S”？S：“W”？W：Return2IfS0:ThenRec(S,M+W:Else Rec(-S,M+W-180: IfEnd:“X=”: X+IX“Y=”: Y+JY：Return3 “X1”? A: “Y1”? B: “X2”? C: “Y2”? D：Return40Z:KN= 180Z: cos(ZP：Return5A+IX: B+JY：Return6Pol(C-A,D-B: JF：Return7 “Y=1, Z=-1”?Q：Return8 “LA”? K: “X1”? A: “Y1”? B: “LN”? N: “X2”? C:“Y2”?D:Return9NK:CA:BD: Return XRULO:L-L(3)(40O2)+L(5)(3456O(4)-L(7)(599040O(6X:L2(6O)-L(4)(336O(3)+L(6)(42240O(5)-L(8)(9676800O(7Y：ReturnHU=0=1U:Prog“4”:Prog“6”:P(V-KL:0H:If(PV).5P(K+N:ThenP(N-VL:180H:F+180Q(N-K-PU)(R)+HF:IfEnd:Rec(T,F: C-IE: D-JG: IfL0:ThenRec(-L,F+180: F+Z+HM: 1MO:Goto2:IfEnd:IfLU:ThenProg“X”:Pol(X,Y:Rec(I,F+PQJcos(H: F+90PQLcos(H)(O)+Z+HM:Goto2:IfEnd: 45(2L-U)(RJ: 45U(RI: F+2PQJcos(H)+Z+HM: UL:Prog“ X”:Rec(2Rsin(J-I,J+I: X+IX: Y+JY:Pol(X,Y:Rec(I,F+PQJcos(H:Lb1 2: E+IX: G+JY：ReturnS1E:AB=-1E:.5RAbs(A-BT: AM:VC= BM:If T Abs(V-C:ThenD+M(V-CH:ElseV-C+TG: D-AT+Gtan(tan-1 (A)+90EG(RH：IfEnd:ReturnFR(N-K)(H-RM: N-K+MU: ML:Prog“X”: XE: YG: UL: Prog“X”:Prog“6”: IZ:Pol(X-E,Y-G: ZIZ: F-QJF:V=K=K+1V:V-K+ML:Prog“X”:Pol(X-E,Y-G:Rec(ZI,F+QJ: F+90QL(OM：Prog“5”：ReturnZProg“4”:Prog“6”:F+ZM: IfVK:ThenRec(V-K,M:Else Rec(K-V,M+180: IfEnd: Prog“5”: 1GO: 1Q：ReturnYProg“6”:Prog“4”:F-90Q(N-K)(R)+ZE:90Q(V-K)(RG:IfQRsin(G) 0:ThenRec(2QRsin(G, E+G:Else Rec(-2QRsin(G, E+G -180: IfEnd: E+2GM:Prog“5”: RO：ReturnBProg“6”:IQ:Pol(M-E,N-G:QIQ:F-JF:Pol(U-E,V-G:Rec(QI,F+J:”X=”: A+IX“Y=”: B+JY：Return匝道 “1-4：WA-WB：5-8：EA-ED：9WN：10WS：11EN：12ES”?T：Lb1 0：Porg“1”:If T=1:Then Porg“P”：Porg“2”：Goto 0: IfEnd: If T=2: Then Porg“WA”：Porg“2”：Goto0: IfEnd: If T=3:Porg“WB”：Porg“2”：Goto0: IfEndP（示例：永蓝主线K99+000-K106+012）第二缓和曲线终点桩号N:IfVN:Then第二缓和曲线起点桩号K:起点x坐标A:起