一、程序功能

本程序由6個主程序、5個次子程序及5個參數子程序組成。主要用于公路測量中坐标正反算，設計任意點高程及橫坡計算,橋涵放樣，路基開挖口及填方坡腳線放樣。程序坐标計算适應于任何線型.

二、源程序

1.主程序1：一般放樣反算程序(①正算坐标、放樣點至置儀點方位角及距離；②反算樁号及距中距離)

程序名:1ZD-XY

Lb1 0:Norm 2

F=1：(正反算判别，F=1正算，F=2反算,也可以改F前加？，改F為變量)

Z[1]=90（與路線右邊夾角）

Prog＂THB＂：F=1=>Goto 1:F=2=>Goto 2

Lb1 1:Ｆix 3:＂ Ｘ=＂：Locate 6,4,Ｘ◢

＂ Ｙ＝＂：Locate 6,4,Ｙ◢

Ｐrog＂3JS”:Goto 0:

Lb1 2:Fix 3:＂ＫＭ＝＂：Locate 6,4,Ｚ◢

＂Ｄ＝＂：Locate 6,4,Ｄ◢

Ｇoto 0

２．主程序２：高程序橫坡程序(設計任意點高程及橫坡)

程序名：2GC

LbI 0:Norm 2

“KM”?Z:?D:

Prog”H”:Fix 3:” H=”:Locate 6,4,H◢

“ I=”: Locate 6,4,I◢

Goto 0

3.主程序3：極坐放樣計算程序(計算放樣點至置儀點方位角及距離)

程序名：3JS

X：Y：

1268．123→K(置儀點X坐标)

2243．545→L（置儀點Y坐标，都是手工輸入,也可以建導線點數據庫子程序,個人認為太麻煩）

Y-L→E：X-K→F：Pol(F,E):IF J<0:Then J 360→J:Int(J) 0.01Int(60Frac(J)) 0.006Frac(60Frac(J)) →J:(不習慣小數點後四位為角度顯示的，也可以用命令J◢DMS◢來直接顯示) Fix 4:” FWJ=”: Locate 6,4,J◢(不習慣小數點後四位為角度顯示的，也可以用命令J◢DMS◢來直接顯示)

Fix 3:” S=”:Locate 6,4,I◢

4．主程序4：涵洞放樣程序（由涵中心樁号計算出各涵角坐标、在主程序3中輸入置儀點坐标後計算放樣點至置儀點方位角及距離)

程序名：4JH-XY

LbI 0:Norm 2

90→Z[1](涵洞中心樁與右邊夾角，手工輸入，也可以修改成前面加？後變為變量)

1→F:Prog”THB”:?L:

Z[2]-Z[1] →E:X Lcos(E) →X:Y Lsin(E) →Y:Fix 3: ＂ Ｘ=＂：Locate 6,4,Ｘ◢

＂ Ｙ＝＂：Locate 6,4,Ｙ◢

Ｐrog＂3JS”:Goto 0:

5．主程序5:路基開挖邊線及填方坡腳線放樣程序（輸入大概樁号及測量坐标、地面标高計算出偏移距離、樁号、距中距離、填挖高度）

程序名：5FBX

LbI 0:Norm 2: 18→DimZ:2→F:90→Z[1]:Prog “THB”:Z:D:”M0”?M:M→Z[4]:D→Z[3]:Prog”6GD”:L→Z[6]:If D<0:Then 0.75-L→D:Goto H:Else L-0.75→D:Goto H:IfEnd

LbI H:Prog”H”:H-0.03-Z[4] →Z[5]:Z[6] →L:If Z[5]<0:Then –z[5] →G:Goto W:Else Z[5] →G:Goto T:

LbI W:Prog “W0”:Z[10] Z[11] →A: If G>A:Then Goto 1:Else If G>Z[10]:Then Goto 2:Else Goto 3:IfEnd:

LbI 1:L Z[12] Z[13] Z[14] (G-A)×Z[9] Z[11]×Z[8] Z[10]×Z[7]:Goto Z:

LbI 2:L Z[12] Z[13] (G-Z[10])×Z[8] Z[10]×Z[7]:Goto Z:

LbI 3:L Z[12] G×Z[7]:Goto z:

LbI T:L 0.5→N:If G>Z[17]:Then (N Z[18] (G-Z[17])×Z[16] Z[17]×Z[15])→S:Goto Z:Else (N G×Z[15])→S:Goto z:

LbI Z:Z[3]→D:Fix 2:Abs(D)-S→T:” L0=”:L Locate 6,4,T◢

＂ＫＭ＝＂：Locate 6,4,Ｚ◢

＂Ｄ＝＂：Locate 6,4,Ｄ◢

“ TW=”: Locate 6,4,Z[5]◢

Goto 0

6．主程序6：路基标準半幅寬度計算程序(對于設計有加寬漸變的有用，如路基寬度無變化，則把此程序直接輸入半幅寬度值至L)

程序名：6GD

Prog “G0”Z-C→E:(B-A)×E/S A→L:L:

7．坐标計算次程序（ＴＨＢ）

程序名：THB

18→DimZ：＂ＫＭ＂?Z：Prog ＂Ｘ０＂

1÷P→Ｃ： (P-R)÷(2HPR) →Ｓ： 180÷π→Ｅ：Ｆ=1=>Goto 1：Ｆ=2=>Goto 2←┘

Lbl 1：?Ｄ： Abs(Ｚ－O) →W：Prog "Ａ"：X：Ｙ：Ｇoto 3

LbI 2:Ｘ：Ｙ：Ｘ→I：Ｙ→J：Prog ＂Ｂ＂：Ｏ＋Ｗ→Z：Ｄ→Ｄ：Ｇoto 3

LbI 3:IF Ｆ＝１Then Ｘ：Ｙ：Else Ｚ：Ｄ

8. 正算子程序(Ａ)

程序名：A

0.1184634425→A： 0.2393143352→B：0.2844444444→N 0.046910077→K： 0.2307653449→L：0.5→M： U W(Acos(G QEKW(C KWＳ)) Bcos(G QELW(C LWＳ)) Ncos(G QEMW(C MWＳ)) Bcos(G QE(1-L)W(C (1-L)WS)) Acos(G QE(1-K)W(C (1-K)WS))) →X：

V W(Asin(G QEKW(C KWＳ)) Bsin(G QELW(C LWＳ)) Nsin(G QEMW(C MWS)) Bsin(G QE(1-L)W(C (1-L)WS)) Asin(G QE(1-K)W(C (1-K)WＳ))) →Y：

G QEW(C WＳ) Ｚ［１］→Ｚ［２］：X Ｄcos（Ｚ［２］）→Ｘ： Y Ｄsin（Ｚ［２］）→Ｙ

9. 反算子程序(Ｂ)

程序名：B

G-90→Ｔ： Abs((Y-V)cos(T)-(X-U)sin(T)) →Ｗ： 0→Ｄ：Lbl 0：Prog "Ａ"： T QEW(C WＳ) →Ｌ： (J-Y)cos(L)-(I-X)sin(L)→Ｄ：IF Abs(Ｄ)<0.01:Then Goto1：Else W Ｄ→W:Goto 0←┘

Lbl 1:0→D：Prog "Ａ"：(J-Y)÷sin(Ｚ［２］) →D:

10．高程計算子程序（H）

程序名：H

Prog “S0”:R:T:C:G:I:C-T→F:Z-F→L:C T→E:G-TI→Q:If T=O:Then Q LI→H:Goto 0:Else If Z<F:Then Q LI→H:Goto 0:Else If Z≤E:Then Q LI L2÷2÷R→H:Goto 0:

LbI 0:H:If D=0:Then Goto I:Else Prog “I”:H V→H:Goto I:LbI I:H:I:

11.高程超高計算程序（I）

程序名：I

Prog”I0”:

W=1=> Goto 0:W=2=>Goto 1:

LbI 0:If L=0:Then Abs(D)×M→V:Goto 2:Else Abs(D)×((N-M)×(Z-C)÷L M)→V:Goto 2:IfEnd:

LbI 1:If L=0:Then Abs(D)×M→V:Goto 2:Else Abs(D)×(((3((Z-C)÷L)2-2((Z-C)÷L)∧(3))×(N-M)) M)→V:Goto 2:IfEnd:

LbI 2:Abs(D)→E:V÷E→I:I(E-K)→V:

12．數據子程序(附後示例)

①程序名：X0（坐标計算要素程序）

If Z≥25900 And Z≤26615.555:Then 25900→O:11587.421→U:1847.983→V:101。09’23.1”→G:715.555→H:1×1045→P: 1×1045→R: 0→Q:Goto 0: IfEnd:

If Z≥26615.555 And Z≤26915.555:Then 26615.555→O:11448.97→U:2550.016→V:101。09’23.1”→G:300→H:1×1045→P: 1800→R:1→Q:Goto 0:IfEnd:

If Z≥26915.555 And Z≤27316.952:Then 26915.555→O:11382.792→U:2842.531→V:105。55’51.83”→G:401.396→H:1800→P: 1800→R:1→Q:Goto 0:IfEnd:

……

LbI 0:O:U:V:G:H:P:R:Q:

程序字母說明：O-線元起點樁号；U-起點X坐标；V-起點Y坐标；G-線元起點樁号切線方位角; Ｈ－線元長度，Ｐ－線元起點曲率半徑，Ｒ－線元終點曲率半徑,Q－線元判别(以道路中線的前進方向（即裡程增大的方向）區分左右；當線元往左偏時，Q=-1；當線元往右偏時，Q=1；當線元為直線時，Q=0)。

②程序名：SO(高程豎曲線要素子程序)

If Z≥25900 And Z≤26157.8:Then -200000→R:157.8→T:26000→C:37.1→G:0.01222→I:Goto 0:IfEnd:

If Z≥26157.8 And Z≤27421.915:Then 300000→R:171.915→T:27250→C:32.65→G:-0.00356→I:Goto 0:IfEnd:

If Z≥27421.915 And Z≤27889.029:Then 1000000→R:134.029→T:27755→C:36.64→G:0.0079→I:Goto 0:IfEnd:

……

LbI 0:R:T:G:C:I:

程序字母說明：C-豎曲線交點樁号；G-交點樁号高程（未豎曲線調整的）；T-豎曲線切線長；R-豎曲線半徑（分正負：凹曲線為正，凸曲線為負）；I－豎曲線前縱坡

③程序名：I0（高程超高參數子程序）

1.5→K:2→W:If D<0:Then Goto L:Else Goto R:IfEnd:

LbI L:If Z≥25900 And Z≤26615.555:Then 0→L:-0.02→M:Goto 0:IfEnd:

If Z≥26615.555 And Z≤26735.555:Then 120→L:26615.555→C:-0.02→M:0.02→N:Goto 0:IfEnd:

If Z≥26735.555 And Z≤26795.555:Then 60→L:26735.555→C:0.02→M:0.04→N:Goto 0:IfEnd

If Z≥26795.555 And Z≤27436.951:Then 0→L:0.04→M: Goto 0:IfEnd

……

LbI R: If Z≥25900 And Z≤26735.555:Then 0→L:-0.02→M: Goto 0:IfEnd

If Z≥26735.555 And Z≤26795.555:Then 60→L:26735.555→C:-0.02→M:-0.04→N:Goto 0:IfEnd

If Z≥26795.555 And Z≤27436.951:Then 0→L:-0.04→M: Goto 0:IfEnd

……

LbI 0 K:W:L:C:M:N:

字程序中字母表示說明：

K－中央分隔帶半幅寬（中樁标高至中樁中線距離，無分隔帶則為0）；W－超高方式參數（W=1為一般直線方式超高，W=2為三次抛物線方式超高）;L－超高漸變段距離（不是漸變段則輸入0）；C-超高漸變段起點樁号（不是漸變段，無需輸入）；M－超高段起點橫坡，N-超高漸變段終點橫坡（不是漸變段，無需輸入）。

④程序名：GO(路基标準半幅寬度參數子程序)

If D<0:Then Goto L:Else Goto R:IfEnd:

LbI L:If Z≥25900 And Z≤28110.727:Then Z→C:17.25→A:Goto 0:IfEnd:

If Z≥28110.727 And Z≤28200.727:Then 28110.727→C:17.25→A:21.25→B:90→S:Goto 0:IfEnd

……

LbI R: If Z≥25900 And Z≤27927.478:Then Z→C:17.25→A: Goto 0:IfEnd

If Z≥27927.478 And Z≤28172:Then 27927.478→C:17.25→A:27.031→B:244.522→S:Goto 0:IfEnd

……

LbI 0 C:A:B:S:

字程序中字母表示說明：

C－寬度漸變段起點樁号（不是漸變段C=Z）；A－寬度漸變段起點寬度; B－寬度漸變段終點寬度（不是漸變段，無需輸入）；S-寬度漸變段距離（不是漸變段，無需輸入）。

⑤程序名：WO(路基填挖邊坡參數子程序)

0.75→Z[7]:1→Z[8]:1.25→Z[9]:10→Z[10]:10→Z[11]:2.6→Z[12]:2→Z[13]:2→Z[14]:

1.5→Z[15]:1.75→Z[16]:8→Z[17]:2→Z[18]

字程序中字母表示說明：

Z[7]－挖方第一階邊坡邊率；Z[8]－挖方第二階邊坡邊率；Z[9]－挖方第三階邊坡邊率；Z[10]－挖方第一階高度；Z[11]－挖方第二階高度；Z[12]－挖方路基碎落台及水溝寬度；Z[13]－挖方第二台階平台寬度；Z[14]－挖方第三台階平台寬度；

Z[15]－填方第一階邊坡邊率；Z[16]－填方第二階邊坡邊率; Z[17]－填方第一階高度；Z[18]-填方第二台階平台寬度。（注：本程序隻做出挖方三個台階，填方二個台階，如需增加，先需在程序THB和5FBZ程序中增加變量。再如一标段有坡率及坡高不一樣時，可以照其它參數程序一樣，用判别語句。）

三、使用說明

1、規定

(1) Z“KM”為所求點樁号，反算時為輸入大概樁号

(2) 當所求點位于中線時，Ｄ=0；當位于中線左鍘時，Ｄ取負值；當位于中線中線右

側時，Ｄ取正值。

(3) 當線元為直線時，其起點、止點的曲率半徑為無窮大，以10的45次代替。

(4) 當線元為圓曲線時，無論其起點、止點與什麼線元相接，其曲率半徑均等于圓

弧的半徑。

(5) 當線元為完整緩和曲線時，起點與直線相接時，曲率半徑為無窮大，以10的45次代替；與圓曲線相接時，曲率半徑等于圓曲線的半徑。止點與直線相接時，曲率半

徑為無窮大，以10的45次代替；與圓曲線相接時，曲率半徑等于圓曲線的半徑。

(6) 當線元為非完整緩和曲線時，起點與直線相接時，曲率半徑等于設計規定的

值；與圓曲線相接時，曲率半徑等于圓曲線的半徑。止點與直線相接時，曲率半徑等于設計規定的值；與圓曲線相接時，曲率半徑等于圓曲線的半徑。

2、輸入與顯示說明

輸入部分：

1. Ｆ＝１：ＺＤ => XY

2.Ｆ＝２： XY = > ＺＤ

Ｆ ? 選擇計算方式，輸入1表示進行由裡程、邊距計算坐标 ；輸入2表示由坐标反算裡程和邊距。Ｚ［１］為線元邊樁與中線右夾角。

ＫＭ？正算時所求點的裡程（反算輸入大概樁号）

Ｄ ？正算時所求點距中線的邊距(左側取負，值右側取正值，在中線上取零)

X0 ？反算時所求點的X坐标（放樣程序中實測X坐标）

Y0 ？反算時所求點的Y坐标（放樣程序中實測Y坐标）

顯示部分：

L? 涵洞放樣程序中涵距中心樁号前後距離,前為正,後為負.

M0? 邊坡放樣程序中原地面标高輸入

X=*** 正算時，計算得出的所求點的X坐标

Y=*** 正算時，計算得出的所求點的Y坐标

FWJ=*** 正算時，計算得出的所求點的至置儀點方位角（在3JS程序中應輸入X、Y坐标值于K、L）

S=*** 正算時，計算得出的所求點的至置儀點距離

ＫＭ=*** 反算時，計算得出的所求點的裡程

Ｄ=*** 反算時，計算得出的所求點的邊距

H=*** 所求點位置設計路面頂标高

I=*** 所求點位置設計路面橫坡

L0=*** 邊樁放樣程序中實測點至設計邊坡點距離，正向内移