Tổng quan nghiên cứu

Trong kỹ thuật công trình hiện đại, việc đánh giá chính xác sức chịu tải giới hạn của kết cấu tấm chịu lực đóng vai trò quyết định đến mức độ an toàn và tối ưu hóa chi phí vật liệu. Phương pháp phân tích giới hạn trực tiếp cho phép xác định tải trọng phá hoại dẻo mà không cần thông qua các bước phân tích đàn hồi - dẻo từng bước phức tạp, giúp tiết kiệm đáng kể thời gian tính toán. Tuy nhiên, khi áp dụng lý thuyết tấm dày Mindlin - Reissner cho các kết cấu có tỷ số kích thước chiều dài trên chiều dày dao động rộng từ 1 đến 500, các phương pháp phần tử hữu hạn truyền thống thường gặp phải hiện tượng khóa cắt nghiêm trọng, làm sai lệch kết quả tải trọng tới hạn lên đến hàng chục phần trăm.

Công trình nghiên cứu này được thực hiện tại Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh trong giai đoạn 2011 – 2012, nhằm mục tiêu phát triển mô hình tính toán số tiên tiến kết hợp giữa phương pháp phần tử hữu hạn trơn dựa trên cạnh và kỹ thuật rời rạc lệch trượt cùng tối ưu hóa nón bậc hai. Nghiên cứu tập trung giải quyết bài toán tấm dày chịu tải trọng phân bố đều và tải trọng tập trung với vật liệu có giới hạn chảy danh định 250 MPa. Kết quả nghiên cứu mang lại ý nghĩa học thuật và thực tiễn sâu sắc khi nâng cao độ chính xác của nghiệm cận trên với mức sai số đối chiếu chỉ khoảng 0,2% đến 4,97% so với các công bố quốc tế, đồng thời loại bỏ hoàn toàn hiện tượng khóa cắt, tạo tiền đề nâng cao hiệu suất tính toán kết cấu tấm trong 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

Nghiên cứu được xây dựng trên nền tảng tích hợp của bốn lý thuyết và mô hình cơ học kết cấu then chốt:

  • Lý thuyết tấm dày Mindlin - Reissner: Khác với giả thiết Kirchhoff bỏ qua biến dạng trượt ngang, lý thuyết Mindlin - Reissner tính đến ảnh hưởng của góc xoay độc lập và biến dạng trượt ngoài mặt phẳng. Mô hình này cho phép phân tích chính xác các tấm có chiều dày thay đổi và tấm dày chịu uốn.
  • Lý thuyết biến dạng cắt bậc nhất: Thiết lập mối quan hệ chuyển vị tuyến tính theo chiều dày tấm, trong đó các đoạn thẳng pháp tuyến với mặt trung hòa trước biến dạng vẫn giữ thẳng nhưng không nhất thiết phải vuông góc với mặt trung hòa sau biến dạng.
  • Tiêu chuẩn chảy dẻo Von Mises: Đóng vai trò là điều kiện giới hạn đàn hồi của vật liệu đẳng hướng, giả định sự chảy dẻo xảy ra khi bất biến thứ hai của tensor ứng suất lệch đạt tới giá trị giới hạn chịu kéo chia cho căn bậc hai của 3.
  • Định lý động phân tích giới hạn cận trên: Xác định hệ số tải trọng phá hoại dẻo thông qua việc cực tiểu hóa hàm tiêu tán năng lượng biến dạng dẻo nội năng với điều kiện cân bằng công ngoại lực bằng 1.

Các khái niệm trọng tâm bao gồm trường vận tốc suy rộng, biến dạng trơn trên miền cạnh, hàm suy biến dẻo nón bậc hai và năng lượng tiêu tán dẻo cục bộ.

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

Mô hình nghiên cứu sử dụng dữ liệu số học mô phỏng cho các kết cấu tấm vuông, chữ nhật và tròn với kích thước nhịp điển hình 10 m và chiều dày biến thiên theo tỷ lệ nhịp trên chiều dày từ 1 đến 500. Về kỹ thuật chia lưới và cỡ mẫu phân tích số, tác giả thiết lập khảo sát trên một phần tư miền tấm đối xứng bằng các lưới phần tử tam giác ba nút có độ mịn tăng dần lần lượt là 7x7, 9x9 và 11x11 phần tử, tương ứng với hàng trăm cạnh trơn độc lập.

Phương pháp phân tích được lựa chọn là sự tích hợp giữa phương pháp phần tử hữu hạn trơn dựa trên cạnh và kỹ thuật rời rạc lệch trượt. Phương pháp này chia miền tính toán thành các miền con dựa trên cạnh chung của phần tử tam giác, giúp làm mềm độ cứng giả tạo của lưới và tính toán chính xác biến dạng uốn, cắt. Tích phân theo chiều dày tấm được thực hiện chính xác với 3 điểm Gauss. Toàn bộ bài toán tối ưu phi tuyến có hàm mục tiêu không khả vi được chuyển đổi tương đương về bài toán quy hoạch nón bậc hai dạng chuẩn, sau đó được giải tự động bằng phần mềm tối ưu Mosek tích hợp trong môi trường lập trình Matlab trên hệ thống vi xử lý 2,2 GHz với bộ nhớ 4 GB RAM. Sự kết hợp này được lựa chọn vì khắc phục triệt để nhược điểm không khả vi của hàm tiêu tán dẻo trong các thuật toán gradient phi tuyến truyền thống.

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 và thực nghiệm số mang lại các kết quả định lượng nổi bật:

  • Triệt tiêu hoàn toàn hiện tượng khóa cắt: Khi tỷ số nhịp trên chiều dày tăng mạnh từ 10 lên 100 và đạt cực hạn 500, giá trị tải trọng giới hạn của tấm vuông biên tựa đơn hội tụ mượt mà về giá trị không đổi 25,07 mà không xuất hiện bất kỳ sự suy giảm độ cứng hay biến động bất thường nào.
  • Độ chính xác cao trong phân tích tấm vuông: Với lưới phân tử 11x11, nghiệm cận trên cho tải trọng phá hoại dẻo của tấm vuông biên tựa đơn chịu tải phân bố đều đạt 25,07, chỉ chênh lệch khoảng 0,2% so với kết quả giải tích 25,02 của các công bố quốc tế kinh điển. Đối với trường hợp tấm vuông liên kết ngàm bốn cạnh, hệ số tải trọng giới hạn đạt 47,54, chênh lệch 4,97% so với mức 45,29 của phương pháp truyền thống.
  • Tính hội tụ đơn điệu theo độ mịn lưới: Khảo sát trên lưới 7x7, 9x9 và 11x11 cho thấy giá trị hàm mục tiêu tiệm cận đều đặn từ trên xuống. Cụ thể ở tấm ngàm có tỷ số nhịp trên chiều dày bằng 10, hệ số tải trọng chuyển biến từ 30,35 ở lưới 7x7 xuống 30,06 ở lưới 9x9 và đạt 29,91 ở lưới 11x11.
  • Hiệu năng xử lý tối ưu của thuật toán nón bậc hai: Thời gian xử lý cho mỗi bài toán tấm phức tạp với hàng nghìn biến số hình nón chỉ mất khoảng vài chục giây trên cấu hình máy tính tiêu chuẩn, chứng minh tốc độ hội tụ vượt trội so với quy hoạch phi tuyến cổ điển.

Thảo luận kết quả

Nguyên nhân chính giúp phương pháp đạt độ chính xác cao là nhờ kỹ thuật làm trơn biến dạng trên từng cạnh phần tử tam giác, giúp giảm độ cứng quá mức vốn có của phương pháp phần tử hữu hạn tiêu chuẩn. Đồng thời, kỹ thuật rời rạc lệch trượt chuyển đổi các ma trận biến dạng cắt từ dạng hàm tuyến tính phức tạp về dạng hằng số ổn định trên mỗi miền trơn, ngăn chặn sự tích tụ năng lượng cắt ảo khi chiều dày tấm tiến về 0.

Dữ liệu tính toán được minh họa trực quan thông qua các bảng đối chuẩn đa điều kiện biên (tựa đơn, ngàm, hai ngàm - hai tựa đơn, ba ngàm - một tự do) và biểu đồ hội tụ năng lượng tiêu tán dẻo. Biểu đồ đường đẳng mức phân bố năng lượng tiêu tán dẻo trên một phần tư miền tấm phản ánh chính xác cơ chế khớp dẻo: năng lượng tập trung dày đặc tại các góc ngàm và tâm tấm, tương thích hoàn hảo với các đường phá hoại dẻo hình học trong thực nghiệm cơ học kết cấu. So với mô hình tấm mỏng Kirchhoff, mô hình Mindlin - Reissner sử dụng phương pháp này luôn cung cấp một chặn trên an toàn và thực tế hơn cho sức chịu tải của cấu kiện công trình.

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

  • Tích hợp mô hình phần tử hữu hạn trơn vào các phần mềm thương mại CAE: Các đơn vị phát triển phần mềm tính toán kết cấu cần lập trình bổ sung module phần tử tam giác trơn dựa trên cạnh để nâng cao độ tin cậy khi mô phỏng bản sàn bê tông cốt thép dày và tấm đáy móng, mục tiêu giảm từ 10% đến 15% lượng thép bố trí dư thừa do tính toán sai lệch trong vòng 12 đến 24 tháng tới.
  • Ứng dụng quy hoạch nón bậc hai cho bài toán kiểm định kết cấu thực tế: Các kỹ sư thiết kế và chuyên gia kiểm định công trình nên áp dụng thuật toán tối ưu nón bậc hai để đánh giá sức chịu tải dư của các kết cấu dầm, sàn chịu tải trọng vượt tải định mức, giúp rút ngắn thời gian xử lý dữ liệu từ vài ngày xuống dưới 60 phút.
  • Chuẩn hóa quy trình tính toán phân tích giới hạn trong tiêu chuẩn xây dựng: Các viện nghiên cứu chuyên ngành và cơ quan quản lý tiêu chuẩn kỹ thuật cần xây dựng bộ hướng dẫn tính toán dẻo cho tấm và vỏ composite, áp dụng thí điểm tại các dự án nhà cao tầng và công trình ngầm trong giai đoạn 2026 – 2028.
  • Nâng cấp cấu hình phần cứng và tối ưu hóa bộ nhớ: Các phòng thí nghiệm cơ học tính toán cần trang bị hệ thống máy chủ xử lý song song với dung lượng bộ nhớ tối thiểu 32 GB RAM để mở rộng quy mô khảo sát cho các kết cấu tấm có mạng lưới phần tử lớn hơn 50x50 nút.

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

  • Kỹ sư thiết kế kết cấu công trình: Nắm vững nguyên lý ứng xử dẻo của tấm dày chịu uốn và cắt để tối ưu hóa thiết kế bản sàn, móng bè, đài cọc và vách hầm trong các công trình dân dụng lẫn công nghiệp.
  • 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: Sử dụng làm tài liệu tham khảo chuyên sâu về kỹ thuật làm trơn hóa biến dạng, cơ sở toán học của thuật toán tối ưu nón bậc hai và phương pháp lập trình mô phỏng trên nền tảng Matlab.
  • Các nhà phát triển phần mềm mô phỏng kỹ thuật: Khai thác chi tiết các công thức giải tích ma trận hàm dạng trơn và thuật toán biến đổi biến nón để tích hợp vào các lõi tính toán phần tử hữu hạn thế hệ mới.
  • Chuyên gia thẩm tra và kiểm định an toàn công trình: Vận dụng phương pháp phân tích cận trên để đánh giá nhanh chóng và an toàn khả năng chịu tải phá hoại của các kết cấu hiện hữu trước nguy cơ quá tải hoặc chịu tác động ngoại lực đặc biệt.

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

Hiện tượng khóa cắt trong phân tích tấm dày là gì và ảnh hưởng thế nào đến kết quả?

Hiện tượng khóa cắt xuất hiện khi chiều dày tấm tiến dần về 0 nhưng ma trận biến dạng cắt của phần tử bậc thấp vẫn giữ năng lượng lớn không triệt tiêu. Điều này làm cho kết cấu trở nên cứng giả tạo, dẫn đến độ võng tính toán nhỏ hơn thực tế và hệ số tải trọng phá hoại bị đánh giá sai lệch nghiêm trọng.

Vì sao phương pháp phần tử hữu hạn trơn dựa trên cạnh lại vượt trội hơn phương pháp truyền thống?

Phương pháp làm trơn trên cạnh thực hiện tích phân trên các miền con nối tâm phần tử và nút chung, giúp làm mềm độ cứng mô hình và triệt tiêu tính phụ thuộc quá mức vào hình dạng lưới. Nhờ đó, phương pháp mang lại độ chính xác cao hơn, hội tụ nhanh hơn khoảng 15% đến 20% so với phần tử tiêu chuẩn trên cùng mật độ lưới.

Thuật toán tối ưu nón bậc hai giải quyết vấn đề gì trong phân tích giới hạn?

Hàm tiêu tán năng lượng dẻo thường chứa các điểm không khả vi tại nơi không có biến dạng dẻo, gây thất bại cho các thuật toán gradient. Quy hoạch nón bậc hai chuyển bài toán về dạng biến nón chuẩn tắc, giúp tìm nghiệm tối ưu toàn cục nhanh chóng, chính xác và ổn định tuyệt đối.

Luận văn đã kiểm chứng phương pháp với những điều kiện biên nào?

Nghiên cứu đã khảo sát bốn trường hợp liên kết biên bao gồm: biên tựa đơn toàn phần, biên ngàm tuyệt đối bốn cạnh, biên hỗn hợp hai ngàm kết hợp hai tựa đơn, và biên ba ngàm kết hợp một cạnh tự do. Tất cả các trường hợp đều cho thấy sự phân bố năng lượng dẻo chính xác theo lý thuyết phá hoại.

Phương pháp này có áp dụng được cho vật liệu phi tuyến phức tạp khác không?

Có. Phương pháp hoàn toàn có khả năng mở rộng cho các loại vật liệu dị hướng, vật liệu composite nhiều lớp hoặc vật liệu có tính năng biến thiên liên tục. Bằng cách thiết lập lại ma trận vật liệu và tiêu chuẩn dẻo tương ứng dưới dạng nón bậc hai, bài toán vẫn được giải quyết hiệu quả.

Kết luận

  • Xây dựng thành công quy trình phân tích giới hạn cận trên cho tấm dày Mindlin - Reissner bằng sự kết hợp giữa phần tử hữu hạn trơn dựa trên cạnh và kỹ thuật rời rạc lệch trượt.
  • Giải quyết triệt để bài toán khóa cắt và bài toán không khả vi của hàm tiêu tán dẻo thông qua ứng dụng hiệu quả thuật toán tối ưu nón bậc hai.
  • Đạt độ chính xác cao và hội tụ đơn điệu qua các ví dụ số tấm vuông, tấm chữ nhật và tấm tròn với sai số chỉ từ 0,2% đến 4,97% so với các công bố quốc tế.
  • Khẳng định tính ưu việt của mô hình tấm dày Mindlin - Reissner như một chặn trên an toàn và sát thực tế hơn so với lý thuyết tấm mỏng cổ điển.
  • Mở ra hướng nghiên cứu tiếp theo trong giai đoạn tới nhằm phát triển mô hình cho kết cấu tấm vỏ composite phức tạp chịu tải trọng động lực học.

Quý độc giả, kỹ sư và nhà nghiên cứu quan tâm đến mô phỏng số và tối ưu hóa kết cấu dẻo hãy ứng dụng ngay mô hình này vào công tác tính toán thực tế để nâng cao chất lượng và độ an toàn cho các công trình xây dựng.