程序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×120000))    ←┘ 

S ≥ 4600 => S ≤ 4900  => H = 124.126 －0.020 ×(4900－S) - ((4900－S )²  ÷ (2×120000))   ←┘ 

S ≥ 4900  => S ≤ 5322  => H = 124.126 + 0.020 ×(S －4900)    ←┘ 

（六）、SM 右线 ：隧道半径计算子程序

{ M }： tan^-1（（ M - H）÷Abs（Z－ 6.7））)≥θ

=> “ R₁= ”：R =√（（ M - H）)²＋（Abs（Z－6.7））²）◢  

 ≠> “R₂  =” ：R =√（（M - H －1.061）²）＋ （（Abs （Z －6.7））-1.061）²

◢  

（七）、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