CASIO 9860G SD线元法隧道三维（坐标正反计算、高程）计算程序
1. A（此为主程序）
Lbl 0：〝1.LC=>XY〝：〝2.XY=>LC〝：〝3.ZHZL=>GC〝：〝PB=>V=1,2,3〞？→V：
If V=1：Then GOTO 1 ：IfEnd ：If V=2：Then GOTO 2 ：IfEnd ：
If V=3：Then GOTO 3 ：Else GOTO 0 ：IfEnd： Lbl 3：〝ZH=H〝？→H ： 〝SDZF=Z〝？→Z： Prog〝ZGCZCX〝： GOTO 0：Lbl 1 ：〝ZH=L〝？→L： If L>173000 And L<174661.96：Then GOTO 4 ： Else GOTO 0 ： IfEnd ： Lbl 4：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〝：GOTO 0： Lbl 2： 〝XO=M〝?→M：〝YO=R〝?→R：173300→L： If M>3845505.273 And M<3846506.099
And R>499371.832 And R<500352.224 ：Then GOTO 5：Else GOTO 2： IfEnd ：Lbl 5：0→Q：0→J：
Prog 〝ZBFS〝：〝LC=L〝：L ▲〝JL=Q〝：Q▲ 〝SDZJ，-Z+Y=Q〝：Q-0.125→Q▲
L→H：Q→Z： Prog 〝ZGCZCX〝：GOTO 0 2.正算坐标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)+Qcos P→X ： B+AbsH/3*(sinC+4(sinG+sinE)+2sinW+sinO)+Qsin P→Y 3.反算坐标： ZBFS
Lbl 0：Prog 〝ZBQXYS〝 ：Prog 〝ZBZS〝：O-90→Z：(R-Y)cosZ-(M-X) sinZ→P ： L+P→L：If Abs P≥0.001：Then GOTO 0 ：Else GOTO 1： IfEnd ： Lbl 1： (R-Y)cosO-(M-X) sinO→Q 4. 曲线元要素数据库：ZBQXYS
If L≥S And L＜K ：Then **→ S：**→ A：**→ B：**→ C：**→ I：**→ K：**→ T IfEnd ：
If L≥S And L＜K ：Then **→ S：**→ A：**→ B：**→ C：**→ I：**→ K：**→ T IfEnd ：
If L≥S And L＜K ：Then **→ S：**→ A：**→ B：**→ C：**→ I：**→ K：**→ T IfEnd ：
If L≥S And L＜K ：Then **→ S：**→ A：**→ B：**→ C：**→ I：**→ K：**→ T IfEnd ：
If L≥S And L＜K ：Then **→ S：**→ A：**→ B：**→ C：**→ I：**→ K：**→ T IfEnd ：
……………………………
If L≥S And L＜K ：Then **→ S：**→ A：**→ B：**→ C：**→ I：**→ K：**→ T IfEnd ┘
(注：如有多个曲线元要素继续添加入数据库ZBQXYS中) 5 高程计算主程序 ZGCZCX (后有修改说明)
Lbi1 ：〝SCGC=X〝？→X：〝R=M〝：5.98→M：〝CS=N〝：1.603→N：Prog＂GCSJK＂：C-D→F：Abs(RF÷2) →T：R AbsF÷F→R：If H≤B-T ： Then 0→K：GOTO 2： IfEnd ：If H>B-T And H 〝YSBSJGC=J〝：G+N+√(M^2-(Z+0.000)^2) →J▲ IfEnd ： Else 〝ZO1XGC=J〝：G+N→J▲If X>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 ：IfEnd 6高程计算主程序子程序：GCSJK If H>起点桩号 And H≤第一个竖曲线起点桩号 Then 第一竖曲线交点高程→A：第一竖曲线交点桩号→B：第一竖曲线前坡度→C：第一竖曲线后坡度→D：第一竖曲线半径→R：IfEnd： If H>第一竖曲线止点桩号 And H≤第二竖曲线起点桩号 Then 第一竖曲线交点高程→A：第一竖曲线交点桩号→B：第一竖曲线前坡度→C：第一竖曲线后坡度→D：第一竖曲线半径→R：IfEnd： If H>第一竖曲线止点桩号 And H≤第二竖曲线起点桩号 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、X〝SCGC〝 测点实测高程 9、M〝R〝 O1点半径（此处如须计算开挖线加上二衬厚度、沉降量、其他结构层的厚度即可） 10、N〝CS〝 线路中桩设计高程至O1点的距离 11、J〝YO1XGC〝 O1点（圆心点）高程 12、 J〝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.125为线路支距。其中Z在左侧取-， 右侧取+。） SQXBG 修改说明 此式中X为实测的该点高程，V=1和3时，请将X输为0，当V=2时，输入为实测的高程数据。一般在测量时，将拱部点测得三维后，反算得到里程和支距， 即用实测高程计算得到该高程点的理论支距，再用理论支距与坐标反算的支距进行对比，超欠挖自明。同时用反算的支距也可计算出该点的理论高程，用实测的高程与计算的理论高程对比，自可得到拱部是高度是否超欠挖。 Lbi1 ：〝SCGC=X〝？→X：： 〝R=M〝？→M：〝CS=N〝？→N：Prog＂GCSJK＂：C-D→F：Abs(RF÷2) →T：R AbsF÷F→R：If H≤B-T ： Then 0→K：GOTO 2： IfEnd ：If H>B-T And H 〝YSBSJGC=J〝：G+N+√(M^2-(Z+0.000)^2) →J： IfEnd ：▲ Else ：〝ZO1XGC=J〝：G+N→J▲If X>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 ：IfEnd 如在4850中请将X〝X=〝 ▲ Y〝Y=〝▲改为〝X=〝：X◢〝Y=〝：Y◢ 将L〝LC=〝 ▲Q〝JL=〝▲改为〝LC=〝：L◢〝JL=〝：Q◢ 将G〝XLZG〝=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 左侧同理。 0403 程序说明 1、程序中有桩号判断，超出计算范围，不会进行计算，要求你重新输入桩号。 2、本程序中有X、Y的最大、最小计算判断，超出X、Y的最大、最小值范围，不会进行计算，要求你重X、Y值。 3、本程序反算中有宽度方向和高度方向的超欠挖。 4、算例部分见第一楼. 提问者：lbs326 人气：1780 赏分：+62 查看楼主其它问题 → 回复内容 只查看楼主与本人回复 按支持人数排序 沙发：lbs326 楼主 算例0403 1、平曲线要素 交点号JD37 交点桩号SK174+133.90 交点坐标X=3845932.100 交点坐标Y=499588.222 右偏角 45°54′41″ 前方位角 297°18′17″ 半径=1200 LS1=LS2=150 ZH点桩号SK173+550.394 HY点桩号SK173+700.394 YH点桩号SK174+511.96 HZ点桩号SK174+661.960 A1=A2=424.264 T1=T2=583.596 SK173+300 QD点坐标X=3845549.528 Y=500329.288 C=297°18′17″ SK173+550.394 ZH点坐标X=3845664.390 Y=500106.7933 C=297°18′17″ SK173+700.394 HY点坐标X=3845735.948 Y=499974.9917 C=300°53′08.5″ SK174+106.177 QZ点坐标X=3845998.621 Y=499668.235 C=320°15′37.47″ K174+511.960 YH点桩号X=3846348.183 Y=499465.9908 C=339°38′06.45″ ZBQXYS If L≥173300 And L＜173550.394：Then 173300→S：3845549.528→A：500329.288→B：297°18′17″→C：1E45→I：173550.394→K：1E45→T： IfEnd： If L≥173550.394 And L＜173700.394：Then 173550.394→S：3845664.390→A：500106.7933→B：297°18′17″→C：1E45→I：173700.394→K：1200→T ：IfEnd： If L≥173700.394 And L＜174106.177：Then 173700.394→S：3845735.948→A：499974.9917→B：00°53′08.55″→C：1200→I：174106.177→K：1200→T： IfEnd： If L≥174106.177 And L＜174511.960：Then 174106.177→S：3845998.621→A ：499668.235→B：320°15′37.47″→C：1200→I：174511.96→K：1200→T ： IfEnd： If L≥174511.96 And L＜174661.960：Then 174511.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。 GCSJK If H>173.145.441 And H≤174187.949 ：Then 1332.762→A：173350→B：0.00761→C：-0.02431→D：24498.033→R：IfEnd ： If H>174187.949 And H≤174900 Then 1307.352→A：174395.429→B：-0.02431→C：-0.0169→D：56000→R：IfEnd 3、实例计算如下 其中： X”SCGC”输入任意实测高程值 M”R”=5.98 N”CS”=1.602 V=1 L输入173806 JL输入0.125 SDZJ,-Z+Y后显示0 X=3845794.175 Y=499886.936 ”FWJ=“305.928 XLZG=1321.677 YO1XG=1323.279 J”ZSBSJGC”=1329.259 V=1 L输入173806 JL输入-4.875 SDZJ,-Z+Y后显示-5 X=3845790.126 Y=499884.002 ”FWJ=“305.928 XLZG=1321.677 YO1XG=1323.279 J”ZSBSJGC”=1326.559 V=1 L输入173806 JL输入5.125 SDZJ,-Z+Y后显示5 X=3845798.224 Y=499889.869 ”FWJ=“305.928 XLZG=1321.677 YO1XG=1323.279 J”ZSBSJGC”=1329.559 J”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”ZGCCQ,C+,Q-”=X-J 均为相应值 V=2 L输入173806 X0=3845794.175 Y=499886.9355 LC=173805.9996 JL=0.125210 SDZJ,-Z+Y=0.000210 XLZG=1321.677 YO1XG=1323.279 J”ZSBSJGC”=1329.259 J”ZGCFSKD” J”ZKDCQ,+C,-Q” J”ZGCCQ,C+,Q-” 均为相应值 V=2 L输入173806 X0=3845790.126 Y=499884.002 LC=173805.9996 JL=-4.875065 SDZJ,-Z+Y=-5.000065 XLZG=1321.677 YO1XG=1323.279 J”ZSBSJGC”=1329.559 J”ZGCFSKD” J”ZKDCQ,+C,-Q” J”ZGCCQ,C+,Q-” 均为相应值 V=2 L输入173806 X0=3845798.224 Y=499889.869 LC=173806.0004 JL=-5.124899 SDZJ,-Z+Y=4.999899 XLZG=1321.677 YO1XG=1323.279 J”ZSBSJGC”=1329.559 J”ZGCFSKD” J”ZKDCQ,+C,-Q” J”ZGCCQ,C+,Q-” 均为相应值 V=3 ZH输入桩号173806 SDZF,-Z+Y =0 XLZJ输入0.125 XLZG=1321.677 YO1XG=1323.279 J”ZSBSJGC”=1329.259 V=3 ZH输入桩号173806 SDZF,-Z+Y =-5 XLZJ输入-4.875 XLZG=1321.677 YO1XG=1323.279 J”ZSBSJGC”=1329.559 V=3 ZH输入173806 SDZF,-Z+Y =5 XLZJ输入5.125 XLZG=1321.677 YO1XG=1323.279 J”ZSBSJGC”=1329.559 因篇幅问题不能全部显示，请点此查看更多更全内容