道路的CASIO 5800P 万能曲线 计算程序
- 格式:doc
- 大小:52.00 KB
- 文档页数:5
道路的CASIO 5800P 万能曲线 计算程序 曲线任意里程中边桩坐标正反算 (CASIO fx-5800P计算器)程序 程序功能 本程序由一个主程序(ZHUCHENXU)和几个子程序——正算子程序(SUB1)、反算子程序(SUB2)、数据库(SUB3)、算方位角程序(FA)、使用方位角算坐标(NE)、边坡放样程序(BIANPO)、隧道圆心放样程序(SDYX)构成,可以根据直线、圆曲线、缓和曲线(完整或非完整型)的线元要素(起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径)及里程边距或坐标,对该曲线段范围内任意里程中边桩坐标进行正反算。 源程序 增加变量20→DimZ 1.主程序(AZFS) (P-R)÷(2HPR) →D:180÷π→E:"1.SZ → XY,2.XY → SZ,3.SDYX,4. BIANPO "?N: Lbl 1:”XYLC=”?S: Prog "SUB3": N≥2=>Goto 2: Abs(S-O)→W:”PYJ=”?→Z:Prog "SUB1": "XS=”:X→X◢"YS=”:Y→Y◢"FS=”:(F-90) →F:F►DMS◢Prog "FA":Goto 1: Lbl 2:”CX=”?X: X→I:”CY=”?Y: Y→J:”DMG=”?→Z[7]:Prog "SUB2": "S=":(O+W) →S◢"Z=":Z→Z◢If N=3:Then Prog”SDYX”: IfEnd :If N=4:Then Prog”BIANPO”:IfEnd:Goto 2 2. 正算子程序(SUB1) 0.1739274226 →Z[1]: 0.3260725774→B: 0.0694318442→K: 0.3300094782→L: (1-L) →F: (1-K)→M: U+W(Z[1]cos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QEFW (C+FWD))+Z[1]cos(G+QEMW(C+MWD))) →X: V+W(Z[1]sin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD)) +Bsin(G+QEFW(C+FWD))+Z[1]sin(G+QEMW(C+MWD))) →Y: G+QEW(C+WD)+90→F: X+ZcosF→X: Y+ZsinF→Y
3. 反算子程序(SUB2) G-90 →Z[2]: Abs((Y-V)cosZ[2]-(X-U)sinZ[2]) →W: 0→Z:Lbl 0:Prog "SUB1": Z[2]+QEW(C+WD) →L: (J-Y)cosL-(I-X)sinL→Z: If Abs(Z)<1E-6:Then Goto1:Else W+Z→W:Goto 0: IfEnd Lbl 1: 0→Z:Prog "SUB1": (J-Y)÷sinF→Z 4. 数据库(SUB3) If S<226100:(曲线起点里程)Then Cls:Stop: IfEnd If S<226255.833:(圆曲线起点里程)Then 226100→O: 99037.9736→U: 4282.3590→V: 201.3266304→G: 1045→P: 1045→R: 253.543→H: 0→Q: Goto 1:IfEnd …… …… If S>236005.651(线路终点里程桩号):Then Cls:Stop: IfEnd Lbl 1: 1÷P→C: 4、SUB3数据库说明: S待求桩号,O=曲线起点桩号,U曲线起点X坐标,V曲线起点Y坐标,G曲线起点切线方位角, P曲线起点半径, R曲线终点半径, H=曲线长度,Q转向(直线0,左转-1,右转1)
5. 算方位角程序(FA) LbI0:“XC=”? A : “YC=”? T :X:Y 0→J: POI(X-A,Y-T)→I:“BC=”:I◢If J<0: Then 360+J→J:ElseJ→J :IfEnd “FW=”:J►DMS◢ Goto0 6. 使用方位角算坐标(NE) LbI0:“X=”? X: “Y=”? Y: “L=”? L:“F0”?F: L×Cos(F)+X→N: “N=”N◢L×Si
n(F)+Y→E: “E=”E◢ Goto0 7. 边坡放样程序(BIANPO) “QPDJL=”? → Z[3]: “QPDGD=”? → Z[4]: “PG=” ? → Z[5]:2→Z[6]: “P1=” ? → Z[8]: “P2=” ? → Z[9]: “P3=” ? → Z[10] (没有的话就删): Z[7]-Z[4] → Z[11]:If Z[11]〈0: Then AbsZ[11] →Z[11]: “T=”:Z[11] ◢Else“W=”:Z[11] ◢IfEnd Z[11]÷Z[5] →Z[12]:If Z[12]〈1 Then Z[8]×Z[11] →Z[13]:Else Z[8]×Z[5]+Z[9]×Z[5]×(Z[12]-1)→Z[13]:IfEnd: AbsZ-(Z[13]+Z[6] ×Int Z[12]+Z[3])→Z[14]:“BP+C-J=”: Z[14] ◢ 说明: Z[3]= QPDJL ……………… 起坡点离中桩距离 Z[4]=“QPDGD” ……………起坡点设计标高 Z[5]=“PG” ……………… 坡高 Z[6]=“PTK” ……………… 平台宽度 Z[7]=“DMG” ………………实测地面高度 Z[8]=“P1” ……………… 一级边坡坡度 Z[9]=“P2” ……………… 二级边坡坡度 Z[10]=“P3” ……………… 三级边坡坡度(没有的话就删) Z[14]=“BP+C-J=” ………………边坡超欠挖宽度(+表超,-表欠) 8. 隧道圆心放样程序(SDYX) “RS=”? →Z[17]:”RZ=” ?→Z[18]: “RX=”? →Z[19]: “RW=”? →Z[20]: (如要储存:把(”RS=”?)换成圆的半径就OK了) IF S<137778.025:THEN 259.468-0.00725×(S-136000) +1.515 →Z[10]:“NGC=”:Z[10] ◢ELSE 248.092-0.00725×(S-137778.025)- (S-137778.025)^2/90000→Z[10]:“YLJ=”:Z[10] ◢IfEnd Z[10]- Z[7] →Z[15]:”GC=”: Z[15] ◢ IF Z[15]≥4.7:√(Z[17]^2-ABS(Z-4.33)^2)- Z[15]→Z[16]: “Z[16]1=”:Z[16] ◢IfEnd IF0≤ Z[15]AND Z[15]<4.7:√(Z[18]^2- Z[15] ^2)-ABS(Z-4.33)→Z[16]: “Z[16] =”:Z[16] ◢IfEnd IF Z[15]< 0AND Z[15]>-2.5:√(Z[19]^2- Z[15] ^2)-ABS(Z-4.33)-2.5→Z[16]: “Z[16] =”:Z[16] ◢IfEnd IF Z[15]<-2.5:√(Z[20]^2-ABS(Z-4.33) ^2)-11.8+ Z[15]→Z[16]: “Z[16]2 =”:Z[1 6] ◢IfEnd 以上为插入的开挖程序(每一条隧道的数据不一样都要休改)
IF S<229000:THEN 243.997+0.0205×(S-226520)-(S-227318)^2/2/16000+1.503 →Z[10]:“NGC=”:Z[10] ◢ELSE 286.385-0.023×(S-230500)- (S-230593.753)^2/2/165000+1.503→Z[10]:“YLJ=”:Z[10] ◢IfEnd Z[10]- Z[7] →Z[15]:“+SG,-JD”:Z[15] ◢√(Z[15]^2+(ABS(Z)-6.8)^2)→Z[16]: “R=”:Z[16] ◢ “XYLC=”:输入所求桩号所在线元的任意桩号 CX、CY:输入实测坐标。DMG输入实测高程 S、C:所测点的实际计算里程、偏距(左-,右+) RS、RZ、RX、RW:输入圆的设计半径. NGC:为程序自算的设计高程,不用算。 GC:为设计高程与实测高程的差值 ,(正值为挖,负值为填) “Z[16]R=”:隧道超欠挖值(正为超,负为欠) 使用说明 1、规定 (1) 以道路中线的前进方向(即里程增大的方向)区分左右;当线元往左偏时, Q=-1;当线元往右偏时,Q=1;当线元为直线时,Q=0。 (2) 当所求点位于中线时,Z=0;当位于中线左铡时,Z取负值;当位于中线中线右侧时,Z取正值。 (3) 当线元为直线时,其起点、止点的曲率半径为无穷大,以10的45次代替。 (4) 当线元为圆曲线时,无论其起点、止点与什么线元相接,其曲率半径均等于圆弧的半径。 (5) 当线元为完整缓和曲线时,起点与直线相接时,曲率半径为无穷大,以10的45次代替;与圆曲线相接时,曲率半径等于圆曲线的半径。止点与直线相接时,曲率半径为无穷大,以10的45次代替;与圆曲线相接时,曲率半径等于圆曲线的半径。 (6) 当线元为非完整缓和曲线时,起点与直线相接时,曲率半径等于设计规定的值;与圆曲线相接时,曲率半径等于圆曲线的半径。止点与直线相接时,曲率半径等于设计规定的值;与圆曲线相接时,曲率半径等于圆曲线的半径。