Tổng quan về luận án
Nghiên cứu của tác giả Phạm Duy Linh (2021) tại Trường Đại học Giao thông Vận tải với đề tài "Nghiên cứu đánh giá một số yếu tố ảnh hưởng đến chất lượng khai thác của mặt đường BTXM đường ô tô và sân bay trong điều kiện Việt Nam" (Chuyên ngành: Xây dựng đường ô tô và đường thành phố, Mã số: 9.05) đặt nền móng tiên phong trong việc giải quyết các xung đột cơ học nhiệt - tải trọng trên kết cấu áo đường cứng nhiều lớp tại Việt Nam. Trong bối cảnh hạ tầng giao thông đường bộ cấp cao và mạng lưới sân bay chịu áp lực nặng nề từ lưu lượng trục xe nặng (ESAL) kết hợp khí hậu nhiệt đới gió mùa nắng nóng khắc nghiệt, mặt đường bê tông xi măng (BTXM) thường xuyên xuất hiện các hư hỏng nghiêm trọng như nứt vỡ tấm, nứt mỏi cạnh khe co và phá hoại chèn kích tại các khe dãn (điển hình tại đường BRT Lê Văn Lương - Hà Nội, sân bay Tân Sơn Nhất và sân bay Đà Nẵng).
Thực trạng thiết kế kết cấu áo đường cứng tại Việt Nam tồn tại khoảng trống nghiên cứu (research gap) cốt lõi: các quy chuẩn hiện hành (QĐ 3230/2012/QĐ-BGTVT và TCVN 10907:2015) chủ yếu tiếp thu nguyên bản từ tiêu chuẩn Trung Quốc (JTG D40-2011) và tiêu chuẩn Liên Xô/Nga (SNiP 2.05.08-85, SNiP 2.05.02-85). Các tiêu chuẩn này giả định gradient nhiệt độ đơn vị ($T_g$) là hằng số tuyến tính ($0{,}86 - 0{,}92^\circ\text{C/cm}$), bỏ qua sự phi tuyến của trường nhiệt theo chiều sâu đối với các tấm có chiều dày lớn ($h = 38 - 40\text{ cm}$). Đồng thời, quy chuẩn coi lớp cách ly (bê tông nhựa chặt - BTNC, đá cát nhựa - SAMI, vữa bitum cát) là lớp phi biến dạng thuần túy cấu tạo, triệt tiêu ma sát mà không lượng hóa hiện tượng suy giảm mô đun đàn hồi của lớp bitum ở nhiệt độ bề mặt $65 - 67^\circ\text{C}$, dẫn đến phân bố lại nội lực bất lợi và gia tăng ứng suất kéo uốn đáy tấm BTXM.
Luận án thiết lập hệ thống câu hỏi nghiên cứu và giả thuyết khoa học:
- Câu hỏi nghiên cứu 1 (RQ1): Chiều dày và sự suy giảm độ cứng đàn nhớt của lớp cách ly ở nhiệt độ cao ảnh hưởng như thế nào đến trạng thái ứng suất - biến dạng trong hệ mặt đường BTXM hai lớp?
- Câu hỏi nghiên cứu 2 (RQ2): Quy luật phân bố trường nhiệt độ và gradient nhiệt thực tế theo chiều sâu tấm BTXM ($h = 22 - 40\text{ cm}$) trong điều kiện bức xạ mặt trời miền Bắc Việt Nam có tuân theo giả định tuyến tính hay không?
- Câu hỏi nghiên cứu 3 (RQ3): Sự cần thiết và khoảng cách bố trí khe dãn dãy tấm BTXM cần được xác lập theo mô hình cơ học giải tích nào để thay thế các quy định cấu tạo thuần túy ($10 - 12\text{ tấm/khe}$ hoặc $150\text{ m}$)?
- Giả thuyết 1 (H1): Khi nhiệt độ bề mặt đạt ngưỡng cực hạn ($65^\circ\text{C}$), lớp cách ly BTNC có chiều dày lớn ($h_{cl} \ge 4\text{ cm}$) bị suy giảm mô đun đàn hồi nghiêm trọng, làm gia tăng độ võng cục bộ và tăng ứng suất kéo uốn đáy tấm BTXM lớp trên so với giả thiết đồng biến dạng.
- Giả thuyết 2 (H2): Gradient nhiệt đơn vị $T_g$ nghịch biến phi tuyến với chiều dày tấm BTXM; tấm càng dày thì $T_g$ càng giảm, khiến các tính toán theo QĐ 3230/2012 đánh giá sai lệch ứng suất uốn vồng ($\sigma_t$).
- Giả thuyết 3 (H3): Áp lực dãn nở nhiệt tích lũy trên dãy tấm vượt quá khả năng hấp thụ biến dạng của ván gỗ chèn khe và cường độ chịu cắt của bê tông mép tấm, đòi hỏi phải tối ưu hóa giải tích khoảng cách khe dãn ($L_{xd}$) theo điều kiện nhiệt độ thi công và nhiệt độ khai thác.
Khung lý thuyết của nghiên cứu tích hợp: Lý thuyết nhiệt đàn hồi (Thermoelasticity) của Bradbury; Cơ học tấm mỏng trên nền đàn hồi nhiều lớp của Westergaard; và Phương pháp sai phân hữu hạn (SPHH). Phạm vi thực nghiệm bao gồm chuỗi quan trắc nhiệt độ liên tục trong 7 ngày nóng nhất năm 2018 tại Hà Nội trên các tấm thực tế có chiều dày $22\text{ cm}, 26\text{ cm}, 30\text{ cm}, 36\text{ cm}$ và $40\text{ cm}$, đem lại giải pháp định lượng đột phá cho công tác thiết kế kết cấu mặt đường ô tô và sân bay cấp cao.
Literature Review và Positioning
Nghiên cứu kết cấu mặt đường BTXM trên thế giới và tại Việt Nam phát triển qua ba dòng học thuật chính:
┌─────────────────────────────────────────────────────────────┐
│ TIẾN TRÌNH PHÁT TRIỂN LÝ THUYẾT MẶT ĐƯỜNG BTXM NHIỀU LỚP │
└──────────────────────────────┬──────────────────────────────┘
│
┌───────────────────────────────────────────┼───────────────────────────────────────────┐
▼ ▼ ▼
┌─────────────────────────────────┐ ┌─────────────────────────────────┐ ┌─────────────────────────────────┐
│ TRƯỜNG PHÁI BÁN THỰC NGHIỆM │ │ TRƯỜNG PHÁI CƠ HỌC GIẢI TÍCH│ │ TRƯỜNG PHÁI CƠ HỌC THỰC NGHIỆM│
│ Westergaard (1926), Bradbury │ │ SNiP (Nga), TCVN 10907:2015, │ │ AASHTO 1993/1998, ME PDG, │
│ - Tấm trên nền Winkler │ │ Nguyễn Văn Liên, Nguyễn H. Minh│ │ FAA (Mỹ), JTG D40 (TQ) │
│ - Giả định gradient tuyến tính │ │ - Tấm hai lớp chịu uốn đồng thời│ │ - Mô phỏng tích lũy hư hỏng │
│ - Chưa xét biến dạng lớp đệm │ │ - Bỏ qua biến dạng nén lớp đệm │ │ - Lớp BTN cách ly dày ≤ 5cm │
└─────────────────────────────────┘ └─────────────────────────────────┘ └─────────────────────────────────┘
│
▼
┌─────────────────────────────────────────────────────────────┐
│ KHOẢNG TRỐNG KHOA HỌC TRONG ĐIỀU KIỆN VIỆT NAM │
│ 1. Nhiệt độ mặt tấm 65-67°C làm mềm lớp cách ly BTN/SAMI │
│ 2. Gradient nhiệt Tg phi tuyến với chiều dày tấm 22-40cm │
│ 3. Khe dãn theo cấu tạo gây nứt vỡ, trồi tấm (blow-up) │
└──────────────────────────────┬──────────────────────────────┘
│
▼
┌─────────────────────────────────────────────────────────────┐
│ ĐÓNG GÓP ĐỘT PHÁ CỦA LUẬN ÁN (PHẠM DUY LINH) │
│ - Thuật toán SPHH (MLCP) xét biến dạng đàn nhớt lớp cách ly│
│ - Mô hình hồi quy phi tuyến tính toán Tg = f(h) và Ttb │
│ - Công thức giải tích xác định chính xác khoảng cách khe dãn│
└─────────────────────────────────────────────────────────────┘
-
Trường phái giải tích và bán thực nghiệm cổ điển: Westergaard (1926) xây dựng nền tảng bài toán tấm đơn vô hạn trên nền đàn hồi Winkler chịu tải trọng trục tĩnh. Bradbury (1938) mở rộng tính toán ứng suất nhiệt uốn vồng dựa trên giả thiết trường nhiệt tuyến tính. Dòng nghiên cứu này tại Việt Nam được tiếp nối bởi GS.TSKH Nguyễn Văn Liên (phương pháp SPHH cho hệ nhiều lớp không ma sát) và PGS.TS Phạm Huy Khang (quy đổi mô đun hệ nhiều lớp dính chặt). Tuy nhiên, các mô hình cổ điển chỉ xét hệ một lớp quy đổi hoặc hệ nhiều lớp dính kết hoàn toàn, chưa phản ánh được sự trượt - lún cơ học của tầng đệm cách ly.
-
Trường phái thiết kế mặt đường nhiều lớp của Liên Xô/Nga: Quy trình thiết kế SNiP 2.05.08-85 và các hiệu chỉnh (2004) đưa ra phương pháp phân chia mô men theo độ cứng kháng uốn tương đối ($D_t, D_d$) với giả định chuyển vị võng hai lớp bằng nhau ($w_1 = w_2$). Lớp cách ly quy định bằng 1-2 lớp giấy dầu hoặc bitum cát mỏng $3 - 5\text{ mm}$. TCVN 10907:2015 kế thừa phương pháp này nhưng bộc lộ nhược điểm khi áp dụng lớp đệm có chiều dày lớn.
-
Trường phái cơ học thực nghiệm Bắc Mỹ và Trung Quốc: AASHTO (1993, 1998), Hướng dẫn thiết kế Cơ học - Thực nghiệm (ME PDG), FAA (AC 150/5320-6D) và Trung Quốc (JTG D40-2011) cho phép sử dụng lớp cách ly bằng BTN nóng dày $30 - 50\text{ mm}$ để bù phụ tạo phẳng. Tuy nhiên, các tiêu chuẩn này được hiệu chuẩn ở vùng khí hậu ôn đới hoặc cận nhiệt đới; khi đưa vào Việt Nam, hiện tượng chảy mềm nhựa bitum dưới bức xạ cực đại chưa được đánh giá định lượng.
Về vấn đề khe dãn, tồn tại tranh luận học thuật sâu sắc giữa trường phái Bắc Mỹ (ACPA, FHWA, FAA) và trường phái khí hậu nhiệt đới:
- Quan điểm 1 (Bắc Mỹ): Triệt tiêu hầu hết khe dãn trên chính tuyến, chỉ bố trí tại vị trí tiếp giáp cầu hoặc giao cắt, vì cho rằng khe dãn làm giảm khả năng truyền lực qua thanh dowel bar, gây lún lệch và tăng chi phí bảo trì. Khoảng cách khe dãn nếu có được đẩy lên tới $150\text{ m}$ (500 ft).
- Quan điểm 2 (Thực tế nhiệt đới/Đề tài KC10-06): Khảo sát các sân bay miền Nam xây dựng trước 1975 theo tiêu chuẩn Mỹ (kích thước tấm $7{,}62 \times 7{,}62\text{ m}$, khe dãn $150\text{ m}$) cho thấy hơn 80% tổng số tấm bị nứt vỡ nhiệt nghiêm trọng. Báo cáo KC10-06 (1994) và các nghiên cứu của Nguyễn Hồng Minh (2007), Nguyễn Duy Đồng (2007), Phạm Đăng Nguyên (2016) khẳng định khoảng cách khe dãn quá lớn trong điều kiện bức xạ nhiệt cao dẫn tới hiện tượng chèn kích đẩy trồi (blow-up) và nứt vỡ mép tấm.
Nghiên cứu của Phạm Duy Linh định vị chính xác tại giao điểm của các tranh luận trên: lấp đầy khoảng trống cơ học về sự làm việc đồng thời của hệ ba cấu phần (Tấm BTXM mặt trên - Lớp cách ly biến dạng - Tầng móng cứng BT nghèo/Gia cố xi măng), đồng thời bác bỏ giả thiết gradient nhiệt tuyến tính bất biến theo chiều dày tấm.
Đóng góp lý thuyết và khung phân tích
Đóng góp cho lý thuyết
Luận án đóng góp 3 bước tiến lý thuyết quan trọng cho ngành cơ học công trình giao thông:
-
Mở rộng lý thuyết tương tác hệ nhiều lớp chịu uốn có tầng đệm đàn nhớt: Luận án phá vỡ giả thiết đồng biến dạng võng ($w_t = w_d$) của mô hình kết cấu áo đường cứng cổ điển (Westergaard, SNiP Nga). Bằng việc mô hình hóa lớp cách ly như một môi trường biến dạng liên tục chịu nén và trượt, nghiên cứu chứng minh rằng khi lớp cách ly có chiều dày $h_{cl} \ge 3\text{ cm}$ và mô đun đàn hồi suy giảm do nhiệt ($E_{cl}$ giảm từ $1200\text{ MPa}$ xuống $200 - 300\text{ MPa}$ khi nhiệt độ tăng từ $30^\circ\text{C}$ lên $60^\circ\text{C}$), phân bố nội lực sẽ chuyển đổi đột ngột: lớp mặt BTXM phải chịu mô men uốn cục bộ lớn hơn $15 - 22%$ so với tính toán theo công thức chia sẻ độ cứng lý thuyết ($D_t / \sum D$).
-
Xác lập quy luật phân bố nhiệt độ phi tuyến và hàm số gradient nhiệt phụ thuộc chiều dày ($h$): Luận án chứng minh quy luật truyền nhiệt trong tấm bê tông tuân theo phương trình vi phân Fourier phi tuyến một chiều với điều kiện biên biến thiên theo chu kỳ bức xạ mặt trời. Gradient nhiệt đơn vị cực đại ($T_g$) không phải là hằng số khu vực như định chuẩn của QĐ 3230/2012 ($T_g = 0{,}86^\circ\text{C/cm}$ miền Bắc) mà là một hàm mũ nghịch biến theo chiều dày tấm:
$$\Delta t(h) = T_g(h) \cdot h$$
Trong đó $T_g$ giảm mạnh từ $0{,}86^\circ\text{C/cm}$ (ở tấm $h = 20\text{ cm}$) xuống còn $0{,}48^\circ\text{C/cm}$ (ở tấm $h = 40\text{ cm}$).
- Mô hình giải tích ứng suất nhiệt và độ giãn nở dãy tấm có xét ma sát đáy và độ cứng ván chèn khe: Luận án phát triển hoàn chỉnh mô hình cân bằng lực co dãn nhiệt trên dãy tấm nhiều khe nối, kết hợp đồng thời ba tham số: hệ số giãn nở nhiệt cốt liệu ($\alpha_c$), lực ma sát đáy tấm theo hàm dịch chuyển trượt phi tuyến, và hệ số kháng ép nén thực tế của ván gỗ chèn khe ngâm tẩm ($E_v, [\sigma_{nen}]$).
Khung phân tích độc đáo
Khung phân tích của luận án tích hợp đa tầng lý thuyết:
- Lý thuyết nhiệt đàn hồi (Bradbury): Đánh giá trạng thái ứng suất uốn vồng ban ngày và uốn võng ban đêm.
- Lý thuyết cơ học tấm mỏng Kirchhoff-Love trên nền đàn hồi Pasternak/Winkler: Thiết lập phương trình vi phân đạo hàm riêng bậc 4 mô tả độ võng của từng lớp tấm.
- Lý thuyết lưu biến vật liệu bitum: Mô tả sự biến đổi độ nhớt - đàn hồi của lớp cách ly gốc nhựa theo gradient nhiệt độ truyền từ bề mặt xuống.
┌────────────────────────────────────────────────────────────────────────────────────────┐
│ KHUNG PHÂN TÍCH ĐỘC ĐÁO CỦA LUẬN ÁN │
├────────────────────────────────────────────────────────────────────────────────────────┤
│ │
│ [ BỨC XẠ MẶT TRỜI & KHÍ HẬU ] ──────► [ TRƯỜNG NHIỆT PHI TUYẾN: T(z, t) ] │
│ │ │ │
│ ▼ ▼ │
│ [ SUY GIẢM MÔ ĐUN LỚP CÁCH LY ] [ GRADIENT NHIỆT ĐỘ: Tg = f(h) ] │
│ - Bitum-cát / SAMI / BTNC - Ứng suất uốn vồng σt,max │
│ - E_cl(T) suy giảm ở 60°C - Nhiệt độ trung bình Ttb │
│ │ │ │
│ └────────────────────┬─────────────────────┘ │
│ │ │
│ ▼ │
│ [ THUẬT TOÁN SAI PHÂN HỮU HẠN CHUYÊN DỤNG: MLCP ] │
│ - Giải hệ phương trình vi phân bậc 4 tương tác │
│ - Phân tích chuyển vị võng w1(x,y) ≠ w2(x,y) │
│ │ │
│ ▼ │
│ [ TỐI ƯU HÓA KẾT CẤU: CHIỀU DÀY LỚP CÁCH LY & KHOẢNG CÁCH KHE DÃN L_xd ] │
│ │
└────────────────────────────────────────────────────────────────────────────────────────┘
Điều kiện biên cơ học được xác định chặt chẽ: biên tự do tại các mép khe dãn, biên liên kết truyền lực đàn hồi (thanh dowel) tại các khe co ngang, và liên kết chống dịch chuyển ngang (thanh tie bar) tại khe dọc.
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 phương pháp nghiên cứu định lượng thực nghiệm - số trị (Empirical-Computational Methodology). Thiết kế nghiên cứu đa tầng bao gồm:
- Thực nghiệm hiện trường đo đạc nhiệt độ tự nhiên quy mô lớn.
- Xây dựng thuật toán và mã nguồn tính toán sai phân hữu hạn (Finite Difference Method - SPHH) chuyên dụng mang tên MLCP (Multi-Layer Concrete Pavement).
- Mô hình hóa thống kê hồi quy phi tuyến trên phần mềm Minitab.
Quy trình nghiên cứu rigorous
Quy trình thực nghiệm được chuẩn hóa nhằm đảm bảo độ tin cậy và tính tái lập cao:
- Gia công mẫu và lắp đặt sensor hiện trường: Chế tạo 5 khối tấm BTXM đại diện với các chiều dày thực tế: $h = 22\text{ cm}, 26\text{ cm}, 30\text{ cm}, 36\text{ cm}$ (kết cấu đường ô tô) và $h = 40\text{ cm}$ (kết cấu sân bay) đặt tại trạm đo khu vực Hà Nội.
- Hệ thống thiết bị đo: Bố trí mạng lưới đầu đo cảm biến nhiệt (thermocouple sensors) độ chính xác $\pm 0{,}1^\circ\text{C}$ theo phương thẳng đứng từ bề mặt đến đáy tấm với mật độ $2 - 5\text{ cm}$/sensor. Hệ thống dây dẫn truyền dữ liệu về hộp ghi tự động (data logger).
- Thu thập dữ liệu cực hạn: Dữ liệu được trích xuất liên tục 24/24h trong đợt nắng nóng đỉnh điểm tháng 6 - 7/2018. Chuỗi 7 ngày nóng nhất có nhiệt độ không khí cực đại $T_{kk} = 40 - 41^\circ\text{C}$ và nhiệt độ bề mặt tấm đạt đỉnh $T_{bm} = 65 - 67^\circ\text{C}$ được lựa chọn làm tập dữ liệu chuẩn để tính toán gradient nhiệt tính toán.
Data và phân tích
┌─────────────────────────────────────────────────────────────────────────────────────────┐
│ THÔNG SỐ MÔ PHỎNG VÀ THỰC NGHIỆM ĐO ĐẠC CỦA LUẬN ÁN │
├──────────────────────────────────────┬──────────────────────────────────────────────────┤
│ Chiều dày tấm BTXM khảo sát ($h$) │ 22 cm, 26 cm, 30 cm, 36 cm, 40 cm │
│ Kích thước tấm mô phỏng │ 4,5 m × 3,75 m (ô tô) ; 5,0 m × 5,0 m (sân bay) │
│ Cường độ chịu nén BTXM ($f_c$) │ ≥ 30 MPa (đường ô tô) ; ≥ 35 MPa (sân bay) │
│ Cường độ kéo uốn BTXM ($f_r$) │ 4,5 - 5,0 MPa │
│ Loại vật liệu lớp cách ly │ Bitum-cát (0,5-2,0cm), SAMI (1,5-3,0cm), BTNC │
│ Chiều dày lớp cách ly BTNC ($h_{cl}$)│ 3,0 cm, 4,0 cm, 5,0 cm, 6,0 cm │
│ Biến thiên nhiệt độ lớp cách ly │ 30°C, 40°C, 50°C, 60°C ($E_{cl}$: 1200 → 250 MPa)│
│ Nhiệt độ bề mặt cực đại ($T_{bm}$) │ 65,0°C - 67,0°C │
│ Chênh lệch nhiệt đáy tấm 40cm │ Thực tế đo: 17,2°C; Tính toán giải tích: 17,3°C │
└──────────────────────────────────────┴──────────────────────────────────────────────────┘
Chương trình tính toán MLCP lập trình trên nền tảng SPHH giải bài toán hệ hai tấm chịu uốn có lớp cách ly trung gian biến dạng:
- Lớp trên: Tấm BTXM ($E_1, h_1, \mu_1$).
- Lớp giữa: Lớp cách ly ($h_{cl}, E_{cl}, G_{cl}$).
- Lớp dưới: Móng cứng BT nghèo/Gia cố xi măng ($E_2, h_2, \mu_2$).
- Nền đất: Mô hình nền đàn hồi tham số biến dạng tương đương ($E_0, K$).
Dữ liệu trường nhiệt được xử lý hồi quy trên phần mềm Minitab, kiểm định phân sai số bình phương nhỏ nhất ($R^2 > 0{,}98$), xác lập các hàm thực nghiệm dự báo $T_g$ và nhiệt độ trung bình tấm ($T_{tb}$).
Phát hiện đột phá và implications
Những phát hiện then chốt
-
Bác bỏ tính bất biến của gradient nhiệt độ đơn vị ($T_g$): Kết quả đo đạc thực nghiệm và mô phỏng số chỉ ra rằng, dưới cùng điều kiện bức xạ mặt trời cực đại tại Hà Nội ($T_{bm} = 65^\circ\text{C}$), gradient nhiệt đơn vị $T_g$ giảm dần khi chiều dày tấm tăng lên:
- Tấm $h = 20\text{ cm} \rightarrow \Delta t = 17{,}3^\circ\text{C} \rightarrow T_g = 0{,}86^\circ\text{C/cm}$ (trùng khớp QĐ 3230).
- Tấm $h = 24\text{ cm} \rightarrow \Delta t = 19{,}0^\circ\text{C} \rightarrow T_g = 0{,}79^\circ\text{C/cm}$.
- Tấm $h = 30\text{ cm} \rightarrow \Delta t = 20{,}5^\circ\text{C} \rightarrow T_g = 0{,}68^\circ\text{C/cm}$.
- Tấm $h = 40\text{ cm} \rightarrow \Delta t = 17{,}2 - 17{,}3^\circ\text{C} \rightarrow T_g = 0{,}43 - 0{,}48^\circ\text{C/cm}$.
Hàm hồi quy phi tuyến đề xuất xác định $T_g$ theo chiều dày tấm $h$ (cm):
$$T_g = a \cdot h^b \quad \text{hoặc} \quad T_g = \frac{\Delta t(h)}{h}$$
Việc áp dụng cứng nhắc $T_g = 0{,}86^\circ\text{C/cm}$ cho tấm sân bay $40\text{ cm}$ theo QĐ 3230 sẽ tính sai lệch $\Delta t = 34{,}4^\circ\text{C}$ (gấp đôi thực tế $17{,}3^\circ\text{C}$), dẫn đến tính toán ứng suất nhiệt uốn vồng $\sigma_{t,\max}$ tăng ảo $80 - 100%$, gây lãng phí vật liệu nghiêm trọng.
-
Ảnh hưởng bất lợi của chiều dày lớp cách ly BTNC ở nhiệt độ cao: Khi sử dụng lớp cách ly bằng BTNC dày $h_{cl} = 4 - 6\text{ cm}$ theo quy định của JTG D40-2011/QĐ 3230/2012, trong điều kiện nắng nóng mùa hè (nhiệt độ lớp cách ly đạt $50 - 60^\circ\text{C}$), mô đun đàn hồi của bê tông nhựa bị suy giảm từ $1200\text{ MPa}$ xuống dưới $300\text{ MPa}$. Sự lún ép biến dạng của lớp cách ly làm giảm khả năng chịu lực liên hợp của lớp móng dưới, làm tăng ứng suất kéo uốn đáy tấm BTXM lớp trên thêm $12 - 18%$ và tăng độ võng $w_{\max}$ thêm $25 - 35%$. Ngược lại, lớp cách ly bitum-cát ($0{,}5 - 1{,}0\text{ cm}$) hoặc màng SAMI ($1{,}5 - 2{,}0\text{ cm}$) duy trì ứng suất kéo uốn đáy tấm ở mức ổn định nhất.
-
Cơ chế phá hoại chèn ép mép tấm tại khe dãn và công thức xác định khoảng cách tối ưu: Dưới tác động của biên độ biến thiên nhiệt độ năm ($\Delta T = T_{tb,\max} - T_{tc}$ lên tới $35 - 40^\circ\text{C}$), dãy tấm BTXM giãn dài gây ép nén cực hạn lên ván gỗ chèn khe. Khi khoảng cách khe dãn lấy theo quy định cấu tạo cũ ($10 - 12\text{ tấm} \approx 45 - 60\text{ m}$ hoặc $150\text{ m}$), ván gỗ chèn khe bị ép vượt quá giới hạn đàn hồi biến dạng ($[\sigma_{nen}] \ge 3\text{ MPa}$), sinh ra lực chèn kích dọc trục cực lớn, gây nứt vỡ cạnh khe dãn (spalling) và đẩy trồi (blow-up) tấm BTXM.
┌────────────────────────────────────────────────────────────────────────────────────────┐
│ SO SÁNH KẾT QUẢ TÍNH TOÁN GRADIENT NHIỆT VÀ ỨNG SUẤT KÉO UỐN │
├───────────────────┬────────────────────────────┬───────────────────────────────────────┤
│ Chiều dày tấm (h) │ Gradient Tg (QĐ 3230/2012) │ Gradient Tg Đề xuất (Phạm Duy Linh) │
├───────────────────┼────────────────────────────┼───────────────────────────────────────┤
│ 20 cm │ 0,86 °C/cm (Δt = 17,2°C) │ 0,86 °C/cm (Δt = 17,2°C) │
│ 24 cm │ 0,86 °C/cm (Δt = 20,6°C) │ 0,79 °C/cm (Δt = 19,0°C) │
│ 30 cm │ 0,86 °C/cm (Δt = 25,8°C) │ 0,68 °C/cm (Δt = 20,5°C) │
│ 40 cm (Sân bay) │ 0,86 °C/cm (Δt = 34,4°C) │ 0,43 °C/cm (Δt = 17,3°C) │
├───────────────────┴────────────────────────────┴───────────────────────────────────────┤
│ => Hệ quả: Với tấm sân bay 40cm, công thức cũ làm tăng ảo ứng suất uốn vồng σt lên 98%│
└────────────────────────────────────────────────────────────────────────────────────────┘
Implications đa chiều
- Về mặt học thuật: Cung cấp bộ dữ liệu thực nghiệm chuẩn và thuật toán SPHH (chương trình MLCP) cho các bài toán tiếp xúc đàn nhớt nhiều lớp.
- Về mặt kỹ thuật công nghệ: Khuyến nghị thay đổi cấu tạo lớp cách ly:
- Tuyệt đối không dùng lớp BTNC dày $\ge 4\text{ cm}$ làm lớp cách ly thuần túy nếu không tính toán liên hợp kết cấu.
- Ưu tiên sử dụng lớp màng hấp thụ ứng suất SAMI ($h_{cl} = 1{,}5 - 2{,}0\text{ cm}$) hoặc bitum-cát ($h_{cl} \le 1{,}0\text{ cm}$).
- Về tiêu chuẩn hóa chính sách: Đề xuất sửa đổi Điều 6.5 và Điều 8 trong QĐ 3230/2012/QĐ-BGTVT và TCVN 10907:2015, đưa biểu thức $T_g = f(h)$ và quy trình kiểm toán giải tích khe dãn $L_{xd}$ vào tiêu chuẩn thiết kế quốc gia.
Limitations và Future Research
Mặc dù đạt được những đột phá học thuật xuất sắc, luận án vẫn ghi nhận 4 giới hạn nghiên cứu:
- Phạm vi địa lý số liệu nhiệt độ: Thực nghiệm hiện trường tập trung chủ yếu tại khu vực Hà Nội (đại diện khí hậu miền Bắc có biên độ nhiệt mùa đông - mùa hè lớn). Các vùng khí hậu bức xạ cao quanh năm như miền Trung và Tây Nam Bộ cần chuỗi đo mở rộng.
- Mô hình ma sát tiếp xúc đáy tấm: Luận án sử dụng hệ số ma sát trung bình giữa tấm và lớp móng; chưa mô hình hóa chi tiết hiện tượng suy giảm dính bám theo chu kỳ mỏi cơ học dưới tác động thủy lực (hiện tượng xói ngầm đáy tấm - pumping).
- Mô phỏng 2D/Tấm mỏng mở rộng: Thuật toán MLCP dựa trên lý thuyết tấm mỏng tương đương; chưa phát triển thành mô hình phần tử hữu hạn phi tuyến 3D đầy đủ để khảo sát tập trung ứng suất tại góc nhọn của khe nối có bố trí thanh dowel lệch trục.
- Vật liệu chèn khe mới: Nghiên cứu tập trung vào ván gỗ ngâm tẩm truyền thống; chưa khảo sát các vật liệu chèn khe đàn hồi cao phân tử thế hệ mới (bọt cao su xốp biến tính polyme).
Định hướng nghiên cứu tương lai:
- Thiết lập bản đồ phân vùng gradient nhiệt độ $T_g(h)$ và biên độ nhiệt tính toán $\Delta T$ cho toàn bộ 63 tỉnh thành Việt Nam.
- Ứng dụng mô hình cơ học rạn nứt phá hủy (Fracture Mechanics) 3D để mô phỏng sự lan truyền vết nứt mỏi cạnh khe co dưới tác động đồng thời của tải trọng bánh xe máy bay siêu nặng (Boeing 777, Airbus A350) và gradient nhiệt độ.
Tác động và ảnh hưởng
- Tác động học thuật: Luận án mở ra hướng tiếp cận chính xác trong cơ học mặt đường nhiều lớp tại các nước nhiệt đới. Các bài báo công bố từ luận án trên Tạp chí Giao thông Vận tải và kỷ yếu hội nghị quốc tế tạo tiền đề cho nhiều đề tài sau đại học.
- Hiệu quả kinh tế - kỹ thuật: Việc hiệu chỉnh gradient nhiệt đơn vị $T_g$ cho mặt đường sân bay ($h = 40\text{ cm}$) giúp tối ưu hóa chiều dày thiết kế, tiết kiệm từ $5 - 8%$ khối lượng bê tông xi măng cường độ cao mà vẫn đảm bảo hệ số độ tin cậy $\gamma_r \ge 1{,}25$.
- Kéo dài tuổi thọ công trình: Áp dụng cấu tạo lớp cách ly SAMI/bitum-cát tối ưu và bố trí chính xác khoảng cách khe dãn giải tích giúp triệt tiêu hiện tượng nứt vỡ mép tấm và phá hoại chèn kích, nâng tuổi thọ khai thác thực tế của mặt đường BTXM từ $15 - 20\text{ năm}$ lên $30\text{ năm}$.
Đối tượng hưởng lợi
- Nghiên cứu sinh & Giới học thuật: Kế thừa phương pháp sai phân hữu hạn giải bài toán hệ nhiều lớp biến dạng và phương pháp luận đo đạc nhiệt độ tấm bê tông khối lớn.
- Kỹ sư tư vấn thiết kế giao thông & sân bay: Sở hữu công thức tính nhanh gradient nhiệt $T_g(h)$, ứng suất uốn vồng $\sigma_t$, và công thức giải tích xác định chính xác khoảng cách bố trí khe dãn dãy tấm.
- Chủ đầu tư & Cơ quan quản lý nhà nước (Bộ GTVT, ACV, Cục Hàng không): Có cơ sở khoa học thực chứng để thẩm định thiết kế, kiểm soát chất lượng thi công và ban hành định mức, tiêu chuẩn kỹ thuật mới.
Câu hỏi chuyên sâu
-
Đóng góp lý thuyết độc đáo nhất của luận án là gì? Trả lời: Luận án đã mở rộng lý thuyết tấm trên nền đàn hồi thông qua việc xây dựng thuật toán sai phân hữu hạn (SPHH) cho phép tính toán hệ hai tấm chịu uốn có xét đến biến dạng nén và trượt của lớp cách ly trung gian bị suy giảm mô đun ở nhiệt độ cao, phá vỡ giả thiết đồng biến dạng võng ($w_1 = w_2$) tồn tại lâu năm trong các tiêu chuẩn SNiP (Nga) và TCVN 10907:2015.
-
Đột phá phương pháp luận so với các nghiên cứu quốc tế trước đây? Trả lời: So với tiêu chuẩn JTG D40-2011 của Trung Quốc và AASHTO của Mỹ vốn coi gradient nhiệt độ $T_g$ là một hằng số bất biến theo chiều dày tấm, luận án đã phối hợp đo đạc thực nghiệm 5 chiều dày tấm ($22 - 40\text{ cm}$) với phân tích hồi quy phi tuyến trên Minitab để chứng minh $T_g$ nghịch biến với chiều dày tấm ($T_g$ giảm từ $0{,}86^\circ\text{C/cm}$ ở tấm $20\text{ cm}$ xuống $0{,}43^\circ\text{C/cm}$ ở tấm $40\text{ cm}$).
-
Phát hiện bất ngờ và phản trực giác nhất từ dữ liệu thực nghiệm? Trả lời: Việc tăng chiều dày lớp cách ly bằng bê tông nhựa chặt ($h_{cl} \ge 4\text{ cm}$) theo khuyến nghị tạo phẳng của nước ngoài không làm tăng khả năng chịu lực của kết cấu mà ngược lại làm tăng ứng suất kéo uốn nguy hiểm tại đáy tấm BTXM lớp trên thêm $12 - 18%$ vào mùa hè, do lớp nhựa bị nung nóng đến $50 - 60^\circ\text{C}$ mất độ cứng, biến thành "đệm mềm" gây lún cục bộ.
-
Nghiên cứu có cung cấp giao thức tái lập (Replication Protocol) không? Trả lời: Có. Luận án cung cấp sơ đồ lưới sai phân chi tiết, hệ phương trình vi phân tương tác tại các nút lưới, thông số vật liệu đầu vào ($E, \mu, h$), sơ đồ bố trí cảm biến nhiệt độ và toàn bộ thuật toán khối của chương trình MLCP.
-
Chương trình nghiên cứu 10 năm tiếp theo được định hình như thế nào? Trả lời: Mở rộng xây dựng hệ thống cơ sở dữ liệu Big Data quan trắc nhiệt ẩm đa điểm trên toàn quốc, tích hợp mô hình biến đổi khí hậu cực đoan và ứng dụng trí tuệ nhân tạo (Machine Learning) để tự động hóa tính toán ứng suất mỏi nhiệt - tải trọng động thời gian thực cho mạng lưới đường cao tốc và sân bay thông minh.
Kết luận
Luận án tiến sĩ của Phạm Duy Linh (2021) là một công trình khoa học mẫu mực với 5 kết luận cốt lõi:
- Xác lập mô hình tính toán hệ nhiều lớp hoàn chỉnh: Xây dựng thành công phương pháp và chương trình tính toán SPHH (MLCP) đánh giá chính xác ảnh hưởng của độ nén ép và suy giảm mô đun đàn nhớt của lớp cách ly đến sự phân bố lại nội lực trong mặt đường BTXM.
- Khuyến nghị chuẩn mực cấu tạo lớp cách ly: Đề xuất chiều dày hợp lý cho lớp cách ly trong điều kiện khí hậu Việt Nam: ưu tiên lớp SAMI ($h_{cl} = 1{,}5 - 2{,}0\text{ cm}$) hoặc cát trộn bitum ($h_{cl} = 0{,}5 - 1{,}0\text{ cm}$); không áp dụng lớp BTNC dày $\ge 4\text{ cm}$ nếu coi là lớp cấu tạo phi biến dạng.
- Quy luật gradient nhiệt độ phi tuyến: Chứng minh gradient nhiệt độ đơn vị $T_g$ thay đổi theo chiều dày tấm, bác bỏ giả thiết $T_g = \text{const}$ của quy chuẩn cũ; đề xuất công thức tính $T_g(h)$ giúp tính toán chính xác ứng suất nhiệt uốn vồng cho cả đường ô tô và sân bay.
- Giải pháp giải tích cho khe dãn: Thiết lập công thức giải tích tính toán kiểm tra sự cần thiết và khoảng cách bố trí khe dãn $L_{xd}$, thay thế hoàn toàn việc bố trí theo cấu tạo cảm tính ($10 - 12\text{ tấm}$), ngăn chặn triệt để hư hỏng nứt vỡ mép tấm và chèn kích đẩy trồi.
- Cơ sở sửa đổi tiêu chuẩn ngành: Cung cấp luận cứ khoa học thực nghiệm vững chắc để sửa đổi, hoàn thiện bộ tiêu chuẩn quốc gia QĐ 3230/2012/QĐ-BGTVT và TCVN 10907:2015, đóng góp giá trị bền vững cho hạ tầng giao thông Việt Nam.