程序清单

AA

Do:Deg:ClrText

“F1 Fangyang”

“F2 Jiance”

“F3 suidao”

‘FN:6→V:3→I:Prog “GN”

‘XL:3→V:Prog “GN”

Prog “FN”:ClrText

“Cont”:1→W:Prog “SC”

LpWhile F=1

ZZ

List 7[1]→I:List 8[I]→I:1→V

Do:I+1→I:List 9[I+1]→Q

LpWhile K>Q-List 12[I+1]

List 10[I]→M:List 11[I]→N

List 16[I]→S:List 17[I]→R

List 9[I]→Q:-List 12[I]→J

6→W

‘QH:Q+J→Q:Prog “SC”:Prog “GN”

K-Q→J:1E-20→B:List 13[I]→O

O<0=>1÷List 17[I-1]→B

AbsO→O:(1÷R-B)÷2÷O→C

ZX:If J<0:Then Prog “SC”

Return :IfEnd

7→W

‘HQ:If J<O:Then Prog “SC”

Return:IfEnd

‘HY:O→J:Prog “SC”:Prog “GN”

Q+O→Q:K-Q→J:List 14[I]→O

8→W

‘YQ:If J<O:Then Prog “SC”

Return:IfEnd

‘YH:O→J:Prog “SC”:Prog “GN”

Q+O→Q:K-Q→J:List 15[I]→O

1÷R→B: 1E20→R

O<0=>List 17[I+1]→R

Abs O→O: (1÷R-B)÷2÷O→C

7→W

‘HQ:If J<O:Then Prog “SC”

Return:IfEnd

‘HZ:O→J: Prog “SC”:Prog “GN”

K-Q-O→J

6→W

‘ZX:Prog “SC”

主程序

设置角度单位为度

F1放样

F2 检测

F3 隧道

接受F1~F3功能键

选择输入线路号

执行上述按键对应的F1~F3功能

是否继续F1是 F2 否

按键为F1时重新开始

中桩计算子程序

获取平曲线数据行号,指定沿曲线按导线递推

读取下一曲线交点里程JD

直到里程小于下一曲线起点(JD-T)

读取交点坐标

读取切线方位角和圆半径

读取交点里程和负的切线长

直线计算功能

由交点推算到起缓

求第一缓曲测线长、默认起点曲率、缓曲长

若第一缓长为负令前曲圆曲率为起点曲率

第一缓长改为正,求缓曲变化系数

K小于曲线起点按直线计算

跳出子程序

缓曲计算功能

K小于缓圆按缓曲计算

跳出子程序

由曲线起点推算到缓圆

求圆曲段测线长及圆曲长

圆曲计算功能

K小于圆缓按圆曲计算

跳出子程序

由缓圆推算到圆缓

求第二缓曲测线长、缓曲长

求起点曲率、默认终点半径

若第二缓长为负令后曲圆半径为终点半径

第二缓长改为正,求缓曲变化系数

缓曲计算功能

K小于曲线终点按缓曲计算

跳出子程序

由圆缓推算到缓终

求直线长

直线计算功能

K大于曲终按直线计算

FN

‘FY:If:F=1:Then Do:ClrText

“Fangyang”

5→V:Prog “GN”:I=0=>Return

Prog “ZZ”:Prog “GN”

Do:ClrText:4→W:0→L

“F1 Daoxian”:“F2 Shexian”:“F3 2Bian”

‘FN:6→V:3→I:Prog “GN”:F→V:F=3=>2→V

“J”?->r:r→J: “A”→θ: θ→A

‘SX.DX:ClrText:Prog “GN”

‘ “XY”:Prog “SC”

‘SX:If F=3:Then -J→r:A→θ

Prog “GN”:Prog “SC”:IfEnd: “OK”▲

ClrText:“F1 Next”: “F2 K?”: “F3 Quit”

‘FN:6→V:3→I:Prog “GN”

LpWhile F=1:LpWhile F=2

Return:IfEnd

‘JC:If F=2:Then Do:ClrText

“Jiance”:0→List 25[1]

“F1 Lumian”:“F2 Bianpo”:“F3 2D”

‘FN:6→V:3→I:Prog “GN”:F→U

‘K?:5→V:Prog “GN”:I=0=>Return

Do:‘FS: “X”?→X:“Y”?→Y

7→V:Prog “GN”: I=0=>Return

‘ “KJ”:ClrText:0→L:5→W:Prog “SC”

If U≠3:Then “Z”?→Z

List 7[1]→I:List 8[I]→A

If List 25[1]=0:Then‘ZG:8→V:Prog “GN”

3→W:Prog “SC”:‘DM:0→W:Prog “SC”

‘LJ:9→W:Prog “SC”:IfEnd

If U=1:Then “-H”?→R:Else Z-G→R:IfEnd

‘PK:0→V:Prog “GN”

If U=1:Then AbsJ→r:Z-R-H-(r-K)C→R:r-K-Q→J

Else Abs J-E-Q→J:List 25[1]→Q:J→List 25[1]

Q≠0=>QJ÷(Q-J)→J:IfEnd

‘ “dKdH”:0→L:2→W:Prog “SC”: “OK” ▲

If U=2:Then ClrText

“Next PJ”:1→W:Prog “SC”:F=1=>0=>List 25[1]

IfEnd:IfEnd:ClrText

“F1 Next”: “F2 Menu”: “F3 Quit”

‘FN:6→V:3→I:Prog “GN”

LpWhile F=1:LpWhile F=2

Return:IfEnd

主功能子程序

放样计算

“放样”

输入里程

中桩坐标计算

指定显示坐标模式

F1 导线,F2 射线,F3 二边

接受F1~F3功能键

读取测线长,转角

执行导线或射线计算

显示坐标

如按键是F3再算对称点

计算,显示图形第二点坐标

F1 下一点,F2 更换中桩,F3 退出放样

接受F1~F3功能键

执行上述按键对应的F1~F2功能

如是F3 退出

检测计算

“检测”

F1 路面,F2 边坡,F3 2维”

接受F1~F3功能键

输入近似里程

输入待检坐标

执行反算求出里程和位置

显示里程和位置

如不是F3则输入测点高程

获取纵横断面数据行号

如第一次测断面则计算中桩高

显示中桩高,获取断面数据行号

计算硬(土)路肩高,宽,

如F1待检层距路面高,F2 计算测点和土路肩高差

求测点边坡宽(相对路肩)

如F1求相对路面边的高宽超欠

如F2 求相对土路肩高差和边坡宽基本超欠

如不是第一点则考虑原地面坡度修正宽超欠

显示高差,宽度超欠

如F2 检测边坡时询问 “下一坡脚 F1是 F2否”

F1则将上次测点距边坡归零

如是F3不执行输入测点高及以上语句

F1 继续原来功能 F2 返回功能选单F3退出检测

接受F1~F3按键

执行上述按键对应的F1~F2功能

如是F3退出

GN

If V=1 Or V=2:Then S+θ→θ:{M,N}+Rec(r, θ)

List Ans[1]→r:List Ans[2]→θ

If V=1:Then r→M: θ→N:S+A→S:IfEnd

Return:IfEnd

If V=3:Then Dim List 6→Q:Q=1=>Stop

Do:ClrText: “”:?→I:Int Abs I+1→I

LpWhile I<2 Or I>Q:I→List 7[I]

Return:IfEnd

If V=4 Or V=5:Then List 7[1]→I

List 6[I]→Q:List 7[I]→M:If V=5:Then ClrText

“XL”: “QDK”: “ZDK”:Locate 4,1,I-1

Locate4,2,Q:Locate4,3,M: “K”?→K:IfEnd

If K<Q Or K>M:Then 0→I: “K over” ▲

Else List 8[I]→I:IfEnd:Return:IfEnd

If V=6:Then Do:Do:8.9-GetKey÷10→F

LpWhile Frac F≠0:LpWhile F>I: Return:IfEnd

If V=7:Then Do:Prog “ZZ”:Prog “GN”

Pol(X-M,Y-N):List Ans[1]→M

Rec(M,List Ans[2]-S:List Ans[2]→J

P*list Ans[1]÷(P-J)→Q:K+Q→K

4→V:Prog “GN”:I=0=>Return

LpWhile Abs Q>1milli:Return:IfEnd

If V=8=:Then List 9[A]→θ

If θ<0:Then ClrText:“ZG”:1→W:Prog “SC”

If F=1:Then –θ→θ:Return:Else 0→θ:IfEnd:IfEnd

If θ=0:Then “ZG”?→θ: -θ→List 9[A] :Return: :IfEnd

θ→I:Do:List 18[I+1]→P:List 21[I]→N

List 19[I]→Q:List 18[I]→O:I+1→I:P-O→P

NAbs P÷P→N:Abs(NP÷2)→P

LpWhile K>Q+P:K-Q→Q:List 20[I]+OQ→θ

Q+P→Q:Q>0=>θ+Q2÷2÷N→θ:Return:IfEnd

If V=9:Then If r<0:Then Do:List 22[I]→P

List 23[I]→Q:I+1→I:List 22[I]→N:List 23[I]→O

LpWhile K>N:IfEnd

Else Do:List 22[I]→P:List 24[I]→Q:I+1→I

List 22[I]→N:List 24[I]→O:LpWhile K>N:IfEnd

(K-P)÷(N-P)→P:3P^(2^S)-2P^(3^S)→P

Q+P(O-Q)→I:Return:IfEnd

If V=0:Then D→I: Do:List 25[I]→P:List 26[I]→Q

I+1→I: List 25[I]→N:List 26[I]→O

LpWhile R<Q Or R>O

(R-P) ÷(N-P)→P: Q+P(O-Q)→Q:Return:IfEnd

功能子程序

计算下一点坐标

射线计算

如果V=1把终点作为起算点(导线)

输线路号,无数据时退出

输入线路号

如有数据,线路所在行号保存到List 7[1]

检查里程

取得起点里程 终点里程

线路起终点里程

显示已知数据,提示输入里程

如果里程不在范围则 “K超范围”

否则输出平曲线数据所在行号

识别按键F1~FI

只对F1~FI有效

反算,先算出中桩坐标

求出测点到中桩距离

求出到切线垂直距离

计算切线方向改正量

检验里程是否超限

切线改正量小于1mm则结束

计算中桩高程

如果上次高程为手输,则询问是保留F1 是F2 否

F1用原数据,F2 重输入

手工输入中桩高

自动计算中桩高,读出后坡,半径

变坡里程,前坡,求坡差

求半径和切线长

如果在竖曲线内则计算直坡各点高

在竖曲线上的加上改正量

直线或3次抛物线内插求值

查左线数据库,

查右线数据库

求出比例系数

内插求值

边坡宽直线内插求值

查询变坡点

找到高程点所在坡面

直线内插求出边坡宽

SC

If W=1:Then “F1 Yes”: “F2 No”

6→V:2→I:Prog “GN” :Return:IfEnd

If W=2:Then “dG=”:”dK=”:L+2→L

Locate 4,L-1,R:Locate 4,L,J:Return:IfEnd

If W=3:Then “ZG=”:L+1→L

Locate 4,L,θ: Return:IfEnd

If W=4:Then “X=”:”Y=”:L+2→L

Locate 3,L-1,r:Locate 3,L, θ:Return:IfEnd

If W=5:Then “K=”:”J=”:L+2→L

Locate 3,L-1,K:Locate 3,L, J:Return:IfEnd

If W=6:Then 0→A:0→θ:J→r:1E20→P:Return:IfEnd

If W=7:Then JrList 2(CJList 2+B)→List 3

Jlist 1cosList 3:Sum List Ans→r

Jlist 1sinList 3:Sum List Ans→θ

Jr(CJ+B)→A:Pol(r, θ):List Ans[1]→r

List Ans[2]→θ:1÷(2CJ+B)→P: Return:IfEnd

If W=8:Then Jr÷R→A:A÷2→θ

2Rsinθ→r:R→P: Return:IfEnd

If W=9:Then J→r:List 10[A]→M

List 14[A]→R:List 11[A]→T

List 12[A]→C:List 15[A]-1→S

If S=-2:Then ClrText

“Next Lm”:1→W:Prog “SC”:F-3→S:IfEnd

If S=-1:Then S→List 15[A]

“CG”?→C:C→List 12[A]

“KD”?→M:M→List 10[A]

“TJ”?→T:T→List 11[A]

“CZ”?→R:R→List 14[A]

Else If S=0 Or S=1:Then 7→V:C→I

Prog “GN”:I→C:T→I:Prog “GN”:I→T

M→I:Prog “GN”:I→M:R→I: Prog “GN”:I→R

If R<0:Then If Abs C>0.02:Then –JC→r

M→I:Prog “GN”:I→R:R(Abs C-0.02)+θ→θ:IfEnd

0→R:IfEnd:θ+MC→H

If C<-0.04:Then H+TC→θ:Else H-0.04T→θ:IfEnd

R+M→K:K+T→r:Return:IfEnd

If W=0:Then List 13[A]→I

If I<0:Then 2→D:If I=-1:Then ClrText

“New Dm”:1→W:Prog “SC”:F-2→I:IfEnd

If I=0:Then “H”?→List 25[D]:“W”? →List 26[D]

Do:D+1→D:“H”?→List 25[D]:“W”? →List 26[D]

ClrText: “Next Bp”:1→W:Prog “SC”:LpWhile F=1

2→D:-1→List 13[A]:IfEnd:Return:IfEnd

7→V:Prog “GN”:O→D:Return:IfEnd

输出等子程序

是否选项 F1 是 F2 否

接受F1, F2键

高度和宽度误差显示

显示数值

高程输出显示

显示数值

坐标输出显示

显示数值

位置输出显示

显示数值

直线延长计算,转角,偏角,弦长,半径

缓曲延长计算,分偏角计算

坐标分量计算

坐标分量计算

转角,弦长

偏角,半径

圆曲延长计算,转角,偏角

弦长,半径

路肩高宽计算,路面保存位置

分隔带,土路肩宽保存位置

超高保存位置,线路等级(内插方式)或数据写入状况

如上次数据是手工临时输入

提示是用原数据还是重输

如要求重输,则先改写数据写入状态为已输入

输入超高

输入路面宽

输入土路肩宽

输入分隔带半宽

如已知线路等级则按内插方式计算

C超高,T,土路肩宽

M路面宽,分隔带宽

如分隔带宽小于0,且超高大于2%

计算曲线内侧路面宽,改正中桩高

计算路面边高

以-4%坡为界,分别计算土路肩高

计算路面宽,土路肩宽

断面号查询

如临时断面数据已输入

提示用原断面还是重输入

如要求重输,先输入第一个边坡点

输入第二个边坡点

提示是否还输入下一点

将数据状态改为已输入,输出2为断面数据位置

否则查表找出断面数据位置.

使用说明书

• 数据准备

1. 高斯公式常数表

高斯公式五节点10个常数保存在LIST第1、2列，共占5行。如要提高运算速度或提高结果精度可改变此两列数据。

1. 线路数据格式

线路数据保存在LIST的6~8列，第一行第第7列保存运行时线路选择结果。

从第二行开始保存各条线路基本情况，第6列为起点里程，第7列为终点里程，第8列为对应线路数据的详细资料的起始行号，数据必须连续存放。

1. 平曲线数据格式

各条线路对应的平曲线数据分段保存在LIST的9-17列，各条线路保存顺序不限，每条线路数据连续存放。

每段的第一行各列为横断面有关数据保存的起始位置或计算参数。第9列为竖曲线，第10列为路面宽度，第11列为土路肩宽度，第12列为超高，第13列为线路等级，第14列为超高旋转方式，第15列为边坡设计形式。

第一行第17列为第一个交点第一缓曲为不完整曲线时的起始半径，如为0表示起点是直线。

第二行起为交点资料，第9列为交点里程，第10列为交点X坐标，第11列为交点Y坐标，第12列为第一缓曲段长度（不完整时长度为负），第13列为圆曲线段长度，第14列为第二缓曲段长度（不完整时长度为负），第15列为圆曲线半径（右转为正），第16列为起点切线方位角，第17列为第一切线长度。

最后一行第9列是终点里程，第15列为最后一个交点第二缓曲为不完整曲线时的终点半径，如为0表示终点是直线，第17列必须为0。

此种数据格式可适应主线交点法及匝道线元法设计资料，主线按设计资料对应输入即可，匝道则先组合成不同的交点，每个交点为直、缓、圆、缓、直的组合，但每段长度都可为0，缓和曲线可为不完整的，令交点里程为第一缓曲起点里程，切线长为0。有一段匝道如下线元组合，线元长略。

第1段：直线

第2段：缓曲

第3段：圆R=200

第4段：缓曲

第5段：圆R=100

第6段：缓曲

第7段：直线

可以分成1、2、3、4和5、6、7两个交点也可分成1、2、3和4、5、6、7两个交点，只要算出每个交点第一个曲线元的起始坐标、方位角就行。连续的数据库中把不完整的缓和曲线线元4的缓曲长改成负的程序就会自动处理。

1. 竖曲线数据格式

各条线路对应的竖曲线数据分段保存在LIST的18-21列，各条线路保存顺序不限，每条线路数据连续存放。

第一行起为竖曲线资料，第18列为前坡度（上坡为正），第19列为变坡点里程，第20列为变坡点高程，第21列为竖曲线半径（下凹为正）。

最后一行18列为后坡度。

1. 横断面数据格式

各条线路对应的路面宽度、路肩宽度、路面超高（向路肩下坡为正）、边坡形式数据（简称数据）分段保存在LIST的22-24列，各条线路保存顺序不限，每条线路数据连续存放。

各行的22列为变化点里程，23列为左侧数据，24列为右侧数据。

1. 边坡详细数据

各种边坡变坡点对应的距路肩宽度、高度数据分段保存在LIST的25-26列，各种边坡形式保存顺序不限，每种边坡数据连续存放。

每种边坡形式数据的第25列为变坡点距路肩高差，第26列为变坡点距路肩宽度。

• 修改程序

为了更方便的使用，GN子程序的线路选择模块暂写了“”，双引号里在正式程序里用序号和线路名称的对照表代替，大致格式如下，每组双引号里内容不超过20字符，每7行加一▲

例如：

“1 ：ZHUXIAN”：“2 ：A-ZADAO”：“3 ：B-ZADAO”：“4 ：C-ZADAO”

“5 ：D-ZADAO”：“6 ：E-ZADAO”：“7 ：F-ZADAO”▲

“8 ：G-ZADAO”：“9 ：H-ZADAO”：“10 ：I-ZADAO”

• 程序使用