公路测量程序

一、程序使用流程

本程序数据和主程序是分开的，编程时将不同的工程数据存放到不同的数据文件里，如A匝道，文件名为A，将匝道A所有的曲线线元参数输入A文件里。运行时只要运行文件名A的程序就可以了，具体运行流程见下图： 二、数据文件的编写 （一）交点法数据文件编辑

交点法编写数据文件必须是对称型的，即直线段→缓和曲线段→圆曲线段→缓和曲线段→直线段，（如果任意一端没有直线段，则把直线段长度看做是0），另外圆曲线两侧缓和曲线的旋转常数必须相等，并且和直线段连接处的半径必须是无穷大。

交点法数据文件编写一般是根据设计图纸提供的平面曲线参数一览表提供的参数来编写，每个弯道包括：弯道起点方位角（C），交点X坐标（D），交点Y坐标（E），缓和曲线长度（F，当没有设缓和曲线时，F=0），交点转交（G，向左转弯，G为负值，向右转弯，G取正值），交点桩号（H），弯道圆曲线半径（R）。

下图是一段市政道路设计参数数据。

根据上图提供的数据，可以编辑成如下的数据文件： 文件名:CHLNR

3→DimZ “X0”?A:”Y0”?B:“Ln”?L:Abs(L)-Int(Abs(1000L))/1000→Z[3]:Lbl 0:If Z[3]≠0.0001: Then ?L: Else “Xp”?X:”Yp”?Y:X→Z[1]：Y→Z[2]：IfEnd:Lbl 1:If L>0 :Then 98°39°35.12°→C：4474.384→D：2415.861→E：140→F：31°17°23°→G：410.007→H:600→R:IfEnd: If L>1060 Then 129°56°58.19°→C:4206.421→D：3093.946→E：70→F：-33°50°48°→G：1285.437→H:600→R:IfEnd:Prog”XLJS”:If Z[3]≠0.0001:Then Goto 0:Else (Z[1]－X)cos(Ｏ)＋(Z[2]－Y)sin(Ｏ)→N:L＋N→L:－(Z[1]－X)sin(O)+(Z[2]－Y)cos（O）→K：If Abs（N）≥0.001：Then Goto 1：Else “L=“:L◢ “K=”:K◢

IfEnd:Goto 0: IfEnd

在面程式中，有两个条件转移语句即If L>0：Then 98°39°35.12°→C:4774.384→D: 2415.861→E:140→F:31°17°23°→G:410.007→H:600→R:IfEnd

If L>1060：Then 129°56°58.19°→C:4206.421→D: 3093.946→E:70→F:-33°50°48°→G:1285.437→H:600→R:IfEnd ……

如果还有其他弯道，可以继续完后加。在这些存放设计参数的语句前后的程序表达式是固定的。 说明：

编辑曲线参数时，每个曲线参数放在一个If L>***.***（两个弯道中间直线段上的任意桩号） :Then ***°**°**.**°→C（弯道起点方位角）：****.***→D（交点X坐标）：****.****→E（交点Y坐标）：***.***→F（缓和曲线长度：**°**°**°→G（转角，向左转为负值，向右转为正值）：***.***

→H（交点里程桩号）:***→R（圆曲线半径）:IfEnd 条件式语句里，如果有多个弯道，一直按上述形式编下去，

变量说明： 1、 L>***.***

弯道参数的起点范围，可以是弯道起点，也可以是弯道前直线段里的任意点桩号

2、 ***°**°**.**°→C 弯道起点方位角 3、****.***→D、****.****→E

交点X坐标、Y坐标 4、***.***→F

缓和曲线长度，如果没设置缓和曲线，则缓和曲线长度看做0 5、**°**°**°→G 转角，向左转为负值，向右转为正值 6、***.***→H、交点里程桩号 7、***→R 圆曲线半径

注意，运行程序时，第一部是输入测站坐标X0，Y0，然后输入近似桩号，主要是为了加快收敛，输入时注意当根据桩号和距中距离计算所求点坐标时，请不要把小数点第四位输入0.0001，如K21+369.0921，请把桩号输入21369.092即可，相反，当根据坐标反算桩号和距中距离时，请在输入时把第四位小数数位0.0001，如K21+200，输入时请输入21200.0001。 （二）线元法数据文件编辑

线元法数据文件编辑适合于各种各样的线型组合，特别是在大型立交区或者是不规则的线型组合地方，如小区道路、重丘区山路等有特别明显的有

点，本程序为了减少使用者设计参数的输入，特别把直线与缓和曲线以及圆曲线连接的地方只要输入缓和曲线参数就可以了，这样可以减少输入很多设计参数，同时节约大量的内存。看下图，

该图是一个半喇叭型的互通立交，一共是5条匝道，我们以D匝道为例，来看看线元法数据文件的建立。 下表是D匝道的线元表：

D匝道线元数据分析：从上表中可以看出，D匝道由三个线元组成，第

半径或缓线元线元名起点桩号 偏向 称 线旋(m) 转常数A 缓和曲线 圆曲线 缓和曲线 DK0+000.00 DK0+039.289 DK0+078.504 右偏 右偏 右偏 55 40 65 39.289 494341.902 39.216 105.625 494347.926 494380.677 478027.825 477989. 873 477971.271 260°38′55.7″ 302°18′26.0″ 358°28′45.9″ 线元起点坐标、方位角 和曲长度X(m) Y(m) 方位角 一段是缓和曲线段，根据上表判断，该缓和曲线属于第一段缓和曲线，其最小半径处桩号为J=K0+039.289，该处坐标X=494347.926,Y=477989. 873，半径为R=40米，方位角W=302°18′26.0″，缓和曲线旋转常数P=55，在缓和曲线终点K0+039.289处，有一圆曲线与其相连接，长度I=39.216。第二段是圆曲线，第三段又是一段缓和曲线，根据上表判断，该曲线属于第二缓和曲线，曲线最小半径位于与园曲线连接处，桩号J=K0+78.504，坐标为X=494380.677，Y=477971.271,方位角W=358°28′45.9″，最小半径R=40，缓和曲线旋转常数P=-65（为了判断是第一段还是第二缓和曲线，我们通过P

的正负来判断，第一缓和段P统一用正数表示，第二缓和段统一用负号表示。起点处相连的圆曲线长度I可以看做0，也可以看做是39.216.由此我们整理下。从整个线型全部是向右转向，所以R全部取正数。具体如下： 第一段缓和曲线线元参数：

最小半径处坐标X=494347.926 Y=477989. 873 最小半径R=40

最小半径处方位角：W=302°18′26.0″ 缓和曲线旋转常数P=55 与之相连的圆曲线长度I=39.216 最小半径处桩号J=39.289

此段数据计算范围可以从K0+000~ K0+78.504 第二段缓和曲线线元参数

最小半径处坐标X=494380.677 Y=477971.271 最小半径R=40

最小半径处方位角：W=358°28′45.9″ 缓和曲线旋转常数P=65 与之相连的圆曲线长度I=39.216 最小半径处桩号J= W=358°28′45.9″

此段数据计算范围可以从K0+78.504到匝道结束。

当桩号L位于K0+000~ K0+78.504时，曲线参数按第一段编写，当L>78.504时，按第二段缓和曲线编写，当然第二段缓和曲线计算范围也可以

从圆曲线起点K0+039.289到K0+78.504中间任意桩号开始。下面是D匝道的数据程序文件（文件名“DZD”）：

3→DimZ：X0”?A:”Y0”?B:“Ln”?L: Abs(L)-Int(Abs(1000L))/1000→Z[3]:Lbl 0:If Z[3]≠0.0001: Then ?L: Else “Xp”?X:”Yp”?Y:X→Z[1]：Y→Z[2]：IfEnd:Lbl 1:If L≥0 :Then 494347.926→X：477989. 873→Y：40→R：302°18°26°→W：55→P:39.216→I: 39.289→J：IfEnd: If L≥78.504：Then 494380.677→X：477971.271→Y：40→R：358°28°45.9°→W：-65→P:39.216→I: 78.504→J：IfEnd:Prog”FDCH”: Prog”XLJS”: If Z[3]≠0.0001:Then Goto 0:Else (Z[1]－X)cos(Ｏ)＋(Z[2]－Y)sin(Ｏ)→N:L＋N→L:－(Z[1]－X)sin(O)+(Z[2]－Y)cos（O）→K：If Abs（N）≥0.001：Then Goto 1：Else “L=“:L◢ “K=”:K◢

IfEnd:Goto 0: IfEnd

在上面程式中，每个缓和曲线的线元参数放在一个If … Then …… IfEnd循环语句中，If 后面紧跟计算的起点范围，终点范围从下段的起点范围结束。

变量定义：

1、X、Y 输入缓和曲线半径最小处的X、Y坐标

2、R 缓和曲线最小半径，当线路左转时，半径数据前面加“-”号，以判断线路转向。

3、W 输入缓和曲线半径最小处的方位角

4、P 缓和曲线旋转常数，第一缓和曲线长度为正数，第二缓和曲线常数为负数。当不设缓和曲线时，P=0

5、I 与缓和曲线最小半径端连接的圆曲线长度。 6、J缓和曲线小半径段的桩号。 再看看E匝道数据文件编辑 E匝道线元设计数据表 在E匝道线元表中，我们可以看出，E匝道式由一段直线段→缓和曲线半径或缓线元线元名称 起点桩号 偏和曲长度向 线旋(m) 转常数A 直线 EK0+000.00 线元起点坐标、方位角 X(m) Y(m) 方位角 120 71.296 494360.619 477706.293 477776.470 477847.870 477907.018 477948.932 478010.658 79°49′57.0″ 缓和曲线 EK0+071.296 右偏 72 494373.205 79°49′57.0″ 圆曲线 EK0+143.296 右偏 200 60.070 494381.631 90°08′44.6″ 缓和曲线 EK0+203.366 右偏 95 45.125 494372.527 107°21′16.3″ 缓和曲线 EK0+248.491 左偏 75 64.841 494355.878 113°49′05.6″ 圆曲线 EK0+313.333 左偏 86.75 17.800 494337.373 92°24′19.0″ 段→圆曲线段→缓和曲线段→缓和曲线段→圆曲线段连接起来的。值得注意的是，在位于K0+248.491处，是两段反向的缓和曲线连接点，在编写数据程序文件时，分三段编写，第一段从K0+000到K0+203.366，第二段从

K0+203.366到K0+248.919，第三段从K0+248.919到线路终点K0+331.13。具体文件如下：（文件名“EZD”）

3→DimZ：X0”?A:”Y0”?B:“Ln”?L:Abs(L)-Int(Abs(1000L))/1000→Z[3]:Lbl 0:If Z[3]≠0.0001: Then ?L: Else “Xp”?X:”Yp”?Y:X→Z[1]：Y→Z[2]：IfEnd:Lbl 1:If L≥0 :Then 494381.6316→X：477847.870→Y：200→R：90°08°44.6°→W：120→P:60.07→I: 143.296→J：IfEnd: If L≥203.366：Then 494372.527→X：477948.932→Y：200→R：107°21°16.3°→W：-95→P:60.07→I: 203.366→J：IfEnd: If L≥248.491：Then 494337.373→X：478010.658→Y：86.75→R：92°24°19°→W：75→P:17.8→I: 313.333→J：IfEnd: Prog”FDCH”: Prog“XLJS”：If Z[3]≠0.0001:Then Goto 0:Else (Z[1]－X)cos(Ｏ)＋(Z[2]?－Y)sin(Ｏ)→N:L＋N→L:Int(1000L)÷1000+0.0001→L:－(Z[1]－X)sin(O)+(Z[2]－Y)cos（O）→K：If Abs（N）≥0.001：Then Goto 1：Else “L=“:L◢ “K=”:K◢ IfEnd:Goto 0: IfEnd

为了让使用者熟悉掌握数据程序文件的编写线元法的数据程序文件，现在将C、F、G等匝道的线元数据列出来，供使用者练习。

C匝道数据

半径线元名称 起点桩号 偏或缓向 和曲线元长度线元起点坐标、方位角 X(m) Y(m) 方位角 (m) 线旋转常半径数A 或缓圆曲线 直线 CK0+000.00 右偏 32 56.876 495005.991 线元478056.677 56°19′42.3″ F

线元起点坐标、方位角 缓和曲线 圆曲线 缓和曲线 线元CK0+056.876起点桩 名称 号 CK0+078.153 CK0+108.332 和曲21.277 494991. 262 478104. 123 158°09′54.4″ 长度 向 线旋X(m) Y(m) 158°09方位角 右偏 65 30.179 494971. 511 478112.036 ′54.4″ (m) 右偏 140 转常30.431 494943.127 478122.242 164°20′25.7″ 65 偏CK0+138.763 右偏 数A 30.179 494913.167 478127.219 176°47′40.6″ 直线 CK0+168.941 右偏 85 缓和曲线 GK0+000.000 右偏 175 圆曲线GK0+133.152 右偏缓和曲线 CK0+590.964 右偏 85 230 缓和曲线 GK0+282.286 右偏 190 圆曲线 CK0+675.964 右偏 85 422.023 494883.008 478126.740 182.°58′133.152 494718.802 478112.963 182°58′11.8″ 11.8″ 494587.600 478093.368 °33′17.585 149.124 494461.552 478104.874 182.199°58′″ 11.8236″ °42′21.1″ 149.72 494472.335 478002.890 72.739 494379.485 478086.681 211°37′04.2″ 匝

道数据

G匝道数据 三、道路计算程序

道路计算程序是核心程序，但是道路计算程序是固定的，使用者不需要做任何修改，也不直接运行道路计算程序分为两部分，第一部分是常规的对

半径或缓线元线元名称 起点桩号 偏和曲长度向 线旋(m) 转常数A 直线 FK0+000.00 线元起点坐标、方位角 成型线

X(m) Y(m) 方位角 型组合

100 110 120 111.525 494591.311 478420.033 252°32′00.2″ 缓和曲线 圆曲线 缓和曲线 FK0+111.505 右偏 FK0+202.434 右偏 FK0+303.548 右偏 90.909 494557.835 101.114 494542.812 130.909 494595.310 478313.650 478224.689 478142.420 252°32′00.2″ 276°12′33.6″ 328°52′35.7″ 计算的程序（文件名：XLJS），利用交点法编辑的数据文件可以直接调用；第二部分是适用于非对称的线型组合计算文件（文件名：FDCH），在数据文件调用完“FDCH”文件后，还要调用“XLJS”文件。 线路计算程序（XLJS）

Deg:Fix 4:G÷Abs(G)→V:F22≤H-T:Then L-H→J:Goto 1:IfEnd:If LLbl 3:180(J-0.5F) ÷πR→O:Q+Rsin (O)-T→M:(P+R(1-cos (O)))V→N:C+VO→O：Goto B

Lbl B:D+Mcos(C)-Nsin(C)→X:E+Msin(C)+Ncos(C)→Y :Goto C:Lbl C:If Z[3]≠0.0001:Then ?K: \"M\"?U:K<0=>180-U →U:O+U×Abs(K)÷(K+10^(-10)) →Z:X+Abs(K)×cos(Z)→X◢ Y+Abs(K)×sin(Z)→Y◢

tan-1((Y-B)÷(X-A+10^(-10))→W

X-A<0=>W+180→W:W<0=>W+360→W:√((X-A)2+(Y-B)2)→S \"DIS\": S◢ \"FW\": W>DMS◢ IfEnd

非对称线型组合计算程序（FDCH）

Abs(P+10^(-10)）÷(P+10^(-10))→S：Abs (R)÷R→V:Abs(R)→R:P2÷R→F:180(I+F)V÷(πR)→G:90FSV÷(πR)→O:If S=1:Then W-O→C:Else If S=-1:Then W-O-G→C :IfEnd:IfEnd:If C<0:Then C=C+360→C:IfEndIf C>360:Then C-360→C:IfEnd:F2÷(24R)-F^(4)÷(2688R^(3))+F^(6)÷(506880R^(5))-F^(8)

÷

→

P:F÷

2-F^(3)

÷(240R2)+F^5

÷

(34560R^(4))-F^(7)÷(8386560R^6)+F^9÷→Q:Q+(R+P)Vtan(G÷2)→T:If S=1:Then J-F+T→H:C→Z:If S=-1:Then J-I-F+T→H:C+G→Z:IfEnd:IfEnd:(F-F^(3)÷(40R2)+F^(5)÷(3456R^(4))-F^(7)÷(599040R^(6))+F^(9)÷→M:(F2÷(6R)-F^(4)÷(336R^(3))+F^(6)÷(42240R^(5))-F^(8)÷(9676800R^(7))+F^(10)÷→N:X-Mcos(Z)+Nsin(Z)→D:Y-Msin(Z)-Ncos(Z)→E 四、程序执行

程序运行分为两种，第一种是根据根据使用者输入的桩号和距中距离计算出所求点的坐标和测站点到该点的方位角和水平距离，可以直接利用计算器提供的距离和方位角进行极坐标放样。另一种是根据使用者输入的实测坐标计算出该点的桩号和距中距离，可以根据桩号和距中距离进行高边坡、隧道断面点的设计高程、以及进行路面高程控制计算方面很方便。 （一）、根据桩号和距中距离计算所求点的坐标、及方位角和水平距离。

要进行程序的线路计算，直接运行数程序文件就可以了，不需要理会计算程序。如果计算某桩号L距离中线K处的坐标X，Y和该点到测站点的距离和方位角，运行提示如下：

X0? 提示输入测站点X坐标

Y0？ 提示输入测站点Y坐标

（测站点的位数和数据程序文件里坐标取位要一致） Ln? 计算点近似桩号

近似桩号输入第四位小数后不能等于0.0001， L？ 输入计算点准确桩号

近似桩号输入第四位小数后不能等于0.0001， K？ 提示输入距离中桩的桩号

如果计算点位于线路的左幅，则在数据前加入“-”号，

如果是计算中桩坐标，则输入0

M？ 输入该点与线路的前进方向的按顺时针转过的角度，常

于斜交的涵洞，桥梁等构造物。

X 显示所求点的X坐标 Y 显示所求点的Y坐标 FW 显示所求点至测站点的方位角 DIS 显示所求点至测站点的距离 要继续计算则从L?开始，计算新点的坐标。 具体流程见程序流程图的左边部分。

例：要计算D匝道的DK0+020中桩坐标、DK0+080左边5.3米的坐标、Dk0+240

斜交75°58′11.8″、距中右侧5米处的坐标，并计算出该点至坐标为X=495005.991，Y=478056.677控制点的水平距离和方位角，运行文件“DZD”提示如 下：

X0? 495005.991 测站X坐标

Y0？ 478056.677 测站Y坐标 Ln？ 0 输入小输点第

坐标坐标测站平距

L？ 20 K? 0 M? 90 X= 494341.482 Y= 478007.903 DIS 666.2959 FW 184L？ 80 K？ -5.3 M？ 90

X= 494382.2290 Y= 477965.9594

四位起不等于0.0001的桩号

输入所求点桩号 输入距中距离 输入斜交角度 显示所求点X

显示所求点Y

显示所求点至

°11′52.42″显示所求点至

测站方位角

输入下点的桩号

DIS 630.324 FW 188°16′29.5″ L？ 240 K？ 5.0

M？ 75°58°11.8° X= 494455.3110 Y= 478101.5491 DIS 552.505 FW 175°20′29.6″

（二）、根据坐标反算桩号和距中

如果是根据某点的坐标，计算出该点在线路中的桩号和距离，程序运行时提示如下：

X0? 提示输入测站点X坐标 Y0？ 提示输入测站点Y坐标

（测站点的位数和数据程序文件里坐标取位要一致） Ln? 计算点近似桩号

近似桩号输入第四位小数后不能等于0.0001， Xp？ 所求点P的X坐标 Yp？ 所求点P的Y坐标 L 显示所求点的桩号

K 显示所求点的距中距离，如果位于线路的左边，则显示

结果为负值，如果K为0时，则该点在线路中线上。

计算新点的桩号和距中距离，则从?Xp 开始。如此循环计算。 具体流程见流程图的右边部分

例：求点X=494382.2290，Y=477965.9594在D匝道中的位置，运行文件“DZD”如下：

X0? 495005.991 测站X坐标(可

输入0)

Y0？ 478056.677 输入0)

Ln？ 60.0001 Xp? 494382.2290 坐标

Yp？ 477965.9594 坐标

L= 79.999999 号

K= -5.2989 测站Y坐标(可

输入小输点第

四位起等于0.0001的桩号

输入所求点X

输入所求点Y

显示所求点桩

显示所求点距

中距离, “-“ 表示该

点位于中桩左侧

五、程序应用

（一）本程序在高边坡放样中的应用

高边坡放样是目前公路道路测量工作普遍要做的一项工作，边坡坡顶线或者填方边坡坡脚线是边坡放样最繁琐的工作，利用本程序进行边坡放样，工作步骤如下：

1、 司镜员走到目的地后，先测出一个点的坐标和高程；

2、 利用本程序求出该点在线路中的桩号和距中距离，然后根据计算

出来的桩号和设计图纸，内插出道路边桩高程，根据道路边桩高程和道路设计宽度、坡度、设计碎落平台的宽度和间隔高度、实测点的高程计算出按实测高程，边坡距中距离。

3、 根据高程算出的距中距离和程序计算出立镜点至距中距离差值

dK，指挥司镜员靠近道路中线或者往外偏移多少米。 4、 重新做1~3工作，直到dK 小于边坡宽度控制精度要求为止。 如上图，司镜员在P点立镜，测出P点坐标和高程后，用本程序dKm 计算出P点到的中桩桩号和距离中桩距离K，根据中桩桩号和有关设计数据内插出边坡坡脚点的设计高程Hs和距离中桩的距离K1，根据m hHc、和Hs的高差计算出K2，按照P点的高程，边坡顶距中桩的理论h宽度K理论，具体计算如下： K1 K2K理论=K1+h*n1+m+h*n2+m+（Hc-Hs-2h）*n3 则dK=K- K理论

dK 求出来后，可以指挥司镜员靠近道路dK米，然后重复以上步骤，直到dK

值能满足边坡控制精度。

（二）本程序在隧道断面欠挖、超挖放样中的应用。 利用本程序进行隧道断面开挖的欠挖和超挖计算，看下图：

用全站仪测出P点的坐标和标高后，根据坐标算出P点的中桩桩号和距中距离K，根据中桩桩号查看P点的路面高程，算出P点至设计路面的高差h，根据P点所在圆弧的圆心O2在断面中的位置（K2，h2）算出P点到O2点距离R2’，用R2’-R2的差距就是P点的超挖或者欠挖的数据。

（三）本程序在路面施工控制放样中的应用

路面施工时，路面的轴线偏位、路面宽度、路面高程是测量的重要控制指标，也是道路施工控制最严格的一道测量控制。利用本程序输入桩号和距中距离，可以直接计算出所求点至测站点的水平距离和方位角，从而直接通过极坐标法直接放出该点在实地的位置。 （四）处理线路设计中的“断链”问题

断链是线路设计中经常出现的问题，断链分为长链与短链，如果是短链，那么只要弯道的计算范围从短链处分界就可以了，如果是长链，则需要把线路分成两个文件，从开始出现断链处分，这样可以避免桩号相同，但位置不同的情况。

（五）本程序在道路竣工验收过程中的应用

道路竣工验收，其主要方面就是轴线偏位和宽度，可以通过现场测量道路中线的坐标，通过坐标反算出该点和设计中线的偏移值K（轴线偏位值）。

附：

附表一：CHLNR计算出来的逐桩坐标表 C匝道坐标一览表

桩号 X(m) Y(m) 桩号 X(m) Y(m) 桩号 X(m) Y(m) CK0+000.00 495005.991 478056.677 CK0+310.00 CK0+020.00 495011.337 478075.612 CK0+330.00 CK0+035.00 495007.349 478089.930 CK0+350.00 CK0+060.00 494988.362 478105.285 CK0+370.00 CK0+145.00 494906.934 478127.439 CK0+390.00 CK0+170.00 494881.950 478126.685 CK0+410.00 CK0+190.00 CK0+210.00 CK0+230.00 CK0+250.00 CK0+270.00 CK0+290.00 494742.138 478119.431 CK0+550.00 494502.461 478106.996 494722.165 478118.395 CK0+570.00 494482.488 478105.960 494702.192 478117.358 CK0+590.00 494462.514 478104.924 494682.219 478116.322 CK0+610.00 494442.551 478103.728 494662.246 478115.286 CK0+630.00 494422.682 478101.484 494642.273 478114.250 CK0+650.00 494403.182 478097.112 494861.977 478125.648 CK0+430.00 494622.299 478113.214 CK0+670.00 494384.666 478089.632 494842.004 478124.612 CK0+450.00 494602.326 478112.177 CK0+690.00 494368.192 478078.372 494822.031 478123.576 CK0+470.00 494582.353 478111.141 CK0+710.00 494354.794 478063.585 494802.058 478122.540 CK0+490.00 494557.387 478109.846 CK0+730.00 494345.212 478046.082 494782.085 478121.503 CK0+510.00 494542.407 478109.069 CK0+748.70 494340.176 478028.110 494762.111 478120.467 CK0+530.00 494522.434 478108.032 D匝道坐标一览表

E匝道坐标一览表 坐标 坐标 坐标 桩号 桩号 桩号 X Y X Y X Y F匝道坐标一览表EK0+000.00 494360.6194 77706.293 EK0+200.00 494373.504 477903.797 EK0+020.00 494364.1494 77725.979E K0+220.00 494366.991 477922.7 G匝道坐标一览表 EK0+040.00 494367.68 477745.665E K0+240.00 494359.2964 77941.159 EK0+060.00 494371.21 477765.351E K0+260.00 494351.2724 77959.479 EK0+080.00 494374.7344 77785.038E K0+280.00 494344.0114 77978.109 六、公路纵断面高程计算 EK0+100.00 494378.0024 77804.768E K0+305.00 494338.104 478002.36 EK0+120.00 494380.4814 77824.612E K0+331.13 494338.4494 78028.394 公路纵断面计算也是采用计算程序和数据程序分开的原则，数据文件EK0+140.00 494381.6134 77844.574 EK0+160.00 494380.8924 77864.553 是专门存放道路纵断面设计参数，对于多个工程，可以建立多个数据文件，EK0+180.00 494378.18 477884.36 这样每个数据文件通过调用计算程序就可以计算出该点的高程。 （一） 数据文件

每个纵断面由纵坡和竖曲线组成，本程序每一个变坡点计算范围是从

竖曲线前的纵坡任意桩号开始到竖曲线后纵坡上的任一桩号结束。

下表是某线路路基纵断面设计的参数一览表： 根据上表编程数据文件如下：文件名BBR 从

上

面

程

序

看

，

这个语句是固定的，“Dh”提醒使用者输入所求点距离路面的厚度，比如路面底基层距离路面厚度为37cm，则输入0.37，表示所求点是路面底基层顶面高程，“L=”，输入所求点桩号,”I”=表示输入路面横坡值，“K=”表示输入路面某点至中桩距离，均为正值。以上数据均在程序运行时输入。

在每一个变坡点处，要求输入的参数有： C；输入变坡点高程 D：输入边坡点桩号

M：输入变坡点小里程段纵坡坡度 N：输入变坡点大里程侧纵坡坡度 R：竖曲线半径

L>表示竖曲线前纵坡任意里程的桩号，本程序每段竖曲线参数计算范

围是从竖曲线前的纵坡上任意点到竖曲线后纵坡上任意点的桩号。 （二） 高程计算程序（文件名：GCJS）

高程计算文件代码如下：

（三） 运行

运行BBR程序显示如下

Dh？ 提示输入至路面厚度 L=? 提示输入桩号 i？ 提示输入横坡度

K？ 提示输入距中距离，均为正数 H= 提示输出所求点高程

例：根据上表计算K38+660中桩高程、K38+900右侧4.5米高程、K39+600底基层面左侧4.5米高程（假设路面沥青厚度为6cm，封釉层1cm，水泥稳定层基层厚度为30cm） 运行BBR文件

Dh？ 0 L=？ 38660 i？ 0 K？ 0 H= 76.9412 L=？ 38900 i？ 0.02 K？ 4.50 H= 85.2200 退出程序，重新运行

Dh？ 0.370 L=？ 39600 i？ 0.02 K？ 4.50 H= 93.0509

附：