推荐 CASIO fx4800Pfx5800P型计算器用于线路施工曲线中线点坐标的计算程序 精品Word下载.docx

1、Q:Q“Z=1；Y=2”：Q=1Q=F-30D2/ （RL）：H=F-90D2/ （RL）：Q=F+30D2/ （RL）：H=F+90D2/（ RL） Goto4：Lbi2:Q=1Q=F-tng-1（U/I）：H=F-O：Q= F+tng-1（U/I）：H=F+OGoto4：Lbi3:Q=1Q=F-tng-1（U/I）:H=F-（A-90（S-（C-K）2/ （RL）: Q= F+tng-1（U/I）:H=F+（A-90（S-（C-K）2/ （RL）: Goto4：Lbi4:B=90+H:H0H“QXFWJ”=B+360 H360H“QXFWJ”=H-360 H“QXFWJ”=HGoto5：L

2、bi5: X“CSD:X”=JcosQ+N+VcosB Y“CSD:Y”=JsinQ+E+VsinBGotoA 3程序说明3.1 输入已知变量 R圆曲线半径，显示R？ L缓和曲线长，显示L？ A曲线转向角，显示A？ E直缓点纵坐标，显示ZH：X=？ N直缓点横坐标，显示ZH：Y=？ F第一切线方位角，即ZH至JD的方位角，显示FWJ=？ K直缓点里程，显示ZH：LC=？3.2 计算待求量 T切线长度，显示T= S曲线全长，显示S= Z3外矢距，即JD到QZ的距离，显示E0= 3.3 输入待求变量 K输入待求（测设）点的里程，显示LC=？ V横向偏距，即测设点左、右侧外移距的偏移量，若为中线点输

3、入0；右侧输入“+”值，左侧输入“-”值。3.4 计算结果和说明X待求（测设）点的纵坐标，显示X= Y待求（测设）点的横坐标，显示Y= 该程序既可对曲线内中线任意点的坐标进行计算，又可以对中线外移桩的坐标进行计算。计算原理如下：先建立整个曲线的坐标系统，确定曲线起点（ZH）的坐标，以ZH到JD的第一切线方位角为计算方位角。 在程序计算过程中，当屏幕显示Z=1；Y=2？的提示时可根据曲线的转向，从下列两种情况中选择一种予以回答。若曲线为左转，则输入数值1；若曲线为右转，则输入数值2。 本程序式输入时，凡遇到、两种符号时，不应将其分开，应视为一个整体的符号。 由于变量的增加，程序式中使用到Z1、Z

4、2、Z3三个计算量，此时应将计算器的内存变量扩大，另行增加36个变量即可。在计算器上的操作过程为：按SHIFT键再按Defm键，然后按数字36即可。3.5 本程序在实际测量工作中的应用 本程序不仅可以计算曲线线路上任意中线点的坐标，亦可以计算线路左、右侧的法线偏距点的坐标，用于曲线中桩和边桩放样。 可以计算曲线线路中心点的切线方位角和法线方位角，该方位角可以用于桥梁墩、台（亦为曲线上桥梁墩、台截面纵、横轴线的方位角）定位。其切线方位角在计算器上显示为“QXFWJ= ”。 本程序具有多方面的计算功能，可根据仪器状况和测设对象以及工作需要而灵活应用，可在测设现场计算，减轻内业的计算时间，提高工作效

5、率。4计算举例 结合某线路的曲线资料，计算曲线线路上中桩或边桩的坐标。4.1已知要素 R=1500， L=150， AZ=30 主点里程：ZHDK8+150.000 HYDK8+300.000 YHDK8+935.398 HZDK9+085.398 4.2 建立平面控制坐标系 建立独立坐标系统，若以ZH为坐标原点，以ZHJD的切线为X轴的正方向，则ZH点的坐标可设为：X=8150.000，Y=1000.000；切线的初始方位角F“FWJ= ”可设为00000。现以计算圆曲线内DK8+400的中线点或左侧15米边桩点的坐标为例。4.3 计算操作步骤 1输入曲线半径R=1500按EXE键 2输入缓

6、和曲线长L=150按EXE键 3输入曲线转向角A=30按EXE键 4输入ZH点X坐标ZH：X=8150.000按EXE键 5. 输入ZH点Y坐标ZH：Y=1000.000按EXE键 6输入切线方位角FWJ=00000按EXE键 7输入ZH点里程ZH：LC=8150按EXE键 显示T=477.0849919（切线长）按EXE键；显示S=935.398164（曲线长）；显示E0=53.56126046（外矢距）；8输入测设点里程LC=8400；9输入中线点或左侧边桩的距离，如计算中桩点坐标输入PJU=0，如需计算左侧15米边桩的坐标则输入PJU=-15米；10屏幕显示Z=1；提示，本例为左偏曲线故

7、输入1，按EXE键,显示FXFWJ=83.31549239，按SHIFT键再按，即DK8+400中桩点的法线方位角为831855.77按EXE键11. 显示CSD：X=8399.5970289，按EXE键； 显示CSD：Y=989.1782961，按EXE键。若需计算DK8+400左侧15米边桩点的坐标，按EXE键，屏幕提示LC=8400；再按EXE键，提示PJU=？，输入-15米，再按EXE键，显示Z=1；输入数字1按EXE键，显示QXFWJ=83.31549239，按EXE键；显示CSD：X=8397.850997，按EXE键；Y=974.2802637，按EXE键。12若需要再计算其他线

8、路中线点或边桩点的坐标，重复811步骤操作，不必再输入各已知数据，仅输入需要计算点的里程或横向偏距及Z=1；的参数即可。ZBFS（坐标正算）A”X1”:B”Y1”:C”X2”:D”Y2”:Fixm：Pol（C-A,D-B）:I”S=” JGoto2L-DT=H=M+ （L-D）F +（F-C）Abs（F-C） （T-L +D）2（2R）：H=M+ （L-D）Goto3LbI2：D-LT=H=M-（D-L）C+（F-C） Abs（F-C）（T-D+ L）2H=M-（D-L）C：Goto3LbI3：HGoto1三、程序执行1、CI1 前进方向输入变坡点后段坡率2、FI2 输入变坡点前段坡率以路线前

9、进方向3、R 输入竖曲线半径。计算出T竖曲线切线长，E竖曲线外矢距。4、L 输入所求点的桩号（里程）5、DBP.LC 输入变坡点桩号（里程）6、MBP.GC 输入变坡点高程。7、H 以路线计算出所求点高程四、说明1、计算凸曲线，所求点桩号（里程）不小于后一个凹曲线临界（终点）桩号，不大于前一个凹曲线临界（起点）桩号。2、计算凹曲线，所求点桩号（里程）不小于后一个凸曲线临界（终点）桩号，不大于前一个凸曲线临界（起点）桩号。曲线任意里程中边桩坐标正反算（CASIO fx-4800P计算器）程序 一、程序功能 本程序由一个主程序（TYQXJS）和两个子程正算子程序（SUB1）、反算子程序（SUB2）

10、序构成，可以根据曲线段直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。另外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲线任意里程中边桩坐标进行正反算。本程序也可以在CASIO fx-4500P计算器及 CASIO fx-4850P计算器上运行。特别申明：（1）. 适用于弧长小于2倍半径的各种线元坐标正反算，精度优于1mm； （2）. 在引用该核心计算部分时，请注明来源。二、源程序1.主程序（TYQXJS）1.SZ = XY：2

11、.XY = SZN：UX0VY0OS0GF0HLSPR0RRNQ：C=1P：D=（P-R）（2HPR）：E=180：N=1=Goto 1：Goto 2Lbl 1：SZ：SZ：W=Abs（S-O）：Prog SUB1XXS=XYYS=YFFS=F-90Goto 1Lbl 2：XY：XY：I=X：J=Y：SUB2S=O+WZ=ZGoto 22. 正算子程序（SUB1）A=0.1739274226：B=0.3260725774：K=0.0694318442：L=0.3300094782：F=1-L：M=1-K：X=U+W（Acos（G+QEKW（C+KWD）+Bcos（G+QELW（C+LWD）+B

12、cos（G+QEFW（C+FWD）+Acos（G+QEMW（C+MWD）：Y=V+W（Asin（G+QEKW（C+KWD）+Bsin（G+QELW（C+LWD）+Bsin（G+QEFW（C+FWD）+Asin（G+QEMW（C+MWD）：F=G+QEW（C+WD）+90：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）cosL-（I-X）sinL：AbsZW=W+Z：Goto 0Z=（J-Y）sinF三、使用说明1、规定 （1） 以道路中线的前进方向

13、（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。 （2） 当所求点位于中线时，Z=0；当位于中线左铡时，Z取负值；当位于中线中线右侧时，Z取正值。 （3） 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。 （4） 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。 （5） 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替； （6） 当线元为非完整缓和曲线时，起点与直

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

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

16、25 16 31.00 269.256 1E45 0769.256 19787.340 28563.378 37.492 221.75 -1806.748 19766.566 28594.574 120 25 54.07 112.779 221.75 221.75919.527 19736.072 28701.893 91 17 30.63 80.285 221.75 9579.228999.812 19744.038 28781.659 80 40 50.00 100.000 0 1、正算 （注意：略去计算方式及线元要素输入，请自行根据所求点所在的线元输入线元要素） S=700 Z=-5 计

17、算得 XS=19831.41785 YS=28509.72590 Z=0 XS=19827.33592 YS=28506.83837 Z= 5 XS=19823.25398 YS=28503.95084 S=780 XS=19785.25749 YS=28575.02270 XS=19781.15561 YS=28572.16358 XS=19777.05373 YS=28569.30446 S=870 XS=19747.53609 YS=28654.13091 XS=19742.68648 YS=28652.91379 XS=19737.83688 YS=28651.69668 S=940

18、Z=-5.123 XS=19741. 59118 YS=28722.05802 XS=19736.47687 YS=28722.35642 Z= 3.009 XS=19733.47298 YS=28722.53168 2、 反算 X=19831.418 Y=28509.726 S=699.9999974 Z= -5 .00018164 X=19827.336 Y=28506.838 S=699.9996493 Z= 0.000145136 X=19823.25398 Y=28503.95084 S=699.9999985 Z= 5.000003137 X=19785.25749 Y=28575

19、.02270 S=780.0000035 Z= -5 .000001663 X=19781.15561 Y=28572.16358 S=780.0000025 Z=- 0.000002979 X=19777.05373 Y=28569.30446 S=780.0000016 Z= 4.99999578 X=19747.536 Y=28654.131 S=870.0001137 Z= -4.99941049 X=19742.686 Y=28652.914 S=870.0003175 Z=- 0.00041814 X=19737.837 Y=28651.697 S=870.0002748 Z= 4

20、.999808656 X=19741.5912 Y=28722.0580 S=939.9999786 Z= -5.123024937 X=19736.4769 Y=28722.3564 S=939.9999862 Z=- 0.000027710 X=19733.4730 Y=28722.5317 S=940.0000238（DMCL）断面测量（fx -4800）ABCDEKIFGJ：B：C”R1” ：D”R2” ：E”R3” ：K”H1” ：I：F”H0” ：G”S0” ：J”LC”Lbl 1：HZ：H：Z”PJ” ：X=F+（J-G）I+K：Y=X+C/2：W=X+A：U=W-Dsin（17

21、5355）：Q=W-（D-E） sin（17HX+C=Goto 3HY=P=（C2-（H-X）2）:Goto 2HU=P=B+（D2-（H-W）2）:Goto 2: P=B+（D-E） cos（175355）+（E2-（H-Q）2） :Goto 2“Z0=“:Z-PGoto 1Lbl 3：“H=“:H-（X+C）JDZB（文件名：曲线坐标计算CASIO fx-5800P）“JDX1”?A:“JDY1”?B: “JDX2”?C:“JDY2”?D:JF：“ZH：LC=”?K：L2/24R-L4/2688R3P:L/2-L3/240R2M：（R+P）*tng（A/2）+MT: “T=”:TRA/180G： G+LS: “S=”:SC+TCOS（F+180） N: “ZH.X=”:ND+TSIN（F+180） E: “ZH.Y=”:E C“DK?+”：V“PJ”：C-KD：If DL:Then D-D5/（40R2L2） I: D3/6RL-D7/（336R3L3）U:（I2+U2）J：IfEnd：IF DG:Then O=90（2D-L）/ （R）： RsinO+MI：R（1