5800計算器程序辛普森公式程序(通用)解析_第1頁
5800計算器程序辛普森公式程序(通用)解析_第2頁
免費預覽已結束,剩余15頁可下載查看

下載本文檔

版權說明:本文檔由用戶提供并上傳,收益歸屬內容提供方,若內容存在侵權,請進行舉報或認領

文檔簡介

1、V=19:FixmJ(與預留的自由變量數有關?。㎡:O“1.ZX 2.BX” JQ=0JLbI 1JOH仁QG :Q”L”:G”XJJD”(顯示“ L”輸入邊樁寬度, 均為正值,顯示“ XJJD”輸入邊樁測點與線路中線上點的連線 與線路中線上點的切線的夾角,如右法線為 90 度,左法線為 270 度或者-90 度,其余同理。)Prog”ZB ”:X二X+Rec(Q,T+G): Y=Y+J:”X= ” X : Pause 0 Y= ” : Y (X ”X= ” Pause 0 : YY二”:)”T= ” : T DMS (T”T二” )Goto 1(注釋:括號內寫法為 Casio4800P 計算

2、器顯示方式)2.ZB ”:坐標計算子程序V=19|FixmJ(與預留的自由變量數有關!)W=21J(與預留的自由變量數有關?。㎏:;LbI 2 J1.”XLZB ”坐標計算主程序F=0:W=27=F=39 903 |: (注釋:方框內數字 27第七段為斷鏈的段號;方框內數字 39.903 為斷鏈值,短鏈輸“ + ”,長鏈輸“-”;如遇斷鏈,依此增加即可!)E= ZW+3a+ZW+F :K ZW+3a =K V=W 1 :工 W=W+1: Goto 2 ;K 24200=K=K F: (方框為斷鏈點的后樁號,即去向樁 號,如遇斷鏈,依此增加即可!)S=K ZV+3a+1:P=(ZV+2a+1 Z

3、V+a+1) - ZV+1:匸P*S:C=ZV+4a+1:D=ZV+a+1:J二C+(l+2D)*S*90/n:L =1:X=ZV+5a+1+S*(CosC+CosJ+4刀(Cos(C+(L+0.5)l/B+2D)(L+0.5)*90*S/B/n(C+(L*l/B+2D)*L*90*S/B/n),L,1,(B-1)/6/B:),L,0,(B-1)+2E(CosY二ZV+6a+1+S*(SinC+SinJ+4刀(Sin(C+(L+0.5)I/B+2D)(L+0.5)*90*S/B/n(C+(L*I/B+2D)*L*90*S/B/n),L,1,(B-1)/6/B:),L,0,(B-1)+2E(Si

4、nT=J:3.SJK.DAT”:數據庫主程序Mcl:Defm(20+7a):Z21=?:Z21+3a=?:Z21+4a=?:Z21+5a=?:Z21+6a=?:Prog”SJK”:(輸入程序時,Defm 后須換成相應數據,a 為曲線元段數)4.SJK”: 數據庫計算子程序a二?:V=20:T=0:B二?(a:曲線元段數;B:等分數,最少為 2, 般為 4 就可以滿足要求)Lbl 0JF=0:T=6=F=39 9031: (注釋:方框內數字 6第七段減一為斷鏈的段號減去 1 所得;方框內數字 39.903 為斷鏈值,短鏈輸“ + ”,長鏈輸“-”;如遇斷鏈,依此增加即可?。〦= ZV+3a+1+

5、ZV+1+F:ZV+3a+2=E :S=ZV+1:匸ZV+2a+1 ZV+a+1:C=ZV+4a+1: D=ZV+a+1:J= C+(I+2D)*S*90/n:ZV+4a+2=J:L=1 :ZV+5a+2=ZV+5a+1+S*(CosC+CosJ+4刀Cos(C+(L+0.5)I/B+2D)(L+0.5)*90*S/B/n), L,0,(B-1)Cos(C+(L*l/B+2D)*L*90*S/B/n),L,1,(B-1)/6/B:ZV+6a+2=ZV+6a+1+S*(SinC+SinJ+4刀Sin(C+(L+0.5)I/B+2D)(L+0.5)*90*S/B/n), L,0,(B-1)Sin(

6、C+(L*I/B+2D)*L*90*S/B/n),L,1,(B-1)/6/B:V=V+1:T=T+1:T a-2 =”END OF THE PROGRAM ” Goto 0為了使程序能夠計算連續多個曲線元的坐標,減少計算過程中 的數據輸入,可利用 CASIOfx-4800/4850P 計算器提供的擴展變量設 置曲線元數據庫, 假定某一段線路包含有 a 個曲線元, 則可設置曲線 元數據庫如下表所示:曲線元數據庫設置表名稱第一段第二段第 i 段第 a-1 段第 a 段曲線兀長度212220+i19+a20+a起點曲率21+a22+a20+a+i19+2a20+2a終點曲率21+2a22+2a20+

7、2a+i19+3a20+3a起點樁號21+3a起點方位角21+4a起點x坐標21+5a起點y坐標21+6a擴展變量 Z1Z20用作自由變量,可以用來彌補字母變量的不足,此項設置可自由選擇!因此,數據庫的起點本程序從Z21開始,表中 Z21Z21+3a以及 Z21+4a、Z21+5a、Z21+6a需預先置入,表中空白部分由程序自動計算生成。CASIOfx-4800/4850P計算器程序源程序該算法共有四個程序: 1.線路坐標計算主程序( XLZB );2.數據庫生成主程序(SJK.DAT); 3坐標計算子程序(ZB); 4數據庫計算子程序(SJK)。程序使用說明不管是整體式還是分離式,主線還是匝

8、道,都能用同一個 計算器就可進行各路段相關坐標計算,只是分別采用了各自路 段的數據庫,如果再結合施工放樣正反算程序及高程程序、隧 道及路基開挖斷面程序,則可形成一個完善實用的施工測量程 序包,特別適合于公路、鐵路現場施工測量人員使用。在計算某一段線路時,首先必須運行相應路段數據庫生成主程序(SJK.DAT),然后運行線路坐標計算主程序(XLZB ),按照程序提示輸入是計算中線還是邊線,若是計算邊線除需輸入 樁號外,還需輸入寬度及斜交角度。因此,根據線路實際情況 可編制相應路段數據庫生成主程序即可完成相應路段坐標計算 工作。算例本文提供聶讓編著全站儀與高等級公路測量 306308 頁算例的 計算

9、結果如下:樁號原著坐標本算法坐標備注XYXY直OO)K0+1161378.2142822.9501378.2142822.950起點樁號:kO+116,起點方位角200,線+1401355.6612814.7411355.6612814.742xq=1378.214,yq=2822.950A + 1501346.2642811.3211346.2642811.321回A +1501346.2642811.3211346.2642811.321pa=0pb=1/124旋+2001300.1422792.1211300.1422792.121線B +2241279.8452779.3631279.

10、8452779.364圓B +2241279.8452779.3631279.8452779.364pb=1/124pc=1/124曲+2801244.1792736.8071244.1792736.808線C +341.8401230.6822677.1131230.6822677.114回C +341.8401230.6822677.1131230.6822677.114pc=1/124pd=1/60旋+3801238.6942640.0771238.6762640.065線D+407.6501254.8242617.8411254.7852617.831圓D+407.6501254.82

11、42617.8411254.7852617.831pd=1/60pe=1/60曲+4601303.3722603.2091303.3172603.148線E+495.8261335.2752618.3081335.2362618.214回E+495.8261335.2752618.3081335.2362618.214pe=1/60pf=0旋+5401357.2722656.0001357.2732655.882線F+577.4931364.6202692.7311364.6592692.605直F+577.4931364.6202692.7311364.6592692.605pf=0pg=0

12、線+6001367.9752714.9871368.0372714.857G +6401373.9382754.5401374.0412754.404名稱第一段第二段第三段第四段第五段第六段第七段曲線兀長度3474117.8465.8188.17681.66762.507起點曲率001/1241/1241/601/600終點曲率01/1241/1241/601/6000起點樁號116起點方位角200起點x坐標1378.214起點y坐標2822.950段數12345678910長度202.5262.324202.5202.5216.244202.5744.7153.6526.26153.6位置Z

13、21222324252627282930起點曲率0-1/1000-1/100001/10001/100000-1/1500-1/1500位置Z31323334353637383940終點曲率-1/1000-1/100001/10001/100000-1/1500-1/15000位置Z41424344454647484950起點樁號22258.882位置Z51525354555657585960起點方位角239-11-15.18位置Z61626364656667686970起點X坐標318841.194位置Z71727374757677787980起點Y坐標493525.938位置Z818283

14、84858687888990說明:1、該線路在第七段線元有斷鏈,具體是 K24160.097二K24200,短鏈 39.903m;2、左轉曲率為“一” ;“ + ”;3、該數據庫對應的程序見前。已知坐標求該點里程并帶入線路高程及主洞開挖放樣的復合計算程序Lbi 7IJZ1=Z5+IcosJ (Z1= I 儀器直接顯示坐標時用)Z2=Z6+ISinJ (Z2= J 儀器直接顯示坐標時用)prog“ZB”Lbi 8Z10=(X Z1)Sin(U 90) (Y Z2)Cos(U 90)AbsZ10 K二K+Z10:Pol(Z1 X,Z2 Y):Z4= : Z4= ICos(J U+90)“ K= ”

15、 : K Z8=1 = Prog “YXGC ”:半Prog “ZXGC”: Prog “ZDKW” :Goto 7:半K=K+Z10: prog “ZB” : Goto 8: 注釋:Z5 、Z6 :測站坐標(直接輸坐標時不用)I、J:實測距離及方位角(或實測坐標)JSZCXZ1 、 Z2 :實測點坐標Z10 :中間變量Z4 :實測點到線路中線的距離,負數為線左,正數為線右“ GC” :設計高程程序(子程序)“ ZDKW ”:開挖放樣子程序X、Y:實測點對應里程的設計坐標,由線路坐標計算子程序“ ZB”帶入U :實測點對應里程的切線方位角三心圓開挖放樣程序“ZDKW”(也可用于同心圓) :請自

16、己畫示意圖理解!CR :C “ H”:Z9+C A=V(Abs(Z4 士 Z8x?)+d)2+(Z9+C H h2)2) RB二Abs(Z4 士 Z8x?)+d V( (R+ r1) 2 (Z9+C H h2)2) V=Z9+C-H-h2-V(R+r1)2(Z4 士 Z8x?)2)工 Z9+C A=V(Abs(Z4 士 Z8x?) 2+(Z9+C H h1)2) RB=Abs(Z4 士 Z8x?) V(R2(Z9+CHh1)2)V=Z9+C-H-h1-V(R2(Z4 士 Z8x?)2)工 V=Z9+C H h1-V(R2 Abs(Z4 士 Z8X?) 2) 注釋:C(顯示 H):實測點至儀器中心的高差(如直接輸入高程,則將C “H”改為 C “SCGC”,且程序中的 Z9+C 就直接用 C 代替)R :開挖半徑(小)Z9 :儀高(直接輸實測高程時不用)A :實測點半徑與設計半徑之差值(徑向) , +超挖,欠挖B :實測點與設計點之差值(水平方向),+超挖,一欠挖V :實測點與設計點之豎向距離r1 :大半徑與小半徑

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯系上傳者。文件的所有權益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網頁內容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
  • 4. 未經權益所有人同意不得將文件中的內容挪作商業或盈利用途。
  • 5. 人人文庫網僅提供信息存儲空間,僅對用戶上傳內容的表現方式做保護處理,對用戶上傳分享的文檔內容本身不做任何修改或編輯,并不能對任何下載內容負責。
  • 6. 下載文件中如有侵權或不適當內容,請與我們聯系,我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論