曲线坐标计算程序

合集下载

CASIO任意曲线坐标计算程序及说明

CASIO任意曲线坐标计算程序及说明

个人测量资料汇集韦习耀2005年8月起缓和曲线(ZH到HY)中桩坐标CASIO4800/4850计算程序A"ZHX":B"ZHY":C"LS":Q:R:K:M"DK ZH"↙LB1 1:{L}:L"DK":X=A+((L-M)-(L-M)5÷(40R2C2))cosQ-((L-M)3÷(6RC))KsinQ◢Y=B+((L-M)-(L-M)5÷(40R2C2))sinQ+((L-M)3÷(6RC))KcosQ ◢Goto 1说明:(A,B)为ZH点坐标,LS为缓和曲线长,Q为ZH点放位角,R 为半径,左转K为-1,右转K为1,M为ZH点里程,L为所求点的里程。

缓和曲线(HY到HZ)中桩坐标CASIO4800/4850计算程序A"HZX":B"HZY":C"LS":Q:R:K:M"DK HZ"↙LB1 1:{L}:L"DK":X=A+((-(L-M))-(-(L-M))5÷(40R2C2))cosQ-((-(L-M))3÷(6RC))KsinQ◢Y=B+((-(L-M))-(-(L-M))5÷(40R2C2))sinQ+((-(L-M))3÷(6RC))KcosQ◢Goto 1说明:(A,B)为HZ点坐标,LS为缓和曲线长,Q为HZ点放位角,R 为半径,左转K为-1,右转K为1,M为HZ点里程,L为所求点的里程。

另外注意,这个程序是把HZ点当作ZH点,倒过来算,方位角Q 和K均应注意其值。

圆曲线(HY到YH)中桩坐标CASIO4800/4850计算程序A"HYX":B"HYY":C"LS":Q:R:K:M"DK HY"↙LB1 1:{L}:L"DK":X=A+2Rsin(90(L-M)÷(πR))cos(Q +90((L-M)+C)÷(πR)K)◢Y=B+2Rsin(90(L-M)÷(πR))sin(Q+90((L-M)+C)÷(πR)K)◢Goto 1说明:(A,B)为HY点坐标,LS为缓和曲线长,Q为ZH点放位角,R 为半径,左转K为-1,右转K为1,M为HY点里程,L为所求点的里程。

坐标计算程序

坐标计算程序

测量实用技术程序一:线路坐标计算程序采用的公式1:曲线要素公式内移值P=Ls xy2/(24R)切线增量q=Ls/2-Ls xy3/(240R xy2)切线长T=(R+p)tg(A/2)+q切线角B0=Ls/(2R)曲线外距E=(R+P)SEC(A/2)-R2:缓和曲线任意一点与ZH点之间的弦长公式C=L-L xy5/(90R xy2Ls xy2)3:缓和曲线任意一点的偏角公式S=(L/LS)xy2S0S0=1/3B04:其他采用坐标增量公式△X=LCOSW△Y=Lsinw然后与JD的坐标相加X=XJD+△XY=YJD+△Y其中L为计算点至交点JD的长度W为JD到计算点的方位角XJD为交点的X坐标YJD为交点的Y坐标圆曲线增设缓和曲线示意图程序功能可计算直线,圆曲线,缓和曲线上任意一点的中桩坐标及其法线边桩坐标A:曲线偏角 T切线长 E外失距 R 园曲半茎ZH 直缓点 HY缓园点 HY 缓园点 QZ曲中 YH 缓园 HZ 缓直线路直曲线坐标计算程序1:A:C”LS=”:D”JD=”:R”R=”2: P=C xy2/24/R-C xy4/2688/R xy3 3:Q=C/2-C xy3/240/R xy24: B=90C/R/∏5:T=(R+P)Tg(absA/2)+Q6: E=(R+P)/COS(A/2)-R7;L=(absA-2B) ∏R/180+2C 8: G”ZH”=D-T▲ H” HY”=G+C 9: I”QZ”=G+L/2 10:K”YH”=G+L-C 11:M”HZ”=G+L 12:N”JDX=”:W”JDY=”:F”FANGWEIJIAO”:J 13:A<0→S=-1: S=1⊿14:U=F+A/2+90S 15:V=E+R 16:B=N+VCOSU 17:O=W+VsinU 18:Lb1 19:{Z} 20:Z≤ G →L=D-Z 21:V=F+18022:U=F23:Goto 2 ⊿24:Z≤H→L=Z-G25:V=L-L xy5/90/(RC) xy226:L=30SL xy2/R/∏/C27:P=F+18028:Q=F+L29:U=F+3L30:Goto5⊿31:Z≤K→L=F+A/2+180+90S+180S(Z-I)/R/∏32:U=L+90S33:Goto5⊿34:Z≤M→L=M-Z35:V=L-L xy5/90(RC) xy236:L=30Sl xy2/(R∏C)37:P=F+A38:Q=P+180-L39:U=Q-2L-18040:Goto4⊿41:Z>M→L=Z-M+T42:V=F+A43:U=V44:Goto2⊿45:Lb1 246:X =N+LcosV▲47:Y =W+LsinV▲48: Goto6⊿49: Lb1 350:{E}51:P”XL”=X+E cos(U-90) ▲52:Q”YL”=Y+E sin(U-90) ▲53:P”XR”=X+E cos(U+90) ▲54:Q”YR”=Y+E sin(U+90) ▲55: Goto 156: Lb1 457:X=N+T cosP+V cosQ▲58:Y=W+T sinP+V sinQ▲59: Goto 660: Lb1 561: X=B+R cosL▲62: Y=O+R sinL▲63; Goto 664: Lb1 665:J=1→ Goto 3⊿66:Goto1 说明:A曲线偏角(左偏为负数)LS缓和曲线长D JD交点里程R园曲线半径N JD的X坐标W JD的Y坐标F是起算方位角ZH→JD Z待算里程 J 输入1时计算中桩及两侧坐标输入0仅计算中桩坐标E输入中桩到两侧的距离XL YL 是中桩左侧的坐标XR LR是右侧坐标已知两点坐标计算方位角程序1:X”XA=”:Y”YH=”:S”XB=”:T”YB=”2:A=(S-X)3:B=(T-Y)4:POL(A,B) ▲5:W≥0→W=W: W=W+360⊿6:F”FWJ”=W▲。

CASIO 5800综合曲线坐标计算程序(使用中)

CASIO 5800综合曲线坐标计算程序(使用中)

C A SI O f x-5800P综合曲线坐标计算程序说明:该程序适用于计算器C ASI O f x-5800P,可计算线路中心的缓和曲线、圆曲线、直线段,中、边桩坐标及切线方位角。

1、α? 输入转角:左转为负,右转为正2、R?输入圆曲线半径3、LS? 输入缓和曲线长度4、J D(DK)?输入交点里程桩号5、X(JD)?输入本交点X 坐标6、Y(JD)?输入本交点Y 坐标7、F W J?输入待求点切线方位角9、J? 输入0 程序计算中桩,输入 1 程序计算边桩10、JSDK? 输入里程桩号主程序名:ZHQ X1. Deg:Fix 4↙2. Lb1 A↙3. "α"?A:?R:"LS"?C:"JD(DK)"?D:"X(JD)"?N:"Y(JD)"?E:"FWJ"?F↙4. Lb1 B↙5. ?J ↙6. C∧2÷24÷R-C∧(4)÷2688÷R∧(3)→P↙7. C÷2-C∧(3)÷240÷R∧2→Q↙8. 90×C÷兀÷R→B↙9. "T=":(R+P)tan(Abs(A)÷2)+Q→T◢10. "E=":(R+P)÷c o s(A÷2)-R→W◢11. "L=":(Abs(A)-2×B)兀R÷180+2×C→L◢12. "ZH=":D-T→G◢13. "HY=":G+C→H◢14. "QZ=":G+L÷2→I◢15. "YH=":G+L-C→K◢16. "HZ=":G+L→M◢17. If A<0:Th en-1→S:E l se1→S:If End↙(本行中0 为数字)18. F+A÷2+90×S→U↙19. W+R→V↙20. N+V×cos(U)→B↙21. E+V×sin(U)→O↙(本行中O 为字母)22. Lbl 1 ↙23. "JS DK"?Z↙24. If Z≤G:Th en T+G-Z→L ↙25. F+180→V ↙26. "FW J=":F→U◢27. Goto 2 ↙28. IfEnd ↙29. If Z≤H:Th en Z-G→L ↙30. L-L∧(5)÷(90×R∧2×C∧2)→V↙31. 30×L∧2×S÷(兀RC)→L ↙32. F+180→P ↙33. F+L→Q↙34. "FW J=":F+3×L→U◢35. Goto 4 ↙36. IfEnd ↙37. If Z≤K:Th en F+A÷2+90×S+180+180(Z-I)×S÷R÷兀→L↙38. "FW J=":L+90×S→U:If U≥360:The n U-360→U◢39. Else U→U◢40. IfEnd:G o to5:If End↙41. If Z≤M:Th en M-Z→L↙42. L-L∧(5)÷(90×R∧2×C∧2)→V↙43. 30×S×L∧2÷(兀RC)→L ↙44. F+A→P ↙45. F+A+180-L→Q↙46. "FW J=":F-3×L+A→U:If U≥360:Th en U-360→U◢47. Else U→U◢48. IfEnd:G o to4:If End↙49. If Z>M:Then Z-M+T→L↙50. F+A→U ↙51. U→V ↙52.G o to2 ↙53. Lbl 2 ↙54. "X=":N+L×cos(V)→X◢55. "Y=":E+L×sin(V)→Y◢56. G o to6 ↙57. If End↙58. Lbl 3 ↙59. ?W60. "XL=":X+W×cos(U-90)→P◢61. "YL=":Y+W×sin(U-90)→Q◢62. "XR=":X+W×cos(U+90)→P◢63. "XL=":Y+W×sin(U+90)→Q◢64. G o to1 ↙65. Lbl 4 ↙66. "X=":N+T×cos(P)+V×cos(Q)→X◢67. "Y=":E+T×sin(P)+V×cos(Q)→Y◢68. G o to6 ↙69.Lbl 5 ↙70. "X=":B+R×cos(L)→X◢71. "Y=":O+R×sin(L)→Y◢(本行中O 为字母)72. G o to6 ↙73. Lbl 6 ↙74. J=1 G o to3 ↙75. G o to1 ↙。

CASIO+FX-5800P+综合曲线坐标计算程序(正确)2010.1.17

CASIO+FX-5800P+综合曲线坐标计算程序(正确)2010.1.17

圆曲线坐标计算交点桩号YK20+734.498;N=2889769.961;E=527002.305 ;方位角314'24'58.3' 转角27'19'31.4' 切线长度T1= 814.345 ;T2= 814.347 ;半径R=3000 ;曲线长度LS1:170;LC=1260.752;LS2= 170 曲线总长=1600.752 求ZH,HY,QZ,YH,HZ以及圆心点的坐标。

1、A?输入转(偏)角;左转为负,右转为正。

转角27'19'31.4'2、R?输入曲线半径。

30003、LS?输入缓和曲线长度。

1704、JD(DK)?输入交点里程桩号。

20734.498结果显示T = 814.347 W=87.776 L=1600.754 ZH=19920.151;HY=20090.151; QZ=20720.528; YH=21350.905; HZ=21520.9055、X(JD)?输入本交点X坐标。

2889769.9616、Y(JD)?输入本交点Y坐标。

527002.3057、FWJ?输入直缓(ZH)点切线方位角。

314'24'58.3'8、J?输入0程序计算中桩,输入1程序计算边桩。

09、Z?输入里程桩号。

ZH=19920.151FWJ:314.416 X=2889200.028 Y= 527583.973完(FWJ:314.416=314'24'57.6' )ZH(2889200.026 527583.975) HY(2889320.140 527463. 682)QZ(2889816.373 527076.808) YH(2890382.383 526801.950)HZ(2890543.311 526747.167)例如2:JD15:4013872.8781;456031.1661 JD16:4015349.321;451808.075ZH:75+137.901 ay:16.45.49.69 R=9000 缓和曲线长:490 T=1571.257 JD(DK) =76709.158 FWJ=289-16-12.85 给我算一下74+400;75+300;75+800;76+500左中线坐标。

曲线坐标计算程序VBA-4.8

曲线坐标计算程序VBA-4.8

单位名称:工程名称:
坐 标 计 算 表
**集团有限公司**路线**标段
单位名称:工程名称:
坐 标 计 算 表
**集团有限公司**路线**标段
单位名称:工程名称:
坐 标 计 算 表
**集团有限公司**路线**标段
单位名称:工程名称:
坐 标 计 算 表
**集团有限公司**路线**标段
单位名称:工程名称:
坐 标 计 算 表
**集团有限公司**路线**标段
单位名称:工程名称:
坐 标 计 算 表
**集团有限公司**路线**标段
单位名称:工程名称:
坐 标 计 算 表
**集团有限公司**路线**标段
单位名称:工程名称:
坐 标 计 算 表
**集团有限公司**路线**标段
单位名称:工程名称:
坐 标 计 算 表
**集团有限公司**路线**标段
单位名称:工程名称:
坐 标 计 算 表
**集团有限公司**路线**标段
单位名称:工程名称:
坐 标 计 算 表
**集团有限公司**路线**标段
单位名称:工程名称:
坐 标 计 算 表
**集团有限公司**路线**标段
单位名称:工程名称:
坐 标 计 算 表
**集团有限公司**路线**标段
单位名称:工程名称:
坐 标 计 算 表
**集团有限公司**路线**标段
单位名称:工程名称:
坐 标 计 算 表
**集团有限公司**路线**标段
单位名称:工程名称:
坐 标 计 算 表
**集团有限公司**路线**标段。

任意曲线坐标放样正反算程序(4800,4850,5800均可以用)

任意曲线坐标放样正反算程序(4800,4850,5800均可以用)

任意曲线坐标放样正反算程序(fx-4800p)正算ZSLbl 0↙{K,L,Q,S}:K:L:Q:Prog“FXL”:Prog“SUB”↙“X=”:X=X:Pause0:Y=Y◢K=S+K↙Goto0↙反算FSLbl 0:{PT}:P“X0”:T“Y0”↙Lbl 1:{K}:L=0:Q=0:Prog“FXL”: Prog“SUB”↙U=I-90↙V=(T-Y)CosU-(P-X)SinU↙AbsV≥0.001=>K=K+V:Goto1:≠>“K0=”:K=K:Pause0:W“L0”=(T-Y)CosI-(P-X)SinI◢Goto0↙子程序SUBNorm↙Deg↙U=(E-D)÷Abs(G-F):V=Abs(K-F):W=UV↙J=C+90V(W+2D)÷π:J<0=>J=J+360⊿J≥360=>J=J-360⊿↙M=C+45V(W÷4+2D)÷(2π):N=C+135V(3W÷4+2D)÷(2π):O=C+45V(W÷2+2D)÷π↙I=J+Q↙Fix3↙X=A+V÷12×(CosC+4(CosM+CosN)+2CosO+CosJ)+L CosI↙Y=B+V÷12×(SinC+4(SinM+SinN)+2SinO+SinJ)+L SinI↙数据库FXLK≥0=>K≤1211.753=>A=9764.547:B=3344.148:C=142,6,26,:D=0:E=0:F=0:G=1211.753⊿⊿↙K≥1211.753=>K≤1406.455=>A=8808.279:B=4088.39:C=142,6,26,:D=1÷700:E=1÷700:F=1211.753:G=1406.455⊿⊿↙使用说明:K-里程桩号; A,B-每段线元的起点X,Y坐标; C-每段线元的起点的切线方位角; F-每段线元的起点桩号;G-每段线元的终点桩号;D-每段线元的起点半径曲率;E-每段线元的终点半径曲率; (左转取-1÷半径曲率,右转取1÷半径曲率,直线为下一个桩号(例: 赋值20确认,即上个桩号加20)D,E取值说明:当线元为直线段,不论其起,止与何线元相接, D,E都取0.当线元为圆曲线,不论其起,止与何线元相接, D,E都取(右偏)1/(左偏)-1÷该圆曲线的半径.当线元为缓和段, 起,止与直线相接, D,E都取0. 与圆曲线相接, D,E都取(右偏)1/(左偏)-1÷该圆曲线的半径. 与缓和段相接, D,E都取0.任意点的切线方位角C说明:编数据库时,若两曲线相接,例:先编ZH至HY段,运行算出HY点的切线方位角,(即求HY点的坐标,然后终止程序调出C值),然后再以计算出的C编HY至YH段.依此类推.反算桩号时需多次确认所输入的X,Y坐标,程序运算至桩号误差小于允许值时,会同时显示结果:K0-桩号,L0-距离.坐标反算程序:ZBFSC“X1”:D“Y1”:E“X2”:F“Y2”:Fixm:Pol(E-C,F-D:I“S=”◢J≤0=>J=J+360⊿J“J=”↙坐标正算程序ZBZSC“X1”:D“Y1”:L:O“J”:Fixm:X=C+Rec(L,O◢Y=D+J↙竖曲线SQXJ“BPD”:B“BPH”:C“I1”:D“I2”:RW=D-C◢A=W÷AbsW:T=Abs(0.5RW)◢L=2T↙E=T2÷2R◢X“QD”=J-T◢Y“ZD”=J+T◢Lbl1:P“K”:H=B+C(P-J)+A(P-J+T)2÷2R◢W“CONT”:W=1=>Goto1:≠>Goto2:⊿↙Lbl2↙说明:BPD-变坡点桩号BPH-变坡点高程I1-前一个坡度(上“+”,下“-”)I2-后一个坡度(上“+”,下“-”)R-半径E-外距QD-该竖曲线起点桩号ZD-该竖曲线终点桩号K-要求的桩号H-要求桩号的高程任意曲线坐标放样正反算程序(fx-4850p)正算ZSLbl 0↙{K,L,Q,S}:K:L:Q:Prog“FXL”:Prog“SUB”↙“X=”:X=X:Pause0:“Y=”:Y=Y◢K=S+K↙Goto0↙反算FSLbl 0:{PT}:P“X0”:T“Y0”↙Lbl 1:{K}:L=0:Q=0:Prog“FXL”: Prog“SUB”↙U=I-90↙V=(T-Y)CosU-(P-X)SinU↙AbsV≥0.001=>“K0=”:K=K+V:Goto1:≠>“K0=”:K=K:Pause0:W=(T-Y)CosI-(P-X)SinI ↙“L0=”:W=W◢Goto0↙子程序SUBNorm↙Deg↙U=(E-D)÷Abs(G-F):V=Abs(K-F):W=UV↙J=C+90V(W+2D)÷π:J<0=>J=J+360⊿J≥360=>J=J-360⊿↙M=C+45V(W÷4+2D)÷(2π):N=C+135V(3W÷4+2D)÷(2π):O=C+45V(W÷2+2D)÷π↙I=J+Q↙Fix3↙X=A+V÷12×(CosC+4(CosM+CosN)+2CosO+CosJ)+L CosI↙Y=B+V÷12×(SinC+4(SinM+SinN)+2SinO+SinJ)+L SinI↙数据库FXLK≥0=>K≤1211.753=>A=9764.547:B=3344.148:C=142,6,26,:D=0:E=0:F=0:G=1211.753⊿⊿↙K≥1211.753=>K≤1406.455=>A=8808.279:B=4088.39:C=142,6,26,:D=1÷700:E=1÷700:F=1211.753:G=1406.455⊿⊿↙使用说明:K-里程桩号; A,B-每段线元的起点X,Y坐标; C-每段线元的起点的切线方位角;D-每段线元的起点半径曲率;E-每段线元的终点半径曲率; (左转取-1÷半径曲率,右转取1÷半径曲率,直线为0. );L-距中桩距离(左负,右正,中桩为0); Q-与线路前进方向的右交角(正交为90); S-赋值加上一个桩号放下一个桩号(例: 赋值20确认,即上个桩号加20)D,E取值说明:当线元为直线段,不论其起,止与何线元相接, D,E都取0.当线元为圆曲线,不论其起,止与何线元相接, D,E都取(右偏)1/(左偏)-1÷该圆曲线的半径.当线元为缓和段, 起,止与直线相接, D,E都取0. 与圆曲线相接, D,E都取(右偏)1/(左偏)-1÷该圆曲线的半径. 与缓和段相接, D,E都取0.任意点的切线方位角C说明:编数据库时,若两曲线相接,例:先编ZH至HY段,运行算出HY点的切线方位角,(即求HY点的坐标,然后终止程序调出C值),然后再以计算出的C编HY至YH段.依此类推.反算桩号时需多次确认所输入的X,Y坐标,程序运算至桩号误差小于允许值时,会同时显示结果:K0-桩号,L0-距离.任意曲线坐标放样正反算程序(fx-4800p)任意曲线坐标放样正反算程序(fx-4850p)。

铁路曲线桥梁支座坐标计算程序

铁路曲线桥梁支座坐标计算程序

铁路曲线桥梁支座坐标计算程序
铁路曲线桥梁支座坐标计算程序是一个用于计算铁路曲线桥梁支座坐标的程序。

它根据给定的曲线参数以及支座的位置和数量,计算出每个支座的坐标。

程序通常包括以下几个步骤:
1. 输入曲线参数:用户需要输入铁路曲线的参数,包括曲线半径、曲线长度、过渡曲线长度等。

2. 输入支座位置和数量:用户需要输入桥梁支座的位置和数量。

3. 计算支座坐标:根据输入的曲线参数和支座位置,程序会计算出每个支座的坐标。

4. 输出结果:程序会将计算结果输出,通常以表格形式呈现。

5. 绘制支座位置示意图:程序可以根据计算结果绘制一个支座位置示意图,以便用户更直观地了解支座的位置。

需要注意的是,铁路曲线桥梁支座坐标计算程序可能会有不同的实现方式和算法,具体的实现细节可能会有所不同。

程序可以使用各种编程语言编写,例如C++、Python等。

CASIO5800CYT综合曲线坐标计算程序V1.0.

CASIO5800CYT综合曲线坐标计算程序V1.0.

[作者原创] CASIO fx-5800P CYT 综合曲线坐标计算程序V1.0说明:本程序适用于卡西欧计算器 CASIO fx-5800P,能进行全线贯通计算。

可计算:第一夹直线段、第一缓和曲线段、圆曲线段、第二缓和曲线段、第二夹直线段的线路中、边桩坐标及切线方位角,该程序主要针对不对称缓和曲线设计以及对现场施工放样减少了计算方位角和距离的繁琐。

主程序名:CYT第1步50→DimZ:Deg:Fix 6第2步"1=>XY":"2=>DK":?J:If J=1:Then Goto A:IfEnd:If J=2:Then Goto B:IfEnd第3步Lb1 A:"JSDK"?Z:"W"?P:"α"?Q:Prog"YT":Goto C:Lb1 B第4步"QJD(X"?U:"QJD(Y"?V:"JD(X"?N:"JD(Y"?E:"JD(DK"?D:?R:"LS1"?S:"LS2"?T:"α"?A 第5步"JSDK"?Z:"W"?P:"α"?Q:Lb1 C第6步 POL(N-U,E-V:J→F:IF J<0:Then J+360→F:IfEnd第7步S^(2÷24÷R-S^(4÷2688÷R^(3→Z[5]:T^(2÷24÷R-T^(4÷2688÷R^(3→Z[6]第8步S÷2-S^(3÷240÷R^(2→Z[7]:T÷2-T^(3÷240÷R^(2→Z[8]第9步90×S÷π÷R→Z[9]:90×T÷π÷R→Z[10]第10步(Abs(A-Z[9]-Z[10]×R×π÷180→Z[2]第11步(R+Z[5]×tan(Abs(A÷2+(Z[6]-Z[5]÷sin(Abs(A+Z[7]→Z[11]第12步(R+Z[6]×tan(Abs(A÷2+(Z[5]-Z[6]÷sin(Abs(A+Z[8]→Z[22]第13步(R+Z[5]÷COS(Abs(A÷2-R→Z[1]第14步Z[2]+S+T→L:D-Z[11]→Z[12]:Z[12]+S→Z[13]:Z[13]+Z[2]→Z[15]:Z[12]+L→Z[16]第15步Cls:"T1=":Locate 4,1,Z[11]:"T2=":Locate 4,2,Z[22]:"E=":Locate4,3,Z[1]:"L=":Locate 4,4,L◢第16步Cls:"ZH=":Locate 4,1,Z[12]:"HY=":Locate 4,2,Z[13]:"YH=":Locate4,3,Z[15]:"HZ=":Locate 4,4,Z[16]◢第17步F+180→Z[17]:If Z[17]>360:Then Z[17]-360→Z[17]:IfEnd第18步N+Z[11]×Cos(Z[17]→Z[18]:E+Z[11]×Sin(Z[17]→Z[19]第19步F+A→Z[20]:If Z[20]>360:Then Z[20]-360→Z[20]:IfEnd第20步N+Z[22]×Cos(Z[20]→Z[21]:E+Z[22]×Sin(Z[20]→Z[23]第21步 If A<0:Then -1→H:Else 1→H:IfEnd:Goto 2:Lb1 1第22步"JSDK"?Z:"W"?P:"α"?Q:Lb1 2第23步If Z≤Z[12]:Then Goto D:IfEnd:If Z≤Z[13]:Then Goto E:IfEnd第24步If Z≤Z[15]:Then Goto F:IfEnd:If Z≤Z[16]:Then Goto G:IfEnd第25步If Z≥Z[16]:Then Goto H:IfEnd:Lb1D:If Z≤Z[12]:Then Z[12]-Z→G:IfEnd第26步Z[18]+G×Cos(Z[17]+P×Cos(F+Q→X:Z[19]+G×Sin(Z[17]+P×Sin(F+Q→Y:F→Z[27]第27步Goto 3:Lb1 E:If Z≤Z[13]:Then Z-Z[12]→W:IfEnd:If A<0:Then -1→H:Else 1→H:IfEnd第28步 (Z-Z[12]^(2÷2÷R÷S×180÷π→Z[26]第29步F+Z[26]×H→Z[27]:W-W^(5÷40÷R^(2÷S^(2→Z[28]:W^(3÷6÷R÷S→Z[29]第30步 Z[18]+Z[28]×Cos(F-H×Z[29]×Sin(F+P×Cos(Z[27]+Q→X第31步Z[19]+Z[28]×Sin(F+H×Z[29]×Cos(F+P×Sin(Z[27]+Q→Y第32步Goto 3:Lb1 F:If Z≤Z[15]:Then Z-Z[12]→K:IfEnd:If A<0:Then -1→H:Else 1→H:IfEnd第33步 ((Z-Z[12]-0.5×S÷R×180÷π→Z[32]第34步F+Z[32]×H→Z[33]:R×Sin(Z[32]+Z[7]→Z[34]:R×(1-Cos(Z[32]+Z[5]→Z[35]第35步 Z[18]+Z[34]×Cos(F-H×Z[35]×Sin(F+P×Cos(Z[33]+Q→X第36步Z[19]+Z[34]×Sin(F+H×Z[35]×Cos(F+P×Sin(Z[33]+Q→Y:Z[33]→Z[27]第37步Goto 3:Lb1 G:If A>0:Then -1→H:Else 1→H:IfEnd:If Z≤Z[16]:ThenZ[16]-Z→M:IfEnd第38步Z[20]+180→Z[42]:If Z[42]>360:Then Z[42]-360→Z[42]:IfEnd第39步M^(2÷2÷R÷T×180÷π→Z[38]:Z[20]+Z[38]×H→Z[39]第40步 M-M^(5÷40÷R^(2÷T^(2→Z[40]:M^(3÷6÷R÷T→Z[41]第41步 Z[21]+Z[40]×Cos(Z[42]-H×Z[41]×Sin(Z[42]+P×Cos(Z[39]+Q→X第42步Z[23]+Z[40]×Sin(Z[42]+H×Z[41]×Cos(Z[42]+P×Sin(Z[39]+Q→Y:Z[39]→Z[27]第43步Goto 3:Lb1 H:If Z≥Z[16]:Then Z-Z[16]→O:IfEnd第44步Z[21]+O×Cos(Z[20]+P×Cos(Z[20]+Q→X:Z[23]+O×Sin(Z[20]+P×Sin(Z[20]+Q→Y:Z[2 0]→Z[27] 第45步 Goto 3:Lbl 3第46步 Cls:"F=":Locate 3,1,Z[27]°:"X=":Locate 3,2,X:"Y=":Locate 3,3,Y◢第47步"Z(X"?B:"Z(Y"?C第48步 POL(X-B,Y-C:If J<0:Then J+360→J:Else J→J:IfEnd第49步 Cls:"F=":Locate 3,1,J°:"S=":Locate 3,2,I◢第50步 Goto 1数据库名:YTGoto 1Lbl 1:If Z<10000:Then3378226.731→U:456053.721→V:3373351.385→N:458628.379→E:9557.84→D:2040→R:260→S:280→T:-14°18ˊ10″→A:Return:IfEndLbl 1:IfZ<…………………………………………………………………………………………………………………………………………………………………………………………Lbl 1:If Z<本交点HZ里程或下个交点ZH里程:Then 前交点坐标X→U:前交点坐标Y→V:本交点坐标X→N:本交点坐标Y→E:交点里程→D:曲线半径→R:第一缓和曲线长→S:第二缓和曲线长→T:转角→A:Return:IfEndLbl 1每增加一行则为增加一个曲线要素,……表示下一个曲线要素。

适合隧道测量用的各种小程序

适合隧道测量用的各种小程序

曲线坐标计算程序:Deg:Norm1↓(设置计算器角度单位为度及小数点显示为小数模式)主程序(程序名称:QXZBJS)R :K:W:C↓(定义变量)K÷2- K3÷240R2→Z[1]:K2÷24R- K4÷2384R3→Z[2] :(Z[2]+R)tan(Abs(W)÷2)+ Z[1]→Z[3] :K+ Abs(W)×πR÷180→Z[4]↓Lbl 1↓Abs(D-C)→I↓If D<C:Then D-C→X:0→Y :0→Q :Goto 4:IfEnd↓If I> Z[4]:Then Z[4]-I→Z[5]:0→Z[6] :0→Q :Goto 3:IfEnd↓If I≤Z[4]-K:ThenI→J:Else Z[4]-I→J:IfEnd↓If J≥K:Then Goto 2:IfEnd↓(J2×180÷(2πRK))°→Q↓(以度数为单位赋值给Q)J- J5÷(40R2K2)→X:J3÷(6RK)- J7÷(336R3K3)→Y↓If I>Z[4]-K:Then X→Z[5]:Y→Z[6]:Goto 3:IfEnd↓Goto 4↓Lbl 2↓((I- 0.5K)×180RπR)°→Q↓(以度数为单位赋值给Q)Rsin(Q)+ Z[1]→X:R(1-cos(Q)) + Z[2]→Y↓Goto 4↓Lbl 3↓Z[3](1+cos(W))- Z[5]cos(W)- Z[6]sin(Abs(W))→X↓(Z[3]- Z[5]) sin(Abs(W)) + Z[6]cos(W)→Y↓Abs(W)-Q→Q↓Goto 4↓Lbl 4↓If W<0:Then -Q→Q:-Y→Y:IfEnd↓X-Esin(Q)→X:Y+Ecos(Q)→Y↓Goto 5↓Lbl 5↓F:G:ZX- Z[3]→I:Y→J↓F+Icos(Z)-Jsin(Z)→X↓G+Isin(Z)+Jcos(Z)→Y↓Goto 6↓Lbl 6↓“X=” :X◢“Y=” :Y◢Pol(X-A, Y-A) ↓“I=” :I◢If J>0:Then“J=” :J◢“J=” :J+360◢IfEnd↓子程序:左线(QL)Lbl 0↓?A:?B↓Goto 1↓Lbl 1↓?D:?E↓If D>220249.506:Then Goto 2:IfEnd↓219215.613→C:1450→R:280→K:-29°47′23″→W↓54117.378→F:94561.972→G: 298°7′23″→Z↓Goto3↓Lbl 2↓If D>222040.807:Then “WU XIAO” ◢Goto 1:IfEnd↓221331.463→C:3000→R:0.00000000000001→K:-13°32′51″→W↓54060.248→F:92598.282→G: 268°20′1″→Z↓Goto3↓Lbl 3↓Prog“QXZBJS” ↓Goto1↓RETURN↓R:圆曲线半径K:缓和曲线长度W:转角(左偏取负值,右偏取正值)C:ZH点里程Z:第一切线方位角A:置镜点X坐标B:置镜点Y坐标F:交点X坐标G:交点点Y坐标D、E:待求点里程和偏距路线为分离式路线时,右线子程序与左线一样,根据右线的直曲转角表输入相应的参数即可导线测量:DXCLLbl 0↓Deg:Norm1↓(设置计算器角度单位为度及小数点显示为小数模式)?A:?B:“A1”→Z[8]:“B1”→Z[9] ↓Pol(A-Z[8], B-Z[9]) ↓If J>0:Then“ J”:J→Z[10]◢Else“J”: J+360→Z[10]◢IfEnd↓Goto 1↓Lbl 1↓?V:?L↓Z[10]+V-180→T↓“X”:A+Lcos(T) →X◢“Y”:B+Lsin(T) →Y◢Goto 1↓A、B:分别为置镜点X坐标和Y坐标“A1” :“B1” 分别为后视点X坐标和Y坐标V:夹角L:边长缓和曲线的已知坐标反算里程和偏距程序主程序(程序名称:HHQXFS)Lbl 0↓Deg:Norm1↓(设置计算器角度单位为度及小数点显示为小数模式)Z[11] :Z[12] :Z[13] :Z[14] :M :N :R :K:C :Z↓(定义变量)Goto 1↓Lbl 1Pol(Z[13]-Z[11] ,Z1[4]-Z[12]) ↓I→Z[15] ↓If J>0:Then J→Z[16]:Else J+360→Z[16]:IfEnd↓Pol(X-Z[11] ,Y-Z[12]) ↓I→Z[17] ↓If J>0:Th en J→Z[18]:Else J+360→Z[18]:IfEnd↓Z[17]×cos(Abs(Z[16]-Z[18])) →Z[19]:Z[17]×sin(Abs(Z[16]-Z[18]))→Z[20]↓Z[19]×K÷Z[15]→SS- S5÷(40R2K2)→Z[21]:S3÷(6RK)- S7÷(336R3K3)→Z[22]↓Goto 2↓Lbl 2↓→Z[23]:tan-1 (Z[22] ÷Z[21])→Z[24]↓If N=1:Then Z+(-1)M×Z[24]→W:Else Z+180- (-1)M×Z[24]→W:IfEnd↓Z[11]+ Z[23] cosW→Z[25] :Z[12]+ Z[23] sinW→Z[26]↓→Z[27]↓→Z[18]↓→Z[29]↓Goto 3↓Lbl 3↓If Z[18]> Z[16]:Then- Z[20]- Z[29]→E:Else Z[20]- Z[29]→E:IfEnd↓“D”:C-(-1) N×S+ →D◢(计算显示反算出的里程)其中为修正系数If M=1:Then“E”:E◢Else “E”:-E→E◢IfEnd↓(计算显示反算出的偏距)子程序:(如当我们想反算慕家梁隧道左线第二缓和曲线时)Lbl 0↓54102.075→Z[11] ↓(HZ点X坐标)94035.965→Z[12] ↓(HZ点Y坐标)54101.208→Z[13] ↓(HY点X坐标)94315.848→Z[14] ↓(HY点Y坐标)1→M:2→N:1450→R:280→K:220249.506→C:268°20′01″→Z↓?X:?Y (待反算点X坐标和Y坐标)Go to 1↓Lbl 1↓Prog“HHQXFS”(引用主程序“HHQXFS”)Goto 0↓Return↓(引用主程序后必须有的结束语句)变量含义:Z[11]、Z[12]、Z[13]、Z[14]:分别为缓和曲线的起终点X坐标和Y坐标;M :缓和曲线左偏M=1,右偏M=2N :计算第一缓和曲线时M=1,第二缓和曲线时M=2R :圆曲线半径K:缓和曲线长C:缓和曲线的起点或终点里程(第一缓和曲线时为起点里程,反之为终点里程);Z:第一切线方位角在实际工作中,当项目变换后,只需根据新项目中的曲线坐标表中给定参数改变子程序中的这些变量(Z[11] :Z[12] :Z[13] :Z[14] :M :N :R :K:C :Z)即可圆曲线反算程序:YQXFSLbl 0↓Deg:Norm1↓(设置计算器角度单位为度及小数点显示为小数模式)Z[30] :Z[31] :Z[32] :Z[33] :R :L:K:Z:C↓(定义变量)Goto 1↓Lbl 1↓(180K÷(2 πR))°→Z[34]↓Z+ ×Z[34]+ ×90→Z[35]↓Z[30]+R cos(Z[35])→Z[36]↓Z[31]+Rsin(Z[35])→Z[37]↓→Z[38]↓Pol(Z[30]- Z[36],Z[31]- Z[37])↓J→Z[39]↓Pol(Z[32]- Z[36],Z[33]- Z[37])↓J→Z[40]↓Pol(Z[X]- Z[36],Y- Z[37])↓J→Z[41]↓Abs(Z[41]- Z[39])×(L-2K)÷Abs(Z[40]- Z[39]) →S↓“D”:C+S →D◢Z[38]-R→E↓If M=1:Then“E”:E◢:Else“E”:-E→E◢IfEnd↓子程序:(如当我们想反算慕家梁隧道左线圆曲线时)Lbl 0↓53993.248→Z[30]:94775.117→Z[31]↓54101.208→Z[32]↓94315.848→Z[33]↓1450→R:1033.893→L:298°7′23″→Z:219495.613→C↓?X:?Y↓Goto 1↓Lbl 1↓Prog“YQXFS”(引用主程序“YQXFS”)Goto 0↓Return↓(引用主程序后必须有的结束语句)Z[30]、Z[31]、Z[32]、Z[33]:分别为圆曲线的起终点X坐标和Y坐标;R:圆曲线半径L:曲线长K:缓和曲线长Z:第一切线方位角高程 GCLbl 0↓Deg:Norm1↓Z[41] :Z[42] :Z[43] :C:R↓(定义变量)Z[41]-Z[42]→Z[44]↓R×Abs(Z[44])÷200→T↓Goto 1↓Lbl 1↓If D<C-T:Then Goto2:IfEnd↓If D≤C:Then Goto3:IfEnd↓If D≤C+T:Then Goto4:IfEnd↓Lbl 2↓“H”:Z[43]-(C-T)×Z[41]÷100→H◢↓Goto 6↓Lbl 3↓(D+T-C) ÷2R→Y↓If Z[44]<0:ThenY→Y:Else-Y→Y:IfEnd↓“H”:Z[43]-(C-T)×Z[41]÷100+Y→H◢↓Goto 6↓Lbl 4↓(C+T-D) ÷2R→Y↓If Z[44]<0:ThenY→Y:Else-Y→Y:IfEnd↓“H”:Z[43]+(D-C)×Z[42]÷100+Y→H◢↓Got o 6↓Lbl 5↓“H”:Z[43]+(D-C)×Z[42]÷100→H◢↓Goto 6↓Lbl 6↓子程序:(左线)Lbl 0↓?D↓Goto 1↓Lbl 1↓If D>221443.412:Then Goto2:IfEnd↓-2.2→Z[41]: -0.5→Z[42]:525.222→Z[43]:219640→C:45000→R↓Goto 3↓Lbl 2↓If D>222000:Then “WU XIAO” ◢Goto 0:IfEnd↓-0.5→Z[41]: -2.473→Z[42]:514.972→Z[43]:221690→C:25000→R↓Goto 3↓Lbl 3↓Prog“GC”(引用主程序“HHQXFS”)Goto 0↓Return↓(引用主程序后必须有的结束语句)Z[41]、Z[42]:W 、WZ[43] :交点高程C:交点里程桩号R:数曲线半径直线反算程序:(ZXFS)“A”?→Z[45]:“B”?→Z[46]↓:?C:?ZLbl 1↓?X:?Y↓Pol(X-Z[45], Y-Z[46])J-Z→Z[47]↓C+Icos(Z[47])→D↓“D”:D◢“E”:Isin(Z[47])→E◢Goto 1↓放样程序:(FY)Deg:Norm1↓(设置计算器角度单位为度及小数点显示为小数模式)?A:?B↓Goto 1↓Lbl 1↓?X:?Y↓Pol(X-A, Y-B)↓“I”:I◢If J>0:Then“ J”:J◢Else“J”: J+360◢IfEnd↓Goto 1↓添加额外变量:例如想添加100个额外变量100→FUNCTION→{PROG}→SHIFT→(Dimz) →EXE当显示器上显示“Done”时,表示添加了指定数量的额外变量。

卡西欧5800曲线坐标计算程序-正确

卡西欧5800曲线坐标计算程序-正确

CASIO4800改5800p坐标计算程序-角度以度为单位1、该程序适用于计算器CASVO fx-4800P,可计算与线路中心成任意夹角的缓和曲线、圆曲线、直线段中、边桩坐标及切线方位角。

2、A?输入转角:左转为负,右转为正3、R?输入圆曲线半径4、LS?输入缓和曲线长度5、JD?输入交点桩号6、J?输入0程序计算中桩,输入1程序计算边桩7、N?E?输入交点X、Y坐标8、FW 待求点切线方位角9、当计算交角为60°时,则公式改为(U+60),(U-120)1 ?A:?R:“LS”?C:“JD”?D2 C∧2/24/R-C∧4/2688/R∧3→P 求缓和曲线内移、值3 C/2-C∧3/240/R∧2→Q 切距差4 90C/兀/R→B 缓和曲线角度β05 (R+P)tan(AbsA/2)+Q→T 切线长6 (R+P)/cos(A/2)-R→W 外矢距7 ((AbsA)-2B)兀R/180+2C→L 缓和曲线长8 D-T→G 直缓点里程9 G+C→H 缓圆点里程10 G+L/2→I 曲中点里程11 G+L-C→K 圆缓点里程12 G+L→M 缓直点里程13 ?N:?E:?F:?Jf为方位角14 If A<0:Then -1→S:Else 1→S:IfEnd15 F+A/2+90S→U16 W+R→V17 N+VcosU→B18 E+VsinU→O19 Lbl 120 ?Z21 If Z≤G:Then T+G-Z→L22 F+180→V23 F→U24 Goto 225 IfEnd26 If Z≤H:Then Z-G→L27 L-L∧5/(90R∧2C∧2)→V28 30L∧2S/(兀RC)→L29 F+180→P30 F+L→Q31 “FW=“:F+3L→U◢32 Goto 433 IfEnd34 If Z≤K:Then F+A/2+90S-180+180(Z-I)S/R/兀→L35 “FW=“:L+90S→U36 Goto 537 IfEnd38 If Z≤M:Then M-Z→L39 L-L∧5/(90R∧2C∧2)→V40 30SL∧2/(兀RC)→L41 F+A→P42 F+A+180-L→Q43 F-3L+A→U44 Goto445 IfEnd46 If Z>M:Then Z-M+T→L47 F+A→U48 U→V49 Goto 250 IfEnd51 Lbl 252 “X=“:N+LcosV→X◢53 “Y=“:E+LsinV→Y◢54 Goto 655 Lbl 356 ?W57 “XL=”:X+Wcos(U-90) ◢58 “YL=”:Y+Wsin(U-90) ◢59 “XR=”:X+Wcos(U+90) ◢60 “YR=”:Y+Wsin (U+90) ◢61 Goto 162 Lbl 463 ”X=“:N+TcosP+VcosQ→X◢64 ”Y=“:E+TsinP+VsinQ→Y◢65 Goto 666 Lbl 567 “X=”:B+R cosL→X◢68 “Y=”:O+RsinL→Y◢69 Goto 670 Lbl 671 J=1=> Goto 372 Goto 1注:◢为输出指令,只要在后面加上,即把前面的结果输出在屏幕上。

Casio fx4800任意曲线程序

Casio fx4800任意曲线程序

Casio fx-4800(4850)计算器程序任意曲线坐标计算程序(HHQX)LbI 0 :{QD}←┘G=I=>Goto 1ΔG=-1=>Goto 2△G=0=>Goto 4 △←┘LbI 1 :L=C÷R :M=L+Q-K :F=A+90(M2-L2)÷(πC):F◢T=A-90L2/πC:Goto 3△←┘LbI 2 :L=C÷R :M=L+K-Q ::T=A+90L2÷πC+180 :F=A-90(M2-L2)÷πC:F◢Goto 3△←┘LbI 3 :N=(M-L)-(M^5-L^5)÷40C2+(M^9-L^9)÷3456C^4-(M^13-L^13)÷599040C^6:E=(M^3-L^3)/6C-(M^7-L^7)÷336C^3+(M^11-L^11)÷42240C^5-(M^15-L^15)÷9676800 C^7 :X=I+NcosT-GEsinT+Dcos(F-90)◢Y=J+NsinT+GEcosT+Dsin(F-90)◢Goto 0←┘LbI 4 :L=Q-K :N=Rsin(180L÷πR):E=R(1-cos(180L÷πR)):F=A+180L÷πR:F◢X=I+NcosA-EsinA+Dcos(F-90)◢Y=J+NsinA+EcosA+Dsin(F-90)◢Goto 0说明:当G=1时,计算第一缓和曲线或由大圆到小圆的卵型曲线;G=-1时,计算第二缓和曲线或由小圆到大圆的卵型曲线;G=0时,计算直线或圆曲线。

R为起算点半径,线路左转时R、C取负值,右转时R、C取正值,计算第一缓和曲线时R取10^45。

计算第二缓和曲线时R取圆曲线半径,计算直线时R取10^45.此计算思路是沿路线的前进方向依次进行分段计算。

不对称平曲线坐标计算程序

不对称平曲线坐标计算程序

平曲线坐标计算程序(CASIO fx-4800P)A子程序:X>M⇒Prog“B”:Goto 0:≠⇒X=M⇒Prog“C”:≠⇒X<M⇒Prog“B”:Q=Q+180:Lb1 0:Q=Q+360:Q≥360⇒Q=Q-360◢≠⇒Q=Q◢B子程序:Q“Q”=tan-1((Y-N)÷(X-M)C子程序:Y≥N⇒Q=90◢≠⇒Q=270◢D子程序:C“SO”=√((X-M)2+(Y-N)2)主程序:R“R”S“LS1”U“LS2”W“ZG”P“WZH”:A=90S÷Rπ:B=90U÷Rπ:C=S2÷24R:D=U2÷24R:E=S÷2-S^3÷240R2:F=U÷2-U^3÷240R2:Q“T1”=E+(R+C)tan(AbsW÷2)+(D-C)÷sin AbsW◢K“T2”=F+(R+C)tan(AbsW÷2)-(D-C)÷tan AbsW◢O“S”=(AbsW-A-B)Rπ÷180+S+U◢J“JD”T“XJD”V“YJD”:A“ZH”=J-Q◢B“HY”=A+S◢E“YH”=A+O-U◢F“HZ”=A+O◢G“XZH”=T+Q cos(P+180)◢H“YZH”=V+Q sin(P+180)◢I“XHZ”=T+K cos(P+W)◢J“YHZ”=V+K sin(P+W)◢M“XO”N“YO”:Prog“A”:Prog“D”◢Lb1 O:{Z}:Z“ZHUANG HAO”:{D}:D:{T}:T“JG”:Lb1 0:Z≤A⇒Goto 1:≠⇒Z >A⇒Z≤B⇒L=Z-A:Goto 2:≠⇒Z>B⇒Z≤E⇒L=S:Goto 2:≠⇒Z>E⇒Z≤F⇒L=F-Z:Goto 2:≠⇒Z>F⇒C=Z-F:Goto 8△Lb1 1:C=A-Z:X=G+C cos(P+180):Y=H+C sin (P+180):V=P:Goto 9△Lb1 2:Z>A⇒Z≤E⇒K=RS:Goto Y:≠⇒Z>E⇒Z≤F⇒K=RU:Goto Y△Lb1 Y:K=0⇒V=P:X=G:Y=H:Goto 6△Goto 3:Lb1 3:X=L-L^5÷40K2+L ^9÷3456K^4-L^13÷599040K^6:Y=L^3÷6K-L^7÷336K^3+L^11÷42240K^5-L ^15÷9676800K^7:C=√(X2+Y2):V=90L2÷πK:Q=tan-1(Y÷X):Goto A△Lb1 A:W<O⇒Q=-Q:V=-V•GotoB△Lb1 B:Z>A⇒Z≤E⇒Goto 5:≠⇒Z>E⇒Z≤F⇒Goto 7△Lb1 5:X=G+C cos(P+Q):Y=H+C sin(P+Q):V=P+V :Goto C△Lb1 C:Z>A ⇒Z≤B⇒Goto 9:≠⇒Z>B⇒Goto 6△Lb1 6:∠L=Z-B:Q=90L÷πR:C=2R sinQ:Goto D△Lb1 D:W<O⇒Q=-Q:Goto E△Lb1 E:X=X+C cos(V+Q):Y=C sin(V+Q)+Y:V=V+2Q:Goto 9△Lb1 7:X=I+C cos(P+W+180-Q):Y=J+C sin(P+W+180-Q):V=P+W-V:Goto 9△Lb1 8:X=I+C cos(P+W):Y=J+C sin(P+W):V=P+W:Goto 9△Lb1 9:X=X+D cos(V+T):“X=”:Pause 0:Y=Y+D sin(V+T)◢Prog“A”:Prog“D”◢Goto O程序执行输入说明R:输入半径Ls1:缓和曲线第一段长度Ls2:缓和曲线第二段长度ZG:转角即路线偏角(左转为负,右转为正)WZH:后段切线方位角(按路线前进方向,在路线方向后为后段切线)T1:第一段切线长度T2:第二段切线长度S:平曲线长度(包括缓和曲线和圆曲线)JD:交点桩号XTD:交点桩号X轴坐标XTD:交点桩号Y轴坐标ZH:直缓点桩号HY:缓圆点桩号YH:圆缓点桩号HZ:缓直点桩号XZH:直缓点X轴坐标YZH:直缓点Y轴坐标XHZ:缓直点X轴坐标YHZ:缓直点Y轴坐标XO:置仪点(测站点)X轴坐标不用输YO:置仪点(测站点)Y轴坐标不用输X:后视点X轴坐标不用输Y:后视点Y轴坐标不用输Q:方位角(置仪点至后视点方位角,当输入桩号“ZHUANG HAO”即所求点桩号时,会计算置仪点至放样点的方位角)SO:距离(置仪点至后视点距离,当输入桩号“ZHUANG HAO”即所求点桩号时,会计算置仪点至放样点的距离,适用于利用经纬仪+测距仪使用)ZHUANG HAO:输入所求点桩号。

计算器曲线计算程序

计算器曲线计算程序

线路测设坐标计算程序(4800-fx Casio )在铁路、公路工程勘测和施工测量中,相对复杂的线路测设工作占了很大比重,且传统线路测设公式已不能满足现代化曲线的需要。

为此,笔者根据工程测量的特点,严格定义了坐标系,推导了线路测设的通用计算公式。

该公式不但可以计算单交点平曲线坐标,同时还可以计算复曲线坐标,如不对称缓和曲线、不完整缓和曲线等。

根据该公式,笔者编写了计算器上的单交点平曲线中桩、边桩坐标计算程序。

该程序可以判断已知数据输入的正确性,可限制输入的测设点里程。

在程序编制过程中,已考虑了零溢出问题。

一.使用说明1. 起动名为“ZXZB ”的计算程序。

2. 输入已知数据:X(QD)—起点的纵坐标(EXE)Y(QD)—起点的横坐标(EXE)X(JD)—起点的纵坐标(EXE)Y(JD)—起点的横坐标(EXE)S(JD)—交点里程(EXE)X(ZD)—终点横坐标(EXE)Y(ZD)—终点纵坐标(EXE)RY —圆曲线半径(EXE)L1 —前缓和曲线长度(EXE)RS —前缓和曲线起始半径(0表示无穷大)(EXE)L2 —后缓和曲线长度(EXE)RE —后半缓和曲线结束半径(0表示无穷大)(EXE)S(HZ1) —起点缓直点里程(EXE)S(ZH1) —交点直缓点里程(EXE)S(HZ2) —交点缓直点里程(EXE)S(ZH2) —终点直缓点里程(EXE)S? —测设点里程(EXE)X=*** —显示测设点纵坐标(EXE)Y=*** —显示测设点横坐标(EXE)I? —边桩距(左负右正)(EXE)XB=*** —边桩纵坐标(EXE)YB=*** —边桩横坐标(EXE)I? —输入0结束边桩计算转入下一测设点计算4800-fx Casio二.计算举例1.已知数据:X(QD)=2175.199 Y(QD)=1015.467X(JD)=2000 Y(JD)=2000S(JD)=6000X(ZD)=2621.088 Y(ZD)=2783.741RY=600L1=220 RS=5800L2=180 RE=3300S(HZ1)=5200 S(ZH1)=5631.047S(HZ2)=6311.047 S(ZH2)=67002.求:各测设点坐标及右边桩20M、左边桩10M坐标3.三.源程序清单ZXZB1. Defm 15:Rad:Norm↙2. A"X(QD)":B"Y(QD)":C"X(JD)":D"Y(JD)":Z"S(JD)":E"X(ZD)":F"Y(ZD)":N"RY"↙3. Q"L1":R"RS":V"L2":W"RE":G"S(HZ1)":H"S(ZH1)":K"S(HZ2)":L"S(ZH2)":Fixm↙4. O=H+Q:P=K-V:M=Pol (C-A,D-B):T=J:Z[3]=J: Z[1]=C+Rec (Z-H,J+π):Z[2]=D+J↙5. M=Pol(E-C,F-D):Z [12]=J:T=J-T:Abs T>π=>T>π=>T=T-2π: ≠>T=T+2π△△↙6. T>0=>Z[13]=1:≠>Z[13]=-1△↙7. R=0=>Z[14]=0:≠>Z[14]=QN÷(R-N)△↙8. W=0=>Z[15]=0:≠>Z[15]=VN÷(W-N)△↙9.Q=0=>I=0:≠>I=(2QZ[14]+Q²) ÷2÷N÷( Z[14]+Q)△↙10.V=0=>J=0:≠>J=(2V(Z[15]+V)-V²)÷2÷N÷( Z[15]+V)△↙11.Z[6]=Z[3]+I Z[13]: Z[9]=Z[12]-JZ[13] ↙12.T=Abs T-I-J:J=TN:Abs(P-O-J)>0.05=>"ERROR"▲Goto 1△↙13.M=Q:Prog "ZHD":Z[4]=X:Z[5]=Y:M=J:Porg "HYD":Z[7]=X:Z[8]=Y:M=V:Prog "YHD": Z[10]=X:Z[11]=Y↙14.Prog "ZXZS"↙15.Lbl 1:0:"":Pause 0ZXZS1.Rad:Norm:Fixm2.Lbl 2:{S}:S<G=>"RANGE OUT"▲Goto 5△↙3.S>L=>"RANGE OUT"▲Goto 5△↙4.S<H=>M=Z[3]:X=Z[1]+Rec(S-H,M):Y=Z[2]+J:Goto 3△↙5.S<O=>M=S-H:Prog "ZHD":Goto 3△↙6.S<P=>M=S-O:Prog "HYD":Goto 3△↙7.S<K=>M=S-P:Prog "YHD":Goto 3△↙8.M=Z[12]:X=Z[10]+Rec(S-K,M):Y=Z[11]+J↙9.Lbl 3:I"X"=Int (1000X+0.5)÷1000▲J"Y"=Int(1000Y+0.5)÷1000▲10.Lb1 4:I=0:{I}:I=0=>Goto 2△↙11.T=X+Rec(I,M+π÷2):U=Y+J:T"XB"=Int (1000T+0.5)÷1000▲U"YB"=Int(1000U+0.5)÷1000▲Goto 4↙12.Lbl 5:0:"":Pause 0ZHD1.Fixm:M<0.001=>T=0:U=0:M=Z[3]:Goto 6△↙2.I=Z[14]:X=N(I+Q):Y=Z[13]:T=M-(20I²+15IM+3M²)M^3÷120÷X²:T=T+((1008I+1680M)I^3+1080I²M²+(315I+35M)M^3)M^5÷120960÷X^4↙3.U=((70I+84M)I²+(35I+5M)M²)M^4÷1680÷X^3:U=(3I+M)M²÷6÷X-U:U=U+((3696I+7920M)I^4+(6930I+3080M)I²M²+(693I+63M)M^4)M^6÷2661120÷X^5:U=UY↙4.M=(2IM+M²)÷2÷X:I=Z[3]:M=I+MY↙5.Lbl 6:X=Z[1]+Tcos I-Usin I:Y=Z[2]+Tsin I+Ucos IHYD1.Fixm:M=M÷N:I=Z[6]:Y=Z[13]:T=Nsin M:U=(N-Ncos M)Y:M=I+MY2.X=Z[4]+Tcos I-Usin I:Y=Z[5]+Tsin I+Ucos IYHD1.Fixm:M<0.001=>T=0:U=0:M=Z[9]:Goto 7△↙2.I=Z[15] +V:X=NI:Y=Z[13]:T=M-(20I²-15IM+3M²)M^3÷120÷X²:T=T+((1008I-1680M)I^3+1080I²M²-(315I-35M)M^3)M^5÷120960÷X^4↙3.U=((70I-84M)I²+(35I-5M)M²)M^4÷1680÷X^3:U=(3I-M)M²÷6÷X-U:U=U+((3696I-7920M)I^4+(6930I-3080M)I²M²+(693I-63M)M^4)M^6÷2661120÷X^5:U=UY ↙4.M=(2IM-M²)÷2÷X:I=Z[9]:M=I+MY↙5.Lbl 7:X=Z[7]+Tcos I-Usin I:Y=Z[8]+Tsin I+UCos I注: 1.由于本计算器程序容量为4500步,本程序只能计算单交点平曲线;2.当已知数据输入2-17步存在错误时,能够自动判断错误后退出;3.当输入的测设点里程超出范围进时,计算自动终止;4.编制程序过程中已经考虑了零溢出问题。

曲线坐标计算程序

曲线坐标计算程序

曲线坐标计算程序关键词:曲线坐标计算EXCEL编程坐标曲线坐标实例摘要:利用EXCEL强大的函数功能通过曲线坐标计算的知识编制成曲线计算坐标的计算程序。

简单的输入曲线的里程桩号,通过坐标旋转、平移结合可以快速的计算完成与线路成任意角度的曲线上各中桩、边桩以及任意点坐标的计算。

1、概述一般计算圆曲线可用坐标正算直接进行计算,具体思路和求解步骤,这里不再阐述。

若计算带有缓和曲线的圆曲线时,将测量中所学的支距法与坐标旋转、平移结合在一起,利用EXCEL表中强大的函数自动计算功能,准确快速的完成对缓和曲线的坐标计算。

比一般的手工计算快10~20倍,比CAD绘图计算快5~10倍。

并可以应用来指导工程施工、施工放样、审核图纸等工作。

2、计算过程分段在计算带有缓和曲线的圆曲线或圆曲线时,只要输入待求点的里程,程序将会自动会计算线路中桩的坐标、与中桩有一定夹角、距离的边桩坐标,与边桩中心线任意夹角的垂直桩基坐标。

若要计算其他的距离和夹角的坐标,相应的修改待求点里程、夹角和距离。

2.1、程序初始化:输入每个曲线所对应交点的半径、缓和曲线长、线路转角、连续三交点的里程和坐标、交点连线的坐标方位角,顺便计算出各个曲线要素以及曲线各主点的里程。

2。

2、初直线H Z i—1~ZH i段:(1)X ZHi—1和Y ZHi-1的计算X ZHi—1= X JDi—1+T i-1×cos(A i-1,i)Y ZHi—1= Y JDi—1+ T i-1×sin(A i-1,i)其中:T iJD i-1曲线的切线长;—1-—A i-1,i——JD i—1与JD i直线的坐标方位角;X JDi-1、Y JDi—1—-JD i-1的坐标;X ZHi—1、Y ZHi—1—-JD i-1对应的ZH点坐标。

(2)中桩计算公式:X中=L A×cos(A i—1,i)+ X ZHi—1Y中= L A×sin(A i-1,i)+ Y ZHi-1其中:L A—-待求点与ZH i的里程差;A i—1,i—-JD i-1与JD i直线的坐标方位角;X中、Y中——待求点里程的中桩坐标;其余符号同上。

利用ExcelVBA程序计算曲线坐标.kdh

利用ExcelVBA程序计算曲线坐标.kdh

利用Excel VBA 程序计算曲线坐标张卫庭高强(中铁济南勘察设计咨询院有限公司山东济南250022)【摘要】Excel 是大家很熟悉的办公软件,相信大家在工作中经常使用。

本文利用Excel VBA 编制程序计算铁路缓和曲线的坐标,并将计算坐标转化为任意平面坐标系下的坐标,最后转化全站仪可识别的文件格式,从而实现了曲线坐标计算和现场放样的一体化。

【关键词】Excel VBA ;曲线要素;坐标转化;全站仪1.曲线坐标计算1.1角度转化为弧度计算曲线坐标,难免要把角度转化为EXCEL所识别的弧度来计算,为此首先要做的是角度转化为弧度。

打开Excel ,按Alt+F11即进入VBAIDE ,在菜单上依次点击[插入]->[模块],然后输入如下代码Public Const pi =3.14159265359Public Function DEG(n As Double)Dim A As Double,B As Double,C As Double,D As Double,E As Double,F As Double,G As Double,KA As DoubleD =Abs(n)+0.000000000000001F =Sgn(n)A =Int(D)B =Int((D -A)*100)C =D -A -B /100DEG =F *(A +B /60+C /0.36)*pi /180End Function1.2曲线坐标计算1.2.1计算曲线坐标的公式加入缓和曲线后的铁路曲线示意图(见右图)坐标系的建立主要取决于控制点的情况。

如果控制点是为测设曲线而布设的,则坐标系一般采用ZH —XY 坐标系统①ZH~HY 段曲线点的坐标根据缓和曲线方程计算:x A =l A -l A 540R 2l 02y A =±l A 36Rl 0-l A7336R 3l 03!"#%%%%%$%%%%%&当曲线右偏时y A 坐标为正,左偏时y A 坐标为负。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

曲线坐标计算程序曲线坐标计算程序关键词:曲线坐标计算 EXCEL编程坐标曲线坐标实例摘要:利用EXCEL强大的函数功能通过曲线坐标计算的知识编制成曲线计算坐标的计算程序。

简单的输入曲线的里程桩号,通过坐标旋转、平移结合可以快速的计算完成与线路成任意角度的曲线上各中桩、边桩以及任意点坐标的计算。

1、概述一般计算圆曲线可用坐标正算直接进行计算,具体思路和求解步骤,这里不再阐述。

若计算带有缓和曲线的圆曲线时,将测量中所学的支距法与坐标旋转、平移结合在一起,利用EXCEL表中强大的函数自动计算功能,准确快速的完成对缓和曲线的坐标计算。

比一般的手工计算快10~20倍,比CAD绘图计算快5~10倍。

并可以应用来指导工程施工、施工放样、审核图纸等工作。

2、计算过程分段在计算带有缓和曲线的圆曲线或圆曲线时,只要输入待求点的里程,程序将会自动会计算线路中桩的坐标、与中桩有一定夹角、距离的边桩坐标,与边桩中心线任意夹角的垂直桩基坐标。

若要计算其他的距离和夹角的坐标,相应的修改待求点里程、夹角和距离。

2.1、程序初始化:输入每个曲线所对应交点的半径、缓和曲线长、线路转角、连续三交点的里程和坐标、交点连线的坐标方位角,顺便计算出各个曲线要素以及曲线各主点的里程。

2.2、初直线H Z i-1~ZH i段:(1)X ZHi-1和Y ZHi-1的计算X ZHi-1= X JDi-1+T i-1×cos(A i-1,i)Y ZHi-1= Y JDi-1+ T i-1×sin(A i-1,i)其中:T i-1——JD i-1曲线的切线长;A i-1,i——JD i-1与JD i直线的坐标方位角;X JDi-1、Y JDi-1——JD i-1的坐标;X ZHi-1、Y ZHi-1——JD i-1对应的ZH点坐标。

(2)中桩计算公式:X中=L A×cos(A i-1,i)+ X ZHi-1Y中= L A× sin(A i-1,i)+ Y ZHi-1其中:L A——待求点与ZH i的里程差;A i-1,i——JD i-1与JD i直线的坐标方位角;X中、Y中——待求点里程的中桩坐标;其余符号同上。

(3)边桩计算公式:X边=L A’×cosα’+ X中Y边= L A’×sinα’+ Y中其中:L A’——待求点与ZH i的里程差;α’——边桩中心线的坐标方位角;X中、Y中——待求点里程的中桩坐标;其余符号同上。

(4)垂直边桩计算公式X垂=L A’’×cosα”+ X边Y垂= L A’’×sinα” + Y边其中:L A——待求点与ZH i的里程差;α”——边桩中心线的坐标方位角;X中、Y中——待求点里程的中桩坐标;其余符号同上。

2.3、缓和曲线Z H i~HY i计算:以下符号的含义类同直线段,不再一一介绍。

(1)支距法计算相对坐标:x= L B- L B5/(40×R2×L02)y= L B3/(6×R×L0)- L B7/(336×R3×L03)(2)支距法转化为测量坐标系中的中桩坐标:X中= x ZHI+ x×cos (A i-1,i)+y×sin(A i-1,i)Y中= y ZHI+ x×cos (A i-1,i)-y×sin(A i-1,i)(3)边桩计算公式:X边= L B’×cos(α’)+ X中Y边= L B’×sin(α’)+ Y中(4)垂直边桩计算公式X垂= X边+ L B’’×cos(α”)Y垂= Y边+ L B’’×sin(α”)2.4、圆曲线HY i~YH i段:(1)支距法计算相对坐标:φc=(Kc-K HY)/R+β0x= R×sin(φc)y= R×(1-cosφc)+p(2)支距法转化为测量坐标系中的中桩坐标:X中= x ZHi+ x×cos (A i-1,i)+ y×sin(A i-1,i)Y中= y ZHi+ x×cos (A i-1,i)- y×sin(A i-1,i)(3)边桩计算公式:X边= L c’×cos(α’)+ X中Y边= L c’×sin(α’)+ Y中(4)垂直边桩计算公式:X垂= X边+ L c’’×cos(α”)Y垂= Y边+ L c’’×sin(α”)2.5、缓和曲线YH i~HZ i段:(1)支距法计算相对坐标:x= L D- L D5/(40×R2×L02)y= L D3/(6×R×L0)- L D7/(336×R3×L03)(2)支距法转化为测量坐标系中的中桩坐标:X中= x ZHI- x×cos (A i,i+1)+ y×sin(A i,i+1)Y中= y ZHI- x×cos (A i,i+1)- y×sin(A i,i+1)(3)边桩计算公式:X边= L D’×cos(α’)+ X中Y边= L D’×sin(α’)+ Y中(4)垂直边桩计算公式:X垂= X边+ L D’’×cos(α”)Y垂= Y边+ L D’’×sin(α”)2.6、HZ i~ZH i+1直线段:(1)中桩计算公式:X中=L E×cos(A i,i+1)+ X HZiY中= L E×sin(A i,i+1)+ Y ZHi(2)边桩计算公式:X边=L E’×cosα’+ X中Y边= L E’×sinα’+ Y中(3)垂直边桩计算公式X垂=L E’’×cosα”+ X边Y垂= L E’’×sinα” + Y边以上公式中的角度在EXCEL中输入时,在每个角度的后面输入:×Pi()/180。

4、工程实例:中铁七局集团郑州工程总公司驻马店至泌阳高速公路K86+961.22跨京广铁路分离式立交桥,线路曲线位于R=4000m的圆曲线上,线路左转偏角为αz=7°14′46.29″,交点里程(JD)为k87+084.664,交点坐标为(x JD=3640304.2268, y JD=522025.9457);直圆点(ZY)里程为k86+831.384,坐标为(x ZY=3640442.5089,y ZY=522238.1426);圆直点(YZ)里程为k87+337.264,坐标为(x YZ=3640140.2841, y YZ=521832.5588);与线路夹角为55°。

计算圆曲线上23墩(里程为k87+111.28)距线路中线为3.938m的桩基坐标,距边桩中线为5m的垂直边桩坐标。

(1)X YZi-1和Y YZi-1的计算缓和曲线为X ZHi-1,Y ZHi-1,而圆曲线时为的X YZi-1和Y YZi-1。

已知数据:JD i-1里程为K85+826.286 (X JDi-1 = 3640991.943 ,Y JDi-1 =523081.2612);A i-1,i = 236°54′32.28″;T i-1 = 310.734m X YZi-1= X JDi-1+T i-1×cos(A i-1,i)=3640991.943+310.734×cos236°54′32.28″= 3640822.291mY YZi-1= Y JDi-1+ T i-1× sin(A i-1,i)= 523081.2612+310.734×sin236°54′32.28″= 522820.9266m(2)X ZYi和Y ZYi的计算L=K ZYi-K YZi-1=k86831.384-k85+826.286=695.609mX ZYi= X YZi-1+ L×cos(A i-1,i)=3640822.291+695.609×cos236°54′32.28″= 3640442.509mY ZYi = Y YZi-1 + L×sin(A i-1,i)=522820.9266+695.609×sin236°54′32.28″=522238.1423m23号墩中桩、边桩和垂直边桩的计算:计算时,只需输入23号墩的里程: 87111.28(1)支距法计算相对坐标:L c = K23-K ZYi=k87+111.28-k86+831.384=279.896m φc=(Kc-K ZY)/R×180/3.1415927=(K87+111.28-K86+831.384)/4000×180/3.1415927 = 4°00′33.17″x= R×sin(φc)=4000×sin4°00′33.17″=279.667my= R×(1-cosφc)=4000×(1-cos4°00′33.17″)= 9.788m(2)支距法转化为测量坐标系中的中桩坐标:X中= x ZYi+ x×cos (A i-1,i) +y×sin(A i-1,i)=3640442.509+279.667×cos236°54′32.28″+9.788×sin236°54′32.28″=3640281.617mY中= y ZYi+ x×cos (A i-1,i) -y×sin(A i-1,i)= 522238.1423+279.667×cos236°54′32.28″-9.788×sin236°54′32.28″=522009.18m同理可计算出23号墩里程k87+111.28加2米(或者更大的,但值越小计算的精度越高)的里程桩号的坐标值,用来计算边桩和垂直边桩的坐标。

(3)边桩坐标计算:边桩计算公式中的角度关系在这里不再进行推导,在每个项目都根据实际情况推导就是。

如果不理解可参考一些测量书籍的相关曲线性质。

已知:L c’= 3.938 m α’=177°53′59.28″X边= L c’×cos(α’)+ X中=3.938×cos(177°53′59.28″)+3640281.617=3640277.682mY边= L c’×sin(α’)+ Y中=3.938×sin(177°53′59.28″)+522009.18=522009.324m(4)垂直边桩坐标计算:已知:L c’’=5m α”=87°53′59.10″X垂= X边+ L c’’×cos(α”)=3640277.682+5×cos87°53′59.10″= 3640277.865mY垂= Y边+ L c’’×sin(α”)=522009.324+5×sin87°53′59.10″= 522014.3206m以上是结合驿阳二标的现场实例,计算出了圆曲线上典型的一个断面的坐标。

相关文档
最新文档