4800公路测量放样计算程序

CASIO fx-4800P计算器程序

DA DI——N E (公路弯道坐标计算)

K“JD”：R：S“LS”：F“PJ——L－R＋”：A“ZH——N”：B“ZH——E”：C“JD——N”：D“JD——E”

M“ZH”=K－（R＋S2÷（24R））tan（Abs F÷2）－S÷2＋S^3÷（240R2）▲

N“HY”=M＋S▲

O“QZ”=M＋（πRAbs F÷180＋S）÷2▲

Z[1]“YH”=M＋πRAbs F÷180▲

Q“HZ”= Z[1]＋S▲

T=K－M▲

L=Q－M▲

E=（R＋S2÷（24R））÷cos（F÷2）－R▲

Lb1 0：{G}：G“ZHUANG HAO”：G＜O＝=＞Z=G－M≠=＞Z=Q－G：△

{U}：U“I——B OUT＋IN－”

Z＜S＋0.005＝=＞H=tan-1（24Z2S2R2－Z^6）÷（48Z^3R^3－6Z^4SR））：V=Z－Z^5÷（40S2R2）＋Usin H：W=Z^3÷（6SR）－Z ^7÷（336Z^3R^3）－Ucos H：≠=＞H=180（Z－S÷2）÷（πR）：V=Rsin H＋S÷2－S^3÷（240R2）＋Usin H：W=R－Rcos H＋ S2÷（24R）－Ucos H：△

F＜0＝=＞G＜O＝=＞X=V：Y=W：≠=＞X=T＋Tcos F－Vcos F－Wsin Abs F：Y=Tsin Abs F－Vsin Abs F＋Wcos F：△≠=＞G＜O＝=＞X=V：Y=－W：≠=＞X=T＋Tcos F－Vcos F－Wsin Abs F：Y=Vsin Abs F－Tsin Abs F－Wcos F：△△

C=A＝=＞D≥B＝=＞P=90：≠=＞P=270：△≠=＞C＜A＝=＞P= tan-1（（D－B）÷（C－A））＋180：≠=＞P= tan-1（（D－B）÷（C－A））：△△

I“COM——N”=A＋√（X2＋Y2）cos（P－tan-1（Y÷X））▲

J“COM——E”=B＋√（X2＋Y2）sin（P－tan-1（Y÷X））▲

G=G＋20

Goto 0

程序运行

输入：JD? ——交点桩号

R? ——圆曲线半径

LS? ——缓和曲线长度

PJ——L－R＋? ——偏角 右偏为正

ZH——N? ——ZH点N坐标

ZH——E? ——ZH点E坐标

JD——N? ——JD点N坐标

JD——E? ——JD点E坐标

输出：ZH= ——ZH点桩号

HY= ——HY点桩号

QZ= ——QZ点桩号

YH= ——YH点桩号

HZ= ——HZ点桩号

T= ——切线长

L= ——曲线长

E= ——外矢距

再输入：ZHUANG HAO? ——输入任一点桩号

I——B OUT＋IN－? ——距中线距离（外侧为正）

再输出：COM——N= ——计算出N坐标

COM——E= ——计算出E坐标

以此循环——再输入——再输出

程序中各字母所代表的含义：

A——ZH点N坐标

B——ZH点E坐标

C——JD点N坐标

D——JD点E坐标

E——外矢距

F——偏角（右偏为正）

G——弯道内任一点桩号

H——弯道内任一点切线角

I——计算出的N坐标

J——计算出的E坐标

K——交点桩号

L——曲线长

M——ZH点桩号

N——HY点桩号

O——QZ点桩号

P——ZH--JD方位角

Q——HZ点桩号

R——圆曲线半径

S——缓和曲线长度

T——切线长

U——距中线距离（外侧为正）

V——支距

W——支距

X——弯道内坐标

Y——弯道内坐标

Z——任一点曲线长度

Z[1]——YH点桩号

另个程序

公路逐桩坐标计算程序

(可以计算对称、不对称缓和曲线)

Lb1 0

Z=?

V=?

W=V+2：Fixm

｛K｝

Lb1 1

K>Z[W+5Z+4]=>W=W+1:Goto 1⊿ (判断桩号在哪个交点范围,就是该交点曲线起点至下一交点曲线起点)

S=K-Z[W+5Z+3] （计算该桩号与曲线起点的距离）

R=Z[W+2Z+2]：L=Z[W+3Z+2]:E=Z[W+4Z+2] （读取该交点曲线要素R、Ls1 、Ls2）

Pol（Z[W]-Z[W-1],Z[W+Z+2]-Z[W+Z+1]） （计算该交点与下一交点直线方位角）

JJ=J+360⊿

A=J

Pol（Z[W-1]-Z[W-2],Z[W+Z+1]-Z[W+Z]） （计算该交点与上一交点直线方位角）

JJ=J+360⊿

C=A-J:A=J （计算偏角）

W=V+2=>Goto2⊿ （如果桩号在起点与第一交点曲线起点之间,则转 Lb1 2 ）

I=Abs（tan（c÷2））

M=L÷2-L^3÷240R^2:N=E÷2-E^3÷240R^2

P=L^2÷6R-L^4÷336R^3-R（1-cos（90L÷πR））

Q=E^2÷6R-E^4÷336R^3-R（1-cos（90E÷πR））

D=（P-Q）I÷2 : F=（P+Q+2R）I÷2

M=F+M-D:Q=F+N+D

N=πRAbsC÷180+（L+E）÷2

X=Z[W-1]-McosA

Y=Z[W+Z+1]-MsinA

M=Z[W-1]+Qcos（A+C）

V=Z[W+Z+1]+Qsin（A+C）

Q=AbsC÷C

S≤L=>P=0:Goto3⊿ （如果桩号在第一缓和曲线内,则转 Lb1 3）

S≤N-E=>S=S-L:Goto4⊿（如果桩号在圆曲线内,则转 Lb1 4）

S≤N=>S=N-S

Q=-Q:A=A+C-180:X=M:Y=V:L=E：P=180:Goto3 ⊿ （如果桩号在第二缓和曲线内,则转 Lb1 3）

P=A+C:S=S-N:D=M+ScosP:F=V+SsinP

Goto6 （如果桩号在直线内,则转 Lb1 6）

Lb1 2

P=A+C

D=Z[W-1]+ScosP

F=Z[W+Z+1]+SsinP:Goto6

Lb1 3

I=S-S^5÷40R^2÷L^2+S^9÷3456R^4÷L^4

J=Q（S^3÷6RL-S^7÷336R^3÷L^3）

P=P+A+90QS^2÷πRL:Goto5

Lb1 4

M=90（2S+L）÷πR

I=RsinM+L÷2-L^3÷240R^2

J=Q（L^2÷24R+R（1-cosM））

P=A+QM

Lb1 5

D=X+IcosA-jsinA:F=Y+JcosA+IsinA

Lb1 6

D〃X=〃◢ （结果显示X坐标）

F〃Y=〃◢ （结果显示Y坐标）

P〃AT=〃◢ （结果显示该桩号方位角）

{BO}：B〃S〃O〃⊿〃 （输入边桩距离,交角）

P=P+O

L〃XB〃=D+BcosP◢ （结果显示边桩X坐标）

M〃YB〃=F+BsinP◢ （结果显示边桩Y坐标）

以上是坐标计算程序,括号内是程序计算的大致原理及说明,中间部分为直线、圆曲线、缓和曲线计算的各种公式,大家也知道,书上也有.

该程序是最基本的,如再加几条语句就可以处理断链.

该程序还需变量输入,把整条线路曲线要素一次输入到4800中,就可以计算线路任意中、边桩坐标.