最新9860G线元法隧道坐标正反计算、高程计算程序汇总
9860G线元法隧道坐标正反计算、高程计
算程序
CASIO 9860G SD线元法隧道三维（坐标正反计算、高程）计算程序悬赏分：+62
作者：刘工
2010-4-3
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＂G CSJK＂：C-D→F：Abs(RF÷2) →T：R AbsF÷F→R：If H≤B-T ： Then 0→K：GOT O 2：
IfEnd ：If H>B-T And H If H≥B+T ： Then 0→K：D→C：GOTO 2：IfEnd ： Lbi 2 ：〝XLZG=G〝：A-（B-H）C-K^2÷2R→G▲If Z≥0： Then If X>100 Then 〝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▲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高程计算主程序子程序：S JK） 1、A：交坡点高程 2、B：交坡点桩号 3、C：前坡度（上坡为正，下坡为负） 4、D：坡度（上坡为正，下坡为负） CASIO 9860G SD线元法隧道三维（坐标正反计算、高程）计算程序 悬赏分：+62 作者：刘工 2010-4-3 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+16 U→Ｏ： 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 ： 公用模块： Option Explicit Public Const PI = 3.14159265358979 '已知A、B两点坐标计算方位角，JSFWJ的中文意思是计算方位角 Public Function JSFWJ(xa As Double, ya As Double, xb As Double, yb As Double) As Double '已知A、B两点坐标计算方位角函数过程Dim vx As Double, vy As Double vx = xb - xa: vy = yb - ya '如果A、B两点坐标相同，出现提示对话框 If vx = 0 And vy = 0 Then MsgBox "您选择的是同一个点！", vbOKOnly + vbExclamation, "提示信息" JSFWJ = 999999999# End If '计算方位角的值 If vx = 0 And vy > 0 Then '与y轴正半轴平行 JSFWJ = RadianToAngle(PI / 2#) ElseIf vx = 0 And vy < 0 Then '与y轴负半轴平行 JSFWJ = RadianToAngle(PI * 3# / 2#) ElseIf vy = 0 And vx > 0 Then '与x轴正半轴平行 JSFWJ = RadianToAngle(0) ElseIf vy = 0 And vx < 0 Then '与x轴负半轴平行 JSFWJ = RadianToAngle(PI) ElseIf vx > 0 And vy > 0 Then '第一象限 JSFWJ = RadianToAngle(Atn(vy / vx)) ElseIf vx < 0 And vy > 0 Then '第二象限 JSFWJ = RadianToAngle(Atn(vy / vx) + PI) ElseIf vx < 0 And vy < 0 Then '第三象限 JSFWJ = RadianToAngle(Atn(vy / vx) + PI) ElseIf vx > 0 And vy < 0 Then '第四象限 JSFWJ = RadianToAngle(Atn(vy / vx) + 2 * PI) End If End Function '已知A、B两点坐标计算距离，JSJLS的中文意思是计算距离S Public Function JSJLS(xa As Double, ya As Double, xb As Double, yb As Double) As Double Dim vx As Double, vy As Double vx = xb - xa: vy = yb - ya '如果A、B两点坐标相同，出现提示对话框 If vx = 0 And vy = 0 Then MsgBox "您选择的是同一个点！", vbOKOnly + vbExclamation, "提示信息" JSJLS = 99999999# End If '计算距离 JSJLS = Sqr(vx * vx + vy * vy) End Function '弧度化角度 Public Function RadianToAngle(ByVal alfa As Double) As Double Dim alfa1 As Double, alfa2 As Double alfa = alfa * 180# / PI 曲线任意里程中边桩坐标正反算(CASIO fx-5800P计算器)程序 一、程序功能及原理 1.功能说明：本程序由一个主程序(TYQXJS)和五个子程——正算子程序(SUB-ZS)、反算子程序(SUB-FS)等构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。本修改版程序既可实现正算全线贯通，亦可实现反算全线贯通。本程序在CASIO fx-5800P计算器运行。 2．计算原理：利用Gauss-Legendre 5点通用公式正算线路中边桩坐标、线外测点至曲线元起点和终点的垂距的符号是否相异（即Dca×Dcb<0=>该测点在其线元内）进行判断并利用该线元要素反算中桩里程、支距，最后计算出放样数据。 二、源程序 1.主程序(TYQXJS)(A) Deg:fix 3 119→DimZ “INPUT(0) Or DATA(Else)”?I Lbl 0:“1.SZ=>XY,2.XY=>SZ，3.TF=>CK,4.SD=>FY，5.TW=>FY”?N If N=1 Or N=5:Then Goto 1 Else If N=2 Or N=3 Or N=4:Then Goto 2 Else Goto 3 IfEnd:IfEnd Lbl 1:“K(m)=”?S If S<0:Then Goto 0:IfEnd “JL(m)=”?Z If Z≠0:Then “ANGLE→R(Deg)=”?M:IfEnd If I=0:Then Prog “DAT1”:Else Prog “DAT2”:IfEnd S-O→W:If W<0:Then Goto 0:Else If W>H:Then Goto 0:IfEnd:IfEnd Prog “SUB-ZS”:Prog “SUB-GC” If Z<0:Then“XL(m)=”:X◢“YL(m)=”:Y◢ If N=5:Then Prog “SUB-TW”:IfEnd Else If Z>0:Then “XR(m)=”:X◢“YR(m)=”:Y◢ If N=5:Then Prog “SUB-TW”:IfEnd Else “X(m)=”:X◢“Y(m)=”:Y◢“Hs(m)=”:L◢“FWJ=”: F?DMS◢ IfEnd:IfEnd CASIO fx5800p全线高程计算程序 GAOCHEN 主程序 Lbl 1 “KM=，<0,Stop”:?K：K<0=>Stop:“PY=”？L：Prog”GK” C-D→E:Abs(RE/2)→T：R(Abs(E)/E)→R If K≤B-T:Then 0→H:Else:If K≥B+T Then 0→H:D→C:Else K-B+T→H:Ifend:Ifend A-(B-K)C-H2/(2R)→G：Cls “KM=”：Locate 4，1，K：Locate 10，1，“PY=”：Locate 13，1，L：Fix 3 “H=”：Locate 4，2，G Prog “PODU”：（E-B）/（D-A）（K-A）+B→I：(F-C)/（D-A）（K-A）+C→J “HL=”：G+IL→X：Locate 4，3，X：Locate 11，3，“I=”：Locate 13，3，I*100 “HR=”：G+JL→Y：Locate 4，4，Y：Locate 11，4，“I=”：Locate 13，4，J*100◢Cls：Norm 2：“BM+HS≤0,Goto 1”？Z：Z≤0=> Goto 1：Cls (输入视线高) “KM=”：Locate 4，1，K：Locate 10，1，“PY=”：Locate 13，1，L：Fix 3 “QSM=”: Locate 6，2，Z-G (显示中桩读数) “QSL=”: Locate 6，3，Z-X (显示左桩读数) “QSR=”: Locate 6，4，Z-Y◢(显示右桩读数) Norm 2：Cls：Goto1 （后面可加已知视线高计算读数部分，不想计算读数则视线高输入0或负数如不想显示麻烦，可将Locate语句去掉） 以下两个子程序不需运行，只是两个独立的数据库赋值程序，字母重复不影响计算结果 GK 数据库子程序 If K≤第二曲线起点桩号：Then 第一曲线交点高程→A：第一曲线交点桩号→B：第一曲线前坡→C：第一曲线后坡→D：第一曲线半径→R：Return：Ifend …………….（有几个变坡点编几个If语句） PODU 计算坡度子程序 If K≤第一变（非变）坡段终点：Then 第一曲线起点桩号→A：第一曲线起点左坡→B：第一曲线起点右坡→C：第一曲终点桩号→D：第一曲终点左坡→E：第一曲终点右坡→F：Return：Ifend ………………（每一个超高变化线元一个If语句） 结果显示： KM=0000.000 PY=0.000 H= 00.000 HL=00.000 I=-1.5 HR=00.000 I=-1.5 KM=0000.000 PY=0.000 QSM= 00.000 QSL= 00.000 QSR= 00.000 在测量岗位工作已经有三个月到时间了，三个月的时间学习和收获了许多，现对这三个月的工作学习做一下总结。 测量工作内容主要有以下两个方面：测量放线（坐标计算），高程控制。 一、测量放线 测量放线到主要技术包括坐标计算和仪器使用。坐标计算包括直线段坐标计算和曲线段坐标计算。 1、直线段坐标计算。直线坐标计算分为中桩坐标计算和边桩坐标计算。 1）中桩坐标计算。根据公式 ααsin ,cos d Y Y d X X +=+=起中起中 d — 所求点到起点距离； α— 该直线坐标方位角。在此顺带详细介绍一下坐标方位角到计算方法： （1）坐标方位角的计算 AB AB A B A B AB x y x x y y ??=--=arctan arctan α当 R y x R y x R y x R y x -360,0,0180,0,0-180,0,0;,0,0?=?+?=??>?αααα；； （2）坐标方位角的推算 ， ， 218021*********βαβααβαβαα-?+=-=+?+=+=B B AB BA B 由此推出：βαα±?+=180后前（“左”→“+”， “右”→“-”），计算中，若α值大于360°，应减去360°；若小于0°，则加上360°。 2）边桩坐标计算 应用公式 )90sin(90cos(?±+=?±+=ααl y y l x x 中边中边）， 进行边桩坐标到计算。北客站为直线车站，坐标计算较简单，现以位于机场线第二段底板的变电所夹层东北角C 点为例进行计算： 以机场线右线为基准来计算中、边桩坐标。已知起点坐标A （22264.4009，11553.2031），终点坐标B （22180.2655，11279.0739），起点里程为YDK0+255.275,C 点里程为YDK0+286.075,偏距为15.33m ，则由以上公式计算C 点坐标： α=arctan(（11279.0739-11553.2031）/（22180.2655-22264.4009）)+180°=252.938°， =中x 22264.4009+（286.075-255.275）*cos252.938°=22255.3640 =中y 11553.2031+（286.075-255.275）*sin252.938°=11523.7586 =c x +15.33*cos （252.938°+90°）=22270.0193 = c y +15.33*sin （252.938°+90°）=11519.2606，则可求出C （22270.0193，11519.2606）。 2、曲线段坐标计算 1）不带缓和曲线的圆曲线中、边桩坐标计算 北 中 x 中 y 公路高程计算公式 ⒈超高 ①超高方式：中央分隔带保持水平， 超高将两侧行车道绕中央分隔带边缘点旋转(包括路肩点)。 ②超高段横断面高程计算 图11 1-1 A 0％ A A A A A A A D D D 3-3 4-4 I％I％ I％ E％ I％ 图12 ⒉ 横坡度计算 外侧横坡度： I L L E I I C X X ?+= )(；(公式中的I 、E 均取正值，下同。) 公式 1 内侧横坡度： I L I E I L L I E I L I E I C C C X X +?+??+? ?= 22)((。 公式 2 1 23 L C I B 左0％ E E +I Z I E +I l c l c I E +I E +I l c I l c E +I I 右 （左 ) 式中：2 I/(E+I)* L C—在L C段内横坡等于I％的长度,m。 X在区间0～2 I/(E+I)*LC时，横坡度为I； 在区间2 I/(E+I)* L C～L C段内时，横坡度为I～E。 I—横坡度设计值， E—超高设计值， L C—缓和曲线长，m。 ⒊竖曲线计算公式： W=I1-I2；当w＞0时，为凸曲线；当w＜0时，为凹曲线。 L=R*W； E=T2/2R； H=l2/2r； T=TA=TB=L/2=R*W/2。 式中：H—切线上任一点至竖曲线上的垂直距离；M． l—曲线上相应于H的P点至切点A或B点的距离，M． R—二次抛物线的参数。(原点处的曲率半径)通常称竖曲率半径，M． I1、I2—切线的斜率，即纵坡度，％．纵坡度(％)，从左向右上坡取“+”，下坡取“-”值．当α很小时， 一个excle 模板的制作 在当今社会，excle的使用已经是越来越来频繁了，几乎涉及所有的行业，路桥施工也 不例外。我在某路桥公司曾经负责过某项目部的测量工作。大家都知道，测量最主要的就 是计算了，如坐标、高程、横坡度等。我现在给大家推荐一款我自己编制的关于测量计算 的excel模板。 首先我会跟大家介绍一下模板的作用，然后再一一讲解此模板的制作过程。 首先给大家看一下此模板的界面如下： 也许大家咋一看，切~ 这算啥，我也会做这张表格，实在是太简单了。不错，如果仅 仅是靠手动输入这样子的数字，也许只要懂一点点excle 的人都会制作出这张表格吧。不过，这张表格并不是你表面所看到的仅仅是几个数字而已，其内在的公式才是它的亮点。也许 这样讲大家还不是很清楚，我继续给大家截个图，看看它里面的公式是什么。 大家注意到上面的公式了吗，并不是仅仅是输入数字就完事的，它是一个自定义函数 zbx（），那么后面的都是一样吗？完全正确，后面的都是自定义函数，它们分别是zby（）、sqx（）、hpz（）、hpy（）。也许大家会问，恩，是不错，但是有什么用呢？那让我先给大家 简述一下这个自定义函数的用法。竟然是一个函数，那么它就必须要有一个自变量，这几个 函数的自变量又是什么呢？其实这个模板里面所有函数的自变量只有一个，就是桩号。什么 意思？就是只要你给出任意一个桩号，都能得到其对应的坐标、中桩高程和横坡度。假设我 们要K38+000~K38+200 段落内每隔20M 一个断面所有点的坐标、中桩高程、以及左右横坡。我就用这个模板给大家演示一下（此模板暂时数据只针对黄祁高速公路六标项目部）。 先在桩号那一列把K38+000~K38+200 输入进去，可不要真的把字母“K”和加号 CASIO5800计算器 公路测量计算程序 程序设计:魏加训 2009.2.28 Casio 5800计算器数据库型万能坐标正反算计算程序 一、主程序：1XY Lbl 0:“1.ZS 2.FS”?→V↙ If V=1: Then “CZ X” ?H: “CZ Y” ?T:Goto 1: Else If V=2: Then Goto 2 :IfEnd: IfEnd ↙ Lbl 1: “K××+×××”?D:?Z: “RJ”?G↙ Prog “P.Z”↙(注:计算另一线路时修改替换此处和FS子程序中的P.Z为对应线路的数据库名称即可) Prog “ZS” ↙ If Z<0: Then Cls:“X（L）=”: “Y（L）=”: Locate 6，1，X : Locate 6，2，Y◢ Pol（X-H，Y-T）: Cls: “S（L）=”: Locate 6，1，I : "F（L）="：360Frac（（J+360）÷360▼DMS◢ Goto 1:IfEnd↙ If Z=0: Then Cls:“X（Z）=”: “Y（Z）=”: Locate 6，1，X : Locate 6，2，Y : “QXFWJ （Z）=”: 360Frac（（J+360）÷360▼DMS◢ Pol（X-H，Y-T）: Cls: “S（Z）=”: Locate 6，1，I : "F（Z）="：360Frac（（J+360）÷360▼DMS◢ Goto 1:IfEnd↙ If Z >0: Then Cls:“X（R）=”: “Y（R）=”: Locate 6，1，X : Locate 6,2，Y◢ Pol（X-H，Y-T）: Cls: “S（R）=”: Locate 6，1，I : "F（R）="：360Frac（（J+360）÷360▼DMS◢ Goto 1:IfEnd↙ Lbl 2: 0→Z:0→G:”X” ?M:”Y” ?I:Prog “FS”: Cls:“K=”：“Z=”：Locate 4，1，D : Locate 4，2，Z◢ Goto 2↙ 二、正算子程序：ZS 5→N: U(E-1-R-1)÷Abs（K-F）→P: Abs（D-F）÷N→Q: 90Q÷π→S: (注:此处5→N是控制计算精度可修改的,一般取值为4～6即可) C+(NPQ+2UR-1)NS→J:1→L↙ A+Q÷6×(Cos (C)+Cos (J) +4∑（Cos (C+((L+0.5)PQ+2UR-1)×(L+0.5)S),L,0,(N-1)）+2∑（Cos (C+((LPQ+2UR-1)LS,L,1,(N-1))）+ZCos（J+G) →X : B+Q÷6×(Sin(C)+Sin( J) +4∑（Sin (C+((L+0.5)PQ+2UR-1)×(L+0.5)S),L,0,(N-1)）+2∑（Sin (C+((LPQ+2UR-1)LS,L,1,(N-1))）+Z Sin（J+G)→Y ↙ 三、反算子程序：FS Lbl 0:Prog “P.Z”:Prog “ZS”↙ (注:计算另一线路时修改替换此处和1XY主程序中的P.Z为对应线路的数据库名称即可) (I-Y)sin(J)+(M-X) cos(J)→ P :D+P→ D ↙ If Abs(P)≥0.001:Then Goto 0 : Else Goto 1 : IfEnd↙ Lbl 1: (I-Y)cos(J)-(M-X) sin(J) →Z ４８００公路中边桩高程计算程序 fx-4800高等级公路中边桩标高程计算程序 （适用超高方式为有分隔带绕左内侧旋转，无分隔绕线路中心旋转） 随着高等级公路建设的不断发展，对路基施工高程控制要求不断提高，经常在现场要对路基中边桩高程进行测量，由于施工线路较长，高程数据较多，而且只提供设计桩号的高程，对非设计断面的横坡及高程计算相当复杂，因此利用fx-4800计算器的编程功能及扩展变量的储存功能，编制一个程序，可以预先输入整条线路的曲线要素，使用时，只要输入任意点桩号就可以计算中边桩高程。 程序由中桩高程计算主程序和横坡计算子程序组成，程序如下： 1、中桩高程计算程序文件名（ZZGC） 程序说明 Lb1 0 Fix 3 （输入变坡点数） N=? （输入扩展变量位置数） V=? M=V （输入待求桩号） {K} Lb1 1 K>Z[2+M]-Z[3N+3+M]=> （确定待求桩号在哪个竖曲线内） M=M+1：Goto 1 （读取该竖曲线半径） R=Z[2N+2+M] （读取该竖曲线切线长） T=Z[3N+2+M] （读取该变坡点前坡率） I=Z[4N+2+M]÷100 （读取该变坡点后坡率） L=Z[4N+3+M]÷100 （计算待求桩号与该竖曲线起点距离） S=K-Z[1+M]+Z[3N+2+M] M=V=> （如果桩号在起点与第1竖曲线起点之间，转Lb1 5 ） Goto 5⊿ （计算坡率差） A=I-L （判断竖曲线凹凸） J=AbsA÷A S 一个建筑设计师应知道的基本数据 一、普通住宅建筑混凝土用量和用钢量： 1、多层砌体住宅： 钢筋：30KG/m2 砼：0.3～0.33m3/m2 2、多层框架： 钢筋：38～42KG/m2 砼：0.33～0.35m3/m2 3、小高层11～12层： 钢筋：50～52KG/m2 砼：0.35m3/m2 4、高层17～18层： 钢筋：54～60KG/m2 砼：0.36m3/m2 5、高层30层H=94米：钢筋：65～75KG/m2 砼：0.42～0.47m3/m2 6、高层酒店式公寓28层H=90米： 钢筋：65～70KG/m2 砼：0.38～0.42m3/m2 7、别墅：混凝土用量和用钢量介于多层砌体住宅和高层11～12层之间； 以上数据按抗震7度区规则结构设计 二、普通多层住宅楼施工预算经济指标 1、室外门窗（不包括单元门、防盗门）面积占建筑面积0.20～0.24 2、模版面积占建筑面积2.2左右 3、室外抹灰面积占建筑面积0.4左右 4、室内抹灰面积占建筑面积3.8 三、施工功效 1、一个抹灰工一天抹灰在35平米 2、一个砖工一天砌红砖1000～1800块 3、一个砖工一天砌空心砖800～1000块 4、瓷砖15平米 5、刮大白第一遍300平米/天，第二遍180平米/天，第三遍压光90平米/天 四、基础数据 1、混凝土重量2500KG/m3 2、钢筋每延米重量0.00617×d×d 3、干砂子重量1500KG/m3，湿砂重量1700KG/m3 4、石子重量2200KG/m3 5、一立方米红砖525块左右（分墙厚） 6、一立方米空心砖175块左右 7、筛一方干净砂需1.3方普通砂 建筑程序歌 要想建设效果好，选择队伍要招标。 小白自学ExcelVBA道路坐标高程计算程序流程 1声明 1.1 本人是个工作时间不长测量施工员。 1.2听说想轻松愉快的干好测量得学门编程，又听说ExcelVBA是个简单的可视化的常用的…… 1.3 本人很懒，面对厚厚的全是“昏天暗地”代码的ExcelVBA教科书籍，实在是没心没力。其实俺还是懂点儿VB编程的（学校学过“=”是赋值，if、while、end的意思，嘿嘿，还有画窗体图框，其他的就都还给老师了）。 1.4 谷歌、度娘，是咱的好帮手。 1.5 仰望的存在，神一般的存在，无私的化身！ 道路中边桩坐标计算程序120424.xls（王中伟”教授”QQ：595077） 曲线坐标计算程序VBA 4.6.xls（陈超”中铁”QQ：295188316） 陈师傅的VBA代码乍一看感觉比王老师的让人头痛，怎么说咱们也是个有理想有目标的人，要啃就啃硬骨头！就选~~~嘿嘿~~~王老师的程序吧！ 1.6 道路中边桩坐标计算程序120424.xls代码获取。 王老师、陈师傅担心不懂编程的人，不小心修改代码导致程序错误，设置了vba密码（个人认为），其次我跟他们不熟，没有正大光明的要的密码，于是求助谷歌度娘，这个关于VBA 密码破解，于是很不厚道的得到了密码，为了方便本次学习，又很不厚道的把道路中边桩坐标计算程序120424.xls的密码贴出来（df750726），神人勿怪呀！ 2 代码理解ing… 我们通过密码打开vba代码窗口，大家可以先把每个表格打开浏览一下代码，是不是感觉很短呀，感觉应该是超级链接的意思，嘿嘿，主要的代码在模块“lx”里，打开下，代码多吧，这就是咱们要肯的骨头。现在咱们冲呀！ CASIO fx-5800P单交点通用型曲线坐标高程计算程序 一、说明： 本程序采用交点法计算道路基本型曲线坐标及高程，在建立好数据库后，能连续计算全线各桩号的中边桩坐标及高程。本程序共包括一个主程序和九个子程序，其中有坐标计算、高程计算、坡口坡脚线放样，锥坡放样坐标计算等子程序。 二、内容： 1．PM5-3 XYZJS（主程序） “SINGLE BASIC TYPE CURVE”◢ “METHOD OF COORDINATE PM5-3”◢ Deg:ClrStat:FreqOn:Fix 3 40→DimZ “INPUT(0) Or DATA(Else)”?N “FUNCTION”?P Prog “SUB5-35” If Z[30]<0:Then -1→Z[20]:Else 1→Z[20]:IfEnd Abs(Z[30])→D Pol(Z[26]-Z[28],Z[27]-Z[29]):Cls If J<0: Then J+360→Z[11]:Else J→Z[11]:IfEnd 计算ZH→JD方位角 Z[11]+Z[30]+180→Z[16] 计算HZ→JD方位角 If Z[16]>360:Then Z[16]-360→Z[12]:Else Z[16]→Z[12]:IfEnd If Z[12]>180:Then Z[12]-180→Z[23]:Else Z[12]+180→Z[23]:IfEnd计算JD→HZ方位角 S2÷（24R）-S4÷(2688R3)→Z[1] 计算第一缓和曲线内移值 0.5S-S3÷(240R2)+S5÷(34560R4)→Z[2] 计算第一缓和曲线切线增长值 T2÷(24R)-T4÷(2688R3)→Z[3] 计算第二缓和曲线内移值 0.5T-T3÷(240R2)+T5÷(34560R4)→Z[4] 计算第二缓和曲线切线增长值 (R+Z[3])÷sin(D)-(R+Z[1])÷tan(D)+Z[2]→Z[5] 计算第一切线长 (R+Z[1])÷sin(D)-(R+Z[3])÷tan(D)+Z[4]→Z[6] 计算第二切线长 90S÷(πR)→Z[7]:90T÷(πR)→Z[8] 计算第一、二缓和曲线偏角πR(D-Z[7]-Z[8])÷180→Z[9] 计算圆曲线长度 S+T+Z[9]→Z[10] 计算曲线总长度 Z[25]-Z[5]→List X[1] 计算直缓点桩号 1→K:Prog “SUB5-37” List X[1]+S→List X[2]:List X[2]+Z[9]→List X[3] 计算缓圆点、圆缓点桩号 Z[26]-Z[5]cos(Z[11])→List Y[1]:Z[27]-Z[5]sin(Z[11])→List Freq[1] If S≠0:Then“ZH PEG(m)=”:List X[1]◢显示直缓点桩号 Else “ZY PEG(m)=”:List X[1]◢显示直圆点桩号 IfEnd “X(m)=”:List Y[1]◢显示直缓(圆)点X坐标 “Y(m)=”:List Freq[1]◢显示直缓(圆)点Y坐标 超高及高程（源程序） L1: R:L“LS”:Z“ZH”:H“HZ”:V“RS”:F“I1”:G“I2”:O“BPH”:N“BP”:T=Vabs(G－F)÷2:X“X(R+1,L-1)”:P “P(R+,L-1)” L2: Lbl 6 : ｛S｝ L3: S+T－N<0 D=O+(S－N)F: S－N<0 D=O+(S－N)F－U(S－N+T)2÷2÷V: S－N≤T D=O+(S－N)G－U(T －S+N)2÷2÷V: D=O+(S－N)G D“H=”◢ L4: R< E= : M= Goto 1: R< E= : M= Goto 1: L5: Lbl 1 : M ≥L C=L: Goto 2: C=M L6: Lbl 2 : P X >0 Goto 3 : Goto 4 L7: Lbl 4 : S ≤Z I=0.02 : Goto 5: S 王天悼卡西欧5800P高程高差现场计算小程序 程序功能目的：输入水准仪现场前视读数及设计快速计算出该点的实测高程及与设计值的高差！ 程序名: 0－GCGC－JS LbI A:Cls:“HS=”?C“BM H（M）=”?G:C+G→S LbI B:Cls: “QS=”?D:If D≥0 :Then Goto D:ElseGoto C:IfEn d LbI C:M→G:Abs（D）→C:Goto A LbI D:S-D→M:Cls:“HS=”:“SXG=”: “QS =”:“H=”:Locate 6,1,C: Locate 6,2,S: Locate 6,3,D: Locate 6,4,M◢ Cls: “1 =>JiSuan,2=>QXGT”?U: U=1=>Goto E:U=2=>Goto F LbI E CLS:“H(SJ)=”？J:“QS=”: “H=”: “H(SJ)=”: “h=” Locate 6,1,D: Locate 6,2,M:Locate 6,3,J: Locate 6,4,J-M－0◢ （0可根据实际情况输入，如：垫层厚度为20cm，要求计算路基标高，则输入0.2） Goto B LbI F CLS:Prog“0-QXGT-GC”:H→J:“H=”: “h=”: “h(Z)=”: “h(Y)=”Locate 6,1,M: Locate 6,2,J-M－0:Locate 6,3,V-M－0: Locate 6,4,W-M －0◢（0可根据实际情况输入，如：垫层厚度为20cm，要求计算路基标高，则输入0.2） Goto B 符号说明： HS? 输入水准点上之后视读数 BG:（标高计算式程序） Lbl 1：Abs(.01(I－J)R÷2)－﹥T： If Z≥L+T：Then Goto 4:IfEnd If Z≥L-T：Then Goto 3:IfEnd Lbl 2：H+.01(Z－L)I－﹥B: Goto 5 Lbl 3:H+.01(Z-L)I+(J-I)(Z+T-L)2÷(2R Abs(J-I)) －﹥B:Goto 5 Lbl 4:H+.01(Z-L)J－﹥B:Goto 5 Lbl 5:B H——竖曲线顶点高程 I——后坡 J——前坡 L——竖曲线顶点桩号 B——设计高程 R——圆曲线半径 GX:（挂线高计算主程序） Lbl 0： ―U=‖?U: ―Z=‖?Z: If U=0: Then ?L:?H:?I:?J:?R: ―Z=‖?Z:Prog―BG‖：―B=‖：B▲―P=‖?P:P:Goto 4: IfEnd If U=1: Then Prog―L1‖：Prog―PY‖：Goto 1: IfEnd If U=2: Then Prog―L2‖：Prog―PY‖：Goto 1: IfEnd If U=3: Then Prog―L1‖：Prog―PZ‖：Goto 1: IfEnd If U=4: Then Prog―L2‖：Prog―PZ‖：Goto 1: IfEnd If U=5: Then Prog―L3‖：Prog―PZ‖：Goto 1: IfEnd If U=6: Then Prog―LA‖：Goto 2: IfEnd If U=7: Then Prog―LB‖：Goto 2: IfEnd If U=8: Then Prog―LC‖：Goto 2: IfEnd If U=9: Then Prog―LD‖：Goto 2: IfEnd 高程计算程序 主程序“1” U“SXG”:Lbl 0:Fixm:｛KBX｝回车 K≤***＝＞K≥***＝＞M=***：L=***回车 终点桩号起点桩号计算公式或者固定数值 。。。逐段输入所有数据。。。 K≤***＝＞K≥***＝＞M=***：L=***回车 Prog“2”: Prog“4”:Goto 0回车 说明：M---左幅超高值计算公式 L----右幅超高值计算公式 按超高方式图输入数据 子程序“2” K≤该段竖曲线终点＝＞K≥上段竖曲线终点＝＞A=***：J=***：C=***：D=***：回车 。。。逐段输入所有数据。。。 K≤该段竖曲线终点＝＞K≥上段竖曲线终点＝＞A=***：J=***：C=***：D=***：回车 Prog“3”回车 说明：A---变坡点高程 J---变坡点桩号 C---第一纵坡 D---第二纵坡 V---竖曲线半径 按纵段面图上的竖曲线数据 子程序“3” E=C-D：T=V×Abs(E÷2):Q=J-T:S=J+T回车 E>0＝＞Y=-(K-Q)2÷(2V)：≠＝＞Y=(K-Q)2÷(2V)回车 K≤Q＝＞H=A-(J-K)C: ≠＝＞K≤S＝＞ K≥S＝＞回车 “H=”：H 说明:本程序计算中桩高程 子程序“4” B≥0＝＞W=H+LB：≠＝＞回车 “W=”：W O=W-（U-X）-0.74回车 “T=”：O×1.25 说明:本程序计算任意宽度的高程 B—计算左半幅输负数,右半幅输正数 W—任意点的路面设计高程 0.74—路面的结构层厚度,需逐层调整 T—仍需填的高度,1.25为松铺系数,根据实验室数据及现场实际情况进行调整 坐标、高程计算公式 一个建筑设计师应知道的基本数据 一、普通住宅建筑混凝土用量和用钢量： 1、多层砌体住宅： 钢筋：30KG/m2 砼：0.3～0.33m3/m2 2、多层框架： 钢筋：38～42KG/m2 砼：0.33～0.35m3/m2 3、小高层11～12层： 钢筋：50～52KG/m2 砼：0.35m3/m2 4、高层17～18层： 钢筋：54～60KG/m2 砼：0.36m3/m2 5、高层30层H=94米：钢筋：65～75KG/m2 砼：0.42～0.47m3/m2 6、高层酒店式公寓28层H=90米： 钢筋：65～70KG/m2 砼：0.38～0.42m3/m2 7、别墅：混凝土用量和用钢量介于多层砌体住宅和高层11～12层之间； 以上数据按抗震7度区规则结构设计 二、普通多层住宅楼施工预算经济指标 1、室外门窗（不包括单元门、防盗门）面积占建筑面积0.20～0.24 2、模版面积占建筑面积2.2左右 3、室外抹灰面积占建筑面积0.4左右 4、室内抹灰面积占建筑面积3.8 三、施工功效 1、一个抹灰工一天抹灰在35平米 2、一个砖工一天砌红砖1000～1800块 3、一个砖工一天砌空心砖800～1000块 4、瓷砖15平米 5、刮大白第一遍300平米/天，第二遍180平米/天，第三遍压光90平米/天 四、基础数据 1、混凝土重量2500KG/m3 2、钢筋每延米重量0.00617×d×d 3、干砂子重量1500KG/m3，湿砂重量1700KG/m3 4、石子重量2200KG/m3 5、一立方米红砖525块左右（分墙厚） 6、一立方米空心砖175块左右 7、筛一方干净砂需1.3方普通砂 FX 5800p市政道路高程计算及放样程序 一、主程序 Cls:Fix 3 “CHA-XUN=0,FANG-YANG≠0”？→J “SHI-XIAN-GAO=”?Z “线路1=1，线路2=2，线路3=3” “XIAN LU”?Q Lb1 1 Q=1=>Prog“线路1” Q=2=>Prog“线路2” Q=3=>Prog“线路3” If K<0：Then“END”：Stop：Else Goto 2 Lb1 2 “PIAN-JU”?N:”JIE-GOU-HOU-DU”?E S-K→L:A-B→W:Abs(R×W÷2)→T:S-T→C:S＋T→D If K≤S:Then Goto 3:Else If K≤D:Then Goto 4:Else G-L×A →H:Goto 5IfEnd:IfEnd:IfEnd Lb1 3:If W＞0:Then G+(C-K)2÷2÷R-L×B→H:Goto 5:Else If W ﹤0:Then G-(C-K)2÷2÷R-L×B→H:Goto 5:IfEnd:IfEnd Lb1 4 If W＞0:Then G+(K-D)2÷2÷R-L×A→H:Goto 5:Else If W ﹤0:Then G-(K-D)2÷2÷R-L×A→H:Goto 5:IfEnd:IfEnd Lb1 5 H-0.015×Abs（N）-E→H “H=”:H◢ If J=0:Then Goto 1：Else Goto 6：IfEnd Lb1 6 “QIAN-SHI=”?X Z-X→S:”DU-SHU=”:S◢ “+T,-W=”:H-S◢ Goto 1 二、线路（竖曲线要素子程序） “LI-CHENG=”?K If K≤曲线段终点桩号：Then [[变坡点后段坡度，变坡点前段坡度，变坡点桩号，变坡点高程，曲线半径]]→Mat F:Prog “竖曲线要素矩阵调用程序“:Return：IfEnd 三、竖曲线要素矩阵调用程序 Mat F[1,1] →A: Mat F[1,2] →B: Mat F[1,3] →S: Mat F[1,4] →G: Mat F[1,5] →R: 注意：0.015为道路横坡，根据需要更改。 一、缓和曲线上的点坐标计算 已知：①缓和曲线上任一点离ZH点的长度：l ②圆曲线的半径：R ③缓和曲线的长度：l0 ④转向角系数：K(1或－1) ⑤过ZH点的切线方位角：α ⑥点ZH的坐标：xZ，yZ 计算过程： 说明：当曲线为左转向时，K=1，为右转向时，K=-1，公式中n的取值如下： 当计算第二缓和曲线上的点坐标时，则： l为到点HZ的长度 α为过点HZ的切线方位角再加上180° K值与计算第一缓和曲线时相反 xZ，yZ为点HZ的坐标 切线角计算公式： 二、圆曲线上的点坐标计算 已知：①圆曲线上任一点离ZH点的长度：l ②圆曲线的半径：R ③缓和曲线的长度：l0 ④转向角系数：K(1或－1) ⑤过ZH点的切线方位角：α ⑥点ZH的坐标：xZ，yZ 计算过程： 说明：当曲线为左转向时，K=1，为右转向时，K=-1，公式中n的取值如下： 当只知道HZ点的坐标时，则： l为到点HZ的长度 α为过点HZ的切线方位角再加上180° K值与知道ZH点坐标时相反 xZ，yZ为点HZ的坐标 三、曲线要素计算公式 公式中各符号说明： l——任意点到起点的曲线长度（或缓曲上任意点到缓曲起点的长度）l1——第一缓和曲线长度 l2——第二缓和曲线长度 l0——对应的缓和曲线长度 R——圆曲线半径 R1——曲线起点处的半径 R2——曲线终点处的半径 P1——曲线起点处的曲率 P2——曲线终点处的曲率 α——曲线转角值 四、竖曲线上高程计算 已知：①第一坡度：i1(上坡为“＋”，下坡为“－”) ②第二坡度：i2(上坡为“＋”，下坡为“－”) ③变坡点桩号：SZ ④变坡点高程：HZ ⑤竖曲线的切线长度：T ⑥待求点桩号：S 计算过程： 五、超高缓和过渡段的横坡计算[整理]9860G线元法隧道(坐标正反计算、高程)计算程序
《坐标方位角及距离计算小程序》代码——Access实现
fx-5800p全线坐标正反算带高程计算程序(线元法)
CASIO fx5800p全线高程计算程序
测量坐标计算及高程计算
公路高程计算公式
怎样计算高速公路路线坐标及高程
5800坐标高程计算程序
4800公路中边桩高程计算程序
坐标、高程计算公式
小白自学ExcelVBA道路坐标高程计算程序流程
[整理]fx-5800P坐标高程计算程序交点法.
超高及高程(源程序)
卡西欧5800现场高程计算程序
5800高程计算(含超高)
路面标高计算程序
最新坐标、高程计算公式
Fx5800p市政道路高程计算及放样程序
坐标高程计算公式