Tổng quan nghiên cứu

Trong kỹ thuật công trình hiện đại, hơn 70% các kết cấu dạng tấm vỏ trong công trình biển, bồn chứa chất lỏng, đập thủy điện và buồng cách âm chịu tác động trực tiếp từ hiện tượng tương tác giữa kết cấu rắn và môi trường lưu chất (Fluid-Structure Interaction - FSI). Khi phân tích động lực học kết cấu tấm đàn hồi, phương pháp phần tử hữu hạn (FEM) truyền thống thường mắc phải hiện tượng quá cứng và hiện tượng khóa cắt nghiêm trọng, gây ra sai số dao động riêng từ 8% đến 15% khi chiều dày tấm giảm dần. Nhằm giải quyết bài toán cấp thiết này, luận văn thạc sĩ chuyên ngành Kỹ thuật Xây dựng Công trình Dân dụng và Công nghiệp của tác giả Nguyễn Quốc Toản tại Trường Đại học Bách Khoa - Đại học Quốc gia Thành phố Hồ Chí Minh đã đề xuất phương pháp kết hợp CS-MIN3/FS-FEM để phân tích dao động tự do của tấm Mindlin trong không gian 3 chiều.

Nghiên cứu được triển khai thực hiện trong giai đoạn từ tháng 2 năm 2012 đến tháng 9 năm 2012 dưới sự hướng dẫn khoa học của Tiến sĩ Lương Văn Hải và Tiến sĩ Nguyễn Thời Trung. Mục tiêu cụ thể của luận văn là thiết lập hệ phương trình ma trận rời rạc ghép nối tương tác rắn - lưu chất, lập trình thuật toán số trên môi trường Matlab và đánh giá chính xác tần số dao động tự do cũng như đáp ứng dao động cưỡng bức. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn khi giúp giảm thiểu sai số tính toán dao động xuống dưới mức 1,5%, tăng tốc độ hội tụ lưới lên gấp 2 lần và tiết kiệm hơn 20% chi phí thử nghiệm mô hình thực tế trong thiết kế công trình.

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 luận văn được xây dựng trên nền tảng cơ học kết cấu tấm Reissner-Mindlin và lý thuyết truyền sóng âm học trong lưu chất. Lý thuyết tấm Mindlin mở rộng phạm vi ứng dụng cho cả tấm mỏng và tấm dày trung bình bằng cách kể đến ảnh hưởng của biến dạng trượt ngang và quán tính quay với hệ số điều chỉnh cắt k = 5/6, loại bỏ hạn chế bỏ qua biến dạng trượt của lý thuyết cổ điển Kirchhoff. Nhằm khắc phục triệt để hiện tượng khóa cắt, phương pháp tấm Mindlin 3 nút trơn hóa dựa trên ô con (CS-MIN3) được áp dụng bằng cách chia mỗi phần tử tam giác thành 3 tam giác con qua trọng tâm và làm trơn trường biến dạng uốn, cắt.

Đối với miền lưu chất 3 chiều, lý thuyết truyền sóng âm trong môi trường không nhớt, không xoáy được mô hình hóa bằng phương pháp phần tử hữu hạn trơn dựa trên mặt (FS-FEM) sử dụng phần tử tứ diện 4 nút. Hai môi trường vật lý được liên kết chặt chẽ thông qua ma trận tương tác bề mặt H, đảm bảo tính liên tục của chuyển vị pháp tuyến và cân bằng áp suất động tại biên tiếp xúc.

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

Nghiên cứu sử dụng phương pháp mô phỏng số thông qua lập trình thuật toán trên phần mềm Matlab, đồng thời đối chiếu trực tiếp với nghiệm từ phần mềm thương mại COMSOL Multiphysics và các phương pháp số MITC4/FEM-H8, MIN3/FEM-T4. Cỡ mẫu khảo sát bao gồm bộ dữ liệu từ 3 mô hình kết cấu 3 chiều điển hình: vách phòng cách âm hình chữ nhật kích thước 3m x 8m x 3m với 1848 bậc tự do, cấu kiện tấm ngàm 4 cạnh trên mặt nước với 3328 bậc tự do, và tấm tròn tương tác bể chứa khí với 3685 bậc tự do.

Phương pháp chọn mẫu có chủ đích được áp dụng nhằm bao quát các dạng hình học phẳng, hình học cong và điều kiện biên tiếp xúc đa dạng trong thực tế công trình. Lý do lựa chọn phương pháp phân tích số CS-MIN3/FS-FEM là khả năng tạo hiệu ứng làm mềm độ cứng giả tạo của mô hình phần tử hữu hạn chuẩn, cho phép sử dụng lưới phần tử tam giác và tứ diện đơn giản nhưng vẫn đạt độ chính xác cao tương đương các phần tử bậc cao. Quá trình nghiên cứu được triển khai chặt chẽ từ ngày 6 tháng 2 năm 2012 và hoàn thành bảo vệ ngày 28 tháng 9 năm 2012.

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

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

Nghiên cứu mang lại 3 phát hiện khoa học quan trọng thông qua các phép thử số nghiệm:

Thứ nhất, phương pháp kết hợp CS-MIN3/FS-FEM triệt tiêu hoàn toàn hiện tượng quá cứng của phần tử hữu hạn chuẩn. Tại bài toán vách phòng kích thước 3m x 8m x 3m, tần số dao động riêng mode 1 đạt giá trị 7,12 Hz với sai số chỉ 0,82% so với nghiệm chuẩn COMSOL, trong khi phương pháp MIN3/FEM-T4 có mức sai lệch lên tới 6,45%.

Thứ hai, tốc độ hội tụ của CS-MIN3/FS-FEM vượt trội trên toàn bộ 9 mode dao động đầu tiên. Ở bài toán tấm đặt trên mặt nước, tần số mode 1 hội tụ chính xác tại giá trị 80,45 Hz với số lượng phần tử giảm 35% so với phương pháp MITC4/FEM-H8. Tại bài toán tấm tròn đáy bể khí, tần số mode 1 hội tụ ổn định ở mức 25,18 Hz ngay trên lưới thô.

Thứ ba, trường áp suất âm học trong miền lưu chất giảm dần phi tuyến khi càng cách xa mặt tiếp xúc, phản ánh chính xác bản chất vật lý truyền sóng với độ tương quan chuyển vị và áp suất vượt trên 98%.

Thảo luận kết quả

Nguyên nhân tạo nên tính ưu việt của phương pháp là nhờ kỹ thuật làm trơn biến dạng trên các ô con tam giác và trơn hóa gradient áp suất trên các mặt tứ diện. Kỹ thuật này giúp khôi phục độ mềm tự nhiên của kết cấu và môi trường âm học, loại bỏ năng lượng biến dạng giả tạo sinh ra do lưới chia thô. Khi so sánh với các nghiên cứu 2 chiều trước đây của Phùng Văn Phúc (2011) và Nguyễn Ngọc Nhân (2011), nghiên cứu này đã tạo bước tiến vượt bậc khi giải quyết trọn vẹn bài toán ghép nối 3 chiều với biên cong phức tạp.

Các kết quả có thể được trình bày trực quan thông qua bảng so sánh sai số tần số riêng giữa 4 phương pháp số và đồ thị đường cong đáp ứng biên độ chuyển vị theo tần số kích thích. Đồ thị cho thấy rõ hiện tượng quán tính gia tăng của lưu chất làm giảm biên độ cộng hưởng dao động cưỡng bức từ 15% đến 25% so với kết cấu dao động trong không khí tự do.

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

Thứ nhất, ứng dụng thuật toán CS-MIN3/FS-FEM vào quy trình tính toán thiết kế kháng chấn cho các bồn chứa chất lỏng và bể nước mái giảm chấn trong các tòa nhà cao tầng. Mục tiêu hướng đến là giảm tối thiểu 12% sai số ứng suất động học trong lộ trình 12 tháng, do các kỹ sư kết cấu tại các viện nghiên cứu và doanh nghiệp tư vấn xây dựng triển khai.

Thứ hai, tích hợp mô-đun phần tử hữu hạn trơn vào các phần mềm mô phỏng thương mại hoặc xây dựng các gói thư viện mã nguồn mở chuyên dụng trên nền tảng Matlab và Python. Nhóm phát triển phần mềm kỹ thuật số cần chủ trì thực hiện nhằm rút ngắn 30% thời gian tạo lưới 3D phức tạp trước quý 4 năm 2026.

Thứ ba, mở rộng mô hình tính toán cho các kết cấu tấm composite nhiều lớp và vật liệu cơ tính biến thiên (FGM) tương tác với chất lỏng có độ nhớt. Mục tiêu là nâng cao độ chính xác dự báo kiểm soát tiếng ồn trong cabin máy bay và khoang xe hơi lên trên 95% trong thời gian 18 tháng, do các nhóm nghiên cứu liên ngành hàng không và cơ học đảm nhiệm.

Thứ tư, ban hành hướng dẫn kỹ thuật chuyên ngành về phân tích tương tác kết cấu - lưu chất cho công trình đập thủy điện, giàn khoan dầu khí và cảng biển. Các hiệp hội kỹ thuật xây dựng và cơ quan quản lý cần phối hợp xây dựng để giúp tiết kiệm khoảng 15% chi phí gia cường kết cấu trước năm 2027.

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

Thứ nhất, học viên cao học và nghiên cứu sinh chuyên ngành Cơ học kỹ thuật, Xây dựng dân dụng và Công nghiệp, Công trình thủy. Đây là tài liệu tham khảo mẫu mực về kỹ thuật trơn hóa S-FEM và phương pháp thiết lập ma trận ghép nối trong bài toán tương tác rắn - lưu chất 3 chiều.

Thứ hai, kỹ sư kết cấu công trình biển, bồn chứa xăng dầu và hệ thống bể giảm chấn chất lỏng. Kỹ sư có thể ứng dụng trực tiếp quy trình tính toán để xác định chính xác tần số dao động riêng và hiện tượng cộng hưởng, đảm bảo an toàn chịu lực công trình với độ tin cậy trên 95%.

Thứ ba, chuyên gia âm học công trình và kỹ sư thiết kế phương tiện cơ giới đường bộ, đường thủy, hàng không. Luận văn cung cấp case study cụ thể về buồng cách âm và tấm vách mỏng, hỗ trợ tối ưu hóa thiết kế giảm thiểu mức ồn từ 3 đến 5 dB trong không gian kín.

Thứ tư, giảng viên đại học và các chuyên gia phát triển phần mềm mô phỏng CAE. Tài liệu cung cấp toàn bộ công thức giải tích rời rạc và thuật toán tối ưu, rất hữu ích để xây dựng giáo trình chuyên đề cơ học tính toán nâng cao hoặc phát triển công cụ phần mềm số mới.

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

Phương pháp CS-MIN3 giải quyết hiện tượng khóa cắt trong tấm mỏng như thế nào? Phương pháp CS-MIN3 chia mỗi phần tử tam giác thành 3 tam giác con và áp dụng kỹ thuật trơn hóa biến dạng trên từng ô con. Nhờ đó, trường biến dạng uốn và cắt được tính trung bình trên miền trơn, giúp loại bỏ hoàn toàn sự quá cứng khi tỷ lệ chiều dày trên nhịp tấm giảm dưới mức 0,001 mà không cần dùng hàm nội suy bậc cao phức tạp.

Vì sao cần kết hợp đồng thời CS-MIN3 cho tấm và FS-FEM cho lưu chất? Mô hình phần tử hữu hạn tiêu chuẩn bị quá cứng ở cả miền rắn lẫn miền âm học 3D. Việc kết hợp CS-MIN3 và FS-FEM tạo ra hiệu ứng làm mềm đồng bộ cho toàn bộ hệ thống ghép nối, giúp cải thiện độ chính xác tần số dao động riêng từ 5% đến hơn 10% ngay trên các lưới phần tử thô, giảm đáng kể tài nguyên tính toán.

Sự tương tác của lưu chất ảnh hưởng như thế nào đến tần số dao động riêng của tấm? Khi tiếp xúc với môi trường lưu chất, áp suất động và quán tính gia tăng từ lưu chất tác động trực tiếp lên bề mặt tấm làm giảm tần số dao động riêng so với trạng thái dao động trong chân không. Trong các ví dụ số của luận văn, tần số mode 1 của tấm khi tương tác chất lỏng giảm từ 20% đến 40% tùy thuộc mật độ môi trường.

Phương pháp CS-MIN3/FS-FEM có ưu thế gì về tốc độ tính toán so với COMSOL? Phương pháp CS-MIN3/FS-FEM sử dụng các phần tử tam giác 3 nút và tứ diện 4 nút dạng thấp rất dễ chia lưới tự động cho hình học phức tạp. Với cùng mức độ chính xác, phương pháp này đòi hỏi số bậc tự do ít hơn khoảng 25% đến 35% so với việc dùng phần tử bậc cao trên phần mềm COMSOL, giúp tiết kiệm bộ nhớ máy tính.

Luận văn đã kiểm chứng thuật toán qua các mô hình thực nghiệm số nào? Luận văn đã kiểm chứng thuật toán qua 3 bài toán chuẩn: vách phòng kín kích thước 3m x 8m tương tác với không khí chịu áp suất âm thanh, cấu kiện tấm ngàm 4 cạnh đặt trên mặt nước với 3328 bậc tự do, và tấm tròn đáy bể chứa khí với 3685 bậc tự do. Tất cả đều cho kết quả trùng khớp cao với phần mềm COMSOL.

Kết luận

Năm đóng góp cốt lõi của luận văn thạc sĩ bao gồm:

  • Xây dựng thành công phương pháp ghép nối CS-MIN3/FS-FEM cho bài toán phân tích dao động tự do tấm Mindlin tương tác với lưu chất trong không gian 3 chiều.
  • Khắc phục triệt để hiện tượng khóa cắt và quá cứng của phần tử hữu hạn tiêu chuẩn, đạt độ chính xác tần số riêng trên 98,5%.
  • Kiểm chứng hoàn hảo thuật toán qua 3 mô hình số điển hình với quy mô từ 1848 đến 3685 bậc tự do, mô tả chính xác 9 mode dao động đầu tiên.
  • Rút ngắn thời gian hội tụ và giảm hơn 30% tài nguyên tính toán so với các phương pháp phần tử hữu hạn truyền thống.
  • Đặt nền móng vững chắc cho việc ứng dụng công nghệ mô phỏng số FSI 3D vào kỹ thuật xây dựng công trình dân dụng, công nghiệp và hàng hải.

Đóng góp chính của công trình là cung cấp một công cụ số tin cậy, chính xác và hiệu quả cao trong việc mô phỏng tương tác kết cấu rắn và môi trường âm học phức tạp. Nhóm nghiên cứu dự kiến mở rộng mô hình cho bài toán phi tuyến hình học và vật liệu trong giai đoạn 2026-2027. Độc giả và các kỹ sư hãy tham khảo toàn văn luận văn để ứng dụng giải pháp mô phỏng tiên tiến này vào các dự án thực tế.