Tổng quan nghiên cứu

Hạ tầng giao thông hiện đại như đường cao tốc và đường cất hạ cánh sân bay đóng vai trò huyết mạch trong việc kết nối kinh tế toàn cầu. Tại các quốc gia phát triển như Mỹ, Nhật Bản, Đức cũng như tại Việt Nam, các kết cấu mặt đường bê tông liên tục phải chịu áp lực từ các phương tiện vận tải có tải trọng lớn và mật độ lưu thông cao. Vấn đề cốt lõi đặt ra là kết cấu mặt đường vừa chịu tải trọng động di chuyển với vận tốc cao, vừa chịu tác động khắc nghiệt từ sự biến thiên nhiệt độ môi trường. Sự kết hợp này gây ra hiện tượng uốn nhiệt, làm phát sinh ứng suất nhiệt lớn cùng biến dạng dư, dẫn đến nứt vỡ và suy giảm tuổi thọ công trình.

Luận văn thạc sĩ chuyên ngành Kỹ thuật xây dựng mã số 8580201 thực hiện tại Trường Đại học Bách khoa, Đại học Quốc gia Thành phố Hồ Chí Minh, hoàn thành vào tháng 01 năm 2023, tập trung giải quyết bài toán động lực học phức tạp này. Mục tiêu nghiên cứu cụ thể là xây dựng mô hình toán học và thuật toán số hoàn chỉnh bằng phương pháp phần tử chuyển động để phân tích ứng xử động của kết cấu tấm Mindlin đặt trên nền đàn nhớt Pasternak. Nghiên cứu thực hiện khảo sát toàn diện qua 13 bài toán mô phỏng số, thiết lập hệ ma trận 27 bậc tự do cho mỗi phần tử tứ giác 9 nút đẳng tham số. Kết quả nghiên cứu cung cấp cơ sở khoa học chính xác, giúp giảm thiểu sai số tính toán võng động xuống dưới 2% so với các phương pháp giải tích truyền thống, đóng góp trực tiếp vào công tác tính toán thiết kế và duy tu 100% các tuyến cao tốc trọng điểm.

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 tích hợp hệ thống 3 nền tảng lý thuyết cơ học kết cấu tiên tiến nhằm phản ánh trung thực trạng thái làm việc của tấm mặt đường trong thực tế:

  • Lý thuyết tấm Mindlin (Reissner-Mindlin Plate Theory): Khắc phục triệt để nhược điểm của lý thuyết tấm mỏng cổ điển Kirchhoff bằng cách tính đến ảnh hưởng của biến dạng trượt ngang với hệ số hiệu chỉnh cắt là 5/6. Lý thuyết này cho phép các đoạn thẳng vuông góc với mặt trung hòa trước biến dạng vẫn giữ nguyên độ thẳng nhưng không nhất thiết vuông góc sau khi chịu lực, đặc biệt chính xác cho kết cấu tấm có tỷ lệ chiều dày trên nhịp lớn hơn 0.05.
  • Mô hình nền đàn nhớt Pasternak hai thông số: Mô hình hóa phản lực của nền đất thông qua hệ số độ cứng nền theo phương đứng và hệ số sức kháng cắt của lớp đệm, kết hợp hệ số cản nhớt nhằm mô tả đầy đủ khả năng tiêu tán năng lượng dao động của nền đường.
  • Lý thuyết nhiệt đàn hồi tổng quát: Mô tả trường nhiệt độ biến thiên tuyến tính theo chiều dày tấm với 2 thành phần chính là nhiệt độ trung bình và độ chênh lệch nhiệt độ giữa hai mặt tấm, từ đó xác định chính xác mô men uốn nhiệt và lực dọc nhiệt trong mặt phẳng.
  • Phương pháp phần tử chuyển động (Moving Element Method - MEM): Chuyển đổi bài toán từ hệ tọa độ cố định sang hệ tọa độ gắn liền và chuyển động cùng tải trọng theo trục di chuyển với vận tốc không đổi, giúp miền tính toán luôn giữ nguyên dạng tĩnh trong không gian khảo sát.

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 việc tự phát triển mã nguồn chương trình tính toán trên nền tảng ngôn ngữ lập trình Matlab.

  • Quy mô mẫu và thiết kế mô phỏng: Cỡ mẫu khảo sát gồm 13 bài toán số chuyên sâu, bao gồm 3 bài toán kiểm chứng độ tin cậy và 10 bài toán phân tích tham số động lực học. Miền tấm được rời rạc hóa với các mức lưới phần tử từ 4x4 đến 16x16 phần tử tứ giác 9 nút đẳng tham số, tạo ra hệ phương trình động lực học tổng thể với hàng trăm bậc tự do động.
  • Phương pháp chọn mẫu tham số: Lựa chọn mẫu tham số có chủ đích dựa trên các tiêu chuẩn thiết kế sân bay và đường cao tốc thực tế, bao gồm dải vận tốc chuyển động từ 10 m/s đến 100 m/s, dải tải trọng tập trung từ 50 kN đến 250 kN, và dải chênh lệch nhiệt độ từ 0 độ C đến 50 độ C.
  • Lý do lựa chọn phương pháp phân tích: Phương pháp tích phân số Newmark từng bước trong miền thời gian với tham số chuẩn không tiêu tán số được lựa chọn nhờ tính ổn định vô điều kiện và tốc độ hội tụ cao. Việc sử dụng phép tích phân Gauss với 9 điểm cầu phương bảo đảm ma trận độ cứng và khối lượng đạt độ chính xác tuyệt đối. Toàn bộ quy trình nghiên cứu được triển khai chặt chẽ từ tháng 07 năm 2022 đến tháng 12 năm 2022.

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

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

Nghiên cứu đã đưa ra các kết luận định lượng rõ ràng về ứng xử động lực học của tấm mặt đường:

  • Ảnh hưởng của gradient nhiệt độ: Khi độ chênh lệch nhiệt độ giữa hai mặt tấm tăng từ 10 độ C lên 40 độ C, độ võng động lớn nhất tại tâm tấm tăng thêm khoảng 22% đến 38%. Hiện tượng này xảy ra do mô men uốn nhiệt gây ra độ cong ban đầu, làm suy giảm đáng kể độ cứng chống uốn tổng thể của kết cấu tấm.
  • Tác động của vận tốc và độ lớn tải trọng di động: Chuyển vị động của tấm tỷ lệ thuận với cường độ tải trọng. Khi vận tốc di chuyển tăng dần từ 20 m/s đến vận tốc tới hạn khoảng 80 m/s, hệ số khuếch đại động lực học của độ võng tăng vọt lên hơn 45%, trước khi có xu hướng ổn định ở dải vận tốc siêu cao.
  • Vai trò củng cố của chiều dày tấm và lớp nền: Khi tăng chiều dày tấm thêm 50% (từ 0.2 m lên 0.3 m), độ võng động cực đại giảm mạnh khoảng 35%. Đồng thời, việc tăng hệ số kháng cắt của móng Pasternak lên gấp 2 lần giúp giảm biên độ dao động động lực học tới 28%.
  • Độ tin cậy của mô hình MEM: Kết quả kiểm chứng 6 mode dao động tự do đầu tiên của tấm Mindlin bằng chương trình tính toán Matlab cho thấy sai số so với các công bố quốc tế uy tín chỉ dao động trong khoảng 0.5% đến 1.4%, khẳng định độ chính xác vượt trội của thuật toán.

Thảo luận kết quả

Các kết quả phân tích phản ánh chính xác bản chất cơ học tương tác động giữa kết cấu và môi trường. Sự thay đổi nhiệt độ tạo ra ứng suất nén nhiệt trong mặt phẳng tấm, làm xuất hiện hiệu ứng suy giảm độ cứng hình học, từ đó khiến tấm dễ bị uốn võng hơn dưới tác dụng của bánh xe hoặc bánh đáp máy bay. So với các nghiên cứu dùng phương pháp phần tử hữu hạn cổ điển, phương pháp phần tử chuyển động MEM loại bỏ hoàn toàn hiện tượng dao động số giả tạo sinh ra do việc chia lại lưới phần tử theo bước thời gian.

Các dữ liệu phân tích số trong luận văn có thể được trình bày một cách trực quan thông qua biểu đồ lịch sử thời gian của chuyển vị tâm tấm theo từng bước lặp thời gian và bảng so sánh chuyển vị cực đại ứng với 5 dải nhiệt độ và 4 cấp vận tốc khác nhau. Việc trực quan hóa dữ liệu qua biểu đồ 2D và bảng tổng hợp ma trận tham số giúp các kỹ sư dễ dàng nhận diện vùng vận tốc nguy hiểm và ngưỡng nhiệt độ giới hạn cần can thiệp kết cấu.

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

Dựa trên các kết quả đạt được, luận văn đưa ra 4 nhóm giải pháp thực tiễn có tính khả thi cao:

  • Tích hợp thuật toán MEM vào phần mềm thiết kế chuyên dụng: Bộ Giao thông Vận tải và các Viện nghiên cứu chuyên ngành cần chủ trì xây dựng bộ công cụ tính toán tự động hóa dựa trên phương pháp phần tử chuyển động trong giai đoạn 12 đến 18 tháng tới. Mục tiêu giảm thiểu 40% thời gian chạy mô phỏng và cắt giảm sai số thiết kế xuống dưới mức 3%.
  • Kiểm soát và giới hạn gradient nhiệt độ mặt đường: Các Ban quản lý dự án và Tổng công ty Cảng hàng không cần ban hành quy chuẩn thi công lớp phủ cách nhiệt hoặc rãnh thoát nhiệt cho đường cất hạ cánh trong vòng 6 tháng tới. Mục tiêu duy trì độ chênh lệch nhiệt độ giữa hai bề mặt tấm luôn dưới 15 độ C nhằm hạn chế 25% nguy cơ nứt gãy do mỏi nhiệt.
  • Nâng cao chỉ tiêu kháng cắt của lớp móng: Các nhà thầu thi công hạ tầng giao thông cần áp dụng các vật liệu gia cố xi măng hoặc lưới địa kỹ thuật nhằm nâng hệ số kháng cắt của lớp nền móng lên tối thiểu 30% so với thiết kế truyền thống. Giải pháp này cần triển khai ngay trong giai đoạn thi công 2024-2025 để giảm 20% độ võng đàn nhớt lâu dài.
  • Bố trí hệ thống cảm biến quan trắc sức khỏe công trình: Các đơn vị vận hành đường cao tốc cần lắp đặt hệ thống quan trắc tải trọng động và nhiệt độ mặt đường theo thời gian thực tại các vị trí xung yếu trong lộ trình 24 tháng, bảo đảm cảnh báo sớm các nguy cơ mất an toàn khi phương tiện vượt tải trọng cho phép.

Đố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ị thiết thực cho 4 nhóm đối tượng chính:

  • Kỹ sư thiết kế kết cấu cầu đường và sân bay: Ứng dụng trực tiếp công thức và bảng tra số liệu để tính toán chiều dày tấm bê tông mặt đường cho các dự án đường cao tốc và đường băng sân bay tiếp nhận các loại máy bay thân rộng có tải trọng trục trên 50 tấn.
  • Nhà thầu thi công và đơn vị tư vấn giám sát: Sử dụng kết quả phân tích ứng xử nhiệt để lập biện pháp thi công đổ bê tông khối lớn hợp lý, kiểm soát nứt do ứng suất nhiệt trong 28 ngày đầu đóng rắn.
  • Học viên cao học và nghiên cứu sinh ngành Xây dựng: Tiếp cận một tài liệu học thuật mẫu mực về phương pháp số nâng cao, nắm vững kỹ thuật lập trình ma trận khối lượng và độ cứng cho phần tử chuyển động 9 nút trên Matlab.
  • Cơ quan quản lý và hoạch định chính sách hạ tầng: Sử dụng các khuyến nghị làm căn cứ cập nhật tiêu chuẩn bảo trì định kỳ cho mạng lưới hơn 1000 km đường cao tốc đang vận hành trên toàn quốc.

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

  • Phương pháp phần tử chuyển động MEM có ưu điểm gì vượt trội so với FEM truyền thống? Phương pháp phần tử chuyển động thiết lập hệ tọa độ di chuyển đồng tốc với tải trọng, giúp tải trọng luôn cố định trong miền tính toán. Nhờ đó, số lượng phần tử cần mô hình hóa giảm hơn 60%, loại bỏ hoàn toàn sai số rời rạc hóa thời gian mà FEM truyền thống thường gặp phải.
  • Tại sao lý thuyết tấm Mindlin lại được ưu tiên lựa chọn thay vì lý thuyết Kirchhoff? Lý thuyết Mindlin có xét đến biến dạng cắt trượt ngang với hệ số hiệu chỉnh 5/6, cực kỳ chính xác cho các kết cấu tấm bê tông mặt đường có bề dày từ 0.2 m đến 0.4 m, giúp tránh hiện tượng khóa trượt và đánh giá đúng ứng suất tại góc mép tấm.
  • Sự chênh lệch nhiệt độ giữa hai mặt tấm ảnh hưởng như thế nào đến độ võng kết cấu? Khi nhiệt độ mặt trên cao hơn mặt dưới từ 20 độ C đến 40 độ C, tấm bị uốn vồng lên ở giữa, làm tăng độ võng động tổng cộng thêm 15% đến 35% khi xe chạy qua, khiến nguy cơ nứt mỏi mặt đường gia tăng gấp 2 lần.
  • Mô hình nền Pasternak phản ánh đất nền tốt hơn nền Winkler ở điểm nào? Nền Pasternak bổ sung thêm lớp truyền cắt giữa các lò xo đàn hồi với hệ số kháng cắt riêng biệt. Điều này khắc phục nhược điểm biến dạng độc lập cục bộ của nền Winkler, mô tả chính xác sự liên tục của nền đất thực tế với sai số dưới 3%.
  • Thuật toán Newmark sử dụng trong luận văn có bảo đảm độ ổn định không? Thuật toán Newmark với các tham số gia tốc trung bình không đổi bảo đảm tính ổn định vô điều kiện trong miền thời gian, cho phép bước lặp tích phân đạt độ chính xác cao và hội tụ mượt mà sau 6 mode dao động kiểm chứng.

Kết luận

  • Xây dựng thành công mô hình toán học và phần tử chuyển động Q9 gồm 27 bậc tự do cho kết cấu tấm Mindlin chịu tải trọng di động có xét trường nhiệt độ.
  • Thiết lập hoàn chỉnh mã nguồn tính toán động lực học trên phần mềm Matlab, đạt độ chính xác cao với sai số kiểm chứng dưới 1.5%.
  • Khảo sát chi tiết 13 kịch bản số, làm rõ quy luật gia tăng độ võng động theo vận tốc tải trọng và gradient nhiệt độ môi trường.
  • Đề xuất định lượng các giải pháp kỹ thuật nhằm tối ưu hóa chiều dày tấm và tăng cường khả năng chịu cắt của lớp móng công trình giao thông.
  • Đóng góp cơ sở lý thuyết và công cụ số vững chắc cho việc thiết kế kết cấu mặt đường chịu tải trọng nặng trong điều kiện biến đổi khí hậu.

Kế hoạch nghiên cứu tiếp theo trong 12 đến 24 tháng tới sẽ mở rộng mô hình cho vật liệu tấm composite biến đổi chức năng FGM và xét đến tải trọng xe nhiều trục ngẫu nhiên. Các cơ quan quản lý và doanh nghiệp xây dựng hạ tầng quan tâm có thể liên hệ hợp tác ứng dụng mô hình số này vào các dự án giao thông trọng điểm ngay hôm nay.