Tổng quan nghiên cứu

Trong bối cảnh hạ tầng giao thông hiện đại phát triển vượt bậc, lưu lượng phương tiện vận tải nặng với tải trọng trục từ 10 đến 15 tấn cùng các dòng máy bay thương mại có tải trọng cất hạ cánh hàng trăm tấn ngày càng gia tăng. Khi các phương tiện này lưu thông với tốc độ từ 80 đến 120 km/h trên đường cao tốc hoặc tiếp đất với vận tốc trên 140 km/h trên đường băng, hiện tượng phanh hãm đột ngột thường xuyên xảy ra. Quá trình hãm khẩn cấp tạo ra lực ma sát tiếp xúc và lực quán tính ngang rất lớn, làm gia tăng độ võng động và ứng suất cắt cục bộ trong kết cấu mặt đường từ 18% đến 28% so với trạng thái chuyển động đều.

Vấn đề đặt ra là các phương pháp tính toán truyền thống phần lớn chỉ khảo sát kết cấu tấm đơn trên nền đàn hồi chịu tải trọng di động với vận tốc không đổi, bỏ qua tương tác phức tạp giữa các lớp kết cấu và hiệu ứng hãm phanh. Mục tiêu cốt lõi của nghiên cứu này là xây dựng mô hình toán học và giải thuật số nhằm phân tích chính xác ứng xử động lực học của kết cấu tấm nhiều lớp trên nền đàn nhớt chịu tải trọng di động có xét đến hiện tượng hãm, thông qua việc phát triển phương pháp phần tử tấm nhiều lớp chuyển động (Multi-layer Moving Plate Method - MMPM).

Phạm vi nghiên cứu tập trung vào mô hình mặt đường sân bay và đường cao tốc gồm tấm bê tông xi măng dày 0.22 m, lớp liên kết nhựa đường dày 0.03 m và lớp móng xi măng đá dày 0.18 m đặt trên nền đất đàn nhớt nhiều lớp, được thực hiện trong giai đoạn từ tháng 09/2019 đến tháng 12/2020 tại Trường Đại học Bách Khoa - ĐHQG TP. Hồ Chí Minh. Nghiên cứu mang lại ý nghĩa khoa học và thực tiễn to lớn khi giúp giảm thiểu hơn 40% thời gian tính toán so với phương pháp phần tử hữu hạn thông thường, đồng thời cung cấp cơ sở dữ liệu tin cậy để tối ưu hóa chiều dày kết cấu, nâng cao tuổi thọ công trình thêm 20% đến 30%.

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 chặt chẽ giữa hai khung lý thuyết cơ học nền tảng:

Lý thuyết tấm dày Reissner-Mindlin: Khác với lý thuyết tấm mỏng Kirchhoff, mô hình Reissner-Mindlin xét đế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 (tương đương khoảng 0.833). Trường chuyển vị tại mỗi điểm trong tấm được biểu diễn qua ba thành phần chuyển vị độc lập gồm độ võng đứng và hai góc xoay của pháp tuyến mặt trung hòa quanh các trục tọa độ. Điều này cho phép mô tả chính xác trạng thái ứng suất - biến dạng ba chiều của các lớp tấm bê tông và tấm xi măng đá có chiều dày lớn.

Mô hình nền đàn nhớt nhiều lớp: Nền đất được mô hình hóa theo trường phái Winkler cải tiến kết hợp các đặc trưng động lực học của Richart và Lysmer (1970) cùng hệ số hiệu chỉnh độ chôn sâu của Whitman (1972). Mô hình này xác định hệ số độ cứng chống nén và hệ số cản nhớt tương đương theo phương đứng dựa trên mô đun kháng cắt và hệ số Poisson trung bình của các lớp đất nền.

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

  • Phương pháp phần tử tấm nhiều lớp chuyển động (MMPM): Kỹ thuật chuyển đổi hệ phương trình vi phân chuyển động từ hệ tọa độ cố định sang hệ tọa độ chuyển động gắn liền với vị trí của tải trọng hãm.
  • Động học tải trọng hãm hai pha: Quá trình tác dụng lực được chia thành pha 1 (chuyển động với vận tốc ban đầu không đổi) và pha 2 (giảm tốc với gia tốc hãm cho đến khi dừng hẳn).
  • Quy đổi tải trọng hãm thành mômen nút: Lực ma sát trượt ngang giữa bánh xe và mặt đường được tích hợp vào các bậc tự do của phần tử tấm dưới dạng các mômen tập trung tương đương.

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

Nguồn dữ liệu đầu vào sử dụng các thông số kỹ thuật tiêu chuẩn của kết cấu mặt đường nhiều lớp: Lớp bê tông trên cùng có mô đun đàn hồi 3.15 x 10^10 N/m2, khối lượng riêng 2400 kg/m3 và hệ số Poisson 0.15; lớp xi măng đá bên dưới có mô đun đàn hồi 1.25 x 10^10 N/m2, khối lượng riêng 2300 kg/m3 và hệ số Poisson 0.20; các lớp đất nền có mô đun biến dạng biến thiên từ 2.35 x 10^7 N/m2 đến 5.0 x 10^7 N/m2.

Phương pháp phân tích số được xây dựng dựa trên việc rời rạc hóa kết cấu thành các phần tử tứ giác 9 nút (Q9) đẳng tham số hai lớp. Mỗi phần tử bao gồm 18 nút với 3 bậc tự do tại mỗi nút, tạo thành ma trận 54 bậc tự do cho mỗi phần tử liên hợp. Tích phân số trên miền phần tử được tính toán chính xác thông qua phép cầu phương Gauss với lưới tích phân 3x3 điểm Gauss. Để khảo sát độ hội tụ và độ chính xác của nghiệm, mô hình thử nghiệm trên các mạng lưới chia từ 20x20 đến 60x60 phần tử. Mạng lưới 60x60 phần tử được lựa chọn nhằm tối ưu hóa giữa độ chính xác tuyệt đối và chi phí bộ nhớ máy tính.

Hệ phương trình vi phân dao động tổng thể bậc hai theo thời gian được giải bằng thuật toán tích phân từng bước Newmark với giả thiết gia tốc trung bình (hệ số gamma = 0.5 và beta = 0.25), bước thời gian được chọn từ 0.001 s đến 0.005 s. Toàn bộ thuật toán được lập trình tự động hóa trên ngôn ngữ MATLAB trong khung thời gian 15 tháng nghiên cứu. Lý do lựa chọn phương pháp MMPM kết hợp Newmark là vì giải pháp này loại bỏ hoàn toàn hiện tượng phản xạ sóng tại biên giả tạo và không cần liên tục tái tạo lưới phần tử khi tải trọng di chuyển, đem lại hiệu năng vượt trội so với các phần mềm phần tử hữu hạn thương mại.

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ố và phân tích độ nhạy của các tham số đầu vào đã chỉ ra 4 phát hiện quan trọng:

Ảnh hưởng của gia tốc hãm: Khi gia tốc hãm tăng từ -5 m/s2 lên -20 m/s2 tại cùng mức vận tốc ban đầu 27.78 m/s (khoảng 100 km/h), độ võng cực đại tại tâm tấm bê tông trên cùng tăng từ 0.385 mm lên 0.467 mm (tương ứng mức tăng 21.3%). Đồng thời, ứng suất uốn cục bộ tại vị trí bánh xe trượt tăng khoảng 26.3% so với trạng thái xe chạy đều.

Ảnh hưởng của vận tốc ban đầu trước khi hãm: Khảo sát dải vận tốc ban đầu từ 20 m/s (72 km/h) đến 40 m/s (144 km/h) dưới gia tốc hãm cố định -10 m/s2 cho thấy độ võng lớn nhất của tấm trên tăng thêm 14.2%, đồng thời vị trí xuất hiện chuyển vị cực đại bị dịch chuyển tịnh tiến dọc theo chiều chuyển động khoảng 3.5 m đến 5.2 m.

Hiệu ứng tương tác giữa hai lớp tấm: Độ cứng liên kết và hệ số cản của lớp dính bám nhựa đường đóng vai trò quyết định đến sự phân phối tải trọng. Khi độ dính bám giữa hai lớp bị suy giảm, độ chênh lệch chuyển vị giữa tấm bê tông trên và tấm xi măng đá dưới tăng vọt lên 35.8%, khiến tấm trên phải chịu tải trọng uốn cục bộ cao hơn 29.4%. Ngược lại, khi liên kết dính bám hoàn hảo, hai lớp tấm cùng làm việc đồng thời, giúp giảm ứng suất kéo khi uốn ở đáy tấm bê tông khoảng 31.2%.

Ảnh hưởng của hệ số ma sát động: Khi hệ số ma sát động giữa lốp xe và mặt đường tăng từ 0.3 lên 0.7 trong điều kiện phanh khẩn cấp, thành phần lực quán tính ngang chuyển đổi thành mômen uốn nút tăng hơn gấp đôi, làm biên độ dao động động lực học của tấm tăng thêm 22.8%.

Thảo luận kết quả

Nguyên nhân vật lý của các hiện tượng trên bắt nguồn từ sự chuyển hóa đột ngột của động năng thành công của lực ma sát và lực quán tính ngang trong thời gian rất ngắn. Lực hãm tác dụng lệch tâm so với mặt trung hòa của tấm tạo ra một mômen xoắn và mômen uốn phụ tác động trực tiếp vào các nút phần tử, kết hợp với tải trọng thẳng đứng gây ra hiệu ứng cộng hưởng cục bộ.

Khi so sánh với các nghiên cứu của Huang và Thambiratnam (2001), kết quả từ phương pháp MMPM cho thấy sự tương đồng tuyệt đối với sai số chuyển vị chỉ dưới 1.85%. So sánh với các công trình nghiên cứu của các tác giả trong nước như Tư (2016) và Nhi (2015), giải thuật MMPM cho kết quả hội tụ mượt mà hơn và phản ánh trung thực trạng thái làm việc đa tầng của kết cấu mặt đường.

Các kết quả nghiên cứu trên được biểu diễn trực quan thông qua các đồ thị lịch sử chuyển vị theo thời gian thực, biểu đồ phân bố đường cong chuyển vị dọc theo trục x và y của tấm, cùng các bảng so sánh sai số tương đối giữa các mật độ lưới phần tử. Các đường cong dao động tắt dần thể hiện rõ vai trò tiêu tán năng lượng của hệ số cản đất nền và lớp cản nhớt liên kết giữa hai tấm.

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

Dựa trên các kết luận rút ra từ phân tích số, luận văn đề xuất 4 nhóm giải pháp kỹ thuật mang tính ứng dụng cao:

Tăng cường độ dính bám kháng cắt của lớp trung gian: Chủ đầu tư và các đơn vị thi công cần áp dụng các loại nhũ tương nhựa đường cải tính polyme hoặc màng chống thấm chuyên dụng nhằm duy trì hệ số độ cứng liên kết tối thiểu 1.5 x 10^8 N/m3, mục tiêu giảm thiểu 25% nguy cơ trượt tách lớp giữa tấm bê tông và lớp móng trong vòng 12 tháng sau khi đưa vào khai thác.

Gia cường cục bộ tại các phân đoạn có xác suất phanh hãm cao: Các đơn vị tư vấn thiết kế cần tăng chiều dày tấm bê tông mặt đường từ 0.22 m lên 0.26 m (tăng khoảng 18%) hoặc bổ sung lưới cốt thép tăng cường tại vùng tiếp đất của đường băng (300 m đến 500 m đầu đường băng) và khu vực trước trạm thu phí, nút giao cao tốc. Giải pháp này giúp triệt tiêu hơn 20% ứng suất uốn phát sinh do lực hãm đột ngột.

Tích hợp thuật toán MMPM vào quy trình tính toán thiết kế chuẩn: Đề xuất Bộ Giao thông Vận tải và Viện Khoa học Công nghệ Giao thông Vận tải xem xét bổ sung phương pháp phần tử tấm nhiều lớp chuyển động vào tiêu chuẩn thiết kế áo đường cứng giai đoạn 2021-2025, thay thế cho các mô hình tính toán tĩnh đơn giản hóa hiện nay.

Kiểm soát chất lượng độ nhám và thoát nước mặt đường: Đơn vị quản lý khai thác cần định kỳ đo đạc và duy trì hệ số ma sát mặt đường trong ngưỡng tối ưu từ 0.5 đến 0.6, kết hợp vệ sinh cao su bám dính trên đường băng để hạn chế hiện tượng trượt lết bánh xe, giảm áp lực quán tính đột ngột lên kết cấu nền đường.

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

Công trình nghiên cứu mang giá trị thực tiễn và học thuật sâu sắc, đặc biệt hữu ích cho 4 nhóm đối tượng sau:

Kỹ sư thiết kế kết cấu hạ tầng giao thông và sân bay: Nắm vững phương pháp tính toán tải trọng động và lực hãm thực tế để tối ưu hóa chiều dày các lớp vật liệu, tiết kiệm từ 15% đến 20% chi phí đầu tư ban đầu mà vẫn đảm bảo an toàn chịu lực.

Các nhà khoa học và nghiên cứu sinh chuyên ngành Cơ học công trình: Sử dụng bộ khung lý thuyết MMPM và phần tử tấm Q9 54 bậc tự do làm nền tảng mở rộng cho các bài toán phức tạp hơn như tấm composite nhiều lớp, nền phi tuyến hoặc tải trọng đoàn xe lưu thông liên tục.

Đơn vị quản lý, khai thác và bảo trì đường cao tốc, sân bay: Ứng dụng mô hình để chẩn đoán nguyên nhân gây nứt vỡ, trồi sụt mặt đường tại các vị trí phanh gấp, từ đó lập kế hoạch duy tu định kỳ chính xác, giảm 30% nguy cơ hỏng hóc đột xuất.

Học viên cao học và sinh viên ngành Xây dựng Cầu đường, Công trình Dân dụng: Tiếp cận tài liệu tham khảo mẫu mực về kỹ thuật lập trình MATLAB giải bài toán động lực học công trình và phương pháp tích phân Newmark trong cơ học kết cấu.

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

Phương pháp MMPM có điểm gì ưu việt hơn so với phương pháp phần tử hữu hạn FEM truyền thống? Phương pháp MMPM thiết lập hệ phương trình chuyển động trong hệ tọa độ di động cùng vận tốc với tải trọng. Điều này giúp vị trí đặt tải luôn cố định trong miền tính toán, loại bỏ nhu cầu chia lại lưới và tránh hiện tượng phản xạ sóng tại biên, giúp tiết kiệm hơn 40% thời gian xử lý dữ liệu.

Tại sao bắt buộc phải xét đến hiện tượng hãm khẩn cấp khi thiết kế mặt đường cao tốc và đường băng? Hiện tượng hãm phanh tạo ra lực ma sát trượt ngang và lực quán tính rất lớn, làm xuất hiện mômen uốn cục bộ làm tăng độ võng và ứng suất uốn của tấm mặt đường thêm 18% đến 28% so với khi xe chạy đều, dễ gây nứt gãy nếu không được tính toán dự phòng.

Lớp nhựa đường mỏng nằm giữa tấm bê tông và lớp xi măng đá có vai trò kết cấu như thế nào? Lớp nhựa đường đóng vai trò là tầng liên kết đàn nhớt giúp truyền và phân phối ứng suất từ tấm trên xuống tấm dưới, đồng thời triệt tiêu dao động nhờ hệ số cản nhớt. Nếu mất liên kết này, độ võng tấm trên sẽ tăng vọt thêm 35.8%.

Mô hình nền đất nhiều lớp trong nghiên cứu được xác định và tính toán ra sao? Nền đất được mô hình hóa bằng hệ lò xo và giảm chấn đàn nhớt Winkler cải tiến. Các thông số độ cứng và độ cản được xác định dựa trên lý thuyết của Richart và Lysmer (1970) kết hợp hệ số hiệu chỉnh độ chôn sâu của Whitman (1972) qua các giá trị mô đun đàn hồi trung bình của từng lớp đất.

Thuật toán Newmark trong nghiên cứu được kiểm chứng độ chính xác và tin cậy như thế nào? Thuật toán Newmark với bước thời gian 0.001 s được đối chiếu trực tiếp với nghiệm giải tích và kết quả công bố của Huang và Thambiratnam (2001). Độ sai lệch ghi nhận luôn duy trì dưới 1.85%, chứng minh thuật toán có độ ổn định và độ tin cậy cơ học tuyệt đối.

Kết luận

  • Luận văn đã phát triển thành công phương pháp phần tử tấm nhiều lớp chuyển động (MMPM) dựa trên lý thuyết tấm dày Reissner-Mindlin để giải quyết toàn diện bài toán động lực học tấm nhiều lớp trên nền đàn nhớt.
  • Thiết lập hoàn chỉnh ma trận khối lượng, ma trận độ cứng, ma trận cản và véc tơ tải trọng có xét đến hiện tượng hãm hai pha và quy đổi lực ma sát ngang thành mômen nút tương đương.
  • Xây dựng chương trình tính toán số hóa tự động trên môi trường MATLAB, chứng minh tốc độ hội tụ nhanh vượt trội và độ chính xác cao với sai số dưới 1.85% so với các tài liệu chuẩn quốc tế.
  • Khẳng định các yếu tố gia tốc hãm, vận tốc ban đầu và chất lượng liên kết dính bám giữa các lớp có ảnh hưởng mang tính quyết định đến tuổi thọ và độ bền mỏi của kết cấu mặt đường.
  • Mở ra hướng nghiên cứu tiếp theo trong giai đoạn 2021-2025 về bài toán tấm trên nền phi tuyến, xét đến ảnh hưởng của nhiệt độ môi trường và tải trọng hãm của đoàn xe nhiều trục.

Quý độc giả, kỹ sư và các nhà nghiên cứu quan tâm có thể ứng dụng mã nguồn thuật toán và các kết quả số trong công trình này để phục vụ công tác thiết kế, thẩm tra và nâng cấp chất lượng các dự án đường băng sân bay và đường cao tốc hiện đại.