5800道路程序终极版附件为隧道版.docx

1、5800道路程序终极版附件为隧道版附件(点击下载)：功能:坐标正反算(含高程),把要素内置化(无需改程序文件,即可更换路线,同时存两条线要素)另含两个边仰(一个横向坡,一个纵向坡)坡放样模块.帮助文件未完成程序浏览:FileName:RESET初始化程序Norm 1:50C:12345JRESET PW?I:I=J=500DimZ 为数据库增加额外变量500个，在SET、SETPFDYS程序根据实际再增减变量50ZC+22PASSWORDS?I:IZC+39重设要素保护密码Cls:StopFileName:DATLOCK要素保护密码确定认程序Cls:Norm 1:50C PASSWORDS?J

2、:Cls:JZC+39=Stop2008-08-07修改FileName:SHELL(外壳程序)50C ;在扩充变量预留前50个给别的程序用。如不够就适量加大。RESET，SHELL，SET，SETPFDYS，这几个程序中C值必需一致ZC+35I:STATION-N?I:IZC+35设置测站N坐标ZC+36I:STATION-E?I:IZC+36设置测站E坐标ZC+44I:STATION-Z?I:IZC+44设置测站Z坐标ZC+45I:STATION-HI?I:IZC+45 设置仪高ZC+41I:GC-DH?I:IZC+41 放样点高差常数Norm 1ZC+23N当前分段要素N坐标ZC+24E

3、当前分段要素E坐标ZC+25M当前分段要素起点桩号ZC+26H当前分段要素起点方位角（正北）单位：弧度ZC+27A当前分段要素起点曲率有左偏负右偏正（注意不半径）ZC+28R当前分段要素终点曲率有左偏负右偏正（注意不半径）ZC+29L当前分段要素长度ZC+3WZC+32I:DAT1 2 3?I:IZC+32平曲线要素数库类型选择1为内置式，2文件式，3实时输入(查看当前要素值) 0IZC+4I:1 2 3?I:IZC+4 选择本程序模式默认为坐标正算，1坐标反算，2横向边仰坡放样，3，纵向边仰坡放样(隧道进出口用到)Lbi 1ZC+4II=1=ProgTURNZHI=2=ProgHXBYPI=

4、3=ProgZXBYPLbi ADeg:Norm 1:ClsZC+1G:ZH?G:GZC+1 ;G桩求桩号ZC+2B:JL?B:BZC+2 ;B横向距离，左正右负Lbi BProg 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+ZC+41；第四行显示Z坐标Locate 10,4,ZC+2 ；第四行显示横向距离0IDo:I+1I:I=1000=Goto

5、1:LpWhile Not(Getkey=57 Or Getkey=27) ;锁定键盘，并在几十秒后自动返回Getkey=57= Goto 1Z+ZC+41+ZC+49-ZC+44-ZC+45ZC+48Pol(ZC+35-X,ZC+36-Y)Cls:DL:Locate 3,1,I:Locate 12,1,ZC+45DH:Locate 3,2,ZC+48:Locate 12,2,ZC+49FWJJJ+360J:JDMS 2008-08-07 修改JZC+46IZC+47Goto 1FileName：ZBJS坐标计算程序ProgREADDATRadG-MQIF AR=0 :Then If A=R

6、:Then 1J:Else 3J :IfEnd:Else If A=R :Then 2J:Else 3J:IfEnd:IfEndIf J=1 :Then HF:N+QCos(H)X:E+QSin(H)Y:IfEnd 直线段直接计算If J=2 :Then H+QRF:Rec(Abs(2Sin(Abs(0.5QR)R),H+QR2):N+IX:E+JY:IfEnd圆弧段直接计算If J=3 :Then Goto 5 :IfEndX+BCos(F+W)XY+BSin(F+W)YDeg:ReturnLbi 5 用五点通用坐标计算计算缓和段0.5(R-A)LKAQIKQ2J0.0469100770P:

7、H+IP+JP2U 0.2307653449P:H+IP+JP2V0.5P :H+IP+JP2D 0.7692346551P:H+IP+JP2F 0.9530899230P:H+IP+JP2T 0.1184634425I0.2393143352J0.2844444444ON+Q(ICos(U)+JCos(V)+OCos(D)+JCos(F)+ICos(T)XE+Q(ISin(U)+JSin(V)+OSin(D)+JSin(F)+ISin(T)YH+AQ+KQ2FX+BCos(F+W)XY+BSin(F+W)YDeg:ReturnFileName:TURNZH(坐标反算) Norm 1ZC+37

8、I:DQD-N?I:IZC+37 输入待求点N坐标I=-1=ProgINFWJDLDHZC+38I:DQD-E?I:IZC+38 输入待求点E坐标0B：M+L2G：ProgZBJSLbi S:RadZC+37-XI:ZC+38-YJIf I=0 And J=0 :Then Goto A:IfEndPol(I,J)Lbi A:Rec(I,J-F):JZC+2G+IG:IF Abs(I)0.0001 :Then ProgZBJS:Goto SGZC+1FileName：INFWJDLDH实测坐标（用方位角，距离，高差）输入程序Deg:Norm 1:ClsZC+46J:DQD-FWJ?J:JZC+4

9、6ZC+47I:DQD-DL?I:IZC+47ZC+48K:DQD-DZ?K:KZC+48ZC+49P:RHT?P:PZC+49Rec(I,J)ZC+35+IZC+37ZC+36+JZC+38ZC+44+ZC+45+K-PZC+43FileName:HXBYP 横向边仰坡放样程序Cls:Norm 1ZC+40I:QPD-DL?I:IZC+40输入起坡点与中桩距离常数ZC+41I:QPD-DZ?I:IZC+41 输入起坡点与中桩高差常数ZC+42I:i?I:IZC+42输入边仰坡坡度，左仰坡(路堑)为正，右仰坡(路堑)为负，左边坡为负，右边坡为正Cls:ZC+43I:DQD-Z?I:IZC+43

10、 输入实测高程ProgTURNZHProg GCJS(ZC+43-Z-ZC+41)ZC+42+ZC+40BFix 3:B-ZC+2 显示与上一次测点的偏差BZC+2FileName:ZXBYP纵向仰坡放样程序(隧道进口使用)Cls:Norm 1:ZC+33I:QPD-ZH?I:IZC+33输入起坡点桩号ZC+34I:QPD-Z?I:IZC+34 输入起坡点高程ZC+42I:i?I:IZC+42仰坡时(隧道进出口仰坡)，进口为正，出口为负)Cls:ZC+43I:DQD-Z?I:IZC+43 输入实测高程ProgTURNZHZC+42(ZC+43-ZC+34)+ZC+33GFix 3:G-ZC+1

11、 显示与上一次测点的偏差GZC+1FileName:SET 设置程序1Lbi S2Norm 1350C450ZC+2250K6ZC+5I:1 Or 2?I:IZC+5;选择1线或者2线7180ZC+3J:XZJJ?J:J180ZC+3;斜桩夹角8If I=1:Then ZC+8ZC+6:ZC+9ZC+79ZC+12ZC+16:ZC+13ZC+17:ZC+18ZC+2010Else ZC+10ZC+6:ZC+11ZC+711ZC+14ZC+16:ZC+15ZC+17:ZC+19ZC+2012IfEnd13ZC+7+3ZC+31140I:1PQX,2SQX,3YSDS?I ;1输入平曲线要素，2输

12、入竖曲线要素，3输入1、2线夹平竖曲线交点(变坡点)数量以让程序规划内存建立要素数据为，其它值退出本程序15I0=ProgDATLOCK 2008-07-28 18：23修改16I=1 = Goto 1:I=2 = Goto 2:I=3 = Goto 32008-07-28 18：23修改17ProgSETPFDYS18Stop19Lbi 1 ;平曲线要素输入部分20ZC+6D21-1ZD+322For 1J To ZC+1623Cls:PQX:Locate 9,1,J24ZD+1I:N?I:IZD+1;请输入N坐标25ZD+2I:E?I:IZD+2 ;请输入E坐标26ZD+3I:R?I:IZ

13、D+3;请输入半径R27ZD+4I:LS1?I:IZD+4 ;请输入缓和曲线1长度28ZD+5I:LS2?I:IZD+5;请输入缓和曲线2长度29D+5D30Next311ZD-232Goto S33Lbi 2 ;竖曲线要素输入部分34ZC+7D35-1ZD+336For 1J To ZC+1737Cls:SQX:Locate 9,1,J38ZD+1I:BPDZH?I:IZD+1 ;请输入变坡点桩号2008-08-07修改39ZD+2I:BPDGC?I:IZD+2 ;请输入变坡点高程 2008-08-07修改40ZD+3I:R?I:IZD+3 ;请输入半径R 2008-08-07修改41D+3

14、D42Next431ZD44Goto S45Lbi 3 ;要素点数输入部分46ZC+12I:1PQXDS?I:IZC+12;输入1线平曲线点数,最小值为347ZC+13I:1SQXDS?I:IZC+13;输入1线竖曲线点数,最小值为348ZC+18I:1JDZH?I:IZC+18 ;输入1线平曲线第一个交点的交点桩号49ZC+14I:2PQXDS?I:IZC+14;输入2线平曲线点数,最小值为350ZC+15I:2SQXDS?I:IZC+15;输入2线竖曲线点数,最小值为351ZC+19I:2JDZH?I:IZC+19;输入2线平曲线第一个交点的交点桩号52C+ZC+22ZC+8 1线路平曲线

15、要素指针偏移基数53ZC+8+5ZC+12ZC+91线路竖曲线要素指针偏移基数54ZC+9+3ZC+13ZC+10 2线路平曲线要素指针偏移基数55ZC+10+5ZC+14ZC+112线路竖曲线要素指针偏移基数56ZC+11+3ZC+15ZC+21 要素数据库结束指针57ZC+21+1ZC+3058ZC+21DimZ 增加额外变量59-1K60Goto SFileName:SETPFDYS 交点要素=分段要素1Norm 1250C3Rad4ZC+6+5D5ZC+21+1F6F+6DimZ7-1ZF设置分段要素开始标志8ZC+20V ; 从数据库读入当前线第一个平曲线交点桩号90G10Pol(Z

16、D+1-ZD-4,ZD+2-ZD-3)；计算第一条直线的距离I，计算方位角J11JH12ZD-4ZF+1直线段13ZD-3ZF+214V-IZF+315For 1G To ZC+16-216Pol(ZD+1-ZD-4,ZD+2-ZD-3) ；计算第本交点与下一交点的距离I，计算方位角J17JH 2008-08-07 修改18Pol(ZD+6-ZD+1,ZD+7-ZD+2)19IW:J-HK ；计算转角K，负值是左转角，正值是右转角 2008-08-07 修改20IF Abs(K):Then If K0:Then K-2K:Else K+2K:IfEnd:IfEnd21ZD+3R;从数据库读入半

17、径R22ZD+4L;从数据库读入缓和曲线1长度23ZD+5M;从数据库读入缓和曲线2长度24L2RB;计算缓和曲线1的B0角25M2RO;计算缓和曲线2的B0角26(Abs(K)-B-O)RP;计算圆曲线长度27L224R-L(4) 2384R(3) Q;内移值计算方法228L2-L(3) 240R2S ;切线增长值计算方法229M224R-M(4) 2384R(3) J;内移值计算方法230M2-M(3) 240R2I ;切线增长值计算方法231Abs(2- Abs(K) A32Rtan(Abs(K2)+S+JCos(A)-QTan(A)T;计算切线133Rtan(Abs(K2)+I+QCo

18、s(A)-JTan(A)U;计算切线234HH+2H351A36K-1A37HZF+4380ZF+5390ZF+640F+6F41F+6DimZ42Rec(T,H)缓和段143ZD+1-IN44ZD+2-JE45NZF+146EZF+247V-TZF+348L=0=Goto A49HZF+4500ZF+551ARZF+652F+6F53F+6DimZ54Lbi A圆弧段550I:0J56L=0=Goto B57Pol(L-L(3)40R2,L26R-L(4)336R(4)58Lbi B59Rec(I,H+AJ)60N+IN61E+JE62NZF+163EZF+264V-T+LZF+365H+A

19、BZF+466ARZF+567ARZF+668F+6F69F+6DimZ70M=0=Goto C缓和段271Rec(2RSin(0.5PR),H+AB+APR2)72N+IN73E+JE74NZF+175EZF+276V-T+L+PZF+377H+AB+APRZF+478ARZF+5790ZF+6 此处删除2008-08-06 修改80F+6F81F+6DimZ82Lbi C此处插入2008-08-06 修改83Rec(U,H+K)84ZD+1+IZF+1此处2008-08-07 修改85ZD+2+JZF+2此处2008-08-07 修改86V-T+L+P+MZF+387V-T+L+P+M-U

20、+WV此处2008-08-07 修改88D+5D89Next90H+KZF+4910ZF+5920ZF+693F+9DimZ941ZF+795VZF+9 要素终点桩号96DegFileName：READDATZC+32II=3=ProgREADDAT3 If GM+L:Then ;验证要素是否可用，否则重设要素I=1=ProgREADDAT1 I=2=ProgREADDAT2 ProgREALDAT3IfEndFileName:READDAT1 平曲线分段要素读取程序1ZC+30F0JLbi 1 搜索要素If GGoto E:F-6F:-1J:Goto 1Else 1=ZF+1=Goto E:

21、-1=J=Goto 2:1J:F+6F:Goto 1IfEndLbi 2FZC+30ZF+1ZC+23当前分段要素N坐标ZF+2ZC+24当前分段要素E坐标ZF+3ZC+25当前分段要素起点桩号ZF+4C+26当前分段要素起点方位角（正北）ZF+5ZC+27当前分段要素起点曲率有左偏负右偏正（注意不半径）ZF+6ZC+28当前分段要素终点曲率有左偏负右偏正（注意不半径）ZF+9-ZF+3ZC+29当前分段要素长度ReturnLbi ECls ERRORStopREADDAT2文件请输入内容为以下内容以便测试结果与本文一致。180JGoto 1Lbi 1If G30008.396:Then Cl

22、s:Stop:IfEnd ;判断待求桩号是否小于本要素起点，如小于就停止运行程序If G30661.038(终点里程):Then 2851422.090(起点X)N:513896.440(起点Y)E:30008.396M:2133252J(起点方位角)H:0(起点曲率)A:0(终点曲率)R:652.642起点至终点长度L:Goto E:IfEndIf G30921.038:Then 2850878.162N:513535.768E:30661.038M:2133252JH:0A:11300R:260L:Goto E:IfEndIf G31426.367:Then 2850666.478N:51

23、3385.004E:30921.038M:2191639JH:11300A:11300R:505.328L:Goto E:IfEndIf G31686.367:Then Cls:Stop:IfEnd ;判断待求桩号是否大于本要素终点，如大于就停止运行程序Lbi ENZC+23:EZC+24:MZC+25:HZC+26:AZC+27:RZC+28:LZC+29FileName:READDAT3 平曲线分段要素读取程序3Norm 1ZC+23N:?N:NZC+23 当前分段要素N坐标ZC+24E:?E:EZC+24 当前分段要素E坐标ZC+25M:SZH?M:MZC+25 当前分段要素起点桩号180ZC+26H:HDMSFWJ?H:H180H:HZC+26 当前分段要素起点方位角（正北）单位：度分秒ZC+27A:SR?A:AZC+27 当前分段要素起点曲率有左偏负右偏正（注意是半径的倒数）ZC+28R:ER?R:RZC+28当前分段要素终点曲率有左偏负右偏正（注意是半径倒数）ZC+29L:?L:LZC+29当前分段要素长度ReturnFileName:GCJS