casio-5800p程序集锦（卡西欧5800编程计算机）

有关路基和隧道的CASIO5800P万能曲线计

算程序

曲线任意里程中边桩坐标正反算

(CASIO仅-5800P计算器)程序

一、程序功能

本程序由一个主程序亿HUCHENXU)和儿个子程序——正算子程序(SUB1)、反

算子程序(SUB2)、数据库(SUB3)、算方位角程序(FA)、使用方位角算坐标(NE)、

边坡放样程序(BIANPO)、隧道圆心放样程序(SDYX)构成，可以根据直线、圆

曲线、缓和曲线(完整或非完整型)的线元要素(起点坐标、起点里程、起点切线

方位角、线元长度、起点曲率半径、止点曲率半径)及里程边距或坐标，对该曲线

段范围内任意里程中边桩坐标进行正反算。

二、源程序

增加变量16-DimZ

1.主程序(ZHUCHENXU)

"1.SZ—XY,2.XY—SZ,3.SDYX,4.B1ANPO"?N：

Lbl1："ZHUANGHAO=''?S：Prog"SUB3"：N>2=>Goto2：Abs(S-O)一W：”

BIANZHU="?—Z:Prog"SUB1"："XS='':XTX/"YS='':Y—Y/"FS='':(F-90)-»F:

FADMS/Prog"FA"：Goto1:

Lbl2："CX="?X:X—I："CY="?Y:Y—J："DMG="?一Z[7]：Prog"SUB2"："S=

":(0+W)-S/"Z="Z—Z/IfN=3:ThenProg"SDYX":IfEnd：IfN=4:ThenProg"

BIANPO”:IfEnd:Goto2

2.正算子程序(SUB1)

0.1739274226->Z[1]：0.3260725774—B：0.0694318442-»K：0.3300094782-L：

(1-L)一F：(1-K).M：U+W(Z[1]cos(G+QEKW(C+KWD))+Bcos(G+QELW(C+

LWD))+Bcos(G+QEFW(C+FWD))+Z[1]cos(G+QEMW(C+MWD)))—X：V+W(Z

[1]sin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))

+Bsin(G+QEFW(C+FWD))+Z[1]sin(G+QEMW(C+MWD)))—Y：G+QEW(C+WD)

+90-F：X+ZcosF-X：Y+ZsinF—Y

3.反算子程序(SUB2)

G-90-Z[2]：Abs((Y-V)cosZ[2]-(X-U)sinZ[2])->W：0—Z：Lbl0：Prog“SUB

1”：Z[2]+QEW(C+WD)-L：(J-Y)cosL-(LX)sinL^Z：

IfAbs(Z)<lE-6：IhenGotol：ElseW+Z—W：Goto0:IfEnd

Lbl1：0一Z：ProgHSUBln：(J-Y)+sinF—Z

4.数据库(SUB3)

IfS<226100:ThenCls:Stop：Iffind

IfS<226255.833:Then226100—0：99037.97361U：4282.3590一V：201.3266

304一G：1(/5一p：I()45一R：253.543—H：0-Q:GotoIJffind

IfS>236005.651:ThenClsiStop:IfEnd

Lbl1:1+P-C：(P-R)+(2HPR)1D：18(RTT->E：

5.算方位角程序(FA)

LblO：“XC="？A:"YC="？T:X：Y

0->J：POI(X-A,Y-T)-I："BC="：I/IfJV0Then360+J—J：ElseJ->J:Iffi

nd

“FW="：JADMS/GotoO

6.使用方位角算坐标(NE)

LblO："X="?X:“Y=”?Y:“L="?L:"F°''?F：LxCos(F)+X—N:"N="N/LxSi

n(F)+Y—E:"E="E/

GotoO

7.边坡放样程序(BIANPO)

“QPDJL="?-Z[3]:“QPDGD="?->Z[4]:"PG="?->Z[5]：2->Z[6]："Pl="?

-Z[8]；"P2=”?-Z[9]：“P3=”？-Z-101(没有的话就删)：Z[7]-Z[旬―

Z[ll]：IfZ[ll]<0ThenAbsZ[ll]-"T="：Z[ll]/Else"W="：Z[ll]

/IfEnd

Z[ll]-Z[5]—Z[12]：IfZ[12]<1ThenZ[8]xZ[ll]->Z[13]：ElseZ[8]xZ[5]+Z[9]

xZ[5]x(Z[12]-l)—Z[13]：IfEnd:AbsZ-(Z[13]+Z[6]xlntZ[12]+Z[3])—Z[14]:

“BP+C-J=":Z[14]/

说明：

Z[3]=QPDJL.......................起坡点离中桩距离

Z[4]="QPDGD”..................起坡点设计标高

Z[5]="PG”.......................坡高

Z[6]="PTK”.......................平台宽度

Z[7]="DMG”.......................实测地面高度

Z[8]="P1”.......................一级边坡坡度

Z[9]="P2”.......................二级边坡坡度

Z[10]="P3”.......................三级边坡坡度(没有的话就删)

Z[14]="BP+C-J="......................边坡超欠挖宽度(+表超，-表欠)

8.隧道圆心放样程序(SDYX)

IFS<229000:THEN243.997+0.0205x(S-226520)-(S-227318)A2/2/l6000+1.503—Z

[10]:“NGC="：Z[10]/ELSE286.385-0.023x(S-230500)-(S-230593.753)A2/2/16500

0+1.503—Z[10]:“YLJ="：Z[10]JlfEnd

Z[10]-Z[7]—Z[15]:"+SG,-JD“：Z[15](Z「15F2+(ABS(Z)-6.8)A2)“R

=”:Z[16]/

三、使用说明

1、规定

(1)以道路中线的前进方向(即里程增大的方向)区分左右；当线元往左偏时，

Q=-l；当线元往右偏时，Q=l；当线元为直线时，Q=0。

(2)当所求点位于中线时，Z=0；当位于中线左锄时，Z取负值；当位于中线中

线右侧时，Z取正值。

(3)当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。

(4)当线元为圆|11|线时，无论其起点、止点与什么线元相接，其曲率半径均等

于圆弧的半径。

(5)当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10

的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时•，

曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆Illi线的

半径。

(6)当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定

的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半

径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。

2、输入与显示说明

SZ一XY,2.XY—SZ,3.SDYX,4.BIANPO

N?选择计算方式，输入1表示进行由里程、边距计算坐标；输入2表示由

坐标反算里程和边距。输入3表示隧道圆心放样和计算实测半径；输入4表示进

行路基边坡放样。

3、盖板涵斜交操作说明：

（以左前点为例）

先把盖板涵斜交角度减90后=人

再盖板涵左长度*CosA=B,盖板涵左长度*SinA=C,盖板涵身宽度/C

osA=D

利用YQXJS程序计算此桩号正交时B宽度的坐标值，及法线方位角。

再转到NE程序，X,Y,F三项不要修改，只改L=C+D/2（计算后点D前

为一）

4、SUB3数据库说明：

S待求桩号，0=曲线起点桩号，U曲线起点X坐标，V曲线起点Y坐标，G曲线

起点切线方位角，

P曲线起点半径，R曲线终点半径，H=曲线长度，Q转向（直线0,左转-1,右

转1）

FX5800计算器公路测量常用程序集2.2版

一、程序功能

主要功能：采用线元法与交点法相结合计算多条线路坐标正反

算，可算任意复杂线型及立交匝道，包括c型，s型、卵型、回头曲

线等；极坐标放样，全线路基边坡开挖口及坡脚放样计算、路基任意

点设计高程、横坡、设计半幅宽度.隧道欠超挖放样计算。

新版本优化：1、调整一次显示结果；2、交点法中考虑了不对称缓和

曲线；3、修改原版本部分地方笔误.（红色为修改处）

二、源程序

1.总主程序（1、坐标计算放样程序2、坐标反算程序；3、高程

计算查阅程序；4、路基半幅标准宽度查阅程序；5、路基边坡及

开挖口放样程序；6、路基标准距离放样；7、桥梁锥坡计算放样

程序；8、极坐标计算程序；9、隧道超欠挖计算程序）运行后输

入1~9数子则选择1至9的程序，返回时，在桩号输入-1,返回

选择选择计算类型。输入-2,返回选择线路。坐标计算中输入-3,

则显示本段曲线要素。

程序名：OZCX

Lbl0:17—DimZ：Norm2:1-A："A:XY=1,ZD=2,GC=3,GD=4,

BP=5,FM=6,ZP=7,JS=8,SD=9"?A:

A=l=>Goto1:A=2=>Goto2:A=3=>Goto3:A=4=>Goto4:

A=5=>Goto5:A=6=>Goto6:A=7=>Goto7:A=8=>Goto8:

A=9=>Goto9：

LbIl:Prog"DX":LbIA:Prog"QX":90一B:"PJ1"?B:B-

C:"PJ2"?C:B­Z[l]:C-*Z[8]:LbI"KM"?Z:Z=

-l=>Goto0:Z=-2=>GotoA:Z=-3=>GotoX:Prog

KM":?D:Porg"TUB":0-L:"LO"?L:Z[2]+Z[1]-Z[8]-

E:X+Lcos(E)-X:Y+Lsin(E)-Y:Prog"XY":Prog"

JS":GotoB：

LbI2:2-F:90-Z[l]:Prog"QX":LbIC:"KM"?Z:

Z=-l=>Goto0:Z=-2=>Goto2:Z=-3=>GotoX:Prog"KM":"X

0"?X:"YO"?Y:Prog"TUB":Porg"ZD":GotoC:

LbI3:Prog"QX":0-B:"H-B"?B:B­Z[9]:LbID:"KM"?

Z:

Z=-l=>Goto0:Z=-2=>Goto3:?D:Porg"H":Fix3:"H=":

Lcoate6,4,H-Z[9]：

"1=":Locate6,4,I：GotoD：

LbI4:Prog"QX":LbIE:"KM"?Z:Z=-l=>Goto0:Z=-2=>Got

o4:?D:Prog"GD":Fix3:"SJGD=":Locate7,4,L：Goto

E：

LbI5:Prog"QX":0.5-B:"TH-GD"?B:B-*Z[17]:LbIF:2-F:

90-Z[l]:"KM"?Z:Z=-1=>Goto0:Z=-2=>Goto5:Prog"K

M""XO"?X:"YO"?Y:OfM:"MO”?M:M-Z[4]:Prog"3F

BX"：

Fix2:P­D：Abs(D)-S-O:"LO=":Locate6,4,0：

Prog"ZD":Z[4]-T:"TW=":Lcoate6,4,T/GotoF:

LbI6:Prog"DX":LbIG:Prog"QX":LbIH:lfF:90-Z[l]:

"KM"?Z:Z=-l=>Goto0:Z=-2=>GotoG:Prog"KM":?D:Pro

g"THB":Prog"XY":Prog"JS":Prog"H":0-M:"MO"?M:

Fix2:H-M^T:"TW=":Locate6,4,T/GotoH:

LbI7:Prog"DX":LbII:Prog"QX":0-*D:"LD:Z-,Y+"?D:

Abs(D)-R:"LR"?R:"ZO"?Z:D—Z[6]:R—Z[7]:Z—Z[8]:

LbIJ:O—Y:"LO:SZ+,DZ-"?Y:Y=-l=>Goto0:Y=-2=>GotoI:

Prog"ZP":GotoJ:

LbI8:Prog"DX":LbIK:?X:X=-l=>Goto0:X=-2=>

Goto8:?Y:Prog"JS":GotoK:

LbI9:Porg"DX":LbIL:Prog"QX”：LbIM:2一F:90-Z[l]:"

KM"?Z:Z=-l=>Goto0:Z=-2=>GotoL:Z=-3=>GotoX:Prog"SD

1":0.5-A:“CQHD"?A:"XO"?X:"YO"?Y:0-M="MO"?M:

M—Z[8]:A-*Z[15]:Prog"5SD"：GotoM:

LbIX:MatFAGoto0

2、次程序:路基开挖边线及填方坡脚线放样程序(输入填方放宽

值、大概桩号及测量坐标、地面标高计算出偏移距离、桩号、距

中距离、填挖高度)

程序名:3FBX

Prog“THB”:DfP:Prog”GD”:IfD<0:Then0.75-L->D:Goto

H:ElseL-0.75-D:GotoH:IfEnd

LbIH:Prog“H”:H-0.03-Z[4]-Z[4]:Prog“Wl”：FfZ[4]

<0:Then-Z[4]-^G:GotoW:ElseZ[4]-G:GotoT:

LbIW:Z[8]+Z[9]-A:IfG>A:ThenGotol:ElseIfG>Z[8]:T

henGoto2:ElseGoto3:IfEnd:

LbIl:L+Z[10]+Z[ll]+Z[12]+(G-A)XZ[7]+Z[9]XZ[6]+Z[8]X

Z[5]-*S:GotoZ:

LbI2:L+Z[10]+Z[ll]+(G-Z[8])XZ[6]+Z[8]XZ[5]->S:Goto

Z:

LbI3:L+Z[10]+GXZ[5]-S:GotoZ:

LbIT:L+Z[17]fT:IfG>Z[15]:Then(T+Z[16]+(G-Z[15])XZ

[14]+Z[15]XZ[13])^S:GotoZ:Else(T+GXZ[13])-S:Goto

Z:

LbIZ

3.次程序3.隧道超欠挖值计算放样程序(输入隧道线路，大概

桩号、输入衬砌厚度、测量三维坐标，计算准确桩号及位置、计

算欠超挖值)

程序名：5SD

Prog”KM”:Prog“THB”：D-*Z[13]:0-D:Prog"H”:Cos(Z[l

0])XZ[4]-E:H+Z[5]-Z[5]:H+Z[7]-Z[7]:E+Z[5]一

Z[10]:IfZ[8]>Z[10]:ThenGotoR:ElseGotoS:IfEnd:

LbIR:V((Z[13]-Z[9])2+(Z[8]-Z[5])2)-Z[4]-Z[15]-L:G

otoL:LbIS:Z[4]-Z[6]-S:V(S-(Z[7]-Z[5])2)XS4-A

bs(S)-*T:Abs(Z[13]-Z[9])-T-*T:V(T2+(Z[8]-Z[7])2)-Z[6]

-Z[15]-L:GotoL:LbIL:Fix2:"L0="：Locate6,4,L：

Prog”ZD”:“H0="：Locate6,4,Z[8]-H/

4：极坐放样计算程序(计算放样点至置仪点方位角及距离)

程序名：JS

X：Y：Z[ll]-K：Z[12]fL：Pol(X-K,Y-L):IFJ<0:ThenJ+3

60fJ:IFEnd:Int(J)+0.Ollnt(60Frac(J))+0.006Frac(60Frac

(J))一J:Fix4:"FWJ="：Locate6,4,J：

Fix3:"S=":Locate6,4,1/

5.路基宽度子程序6：路基标准半幅宽度计算程序(对于设计有

加宽渐变的有用，如路基宽度无变化，则把此程序直接输入半幅

宽度值至L)

程序名：GD

1-S:Prog“G1”：Z-C-E:(B-A)*E/S+A-L:

6.坐标计算次程序(THB)

程序名：THB

IfQ>l:ThenGotoJ:Elsel+P-C：(P-R)4-(2HPR)fS：1

80+n-E：IfF=l:ThenAbs(Z-O)-W：Prog〃A〃：Goto

2:ElseX—I：Y—J：Prog"B"：O+W-Z：D—D：G

oto2:IfEnd:

LbIJ:IfF=l:ThenProg〃Z“：Gotol:ElseProg"ZX"：

Goto2:IfEnd:

LbIl:I+DXcos(Z[2]+Z[l])-X:J+DXsin(Z[2]+Z[l])一Y:

Lbl2

7.线元法正算子程序(A)

程序名：A

0.1184634425-A：0.2393143352-B：0.2844444444-N：0.0

46910077-K：0.2307653449-L：0.5-M：U+W(Acos(G+QEKW

(C+KWS))+Bcos(G+QELW(C+LWS))+Ncos(G+QEMW(C+MWS))+Bco

s(G+QE(1-L)W(C+(l-L)WS))+Acos(G+QE(1-K)W(C+(1-K)WS)))

-X：

V+W(Asin(G+QEKW(C+KWS))+Bsin(G+QELW(C+LWS))+Nsin(G+QE

MW(C+MWS))+Bsin(G+QE(l-L)W(C+(l-L)WS))+Asin(G+QE(l-K)W

(C+(1-K)WS)))-Y：

G+QEW(C+WS)-Z[2]：X+Dcos(Z[2]+Z[l])一X：

Y+Dsin(Z[2]+Z[l])一丫

8.线元法反算子程序(B)

程序名：B

G-90-T：Abs((Y-V)cos(T)-(X-U)sin(T))-W：0-D：L

bl0：Prog〃A〃：T+QEW(C+WS)-L：(J-Y)cos(L)-(I-X)

sin(L)-D：IFAbs(D)<0.01:ThenGoto1：ElseW+D-W:G

oto0:IfEnd:"*-1

Lbll:0-D：Prog〃A〃：(J-Y)4-sin(Z[2]+90)-D

9.交点法正算子程序(Z)

程序名：Z

H24-R4-24-HA(4)4-2688+RA(3)fA:H4-2TA⑶4-240+R,-B:

((H-N2)H-244-R)4-Sin(Abs(P))-((HA(4)-NA(4))4-26884-R

A(3))+sin(Abs(P))-E:(R+A)tan(Abs(P)4-2)+B-E-T:P+Abs

(P)->W:0-M:H-C:

IfZW0-T:ThenZ-0-*S:G->Z[2]:Goto2:IfEnd:

IfZW0-T+H:ThenZ-0+T->S:Prog“HX”:G+WK—Z[2]:Goto

4:IfEnd:

IfZ^O-T+HRXAbs(P)4-180+H-4-2-N4-2:Then180(Z-0+T-0.

5H)+R+r|fS:A+R(l-Cos(S))^B：II4-2-HA(3)4-2404-R2+Rs

in(S)-*A:R->M:G+WS^Z[2]:Goto4:IfEnd:

O-T+HRXAbs(P)4-180+H4-2+N4-2-Z-*S:(R+N2^-R4-24-NA(4)

4-26884-RA(3))tan(Abs(P)+2)+N+2-NA(3)+2404-R2+E->T:N

-II:ProgWHXW:G+P-*S:S-WK-*Z[2]:U+(T-A)Cos(S)-WBSin(S)

-I:V+(T-A)Sin(S)+WBcos(S)-J:Goto3:

LbI4:U+(A-T)cos(G)-WBsin(G)f]:V+(A-T)Sin(G)+WBcos(G)

—J:Goto3:

LbI2:U+Scos(Z[2])I:V+Ssin(Z[2])-J:LbI3:C->H

10.交点法缓和段转化子程序(HX)

程序名：HX

S-SA(5)4-404-R24-H2+SA(9)+3456+-(4)4-HA(4)-A:SA(3)+

6+R+H-SA⑺4-3364-RA⑶9HA⑶+SA(11)4-422404-RA(5)4-

HA⑸-B:90s29n+R9H-K:RH+S-M

11.交点法反算子程序(zx)

程序名：zx

Z:0-*D:LbI0:Prog“Z”：Pol(X-I,Y-J):J-Z[2]-*J:Isin(J)

一S:Icos(J)-I:IfAbs(I)<0.l:ThenZ+I—Z:S-D:Goto2:E

IseGoto1:LbI1:IfM=0:ThenZ+I^Z:Goto0:ElesPol(M-

WS,I):(_JMn)/180-I:Z+I-Z:Goto0:IfEnd

LbI2

12.高程计算子程序(H)

程序名：H

Prog“SI”：C-TfF:Z-F—S:C+TfE:G-TIfQ:IfT=0:ThenQ+

SI^H:Goto0:ElseIfZ<F:ThenQ+SI^H:Goto0:ElseIfZ

WE:ThenQ+SI+S24-24-R^H:Goto0:

LbI0:H:IfD=0:ThenGotoI:ElseProg“I”：H+V-H:Goto

I:LbII

13.高程超高计算程序(I)

程序名：I

IfZ[3]=l:ThenProg“II”:Goto1:ElseProg"12”:Goto

1:IfEnd:

LbI1:IfW=l:ThenGotoZ:ElseGotoX:IfEnd

LbIZ:IfS=0:ThenAbs(D)XM^V:Goto2:ElseAbs(D)X((N-

M)X(Z-C)4-S+M)^V:Goto2:IfEnd:

LbIX:IfS=0:ThenAbs(D)XM-V:Goto2:ElseAbs(D)X(((3

((Z-C)4-S)-2((Z-C)(S)A⑶)X(N-M))+M)-V:Goto2:IfEn

d:

LbI2:Abs(D)fE:V+E-I:I(E-K)-V

14.导线点子程序(DX)

程序名：DX

Z[ll]—K:Z[12]一L:“XZ”？K:“YZ”？L:K—Z[U]:L—Z[12]

15.线路选择子程序(线路选择输0时。则输曲线参数(QX)

程序名：QX

1—A:"QX:X1=1,X2=2W?A:A-*Z[3]:

IfA=O:ThenGoto1:ElseGoto2:IfEnd:

LbI1:“JDO”？0:“XO”?U:"YO”?V:"FW"？G:“LS”?H:?P:?

R:"J=2,3;X=(0,1,-1)w?Q:IfQ=2:ThenH-*N:Goto2:ElseI

fQ=3:Then“LS”?N:LbI2

16.线路中线元段判别子程序(KM)

程序名：KM

IfZ[3]=0:ThenGoto0:ElseIfZ[3]=l:ThenProg“XI”：

Goto0:ElseProg"X2”:Goto0:IfEnd:LbI0

17.坐标计算参数转化子程序(F)

程序名：F

Mat-0：MatF[l,2]-U:MatF[l,3]^V:MatF[l,4]

—G:MatF[l,5]-H:MatF[l,6]-P:MatF[l,7]-R:MatF

Q=2=>H^N:Q=3=>MatF[1,9]-*N

18.桥梁锥坡放样子程序(ZP)

程序名:ZP

V(Z[6]2(1-Y7Z[7]2))-X：Z[8]—Z:Prog”GD”:IfD<0:T

henL+X-D:Else-L-X—D:IfEnd:Z[8]+Y-Z:D:lfF:90-Z

[1]:Prog”KM”:Prog"THB”:Prog”XY”:Prog”JS”

19.显示子程序(XY)

程序名：XY

Fix3:"X="：Locate6,4,X：

tlY="：Locate6,4,Y：

程序名：ZD

Fix3:"KM="：Locate6,4,Z：

"D=":Locate6,4,D：

20.数据子程序（附后示例）

①程序名：XI（线路1坐标计算要素程序）

IfZW26615.555:Then[[25900,11587.421,1847.983,10V0

9'23.1”，715.555,IXIO'、ixiO45,0]]->MatF:Goto0:If

End:

IfZW26915.555:Then[[26615.555,11448.97,2550.016,10?

09'23.1v,300,1X1045,1800,1]]-*MatF:GotoO:IfEnd:

LbI0:Prog"F”

①-2程序名：X2（线路2坐标计算要素程序）

IfZW122269.856:Then[[121656.952,2603848.019,50835

1.147,98°53,22.6〃，660,41°47'44.1”,1415.097,3,3

00]]—MatF:Goto0:IfEnd:

IfZW125115.9:Then[[123984.982,2601974.464,509885.4

41,140°41'06.8〃，700,-39°3T43.1",2800,3,400]]

~*MatF:Goto0:IfEnd:

LbI0:Prog"F”

说明：1、交点法：IfZW交点HZ点桩号：Then[[交点桩号，交

点X坐标，交点丫坐标，前交点至止交点方位角，交点前段缓和

曲线长度，交点总偏角，圆半径，3,交点后段缓和曲线长度]]

-MatF:Goto0:IfEnd:（如对称缓和曲线，圆半径后输，2后

不用输后段缓和长，总偏角分正负，左偏为负，右偏为正）

2、线元法：IfZW本段线元终点桩号:Then[[本段线元起点桩

号，线元起点X坐标，线元起点丫坐标,线元起点切线方位角，线

元长度，线元起点半径，线元终点半径，线元判别]]一MatF:Go

to0:IfEnd:

（线元判别及线元半径见后说明）

②-1程序名：S1（线路1高程竖曲线要素子程序）

IfZ[3]=l:ThenGoto1:ElseProg"S2”:Goto0:IfEnd:

LbIl:IfZW26157.8:Then-200000-R:157.8-「26000-C:3

7.IfG:0.01222—1:Goto0:IfEnd:

IfZW2742L915:Then300000-R:171.915一「27250-C:32.6

5-*G:-0.00356^1:Goto0:IfEnd:

IfZW27889.029:Then1000000-R:134.029一「27755-C:36.

64-*G:0.0079^1:Goto0:IfEnd:

LbI0

②-2程序名：S2（线路2高程竖曲线要素子程序）

IfZW183.007:Then1365-R:20.475-T:162.532-C:46.1一G:

0.003—I:Goto0:IfEnd:

IfZW321.264:Then-1200-R:38.732-T:282.532-C:50.175

-*G:0.033-*I:Goto0:IfEnd:

LbI0

程序字母说明：C-竖曲线交点桩号；G-交点桩号高程（未竖曲线

调整的）；T-竖曲线切线长（标段结尾或开头没有竖曲线的输0）;

R-竖曲线半径（分正负：凹曲线为正，凸曲线为负）；I—竖曲

线前纵坡

③-1程序名：II（线路1高程超高参数子程序）

1.5fK:2-W:IfD<0:ThenGotoL:ElseGotoR:IfEnd:

LbIL:IfZW26615.555:Then0-S：-0.02-M:Goto0:IfEnd:

IfZW26735.555:Then120-S:26615.555-C：-0.02-M:0.02

—N:Goto0:IfEnd:

IfZW26795.555:Then60-S:26735.555-C:0.02-M:0.04-N:

Goto0:IfEnd

LbIR:IfZW26735.555:Then0-S：-0.02-M:Goto0:IfEnd

IfZW26795.555:Then60-S:26735.555-C：-0.02-M：-0.04

—N:Goto0:IfEnd

IfZW27436.951:Then0-S：-0.04-M:Goto0:IfEnd

LbI0

③-2程序名：12（线路2高程超高参数子程序）

0.5-*K:1-^W:IfD<0:ThenGotoL:ElseGotoR:IfEnd:

LbIL:If100:Then0-S：-0.015一业Goto0:IfEnd:

IfZW110:Then10^S:100-*C：-0.015^M：-0.02-*N:Goto0:1

fEnd:IfZW343.665:Then0-S:-0.02-M:Goto0:IfEnd

LbIR:IfZW100:Then0-S：-0.015-M:Goto0:IfEnd

IfZW110:Then10—S:100-C：-0.015-M：-0.02-N:Goto0:1

fEndIfZW375.665:Then0-S：-0.02-M:Goto0:IfEnd

LbI0

字程序中字母表示说明：

K—中央分隔带半幅宽（中桩标高至中桩中线距离，无分隔带则

为0）;W—超高方式参数（W=l为一般直线方式超高，W=2为三

次抛物线方式超高）;S—超高渐变段距离（不是渐变段则输入0）；

C-超高渐变段起点桩号（不是渐变段，无需输入）；M—超高段

起点横坡，N-超高渐变段终点横坡（不是渐变段，无需输入）。

@-1程序名：G1（线路1路基标准半幅宽度参数子程序）

IfZ[3]=l:ThenGotol:ElseProg"G2”:Goto0:IfEnd

LbI1:IfD<0:ThenGotoL:ElseGotoR:IfEnd:

LbIL:IfZW28110.727:ThenZ-C:17.25fA:Goto0:IfEnd:

IfZW28200.727:Then28110.727-C:17.25-A:21.25一B:90

-*S:Goto0:IfEnd

LbIR:ZW27927.478:ThenZfC:17.25-A:Goto0:IfEnd

IfZW28172:Then27927.478-C:17.25-A:27.031-B:244.52

2^S:Goto0:IfEnd

LbI0

@-2程序名：G2（线路2路基标准半幅宽度参数子程序）

IfD<O:ThenGotoL:ElseGotoR:IfEnd:

LbIL:IfZ^69.998:ThenZ-C:28.25-A:Goto0:IfEnd:

IfZ^82.402:Then69.998-C:28.25-A:25.545-B:12.404-

S:Goto0:IfEnd

LbIR:Z^69.998:ThenZfC:17.45-A:Goto0:IfEnd

IfZ^80.367:Then69.998-C:17.45-A:16.787-B:10.369-

S:Goto0:IfEnd:LbI0

字程序中字母表示说明：

C—宽度渐变段起点桩号（不是渐变段C=Z）；A—宽度渐变段起

点宽度；B—宽度渐变段终点宽度（不是渐变段，无需输入）；S

-宽度渐变段距离（不是渐变段，无需输入）。

⑤-1程序名：SD1（线路1隧道参数子程序）

IfZ[3]=1:ThenGoto1:ElseProg"SD2":Goto0:IfEnd:

LbI1:5.1^Z[4]:1.8fz[5]:1.4-Z[6]:0.59-Z[7]:0-Z[9]:

109OP30”-Z[10]:

LbI0

⑤-2程序名：SD2（线路2隧道参数子程序）

5.43fz[4]：1.634fz[5]:7.93fz[6]:1.634fz[7]:0-Z[9]:

90-Z[10]

字程序中字母表示说明：

Z[4]一隧道上拱半径；Z[5]一隧道上拱圆心至设计面高度；Z[6]

一隧道中腰部半径；Z[7]一中腰圆心至设计面高度；Z[9]一设计

隧道中线与设计路基中线偏移值，左为负数，右为正数；Z[10]

一上部拱跨隧道半幅角度；

@-1程序名：W1（线路1路基填挖边坡参数子程序）

IfZ[3]=1:ThenGoto1:ElseProg"W2":Goto0:IfEnd:

LbI1:0.75fz[5]:1-Z[6]:1.25fz[7]:10-*Z[8]:10^Z[9]:2.

6fz[10]:2一Z[ll]:2-Z[12]:1.5->Z[13]:l.75-Z[14]:8一Z

[15]:2-Z[16]：LbI0

⑥-2程序名：W2（线路2路基填挖边坡参数子程序）

1-*Z[5]:1-Z[6]:1.25fz[7]:10-Z[8]:10-*Z[9]:2.6fz[10]:

1.5-Z[ll]:L5-Z[12]:L5-Z[13]:1.75-Z[14]:8-Z[15]:2

-Z[16]

字程序中字母表示说明：

Z[5]一挖方第一阶边坡边率；Z[6]一挖方第二阶边坡边率；Z[7]

一挖方第三阶边坡边率；Z[8]一挖方第一阶高度；Z[9]一挖方第

二阶高度；Z[10]一挖方路基碎落台及水沟宽度；Z[ll]一挖方第

二台阶平台宽度；Z[12]一挖方第三台阶平台宽度；

Z[13]一填方第一阶边坡边率；Z[14]一填方第二阶边坡边率；Z

[15]一填方第一阶高度；Z[16]-填方第二台阶平台宽度。（注：

本程序只做出挖方三个台阶，填方二个台阶，如需增加，先需在

程序TUB和4FBZ程序中增加变量。再如一标段有坡率及坡高不

一样时，可以照其它参数程序一样，用判别语句。)

三、使用说明

1、规定

(1)“KM”为所求点桩号，反算时为输入大概桩号

(2)当所求点位于中线时，D=0；当位于中线左锄时，

D取负值；当位于中线中线右侧时，D取正值。

(3)当线元为直线时，其起点、止点的曲率半径为无穷

大，以10的45次代替。(在用交点法时，P子母代表偏角)

(4)当线元为圆曲线时，无论其起点、止点与什么线元

相接，其曲率半径均等于圆弧的半径。(交点法为圆半径)

(5)当线元为完整缓和曲线时，起点与直线相接时，曲

率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半

径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，

以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半

径。

(6)当线元为非完整缓和曲线时，起点与直线相接时，曲率

半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线

的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆

曲线相接时，曲率半径等于圆曲线的半径。

(7)Q值是判别用采用线元法还是交点法。交点法值为2或

3,2为对称缓利曲线段交点法，3为不对称缓和段交点法。线元

法，根据线元段取值，直线段取0,缓和段和圆曲线段，根据左

偏取T,右偏取1.

（8）为了区别字母0及数子0,特把字母0加粗。

2、输入与显示说明

输入部分：

"A:XY=1,ZD=2,GC=3,GD=4,BP=5,FM=6,ZP=7,

JS=8,SD=9"?输入：1、坐标计算放样程序2、坐标反算程

序；3、高程计算查阅程序；4、路基半幅标准宽度查阅程序；5、

路基边坡及开挖口放样程序；6、路基标准距离放样；7、桥梁锥

坡计算放样程序；8、极坐标计算程序；9、隧道超欠挖计算程序

XZ?YZ?为置仪点坐标，第一次运行输一次，以下运行不

需再输入，下次再运行此程序调上次输值，如需改，则输入，不

需改，则按确认键。

“QX:X1=1,X2=2"？选择经路，坐标计算中输入。值，变为

输入线元段曲线要素。输1为第一条线路。输2为第二条线路，

本程序暂设计实例为2条线路。

KM?正算时所求点的里程（反算输入大概桩号）输入-1,

返回输入选择计算类型，输入-2,返回选择线路，输入-3,调出

本段线元参数，以矩阵显示，再运行从头开始。

D?正算时所求点距中线的边距（左侧取负，值右侧取正

值，在中线上取零）

xo?反算时所求点的X坐标（放样程序中实测X坐标）

Y0?反算时所求点的丫坐标（放样程序中实测丫坐标）

显示部分：

L?涵洞放样程序中涵距中心桩号前后距离,前为正,后为

负.

M0?边坡放样程序中原地面标高输入

H-B?计算面至路面高差（默认路面为0）

PJ1?PJ2?涵洞计算中第一个夹角同PJo第二夹角为边线

至路前进方向夹角。（同理适合在桥梁中斜交桥，平行布置中，桥

台桩基坐标计算）默认值均为90度.

TH-GH?边坡放样中填方超放宽值，默认值为0.5m.

"LD:Z-Y+"?桥梁锥坡向左右方向边长，左为负数，右为

正数。

"LR"?锥坡向路线前进或向后退方向边长。

"Z0"?锥坡顶边桩号

输出部分

X=***正算时，计算得出的所求点的X坐标

丫二***正算时，计算得出的所求点的丫坐标

FWJ=***正算时，计算得出的所求点的至置仪点方位角

S=***正算时，计算得出的所求点的至置仪点距离

KM=***反算时，计算得出的所求点的里程

D=***反算时，计算得出的所求点的边距

H=***所求点位置设计路面顶标高

1=***所求点位置设计路面横坡（D输0时为，本段

路基纵坡）

SJLG=***所求点位置路基设计半幅标准宽度

L0=***边桩放样程序中实测点至设计边坡点距离，正

向内移，负向外移。（隧道程序中为超欠挖值，正为超挖，负为

欠挖）

TW=***边桩放样程序中，实测放样边坡点至止桩号设

计顶填挖高度值

H0=***隧道放样程序中测量点至设计面高度

四、结束语

1、本套程序实例线路给二条，如标段像匝道有多条线时，同理

增加。

2、由于本人没做过隧道，隧道程序只是参照本项目另一个标的

图纸所做，如有不同，可以调整参数。原理应是一样的。

3、为了更好的开发计算器程序，本人特建一QQ群，欢迎各位测

友加入，QQ群号是：24817026c

4、所有对程序有建议或问题的，加本人QQ后可以直接留言，由

于工作原因，可能不能马上给予答复，敬请谅解！（加QQ为好友

时请注明：计算器程序）

主要功能：采用线元法与交点法相结合计算多条线路坐标正反算，可算任意复杂线型

及立交匝道，包括c型，s型、卵型、回头曲线等；极坐标放样，全线路基边坡开挖口

及坡脚放样计算、路基任意点设计高程、横坡、设计半幅宽度计算.桥梁锥坡放样计算、

隧道欠超挖放样计算。

2009年11月16日.根据网友调试成功;更改THB程序中一处漏错.

2009年11月18日完成2.1版，经群友调试后，整理上传。新版本建主程序合并原所有

计算类型，在主程序中可选择操作类型。隧道欠超挖增加变量衬砌厚度，因有设计衬砌

厚度根据石岩来定的。增加锥坡放样计算。I/B]

[B]2009年12月4日更新为2.2版,新版本:1、调整一次显示结果；2、交点法中考虑了

不对称缓和曲线；3、修改原版本部分地方笔误.（红色为修改处）

2009年12月5日，修改F程序中最后一句为：Q=3=>MatF[l,9]—N.解决了线元法和

对称缓和交点法只要输八个参数的问题.不对称缓和段输九个参数.

2009年12月8日，经网友调试成功，修改F子程序中后面一句:Q=2=>H—N:Q=3=>Mat

F[l,9]—N.

曲^任意里程中遏格坐檄正反算5800p

TYQXJS（主程序名）

Lbl0:4—>Dimz:Cls

“l.SZ—XY”

”2.XY一SZ”?一NCs

“XO“？U:“YO”?V:“SO"？O:“FO"？G:“LS”？H:,,RO,,?P:,,RN,,?R:?Q

l+P-C：（P-R）+（2HPR）-D：180「一E

N=2=>Goto2

Lbll:Cls:?S

S=0=>Goto0（注:常S输入0畤，程序畲由Lbl0始重新建行）

?Z:Abs（S・O）—W

PrognSUBl":Cls

“xs="：x/

“YS=”:Y/

“FS=”:F-90ADMS，

Goto1

LbI2:Cls:?X

X=0=>Goto0（注:常X输入0畤，程序由畲Lbl0始重新遵行）

?Y:X—I:Y—J

Prog”SUB2”:O+W—S:Cls

“S=":S/

“z="：z/

Goto2

注:缸色部可以不要，亦可以要,它的作用是揩S或X输入0畤,程序畲重新由始

行,（即是可以言十算下段^元）

SUB1（子程序名）

0.1739274226—AO3260725774—B

0.0694318442—K03300094782一L

1・L—F:1・K—M

G+QEKW(C+KWD)—Z[1]

G+QELW(C+LWD)—Z[2]

G+QEFW(C+FWD)—Z[3]

G+QEMW(C+MWD)-Z[4]

Axcos(Z[l]HX(剩虢可省略)

X+Bcos亿[2])-X

X+Bcos(Z[3])一X

X+Acos(Z[4])->X

u+wx—x

Asin(Z[l])-Y

Y+BsiMZ[2])一Y

Y+Bsin(Z[3])->Y

Y+Asin(Z[4])->Y

V+WY—Y

G+QEW(C+WD)+90-F

X+Zcos(F)一X

Y+Zsin(F)->Y

SUB2(子程序名)

G・90一T

(Y-V)cos(T)-(X-U)sin(T)^W

Abs(W)-W:0一Z

Lbl4:ProgMSUBr,

T+QEW(C+WD)-L

(J-Y)cos(L)-(I-X)sin(L)->Z

IfAbs(Z)<lE-6:Then0-Z:ProgMSUBl"

(J・Y)+sin(F)—Z

ElseW+Z—W:Goto4:IfEnd

注:IE-6即是10的负6次方，亦即是等於0.000001

输入畤按SHIFT—log—(-)—6

亦可直接酹IfAbs(Z)<lE-6离成IfAbs(Z)<0.000001

注意:本程序是根摞yhhf编瘾的4800P改嘉而成,以下是源程序及算例

大家有冏题亦可在以卜^址查看及向yshf赞冏

曲线任意里程中边桩坐标正反算(CASIOfx-4800P计算器)程序

一、程序功能

本程序由一个主程序(TYQXJS)和两个子程——正算子程序(SUB1)、反算子程序(

SUB2)序构成，可以根据曲线段——直线、圆曲线、缓和曲线(完整或非完整型)的线

元要素(起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲

率半径)及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。另

外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲

线任意里程中边桩坐标进行正反算。本程序也可以在CASIOfx4500P计算器及CASIO

仅-4850P计算器上运行。

特别申明：(1).适用于弧长小于2倍半径的各种线元坐标正反算，精度优

于1mm；

(2).在引用该核心计算部分时，请注明来源。

二、源程序

1.主程序(TYQXJS)

"l.SZ=>XY"："2.XY=>SZ"：N：U"X0"：V"Y0"：O"S0"：G"F0"：H"LS"：P"R0"：

R"

RN"：Q：C=l-P：D=(P-R)-(2HPR)：E=180^7t：N=l=>Goto1：/>Goto2A—」

Lbl1：{SZ}：SZ：W=Abs(S-O)：Prog"SUB1"：X"XS"=XJ

Y"YS"=Y/

F"FS"=F-90J

Goto1<一1

Lbl2：{XY}：XY：I=X：J=Y：Prog"SUB2"：S"S『O+W/

Z"Z"=Z/

Goto2

2.正算子程序(SUB1)

A=0.1739274226：B=0.3260725774：K=0.0694318442：L=0.3300094782：F=l-L：

M=l-K：X=U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QEFW

(C+FWD))+Acos(G+QEMW(C+MWD)))：Y=V+W(Asin(G+QEKW(C+KWD))+Bsin(G+

QELW(C+LWD))+Bsin(G+QEFW(C+FWD))+Asin(G+QEMW(C+MWD)))：F=G+QEW(C+

WD)+90：X=X+ZcosF：Y=Y+ZsinF

3.反算子程序(SUB2)

T=G-90：W=Abs((Y-V)cosT-(X-U)sinT)：Z=0：Lbl0：Prog"SUB1"：L=T+QEW(C+

WD)：Z=(J-Y)cosL-(I-X)sinL：AbsZ<lE-6=>Gotol：/>W=W+Z：Goto0A」

Lbl1：Z=0：Prog"SUB1"：Z=(J-Y)-sinF

三、使用说明

1、规定

(1)以道路中线的前进方向(即里程增大的方向)区分左右；当线元往左偏时，

Q=-l；当线元往右偏时，Q=l；当线元为直线时，Q=0o

(2)当所求点位于中线时，Z=0；当位于中线左恻时，Z取负值；当位于中线中线右

侧时，Z取正值。

(3)当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。

(4)当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆

弧的半径。

(5)当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的

45

次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半

径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。

（6）当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的

值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等

于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。

2、输入与显示说明

输入部分：

1.SZ=>XY

2.XY=>SZ

N?选择计算方式，输入1表示进行由里程、边距计算坐标；输入2表示由坐标反

算

里程和边距。

X0?线元起点的X坐标

Y0?线元起点的丫坐标

S0?线元起点里程

F0?线元起点切线方位角

LS?线元长度

R0?线元起点曲率半径

RN?线元止点曲率半径

Q?线元左右偏标志（左偏Q=-l,右偏Q=l,直线段Q=0）

S?正算时所求点的里程

Z?正算时所求点距中线的边距（左侧取负，值右侧取正值，在中线上取零）

X?反算时所求点的X坐标

Y?反算时所求点的丫坐标

显示部分：

XS=xxx正算时，计算得出的所求点的X坐标

YS=xxx正算时，计算得出的所求点的丫坐标

FS=xxx正算时，所求点对应的中线点的切线方位角

S=xxx反算时，计算得出的所求点的里程

Z=xxx反算时，计算得出的所求点的边距

四、算例

某匝道的由五段线元（直线+完整缓和曲线+圆曲线+非完整缓和曲线+直线）组

成，各段线元的要素（起点里程S0、起点坐标X0Y0、起点切线方位角F0、线元长度

LS、起点曲率半径R0、止点曲率半径RN、线元左右偏标志Q）如卜.：

soX0Y0F0LSRORNQ

500.00019942.83728343.5611251631.00269.2561E451E450

769.25619787.34028563.3781251631.0037.4921E45221.75-1

806.74819766.56628594.5741202554.07112.779221.75221.75-1

919.52719736.07228701.893911730.6380.285221.759579.228-1

999.81219744.03828781.659804050.00100.0001E451E450

1、正算

（注意：略去计算方式及线元要素输入，请自行根据所求点所在的线元输入线元

要素）

S=700Z=-5计算得XS=19831.41785YS=28509.72590

S=700Z=0计算得XS=19827.33592YS=28506.83837

S=700Z=5计算得XS=19823.25398YS=28503.95084

S=780Z=-5计算得XS=19785.25749YS=28575.02270

S=780Z=0计算得XS=19781.15561YS=28572.16358

S=780Z=5计算得XS=19777.05373YS=28569.30446

S=870Z=-5计算得XS=19747.53609YS=28654.13091

S=870Z=0计算得XS=19742.68648YS=28652.91379

S=870Z=5计算得XS=19737.83688YS=28651.69668

S=940Z=-5.123计算得XS=19741.59118YS=28722.05802

S=940Z=0计算得XS=19736.47687YS=28722.35642

S=940Z=3.009计算得XS=19733.47298YS=28722.53168

2、反算

X=19831.418Y=28509.726计算得S=699.9999974Z=-5.00018164

X=19827.336Y=28506.838计算得S=699.9996493Z=0.000145136

X=19823.25398Y=28503.95084计算得S=699.9999985Z=5.000003137

X=19785.25749Y=28575.02270计算得S=780.0000035Z=-5.000001663

X=19781.15561Y=28572.16358计算得S=780.0000025Z=-0.000002