程序2020.12.20
4850综合程序 更新版
(一)、SUB3520 右线:
主程序 ( 综合曲线 里程、边距、坐标计算 )
“1.SZ => XY” :“XY =>2. SZ” : { N } : N=1=>Goto 1 : ≠> Go to 2 ◢
Lbl 1: {SZ} : SZ : Prog “ SUBθ-Y ” : D=(P-C) ÷(2PC(B-O )) :W=Abs(S-O) :Prog “ SUB1 ” :“ XS= ”: X=X ◢ “ YS= ” : Y=Y ◢ “FS= ” : F-90 ◢
Goto 1 ←┘
Lbl 2 : {SXY} :SXY : Prog “??? ” : D=(P-C) ÷(2PC(B-O )) : I=X : J=Y
: Prog“SUB2 ” : S = O+W :“ S = ” : S = S ◢
“ Z = ”:Z = Z ◢
Prog “ SH3403-Z ” ◢ (竖曲线)
Prog “ HAA ” ◢ (路面横坡)
Prog "HB" ◢ (电缆沟盖板顶高程)
Prog "SM-Z " ◢ (隧道半径计算)
Go to 2 ←┘
(二)、SUB1 子程序 (正算)
Z【99】= 0.1739274226 :Z【98】= 0.3260725774:
Z【97】= 0.0694318442 :Z【96】= 0.3300094782
:X=U十W(Z【99】 cos(G+(180÷P)Z【97】QW((1÷P)十 Z【97】WD))
+ Z【98】 cos(G+(180÷P)Q W Z【96】((1÷P)+ WD Z【96】))
十 Z【98】 cos(G+(180÷P)(1-Z【96】)Q W((1÷P)+(1-Z【96】)WD))
十 Z【99】 cos(G+(180÷∏)(1-Z【97】)Q W((1÷P)+(1-Z【97】)WD))) :Y= V十W(Z【99】 sin(G+(180÷∏)Z【97】Q W((1÷P)+ Z【97】 WD))
+ Z【98】 sin(G+(180÷P) QW Z【96】((1÷P)十WD Z【96】))
+ Z【98】 sin(G+(180÷∏)(1-Z【96】)QW((1÷P)十(1-Z【96】)WD))
+ Z【99】 sin(G+(180÷∏)(1-Z【97】) QW((1÷P)十(1-Z【97】)WD)))
:F =(G 十(180÷∏)Q W((1÷P)十WD))十90 :X =X十ZcosF :Y=Y十ZsinF←┘
(三)、SUB2子程序 ( 反算 )
T =G-90 :W =Abs((Y-V)cos(T)-(X-U)sin(T) :Z =θ:
Lblθ:Prog “SUB1” :L= T十(180÷∏)Q W((1÷P)+W D):Z= (J-Y)cos(L)-(I-X)sin(L) : Abs Z < 1÷0.000001 => Goto 1 :≠> W = W十Z :Goto θ
Lbl 1:Z =θ:Prog“ SUB1” :Z =( J-Y ) ÷sin F ←┘
(四)、SUBθ-Y 子程序(右线:数据库)
S ≥ 2378.450 => S ≤3564.857 => U = 2567960.787
: V = 478640.063 : G = 79° 37′10.4 ″: P = 3301: C = 3301
: O = 2378.45 :B = 3564.857 : Q =-1 ←┘(弯箭头)
S ≥3564.857 => S ≤ 3709.783 => U = 2568377.450
: V = 479744.086: G = 59° 01′37.1 ″: P = 3301:C = 3600
: O = 3564.857: B = 3709.783: Q = -1 ←┘(弯箭头)
S ≥3709.783 => S ≤ 5480.327 => U = 2568454.663
: V = 479866.718: G = 56° 36′57.39″: P = 3600:C = 3600
: O = 3709.783: B = 5480.327: Q = -1 ←┘(弯箭头)
S ≥5480.327 => S ≤ 5880.327 => U = 2569746.383
: V = 481051.456: G =28° 26′12.69″: P = 3600:C =10∧45
: O = 5480.327: B = 5880.327: Q = -1 ←┘(弯箭头)
(五)、SH 3520 - Y ( 竖曲线 程序 )
S ≥ 3400 => S ≤ 4300 => H = 88.126 +0.025 ×(S-3400 ) ←┘
S ≥ 4300 => S ≤ 4600 => H = 110.626 +0.0025 ×(S-4300 ) - (( S-4300 )2 ÷ (2×120000)) ←┘
S ≥ 4600 => S ≤ 4900 => H = 124.126 -0.020 ×(4900-S) - ((4900-S )2 ÷ (2×120000)) ←┘
S ≥ 4900 => S ≤ 5322 => H = 124.126 + 0.020 ×(S -4900) ←┘
(六)、SM 右线 :隧道半径计算子程序
{ M }: tan-1(( M - H)÷Abs(Z- 6.7)))≥θ
=> “ R1= ”:R =√(( M - H))2+(Abs(Z-6.7))2)◢
≠> “R2 =” :R =√((M - H -1.061)2)+ ((Abs (Z -6.7))-1.061)2
◢
(七)、HAA: 子程序( 路面横坡计算程序 )
“ HA = ” :(H-0.300)+Z×0.020 ◢
(八)、HB: 子程序( 电缆沟盖板顶高程 )
“ HB = ” :H + 0.64 ◢
输入与显示说明 ( 待修改)
(一)数据库输入说明
X0 :线元起点(切点)的坐标(X); Y0 :线元起点(切点)的坐标(Y );
0 :线元起点(切点)里程 ; B?:所求线元的终点桩号;
F0 :线元起点(切线)的方位角;
P :线元(曲线)起点的曲率半径 ,(当线元为完整缓和曲线时)。
①起点与直线相接时,曲率半径为无穷大,以10的45次代替;
②起点与圆曲线相接时,曲率半径等于圆曲线的半径。
C :线元终点的曲率半径。(当线元为完整缓和曲线时)。
①终点与直线相接时,曲率半径为无穷大,以10的45次代替;
②终点与圆曲线相接时,曲率半径等于圆曲线的半径。
Q : 线 元转角方向(左偏Q= - 1,右偏Q= 1,直线段Q= 0);
(二) 正算输入说明
N? : 输入1 ( 用边距和桩号 计算 ____ 坐标X,Y ).
S ?: 所求点的桩号.
Z ?:所求点距线路中线的边距(左侧取负值, 右侧取正值,在中线上取零)
(三) 反算输入说明
N ?: 输入2 . (用X,Y的坐标反算桩号和边距)。
S ?: 所求点的接近桩号.
X ?:反算时所求点的实测(X)坐标;
Y ?:反算时所求点的实测(Y)坐标. FS: 计算得出的(所求点的)切线方位角.
M ?:为所求点的实测高程。
(四)显示部分:
① 正算时: XS= : 计算得出的(所求点)X坐标;
YS= : 计算得出的(所求点)Y坐标;
② 反算时: S= : ,所求点的里程 ;
Z= : 反算时,所求点至线路中桩的边距。
H=: 为设计高程。 R1和R2为隧道道洞身半径. 2020.12.20