Tổng quan về luận án

Trong kỹ thuật xây dựng công trình giao thông hiện đại, hiện tượ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 nặng, lưu lượng xe gia tăng và nền nhiệt độ cao là thách thức kỹ thuật mang tính toàn cầu. Luận án tiến sĩ kỹ thuật với tiêu đề "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 tại Trường Đại học Giao thông Vận tải (2021) dưới sự hướng dẫn của PGS.TS Trần Thị Kim Đăng và PGS.TS Lã Văn Chăm thuộc chuyên ngành Xây dựng đường ôtô và đường thành phố (Mã số: 9.05), là công trình tiên phong thiết lập cơ sở khoa học và thực nghiệm kết nối trạng thái ứng suất phá hoại cắt với động học tích lũy biến dạng dẻo tại Việt Nam.

                  +----------------------------------------------+
                  | Tải trọng trục nặng + Nhiệt độ cao (60°C)   |
                  +----------------------+-----------------------+
                                         |
                                         v
                  +----------------------------------------------+
                  |  Ứng suất cắt cục bộ vượt ngưỡng kháng cắt   |
                  |          tau_max >= tau_f = c + sigma*tan(phi)|
                  +----------------------+-----------------------+
                                         |
                                         v
                  +----------------------------------------------+
                  |       Trượt cắt dẻo & Đẩy trồi hai bên       |
                  |     (Chiếm ~80% tổng chiều sâu vệt hằn lún)   |
                  +----------------------+-----------------------+
                                         |
                                         v
                  +----------------------------------------------+
                  | Mô hình dự báo biến dạng không hồi phục      |
                  |     epsilon_p = f(N, tau/tau_f, T, sigma_3)  |
                  +----------------------------------------------+

Bối cảnh khoa học và tính tiên phong

Trước đây, các nghiên cứu trong nước chủ yếu tiếp cận vấn đề chống lún bằng giải pháp cải thiện thành phần vật liệu đơn thuần (sử dụng nhựa cải tiến polime SBS, tối ưu hóa cấp phối đá dăm theo Quyết định 858/QĐ-BGTVT hoặc TCVN 8819:2011) hoặc đánh giá khả năng kháng hằn lún thông qua các thí nghiệm mô phỏng quy ước như thiết bị hằn lún bánh xe Hamburg (Hamburg Wheel-Tracking). Tuy nhiên, các phương pháp này chưa giải quyết được bản chất cơ học của hiện tượng phá hoại cắt trượt và thiếu vắng mô hình dự báo định lượng biến dạng dẻo theo thời gian khai thác. Công trình của tác giả là nghiên cứu đầu tiên tại Việt Nam cải tiến thành công hệ thống thiết bị thí nghiệm nén ba trục (triaxial shear test) và hệ thống nén dọc trục tải trọng lặp có áp lực hông (Repeated Load Permanent Deformation - RLPD) để giải mã mối quan hệ giữa sức kháng cắt và biến dạng dẻo tích lũy.

Khoảng trống nghiên cứu (Research Gap)

  1. Thiếu vắng thông số cơ học cơ bản: Hệ thống tiêu chuẩn thiết kế áo đường mềm hiện hành ở Việt Nam (22 TCN 211-06) chưa đưa trạng thái giới hạn về ứng suất cắt trượt vào quy trình tính toán kiểm tra kết cấu mặt đường BTA.
  2. Chưa xác lập ngưỡng giá trị $c$ và $\phi$ chuẩn hóa: Chưa có cơ sở dữ liệu thực nghiệm về lực dính đơn vị ($c$) và góc nội ma sát ($\phi$) của các cấp phối BTA thông dụng (BTAC19, BTAC12.5) sử dụng các nguồn cốt liệu nội địa (đá mỏ Đồng Ao, đá mỏ Tân Đông Hiệp) và các loại chất kết dính khác nhau (nhựa 60/70, 40/50, nhựa polime PMB III).
  3. Thiếu mô hình toán học tương quan: Chưa từng có phương trình thực nghiệm nào mô tả định lượng mối quan hệ giữa tỷ số ứng suất cắt trên cường độ kháng cắt ($\tau/\tau_f$) với tốc độ tích lũy biến dạng không hồi phục ($\varepsilon_p$) dưới điều kiện nhiệt độ cao bất lợi ($60^\circ\text{C}$).

Câu hỏi nghiên cứu và Giả thuyết khoa học

  • Câu hỏi nghiên cứu 1 (Q1): Bản chất cơ chế hình thành lún vệt bánh xe liên quan như thế nào đến sự suy giảm cường độ kháng cắt của BTA ở nhiệt độ cao và tốc độ gia tải chậm?
  • Câu hỏi nghiên cứu 2 (Q2): Các thông số kháng cắt $(c, \phi)$ biến thiên ra sao dưới ảnh hưởng tương tác đa yếu tố (nguồn gốc khoáng sét/đá dăm, loại bitum, nhiệt độ, tốc độ biến dạng)?
  • Câu hỏi nghiên cứu 3 (Q3): Liệu có thể mô hình hóa biến dạng không hồi phục tích lũy $\varepsilon_p$ thông qua thông số không thứ nguyên $\tau/\tau_f$ kết hợp với số chu kỳ tác dụng tải trọng $N$?
  • Giả thuyết 1 (H1): Biến dạng dẻo vệt hằn bánh xe không phải do quá trình nén chặt thể tích sau giai đoạn cố kết ban đầu, mà chủ yếu do biến dạng trượt cắt dẻo (plastic shear flow) gây đẩy trồi sang hai bên.
  • Giả thuyết 2 (H2): Tồn tại hàm phụ thuộc phi tuyến đồng nhất giữa biến dạng không hồi phục tích lũy $\varepsilon_p$ và tỷ số ứng suất cắt $\tau/\tau_f$, cho phép dự báo biến dạng dẻo của các loại BTA khác nhau trên cùng một quy luật hàm số.

Khung lý thuyết và Quy mô nghiên cứu

Nghiên cứu tích hợp lý thuyết bền phá hoại Mohr-Coulomb, động học lưu biến vật liệu nhớt-đàn-dẻo (visco-elasto-plasticity), và nguyên lý thiết kế cơ học - thực nghiệm (Mechanistic-Empirical Pavement Design - M-E). Quy mô thực nghiệm bao gồm hàng trăm mẫu nén Marshall, hàng chục tổ mẫu nén ba trục tĩnh và nén ba trục tải trọng lặp trên các cấp phối BTAC19 mỏ Đồng Ao, BTAC12.5 (mỏ Đồng Ao và Tân Đông Hiệp) kết hợp 3 loại bitum (60/70, 40/50, PMB III) tại các dải nhiệt độ $30^\circ\text{C}, 45^\circ\text{C}, 60^\circ\text{C}$ và tốc độ nén từ $0.05\text{ mm/phút}$ đến $50\text{ mm/phút}$.

Literature Review và Positioning

Tổng quan các trường phái nghiên cứu

Trong cơ học mặt đường asphalt, lịch sử nghiên cứu lún vệt bánh xe chia thành hai trường phái học thuật rõ rệt:

  1. Trường phái Cải thiện Thành phần Thực nghiệm (Empirical Material Optimization):

    • Đại diện bởi các nghiên cứu thuộc Chương trình Nghiên cứu Chiến lược Đường bộ Hoa Kỳ (Strategic Highway Research Program - SHRP) với phương pháp thiết kế Superpave và phân loại chất kết dính theo cấp tính năng PG (Performance Grade). Các tác giả như Roberts et al. (1996), Kandhal và Mallick (2001) tập trung vào độ góc cạnh cốt liệu thô/mịn, hàm lượng bột khoáng lọt sàng 0.075mm, và độ nhớt động học của bitum.
    • Tại Việt Nam, nhóm tác giả Trường Đại học Xây dựng (2015, 2017) và Học viện Kỹ thuật Quân sự (2018) đã thử nghiệm đặc trưng kháng cắt bằng mô hình cắt phẳng Leutner hoặc phối hợp nén dọc trục và ép chẻ theo đường sinh theo tiêu chuẩn Nga ГОСТ 12801-98. Tuy nhiên, các mô hình này không thiết lập được trạng thái ứng suất khối 3 chiều.

  2. Trường phái Cơ học - Thực nghiệm và Phân tích Trạng thái Ứng suất Khối (Mechanistic-Empirical & Triaxial Stress Modeling):

    • Xuất phát từ lý thuyết cơ học môi trường liên tục, Van der Loo (1978) đã chỉ ra rằng thí nghiệm từ biến nén một trục không nở hông tự do không thể mô phỏng đúng trường biến dạng thực tế ngoài hiện trường.
    • Smith (1990) tại Đại học Quốc gia Singapore đề xuất đường bao giới hạn phá hoại $c-\phi$ dựa trên phân tích phần tử hữu hạn (FEM) để tích hợp thiết kế hỗn hợp và kết cấu áo đường.
    • Monismith et al. (1994) và Sousa et al. (1994) chứng minh rằng lún vệt bánh xe hình thành chủ yếu do biến dạng trượt cắt khi ứng suất cắt vượt quá giới hạn dẻo của vật liệu.
                           KHUNG LITERATURE POSITIONING
===================================================================================
Trường phái Cải thiện Vật liệu              Trường phái Cơ học - Thực nghiệm (M-E)
(SHRP, Superpave, Marshall, Hveem)          (Shell, Monismith, Sousa, M-EPDG)
                \                                        /
                 \                                      /
                  v                                    v
       Tranh luận học thuật: Lún do nén chặt (Densification) vs. Biến dạng trượt cắt (Plastic Shear)
       -> Thực nghiệm APT & Shell chứng minh: Trượt cắt dẻo chiếm ~80% chiều sâu lún.
                                      |
                                      v
       Nghiên cứu quốc tế: Baek et al. (Hàn Quốc), Zheng et al. (Trung Quốc), Smith (Singapore)
                                      |
                                      v
     [ KHOẢNG TRỐNG TẠI VN: Thiếu thiết bị 3 trục, chưa có hàm liên hệ tau/tau_f & epsilon_p ]
                                      |
                                      v
      CÔNG TRÌNH LUẬN ÁN (NGÔ NGỌC QUÝ, 2021): Cải tiến thiết bị 3 trục + RLPD -> Mô hình epsilon_p
===================================================================================

Tranh luận học thuật cốt lõi (Scientific Debates)

Một cuộc tranh luận kéo dài trong công nghệ áo đường xoay quanh cơ chế biến dạng chính gây hằn lún: Lún do nén chặt thể tích (Densification) hay Lún do biến dạng trượt cắt dẻo (Plastic Shear Deformation)?

  • Quan điểm ban đầu dựa trên mô hình nén chặt cho rằng vệt lún xuất hiện do lỗ rỗng dư trong bê tông giảm dần dưới tác động đầm lèn của bánh xe.
  • Tuy nhiên, các thực nghiệm hiện trường và thử nghiệm gia tốc mặt đường (Accelerated Pavement Testing - APT) của AASHTO, Shell và nghiên cứu của Zheng et al. (2007) tại Đại học Giao thông Tây Nam (Trung Quốc) đã bác bỏ giả thuyết trên. Kết quả khẳng định: giai đoạn nén chặt thể tích chỉ diễn ra trong thời gian ngắn ban đầu, trong khi biến dạng trượt cắt kèm theo hiện tượng trồi hai bên vệt bánh xe chiếm tới 80% tổng chiều sâu lún. Khi trích dẫn nguyên văn từ các kết quả cơ học biến dạng: "lún vệt bánh xe phần lớn hình thành do biến dạng cắt của hỗn hợp bê tông nhựa... xảy ra khi ứng suất cắt vượt quá cường độ kháng cắt của bê tông asphalt hoặc ứng suất cắt đủ để biến dạng từ biến tích lũy lại".

So sánh với các công trình quốc tế tiêu biểu

  1. So sánh với nghiên cứu của Zheng, Xiao et al. (2007 - Trung Quốc): Zheng et al. sử dụng thí nghiệm nén ba trục trên các hỗn hợp AC13, SMA13, Sup13 để khảo sát quan hệ giữa $(c, \phi)$ với độ sâu lún vệt bánh xe Wheel-Tracking. Công trình của Ngô Ngọc Quý phát triển xa hơn khi không chỉ dừng lại ở chỉ tiêu hằn lún quy ước mà đi sâu vào động học biến dạng không hồi phục qua thí nghiệm tải trọng lặp RLPD với áp lực hông.
  2. So sánh với nghiên cứu của Baek, Park et al. (2008, 2009 - Đại học Sejong, Hàn Quốc): Baek et al. xây dựng tương quan giữa tỷ số $\tau/\tau_f$ và biến dạng không hồi phục $\varepsilon_p$ cho vật liệu BTA tại Hàn Quốc ở $50^\circ\text{C}$. Luận án của Ngô Ngọc Quý đã mở rộng và tái định chuẩn mô hình này trong điều kiện khắc nghiệt của khí hậu nhiệt đới gió mùa Việt Nam ở mức nhiệt độ bất lợi $60^\circ\text{C}$, đồng thời kiểm chứng trên cả cấp phối BTAC19 thông thường và BTAC12.5 sử dụng bitum cải tiến polime PMB III.

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

                     KHUNG TÍCH HỢP LÝ THUYẾT VÀ PHÂN TÍCH
+-----------------------------------------------------------------------------------+
| 1. LÝ THUYẾT BỀN MOHR-COULOMB:  tau_f = c + sigma * tan(phi)                      |
|    Xác định khả năng chịu cắt cực hạn từ lực dính đơn vị (c) và góc nội ma sát    |
+------------------------------------------+----------------------------------------+
                                           |
                                           v
+-----------------------------------------------------------------------------------+
| 2. LÝ THUYẾT NHỚT-ĐÀN-DẺO (VISCOELASTICITY - CREEP & FLOW THEORY):                 |
|    Xác định các pha biến dạng: Cố kết ban đầu -> Từ biến ổn định -> Phá hoại dẻo  |
|    Thông số kiểm soát: Flow Number (FN), Flow Time (FT), Modulus Estiff           |
+------------------------------------------+----------------------------------------+
                                           |
                                           v
+-----------------------------------------------------------------------------------+
| 3. NGUYÊN LÝ THIẾT KẾ CƠ HỌC - THỰC NGHIỆM (M-E DESIGN):                          |
|    Chuẩn hóa trạng thái ứng suất qua Tỷ số ứng suất cắt: SSR = tau / tau_f        |
+------------------------------------------+----------------------------------------+
                                           |
                                           v
+-----------------------------------------------------------------------------------+
| MÔ HÌNH TOÁN HỌC DỰ BÁO:  epsilon_p = A * (N)^B * (tau / tau_f)^C                 |
| R^2 > 0.95 với điều kiện biên: T = 60°C; sigma_3 = 69-138 kPa; f = 1Hz (0.1s/0.9s)|
+-----------------------------------------------------------------------------------+

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

  1. Mở rộng phạm vi ứng dụng của Định luật bền Mohr-Coulomb cho vật liệu gốc bitum: Phương trình kinh điển: $$\tau_f = c + \sigma \tan \phi$$ được luận án chứng minh là hoàn toàn áp dụng được để mô tả trạng thái giới hạn cắt của BTA chặt. Trong đó, tác giả đã tách biệt rõ rệt nguồn gốc của hai thành phần: lực dính đơn vị $c$ đặc trưng cho độ quánh và tính dính kết của màng vữa bitum - bột khoáng; góc nội ma sát $\phi$ đặc trưng cho sự chèn móc, đan cài của khung xương cốt liệu thô.
  2. Xác lập tham số trạng thái không thứ nguyên $\tau/\tau_f$ làm cầu nối chuyển tiếp pha biến dạng: Luận án làm sáng tỏ rằng mức độ tích lũy biến dạng dẻo không đơn thuần phụ thuộc vào độ lớn tuyệt đối của tải trọng trục $\sigma_1$, mà bị chi phối trực tiếp bởi tỷ số ứng suất cắt tác dụng trên cường độ kháng cắt cực hạn $\tau/\tau_f$ (Shear Stress Ratio - SSR). Khi tỷ số này vượt qua một ngưỡng tới hạn, vật liệu sẽ chuyển nhanh từ pha biến dạng ổn định (pha 2) sang pha phá hoại trượt dẻo dòng chảy (pha 3 - Tertiary Flow).
  3. Phát triển bước tiến mô hình hóa cho Thiết kế Cơ học - Thực nghiệm (M-E): Thay thế các công thức kinh nghiệm thuần túy bằng phương trình tích lũy biến dạng không hồi phục có xét đến thông số kháng cắt thực tế của vật liệu, tạo tiền đề toán học cho việc kiểm soát lún ngay từ khâu thiết kế thành phần hỗn hợp.

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

Khung phân tích của luận án tích hợp đồng thời 3 nền tảng lý thuyết:

  • Lý thuyết phá hoại Mohr-Coulomb (xác định ngưỡng kháng cắt $\tau_f$).
  • Lý thuyết lưu biến và từ biến dẻo (xác định tốc độ biến dạng $\dot{\varepsilon}$ và chu kỳ chảy $FN$).
  • Phương pháp tương quan thống kê phi tuyến đa biến (xây dựng hàm hồi quy liên kết giữa $\varepsilon_p$, $N$, và $\tau/\tau_f$).

Điều kiện biên xác định (Boundary Conditions):

  • Chế độ nhiệt độ chuẩn: $T = 60^\circ\text{C}$ (nhiệt độ thiết kế mặt đường mùa nóng tại Việt Nam).
  • Áp lực hông hạn chế nở ngang: $\sigma_3 = 69 \text{ kPa} \div 138 \text{ kPa}$.
  • Tần số tác dụng tải lặp: $f = 1 \text{ Hz}$ (thời gian tác dụng tải xung $t_{load} = 0.1\text{s}$, thời gian nghỉ $t_{rest} = 0.9\text{s}$).

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

Thiết kế nghiên cứu

  • Triết lý nghiên cứu: Thực chứng luận (Positivism) kết hợp phương pháp luận Cơ học - Thực nghiệm (Mechanistic-Empirical).
  • Thiết kế giai tầng đa mức (Multi-level Factorial Design): Đánh giá tương tác đồng thời giữa 4 nhóm biến số độc lập:
    1. Biến số vật liệu: 2 nguồn cốt liệu (Mỏ Đồng Ao - đá vôi/dolomit cứng; Mỏ Tân Đông Hiệp - đá andesit/granit), 3 loại mác nhựa (60/70, 40/50, PMB III).
    2. Biến số cấp phối: Cấp phối BTAC19 (Dmax = 19 mm) và BTAC12.5 / AC13 (Dmax = 12.5 mm).
    3. Biến số nhiệt độ: $30^\circ\text{C}, 45^\circ\text{C}, 60^\circ\text{C}$.
    4. Biến số cơ học: Tốc độ nén biến dạng ($0.05, 0.5, 5.0, 7.5, 50\text{ mm/phút}$), áp lực hông $\sigma_3$ ($0, 50, 100, 200\text{ kPa}$), mức tải trọng lặp $\tau/\tau_f$ ($0.175, 0.269, 0.425$).
                    MA TRẬN THỰC NGHIỆM ĐA YẾU TỐ
===============================================================================
[Vật liệu đầu vào]       [Thí nghiệm Kháng cắt]          [Thí nghiệm Biến dạng dẻo]
- Cốt liệu: Đồng Ao,      - Nén 3 trục tĩnh (cải tiến)    - Creep tĩnh (SCT):
  Tân Đông Hiệp             (sigma_3 = 0, 50, 100 kPa)      (T = 30, 45, 60°C)
- Nhựa: 60/70, 40/50,     - Vận tốc: 0.05 - 50 mm/phút    - Tải lặp 3 trục (RLPD):
  PMB III                 - T = 30°C, 45°C, 60°C            (tau/tau_f = 0.175, 0.269, 0.425)
- Cấp phối: AC19, AC13    -> Xác định c, phi                -> Xác định epsilon_p, FN, Estiff
===============================================================================

Quy trình nghiên cứu thực nghiệm và Cải tiến thiết bị

  1. Thiết kế hỗn hợp và chế tạo mẫu:
    • Thiết kế thành phần cấp phối cốt liệu và hàm lượng nhựa tối ưu theo phương pháp Marshall kết hợp kiểm tra độ rỗng dư ($VA = 4.0% \pm 0.2%$).
    • Chế bị mẫu hình trụ tiêu chuẩn ($D = 100\text{ mm}, H = 150\text{ mm}$ đến $200\text{ mm}$) đầm nén bằng thiết bị đầm xoay Superpave (SGC) và khoan cắt gọt phẳng hai mặt đáy với độ sai lệch song song $< 0.05\text{ mm}$.
  2. Cải tiến thiết bị thí nghiệm nén ba trục:
    • Luận án đã có đóng góp kỹ thuật quan trọng khi thiết kế cải tiến buồng nén ba trục của đất trong phòng thí nghiệm Vilas 047 và LasXD 1256 sang thiết bị nén ba trục cho BTA: bổ sung áo gia nhiệt kiểm soát nhiệt độ chính xác $\pm 0.5^\circ\text{C}$, lắp đặt màng cao su chịu nhiệt, hệ thống tạo áp lực hông khí nén ổn định.
    • Luận chứng và đánh giá thực nghiệm ảnh hưởng của màng bao mẫu (membrane effect) đến ứng suất chính $\sigma_1$ và biến dạng dọc trục $\varepsilon_1$, chứng minh sai số do màng bao nằm trong khoảng tin cậy 95% cho phép.
  3. Lắp dựng hệ thống nén dọc trục tải trọng lặp (RLPD):
    • Tích hợp Loadcell cảm biến lực độ chính xác cao, 02 cảm biến biến dạng vi sai tuyến tính LVDT (model CDP-25), bộ chuyển đổi tín hiệu tương tự - số ADC tốc độ cao kết nối máy tính ghi nhận dữ liệu liên tục theo thời gian thực.
                      SƠ ĐỒ BUỒNG NÉN BA TRỤC CẢI TIẾN
                    +--------------------------------+
                    |        Piston gia tải          |
                    |              |                 |
                    |              v                 |
                    |    +----------------------+    |
                    |    |    Thớt nén trên     |    |
                    |    +----------------------+    |
                    |    |                      |    |
   Khí nén áp lực --+--->|   Mẫu BTA (D100xH150)|<---+--- Áp lực hông sigma_3
   hông sigma_3     |    |   Màng cao su bao    |    |
                    |    |   LVDT đo biến dạng  |    |
                    |    |                      |    |
                    |    +----------------------+    |
                    |    |    Thớt nén dưới     |    |
                    |    +----------------------+    |
                    |    | Áo gia nhiệt (±0.5°C)|    |
                    +----+----------------------+----+

Xử lý dữ liệu và Phân tích thống kê

  • Phân tích phương sai (ANOVA): Kiểm tra mức độ ảnh hưởng của các yếu tố nguồn gốc đá, loại nhựa, nhiệt độ và tốc độ gia tải lên các biến phụ thuộc $(c, \phi, \varepsilon_p)$ với mức ý nghĩa $\alpha = 0.05$.
  • Kiểm định so sánh bội Tukey HSD: Phân tách nhóm có ý nghĩa thống kê giữa các cặp xử lý thực nghiệm.
  • Phân tích hồi quy phi tuyến: Xây dựng mô hình toán học tương quan sử dụng thuật toán bình phương tối thiểu Gauss-Newton, kiểm soát sai số tương đối và đánh giá độ tương thích qua hệ số xác định $R^2$, sai số chuẩn hồi quy ($S_e/S_y$).

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

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

               SO SÁNH CÁC THÔNG SỐ KHÁNG CẮT (TẠI 60°C, 7.5 mm/phút)
+-------------------------+--------------------+---------------------+
| Hỗn hợp BTA             | Lực dính c (MPa)   | Góc nội ma sát phi  |
+-------------------------+--------------------+---------------------+
| BTAC19 - Nhựa 60/70     |    0.20 - 0.23     |     41.2° - 42.5°   |
| BTAC13 - Nhựa 60/70     |    0.22 - 0.25     |     38.5° - 39.8°   |
| BTAC13 - Nhựa 40/50     |    0.28 - 0.31     |     39.2° - 40.1°   |
| BTAC13 - Nhựa PMB III   |    0.38 - 0.44     |     42.8° - 44.5°   |
+-------------------------+--------------------+---------------------+

  1. Hiệu ứng vượt trội của bitum polyme cải tiến PMB III lên lực dính đơn vị: Phân tích ANOVA và kiểm định Tukey khẳng định loại chất kết dính là yếu tố quyết định hàng đầu đến lực dính đơn vị $c$. Ở $60^\circ\text{C}$, việc chuyển đổi từ nhựa thông thường 60/70 sang nhựa cải tiến polime PMB III làm lực dính đơn vị $c$ tăng từ $0.22\text{ MPa}$ lên $0.41\text{ MPa}$ (tăng gần 86%), đồng thời giúp góc nội ma sát $\phi$ tăng thêm từ $2^\circ$ đến $4^\circ$ do cấu trúc mạng lưới polyme liên kết giữ chặt các hạt cốt liệu mịn.

  2. Quy luật suy giảm phi tuyến của sức kháng cắt theo nhiệt độ: Khi nhiệt độ tăng từ $30^\circ\text{C}$ lên $60^\circ\text{C}$, lực dính đơn vị $c$ của hỗn hợp BTAC19 nhựa 60/70 suy giảm nghiêm trọng tới $65% \div 70%$ (từ xấp xỉ $0.72\text{ MPa}$ xuống còn $0.21\text{ MPa}$). Trái lại, góc nội ma sát $\phi$ có độ nhạy cảm với nhiệt độ thấp hơn nhiều (chỉ giảm nhẹ từ $44.5^\circ$ xuống $41.8^\circ$), chứng minh rằng khả năng kháng tải ở nhiệt độ cao phụ thuộc chủ yếu vào sự tiếp xúc và chèn móc của bộ khung cốt liệu thô.

  3. Độ nhạy cảm của sức kháng cắt đối với tốc độ gia tải: Khi tốc độ gia tải tăng từ $0.05\text{ mm/phút}$ (mô phỏng xe dừng đỗ hoặc leo dốc tốc độ cực chậm) lên $50\text{ mm/phút}$ (mô phỏng dòng xe chạy tốc độ cao), lực dính $c$ tăng hơn 2.5 lần. Điều này giải thích tại sao hiện tượng hằn lún vệt bánh xe luôn diễn ra trầm trọng nhất ở các đoạn đường dốc, trạm thu phí, nút giao đèn tín hiệu giao thông.

  4. Xác lập Phương trình Tương quan Thực nghiệm giữa $\tau/\tau_f$ và $\varepsilon_p$: Luận án đã xây dựng thành công phương trình hồi quy dự báo biến dạng không hồi phục tích lũy cho hai loại hỗn hợp chủ lực tại Việt Nam ở $60^\circ\text{C}$:

    • Đối với cấp phối BTAC19 sử dụng nhựa 60/70 (lớp mặt dưới): $$\varepsilon_p = 0.00142 \cdot N^{0.284} \cdot \left(\frac{\tau}{\tau_f}\right)^{1.452} \quad (R^2 = 0.962)$$
    • Đối với cấp phối BTAC12.5 sử dụng nhựa PMB III (lớp mặt trên): $$\varepsilon_p = 0.00085 \cdot N^{0.241} \cdot \left(\frac{\tau}{\tau_f}\right)^{1.318} \quad (R^2 = 0.978)$$

    Kiểm chứng thống kê cho thấy các mô hình đạt hệ số tương quan $R^2 > 0.95$, sai số dự báo so với số liệu thực nghiệm RLPD nhỏ hơn $8.5%$.

                      ĐỒ THỊ TÍCH LŨY BIẾN DẠNG DẺO epsilon_p
    epsilon_p (%)
         ^
         |                                          / Mức tải cao (tau/tau_f = 0.425)
     4.0 |                                         /  [Giai đoạn 3: Phá hoại dẻo]
         |                                        /
     3.0 |                                      /--- Mức tải TB (tau/tau_f = 0.269)
         |                        -------------/     [Giai đoạn 2: Tốc độ ổn định]
     2.0 |            -----------/
         |       ----/------------------------------ Mức tải thấp (tau/tau_f = 0.175)
     1.0 |  ----/
         | / [Giai đoạn 1: Cố kết ban đầu]
       0 +-------------------------------------------------->
         0        10,000      20,000      30,000     40,000   Chu kỳ tải (N)

Implications đa chiều

  • Ý nghĩa khoa học: Chuyển biến phương pháp luận đánh giá chất lượng BTA từ các chỉ tiêu tĩnh quy ước (độ ổn định Marshall, độ dẻo) sang các tham số cơ học thực chất $(c, \phi, \tau/\tau_f, FN)$.
  • Phương pháp luận ứng dụng: Cho phép các phòng thí nghiệm chuyên ngành tại Việt Nam sử dụng các thiết bị nén ba trục cải tiến để đánh giá chất lượng hỗn hợp với chi phí thấp hơn nhiều so với việc nhập khẩu các hệ thống thí nghiệm phức tạp của nước ngoài.
  • Thực tiễn thiết kế và thi công: Cung cấp công cụ cho các kỹ sư tư vấn thiết kế tính toán kiểm tra chiều sâu lún vệt bánh xe dự kiến theo tuổi thọ khai thác $N$ trục xe tiêu chuẩn, từ đó lựa chọn chiều dày và loại vật liệu phù hợp cho từng tầng mặt đường.

Limitations và Future Research

Hạn chế học thuật (Limitations)

  1. Phạm vi nhiệt độ mô hình tương quan: Phương trình toán học liên hệ giữa $\tau/\tau_f$ và $\varepsilon_p$ trong luận án mới được thiết lập chi tiết tại một mức nhiệt độ tham chiếu chuẩn là $60^\circ\text{C}$. Dù đây là nhiệt độ tới hạn tính toán lún, mô hình cần được mở rộng sang dải nhiệt độ biến thiên liên tục ($20^\circ\text{C} \div 70^\circ\text{C}$).
  2. Quy mô loại hình hỗn hợp: Nghiên cứu mới tập trung giải quyết trên 2 loại cấp phối chặt truyền thống (BTAC19 nhựa 60/70 và BTAC12.5 nhựa PMB III), chưa bao phủ hết các hỗn hợp đặc biệt như bê tông nhựa rỗng thoát nước (OGFC), đá - vữa nhựa (SMA), hay bê tông nhựa tái sinh cào bóc nguội (RAP).
  3. Mô phỏng điều kiện áp lực nước lỗ rỗng và độ ẩm: Thí nghiệm nén ba trục được thực hiện trong điều kiện mẫu khô, chưa xét đến tác động phá hủy kết hợp của độ ẩm và áp lực nước thủy động (moisture-induced damage).

Hướng nghiên cứu tiếp theo (Future Research Agenda)

  • Mở rộng mô hình tương quan đa biến tổng quát (Level 3 Model): Tích hợp đồng thời nhiệt độ $T$, độ ẩm, hàm lượng rỗng dư $VA$, và thời gian tác dụng tải $t_{load}$ vào một hàm trạng thái duy nhất.
  • Phát triển chương trình phần tử hữu hạn (FEM) chuyên dụng: Lập trình chương trình tính toán kết cấu áo đường mềm nhiều lớp phi tuyến dựa trên đường cong quan hệ ứng suất - biến dạng dẻo thực nghiệm xác định từ thí nghiệm nén 3 trục.
  • Kiểm chuẩn trên quy mô hiện trường thực tế (Full-Scale Accelerated Pavement Testing): Xây dựng các đoạn đường thử nghiệm có đặt các đầu đo cảm biến biến dạng sâu (Multi-Depth Deflectometers - MDD) để hiệu chỉnh hệ số tương quan giữa dự báo trong phòng và thực tế khai thác.

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

                                HỆ SINH THÁI TÁC ĐỘNG
+-----------------------------------------------------------------------------------+
| 1. GIỚI HỌC THUẬT & ĐÀO TẠO:                                                      |
|    - Cung cấp dữ liệu chuẩn c, phi cho các nghiên cứu sinh, học viên cao học      |
|    - Tài liệu tham khảo nền tảng cho môn học Cơ học mặt đường nâng cao            |
+------------------------------------------+----------------------------------------+
                                           |
                                           v
+-----------------------------------------------------------------------------------+
| 2. CƠ QUAN QUẢN LÝ NHÀ NƯỚC (BỘ GIAO THÔNG VẬN TẢI):                              |
|    - Cơ sở khoa học để cập nhật tiêu chuẩn thiết kế áo đường mềm thay thế 22TCN211|
|    - Hoàn thiện khung kỹ thuật kiểm soát lún theo Quyết định 858/QĐ-BGTVT         |
+------------------------------------------+----------------------------------------+
                                           |
                                           v
+-----------------------------------------------------------------------------------+
| 3. DOANH NGHIỆP TƯ VẤN, THI CÔNG & CHỦ ĐẦU TƯ CAO TỐC:                            |
|    - Giảm thiểu rủi ro hằn lún vệt bánh xe trên các tuyến huyết mạch              |
|    - Tiết kiệm hàng trăm tỷ đồng chi phí cào bóc, sửa chữa mặt đường định kỳ      |
+-----------------------------------------------------------------------------------+
  • Tác động học thuật: Công trình mở ra hướng nghiên cứu chuyên sâu về cơ học cắt trượt trong vật liệu áo đường tại Việt Nam, kích thích hàng loạt các đề tài nghiên cứu tiếp nối về ứng xử mỏi và trượt dẻo của hỗn hợp asphalt cải tiến.
  • Ảnh hưởng chính sách và tiêu chuẩn hóa: Đóng góp luận cứ khoa học thực nghiệm trực tiếp cho Bộ Giao thông Vận tải trong việc rà soát, điều chỉnh Tiêu chuẩn thiết kế mặt đường mềm, định hướng chuyển dịch từ phương pháp thiết kế tĩnh truyền thống sang phương pháp thiết kế Cơ học - Thực nghiệm (M-E).
  • Lợi ích kinh tế - xã hội: Ứng dụng kết quả nghiên cứu giúp tối ưu hóa chiều dày kết cấu mặt đường, lựa chọn đúng loại bitum polyme cho tầng mặt trên chịu ứng suất cắt lớn, từ đó hạn chế triệt để hư hỏng lún vệt bánh xe, nâng cao an toàn giao thông, kéo dài tuổi thọ công trình và tiết kiệm chi phí duy tu sửa chữa hàng năm cho ngân sách quốc gia.

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

  • Nghiên cứu sinh, Học viên Cao học ngành Công trình Giao thông: Tiếp cận phương pháp luận nghiên cứu thực nghiệm chuẩn mực, quy trình cải tiến thiết bị nén ba trục, và thuật toán phân tích phương sai ANOVA/Tukey xử lý số liệu vật liệu xây dựng.
  • Các Viện nghiên cứu và Trường Đại học Kỹ thuật: Sở hữu tài liệu tham khảo chuyên sâu và bộ số liệu thực nghiệm gốc về đặc tính kháng cắt của vật liệu BTA tại Việt Nam.
  • Kỹ sư Tư vấn Thiết kế đường bộ: Nắm bắt phương pháp tính toán dự báo biến dạng dẻo $\varepsilon_p$ theo tỷ số $\tau/\tau_f$ để chủ động kiểm soát trạng thái hằn lún ngay từ giai đoạn lập dự án kỹ thuật.
  • Chủ đầu tư và Cơ quan Quản lý Đường cao tốc: Có căn cứ khoa học chuẩn xác để nghiệm thu, lựa chọn vật liệu asphalt chất lượng cao, giám sát chặt chẽ độ ổn định nhiệt độ và chất lượng thi công mặt đường.

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?

Đóng góp lý thuyết độc đáo nhất là việc chứng minh và lượng hóa thành công sự phụ thuộc của biến dạng không hồi phục tích lũy $\varepsilon_p$ vào Tỷ số ứng suất cắt trên sức kháng cắt $\tau/\tau_f$ (SSR) dưới dạng hàm số mũ phi tuyến đồng nhất. Nghiên cứu đã mở rộng Định luật phá hoại Mohr-Coulomb từ một mô hình trạng thái giới hạn phá hủy tĩnh sang môi trường động học lưu biến nhớt-đàn-dẻo dưới tải trọng lặp, thiết lập mối liên kết toán học trực tiếp giữa khả năng chống cắt trượt cực hạn và tốc độ phát triển biến dạng dẻo theo thời gian.

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

So với các nghiên cứu trong nước trước đây (chỉ dùng mô hình cắt phẳng Leutner hoặc ép chẻ gián tiếp không có áp lực hông) và các nghiên cứu quốc tế dùng thiết bị đắt tiền, tác giả đã:

  1. Thiết kế, chế tạo cải tiến thành công buồng thí nghiệm nén 3 trục và nén dọc trục tải trọng lặp (RLPD) có kiểm soát nhiệt độ chính xác $\pm 0.5^\circ\text{C}$ từ thiết bị cơ học đất sẵn có, giúp chuẩn hóa điều kiện biên ứng suất khối $(\sigma_1, \sigma_3)$ cho mẫu asphalt đường kính lớn ($100\text{ mm}$).
  2. Loại trừ và kiểm chuẩn định lượng thành công ảnh hưởng của màng bao cao su (membrane effect) lên số liệu đo biến dạng và ứng suất chính, đảm bảo tính khách quan khoa học cao nhất.

3. Phát hiện thực nghiệm nào gây bất ngờ nhất trong quá trình nghiên cứu?

Phát hiện bất ngờ và có giá trị nhất là: Góc nội ma sát $\phi$ của hỗn hợp BTA hầu như không bị suy giảm đáng kể theo nhiệt độ (chỉ giảm khoảng $2^\circ \div 3^\circ$ khi nhiệt độ tăng từ $30^\circ\text{C}$ lên $60^\circ\text{C}$), trong khi lực dính đơn vị $c$ lại sụt giảm tới gần $70%$. Điều này chứng minh một cách tường minh rằng ở nhiệt độ khai thác mùa hè khắc nghiệt ($60^\circ\text{C}$), màng nhựa bitum bị hóa mềm và mất hầu hết khả năng liên kết dính bám; toàn bộ khả năng kháng cắt và chống lún của mặt đường khi đó phải dựa vào sự chèn móc cơ học của bộ khung cốt liệu đá dăm thô có góc cạnh cao.

4. Quy trình thí nghiệm của luận án có khả năng lặp lại và chuyển giao (Replication Protocol) không?

Quy trình hoàn toàn có khả năng tái lập và chuyển giao công nghệ 100%. Luận án đã mô tả chi tiết:

  • Kích thước chế tạo mẫu, số vòng đầm xoay SGC, quy trình cắt gọt mẫu phẳng đáy.
  • Sơ đồ mạch đo LVDT, hiệu chuẩn cảm biến Loadcell và CDP-25.
  • Trình tự thiết lập áp lực buồng khí nén $\sigma_3$, tốc độ gia tải chuẩn $7.5\text{ mm/phút}$ cho nén tĩnh và xung tải $1\text{ Hz}$ ($0.1\text{s}$ nén / $0.9\text{s}$ nghỉ) cho nén lặp, hoàn toàn tương thích với các phòng thí nghiệm LasXD tại Việt Nam.

5. Khung chương trình nghiên cứu 10 năm tới (10-Year Research Agenda) được định hình ra sao?

  1. Năm 1-3: Xây dựng phần mềm tích hợp mô hình $\varepsilon_p = f(N, \tau/\tau_f)$ vào bài toán kết cấu áo đường mềm nhiều lớp phục vụ thiết kế thực tiễn.
  2. Năm 4-6: Mở rộng ngân hàng dữ liệu $c, \phi$ cho toàn bộ các loại vật liệu mới (SMA, Asphalt cao su hóa CR-PMA, RAP tái sinh nguội).
  3. Năm 7-10: Thiết lập hệ số chuyển đổi thực nghiệm từ phòng thí nghiệm ra hiện trường thông qua quan trắc trạm thử nghiệm gia tốc mặt đường APT và hệ thống trạm cân tải trọng tự động WIM (Weigh-In-Motion), hoàn thiện trọn vẹn Tiêu chuẩn Thiết kế Áo đường M-E của Việt Nam.

Kết luận

  1. Khẳng định bản chất cơ chế hằn lún: Luận án đã chứng minh một cách khoa học rằng hằn lún vệt bánh xe ở nhiệt độ cao chủ yếu do biến dạng trượt cắt dẻo của hỗn hợp bê tông asphalt gây nên, trong đó trạng thái ứng suất cắt cục bộ vượt quá ngưỡng cường độ kháng cắt cực hạn của vật liệu.
  2. Cải tiến thành công hệ thống thiết bị thí nghiệm chuyên dụng: Chế tạo và chuẩn hóa thành công buồng nén ba trục tĩnh và nén ba trục tải trọng lặp (RLPD) kiểm soát nhiệt độ từ thiết bị cơ học đất, mở ra giải pháp thực nghiệm cơ học khối ba chiều chi phí thấp, độ chính xác cao cho ngành giao thông Việt Nam.
  3. Bộ cơ sở dữ liệu thực nghiệm chuẩn hóa đầu tiên: Xác lập bộ thông số toàn diện về lực dính đơn vị $c$ và góc nội ma sát $\phi$ của các cấp phối BTAC19, BTAC12.5 kết hợp nhựa 60/70, 40/50 và PMB III dưới các dải nhiệt độ ($30^\circ\text{C} \div 60^\circ\text{C}$) và vận tốc biến dạng khác nhau.
  4. Thiết lập mô hình toán học tương quan có giá trị ứng dụng cao: Xây dựng thành công phương trình hồi quy phi tuyến dự báo biến dạng không hồi phục tích lũy $\varepsilon_p$ theo chu kỳ tải $N$ và tỷ số ứng suất cắt $\tau/\tau_f$ với độ tin cậy $R^2 > 0.95$.
  5. Định hướng ứng dụng đột phá trong thiết kế kết cấu: Tạo tiền đề chuyển giao công nghệ tính toán kết cấu áo đường từ phương pháp đàn hồi quy ước sang phương pháp Cơ học - Thực nghiệm (M-E), trực tiếp phục vụ mục tiêu nâng cao chất lượng khai thác và tuổi thọ của hệ thống đường cao tốc tại Việt Nam.