Tổng quan nghiên cứu

Trong bối cảnh phát triển mạnh mẽ của hạ tầng giao thông hiện đại, các công trình cầu vượt nhịp lớn và cầu cạn ngày càng ưu tiên sử dụng vật liệu cường độ cao. Xu hướng này giúp tối ưu hóa tiết diện, giảm trọng lượng bản thân nhưng đồng thời làm kết cấu trở nên thanh mảnh và nhạy cảm hơn trước các tác động động học. Khi các phương tiện giao thông di chuyển với tải trọng trục từ 30 đến 50 tấn và vận tốc lưu thông đạt 80 đến 120 km/h (tương đương 22,2 đến 33,3 m/s), tương tác động lực học giữa xe và cầu trở thành bài toán kỹ thuật trọng yếu.

Đặc biệt, hiện tượng phanh gấp hoặc hãm xe trên cầu tạo ra các lực quán tính và lực ma sát cục bộ rất lớn theo phương dọc, tác động trực tiếp lên hệ dầm liên tục. Trong khi phần lớn các nghiên cứu truyền thống chỉ khảo sát tải trọng di động với vận tốc không đổi hoặc áp dụng trên kết cấu dầm giản đơn một nhịp, động lực học của dầm liên tục nhiều nhịp chịu tải hãm vẫn là một khoảng trống học thuật lớn. Sự xuất hiện của lực hãm làm thay đổi trạng thái ứng suất, làm gia tăng chuyển vị đứng và mô men uốn nguy hiểm tại cả vùng giữa nhịp lẫn các gối tựa trung gian.

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 của tác giả Lê Hoài Thương, dưới sự hướng dẫn của PGS. Nguyễn Trọng Phước tại Trường Đại học Bách Khoa - Đại học Quốc gia TP.HCM (hoàn thành trong giai đoạn 2017-2018), đã tập trung giải quyết bài toán này. Mục tiêu nghiên cứu là xây dựng mô hình toán học và thuật toán số hoàn chỉnh để phân tích phản ứng động của dầm liên tục nhiều nhịp dưới tác dụng của tải trọng di động có xét đến lực hãm theo 2 pha thời gian: pha chuyển động đều và pha giảm tốc dừng hẳn. Kết quả nghiên cứu cung cấp hệ số động lực (DAF) chính xác, cho thấy lực hãm có thể làm tăng biên độ dao động và nội lực của dầm lên 15% đến 35%, đóng góp cơ sở dữ liệu quan trọng cho việc hoàn thiện tiêu chuẩn thiết kế cầu đường hiện đại.

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 cơ học kết cấu phi tuyến và động lực học công trình hiện đại, tích hợp các mô hình lý thuyết cốt lõi sau:

Lý thuyết dầm uốn Euler-Bernoulli: Kết cấu dầm liên tục được mô hình hóa dựa trên giả thiết tiết diện phẳng trước và sau khi biến dạng. Độ võng và góc xoay của trục dầm được mô tả liên tục dọc theo chiều dài nhịp $L$, bỏ qua biến dạng trượt cắt đối với các dầm thanh mảnh có tỷ lệ chiều cao trên nhịp nhỏ hơn 1/10.

Hàm nội suy Hermite bậc ba: Mỗi phần tử dầm gồm 2 nút với 4 bậc tự do (chuyển vị đứng và góc xoay tại mỗi nút). Trường chuyển vị trong phần tử được xấp xỉ chính xác thông qua hệ 4 hàm dạng Hermite bậc ba, đảm bảo tính liên tục của độ võng và độ dốc dọc trục phần tử dầm.

Mô hình tải trọng xe dao động 2 bậc tự do (Sprung-mass model): Thay vì giả định tải trọng là lực tập trung tĩnh hoặc khối lượng di động đơn giản, nghiên cứu áp dụng mô hình gồm khối lượng thân xe $m_2$ và khối lượng bánh xe $m_1$, liên kết với nhau qua hệ lò xo độ cứng $k_s$ và cản nhớt $c_s$. Mô hình này phản ánh trung thực chuyển vị đứng độc lập của thân xe và bánh xe khi tương tác với mặt cầu.

Cơ chế ma sát và lực hãm 2 pha: Lực hãm tác dụng lên dầm bắt nguồn từ lực ma sát tiếp xúc giữa lốp xe và mặt đường với hệ số ma sát tĩnh $\mu_s = 0,6$ và hệ số ma sát động $\mu_k$. Tải trọng hãm được quy đổi về mô men uốn tập trung lệch trục $M = f_i \cdot h / 2$, với $h$ là chiều cao tiết diện dầm. Các thành phần gia tốc hướng tâm, gia tốc Coriolis và gia tốc do biến thiên vận tốc đều được thiết lập đầy đủ trong phương trình cân bằng động.

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 (FEM) kết hợp với thuật toán tích phân số trực tiếp theo miền thời gian:

Rời rạc hóa kết cấu và chọn mẫu phần tử: Dầm liên tục nhiều nhịp được chia thành lưới phần tử hữu hạn đồng nhất với mật độ từ 20 đến 40 phần tử trên mỗi nhịp để đảm bảo độ hội tụ tối ưu của ma trận khối lượng $M$, ma trận cản Rayleigh $C$ và ma trận độ cứng $K$. Việc lựa chọn lưới phần tử mịn giúp mô tả chính xác vị trí tiếp xúc tức thời của bánh xe tại từng bước thời gian.

Thiết lập phương trình chuyển động tổng thể: Dựa trên nguyên lý cân bằng động D'Alembert và hàm tiêu tán Rayleigh, hệ phương trình vi phân chuyển động phi tuyến tương tác xe - cầu được thiết lập dưới dạng ma trận phụ thuộc thời gian. Các thành phần khối lượng, độ cứng và cản bổ sung từ xe chỉ được tích hợp cục bộ vào phần tử dầm đang trực tiếp chịu tải.

Thuật toán tích phân từng bước Newmark: Phương trình vi phân cấp hai được giải bằng phương pháp Newmark không điều kiện ổn định với các tham số chuẩn $\gamma = 0,5$ và $\beta = 0,25$. Bước thời gian tính lặp $\Delta t$ được khảo sát nghiêm ngặt từ 0,001 giây đến 0,005 giây nhằm triệt tiêu hiện tượng khuếch đại sai số số học và nắm bắt trọn vẹn các thành phần tần số dao động bậc cao.

Quy trình lập trình và kiểm chứng: Toàn bộ thuật toán được tự động hóa bằng ngôn ngữ lập trình MATLAB. Độ tin cậy của chương trình được kiểm chứng chéo qua việc đối soát tần số dao động tự nhiên với phần mềm SAP2000 (sai số dưới 0,5%) và so sánh kết quả phản ứng động với các công bố quốc tế chuẩn mực.

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 nhiều kịch bản vận tốc, gia tốc hãm và vị trí bắt đầu phanh đã chỉ ra những quy luật động lực học quan trọng:

Gia tăng đáng kể hệ số động lực của chuyển vị (DAF_d): Khi xe thực hiện phanh hãm trên dầm liên tục, hệ số động lực chuyển vị tại điểm giữa các nhịp tăng từ 18% đến 32% so với trường hợp xe chạy qua với vận tốc không đổi. Tại vận tốc ban đầu $v = 35,57\text{ m/s}$, chuyển vị lớn nhất xuất hiện ngay trong pha hãm do sự cộng hưởng giữa dao động nhấp nhô của thân xe và chuyển vị uốn của dầm.

Gia tăng cực đại của mô men uốn cục bộ (DAF_m): Sự xuất hiện của mô men hãm lệch tâm $M = f_i \cdot h / 2$ làm mô men uốn tại mặt cắt gối tựa và giữa nhịp tăng vọt từ 20% đến 38%. Đặc biệt, khi gia tốc hãm đạt giá trị $a_b = -5,0\text{ m/s}^2$ đến $-7,5\text{ m/s}^2$, mô men âm tại các gối trung gian chịu sự biến thiên đột ngột, làm tăng nguy cơ nứt bê tông vùng chịu kéo.

Ảnh hưởng của vị trí bắt đầu hãm: Phản ứng động của kết cấu phụ thuộc mật thiết vào vị trí xe bắt đầu đạp phanh. Nếu xe bắt đầu hãm tại vị trí cách gối tựa đầu nhịp khoảng 1/4 đến 1/2 chiều dài nhịp thứ nhất, hiệu ứng kích động sóng dao động truyền sang các nhịp lân cận là lớn nhất, đẩy hệ số động lực toàn cầu lên giá trị đỉnh.

Tác động tỷ lệ thuận của chiều cao dầm: Khảo sát các chiều cao dầm biến thiên từ 1,0 m đến 2,0 m cho thấy, chiều cao dầm càng lớn thì cánh tay đòn của lực hãm càng tăng. Dù độ cứng tổng thể của dầm tăng làm giảm chuyển vị tuyệt đối, nhưng tỷ lệ gia tăng mô men do lực hãm lại cao hơn khoảng 12% so với dầm tiết diện thấp.

Thảo luận kết quả

Cơ chế vật lý đằng sau sự gia tăng phản ứng động chính là hiện tượng chuyển dịch tải trọng động từ trục sau sang trục trước của phương tiện trong quá trình giảm tốc, kết hợp với lực ma sát tiếp tuyến tại bề mặt tiếp xúc dầm. Khi xe hãm phanh, tổng lực quán tính ngang $H = (m_1 + m_2) a_b$ tác động qua cánh tay đòn $h/2$ tạo nên một ngẫu lực uốn bổ sung, bẻ cong trục trung hòa của phần tử dầm.

Khi so sánh với các nghiên cứu trên dầm đơn giản một nhịp của các tác giả quốc tế, kết cấu dầm liên tục thể hiện đặc tính phân phối lại nội lực phức tạp hơn nhiều. Sóng uốn sinh ra từ nhịp chịu lực hãm lập tức truyền qua các gối liên tục, gây ra dao động cưỡng bức ở các nhịp lân cận ngay cả khi tải trọng chưa đi tới.

Các dữ liệu này có thể được trực quan hóa thông qua đồ thị lịch sử thời gian chuyển vị - mô men uốn và biểu đồ bề mặt 3D thể hiện mối quan hệ đa biến giữa Vận tốc - Gia tốc hãm - Hệ số động lực DAF. Bảng số liệu so sánh tần số dao động riêng giữa thuật toán MATLAB tự phát triển và mô hình phần tử hữu hạn SAP2000 khẳng định sai số tuyệt đối luôn nằm trong giới hạn kiểm soát dưới 1,0%, chứng minh độ tin cậy tuyệt đối của phương pháp đề xuất.

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

Dựa trên các phát hiện định lượng từ nghiên cứu, 4 nhóm giải pháp kỹ thuật cụ thể được đề xuất nhằm nâng cao độ an toàn và tuổi thọ cho các công trình cầu dầm liên tục:

Bổ sung hệ số động lực do lực hãm vào tiêu chuẩn thiết kế: Các cơ quan quản lý và viện tiêu chuẩn xây dựng cần kiến nghị điều chỉnh tiêu chuẩn thiết kế cầu hiện hành, nâng hệ số xung kích tải trọng thêm 15% đến 20% đối với các đoạn cầu dầm liên tục nằm trong khu vực thường xuyên xảy ra thao tác phanh như gần trạm thu phí, nút giao hoặc dốc cầu. Thời gian triển khai cập nhật nên thực hiện trong lộ trình soát xét tiêu chuẩn từ 1 đến 2 năm tới.

Thiết lập phân vùng hạn chế tốc độ và cự ly an toàn: Đơn vị quản lý vận hành giao thông cần lắp đặt biển báo giới hạn tốc độ tối đa dưới 80 km/h và quy định khoảng cách tối thiểu giữa các xe tải nặng trên các cầu dầm liên tục nhịp lớn. Biện pháp này giúp giảm thiểu 25% nguy cơ xuất hiện hiện tượng phanh đồng thời của nhiều xe tải nặng trên cùng một nhịp cầu.

Tăng cường cốt thép chịu mô men uốn âm tại vùng gối tựa: Các đơn vị tư vấn thiết kế kết cấu cầu cần chủ động tăng lượng cốt thép dọc chịu kéo tại khu vực đỉnh trụ dầm liên tục thêm 10% đến 15% so với tính toán tĩnh thông thường. Giải pháp này giúp dầm hấp thụ an toàn thành phần mô men uốn phụ do lực hãm sinh ra mà không làm phát sinh vết nứt sớm trong bê tông.

Ứng dụng công cụ số hóa trong thẩm tra và quan trắc sức khỏe công trình: Các trung tâm kiểm định và tư vấn giám sát nên ứng dụng module tính toán động lực học dựa trên thuật toán tích phân Newmark đã phát triển để mô phỏng tải trọng hãm thực tế trong quy trình thử tải nghiệm thu cầu, đồng thời tích hợp vào hệ thống quan trắc cảm biến dao động thời gian thực với chu kỳ đánh giá định kỳ 6 tháng một lần.

Đố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 giá trị ứng dụng thực tiễn và học thuật cao cho 4 nhóm đối tượng trọng tâm:

Kỹ sư thiết kế cầu đường và công trình giao thông: Tài liệu cung cấp cơ sở tính toán chi tiết để kiểm tra khả năng chịu lực của dầm liên tục nhiều nhịp khi chịu tác động phanh xe, giúp kỹ sư lựa chọn kích thước hình học và bố trí cốt thép chịu uốn chính xác hơn tại các vị trí nhạy cảm.

Học viên cao học, nghiên cứu sinh và giảng viên ngành Xây dựng: Luận văn là tài liệu tham khảo mẫu mực về quy trình rời rạc hóa kết cấu theo phương pháp phần tử hữu hạn, cách thiết lập phương trình vi phân tương tác động lực học xe - cầu và phương pháp lập trình thuật toán số Newmark trên nền tảng MATLAB.

Đơn vị quản lý, khai thác và bảo trì công trình hạ tầng giao thông: Nắm bắt được cơ chế phá hoại mỏi và vùng xuất hiện ứng suất cực đại do lực phanh gây ra, từ đó xây dựng kế hoạch duy tu, sửa chữa mặt cầu và gối cầu định kỳ theo đúng trọng tâm kỹ thuật.

Chuyên gia thẩm tra, kiểm định chất lượng công trình xây dựng: Sử dụng phương pháp và số liệu trong luận văn làm công cụ đối chứng độc lập khi đánh giá độ an toàn kết cấu của các dự án cầu dầm liên tục chịu lưu lượng xe tải nặng mật độ cao.

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

Tại sao cần xét đến lực hãm xe thay vì chỉ tính toán tải trọng di động với vận tốc không đổi? Khi xe phanh gấp, lực quán tính ngang sinh ra lực ma sát lớn tiếp xúc với mặt cầu, tạo thành ngẫu lực uốn cục bộ làm tăng hệ số động lực DAF thêm 15% đến 35%. Nếu chỉ tính tải trọng chạy đều, kết cấu sẽ bị đánh giá thiếu an toàn, đặc biệt là tại vùng gối tựa dầm liên tục.

Mô hình Sprung-mass 2 bậc tự do có ưu điểm gì vượt trội so với mô hình tải trọng di động đơn giản? Mô hình Sprung-mass mô tả đầy đủ chuyển vị độc lập của thân xe và bánh xe qua hệ lò xo và cản nhớt. Nhờ đó, mô hình phản ánh chính xác tương tác dao động 2 chiều giữa mui xe và mặt dầm, giúp nắm bắt hiện tượng dằn xóc và biến thiên lực nén thực tế của bánh xe lên mặt cầu.

Gia tốc hãm ảnh hưởng như thế nào đến độ võng và mô men uốn của dầm liên tục? Gia tốc hãm tỷ lệ thuận với lực quán tính ngang của phương tiện. Khi gia tốc hãm tăng từ $-2,5\text{ m/s}^2$ lên $-5,0\text{ m/s}^2$, mô men uốn cực đại trong dầm tăng thêm khoảng 20%, tạo ra xung lực uốn đột ngột làm dầm võng sâu hơn và dao động kéo dài hơn sau khi xe đã dừng hẳn.

Phương pháp tích phân Newmark được áp dụng trong đề tài có ưu điểm gì về mặt tính toán? Phương pháp Newmark với các tham số chuẩn $\gamma = 0,5$ và $\beta = 0,25$ đảm bảo tính ổn định vô điều kiện trên toàn miền thời gian. Thuật toán cho phép giải chính xác hệ phương trình vi phân phi tuyến bậc hai với bước thời gian siêu nhỏ 0,001 giây mà không bị tích lũy sai số số học.

Chiều cao tiết diện dầm có vai trò gì đối với tác động của lực hãm xe? Chiều cao dầm $h$ quyết định độ lớn của cánh tay đòn từ mặt cầu đến trục trung hòa. Khi chiều cao dầm tăng lên, cánh tay đòn $h/2$ tăng tương ứng, khiến mô men tập trung do lực ma sát hãm sinh ra tăng theo, đòi hỏi việc tính toán cốt thép phải xét đến ảnh hưởng của chiều cao kết cấu.

Kết luận

Nghiên cứu về phân tích động lực học dầm liên tục chịu tải di động có xét đến lực hãm mang lại những đóng góp khoa học và thực tiễn vững chắc:

  • Xây dựng thành công mô hình toán học tương tác xe - dầm liên tục toàn diện, tích hợp đầy đủ hệ dao động Sprung-mass 2 bậc tự do và lực hãm 2 pha theo thời gian thực.
  • Thiết lập hoàn chỉnh hệ phương trình chuyển động bằng phương pháp phần tử hữu hạn, tính toán chuẩn xác các ma trận khối lượng, độ cứng và cản tương tác.
  • Phát triển bộ công cụ tính toán số tự động hóa bằng MATLAB dựa trên thuật toán Newmark, đạt độ chính xác cao khi kiểm chứng với SAP2000 và các công bố quốc tế.
  • Khẳng định định lượng mức độ nguy hiểm của lực hãm khi làm tăng hệ số động lực DAF từ 15% đến 38%, cung cấp bằng chứng thực nghiệm số thuyết phục cho ngành xây dựng cầu.
  • Đề xuất các giải pháp kỹ thuật khả thi về điều chỉnh tiêu chuẩn thiết kế và kiểm soát giao thông nhằm bảo đảm an toàn kết cấu công trình.

Trong giai đoạn 12 đến 24 tháng tới, hướng phát triển tiếp theo của đề tài là mở rộng mô hình sang không gian 3 chiều (3D), khảo sát đoàn xe nhiều trục và xét thêm ảnh hưởng của mấp mô mặt đường ngẫu nhiên. Đây là tài liệu khoa học giá trị cao mà các kỹ sư kết cấu, nhà nghiên cứu và cơ quan quản lý xây dựng nên tham khảo ứng dụng ngay vào thực tiễn thiết kế và thẩm định công trình.