Tổng quan về luận án

Mặt đường bê tông xi măng (BTXM) tại các cảng hàng không là công trình hạ tầng chịu tác động đồng thời của tải trọng bánh máy bay siêu nặng (tải trọng khai thác vượt trội so với đường ô tô thông thường) và luồng khí phụt nhiệt độ cao từ động cơ phản lực, kết hợp với điều kiện môi trường khí hậu khắc nghiệt. Nhằm duy trì an toàn khai thác bay theo quy định của Tổ chức Hàng không Dân dụng Quốc tế (ICAO), phương pháp thí nghiệm không phá hủy (NDT) bằng thiết bị đo võng xung lực quả nặng rơi hạng nặng (Heavy Weight Deflectometer - HWD PRIMAX 2500, tải trọng xung động lên đến trên 150 kN / 34.000 lb) đã trở thành tiêu chuẩn kỹ thuật toàn cầu. Tuy nhiên, kết quả đo đạc độ võng thực tế và các thông số tính toán ngược (backcalculation) thường xuyên xảy ra sai số mang tính hệ thống do tác động của biến thiên trường nhiệt độ môi trường. Luận án tiến sĩ của nghiên cứu sinh Nguyễn Thị Ngân (2022) thuộc chuyên ngành Xây dựng đường ô tô và đường thành phố (Mã số: 9580205-2), Trường Đại học Xây dựng Hà Nội, dưới sự hướng dẫn khoa học của TS. Ngô Việt Đức, đã giải quyết căn bản bài toán khoa học mang tính cấp thiết này.

Khoảng trống nghiên cứu (research gap) cốt lõi được xác định rõ: Các quy trình kiểm định và tiêu chuẩn thiết kế hiện hành (TCVN 10907:2015, TCVN 8753:2011, FAA AC 150/5370-11B) mới dừng lại ở việc chuẩn hóa quy trình đo võng bằng HWD nhưng hoàn toàn bỏ qua việc định lượng và hiệu chỉnh tác động của trường nhiệt độ phi tuyến theo chiều sâu tấm đến sự làm việc của kết cấu mặt đường BTXM sân bay. Cụ thể, hiện tượng uốn vồng do chênh lệch nhiệt độ mặt trên - mặt dưới ($\Delta t = t_1 - t_2$) làm thay đổi trạng thái tiếp xúc đáy tấm với lớp móng, đồng thời sự giãn nở nhiệt làm biến thiên cơ chế chèn móc cốt liệu (aggregate interlock) tại khe nối, dẫn đến sự dao động giả tạo của các chỉ số kết cấu được tính toán ngược.

Nghiên cứu tập trung giải quyết 3 câu hỏi nghiên cứu ($RQ$) và kiểm chứng 3 giả thuyết khoa học ($H$):

  • $RQ_1$: Quy luật truyền nhiệt thực tế theo chiều sâu tấm BTXM sân bay dày 30–40 cm trong điều kiện khí hậu miền Bắc Việt Nam tuân theo mô hình giải tích nào? Giả thuyết $H_1$: Trường nhiệt độ trong tấm BTXM phân bố phi tuyến và có thể mô phỏng chính xác thông qua phương trình vi phân truyền nhiệt Fourier một chiều được hiệu chỉnh tham số thực nghiệm ($a, K$).
  • $RQ_2$: Sự biến thiên nhiệt độ bề mặt ($t_m$), nhiệt độ trung bình tấm ($t_{tb}$) và chênh lệch nhiệt độ ($\Delta t$) chi phối định lượng như thế nào đến độ võng ($d_i$), hệ số phản lực nền ($k$), mô đun đàn hồi tấm ($E_{bt}$), hiệu quả truyền lực khe nối ($LTE$) và chỉ số phân loại mặt đường ($PCN$)? Giả thuyết $H_2$: Độ võng tại tâm tấm phụ thuộc chủ yếu vào $\Delta t$, trong khi độ võng tại mép tấm và $LTE$ chịu sự tương tác đồng thời của cả $t_m$ và $\Delta t$.
  • $RQ_3$: Làm thế nào để chuẩn hóa số liệu đánh giá HWD tại hiện trường về trạng thái nhiệt độ tiêu chuẩn nhằm đảm bảo tính khách quan của chỉ số $PCN$? Giả thuyết $H_3$: Tồn tại một khung hàm hiệu chỉnh toán học và khung giờ đo tối ưu cho phép quy đổi các thông số đo động tại nhiệt độ thực địa về điều kiện chuẩn ($20^\circ\text{C}$).

Khung lý thuyết của luận án được xây dựng dựa trên sự tích hợp giữa: (1) Lý thuyết truyền nhiệt Fourier trong môi trường bán không gian vô hạn; (2) Lý thuyết uốn tấm đàn hồi trên nền đàn hồi Winkler của Westergaard; (3) Lý thuyết tiếp xúc và truyền lực qua khe nối (aggregate interlock và dowel action); và (4) Khung phân loại sức chịu tải kết cấu ACN-PCN của ICAO. Nghiên cứu thực nghiệm được triển khai với chuỗi quan trắc nhiệt độ liên tục 365 ngày trên các tấm mẫu thực nghiệm kết hợp thí nghiệm HWD trên đường cất hạ cánh (CHC), đường lăn tại Cảng hàng không Quốc tế Vân Đồn và Nội Bài.

Literature Review và Positioning

Lịch sử nghiên cứu cơ học mặt đường cứng ghi nhận ba dòng lý thuyết chủ lưu trong phân tích ứng suất nhiệt và đánh giá không phá hủy:

  1. Dòng lý thuyết phân bố nhiệt độ trong tấm bê tông: Bắt đầu từ giả định phân bố nhiệt tuyến tính cổ điển của Westergaard (1926), Bradbury (1938), phát triển qua các mô hình truyền nhiệt bán giải tích của Kaufman, Barber (1957), Choubane & Tia (1992), Qin & Hiller (2011), và các nghiên cứu của Dương Học Hải, Hoàng Tùng tại Việt Nam.
  2. Dòng lý thuyết tính toán ngược thông số mặt đường từ thiết bị FWD/HWD: Được định hình bởi Foxworthy (1984), Ioannides và cộng sự (1989) với phương pháp diện tích chậu võng (AREA Method), cùng hệ thống phần mềm COMFAA, BAKFAA của Cục Hàng không Liên bang Hoa Kỳ (FAA).
  3. Dòng lý thuyết truyền lực và ứng suất tiếp xúc khe nối: Được xây dựng bởi Teller & Sutherland (1935), Eisenmann (1979), Zollinger và cộng sự (1998) về cơ chế chèn móc cốt liệu và thanh truyền lực.

Văn bản khoa học tồn tại những tranh luận học thuật sâu sắc giữa hai quan điểm đối lập: Quan điểm truyền thống xem xét gradien nhiệt độ trong tấm là hàm tuyến tính đơn giản nhằm giảm thiểu độ phức tạp tính toán thiết kế; ngược lại, các công trình hiện đại chứng minh sự phân bố nhiệt độ thực tế mang tính phi tuyến cao, trong đó sự biến thiên mạnh nhất tập trung ở nửa trên chiều dày tấm tiếp xúc trực tiếp với bức xạ mặt trời.

                  ┌─────────────────────────────────────────────────────────┐
                  │              LITERATURE POSITIONING MAP                 │
                  └─────────────────────────────────────────────────────────┘
                                               │
             ┌─────────────────────────────────┴─────────────────────────────────┐
             ▼                                                                   ▼
┌─────────────────────────┐                                         ┌─────────────────────────┐
│ International Studies   │                                         │ Thesis Positioning     │
│ • USA (FHWA/FAA):       │                                         │ (Nguyen Thi Ngan, 2022) │
│   k drops 80->20 kPa/mm │                                         │ • Calibrated 1D Fourier │
│   Reference Ts = 20°C   │                                         │   with dynamic K, a     │
│ • Brazil (Balbo et al): │ ──[Identified Gap: Tropical Climate]──> │ • Quantitative links    │
│   Delta t = 15°C summer │                                         │   Delta t, Tm to HWD    │
│ • Czech (Stehlik):      │                                         │ • Multi-parameter       │
│   Deflection ~ Gradient │                                         │   calibration model     │
│ • China (Zhao et al):   │                                         │   (d, k, LTE, PCN,      │
│   Tri-phasic daily cycle│                                         │    x, Mu, sigma_s)      │
└─────────────────────────┘                                         └─────────────────────────┘

Vị thế học thuật của luận án được xác lập rõ khi so sánh trực tiếp với các nghiên cứu quốc tế tiêu biểu:

  • So với nghiên cứu tại Hoa Kỳ (FHWA/FAA, Foxworthy, 2000s): Các nghiên cứu tại Mỹ khuyến cáo nhiệt độ đo FWD không vượt quá $27^\circ\text{C}$ và chọn nhiệt độ chuẩn là $20^\circ\text{C}$, ghi nhận hệ số nền $k$ sụt giảm từ $80\text{ kPa/mm}$ xuống $20\text{ kPa/mm}$ khi nhiệt độ mặt đường tăng lên $32,2^\circ\text{C}$. Luận án đã mở rộng phạm vi kiểm chứng này trên kết cấu tấm sân bay dày $30\text{--}40\text{ cm}$ chịu tải trọng HWD siêu nặng $150\text{ kN}$.
  • So với nghiên cứu tại Brazil (Balbo & Severi, 2002; [65][77][84]): Nghiên cứu tại Brazil ghi nhận trong điều kiện nhiệt đới, chênh lệch nhiệt độ cực đại $\Delta t$ đạt $15^\circ\text{C}$ vào mùa hè và $10^\circ\text{C}$ vào mùa đông đối với tấm dày $25\text{ cm}$. Luận án đã phát hiện trong điều kiện miền Bắc Việt Nam, $\Delta t$ trên tấm $40\text{ cm}$ biến thiên phức tạp hơn do ảnh hưởng của độ ẩm gió mùa, với biên độ nhiệt dao động ngày đêm sâu hơn.
  • So với nghiên cứu tại Cộng hòa Séc (Stehlík et al. [62]) và Trung Quốc (Zhao et al. [61][88]): Stehlík chứng minh tương quan bậc nhất giữa gradien nhiệt và độ võng, trong khi Zhao chỉ ra quy luật 3 giai đoạn của $LTE$ theo nhiệt độ bề mặt (tăng - ổn định từ 13h-17h - giảm). Luận án của Nguyễn Thị Ngân đã nâng cao mức độ đóng góp bằng cách mô hình hóa đồng thời cả 2 nhóm thông số: nhóm làm việc ($d_i, k, LTE, PCN$) và nhóm thiết kế kết cấu ($x, M_u, \sigma_s$).

Đó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 quan trọng vào lý thuyết cơ học mặt đường thông qua việc hiệu chỉnh phương trình vi phân truyền nhiệt Fourier một chiều trong môi trường bán không gian vô hạn:

$$\frac{\partial t}{\partial T} = a \frac{\partial^2 t}{\partial z^2}$$

Trong đó $t$ là nhiệt độ ($^\circ\text{C}$), $T$ là thời gian (giờ), $z$ là độ sâu tính từ bề mặt tấm (cm), và $a$ là hệ số truyền nhiệt độ ($a = \lambda / (c_v \cdot \rho)$, đơn vị $\text{m}^2/\text{h}$). Dạng nghiệm tổng quát phục vụ tính toán phân bố nhiệt được hiệu chỉnh:

$$t(z, T) = t_{tb.mat} - K \cdot z + t_{n\max} \cdot e^{-z \sqrt{\frac{\omega}{2a}}} \cdot \sin\left(\omega T - z \sqrt{\frac{\omega}{2a}} - \frac{\pi}{2}\right)$$

Nghiên cứu đã xác lập 3 mệnh đề lý thuyết nền tảng:

  • Mệnh đề 1: Hệ số truyền nhiệt độ thực tế $a$ và thông số phân bố trường nhiệt độ dừng $K$ không phải là hằng số cố định mà biến thiên theo mùa và theo độ sâu lớp vật liệu do tương tác nhiệt ẩm nội tại.
  • Mệnh đề 2: Cơ chế uốn vồng nhiệt ($\Delta t > 0$ uốn vồng dương vào ban ngày, $\Delta t < 0$ uốn vồng âm vào ban đêm) làm thay đổi căn bản điều kiện biên tiếp xúc giữa tấm BTXM và lớp móng, dẫn đến sự suy giảm độ cứng chống uốn biểu kiến khi chịu xung tải HWD.
  • Mệnh đề 3: Chỉ số $LTE$ tại khe nối chịu sự chi phối phi tuyến của hiệu ứng chèn móc cốt liệu: nhiệt độ trung bình $t_{tb}$ kích hoạt sự giãn nở hình học làm tăng diện tích tiếp xúc khe, nâng cao khả năng truyền tải trọng cắt ngang qua mối nối.
┌────────────────────────────────────────────────────────────────────────────────────────┐
│                        CONCEPTUAL & THEORETICAL FRAMEWORK                              │
└────────────────────────────────────────────────────────────────────────────────────────┘
                                             │
      ┌──────────────────────────────────────┼──────────────────────────────────────┐
      ▼                                      ▼                                      ▼
┌───────────────────────────┐  ┌───────────────────────────┐  ┌───────────────────────────┐
│ 1. 1D Fourier Heat Flow   │  │ 2. Westergaard-Winkler    │  │ 3. Aggregate Interlock    │
│    Calibration            │  │    Curling Contact Model  │  │    Joint Transfer Model   │
│ • a = lambda / (c_v * rho)│  │ • Day: Delta t > 0        │  │ • Expansion -> tight joint│
│ • Dynamic parameter K     │  │   (Upward curling)        │  │   -> LTE increases        │
│ • Damping depth z_max     │  │ • Night: Delta t < 0      │  │ • Contraction -> open gap │
│                           │  │   (Downward curling)      │  │   -> LTE drops            │
└───────────────────────────┘  └───────────────────────────┘  └───────────────────────────┘
      │                                      │                                      │
      └──────────────────────────────────────┼──────────────────────────────────────┘
                                             ▼
                 ┌─────────────────────────────────────────────────────────┐
                 │ Coupled Analytical Model & Multi-parameter Calibration  │
                 │ Outputs: Deflection (d0), Subgrade (k), LTE, PCN,       │
                 │          Neutral axis (x), Moment (Mu), Stress (sigma_s)│
                 └─────────────────────────────────────────────────────────┘

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

Khung phân tích của luận án tích hợp thành công ba lý thuyết trụ cột: Lý thuyết nhiệt vật lý công trình, Lý thuyết bản đàn hồi tiếp xúc phi tuyến, và Lý thuyết động lực học kết cấu áo đường chịu tải trọng va đập. Tiếp cận phân tích mới này thiết lập chuỗi chuyển tiếp dữ liệu: từ số liệu quan trắc nhiệt độ bề mặt $t_m$ $\rightarrow$ suy diễn trường nhiệt theo chiều sâu $t(z, T)$ $\rightarrow$ xác định $\Delta t$ $\rightarrow$ hiệu chỉnh chậu võng động HWD $\rightarrow$ tính toán ngược các thông số đàn hồi ($k, E_{bt}$) và thông số khai thác ($LTE, PCN, M_u$). Điều kiện biên được xác định rõ ràng cho kết cấu mặt đường BTXM sân bay có chiều dày $30\text{--}40\text{ cm}$, móng cấp phối đá dăm gia cố xi măng hoặc bê tông nhựa, trong dải nhiệt độ bề mặt thực tế từ $10^\circ\text{C}$ đến $60^\circ\text{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) với phương pháp tiếp cận định lượng đa tầng (multi-level quantitative experimental design), kết hợp chặt chẽ giữa nghiên cứu lý thuyết giải tích, thực nghiệm hiện trường diện rộng và mô hình hóa xác suất thống kê. Thiết kế thí nghiệm đa tầng bao gồm:

  • Tầng 1 (Quan trắc nhiệt vật lý chuyên sâu): Thiết lập hệ thống cảm biến nhiệt chôn sẵn trong tấm BTXM thực nghiệm dày $40\text{ cm}$ để thu thập chuỗi thời gian liên tục 24/7 trong 12 tháng.
  • Tầng 2 (Thực nghiệm động lực học kết cấu HWD): Tiến hành thử tải động bằng thiết bị HWD trên đường CHC và đường lăn sân bay thực tế (Vân Đồn, Nội Bài) tại các thời điểm nhiệt độ khác nhau trong ngày và các mùa trong năm.
┌────────────────────────────────────────────────────────────────────────────────────────┐
│                        EXPERIMENTAL ARCHITECTURE & FLOW                                │
└────────────────────────────────────────────────────────────────────────────────────────┘
                                             │
             ┌───────────────────────────────┴───────────────────────────────┐
             ▼                                                               ▼
┌─────────────────────────────────────────┐     ┌─────────────────────────────────────────┐
│ Multi-depth Thermal Acquisition         │     │ Dynamic In-situ HWD Testing             │
│ • EXTECH SDL200 datalogger              │     │ • HWD PRIMAX 2500 (Sweco, 150 kN load)  │
│ • JMZX-300L contact thermal sensors     │     │ • 10 geophones (0 to 210 cm)            │
│ • Depths: z = 1.5 cm, 20 cm, 40 cm      │     │ • Center slab & Joint edge loading      │
│ • 365-day continuous recording          │     │ • Runway / Taxiway (Van Don & Noi Bai)  │
└─────────────────────────────────────────┘     └─────────────────────────────────────────┘
             │                                                               │
             └───────────────────────────────┬───────────────────────────────┘
                                             ▼
┌────────────────────────────────────────────────────────────────────────────────────────┐
│ Statistical Processing & Backcalculation Engine                                        │
│ • Minitab regression analysis (R^2, Se, Sy, p < 0.05)                                  │
│ • AREA Method for relative stiffness radius (l_est) and subgrade reaction (k_est)      │
│ • COMFAA software for standardized PCN determination (ICAO compliance)                 │
└────────────────────────────────────────────────────────────────────────────────────────┘

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu được kiểm soát nghiêm ngặt theo các tiêu chuẩn quốc tế:

  • Thiết bị đo nhiệt: Sử dụng thiết bị ghi nhiệt EXTECH SDL200 độ chính xác cao kết hợp các đầu đo nhiệt kiểu tiếp xúc JMZX-300L và thiết bị hồng ngoại cầm tay UNI-T. Các cảm biến được bố trí tại các độ sâu $z = 1,5\text{ cm}$ (sát bề mặt), $z = 20\text{ cm}$ (giữa tấm), và $z = 40\text{ cm}$ (đáy tấm).
  • Thiết bị đo tải trọng động: Thiết bị HWD PRIMAX 2500 do hãng Sweco chế tạo, được hiệu chuẩn kiểm định quốc tế. Thiết bị tạo xung tải trọng tác dụng từ $70\text{ kN}$ đến trên $150\text{ kN}$ thông qua bộ giảm chấn truyền lực qua tấm nén đường kính $D = 30\text{ cm}$ trong khoảng thời gian từ $20\text{--}30\text{ ms}$.
  • Hệ thống cảm biến đo võng: Bố trí 10 đầu đo vận tốc địa chấn (geophones) hướng tâm tại các khoảng cách cố định: $0, 20, 30, 45, 60, 90, 120, 150, 180, 210\text{ cm}$ từ tâm tấm ép.
  • Quy trình tam giác hóa dữ liệu (Triangulation): Kết hợp đối chứng chéo giữa số liệu đo đạc nhiệt độ thực tế, kết quả giải tích phương trình Fourier, và phân tích hồi quy thống kê đa biến.

Data và phân tích

Phân tích dữ liệu được thực hiện thông qua công cụ chuyên dụng Minitab để xây dựng phương trình tương quan hồi quy và phần mềm COMFAA của FAA để tính toán chỉ số PCN.

  • Phương pháp diện tích chậu võng (AREA method) được ứng dụng để tính toán ngược bán kính độ cứng tương đối $l_{est}$ và mô đun phản lực nền $k_{est}$:

$$AREA = 6 \cdot \left[1 + 2 \cdot \left(\frac{d_{30}}{d_0}\right) + 2 \cdot \left(\frac{d_{60}}{d_0}\right) + \left(\frac{d_{90}}{d_0}\right)\right]$$

  • Hiệu quả truyền lực tại khe nối ($LTE$) được xác định chính xác từ độ võng tấm chịu tải ($d_{\text{loaded}}$) và tấm không chịu tải ($d_{\text{unloaded}}$):

$$LTE = \frac{d_{\text{unloaded}}}{d_{\text{loaded}}} \cdot 100% \quad \text{hoặc} \quad LTE = \frac{2 \cdot d_{\text{unloaded}}}{d_{\text{loaded}} + d_{\text{unloaded}}} \cdot 100%$$

Mức độ phù hợp của các mô hình hồi quy được kiểm soát thông qua hệ số tương quan $R^2$, sai số chuẩn ước tính $S_e$, và độ lệch chuẩn $S_y$ với mức ý nghĩa thống kê $p < 0,05$.

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

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

  1. Quy luật phân bố nhiệt độ phi tuyến và chiều sâu tắt dao động nhiệt: Số liệu thực nghiệm 365 ngày khẳng định nhiệt độ phân bố phi tuyến theo chiều sâu tấm BTXM sân bay. Giá trị thông số $K$ trung bình biến thiên mạnh theo các tháng trong năm (tại $z = 1,5\text{ cm}$, $z = 20\text{ cm}$ và $z = 40\text{ cm}$). Chiều sâu tắt biên độ dao động nhiệt ngày đêm $z_{\max}$ trên tấm bê tông sân bay miền Bắc Việt Nam được xác định chính xác, tạo cơ sở thực nghiệm hiệu chỉnh điều kiện biên cho phương trình Fourier.
  2. Tác động phân hóa của nhiệt độ tới độ võng tại tâm và mép tấm: Độ võng tại tâm tấm ($d_0$) có mối tương quan thuận chặt chẽ với chênh lệch nhiệt độ $\Delta t$ ($R^2 > 0,85$). Khi $\Delta t > 0$ tăng cao vào buổi trưa, hiện tượng uốn vồng dương làm tâm tấm tách khỏi móng, dẫn đến độ võng đo được dưới tải trọng HWD tăng mạnh. Ngược lại, độ võng tại mép tấm phụ thuộc trực tiếp vào nhiệt độ bề mặt $t_m$ do sự co giãn khe nối.
  3. Sự biến động phi tuyến của hiệu quả truyền lực khe nối ($LTE$): Tại các khe nối không có thanh truyền lực hoặc thanh truyền lực giảm hiệu quả, $LTE$ dao động cực lớn từ dưới $20%$ vào sáng sớm (khi nhiệt độ thấp, tấm co lại, khe hở mở rộng làm mất hiệu ứng aggregate interlock) lên đến xấp xỉ $80\text{--}90%$ vào khoảng thời gian từ 13h đến 17h (khi nhiệt độ tấm đạt cực đại, tấm giãn dài chèn chặt khe nối).
  4. Sai lệch đáng kể của hệ số phản lực nền biểu kiến ($k$): Giá trị $k$ xác định bằng HWD vào thời điểm giữa trưa (nhiệt độ bề mặt cao, uốn vồng dương lớn) bị đánh giá thấp hơn từ $40\text{--}60%$ so với giá trị đo vào sáng sớm trên cùng một vị trí kết cấu, dù cường độ nền đất không hề thay đổi.
  5. Biến thiên sai lệch của chỉ số phân loại mặt đường ($PCN$): Chỉ số $PCN$ tính toán từ phần mềm COMFAA nếu không được hiệu chỉnh nhiệt độ sẽ dao động từ $15\text{--}25%$ trong ngày. Điều này có thể dẫn đến kết luận sai lầm về năng lực tiếp nhận dòng máy bay khai thác (như Boeing 777, 787, Airbus A350).
┌────────────────────────────────────────────────────────────────────────────────────────┐
│                 SUMMARY OF KEY QUANTITATIVE RELATIONSHIPS & FINDINGS                   │
├───────────────────────┬─────────────────────────────────┬──────────────────────────────┤
│ Parameter Evaluated   │ Primary Temperature Driver      │ Physical & Structural Effect │
├───────────────────────┼─────────────────────────────────┼──────────────────────────────┤
│ Center Deflection (d0)│ Temperature Difference (Delta t)│ Daytime upward curling lifts │
│                       │                                 │ slab center -> d0 increases  │
├───────────────────────┼─────────────────────────────────┼──────────────────────────────┤
│ Joint Edge Deflection │ Surface Temp (Tm) & (Delta t)   │ Joint expansion closes gaps; │
│                       │                                 │ subgrade contact variation   │
├───────────────────────┼─────────────────────────────────┼──────────────────────────────┤
│ Load Transfer (LTE)   │ Mean Slab Temp (T_tb) & Tm      │ LTE swings from <20% (dawn)  │
│                       │                                 │ to >80% (13h-17h, interlock) │
├───────────────────────┼─────────────────────────────────┼──────────────────────────────┤
│ Subgrade Reaction (k) │ Thermal Curling (Delta t)       │ Apparent k drops up to 40-60%│
│                       │                                 │ at peak noon temperatures    │
├───────────────────────┼─────────────────────────────────┼──────────────────────────────┤
│ Pavement Index (PCN)  │ Coupled Thermal Basin (d_i)     │ Uncorrected PCN fluctuates   │
│                       │                                 │ 15-25% across daily cycles   │
└───────────────────────┴─────────────────────────────────┴──────────────────────────────┘

Implications đa chiều

  • Về mặt lý thuyết: Mở rộng cơ học mặt đường sân bay từ phân tích tĩnh sang phân tích tương tác nhiệt động lực học đa trường (thermo-mechanical coupled analysis).
  • Về mặt phương pháp luận: Cung cấp thuật toán và bảng tra hệ số hiệu chỉnh đưa độ võng $d_i$, hệ số $k$, $LTE$, và $PCN$ đo ở nhiệt độ bất kỳ về điều kiện chuẩn ($T_0 = 20^\circ\text{C}$).
  • Về mặt ứng dụng thực tiễn: Kiến nghị khung thời gian tối ưu thực hiện thí nghiệm HWD tại các sân bay Việt Nam (khuyến cáo đo vào thời điểm chênh lệch nhiệt độ ổn định, sáng sớm hoặc chiều mát) nhằm giảm thiểu sai số uốn vồng.
  • Về mặt chính sách và quản lý bay: Giúp Cục Hàng không Việt Nam và Tổng công ty Cảng hàng không Việt Nam (ACV) đánh giá chính xác chỉ số PCN thực tế của đường cất hạ cánh, ngăn ngừa nguy cơ cấm tải nhầm hoặc khai thác quá tải gây phá hoại kết cấu mặt đường.

Limitations và Future Research

Nghiên cứu có các giới hạn khoa học cần được ghi nhận:

  1. Mô hình truyền nhiệt một chiều: Luận án tiếp cận bài toán truyền nhiệt một chiều theo phương đứng ($z$), chưa mô phỏng trường nhiệt độ 3 chiều tại các vị trí góc tấm chịu bức xạ biên đa hướng.
  2. Phạm vi khí hậu khu vực: Số liệu thực nghiệm tập trung tại khu vực miền Bắc Việt Nam (Cảng hàng không Quốc tế Vân Đồn, Nội Bài); cần mở rộng kiểm chứng tại khu vực miền Trung và miền Nam với đặc thù bức xạ nhiệt và mùa mưa khác biệt.
  3. Đặc tính ẩm kết hợp (hydro-thermal): Nghiên cứu tập trung vào gradient nhiệt độ mà chưa tích hợp đầy đủ biến thiên gradient độ ẩm theo chiều sâu tấm bê tông trong mùa mưa bão.

Chương trình nghiên cứu tiếp theo (Future Research Agenda):

  • Phát triển mô hình phần tử hữu hạn 3D mô phỏng tương tác động lực học giữa bánh máy bay - tấm BTXM - nền móng có xét đến gradient nhiệt ẩm phi tuyến.
  • Ứng dụng trí tuệ nhân tạo (Machine Learning/Deep Learning) để tự động hóa quá trình nhận dạng chậu võng và hiệu chỉnh nhiệt độ thời gian thực từ dữ liệu cảm biến IoT gắn trong đường CHC.
  • Nghiên cứu ảnh hưởng của hiện tượng mỏi nhiệt kết hợp mỏi tải trọng lặp (thermo-mechanical fatigue interaction) đến tuổi thọ còn lại của tấm BTXM sân bay.
  • Mở rộng thuật toán hiệu chỉnh cho các kết cấu mặt đường composite (bê tông nhựa phủ trên mặt đường bê tông xi măng cũ).

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

  • Tác động học thuật: Cung cấp tài liệu tham khảo mẫu mực cho các chương trình đào tạo sau đại học ngành Kỹ thuật Xây dựng Công trình Giao thông; mở ra hướng nghiên cứu chuyên sâu về NDT kết cấu sân bay trong điều kiện biến đổi khí hậu.
  • Chuyển đổi ngành công nghiệp: Định hình lại quy trình kiểm định, bảo trì định kỳ đường CHC, đường lăn tại các cảng hàng không Việt Nam; giúp các đơn vị tư vấn thiết kế, kiểm định áp dụng quy trình chuẩn hóa số liệu HWD theo tiêu chuẩn ICAO/FAA.
  • Hiệu quả kinh tế - xã hội: Việc đánh giá chính xác chỉ số PCN giúp tối ưu hóa kế hoạch bảo trì, tránh việc sửa chữa lớn đường CHC sớm hơn thực tế (tiết kiệm hàng trăm tỷ đồng ngân sách hạ tầng hàng không), đồng thời loại bỏ nguy cơ gián đoạn khai thác các đường bay quốc tế.

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

  • Nghiên cứu sinh và học viên sau đại học: Tiếp cận phương pháp luận nghiên cứu kết hợp giữa mô hình hóa giải tích, thực nghiệm hiện trường và hồi quy thống kê chuyên sâu.
  • Các nhà khoa học và giảng viên chuyên ngành: Kế thừa dữ liệu thực nghiệm nhiệt độ 365 ngày và mô hình truyền nhiệt hiệu chỉnh để phát triển các bài toán cơ học mặt đường nâng cao.
  • Kỹ sư R&D và đơn vị kiểm định sân bay: Sử dụng trực tiếp công thức và bảng tra hiệu chỉnh độ võng HWD về điều kiện chuẩn để lập báo cáo PCN chính xác.
  • Cơ quan quản lý hàng không (CAAV, ACV): Có cơ sở khoa học vững chắc để ban hành quy định kỹ thuật, quy trình kiểm tra sức chịu tải đường CHC bằng thiết bị gia tải động HWD.

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ì? Trả lời: Đóng góp lý thuyết độc đáo nhất là việc thiết lập và giải thành công phương trình vi phân truyền nhiệt Fourier một chiều có tham số hiệu chỉnh động ($a, K$) dựa trên chuỗi số liệu thực nghiệm liên tục 1 năm tại Việt Nam, từ đó tích hợp định lượng sự phân bố nhiệt độ phi tuyến vào phương pháp tính toán ngược của Westergaard-Winkler cho mặt đường BTXM sân bay dày $30\text{--}40\text{ cm}$.

2. Điểm đổi mới phương pháp luận so với các nghiên cứu trước đây? Trả lời: Trước đây, các nghiên cứu trong nước chủ yếu khảo sát bằng thiết bị FWD với tải trọng nhỏ dành cho đường ô tô và chỉ dừng lại ở chỉ số PCN đơn lẻ. Luận án đã sử dụng thiết bị HWD PRIMAX 2500 công suất lớn (>150 kN), đồng thời phát triển khung đánh giá toàn diện hai nhóm thông số: nhóm thông số làm việc ($d_i, k, E_{bt}, LTE, PCN$) và nhóm thông số tính toán thiết kế ($x, M_u, \sigma_s$).

3. Phát hiện nào bất ngờ nhất từ dữ liệu thực nghiệm? Trả lời: Phát hiện bất ngờ nhất là sự biến động cực lớn của chỉ số truyền lực khe nối $LTE$ (từ <20% lên >80%) và sự sụt giảm biểu kiến của hệ số phản lực nền $k$ (tới 40-60%) chỉ trong chu kỳ 24 giờ do sự kết hợp giữa hiệu ứng uốn vồng nhiệt ($\Delta t$) và sự giãn nở thể tích đóng mở khe nối (aggregate interlock), dù cấu trúc vật liệu nền móng không thay đổi.

4. Nghiên cứu có cung cấp quy trình tái lặp (Replication Protocol) không? Trả lời: Có. Luận án cung cấp đầy đủ thông số kỹ thuật thiết bị (EXTECH SDL200, Sweco HWD PRIMAX 2500), tọa độ khoảng cách 10 đầu đo geophone ($0\text{--}210\text{ cm}$), cấu trúc khai báo đầu vào phần mềm COMFAA, các bảng số liệu nhiệt độ $K_{tb}$, và thuật toán giải phương trình vi phân, cho phép tái lập quy trình trên bất kỳ cảng hàng không nào.

5. Định hướng nghiên cứu 10 năm tới được phác thảo như thế nào? Trả lời: Xây dựng hệ thống bản sao số (Digital Twin) kết cấu mặt đường sân bay tích hợp cảm biến nhiệt-ẩm-biến dạng sợi quang (FBG) theo dõi sức khỏe công trình theo thời gian thực và tự động cập nhật chỉ số PCN thông báo bay.

Kết luận

  1. Hoàn thiện thành công phương trình truyền nhiệt Fourier phi tuyến theo chiều sâu tấm BTXM sân bay, xác định chính xác các tham số nhiệt vật lý ($a, K, z_{\max}$) đặc thù cho khí hậu miền Bắc Việt Nam.
  2. Đo đạc và phân tích toàn diện chuỗi dữ liệu độ võng động HWD PRIMAX 2500 trên các sân bay quốc tế (Vân Đồn, Nội Bài), làm rõ cơ chế uốn vồng nhiệt và cơ chế chèn móc cốt liệu tác động đến chậu võng.
  3. Định lượng hóa ảnh hưởng của nhiệt độ bề mặt ($t_m$), nhiệt độ trung bình ($t_{tb}$) và chênh lệch nhiệt độ ($\Delta t$) đến hệ thống thông số làm việc ($d_0, k, E_{bt}, LTE, PCN$) và thông số thiết kế kết cấu ($x, M_u, \sigma_s$).
  4. Thiết lập thành công mô hình toán học và khung bảng tra hiệu chỉnh kết quả thí nghiệm HWD từ điều kiện nhiệt độ thực tế về điều kiện chuẩn ($20^\circ\text{C}$).
  5. Đề xuất khung thời gian đo đạc tối ưu tại hiện trường, góp phần hoàn thiện tiêu chuẩn kỹ thuật kiểm định mặt đường sân bay theo chuẩn mực ICAO.
  6. Mở ra ba hướng nghiên cứu học thuật mới: Mô hình hóa tương tác nhiệt - ẩm - cơ học 3D; Đánh giá mỏi kết cấu dưới tác động nhiệt động lực học; và Ứng dụng AI/IoT trong giám sát sức khỏe mặt đường sân bay thông minh.