9860線元法隧道測量(坐標正反計算、高程)計算程序

CASIO9860GSD線元法隧道三維〔坐標正反計算、高程〕計算程序1.A〔此為主程序〕Lbl0：〝1.LC=>XY〝：〝2.XY=>LC〝：〝3.ZHZL=>GC〝：〝PB=>V=1,2,3〞？→V：IfV=1：ThenGOTO1：IfEnd：IfV=2：ThenGOTO2：IfEnd：IfV=3：ThenGOTO3：ElseGOTO0：IfEnd：Lbl3：〝ZH=H〝？→H：〝SDZF=Z〝？→Z：Prog〝ZGCZCX〝：GOTO0：Lbl1：〝ZH=L〝？→L：IfL>173000AndL<174661.96：ThenGOTO4：ElseGOTO0：IfEnd：Lbl4：L→L：〝SDZF=Q〝？→Q：〝XLZJ,-Z+Y=Q〝：Q+0.125→Q：Prog〝ZBQXYS〝：〝JSJD=J〝：90→J▲Prog〝ZSZB〝：〝X=〝：X▲〝Y=〝：Y▲〝FWJ=O〝：O▼DMS▲L→H：Q-0.125→Z：Prog〝ZGCZCX〝：GOTO0：Lbl2XO=M〝?→M：〝YO=R〝?→R：173300→L：IfM>3845505.273AndM<3846506.099AndR>499371.832AndR<500352.224：ThenGOTO5：ElseGOTO2：IfEnd：Lbl5：0→Q：0→J：Prog〝ZBFS〝：〝LC=L〝：L▲〝JL=Q〝：Q▲〝SDZJ，-Z+Y=Q〝：Q-0.125→Q▲L→H：Q→ZProg〝ZGCZCX〝：GOTO02.正算坐標ZBZS(L-S)/4→H：90/π→F：HHF〔1/T-1/I〕/〔K-S〕→U：2HF/I→D：C+4D+16U→Ｏ：O+J→P：C+D+U→E：C+2D+4U→W：C+3D+9U→G：A+AbsH/3*(cosC+4(cosG+cosE)+2cosW+cosO)+QcosP→X：B+AbsH/3*(sinC+4(sinG+sinE)+2sinW+sinO)+QsinP→Y反算坐標：ZBFSLbl0：Prog〝ZBQXYSProg〝ZBZS〝：O-90→Z：(R-Y)cosZ-(M-X)sinZ→P：L+P→L：IfAbsP≥0.001：ThenGOTO0：ElseGOTO1：IfEnd：Lbl1：(R-Y)cosO-(M-X)sinO→Q曲線元要素數(shù)據(jù)庫：ZBQXYSIfL≥SAndL＜K*→TIfEnd：：Then**→S：**→A：**→B：**→C：**→I：**→K：*IfL≥SAndL＜K*→TIfEnd：：Then**→S：**→A：**→B：**→C：**→I：**→K：*IfL≥SAndL＜K*→TIfEnd：：Then**→S：**→A：**→B：**→C：**→I：**→K：*IfL≥SAndL＜K*→TIfEnd：：Then**→S：**→A：**→B：**→C：**→I：**→K：*IfL≥SAndL＜K*→TIfEnd：：Then**→S：**→A：**→B：**→C：**→I：**→K：*……………IfL≥SAndL＜K：Then**→S：**→*→TIfEnd┘A：**→B：**→C：**→I：**→K：*(注：如有多個曲線元要素連續(xù)添參與數(shù)據(jù)庫ZBQXYS中)5ZGCZCX(后有修改說明)Lbi1：〝SCGC=X〝？→X：〝R=M〝：5.98→M：〝CS=N〝：1.603→N：Prog＂GCIfEnd：IfH>B-TAndH<B+TThenH-B+T→K：GOTO2：IfEnd：IfH≥B+T：Then0→K：D→C：GOTO2：IfEnd：Lbi2：〝XLZG=G〝：A-〔B-H〕C-K^2÷2R→G▲IfZ≥0：ThenIf X>100Then〝YO1XGC=J〝：G+N→J▲〝YGCFSKD=J〝：Abs√(M^2-(X-〔G+N〕)^2)→J▲〝YKDCQ,+C,-Q=J〝：J-Abs(Z+0.000)→J▲〝YSBSJGC=J〝：G+N+√(M^2-(Z+0.000)^2)→J▲〝YGCCQ,C+,Q-=J〝：X-J→J▲Else〝YO1XGC=J〝：G+N→J▲〝YSBSJGC=J〝：G+N+√(M^2-(Z+0.000)^2)→J▲IfEnd：Else〝ZO1XGC=J〝：G+N→J▲IfX>100：Then〝ZGCFSKD=J〝：Abs√(M^2-(X-〔G+N〕)^2)→J▲〝ZKDCQ,+C,-Q=J〝：J-Abs(Z+0.000)→J▲〝ZSBSJGC=J〝：G+N+√(M^2-(Z+0.000)^2)→J▲〝ZGCCQ,C+,Q-=J〝：X-J→J▲Else〝ZSBSJGC=J〝：G+N+√(M^2-(Z+0.000)^2)→J▲IfEnd：IfEnd6高程計算主程序子程序：GCSJKIfH>起點樁號AndH≤Then第一豎曲線交點高程→A：第一豎曲線交點樁號→B：第一豎曲線前坡度→C：第一豎曲線后坡度→D：第一豎曲線半徑→R：IfEnd：IfH>AndH≤Then第一豎曲線交點高程→A：第一豎曲線交點樁號→B：第一豎曲線前坡度→C：第一豎曲線后坡度→D：第一豎曲線半徑→R：IfEnd：IfH>AndH≤Then第一豎曲線交點高程→A：第一豎曲線交點樁號→B：第一豎曲線前坡度→C：第一豎曲線后坡度→D：第一豎曲線半徑→R：IfEnd………(連續(xù)添加要素)說明：第一局部坐標局部〔1、2、3、4〕V=1進入坐標正算V=2進入坐標反算V=3進入單獨的高程計算V1、2、3時，則返回程序，要求再次輸入V值。變量說明：S……..起點里程A…..起點X坐標B….起點Y坐標C….起點方位角I….….起點半徑K…..終點里程 T.…終點半徑 L….計算點里程J…輸入左右〔左-，零，右+〕角度Q….中樁到邊樁的距離〝SDZJ-Z+Y〝為隧中支距，左〝-+〝坐標正算中輸入的Z值為線路支距坐標反算中得出的Z值也為線路支距留意：1JD,JL等于零。2例如：〔ZY,YZ,ZH,QZ,HY,YZ〕每個點都是起點所以特別留意3、規(guī)定以道路中線的前進方向〔即里程增大的方向〕區(qū)分左右；當線元往左偏時，=1045次代替。當所求點位于中線時，Q,J=0；當位于中線左側時，J取負值；當位于中線右側時，J取正值。當線元為直線時，其起點、止點的曲率半徑為無窮大，以1045次代替。當線元為圓曲線時，無論其起點、止點與什么線元相接，其曲率半徑均等于圓弧的半徑。當線元為完整緩和曲線時，起點與直線相接時，曲率半徑為無窮大，以1045次代替；與圓曲線相接時，曲率半徑等于圓曲線的半徑。止點與直線相接時，曲率半1045次代替；與圓曲線相接時，曲率半徑等于圓曲線的半徑。與圓曲線相接時，曲率半徑等于圓曲線的半徑。止點與直線相接時，曲率半徑等于設計規(guī)定的值；與圓曲線相接時，曲率半徑等于圓曲線的半徑。4、輸入與顯示說明V=1的時候JL(在中線上取零)J？邊樁時左右邊樁連線與線路中線的交角線元左右偏標志(左偏J=-1*角度，右偏J=角度，直線段J=0)角度X=××× 正算時，計算得出的所求點的X坐標Y=××× 正算時，計算得出的所求點的Y坐標FWJ=*** 計算得出所求點的方位角°′″V=2的時候XO=×××所求點的X坐標YO=××× 所求點的Y坐標LC=**** 求點所對應的里程JL=**** -5就是左邊相反就是右邊V=1看看坐標是否全都留意：驗算的時候不需要改動JL？J？[/center]◢說明：其次局部高程局部〔5高程計算主程序：SQXBG6高程計算主程序子程序：SJK〕1、A：交坡點高程2、B：交坡點樁號3、C：前坡度〔上坡為正，下坡為負〕4、D：坡度〔上坡為正，下坡為負〕5、R：豎曲線半徑6、H：待求點樁號7、E：下一豎曲起點樁號8、X〝SCGC〝測點實測高程9、M〝R〝O1點半徑〔此處如須計算開挖線加上二襯厚度、沉降量、其他構造層的厚度即可〕10、N〝CS〝線路中樁設計高程至O1點的距離11、J〝YO1XGC〝O1點〔圓心點〕高程12J〝YSBSJGC〝依據(jù)樁號及支距〔正反算程序中〕計算得到的測點〔設計〕高程J〝YGCFSKD〝反算中依據(jù)高程反算得到的支距J〝YKDCQ,+C,-Q〝實測高程計算得支距與坐標反算得支距進展比照,得到超欠挖J〝YGCCQ,C+,Q-〝實測高程與坐標反算得到高程進展比照,得到超欠挖13、J〝ZO1XGC〝O1點〔圓心點〕高程14、J〝YSBSJGC〝依據(jù)樁號及支距〔正反算程序中〕計算得到的測點〔設計〕高程J〝ZGCFSKD〝反算中依據(jù)高程反算得到的支距J〝ZKDCQ,+C,-Q〝實測高程計算得支距與坐標反算得支距進展比照,得到超欠挖J〝ZGCCQ,C+,Q-〝實測高程與坐標反算得到高程進展比照,得到超欠挖〔本人施測的是一大路隧道的上行線，線路中心線至隧道中心線的距離為0.125米，Z〝SDZJ-Z+Y〝為隧中支距，Z+0.125為線路支距。Z-0.125Z在左側取-+?！砈QXBG修改說明此式中X為實測的該點高程，V=13時，請將X0，當V=2時，輸入為實測的高程數(shù)據(jù)。一般在測量時，將拱部點測得三維后，反算得到里程和支距，即用實測高程計算得到該高程點的理論支距，再用理論支距與坐標反算的支距進展比照，超欠挖自明。同時用反算的支距也可計算出該點的理論高程，用實測的高程與計算的理論高程比照，自可得到拱部是高度是否超欠挖。C-D→F：Abs(RF÷2)→T：RAbsF÷F→R：IfH≤B-T：Then0→K：GOTO2：IfEnd：IfH>B-TAndH<B+TThenH-B+T→K：GOTO2：IfEnd：IfH≥B+T：Then0→K：D→C：GOTO2：IfEnd：Lbi2：〝XLZG=G〝：A-〔B-H〕C-K^2÷2R→G▲IfZ≥0：〝YO1XGC=J〝：G+N→J▲ThenIf X>100Then〝YGCFSKD=J〝：Abs√(M^2-(X-〔G+N〕)^2)→J▲〝YKDCQ,+C,-Q=J〝：J-Abs(Z+0.000)→J▲〝YSBSJGC=J〝：G+N+√(M^2-(Z+0.000)^2)→J▲〝YGCCQ,C+,Q-=J〝：X-J→J▲Else〝YSBSJGC=J〝：G+N+√(M^2-(Z+0.000)^2)→J：IfEnd：▲Else：〝ZO1XGC=J〝：G+N→J▲IfX>100：Then〝ZGCFSKD=J〝：Abs√(M^2-(X-〔G+N〕)^2)→J▲〝ZKDCQ,+C,-Q=J〝：J-Abs(Z+0.000)→J▲〝ZSBSJGC=J〝：G+N+√(M^2-(Z+0.000)^2)→J▲〝ZGCCQ,C+,Q-=J〝：X-J→J▲Else〝ZSBSJGC=J〝：G+N+√(M^2-(Z+0.000)^2)→J▲IfEnd：IfEnd4850中請將X〝X=〝▲Y〝Y=〝▲改為〝X=〝：X◢〝Y=〝：Y◢L〝LC=▲Q〝JL=〝▲改為〝LC=〝：L◢〝JL=〝：Q◢〝=A-〔B-H〕C-K2÷2R▲改為〝XLZG=〝：G=A-〔B-H〕C-K2÷2R▲將J〝YO1XGC〝=G+N▲J〝YO1XGC=〝：G+N▲將J〝YKDCQ,+C,-Q〝=J-Abs(Z+0.125)▲改為J〝YKDCQ,+C,-Q=〝：J-Abs(Z+0.125)▲將J〝YSBSJGC〝=G+N+√(M^2-(Z+0.125)^2)▲改為J〝YSBSJGC〝:G+N+√(M^2-(Z+0.125)^2)將J〝YGCCQ,C+,Q-〝=X-J▲改為J〝YGCCQ,C+,Q-〝:X-J左側同理。04030403程序說明1、程序中有樁號推斷，超出計算范圍，不會進展計算，要求你重輸入樁號。2、本程序中有X、Y的最大、最小計算推斷，超出X、Y的最大、最小值范圍，不會進展計算，要求你重X、Y值。3、本程序反算中有寬度方向和高度方向的超欠挖。4、算例局部見第一樓.提問者：lbs326人氣：1780賞分：+62查看樓主其它問題-回復內(nèi)容 只查看樓主與本人回復 按支持人數(shù)排序沙發(fā)：lbs326樓主04031、平曲線要素交點號JD37交點樁號SK174+133.90 交點坐標X=3845932.100交點坐標Y=499588.222右偏角45°54′41″ 前方位角297°18′17″ 半徑=1200 LS1=LS2=150ZH點樁號SK173+550.394 HY點樁號SK173+700.394YH點樁號SK174+511.96 樁號SK174+661.960A1=A2=424.264 T1=T2=583.596SK173+300 QD點坐標X=3845549.528Y=500329.288 C=297°18′17″SK173+550.394ZH點坐標X=3845664.390Y=500106.7933C=297°18′17″SK173+700.394HY點坐標X=3845735.948Y=499974.9917C=300°53′08.5″SK174+106.177QZ點坐標X=3845998.621Y=499668.235C=320°15′37.47″K174+511.960YH點樁號X=3846348.183Y=499465.9908C=339°38′06.45″ZBQXYSIfL≥173300AndL＜173550.394：Then173300→S：3845549.528→A：500329.288→B：297°18′17″→C：1E45→I：173550.394→K：1E45→T：IfEnd：IfL≥173550.394AndL＜173700.394：Then173550.394→S：3845664.390→A：500106.7933→B：297°18′17″→C：1E45→I：173700.394→K：1200→T：IfEnd：IfL≥173700.394AndL＜174106.177：Then173700.394→S：3845735.948→A：499974.9917→B：00°53′08.55″→C：1200→I：174106.177→K：1200→T：IfEnd：IfL≥174106.177AndL＜174511.960：Then174106.177→S：3845998.621→A：499668.235→B：320°15′37.47″→C：1200→I：174511.96→K：1200→T：IfEnd：IfL≥174511.96AndL＜174661.960：Then174511.960→S：3846348.183→A：499465.9908→B：339°38′06.45″→C：1200→I：174661.96→K：1E45→T：IfEnd注：1E45，意即無窮大。2、豎曲線要素第一變坡點樁號SK173+350 變坡點高程=1332.762 豎曲線半徑=24498.033前坡度-0.761%后坡度-2.431%豎曲線起點：SK173+145.441豎曲線止點：SK173+554.559其次變坡點樁號SK174+395.249 變坡點高程=1307.352 曲線半徑=56000前坡度-2.431%后坡度-1.690%豎曲線起點：SK174+187.949豎曲線止點：SK174+602.909第三變坡點的豎曲線起點樁號：SK174+900。GCSJKIfH>173.145.441AndH≤174187.949：Then1332.762→A：173350→B：0.00761→C：-0.02431→D：24498.033→R：IfEnd：IfH>174187.949AndH≤174900Then1307.352→A：174395.429→B：-0.02431→C：-0.0169→D：56000→R：IfEnd3、實例計算如下其中：X”SCGC”輸入任意實測高程值M”R”=5.98 N”CS”=1.602V=1 L輸入173806 JL輸入0.125SDZJ,-Z+Y后顯示0X=3845794.175Y=499886.936”FWJ=“305.928XLZG=1321.677YO1XG=1323.279J”ZSBSJGC”=1329.259V=1 L輸入173806 JL輸入-4.875SDZJ,-Z+Y后顯示-5X=3845790.126Y=499884.002 ”FWJ=“305.928XLZG=1321.677 YO1XG=1323.279J”ZSBSJGC”=1326.559V=1 L輸入173806 JL輸入5.125SDZJ,-Z+Y后顯示5X=3845798.224Y=499889.869 ”FWJ=“305.928XLZG=1321.677 YO1XG=1323.279J”ZSBSJGC”=1329.559J”ZGCFSKD”=Abs√(M^2-(X-〔G+N〕)^2)▲J”ZKDCQ,+C,-Q”=J-Abs(Z+0.125)▲J”ZSBSJGC”=G+N+√(M^2-(Z+0.125)^2)▲J”