竖曲线高程计算程序

“SQXJSCX”↙
Lb1 1↙
CLs:Fix 3↙
“K=”?k◢ (计算点里程输入)
If k<67549.755 AND K≥66894.3 :Then -0.00052→A : 0→B : 67394.3→S : 67.37→G : 600000→R : Goto 2 : IfEnd↙
If k<68708.391 AND K≥67549.755 :Then -0.0048→A : -0.00052→B : 68494.3→S :66.8→G : 100000→R : Goto 2 : IfEnd↙
If k<本段竖曲线终点里程 AND K≥前一竖曲线终点里程 :Then -0.0048→前坡（大里程向）A : -0.00052→后坡（小里程向）B : 68494.3→竖曲线交点里程S :66.8→交点高程G : 100000→曲率半径R : Goto 2 : IfEnd↙
………
依次类推，计算原始数据完成输入，坡度换算成小数。

Lb1 2↙
S-K→L ： A-B→W : Abs（R*W/2）→T : S-T→C : S+T→D ↙
If K≤C : Then G-L*B →H : Goto 5 : Else If K≤S : Then Goto3 : Else If K≤D : Then Goto4 : IfEnd : IfEnd : IfEnd↙
Lb1 3↙
If W>0 : Then G+(C-K)^ 2/2/R-L*B→H : Goto 5 : Else If W<0 : Then G-(C-K) ^2/2/R-L*B→H : Goto 5 : IfEnd : IfEnd↙

Lb1 4↙
If W>0 : Then G+(K-D) ^2/2/R-L*A→H : Goto 5 : Else If W<0 : Then G-(K-D)^ 2/2/R-L*A→H : Goto 5 : IfEnd : IfEnd↙
1 5↙
H-0.304→H↙ (0.304为中心~路肩高差)
“H=”：H◢
Goto 1↙
Norm 1:50→C:12345→J
"RESET PW"?I:I=J=>500→DimZ 为数据库增加额外变量500个，在SET、SETPFDYS程序根据实际再增减变量
50→Z[C+22]
"PASSWORDS"?I:I→Z[C+39] 重设要素保护密码
Cls:Stop
FileName:DATLOCK 要素保护密码确定认程序
Cls:Norm 1:50→C
"PASSWORDS"?J:Cls:J≠Z[C+39]=>Stop 2008-08-07修改
FileName:SHELL(外壳程序)
50→C ;在扩充变量预留前50个给别的程序用。如不够就适量加大。RESET，SHELL，SET，SETPFDYS，这几个程序中C值必需一致
Z[C+35]→I:"STATION-N"?I:I→Z[C+35] 设置测站N坐标
Z[C+36]→I:"STATION-E"?I:I→Z[C+36] 设置测站E坐标
Z[C+44]→I:"STATION-Z"?I:I→Z[C+44] 设置测站Z坐标
Z[C+45]→I:"STATION-HI"?I:I→Z[C+45] 设置仪高
Z[C+41]→I:"GC-DH"?I:I→Z[C+41] 放样点高差常数
Norm 1
Z[C+23]→N　当前分段要素N坐标
Z[C+24]→E　当前分段要素E坐标
Z[C+25]→M　当前分段要素起点桩号
Z[C+26]→H　当前分段要素起点方位角（正北）单位：弧度
Z[C+27]→A　当前分段要素起点曲率有左偏负右偏正（注意不半径）
Z[C+28]→R　当前分段要素终点曲率有左偏负右偏正（注意不半径）
Z[C+29]→L　当前分段要素长度
Z[C+3]→W
Z[C+32]→I:"DAT1 2 3"?I:I→Z[C+32] 平曲线要素数库类型选择1为内置式，2文件式，3实时输入(查看当前要素值)
0→I
Z[C+4]→I:"1 2 3"?I:I→Z[C+4] 选择本程序模式默认为坐标正算，1坐标反算，2横向边仰坡放样，3，纵向边仰坡放样(隧道进出口用到)
Lbi 1
Z[C+4]→I
I=1=>Prog"TURNZH"
I=2=>Prog"HXBYP"
I=3=>Prog"ZXBYP"
Lbi A

Deg:Norm 1:Cls
Z[C+1]→G:"ZH"?G:G→Z[C+1] ;G桩求桩号
Z[C+2]→B:"JL"?B:B→Z[C+2] ;B横向距离，左正右负
Lbi B
Prog "ZBJS" ;坐标正算
Prog "GCJS" ;高程计算
Fix 3:Cls ;设置三位小数
"ZH＝":Locate 4,1,G ;第一行显示桩号
"X＝":Locate 3,2,X ;第二行显示X(N)坐标
"Y＝":Locate 3,3,Y ；第三行显示Y(E)坐标
Prog "GCJS" ;高程计算
"Z":Locate 2,4,Z+Z[C+41] ；第四行显示Z坐标
Locate 10,4,Z[C+2] ；第四行显示横向距离
0→I
Do:I+1→I:I=1000=>Goto 1:LpWhile Not(Getkey=57 Or Getkey=27) ;锁定键盘，并在几十秒后自动返回
Getkey=57=> Goto 1
Z+Z[C+41]+Z[C+49]-Z[C+44]-Z[C+45]→Z[C+48]
Pol(Z[C+35]-X,Z[C+36]-Y)
Cls:"DL":Locate 3,1,I:Locate 12,1,Z[C+45]
"DH":Locate 3,2,Z[C+48]:Locate 12,2,Z[C+49]
"FWJ"
J<0=>J+360→J:J◤DMS◢ 2008-08-07 修改
J→Z[C+46]
I→Z[C+47]
Goto 1
FileName：ZBJS 坐标计算程序
Prog"READDAT"
Rad
G-M→Q
IF AR=0 :Then If A=R :Then 1→J:Else 3→J :IfEnd:Else If A=R :Then 2→J:Else 3→J:IfEnd:IfEnd
If J=1 :Then H→F:N+QCos(H)→X:E+QSin(H)→Y:IfEnd 直线段直接计算
If J=2 :Then H+QR→F:Rec(Abs(2Sin(Abs(0.5QR))÷R),H+QR÷2):N+I→X:E+J→Y:IfEnd 圆弧段直接计算
If J=3 :Then Goto 5 :IfEnd
X+BCos(F+W)→X
Y+BSin(F+W)→Y
Deg:Return
Lbi 5 用五点通用坐标计算计算缓和段
0.5(R-A)÷L→K
AQ→I
KQ2→J
0.0469100770→P:H+IP+JP2→U
0.2307653449→P:H+IP+JP2→V
0.5→P :H+IP+JP2→D
0.7692346551→P:H+IP+JP2→F
0.9530899230→P:H+IP+JP2→T
0.1184634425→I
0.2393143352→J
0.2844444444→O
N+Q(ICos(U)+JCos(V)+OCos(D)+JCos(F)+ICos(T))→X
E+Q(ISin(U)+JSin(V)+OSin(D)+JSin(F)+ISin(T))→Y
H+AQ+KQ2→F
X+BCos(F+W)→X
Y+BSin(F+W)→Y
Deg:Return
FileName:TURNZH(坐标反算)
Norm 1
Z[C+37]→I:"DQD-N"?I:I→Z[C+37] 输入待求点N坐标
I=-1=>Prog"INFWJDLDH"
Z[C+38]→I:"DQD-E"?I:I→Z[C+38] 输入待求点E坐标
0→B：M+L÷2→G：Prog"ZBJS"
Lbi S:Rad
Z[C+37]-X→I:Z[C+38]-Y→J
If I=0 And J=0 :Then Goto A:IfEnd
Pol(I,J)
Lbi A:Rec(I,J-F):J→Z[C+2]
G+I→G:IF Abs(I)>0.0001 :Then Prog"ZBJS":Goto S↙
G→Z[C+1]
FileName：INFWJDLDH 实测坐标（用方位角，距离，高差）输入程序
Deg:Norm 1:Cls
Z[C+46]→J:"DQD-FWJ"?J:J→Z[C+46]
Z[C+47]→I:"DQD-DL"?I:I→Z[C+47]
Z[C+48]→K:"DQD-DZ"?K:K→Z[C+48]
Z[C+49]→P:"RHT"?P:P→Z[C+49]
Rec(I,J)
Z[C+35]+I→Z[C+37]
Z[C+36]+J→Z[C+38]
Z[C+44]+Z[C+45]+K-P→Z[C+43]
FileName:HXBYP 横向边仰坡放样程序
Cls:Norm 1
Z[C+40]→I:"QPD-DL"?I:I→Z[C+40] 输入起坡点与中桩距离常数
Z[C+41]→I:"QPD-DZ"?I:I→Z[C+41] 输入起坡点与中桩高差常数
Z[C+42]→I:"i"?I:I→Z[C+42] 输入边仰坡坡度，左仰坡(路堑)为正，右仰坡(路堑)为负，左边坡为负，右边坡为正
Cls:Z[C+43]→I:"DQD-Z"?I:I→Z[C+43] 输入

实测高程
Prog"TURNZH"
Prog "GCJS"
(Z[C+43]-Z-Z[C+41])Z[C+42]+Z[C+40]→B
Fix 3:B-Z[C+2]◢ 显示与上一次测点的偏差
B→Z[C+2]
FileName:ZXBYP 纵向仰坡放样程序(隧道进口使用)
Cls:Norm 1:Z[C+33]→I:"QPD-ZH"?I:I→Z[C+33] 输入起坡点桩号
Z[C+34]→I:"QPD-Z"?I:I→Z[C+34] 输入起坡点高程
Z[C+42]→I:"i"?I:I→Z[C+42] 仰坡时(隧道进出口仰坡)，进口为正，出口为负)
Cls:Z[C+43]→I:"DQD-Z"?I:I→Z[C+43] 输入实测高程
Prog"TURNZH"
Z[C+42](Z[C+43]-Z[C+34])+Z[C+33]→G
Fix 3:G-Z[C+1]◢ 显示与上一次测点的偏差
G→Z[C+1]
FileName:SET 设置程序
1\Lbi S
2\Norm 1
3\50→C
4\50→Z[C+22]
5\0→K
6\Z[C+5]→I:"1 Or 2"?I:I→Z[C+5] ;选择1线或者2线
7\180Z[C+3]÷π→J:"XZJJ"?J:Jπ÷180→Z[C+3] ;斜桩夹角
8\If I=1:Then Z[C+8]→Z[C+6]:Z[C+9]→Z[C+7]
9\Z[C+12]→Z[C+16]:Z[C+13]→Z[C+17]:Z[C+18]→Z[C+20]
10\Else Z[C+10]→Z[C+6]:Z[C+11]→Z[C+7]
11\Z[C+14]→Z[C+16]:Z[C+15]→Z[C+17]:Z[C+19]→Z[C+20]
12\IfEnd
13\Z[C+7]+3→Z[C+31]
14\0→I:"1PQX,2SQX,3YSDS"?I ;1输入平曲线要素，2输入竖曲线要素，3输入1、2线夹平竖曲线交点(变坡点)数量以让程序规划内存建立要素数据为，其它值退出本程序
15\I≠0=>Prog"DATLOCK" 2008-07-28 18：23修改
16\I=1 => Goto 1:I=2 => Goto 2:I=3 => Goto 3 2008-07-28 18：23修改
17\Prog"SETPFDYS"
18\Stop
19\Lbi 1 ;平曲线要素输入部分
20\Z[C+6]→D
21\-1→Z[D+3]
22\For 1→J To Z[C+16]
23\Cls:"PQX":Locate 9,1,J
24\Z[D+1]→I:"N"?I:I→Z[D+1] ;请输入N坐标
25\Z[D+2]→I:"E"?I:I→Z[D+2] ;请输入E坐标
26\Z[D+3]→I:"R"?I:I→Z[D+3] ;请输入半径R
27\Z[D+4]→I:"LS1"?I:I→Z[D+4] ;请输入缓和曲线1长度
28\Z[D+5]→I:"LS2"?I:I→Z[D+5] ;请输入缓和曲线2长度
29\D+5→D
30\Next
31\1→Z[D-2]
32\Goto S
33\Lbi 2 ;竖曲线要素输入部分
34\Z[C+7]→D
35\-1→Z[D+3]
36\For 1→J To Z[C+17]
37\Cls:"SQX":Locate 9,1,J
38\Z[D+1]→I:"BPDZH"?I:I→Z[D+1] ;请输入变坡点桩号 2008-08-07修改
39\Z[D+2]→I:"BPDGC"?I:I→Z[D+2] ;请输入变坡点高程 2008-08-07修改
40\Z[D+3]→I:"R"?I:I→Z[D+3] ;请输入半径R 2008-08-07修改
41\D+3→D
42\Next
43\1→Z[D]
44\Goto S
45\Lbi 3 ;要素点数输入部分
46\Z[C+12]→I:"1PQXDS"?I:I→Z[C+12] ;输入1线平曲线点数,最小值为3
47\Z[C+13]→I:"1SQXDS"?I:I→Z[C+13] ;输入1线竖曲线点数,最小值为3
48\Z[C+18]→I:"1JDZH"?I:I→Z[C+18] ;输入1线平曲线第一个交点的交点桩号
49\Z[C+14]→I:"2PQXDS"?I:I→Z[C+14] ;输入2线平曲线点数,最小值为3
50\Z[C+15]→I:"2SQXDS"?I:I→Z[C+15] ;输入2线竖曲线点数,最小值为3
51\Z[C+19]→I:"2JDZH"?I:I→Z[C+19] ;输入2线平曲线第一个交点的交点桩号
52\C+

Z[C+22]→Z[C+8] 1线路平曲线要素指针偏移基数
53\Z[C+8]+5Z[C+12]→Z[C+9] 1线路竖曲线要素指针偏移基数
54\Z[C+9]+3Z[C+13]→Z[C+10] 2线路平曲线要素指针偏移基数
55\Z[C+10]+5Z[C+14]→Z[C+11] 2线路竖曲线要素指针偏移基数
56\Z[C+11]+3Z[C+15]→Z[C+21] 要素数据库结束指针
57\Z[C+21]+1→Z[C+30]
58\Z[C+21]→DimZ 增加额外变量
59\-1→K
60\Goto S
FileName:SETPFDYS 交点要素=>分段要素
1\Norm 1
2\50→C
3\Rad
4\Z[C+6]+5→D
5\Z[C+21]+1→F
6\F+6→DimZ
7\-1→Z[F] 设置分段要素开始标志
8\Z[C+20]→V ; 从数据库读入当前线第一个平曲线交点桩号
9\0→G
10\Pol(Z[D+1]-Z[D-4],Z[D+2]-Z[D-3]) ；计算第一条直线的距离I，计算方位角J
11\J→H
12\Z[D-4]→Z[F+1] 直线段
13\Z[D-3]→Z[F+2]
14\V-I→Z[F+3]
15\For 1→G To Z[C+16]-2
16\Pol(Z[D+1]-Z[D-4],Z[D+2]-Z[D-3]) ；计算第本交点与下一交点的距离I，计算方位角J
17\J→H 2008-08-07 修改
18\Pol(Z[D+6]-Z[D+1],Z[D+7]-Z[D+2])
19\I→W:J-H→K ；计算转角K，负值是左转角，正值是右转角 2008-08-07 修改
20\IF Abs(K)>π:Then If K>0:Then K-2π→K:Else K+2π→K:IfEnd:IfEnd
21\Z[D+3]→R ;从数据库读入半径R
22\Z[D+4]→L ;从数据库读入缓和曲线1长度
23\Z[D+5]→M ;从数据库读入缓和曲线2长度
24\L÷2÷R→B ;计算缓和曲线1的B0角
25\M÷2÷R→O ;计算缓和曲线2的B0角
26\(Abs(K)-B-O)R→P ;计算圆曲线长度
27\L2÷24÷R-L^(4) ÷2384÷R^(3) →Q ;内移值计算方法2
28\L÷2-L^(3) ÷240÷R2→S ;切线增长值计算方法2
29\M2÷24÷R-M^(4) ÷2384÷R^(3) →J ;内移值计算方法2
30\M÷2-M^(3) ÷240÷R2→I ;切线增长值计算方法2
31\Abs(π÷2- Abs(K)) →A
32\Rtan(Abs(K÷2))+S+J÷Cos(A)-QTan(A)→T ;计算切线1
33\Rtan(Abs(K÷2))+I+Q÷Cos(A)-JTan(A)→U ;计算切线2
34\H<0=>H+2π→H
35\1→A
36\K<0 =>-1→A
37\H→Z[F+4]
38\0→Z[F+5]
39\0→Z[F+6]
40\F+6→F
41\F+6→DimZ
42\Rec(T,H) 缓和段1
43\Z[D+1]-I→N
44\Z[D+2]-J→E
45\N→Z[F+1]
46\E→Z[F+2]
47\V-T→Z[F+3]
48\L=0=>Goto A
49\H→Z[F+4]
50\0→Z[F+5]
51\A÷R→Z[F+6]
52\F+6→F
53\F+6→DimZ
54\Lbi A 圆弧段
55\0→I:0→J
56\L=0=>Goto B
57\Pol(L-L^(3)÷40÷R2,L2÷6÷R-L^(4)÷336÷R^(4))
58\Lbi B
59\Rec(I,H+AJ)
60\N+I→N
61\E+J→E
62\N→Z[F+1]
63\E→Z[F+2]
64\V-T+L→Z[F+3]
65\H+AB→Z[F+4]
66\A÷R→Z[F+5]
67\A÷R→Z[F+6]
68\F+6→F
69\F+6→DimZ
70\M=0=>Goto C 缓和段2
71\Rec(2RSin(0.5P÷R),H+AB+AP÷R÷2)
72\N+I→N
73\E+J→E
74\N→Z[F+1]
75\E→Z[F+2]
76\V-T+L+P→Z[F+3]
77\H+AB+AP÷R→Z[F+4]
78\A÷R→Z[F+5]
79\0→Z[F+6]
此处删除2008-08-06 修改
80\F+6→F
81\F+6→DimZ
82\Lbi C 此处插入2008-08-06 修改
83

\Rec(U,H+K)
84\Z[D+1]+I→Z[F+1] 此处2008-08-07 修改
85\Z[D+2]+J→Z[F+2] 此处2008-08-07 修改
86\V-T+L+P+M→Z[F+3]
87\V-T+L+P+M-U+W→V 此处2008-08-07 修改
88\D+5→D
89\Next
90\H+K→Z[F+4]
91\0→Z[F+5]
92\0→Z[F+6]
93\F+9→DimZ
94\1→Z[F+7]
95\V→Z[F+9] 要素终点桩号
96\Deg
FileName：READDAT
Z[C+32]→I
I=3=>Prog"READDAT3"
If GM+L :Then ;验证要素是否可用，否则重设要素
I=1=>Prog"READDAT1"
I=2=>Prog"READDAT2"
Prog"REALDAT3"
IfEnd
FileName:READDAT1 平曲线分段要素读取程序1
Z[C+30]→F
0→J
Lbi 1 搜索要素
If G-1=Z[F]=>Goto E:F-6→F:-1→J:Goto 1
Else 1=Z[F+1]=>Goto E:-1=J=>Goto 2:1→J:F+6→F:Goto 1
IfEnd
Lbi 2
F→Z[C+30]
Z[F+1]→Z[C+23]　当前分段要素N坐标
Z[F+2]→Z[C+24]　当前分段要素E坐标
Z[F+3]→Z[C+25]　当前分段要素起点桩号
Z[F+4]→[C+26]　当前分段要素起点方位角（正北）
Z[F+5]→Z[C+27]　当前分段要素起点曲率有左偏负右偏正（注意不半径）
Z[F+6]→Z[C+28] 当前分段要素终点曲率有左偏负右偏正（注意不半径）
Z[F+9]-Z[F+3]→Z[C+29] 当前分段要素长度
Return
Lbi E
Cls
"ERROR"
Stop
READDAT2　文件请输入内容为以下内容以便测试结果与本文一致。π÷180→J
Goto 1
Lbi 1
If G<30008.396:Then Cls:Stop:IfEnd↙ ;判断待求桩号是否小于本要素起点，如小于就停止运行程序
If G<30661.038(终点里程):Then 2851422.090(起点X)→N:513896.440(起点Y)→E:30008.396→M:213°32′52″J(起点方位角)→H:0(起点曲率)→A:0(终点曲率)→R:652.642起点至终点长度→L:Goto E:IfEnd↙
If G<30921.038:Then 2850878.162→N:513535.768→E:30661.038→M:213°32′52″J→H:0→A:1÷1300→R:260→L:Goto E:IfEnd↙
If G<31426.367:Then 2850666.478→N:513385.004→E:30921.038→M:219°16′39″J→H:1÷1300→A:1÷1300→R:505.328→L:Goto E:IfEnd↙
If G<31686.367:Then 2850346.483→N:512998.015→E:31426.367→M:241°32′56″J→H:1÷1300→A:0→R:260→L:Goto E:IfEnd↙If G>31686.367:Then Cls:Stop:IfEnd↙ ;判断待求桩号是否大于本要素终点，如大于就停止运行程序
Lbi E
N→Z[C+23]:E→Z[C+24]:M→Z[C+25]:H→Z[C+26]:A→Z[C+27]:R→Z[C+28]:L→Z[C+29]
FileName:READDAT3 平曲线分段要素读取程序3
Norm 1
Z[C+23]→N:?N:N→Z[C+23] 当前分段要素N坐标
Z[C+24]→E:?E:E→Z[C+24] 当前分段要素E坐标
Z[C+25]→M:"SZH"?M:M→Z[C+25] 当前分段要素起点桩号
180Z[C+26]÷π→H:H◤DMS◢
"FWJ"?H:Hπ÷180→H:H→Z[C+26] 当前分段要素起点方位角（正北）单位：度分秒
Z[C+27]→A:"SR"?A:A→Z[C+27] 当前分段要素起点曲率有左偏负右偏正（注意是半径的倒数）
Z[C+28]→R:"ER"?R:R→Z[C+28] 当前分段要素终点曲率有左偏负右偏正（注意是半径倒数）
Z[C+29]→L:?L:L→Z[C+2

9] 当前分段要素长度
Return
FileName:GCJS　高程计算
Z[C+31]→D
0→J
Lbi S
D→Z[C+31]
Z[D-2]→I
Z[D-1]→V
Z[D+1]→K
Z[D+2]→O
(O-V)÷(K-I)→U
Z[D+4]→I
Z[D+5]→V
(V-O)÷(J-K)→V
U-V→Z ;为正时是凸型竖曲线,负是时为凹竖曲线
Z[D+3]→I ;读入半径
Abs(IZ)÷2→T ;切线长
If GIf J=1 Or Z[D]=-1 :Then UG-UK+O→Z :Return:IfEnd ;计算高程Z
-1→J:D-3→D:Goto S
Else If GG-K+T→Q
Q2÷2÷I→I
Z>0=>-I→I
O-UT+UQ+I→Z ;计算高程Z
Return
Else If GG-K-T→Q
Q2÷2÷I→I
Z>0=>-I→I
O+VT+VQ+I→Z ;计算高程Z
Return
Else If J=-1 Or Z[D+3]=1 :Then
O+VG-VK→Z ;计算高程Z
Return
Else 1→J:D+3→D:Goto S
IfEnd :IfEnd:IfEnd:IfEnd