坐标计算4800源程序.doc

公路通用复化辛普森公式匝道点位坐标计算4800源程序 公路通用复化辛普森公式匝道点位坐标计算4800源程序-杭浦高速临平互通- 本文利用的是计算公路匝道点位坐标的复化辛普森通用公式数学模型,集直线、圆曲线、回旋线通用，占字符内存较小，计算精度不限的程序一、运行变量名称说明： V=1、2分别进入坐标计算、桩号反算 K1、K2-曲线起点、终点里程 F0-曲线起点方位角 R1、 R2-曲线起点、终点半径 (左-右+，0为直线) X0、Y0-曲线起点、终点坐标 M-求和累积次数n的2倍(偶数)，精度迭代次数 K-曲线待求点里程 BP-求点左右偏距(左-右+) ANG- -求点的右斜交角 X、Y-曲线求得坐标 FW-待求点的即时切线方位角 XF、YF-为需求桩号的点坐标 DL、K+O、LP分别为桩号误差、求得桩号、左右偏距(左-右+)当曲线的设计半径较小时，为保证点位计算精度，M(即程序中n的2倍)的取值可适当的大些。M为偶数，直线时M=2即可，经计算M=16即可满足半径为60的小半径曲线精度。二、曲线计算程序名: Prog CURVEDefm 4VV=1 2Lbl 0:KLWLbl 4QOPT:M0AB1C2D3E4FH5G6I7J8CR9=0=Prog MQ=1=Prog ABQ=2=Prog CQ=3=Prog DQ=4=Prog EQ=5=Prog FHQ=6=Prog GQ=7=Prog IQ=8=Prog JQ=9=Prog CRAK1BK2CF0DR1ER2FX0GY0D0=I=1/D:=I=DE0=J=1/E:=J=EAbsD+AbsE=0=M=2:=M=16V=2=L=0:W=90KLBPWANGN=0：Z1=0：Z2=0：Z3=0：Z4=0 进入坐标迭代计算Lbl2N=N+1：H=2(K-A)/M：R=NH/2+A：R=C+180/*(I+(J-I)/2(B-A)*(R-A)*(R-A)Int(N/2)=N/2=Z1=Z1+cosR：Z2=Z2+sinR：=Z3=Z3+cosR：Z4=Z4+ sinRN=M=Goto3：=Goto 2Lbl3X=F+H/6*(cosC+4Z3+2Z1-cosR)+Lcos(R+W)Y=G+H/6*(sinC+4Z4+2Z2-sinR)+Lsin(R+W)V=2=Goto 6XX=YY= RFW=R-360Intg(R/360Goto 0 Lbl 6 进入桩号求算Pol(TXF-X,UYF-YO=Icos(J-RAbsO1e-4=ODLK=K+O OLP=Isin(J-RTUGoto 6:=K=K+O:L=0:Goto 4三、数据文件： 线元要素数据文件每行为一个线元段，逐句执行赋值，直至不满足、运行完成。主线M:K=5865=A=5865:B=586.526:C=1014958:D=0:E=0:F=3361488.152:G=449202.642KB=A=B:B=6235.526:C=1014958:D=0:E=-3250:F=3361488.044:G=449203.157KB=A=B:B=7674.84:C=983416.9:D=-3250:E=-3250:F=3361419.068:G=449566.617KB=A=B:B=8044.84:C=731149.3:D=-3250:E=0:F=3361521.529:G=450990.516KB=A=B:B=8273.243:C=69567.9:D=0:E=0:F=3361641.834:G=451340.3558+273=8+100 S=1 2:SS=2=A=8100:B=9610:C=C:D=0:E=0:F=3361720.194:G=451554.895注：主线中S=1、2为断链1、2段(如K8+273.243=K8+100)匝ABK0:A= 0:B=716.396:C=854615:D=800:E=800:F=3361721.181:G=450482.643KB=A=B:B=776.396:C=137444:D=0:E=-60:F=3361468.151:G=451127.474KB=A=B:B=842.929:C=1082551.6:D=E:E=-60:F=3361431.988:G=451174.516KB=A=B:B=902.929:C=445350:D=E:E=0:F=3361446.56:G=451235.989KB=A=B:B=947.929:C=161458:D=E:E=180:F=3361499.992:G=451261.793KB=A=B:B=1047.145:C=232441.1:D=E:E=180:F=3361542.602:G=451276.164KB=A=B:B=1141.034:C=545934.26:D=E:E=0:F=3361618.517:G=451338.083匝CK842.929:A= 842.929:B=932.843:C=445350:D=-60:E=-60:F=3361446.56:G=451235.989KB=A=B:B=1067.843:C=3190208.44:D=E:E=0:F=3361528.246:G=451238.793KB=A=B:B=1204.663:C=2543440:D=E:E=0:F=3361541.195:G=451111.894匝DK0:A= 0:B=46.596:C=115737:D=791.747:E=60:F=3361640.231:G=450877.857KB=A=B:B=101.697:C=139339.2:D=E:E=60:F=3361614.603:G=450916.347KB=A=B:B=234.195:C=1914042.5:D=E:E=3237.41:F=3361563.143:G=450929.781匝EK0:A= 0:B=53.571:C=1153028:D=0:E=-420:F=3361652.781:G=450889.678KB=A=B:B=303.928:C=1115113.6:D=E:E=-420:F=3361630.748:G=450938.498KB=A=B:B=385.416:C=77421.6:D=E:E=0:F=3361610.206:G=451184.308KB=A=B:B=516.85:C=72831.93:D=E:E=0:F=3361632.665:G=451262.605匝FHK0:A= 0:B=687.778:C=2654615:D=0:E=0:F=3361721.183:G=450482.658KB=A=B:B=747.778:C=2654615:D=E:E=-60:F=3361670.463:G=449796.752KB=A=B:B=956.691:C=237722.6:D=E:E=-60:F=3361656.351:G=449739.118KB=A=B:B=1038.357:C=373732.61:D=E:E=0:F=3361569.331:G=449819.209KB=A=B:B=1242.528:C=3583759:D=E:E=0:F=3361647.7:G=449835.267匝GK838.838:A= 838.838:B=898.838:C=150102:D=-60:E=0:F=3361576.113:G=449719.64KB=A=B:B=949.463:C=121319.6:D=E:E=160:F=3361537.147:G=449764.388KB=A=B:B=1033.468:C=130351.34:D=E:E=160:F=3361508.476:G=449806.043KB=A=B:B=1121.958:C=1603956.8:D=E:E=986.28538:F=3361439.935:G=449852.931匝IK0:A= 0:B=66.176:C=2654615:D=0:E=85:F=3361694.914:G=450016.224KB=A=B:B=105.294:C=288427:D=E:E=85:F=3361698.579:G=449950.594KB=A=B:B=157.408:C=3142632.5:D=E:E=400:F=3361718.698:G=449917.45KB=A=B:B=256.758:C=3354420:D=E:E=400:F=3361762.343:G=449889.539KB=A=B:B=337.758:C=3495811:D=E:E=0:F=3361857.034:G=449860.326匝JK0:A= 0:B=66.176:C=2654615:D=0:E=-85:F=3361678.559:G=450017.433KB=A=B:B=104.77:C=243282.8:D=E:E=-85:F=3361665.281:G=449953.056KB=A=B:B=161.779:C=217279:D=E:E=-350:F=3361640.922:G=449923.549KB=A=B:B=215.557:C=1933420.8:D=E:E=-350:F=3361588.865:G=449901.365KB=A=B:B=279.843:C=184468:D=E:E=0:F=3361535.825:G=449892.803KB=A=B:B=387.002:C=1793025:D=E:E=0:F=3361471.579:G=449891.389匝CRK3080:A= 3080:B= 3666.041:C=1754615:D=0:E=0:F=3362052.282:G=449829.361KB=A=B:B= 3786.041:C=1754615:D=E:E=900:F=3361467.83675:G=449872.57911KB=A=B:B=3920:C=1793525.95:D=E:E=900:F=3361348.02009:G=449878.76613其他说明：数字前的乘号不能省。如A5、6不当，应改为A*5、*6或5A、6。括号前的乘号可省,三角函数前的乘号可省。如H(A+B)、LcosR正确以下形式为连续除法：1/AB=1/A/B,1/5(B+C)=1/5/(B+C),1/A2=1/A/A。句末尾的括号可省。如Pol(TXF-X,UYF-Y。=前不用冒号=前有冒号前后不用冒号，且可不换行。李甲生 交点法坐标计算fx-4800程序平曲线坐标、反算桩号计算程序prog PQXmn:MX0:NY0:Defm 1:Z1=0Lbl 1:ABGHVWK:KJL:GJX:HJY:VI0:WJ0:AR:BLS:FixmP=B2/24/A:Q=B/2-B3/240/A2:T=(A+P)tan(Abs W/2)+Q:C=K-T:D=C+B:F=D+*A*AbsW/180:E=F-Bws=-1:=s=1 Lbl 2:L:L:L=0=GOTO 1Lbl 6:LC=O=K-L:R=G-OcosV:U=H-OsinV:Z=V:GOTO 3LF=O=L-F+T:Z=V+W:R=G+OcosZ:U=H+OsinZ:GOTO 3LO=L-C:prog XY:Z=V+SR:R=G+(I-T)cosV-SJsinV:U=H+(I-T)sinV+JScosV:GOTO 3LE=O=F-L:prog XY:O=V+W:Z=O-SR:R=G+(T-I)cosO-SJsinO:U=H+(T-I)sinO+SJcosO:GOTO 3O=180(L-D+B/2)/A:I=AsinO+Q:J=P+A(1-cosO):Z=V+SO:R=G-(T-I)cosV-SJsinV:U=H-(T-I)sinV+SJcosVLbl 3:Z1=0.002=GOTO 5:=O=0RXLUYLZZLLbl 4:O:OYC:O=0=GOTO 2O=-1=GOTO 5X=R-OsinZ:Y=U+OcosZ:pol(X-M,Y-N:JJ=J+360XXZYYZJA0IDDGOTO 4Lbl 5:XY:XXF:X=-1=Z1=0:GOTO 4YYF:pol(X-R,Y-U:O=Icos(J-Z):AbsO0.001=ODLL+OO=Isin(J-Z):OYCGOTO 5:=Z1=0.002:L=L+O:GOTO 6prog XYU=AB:I=O-O5/40/U2:J=O3/6/U-O7/336/U3:R=90O2/U变量名称说明：MX0:NY0为测站点坐标 KJL:GJX:HJY:VI0:WJ0:AR:BLS分别为交点桩号、坐标、起始方位角、转角(左-右+)、半径、缓长L、RXL、UYL、ZZL 分别为输入桩号、求出中桩坐标、即时方位角 OYC：XXZYYZJA0IDD 分别为输入边桩距离(左-右+)、求出边桩坐标、边桩点到测站点的方位角和平距OYC=0返回计算中桩坐标OYC=-1进入计算桩号过程(输入坐标XXF 、YYF分别为求出桩号L+O,及至中桩边距OYC，ODL为桩号计算的误差距)OYC=约0.00001可求中桩点到测站点方位角和平距 XXF=-1返回计算坐标过程 计算桩号过程时，需多次确认输入的坐标XXF、YYF直到出现ODL桩号误差距较小时即可得出准确的桩号=前不用冒号=前有冒号前后不用冒号fx-4800公路超高及加宽计算作者：yangchy1文章来源：本站原创点击数： 1762更新时间：2007-5-21一、平曲线加宽、超高程序符号A路肩宽度 B路面宽度 I0路肩坡度 I1路拱坡度IB超高横坡 IX路线纵坡，上坡为正，下坡为负 BJ路面加宽值BJXX距离处路面加宽值 W未加宽前路面宽度 WS待求定的路面加宽值M加宽缓和段起点桩号 N加宽缓和段终点桩号 HN路基内缘与设计高之高差HW路基外缘与设计高之高差超高及加宽计算.功能及适用范围.绕中轴旋转的超高及加宽计算。.绕中轴旋转的超高及加宽计算。.有关规定.当线形设计须采用较长的回旋线时，横坡坡度由2%（或1.5%）过渡到0%路段的超高渐变率不小于1/330。超高过渡应在回旋线全长范围内进行。当超高渐变率太小时，超高的过渡可设在回旋线的某一区段范围之内。.路肩横坡：A、直线路段的路肩横坡可与行车道横坡相同。B、曲线段的路肩横坡：当硬路肩宽度2.25米时，曲线外侧路肩横坡的方向及其坡度值见表1。 路肩横坡方向及其坡度值 表1行车道超高值（%）2、3、4、56、78、9、10曲线外侧路肩横坡方向向外侧倾斜向外侧倾斜向内侧倾斜曲线外侧路肩坡度（%）-2-1与行车道横坡相同 当硬路肩宽度2.25米时，曲线外侧路肩横坡的方向及其坡度值与行车道相同。 路肩横坡方向与坡度变化处应设过渡段，过渡段的渐变率规定见表2。 路肩横坡过渡的渐变率 表2与行车道横坡方向相同时1/150与行车道横坡方向相反时1/100 路肩的横坡改变倾斜方向的旋转轴为路缘带外侧边缘。.程序清单：.中轴旋转ZHONGK“I1”Z=0Goto 1：Goto 2Lb1 1I“I0”H“HP”=-（AI/100+BK/200）Goto 9Lb1 2C“ZH”：J“IB”：S“LS”X=Abs（Q-C）F=0Goto 3：Goto 4Lb1 3D=-K/100+（K+J）*X/100*XH“HW”=D（A+B/2）Goto 9Lb1 4W“BJ”L=2KS/100（K+J）/100）M=0Goto 5：Goto 6Lb1 5G“BJX”=XW/S XLGoto 7：Goto 8Lb1 6P=X/SG“BJX”=（4P3-3P4）W XLGoto 7：Goto 8Lb1 7H“HN”=-（A+B/2+G）K/100 Goto 9Lb1 8H“HN”