卡西欧fx-4800导线坐标计算程序2005-8-14.doc

卡西欧fx-4800计算器（2007年07月08日）一.计算线路平曲线坐标计算主程序XYZ2及子程序XYZ1： 子程序名XYZ1X=S(1-S440A4+S83456A8-S12599040A12)Y=S36A2(1-S456A4+S87040A8-S121612800A12)主程序名XYZ2BOT: ZJ K E XH YC AARL:Z=1 F=K-J: G=: D=Z=2 L=K-J:D=(AbsLR)r: F=Rsin(LR)r:G=Ftan(D2)Z=3 S=K-J:D=9RS2(A2AbsR):progXYZ1:F=X:G=YAbsRR Z=4 S=A2AbsR:P=K-J:L=S-AbsP:C=C+9A2(RAbsR):D=-9L2AbsR(RA2):progXYZ1:F=X:G=Y:S=L: progXYZ1:F=P(F-X)AbsP:G=P(Y-G)AbsR(RAbsP)C+DX=E+FcosC+Gcos(C+9)+Bcos(C+D+O)+Tcos(C+D) Y=H+FsinC+Gsin(C+9)+Bsin(C+D+O)+Tsin(C+D)说明：、 Z表示平曲线的线型：Z=1表示直线，Z=2表示圆曲线，Z=3表示缓和曲线，Z=4表示不完全缓和曲线。、 K（J）表示起始点的里程，X（E）表示其始点的X坐标，Y（H）表示其始点的Y坐标。注意：在计算器中0123456789中的0显示为，以与英语字母的OPQ中O区分开。、 A（C）表示起始点的方位角；RL（A）表示平曲线的曲线半径与缓和曲线长度乘积的开平方，K表示待求坐标点的里程。、 R表示平曲线的曲线半径R，左偏输入负值-R，右偏输入正值+R。、 C+D为待求坐标点处的切线方位角，此值为程序自动输出的中间变量，不用人工输入此值。、 B表示计算第一偏距；O表示计算第一夹角（按顺时针）；T表示计算第二偏距，前进大里程方向输入负值+T，后退小里程输入正值-T。详见下图所示： 线路前进方向 坐标放样待求点 T O B该程序可计算平曲线线路上任意一点坐标和及其线路方位角。以及与线路任意夹角任意距离处的任意一点的边桩坐标。举例如下：某曲线资料为：左偏=29-15-30，R=1600m,L=140m,T=487.35ZH点里程:DK12+343.93,ZH点坐标:(12343.930,500.000)ZH点方位角:00-00-00,HY点里程:DK12+483.930。1）.缓和曲线上点的坐标具体计算见下：例如计算DK12+441.63处(位于缓和曲线上)的线路上的坐标步骤依次为：依次按键MODE 5 2 选择程序名XYZ2 按 EXE 计算器显示 Z？输入 3按EXE计算器显示K？输入12343.93按EXE计算器显示X？输入12343.93按EXE计算器显示Y？输入500按EXE 计算器显示 A？输入0按EXE 计算器显示RL？输入(1600140)按 EXE计算器显示K？输入12441.63 按 EXE 计算器显示 R？输入(-)1600按EXE 计算器显示 C+D 1.2207699713按SHIFT 计算器显示11314.77按EXE计算器显示B？输入0按EXE 计算器显示 O？输入0 按EXE 计算器显示 T？输入0按EXE 计算器显示 X=12441.62556 按 EXE计算器显示Y=499.3061424。因此 DK12+441.63处(位于缓和曲线上)的线路上的坐标为(12441.626，499.306)。同样可计算出HY点（里程DK12+483.930）处的线路上的坐标为(12483.9032，497.9586)以及HY点线路上该处的切线方位角-23024.09（即C+D）。若计算DK12+441.63处（位于缓和曲线上）右偏0.03米法线上的墩中心坐标步骤依次为：依次按键MODE 5 2 选择程序名：XYZ2 按 EXE 计算器显示 Z？输入 3按EXE计算器显示K？输入12343.93按EXE计算器显示X？输入12343.93按EXE计算器显示Y？输入500按EXE 计算器显示 A？输入0按EXE 计算器显示RL？输入(1600140)按 EXE计算器显示K？输入12441.63 按 EXE 计算器显示 R？输入(-)1600按EXE 计算器显示 C+D1.2207699713按SHIFT 计算器显示11314.77按EXE计算器显示B？输入0.03按EXE 计算器显示 O？输入90 按EXE 计算器显示 T？输入0按EXE 计算器显示 X=12441.6262 按 EXE计算器显示Y=499.3361356。因此 DK12+441.63处(位于缓和曲线上) 右偏0.03米法线上的墩中心坐标为(12441.626，499.336)。若要计算DK12+441.63 (位于缓和曲线上) 在ZH-xy小坐标系下的坐标步骤依次为：依次按键MODE 5 2 选择程序名XYZ2 按 EXE 计算器显示 Z？输入 3按EXE计算器显示K？输入12343.93按EXE计算器显示X？输入0.000按EXE计算器显示Y？输入0.000按EXE 计算器显示 A？输入0按EXE 计算器显示RL？输入(1600140)按 EXE计算器显示K？输入12441.63 按 EXE 计算器显示 R？输入(-)1600按EXE 计算器显示 C+D 1.2207699713按SHIFT 计算器显示11314.77按EXE计算器显示B？输入0按EXE 计算器显示 O？输入0 按EXE 计算器显示 T？输入0按EXE 计算器显示 X=97.69556486 按 EXE计算器显示Y=-0.6938575846。因此DK12+441.63(位于缓和曲线上)在ZH-xy小坐标系下的坐标为(97.6956，-0.6939)。这里的计算结果与用程序名为HHQXFY的计算结果相同(详见第8页中五.缓和曲线上任一点的放样程序)。同样可计算出HY点（里程DK12+483.930）在ZH-xy小坐标系下的坐标为(139.9732，-2.0414)以及HY点线路上该处的切线方位角-23024.09（即C+D）。2）.圆曲线上点的坐标具体计算见下：例如计算DK12+539.84处(位于圆曲线上)的线路上的坐标步骤依次为：依次按键MODE 5 2 选择程序名XYZ2 按 EXE 计算器显示 Z？输入2按EXE计算器显示K？输入12483.93按EXE计算器显示X？输入12483.9032(此值来源于1)中的的计算结果)按EXE计算器显示Y？输入497.9586(此值来源于1)中的的计算结果)按EXE 计算器显示 A？输入(-) 2 30 24.09 (此值来源于1)中的的计算结果)按EXE 计算器显示RL？输入(1600140)按 EXE计算器显示K？输入12539.84 按 EXE 计算器显示R？输入(-)1600按EXE 计算器显示 C+D4.508821062按SHIFT 计算器显示43031.76按EXE计算器显示B？输入0按EXE 计算器显示 O？输入0 按EXE 计算器显示 T？输入0按EXE 计算器显示 X=12539.70561 按 EXE计算器显示Y=494.5379956。因此 DK12+539.84处(位于圆曲线上)的线路上的坐标为(12539.706，494.538)。同样若计算DK12+539.84处(位于圆曲线上)右偏0.04米法线上的墩中心的坐标步骤依次为：依次按键MODE 5 2 选择程序名XYZ2 按 EXE 计算器显示 Z？输入2按EXE计算器显示K？输入12483.93按EXE计算器显示X？输入12483.9032(此值来源于1)中的的计算结果)按EXE计算器显示Y？输入497.9586(此值来源于1)中的的计算结果)按EXE 计算器显示 A？输入(-) 2 30 24.09 (此值来源于1)中的的计算结果)按EXE 计算器显示RL？输入(1600140)按 EXE计算器显示K？输入12539.84 按 EXE 计算器显示R？输入(-)1600按EXE 计算器显示 C+D4.508821062按SHIFT 计算器显示43031.76按EXE计算器显示B？输入0.04按EXE 计算器显示 O？输入90 按EXE 计算器显示 T？输入0按EXE 计算器显示 X=12539.70876 按 EXE计算器显示Y=494.5778718。因此 DK12+539.84处(位于圆曲线上)右偏0.04米法线上的墩中心坐标为(12539.709，494.578)。若要计算DK12+539.84 (位于圆曲线上) 在ZH-xy小坐标系下的坐标步骤依次为：依次按键MODE 5 2 选择程序名XYZ2 按 EXE 计算器显示 Z？输入 2按EXE计算器显示K？输入12483.93按EXE计算器显示X？输入139.9732(此值来源于1)中的的计算结果)按EXE计算器显示Y？输入(-)2.0414(此值来源于1)中的计算结果)按EXE 计算器显示 A？输入输入(-) 2 30 24.09 (此值来源于1)中的的计算结果)按EXE 计算器显示RL？输入(1600140)按 EXE计算器显示K？输入12539.84 按 EXE 计算器显示 R？输入(-)1600按EXE 计算器显示 C+D 4.508821062按SHIFT 计算器显示43031.76按EXE计算器显示B？输入0按EXE 计算器显示 O？输入0 按EXE 计算器显示 T？输入0按EXE 计算器显示 X=195.7756148 按 EXE计算器显示Y=-5.462004424。因此DK12+539.84(位于圆曲线上)在ZH-xy小坐标系下的坐标为(195.7756，-5.4620)。这里的计算结果与用程序名为YQXFY的计算结果相同(详见第9页中六. 圆曲线上任一点的放样程序)。二反算方位角和边长程序：程序名XY-SAM X1 N Y1P X2Q Y2S= (M-P)2+(N-Q)2 ) A=tan-1(Q-N)(P-M)+(1-Abs(P-M)(P-M)90+( 1+Abs (P-M)(P-M)(1- Abs(Q-N)(Q-N)90示意图：X (X2，Y2) A S (X1，Y1) Y 举例如下：某曲线资料为：左偏=29-15-30，R=1600m,L=140m,T=487.35ZH点里程:DK12+343.93,坐标:(12343.930,500.000)。HY点里程:DK12+483.93,坐标:(12483.903,497.986)。则计算ZH点与HY点这两点间的距离S及ZH点与HY点这两点连线的方位角A的具体步骤依次为：依次按键MODE 5 2 选择程序名XY-SA 按 EXE 计算器显示 X1？输入 12343.930按EXE计算器显示Y1？输入500.000按EXE计算器显示X2？输入12483.903按EXE计算器显示Y2？输入497.986按EXE计算器显示 S=139.9874885 按 EXE计算器显示A=359.1756572按SHIFT 计算器显示3591032.37。三已知方位角和边长计算坐标程序： 程序名SA-XYM X N Y X=M+ ScosAY=N+SsinA示意图：X (X，Y) A S (X，Y) Y 举例如下：某曲线资料为：左偏=29-15-30，R=1600m,L=140m,T=487.35ZH点里程:DK12+343.93,坐标:(12343.930,500.000)。ZH点与HY点这两点间的距离S为:139.9875, ZH点与HY点这两点间连线的方位角A为: 3591032.37。则计算HY点坐标（X，Y）的具体步骤依次为：依次按键MODE 5 2 选择程序名SA-XY 按 EXE 计算器显示 X？输入 12343.93按 EXE计算器显示 Y？输入 500按EXE计算器显示 S？输入 139.9875按EXE计算器显示A？输入359 10 32.37 按EXE计算器显示 X=12483.90301 按 EXE计算器显示Y=497.9860027。四.高程计算主程序：程序名GCJKYHRE I1 F I2 K：W=Abs(E-F)：Q=WR/2：M=W/(E-F)：KJ-Q H=Y-E(J-K): Goto1KJ H=Y-E(J-K)-M(K-J+Q)2/(2R):Goto1KJ+Q H=Y+F(K-J)M(J+Q-K)2/(2R):Goto1H=Y+F(K-J): Lb11:H=H 说明：K (J)表示变坡点的桩号；H (Y)表示变坡点的高程；R表示竖曲线半径；I1(E)表示第一坡度；I2(F)表示第二坡度；K表示欲求高程点的桩号。该程序可计算线路纵坡任一点的高程。举例如下： DK12+500某线路纵坡及竖曲线资料为： +2.5% (112.50m) +1.0% DK12+000(变坡点) DK11+500 (100.00m) (95.00m)其中，DK12+000处的竖曲线资料为：R=4000m，T=30m，E=0.1125m。竖曲线的里程范围: DK11+970DK12+030。则计算DK11+980处纵坡高程的具体步骤依次为：依次按键MODE 5 2 选择程序名GC 按 EXE 计算器显示 K？输入 12000按EXE计算器显示H？输入100.00按EXE计算器显示R？输入4000按EXE计算器显示I1？输入0.01按EXE计算器显示I2？输入0.025按 EXE 计算器显示 K？输入 11980按EXE计算器显示 H=99.8125。若计算DK12+020处纵坡高程的具体步骤依次为：依次按键MODE 5 2 选择程序名GC 按 EXE 计算器显示 K？输入 12000按EXE计算器显示H？输入100.00按EXE计算器显示R？输入4000按EXE计算器显示I1？输入0.01按EXE计算器显示I2？输入0.025按 EXE 计算器显示 K？输入 12020按EXE计算器显示 H=100.5125。同样，可计算出DK11+900的纵坡高程为：99.00，DK11+990的纵坡高程为：99.95，DK12+000的纵坡高程为：100.1125，DK12+010的纵坡高程为：100.30，DK12+100的纵坡高程为：102.50。依次类推。五.缓和曲线上任一点的放样程序：程序名HHQXFYKJKRLS=K-J： X=S-S5(40(RL)2 )+ S9(3456(RL) 4 )- S13(599040(RL) 6) Y= S3(6(RL) )- S7(336(RL) 3 )+ S11(42240(RL) 5)- S15(1612800(RL) 7) A=tan-1 (YX) H= (X2+Y2)其中：K为缓和曲线上待求放样点的里程；K (J)为ZH点的里程；(S为缓和曲线上待求点与ZH点里程的差值，即S=K-J)；L为缓和曲线长；R为曲线半径，左偏输入负值-R，右偏输入正值+R。该程序只能计算缓和曲线上任一点在ZH-xy自定小坐标系下的坐标。示意图如下 ： X HY ZH Y 举例如下：某曲线资料为：左偏=29-15-30，R=1600m,L=140m,T=487.35ZH点里程:DK12+343.93,ZH点坐标:(0.000,0.000)缓和曲线上任一点在ZH-xy小坐标系下的坐标具体计算见下：例如计算DK12+441.63(位于缓和曲线上)在ZH-xy小坐标系下的坐标步骤依次为：依次按键MODE 5 2 选择程序名HHQXFY 按 EXE 计算器显示 K？输入12441.63按EXE计算器显示K？输入12343.93按 EXE 计算器显示 R？输入(-)1600按EXE计算器显示L？输入140按EXE 计算器显示 X=97.695564