千赢首页

Casio5800交点法与线元法(积木法)匝道坐标正反算放样程序

Casio5800交点法与线元法(积木法)匝道坐标正反算放样程序(XUFENG 2011.2.14)本人一直以来想找一个交点法与线元法相结合的坐标正反算程序,在网上找了很久很久,没能找到一个较为满意的,有...

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↙
Cls:"1.ZS  2.FS"? I: I=2=>Goto 3↙
LbI 1 : Cls: If Z[61]=1: Then"JD ZS KX+XXX"?K :Prog"4": Else "ZHADAO ZS KX+XXX"?K :IfEnd↙
LbI  2: Cls:90→B: Cls:"RJ Or 0  To K"?B:B=0 =>Goto 1:"Z"?T↙
Prog "XY-A"↙
X+Tcos(M+B)→X↙
Y+Tsin(M+B)→Y↙
360Frac((M+360)÷360→M↙
Pol(X-H,Y-Z : 360Frac((J+360)÷360→J↙
2→O: Prog "XY-B":Goto 2↙
LbI 3 : Cls: If Z[61]=1: Then"JD FS KN+"?K:"X"?C:"Y"?D:Prog"4":Else Cls: " ZHADAO FS":"X"?C:"Y"?D:IfEnd↙
LbI 4 :Prog "XY-A"↙
(D-Y)sin(M)+(C-X)cos(M)→H↙
If Abs(H)>X10-3 :Then K+H→K:Goto 4:IfEnd↙
(D-Y)÷cos(M)→T↙
3→O: Prog "XY-B":Goto 3↙
3子程序名: XY-A(功能:坐标计算程序)
5→N: G(Q-1-P-1)÷Abs(L-A)→F: Abs(K-A)÷N→R: 90R÷π→S:
W+(FNR+2GP-1)NS→M:1→E↙
U+R÷6×(Cos (W)+Cos (M) +4∑(Cos (W+((E+0.5)FR+2GP-1)×(E+0.5)S),E,0,(N-1))+2∑(Cos (W+((EFR+2GP-1)ES,E,1,(N-1)))→X ↙
V+R÷6×(sin (W)+sin (M) +4∑(sin (W+((E+0.5)FR+2GP-1)×(E+0.5)S),E,0,(N-1))+2∑(sin (W+((EFR+2GP-1)ES,E,1,(N-1)))→Y↙
4子程序名: XY-B(功能:显示正算或反算结果)
If O=2:Then↙
Cls :"K×××=":"Z=":"X=":"Y=": Locate 6,1, K : Locate 4, 2, T : Locate 4,3, X : Locate 4,4, Y◢
If T=0 :Then Cls :"QF(Z)=": Locate 8,1, M:M▼DMS◢
IfEnd↙
Cls :"K×××=":"S=": Locate 6,1, K : Locate 4, 2, I :
"F=":J:J▼DMS◢
IfEnd↙
If O=3:Then "X=":"Y=":"K×××=":"Z=": Locate 4,1,C: Locate 4, 2, D : Locate 6,3,K :Locate 4,4,T◢
IfEnd:Cls↙
5子程序名:4(功能:将交点参数转为线元计算参数)
LbI 1: IF Z[48]<0 :Then -1→Z[62] : Else:1→Z[62]:IfEnd
LbI 2: If K≥Z[57]:Then Z[57]→A  :Z[1]→L  :Z[23]→U  :Z[24]→V : Z[31]→W : 10^45→P   :10^45→Q : 0→G:IfEnd↙
LbI 3:If K≥Z[1]:Then Z[1]→A : Z[2]→L : Z[19]→U : Z[20]→V  :Z[29]→W : 10^45→P   :Z[46]→Q : Z[62]→G: IfEnd↙
LbI 4:If K≥Z[2]:Then Z[2]→A : Z[4]→L  :Z[25]→U : Z[26]→V  :Z[32]→W : Z[46]→P  : Z[46]→Q : Z[62]→G: IfEnd↙
LbI 5:f K≥Z[4]:Then Z[4]→A : Z[5]→L : Z[27]→U  :Z[28]→V : Z[33]→W : Z[46]→P  : 10^45→Q : Z[62]→G: IfEnd↙
LbI 6:If K≥Z[5]:Then Z[5]→A : Z[5]+1000→L  :Z[21]→U : Z[22]→V : Z[30]→W  :10^45→P :  10^45→Q : 0→G : IfEnd↙
6子程序名:JDYS(功能:输入交点要素、显示交点要素及主点坐标)
Cls : " BP"?H:H→Z[57]:"K(JD)"?K:K→Z[41] :"X(JD)"?X :X→Z[42]:"Y(JD)"?Y:Y →Z[43]:"LS1"?B:B→Z[44] :"LS2"?C:C →Z[45]: ?R:R →Z[46]:"(ZH)FWJ°"?M:M→Z[47] : "α(Z-,Y+)°"?O:O→Z[48] : Z[47]+Z[48]→Z[49]: Prog "1":Prog "2"↙
Cls :"T1=":"T2=":"L=":"LY=": Locate 4,1, Z[50] : Locate 4,2, Z[51]: Locate 4,3, Z[52] : Locate 4,4, Z[53]◢
Cls :"E=": Locate 7,1, Z[54]
Cls :"K(QD)=": "X=": "Y=": "FWJ="Locate 7,1,Z[57] :Locate 7,2, Z[23] :Locate 7,3, Z[24] :Locate 7,4, Z[31] ◢
Cls :"K(ZH)=": "X=": "Y=": "FWJ=":Locate 7,1,Z[1] : Locate 7,2, Z[19] :Locate 7,3, Z[20] :Locate 7,4, Z[29]◢
Cls : "K(HY)=": "X=": "Y=": "FWJ=": Locate 7,1,Z[2] : Locate 7,2, Z[25] :Locate 7,3, Z[26] :Locate 7,4, Z[32]◢
Cls :"K(QZ)=": Locate 7,1,Z[3]◢
Cls :"K(YH)=": "X=": "Y=": "FWJ=": Locate 7,1,Z[4] : Locate 7,2, Z[27] :Locate 7,3, Z[28] :Locate 7,4, Z[33]◢
Cls :"K(HZ)=": "X=": "Y=": "FWJ=": Locate 7,1,Z[5] : Locate 7,2, Z[21] :Locate 7,3, Z[22] :Locate 7,4, Z[30]◢
7子程序名: 1(功能:计算交点要素)
If Z[48]<0 :Then  -1→Z[55] : Else 1→Z[55] : IfEnd : Z[55]* Z[48]→Z[56] ↙
Z[44] 2 ÷24÷Z[46]- Z[44]^(4)÷2688÷Z[46] ^(3) →Z[6] ↙
Z[45] 2 ÷24÷Z[46]- Z[45]^(4)÷2688÷Z[46] ^(3) →Z[7] ↙
Z[44]÷2-Z[44]^(3)÷240÷Z[46]2 →Z[8] ↙
Z[45]÷2-Z[45]^(3)÷240÷Z[46]2 →Z[9] ↙
Z[8]+(( Z[46]+Z[7]-( Z[46]+Z[6])cos(Z[56]))÷sin(Z[56]))→Z[50]↙
Z[9]+(( Z[46]+Z[6]-( Z[46]+Z[7])cos(Z[56]))÷sin(Z[56]))→Z[51]↙
Z[46]* Z[56]π÷180+( Z[44]+ Z[45]) ÷2→Z[52]↙
Z[46]* Z[56]π÷180-( Z[44]+ Z[45]) ÷2→Z[53]↙
(Z[46]+(Z[6]+Z[7])÷2)÷cos(Z[56]÷2)- Z[46]→Z[54]↙
Z[41]-Z[50]→Z[1] ↙↙
Z[1]+Z[44]→Z[2] ↙↙
Z[2]+Z[53]÷2→Z[3]↙
Z[1]+Z[52]-Z[45]→Z[4]↙
Z[4]+Z[45]→Z[5]↙
8子程序名: 2(功能:计算主点坐标及切线方位角)
Z[42]-Z[50]cos(Z[47])→Z[19]: (直缓坐标)
Z[43]-Z[50]sin(Z[47])→Z[20]↙
Z[47]→Z : 360Frac((Z+360)÷360→Z[29] (方位角)
Z[42]+Z[51]cos(Z[49])→Z[21]: (缓直坐标)
Z[43]+Z[51]sin(Z[49])→Z[22]↙
Z[49]→Z: 360Frac((Z+360)÷360→Z[30] (方位角)
Z[1]-Z[57]→L↙     (H→Z[57]为前直线起点桩号)
Z[42]-( Z[50]+L)cos(Z[47])→Z[23]↙       (前直线起点坐标)
Z[43]-( Z[50]+L)sin(Z[47])→Z[24]↙
Z[47]→Z : 360Frac((Z+360)÷360→Z[31]↙   (方位角)
Z[44]→Z[12]:Z[44]→Z[13]:Prog"3"↙
Z[4]-Z[1]→L:90(2L-Z[44])÷Z[46]÷π→Z[11]↙
Z[46]sin(Z[11])+Z[8]→Z[14]:Z[46](1-cos(Z[11]))+Z[6]→Z[15]↙
Z[19]+Z[14]cos(Z[47])-Z[55]Z[15]sin(Z[47])]→Z[27]↙  (圆缓点坐标)
Z[20]+Z[14]sin(Z[47])+Z[55]Z[15]cos(Z[47])]→Z[28]↙
Z[47]+Z[55]Z[11]→Z: 360Frac((Z+360)÷360→Z[33]↙(方位角)
Z[2]-Z[1]→L:90(2L-Z[44])÷Z[46]÷π→Z[58]↙
Z[46]sin(Z[58])+Z[8]→Z[14]:Z[46](1-cos(Z[58]))+Z[6]→Z[15]↙
Z[19]+Z[14]cos(Z[47])-Z[55]Z[15]sin(Z[47])]→Z[25]↙  (缓圆点坐标)
Z[20]+Z[14]sin(Z[47])+Z[55]Z[15]cos(Z[47])]→Z[26]↙
Z[47]+Z[55]Z[58]→Z: 360Frac((Z+360)÷360→Z[32]↙(方位角)
9子程序名: 3(主点坐标计算辅助程序)
If Z[12]=0 :Then 0→Z[14]: 0→Z[15]:Else↙
Z[12]- Z[12]^(5)÷40÷(Z[46]*Z[13])2+ Z[12]^(9)÷3456÷(Z[46]*Z[13])^(4) →Z[14]↙
Z[12]^(3)÷6÷(Z[46]*Z[13])-Z[12]^(7)÷336÷(Z[46]*Z[13])^(3)+ Z[12]^(11) ÷42240÷(Z[46]*Z[13])^(5)→Z[15] ↙
IfEnd↙

程序说明:
1、进入程序:1.JD ZFS  2. ZHADAO ZFS? 选1为交点法正反算(以后操作均为交点法计算),选2为线元法正反算(以后操作均为线元法计算)
2、ZHONG SHU JS 2. JS?选1重输参数,选2直接进入交点法或线元法正反算(参数为已输过的参数)
3、参数输入:
一、交点法已知数据输入:
BP?上一交点ZH桩号
K(JD)?交点桩号
X(JD)?交点X坐标
Y(JD)?交点Y坐标
LS1  ?第一缓和曲线长度
LS2  ?第二缓和曲线长度
R  ? 圆曲线半径
(ZH)FWJ°?交点前(即前交点至本交点也即ZH点)的正切线方位角
α(Z-,Y+)?本交点处线路转角(左转为负,右转为正,度分秒输入)
交点法计算要素显示:
T1=第一切线长
T2=第二切线长
L=曲线总长
LY=圆曲线长
E=曲线外距
K(ZH)=直缓点桩号
K(HY)=缓圆点桩号
K(QZ)=曲中点桩号
K(YH)=圆缓点桩号
K(HZ)=缓直点桩号
二、线元法已知数据输入:
K0? KN? R0? RN?F0?X0? Y0?ZX? 分别为线元起点桩号 、终点桩号、起点半径、终点半径、起点切线方位角、起点X坐标、起点Y坐标、线元转向。
4、XC ?    YC? 输入置镜点即测站的X,Y坐标
5、1.ZS 2.FS选1正算 选2反算
程序线元判断原则:
(1) 以道路中线的前进方向(即里程增大的方向)区分左右;
(2) 当所求点位于中线时,Z=0;当位于中线左侧时,Z取负值;当位于中线右侧时,Z取正值。
(3) 当线元为直线时,其起点、终点的曲率半径为无穷大,以10的45次代替输入×10x45。
(4) 当线元为圆曲线时,无论其起点、终点与什么线元相接,其曲率半径均等于圆弧的半径。
(5) 当线元为完整缓和曲线时,起点与直线相接时,曲率半径为无穷大,以10的45次代替输入×10x45;与圆曲线相接时,曲率半径等于圆曲线的半径。终点与直线相接时,曲率半径为无穷大,以10的45次代替输入×10x45;与圆曲线相接时,曲率半径等于圆曲线的半径。
(6) 当线元为非完整缓和曲线时,起点与直线相接时,曲率半径等于设计规定的值;与圆曲线相接时,曲率半径等于圆曲线的半径。终点与直线相接时,曲率半径等于设计规定的值;与圆曲线相接时,曲率半径等于圆曲线的半径。
6、循环计算部份:
输入与显示简单说明
JD ZS KX+XXX? 或 ZHADAO ZS KX+XXX?输入所求的桩号(交点法或线元法 )
Z ?              输入所求点距中线的边距(在中线输零,左负右正)
RJ?              输入边桩时左右边桩连线与线路前进方向中桩切线的右交角(当输入数字0时进入下一个桩号计算输入)
X =、Y =            计算得出的所求点的左、中、右 的X Y坐标
QF(Z)= ×××            计算得出所求点的中桩切线方位角
F= ×××    计算得出置镜点到测点的方位角
S= ×××   计算得出置镜点到测点的水平距离
JD FS KN+?  交点法反算时需输入反算点所在(前直线,前缓和段、圆曲线、后缓和段、后直线)线元中的任一点桩号(定位线元用)
X=×××    反算输入所求点的X坐标
Y=×××     反算输入所求点的Y坐标
K=×××      计算得出求点所对应的里程
Z=×××       计算得出求点到所对应的里程的垂直距离 (负就是左边,正就是右边)

  • 发表于 2018-11-05 14:18
  • 阅读 ( 194 )
  • 分类:5800

你可能感兴趣的文章

相关问题

0 条评论

请先 登录 后评论
admin
admin

117 篇文章

作家榜 »

  1. admin 117 文章
  2. 胡亮 3 文章
  3. 梁雄 2 文章
  4. 肖肖 2 文章
  5. 12 1 文章
  6. boyving 1 文章
  7. 张峰 0 文章
  8. 陈建业 0 文章