Tổng quan nghiên cứu

Trong kỹ thuật xây dựng công trình dân dụng và công nghiệp, thiết kế kết cấu bê tông cốt thép theo trạng thái giới hạn dẻo mang lại hiệu quả kinh tế vượt trội so với phương pháp đàn hồi truyền thống. Thực nghiệm từ công trình nhà cao tầng Cardington tại Vương quốc Anh đã chứng minh việc áp dụng lý thuyết đường chảy dẻo giúp giảm tới 37% khối lượng cốt thép chịu uốn trong sàn, cụ thể từ 23.2 tấn/tầng theo thiết kế đàn hồi xuống còn 14.5 tấn/tầng khi tính toán dẻo. Tuy nhiên, phương pháp phân tích đường chảy dẻo cổ điển bộc lộ nhiều hạn chế khi giải quyết các tấm sàn có hình học phức tạp, bản khoét lỗ hoặc có cột chống cục bộ, do kỹ sư bắt buộc phải phỏng đoán trước cơ cấu phá hoại dẻo.

Mặt khác, các phần mềm thương mại hiện nay thường dựa trên phương pháp phân tích phi tuyến từng bước với khối lượng tính toán lớn, đòi hỏi hàng trăm bước gia tải nhỏ và phụ thuộc quá nhiều vào thông số vật liệu đầu vào. Nhằm khắc phục hạn chế này, lý thuyết phân tích trực tiếp tải trọng giới hạn cận trên kết hợp với tiêu chuẩn chảy dẻo Nielsen mở ra giải pháp xác định trực tiếp khả năng chịu lực cực hạn của kết cấu mà không cần trải qua các bước trung gian. Tiêu chuẩn Nielsen phản ánh chính xác trạng thái giới hạn của tấm bê tông cốt thép chịu uốn thuần túy, nhưng lại tạo ra trở ngại toán học lớn do hàm chảy dẻo cấu thành từ giao thoa của hai hình nón bậc hai xoay, không thể áp dụng trực tiếp luật chảy dẻo kết hợp thông thường.

Được thực hiện trong giai đoạn từ tháng 07 năm 2012 đến tháng 01 năm 2013 tại Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh, nghiên cứu tập trung thiết lập công thức tính năng lượng tiêu tán dẻo qua kỹ thuật đối ngẫu điểm nội, chuyển đổi bài toán phân tích giới hạn cận trên về dạng quy hoạch nón bậc hai và giải quyết bằng phần tử hữu hạn tương thích HCT. Kết quả nghiên cứu đem lại độ chính xác cao với sai số chỉ 0.27% trên lưới 512 phần tử cho bản tựa đơn và thời gian xử lý số tính bằng giây, tạo lập công cụ hỗ trợ thiết kế tối ưu hóa kết cấu sàn thực tế.

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 xây dựng trên nền tảng lý thuyết dẻo và cơ học kết cấu tấm mỏng cổ điển Kirchhoff. Vật liệu bê tông cốt thép được giả thiết ứng xử theo mô hình cứng - dẻo lý tưởng của von Mises, bỏ qua biến dạng đàn hồi và hiện tượng tái bền vật liệu khi mô đun đàn hồi tiến tới vô cực. Tính ổn định của vật liệu được đảm bảo thông qua định đề Drucker với biểu thức tích phân công gia số biến dạng không âm trong một chu trình tải trọng khép kín.

Tiêu chuẩn chảy dẻo Nielsen cho tấm bê tông cốt thép hai lớp chịu uốn thuần túy được biểu diễn thông qua hệ phương trình mặt chảy giới hạn bởi các mô men kháng uốn giới hạn theo phương x và phương y. Do mặt chảy Nielsen tạo bởi hai hình nón bậc hai xoay trong không gian ứng suất ba chiều của mô men uốn và mô men xoắn, nghiên cứu áp dụng thuật toán đối ngẫu điểm nội để biểu diễn hàm tiêu tán năng lượng dẻo dưới dạng tổng các thành phần biến dạng độ cong khả dĩ động thuộc nón bậc hai tự đối ngẫu. Định lý động học cận trên chỉ ra rằng tải trọng giới hạn tìm được là giá trị nhỏ nhất của công suất tiêu tán năng lượng dẻo nội tại khi công suất ngoại lực được trực giao hóa về đơn vị.

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 phần tử hữu hạn tương thích C1 Hsieh-Clough-Tocher nhằm xấp xỉ trường vận tốc chuyển vị võng và trường góc xoay của tấm sàn. Mỗi phần tử tam giác HCT gồm 12 bậc tự do chính (gồm độ võng và hai góc xoay tại 3 nút đỉnh, cùng các góc xoay pháp tuyến tại 3 điểm giữa cạnh), được phân chia nội bộ thành 3 tam giác con đồng quy tại trọng tâm. Ba bậc tự do tại nút trọng tâm được khử tĩnh thông qua ma trận tương thích liên tục, giúp tăng độ chính xác tương đương 3 phần tử thông thường mà không làm tăng bậc tự do tổng thể của hệ thống.

Quá trình tích phân số trong từng tam giác con được rời rạc hóa qua 3 điểm Gauss. Toàn bộ thuật toán được lập trình trên môi trường MATLAB phiên bản R2009a, tích hợp công cụ giải quy hoạch toán học MOSEK phiên bản 5.0 trên hệ thống phần cứng vi xử lý 2.0 GHz. Cỡ mẫu khảo sát bao gồm 4 cấu hình hình học kết cấu tiêu biểu: tấm vuông, tấm chữ nhật, tấm tròn và tấm chữ L với các cấp độ chia lưới phần tử tăng dần từ 8, 32, 128 đến 512 phần tử HCT nhằm kiểm tra sự hội tụ và tính ổn định của thuật toán. Phương pháp phân tích trực tiếp được lựa chọn vì loại bỏ hoàn toàn sự tích lũy sai số qua các bước lặp thời gian trong phân tích phi tuyến từng bước 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 số trên các bài toán chuẩn đã chứng minh tính ưu việt vượt trội của việc kết hợp phần tử HCT và tiêu chuẩn Nielsen:

Thứ nhất, đối với bài toán tấm sàn hình vuông kích thước 5m x 5m, chiều dày 100mm, bê tông có cường độ nén 28 MPa và cốt thép D12 khoảng cách 100mm có giới hạn chảy 280 MPa với điều kiện bốn biên tựa đơn, hệ số tải trọng giới hạn hội tụ nhanh chóng về nghiệm giải tích chính xác của Prager là 24.00. Khi tăng số lượng phần tử HCT từ 8, 32, 128 lên 512 phần tử, hệ số tải trọng giảm dần từ 27.12 (sai số 13.00%), 24.80 (sai số 3.35%), 24.22 (sai số 0.91%) xuống mức 24.06 với sai số chỉ còn 0.27%. Thời gian tính toán cho lưới 512 phần tử chỉ mất 19 giây.

Thứ hai, đối với bài toán tấm hình vuông bốn biên ngàm chịu tải trọng phân bố đều, hệ số tải trọng giới hạn thu được ở lưới 512 phần tử đạt giá trị 46.49 với sai số 8.49% so với nghiệm giải tích của Fox là 42.851. Giá trị này chặt chẽ hơn đáng kể so với các kết quả phân tích cận trên sử dụng tiêu chuẩn von Mises của Hodge và Belytschko (49.25) hay nghiên cứu của nhóm tác giả Canh và cộng sự năm 2010 (46.71).

Thứ ba, về hiệu năng tính toán số, phương pháp đề xuất vượt trội so với các phương pháp không lưới. Ở cùng điều kiện bài toán tấm biên tựa, phương pháp phần tử không lưới EFG trên hệ lưới 14x14 điểm đòi hỏi 42 giây để đạt hệ số 25.01, trong khi phần tử HCT với 128 phần tử chỉ mất đúng 1 giây để đạt hệ số 24.22.

Thứ tư, biểu đồ phân bố trường chuyển vị độ võng và các trường góc xoay thu được từ nghiệm cận trên phản ánh chính xác vị trí phát triển các đường khớp dẻo tại giao điểm hai đường chéo tấm sàn, hoàn toàn trùng khớp với cơ cấu phá hoại thực nghiệm của Johansen trên tấm mẫu 2m x 2m và kết quả mô phỏng bằng chương trình DLO.

Thảo luận kết quả

Các kết quả tính toán có thể được biểu diễn trực quan thông qua bảng đối chiếu độ hội tụ và các đồ thị quan hệ giữa hệ số tải trọng với mật độ lưới chia. Biểu đồ đường cong hội tụ thể hiện rõ sai số giảm theo hàm mũ khi số lượng điểm Gauss và số phần tử HCT tăng lên.

Sự khác biệt về độ chính xác giữa trường hợp biên tựa đơn (sai số 0.27%) và biên ngàm (sai số 8.49%) bắt nguồn từ cơ chế tiêu tán năng lượng dẻo. Ở điều kiện biên ngàm, năng lượng tiêu tán tập trung với mật độ rất cao dọc theo suốt chiều dài các cạnh ngàm do sự xuất hiện của mô men uốn âm cực đại. Khi xấp xỉ chuyển vị liên tục bậc ba, việc rời rạc hóa vùng biên ngàm làm tích lũy sai số cục bộ lớn hơn so với vùng giữa nhịp của bản tựa đơn. Việc áp dụng tiêu chuẩn Nielsen trong mô hình quy hoạch nón bậc hai tránh được hiện tượng tuyến tính hóa mặt chảy dẻo vốn làm sai lệch kết quả từ 5% đến 15% trong các nghiên cứu cổ điển. Mô hình chuyển vị 3D và các bản đồ đẳng mức góc xoay cung cấp bức tranh toàn diện về động học phá hủy, cho phép kỹ sư xác định chính xác vị trí cần gia cường cốt thép chống uốn và chống nứt trên mặt sàn.

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

Dựa trên các kết quả phân tích số và cơ sở lý thuyết đã được kiểm chứng, nghiên cứu đưa ra 4 khuyến nghị then chốt nhằm ứng dụng và mở rộng phương pháp trong thực tiễn:

  1. Tích hợp thuật toán quy hoạch nón bậc hai vào các phần mềm thiết kế kết cấu chuyên dụng: Các đơn vị phát triển phần mềm xây dựng trong nước cần nghiên cứu tích hợp module giải SOCP kết hợp phần tử HCT vào các công cụ CAD/CAE hiện hành nhằm rút ngắn từ 60% đến 80% thời gian tính toán kết cấu sàn phức tạp trong lộ trình 12 đến 24 tháng tới.

  2. Ứng dụng thiết kế dẻo cho các công trình nhà nhiều tầng và sàn khẩu độ lớn: Các kỹ sư thiết kế kết cấu tại các viện quy hoạch và công ty tư vấn xây dựng cần chủ động áp dụng phương pháp phân tích giới hạn trực tiếp cho hệ sàn nấm, sàn phẳng không dầm hoặc sàn có khoét lỗ kỹ thuật nhằm mục tiêu cắt giảm từ 20% đến 35% hàm lượng cốt thép chịu uốn mà vẫn đảm bảo tuyệt đối hệ số an toàn chịu lực.

  3. Phát triển kỹ thuật chia lưới tự thích nghi (Adaptive Mesh Refinement): Các nhóm nghiên cứu chuyên sâu về cơ học tính toán cần tiếp tục hoàn thiện giải thuật tự động làm mịn lưới phần tử cục bộ dọc theo các đường dẻo và biên ngàm trong vòng 6 đến 12 tháng, nhằm mục tiêu giảm sai số tính toán của bài toán biên ngàm từ mức 8.49% xuống dưới 3.0%.

  4. Cập nhật và hoàn thiện các tiêu chuẩn thiết kế quốc gia: Đề nghị các hội đồng chuyên ngành thuộc Bộ Xây dựng xem xét bổ sung các hướng dẫn tính toán trạng thái giới hạn dẻo theo tiêu chuẩn Nielsen vào hệ thống Tiêu chuẩn Xây dựng Việt Nam (TCVN) trong chu kỳ soát xét 3 năm tới, tạo hành lang pháp lý vững chắc cho việc tối ưu hóa vật liệu xây dựng.

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

Nội dung và phương pháp luận của luận văn mang lại giá trị học thuật và ứng dụng thực tiễn cao cho 4 nhóm đối tượng cụ thể:

  • Kỹ sư thiết kế kết cấu công trình dân dụng và công nghiệp: Tài liệu là cẩm nang hữu ích giúp kỹ sư nắm bắt phương pháp đánh giá chính xác sức chịu tải giới hạn của bản sàn bê tông cốt thép có hình dạng bất kỳ, sàn khoét lỗ hoặc chịu tải tập trung, từ đó tối ưu hóa đường kính và khoảng cách bố trí thép sàn trong các dự án thực tế.

  • Học viên cao học, nghiên cứu sinh và giảng viên ngành Kỹ thuật Xây dựng - Cơ học kết cấu: Cung cấp tài liệu tham khảo chuẩn mực về cách thiết lập công thức biến phân, kỹ thuật đối ngẫu điểm nội trong tối ưu hóa phi tuyến và quy trình lập trình phương pháp phần tử hữu hạn bậc cao trên môi trường MATLAB.

  • Các chuyên gia phát triển phần mềm tính toán kết cấu (CAE): Nguồn tham khảo giá trị về thuật toán chuyển đổi bài toán cơ học phi tuyến phức tạp về dạng tối ưu hình nón bậc hai (SOCP) kết hợp bộ giải tối ưu MOSEK để tăng tốc độ tính toán cho các phần mềm thương mại.

  • Chủ đầu tư, nhà thầu thi công và các đơn vị thẩm tra thiết kế: Cung cấp cơ sở khoa học để đánh giá phương án kinh tế - kỹ thuật, hỗ trợ thẩm định khả năng chịu lực thực tế của kết cấu sàn nhằm tiết kiệm chi phí vật tư thép xây dựng từ 20% đến 37% trên tổng mức đầu tư phần thô.

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

Tại sao tiêu chuẩn Nielsen lại phù hợp hơn tiêu chuẩn von Mises cho tấm sàn bê tông cốt thép?

Tiêu chuẩn von Mises vốn được xây dựng cho vật liệu kim loại đẳng hướng chịu dẻo thuần túy, không phản ánh đúng tính chất làm việc dị hướng và khả năng chịu kéo nén khác biệt của bê tông cốt thép. Tiêu chuẩn Nielsen thiết lập mặt chảy dẻo dựa trên mô men uốn kháng nứt thực tế theo hai phương đặt thép độc lập, thể hiện chính xác mối quan hệ tương tác giữa mô men uốn và mô men xoắn trong tấm bản chịu lực phẳng.

Phương pháp phân tích trực tiếp tải trọng giới hạn có ưu điểm gì so với phân tích phi tuyến từng bước?

Phân tích từng bước đòi hỏi chia tải trọng thành hàng trăm gia số nhỏ và giải lặp Newton-Raphson tại từng bước, gây tốn kém tài nguyên tính toán và dễ mất hội tụ khi kết cấu bước vào giai đoạn chảy dẻo cục bộ. Phương pháp phân tích giới hạn trực tiếp giải một bài toán quy hoạch toán học duy nhất để xác định ngay hệ số tải trọng phá hoại mà không cần theo dõi lịch sử biến dạng trung gian, rút ngắn thời gian xử lý từ hàng chục phút xuống vài giây.

Vì sao phần tử HCT lại được ưu tiên lựa chọn để xấp xỉ trường chuyển vị của tấm?

Phần tử tam giác Hsieh-Clough-Tocher là phần tử tương thích liên tục C1 sử dụng đa thức bậc ba, đảm bảo tính liên tục của cả độ võng lẫn góc xoay trên toàn bộ các biên tiếp giáp phần tử. Với 12 bậc tự do sau khi đã khử nút trong qua ma trận tương thích, phần tử HCT mang lại độ chính xác cao tương đương 3 phần tử tam giác thông thường nhưng giảm tải đáng kể số lượng biến số cần giải trong bộ tối ưu.

Kỹ thuật quy hoạch nón bậc hai (SOCP) giải quyết trở ngại toán học nào trong bài toán cận trên?

Mặt chảy dẻo Nielsen có các điểm kỳ dị và đường giao thoa không khả vi liên tục, khiến việc lấy đạo hàm theo luật chảy dẻo kết hợp bị bế tắc. Kỹ thuật đối ngẫu điểm nội cho phép biến đổi hàm năng lượng tiêu tán dẻo phi tuyến về dạng điều kiện ràng buộc nón bậc hai xoay chuẩn tắc, giúp các thuật toán điểm nội giải quyết bài toán với độ hội tụ chắc chắn và tốc độ tối ưu.

Nguyên nhân nào khiến sai số của bài toán tấm biên ngàm (8.49%) cao hơn bài toán tấm biên tựa (0.27%)?

Tại biên ngàm, mật độ tiêu tán năng lượng dẻo tập trung đột biến dọc theo các cạnh cố định do mô men uốn âm đạt giá trị cực hạn. Lưới phần tử đồng đều thông thường chưa phản ánh trọn vẹn sự biến thiên độ cong rất lớn tại sát mép ngàm, dẫn đến sai số tích lũy cao hơn so với bản biên tựa đơn nơi năng lượng dẻo phân bố đều hơn về phía trung tâm tấm sàn.

Kết luận

  • Luận văn đã xây dựng thành công quy trình số hoàn chỉnh phân tích giới hạn cận trên cho tấm sàn bê tông cốt thép dựa trên tiêu chuẩn chảy dẻo Nielsen và phương pháp phần tử hữu hạn tương thích HCT.
  • Thiết lập công thức tính toán năng lượng tiêu tán dẻo thông qua kỹ thuật đối ngẫu điểm nội, chuyển đổi bài toán cơ học phi tuyến về dạng quy hoạch nón bậc hai giải bằng công cụ MOSEK.
  • Kết quả mô phỏng số đạt độ chính xác cao với sai số hội tụ chỉ 0.27% trên bài toán tấm vuông biên tựa và 8.49% trên bài toán tấm vuông biên ngàm, vượt trội so với các mô hình tiêu chuẩn von Mises trước đây.
  • Tối ưu hóa triệt để thời gian tính toán số, đạt lời giải chính xác chỉ trong khoảng thời gian từ 1 đến 19 giây trên máy tính cá nhân tiêu chuẩn, chứng minh tính khả thi trong ứng dụng thực tế.
  • Khẳng định tiềm năng ứng dụng lý thuyết dẻo vào thực tiễn thiết kế xây dựng nhằm tiết kiệm tới 37% khối lượng cốt thép chịu uốn cho kết cấu sàn bê tông cốt thép nhà cao tầng.

Đóng góp then chốt của công trình là đã tháo gỡ điểm nghẽn toán học trong việc tính toán năng lượng tiêu tán dẻo của tiêu chuẩn Nielsen, cung cấp công cụ tính toán cận trên đáng tin cậy cho kỹ thuật xây dựng. Để tiếp tục phát triển nghiên cứu, lộ trình tiếp theo cần tập trung tích hợp kỹ thuật lưới tự thích nghi (Adaptive FEM) và mở rộng cho kết cấu sàn bê tông ứng suất trước trong vòng 12 đến 24 tháng tới. Các kỹ sư và nhà nghiên cứu quan tâm có thể ứng dụng ngay khung thuật toán này để tối ưu hóa thiết kế kết cấu bê tông cốt thép cho các dự án xây dựng hiện đại.