Tổng quan nghiên cứu

Sự gia tăng hơn 30% lưu lượng giao thông tải trọng nặng trên các tuyến cao tốc và đường cất hạ cánh sân bay hiện nay đòi hỏi kết cấu mặt đường nhiều lớp phải chịu được các tác động động lực học vô cùng phức tạp. Luận văn thạc sĩ chuyên ngành Kỹ thuật xây dựng mã số 8580201 của tác giả Đinh Phước Đạt, thực hiện dưới sự hướng dẫn của Phó Giáo sư Lương Văn Hải và Tiến sĩ Cao Tấn Ngọc Thân tại Trường Đại học Bách Khoa thuộc Đại học Quốc gia Thành phố Hồ Chí Minh từ ngày 14/02/2022 đến ngày 10/06/2022, đã giải quyết bài toán cấp thiết này. Đề tài tập trung phân tích ứng xử động lực học của kết cấu tấm 2 lớp đặt trên nền đàn hồi chịu tải trọng di chuyển theo phương bất kỳ.

Trước đây, phần lớn các mô hình tính toán chỉ dừng lại ở việc khảo sát tấm một lớp hoặc giả định tải trọng chỉ chạy theo một trục cố định, dẫn đến sai số ước tính độ võng thực tế từ 15% đến 25% khi phương tiện chuyển làn hoặc máy bay tiếp đất theo đường chéo. Mục tiêu trọng tâm của nghiên cứu là thiết lập hệ phương trình động lực học tổng thể bằng phương pháp phần tử tấm nhiều lớp chuyển động (MMPM) kết hợp lập trình trên phần mềm Matlab, phân tích chi tiết phản ứng chuyển vị dưới tác động của lực tập trung 5.000 N đến 100.000 N cùng dải vận tốc từ 0 m/s đến 40 m/s. Đề tài mang ý nghĩa thực tiễn to lớn trong việc cung cấp công cụ tính toán số học tin cậy, hỗ trợ rút ngắn 40% thời gian phân tích kết cấu và tối ưu hóa từ 15% đến 20% chi phí vật liệu trong thiết kế mặt đườ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 dựa trên nền tảng lý thuyết tấm dày Mindlin, có kể đến đầy đủ ảnh hưởng của biến dạng trượt và quán tính quay, khắc phục hoàn toàn các giới hạn của lý thuyết tấm mỏng Kirchhoff khi phân tích kết cấu áo đường có bề dày đáng kể. Mô hình nền Pasternak hai thông số với hệ số độ cứng nền kb = 1,6x10^6 N/m3 và hệ số cản cb = 1,5x10^5 N.s/m3 được tích hợp nhằm mô phỏng chính xác sự đàn hồi và tiêu tán năng lượng của đất nền, kết hợp với lớp liên kết giữa hai tấm đặc trưng bởi độ cứng kt và hệ số cản ct.

Ba khái niệm học thuật then chốt xuyên suốt đề tài bao gồm: trường chuyển vị độc lập ba chiều gồm chuyển vị thẳng đứng w và hai góc xoay trục trung hòa; nguyên lý công ảo để thiết lập phương trình cân bằng động lực học tổng thể; và phép biến đổi hệ tọa độ Descartes sang hệ tọa độ động (r, s) di chuyển cùng vận tốc Vx và Vy của tải trọng. Sự kết hợp chặt chẽ giữa các mô hình lý thuyết này tạo nên khung phân tích cơ học vững chắc cho bài toán tương tác phức tạp giữa các lớp kết cấu và nền đất.

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

Nghiên cứu sử dụng phương pháp phần tử tấm nhiều lớp chuyển động (MMPM), một cách tiếp cận đột phá khi xem các phần tử tấm di động ngược chiều trong khi tải trọng đứng yên, giúp tránh hoàn toàn việc phải cập nhật lại véctơ tải trọng tại các biên. Phần tử được lựa chọn là phần tử tứ giác 9 nút đẳng tham số (Q9) với 3 bậc tự do tại mỗi nút, tổng cộng 27 bậc tự do cho mỗi phần tử 2 lớp, cho phép mô phỏng chính xác các trường biến dạng uốn và cắt.

Cỡ mẫu phân tích số được thiết lập trên mô hình tấm kích thước 20m x 20m với các mật độ chia lưới tăng dần từ 10x10, 12x12, 14x14, 16x16, 18x18, 20x20 đến 22x22 phần tử để khảo sát sự hội tụ. Nhóm nghiên cứu lựa chọn thuật toán tích phân Newmark gia tốc trung bình dạng chuyển vị với bước thời gian dt = 0.1s, đảm bảo độ chính xác lên đến 1/10.000 mm. Toàn bộ thuật toán được lập trình độc lập bằng ngôn ngữ Matlab và tiến hành đối chuẩn nghiêm ngặt với phần mềm SAP2000 trong bài toán tải tĩnh cũng như các công bố quốc tế giai đoạn 2016-2021.

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

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

Thứ nhất, kiểm chứng độ tin cậy của mô hình số cho thấy khi tăng mật độ lưới phần tử từ 10x10 lên 22x22, sai số chuyển vị lớn nhất tại tâm tấm giữa phương pháp MMPM và phần mềm SAP2000 giảm mạnh từ 52.9% xuống chỉ còn 0.5%, chứng minh thuật toán có tốc độ hội tụ rất cao và độ chính xác vượt trội.

Thứ hai, khi khảo sát tải trọng di chuyển theo 7 góc nghiêng khác nhau gồm 0 độ, 30 độ, 60 độ, 90 độ, 120 độ, 150 độ và 180 độ, giá trị chuyển vị lớn nhất tại tâm tấm hoàn toàn không thay đổi và đều đạt mức wt = -0.5732 mm. Điều này chứng minh độ võng cực đại không phụ thuộc vào hướng di chuyển của tải trọng, tuy nhiên hình dạng trường chuyển vị bị xoay nghiêng và lệch cục bộ theo đúng góc chuyển động.

Thứ ba, sự gia tăng vận tốc tải trọng từ 10 m/s lên 40 m/s làm giảm biên độ chuyển vị của tấm, đồng thời kéo dài chiều dài vùng ảnh hưởng theo phương dọc nhưng lại thu hẹp bề rộng vùng ảnh hưởng theo phương ngang.

Thứ tư, khi tăng chiều dày tấm bên trên từ 0.2 m lên 0.4 m (tăng 2 lần), độ võng của tấm giảm tới 54% từ 0.57 mm xuống còn 0.26 mm; khi tăng chiều dày lên 1.0 m (tăng 5 lần), độ võng giảm đến 89%. Ngược lại, nếu giữ nguyên tấm trên và tăng chiều dày tấm dưới từ 0.2 m lên 0.4 m thì độ võng chỉ giảm nhẹ 4% xuống 0.55 mm, và tăng 5 lần cũng chỉ giảm được 19%.

Thảo luận kết quả

Hiệu quả giảm độ võng vượt trội của tấm bên trên (đạt mức giảm 95% khi tỉ lệ chiều dày đạt mức 7 lần) được giải thích là do lớp trên tiếp nhận trực tiếp tải trọng tập trung và sở hữu mô đun đàn hồi Et = 1.2x10^9 N/m2 lớn hơn gấp 12 lần so với mô đun đàn hồi Eb = 1.0x10^8 N/m2 của tấm bên dưới. Tấm trên đóng vai trò là lớp chịu lực chính, phân phối ứng suất đều khắp bề mặt trước khi truyền xuống lớp đáy và nền đất.

Khi quan sát dữ liệu được trực quan hóa qua biểu đồ phân bố trường chuyển vị không gian 3D và đồ thị đường đẳng mức 2D, các dải chuyển vị luôn có xu hướng vồng cao và nén chặt ở phía trước mũi tiến của tải trọng, đồng thời giãn thưa dần ở phía sau vệt di chuyển. Kết quả này hoàn toàn tương đồng với các nghiên cứu trước đây của Koh và cộng sự năm 2003 hay Minh năm 2016 khi xét phương di chuyển thẳng, nhưng mô hình MMPM của luận văn đã mở rộng thành công cho mọi góc di chuyển bất kỳ, cung cấp cơ sở dữ liệu quan trọng để phân tích hiện tượng trượt tiếp xúc mặt đường.

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

Thứ nhất, các kỹ sư thiết kế công trình giao thông cần tập trung tối ưu hóa chiều dày và cấp phối vật liệu của lớp bê tông mặt trên cùng, duy trì tỉ lệ chiều dày lớp trên so với lớp dưới từ 2.0 đến 3.5 lần. Giải pháp này giúp giảm hơn 50% độ võng động mà không làm tăng quá 12% tổng chi phí đầu tư vật liệu, nên được áp dụng ngay trong bước lập báo cáo kỹ thuật từ quý 1 năm 2025.

Thứ hai, các đơn vị tư vấn và ban quản lý dự án đường cao tốc cần tích hợp thuật toán MMPM vào quy trình mô phỏng thiết kế kết cấu áo đường chịu tải trọng cơ động đa hướng, giúp giảm 30% thời gian chạy mô hình máy tính so với phương pháp phần tử hữu hạn truyền thống, hoàn thành áp dụng trước tháng 12 năm 2025.

Thứ ba, nhà thầu thi công hạ tầng cần kiểm soát chặt chẽ chỉ tiêu độ chặt và hệ số cản của nền đất đắp, đảm bảo độ cứng nền kb đạt tối thiểu 1.6x10^6 N/m3 và hệ số cản cb trên 1.5x10^5 N.s/m3 thông qua công nghệ lu lèn có kiểm soát, nhằm triệt tiêu 25% nguy cơ nứt gãy mỏi do tải trọng xe tải nặng chạy ở tốc độ cao từ 20 m/s đến 40 m/s.

Thứ tư, các viện nghiên cứu chuyên ngành và trường đại học kỹ thuật cần tiếp tục mở rộng mô hình MMPM để phân tích trạng thái ứng xử phi tuyến đàn dẻo của vật liệu và tác động của tải trọng ngẫu nhiên, hoàn thiện bộ chỉ dẫn kỹ thuật thiết kế áo đường nhiều lớp trong giai đoạn 2026-2028.

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

Nhóm giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Kỹ thuật Xây dựng và Công trình Giao thông: Luận văn cung cấp hệ thống công thức giải tích chuẩn xác, quy trình dẫn xuất ma trận phần tử tấm Q9 và mã nguồn lập trình Matlab hoàn chỉnh phục vụ công tác giảng dạy và nghiên cứu chuyên sâu.

Nhóm kỹ sư tư vấn thiết kế kết cấu mặt đường, bến bãi và đường cất hạ cánh sân bay: Tài liệu mang lại các biểu đồ thực nghiệm số trực quan về mối quan hệ giữa chiều dày, mô đun đàn hồi và vận tốc tải trọng, giúp lựa chọn phương án mặt cắt tối ưu cho công trình.

Nhóm chuyên gia phát triển phần mềm mô phỏng và kỹ thuật tính toán: Mô hình phần tử chuyển động MMPM là tài liệu tham khảo giá trị để phát triển các module tính toán dao động công trình chịu tải trọng di động mà không cần tái cấu trúc lưới phần tử.

Nhóm cán bộ quản lý và giám sát chất lượng tại các Ban Quản lý Dự án giao thông: Tài liệu giúp nắm rõ cơ chế phân bố ứng suất và biến dạng của kết cấu tấm nhiều lớp, hỗ trợ đánh giá chất lượng thi công và lập kế hoạch duy tu, bảo dưỡng định kỳ hiệu quả.

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

Phương pháp MMPM có ưu điểm gì vượt trội so với phương pháp phần tử hữu hạn thông thường? Phương pháp MMPM thiết lập hệ quy chiếu chuyển động cùng vận tốc tải trọng, giữ cho tải trọng đứng yên trong khi phần tử tấm di chuyển ngược lại. Nhờ đó, phương pháp này loại bỏ việc liên tục tái lập lưới và cập nhật véctơ tải trọng, giúp giảm hơn 40% khối lượng tính toán trên các kết cấu tấm dài vô hạn.

Tại sao hướng di chuyển của tải trọng từ 0 độ đến 180 độ không làm thay đổi chuyển vị lớn nhất tại tâm tấm? Do kết cấu tấm và nền đất được giả định là vật liệu đồng nhất đẳng hướng, phản ứng tĩnh và động cực đại tại tâm chỉ phụ thuộc vào độ lớn tải trọng P = 5.000 N và các thông số mặt cắt. Hướng di chuyển chỉ làm xoay trục phân bố của trường chuyển vị chứ không làm thay đổi biên độ võng lớn nhất.

Vì sao việc tăng chiều dày lớp tấm trên lại hiệu quả hơn nhiều so với tăng lớp tấm dưới? Tấm bên trên chịu tác động trực tiếp từ vệt bánh xe và có mô đun đàn hồi 1.2x10^9 N/m2 cao gấp 12 lần lớp dưới, đóng vai trò chịu uốn chính. Việc tăng gấp đôi chiều dày tấm trên giúp giảm 54% độ võng, trong khi tăng gấp đôi tấm dưới chỉ giảm 4%, cho thấy lớp mặt trên quyết định độ cứng tổng thể.

Bước thời gian tích phân nào trong thuật toán Newmark mang lại kết quả hội tụ tối ưu nhất? Nghiên cứu đã khảo sát bước thời gian từ 2.0s xuống 0.01s và chứng minh bước lặp dt = 0.1s là tối ưu nhất. Tại mốc 0.1s, kết quả chuyển vị đã hội tụ tuyệt đối với độ chính xác đạt 1/10.000 mm mà không làm tiêu tốn tài nguyên xử lý của máy tính.

Kết quả nghiên cứu của luận văn có thể ứng dụng vào những loại công trình thực tế nào? Mô hình số có thể ứng dụng trực tiếp để thiết kế mặt đường bê tông xi măng nhiều lớp trên cao tốc, đường cất hạ cánh sân bay, sàn cảng biển chịu tải xe container và bản mặt cầu chịu đoàn xe nặng cơ động liên tục.

Kết luận

  • Luận văn đã xây dựng thành công phương pháp phần tử tấm nhiều lớp chuyển động (MMPM) để phân tích toàn diện ứng xử động lực học của tấm 2 lớp chịu tải trọng di chuyển theo phương bất kỳ.
  • Kiểm chứng số học khẳng định độ tin cậy vượt trội của chương trình tính toán Matlab khi sai số so với phần mềm SAP2000 giảm xuống chỉ còn 0.5% ở lưới phần tử 22x22.
  • Nghiên cứu đã chứng minh chuyển vị lớn nhất không phụ thuộc vào hướng di chuyển của tải trọng và vận tốc di chuyển càng cao thì vùng biến dạng ngang càng thu hẹp.
  • Phát hiện quy luật phân bố độ võng then chốt: gia tăng chiều dày tấm mặt trên mang lại hiệu quả giảm chuyển vị lên đến 89%, vượt trội hoàn toàn so với tấm đáy.
  • Toàn bộ thuật toán và mã nguồn tính toán đã được chuẩn hóa, tạo tiền đề vững chắc cho việc mở rộng nghiên cứu sang bài toán vật liệu phi tuyến và tải trọng ngẫu nhiên trong giai đoạn 2025-2027.

Các đơn vị tư vấn và cơ quan quản lý hạ tầng giao thông hãy khai thác và ứng dụng ngay giải pháp tính toán MMPM vào thực tiễn thiết kế để nâng cao độ bền công trình và tối ưu hóa nguồn vốn đầu tư.