公路中、边桩坐标计算及放样程序.doc

公路中、边桩坐标计算及放样程序适用范围：单圆曲线、基本型缓和曲线、非对称基本型缓和曲线、卵型曲线、回头曲线（转角大于180度）特点：1．任意半径曲线、任意交角边桩2．子程序可独立执行源程序：F 2（文件名：2）“O=0 X Y Q SO O≠0 X Y ”：ORS“LS1”U“LS2”P“ZH-Q”W“Z-G”T“J-G”A“ZH”：B=A+S：F“HZ”：E=F-U：G“X-ZH”H“Y-ZH”I“X：HZ”J“Y：HZ”M“XO”N“YO”X“X-HS”Y“Y-HS”：Prog“A”：Prog“D”：Lbi0：{Z}：Z：D=0=>Goto1：≠>Goto2：Lbi2：{D}：D：Lbi1：Z≤A=>Prog“3”：Prog“0”：Goto0：≠>Z＞A=>Z≤B=>L=Z-A：K=RS：Prog“4”：Prog“5”：Prog“0”：Goto0：≠>Z＞B=>Z≤E=>L=S：K=RS：Prog“4”：Prog“5”：Prog“6”：Prog“0”：Goto0：≠>Z＞E=>Z≤F=>L=F-Z：K=RU：Prog“4”：Prog“9”：Prog“7”：Prog“0”：Goto0：≠>Z≥F=>C=Z-F：Prog“8”：Prog“0”：Goto0F 3C=A-Z：X=G+Ccos（P+180）：Y=H+Csin（P+180）：V=PF 4S=0=>X=0：Y=0：C=0：Q=0：V=0：≠>Goto1：Lbi1：X=L-L5÷40K2+L9÷3456K4-L13÷599040K6+L17÷175472640K8-L21÷7.80337152E10K10：Y=L3÷6K-L7÷336K3+L11÷42240K5-L15÷9676800K7+L19÷3530096640K9-L23÷1.880240947E12K11：C=√（X2+Y2）：V=90L2÷лK：X≠0=>Q=tan-1Y┛X：≠>Q=0：⊿F 5Prog“9”：X=G+Ccos（P+Q）：Y=H+Csin（P+Q）： V=P+VF 6L=Z-B：Q=90L÷лR：C=2RsinQ：Prog“9”：X=X+Ccos（AbsV+Q）：Y=Y+Csin（AbsV+Q）：V=AbsV+2QF 7X=I+Ccos（P+W-180-Q）：Y=J+Csin（P+W-180-Q）：V=P+W-VF 8X=I+Ccos（P+W）：Y=J+Csin（P+W）：V=P+WF 9W≥0=>Q：V：≠>Q= -Q：V= -V：⊿F 0X=X+Dcos（V+T）：Y=Y+Dsin（V+T）：O=0=>“X=”：X：Pause0：Y=Y◢Prog“A”：Prog“D”：≠>O≠0=>“X=”：X：Pause0：Y=Y◢V=V◢F AX＞M=>Prog“B”：Goto0：≠>X=M=>Prog“C”：≠>X＜M=>Prog“B”：Q=Q+180：Lbi0：Q=Q+360：Q≥360=>Q=Q-360◢≠>Q=Q◢F BQ=tan-1（Y-N）┛（X-M）F C Y≥N=>Q=90◢≠>Q=270◢F DC=√（（X-M）2+（Y-N）2））：“SO=”：Pause0说明：一、程序中字母及符号意义：R：圆曲线半径；LS1（S）：第一缓和曲线长LS2（U）：第二缓和曲线长ZH-Q（P）：直缓（直圆）点切线方位角Z-G（W）弯道转角（左转为负，右转为正）J-G（T）：中桩至右侧某点方向与中桩切线方位角的夹角（大于等于0度且小于等于180度，当正交时为90度）ZH（A）：直缓或直圆点桩号HZ（F）：缓直或圆直点桩号X-ZH（G）：直缓或直圆点X坐标；Y-ZH（H）：直缓或直圆点Y坐标X-HZ（I）：缓直或圆直点X坐标Y-HZ（J）：缓直或圆直点Y坐标XO（M）：测站X坐标； YO（N）：测站X坐标 X-HS：后视点X坐标； Y-HS：后视点Y坐标X、Y：计算或放样点坐标Q（Q）：计算或放样方位角 SO（C）：计算或放样距离Z（Z）：计算点桩号D（D）：边桩距中桩宽度（左为负值，右为正值）V（V）：中桩切线方位角二、输入、计算要点1．该程序一次只可输入一个弯道的参数，计算段落为上一弯道终点（HZ或YZ）至下一弯道起点（ZH或ZY）。

2．计算单圆曲线时LS1、LS2输入时输03．当只计算第一缓和曲线及圆曲线，不计算第二缓和曲线时，弯道转角只需输入正或负值（左转为正、右转为负）即可，可不输入准确的角度4．如只计算上一弯道曲线（HZ或YZ）终点至ZH点间的桩号时，除ZH点桩号、方位角和坐标须输入正确外，其它参数可任意输入（但不能有逻辑错误，如HZ点桩号小于了ZH点的桩号）5．后视点X、Y坐标在每次重启程序时需重新输入，建议第一次输入计算出后视角度及距离后用笔将计算的后视角度和距离在纸上记录下来，便于中途后视时核对6．用于卵型曲线时坐标计算时，需先计算出卵型曲线的相关参数后再用本程序计算即可7．为保证计算精度，本程序所用切线支距公式（F4子程序中）取至6项，计算非卵型曲线时为节省存储空间可只取前4项，建议至少取3项8．程序启动时，当问O？时，如输入0则计算结果显示坐标X、Y及放样角度Q和放样距离SO；如输入E3（1000）时显示坐标X、Y及中桩切线方位角9．当问距中桩宽度D？时，如只放中桩可输入0，程序默认为中桩后以后不再提示输入宽度；如要同时计算中、边桩，在计算中桩时可输入0.0001或E- 4等更小的数即可10．如已知两点坐标只计算其方位角或距离时，直接运行F A或F D子程序即可。

11．如需将多个曲线参数一次置入程序，减少每次输入的繁锁，可对F 2程序进行修改即可，其它子程序不用改动公路路线坐标反算适用范围：单圆曲线、基本型缓和曲线、非对称基本型缓和曲线、卵型曲线、回头曲线（转角大于180度）特点：1．已知任意点坐标求对应于中桩任意方向其对应中桩的桩号及到中桩距离，斜交或正交2．无需人为判断即可求得其位于路线的关系3．本程序公用了前面《公路中、边桩坐标计算及放样程序》中F 3～F9子程序,因此节约了空间源程序一与源程序二任取一个均可源程序一：F E ZUO BIAO FAN SUANZ=0：O=0：D=0：RS“LS1”U“LS2”P“ZH-Q”W“Z-G”T“J-G”A“ZH”：B=A+S：F“HZ”：E=F-U：G“X-ZH”H“Y-ZH”I“X：HZ”J“Y：HZ”M“XO”N“YO”：{D}：D：Lbi0：Prog“F”： Goto0F FZ=Z+O÷2：Lbi1：Z≤A=> Prog“3”：Prog“G”： Goto1：≠>Z＞A=>Z≤B=> L=Z-A：K=RS：Prog“4”：Prog“5”：Prog“G”：≠>Z＞B=>Z≤E=> L=S：K=RS：Prog“4”：Prog“5”：Prog“6”：Prog“G”：≠>Z＞E=>Z≤F=> L=F-Z：K=RU：Prog“4”：Prog“9”：Prog“7”：Prog“G”：≠>Z＞F=> C=Z-F：Prog“8”：Prog“G” F G C=√（（X-M）2+（Y-N）2））：D＜0=>D= -C：Goto0：≠>D≥0=>D=C：Goto0：Lbi0：X=X+Dcos（V+T）：“X=”：Pause0：Y=Y+Dsin（V+T）：“Y=”：Pause0：O=√（（X-M）2+（Y-N）2））：O＜E9=> D=D：Goto1：≠>O≥E9=>Z=0：O=0：D= -D：Goto1：Lbi1： Abs（X-M）≤E-4=> Abs（Y-N）≤E-4=> X ：“X=”：Pause0： Y ：“Y=” ◢ Z ：“Z=”： Pause0： D ：“D=” ◢{MND}：M“XO”N“YO”D“D”：Z=0： O=0：≠>Abs（X-M）＞E-4=> Abs（Y-N）＞E-4=> Z源程序二：F E ZUO BIAO FAN SUANZ=0：O=0：D=0：RS“LS1”U“LS2”P“ZH-Q”W“Z-G”T“J-G”A“ZH”：B=A+S：F“HZ”：E=F-U：G“X-ZH”H“Y-ZH”I“X：HZ”J“Y：HZ”M“XO”N“YO”：{D}：D：Lbi0：Z=Z+O÷2：Z≤A=> Prog“3”：Prog“F”： Goto0：≠>Z＞A=>Z≤B=> L=Z-A：K=RS：Prog“4”：Prog“5”：Prog“F”：Goto0：≠>Z＞B=>Z≤E=> L=S：K=RS：Prog“4”：Prog“5”：Prog“6”：Prog“F”：Goto0：≠>Z＞E=>Z≤F=> L=F-Z：K=RU：Prog“4”：Prog“9”：Prog“7”：Prog“F”：Goto0：≠>Z＞F=> C=Z-F：Prog“8”：Prog“F” ：Goto0F F AbsX≥E9=>X=0：Y=0：Z=0：O=0：D= -D：≠> Goto0：Lbi0： C=√（（X-M）2+（Y-N）2））：D＜0=>D= -C：Goto1：≠> D=C： Lbi1：X=X+Dcos（V+T）：“X=”：Pause0：Y=Y+Dsin（V+T）：“Y=”：Pause0：O=√（（X-M）2+（Y-N）2））：Abs（X-M）≤E-4=> Abs（Y-N）≤E-4=> X ：“X=”：Pause0： Y ：“Y=” ◢ Z ：“Z=”： Pause0： D ：“D=” ◢{MND}：M“XO”N“YO”D“D”：Z=0说明：一、程序中字母及符号意义：D：待求点距中桩宽度（左负右正）XO（M）：待求点已知X坐标YO（N）：待求点已知Y坐标Z：待求点对应中桩桩号其它与公路中、边桩坐标计算及放样程序中符号意义相同。

二、输入、计算要点1．D值的输入：如能判断其位于路线的左或右侧（其值为左负右正），则会缩短程序运行时间；也可不输入由程序自行计算2．参数的输入：如能确定所求点坐标所处的曲线位置时则输入该曲线参数，也会缩短程序运行时间3．假如事先不知道已知点位于哪一曲线上又怎样处理呢？输入你认为已知点可能位于的那个曲线参数，然后看计算结果的桩号是否处于你输入的那个曲线，如在则结果可用，如不在则输入计算结果显示桩号所处的那个曲线参数即可4．程序计算结果精度设定小于0.1mm，为缩短运行时间可将精度重设，方法是在子程序F G中将有E-4的改为E-3（精度设定小于1mm）5．当已知点所处的曲线半径较小时，如R为20、30m时，程序可能无法得出正确结果或出现死循环时，将程序中止进入程序F F中将表达式Z=Z+O÷2改为Z=Z+O÷4即可，改后的程序执行所需时间较修改前长，建议一般不改，除非出现上述情况时再作修改6．此程序必须与《公路中、边桩坐标计算及放样程序》中F 3～F 9子程序结合使用，不能独立运行已知一点A求其对应于中桩的桩号及到中桩的宽席计算原理（试算渐近法）：此分析是基于用于编程而考虑的1. 先假定从路线起点K0+000处为计算起点(如果不是编程序可不从0桩号开始，可根据经验或实际进行估计)，求得0桩号的中桩坐标，再求得其与已知点的距离C1，利用中桩算得C1宽度的边桩坐标点B1，再算B1点与已知点A的距离D1。

2. 算得Z1的桩号（为Z0的桩号加上前一个桩号对应边桩。