CASIO 9860G SD线元法隧道三维（坐标正反计算、高程）计算程序.doc

CASIO 9860G SD线元法隧道三维（坐标正反计算、高程）计算程序悬赏分：+62 作者：刘工2010-4-31. A（此为主程序）Lbl 0：1.LC=XY：2.XY=LC：3.ZHZL=GC：PB=V=1,2,3？V：If V=1：Then GOTO1 ：IfEnd ：If V=2：Then GOTO 2 ：IfEnd ：If V=3：Then GOTO3 ：Else GOTO 0 ：IfEnd： Lbl 3：ZH=H？H ：SDZF=Z？Z： ProgZGCZCX： GOTO 0：Lbl1 ：ZH=L？L：IfL173000AndL3845505.273AndM499371.832AndRB-T And H100Then YO1XGC=J：G+NJYGCFSKD=J：Abs(M2-(X-（G+N）)2) JYKDCQ,+C,-Q=J：J- Abs (Z+0.000) JYSBSJGC=J：G+N+(M2-(Z+0.000)2) JYGCCQ,C+,Q-=J：X-JJElseYO1XGC=J：G+NJYSBSJGC=J：G+N+(M2-(Z+0.000)2) J IfEnd ：Else ZO1XGC=J：G+NJIfX100：Then ZGCFSKD=J：Abs(M2-(X-（G+N）)2) JZKDCQ,+C,-Q=J：J- Abs (Z+0.000) JZSBSJGC=J：G+N+(M2-(Z+0.000)2) JZGCCQ,C+,Q-=J：X-JJElseZSBSJGC=J：G+N+(M2-(Z+0.000)2) J IfEnd ：IfEnd6高程计算主程序子程序：GCSJKIf H起点桩号AndH第一个竖曲线起点桩号 Then第一竖曲线交点高程A：第一竖曲线交点桩号B：第一竖曲线前坡度C：第一竖曲线后坡度D：第一竖曲线半径R：IfEnd：If H第一竖曲线止点桩号AndH第二竖曲线起点桩号 Then第一竖曲线交点高程A：第一竖曲线交点桩号B：第一竖曲线前坡度C：第一竖曲线后坡度D：第一竖曲线半径R：IfEnd：If H第一竖曲线止点桩号AndH第二竖曲线起点桩号 Then第一竖曲线交点高程A：第一竖曲线交点桩号B：第一竖曲线前坡度C：第一竖曲线后坡度D：第一竖曲线半径R：IfEnd(继续添加要素)说明：第一部分坐标部分（1、2、3、4）V=1进入坐标正算V=2进入坐标反算 V=3进入单独的高程计算当V不等于1、2、3时，则返回程序，要求再次输入V值。变量说明：S.起点里程 A.起点X坐标 B.起点Y坐标 C.起点方位角I .起点半径 K.终点里程 T .终点半径 L.计算点里程J输入左右（左-，零，右+）角度 Q.中桩到边桩的距离SDZJ-Z+Y 为隧中支距，左-右+坐标正算中输入的Z值为线路支距坐标反算中得出的Z值也为线路支距注意：1：计算中桩坐标 J D,JL 等于零。2：这里的方位角是指 主点作为起点 的方位角例如：（ZY , YZ , ZH ,QZ, HY , YZ）每个点都是 起点 所以特别注意3、规定(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时， 半径半径*-1；当线元往右偏时，半径=半径；当线元为直线时，半径=以10的45次代替。(2) 当所求点位于中线时，Q,J=0；当位于中线左侧时，J取负值；当位于中线右侧时，J取正值。(3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。 (4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。 (5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。 (6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。 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高程计算主程序：SQXBG 6高程计算主程序子程序：SJK）1、A：交坡点高程2、B：交坡点桩号3、C：前坡度（上坡为正，下坡为负）4、D：坡度（上坡为正，下坡为负）5、R：竖曲线半径6、H：待求点桩号7、E：下一竖曲起点桩号8、XSCGC 测点实测高程9、MRO1点半径（此处如须计算开挖线加上二衬厚度、沉降量、其他结构层的厚度即可）10、NCS线路中桩设计高程至O1点的距离11、JYO1XGCO1点（圆心点）高程12、 JYSBSJGC 根据桩号及支距（正反算程序中）计算得到的测点（设计）高程JYGCFSKD 反算中根据高程反算得到的支距JYKDCQ,+C,-Q 实测高程计算得支距与坐标反算得支距进行对比,得到超欠挖JYGCCQ,C+,Q-实测高程与坐标反算得到高程进行对比,得到超欠挖13、JZO1XGCO1点（圆心点）高程14、J YSBSJGC 根据桩号及支距（正反算程序中）计算得到的测点（设计）高程JZGCFSKD 反算中根据高程反算得到的支距JZKDCQ,+C,-Q 实测高程计算得支距与坐标反算得支距进行对比,得到超欠挖JZGCCQ,C+,Q-实测高程与坐标反算得到高程进行对比,得到超欠挖（本人施测的是一公路隧道的上行线，线路中心线至隧道中心线的距离为0.125米，线路中心线在隧道中心线的左侧。ZSDZJ-Z+Y为隧中支距，Z+0.125为线路支距。如线路中线在隧道中线的右侧，则为Z-0.125为线路支距。其中Z在左侧取-， 右侧取+。）SQXBG 修改说明 此式中X为实测的该点高程，V=1和3时，请将X输为0，当V=2时，输入为实测的高程数据。一般在测量时，将拱部点测得三维后，反算得到里程和支距，即用实测高程计算得到该高程点的理论支距，再用理论支距与坐标反算的支距进行对比，超欠挖自明。同时用反算的支距也可计算出该点的理论高程，用实测的高程与计算的理论高程对比，自可得到拱部是高度是否超欠挖。Lbi1 ：SCGC=X？X： R=M？M：CS=N？N：ProgGCSJK：C-DF：Abs(RF2) T：R AbsFFR：IfHB-T ： Then0K：GOTO2：IfEnd ：If HB-T And H100 Then YGCFSKD=J：Abs(M2-(X-（G+N）)2) JYKDCQ,+C,-Q=J：J- Abs (Z+0.000) JYSBSJGC=J：G+N+(M2-(Z+0.000)2) JYGCCQ,C+,Q-=J：X-JJElseYSBSJGC=J：G+N+(M2-(Z+0.000)2) J： IfEnd ：Else ：ZO1XGC=J：G+NJIfX100：Then ZGCFSKD=J：Abs(M2-(X-（G+N）)2) JZKDCQ,+C,-Q=J：J- Abs (Z+0.000) JZSBSJGC=J：G+N+(M2-(Z+0.000)2) JZGCCQ,C+,Q-=J：X-JJElseZSBSJGC=J：G+N+(M2-(Z+0.000)2) J IfEnd ：IfEnd如在4850中请将XX= YY=改为X=：XY=：Y将LLC= QJL=改为LC=：LJL=：Q将GXLZG=A-（B-H）C-K22R改为XLZG=： G=A-（B-H）C-K22R将JYO1XGC=G+N改为JYO1XGC=：G+N将JYKDCQ,+C,-Q= J- Abs (Z+0.125) 改为JYKDCQ,+C,-Q= ：J- Abs (Z+0.125) 将JYSBSJGC=G+N+(M2-(Z+0.125)2) 改为JYSBSJGC: G+N+(M2-(Z+0.125)2) 将JYGCCQ,C+,Q-=X-J 改为JYGCCQ,C+,Q- : X-J左侧同理。0403 程序说明1、程序中有桩号判断，超出计算范围，不会进行计算，要求你重新输入桩号。2、本程序中有X、Y的最大、最小计算判断，超出X、Y的最大、最小值范围，不会进行计算，要求你重X、Y值。3、本程序反算中有宽度方向和高度方向的超欠挖。4、算例部分见第一楼.提问者：lbs326 人气：1780 赏分：+62 查看楼主其它问题回复内容只查看楼主与本人回复按支持人数排序沙发：lbs326 楼主 10-4-4 0:29 算例04031、平曲线要素交点号JD37交点桩号SK174+133.90 交点坐标X=3845932.100 交点坐标Y=499588.222右偏角 455441前方位角 2971817 半径=1200 LS1=LS2=150ZH点桩号SK173+550.394 HY点桩号SK173+700.394YH点桩号SK174+511.96HZ点桩号SK174+661.960A1=A2=424.264T1=T2=583.596SK173+300QD点坐标X=3845549.528Y=500329.288 C=2971817SK173+550.394ZH点坐标X=3845664.390Y=500106.7933C=2971817SK173+700.394HY点坐标X=3845735.948Y=499974.9917C=3005308.5SK174+106.177 QZ点坐标X=3845998.621Y=499668.235C=3201537.47K174+511.960YH点桩号X=3846348.183Y=499465.9908C=3393806.45ZBQXYSIfL173300 And L173550.394：Then 173300S：3845549.528A：500329.288B：2971817C：1E45I：173550.394K：1E45T： IfEnd：IfL173550.394 And L173700.394：Then 173550.394S：3845664.390A：500106.7933B：2971817C：1E45I：173700.394K：1200T ：IfEnd：If L173700.394 And L174106.177：Then 173700.394S：3845735.948A：499974.9917B：005308.55C：1200I：174106.177K：1200T：IfEnd：If L174106.177 And L174511.960：Then174106.177S：3845998.621A ：499668.235B：3201537.47C：1200I：174511.96K：1200T ： IfEnd：If L174511.96 And L174661.960：Then 174511.960S：3846348.183A：499465.9908B：3393806.45C：1200I：174661.96K：1E45T ：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。GCSJKIfH173.145.441 AndH174187.949 ：Then1332.762A：173350B：0.00761C：-0.02431D：24498.033R：IfEnd ：If H174187.949 AndH174900Then1307.352A：174395.429B：-0.02431C：-0.0169D：56000R：IfEnd 3、实例计算如下其中：X”SCGC”输入任意实测高程值M”R”=5.98 N”CS”=1.602V=1 L输入173806JL输入0.125SDZJ,-Z+Y后显示0 X=3845794.175Y=499886.936”FWJ=“305.928 XLZG=1321.677YO1XG=1323.279J”ZSBSJGC”=1329.259V=1 L输入173806JL输入-4.875SDZJ,-Z+Y后显示-5 X=3845790.126Y=499884.002”FWJ=“305.928 XLZG=1321.677YO1XG=1323.279J”ZSBSJGC”=1326.559V=1 L输入173806JL输入5.125SDZJ,-Z+Y后显示5 X=3845798.224Y=499889.869”FWJ=“305.928 XLZG=1321.677YO1XG=1323.279J”ZSBSJGC”=1329.559J”ZGCFSKD”= Abs(M2-(X-（G+N）)2) J”ZKDCQ,+C,-Q”=J-Abs(Z+0.125) J”ZSBSJGC”=G+N+(M2-(Z+0.125)2) J”ZGCCQ,C+,Q-”=X-J均为相应值V=2 L输入173806X0=3845794.175Y=499886.9355 LC=173805.9996 JL=0.125210SDZJ,-Z+Y=0.000210 XLZG=1321.677YO1XG=1323.279J”ZSBSJGC”=1329.25