Tổng quan luận án

Luận án tiến sĩ thuộc chuyên ngành Kỹ thuật xây dựng công trình giao thông (kỹ thuật đường bộ), tiếp cận theo hướng thực nghiệm phòng thí nghiệm kết hợp mô hình hóa toán học - cơ học và mô phỏng cơ học thực nghiệm (Mechanistic - Empirical, viết tắt là ME).

Về tính cấp thiết và khoảng trống nghiên cứu, văn bản luận án nêu rõ phương pháp thiết kế kết cấu mặt đường mềm tại Việt Nam hiện vẫn dựa trên các quy trình truyền thống sử dụng mô đun đàn hồi tĩnh hoặc mô đun đàn hồi động đo theo phương pháp kéo gián tiếp. Cách tiếp cận này chưa phản ánh đúng bản chất làm việc thực tế của vật liệu bitum và bê tông nhựa (BTN) là vật liệu có tính đàn nhớt biến thiên theo nhiệt độ và thời gian tác dụng của tải trọng. Sự thiếu hụt này dẫn tới việc xác định chiều dày kết cấu chưa tối ưu, mặt đường xuất hiện hư hỏng sớm, đặc biệt là hiện tượng hằn lún vệt bánh xe và nứt mỏi, làm suy giảm tuổi thọ công trình. Mặc dù Đề tài cấp Bộ Giao thông vận tải năm 2017 (mã số DT174055) đã kết luận việc áp dụng phương pháp thiết kế mặt đường cơ học thực nghiệm (ME) vào điều kiện Việt Nam là khả thi, song để vận hành được mô hình ME bắt buộc phải "địa phương hóa các thông số đầu vào". Hai chỉ tiêu cốt lõi của phương pháp này là mô đun cắt động của bitum ($|G^|$) và mô đun động của bê tông nhựa ($|E^|$). Tại thời điểm thực hiện đề tài, các nghiên cứu trong nước về $|G^|$, góc pha ($\delta_b$) và $|E^|$ còn rất hạn chế, chưa có cơ sở dữ liệu thực nghiệm toàn diện và chưa từng có công trình nào tiến hành hiệu chỉnh các hệ số dự báo $|E^*|$ theo điều kiện vật liệu cụ thể tại Việt Nam.

Mục tiêu nghiên cứu tổng quát của luận án là xác định mô hình dự báo $|E^*|$ theo tính chất của loại bitum sử dụng phù hợp với điều kiện vật liệu của Việt Nam. Để đạt được mục tiêu này, luận án xác định 5 mục tiêu cụ thể được đánh số theo văn bản:

  1. Hiểu rõ tính chất của các vật liệu thành phần của hỗn hợp BTN gồm mô đun cắt động $|G^|$, góc pha ($\delta_b$), mô đun động của BTN ($|E^|$).
  2. Hiểu rõ được các yếu tố ảnh hưởng tới mô đun cắt động của bitum ($|G^|$), góc pha ($\delta_b$), mô đun động của BTN ($|E^|$) và yếu tố nào có ảnh hưởng nhiều nhất tới giá trị mô đun động của BTN.
  3. Hiểu rõ được phương pháp xây dựng đường cong chủ của $|G^|$, góc pha ($\delta_b$), mô đun động của BTN ($|E^|$) và việc đánh giá mức độ phù hợp của các mô hình xây dựng đường cong chủ ($|G^|$, góc pha ($\delta_b$), $|E^|$) cho các loại bitum và BTN.
  4. Hiểu rõ được các mô hình dự báo $|E^*|$ của Hoa Kỳ theo các mô hình Witczak, mô hình Hirsch và cách đánh giá mức độ phù hợp của các mô hình dự báo này theo tiêu chuẩn thống kê.
  5. Xác định được mô hình dự báo $|E^|$ trong ba mô hình của Hoa Kỳ (các mô hình Witczak và mô hình Hirsch) có khả năng dự báo $|E^|$ với độ chính xác cao nhất theo tiêu chuẩn thống kê theo điều kiện vật liệu của Việt Nam.

Đối tượng và phạm vi nghiên cứu của luận án được giới hạn cụ thể như sau:

  • Đối tượng nghiên cứu: Mô đun cắt động $|G^|$ và góc pha ($\delta_b$) của chất kết dính bitum; mô đun động $|E^|$ của hỗn hợp bê tông nhựa chặt (BTNC); các mối quan hệ tương quan thực nghiệm dự báo $|E^*|$; đặc trưng khai thác của kết cấu mặt đường mềm phân tích theo phương pháp ME.
  • Phạm vi vật liệu: Nghiên cứu 4 loại bitum gồm bitum thông thường (60/70, 35/50, 40/50) và bitum cải tiến polymer (PMBIII); 2 loại hỗn hợp bê tông nhựa chặt gồm BTNC 12.5 và BTNC 19; cốt liệu khoáng sản xuất từ 3 nguồn đá gồm đá bazan, đá granit, đá vôi; bột khoáng sản xuất từ đá gốc đá vôi. Các yêu cầu kỹ thuật và cấp phối tuân thủ theo TCVN 8819:2011.
  • Phạm vi không gian và thời gian: Điều kiện khí hậu, nhiệt độ tính toán lấy theo khu vực miền Bắc Việt Nam; dữ liệu giao thông và kết cấu phục vụ mô phỏng kiểm chứng thu thập từ một dự án giao thông thực tế tại Việt Nam.

Tổng quan tài liệu và vị trí của luận án

Luận án hệ thống hóa các công trình nghiên cứu trên thế giới và trong nước liên quan đến chỉ tiêu đàn nhớt của bitum và mô đun động của bê tông nhựa:

  • Các hướng nghiên cứu quốc tế về tính chất đàn nhớt và độ cứng của bitum: Van der Poel (1954) đưa ra khái niệm mô đun độ cứng của bitum $S_{t,T}$ theo thời gian tác dụng tải trọng $t$ và nhiệt độ $T$, xây dựng toán đồ (nomograph) dự báo độ cứng dựa trên độ kim lún và điểm hóa mềm. Chương trình Nghiên cứu Chiến lược Đường bộ Hoa Kỳ (SHRP) xây dựng hệ thống phân loại bitum theo cấp đặc tính khai thác PG (AASHTO M320), trong đó $|G^|$ và góc pha $\delta_b$ đo bằng thiết bị DSR là các chỉ tiêu then chốt. NCHRP (2002) thiết lập mối quan hệ thực nghiệm chuyển đổi $|G^|$ và $\delta_b$ tại tần số 10 rad/s sang độ nhớt bitum để làm đầu vào cho mô hình dự báo. Javed Bari và Matthew W. Witczak (2006) phát triển mô hình dự báo góc pha và $|G^*|$ cho cả bitum thông thường lẫn bitum cải tiến polymer.
  • Các hướng nghiên cứu quốc tế về mô hình dự báo mô đun động $|E^|$: Tập đoàn Shell phát triển phần mềm BANDS 2.0 và toán đồ dự báo độ cứng BTN dựa trên độ cứng bitum phục hồi, thể tích bitum $V_b$ và thể tích cốt liệu $V_g$, nhưng phương pháp này bị hạn chế khi độ cứng bitum nhỏ hơn 5 MPa. Viện Asphalt Hoa Kỳ (Hwang và Witczak, 1979) bổ sung thêm hàm lượng lọt sàng số 200 ($P_{200}$), tần số tải trọng $f$ và độ nhớt bitum tại $70^\circ\text{F}$. Nhóm nghiên cứu của Witczak sau đó phát triển mô hình Witczak truyền thống (Original Witczak) dựa trên 2.750 bộ dữ liệu thí nghiệm của 205 hỗn hợp, tiếp theo là mô hình Witczak cải tiến (Modified Witczak, Bari và Witczak, 2006) dựa trên 7.400 giá trị $|E^|$ của 346 hỗn hợp BTN, sử dụng trực tiếp $|G^*|$ và $\delta_b$. Christensen và cộng sự phát triển mô hình Hirsch bán thực nghiệm dựa trên lý thuyết vật liệu tổ hợp phần tử nối tiếp và song song (với 18 hỗn hợp BTN, 8 loại bitum, 206 điểm dữ liệu). Al-Khateeb và cộng sự đề xuất mô hình định luật hỗn hợp tương đương (Law of mixtures parallel model). Ahmad M. Abu Abdo (Đại học Dhofar, Oman; 2009) áp dụng phân tích độ nhạy cho mô hình Idaho. Tại Úc và Hàn Quốc, các nghiên cứu thực nghiệm cũng được triển khai nhằm đánh giá và hiệu chỉnh các mô hình Witczak và Hirsch theo nguồn vật liệu bản địa (chẳng hạn Hàn Quốc nghiên cứu trên 3 cấp phối chặt và cấp phối SMA với bitum PG58-22 và PG64-22).
  • Các nghiên cứu tại Việt Nam: Đề tài cấp Bộ năm 2015 (mã số DT154015) bước đầu nghiên cứu phân loại bitum theo phương pháp Superpave. Đề tài cấp Bộ năm 2017 (mã số DT174055) đánh giá khả năng áp dụng phương pháp ME nhưng chỉ dừng lại ở phân tích khung lý thuyết và khuyến nghị địa phương hóa thông số. Một số nghiên cứu của TS. Nguyễn Mai Lân sử dụng thiết bị DMA khảo sát mô đun phức của bitum 60/70; Viện Khoa học và Công nghệ GTVT thí nghiệm phân loại AASHTO M320 cho các mác nhựa 40/50, 60/70, PMBI, PMBII, PMBIII.

Khoảng trống nghiên cứu mà luận án lựa chọn giải quyết: Các công trình trong nước trước đó chỉ thực hiện thí nghiệm rời rạc ở số lượng mẫu nhỏ, phục vụ mục đích phân loại mác nhựa hoặc làm luận văn thạc sĩ, chưa có nghiên cứu hệ thống nào xây dựng đường cong chủ đàn nhớt toàn diện cho các loại bitum ở Việt Nam bằng mô hình cơ học tiên tiến, chưa áp dụng quy hoạch thực nghiệm Taguchi để tối ưu hóa việc xác định $|E^|$, và chưa từng thực hiện phân tích độ nhạy thống kê kết hợp hiệu chỉnh bộ hệ số cho các mô hình dự báo $|E^|$ (Witczak ban đầu, Witczak cải tiến, Hirsch) theo nguồn vật liệu cát, đá, bitum tại Việt Nam.

Cơ sở lý thuyết và phương pháp nghiên cứu

Cơ sở lý thuyết và khung phân tích

  • Lý thuyết đàn nhớt tuyến tính (Linear Viscoelasticity): Ứng xử cơ học của bitum và BTN thể hiện tính chất trễ pha giữa ứng suất hình sin tác dụng và biến dạng sinh ra. Mô đun phức cắt của bitum được biểu diễn qua hai thành phần: mô đun dự trữ đàn hồi ($G'$) và mô đun tổn thất nhớt ($G''$), với $|G^| = \sqrt{(G')^2 + (G'')^2}$. Tương tự, mô đun phức của bê tông nhựa $E^$ gồm phần thực $E_1$ và phần ảo $E_2$, với giá trị tuyệt đối $|E^*| = \sigma_0 / \varepsilon_0$.
  • Các mô hình xây dựng đường cong chủ (Master curve): Luận án xem xét các mô hình cơ học kinh điển gồm Maxwell, Voigt-Kelvin, Burgers, Huet, Huet-Sayegh, và tập trung vào mô hình cơ học 2S2P1D (gồm 2 lò xo đàn hồi tuyến tính, 2 phần tử biến dạng parabol và 1 giảm chấn tuyến tính).
  • Nguyên lý tương đương thời gian - nhiệt độ (TTSP): Chuyển đổi các tần số đo ở nhiệt độ khác nhau về tần số quy đổi (reduced frequency, $f_r$) tại nhiệt độ tham chiếu $T_{ref}$ qua hệ số dịch chuyển $a_T$: $\log(f_r) = \log(f) + \log(a_T)$. Hệ số dịch chuyển $a_T$ được xác định theo phương trình Williams - Landel - Ferry (WLF): $$\log(a_T) = \frac{-C_1 (T - T_{ref})}{C_2 + (T - T_{ref})}$$ Các hằng số $C_1, C_2$ được xác định bằng phương pháp bình phương tối thiểu trên MS Excel.

Phương pháp nghiên cứu thực nghiệm và công cụ phân tích

  • Thiết bị và quy trình thí nghiệm bitum:
    • Thiết bị cắt động lưu biến DSR (Dynamic Shear Rheometer) theo tiêu chuẩn AASHTO T 315-12, sử dụng tần số góc 10 rad/s (tương đương 1,59 Hz), sử dụng hệ đĩa song song đường kính 25 mm (cho nhiệt độ $\ge 30^\circ\text{C}$) và đĩa 8 mm (cho nhiệt độ $\le 30^\circ\text{C}$).
    • Thiết bị phân tích cơ học động DMA (Dynamic Mechanical Analyzer) MetraviB (IFSTTAR) với đầu đo lực đến 450 N, đầu đo chuyển vị đến 500 $\mu\text{m}$, tần số làm việc từ 1 Hz đến 200 Hz theo tiêu chuẩn Pháp XP-T-66-065. Phương thức tác dụng lực gồm: kéo nén (K/N) cho nhiệt độ $\le 20^\circ\text{C}$ và cắt góc (C/G) cho nhiệt độ $\ge 20^\circ\text{C}$.
  • Thiết bị và quy trình thí nghiệm bê tông nhựa: Thiết bị Cooper của hãng Keith Cooper & Steven Brown (Đại học Nottingham) đặt tại Phòng thí nghiệm vật liệu Trường Đại học Giao thông vận tải, đáp ứng tiêu chuẩn AASHTO TP-62 và các tiêu chuẩn châu Âu EN 12697-24, EN 12697-25, EN 12697-26. Dải tần số thí nghiệm mặc định gồm: 0.5 Hz, 1 Hz, 5 Hz, 10 Hz và 25 Hz.
  • Kỹ thuật quy hoạch thực nghiệm: Áp dụng phương pháp Taguchi với mảng trực giao $L_8(2^7)$ nhằm tối ưu hóa số lượng tổ mẫu chế tạo và thí nghiệm mà vẫn bảo đảm độ tin cậy khoa học và phân tích đa biến cho hỗn hợp BTNC 12.5 và BTNC 19.
  • Kỹ thuật phân tích độ nhạy và mô phỏng: Kết hợp phân tích độ nhạy cục bộ (One-At-a-Time - OAT) và phân tích độ nhạy tổng thể (Global Sensitivity Analysis - GSA) thông qua mô phỏng Monte Carlo, sử dụng phần mềm Oracle Crystal Ball tích hợp trong MS Excel 2016.
  • Thiết kế kết cấu mặt đường: Sử dụng phần mềm AASHTOWare Pavement ME Design (DARWin-ME) để phân tích ứng xử và dự báo hư hỏng kết cấu mặt đường theo phương pháp cơ học thực nghiệm.

Nguồn tư liệu và dữ liệu

Dữ liệu đầu vào gồm các kết quả thí nghiệm trực tiếp chỉ tiêu cơ lý của 4 mác bitum, các chỉ tiêu cơ lý và độ dính bám của 3 loại đá (bazan, granit, đá vôi), kết quả thí nghiệm đúc mẫu nén nén dọc trục đo biến dạng phục hồi; số liệu khí hậu thực tế khu vực miền Bắc; dữ liệu đếm xe quy đổi lưu lượng giao thông trung bình hàng ngày hàng năm (AADT); các chỉ tiêu mô đun đàn hồi nền đất subgrade và các lớp móng cấp phối đá dăm từ một dự án đường ô tô cụ thể tại Việt Nam.

Nội dung chính theo từng chương

CHƯƠNG 1: Tổng quan về ảnh hưởng tính chất của bitum đến mô đun động của bê tông nhựa

Chương này tập trung phân tích cơ sở lý luận về tính đàn nhớt của bitum và BTN, khẳng định bitum dù chỉ chiếm từ 4% đến 6% khối lượng hỗn hợp nhưng là pha kết dính quyết định đặc tính lưu biến và khả năng chống biến dạng của BTN. Tác giả phân tích cấu tạo và nguyên lý hoạt động của thiết bị DSR và DMA, trong đó chỉ rõ ưu nhược điểm của từng phương pháp: thiết bị DSR gặp bước nhảy kết quả tại vùng nhiệt độ chuyển tiếp $30^\circ\text{C}$ do thay đổi kích thước đĩa đo (từ đĩa 8 mm sang 25 mm); trong khi thiết bị DMA cho dải tần số rộng hơn, đường cong thực nghiệm trơn và không bị biến động tại nhiệt độ chuyển tiếp.

Chương 1 tổng hợp các yếu tố ảnh hưởng tới $|E^|$, khẳng định rằng tính chất của bitum ($|G^|, \delta_b, \eta$) giữ vai trò chi phối áp đảo so với các yếu tố như độ rỗng dư ($V_a$), tỷ lệ bột khoáng trên bitum có hiệu ($filler/P_{be}$), hình dạng hay thành phần hạt. Tác giả khảo sát chi tiết cấu trúc toán học của các mô hình dự báo $|E^*|$ kinh điển từ Shell, Asphalt Institute, mô hình Witczak ban đầu, mô hình Witczak cải tiến, mô hình Hirsch, và mô hình Al-Khateeb.

Tiêu chí Mô hình Shell Mô hình Asphalt Institute Mô hình Witczak ban đầu Mô hình Witczak cải tiến Mô hình Hirsch
Năm công bố 1977 1979 1999 (NCHRP 1-37A) 2006 (Bari & Witczak) 2003 (Christensen)
Cơ sở dữ liệu Toán đồ thực nghiệm Thực nghiệm phòng 2.750 điểm dữ liệu 7.400 điểm dữ liệu 206 điểm dữ liệu
Thông số bitum Mô đun độ cứng phục hồi Độ nhớt $\eta$ tại $70^\circ\text{F}$ Độ nhớt $\eta$ theo nhiệt độ Mô đun cắt $ G^*
Thông số hỗn hợp Thể tích $V_b, V_g$ $V_a, V_b, P_{200}$, tần số $f$ $P_{200}, P_4, P_{3/8}, P_{3/4}, V_a, V_{beff}, f$ $P_{200}, P_4, P_{3/8}, P_{3/4}, V_a, V_{beff}, f$ Thể tích $VMA, VFA$, hệ số tiếp xúc $P_c$
Hạn chế chính Giới hạn $S_{t,T} > 5\text{ MPa}$ Độ nhớt quy đổi kém chính xác Độ nhớt không đo trực tiếp biến dạng cắt Cần thí nghiệm DSR phức tạp Mô hình bán thực nghiệm đơn giản hóa

CHƯƠNG 2: Nghiên cứu thực nghiệm và đề xuất sử dụng mô hình 2S2P1D để xây dựng đường cong chủ của mô đun cắt động và góc pha của một số loại bitum ở Việt Nam

Chương 2 trình bày kết quả thí nghiệm xác định các chỉ tiêu vật lý và cơ học lưu biến của 4 loại bitum (60/70, 35/50, 40/50 và bitum polyme PMBIII). Tác giả so sánh các mô hình toán học và cơ học (Maxwell, Voigt-Kelvin, Huet, Huet-Sayegh) và đề xuất lựa chọn mô hình 2S2P1D để xây dựng đường cong chủ của $|G^*|$ và góc pha $\delta_b$ tại nhiệt độ tham chiếu $T_{ref} = 30^\circ\text{C}$ và $40^\circ\text{C}$.

Mô hình 2S2P1D cho thấy khả năng mô tả chính xác ứng xử lưu biến của cả bitum thông thường lẫn bitum polyme cải tiến trên toàn bộ dải tần số quy đổi. Kết quả tính toán các hằng số $C_1, C_2$ của phương trình WLF cùng các thông số của mô hình 2S2P1D được kiểm định qua tiêu chuẩn thống kê nghiêm ngặt (hệ số xác định $R^2$ và tỷ số sai số chuẩn trên độ lệch chuẩn $S_e/S_y$), chứng minh mức độ dự báo rất tốt. Tác giả tiến hành phân loại và đánh giá cấp chất lượng của các loại bitum theo tiêu chuẩn Superpave AASHTO M320 thông qua chỉ tiêu chống lún hằn vệt bánh xe ban đầu $G^/\sin\delta$, đồng thời trích xuất các thông số $|G^|$, $\delta_b$ và độ nhớt $\eta$ tương ứng phục vụ làm đầu vào cho các mô hình dự báo $|E^*|$ ở chương kế tiếp.

CHƯƠNG 3: Nghiên cứu mối quan hệ thực nghiệm giữa tính chất của bitum và mô đun động của bê tông nhựa chặt ở Việt Nam

Chương 3 là phần thực nghiệm và mô hình hóa trung tâm của luận án. Đầu tiên, tác giả thực hiện thiết kế hỗn hợp cho BTNC 12.5 và BTNC 19 sử dụng cốt liệu đá bazan, đá granit, đá vôi và bột khoáng đá vôi theo TCVN 8819:2011. Tiếp đó, tác giả áp dụng phương pháp thiết kế thực nghiệm Taguchi với mảng trực giao $L_8(2^7)$ để giảm thiểu số lượng mẫu thử nhưng vẫn khảo sát đầy đủ các yếu tố ảnh hưởng.

Tác giả sử dụng kỹ thuật mô phỏng Monte Carlo trên phần mềm Oracle Crystal Ball để phân tích độ nhạy tổng thể (GSA) và phân tích OAT đối với các biến đầu vào trong các mô hình dự báo $|E^|$ của Hoa Kỳ. Kết quả biểu đồ Tornado khẳng định thông số tính chất bitum ($|G^|$, $\delta_b$ hoặc $\eta$) có độ nhạy lớn nhất đối với giá trị $|E^*|$, vượt trội so với các chỉ tiêu thành phần hạt và độ rỗng thể tích.

Sau khi tiến hành thí nghiệm đo $|E^|$ trên thiết bị Cooper tại Trường Đại học Giao thông vận tải theo tiêu chuẩn AASHTO TP-62 qua các cấp nhiệt độ và tần số, luận án xây dựng đường cong chủ của $|E^|$ (lấy nhiệt độ tham chiếu $10^\circ\text{C}$). Khi áp dụng trực tiếp 3 mô hình dự báo nguyên bản của Hoa Kỳ (Witczak ban đầu, Witczak cải tiến, Hirsch) cho số liệu thực nghiệm Việt Nam, sai lệch giữa giá trị dự báo và thực tế là rất lớn. Tác giả tiến hành hồi quy phi tuyến để hiệu chỉnh lại toàn bộ hệ số cho cả 3 mô hình. Đánh giá thống kê mức độ phù hợp (Goodness-of-fit statistic) trước và sau khi hiệu chỉnh chứng minh sự vượt trội của mô hình sau địa phương hóa:

| Mô hình dự báo $|E^|$ | Trạng thái | Đánh giá thống kê ($R^2$ và $S_e/S_y$) | Khả năng dự báo | | :--- | :--- | :--- | :--- | | Witczak ban đầu | Trước hiệu chỉnh | $R^2$ thấp, $S_e/S_y$ lớn, độ lệch cao | Kém phù hợp với BTNC Việt Nam | | | Sau hiệu chỉnh | $R^2$ tăng đáng kể, sai số giảm rõ rệt | Chấp nhận được | | Witczak cải tiến | Trước hiệu chỉnh | Tồn tại độ lệch hệ thống so với số liệu đo | Chưa phản ánh đúng vật liệu địa phương | | | Sau hiệu chỉnh | Đạt $R^2$ cao nhất, $S_e/S_y$ nhỏ nhất | Độ chính xác cao nhất | | Hirsch | Trước hiệu chỉnh | Đánh giá thấp giá trị $|E^|$ ở tần số cao | Độ tin cậy trung bình | | | Sau hiệu chỉnh | Mức độ phù hợp thống kê cải thiện tốt | Khá tốt, sau Witczak cải tiến |

Kết quả so sánh khẳng định mô hình Witczak cải tiến sau khi được hiệu chỉnh hệ số theo điều kiện Việt Nam cho kết quả dự báo $|E^*|$ chính xác nhất, phù hợp nhất để đưa vào phần mềm thiết kế ME.

CHƯƠNG 4: Nghiên cứu ảnh hưởng của mô đun động của bê tông nhựa chặt tới đặc trưng khai thác của kết cấu mặt đường mềm ở Việt Nam

Chương 4 đưa các kết quả nghiên cứu vật liệu vào ứng dụng tính toán thực tế trên phần mềm cơ học thực nghiệm AASHTOWare Pavement ME Design cho một dự án đường bộ tại Việt Nam. Tác giả thu thập đầy đủ các tập dữ liệu: giao thông phân bố theo phổ tải trọng trục và AADT; dữ liệu thời tiết khí hậu miền Bắc; chỉ tiêu mô đun của lớp nền thượng (subgrade layer), lớp móng cấp phối đá dăm (CPĐD), và các thông số vật liệu BTNC với giá trị $|E^*|$ đã xác lập ở Chương 3.

Nghiên cứu thiết lập 6 phương án kết cấu mặt đường thử nghiệm (kết cấu 1 đến kết cấu 6) và phân tích ứng xử theo các dạng hư hỏng quy định bởi ME:

  • Hư hỏng hằn lún vệt bánh xe (tổng biến dạng vĩnh cửu của các lớp kết cấu và nền đất).
  • Hư hỏng nứt mỏi (nứt phân bố dạng đáy - đỉnh / bottom-up cracking).
  • Hư hỏng nứt nhiệt theo thời gian.
  • Chỉ số độ gồ ghề quốc tế (International Roughness Index - IRI).

Các kết quả dự báo cho thấy phương án kết cấu ban đầu có thể không thỏa mãn chỉ tiêu giới hạn hằn lún vệt bánh xe hoặc nứt mỏi trong thời gian khai thác thiết kế. Tác giả tiến hành điều chỉnh chiều dày các lớp bê tông nhựa và lớp móng để đạt được kết cấu khả thi, đồng thời phân tích độ nhạy của các biến số khai thác theo lưu lượng AADT. Phân tích khẳng định việc sử dụng BTNC có mác bitum độ quánh cao hoặc bitum polymer PMBIII giúp tăng đáng kể giá trị $|E^*|$ ở nhiệt độ cao, từ đó kiểm soát hiệu quả biến dạng vệt bánh xe và duy trì chỉ số IRI trong ngưỡng cho phép.

Kết quả và những đóng góp mới

Đóng góp về mặt khoa học

  • Luận án đã luận giải đầy đủ cơ sở khoa học của các thông số đàn nhớt ($|G^|, \delta_b, |E^|$) sử dụng trong phương pháp thiết kế mặt đường cơ học thực nghiệm (ME), chứng minh tính khả thi của việc ứng dụng phương pháp ME vào điều kiện thực tế ở Việt Nam.
  • Thiết lập thành công mối quan hệ tương quan thực nghiệm giữa tính chất lưu biến của bitum và mô đun động $|E^*|$ của các loại bê tông nhựa chặt điển hình sản xuất từ vật liệu địa phương (đá bazan, đá granit, đá vôi).
  • Chứng minh bằng các công cụ phân tích độ nhạy OAT, GSA và mô phỏng Monte Carlo rằng tính chất của bitum ($|G^|, \delta_b, \eta$) là nhóm thông số ảnh hưởng lớn nhất đến giá trị mô đun động $|E^|$ của BTNC; đồng thời giá trị $|E^*|$ có vai trò chi phối hàng đầu đối với sự hình thành hư hỏng hằn lún vệt bánh xe và chỉ số IRI của kết cấu mặt đường mềm.
  • Đề xuất và ứng dụng thành công mô hình cơ học 2S2P1D kết hợp phương trình WLF để xây dựng đường cong chủ của bitum, cùng với việc ứng dụng mảng trực giao Taguchi $L_8(2^7)$ trong quy hoạch thực nghiệm vật liệu đường bộ tại Việt Nam.

Đóng góp về mặt thực tiễn

  • Xác định và công bố các chỉ tiêu kỹ thuật lưu biến theo cấp đặc tính khai thác (PG) theo tiêu chuẩn AASHTO M320 cho các loại bitum đang lưu hành phổ biến tại Việt Nam (60/70, 35/50, 40/50, PMBIII).
  • Xây dựng hệ thống cơ sở dữ liệu thực nghiệm về mô đun động $|E^*|$ cho hai loại hỗn hợp chủ lực BTNC 12.5 và BTNC 19 kết hợp với các nguồn cốt liệu điển hình tại các cấp nhiệt độ và tần số tải trọng khác nhau.
  • Đề xuất bộ hệ số hiệu chỉnh hoàn chỉnh cho 3 mô hình dự báo $|E^*|$ (Witczak ban đầu, Witczak cải tiến và Hirsch) phù hợp với nguồn vật liệu Việt Nam; trong đó khẳng định mô hình Witczak cải tiến sau hiệu chỉnh có độ chính xác cao nhất và khuyến nghị tích hợp vào quy trình tính toán thiết kế.
  • Đề xuất quy trình tính toán, lựa chọn và điều chỉnh chiều dày kết cấu mặt đường mềm theo phương pháp ME, phục vụ việc kiểm soát hư hỏng hằn lún vệt bánh xe tại các dự án đường ô tô chịu tải trọng nặng ở Việt Nam.

Hạn chế và hướng nghiên cứu tiếp

Văn bản luận án nêu rõ các giới hạn phạm vi nghiên cứu và đề xuất định hướng tiếp theo:

  • Giới hạn nghiên cứu: Đề tài mới chỉ tiến hành nghiên cứu trên 4 loại bitum (60/70, 35/50, 40/50 và PMBIII); 2 loại cấp phối bê tông nhựa chặt (BTNC 12.5 và BTNC 19); 3 loại cốt liệu đá (bazan, granit, đá vôi) và bột khoáng từ đá vôi; điều kiện nhiệt độ và dữ liệu thời tiết chỉ mới lấy đại diện cho khu vực miền Bắc Việt Nam.
  • Hướng nghiên cứu tiếp theo:
    • Mở rộng phạm vi nghiên cứu thực nghiệm xác định $|E^*|$ cho các loại hỗn hợp bê tông nhựa khác như bê tông nhựa rỗng (BTNR), bê tông nhựa có độ nhám cao (VTO), bê tông nhựa đá - mastix (SMA), và bê tông nhựa sử dụng bitum cải tiến khác.
    • Mở rộng việc thu thập và phân tích dữ liệu khí hậu cho các vùng miền khác tại Việt Nam (miền Trung và miền Nam) để hoàn thiện cơ sở dữ liệu quốc gia cho phương pháp ME.
    • Tiếp tục theo dõi, quan trắc hiện trường các đoạn đường thử nghiệm thực tế để kiểm chứng và hiệu chỉnh thêm các phương trình dự báo hư hỏng (nứt mỏi, lún vệt bánh, nứt nhiệt) theo thời gian khai thác.

Giá trị tham khảo

Luận án là tài liệu tham khảo có giá trị chuyên môn cao cho nhiều đối tượng trong ngành xây dựng công trình giao thông:

  • Đối với nghiên cứu sinh, học viên cao học và giảng viên chuyên ngành Đường ô tô và Công trình giao thông: Luận án cung cấp hệ thống lý thuyết chuẩn mực về cơ học đàn nhớt vật liệu đường, quy trình sử dụng thiết bị DSR, DMA và Cooper, phương pháp xây dựng đường cong chủ bằng mô hình 2S2P1D, kỹ thuật thiết kế thực nghiệm Taguchi $L_8(2^7)$ và phân tích độ nhạy tổng thể bằng mô phỏng ngẫu nhiên Monte Carlo.
  • Đối với các kỹ sư tư vấn thiết kế và cán bộ kỹ thuật: Có thể sử dụng trực tiếp các phương trình dự báo $|E^*|$ (đặc biệt là phương trình Witczak cải tiến đã hiệu chỉnh) để xác định thông số đầu vào cho bê tông nhựa chặt trong điều kiện thiếu thiết bị thí nghiệm động lực học phức tạp, giúp việc tính toán kết cấu mặt đường mềm theo xu hướng ME đạt độ chính xác cao.
  • Đối với các nhà quản lý kỹ thuật, cơ quan quản lý chuyên ngành (Bộ Giao thông vận tải, Cục Đường bộ Việt Nam): Cung cấp luận cứ khoa học và bộ thông số thực nghiệm để xây dựng, sửa đổi và ban hành các tiêu chuẩn thiết kế mặt đường mới thay thế các quy trình truyền thống vốn bộc lộ nhiều hạn chế về tuổi thọ công trình.

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

1. Trong ba mô hình dự báo mô đun động $|E^*|$ của Hoa Kỳ được luận án hiệu chỉnh, mô hình nào cho kết quả phù hợp nhất với điều kiện Việt Nam?
Dựa trên các tiêu chí đánh giá mức độ phù hợp thống kê (hệ số xác định $R^2$ cao nhất và tỷ số sai số chuẩn trên độ lệch chuẩn $S_e/S_y$ nhỏ nhất), mô hình Witczak cải tiến (Modified Witczak) sau khi hiệu chỉnh lại các hệ số thực nghiệm cho khả năng dự báo $|E^*|$ với độ chính xác cao nhất đối với bê tông nhựa chặt ở Việt Nam.

2. Thiết bị nào đã được luận án sử dụng để xác định mô đun động $|E^*|$ của hỗn hợp bê tông nhựa chặt trong phòng thí nghiệm?
Thí nghiệm xác định $|E^*|$ được thực hiện trên thiết bị Cooper (phát triển bởi Keith Cooper và Giáo sư Steven Brown thuộc Đại học Nottingham) đặt tại Phòng thí nghiệm vật liệu của Trường Đại học Giao thông vận tải, vận hành theo tiêu chuẩn AASHTO TP-62 với các tần số đo 0.5 Hz, 1 Hz, 5 Hz, 10 Hz và 25 Hz.

3. Tại sao luận án lựa chọn mô hình 2S2P1D thay vì các mô hình cơ học cổ điển để xây dựng đường cong chủ của bitum?
Mô hình 2S2P1D (gồm 2 lò xo, 2 phần tử biến dạng parabol và 1 giảm chấn) có khả năng mô tả chính xác cả mô đun cắt động $|G^*|$ và góc pha $\delta_b$ trên toàn bộ dải tần số rộng mà không bị sai lệch ở vùng nhiệt độ chuyển tiếp, thích hợp cho cả bitum thông thường lẫn bitum polyme cải tiến (PMBIII).

4. Kỹ thuật nào được luận án áp dụng để phân tích mức độ ảnh hưởng của các thông số đầu vào tới mô đun động $|E^*|$?
Luận án kết hợp phân tích độ nhạy cục bộ (OAT) và phân tích độ nhạy tổng thể (GSA) thông qua phương pháp mô phỏng ngẫu nhiên Monte Carlo, sử dụng phần mềm Oracle Crystal Ball tích hợp trong Microsoft Excel 2016.

5. Kết quả mô phỏng trên phần mềm AASHTOWare Pavement ME Design cho thấy thông số $|E^*|$ ảnh hưởng mạnh nhất tới dạng hư hỏng nào của mặt đường mềm?
Kết quả phân tích ứng xử và độ nhạy kết cấu chỉ ra rằng giá trị $|E^*|$ có ảnh hưởng nhiều nhất đến mức độ phát triển hư hỏng hằn lún vệt bánh xe (tổng biến dạng vĩnh cửu) và chỉ số độ gồ ghề quốc tế (IRI) của mặt đường mềm.

Kết luận

Luận án đã giải quyết hoàn chỉnh vấn đề thực nghiệm và mô hình hóa các chỉ tiêu đàn nhớt của bitum và hỗn hợp bê tông nhựa chặt tại Việt Nam, phục vụ trực tiếp cho việc địa phương hóa phương pháp thiết kế mặt đường cơ học thực nghiệm (ME). Công trình đã xác lập bộ dữ liệu thực nghiệm lưu biến của 4 loại bitum và mô đun động $|E^*|$ của BTNC 12.5 và BTNC 19, đồng thời chứng minh mô hình Witczak cải tiến sau hiệu chỉnh đạt độ chính xác dự báo cao nhất. Những kết quả nghiên cứu và quy trình phân tích kết cấu đề xuất trong luận án tạo cơ sở khoa học và thực tiễn vững chắc cho việc nâng cao chất lượng và tuổi thọ công trình đường bộ tại Việt Nam.