4850隧道测量超欠挖程序

1. WGZZ（计算线路中桩坐标及高程）：

X=2808751.636+(K-1906717.28)cos133.5483161+2.5cos223.5483161▲

Y=441813.404+(K-1906717.28）cos133.5483161+2.5cos223.5483161▲

K≤1917014 => prog“WGA”：≠＝> K≤1917200 => prog“WGB”： ≠＝> K≤1917386 => prog“WGC”： ≠＝> K≤1918293.25 => prog“WGD”： ≠＝>K≤1918400 => prog“WGE”： ≠＝> K≤1918506.75 => prog“WGF”： ≠＝> K≤1920417.5=> prog“WGD”： ≠＝> “END” ▲

（2808751.636，441813.404）是直线左线起点坐标，1906717.28为直线起点里程，133.5483161为直线方位角。2.5为线路中线和隧道中线的距离，223.5483161为直线法线方位角。 加下划线的是直线.段的高程程序 2. WGZk（直线坐标反算里程）： A“X1”：B“Y1”：I=0：J=0

POL（（A-2808451.636），（B-441813.404）） K=ICOS（J-133°32′54″）+1906717.28▲ X1，Y1为实测坐标。133°32′54″直线方位角 (2808451.636, 441813.404) 是直线左线起点坐标 3. WGZC（坐标反算宽度）：（主程序）

A“X1”：B“Y1”：C“Z1”（实测坐标与高程） Prog “WGZK ” Prog “WGZZ”

O=C-Z▲(O表示实测高程与内轨顶高差，Z 表示内轨顶高) J=Z+2.43▲（圆心高程）

D=√((A-X)2+(B-Y)2)▲（计算出实测值和隧道中线的距离）。 R=√ (D2+(C-Z-2.43)2)-7.50▲7.5表示半径。（超欠挖，“+”表示超，“-”表示欠） 4. WGA（子程序）：

Z=194.19+（1917014-K）* 0.0155▲ 坡率为0.0155的高程程序 5. WGB（子程序）：

Z=194.19-0.0155（K-1917014）+（K-1914014）2/30000*2▲ 30000为竖曲线半径，194.19为起坡点高程, 1917014起坡点里程 6. WGC（子程序）：

Z=190.7304+0.0031（1917386- K）+（1917386- K）2/60000▲ 60000竖曲线半径的2倍, 190.7304竖曲线终点高程, 1917386竖曲线终点里程 7.WGD （子程序）：

Z=190.7304-0.0031（K-1917386）▲ 8. WGE （子程序）：

Z=187.9179-0.0031（K-1918293.25）+（K-1918293.25）2/70000

▲

9. WGF （子程序）：

Z=187.90725-0.003（1918506.75-K）+（1918506.75-K）2/70000▲

10. WGG （子程序）：

Z=187.90725-0.003（K -1918506.75）▲

Z表示实测高程，其中用竖曲线计算公式和相应的里程 缓和曲线计算；

11，GHZ：（缓和曲线中线坐标及高程）

L=K-1919885.166 （1919885.166为直缓点里程）。

X=2799679.41+(L-L5/7396000/108+L9/1016/3456/4304)cos133.5483161-(L3/2580/104-L7/10/336/4303)sin133.5483161+2.5cos(133.5483161+9

120+L2/860/10000180/∏)▲

Y=451357.404+(L-L5/108/7396000+L9*3456*1016*430)sin133.548316

41-(L3/104/2580-L7/336/1012/4303)

cos133.5483161+2.5sin

(133.5483161+90+L2/860/10000*180/∏)▲

（2799679.41，451357.404）为直缓点左线坐标，里面关系到坐标系的转换。

Z=187.90725+0.003*(k-1918506.75) ▲ 12.GHK(坐标反算里程) A“X1”：B“Y1” Prog“GHZ”

Pol（（A-X），（B-Y）：J≤0＝> J=J+360

Q=133.5483161+L2/20000/430/∏*180：Q≥360Q=Q-360

133.5483161为直线段方位角，20000为曲线半径的2倍，430为缓和曲线长，线路左偏为 “-“,右偏为 “+”. U=ICOS（Q-J）▲ U≥3＝> “ERROR”  U≤-3＝> “ERROR”  K=u+k▲

13．GHC；坐标反算宽度断面及超欠情况。（主程序） A“X1”：B“X2”：C“Z1” Prog“GHK” Prog“GHZ”

J=Z+2.43▲Z 表示内轨顶高 H=Z-0.2▲

D=√((A-X)2+(B-Y)2)▲反算的实测点到中线的距离 R=(D2+(C-Z-2.43)2)-7.50▲超欠挖

(D2+(Z-C+13.670)2)-16.22▲≠＝>C≤Z+0.125

▲

≠

＝

>

C≤Z-1.809＝﹥R=＝R=﹥

R=

((D-3.932)2+(C-Z-1.101)2)-3.05

(D2+((C-Z-2.43)2)-7.10▲≠＝> “NED” ▲

圆曲线段计算坐标及里程； 15， GYZ：线路中桩坐标和高程：

X=2792282.9881+9997.5COS（44.78016+（K-1920315.166）*180/10000/

∏）▲

Y=444623.0433+9997.5SIN（44.78016+（K-1920315.166）*180/10000/∏）▲

（2792282.9881，444623.0433）为圆心坐标，44.78016是连接圆心与缓圆点直线的方位角。9997.5为圆心到隧道中线的距离. 10000为设计半径.

K≤1920417.5=> Prog“WGG”： ≠＝> Prog“WGG” 14，子程序：WGH

Z=193.6395+（K-1920417.5）* 0.003-（K-1920417.5）2/60000▲ 圆曲线坐标反算里程； 16．GYK A“X1”：B“Y1”

Pol（（A-2792282.7881），（B-444623.0433）） K=（J-44.78016）*∏*10000/180+1920315.166▲ 17 ，GYC， 圆曲线计算主程序 A“X1”：B“Y1”：C“Z1” Prog“GYK” Prog“GYZ” O=Z+2.43▲

D=√((A-X)2+(B-Y)2)▲ R=

(D2+(C-Z-2.43)2)-7.5▲