Tối ưu hóa dựa trên độ tin cậy kết cấu bồn gia cường chứa dầu bằng phương pháp lai pso sqp

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC MỞ TP HỒ CHÍ MINH LÊ ĐÌNH VĂN TỐI ƯU HĨA DỰA TRÊN ĐỘ TIN CẬY KẾT CẤU BỒN GIA CƯỜNG CHỨA DẦU BẰNG PHƯƠNG PHÁP LAI PSO - SQP Tai Lieu Chat Luong LUẬN VĂN THẠC SỸ XÂY DỰNG CÔNG TRÌNH DÂN DỤNG VÀ CƠNG NGHIỆP TP HỒ CHÍ MINH, NĂM 2015 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC MỞ TP HỒ CHÍ MINH LÊ ĐÌNH VĂN TỐI ƯU HÓA DỰA TRÊN ĐỘ TIN CẬY KẾT CẤU BỒN GIA CƯỜNG CHỨA DẦU BẰNG PHƯƠNG PHÁP LAI PSO - SQP Chun ngành: Xây dựng cơng trình dân dụng công nghiệp Mã số chuyên ngành: 60 58 02 08 LUẬN VĂN THẠC SỸ Người hướng dẫn khoa học PGS.TS NGUYỄN THỜI TRUNG Tp Hồ Chí Minh, năm 2015 -i- LỜI CAM ĐOAN Luận văn “Tối ưu hóa dựa độ tin cậy kết cấu bồn gia cường chứa dầu phương pháp lai PSO – SQP” nghiên cứu tơi thực hướng dẫn PGS.TS Nguyễn Thời Trung Tôi xin cam đoan 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 tất công việc mà thực luận văn Tp.HCM, ngày … tháng … năm 2015 Học viên thực LÊ ĐÌNH VĂN - iii - TÓM TẮT LUẬN VĂN TÊN ĐỀ TÀI “TỐI ƯU HÓA DỰA TRÊN ĐỘ TIN CẬY KẾT CẤU BỒN GIA CƯỜNG CHỨA DẦU BẰNG PHƯƠNG PHÁP LAI PSO – SQP” Luận văn nhằm thực đề tài “Tối ưu hóa dựa độ tin cậy kết cấu bồn gia cường chứa dầu phương pháp lai PSO – SQP” Mơ hình tính tốn kết cấu bồn gia cường chứa dầu vỏ có gia cường dầm với ngàm biên đáy bồn Phương pháp phần tử hữu hạn (Finite Element Method – FEM) với phần tử CS-DSG3 (Cellbased Smoothed Discrete Shear Gap) sử dụng để phân tích ứng xử cho tốn vỏ gia cường dầm Bài toán thiết kế tối ưu dựa độ tin cậy (Reliability-Based Design Optimization – RBDO) thiết lập với hàm mục tiêu trọng lượng thành bồn, biến thiết kế chiều dày thành bồn kích thước mặt cắt ngang dầm Trong đó, phương pháp lai PSO – SQP (là kết hợp phương pháp: tối ưu hóa bầy đàn (Particle Swarm Optimization - PSO) giải thuật bình phương (Sequential Quadratic Programming – SQP)) Phương pháp độ tin cậy bậc ngược (Inverse First Order Reliability Method – Inv-FORM) sử dụng để giải toán tối ưu hóa tốn đánh giá độ tin cậy vịng lặp giải thuật tối ưu hóa dựa độ tin cậy (Sequential Optimization and Reliability Assessement – SORA) Việc giải tốn tối ưu hóa dựa độ tin cậy kết cấu bồn gia cường thực qua bước sau: Bước 1: Phân tích ứng xử kết cấu bồn gia cường dầm phần tử CS-DGG3 Bước 2: Thiết lập giải toán thiết kế tối ưu kết cấu bồn gia cường chứa dầu phương pháp lai PSO - SQP Bước 3: Thiết lập giải toán thiết kế tối ưu dựa độ tin cậy (RBDO) kết cấu bồn gia cường chứa dầu hai phương pháp PSO – SQP phương pháp độ tin cậy bậc ngược (Inv-FORM) - iv - Các kết phân tích ứng xử kiểm chứng kết thu từ phần mềm ANSYS lời giải giải tích tốn vỏ trụ trịn, trục thẳng đứng chịu uốn lý thuyết vỏ Đối với toán tối ưu hóa bước 2, để đánh giá hiệu phương pháp lai PSO – SQP nghiệm tối ưu đạt toán so sánh với nghiệm tối ưu phương pháp tiến hóa cải tiến (Improved Constrained Diffirential Evolution – ICDE) phương pháp PSO Nghiệm tối ưu thu sau giải toán thiết kế tối ưu dựa độ tin cậy phương pháp lai PSO – SQP bước kiểm tra, so sánh với thiết kế thực tế Từ khóa: Phần tử hữu hạn (PTHH); phương pháp lai tối ưu hóa bầy đàn - giải thuật bình phương (PSO – SQP); phương pháp thiết kế tối ưu dựa độ tin cậy (RBDO); thuật giải tối ưu dựa độ tin cậy (SORA) -1- MỤC LỤC LỜI CAM ĐOAN i LỜI CẢM ƠN .ii TÓM TẮT LUẬN VĂN iii MỤC LỤC DANH MỤC BẢNG BIỂU DANH MỤC HÌNH ẢNH DANH MỤC CÁC TỪ VIẾT TẮT TỔNG QUAN 1.1 Giới thiệu chung 1.2 Tình hình nghiên cứu giới 11 1.3 Tình hình nghiên cứu nước .12 1.4 Mục tiêu đề tài 13 1.5 Phạm vi nghiên cứu 13 CƠ SỞ LÝ THUYẾT 14 2.1 Lý thuyết vỏ Reissner – Mindlin gia cường dầm Timoshenko 14 2.1.1 Tấm Reissner – Mindlin 14 2.1.2 Vỏ Reissner - Mindlin 20 2.1.3 Dầm Timoshenko .22 2.1.4 Vỏ Reissner-Mindlin gia cường dầm Timoshenko 28 2.2 Phương pháp PTHH cho toán vỏ gia cường dầm .29 2.2.1 Phương pháp PTHH cho toán vỏ 30 2.2.2 Phương pháp PTHH cho toán dầm 37 2.2.3 Phương pháp PTHH cho toán vỏ gia cường dầm 40 2.3 Lý thuyết tối ưu hóa sử dụng phương pháp lai PSO – SQP 41 2.3.1 Giới thiệu 41 2.3.2 Phương pháp PSO 42 2.3.3 Phương pháp SQP 45 2.3.4 Phương pháp lai PSO – SQP .46 -2- 2.4 Lý thuyết thiết kế tối ưu dựa độ tin cậy (RBDO) .47 2.4.1 Phân tích độ nhạy .50 2.4.2 Phương pháp độ tin cậy bậc ngược Inv-FORM 50 2.4.3 Giải thuật RBDO dựa phương pháp lai PSO – SQP phương pháp Inv-FORM 53 KẾT QUẢ SỐ 57 3.1 Bài toán bồn gia cường chứa dầu có dung tích 5.000 m3 59 3.2 Bài toán bồn gia cường chứa dầu có dung tích 20.000 m3 66 3.3 Nhận xét kết 73 3.4 Mối quan hệ hệ số an toàn độ tin cậy kết cấu 74 3.5 So sánh kết luận văn với kết thiết kế thực tế 75 3.6 Khảo sát tốn tối ưu hóa dựa độ tin cậy ứng với tỉ lệ thay đổi khác biến ngẫu nhiên 76 KẾT LUẬN VÀ KHUYẾN NGHỊ 79 4.1 Kết luận 79 4.2 Khuyến nghị .80 4.2.1 Những hạn chế nghiên cứu 80 4.2.2 Hướng phát triển đề tài 81 TÀI LIỆU THAM KHẢO 82 PHỤ LỤC 87 -3- DANH MỤC BẢNG BIỂU Bảng 3.1 So sánh kết phân tích ứng xử vỏ không gia cường 60 Bảng 3.2 So sánh kết phân tích ứng xử vỏ gia cường dầm 61 Bảng 3.3 Kết phân tích độ nhạy hàm chuyển vị hàm ứng suất 62 Bảng 3.4 Nghiệm tốn tối ưu (khơng xét độ tin cậy) ba phương pháp PSO – SQP, ICDE PSO .62 Bảng 3.5 Kết tối ưu hóa tốn RBDO ứng với trường hợp xác suất an toàn khác 63 Bảng 3.6 Chênh lệch giá trị hàm mục tiêu hai trường hợp tính tốn 64 Bảng 3.7 So sánh nghiệm tối ưu hai trường hợp 65 Bảng 3.8 So sánh kết phân tích ứng xử vỏ không gia cường 67 Bảng 3.9 So sánh kết phân tích ứng xử vỏ có gia cường 68 Bảng 3.10 Kết phân tích độ nhạy hàm chuyển vị hàm ứng suất 69 Bảng 3.11 Nghiệm tốn tối ưu (khơng xét độ tin cậy) ba phương pháp PSO – SQP, ICDE PSO .70 Bảng 3.12 Kết tối ưu hóa toán RBDO ứng với trường hợp xác suất an toàn khác 71 Bảng 3.13 Chênh lệch giá trị hàm mục tiêu hai trường hợp tính toán 72 Bảng 3.14 So sánh nghiệm tối ưu hai trường hợp 73 Bảng 3.15 Mối liên quan hệ số an toàn độ tin cậy kết cấu 74 Bảng 3.16 Kết thu từ luận văn thiết kế thực tế 75 Bảng 3.17 Khảo sát nghiệm toán RBDO (bồn 5.000 m3) thay đổi tỉ lệ thay đổi V biến ngẫu nhiên ν P 77 Bảng 3.18 Khảo sát nghiệm toán RBDO (bồn 20.000 m3) thay đổi tỉ lệ thay đổi V biến ngẫu nhiên ν P 77 -4- DANH MỤC HÌNH ẢNH Hình 1.1 Bồn chứa dầu nhà máy lọc dầu Dung Quất Hình 1.2 Dự án kho xăng dầu hàng không miền Nam .8 Hình 1.3 Hình dáng chung bồn chứa dầu với mái kết cấu trụ cầu Hình 1.4 Mặt cắt đứng bồn chứa dầu Hình 2.1 Phần tử vỏ thoải 14 Hình 2.2 Quy ước dấu 15 Hình 2.3 a) Quy ước dấu mặt phẳng Oxz; b) Quy ước dấu mặt phẳng Oyz 16 Hình 2.4 Biến dạng 17 Hình 2.5 Chuyển hệ tọa độ cho phần tử tam giác 20 Hình 2.6 Chuyển hệ tọa độ cho phần tử dầm 22 Hình 2.7 Quy ước dấu cho dầm .23 Hình 2.8 Chuyển vị biến dạng dầm 24 Hình 2.9 Đổi biến dầm 26 Hình 2.10 Phần tử vỏ gia cường dầm 29 Hình 2.11 Phần tử tam giác miền tham chiếu 32 Hình 2.12 Tọa độ nút hệ tọa độ tổng thể 33 Hình 2.13 Phần tử CS-DSG3 chia nhỏ thành tam giác 35 Hình 2.14 Phần tử dầm hai nút 37 Hình 2.15 Phần tử dầm miền tham chiếu 39 Hình 2.16 Sơ đồ thuật giải PSO 43 Hình 2.17 Khởi tạo quần thể ban đầu 44 Hình 2.18 Cập nhật vận tốc vị trí cá thể 45 Hình 2.19 Sơ đồ khối phương pháp lai PSO – SQP 47 Hình 2.20 Kết tối ưu xét đến độ tin cậy không xét đến độ tin cậy 48 Hình 2.21 Hai vòng lặp riêng biệt giải thuật SORA 49 Hình 2.22 Điểm thiết kế MPP không gian vật lý không gian chuẩn hóa 51 Hình 2.23 Sơ đồ thuật tốn Inv-FORM 53 -5- Hình 2.24 Mơ hình chuyển đổi hàm trạng thái giới hạn 54 Hình 2.25 Biến đổi tương đương hai hàm ràng buộc xác suất 54 Hình 2.26 Hiệu chỉnh giá trị MPP 55 Hình 2.27 Sơ đồ giải thuật RBDO dựa phương pháp lai PSO – SQP 56 Hình 3.1 Mơ hình phân tích ứng xử ¼ thành bồn phần tử CS-DSG3 với lưới 10x10 .57 Hình 3.2 Mơ hình thành bồn gia cường dầm 5.000 m3 59 Hình 3.3 Chuyển vị vỏ không gia cường a) Lý thuyết vỏ b) ANSYS 60 Hình 3.4 Kết biến dạng vỏ khơng gia cường a) trước có tải b) sau có tải .60 Hình 3.5 Chuyển vị vỏ gia cường chạy từ phần mềm ANSYS 61 Hình 3.6 Kết biến dạng vỏ gia cường a) trước có tải b) sau có tải 61 Hình 3.7 Thay đổi hàm mục tiêu ứng với thay đổi xác suất an toàn 64 Hình 3.8 So sánh giá trị hàm mục tiêu hai trường hợp tính tốn 65 Hình 3.9 Mơ hình thành bồn gia cường dầm 20.000 m3 66 Hình 3.10 Chuyển vị vỏ không gia cường a) lý thuyết vỏ b) ANSYS 67 Hình 3.11 Kết biến dạng vỏ không gia cường a) trước có tải b) sau có tải 68 Hình 3.12 Chuyển vị vỏ có gia cường chạy từ phần mềm ANSYS 68 Hình 3.13 Kết biến dạng vỏ gia cường a) trước có tải b) sau có tải 69 Hình 3.14 Thay đổi hàm mục tiêu ứng với thay đổi xác suất an toàn 71 Hình 3.15 So sánh giá trị hàm mục tiêu hai trường hợp tính tốn 72 Hình 3.16 Khảo sát mối liên quan hệ số an toàn độ tin cậy a) Trường hợp bồn 5.000 m3 b) Trường hợp bồn 20.000 m3 75 - 72 - tin cậy Tỉ lệ tăng dao động từ 12,40% đến 30,78% tùy theo giá trị xác suất an toàn Ngoài ra, ta thấy rõ xác suất an toàn cao tỉ lệ tăng cao Điều cho thấy việc kết hợp tốn phân tích độ tin cậy với tốn tối ưu hóa quan trọng Nhất xác suất an toàn ấn định phạm vi lớn 99,98% Bảng 3.13 Chênh lệch giá trị hàm mục tiêu hai trường hợp tính tốn Hình 3.15 So sánh giá trị hàm mục tiêu hai trường hợp tính tốn Tương ứng với giá trị hàm mục tiêu, giá trị nghiệm tối ưu toán thiết kế tối ưu dựa độ tin cậy tăng lên so với giá trị nghiệm tối ưu tốn tối ưu hóa không xét độ tin cậy Bảng 3.13 so sánh giá trị nghiệm tối ưu trường hợp thiết kế tối ưu không xét độ tin cậy với giá trị nghiệm tối ưu tối ưu có ràng buộc độ tin cậy (ứng với ràng buộc xác suất an toàn Pr(  [])  99,38% (tương ứng số độ tin cậy β = 2,5) kết cấu) Kết cho thấy rằng, giá trị biến chiều dày toán - 73 - tối ưu có xét độ tin cậy (t = 15,674 mm) tăng đến 24,26% so với giá trị biến chiều dày toán tối ưu không xét độ tin cậy (t = 12,614 mm) Bảng 3.14 So sánh nghiệm tối ưu hai trường hợp 3.3 Nhận xét kết Dựa kết tốn tối ưu 3.1 3.2, ta có nhận xét sau:  Kết phân tích ứng xử kết cấu bồn gia cường phần tử CSDSG3 đáng tin cậy sử dụng để tính tốn cho tốn tối ưu toán thiết kế tối ưu dựa độ tin cậy;  Trong toán thiết kế tối ưu hóa khơng xét đến độ tin cậy kết cấu bồn gia cường, phương pháp lai PSO – SQP sử dụng để tìm kiếm nghiệm tối ưu Khi so sánh kết với hai phương pháp khác ICDE PSO, phương pháp lai PSO – SQP cho kết xác tiết kiệm chi phí tính tốn nhiều Phương pháp ICDE cho kết xác tương tự tốn nhiều thời gian tính tốn hơn, phương pháp PSO có thời gian tính tốn nhỏ phương pháp ICDE lớn phương pháp PSO – SQP kết tối ưu lại khơng xác hai phương pháp lại Như vậy, việc sử dụng phương pháp lai PSO – SQP hợp lý hiệu toán tối ưu hóa bồn gia cường;  Đối với tốn thiết kế tối ưu xét đến độ tin cậy, ta sử dụng giải thuật SORA cách kết hợp phương pháp lai PSO – SQP cho vòng lặp tối ưu hóa phương pháp Inv-FORM cho vịng lặp độ tin cậy Từ kết thu được, ta thấy giá trị hàm mục tiêu tăng so với không xét đến độ tin cậy Giá trị hàm mục tiêu tỉ lệ thuận với xác suất an toàn ấn định trước Đặc biệt, giá trị hàm mục tiêu tăng mạnh ứng với xác suất an toàn kết cấu ấn định phạm vi lớn 99,98% - 74 - 3.4 Mối quan hệ hệ số an toàn độ tin cậy kết cấu Trong tiêu chuẩn thiết kế, để nâng cao xác suất an tồn cơng trình, người thiết kế sử dụng tĩnh tải tính tốn tải tiêu chuẩn nhân thêm hệ số an tồn ( = 1,1 ÷ 1,2) Tuy nhiên, điều có ý nghĩa mặt định tính, người thiết kế không xác định mặt định lượng xác suất an tồn cơng trình cụ thể Với toán thiết kế tối ưu dựa độ tin cậy, người kỹ sư đảm bảo thiết kế vừa tối ưu mặt kinh tế vừa đảm bảo xác suất an toàn cho phép cơng trình Như vậy, ứng với giá trị thiết kế từ hai toán trên, ta thấy có liên hệ tương ứng hệ số an toàn với xác suất an toàn (hay số độ tin cậy β) Trong ngữ cảnh luận văn này, ta định nghĩa hệ số an toàn (Safety Factor – SF) tỉ số ứng suất cho phép ứng suất tính tốn, cách sử dụng phương pháp phân tích độ tin cậy bậc ngược Inv-FORM, ta xác định mối quan hệ hệ số an toàn độ tin cậy kết cấu toán 3.1 3.2, minh họa Bảng 3.15 Hình 3.16 Bảng 3.15 Mối liên quan hệ số an toàn độ tin cậy kết cấu - 75 - Kết Bảng 3.15 thiết thực cho phép người thiết kế chọn hệ số an toàn phù hợp để thiết kế (tức khơng cần giải tốn tối ưu) mà đảm bảo xác suất an toàn ấn định trước chủ đầu tư a) b) Hình 3.16 Khảo sát mối liên quan hệ số an toàn độ tin cậy a) Trường hợp bồn 5.000 m3 b) Trường hợp bồn 20.000 m3 3.5 So sánh kết luận văn với kết thiết kế thực tế Các cơng trình bồn chứa dầu thực tế như: Tổng kho xăng dầu Nhà Bè, Kho xăng dầu Chân Mây, Kho xăng dầu hàng không Miền Nam, v.v có phần kết cấu thép bồn thiết kế theo tiêu chuẩn API 650 11th, ADD.2008 Viện Dầu khí Quốc gia Mỹ Từ kết thiết kế thực tế bồn chứa 5.000 m3 bồn chứa 20.000 m3 sử dụng kết Bảng 3.15 sử dụng phương pháp phân tích độ tin cậy bậc FORM, ta thu kết từ luận văn (gồm xác suất an toàn số độ tin cậy β) ứng với hệ số an toàn thiết kế thực tế Kết thu từ nghiên cứu thiết kế thực tế trình bày Bảng 3.16 Bảng 3.16 Kết thu từ luận văn thiết kế thực tế Từ mục 3.4 Bảng 3.16, ta có nhận xét kết hàm mục tiêu thu nghiên cứu thiết kế thực tế sau: - 76 -  Đối với bồn chứa dầu dung tích 5.000 m3, với hệ số an tồn SF = 1,225 thiết kế thực tế, ta tính xác suất an toàn kết cấu bồn 99,754% giá trị hàm mục tiêu hay trọng lượng thành bồn 62.795 kg Giá trị hàm mục tiêu có mức chênh lệch nhỏ so với giá trị hàm mục tiêu thiết kế thực tế 0,29%;  Đối với bồn chứa dầu dung tích 20.000 m3, với hệ số an toàn SF = 1,214 thiết kế thực tế, ta tính xác suất an tồn kết cấu bồn 99,629% giá trị hàm mục tiêu hay trọng lượng thành bồn 334.729 kg Giá trị hàm mục tiêu có mức chênh lệch nhỏ so với thiết kế thực tế 2,28%;  Như vậy, giá trị hàm mục tiêu thu toán tối ưu dựa độ tin cậy từ nghiên cứu chênh lệch nhỏ so với thiết kế thực tế Điểm vượt trội từ kết nghiên cứu việc xác định cụ thể xác suất an toàn kết cấu ứng với hệ số an toàn sử dụng thiết kế thực tế Trong khi, với cách thiết kế sử dụng hệ số an tồn truyền thống ta khơng thể xác định xác suất an toàn kết cấu;  Cũng cần lưu ý rằng, kết đạt tốn tối ưu hóa dựa độ tin cậy mục tương ứng với tỉ lệ thay đổi V hai biến ngẫu nhiên ν P 10% 8%, tức ứng với độ lệch phân bố ngẫu nhiên lớn Nếu độ lệch phân bố ngẫu nhiên nhỏ (tức tỉ lệ thay đổi V nhỏ hơn) tăng tỉ lệ thay đổi biến ngẫu nhiên P kết đạt tốn tối ưu hóa dựa độ tin cậy khác biệt rõ so với kết thiết kế thực tế Để làm rõ ý phân tích này, học viên tiến hành khảo sát nghiệm tốn tối ưu hóa dựa độ tin cậy ứng với tỉ lệ thay đổi khác biến ngẫu nhiên 3.6 Khảo sát toán tối ưu hóa dựa độ tin cậy ứng với tỉ lệ thay đổi khác biến ngẫu nhiên Ta khảo sát nghiệm toán tối ưu dựa độ tin cậy ứng với trường hợp tỉ lệ thay đổi hai biến ngẫu nhiên ν P sau: - Trường hợp 1: Vν = 5%; VP = 5% - 77 - - Trường hợp 2: Vν = 5%; VP = 8% - Trường hợp 3: Vν = 10%; VP = 8% - Trường hợp 4: Vν = 10%; VP = 12% Các trường hợp tính tốn ứng với hệ số an toàn thiết kế thực tế thể Bảng 3.16 Nghiệm toán tối ưu dựa độ tin cậy ứng với trường hợp trình bày Bảng 3.17 (bồn chứa dầu 5.000 m3) Bảng 3.18 (bồn chứa dầu 20.000 m3) Bảng 3.17 Khảo sát nghiệm toán RBDO (bồn 5.000 m3) thay đổi tỉ lệ thay đổi V biến ngẫu nhiên ν P Bảng 3.18 Khảo sát nghiệm toán RBDO (bồn 20.000 m3) thay đổi tỉ lệ thay đổi V biến ngẫu nhiên ν P Hai Bảng 3.17 3.18 cho ta thấy nghiệm tốn RBDO khơng thay đổi nhiều tỉ lệ thay đổi biến ngẫu nhiên P VP = 8% biến ngẫu nhiên ν Vν = 5% ÷ 10% giá trị hàm mục tiêu từ nghiên cứu không khác biệt nhiều so với giá trị hàm mục tiêu thiết kế thực tế Khi tỉ lệ thay đổi V biến ngẫu nhiên nhỏ (trường hợp 1), giá trị hàm mục tiêu tính tốn nhỏ so với giá trị hàm mục tiêu thiết kế thực tế 7,08% (bồn 5.000m3) 8,94% (bồn 20.000 m3) Khi tỉ lệ thay đổi biến ngẫu nhiên lớn (trường hợp 4), giá trị hàm mục tiêu tính tốn lớn so với giá trị hàm mục tiêu thiết kế thực tế 9,01% (bồn 5.000m3) 8,07% (bồn 20.000 m3) Như vậy, hai trường hợp 4, kết đạt - 78 - tốn RBDO có khác biệt đáng kể so với kết thiết kế thực tế Kết cho thấy rằng, việc thiết kế tối ưu dựa độ tin cậy với đầy đủ liệu phân bố ngẫu nhiên thực tế cho kết thiết kế sát với điều kiện sử dụng thực tế mà đáp ứng xác suất an toàn cho trước - 79 - CHƯƠNG KẾT LUẬN VÀ KHUYẾN NGHỊ 4.1 Kết luận Luận văn thực nghiên cứu toán thiết kế tối ưu dựa độ tin cậy RBDO cho kết cấu bồn gia cường Trên sở kết hợp phương pháp phân tích kết cấu CS-DSG3 cho tốn phân tích ứng xử vỏ gia cường, phương pháp lai PSO – SQP cho toán tối ưu hóa phương pháp Inv-FORM cho phân tích độ tin cậy giải thuật tối ưu hóa dựa độ tin cậy SORA, học viên thiết lập giải thành cơng tốn thiết kế tối ưu dựa độ tin cậy kết cấu Bài toán xây dựng với hàm mục tiêu trọng lượng thành bồn, biến ngẫu nhiên hệ số poison ν, tải thủy tĩnh P biến thiết kế chiều dày thành bồn t, chiều rộng b, chiều cao h mặt cắt ngang dầm Các ràng buộc toán gồm: ràng buộc ứng suất   [], ràng buộc xác suất an toàn Pr (  []) ≥ [P] Đối với kết cấu phức tạp, hàm mục tiêu ràng buộc có tính phi tuyến cao, việc tìm lời giải tối ưu tồn cục khó khăn Việc sử dụng phương pháp kết hợp PSO SQP giúp giải hiệu vấn đề PSO mạnh tìm điểm tối ưu tồn cục SQP có ưu điểm hội tụ nhanh nghiệm tối ưu địa phương Đánh giá độ tin cậy với phương pháp Inv-FORM đảm bảo nghiệm tối ưu thỏa mãn xác suất an toàn ấn định cho trước Trong tốn tối ưu hóa kết cấu bồn gia cường, kết hợp phương pháp tối ưu toàn cục PSO phương pháp tối ưu địa phương SQP cho thấy tính hiệu cao Các kết tối ưu đạt từ phương pháp lai PSO – SQP đáng tin cậy hiệu số phương pháp so sánh ICDE PSO Kết tốn tối ưu khơng xét độ tin cậy thu cho thấy nghiệm tối ưu nằm ranh giới vùng an tồn khơng an tồn (qua kiểm tra, ta tính tốn xác suất an tồn 50,001% hay xác suất phá hủy 49,999%) Khi biến ngẫu nhiên thay đổi theo quy luật phân bố định (ảnh hưởng từ trình sản xuất, từ mơi trường xung quanh, từ q trình vận hành, v.v.) nghiệm tối ưu rơi vào vùng khơng an tồn Bài tốn thiết kế tối ưu dựa độ tin cậy áp - 80 - dụng đảm bảo nghiệm tối ưu nằm vùng an toàn ứng với xác suất an toàn cho trước Cùng hệ số an toàn với tỉ lệ thay đổi biến ngẫu nhiên Vν = 10%, VP = 8%, giá trị hàm mục tiêu thu từ luận văn tương đương với thiết kế thực tế nghiên cứu xác suất an toàn kết cấu Vì vậy, áp dụng tốn thiết kế tối ưu dựa độ tin cậy cho kết cấu, người thiết kế xác định nghiệm thiết kế tối ưu ứng với hệ số an tồn mà chủ đầu tư u cầu đồng thời xác định xác suất an toàn kết cấu Ngoài ra, độ lệch phân bố ngẫu nhiên (tức tỉ lệ thay đổi V) biến ngẫu nhiên thay đổi kết tốn tối ưu hóa dựa độ tin cậy thay đổi Kết nhỏ hay lớn so với kết thiết kế thực tế phụ thuộc vào thay đổi nhỏ hay lớn tỉ lệ thay đổi V Kết cho thấy rằng, việc thiết kế tối ưu dựa độ tin cậy với đầy đủ liệu phân bố ngẫu nhiên thực tế cho kết thiết kế sát với điều kiện sử dụng thực tế mà đáp ứng xác suất an toàn cho trước Kết số cho thấy, việc thiết kế tối ưu dựa độ tin cậy RBDO sử dụng phương pháp lai PSO – SQP cho vịng lặp tối ưu hóa giải thuật tách rời SORA (trong hai vịng lặp tối ưu đánh giá độ tin cậy thực độc lập), lựa chọn tối ưu để giải toán thiết kế tối ưu dựa độ tin cậy cho kết cấu bồn gia cường chứa dầu 4.2 Khuyến nghị 4.2.1 Những hạn chế nghiên cứu Mặc dù học viên cố gắng tiếp cận nhiều phương pháp tính tốn khác q trình nghiên cứu thời gian thực hạn chế nên đề tài dừng lại toán thiết kế tối ưu dựa độ tin cậy phần kết cấu thép thành bồn chứa Để nghiên cứu đưa vào ứng dụng thực tế, người thiết kế phải giải tiếp toán thiết kế tối ưu dựa độ tin cậy cho toàn kết cấu bồn chứa dầu điều kiện tải tác dụng khác như: động đất, gió, v.v Sau đó, người thiết kế phải thuyết phục khách hàng sử dụng sản phẩm thiết kế với nhiều ưu điểm so với thiết kế thực tế - 81 - Ngồi ra, học viên cần phải có số liệu thống kê thực tế của biến ngẫu nhiên như: giá trị trung bình, độ lệch chuẩn qui luật phân bố Điều đòi hỏi nhiều nỗ lực, thời gian đo đạc thống kê liệu 4.2.2 Hướng phát triển đề tài  Mở rộng toán thiết kế tối ưu dựa độ tin cậy RBDO cho toàn kết cấu thép bồn chứa dầu (thêm phần đáy mái bồn) Trong luận văn này, khối lượng công việc nhiều thời gian thực có giới hạn nên học viên chọn biến thiết kế tối ưu chiều dày thành bồn (t), kích thước mặt cắt ngang dầm (b,h) số lượng dầm chọn cố định bố trí theo chiều cao bồn Mặt khác, kết cấu bồn chứa dầu nên áp lực giảm dần theo chiều cao, việc tối ưu số lượng dầm gia cường tối ưu cách bố trí dầm gia cường đề tài nghiên cứu cần thiết  Phạm vi nghiên cứu tải trọng tác dụng áp lực thủy tĩnh Do cần nghiên cứu thêm tốn với loại tải tác dụng khác như: động đất, gió, áp lực động lực học, v.v  Luận văn áp dụng giải thuật SORA, vịng lặp tối ưu phương pháp lai PSO – SQP Phương pháp mở rộng để giải toán thiết kế tối ưu dựa độ tin cậy cho loại kết cấu, vật liệu loại chất lỏng khác  Trong luận văn này, thiết kế tối ưu dựa độ tin cậy nghiên cứu tốn bồn gia cường Vì vậy, ta sử dụng thuật giải khác để tính tốn so sánh với kết để có lời giải tối ưu Kết hợp với việc xác định thông số biến ngẫu nhiên qua thực tiễn, thí nghiệm, ta đưa nghiên cứu vào ứng dụng thực tế sản phẩm thiết kế đáng tin cậy - 82 - TÀI LIỆU THAM KHẢO [1] S Timoshenko S Woinowsky Krieger, Theory of Plates and Shells McGraw – Hill Book Company, Second Edition, 1987 [2] Lương Văn Hải, Bài giảng môn học kết cấu vỏ Trường Đại học Mở TP HCM, 2013 [3] Ugural, A.C.Tresses, in plates and shell, 1999 [4] Nguyễn Thời Trung Nguyễn Xuân Hùng, Phương pháp phần tử hữu hạn sử dụng Matlab Nhà xuất xây dựng, 2015 [5] Chu Quốc Thắng, Giáo trình Phương pháp phần tử hữu hạn Nhà xuất khoa học kỹ thuật, 1997 [6] Jacob Fish Ted Belytschko, A first course in Finite elements Printed and bound in Great Britain by Antony Rowe Ltd, Chippenham, Wiltshire, 2007 [7] Trần Ích Thịnh Ngơ Như Khoa, Phương pháp phần tử hữu hạn Nhà Xuất Hà Nội, 2007 [8] T Nguyen-Thoi, P Phung-Van, H Nguyen-Xuan, C Thai-Hoang, A cellbased smoothed discrete shear gap method (CS-DSG3) using triangular elements for static and free vibration analyses of Reissner-Mindlin plates International Journal for Numerical Methods in Engineering, 91(7), 705-741, 2012 [9] Andrzej S Nowak Kevin R Collins, Reliability of Structures The McGrawHill Companies, Inc, 2000 [10] Nguyễn Thời Trung, Bài giảng môn học Độ tin cậy kết cấu Trường Đại học Mở TP.HCM, 2013 [11] Russel Ebenhart James Kennedy, Particle Swarm Optimization Proc of IEEE International Conf on Neural Networks IV, 1942–1948, 1995 - 83 - [12] Riccardo Poli, James Kennedy, Tim Blackwell, Particle swarm optimization – An overview Swarm Intell 1, 33–57, 2007 [13] Qinghai Bai, Analysis of Particle Swarm Optimization Algorithm Computer and Information Science, 3(1), 2010 [14] Dian Palupi Rini, Siti Mariyam Shamsuddin, Siti Sophiyati Yuhaniz, Particle Swarm Optimization: Technique, System and Challenges International Journal of Computer Applications (0975-8887), 14(1), 2011 [15] Philip E Gill, Walter Murray, Michael A Saunders, An SQP Algorithm for Large scale constrained Optimization Society for Industrial and Applied Mathematics, 47(1), 99–131, 2005 [16] Paul T Boggs Jon W Tolle, Sequential quadratic programming Contribution of the National Institute of Standards and Technology, 1996 [17] Craig Travers Lawrence, A Computationally Efficient Feasible Sequential Quadratic Programming Algorithm Doctoral Dissertation, Copyright by Craig Travers Lawrence, 1998 [18] Soren Ebbesen, Pascal Kiwitz, Lino Guzzella, A Generic Particle Swarm Optimization Matlab Function American Control Conference Fairmont Queen Elizabeth, Montréal, Canada, 2012 [19] Zach D Richards, Constrained Particle Swarm Optimization for Sequential quadratic programming International Journal of Modelling, Identification and Control, 8(4), 2009 [20] Vagelis Plevris Manolis Papadrakakis, A Hybrid Particle Swarm – Gradient Algorithm for Global Structural Optimization Computer-Aided Civil and Infrastructure Engineering, 26, 48-68, 2011 [21] Shashank Gupta Nitin Narang, Integrated PSO-SQP technique for Shortterm Hydrothermal Scheduling International Journal of Advanced Research in Computer Engineering & Technology (IJARCET), 4(4), 2015 - 84 - [22] S Rahman D Wei, Design sensitivity and Reliability-based structural Optimization by Univariate decomposition Struct Multidisc Optim 35: 245–261, 2008 [23] Jian Tu, Kyung K Choi, Young H Park, A new study on Reliability based Design Optimization ASME Journal of Mechanical Design, 121(4), 557-564, 1999 [24] K Kolanek S Jendo, Reliability - Based Optimization of truss structures AMAS Course on Reliability - Based Optimization, 237-246, 2002 [25] Harish Agarwal, Reliability Based Design Optimization: Formulations and Methodologies Graduate Program in Aerospace and Mechanical Engineering, 2004 [26] Nikos D.Lagaros Vissarion Papadopoulos, Optimum design of shell structure with random geometric, meterial and thicness imperfections International Journal of Solids and Structure 43, 6948-6964, 2006 [27] Leandro Luis Corso Herbert Martins Gomes, Reliability based design optimization of trusses by heuristic algorithms 20th International Congress of Mechanical Engineering, November 15-20, Gramado, RS, Brazil, 2009 [28] Rafael Holdorf Lopeza, Didier Lemossea, Jos´e Eduardo Souza de Cursia, Jhojan Rojasb, Abdelkhalak El-Hamia, An Approach for the Reliability Based Design Optimization of Laminated Composites Engineering Optimization 00(00), 1–17, 2010 [29] Lê Trung Kiên, Tính tối ưu dàn phẳng sử dụng giải thuật di truyền Luận văn thạc sỹ Trường Đại học Bách Khoa Tp.HCM, 2000 [30] Nguyễn Tấn Dũng, Tính tối ưu dàn phẳng Luận văn thạc sĩ Trường Đại học Kiến Trúc Hà Nội, 2002 - 85 - [31] Hồ Hữu Vịnh, Trần Văn Dần, Phùng Văn Phúc, Lâm Phát Thuận, Lương Văn Hải, Nguyễn Thời Trung, Tối ưu hóa dựa độ tin cậy composite nhiều lớp chịu uốn Hội nghị học toàn quốc lần thứ IX, 2012 [32] Trần Văn Dần, Tối ưu hóa dựa độ tin cậy composite laminate giải thuật di truyền phần tử CS-DSG3 Luận văn thạc sĩ Trường đại học Bách Khoa TP.HCM, 2013 [33] Huỳnh Thanh Phương, Tối ưu hóa dựa độ tin cậy kết cấu dàn Luận văn thạc sỹ Trường Đại học Mở TP.HCM, 2013 [34] Nguyễn Viết Cường, Tối ưu hóa dựa độ tin cậy composite nhiều lớp sử dụng giải thuật lặp SORA Luận văn thạc sỹ Trường Đại học Công nghệ TP.HCM, 2015 [35] Xiaoping Du Wei Chen, Sequential Optimization and Reliability Asessement Method for Efficient Probabilistic Design ASME Journal of Mechanical Design Revision Manuscript Number 20302, 2004 [36] Guanbo Jia, Yong Wang, Zixing Cai, Yaochu Jin, An improved (µ+) – constrained differential evolution for constrained optimization Information Sciences 222, 302-322, 2013 [37] Xiaoping Du, First Order and Second Reliability Methods – Chapter University of Missouri Rolla, 2005 [38] Yibing Xiang Yongming Liu, Inverse FORM method for probabilistic fatigue prognosis 51st AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference, AIAA 2010-2684, 2010 [39] Nguyễn Thời Trung, Phân tích độ tin cậy cơng trình xây dựng: tổng quan, thách thức triển vọng Hội thảo khoa học lần V khoa Xây dựng Điện trường Đại học Mở TP Hồ Chí Minh, 2014 - 86 - [40] Kai-Uwe Bletzinger, Manfred Bischoff, Ekkehard Ramm, a unified approach for shear-locking-free triangular and rectangular shell finite elements COMPUTATIONAL MECHANICS New Trends and Applications S Idelsohn, E Onate and E Dvorkin (Eds.) ©CIMNE, Barcelona, Spain 1998 [41] P Phung-Van, T Nguyen-Thoi, H Luong-Van, C Cao-Hoang, H NguyenXuan, A cell-based smoothed discrete shear gap method (CS-FEM-DSG3) using layerwise deformation theory for dynamic response of composite plates resting on viscoelastic foundation Computer Methods in Applied Mechanics and Engineering 272, 138-159, 2014