隧道三心圆超欠挖程序修改

正算主程序(ZS)：

3→DimZ

Lb1 0：CLS:“K=”?S：“BJ=”?Z：Prog “PM-SJ”： Abs(S-O) → W：Prog \"SUB1\"：CLS: Locate 1,1,\"XS=\"：Locate 4,1,X: Locate 1,2,\"YS=\"：Locate 4,2,Y:Fix 4: Locate 1,3,”FWJ=“：F-90→F：Locate 5,3,Int(F)+int(Frec(F)*60)/100+

Frec((Frec(F)*60)*60)/10000:Prog“SQX”： Locate 1,4,“H=”： Locate 3,4,H◢ Goto 0

反算主程序(FS)录入X,Y及实测高程

3→DimZ “JS-K”?S

Lb1 0:CLS:“X=”?X：“Y=”？Y：Z[1] →Z:“GC=”?Z:Z→ Z[1]:Prog“PM-SJ”： X→ I： Y→ J：Prog \"SUB2\"：CLS: Locate 1,1, \"K=\"：O+W→S: Locate 3,1,S: Locate 1,2,\"Z=\"：Locate 3,2,Z :Prog“SQX” ：Locate 1,3,“H+2.95=”:Locate 8,3,H+2.95:Locate 1,4,“H=”：Locate 3,4,H◢ Prog “CQW”:Goto 0 2.95上导起拱线1米位置

正算子程序(SUB1)

1÷P→ C： (P-R)÷(2HPR) → D： 180÷π→ E： 0.1739274226→ A： 0.3260725774→ B： 0.0694318442→ K： 0.3300094782→ L： 1-L→ F：1-K→ M：

U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QEFW(C+FWD))+Acos(G+QEMW(C+MWD))) → X：

V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEFW(C+FWD))+Asin(G+QEMW(C+MWD))) → Y： G+QEW(C+WD)+90→ F： X+ZcosF→ X： Y+ZsinF→ Y 反算子程序(SUB2)

Abs((Y-V)cos(G-90)-（X-U）sin（G-90）) → W：Lbl 0：0→ Z :Prog \"SUB1\"： (J-Y)cos(F-180)-(I-X)sin(F-180) → Z：if AbsZ<0.0001：then Goto 1：Else W+Z → W：Goto 0：IfEnd

Lbl 1：(J-Y)÷sin(F)→ Z 子程序（平面线形数据库）PM

If (S -45798.226线元起点里程)(S-45798.226-200终点里程)≤0 ：Then

[[起点X坐标，起点Y坐标，起点里程，起点方位角，线元长度,起点半径,终点半径，左右转标志]] →Mat A Ifend

说明：若是直线或者由变径变至直线，则直线处的半径为10^45 有多少个线元输入多少个,直到线元全部录入

隧道3心圆放样数据库（CQW-SJ）

If (H-1.95)（H-10（^45））≤0 有多少段就写多少段，以右半部分为基础计算 Then 1.95 → I:0→ J:4.780→ R I，J，R及分段高度值。注意I是圆心到 Ifend 基线的高度，J是相对隧道中线的偏距， If (H-1.95)(H+10^(45)) ≤0 右侧为正值，R为隧道内设计半径

Then 1.95→I:4.-1.5→J:6.28→R H-1.95表示圆心到基线高度。 H+10的45次方

三心圆（CQW-SJ1）

If (H-5.33)(H-10^（45）) ≤0 Then 0.323→I：0→J：7.078→R

Ifend 有多少段就写多少段，以右半部分为基础计算 If （H-1.95）（H-5.33）≤0 I，J，R及分段高度值。注意I是圆心到 Then 1.95→I:1.625→J:4.78→R 基线的高度，J是相对隧道中线的偏距， Ifend 右侧为正值，R为隧道内设计半径

If (H-1.95)(H+10^(45)) ≤0 H-1.95表示圆心到基线高度。 H+10的45次方 Then 1.95→I:0.126→J:6.28→R Ifend

隧道3心圆放样主程序（CQW）

Z[1]-H→ H:prog “COW-SJ” Cls

Z[2] →A: “JGC=” ?A :A→Z[2] 录入结构层厚度

Z[3] →A: “C-Z=” ?A :A →Z[3] 录入隧道中线到测量线偏距，左侧时为负中线时为0 CLS:Locate 1,1,”K=”:Locate 3,1,S

Locate 1,2,”BJ=” :Locate 4,2,Z: Locate 1,3,”H-?GC=” :Locate 7,3,H: Locate 1,4,”R0=”:Locate 4,4,R◢

CLS:Locate 1,1,”OH=”:Locate 4,1,I：Locate 1,2,”OD=” :Locate 4,2,J: Locate 1,3,”R=” :Locate 3,3,((H-I)2+（Abs(Z -Z[3]) -J)2)^(0.5） Locate 1,4,”JGC=” :Locate 6,4,Z[2]◢

CLS:Locate 1,1,”+C-Q=”:Locate 7,1, ((H-I)2+（Abs(Z -Z[3]) -J)2)^(0.5）-R-Z[2]：Locate 1,2,”BJ=” :Locate 4,2,Z: Locate 1,3,”J-?GC=” :Locate 7,3,H: Locate 1,4,”R=” :Locate 3,4,((H-I)2+（Abs(Z -Z[3]) -J)2)^(0.5）◢

K=表桩号 BJ=边距 DH=表测点到设计线高差 R0为设计半径 OH 所在圆弧圆心高度，OD 所在圆弧相对圆心高度 JGC 表标结构 层厚度 R=实测半径 ，+CQ结果为正时为超，负时为欠，本程序可以计算多心圆

子程序（平面线型要素转换）PM-SJ Prog “PM”

Mat A[1,1] →U: Mat A[1,2] →V: Mat A[1,3] →O: Mat A[1,4] →G: Mat A[1,5] →H: Mat A[1,6] →P: Mat A[1,7] →R: Mat A[1,8] →Q 子程序（竖曲线要素库）SQX-SJ

If (S -竖曲线能计算起点里程)(S-竖曲线能计算终点里程)≤0 ：Then [[变坡点里程，变坡点高程，变坡半径，前坡，后坡]] →Mat B Ifend

If (S-0)(S-28800) ≤0:Then

[[27490,1783.002,4500,-0.01665,0.00795]] →Mat B Ifend

If (S-2800)(S-29000) ≤0 Then

[[28340,1789.758,20000,0.00795,-0.0075]] →Mat B Ifend

子程序（竖曲线要素库）SQX

Prog “SQX-SJ”

Mat B [1,2] +(S-Mat B [1,1])* Mat B [1,4] →H

If S> Mat B [1,1] :then Mat B [1,2] +(S-Mat B [1,1])* Mat B [1,5] →H:Ifend Abs(Mat B [1,3]*( Mat B [1,4]- Mat B [1,5]))/2 →T If Abs(S- Mat B [1,1])H-0.5*(Abs(S- Mat B [1,1])-T)^(2)*Abs（Mat B [1,4]- Mat B [1,5]-10^(-9)）/ （Mat B [1,3]*（Mat B [1,4] -Mat B [1,5]）-10^(-9)）→H Ifend 字母表示意思

”FWJ=“：F-90→F：Locate 5,3,Locate表示定位显示， 5表示方位角屏幕由左到右第5个格显示，3表示右上到下第3格显示