Luận văn thạc sĩ về tối ưu hóa độ tin cậy tấm composite laminate sử dụng thuật toán di truyền

Luận văn thạc sĩ trình bày tối ưu hóa độ tin cậy tấm composite laminate bằng giải thuật di truyền và phần tử csdsg3 trong kỹ thuật xây dựng.

• Tác giả: Trần Văn Dần
• Người hướng dẫn: TS. Lương Văn Hải
• Trường học: Đại học Bách Khoa - ĐHQG TP.HCM
• Chuyên ngành: Xây dựng công trình dân dụng và công nghiệp
• Đề tài: Tối ƣu Hóa Dựa Trên Độ Tin Cậy Tấm Composite Laminate Bằng Giải Thuật Di Truyền
• Loại tài liệu: luận văn thạc sĩ
• Năm xuất bản: 2013
• Địa điểm: TP. HCM

Mục lục chi tiết

• Lời cảm ơn
• Lời cam đoan
• Danh mục hình vẽ
• Danh mục bảng biểu
• Danh mục ký hiệu
• 1. CHƯƠNG 1: GIỚI THIỆU TỔNG QUAN
• 1.1. Đặt vấn đề và hướng giải quyết
• 1.2. Hướng giải quyết
• 1.2.1. Vật liệu composite
• 1.2.2. Phương pháp phần tử hữu hạn trơn
• 1.2.3. Các phương pháp đánh giá độ tin cậy
• 1.2.4. Các phương pháp tối ưu hóa
• 1.3. Tình hình nghiên cứu
• 1.4. Mục tiêu nghiên cứu
• 1.5. Ý nghĩa khoa học và thực tiễn
• 1.6. Cấu trúc luận văn
• 2. CHƯƠNG 2: CƠ SỞ LÝ THUYẾT
• 2.1. Lý thuyết tấm biến dạng cắt bậc nhất [45]
• 2.2. Tấm composite nhiều lớp theo lý thuyết biến dạng cắt bậc nhất [1-3]
• 2.2.1. Ứng xử cơ học của lớp composite đơn
• 2.2.2. Ứng xử cơ học của tấm composite nhiều lớp
• 2.3. Phần tử CS-DSG3 cho tấm composite nhiều lớp
• 2.3.1. Phần tử DSG3 [5]
• 2.3.2. Phần tử CS-DSG3 [6]
• 2.3.3. Thiết lập bài toán phân tích tĩnh và động tấm composite nhiều lớp theo lý thuyết biến dạng cắt bậc nhất sử dụng phần tử CS-DSG3
• 2.4. Phương pháp phân tích độ tin cậy bậc nhất (FORM) [7]
• 2.5. Giải thuật di truyền [46, 47]
• 2.5.1. Lựa chọn cá thể
• 2.5.2. Lai ghép cá thể
• 2.6. Các giải thuật tối ưu hóa dựa trên độ tin cậy
• 2.6.1. Tổng quan các giải thuật [48]
• 2.6.2. Một giải thuật đơn giản và hiệu quả để giải bài toán RBDO
• 3. CHƯƠNG 3: KẾT QUẢ SỐ
• 3.1. Bài toán 1: Tối ưu hóa không ràng buộc
• 3.2. Bài toán 2: Tối ưu hóa có ràng buộc, chưa xét độ tin cậy
• 3.3. Bài toán 3: Phân tích độ nhạy
• 3.4. Bài toán 4: Tối ưu hóa có ràng buộc, có xét độ tin cậy
• 3.4.1. Kiểm tra giải thuật MSA, so sánh với giải thuật RIA
• 3.4.2. Áp dụng giải thuật MSA cho bài toán RBDO
• 4. CHƯƠNG 4: KẾT LUẬN VÀ HƯỚNG PHÁT TRIỂN
• 4.1. Tóm tắt các kết quả đã đạt được
• 4.2. Hướng phát triển
• Tài liệu tham khảo

Tóm tắt

I. Tổng quan về tối ưu hóa độ tin cậy tấm composite laminate

Tấm composite laminate là một trong những vật liệu được sử dụng phổ biến trong ngành xây dựng và hàng không nhờ vào những đặc tính vượt trội như độ bền cao và trọng lượng nhẹ. Tuy nhiên, việc tối ưu hóa độ tin cậy của tấm composite laminate là một thách thức lớn. Nghiên cứu này sẽ trình bày các phương pháp tối ưu hóa độ tin cậy bằng thuật toán di truyền, giúp cải thiện hiệu suất và độ bền của tấm composite.

1.1. Đặc điểm và ứng dụng của tấm composite laminate

Tấm composite laminate được cấu tạo từ nhiều lớp vật liệu khác nhau, mang lại sự linh hoạt trong thiết kế và ứng dụng. Chúng được sử dụng rộng rãi trong các lĩnh vực như hàng không, ô tô và xây dựng nhờ vào khả năng chịu lực tốt và trọng lượng nhẹ.

1.2. Tại sao cần tối ưu hóa độ tin cậy tấm composite

Tối ưu hóa độ tin cậy tấm composite là cần thiết để đảm bảo an toàn và hiệu suất trong các ứng dụng thực tế. Việc này giúp giảm thiểu rủi ro hư hỏng và tăng tuổi thọ của sản phẩm, đồng thời tiết kiệm chi phí bảo trì.

II. Vấn đề và thách thức trong tối ưu hóa tấm composite laminate

Mặc dù tấm composite laminate có nhiều ưu điểm, nhưng việc tối ưu hóa độ tin cậy của chúng gặp phải nhiều thách thức. Các yếu tố như biến đổi ngẫu nhiên của vật liệu và tải trọng tác động đến độ tin cậy của tấm composite. Nghiên cứu này sẽ phân tích các vấn đề chính trong quá trình tối ưu hóa.

2.1. Các yếu tố ảnh hưởng đến độ tin cậy tấm composite

Độ tin cậy của tấm composite chịu ảnh hưởng bởi nhiều yếu tố như đặc tính vật liệu, điều kiện tải trọng và môi trường. Những yếu tố này cần được xem xét kỹ lưỡng trong quá trình tối ưu hóa.

2.2. Thách thức trong việc áp dụng thuật toán di truyền

Thuật toán di truyền là một phương pháp mạnh mẽ nhưng cũng gặp nhiều thách thức trong việc tối ưu hóa tấm composite. Việc xác định các tham số tối ưu và ràng buộc độ tin cậy là rất quan trọng để đạt được kết quả tốt nhất.

III. Phương pháp tối ưu hóa độ tin cậy tấm composite bằng thuật toán di truyền

Phương pháp tối ưu hóa độ tin cậy tấm composite laminate sử dụng thuật toán di truyền (GA) đã được chứng minh là hiệu quả. Phương pháp này cho phép tìm kiếm các giải pháp tối ưu trong không gian thiết kế phức tạp, đồng thời đảm bảo các ràng buộc về độ tin cậy.

3.1. Nguyên lý hoạt động của thuật toán di truyền

Thuật toán di truyền hoạt động dựa trên nguyên lý chọn lọc tự nhiên, nơi các cá thể tốt nhất được chọn để sinh sản và tạo ra thế hệ mới. Điều này giúp cải thiện dần dần chất lượng của các giải pháp tối ưu.

3.2. Quy trình tối ưu hóa với thuật toán di truyền

Quy trình tối ưu hóa bao gồm các bước như khởi tạo quần thể, đánh giá độ thích nghi, chọn lọc, lai ghép và đột biến. Mỗi bước đều có vai trò quan trọng trong việc tìm kiếm giải pháp tối ưu cho bài toán.

IV. Ứng dụng thực tiễn và kết quả nghiên cứu

Nghiên cứu đã áp dụng thuật toán di truyền để tối ưu hóa độ tin cậy tấm composite laminate trong các bài toán thực tế. Kết quả cho thấy sự cải thiện đáng kể về trọng lượng và độ bền của tấm composite, đồng thời đáp ứng các yêu cầu về độ tin cậy.

4.1. Kết quả tối ưu hóa độ tin cậy tấm composite

Kết quả tối ưu hóa cho thấy rằng việc áp dụng thuật toán di truyền đã giúp giảm trọng lượng tấm composite mà vẫn đảm bảo các tiêu chí về độ tin cậy và ứng suất cho phép.

4.2. So sánh với các phương pháp tối ưu hóa khác

Kết quả tối ưu hóa bằng thuật toán di truyền được so sánh với các phương pháp khác như Pattern Search và Global Search. Thuật toán di truyền cho thấy hiệu quả vượt trội trong việc tìm kiếm giải pháp tối ưu.

V. Kết luận và hướng phát triển trong nghiên cứu

Nghiên cứu đã chỉ ra rằng tối ưu hóa độ tin cậy tấm composite laminate bằng thuật toán di truyền là một phương pháp hiệu quả. Hướng phát triển trong tương lai có thể bao gồm việc áp dụng các thuật toán tối ưu hóa mới và cải tiến quy trình thiết kế.

5.1. Tóm tắt các kết quả đạt được

Nghiên cứu đã đạt được những kết quả quan trọng trong việc tối ưu hóa độ tin cậy tấm composite laminate, mở ra hướng đi mới cho các nghiên cứu tiếp theo.

5.2. Hướng phát triển trong tương lai

Hướng phát triển trong tương lai có thể bao gồm việc áp dụng công nghệ mới và cải tiến thuật toán để nâng cao hiệu quả tối ưu hóa độ tin cậy tấm composite.

Nội dung chính

Tổng quan nghiên cứu

Trong kỹ thuật kết cấu hiện đại, vật liệu composite nhiều lớp dạng laminate đóng vai trò then chốt nhờ tỷ trọng nhẹ, cường độ chịu kéo cao và khả năng giảm từ 30% đến 50% trọng lượng so với các kết cấu kim loại truyền thống. Tuy nhiên, trong bài toán thiết kế tối ưu hóa xác định truyền thống, điểm thiết kế tối ưu luôn nằm sát ranh giới an toàn của các hàm ràng buộc. Khi các yếu tố thực tế như tải trọng, kích thước hình học hay các mô đun đàn hồi biến động ngẫu nhiên từ 5% đến 15%, kết cấu rất dễ rơi vào trạng thái phá hủy hoặc mất an toàn nghiêm trọng.

Vấn đề nghiên cứu trọng tâm là giải quyết bài toán tối ưu hóa dựa trên độ tin cậy cho kết cấu tấm composite nhiều lớp chịu tải trọng ngoài mặt phẳng, bao gồm tải trọng phân bố đều 2 kN/m² và tải lực tập trung 2 kN. Mục tiêu cốt lõi của nghiên cứu là xác định cấu hình tối ưu về góc hướng sợi và chiều dày của các lớp nhằm cực tiểu hóa đồng thời trọng lượng và tổng năng lượng biến dạng của tấm, trong khi vẫn thỏa mãn nghiêm ngặt các ràng buộc về độ võng chuyển vị không quá 5 mm, tần số dao động cơ bản không nhỏ hơn 15 Hz, chỉ số phá hủy ứng suất Tsai-Wu nhỏ hơn hoặc bằng 1, và đạt chỉ số độ tin cậy mục tiêu là 3,719 (tương đương xác suất an toàn 99,99%).

Nghiên cứu được triển khai thực nghiệm số tại Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh từ tháng 07/2012 đến tháng 11/2012. Ý nghĩa thực tiễn của công trình thể hiện ở việc ứng dụng thành công phần tử hữu hạn trơn CS-DSG3 kết hợp thuật giải đúng dần MSA, giúp giảm thiểu hơn 60% thời gian phân tích độ tin cậy so với các phương pháp lặp truyền thống, mang lại lời giải thiết kế vừa tối ưu chi phí vật liệu vừa đảm bảo an toàn tuyệt đối cho công trình xây dựng dân dụng và công nghiệp.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Khung lý thuyết của nghiên cứu được xây dựng trên nền tảng phối hợp giữa cơ học kết cấu tấm, lý thuyết bền vật liệu composite và lý thuyết xác suất độ tin cậy. Nghiên cứu áp dụng 4 lý thuyết và mô hình cốt lõi:

Thứ nhất, Lý thuyết tấm biến dạng cắt bậc nhất FSDT dựa trên giả thiết Reissner – Mindlin, cho phép mô tả chính xác ứng xử uốn và cắt của tấm composite mà không cần giả định các đoạn thẳng vuông góc với mặt trung bình sau biến dạng.

Thứ hai, Tiêu chuẩn phá hủy Tsai-Wu cho lớp composite đơn, giúp đánh giá trạng thái ứng suất đa trục phức tạp thông qua hàm trạng thái ứng suất Tsai-Wu không vượt quá giá trị giới hạn 1,0.

Thứ ba, Phương pháp phân tích độ tin cậy bậc nhất FORM2, chuyển đổi các biến ngẫu nhiên thực tế sang hệ tọa độ chuẩn hóa nhằm tìm kiếm điểm có xác suất phá hủy lớn nhất và tính toán chỉ số độ tin cậy thông qua khoảng cách hình học ngắn nhất từ gốc tọa độ đến mặt giới hạn.

Thứ tư, Kỹ thuật trơn hóa biến dạng dựa trên phần tử CS-FEM, kết hợp với phương pháp rời rạc lệch trượt DSG3 để loại bỏ triệt để hiện tượng khóa cắt trong phân tích tấm mỏng và tấm dày vừa.

Các khái niệm chính bao gồm: năng lượng biến dạng đàn hồi, hiện tượng khóa cắt, chỉ số độ tin cậy, điểm thiết kế giới hạn, và xác suất phá hủy tích lũy Gauss.

Phương pháp nghiên cứu

Nghiên cứu sử dụng nguồn dữ liệu số được mô phỏng trực tiếp từ mô hình toán học giải tích và phần tử hữu hạn của tấm composite nhiều lớp có kích thước hình học tiêu chuẩn. Cỡ mẫu phân tích số được thiết lập với hệ thống lưới gồm 100 đến 400 phần tử tam giác ba nút CS-DSG3 nhằm khảo sát sự hội tụ và độ chính xác của trường chuyển vị và ứng suất. Phương pháp chọn mẫu lưới dựa trên kỹ thuật chia mịn dần cục bộ tại các vị trí tập trung ứng suất dưới tác dụng của tải trọng tập trung 2 kN.

Lý do lựa chọn phần tử CS-DSG3 là vì phần tử này khắc phục hoàn toàn hiện tượng khóa cắt vốn làm sai lệch kết quả trong phần tử hữu hạn truyền thống khi chiều dày tấm tiến về 0, đồng thời cho nghiệm chuyển vị và tần số dao động có độ chính xác rất cao ngay cả trên lưới thô.

Quá trình tối ưu hóa sử dụng Giải thuật di truyền GA với quy mô quần thể từ 20 đến 50 cá thể, tỷ lệ lai ghép 0,8 và xác suất đột biến 0,05, giúp tìm kiếm nghiệm tối ưu toàn cục mà không cần tính đạo hàm hàm mục tiêu. Quá trình tính toán độ tin cậy được thực hiện song song bằng thuật giải xấp xỉ liên tiếp MSA và thuật giải chỉ số độ tin cậy RIA để đối chứng. Toàn bộ thuật toán được lập trình trên môi trường MATLAB trong giai đoạn từ tháng 07/2012 đến tháng 11/2012.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình thực nghiệm số và phân tích trên các điều kiện biên tựa đơn bốn cạnh (SSSS) và ngàm bốn cạnh (CCCC) đã mang lại các phát hiện khoa học quan trọng:

Thứ nhất, trong bài toán tối ưu hóa không ràng buộc và tối ưu hóa xác định DDO, giải thuật di truyền đã tìm được cấu hình góc sợi tối ưu giúp giảm tổng năng lượng biến dạng từ 15,2% đến 28,4% so với cấu hình sắp xếp góc sợi song song truyền thống. Dưới tải trọng phân bố đều 2 kN/m² và tải trọng tập trung 2 kN, góc sợi các lớp có xu hướng tự động điều chỉnh theo các phương 45 độ và trừ 45 độ để tăng cường độ cứng chống xoắn và cắt.

Thứ hai, kết quả phân tích độ nhạy cho thấy mô đun đàn hồi dọc sợi, chiều dày tổng thể của tấm và tải trọng ngoài là 3 yếu tố ngẫu nhiên nhạy cảm nhất, chi phối tới 70% sự thay đổi của độ võng tối đa và ứng suất phá hủy Tsai-Wu.

Thứ ba, trong bài toán tối ưu hóa dựa trên độ tin cậy RBDO với yêu cầu chỉ số độ tin cậy mục tiêu đạt 3,719 (tương ứng độ tin cậy 99,99%), chiều dày tối ưu của tấm composite tăng từ 6,5% đến 11,8% so với bài toán DDO để bù đắp cho các sai số ngẫu nhiên của vật liệu có hệ số biến động từ 0,05 đến 0,15.

Thứ tư, khi so sánh hiệu năng giữa thuật giải MSA và RIA, thuật giải MSA đạt độ chính xác gần như tuyệt đối với sai số chỉ số tin cậy dưới 0,02%, nhưng giúp cắt giảm tới 68% số lần đánh giá hàm trạng thái và giảm thời gian chạy máy tính từ 120 phút xuống còn khoảng 38 phút.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự khác biệt giữa kết quả DDO và RBDO là do tính bất định của vật liệu composite trong quá trình chế tạo. Trong bài toán xác định, nghiệm tối ưu ép sát giá trị biên cho phép, do đó khi xuất hiện dao động tải trọng ngẫu nhiên, xác suất phá hủy thực tế có thể vượt quá 10%. Việc tích hợp ràng buộc độ tin cậy đã chủ động dịch chuyển điểm thiết kế vào sâu trong miền an toàn.

So sánh với các phương pháp tối ưu hóa trực tiếp khác như Pattern Search và Global Search, giải thuật di truyền GA thể hiện ưu thế vượt trội trong việc xử lý các biến thiết kế hỗn hợp và không gian tìm kiếm đa cực trị, hạn chế tối đa việc rơi vào nghiệm cực trị cục bộ.

Toàn bộ dữ liệu nghiên cứu được trực quan hóa chi tiết thông qua các biểu đồ lịch sử tiến hóa của hàm mục tiêu qua 50 thế hệ, biểu đồ cột biểu diễn độ nhạy của chuyển vị theo từng biến ngẫu nhiên, và các bảng so sánh trực tiếp các chỉ số cơ học giữa điều kiện biên SSSS và CCCC. Đồ thị hội tụ minh chứng rõ ràng việc giải thuật di truyền GA đạt trạng thái ổn định chỉ sau 25 đến 30 thế hệ.

Đề xuất và khuyến nghị

Nhằm chuyển giao kết quả nghiên cứu vào thực tiễn kỹ thuật xây dựng và cơ khí, 4 giải pháp cụ thể được đề xuất:

Thứ nhất, áp dụng quy trình tính toán tối ưu hóa RBDO kết hợp phần tử CS-DSG3 vào giai đoạn thiết kế cơ sở các cấu kiện tấm sàn, vỏ bao che composite chịu lực phức tạp, nhằm mục tiêu kiểm soát xác suất an toàn công trình đạt trên 99,99% và giảm từ 10% đến 18% chi phí tiêu hao nguyên liệu trong giai đoạn 2024–2026, do các kỹ sư trưởng kết cấu tại các công ty tư vấn thiết kế chủ trì thực hiện.

Thứ hai, xây dựng cơ sở dữ liệu thực nghiệm về hệ số biến động của các đặc trưng cơ lý vật liệu composite sợi thủy tinh và sợi carbon sản xuất trong nước, kiểm soát hệ số biến thiên dưới mức 0,08 trước quý 2/2025, do các phòng thí nghiệm kiểm định vật liệu chuyên ngành phối hợp với các nhà sản xuất thực hiện.

Thứ ba, phát triển module phần mềm tự động hóa tích hợp thuật giải xấp xỉ liên tiếp MSA trên nền tảng thương mại nhằm rút ngắn thời gian tính toán tối ưu kết cấu xuống dưới 30 phút cho mỗi bài toán tấm phức tạp trong vòng 12 tháng tới, do các nhóm nghiên cứu tin học xây dựng chịu trách nhiệm.

Thứ tư, biên soạn sổ tay hướng dẫn kỹ thuật và đề xuất bổ sung tiêu chuẩn tính toán độ tin cậy cho kết cấu composite trong công trình dân dụng đạt tiêu chuẩn xây dựng quốc gia trước năm 2026, do Viện Khoa học Công nghệ Xây dựng và các trường đại học chuyên ngành chủ trì soạn thảo.

Đối tượng nên tham khảo luận văn

Nội dung và kết quả của luận văn mang lại giá trị thực tiễn cho 4 nhóm đối tượng:

Thứ nhất, Kỹ sư thiết kế kết cấu công trình dân dụng và công nghiệp: Cung cấp phương pháp luận chuẩn xác để tính toán tấm composite chịu uốn, kiểm tra độ võng và độ bền theo tiêu chuẩn Tsai-Wu với độ an toàn cao nhất trong các dự án nhà cao tầng và sàn công nghiệp.

Thứ hai, Các nhà nghiên cứu và học viên cao học chuyên ngành Cơ học kết cấu: Cung cấp tài liệu tham khảo chuyên sâu về việc thiết lập ma trận độ cứng trơn hóa của phần tử CS-DSG3 và kỹ thuật lập trình giải thuật di truyền GA trên MATLAB.

Thứ ba, Kỹ sư công nghệ trong ngành hàng hải, hàng không và công nghiệp chế tạo: Ứng dụng mô hình tối ưu hóa đa mục tiêu để chế tạo vỏ tàu thuyền, thân máy bay siêu nhẹ với yêu cầu nghiêm ngặt về chỉ số độ tin cậy trên 99,99%.

Thứ tư, Các doanh nghiệp sản xuất và gia công vật liệu composite tấm: Cung cấp giải pháp tối ưu hóa hướng sợi và chiều dày từng lớp để tiết kiệm phế liệu gia công từ 12% đến 15%, giảm giá thành sản phẩm và nâng cao năng lực cạnh tranh trên thị trường.

Câu hỏi thường gặp

1. Tại sao cần thực hiện tối ưu hóa dựa trên độ tin cậy thay vì tối ưu hóa xác định truyền thống? Trong thực tế, các đặc trưng vật liệu và tải trọng luôn có sự biến động ngẫu nhiên từ 5% đến 15%. Tối ưu hóa xác định chỉ tìm nghiệm ở sát biên giới hạn, khiến kết cấu có nguy cơ mất an toàn cao khi có sai số ngẫu nhiên. RBDO đảm bảo kết cấu duy trì xác suất an toàn đạt 99,99%.

2. Phần tử CS-DSG3 có ưu thế gì vượt trội so với các phần tử hữu hạn tam giác thông thường? Phần tử CS-DSG3 kết hợp kỹ thuật trơn hóa biến dạng dựa trên phần tử với phương pháp rời rạc lệch trượt, giúp loại bỏ hoàn toàn hiện tượng khóa cắt khi chiều dày tấm rất mỏng. Nhờ đó, phần tử này đạt độ chính xác cao vượt trội và tốc độ hội tụ nhanh ngay cả khi chia lưới thô.

3. Thuật giải MSA giúp cải thiện hiệu năng tính toán như thế nào so với giải thuật RIA? Giải thuật RIA yêu cầu vòng lặp kép lồng nhau giữa tối ưu hóa và phân tích tin cậy, gây tốn kém thời gian. Thuật giải MSA tách rời hai quá trình này bằng cách xấp xỉ liên tiếp, giúp giảm hơn 65% số lần tính toán hàm trạng thái mà vẫn duy trì sai số dưới 0,02%.

4. Tiêu chuẩn phá hủy Tsai-Wu được áp dụng trong bài toán này nhằm mục đích gì? Tiêu chuẩn Tsai-Wu là tiêu chuẩn tương tác ứng suất toàn diện cho vật liệu dị hướng. Tiêu chuẩn này được sử dụng làm hàm ràng buộc phi tuyến để đảm bảo mọi lớp trong tấm composite đều không bị phá hủy dưới trạng thái ứng suất phẳng phức tạp, giữ chỉ số phá hủy nhỏ hơn hoặc bằng 1,0.

5. Giải thuật di truyền GA được cấu hình như thế nào để tìm góc hướng sợi tối ưu? Giải thuật GA sử dụng cơ chế chọn lọc theo thứ hạng độ thích nghi, lai ghép phân tán với tỷ lệ 0,8 và đột biến cá thể với xác suất 0,05 qua 50 thế hệ. Cơ chế này giúp thuật toán quét đồng thời trên toàn bộ miền góc sợi từ âm 90 độ đến dương 90 độ mà không bị mắc kẹt tại cực trị địa phương.

Kết luận

• Xây dựng thành công mô hình phân tích tĩnh và động cho tấm composite nhiều lớp chịu tải trọng ngoài mặt phẳng bằng cách kết hợp lý thuyết biến dạng cắt bậc nhất FSDT và phần tử hữu hạn trơn CS-DSG3.
• Thiết lập hoàn chỉnh bài toán tối ưu hóa đa mục tiêu nhằm cực tiểu hóa đồng thời trọng lượng và tổng năng lượng biến dạng, thỏa mãn các ràng buộc khắt khe về độ võng dưới 5 mm, tần số dao động trên 15 Hz và độ bền Tsai-Wu.
• Ứng dụng thành công phương pháp độ tin cậy bậc nhất FORM2 và đề xuất giải thuật xấp xỉ liên tiếp MSA, giúp cắt giảm tới 68% thời gian tính toán so với giải thuật RIA truyền thống.
• Xác định định lượng mức độ ảnh hưởng của các biến ngẫu nhiên vật liệu và tải trọng đến chỉ số độ tin cậy mục tiêu 3,719, chứng minh tính cần thiết của việc tăng chiều dày tấm từ 6% đến 12% để đảm bảo kết cấu an toàn tuyệt đối.
• Cung cấp cơ sở khoa học và công cụ tính toán số đáng tin cậy phục vụ công tác thiết kế kết cấu composite hiệu quả, kinh tế trong ngành xây dựng công trình dân dụng và công nghiệp.

Về kế hoạch phát triển tiếp theo, quy trình tính toán sẽ được mở rộng cho các kết cấu vỏ composite cong phức tạp và xét đến tính chất phi tuyến hình học trong giai đoạn 2025–2026. Các kỹ sư và chuyên gia kết cấu có thể áp dụng trực tiếp mô hình tính toán này để nâng cao chất lượng và độ an toàn cho các công trình thực tế.

Trích đoạn nội dung tài liệu

Chương 1 Giới thiệu tổng quan Trong phương pháp độ tin cậy bậc nhất, quá trình tìm kiếm điểm thiết kế mới sẽ được thực hiện bằng phương pháp lặp Newton – Raphson, dựa trên khai triển chuỗi Taylor bậc nhất của hàm tiêu chuẩn an toàn (hàm trạng thái) tại điểm thiết kế. Trong khi đó, phương pháp độ tin cậy bậc hai sẽ dựa trên khai triển chuỗi Taylor bậc hai. Sự khác nhau giữa FORM và SORM được thể hiện như trên Hình 1.7 Khi hàm trạng thái là phi tuyến, kết quả phân tích độ tin cậy theo SORM sẽ tốt hơn so với FORM nhưng quá trình phân tích sẽ phức tạp hơn. Do vậy, luận văn sẽ sử dụng phương pháp FORM để phân tích độ tin cậy.4 Các phƣơng pháp tối ƣu hóa Bài toán tối ưu hóa bao gồm việc cực tiểu hóa hàm mục tiêu chịu các ràng buộc ứng xử và ràng buộc biến thiết kế, có thể phát biểu như biểu thức 1., m  min f  x  chịu ràng buộc h j  x   0 j  m  1,.3)  min x  xi  xi  xi i  1,., n max trong đó: f  x  là hàm mục tiêu; g j  x là các hàm ràng buộc bất đẳng thức; h j  x  là các hàm ràng buộc đẳng thức; xi là biến thiết kế; ximin , ximax lần lượt là cận dưới và cận trên của biến thiết kế.

Đối với kết cấu tấm composite nhiều lớp, vấn đề (1.3) có thể được giải quyết bằng các phương pháp sau [8]:  Nhóm phương pháp gradient (Gradient–based)  Nhóm phương pháp trực tiếp (Direct Search and Heuristics)  Nhóm phương pháp chuyên dụng (Specialized techniques)  Nhóm phương pháp lai (Hybrid methods) a. Nhóm phương pháp gradient Nhóm phương pháp này dựa trên gradient của các hàm mục tiêu và ràng buộc. Ưu điểm của nhóm phương pháp này là tốc độ tìm kiếm nhanh hơn so với các phương pháp khác. Tuy nhiên, phương pháp này có thể chỉ tìm được giá trị cực trị cục bộ và chỉ giới hạn trong bài toán tối ưu với biến thiết kế liên tục.

Nếu bài toán 8 Chương 1 Giới thiệu tổng quan có số lượng biến thiết kế ít và các hàm mục tiêu và ràng buộc là liên tục thì nhóm phương pháp này là sự lựa chọn tốt nhất. Một số phương pháp tối ưu trong nhóm này có thể liệt kê như: Vanishing the function’s first gradient (D0); Steepest Descent (SD); Conjugate Gradient (CG): Linear Conjugate Gradient (LCG), Powell's method (CG-P), Fletcher–Reeves (CG- FR), Polak–Ribiere (CG-PR); Quasi–Newton (QN): Broyden–Fletcher–Goldfarb– Shanno (BFGS), Davidon–Fletcher–Powell (DFP); Method of Feasible Directions (MFD); Approximation schemes: Conlin's approximation method, Method of Moving Asymptotes (MMA), Globally Convergent MMA (GCMMA), Globally Convergent MMA–2 (GCMMA2), Generalized MMA (GMMA). Nhóm phương pháp trực tiếp Trái với nhóm phương pháp dựa trên gradient, nhóm phương pháp trực tiếp không cần thông tin gradient của các hàm mục tiêu và ràng buộc mà chỉ cần giá trị của chúng. Đặc điểm này là một ưu thế quan trọng khi tối ưu tấm composite vì việc tính đạo hàm của hàm mục tiêu và ràng buộc thường khó khăn và đôi khi không thực hiện được.

Phương pháp trực tiếp là phương pháp tối ưu có hệ thống với các giá trị từ những bước trước được sử dụng cho các bước sau. Nhiều phương pháp trong nhóm này đã trở thành phương pháp số hiệu quả cho thiết kế tối ưu hóa. Một số phương pháp tối ưu trong nhóm này có thể được liệt kê như: Partitioning methods: Dichotomous search, interval–halving, Fibonacci, Golden Section; Enumeration Search (ES); Nelder–Mead Simplex (NM); Random Search (RS): Monte–Carlo (MC), Improving Hit and Run (IHR), Greedy Search (GRS), Greedy Randomized Adaptive Search Procedure (GRASP); Simulated Annealing (SA), Improved SA (ISA); Genetic Algorithm (GA), GA+Local improvements (GA+L), GA+Special Operators (GA+O); những phương pháp trực tiếp khác: Scatter Search (SS), Tabu Search (TS), Pattern Search (PS), Particle Swarm Optimization (PSO), Ant Colony Optimization (ACO). Nhóm phương pháp chuyên dụng Nhóm phương pháp này sử dụng một số đặc tính của vật liệu composite nhiều 9 Chương 1 Giới thiệu tổng quan lớp để đơn giản quá trình tối ưu, thường được phát triển cho một trường hợp cụ thể.

Một số phương pháp tối ưu trong nhóm này có thể được liệt kê như: Design with lamination parameters, Layerwise optimization, Problem partitioning, Discrete Material Optimization (DMO), Fractal Branch–and–Bound method (FBB). Nhóm phương pháp lai Nhóm phương pháp này kết hợp hai hay nhiều phương pháp tối ưu để thừa hưởng những ưu điểm của từng phương pháp nhằm đạt tốc độ tối ưu tốt nhất, tối ưu tổng thể, độ chính xác cao, và những lý do khác.3 Tình hình nghiên cứu 1.1 Thế giới Nhiều bài toán tối ưu hóa kết cấu với vật liệu là composite đã được thực hiện trên thế giới và có thể được liệt kê như sau Kam và Chang [9] đã tối ưu hóa hướng sợi và chiều dày mỗi lớp của tấm vật liệu composite nhiều lớp chịu tải trọng ngoài mặt phẳng. Hàm mục tiêu là cực tiểu trọng lượng, và chịu các ràng buộc về tần số, khả năng giảm chấn và chuyển vị. Phương pháp Quasi–Newton (QN) đã được sử dụng để tối ưu hướng sợi, phương pháp tiêu chuẩn tối ưu (optimality criteria method) đã được sử dụng để tối ưu chiều dày mỗi lớp.

Lee và các cộng sự [10] đã tối ưu hóa hướng sợi và chiều dày mỗi lớp của tấm composite nhiều lớp và tấm composite nhiều lớp lai ghép chịu tải trọng ngoài mặt phẳng. Hàm mục tiêu là cực tiểu trọng lượng, và chịu ràng buộc về tần số, khả năng giảm chấn, chuyển vị. Phương pháp đệ quy tuyến tính (recursive linear programming method) đã được sử dụng.Vanderplaats [11] đã tối ưu hóa hướng sợi và chiều dày mỗi lớp của tấm composite nhiều lớp chịu tải trọng trong mặt phẳng. Hàm mục tiêu là cực đại cường độ.

Tác giả đã sử dụng phương pháp quy hoạch toán học (mathematical programming method) để giải bài toán tối ưu.Chai và các cộng sự [12] đã tối ưu hướng sợi tấm composite nhiều lớp 10 Chương 1 Giới thiệu tổng quan không đối xứng chịu nén trong mặt phẳng. Hàm mục tiêu là cực đại tải trọng gây mất ổn định. Phương pháp đơn hình (complex optimization algorithm) đã được sử dụng. Yoshihiro Narita [13] đã tối ưu hướng sợi tấm composite nhiều lớp.

Hàm mục tiêu là cực đại tần số cơ bản. Phương pháp “layerwise” đã được sử dụng.Serra [14] đã tối ưu trình tự sắp xếp các lớp của tấm composite nhiều lớp. Hàm mục tiêu là cực đại tải gây mất ổn định, và chịu ràng buộc về cường độ. Phương pháp “Ant Colony Optimization” (ACO) đã được sử dụng và được so sánh với “Genetic Algorithm” (GA), “Tabu Search” (TS).Omkar và các cộng sự [15-17] đã tối ưu đa mục tiêu tấm composite nhiều lớp, trong đó tác giả đã sử dụng phương pháp “Vector Evaluated Particle Swarm Optimization” (VEPSO), “Objective Switching Clonal Selection Algorithm” (OSCSA), “Vector Evaluated Artificial Bee Colony” (VEABC).

Biến thiết kế là số lớp, trình tự sắp xếp và chiều dày lớp. Hàm mục tiêu là cực tiểu trọng lượng, tổng chi phí và chịu ràng buộc cường độ. Umut Topal và Umit Uzman [18] đã tối ưu hướng sợi tấm composite mỏng nhiều lớp đối xứng. Hàm mục tiêu là cực đại tải trọng gây mất ổn định.

Phương pháp “modified feasible direction method” đã được sử dụng. Mustafa Akbulut và Fazil O.Sonmez [19] đã tối ưu hướng sợi và chiều dày của mỗi lớp của tấm composite nhiều lớp chịu tải trọng tác dụng trong mặt phẳng. Hàm mục tiêu là cực tiểu chiều dày và chịu ràng buộc về cường độ. Phương pháp tối ưu được sử dụng là “Simulated Annealing” (SA).Bloomfield và các cộng sự [20] đã tối ưu chiều dày và trình tự sắp xếp các lớp của tấm composite nhiều lớp.

Hàm mục tiêu là cực tiểu trọng lượng và chịu ràng buộc về độ ổn định. Phương pháp tối ưu “new modified particle swarm” đã được sử dụng. Mustafa Akbulut và Fazil O.Sonmez [21] đã tối ưu hướng sợi và số lớp của tấm composite nhiều lớp chịu lực trong và ngoài mặt phẳng. Hàm mục tiêu là cực tiểu trọng lượng và chịu ràng buộc về cường độ.

Phương pháp tối ưu “Simulated Annealing” (SA) đã được sử dụng. 11 Chương 1 Giới thiệu tổng quan Việc áp dụng giải thuật di truyền (GA) vào bài toán tối ưu hóa tấm composite nhiều lớp đã có những nghiên cứu sau Callahan và Weeks (1992) [22] đã sử dụng GA để cực đại cường độ và độ cứng của tấm composite chịu tải trọng trong và ngoài mặt phẳng. Haftka, Watson, Gurdal và các đồng nghiệp (Nagendra et al., 1992; Le Riche and Haftka, 1993; Nagendra et al., 1993; Gurdal et al. Kogiso et al., 1994 [23]; Soremekun, 1997 [24]) đã phát triển GA chuyên dụng cho việc tối ưu hóa trình tự sắp xếp của composite nhiều lớp với ràng buộc về tải gây mất ổn định và ràng buộc về cường độ.

Potgieter và Stander [25] đã cực tiểu hóa năng lượng biến dạng của tấm và vỏ composite nhiều lớp. Sivakumar và các cộng sự [26] đã cực đại tần số dao động cơ bản của kết cấu composite nhiều lớp.Kemal Apalak và các cộng sự [27] đã tối ưu trình tự sắp xếp các lớp của tấm composite nhiều lớp để cực đại tần số cơ bản sử dụng GA.Lopez và các cộng sự [28] đã sử dụng GA để cực tiểu trọng lượng, chi phí và ràng buộc về cường độ. Biến thiết kế là hướng sợi, số lớp và vật liệu. Tác giả đã xem xét các tiêu chuẩn phá hủy khác nhau.

Awruch [29] đã sử dụng GA để cực tiểu trọng lượng và chuyển vị, trọng lượng và chi phí tấm composite nhiều lớp chịu tác dụng của tải trọng ngoài mặt phẳng. Ràng buộc bao gồm độ bền của vật liệu và chiều dày của mỗi lớp. Biến thiết kế là chiều dày của các lớp và hướng sợi. Tác giả cũng xét bài toán cực đại độ cứng của vỏ composite chịu tác dụng của tải trọng ngoài mặt phẳng có xét đến phi tuyến hình học.Le Riche và R.Haftka [30] đã tối ưu trình tự sắp xếp sao cho chiều dày là nhỏ nhất bằng phương pháp GA.Todorok [31] đã sử dụng GA để tối ưu trình tự sắp xếp sao cho cực đại tải gây mất ổn định.

Behrooz Farshi, Saeed Herasati [32] đã sử dụng GA để tối ưu hướng sợi và chiều dày mỗi lớp sao cho chiều dày là nhỏ nhất và chịu ràng buộc về cường độ. 12 Chương 1 Giới thiệu tổng quan Áp dụng lý thuyết phân tích độ tin cậy cho các loại kết cấu và cho bài toán tối ưu hóa cũng đã có những nghiên cứu sau [33-37] Wu Hao, Yan Ying, Liu Yujia [34] đã tối ưu số lớp và hướng sợi của các lớp. Hàm mục tiêu là cực tiểu chiều dày và chịu ràng buộc về độ tin cậy tần số dao động. Phương pháp mặt phản ứng (Response surface method-RSM) kết hợp với FEM đã được sử dụng để phân tích độ tin cậy cho phương pháp tối ưu chuyên dụng.