Tổng quan nghiên cứu

Sự gia tăng mật độ phương tiện giao thông tải trọng lớn đang đặt ra thách thức nghiêm trọng đối với độ an toàn của các công trình xây dựng. Trong thực tế, nhiều vụ tai nạn thảm khốc đã xảy ra khi phương tiện mất lái đâm trực diện vào kết cấu chịu lực. Điển hình vào ngày 06/05/2018 tại Cửa Lò, Nghệ An, một xe rác va chạm làm gãy cột trụ nhà dân khiến góc công trình sụp đổ và làm 3 người thương vong. Tương tự, vào ngày 03/10 trên Quốc lộ 1A thuộc tỉnh Quảng Ngãi, một xe tải chở gỗ mất lái đã tông sập 4 căn nhà liền kề. Những sự cố này minh chứng rằng tải trọng va đập đột ngột gây ra hiểm họa phá hủy dây chuyền vô cùng nặng nề.

Hiện nay, các tiêu chuẩn thiết kế như Eurocode 1 Phần 1-1 chỉ giới hạn vận tốc tính toán va chạm dưới 16.2 km/h, trong khi tiêu chuẩn Tây Ban Nha CTE quy định lực tĩnh tương đương cố định ở mức 50 kN. Những quy định này chưa phản ánh đầy đủ bản chất phi tuyến phức tạp khi xảy ra va chạm ở dải vận tốc cao từ 33 km/h đến 96.6 km/h.

Luận văn thạc sĩ chuyên ngành Kỹ thuật Xây dựng Công trình Dân dụng và Công nghiệp, thực hiện tại Trường Đại học Bách Khoa – ĐHQG TP.HCM vào tháng 08/2019, tập trung mô phỏng ứng xử phi tuyến của dầm và cột bê tông cốt thép chịu tải trọng va đập bằng phương pháp phần tử hữu hạn trên phần mềm LS-DYNA. Nghiên cứu phân tích các đại lượng động lực học quan trọng gồm ứng suất Von Mises, biến dạng dẻo hiệu quả, chuyển vị và phổ lực va đập trong khoảng thời gian từ 0.01 đến 0.2 giây. Kết quả nghiên cứu đóng góp cơ sở dữ liệu chuẩn xác, hỗ trợ kỹ sư tối ưu hóa giải pháp kháng chấn và gia cường, giảm thiểu tới hơn 40% nguy cơ sụp đổ công trình khi chịu va chạm bất khả kháng.

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 ứng dụng hệ phương trình vi phân cân bằng động lực học kết cấu, được giải bằng phương pháp tích phân trực tiếp Newmark với các hệ số ổn định gamma bằng 0.5 và beta bằng 0.25, kết hợp phương pháp sai phân trung tâm với bước thời gian tới hạn nghiêm ngặt. Phản ứng phi tuyến của vật liệu được xây dựng trên nền tảng cơ học phá hủy liên tục, trong đó hàm suy giảm mô đun đàn hồi Young được xác định qua tham số tổn thương d biến thiên từ 0 đến 1.

Dưới tác động của tải trọng động tốc độ cao, hiệu ứng tốc độ biến dạng đóng vai trò quyết định đến độ bền vật liệu. Khi tốc độ biến dạng gia tăng từ 10^-4 s^-1 lên đến ngưỡng 10^4 s^-1, hệ số tăng động DIF làm gia tăng đáng kể cường độ nén và kéo của bê tông. Quá trình này được thúc đẩy bởi lực quán tính ngăn cản sự giãn nở ngang, tạo nên trạng thái ứng suất nén ba trục cục bộ và buộc các vết nứt phải xuyên qua cốt liệu thay vì men theo pha xi măng. Tốc độ truyền sóng ứng suất trong bê tông đạt khoảng 3500 m/s, chi phối cơ chế phá hoại cục bộ trước khi kết cấu kịp phản ứng tổng thể. Nghiên cứu áp dụng mô hình bề mặt nắp liên tục CSCM cho bê tông và mô hình dẻo đẳng hướng tuyến tính từng đoạn Piecewise cho cốt thép nhằm mô phỏng chính xác sự tương tác đàn dẻo.

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

Luận văn thu thập nguồn dữ liệu thực nghiệm chuẩn quốc tế để xây dựng và hiệu chỉnh mô hình số. Cỡ mẫu nghiên cứu bao gồm 4 trường hợp mô phỏng kiểm chứng độc lập và 2 bài toán mô phỏng va chạm quy mô lớn. Các bộ dữ liệu được chọn mẫu có chủ đích dựa trên các thí nghiệm điển hình của Zhao và cộng sự năm 2017 về dầm bê tông cốt thép, thí nghiệm của Feyeabend về cột chịu tải trọng rơi, dữ liệu va chạm của Murray ở vận tốc 33 km/h và thử nghiệm va chạm thực tế của xe tải Ford F800 nặng 2.89 tấn ở vận tốc 96.6 km/h do Buth thực hiện năm 2009.

Phương pháp phân tích số được thực hiện qua bộ giải tường minh của phần mềm LS-DYNA. Bê tông được chia lưới bằng phần tử khối rắn 8 nút với 1 điểm tích phân Gauss, kiểm soát tỷ lệ năng lượng Hourglass nghiêm ngặt dưới 5% tổng nội năng. Cốt thép được mô hình hóa bằng phần tử dầm Hughes-Liu liên kết chặt chẽ với bê tông qua kỹ thuật chia sẻ nút hoặc ràng buộc tương hỗ. Trọng lượng bản thân kết cấu được thiết lập trạng thái ứng suất ban đầu thông qua phương pháp thư giãn động Dynamic Relaxation. Toàn bộ quy trình mô phỏng và phân tích dữ liệu được triển khai hoàn chỉnh trong giai đoạn từ năm 2018 đến năm 2019.

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

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

Phổ lực va đập theo thời gian ghi nhận quy luật xuất hiện 2 đỉnh xung lực rõ rệt. Đỉnh lực thứ nhất phát sinh trong khoảng 0.01 đến 0.02 giây đầu tiên khi mui xe tiếp xúc và biến dạng; đỉnh lực thứ hai có giá trị cực đại từ 1500 kN đến hơn 2200 kN xuất hiện ở thời điểm 0.03 đến 0.05 giây khi khối động cơ có độ cứng lớn đâm trực tiếp vào cột bê tông cốt thép.

Tác động của hiệu ứng tốc độ biến dạng trên 10^2 s^-1 làm tăng khả năng kháng nén động của bê tông từ 25% đến 45% so với giá trị cường độ tĩnh ban đầu. Biến dạng dẻo hiệu quả tập trung cao độ tại vùng chân cột trong phạm vi chiều cao từ 0.3 m đến 1.2 m. Chuyển vị ngang cực đại của cột dao động từ 35 mm đến 120 mm tùy thuộc vào vận tốc va chạm của phương tiện.

Cơ chế phá hoại của cấu kiện thể hiện rõ sự hình thành vết nứt trượt cắt chéo trước khi đạt tới giới hạn xoay dẻo uốn. Năng lượng tiêu tán của hệ thống được kiểm soát ổn định, trong đó tỷ lệ năng lượng Hourglass duy trì ở mức dưới 4.2%, đảm bảo độ tin cậy và tính hội tụ cao của mô hình số.

Thảo luận kết quả

Sự xuất hiện của phá hoại cắt cục bộ bắt nguồn từ vận tốc truyền sóng 3500 m/s. Khi phương tiện va chạm ở tốc độ cao, sóng ứng suất chưa kịp phân tán đều khắp chiều dài cấu kiện đã tạo ra ứng suất tiếp tuyến kéo vượt quá khả năng chịu lực của bê tông tại vùng tiếp xúc.

So sánh với các tiêu chuẩn thiết kế quốc tế, giá trị lực va đập tức thời từ mô phỏng FEM lớn hơn từ 35% đến 60% so với giá trị lực tĩnh tương đương quy định trong Eurocode 1 Phần 1-7 (thường dao động từ 500 kN đến 1000 kN). Tiêu chuẩn Trung Quốc JTG D81-2006 xem xét va chạm dưới dạng dao động điều hòa đơn giản đã bỏ qua xung lực thứ hai từ động cơ, dẫn đến đánh giá thiếu an toàn. Đối chiếu với tiêu chuẩn thiết kế kháng nổ ASCE, góc xoay gối tại chân cột vượt quá ngưỡng 2 độ, xếp cấu kiện vào mức độ hư hỏng trung bình đến nặng.

Để phục vụ công tác thiết kế thực tế, các dữ liệu này nên được thể hiện qua biểu đồ lịch sử thời gian của lực va chạm, biểu đồ phân bố ứng suất Von Mises và bảng tổng hợp chuyển vị ngang theo chiều cao cột nhằm giúp kỹ sư dễ dàng nhận diện vị trí xung yếu.

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

Cập nhật tiêu chuẩn thiết kế kết cấu: Đề xuất Bộ Xây dựng và các viện nghiên cứu bổ sung quy định tính toán tải trọng va đập động lực học phi tuyến cho các công trình nằm trong phạm vi 5.0 m tính từ mép đường giao thông. Mục tiêu giảm thiểu 40% sai số thiết kế so với phương pháp tĩnh truyền thống, áp dụng trong lộ trình quy chuẩn giai đoạn 2026 - 2028.

Tăng cường cốt thép kháng cắt cho cấu kiện: Yêu cầu các đơn vị tư vấn thiết kế thu hẹp khoảng cách cốt thép đai xuống mức 80 mm đến 100 mm trong đoạn 1.5 m tính từ chân cột tầng trệt. Kết hợp giải pháp bọc áo thép dày 6 mm đến 10 mm hoặc bọc sợi composite carbon CFRP để tăng khả năng tiêu tán năng lượng va đập thêm 35% đến 50%, thực hiện ngay từ năm 2026.

Lắp đặt dải phân cách và trụ bảo vệ thụ động: Các ban quản lý dự án hạ tầng đô thị cần triển khai hệ thống rào chắn bê tông gờ chữ F cao từ 0.81 m đến 1.1 m hoặc lắp đặt trụ thép nhồi bê tông trước các cột trụ cầu và tòa nhà trọng yếu. Giải pháp này giúp triệt tiêu từ 50% đến 70% động năng ban đầu của phương tiện có khối lượng từ 3.0 đến 8.0 tấn di chuyển ở vận tốc 60 km/h đến 80 km/h trước khi va chạm vào kết cấu chính.

Chuẩn hóa quy trình mô phỏng số chuyên ngành: Các viện nghiên cứu và doanh nghiệp tư vấn xây dựng cần áp dụng cẩm nang mô phỏng số chuẩn hóa trên phần mềm LS-DYNA, sử dụng mô hình vật liệu CSCM và giới hạn năng lượng Hourglass dưới 5% để đạt độ chính xác dự báo phản ứng kết cấu trên 90%, hoàn thành ban hành nội bộ trong năm 2026.

Đố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 dân dụng và công nghiệp: Ứng dụng kết quả nghiên cứu để tính toán cốt thép đai chịu cắt, kiểm tra khả năng xoay dẻo của cột tầng hầm và tầng trệt khi chịu tải trọng xe tải va chạm từ 2.5 đến 10 tấn, ngăn ngừa nguy cơ sụp đổ dây chuyền.

Kỹ sư quy hoạch và thiết kế hạ tầng giao thông: Khai thác mô hình tương tác va chạm giữa xe và kết cấu để thiết kế các kết cấu mố trụ cầu vượt, hệ thống rào chắn phân cách trên đường cao tốc chịu được vận tốc va chạm trên 80 km/h.

Giảng viên, nghiên cứu sinh và học viên cao học xây dựng: Sử dụng luận văn làm tài liệu chuyên khảo sâu về cơ học động lực học phi tuyến, kỹ thuật chia lưới phần tử hữu hạn explicit và phương pháp gán mô hình vật liệu CSCM trong LS-DYNA.

Đơn vị thẩm tra kỹ thuật và công ty bảo hiểm công trình: Sử dụng bộ tiêu chí biến dạng dẻo và góc xoay ASCE trong luận văn làm căn cứ khoa học để thẩm định hồ sơ an toàn kết cấu và giám định tổn thất sau tai nạn đâm va phương tiện.

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

Tại sao phương pháp lực tĩnh tương đương trong tiêu chuẩn không đủ an toàn khi tính toán va đập? Phương pháp lực tĩnh tương đương chỉ áp dụng mức tải cố định từ 50 kN đến 1000 kN, bỏ qua hiệu ứng tốc độ biến dạng và sự lan truyền sóng ứng suất 3500 m/s. Thực tế mô phỏng cho thấy xe tải 2.89 tấn chạy 96.6 km/h tạo ra xung lực đỉnh vượt 1800 kN trong 0.02 giây, gây phá hoại cắt tức thời mà tính toán tĩnh không thể bao quát hết.

Mô hình vật liệu CSCM trong LS-DYNA mang lại lợi thế gì khi mô phỏng bê tông? Mô hình CSCM tích hợp mặt chảy dẻo liên tục kết hợp nắp mặt cầu, cho phép mô tả chính xác phản ứng nén ba trục, hóa mềm khi kéo và sự suy giảm mô đun đàn hồi với tham số tổn thương d từ 0 đến 1. Mô hình tự động điều chỉnh hệ số tăng động DIF khi tốc độ biến dạng thay đổi từ 10^-4 s^-1 đến 10^4 s^-1.

Năng lượng Hourglass là gì và làm thế nào để đảm bảo mô phỏng FEM không bị sai lệch? Năng lượng Hourglass phát sinh do biến dạng giả không sinh ứng suất khi sử dụng phần tử khối rắn 8 nút tích phân 1 điểm Gauss. Để đảm bảo kết quả mô phỏng đáng tin cậy, kỹ sư cần thiết lập lực kháng Hourglass phù hợp và kiểm tra biểu đồ năng lượng sao cho tỷ lệ năng lượng Hourglass luôn duy trì dưới 5% tổng nội năng hệ thống.

Tốc độ biến dạng cao tác động như thế nào đến ứng xử phá hoại của cột bê tông cốt thép? Ở tốc độ biến dạng cao trên 10^2 s^-1, lực quán tính cản trở vật liệu giãn nở hông, tạo hiệu ứng bó hông ba trục làm tăng cường độ nén động của bê tông thêm 25% đến 45%. Đồng thời, vết nứt buộc phải truyền xuyên qua cốt liệu, làm gia tăng tính giòn và khiến cấu kiện dễ bị phá hoại dạng trượt cắt cục bộ.

Biện pháp nào hiệu quả nhất để gia cường kháng va đập cho các cột bê tông cốt thép hiện hữu? Giải pháp tối ưu là bọc áo thép dày 6 mm đến 10 mm hoặc dán 2 đến 3 lớp sợi composite carbon CFRP quanh vùng chân cột cao 1.5 m kết hợp chèn vữa không co ngót. Biện pháp này giúp tăng khả năng kháng cắt lên trên 40%, hạn chế góc xoay gối dưới 1.5 độ nhằm đáp ứng cấp phá hoại nhẹ theo tiêu chuẩn ASCE.

Kết luận

  • Xây dựng thành công quy trình mô phỏng số phi tuyến explicit bằng phần mềm LS-DYNA nhằm khảo sát chi tiết phản ứng động học của dầm và cột bê tông cốt thép dưới tác động va đập của phương tiện giao thông.
  • Xác định phổ lực va chạm thực tế gồm 2 đỉnh xung lực xuất hiện trong thời gian cực ngắn dưới 0.05 giây, với xung lực đỉnh vượt tải trọng tĩnh quy ước từ 35% đến 60%.
  • Chứng minh cơ chế phá hoại trượt cắt cục bộ chiếm ưu thế tại vùng chân cột trong phạm vi từ 0.3 m đến 1.2 m do ảnh hưởng của sóng ứng suất truyền với tốc độ 3500 m/s.
  • Làm rõ vai trò của hệ số tăng động DIF và hiệu ứng bó hông quán tính giúp gia tăng cường độ nén tức thời của bê tông lên 1.25 đến 1.45 lần ở tốc độ biến dạng cao.
  • Đề xuất các giải pháp cấu tạo cốt thép đai dày với khoảng cách 80 mm đến 100 mm kết hợp rào chắn tiêu năng giúp giảm trên 50% nguy cơ sụp đổ công trình.

Luận văn đã đóng góp cơ sở dữ liệu số hóa tin cậy và phương pháp luận khoa học vững chắc cho chuyên ngành Kỹ thuật Xây dựng tại Việt Nam trong bài toán phòng chống tải trọng đặc biệt. Trong giai đoạn 2026 - 2028, hướng nghiên cứu tiếp theo cần mở rộng khảo sát cho các kịch bản va chạm góc nghiêng từ 30 đến 60 độ và tổ hợp đồng thời tải trọng nổ - va đập. Các kỹ sư và nhà quản lý công trình cần chủ động ứng dụng phương pháp mô phỏng explicit và tiêu chuẩn biến dạng ASCE vào thực tiễn thiết kế để bảo vệ tối đa tính mạng con người và sự an toàn bền vững của công trình.