隧道断面测量程序(ls1=ls2)

Fixm:

{E，F，K，B，R，L， U， A， V ，W }：E”JDX”F”JDY”K”JD”B”PJ”R

L”LS””U”T” A”A0” V”XZZ”W”WIDE”输入交点x、y坐标，交点桩号、转角值、圆曲线半径、缓和曲线长、切线长、起始方位角、校正值、设计线到隧道中线的距离。

Lb1 0:{X,Y,Z,S}:X”X”Y”Y”Z”Z”S”KWBJ” 输入测量三维坐标及设计开挖半径

PrOg”SUB”……………………………………平曲线文件

PrOg”SQX”……………………………………竖曲线文件

T≤48053.809=>N=460.88-(48200-T)×0.0062: GOtO 1⊿T≤48346.191

=>N=460.88-146.191×0.0062+(T-48053.809)2÷44000: GOtO 1⊿N=

460.88-(T-48200) ×0.00709: Lb1 1:（竖曲线计算高程）

”SJGC”:N◢

“⊿O”：O=Z-(N+1.545) ◢

“⊿R”：√(D2+O2)-R ◢

“⊿Y”：Abs D-√(R2-O2)◢

“⊿H”：Abs O-√(R2-D2)◢

GOtO 0

SUB (数据输入子程序)

G=E-UcosA:H=F-UsinA……………………计算直缓点坐标

C=（X-G）cosA+(Y-H)sinA:C≤0 => PrOg”QZX”

: ≠=>D=L:Q=A:H=U: PrOg”ZBJS”:C≤0 => PrOg”QHHQX”

: ≠=>D=-L:Q=A+B:H=-U: PrOg”ZBJS”:C≤0 => PrOg”YQX”

: ≠=>G=E+Ucos(A+B):H=F+Usin(A+B): C=（X-G）cos（A＋Ｂ）+(Y-H)

sin（A＋Ｂ）:C≤0 => H=-U:PrOg”QHHQX”

: ≠=> PrOg”HZX” ⊿⊿⊿⊿

QZX (前段直线计算)

“LC”:T=K-U+C◢ 输出测点里程

“PC”:D= W -(Y-H)cosA+(X-G)sinA◢ 输出测点到中线距离

QHHQX (缓和曲线段计算)

Defm 8:

Z[1]=(X-E+HcosQ)cosQ+(Y-F+HsinQ)sinQ 将X转换到施工坐标系

Z[2]=-(X-E+HcosQ)sinQ+(Y-F+HsinQ)cosQ 将Y转换到施工坐标系

Z[3]=Z[1]: Lb1 1:Z[4]=Z[1]2÷2RD×180÷π:Z[5]=Z[1]-Z[1]^5÷40(RD)2：Z[6]=Z[1]^3÷6RD-Z[1]^7÷366÷(RD)^3:

B＞0 => Z[6]= - Z[6] ⊿ 切线坐标投影

Z[7]=(Z[3]-Z[5])cosZ[4]-(Z[2]+Z[6])sinZ[4]: 判断投影里程差

Abs Z[7]≥0.001=>Z[1]=Z[1]+Z[7]: GOtO 1⊿

“HLC”:T=K-H+Z[1]：T≥K=>T=T-V⊿T=T◢ 输出测点里程

“HPC”:D=(Z[5]-Z[3])sinZ[4]-(Z[2]+Z[6])cosZ[4]＋W◢ 到中线距离

YQX (圆曲线段计算)

Defm 8:

Z[5]=A+B+P：

Z[1]=(X-G)cos Z[5] +(Y-I)sin Z[5]

Z[2]=-(X-G)sin Z[5] +(Y-I)cos Z[5]

Z[3]=R:B<0=>Z[3]=-R:pol(Z[1],(Z[2]- Z[3])):B<0=>J=J-90:≠ =>

J=J+90 ⊿

“YLC”:T=K+V-L-V+πZ[3] J÷180◢ 输出测点里程

“YPC”:D=√(Z[1]2+(Z[2]-Z[3])2)－Ｒ＋W◢ 输出测点到中线距离

HZX (后段直线计算)

“HZXLC”:T=K+U+C-V◢ 输出测点里程

“HZXPC”:D=W -(Y-H)cos(A+B)+(X-G)sin(A+B)◢ 输出测点到中线距离

ZBJS (缓和曲线段坐标计算)

N=D-D^3÷40R2:O=D2÷6R-D^4÷336R^3:P=90D÷πR:B≤0=>O=-O:P=-P⊿

G=E+(N-H)cosQ-OsinQ:I=F+(N-H)sinQ+OcosQ:…… 缓圆（圆缓）点坐标

C=(X-G)cos(Q+P)+(Y-I)sin(Q+P)………………… 将测量点坐标投影