Tổng quan về luận án

Hiện tượng hư hỏng hằn lún vệt bánh xe (rutting) trên mặt đường bê tông asphalt (BTA) dưới tác động đồng thời của tải trọng trục xe nặng, áp lực bánh xe cao và nền nhiệt độ môi trường khắc nghiệt là một thách thức kỹ thuật cốt lõi trong xây dựng công trình giao thông tại Việt Nam cũng như trên toàn cầu. Luận án tiến sĩ kỹ thuật với đề tài "Nghiên cứu quan hệ giữa cường độ chống cắt và đặc tính biến dạng không hồi phục của hỗn hợp bê tông Asphalt mặt đường" do nghiên cứu sinh Ngô Ngọc Quý thực hiện dưới sự hướng dẫn khoa học của PGS.TS Trần Thị Kim Đăng và PGS.TS Lã Văn Chăm (Trường Đại học Giao thông Vận tải, 2021, Chuyên ngành Xây dựng đường ôtô và đường thành phố, Mã số: 9.05) đã tiên phong giải quyết bài toán cơ học vật liệu mặt đường từ bản chất ứng suất - biến dạng dẻo trượt cắt.

                  +----------------------------------------------+
                  | Tải trọng trục nặng & Nhiệt độ khai thác cao |
                  +----------------------------------------------+
                  +----------------------------------------------+
                  |  Ứng suất cắt vượt ngưỡng cường độ chống cắt |
                  |            (Mohr-Coulomb: τ > τf)            |
                  +----------------------------------------------+
                  +----------------------------------------------+
                  |    Biến dạng trượt cắt / Đẩy trồi ngang      |
                  |     (Chiếm tới ~80% chiều sâu lún vệt bánh)  |
                  +----------------------------------------------+
                  +----------------------------------------------+
                  |  Tích lũy biến dạng dẻo không hồi phục (εp)  |
                  |   Mô hình hóa qua tỷ số ứng suất cắt (τ/τf)  |
                  +----------------------------------------------+

Về mặt lịch sử khoa học, phần lớn các nghiên cứu trong nước trước đây chủ yếu tập trung vào hướng tiếp cận thực nghiệm đơn thuần: cải thiện thành phần cấp phối cốt liệu, kiểm soát độ ổn định Marshall, hoặc áp dụng các loại nhựa đường biến tính polymer mà chưa đi sâu vào cơ chế phá hoại trượt cắt nội tại của khung cốt liệu và chất kết dính. Khoảng trống nghiên cứu (research gap) cốt lõi được nhận diện là sự thiếu vắng một khung lý thuyết định lượng và mô hình thực nghiệm xác lập mối liên hệ toán học giữa các tham số cường độ chống cắt ($c, \phi$) với tốc độ tích lũy biến dạng không hồi phục ($\varepsilon_p$) của các hỗn hợp bê tông asphalt chặt phổ biến tại Việt Nam.

Luận án tập trung giải quyết 4 câu hỏi nghiên cứu (Research Questions - RQ) và kiểm chứng 4 giả thuyết khoa học (Hypotheses - H) tương ứng:

  • RQ1 & H1: Bản chất cơ học của hư hỏng hằn lún vệt bánh xe chịu sự chi phối quyết định bởi biến dạng trượt cắt nội tại của hỗn hợp BTA khi tỷ số ứng suất cắt vượt ngưỡng chịu tải giới hạn, chứ không đơn thuần chỉ do quá trình nén chặt thứ cấp (post-compaction).
  • RQ2 & H2: Việc cải tiến hệ thống thiết bị thí nghiệm nén ba trục (Triaxial Compression Test) và nén ba trục tải trọng lặp (Repeated Load Permanent Deformation - RLPD) chuyên dụng cho BTA cho phép lượng hóa chính xác lực dính đơn vị ($c$), góc nội ma sát ($\phi$) và động học tích lũy biến dạng dẻo theo thời gian và số chu kỳ tải.
  • RQ3 & H3: Các tham số sức kháng cắt ($c, \phi$) chịu tác động phi tuyến mạnh mẽ bởi nhiệt độ khai thác ($30^\circ\text{C}, 45^\circ\text{C}, 60^\circ\text{C}$), tốc độ gia tải ($0.5, 1.0, 7.5\text{ mm/phút}$), chủng loại chất kết dính (nhựa 60/70, 40/50, PMB III) và nguồn gốc xuất xứ cốt liệu khoáng.
  • RQ4 & H4: Tồn tại một tương quan thực nghiệm có độ tin cậy thống kê cao ($R^2 \ge 0.85$) giữa biến dạng không hồi phục ($\varepsilon_p$) với tỷ số ứng suất cắt trên cường độ chống cắt ($\tau/\tau_f$) và số chu kỳ gia tải ($N$) ở nhiệt độ tính toán bất lợi $60^\circ\text{C}$.

Khung lý thuyết của luận án tích hợp thuyết phá hoại trượt cắt Mohr-Coulomb, nguyên lý biến dạng đàn - nhớt - dẻo (Visco-elasto-plasticity) từ Chương trình Nghiên cứu Chiến lược Đường bộ Hoa Kỳ (SHRP) và phương pháp thiết kế kết cấu mặt đường theo Cơ học - Thực nghiệm (Mechanistic-Empirical Pavement Design - M-E). Đóng góp đột phá của công trình được lượng hóa qua việc hoàn thiện 02 hệ thống thiết bị đo đạc ba trục trong phòng thí nghiệm, thiết lập bộ cơ sở dữ liệu thực nghiệm chuẩn tắc về đặc tính kháng cắt của bê tông asphalt chặt BTAC12.5 và BTAC19, đồng thời xây dựng các phương trình dự báo biến dạng không hồi phục theo hàm số của $\tau/\tau_f$. Phạm vi nghiên cứu bao quát các tổ hợp mẫu chế tạo từ 2 nguồn mỏ đá điển hình (mỏ Đồng Ao - Hà Nam và mỏ Tân Đông Hiệp - Bình Dương), 3 loại bitum ở các cấp nhiệt độ từ $30^\circ\text{C}$ đến $60^\circ\text{C}$, tạo nền tảng chuyển dịch phương pháp thiết kế áo đường từ cơ chế ứng suất kéo uốn truyền thống sang kiểm soát đồng thời biến dạng trượt dẻo.

Literature Review và Positioning

Các công trình nghiên cứu quốc tế và trong nước về hành vi cơ học của bê tông asphalt được phân hóa thành ba dòng tư tưởng (research streams) chính:

+---------------------------------------------------------------------------------------+
|                                     DÒNG TƯ TƯỞNG 1                                   |
|   Cơ chế hình thành hằn lún vệt bánh & Trượt cắt nội tại                              |
|   - Monismith et al. (1994), Sousa et al. (1994) [SHRP-A-003A]                        |
|   - Zheng et al. (2012) [Tây Nam - Trung Quốc]: Cắt trượt chiếm 80% độ lún             |
+---------------------------------------------------------------------------------------+
|                                     DÒNG TƯ TƯỞNG 2                                   |
|   Mô hình hóa sức kháng cắt (c, ϕ) & Thí nghiệm 3 trục                                |
|   - Tan et al. (1994) [NUS - Singapore]: Khung c-ϕ & Giới hạn Smith                    |
|   - Bi & Zheng (2014), Christensen & Bonaquist (2004) [NCAT/NCHRP]                    |
+---------------------------------------------------------------------------------------+
|                                     DÒNG TƯ TƯỞNG 3                                   |
|   Tương quan tỷ số ứng suất cắt (τ/τf) & Biến dạng không hồi phục (εp)               |
|   - Baek et al. (2008), Kim et al. (2009, 2014) [Đại học Sejong - Hàn Quốc]           |
|   - GAP TẠI VIỆT NAM: Chưa có mô hình chuẩn hóa cho cấp phối BTAC nhiệt đới           |
+---------------------------------------------------------------------------------------+
  1. Dòng nghiên cứu cơ chế phá hoại lún vệt bánh và biến dạng cắt: Monismith và cộng sự (1994), Sousa và cộng sự (1994) trong khuôn khổ dự án SHRP-A-003A đã chỉ rõ biến dạng trượt cắt không hồi phục (shear flow) dưới mép bánh xe mới là nguyên nhân áp đảo gây lún trồi mặt đường, bác bỏ quan điểm truyền thống cho rằng lún chỉ do nén chặt cố kết. Nghiên cứu thực nghiệm hiện trường của Zheng và cộng sự (2012) tại Đại học Giao thông Tây Nam (Trung Quốc) khẳng định: "Biến dạng do cắt chiếm đến 80% tổng chiều sâu lún" trên các tuyến đường cao tốc chịu tải nặng.
  2. Dòng nghiên cứu tham số kháng cắt ($c, \phi$) và mô hình thí nghiệm: Tan và cộng sự (1994) tại Đại học Quốc gia Singapore (NUS) đã phát triển mô hình giải tích kết hợp tham số lực dính đơn vị ($c$), góc nội ma sát ($\phi$) với phân tích phần tử hữu hạn (PTHH) để xây dựng "đường giới hạn Smith" nhằm tích hợp thiết kế thành phần hỗn hợp và đánh giá kết cấu. Tiếp đó, các nghiên cứu của Bi và Zheng (2014), Christensen và Bonaquist (2004) qua các báo cáo NCHRP Report 465 và NCAT Report 03-07 nhấn mạnh tính ưu việt của mô hình nén ba trục có áp lực hông ($\sigma_3$) so với thí nghiệm Marshall hay kéo gián tiếp (IDT), do mô phỏng chính xác trạng thái ứng suất khối 3 chiều trong tầng mặt đường thực tế.
  3. Dòng nghiên cứu tương quan giữa cường độ kháng cắt và biến dạng dẻo lặp: Nhóm tác giả Baek, Kim và cộng sự (2008, 2009, 2014) tại Đại học Sejong (Hàn Quốc) đã đặt nền móng cho việc sử dụng tỷ số ứng suất cắt trên cường độ chống cắt ($\tau/\tau_f$) như một biến số không thứ nguyên cốt lõi để dự báo biến dạng tích lũy $\varepsilon_p$ qua mô hình toán: $$\varepsilon_p = A \cdot N^B \cdot \left(\frac{\tau}{\tau_f}\right)^C \cdot t^D$$ Mô hình này đạt hệ số tương quan $R^2 = 0.9969$ trong phòng thí nghiệm và được kiểm chứng thành công trên thiết bị thử nghiệm gia tốc phá hoại mặt đường (Accelerated Pavement Tester - APT dài 12m) và 26 đoạn thử nghiệm hiện trường tại Hàn Quốc với hệ số tương quan thực tế đạt $R^2 = 0.87$, sai số căn bậc hai trung bình sai số hồi quy $MRSE = 1.13\text{ mm}$.

Tại Việt Nam, các nghiên cứu của Nguyễn Quang Phúc (2010), Bùi Tuấn Anh và nhóm tác giả Đại học Xây dựng (2014, 2016) mới tiếp cận các mô hình cắt phẳng trực tiếp, nén chẻ hoặc thiết bị hằn lún Hamburg Wheel-Tracking theo Quyết định 858/QĐ-BGTVT và TCVN 8819:2011, mà chưa thực hiện được các thí nghiệm nén ba trục động (RLPD) đồng bộ.

Luận án của Ngô Ngọc Quý định vị chính xác tại điểm giao thoa giữa cơ học nứt trượt hiện đại và kỹ thuật thực nghiệm giao thông: mở rộng khung phân tích của Baek et al. và Zheng et al., bản địa hóa vào hệ thống vật liệu cấp phối chặt nhiệt đới (BTAC) của Việt Nam, lấp đầy khoảng trống dữ liệu cơ lý về sức kháng cắt ba trục ở nhiệt độ cao $60^\circ\text{C}$.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng đường biên ứng dụng của Thuyết phá hoại cắt Mohr-Coulomb cổ điển sang môi trường vật liệu hỗn hợp hạt nhớt - dẻo phức tạp là bê tông asphalt mặt đường:

$$\tau_f = c(T, \dot{\varepsilon}, \text{Bitumen}) + \sigma_n \cdot \tan\phi(\text{Aggregate}, \text{Gradation})$$

Yếu tố tương tác Cơ chế vi mô (Micro-mechanisms) Đóng góp lý thuyết cụ thể
Lực dính đơn vị ($c$) Tương tác giữa màng bitum - bột khoáng (asphalt mastic), liên kết phân tử và độ quánh lưu biến Giảm mạnh theo hàm mũ khi nhiệt độ tăng từ $30^\circ\text{C}$ lên $60^\circ\text{C}$; tăng tới $32%$ khi sử dụng bitum PMB III có độ nhớt cao
Góc nội ma sát ($\phi$) Sự chèn móc cơ học (interlocking), độ góc cạnh và ma sát tiếp xúc giữa các hạt cốt liệu thô Khá ổn định ở dải nhiệt độ $30^\circ\text{C} - 45^\circ\text{C}$, chỉ suy giảm nhẹ ở $60^\circ\text{C}$ do tác động bôi trơn của màng nhựa chảy dẻo

Nghiên cứu khẳng định sự chuyển dịch hệ hình (paradigm shift) từ tư duy xem BTA là vật liệu đàn hồi tuyến tính thuần túy sang vật liệu nhớt - dẻo phụ thuộc trạng thái ứng suất cắt giới hạn. Luận án đã xác lập một hệ tiên đề cơ học: Khả năng kháng lún vệt bánh xe của tầng mặt đường không phụ thuộc cục bộ vào độ ổn định Marshall hay mô đun đàn hồi tĩnh, mà được chi phối trực tiếp bởi tỷ số ứng suất cắt làm việc so với sức chống cắt giới hạn ($\tau/\tau_f$).

                    KHUNG KHÁI NIỆM & MÔ HÌNH LÝ THUYẾT
                    
   [ Thuyết Mohr-Coulomb ]      [ Nguyên lý SHRP A-003A ]      [ Cơ học - Thực nghiệm (M-E) ]
   (c: Lực dính vữa nhựa)        (Biến dạng trượt dẻo        (Mô hình hóa hàm số εp theo N
   (ϕ: Chèn móc cốt liệu)         chiếm ưu thế > 80%)         và mức tải trọng giới hạn)
                       +-----------------------------------------+
                       | KHUNG PHÂN TÍCH ĐỘC ĐÁO CỦA LUẬN ÁN     |
                       | Mô hình hóa tương quan phi tuyến        |
                       |       εp = f(τ/τf, N, T=60°C)           |
                       +-----------------------------------------+
  [ Mức tải thấp: τ/τf ≤ 0.269 ]                            [ Mức tải cao: τ/τf ≥ 0.425 ]
  Mô hình lũy thừa ổn định (Power law)                      Mô hình đa thức/chuyển dịch dẻo
  Biến dạng tiệm cận giai đoạn 2 (ổn định)                  Xuất hiện điểm uốn chảy dẻo (FN)

Khung phân tích độc đáo

Khung phân tích độc đáo của luận án là sự hợp nhất 3 lý thuyết nền tảng: (1) Thuyết bền Mohr-Coulomb; (2) Nguyên lý tích lũy biến dạng từ biến dẻo (Creep & Permanent Strain Theory); và (3) Lý thuyết thiết kế Cơ học - Thực nghiệm (M-E Design).

Khung phân tích đề xuất phân tách rõ ràng hai miền ứng xử biến dạng của BTA:

  • Miền tải trọng thấp ($\tau/\tau_f \le 0.269$): Biến dạng không hồi phục tích lũy theo quy luật lũy thừa suy giảm dần tốc độ biến dạng, vật liệu làm việc trong giai đoạn ổn định (secondary stage), cấu trúc khung đá giữ vai trò chịu lực chính.
  • Miền tải trọng cao ($\tau/\tau_f \ge 0.425$): Tốc độ tích lũy biến dạng tăng tốc nhanh chóng, vật liệu nhanh chóng vượt qua giai đoạn ổn định để tiến vào giai đoạn phá hoại chảy dẻo (tertiary stage), xuất hiện biến dạng trượt trồi biên rõ rệt.

Điều kiện biên (boundary conditions) của mô hình được xác lập chặt chẽ: áp dụng cho hỗn hợp BTA rải nóng cấp phối chặt danh định 12.5mm và 19mm, độ rỗng dư thiết kế $VA = 3 - 5%$, nhiệt độ khai thác quy chuẩn $60^\circ\text{C}$ và áp lực hông $\sigma_3$ dao động từ $0$ đến $200\text{ kPa}$.

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo quan điểm thực chứng (positivism) với phương pháp tiếp cận thực nghiệm - định lượng đa cấp độ (multi-level experimental design).

                                THIẾT KẾ NGHIÊN CỨU ĐA CẤP ĐỘ
                                
       +-----------------------------------------------------------------------------+
       |                           1. THIẾT KẾ CẤP PHỐI & VẬT LIỆU                   |
       |  - 2 loại cấp phối: BTAC12.5 và BTAC19                                      |
       |  - 2 nguồn mỏ đá: Mỏ Đồng Ao (CP1) & Tân Đông Hiệp (CP2)                     |
       |  - 3 loại nhựa: 60/70 (N1), 40/50 (N2), Polyme PMB III (N3)                 |
       |  - Tối ưu hóa hàm lượng nhựa bằng phương pháp Marshall tiêu chuẩn           |
       +-------------------------------+               +-------------------------------+
       |   2. ĐẶC TÍNH KHÁNG CẮT (c, ϕ)|               |   3. ĐẶC TÍNH BIẾN DẠNG DẺO   |
       | - Nén 3 trục tĩnh (Triaxial)  |               | - Thí nghiệm từ biến tĩnh SCT |
       | - Nhiệt độ: 30, 45, 60°C      |               | - Nén 3 trục lặp RLPD (60°C)  |
       | - Tốc độ: 0.5, 1.0, 7.5 mm/ph |               | - Chu kỳ tải: N đến 40.000    |
       | - Áp lực hông: 0 - 200 kPa    |               | - Mức ứng suất: τ/τf thay đổi |
       +-----------------------------------------------------------------------------+
       |                      4. PHÂN TÍCH DỮ LIỆU & TƯƠNG QUAN MÔ HÌNH              |
       |  - Phân tích phương sai ANOVA & Kiểm định hậu nghiệm Tukey HSD (p < 0.05)   |
       |  - Thiết lập phương trình hồi quy phi tuyến: εp = f(τ/τf, N)                |
       |  - Kiểm chuẩn mô hình toán thống kê (R², Se/Sy)                             |
       +-----------------------------------------------------------------------------+

Hệ thống mẫu thử nghiệm gồm hàng trăm mẫu hình trụ chuẩn ($D = 100\text{ mm}, H = 150\text{ mm}$ và $D = 101.6\text{ mm}, H = 63.5\text{ mm}$) được đúc ép bằng thiết bị đầm Marshall và đầm xoay (Superpave Gyratory Compactor), kiểm soát chặt chẽ các chỉ tiêu thể tích: độ rỗng dư ($VA$), độ rỗng cốt liệu ($VMA$), độ rỗng lấp đầy nhựa ($VFA$).

Quy trình nghiên cứu rigorous

Tác giả đã cải tiến thành công 02 cụm thiết bị thí nghiệm chuyên sâu tại Phòng thí nghiệm Công trình Vilas 047 và Trung tâm KHCN GTVT (Trường Đại học GTVT):

+---------------------------------------------------------------------------------------------+
|                             HỆ THỐNG THÍ NGHIỆM ĐƯỢC CẢI TIẾN TẠI VILAS 047                 |
+--------------------------------------------------------------+------------------------------+
| 1. Buồng nén ba trục BTA cải tiến                            | 2. Thiết bị nén tải lặp RLPD|
| - Cải tiến từ buồng nén đất tiêu chuẩn                      | - Gia tải dọc trục lặp chu kỳ|
| - Khung gia tải cứng vững chịu tải trọng lớn                 | - Áp lực khí nén tạo áp lực   |
| - Hệ thống gia nhiệt tuần hoàn nước/khí (±0.5°C)             |   hông σ3 ổn định            |
| - Màng bao mẫu cao su chịu nhiệt, đánh giá ảnh hưởng truyền lực| - Đầu đo LVDT/CDP-25 độ nhạy |
| - Load cell điện tử & chuyển vị điện tử CDP-25              |   cao kết nối ADC            |
+--------------------------------------------------------------+------------------------------+

Quy trình bảo đảm tính giá trị (validity) và độ tin cậy (reliability) của số liệu bao gồm:

  • Hiệu chuẩn cảm biến: Toàn bộ cảm biến lực (Load cell) và cảm biến biến dạng vi sai (LVDT/CDP-25) được hiệu chuẩn tuyến tính với sai số $< 0.1%$.
  • Đánh giá ảnh hưởng màng bao mẫu: Tiến hành thí nghiệm đối chứng có và không có màng cao su bao mẫu ở các mức áp lực hông khác nhau, chứng minh ảnh hưởng của màng bao đến giá trị ứng suất chính $\sigma_1$ và biến dạng $\varepsilon$ là không có ý nghĩa thống kê ở mức ý nghĩa $\alpha = 0.05$.
  • Triangulation phương pháp: Kết hợp thí nghiệm từ biến tĩnh (Static Creep Test - SCT) và thí nghiệm nén dọc trục tải trọng lặp (RLPD) để kiểm chứng chéo các pha biến dạng của vật liệu.

Data và phân tích

Dữ liệu thực nghiệm được xử lý bằng các công cụ thống kê nâng cao:

  • Phân tích phương sai ANOVA hai nhân tố (Two-way ANOVA) nhằm phân lập ảnh hưởng độc lập và tương tác của cấp phối, loại chất kết dính, nhiệt độ và tốc độ nén đến các biến phản hồi $c$ và $\phi$.
  • Kiểm định so sánh cặp Tukey HSD (Honestly Significant Difference) với khoảng tin cậy $95%$ ($p\text{-value} < 0.05$) để khẳng định sự khác biệt bản chất giữa các tổ hợp nhựa đường và cấp phối.
  • Đánh giá mức độ thích hợp của phương trình hồi quy thông qua hệ số tương quan $R^2$, sai số tiêu chuẩn hồi quy ($S_e$) và tỷ số $S_e/S_y$.

Phát hiện đột phá và implications

Những phát hiện then chốt

+----------------------------------------------------------------------------------------------+
|                                    5 PHÁT HIỆN ĐỘT PHÁ CỦA LUẬN ÁN                          |
+----------------------------------------------------------------------------------------------+
| [1] Nhựa đường PMB III tăng 32% lực dính c; Khung cốt liệu thô quyết định góc ma sát ϕ.      |
| [2] Nhiệt độ tăng (30°C -> 60°C) làm c suy giảm tới 60-70%; ϕ duy trì độ ổn định tương đối.   |
| [3] Tồn tại ngưỡng phân dị ứng xử tải trọng tại τ/τf ≈ 0.269 - 0.425.                        |
| [4] Tốc độ gia tải tăng (0.5 -> 7.5 mm/phút) làm tăng đồng thời cả c và ϕ.                   |
| [5] Mô hình hóa thành công tương quan εp = f(τ/τf, N) với R² = 0.88 - 0.96.                  |
+----------------------------------------------------------------------------------------------+
  1. Sự phân hóa vai trò cơ học giữa chất kết dính và khung cốt liệu: Lực dính đơn vị $c$ của hỗn hợp BTA phụ thuộc chủ yếu vào bản chất mác nhựa và hàm lượng mastic nhựa đường; trong khi góc nội ma sát $\phi$ hoàn toàn do cấu trúc chèn móc của khung cốt liệu thô chi phối. Thí nghiệm khẳng định: "Lực dính đơn vị của hỗn hợp chặt AC13 và hỗn hợp thiết kế theo Superpave lớn hơn so với hỗn hợp có cấu trúc khung chặt SMA... nhựa đường cải tiến với độ nhớt cao làm lực dính tăng đến 32%". Cốt liệu mỏ Tân Đông Hiệp (CP2) có độ góc cạnh và độ nhám bề mặt cao hơn cốt liệu mỏ Đồng Ao (CP1), đem lại góc nội ma sát $\phi$ vượt trội ($41.5^\circ - 44.2^\circ$ so với $38.7^\circ - 41.0^\circ$).
  2. Quy luật suy giảm nhiệt phi đối xứng: Khi nhiệt độ tăng từ $30^\circ\text{C}$ lên $60^\circ\text{C}$, lực dính đơn vị $c$ suy giảm nghiêm trọng từ $60%$ đến hơn $70%$ (đặc biệt đối với nhựa thông thường 60/70), trong khi góc nội ma sát $\phi$ chỉ suy giảm nhẹ ($1.5^\circ - 3.0^\circ$). Điều này giải thích hiện tượng mặt đường thường bị phá hoại cắt trượt và lún vệt bánh tập trung vào các tháng mùa hè khi nhiệt độ mặt đường đạt ngưỡng $60^\circ\text{C}$.
  3. Hiệu ứng tốc độ gia tải (Strain Rate Sensitivity): Khi tăng tốc độ nén từ $0.5\text{ mm/phút}$ lên $7.5\text{ mm/phút}$, cả $c$ và $\phi$ đều có xu hướng tăng do tính chất nhạy tốc độ biến dạng của chất liên kết bitum nhớt - dẻo. Phân tích phương sai ANOVA khẳng định tác động của tốc độ gia tải đạt mức ý nghĩa thống kê cao ($p < 0.001$).
  4. Phát hiện ngưỡng tỷ số ứng suất cắt chuyển pha biến dạng: Tỷ số $\tau/\tau_f$ là biến số quyết định trạng thái tích lũy biến dạng dẻo:
    • Ở mức $\tau/\tau_f \le 0.269$, đường cong $\varepsilon_p - N$ thuần túy nằm ở pha 2 ổn định, tốc độ biến dạng dẻo tiệm cận hằng số nhỏ.
    • Khi $\tau/\tau_f \ge 0.425$, biến dạng không hồi phục tăng vọt, mẫu nhanh chóng chuyển sang pha 3 (pha phá hoại cắt dẻo với điểm uốn chu kỳ chảy $FN$).
  5. Xác lập mô hình tương quan thực nghiệm $\varepsilon_p - (\tau/\tau_f) - N$: Luận án đã xây dựng thành công bộ phương trình hồi quy phi tuyến mô tả biến dạng không hồi phục tích lũy ở nhiệt độ $60^\circ\text{C}$:
    • Đối với cấp phối BTAC13 sử dụng nhựa PMB III (AC13_CP1_N3): $$\varepsilon_p = 0.0124 \cdot N^{0.218} \cdot \left(\frac{\tau}{\tau_f}\right)^{1.452} \quad (R^2 = 0.942)$$
    • Đối với cấp phối BTAC19 sử dụng nhựa 60/70 (AC19_CP1_N1): $$\varepsilon_p = 0.0187 \cdot N^{0.245} \cdot \left(\frac{\tau}{\tau_f}\right)^{1.618} \quad (R^2 = 0.915)$$
                                  SO SÁNH CÁC NGHIÊN CỨU QUỐC TẾ
                                  
+----------------------+-----------------------------+-----------------------------+-----------------------------+
| Tiêu chí so sánh     | Baek et al. (2008, 2014)    | Zheng et al. (2012, 2016)   | Luận án Ngô Ngọc Quý (2021) |
|                      | (Đại học Sejong - Hàn Quốc) | (ĐH GTVT Tây Nam - TQ)      | (ĐH GTVT - Việt Nam)        |
+----------------------+-----------------------------+-----------------------------+-----------------------------+
| Chủng loại hỗn hợp   | Superpave 19mm, SMA13       | AC13, SMA13, Sup13, PA13    | BTAC12.5, BTAC19            |
| Loại chất kết dính   | PG 64-22, PG 76-22 (SBS)    | Bitum AH-70, SBS biến tính  | Nhựa 60/70, 40/50, PMB III  |
| Nhiệt độ thí nghiệm  | 40°C, 50°C, 60°C            | 25°C, 40°C, 60°C            | 30°C, 45°C, 60°C            |
| Mô hình nén 3 trục   | Triaxial tĩnh & RLPD        | Triaxial tĩnh               | Cải tiến Triaxial & RLPD    |
| Đóng góp mô hình hóa | εp = f(τ/τf, N, t, T)        | Tương quan c, ϕ với lún     | εp = f(τ/τf, N) ở 60°C      |
|                      | R² = 0.87 (kiểm chuẩn APT)  | Biến dạng cắt chiếm 80% lún | ANOVA & Tukey CI 95%        |
+----------------------+-----------------------------+-----------------------------+-----------------------------+

Implications đa chiều

  • Về mặt lý thuyết: Bổ sung cơ sở khoa học cho cơ học vật liệu mặt đường tại Việt Nam, kết nối liền mạch giữa lý thuyết bền chống cắt địa kỹ thuật và động học tích lũy biến dạng dẻo áo đường mềm.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình chế bị mẫu, lắp dựng và vận hành hệ thống nén ba trục tĩnh và động, cung cấp quy trình phân tích thống kê ANOVA - Tukey HSD mẫu mực cho nghiên cứu vật liệu giao thông.
  • Về mặt ứng dụng thực tiễn: Khuyến cáo bắt buộc sử dụng nhựa đường polymer (như PMB III) cho các lớp bê tông asphalt mặt trên (BTAC12.5) tại các đoạn tuyến có mật độ xe nặng cao, xe chạy chậm hoặc dừng đỗ kéo dài (dốc cao, trạm thu phí, nút giao ngã tư), nhằm nâng cao lực dính $c$ và khống chế tỷ số $\tau/\tau_f < 0.25$.
  • Về mặt chính sách và tiêu chuẩn: Cung cấp luận cứ kỹ thuật để Bộ Giao thông Vận tải bổ sung tiêu chí kiểm soát trạng thái giới hạn cắt trượt ($\tau \le [\tau_f]$) vào Tiêu chuẩn thiết kế mặt đường mềm thay thế tiêu chuẩn 22 TCN 211-06 hiện hành.

Limitations và Future Research

  1. Giới hạn phạm vi nhiệt độ trong mô hình RLPD: Mô hình tương quan $\varepsilon_p - (\tau/\tau_f) - N$ mới được xây dựng chi tiết ở nhiệt độ bất lợi $60^\circ\text{C}$; chưa bao quát được phương trình liên tục cho toàn bộ dải nhiệt độ dao động trong năm ($10^\circ\text{C} - 70^\circ\text{C}$).
  2. Giới hạn về chủng loại cấp phối và vật liệu: Nghiên cứu tập trung vào 02 loại cấp phối chặt BTAC12.5 và BTAC19; chưa mở rộng khảo sát cho các hỗn hợp đá - vữa nhựa (SMA), bê tông nhựa rỗng thoát nước (OGFC) hoặc hỗn hợp có sử dụng tỷ lệ cao bê tông nhựa tái chế (RAP).
  3. Mô phỏng cơ chế gia tải thực tế: Thí nghiệm nén ba trục trong phòng thí nghiệm sử dụng sóng tải trọng dọc trục xung đơn phương có áp lực hông không đổi, chưa mô phỏng trọn vẹn sự đổi chiều của vector ứng suất cắt (shear stress reversal) và sự xoay trục ứng suất chính (principal stress rotation) khi bánh xe lăn thực tế trên mặt đường.

Định hướng nghiên cứu tiếp theo (Future Research Agenda):

  • Tích hợp thêm tham số thời gian gia tải ($t$) và dải nhiệt độ ($T$) rộng để phát triển mô hình tổng quát bậc cao dạng Master-Curve cho biến dạng không hồi phục.
  • Ứng dụng phần mềm phần tử hữu hạn phi tuyến 3D (như ABAQUS hoặc ANSYS) kết hợp mô hình vật liệu Mohr-Coulomb/Drucker-Prager cải tiến để mô phỏng chính xác trường ứng suất $\tau/\tau_f$ trong kết cấu áo đường nhiều lớp chịu tải trọng động.
  • Triển khai thí nghiệm kiểm chứng trên mô hình thử nghiệm tăng tốc quy mô lớn (Accelerated Pavement Testing - APT) hoặc theo dõi dài hạn trên các phân đoạn đường cao tốc thực tế tại Việt Nam.

Tác động và ảnh hưởng

                                    MA TRẬN TÁC ĐỘNG ĐA LĨNH VỰC
                                    
  +----------------------+------------------------------------------------------------------+
  | HỌC THUẬT &          | - Tiên phong xác lập hệ dữ liệu sức kháng cắt c, ϕ ba trục BTA.  |
  | NGHIÊN CỨU           | - Tiềm năng trích dẫn cao trong các nghiên cứu M-E Design.        |
  +----------------------+------------------------------------------------------------------+
  | CHUYỂN GIAO          | - Làm chủ công nghệ và quy trình thí nghiệm 3 trục BTA.          |
  | CÔNG NGHỆ            | - Tối ưu hóa thiết kế cấp phối chống hằn lún cho nhà thầu R&D.   |
  +----------------------+------------------------------------------------------------------+
  | QUẢN LÝ &            | - Đề xuất sửa đổi tiêu chuẩn thiết kế áo đường mềm 22 TCN 211-06. |
  | CHÍNH SÁCH           | - Quy định bắt buộc kiểm tra sức kháng cắt tầng mặt đường.       |
  +----------------------+------------------------------------------------------------------+
  | KINH TẾ &            | - Giảm thiểu chi phí sửa chữa hằn lún định kỳ hàng nghìn tỷ VNĐ. |
  | XÃ HỘI               | - Nâng cao an toàn giao thông, hạn chế hiện tượng xe trượt vệt.   |
  +----------------------+------------------------------------------------------------------+
  • Tác động học thuật: Đặt nền móng khoa học vững chắc cho các đề tài sau đại học và các nhóm nghiên cứu chuyên sâu về cơ học mặt đường tại các trường đại học khối kỹ thuật giao thông.
  • Tác động ngành và công nghiệp: Cung cấp công cụ đánh giá chính xác chất lượng nhựa đường biến tính polymer và phụ gia kháng hằn lún cho các doanh nghiệp xây dựng giao thông, các phòng thí nghiệm trọng điểm (LAS-XD).
  • Lợi ích kinh tế - xã hội: Việc kiểm soát và ngăn ngừa hiện tượng lún vệt bánh xe ngay từ khâu thiết kế cấp phối và kết cấu giúp kéo dài tuổi thọ khai thác của mạng lưới đường ôtô cao tốc, tiết kiệm hàng nghìn tỷ đồng kinh phí bảo trì, cào bóc thảm lại hàng năm, đồng thời nâng cao an toàn giao thông phương tiện.

Đối tượng hưởng lợi

  • Nghiên cứu sinh và học viên cao học: Tiếp cận phương pháp luận nghiên cứu thực nghiệm nén ba trục chuẩn tắc, khai thác bộ dữ liệu tham số $c, \phi$ và mô hình $\tau/\tau_f$ để phát triển các bài toán mô phỏng số.
  • Các nhà khoa học và giảng viên chuyên ngành: Sử dụng kết quả luận án làm tài liệu tham khảo chuyên sâu phục vụ công tác giảng dạy môn học Thiết kế mặt đường nâng cao, Cơ học mặt đường và Vật liệu xây dựng đường.
  • Kỹ sư thiết kế và tư vấn giao thông: Ứng dụng quy trình tính toán kiểm tra điều kiện chống cắt trượt để lựa chọn chiều dày hợp lý và chỉ định chủng loại vật liệu (BTAC19 nhựa thường cho lớp dưới, BTAC12.5 nhựa PMB III cho lớp trên).
  • Các nhà quản lý và hoạch định chính sách: Ban hành các quy chuẩn, định mức kinh tế kỹ thuật và hướng dẫn thiết kế hỗn hợp bê tông asphalt theo định hướng hiệu năng thực tế (performance-based specifications).

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Luận án đã mở rộng Thuyết phá hoại cắt Mohr-Coulomb vào lĩnh vực cơ học vật liệu asphalt nhớt - dẻo tại Việt Nam, thiết lập mối liên hệ định lượng giữa độ bền cắt tĩnh ($c, \phi$) với động học tích lũy biến dạng dẻo ($\varepsilon_p$) thông qua tỷ số ứng suất cắt $\tau/\tau_f$.

2. Điểm cải tiến mang tính đột phá về phương pháp luận so với các nghiên cứu trước đây?

So với các nghiên cứu của Nguyễn Quang Phúc (2010) và nhóm tác giả ĐH Xây dựng (2014) chỉ dùng thí nghiệm cắt trực tiếp hoặc nén chẻ, luận án đã thiết kế, cải tiến thành công buồng nén ba trục có kiểm soát nhiệt độ chính xác và thiết bị nén ba trục tải trọng lặp (RLPD), thu thập dữ liệu biến dạng động với tần số lấy mẫu cao qua hệ thống ADC và đầu đo LVDT chuyên dụng.

3. Phát hiện thực nghiệm nào gây bất ngờ nhất và bằng chứng số liệu đi kèm?

Góc nội ma sát $\phi$ của hỗn hợp BTA hầu như không bị suy giảm theo nhiệt độ ở dải $30^\circ\text{C} - 45^\circ\text{C}$ ($38.7^\circ - 44.2^\circ$), chứng minh rằng khả năng chịu lực của khung cốt liệu thô vẫn được bảo toàn nguyên vẹn; hiện tượng mất ổn định trượt cắt mặt đường ở nhiệt độ cao $60^\circ\text{C}$ hoàn toàn bắt nguồn từ sự suy giảm đột ngột (giảm $60-70%$) của lực dính đơn vị $c$ do vữa nhựa bị hóa mềm.

                      BIẾN THIÊN CÁC THAM SỐ KHÁNG CẮT THEO NHIỆT ĐỘ
                      
     Lực dính đơn vị c (MPa)                  Góc nội ma sát ϕ (độ)
      30°C    45°C    60°C (Nhiệt độ)         30°C    45°C    60°C (Nhiệt độ)

4. Luận án có cung cấp quy trình tái lập thực nghiệm (Replication Protocol) hoàn chỉnh không?

Luận án cung cấp đầy đủ và chi tiết toàn bộ quy trình công nghệ: từ thông số thiết kế buồng nén 3 trục, quy cách chế tạo mẫu ép Marshall và đầm xoay, quy trình kiểm soát gia nhiệt màng bao mẫu cao su, chế độ hiệu chuẩn Load cell/LVDT, đến các đoạn mã xử lý phương sai ANOVA và các bảng số liệu gốc tại phần Phụ lục.

5. Chương trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?

Xây dựng một hệ thống phần mềm tính toán kết cấu áo đường mềm hoàn chỉnh theo phương pháp M-E cho điều kiện khí hậu và tải trọng Việt Nam (tương đương AASHTOWare Pavement ME Design), trong đó mô-đun kiểm soát lún cắt trượt dựa trên ma trận dữ liệu $c, \phi, \tau/\tau_f$ là hạt nhân tính toán chính.

Kết luận

Công trình luận án tiến sĩ của tác giả Ngô Ngọc Quý đã hoàn thành xuất sắc các mục tiêu đề ra, mang lại 5 đóng góp cốt lõi:

  1. Hệ thống hóa và làm sáng tỏ bản chất cơ học của hiện tượng hư hỏng hằn lún vệt bánh xe, khẳng định biến dạng trượt cắt nội tại dưới tác dụng của tải trọng trục nặng ở nhiệt độ cao là nguyên nhân cốt tử.
  2. Hoàn thiện thành công hệ thống thiết bị và quy trình thí nghiệm nén ba trục tĩnh và nén ba trục tải trọng lặp (RLPD) chuyên dụng cho bê tông asphalt trong điều kiện nghiên cứu tại Việt Nam.
  3. Xác lập bộ cơ sở dữ liệu thực nghiệm chuẩn xác về lực dính đơn vị ($c$) và góc nội ma sát ($\phi$) của các hỗn hợp BTAC12.5 và BTAC19 với nhiều tổ hợp chất kết dính (60/70, 40/50, PMB III) và nguồn cốt liệu khoáng.
  4. Lượng hóa quy luật suy giảm nhiệt, ảnh hưởng của tốc độ nén và cấu trúc khung đá đến các tham số cường độ chống cắt thông qua phân tích phương sai ANOVA và kiểm định Tukey HSD ở độ tin cậy $95%$.
  5. Thiết lập thành công các mô hình tương quan thực nghiệm $\varepsilon_p = f(\tau/\tau_f, N)$ ở nhiệt độ tính toán $60^\circ\text{C}$, đặt nền móng lý thuyết và thực tiễn vững chắc cho việc chuyển dịch sang phương pháp thiết kế kết cấu áo đường mềm có kiểm soát trạng thái giới hạn biến dạng trượt cắt tại Việt Nam.