LUẬN VĂN Đề tài : “Khảo sát phương pháp bình sai lới trắc địa tự ứng dụng x lý li thi cụng cụng trỡnh Khoa Trắc địa §å ¸n tèt nghiƯp mơc lơc Trang Mơc lơc Mở đầu Ch­¬ng Tổng quan lưới khống chế trắc địa công trình 1.1 Mét số đặc điểm lưới khống chế trắc địa công trình 1.2 Nguyên tắc ước tính lưới trắc địa công trình 1.3 Lựa chọn hệ toạ độ mặt chiếu cho lưới trắc địa công trình 13 1.4 Tính chuyển toạ độ 16 1.5 Lưới khống chế thi công công trình 22 Ch­¬ng Khảo sát phương pháp bình sai lưới tự 28 2.1 Khái niệm chung lưới trắc ®Þa tù 28 2.2 Định vị lưới 32 2.3 Một số tính chất kết bình sai lưới tự 33 2.4 NhËn xÐt vỊ b×nh sai l­íi tù 37 Ch­¬ng – øng dơng ph­¬ng pháp bình sai lưới trắc địa tự để xử lý lưới thi công công trình 39 3.1 Thuật toán xử lý số liệu lưới thi công 39 57 62 3.2 Lập chương trình bình sai 3.3 Tính toán thực nghiệm Kết luận phụ lôc 69 Tµi liƯu tham khảo Cao Bá Hạ -1- 70 87 Lớp Trắc địa B-K48 Khoa Trắc địa Đồ án tốt nghiệp Mở đầu Trong năm gần công trình xây dựng như: công trình cầu, công trình công nghiệp - dân dụng, công trình thủy lợi, thủy điện xây dựng ngày rộng rÃi Để đáp ứng yêu cầu công trình công tác trắc địa đóng vai trò quan trọng từ giai đoạn khảo sát, thiết kế, thi công công trình đến công trình bắt đầu vào sử dụng ổn định Các mạng lưới trắc địa công trình thường có tính đặc thù cao, với việc đo đạc xác (lựa chọn dụng cụ máy móc phương pháp đo) việc nghiên cứu phương pháp tính toán, quy trình xử lý số liệu cách hợp lý, phù hợp với đặc điểm chất lưới trắc địa công trình cần thiết Được hướng dẫn thầy giáo TS Nguyễn Quang Phúc đà lựa chọn đề tài tốt nghiệp: Khảo sát phương pháp bình sai lưới trắc địa tù vµ øng dơng xư lý l­íi thi công công trình Trong đề tài đà đặt khảo sát, nghiên cứu nội dung sau: Chương 1: Tổng quan lưới khống chế trắc địa công trình Chương 2: Khảo sát phương pháp bình sai lưới tự Chương 3: ứng dụng phương pháp bình sai lưới trắc địa tự để xử lý lưới thi công công trình Kết luận Hà Nội, tháng - 2006 Sinh viên Cao Bá Hạ Cao Bá Hạ -2- Lớp Trắc địa B-K48 Khoa Trắc địa Đồ ¸n tèt nghiƯp Ch­¬ng Tỉng quan vỊ l­íi khèng chế trắc địa công trình 1.1 Một số đặc điểm lưới khống chế trắc địa công trình 1.1.1 Lưới khống chế mặt Lưới khống chế lập khu vực xây dựng công trình công nghiệp, thành phố, khu vực xây dựng cầu cảng, đường hầm sở trắc địa phục vụ cho khảo sát thiết kế, thi công xây dựng cho khai thác sử dụng công trình Theo mục đích thành lập lưới trắc địa công trình phân thành nhóm: Lưới khảo sát công trình, lưới thi công công trình, lưới quan trắc chuyển dịch biến dạng công trình Trong giai đoạn khảo sát , người ta phải nghiên cứu tổng hợp điều kiện tự nhiên vùng xây dựng, thu thập số liệu địa hình, địa mạo, địa chất công trình, địa chất thủy văn, khí tượng,để có giải pháp kỹ thuật hợp lý thiết kế xây dựng công trình Đề xuất yêu cầu tiến hành thành lập lưới khống chế trắc địa khu vực, đo vẽ đồ địa hình tỉ lệ trung bình, đo vẽ mặt cắt địa hình theo tuyến công trình, phục vụ cho giai đoạn thiết kế giai đoạn thiết kế cần lập đồ địa hình, mặt cắt tỷ lệ lớn để thiết kế chi tiết Chuẩn bị phương án trắc địa để chuyển thiết kế thực địa Trong giai đoạn thi công, phải tiến hành công tác xây dựng lưới sở trắc địa phục vụ cắm công trình, chuyển thiết kế thực địa vị trí, kích thước đà thiết kế Mặt khác phải theo dõi thi công hàng ngày để đảm bảo công trình có kết cấu thiết kế Sau hoàn thành công trình cần đo vẽ hoàn công để kiểm tra vị trí, kích thước công trình đà xây dựng Công tác trắc địa phục vụ quan trắc biến dạng công trình gồm đo độ lún móng, xác định độ xê dịch vị trí mặt độ nghiêng công trình nước ta lưới trắc địa công trình phận hệ toạ độ quốc gia VIệT NAM mô tả (Hình1.1) Cao Bá Hạ -3- Lớp Trắc địa B-K48 Khoa Trắc địa Đồ án tốt nghiệp Lưới trắc địa mặt Hình 1.1 Phân loại lưới khống chế trắc địa mặt Cao Bá Hạ -4- Lớp Trắc địa B-K48 Lưới quan trắc biến dạng Lưới thi công CT Lưói TĐCT Lưới khảo sát CT Địa II Địa I Giao hội Tam giác nhỏ Đường chuyền Địa sở Lưói địa Lưói đo vẽ Đa giăc Giải tích II Hạng IV Hạng III Hạng II Hạng I Giải tích I Lưói khu vực Lưói nhà nước Khoa Trắc địa Đồ án tốt nghiệp Tuy mục đích thành lập có khác nhìn chung, lưới trắc địa công trình chủ yếu lập theo phương pháp truyền thống đà biết như: phương pháp tam giác, đa giác hay giao hội Ngoài lưới trắc địa công trình thành lập theo phương pháp đặc biệt lưới tứ giác không đường chéo, lưới ô vuông xây dựng, lưới tam giác nhỏ đo oàn cạnh độ xác cao, lưới tam giác bẹt Hiện nay, công nghệ GPS đà bước ứng dụng trắc địa công trình nước ta nhìn chung, việc lập lưới trắc địa công trình trị đo mặt đất chiếm vị trí chủ yếu Độ xác lưới trắc địa công trình quy định Tuỳ thuộc vào yêu cầu nhiệm vụ cụ thể mà phải giải tuỳ theo giai đoạn khảo sát thiết kế, thi công xây dựng khai thác sử dụng Vị trí mật độ số lượng điểm khống chế tuỳ thuộc mục đích thành lập đặc điểm kỹ thuật loại công trình Theo độ xác lưới trắc địa công trình phân thành hai trường hợp: + Trường hợp 1: Lưới trắc địa công trình có độ xác tương đương lưới đo vẽ đồ (thường gặp thời kỳ khảo sát công trình để đo vẽ địa hình công trình) Trong trường hợp sở để ước tính sai số trung phương bậc lưới cấp cuối (lưới đo vẽ) so với lưới khu vực nhà nước không vượt 0, mm.M Lưới trắc địa công trình lúc phát triển dựa điểm lưới nhà nước theo nguyên tác chung từ tổng quát đến chi tiết + Trường hợp 2: Lưới trắc địa công trình có yêu cầu độ xác cao hẳn so với lưới đo vẽ đồ (thường gặp giai đoạn thi công, sử dụng công trình) Trong trường hợp cần phải lập lưới chuyên dùng cho công trình Trong trường hợp thứ lưới trắc địa công trình phát triển theo nguyên tắc chung từ tổng quát đến chi tiết lấy điểm khống chế nhà nước làm sở (coi điểm sai số) Cơ sở ước tính độ xác cho lưới trắc địa công trình lúc sai số trung phương vị trí điểm yếu bậc lưới cấp cuối so với điểm lưới bậc cao không 0, mm.M Cao Bá Hạ -5- Lớp Trắc địa B-K48 Khoa Trắc địa Đồ án tốt nghiệp Trong trường hợp thứ hai điểm lưới nhà nước có khu vực sử dụng số liệu để định hướng định vị cho lưới trắc địa công trình 1.1.2 L­íi khèng chÕ ®é cao L­íi khèng chÕ ®é cao lập khu vực xây dựng công trình sở trắc địa phục vụ cho đo vẽ địa hình công trình, cho thi công công trình cho quan trắc chuyển dịch biến dạng công trình Lưới độ cao trắc địa công trình cố thể thành lập theo dạng sau: Phương pháp thủy chuẩn hình học tia ngắm ngắn, phương pháp đo cao lượng giác tia ngắm ngắn, phương pháp thủy chuẩn thủy tĩnh Chính mục đích hành lập nên lưới độ cao trắc địa công trình có đặc điểm khác so với lưới độ cao nhà nước: Thứ nhất: Cấp hạng lưới khống chế độ cao quy định tuỳ thuộc vào diện tích khu vực xây dựng công trình: Bảng 1.1 Diện tích xây dựng Cấp hạng thủy chuÈn > 500 km  500 km 10  50 km 50 I, II, III, IV II, III, IV III, IV < 10 km IV Thø hai: để phục vụ cho đo vẽ địa hình công trình lưới độ cao trắc địa công trình phát triển dựa điểm lưới độ cao nhà nước theo nguyên tắc từ tổng quát đến chi tiết Thứ ba: để thi công công trình, lưới độ cao cần phải xây dựng tuỳ thuộc vào đặc điểm yêu cầu kỹ thuật loại công trình, điểm quan trọng phải đo nối với lưới độ cao nhà nước Thứ tư: so với lưới nhà nước mật độ điểm lưới trắc địa công trình dày hơn, chiều dài rút ngắn Để thấy rõ ta tìm hiểu số tiêu lưới độ cao trắc địa công trình: Cao Bá Hạ -6- Lớp Trắc địa B-K48 Khoa Trắc địa Đồ án tốt nghiệp Bảng 1.2 Chỉ tiêu kỹ tht CÊp h¹ng thđy chn H¹ng II H¹ng III H¹ng IV 40/270 15/65 4/15 10/140 3/25 2/6 - Khu vùc x©y dùng: 0.2 0.2  0.5 - Khu vùc ch­a x©y dùng: 0.8 0.5  2.0 L 10 L 20 L (mm) (mm) (mm) ChiỊu dµi lín (km): - Giữa điểm gốc: (công trình/nhà nước) - Giữa điểm nút: (công trình/nhà nước) Khoảng cách mốc thuỷ chuẩn (km): Sai số khép giới hạn tuyến: ( L tính km) 1.2 Nguyên tắc ước tính lưới trắc địa công trình Trên sở tài liệu thu thập yêu cầu độ xác bố trí công trình ta dự tính độ xác cho lưới tiến hành ước tính độ xác lưới Có hai phương pháp ước tính độ xác yếu tố mạng lưới phương pháp ước tính chặt chẽ phương pháp gần 1.2.1 Phương pháp ước tính gần Trước kia, phương tiện phục vụ công tác đo đạc tính toán lưới chưa nhiều người ta thường thiết kế lưới theo số dạng đồ hình mẫu chuỗi tam giác, lưới đa giác trung tâm, lưới đường chuyền sử dụng công thức lập sẵn sở toán bình sai điều kiện tập trung xem xét số yếu tố đặc trưng vị trí yếu mạng lưới Ví dụ ước tính độ xác chuỗi tam giác, đa giác trung tâm, người ta thiết kế lưới theo dạng đồ hình mẫu đơn giản sủ dụng công thức lập sẵn theo dạng đồ hình Sau số ví dụ Cao Bá Hạ -7- Lớp Trắc địa B-K48 Khoa Trắc địa Đồ án tốt nghiệp 1.2.1.1 ước tính độ xác chuỗi tam giác E S F Hình 1.2 Đồ hình chuỗi tam giác Với chuỗi tam giác có dạng Hình 2.1 độ xác yếu tố vị trí yếu đặc trưng sai số: Sai số trung phương chiều dài cạnh yếu: Ký hiệu: m S sai số trung phương chiều dài cạnh EF tính theo đường (1) m S sai số trung phương chiều dài cạnh EF tÝnh theo ®­êng (2) Lóc ®ã träng sè chiều dài cạnh EF tính theo nguyên tắc cộng träng sè m s21  m s22 2 2 P EF = P + P = + =  2 m s1 m s2 m s1 m s2 Sai sè trung ph­¬ng chiỊu dài cạnh EF: mS = EF mlg s1 =  PEF =± m s2 m s2 (1.1) m s21  m s22 2 2 2 m  ( A   B   C ) = m   R 3 m S , m S xác định thông qua sai số trung phương loga chiều dài cạnh: m Si  mlg s S (1.2)  10 μ = m  víi m  lµ sai sè trung phương đo góc dự kiến Sai số trung phương phương vị cạnh yếu: đo hai đầu phương vị yếu phương vị cạnh Tương tự trên, gọi m , m sai số trung phương phương vị cạnh EF tính theo đường (1) (2) ta tính sai số trung phương phương vị cạnh yếu m là: EF Cao Bá Hạ -8- Lớp Trắc địa B-K48 Khoa Trắc địa Đồ án tốt nghiệp m EF = m1 m (1.3) 2 m1  m Sai số trung phương tương hỗ vị trí điểm: m 2 ) s  mth   m s2 ( (1.4) 1.2.1.1 ước tính độ xác đa giác trung tâm D E S c F B A Hình 1.3 Đồ hình đa giác trung tâm Đối với đa giác trung tâm hình 1.3, việc ước tính ®é chÝnh x¸c cã thĨ thùc hiƯn theo nh­ ­íc tính độ xác chuỗi tam giác đơn trường hợp xuất phát từ cạnh ®¸y AB = b , tøc: mlg b1  mlg b2  mlg b0 m d  m C m Do cạnh DE = S, ta cã: 2 mlg s1 mlg s2 S mS   mlg b0 2 M 10 mlg s1  mlg s2 Vµ m  Cao Bá Hạ 2 m1 m m m 1 2 (1.5) (1.6)  m -9- Lớp Trắc địa B-K48 Khoa Trắc địa 1130 1118 1117 1116 Đồ án tốt nghiệp Do 1116 I=1,Ian Do 1117 K=I,Ian J=(K-1)*K/2+I Do 1118 L=1,3 If(Ng.EQ.Ndi) Goto 1130 CC(J)=CC(J)+C(I,L)*C(K,L) BB(J)=BB(J)+B(I,L)*B(K,L) Continue Continue Continue KL=(Ian-1)*Ian/2+Ian Do 1131 L=1,KL If(Ng.EQ.Ndi) Goto 132 PTC(L)=PTC(L)+CC(L)*10000000.d0 Goto 1131 132 PTC(L)=PTC(L)+BB(L)*10000000.d0 1131 Continue Ls=Ian-1 Call DTDNG(Ls) Call DNDAO(Ian) 800 805 806 802 801 889 701 Do 800 I2=1,Ian RLV(I2)=0.d0 Continue Do 801 I3=1,Ian Do 802 J3=1,Ian If(J3.GT.I3) Goto 805 K=(I3-1)*I3/2+J3 Goto 806 K=(J3-1)*J3/2+I3 RLV(I3)=RLV(I3)-PTC(K)*PTC(KL+J3) Continue Continue Write(4, 701) Format(//14x, 'III.KET QUA PHAN TICH DO ON DINH CUA CAC MOC '/ * 15x, '========================================= '// *10x,69('=')/ *10x,'| S | TEN | TOA DO | lech k/c * |'/ *10x,'| T | | -| -* -|'/ *10x,'| T | DIEM | dX(m) | dY(m) | ds(cm) * |'/ *10x,'| -| -| -| -| -* -|') sdy=0.d0 DRmax=0d0 Do 989 K=(nxd+1),Ndi RX= RLV(K*2-1) RY= RLV(K*2) DR= dsqrt(RX**2+RY**2) If(DR.le.DRmax) Goto 710 DRmax=DR Tendy=Khd(K) sttmax=K Xmax=X(K) Cao Bá Hạ - 73 - Lớp Trắc địa B-K48 Khoa Trắc địa 710 711 989 712 820 513 812 832 833 830 831 Đồ án tốt nghiệp Ymax=Y(K) Write (4, 711) k,Khd(k),RX,RY,DR Format (10x,'|', I2,' |' , A7,'|', F15.3,'|', F15.3, '| ',F15.3,' * |') Continue write(4,712) Format(10x,69('=')) Write(4, 820) Format(//15x, 'IV.KET QUA TRI DO SAU BINH SAI'/ * 15X, '==============================='//) If(Ngoc.eq.0.d0) Goto 514 Write(4,513) Format(15X, 'A.TRI DO GOC NGANG SAU BINH SAI'// *10x,69('=')/ *10x,'| SO | TEN-DINH ' *'| GIA TRI GOC | HIEU | TRI GOC SAU |'/ *10x,'| THU| -' *'| DO | CHINH | BINH SAI |'/ *10x,'| TU | Trai-Giua-Phai ' *'| Do-Phut-Giay| Giay | Do-Phut-Giay|'/ *10x'| | -| -| ' *'| -|') Pvv1=0.d0 Do 810 M= 1,Ngoc Pg=1.d0/(Ssg**2) I=Ig(M) J=It(M) K=Ip(M) Call dSHCG(Ro,Pi2,I,J,K,M,Ian) V= 0.d0 Do 812 L= 1, Ian V=V + HV(L)*RLV(L) Continue V=V+HV(Ian+1) Pvv1=Pvv1+V*V*Pg Call Drado(Goc(M),Id(m),Iph(m),giay(m),Ro) Id1(m) = Id(m) Iph1(m) = Iph(m) Giay1(m)= Giay(m)+V If(Giay1(m).ge.60.d0) Goto 830 If(Giay1(m).lt.0.d0) Goto 832 Goto 860 II = Giay1(m)/60.d0-1 Iph1(m) = Iph1(m)+II Giay1(m)= Giay1(m)+Abs(II*60.d0) If(Iph1(m).le.0.d0) Goto 833 Goto 860 II = Iph1(m)/60.d0-1 Id1(m) = Id1(m)+II Iph1(m)= Iph1(m)+Abs(II*60.d0) Goto 860 II = Giay1(m)/60.d0 Iph1(m) = Iph1(m)+II Giay1(m)= Giay1(m)-II*60.d0 If(Iph1(m).ge.60) Goto 831 Goto 860 II = Iph1(m)/60.d0 Id1(m) = Id1(m)+II Iph1(m)= Iph1(m)-II*60.d0 Goto 860 Cao Bá Hạ - 74 - Lớp Trắc địa B-K48 Khoa Trắc địa 860 Đồ án tốt nghiệp Write (4,813) M,Khd(It(M)),Khd(Ig(M)),Khd(Ip(M)) ,Id(m), Iph(m) , Giay(m), V, Id1(m), Iph1(m), Giay1(m) Format ( 10x,'|', I2, ' |', 3A7, '|', I4, I3, F6.2, '|', *F8.2,4X,'|', I4, I3, F6.2, '|' ) Continue Write(4, 890) Format(10x,69('=')) If (Ncanh.eq.0.d0) Goto 515 Write(4,920) Format(//15X, 'B.TRI DO CANH SAU BINH SAI'// *10x, 59('=')/ *10x,'| SO | TEN-CANH | CHIEU DAI | HIEU | CHIEU ' *'DAI |'/ *10x,'| THU| | DO | CHINH | SAU ' *'BINH SAI |'/ *10x,'| TU | DAU-CUOI | (m) | (m) | (m)' *' |'/ *10x,'| | | | -| ' *' |') Pvv2=0.d0 Do 910 M=1,(Ncanh) I=Kd(M) J=Kc(M) Pc=1/((a+b1*Canh(m))**2) Call dSHCC(M,I,J,Ian) V=0.d0 Do 912 L=1,Ian V=V+HV(L)*RLV(L) Continue V=V+HV(Ian+1) Pvv2=Pvv2+V*V*Pc Canhbs=Canh(M)+V Write (4, 913) M,Khd(Kd(M)),Khd(Kc(M)),Canh(M),V,Canhbs Format (10x,'|',I2,' |' , A7, ' ', A7, ' |', F10.3, ' |', F6.3 * ' |', F10.3, ' |') Continue Write(4, 914) Format(10x,59('=')) If(Nfv.eq.0.d0) Goto 1222 Write(4,516) Format(//15X, 'C.TRI DO PHUONG VI SAU BINH SAI'// *10x, 64('=')/ *10x,'| SO | TEN-HUONG | GIA TRI GOC | HIEU |' *' TRI GOC SAU |'/ *10x,'| THU| | DO | CHINH |' *' BINH SAI |'/ *10x,'| TU | DAU-CUOI | Do-Phut-Giay| Giay |' *' Do-Phut-Giay|'/ *10x,'| | | -| |' *' -|') Pvv3=0.d0 Do 310 M= 1,Nfv Pfv=1/(Ssfv**2) I=Idfv(M) J=Icfv(M) Call dSHCFV(Ro,Pi2,I,J,M,Ian) V= 0.d0 Do 312 L= 1, Ian V=V + HV(L)*RLV(L) Continue V=V+HV(Ian+1) * 813 810 890 514 919 920 912 913 910 914 515 516 312 Cao B¸ Hạ - 75 - Lớp Trắc địa B-K48 Khoa Trắc ®Þa 332 333 330 331 360 313 310 390 1222 1939 1989 1889 1701 Đồ án tốt nghiệp Pvv3=Pvv3+V*V*Pfv Call Drado(Gocfv(M),Id(m),Iph(m),giay(m),Ro) Id1(m) = Id(m) Iph1(m) = Iph(m) Giay1(m)= Giay(m)+V If(Giay1(m).ge.60.d0) Goto 330 If(Giay1(m).lt.0.d0) Goto 332 Goto 360 II = Giay1(m)/60.d0-1 Iph1(m) = Iph1(m)+II Giay1(m)= Giay1(m)+Abs(II*60.d0) If(Iph1(m).le.0.d0) Goto 333 Goto 360 II = Iph1(m)/60.d0-1 Id1(m) = Id1(m)+II Iph1(m)= Iph1(m)+Abs(II*60.d0) Goto 360 II = Giay1(m)/60.d0 Iph1(m) = Iph1(m)+II Giay1(m)= Giay1(m)-II*60.d0 If(Iph1(m).ge.60) Goto 331 Goto 360 II = Iph1(m)/60.d0 Id1(m) = Id1(m)+II Iph1(m)= Iph1(m)-II*60.d0 Goto 360 Write (4,313) M,Khd(Idfv(M)),Khd(Icfv(M)) * ,Id(m), Iph(m) , Giay(m), V, Id1(m), Iph1(m), Giay1(m) Format ( 10x,'|', I2, ' |', A7, ' ', A7, ' |', I4, I3, * F6.2, '|', F8.2,4X,'|', I4, I3, F6.2, '|' ) Continue Write(4, 390) Format(64('=')) PVV=Pvv1+Pvv2+Pvv3 Goto 1939 PVV=Pvv1+Pvv2 SSDV=dsqrt(PVV/(Ngoc+Ncanh+Nfv-2*Ndi+3)) Do 1989 K3=1,Ndi X(K3)=X(K3)+RLV(K3*2-1) Y(K3)=Y(K3)+RLV(K3*2) Continue Write(4, 1701) Format(//15x, 'V.KET QUA TOA DO DIEM SAU BINH SAI '/ * 15x, '==================================== '// *10x,91('=')/ *10x,'| S | TEN | TOA DO | SAI SO VI' *' TRI DIEM | Elip sai so |'/ *10x,'| T | | -| ' *' | |'/ *10x,'| T | DIEM | X(m) | Y(m) |Mx(cm)|' *'My(cm)|Md(cm)| E | F | goc dinh huong |'/ *10x,'| -| -| | | |' *' | | |') Do 1710 I=1,Ndi Kx =(I*2-1)*(I*2-2)/2+I*2-1 Ky =(I*2)*(I*2-1)/2+I*2 Kxy=Ky-1 Qx =PTC(Kx) Qy =PTC(Ky) Qxy=PTC(Kxy) Cao Bá Hạ - 76 - Lớp Trắc địa B-K48 Khoa Trắc địa Đồ án tốt nghiÖp SSx=100.d0*SSDV*Dsqrt(ABS(Qx)) SSy=100.d0*SSDV*Dsqrt(ABS(Qy)) SSd(I)=Dsqrt(SSx*SSx+SSy*SSy) ZM(I)=Qx-Qy ZN(I)=2*Qxy ZZ=Qx+Qy ZR=Dsqrt(ZM(I)*ZM(I)+ZN(I)*ZN(I)) E(I)=SSDV*Dsqrt((ZZ+ZR)/2) F(I)=SSDV*Dsqrt((ZZ-ZR)/2) E1=1000.d0*E(I) F1=1000.d0*F(I) Ph(I)=dMod(dAtan2(ZN(I),ZM(I))+Pi2,Pi2) If(Ph(I).Lt.0.d0) Ph(I)=Ph(I)+Pi2 If(Ph(I).GT.Pi) Ph(I)=Ph(I)-Pi Ph1=Ph(I) Call Drado(Ph1,Id,Iph,giay,Ro) 1711 1710 1712 1021 1020 80 Write (4, 1711) I,Khd(I),X(I),Y(I),SSx,SSy,SSd(I),E1,F1,Id(I) *,Iph(I) Format (10x,'|', I2,' |' , A7,'|', F12.3,'|', F12.3, '| ', *F4.2, ' | ', F4.2, ' | ', F4.2, ' |',F5.2,' |',F5.2,' | ',I4,' *| ',I3,' |') Continue Write(4, 1712) Format(10x,91('=')) Do 1020 M=1, Ncanh I=Kd(M) J=Kc(M) Dx1= X(J)-X(I) Dy1= Y(J)-Y(I) Ds(m)= Dsqrt(Dx1**2+Dy1**2) Do 1021 L=1,Ian HV(L)= 0.d0 Continue HV(I*2-1)= -Dx1/Ds(m) HV(I*2) = -Dy1/Ds(m) HV(J*2-1)= Dx1/Ds(m) HV(J*2) = Dy1/Ds(m) QF= 0.d0 Call Dhts(Ian, QF) Sm(m)= SSDV*Dsqrt(ABS(QF)) Ktt(m)=(Ds(m)/Sm(m)) Ktt(m)=(Ktt(m)/1000.d0) Ktt(m)=Ktt(m)*1000.d0 Continue KI=1000000000 Do 80 I=1,Ncanh If(Ktt(I).GE.KI) Goto 80 KI=Ktt(I) Kj=Kd(I) Kn=Kc(I) Continue Do 760 m=1,Ncanh I=Kd(M) J=Kc(M) Dx1= X(J)-X(I) Dy1= Y(J)-Y(I) Ds(m)= Dsqrt(Dx1**2+Dy1**2) Do 1028 L=1,Ian HV(L)= 0.d0 Cao B¸ Hạ - 77 - Lớp Trắc địa B-K48 Khoa Trắc địa Đồ án tốt nghiệp 1028 Continue HV(I*2-1)= Ro*dy1/(Ds(m)**2) HV(I*2) = -Ro*dx1/(Ds(m)**2) HV(J*2-1)= -Ro*dy1/(Ds(m)**2) HV(J*2) = Ro*dx1/(Ds(m)**2) QF=0.d0 Call Dhts(Ian,QF) Sa(m)=SSDV*Dsqrt(ABS(QF)) 760 Continue ss=0.d0 Do 341 I=1,Ncanh If(Sa(I).LE.ss) Goto 341 ss=Sa(I) Ko=Kd(I) Kz=Kc(I) 341 Continue Sd=0.d0 Do 342 I=1,Ndi If(SSd(I).LE.Sd) Goto 342 Sd=SSd(I) Km=I 342 Continue Write(4,406) 406 Format(//15x,'VI CAC YEU TO TUONG HO CUA LUOI'/ * 15x,'================================'// *10x,82('=')/ *10x,'| No | Ten canh | Chieu dai | Phuong vi |Ma(")|' *' Ms(m)| Ms/S | Mth |'/ *10x,'| | | -| -| -|' *' | | |') Do 400 M=1,Ncanh I=Kd(m) J=Kc(m) DX1=X(J)-X(I) DY1=Y(J)-Y(I) S(m)=dsqrt(DX1**2+DY1**2) tu=(S(m)*Sa(m))/Ro Sth(m)=dsqrt(Sm(m)*Sm(m)+Tu**2) afa=dMod(dAtan2(DY1,DX1)+Pi2,Pi2) If(afa.Lt.0.d0) afa=afa+Pi2 Call Drado (afa,Id,Iph,giay,Ro) Write(4,405) M,Khd(Kd(m)),Khd(Kc(m)),S(m),Id(m),Iph(m),giay(m), *Sa(m),Sm(m),Ktt(m),Sth(m) 405 Format(10x,'|',I3,' |',A7,A7,'|',F11.3,'|',I4,I3,F6.2,'|', *F5.2,'|',F6.3,'|',' 1:',I8,' |',F6.3,'|') 400 Continue Write(4,407) 407 Format(10x,82('=')) Write(4, 2001) 2001 Format(/10x,'VII.KET LUAN'/ * 10x,'============') Write(4,1027) SSDV 1027 Format(/15x, * '1 Sai so trung phuong so don vi mo =', * F6.2, ' " ') Write (4,348) Khd(Ko),Khd(Kz),ss Format(15x,'2 Phuong vi yeu nhat * A7,'- ',A7,' ma =',F6.2,' " ') Write (4,1119) Khd(Kj),Khd(Kn),KI 1119 Format(15x,'3 Canh yeu nhat 348 Cao Bá Hạ - 78 - ', ', Lớp Trắc địa B-K48 Khoa Trắc ®Þa 453 566 69 93 83 162 131 137 138 821 839 721 739 957 961 Đồ án tốt nghiệp *A7,'- ',A7,'ms/S = 1/',I6) Write(4,453) Khd(Km),Sd Format(15x,'4 Diem yeu nhat ',A7,' * F6.2,'(cm)') Write(4,566) Format(///40x,' Nguoi thuc hien: Cao Ba Ha '/ * 40x,' Lop : Trac dia B - K48'/ * 40x,' TRUONG DAI HOC MO - DIA CHAT '/) M =', Write(6,69) Format(' 0'/'SECTION'/' 2'/'ENTITIES') Do 83 I=1,Ncanh Write(6,93) Y(Kd(I)),X(Kd(I)),Y(Kc(I)),X(Kc(I)) Format( *' 0'/'LINE'/' 8'/'CANH'/' 62'/'131'/' 10'/f18.3/' 20'/f18.3/ *' 30'/'0.0'/' 11'/f18.3/' 21'/f18.3/' 31'/'0.0') Continue Do 131 I=1,Ngoc Write(6,162) Y(It(I)),X(It(I)),Y(Ig(I)),X(Ig(I)) Format( *' 0'/'LINE'/' 8'/'CANH'/' 62'/'131'/' 10'/f18.3/' 20'/f18.3/ *' 30'/'0.0'/' 11'/f18.3/' 21'/f18.3/' 31'/'0.0') Continue Do 138 I=1,Ngoc Write(6,137) Y(Ig(I)),X(Ig(I)),Y(Ip(I)),X(Ip(I)) Format( *' 0'/'LINE'/' 8'/'CANH'/' 62'/'131'/' 10'/f18.3/' 20'/f18.3/ *' 30'/'0.0'/' 11'/f18.3/' 21'/f18.3/' 31'/'0.0') Continue Do 839 I=(Nxd+1),Ndi Write(6,821) Y(I),X(I),Cch,Khd(I) Format(' 0'/'TEXT'/' 8'/'TEN_DIEM'/' 62'/'30'/' 10'/F18.3/ *' 20'/F18.3/' 30'/'0.0'/' 40'/F8.1/' 1'/A7/' 41'/'0.7'/ *' 7'/'STANDARD') Continue Do 739 k=1,Nxd Write(6,721) Y(k),X(k),Cch,Khd(k) Format(' 0'/'TEXT'/' 8'/'TEN_DIEM'/' 62'/'50'/' 10'/F18.3/ *' 20'/F18.3/' 30'/'0.0'/' 40'/F8.1/' 1'/A7/' 41'/'0.7'/ *' 7'/'STANDARD') Continue Cd=10.d0 Do 961 I=1,Ncanh I1=Kd(I) J1=Kc(I) dx11 =X(J1)-X(I1) dy11 =Y(J1)-Y(I1) afS = dMod(dAtan2(dy11,dx11)+Pi2,Pi2) If(afS.LT.0) afS = afS + Pi2 Xg(I)=(X(I1)+X(J1))/2 Yg(I)=(Y(I1)+Y(J1))/2 Xq1=Xg(I)+cd*dcos(afS+Pi/2) Yq1=Yg(I)+cd*dsin(afS+Pi/2) Xq2=Xg(I)+cd*dcos(afS+(3*Pi)/2) Yq2=Yg(I)+cd*dsin(afS+(3*Pi)/2) Write(6,957) Yq1,Xq1,Yq2,Xq2 Format( *' 0'/'LINE'/' 8'/'KH-CANHDO'/' 62'/'120'/' 10'/f14.3/ *' 20'/f14.3/' 30'/'0.0'/' 11'/f14.3/' 21'/f14.3/' 31'/'0.0') Continue Cao Bá Hạ - 79 - Lớp Trắc địa B-K48 Khoa Trắc địa 849 847 201 654 767 768 769 770 §å ¸n tèt nghiƯp dl =25.d0 Do 847 I=1,Ngoc dxt=X(It(I))-X(Ig(I)) dyt=Y(It(I))-Y(Ig(I)) dxp=X(Ip(I))-X(Ig(I)) dyp=Y(Ip(I))-Y(Ig(I)) afat= dMod(dAtan2(dxt,dyt)+Pi2,Pi2) afap= dMod(dAtan2(dxp,dyp)+Pi2,Pi2) afat= afat*180.d0/pi afap= afap*180.d0/pi Write(6,849) Y(Ig(I)),X(Ig(I)),dl,afap,afat Format(' 0'/'ARC'/' 8'/'KH-GOC'/' 62'/'6'/' 10'/F18.3/' 20' * /F18.3/' 30'/' 0.0'/' 40'/F8.1/' 50'/F11.6/' 51'/F11.6) Continue Do 654 I=1,Ndi Write(6,201) Y(I),X(I),bk FORMAT(' 0'/'CIRCLE'/' 8'/'VT-diem'/' 62'/'13'/' 10'/F12.3/ *' 20'/F12.3/' 30'/'0.0'/' 40'/F10.3) Continue Do 766 I=1,Ndi Yf=F(I)/50 Do 767 K=2,50 K1=K-1 Ye(K)=-K1*Yf Ca1=E(I)*dsqrt(F(I)*F(I)-Ye(K)*Ye(K)) Xe(K)=Ca1/F(I) Se(K)=dsqrt(Xe(K)*Xe(K)+Ye(K)*Ye(K)) Continue Do 768 J=52,100 J1=J-51 Ye(J)=-(F(I)-J1*Yf) Ca2=E(I)*dsqrt(F(I)*F(I)-Ye(J)*Ye(J)) Xe(J)=-Ca2/F(I) Se(J)=dsqrt(Xe(J)*Xe(J)+Ye(J)*Ye(J)) Continue Do 769 J2=102,150 J3=J2-101 Ye(J2)=J3*Yf Ca3=E(I)*dsqrt(F(I)*F(I)-Ye(J2)*Ye(J2)) Xe(J2)=-Ca3/F(I) Se(J2)=dsqrt(Xe(J2)*Xe(J2)+Ye(J2)*Ye(J2)) Continue Do 770 K2=152,200 K3=K2-152 Ye(K2)=(F(I)-K3*Yf) Ca4=E(I)*dsqrt(F(I)*F(I)-Ye(K2)*Ye(K2)) Xe(K2)=Ca4/F(I) Se(K2)=dsqrt(Xe(K2)*Xe(K2)+Ye(K2)*Ye(K2)) Continue Xe(1) = E(I) Ye(1) = 0.d0 Se(1) = E(I) Xe(51)= 0.d0 Ye(51)=-F(I) Se(51)= F(I) Xe(101)=-E(I) Ye(101)= 0.d0 Se(101)= E(I) Xe(151)= 0.d0 Cao Bá Hạ - 80 - Lớp Trắc địa B-K48 Khoa Trắc địa 771 773 772 766 209 Đồ án tốt nghiệp Ye(151)= F(I) Se(151)= F(I) Aft=dmod(dAtan2(Ye(1),Xe(1))+Pi2,Pi2) Do 771 Im=1,200 Afp=dmod(dAtan2(Ye(Im),Xe(Im))+Pi2,Pi2) Beta(Im)= Afp-Aft If(Beta(Im).Lt.0.d0) Beta(Im)=Beta(Im)+Pi2 Se(Im)=15000.d0*Se(Im) Af1=Ph(I)+Beta(Im) X2(Im)=X(I)+Se(Im)*dcos(Af1) Y2(Im)=Y(I)+Se(Im)*dsin(Af1) Continue Do 772 I6=1,200 X2(201)=X2(1) Y2(201)=Y2(1) Write(6,773) Y2(I6),X2(I6),Y2(I6+1),X2(I6+1) Format( *' 0'/'LINE'/' 8'/'ELIP-SAISO'/' 62'/'10'/' 10'/f14.3/ *' 20'/f14.3/' 30'/'0.0'/' 11'/f14.3/' 21'/f14.3/' 31'/'0.0') Continue Continue Write(6,209) Format(' 0'/'ENDSEC'/' 0'/'EOF') Close(6) Stop End Chương trình bình sai lưới độ cao thi c«ng $debug $large c Binh sai luoi cao thi cong c Ho va ten : Cao Ba Ha c Lop : Trac dia BK48 Implicit Real *8(a-h, o-z) Character ff*50, Tenct*50, Tr*7 Character Khd(200)*8 Dimension H(200),hh(200),S(200),Kd(500),Kc(500),RLV(100),RH(100) Common /v1/ HV(200), PTC(5000) Write (*,*) 'Nhap ten file so lieu : ' Read (*,'(a)') ff Open (1, file=ff, status='old') Write (*,*) 'Nhap ten file ket qua : ' Read (*,'(a)') ff Open (4, file=ff, status='new') Read (1,'(a)') Tenct Read (1,*) Ng, Nxd, Ncc,Gh a= -10000000000.d0 Ndi= (Ng+Nxd) Do I= 1, Ndi H(I)= a Continue Read(1,89) (Khd(m),m=1,Ndi) 89 Format(10A6) Read(1,*) (H(k),K=(Nxd+1),Ndi) Do I= 1, Ncc Read (1,*) Kd(I), Kc(I), hh(I), S(I) Continue Write (4,5) Ng, Nxd, Ncc Format (//15x, 'KET QUA BINH SAI LUOI DO CAO '/ * 15x, '===========* * *============'// * 15x, 'I.SO LIEU KHOI TINH'/ * 15x, '==================='/ * 15x, 'So diem cao goc : ',I3/ * 15x, 'So diem can xac dinh : ',I3/ Cao Bá Hạ - 81 - Lớp Trắc địa B-K48 Khoa Trắc địa * 14 11 12 13 10 20 * * * * * * 22 21 23 199 101 100 66 65 71 70 80 §å ¸n tèt nghiÖp 15x, 'So luong chenh cao : ',I3// ) LT= Do 10 I= 1, Ncc If (H(Kd(I)).ne.a.and.H(Kc(I)).eq.a) Goto 11 If (H(Kd(I)).eq.a.and.H(Kc(I)).ne.a) Goto 12 Goto 10 H(Kc(I))= H(Kd(I))+hh(I) Goto 13 H(Kd(I))= H(Kc(I))-hh(I) LT= LT+1 If (LT.eq.Nxd) Goto 20 Continue Goto 14 Write (4,6) 'II DO CAO DIEM DINH VI' Format (15x, A23/ 15X, 23('=')// 10x, 44('-')/ 10x, '| S | TEN | DO CAO |'/ 10x, '| T | | |'/ 10x, '| T | DIEM | H(mm) Ghi chu |'/ 10x, '| |') Do 21 J= (Nxd+1), Ndi Write (4,22) J, Khd(J), H(J) Format (10x,'|', I3,' |', A7,'|', F9.2,' Continue Write(4,23) Format(10x,44('-')) Ian=Ndi KL=(Ian-1)*Ian/2+Ian Ms=Ian*(Ian+1)/2 Do 199 K=1,Ms PTC(K)=0.d0 Continue Do 100 M= 1, Ncc Do 101 K= 1, (Ian+1) HV(K)=0d0 continue I= Kd(M) J= Kc(M) HV(I)= -1d0 HV(J)= 1d0 HV(Ian+1)= H(J)-H(I)-hh(M) P= 1/S(M) Call dPTC (Ian, P) Continue Do 65 I=(Nxd+1),Ian Do 66 J=I,Ian K=J*(J-1)/2+I PTC(K)=PTC(K)+1.d0 Continue Continue Ls=Ian-1 Call DTDNG(Ls) Call DNDAO(Ian) Do 70 I=1,Ian Do 71 J=I,Ian L1=J*(J-1)/2+I PTC(L1)=PTC(L1)-1.d0/(Ng*1.d0)**2 Continue Continue Do 80 I1=1,Ian RLV(I1)=0.d0 Continue Do 82 I=1,Ian Do 83 J=1,Ian If(J.GT.I) Goto 84 K=(I-1)*I/2+J Cao Bá Hạ - 82 - |',14X,'|') Lớp Trắc địa B-K48 Khoa Trắc địa 84 85 83 82 60 711 33 Đồ án tèt nghiÖp Goto 85 K=J*(J-1)/2+I RLV(I)=RLV(I)-PTC(K)*PTC(Ms+J) Continue Continue Write (4,60) 'III DO LECH CAO DO CUA CAC DIEM DINH VI' Format (12x, A40/ * 12X, 40('=')// * 10x, 44('-')/ * 10x, '| S | TEN | DO lech |'/ * 10x, '| T | | |'/ * 10x, '| T | DIEM | dH(mm) Ghi chu |'/ * 10x, '| |') Write (4, 711) (K,Khd(K),RLV(K),K=(nxd+1),Ndi) Format (10x,'|', I3,' |', A7,'|', F9.2,' |',14X,'|') Write(4,33) Format(10x,44('-')) Write (4,7) 'IV.KET QUA TRI DO SAU BINH SAI' Format (//15x, A30/ * 15x, 30('=')/ * 15x, 'TRI DO CHENH CAO SAU BINH SAI'/ * 15x, ' -'/ * 67('=')/ *'| SO | TEN - CANH | GIA TRI CHENH | HIEU | TRI CHENH CAO | Mh 501 502 503 500 504 259 * |'/ *'| THU| - | CAO DO | CHINH | BINH SAI | * |'/ *'| TU | DAU – CUOI | (mm) | (mm) | (mm) |(mm *)|'/ * 67('-')) PVV= 0.d0 Do 500 m= 1, ncc Do 501 K= 1, (Ian+1) HV(K)=0.d0 continue I= Kd(M) J= Kc(M) HV(I)= -1d0 HV(J)= 1d0 QF=0.d0 Call Dhts(Ian,QF) HV(Ian+1)= H(J)-H(I)-hh(M) P= 1/S(M) V= 0.d0 Do 502 K=1, Ian V= V+HV(K)*RLV(K) Continue V= V+HV(Ian+1) PVV= PVV+ V*V*P PVV= abs(PVV) SSDV= dSQRT(PVV/(Ncc-Ndi+1)) Sa=SSDV*Dsqrt(ABS(QF)) hhbs= hh(m)+V Write(4, 503) m, Khd(I), Khd(J), hh(m), V, hhbs, sa Format ('|', I3, ' | ', A6,'- ',A6, '|' * F10.2, 5X, '|', F7.2, '| ', F9.2, 4X,'|', F4.2, '|') Continue Write(4,504) Format(67('-')) Write (4, 259) 'V.DO CAO DIEM SAU BINH SAI' Format (//15x, A28/ * 17x, 26('=')// * 57('=')/ *'| S | TEN | DO CAO | SAI SO | | '/ *'| T | | | GHI CHU | '/ *'| T | DIEM | (mm) | (mm) | '/ *' | ') Cao Bá Hạ - 83 - Lớp Trắc địa B-K48 Khoa Trắc địa Đồ ¸n tèt nghiÖp Do 400 I= 1, Ian K= I*(I-1)/2 + I QH= ptc(k) SSH= SSDV*dSqrt(ABS(QH)) HBS= H(I)+RLV(I) Write (4, 401) I, Khd(I), HBS, SSH Format ('|', I3, ' | ', A7, '|', F10.2, ' 401 * 400 402 505 | ', F7.2,' 12x, '|') |' Continue Write(4,402) Format(57('-')) Write (4,505) ssdv Format (//15x, 'DO CHINH XAC DO DAC'/ * 15x, ' -'// * 15x,'Sai so trung phuong so don vi:',F5.2, 'mm'// * ' GHI CHU : '/ * 10x, '- Nguoi tinh : Cao Ba Ha'/ * 10x, '- Lop trac dia BK48'/ * 10x, '- Truong dai hoc Mo Dia chat'/) Stop End Chương trình tính chuyển toạ ®é $debug $large c CHUONG TRINH TINH CHUYEN TU XY SANG BL c Ho va ten :Cao Ba Ha c Lop : Trac dia B-K48 Implicit Real *8(a-h, o-z) Character ff*50,Tenct*100, td(100)*100 Dimension X(1000),Y(1000) b1= 1.d0 Pi= 4.d0*dAtan(b1) Pi2= 2.d0*Pi Ro= 3600.d0*180.d0/Pi Write (*,'(a\)') ' Nhap ten file so lieu : ' Read (*,'(a)') ff Open (1, file=ff, status='old') Write (*,'(a\)') ' Nhap ten file ket qua : ' Read (*,'(a)') ff Open (4, file=ff, status='new') 50 Read (1,'(a)') Tenct Read(1,*) Nd,a,b2,dMo Read(1,50) (td(I),I=1,Nd) Format(10A7) Read(1,*) Id1,Iph1,G1 GL0=(Id1*3600.d0+Iph1*60.d0+G1)/Ro Read (1,*) I0do, I0ph, g0 dl0= (I0do*3600.d0 + I0ph*60.d0 +g0)/Ro Do I=1,Nd Read(1,*) I,X(I),Y(I) Continue Gh=0.0000000000000001d0 a22=a*a b22=b2*b2 e2=(a22-b22)/a22 ep2=(a22-b22)/b22 c=a/dsqrt(1-e2) Ao1=1.d0-3.d0*ep2/4.d0+ 45.d0*(ep2**2)/64.d0 Ao2=-175.d0*(ep2**3)/256.d0+ 11025.d0*(ep2**4)/16384.d0 Cao Bá Hạ - 84 - Lớp Trắc địa B-K48 Khoa Trắc địa Đồ án tốt nghiệp Ao=Ao1+Ao2 A2=Ao-1.d0 A4=15.d0*(ep2**2)/32.d0-175.d0*(ep2**3)/384.d0+3675.d0*(ep2**4) */8192.d0 A6=-35.d0*(ep2**3)/96.d0+735.d0*(ep2**4)/2048.d0 A8=315.d0*(ep2**4)/1024.d0 Write(4,6) tenct Format(///32x,'KET QUA TINH CHUYEN TOA DO PHANG'/ *36x,'TRONG CUNG HE QUY CHIEU'// *27x,A100/ *32x,' '/ *113('=')/ *'| | | TOA DO TRONG HE CU |' *' TOA DO TRONG HE MOI | TOA DO TRAC DIA |'/ *'| STT | TEN DIEM | -|' *' |'/ *'| | | X(m) | Y(m) |' *' X(m) | Y(m) | B | L |'/ *'| -| | -| -|' *' -| -| -| |') Do I=1,Nd Y(I)=Y(I)-500000.d0 Bo1=20.d0 t=dcos(Bo1) t1=dsin(Bo1) t2=t*t t4=t2**2 t6=t2**3 Bo=(X(I)/(dMo*c)-(A2+A4*t2+A6*t4+A8*t6)*t*t1)/Ao ep1=ABS(Bo-Bo1) If(ep1.LE.Gh) Goto Bo1=Bo Goto to=dtan(Bo) tN=a/dsqrt(1.d0-e2*dsin(Bo)*dsin(Bo)) tM=(1.d0-e2)*tN/(1.d0-e2*dsin(Bo)*dsin(Bo)) P=tN/tM y2=Y(I)*Y(I) y3=Y(I)*Y(I)*Y(I) y4=y2*y2 y5=y2*y3 y6=y2**3 y7=y3*y4 y8=y2**4 dM2=dMo**2 dM3=dMo**3 dM4=dMo**4 dM5=dMo**5 dM6=dMo**6 dM7=dMo**7 dM8=dMo**8 tN3=tN**3 tN5=tN**5 tN7=tN**7 to2=to**2 to4=to**4 to6=to**6 P2=P**2 P3=P**3 P4=P**4 Cao Bá Hạ - 85 - Lớp Trắc địa B-K48 Khoa Trắc địa Đồ án tốt nghiệp R1=(to*y2)/(2.d0*dM2*tM*tN) R2=(to*y4)/(24.d0*dM4*tM*tN3) R21=(-4.d0*P2+9.d0*P*(1.d0-to2)+12.d0*to2) R3=(to*y6)/(720.d0*dM6*tM*tN5) R31=8.d0*P4*(11.d0-24.d0*to2)-12.d0*P*P2*(21.d0-71.d0*to2) R32=15.d0*P2*(15.d0-98.d0*to2+15.d0*to4)+180.d0*P*(5.d0*to2*3.d0*to4)+360.d0*to4 R4=(to*y8)/(40320.d0*dM8*tM*tN7) R41=(1385.d0+3633.d0*to2+4095.d0*to4+1575.d0*to6) GB=Bo-R1+R2*R21-R3*(R31+R32)+R4*R41 SB=1.d0/dcos(Bo) q1=Y(I)/(dMo*tN) q2=y3*(P+2.d0*to2)/(6.d0*dM3*tN3) q3=y5/(120.d0*dM5*tN5) q31=(-4.d0*P3*(1.d0-6.d0*to2)+P2*(9.d0-68.d0*to2)+72.d0*P*to2 *+24.d0*to4) q4=y7/(5040.d0*dM7*tN7) q41=(61.d0+662.d0*to2+1320.d0*to4+720.d0*to6) dl=sb*(q1-q2+q3*q31-q4*q41) GL=GL0+dl e=dsqrt((a22-b22)/a22) dl1=Gl-dl0 dt=dtan(GB) fi=(1.d0-e**2*(dsin(GB))**2)/(1.d0-e**2) b0=1.d0+0.75d0*e**2+(45.d0/64.d0)*e**4+(175.d0/256.d0)*e**6 b222=0.75d0*e**2+(15.d0/16.d0)*e**4+(525.d0/512.d0)*e**6 b4=(15.d0/64.d0)*e**4+(105.d0/256.d0)*e**6 b6=(35.d0/512.d0)*e**6 X0=a*(1.d0-e**2)*(b0*GB-(b222/2.d0)*dsin(2.d0*GB)+(b4/4.d0)* *dsin(4.d0*GB)-(b6/6.d0)*dsin(6.d0*GB)) dN=a/dsqrt(1.d0-e**2*(dsin(GB))**2) dc=dcos(GB) ds=dsin(GB) 16 dx=dmo*(X0+dN*ds*dc*dl1**2/2.d0+dN*ds*dc**3*(4.d0*fi**2+fi*dt**2)*dl1**4/24.d0+dN*ds*dc**5.d0*(8*fi**4*(11.d0-24.d0*dt**2) *-28.d0*fi**3*(1.d0-6*dt**2)+fi**2*(1.d0-32.d0*dt**2)-fi*(2.d0* *dt**2)+dt**4)*dl1**6/720.d0+dN*ds*dc**7*(1385.d0-3111.d0*dt**2+ *543.d0*dt**4-dt**6)*dl1**8/40320.d0) dy=dmo*(dN*dl1*dc+dN*dc**3*(fi-dt**2)*dl1**3/6.d0+dN*dc**5* *(4.d0*fi**3*(1.d0-6.d0*dt**2)+fi**2*(1+8.d0*dt**2)-2.d0*fi* *dt**2+dt**4)*dl1**5/120.d0+dN*dc**7*(61.d0-479.d0*dt**2+179.d0* *dt**4-dt**6)*dl1**7/5040.d0) dy=dy+500000.d0 Y(I)=Y(I)+500000.d0 Call Drado(GB,Id,Iph,g,Ro) Call Drado(GL,Id2,Iph2,g2,Ro) Write (4,7) I,td(I),x(I),Y(I), dx, dy,Id,Iph,g,Id2,Iph2,g2 Format ( '|',I5,'| ',A7,' |',F15.3,'|',F15.3,'|',F15.3,'|',F15.3, *'|',I3 , I3 , F8.4 ,' |', I3 , I3 , F8.4,'|') Continue Write(4,8) Format (113('=')) Write(4,16) Format(///60x,' Nguoi thuc hien: Cao Ba Ha '/ * 60x,' Lop : Trac dia B - K48'/ * 60x,' TRUONG DAI HOC MO - DIA CHAT '/) Stop End Cao Bá Hạ - 86 - Lớp Trắc địa B-K48 Khoa Trắc địa Đồ án tốt nghiệp tài liệu tham khảo [1] GS.TSKH Hoàng Ngọc Hà - TS Trương Quang Hiếu (1999), Cơ sở toán học xử lý số liệu trắc địa, Nhà xuất giao thông vận tải [2] Phan Văn Hiến - Ngô Văn Hợi - Trần Khánh - Nguyễn Quang Phúc Nguyễn Quang Thắng - Phan Hồng Tiến - Trần Viết Tuấn (2001), Trắc địa công trình, Nhà xuất giao thông vận tải [3] Nguyễn Quang Phúc (2004), Đặc điểm công tác thiết kế lưới khống chế mặt phẳng trắc địa công trình, Tạp chí KHKT Mỏ - Địa chất, Trường Đại học Mỏ - Địa chất [4] Trần Khánh (1996), Nghiên cứu øng dơng b×nh sai tù lÜnh vùc xư lý số liệu trắc địa công trình, Luận án PTS khoa học kỹ thuật, Thư viện trường Đại học Mỏ-Địa chất, Hà Nội [5] Nguyễn Quang Phúc (2002), Tiêu chuẩn ổn định điểm khống chế sở quan trắc chuyển dịch ngang công trình, Tuyển tập báo cáo hội nghị khoa học lần thứ 15, Đại học Mỏ - Địa chất, Hà Nội [6] Nguyễn Quang Phúc (2001), Nghiên cứu phương pháp phân tích độ ổn định mốc chuẩn xử lý số liệu đo lún công trình, Luận văn thạc sỹ kỹ thuật,Thư viện trường Đại học Mỏ - Địa chất Hà Nội [7] Phan Văn Hiến - Vi Trường - Trương Quang Hiếu (1985), Lý thuyết sai số phương pháp số bình phương nhỏ nhất, Nhà xuất giao thông vận tải [8] Nguyễn Quang Phúc ( ), Bàn thêm vấn đề định vị lưới tự trắc địa công trình, Tạp chí KHKT Mỏ- Địa Chất số 19, Trường Đại Học Mỏ- Địa Chất [9] Đỗ Ngọc Đường - Đặng Nam Chinh (2007), Bài giảng công nghệ GPS, Trường Đại Học Mỏ - Địa Chất [10] Nguyễn Quang Phúc (2004), Nghiên cứu ảnh hưởng đại lượng đo đến hình dạng kích thước ellipse sai số, Tạp chí KHKT Mỏ - Địa Chất số 7, Trường Đại Học Mỏ - Địa Chất Cao Bá Hạ - 87 - Lớp Trắc địa B-K48 ... bình sai để xử lý mạng lưới thi công, có số phương pháp bình sai thông dụng sau: phương pháp bình sai lưới phụ thuộc, phương pháp bình sai lưới tự Phương pháp bình sai lưới phụ thuộc: Phương pháp. .. án tốt nghiệp Chương ứng dụng phương pháp bình sai lưới trắc địa tự để xử lý lưới thi công công trình 3.1 thuật toán xử lý số liệu lưới thi công Lưới thi công thành lập giai đoạn xây dựng công. .. quan vỊ l­íi khèng chế trắc địa công trình Chương 2: Khảo sát phương pháp bình sai lưới tự Chương 3: ứng dụng phương pháp bình sai lưới trắc địa tự để xử lý lưới thi công công trình Kết luận Hà