Tổng quan về luận án

Nghiên cứu của tác giả Bùi Văn Phú (2024) tại Trường Đại học Giao thông Vận tải với đề tài "Nghiên cứu ứng xử cơ học của vật liệu và kết cấu áo đường mềm dưới tác dụng của tải trọng động trong điều kiện Việt Nam" (Mã số: 9580205, Chuyên ngành: Kỹ thuật xây dựng công trình giao thông, dưới sự hướng dẫn của PGS. Nguyễn Quang Tuấn và PGS. Nguyễn Quang Phúc) là một công trình khoa học mang tính đột phá và hệ thống. Đề tài giải quyết trực diện điểm nghẽn tồn tại nhiều thập kỷ trong ngành công trình giao thông tại Việt Nam: sự thiếu tương thích giữa các mô hình thiết kế áo đường truyền thống dựa trên giả định đàn hồi tĩnh tuyến tính với ứng xử cơ học thực tế phức tạp của vật liệu dưới tác dụng của tải trọng trục xe động và điều kiện khí hậu nhiệt đới gió mùa.

Về khoảng trống nghiên cứu (research gap), thực tiễn thiết kế kết cấu áo đường mềm tại Việt Nam hiện nay chủ yếu áp dụng tiêu chuẩn TCCS 38:2022/TCĐBVN. Tiêu chuẩn này tính toán kết cấu dựa trên các giá trị mô đun đàn hồi tĩnh, vốn quy ước tải trọng có phương, chiều và độ lớn không đổi trong suốt quá trình tác dụng. Tuy nhiên, thực tế khai thác cho thấy tải trọng bánh xe chạy qua một phân tố kết cấu là một hàm biến thiên theo thời gian có bản chất sóng chu kỳ (dynamic load pulse), đồng thời vật liệu bê tông nhựa (BTN) là vật liệu đàn nhớt tuyến tính (Linear Viscoelasticity - LVE) cực kỳ nhạy cảm với nhiệt độ và thời gian gia tải. Mặc dù tiêu chuẩn cơ sở TCCS 37:2022/TCĐBVN đã từng bước tiếp cận phương pháp cơ học - thực nghiệm (Mechanistic-Empirical Pavement Design Guide - MEPDG), nguồn dữ liệu thực nghiệm về mô đun phức động ($E^$, $G^$), góc trễ pha ($\phi$), hệ số Poát-xông phức động 3D ($\nu^*$) và mô đun đàn hồi động phục hồi ($M_R$) của vật liệu bản địa vẫn còn thiếu hụt trầm trọng.

Nghiên cứu tập trung giải quyết 4 câu hỏi nghiên cứu (Research Questions - RQ) và kiểm định 4 giả thuyết (Hypotheses - H) cụ thể:

  • RQ1 & H1: Đặc tính đàn nhớt của nhựa đường chuyển tiếp sang pha ma tít và bê tông nhựa tuân theo quy luật lưu biến nào? Giả thuyết $H_1$ khẳng định nguyên tắc Tương quan Nhiệt độ - Tần số (Time-Temperature Superposition Principle - TTSP) áp dụng hoàn toàn cho nhựa thông thường và bán toàn phần (semi-TTSP) cho nhựa cải tiến polyme (PMB).
  • RQ2 & H2: Nhiệt độ hóa mềm thực nghiệm (Vòng và Bi) có thể quy đổi sang các tham số lưu biến cơ học cơ bản không? Giả thuyết $H_2$ khẳng định tồn tại tương quan giải tích xác định giữa nhiệt độ hóa mềm với mô đun cắt động $|G^*|$ và góc pha $\phi$.
  • RQ3 & H3: Ứng xử biến dạng không gian 3 chiều (3D) của BTN trong miền thời gian (thí nghiệm từ biến) có thể dự đoán chính xác từ miền tần số (thí nghiệm mô đun phức động) hay không? Giả thuyết $H_3$ cho rằng mô hình Kelvin-Voigt tổng quát 3 chiều (GKV 3D) có khả năng chuyển đổi trực tiếp các tham số biến dạng thể tích.
  • RQ4 & H4: Mức độ sai lệch ứng xử ứng suất - biến dạng giữa mô phỏng đàn nhớt tuyến tính và mô phỏng đàn hồi tĩnh lớn đến đâu? Giả thuyết $H_4$ khẳng định mô phỏng đàn hồi tĩnh bỏ qua hoàn toàn độ trễ pha và phân bố ứng suất cắt trượt nguy hiểm ở nhiệt độ cao.

Khung lý thuyết của luận án được xây dựng dựa trên cơ học môi trường liên tục đàn nhớt tuyến tính, kết hợp các mô hình lưu biến học kinh điển và hiện đại: mô hình Huet-Sayegh, mô hình 2S2P1D (Olard & Di Benedetto), mô hình Kelvin-Voigt tổng quát 3D, phương trình dịch chuyển Williams-Landel-Ferry (WLF) và mô hình phi tuyến phụ thuộc trạng thái ứng suất $K-\theta$ (Hicks & Monismith). Phạm vi thực nghiệm của luận án khảo sát toàn diện 3 chủng loại chất kết dính (nhựa đường 60/70, nhựa đường 35/50, nhựa polyme PMB3), hệ ma tít tương ứng, 3 loại hỗn hợp BTN (BTN C12.5, BTN C19, BTN P12.5) cùng các lớp móng cấp phối đá dăm (CPĐD) và đất nền trong dải nhiệt độ $15^\circ\text{C} - 60^\circ\text{C}$ và tần số $0.01\text{ Hz} - 25\text{ Hz}$.

Literature Review và Positioning

Tổng quan y văn quốc tế cho thấy lịch sử phát triển của cơ học áo đường là quá trình chuyển dịch mô hình (paradigm shift) mạnh mẽ từ đàn hồi tĩnh sang đàn nhớt động lực học:

[Mô hình Đàn hồi Tĩnh] (Boussinesq, 1885; Burmister, 1945; TCCS 38:2022)
[Mô hình Cơ học Lưu biến 1D] (Huet, 1963; Sayegh, 1965; Di Benedetto, 1990)
[Mô hình Đàn nhớt Đa chiều 3D & 2S2P1D] (Olard & Di Benedetto, 2003; Pouget et al., 2010)
[Khung Đánh giá Toàn diện Đa cấp độ & Tải trọng Động] (Bùi Văn Phú, 2024 - Điều kiện VN)

Luận án đã tổng hợp và phân loại các luồng nghiên cứu chính trên thế giới:

  1. Luồng nghiên cứu lưu biến chất kết dính và ma tít: Dẫn đầu bởi Di Benedetto et al. (1990, 2004) [31, 36], Anderson et al. (1994), Airey (2002) [25]. Các tác giả đã khẳng định chất kết dính bitum kiểm soát hành vi phụ thuộc nhiệt độ và thời gian của hỗn hợp. Ma tít đóng vai trò là "chất keo thực thụ" bao bọc cốt liệu thô, trong đó tương tác bề mặt bột khoáng - bitum quyết định khả năng kháng nứt nhiệt và kháng biến dạng vệt bánh xe (hằn lún).
  2. Luồng mô hình hóa đàn nhớt hỗn hợp BTN: Huet (1963) và Sayegh (1965) [98] khởi xướng mô hình kết hợp lò xo và phần tử nhớt parabolic. Sau đó, Olard và Di Benedetto (2003, 2005) [85, 86] hoàn thiện mô hình 2S2P1D (gồm 2 lò xo, 2 phần tử parabolic và 1 piston nhớt tuyến tính), khắc phục triệt để nhược điểm mất kiểm soát của mô hình Huet-Sayegh tại vùng nhiệt độ cao/tần số thấp.
  3. Luồng vật liệu rời và đất nền: Hicks và Monismith (1971) [58], Uzan (1985) [102], Witczak et al. (2004) [93] đã chứng minh cấp phối đá dăm và đất nền có ứng xử phi tuyến phụ thuộc trạng thái ứng suất, đặc trưng bởi mô đun đàn hồi phục hồi $M_R$ biến thiên theo ứng suất khối $\theta$ và ứng suất cắt bát diện $\tau_{oct}$.

Trong y văn tồn tại hai cuộc tranh luận học thuật lớn:

  • Tranh luận 1: Giả định hệ số Poát-xông ($\nu$) hằng số so với hệ số Poát-xông phức động biến thiên $\nu^(\omega, T)$.* Đại đa số các tiêu chuẩn thiết kế (kể cả MEPDG của AASHTO) thường giản lược $\nu$ thành một giá trị cố định (0.30 - 0.35). Tuy nhiên, các nghiên cứu chuyên sâu của Maher et al. (2000) [68] và Schwartz (2002) [99] chỉ ra rằng sự biến thiên của $\nu$ theo nhiệt độ và tần số gây tác động nghiêm trọng đến sai số dự báo lún và nứt mỏi đáy áo đường.
  • Tranh luận 2: Tính duy nhất của nguyên tắc TTSP trên nhựa biến tính polyme. Nhiều học giả cho rằng nhựa polyme (như SBS) làm phá vỡ biểu đồ pha liên tục trong mặt phẳng Black. Di Benedetto và Olard (2005) [86] đề xuất khái niệm "nguyên tắc tương quan bán toàn phần", thừa nhận sự không duy nhất trên đồ thị Black nhưng vẫn cho phép xây dựng đường cong đặc trưng mô đun động $|E^*|$.

So sánh với hai nghiên cứu thực nghiệm quốc tế tiêu biểu:

  • Nghiên cứu của Juliette Blanc et al. (2017) [40] (Pháp): Sử dụng mạng lưới cảm biến đo đạc biến dạng đáy áo đường thực tế trong nhiều năm, chứng minh rằng biến dạng thẳng đứng tại đỉnh lớp đất nền và biến dạng kéo ngang đáy lớp mặt do các trục xe tác dụng có tính chất trễ pha rõ rệt và độ lớn biến dạng phụ thuộc nghiêm ngặt vào vận tốc xe chạy.
  • Nghiên cứu của Lofft / Loulizi et al. (2005) [65] (Mỹ): So sánh biến dạng ngang $\epsilon_{yy}$ giữa mô phỏng đàn hồi tuyến tính, mô phỏng đàn nhớt và số liệu đo đạc thực tế tại đường thử nghiệm thông minh Virginia Smart Road. Kết quả cho thấy mô hình đàn hồi tạo ra đường đáp ứng biến dạng đối xứng không thực tế, trong khi mô hình đàn nhớt tái hiện chính xác hình dạng bất đối xứng và độ phục hồi biến dạng trễ của kết cấu.

Vị trí của luận án Bùi Văn Phú (2024) được định vị là công trình tiên phong tại Việt Nam tích hợp chuỗi chuyển tiếp đa cấp độ (multi-scale): từ bitum $\rightarrow$ ma tít $\rightarrow$ hỗn hợp BTN (1D & 3D) $\rightarrow$ kết cấu áo đường nhiều lớp chịu tải trọng xe chạy động học, xây dựng cơ sở thực nghiệm vững chắc phù hợp điều kiện vật liệu và khí hậu nhiệt đới Việt Nam.

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

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

Luận án mang lại những đóng góp nền tảng cho lý thuyết lưu biến và cơ học kết cấu áo đường:

  1. Mở rộng và ứng dụng mô hình 2S2P1D trong dự đoán đa pha: Luận án đã mở rộng ứng dụng mô hình 2S2P1D của Olard & Di Benedetto (2003) để dự đoán chính xác mô đun cắt động $|G^*|$ của ma tít từ các thông số đàn nhớt của nhựa đường nguyên bản kết hợp với hàm lượng thể tích bột khoáng $V(%)$, thông qua phương trình hàm mũ hồi quy giữa hệ số dịch chuyển logarit $\log(\alpha)$ và $V(%)$.
  2. Thiết lập mối quan hệ giải tích giữa chỉ tiêu kinh nghiệm và cơ học lưu biến: Luận án đã giải mã thành công mối quan hệ giữa nhiệt độ hóa mềm kinh nghiệm (phương pháp Vòng và Bi) với các thông số vật lý cơ bản thu được từ thiết bị đo lưu biến cắt động (DSR) bao gồm $|G^*|$, góc pha $\phi$ và tham số độ cứng nhiệt $K$. Luận án mở rộng phương pháp xác định nhiệt độ hóa mềm tương đương $T_V$ của Alisov và đề xuất quy trình xác định nhiệt độ hóa mềm tương đương $T_K$ cho cả nhựa thông thường và nhựa polyme PMB3.
  3. Chuyển đổi lý thuyết 3D giữa miền tần số và miền thời gian: Đóng góp lý thuyết nổi bật là việc chứng minh sự tương đương giữa thí nghiệm mô đun phức động 3D ($|E^|, |\nu^|$) và thí nghiệm từ biến 3D ($J(t), \nu(t)$) thông qua mô hình Kelvin-Voigt tổng quát 3D (GKV 3D), cho phép dự đoán chính xác sự phát triển của biến dạng dọc trục và biến dạng nở hông theo thời gian mà không cần tiến hành các thí nghiệm từ biến phức tạp.
                  CHUYỂN TIẾP ĐA CẤP ĐỘ CỦA VẬT LIỆU

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

Khung phân tích của luận án kết nối chặt chẽ ba lý thuyết trụ cột:

  • Lý thuyết Đàn nhớt tuyến tính môi trường liên tục (Linear Viscoelasticity Continuum Theory): Xác định quan hệ ứng suất - biến dạng dưới dạng tích phân chập Volterra và hàm phức số: $$\sigma^(t) = E^(\omega)\cdot\epsilon^(t) = (|E^|\cos\phi + j|E^|\sin\phi)\epsilon^(t)$$
  • Nguyên tắc Tương quan Nhiệt độ - Tần số (TTSP): Chuẩn hóa dữ liệu đo ở các nhiệt độ khác nhau về một đường cong đặc trưng (master curve) duy nhất tại nhiệt độ tham chiếu ($T_{ref} = 25^\circ\text{C}$ cho bitum/ma tít và $T_{ref} = 40^\circ\text{C}$ cho BTN) qua hàm dịch chuyển WLF: $$\log(a_T) = \frac{-C_1(T - T_{ref})}{C_2 + T - T_{ref}}$$
  • Lý thuyết Môi trường Rời phi tuyến ứng suất: Mô tả lớp cấp phối đá dăm và đất nền bằng mô hình $K-\theta$ của Hicks & Monismith: $$M_R = k_1(\theta)^{k_2} = k_1(\sigma_1 + \sigma_2 + \sigma_3)^{k_2}$$

Điều kiện biên xác định: Miền ứng xử đàn nhớt tuyến tính của BTN được khống chế nghiêm ngặt ở mức biến dạng cực nhỏ $\epsilon < 10^{-4}\text{ m/m}$ (dưới 100 microstrain) để đảm bảo cấu trúc khung cốt liệu không bị suy thoái phi tuyến hoặc phá hủy mỏi trong suốt quá trình thử nghiệm.

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

Thiết kế nghiên cứu

Luận án tuân thủ thế giới quan thực chứng (positivism) kết hợp chủ nghĩa hiện thực phản biện (critical realism). Thiết kế nghiên cứu thực nghiệm đa mức độ (multi-scale experimental design) kết hợp mô phỏng số giải tích động lực học kết cấu:

Cấp độ khảo sát Đối tượng vật liệu Thiết bị & Tiêu chuẩn Thông số đo đạc chính
Cấp 1: Vi mô Nhựa đường (60/70, 35/50, PMB3) DSR Rheotest RN 4.3 (ASTM D7175) $|G^*|, \phi, T_V, T_K$, Master curve ($T_{ref}=25^\circ\text{C}$)
Cấp 2: Trung gian Ma tít át phan (Bitum + Bột khoáng) DSR Rheotest RN 4.3 (Cắt góc / Kéo nén) $|G^*|_{mastic}, \phi, a_T$, Hệ số tương quan $\alpha$
Cấp 3: Vĩ mô Hỗn hợp BTN (C12.5, C19, P12.5) Cooper Test cải tiến 3D LVDTs $|E^|, \phi_E, |\nu^|, \phi_\nu, J(t), \nu(t)$ ($15, 30, 40, 60^\circ\text{C}$)
Cấp 4: Nền móng Cấp phối đá dăm & Đất nền Hệ thống nén 3 trục chu kỳ (AASHTO T307) Mô đun đàn hồi phục hồi $M_R$, hệ số hồi quy $k_1, k_2$
Cấp 5: Kết cấu Áo đường mềm điển hình Việt Nam Phần mềm chuyên dụng Viscoroute 2.0 Trường ứng suất $\sigma_{xx}, \tau_{xz}$, độ võng, biến dạng $\epsilon_{yy}$

Quy trình nghiên cứu rigorous

Tác giả đã thực hiện các cải tiến kỹ thuật quan trọng trên thiết bị thí nghiệm gia tải Cooper để thực hiện thành công phép đo đồng thời biến dạng dọc trục và biến dạng nở hông:

  • Cải tiến đo đạc 3D: Lắp đặt hệ thống vòng đo cảm biến biến dạng ngoài (external radial LVDT array) gồm các đầu đo dịch chuyển vi sai đối xứng $120^\circ$ quanh thân mẫu hình trụ ($\phi 100\text{mm} \times h 150\text{mm}$), cho phép ghi nhận biến dạng nở hông cực nhỏ với độ phân giải micromet.
  • Giao thức kiểm định AASHTO T307: Mẫu đất nền chịu quá trình tiền gia tải (conditioning) 1000 chu kỳ ở áp lực buồng $\sigma_3 = 41.4\text{ kPa}$ và ứng suất lệch $\sigma_d = 24.8\text{ kPa}$ nhằm triệt tiêu biến dạng dư ban đầu. Tiếp đó, mẫu được thí nghiệm qua 3 cấp áp lực buồng ($\sigma_3 = 41.4, 27.6, 13.8\text{ kPa}$), mỗi cấp buồng chịu 5 cấp ứng suất lệch ($\sigma_d = 12.4, 24.8, 37.3, 49.7, 62.0\text{ kPa}$), mỗi cấp 100 chu kỳ tải dạng haversine ($0.1\text{s}$ tác dụng, $0.9\text{s}$ nghỉ). Tổng cộng 16 giai đoạn gia tải cho một mẫu thí nghiệm.
XUNG TẢI TRỌNG HAIVERSINE TRONG THÍ NGHIỆM M_R (AASHTO T307)
 Ứng suất σ
        0.1s    1.0s      1.1s    2.0s

Data và phân tích

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

  • Kiểm định Goodness-of-Fit: Đánh giá độ phù hợp của mô hình 2S2P1D và GKV thông qua hệ số xác định hiệu chỉnh $R^2$, sai số bình phương trung bình chuẩn hóa ($NRMSE$) và kiểm định sai số tỉ đối Se/Sy theo tiêu chuẩn MEPDG.
  • Phân tích phương sai (ANOVA): Đánh giá mức độ ảnh hưởng và ý nghĩa thống kê của các nhân tố nhiệt độ, tần số và loại nhựa đường đối với $|E^|$, $|G^|$ và hệ số Poát-xông $|\nu^*|$ ở mức ý nghĩa $\alpha = 0.05$ ($p < 0.001$).
  • Mô phỏng cơ học kết cấu: Ứng dụng phần mềm phần tử hữu hạn chuyển vị đàn nhớt nhiều lớp Viscoroute 2.0 do Viện LCPC/IFSTTAR (Pháp) phát triển, mô phỏng tải trọng trục đơn bánh đôi ($P = 100\text{ kN}$, áp lực vệt bánh $p = 0.707\text{ MPa}$) di động với các dải vận tốc $V = 20, 40, 60, 80\text{ km/h}$ trên mặt đường ở các dải nhiệt độ phân tầng thực tế từ $15^\circ\text{C}$ đến $60^\circ\text{C}$.

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

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

  1. Tính biến thiên của hệ số Poát-xông phức động 3D ($|\nu^*|$) và hệ số Poát-xông từ biến ($\nu(t)$): Khác biệt hoàn toàn với giả định truyền thống xem $\nu$ là hằng số, thực nghiệm 3D chứng minh $|\nu^|$ của BTN C12.5, BTN C19 và BTN P12.5 tăng mạnh khi nhiệt độ tăng và tần số giảm. Cụ thể, $|\nu^|$ biến thiên từ $0.22 - 0.25$ ở nhiệt độ thấp ($15^\circ\text{C}$, tần số cao $10\text{ Hz}$) lên tới $0.45 - 0.48$ ở nhiệt độ cao ($60^\circ\text{C}$, tần số thấp $0.01\text{ Hz}$). Trong thí nghiệm từ biến, $\nu(t)$ có xu hướng giảm nhẹ trong giai đoạn đầu gia tải và duy trì ổn định khi dỡ tải.

  2. Sự tập trung ứng suất cắt trượt cực đại ($\tau_{xz}$) tại đáy lớp mặt BTN ở nhiệt độ cao: Tại nhiệt độ mặt đường mùa hè ($60^\circ\text{C}$), ứng suất cắt trượt $\tau_{xz}$ tại đáy lớp BTN C19 tăng vọt lên mức nghiêm trọng, vượt qua khả năng chịu cắt của vật liệu. Đây là phát hiện then chốt giải thích nguyên nhân gốc rễ gây ra hư hỏng trượt tầng và hằn lún vệt bánh xe tại các tuyến quốc lộ và cao tốc ở Việt Nam.

PHÂN BỐ ỨNG SUẤT CẮT TRƯỢT τ_xz THEO CHIỀU SÂU Ở NHIỆT ĐỘ 60°C
 Chiều sâu z (cm)

  1. Hiện tượng bất đối xứng và độ trễ pha của trường biến dạng ngang $\epsilon_{yy}$: Kết quả mô phỏng đàn nhớt bằng Viscoroute 2.0 chỉ ra rằng sóng biến dạng tại đáy lớp BTN khi bánh xe chạy qua hoàn toàn bất đối xứng: biến dạng nén xuất hiện trước khi bánh xe tới đỉnh, chuyển sang biến dạng kéo cực đại khi bánh xe nằm ngay trên điểm khảo sát, và duy trì biến dạng dư trễ pha kéo dài sau khi bánh xe đã đi qua. Tính toán đàn hồi cổ điển hoàn toàn thất bại trong việc mô tả hiện tượng vật lý này.

  2. Quy luật chuyển tiếp ma tít và độ tin cậy của mô hình 2S2P1D: Hệ số dịch chuyển đường cong đặc trưng $a_T$ của ma tít hoàn toàn đồng nhất với $a_T$ của nhựa đường gốc. Mức độ tăng cứng của ma tít tỉ lệ thuận với hàm lượng thể tích bột khoáng $V(%)$ theo quan hệ tuyến tính chặt chẽ $\log(\alpha) = A\cdot V(%) + B$ với hệ số tương quan $R^2 > 0.98$, cho phép dự báo trực tiếp mô đun cắt $|G^*|$ của hỗn hợp ma tít mà không cần thí nghiệm DSR lặp lại.

  3. Nhiệt độ hóa mềm tương đương $T_K$ và $T_V$: Thiết lập mối tương quan hồi quy có độ tin cậy rất cao ($R^2 > 0.96$) giữa nhiệt độ hóa mềm thực nghiệm với các chỉ tiêu lưu biến học DSR, chứng minh rằng tại nhiệt độ hóa mềm, góc lệch pha $\phi$ của bitum nguyên gốc 60/70 và 35/50 đạt trạng thái chuyển pha nhớt xấp xỉ $80^\circ - 85^\circ$, trong khi bitum polyme PMB3 giữ được độ đàn hồi vượt trội nhờ mạng lưới polyme SBS ($\phi \approx 65^\circ - 70^\circ$).

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp bức tranh toàn diện về cơ học đàn nhớt tuyến tính 3D cho vật liệu đường bộ Việt Nam; hoàn thiện phương pháp chuyển đổi giải tích giữa miền tần số ($E^*$) và miền thời gian ($J(t)$).
  • Về mặt phương pháp luận: Thiết lập quy trình chuẩn thí nghiệm mô đun phức động 3D trên thiết bị Cooper và giao thức nén 3 trục chu kỳ xác định $M_R$ theo AASHTO T307, có thể chuyển giao cho các phòng thí nghiệm LAS-XD trọng điểm.
  • Về mặt thực tiễn kỹ thuật: Cảnh báo nguy cơ phá hoại do cắt trượt ở tầng đáy lớp BTN mặt đường tại vùng khí hậu nắng nóng; kiến nghị sử dụng BTN polyme (PMB) cho lớp mặt trên và tối ưu hóa cấp phối đá dăm kháng trượt cho lớp dưới.
  • Về mặt chính sách: Đóng vai trò là luận cứ khoa học thực nghiệm trực tiếp thúc đẩy Bộ Giao thông Vận tải hoàn thiện lộ trình chuyển đổi từ tiêu chuẩn đàn hồi tĩnh TCCS 38:2022 sang tiêu chuẩn cơ học - thực nghiệm hiện đại TCCS 37:2022.

Limitations và Future Research

Luận án thẳng thắn thừa nhận các giới hạn nghiên cứu:

  1. Phạm vi vật liệu: Chưa khảo sát hết các loại phụ gia tái sinh nhựa đường (RAP), nhựa đường bọt hay bê tông nhựa ấm (WMA).
  2. Miền biến dạng: Giới hạn trong miền đàn nhớt tuyến tính ($\epsilon < 10^{-4}\text{ m/m}$), chưa mô hình hóa phi tuyến phá hủy mỏi lớn (fatigue cracking damage) và biến dạng dẻo tích lũy vĩnh cửu (rutting creep damage) ở giai đoạn phá hoại bậc 3.
  3. Mô hình hóa nền đất: Lớp cấp phối đá dăm và đất nền được gán mô đun tương đương phân lớp theo trạng thái ứng suất tĩnh, chưa tích hợp phần tử hữu hạn phi tuyến 3D hoàn toàn (Full 3D Non-linear Dynamic FEM).

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Mở rộng mô hình 2S2P1D sang miền phi tuyến đàn - nhớt - dẻo (Visco-Elasto-Plastic DVEP Model).
  • Nghiên cứu ứng xử tự hồi phục (self-healing) của nhựa đường và ma tít dưới tác động của nhiệt độ và chu kỳ nghỉ giữa các lượt xe.
  • Ứng dụng công nghệ cảm biến sợi quang FBG chôn ngầm trong áo đường thực tế để hiệu chỉnh mô hình số liệu đo đạc hiện trường tại Việt Nam.

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

  • Học thuật: Công trình tạo nền móng vững chắc cho các nghiên cứu sinh, học viên cao học chuyên ngành Cầu đường; tiềm năng trích dẫn cao trong các tạp chí thuộc danh mục ISI/Scopus (Q1/Q2) như Construction and Building Materials, Materials and Structures, Road Materials and Pavement Design.
  • Chuyển đổi ngành công nghiệp: Cung cấp cơ sở dữ liệu vật liệu đầu vào chính xác cho các đơn vị tư vấn thiết kế giao thông hàng đầu (TEDI, TEDIS, các Viện KHCN GTVT), giúp tối ưu hóa chiều dày kết cấu, giảm thiểu tình trạng thiết kế thừa hoặc phá hoại sớm.
  • Lợi ích kinh tế - xã hội: Giảm thiểu hư hỏng hằn lún vệt bánh xe giúp tiết kiệm hàng trăm tỷ đồng ngân sách bảo trì mặt đường hàng năm, nâng cao an toàn giao thông và giảm thời gian ách tắc trên các hành lang kinh tế trọng điểm.

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

  • Nghiên cứu sinh & Giới học thuật: Tiếp cận bộ thông số thực nghiệm toàn diện về mô đun động 1D/3D, mô hình toán 2S2P1D và quy trình chuyển đổi miền tần số - thời gian.
  • Kỹ sư thiết kế & Tư vấn giao thông: Sở hữu phương pháp luận và dữ liệu thực tế để áp dụng tiêu chuẩn tiên tiến TCCS 37:2022/TCĐBVN và AASHTO MEPDG vào dự án đường cao tốc.
  • Cơ quan quản lý nhà nước (Cục Đường bộ VN, Bộ GTVT): Có căn cứ khoa học chính xác để ban hành các định mức kinh tế kỹ thuật, tiêu chuẩn vật liệu và chỉ dẫn kỹ thuật thi công kiểm soát chất lượng mặt đường.
  • Doanh nghiệp sản xuất nhựa đường & Bê tông nhựa: Hiểu rõ cơ chế tương tác vi mô bột khoáng - bitum để thiết kế cấp phối ma tít tối ưu.

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 độc đáo nhất là việc tích hợp thành công mô hình lưu biến 2S2P1D (Olard & Di Benedetto) và mô hình Kelvin-Voigt tổng quát 3D (GKV 3D) để thiết lập chuỗi chuyển tiếp giải tích liên tục từ tính chất đàn nhớt vi mô của chất kết dính (nhựa đường) sang ma tít, và từ mô đun phức động đa chiều 3D ($|E^|, |\nu^|$) sang hàm từ biến theo thời gian ($J(t), \nu(t)$) của bê tông nhựa trong điều kiện vật liệu Việt Nam.

2. Tính đổi mới trong phương pháp nghiên cứu khi so sánh với các công trình trước đây?

So với các nghiên cứu trước đây tại Việt Nam (vốn chỉ đo mô đun kéo nén tĩnh 1D hoặc mô đun đàn hồi uốn theo TCCS 38:2022), luận án đã:

  • Cải tiến thành công hệ thống gá đo và cảm biến dịch chuyển vi sai ngoài (3D Radial LVDTs) trên máy nén Cooper để thu thập đồng thời biến dạng thể tích hông, xác định hệ số Poát-xông phức động $\nu^*(\omega, T)$.
  • Áp dụng triệt để quy trình nén 3 trục chu kỳ 16 cấp gia tải theo AASHTO T307 trên đất nền và CPĐD thay vì thí nghiệm ép lún tấm ép tĩnh hoặc CBR truyền thống.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì?

Phát hiện bất ngờ nhất là tại nhiệt độ cao mùa hè ($60^\circ\text{C}$), ứng suất kéo uốn tại đáy lớp BTN giảm mạnh do vật liệu bị mềm hóa (chuyển sang trạng thái nhớt chảy), nhưng ngược lại ứng suất cắt trượt $\tau_{xz}$ lại tập trung cực đại tại mặt tiếp xúc giữa lớp mặt dưới BTN C19 và lớp móng CPĐD. Điều này chứng minh nứt mỏi đáy áo đường không phải là nguy cơ lớn nhất trong mùa hè, mà chính hiện tượng trượt tầng do ứng suất cắt kết hợp suy giảm mô đun mới là nguyên nhân kích hoạt hằn lún vệt bánh xe.

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

Có. Luận án mô tả chi tiết từng bước chuẩn bị mẫu hình trụ Marshall/Superpave khoan cắt, tỷ lệ bột khoáng cấp phối, thông số cài đặt buồng nhiệt, tần số quét ($0.01 - 25\text{ Hz}$), biên độ biến dạng khống chế ($\epsilon \le 50\mu\epsilon$), quy trình tiền gia tải buồng nén 3 trục, và mã lệnh thiết lập phần mềm Viscoroute 2.0.

5. Định hướng chương trình nghiên cứu 10 năm tới được vạch ra như thế nào?

Xây dựng cơ sở dữ liệu số quốc gia về mô đun phức động ($E^*$) cho toàn bộ các mỏ đá và nguồn nhựa đường tại 3 miền Bắc - Trung - Nam; tích hợp mô hình phá hủy cơ - nhiệt - ẩm liên hợp (coupled thermo-hydro-mechanical damage model); và ứng dụng trí tuệ nhân tạo (Machine Learning/PINNs) để dự báo suy thoái mặt đường thời gian thực.

Kết luận

  1. Chuẩn hóa cơ sở thực nghiệm lưu biến: Luận án đã xây dựng bộ cơ sở dữ liệu thực nghiệm hoàn chỉnh đầu tiên tại Việt Nam về mô đun cắt động $|G^|$, mô đun phức động 1D/3D $|E^|$, góc lệch pha $\phi$, hệ số Poát-xông động $|\nu^*|$ và mô đun phục hồi $M_R$ cho các tầng vật liệu áo đường mềm.
  2. Xác thực nguyên tắc TTSP và mô hình 2S2P1D: Chứng minh tính quy luật của nguyên tắc Tương quan Nhiệt độ - Tần số trên đường cong đặc trưng master curve, khẳng định sự phù hợp tuyệt đối của mô hình 2S2P1D và mô hình GKV 3D trong việc mô phỏng ứng xử vật liệu.
  3. Giải mã nhiệt độ hóa mềm tương đương: Đề xuất thành công công thức hồi quy tương quan chuyển đổi giữa chỉ tiêu kinh nghiệm Vòng và Bi với các tham số độ cứng lưu biến học DSR ($T_V, T_K$), nâng tầm giá trị khoa học cho các phép thử truyền thống.
  4. Đột phá trong mô phỏng kết cấu động học: Chứng minh tính vượt trội của mô hình đàn nhớt tuyến tính nhiều lớp (Viscoroute 2.0) so với mô hình đàn hồi tĩnh cổ điển (TCCS 38:2022), vạch rõ cơ chế phá hoại do tích tụ biến dạng trễ pha và ứng suất cắt trượt nguy hiểm ở nhiệt độ cao ($60^\circ\text{C}$).
  5. Mở ra 3 hướng nghiên cứu mới: Mở đường cho các nghiên cứu sâu về mô hình đàn nhớt phi tuyến 3D, cơ học phá hủy mỏi liên hợp, và tối ưu hóa thiết kế áo đường cao tốc thông minh thích ứng với biến đổi khí hậu tại Việt Nam.