施工放样 计算器应用程序

一、数学模型：

坐标系统分：以ZH、HZ为原点交点方向为X轴。

1) 测点里程处在第一缓和曲线时，测点坐标计算：

Xi=K i-（Ki）5/40R2·C2

Yi=（Ki）3/6RC-（Ki）7/336R3·C3

C=L0、K=L、K i=桩点里程-ZH里程

2) 测点里程处于曲线时、测点坐标及偏角采用：Ai=Ki/R·180/π＋90C/πR

Xi=R·SinAi＋（C/2－C3/240R3）

Yi=R（1－CosAi）＋C2/24R

Ki=桩点里程－HY里程

B0=90C/πr、M=（C/2—C3/240R2）

ρ=C2/24R

3) 测点里程处于第二缓和曲线时测点坐标采用：Xi=（T－Xt）Cos A－YtSinA＋T

Yi=（T－Xt）Sin A＋YtCosA

式中：

Xt=（Ki）－（Ki）5/40R2·C2

Yt=（Ki）3/6RC－（Ki）7/336R3·C3

B0=L0/2R·180/π=90C/πr

m=L0/2－L03/240R2

ρ=L02/24R

C=L0缓和曲线长

Ai= 圆弧所对的圆心角

A= 路线的转折角

4)曲线要素计算公式：

T=(R+P)tag(A/2)+M

L=A·R·π/180+C

E0=(R+P) ·Scs(A/2)-R

g=2T-L

P=C2/24R

M=C/2-C3/240R2

注：L=( R·∏/180)(a-2B0)+2L0

M=m0

5）里程计算：

HY里程=ZH里程+C

QZ里程=ZH里程+L/2

YH里程=ZH里程+L-C

HZ里程=ZH里程+L

6)测站点坐标计算

X ＝X ＋S ·CosI

其中X 、Y 初值为0，I=I+180-B

Y=Y+S ·SinI

7)欲测点（里程为K ）的拨角和距离计算

II=arcTg(YI/XI)-AI=arcTg((YI 中-Y)/(XI 中-X))-AI DI=

))Y YI ()X XI 22中－中－（+

COSIO FX-4500 缓和曲线计算程序

A:L1:DEG:P ”A0”:ProgN:I=P:L ”I ”:R:F:

A=L 2/24/R: M=L/2-L 3/240/R 2:B=90L/π/R L2:T ”T ”=(R+A)Tan(I/2)+M ▲D ”L ”=R π/180×(I-2B)+2L ▲ E ”E ”=(R+A)/Cos(I/2)-R ▲S ”Q ”=2T-D ▲H=D-L

L3:C ”X ”:O ”Y ”{P}:p”

L4:LbiO:{G}G ”L ”:G> D=>GOTOJ:△G>H=>GOTO1: △G>L=>GOTO2: △N=G:ProgI:GOTO3

L5:Lbi1:N=D-G:ProgI:V=X: W=Y: X=(T-V)ConI-W SinI+T: Y=(T-V)Sin+W CosI:GOTO3

L6:Lbi2:N=G-L:J=N/R ×180/π+B: X=R SinJ+M: Y=R(1-CosJ)+A: GOTO3

L7:Lbi3:X=0=>K=0:≠K=1/Tag(Y/X): △J=Z+F ×K: S=

(X2+Y2):V ”X ”=C+S CosJ ▲ W ”Y ”=O+S SinJ ▲ GOTO0 L8:Lbi:

(-2) ▲

I:L1:X=N-N 5/40/R 2/L 2: Y=N 3/6/R/L-N 7/336/R 3/L

3

N:L1:P: Q=FracP 102: P=IntP+IntQ/60+FracQ/36

A:坐标传递计算

L1:DEG:A ”XO ”:B ”Y0”:P ”<0”:ProgN:C=P

L2:LbiO:{P}:P ”<”:ProgN:J=P:{S}

L3:C=C+J-180:C<0=>C ”<”=C+360▲ ≠>C:C ≥360=>C ”<”=C-360▲≠>C ”<”=C ▲△△

L4:A ”X ”=A+Cos C ×S ▲

L5:B ”Y ”=B+SinC ×S ▲

L6:GOTO 0

α=θ+α1+180°

Xn=X+L ×Cos α

Yn=Y+L ×Sin α

B:直角坐标求极坐标

L1:DEG: A ”X0”: B ”Y0”

L2:Lbi:{C }:C”X1”:{D}:D ”Y1”

L3:E=C-A:F=D-B

L4:S ”S ”=Poi(E,F)▲W ＜0=＞H ”＜”=360+W ▲≠＞H ”＜”= W ▲△

L5:GOTO 0

C:坐标变换施工0坐标（旧变新）

L1:DEG: Lbi 0

L2:｛A ｝:A ”XA ”: ｛B ｝:B ”YA ”:｛P ｝:P ”＜A ”:Prog N:C=P: C=1(TanC)

L3:｛D ｝:D ”XB ”:｛E ｝:E ”YB ”:｛P ｝:P ” ＜B ”: ProgN: F=P: F=1/(TanF)

L4:G ”XP ”=(F*A+D*C-B+E)-(C+F) ▲

L5:H ”YP ”=(B*F+E*C+A-D)/(C+F) ▲GOTO 0

说明：计算公式为

Xp=(Xa*TanB+ Xb*TanA-Ya+Yb)/(TanA+ TanB)

Yp=(Ya* TanB+Yb*TanA+ Xa-Xb)/ (TanA+ TanB)

点位编号应为逆时针旋转方向编号。

E:多点支导线

L1:A ”X0”:B”Y0”:P” ＜0”: ProgN: C=P

L2:Lbi0:｛p ｝p ” ＜0”: ProgN: J=P:｛S ｝

L3:D=C=J-180:D ＜0=＞D ”＜”=D+360▲≠＞D:D ＞360=＞D ”＜”= D-360▲≠＞D ”＜”=D ▲△△

L4:E ”X ”=A+CosD*S ▲

L5:”Y ”=B+SinD*S ▲GOTO 0

F:极坐标放样，计算边长，拨角

L1: DEG:A ”X0”: B ”Y0” : C ”X1”: D ”Y1” :E=A-C:F=B-D:S ”S0”=Poi(E,F) ▲G=W

L2: Lbi0: ｛X ｝: ｛Y ｝:H=X-C:I=Y-D:S ”S ”=Poi(H,I) ▲W=W-G:W<0=＞J ”<”=360+W ▲≠＞H ”<”=W ▲△GOTO 0

G:联系三角形算

L1:Lbi0:p ”A ＜”: ProgN: C=P

L2:A ”SA ”:B ”SB ”:E ”SC ” L3:D=((A-CosC) *E *B *2-E2(B2 )/3:A ”SA ”=A-D ▲B ”SB ”=B-D ▲B ”SC ”=E+D ▲

L4:C ”＜A ”=C ▲G=A/SinC:H ”＜B ”=1/Sin(B/G) ▲F ” ＜Y ”=180-1/Sin(E/G) ▲GOTO 0

计算：开计算机，用点度形式输入a 角值。再依次输入a 、b 、c 各丈量连长。A 角值为仪器观测值。输出a 、b 、C 边长平差数值后，再输出a 、β、γ的平差后角值。（单位：度）