1÷P→C:(P-R)÷(2HPR) →D: 180÷π→E←┚
0.1739274226→A: 0.3260725774→B: 0.0694318442→滑蠢K: 0.3300094782→L: 1-L→F: 1-K→M:U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QEFW(C+FWD))+Acos(G+QEMW(C+MWD)))→X: V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEFW(C+FWD))+Asin(G+QEMW(C+MWD)))→Y: G+QEW(C+WD)+90→F: X+ZcosF→X: Y+ZsinF→Y
3反算子程序(FS)(不运行)
G-90→T: Abs((Y-V))cos(T)-(X-U)sin(T)→W: 0→Z←┚
Lbl 0←┚
prog"ZS":T+QEW(C+WD)→L:(J-Y)cosL-(I-X)sinL→Z:If Abs(Z)<0.001: Then Goto1:Else W+Z→W:IfEnd :Goto 0←┚
Lbl 1: 0→Z:Prog"ZS":(J-Y)÷sinF→Z
ZX-WN(运行主程序,该程序适用于匝道,对全段完整非对称曲线计算太过麻烦)
"CZX"?S(仪站X坐标):"CZY"?W(仪站Y坐标)
LbI 0: "K"?K: If K≤126.9(第1条单一曲线的终点下面同上): Then Goto A:Else If K≤176.9 And K>126.9: Then Goto B:Else If K≤259.562 And K>176.9: Then Goto C:Else If K≤309.562 And K>259.562: Then Goto D: IfEnd:IfEnd:IfEnd:IfEnd←┚
LbI A:92094.336(起点X坐标)→U:8287.962(起点Y坐标)→V:0起点桩号→O:145°11'54"(方位角)→G:126.9(曲线长度)→H:10^(45) (起点半径)→P:10^(45) (止点半径)→R:+1(路线左右偏,左负右正,直线输0)→Q:Prog "WN":Goto W
LbI B:91990.134→U:8360.388→V:126.9→O:145°11'54"→G:50→H:10^(45) →P:150 →R:+1→Q:Prog "WN":Goto W
LbI C:91947.609→U:8386.570→V:176.9→O:154°44'51"→G:82.662→H:150 →P:150 →R:+1→Q:Prog "WN":Goto W
LbI D:91867.101→U:8399.992→V:259.562→O:186°19'19"→G:50→H:150 →P: 10^(45) →R:+1→Q:Prog "WN":Goto W
LbI W:"PJ1"?T←┚(该项是只左右边桩再左右偏距,输0为不再左右偏,可以直接出方位角、距离,适用于桥梁桥台、盖梁、涵洞八字墙等放样,以路线前进方向,前-,后+)
If T=0:Then Z[3]→Z: Z[4]→M: Goto P: Else Goto V: IfEnd←核激┚
LbI V:"JJ1"? O←┚(左右边桩再偏夹角,用180°减去路线夹角)
"X1=": Z[3]+Tcos(F+N+ O)→U◢
"Y1=" :Z[4]+Tsin(F+N+ O)→V◢
U→Z: V→M: GotoP←┚
LbI P:Pol(Z-S,M-W): If J≤0:Then J+360°→J: IfEnd←┚
"JL=":I◢(仪站至测点的距离)
"FWJ=":J DMS◢(仪站至测点的方位角)
Goto0
备注:
1. 规定以道路中线前进方向(即里程增大的方向)区分左右;当线元往左偏时Q=-1;当元线往右偏时,Q=1:当元线为直线时,Q=0。
2. 当线元为信氏陪直线时,其起点、止点的曲率半径为无穷大,以10的45次方代替。
3. 当线元为圆曲线时,无论其起点、止点与什么线元相接,其曲率半径均等于
坐标正反算程序
正算(ZS1)
Lbl0:"DK="?K:"W="?Z:"A="?B↙
Prog"data1"↙
"X=":N+Zcos(F+B)◢
"Y=":E+Zsin(F+B)◢
"FWJ=":F◢
Goto 0
K=?:待求桩号
W=?:偏距(左负,右正)
A=?:偏角
反算(FS1)
"X="?C:"Y="?D:"DK="?K↙
Lbl 0:Prog"data1"↙
Pol(C-N,D-E)↙
ListAns[1]→I↙
ListAns[2]→J↙
Icos(F-J)→S:K+S→K↙
Abs(S)>0.0001=>Goto0↙
"DK=":K◢
"W=":Isin(J-F)→Z◢
X?=:反算X坐标
Y?=:反算Y坐标
DK=:近似里程(必须在线元内)
W=:偏距(左负右正)
3.子程序(XY1)
K-A→S:(Q-P)÷L→I↙
N+∫(cos(F+X(2P+XI)×90÷π),0,S)→N↙
E+∫(sin(F+X(2P+XI)×90÷π),0,S)→E↙
F+S(2P+SI)×90÷π→F
4.数据库(Data1)
K≤7586.707(判敏线元起点)=>Stop↙
7586.707→(起点桩号)A:3378605.445→(起点X)N:453648.704→(起点Y)E:98°56′55.62″→(起渣野点方位角)F:0→(起点曲率半径、转角左负右正)P:1/4500→(止点曲率半掘梁枝径、转角左负右正)Q:360→(线元长度)L:K≤A+L=>GoTo 1↙
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)