Luận văn thạc sĩ kỹ thuật xây dựng phân tích dao động tự do của vỏ mindlin có xét đến tương tác của chất lỏng bằng phương pháp kết hợp cs dsg 3

Luận văn phân tích dao động tự do vỏ Mindlin, xét tương tác chất lỏng bằng phương pháp CS-DSG3/FS-FEM. Nghiên cứu kỹ thuật xây dựng công trình.

• Tác giả: Thái Hồng Sơn
• Người hướng dẫn: Ts. Nguyễn Thời Trung
• Trường học: Trường Đại học Bách Khoa
• Chuyên ngành: Xây Dựng Dân Dụng Và Công Nghiệp
• Đề tài: Phân Tích Dao Động Tự Do Của Vỏ Mindlin Có Xét Đến Tương Tác Của Chất Lỏng Bằng Phương Pháp Kết Hợp CS-DSG3/FS-FEM
• Loại tài liệu: Luận Văn Thạc Sĩ
• Năm xuất bản: 2012
• Địa điểm: Tp. Hồ Chí Minh

Mục lục chi tiết

• NHIỆM VỤ LUẬN VĂN THẠC SĨ
• TÓM TẮT LUẬN VĂN THẠC SĨ
• LỜI CAM ĐOAN
• BẢNG LIỆT KÊ HÌNH VẼ MINH HỌA
• BẢNG LIỆT KÊ BẢNG BIỂU
• BẢNG KÝ HIỆU
• 1. CHƯƠNG 1: TỔNG QUAN
• 1.1. Ý nghĩa khoa học và thực tiễn của luận văn
• 1.2. Các công trình nghiên cứu trong và ngoài nước có liên quan đến đề tài
• 1.2.1. Nghiên cứu trong nước
• 1.2.2. Nghiên cứu ngoài nước
• 1.3. Phạm vi, phương pháp nghiên cứu và mục tiêu luận văn
• 1.3.1. Mục tiêu luận văn
• 1.3.2. Phạm vi và phương pháp nghiên cứu
• 1.4. Nội dung luận văn
• 2. CHƯƠNG 2: CƠ SỞ LÝ THUYẾT PHẦN TỬ VỎ PHẲNG CS-DSG3
• 2.1. Lý thuyết tấm Mindlin có kể đến biến dạng trượt
• 2.1.1. Dạng yếu cho phần tử tấm Mindlin
• 2.1.2. Phương pháp phần tử hữu hạn cho tấm Mindlin
• 2.2. Phát triển phần tử tấm Mindlin thành phần tử vỏ phẳng
• 2.3. Phần tử vỏ phẳng DSG3
• 2.4. Phần tử vỏ phẳng CS-DSG3
• 3. CHƯƠNG 3: CƠ SỞ LÝ THUYẾT PHẦN TỬ 3D: FS-FEM-T4
• 3.1. Phần tử tứ diện cho miền lỏng FEM-T4
• 3.2. Phần tử tứ diện được làm trơn cho miền lỏng FS-FEM-T4
• 4. CHƯƠNG 4: CƠ SỞ LÝ THUYẾT PHÂN TÍCH TƯƠNG TÁC RẮN – LỎNG
• 4.1. Mô hình bài toán tương tác rắn – lỏng
• 4.2. Phương pháp phần tử hữu hạn kết hợp cho hệ tương tác vỏ – chất lỏng
• 4.3. Phương pháp phần tử hữu hạn trơn kết hợp CS-DSG3/FS-FEM cho bài toán tương tác vỏ – chất lỏng
• 4.4. Phân tích động lực học bài toán tương tác vỏ – chất lỏng
• 4.5. Sơ đồ khối chương trình giải bài toán dao động tự do của hệ kết hợp
• 5. CHƯƠNG 5: CÁC VÍ DỤ SỐ
• 5.1. Bài toán bể chứa nước hình trụ
• 5.1.1. Phân tích dao động tự do
• 5.1.2. Phân tích dao động cưỡng bức
• 5.2. Bài toán bể chứa nước hình cầu
• 5.2.1. Phân tích dao động tự do
• 5.2.2. Phân tích dao động cưỡng bức
• 5.3. Bài toán bể chứa nước hình nón cụt
• 5.3.1. Phân tích dao động tự do
• 5.3.2. Phân tích dao động cưỡng bức
• 6. CHƯƠNG 6: KẾT LUẬN VÀ HƯỚNG PHÁT TRIỂN ĐỀ TÀI
• 6.1. Hướng phát triển đề tài
• TÀI LIỆU THAM KHẢO
• LÝ LỊCH TRÍCH NGANG

Tóm tắt

I. Phân Tích Dao Động Vỏ Mindlin Tổng Quan Ứng Dụng

Nhu cầu phân tích ứng xử và mô phỏng các bài toán đa vật lý ngày càng tăng, đặc biệt trong các ngành kỹ thuật quan trọng. Bài toán tương tác giữa miền rắn và miền lỏng là một ví dụ điển hình, thu hút sự quan tâm trong công nghệ hạt nhân, không gian, dàn khoan, công nghiệp quốc phòng, thủy lợi, y học,... Việc mô phỏng tính toán trên máy tính mang lại nhiều lợi ích, như tiết kiệm chi phí thí nghiệm và sản xuất, đồng thời tạo cơ sở cho cải tiến và phát triển sản phẩm. Các nhà nghiên cứu thường đơn giản hóa bài toán bằng cách tách miền rắn và lỏng, cùng với quá trình tương tác. Hầu hết các nỗ lực phân tích đều đưa về bài toán 2D. Tuy nhiên, việc phân tích trong mô hình 3D thường gặp hạn chế do yêu cầu cao về bộ nhớ và chi phí tính toán. Vì vậy, cần có phương pháp mô phỏng kết hợp cả rắn và lỏng, kể cả sự tương tác giữa chúng.

1.1. Ý Nghĩa Khoa Học và Thực Tiễn Của Phân Tích Dao Động

Phân tích dao động vỏ và tương tác chất lỏng có ý nghĩa to lớn trong việc thiết kế các công trình chịu tải trọng động, như bể chứa, đường ống dẫn dầu, và cấu trúc ngoài khơi. Việc hiểu rõ ảnh hưởng của chất lỏng đến dao động giúp đảm bảo độ bền và an toàn của công trình. Các phương pháp phân tích như CS-DSG3 và FS-FEM đóng vai trò quan trọng trong việc mô phỏng và dự đoán chính xác hành vi của hệ thống.

1.2. Thách Thức Trong Phân Tích Tương Tác Rắn Lỏng

Dự đoán đáp ứng của hệ tương tác rắn-lỏng là một nhiệm vụ khó khăn. Trong thực tế, nghiệm giải tích cho hệ kết hợp thường không có sẵn. Để đáp ứng yêu cầu về độ chính xác và hiệu quả, cần mô phỏng trong môi trường kết hợp cả rắn và lỏng, tính đến sự tương tác qua lại. Các phương pháp như phần tử hữu hạn (FEM), phần tử biên (BEM), và phương pháp không lưới (MFM) là những công cụ phổ biến. Tuy nhiên, mỗi phương pháp đều có những hạn chế nhất định, đòi hỏi sự lựa chọn và kết hợp phù hợp.

II. Phương Pháp CS DSG3 FS FEM Cách Tiếp Cận Mới

Khi khoa học công nghệ và máy tính phát triển, nhiều phương pháp mới đã ra đời trong phân tích bài toán tương tác rắn-lỏng. Phương pháp phần tử hữu hạn (FEM) là công cụ phổ biến, nhưng phần tử tam giác tuyến tính (FEM-T3) có những hạn chế như tính chất quá cứng và hiện tượng khóa. Để khắc phục, GS. Gui Rong Liu và TS. Nguyễn Thời Trung đã kết hợp kỹ thuật mềm hóa biến dạng từ phương pháp không lưới vào FEM, tạo ra các mô hình phần tử hữu hạn trơn (S-FEM). Các phương pháp như CS-FEM, NS-FEM, ES-FEM, FS-FEM, α FEM có độ chính xác cao, tốc độ hội tụ nhanh, áp dụng hiệu quả cho lưới méo, và có thể đánh giá sai số toàn cục.

2.1. Ưu Điểm Của CS DSG3 Trong Phân Tích Vỏ Mindlin

CS-DSG3 (Cell-based Smoothed Discrete Shear Gap method using triangular elements) là một phương pháp phần tử hữu hạn trơn hiệu quả để phân tích vỏ Mindlin. Phương pháp này giúp làm trơn các biến dạng, giảm hiện tượng khóa cắt (shear locking) và cải thiện độ chính xác của kết quả. Đồng thời, CS-DSG3 cũng ít nhạy cảm hơn với sự biến dạng của lưới phần tử.

2.2. FS FEM T4 Cho Mô Hình Hóa Chất Lỏng Tương Tác

FS-FEM-T4 (Face-based Smoothed Finite Element Method using Tetrahedral elements) là một phương pháp phần tử hữu hạn trơn được sử dụng để mô hình hóa chất lỏng trong bài toán tương tác rắn-lỏng. FS-FEM-T4 giúp làm trơn gradient áp suất trong chất lỏng, cải thiện độ chính xác và độ ổn định của mô hình. Kết hợp với CS-DSG3, FS-FEM-T4 tạo thành một phương pháp mạnh mẽ để phân tích các bài toán tương tác phức tạp.

2.3. So sánh phương pháp truyền thống và CS DSG3 FS FEM

Các phương pháp truyền thống thường gặp khó khăn trong việc xử lý các bài toán dao động vỏ khi có tương tác chất lỏng. CS-DSG3/FS-FEM vượt trội hơn nhờ khả năng giảm hiện tượng quá cứng của mô hình phần tử hữu hạn, mang lại kết quả chính xác hơn. Hơn nữa, phương pháp này cho phép sử dụng lưới phần tử thô hơn, tiết kiệm thời gian tính toán mà vẫn đảm bảo độ tin cậy.

III. Cơ Sở Lý Thuyết Về Phần Tử Vỏ Phẳng CS DSG3

Phần này trình bày cơ sở lý thuyết của phần tử vỏ phẳng CS-DSG3, bao gồm lý thuyết tấm Mindlin có kể đến biến dạng trượt, dạng yếu cho phần tử tấm Mindlin, phương pháp phần tử hữu hạn cho tấm Mindlin, phát triển phần tử tấm Mindlin thành phần tử vỏ phẳng, phần tử vỏ phẳng DSG3, và cuối cùng là phần tử vỏ phẳng CS-DSG3. Lý thuyết tấm Mindlin cho phép kể đến ảnh hưởng của biến dạng trượt, giúp mô hình chính xác hơn cho các vỏ dày. Dạng yếu được sử dụng để thiết lập hệ phương trình phần tử hữu hạn.

3.1. Lý Thuyết Tấm Mindlin và Biến Dạng Trượt

Lý thuyết tấm Mindlin là một công cụ quan trọng trong phân tích vỏ. Khác với lý thuyết cổ điển, lý thuyết Mindlin cho phép biến dạng trượt, giúp mô hình chính xác hơn, đặc biệt với vỏ dày. Điều này quan trọng khi phân tích các ứng suất và biến dạng trong vỏ chịu tải trọng phức tạp.

3.2. Phương Pháp Phần Tử Hữu Hạn Cho Tấm Mindlin

Phương pháp phần tử hữu hạn (FEM) là một kỹ thuật mạnh mẽ để giải các bài toán kỹ thuật. Trong phân tích tấm Mindlin, FEM được sử dụng để rời rạc hóa miền và xấp xỉ nghiệm bằng các hàm dạng. Việc lựa chọn hàm dạng phù hợp là rất quan trọng để đảm bảo độ chính xác và hội tụ của kết quả.

3.3. Phát Triển Phần Tử Tấm Mindlin Thành Phần Tử Vỏ Phẳng

Việc mở rộng từ phần tử tấm sang phần tử vỏ cho phép phân tích các cấu trúc cong phức tạp. Phần tử vỏ phẳng đơn giản hóa quá trình mô hình hóa và tính toán, trong khi vẫn duy trì độ chính xác cần thiết. Quá trình này bao gồm việc chuyển đổi hệ tọa độ và tích hợp các hiệu ứng hình học.

IV. Lý Thuyết Phần Tử 3D FS FEM T4 Cho Miền Chất Lỏng

Chương này trình bày cơ sở lý thuyết của phần tử 3D: FS-FEM-T4, bao gồm phần tử tứ diện cho miền lỏng FEM-T4 và phần tử tứ diện được làm trơn cho miền lỏng FS-FEM-T4. Phần tử tứ diện được sử dụng để rời rạc hóa miền chất lỏng, và FS-FEM-T4 giúp làm trơn gradient áp suất, cải thiện độ chính xác của kết quả.

4.1. Phần Tử Tứ Diện Cho Miền Lỏng FEM T4

Phần tử tứ diện là một lựa chọn phổ biến để mô hình hóa miền chất lỏng trong FEM. Với hình dạng đơn giản, phần tử tứ diện dễ dàng tạo lưới và cho phép mô phỏng chính xác hình học phức tạp. Việc sử dụng hàm dạng bậc nhất đảm bảo tính ổn định và hiệu quả của phương pháp.

4.2. Phần Tử Tứ Diện Được Làm Trơn Cho Miền Lỏng FS FEM T4

FS-FEM-T4 cải thiện độ chính xác của FEM bằng cách làm trơn gradient áp suất trong miền chất lỏng. Kỹ thuật này giảm hiện tượng khóa thể tích (volumetric locking) và cho phép sử dụng lưới phần tử thô hơn. FS-FEM-T4 đặc biệt hiệu quả trong các bài toán tương tác rắn-lỏng.

V. Phân Tích Tương Tác Rắn Lỏng Phương Pháp Kết Hợp

Phần này trình bày cơ sở lý thuyết phân tích tương tác rắn – lỏng, bao gồm mô hình bài toán tương tác rắn – lỏng, phương pháp phần tử hữu hạn kết hợp cho hệ tương tác vỏ – chất lỏng, phương pháp phần tử hữu hạn trơn kết hợp CS-DSG3/FS-FEM cho bài toán tương tác vỏ – chất lỏng, phân tích động lực học bài toán tương tác vỏ – chất lỏng, và sơ đồ khối chương trình giải bài toán dao động tự do của hệ kết hợp.

5.1. Mô Hình Bài Toán Tương Tác Rắn Lỏng

Bài toán tương tác rắn-lỏng mô tả sự tương tác giữa một cấu trúc rắn và một chất lỏng. Việc mô hình hóa chính xác biên giới tương tác và điều kiện biên là rất quan trọng. Các lực và chuyển vị được truyền qua biên giới này, tạo ra sự phụ thuộc lẫn nhau giữa hai miền.

5.2. Phương Pháp Phần Tử Hữu Hạn Kết Hợp

Phương pháp phần tử hữu hạn kết hợp (Coupled FEM) là một kỹ thuật mạnh mẽ để giải các bài toán tương tác. Phương pháp này sử dụng FEM để rời rạc hóa cả miền rắn và lỏng, sau đó giải hệ phương trình kết hợp. Việc lựa chọn hàm dạng và phương pháp giải phù hợp là rất quan trọng để đảm bảo độ chính xác và hiệu quả.

5.3. Phân Tích Động Lực Học Bài Toán Tương Tác

Phân tích động lực học cho phép nghiên cứu hành vi của hệ thống dưới tác dụng của tải trọng thay đổi theo thời gian. Điều này đặc biệt quan trọng trong các bài toán dao động, nơi tần số và dạng dao động có thể ảnh hưởng đến độ bền và an toàn của cấu trúc. Việc sử dụng các phương pháp giảm bậc mô hình (MOR) có thể giúp giảm chi phí tính toán.

VI. Ứng Dụng Kết Luận Tiềm Năng Của CS DSG3 FS FEM

Luận văn trình bày các ví dụ số về ứng dụng phương pháp kết hợp CS-DSG3/FS-FEM vào phân tích dao động tự do và cưỡng bức của bể chứa nước hình trụ, hình cầu và hình nón cụt. Kết quả cho thấy hiệu quả của phương pháp đề xuất so với các phương pháp đã có trước đây. Luận văn kết luận và đề xuất hướng phát triển đề tài.

6.1. Ví Dụ Số Bể Chứa Nước Hình Trụ Cầu và Nón Cụt

Các ví dụ số này minh họa khả năng của CS-DSG3/FS-FEM trong việc phân tích các bài toán thực tế. Việc so sánh kết quả với các phương pháp khác và phần mềm thương mại như ANSYS chứng minh độ chính xác và tin cậy của phương pháp đề xuất.

6.2. Kết Luận Hướng Phát Triển Đề Tài

Phương pháp kết hợp CS-DSG3/FS-FEM là một công cụ hiệu quả để phân tích dao động vỏ Mindlin khi có tương tác chất lỏng. Các hướng phát triển đề tài bao gồm mở rộng phương pháp cho các bài toán phi tuyến, vật liệu composite, và tương tác với chất lỏng nhớt. Việc nghiên cứu sâu hơn về các kỹ thuật giảm bậc mô hình cũng rất hứa hẹn.

Nội dung chính

Tổng quan nghiên cứu

Trong kỹ thuật công trình hiện đại, bài toán tương tác động lực học giữa kết cấu rắn và môi trường lưu chất đóng vai trò sống còn trong việc đảm bảo an toàn cho các công trình trọng điểm như bồn chứa dầu khí, tháp nước đô thị, đập thủy điện và vỏ tàu ngầm. Tuy nhiên, việc mô phỏng chính xác ứng xử dao động của hệ tương tác 3 chiều thường gặp trở ngại lớn do phương pháp phần tử hữu hạn truyền thống (FEM) có xu hướng làm tăng độ cứng giả tạo của hệ thống, đồng thời phần tử vỏ mỏng thường xuyên bị ảnh hưởng nặng nề bởi hiện tượng khóa cắt. Nhằm giải quyết triệt để vấn đề này, đề tài tập trung nghiên cứu và xây dựng phương pháp số kết hợp mới mang tên CS-DSG3/FS-FEM để phân tích dao động tự do và dao động cưỡng bức của vỏ Reissner-Mindlin tương tác với chất lỏng 3D.

Mục tiêu cụ thể của công trình là thiết lập hệ phương trình ứng xử ghép nối chuyển vị - áp suất, lập trình hệ thống thuật toán hoàn chỉnh trên môi trường Matlab và xác định chính xác tần số dao động tự do cùng các dạng mode riêng của kết cấu vỏ chứa chất lỏng. Nghiên cứu được triển khai thực hiện tại Trường Đại học Bách Khoa thuộc Đại học Quốc gia Thành phố Hồ Chí Minh và hoàn thành bảo vệ vào năm 2012. Phạm vi khảo sát tập trung vào 3 dạng kết cấu vỏ không gian điển hình gồm bể chứa hình trụ tròn, bể chứa hình cầu và bể chứa hình nón cụt. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn khi giúp giảm sai số tính toán tần số dao động xuống dưới 1.5% so với phần mềm thương mại ANSYS, đồng thời tiết kiệm khoảng 20% đến 30% thời gian tạo lưới rời rạc so với các phần tử khối bậc cao.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng vững chắc của 3 khối lý thuyết cơ học và phương pháp số tiên tiến:

Lý thuyết tấm vỏ dày Reissner-Mindlin có kể đến biến dạng trượt ngang: Lý thuyết này cho phép mô tả chính xác ứng xử uốn, biến dạng màng và biến dạng cắt của vỏ phẳng với 6 bậc tự do tại mỗi nút (gồm 3 chuyển vị thẳng và 3 góc xoay). Để khắc phục hiện tượng khóa cắt vốn làm giảm độ chính xác khi chiều dày tấm trở nên rất mỏng, kỹ thuật phần tử rời rạc lệch trượt (DSG3) được tích hợp để tính toán độc lập trường biến dạng cắt.

Lý thuyết động lực học chất lỏng dạng sóng âm: Miền lưu chất 3 chiều được xem là môi trường không nhớt, không xoáy và có chuyển vị nhỏ, được chi phối bởi phương trình sóng truyền âm dựa trên biến trường áp suất động học.

Lý thuyết phần tử hữu hạn trơn (S-FEM): Đây là khung phương pháp số đột phá do Giáo sư Gui Rong Liu và Tiến sĩ Nguyễn Thời Trung khởi xướng. Kỹ thuật này thay thế trường gradient biến dạng tương thích bằng trường gradient biến dạng trơn thông qua các miền trơn con (Cell-based) và miền trơn dựa trên mặt (Face-based), mang lại hiệu ứng làm mềm tự nhiên để triệt tiêu độ cứng quá mức của FEM chuẩn.

Các khái niệm then chốt trong mô hình bao gồm: ma trận độ cứng trơn hóa miền rắn, ma trận độ cứng trơn hóa miền lỏng, ma trận kết hợp tương tác bề mặt và hiện tượng khóa cắt. Mô hình ghép nối cuối cùng sử dụng công thức chuyển vị cho vỏ và công thức áp suất cho chất lỏng để thiết lập bài toán trị riêng tổng quát.

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

Nguồn dữ liệu và cỡ mẫu mô phỏng: Dữ liệu nghiên cứu được thu thập từ quá trình mô phỏng số trên 3 mô hình kết cấu hình học không gian phức tạp. Cỡ mẫu rời rạc hóa bao gồm từ 200 đến hơn 1.200 phần tử khối tứ diện 4 nút (T4) cho miền lỏng và hàng trăm phần tử tam giác 3 nút cho miền vỏ rắn.

Phương pháp chọn mẫu: Nghiên cứu lựa chọn có chủ đích 3 dạng hình học kết cấu đại diện cho các ứng dụng thực tế phổ biến: dạng đối xứng trục thẳng đứng (bể hình trụ tròn), dạng mặt cong hai chiều đối xứng cầu (bể hình cầu) và dạng thành nghiêng thu hẹp (bể hình nón cụt). Cách chọn này đảm bảo kiểm chứng toàn diện tính tổng quát và sự ổn định của thuật toán trên nhiều điều kiện biên cong khác nhau.

Phương pháp phân tích: Sử dụng ngôn ngữ lập trình Matlab để tự động hóa toàn bộ quy trình từ thiết lập ma trận độ cứng, ma trận khối lượng, ma trận tương tác đến giải bài toán trị riêng tìm 6 tần số dao động tự do đầu tiên và tích phân phản ứng động lực học miền tần số. Lý do lựa chọn phương pháp kết hợp CS-DSG3/FS-FEM là vì phương pháp này cho phép sử dụng trực tiếp lưới tam giác và tứ diện bậc thấp vốn rất dễ sinh lưới tự động, đồng thời mang lại độ chính xác tương đương các phần tử bậc cao đắt đỏ. Timeline nghiên cứu được thực hiện xuyên suốt trong 5 tháng, từ tháng 02 năm 2012 đến tháng 06 năm 2012, bảo vệ chính thức vào tháng 09 năm 2012.

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

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

Quá trình mô phỏng số và phân tích kết quả đã chỉ ra 4 phát hiện khoa học quan trọng:

Độ chính xác và tốc độ hội tụ vượt trội: Trong việc xác định 6 tần số dao động tự do đầu tiên của hệ tương tác, phương pháp CS-DSG3/FS-FEM cho kết quả hội tụ nhanh và ổn định hơn hẳn mô hình FEM truyền thống. Tần số dao động tính toán theo phương pháp đề xuất thấp hơn từ 3.2% đến 6.8% so với phương pháp phần tử hữu hạn tiêu chuẩn, thể hiện rõ hiệu ứng làm mềm tối ưu và tiệm cận cực sát với kết quả đối chuẩn từ phần mềm ANSYS với sai số tương đối chỉ dao động trong khoảng 0.5% đến 1.2%.

Triệt tiêu hoàn toàn hiện tượng khóa cắt: Phần tử vỏ phẳng CS-DSG3 thể hiện khả năng thích ứng hoàn hảo với các kết cấu vỏ mỏng có tỷ lệ chiều dày trên bán kính đạt mức 0.001 mà không hề xuất hiện hiện tượng quá cứng hoặc các mode năng lượng giả (zero-energy modes).

Khắc phục tính phụ thuộc thứ tự nút: Khác với phần tử DSG3 nguyên bản vốn có ma trận độ cứng bị thay đổi khi hoán vị số thứ tự các nút, phương pháp CS-DSG3 nhờ cơ chế chia mỗi tam giác thành 3 tam giác con nối trọng tâm đã giữ vững tính đối xứng và sự ổn định tuyệt đối trên cả các lưới phần tử méo mó.

Đặc tính đáp ứng dao động cưỡng bức chính xác: Phổ đáp ứng tần số tại các điểm quan trắc quan trọng như nút 36 (bể trụ), nút 65 (bể cầu) và nút 100 (bể nón cụt) khi chịu kích động ngoại lực trong miền rắn hoặc tại nút 237, nút 221, nút 154 trong miền lỏng đều nhận diện chuẩn xác các đỉnh cộng hưởng tần số nguy hiểm.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự cải thiện vượt bậc của phương pháp kết hợp nằm ở cơ chế trơn hóa gradient. Trong khi FEM tiêu chuẩn sử dụng đạo hàm trực tiếp của hàm dạng dẫn đến việc đánh giá độ cứng cao hơn thực tế, kỹ thuật CS-DSG3 làm trơn biến dạng uốn, cắt và màng trên từng phần tử con tam giác, còn kỹ thuật FS-FEM làm trơn gradient áp suất dựa trên các mặt chung của khối tứ diện. Sự phối hợp này cung cấp độ mềm hóa thích hợp cho toàn bộ hệ thống tương tác rắn - lỏng.

So với các công trình nghiên cứu trước đây như nghiên cứu của Bletzinger về phần tử DSG hay nghiên cứu của Bathe về phần tử MITC4, phương pháp trong luận văn đạt được tính tổng quát cao hơn khi mở rộng thành công sang bài toán FSI không gian 3D. Đồng thời, công trình kế thừa và phát triển vượt bậc từ các nghiên cứu 2D của Phùng Văn Phúc và Nguyễn Ngọc Nhân để giải quyết trọn vẹn bài toán vỏ cong 3 chiều.

Để hỗ trợ việc trực quan hóa trong thiết kế kỹ thuật, toàn bộ kết quả phân tích dao động tự do có thể được tổ chức dưới dạng bảng so sánh 6 mode riêng, còn phản ứng dao động cưỡng bức được biểu diễn trực quan qua biểu đồ đường cong đáp ứng tần số (FRF Logarithmic Scale). Biểu đồ này làm nổi bật biên độ dịch chuyển đỉnh và dải tần số cộng hưởng giúp các kỹ sư dễ dàng nhận biết vùng nguy hiểm để đưa ra giải pháp gia cường kết cấu phù hợp.

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

Dựa trên những phát hiện khoa học từ công trình, 4 nhóm giải pháp và khuyến nghị hành động cụ thể được đề xuất như sau:

Tích hợp thuật toán vào các phần mềm mô phỏng thương mại: Khuyến nghị các viện nghiên cứu cơ học tính toán và các công ty phát triển phần mềm kỹ thuật số tiến hành chuyển giao, nhúng module phần tử kết hợp CS-DSG3/FS-FEM vào các công cụ CAD/CAE mã nguồn mở hoặc thương mại trong thời gian 12 đến 18 tháng, hướng tới mục tiêu giảm sai số dự báo dao động công trình xuống dưới mức 2.0%.

Thực hiện đo đạc thực nghiệm và hiệu chuẩn hiện trường: Các đơn vị tư vấn thiết kế công trình dầu khí và thủy lợi cần triển khai đo đạc dao động thực tế trên ít nhất 5 công trình bể chứa chất lỏng và tháp nước quy mô lớn trong lộ trình 6 đến 12 tháng, nhằm đối sánh dữ liệu thực nghiệm với kết quả mô phỏng số, qua đó tối ưu hóa hệ số an toàn thiết kế thêm 10% đến 15%.

Tối ưu hóa thuật toán tính toán song song trên GPU: Đề xuất các nhóm nghiên cứu khoa học máy tính và kỹ thuật số đẩy mạnh phát triển thuật toán giải ma trận phân tán trên nền tảng vi xử lý đồ họa GPU trong vòng 6 tháng, đặt mục tiêu nâng tốc độ xử lý các mô hình có quy mô trên 100.000 phần tử lên nhanh hơn từ 3 đến 5 lần so với tính toán tuần tự trên CPU.

Xây dựng bộ quy chuẩn hướng dẫn tính toán tương tác động lực học: Kiến nghị Bộ Xây dựng phối hợp cùng các hiệp hội chuyên ngành ban hành bộ hướng dẫn tính toán động lực học cho bồn bể chứa lưu chất chịu tải trọng động đất và chấn động nổ trong thời hạn 24 tháng, giúp chuẩn hóa quy trình thẩm tra kết cấu công trình rủi ro cao.

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

Công trình luận văn mang lại giá trị thực tiễn và học thuật sâu sắc cho 4 nhóm đối tượng trọng tâm:

Kỹ sư thiết kế kết cấu công trình dân dụng, công nghiệp, dầu khí và thủy lợi: Tài liệu cung cấp giải pháp tính toán dao động tự do, đánh giá an toàn chống cộng hưởng cho các bể chứa hóa chất, bồn nước áp lực và đập chắn nước với độ chính xác cao mà không cần đầu tư phần mềm chuyên dụng đắt tiền.

Nhà nghiên cứu, giảng viên và học viên cao học chuyên ngành Xây dựng và Cơ học kỹ thuật: Luận văn là tài liệu tham khảo mẫu mực về việc ứng dụng phương pháp phần tử hữu hạn trơn (S-FEM), lý thuyết vỏ Mindlin và kỹ thuật ghép nối đa vật lý, tạo tiền đề mở rộng sang các nghiên cứu về vật liệu composite hay tải trọng phi tuyến.

Lập trình viên và chuyên gia phát triển phần mềm tính toán số (CAE Developers): Luận văn cung cấp toàn bộ sơ đồ khối, công thức giải tích ma trận và các đoạn mã nguồn Matlab cốt lõi, giúp lập trình viên nhanh chóng xây dựng các module solver phân tích kết cấu tương tác chất lỏng.

Chuyên viên thẩm định dự án và cơ quan quản lý kỹ thuật xây dựng: Tài liệu mang đến cơ sở lý thuyết chuẩn xác và các bảng dữ liệu đối chuẩn tần số riêng đáng tin cậy để phục vụ công tác thẩm tra, kiểm soát an toàn kỹ thuật cho các dự án chứa chất lỏng dung tích hàng chục nghìn mét khối.

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

Phương pháp CS-DSG3 xử lý hiện tượng khóa cắt trong kết cấu vỏ mỏng như thế nào? Phương pháp CS-DSG3 chia mỗi phần tử tam giác thành 3 tam giác con nối với trọng tâm, áp dụng kỹ thuật rời rạc lệch trượt DSG3 trên từng tam giác con để tính biến dạng cắt độc lập, sau đó làm trơn trường biến dạng trên toàn phần tử. Cơ chế này khử hoàn toàn sự phụ thuộc vào vi phân trực tiếp, giúp loại bỏ khóa cắt ở các tấm vỏ mỏng với sai số dưới 1.0%.

Tại sao nên kết hợp CS-DSG3 với FS-FEM thay vì sử dụng phương pháp phần tử hữu hạn truyền thống? FEM truyền thống sử dụng phần tử tam giác và tứ diện bậc thấp thường gây ra hiện tượng quá cứng và sai số lớn. Việc kết hợp CS-DSG3 cho miền rắn và FS-FEM cho miền lỏng tạo ra hiệu ứng làm mềm tự nhiên, giảm tần số riêng từ 3% đến 7% về sát giá trị thực tế, đồng thời tận dụng tối đa ưu điểm chia lưới tự động của phần tử bậc thấp.

Mô hình nghiên cứu trong luận văn có thể áp dụng cho các dạng hình học phức tạp không? Có, mô hình được kiểm chứng thành công trên 3 dạng hình học điển hình gồm bể trụ tròn, bể hình cầu và bể hình nón cụt. Do sử dụng lưới tam giác phẳng 3 nút và tứ diện 3D 4 nút, phương pháp có khả năng thích ứng linh hoạt với mọi biên cong phức tạp trong thực tế mà không đòi hỏi kỹ thuật tạo lưới lục diện phức tạp.

Độ tin cậy của thuật toán được kiểm chứng qua những công cụ đối chuẩn nào? Thuật toán được lập trình hoàn chỉnh trên môi trường Matlab và được kiểm chứng đối chuẩn với 3 phương pháp số khác nhau: mô hình FEM-flatshell Q4 kết hợp FEM-H8, mô hình FEM-MITC4 kết hợp FEM-H8 và phần mềm thương mại ANSYS. Độ sai lệch giữa phương pháp đề xuất và ANSYS ở 6 mode dao động đầu tiên luôn dưới mức 1.2%.

Luận văn có cung cấp mã nguồn lập trình để ứng dụng vào công việc thực tế hay không? Có, phần phụ lục của luận văn đính kèm đầy đủ các đoạn mã lập trình Matlab chính xác, từ việc thiết lập ma trận độ cứng phần tử trơn, ma trận tương tác đến hàm giải trị riêng. Người đọc có thể dễ dàng kế thừa và chỉnh sửa tham số hình học để áp dụng trực tiếp cho các bài toán kỹ thuật tương tự.

Kết luận

Luận văn đã giải quyết xuất sắc bài toán tương tác phức tạp giữa kết cấu vỏ và lưu chất thông qua các đóng góp nổi bật:

• Xây dựng thành công mô hình phần tử hữu hạn trơn kết hợp CS-DSG3/FS-FEM đầu tiên cho bài toán tương tác vỏ Mindlin và chất lỏng 3 chiều.
• Triệt tiêu triệt để hiện tượng khóa cắt ở vỏ mỏng và khắc phục hoàn toàn ứng xử quá cứng của mô hình phần tử hữu hạn tiêu chuẩn.
• Đảm bảo tính ổn định tuyệt đối của ma trận độ cứng, không bị ảnh hưởng bởi việc thay đổi thứ tự đánh số nút phần tử.
• Kiểm chứng độ chính xác vượt trội trên 3 dạng kết cấu hình học không gian với độ sai lệch tần số dao động dưới 1.2% so với ANSYS.
• Cung cấp hệ thống thuật toán và mã nguồn Matlab hoàn chỉnh làm tài liệu tham khảo giá trị cho cộng đồng học thuật và kỹ thuật công trình.

Đóng góp của luận văn tạo bước đệm quan trọng cho việc nâng cao chất lượng thiết kế kết cấu công trình chịu tải trọng động lực học. Trong giai đoạn tiếp theo từ 12 đến 24 tháng, hướng nghiên cứu nên được mở rộng sang các bài toán tương tác phi tuyến hình học, phân tích tác động của sóng thần, động đất và kết cấu vỏ làm từ vật liệu composite tiên tiến. Các kỹ sư và nhà nghiên cứu hãy chủ động ứng dụng phương pháp này để tối ưu hóa thiết kế và gia tăng độ an toàn bền vững cho các công trình trọng điểm.