Tổng quan nghiên cứu

Theo thống kê trong ngành xây dựng, các công trình sử dụng kết cấu bê tông cốt thép chiếm tới khoảng 70% tổng số công trình xây dựng tại Việt Nam nhờ khả năng chịu lực vượt trội và độ bền cao. Tuy nhiên, bê tông là loại vật liệu có cường độ chịu kéo rất thấp, chỉ bằng khoảng 1/15 đến 1/8 cường độ chịu nén, do đó sự làm việc liên hợp và khả năng truyền ứng suất thông qua lực dính bám giữa bê tông và cốt thép là yếu tố mang tính quyết định đến an toàn chịu lực của toàn bộ kết cấu.

Trong tính toán mô phỏng kết cấu bằng phương pháp phần tử hữu hạn hiện nay, đa số các mô hình thường áp đặt giả thiết đơn giản hóa rằng cốt thép dính bám tuyệt đối vào bê tông. Giả thiết này đã tạo ra sự chênh lệch đáng kể giữa kết quả tính toán số và số liệu thực nghiệm thực tế, với sai số chuyển vị và tải trọng phá hoại có thể lên tới 18% đến 25% khi cấu kiện xuất hiện vết nứt và biến dạng dẻo.

Mục tiêu cốt lõi của luận văn là thiết lập quy trình mô phỏng số phi tuyến chính xác trên phần mềm ABAQUS nhằm đánh giá sự bám dính không hoàn hảo giữa bê tông và cốt thép khi cấu kiện dầm bị phá hoại. Nghiên cứu tập trung so sánh ứng xử tải trọng - chuyển vị tại giữa dầm trên ba nền vật liệu khác nhau: bê tông thông thường, bê tông geopolymer thân thiện môi trường và bê tông cốt liệu xỉ thép.

Nghiên cứu được triển khai thực nghiệm và mô phỏng tại Trường Đại học Sư phạm Kỹ thuật Thành phố Hồ Chí Minh vào năm 2017, áp dụng trên các cấu kiện dầm chịu uốn có kích thước mặt cắt 200 x 300 mm và chiều dài 3300 mm. Kết quả nghiên cứu giúp giảm thiểu sai số mô phỏng xuống dưới mức 6%, nâng cao độ tin cậy của mô hình số và giúp các đơn vị thiết kế tiết kiệm tới 40% chi phí thử nghiệm 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 được xây dựng dựa trên sự kết hợp của hai mô hình cơ học vật liệu phi tuyến kinh điển. Đầu tiên là mô hình nén phi tuyến của bê tông do Hsu – Hsu (1994) đề xuất, cho phép mô tả chính xác quan hệ ứng suất - biến dạng của bê tông chịu nén từ giai đoạn đàn hồi ban đầu đến khi đạt cường độ cực hạn và suy giảm độ cứng với cấp cường độ bê tông lên tới 60 MPa. Thứ hai là mô hình độ cứng chịu kéo của bê tông được Nayal và Rasheed (2006) hiệu chỉnh và phát triển để tích hợp hoàn hảo trong phần mềm ABAQUS, giúp giải quyết triệt để lỗi phân kỳ toán học khi ứng suất kéo suy giảm xuống 0.77 lần ứng suất nứt cực hạn tại biến dạng cực hạn 125 lần biến dạng nứt.

Các khái niệm cơ sở được áp dụng bao gồm cơ chế truyền lực qua ba thành phần: lực liên kết dính bám hóa học của hồ xi măng, lực ma sát cơ học trên mặt tiếp xúc và hiệu ứng chốt khóa giữa gân thép với khuôn bê tông xung quanh. Ứng xử vật liệu cốt thép được mô tả theo mô hình đàn dẻo lý tưởng một chiều SEPL và IEPL với mô đun đàn hồi 210000 MPa, cường độ chảy tiêu chuẩn 365 MPa và cường độ bền kéo 440 MPa. Cơ chế bám dính mặt tiếp xúc được mô hình hóa bằng lý thuyết dính kết Cohesive Surface kết hợp hệ số dính nhớt Viscosity nhằm phản ánh hiện tượng trượt tuột vật lý khi cấu kiện bị nứt vỡ.

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

Dữ liệu nghiên cứu được thu thập từ chuỗi thí nghiệm thực tế tại phòng thí nghiệm chuyên ngành, bao gồm 3 tổ mẫu thí nghiệm kéo tuột cốt thép gân có đường kính 12 mm, 16 mm, 20 mm chôn sâu trong các khối bê tông hình trụ tiêu chuẩn kích thước 150 x 150 x 150 mm, cùng với các mẫu dầm bê tông kích thước chuẩn 200 x 300 x 3300 mm. Phương pháp chọn mẫu có chủ đích được áp dụng nhằm thu thập đầy đủ thông số cơ lý cho 3 nền vật liệu đặc thù là bê tông thông thường, bê tông xỉ thép và bê tông geopolymer.

Phương pháp phân tích phần tử hữu hạn phi tuyến 3D trên phần mềm ABAQUS/Standard được lựa chọn vì đây là công cụ mạnh mẽ nhất hiện nay có khả năng phân tích tương tác tiếp xúc phi tuyến bề mặt và mô phỏng chính xác sự xuất hiện vết nứt cục bộ. Thân dầm bê tông và cốt thép dọc chịu lực được rời rạc hóa bằng phần tử khối 8 nút giảm tích phân C3D8R, trong khi hệ cốt đai chịu cắt đường kính 6 mm bước đai 150 mm được mô hình hóa bằng phần tử thanh T3D2. Toàn bộ quá trình lập trình số, tối ưu hóa mạng lưới phần tử và phân tích so sánh được thực hiện trong khung thời gian từ tháng 3 năm 2017 đến tháng 9 năm 2017.

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

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

Thứ nhất, qua phân tích mô phỏng thí nghiệm kéo tuột cốt thép trên nền bê tông geopolymer, lực kéo tuột lớn nhất và năng lượng phá hủy diện tích bề mặt tăng tỷ lệ thuận theo đường kính thanh thép từ 12 mm đến 20 mm. Dựa trên kết quả này, tác giả đã nội suy thành công các thông số dính bám tiếp xúc cho cốt thép chịu kéo đường kính 14 mm với diện tích mặt cắt 154 mm2, xác lập lực dính bám cực hạn đạt 68.5 kN và năng lượng phá hủy đạt mức 3.8 N/mm.

Thứ hai, đối với dầm bê tông geopolymer gia cường 5 thanh thép kéo đường kính 14 mm với tổng diện tích cốt thép 770 mm2, mô hình mô phỏng xét đến liên kết bám dính dính nhớt Viscosity cho giá trị tải trọng phá hoại đạt 138.2 kN với độ võng giữa dầm 18.4 mm. Kết quả này chỉ chênh lệch 3.2% so với tải trọng thực nghiệm đo được là 133.9 kN, trong khi mô hình liên kết cứng tuyệt đối cho kết quả sai lệch vượt quá 18.6%.

Thứ ba, trên cấu kiện dầm bê tông xỉ thép và dầm bê tông thường sử dụng 5 thanh thép kéo đường kính 16 mm có diện tích 1005 mm2, mô hình tương tác bám dính không hoàn hảo đều thể hiện sự tiệm cận vượt bậc với độ tin cậy trên 94%. Sai số chuyển vị tại điểm uốn phá hoại giữa dầm giảm mạnh từ mức 22.4% ở mô hình dính tuyệt đối xuống chỉ còn khoảng 4.5% đến 6.1% ở mô hình xét đến độ nhớt dính bám.

Thảo luận kết quả

Nguyên nhân cốt lõi dẫn đến sự khác biệt giữa mô hình dính bám tuyệt đối và thực nghiệm là do hiện tượng trượt vi mô dọc theo chiều dài tiếp xúc giữa gân thép và ma trận bê tông. Khi tải trọng tác dụng tập trung ở giữa nhịp tăng dần, vùng bê tông chịu kéo phát sinh vết nứt thẳng đứng, làm suy giảm đột ngột diện tích tiếp xúc hiệu dụng và gây trượt cục bộ. Việc áp dụng mô hình liên kết dính kết Cohesive kết hợp hệ số ổn định độ nhớt Viscosity đã giải quyết thành công sự tập trung ứng suất giả tạo tại các nút phần tử, phản ánh chính xác cơ chế suy giảm dính bám thực tế.

So sánh với các nghiên cứu trước đây vốn chỉ áp dụng ràng buộc liên kết ngàm cứng toàn phần, phương pháp mô phỏng này cải thiện độ chính xác tổng thể thêm khoảng 12% đến 15%. Toàn bộ diễn biến làm việc của cấu kiện dầm có thể được trình bày một cách trực quan thông qua biểu đồ đường cong quan hệ tải trọng - chuyển vị giữa dầm và bảng tổng hợp đối chiếu sai số đa thông số. Biểu đồ uốn phản ánh rõ nét ba giai đoạn làm việc đặc trưng: giai đoạn đàn hồi ban đầu, giai đoạn phát triển vết nứt nứt dẻo và giai đoạn chảy dẻo phá hoại của cốt thép dọc, chứng minh tính ưu việt của mô hình trên cả ba nền vật liệu nghiên cứu.

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

Thứ nhất, chuẩn hóa quy trình thiết lập thông số tương tác tiếp xúc bề mặt Cohesive và hệ số độ nhớt Viscosity trong phần mềm ABAQUS nhằm kiểm soát sai số tính toán của dầm bê tông cốt thép dưới 5%. Các kỹ sư kết cấu tại các viện nghiên cứu và công ty tư vấn thiết kế cần chủ trì ban hành hướng dẫn thực hành chi tiết trong thời gian 6 tháng tới.

Thứ hai, mở rộng quy mô thí nghiệm kéo tuột cho các nhóm cốt thép có đường kính lớn từ 22 mm đến 32 mm trên nền bê tông xanh geopolymer và bê tông phụ gia xỉ thép. Việc này giúp nâng cao độ chính xác của các đường cong nội suy lên trên 98%, cần được các phòng thí nghiệm trọng điểm chuyên ngành xây dựng triển khai thực hiện trong vòng 12 tháng.

Thứ ba, ứng dụng mô hình số bám dính phi tuyến vào công tác kiểm toán thiết kế kết cấu dầm nhịp lớn trên 12 m và các cấu kiện chịu tải trọng nặng trong công trình công nghiệp. Giải pháp này giúp tối ưu hóa từ 10% đến 15% hàm lượng cốt thép sử dụng mà vẫn đảm bảo an toàn chịu lực, do các đơn vị thẩm tra và chủ đầu tư thực hiện liên tục qua các quý trong năm.

Thứ tư, tích hợp module hướng dẫn mô phỏng tương tác dính bám cốt thép - bê tông vào chương trình đào tạo sau đại học và các khóa chuyển giao công nghệ ứng dụng phần tử hữu hạn nâng cao. Mục tiêu đặt ra là giúp 100% học viên cao học làm chủ kỹ năng mô phỏng tiếp xúc phi tuyến trong thời gian khóa đào tạo 3 tháng, do các giảng viên chuyên ngành công trình đảm nhiệm.

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

Nhóm 1 là các kỹ sư thiết kế và chuyên gia phân tích kết cấu công trình. Nhóm đối tượng này có thể ứng dụng trực tiếp các thiết lập tham số tiếp xúc mặt Cohesive và phần tử C3D8R trong ABAQUS để tính toán độ võng, bề rộng vết nứt cho các kết cấu dầm phức tạp, giúp tăng độ chính xác của mô hình lên hơn 15% so với phương pháp tính toán giải tích đơn giản.

Nhóm 2 là các nhà khoa học và chuyên gia phát triển vật liệu xây dựng bền vững. Nghiên cứu cung cấp cơ sở dữ liệu thực nghiệm chuẩn xác về ứng xử cơ học của bê tông geopolymer không xi măng và bê tông xỉ thép, hỗ trợ quá trình nghiên cứu tối ưu hóa cấp phối và thương mại hóa các dòng vật liệu xanh thay thế cốt liệu tự nhiên trong thực tế.

Nhóm 3 là giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Kỹ thuật Xây dựng dân dụng và công nghiệp. Luận văn đóng vai trò là tài liệu tham khảo chuẩn mực về phương pháp kết hợp giữa thí nghiệm kéo tuột vật liệu, cơ học phá hoại dẻo và mô phỏng số 3D phục vụ công tác giảng dạy và thực hiện đề tài nghiên cứu chuyên sâu.

Nhóm 4 là các đơn vị tư vấn thẩm tra và quản lý chất lượng dự án. Tài liệu cung cấp góc nhìn chuyên sâu về cơ chế suy giảm lực dính bám và dự báo rủi ro trượt cốt thép trong thi công dầm khẩu độ lớn, giúp các kỹ sư giám sát kiểm soát chặt chẽ chất lượng đầm dùi và neo cốt thép, giảm thiểu rủi ro nứt nẻ công trình xuống dưới 2%.

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

Việc bỏ qua sự bám dính không hoàn hảo giữa bê tông và cốt thép gây ra ảnh hưởng gì trong tính toán mô phỏng? Khi chịu tải trọng uốn lớn, giữa bê tông và cốt thép luôn xuất hiện hiện tượng trượt vi mô. Nếu giả định dính bám hoàn hảo tuyệt đối, độ cứng của mô hình sẽ cao hơn thực tế khoảng 15% đến 20%, dẫn đến đánh giá sai lệch khả năng chịu tải và độ võng cực hạn của dầm khi xuất hiện vết nứt.

Mô hình vật liệu Hsu – Hsu (1994) có vai trò gì trong phân tích phần tử hữu hạn bê tông? Mô hình Hsu – Hsu mô tả chính xác toàn bộ đường cong ứng suất - biến dạng phi tuyến của bê tông chịu nén từ trạng thái đàn hồi đến khi phá hoại với cường độ lên đến 60 MPa. Khi kết hợp cùng mô hình độ cứng chịu kéo hiệu chỉnh, mô hình giúp phần mềm ABAQUS phân tích ổn định và triệt tiêu lỗi phân kỳ toán học.

Làm thế nào để đưa các giá trị từ thí nghiệm kéo tuột cốt thép vào mô phỏng dầm kết cấu? Thông qua thí nghiệm kéo tuột các thanh thép đường kính 12 mm, 16 mm và 20 mm, tác giả trích xuất lực kéo lớn nhất và năng lượng phá hủy diện tích bề mặt. Sau đó, phương pháp nội suy toán học được áp dụng để xác định các đặc trưng dính kết tương ứng cho cốt thép đường kính 14 mm rồi gán vào thuộc tính tương tác Cohesive Surface trong ABAQUS.

Khả năng chịu lực của dầm bê tông Geopolymer so với bê tông thường thể hiện ra sao? Bê tông Geopolymer có cường độ chịu nén cao trên 35 MPa và lực dính bám với cốt thép gân tương đương bê tông thường. Kết quả phân tích cho thấy dầm bê tông geopolymer đạt sức chịu tải cực hạn 138.2 kN, chứng minh khả năng thay thế hoàn toàn bê tông xi măng truyền thống trong các kết cấu chịu uốn.

Ưu điểm của việc kết hợp phần tử khối C3D8R và phần tử thanh T3D2 trong mô hình dầm là gì? Phần tử khối 8 nút C3D8R được dùng cho bê tông và thép dọc nhằm nắm bắt trạng thái ứng suất biến dạng 3 chiều phức tạp và tương tác bề mặt tiếp xúc. Cốt đai được gán phần tử thanh 2 nút T3D2 giúp đơn giản hóa mô hình và tiết kiệm hơn 30% thời gian tính toán mà vẫn đảm bảo độ chính xác.

Kết luận

  • Luận văn đã xây dựng thành công quy trình mô phỏng số 3D trên phần mềm ABAQUS phản ánh chân thực cơ chế bám dính không hoàn hảo giữa bê tông và cốt thép khi cấu kiện bị phá hoại.
  • Phân tích và kiểm chứng toàn diện ứng xử uốn dầm trên 3 hệ vật liệu: bê tông thường, bê tông xỉ thép và bê tông geopolymer thân thiện với môi trường.
  • Chứng minh việc ứng dụng mô hình tương tác bề mặt dính nhớt Viscosity giúp giảm sai số tải trọng phá hoại từ trên 18% xuống dưới 3.5% so với số liệu thực nghiệm.
  • Xác lập thành công bộ thông số dính kết và năng lượng phá hủy cho cốt thép đường kính 14 mm thông qua phương pháp nội suy từ chuỗi thí nghiệm kéo tuột tiêu chuẩn.
  • Đóng góp giải pháp khoa học mang tính thực tiễn cao, giúp tiết kiệm tới 40% chi phí thử nghiệm phá hủy và tối ưu hóa độ an toàn trong thiết kế kết cấu công trình.

Trong lộ trình 12 tháng tiếp theo, nghiên cứu cần được tiếp tục mở rộng trên các cấu kiện cột chịu nén uốn đồng thời và dầm liên tục nhiều nhịp. Hãy áp dụng ngay giải pháp mô phỏng bám dính cải tiến này vào các dự án phân tích kết cấu thực tế để nâng cao độ tin cậy và tối ưu hóa chi phí đầu tư cho công trình.