5800隧道专版(重点)

1、5800隧道专版本版本是道路版的升级版，程序只改变了SHELL程序，从SHELL中拆分出WORK-SET(工作设置程序)，加入SDPY(隧道放样)。 FileName:RESET    初始化程序 Norm 1:50C:12345J "RESET PW"?I:I=J=>500DimZ   为数据库增加额外变量500个，在SET、SETPFDYS程序根据实际再增减变量以保证不浪费内存 50ZC+22：本程序设置变量个数目前50个刚刚够用。 "PASSWORDS"?I:IZC+39 

2、     重设要素保护密码 Cls:Stop FileName:DATLOCK      要素保护密码确定认程序 Cls:Norm 1:50C "PASSWORDS"?J:Cls:JZC+39=>Stop FileName:SHELL(外壳程序) 50C    在扩充变量预留前50个给别的程序用。如不够就适量加大。RESET，SHELL，SET，SETPFDYS，这几个程序中C值必需一致 Norm 1 ZC+23N当前分

3、段要素N坐标 ZC+24E当前分段要素E坐标 ZC+25M当前分段要素起点桩号 ZC+26H当前分段要素起点方位角（正北）单位：弧度 ZC+27A当前分段要素起点曲率有左偏负右偏正（注意不是半径） ZC+28R当前分段要素终点曲率有左偏负右偏正（注意不是半径） ZC+29L当前分段要素长度 ZC+3W   斜桩角度 Lbi S Prog"PROGMODE" ;进入模式功能选择 Lbi A Deg:Norm 1:Cls ZC+1G:"ZH"?G    桩求桩号 If G=-1:Then Prog&

4、quot;WORK-SET":Goto S:IfEnd    进入工作设置 GZC+1 ZC+2B:"JL"?B:BZC+2 ;B横向距离，左正右负 Lbi B Prog "ZBJS"     ;坐标正算 Fix 3:Cls ;设置三位小数 "ZH":Locate 4,1,G ;第一行显示桩号 "X":Locate 3,2,X ;第二行显示X(N)坐标 "Y":Locate 3,3,Y；第三行显示Y(E)坐

5、标 Prog "GCJS"     ;高程计算 "Z":Locate 2,4,Z+ZC+41 ；第四行显示Z坐标 Locate 10,4,B         ；第四行显示横向距离 0I Do:I+1I:I=1000=>Goto 1:LpWhile Not(Getkey=57 Or Getkey=27 Or Getkey=26) ;锁定键盘，并在几十秒后自动返回 Getkey=57=> Goto 1 If Getkey=26=

6、> Prog "SDFY":Goto B:EndIf Pol(ZC+35-X,ZC+36-Y) Cls:"DL":Locate 3,1,I:Locate 12,1,ZC+45 "":Locate 3,2,B:Locate 12,2,ZC+49 "FWJ" J<0=>J+360J:JDMS JZC+46 IZC+47 Goto 1 FileName：SDFY 隧道放样 Cls ZC+2V ZC+41D ZC+43I:"DQD-Z"?I:IZC+43 ZC+42J:"R&

7、quot;?J:JZC+42 ZC+40K:"R-DL"?K:KZC+40 Abs(V-K)P (P2+(I-Z-D)2)U     计算实际半径 Fix 3:Cls ;设置三位小数 "ZH":Locate 4,1,G "R":Locate 2,2,U Locate 10,2,U-J If I-Z-D-J>0 :Then I-Z-D-(J2-P2)P:0O Else If P-J>0:Then P-(J2-(I-Z-D)2)O:0P Else P-(J2-(I-Z-D)2)O:I-Z

8、-D-(J2-P2)P IfEnd:IfEnd "H":Locate 2,3,O Locate 8,3,"V" Locate 9,3,P Locate 1,4,Z+D Locate 9,4,I 0T Do:T+1T:T=1000=>Return:LpWhile Not(Getkey=57 Or Getkey=27) ;锁定键盘，并在几十秒后自动返回 Getkey=57=>Return If V-K>0 :Then V-OB:Else V+OB:IfEnd FileName：ZBJS  坐标计算程序 Prog"

9、;READDAT" Rad G-MQ IF AR=0 :Then If A=R :Then 1J:Else 3J :IfEnd:Else If A=R :Then 2J:Else 3J:IfEnd:IfEnd 判断线元类型 If 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+QR÷2):N+IX:E+JY:IfEnd  圆弧段直接计算 If J=3 :Then Goto 5

10、:IfEnd X+BCos(F+W)X Y+BSin(F+W)Y Deg:Return Lbi 5     用五点通用坐标计算计算缓和段 0.5(R-A)÷LK AQI KQ2J 0.0469100770P:H+IP+JP2O 0.2307653449P:H+IP+JP2T 0.5P         :H+IP+JP2D 0.7692346551P:H+IP+JP2F 0.9530899230P:H+IP+JP2P 0.1184634425I 0.239314335

11、2J 0.2844444444Y N+Q(ICos(O)+JCos(T)+YCos(D)+JCos(F)+ICos(P)X E+Q(ISin(O)+JSin(T)+YSin(D)+JSin(F)+ISin(P)Y H+AQ+KQ2F X+BCos(F+W)X Y+BSin(F+W)Y Deg:Return FileName:TURNZH(坐标反算) Prog"INNEZ" ZC+37U ZC+38V 0B：M+L÷2G：Prog"ZBJS" Lbi S:Rad U-XI:V-YJ If I=0 And J=0 :Then Goto A:IfEn

12、d Pol(I,J) Lbi A:Rec(I,J-F) G+IG:IF Abs(I)>0.0001 :Then Prog"ZBJS":Goto S JB GZC+1 JZC+2 Prog"GCJS" FileName：INNEZ 实测坐标输入程序 Lbl S:Norm 1:Cls ZC+50=2=>Goto 1 ZC+37I I<0=>Goto 1 "DQD-N"?I     输入待求点N坐标 I<0=>Goto 1 IZC+37 ZC+38I:"D

13、QD-E"?I    输入待求点E坐标 I<0=>Goto 1 IZC+38 ZC+43I:"DQD-Z"?I 输入实测高程 I<0=>Goto 1 IZC+43     Return Lbi 1 I=-1=>Then 2ZC+50 If I=-2:Then Prog"WORK-SET":Goto S:IfEnd    进入工作设置 Prog"INFWJDLDH" ZC+50=

14、1=>Goto S Return FileName：INFWJDLDH实测坐标（用方位角，距离，高差）输入程序 Lbi S:Deg:Norm 1:Cls ZC+46J:"DQD-FWJ"?J J<0=>Goto 1 JZC+46 ZC+47I:"DQD-DL"?I:IZC+47 ZC+48K:"DQD-DZ"?K:KZC+48 ZC+49P:"RHT"?P:PZC+49 Rec(I,J) ZC+35+IZC+37 ZC+36+JZC+38 ZC+44+ZC+45+K-PZC+43 Return L

15、bi 1 If J=-2:Then Prog"WORK-SET":Goto S:IfEnd      进入工作设置 1ZC+50 Return FileName: WORK-SET 测站设置程序 Norm 1:Cls 50C ZC+32I:"DAT1 2 3"?I:IZC+32  平曲线要素数库类型选择1为内置式，2文件式，3实时输入(查看当前要素值) ZC+4I:"0 1 2 3"?I:IZC+4   选择本程序模式0为坐标正算，1坐

16、标反算，2横向边仰坡放样，3，纵向边仰坡放样(隧道进出口用到) ZC+41I:"GC-DH"?I:IZC+41     放样点高差常数 ZC+35I:"STATION-N"?I:IZC+35      设置测站N坐标 ZC+36I:"STATION-E"?I:IZC+36      设置测站E坐标 ZC+44I:"STATION-Z"?I:IZC+44设置测

17、站Z坐标 ZC+45I:"STATION-HI"?I:IZC+45     设置仪高 Cls FileName: PROGMODE 程序功能模式选择 ZC+4I I=0=>Return ;正算模式 I=1=>Prog"TURNZH"      反算模式 I=2=>Prog"HXBYP"     ;横向边仰坡放样模式 I=3=>Prog"ZXBYP" ;纵

18、向边仰坡放样模式 FileName:HXBYP   横向边仰坡放样程序 Cls:Norm 1 ZC+40I:"QPD-DL"?I:IZC+40    输入起坡点与中桩距离常数 ZC+41I:"QPD-DZ"?I:IZC+41       输入起坡点与中桩高差常数 ZC+42I:"i"?I:IZC+42          

19、;  输入边仰坡坡度，左仰坡(路堑)为正，右仰坡(路堑)为负，左边坡为负，右边坡为正 Cls:ZC+43I:"DQD-Z"?I:IZC+43     输入实测高程 Prog"TURNZH" Prog "GCJS" (ZC+43-Z-ZC+41)ZC+42+ZC+40B Fix 3:B-ZC+2 显示与设计位置的偏差 BZC+2 FileName:ZXBYP  纵向仰坡放样程序(隧道进口使用) Cls:Norm 1:ZC+33I:"QPD-ZH&q

20、uot;?I:IZC+33    输入起坡点桩号 ZC+34I:"QPD-Z"?I:IZC+34       输入起坡点高程 ZC+42I:"i"?I:IZC+42            仰坡时(隧道进出口仰坡，进口为正，出口为负) Cls:ZC+43I:"DQD-Z"?I:IZC+43  

21、0;  输入实测高程 Prog"TURNZH" ZC+42(ZC+43-ZC+34)+ZC+33G Fix 3:G-ZC+1 显示与设计位置的偏差 GZC+1 FileName:SET 设置程序 Lbi S Norm 1 50C 50ZC+22 本程序设置变量个数目前50个刚刚够用。 0K ZC+5I:"1 Or 2"?I:IZC+5 ;选择1线或者2线 180ZC+3÷J:"XZJJ"?J:J÷180ZC+3;斜桩夹角 If I=1:Then ZC+8ZC+6:ZC+9ZC+7 ZC+12ZC+16:Z

22、C+13ZC+17:ZC+18ZC+20 Else ZC+10ZC+6:ZC+11ZC+7 ZC+14ZC+16:ZC+15ZC+17:ZC+19ZC+20 IfEnd   ZC+7+3ZC+31 0I:"1PQX,2SQX,3YSDS"?I   1输入平曲线要素，2输入竖曲线要素，3输入1、2线夹平竖曲线交点(变坡点)数量以让程序规划内存建立要素数据为，其它值退出本程序 I0=>Prog"DATLOCK" I=1 => Goto 1:I=2 => Goto 2:I=3 => Goto 3

23、Prog"SETPFDYS" Stop Lbi 1 ;平曲线要素输入部分 ZC+6D -1ZD+3 For 1J To ZC+16 Cls:"PQX":Locate 9,1,J ZD+1I:"N"?I:IZD+1     请输入N坐标 ZD+2I:"E"?I:IZD+2 ;请输入E坐标 ZD+3I:"R"?I:IZD+3      ;请输入半径R ZD+4I:"LS1"?I:IZD+4 ;请输入缓

24、和曲线1长度 ZD+5I:"LS2"?I:IZD+5        请输入缓和曲线2长度 D+5D Next 1ZD-2 Goto S Lbi 2                             

25、60; ;竖曲线要素输入部分 ZC+7D -1ZD+3 For 1J To ZC+17 Cls:"SQX":Locate 9,1,J ZD+1I:"BPDZH"?I:IZD+1     请输入变坡点桩号 ZD+2I:"BPDGC"?I:IZD+2 ;请输入变坡点高程 ZD+3I:"R"?I:IZD+3         ;请输入半径R D+3D Next 1ZD Goto S Lbi 3 

26、                            ;要素点数输入部分 ZC+12I:"1PQXDS"?I:IZC+12 ;输入1线平曲线点数,最小值为3 ZC+13I:"1SQXDS"?I:IZC+13 ;输入1线竖曲线点数,最小值为3 ZC+18I:"1JDZH"

27、;?I:IZC+18     输入1线平曲线第一个交点的交点桩号 ZC+14I:"2PQXDS"?I:IZC+14 ;输入2线平曲线点数,最小值为3 ZC+15I:"2SQXDS"?I:IZC+15 ;输入2线竖曲线点数,最小值为3 ZC+19I:"2JDZH"?I:IZC+19 ;输入2线平曲线第一个交点的交点桩号 C+ZC+22ZC+8             &

28、#160; 1线路平曲线要素指针偏移基数 ZC+8+5ZC+12ZC+9              1线路竖曲线要素指针偏移基数 ZC+9+3ZC+13ZC+10         2线路平曲线要素指针偏移基数         ZC+10+5ZC+14ZC+11   &#

29、160;  2线路竖曲线要素指针偏移基数 ZC+11+3ZC+15ZC+21     要素数据库结束指针 ZC+21+1ZC+30 ZC+21DimZ                 增加额外变量 -1K Goto S FileName:SETPFDYS   交点要素=>分段要素 Norm 1 50C Rad ZC+6+5D ZC+21+1F F+6D

30、imZ -1ZF          设置分段要素开始标志 ZC+20V         ; 从数据库读入当前线第一个平曲线交点桩号 0G Pol(ZD+1-ZD-4,ZD+2-ZD-3) ；计算第一条直线的距离I，计算方位角J JH ZD-4ZF+1          直线段 ZD-3ZF+2 V-IZF+3 For

31、1G To ZC+16-2   Pol(ZD+1-ZD-4,ZD+2-ZD-3)    ；计算第本交点与下一交点的距离I，计算方位角J JH Pol(ZD+6-ZD+1,ZD+7-ZD+2) IW:J-HK                             ；计算转角K，负

32、值是左转角，正值是右转角 IF Abs(K)>:Then If K>0:Then K-2K:Else K+2K:IfEnd:IfEnd ZD+3R ;从数据库读入半径R ZD+4L ;从数据库读入缓和曲线1长度 ZD+5M ;从数据库读入缓和曲线2长度 L÷2÷RB ;计算缓和曲线1的B0角 M÷2÷RO ;计算缓和曲线2的B0角 (Abs(K)-B-O)RP ;计算圆曲线长度 L2÷24÷R-L(4) ÷2384÷R(3) Q      内移

33、值计算方法2 L÷2-L(3) ÷240÷R2S     ;切线增长值计算方法2 M2÷24÷R-M(4) ÷2384÷R(3) J      内移值计算方法2 M÷2-M(3) ÷240÷R2I     ;切线增长值计算方法2 Abs(÷2- Abs(K) A Rtan(Abs(K÷2)+S+J÷Cos(A)-QTan(A)T

34、 ;计算切线1 Rtan(Abs(K÷2)+I+Q÷Cos(A)-JTan(A)U ;计算切线2 H<0=>H+2H 1A K<0 =>-1A HZF+4 0ZF+5 0ZF+6 F+6F F+6DimZ Rec(T,H)缓和段1 ZD+1-IN ZD+2-JE NZF+1 EZF+2 V-TZF+3 L=0=>Goto A HZF+4 0ZF+5 A÷RZF+6 F+6F F+6DimZ Lbi A圆弧段 0I:0J L=0=>Goto B Pol(L-L(3)÷40÷R2,L2÷6÷R

35、-L(4)÷336÷R(4) Lbi B Rec(I,H+AJ) N+IN E+JE NZF+1 EZF+2 V-T+LZF+3 H+ABZF+4 A÷RZF+5 A÷RZF+6 F+6F F+6DimZ M=0=>Goto C        缓和段2 Rec(2RSin(0.5P÷R),H+AB+AP÷R÷2) N+IN E+JE NZF+1 EZF+2 V-T+L+PZF+3 H+AB+AP÷RZF+4 A÷RZF+

36、5 0ZF+6 F+6F F+6DimZ Lbi C Rec(U,H+K) ZD+1+IZF+1 ZD+2+JZF+2 V-T+L+P+MZF+3 V-T+L+P+M-U+WV D+5D Next H+KZF+4 0ZF+5 0ZF+6 F+9DimZ 1ZF+7 VZF+9   要素终点桩号 Deg FileName：READDAT ZC+32I I=3=>Prog"READDAT3" If G<M Or G>M+L  :Then ;验证要素是否可用，否则重设要素 I=1=>Prog"READDAT

37、1" I=2=>Prog"READDAT2" Prog"READDAT3"   IfEnd FileName:READDAT1 平曲线分段要素读取程序1 ZC+30F 0J Lbi 1           搜索要素 If G<ZF+3:Then -1=ZF=>Goto E:F-6F:-1J:Goto 1 Else 1=ZF+1=>Goto E:-1=J=>Goto 2:1J:F+6F:Goto 1

38、 IfEnd Lbi 2 FZC+30 ZF+1ZC+23当前分段要素N坐标 ZF+2ZC+24当前分段要素E坐标 ZF+3ZC+25当前分段要素起点桩号 ZF+4C+26当前分段要素起点方位角（正北） ZF+5ZC+27当前分段要素起点曲率有左偏负右偏正（注意不半径） ZF+6ZC+28  当前分段要素终点曲率有左偏负右偏正（注意不半径） ZF+9-ZF+3ZC+29  当前分段要素长度 Return Lbi E Cls "ERROR" Stop READDAT2文件请输入内容为以下内容以便测试结果与本文一致。 ÷180J

39、 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.396M:213°3252J(起点方位角)H:0(起点曲率)A:0(终点曲率)R:652.642起点至终点长度L:Goto E:IfEnd If G<30921.038:Then 2850878.162N:513535.768E:30661.038M:213°3252JH:0A:1÷1300R:260L:Goto E:IfEnd If G<31426.367:Then 2850666.478N:513385.004E:30921.038M:219°1639JH:1÷1300A:1÷1300R:505.328L:Goto E:IfEnd If G<31686.367:Then 2850346.483N:512998.015E