4500程序.docVIP

4500程序 1：以知线外任意点坐标，求对应线路里程 在缓和曲线上，要计算任意里程的法线方向及任意宽度的边线坐标，非常简单。但要计算任意一个已知坐标点，是对应哪一个里程法线方向上的点，就有一些困难。很难推导一个这样的计算公式。唯一的方法“渐进”，如果手工计算这可不是一个好方法。但在有CASIO系列可编程计算器，如：FX-4500的情况下就变的非常简单了。亦可用于直线和圆曲线的计算。 首先在缓和曲线上任选一点A为起始点，计算该点的坐标和切线方位角，通过坐标反算求起始点A与计算点B的方位角和距离，B点肯定对应A点切线方向上有一个垂足C点，把三点看成一个直角三角形，通过解直角三角形计算AC的距离，当该距离大于某一数值，如0。001m，A点里程加AC的距离等于C点的里程，回到开始重新进入新一轮的计算，如果AC的距离小于某一规定值，则计算C点的里程与BC的距离即可。 求对应线路里程程序： 主程序QLC (已知坐标求里程) Lb1 0：｛LＤE｝：Prog XH：Goto 0 子程序：XH (循环) L1 Lb1 1 L2 Norm： Prog LYYD： L3 PO1(D-X,E-Y)：W≤0= W=W+360⊿ L4 Z=W-I： A=V×cos Z：L=L+A L5 Abs A≥0.001=Goto 1：≠=B=V×sinZ：Fix 3：“FXJL=” L6 L:Fix3:“DYLC=” ◢ 程序中字母代表 D 任意点X坐标 ， E 任意点Y坐标，DYLC 对应里程， FXJL 中线法线距离。程序中有坐标反算功能。 使用方法：只需输入计算点坐标、和较为接近的桩号。桩号越接近计算速度越快 2：逐桩坐标计算 2.1编制方法：线路坐标程序是按照平曲线为单元，直线部分归属在曲线两端的方法，把整段路线分段装进数据库，根据桩号判断采用数据通过共用程序，进行任意点的坐标计算，在坐标转换示意土，第一直线段，是通过方位角和距离直接计算大地坐标，第一缓和曲线和圆曲线段，是先计算任意点切线支距和方位角然后转换大地坐标，第二缓和曲线段和直线段是先计算任意点切线支距和方位角。然后转换为ZH坐标系的坐标，通过ZH坐标系的坐标再转换为大地坐标。 2.2使用方法 2.2.1准备工作：室内把已知曲线条件，装进数据库，曲线划分界线、判断条件装进子程序LYYD (路由引导)。 2.2.2现场使用：根据计算机提示输入相关数据即可。提示情况如下： K 公里桩号如 312，启动程序出现一次。 L 细部里程桩号如 518.如采用渐进只出现一次，否则逐桩输入。过千米桩时需输入1000确认。 O 渐进长度，如20米一点,取O =20，公里桩号也自动渐进。否则O=0，启动程序出现一次。 Y 断链条件，执行输0不执行默认 ，不输入【 】以内的程序，Y不出现。 E 边线角度，法线为90度，分正负值，输E=0此后则不在出现计算边线的过程。 D 边线点至中线点的距离 V W 输出的边线1的大地坐标 X Y 输出的边线2和中线的大地坐标, 2.3逐桩坐标计算程序 主程序：XLZB(线路坐标) L1 Lbl 0：L≥1000=P=P+1：L=L-1000 L2 O=0 = prog FJJ⊿ L=L+O： prog LYYD： progXSZB： E≠0 = progBX Goto 0 子程序：FJJ (非渐进) {L}：L=L 子程序：LYYD (路由引导) N=（P“K”+L/1000）×1000：【X=0：Y=6】 N≥***.*** =prog 1： prog PQX： prog ZJ N≤***.*** =prog 2： prog PQX： prog YJ N≥***.*** =prog 3： prog PQX： prog ZJ N≥***.*** =prog 4： prog PQX： prog YJ N≤***.***= prog 5： prog PQX： prog YJ N≥*** …… ……progZB 子程序：PQX（单圆曲线及带缓和曲线的平曲线） L1 B=Z+Q-S：H=Z+Q：A=S/2-S3/240 R2： T=A+（R+ S2/24 R）tan（F/2） L2 N≤Z =V=N-Z：W=0：I=0：Goto 1 L3 U=N-Z：N≤(Z+S)=V=U-U5/40R2S2：W=U7/336 R3S3-U3/6RS：I=90U2/πRS：Goto 1 【L*{Y}：Y=0=U=U+X】 L4 N≤B =I=90（2U-S）/πR：V=RsinI+A：W=RcosI-R- S2/24R：Goto 1 L5 N≤H =U=H-N：C= U-U5/40R2S2：G= U3/