搜索
5800公路测量程序使用说明一、程序使用流程本程序数据和主程序是分开的，编程时将不同的工程数据存放到不同的数据文件里，如A匝道，文件名为A，将匝道A所有的曲线线元参数输入A文件里。运行时只要运行文件名A的程序就可以了，具体运行流程见下图：二、数据文件的编写（一）交点法数据文件编辑交点法编写数据文件必须是对称型的，即直线段→缓和曲线段→圆曲线段→缓和曲线段→直线段，（如果任意一端没有直线段，则把直线段长度看做是0），另外圆曲线两侧缓和曲线的旋转常数必须相等，并且和直线段连接处的半径必须是无穷大。交点法数据文件编写一般是根据设计图纸提供的平面曲线参数一览表提供的参数来编写，每个弯道包括：弯道起点方位角（C），交点X坐标（D），交点Y坐标（E），缓和曲线长度（F，当没有设缓和曲线时，F=0），交点转交（G，向左转弯，G为负值，向右转弯，G取正值），交点桩号（H），弯道圆曲线半径（R）。下图是一段市政道路设计参数数据。根据上图提供的数据，可以编辑成如下的数据文件：文件名:CHLNR3→DimZ“X0”?A:”Y0”?B:“Ln”?L:Abs(L)-Int(Abs(1000L))/1000→Z[3]:Lbl0:IfZ[3]≠0.0001:Then?L:Else“Xp”?X:”Yp”?Y:X→Z[1]：Y→Z[2]：IfEnd:Lbl1:IfL0:Then98°39°35.12°→C：4474.384→D：2415.861→E：140→F：31°17°23°→G：410.007→H:600→R:IfEnd:IfL1060Then129°56°58.19°→C:4206.421→D：3093.946→E：70→F：-33°50°48°→G：1285.437→H:600→R:IfEnd:Prog”XLJS”:IfZ[3]≠0.0001:ThenGoto0:Else(Z[1]－X)cos(Ｏ)＋(Z[2]－Y)sin(Ｏ)→N:L＋N→L:－(Z[1]－X)sin(O)+(Z[2]－Y)cos（O）→K：IfAbs（N）≥0.001：ThenGoto1：Else“L=“:L◢“K=”:K◢IfEnd:Goto0:IfEnd在面程式中，有两个条件转移语句即IfL0：Then98°39°35.12°→C:4774.384→D:2415.861→E:140→F:31°17°23°→G:410.007→H:600→R:IfEndIfL1060：Then129°56°58.19°→C:4206.421→D:3093.946→E:70→F:-33°50°48°→G:1285.437→H:600→R:IfEnd……如果还有其他弯道，可以继续完后加。在这些存放设计参数的语句前后的程序表达式是固定的。说明：编辑曲线参数时，每个曲线参数放在一个IfL***.***（两个弯道中间直线段上的任意桩号）:Then***°**°**.**°→C（弯道起点方位角）：****.***→D（交点X坐标）：****.****→E（交点Y坐标）：***.***→F（缓和曲线长度：**°**°**°→G（转角，向左转为负值，向右转为正值）：***.***→H（交点里程桩号）:***→R（圆曲线半径）:IfEnd条件式语句里，如果有多个弯道，一直按上述形式编下去，变量说明：1、L***.***弯道参数的起点范围，可以是弯道起点，也可以是弯道前直线段里的任意点桩号2、***°**°**.**°→C弯道起点方位角3、****.***→D、****.****→E交点X坐标、Y坐标4、***.***→F缓和曲线长度，如果没设置缓和曲线，则缓和曲线长度看做05、**°**°**°→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匝道由三个线元组成，第一段是缓和曲线段，根据上表判断，该缓和曲线属于第一段缓和曲线，其最小半径处桩号为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.926Y=477989.873最小半径R=40最小半径处方位角：W=302°18′26.0″缓和曲线旋转常数P=55与之相连的圆曲线长度I=39.216最小半径处桩号J=39.289此段数据计算范围可以从K0+000~K0+78.504第二段缓和曲线线元参数线元名称起点桩号偏向半径或缓和曲线旋转常数A线元长度(m)线元起点坐标、方位角X(m)Y(m)方位角缓和曲线DK0+000.00右偏5539.289494341.902478027.825260°38′55.7″圆曲线DK0+039.289右偏4039.216494347.926477989.873302°18′26.0″缓和曲线DK0+078.504右偏65105.625494380.677477971.271358°28′45.9″最小半径处坐标X=494380.677Y=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时，曲线参数按第一段编写，当L78.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]:Lbl0:IfZ[3]≠0.0001:Then?L:Else“Xp”?X:”Yp”?Y:X→Z[1]：Y→Z[2]：IfEnd:Lbl1:IfL≥0:Then494347.926→X：477989.873→Y：40→R：302°18°26°→W：55→P:39.216→I:39.289→J：IfEnd:IfL≥78.504：Then494380.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”:IfZ[3]≠0.0001:ThenGoto0:Else(Z[1]－X)cos(Ｏ)＋(Z[2]－Y)sin(Ｏ)→N:L＋N→L:－(Z[1]－X)sin(O)+(Z[2]－Y)cos（O）→K：IfAbs（N）≥0.001：ThenGoto1：Else“L=“:L◢“K=”:K◢IfEnd:Goto0:IfEnd在上面程式中，每个缓和曲线的线元参数放在一个If…Then……IfEnd循环语句中，If后面紧跟计算的起点范围，终点范围从下段的起点范围结束。变量定义：1、X、Y输入缓和曲线半径最小处的X、Y坐标2、R缓和曲线最小半径，当线路左转时，半径数据前面加“-”号，以判断线路转向。3、W输入缓和曲线半径最小处的方位角4、P缓和曲线旋转常数，第一缓和曲线长度为正数，第二缓和曲线常数为负数。当不设缓和曲线时，P=05、I与缓和曲线最小半径端连接的圆曲线长度。6、J缓和曲线小半径段的桩号。再看看E匝道数据文件编辑E匝道线元设计数据表在E匝道线元表中，我们可以看出，E匝道式由一段直线段→缓和曲线段→圆曲线段→缓和曲线段→缓和曲线段→圆曲线段连接起来的。值得注意的是，在位于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]:Lbl0:IfZ[3]≠0.0001:Then?L:Else“Xp”?X:”Yp”?Y:X→Z[1]：Y→Z[2]：IfEnd:Lbl1:IfL≥0:Then494381.6316→X：477847.870→Y：200→R：90°08°44.6°→W：120→P:60.07→I:143.296→J：IfEnd:IfL≥203.366：Then494372.527→X：477948.932→Y：200→R：107°21°16.3°→W：-95→P:60.07→I:203.366→J：IfEnd:IfL≥248.491：Then494337.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”：IfZ[3]≠0.0001:ThenGoto0: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：IfAbs（N）≥0.001：ThenGoto1：Else“L=“:L◢线元名称起点桩号偏向半径或缓和曲线旋转常数A线元长度(m)线元起点坐标、方位角X(m)Y(m)方位角直线EK0+000.0071.296494360.619477706.29379°49′57.0″缓和曲线EK0+071.296右偏12072494373.205477776.47079°49′57.0″圆曲线EK0+143.296右偏20060.070494381.631477847.87090°08′44.6″缓和曲线EK0+203.366右偏9545.125494372.527477907.018107°21′16.3″缓和曲线EK0+248.491左偏7564.841494355.878477948.932113°49′05.6″圆曲线EK0+313.333左偏86.7517.800494337.373478010.65892