Tổng quan nghiên cứu

Bài toán tương tác giữa kết cấu vỏ và môi trường lưu chất là một trong những thách thức cơ học đa vật lý trọng tâm trong kỹ thuật xây dựng công trình dân dụng, công nghiệp và quốc phòng. Thực tế cho thấy, khoảng 70% các sự cố hư hỏng bồn chứa chất lỏng công nghiệp hoặc tháp nước chịu rung chấn bắt nguồn từ việc đánh giá không đầy đủ phản ứng động lực học của hệ kết hợp. Khi chất lỏng chuyển động, áp suất động tác động trực tiếp lên thành vỏ và làm thay đổi đáng kể các đặc trưng tần số tự nhiên. Tuy nhiên, việc áp dụng phương pháp phần tử hữu hạn truyền thống thường dẫn đến hiện tượng quá cứng của mô hình và sai số tần số dao động có thể vượt mức 20% ở các kết cấu mỏng.

Luận văn thạc sĩ chuyên ngành Xây dựng Dân dụng và Công nghiệp được 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, hoàn thành vào tháng 10 năm 2012. Mục tiêu trọng tâm của nghiên cứu là xây dựng, cài đặt thuật toán và kiểm chứng phương pháp số kết hợp mới mang tên CS-DSG3 và FS-FEM nhằm phân tích chính xác dao động tự do và dao động cưỡng bức của vỏ Mindlin có xét tương tác chất lỏng 3 chiều. Phạm vi nghiên cứu khảo sát chi tiết trên 3 dạng hình học bồn chứa phổ biến gồm hình trụ tròn, hình cầu và hình nón cụt. Giá trị thực tiễn của công trình thể hiện ở việc nâng cao độ chính xác dự báo tần số dao động lên trên 98,5% so với phần mềm thương mại ANSYS, đồng thời giảm khoảng 40% chi phí tính toán và bộ nhớ lưu trữ cho hệ phương trình tương tác.

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 kết hợp giữa cơ học kết cấu vỏ tấm và động lực học chất lỏng sóng âm tuyến tính. Về miền rắn, lý thuyết tấm vỏ Reissner-Mindlin được áp dụng có kể đến biến dạng trượt ngang với hệ số hiệu chỉnh cắt k bằng 5/6, cho phép mô tả chính xác cả vỏ mỏng và vỏ dày trung bình. Phần tử vỏ phẳng tam giác 3 nút được phát triển dựa trên sự kết hợp giữa biến dạng màng và biến dạng uốn 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 suy biến ma trận, độ cứng ảo tương ứng với bậc tự do xoay quanh trục pháp tuyến được gán xấp xỉ bằng 0,001 lần giá trị đường chéo lớn nhất trong ma trận độ cứng phần tử.

Phương pháp CS-DSG3 đóng vai trò làm trơn biến dạng uốn, cắt và màng bằng cách chia mỗi phần tử tam giác thành 3 tam giác con liên kết qua trọng tâm, từ đó loại bỏ hoàn toàn hiện tượng khóa cắt mà không làm phát sinh thêm bậc tự do giải tích. Đối với miền lỏng 3 chiều, phương trình truyền sóng Helmholtz cho lưu chất không nhớt, không xoáy và biến dạng bé được sử dụng với vận tốc âm thanh c0 và khối lượng riêng của chất lỏng. Kỹ thuật FS-FEM-T4 làm trơn gradient áp suất dựa trên các miền trơn liên kết với mặt chung của hai phần tử tứ diện 4 nút liền kề. Biên tương tác giữa chất lỏng và kết cấu vỏ tuân thủ điều kiện cân bằng áp suất động học và tính liên tục của gia tốc pháp tuyến.

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

Nguồn dữ liệu phân tích của luận văn xuất phát từ việc tính toán số trị trên 3 cấu kiện bồn chứa chuẩn mực với dải kích thước hình học đa dạng và mật độ lưới từ 100 nút đến hơn 1.200 phần tử khối. Phương pháp chọn mẫu cấu kiện tập trung vào các dạng hình học đối xứng trục đại diện cho các trạng thái ứng suất màng và ứng suất uốn phức tạp nhất trong thực tế xây dựng công nghiệp. Lý do lựa chọn mô hình phần tử hữu hạn trơn kết hợp CS-DSG3/FS-FEM là vì phương pháp này cho phép tính toán trực tiếp ma trận độ cứng trên tọa độ thực thông qua tích phân biên dọc theo các miền trơn, không phụ thuộc vào đạo hàm hàm dạng và không cần sử dụng phần tử đẳng tham số phức tạp.

Toàn bộ quy trình giải thuật dạng yếu Galerkin cho hệ phương trình tương tác rắn - lỏng được lập trình bằng ngôn ngữ tính toán kỹ thuật MATLAB. Để đánh giá độ tin cậy, các kết quả thu được được so sánh đối chứng đa chiều với các phần tử hữu hạn chuẩn như FEM-flatshell Q4 kết hợp FEM-H8, phần tử FEM-MITC4 kết hợp FEM-H8 và phần mềm ANSYS. Tiến độ nghiên cứu được triển khai đồng bộ từ khâu thiết lập công thức toán học vào tháng 02/2012, hoàn thiện thuật toán lập trình vào tháng 06/2012 và nghiệm thu bảo vệ chính thức vào cuối 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 phân tích số trị trên các mô hình bồn chứa mang lại các kết quả nổi bật sau:

Thứ nhất, về phân tích dao động tự do, phương pháp kết hợp CS-DSG3/FS-FEM cho độ chính xác vượt trội đối với 6 dạng dao động riêng đầu tiên. Tại mô hình bồn chứa nước hình trụ tròn, tần số dao động ở mode cơ bản thứ nhất đạt mức độ sai lệch chỉ 0,85% so với lời giải ANSYS, trong khi phương pháp FEM truyền thống có độ lệch lên tới 8,42%.

Thứ hai, tốc độ hội tụ của nghiệm tăng trưởng nhanh hơn khoảng 50% khi tăng mật độ lưới chia. Ứng xử quá cứng đặc trưng của phần tử tứ diện tiêu chuẩn FEM-T4 được giảm thiểu hơn 30% nhờ kỹ thuật làm trơn gradient áp suất trên mặt trong FS-FEM, giúp đường cong hội tụ tiệm cận nghiệm chính xác ngay từ các mức lưới phân chia thô.

Thứ ba, hiện tượng khóa cắt trong miền vỏ mỏng được loại bỏ hoàn toàn. Khi khảo sát tỷ lệ độ dày trên kích thước hình học nhỏ hơn 0,01, mô hình CS-DSG3 vẫn bảo đảm nghiệm ổn định tuyệt đối, không xuất hiện các mode dao động giả có năng lượng bằng không và không bị ảnh hưởng bởi thứ tự đánh số nút trên lưới tam giác méo.

Thứ tư, trong bài toán dao động cưỡng bức, hàm đáp ứng tần số tại nút 36 ở bồn trụ, nút 65 ở bồn cầu và nút 100 ở bồn nón cụt ghi nhận các đỉnh cộng hưởng trùng khớp hoàn toàn với các tần số tự nhiên đã xác định từ bài toán trị riêng, phản ánh chính xác tác động kích động từ các vị trí nguồn tải tại nút 237, 221 và 154 trong miền chất lỏng.

Thảo luận kết quả

Sở dĩ phương pháp CS-DSG3/FS-FEM đạt được độ chính xác và tính ổn định cao là nhờ cơ chế làm mềm độ cứng liên tục của các toán tử làm trơn biến dạng. Việc tích phân biến dạng trên 3 tam giác con trong CS-DSG3 giúp khôi phục các thành phần biến dạng cắt bị triệt tiêu sai lệch trong lý thuyết dầm Mindlin cổ điển. Đồng thời, miền trơn đa giác FS-FEM nối các trọng tâm khối tứ diện với mặt chung giúp phân phối trường gradient áp suất mượt mà hơn hẳn so với phép tính đạo hàm bậc nhất của phần tử tứ diện 4 nút chuẩn.

Dữ liệu tính toán trong luận văn được hệ thống hóa qua 3 bảng so sánh tần số chi tiết cho 6 mode dao động và hơn 20 đồ thị biểu diễn đường cong hội tụ cũng như phổ đáp ứng tần số biên độ. Các hình dạng mode dao động không gian 3 chiều minh họa trực quan sự dao động uốn lượn sóng của thành vỏ kết hợp với sự biến thiên áp suất thủy động trong lòng chất lỏng. So với mô hình phần tử vỏ tứ giác FEM-flatshell Q4 kết hợp phần tử lục diện 8 nút FEM-H8, phương pháp đề xuất thích ứng linh hoạt hơn trên các bề mặt cong phức tạp nhờ khả năng tự động chia lưới tam giác và tứ diện mà vẫn đảm bảo độ hội tụ vượt trội.

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

Nhằm phát huy giá trị ứng dụng của phương pháp tính toán trong thực tiễn xây dựng và thiết kế công trình, các giải pháp cụ thể được kiến nghị như sau:

Thứ nhất, tích hợp mô hình CS-DSG3/FS-FEM vào các công cụ tính toán nội bộ tại các viện thiết kế chuyên ngành xây dựng và giao thông thủy lợi. Mục tiêu là rút ngắn khoảng 30% thời gian mô phỏng kết cấu bồn bể chứa nước và xăng dầu, hoàn thành lộ trình chuẩn hóa quy trình phân tích trong vòng 12 tháng.

Thứ hai, áp dụng phương pháp phân tích tương tác vỏ và lưu chất để đánh giá lại mức độ an toàn chịu chấn động cho ít nhất 40 công trình bồn chứa công nghiệp trọng điểm tại các khu công nghiệp ven biển trong giai đoạn 2026 đến 2028 do các tổng công ty xây lắp thực hiện.

Thứ ba, xây dựng thư viện phần mềm mã nguồn mở trên nền tảng Python hoặc C++ mở rộng từ mã nguồn MATLAB của luận văn. Nhóm nghiên cứu tại các trường đại học kỹ thuật cần chủ trì thực hiện nhằm tối ưu hóa thuật toán giải ma trận thưa, hướng tới nâng cao hiệu suất xử lý lên 40% trong thời gian 18 tháng.

Thứ tư, mở rộng phạm vi nghiên cứu sang các bài toán phi tuyến hình học, xét đến hiện tượng sóng vỗ mặt tự do và tương tác kết cấu vỏ composite nhiều lớp. Cơ quan quản lý khoa học công nghệ cấp bộ cần xem xét cấp kinh phí hỗ trợ đề tài nghiên cứu chuyên sâu với tiến độ thực hiện 24 tháng.

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

Luận văn là nguồn tài liệu học thuật và kỹ thuật giá trị dành cho các nhóm độc giả sau:

  1. Học viên cao học và nghiên cứu sinh chuyên ngành Kỹ thuật Xây dựng, Cơ kỹ thuật và Cơ học tính toán. Tài liệu cung cấp tường minh cơ sở lý thuyết phần tử hữu hạn trơn, giúp người đọc nắm vững phương pháp thiết lập ma trận độ cứng cho bài toán đa vật lý tương tác rắn - lỏng.

  2. Kỹ sư kết cấu công trình dân dụng, công nghiệp và dầu khí. Luận văn cung cấp phương pháp luận chuẩn xác để tính toán dao động, kiểm tra an toàn động học cho các bồn chứa dung tích lớn, tháp làm mát và si-lô chứa chất lỏng chịu tải trọng động đất hoặc gió bão.

  3. Chuyên gia phát triển phần mềm mô phỏng kỹ thuật và phân tích phần tử hữu hạn. Người đọc có thể khai thác thuật toán CS-DSG3 và FS-FEM-T4 để tích hợp các phần tử mới vào phần mềm CAE thương mại hoặc mã nguồn mở, nâng cao năng lực xử lý lưới tự động.

  4. Giảng viên và nhà nghiên cứu tại các viện và trường đại học khối kỹ thuật. Toàn bộ sơ đồ khối, công thức giải tích và hệ thống mã lệnh MATLAB trong luận văn có thể sử dụng làm học liệu tham khảo trực tiếp cho các môn học Động lực học công trình nâng cao và Phương pháp phần tử hữu hạn.

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

Phương pháp CS-DSG3 giải quyết hiện tượng khóa cắt trong kết cấu vỏ mỏng như thế nào? Phương pháp chia mỗi phần tử tam giác thành 3 tam giác con và sử dụng kỹ thuật lệch trượt rời rạc DSG để xấp xỉ biến dạng cắt. Sau đó, toán tử làm trơn biến dạng trên toàn phần tử giúp triệt tiêu hoàn toàn độ cứng giả tạo, bảo đảm kết quả tần số chính xác ngay cả khi bề dày tấm vỏ rất nhỏ.

Tại sao phần tử tứ diện FS-FEM-T4 lại ưu việt hơn phần tử tứ diện FEM-T4 truyền thống? Phần tử FEM-T4 thông thường có độ cứng quá mức và dễ phát sinh mode năng lượng thừa. FS-FEM-T4 khắc phục triệt để nhược điểm này bằng cách thực hiện tích phân trên các miền trơn liên kết với các mặt tứ diện, giảm hơn 30% độ cứng nhân tạo và tăng đáng kể tốc độ hội tụ nghiệm áp suất.

Chất lỏng trong mô hình nghiên cứu có xét đến độ nhớt và dòng xoáy hay không? Trong phạm vi luận văn, chất lỏng được giả thiết là môi trường lý tưởng không nhớt, không xoáy, có tính nén nhẹ và dịch chuyển bé. Mô hình này hoàn toàn phù hợp và mang lại độ chính xác cao cho các bài toán phân tích dao động âm học và chấn động âm thanh của bể chứa trong dải tần số từ 0 đến 100 Hz.

Kết quả tính toán của luận văn được kiểm chứng với những phương pháp nào? Tác giả đã lập trình thuật toán trên MATLAB và so sánh trực tiếp kết quả với phần mềm ANSYS cùng các phương pháp phần tử hữu hạn chuẩn như FEM-flatshell Q4/FEM-H8 và FEM-MITC4/FEM-H8. Sai số tần số dao động ở 6 mode đầu tiên luôn duy trì ở mức dưới 1,5%.

Thuật toán trong luận văn có thể áp dụng cho các dạng hình học bồn chứa phức tạp không? Thuật toán hoàn toàn áp dụng hiệu quả cho mọi kết cấu có hình dạng bất kỳ nhờ khả năng chia lưới tự động bằng phần tử tam giác cho vỏ và phần tử tứ diện cho chất lỏng. Kỹ thuật làm trơn biến dạng giúp phương pháp không bị nhạy cảm với hiện tượng méo lưới.

Kết luận

  • Luận văn đã xây dựng thành công phương pháp số kết hợp CS-DSG3/FS-FEM cho bài toán tương tác động lực học giữa kết cấu vỏ Mindlin và chất lỏng 3 chiều.
  • Khử triệt để hiện tượng khóa cắt trong kết cấu vỏ và hiện tượng quá cứng trong miền chất lỏng mà không làm phát sinh thêm bậc tự do.
  • Đạt độ chính xác cao tương đương phần mềm ANSYS với độ sai lệch tần số dao động chỉ từ 0,5% đến 1,2% trên 6 dạng dao động đầu tiên.
  • Chứng minh tính hiệu quả vượt trội qua 3 bài toán bồn chứa điển hình gồm hình trụ tròn, hình cầu và hình nón cụt.
  • Cung cấp mã nguồn lập trình MATLAB hoàn chỉnh, tạo tiền đề vững chắc cho việc phát triển các phần mềm tính toán kết cấu chuyên dụng.

Công trình là đóng góp học thuật quan trọng, khẳng định năng lực ứng dụng phương pháp phần tử hữu hạn trơn vào các bài toán cơ học đa vật lý phức tạp. Kế hoạch nghiên cứu tiếp theo sẽ tập trung mở rộng mô hình sang phân tích phi tuyến hình học và dao động chịu tải trọng va đập trong vòng 12 tháng tới. Độc giả và các kỹ sư kết cấu hãy tham khảo chi tiết toàn văn luận văn để ứng dụng hiệu quả giải thuật tiên tiến này vào công tác tính toán và thiết kế công trình thực tế.