Tổng quan về luận án

Công nghệ bê tông đầm lăn (Roller Compacted Concrete - RCC/BTĐL) đại diện cho bước đột phá mang tính cách mạng trong ngành xây dựng đập trọng lực quy mô lớn trên thế giới từ thập niên 1970 đến nay. Bản chất của BTĐL là hỗn hợp bê tông không độ sụt, sử dụng hàm lượng xi măng thấp (dao động từ 60 đến 100 kg/m³), được đầm chặt bằng lu rung tương tự công nghệ thi công đắp đập đất đá. Nhờ tốc độ thi công vượt trội, tiết giảm chi phí cốp pha, giảm quy mô công trình phụ trợ dẫn dòng và hạ giá thành xây dựng từ 15% đến 30%, BTĐL đã được triển khai rộng rãi trên hơn 20 công trình thủy lợi - thủy điện trọng điểm tại Việt Nam như Pleikrông (71 m), Sơn La (138 m), Lai Châu (137 m), Bản Vẽ (138 m), Đồng Nai 3 (110 m), Đồng Nai 4 (129 m) và Trung Sơn (84 m).

Tuy nhiên, đặc thù khí hậu nhiệt đới gió mùa tại Việt Nam với nền nhiệt độ cao, bức xạ mặt trời gay gắt, độ ẩm biến động mạnh và biên độ nhiệt ngày - đêm cũng như mùa hè - mùa đông chênh lệch lớn đã tạo ra thách thức kỹ thuật nghiêm trọng: hiện tượng nứt khối bê tông do ứng suất nhiệt trong quá trình thủy hóa chất kết dính (CKD). Thực tiễn cho thấy nhiều đập BTĐL quy mô quốc gia tại Việt Nam đã xuất hiện các dạng nứt nguy hiểm (nứt bề mặt, nứt xuyên từ nền, nứt song song trục đập và nứt quanh hành lang kiểm tra), đe dọa trực tiếp đến tính chỉnh thể, khả năng chống thấm và tuổi thọ công trình.

Khoảng trống nghiên cứu (Research Gap) cốt lõi xuất phát từ việc các tiêu chuẩn thiết kế và giải pháp kiểm soát nhiệt hiện hành tại Việt Nam chủ yếu sao chép nguyên bản từ quy phạm của Hoa Kỳ (USACE, USBR, ACI 207.5R) hoặc Trung Quốc (SL314-2004, SL319-2005) – vốn được xây dựng trên điều kiện khí hậu ôn đới hoặc bán khô hạn và nguồn vật liệu bản địa khác biệt. Chưa có công trình nào giải quyết một cách toàn diện mối quan hệ tương tác giữa: (1) thành phần khoáng hoạt tính của xi măng nội địa Việt Nam ($C_3S, C_2S, C_3A, C_4AF$), (2) tỷ lệ thay thế của phụ gia khoáng (PGK) như tro bay nhiệt điện Phả Lại, Puzơlan tự nhiên Long Phước, và (3) đặc thù khí hậu phân vùng lãnh thổ để tối ưu hóa trường nhiệt độ và trường ứng suất nhiệt theo thời gian thực thi công.

Luận án tiến sĩ kỹ thuật của NCS. Nguyễn Minh Việt với đề tài "Nghiên cứu giải pháp giảm ứng suất nhiệt của bê tông đầm lăn trong xây dựng công trình thủy lợi thủy điện tại Việt Nam" (Viện Khoa học Thủy lợi Việt Nam, 2017) đã giải quyết triệt để bài toán khoa học này. Luận án đặt ra 03 câu hỏi nghiên cứu và giả thuyết tương ứng:

  • Câu hỏi nghiên cứu 1 (RQ1): Quy luật phát triển nhiệt thủy hóa thực tế của chất kết dính BTĐL phụ thuộc như thế nào vào hàm lượng các khoáng vật xi măng sản xuất tại Việt Nam ($C_3A + C_3S$) và tỷ lệ pha trộn phụ gia khoáng PGK/CKD?
    • Giả thuyết 1 (H1): Tích hợp thành phần khoáng thực tế vào phương trình nguồn nhiệt thủy hóa phi tuyến sẽ phản ánh chính xác trường nhiệt độ nội tại khối đập so với các công thức kinh nghiệm chuẩn hóa quốc tế.
  • Câu hỏi nghiên cứu 2 (RQ2): Mức độ tích tụ ứng suất kéo chính tại các vùng xung yếu (biên thượng lưu, hạ lưu, đáy tiếp xúc nền đá) diễn biến ra sao dưới tác động đồng thời của nhiệt độ môi trường chu kỳ, nhiệt độ đổ ban đầu ($T_p$) và quá trình phát triển mô đun đàn hồi phi tuyến theo thời gian?
    • Giả thuyết 2 (H2): Trạng thái ứng suất kéo nguy hiểm nhất xuất hiện tại bề mặt bóc lộ vào mùa lạnh do gradient nhiệt độ trong - ngoài ($\Delta T_{in-out}$) và tại đáy khối đổ do kiềm chế biên nền đá cứng ($R_n$).
  • Câu hỏi nghiên cứu 3 (RQ3): Bộ giải pháp công nghệ phân vùng nào cho phép kiểm soát ứng suất kéo nhỏ hơn cường độ kháng kéo giới hạn ($[\sigma_t] \le f_t(t)$) mà vẫn bảo đảm tính khả thi kinh tế - kỹ thuật tại từng vùng khí hậu Việt Nam?
    • Giả thuyết 3 (H3): Thiết lập ngưỡng tối ưu hóa cặp thông số $(C_3A + C_3S)$ xi măng và tỷ lệ PGK cho từng phân vùng địa lý sẽ loại trừ nguy cơ nứt nhiệt mà không cần hệ thống làm lạnh nhân tạo đắt đỏ.

Khung lý thuyết của luận án tích hợp chặt chẽ: Lý thuyết dẫn nhiệt vi phân Fourier-Laplace, Lý thuyết đàn hồi nhiệt môi trường liên tục (Thermo-elasticity Theory), Động học thủy hóa xi măng phi tuyến và Cơ học phá hủy bê tông khối lớn. Nghiên cứu thực hiện mô phỏng số chi tiết trên 02 công trình kiểm chứng thực tế: Thủy điện Sơn La (đập cao 138 m, khối lượng BTĐL xấp xỉ 3,0 triệu m³) và Thủy điện Trung Sơn (đập cao 84 m, khối lượng BTĐL trên 700.000 m³), đối chiếu trực tiếp với chuỗi số liệu quan trắc nhiệt kế nhiệt điện trở đa điểm theo thời gian thực.

Literature Review và Positioning

Nghiên cứu về bê tông đầm lăn và bài toán nhiệt khối lớn đã trải qua hơn 5 thập kỷ phát triển với các trường phái lý thuyết và công nghệ rõ rệt:

  1. Giai đoạn khởi xướng và hình thành công nghệ (1960 - 1980): Khởi đầu từ việc thử nghiệm bê tông không độ sụt rải bằng máy ủi tại đập Alpe Gera (Italia, 1961) và Manicouagan (Canada, 1961), bước ngoặt lý thuyết được thiết lập khi Raphael (1970) công bố báo cáo kinh điển "The Optimum Gravity Dam", đề xuất thay thế toàn bộ quy trình đổ đập truyền thống bằng công nghệ rải và đầm lèn liên tục. Tiếp đó, Cannon (1972, 1974) tại Tennessee Valley Authority (Hoa Kỳ) đã chuẩn hóa khái niệm bê tông nghèo xi măng đầm bằng lu rung, dẫn tới việc Công binh Lục quân Hoa Kỳ (USACE) thi công thực nghiệm tại đập Lost Creek (1977) và xây dựng đập BTĐL đầu tiên trên thế giới Willow Creek (Oregon, 1980, cao 52 m, dài 543 m, thể tích 331.000 m³ BTĐL). Tại Anh, Dunstan (1970) cùng Hiệp hội Nghiên cứu và Thông tin Công nghiệp Xây dựng (CIRIA) đã tiên phong phát triển trường phái BTĐL hàm lượng tro bay cao (High Fly-Ash RCC), làm tiền đề cho Cục Khai hoang Hoa Kỳ (USBR) thiết kế đập Upper Stillwater (1987, cao 90 m, 1,15 triệu m³ BTĐL).
  2. Sự phân hóa trường phái công nghệ quốc tế (1980 - nay): Thế giới hình thành 3 trường phái chính:
    • Trường phái Mỹ (RCC): Ưu tiên BTĐL nghèo chất kết dính ($CKD < 100\text{ kg/m}^3$), sử dụng màng chống thấm thượng lưu hoặc lớp bê tông thường (CVC) đổ đồng thời.
    • Trường phái Nhật Bản (RCD - Roller Compacted Dam-concrete): Phát triển từ đập Shimajigawa (1980, cao 89 m), áp dụng nguyên lý kết cấu "vàng bọc bạc" (lõi BTĐL, vỏ CVC dày 2–3 m chống thấm cao) và kiểm soát khe co giãn nghiêm ngặt.
    • Trường phái Trung Quốc (RCCD): Bắt đầu từ đập Khanh Khẩu (1986), phát triển mạnh mẽ về đập BTĐL siêu cao (Longtan cao 216 m, Baise cao 130 m) sử dụng BTĐL giàu chất kết dính kết hợp tro bay nghiền mịn và công nghệ cắt khe nhiệt cưỡng bức.
  3. Tiến trình mô hình hóa trường nhiệt và ứng suất nhiệt: Barret và Tatro (1992) là những người tiên phong ứng dụng phương pháp phần tử hữu hạn (FEM) để giải bài toán nhiệt - ứng suất có xét đến hiệu ứng từ biến (creep) và co ngót (shrinkage) theo thời gian. Các nghiên cứu tiếp theo của Zhu Bofang (1999), Cervera et al. (2000), và ACI Committee 207 (2007) đã hoàn thiện mô hình tương tác nhiệt - cơ học phi tuyến trong bê tông khối lớn.
  4. Các tranh luận khoa học đối lập (Scientific Debates):
    • Tranh luận 1 (Hàm lượng chất kết dính vs. Nhiệt phát sinh): Nhóm nghiên cứu theo trường phái USACE bảo vệ quan điểm giữ hàm lượng CKD cực thấp ($< 100\text{ kg/m}^3$) để khống chế triệt để nhiệt thủy hóa, chấp nhận cường độ chịu kéo thấp ($f_t < 1.0\text{ MPa}$). Ngược lại, trường phái Trung Quốc và Tây Âu chủ trương nâng cao CKD ($150 - 220\text{ kg/m}^3$) kết hợp tỷ lệ PGK lớn ($> 50%$) nhằm tăng tính chống thấm và cường độ dán dính lớp rải ($f_t > 1.8\text{ MPa}$), kiểm soát nhiệt bằng phụ gia khoáng hoạt tính.
    • Tranh luận 2 (Phương pháp kiểm soát nứt do nền kiềm chế): Một luồng quan điểm (Dunstan, USBR) đề xuất hạ thấp nhiệt độ đổ ($T_p \le 15 - 18^\circ\text{C}$) bằng trạm làm lạnh nước và đá vẩy. Luồng quan điểm khác (Lê Quốc Toàn 2012, Đỗ Văn Lượng 2008) cho rằng trong điều kiện nhiệt đới, việc hạ thấp $T_p$ quá sâu gây tốn kém năng lượng khổng lồ và làm tăng biên độ sốc nhiệt bề mặt khi tiếp xúc với không khí môi trường, thay vào đó cần tập trung kiểm soát tốc độ tỏa nhiệt của khoáng xi măng.

Vị thế học thuật của luận án được xác lập rõ nét: Luận án không dừng lại ở việc áp dụng đơn thuần các gói phần mềm thương mại để tính toán kiểm tra cho một dự án cá biệt (như các nghiên cứu trước đây tại Pleikrông của Lê Anh Vân 2005, hay Bản Chát của TVXD1 2011), mà đã tiên phong xây dựng khung phương pháp luận hoàn chỉnh từ vi mô (thành phần khoáng xi măng Việt Nam $C_3A, C_3S$) đến vĩ mô (trường ứng suất nhiệt kết cấu đập - nền đá tương tác), lượng hóa chính xác ảnh hưởng của các yếu tố khí hậu bản địa để đề xuất giải pháp kỹ thuật theo phân vùng địa lý đầu tiên tại Việt Nam.

Đó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 và làm sâu sắc thêm Lý thuyết truyền nhiệt Fourier trong môi trường sinh nhiệt nội vi biến thiên và Lý thuyết đàn hồi nhiệt phi tuyến (Thermo-viscoelasticity) trong vật liệu bê tông trẻ:

  1. Hiệu chỉnh mô hình động học nhiệt thủy hóa xi măng theo thành phần khoáng thực tế: Kế thừa các nghiên cứu kinh điển về enthalpy phân rã của Woods (1932), McCarthy (1985) và Verbeck (1968), luận án đã lượng hóa chính xác phương trình nhiệt sinh tổng cộng của xi măng Poóc-lăng sản xuất tại Việt Nam theo tỷ lệ 4 khoáng chất chủ đạo: $$Q(t) = Q_0 \cdot \left(1 - e^{-r t}\right)$$ Trong đó công suất sinh nhiệt cực đại $Q_0$ được xác định dựa trên hàm lượng phần trăm khoáng chất: $$Q_0 = a \cdot C_3S(%) + b \cdot C_2S(%) + c \cdot C_3A(%) + d \cdot C_4AF(%)\quad (\text{J/g})$$ với các trọng số thực nghiệm trung bình: $a = 520\text{ J/g}$, $b = 247\text{ J/g}$, $c = 1023\text{ J/g}$, $d = 337\text{ J/g}$. Luận án chứng minh rằng hai khoáng chất $C_3A$ và $C_3S$ chiếm hơn 75% tổng lượng nhiệt tỏa ra trong 72 giờ đầu, định hình tốc độ dâng nhiệt đột biến trong lòng khối đập.
  2. Mô hình hóa hàm suy giảm nhiệt lượng khi tích hợp phụ gia khoáng: Thiết lập hàm tương quan phi tuyến giữa tỷ lệ thay thế PGK/CKD (tro bay, Puzơlan) và độ suy giảm nhiệt dung tích, chứng minh cơ chế pha loãng nồng độ clinker kết hợp phản ứng Puzolanic thứ cấp làm dịch chuyển đỉnh nhiệt độ (peak temperature) chậm lại từ 3 đến 7 ngày, tạo điều kiện thuận lợi cho quá trình tản nhiệt tự nhiên qua bề mặt.
  3. Mở rộng mô hình phát triển cơ lý tính theo thời gian thực: Luận án đã chuẩn hóa hệ phương trình thực nghiệm phản ánh sự phát triển đồng thời của cường độ chịu nén ($f_c$), cường độ kháng kéo ($f_t$) và mô đun đàn hồi ($E_c$) theo tỷ lệ nước/chất kết dính ($N/CKD$) và hàm lượng $PGK/CKD$ ở các mốc tuổi kiểm soát 28 ngày và 90 ngày: $$f_c^{28} = 19{,}459\left(\frac{CKD}{N}\right) - 0{,}147\left(\frac{PGK}{CKD}\right) - 11{,}681\quad (\text{MPa})$$ $$f_c^{90} = 19{,}326\left(\frac{CKD}{N}\right) - 0{,}333\left(\frac{PGK}{CKD}\right) + 5{,}968\quad (\text{MPa})$$ $$E_c^{90} = 21{,}217\left(\frac{CKD}{N}\right) - 0{,}197\left(\frac{PGK}{CKD}\right) - 1{,}376\quad (\text{GPa})$$
+-----------------------------------------------------------------------------+
|                     KHUNG PHÂN TÍCH NHIỆT - ỨNG SUẤT BTĐL                   |
+-----------------------------------------------------------------------------+
|  [ĐẦU VÀO VẬT LIỆU & KHÍ HẬU]                                               |
|  - Khoáng Xi măng: C3A, C3S, C2S, C4AF                                      |
|  - Tỷ lệ Phụ gia khoáng: PGK/CKD (Tro bay / Puzơlan)                        |
|  - Khí hậu 3 vùng: Tm (trung bình), Tb (biên độ), Độ ẩm H%                  |
|  - Nhiệt độ đổ ban đầu: Tp                                                  |
+-----------------------------------------------------------------------------+
+-----------------------------------------------------------------------------+
|  [BÀI TOÁN TRUYỀN NHIỆT TRANSIENT 2D (PLANE55)]                             |
|  Phương trình Fourier-Laplace: c.rho.(dT/dt) = div(lambda.grad T) + qv      |
|  - Biên loại 1: Nhiệt độ tiếp xúc nền đá                                    |
|  - Biên loại 3: Tỏa nhiệt đối lưu beta.(Tw - Tf) + Bức xạ mặt trời Q0       |
|  ==> Xuất Trường Nhiệt độ T(x, y, t) theo tiến độ rải lớp                   |
+-----------------------------------------------------------------------------+
+-----------------------------------------------------------------------------+
|  [BÀI TOÁN ĐÀN HỒI NHIỆT KẾT CẤU TƯƠNG TÁC (PLANE182)]                      |
|  Biến dạng tổng: epsilon = epsilon_e + epsilon_th + epsilon_cr              |
|  - Module biến dạng phi tuyến: Ec(t), ft(t)                                 |
|  - Kiềm chế biên nền đá & khối đổ cũ (Rn)                                   |
|  ==> Xuất Trường Ứng suất chính lớn nhất sigma_1(x, y, t)                   |
+-----------------------------------------------------------------------------+
+-----------------------------------------------------------------------------+
|  [TIÊU CHÍ AN TOÀN KHỐNG CHẾ NỨT NHIỆT]                                     |
|  - Delta T_đáy = Tp + Tr - Tf <= [Delta T_đáy] (Quy phạm SL314-2004)        |
|  - TM <= [15 - 20 oC] + Tm - Tb (Quy phạm SL319-2005)                       |
|  - Điều kiện bền phá hủy: sigma_1(t) <= ft(t) (Hệ số an toàn K_t >= 1.2)    |
+-----------------------------------------------------------------------------+

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

Khung phân tích của luận án tích hợp liên môn giữa Hóa học Silicate, Cơ học Môi trường liên tục và Phương pháp Phần tử hữu hạn chuyên sâu:

  • Tiếp cận phân tích đa quy mô (Multi-scale): Kết nối từ cấp độ cấu trúc vi mô của phản ứng thủy hóa khoáng xi măng sang bài toán ứng xử cơ - nhiệt vĩ mô của toàn bộ khối đập làm việc đồng thời với nền đá.
  • Mô hình liên kết tuần tự (Sequentially Coupled Thermal-Stress Formulation): Trường nhiệt độ $T(x,y,t)$ tại từng nút phần tử được tính toán từ phân tích nhiệt quá độ (Transient Thermal Analysis), sau đó được nạp trực tiếp làm véc-tơ tải trọng nhiệt (Thermal Body Load) để giải bài toán biến dạng và ứng suất phẳng (Plane Strain) tại cùng bước thời gian.
  • Điều kiện biên xác định (Boundary Conditions):
    • Biên bức xạ và đối lưu không khí - bê tông: $q = \beta (T_w - T_f) - A_s Q_0$, với hệ số trao đổi nhiệt đối lưu $\beta$ phụ thuộc trực tiếp vào vận tốc gió và độ ẩm không khí vùng lòng hồ thủy điện.
    • Biên truyền nhiệt dẫn qua nền đá: Xét chiều sâu vùng ảnh hưởng nhiệt của nền đá bằng 1,0 đến 1,5 lần chiều cao đập.

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo trường phái Thực chứng luận (Positivism) kết hợp Thực nghiệm số (Numerical Experimentation) và kiểm chứng thực địa (Field Instrumentation Validation):

  • Triết lý nghiên cứu: Mọi quy luật toán học về phân bố ứng suất nhiệt đều phải được suy diễn từ định luật nhiệt động lực học bảo toàn năng lượng và được xác nhận thông qua dữ liệu đo lường thực tế tại công trình cấp đặc biệt.
  • Quy mô mẫu và thông số mô phỏng: Khảo sát toàn diện 21 đập BTĐL tại Việt Nam, thu thập mẫu thí nghiệm thành phần khoáng của 6 nhà máy xi măng lớn (Bút Sơn, Hoàng Thạch, Bỉm Sơn, Nghi Sơn, Tam Điệp, Thăng Long), kiểm tra cơ lý 180 tổ mẫu bê tông BTĐL sử dụng tro bay Phả Lại và Puzơlan Long Phước ở các tỷ lệ thay thế từ 0% đến 50%.
  • Mô hình hóa hình học: Thiết lập lưới phần tử hữu hạn 2D mặt cắt ngang đập đại diện lớn nhất (Mặt cắt tràn xả lũ và tuyến đập không tràn), mô phỏng tuần tự từng lớp rải $0{,}3\text{ m}$, chia đợt đổ 3–5 lớp ($0{,}9 - 1{,}5\text{ m}$) kèm thời gian nghỉ công nghệ từ 2 đến 7 ngày.

Quy trình nghiên cứu rigorous

Quy trình tính toán số được tự động hóa hoàn toàn bằng Ngôn ngữ lập trình tham số APDL (ANSYS Parametric Design Language) trong phần mềm ANSYS:

[BƯỚC 1: TIỀN XỬ LÝ NHIỆT]
[BƯỚC 2: MÔ PHỎNG QUÁ ĐỘ NHIỆT THỦY HÓA & THI CÔNG]
[BƯỚC 3: CHUYỂN ĐỔI BÀI TOÁN NHIỆT - CƠ HỌC]
[BƯỚC 4: TÍNH TOÁN TRƯỜNG ỨNG SUẤT NHIỆT & KIỂM ĐỊNH NỨT]

Tính hợp thức và độ tin cậy (Validity & Reliability) của mô hình được kiểm chứng nghiêm ngặt qua độ trùng khớp giữa trường nhiệt độ mô phỏng số và số liệu nhiệt kế điện trở PT64 tại đập Sơn La (đoạn đập cao trình +114,0 m) và đập Trung Sơn: Sai số tuyệt đối $|\Delta T| \le 1{,}2^\circ\text{C}$, sai số tương đối dưới 3,5%, chứng minh thuật toán mô phỏng có độ tin cậy tuyệt đối.

Data và phân tích

Các chỉ tiêu cơ lý và thông số nhiệt nhiệt động lực học đầu vào chuẩn hóa trong mô hình tính toán:

  • Bê tông đầm lăn (BTĐL): Khối lượng thể tích $\rho = 2400\text{ kg/m}^3$; Hệ số dẫn nhiệt $\lambda = 1{,}75\text{ kcal/(m}\cdot\text{h}\cdot^\circ\text{C)}$ ($2{,}035\text{ W/m}\cdot\text{K}$); Nhiệt dung riêng $C = 0{,}24\text{ kcal/(kg}\cdot^\circ\text{C)}$ ($1004{,}8\text{ J/kg}\cdot\text{K}$); Hệ số khuếch tán nhiệt $a = 0{,}00304\text{ m}^2/\text{h}$; Hệ số dãn nở nhiệt $\alpha = 0{,}85 \times 10^{-5}\text{ /}^\circ\text{C}$; Hệ số Poisson $\nu = 0{,}18 - 0{,}20$.
  • Nền đá móng (Granite/Gneiss cường độ cao): $\rho_r = 2650\text{ kg/m}^3$; $\lambda_r = 2{,}10\text{ kcal/(m}\cdot\text{h}\cdot^\circ\text{C)}$; $C_r = 0{,}22\text{ kcal/(kg}\cdot^\circ\text{C)}$; Mô đun đàn hồi nền $E_r = 18 - 25\text{ GPa}$.
  • Phân tích độ nhạy (Sensitivity Analysis): Thực hiện tính toán lặp trên 54 kịch bản tổ hợp tham số biến thiên:
    • Hàm lượng khoáng $(C_3A + C_3S)$ xi măng: $48%, 52%, 56%, 60%, 65%$.
    • Tỷ lệ thay thế PGK: $0%, 20%, 30%, 35%, 40%, 45%, 50%$.
    • Nhiệt độ hỗn hợp bê tông khi đổ $T_p$: $18^\circ\text{C}, 20^\circ\text{C}, 22^\circ\text{C}, 25^\circ\text{C}, 28^\circ\text{C}$.
    • Nhiệt độ môi trường trung bình $T_m$: $18^\circ\text{C}$ (Miền Bắc mùa đông) đến $32^\circ\text{C}$ (Nam Trung Bộ mùa khô).

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

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

Kết quả mô phỏng số và phân tích dữ liệu thực nghiệm đã mang lại 5 phát hiện mang tính đột phá cho ngành xây dựng đập thủy điện:

  1. Hiệu ứng chi phối của tổ hợp khoáng $(C_3A + C_3S)$ đến ứng suất kéo bề mặt: Hàm lượng tổng khoáng $(C_3A + C_3S)$ trong xi măng tăng từ 50% lên 65% làm nhiệt độ đỉnh trong thân đập tăng thêm $6{,}8^\circ\text{C}$ đến $9{,}2^\circ\text{C}$, kéo theo ứng suất kéo chính lớn nhất ($\sigma_1$) tại mép thượng lưu và hạ lưu tăng vọt từ $1{,}24\text{ MPa}$ lên $2{,}15\text{ MPa}$ (vượt quá cường độ kháng kéo giới hạn $f_t^{28} \approx 1{,}45\text{ MPa}$, dẫn đến phá hoại nứt bề mặt).
  2. Quy luật giải phóng ứng suất nhiệt thông qua tỷ lệ PGK tối ưu: Tăng hàm lượng tro bay hoặc Puzơlan từ 20% lên 40% trong tổng lượng CKD giúp giảm nhiệt độ cực đại trong thân đập từ $46{,}5^\circ\text{C}$ xuống $38{,}2^\circ\text{C}$ ($\Delta T_{max} = -8{,}3^\circ\text{C}$), đồng thời làm giảm mô đun đàn hồi sớm $E_c(t)$, tạo điều kiện cho hiện tượng từ biến (creep) thư giãn ứng suất kéo, đưa $\sigma_1$ tại vùng đáy đập từ trạng thái nguy hiểm $1{,}85\text{ MPa}$ về ngưỡng an toàn $0{,}92\text{ MPa}$ ($\text{giảm } 50{,}3%$).
  3. Phát hiện nghịch lý về nhiệt độ đổ bê tông ($T_p$) tại vùng nhiệt đới: Hạ thấp nhiệt độ đổ $T_p$ quá sâu ($< 18^\circ\text{C}$) trong mùa hè nắng nóng không làm giảm đáng kể nguy cơ nứt mà trái lại làm tăng gradient nhiệt độ cục bộ tại lớp bề mặt tiếp xúc với không khí nóng ($35 - 38^\circ\text{C}$), tạo ra ứng suất kéo bề mặt tức thời rất cao ngay trong 48 giờ đầu sau khi rải. Ngưỡng $T_p$ an toàn và kinh tế nhất tại Việt Nam được xác định là $22^\circ\text{C} \le T_p \le 25^\circ\text{C}$.
  4. Cơ chế nứt tiếp xúc nền đá do kiềm chế biên: Ứng suất kéo nguy hiểm nhất tập trung tại dải chiều cao $(0{,}0 \div 0{,}2)H$ tính từ đáy đập tiếp xúc nền đá cứng, phát sinh chủ yếu trong giai đoạn hạ nhiệt ổn định (từ ngày thứ 60 đến 180). Chênh lệch nhiệt độ đáy đập $\Delta T_{đáy} = T_p + T_r - T_f$ bắt buộc phải khống chế không vượt quá $16 - 19^\circ\text{C}$ đối với khối đổ toàn mặt cắt rộng trên $30\text{ m}$.
  5. Bản đồ phân vùng giải pháp nhiệt công nghệ đập BTĐL tại Việt Nam:
Phân vùng địa lý Đặc trưng khí hậu & Nguồn vật liệu Yêu cầu khoáng XM $(C_3A + C_3S)$ Tỷ lệ PGK / CKD Nhiệt độ đổ khống chế $T_p$
Vùng 1: Miền núi phía Bắc $T_m = 22 - 24^\circ\text{C}$, mùa đông lạnh ($< 12^\circ\text{C}$), nguồn tro bay Phả Lại, Hải Phòng dồi dào $\le 55%$ (Ưu tiên XM ít tỏa nhiệt) $35% \div 45%$ (Tro bay tuyển) $\le 22^\circ\text{C}$ (Mùa hè làm mát nước trộn)
Vùng 2: Bắc Trung Bộ $T_m = 25 - 26^\circ\text{C}$, gió phơn Tây Nam khô nóng, mùa đông rét buốt $\le 52%$ $30% \div 40%$ (Tro bay Nghi Sơn/Vũng Áng) $\le 23^\circ\text{C}$
Vùng 3: Nam Trung Bộ & Tây Nguyên $T_m = 26 - 28^\circ\text{C}$, nền nhiệt cao quanh năm, Puzơlan tự nhiên phong phú $\le 58%$ $35% \div 45%$ (Puzơlan mịn) $\le 25^\circ\text{C}$ (Che chắn bức xạ mặt trời)

Implications đa chiều

  • Về mặt học thuật: Cung cấp phương pháp luận hoàn chỉnh xác định trường nhiệt - ứng suất liên kết cho kết cấu bê tông khối lớn trong điều kiện khí hậu nhiệt đới gió mùa, đóng góp dữ liệu thực nghiệm quý giá cho kho tàng tài liệu của Hội Đập lớn Thế giới (ICOLD).
  • Về mặt công nghệ thi công: Loại bỏ sự phụ thuộc vào các giải pháp đắt tiền như đặt ống nước làm lạnh ngầm trong thân đập (Post-cooling pipes), thay thế hoàn toàn bằng giải pháp điều phối cấp phối vật liệu (Pre-cooling và điều chỉnh khoáng chất).
  • Về mặt chính sách và tiêu chuẩn hóa: Là cơ sở khoa học trực tiếp để Hội đồng Khoa học Bộ Nông nghiệp & PTNT, Bộ Công Thương và Bộ Xây dựng soát xét, cập nhật Tiêu chuẩn Quốc gia TCVN về "Thiết kế và thi công đập bê tông đầm lăn".

Limitations và Future Research

Mặc dù đạt được những đóng góp học thuật xuất sắc, luận án thẳng thắn chỉ rõ 03 giới hạn nghiên cứu:

  1. Giới hạn mô hình hình học phẳng (2D Plane Strain): Chưa xét đến hiệu ứng không gian 3 chiều tại các khối đổ cong quanh vai đập, vùng giao cắt giữa khối tràn xả lũ và trụ pin, nơi xuất hiện trạng thái ứng suất tập trung 3 trục phức tạp ($\sigma_1, \sigma_2, \sigma_3$).
  2. Giả thiết tính chất nhiệt bất biến: Hệ số dẫn nhiệt ($\lambda$) và nhiệt dung riêng ($C$) được giả thiết là hằng số không phụ thuộc vào độ ẩm nội tại và mức độ thủy hóa của bê tông trẻ.
  3. Mô hình từ biến tuyến tính: Chưa tích hợp mô hình từ biến phi tuyến hóa dẻo bậc cao (Nonlinear Viscoelastic-plastic Damage Model) dưới tác động của tải trọng thủy tĩnh động khi tích nước hồ chứa sớm.

Chương trình nghiên cứu tiếp nối trong giai đoạn 5 - 10 năm:

  • Hướng 1: Phát triển phần mềm tính toán nhiệt - ứng suất 3D tích hợp công nghệ trí tuệ nhân tạo (AI/Machine Learning) để dự báo thời gian thực nguy cơ nứt đập dựa trên mạng lưới cảm biến IoT không dây.
  • Hướng 2: Nghiên cứu ứng dụng tro xỉ lò cao nghiền mịn (GGBS) và tro xỉ đáy kết hợp Nano-silica trong BTĐL để giảm hơn nữa lượng nhiệt thủy hóa ban đầu.
  • Hướng 3: Đánh giá ảnh hưởng của biến đổi khí hậu cực đoan (nhiệt độ tăng cao, hạn hán kéo dài) đến độ bền lâu và trạng thái ứng suất nhiệt vận hành dài hạn của hệ thống đập BTĐL tại Việt Nam.

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

  • Tác động học thuật: Luận án được đánh giá là công trình tiên phong có tính hệ thống cao nhất tại Việt Nam về ứng suất nhiệt BTĐL, cung cấp nền tảng 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, Địa kỹ thuật và Vật liệu Xây dựng.
  • Tác động kinh tế - kỹ thuật: Ứng dụng giải pháp tối ưu hóa cấp phối khoáng và tỷ lệ PGK tại công trình Thủy điện Trung Sơn (Thanh Hóa) đã giúp chủ đầu tư và nhà thầu tiết kiệm hàng chục tỷ đồng chi phí trạm lạnh làm mát cốt liệu, rút ngắn tiến độ thi công toàn dự án thêm 3 tháng mà không xuất hiện bất kỳ vết nứt nhiệt nguy hiểm nào qua 5 năm tích nước vận hành.
  • Lợi ích xã hội và môi trường: Khuyến khích tận dụng tối đa hàng triệu tấn tro bay phế thải từ các nhà máy nhiệt điện than (Phả Lại, Nghi Sơn, Vũng Áng) làm chất kết dính xây đập, giảm thiểu phát thải khí nhà kính $CO_2$ từ ngành sản xuất xi măng và giải phóng quỹ đất chứa xỉ thải.

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

  • Nghiên cứu sinh & Nhà khoa học: Tiếp cận phương pháp lập trình APDL ANSYS chuyên sâu mô phỏng liên kết nhiệt - cơ học và phương trình động học nhiệt khoáng xi măng.
  • Kỹ sư Tư vấn Thiết kế (TVTK): Nắm vững bộ công cụ tính toán kiểm tra nứt nhiệt, lựa chọn chính xác các chỉ tiêu chênh lệch nhiệt độ cho phép ($\Delta T_{đáy}$, $\Delta T_{trong-ngoài}$) theo tiêu chuẩn quốc tế.
  • Nhà thầu Thi công đập: Ứng dụng quy trình kiểm soát nhiệt độ đổ bê tông ($T_p$), chiều dày lớp rải, khoảng thời gian nghỉ giữa các lớp đầm lèn phù hợp với thời tiết thực tế từng mùa.
  • Chủ đầu tư & Cơ quan Quản lý Nhà nước: Có cơ sở khoa học để nghiệm thu chất lượng công trình, bảo đảm an toàn tuyệt đối cho các hồ chứa thủy điện bậc thang phòng chống thảm họa vỡ đập.

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?

Trả lời: Đóng góp độc đáo nhất là việc mở rộng Lý thuyết nhiệt thủy hóa của Woods (1932) và Verbeck (1968) bằng cách thiết lập phương trình nguồn nhiệt nội sinh $Q_v(t)$ hiệu chỉnh riêng cho xi măng Poóc-lăng Việt Nam có xét đến tỷ lệ hoạt tính của cặp khoáng chất sinh nhiệt nhanh $(C_3A + C_3S)$ kết hợp hàm suy giảm nhiệt lượng phi tuyến do tỷ lệ phụ gia khoáng PGK/CKD chi phối, giải quyết triệt để sự thiếu chính xác của các mô hình truyền nhiệt kinh nghiệm truyền thống.

2. Đột phá về phương pháp luận nghiên cứu so với các công trình quốc tế trước đây?

Trả lời: So với Barret & Tatro (1992) hay Zhu Bofang (1999), luận án đã tự động hóa quy trình phân tích nhiệt - kết cấu liên kết tuần tự trên ANSYS APDL bằng thuật toán Sinh - Diệt phần tử (Element Birth & Death) bám sát tuyệt đối tiến độ thi công rải từng lớp mỏng $0{,}3\text{ m}$, đồng thời tích hợp hàm điều kiện biên khí hậu thực tế biến thiên liên tục theo chu kỳ ngày - đêm và mùa năm tại tọa độ công trình.

3. Phát hiện bất ngờ nhất có số liệu thực nghiệm chứng minh là gì?

Trả lời: Việc hạ thấp nhiệt độ đổ bê tông $T_p$ xuống dưới $18^\circ\text{C}$ vào mùa hè không giúp loại trừ nứt mà ngược lại làm tăng đột biến ứng suất kéo bề mặt $\sigma_1 > 1{,}8\text{ MPa}$ trong 48 giờ đầu do hiện tượng "sốc nhiệt" (Thermal Shock) giữa bề mặt ngoài tiếp xúc không khí nóng ($> 36^\circ\text{C}$) và khối lõi bên trong, từ đó chứng minh nhiệt độ đổ tối ưu về mặt cơ học nhiệt là $22 - 25^\circ\text{C}$.

4. Luận án có cung cấp quy trình tái lập nghiên cứu (Replication Protocol) không?

Trả lời: Toàn bộ mã nguồn chương trình tính toán APDL, bảng tra cứu thành phần khoáng xi măng tại Phụ lục I, thông số nhiệt động lực học vật liệu và quy trình thiết lập điều kiện biên trong ANSYS được trình bày chi tiết và minh bạch tại Chương 3, Chương 4 và Phụ lục II, cho phép các kỹ sư và nhà nghiên cứu tái lập 100% kết quả mô phỏng.

5. Lộ trình nghiên cứu phát triển trong 10 năm tới được vạch ra như thế nào?

Trả lời: Tập trung nâng cấp mô hình lên bài toán 3D phi tuyến đầy đủ, tích hợp ảnh hưởng của từ biến phi tuyến dẻo - phá hủy, nghiên cứu ứng dụng vật liệu xanh (Geopolymer, tro xỉ hạt lò cao nghiền mịn GGBS) và phát triển hệ thống cảnh báo sớm rủi ro nứt nhiệt đập thông minh dựa trên bản sao số (Digital Twin) và dữ liệu cảm biến quang Fiber Optic theo thời gian thực.

Kết luận

  1. Luận án đã giải quyết thành công và trọn vẹn bài toán cấp bách về kiểm soát nhiệt và giảm ứng suất nhiệt chống nứt cho đập bê tông đầm lăn trong điều kiện đặc thù của Việt Nam.
  2. Thiết lập thành công phương trình nguồn nhiệt thủy hóa chính xác dựa trên hàm lượng khoáng xi măng nội địa ($C_3A, C_3S, C_2S, C_4AF$) và tỷ lệ phụ gia khoáng hoạt tính (Tro bay, Puzơlan).
  3. Xây dựng chương trình tính toán số tin cậy trên nền tảng ANSYS APDL, mô phỏng chính xác trạng thái ứng suất - biến dạng không gian 2 chiều tương tác đập - nền đá theo thời gian thực thi công.
  4. Đề xuất bản đồ phân vùng giải pháp công nghệ giảm ứng suất nhiệt khoa học cho 3 miền địa lý của Việt Nam, kiểm chứng thành công tại đập Sơn La và đập Trung Sơn.
  5. Tạo bước chuyển dịch quan trọng từ phương pháp làm mát cưỡng bức thụ động sang phương pháp chủ động điều phối vật liệu và kiểm soát tiến độ công nghệ, mang lại giá trị kinh tế hàng trăm tỷ đồng và bảo đảm an toàn vĩnh cửu cho hệ thống công trình thủy lợi - thủy điện quốc gia.