连续运算的公路测量程序(含数据库).doc

一、公路基本型非对称曲线中边桩坐标计算放样程序（连续运算）1、子程序 SUB.B （前后交点坐标反求偏角）参数意义：ZJD1X，UJD1Y，NJDX，EJDY，BJD2X，CJD2Y， OJD交点桩号，RR，SLS1，PLS2Z1B前方位角值（角度），Z2C后方位角值，Z3I判断左偏或右偏（I=-1左偏，I=1右偏），A夹角（角度单位）程序清单：Pol（NZ，EU）J0Z1=360+J：Z1=J：Pol（BN，CE）：J0Z2=360+JZ2=J A=Z2 Z1：M= ABS AM180 GOTO 1GOTO 2Lb1 1A0“I=”：Z3=-1“A=”：A= ABS ADMS GOTO 3“I=”：Z3= 1“A=”：A= ABS A DMSGOTO 3Lb1 2A0“I=”：Z3=1“A=”：A=360-Abs ADMSGOTO 3 “I=”：Z3= -1“A=”：A=360-Abs ADMSGOTO 3Lb1 32、主程序 ZBZZB2L (中边桩坐标计算主程序)参数意义：K中桩桩号，D边桩距离（恒为正），Q边桩与中桩切线的夹角（左负右正），G测站X坐标，H测站Y坐标，“N=”：X中桩X坐标，“E=”：Y中桩Y坐标，Z1B前方位角值（角度），Z2C后方位角值，Z3I判断左偏或右偏（I=-1左偏，I=1右偏），NJDX，EJDY程序清单：G“che X”+H“che Y”Lb1 1FixmKQDProg“AA”：Z21=S2(24R)S4(2384R3)：Z22=P2(24R)P4(2384R3)Z23=S2S3240R2：Z24=P2P3240R2Z25=S2R180：Z26=P2R180“T1=”：Z27=（R+Z22 （R+Z21）cosA）sinA+Z23 “T2=”：Z28=（R+Z21 （R+Z22）cosA）sinA+Z24 “LY=”：M=（AZ25Z26）R180“L=”：L=M+S+P“E=”：Z9=（R+（Z21+Z22）2）cos (A2)-RZ4=OZ27：Z5=Z4+S：Z6=Z4+S+M2：Z7=Z5+M：Z8=Z4+LKZ8 “N=”：X=N+（KZ8+Z28）cos Z2 “E=”：Y=E+（KZ8+ Z28）sin Z2 GOTO 6KZ4 “N=”：X=N（O-K）cos Z1 “E=”：Y=E（OK）sin Z1 GOTO 6KZ6 V=Z8 K：Z11=N+Z28cos Z2 ：Z12=E+ Z28sin Z2 ：Z9=Z2GOTO 2V=KZ4：Z11=NZ27cos Z1 ：Z12=EZ27sin Z1 ：Z9=Z1GOTO 3Lb1 2VP“LS2X=”：X=V-V5(40R2P2) “LS2Y=”：Y=V3(6RP)-V7(336R3P3) GOTO 4W=(V-P2)180（R）：“LYX=”：X=Rsin W+Z24 “LYY=”：Y=R（1-cos W）+Z22 GOTO 4Lb1 3VS“LS1X=”：X=V-V5(40R2S2) “LS1Y=”：Y=V3(6RS)-V7(336R3S3) GOTO 4W=(V-S2)180（R）：“LYX=”：X=Rsin W+Z23 “LYY=”：Y=R（1-cos W）+Z21 GOTO 4Lb1 4KZ6 Z14= X：GOTO 5：Z14= X：GOTO 5Lb1 5“N=”：X=Z11+Z14cos Z9- Z3Ysin Z9 “E=”：Y=Z12+Z14sin Z9+Z3Ycos Z9 Lb1 6D=0Prog“SUBRD”：GOTO 1：Prog“F”：Prog“SUBRD”：GOTO 13、子程序 SUBRD （已知两点坐标求方位角及距离即极坐标方法放样程序）参数意义：G测站X坐标，H测站Y坐标，X待放样点X坐标，Y待放样点Y坐标，Z13FWJ方位角，IZD L放样点与测站点的距离程序清单：G+H+X+YPol（X-G，Y-H）J0 “FWJ=”：Z13= 360+JDMS“FWJ=”：Z13= JDMS“ZL D=”：I4、子程序 F (切向方位角及边桩坐标计算子程序)参数意义：Z17（F）切向方位角，Q边桩与中桩切线的夹角（左负右正），D边桩距离（恒为正），NB边桩X坐标，EB边桩Y坐标，Z18（A）边桩方位角程序清单：KZ4 “F=”：Z17=Z1 DMSGOTO 5 KZ5 “F=”：Z17=Z1 +180Z3 （K-Z4）2（2RS）DMSGOTO 5 KZ7 “F=”：Z17=Z1 +180Z3 （2（K-Z4）-S）（2R）DMSGOTO 5 KZ8 “F=”：Z17=Z2-180Z3 （Z8-K）2（2RP）DMSGOTO 5“F=”：Z17=Z2 DMSLb1 5Z18= Z17+Q“NB=”：X=X+Dcos Z18“EB=”：Y=Y+Dsin Z185、数据库调用子程序 AA 程序清单：KZH（JD2处）prog“A1”：GOTO 1KZH（JD3处）prog“A2”：KZH（JD4处）prog“A3” ：GOTO 1 K终点prog“AN”：GOTO 1Lb0 16、数据库子程序 A1AN （N=1，2，3N）程序清单：程序A1Z=起点X：U=起点Y：N=JD1X：E=JD1Y：B=JD2X：C=JD2Y：O=JD1桩号：R=JD1半径：S=JD1第一缓和曲线长：P=JD1第二缓和曲线长Prog“SUB.B”程序A2Z=JD1X：U=JD1Y：N=JD2X：E=JD2Y：B=JD3X：C=JD3Y：O=JD2桩号：R=JD2半径：S=JD2第一缓和曲线长：P=JD2第二缓和曲线长Prog“SUB.B”程序ANZ=JDN-1X：U=JD N-1Y：N=JD N-X：E=JD N-Y：B=JD N+1X：C=JD N+1Y：O=JD N2桩号：R=JD N半径：S=JD N第一缓和曲线长：P=JD N第二缓和曲线长Prog“SUB.B”说明：本程序主要运行主程序 ZBZZB2，它调用其它相关的子程序来进行计算，能计算及放样出数据库内任意桩号及任意角度的中边桩坐标计算及极坐标放样。输入的基本参数为中桩桩号，边桩夹角，边桩距离K，Q，D二、公路匝道坐标计算及放样主程序（连续运算）1、主程序 QQQXZD （匝道坐标计算及放样）参数意义：A（XA）曲线元起点X坐标，B（YA）曲线元起点Y坐标，C（CA）曲线元起点切线方位角；D（+-RA）曲线元起点曲率半径（左偏负、右偏正），E（+-RB）曲线元终点曲率半径（左负右正）；F（DKA）曲线元起点桩号，G（DKB）曲线元终点桩号；M（che X）测站X坐标，N（che Y）测站Y坐标；H待求中桩桩号，X待求点中边桩的X坐标，Y待求点中边桩的Y坐标；J（FWJ）放样点至测站点的方位角，I（DL）放样点至测站点的距离；Z（BJ +0+）边桩距中线平距，恒为正；Q（JIA）中点切线与与两边桩连线夹角（左负中0右正），Z13FWJ方位角，IZL放样点与测站点的距离O切线方位角。关于曲率半径，直线段取1045；缓和曲线段ZH取1045，HY取R，YH取R，HZ取1045。并注意左偏或右偏。程序清单：M“che X”：N“che Y”Lb1 1FixmHProg“ZZ”：L=（H-F）4：I=180LD：U=90LL（D-E）DE（G-F）：T=C+I+U：W=C+2I+4U：K=C+3I+9U：O=C+4I+16U：ZQ：Z“BJ（+0+）”：P=O+Q“JIA（-0+）”：X=A+L（cos C+4（cos K+cos T）+2cos W+cos O）3+Zcos P：Y=B+L（sin C+4(sin K+sin T)+2sin W+sin O）3+Zsin P“X=”：X“Y=”：YI=0：J=0：Pol（X-M，Y-N）J0 “FWJ=”：Z13= 360+JDMS“FWJ=”：Z13= JDMS “ZL =”：IGOTO 12、 数据库调用子程序 ZZ 程序清单：K第1曲线元终点桩号prog“Z1”：GOTO 1K第2曲线元终点桩号：prog“Z2”：GOTO 1K第3曲线元终点桩号：prog“Z3”：GOTO 1 K终点prog“ZN”：GOTO 1Lb0 13、 数据库子程序 Z1ZN （N=1，2，3N）程序Z1A=第1曲线元起点X坐标：B=第1曲线元起点Y坐标：C=第1曲线元起点切线方位角：D=第1曲线元起点半径（左偏-，右偏+）：E=第1曲线元终点半径（左偏-，右偏+）：F=第1曲线元起点桩号：G=第1曲线元终点桩号程序Z2A=第2曲线元起点X坐标：B=第2曲线元起点Y坐标：C=第2曲线元起点切线方位角：D=第2曲线元起点半径（左偏-，右偏+）：E=第2曲线元终点半径（左偏