ĐẠI HỌC QUỐC GIA TP.HCM TRƯỜNG ĐẠI HỌC BÁCH KHOA ĐOÀN TẤN VIỆT CƯỜNG ỨNG XỬ CỦA TƯỜNG VÂY HỐ ĐÀO CƠNG TRÌNH XÉT ĐẾN TẢI TRỌNG ĐỘNG DO XE VẬN CHUYỂN TRONG Q TRÌNH THI CƠNG Chun ngành: Kỹ Thuật Xây Dựng Mã số: 8580201 LUẬN VĂN THẠC SĨ TP HỒ CHÍ MINH, tháng 07 năm 2022 CƠNG TRÌNH ĐƯỢC HOÀN THÀNH TẠI: TRƯỜNG ĐẠI HỌC BÁCH KHOA –ĐHQG -HCM Cán hướng dẫn khoa học: TS Đỗ Thanh Hải PSG.TS Lương Văn Hải Cán chấm nhận xét 1: TS Nguyễn Đình Hùng Cán chấm nhận xét 2: TS Nguyễn Hồng Ân Luận văn thạc sĩ bảo vệ Trường Đại học Bách Khoa, ĐHQG Tp HCM ngày 22 tháng 07 năm 2022 Thành phần Hội đồng đánh giá luận văn thạc sĩ gồm: PGS.TS Đào Đình Nhân – Chủ tịch TS Liêu Xuân Q – Thư ký TS Nguyễn Đình Hùng – Phản biện TS Nguyễn Hồng Ân – Phản biện TS Nguyễn Thái Bình – Ủy viên Xác nhận Chủ tịch Hội đồng đánh giá LV Trưởng Khoa quản lý chuyên ngành sau luận văn sửa chữa (nếu có) CHỦ TỊCH HỘI ĐỒNG TRƯỞNG KHOA KỸ THUẬT XÂY DỰNG i ĐẠI HỌC QUỐC GIA TP.HCM CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM TRƯỜNG ĐẠI HỌC BÁCH KHOA Độc lập - Tự - Hạnh phúc NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên: ĐOÀN TẤN VIỆT CƯỜNG Ngày, tháng, năm sinh: 22/08/1996 Chuyên ngành: Kỹ thuật xây dựng MSHV: 2170197 Nơi sinh: Bình Thuận Mã số: 8580201 I TÊN ĐỀ TÀI: “Ứng xử tường vây hố đào cơng trình xét đến tải trọng động xe vận chuyển q trình thi cơng” (The behavior of the excavation diaphragm in consideration of the dynamic load carried) by the vehicle during the construction process II NHIỆM VỤ VÀ NỘI DUNG: Nhiệm vụ: Phân tích ứng xử tường vây hố đào cơng trình xét đến tải trọng động xe vận chuyển q trình thi cơng Nội dung: Mở đầu Chương 1: Tổng quan Chương 2: Cơ sở lý thuyết Chương 3: Tính tốn mơ phần mềm Matlab Chương 4: Ứng dụng khảo sát chuyển vị tường vây xét đến tải trọng xe vận chuyển Kết luận kiến nghị Tài liệu tham khảo Phụ lục III NGÀY GIAO NHIỆM VỤ: 14/02/2022 IV NGÀY HOÀN THÀNH NHIỆM VỤ: 17/06/2022 V CÁN BỘ HƯỚNG DẪN: TS ĐỖ THANH HẢI PSG.TS LƯƠNG VĂN HẢI Tp Hồ Chí Minh, ngày tháng năm 20 CÁN BỘ HƯỚNG DẪN CHỦ NHIỆM BỘ MÔN ĐÀO TẠO TS ĐỖ THANH HẢI PGS TS LƯƠNG VĂN HẢI TRƯỞNG KHOA KỸ THUẬT XÂY DỰNG ii LỜI CẢM ƠN Lời đầu tiên, xin chân thành cảm ơn quý thầy cô khoa Kỹ thuật Xây dựng Ban giám hiệu trường Đại học Bách Khoa TP.HCM tạo điều kiện cho thực luận văn Tơi xin bày tỏ lịng biết ơn chân thành sâu sắc tới TS Đỗ Thanh Hải PSG.TS Lương Văn Hải tận tình hướng dẫn tơi hồn thành luận văn Bên cạnh đó, tơi xin gửi lời tri ân tới Công ty CP Địa kỹ thuật Nền móng Thái Dương Hệ giúp tơi số liệu quan trắc để thực luận văn Cuối cùng, xin cảm ơn gia đình bạn bè ủng hộ suốt thời gian vừa qua Tp Hồ Chí Minh, ngày 17 tháng 06 năm 2022 Học viên ĐỒN TẤN VIỆT CƯỜNG ii TĨM TẮT Tường vây hố đào cơng trình chịu tác động gây tải trọng tải trọng khối đất, áp lực nước ngầm, tải trọng cơng trình lân cận, tải trọng phương tiện giao thông di chuyển đường gần với tầng hầm, tải trọng thi công bên công trường xung quanh hố đào tầng hầm Luận văn trình bày phân tích, đánh giá nhân tố khoảng cách đặt tải trọng động, chiều sâu hố đào, vận tốc xe chạy ảnh hưởng tới áp lực tải trọng động gây lên tường vây hố đào cơng trình Q trình nghiên cứu sử dụng phương pháp giải tích, sử dụng phần mềm MATLAB PLAXIS để tính tốn áp lực tường vây hố đào bị lực tải trọng xe tác động lên, từ phân tích ứng xử tường vây tác động trường hợp hoạt tải xe tác động khác Công trình Cityland Park Hills nằm số 18 Phan Văn Trị, phường 10, quận Gò Vấp TPHCM chọn để nghiên cứu Tải tác động tính tốn mơ với trường hợp Kết phân tích cho thấy vận tốc xe chạy khoảng 5-30 km/h, gây áp lực động lớn lên tường vây Ngoài khoảng cách tải trọng động xe gây gần cơng trình gây áp lực cho tường vây lớn Tuy nhiên, khoảng cách với tường vây từ 5m trở ảnh hưởng không tuân theo quy luật tuyến tính mà phụ thuộc vào cộng hưởng tần số xe tạo tần số riêng đất Khi tính tốn thiết kế tường vây không xét đến ảnh hưởng tải trọng xe, chuyển vị tường vây tính nhỏ giá trị quan trắc, kết thiết kế tính khơng an tồn cho tường vây Nghiên cứu toán động tải xe thiết kế để áp dụng vào mơ hình tính tốn tường vây an toàn cho tường vây, tránh xảy rủi ro thi cơng cơng trình hầm, nằm sát đường có xe tải nặng chạy sát đường nội xe vận chuyển vật tư chạy nhiều iii ABSTRACT The diaphragm wall of the excavation pit is affected by loads such as earth mass load, groundwater pressure, neighboring structures, and traffic loads when traveling on roads close to the basement, construction loads inside the construction site and around the excavation pit of the basement This thesis presents the analysis and evaluation of factors such as the distance of placing dynamic loads, depth of excavation, vehicle speed affecting the pressure caused by dynamic loads on the diaphragm wall of the excavation pit The research process used analytical methods, using MATLAB and PLAXIS software to calculate the pressure of the diaphragm wall of the excavation when being acted upon by the vehicle load, thereby analyzing the behavior of the diaphragm wall under the influence of different driving loads Cityland Park Hills project located at 18 Phan Van Tri, Ward 10, Go Vap District, Ho Chi Minh City was selected for research Impact loads are calculated and simulated with cases The analysis results show that the vehicle speed is in the range of 5-30 km/h, which will cause great dynamic pressure on the diaphragm wall In addition, the closer the dynamic load caused by the vehicle is to the construction site, the greater the pressure exerted on the diaphragm wall However, the distance from the diaphragm wall from 5m onwards, this influence does not follow the linear law but depends on the resonance of the frequency generated by the vehicle and the natural frequency of the ground When calculating the diaphragm wall design without considering the influence of vehicle load, the calculated diaphragm wall displacement is smaller than the observed value, the calculated design results will not be safe for the diaphragm wall Research on the dynamic problem of vehicle loading when designing to apply to the diaphragm wall calculation model, it will be safer for the diaphragm wall, avoiding risks when constructing tunnels, located close to the road with heavy trucks running or close to the internal road, the vehicle transporting materials runs a lot iv LỜI CAM ĐOAN Tôi xin cam đoan luận văn đề tài “Ứng xử tường vây hố đào cơng trình xét đến tải trọng động xe vận chuyển q trình thi cơng” cơng trình nghiên cứu cá nhân tơi thời gian qua Các kết quả, số liệu luận án trung thực chưa công bố cơng trình khoa học khác Tơi xin chịu trách nhiệm cơng việc thực Tp Hồ Chí Minh, ngày 17 tháng 06 năm 2022 Học viên ĐOÀN TẤN VIỆT CƯỜNG v MỤC LỤC MỞ ĐẦU Sự cần thiết đề tài Mục tiêu đề tài .1 Phương pháp nghiên cứu .1 CHƯƠNG TỔNG QUAN 1.1 Thực trạng nghiên cứu nước .2 1.2 Thực trạng nghiên cứu nước .2 CHƯƠNG CƠ SỞ LÝ THUYẾT 2.1 Lý thuyết đàn hồi tuyến tính 2.2 Các mơ hình đất .9 2.2.1 Mơ hình Mohr – Coulomb 2.2.2 Mơ hình Hyperbol .11 2.2.3 Mơ hình Cam-Clay cải tiến 12 2.2.4 Mô hình Hardening Soil 13 2.3 Mơ hình tải động 17 2.3.1 Mơ hình hệ dao động .17 2.3.2 Mô hình hệ dao động hoạt tải xe 18 2.4 Sự truyền tải động môi trường liên tục đàn hồi đẳng hướng 21 CHƯƠNG TÍNH TỐN MƠ PHỎNG BẰNG PHẦN MỀM MATLAB 22 3.1 Các toán đưa .22 3.2 Kiểm chứng mơ hình tính toán 22 3.2.1 Mục đích 22 3.2.2 Kết tính tốn 23 3.3 Đánh giá ảnh hưởng khoảng cách từ tải trọng xe đến tường vây lên áp lực đất 24 3.4 Đánh giá áp lực động tải trọng xe lên tường vây 26 3.5 Ảnh hưởng tải trọng động tác động lên hệ giằng tường vây .28 CHƯƠNG ỨNG DỤNG KHẢO SÁT CHUYỂN VỊ TƯỜNG VÂY KHI XÉT ĐẾN TẢI TRỌNG XE VẬN CHUYỂN 30 4.1 Giới thiệu cơng trình 30 vi 4.1.1 Tổng quan cơng trình 30 4.1.2 Khảo sát địa chất 31 4.1.3 Quan trắc chuyển vị tường vây 33 4.2 Tính tốn mơ 38 4.2.1 Thông số địa chất sử dụng cho mơ hình 38 4.2.2 Cừ Larsen 39 4.2.3 Sàn hầm .39 4.2.4 Phụ tải mặt đất .40 4.2.5 Trình tự thi cơng 41 4.2.6 Kết tính tốn 43 4.2.7 Tổng hợp kết 51 KẾT LUẬN VÀ KIẾN NGHỊ .54 Kết luận 54 Kết nghị .54 TÀI LIỆU THAM KHẢO .55 PHỤ LỤC 57 Phụ lục A.1: Kết quan trắc chuyển vị I2 (đơn vị mm) .57 Phụ lục A.2: Kết quan trắc chuyển vị I3 (đơn vị mm) .61 Phụ lục A.3: Kết quan trắc chuyển vị I4 (đơn vị mm) .65 Phụ lục B.1: Code Matlab “So sánh với kết Veletsos Younan (1994)” 69 Phụ lục B.2: Code Matlab “Ảnh hưởng khoảng cách từ tải trọng đến tường vây” 70 Phụ lục B.3: Code Matlab “Áp lực động tải trọng xe lên tường vây” 71 Phụ lục B.4: Code Matlab “Lực hệ giằng tải trọng xe tác động lên tường vây- trường hợp hệ giằng” 72 Phụ lục B.5: Code Matlab “Lực hệ giằng tải trọng xe tác động lên tường vây- trường hợp hệ giằng” 73 vii DANH MỤC HÌNH ẢNH Hình 1.1 Minh họa tác động tải trọng động lên tường chắn .3 Hình 1.2 Đồ thị thể thay đổi áp lực đất đỉnh tường chắn đất phụ thuộc vào tương quan tần số dao động điều hòa tần số riêng, với hệ số nhớt 𝛿 = 0.1 .6 Hình 1.3 Đồ thị thể thay đổi áp lực đất lên tường chắn đất theo chiều sâu phụ thuộc vào tương quan tần số dao động điều hòa tần số riêng, với hệ số nhớt 𝛿 = 0.1 Hình 1.4 Đồ thị thể thay đổi theo phương ngang ứng suất pháp phương ngang gây chuyển động điều hòa bề mặt lớp đất; 𝑣= 0.3; 𝛿 = 0.1 .7 Hình 1.5 Đồ thị thể ảnh hưởng cản nhớt vật liệu tới độ lớn lực cắt đáy gây dao động điều hòa; 𝑣 = 0.3 .8 Hình 1.6 Đồ thị thể so sánh đường cong đáp ứng tần số lực cắt tường tính tốn cách xem xét chế độ rung; .8 Hình 2.1 Vịng trịn Mohr ngưỡng dẻo .9 Hình 2.2 Các mặt bao phá hoại theo tiêu chuẩn Mohr – Coulomb [8] 10 Hình 2.3 Quan hệ ứng suất – biến dạng theo mơ hình đàn hồi – dẻo lý tưởng .11 Hình 2.4 Quan hệ ứng suất – biến dạng mô đun mơ hình Hypebol 12 Hình 2.5 Mặt dẻo mơ hình Cam-clay cải tiến 13 Hình 2.6 Mặt dẻo mơ hình HS khơng gian ứng suất (p-q) 14 Hình 2.7 Định nghĩa E50 Eur thí nghiệm nén ba trục nước .15 Hình 2.8 Định nghĩa 𝐸oed 𝑟𝑒𝑓 theo kết thí nghiệm nén trục 16 Hình 2.9 Mơ hình hệ dao động bậc tự .17 Hình 2.10 Ảnh hưởng trọng lực vào cân hệ dao động bậc tự do.18 Hình 2.11 Mơ hình thành phần toán ảnh hưởng tải trọng động 18 Hình 2.12 Phân bố tải trọng xe xuống mặt đường Huang (1993) 19 Hình 2.13 Phân bố tải trọng xe theo vệt bánh xe Tran (2004) 19 Hình 2.14 Giả thiết mặt đường thí nghiệm Perret (2003) .20 Hình 2.15 Tính chất biến dạng mặt đường chịu tải trọng xe (Perret, 2003) .21 Hình 3.1 Giao diện tổng 22 Hình 3.2 Chọn vào “So sánh với kết Veletsos Younan (1994)” để kết .23 viii LUẬN VĂN THẠC SĨ I2 HVTH: ĐOÀN TẤN VIỆT CƯỜNG Ngày đo 29/06/ 2018 03/07/ 2018 07/07/ 2018 11/07/ 2018 14/07/ 2018 17/07/ 2018 21/07/ 2018 24/07/ 2018 28/07/ 2018 Chiều sâu đo Chu kỳ 31 Chu kỳ 32 Chu kỳ 33 Chu kỳ 34 Chu kỳ 35 Chu kỳ 36 Chu kỳ 37 Chu kỳ 38 Chu kỳ 39 -0.5 61.9875 61.875 62.1375 61.05 62.025 63.45 62.7375 63.8625 62.55 -1 52.2375 52.3125 52.5375 51.5625 52.8 53.625 53.175 53.4375 52.725 -1.5 42.4875 42.7875 42.9375 42.075 43.575 43.8 43.6125 43.7625 42.9 -2 34.6125 35.0625 34.9875 34.4625 35.8875 35.7375 35.925 35.7375 34.9125 -2.5 29.475 30.0375 29.925 29.2875 30.825 30.5625 30.75 30.7875 29.925 -3 21 21.6 21.4125 20.8875 22.275 21.975 22.05 22.1625 21.2625 -3.5 7.5 7.95 7.8375 7.3875 8.5875 8.325 8.3625 8.475 7.575 -4 2.25 2.7 2.5125 2.1375 3.3 2.925 3.075 3.15 2.2875 -4.5 1.2 1.575 1.3875 0.975 2.025 1.6875 1.875 1.9875 1.0875 -5 0.7875 1.125 0.9375 0.45 1.5 1.0875 1.3125 1.35 0.525 -5.5 0.1875 0.6 0.375 -0.1875 0.825 0.4125 0.675 0.7125 -0.1125 -6 -0.225 0.1125 -0.1125 -0.675 0.225 -0.075 0.075 0.075 -0.7125 -6.5 -0.7875 -0.4875 -0.7125 -1.2375 -0.4125 -0.7125 -0.6 -0.5625 -1.35 -7 -1.0875 -0.8625 -1.05 -1.425 -0.7125 -0.975 -0.975 -0.9375 -1.6875 -7.5 -1.3875 -1.0875 -1.3125 -1.65 -1.2 -1.3875 -1.425 -1.3125 -1.95 -8 -1.725 -1.575 -1.7625 -1.9125 -1.6125 -1.725 -1.7625 -1.6875 -2.2125 -8.5 -1.725 -1.5375 -1.6875 -1.8 -1.65 -1.725 -1.725 -1.65 -1.9875 -9 -1.5375 -1.5 -1.5375 -1.5375 -1.5 -1.5 -1.5375 -1.4625 -1.6875 -9.5 0 0 0 0 60 LUẬN VĂN THẠC SĨ HVTH: ĐOÀN TẤN VIỆT CƯỜNG Phụ lục A.2: Kết quan trắc chuyển vị I3 (đơn vị mm) I3 Ngày đo 27/03/ 2018 30/03/ 2018 03/04/ 2018 05/04/ 2018 09/04/ 2018 13/04/ 2018 17/04/ 2018 20/04/ 2018 24/04/ 2018 27/04/ 2018 Chiều sâu đo Chu kỳ Chu kỳ Chu kỳ Chu kỳ Chu kỳ Chu kỳ Chu kỳ Chu kỳ Chu kỳ Chu kỳ 10 -0.5 -0.4875 1.1625 -0.675 0.8625 3.3 0.4875 1.0875 3.2625 -2.5125 -1.65 -1 0.4875 2.2125 1.1625 2.3625 4.425 1.65 2.2125 4.3125 -0.6375 1.3125 -1.5 1.4625 3.2625 3.9 5.55 2.8125 3.3375 5.3625 1.2375 4.275 -2 0.8625 2.6625 1.6875 1.9875 3.5625 1.2 1.425 2.925 -0.4875 2.4375 -2.5 -0.675 0.7875 -0.225 -0.3 1.2375 -0.9 -1.5 0.15 -2.8125 -0.4875 -3 -0.7125 0.975 -0.4875 -0.0375 1.6875 -0.15 -1.2 0.5625 -2.5125 -0.225 -3.5 -1.125 0.675 -0.8625 -0.525 0.9375 -0.7875 -1.9875 -0.4125 -1.95 -0.8625 -4 -1.35 -1.0125 -1.05 -0.9375 -0.9375 -1.2375 -1.8375 -1.5375 -2.925 -2.175 -4.5 -1.35 -0.9375 -1.0125 -0.8625 -0.75 -0.9 -1.4625 -1.1625 -2.4 -1.575 -5 -1.35 -1.0875 -1.125 -1.05 -0.75 -0.7875 -1.1625 -1.1625 -2.1 -1.35 -5.5 -1.05 -0.8625 -0.975 -0.6375 -1.8375 -1.0875 -1.4625 -1.5375 -2.3625 -1.7625 -6 -0.75 -0.6 -0.675 -0.1875 -0.9375 -0.6375 -0.3375 -0.825 -1.425 -0.7125 -6.5 -0.3 -0.1125 -0.375 -0.1125 -0.1125 -0.525 -0.2625 -0.7875 -1.0875 -0.4875 -7 -0.2625 -0.4125 -0.375 -0.1875 -0.1125 -0.3 -0.5625 -0.675 -0.7125 -0.4875 -7.5 -0.4875 -0.5625 -0.6375 -0.4875 -0.375 -0.4875 -0.7125 -1.0125 -0.975 -0.525 -8 -0.225 -0.375 -0.3 -0.3 -0.3 -0.5625 -0.9375 -0.825 -1.4625 -1.1625 -8.5 0.525 0.4875 0.375 0.3375 0.2625 -0.225 -0.825 -0.975 -1.575 -1.575 -9 0.4125 0.4125 0.4125 0.3 0.1875 -0.1875 -0.7125 -0.825 -1.35 -1.35 -9.5 0 0 0 0 0 61 LUẬN VĂN THẠC SĨ I3 HVTH: ĐOÀN TẤN VIỆT CƯỜNG Ngày đo 02/05/ 2018 04/05/ 2018 09/05/ 2018 11/05/ 2018 13/05/ 2018 14/05/ 2018 16/05/ 2018 18/05/ 2018 21/05/ 2018 23/05/ 2018 Chiều sâu đo Chu kỳ 11 Chu kỳ 12 Chu kỳ 13 Chu kỳ 14 Chu kỳ 15 Chu kỳ 16 Chu kỳ 17 Chu kỳ 18 Chu kỳ 19 Chu kỳ 20 -0.5 -2.0625 -2.55 -1.05 0.2625 0.7125 1.8 6.075 7.575 8.325 10.0875 -1 0.9375 1.275 1.65 3.3375 3.75 4.1625 8.175 10.5375 8.8125 11.2125 -1.5 3.9375 5.175 4.3125 6.4125 6.7875 6.525 10.2375 13.5 9.3 12.3375 -2 1.875 3.15 1.95 3.825 4.1625 3.75 7.65 10.9125 4.8 7.575 -2.5 -1.9125 -0.525 -2.7375 -0.6375 -0.675 -1.0125 2.925 6.525 -0.1875 2.9625 -3 -2.925 -1.1625 -4.9125 -1.9875 -2.2125 -2.85 1.6125 5.5125 -1.125 2.175 -3.5 -3.75 -1.9125 -5.25 -2.5125 -2.925 -3.2625 0.975 4.9125 -1.7625 1.3125 -4 -4.5375 -2.7375 -4.6125 -3.3375 -3.3 -2.6625 -1.1625 2.8875 -2.475 0.525 -4.5 -4.6875 -2.475 -4.3125 -3.0375 -2.925 -2.55 -1.0875 2.9625 -2.0625 0.7125 -5 -4.5375 -2.625 -4.3125 -2.775 -3.1125 -2.925 -1.725 2.2125 -2.4 -0.075 -5.5 -4.275 -2.8875 -3.525 -3 -2.925 -3.075 -1.9125 2.1375 -2.55 -0.3375 -6 -2.475 -1.0875 -1.5375 -0.825 -0.825 -1.125 0.0375 3.8625 -0.7125 1.35 -6.5 -1.95 -0.15 -0.4125 -0.7125 0.15 0.5625 1.65 5.2125 0.75 2.4 -7 -1.0875 -0.45 -0.6375 -0.825 -0.0375 0.6375 1.425 5.1375 0.75 2.2125 -7.5 -0.7875 -0.4125 -0.5625 -0.75 -0.3 0.1875 0.8625 4.5375 0.4125 1.425 -8 -1.1625 -0.7875 -0.8625 -1.2 -0.9375 -0.6375 -0.15 3.3375 -0.6 0.2625 -8.5 -1.6125 -1.35 -1.3875 -1.575 -1.725 -1.35 -1.1625 2.25 -1.6125 -1.1625 -9 -1.4625 -1.6125 -1.5 -1.875 -1.8375 -1.725 -1.7625 2.4375 -2.1 -1.8 -9.5 0 0 0 0 0 62 LUẬN VĂN THẠC SĨ I3 HVTH: ĐOÀN TẤN VIỆT CƯỜNG Ngày đo 25/05/ 2018 30/05/ 2018 02/06/ 2018 05/06/ 2018 08/06/ 2018 12/06/ 2018 15/06/ 2018 20/06/ 2018 22/06/ 2018 26/06/ 2018 Chiều sâu đo Chu kỳ 21 Chu kỳ 22 Chu kỳ 23 Chu kỳ 24 Chu kỳ 25 Chu kỳ 26 Chu kỳ 27 Chu kỳ 28 Chu kỳ 29 Chu kỳ 30 -0.5 15.6 12.9 19.2 19.5375 15.7875 16.0875 40.9125 38.5125 39.375 41.4 -1 13.6125 12.225 16.9125 17.25 14.175 14.7 35.1375 33.75 33.675 36 -1.5 11.5875 11.5875 14.625 14.9625 12.5625 13.275 29.3625 28.9125 27.9375 30.5625 -2 6.8625 6.75 9.3 9.4125 7.1625 7.6125 20.1 19.7625 18.5625 21.1125 -2.5 2.025 1.6125 4.425 4.5 2.1375 2.4 12.9375 12.525 10.9875 13.425 -3 1.3125 1.125 3.675 3.9 1.4625 1.8 9.4875 9.15 7.3125 9.75 -3.5 0.3375 -0.0375 2.5125 2.7375 0.3 0.525 5.2875 5.025 3.1125 5.5875 -4 -0.6 -2.475 1.4625 1.8 -0.6375 -0.6 3.225 2.8875 0.8625 3.15 -4.5 -0.45 -2.0625 1.65 1.95 -0.3375 -0.225 3.3 2.925 0.9 3.0375 -5 -1.2 -2.6625 0.7875 1.125 -1.0875 -0.9 2.025 1.6875 -0.225 1.7625 -5.5 -0.975 -2.7375 0.15 0.375 -1.425 -1.35 1.2375 0.8625 -0.1125 1.1625 -6 0.8625 -0.525 1.9125 2.175 0.75 0.75 3.2625 2.9625 2.3625 3.2625 -6.5 2.025 0.8625 3.1125 2.3625 2.2125 4.5375 4.6875 4.7625 4.65 -7 1.9875 2.025 3.0375 2.9625 2.7 2.85 4.35 4.65 4.9125 4.65 -7.5 1.2375 1.5375 2.025 1.9125 1.725 1.8375 3.075 3.1875 3.5625 3.225 -8 -0.15 0.225 0.4125 0.2625 0.15 0.225 1.0125 0.975 1.0875 0.7875 -8.5 -1.425 -1.0875 -1.125 -1.2 -1.2375 -1.2 -1.125 -1.05 -1.275 -1.275 -9 -1.9875 -1.9125 -1.9125 -1.95 -1.9875 -1.95 -2.4375 -2.2875 -2.4375 -2.4 -9.5 0 0 0 0 0 63 LUẬN VĂN THẠC SĨ I3 HVTH: ĐOÀN TẤN VIỆT CƯỜNG Ngày đo 29/06/ 2018 03/07/ 2018 07/07/ 2018 11/07/ 2018 14/07/ 2018 17/07/ 2018 21/07/ 2018 24/07/ 2018 Chiều sâu đo Chu kỳ 31 Chu kỳ 32 Chu kỳ 33 Chu kỳ 34 Chu kỳ 35 Chu kỳ 36 Chu kỳ 37 Chu kỳ 38 -0.5 40.8375 42.1875 42 42.4875 41.0625 42.975 43.8 45.225 -1 36.225 37.6875 37.0875 37.5 36.225 37.8 38.475 39.6375 -1.5 31.6125 33.1875 32.175 32.5125 31.3875 32.625 33.15 34.0875 -2 22.0125 23.55 22.275 22.425 21.6 22.65 23.2875 24.1125 -2.5 13.95 15.5625 13.65 13.725 13.3125 14.0625 14.775 15.4125 -3 10.2 11.7 9.6 9.675 9.525 9.975 10.725 11.1375 -3.5 5.7 7.4625 4.725 4.95 5.175 5.2125 5.9625 6.375 -4 2.9625 4.8 2.025 2.0625 2.5875 2.7375 3.45 3.45 -4.5 2.925 4.8375 1.8375 1.95 2.6625 2.775 3.375 3.3375 -5 1.9875 3.6375 0.45 0.675 1.425 1.425 2.175 2.0625 -5.5 1.3875 3.0375 0.0375 0.2625 1.2 0.9375 1.6875 1.575 -6 3.5625 5.25 2.2125 2.6625 3.375 3.3375 3.975 3.9 -6.5 5.0625 5.925 3.825 4.35 4.9125 4.95 5.2875 5.4375 -7 5.7375 5.775 3.8625 4.575 4.875 4.8375 5.2125 5.2125 -7.5 4.2 4.0875 2.625 3.15 3.3375 3.3375 3.6 3.6 -8 1.5375 1.425 0.2625 0.75 0.825 0.825 1.0125 0.9375 -8.5 -1.0875 -1.3125 -1.7625 -1.4625 -1.4625 -1.5375 -1.3125 -1.425 -9 -2.475 -2.55 -2.7375 -2.55 -2.475 -2.475 -2.475 -2.5125 -9.5 0 0 0 0 64 LUẬN VĂN THẠC SĨ HVTH: ĐOÀN TẤN VIỆT CƯỜNG Phụ lục A.3: Kết quan trắc chuyển vị I4 (đơn vị mm) Chiều sâu đo 27/03/ 2018 Chu kỳ 30/03/ 2018 Chu kỳ 03/04/ 2018 Chu kỳ 05/04/ 2018 Chu kỳ 09/04/ 2018 Chu kỳ 13/04/ 2018 Chu kỳ 17/04/ 2018 Chu kỳ 20/04/ 2018 Chu kỳ 24/04/ 2018 Chu kỳ 27/04/ 2018 Chu kỳ 10 -0.5 2.30625 2.353125 -2.371875 -1.246875 6.159375 13.37813 113.48438 115.05938 115.8 -1 1.359375 1.275 -1.7625 -1.284375 1.678125 5.86875 107.16563 110.4375 111.79688 -1.5 0.375 0.24375 -1.153125 -1.340625 -2.803125 -1.6875 100.84688 105.81563 107.8125 -2 -0.20625 -0.525 -1.040625 -1.565625 -4.171875 -4.603125 94.78125 100.96875 102.92813 -2.5 -0.31875 -0.76875 -0.890625 -1.396875 -3.309375 -4.125 88.55625 95.184375 97.040625 -3 -0.29063 -0.75 -0.80625 -1.2 -1.7625 -2.315625 81.4125 87.20625 88.809375 -3.5 -0.1875 -0.609375 -0.665625 -0.946875 -0.796875 -0.7875 71.615625 75.99375 77.23125 -4 -0.17813 -0.553125 -0.609375 -0.815625 -0.609375 -0.290625 56.465625 59.446875 60.459375 -4.5 -0.15938 -0.45 -0.515625 -0.665625 -0.58125 -0.365625 33.628125 36.28125 37.471875 -5 -0.04688 -0.290625 -0.3 -0.4125 -0.39375 -0.5625 10.603125 13.44375 14.615625 -5.5 0.009375 -0.1875 -0.196875 -0.290625 -0.24375 -0.778125 -1.2375 0.43125 1.275 -6 0.01875 -0.16875 -0.121875 -0.225 -0.178125 -1.03125 -3.80625 -3.534375 -2.878125 -6.5 0.028125 -0.16875 -0.09375 -0.215625 -0.084375 -1.171875 -1.996875 -2.30625 -1.65 -7 0.0375 -0.15 -0.0375 -0.16875 -0.075 -0.9 0.7875 0.909375 1.584375 -7.5 0.046875 -0.140625 -0.009375 -0.15 -0.009375 -0.0375 2.615625 3.478125 4.05 -8 0.046875 -0.103125 0.084375 -0.09375 0.01875 1.0125 3.028125 4.0875 4.565625 -8.5 0.075 -0.075 0.13125 -0.028125 0.065625 0.778125 2.2125 3.075 3.440625 -9 0.075 -0.084375 0.09375 -0.0375 0.075 -0.196875 0.73125 1.471875 1.753125 -9.5 -0.00938 -0.15 0.0375 -0.13125 0.028125 -0.140625 0.1125 0.5625 0.73125 -10 -0.01875 -0.140625 0.009375 -0.1125 0.075 -0.20625 -0.0375 0.28125 0.4125 -10.5 -0.0375 -0.159375 -0.0375 -0.13125 0.0375 -0.309375 -0.121875 0.178125 0.3375 -11 -0.01875 -0.178125 -0.065625 -0.16875 -0.046875 -0.384375 -0.665625 -0.534375 -0.43125 -11.5 -0.00938 -0.1875 -0.075 -0.15 -0.046875 -0.375 -0.703125 -0.65625 -0.628125 -12 0.028125 -0.16875 -0.065625 -0.121875 -0.309375 -0.496875 -0.478125 -0.4875 -12.5 -0.01875 -0.215625 -0.1125 -0.196875 -0.065625 -0.28125 -0.4875 -0.35625 -0.43125 -13 -0.02813 -0.178125 -0.121875 -0.13125 -0.09375 -0.225 -0.328125 -0.1875 -0.24375 -13.5 0.05625 -0.103125 -0.009375 -0.0375 0.05625 -0.1125 -0.09375 -0.0375 -0.084375 -14 0.009375 -0.046875 -0.009375 0.01875 0.075 -0.05625 0.046875 0.0375 0.046875 -14.5 0.01875 -0.028125 0.028125 0.103125 -0.009375 0.09375 0.065625 0.1125 -15 0 0 0 0 0 Ngày đo I4 65 LUẬN VĂN THẠC SĨ HVTH: ĐOÀN TẤN VIỆT CƯỜNG Chiều sâu đo 02/05/ 2018 Chu kỳ 11 04/05/ 2018 Chu kỳ 12 09/05/ 2018 Chu kỳ 13 11/05/ 2018 Chu kỳ 14 13/05/ 2018 Chu kỳ 15 14/05/ 2018 Chu kỳ 16 16/05/ 2018 Chu kỳ 17 18/05/ 2018 Chu kỳ 18 21/05/ 2018 Chu kỳ 19 23/05/ 2018 Chu kỳ 20 -0.5 134.08125 129.34688 147.62813 146.775 149.91563 150.77813 150.80625 149.91563 150.3375 150.45 -1 128.29688 125.925 141.90938 142.875 144.01875 144.74063 144.50625 144.01875 150.8625 151.15313 -1.5 122.50313 122.41875 136.19063 138.96563 138.12188 138.70313 138.21563 138.12188 151.3875 151.86563 -2 116.56875 117.2625 129.9375 133.07813 131.85938 132.4125 131.9625 131.85938 149.23125 149.67188 -2.5 109.95 110.67188 122.75625 125.55938 124.71563 125.29688 124.9125 124.71563 145.09688 145.48125 -3 100.99688 101.61563 113.30625 115.92188 115.20938 115.78125 115.4625 115.20938 136.77188 137.04375 -3.5 88.003125 88.734375 99.5625 101.67188 101.0625 101.5875 101.35313 101.0625 120.36563 120.51563 -4 69.0375 69.590625 78.20625 79.228125 78.675 79.25625 79.275 78.675 92.30625 92.45625 -4.5 42.253125 42.825 47.2125 47.634375 47.175 47.915625 47.896875 47.175 52.725 52.89375 -5 16.171875 16.884375 17.840625 18.3375 17.953125 18.43125 18.553125 17.953125 18.75 19.125 -5.5 1.809375 2.484375 2.653125 2.971875 2.971875 2.953125 3.159375 2.971875 3.234375 3.628125 -6 -2.5875 -1.9125 -1.63125 -1.153125 -1.05 -1.228125 -0.965625 -1.05 -0.825 -0.346875 -6.5 -1.340625 -0.703125 0.196875 0.6375 0.75 0.571875 0.75 0.75 0.984375 1.4625 -7 1.9875 2.625 4.06875 4.584375 4.659375 4.4625 4.66875 4.659375 4.93125 5.38125 -7.5 4.575 5.175 6.965625 7.4625 7.5375 7.35 7.5375 7.5375 7.753125 8.15625 -8 5.053125 5.559375 7.021875 7.48125 7.490625 7.359375 7.471875 7.490625 7.7625 8.090625 -8.5 3.9375 4.359375 5.11875 5.41875 5.4375 5.334375 5.503125 5.4375 5.68125 -9 2.15625 2.578125 2.878125 3.159375 3.1875 3.084375 3.215625 3.1875 3.365625 3.665625 -9.5 0.99375 1.359375 1.55625 1.8375 1.89375 1.725 1.828125 1.89375 1.921875 2.221875 -10 0.65625 1.003125 1.059375 1.35 1.453125 1.25625 1.303125 1.453125 1.340625 1.659375 -10.5 0.478125 0.815625 0.825 1.115625 1.2 0.99375 1.0875 1.2 1.03125 1.340625 -11 -0.35625 -0.0375 -0.103125 0.178125 0.24375 0.028125 0.121875 0.24375 -0.01875 0.2625 -11.5 -0.628125 -0.35625 -0.496875 -0.1875 -0.084375 -0.309375 -0.253125 -0.084375 -0.384375 -0.103125 -12 -0.440625 -0.215625 -0.3 -0.046875 0.009375 -0.215625 -0.05625 0.009375 -0.31875 -0.046875 -12.5 -0.403125 -0.225 -0.328125 -0.046875 0.05625 -0.159375 -0.028125 0.05625 -0.253125 0.01875 -13 -0.20625 -0.05625 -0.140625 0.103125 0.225 0.0375 0.13125 0.225 -0.065625 0.20625 -13.5 -0.01875 0.05625 0.05625 0.140625 0.215625 0.159375 0.196875 0.215625 0.196875 -14 0.140625 0.140625 0.15 0.196875 0.31875 0.24375 0.225 0.31875 0.084375 0.20625 -14.5 0.178125 0.1875 0.16875 0.20625 0.271875 0.2625 0.234375 0.271875 0.103125 0.16875 -15 0 0 0 0 0 Ngày đo I4 66 LUẬN VĂN THẠC SĨ HVTH: ĐOÀN TẤN VIỆT CƯỜNG Chiều sâu đo 25/05/ 2018 Chu kỳ 21 30/05/ 2018 Chu kỳ 22 02/06/ 2018 Chu kỳ 23 05/06/ 2018 Chu kỳ 24 08/06/ 2018 Chu kỳ 25 12/06/ 2018 Chu kỳ 26 15/06/ 2018 Chu kỳ 27 20/06/ 2018 Chu kỳ 28 22/06/ 2018 Chu kỳ 29 26/06/ 2018 Chu kỳ 30 -0.5 153.34688 159.83438 152.69063 150.39375 149.59688 154.8375 152.83125 152.5875 155.57813 158.99063 -1 152.46563 157.275 153.12188 151.80938 151.40625 154.24688 152.80313 153.06563 154.82813 157.11563 -1.5 151.58438 154.60313 153.60938 153.225 153.21563 153.5625 152.89688 153.51563 153.9375 155.25 -2 149.40938 150.51563 151.82813 152.175 152.14688 152.27813 151.8 152.87813 152.39063 152.7 -2.5 145.3125 145.575 147.72188 148.45313 148.41563 148.575 148.2 149.50313 148.92188 149.05313 -3 136.89375 136.71563 138.64688 139.35938 139.3125 139.425 139.125 140.33438 139.94063 140.16563 -3.5 120.61875 120.525 121.3125 121.86563 121.96875 121.70625 121.52813 122.5875 122.19375 122.625 -4 92.821875 92.5875 92.446875 92.7375 93.01875 92.446875 92.475 93.309375 92.878125 93.459375 -4.5 53.446875 53.165625 53.08125 53.38125 53.728125 53.11875 53.23125 54.05625 53.634375 54.1875 -5 19.3875 19.425 19.471875 19.696875 19.96875 19.55625 19.509375 20.30625 19.978125 20.446875 -5.5 3.609375 3.665625 3.95625 4.05 4.21875 3.984375 3.740625 4.528125 4.246875 4.621875 -6 -0.571875 -0.54375 -0.028125 0.01875 0.075 -0.009375 -0.346875 0.45 0.1875 0.54375 -6.5 1.18125 1.06875 1.621875 1.65 1.70625 1.659375 1.303125 2.053125 1.809375 2.128125 -7 5.071875 4.89375 5.41875 5.4375 5.503125 5.49375 5.11875 5.821875 5.5875 5.94375 -7.5 7.921875 7.78125 8.23125 8.23125 8.325 8.30625 7.93125 8.634375 8.4 8.746875 -8 7.9875 7.8375 8.203125 8.25 8.353125 8.34375 8.015625 8.68125 8.4375 8.765625 -8.5 5.971875 5.784375 6.15 6.196875 6.328125 6.3 6.628125 6.43125 6.740625 -9 3.6375 3.46875 3.84375 3.834375 3.975 3.95625 3.646875 4.265625 4.06875 4.378125 -9.5 2.146875 2.08125 2.409375 2.38125 2.446875 2.465625 2.165625 2.746875 2.53125 2.85 -10 1.59375 1.5375 1.8375 1.78125 1.8375 1.81875 1.584375 2.128125 1.9125 2.2125 -10.5 1.321875 1.2 1.509375 1.4625 1.51875 1.471875 1.275 1.771875 1.6125 1.903125 -11 0.24375 0.140625 0.4125 0.346875 0.375 0.346875 0.159375 0.590625 0.440625 0.721875 -11.5 -0.159375 -0.2625 0.01875 -0.065625 -0.05625 -0.084375 -0.271875 0.103125 -0.05625 0.215625 -12 -0.046875 -0.1875 0.065625 -0.01875 0.01875 -0.009375 -0.1875 0.140625 0.009375 0.2625 -12.5 -0.09375 -0.1125 0.103125 -0.009375 0 -0.1875 0.13125 0.01875 0.225 -13 0.09375 0.121875 0.309375 0.178125 0.159375 0.196875 -0.009375 0.2625 0.140625 0.375 -13.5 0.16875 0.159375 0.253125 0.1875 0.24375 0.225 0.075 0.253125 0.20625 0.384375 -14 0.20625 0.159375 0.28125 0.15 0.271875 0.215625 0.140625 0.253125 0.234375 0.35625 -14.5 0.215625 0.196875 0.253125 0.159375 0.253125 0.16875 0.215625 0.2625 0.225 0.271875 -15 0 0 0 0 0 Ngày đo I4 67 LUẬN VĂN THẠC SĨ HVTH: ĐOÀN TẤN VIỆT CƯỜNG Chiều sâu đo 29/06/ 2018 Chu kỳ 31 03/07/ 2018 Chu kỳ 32 07/07/ 2018 Chu kỳ 33 11/07/ 2018 Chu kỳ 34 14/07/ 2018 Chu kỳ 35 17/07/ 2018 Chu kỳ 36 21/07/ 2018 Chu kỳ 37 24/07/ 2018 Chu kỳ 38 28/07/ 2018 Chu kỳ 39 -0.5 220.54688 220.56563 220.26563 219.8625 195.35625 177.675 185.625 188.86875 188.65313 -1 193.10625 193.28438 192.975 192.54375 180.31875 169.18125 174.73125 176.72813 176.42813 -1.5 165.6375 166.00313 165.68438 165.225 165.29063 160.32188 163.875 164.5875 164.20313 -2 148.79063 149.30625 149.12813 148.575 152.84063 151.9875 152.98125 153.09375 152.87813 -2.5 142.47188 142.9875 142.82813 142.125 145.54688 144.9 144.7125 144.87188 144.6375 -3 136.69688 137.14688 136.9875 136.20938 137.48438 138.3 138.45 138.6 138.38438 -3.5 122.19375 122.42813 122.30625 121.56563 121.59375 123 123.1125 123.23438 122.98125 -4 93.328125 93.328125 93.290625 92.784375 92.840625 93.271875 93.421875 93.525 93.2625 -4.5 54.13125 54.01875 54 53.634375 53.653125 53.840625 53.8875 53.9625 53.728125 -5 20.30625 20.428125 20.390625 20.034375 20.0625 20.34375 20.334375 20.475 20.23125 -5.5 4.471875 4.790625 4.678125 4.36875 4.471875 4.715625 4.659375 4.74375 4.66875 -6 0.328125 0.69375 0.6375 0.346875 0.39375 0.65625 0.6 0.73125 0.646875 -6.5 1.875 2.259375 2.23125 1.9125 1.89375 2.15625 2.109375 2.2125 2.175 -7 5.653125 6.0375 5.990625 5.709375 5.709375 5.9625 5.925 5.98125 -7.5 8.49375 8.86875 8.83125 8.55 8.521875 8.765625 8.728125 8.75625 8.803125 -8 8.596875 8.878125 8.85 8.60625 8.55 8.765625 8.728125 8.71875 8.775 -8.5 6.6 6.834375 6.825 6.609375 6.54375 6.721875 6.69375 6.675 6.7125 -9 4.275 4.48125 4.490625 4.284375 4.2375 4.396875 4.36875 4.35 4.359375 -9.5 2.7375 2.94375 2.9625 2.728125 2.709375 2.840625 2.83125 2.83125 2.83125 -10 2.090625 2.34375 2.371875 2.11875 2.071875 2.175 2.203125 2.221875 2.175 -10.5 1.809375 2.04375 2.090625 1.81875 1.8 1.85625 1.9125 1.921875 1.85625 -11 0.61875 0.871875 0.909375 0.61875 0.61875 0.646875 0.7125 0.721875 0.65625 -11.5 0.103125 0.3375 0.365625 0.084375 0.075 0.084375 0.15 0.178125 0.09375 -12 0.159375 0.35625 0.403125 0.09375 0.103125 0.1125 0.196875 0.20625 0.1125 -12.5 0.103125 0.28125 0.31875 0.0375 0.075 0.028125 0.140625 0.178125 0.084375 -13 0.234375 0.384375 0.375 0.20625 0.253125 0.1875 0.290625 0.365625 0.2625 -13.5 0.24375 0.384375 0.4125 0.253125 0.3 0.20625 0.31875 0.384375 0.28125 -14 0.2625 0.346875 0.375 0.290625 0.3 0.196875 0.31875 0.365625 0.290625 -14.5 0.20625 0.271875 0.328125 0.253125 0.20625 0.13125 0.253125 0.271875 0.28125 -15 0 0 0 0 Ngày đo I4 68 LUẬN VĂN THẠC SĨ HVTH: ĐOÀN TẤN VIỆT CƯỜNG Phụ lục B.1: Code Matlab “So sánh với kết Veletsos Younan (1994)” function f1 = baitoan1 nu = 0.3; psisig = 2/sqrt((1-nu)*(2-nu)); psie = sqrt((2-nu)/(1-nu)); delta = 0.1; eta = 1; omega_omega1 = [0 0.9 4]; xi = linspace(0,10,101); vectK = linspace(0,100,101); vectN = 2*vectK+1; vect1N = 1./(2*vectK+1); matNN = ones(length(omega_omega1),1)*vectN; phiN = omega_omega1'*vect1N; LambdaN = 1-(phiN.^2/(1+1i*delta)); alphaN = pi/2/psie*( sqrt(LambdaN).*matNN ); for jomega = 1:length(omega_omega1) for jxi = 1:length(xi) vectSinN= sin(pi/2*eta*vectN); sigmaN(jomega,jxi) = -(8*psisig/pi^2) * (1+1i*delta)*sum (vectSinN.*(vect1N.^2) / (1phiN(jomega,:).^2+1i*delta).*((sqrt(LambdaN(jomega,:))).*e xp(-alphaN(jomega,:).*xi(jxi)))); end end f1 = sigmaN; 69 LUẬN VĂN THẠC SĨ HVTH: ĐOÀN TẤN VIỆT CƯỜNG Phụ lục B.2: Code Matlab “Ảnh hưởng khoảng cách từ tải trọng đến tường vây” function f2= baitoan2(E,nu,rho,delta,H,L) G = E/(2*(1+nu)); psisig = 2/sqrt((1-nu)*(2-nu)); psie = sqrt((2-nu)/(1-nu)); omega1 = pi/2/H*sqrt(G/rho); xi = linspace(0,L/H,101); eta = 1; vitkmh = linspace(0,100,101); vitms = vitkmh/3.6; f = 0.45*vitkmh; omega = 2*pi*f; vectK = linspace(0,100,101); vectN = 2*vectK+1; vect1N = 1./(2*vectK+1); matNN = ones(length(omega),1)*vectN; phiN = omega'*vect1N./omega1; LambdaN = 1-(phiN.^2/(1+1i*delta)); alphaN = pi/2/psie*( sqrt(LambdaN).*matNN ); vectSinN= sin(pi/2*eta'*vectN); for jomega = 1:length(omega) for jxi = 1:length(xi) sigma(jomega,jxi) = -(8*psisig/pi^2) *(rho*H*vitms(jomega)*omega(jomega))* (1+1i*delta)* sum ( vectSinN.*(vect1N.^2) / (1phiN(jomega,:).^2+1i*delta).* ( (sqrt(LambdaN(jomega,:))).* exp(-alphaN(jomega,:).*xi(jxi)) ) ) ; end end f2 = sigma; 70 LUẬN VĂN THẠC SĨ HVTH: ĐOÀN TẤN VIỆT CƯỜNG Phụ lục B.3: Code Matlab “Áp lực động tải trọng xe lên tường vây” function f3= baitoan3(E,nu,rho,delta,H,L) G = E/(2*(1+nu)); psisig = 2/sqrt((1-nu)*(2-nu)); psie = sqrt((2-nu)/(1-nu)); omega1 = pi/2/H*sqrt(G/rho); xi = L/H; eta = linspace(0,1,101); vitkmh = linspace(0,100,101); f = 0.45*vitkmh; omega = 2*pi*f; vitms = vitkmh/3.6; vectK = linspace(0,100,101); vectN = 2*vectK+1; vect1N = 1./(2*vectK+1); matNN = ones(length(omega),1)*vectN; phiN = omega'*vect1N./omega1; LambdaN = 1-(phiN.^2/(1+1i*delta)); alphaN = pi/2/psie*( sqrt(LambdaN).*matNN ); for jomega = 1:length(omega) for jeta = 1:length(eta) vectSinN= sin(pi/2*eta(jeta)*vectN); sigma(jomega,jeta) = -(8*psisig/pi^2) *(rho*H*vitms(jomega)*omega(jomega))* (1+1i*delta)* sum ( vectSinN.*(vect1N.^2) / (1phiN(jomega,:).^2+1i*delta).* ( (sqrt(LambdaN(jomega,:))).* exp(-alphaN(jomega,:).*xi) ) ) ; end end f3 = sigma; 71 LUẬN VĂN THẠC SĨ HVTH: ĐOÀN TẤN VIỆT CƯỜNG Phụ lục B.4: Code Matlab “Lực hệ giằng tải trọng xe tác động lên tường vây- trường hợp hệ giằng” function [F1]= baitoan4(E,nu,rho,delta,H,L,H1) G = E/(2*(1+nu)); psisig = 2/sqrt((1-nu)*(2-nu)); psie = sqrt((2-nu)/(1-nu)); omega1 = pi/2/H*sqrt(G/rho); xi = L/H; eta = linspace(0,H/H,101); vitkmh = linspace(0,100,101); f = 0.45*vitkmh; omega = 2*pi*f; vitms = vitkmh/3.6; vectK = linspace(0,100,101); vectN = 2*vectK+1; vect1N = 1./(2*vectK+1); matNN = ones(length(omega),1)*vectN; phiN = omega'*vect1N./omega1; LambdaN = 1-(phiN.^2/(1+1i*delta)); alphaN = pi/2/psie*( sqrt(LambdaN).*matNN ); for jomega = 1:length(omega) M(jomega)=0; for jeta = 1:length(eta) vectSinN= sin(pi/2*eta(jeta)*vectN); sigma(jomega,jeta) = -(8*psisig/pi^2) *(rho*H*vitms(jomega)*omega(jomega)/10)* (1+1i*delta)* sum ( vectSinN.*(vect1N.^2) / (1phiN(jomega,:).^2+1i*delta).* ( (sqrt(LambdaN(jomega,:))).* exp(-alphaN(jomega,:).*xi) ) ) ; M(jomega)=M(jomega)+sigma(jomega,jeta)*H/100.*eta(jeta)*H; end F1(jomega)=abs(M(jomega)/H1)/1e3; end 72 LUẬN VĂN THẠC SĨ HVTH: ĐOÀN TẤN VIỆT CƯỜNG Phụ lục B.5: Code Matlab “Lực hệ giằng tải trọng xe tác động lên tường vây- trường hợp hệ giằng” function [F1,F2]= baitoan5(E,nu,rho,delta,H,L,H1,H2) G = E/(2*(1+nu)); psisig = 2/sqrt((1-nu)*(2-nu)); psie = sqrt((2-nu)/(1-nu)); omega1 = pi/2/H*sqrt(G/rho); xi = L/H; eta = linspace(0,H/H,101); vitkmh = linspace(0,100,101); f = 0.45*vitkmh; omega = 2*pi*f; vitms = vitkmh/3.6; vectK = linspace(0,100,101); vectN = 2*vectK+1; vect1N = 1./(2*vectK+1); matNN = ones(length(omega),1)*vectN; phiN = omega'*vect1N./omega1; LambdaN = 1-(phiN.^2/(1+1i*delta)); alphaN = pi/2/psie*( sqrt(LambdaN).*matNN ); for jomega = 1:length(omega) F3(jomega)=0; Fq(jomega)=0; M(jomega)=0; for jeta = 1:length(eta) vectSinN= sin(pi/2*eta(jeta)*vectN); sigma(jomega,jeta) = -(8*psisig/pi^2) *(rho*H*vitms(jomega)*omega(jomega)/10)* (1+1i*delta)* sum ( vectSinN.*(vect1N.^2) / (1phiN(jomega,:).^2+1i*delta).* ( (sqrt(LambdaN(jomega,:))).* exp(-alphaN(jomega,:).*xi) ) ) ; if jeta