《fx-5800公路综合放样及常用程序》由会员分享，可在线阅读，更多相关《fx-5800公路综合放样及常用程序（6页珍藏版）》请在金锄头文库上搜索。
1、曲线任意里程中边桩坐标正反算、放样及公路常用（fx-5800）程序 源程序 1.对称型曲线(对称 缓+园+缓 含两侧直线,不能反算)1.1 FYZ (放样主程序) “9=K“?J: If J=9:Then Goto F:IfEnd:“(+,-)“?A:A/ Abs (AZ: Abs (AA:?R:“Lsci“?S:R+S/24/R-S(4) / 2384/R（3）Q: “T=“:S/2- S(3)/240/R+Anstan .5AT “L=“:RA/1r+SL“E=“:(R+S /6/R-S(4) /336/R(3)-R(1-cos (90S/R)/cos (.5A)-RE “JD“?W:“Z
2、H=“:W-TH(“HY=“:H+S “QZ=“:H+.5L “YH=“:H+L-S)括号内可不要 “HZ=“:H+LL “XJ“?D:“YJ“?K:“FWJ“？B:“XCe“?O:“YCe“?U:“Xh“?G:“Yh“?M: Pol(G-O,M-U: If JA=转角（左-，右+）: S=缓和曲线长: R=半径: W=交点桩号:D=交点X坐标:K=交点Y坐标:B=计算方位角 : R=另一交点编号；有几个=就有几个）1.4 子程序 FY （1、2、3程序公用） Pol(X-O,Y-U: If JK? 输入9进入桩号 (+,-)? 输入转角左负右正“ Lsci“? 输入缓和曲线长 “JD“? 输
3、入交点桩号“XJ“YJ“? 输入交点坐标 “ FWJ“? 输入计算方位角 “Xce“Yce “? 输入测站坐标 “Xh“Yh “? 输入后视坐标K？ 输入桩号 L(-Z,+Y) 输入边距左负右正 (0,1)？ 输入0夹角放样（即后视归0）；输入1方位角放样1.8显示部分： T= 切线长； L= 曲线长； E= 外距；ZH= 直缓桩； HY= 缓圆桩； QZ= 曲中桩；YH= 圆缓桩； HZ= 缓直桩； X= X坐标；Y= Y坐标= 1、测站到后视方位角；2、放样夹角或方位角P= 1、测站到后视距离； 2、放样距离2正反算程序2.1正算主程序 GSZS “X0“？Z:“Y0“？S:“K0“？W:
4、“F0“？G:“KN“？H:“R0“？N:“RN“？R:“Q(-Z +Y) “？Q:“Xce“？O:“Yce“？U:“Xh“？A:“Yh“？B:Pol(A-O,B-U: If J就有几个) KH=Z=*:S=*:W=*:G=*:H=*:N=*:R=*:Q=*: (注：如有多个曲线元要素继续添加入数据库DAT-M中) 2.5、M线（坐标正算）组合程序 Z-S Prog”DAT-M”:Prog”GSZS”2.6、M线（坐标反算）组合程序 F-S Prog”DAT-M”:Prog”GSFS” 说明： 一、俩个程序优缺点：1.优点：1程序一个交点内只须输一组数据（直+缓+园+缓+直）、并且均是 曲线要
5、素表中的已知数据；2程序可以正反算、任意曲线均可用。角度均为度分 秒2.缺点：1程序不能反算、要求对称；2程序输入数据库较大、曲线起点坐 标及切线方位角在曲线要素表中无须计算。3.建议都输入计算器、能用1程序就尽量用1程序、特殊曲线才用2程序。子 程序1.4“ FY”公用子程序；号会与+号混淆所以在程序中除号均以/号写。二、程序功能及原理 1.功能说明： 本程序由两个主程序正算主程序(GSZS)、反算主程序(GSFS)和两个子程正算子程序(SUB1)、线元数据库(DAT-M)构成，可以根据曲线段直线、圆 曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点 切线方位角、终点里程
6、、起点曲率半径、止点曲率半径）及里程边距或坐标， 对该曲线段范围内任意里程中边桩坐标进行正反算。本程序在CASIO fx-5800P 计算器上运行。由于加入了数据库(DAT-M)，可实现坐标正反算的全线贯通。2计算原理： 利用Gauss-Legendre 5点通用公式计算线路中边桩坐标并计算放样数据。利用待求点至线元起点切线作垂线，逐次迭代趋近原理反算里程及边距。三、使用说明 1、规定(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往 左偏时， Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。(2) 当所求点位于中线时，L=0；当位于中线左侧时，L取负值；当位于 中线
7、右侧时，L取正值。(3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次 代替。(4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径 均等于圆弧的半径。 (5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10 的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接 时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆 曲线的半径。(6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设 计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接 时，曲率半径等于设计规定的值；与圆曲线相接
8、时，曲率半径等于圆曲线的半 径。 (7)曲线元要素数据库（DAT-M）可根据线型不同分为各个线元段输入到DAT-M中， 即分为直线段、缓和曲线、圆曲线等。 (8)正算时可仅输入里程和边距及右交角可实现全线计算，但反算时只能通过首 先输入里程K值读取数据库DAT-M，计算器自动将里程K所在线元数据赋给反算主 程序GSFS进行试算，试算出的里程和边距须带入正算主程序GSZS中计算坐标， 若坐标吻合则反算正确。 2、输入与显示说明（1）输入部分：X0 ？线元起点的X坐标 Y0 ？线元起点的Y坐标 K0 ？线元起点里程F0 ？线元起点切线方位角KN ？线元终点里程R0 ？线元起点曲率半径RN ？线元止
9、点曲率半径Q ？ 线 元左右偏标志(左偏Q=-1，右偏Q=1，直线段Q=0)K ？ 正算时所求点的里程L ？ 正算时所求点距中线的边距(左侧取负值，右侧取正值，在中线 上取零)X ？ 反算时所求点的X坐标Y ？ 反算时所求点的Y坐标A、 B、Z4 是Gauss-Legendre求积公式中的插值系数C 、E、0.5 是Gauss-Legendre求积公式中的求积节点 （2）显示部分： X= 正算时，计算得出的所求点的X坐标Y= 正算时，计算得出的所求点的Y坐标K= 反算时，计算得出的所求点的里程L= 反算时，计算得出的所求点的边距3.直线程序 ZX (已改为5800) “X0“?A: “Y0“?
10、B: “T0“?T: “K0“?E：Lb1 1：?K: ?L: “X“:A+ Rec(Pol(K- E,L,T+JX“Y“:B+JYProg“FY“:Goto 1 (说明：已知某一直线线路上一点的坐标为X0，Y0，该点里程为K0以及该直线沿 前进方向的坐标方位角T0，按照提示输入以上数据，如果输入某一里程K和边距 L（向右移为正值，向左移为负值）即可计算出该里程点的法线外移桩坐标 X、Y，当L=0时计算结果为中桩坐标。) 4.超高加宽程序 CG JK “ZH“? Q: “LSci“? L: “JA K“? C:“cao g“？I:“Q=0，Z=1“?F:Lbl 1:Abs(Q - A“K“)/
11、L K: If F=0:Then 4.25+(4K3-3K4)C B:Else KC+4.25B:IfEnd: “B=“:B?B: “E=“: (Abs (A - Q) /L(I%+2%) 2%)BE Goto 1 说明: 4.25为 路基宽,(路基是多宽改为多宽,0.02为路基横坡，若不是2%须改) 输入部分： ZH? 直缓桩号; LSci? 缓和曲线长; JA K? 全加宽值 cao g ? 超高横坡; K ? 桩号; B ? 路基宽 Q=0，Z=1 ? 0为抛物线加宽, 1为直线加宽; 显示部分：B 加宽后的路基宽； E= 超高值；5.竖曲线程序 SHU Q X(已改为5800) “K0
12、+x“? A : “H0“? G:?T:?R: “1i“?C: “2i“？D: Lbl 1: ? K: Abs (A-K) B: If BT:Then 0E :Else (T-B)/2/RE:IfEnd: If DC:Then -EE:IfEnd: “E=“:E If KA:Then G-BC%+EH:Else G+BD%+EH:IfEnd: “H=“:H?L:“H1=“: H- 2%LGoto 1 说明: 输入部分： K0+x? 变坡点桩号; H0? 变坡点高程; T? 切线长 R ? 半径; 1i ? 前纵坡; 2i ? 后纵坡 K ? 要算点桩号; L? 边距显示部分：E 竖曲线改正值；
13、 H= 要算点高程；H1=边桩高程（若横坡不是2%须改）后方交会测量程序后方交会测量程序 HFJH（5800）“X1“？C：“Y1“？D：“X2“?G：“Y2“?M：“P1“?A：“P2“?B： I=0：J=0:POL（G-C，M-D: COS-1（I2+A2-B2）/（2IA） ）: “X=“:C+Rec(A,J+AnsO”Y=”:D+JU注：输入部分：X1? 左边点 X 坐标，Y1? 左边点 Y 坐标X2？ 右边点 X 坐标，Y2？ 右边点 Y 坐标 P1？ 左边边长， P2？ 右边边长显示部分：X= 置仪点 X 坐标； Y= 置仪点 Y 坐标后方交会测量程序后方交会测量程序 HFJH（4800）C“X1“D“Y1“G“X2“M“Y2“A“P1“B“P2“：I=0：J=0:POL（G-C，M-D: COS-1（I2+A2-B2）/（2IA） ）:O“X“=C+Rec(A,J+AnsU“Y“=D+J
