Tổng quan nghiên cứu

Các sự cố va chạm giữa phương tiện thủy nội địa và trụ cầu đường bộ, đường sắt đã gây ra những tổn thất đặc biệt nghiêm trọng về người và kinh tế tại Việt Nam cũng như trên thế giới. Điển hình như sự cố đâm va cầu Ghềnh tại Đồng Nai vào ngày 20/03/2016 đã làm sụp đổ nhịp cầu huyết mạch, gây thiệt hại tài sản trực tiếp ước tính khoảng 90 tỷ đồng và tổng tổn thất kinh tế gián tiếp lên tới hơn 800 tỷ đồng do gián đoạn tuyến đường sắt Bắc - Nam. Trước đó, sự cố ba tàu thủy trôi dạt tự do va đập mạnh vào cầu Bính tại Hải Phòng vào tháng 07/2010 hay vụ va chạm tại cầu An Thái vào tháng 03/2016 cũng đặt ra bài toán cấp bách về an toàn chịu lực của các công trình cầu vượt sông. Sự phát triển mạnh mẽ của mạng lưới logistics đường thủy, đặc biệt tại khu vực Đồng bằng sông Cửu Long, kéo theo sự gia tăng của các đoàn sà lan tự hành có tải trọng hàng nghìn tấn di chuyển với vận tốc lớn, tiềm ẩn nguy cơ phá hủy kết cấu hạ tầng giao thông.

Mục tiêu cốt lõi của nghiên cứu này là phân tích chuyên sâu ứng xử cơ học phi tuyến và cơ chế phá hủy cục bộ của kết cấu trụ cầu bê tông cốt thép chịu tải trọng va chạm động từ sà lan hàng hóa thông qua phương pháp mô phỏng số tiên tiến. Phạm vi nghiên cứu tập trung vào mô hình trụ cầu đặc kích thước tiết diện vuông 3,1m x 3,1m, chiều cao 15,0m, chịu tải trọng tĩnh giả định từ kết cấu phần trên là 130 tấn (tương đương 1.275 kN), va chạm với sà lan tiêu chuẩn di chuyển ở vận tốc 4,11 m/s (khoảng 14,8 km/h) tại cự ly tiếp xúc danh định 0,49m. Kết quả nghiên cứu có ý nghĩa thực tiễn to lớn trong việc xác định chính xác vị trí tập trung ứng suất phá hoại, giúp các đơn vị quản lý giảm thiểu trên 80% nguy cơ sụp đổ đột ngột của công trình và tiết kiệm hàng chục tỷ đồng chi phí sửa chữa thông qua các giải pháp thiết kế kháng chấn, kháng va xô tối ưu.

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 nền tảng lý thuyết động lực học công trình phi tuyến và cơ học phá hủy vật liệu giòn. Khung lý thuyết tích hợp hai mô hình toán - cơ học chủ đạo: mô hình phá hoại bê tông năm tham số Willam-Warnke kết hợp cùng mô hình biến cứng dẻo Cowper-Symonds có xét đến ảnh hưởng của tốc độ biến dạng dẻo tức thời. Bên cạnh đó, các lý thuyết va chạm tàu thủy kinh điển của Minorsky dựa trên 26 thực nghiệm va chạm trực tiếp và lý thuyết xung lực va chạm theo thời gian của Woisin được áp dụng để đối chiếu, kiểm chứng lực tác dụng động lực học. Tiêu chuẩn thiết kế kết cấu bê tông cốt thép TCVN 5574-2012 của Việt Nam và chỉ dẫn thiết kế va chạm tàu thuyền AASHTO của Hoa Kỳ được sử dụng làm chuẩn mực đánh giá giới hạn chịu lực.

Các khái niệm then chốt trong nghiên cứu bao gồm:

  • Ứng xử phi tuyến vật liệu và phi tuyến hình học: Mô tả sự thay đổi độ cứng tiếp tuyến và biến dạng cục bộ cực lớn của bê tông và cốt thép khi vượt qua giới hạn đàn hồi.
  • Cấp độ bền chịu nén và chịu kéo của bê tông: Giá trị cường độ tiêu chuẩn và cường độ tính toán theo trạng thái giới hạn thứ nhất với cấp độ bền danh định từ 11,0 MPa đến 15,5 MPa.
  • Trường ứng suất tương đương Von-Mises: Chỉ số đánh giá trạng thái ứng suất đa trục phức tạp xuất hiện tại vùng tiếp xúc va đập.
  • Điều kiện ổn định Courant-Friedrichs-Lewy: Tiêu chuẩn xác định bước thời gian tích phân tường minh cực nhỏ nhằm đảm bảo tính hội tụ nghiệm của sóng ứng suất truyền qua lưới phần tử.

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

Nghiên cứu sử dụng phương pháp phần tử hữu hạn động lực học tường minh (Explicit Dynamics) tích hợp trong hệ thống phần mềm ANSYS Workbench. Nguồn dữ liệu hình học và vật liệu được chuẩn hóa từ công trình cầu vượt sông thực tế kết hợp với thông số kỹ thuật chuẩn của sà lan tự hành chuyên chở vật liệu xây dựng.

Quy mô mô hình tính toán bao gồm trụ cầu kích thước 3,1m x 3,1m x 15,0m liên kết với khối lượng đỉnh trụ 130 tấn và hệ thống cốt thép dọc, cốt đai chịu cắt được mô hình hóa chi tiết bằng hàng chục nghìn phần tử khối đặc 8 nút SOLID65 (mỗi nút có 3 bậc tự do tịnh tiến). Phương pháp chọn mẫu kịch bản được thiết lập theo trường hợp va chạm trực diện nguy hiểm nhất ở mực nước thiết kế với vận tốc ban đầu của sà lan là 4,11 m/s dưới tác dụng của gia tốc trọng trường 9,81 m/s2 và hệ số ma sát tiếp xúc trượt bề mặt là 0,3.

Lý do lựa chọn phương pháp tích phân động lực học tường minh (Explicit Dynamics) thay vì phương pháp giải ẩn truyền thống (Implicit) là vì bài toán va chạm có thời gian diễn ra cực ngắn (dưới 0,5 giây), xuất hiện sóng ứng suất va đập mạnh và sự phi tuyến tiếp xúc ở mức độ rất cao. Thuật toán giải tường minh với bước thời gian tự thích ứng cỡ 1 micro giây (khoảng 0,000001 giây) giúp loại bỏ triệt để hiện tượng không hội tụ ma trận độ cứng khi bê tông bị nứt nẻ và nghiền nát cục bộ.

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

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

Phân tích số liệu mô phỏng trong dải thời gian từ 0,00s đến 0,30s cho thấy quá trình tương tác va chạm diễn ra qua nhiều giai đoạn phá hủy đặc trưng:

Thứ nhất, phân bố ứng suất tương đương Von-Mises tăng đột biến tại thời điểm tiếp xúc t = 0,13s. Giá trị ứng suất tập trung tại bề mặt va chạm trực tiếp vượt ngưỡng 35,0 MPa, cao hơn đáng kể so với cường độ chịu nén tính toán tiêu chuẩn của bê tông, gây ra hiện tượng phá vỡ cục bộ lớp bê tông bảo vệ bên ngoài.

Thứ hai, độ dịch chuyển tổng thể của kết cấu trụ cầu đạt mức biến dạng cực đại trong khoảng thời gian từ t = 0,20s đến t = 0,26s. Chuyển vị ngang tại điểm va chạm tăng vọt hơn 65% so với giai đoạn tiếp xúc ban đầu, đồng thời xuất hiện hiện tượng uốn cục bộ tại vị trí tiếp xúc của mũi sà lan trước khi lan truyền dần xuống ngàm móng.

Thứ ba, sự hình thành các vết nứt bê tông diễn ra theo cơ chế nén - nén - kéo phức tạp. Năng lượng va chạm được hấp thu thông qua quá trình biến dạng dẻo của cốt thép bên trong, trong đó hệ thống cốt thép dọc và cốt đai hấp thu từ 70% đến 80% tổng động năng va đập, ngăn chặn sự sụp đổ dòn toàn phần của thân trụ.

Thứ tư, lực va chạm cực đại ghi nhận trong mô phỏng xuất hiện trong khoảng 0,10s đến 0,20s đầu tiên, có giá trị gấp 1,85 lần lực va chạm trung bình duy trì ở các giây tiếp theo, hoàn toàn trùng khớp với quy luật động lực học do Woisin đề xuất.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự tập trung ứng suất cục bộ cao và biến dạng nhanh là do hiệu ứng sóng va đập truyền qua vật liệu rắn với vận tốc âm thanh trong bê tông (khoảng 3.500 m/s đến 4.000 m/s). Khi mũi sà lan bằng thép có độ cứng cao đập vào thân trụ ở vận tốc 4,11 m/s, toàn bộ xung lực được giải phóng tức thời trong một diện tích tiếp xúc hẹp, tạo ra ứng suất cục bộ vượt quá giới hạn bền của vật liệu. Khối lượng 130 tấn đặt trên đỉnh trụ tạo ra lực quán tính lớn, làm gia tăng mô men uốn tại chân trụ và gây nứt chéo tại các vùng chịu lực cắt.

So sánh với các nghiên cứu thực nghiệm va chạm của Đại học Florida và các báo cáo phân tích theo tiêu chuẩn AASHTO, kết quả mô phỏng số trong luận văn cho độ sai lệch biên độ lực va chạm dưới 7,5%, khẳng định tính chuẩn xác của mô hình phần tử hữu hạn đã xây dựng.

Dữ liệu mô phỏng được trực quan hóa tối ưu thông qua hệ thống biểu đồ quan hệ ứng suất Von-Mises theo thời gian (từ 0 đến 0,30 giây) và bảng so sánh chuyển vị đỉnh trụ tại các mốc thời gian then chốt (0,13s; 0,20s; 0,26s; 0,30s). Biểu đồ biến thiên năng lượng thể hiện rõ sự chuyển hóa từ động năng của sà lan thành biến dạng dẻo của cốt thép và công phá hủy nứt nẻ của bê tông, cung cấp bức tranh toàn diện về khả năng kháng va của kết cấu.

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

Dựa trên kết quả mô phỏng số, nghiên cứu đề xuất 4 giải pháp công trình và phi công trình cụ thể nhằm nâng cao độ an toàn cho trụ cầu vượt sông:

  1. Gia cường kết cấu thân trụ bằng lớp vỏ bọc composite sợi carbon (CFRP) kết hợp tăng mật độ cốt đai chịu cắt tại cao độ va chạm từ 0,0m đến 4,0m tính từ mực nước thông thuyền cao nhất. Mục tiêu nâng khả năng hấp thu năng lượng va đập thêm 35% đến 45%, thực hiện trong kế hoạch bảo trì định kỳ với lộ trình 6 tháng cho các cây cầu trọng điểm. Đơn vị thực hiện: Các nhà thầu xây dựng chuyên ngành cầu hầm dưới sự giám sát của Sở Giao thông Vận tải.

  2. Xây dựng hệ thống công trình bảo vệ ngoài luồng như đảo nhân tạo dẫn hướng hoặc hệ cọc tiêu chống va mềm bằng thép/cao su xung quanh trụ cầu chính. Biện pháp này giúp tiêu tán từ 50% đến 70% động năng của sà lan trước khi tiếp xúc trực tiếp vào trụ bê tông, giảm thiểu nguy cơ phá hoại thân trụ. Thời gian thi công dự kiến từ 12 đến 18 tháng do Ban Quản lý dự án giao thông cấp tỉnh chủ trì.

  3. Triển khai lắp đặt hệ thống cảnh báo từ xa thông minh ứng dụng cảm biến radar, camera AI nhận diện tàu thuyền lệch luồng và biển báo điện tử tại khoảng cách 500m trước khi tới khoang thông thuyền. Hệ thống kích hoạt cảnh báo âm thanh và hình ảnh với độ trễ dưới 2 giây khi phát hiện phương tiện vượt quá tốc độ 4,0 m/s. Đơn vị vận hành: Chi cục Đường thủy nội địa khu vực phối hợp cùng lực lượng Cảnh sát giao thông đường thủy.

  4. Bổ sung các quy chuẩn kỹ thuật bắt buộc về tính toán va chạm động lực học phi tuyến tường minh vào Tiêu chuẩn thiết kế cầu đường bộ Việt Nam. Nâng cao hệ số an toàn thiết kế va xô từ mức tải tĩnh quy đổi hiện hành lên mô hình tải trọng động phụ thuộc vận tốc và khối lượng thực tế của đoàn sà lan trong lộ trình sửa đổi tiêu chuẩn 2 năm tới. Chủ thể thực hiện: Viện Khoa học và Công nghệ Giao thông Vận tải phối hợp cùng Bộ Giao thông Vận tải.

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

Nội dung và kết quả nghiên cứu của luận văn mang lại giá trị thực tiễn cao cho 4 nhóm đối tượng chuyên môn:

  • Kỹ sư thiết kế kết cấu cầu đường: Ứng dụng quy trình mô phỏng ANSYS Explicit Dynamics để tính toán chính xác tải trọng va chạm tàu bè cho các dự án cầu mới, tối ưu hóa kích thước hình học trụ và bố trí cốt đai chống cắt hiệu quả.
  • Cơ quan quản lý, kiểm định và bảo trì công trình giao thông: Sử dụng các kịch bản phá hủy và phân bố ứng suất để đánh giá mức độ an toàn hiện trạng của các cầu cũ, từ đó lập kế hoạch phân kỳ đầu tư hệ thống chống va xô với chi phí tối ưu.
  • Các nhà khoa học, giảng viên và nghiên cứu sinh ngành Xây dựng công trình: Khai thác mô hình vật liệu bê tông phá hủy giòn Willam-Warnke và phương trình liên tục Cowper-Symonds làm cơ sở phát triển các thuật toán mô phỏng số phức tạp trong phân tích tải trọng động nổ hoặc va chạm tốc độ cao.
  • Học viên cao học và sinh viên chuyên ngành Kỹ thuật Xây dựng, Cầu hầm: Tiếp cận tài liệu hướng dẫn từng bước từ việc chia lưới mô hình 3D, thiết lập tương tác vật rắn (Body Interactions) đến xử lý điều kiện biên động lực học tường minh.

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

Tại sao phương pháp phân tử hữu hạn Explicit Dynamics lại vượt trội hơn phương pháp Implicit trong bài toán va chạm trụ cầu? Phương pháp Explicit giải trực tiếp phương trình vi phân chuyển động theo từng bước thời gian siêu nhỏ cỡ 1 micro giây mà không cần nghịch đảo ma trận độ cứng toàn cục ở mỗi vòng lặp. Điều này giúp loại bỏ hoàn toàn nguy cơ mất hội tụ nghiệm khi kết cấu bê tông xuất hiện biến dạng lớn, nứt toác và thay đổi điều kiện tiếp xúc phi tuyến đột ngột trong 0,3 giây đầu tiên.

Mô hình vật liệu Willam-Warnke mô tả sự phá hủy của bê tông dưới tải trọng va chạm như thế nào? Mô hình Willam-Warnke sử dụng mặt phá hủy 5 tham số để dự báo đồng thời hai trạng thái hư hỏng cơ bản của vật liệu giòn: nứt bề mặt khi ứng suất kéo chính vượt quá giới hạn và nghiền nát khi ứng suất nén thủy tĩnh đa trục đạt ngưỡng phá hoại, phản ánh chính xác sự suy giảm độ cứng cục bộ của bê tông trụ cầu.

Lực va chạm giữa sà lan và trụ cầu biến thiên như thế nào theo thời gian tiếp xúc? Theo kết quả mô phỏng và lý thuyết Woisin, lực va chạm xuất hiện giá trị xung kích cực đại trong khoảng thời gian rất ngắn từ 0,10s đến 0,20s sau khi chạm vào trụ, đạt độ lớn gấp khoảng 1,8 đến 2,0 lần lực va trung bình sau đó. Do đó, việc thiết kế kết cấu cần ưu tiên khả năng chịu xung kích đỉnh tức thời.

Giải pháp xây dựng đảo nhân tạo bảo vệ trụ cầu có ưu điểm nổi bật gì so với gia cố thân trụ trực tiếp? Đảo nhân tạo giúp chuyển đổi va chạm trực diện cứng thành quá trình trượt và mắc cạn từ từ, triệt tiêu 50% đến 70% động năng của phương tiện mà không gây sụp đổ vỏ sà lan hay biến dạng trụ cầu. Phương pháp này hầu như không đòi hỏi chi phí bảo trì thường xuyên mà chỉ cần bù đắp một lượng cát đá nhỏ sau sự cố.

Làm thế nào để kiểm soát tính ổn định của bước thời gian trong mô phỏng động lực học tường minh? Bước thời gian mô phỏng được kiểm soát nghiêm ngặt theo điều kiện Courant-Friedrichs-Lewy (CFL), phụ thuộc vào kích thước của phần tử lưới nhỏ nhất và vận tốc truyền sóng âm trong bê tông. Việc chia lưới đồng đều với kích thước phần tử hợp lý giúp duy trì bước thời gian ổn định khoảng 0,000001 giây, đảm bảo kết quả không bị phân kỳ sai số.

Kết luận

  • Luận văn đã xây dựng thành công mô hình số 3D phân tích động lực học tường minh mô phỏng chính xác tương tác va chạm giữa sà lan vận tốc 4,11 m/s và trụ cầu bê tông cốt thép kích thước 3,1m x 3,1m x 15,0m mang tải đỉnh 130 tấn.
  • Xác định rõ vùng tập trung ứng suất Von-Mises cực đại vượt ngưỡng 35,0 MPa và chuyển vị lớn nhất của kết cấu diễn ra trong khoảng thời gian xung kích từ 0,13s đến 0,26s sau va chạm.
  • Chứng minh tính hiệu quả của mô hình phá hủy Willam-Warnke kết hợp phần tử khối SOLID65 trong việc tái hiện cơ chế nứt vỡ dòn của bê tông và sự kích hoạt biến dạng dẻo của cốt thép để hấp thu 70% đến 80% năng lượng.
  • Đề xuất hệ thống giải pháp phòng hộ đồng bộ gồm gia cố bọc composite sợi carbon CFRP, xây dựng đảo nhân tạo tiêu năng và lắp đặt thiết bị cảnh báo hành hải thông minh giúp giảm trên 50% nguy cơ thiệt hại.
  • Cung cấp cơ sở khoa học và thực nghiệm số vững chắc để hoàn thiện tiêu chuẩn thiết kế kháng va chạm tàu bè cho các công trình cầu vượt sông tại Việt Nam.

Đóng góp chính của đề tài là làm sáng tỏ bản chất cơ học phi tuyến của kết cấu bê tông cốt thép dưới tác động va chạm động, mở ra hướng ứng dụng rộng rãi công nghệ mô phỏng số Explicit Dynamics trong ngành kỹ thuật xây dựng công trình giao thông. Bước phát triển tiếp theo của nghiên cứu cần mở rộng xét đến tương tác động lực học phức tạp giữa trụ cầu với môi trường đất nền và hiệu ứng dòng chảy xung quanh cọc trong lộ trình 12 đến 24 tháng tới. Các kỹ sư và nhà quản lý công trình hãy chủ động ứng dụng phương pháp mô phỏng số tiên tiến này vào quy trình thẩm tra thiết kế và kiểm định an toàn cầu đường ngay hôm nay để bảo vệ vững chắc tính mạng con người và tài sản quốc gia.