任意坐标正反算及隧道超欠挖程序
1-JIN YU(主程序)
Lbl 4:"1→ZS,2→FS,3→CQW"?N
N=1=>Goto 1:N=2=>Goto 2:N=3=>Goto 3
Lbl 1: “K=”?S：“P=”?Z：Prog “SJ?-PM”：Abs(S-O) → W：Prog "SUB1-ZS"：“X="：Locate4,4,X:"Y="：Locate4,4,Y:F-90→F:S→ K:Prog“SJ?-GC”：“H=”：Locat e4,4,H◢
Goto 4
Lbl 2:“X=”?X：“Y=”？Y:Prog“SJ?-PM”：X→ I：Y→ J：Prog "SUB2-FS"：O+ W→S: “K="：S ◢
“P="：Z◢
S→ K:Prog“SJ?-GC” ：“H=”：H◢
Goto 4
Lbl 3: “X=”?X：“Y=”？Y:Prog“SJ?-PM”：X→ I：Y→ J：Prog "SUB2-FS"：O +W→S: “K="：S ◢
“P="：Z◢
S→ K:Prog“SJ?-GC” ：“H=”：H◢
Prog“SJ-CQW” ◢
Goto4
SJ?-PM（子程序名-平面线形数据库）
If S ≥59227.681（线元起点里程）:Then 2599818.013→U（线元起点X坐标）：49 6887.918→V（线元起点Y坐标）：59227.681→O（线元起点里程）：208028’39”→G（线元起点方位角）：100 →H（线元长度）：1×1045→P（线元起点曲率半径）：1×1045→R（线元终点曲率半径）：0 →Q（线元左右偏标志：左-1右1）：IfEnd If S ≥59327.681（线元起点里程）:Then 2599730.112→U（线元起点X坐标）：49 6840.237→V（线元起点Y坐标）：59327.681→O（线元起点里程）：208028’39”→G（线元起点方位角）：90 →H（线元长度）：1×1045→P（线元起点曲率半径）：750→R（线元终点曲率半径）：-1 →Q（线元左右偏标志：左-1右1）：IfEnd
……
……
SJ-GC?（子程序名-竖数据库）
If S＜变坡终点里程And S≥变坡起点里程：Then大里程坡度→A：小里程坡度→B：变坡点里程→O：变坡点高度→G：半径→R：Prog“SUB3-GC”:IfEnd
(注：1.有多个竖曲线，依照上面的依次变更，每多一个，就增加一个。每次只需要修改以上的数据或增加一个判断，子程序不用变动。2.如整条线只有一个纵坡比喻为2﹪，那么程序应为If S＜终点里程And S≥起点里程：Then 0.02→A：0.02→B：起点里程→O：起点高程→G：1×1045→R：Prog“SUB3-GC”:IfEnd)
SJ-CQW (子程序名-隧道数据库)
“H1=”?F:”P=”?Z
If F-H＞5.845:Then √((Abs (Z))2+(F-H-0.45)2 )-7.45→W: IfEnd:
If F-H≤5.845 And F-H＞-1.169: Then √((Abs (Z)-1)2+( Abs (F-H-1.5)) 2)-6→W : IfEnd:
If F-H≤-1.169: Then √((Abs (Z)) 2+(H-F+14.05) 2)-16.5→W: IfEnd:
…
…
“+CQ,-CQ=”: W
SUB1-ZS(正算子程序)?
1÷P→ C：(P-R)÷(2HPR) → D：180÷π→ E：0.1739274226→ A：0.32607257 74→ B：0.0694318442→ K：0.3300094782→ L：1-L→ F：1-K→ M：U+W(A cos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QEFW(C+FWD))+Aco s(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+Zcos(F)→ X：Y+Zsin(F)→ Y
SUB2-FS(反算子程序)?
G-90→T：Abs((Y-V)cos(T)-(X-U)sin（T）) → W：0→ Z：Lbl 0：Prog "SUB1-ZS"：T+QEW(C+WD) →L：(J-Y)cos(L)-(I-X)sin(L) → Z：If Abs(Z)<1E-6(注：此处1E-6可输入0.000001)：then Goto1：Else W+Z → W：Goto 0：IfEnd
Lbl 1：0→ Z：Prog "SUB1-ZS"：(J-Y)÷sin(F)→ Z
SUB3-GC（高程子程序）
O-S→L:A-B→W:Abs(R*W÷2) →T:O-T→M:O+T→P
If S≤M:ThenG-L*B→H:Goto5:Else If S≤O:Then Goto3:Else If S≤P:Then Goto4:If End: IfEnd: IfEnd
Lbi3
If W＞O：Then G+(M-S)2÷2÷R-L*B→H:Goto5: Else If W＜O: Then G-(M-S)2÷2÷R-L*B→H:Goto5: If End: IfEnd
Lbi4
If W＞O：Then G+( S-P)2÷2÷R-L*A→H:Goto5: Else If W＜O: Then G-( S-P)2÷2÷R-L*A→H:Goto5: IfEnd: IfEnd
Lbi5
H→H: Return
说明：(正算1秒，反算和超欠挖5秒完成)
所有程序名。不同线路，只需改动SJ-PM，SJ-GC，SJ-CQW三个子程序里的内容，其它不变。
1为正算，2为反算，3为超欠挖.
K= 正算时，输入所求点里程。反算时得出里程结果
X=，Y= 正算时得出结果。反算时输入实测坐标
P= 正算时输入偏距。反算时得出偏距
H= 正算时得出高程。反算时得出高程
H1= 超欠挖时输入实测高程
+CQ,-CQ= 得出超欠挖值
注，超欠挖输入时显示里程，编距，设计高程，不需输入。只输入实测量高程H1 谢谢发贴：+302 分【测量空间】https://www.doczj.com/doc/971062468.html,[0
复化辛卜生公式及公路中线坐标计算中的应用
时间:2010-06-05 09:12:33 来源:本站作者:未知我要投稿我要收藏投稿指南
在公路中线坐标计算中，我们通常采用切线支距公式来计算曲线上各点的坐标。但当在不同的曲线上计算时就需用不同的计算公式，这为计算也带来不便。在设有缓和曲线的圆曲线半径较小或是卵形曲线上的坐标计算时，如公式选用不当就会出现较大计算误差，即便是能对切线支距公式进行多项展开，也会增加计算的难度。而用复化辛卜生公式不仅能解决不同曲线线型或直线上的坐标计算问题，而且用复化辛卜生公式计算完全是可逆的（即：可顺前进方向也可逆向计算），尤其在计算第二缓和曲线和卵形曲线时显得尤为方便。
用辛卜生公式计算坐标的精度可由人为或程序自行判断，其计算结果完全能保证坐标计算的精度要求。因此，可以说复化辛卜生公式是一个计算公路中线坐标的万能公式。下面本人就该公式在公路中线坐标计算中的具体应用进行实例解析。
一、复化辛卜生公式
式中：
H=（Z
i －Z
A
）/n
(公式2)
(公式3)
Zi—待求点桩号
Z
A
—曲线元起点桩号
Z
B
—曲线元终点桩号
ρ
A
—曲线元起点曲率
ρ
B
—曲线元终点曲率
a
i
曲线上任意一点处切线方位角的计算方法有以下三种方法：
1．利用公式（3）求得曲率代入公式（2）计算
2．利用曲线元上已知起点和终点曲率用内插法求得曲率代入公式（2）计算3．利用切线角公式计算
二、算例
例：已知雅（安）攀（枝花）高速公路西昌西宁立交A匝道一卵形曲线（卵形曲线相关参数见图一，其计算略。），相关设计数据见下表。现用辛卜生公式来计算卵形曲线中桩坐标。
图一
已知相关设计数据见下表：
主点桩号坐标
（m）
切线方位角
（θ）
X Y ° ’ ”
ZH
AK0+090
9987.403 10059.378 92 17 26.2 HY1
AK0+160
9968.981 10125.341 132 23 51.6 YH1
AK0+223.715
9910.603 10136.791 205 24 33.6 HY2
AK0+271.881
9880.438 10100.904 251 24 18.5 YH2
AK0+384.032
9922.316 10007.909 337 04 54.2 HZ
AK0+444.032
9981.363 10000.000 0 00 00
（一）由+271.881推算Zi=+223.715的坐标，n取2等分
用公式(3)、公式（2）计算+247.798处曲线及方位角:
ρ
=1÷75+(1÷50-1÷75)(247.798-271.881) ÷(223.715-271.881) +247.798
=0.01666666666666667
a
=71°24’18.5”
+247.798
+(0.016666667+1÷75)(247.798-271.881)×180÷π÷2
=50°42’26.37”
其它各点依次代入公式计算，结果见下表:
曲率及切线方位角计算表