FX5800计算器公路测量常用程序集.doc

fx5800计算器公路测量常用程序集2.3版一、程序功能 主要功能：采用交点法方式计算多条线路坐标正反算，可算任意复杂线型及立交匝道，包括c型，s型、卵型、回头曲线等；极坐标放样，全线路基边坡开挖口及坡脚放样计算、路基任意点设计高程、横坡、设计半幅宽度、隧道欠超挖放样计算。 新版本优化：1、优化程序语句、2、取消原线元法计算，化线元法为交点法。3、高程计算修改,4、附计算坐标及高程参数设计要素输入实例二、源程序（绿色为程序名；蓝色为输入计算器内容）红色为网友调试后修改笔误。0.总主程序（1、坐标计算放样程序2、坐标反算程序；3、高程计算查阅程序；4、路基半幅标准宽度查阅程序；5、路基边坡及开挖口放样程序；6、路基标准距离放样；7、桥梁锥坡计算放样程序；8、极坐标计算程序；9、隧道超欠挖计算程序）运行后按19数子约半秒，则选择1至9的程序，返回时，在桩号输入-1,返回选择计算类型。输入-2,返回选择线路。程序名：0（数子0）1.主程序：一般坐标计算及放样程序程序名:1xyprog dx:progqx:90a: pj1?a:ab: pj2?b:az1:bz9: do:km?z:z=-1=stop:z=-2=break: progk: ?d: 0l: l0?l:lz10: porga:z2+z1 a :a-z9e:i+dcos(a)+z10cos(ex:j+dsin(a)+z10sin(ey:progxy:progjs:lpwhile z-12.主程序：由大概桩号及坐标反算桩号及距离程序名:2zdprogqx:do: km?z:z=-1=stop:z=-2=break: progk: xo?x: y0?y: porgb: progzd:lpwhile z-13.主程序：任意点高程计算及横坡程序名:3gcprogqx:0b: h-b?b: bz9: do: km?z:z=-1=stop:z=-2=break:?d:progh:fix 3: h=:lcoate 6,4,h-z9： i=:locate 6,4,i：lpwhile z-14.主程序:任意桩号位置的标准路基宽度计算程序名:4gdprogqx: do: km?z:z=-1=stop:z=-2=break:-1d:progc:lp：1d: progc:fix 3: lgd=:lcoate 6,4,p： rgd=:locate 6,4,l：lpwhile z-15、主程序:路基开挖边线及填方坡脚线放样程序程序名：4bpprog dx:progqx:0.5b: th-gd?b: bz7: do:km?z:z=-1=stop:z=-2=break: progk: “x0?x:y0?y:z8m：m0?m:mz8: porgb:dp :progc:if d0:then goto t: ifend:-gh:if he+f:then l+m+n+u+c(h-e-f)+bf+aes:goto 0:ifend:if he:then l+m+n+ b(h-e)+aes:goto 0:ifend:if hw then l+z7+v+j(g-w)+iws:else l+z7+gis:ifend:lbi 0: fix 2: pd ：abs（d）-st: lx=:locate 6,4,t：progzd: tw=:lcoate 6,4,glpwhile z-16、主程序：定值放样路基桩计算填挖值程序名：6fmprog dx:progqx: do:km?z:z=-1=stop:z=-2=break: progk: ?d: porga:z2+90a :i+dcos(ax:j+dsin(ay:progxy:progjs:z8m:”m0”?m:mz8:prog “h”:” tw=”:locate 6,4,h-z8: lpwhile z-17、主程序：桥梁锥坡计算放样程序名：7zpprogdx:progqx:z0?z:ld:z-,y+? d:abs(d) r: lr?r: dz7:rz10:zz9:prog”c”:lz10:prog”k”:do:z10y:l0:sz+,dz-?y: y=-1=stop:y=-2=break: (z72（1-y2/z102）)x：z9+yz: prog”c”: if d0:then l+xd:else l-xd:ifend:prog”a”: z2+90a :i+dcos(ax:j+dsin(ay:progxy:progjs: lpwhile y-18、主程序：由坐标计算方位角及距离程序名：8jsprog dx:do:?x: x=-1=stop:x=-2=break:?y:prog”js”: lpwhile x-19、主程序：隧道超欠挖放样计算程序名：9sdprogdx:0.65a:”cqhd”?a:az7:do:km?z:z=-1=stop:z=-2=break: progk: “x0?x:y0?y:z8m：m0?m:mz8: porgb:dz9:0d:prog”h”:prog”sd1”:z9 d:z8m:z7x:rcos(pc:h+aa:h+bb:c+ae:if mgoto 1: (s2-(b-a)2)sabs(st:abs(d-q)-tt:(t2+(m-b)2)-v-xl:goto 0:ifend:lbi 1:(d-q)2+(m-a)2)-r-xl:lbi 0:fix 2: ” lx=”: locate 6,4,l:prog”zd”:“ h0=”: locate 6,4,m-hlpwhile z-110. 交点法正算子程序()程序名：apabs(pw:o-t+rabs(p)180+h2-n2q:hc:0m:if zo-t:then z-os:gz2:goto 1: ifend: if zo-t+h:then z-o+ts:prog “hx”:g+wkz2:a-ta:ge:goto 2:ifend:if zq: then 180(z-o+t-0.5h)rs : prog”ab”:a+r(1-cos(se：b+rsin(sa:eb:rm:g+wsz2:a-ta:ge:goto 2: ifend:if zq+n:then q+n-zs:nh:prog “hx”:g+pe:e-wkz2:f-aa:ch:goto 2:ifend:z-q-n+fs:g+pz2: lbi 1:u+scos(z2i:v+ssin(z2j:return:lbi 2:u+acos(e)-wbsin(ei:v+asin(e)+wbcos(ej11. 交点法反算子程序()程序名：b0d:do:z+dd:prog”a”:pol(x-i,y-j+10-9:j-z2 j:isin(js:icos(jd:if m0:then pol(m-ws,d:jm180d:ifend: lpwhile abs(d.001:z+dz:sd12.缓和曲线段直角坐标计算子程序(b)程序名：abh2r24-h(4)2688r(3)+h (6)506880r(5)- h (8)154828800r(7)a:h2-h(3)240r2+ h (5)34560r(4)- h (7)8386560r(6)+ h (9)3158507520r(8)b13. 计算缓和曲线段坐标增量及偏角子程序(hx)程序名：hxs-s(5)40r2h2+s(9)3456r(4)h(4)-s(13)599040r(6)h(6)+s(17)175472640r(8)h(8)a:s(3)6rh-s(7)336r(3)h(3)+s(11)42240r(5)h(5)-s(15)9676800r(7)h(7)+s(19)3530097000r(9)h(9)b:90s2rhk:if s=0:then 0m:else rhsm:ifend14高程计算子程序（h）程序名：hprog “s1”:n-ia:aabs(aw:rabs(a2t:if zc-t:then g+i(z-ch:goto 0:ifend:if zc+t:then g-ti+(i+(z-c+t)w2r)(z-c+th:goto 0: ifend:ni: g+i(z-ch:ilbi 0:if d=0:then hh:return:else prog “i”:h+i(abs(d)-kh:ifend15.高程超高计算横坡子程序（i）程序名：iprog “i1”：if s=0:then mi: return:else if w=1:then (n-m)(z-c)s+mi:return:else (3(z-c)s)2-2(z-c)s)(3)(n-m)+mi:ifend:16路基宽度子程序程序名： c1s: prog “g1”:z-ce:(b-a)es+al17：极坐放样计算程序(计算放样点至置仪点方位角及距离)程序名：jsz5k：z6l:pol(x-k, y-l:if j0:then j+360j:ifend:int(j)+0.01int(60frac(j)+0.006frac(60frac(jj: fix 4:” fwj=”: locate 6,4,j：fix 3:” s=”:locate 6,4,i18导线点子程序（dx）程序名：dxz5k:z6l:“xz”?k:”yz”?l:kz5:lz619线路选择子程序(线路选择输0时。则输曲线参数（qx）程序名：qx1a:”qx:x1=1,x2=2”?a:az4:if a=0:then ”jdkm”?o:”jdx”?u:”jdy”?v:”fwj”?g:”zj”?p:”ls1”?h:”ls2”?n:”t1”?t:”t2”?f:?r:ifend20线路中线元段判别子程序（k）程序名：kif z3=0:then return: ifend:if z3=1:then prog “x1”:return:ifend：if z3=2:then prog “x2”:return:ifend21坐标参数矩阵调用程序（f）程序名：fmat f1,1o: mat f1,2u: mat f1,3v: mat f1,4g: mat f1,5p: mat f1,6h: mat f1,7n: mat f1,8t:mat f1,9f: mat f1,10r22显示子程序程序名：xyix 3: =：locate 6,4,： ：locate 6,4,：程序名：zdfix 3:：locate 6,4,：locate 6,4,23数据子程序(附后示例)匝道线元转交点数据程序（t）程序名：tnorm 2:“qkm”？z：“xq”？u：“yq”？v：“ls1”？m：“ls2”？o：“a1”？a：“a2”？b：“fwj“？g：”zj“？p：？r：pabs(pq:a2rh:hmh:h-ms:b2rn:n-oc:o n=on:h=0=10-9h:hx:90s2(rhk:g-qkg:n=0=10-9n: 90c2(rnd:p+qk+qdp:prog”ab”: (h2-h2)24rsin(abs(p)-(h(4)-n(4)2688r(3)sin(abs(p)+(h (6)-n(6)506880r(5)sin(abs(p)-(h (8)-n(8)154828800r(7)sin(abs(pe:(r+a)tan(abs(p)2)+b-et:nh:prog”ab”: (r+a)tan(abs(p)2)+b+ef: xh: z-s+to:prog”hx”:u-(a-t)cos(g)+qbsin(gu:v-(a-t)sin(g)-qbcos(gv:fix 3:”jdkm=”:locate 6,4,o:”jdx=”:locate 6,4,u: ”jdy=”:locate 6,4,v: ”fwj=”:gdms”zj=”pdms”ls1=”: locate 6,4,h: ”ls2=”: locate 6,4,n: ”t1=”: locate 6,4,t: ”t2=”: locate 6,4,f”r=”:r说明：运行后“qkm”？输入起点桩号；“xq”？ “yq”？输入起点x、y坐标；“ls1”？ “ls2”？输入前段缓和曲线长度，后段缓和曲线长度；“a1”？ “a2”？输入前段缓和曲线参数，后段缓和曲线参数；“fwj“？输入起点切线方位角；”zj“？输入转角（终点方位角减起点方位角之差，左转为负，右转为正）；计算结果要素同qx子程序中要素字母。-1程序名：x1（线路1坐标计算要素程序）if z241.086:then -336.478，2600441.639，514392.519， 280。2028.78”， 203。0828.32”，106.667，81.667，-270.668，-283.056，60mat f:prog”f”: ifend-2程序名：x2（线路2坐标计算要素程序）if z273.006:then 171.674，2600784.547，514635.863，58。4917.2”， 68。0926.91”，93.889，125，674.184，771，180mat f:prog”f”: return: ifendif z466.227:then 341.96，2600664.76，514799.865，123。4325.9”， 15。5849.5”,0，138.211，68.954，125.273，445mat f:prog”f”: return: ifend说明：1、交点法：if z计算交点段终点桩号：then交点桩号，交点x坐标，交点y坐标，起点至交点方位角，交点转角，交点前段缓和曲线长度，交点后段缓和曲线长度，交点前段切线长度，交点后段切线长度，圆半径 mat f: prog”f”: return: ifend（转角则为交点段终点方位角减起点方位角，分正负，左偏为负，右偏为正）如整条线无曲线，全为直线段，则交点为终点桩号，xy为终点坐标，方位角为起点方位角，转角和圆半径为1，其它全为0。-1程序名：s1(线路1高程竖曲线要素子程序)if z3=1:then goto 1: else prog “s2”: return: ifend: lbi 1:if z241.09:then 95.356c:40.742g:2000r:-.03155i:.00812n:ifend:-2程序名：s2(线路2高程竖曲线要素子程序)if z133.601:then 106c:45.12g:1525r:-.0064i:.0298n:return:ifend:if z466.23:then 298c:50.841g:1945r:.0298i:-.01968n: return:ifend:程序字母说明：c-竖曲线交点桩号；g-交点桩号高程（未竖曲线调整的）；r-竖曲线半径；i竖曲线前纵坡; n竖曲线后纵坡-1程序名：i1（线路1高程超高参数子程序）0k:1w: if z4=1:then goto 1:else prog “i2”:return: ifend:lbi 1:if d0:then goto l:else goto r:ifend:lbi l:if z159.419:then 0s:.06m: return:ifend:if z187.987:then 28.568s:159.419c:.06m:.02n: return:ifend:if z241.086:then 0s:.02m: return:ifend:lbi r:if z159.419:then 0s:-.06m: return:ifend:if z187.987:then 28.568s:159.419c:-.06m:-.02n: return:ifend:if z241.086:then 0s:.-02m: return:ifend-2程序名：i2（线路2高程超高参数子程序）if d0:then goto l:else goto r:ifend:lbi l:if z77.872:then 0s:.02m:return:ifend:if z93.889:then 16.017s:77.872c:.02m:.05n:return:ifend:if z198.567:then 0s: .05m:return:ifend: if z273.006:then 74.439s:198.567c:.05m:.02n:return:ifend:if z335:then 0s: .02m:return:ifend: z346.235:then 11.235s:335c:.02m:.015n:return:ifend:lbi r:if z77.872:then 0s:-.02m:return:ifend:if z93.889:then 16.017s:77.872c:-.02m:-.05n: return:ifend:if z198.567:then 0s: -.05m: return:ifend: if z273.006:then 74.439s:198.567c:-.05m:-.02n: return:ifend:if z335:then 0s: -.02m: return:ifend: z346.235:then 11.235s:335c:-.02m:-.015n: return:ifend:子程序中字母表示说明：k中央分隔带半幅宽（中桩标高至中桩中线距离，无分隔带则为0）；w超高方式参数（w=1为一般直线方式超高，w=2为三次抛物线方式超高）;s超高渐变段距离（不是渐变段则输入0）；c-超高渐变段起点桩号（不是渐变段，无需输入）；m超高段起点横坡，n-超高渐变段终点横坡（不是渐变段，无需输入）。-1程序名：g1(线路1路基标准半幅宽度参数子程序)if z4=1:then goto 1:else prog “g2”:return: ifendlbi 1:if d0:then goto l:else goto r:ifend:lbi l:if z241.086:then zc:3.5a:return:ifend:lbi r: z40.857:then 0c:5a:5.25b:40.857s:return:ifend:if z159.419:then zc:5.25a:return:ifend: z187.987:then 159.419c:5.25a:5.5b:28.568s:return:ifend:if z241.086:then zc:5.5a:return:ifend-2程序名：g2(线路2路基标准半幅宽度参数子程序)if d0:then goto l:else goto r:ifend:lbi l:if z466.227:then zc:3.5a:return:ifend:lbi r: z321.611:then zc:5a:return:ifend:if z346.235:then 321.611c:5a:5.75b:24.624s:return:ifend:if z466.227:then zc:5.75a:return:ifend子程序中字母表示说明：c宽度渐变段起点桩号（不是渐变段c=z）；a宽度渐变段起点宽度; b宽度渐变段终点宽度（不是渐变段，无需输入）；s-宽度渐变段距离（不是渐变段，无需输入）。-1程序名：sd1(隧道参数子程序)if z4=1:then 5.1r:1.8a:1.4v:1.8b:0q: 109。0130”p: return:ifend:if z4=2:then 8.75r:-1.1a:4.77v:1.47b:-.25q: 51。2857.1”p: return:ifend:子程序中字母表示说明：r隧道上拱半径；a隧道上拱圆心至设计面高度；v隧道中腰部半径；b中腰圆心至设计面高度；q设计隧道中线与设计路基中线偏移值，左为负数，右为正数；p上部拱跨隧道半幅角度；-1程序名：w1(线路1路基填挖边坡参数子程序)if z4=1:then goto 1: else prog “w2”: return:ifend:lbi 1: 10e:10f:2.6m:2n:2s:1.5i:1.75j:8w:2v:if d0:then goto l:else goto r:ifend:lbi l:1a:1.25b:1.25c:return:lbi r:if z110:then 1.25a:1.25b:1.5c: return:else goto l: ifend-2程序名：w2(线路2路基填挖边坡参数子程序)10e:10f:2.6m:2n:2s:1.5i:1.75j:8w:2v:if d0:then goto l:else goto r:ifend:lbi l:0.75a:0.75b:0.75c:return:lbi r:if z445:then goto l:else if z506.9:then 0.25(z-445)21.2q:.75+qa:ab:ac:return:ifend子程序中字母表示说明：a挖方第一阶边坡边率；b挖方第二阶边坡边率；c挖方第三阶边坡边率；e挖方第一阶高度；f挖方第二阶高度；m挖方路基碎落台及水沟宽度；n挖方第二台阶平台宽度；s挖方第三台阶平台宽度；i填方第一阶边坡边率；j填方第二阶边坡边率; w填方第一阶高度；v-填方第二台阶平台宽度。（注：本程序只做出挖方三个台阶，填方二个台阶，如需增加，先需在程序5bp程序相应增加。再如一标段有坡率及坡高不一样时，可以照其它参数程序一样，用判别语句。）三、使用说明1、规定 (1) 坐标计算以交点方式计算，一般情况下一个交点计算一段线是：直线+缓和曲线+圆曲线+缓和曲线+直线。在匝道中先根据实际线来调整，可以单独计算一段不完整的缓和曲线+圆曲线。也可以单独计算一段缓和曲线或一段圆曲线，需先运行t程序整理成交点数据。无缓和曲线时，缓和长度输0。(2) 程序中：10-9 为10的负9次方即0.000000001; 为圆周率即3.1415926; 为根号； =：x前有空格二个，为了显示效果 (3) 为了区别字母o及数子0，特把字母0加粗。2、输入与显示说明 输入提示部分：a:xy=1， zd=2 ，gc=3， gd=4， bp=5， fm=6， zp=7， js=8，sd=9? 按对应数子键约半秒:1、坐标计算放样程序2、坐标反算程序；3、高程计算查阅程序；4、路基半幅标准宽度查阅程序；5、路基边坡及开挖口放样程序；6、路基标准距离放样；7、桥梁锥坡计算放样程序；8、极坐标计算程序；9、隧道超欠挖计算程序xz？yz？ 为置仪点坐标，第一次运行输一次，以下运行不需再输入，下次再运行此程序调上次输值，如需改，则输入，不需改，则按确认键。(后附有导线点坐标存储子程序,有兴趣的可加上)“qx:x1=1,x2=2”?选择经路，坐标计算中输入0值，变为输入线元段曲线要素。输1为第一条线路。输2为第二条线路，本程序暂设计实例为2条线路。？正算时所求点的里程（反算输入大概桩号）输入-1，返回输入选择计算类型，输入-2,返回选择线路。 ？正算时所求点距中线的边距(左侧取负值,右侧取正值，在中线上取零) x0 ？反算时所求点的x坐标（放样程序中实测x坐标） y0 ？反算时所求点的y坐标（放样程序中实测y坐标）显示部分：l0? 涵洞放样程序中涵距中心桩号前后距离,前为正,后为负.(即涵洞半幅宽度)m0? 放样程序中原地面标高输入h-b？计算面至路面高差（默认路面为0）pj1? pj2? pj1为路线中线与边线的右夹角。第二夹角为边线至路前进方向夹角,默认值均为90度.th-gh？边坡放样中填方超放宽值，默认值为0.5m.ld:z-,y+? 桥梁锥坡向左右方向边长，左为负数，右为正数。(即为锥坡高程减地面高程后乘以填方坡比)lr? 锥坡向路线前进或向后退方向边长。z0? 锥坡顶边桩号l0:sz+,dz-?锥坡放样中，椭圆上分段放样距离，小桩号锥坡输正数，大桩号锥坡输负数。“cqhd”? 隧道超欠挖程序中衬砌厚度, 默认值为0.5m.qx子程序中和t程序结果中：“jdkm”？输入交点桩号；“jdx”？ “jdy”？输入交点x、y坐标；“fwj“？输入起点至交点方位角”zj“？输入转角（终点方位角-起点方位角之差，左转为负，右转为正）;“ls1”？“ls2”？输入交点前段缓和曲线长度，交点后段缓和曲线长度；“t1”？ “t2”？ 输入交点前段切线长度，交点后段切线长度；”r”?圆半径输出部分 x=* 正算时，计算得出的所求点的x坐标 y=* 正算时，计算得出的所求点的y坐标fwj=* 正算时，计算得出的所求点的至置仪点方位角s=* 正算时，计算得出的所求点的至置仪点距离=* 反算时，计算得出的所求点的里程 =* 反算时，计算得出的所求点的边距h=* 所求点位置设计计算面顶标高i=* 所求点位置设计路面横坡（d输0时为，本段路基纵坡）lgd=* 所求点位置设计路面顶左半幅标准宽度lgd=* 所求点位置设计路面顶右半幅标准宽度lx=* 边桩放