Phương pháp kết hợp giữa phần tử biên và phần tử hữu hạn trong phân tích tấm nổi trực hướng chịu tải trọng di động

ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA NGUYỄN PHƯỢNG KIỀU PHƯƠNG PHÁP KẾT HỢP GIỮA PHẦN TỬ BIÊN VÀ PHẦN TỬ HỮU HẠN TRONG PHÂN TÍCH TẤM NỔI TRỰC HƯỚNG CHỊU TẢI TRỌNG DI ĐỘNG Chun ngành: Kỹ Thuật Xây Dựng Cơng Trình Dân dụng Và Công nghiệp Mã số ngành: 60580208 LUẬN VĂN THẠC SĨ Tp.HCM, năm 2018 CƠNG TRÌNH ĐƯỢC HỒN THÀNH TẠI TRƯỜNG ĐẠI HỌC BÁCH KHOA ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH Cán hướng dẫn khoa học: Cán hướng dẫn: PGS.TS LƯƠNG VĂN HẢI Cán chấm nhận xét 1: TS Nguyễn Hồng Ân Cán chấm nhận xét 2: PGS TS Nguyễn Trọng Phước Luận văn thạc sĩ bảo vệ Trường Đại học Bách Khoa, ĐHQG Tp HCM, ngày 07 tháng 02 năm 2018 Thành phần Hội đồng đánh giá Luận văn thạc sĩ gồm: PGS TS Bùi Công Thành - Chủ tịch Hội đồng PGS TS Ngô Hữu Cường - Thư ký TS Nguyễn Hồng Ân - Ủy viên (Phản biện 1) PGS TS Nguyễn Trọng Phước - Ủy viên (Phản biện 2) PGS TS Nguyễn Văn Hiếu Ủy viên 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ự Do - Hạnh Phúc NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên: NGUYỄN PHƯỢNG KIỀU Ngày, tháng, năm sinh: 03/12/1993 MSHV: 1570647 Nơi sinh: Bạc Liêu Chuyên ngành: Kỹ Thuật Xây Dựng Dân Dụng Công Nghiệp Mã số: 60580208 I TÊN ĐỀ TÀI: PHƯƠNG PHÁP KẾT HỢP GIỮA PHẦN TỬ BIÊN VÀ PHẦN TỬ HỮU HẠN TRONG PHÂN TÍCH TẤM NỔI TRỰC HƯỚNG CHỊU TẢI TRỌNG DI ĐỘNG II NHIỆM VỤ VÀ NỘI DUNG Trình bày sở lý thuyết trực hướng, chất lỏng lý tưởng, phương pháp phần tử hữu hạn phương pháp phần tử biên Phát triển thuật toán giải hệ phương trình tương tác chất lỏng miền thời gian xây dựng chương trình Matlab mơ ứng xử hệ tấmnước Kiểm chứng chương trình thực toán khảo sát để đánh giá ảnh hưởng tính trực hướng lên đại lượng động lực học đặc trưng hệ Đưa kết luận kiến nghị từ kết nghiên cứu số III NGÀY GIAO NHIỆM VỤ : 04/09/2017 IV NGÀY HOÀN THÀNH NHIỆM VỤ : 08/01/2018 V HỌ VÀ TÊN CÁN BỘ HƯỚNG DẪN: PGS TS LƯƠNG VĂN HẢI CÁN BỘ HƯỚNG DẪN (Họ tên chữ ký) Tp HCM, ngày tháng năm 2018 BAN QUẢN LÝ CHUYÊN NGÀNH (Họ tên chữ ký) Lương Văn Hải TRƯỞNG KHOA KỸ THUẬT XÂY DỰNG (Họ tên chữ ký) ii LỜI CẢM ƠN Luận văn thạc sĩ Xây dựng cơng trình dân dụng cơng nghiệp nằm hệ thống luận cuối khóa nhằm trang bị cho học viên cao học khả tự nghiên cứu, biết cách giải vấn đề cụ thể đặt thực tế xây dựng v.v… Đó trách nhiệm niềm tự hào học viên cao học Để hoàn thành luận văn này, cố gắng nỗ lực thân, nhận giúp đỡ nhiều từ tập thể cá nhân Tơi xin ghi nhận tỏ lịng biết ơn tới tập thể cá nhân dành cho tơi giúp đỡ q báu Đầu tiên tơi xin bày tỏ lòng biết ơn sâu sắc đến Thầy PGS.TS Lương Văn Hải Thầy TS Trần Minh Thi Các thầy đưa gợi ý để hình thành nên ý tưởng đề tài, góp ý cho nhiều cách nhận định đắn vấn đề nghiên cứu, cách tiếp cận nghiên cứu hiệu Tôi xin chân thành cảm ơn quý Thầy Cô Khoa Kỹ thuật Xây dựng, trường Đại học Bách Khoa Tp.HCM truyền dạy kiến thức quý giá cho tơi, kiến thức khơng thể thiếu đường nghiên cứu khoa học nghiệp sau Tôi xin gửi lời cảm ơn đến anh NCS Nguyễn Xuân Vũ giúp đỡ tơi nhiều q trình thực luận văn Luận văn thạc sĩ hoàn thành thời gian quy định với nỗ lực thân, nhiên khơng thể khơng có thiếu sót Kính mong q Thầy Cơ dẫn thêm để tơi bổ sung kiến thức hồn thiện thân Xin trân trọng cảm ơn Tp HCM, ngày 08 tháng 01 năm 2018 NGUYỄN PHƯỢNG KIỀU iii TÓM TẮT LUẬN VĂN THẠC SĨ Đại dương bao phủ 70% Trái Đất nơi làm giảm nhu cầu cấp bách vấn đề đất đai ven biển Gần 50% giới cơng nghiệp hóa sống vịng km bờ biển, nhu cầu tài nguyên đất không gian bắt đầu tiến tới giai đoạn quan trọng dân số giới tiếp tục tăng với tốc độ đáng báo động Cần có phát triển bền vững thân thiện với môi trường Một sáng kiến thân thiện với môi trường xuất thời gian gần khái niệm cơng trình siêu lớn (VLFS) - công nghệ cho phép tạo vùng đất nhân tạo từ biển mà không phá hủy môi trường biển, không làm ô nhiễm nguồn nước ven biển khơng thay đổi dịng chảy thủy triều dịng chảy tự nhiên Cho đến nhiều phương pháp số phát triển để phân tích VLFS mơi trường sóng biển, luận văn sử dụng phương pháp kết hợp phần tử biên (BEM) phương pháp phần tử hữu hạn (FEM) để phân tích kết cấu trực hướng chịu tải trọng di động để giải toán kết cấu siêu lớn Chương phần giới thiệu tổng quan kết cấu siêu lớn, lịch phát triển, ưu điểm ứng dụng kết cấu siêu lớn thực tế, mục tiêu hướng nghiên cứu đề tài Trong chương trình bày mơ hình kết cấu chất lỏng với điều kiện biên thiết lập q trình tính tốn, trình bày sở lý thuyết trực hướng, chất lỏng lý tưởng, phương pháp phần tử hữu hạn phương pháp phần tử biên sử dụng q trình tính tốn Chương sử dụng thuật tốn giải hệ phương trình tương tác chất lỏng miền thời gian xây dựng chương trình Matlab mơ ứng xử hệ tấmnước, tiến hành phân tích với giá trị độ cứng uốn phương thay đổi với thay đổi vận tốc, khối lượng tải trọng di động, chiều dày chiều rộng Từ kết chương 3, nhận xét kiến nghị trình bày chương cuối danh mục tài liệu tham khảo sử dụng luận văn iv LỜI CAM ĐOAN Tôi xin cam đoan cơng việc tơi thực hướng dẫn Thầy Lương Văn Hải Các kết Luận văn thật chưa công bố nghiên cứu khác Tôi xin chịu trách nhiệm cơng việc thực Tp HCM, ngày 08 tháng 01 năm 2018 NGUYỄN PHƯỢNG KIỀU v MỤC LỤC NHIỆM VỤ LUẬN VĂN THẠC SĨ i LỜI CẢM ƠN .ii TÓM TẮT LUẬN VĂN THẠC SĨ iii LỜI CAM ĐOAN iv MỤC LỤC v DANH MỤC CÁC HÌNH VẼ vi DANH MỤC CÁC BẢNG BIỂU .vii MỘT SỐ KÝ HIỆU VIẾT TẮT viii CHƯƠNG TỔNG QUAN 1.1 Giới thiệu 1.1.1 Lịch sử phát triển 1.1.2 Ưu điểm 1.2 Ứng dụng 1.3 Tổng quan nghiên cứu cần thiết tiến hành nghiên cứu 1.3.1 Lý thuyết Hydroelastic .6 1.3.2 Các phương pháp phân tích 1.4 Mục tiêu hướng nghiên cứu 10 1.4.1 Mục tiêu 10 1.4.2 Hướng nghiên cứu 10 1.5 Cấu trúc luận văn 10 CHƯƠNG CƠ SỞ LÝ THUYẾT 11 2.1 Mơ hình cấu trúc chất lỏng điều kiện biên 11 2.1.1 Mô hình cấu trúc chất lỏng .11 2.1.2 Điều kiện biên 12 2.2 Lý thuyết trực hướng 13 2.3 Phương pháp phần tử hữu hạn áp dụng cho trực hướng 15 2.4 Phương pháp phần tử biên 17 v 2.4.1 Nghiệm phương trình Laplace 17 2.4.2 Phương trình động lực học kết cấu 20 2.5 Lưu đồ thuật toán 22 CHƯƠNG VÍ DỤ SỐ 23 3.1 Bài toán 1: Kiếm chứng chương trình phân tích hydroelastic đẳng hướng so với kết tác giả Ismail 24 3.2 Kiểm chứng chương trình Matlab phân tích trực hướng 25 3.2.1 Bài tốn 2: Phân tích ứng xử trực hướng chịu tác dụng tải trọng tĩnh 25 3.2.2 Bài tốn 3: Phân tích dao động tự trực hướng .28 3.3 Bài toán 4: Phân tích ứng xử động trực hướng với giá trị độ cứng uốn theo phương thay đổi ứng với thay đổi vận tốc v 32 3.3.1 Khảo sát chuyển vị lớn thay đổi vận tốc tải trọng di động 32 3.3.2 Khảo sát chuyển vị theo phương x y thay đổi vận tốc tải trọng di động 35 3.3.3 Khảo sát chuyển vị điểm theo thời gian thay đổi vận tốc tải trọng di động .41 3.4 Bài tốn 5: Phân tích ứng xử động trực hướng với giá trị độ cứng uốn theo phương thay đổi ứng với thay đổi tải trọng tác dụng P 43 3.5 Bài tốn 6: Phân tích ứng xử động trực hướng với giá trị độ cứng uốn theo phương thay đổi ứng với thay đổi chiều dày hs 46 3.6 Bài tốn 7: Phân tích ứng xử động trực hướng với giá trị độ cứng uốn theo phương thay đổi ứng với thay đổi kích thước L B 48 CHƯƠNG KẾT LUẬN VÀ KIẾN NGHỊ 53 4.1 Kết luận 53 v 4.2 Kiến nghị 54 TÀI LIỆU THAM KHẢO 55 vi DANH MỤC CÁC HÌNH VẼ Hình 1 Các thành phần hệ thống kết cấu Hình Mega-Float vịnh Tokyo, Nhật Bản Hình Sân bay Quốc tế Tokyo, (Haneda) Hình Bến tàu (Valdez, Alaska) Hình Bến container di động Hình Căn cứu hộ khẩn cấp vịnh Tokyo vịnh Osaka Hình Bãi đáp trực thăng nổi, Vancouver, Canada Hình Căn lưu trữ dầu Shirashima, Japan Hình Căn lưu trữ dầu Kamigoto, tỉnh Nagasaki , Japan Hình 10 Phản ứng toàn phần tác dụng tải trọng tĩnh Hình 11 Sơ đồ phản ứng toàn phần kết cấu Hình Hệ cấu trúc chất lỏng hệ tọa độ 11 Hình 2 Kết cấu siêu lớn VLFS 13 Hình Phần tử chữ nhật với điểm nút 16 Hình Hàm Green G  P, Q  17 Hình Miền khảo sát, mặt biên tròn mặt biên tròn S ,   r 18 Hình Mơ hình kết cấu VLFS 23 Hình Mơ hình kích thước 9, 75  1, 95m 24 Hình Kết chuyển vị Ismail luận văn 25 Hình 4 Chuyển vị vị trí đặt lực ứng với lưới chia phần tử 26 Hình Chuyển vị vị trí đặt lực với lưới chia phần tử 30  30 27 Hình Tần số dao động tự nhiên mode dao động thứ với thay đổi module đàn hồi 28 Hình So sánh tần số dao động tự nhiên mode dao động thứ với phần mềm SAP2000 30 Hình Tần số dao động tự nhiên sáu mode dao động 31 Chương 3: Ví dụ số 3.4 Bài tốn 5: Phân tích ứng xử động trực hướng với giá trị độ cứng uốn theo phương thay đổi ứng với thay đổi tải trọng tác dụng P Thông số kích thước tấm, thơng số vật liệu cho Bảng Trong toán này, khảo sát ảnh hưởng lực di chuyển với vận tốc v  9, m / s đến ứng xử động kết cấu xem xét trường hợp tải trọng di động P  1, kN ; P  2, kN ; P  3, kN ; P  5, kN P  7, kN Chuyển vị dọc theo trục x  y  bs /  thể Hình 19, tương ứng trọng tâm xe vị trí x  30, m a Lực di chuyển P  1, kN b Lực di chuyển P  2, kN 43 Chương 3: Ví dụ số c Lực di chuyển P  3, kN d Lực di chuyển P  5, kN e Lực di chuyển P  7, kN Hình 19 Chuyển vị theo phương x  y  bs /  thay đổi giá trị tải trọng 44 Chương 3: Ví dụ số Bảng 10 Chuyển vị lớn theo phương x  y  bs /  thay đổi giá trị tải trọng Tải trọng (kN) Module đàn hồi E1  0, E2 E1  0, E2 E1  1, E2 E1  1, 5E2 E1  2, E2 P  1, 5kN -6,04 -5,81 -5,40 -4,74 -3,94 P  2, 5kN -10,06 -9,69 -9,01 -7,91 -6,56 P  3, 5kN -14,09 -13,56 -12,61 -11,07 -9,19 P  5, 0kN -20,13 -19,37 -18,02 -15,82 -13,13 P  7, 5kN -30,19 -29,06 -27,02 -23,72 -19,69 Tỷ lệ trường hợp so với P  1, 5kN P  1, 5kN 1,00 1,00 1,00 1,00 1,00 P  2, 5kN 1,67 1,67 1,67 1,67 1,66 P  3, 5kN 2,33 2,33 2,34 2,34 2,33 P  5, 0kN 3,33 3,33 3,34 3,34 3,33 P  7, 5kN 5,00 5,00 5,00 5,00 5,00 Từ kết cho Hình 19 Bảng 10, module đàn hồi E1 tăng dần chuyển vị giảm dần, cụ thể trường hợp P  7, kN E1  0, E2 chuyển vị 30,19 mm , module đàn hồi E1  1, E2 chuyển vị 27, 02 mm (giảm 10, 5% ) module đàn hồi E1  2, E2 chuyển vị 19, 69 mm (giảm 34, 8% ) Từ kết trên, ta nhận thấy giá trị lực di chuyển P tăng dần chuyển vị tăng dần, cụ thể P tăng 3, lần chuyển vị tăng 3, lần P tăng lần chuyển vị tăng lần Tỷ lệ so với P  1, kN không thay đổi module đàn hồi thay đổi 45 Chương 3: Ví dụ số 3.5 Bài tốn 6: Phân tích ứng xử động trực hướng với giá trị độ cứng uốn theo phương thay đổi ứng với thay đổi chiều dày hs Bài toán khảo sát ảnh hưởng chiều dày đến ứng xử động lực học kết cấu xem xét trường hợp hs1  0,17 m , hs  0, 25 m , hs  0, m hs  1, m Với tải trọng di động P  kN di chuyển với vận tốc v  9, m / s Chuyển vị dọc theo trục lực di chuyển thể Hình 20 tương ứng với trường hợp module đàn hồi, trọng tâm xe vị trí x  30, m a Chuyển vị theo phương x mơ hình 3D hs  0,17 m b Chuyển vị theo phương x mô hình 3D hs  0, 25 m c Chuyển vị theo phương x mơ hình 3D hs  0, m 46 Chương 3: Ví dụ số d Chuyển vị theo phương x mơ hình 3D hs  1, m Hình 20 Chuyển vị theo phương x chiều dày thay đổi Bảng 11 Chuyển vị lớn theo phương x  y  bs /  chiều dày thay đổi Module đàn hồi Chiều dày (m) hs  0,17 m hs  0, 25 m hs  0, 50 m hs  1, 00 m E1  0, E2 -20,13 -12,37 -4,59 -1,59 E1  0, E2 -19,37 -11,46 -4,42 -1,53 E1  1, E2 -18,02 -10,51 -4,09 -1,43 E1  1, 5E2 -15,82 -9,47 -3,16 -1,19 E1  2, E2 -13,13 -7,89 -2,69 -1,00 So sánh với trường hợp E1  0, E2 E1  0, E2 X X X X E1  0, E2 -3,7% -7,4% -3,8% -3,9% E1  1, E2 -10,5% -15,0% -10,9% -10,1% E1  1, 5E2 -21,4% -23,4% -31,3% -25,6% E1  2, E2 -34,8% -36,2% -41,4% -37,3% 47 Chương 3: Ví dụ số Cũng giống toán khảo sát trước, module đàn hồi tăng dần chuyển vị giảm dần chiều dày nhỏ số đỉnh sóng nhiều Cụ thể trường hợp E1  0, E2 , chiều dày hs  0,17 m số đỉnh sóng 4, hs  0, 25 m số đỉnh sóng hs  1, m số đỉnh sóng Cũng trường hợp chuyển vị lớn hs  1, m  1, 59 mm  nhỏ 12, 66 lần so với hs  0,17 m  20,13 mm  Hình 21 thể so sánh chuyển vị theo phương x tương ứng với trường hợp chiều dày hs E1  0, E2 Hình 21 Chuyển vị theo phương x  y  bs /  chiều dày thay đổi trường hợp E1  0, E2 3.6 Bài toán 7: Phân tích ứng xử động trực hướng với giá trị độ cứng uốn theo phương thay đổi ứng với thay đổi kích thước L B Bài toán khảo sát ảnh hưởng tỷ lệ chiều dài chiều rộng đến ứng xử động lực học kết cấu Chiều dài kết cấu Ls  60 m chiều rộng thay đổi bs1  60 m , bs  48 m , bs  36 m , bs  24 m bs  12 m Với tải trọng di động P  2, kN di chuyển với vận tốc v  9, m / s , chuyển vị lớn thể Hình 22 tương ứng với trường hợp module đàn hồi phương thay đổi trọng tậm xe vị trí x  30, m Theo kết khảo sát 48 Chương 3: Ví dụ số chuyển vị lớn đạt giá trị nhỏ tỷ lệ Ls / bs 1, 67 2, 50 tương ứng với bs  36 m đến bs  24 m Kết giúp cho việc chọn lựa kích thước phù hợp để chuyển vị kết cấu nhỏ Hình 22 Chuyển vị lớn thay đổi chiều rộng a Chuyển vị theo phương x mô hình 3D Ls  60 m , bs  60 m b Chuyển vị theo phương x mơ hình 3D Ls  60 m , bs  48 m 49 Chương 3: Ví dụ số c Chuyển vị theo phương x mô hình 3D Ls  60 m , bs  36 m d Chuyển vị theo phương x mơ hình 3D Ls  60 m , bs  24 m e Chuyển vị theo phương x mơ hình 3D Ls  60 m , bs  12 m Hình 23 Chuyển vị theo phương x  y  bs /  thay đổi chiều rộng 50 Chương 3: Ví dụ số Bảng 12 Chuyển vị lớn theo phương x  y  bs /  thay đổi chiều rộng Module đàn hồi Chiều rộng kết cấu bs  m  bs  60 m bs  48 m bs  36 m bs  24 m bs  12 m E1  0, E2 -10,06 -9,40 -8,48 -7,69 -7,73 E1  0, E2 -9,69 -9,37 -8,78 -8,04 -6.44 E1  1, E2 -9,01 -8,51 -8,54 -7,95 -7,27 E1  1, 5E2 -7,91 -8,25 -7,55 -7,32 -6.60 E1  2, E2 -6,56 -6,88 -6,47 -6,52 -5.53 So sánh chuyển vị lớn theo phương x y  bs / với trường hợp E1  0, E2 E1  0, E2 X X X X X E1  0, E2 -3.7% -0.3% 3.5% 4.5% -16,8% E1  1, E2 -10.5% -9.5% 0.7% 3.3% -6,0% E1  1, 5E2 -21.4% -12.2% -10.9% -4.9% -14,7% E1  2, E2 -34.8% -26.8% -23.7% -15.3% -28,6% Từ Hình 23 ta thấy dạng đồ thị chuyển vị theo phương x đường thẳng y  bs / gần không đổi tỉ lệ Ls / bs thay đổi Theo kết Bảng 12, trường hợp bs  60 m bs  48 m module đàn hồi tăng dần chuyển vị lớn theo phương x giảm dần, phù hợp với toán khảo sát bên Nhưng trường hợp lại module đàn hồi tăng dần giá trị chuyển vị lớn theo phương x lại khơng ổn định Vì lúc 51 Chương 3: Ví dụ số kích thước nhỏ, chiều rộng bs  12 m , bs  24 m , bs  36 m tải trọng di động tác dụng lên lan truyền nhanh đến biên sau tác dụng ngược lại nên làm cho chuyển động ổn định Bảng 13 kết khảo sát chuyển vị lớn biên theo phương x y  Chuyển vị lúc giảm dần giá trị module đàn hồi tăng dần Bảng 13 Chuyển vị lớn theo phương x  y   thay đổi chiều rộng Module đàn hồi Chiều rộng kết cấu bs  m  bs  48 m bs  36 m bs  24 m bs  12 m E1  0, E2 2.42 3.58 6.41 13.90 E1  0, E2 1.65 2.67 6.02 11.31 E1  1, E2 1.36 1.92 4.87 10.79 E1  1, 5E2 1.27 1.57 3.51 8.65 E1  2, E2 0.99 1.52 3.18 5.19 So sánh chuyển vị lớn theo phương x y  với trường hợp E1  0, E2 E1  0, E2 X X X X E1  0, E2 -31.6% -25.4% -6.2% -18.6% E1  1, E2 -43.7% -46.3% -24.1% -22.4% E1  1, 5E2 -47.3% -56.1% -45.3% -37.8% E1  2, E2 -59.2% -57.5% -50.5% -62.7% 52 Chương 4: Kết luận kiến nghị CHƯƠNG KẾT LUẬN VÀ KIẾN NGHỊ Luận văn sử dụng phương pháp kết hợp phần tử biên (BEM) phương pháp phần tử hữu hạn (FEM) phân tích trực hướng chịu tải trọng di động để giải toán VLFS Phân tích ứng xử động trực hướng với giá trị độ cứng uốn theo phương thay đổi ứng với thay đổi vận tốc lực di chuyển, tải trọng lực di chuyển, chiều dày chiều rộng Qua kết phân tích số đạt trình bày Chương 4, số kết luận quan trọng kiến nghị hướng phát triển đề tài tương lai rút 4.1 Kết luận Mơ hình phản ánh làm việc hợp lý kết cấu trực hướng nước Mơ hình đảm bảo độ tin cậy, độ xác xu hướng hợp lý việc xác định ứng xử động lực học trực hướng tác dụng tải trọng xe di chuyển Đối với tải trọng không đổi, vận tốc thay đổi kết cấu tác dụng tải trọng di động có chuyển vị lớn giảm dần giá trị module đàn hồi tăng dần Kết khảo sát với trường hợp vận tốc cho thấy tồn vận tốc tới hạn Ảnh hưởng thông số đến ứng xử động trực hướng chịu tải trọng di động như: chiều dày chiều rộng có tính chất chung độ lớn số tăng chuyển vị giảm Ngược lại độ lớn tải trọng tăng chuyển vị tăng Đồng thời thay đổi thông số chu kì chuyển vị khơng thay đổi 53 Chương 4: Kết luận kiến nghị Kiến nghị 4.2 Mặc dù luận văn đạt số kết định trình bày số vấn đề chưa khám phá cần nghiên cứu thêm tương lai Những vấn đề bao gồm: - Mở rộng mơ hình trực hướng sang nhiều lớp phân lớp chức - Xét thêm trường hợp tải trọng hướng nghiên cứu cần khai thác 54 Tài liệu tham khảo TÀI LIỆU THAM KHẢO [1] Suzuki, H., & Yoshida, K (1996, November) Design flow and strategy for safety of very large floating structure In Proceedings of Int Workshop on Very Large Floating Structures, VLFS (Vol 96, pp 21-27) [2] Khabakhpasheva, T I., & Korobkin, A.A (2002) Hydroelastic behaviour of compound floating plate in waves Journal of engineering mathematics, 44(1), 21-40 [3] Kashiwagi, M (2004) Transient responses of a VLFS during landing and take-off of an airplane Journal of Marine Science and Technology, 9(1), 1423 [4] Watanabe E, Utsunomiya T, Wang CM and Le Thi Thu Hang Benchmark hydroelastic responses of circular VLFS underwave action Engrg Struct.; 2000, 28(6), pp 423–30 [5] Fujikubo, M., & Yao, T (2001) Structural modeling for global response analysis of VLFS Marine structures, 14(3), 295-310 [6] Jin, J (2008) A mixed mode function: boundary element method for very large floating structure: water interaction systems excited by airplane landing impacts (Doctoral dissertation, University of Southampton) [7] C.M Wang, E Watanabe, T Utsunomiya, (Eds.) (2007) Very large floating structures CRC Press [8] R Szilard, (2004) Theories and applications of plate analysis: classical numerical and engineering methods John Wiley & Sons [9] H Suzuki, H.R Riggs, M Fujikubo, T.A Shugar, H Seto, Y Yasuzawa, &H Shin, (2007, January) Very large floating structures In ASME 2007 26th International Conference on Offshore Mechanics and Arctic Engineering (pp 597-608) American Society of Mechanical Engineers 55 Tài liệu tham khảo [10] Kashiwagi, M (2000a) “Research on hydroelastic responses of VLFS: Recent progress and future work,” International Journal of Offshore and Polar Engineering, 10, pp 81–90 [11] Watanabe, E., Utsunomiya, T., & Wang, C M (2004) Hydroelastic analysis of pontoon-type VLFS: a literature survey Engineering structures, 26(2), 245-256 [12] Newman, J N (2005) Efficient hydrodynamic analysis of very large floating structures Marine structures, 18(2), 169-180 [13] Ohmatsu, S (2005) Overview: Research on wave loading and responses of VLFS Marine Structures, 18(2), 149-168 [14] Suzuki, H., Bhattacharya, B., Fujikubo, M., Hudson, D A., Riggs, H R., Seto, H., & Zong, Z (2006) ISSC committee VI 2: Very large floating structures In Sixteenth International Ship and Offshore Structures Congress, Southampton, UK (pp 394-442) [15] Yasuzawa, Y., Kagawa, K., Kawano, D and Kitabayashi, K 1997 Dynamic response of alarge flexible floating structure in regular waves Proceedings of the 16th Conference on Off shore Mechanics and Arctic Engineering Yokohama, Japan, 187-194 [16] Yago, K and Endo, H 1996 On the hydroelastic response of box-shaped floating structure with shallow draft Journal of the Society of Naval Architects of Japan 180 341-352 [17] Yasuzawa, Y., Kagawa, K., Kawano, D and Kitabayashi, K 1996 Wave response analysis of a flexible large floating structure Proceedings of the 2nd International Workshop onVery Large Floating Structures Hayama, Japan, 221-228 [18] Hamamoto, T., Suzuki, A and Fujita, K 1997 Hybrid dynamic analysis of large tension legfloating structures using plate elements Proceedings of the 7th International Offshore and Polar Engineering Conference Honolulu, HI, 285-292 56 Tài liệu tham khảo [19] Ohkusu, M., & Namba, Y (1996) Analysis of hydroelastic behavior of a large floating platform of thin plate configuration in waves VLFS, 96, 143-148 [20] Namba, Y., & Ohkusu, M (1999) Hydroelastic behavior of floating artificial islands in waves International Journal of Offshore and Polar Engineering, 9(01) [21] Mamidipudi, P (1994) The Motion Performance of a Mat-like Floating Airports In Proc of Int Conf on Hydroelasticity in Marine Technology, 1994 (pp 363-375) [22] Wu, C., Watanabe, E., & Utsunomiya, T (1995) An eigenfunction expansionmatching method for analyzing the wave-induced responses of an elastic floating plate Applied Ocean Research, 17(5), 301-310 [23] Kim, J W., & Webster, W C (1998) The drag on an airplane taking off from a floating runway Journal of Marine science and Technology, 3(2), 76-81 [24] Watanabe, E., Utsunomiya, T., & Tanigaki, S (1998) A transient response analysis of a very large floating structure by finite element method Doboku Gakkai Ronbunshu, 1998(598), 1-9 [25] Watanabe, E (1996) Transient response analysis of VLFS at airplane landing In Proceedings of the 2nd International Workshop on Very Large Floating Structures, Hayama (pp 243-247) [26] Endo, H (2000) The behavior of a VLFS and an airplane during takeoff/landing run in wave condition Marine structures, 13(4), 477-491 [27] Kashiwagi, M (2000b) A time-domain mode-expansion method for calculating transient elastic responses of a pontoon-type VLFS Journal of Marine Science and Technology, 5(2), 89-100 [28] Ismail, R E S (2016) Time-Domain Three Dimensional BE-FE Method for Transient Response of Floating Structures Under Unsteady Loads Latin American Journal of Solids and Structures, 13(7), 1340-1359 57 ... TÊN ĐỀ TÀI: PHƯƠNG PHÁP KẾT HỢP GIỮA PHẦN TỬ BIÊN VÀ PHẦN TỬ HỮU HẠN TRONG PHÂN TÍCH TẤM NỔI TRỰC HƯỚNG CHỊU TẢI TRỌNG DI ĐỘNG II NHIỆM VỤ VÀ NỘI DUNG Trình bày sở lý thuyết trực hướng, chất... Trong luận văn sử dụng phương pháp kết hợp phần tử biên (BEM) phương pháp phần tử hữu hạn (FEM) phân tích động lực học trực hướng chịu tải trọng di động để giải toán VLFS Hướng nghiên cứu 1.4.2... nhiều phương pháp số phát triển để phân tích VLFS mơi trường sóng biển, luận văn sử dụng phương pháp kết hợp phần tử biên (BEM) phương pháp phần tử hữu hạn (FEM) để phân tích kết cấu trực hướng chịu