测量计算程序汇篇.docVIP

测量计算程序 一，平曲线计算程序 1，主程序【WXA】 “1,SZ=＞XY” ： “2,XY=＞SZ” ：N ：U“X0” ：V“Y0” ： O“S0” ：G “F0” ： H“LS” ：P“R0” ：R“RN” ：Q ：C=1÷P:D=(P-R)÷(2HPR)：E=180÷∏：N=1=＞Goto 1：≠＞Goto 2 ⊿ LbI 1：｛SZ｝：S“DK”：Z“ZJ”：W=Abs（S-O）:Prog“SUB1” ： X“XP”=X◢ Y“YP”=Y◢ Goto 1 LbI 2：｛XY｝：X“XP ” ：Y“YP” ：I=X ： J=Y ： Prog“SUB2” ：S“DK”=O+W◢ Z“ZJ”=Z◢ Prog“SUB3” Goto 2 2,子程序【SUB1】 A=0.1739274226 ： B=0.3260725774 ：K=0.0694318442 ： L=0.3300094782 ： F=1-L ： M=1-K ： X=U+W(Acos(G+QEKW(C+KWD)) +Bcos(G+QELW(C+LWD))+Bcos(G+QEFW(C+FWD))+Acos(G+QEMW(C+MWD)))： Y=V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEFW(C+FWD))+Asin(G+QEMW(C+MWD)))：F=G+QEW(C+WD)+90 ： X=X+ZcosF ： Y=Y+ZsinF 3,子程序【SUB2】 T=G-90 ： W=Abs((Y-V)cosT-(X-U)sinT) ： Z=0 ： LbI 0 ： Prog“SUB1” ： L=T+QEW(C+WD) ： Z=(J-Y)cosL-(I-X)sinL ： AbsZ＜1E-6=＞Goto 1：≠＞W=W+Z：Goto 0 ⊿ LbI 1：Z=0：Prog“SUB1” ：Z=（J-Y）÷sinF 4, 子程序【SUB3】 K=AbsZ : ｛Z｝：Z“ZP” S≥6009.2 =＞A=36.44-0.014(6250-S)-(S-6009.2)∧2÷(2*20000): ≠＞A=36.44-0.014(6250-S) ⊿ Z≥A+5.421=＞B“R1”=√(K∧2+(Z-A+1.2) ∧2) ◢≠＞B“R2” =√((K-3.31) ∧2+(Z-A-1.56) ∧2) ◢ 使用说明 线路前进方向分左右，线路左偏转时Q=-1右偏转时Q=1；当线路为直线时Q=0。 所求点位于线路左侧时ZJ取负值右侧取正值中线上取零。 直线上的点的曲率半径为无穷大以10的45次方代替。曲线上的点的曲率半径为该点半径。 N=1坐标正算，N=2坐标反算；X0，Y0起点坐标；S0起点里程；F0起点切线方位角；Ls线元长度；R0起点曲率半径，RN止点曲率半径；DK所求点里程；ZJ所求点的支距；XP，YP所求点的坐标；F测点法线方位角；（F-90）所求点切线方位角。 竖曲线计算程序 【SQX】 ｛D｝：D“SJD” ：｛E｝：E“I1” ：｛F｝：F“I2” ：｛R｝：R ：｛H｝：H：T=R/2*(E-F) LbI 0 ｛K｝：K≤D＝＞Prog“AB” ：≠＞Prog“CD” Goto 0 【AB】 K≤D-T＝＞I “H0”=H+E(K-D) ◢≠＞I “H0”=H-E(D-K)-(K-D+T) ∧2/(2R) ◢ 【CD】 K≥D+T＝＞I “H0”=H+F(K-D) ◢≠＞I “H0”=H+F(K-D)-(D+T-K) ∧2/(2R) ◢ 说明：SJD竖交点；I1，I2纵坡设计坡度；R竖曲线半径；H变坡点设计高程；H0所求点高程。 缓和曲线计算公式 β=L*180/(2R∏Ls) 切线角 δ=L*180/(6R∏Ls) 弦偏角 x=L-L/(40RLs) y=L/6RLs -L/(336RLs) 缓和曲线参数方程。 P=Ls/24R 内移距 q=Ls/2 -Ls/(240R) 切线增长值 T=（R+p）tg a/2 + q 切线长 Lh=(a-2β)∏R/180 + 2Ls 曲线长 Ly=Lh-2Ls E=(R+p)sec a/2 –R 外矢距 D=2T-LQ 切曲差 直缓点里程 ZH=JD-T 缓圆点里程 HY=ZH+Ls 圆缓点里程 YH=HY+Ly 缓直点里程