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曲线元要素数据库：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：*(注：如有多个曲线元要素连续添参与数据库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………(连续添加要素)说明：第一局部坐标局部〔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、规定以道路中线的前进方向〔即里程增大的方向〕区分左右；当线元往左偏时，=1045次代替。当所求点位于中线时，Q,J=0；当位于中线左侧时，J取负值；当位于中线右侧时，J取正值。当线元为直线时，其起点、止点的曲率半径为无穷大，以1045次代替。当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以1045次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半1045次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。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〝YGCFSKD〝反算中依据高程反算得到的支距J〝YKDCQ,+C,-Q〝实测高程计算得支距与坐标反算得支距进展比照,得到超欠挖J〝YGCCQ,C+,Q-〝实测高程与坐标反算得到高程进展比照,得到超欠挖13、J〝ZO1XGC〝O1点〔圆心点〕高程14、J〝YSBSJGC〝依据桩号及支距〔正反算程序中〕计算得到的测点〔设计〕高程J〝ZGCFSKD〝反算中依据高程反算得到的支距J〝ZKDCQ,+C,-Q〝实测高程计算得支距与坐标反算得支距进展比照,得到超欠挖J〝ZGCCQ,C+,Q-〝实测高程与坐标反算得到高程进展比照,得到超欠挖〔本人施测的是一大路隧道的上行线，线路中心线至隧道中心线的距离为0.125米，Z〝SDZJ-Z+Y〝为隧中支距，Z+0.125为线路支距。Z-0.125Z在左侧取-+。〕SQXBG修改说明此式中X为实测的该点高程，V=13时，请将X0，当V=2时，输入为实测的高程数据。一般在测量时，将拱部点测得三维后，反算得到里程和支距，即用实测高程计算得到该高程点的理论支距，再用理论支距与坐标反算的支距进展比照，超欠挖自明。同时用反算的支距也可计算出该点的理论高程，用实测的高程与计算的理论高程比照，自可得到拱部是高度是否超欠挖。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查看楼主其它问题-回复内容 只查看楼主与本人回复 按支持人数排序沙发：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”