Khảo sát phương pháp bình sai lưới trắc địa tự do và ứng dụng trong xử lý số liệu quan trắc lún công trình

LUẬN VĂN Đề tài: “Khảo sát phương pháp bình sai lới trắc địa tự vµ ứng dụng xử lý số liệu quan trắc lún cơng trình” LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Khoa Trắc địa Đồ án tốt nghiƯp mơc lơc Trang Mơc lơc Mở đầu Ch­¬ng - quan trắc lún công trình 1.1 Những vấn đề chung quan trắc chuyển dịch biến dạng công trình 1.2 Quan trắc lún công tr×nh 1.3 Thực trạng công tác quan trắc công trình nước ta 20 Chương - khảo sát phương pháp bình sai lưới trắc địa tự 22 2.1 Mét sè khái niệm lưới trắc địa tự 22 2.2 Mô hình toán học phương pháp bình sai lưới trắc địa tự 23 2.3 Tính chất kết bình sai lưới tự 28 2.4 Vấn đề định vị hƯ thèng l­íi ®é cao ®o lón 30 Ch­¬ng – øng dơng phương pháp bình sai lưới tự để xử lý số liệu quan trắc lún công trình 32 3.1 Thuật to¸n 32 3.2 S¬ đồ khối quy trình xử lý lưới độ cao ®o lón 35 3.3 Lập trình bình sai lưới quan trắc độ lún 38 3.4 Ch­¬ng trình nguồn tệp liệu 41 3.5 Sư dơng chương trình 49 3.6 TÝnh to¸n thùc nghiÖm 51 KÕt luËn 58 Tài liệu tham khảo 59 Phô lôc 60 Phô lôc 63 Phô lôc 66 Lương Anh Tuấn -1- Lớp Trắc ®Þa B-K48 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Khoa Trắc địa Đồ án tốt nghiệp Mở đầu Đối với công tác quan trắc lún công trình, tính đắn trình lún công trình phụ thuộc vào độ xác quan trắc, mà chịu ảnh hưởng lớn phương pháp xử lý số liệu Tuy nhiên, phương pháp xử lý số liệu quan trắc lún công trình thực tế chưa trọng mức Vì vậy, việc nghiên cứu đề biện pháp quy trình xử lý số liệu quan trắc lún công trình cách hợp lý, phù hợp với đặc điểm chất lưới quan trắc biến dạng vấn đề thời cần thiết Lưới trắc địa công trình nói chung lưới quan chuyển dịch công trình nói riêng xây dựng theo quy trình tiêu kỹ thuật riêng nhằm giải nhiệm vụ đa dạng chuyên nghành Do không giống lưới đo vẽ đồ, mà tính đặc thù cao, đòi hỏi cao độ xác, hệ thống điểm gốc khởi tính không ổn định Với đặc thù lưới trắc địa công trình nêu đòi hỏi phải có kỹ thuật xử lý số liệu riêng phù hợp với đặc điểm chất lưới Nhận thấy tầm quan trọng công tác xử lý số liệu quan trắc lún công trình, thời gian làm đồ án tốt nghiệp khoá học, em đà chọn nghiên cứu đề tài với nội dung: Khảo sát phương pháp bình sai lưới trắc địa tự ứng dụng xử lý số liệu quan trắc lún công trình Nội dung đồ án em trình bày chương: Chương 1: Quan trắc lún công trình Chương 2: Khảo sát phương pháp bình sai lưới trắc địa tự Chương 3: ứng dụng phương pháp bình sai lưới trắc địa tự để xử lý số liệu quan trắc lún công trình Em xin chân thành cảm ơn giúp đỡ bảo tận tình thầy giáo TS NGUYễN QUANG PHúC suốt trình em làm đồ án Do thời gian chuyên môn có hạn nên đồ án không tránh khỏi khiếm khuyết Em mong nhận đóng góp ý kiến thầy cô giáo bạn đồng nghiệp đề đồ án hoàn thiện Em xin chân thành cảm ơn ! Lương Anh Tuấn -2- Lớp Trắc địa B-K48 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Khoa Trắc địa Đồ án tốt nghiệp Chương QUAN TRắC LúN CÔNG TRìNH 1.1 NHữNG VấN Đề chung Về quan trắc chuyển dịch biến dạng công trình 1.1.1 Hiện tượng chuyển dịch biến dạng công trình a Hiện tượng chuyển dịch Là thay đổi vị trí công trình không gian theo thời gian so với vị trí ban đầu Có thể chia chuyển dịch công trình thành hai loại: - Chuyển dịch thẳng đứng: thay đổi vị trí công trình theo phương dây dọi Chuyển dịch theo hướng xuống gọi lún Chuyển dịch theo hướng lên gọi trồi - Chuyển dịch ngang: thay vị trí công trình mặt phẳng nằm ngang Chuyển dịch ngang theo hướng theo hướng xác định (hướng áp lực lớn nhất) b Hiện tượng biến dạng Là thay đổi hình dạng kích thước công trình không gian theo thời gian Biến dạng hậu tất yếu chuyển dịch không công trình biểu thường gặp sự: cong, vênh, vặn xoắn, vết rạn nứt 1.1.2 Nguyên nhân gây nên chuyển dịch biến dạng công trình Có nhiều nguyên nhân gây tượng chuyển dịch biến dạng công trình, quy nạp lại có hai nhóm nguyên nhân Cụ thể: a Nhóm nguyên nhân liên quan đến điều kiện tự nhiên Nhóm nguyên nhân gây : Tính chất lý lớp đất đá móng công trình, ảnh hưởng yếu tố khí tượng (như nhiệt độ, ®é Èm, h­íng chiÕu s¸ng ), sù thay ®ỉi chÕ độ nước mặt, nước ngầm vận động nội sinh lòng trái đất gây nên chuyển dịch biến dạng công trình (tuy nhiên mức độ chuyển dịch nguyên nhân gây thường bé) Lương Anh Tuấn -3- Lớp Trắc địa B-K48 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Khoa Tr¾c địa Đồ án tốt nghiệp b Nhóm nguyên nhân có liên quan đến trình xây dựng vận hành công trình Trong trình xây dựng vận hành công trình gia tăng tải trọng công trình, sai sót trình khảo sát địa chất công trình, việc khai thác nước ngầm gây nên tượng sụt lún lòng đất việc xây dựng công trình ngầm, công trình xây chen .đà gây nên chuyển dịch biến dạng công trình 1.1.3 Công tác quan trắc chuyển dịch biến dạng công trình a Mục đích quan trắc Công tác quan trắc chuyển dịch biến dạng công trình cần tiến hành theo phương án kỹ thuật nhằm đạt mục đích sau: - Thứ xác định giá trị độ lún, độ chuyển dịch tuyệt đối tương đối nhà công trình so với giá trị tính toán theo thiết kế chúng Từ tìm nguyên nhân gây lún, chuyển dịch mức độ nguy hiểm chúng trình làm việc bình thường nhà công trình sở đưa giải pháp nhù hợp nhằm phòng ngưà cố xảy ra; - Thứ hai xác định thông số đặc trưng cần thiết độ ổn định công trình, làm xác thêm số liệu đặc trưng cho tính chất lý đất; Dùng làm số liệu kiểm tra phương pháp tính toán, xác định giá trị độ lún, độ chuyển dịch giới hạn cho phép loại đất công trình khác b Nguyên tắc thực công tác quan trắc Công tác quan trắc chuyển dịch biến dạng công trình tiến hành theo nguyên tắc sau: - Việc quan trắc chuyển dịch biến dạng phải thực theo nhiều thời điểm, thời điểm gọi chu kỳ Chu kỳ đầu gọi chu kỳ - Chuyển dịch biến dạng công trình so sánh tương đối tượng khác xem ổn định - Chuyển dịch biến dạng công trình thường có trị số nhỏ phải có Lương Anh Tuấn -4- Lớp Trắc địa B-K48 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Khoa Trắc địa Đồ án tốt nghiệp phương pháp phương tiễn có độ xác cao - Cần phải có kỹ thuật xử lý riêng phù hợp với đặc điểm chất mạng lưới quan trắc chuyển dịch biến dạng công trình c Yêu cầu độ xác quan trắc Yêu cầu độ xác quan trắc chuyển dịch độ xác cần thiết xác định chuyển dịch công trình, tiêu định lượng đại lượng phụ thuộc chủ yếu vào tính chất lý đất đá móng, đặc điểm kết cấu vận hành công trình Yêu cầu độ xác xác định theo giá trị chuyển dịch dự báo (cho thiết kế) xác định theo tiêu chuẩn xây dựng, vận hành công trình (quy định tiêu chuẩn ngành) - Nếu theo độ chuyển dịch dự báo (cho thiết kế xác định theo số chu kỳ đà quan trắc), yêu cầu độ xác quan trắc xác định theo công thức: mQ Q Với mQ yêu cầu độ xác quan trắc thời điểm t Q giá trị chuyển dịch dự báo chu kỳ quan trắc hệ số đặc trưng cho độ tin cậy kết quan trắc Thường chọn lúc có mQ 0.17 Q - Nếu yêu cầu độ xác quan trắc xác định dựa vào quy mô công trình tính chất đất móng công trình yêu cầu độ xác quan trắc quy định theo bảng 1.1 Bảng 1.1 Độ xác quan trắc Công trình xây dựng đá gốc nửa đá gốc Công trình sét chịu lực Các loại đập đất, đá chịu lực cao Các công trình xây dựng trượt Các loại công trình đất đắp Độ xác quan trắc(mm) 1.0 3.0 5.0 10.0 15.0 Lương Anh Tuấn Lớp Trắc địa B-K48 Loại công trình móng -5- LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Khoa Trắc địa Đồ án tốt nghiệp d Chu kỳ quan trắc Nhìn chung chu kỳ quan trắc chuyển dịch biến dạng công trình quy định phụ thuộc vào yếu tố sau: - Tùy thuộc vào loại công trình tính chất đất đá móng công trình - Tùy thuộc vào giai đoạn xây dựng vận hành công trình Đối với quan trắc chuyển dịch thẳng đứng công trình Việc xác định thời gian đo (chu kỳ đo) chiếm vai trò quan trọng Theo kinh nghiệm quan trắc công trình ngời ta chia làm giai đoạn: - Quan trắc lún giai đoạn thi công; - Quan trắc lún công trình đa vào sử dụng; Giai đoạn thi công, quan trắc lún thờng đợc xác định theo tiến độ thi công mức độ phức tạp công trình Để dễ dàng cho việc theo dõi, ngời ta đo theo tải trọng hoàn thành trình xây dựng cụ thể là: - Công trình hoàn thành xong phần móng - Công trình đạt tới 20% tải trọng - Công trình đạt tới 50% tải trọng - Công trình đạt tới 75% tải trọng - Công trình đạt tới 100% tải trọng Đối với công trình phức tạp, việc theo dõi chuyển dịch biến dạng móng (khi hoàn thành xây xong phần móng) đạt 10% tải trọng cần phải quan trắc lần Tại lần quan trắc, kết so sánh với lần đo trớc gần sau xem xét hiệu chênh lệch cao hai lần đo kề h (độ lún) sở để định việc tăng dầy lần đo hay tiến hành đo theo tiến độ đà ấn định từ đầu - giai đoạn thứ hai công trình đ đa vào sử dụng Việc phân định số lần đo phụ thuộc hoàn toàn vào yêu cầu độ xác đo lún công trình nh đ trình bày Nếu sai số cho phép đo cấp xác nhỏ chu kỳ (thời gian) cách lớn ngợc lại sai số cho phép đo độ xác lớn chu kỳ đo cách Khi công trình Lương Anh Tuấn -6- Lớp Trắc địa B-K48 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Khoa Trắc địa Đồ án tốt nghiệp có dấu hiệu biến dạng lớn chu kỳ đo với số yêu cầu đặc biệt ngời t vấn thiết kế quy định Thời kỳ công trình vào ổn định (tốc độ chuyển dịch công trình đạt từ 1mm /năm 2mm/năm), thời kỳ chu kỳ quan trắc tháng hoăc năm năm Đối với quan trắc chuyển dịch ngang công trình Thời gian thực chu kỳ quan trắc chuyển dịch tiến hành dựa vào yếu tố : - Loại nhà công trình; - Loại đất xây dựng nhà công trình; - Đặc điểm áp lực ngang; - Mức độ chuyển dịch ngang; - Tiến độ thi công xây dựng công trình Chu kỳ quan trắc thực sau xây dựng móng công trình trước có áp lực ngang tác động đến công trình Các chu kỳ thực tuỳ thuộc vào mức tăng giảm áp lực ngang tác động vào công trình quan trắc tháng lần thời gian xây dựng công trình Trong thời gian sử dụng công trình, số lượng chu kỳ quan trắc tiÕn hµnh tõ  chu kú mét năm, vào thời điểm mà điều kiện ngoại cảnh khác biệt Ngoài cần phải quan trắc bổ sung công trình có độ chuyển dịch ngang lớn, quan trắc bổ sung để tìm nguyên nhân gây nên cố công trình 1.2 quan trắc lún công trình 1.2.1 Các phương pháp quan trắc lún công trình a Đo cao hình học Phương pháp dựa nguyên lý sử dụng tia ngắm nằm ngang xác định chênh cao hai điểm (Hình 1.1) Lương Anh Tuấn -7- Lớp Trắc địa B-K48 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Khoa Trắc địa Đồ án tốt nghiệp Hình 1.1 Trạm đo cao hình học Nếu máy thủy chuẩn đặt khoảng A, B, ký hiệu (a), (b) số đọc tương ứng mia sau (đặt A) mia trước (đặt B), chênh cao hai điểm A, B tính theo công thức: hAB = (a) (b) Việc quan trắc để xác định độ lún công trình phải đợc tiến hành theo quy định đo cao hình học xác đặc biệt hay gọi đo cao hình học tia ngắm ngắn Những tiêu kỹ thuật đo cao hình học tia ngắm ngắn quy định bảng 1.3 Bảng 1.3 Chỉ tiêu kỹ thuật thủy chuẩn hình học tia ngắm ngắn Hạng thủy chuẩn Chỉ tiêu Chiều dài tia ngắm (m) Chiều cao tia ngắm (m) Chênh lệch khoảng ngắm (m): - Trên trạm - Trên toàn tuyến Sai số khÐp cho phÐp (mm) I ≤ 25 ≥ 0,8 II ≤ 25 ≥ 0,5 III ≤ 40 ≥ 0,3 0,4 2,0 1,0 4,0 2,0 5,0 ≤  0,3 n ≤  1,0 n ≤  2,0 n b §o cao thuỷ tĩnh Phương pháp đo cao thủy tĩnh dựa nguyên lý bình thông nhau: Bề mặt chất lỏng bình thông có vị trí nằm ngang (vuông góc phương dây dọi) có độ cao, không phụ thuộc vào hình dạng mặt cắt khối lượng chất lỏng bình Lương Anh Tuấn -8- Lớp Trắc địa B-K48 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Khoa Trắc địa Đồ án tốt nghiÖp File sè liÖu : 4 10 11 12 13 14 15 16 17 20 20 18 19 20 10 23 20 22 21 24 25 28 29 27 26 28 26 26 29 28 10 11 12 13 14 15 16 17 18 19 20 25 10 11 12 13 14 15 16 17 1 18 19 13 10 23 22 21 24 25 10 28 29 29 27 27 27 N1 N2 N3 N4 N5 N6 N7 N8 N9 N10 N11 N12 N13 N14 N15 N16 N17 N18 N19 N20 38 0.95 29.07 -70.18 17.62 63.59 -87.44 95.13 0.06 -100.65 -28.69 29.08 45.63 -44.86 -31.49 43.9 -78.92 -2.3 120.32 91.49 -21.77 193.36 -131.81 9.87 -56.43 84.45 63.75 61.74 -114.7 178.66 -146.43 537.41 537.08 206.51 -186.26 -39.46 -226.03 22.47 248.81 209.58 L­¬ng Anh TuÊn 1 1 1 1 1 1 1 2 1 1 3 3 - 54 - Lớp Trắc địa B-K48 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Khoa Trắc địa 21 22 23 24 25 26 27 28 29 N21 N22 N23 N24 N25 M1 TC M2 M3 §å ¸n tèt nghiƯp 8618.26 8641.7 8429.74 8389.9 Sư dơng chương trình xử lý số liệu viết để tính toán bình sai, kết sau: KET QUA TINH TOAN BINH SAI LUOI DO CAO DO LUN ************************************************* CONG TRINH: DON NGUYEN II - THUOC LO - KDTM DINH CONG CHU KY: 11 - NGAY DO: 07/09/2006 CAC CHI TIEU KY THUAT CUA LUOI : ================================ 1.So luong diem quan trac: 25 2.So luong diem DINH VI : 3.So luong chenh cao do: 38 PHAN TICH DO ON DINH CUA CAC DIEM CO SO: ======================================== Phan tich lan : =================== BANG: ============================================================== | TEN | DO CAO | DO LECH | KET QUA | | DIEM | (mm) | CAO DO | PHAN TICH | ============================================================== | M1 | 8618.26 | -.88 | On dinh | | TC | 8641.70 | -1.56 | Khong on dinh | | M2 | 8429.74 | +1.05 | Khong on dinh | | M3 | 8389.90 | +1.39 | Khong on dinh | ============================================================== Loai diem TC khoi danh sach cac diem DINH VI Phan tich lan : =================== BANG: ============================================================== | TEN | DO CAO | DO LECH | KET QUA | | DIEM | (mm) | CAO DO | PHAN TICH | ============================================================== | M1 | 8618.26 | -1.40 | Khong on dinh | | M2 | 8429.74 | +.53 | On dinh | | M3 | 8389.90 | +.87 | On dinh | ============================================================== Loai diem L­¬ng Anh TuÊn M1 khoi danh sach cac diem DINH VI - 55 - Lớp Trắc địa B-K48 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Khoa Trắc địa Đồ án tốt nghiÖp Phan tich lan : =================== BANG: ============================================================== | TEN | DO CAO | DO LECH | KET QUA | | DIEM | (mm) | CAO DO | PHAN TICH | ============================================================== | M2 | 8429.74 | -.17 | On dinh | | M3 | 8389.90 | +.17 | On dinh | ============================================================== Lay diem o bang de dinh vi luoi! ==================================== BINH SAI VA VA DANH GIA DO CHINH XAC: ===================================== Tri sau binh sai va chinh xac cua tri : ================================================ BANG: ============================================================== | T | TUYEN DO |CH CAO DO| Vh | SO |CH CAO SAU| Mh | | T | DAU-CUOI | (mm) | (mm) |TRAM | BINH SAI | (mm) | ============================================================== | | N1 N2 | +29.07 | -.06 | | +29.01 | 16 | | | N2 N3 | -70.18 | -.06 | | -70.24 | 16 | | | N3 N4 | +17.62 | -.06 | | +17.56 | 16 | | | N4 N5 | +63.59 | -.06 | | +63.53 | 16 | | | N5 N6 | -87.44 | -.07 | | -87.51 | 16 | | | N6 N7 | +95.13 | +.01 | | +95.14 | 16 | | | N7 N8 | +.06 | +.01 | | +.07 | 16 | | | N8 N9 | -100.65 | +.01 | | -100.64 | 16 | | | N9 N10 | -28.69 | +.23 | | -28.46 | 14 | | 10 | N10 N11 | +29.08 | +.02 | | +29.10 | 16 | | 11 | N11 N12 | +45.63 | +.02 | | +45.65 | 16 | | 12 | N12 N13 | -44.86 | +.02 | | -44.84 | 16 | | 13 | N13 N14 | -31.49 | +.02 | | -31.47 | 17 | | 14 | N14 N15 | +43.90 | +.02 | | +43.92 | 17 | | 15 | N15 N16 | -78.92 | +.02 | | -78.90 | 17 | | 16 | N16 N17 | -2.30 | +.02 | | -2.28 | 17 | | 17 | N17 N1 | +120.32 | +.02 | | +120.34 | 17 | | 18 | N20 N1 | +91.49 | -.12 | | +91.37 | 21 | | 19 | N20 N18 | -21.77 | 00 | | -21.77 | 17 | | 20 | N18 N19 | +193.36 | -.01 | | +193.35 | 22 | | 21 | N19 N13 | -131.81 | -.01 | | -131.82 | 22 | | 22 | N20 N10 | +9.87 | -.04 | | +9.83 | 15 | | 23 | N10 N23 | -56.43 | +.22 | | -56.21 | 14 | | 24 | N23 N9 | +84.45 | +.22 | | +84.67 | 14 | | 25 | N20 N22 | +63.75 | +.09 | | +63.84 | 17 | | 26 | N22 N21 | +61.74 | +.09 | | +61.83 | 17 | | 27 | N21 N24 | -114.70 | +.17 | | -114.53 | 23 | | 28 | N24 N25 | +178.66 | +.09 | | +178.75 | 17 | | 29 | N25 N6 | -146.43 | +.26 | | -146.17 | 26 | | 30 | M2 N5 | +537.41 | -.02 | | +537.39 | 20 | | 31 | M3 N1 | +537.08 | -.04 | | +537.04 | 16 | | 32 | TC N10 | +206.51 | +.15 | | +206.66 | 23 | ============================================================== L­¬ng Anh Tuấn - 56 - Lớp Trắc địa B-K48 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Khoa Trắc địa Đồ ¸n tèt nghiÖp BANG: (Tiep) ============================================================== | T | TUYEN DO |CH CAO DO| Vh | SO |CH CAO SAU| Mh | | T | DAU-CUOI | (mm) | (mm) |TRAM | BINH SAI | (mm) | ============================================================== | 33 | M1 M2 | -186.26 | -.33 | | -186.59 | 23 | | 34 | M2 M3 | -39.46 | -.04 | | -39.50 | 18 | | 35 | M1 M3 | -226.03 | -.06 | | -226.09 | 24 | | 36 | M1 TC | +22.47 | +.29 | | +22.76 | 23 | | 37 | M3 TC | +248.81 | +.04 | | +248.85 | 20 | | 38 | M2 TC | +209.58 | -.23 | | +209.35 | 21 | ============================================================== Do cao binh sai va sai so cua chung : ===================================== BANG: ============================================================== | T | TEN DIEM | DO CAO | Mx | GHI | | T | | (mm) | (mm) | CHU | ============================================================== | | N1 | 8927.11 | 17 | | | | N2 | 8956.12 | 21 | | | | N3 | 8885.88 | 23 | | | | N4 | 8903.43 | 22 | | | | N5 | 8966.96 | 20 | | | | N6 | 8879.45 | 23 | | | | N7 | 8974.59 | 25 | | | | N8 | 8974.67 | 26 | | | | N9 | 8874.03 | 24 | | | 10 | N10 | 8845.58 | 22 | | | 11 | N11 | 8874.68 | 26 | | | 12 | N12 | 8920.34 | 28 | | | 13 | N13 | 8875.50 | 27 | | | 14 | N14 | 8844.03 | 28 | | | 15 | N15 | 8887.95 | 28 | | | 16 | N16 | 8809.05 | 26 | | | 17 | N17 | 8806.77 | 23 | | | 18 | N18 | 8813.97 | 27 | | | 19 | N19 | 9007.32 | 30 | | | 20 | N20 | 8835.75 | 23 | | | 21 | N21 | 8961.41 | 30 | | | 22 | N22 | 8899.58 | 27 | | | 23 | N23 | 8789.36 | 25 | | | 24 | N24 | 8846.88 | 32 | | | 25 | N25 | 9025.62 | 31 | | | 26 | TC | 8638.92 | 18 | | | 27 | M1 | 8616.16 | 22 | | | 28 | M2 | 8429.74 | 09 | | | 29 | M3 | 8389.90 | 09 | | ============================================================== Sai so trung phuong so don vi : 18mm Chuong trinh duoc viet bang ngon ngu lap trinh VISUALBASIC6.0 L­¬ng Anh TuÊn - 57 - Lớp Trắc địa B-K48 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Khoa Trắc địa Đồ án tốt nghiệp Kết luận Từ kết nghiên cứu lý thuyết tính toán ứng dụng thực tiễn đà trình bày nội dung đồ án đà cho phép rút kết luận sau đây: Khi xử lý số liệu đo lún công trình cần kết hợp phân tích độ ổn định mốc độ cao sở Chỉ điểm độ cao sở ổn định tham gia vào việc định vị mạng lưới quan trắc Như độ lún thực tế công trình phản ánh cách khách quan 2.Việc sử dụng phương pháp bình sai tù ®Ĩ xư lý hƯ thèng l­íi ®é cao đo lún giải pháp toàn diện, phù hợp với đặc điểm chất dạng lưới Các kết nghiên cứu lý thuyết đà cụ thể hoá việc xây dựng chương trình phần mềm để giải nhiệm vụ đặt máy tính So sánh kết bình sai chương trình viết kết thầy cô môn hoàn toàn giống Điều khẳng định thuật toán, quy trình xử lý chương trình bình sai mà đà trình bày đồ án đắn, đáp ứng đầy đủ yêu cầu kü tht vỊ xư lý sè liƯu quan tr¾c lón công trình Xử lý số liệu đo biến dạng công trình nói chung đo lún công trình nói riêng có vai trò quan trọng, nên vấn đề nhiều người quan tâm Nhận thức tầm quan trọng vấn đề này, thời gian làm đồ án tốt nghiệp khoá học em đà cố gắng nhiều để hoàn thiện quy trình công nghệ công tác trắc địa xử lý số liệu quan trắc lún công trình Hy vọng nhận góp ý thầy cô bạn đồng nghiêp để chương trình hoàn thiện Một lần em xin chân thành cảm ơn thầy giáo TS Nguyễn Quang Phúc, thầy cô giáo khoa trắc địa bạn đồng nghiệp ! Lương Anh Tuấn - 58 - Lớp Trắc địa B-K48 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Khoa Trắc địa Đồ án tốt nghiệp Tài liệu tham khảo [1] Phan Văn Hiến (1997), Quan trắc chuyển dịch biến dạng công trình, Trường đại học Mỏ-Địa chất, Hà Nội [2] 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, Trường Đại học Mỏ-Địa chất, Hà Nội [3] Nguyễn Quang Phúc (2001), Nghiên cứu phương pháp xử lý số liệu đo lún công trình xây dựng , Tuyển tập công trình khoa học- Tập 32, Trường đại học Mỏ-Địa chất, Hà Nội [4] Nguyễn Quang Phúc (2001), Tiêu chuẩn ổn định điểm độ cao sở đo lún công trình , Tuyển tập công trình khoa học- Tập 33, Trường đại học Mỏ-Địa chất, Hà Nội [5] Nguyễn Quang Phúc (2007), Quan trắc phân tích biến dạng công trình, Bài giảng dùng cho học viên cao học, Trường đại học Mỏ-Địa chất, Hà Nội [6] Trần Khánh (1996), Nghiên cứu ứng dụng bình sai tự xử lý số liệu trắc địa công trình, Luận án PTS khoa học kỹ thuật, Trường đại học Mỏ-Địa chất, Hà Nội [7] Tiêu chuẩn xây dựng Việt Nam 203 - 1997 Bộ Xây Dựng Lương Anh Tuấn - 59 - Lớp Trắc địa B-K48 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Khoa Trắc địa Mà nguồn modul Dự báo lún Đồ án tốt nghiệp Phô lôc Private Sub mnuPhuthuoc_Click() Dim tencongtrinh As String tencongtrinh = txtTencongtrinh Dim Chukydo As String Chukydo = txtChukydo Dim Ngaydo As String Ngaydo = txtngaydo Dim Dg As Long, nDqt As Long, nh As Long Dg = VSGrid.TextMatrix(0, 0) nDqt = VSGrid.TextMatrix(0, 1) nh = VSGrid.TextMatrix(0, 2) Dim id() As Long Dim ic() As Long Dim h() As Double Dim nt() As Long Dim i, j As Long ReDim id(1 To nh), ic(1 To nh), h(1 To nh), nt(1 To nh) For i = To nh id(i) = VSGrid.TextMatrix(i, 0) ic(i) = VSGrid.TextMatrix(i, 1) h(i) = VSGrid.TextMatrix(i, 2) nt(i) = VSGrid.TextMatrix(i, 3) Next Dim Tendiem() As String ReDim Tendiem(1 To Dg + nDqt) For i = To (Dg + nDqt) Tendiem(i) = VSGrid.TextMatrix(nh + i, 1) Next Dim hgd() As Double ReDim hgd(1 To Dg + nDqt) Dim A As Double A = -99999999 For i = To nDqt hgd(i) = A Next j = For i = (nDqt + 1) To (nDqt + Dg) j = j + hgd(i) = VSGrid.TextMatrix(nh + nDqt + j, 2) Next Dim dem As Long dem = 1: For i = To nh If (hgd(id(i)) A) And (hgd(ic(i)) = A) Then hgd(ic(i)) = hgd(id(i)) + h(i) dem = dem + End If If (hgd(id(i)) = A) And (hgd(ic(i)) A) Then hgd(id(i)) = hgd(ic(i)) - h(i) dem = dem + End If Lương Anh Tuấn - 60 - Lớp Trắc địa B-K48 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Khoa Trắc địa Đồ án tốt nghiệp Next If dem < nDqt Then GoTo End If Dim mtA() As Double ReDim mtA(1 To nh, To nDqt + Dg) For i = To nh For j = To nDqt + Dg mtA(i, j) = Next mtA(i, id(i)) = -1 mtA(i, ic(i)) = Next Dim mtB() As Double ReDim mtB(1 To nh, To nDqt) For i = To nh For j = To nDqt mtB(i, j) = mtA(i, j) Next Next Dim l() As Double ReDim l(1 To nh, To 1) For i = To nh l(i, 1) = (hgd(ic(i)) - hgd(id(i))) - h(i) Next Dim p() As Double ReDim p(1 To nh, To nh) For i = To nh p(i, i) = / nt(i) Next Dim N() As Double ReDim N(1 To nDqt, To nDqt) N = Nhan_MaTran(chuyenvi_matran(mtB(), nh, nDqt), _ Nhan_MaTran(p(), mtB(), nh, nh, nDqt), nDqt, nh, nDqt) Dim M() As Double ReDim M(1 To nh, To 1) M = Nhan_MaTran(chuyenvi_matran(mtB(), nh, nDqt), _ Nhan_MaTran(p(), l(), nh, nh, 1), nDqt, nh, 1) Dim Nnd() As Double ReDim Nnd(1 To nDqt, To nDqt) Nnd = nghichdao_matran(N(), nDqt) Dim x() As Double ReDim x(1 To nDqt, To 1) x = DoiDau_MaTran(Nhan_MaTran(Nnd(), M(), nDqt, nDqt, 1), nDqt, 1) Dim Hbs() As Double ReDim Hbs(1 To nDqt, To 1) Hbs = cong_matran(ho(), x(), nDqt, 1) Dim V() As Double ReDim V(1 To nh, To 1) V = cong_matran(Nhan_MaTran(mtA(), x(), nh, nDqt, 1), l(), nh, 1) Dim hhbs() As Double ReDim hhbs(1 To nh, To 1) hhbs = cong_matran(hh(), V(), nh, 1) Dim vpv() As Double Dim mo() As Double ReDim vpv(1 To 1, To 1) Lương Anh Tuấn - 61 - Lớp Trắc ®Þa B-K48 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Khoa Trắc địa Đồ án tốt nghiệp ReDim mo(1 To 1, To 1) vpv = Nhan_MaTran(Nhan_MaTran(chuyenvi_matran(V(), nh, 1), p(), 1, nh, nh), V(), 1, nh, 1) For i = To For j = To mo(i, j) = Sqr(vpv(i, j) / (nh - nDqt)) Next Next Dim Rxx() As Double ReDim Rxx(1 To nDqt, To 1) For i = To nDqt Rxx(i, 1) = Sqr(Nnd(i, i)) Next Dim mx() As Double ReDim mx(1 To nDqt, To 1) mx = Nhan_MaTran(Rxx(), mo(), nDqt, 1, 1) Dim fRf() As Double ReDim fRf(1 To nh, To nh) fRf = Nhan_MaTran(Nhan_MaTran(mtA(), Nnd(), nh, nDqt, nDqt), _ chuyenvi_matran(mtA(), nh, nDqt), nh, nDqt, nh) Dim cfRf() As Double ReDim cfRf(1 To nh, To 1) For i = To nh cfRf(i, 1) = Sqr(fRf(i, i)) Next Dim mf() As Double ReDim mf(1 To nh, To 1) mf = Nhan_MaTran(cfRf(), mo(), nh, 1, 1) With CommonDialog1 DialogTitle = "In file nao" CancelError = True Filter = "Text file *.txt|*.txt|All Files | *.*" ShowSave End With str = CommonDialog1.FileName Dim FileNum2 As Long FileNum2 = FreeFile Open str For Output As FileNum2 Print #FileNum2, String(5, " ") & "Chuong trinh duoc viet bang ngon ngu lap trinh VISUALBASIC 6.0" Close FileNum2 RichTextBox1.LoadFile CommonDialog1.FileName End Sub Lương Anh Tuấn - 62 - Lớp Trắc địa B-K48 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Khoa Trắc địa Đồ án tốt nghiệp Phụ lục Mà nguồn modul Tính thông số chuyển dịch Private Sub mnuCacthongsochuyendich_Click() Dim tencongtrinh As String tencongtrinh = txtTencongtrinh Dim Chukydo As String Chukydo = txtChukydo Dim Ngaydo As String Ngaydo = txtngaydo Dim sodiemquantrac As Long Dim sochukyquantrac As Long Dim i, j As Long i = Do While VSGrid.TextMatrix(1, i) "" i = i + Loop sochukyquantrac = i - i = Do While VSGrid.TextMatrix(i, 1) "" i = i + Loop sodiemquantrac = i - Dim khoangthoigiantinhtuCkdau() As Double ReDim khoangthoigiantinhtuCkdau(1 To 1, To sochukyquantrac) Dim khoangthoigiantinh2Ck() As Double ReDim khoangthoigiantinh2Ck(1 To 1, To sochukyquantrac) For i = To sochukyquantrac khoangthoigiantinhtuCkdau(1, i) = VSGrid.TextMatrix(0, i) Next For i = To sochukyquantrac - khoangthoigiantinh2Ck(1, i) = khoangthoigiantinhtuCkdau(1, i + 1) - khoangthoigiantinhtuCkdau(1, i) Next Dim Tendiem() As String ReDim Tendiem(1 To sodiemquantrac) For i = To sodiemquantrac Tendiem(i) = VSGrid.TextMatrix(i, 0) Next Dim docao() As Double ReDim docao(1 To sodiemquantrac, To sochukyquantrac) For i = To sodiemquantrac For j = To sochukyquantrac docao(i, j) = VSGrid.TextMatrix(i, j) Next Next Dim DoLun() As Double ReDim DoLun(1 To sodiemquantrac, To 1) Dim Dlgiua2chukylientiep() As Double ReDim Dlgiua2chukylientiep(1 To sodiemquantrac, To sochukyquantrac - 1) Dim Dlsovoickdautien() As Double ReDim Dlsovoickdautien(1 To sodiemquantrac, To sochukyquantrac - 1) For i = To sodiemquantrac For j = To sochukyquantrac - L­¬ng Anh Tuấn - 63 - Lớp Trắc địa B-K48 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Khoa Trắc địa Đồ ¸n tèt nghiÖp DoLun(i, 1) = docao(i, j + 1) - docao(i, j) Dlgiua2chukylientiep(i, j) = DoLun(i, 1) Dlsovoickdautien(i, j) = docao(i, j + 1) - docao(i, 1) Next Next Dim TongDocaotungck() As Double ReDim TongDocaotungck(1 To 1, To sochukyquantrac) Dim TG() As Double ReDim TG(1 To sochukyquantrac, To sodiemquantrac) TG() = chuyenvi_matran(docao(), sodiemquantrac, sochukyquantrac) Dim A, b, c As Double For i = To sochukyquantrac b = For j = To sodiemquantrac c = b A = c + TG(i, j) b = A Next TongDocaotungck(1, i) = b NextDim TongDltheo2cklientiep() As Double ReDim TongDltheo2cklientiep(1 To 1, To sochukyquantrac - 1) Dim TongDltheoSovoick0() As Double ReDim TongDltheoSovoick0(1 To 1, To sochukyquantrac - 1) For i = To sochukyquantrac - TongDltheo2cklientiep(1, i) = TongDocaotungck(1, i + 1) TongDocaotungck(1, i) TongDltheoSovoick0(1, i) = TongDocaotungck(1, i + 1) TongDocaotungck(1, 1) Next Dim DlTrungbinhtheo2cklientiep() As Double ReDim DlTrungbinhtheo2cklientiep(1 To 1, To sochukyquantrac - 1) Dim DlTrungbinhtheoSovoick0() As Double ReDim DlTrungbinhtheoSovoick0(1 To 1, To sochukyquantrac - 1) For i = To sochukyquantrac - DlTrungbinhtheo2cklientiep(1, i) = TongDltheo2cklientiep(1, i) / sodiemquantrac DlTrungbinhtheoSovoick0(1, i) = TongDltheoSovoick0(1, i) / sodiemquantrac Next Dim Tocdoluntheo2Cklientiep() As Double ReDim Tocdoluntheo2Cklientiep(1 To 1, To sochukyquantrac - 1) Dim TocdoluntheosovoiCk0() As Double ReDim TocdoluntheosovoiCk0(1 To 1, To sochukyquantrac - 1) For i = To sochukyquantrac - Tocdoluntheo2Cklientiep(1, i) = DlTrungbinhtheo2cklientiep(1, i) / khoangthoigiantinh2Ck(1, i) TocdoluntheosovoiCk0(1, i) = DlTrungbinhtheoSovoick0(1, i) / khoangthoigiantinhtuCkdau(1, i + 1) Next Dim Dllechtheo2Cklientiep() As Double ReDim Dllechtheo2Cklientiep(1 To 1, To sochukyquantrac - 1) Dim DllechtheoSovoiCk0() As Double ReDim DllechtheoSovoiCk0(1 To 1, To sochukyquantrac - 1) Dim Smax1() As Double Dim Smax2() As Double L­¬ng Anh TuÊn - 64 - Lớp Trắc địa B-K48 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Khoa Trắc địa Đồ án tốt nghiệp Dim Smin1() As Double Dim Smin2() As Double ReDim Smax1(1 To 1, To sochukyquantrac - 1) ReDim Smax2(1 To 1, To sochukyquantrac - 1) ReDim Smin1(1 To 1, To sochukyquantrac - 1) ReDim Smin2(1 To 1, To sochukyquantrac - 1) For i = To sochukyquantrac - Smax1(1, i) = CviDl2Cklt(i, 1) Smax2(1, i) = CviDlSvCk0(i, 1) For j = To sodiemquantrac If CviDl2Cklt(i, j) >= Smax1(1, i) Then Smax1(1, i) = CviDl2Cklt(i, j) Else Smax1(1, i) = Smax1(1, i) End If If CviDlSvCk0(i, j) >= Smax2(1, i) Then Smax2(1, i) = CviDlSvCk0(i, j) Else Smax2(1, i) = Smax2(1, i) End If Next Next For i = To sochukyquantrac - Smin1(1, i) = CviDl2Cklt(i, 1) Smin2(1, i) = CviDlSvCk0(i, 1) For j = To sodiemquantrac If CviDl2Cklt(i, j)