BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP & PTNN TRƯỜNG ĐẠI HỌC THỦY LỢI NGUYỄN VĂN DŨNG NGHIÊN CỨU CƠ SỞ KHOA HỌC CỦA GIẢI PHÁP TƯỜNG ĐỈNH GIẢM SÓNG TRÀN TRÊN ĐÊ BIỂN LUẬN ÁN TIẾN SĨ KỸ THUẬT HÀ NỘI, NĂM 2017 BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP & PTNN TRƯỜNG ĐẠI HỌC THỦY LỢI NGUYỄN VĂN DŨNG NGHIÊN CỨU CƠ SỞ KHOA HỌC CỦA GIẢI PHÁP TƯỜNG ĐỈNH GIẢM SÓNG TRÀN TRÊN ĐÊ BIỂN Chuyên ngành: Xây dựng công trình thủy Mã số:62-58-40-01 NGƯỜI HƯỚNG DẪN KHOA HỌC: PGS.TS Lê Xuân Roanh GS.TS Thiều Quang Tuấn HÀ NỘI, NĂM 2017 i LỜI CAM ĐOAN Tác giả xin cam đoan công trình nghiên cứu thân tác giả Các kết nghiên cứu kết luận luận án trung thực, không chép từ nguồn hình thức Việc tham khảo nguồn tài liệu thực trích dẫn ghi nguồn tài liệu tham khảo theo quy định Tác giả luận án Nguyễn Văn Dũng i MỤC LỤC MỞ ĐẦU 1 Tính cấp thiết đề tài .1 Mục tiêu nghiên cứu Đối tượng phạm vi nghiên cứu .2 Nội dung nghiên cứu Cách tiếp cận phương pháp nghiên cứu Ý nghĩa khoa học thực tiễn .3 Cấu trúc luận án CHƯƠNG 1:TỔNG QUAN TÌNH HÌNH NGHIÊN CỨU VỀ SÓNG TRÀN QUA ĐÊ BIỂN 1.1 Tổng quan đê có tường đỉnh 1.1.1 Đê có tường đỉnh Hà Tĩnh 1.1.2 Đê biển có tường đỉnh Thanh Hóa .7 1.1.3 Đê biển có tường đỉnh Nam Định .7 1.1.4 Đê biển có tường đỉnh Hải Phòng 1.2 Tổng quan nghiên cứu sóng tràn qua đê biển 1.3 Ảnh hưởng tường đỉnh đê đến sóng tràn 12 1.3.1 Nghiên cứu TAW (2002) 12 1.3.2 Nghiên cứu Thiều Quang Tuấn (2009, 2013) 14 1.3.3 Nghiên cứu Nguyễn Văn Thìn (2014) 16 1.3.4 Nghiên cứu Koen Van Doorslear nnk (2015) 16 1.4 Áp lực sóng lên tường đỉnh đê 18 1.4.1 Nghiên cứu áp lực sóng lên tường đỉnh giới 18 1.4.2 Nghiên cứu áp lực sóng lên tường đỉnh Việt Nam 28 Kết luận chương 29 CHƯƠNG 2: XÁC ĐỊNH LƯU LƯỢNG SÓNG TRÀN TRUNG BÌNH QUA ĐÊ BIỂN CÓ TƯỜNG ĐỈNH VỚI MŨI HẮT SÓNG BẰNG THÍ NGHIỆM MÔ HÌNH VẬT LÝ .30 2.1 Cơ sở lý thuyết chung .30 2.1.1 Tương tự hình học 30 ii 2.1.2 Tương tự động học 30 2.1.3 Tương tự động lực học .31 2.2 Thí nghiệm nghiên cứu sóng tràn qua đê biển .31 2.2.1 Thiết bị thí nghiệm 31 2.2.2 Mô hình đê tham số thí nghiệm 31 2.2.3 Trình tự thí nghiệm 33 2.3 Phân tích đánh giá kết thí nghiệm .34 2.3.1 Kiểm định nghiên cứu Thiều Quang Tuấn (2013) trường hợp tường đỉnh mũi hắt sóng 34 2.3.2 Ảnh hưởng tham số tường đỉnh có mũi hắt sóng đến sóng tràn qua đê biển .35 2.4 Xây dựng phương pháp tính toán lưu lượng sóng tràn trung bình qua đê biển có tường đỉnh 41 2.4.1 Cơ sở lý thuyết hệ số ảnh hưởng thành phần mũi hắt 41 2.4.2 Mối quan hệ tham số a 41 2.4.3 Xây dựng công thức thực nghiệm .44 Kết luận chương 49 CHƯƠNG 3: NGHIÊN CỨU ÁP LỰC SÓNG LỚN NHẤT LÊN TƯỜNG ĐỈNH BẰNG MÔ HÌNH VẬT LÝ MÁNG SÓNG 50 3.1 Giới thiệu thiết bị đo áp lực sóng 50 3.1.1 Đầu đo áp lực 50 3.1.2 Máy đo áp lực (máy đa kênh) .51 3.2 Các kịch thí nghiệm trình tự thí nghiệm .51 3.2.1 Các kịch thí nghiệm .51 3.2.2 Công tác chuẩn bị 52 3.2.3 Trình tự thí nghiệm 53 3.3 Phân tích đánh giá kết thí nghiệm .54 3.3.1 Kết thí nghiệm .54 3.3.2 Phân tích áp lực lớn 57 3.4 Xây dựng phương pháp tính toán áp lực sóng lên tường đỉnh 59 3.4.1 Áp lực sóng vị trí Y = cm .59 3.4.2 Xây dựng biểu đồ đường bao phân bố áp lực sóng lớn 71 iii 3.4.3 Xây dựng biể đồ phân bố áp lực sóng thực tế lớn 76 Kết luận chương 77 CHƯƠNG 4: ỨNG DỤNG KẾT QUẢ NGHIÊN CỨU TRONG THIẾT KẾ ĐÊ BIỂN BIỂN HẢI NINH – THANH HÓA .79 4.1 Chọn địa điểm áp dụng 79 4.1.1 Sơ lược đê Hải Ninh 79 4.1.2 Hiện trạng đê biển Hải Ninh 81 4.2 Ứng dụng kết nghiên cứu thiết kế mặt cắt ngang đê biển Hải Ninh 83 4.2.1 Cấp công trình tần suất thiết kế .83 4.2.2 Điều kiện biên thủy lực .83 4.2.3 Xác định cao trình đỉnh đê 84 4.2.4 Đánh giá hiệu phương án thiết kế 86 4.3 Kiểm tra ổn định tường đỉnh theo kết nghiên cứu .87 4.3.1 So sánh kết xác định áp lực sóng 87 4.3.2 Kiểm tra ổn định tường đỉnh đê .88 Kết luận chương 4: .91 KẾT LUẬN VÀ KIẾN NGHỊ .92 I Kết đạt luận án .92 Nghiên cứu tổng quan 92 Nghiên cứu thực nghiệm mô hình vật lý 92 Nghiên cứu ứng dụng vào công trình thực tế 92 II Những đóng góp luận án 93 III Tồn hướng phát triển 93 Những tồn 93 Hướng phát triển .93 IV Kiến nghị 94 DANH MỤC CÁC CÔNG TRÌNH ĐÃ CÔNG BỐ 95 TÀI LIỆU THAM KHẢO 96 PHỤ LỤC 101 iv DANH MỤC HÌNH ẢNH Hình 1.1 Một số dạng điển hình kết cấu tường đỉnh đê biển Hình 1.2 Đê biển Hội Thống, Nghi Xuân, Hà Tĩnh Hình 1.3 Đê biển Phúc Long Nhượng, Cẩm Xuyên, Hà Tĩnh Hình 1.4 Đê biển Hậu Lộc, Thanh Hóa Hình 1.5 Đê biển Giao Thủy, Nam Định Hình 1.6 Đê biển Quất Lâm, Nam Định Hình 1.7 Đê biển Cát Hải, Hải Phòng Hình 1.8 Đê biển Afsluitdijk Hà Lan .9 Hình 1.9 Sóng tràn gây vỡ đê biển Hậu lộc, Thanh Hóa (bão số 7/2005) 10 Hình 1.10 Xác định độ dốc mái đê quy đổi có tường đỉnh 13 Hình 1.11 Biểu đồ áp lực sóng Goda 19 Hình 1.12 Sơ đồ xác định áp lực sóng Oumeraci, H (2001) 20 Hình 1.13 Hình minh họa áp lực sóng tác dụng lên tường 21 Hình 1.14 Sơ đồ áp lực sóng tác dụng lên tường đỉnh 22 Hình 1.15 Mặt cắt dọc đê có tường đỉnh máng sóng 23 Hình 1.16 Sơ đồ bố trí đầu đo áp lực .24 Hình 1.17 Áp lực sóng đo đầu áp lực theo thời gian 24 Hình 1.18 Áp lực sóng chu kỳ sóng 25 Hình 1.19 Áp lực sóng phân bố theo chiều cao tường 25 Hình 1.20 Đường hồi quy xác định áp lực sóng theo thời gian gia tăng áp lực .26 Hình 1.21 Mũi hắt sóng làm tăng áp lực so với mũi hắt sóng 27 Hình 1.22 Độ thị quan hệ lực sóng theo Hmo Rc 27 Hình 1.23 Phân bố áp lực sóng lên tường xung quanh thời điểm t* .28 Hình 2.1 Sơ đồ bố trí thí nghiệm 32 Hình 2.2 Sơ đồ mô hình tườngđỉnh có mũi hắt sóng nghiên cứu 33 Hình 2.3 Kiểm định ảnh hưởng tường đỉnh đến sóng tràn qua đê theo Thiều Quang Tuấn (2013), tường mũi hắt sóng, sóng vỡ 34 Hình 2.4 Ảnh hưởng mũi hắt tường đỉnh đến sóng tràn qua đê biển .35 Hình 2.5 Ảnh hưởng bề rộng thềm trước đến lưu lượng tràn tường đỉnh mũi hắt sóng 36 v Hình 2.6 Ảnh hưởng bề rộng thềm trước đến lưu lượng tràn tường đỉnh có mũi hắt sóng (β = 450) 36 Hình 2.7 Ảnh hưởng bề rộng thềm trước đến lưu lượng tràn tường đỉnh có mũi hắt sóng (β = 900) 37 Hình 2.8 Ảnh hưởng bề rộng thềm trước đến sóng tràn qua đê có mũi hắt sóng 38 Hình 2.9 Ảnh hưởng góc mũi hắt sóng đến lưu lượng sóng tràn qua đê 39 Hình 2.10: Ảnh hưởng chiều cao mũi hắt tương đối hn/W đến lưu lượng sóng tràn .40 Hình 2.11 Sự phụ thuộc a với bề rộng thềm góc mũi hắt sóng 45 Hình 2.12 Sự phụ thuộc a với bề rộng thềm chiều cao mũi tương đối .46 Hình 2.13 Đường hồi quy hàm số thực nghiệm xác định hệ số ảnh hưởng thềm mũi hắt a .47 Hình 2.14 Hệ số ảnh hưởng tổng hợp tường đỉnh có mũi hắt sóng v thực đo tính toán 48 Hình 2.15 Sóng tràn qua đê biển với hệ số ảnh hưởng tổng hợp cho tường đỉnh có mũi hắt sóng 48 Hình 3.1 Kết nối đầu đo PDB với máy đo áp lực truyền liệu vào máy tính .50 Hình 3.2 Sơ đồ bố trí thí nghiệm đo áp lực sóng 52 Hình 3.3 Kiểm tra đầu đo áp lực 53 Hình 3.4 Lắp đặt đầu đo áp lực PDB-200KPa lên tường đỉnh 53 Hình 3.5 Tín hiệu áp lực ghi lại dạng hiệu điện (mV) 55 Hình 3.6 Áp lựcsóng theo đợt sóng .55 Hình 3.7 Biểu đồ phân bố áp lực sóng số thời điểm 56 Hình 3.8 Áp lực sóng tác dụng lên tường theo thời gian vị trí đo 57 Hình 3.9 Xác định đỉnh áp lực sóng tác dụng lên tường theo thời gian 58 Hình 3.10 Áp lực sóng lớn kịch thí nghiệm Test_196 58 Hình 3.11 Ảnh hưởng Hm0 đến áp lực sóng lên tường (Đường mầu đỏ - mũi tên thể xu quan hệ) 59 Hình 3.12 Ảnh hưởng Tp đến áp lực sóng lên tường 60 Hình 3.13 Ảnh hưởng Rc đến áp lực sóng lên tường 60 Hình 3.14 Ảnh hưởng bề rộng thềm trước tường (S) đến áp lực sóng 61 Hình 3.15 Ảnh hưởng độ dốc mái đê phía biển đến áp lực sóng lên tường 61 Hình 3.16 Ảnh hưởng W đến áp lực sóng lên tường 62 vi Hình 3.17 Sự phụ thuộc p/(gρHm0) với độ cao lưu không tương đối 65 Hình 3.18 Sự phụ thuộc p/(gρHm0) với bề rộng thềm tương đối 66 Hình 3.19 Sự phụ thuộc p/(gρHm0) với số Iribarren ξm 66 Hình 3.20 Sự phụ thuộc p/(gρHm0) với chiều cao tường tương đối 67 Hình 3.21Đường hồi quy áp lực sóng Y = 1cm, tường mũi hắt sóng 68 Hình 3.22 Đường hồi quy áp lực sóng có mũi hắt sóng .68 Hình 3.23 Sự phụ thuộc áp lực sóng vào góc mũi hắt sóng 69 Hình 3.24 Sự phụ thuộc áp lực vào chiều cao mũi tương đối .70 Hình 3.25 Đường hồi quy hàm số thực nghiệm xác định áp lực sóng vị trí Y = 1cm xét đến mũi hắt sóng 71 Hình 3.26 Đường hồi quy hàm số thực nghiệm áp lực sóng vị trí Y = 3cm 72 Hình 3.27 Đường hồi quy hàm số thực nghiệm áp lực sóng vị trí Y = 4,5cm 72 Hình 3.28 Đường hồi quy hàm số thực nghiệm lực sóng vị trí Y = 6cm 73 Hình 3.29 Đường hồi quy hàm số thực nghiệm lực sóng vị trí Y = 6,5 cm .73 Hình 3.30 Đường hồi quy hàm số thực nghiệm áp lực sóng vị trí Y = 8,5cm 74 Hình 3.31 Đường hồi quy hàm số thực nghiệm lực sóng vị trí Y = 9,5cm 74 Hình 3.32 Đường hồi quy hàm số thực nghiệm áp lực sóng vị trí Y = 11,5cm .75 Hình 3.33 Đường hồi quy hàm số thực nghiệm xác định hệ số Cp theo chiều cao tương đối vị trí áp lực chiều cao tường .76 Hình 3.34 Đường hồi Đường hồi quy xác định tương quan lực lớn (Fmax) lực đồng thời xảy áp lực p1/250 (F(p1/250)) 77 Hình 4.1 Vị trí tuyến đê Hải Ninh 78 Hình 4.2 Bản đồ vị trí đê thiết kế 82 Hình 4.3 Tuyến đê biển Hải Ninh xây dựng 82 Hình 4.4 Mặt cắt ngang thiết kế 83 Hình 4.5 Đường tần suất mực nước tổng hợp điểm MC20 .84 Hình 4.6 Mặt cắt ngang đê biển Hải Ninh (phương án đề xuất) 86 Hình 4.7 Chi tiết tường đỉnh đê .86 Hình 4.8 Biểu đồ áp lực sóng lên tường đỉnh đê Hải Ninh 88 Hình 4.9 Sơ đồ lực tác dụng lên tường đỉnh 88 vii DANH MỤC BẢNG BIỂU Bảng 1.1 Các hệ số thực nghiệm công thức Owen (1980) 11 Bảng 2.1 Tổng hợp chương trình thí nghiệm 33 Bảng 3.1 Số liệu thí nghiệm kịch Test_196 54 Bảng 3.2 Giá trị hệ số Cp theo vị trí tường 75 Bảng 4.1 Kết tính cao trình đỉnh đê .85 Bảng 4.2 Kết so sánh cao trình đỉnh đê .87 Bảng 4.3 Áp lực sóng lên tường theo phương páp 87 viii 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 0.09 0.09 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.09 0.09 0.09 0.09 0.09 0.09 0.09 0.09 0.09 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.09 0.25 0.25 0 0.1 0.1 0.1 0.25 0.25 0.25 0 0.1 0.1 0.25 0.25 0.25 0 0.1 0.1 0.1 0.25 0.25 0.25 0 0.1 0.1 0.1 0.25 0.25 0.25 0 0.1 0.1 0.1 0.25 0.25 0.25 0 0 0 0 0 0 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 2.61 2.69 1.903 2.078 2.209 2.12 2.289 2.406 2.557 2.7 2.778 1.7 1.876 2.015 2.05 2.182 2.22 2.382 2.49 1.779 1.956 2.094 1.957 2.133 2.263 2.303 2.468 2.574 1.854 2.033 2.169 2.034 2.211 2.34 2.381 2.548 2.653 1.728 1.901 2.037 1.939 2.106 2.229 2.375 2.514 2.596 1.819 113 0 0 0 0 0 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.19 0.19 0.22 0.22 0.22 0.22 0.22 0.22 0.22 0.22 0.22 0.21 0.21 0.21 0.21 0.21 0.21 0.21 0.21 0.24 0.24 0.24 0.24 0.24 0.24 0.24 0.24 0.24 0.27 0.27 0.27 0.27 0.27 0.27 0.27 0.27 0.27 0.16 0.16 0.16 0.16 0.16 0.16 0.16 0.16 0.16 0.19 5.44 8.78 4.36 5.66 7.77 5.14 5.99 5.18 5.26 7.67 12.05 7.18 8.62 8.74 3.55 6.03 14.76 1.57 3.47 12.04 4.64 2.31 4.01 2.09 12 23.97 1.61 6.09 11.38 2.93 10.99 14.62 9.8 15.37 23.88 12.51 19.39 31.98 5.903 22.81 18.33 4.68 13.52 11.83 16.54 26.88 10.35 12.39 25.15 6.16 6.16 14.61 4.91 14.83 15.23 2.94 12.34 7.91 13.09 12.66 24.56 3.11 12.33 17.8 8.24 13.33 17.48 16.34 12.09 11.12 12.2 20.95 27.11 16.56 13 11.51 10.44 14.2 10.8 8.78 14.6 14.6 5.75 8.17 12.2 4.49 10.6 10.4 6.63 11.4 16.4 8.57 13.4 9.43 9.81 11 6.13 9.65 13.5 4.76 5.9 10.8 5.46 8.4 12.5 8.68 9.49 15.8 7.39 10 13.7 7.29 9.52 11.7 14.1 16.1 26.5 12.5 14.3 15.1 7.53 8.49 12.1 9.08 30.777 36.514 34.156 53.103 73.295 37.675 55.478 85.056 22.672 66.938 62.148 13.939 16.979 29.017 10.153 22.855 30.217 7.684 20.842 22.267 34.181 32.629 60.359 14.140 39.361 60.739 22.425 36.444 49.382 19.012 28.280 36.219 14.914 23.985 41.088 11.872 14.595 21.162 22.723 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 0.09 0.09 0.09 0.09 0.09 0.09 0.09 0.09 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.09 0.09 0.09 0.09 0.09 0.09 0.09 0.09 0.09 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.06 0.06 0.06 0.06 0.06 0 0.1 0.1 0.1 0.25 0.25 0.25 0 0.1 0.1 0.1 0.25 0.25 0.25 0 0.1 0.25 0.25 0.25 0 0.1 0.1 0.1 0.25 0.25 0.25 0 0.1 0.1 0.1 0.25 0.25 0.25 0 0.1 0.1 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 1.992 2.125 2.033 2.201 2.32 2.47 2.61 2.69 1.903 2.078 2.209 2.12 2.289 2.406 2.557 2.7 2.778 1.7 1.876 2.015 1.876 2.22 2.382 2.49 1.779 1.956 2.094 1.957 2.133 2.263 2.303 2.468 2.574 1.854 2.033 2.169 2.034 2.211 2.34 2.381 2.548 2.653 1.728 1.901 2.037 1.939 2.106 114 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.6 0.6 0.6 0.6 0.6 0.19 0.19 0.19 0.19 0.19 0.19 0.19 0.19 0.22 0.22 0.22 0.22 0.22 0.22 0.22 0.22 0.22 0.21 0.21 0.21 0.21 0.21 0.21 0.21 0.24 0.24 0.24 0.24 0.24 0.24 0.24 0.24 0.24 0.27 0.27 0.27 0.27 0.27 0.27 0.27 0.27 0.27 0.16 0.16 0.16 0.16 0.16 9.27 10.40 9.56 9.36 13.45 5.08 6.32 10.88 4.07 3.53 5.72 8.22 14.07 25.92 1.18 2.6 16.22 6.01 6.71 7.25 4.19 9.50 13.37 3.73 5.81 5.26 3.3 15.27 23 6.1 4.8 7.42 1.97 7.69 19.03 15.72 18.39 6.13 11.91 24.38 5.14 6.07 7.9 5.22 23.27 36.4 16.24 14.1 32.97 20.89 20.66 29.03 18.35 17.75 19.81 5.74 9.19 13.03 10.98 3.59 7.81 15.76 6.18 16.67 9.59 9.71 6.67 9.75 9.893 24.74 30.53 10.31 24.8 25.04 5.285 13.69 19.76 11.52 20.19 18.17 12.23 16.51 10.4 16.9 6.21 10.2 12.7 5.98 7.7 11.6 11 12.3 14.5 8.16 12.3 13.3 8.95 10.6 12 5.96 13.8 18 11 5.73 9.03 13.8 5.09 9.42 20.1 4.87 9.34 11 10.1 9.55 4.72 12.4 25.4 9.22 14.5 13.3 6.43 8.25 13.2 16 18.3 29.5 9.47 6.11 54.330 80.384 41.522 46.119 87.493 45.289 47.542 72.705 13.939 16.979 29.017 10.153 22.855 30.217 7.684 20.842 22.267 9.916 17.869 34.306 11.883 28.675 23.597 19.243 15.961 19.514 23.848 65.147 93.401 31.017 60.020 61.062 16.794 36.906 60.113 15.816 24.214 44.118 15.149 20.054 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 3 3 3 3 3 3 3 3 3 3 4 4 4 4 4 4 4 4 4 4 4 4 4 0.06 0.06 0.06 0.09 0.09 0.09 0.09 0.09 0.09 0.09 0.09 0.09 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.09 0.09 0.09 0.09 0.09 0.09 0.09 0.09 0.09 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.25 0.25 0.25 0 0.1 0.1 0.1 0.25 0.25 0.25 0 0.1 0.1 0.1 0.25 0.25 0.25 0 0.1 0.1 0.1 0.25 0.25 0.25 0 0.1 0.1 0.1 0.25 0.25 0.25 0 0.1 0.1 0.1 0.25 0.25 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 0 0 0 0 0 0 0 0 0 0 0 0 0 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.375 2.514 2.596 1.819 1.992 2.125 2.033 2.201 2.32 2.47 2.61 2.69 1.903 2.078 2.209 2.12 2.289 2.406 2.557 2.7 2.778 1.293 1.425 1.529 1.392 1.523 1.623 1.573 1.699 1.788 1.362 1.495 1.597 1.464 1.596 1.694 1.65 1.776 1.862 1.428 1.562 1.663 1.533 1.666 1.762 1.722 1.85 115 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0 0 0 0 0 0 0 0 0 0 0 0 0 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.16 0.16 0.16 0.19 0.19 0.19 0.19 0.19 0.19 0.19 0.19 0.19 0.22 0.22 0.22 0.22 0.22 0.22 0.22 0.22 0.22 0.21 0.21 0.21 0.21 0.21 0.21 0.21 0.21 0.21 0.24 0.24 0.24 0.24 0.24 0.24 0.24 0.24 0.24 0.27 0.27 0.27 0.27 0.27 0.27 0.27 0.27 5.04 4.43 6.61 2.91 7.69 7.33 4.69 5.30 5.33 11.84 15.03 27.71 7.09 16.78 24.31 9.88 12.57 4.35 2.06 1.15 1.5 0.16 0.9 2.64 1.02 1.33 3.51 1.13 1.83 2.78 1.09 2.63 4.14 1.04 1.88 3.19 3.02 1.15 1.85 4.8 0.72 1.11 0.2 2.91 16.51 19.98 12.68 7.41 12.32 25.04 7.75 8.88 9.86 10.5 7.72 8.02 12.6 27.54 35.97 6.21 10.97 30.14 14.79 21.86 21.94 2.56 2.49 4.34 3.46 4.56 8.77 3.72 4.27 7.88 1.65 1.67 3.11 1.41 2.92 4.33 2.01 2.45 3.73 1.87 1.09 2.92 2.04 3.95 3.43 0.43 1.77 6.11 12.9 13.2 8.07 12.7 19.9 7.8 10.3 13.8 8.98 7.66 12.5 8.64 17.9 15.4 7.63 11 14.2 6.38 9.7 12.4 4.07 4.55 10.1 4.53 4.55 5.67 2.9 5.15 7.15 4.4 5.73 9.18 1.82 2.33 4.09 3.95 4.24 4.4 2.5 2.17 4.59 5.38 4.93 12.4 0.63 2.9 23.438 19.712 15.974 19.704 38.348 75.239 95.946 23.150 41.395 82.455 38.439 55.661 52.743 4.038 3.908 7.528 4.282 5.282 9.199 4.070 5.448 9.429 6.459 7.980 13.315 4.261 8.115 12.670 6.612 8.147 11.227 9.749 5.724 12.757 14.944 14.017 20.824 1.753 9.818 25.954 5.054 7.188 11.830 7.148 8.968 12.305 4.796 185 186 187 188 189 190 191 192 193 194 195 185 186 187 188 189 190 191 192 193 194 195 4 4 4 4 4 0.12 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.09 0.25 0 0.1 0.1 0.1 0.25 0.25 0.25 0 0 0 0 0 0 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 2.058 1.458 1.73 2.058 1.458 1.73 2.058 1.458 1.73 2.058 1.458 1.934 1.351 1.447 1.548 1.473 1.563 1.656 1.705 1.776 1.851 1.432 0 0 0 0 0 0.55 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.27 0.16 0.16 0.16 0.16 0.16 0.16 0.16 0.16 0.16 0.19 0.83 1.95 1.29 4.19 1.31 1.8 2.18 2.05 1.55 8.45 1.64 8.71 3.71 5.55 7.39 5.89 8.36 11.37 4.19 7.29 9.32 2.69 6.75 4.68 7.29 13.5 6.58 6.44 9.23 3.12 6.21 7.81 5.83 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 0.09 0.09 0.09 0.09 0.09 0.09 0.09 0.09 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.09 0.09 0.09 0.09 0.09 0.09 0.09 0.09 0.09 0.12 0 0.1 0.1 0.1 0.25 0.25 0.25 0 0.1 0.1 0.1 0.25 0.25 0.25 0 0.1 0.1 0.1 0.25 0.25 0.25 0 0.1 0.1 0.1 0.25 0.25 0.25 0 0 0 0 0 0 0 0 0 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 1.73 2.058 1.458 1.73 2.058 1.458 1.73 2.058 1.458 1.73 2.058 1.458 1.73 2.058 1.458 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.527 1.625 1.559 1.646 1.736 1.796 1.865 1.935 1.51 1.603 1.699 1.64 1.725 1.813 1.882 1.949 2.016 1.293 1.425 1.529 1.392 1.523 1.623 1.573 1.699 1.788 1.362 1.495 1.597 1.464 1.596 1.694 1.65 1.776 1.862 1.428 0 0 0 0 0 0 0 0 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.19 0.19 0.19 0.19 0.19 0.19 0.19 0.19 0.22 0.22 0.22 0.22 0.22 0.22 0.22 0.22 0.22 0.21 0.21 0.21 0.21 0.21 0.21 0.21 0.21 0.21 0.24 0.24 0.24 0.24 0.24 0.24 0.24 0.24 0.24 0.27 2.46 3.65 1.62 2.54 3.61 1.21 2.00 3.27 3.3 5.11 9.03 0.97 1.07 6.84 2.53 0.67 4.36 3.09 5.98 3.3 6.09 10.48 1.64 5.15 8.12 12.33 14.58 30.53 1.29 4.78 2.48 4.54 2.31 5.52 2.17 4.61 11.17 5.4 7.86 8.65 5.35 3.61 4.5 1.94 3.07 3.93 2.24 5.11 7.03 1.66 3.39 5.19 1.48 8.16 14.5 8.15 7.51 3.91 10.5 7.7 4.72 12.4 9.66 3.41 7.37 12.1 6.67 5.6 9.26 3.18 6.26 12.3 3.3 4.31 7.47 14 4.99 5.22 3.65 3.83 7.63 3.12 3.69 4.54 2.08 4.6 5.77 1.52 116 5.12 5.13 8.63 3.76 6.62 10.30 0.83 4.41 6.04 0.64 8.052 12.183 9.428 12.315 22.040 11.951 17.467 23.865 6.165 16.675 20.186 34.243 48.639 83.706 7.186 18.585 24.832 18.987 11.188 25.648 9.916 17.869 34.306 11.883 28.675 23.597 19.243 15.961 6.737 8.110 13.116 6.152 10.852 15.031 3.614 8.713 12.164 4.484 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.09 0.09 0.09 0.09 0.09 0.09 0.09 0.09 0.09 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.09 0.09 0.09 0 0.1 0.1 0.1 0.25 0.25 0.25 0 0.1 0.1 0.1 0.25 0.25 0.25 0 0.1 0.1 0.1 0.25 0.25 0.25 0 0.1 0.1 0.1 0.25 0.25 0.25 0 0.1 0.1 0.1 0.25 0.25 0.25 0 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 90 90 90 90 90 90 90 90 90 90 90 90 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.458 1.73 2.058 1.458 1.73 2.058 1.458 1.73 2.058 1.458 1.73 2.058 1.458 1.73 2.058 1.458 1.73 2.058 1.458 1.73 2.058 1.458 1.73 2.058 1.458 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.562 1.663 1.533 1.666 1.762 1.722 1.85 1.934 1.351 1.447 1.548 1.473 1.563 1.656 1.705 1.776 1.851 1.432 1.527 1.625 1.559 1.646 1.736 1.796 1.865 1.935 1.51 1.603 1.699 1.64 1.725 1.813 1.882 1.949 2.016 1.293 1.425 1.529 1.392 1.523 1.623 1.573 1.699 1.788 1.362 1.495 1.597 117 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.27 0.27 0.27 0.27 0.27 0.27 0.27 0.27 0.16 0.16 0.16 0.16 0.16 0.16 0.16 0.16 0.16 0.19 0.19 0.19 0.19 0.19 0.19 0.19 0.19 0.19 0.22 0.22 0.22 0.22 0.22 0.22 0.22 0.22 0.22 0.21 0.21 0.21 0.21 0.21 0.21 0.21 0.21 0.21 0.24 0.24 0.24 0.7 1.38 1.38 0.34 7.57 1.25 3.47 5.68 5.24 6.74 8.79 4.00 8.28 3.78 1.10 1.97 6.71 1.96 1.23 6.74 2.2 5.7 10.55 0.28 1.00 1.80 2.68 5.27 3.61 2.98 4.48 1.47 2.29 5.1 6.787 12.08 6.37 5.4 7.47 13.82 6.27 4.33 10.87 3.67 5.09 9.36 2.24 5.01 7.66 4.48 3.47 9.6 3.45 7.06 8.16 4.43 4.73 8.24 16.45 14.67 27.37 12.56 16.11 13.51 3.12 7.29 10.55 1.77 4.28 6.69 1.86 4.74 5.17 1.73 2.63 3.95 2.72 3.15 2.41 4.42 7.86 2.11 3.08 4.68 5.82 8.9 21.1 7.59 4.48 14.5 4.77 5.23 8.31 5.55 8.02 12.2 4.04 6.24 6.53 4.27 5.85 10.4 2.62 4.04 11.6 4.04 5.33 5.35 5.06 10.9 3.22 6.77 10.1 1.94 3.24 7.68 3.21 4.06 5.78 3.52 3.61 8.55 8.276 11.732 5.875 8.551 21.829 14.523 26.045 19.485 8.445 13.821 20.482 9.356 11.497 16.643 7.224 12.649 18.940 11.523 12.942 30.326 33.545 34.259 59.264 26.608 37.087 35.127 5.439 7.528 11.774 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 0.09 0.09 0.09 0.09 0.09 0.09 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.09 0.09 0.09 0.09 0.09 0.09 0.09 0.09 0.09 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.1 0.1 0.1 0.25 0.25 0.25 0 0.1 0.1 0.1 0.25 0.25 0.25 0 0.1 0.1 0.1 0.25 0.25 0.25 0 0.1 0.1 0.1 0.25 0.25 0.25 0 0.1 0.1 0.1 0.25 0.25 0.25 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 0.124 0.148 0.178 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.422 1.73 2.058 1.458 1.73 2.058 1.458 1.73 2.058 1.458 1.73 2.058 1.458 1.73 2.058 1.458 1.73 2.058 1.458 1.73 2.058 1.458 1.73 2.058 1.458 1.73 2.058 1.458 1.73 2.058 1.464 1.596 1.694 1.65 1.776 1.862 1.428 1.562 1.663 1.533 1.666 1.762 1.722 1.85 1.934 1.351 1.447 1.548 1.473 1.563 1.656 1.705 1.776 1.851 1.432 1.527 1.625 1.559 1.646 1.736 1.796 1.865 1.935 1.51 1.603 1.699 1.64 1.725 1.813 1.882 1.949 2.016 118 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.55 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.24 0.24 0.24 0.24 0.24 0.24 0.27 0.27 0.27 0.27 0.27 0.27 0.27 0.27 0.27 0.16 0.16 0.16 0.16 0.16 0.16 0.16 0.16 0.16 0.19 0.19 0.19 0.19 0.19 0.19 0.19 0.19 0.19 0.22 0.22 0.22 0.22 0.22 0.22 0.22 0.22 0.22 1.41 5.31 5.54 3.25 4.38 9.48 3.16 0.82 1.52 1.51 1.52 2.86 3.78 0.25 0.68 3.69 3.71 5.50 4.87 3.83 5.45 3.51 5.29 7.39 0.68 4.35 2.35 0.21 5.35 2.54 1.52 0.97 6.98 1.94 1.94 3.63 2.74 0.7 4.66 1.67 3.87 6.1 3.794 16.91 8.92 0.96 7.24 1.63 7.96 12.64 23.69 11.41 6.28 9.52 5.16 7.19 8.15 5.21 5.16 8.27 4.14 4.56 6.37 5.02 4.43 9.26 3.46 6.47 13.23 3.46 6.97 14.39 3.75 5.02 22.02 2.74 2.74 3.78 2.95 1.32 5.09 2.23 1.87 3.19 2.59 3.79 5.38 2.12 4.84 4.59 7.52 10.5 21.9 9.73 4.89 8.91 5.21 6.5 11.9 6.75 10 12.4 2.56 4.08 9.44 7.13 8.63 3.52 5.99 7.77 6.4 4.1 6.62 6.57 6.34 8.68 4.616 6.176 8.855 6.379 3.446 12.481 7.281 9.023 14.624 10.097 33.601 22.546 6.398 17.323 7.837 11.932 14.160 19.977 8.294 9.404 15.795 11.126 11.980 19.183 9.878 20.034 31.906 12.429 19.694 32.861 14.009 15.553 51.000 Phụ lục 4: Code matlap phân tích áp lực sóng clearall; clc; [inname, inpath] = uigetfile('*.DAT','file to read'); infname=[inpath inname]; if infname~0 a=load(infname); np=length(a(1,:))-1; linstyl={':r','-.k','-g',' b','-m','-c','-y','-r','-b'}; leg={'Probe No.1','Probe No.2','Probe No.3', 'Probe No.4', 'Probe No.5'}; %% Overview of all channel figure(1); for i=1:np plot(a(:,1),a(:,i+1),linstyl{i}); holdon end legend(leg); prompt = {'Time step to start','Time step to end','Time conversion factor(-)','Pressure conversion factor(KPa/mV)','Reference channel','Peak interval (s)','Positve pressure analysis (1-YES/0-NO)','Write output (1YES/0-NO)','HPF frequency (Hz)','LPF frequency'}; dlg_title = 'Wave on wall pressure analysis'; def={'-1','-1','0.001','200','3','2.5','1','1','0.1','49'}; Inparm = inputdlg(prompt,dlg_title,1,def,'off'); uiwait; if length(Inparm)>0 start= str2num(Inparm{1}); endt= str2num(Inparm{2}); tscl= str2num(Inparm{3}); pscl= str2num(Inparm{4}); n_r= str2num(Inparm{5}); %reference channel for finding peak pressure T1=str2num(Inparm{6}); T2=T1; neg=str2num(Inparm{7}); writeans=str2num(Inparm{8}); hz=str2num(Inparm{9});%low-passed filter hzl=str2num(Inparm{10}); %high-pass filter end if start~=-1 a(1:start,:)=[]; end if endt~=-1 a(endt-start:end,:)=[]; end t=a(:,1)*tscl; delt=t(2)-t(1); p=ones(length(t),np); p=a(:,2:np+1); cleara; %% High-passed filter Wn_h=hz/(100/2); %cutoff frequency or interested signal above hz, sampling freq 100 Hz/2 [b,a] = butter(3,Wn_h,'high'); for i=1:np p_hf(:,i) = filter(b,a,p(:,i)); %high-pass filter end %% Low-passed filter Wn_l=hzl/(100/2); %cutoff frequency or interested signal above hz, sampling freq 100 Hz/2 [b,a] = butter(3,Wn_l,'low'); 119 for i=1:np p_lf(:,i) = filter(b,a,p_hf(:,i)); %low-pass filter end %% figure(1); plot(t,p(:,n_r),'b ',t,p_hf(:,n_r),'r-',t,p_lf(:,n_r),'g-','linewidth',1); legend({'Orginal','High-pass filter','Low-pass filter'}); p=p_lf; clearp_hf; clearp_lf; figure(2); for i=1:np plot(t,p(:,i),linstyl{i},'linewidth',1); holdon end figure(3); for i=1:np if neg==1 p(p(:,i)