任意桩号坐标及高程
一、程序功能
本程序由一个主程序(ZHUCHENXU)和几个子程序——正算子程序(SUB1)、反算子程序(SUB2)、数据库(SUB3)、曲线要数（SUB4,4-1）、边坡放样程序(BIANPO)，可以根据直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标及中桩高程，对该曲线段范围内任意里程中边桩坐标进行正反算。 二、源程序 增加变量18→DimZ 1.主程序(ZHUCHENXU)
\"1.SZ → XY，2.XY → SZ，3.SDYX，4. BIANPO \"？N：
Lbl 1：”ZHUANGHAO=”?S： Prog \"SUB3\"： N≥2=>Goto 2： Abs(S-O)→W：”BIANZHU=”?→Z:Prog \"SUB1\"：Prog\"SUB4\"：“H=”H◢ \"XS=”:X→X◢\"YS=”:Y→Y◢\"FS=”:(F-90) →F:F►DMS◢ Goto1 Lbl 2：”CX=”?X: X→I：”CY=”?Y: Y→J：”DMG=”?→Z[7]：Prog \"SUB2\"： \"S=\":(O+W) →S◢\"Z=\"Z→Z◢ Z→Z[17] :If N=3:Then Prog”SDYX”: IfEnd ：If N=4:Then Prog”BIANPO”:IfEnd:Goto 2
2.SUB1(正算子程序名)
0.1739274226→A:0.3260725774→B
0.0694318442→K:0.330009482→L 1-L→F:1-K→M G+QEKW(C+KWD)→Z[1] G+QELW(C+LWD)→Z[2] G+QEFW(C+FWD)→Z[3] G+QEMW(C+MWD)→Z[4]
A×cos(Z[1])→X (剩號可省略) X+Bcos(Z[2])→X X+Bcos(Z[3])→X X+Acos(Z[4])→X U+WX→X Asin(Z[1])→Y Y+Bsin(Z[2])→Y Y+Bsin(Z[3])→Y Y+Asin(Z[4])→Y V+WY→Y
G+QEW(C+WD)+90→F X+Zcos(F)→X Y+Zsin(F)→Y
3.SUB2(反算子程序名) G-90→T
(Y-V)cos(T)-(X-U)sin(T)→W Abs(W)→W:0→Z LbI 4:Prog\"SUB1\" T+QEW(C+WD)→L
(J-Y)cos(L)-(I-X)sin(L)→Z
If Abs(Z)＜1E-6:Then 0→Z :Prog\"SUB1\" (J-Y)÷sin(F)→Z Else W+Z→W:Goto 4:IfEnd
注:1E-6即是10的負6次方,亦即是等於0.000001 輸入時按SHIFT---log---(-)---6
4. 数据库(SUB3)
If S<226100:Then Cls:Stop: IfEnd
If S<226255.833:Then 226100→O： 99037.9736→U： 4282.3590→V： 201.3266304→G： 1045→P： 1045→R： 253.543→H： 0→Q: Goto 1:IfEnd If S>And S<„„ „„
If S>236005.651:Then Cls:Stop: IfEnd
Lbl 1: 1÷P→C：(P-R)÷(2HPR) →D：180÷π→E：
5、曲线主程序: 1-E.BG(SUB4) Fix 3:Prog”Z.EBG”: R Abs(D-C)/2→T:
if D-C>θ :Then 1→W:Else -1→W:ifend if S 曲线数据库: Z.EBG(SUB4-1) if S<332.599:Then 0.02793→C:-0.03309→D:1800→R:235→A:32.7→Z: goto 1:ifend if S>322.599 And K<734.085:Then -0.03309→C:-0.00485→D:6000→R:460→A:25.254→Z: goto 1:ifend -- -- Lbl 1 C—前坡:D—后坡:R—呵呵,都知道是什么!A—变坡点:Z—变坡点设计高 O—从设计高程点处垂直向下的距离:N—距设计高程处的平距(若你的设计高程线是路中线的话，请把0.75改为0):J—所求桩号处的横坡. 6. 边坡放样程序(BIANPO) （SUB6） “QPDJL=”? → Z[3]: “QPDGD=” : Prog “（SUB4）”：H→ Z[4]◢“PG=” ? → Z[5]：PTK→Z[6]： “P1=” ? → Z[8]： “P2=” ? → Z[9]： “P3=” ? → Z[10] (没有的话就删)： Z[7]-Z[4] → Z[11]：If Z[11]〈0 Then AbsZ[11] →Z[11]: “T=”：Z[11] ◢Else“W=”：Z[11] ◢IfEn Z[11]÷Z[5] →Z[12]：If Z[12]〈1 Then Z[8]×Z[11] →Z[13]：Else Z[8]×Z[5]+Z[9]×Z[5]×（Z[12]-1）→Z[13]：IfEnd: AbsZ[17]-（Z[13]+Z[6] ×Int (Z[12])+Z[3]）→Z[14]:“BP+C-J=”: Z[14] ◢ 说明： Z[3]= QPDJL „„„„„„ 起坡点离中桩距离 Z[4]=“QPDGD” „„„„„起坡点设计标高 Z[5]=“PG” „„„„„„ 坡高 Z[6]=“PTK” „„„„„„ 平台宽度 Z[7]=“DMG” „„„„„„实测地面高度 Z[8]=“P1” „„„„„„ 一级边坡坡度 Z[9]=“P2” „„„„„„ 二级边坡坡度 Z[10]=“P3” „„„„„„ 三级边坡坡度(没有的话就删) Z[14]=“BP+C-J=” „„„„„„边坡超欠挖宽度（+表超，-表欠） 三、使用说明 1、规定 (1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时， Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。 (2) 当所求点位于中线时，Z=0；当位于中线左铡时，Z取负值；当位于中线中线右侧时，Z取正值。 (3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。 (4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。 (5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。 (6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。 2、输入与显示说明 SZ → XY，2.XY → SZ，3.SDYX，4. BIANPO N ? 选择计算方式，输入1表示进行由里程、边距计算坐标 ； 输入2表示由坐标反算 里程和边距。输入3表示隧道圆心放样和计算实测半径；输入4表示进行路基边坡放样。 4、SUB3数据库说明： S待求桩号，O=曲线起点桩号，U曲线起点X坐标，V曲线起点Y坐标，G曲线起点切线方位角， P曲线起点半径， R曲线终点半径， H=曲线长度，Q转向（直线0，左转-1，右转1） 注：本程序是经过本人改编后的，程序起头3.SDYX可以输入也可不输 5800计算器全线坐标计算放样程序 “XLZBJSCX” ◢ LB1 0 ↙ CLS : FIX 4 : 30→DIM Z ↙ “XHS=\"?G(后视点X):\"YHS=\"?L(后视点Y):\"XZJ=\"?M(置镜点X):\"YZJ=\"?N(置镜点 Y):Pol(G-M,L-N):\"DH=\":I(后视 距)◢J<0=>J+360→J:\"FH=\":J→DMS◢(后视方位角) LB1 1 ↙ “K=”?K ◢（计算里程） IF K<本段曲线终点里程 AND K≥上段曲线终点里程 ：THEN 本段终点里程→Z[1] : 上段曲线终点里程→Z[2] ：1→O （注：左偏曲线输入-1→O,右偏曲线输入1→O）: 偏角→A：半径→R : 第一缓和曲线→Z[6] : 第二缓和曲线→Z[7] : 交点 X→B :交点 Y→C : 小里程向交点方位角→E : 交点向大里程方位角→F : GOTO 2 : IFEND↙ „„„„（曲线段分段输入） 补充直线段输入如下 IF K<本段直线终点里程 AND K≥本段直线起点里程:THEN 1→O:本段直线终点里程→Z[3]:终点X→Z[16]:终点Y→Z[17]:方位角→E:GOTO 4:IFEND LB1 2 ↙（曲线要素计算） Z[6]/2- Z[6]^3/(240*R^2)+ Z[6]^5/(34560*R^4) →Z[8] ↙ （M1） Z[7]/2- Z[7]^3/(240*R^2)+ Z[7]^5/(34560*R^4) →Z[9] ↙ （M2） Z[6]^2/(24*R)- Z[6]^4/(2688*R^3) →Z[10] ↙ （P1） Z[7]^2/(24*R)- Z[7]^4/(2688*R^3) →Z[11] ↙ （P2） π*A*R/180+0.5*( Z[6]+ Z[7])→S ↙ （曲线总长） 90* Z[6]/(R*π) →Z[14] ↙ （第一缓和曲线总偏角） 90* Z[7]/(R*π) →Z[15] ↙（第二缓和曲线总偏角,可以省略） Z[8]＋（R+Z[10])TAN(A/2)-(Z[10]-Z[11] )/SIN A→Z[12]↙ (切线T1) Z[9]＋（R+Z[11])TAN(A/2)+(Z[10]-Z[11] )/SIN A→Z[13]↙ (切线T2) B+ Z[12]*COS (E+180)→ Z[16] ↙ （ZH点X） C+ Z[12]*SIN(E+180)→ Z[17] ↙ （ZH点Y） Z[1]-S→Z[3] ↙ (ZH点里程) Z[3]+ Z[6]→Z[4] ↙ (HY点里程) Z[1]- Z[7]→Z[5] ↙ (YH点里程) GOTO 3 ↙ LB1 3 ↙(判断里程点与曲线关系) IF K≤Z[3] AND K> Z[2] : THEN GOTO 4 : IFEND ↙ IF K≤Z[4] AND K> Z[3] : THEN GOTO 5 : IFEND ↙ IF K≤Z[5] AND K> Z[4] : THEN GOTO 6 : IFEND ↙ IF K≤Z[1] AND K> Z[5] : THEN GOTO 7 : IFEND ↙ LB1 4 ↙（里程小于直缓点直线独立坐标） K- Z[3] →X : 0→Y : E→T : PROG“TYZBCX” ：GOTO 1 ↙ LB1 5 ↙（第一缓和曲线独立坐标） K- Z[3] →H ↙ H-H^5/(40*R^2* Z[6]^2)+H^9/(3456*R^4* Z[6]^4) →X ↙ H^3/(6*R* Z[6])-H^7/(336*R^3* Z[6]^3) →Y ↙ 90*H^2/( R*π* Z[6]) →T ↙ IF O>0 ：THEN T +E→T : ELSE E-T →T : T<0=>360+T→T : IFEND ↙ PROG“TYZBCX” ：GOTO 1 ↙ LB1 6 ↙（圆曲线独立坐标） K- Z[4] →H ↙ H*180/( R*π)+ Z[14]→T ↙ R*SIN T+ Z[8]→X ↙ R*(1-COS T)+ Z[10]→Y ↙ IF O>0 ：THEN T +E→T : ELSE E-T →T : T<0=>360+T→T : IFEND ↙ PROG“TYZBCX” ：GOTO 1 ↙ LB1 7 ↙（第二缓和曲线独立坐标） Z[1] -K →H ↙ H-H^5/(40*R^2* Z[7]^2)+H^9/(3456*R^4* Z[7]^4) →U↙ H^3/(6*R* Z[7])-H^7/(336*R^3* Z[7]^3) →V ↙ 90*H^2/( R*π* Z[7]) →T ↙ Z[13]COS A+ Z[12]-U*COS A-V*SIN A→X ↙ Z[13]*SIN A-U*SIN A+V*COS A→Y ↙ IF O>0 ：THEN F-T→T : T<0=>360+T→T : ELSE F+T →T : IFEND ↙ PROG“TYZBCX” ：GOTO 1 ↙ 子程序： “TYZBCX” ↙ （统一坐标计算） IF O<0 : THEN -Y→Y : IFEND ↙ “QXJ=” :T◢ （计算里程点切线方位角，可以不显示） Z[16]+X*COS E-Y*SIN E→Z[18] ↙ Z[17]+X*SIN E+Y*COS E→Z[19] ↙ “XI=” : Z[18] ◢ （中线X） “YI=” : Z[19] ◢ （中线Y） Pol(Z[18]-M,Z[19]-N):\"DI=\":I （ 中 桩 放 样 距 ） ◢J<0=>J+360→J:\"FI=\":J→DMS◢（中桩放样方位角） “PJ=”?P◢ (输入边桩与线路夹角，左-右+) “PD=”?D◢ （输入边桩距） Z[18]+D*COS(T+P) →Z[20] ↙ Z[19]+D*SIN(T+P) →Z[21] ↙ “XP=”: Z[20] ◢ （边桩X） “YP=”: Z[21] ◢ （边桩Y） Pol(Z[20]-M,Z[21]-N):\"DP=\":I◢ （ 边 桩 放样 距 ） J<0=>J+360→J:\"FP=\":J→DMS◢（边桩放样方位角） RETURN↙ 经演算可以放心使用 因篇幅问题不能全部显示，请点此查看更多更全内容