程序计算在道路测量中的运用.docVIP

程序计算在道路测量中的运用 二、 曲线坐标计算： 传统计算坐标方法如下： 切线支距法： 是以直缓点ZH或缓直点HZ为坐标原点，以过原点的切线为X轴，过原点的曲率半径为Y轴，利用缓和曲线和圆曲线上各点的X、Y坐标测设曲线。 在缓和曲线上各点的坐标可按以下公式计算，即： X=L-L5/40R2Ls2 Y=L3/6RLs 圆曲线上各点坐标的计算公式如下： X=Rsinφ+q Y=R(1-cosφ)+p 式中：φ=L180o/R*π+βo； L：该点到HY或YH的曲线长。 2.2偏角法 用上述切线法算出缓和曲线和圆曲线上各点的坐标X、Y后，偏角法测设平曲线时，可将经纬仪置于ZH或HZ点进行测设。设平曲线上任意一点P的支坐标为X、Y，则其偏角δ和弦长c的计算可由三角形得： δ=arcotY/X c=√（X2+Y2） 以上两种方法为传统计算法，图示及计算实例略。 此上方法计算过程繁琐,而且高速公路坐标为国家统一坐标系,要进行坐 标转换容易产生错误,不利于施工过程的快速、准确开展,笔者在道路施工测量过程中,总结出了一套运用CASIO（fx-4800P）计算器程序计算缓和曲线上各点坐标,其特点是灵活、速度快、精度高。 .1 缓和曲线计中边桩推算思路： 如上图： 根据切线支距法公式先计算P点坐标，由P点坐标算出P-ZH距离，由起始方位角减P点的偏角δ得出ZH-P点方位角，后计算P点统一坐标；计算P点切线方位角计算左右边桩。 Xp=X0+Cos(F-δ)*√(X2+Y2); (δ=β0*L2*180/3Ls2π) Yp=Y0+Sin(F-δ)*√(X2+Y2)； （β0=Ls/2R*180/π） X边=Xp+Cos(90±F-β)*S; (β=L2/2RLs) Y边=Yp+Sin(90±F-β)*S。 式中X0、Y0，为ZH点线路统一坐标；X、Y为P点直角坐标；F为起始方位角；β为P点切线角；δ为P点偏角；L为P至ZH的缓和曲线长；S为P点左右边桩距离。 如上图： 根据圆曲线切线支距法计算公式，计算P点坐标（X，Y）由ZH点统一坐 标推算P’坐标再由P’计算P点统一坐标。 Xp’=X0+X*CosF Yp’=Y0+X*SinF; Xp=Xp’+Y*Cos(90-F) Yp=Yp’+Y*Sin(90-F); X边=Xp+Cos(90±β0-φ)*S Y边=Yp+Sin(90±β0-φ)*S。 式中F为起始方位角；S为P点左右边桩距离；φ为P点至HY点弧长所对圆心角；β0为缓和曲线切线角。 3.3 线路坐标计算及极坐标放样程序（适用于CASIO fx-4800P） 主 程 序 (QXZB) Lb1 5: Prog“ZJS”:｛LP｝: O=V:M=L-A:F=Abs(D-E)÷（D-E）:D-E(0(P=-P( T=D:U=B:Q=C:L≤A(Goto4( L≤A+V(Goto 2 ( L(G-V(Goto 1( M=G-L:T=180+E:U=H:Q=N:F=-F:LG(Goto4( Lb1 1: Z=90(M-O)÷πR:W=D-90FO÷πR-FZ:X=U+Rec(R*Sin(180(M-O) ÷2)÷πR)+O÷2-O3240R2 ,T):Y=Q+J:X=X+Rec(R-Cos(180(M-O÷2)÷πR)+O2÷ 24R,T-90F):Y=Y+J( Goto3: Lb1 2: Z=30M2÷πRO:W=T-FZ:X=U+Rec(√((M-M540R2O2)2÷(M3÷6RO)2),W):Y=Q+J:Z=2Z:S=W-FZ:S360(S=S-180:( Lb1 3: FP(0(X=X+Rec(AbsFP,270+W-FZ):Y=Y+J:( FP0(X=X+Rec(FP,90+W-FZ):Y=Y+J:( Goto 0: Lb1 4: X=U+Rec(AbsM,T+180):Y=Y+J:P(0(T=T+180( P(0(P=AbsP( X=X+Rec(P,T+90F):Y=Y+J:L≤A(S=E:LG(S=D:( Lb1 0: X ( Y ( S ( L0(Prog“JDJL”: Goto 5: 子 程 序 (ZJS) A=:G=:L≤G(B=:C=:H=:N=:D=:E=:R=:V=:Goto 1( Lb1 0 (JDJL) Pol(Z(3(-Z(1( ,Z(4(-Z(2():Z(5(=J:Z(5((0(Z(5(=Z(5(+360:Pol(X-Z(1( ,Y-Z(2():Z(6(=I: Z(6(=( J(0(J=J+360:Z(7(=J-Z(5(: Z(7(( 3.4 程序说明: 子程序(ZJS)须输入： A为ZH 点里程；G为HZ里程点