Tổng quan về luận án

Sự phát triển của hạ tầng giao thông hiện đại đặt ra yêu cầu ngày càng cao đối với các kết cấu bê tông cốt thép (BTCT) khối lớn trong công trình cầu như mố, trụ và bệ móng. Khác với các công trình thủy công hay đập bê tông đầm lăn sử dụng bê tông mác thấp với lượng xi măng khiêm tốn, các bộ phận công trình cầu hiện nay thường đòi hỏi bê tông cường độ cao (từ cấp C25, C30 đến C40) với hàm lượng xi măng poóclăng vượt ngưỡng $400\text{ kg/m}^3$. Hệ quả tất yếu của quá trình này là sự phát sinh một lượng nhiệt thủy hóa khổng lồ trong 72 giờ đầu sau khi đổ, dẫn đến hiện tượng chênh lệch nhiệt độ lớn giữa tâm khối đổ và bề mặt biên, tạo ra ứng suất kéo vượt quá cường độ chịu kéo non tuổi của bê tông và gây ra nứt phi kết cấu nghiêm trọng.

Trong thực tiễn kỹ thuật, khoảng trống nghiên cứu (research gap) cốt lõi tồn tại ở việc hầu hết các mô hình tính toán truyền thống đều xem kết cấu bê tông khối lớn là một môi trường đẳng hướng, thuần túy đồng nhất mà bỏ qua sự hiện diện của mật độ cốt thép dày đặc tại lớp vỏ bảo vệ biên. Lớp cốt thép này không chỉ làm thay đổi cục bộ hệ số dẫn nhiệt mà còn tạo ra mức độ kiềm chế biến dạng cơ học phức tạp. Luận án tiến sĩ kỹ thuật với đề tài "Ứng dụng lý thuyết đồng nhất hóa để phân tích trạng thái phân bố nhiệt độ và ứng suất do nhiệt thủy hóa xi măng trong bê tông cốt thép công trình cầu" do Nghiên cứu sinh Nguyễn Xuân Lam thực hiện tại Trường Đại học Giao thông Vận tải dưới sự hướng dẫn của PGS.TS Nguyễn Ngọc Long và PGS.TS Nguyễn Duy Tiến đã giải quyết trực tiếp khoảng trống học thuật này.

Nghiên cứu tập trung vào các câu hỏi nghiên cứu và giả thuyết khoa học cụ thể:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để mô hình hóa chính xác ảnh hưởng của cốt thép biên đến quá trình truyền nhiệt và phân bố ứng suất mà không làm bùng nổ khối lượng tính toán trong phương pháp phần tử hữu hạn (FEM)?
  • Câu hỏi nghiên cứu 2 (RQ2): Quy luật phát triển nhiệt lượng đoạn nhiệt của các cấp phối bê tông C30 và C35 điển hình dùng cho công trình cầu tại Việt Nam diễn ra như thế nào theo thời gian?
  • Câu hỏi nghiên cứu 3 (RQ3): Mức độ tương quan giữa kết quả mô phỏng số theo lý thuyết đồng nhất hóa đa cấp độ với số liệu quan trắc thực nghiệm nhiệt độ và ứng suất ngoài hiện trường của thân trụ cầu đạt độ tin cậy bao nhiêu?
  • Giả thuyết khoa học 1 (H1): Việc áp dụng lý thuyết đồng nhất hóa trên phần tử thể tích đại diện (Representative Elementary Volume - REV) cho phép quy đổi chính xác lớp vỏ BTCT thành một lớp vật liệu tương đương có hệ số dẫn nhiệt $k_{eq}$ và mô đun đàn hồi $E_{eq}$ dị hướng biến thiên theo thời gian.
  • Giả thuyết khoa học 2 (H2): Đỉnh ứng suất kéo nhiệt trên bề mặt thân trụ cầu xuất hiện sớm (trong khoảng 60 giờ đầu), trước khi bê tông đạt đủ cường độ thiết kế, và là nguyên nhân chính khởi phát các vết nứt bề mặt.

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 lý thuyết đồng nhất hóa vi-vĩ mô (Multi-scale Homogenization Theory), phương trình vi phân truyền nhiệt Fourier phi tuyến có xét nguồn nhiệt đoạn nhiệt theo mô hình của Hiệp hội Kỹ sư Dân dụng Hoa Kỳ (ASCE), và cơ học phá hủy bê tông non tuổi có xét đến từ biến và co ngót theo Fib Model Code 2010. Phạm vi nghiên cứu bao gồm thực nghiệm nhiệt đoạn nhiệt trong phòng thí nghiệm trên các mẫu bê tông C30 và C35, mô phỏng số thông qua việc tự phát triển bộ công cụ phần mềm TCon1 và TCon2, kết hợp đối chứng với số liệu quan trắc thực tế tại thân trụ cầu kích thước lớn ngoài hiện trường.

Literature Review và Positioning

Vấn đề kiểm soát nhiệt thủy hóa và nứt do nhiệt trong bê tông khối lớn đã thu hút sự quan tâm của nhiều thế hệ nhà khoa học trên thế giới. Dòng nghiên cứu nền tảng về nhiệt hóa học và động học thủy hóa bắt nguồn từ các công trình kinh điển của Baikov và Rebinder, phân chia quá trình thủy hóa xi măng thành 5 giai đoạn rõ rệt: hòa tan, ngủ đông, gia tốc, giảm tốc và ổn định. Trong giai đoạn gia tốc (thường từ 2 đến 10 giờ sau khi trộn), khoáng alit ($C_3S$) và aluminat tricanxit ($C_3A$) phản ứng mãnh liệt với nước, tạo ra các tinh thể hydrosilicat canxi (C-S-H) và ettringite ($3\text{CaO}\cdot\text{Al}_2\text{O}_3\cdot 3\text{CaSO}_4\cdot 31\text{-}32\text{H}_2\text{O}$), giải phóng một lượng nhiệt đơn vị cực lớn.

Ở bình diện quốc tế, các mô hình nhiệt - cơ cho bê tông tuổi sớm đã được mở rộng mạnh mẽ:

  • Ulm và Coussy (2001) đã đề xuất phương pháp đánh giá đặc trưng khuếch tán nhiệt thủy hóa dựa trên nhiệt động lực học môi trường liên tục, gắn kết mức độ thủy hóa với sự tiến triển của cơ chế sinh nhiệt nội tại.
  • De Schutter và Taerwe (1996) phát triển khái niệm độ dày tương đương ($d_{eq} = \gamma_a \cdot M$) để chuẩn hóa quan hệ giữa hình học kết cấu và gradient nhiệt độ trong các cấu kiện bê tông khối lớn có hình dạng phức tạp.
  • Tiêu chuẩn ACI 207 và ACI 116R đưa ra các chỉ dẫn thực hành về bê tông khối lớn, định nghĩa khối đổ dựa trên kích thước đủ lớn gây tích tụ nhiệt, trong khi FIB Model Code 2010 cung cấp khung phân tích phi tuyến cơ bản để dự báo chỉ số nứt do nhiệt kết hợp co ngót và từ biến.

Tuy nhiên, trong y văn tồn tại hai luồng quan điểm tranh luận sâu sắc:

  1. Luồng quan điểm thứ nhất cho rằng việc kiểm soát nứt chỉ cần tối ưu hóa cấp phối (giảm lượng xi măng xuống dưới $350\text{ kg/m}^3$, thay thế tro bay từ 25-30%) và áp dụng các biện pháp làm mát thông thường.
  2. Luồng quan điểm thứ hai chứng minh rằng đối với các kết cấu trụ cầu hiện đại chịu tải trọng nặng, việc tăng cường cốt thép biên tạo ra độ kiềm chế cục bộ rất lớn; nếu không xét đến sự tương tác nhiệt - cơ giữa lớp cốt thép này và ma trận bê tông thì các biện pháp giảm nhiệt cấp phối vẫn không ngăn ngừa được nứt bề mặt.

Để định vị nghiên cứu, luận án đã phân tích và đối chiếu với nhiều trường hợp sự cố nứt kết cấu bê tông khối lớn kinh điển trên thế giới:

  • Đập tràn Nam Mỹ: Kích thước khối đập $35 \times 15 \times 7{,}5\text{ m}$, sử dụng $350\text{ kg xi măng/m}^3$ ($W/C = 0{,}52$). Nhiệt độ môi trường thi công $32^\circ\text{C}$, nhiệt độ lõi đập đạt đỉnh $62{,}5^\circ\text{C}$, tạo ra chênh lệch nhiệt độ bề mặt lên tới $42^\circ\text{C}$, gây ra các vết nứt rộng từ $0{,}2$ đến $1{,}2\text{ mm}$ lan sâu vào thân đập từ 10 đến 100 ngày sau khi đúc.
  • Đập Itaipu (Brazil/Paraguay): Khối đổ thể tích $13.000\text{ m}^3$ gồm 83 phân đoạn kích thước $17 \times 35 \times 85\text{ m}$. Dù nhiệt độ bê tông tươi được khống chế ở $7^\circ\text{C}$ và hàm lượng xi măng chỉ từ $108 - 169\text{ kg/m}^3$, các vết nứt nông dài $10 - 20\text{ m}$, rộng $0{,}3 - 0{,}9\text{ mm}$ vẫn xuất hiện tại 34/47 khối do biến dạng kiềm chế và từ biến non tuổi.
  • Trụ cầu tại Cộng hòa Séc: Kích thước $5 \times 5 \times 4{,}57\text{ m}$, cấp bê tông C35/45 ($400\text{ kg xi măng/m}^3$, $40\text{ kg tro bay/m}^3$). Bố trí 3 lớp cốt thép dọc đường kính $32\text{ mm}$ (tỷ lệ $0{,}8%$). Nhiệt độ lõi tăng từ $20^\circ\text{C}$ lên $61^\circ\text{C}$ sau 9 ngày. Sau 1 năm, xuất hiện các vết nứt dọc rộng $0{,}1 - 0{,}5\text{ mm}$ xuyên qua cả lớp cốt thép gia cố do mật độ thép đai không đủ kiềm chế ứng suất nhiệt cục bộ.
  • Hệ thống móng trụ cầu tại Bắc Carolina (Mỹ): Khảo sát các cầu Oak Island (OIB), Sunset Beach (SBB) và Wilmington Bypass (WBB) với lượng chất kết dính $410\text{ kg/m}^3$. Tại cầu WBB, chênh lệch nhiệt độ cực đại đạt $55^\circ\text{C}$ tại thời điểm 156 giờ, gây nứt nghiêm trọng theo cả phương dọc và ngang bất chấp việc đã dùng ván khuôn thép cách nhiệt.
  • Trụ cầu cạn tại Ý: Thân trụ đường kính $3{,}5 - 4{,}5\text{ m}$, dùng $330\text{ kg xi măng/m}^3$. Tháo ván khuôn sau 1 ngày dẫn đến nhiệt độ lõi $65^\circ\text{C}$ tạo gradient nhiệt độ đột ngột, làm nứt bề mặt với độ mở rộng $0{,}05 - 0{,}3\text{ mm}$.

Tại Việt Nam, các nghiên cứu của các tác giả trong nước trước đây thường chỉ dừng lại ở việc khảo sát các khối bê tông thuần túy không có cốt thép (khối lập phương kích thước từ $2\text{ m}$ đến $5\text{ m}$) hoặc chỉ xét ảnh hưởng của phụ gia khoáng tro bay trên mẫu phòng thí nghiệm. Việc luận án định vị trọng tâm vào kết cấu trụ cầu BTCT thực tế, kết hợp lý thuyết đồng nhất hóa đa cấp độ để giải bài toán nhiệt thủy hóa xi măng, là một bước tiến học thuật đột phá, lấp đầy khoảng trống giữa lý thuyết truyền nhiệt thuần túy và ứng xử cơ học phi tuyến của kết cấu BTCT non tuổi.

Công trình / Nghiên cứu Kích thước cấu kiện Cấp phối & Hàm lượng xi măng Chênh lệch nhiệt độ cực đại ($\Delta T_{max}$) Dạng hư hỏng / Bề rộng vết nứt ghi nhận
Đập tràn Nam Mỹ $35 \times 15 \times 7{,}5\text{ m}$ $350\text{ kg/m}^3$, $W/C = 0{,}52$ $42^\circ\text{C}$ (Lõi đạt $62{,}5^\circ\text{C}$) Vết nứt bề mặt rộng $0{,}2 - 1{,}2\text{ mm}$
Đập Itaipu (Brazil) $17 \times 35 \times 85\text{ m}$ $108 - 169\text{ kg/m}^3$ + Tro bay Lõi $36^\circ\text{C}$, mặt $22^\circ\text{C}$ Nứt dọc tại 34/47 khối, rộng $0{,}3 - 0{,}9\text{ mm}$
Trụ cầu CH Séc $5 \times 5 \times 4{,}57\text{ m}$ $400\text{ kg/m}^3$ (Bê tông C35/45) Tăng từ $20^\circ\text{C}$ lên $61^\circ\text{C}$ sau 9 ngày Vết nứt dọc xuyên qua 3 lớp thép D32 ($0{,}1 - 0{,}5\text{ mm}$)
Móng cầu WBB (Mỹ) $28{,}8 \times 18{,}5\text{ m}$ $410\text{ kg/m}^3$ chất kết dính $55^\circ\text{C}$ (Lõi đạt $78^\circ\text{C}$ sau 120h) Nứt dọc và ngang nghiêm trọng
Trụ cầu Vĩnh Tuy (VN) Trụ đặc kích thước lớn Bê tông mác cao $> 400\text{ kg/m}^3$ Vượt ngưỡng $20^\circ\text{C}$ (TCVN 9341) Nứt nhiệt bề mặt trong giai đoạn non tuổi
Nghiên cứu của Luận án Thân trụ cầu BTCT thực tế Bê tông C30 và C35 Đo đạc & mô phỏng theo TCon2 Xác định đỉnh ứng suất kéo tại 60 giờ

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

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

Luận án đã mở rộng lý thuyết đồng nhất hóa môi trường liên tục (Homogenization Theory for Continuum Media) từ bài toán cơ học thuần túy sang bài toán ghép nối nhiệt - cơ (coupled thermo-mechanical problem) trong môi trường bê tông cốt thép non tuổi biến thiên theo thời gian. Đóng góp lý thuyết cốt lõi bao gồm việc xây dựng mô hình phần tử thể tích đặc trưng (Representative Elementary Volume - REV) chứa ma trận bê tông và các thanh cốt thép biên được chia lưới tam giác hữu hạn.

Thông qua việc thiết lập các điều kiện biên chuyển vị tuần hoàn và điều kiện biên biến dạng đơn vị: $$\varepsilon = \varepsilon_{11} e_1 \otimes e_1, \quad \varepsilon = \varepsilon_{22} e_2 \otimes e_2, \quad \varepsilon = \varepsilon_{12} (e_1 \otimes e_2 + e_2 \otimes e_1)$$

Mô hình lý thuyết này chứng minh rằng lớp vỏ BTCT dày từ 50 mm đến hàng trăm milimét có thể được chuyển đổi toán học tương đương thành một lớp vật liệu liên tục đồng nhất có các thông số dẫn nhiệt ($k_{xx}^{eq}, k_{yy}^{eq}$), nhiệt dung riêng ($C_{eq}$) và ma trận độ cứng đàn hồi ($E_x^{eq}, E_y^{eq}, \nu_{yx}^{eq}$) dị hướng. Điều này thay đổi hoàn toàn giả định trước đây vốn coi lớp bảo vệ bê tông có tính chất nhiệt tương đồng với phần lõi không gia cốt thép.

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

Khung phân tích của luận án tích hợp đồng thời ba hệ thống lý thuyết bổ trợ:

  1. Lý thuyết truyền nhiệt phi tuyến trong môi trường xốp: Phương trình truyền nhiệt tổng quát dạng đạo hàm riêng: $$\frac{\partial}{\partial x}\left(k_x \frac{\partial T}{\partial x}\right) + \frac{\partial}{\partial y}\left(k_y \frac{\partial T}{\partial y}\right) + q(t) = \rho C \frac{\partial T}{\partial t}$$ Trong đó nguồn nhiệt phát sinh $q(t)$ được biểu diễn theo hàm mức độ thủy hóa $\alpha(t)$ chuẩn hóa bởi ASCE: $$q(t) = \rho \cdot C \cdot T_{\max} \cdot \alpha(t) \cdot \frac{1}{24} e^{-\alpha \frac{t}{24}}$$
  2. Lý thuyết biến dạng nhiệt và tự do kiềm chế: Kết hợp tác động của gradient nhiệt độ nội tại với các điều kiện biên ngăn cản chuyển vị pháp tuyến của thân trụ mố cầu.
  3. Mô hình đàn-nhớt phi tuyến (Viscoelasticity) cho bê tông non tuổi: Tích hợp các hàm ảnh hưởng của từ biến (creep compliance) và co ngót tự sinh (autogenous shrinkage) theo thời gian, giúp phản ánh hiện tượng chùng ứng suất (stress relaxation) khi khối bê tông hạ nhiệt.

Điều kiện biên của khung phân tích được xác định tường minh: biên tiếp xúc ván khuôn và môi trường không khí được mô hình hóa qua điều kiện biên loại 3 (Robin boundary condition) với hệ số trao đổi nhiệt đối lưu bề mặt biến thiên theo vận tốc gió và nhiệt độ khí quyển thực đo.

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ lập trường triết học thực chứng (positivism) kết hợp hiện thực phản biện (critical realism), xây dựng một thiết kế nghiên cứu đa cấp độ (Multi-level Research Design) từ cấp độ vi mô (thành phần khoáng clinker xi măng), cấp độ trung mô (phần tử thể tích đặc trưng REV của bê tông và cốt thép) đến cấp độ vĩ mô (kết cấu thân trụ cầu kích thước thật).

Phương pháp nghiên cứu kết hợp chặt chẽ giữa:

  • Nghiên cứu thực nghiệm trong phòng: Xác định thành phần hóa học, thành phần khoáng vật của xi măng và đo nhiệt độ tăng đoạn nhiệt của hỗn hợp bê tông.
  • Mô phỏng số và phát triển phần mềm: Lập trình thuật toán đồng nhất hóa và thuật toán PTHH phân tích trường nhiệt - ứng suất.
  • Thực nghiệm hiện trường: Đo đạc nhiệt độ trực tiếp tại hiện trường bằng hệ thống cảm biến nhiệt độ gắn sâu trong các phân đoạn thân trụ cầu.

Quy trình nghiên cứu rigorous

Quy trình thực nghiệm và thu thập số liệu được chuẩn hóa nghiêm ngặt qua các giai đoạn:

  1. Phân tích hóa lý xi măng: Xác định hàm lượng các oxit chủ yếu bằng phương pháp huỳnh quang tia X: $\text{CaO}$ ($63-66%$), $\text{SiO}_2$ ($21-24%$), $\text{Al}_2\text{O}_3$ ($4-8%$), $\text{Fe}_2\text{O}_3$ ($2-4%$), cùng hàm lượng khoáng vật clinker chính: Alit ($C_3S$: $45-60%$), Belit ($C_2S$: $20-30%$), Aluminat tricanxit ($C_3A$: $4-12%$), Feroaluminat tetracanxit ($C_4AF$: $10-12%$).
  2. Thí nghiệm nhiệt lượng đoạn nhiệt (Adiabatic Calorimetry): Chế tạo mẫu bê tông cấp C30 và C35 theo đúng cấp phối thực tế công trình cầu. Mẫu sau khi trộn được đặt ngay vào thùng chứa đoạn nhiệt chuyên dụng có kết nối hệ thống cảm biến nhiệt và máy tính theo dõi tự động để ghi nhận đường cong gia tăng nhiệt độ $T_{ad}(t)$ liên tục từ 0 đến 168 giờ.
  3. Phát triển chương trình phần mềm TCon1 và TCon2:
    • Chương trình TCon1: Nhập dữ liệu hình học REV, đường kính cốt thép ($D18, D20, D25, D32\text{ mm}$), bề dày lớp bảo vệ ($50\text{ mm}$), giải bài toán phần tử hữu hạn lưới tam giác để xuất ra các bảng hệ số dẫn nhiệt tương đương $k_{eq}$, mô đun đàn hồi tương đương $E_x, E_y$, hệ số Poisson $\nu_{yx}$, và nhiệt dung riêng tương đương $C_{eq}$.
    • Chương trình TCon2: Đọc kết quả từ TCon1 và đường cong nhiệt đoạn nhiệt thực nghiệm, gán điều kiện biên nhiệt môi trường, hàm đối lưu, hàm từ biến và co ngót để tính toán phân bố nhiệt độ và ứng suất nhiệt 3D cho toàn khối thân trụ theo các bước thời gian.
  4. Quan trắc thực tế ngoài hiện trường: Thiết lập mạng lưới đầu đo nhiệt kế điện tử (Thermometer) trên mặt cắt ngang thân trụ cầu thực tế. Các cảm biến được bố trí dày đặc theo cả hai trục cạnh dài và cạnh ngắn, từ tâm trụ ra tới sát bề mặt lớp vỏ BTCT.

Data và phân tích

Dữ liệu thực nghiệm được xử lý qua mô hình hồi quy phi tuyến để khớp đường cong mức độ thủy hóa $\alpha(t)$. Kiểm tra độ tin cậy và tính vững (robustness checks) được thực hiện bằng cách so sánh đối chuẩn chéo (cross-validation) giữa trường nhiệt độ tính toán bởi mô hình số thuần túy bê tông (mô hình truyền thống), mô hình BTCT đồng nhất hóa (TCon2), và dữ liệu quan trắc thực tế tại các mốc thời gian quan trọng: 12 giờ, 24 giờ, 48 giờ, 60 giờ, 120 giờ và 1000 giờ.

      Nhiệt độ (°C)
          0     20     40     60     80    100    120

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

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

  1. Quy luật dẫn nhiệt tương đương phụ thuộc mạnh vào đường kính cốt thép và lớp bảo vệ: Kết quả tính toán từ TCon1 chứng minh rằng hệ số dẫn nhiệt tương đương $k_{eq}$ của lớp vỏ BTCT tăng đáng kể khi đường kính cốt thép tăng từ D18 đến D32. Lớp cốt thép biên đóng vai trò như một "cầu dẫn nhiệt" cục bộ làm tăng tốc độ thất thoát nhiệt từ vùng biên ra môi trường ngoài, khiến gradient nhiệt độ tại vùng sát bề mặt trở nên dốc hơn nhiều so với dự đoán của mô hình bê tông trơn.
  2. Đỉnh ứng suất kéo nhiệt nguy hiểm nhất xuất hiện tại thời điểm 60 giờ: Mô phỏng TCon2 và dữ liệu hiện trường khẳng định rằng trong khi nhiệt độ tâm khối đổ có thể đạt đỉnh muộn hơn, thì ứng suất kéo lớn nhất do nhiệt thủy hóa tại bề mặt bê tông lại đạt giá trị cực đại ở mốc xấp xỉ 60 giờ sau khi đổ. Tại thời điểm này, cường độ chịu kéo của bê tông non tuổi phát triển chưa hoàn thiện, khiến chỉ số an toàn chống nứt ($I_{cr} = f_{ct}(t) / \sigma_t$) tụt xuống mức thấp nhất, giải thích nguyên nhân tại sao các vết nứt chân chim và nứt dọc thân trụ thường hình thành ngay trong 3 ngày đầu.
  3. Hiệu ứng giải tỏa ứng suất nhờ từ biến sau 1000 giờ: Phân tích trạng thái ứng suất dài hạn đến 1000 giờ cho thấy hiện tượng chùng ứng suất do từ biến đã giúp làm giảm đáng kể biên độ ứng suất kéo dư. Tuy nhiên, nếu các vết nứt tế vi đã bị kích hoạt ở giai đoạn 60 giờ, chúng sẽ đóng vai trò như các tập trung ứng suất (stress concentrators) làm suy giảm tính toàn vẹn kết cấu khi chịu tải trọng khai thác sau này.
  4. Sự không phù hợp của quy chuẩn chung khi áp dụng cho trụ cầu bê tông mác cao: Tiêu chuẩn hiện hành TCVN 9341:2012 quy định độ chênh lệch nhiệt độ $\Delta T \le 20^\circ\text{C}$ và mô đun độ chênh lệch nhiệt độ $M_T \le 50^\circ\text{C/m}$. Luận án chứng minh rằng với bê tông cấp C30-C35 có hàm lượng xi măng $> 400\text{ kg/m}^3$, nếu chiều dày thân trụ lớn hơn $1{,}6\text{ m}$, gradient nhiệt độ tự nhiên sẽ dễ dàng vượt ngưỡng $50^\circ\text{C/m}$ nếu chỉ áp dụng các biện pháp bảo dưỡng thông thường không có kiểm soát nhiệt chủ động.

Implications đa chiều

  • Về mặt học thuật và lý thuyết: Xác lập cơ sở khoa học hoàn chỉnh cho phương pháp mô hình hóa phân lớp đồng nhất (homogenized layered approach), loại bỏ sự cần thiết phải chia lưới phần tử hữu hạn quá mịn quanh từng thanh cốt thép, tiết kiệm tới 80% tài nguyên tính toán nhưng vẫn bảo toàn độ chính xác về ứng xử nhiệt - cơ vi mô.
  • Về mặt công nghệ và thực tiễn thi công:
    • Cung cấp cơ sở định lượng để tối ưu hóa chiều dày khối đổ và lượng xi măng theo từng mác bê tông (M200, M250, M300, M350, M400) nhằm đảm bảo thân trụ không bị nứt.
    • Đưa ra khuyến nghị bắt buộc về thời gian giữ ván khuôn cách nhiệt: không được tháo ván khuôn trong vòng 3 đến 5 ngày đầu đối với các thân trụ kích thước $> 1{,}5\text{ m}$ để tránh sốc nhiệt bề mặt.
    • Ứng dụng phụ gia khoáng tro bay hoạt tính cao nhằm kéo dài giai đoạn ngủ đông và giảm tốc độ tỏa nhiệt cực đại của khoáng alit trong giai đoạn gia tốc.
  • Về mặt chính sách và tiêu chuẩn xây dựng: Đóng góp dữ liệu thực chứng quan trọng phục vụ công tác sửa đổi, bổ sung Tiêu chuẩn thiết kế cầu đường bộ TCVN 11823:2017 và Tiêu chuẩn thi công bê tông khối lớn TCVN 9341:2012, bổ sung các điều khoản riêng biệt cho kết cấu phần dưới công trình cầu giao thông.

Limitations và Future Research

Mặc dù đạt được những đóng góp đột phá, luận án thẳng thắn ghi nhận một số giới hạn nghiên cứu:

  1. Giả định dính bám hoàn hảo (Perfect Bond): Mô hình REV đồng nhất hóa giả định dính bám hoàn toàn giữa bê tông và cốt thép, chưa xét đến hiện tượng trượt dính bám (bond-slip) và vi nứt cục bộ tại đới chuyển tiếp tiếp xúc (interfacial transition zone - ITZ) dưới tác động của chu trình nhiệt độ biến thiên liên tục.
  2. Phạm vi cấp phối và nguồn vật liệu: Thực nghiệm nhiệt đoạn nhiệt mới tập trung khảo sát trên các cấp phối C30 và C35 sử dụng xi măng poóclăng phổ biến tại thị trường miền Bắc Việt Nam, chưa bao quát toàn diện các loại xi măng xỉ lò cao hoặc xi măng bền sunfat dùng cho công trình biển.
  3. Bỏ qua ghép nối với tải trọng cơ học thi công: Nghiên cứu tập trung giải quyết trường ứng suất thuần túy do nhiệt thủy hóa, từ biến và co ngót, chưa xét đến sự cộng hưởng đồng thời của tải trọng cơ học phát sinh trong quá trình thi công đốt trụ tiếp theo hoặc tải trọng gió bão.

Các hướng nghiên cứu tương lai được đề xuất bao gồm:

  • Mở rộng ứng dụng lý thuyết đồng nhất hóa cho bê tông siêu tính năng (UHPC) và bê tông tự lèn (SCC) có hàm lượng chất kết dính siêu cao dùng trong các cầu dây văng nhịp lớn.
  • Phát triển mô hình REV ngẫu nhiên (Stochastic REV) có xét đến sự phân bố cốt thép không đều và các khuyết tật rỗng tế vi.
  • Ứng dụng cảm biến sợi quang FBG (Fiber Bragg Grating) kết hợp thuật toán học máy để quan trắc và dự báo trực tuyến theo thời gian thực (Real-time Digital Twin) trường nhiệt - ứng suất của mố trụ cầu ngay khi đổ bê tông.

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

Nghiên cứu của luận án tạo ra chuỗi giá trị tác động sâu rộng:

  • Tác động học thuật: Mở ra hướng nghiên cứu liên ngành cơ học vật liệu tính toán và kỹ thuật xây dựng cầu hầm tại Việt Nam; cung cấp công cụ mã nguồn TCon1 và TCon2 làm tiền đề cho các nghiên cứu sinh và học viên cao học tiếp tục phát triển mô phỏng truyền nhiệt phi tuyến đa tỷ lệ.
  • Chuyển đổi thực tiễn ngành xây dựng cầu: Giúp các nhà thầu thi công lớn (như Cienco, Đèo Cả, Tổng công ty Thăng Long...) kiểm soát hoàn toàn rủi ro nứt nhiệt tại các dự án cầu vượt sông, cầu cạn cao tốc, giảm thiểu chi phí sửa chữa vết nứt sau thi công vốn có thể chiếm từ $3 - 5%$ giá trị gói thầu kết cấu phần dưới.
  • Lợi ích kinh tế - xã hội: Kéo dài tuổi thọ khai thác của công trình cầu từ 50 lên 100 năm thông qua việc triệt tiêu các vết nứt nhiệt non tuổi – con đường chính dẫn truyền ion clorua gây rỉ cốt thép sớm trong môi trường nhiệt đới ẩm gió mùa và xâm thực ven biển Việt Nam.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận phương pháp luận đồng nhất hóa vật liệu hoàn chỉnh, có thể kế thừa thuật toán và số liệu thực nghiệm để công bố các công trình quốc tế về cơ học công trình.
  • Kỹ sư Tư vấn Thiết kế: Sử dụng các bảng tra tương quan giữa mác bê tông ($M$), chiều dày thân trụ ($a$) và lượng xi măng tối ưu do luận án thiết lập để lựa chọn kích thước hình học và hàm lượng cốt thép hợp lý ngay từ giai đoạn thiết kế cơ sở.
  • Kỹ sư Giám sát và Nhà thầu Thi công: Có được cẩm nang chỉ dẫn kỹ thuật định lượng về thời điểm nguy hiểm nhất (mốc 60 giờ), kiểm soát chặt chẽ nhiệt độ bê tông tươi ban đầu, áp dụng hệ thống bảo dưỡng cách nhiệt bề mặt đúng thời điểm nhằm đáp ứng các tiêu chuẩn kỹ thuật nghiêm ngặt.
  • Cơ quan Quản lý Nhà nước (Bộ Giao thông Vận tải, Bộ Xây dựng): Có cơ sở khoa học thực chứng để ban hành các quy chuẩn, tiêu chuẩn kỹ thuật quốc gia về kiểm soát nứt phi kết cấu cho hạ tầng giao thông trọng điểm.

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à nó mở rộng lý thuyết nào?

Luận án đã mở rộng lý thuyết đồng nhất hóa đa cấp độ (Multi-scale Homogenization Theory) sang bài toán truyền nhiệt và ứng suất phi tuyến của bê tông cốt thép non tuổi. Đột phá độc đáo nhất là việc lượng hóa thành công lớp vỏ BTCT thành một môi trường đồng nhất tương đương có hệ số dẫn nhiệt $k_{eq}$ và mô đun đàn hồi $E_{eq}$ dị hướng biến thiên theo thời gian thực dựa trên cấu trúc phần tử thể tích đặc trưng (REV), giải quyết triệt để bài toán bùng nổ khối lượng tính toán khi mô phỏng mạng lưới cốt thép dày đặc trong kết cấu cầu lớn.

2. Đổi mới phương pháp luận của luận án so với các nghiên cứu quốc tế trước đây thể hiện ở điểm nào?

So với nghiên cứu của Ulm & Coussy (2001) chỉ xét khối bê tông trơn đồng nhất, hay De Schutter & Taerwe (1996) chỉ dùng hệ số hình học đơn giản $d_{eq}$, phương pháp luận của luận án tích hợp chu trình khép kín ba bước: (1) Thí nghiệm đo nhiệt lượng đoạn nhiệt thực tế cho cấp phối C30/C35; (2) Đồng nhất hóa REV lớp vỏ BTCT qua phần mềm tự lập trình TCon1; (3) Giải bài toán trường nhiệt - ứng suất 3D phi tuyến có xét từ biến/co ngót qua phần mềm TCon2 và kiểm chứng trực tiếp bằng chuỗi số liệu quan trắc nhiệt kế Thermometer tại hiện trường thân trụ cầu thật.

3. Phát hiện nào trong dữ liệu mang tính bất ngờ nhất và lời giải thích cơ học là gì?

Phát hiện bất ngờ nhất là đỉnh ứng suất kéo nhiệt nguy hiểm nhất trên bề mặt thân trụ cầu không xuất hiện ở thời điểm nhiệt độ lõi cao nhất, mà xuất hiện cục bộ ở mốc khoảng 60 giờ sau khi đổ. Lời giải thích cơ học nằm ở sự trễ pha giữa tốc độ tích nhiệt bên trong và tốc độ thất thoát nhiệt qua lớp vỏ BTCT dẫn nhiệt cao ra môi trường, tạo ra gradient nhiệt độ bề mặt dốc nhất đúng vào thời điểm bê tông chưa phát triển đủ cường độ chịu kéo non tuổi ($f_{ct}(t)$), dẫn đến nguy cơ nứt cao nhất.

4. Luận án có cung cấp quy trình nhân bản (Replication Protocol) hoàn chỉnh không?

Luận án cung cấp quy trình nhân bản vô cùng chi tiết: từ thông số cấp phối $1\text{ m}^3$ bê tông, bảng thành phần oxit hóa học và khoáng vật clinker xi măng, sơ đồ đấu nối thí nghiệm thùng đoạn nhiệt, hệ phương trình toán học giải tích, sơ đồ khối thuật toán TCon1, TCon2 cho đến tọa độ bố trí các đầu đo Thermometer trên mặt cắt ngang thân trụ, cho phép các phòng thí nghiệm độc lập tái lập hoàn toàn quy trình nghiên cứu.

5. Lộ trình nghiên cứu 10 năm tiếp theo từ nền tảng luận án này nên tập trung vào đâu?

Lộ trình 10 năm tới cần tập trung vào ba trụ cột: (1) Nâng cấp mô hình đồng nhất hóa cho vật liệu bê tông thế hệ mới (UHPC, Geopolymer, bê tông hấp thụ carbon); (2) Tích hợp mô hình cơ học vi nứt dính bám trượt (bond-slip damage mechanics) vào phần tử REV; (3) Tự động hóa thiết kế tối ưu hệ thống ống giải nhiệt cưỡng bức (cooling pipes) thông qua trí tuệ nhân tạo (AI) và tối ưu hóa đa mục tiêu dựa trên nền tảng tính toán nhiệt TCon2.

Kết luận

Luận án tiến sĩ kỹ thuật của Nghiên cứu sinh Nguyễn Xuân Lam đã hoàn thành xuất sắc các mục tiêu nghiên cứu đặt ra với các đóng góp cốt lõi:

  1. Xây dựng thành công cơ sở lý thuyết đồng nhất hóa đa cấp độ cho bài toán nhiệt - cơ trong kết cấu bê tông cốt thép non tuổi công trình cầu.
  2. Phát triển bộ công cụ phần mềm chuyên dụng TCon1 và TCon2, cho phép xác định chính xác các đặc trưng vật liệu tương đương ($k_{eq}, C_{eq}, E_x, E_y, \nu_{yx}$) và phân tích phân bố trường nhiệt độ, ứng suất 3D theo thời gian.
  3. Thiết lập hệ thống dữ liệu thực nghiệm nhiệt lượng đoạn nhiệt chuẩn xác cho các cấp phối bê tông C30 và C35 dùng phổ biến trong kết cấu phần dưới công trình cầu tại Việt Nam.
  4. Xác định tường minh cơ chế hình thành ứng suất nhiệt và đỉnh nguy cơ nứt tại mốc 60 giờ, chứng minh vai trò cầu dẫn nhiệt và kiềm chế cơ học của lớp cốt thép biên.
  5. Cung cấp hệ thống giải pháp công nghệ và biểu đồ tương quan kỹ thuật, giúp các kỹ sư thiết kế và thi công kiểm soát triệt để hiện tượng nứt do nhiệt thủy hóa, nâng cao chất lượng và độ bền vững lâu dài cho các công trình cầu giao thông trọng điểm.