线路计算程序

1、公 路 工 程线路计算程序汇编Fx 一 4500 、 4800 、 4850PA （卡西欧计算器）Fx4500 、4800 、4850PA（卡西欧计算器）一 坐标计算（中、边桩）、 文件名 ZBL1 DJD：A：R：CLs L2 P=C2/24R-C4/2688R3L3 Q=C/2-C3/240R2 L4 B=90C/RL5 T=（R+P）tan（Abs A/2）+Q L6 WE=（R+P）/cos（A/2）-RL7 L=(Abs（A-2B)）R /180+2CL8 GZH=D-T HHY=G+C IQZ=G+L/2 KYH=G+L-CL9 MHZ=G+L L10 N ：E ：F ：JL11

2、A<0=>S=-1：=>S=1 L12 U=F+A/2+90SL13 V=W+R L14 B=N+VcosUL15 O=E+VsinU L16 W=6.5L17 Lb1 1 L18 ZL19 ZG=>L=T+G-Z L20 V=F+180L21 U=F L22 Goto 2L23 ZH=>L=Z-G L24 V=L-L5/90（RC）2L25 L=30SL2/(RC) L26 P=F+180L27 Q=F+L L28 U=F+3LL29 Goto 4 L30 ZK=>L=F+A/2+180+90S+180（Z-I）S/RL31 U=L+90S L32 Got

3、o 5 L33 ZM=>L=M-Z L34 V=L-L5/90（RC）2L35 L=30SL2/(CR ) L36 P=F+AL37 Q=F+A+180-L L38 U=F+A-3LL39 Goto 4 L40 Z>M=>L=Z-M+T L41 U=F+A L42 V=U L43 Goto 2 L44 Lb1 2L45 X=N+LcosV L46 Y=E+LsinV L47 Goto 6 L48 Lb1 3 L49 WL50 PXL=X+Wcos（U-90） QYL=Y+Wsin（U-90） PXR=X+Wcos（U+90） QYR=Y+Wsin（U+90）L51 Goto

4、1 L52 Lb1 4L53 X=N+T cosP+VcosQ Y=E+TsinP+VsinQ L54 Goto 6L55 Lb1 5 L56 X=B+RcosL Y=O+RsinL L57 Goto 6 L58 Lb1 6 L59 J=1=>Goto 3 L60 Goto 1 （共计667步，计算器 “”为零， “O”为英文字母O） ZB坐标 JD交点桩号 Ls缓和曲线长度（单位m）A 偏角值（转角值），左偏为“-”，右偏为“+”R曲线半径 T切线长 E外矢距L曲线全长 ZH HY QZ YH HZ要素点里程桩号N交点X坐标 E交点Y坐标F交点i-1至交点i的方位角J输入“零”（）为中

5、桩坐标计算 输入“1”为边桩坐标计算X所求桩号X坐标 Y所求桩号Y坐标W输入所求边桩到中桩的垂距（单位为m） XL左边桩X坐标 YL左边桩Y坐标 XR右边桩X坐标 YR右边桩Y坐标 Z待求点桩号JD10JD9JD11ZHHZZYYZHZ 注： 以上示意图，利用JD10的数据，可以自JD9的HZ点计算到JD11的ZY点，其它类推。二 用坐标法计算面积 文件名 S L1 N:I=1 L2 X:YL3 A=X:B=Y:C=X:D=Y L4 I=Z:S=0L5 LbI 0:I:XY:X:Y L6 S=S+(CY-DX)/2L7 C=X：D=Y L8 I=N => Goto 1L9 I=I+1：G

6、oto 0 L10 LbI 1:S=S+(CB-AD)/2说明：N为输入转折点个数三.方位角距离计算程序已知两点坐标，求方位和距离） 文件名：YZ LD ZB FJL1 A ：B ：C ：D L2 X=A-C L3 y=B-D L4 L=Po1（X ，Y） L5 W=J+180 注：（W4500，J4800、4850）已知两点坐标求方位角、距离 文件名：YZ ZB WDL1 Lb1 0 L2 X ,Y L3 H=X-A:U=Y-B L4 Poi(H ,U) L5 U0=>J L6 U<0=>J+360 L7 Goto 0 注：（W4500，J4800、4850） 已知测站和待

7、求点坐标，求方位角、距离 文件名： CZ DQ ZB FLL1 Lb1 1 L2 A ：B ：N ，E L3 X=N-A Y=E-BL4 D=(X2+Y2) L5 X0=>F=tan-1（Y/X） => Goto 2L6 Lb1 2 L7 X0=> F=360+tan-1（Y/X） => Goto 3L8 Lb1 3 L9 X<0=> F=180+tan-1（Y/X） Goto 1 四.设计高程计算程序 文件名：BG1、 P：H：I：J：R 2、 W=Abs（I-J） T=R W/2 E=T2/2R 3、 Lb1 1 4、 Z,D 5、 ZP=>X=T

8、-(P-Z) =>Goto 26、 X0=>KHI= H -（P-Z）I B=K-DV=>Goto 4 7、 Goto 1 8、 Lb1 4 9、 Y=X2/2R 10. J>I=>KHI= H -（P-Z）×I+Y B=K-DV =>KHI= H -（P-Z）×I-Y 11、 Goto 1 12、 Lb1 2 13、 X=T-(Z-P) 14 X0=>KHI= H +(Z-P)×J B=K-DV =>Goto 3 15、 Goto1 16、 Lb1 3 17、 Y=X2/2R 18、 J>I=>KHI=

9、 H+(Z-P)×J+Y B=K-DV =>KHI=H+(Z-P)×J-Y B=K-DV19、 Goto1 （ 共计253步）说明： P竖曲线变坡点桩号 H 竖曲线变坡点高程（即SJD之高程） I前坡（起点方向纵坡，上为“+”，下为“-”） J后坡（终点方向纵坡，上为“+”，下为“-”） R竖曲线半径 Z所求点桩号 T-竖曲线切线长 B边桩高程 D为中到边的垂距 V为横坡度i%五、单纯圆曲线点偏角法，“付偏角法”，弦线支距法、切线支距:文件名：YQX PJ XT (计算区间：ZY YZ)1、 MZY:NYZ:R 2、 LbI 1 3、 P4、 A=180(P-M)/2

10、R 5、 =180(N-P)/2R6、 Z=(P-M)-(P-M)3/24R2 7、 Y=Zsin8、 XX1=Zcos 9、 TTX=ZcosA 10、 HTZ=ZsinA 11、 Goto 1 (共计142步)六、 单纯缓和曲线的偏角、付偏角、弦线支距、切线支距： 文件名：HQX PJF XT （计算区间：ZHHY）1、 MZH:NHY:R:FLs 2、 B=RF 3、 Lbi 1 4、 P 5、 A=30(P-M)2/B 6、 Q=30F2-(P-M)2/B7、 D=(P-M)-（P-M）5/90/R2/F2 8、 Y=DsinQ 9、 XX1=DcosQ 10、 TTX=DcosA 1

11、1、 HTZ=DsinA 12、 Goto 1 （共计144步）七、缓和曲线与圆曲线组合时的切线支距法（程序之一）计算区间ZHHYQZ文件名 LS-QXZJ(意为：切线支距)1、 AZH:KHY:CLS:R 2、 LbI 1 3、 P4、 PK=>X=(P-A)-(P-A）5/(40R2C2)Y=（P-A）3/6/R/C=>Goto 25、 Goto 1 6、 LbI 2 7、 B=90C/R 8、 F=180(P-K)/R+B 9、 D=C2/24R 10、Q=C/2-C3/240R211、 X=RsinF+Q Y=R(1-cosF)+ D 12、 Goto 1 (共计177步)

12、注：后半部分曲线需倒算，桩号需另外假设，即HZQE八、缓和曲线与圆曲线组合式的切线支距法（程序之二）（此可取代第七套程序）文件名：ZH-QE,HZ-QZ （亦即计算区间）1、 AZH:KHY:MYH:NHZ:CLS:R2、 Lb1 1 3. P4、 PK=>X=(P-A)-(P-A)5/40(RC)2 Y=(P-A)3/(6RC）=>Goto 25、 Goto 1 6、 Lb1 2 7、 B=90C/R 8、 F=180(P-K)/R+B 9、 D=C2/24R 10、 Q=C/2-C3 /240R2 11、 X=RSinF+Q Y=R(1-cosF)+D12、 PM=>XT

13、X=(N-P)-(N-P)5/40 (RC)2 YTY=(N-P)3/(6RC) =Goto 313、 Goto 1 14、 Lbi1 3 15、 V=180(M-P)/R+B 16、 XTX=RSinV+Q YTY=R(1-CosV)+D17、 Goto 1 （计300步）说明：X、Y为正(顺)算之值；TX、TY为倒(逆)算之值。即ZHQZ 为X、Y HZQZ 为TX、TY 不用改变(假设)桩号，直接输入实际里程(桩号)。也可使用纯圆曲线（P3 、P4 直接输入即可）九、直线某点中桩、（左、右）边桩坐标计算文件名：ZBJS(中、边桩坐标计算)AXO：BYO：CK0：F 1、 Lb1 1 2、

14、 P，D 3、 ZXP=A+（P-C）CosF 4、 HYP=B+（P-C）SinF 6、 JXL=Z+DCos(F-90） 7、 KYL=H+DSin(F-90）8、 UXR=Z+DCos(F+90） 9、 MYR=H+DSin(F+90）10、 Goto 1 （计137步）备注： X0为起算点X坐标； Y0为起算点Y坐标； K0为起算点桩号； P为待求点桩号； D为中到边的垂距。十、缓和曲线中、边桩坐标计算程序（ZHHY）文件名：Ls ZH HY1 DZH：TXZH：BYZH：FFZH：R：CLs1：A±2 A>O=>S=1 ：=>S=-1 3 Lb1 1 4

15、P5 偏角 E=30(P-D)2/(RC)6 弦长 Z=(P-D)-(P-D)5/90(RC)27 弦方位角 GFZHP=F+ES8 BP H=90(P-D)2 /(RC)9 P点切线方位 IFP =F+HS10 中桩坐标 X=T+ZCosG Y=B+ZSinG 11 中线到边桩的垂距 W12 左边 JXL=X+WCos(I-90) KYL=Y+WSin(I-90)13 右边LXR=X+WCos(I+90) MYR=Y+WSin(I+90)14 Goto 1 大约230步注：A路线前进方向ZHYH，左“-”、右“+”输入一个代正负号的任意数十一、缓和曲线中、边桩坐标计算程序(YHHZ)文件名：

16、Ls YH HZ DHZ：TXHZ：BYHZ：FFHZYH：R：CLs2：A±A>O=>S=1：=>S=-1 1、 Lb1 1 2、 P3、 偏角 E=30(D-P)2/(RC) 4、 弦长 Z=(D-P)-(D-P)5/90(RC)25、 弦方位角 GFFPYH=F + ES6、 BP H=90(D-P)2 /(RC)7、 P点切线倒方位 IFFPYH=F + HS8、 中桩坐标 X=T+ZcosG Y=B+ZsinG 9、 中线到边桩的垂距 W10、 左边 JXL=X+Wcos(I+90) KYL=Y+Wsin(I+90)11、 右边 LXR=X+Wcos(I-

17、90) MYR=Y+Wsin(I-90)12、 Goto 1 大约238步注：1、第二缓和曲线需倒算，倒算时显示的左侧边桩，实为路线前进方向的右侧边桩；2、“FHZYH”=图纸上FHZ±180，即HZ点切线去YH的方位。 3、 A路线前进方向YH-HZ，左“+”、右“-”输入一个代正负号的任意数十二、复曲线中间插入缓和曲线(卵形曲线)中、边桩坐标计算程序文件名：LXQXDX (由大D小X) BXR:DXA1:GYH:HXYH:IYYH:JFYH:KHY:A± A>0=>S=1:= S=-1 NXLs1=D2/B MZH=K-N1、Lb1 1 2、 P 3、L&q

18、uot;X1"=(PG)20(GM)2(PG)3+15(GM)(PG)4+3(PG)5120(BN)24、O"Y1"=(PG)3/(6BN)+(GM)(PG)2/（2BN）5、Q=30(PG)PM+2(GM)/(BN) 6、F"FC"=J+QS 大到小，右用“+”，左用“-”，易程序 7、X=H+(L2+O2) cosF Y=I+(L2+O2) sinF 8、T"FP"=J+S(P-G)(P-G)+2(G-M)/(2BN)×（180 /)9、W 注：S有正负之分，大小，左偏为“-“,右偏为“+”（易程序） 10、U

19、XL=X+WCos(T-90) 11、VYL=Y+Wsin(T-90)12、CXR=X+Wcos(T+90) 13、ZYR=Y+Wsin(T+90) 14、Goto 1 大约380步左右 显示注明：DR-大圆半径 XR-小圆半径 XA1-小圆缓和曲线参数A-大到小（路线前进方向），左偏为“-”，右偏为“+”（任意正、负数）X,Y-任意待求点P的中桩坐标， W-中到边之垂距。XL YL-左边坐标 XR YR-右边坐标 对应于P点小圆第一缓和曲线长度 XLs1=XA12/XR小圆半径ZH点里程： ZH=HY-XLs1 注明：该程序计算方向为：从大圆到小圆；反之，倒算。倒算时所显示的左右边桩坐标，实

20、质为路线前进方向的左右边桩坐标。十三、中间缓和曲线（卵形曲线）中边桩坐标计算程序文件名：LXQX-X-D (意为：卵形曲线由小大) BXR:DXA2:GHY:HXHY:IYHY:JFHY:KYH:A±A>0=>S=1:=>S=-1 NXLs2=D2/BMHZ=K+N 注：FHY=图纸上FHY±180，即HY点去小圆方向的切线方位1、 Lb1 1 2、P3、 LX1=(G-P) 20(MG)2(GP)3+15(MG)(GP)4+3(GP)5/120(BN)24、 O"Y1"=(PG)3/(6BN)+(MG)(GP)2/（2BN）5、 Q=

21、30(GP)M-P+2(MG)/(BN) 6、 F"FC"=J+QS 小到大，左用“+”，右用“-”，易程序 7、 X=H+(L2+O2）cosF Y=I+ (L2+O2）sinF 8、 T"FP"=J+S(G-P)(G-P)+2(M-G)/(2BN)×(180/)9、 W 注：S有正负之分，小 大，左偏为“+”,右偏为“-”（序）易程 10、 UXL=X+WCos(T+90) 11、 VYL=Y+Wsin(T+90)12、 CXR=X+Wcos(T-90) 13、 ZYR=Y+Wsin(T-90) 14、 Goto 1 大约380步左右 显示

22、的左右边桩，乃实际路线前进方向之左右A小倒大路线前进方向，左偏为“+”、右偏为“-”中间缓和曲线（大圆小圆）计算理论及过程：小圆自HY或YH向ZH或HZ方向以回旋线的形式渐变，必有一点要与大圆相接于大圆的YH或HY，就形成YHHY一段曲线，称其为中间缓和曲线或称卵形曲线。中、边桩任意点坐标计算步骤：1、 从图上查得：Xr小圆半径, Xa1、Xa2-小半径圆曲线两端缓和曲线参数 YH、HY里程桩号, XYH 、YYH YH点坐标, FYHYH点的切线方位 A视左、右偏，大小，小大而确定正、负值任意一数；2、以YH为原点、其切线为X1轴，垂直切线方向为Y1轴，建立直角坐标系，求出待求点P的切线支距

23、（X1，Y1）;3、 求得YH到P点偏角,YH到P点的弦长,YH到P点的弦线方位角;4、 计算P点的中桩坐标X,Y;5、 求出P6、 求出P点的切线方位FP = FYH±1直线外一点P(X ，Y)到直线和直线上某两点1、2的距离文件名：PZXLDJL A ：B ：C ：D 1、 LbI 1 2、 X ，Y3、 L=(A-C)2+(B-D)2 4、 K=(A-X)2+(B-Y)25、 M=(C-X)2+(D-Y)2 6、 H=Cos-1(K2+L2-M2)/2KL 7、 Q "X1"=Cos H.×K 8、 T"Y1"=Sin H.&#

24、215;K 10、 Goto 1LY1X1HMK2#(C,D)1#(A,B)P(X,Y)十四、单圆曲线中桩、边桩坐标计算程序文件名: ZYYZ ZB kzy:AXzy:BYzy:FFzy:R:O2. O>0S=1:S= 1 备注: Q为 ±号3. Lbi 1 P 备注: P为待求点桩号 Q=180(Pk)R 备注: Q为P点的圆心角 IFp=F+QS 备注: I为P点切线方位角 N=RSinQ E=R(1CosQ) 备注: N 、E为P点支距 X=A+NCosFES×SinF Y=B+NSinF+ESCosF 备注: X 、Y为P点坐标 边桩坐标 D D为中桩到边桩的

25、垂距 ZXL=X+DCos(I90) LYL=Y+DSin(I90) 备注: Z 、L为左边桩坐标 HXR=X+DCos(I+90) MYR=Y+DSin(I+90) 备注: H 、M为右边桩坐标16、Goto 1十五、隧道掌子面计算程序1文件名：SDZZM1L1 ZR1: YR2 : XR3: O : B : C L2 VH1=Z(1-Cos45°)L3 UH2=Z-C L4 TH3=UYSin（O-45°）L5 SH4=U(Y-X)Sin（O-45°）XCos（135°-O -B）L6 QH5 = (Y-X)Sin(O -45°) L7 PH6 = (Y-X)Cos(O -45°)L8 LbI 1 L9 H L10 H<0=>L=0 =>Goto 2 L11 Goto 1L12 LbI 2 L13 HV=>L=(Z2 -(Z- H ) 2) =>Goto 3L14 Goto 1L15 LbI 3 L16 H U =>L=C（Y2 -（Z - H -C) 2） =>Goto 4L17 Goto 1 L18 LbI 4L1