Tổng quan về luận án

Công nghệ bê tông đầm lăn (Roller Compacted Concrete - BTĐL/RCC) là giải pháp đột phá trong xây dựng các công trình thủy công quy mô lớn nhờ tốc độ thi công nhanh, giảm thiểu chi phí kết cấu phụ trợ và tối ưu hóa hiệu quả kinh tế. Tuy nhiên, bản chất vật liệu bê tông khối lớn với hàm lượng xi măng thấp (60–100 kg/m³) kết hợp phụ gia khoáng hoạt tính (PGK) vẫn đối mặt với hiểm họa nứt do ứng suất nhiệt phát sinh từ quá trình thủy hóa tỏa nhiệt chậm nhưng kéo dài. Tại Việt Nam, nhiều công trình trọng điểm quy mô lớn như đập thủy điện Sơn La (cao 138 m), Pleikrông (71 m), Bản Vẽ (138 m) đã ghi nhận hiện tượng nứt bề mặt và nứt xuyên do ứng suất nhiệt vượt quá cường độ kháng kéo của bê tông dưới điều kiện khí hậu nhiệt đới gió mùa khắc nghiệt.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ việc các dự án tại Việt Nam chủ yếu áp dụng quy chuẩn, tiêu chuẩn từ Hoa Kỳ (USACE, USBR, ACI), Trung Quốc (SL314-2004, SL319-2005) hoặc Nhật Bản (RCD) vốn được xây dựng trên điều kiện biên khí hậu ôn đới hoặc cận nhiệt đới và nguồn vật liệu bản địa khác biệt. Chưa có nghiên cứu hệ thống nào tích hợp đồng thời các yếu tố thành phần khoáng clinker xi măng sản xuất trong nước ($C_3A, C_3S$), đặc tính phụ gia khoáng tro bay/puzơlan, nhiệt độ hỗn hợp đổ ban đầu ($T_p$) và dao động vi khí hậu ba miền để đưa ra giải pháp giảm ứng suất nhiệt phân vùng tối ưu.

Luận án tiến sĩ kỹ thuật của tác giả Nguyễn Minh Việt (2017) tại Viện Khoa học Thủy lợi Việt Nam, dưới sự hướng dẫn của PGS. Hoàng Phó Uyên và GS. Phạm Ngọc Khánh, đã giải quyết căn cơ bài toán này với hệ thống câu hỏi và giả thuyết nghiên cứu rõ ràng:

  1. Câu hỏi nghiên cứu 1 (RQ1): Quy luật phát triển trường nhiệt độ và ứng suất nhiệt trong thân đập BTĐL thay đổi như thế nào dưới tác động tương hỗ giữa nhiệt thủy hóa của xi măng nội địa và điều kiện khí tượng từng vùng lãnh thổ?
  2. Câu hỏi nghiên cứu 2 (RQ2): Tỷ lệ phối trộn phụ gia khoáng (PGK/CKD) và giới hạn hàm lượng khoáng tỏa nhiệt ($C_3A + C_3S$) bao nhiêu là tối ưu để triệt tiêu nguy cơ nứt do ứng suất nhiệt mà vẫn bảo đảm độ bền cơ lý?
  3. Giả thuyết nghiên cứu 1 (H1): Khống chế tổng hàm lượng khoáng $(C_3A + C_3S) \le 58%$ trong xi măng kết hợp điều chỉnh tỷ lệ PGK từ 35–50% sẽ làm giảm nhiệt độ đỉnh khối đổ ($\Delta T$) từ $25^\circ\text{C}$ xuống dưới $17^\circ\text{C}$, qua đó triệt tiêu ứng suất kéo nhiệt nguy hiểm tại các biên thượng/hạ lưu.
  4. Giả thuyết nghiên cứu 2 (H2): Một thuật toán mô phỏng phần tử hữu hạn (PTHH) sử dụng ngôn ngữ APDL trong ANSYS kết hợp kỹ thuật phần tử “sinh - tử” có thể dự báo chính xác trường nhiệt - ứng suất nhiệt theo từng lớp đổ (lift 0,3 m) với sai số quan trắc dưới 5%.

Khung lý thuyết của luận án tích hợp phương trình vi phân dẫn nhiệt Fourier phi dừng, định luật khuếch tán nhiệt độ không gian ba chiều, lý thuyết đàn hồi nhiệt Hooke mở rộng có xét đến từ biến (creep) và co ngót (shrinkage). Phạm vi nghiên cứu bao phủ dữ liệu thực nghiệm và quan trắc tại hơn 20 đập BTĐL lớn tại Việt Nam, kiểm chứng thực tế tại thủy điện Sơn La và chuyển giao tính toán ứng dụng cho đập thủy điện Trung Sơn (cao 84 m, Thanh Hóa).

Literature Review và Positioning

Lịch sử phát triển công nghệ BTĐL khởi nguồn từ những năm 1960–1970 với các công trình tiên phong như đập Alpe Gera (Italia, 1961), Manicouagan (Canada, 1961) và báo cáo kinh điển "Đập trọng lực tối ưu" của Giáo sư Jerome M. Raphael (1970). Nghiên cứu của Cannon R.W. (1972–1974) cùng Công binh Lục quân Hoa Kỳ (USACE) tại đập Lost Creek (1977) và Willow Creek (1980, cao 52 m, 331.000 m³ BTĐL) đã xác lập trường phái RCC nghèo xi măng. Tại Anh, Dunstan (1970) và Hiệp hội Nghiên cứu Thông tin Công nghiệp Xây dựng (CIRIA, 1976–1979) phát triển trường phái BTĐL giàu tro bay, tiền đề cho Cục Khai hoang Hoa Kỳ (USBR) xây dựng đập Upper Stillwater (1987, cao 90 m, 1,126 triệu m³ BTĐL). Đồng thời, Nhật Bản thiết lập trường phái RCD (Roller Compacted Dams) với kết cấu “vàng bọc bạc” (lõi BTĐL, vỏ bê tông thường chống thấm) tại đập Shimajigawa (1980, cao 89 m), trong khi Trung Quốc hoàn thiện hệ tiêu chuẩn RCCD qua hơn 165 công trình (đập Khanh Khẩu 1986, Tam Hiệp, Long Than).

Trong bức tranh học thuật quốc tế, Barrett và Tatro (1992) đặt nền móng giải pháp phân tích nhiệt - ứng suất nhiệt bằng phương pháp PTHH có xét đến biến dạng từ biến và co ngót. Về mặt nhiệt học clinker, các công trình của Woods (1932), McCarthy (1984) và Verbeck (1968) đã lượng hóa nhiệt thủy hóa của các khoáng đơn: $C_3A$ (1023 J/g), $C_3S$ (520 J/g), $C_4AF$ (337 J/g), $C_2S$ (247 J/g).

Tại Việt Nam, các tác giả Nguyễn Như Quý (2000), Nguyễn Tiến Đích (2001), Đỗ Hồng Hải (2006) và Lê Quang Hùng (2007) đã phân tích khả năng giảm nhiệt cách nhiệt của tro bay Phả Lại và puzơlan Long Phước. Lê Anh Vân, Nguyễn Văn Mạo (2006) khảo sát phát triển nhiệt tại đập Pleikrông; Đinh Hữu Dụng (2009) phân tích nhiệt tải thi công đập Sơn La; Đỗ Văn Lượng (2010) nghiên cứu ứng suất nhiệt đập bê tông khối lớn; Lê Quốc Toàn (2014) khảo sát tốc độ lên đập Đồng Nai 4 bằng ANSYS nhưng chỉ dừng ở các thông số cơ lý đơn lẻ, thiếu tính tổng thể theo vùng địa lý tự nhiên.

Tồn tại một tranh luận học thuật lớn giữa hai trường phái:

  • Trường phái 1 (USACE, Trung Quốc): Ưu tiên các biện pháp cưỡng bức cơ học như làm lạnh cốt liệu, dùng đá vảy, tưới nước mặt khối đổ và phân khe co giãn dày đặc.
  • Trường phái 2 (Dunstan, CIRIA, luận án này): Chủ động điều chỉnh từ gốc phản ứng nhiệt hóa sinh bằng cách tối ưu hóa thành phần khoáng clinker xi măng và tăng tối đa hàm lượng PGK hoạt tính cao, kết hợp kiểm soát vi khí hậu khu vực.

Vị thế học thuật (positioning) của luận án được xác lập rõ nét: Đây là công trình đầu tiên tại Việt Nam giải mã bản chất tương tác giữa thành phần khoáng xi măng nội địa ($C_3A, C_3S$) và tỷ lệ phối trộn PGK trong môi trường khí hậu nhiệt đới gió mùa phân vùng (Miền núi phía Bắc, Bắc Trung Bộ, Nam Trung Bộ & Tây Nguyên), khắc phục triệt để hạn chế của các nghiên cứu đơn lẻ trước đó.

Đó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 nhiệt thủy hóa bê tông khối lớn thông qua việc thiết lập hàm toán học liên tục mô tả động học tỏa nhiệt của chất kết dính hỗn hợp (Xi măng + PGK) có điều chỉnh theo tỷ lệ thành phần khoáng clinker Việt Nam: $$Q = a \cdot C_3S(%) + b \cdot C_2S(%) + c \cdot C_3A(%) + d \cdot C_4AF(%) \quad (\text{J/g})$$ với hệ số thực nghiệm hiệu chỉnh trung bình: $a = 520\text{ J/g}$, $b = 247\text{ J/g}$, $c = 1023\text{ J/g}$, $d = 337\text{ J/g}$.

Đồng thời, luận án đã xác lập mô hình giải tích thực nghiệm dự báo sự phát triển cường độ kháng nén ($f_c$), mô đun đàn hồi ($E_c$) và độ dãn dài giới hạn ($\varepsilon_p$) của BTĐL theo thời gian và hàm lượng chất kết dính: $$f_{c28} = 19,459(CKD/N) - 0,147(PGK/CKD) - 11,681 \quad (\text{MPa})$$ $$f_{c90} = 19,326(CKD/N) - 0,333(PGK/CKD) + 5,968 \quad (\text{MPa})$$ $$E_{c90} = 21,217(CKD/N) - 0,197(PGK/CKD) - 1,376 \quad (\text{GPa})$$

Mô hình lý thuyết này chứng minh rằng việc tăng tỷ lệ PGK/CKD trực tiếp làm giảm mô đun đàn hồi $E_c$ ở giai đoạn đầu, tạo điều kiện thuận lợi cho sự xả ứng suất thông qua từ biến dẻo, qua đó làm giảm đáng kể nguy cơ nứt do nền kiềm chế.

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 khối lý thuyết trụ cột:

  1. Lý thuyết truyền nhiệt phi tuyến (Fourier - Laplace): Giải phương trình vi phân dẫn nhiệt không dừng kết hợp tỏa nhiệt nội sinh $q_v(t)$, trao đổi nhiệt đối lưu bề mặt $\alpha(T_w - T_f)$ và hấp thụ bức xạ mặt trời $T_{ow} = \frac{Q_o}{\alpha_w}$.
  2. Lý thuyết cơ học đàn hồi - từ biến vật liệu composite hạt: Tính toán tương tác giữa pha nền xi măng thủy hóa chậm và cốt liệu chịu kiềm chế ngoại vi.
  3. Lý thuyết kiểm soát chênh lệch nhiệt độ biên quy chuẩn: Định lượng chênh lệch nhiệt độ đáy móng ($T = T_p + T_r - T_f$), chênh lệch lớp trên - dưới ($T_M \le 15 \div 20^\circ\text{C} + T_m - \varepsilon T_b$) và chênh lệch trong - ngoài ($T_m = \Delta T_0 + \Delta T_m + T_S$).

Điều kiện biên xác định rõ: Khung phân tích áp dụng cho các đập bê tông đầm lăn trọng lực có chiều cao từ 50 m đến trên 140 m đặt trên nền đá cứng có mô đun biến dạng tương thích, thi công theo công nghệ rải lớp đầm lăn liên tục trong điều kiện nhiệt đới gió mùa.

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

Thiết kế nghiên cứu

Nghiên cứu được định vị theo bản thể luận khách quan (positivism) và nhận thức luận diễn dịch thực nghiệm (deductive quantitative paradigm). Thiết kế nghiên cứu kết hợp chặt chẽ giữa mô hình toán số trị giải tích phần tử hữu hạn đa trường (Coupled Transient Thermo-Mechanical Analysis) và đối soát thực chứng quan trắc thực địa.

Quy mô dữ liệu bao gồm:

  • Toàn bộ chuỗi số liệu quan trắc nhiệt độ của đập thủy điện Sơn La (trạm đo nhiệt độ tại cao trình +114,0 m với các cảm biến nhiệt điện trở chôn sẵn trong thân đập).
  • Tập mẫu thí nghiệm thành phần khoáng clinker từ 12 nhà máy xi măng lớn tại Việt Nam (Phụ lục I: Hoàng Thạch, Bút Sơn, Bỉm Sơn, Nghi Sơn, Chinfon, Holcim...).
  • Dữ liệu khí tượng thủy văn hơn 30 năm từ Trung tâm Dự báo Khí tượng Thủy văn Quốc gia đại diện cho 3 phân vùng sinh khí hậu.

Quy trình nghiên cứu rigorous

Quy trình mô phỏng số được thiết lập tự động hóa bằng ngôn ngữ APDL trong ANSYS qua hai giai đoạn liên hoàn:

  1. Phân tích truyền nhiệt phi định thường (Thermal Analysis): Sử dụng phần tử phẳng 4 nút PLANE55 (có bậc tự do nhiệt độ tại mỗi nút) để mô phỏng quá trình truyền nhiệt dẫn nhiệt nội tại, đối lưu qua mặt thoáng và bức xạ mặt trời. Kỹ thuật "Element Birth and Death" được lập trình để kích hoạt từng lớp đổ bê tông dày 0,30 m theo đúng tiến độ rải và đầm lèn thực tế.
  2. Phân tích trường ứng suất đàn dẻo kết hợp từ biến (Structural Analysis): Chuyển đổi toàn bộ phần tử PLANE55 sang phần tử cơ học tương thích PLANE182 (2 bậc tự do chuyển vị $u_x, u_y$ tại mỗi nút). Trường nhiệt độ tính được từ bước 1 được gán làm tải trọng thể tích nút (body load) cho bước 2 để tính toán trường ứng suất chính $\sigma_1, \sigma_3$ và ứng suất cắt $\tau_{xy}$.

Tính tam giác hóa (triangulation) được bảo đảm nghiêm ngặt: Kết quả mô hình toán được đối chiếu đồng thời với (1) Lời giải giải tích cổ điển Sitmit sai phân, (2) Dữ liệu đo đạc thực tế tại đập Sơn La và (3) Kết quả thí nghiệm kiểm chứng tại đập Trung Sơn. Độ tin cậy của thuật toán đạt độ chính xác cao với hệ số tương quan $R^2 > 0,96$ giữa đường cong nhiệt mô phỏng và thực đo.

       QUY TRÌNH MÔ PHỎNG SỐ TRONG ANSYS (APDL)
       
+-------------------------------------------------------+
|  BƯỚC 1: MÔ HÌNH NHIỆT (Thermal Module: PLANE55)       |
|  - Gán nhiệt thủy hóa Q(t) theo khoáng C3A, C3S       |
|  - Gán biên đối lưu alpha(t), bức xạ Tow, T_kk(t)     |
|  - Kích hoạt lớp đổ 0,3m bằng Element Birth/Death    |
|  - Giải trường nhiệt độ phi tuyến: T(x, y, z, t)      |
+---------------------------+---------------------------+
                            |
                     Chuyển đổi phần tử
                 (ETCHG, TTS: PLANE55 -> PLANE182)
                            |
                            v
+-------------------------------------------------------+
|  BƯỚC 2: MÔ HÌNH CƠ HỌC (Structural: PLANE182)        |
|  - Gán T(x, y, z, t) làm tải trọng thể tích Body Load |
|  - Gán điều kiện kiềm chế đàn hồi của nền đá          |
|  - Cập nhật hàm phát triển Ec(t), Creep, fc(t), ft(t) |
|  - Tính toán trường ứng suất kéo chính Sigma_1(x, y)  |
+---------------------------+---------------------------+
                            |
                            v
+-------------------------------------------------------+
|  BƯỚC 3: ĐÁNH GIÁ NGUY CƠ NỨT & TỐI ƯU PHÂN VÙNG      |
|  - Kiểm tra tiêu chí an toàn: Sigma_1 <= [ft] / K_at  |
|  - Xuất biểu đồ tương quan PGK, CX, Tp theo 3 vùng   |
+-------------------------------------------------------+

Data và phân tích

Mô hình PTHH làm việc đồng thời giữa thân đập và nền đá: Nền đá được mô phỏng kéo dài về phía thượng lưu, hạ lưu và đáy sâu một khoảng bằng $1,5 \div 2,0$ lần chiều cao đập nhằm triệt tiêu ảnh hưởng biên giả tạo.

  • Thông số nhiệt - cơ nền đá: Mô đun đàn hồi $E_r = 15 \div 25\text{ GPa}$, hệ số Poisson $\nu = 0,22$, hệ số dẫn nhiệt $\lambda = 2,0\text{ kcal/m}\cdot\text{h}\cdot^\circ\text{C}$.
  • Thông số BTĐL: Cấp độ bền nén tuổi 90 ngày đạt 15–20 MPa, cường độ kháng kéo $f_t = 1,2 \div 1,8\text{ MPa}$, hệ số giãn nở nhiệt $\alpha_T = 1,0 \times 10^{-5}\text{ /}^\circ\text{C}$, khối lượng thể tích $\gamma = 2400\text{ kg/m}^3$.
  • Kiểm tra độ vững chắc (Robustness checks): Mô hình được thử nghiệm chạy qua hàng loạt kịch bản thay đổi bước thời gian giải tích ($\Delta t = 1\text{ h}, 6\text{ h}, 12\text{ h}, 24\text{ h}$) và kích thước lưới chia phần tử hữu hạn ($0,15\text{ m} \div 0,60\text{ m}$), xác nhận kết quả hội tụ ổn định tuyệt đối.

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

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

Luận án rút ra các phát hiện khoa học mang tính quy luật có ý nghĩa quyết định đối với an toàn kết cấu đập BTĐL:

  1. Quy luật chi phối của tổ hợp khoáng $(C_3A + C_3S)$ đối với đỉnh nhiệt và ứng suất nhiệt: Hàm lượng khoáng $(C_3A + C_3S)$ trong xi măng tỷ lệ thuận tuyến tính với nhiệt độ cực đại và ứng suất kéo chính cực đại tại biên thượng/hạ lưu. Khi tổng hàm lượng $(C_3A + C_3S)$ tăng từ 50% lên 66%, nhiệt độ đỉnh trong thân đập tăng từ $3,2^\circ\text{C}$ đến $4,8^\circ\text{C}$, kéo theo ứng suất kéo chính $\sigma_1$ tại mép biên tăng thêm $0,35 \div 0,52\text{ MPa}$, trực tiếp vượt ngưỡng kháng kéo tức thời của bê tông non tuổi và gây nứt bề mặt.

  2. Hiệu ứng xả ứng suất của Phụ gia khoáng (PGK): Sử dụng PGK (tro bay hoặc puzơlan) thay thế xi măng với tỷ lệ từ 35% đến 50% tổng lượng chất kết dính có tác dụng kép: Giảm lượng nhiệt sinh ra trong 72 giờ đầu từ $30 \div 45%$ và làm chậm tốc độ tăng mô đun đàn hồi $E_c(t)$, giúp khối bê tông hấp thụ biến dạng thông qua từ biến mà không tích tụ ứng suất kéo nguy hiểm.

  3. Tác động phi đối xứng của điều kiện nhiệt ẩm môi trường theo vùng miền:

  • Tại Miền núi phía Bắc (biên độ nhiệt năm $T_b$ lớn từ $10 \div 14^\circ\text{C}$, mùa đông lạnh khô): Ứng suất kéo nguy hiểm nhất xuất hiện ở bề mặt thượng - hạ lưu vào mùa đông do chênh lệch nhiệt độ trong - ngoài lớn ($T_m > 22^\circ\text{C}$).
  • Tại Bắc Trung Bộ (chịu ảnh hưởng gió phơn Tây Nam khô nóng): Hiện tượng thoát ẩm bề mặt cực nhanh làm gia tăng ứng suất co ngót nhiệt kết hợp co khô.
  • Tại Nam Trung Bộ & Tây Nguyên (nhiệt độ nền cao quanh năm): Nhiệt độ hỗn hợp bê tông ban đầu ($T_p$) là nhân tố nhạy cảm nhất quyết định đỉnh nhiệt thủy hóa.
  1. Ngưỡng nhiệt độ đổ khống chế ($T_p$): Mỗi mức giảm $1^\circ\text{C}$ của nhiệt độ bê tông khi đổ ($T_p$) giúp giảm đỉnh nhiệt thân đập xấp xỉ $0,85 \div 0,92^\circ\text{C}$ và giảm ứng suất kéo bề mặt $0,08\text{ MPa}$. Giới hạn an toàn của $T_p$ cho các vùng tại Việt Nam được xác định dao động từ $20^\circ\text{C} \div 24^\circ\text{C}$.
Thông số Phân vùng / Giải pháp Miền núi phía Bắc Bắc Trung Bộ Nam Trung Bộ & Tây Nguyên
Đặc thù khí hậu Mùa đông rất lạnh, $T_b = 10 \div 14^\circ\text{C}$ Nắng nóng gay gắt, gió phơn khô nóng Nhiệt độ cao ổn định quanh năm
Hàm lượng khoáng $(C_3A + C_3S)$ $\le 52 \div 54%$ $\le 54 \div 56%$ $\le 56 \div 58%$
Tỷ lệ phối trộn PGK/CKD $35 \div 45%$ (Tro bay Phả Lại/Vũng Áng) $40 \div 50%$ (Tro bay Nghi Sơn/Vũng Áng) $35 \div 45%$ (Puzơlan Long Phước/Gia Lai)
Nhiệt độ đổ bê tông khống chế ($T_p$) $\le 21 \div 23^\circ\text{C}$ $\le 20 \div 22^\circ\text{C}$ $\le 22 \div 24^\circ\text{C}$
Biện pháp dưỡng hộ bề mặt bổ sung Che chắn nhiệt biên mùa đông, chống sốc nhiệt Phun sương ẩm liên tục chống co ngót khô Che phủ bạt cách nhiệt, làm mát cốt liệu

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp phương pháp luận định lượng liên kết giữa hóa học tinh thể clinker xi măng với cơ học môi trường liên tục trong công trình thủy công khối lớn.
  • Về mặt phương pháp: Xây dựng bộ module macro APDL hoàn chỉnh cho phần mềm ANSYS, cho phép các kỹ sư thiết kế mô phỏng tức thì bài toán nhiệt - ứng suất nhiệt cho bất kỳ hình học mặt cắt đập nào tại Việt Nam.
  • Về mặt thực tiễn và chính sách: Đề xuất sửa đổi, bổ sung vào các Tiêu chuẩn Quốc gia (TCVN) và Quy chuẩn Kỹ thuật Quốc gia (QCVN) về thiết kế, thi công đập BTĐL các điều khoản bắt buộc về giới hạn thành phần khoáng xi măng và kiểm soát nhiệt phân vùng sinh khí hậu.

Limitations và Future Research

Luận án thừa nhận một số giới hạn nghiên cứu khách quan:

  1. Mô hình hình học: Phân tích ứng suất nhiệt dừng lại ở bài toán biến dạng phẳng 2D (mặt cắt đại diện đập - nền), chưa phản ánh toàn diện hiệu ứng tập trung ứng suất 3D cục bộ tại các khu vực tiếp giáp vai đập có địa hình thung lũng dốc hoặc xung quanh cụm công trình ngầm (hành lang tiêu nước, ống luồn tuabin).
  2. Đặc trưng phi tuyến vật liệu: Mô hình đàn dẻo nhiệt xem nền đá là môi trường liên tục đồng nhất, chưa xét đến tính dị hướng của các hệ thống khe nứt, đứt gãy kiến tạo phức tạp trong khối đá nền.
  3. Đặc trưng độ ẩm ghép đôi: Ảnh hưởng của độ ẩm môi trường chủ yếu được tích hợp thông qua hệ số trao đổi nhiệt đối lưu bề mặt $\alpha$, chưa giải bài toán truyền ẩm phi tuyến liên kết hoàn toàn (fully-coupled hygro-thermo-mechanical analysis).

Chương trình nghiên cứu tiếp theo (Future Research Agenda) cần tập trung vào 4 hướng chiến lược:

  • Hướng 1: Phát triển mô hình mô phỏng số 3D toàn khối đập kết hợp tương tác động đất và áp lực nước thấm trong điều kiện hồ chứa tích nước.
  • Hướng 2: Ứng dụng vật liệu nano-silica và tro bay siêu mịn nghiền cơ hóa để nâng cao độ đặc chắc và khả năng kháng nứt sớm của BTĐL.
  • Hướng 3: Tích hợp trí tuệ nhân tạo (AI/Machine Learning) và mạng nơ-ron vật lý (Physics-Informed Neural Networks - PINNs) để hiệu chỉnh thời gian thực trường nhiệt độ từ dữ liệu cảm biến IoT tại hiện trường.
  • Hướng 4: Đánh giá suy thoái độ bền lâu dài của đập BTĐL dưới tác động của biến đổi khí hậu cực đoan và xâm thực hóa học.

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

Nghiên cứu mang lại giá trị gia tăng to lớn trên nhiều phương diện:

  • Giá trị học thuật: Đặt nền tảng cho trường phái tính toán nhiệt - cơ đập BTĐL nhiệt đới tại Việt Nam, đóng góp nguồn tài liệu tham khảo cốt lõi cho các chương trình đào tạo sau đại học chuyên ngành Xây dựng Công trình Thủy và Địa kỹ thuật Công trình.
  • Chuyển đổi ngành công nghiệp xây dựng thủy điện: Giúp các nhà thầu lớn (Tập đoàn Điện lực Việt Nam - EVN, Tổng Công ty Sông Đà, LILAMA...) tối ưu hóa dây chuyền sản xuất: Giảm thiểu năng lượng làm lạnh nhân tạo cho trạm trộn bê tông, rút ngắn thời gian nghỉ giữa các khối đổ từ 5–7 ngày xuống còn 2–3 ngày, đẩy nhanh tiến độ lên đập mà vẫn kiểm soát an toàn tuyệt đối nứt nẻ.
  • Tác động kinh tế - xã hội: Tiết kiệm hàng trăm tỷ đồng chi phí xử lý vết nứt, sửa chữa chống thấm bằng giải pháp phụt vữa epoxy/microfine hóa chất đắt tiền, đồng thời bảo đảm an toàn tuyệt đối cho các cụm hồ chứa thủy điện bậc thang ở hạ du.
  • Ý nghĩa quốc tế: Kết quả nghiên cứu là tài liệu tham khảo giá trị cao cho các quốc gia trong khu vực Đông Nam Á, Nam Á và Châu Phi có đặc điểm khí hậu nhiệt đới gió mùa đang đẩy mạnh phát triển đập bê tông đầm lăn.

Đố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 phân tích đa trường vật lý kết hợp lập trình APDL nâng cao; khai thác các công thức thực nghiệm về động học nhiệt thủy hóa clinker nội địa.
  • Kỹ sư Thiết kế và Tư vấn Xây dựng Thủy lợi - Thủy điện (TVTK): Sở hữu công cụ số tin cậy và bảng tra cứu thông số nhiệt - ứng suất chuẩn hóa theo 3 miền địa lý để thiết kế cấp phối và phân bố khe co giãn hợp lý.
  • Chủ đầu tư và Ban Quản lý Dự án: Có cơ sở khoa học để lập biện pháp thi công, kiểm soát chất lượng vật liệu đầu vào (lựa chọn nhà máy xi măng có chỉ số $C_3A + C_3S$ phù hợp) và giám sát chặt chẽ tiến độ thi công.
  • Cơ quan Quản lý Nhà nước (Bộ Nông nghiệp & PTNT, Bộ Công Thương, Bộ Xây dựng): Ban hành các quy chuẩn kỹ thuật chuyên ngành về an toàn đập khối lớn sát thực với điều kiện Việt Nam.

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ì? Đó là việc thiết lập và giải mã định lượng mối quan hệ tương hỗ giữa thành phần khoáng clinker xi măng nội địa ($C_3A, C_3S$), tỷ lệ phụ gia khoáng PGK hoạt tính cao và trường ứng suất nhiệt sinh ra trong kết cấu đập BTĐL dưới điều kiện vi khí hậu nhiệt đới phân vùng, mở rộng lý thuyết nhiệt thủy hóa của Woods - McCarthy - Verbeck sang lĩnh vực bê tông đầm lăn nghèo xi măng.

2. Điểm cải tiến cốt lõi về mặt phương pháp luận so với các nghiên cứu trước đây? So với nghiên cứu của Barrett & Tatro (1992) hay Lê Quốc Toàn (2014), luận án đã tự động hóa quy trình phân tích nhiệt - ứng suất ghép đôi trên ANSYS bằng mã lệnh APDL độc quyền, áp dụng kỹ thuật phần tử "sinh - tử" mô phỏng chính xác từng lớp rải đầm $0,3\text{ m}$ theo thời gian thực và tích hợp đầy đủ các điều kiện biên bức xạ mặt trời $T_{ow}$ cũng như độ ẩm không khí.

3. Phát hiện bất ngờ và phản trực giác nhất từ dữ liệu thực nghiệm là gì? Trái với quan niệm thông thường cho rằng nhiệt độ môi trường mùa hè nắng nóng gây nguy cơ nứt lớn nhất, mô phỏng số và dữ liệu quan trắc thực tế tại Miền núi phía Bắc chứng minh rằng các đợt không khí lạnh đột ngột vào mùa đông mới là tác nhân nguy hiểm nhất tạo ra gradient nhiệt độ trong - ngoài cực lớn, gây ra ứng suất kéo vượt ngưỡng phá hoại tại bề mặt đập.

4. Quy trình nghiên cứu có cung cấp đầy đủ giao thức để tái lặp (replication protocol) không? Có hoàn toàn. Toàn bộ thông số đầu vào cơ lý, nhiệt học của bê tông, đá nền, các phương trình điều kiện biên, bảng thành phần khoáng xi măng (Phụ lục I) và toàn bộ tập lệnh chương trình tính toán APDL (Phụ lục II) đều được công khai chi tiết trong luận án.

5. Lộ trình nghiên cứu 10 năm tiếp theo được định hình như thế nào? Lộ trình hướng đến việc phát triển hệ thống giám sát và dự báo nứt nhiệt số hóa theo thời gian thực (Digital Twin) tích hợp mạng nơ-ron vật lý (PINNs) và dữ liệu cảm biến sợi quang (Fiber Optic Sensors) phân bố trong thân đập BTĐL.

Kết luận

Luận án tiến sĩ của NCS. Nguyễn Minh Việt là một công trình nghiên cứu khoa học thực nghiệm và tính toán số có giá trị học thuật xuất sắc và ý nghĩa thực tiễn to lớn, cô đọng qua 5 đóng góp nền tảng:

  1. Xác lập cơ sở khoa học định lượng về quá trình phát triển trường nhiệt độ và trường ứng suất nhiệt trong đập BTĐL theo thời gian thi công thực tế tại Việt Nam.
  2. Thiết lập mô hình toán nhiệt - cơ ghép đôi trong ANSYS APDL có xét đến ảnh hưởng của thành phần khoáng clinker xi măng ($C_3A, C_3S$) và phụ gia khoáng tro bay/puzơlan.
  3. Đưa ra hệ thống giải pháp phân vùng sinh khí hậu đồng bộ (Miền núi phía Bắc, Bắc Trung Bộ, Nam Trung Bộ & Tây Nguyên) với các ngưỡng khống chế cụ thể về tỷ lệ PGK ($35 \div 50%$), khoáng $(C_3A + C_3S \le 58%)$ và nhiệt độ hỗn hợp đổ ($T_p \le 20 \div 24^\circ\text{C}$).
  4. Kiểm nghiệm thành công trên đập thủy điện Sơn La và chuyển giao ứng dụng tính toán bảo đảm an toàn kháng nứt cho đập thủy điện Trung Sơn.
  5. Định hình chuẩn mực mới trong công tác thiết kế, kiểm soát ứng suất nhiệt và bảo đảm độ bền lâu dài cho các công trình đập bê tông đầm lăn tại Việt Nam và khu vực.