Tổng quan về luận án

Sự phát triển mạnh mẽ của bê tông tính năng cao (High-Performance Concrete - HPC) với tỷ lệ nước/xi măng ($w/c$) thấp mang lại nhiều lợi ích vượt trội về cường độ chịu lực, độ đặc chắc và độ bền kết cấu. Tuy nhiên, các vật liệu này lại đối mặt với thách thức lớn về nứt sớm do co ngót tự sinh (autogeneous shrinkage). Hiện tượng co ngót tự sinh khởi nguồn từ sự tự sấy khô (self-desiccation) và co ngót hóa học (chemical shrinkage / co ngót Le Châtelier với độ giảm thể tích từ 6% đến 8% ở mức thủy hóa hoàn toàn). Khi khung rắn bắt đầu liên kết (percolation), các lỗ rỗng mao dẫn chưa bão hòa hình thành giao diện khum chất lỏng - khí, tạo ra ứng suất căng mao dẫn (capillary tension) cục bộ và gây biến dạng thể tích vĩ mô. Việc dự báo định lượng chính xác co ngót tự sinh ở tuổi sớm vẫn là thách thức lớn đối với cơ học vật liệu xây dựng.

Luận án tiến sĩ số 5881 (2013) tại Trường Bách khoa Liên bang Lausanne (EPFL) của tác giả Đỗ Quang Huy, dưới sự hướng dẫn của GS. Karen Scrivener và GS. Shashank Bishnoi tại Laboratoire des Matériaux de Construction (LMC), mang tên: "Modelling Properties of Cement Paste from Microstructure: Porosity, Mechanical Properties, Creep and Shrinkage". Nghiên cứu đặt mục tiêu xây dựng một khung mô hình vi cơ học số thống nhất (unified micromechanical framework) dựa trên nền tảng vi cấu trúc vector 3D $\mu$ic để mô phỏng đồng thời quá trình thủy hóa, tiến hóa độ xốp, đặc tính đàn hồi, từ biến cơ bản (basic creep) và co ngót tự sinh của hồ xi măng Portland.

Hệ thống câu hỏi và giả thuyết nghiên cứu được xác lập chặt chẽ:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để mô hình hóa chính xác sự tiến hóa kích thước lỗ rỗng và mô phỏng phép đo đo độ xốp xâm nhập thủy ngân (MIP) trên cấu trúc 3D khi xét đến hình thái thực của pha Calcium Silicate Hydrate (C-S-H)?
  • Câu hỏi nghiên cứu 2 (RQ2): Làm thế nào để triệt tiêu các liên kết nhân tạo (artefacts) do quá trình chia lưới phần tử hữu hạn (FEM) gây ra hiện tượng kết nối mạng rắn giả tạo trước thời ninh kết?
  • Câu hỏi nghiên cứu 3 (RQ3): Cơ chế từ biến cơ bản của pha C-S-H biến thiên theo tuổi và mật độ đóng gói ($\eta$) đóng góp như thế nào vào biến dạng co ngót tự sinh vĩ mô?
  • Giả thuyết nghiên cứu 1 (H1): Sự biến thiên mật độ khối của C-S-H ($\rho$ tăng dần theo thời gian đến giá trị tới hạn $2.0\text{ g/cm}^3$, tương ứng $\eta = 0.7$) kết hợp với hiện tượng kết bông (flocculation) hạt ban đầu giải thích chính xác ngưỡng percolation cơ học và mô đun đàn hồi tuổi sớm.
  • Giả thuyết nghiên cứu 2 (H2): Phương pháp chồng chất từ biến (creep-superposition approach) phản ánh bản chất tương tác vi cơ học chính xác hơn mô hình xốp đàn hồi kinh điển (poro-elasticity approach) trong việc dự báo co ngót tự sinh ở các tỷ lệ $w/c$ thấp.

Luận án bao quát phạm vi mô phỏng vi cấu trúc 3D với thể tích tính toán chứa hàng triệu hạt clinker (alite $C_3S$, belite $C_2S$, aluminate $C_3A$, ferrite $C_4AF$), kiểm chứng thực nghiệm tại điều kiện đẳng nhiệt 20°C từ 18 giờ đến giai đoạn trưởng thành.

Literature Review và Positioning

Các nghiên cứu vi cấu trúc xi măng trong hơn ba thập kỷ qua phân hóa thành nhiều trường phái mô hình hóa với các giả định khác nhau:

  • Mô hình HymoStruc (van Breugel, 1995a, 1995b; Koenders & van Breugel, 1997): Dựa trên tiếp cận vector biểu diễn các hạt xi măng hình cầu phát triển các lớp đồng tâm. Nhược điểm cốt lõi của HymoStruc là giả định chỉ có một pha hydrate duy nhất bọc quanh hạt clinker, không phản ánh sự mầm hóa (nucleation) độc lập trong không gian lỗ rỗng và bỏ qua sự tương tác va chạm hạt (impingements) phức tạp.
  • Mô hình CemHyd3D (Bentz & Garboczi, 1991; Bentz, 1997, 2000): Ứng dụng phương pháp rời rạc hóa trên lưới voxel 3D. Dù mô tả được hình học bất kỳ, CemHyd3D chịu giới hạn độ phân giải voxel ($1\text{ }\mu\text{m}$), gây hiện tượng nghẽn mạng (voxel locking) khi chia lưới mịn và không thể mô tả các lỗ rỗng mao dẫn siêu nhỏ hay vi cấu trúc dưới micron.
  • Mô hình Động học hạt tích hợp IPKM (Navi & Pignat, 1996; Pignat et al., 2005): Khởi xướng tiếp cận vector tính toán động học mầm hóa và phát triển pha theo phương trình Avrami, nhưng bị hạn chế nghiêm trọng về thuật toán cấu trúc dữ liệu, chỉ xử lý được khoảng 20.000 hạt, không đủ đại diện cho thể tích vi mô thực tế.

Trong lý thuyết về động lực co ngót tự sinh, các cuộc tranh luận học thuật tập trung vào ba cơ chế:

  1. Áp suất tách (Disjoining pressure) (Ferraris, 1986; Visser, 1998): Chi phối bởi lực van der Waals giữa các bề mặt hydrate cách nhau dưới hai lớp phân tử nước hấp phụ. Tuy nhiên, áp suất tách gần như không đổi trong khoảng độ ẩm tương đối ($RH$) từ 80% đến 100%.
  2. Sức căng bề mặt hạt keo (Bangham & Fakhoury, 1930; Powers, 1965, 1968; Wittmann, 1968, 1982): Chỉ phát huy tác dụng ở dải $RH$ rất thấp (5% – 50%), trong khi quá trình tự sấy khô do thủy hóa chỉ làm giảm $RH$ nội tại xuống khoảng 75% – 100%.
  3. Ứng suất căng mao dẫn (Capillary tension) (Hua et al., 1995; Lura et al., 2003; Coussy et al., 2004): Được thừa nhận rộng rãi là động lực chủ đạo tại vùng $RH > 45%$, liên kết với bán kính mặt khum chất lỏng qua phương trình Laplace và Kelvin.

Luận án của Đỗ Quang Huy định vị đột phá bằng cách khắc phục khoảng trống phương pháp luận của các nghiên cứu quốc tế:

  • So với mô hình vi cơ học đa tỷ lệ thuần giải tích của Bernard et al. (2003) và Ulm et al. (2004) vốn giả định pha C-S-H đồng nhất không đổi, luận án tích hợp trực tiếp sự biến thiên mật độ C-S-H theo thời gian trên nền tảng vi cấu trúc thực.
  • So với các mô hình phần tử hữu hạn vĩ mô của Pichler et al. (2007) hay mô hình bán thực nghiệm của Lura (2003), nghiên cứu này thực hiện bước nhảy vọt khi sử dụng cùng một vi cấu trúc vector 3D từ nền tảng $\mu$ic (Bishnoi & Scrivener, 2009) để giải quyết chuỗi liên hoàn: Thủy hóa $\rightarrow$ Độ xốp $\rightarrow$ Đàn hồi $\rightarrow$ Từ biến lão hóa $\rightarrow$ Co ngót tự sinh.

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

graph TD
    A[Mô phỏng Thủy hóa Vector 3D mic] --> B[Tiến hóa Vi cấu trúc & Độ xốp]
    B --> C[Phân tích Độ xốp & Mô phỏng MIP]
    B --> D[Đồng nhất hóa Đàn hồi FEM & SCS]
    D --> E[Loại bỏ Liên kết giả tạo - Double Burning Algorithm]
    B --> F[Mô hình Từ biến Viscoelastic Maxwell C-S-H]
    C --> G[Xác định Áp lực Mao dẫn Capillary Stress]
    F --> H[Phương pháp Chồng chất Từ biến Creep-Superposition]
    D --> H
    G --> H
    H --> I[Dự báo Co ngót Tự sinh Vĩ mô]

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

Luận án mang lại bốn đóng góp nền tảng cho lý thuyết hóa - cơ vật liệu xi măng:

  1. Mở rộng lý thuyết vi cơ học về pha C-S-H: Bác bỏ giả định truyền thống xem C-S-H là pha đồng nhất có khối lượng riêng cố định. Nghiên cứu chứng minh C-S-H hình thành ở tuổi sớm có độ đặc chắc thấp và trải qua quá trình nén chặt nội tại (densification) liên tục, tiến dần từ mật độ thấp đến khối lượng riêng thể tích $\rho = 2.0\text{ g/cm}^3$ và mật độ đóng gói $\eta = 0.7$.
  2. Làm rõ bản chất hình thái C-S-H trong không gian rỗng: Chứng minh sự sai lệch giữa mô phỏng số và thực nghiệm MIP bắt nguồn từ việc C-S-H phát triển dạng sợi kim (needles/diffuse growth) đâm xuyên vào không gian mao dẫn, tạo ra các lỗ rỗng siêu mịn bị giữ nước ngay từ giai đoạn đầu, thay vì tạo các mặt cầu nhẵn trơn như các giả định mô hình trước đây.
  3. Xác lập vai trò thống trị của từ biến trong co ngót tự sinh: Khẳng định co ngót tự sinh không đơn thuần là biến dạng xốp đàn hồi tức thời dưới áp suất mao dẫn mà là kết quả của sự tích lũy biến dạng từ biến lão hóa (aging creep) theo thời gian. Sự đóng góp của từ biến giải thích nguyên nhân tại sao hồ xi măng có tỷ lệ $w/c$ thấp lại chịu co ngót vượt bậc.
  4. Chuẩn hóa lý thuyết thẩm thấu cơ học (percolation threshold): Chứng minh rằng trạng thái đông kết cơ học của hồ xi măng được kiểm soát bởi sự kết hợp giữa hiệu ứng kết bông ban đầu của các hạt khoáng $C_3S$ và tốc độ gia tăng mật độ của gel C-S-H.

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

Khung phân tích tích hợp ba lý thuyết vật lý cốt lõi:

  • Lý thuyết Poro-mechanics và Ứng suất Mao dẫn (Young-Laplace & Kelvin):

$$p_c = -\frac{2\sigma\cos\theta}{r} = \frac{RT}{V_m}\ln(RH)$$

Trong đó $p_c$ là áp suất mao dẫn, $\sigma$ là sức căng bề mặt chất lỏng, $\theta$ là góc thấm ướt, $r$ là bán kính mặt khum Kelvin, $R$ là hằng số khí ($8.314\text{ J/mol}\cdot\text{K}$), $T$ là nhiệt độ tuyệt đối và $V_m$ là thể tích mol của nước lỗ rỗng.

  • Lý thuyết Đàn nhớt tuyến tính (Linear Viscoelasticity) với mô hình Maxwell tổng quát: Hàm từ biến nội tại của C-S-H được mô tả qua chuỗi Maxwell nhiều nhánh kết hợp phần tử lò xo đàn hồi và phần tử giảm chấn nhớt, tham số hóa từ dữ liệu thí nghiệm ấn nano (nano-indentation).
  • Lý thuyết Đồng nhất hóa Cơ học Vi mô (Micromechanical Homogenization): Kết hợp lược đồ tự nhất quán (Self-Consistent Scheme - SCS) và Phương pháp Phần tử Hữu hạn (FEM) trên thể tích đại diện (Representative Element Volume - REV) để chuyển đổi đặc tính cơ học từ cấp độ nano/micro lên cấp độ vĩ mô.

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

Thiết kế nghiên cứu

Nghiên cứu vận hành theo thế giới quan thực chứng (positivism) và chủ nghĩa hiện thực phản biện (critical realism), kết hợp mô phỏng số học chuyên sâu với dữ liệu thực nghiệm đối chứng đa phương thức.

Thiết kế nghiên cứu đa tỷ lệ bao gồm:

  • Cấp độ Nano: Ứng xử đàn nhớt nội tại của các lá C-S-H đơn lẻ ($10^{-9}\text{ m}$).
  • Cấp độ Micro: Tương tác pha khoáng clinker ($C_3S, C_2S, C_3A, C_4AF$), Portlandite ($CH$), C-S-H và hệ thống lỗ rỗng trong khối thể tích đại diện REV 3D ($10^{-6}\text{ m}$ đến $10^{-4}\text{ m}$).
  • Cấp độ Macro: Biến dạng co ngót thể tích và cơ tính vĩ mô của mẫu hồ xi măng ($10^{-2}\text{ m}$).
[Nano-scale: C-S-H Viscoelasticity (Nano-indentation)]
                        │
                        ▼
[Micro-scale: 3D Vector mic Platform + FEM Homogenization + Double Burning]
                        │
                        ▼
[Macro-scale: Creep-Superposition Shrinkage Formulation]

Quy trình nghiên cứu rigorous

  1. Khởi tạo vi cấu trúc 3D: Sử dụng nền tảng $\mu$ic với phân bố kích thước hạt (particle size distribution - psd) thực tế của xi măng thương mại, chứa tới hơn 3 triệu hạt trong thể tích tính toán, triệt tiêu hoàn toàn rào cản tính toán của các mô hình vector trước đây bằng thuật toán chia nhỏ lưới phụ và lấy mẫu điểm (point sampling).
  2. Kỹ thuật chia lưới và giải thuật "Double Burning": Khi chuyển đổi cấu trúc hình học vector sang lưới phần tử hữu hạn hexahedral, các hạt tiếp xúc gần tạo ra các "nút thắt liên kết giả tạo" (meshing artefacts), khiến khung chất rắn bị kết nối mạng (percolate) ngay tại thời điểm $t = 0$ (trước khi thủy hóa bắt đầu). Tác giả đã thiết kế thuật toán loại bỏ kết nối nhân tạo (double burning algorithm) dựa trên cơ sở dữ liệu vector hình học gốc, khôi phục chính xác thời điểm bắt đầu đông kết thực tế của vật liệu.
  3. Mô phỏng phép đo độ xốp: Tích hợp thuật toán đo kích thước lỗ rỗng và xâm nhập thủy ngân số (numerical MIP), so sánh đối chứng với phép đo thực nghiệm đo độ xốp xâm nhập thủy ngân và phổ cộng hưởng từ hạt nhân (Nuclear Magnetic Resonance - NMR).
  4. Đồng nhất hóa cơ học và tính toán từ biến: Lập trình giải thuật phần tử hữu hạn 3D đàn nhớt đẳng hướng với các biến trạng thái nội (internal variables), áp dụng điều kiện biên tuần hoàn và biến dạng thuần nhất để tính toán các tensor đàn hồi và tensor từ biến theo tuổi chất tải.

Data và phân tích

  • Mẫu vật liệu mô phỏng: Hệ hồ xi măng tricalcium silicate tinh khiết ($C_3S$) và xi măng Portland thương phẩm với các tỷ lệ $w/c$ biến thiên từ $0.30, 0.35, 0.40$ đến $0.475$.
  • Phần mềm và công cụ: Nền tảng vi cấu trúc $\mu$ic, bộ mã phần tử hữu hạn FE chuyên dụng phát triển trên ngôn ngữ C++/Fortran, công cụ tính toán đồng nhất hóa SCS.
  • Dữ liệu đối chứng thực nghiệm: Thí nghiệm đo biến dạng tuyến tính bằng ống lượn sóng (corrugated tube method theo tiêu chuẩn ASTM C1698), giàn đo biến dạng tự do (free deformation rig $100\times 100\times 500\text{ mm}$), số liệu ấn nano độ sâu xung lực và số liệu từ biến cơ bản tuổi sớm tại phòng thí nghiệm NAVIER (Université Paris-Est) và LMC (EPFL).

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

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

  1. Đột phá về cơ chế nén chặt C-S-H: Nếu giả định C-S-H có mật độ khối không đổi $\rho = 2.0\text{ g/cm}^3$ ($\eta = 0.7$), mô hình phần tử hữu hạn dự báo chính xác từ biến thực nghiệm ở các tuổi chất tải 24 giờ và 30 giờ, nhưng đánh giá quá thấp biến dạng từ biến ở tuổi chất tải sớm 18 giờ. Khi đưa hàm mật độ C-S-H tăng dần theo độ thủy hóa ($\rho(t)$ biến thiên), mô hình tái lập hoàn hảo ứng xử từ biến lão hóa tuổi cực sớm.
  2. Giải mã nghịch lý đường kính ngưỡng (break-through diameter) trong MIP: Đường kính cổ lỗ rỗng tới hạn thu được từ thuật toán số trên vi cấu trúc kỹ thuật số luôn lớn hơn đáng kể so với thực nghiệm MIP. Luận án chỉ ra nguyên nhân vật lý cốt lõi: C-S-H không phát triển lớp vỏ trơn nhẵn như các mô hình hình học cổ điển giả định, mà phát triển theo dạng tinh thể hình kim phân tán vào không gian rỗng, tạo ra mạng lưới vi mao dẫn nano giữ nước từ rất sớm.
  3. Ưu thế vượt trội của phương pháp Chồng chất Từ biến (Creep-Superposition): So sánh trực tiếp giữa hai hướng tiếp cận chỉ ra rằng mô hình xốp đàn hồi (poro-elasticity) đánh giá thấp nghiêm trọng độ co ngót tự sinh vĩ mô. Ngược lại, phương pháp chồng chất từ biến kết hợp mật độ C-S-H tiến hóa tái tạo chính xác độ dốc và biên độ co ngót ở các tỷ lệ $w/c$ khác nhau.
  4. Hiệu ứng hiệp đồng giữa kết bông hạt và tăng mật độ hydrate: Việc đưa hiệu ứng kết bông (flocculation) khi rải hạt ban đầu và sự nén chặt của C-S-H vào mô hình là điều kiện tiên quyết để đạt được sự trùng khớp tuyệt đối về ngưỡng percolation cơ học và động học phát triển mô đun đàn hồi $E(t)$.
So sánh biến dạng co ngót theo thời gian:
Biến dạng Co ngót (µm/m)
 │                                        [Creep-Superposition + Densification] ── Thực nghiệm
 │                                      ╱
 │                                    ╱
 │                                  ╱  [Poro-elasticity thuần túy] (Dự báo thấp)
 │                        ─────────
 └────────────────────────────────────────────► Thời gian thủy hóa (giờ)

Implications đa chiều

  • Về mặt lý thuyết: Thống nhất các nhánh nghiên cứu riêng rẽ về động học hóa học, vi cơ học đàn nhớt và nhiệt động học mao dẫn thành một hệ thống tính toán khép kín.
  • Về phương pháp luận: Cung cấp thuật toán Double Burning giải quyết triệt để lỗi chia lưới trong vi cơ học tính toán vật liệu composite hạt, có thể chuyển giao cho nghiên cứu gốm sứ, địa chất và vật liệu polyme xốp.
  • Về ứng dụng thực tiễn: Cho phép các kỹ sư vật liệu dự báo chính xác nguy cơ nứt sớm của các trạm trộn bê tông HPC/UHPC dựa trên thành phần khoáng và độ mịn xi măng mà không cần chờ đợi các quy trình đo thực nghiệm kéo dài nhiều tuần.

Limitations và Future Research

Luận án thẳng thắn chỉ ra bốn giới hạn phương pháp luận chính:

  1. Rào cản quy mô mô phỏng cấu trúc rỗng: Do hình thái phát triển dạng sợi kim phân tán cực kỳ phức tạp của C-S-H, việc tái tạo trực tiếp toàn bộ chi tiết không gian rỗng ở quy mô hàng trăm micromet bằng mô hình số thuần túy vẫn chưa thể thay thế hoàn toàn dữ liệu đầu vào thực nghiệm về áp suất mao dẫn.
  2. Giả định đẳng nhiệt và đẳng ẩm lý tưởng: Mô hình tập trung vào điều kiện kín đẳng nhiệt 20°C, chưa ghép nối đầy đủ trường nhiệt phi ổn định sinh ra từ phản ứng tỏa nhiệt thủy hóa khối lớn và ảnh hưởng của phụ gia giảm co ngót (SRA) hay tro bay/xỉ lò cao.
  3. Giới hạn mô tả cơ chế trương nở ban đầu: Sự xuất hiện của áp suất kết tinh (crystallization pressure) do tạo mầm ettringite hoặc portlandite gây trương nở nhẹ ở giai đoạn đầu ninh kết chưa được mô hình hóa tự động hoàn toàn mà mới chỉ được xem xét như một thành phần bù trừ biến dạng.

Các hướng nghiên cứu tiếp nối trong chương trình 10 năm:

  • Phát triển mô hình tăng trưởng C-S-H phân nhánh fractal để mô phỏng chính xác cấu trúc vi mao dẫn mà không cần hiệu chỉnh bán thực nghiệm.
  • Ghép nối mô hình nhiệt - ẩm - cơ (Thermo-Hygro-Mechanical - THM) đa chiều cho kết cấu bê tông khối lớn.
  • Mở rộng thư viện vật liệu $\mu$ic cho các hệ chất kết dính xanh thế hệ mới (Alkali-activated materials, LC3 - Limestone Calcined Clay Cement).

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

  • Tác động học thuật: Luận án là đỉnh cao nghiên cứu suốt gần 20 năm tại phòng thí nghiệm LMC/EPFL về mô hình hóa vi cấu trúc xi măng, đóng vai trò nền tảng cho hàng trăm công trình trích dẫn quốc tế trong lĩnh vực cơ học vật liệu xây dựng tính toán.
  • Chuyển đổi công nghiệp: Cung cấp công cụ mô phỏng vi cấu trúc phục vụ các tập đoàn vật liệu xây dựng toàn cầu (như LafargeHolcim, Sika) tối ưu hóa cấp phối bê tông cường độ siêu cao (UHPC), cắt giảm hàm lượng clinker phát thải $CO_2$ mà vẫn bảo đảm độ bền nứt.
  • Lợi ích kinh tế - xã hội: Giảm thiểu chi phí bảo trì và rủi ro sụp đổ cầu cảng, công trình ngầm và nhà cao tầng do nứt nẻ tuổi sớm gây ăn mòn cốt thép.

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

  • Nghiên cứu sinh & Nhà khoa học trẻ: Tiếp cận mã nguồn phương pháp luận đồng nhất hóa cơ học số và giải thuật xử lý lưới vector-to-FEM không lỗi.
  • Giáo sư & Chuyên gia vi cơ học: Thừa hưởng khung lý thuyết đã chứng minh về từ biến lão hóa C-S-H và phương pháp chồng chất từ biến trong vật liệu xốp.
  • Phòng R&D Doanh nghiệp Xây dựng: Ứng dụng các thông số tối ưu về tỷ lệ $w/c$, độ mịn và thành phần khoáng để thiết kế hỗn hợp bê tông chống nứt tin cậy.
  • Cơ quan quản lý tiêu chuẩn kỹ thuật: Bổ sung cơ sở khoa học để hoàn thiện các tiêu chuẩn dự báo co ngót và nứt kết cấu bê tông hiện đại.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp độc đáo nhất là việc tích hợp thành công lý thuyết nén chặt mật độ C-S-H theo thời gian vào phương pháp chồng chất từ biến (creep-superposition approach) trên nền tảng vi cơ học. Công trình đã mở rộng lý thuyết xốp đàn nhớt (poro-viscoelasticity) của Biot và Coussy, chứng minh rằng biến dạng co ngót tự sinh ở tuổi sớm bị chi phối áp đảo bởi hiện tượng từ biến lão hóa của gel C-S-H dưới tác dụng của ứng suất căng mao dẫn biến thiên.

2. Đột phá phương pháp luận của nghiên cứu khi so sánh với các công trình quốc tế tiền nhiệm?

So với mô hình rời rạc CemHyd3D (Bentz, 1997) bị giới hạn độ phân giải $1\text{ }\mu\text{m}$ và mô hình HymoStruc (van Breugel, 1995) chỉ xét một pha hydrate duy nhất, luận án đã:

  • Ứng dụng nền tảng vector 3D $\mu$ic xử lý hàng triệu hạt clinker đa khoáng với động học mầm hóa chân thực.
  • Sáng tạo thuật toán Double Burning Algorithm triệt tiêu hoàn toàn các liên kết giả tạo của lưới phần tử hữu hạn, loại bỏ hiện tượng percolation rắn nhân tạo trước ninh kết – một vấn đề nan giải tồn tại suốt nhiều thập kỷ trong mô hình hóa vi cơ học xi măng.

3. Phát hiện thực nghiệm/mô phỏng nào gây bất ngờ nhất và bằng chứng dữ liệu đi kèm?

Phát hiện bất ngờ nhất là đường kính xuyên phá (break-through diameter) tính toán từ mô phỏng MIP trên cấu trúc số hóa luôn lớn hơn nhiều so với thực nghiệm đo MIP thực tế. Dữ liệu chứng minh rằng C-S-H không tạo thành các mặt cầu nhẵn bao bọc hạt mà mọc dạng kim/sợi phân tán (diffuse/needle-like growth) lấp đầy không gian rỗng, tạo ra mạng vi mao dẫn nano giữ nước ngay từ 18–24 giờ đầu.

Mô hình giả định cũ:                 Thực tế chỉ ra trong luận án:
   [ Hạt Clinker ]                     [ Hạt Clinker ]
    ╭───────────╮                       ╭───*─*──*──╮
    │  C-S-H    │ (Vỏ nhẵn trơn)        │ ╲╱ ╲╱ ╲╱  │ (C-S-H dạng sợi kim/diffuse)
    │           │                       │  *  *  *  │ -> Giữ nước vi mao dẫn
    ╰───────────╯                       ╰───*─*──*──╯

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) rõ ràng không?

Có. Luận án chi tiết hóa toàn bộ hệ phương trình điều khiển đàn nhớt 3D (Phụ lục C, E), thuật toán chia lưới vector, bộ thông số đàn hồi của từng pha clinker ($C_3S, C_2S, C_3A, C_4AF, CH, C\text{-}S\text{-}H$) và quy trình khớp mô hình Maxwell tổng quát từ dữ liệu ấn nano, cho phép cộng đồng khoa học lập trình tái lập hoàn toàn kết quả mô phỏng.

5. Lộ trình nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?

Lộ trình tập trung vào việc:

  • Mô hình hóa trực tiếp hình thái phân tán sợi kim nano của C-S-H để lượng hóa áp suất mao dẫn mà không cần dữ liệu phụ trợ thực nghiệm.
  • Mở rộng sang bài toán ghép nối đa trường Thermo-Hygro-Chemo-Mechanical (THCM) cho kết cấu bê tông đại khối trong điều kiện môi trường thực tế.
  • Tích hợp mô phỏng vi cấu trúc cho các hệ chất kết dính bền vững chứa hàm lượng phụ gia khoáng hoạt tính cao.

Kết luận

Luận án tiến sĩ của Đỗ Quang Huy đã xác lập một mốc son học thuật quan trọng trong lĩnh vực cơ học vật liệu xi măng tính toán thông qua 5 đóng góp cốt lõi:

  1. Xây dựng khung mô hình vi cơ học 3D hoàn chỉnh: Liên kết liền mạch từ động học thủy hóa hạt vector đến đặc tính cơ học đàn nhớt và biến dạng thể tích vĩ mo.
  2. Phát minh thuật toán Double Burning: Loại bỏ triệt để lỗi liên kết giả tạo của phương pháp phần tử hữu hạn, xác định chính xác thời điểm đông kết cơ học.
  3. Chứng minh quy luật nén chặt của C-S-H: Khẳng định mật độ khối và mật độ đóng gói của C-S-H biến thiên theo tuổi thủy hóa, giải thích chính xác từ biến lão hóa từ 18 giờ đến 30 giờ.
  4. Làm rõ bản chất cơ chế co ngót tự sinh: Khẳng định phương pháp chồng chất từ biến (creep-superposition) là chìa khóa duy nhất mô tả chính xác ảnh hưởng của tỷ lệ $w/c$ thấp lên co ngót tự sinh.
  5. Định danh hình thái sợi phân tán của hydrate: Mở đường cho các lý thuyết mới về truyền dẫn ion và độ bền mao dẫn trong vật liệu gốc xi măng.

Công trình không chỉ giải quyết các câu hỏi cơ bản về vật lý - hóa học vật liệu mà còn cung cấp nền tảng khoa học vững chắc giúp ngành công nghiệp xây dựng toàn cầu kiểm soát hiện tượng nứt sớm, nâng cao tuổi thọ và độ an toàn cho các công trình hạ tầng thế kỷ.