曲线任意里程中边桩坐标正反算文档格式.docx

1、SUB1XS=:X=XYS=Y=YGoto 3Lbl 2：XY：XY：I=X：J=Y：SUB2S=S=O+WZ=Z=ZGoto 2Lbl 3:J=0:I=Pol（X-Z1,Y-Z2）:F=J:FF=F+360:F360=F=F-36”DIST=”:I”FW=”:FDMS Goto 12. 正算子程序（SUB1）A=0.1184634425:B=0.2393143352:N=0.2844444444:K=0.0469100770:L=0.2307653449:M=0.5:F=1-L:M=1-K:X=U+W（Acos（G+QEKW（C+KWD）+Bcos（G+QELW（C+LWD）+Ncos（G+

2、QEMW（C+MWD）+Bcos（G+QE（1-L）W（C+（1-L）WD）+Acos（G+QE（1-K）W（C+（1-K）WD）: Y=V+W（Asin（G+QEKW（C+KWD）+Bsin（G+QELW（C+LWD）+Nsin（G+QEMW（C+MWD）+Bsin （G+QE（1-L）W（C+（1-L）WD）+Asin （G+QE（1-K）W（C+（1-K）WD）：F=G+QEW（C+WD）+T：X=X+ZcosF：Y=Y+ZsinF3. 反算子程序（SUB2）T=G-90：W=Abs（Y-V）cosT-（X-U）sinT）：Z=0：Lbl 0：L=T+QEW（C+WD）：Z=（J-Y）c

3、osL-（I-X）sinL：AbsZGoto1：W=W+Z：Goto 0Z=（J-Y）sinF4.曲线元要素数据库：DAT1J=1=U=xxxx.xxx:V=xxxx.xxx:O=xxxx.xxx:G=xxxx.xxx:H=xxxx.xxx:P=xxxx.xxx:R=xxxx.xxx:Q=xxJ=2=Q=xx.J=n=U. （注：如有多个曲线元要素继续添加入数据库DAT1中）三、使用说明1、规定（1） 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。（2） 当所求点位于中线时，Z=0；当位于中线左铡时，Z取负值；当位于

4、中线中线右侧时，Z取正值。（3） 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。（4） 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。 （5） 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替； （6） 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；止点与直线相接时，曲率半径等于设计规定的值；（7）曲线元要素数据库（DAT1）可根据线型不同分为各个线元段输入到DAT1中，即分为直线段、缓和曲

5、线、圆曲线等。2、输入与显示说明 输入部分：1. SZ = XY2. XY = SZ N ?选择计算方式，输入1表示进行由里程、边距计算坐标 ；输入2表示由坐标反算里程和边距。 X0 ？线元起点的X坐标 Y0 ？线元起点的Y坐标 S0 ？线元起点里程 F0 ？线元起点切线方位角 LS ？线元长度 R0 ？线元起点曲率半径 RN ？线元止点曲率半径 Q ？线 元左右偏标志（左偏Q=-1，右偏Q=1，直线段Q=0）S ？ 正算时所求点的里程Z ？正算时所求点距中线的边距（左侧取负，值右侧取正值，在中线上取零）ANG？正算边桩时左右边桩连线与线路中线的右交角J？曲线元数据库曲线段判断系数（J=1、2

6、.n）X ？反算时所求点的X坐标Y ？反算时所求点的Y坐标 显示部分： XS=正算时，计算得出的所求点的X坐标 YS=正算时，计算得出的所求点的Y坐标 S=反算时，计算得出的所求点的里程 Z=反算时，计算得出的所求点的边距四、算例某匝道的由五段线元（直线+完整缓和曲线+圆曲线+非完整缓和曲线+直线）组成，各段线元的要素（起点里程S0、起点坐标X0 Y0、起点切线方位角F0、线元长度LS、起点曲率半径R0、止点曲率半径RN、线 元左右偏标志Q）如下：S0 X0Y0 F0 LSR0 RN Q500.00019942.83728343.561 125 16 31.00269.256 1E45 076

7、9.25619787.34028563.37837.492 221.75-1806.74819766.56628594.574 120 25 54.07 112.779221.75 -1919.52719736.07228701.89391 17 30.6380.2859579.228999.81219744.03828781.65980 40 50.00 100.0001E450 1、正算（注意：略去计算方式及线元要素输入，请自行根据所求点所在的线元输入线元要素）S=700Z=-5计算得XS=19831.41785 YS=28509.72590Z=0 计算得XS=19827.33592 YS

8、=28506.83837Z= 5XS=19823.25398 YS=28503.95084S=780XS=19785.25749 YS=28575.02270XS=19781.15561 YS=28572.16358XS=19777.05373 YS=28569.30446 S=870XS=19747.53609 YS=28654.13091XS=19742.68648 YS=28652.91379XS=19737.83688 YS=28651.69668S=940Z=-5.123XS=19741. 59118YS=28722.05802XS=19736.47687 YS=28722.3564

9、2Z= 3.009XS=19733.47298 YS=28722.53168 2、 反算 X=19831.418 Y=28509.726S=699.9999974Z= -5 .00018164 X=19827.336 Y=28506.838S=699.9996493Z= 0.000145136 X=19823.25398 Y=28503.95084S=699.9999985Z= 5.000003137 X=19785.25749 Y=28575.02270S=780.0000035Z= -5 .000001663 X=19781.15561 Y=28572.16358S=780.0000025

10、Z=- 0.000002979 X=19777.05373 Y=28569.30446S=780.0000016Z= 4.99999578X=19747.536 Y=28654.131S=870.0001137Z= -4.99941049 X=19742.686 Y=28652.914S=870.0003175Z=- 0.00041814 X=19737.837 Y=28651.697S=870.0002748Z= 4.999808656X=19741.5912 Y=28722.0580S=939.9999786Z= -5.123024937X=19736.4769 Y=28722.3564S

11、=939.9999862Z=- 0.000027710X=19733.4730 Y=28722.5317S=940.0000238Z= 3.00898694利用Simpson公式和垂距趋近正反算曲线程序利用Simpson公式和垂距趋近正反算曲线程序 一、程序编制原理和功能：1、原理：利用Simpson公式和点至主点法线垂距逐步迭代趋近单向收敛原理2、功能：a、正算任意角度、偏距点坐标；b、反算任意角度的中桩里程和偏距。二、程序主体：1、主程序：ZFST1.KDF=XY,2.XY=KDAQDBZD1/R01/RNCDET=1=Goto 1:Goto 2Lbl 1:LIF:LK=IDP=FANG=

12、SIM2XI=X=X:Pause 0:YI=Y=YFW=WDMSGoto 1Lbl 2:XY:XY:O=X:U=Y:L=A+Abs （Y-D）cos（E+F-180）-（X-C）sin（E+F-180）I=0:Z=（U-Y）cos（W+F-180）-（O-X）sin（W+F-180）:1E-4=Goto 4:L=L+Z:Goto 3Lbl 4:Z=（U-Y）/sin（W+F）Lbl 5:L:ZGoto 22、子程序：SIM2N=6：T=（VR）（BA）：H=（LA）P=TNH+R：S=9HW=E+（P+R）NSG=1： X =C+H6（cosE+cosW+4（cos（E+（G+.5）HT+2R

13、）（G+.5）S），G，（N1）+2（cos（E+（THG+2R）GS），G，1，（N1）+Icos（W+F）：Y =D+H（sinE+sinW+4（sin（E+（G+.5）HT+2R）（G+.5）S，G，（N1）+2（sin（E+（THG+2R）GS，G，1,（N1）+Isin（W+F）三、使用说明： QD起点里程； ZD终点里程； 1/R、1/RN起、终点曲率（左负右正） X、Y起点坐标； F起点方位角；N可取412，视线型而定，亦可改为变量； ANG线外左右测点连线与路线方向斜交右角； X、Y待求点中桩坐标 XI、YI计算点中边桩坐标。 FW切线方位角 K= 正算时为计算点里程；反算时为

14、测点对应中桩里程 DP正算时为支距（左负右正）；反算时为测点对应于中线某一角度ANG的支距（若ANG=90则为垂支距，左负右正）。利用真数（级数展开）公式和垂距趋近正反算曲线程序（for 4850P）利用真数（级数展开）公式和点至法线垂距逐步迭代趋近单向收敛原理；功能：1、正算任意角度、偏距点坐标；2、反算任意角度的中桩里程和偏距。ZFS-ZHSHQ-1,+1H=L-A:T=（V-R） / 2 /（B-A）:W=E+Q（RH+TH2）180/P=H-R2H3/6-TRH4/4+（R4-12T2）H5/120+TR3H6/36+T2R2H7/28+T3RH8/48+T4H9/216S=RH2/2

15、+TH3/3-R3H4/24-TR2H5/10+（R5-60T2R）H6/720+（TR4-4T3）H7/168+T2R3H8/96+T3R2H9/108+T4RH10/240+T5H11/1320X=C+PcosE-QS sinE+Icos （W+F）Y=D+PsinE+QS cosE+Isin （W+F） 1/R、1/RN起、终点曲率（直线输入零）Q曲线偏向（左-1、右+1）；线路中边桩任意里程坐标正反算改编程序CASIO4850 含数据库、测站放样计算（第五次修改-正式版）1.加编数据库及计算总调度程序,计算中不必逐项输入线元要素,提高运算速度，避免现场忙中出错2.可加入多条线路的数据库

16、,内业输入数据,外业一目了然3.计算直观,人性化4.正算直接输入里程和边距,反算输入近似里程便可5.增加了“计算点与测站点”的距离和方位角计算语句，方便直接放样6.愿收获与大家共享,同时也希望大家提出心得和宝贵建议7.核心计算程序摘自“yshf”一.改动后的程序清单，增设数据库程序,可输入n条线路的数据库1.（QXZDJS 计算总调度程序）M=1=SJK1SUBSJK.M=n=SJKn2.（SJK1 数据库程序）N,S：SDKIS本线元终点里程U=本线元起点X坐标：V=本线元起点Y坐标：G=本线元起点正切线方位角：P=本线元起点曲率半径：R=本线元终点曲率半径：O=本线元起点里程：H=本线元终

17、点里程：Q=本线元左直右偏向（0或1、-1）:Goto0S下一线元终点里程=O=。 Goto0LB1 03.（SUBSJK 运算主程序）UVGPROHQ：XAYAFARARBDKADKBQDeg:SH=Goto AD=（P-R）2PR（H-O）：Goto 2Z：Z：XYQX-FWJ=F=F-90DMSI=0：J=0：Pol（X-CXC，Y-EYC）：J 二. 使用说明：1.把所有相关的”线元要素“依次输入”SJK“，如果有多条线路，可把“SJK”分成SJK1.SJKn等n条线路数据库，然后对应“计算总调度程序”中的M1.Mn备用2.运算时直接调用“SJK”运行，程序提示输入里程“S”？时，正算直接输入待求点里程，反算输入所求点“近似“里程3.如果在正算中途想进行反算，重新给”N“？赋值2便可4.程序中“XC、YC”为测站坐标，“I”、“FWJ”为放样距离和方位角5.程序中“QX-FWJ=”为计算点在中线上的正切线方位角，此方位角在新线路输入数据库时非常重要，我们在每输入完一段曲线元要素后，立即退出返回到“计算总调度程序（SUBSJK）”中，然后输入该段曲线元的终点里程，其计算出来的X，Y中桩坐标应该与url=javascr