道路曲线计算程序-4850

CZB(测站准备)I=Z[1]：J=Z[2]：X=Z[3]：Y=Z[4]Lb1 1:{VW}:V”YG”:W”HC”V>10=>GOTO 1△Lb1 2:{IJXY}:I”XC”:J”YC”:X”HSX”:Y”HSX”Z[1]=I:Z[2]=J:Z[3]=X:Z[4]=YPROG”FWJ”“A”:J▲“D”:I▲FWJ（方位角计算）POL（X-I，Y-J）：J〈0=>J=J+360FX ALL(全线放线)Lb1 1:Z=Z[5]:{Z}:Z”FX Z”:Z[5]=ZFixm:Z<Z[6] =>PROG”JD ALL”: ≠>Z>N =>PROG”JD ALL”△△Z[8]=0：PROG“JKCG”：E=Z[7]：{E}：E“Z-Y+”：Z[7]=EPROG“ZB”：PROG“FWJ”：“A”：J->DMS▲“D”:I▲Lb1 2:X=0:Y=1.3:{XY}:X”HV”:Y”JG”“CD HH”V+W+X-Y▲ZB”坐标计算“M=Z[17]：N=Z[19]：H=Z[21]：Z>N-H=>L=N-Z:T=Z[15]:K=H: ≠>L=Z-M:T=Z[13]: K=S: △K=0=>GOTO 1△A=(L^2/2RK)r（弧度转化为度）X=L-L^5/40R^2K^2Y=L^3/6RK-L^7/336R^3K^3Lb1 1:Z>M+S=>Z≤N-H=>A=（（L-K/2）/R）r：X=Q+R Sin A:Y=P+R(1-Cos A) △△L<0=>Y=0:A=0:X=L△POL(T-X,Y)Z>N-H=>J=180-J-Abs U:A=ＡbsU-A△X=B+Rec(-I,F-JO(字母))Y=C+JX=X+Rec(E,F+AO(字母)+90):Y=Y+J;I=Z[1]:J=Z[2]:”X”:X▲“Y”:Y▲PQXJS（平曲线计算）O(字母)=U/Abs U:U=Abs U:K=S:PROG”Q1”P=X:Q=Y:I=JK=H:PROG”Q1”Z[13]=(R+X-(R+P)Cos u)/Sin u+QZ[15]=(R+P-(R+X)Cos u)/Sin u+YL=S+H+(Uπ/180-I-J)RM=D-Z[13]：N=M+LJKCG（加宽超高计算）M=Z[17]：N=Z[19]：H=Z[21]：Z>N-H=>L=N-Z: ≠>L=Z-M △R>600=>I=-2(路弓横坡): ≠>I=2(最小超高值)：R〈390=>I=3……(一直到把所有要求输完) △△……(有几个=> 就有几个△)R>250=>X=0(加宽值): ≠>X=0.8：R〈200=>X=1……(一直到把所有要求输完) △△……(有几个=> 就有几个△)K=S:T=HK+T=0=>K=1:T=1:≠>KT≠0=>GOTO 1: ≠>K=1: ≠>T=1 △△△Lb1 1:Z>N-H=>Y=XL/T:J=(I+2)L/T-2:≠> Y=XL/K:J=(I+2)L/K-2△Y<0=>Y=0:J=-2△J>I=>Y=X:J=I△O(字母)<0=>”ZJK”:Y▲△“CG PD”:J▲Z[8] ≠0=> E=Z[7]：{E}：E“Z-Y+”：Z[7]=E:E=Abs E:I=0.01EJ:J<2=>J=-0.02E: ≠>J =-0.01EJ△O(字母)<0=>H=I:I=J:J=H△Z[7]>0=>”ZH”Z[25]+J▲≠>”HZ”:Z[25]+I▲JD ALL（平曲线）N=43433（起点桩号）I=23232.3232:J=4343433.433(起点坐标)X=22322.323:Y=434322.242(交点坐标)D=54435.899（交点桩号）U=-63(度)32…23“（左偏角）R=323.121S=35（缓和曲线长度1）:H=50（缓和曲线长度2）PROG“QX2”：Z〈N=>GOTO 2△X=23322.323:Y=424322.242(交点坐标)D=55435.899（交点桩号）U=33(度)32…23“（右偏角）R=323.121S=35（缓和曲线长度1）:H=50（缓和曲线长度2）PROG“QX2”：Z〈N=>GOTO 2△……Lb1 2PQXYS(平曲线要素)“A”：U-〉DMS▲{DRSH}：D“JD”：R：S“LS1”H“LS2”：Z[21]=H：PROG“PQXJS”“HZ“：Z[17]=M▲“HZ“：Z[19]=N▲Q1（平曲线计算子程序）X=K^2/24R-K^4/2384R^3Y=K/2-K^3/240R^2J=K/2RQ2(JD ALL子程序)Z[6]=N：B=X：C=Y：PROG“FWJ”：F=J：PROG“PQXJS”I=B：J=C：Z[23]=D：Z[1 7]=M：Z[19]=N：Z[21]=HSQX ALL（竖曲线）Z>49020(起点桩号) =>GOTO A △Lb1 1:D=45100(交点桩号):J=-6.289(第1个纵坡):H=528.88D=45330:I=-6.289(第2个纵坡):H=434.43:T（竖曲线半劲）=5000：PROG”SQ 2”:Z <N=>GOTO 2△D=46546:I=3.232:H=323:T=4000：PROG”SQ 2”:Z<N=>GOTO 2△……(到曲线资料完)Lb1 2SQ2（SQ ALL子程序）Z[26]=N:X=0.01(J-I):L=Abs（绝对值） TX/2:Z[14]=D-L:N=D+L:Z[16]=N:Z[18]=D: Z[20]=H:Z[22]=T:Z[24]=-X/Abs X:Z[28]=J:Z[30]=I:J=ISQJS(竖曲线计算)I=Z[30]:J=Z[28]:T=Z[22]:X=Z[24]:D=Z[18]:H=Z[20]:M=Z[14]:N=Z[16]Z《M=>X=0: ≠>Z≥N=>X=0△△Z>D=>L=H+0.01(Z-D)I+X(N-Z)^2/2T≠> L=H+0.01(D-Z)J+X(Z-M)^2/2TGCJS（全线高程计算）Z=Z[5]：{Z}：Z“JSZ”：Fixm:Z>Z[16] =>PRONG”SQX ALL”: ≠>Z〈Z[26] =>PRO NG”SQX ALL” △△Z[8]=1：PROG“SQXJS”：Z[25]=LZ<Z[6] =>PROG”JD ALL”: ≠>Z>N =>PROG”JD ALL”△△PROG“JKCG”。

曲线任意里程中边桩坐标正反算(CASIO fx-4850P计算器)程序一、程序功能本程序由一个主程序(1)和九个子程构成可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。

另外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲线任意里程中边桩坐标进行正反算。

本程序也可以在CASIO fx-4800P计算器计算器上运行。

特别申明：适用于弧长小于2倍半径的各种线元坐标正反算，精度优于1mm；二、源程序F1 (1)D“X1”:E“Y2”:F“X2”:G“Y2”:H“X3”:Y“Y3”:RL“LS”:X“JD”←V=D-F:W=E-G: Prog “C”:O=J:V=H-F:W=0-G:Prog“C”:U=J←M=O-U:M<0=>M=M+360△M<180 =>M=180-M:C=1: ≠>M=M-180:C=-1←P=L^2/(24R)-L^4/(2688R^3):Q=L/2-L^3/(240R^2): “T=”:T=Q+(P+R)tan(M/2)◢“S=”:S=RMπ /180+L：S◢ Prog“2”←F2 (2)Lb1 0: Z:Fixm:{KB}：V=K-X+T：W=Z：B←V＜S-L=>V≤0=>M=V+BcosZ:N=BsinZ:Goto 2△C=1=>W=180-Z△V≤L=> Prog“D”：≠> Prog“F”△Lb1 2：W=O：Prog“G”：Goto 0△V=S-V←V≤0=>M=V-BcosZ:N=-BsinZ◢≠>C≠1=>W=180-Z△Prog“D”:N=-N△W=U:Prog“G”:Goto 0←F3 (D)A=RL:M=V-V^5/(40A^2):N=V^3/(6A)-V^7/(336A^3):A=90V^2/πA：Prog“E”F4 (E)A=A-W:M=M-CBcosA:N=CN-BsinAF5 (F)A=(180V-90L)/( πR):M=RsinA+Q:N=P+R(1-cosA):Prog“E”F6 (G)A=cosW:W=sinW:V=F+TA-MA+NW:W=G+TW-MW-NA: “X0=”:V◢“Y0=”:W◢Po1(V-0,W-0):J≤0=>J=J+360△”J=”:J◢”I=”:I◢F7 (C)Po1(V,W):J＜0=>J=J+360△F8 (ZB)]A“X0”:B“Y0”DQ←A+DcosQ: “X”◢B+DsinQ: “Y”F9 (AB)A“X0”:B“Y0”XY ←L=Po1(V-A,Y-B):L◢Q=W←Q≥0=>Q◢≠>Q=Q+360:Q◢F10 (00)QRZ“LS”:X“JD”←P=Z^2/(24R) ←V=Z/2-Z^3/(240R^2) ←A=(P+R)tan(Q/2)+V: “TH”◢B=RQπ/180+Z: “LH”◢(R+P)/cos(Q/2)-R: “EH”◢2A-B: “DH”◢C=J-A: “ZH”◢C+Z: “HY”◢C+B/2: “QZ”◢C+B-Z: “YH”◢C+B: “HZ”◢。

1．中边桩计算程序<1>主程序名—XYLbl 1:{KM}：K：M：Prog”3”:X=C+Mcos(W-90:Y=D+Msin(W-90:”X=”:X:Pause 0:”Y=”:Y▲Goto1<2>子程序—3PrOg”QX”：Prog”20”<3>曲线要素—QX数据实例：某高速公路中线设计如下图所示K<2000=>B:E:F:O:P:Q:R: K≤4570=>B=25:E=1000:F=1000:O=100P=66:Q=4000:R=800: K≤5040=>B=-17:E=1086.4378:F=1555.1451:O=90:P=91:Q=4570:R=800: K≤5670=>B=35:E=1168.2556:F=2016.0909:O=120:P=74:Q=5040:R=1000: B=-25:E=1180.9927:F=2636.2334:O=150:P=109:Q=5670:R=1500注：无缓和曲线时，O=0.00001当输入里程小于2000时，为新曲线要素输入，可修改值域范围<4> 计算过程—20Z=O：Prog”8”：A=B/AbsB:L=ABπR/180+O:T=(D+RcosN)tan(AB/2)+C-RsinN：K≤Q-T=>C=(K-Q)cosP+E：D=(K-Q)sinP+F:W=P： K≤Q-T+O=>Z=K-Q+T：Prog”9”：K≤Q-T+L-O=>Z=O: Prog”9”:V=90(K-Q+T-O)/ πR:S=2R sinV:N=W+AV:C= C+S cosN:D=D+S sinN:W=W+2AV: Z=Q-T+L-K: Prog”8”:V=P+B-180-AV:W=P+B-AN :C=S cosV+E+T cos(P+B:D=S sinV+F+T sin(P+B: Z=K-Q+2T-L:W=P+B:C=E+ZcosW:D=F+Z sinW△△△△N=Z[3]<5> 缓和曲线方程式—8V=RO:C=Z-Z^5/40V2+Z^9/3456V^4:D=Z^3/6V-Z^7/336V^3+Z^11/42240V^5:N= 90Z^2/πV:V=tan-1(D/C:S=(D^2+C^2)^0.5<6> 子程序—9Prog”8”:C=Scos(P+AV)+E-TcosP:D=Ssin(P+AV)+F-TsinP:W=P+AN2.桩号、边桩距返算程序—KM方式一：试算法（速度慢，程序简单）Fixm:K=Q-T+L-O:Prog”20”Lbl 1:{XY}Lbl 2:POL(X-C,Y-D:K=Isin(J-W+90)+K:Abs sin(J-W+90)<0.0001＝>”K=”:K:Pause0:”M=”:M=Icos(J-W+90▲≠>Prog”20”:Goto2△方式二：解算法（程序复杂，计算速度快）Defm5: Lbl 1:Fixm:{XY}:Pol(X-E,Y-F:Z[1]=I:Z[2]=J:cos(J-P<0＝>K=Q+Icos(J-P: ≠>K=Q+Icos(J-P-B)+2T+L△K≤Q-T＝>M=Icos(J-P+90:Goto5≠>Goto2△△Lbl 2:Z=O:Porg”8”:U=(D+RcosN)+cos(B/2):Z=P+B/2+90A:C=UcosZ+E:D=UsinZ+F:Pol(X-C,Y-D:Asin(J-P-AN+90 A<0＝>K=K-AZ[1]cos(Z[2]-P+90)tan(90 (K-Q+T)^2/πRO):Goto3:≠>Asin(J-P-B+AN+90A>0 ＝>K=K-AZ[1]cos(Z[2]-P-B-90 )tan(90(K-Q+T-L)^2/πRO):Goto4:≠>K=Q-T+O+Asin-1(sin(J-P-AN+90 A))×πR/180 :M=A(I-R:Goto5△△Lbl 3:Z=K-Q+T:Prog”9”Pol(X-C,Y-D:K=K+Isin(J-W+90 :Abs sin(J-W+90<.001＝>M=Icos(J-W+90 :Goto5: ≠>Goto3△Lbl 4:Z=Q-T+L-K:Prog”8”:V=P+B-180-AV:W=P+B-AN:C=cosV+E+Tcos(P+B:Pol(X-C,Y-D:K=K+Isin(J-W+90 :Abs sin(J-W+90 <.001＝>M=Icos(J-W+90:Goto5: ≠>Goto4△Lbl 5:”K=”:K:Pause0:”M=”:M▲注：程序“20”为中桩计算程序Q—直缓点里程；T—切线长，L—曲线长；O—缓和曲线长；W—切线方位角Pol—极坐标转换，为涵数特别注意数字0和字母O的区别一、程序功能本程序由一个主程序(TYQXjs)和两个子程——正算子程序(SUB1)、反算子程序(SUB2)序构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。

FX-4850P公路放样程序单位:云南第二公路桥梁工程有限公司作者:杨艳刚文件名ZBJSDefm 12 :K“X0”:V“Y0”:L“Ls”: Z【8】=180÷πR ↙Q=0.5L-L∧3÷240R2:T=（R+ L2÷24R）tan.5A+ Q↙“T=”:T▲B= Z【8】L÷2: Z【10】“QX”=A÷Z【8】+ L↙“QX=”: Z【10】▲E=(R+ L2÷24R）÷cos.5 A-R↙“E=”:E▲E“ZH”=○“JD”-T↙“ZH=”: E▲Z【1】“HY”=E+ L↙“HY=”: Z【1】▲Z【9】“QZ”= E+.5 Z【10】↙“QZ=”: Z【9】▲Z【2】“YH”= E+Z【10】- L↙“YH=”: Z【2】▲Z【3】“HZ” = E+Z【10】↙“HZ=”: Z【3】▲P“1，2”：u “X(JD)”: Z“Y(JD)”↙Lb1 0: ｛m｝:M“CD”M≤△M≤Z【1】△↙M＜Z【2】△↙M≥Z【3】△↙M≥Z【2△↙Lb1 1:G=M-E:D=√((G-G∧5÷40R2 L2)2+( G∧3÷6RL- G∧5÷336 R∧3L)2) ↙P＜【8】G2÷6L:H=F+BG2÷L2: 【8】G2÷6L:H=F-BG2÷L2△Goto3↙Lb1 2:G=M- Z【1】:W=R(1-Cos(B+ Z【8】G))+ L2÷24R：D=√（W2+（R Sin(B+ Z【8】G)+Q）2）↙P＜tan-1(W÷（R Sin(B+ Z【8】G)+ Q)）:H=F+B+ Z【8】G: ≠F- tan-1(W÷（R Sin(B+ Z【8】G)+ Q)）:H= F-B- Z【8】G△↙Lb1 3:X=U+T Cos(F-180)+D Cos C↙“X=”: X▲Y=Z+T Sin (F-180)+D Sin C↙“Y=”: Y▲Prog “FY”↙Lb1 4: Z【6】=X: Z【7】=Y:｛S｝:X“ZX”= Z【6】+S Cos(H-90)“ZX=”:X ▲Y“ZY”= Z【7】+S Sin (H-90) ↙“ZY=”:Y▲Prog “FY”: ｛N｝:X“YX”= Z【6】+N Cos(H+90) ↙“YX=”:X▲Y“YY”= Z【7】+N Sin (H+90) ↙“YY=”: Y ▲Prog “FY”: Goto 0Lb1 5:G= Z【3】-M: Z【4】=(T-G+ G∧5÷40 R2 L2) SinA+( G∧3÷6RL- G∧5÷336 R∧3L) CosA: Z【5】= T+( T -G+ G∧5÷40 R2 L2) CosA-( G∧3÷6RL- G∧5÷336 R∧3L) SinA: ↙D=√(Z【4】2+ Z【5】2): P＜-1(Z【4】÷Z【5】)：H=F+A-BG2÷L2: ≠-1(Z【4】÷Z【5】)：H=F-A+BG2÷L2△Goto3↙Lb1 6:D=○-M: C= F+180:H=F: Goto8↙Lb1 7: D=T+M- Z【3】: P＜≠△H=C↙Lb1 8:X=U+D CosC“X=”:X▲Y=Z+ D SinC↙“Y=”:Y▲Prog “FY”: Goto 4↙文件名FYI=0:J=O:PoL(X-K,Y-V): Z【11】“D”=I↙“D=”: Z【11】▲Y-V≥【12】“F”=J↙≠【12】“F”=360+J↙“F=”： Z【12】▲使用程序时依次输入：1、置仪点坐标：X0？ Y0？（分别为X坐标Y坐标）2、缓和曲线长度Ls?3、曲线半径R？4、导线转角A？5、交点里程桩号JD？6、导线转向1,2？（右转输1，左转输2）7、交点坐标X（JD）?Y(JD)?8、待测点桩号CD?9、导线方位角F？10、所测中桩的左右距离S？N?（S-左边距 N-右边距）注：文件ZBJS 和FY都须输入才能运行D=-待放样点的距离F=-待放样点的角度在FY中运行，为坐标放样X？为待放点的X坐标 Y？为待放点的Y坐标K？为置仪点的X坐标 V？为置仪点的Y坐标D=-置仪点到待放点的距离F=-置仪点到待放点的角度。