单位：中铁二十局集团第二工程有限公司           作者：辜建君  
  1、主程序  TYQXJSDefm10：
Lbl 0:  “1.LC=>XY”: “2.XY=>LC”: {V}:V=1=>GOTO 1: ≠> GOTO 2 △
Lbl 1: {DZG}:D:Z”JL=”:G”J=”:Z=0=>G=0：△Prog “QXYS”;Prog “ZSZB” : ”X=”：X ▲ ”Y=”：Y▲"FWJ="：J→DMS▲
{HU}：H“XC”：U“YC”：I=0：J=0：Pol（X-H，Y-U）：J<0＝＞J=J+360：Δ“FW=”：J→DMS◢“S=”：I◢Prog “SQX” GOTO 0
Lbl 2: D=***:Z=0:G=0:{MI}:M”XO”:I”YO”: Prog “FSZB”:”D=”：D ▲”JL=”：Z▲Prog “SQX” GOTO 0    D=有效里程 即线路起点 里程
    2. 正算坐标ZSZB          
N=5:P＝（1\E-1\R）÷Abs（K-F）：Q＝Abs（D-F）÷N：S=90Q÷π： J＝C+(NPQ+2\R)NS：L=1←┘X＝A+Q÷6×(Cos C+Cos J +4∑（Cos (C+((L+0.5)PQ+2\R)×(L+.5)S),L,0,(N-1)）+2∑（Cos (C+((LPQ+2\R)LS,L,1,(N-1))）+ZCos（J+G)：←┘     Y＝B+Q÷6×(Sin C+Sin J +4∑（Sin (C+((L+0.5)PQ+2\R)×(L+.5)S),L,0,(N-1)）+2∑（Sin (C+((LPQ+2\R)LS,L,1,(N-1))）+Z Sin（J+G)： ←┘
    3. 反算坐标：FSZBLbl 0:Prog “QXYS” :Prog “ZSZB”:O=J-90:P= (I-Y)cosO-(M-X) sinO :D=D+P:AbsP≥0.001=> GOTO 0 :≠> GOTO 1 △ Lbl 1: Z= (I-Y)cosJ-(M-X) sinJ
    4. 《SQX》高程计算子程序
    Prog“SB”：W=Z[1]－Z[2]:W＞0=>Z[6]＝－1：≠>Z[6]＝1△Z[3]：T“T”＝Abs(Z[3]W/2) ：E“E”=T＾2÷2Z[3]：C=D-Z[4]：D≦Z[4]=>Z[7]=Z[1]：≠>Z[7]=Z[2]△Z[5]：AbsC≦T=> “H0”: Z[8]=Z[5]+CZ[7]+Z[6](T- Abs C)＾2÷2Z[3] ▲≠> “H0”: Z[8]=Z[5]+CZ[7] ▲
     5、曲线元要素数据库：QXYS
D≥F=>D＜K=> F=**:A=**:B=**:C=**:R=**:K=**E=**⊿⊿←┘
D≥F=>D＜K=>     F=**:A=**:B=**:C=**:R=**:K=**E=**⊿⊿←┘
D≥F=>D＜K=> F=**:A=**:B=**:C=**:R=**:K=**E=**⊿⊿←┘
D≥F=>D＜K=> F=**:A=**:B=**:C=**:R=**:K=**E=**⊿⊿←┘
D≥F=>D＜K=> F=**:A=**:B=**:C=**:R=**:K=**E=**⊿⊿←┘                     ……………………………D≥F=>D＜K=> F=**:A=**:B=**:C=**:R=**:K=**E=**⊿⊿←┘                
(注：如有多个曲线元要素继续添加入数据库QXYS中)
     6、《SB》竖曲线数据库
D≥竖曲线起点里程=>Z[1]=？：Z[2]=？：Z[3]=？：Z[4]=？：Z[5]=？△
D≥？=>Z[1]=？：Z[2]=？：Z[3]=？    ：…………△
     附注：  
     1、D所求的桩号：Z[1]=坡度1:Z[2]=坡度2：Z[3]=半径：Z[4]=交点桩号：Z[5]=变坡点控制高程：
     2、注意Z[1]、Z[2]的正负。        
     3、本程序可以计算上百公里的路线中桩高程。
     4、若在直线上时半径为10的45次方，竖曲线起点里程为直线上的起点里程，交点里程为直线中点里程。
    变量说明：
     F……..起点里程   A…..起点X坐标   B….起点Y坐标   C….起点方位角R….….起点半径   K…..终点里程     E.…终点半径     D….计算点里程G…输入左右（左-，零，右+）角度   Z….中桩到边桩的距离  
   注意：
     1:计算中桩坐标 J,JL 等于零。        
     2：这里的方位角从C是指 线元起点的方位角例如：（ZY , YZ , ZH , HY , YZ）每个点都是 起点
   所以特别注意
      1、规定      
     (1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，半径＝负半径；当线元往右偏时,半径=半径；当线元为直线时，半径=以10的45次代替。      
     (2) 当所求点位于中线时，Z,J=0；当位于中线左侧时，J取负值；当位于中线右侧时，J取正值。      
     (3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。    
     (4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。  
     (5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。      
     (6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。
     2、输入与显示说明      
     V=1的时候      JL ？  正算时所求点距中线的边距(在中线上取零)      J？ 边桩时左右边桩连线与线路前进中线的交角线元左右偏标志       (左偏J= -1*角度，右偏J=角度，直线段J=0) 角度         X=×××    正算时，计算得出的所求点的X坐标       Y=×××    正算时，计算得出的所求点的Y坐标      FWJ=***     计算得出所求点的切线方位角 °′″       FW=***    计算得出置镜点到测点方位角    H0=***    计算得出所求点的线路中线设计标高 V=2的时候      XO=×××  所求点的X坐标       YO=×××    所求点的Y坐标      D=****   求点所对应的里程      JL=****    求点到所对应的里程的距离  ：负就是左边，正就是右边
      1.      若果采用4800则把程序中Lbl 1: {DZG}:D:Z”JL=”:G”J=”:Z=0=>G=0：△Prog “QXYS”;Prog “ZSZB” : ”X=”：X ▲ ”Y=”：Y▲"FWJ="：J→DMS▲{HU}：H“XC”：U“YC”：I=0：J=0：Pol（X-H，Y-U）：J<0＝＞J=J+360：Δ“FW=”：J→DMS◢“S=”：I◢Prog “SQX” GOTO 0Lbl 2: D=***:Z=0:G=0:{MI}:M”XO”:I”YO”: Prog “FSZB”:”D=”：D ▲”JL=”：Z▲Prog “SQX” GOTO 04. 《SQX》高程计算子程序Prog“SB”：W=Z[1]－Z[2]:W＞0=>Z[6]＝－1：≠>Z[6]＝1△Z[3]：T“T”＝Abs(Z[3]W/2) ：E“E”=T＾2÷2Z[3]：C=D-Z[4]：D≦Z[4]=>Z[7]=Z[1]：≠>Z[7]=Z[2]△Z[5]：AbsC≦T=> “H0=”: Z[8]=Z[5]+CZ[7]+Z[6](T- Abs C)＾2÷2Z[3] ▲≠> “H0=”: Z[8]=Z[5]+CZ[7] ▲
      修改为Lbl 1: {DZG}:D:Z”JL=”:G”J=”:Z=0=>G=0：△Prog “QXYS”;Prog “ZSZB” :X ”X=”▲Y ”Y=”▲J"FWJ="▲{HU}：H“XC”：U“YC”：I=0：J=0：Pol（X-H，Y-U）：J<0＝＞J=J+360：ΔJ“FW=”◢I“S=”◢Prog “SQX” GOTO 0Lbl 2: D=***:Z=0:G=0:{MI}:M”XO”:I”YO”: Prog “FSZB”:D”D=”▲Z”JL=”▲Prog “SQX” GOTO 04. 《SQX》高程计算子程序Prog“SB”：W=Z[1]－Z[2]:W＞0=>Z[6]＝－1：≠>Z[6]＝1△Z[3]：T“T”＝Abs(Z[3]W/2) ：E“E”=T＾2÷2Z[3]：C=D-Z[4]：D≦Z[4]=>Z[7]=Z[1]：≠>Z[7]=Z[2]△Z[5]：AbsC≦T=> Z[8] “H0=” =Z[5]+CZ[7]+Z[6](T- Abs C)＾2÷2Z[3] ▲≠> Z[8] “H0=”=Z[5]+CZ[7] ▲



由本站论坛整理，未验证其正确性