5800道路全程序
一、程序功能
本程序由6个主程序、5个次子程序及5个参数子程序组成。主要用于公路测量中坐标正反算，设计任意点高程及横坡计算,桥涵放样，路基开挖口及填方坡脚线放样。程序
坐标计算适应于任何线型.
二、源程序
1.主程序1：一般放样反算程序(①正算坐标、放样点至置仪点方位角及距离；②反算桩号
及距中距离)
程序名:1ZD-XY
Lb1 0:Norm 2↙
F=1：(正反算判别，F=1正算，F=2反算,也可以改F前加？，改F为变量)
Z[1]=90（与路线右边夹角）↙
Prog＂THB＂：F=1=>Goto 1:F=2=>Goto 2↙
Lb1 1:Ｆix 3:＂Ｘ=＂：Locate 6,4,Ｘ◢↙
＂Ｙ＝＂：Locate 6,4,Ｙ◢↙
Ｐrog＂3JS”:Goto 0:↙
Lb1 2:Fix 3:＂ＫＭ＝＂：Locate 6,4,Ｚ◢↙
＂Ｄ＝＂：Locate 6,4,Ｄ◢↙
Ｇoto 0↙
２．主程序２：高程序横坡程序(设计任意点高程及横坡)
程序名：2GC
LbI 0:Norm 2↙
“KM”?Z:?D:↙
Prog”H”:Fix 3:”H=”:Locate 6,4,H◢↙
“I=”: Locate 6,4,I◢↙
Goto 0↙
3.主程序3：极坐放样计算程序(计算放样点至置仪点方位角及距离)
程序名：3JS
X：Y：↙
1268．123→K(置仪点X坐标)
2243．545→L（置仪点Y坐标，都是手工输入,也可以建导线点数据库子程序,个人认为太麻
烦）↙
Y-L→E：X-K→F：Pol(F,E):IF J<0:Then
J+360→J:Int(J)+0.01Int(60Frac(J))+0.006Frac(60Frac(J)) →J:(不习惯小数点后四位为角度显示的，也可以用命令J◢DMS◢来直接显示) Fix 4:” FWJ=”: Locate 6,4,J◢(不习惯小数点后四位为角度显示的，也可以用命令J◢DMS◢来直接显示) ↙
Fix 3:”S=”:Locate 6,4,I◢
4．主程序4：涵洞放样程序（由涵中心桩号计算出各涵角坐标、在主程序3中输入置仪点
坐标后计算放样点至置仪点方位角及距离)
程序名：4JH-XY
LbI 0:Norm 2↙
90→Z[1](涵洞中心桩与右边夹角，手工输入，也可以修改成前面加？后变为变量)
1→F:Prog”THB”:?L: ↙
Z[2]-Z[1] →E:X+Lcos(E) →X:Y+Lsin(E) →Y:Fix 3: ＂Ｘ=＂：Locate 6,4,Ｘ◢
＂Ｙ＝＂：Locate 6,4,Ｙ◢
Ｐrog＂3JS”:Goto 0:↙
5．主程序5:路基开挖边线及填方坡脚线放样程序（输入大概桩号及测量坐标、地面标高计
算出偏移距离、桩号、距中距离、填挖高度）
程序名：5FBX
LbI 0:Norm 2: 18→DimZ:2→F:90→Z[1]:Prog
“THB”:Z:D:”M0”?M:M→Z[4]:D→Z[3]:Prog”6GD”:L→Z[6]:If D<0:Then
0.75-L→D:Goto H:Else L-0.75→D:Goto H:IfEnd↙
LbI H:Prog”H”:H-0.03-Z[4] →Z[5]:Z[6] →L:If Z[5]<0:Then –z[5] →G:Goto W:Else
Z[5] →G:Goto T:↙
LbI W:Prog “W0”:Z[10]+Z[11] →A: If G>A:Then Goto 1:Else If G>Z[10]:Then Goto
2:Else Goto 3:IfEnd: ↙
LbI 1:L+Z[12]+Z[13]+Z[14]+(G-A)×Z[9]+Z[11]×Z[8]+Z[10]×Z[7]:Goto Z:↙LbI 2:L+Z[12]+Z[13]+(G-Z[10])×Z[8]+Z[10]×Z[7]:Goto Z:↙
LbI 3:L+Z[12]+G×Z[7]:Goto z:↙
LbI T:L+0.5→N:If G>Z[17]:Then (N+Z[18]+(G-Z[17])×Z[16]+Z[17]×Z[15])→S:Goto
Z:Else (N+G×Z[15])→S:Goto z:↙
LbI Z:Z[3]→D:Fix 2:Abs(D)-S→T:”L0=”:L Locate 6,4,T◢
＂ＫＭ＝＂：Locate 6,4,Ｚ◢
＂Ｄ＝＂：Locate 6,4,Ｄ◢
“TW=”: Locate 6,4,Z[5]◢
Goto 0↙
6．主程序6：路基标准半幅宽度计算程序(对于设计有加宽渐变的有用，如路基宽度无变化，
则把此程序直接输入半幅宽度值至L)
程序名：6GD
Prog “G0”Z-C→E:(B-A)×E/S+A→L:L:↙
7．坐标计算次程序（ＴＨＢ）
程序名：THB
18→DimZ：＂ＫＭ＂?Z：Prog ＂Ｘ０＂↙
1÷P→Ｃ： (P-R)÷(2HPR) →Ｓ：180÷π→Ｅ：Ｆ=1=>Goto 1：Ｆ=2=>Goto 2←┘↙Lbl 1：?Ｄ： Abs(Ｚ－O) →W：Prog "Ａ"：X：Ｙ：Ｇoto 3↙LbI 2:Ｘ：Ｙ：Ｘ→I：Ｙ→J：Prog ＂Ｂ＂：Ｏ＋Ｗ→Z：Ｄ→Ｄ：Ｇoto 3↙
LbI 3:IF Ｆ＝１Then Ｘ：Ｙ：Else Ｚ：Ｄ↙
8. 正算子程序(Ａ)
程序名：A
0.1184634425→A： 0.2393143352→B：0.2844444444→N 0.046910077→K：
0.2307653449→L：0.5→M： U+W(Acos(G+QEKW(C+KWＳ))+Bcos(G+QELW(C+LW
Ｓ))+Ncos(G+QEMW(C+MW
Ｓ))+Bcos(G+QE(1-L)W(C+(1-L)WS))+Acos(G+QE(1-K)W(C+(1-K)WS))) →X：↙
V+W(Asin(G+QEKW(C+KWＳ))+Bsin(G+QELW(C+LW
Ｓ))+Nsin(G+QEMW(C+MWS))+Bsin(G+QE(1-L)W(C+(1-L)WS))+Asin(G+QE(1-K)W(C+(1-K)W
Ｓ))) →Y：↙
G+QEW(C+WＳ)+Ｚ［１］→Ｚ［２］：X+Ｄcos（Ｚ［２］）→Ｘ： Y+Ｄsin（Ｚ［２］）→
Ｙ↙
9. 反算子程序(Ｂ)
程序名：B
G-90→Ｔ： Abs((Y-V)cos(T)-(X-U)sin(T)) →Ｗ：0→Ｄ：Lbl 0：Prog "Ａ"： T+QEW(C+W Ｓ) →Ｌ： (J-Y)cos(L)-(I-X)sin(L)→Ｄ：IF Abs(Ｄ)<0.01:Then Goto1：Else W+Ｄ
→W:Goto 0←┘↙
Lbl 1:0→D：Prog "Ａ"：(J-Y)÷sin(Ｚ［２］) →D:↙
10．高程计算子程序（H）
程序名：H
Prog “S0”:R:T:C:G:I:C-T→F:Z-F→L:C+T→E:G-TI→Q:If T=O:Then Q+LI→H:Goto 0:Else If Z LbI 0:H:If D=0:Then Goto I:Else Prog “I”:H+V→H:Goto I:LbI I:H:I:↙ 11.高程超高计算程序（I） 程序名：I Prog”I0”:↙ W=1=> Goto 0:W=2=>Goto 1: ↙ LbI 0:If L=0:Then Abs(D)×M→V:Goto 2:Else Abs(D)×((N-M)×(Z-C)÷L+M)→V:Goto 2:IfEnd: ↙ LbI 1:If L=0:Then Abs(D)×M→V:Goto 2:Else Abs(D)×(((3((Z-C)÷L)2 -2((Z-C)÷L)∧(3))×(N-M))+M)→V:Goto 2:IfEnd:↙ LbI 2:Abs(D)→E:V÷E→I:I(E-K)→V:↙ 12．数据子程序(附后示例) ①程序名：X0（坐标计算要素程序） If Z≥25900 And Z≤26615.555:Then 25900→O:11587.421→U:1847.983→V:101。 09’23.1”→G:715.555→H:1×1045→P: 1×1045→R: 0→Q:Goto 0: IfEnd:↙ If Z≥26615.555 And Z≤26915.555:Then 26615.555→O:11448.97→U:2550.016→V:101。 09’23.1”→G:300→H:1×1045→P: 1800→R:1→Q:Goto 0:IfEnd:↙ If Z≥26915.555 And Z≤27316.952:Then 26915.555→O:11382.792→U:2842.531→V:105。 55’51.83”→G:401.396→H:1800→P: 1800→R:1→Q:Goto 0:IfEnd:↙ ……↙ LbI 0:O:U:V:G:H:P:R:Q: ↙ 程序字母说明：O-线元起点桩号；U-起点X坐标；V-起点Y坐标；G-线元起点桩号切线方位角; Ｈ－线元长度，Ｐ－线元起点曲率半径，Ｒ－线元终点曲率半径,Q－线元判别(以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时， Q=1；当线元为直线时，Q=0)。 ②程序名：SO(高程竖曲线要素子程序) If Z≥25900 And Z≤26157.8:Then -200000→R:157.8→T:26000→C:37.1→G:0.01222→I:Goto 0:IfEnd: ↙ If Z≥26157.8 And Z≤27421.915:Then 300000→R:171.915→T:27250→C:32.65→G:-0.00356→I:Goto 0:IfEnd:↙ If Z≥27421.915 And Z≤27889.029:Then 1000000→R:134.029→T:27755→C:36.64→G:0.0079→I:Goto 0:IfEnd:↙ ……↙ LbI 0:R:T:G:C:I: ↙ 程序字母说明：C-竖曲线交点桩号；G-交点桩号高程（未竖曲线调整的）；T-竖曲线切线长； R-竖曲线半径（分正负：凹曲线为正，凸曲线为负）；I－竖曲线前纵坡 ③程序名：I0（高程超高参数子程序） 1.5→K:2→W:If D<0:Then Goto L:Else Goto R:IfEnd:↙ LbI L:If Z≥25900 And Z≤26615.555:Then 0→L:-0.02→M:Goto 0:IfEnd:↙ If Z≥26615.555 And Z≤26735.555:Then 120→L:26615.555→C:-0.02→M:0.02→N:Got o 0:IfEnd: ↙ If Z≥26735.555 And Z≤26795.555:Then 60→L:26735.555→C:0.02→M:0.04→N:Goto 0:IfEnd↙ If Z≥26795.555 And Z≤27436.951:Then 0→L:0.04→M: Goto 0:IfEnd↙ ……↙ LbI R: If Z≥25900 And Z≤26735.555:Then 0→L:-0.02→M: Goto 0:IfEnd↙ If Z≥26735.555 And Z≤26795.555:Then 60→L:26735.555→C:-0.02→M:-0.04→N:Goto 0:IfEnd↙ If Z≥26795.555 And Z≤27436.951:Then 0→L:-0.04→M: Goto 0:IfEnd ……↙ LbI 0 K:W:L:C:M:N: ↙ 字程序中字母表示说明： K－中央分隔带半幅宽（中桩标高至中桩中线距离，无分隔带则为0）；W－超高方式参数（W=1为一般直线方式超高，W=2为三次抛物线方式超高）;L－超高渐变段距离（不是渐变段则输入0）；C-超高渐变段起点桩号（不是渐变段，无需输入）；M－超高段起点横坡，N-超高 渐变段终点横坡（不是渐变段，无需输入）。 ④程序名：GO(路基标准半幅宽度参数子程序) If D<0:Then Goto L:Else Goto R:IfEnd: ↙ LbI L:If Z≥25900 And Z≤28110.727:Then Z→C:17.25→A:Goto 0:IfEnd:↙ If Z≥28110.727 And Z≤28200.727:Then 28110.727→C:17.25→A:21.25→B:90→S:Goto 0:IfEnd↙ ……↙ LbI R: If Z≥25900 And Z≤27927.478:Then Z→C:17.25→A: Goto 0:IfEnd↙ If Z≥27927.478 And Z≤28172:Then 27927.478→C:17.25→A:27.031→B:244.522→S:Goto 0:IfEnd↙ ……↙ LbI 0 C:A:B:S: ↙ 字程序中字母表示说明： C－宽度渐变段起点桩号（不是渐变段C=Z）；A－宽度渐变段起点宽度; B－宽度渐变段终点宽度（不是渐变段，无需输入）；S-宽度渐变段距离（不是渐变段，无需输入）。 ⑤程序名：WO(路基填挖边坡参数子程序) 0.75→Z[7]:1→Z[8]:1.25→Z[9]:10→Z[10]:10→Z[11]:2.6→Z[12]:2→Z[13]:2→Z[14]: ↙ 1.5→Z[15]:1.75→Z[16]:8→Z[17]:2→Z[18]↙ 字程序中字母表示说明： Z[7]－挖方第一阶边坡边率；Z[8]－挖方第二阶边坡边率；Z[9]－挖方第三阶边坡边率；Z[10]－挖方第一阶高度；Z[11]－挖方第二阶高度；Z[12]－挖方路基碎落台及水沟宽度；Z[13]－挖方第二台阶平台宽度；Z[14]－挖方第三台阶平台宽度； Z[15]－填方第一阶边坡边率；Z[16]－填方第二阶边坡边率; Z[17]－填方第一阶高度； Z[18]-填方第二台阶平台宽度。（注：本程序只做出挖方三个台阶，填方二个台阶，如需增加，先需在程序THB和5FBZ程序中增加变量。再如一标段有坡率及坡高不一样时，可以照 其它参数程序一样，用判别语句。） 三、使用说明 1、规定 (1) Z“KM”为所求点桩号，反算时为输入大概桩号 (2) 当所求点位于中线时，Ｄ=0；当位于中线左铡时，Ｄ取负值；当位于中 线中线右 侧时，Ｄ取正值。 (3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次 代替。 (4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径 均等于圆 弧的半径。 (5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲 率半 径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。 (6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设 计规定的 值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。 2、输入与显示说明 输入部分： 1. Ｆ＝１：ＺＤ => XY 2.Ｆ＝２： XY = > ＺＤ Ｆ ? 选择计算方式，输入1表示进行由里程、边距计算坐标；输入2表示 由坐标反算里程和边距。Ｚ［１］为线元边桩与中线右夹角。 ＫＭ？正算时所求点的里程（反算输入大概桩号） Ｄ？正算时所求点距中线的边距(左侧取负，值右侧取正值，在中线上取 零) X0 ？反算时所求点的X坐标（放样程序中实测X坐标） Y0 ？反算时所求点的Y坐标（放样程序中实测Y坐标） 显示部分： L? 涵洞放样程序中涵距中心桩号前后距离,前为正,后为负. M0? 边坡放样程序中原地面标高输入 X=*** 正算时，计算得出的所求点的X坐标 Y=*** 正算时，计算得出的所求点的Y坐标 FWJ=*** 正算时，计算得出的所求点的至置仪点方位角（在3JS程序中应输入X、Y坐 标值于K、L） S=*** 正算时，计算得出的所求点的至置仪点距离 ＫＭ=*** 反算时，计算得出的所求点的里程 Ｄ=*** 反算时，计算得出的所求点的边距 H=*** 所求点位置设计路面顶标高 I=*** 所求点位置设计路面横坡 L0=*** 边桩放样程序中实测点至设计边坡点距离，正向内移， 直线段边坡超欠挖检查及开口线放样程序程序名：BPFY 程序：A“QX”：B“QY”：C“QZ”：E“ZX”：F“ZY”：G“ZZ”：P“PB”：K“XP”：L“YP”：M“ZP”：Fixm：Pol(E-A，F-B：J≤0=>J=J 360⊿D=Abs((K-A)*CosJ (L-B)*SinJ)：S=Abs((L-F)*CosJ-(K-E)*SinJ)：X=A D*CosJ：Y=B D*SinJ：Z=(G-C)/I*D C：W=AbS(M-Z)：V=S/P-W：V<0=>O“Hcw”=V◢⊿V≥0=>O“Hqw”=V◢⊿U=W*P-S：U<0=>N“Scw”=U◢⊿U≥0=>N“Sqw”=U◢⊿“END”说明：A、B、C为边坡底线的起点，显示为: QX?QY?QZ? 依次输入地线起点的X,Y,Z；E、F、G为边坡底线的终点，显示为: ZX?ZY?ZZ? 依次输入地线终点的X,Y,Z；P为坡比，无正负条件；K、L、M 为测量点，显示为: XP?YP?PZ? 依次输入测量点的X,Y,Z；中间计算测量点至边坡底线的垂足点，保存变量为X、Y、Z。O为边坡的高程超欠挖值，显示为Hcw或Hqw，分别表示超挖或欠挖，即测点高程设计高程减实测高程；N为边坡的距离超欠挖值，显示为Scw或Sqw，分别表示超挖或欠挖，按高差计算的设计距离减实测点到垂足点距离。圆弧段边坡超欠挖检查及开口线放样程序程序名：YFFY 程序：A“YX”：B“YY”：X“XP”：Y“YP”：Z“ZP”：Fixm：Pol(A-X，B-Y)：T=Abs(I-R)：W=Abs(H-Z)：S=W*P-T：V=T/P-W：S<0=>M“Scw”=S◢⊿S≥0=>M “Sqw”=S◢⊿V<0=>K“Hcw”=V◢⊿V≥0=>K“Hqw”=V◢⊿“END”说明：A、B为圆心坐标，显示为YX？YY？，依次输入圆心点的X、Y；X、Y、Z为测点坐标，显示为XP？YP？ZP？依次输入测点的X、Y，Z；R为到边坡底线的距离（半径）；H为边坡底线的高程（马道或平台的高程）；I不用管它，直接按“EXE”，显示下一个输入；P为坡比，无须输入正负号；M为计算的边坡距离超欠值，显示为Scw或Sqw，分别表示超挖或欠挖，距离超欠值为：设计距离减测量距离；K为计算的边坡高程超欠值，显示为Hcw或Hqw，分别表示超挖或欠挖，高程超欠值为：设计高程减测量高程。从网上弄来的，自己没有用过，你用的时候最好多做检核。 边坡测量（由坐标反算桩号及填挖高度方法的运用） 在施工测量放样中，最常用的就是正算放样（已知该点坐标进行定点放样），但在某些特殊情况下要进行的放样工作，这种办法就显得黔驴技穷了，比如在南方的山区道路施工中，往往会有高填高挖路段，有的高达几十米，那么这些段落的填土边线及挖方开口线放样就给我们施工测量带来了麻烦，如果采用常规的办法一般会有以下几个步骤： 1、首先放出中线 2、然后对填挖断面进行测量 3、再估算填挖边线距中桩的距离 4、试放出该点位置 5、对该点进行高程测量，比较该点实测高程与设计高程的差距，重新计算距中桩距离 6、重复4、5步，试放该点，测量高程，至到与设计高程相符为止 线元法万能曲线正反算简介 我的线元法是把线形分为直线和曲线，直线就不用说了，起止点桩号，坐标和方位角就可以算了；曲线最基本的组合：是由一段缓和曲线+一段圆曲线组成，任意复杂的曲线都可以分解成缓和曲线+圆曲线或者其中之一就可以。 分析最复杂的曲线可以看到: 一般复杂线形由Ls1 ,R1,Ls2, R2组成，相邻的Ls1+R1,一般满足A*A=Ls1*R1,这就是一个线元法单元，即使不满足也可以作为一个线元： 当Ls1= Ls2，且R1= R2时，为单曲线 当Ls1≠ Ls2，或者R1≠R2时，为复合曲线 当Ls1= Ls2=0时，线性为圆曲线， 当圆曲线长度为0时，线性为缓和曲线+缓和曲线， 当A*A≠Ls1*R1时，为卵形曲线，需要计算虚拟起点坐标 综合以上线形，本程序正反算计算全部可以处理。结合目前流行的线元法，本程序也可以，分为缓和曲线和圆曲线录入，方法是一样的，所不同的是起点要注意，复杂曲线，是两边向中间定义数据库，缓和曲线永远是ZH点或HZ点为起点。 曲线要素说明（有9个）： 1、起点桩号：（一般为ZH点或HZ点，或ZY点或YZ点，或者卵形公切点GQ） 2~3、起点坐标：（X，Y） 4、起点方位角：FWJ 114°15′24.33″写成：114.152433 5、线性特征：直线，左偏，右偏；三个选一个 6、终点桩号：如果起点为ZH点，终点一边为YH点，QZ点，HY点，都可以，一般为YH点，缓和曲线+圆曲线。如果缓和曲线Ls=0，就是YZ点；大小不一定按路线顺序，如果起点为HZ点，终点根据缓和曲线+圆曲线的特点，和上个线元对接上就可以了。 7、缓和曲线长度Ls： 8、圆曲线半径R: 9、回旋参数A: 一般满足A*A=Ls1*R1，不满足条件的是卵形曲线。 可以处理任意数量断链。 操作流程：1、先编辑线元数据，保存后推出。 2、如果有线元断链的输以下线元断链数据 3、打开线元万能曲线计算单点计算就可以了。 目前，已有一个例子文件在里面，在安装文件目录下“ \dmfx4.0\demo\左线”，有个CAD文件，里面有校核数据，可以看到本软件处理的逐桩表和要素表，可以验证软件的数据，任意数据坐标反算可以得到桩号和距中，任意输入桩号和距中可以正算得坐标。 授权版用户，可以通过运行交点文件编辑，保存后，退出；打开线元法数据编辑，浏览正在使用的主项目文件，就可以看到一个线元数据，点击这个文件确定，保存退出。就完成交点法数据转换线元法数据过程。 阅览室馆友我的图书馆 kaixin100 | | 分享 U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos( G+QEFW (C+FWD))+Acos(G+QEMW(C+MWD)))→X： V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+ QEFW(C+FWD))+Asin(G+QEMW(C+MWD))) →Y： G+QEW(C+WD)+90→F：X+Zcos（F）→X：Y+Zsin（F）→Y 反算子程序SUB2 G-90→T (Y-V)cosT-(X-U)sin(T) →W Abs(W)→W:0→Z Lbl6：Prog "SUB1" T+QEW(C+WD) →L：(J-Y)cos（L）-(I-X)sin（L）→Z IF Abs（Z）<1E-6：Then0→Z:Prog "SUB1":(J-Y)÷sin（F）→Z:Else W+Z→W：Goto6:IfEnd 数据库子程序SUB0 Goto 1(线元可输入多条，分离式可在前多加一位，匝道一样。例：左幅为K129+500,右幅输线元参数里程为1129+500，其他不变，前面 1为任意数字，计算机便于区分) Lbl 1：IF S<线元终点里程：Then@@@→O(线元起点里程) ：@@@ →U(线元起点X坐标)：@@@→V(线元起点Y坐标)：@@@→G(线元起点计算方位角)：@@@→P(线元起点半径)：@@@→R(线元止点半径)：@@@→H(线元长度)：@@@→Q(线元左、右偏标志，左偏-1，右偏1，直线为0)：Return：IfEnd IF S<线元终点里程：Then@@@→O(线元起点里程) ：@@@ →U(线元起点X坐标)：@@@→V(线元起点Y坐标)：@@@→G线元(起点计算方位角)：@@@→P(线元起点半径)：@@@→R(线元止点半径)：@@@→H(线元长度)：@@@→Q(线元左、右偏标志，左偏-1，右偏1，直线为0)：Return：IfEnd 。。。。。。。。。。。。。。。。。。。。。。。。。。。 一程序功能 本程序由一个主程序(ZBJS)和3个子程——正算子程序(SUB1)、反 算子程序( SUB2) 、数据库子程序（SUB0）构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线 元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲 率半径、止点曲 率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标 进行正反算。另 外也可以将本程序中核心算法部分的两个子程序移植到其它相关的 程序中，用于对曲 线任意里程中边桩坐标进行正反算。本程序也可以在CASIO fx-4500P计算器及CASIO fx-4850P计算器上运行。 二、使用说明 道路中边桩坐标放样正反算CASIO fx-5800P程序（全线贯通） 编辑 | 删除 | 权限设置 | 更多▼ 设置置顶推荐日志转到私密记事本 转载自王中伟转载于2009年08月12日 17:34 阅读(1) 评论(0) 分类：技术交流权限: 公开 一、前言 本程序是《CASIO fx-5800P计算与道路坐标放样计算》中道路坐标放样计算程序的升级改进版本。原道路坐标放样计算程序只基于道路的单个基本型曲线，有效计算范围仅包括平曲线部分和前后的两条直线段，使用时需要输入平曲线设计参数，无坐标反算桩号功能。 改进后的程序名称为：道路中边桩坐标放样正反算程序（全线贯通），增加了可实现全线贯通的数据库功能和坐标反算桩号功能，主要是： 1．使用道路平面数据库子程序，可将一段或若干段道路的交点法格式平面参数（可容易从直线、曲线及转角表中获得）以数据库子程序形式输入计算器，程序在计算时省却了输入原始数据的麻烦； 2．坐标正算方面，输入桩号即可进行道路的中、边桩坐标计算，若输入了测站坐标，还可同时计算全站仪极坐标放样数据（拨角和平距）； 3．坐标反算方面，输入平面坐标，即可计算对应的桩号和距中距离（含左右信息）； 4．对于存在断链的道路，可分段分别编写数据库子程序，然后在主程序中添加一个路段选择的功能即可实现（可参照立交匝道程序中匝道的选择）。 程序的特点： 1．可进行中桩坐标的正、反算，程序代码简洁，便于阅读和改写； 2．主程序通过调用数据库子程序，省却了使用时输入平面参数的繁琐； 3．使用数据库子程序，换项目只需改写数据库子程序，程序通用性强。 二、道路示例项目基本资料 基本资料同《CASIO fx-5800P计算与道路坐标放样计算》第6章HY高速公路第2合同段（合同段起止桩号：K4+800~K9+600）。这里摘取直线、曲线及转角表资料如下（若图片不清晰，请参见参见教材P161附录1）： . 颜色文字Casio5800交点法与线元法（积木法）匝道坐标正反算放样程序（XUFENG 2011.2.14） 本人一直以来想找一个交点法与线元法相结合的坐标正反算程序，在网上找了很久很久，没能找到一个较为满意的，有幸在测量空间看到大歪哥的《Casio5800交点法程序》与《线元法（积木法）匝道坐标正反算放样程序》，根据歪哥意见“需要的自行修改结合XY框架自己修改为数据库反算程序等”，本人不才，采用最笨的办法将两个程序综合了一下，使之能既能进行交点法正反算，又能进行线元法正反算。在此特别感谢大歪哥！ 将程序发上来，愿与大家一同交流学习欢迎大家吐口水，只要能进步就行！ 程序由一个主程序ZBZFS和8个子程序（JS、XY-A、XY-B、JDYS、1、2、3、4）构成，运行时只需运行主程序即可！ 本程序适用于单交点对称型、不对称型、无缓和曲线单圆曲线型一个交点范围内（含交点前后有直线段时）的曲线要素核对和坐标正反算，手工输入要素，对设计图纸的“直线、曲线转角表”中交点数据进行复核验证，并能对单一线元进行坐标正反算。 1主程序名：ZBZFS（功能：进入计算主程序） 65→Dimz↙ Deg:Fix 3↙ ＂1.JD ZFS 2. ZHADAO ZFS＂? I: I→Z[61]: ＂1.ZHONG SHU JS 2. JS＂? I ↙ If I=1: Then Goto1: Else Goto2:IfEnd↙ LbI 1 :If Z[61]=1: Then Prog＂JDYS＂:Else Cls:＂K0＂?A:＂KN＂?L :＂X0＂?U :＂Y0＂?V :＂F0＂?W :＂R0＂?P :＂RN＂?Q:＂ZX:-1,+1,0＂?G:IfEnd↙ LbI 2 :Prog＂JS＂ 2子程序名：JS（功能：选择正算或反算模式） Cls:＂XC＂?H:＂YC＂?Z↙ 1fx－5800P计算器编程 缓和曲线程序 14→DimZ ：“ZHK”：?K：?R ：? A：?L：? T：? F：“Y＋1，Z－1”：?P：“JDX”：?Q：“JDY”：?W：180÷（πR）→Z：0.5L－L^3÷(240R2)→Z[8]：“LZ”：ZL÷2→B：A÷Z＋L→Z[13]◢ “ZHK=”：K◢ “HYK=”：K+L→Z[1]◢ “YHK=”：K+ Z[13] －L→Z[2]◢ “HZK=”：K+ Z[13]→Z[3]◢ “ZHX=”：Q+Tcos(F+180)→U◢ “ZHY=”：W+Tsin(F+180)→V◢ “HZX=”：Q+ Tcos(F+PA)→Z[6]◢ “HZY=”：W+ Tsin (F+PA)→Z[7]◢ Lb1 0：“CDZH”：?M If M≤K：Then Goto 6：Else If M≤Z[1] ：Then Goto 1：Else If M≤Z[2] ：Then Goto 2：Else If M≥Z[3] ：Then Goto 7：Else If M≥Z[2] ：Then Goto 5：If End：If End：If End：If End：If End Lb1 1：M－K→G：√（（G－G^(5)÷(40L2R2)＋G^(9)÷(3456R^(4)L^(4)))2＋(G^(3)÷(6RL)－G^(7)÷(336R^(3)L^(3)) ＋G^(11)÷(42240R^(5)L^(5)))2）→D F＋PZG2÷(6L)→C：F＋PBG2÷L2→H：Goto 3 Lb1 2：M－Z[1]→G：R（1－cos（B＋ZG））＋L2÷（24R）→Z[9]：√（Z[9]2＋（Rsin（B＋ZG）＋Z[8]）2）→D F＋Ptg-1(Z[9]÷(Rsin(B＋ZG)＋Z[8]))→C：F＋P(B＋ZG)→H Lb1 3 “X=”:U＋Dcos(C)→X◢ “Y=”:V＋Dsin(C)→Y◢ Lb1 B：“ZBJS1，FY－1”：?J If J>0：Then Goto 4：Else If J〈0：Then Goto A：If End：If End Lb1 4：“ZB”：?S “ZBX”:X＋Scos（H－90）◢ “ZBY”:Y＋Ssin（H－90）◢ “YB”：?N “YBX”:X＋Ncos（H＋90）◢ 一、前言本程序是《CASIO fx-5800P计算与道路坐标放样计算》中道路坐标放样计算程序的升级改进版本。原道路坐标放样计算程序只基于道路的单个基本型曲线，有效计算范围仅包括平曲线部分和前后的两条直线段，使用时需要输入平曲线设计参数，无坐标反算桩号功能。改进后的程序名称为：道路中边桩坐标放样正反算程序（全线贯通），增加了可实现全线贯通的数据库功能和坐标反算桩号功能，主要是： 1．使用道路平面数据库子程序，可将一段或若干段道路的交点法格式平面参数（可容易从直线、曲线及转角表中获得）以数据库子程序形式输入计算器，程序在计算时省却了输入原始数据的麻烦； 2．坐标正算方面，输入桩号即可进行道路的中、边桩坐标计算，若输入了测站坐标，还可同时计算全站仪极坐标放样数据（拨角和平距）； 3．坐标反算方面，输入平面坐标，即可计算对应的桩号和距中距离（含左右信息）； 4．对于存在断链的道路，可分段分别编写数据库子程序，然后在主程序中添加一个路段选择的功能即可实现（可参照立交匝道程序中匝道的选择）。程序的特点： 1．可进行中桩坐标的正、反算，程序代码简洁，便于阅读和改写； 2．主程序通过调用数据库子程序，省却了使用时输入平面参数的繁琐； 3．使用数据库子程序，换项目只需改写数据库子程序，程序通用性强。二、道路示例项目基本资料基本资料同《CASIO fx-5800P计算与道路坐标放样计算》第6章HY高速公路第2合同段（合同段起止桩号： K4+800~K9+600）。这里摘取直线、曲线及转角表资料如下. . . 三、程序代码 . . . . . 5800全线任意坐标计算程序 1. 正算主程序(ZHCX) （不运行） 8→DimZ 1÷P→Z[4 ]：(P-R)÷(2HPR)→D: 180÷π→E “Z=”?Z:”YJJ=”?A:Abs(S-O)→W 0.26→Z[1 ]: 0.74→B: 0.02→K: 0.82→Z[3 ]: 1-Z[3 ]→F:1-K→Z[2 ] U+W(Z[1 ]cos(G+QEKW(Z[4 ]+KWD))+Bcos(G+Z[3 ]QEW(Z[4 ]+ Z[3 ]WD))+Bcos(G+QEFW (Z[4 ]+FWD))+ Z[1 ]cos(G+ Z[2 ]QEW(Z[4 ]+ Z[2 ]WD)))→X: V+W(Z[1 ] sin (G+QEKW(Z[4 ]+KWD))+B sin(G+ Z[3 ]QEW(Z[4 ]+ Z[3 ]WD))+B sin(G+QEFW (Z[4 ]+FWD))+ Z[1 ] sin(G+ Z[2 ]QEW(Z[4 ]+ Z[2 ]WD)))→Y： G+QEW(Z[4 ]+WD)→F：X+Zcos(F+A)→X：Y+Zsin(F+A)→Y:If F≧360:Then F-360→F:IfEnd ”X=”:X→X◢ ”Y=”:Y→Y◢ If F﹤0:Then F+360→F:IfEnd ”QX FWJ=”:F▼DMS◢ “C=1=>XX: C=2=>XZ”: ”C=”?C: ”QHJU=”?L: If C=1：Then Goto 1:Else Goto 2: IfEnd 可以计算斜交斜做或斜交正做的桥涵坐标 Lbi 1 X+L cos(F)→X：Y+Lsin(F)→Y: Goto 3 Lbi 2 X+L cos(F+A-90)→X：Y+Lsin(F+A-90)→Y: Goto 3 Lbi 3 “QH-X=”: X →X◢ “QH-Y=”: Y →Y◢ Prog “FY” 2 . 参数子程序（直接运行） M(主线) 一条线路一个名称 “S=”?S If S≦线元终点：Then 线元起点X值→U: 线元起点Y值→V:线元起点切线方位角→G:线元起点桩号→O:线元长度→H：线元起点半径→P：线元终点半径→R：（左偏-1，或右偏 1）→Q:Goto 1:IfEnd … … If S≦线元终点：Then 线元起点X值→U: 线元起点Y值→V:线元起点切线方位角→G:线元起点桩号→O:线元长度→H：线元起点半径→P：线元终点半径→R：（左偏-1，或右偏 1）→Q:Goto 1:IfEnd Lbi 1 Prog “ZBJS” 3. 放样程序（FY）（不运行） “X0=”?M:“Y0=”？N Pol((X-M, Y-N) FX5800计算器测量程序集版 一、程序功能 主要功能：采用交点法方式计算多条线路坐标正反算，可算任意复杂线型及立交匝道，包括C型，S型、卵型、回头曲线等；极坐标放样，全线路基边坡开挖口及坡脚放样计算、路基任意点设计高程、横坡、设计半幅宽度.隧道欠超挖放样计算。 新版本优化：1、优化程序语句、2、以复数形式输入变量及做数据库，取消原矩阵数据库;3、修改隧道超欠挖程序为通用形，不受圆心个数限制、4、新增测量资料表计算 二、源程序（绿色为程序名；蓝色为输入计算器内容）紫色为新版改动处（可以根据自己标段情况用相关主程序及子程序，再在0程序中汇总）0.汇总程序（1、坐标计算放样程序(1XY、A、AB、HX、JS、DX、QX、F、XY、X1)；2、坐标反算程序（2ZD、A、B、AB、HX、QX、F、ZD、X1)；3、高程计算查阅程序（3GC、H、I、QX、S1、I1）；4、路基半幅标准宽度查阅程序（4GD、C、QX、G1）；5、路基边坡及开挖口放样程序（5BP、 A、B、AB、HX、H、I、C、JS、DX、QX、F、ZD、X1、S1、I1、G1、W1）； 6、路基标准距离放样（6FM、A、AB、HX、H、I、JS、DX、QX、F、XY、X1、S1、I1)； 7、桥梁锥坡计算放样程序（7ZP、A、AB、HX、C、JS、DX、QX、F、XY、X1、G1)； 8、极坐标计算程序（8JS、JS、DS）； 9、隧道超欠挖计算程序（A、B、AB、HX、H、I、QX、S、SD、F、ZD、X1、S1、I1、SD1）运行后按1~9数子约半秒，则选择1至9的程序，返回时，在桩号输入-1,返回选择计算类型。输入-2,返回选择线路。 程序名：0（数子0） ClrMat:ClrVar:12→DimZ：Norm 2:Do:＂(XY=1，ZD=2 ，GC=3，GD=4，BP=5，FM=6，ZP=7，JS=8，SD=9)===>QING AN 1-9＂：Getkey→Z[3]:While Z[3]=35:Prog＂1XY＂:WhileEnd:While Z[3]=36:Prog＂2ZD＂：WhileEnd: While Z[3]=37:Prog＂3GC＂：WhileEnd: While Z[3]=21:Prog＂4GD＂：WhileEnd: While Z[3]=22:Prog＂5BP＂：WhileEnd: While Z[3]=23:Prog＂6FM＂：WhileEnd: While Z[3]=31:Prog＂7ZP＂：WhileEnd: While Z[3]=32:Prog＂8JS＂： 任意里程中边桩坐标正反算(CASIO fx-4800P计算器)程序 时间:2007-12-18 13:13:15 来源:工测员网作者:未知 1. 加编数据库作为主程序 , 计算中不必逐项输入 " 线元要素 ", 提高运算速度，避免现场忙中出错 2. 将原来的主程序并入数据库 3. 计算直观 , 人性化 4. 正算直接输入里程和边距 , 反算输入近似里程便可 5. 增加了“ 计算点与测站点” 的距离和方位角计算语句，方便直接放样 6. 愿收获与大家共享 7. 核心计算程序摘自“yshf” 1. 正算子程序 (SUB1) A=0.1739274226 ： B=0.3260725774 ： K=0.0694318442 ： L=0.3300094782 ：F=1-L ： M=1-K ： X=U+W(Acos(G+57.2958QKW(1/P+KWD))+Bcos(G+57.2958QLW(1/P+LWD))+Bcos(G+ 57.2958QFW (1/P+FWD))+Acos(G+57.2958QMW(1/P+MWD))) ： Y=V+W(Asin(G+57.2958QKW(1/P+KWD))+Bsin(G+ 57.2958QLW(1/P+LWD))+Bsin(G+57.2958QFW(1/P+FWD))+Asin(G+57.2958QMW(1/ P+MWD))) ： F=G+57.2958QW(1/P+ WD)+90 ： X=X+ZcosF ： Y=Y+ZsinF 2. 反算子程序 (SUB2) T=G-90 ： W=Abs((Y-V)cosT-(X-U)sinT) ： Z=0 ： Lbl 0 ： Prog "SUB1" ：L=T+57.2958QW(1/P+ WD) ： Z=(J-Y)cosL-(I-X)sinL ： AbsZ<1E-6=>Goto1 ：≠>W=W+Z ： Goto 0Δ←┘ Lbl 1 ： Z=0 ： Prog "SUB1" ： Z=(J-Y)÷sinF 3、 . 增设数据库程序（ＳＪＫ主程序） Lb1 4 ： "1.SZ => XY" ： "2.XY => SZ" ：｛ NS ｝：S ∠ 下一线元起点里程＝＞O = 本线元起点里程： U= 本线元起点 X ： V= 本线元起点 Y ： G= 本线元起算方位角： H= 本线元长度： P= 起点曲率半径： R= 终点曲率半径： Q=0 或 1 、 -1 ： Goto0Δ←┘( 第一线元数据要素) S ∠ 下一线元起点里程＝＞O = 本线元起点里程： U= 本线元起点 X ： V= 本线元起点 Y ： G= 本线元起算方位角： H= 本线元长度： P= 起点曲率半径： R= 终点曲率半径： Q=0 或 1 、 -1 ： Goto0Δ←┘( 第二线元数据要素) S ∠ 下一线元起点里程＝＞O = 本线元起点里程： U= 本线元起点 X ： V= 本线元起点 Y ： G= 本线元起算方位角： H= 本线元长度： P= 起点曲率半径： R= 终点曲率半径： Q=0 或 1 、 -1 ： Goto0Δ←┘( 第三线元数据要素) 5800-平曲线程序 程序目的：依平曲线要素计算直线、圆曲线、缓和曲线的任意中桩、左、右桩坐标。 程序说明：K0：起始桩号 X0：起始X坐标 Y0：起始Y坐标 ALF：起始方位角 R：半径 LS：缓和曲线长 N：曲线左转N=1，右转N=2 K：待求桩号 LL、LR：左、右桩距离 Q：左、右桩与中线斜交角求得XZ、YZ、XL、YL、XR、YR分别为中桩、左、右桩坐标。 一、直线段 文件名：ZX (COMP) “K0=”? L：“X0=”? O：“Y0=”? P：“ALF=”? W： Lbl 0：“K=”?：O+(K-L)cosW→X： P+(K-L)sinW→Y：“XZ=”：X◢“YZ=”：Y◢“LL=”?B：“XJJ=”?Q： X-Bcos(W +Q) →E：Y-Bsin(W+Q) →F：“XL=”：E◢“YL=”：F◢“RL=”？C： X+Ccos(W+Q) →I：Y+Csin(W+Q) →J：“XR=”：I ◢“YR=”：J◢Goto 0 注：在程序执行过程中，赋给的要素变数的值被固定不变，可对变数（K、LL、LR）赋予不同值，迅速求得所需坐标。 二、圆曲线段 文件名：YQX (COMP) “K0=”？ L：“X0=”？ O：“Y0=”？ P：“ZH—JD FWJ=”？ G：“LS=”？L：“R=”？R：G+L÷2÷R×180÷π→W：Lbl 0：“K=”？K：(-1)^N(K-L)÷R×180÷π→J： 2Rsin((-1)^N J÷2)→D：O+Dcos(W+J÷2)→X：P+Dsin(W+J÷2)→Y：“XZ=”：X◢“YZ=”：Y◢“LL=”？B：“XJJ=”？Q：X-Bcos(W+J+Q)：Y-Bsin(W+J+Q) →F：“XL=”：E◢“YL =”：F◢ “RL=”?C：“XJJ=”?Q：X+Ccos(W+J+Q) →U：Y+Csin(W+J+Q) →V：“XR=”：U◢“YR=”：V◢ Goto 0 注：公式里面W为HY点方位角，通过ZH-ＪＤ方位角计算，公式为Ｇ±Ls/2/R×180/π,(曲线左转-，曲线右转+)。 三、缓和曲线 文件名：HHQX (COMP) “K0＝”？K: “X0=”？O：“Y0=”？“ALF=”？W：“LS=”？M：“R=”？R： Lbl 0：“K=”？K：“N=”？N：(-1)^N ×(K-L)^2÷M÷R÷6×180÷π→I：(K-L)-(K-L)^5÷90÷(RM)^2→D：O+Dcos(W+I) →X：P+Dsin(W+I) →Y：“XZ=”：X ◢“YZ=”：Y◢“LL=”？B：“XJJ=”？Q：X-Bcos(W+3I+Q) →E：Y-Bsin(W+3I+Q) →F：“XL=”：E◢“YL=”：F◢“R L=”？C：X+Ccos(W+3I+Q) →I：Y+Csin(W+3I+Q) →J：“XR=”：I◢“YR=”：J◢Goto 0 注：1、坐标计算方法是根据偏角法原理； 2、缓和曲线（ZH~HY或YH~HZ）以ZH（或HZ）为起始点； 3、平曲线左转（ZH~HY段N=1，YH~HZ段N=2），曲线右转（ZH~HY段N=2，YH~HZ段N=1）。 2008/12/8 网上搜索来的 一、程序功能 本程序由6个主程序、5个次子程序及5个参数子程序组成。主要用于公路测量中坐标正反算，设计任意点高程及横坡计算,桥涵放样，路基开挖口及填方坡脚线放样。程序坐标计算适应于任何线型. 二、源程序 1.主程序1：一般放样反算程序(①正算坐标、放样点至置仪点方位角及距离；②反算桩号及距中距离) 程序名:1ZD-XY Lb1 0:Norm 2 F=1：(正反算判别，F=1正算，F=2反算,也可以改F前加？，改F为变量) Z[1]=90（与路线右边夹角） Prog＂THB＂：F=1=>Goto 1:F=2=>Goto 2 Lb1 1:Ｆix 3:＂Ｘ=＂：Locate 6,4,Ｘ◢ ＂Ｙ＝＂：Locate 6,4,Ｙ◢ Ｐrog＂3JS”:Goto 0: Lb1 2:Fix 3:＂ＫＭ＝＂：Locate 6,4,Ｚ◢ ＂Ｄ＝＂：Locate 6,4,Ｄ◢ Ｇoto 0 ２．主程序２：高程序横坡程序(设计任意点高程及横坡) 程序名：2GC LbI 0:Norm 2 “KM”?Z:?D: Prog”H”:Fix 3:”H=”:Locate 6,4,H◢ “I=”: Locate 6,4,I◢ Goto 0 3.主程序3：极坐放样计算程序(计算放样点至置仪点方位角及距离) 程序名：3JS X：Y： 1268．123→K(置仪点X坐标) 2243．545→L（置仪点Y坐标，都是手工输入,也可以建导线点数据库子程序,个人认为太麻烦） Y-L→E：X-K→F：Pol(F,E):IF J<0:Then J+360→J:Int(J)+0.01Int(60Frac(J))+0.006Frac(60Frac(J)) →J:(不习惯小数点后四位为角度显示的，也可以用命令J◢DMS◢来直接显示) Fix 4:” FWJ=”: Locate 6,4,J◢(不习惯小数点后四位为角度显示的，也可以用命令J◢DMS◢来直接显示) Fix 3:”S=”:Locate 6,4,I◢ 4．主程序4：涵洞放样程序（由涵中心桩号计算出各涵角坐标、在主程序3中输入置仪点坐标后计算放样点至置仪点方位角及距离) 程序名：4JH-XY LbI 0:Norm 2 90→Z[1](涵洞中心桩与右边夹角，手工输入，也可以修改成前面加？后变为变量) 1→F:Prog”THB”:?L: Z[2]-Z[1] →E:X+Lcos(E) →X:Y+Lsin(E) →Y:Fix 3: ＂Ｘ=＂：Locate 6,4,Ｘ◢ ＂Ｙ＝＂：Locate 6,4,Ｙ◢ Ｐrog＂3JS”:Goto 0: 5．主程序5:路基开挖边线及填方坡脚线放样程序（输入大概桩号及测量坐标、地面标高计算出偏移距离、桩号、距中距离、填挖高度） 程序名：5FBX LbI 0:Norm 2: 18→DimZ:2→F:90→Z[1]:Prog “THB”:Z:D:”M0”?M:M→Z[4]:D→Z[3]:Prog”6GD”:L→Z[6]:If D<0:Then 0.75-L→D:Goto H:Else L-0.75→D:Goto H:IfEnd LbI H:Prog”H”:H-0.03-Z[4] →Z[5]:Z[6] →L:If Z[5]<0:Then –z[5] →G:Goto W:Else Z[5] →G:Goto T: LbI W:Prog “W0”:Z[10]+Z[11] →A: If G>A:Then Goto 1:Else If G>Z[10]:Then Goto 2:Else Goto 3:IfEnd: LbI 1:L+Z[12]+Z[13]+Z[14]+(G-A)×Z[9]+Z[11]×Z[8]+Z[10]×Z[7]:Goto Z: LbI 2:L+Z[12]+Z[13]+(G-Z[10])×Z[8]+Z[10]×Z[7]:Goto Z: LbI 3:L+Z[12]+G×Z[7]:Goto z: LbI T:L+0.5→N:If G>Z[17]:Then (N+Z[18]+(G-Z[17])×Z[16]+Z[17]×Z[15])→S:Goto Z:Else (N+G×Z[15])→S:Goto z: LbI Z:Z[3]→D:Fix 2:Abs(D)-S→T:”L0=”:L Locate 6,4,T◢ ＂ＫＭ＝＂：Locate 6,4,Ｚ◢ ＂Ｄ＝＂：Locate 6,4,Ｄ◢ “TW=”: Locate 6,4,Z[5]◢ 曲线任意里程中边桩坐标正反算(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 5800计算器测量计算程序 上上月做这个东西的时候没仔细检查，有好几处输错了的地方，今天把它修改过来。 简要介绍： 1. 新版程序把线元法和交点法已经集成在一个模块中了，用户只需修改JD 程序和ZA程序中的数据部分即可,其余不需作任何的改动。 2. 因为每条路高程计算不尽相同，且比较复杂，现在可利用PC机EXCEL计算好打印成表格带到工地上使用，所以本版程序未对线路高程序进行专门的编程计算，而是利用统计计算模式中来输入桩号（第一列X）及左、右高程（第二、三列Y，Freq），这种输入数据的方式最为直观，易发现错误，也易修改，输入完毕后运行S程序对数据按桩号进行排序，在程序中通过调用GG程序来进行内插计算，SG=-1得左标高，SG=1得右标高（若SG输入0,则可进行一般的线性内插计算）。 3. 在JD程序和XY程序中，先将一个计算单元的数据置入矩阵F中（1行8列或1行9列），这样程序可读性极好。 4．相比原CASIO4850程序操作习惯，作了一点小小的改动，测站坐标存在Z[10]，N中，X坐标原存在M中容易被误操作修改，而设计标高存在M中，这样易于修改，因为CASIO5800没有IN，OUT功能，很不方便。 4. 程序利用Z[2]变量值来判断是采用交点法还是线元法模型计算，Z[2]＝0为线元法，否则为交点法。 一、PQX程序：计算中边桩坐标及近似的桩号反算，在运行模式直接调用。 ①Z[10]→S:”XO”?S:S→Z[10]:”YO”?N:Prog “AU” ②Lbl 2:?L:Prog “Z”:Prog “E”:1n→O:90→S ③Lbl 4:”JJ”?S:”YC”?O:SO=0 =>Goto 2…原来lbl 后没有标号４的。 ④O=-1 =>Goto 6 ⑤“X,Y”:R+OCos(Z+S)→X▲U+OSin(Z+S)→Y▲Prog “D”:Goto 4 ⑥Lbl 6:Z[7]→X:Z[8]→Y:”XF”?X:”YF”?Y:XY=0 =>Goto 4 ⑦X→Z[7]:Y→Z[8]:Pol(X-R,Y-U+1p):Z+S-J→J:”YC,DL,L”:ICos(J)→O▲ISin( J)→I▲L+I▲Goto 6 二、P程序：在程序中提供一个自由运算的模式。 ①Lbl 1:”TMP”?I:If I≠0:Then “RST”:I▲Goto 1:IfEnd 二、LYC程序：进行桩号反算及边坡放样，在运行模式直接调用。 ①Prog “AU” ②Lbl 1:Z[7]→X: Z[8]→Y: Z[6]→S: ”XF”?X :X→Z[7]:”YF”?Y:Y→Z[8]: ”ZF”?S: S→Z[6] ③Lbl 2:Prog “Z”:Y=U =>Y+1p→Y ④Pol(X-R,Y-U):J-Z→J:Isin(J)→O:Icos(J)→I ⑤If Abs(I)≤0.1:Then Prog “E”:”L,YC”:L+I→L▲O▲Goto 3:IfEnd ⑥If Z[9]≠0:Then Pol(Z[9]-SO,I):πJZ[9]÷180→I:IfEnd ⑦”DL”:I▲L+I→L:Goto 2 ⑧Lbl 3: Z[6]→S:If S=0:Then Goto 1:IfEnd ⑧M→Z ⑨Lbl 4:”SG”?Z:Z→M:If Abs(Z)=1: Then Prog “GG”:Y→Z:If X=1:Then 公路边坡（开挖点或坡脚点）精准放样通用程序 抛弃渐进法、放坡脚线开挖点，一次准确定位。开挖点或坡脚点无论填方挖方通用，再也无需逐级边坡主子程搞边坡放样程序，根据数学模型推理计算编写的实用的5800程序，用此法可以解决边坡放样，卡西欧5800代码如下： 1.程序名：BP LbI S↙ Cls↙ “K×+×××”?V:“L（-1） Or R（+1）”？Q:“Z0(SJ)”？E: “H0(SJ)”？F:“Z1(W)”?A:“H1(W)”？B: “Z2(N)”？C↙ If C≥A:Then Cls:Locate 6 ,2 ,“ERR!”▲ Stop:Goto S:Else↙ Cls: “H2(N)”？U: “BP(1:m)”?N↙ If U=B:Then 0→M:Else Abs（（A-C）÷（B-U））→M:IfEnd↙Abs（U-F）→G:C-E→S:(MG-S)N÷（M-N）→L↙ Q=-1=>-(E+L) →D↙ Q=1=>E+L→D↙ Cls:“Z=”:Locate 4，1，D▲ Cls:“HC”？H↙ L÷Abs（H-F）→K↙ Cls:“K×××=”: “1:M=”:“1:m=” :“M- m=”:Locate 5，1， V:Locate 5，2，K:Locate 5，3，N:Locate 5，4，K-N▲ Goto S↙ IfEnd↙ 2.程序设计原理： 公路挖方段的开口线或填方段的坡脚线是沿着原地面随高程变 化的一条曲线，放样坡口开挖点或坡脚的坡脚点实际上是放样不同里程桩号断面地面线上的点。通常的放样法是：先计算断面上坡口（角）点附近一点的坐标进行放样并测出该点的高程，再根据设计边坡坡度值用逐渐趋近的方法放出开挖点或坡脚点，弊端是当放出的点高程满足边坡坡度设计要求时候却又很难 保证该点是否在该断面上了的，并且外业作业繁琐，反复试放，点位精度也低。 本程序是根据先放出断面上坡口（角）点附近的两点（不打桩）直接测量两点的实测高程，由程序内部计算出该坡口点所在的地面线的自然坡度，并与该横断面的设计边坡值联立方程，从而直接推导出计算真正坡口（角）点距该断面中桩的水平距离值，如此，代入一般的线路中边桩坐标计算程序中从而相当精确的计算出坡口（角）点的坐标利用全站仪直接一次放样的思路原理。同时程序加入了实际所在点位所形成的实际边坡坡比值与设计坡比值比较校核的功能，若与设计边坡吻合则就可一次直接打桩定位，若相差不大吻合，则是断面上坡口（角）点附近的两点选择的不合适，将两点重新缩进重新找 边界元与有限元 边界元法boundary element method 定义：将力学中的微分方程的定解问题化为边界积分方程的定解问题，再通过边界的离散化与待定函数的分片插值求解的数值方法。 所属学科：水利科技(一级学科) ；工程力学、工程结构、建筑材料(二级学科) ；工程力学（水利）(三级学科) 边界元法（boundary element method）是一种继有限元法之后发展起来的一种新数值方法，与有限元法在连续体域内划分单元的基本思想不同，边界元法是只在定义域的边界上划分单元,用满足控制方程的函数去逼近边界条件。所以边界元法与有限元相比，具有单元个数少,数据准备简单等优点.但用边界元法解非线性问题时，遇到同非线性项相对应的区域积分，这种积分在奇异点附近有强烈的奇异性，使求解遇到困难。 简介 边界元法是在有限元法之后发展起来的一种较精确有效的工程数值分析方法。又称边界积分方程-边界元法。它以定义在边界上的边界积分方程为控制方程，通过对边界分元插值离散，化为代数方程组求解。它与基于偏微分方程的区域解法相比，由于降低了问题的维数，而显著降低了自由度数，边界的离散也比区域的离散方便得多，可用较简单的单元准确地模拟边界形状，最终得到阶数较低的线性代数方程组。又由于它利用微分算子的解析的基本解作为边界积分方程的核函数，而具有解析与数值相结合的特点，通常具有较高的精度。特别是对于边界变量变化梯度较大的问题，如应力集中问题，或边界变量出现奇异性的裂纹问题，边界元法被公认为比有限元法更加精确高效。由于边界元法所利用的微分算子基本解能自动满足无限远处的条件，因而边界元法特别便于处理无限域以及半无限域问题。边界元法的主要缺点是它的应用范围以存在相应微分算子的基本解为前提，对于非均匀介质等问题难以应用，故其适用范围远不如有限元法广泛， FX5800计算器公路测量常用程序集 一、程序功能 本程序由6个主程序、5个次子程序及5个参数子程序组成。主要用于公路测量中坐标正反算，设计任意点高程及横坡计算,桥涵放样，路基开挖口及填方坡脚线放样。程序坐标计算适应于任何线型. 二、源程序 1.主程序1：一般放样反算程序(①正算坐标、放样点至置仪点方位角及距离；②反算桩号及距中距离) 程序名:1ZD-XY Lb1 0:Norm 2 F=1：(正反算判别，F=1正算，F=2反算,也可以改F前加？，改F为变量) Z[1]=90（与路线右边夹角） Prog＂THB＂：F=1=>Goto 1:F=2=>Goto 2 Lb1 1:Ｆix 3:＂Ｘ=＂：Locate 6,4,Ｘ◢ ＂Ｙ＝＂：Locate 6,4,Ｙ◢ Ｐrog＂3JS”:Goto 0: Lb1 2:Fix 3:＂ＫＭ＝＂：Locate 6,4,Ｚ◢ ＂Ｄ＝＂：Locate 6,4,Ｄ◢ Ｇoto 0 ２．主程序２：高程序横坡程序(设计任意点高程及横坡) 程序名：2GC LbI 0:Norm 2 “KM”?Z:?D: Prog”H”:Fix 3:”H=”:Locate 6,4,H◢ “I=”: Locate 6,4,I◢ Goto 0 3.主程序3：极坐放样计算程序(计算放样点至置仪点方位角及距离) 程序名：3JS X：Y： 1268．123→K(置仪点X坐标) 2243．545→L（置仪点Y坐标，都是手工输入,也可以建导线点数据库子程序,个人认为太麻烦） Y-L→E：X-K→F：Pol(F,E):IF J<0:Then J+360→J:Int(J)+0.01Int(60Frac(J))+0.006Frac(60Frac(J)) →J:(不习惯小数点后四位为角度显示的，也可以用命令J◢DMS◢来直接显示) Fix 4:” FWJ=”: Locate 6,4,J◢(不习惯小数点后四位为角度显示的，也可以用命令J◢DMS◢来直接显示) Fix 3:”S=”:Locate 6,4,I◢ 4．主程序4：涵洞放样程序（由涵中心桩号计算出各涵角坐标、在主程序3中输入置仪点坐标后计算放样点至置仪点方位角及距离) 程序名：4JH-XY LbI 0:Norm 2 90→Z[1](涵洞中心桩与右边夹角，手工输入，也可以修改成前面加？后变为变量) 1→F:Prog”THB”:?L: Z[2]-Z[1] →E:X+Lcos(E) →X:Y+Lsin(E) →Y:Fix 3: ＂Ｘ=＂：Locate 6,4,Ｘ◢5800边坡放样程序
线元法简介
5800计算器程序下载
FX5800道路路线测量程序
5800坐标程序
卡西欧fx5800计算器缓和曲线编程
公路测量卡西欧5800万能程序
5800简单全线坐标计算程序
Fx5800计算器公路测量程序设计
匝道线元法
5800平曲线放样程序(分段)
5800p全功能程序
fx-5800p全线坐标正反算带高程计算程序(线元法)
CASIO 5800计算器测量计算程序
公路边坡(开挖点或坡脚点)精准放样通用程序
边界元与有限元
5800卡西欧计算器程序大全