新曲线任意里程中边桩坐标反正算.doc

曲线任意里程中边桩坐标反正算（CASIO fx-5800P计算器）程序一， 程序功能 本程序由一个主程序（ZHUCHENXU）和几个子程序-正算子程序（SUB1）反算子程序（SBU2）数据库（SBU3）,算方位角程序（FA），使用方位角算坐标（NE），边坡放样程序（BIANPO）,隧道圆心放样程序（SDYX）构成，可以根据直线，圆曲线，缓和曲线（完整或非完整型）的线元要素（起点坐标，起点里程，起点切线方位角，线元长度，起点曲率半经，止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩进行正反算。二， 源程序增加边量22DimZ1,主程序（ZHUCHENXU） 1.SZXY , 2.XYSZ , 3.SDYX, 4BIANPO? N: Lb1 1: ZHUANGHAO= ?S: Prog SUB3 : N2=Goto 2: Abs(S-0) W: BIANZHU=? Z:ProgSUB1:XS=:XXYS=:YYFS=:(F-90) F:FDMSProgFA:Goto 1 :Lb1 2: CX=?X: XI:CY=?Y:YJ: DMG=? Z7:ProgSUB2:S=:(O+W) S Z:ZZIf N=3:Then ProgSDYX:IfEnd:IfN=4:Then ProgBIANPO”:IfEnd:Goto 2 2，正算子程序（SUB1） 0.1739274226Z1: 0.3260725774B: 0.0694318442K:0.3300094782L: (1-L) F: (1-K) M:U+W(Z1cos(G+QEKW(C+KWD)+Bcos(G+QELW(C+LWD)+Bcos(G+QEFW(C+FWD)+Z1cos(G+QEMW(C+MWD) X:V+W(Z1sin(G+QEKW(C+KWD)+Bsin(G+QELW(C+LWD)+Bsin(G+QEFW(C+FWD)+Z1sin(G+QEMW(C+MWD) )Y:G+QEW(C+WD)+90F: X+ZcosFX: Y+ZsinFY 3,反算子程序（SUB2） G-90Z2: Abs(Y-V)cosZ2-(X-U)sinZ2) W: 0Z:Lb1 0:ProgSUB1: Z2+QEW(C+WD) L: (J-Y)cosL-(I-X)sinLZ:If Abs(z)1E-6:Then Goto1:Else W+ZW:Goto 0:IfEndLb1 1:0Z: ProgSUB1: (J-Y)sinFZ 4,数据库（SUB3）If S226100:Then Cls:Stop:IfEndIf S226255.833:Then 2261000: 99037.9736U: 4282.3590V: 201.3266304G: 1045P: 1045R: 253.543H: 0Q:Goto1:IfEnd.If S236005.651:ThenCls:IfEndLb1 1:1PC: (P-R) (2HPR) D:180E: 5算方位角程序（FA） LbI 0: XC=?A: YC=? T:X:Y0J:POI(X-A,Y-T) I: BC=:IIf J0:Then360+JJ:Else JJ:IfEndFW=:JDMS Goto0 6.使用方位角算坐标（NE）Lbl 0： X=? X: Y=? Y: L=? L: F0 ?F: LCos(F)+XN:N=NLSin(F)+YE: E=EGoto0 7.边坡放样程序（BIANPO）ProgTCD:ProgBHD:IF S226255.731 :then 232.257Z17:225844.269Z18: - 0.69Z19:11789.745 Z22:Goto 3:IfEndIF S226351.25:then 230.837Z17:226050Z18:2.8Z19:1045Z22:Goto 3:IfEndIF S226688.75:then 239.272 Z17:226351.25Z18:2.8Z19:- 45000Z22:Goto3:IfEndIF S227318:then 243.997Z17:226520Z18:2.05Z19:1045Z22:Goto3: IfEndIF S227982:then 260.356Z17:227318Z18:2.05Z19:-16000Z22:Goto3:IfEndLb1 3:Z17+(S-Z18) Z19/100+(S-Z182/2/Z22-(13+Z20) 0.02Z4: QPDGD=:Z4 Z7-Z4 Z11:If Z110 Then Abs(Z11) Z11: T=:Z11 13+Z20+Z21Z3:8Z5:0Z6:1.5Z8: 1.75Z9:ElseW=:Z11 15.6+Z20+Z21 Z3:10Z5:2Z6:1.0Z8:1.25Z9:IfEndZ11 Z5 Z12:If Z121 Then Z8Z11 Z13:Else Z8 Z5+Z9Z5 (Z12-1) Z13:IfEnd: AbsZ-(Z13+Z6Int Z12+Z3) Z14: BP+C-J=:Z14 说明： 只适用二级边坡Z14= BP+C-J.边坡超欠挖宽度（+表超，-表欠） 8，隧道圆心放样程序（SDYX）IF S229000:THEN 243.997+0.0205(S-226520)-(S-227318)2/2/16000+1.503Z10: NGC=:Z10 ELSE 286.385-0.023(S-230500)-(S-230593.753)2/2/165000+1.503Z10: Y L J= :Z10 IfEndZ10-Z7 Z15: +SG,-JD : Z15 (Z152+(ABS(Z)-6.8)2) Z16: R= :Z16 9,TCD(停车带)IF S227020: then 0Z20 :reture: IfEndIF S227035:then(S-227020) 3.5/15Z20 :reture: IfEndIF S227065: then 3.5Z20 : reture: IfEndIF S227080: then(S-227065) 3.5/15Z20 :Else 0Z20 : IfEnd 10,BHD(变化段)IF S227300: then 0Z21 :reture: IfEndIF S227450: then(S-227300) 1.85/300Z21 :Else 0.925Z21 : IfEnd三，使用说明1，规定（1） 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1; 当线元往右偏时，Q=1； 当线元为直线时，Q=0。（2） 当所求点位于中线时，Z=0; 当位于中线左侧时，Z取负值；当位于中线右侧时，Z取正值 （3）当线元为直线时，其起点,止点的曲线半径为无穷大，以10的45次代替。 （4）当线元为圆曲线时，无论其起点,止点于什么线元相接，其曲率半经均等于圆弧的半径。 （5）当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。 （6）当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。 2，输入于显示说明 SZXY, 2. XYSZ, 3. SDYX, 4. BIANPO N? 选择计算机方式，输入1表示进行由里程,边距计算坐标；输入2表示由做坐标反算里程和边距。输入3表示隧道圆心放样和计算实测半径；输入4表示进行路基边坡放样。 3，盖板涵斜交操作说明： （以左前点为列）先把盖板涵斜交角度减90后=A再盖板涵左长度*C o s A=B, 盖板涵左长度*SinA=C,盖板涵身宽度/