不等长缓和曲线计算.docVIP

5800不等长缓和曲线计算 “XLZBJSCX” ◢LB1??1 ↙CLS : FIX 4 : 30→DIM Z ↙“K=”?K ◢（计算里程）IF??K本段曲线终点里程??AND??K≥上段曲线终点里程 ：THEN??本段终点里程→Z[1] : 上段曲线终点里程→Z[2] ：曲线偏角→A : 1→I （注：左偏曲线输入-1→I,右偏曲线输入1→I）: 半径→R : 第一缓和曲线→Z[6] : 第二缓和曲线→Z[7] : T1切线→Z[12]: T2切线→Z[13]: 交点 X→B :交点 Y→C : 小里程向交点方位角→E : 交点向大里程方位角→F : GOTO??2 : IFEND…………LB1??2 ↙（曲线要素计算）Z[6]/2- Z[6]^3/(240*R^2)+ Z[6]^5/(34560*R^4) →Z[8] ??（M1）Z[7]/2- Z[7]^3/(240*R^2)+ Z[7]^5/(34560*R^4) →Z[9] ??（M2）Z[6]^2/(24*R)- Z[6]^4/(2688*R^3) →Z[10] ??（P1）Z[7]^2/(24*R)- Z[7]^4/(2688*R^3) →Z[11] ??（M2）π*A*R/180+0.5*( Z[6]+ Z[7])→S （曲线长）90* Z[6]/(R*π) →Z[14] （第一缓和曲线总偏角）90* Z[7]/(R*π) →Z[15] （第二缓和曲线总偏角）B+ Z[12]*COS (E+180)→ Z[16] （ZH点X）C+ Z[12]*SIN(E+180)→ Z[17] （ZH点Y）Z[1]-S→Z[3] (ZH点里程)Z[1]-S+ Z[6]→Z[4] (HY点里程)Z[1]- Z[7]→Z[5] (YH点里程)GOTO 3 LB1??3 ↙(判断里程点与曲线关系)IF??K≤Z[3]??AND??K Z[2] : THEN??GOTO??4 : IFEND IF??K≤Z[4]??AND??K Z[3] : THEN??GOTO??5 : IFEND ↙IF??K≤Z[5]??AND??K Z[4] : THEN??GOTO??6 : IFEND ↙IF??K≤Z[1]??AND??K Z[5] : THEN??GOTO??7 : IFEND ↙LB1??4 ↙（里程小于直缓点直线独立坐标）K- Z[3] →X : 0→Y : E→T : PROG“TYZBCX” ：GOTO??1 LB1??5 ↙（第一缓和曲线独立坐标）K- Z[3] →H H-H^5/(40*R^2* Z[6]^2)+H^9/(3456*R^4* Z[6]^4) →X ↙H^3/(6*R* Z[6])-H^7/(336*R^3* Z[6]^3) →Y ↙90*H^2/( R*π* Z[6]) →T ↙IF??I0??：THEN??T +E→T : ELSE??E-T →T : T0=360+T→T : IFEND PROG“TYZBCX” ：GOTO??1 LB1 6 ↙（圆曲线独立坐标）K- Z[4] →H H*180/( R*π)+ Z[14]→T ↙R*SIN T+ Z[8]→X ↙R*(1-COS T)+ Z[10]→Y ↙IF??I0??：THEN??T +E→T : ELSE??E-T →T : T0=360+T→T : IFEND PROG“TYZBCX” ：GOTO??1 LB1 7 ↙（第二缓和曲线独立坐标）Z[1] -K??→H H-H^5/(40*R^2* Z[7]^2)+H^9/(3456*R^4* Z[7]^4) →U↙H^3/(6*R* Z[7])-H^7/(336*R^3* Z[7]^3) →V ↙90*H^2/( R*π* Z[7]) →T ↙Z[13]COS A+ Z[12]-U*COS A-V*SIN A→X ↙Z[13]*SIN A-U*SIN A+V*COS A→Y ↙IF??I0??：THEN??F-T→T : T0=360+T→T : ELSE??F+T →T : IFEND PROG“TYZBCX” ：GOTO??1 子程序：“TYZBCX” （统一坐标计算）IF??I0??: THEN??-Y→Y : IFEND “QXJ=” :T◢（切线方位角）Z[16]+X*COS E-Y*SIN E→Z[18] Z[17]+X*SIN E+Y*COS E→Z[19] ↙“XI=” : Z[18] ◢（中线X）“YI=” : Z[19] （中线Y）“PJ=”?J(输入边桩与线路夹角，左-右+)