Tổng quan về luận án
Công nghệ bê tông đầm lăn (Roller Compacted Concrete - BTĐL/RCC) đã tạo ra bước đột phá trong ngành xây dựng thủy lợi - thủy điện toàn cầu nhờ tốc độ thi công vượt trội, tối ưu hóa chi phí cốp pha và giảm thiểu hệ thống phụ trợ dẫn dòng. Tại Việt Nam, công nghệ này đã được triển khai trên hơn 20 công trình thủy điện trọng điểm quy mô lớn như Sơn La (cao 138 m), Lai Châu (137 m), Bản Vẽ (138 m), Đồng Nai 4 (129 m) và Pleikrông (71 m). Tuy nhiên, đặc thù khí hậu nhiệt đới gió mùa với biên độ dao động nhiệt độ ngày - đêm và mùa đông - mùa hè lớn, kết hợp bức xạ mặt trời gay gắt đã dẫn đến hiện tượng nứt do ứng suất nhiệt nghiêm trọng tại các khối đổ bê tông khối lớn, điển hình như sự cố nứt tại đập thủy điện Sơn La.
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 đập BTĐL tại Việt Nam chủ yếu áp dụng trực tiếp các tiêu chuẩn ngoại nhập từ Hoa Kỳ (USACE, USBR, ACI) hoặc Trung Quốc (SL314-2004, SL319-2005). Các hướng dẫn thiết kế này giả định điều kiện biên nhiệt tĩnh và đặc tính nhiệt thủy hóa của xi măng ngoại quốc, hoàn toàn không tương thích với thành phần khoáng học của các chủng loại xi măng nội địa cũng như điều kiện khí hậu đặc thù từng vùng miền Việt Nam. Luận án tiến sĩ 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" đã giải quyết triệt để bài toán này thông qua việc thiết lập mô hình tính toán nhiệt - ứng suất nhiệt động lực học không ổn định.
Nghiên cứu tập trung giải quyết hai câu hỏi nghiên cứu (Research Questions) và kiểm chứng hai giả thuyết khoa học (Hypotheses):
- RQ1: Sự biến thiên thành phần khoáng học ($C_3S, C_2S, C_3A, C_4AF$) trong xi măng nội địa và hàm lượng phụ gia khoáng (PGK) tro bay/Puzơlan ảnh hưởng như thế nào đến hàm nhiệt thủy hóa thực tế của BTĐL theo thời gian?
- RQ2: Quy luật phân bố trường nhiệt độ và trường ứng suất nhiệt không ổn định trong thân đập thay đổi ra sao dưới tác động đồng thời của điều kiện biên nhiệt biến thiên theo chu kỳ thời tiết và tiến trình thi công nhiều lớp?
- H1: Việc hiệu chỉnh phương trình phát nhiệt thủy hóa dựa trên tổng hàm lượng khoáng sinh nhiệt $(C_3S + C_3A)$ của xi măng sản xuất tại Việt Nam sẽ cho phép mô phỏng chính xác trường nhiệt độ khối đổ với sai số dưới $5%$ so với quan trắc thực tế.
- H2: Tối ưu hóa tỷ lệ thay thế phụ gia khoáng (PGK/CKD từ $30%$ đến $50%$) kết hợp kiểm soát công nghệ đầm lèn phân lớp và nhiệt độ đổ hỗn hợp ($T_p \le 23^\circ\text{C}$) sẽ triệt tiêu ứng suất kéo cục bộ vượt giới hạn bền chịu kéo ($f_t$), ngăn chặn hoàn toàn nứt xuyên và nứt bề mặt.
Khung lý thuyết của nghiên cứu được xây dựng trên sự tích hợp giữa Lý thuyết truyền nhiệt Fourier phi tuyến, Lý thuyết đàn hồi nhiệt Duhamel-Neumann và Động học thủy hóa xi măng Poóc-lăng. Phạm vi nghiên cứu bao quát dữ liệu cơ lý - nhiệt của 12 nhà máy xi măng lớn tại Việt Nam, phân tích chuyên sâu cho các kịch bản thi công đập trọng lực, và ứng dụng kiểm chứng thực chứng tại đập BTĐL Trung Sơn (tỉnh Thanh Hóa, chiều cao 84 m).
Literature Review và Positioning
Lịch sử phát triển BTĐL thế giới ghi nhận các cột mốc quan trọng từ báo cáo "Đập trọng lực tối ưu" của Giáo sư Raphael (1970), các nghiên cứu bê tông nghèo xi măng của Cannon R.W (1972–1974), đến các công trình thực nghiệm của Dunstan (Anh, thập niên 1970) và Hiệp hội Công binh Hoa Kỳ (USACE) tại đập Lost Creek (1977) và Willow Creek (Oregon, Mỹ, 1980 - cao 52 m). Qua hơn 4 thập kỷ, ba trường phái công nghệ BTĐL chính đã định hình:
- Trường phái Mỹ (RCC nghèo xi măng - USACE/USBR): Sử dụng hàm lượng chất kết dính (CKD) dưới $100\text{ kg/m}^3$, đập Willow Creek và Upper Stillwater (cao 90 m, khối lượng 1,17 triệu $\text{m}^3$) làm điển hình, chống thấm bằng màng thượng lưu hoặc bê tông CVC truyền thống.
- Trường phái Nhật Bản (RCD - Roller Compacted Dam): Triển khai lần đầu tại đập Shimajigawa (1974, cao 89 m), áp dụng nguyên lý kết cấu "vàng bọc bạc" với lõi là BTĐL và lớp vỏ bao bọc bằng bê tông thường CVC chống thấm cao dày $2\div3\text{ m}$.
- Trường phái Trung Quốc (RCCD): Bắt đầu từ đập Khanh Khẩu (1986) đến các đại công trình như Longtan (cao 216,5 m), thiên về sử dụng hàm lượng CKD trung bình và giàu, kết hợp hàm lượng tro bay cao ($>30%$) thi công trên toàn mặt cắt.
Trong lĩnh vực phân tích nhiệt và ứng suất nhiệt, Barret và Tatro (1992) đã tiên phong đặt nền móng giải bài toán phần tử hữu hạn (PTHH) mô phỏng tương tác giữa nhiệt thủy hóa, từ biến (creep) và co ngót (shrinkage). Tại Việt Nam, nhiều nhà khoa học đã tiếp cận từng khía cạnh đơn lẻ: Nguyễn Như Quý (2006) đo đạc nhiệt độ cách nhiệt của BTĐL sử dụng tro bay nhiệt điện ($\Delta T$ giảm từ $24,9^\circ\text{C}$ xuống $8,1\div12,7^\circ\text{C}$); Đỗ Hồng Hải (2007) khảo sát thay thế $35\div40%$ Puzơlan Long Phước cho đập Lòng Sông; Lê Quang Hùng và Nguyễn Quang Hiệp (2008) kiến nghị tỷ lệ tro bay $30\div50%$ cho BTĐL đập thủy điện; Lê Anh Vân, Nguyễn Văn Mạo (2006) và Đinh Hữu Dụng (2010) phân tích nứt và trường nhiệt đập Pleikrông, Sơn La; Lê Quốc Toàn (2014) khảo sát ảnh hưởng tốc độ thi công và mô-đun hóa đầu vào trên ANSYS.
Mặc dù có nhiều công trình nghiên cứu, y văn tồn tại tranh luận gay gắt giữa hai luồng quan điểm: (1) Nhóm quan điểm cho rằng việc kiểm soát nhiệt độ chỉ cần tối ưu hóa cấp phối tĩnh và áp dụng tiêu chuẩn nhiệt độ đổ cố định $T_p \le 23^\circ\text{C}$; (2) Nhóm quan điểm cho rằng hiện tượng nứt phụ thuộc chủ yếu vào động học giải phóng nhiệt và độ chênh gradient nhiệt độ tức thời giữa các lớp khối đổ với môi trường ngoài. Luận án đã định vị chính xác khoảng trống: Chưa có công trình nào lượng hóa sự biến thiên thành phần khoáng học đặc thù của xi măng sản xuất tại Việt Nam ($C_3S$ dao động $34,95\div81,70%$, $C_3A$ dao động $6,14\div15,00%$) vào phương trình động học phát nhiệt thủy hóa, đồng thời tích hợp bài toán biên nhiệt thay đổi theo chu kỳ tự nhiên trong môi trường phần tử hữu hạn đa trường vật lý.
Đó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 vật liệu chất kết dính của Lưu Thụ Hoa (Lu Shuhua) vốn xây dựng từ thực nghiệm xi măng Trung Quốc. Phương trình gốc của Lưu Thụ Hoa:
$$Q = c \cdot k_F \cdot Q_t$$
với $c = 0,008$, $k_F = -1,89F + 238,07$ ($F$ là tỷ lệ PGK/CKD tính bằng $%$), và $Q_t = 36,75\ln(t) + 101,15\text{ (J/g)}$ ($R^2 = 0,9899$).
Nhận thấy công thức trên đánh đồng đặc tính sinh nhiệt của mọi loại xi măng, luận án đã thiết lập đột phá lý thuyết bằng việc bổ sung hệ số ảnh hưởng khoáng học $k_C$, chuyển đổi phương trình phát nhiệt thủy hóa thành:
$$Q = c \cdot k_C \cdot k_F \cdot Q_t$$
Hệ số $k_C$ được phát triển dựa trên mối tương quan tuyến tính giữa tổng hàm lượng khoáng phát nhiệt chủ đạo $C_X = %(C_3S + C_3A)$ của xi măng nội địa với tỷ lệ phát nhiệt chuẩn:
$$k_C = 0,0073 \cdot C_X + 0,4842 \quad (R^2 = 0,9893)$$
kC (Hệ số khoáng học)
40 60 80 100
Mô hình lý thuyết này chứng minh rằng cùng một hàm lượng xi măng $80\text{ kg/m}^3$, việc sử dụng xi măng có hàm lượng khoáng $C_3S+C_3A$ cao như Quang Sơn ($81,70% + 7,37% = 89,07%$) sẽ sinh ra tổng nhiệt lượng thủy hóa lớn hơn $41,25%$ so với xi măng Nghi Sơn ($34,95% + 7,61% = 42,56%$). Đây là luận cứ khoa học giải thích nguyên nhân các công trình áp dụng cùng cấp phối định mức nhưng vẫn xuất hiện mức độ nứt nhiệt rất khác nhau.
Khung phân tích độc đáo
Khung phân tích của luận án kết hợp 3 lý thuyết nền tảng:
-
Phương trình vi phân dẫn nhiệt không ổn định tổng quát (Fourier - Laplace):
$$\frac{\partial T}{\partial t} = a \left( \frac{\partial^2 T}{\partial x^2} + \frac{\partial^2 T}{\partial y^2} + \frac{\partial^2 T}{\partial z^2} \right) + \frac{q_v}{C \rho}$$
Trong đó $a = \frac{\lambda}{C \rho}$ là hệ số khuếch tán nhiệt; $\lambda = 1,1 \div 2,3\text{ kcal/m}\cdot^\circ\text{C}$; $C = 0,22 \div 0,27\text{ kcal/kg}\cdot^\circ\text{C}$; $q_v$ là công suất nguồn nhiệt nội tại phụ thuộc vào đạo hàm thời gian của hàm nhiệt thủy hóa $Q(t)$.
-
Hệ điều kiện biên phi tuyến đa thông số:
- Biên đối lưu nhiệt (Loại 3): $dQ = \alpha (T_w - T_f) dt$ mô tả trao đổi nhiệt giữa bề mặt hở của bê tông với không khí hoặc nước.
- Biên bức xạ mặt trời: Gia số nhiệt độ bề mặt $T_{ow} = \frac{Q_0}{\beta_w}$, phản ánh bức xạ nhiệt tại công trình.
- Biên tiếp xúc lớp đổ động (Layer-by-Layer Boundary): Biên nhiệt thay đổi liên tục theo thời gian đổ lớp thứ $I$ ($T_D$) và thời gian nghỉ giữa các lớp ($T_N$).
- Mô hình đàn hồi nhiệt liên hợp (Thermo-mechanical Coupling Model):
Mối quan hệ ứng suất - biến dạng phẳng được xác định thông qua ma trận đàn hồi $[D]$ và véctơ tải trọng nút do biến thiên nhiệt độ:
$${\sigma} = [D]{\varepsilon} - [D]{\varepsilon_0}$$
$${F_3} = \sum_{e=1}^n \iint_V [B]^T [D] {\varepsilon_0} dV$$
với ${\varepsilon_0} = {\alpha \Delta T, \alpha \Delta T, 0}^T$, $\Delta T = \frac{T_i + T_j + T_k}{3} - T_{ref}$.
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ủ thế giới quan thực chứng (Positivism) với phương pháp định lượng mô phỏng số kết hợp kiểm chứng thực tế hiện trường. Thiết kế nghiên cứu đa tầng (Multi-level design) bao gồm:
- Tầng 1 (Vật liệu vi mô): Phân tích hóa học định lượng thành phần 4 khoáng chất chính ($C_3S, C_2S, C_3A, C_4AF$) của 12 mẫu xi măng Poóc-lăng thương phẩm theo tiêu chuẩn ASTM C150 / TCVN 6069:2007.
- Tầng 2 (Cơ lý vĩ mô theo thời gian): Mô hình hóa hàm cường độ chịu nén ($f_c$), cường độ kéo ($f_t$), mô-đun đàn hồi ($E_c$) và biến dạng kéo giới hạn ($\varepsilon_p$) ở các mốc 28 và 90 ngày tuổi:
$$f_c^{28} = 19,459(CKD/N) - 0,147(PGK/CKD) - 11,681 \quad (\text{MPa})$$
$$f_c^{90} = 19,326(CKD/N) - 0,333(PGK/CKD) + 5,968 \quad (\text{MPa})$$
$$E_c^{90} = 21,217(CKD/N) - 0,197(PGK/CKD) - 1,376 \quad (\text{GPa})$$
- Tầng 3 (Kết cấu công trình): Giải bài toán biến dạng phẳng 2 chiều trên mô hình mặt cắt đập thực tế chiều cao $H$, bề rộng đáy $B_D$, chiều cao cổ đập $H_D$, bề rộng chân đập thượng lưu $B_C$ đặt trên nền đá đàn hồi.
Quy trình nghiên cứu rigorous
Quy trình nghiên cứu được số hóa hoàn toàn bằng ngôn ngữ lập trình tham số APDL (ANSYS Parametric Design Language) theo sơ đồ thuật toán tự động hóa:
Data và phân tích
Mô hình phân tích sử dụng cặp phần tử tương đương hữu hạn:
- Phần tử nhiệt PLANE55: Phần tử 2D phẳng 4 điểm nút với bậc tự do duy nhất là nhiệt độ ($T$), có khả năng dẫn nhiệt đẳng hướng/bất đẳng hướng và mô phỏng nguồn nhiệt sinh bên trong phụ thuộc thời gian.
- Phần tử kết cấu PLANE182: Phần tử kết cấu 2D 4 điểm nút với 2 bậc tự do chuyển vị tại mỗi nút ($u_x, u_y$), hỗ trợ tính dẻo, từ biến, giãn nở nhiệt và siêu đàn hồi.
Dữ liệu đầu vào thực nghiệm cho 12 loại xi măng được chuẩn hóa với các chỉ số khoáng học:
| Chủng loại Xi măng |
$C_3S$ (%) |
$C_2S$ (%) |
$C_3A$ (%) |
$C_4AF$ (%) |
Tổng $C_3S+C_3A$ (%) |
Nhiệt phát sinh (J/g) |
| Nghi Sơn |
34,95 |
39,39 |
7,61 |
9,00 |
42,56 |
197 |
| La Hiên |
41,84 |
29,29 |
7,98 |
10,67 |
49,82 |
256 |
| Bút Sơn |
46,27 |
30,24 |
15,00 |
7,39 |
61,27 |
257 |
| Quán Triều |
48,21 |
27,61 |
7,10 |
11,16 |
55,31 |
259 |
| Tam Điệp |
50,43 |
27,56 |
7,64 |
10,18 |
58,07 |
266 |
| Hà Tiên |
52,78 |
18,48 |
10,00 |
8,00 |
62,78 |
286 |
| Xuân Thành |
51,57 |
25,84 |
6,97 |
8,72 |
58,54 |
292 |
| Điện Biên |
59,90 |
18,33 |
7,00 |
10,67 |
66,90 |
304 |
| Công Thanh |
62,92 |
16,18 |
6,14 |
9,45 |
69,06 |
319 |
| Tân Quang |
61,29 |
14,25 |
8,05 |
9,45 |
69,34 |
321 |
| Holcim |
65,27 |
7,65 |
6,69 |
10,61 |
71,96 |
355 |
| Quang Sơn |
81,70 |
0,26 |
7,37 |
10,67 |
89,07 |
420 |
Phát hiện đột phá và implications
Những phát hiện then chốt
-
Quy luật chi phối của khoáng học xi măng đến đỉnh nhiệt độ khối đập:
Sự chênh lệch hàm lượng khoáng sinh nhiệt $(C_3S + C_3A)$ giữa các nhà máy xi măng nội địa tạo ra sự chênh lệch đỉnh nhiệt độ thủy hóa trong thân đập lên tới $12,4^\circ\text{C}$ dưới cùng một điều kiện thi công. Cụ thể, xi măng Quang Sơn đạt nhiệt lượng $420\text{ J/g}$, trong khi xi măng Nghi Sơn chỉ sinh nhiệt $197\text{ J/g}$.
-
Hiệu ứng giải tỏa ứng suất nhiệt của Phụ gia khoáng (PGK):
Khi tăng tỷ lệ tro bay từ $30%$ lên $50%$ trong thành phần chất kết dính, nhiệt độ đoạn nhiệt của khối đổ BTĐL giảm từ $\Delta T = 25^\circ\text{C}$ xuống còn $\Delta T = 17^\circ\text{C}$, tương ứng mức giảm $32%$. Mức giảm này kéo theo sự suy giảm $38,5%$ giá trị ứng suất kéo cực đại tại vùng biên tiếp giáp nền đá ($\sigma_{max}$ giảm từ $2,18\text{ MPa}$ xuống $1,34\text{ MPa}$), đưa trạng thái ứng suất về dưới ngưỡng kháng kéo an toàn ($f_t \approx 1,5\div1,8\text{ MPa}$).
-
Hiện tượng nứt do chênh lệch nhiệt độ trong - ngoài ($\Delta T_{tn}$) và trên - dưới ($\Delta T_{td}$):
Nghiên cứu chứng minh rằng nứt bề mặt và nứt song song mặt đập không xuất hiện ở giai đoạn tăng nhiệt ban đầu mà phát triển mạnh nhất vào mùa lạnh đầu tiên sau khi đổ bê tông. Lúc này, nhiệt độ bên ngoài giảm mạnh làm bề mặt co lại trong khi lõi đập còn duy trì nhiệt độ cao, gây ra ứng suất kéo bề mặt vượt giới hạn kéo dãn cực hạn $\varepsilon_p = (0,85 \div 1,1) \times 10^{-4}$.
-
Kịch bản thi công tối ưu cho vùng khí hậu miền Trung và Bắc Bộ:
Kết quả mô phỏng số khẳng định đối với các đập có quy mô như Trung Sơn, việc khống chế nhiệt độ hỗn hợp BTĐL trước khi đổ ở mức $T_p \le 23^\circ\text{C}$, chiều dày mỗi lớp đổ sau đầm lèn đạt $0,30\text{ m}$, thi công liên tục 3 lớp trong 1 đợt đổ ($0,90\text{ m}$) kết hợp thời gian nghỉ giãn cách giữa các đợt từ $3\div5$ ngày là điều kiện biên động lực tối ưu để không vượt quá tiêu chuẩn nhiệt độ đỉnh $[T_M] \le 15\div20^\circ\text{C} + T_m + T_b$.
Implications đa chiều
- Về mặt lý thuyết: Bổ sung tham số khoáng học động $k_C$ vào phương trình phát nhiệt phi tuyến của vật liệu liên hợp polymer vô cơ / bê tông đầm lăn tro bay cao, mở rộng trường ứng dụng của cơ học phá hủy nhiệt trong bê tông khối lớn.
- Về phương pháp luận: Chuẩn hóa quy trình mô phỏng liên hợp Nhiệt - Cơ học bằng thuật toán APDL trên ANSYS, cho phép các kỹ sư tính toán chính xác ứng suất nhiệt công trình nhiều lớp đổ mà không cần đơn giản hóa điều kiện biên thành hằng số tĩnh.
- Về ứng dụng thực tiễn: Cung cấp giải pháp định lượng cụ thể cho các nhà thầu và ban quản lý dự án thủy điện trong việc lựa chọn thương hiệu xi măng (ưu tiên các dòng xi măng ít tỏa nhiệt có $C_3S < 50%, C_3A < 8%$) và tỷ lệ phối trộn tro bay/Puzơlan thích hợp cho từng vị trí khối đổ (vùng đáy kiềm chế mạnh cần PGK $>40%$).
- Về chính sách và quy chuẩn: Đặt cơ sở dữ liệu thực nghiệm để Bộ Xây dựng và Bộ Nông nghiệp & Phát triển Nông thôn xây dựng Tiêu chuẩn Quốc gia (TCVN) chuyên biệt về Thiết kế và Kiểm soát nhiệt cho Đập bê tông đầm lăn, thay thế việc lệ thuộc vào các tiêu chuẩn ngoại nhập SL314/SL319 của Trung Quốc.
Limitations và Future Research
Mặc dù đạt được những kết quả đột phá, luận án có các giới hạn học thuật cần ghi nhận:
- Giới hạn bài toán phẳng 2 chiều: Mô hình tính toán được thiết lập dưới giả thiết biến dạng phẳng (2D Plane Strain). Giả thiết này phản ánh tốt mặt cắt ngang giữa đập nhưng chưa mô phỏng trọn vẹn hiệu ứng kiềm chế 3 chiều (3D) tại các khối đổ áp sát hai vai đập địa hình hẻm vực chữ V.
- Giả thiết hằng số về hệ số dẫn nhiệt: Các thông số vật lý nhiệt như hệ số dẫn nhiệt ($\lambda$) và nhiệt dung riêng ($C$) được giả định là hằng số không thay đổi theo ngày tuổi đông kết, dù thực tế chúng có biến thiên nhỏ trong 72 giờ đầu.
- Cơ chế từ biến phi tuyến: Luận án chủ yếu xem xét ứng suất đàn hồi nhiệt tuyến tính có xét hệ số từ biến kinh nghiệm, chưa tích hợp đầy đủ mô hình lưu biến học phi tuyến (Nonlinear Viscoelastic-Plastic Creep) theo thời gian thực.
Định hướng nghiên cứu tiếp theo (Future Agenda):
- Phát triển mô hình mô phỏng không gian 3 chiều (3D FEA) tích hợp đầy đủ tương tác nhiệt - ứng suất - thủy lực (THM Coupling) có xét đến ảnh hưởng của hệ thống ống làm lạnh cưỡng bức đặt trong thân đập.
- Nghiên cứu cơ chế tự liền sẹo (self-healing) của các vi vết nứt nhiệt trong môi trường ngập nước của BTĐL sử dụng tro bay hoạt tính cao.
- Xây dựng thuật toán trí tuệ nhân tạo (Machine Learning/Deep Neural Networks) dự báo thời gian thực trường ứng suất nhiệt dựa trên dữ liệu chuỗi thời gian quan trắc từ các cảm biến sợi quang (Fiber Optic Sensors) gắn trong đập.
Tác động và ảnh hưởng
- Ảnh hưởng học thuật: Các công thức hiệu chỉnh $k_C$ và quy trình macro APDL mở ra hướng nghiên cứu mới về vật liệu nhiệt đầm lăn, ước tính tạo tiền đề cho hàng chục công bố khoa học chuyên ngành trên các tạp chí thuộc danh mục Scopus/Web of Science (như Construction and Building Materials, Materials and Structures).
- Chuyển đổi ngành công nghiệp xây dựng thủy điện: Giúp tối ưu hóa tỷ lệ cấp phối BTĐL tại các dự án thủy điện lớn, giảm thiểu chi phí xử lý sự cố nứt hàng chục tỷ đồng cho mỗi công trình, đồng thời rút ngắn tiến độ thi công từ $2\div4$ tháng nhờ khả năng kiểm soát đợt đổ an toàn.
- Tác động môi trường và xã hội: Tăng cường tỷ lệ sử dụng tro bay phế thải từ các nhà máy nhiệt điện than (như Phả Lại, Nghi Sơn, Vũng Áng) lên đến $40\div50%$ lượng chất kết dính, trực tiếp giảm phát thải khí nhà kính $\text{CO}_2$ từ công nghiệp sản xuất clinker xi măng, thúc đẩy kinh tế tuần hoàn trong xây dựng hạ tầng năng lượng.
Đối tượng hưởng lợi
- Nghiên cứu sinh & Giới học thuật: Tiếp cận phương pháp luận liên hợp nhiệt - cơ học chuẩn xác, sở hữu bộ số liệu khoáng học xi măng Việt Nam được phân tích hệ thống phục vụ các nghiên cứu mở rộng.
- Kỹ sư tư vấn thiết kế (TVTK): Nắm vững công cụ tính toán kiểm tra an toàn nhiệt theo quy chuẩn, tự tin đưa ra các thông số khe co dãn nhiệt, kích thước khối đổ và định mức vật liệu tối ưu.
- Nhà thầu thi công: Làm chủ công nghệ kiểm soát nhiệt độ hỗn hợp bê tông trước khi đổ ($T_p$), tốc độ nâng đập, và quy trình bảo dưỡng ẩm bề mặt ngăn ngừa nứt nẻ.
- Chủ đầu tư & Cơ quan Quản lý Nhà nước (EVN, Bộ NN&PTNT, Bộ Công Thương): Đảm bảo an toàn tuyệt đối cho các công trình đập dâng trọng lực cấp đặc biệt, kéo dài tuổi thọ vận hành công trình trên 100 nă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à đã mở rộng lý thuyết nào?
Đóng góp độc đáo nhất là việc thiết lập Hệ số ảnh hưởng khoáng học $k_C = 0,0073(C_3S + C_3A) + 0,4842$, mở rộng trực tiếp Lý thuyết động học nhiệt thủy hóa chất kết dính của Lưu Thụ Hoa (Lu Shuhua). Đóng góp này chuyển đổi mô hình tính toán nhiệt từ dạng "định mức đồng nhất" sang mô hình "động học thích ứng khoáng vật địa phương".
2. Điểm mới về mặt phương pháp luận so với các nghiên cứu trước đây là gì?
So với nghiên cứu của Barret & Tatro (1992) hay Lê Quốc Toàn (2014) vốn cố định điều kiện biên nhiệt môi trường, luận án đã xây dựng chu trình thuật toán APDL tự động hóa trên ANSYS mô phỏng biên nhiệt thay đổi động lực theo từng bước đổ (Layer-by-Layer), liên kết trực tiếp giữa phân tích truyền nhiệt không ổn định (phần tử PLANE55) sang phân tích ứng suất cơ học phẳng (phần tử PLANE182).
3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì?
Sự phân hóa cực lớn về khả năng sinh nhiệt giữa các dòng xi măng thương phẩm tại Việt Nam: Xi măng Quang Sơn có tổng lượng khoáng $C_3S+C_3A$ lên tới $89,07%$, sinh nhiệt $420\text{ J/g}$ (cao gấp 2,13 lần so với xi măng Nghi Sơn có tổng khoáng $42,56%$ và sinh nhiệt $197\text{ J/g}$). Điều này lý giải vì sao các tiêu chuẩn thiết kế cố định nhiệt độ đổ $T_p$ toàn quốc luôn thất bại trong việc ngăn ngừa nứt tại một số vùng miền.
Nhiệt thủy hóa (J/g)
4. Luận án có cung cấp quy trình tái lập (Replication Protocol) không?
Có. Luận án cung cấp chi tiết:
- Toàn bộ bảng thành phần khoáng và nhiệt năng 12 loại xi măng (Bảng 3.1).
- Các hệ phương trình cơ lý phụ thuộc tỷ lệ $N/CKD$ và $PGK/CKD$ (Công thức 1.1 đến 1.10).
- Cấu trúc sơ đồ khối thuật toán APDL, câu lệnh phần tử tương đương và các bước đặt tải biên nhiệt theo thời gian thực.
5. Lộ trình nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?
- Giai đoạn 1 (1-3 năm): Mở rộng hệ thống hóa dữ liệu khoáng học cho $100%$ các nhà máy xi măng trên toàn quốc và tích hợp mô-đun hóa 3D cho hẻm vực dốc.
- Giai đoạn 2 (4-6 năm): Nghiên cứu ứng dụng bê tông đầm lăn tro bay siêu mịn kết hợp sợi gia cường phân tán nhằm tăng khả năng chịu biến dạng kéo giới hạn $\varepsilon_p$.
- Giai đoạn 3 (7-10 năm): Chuyển giao công nghệ xây dựng phần mềm chuyên dụng tính toán nhiệt đập BTĐL mang thương hiệu Việt Nam (VN-RCC Thermal Solver) và hoàn thiện hệ thống Tiêu chuẩn Quốc gia TCVN.
Kết luận
- Luận án đã giải quyết thành công bài toán cấp bách về nứt do ứng suất nhiệt trong thi công đập BTĐL tại Việt Nam – vấn đề từng gây tổn hại nghiêm trọng đến độ bền và an toàn công trình thủy lợi - thủy điện.
- Thiết lập thành công mô hình toán học hiệu chỉnh nhiệt thủy hóa chất kết dính $Q = c \cdot k_C \cdot k_F \cdot Q_t$ có xét đến đặc thù khoáng học $C_3S + C_3A$ của xi măng nội địa với độ chính xác và độ tương quan rất cao ($R^2 = 0,9893$).
- Xây dựng hoàn chỉnh bộ chương trình tính toán số tự động hóa bằng ngôn ngữ APDL trong môi trường ANSYS, giải quyết trọn vẹn bài toán truyền nhiệt phi tuyến liên hợp ứng suất phẳng có xét đến sự thay đổi điều kiện biên và tiến độ phân lớp thi công thực tế.
- Xác định định lượng vai trò của phụ gia khoáng tro bay (tỷ lệ tối ưu $30\div50%$) trong việc hạ thấp nhiệt độ đoạn nhiệt $\Delta T$ xuống $17^\circ\text{C}$ và giảm ứng suất kéo tức thời xuống dưới ngưỡng nứt cho phép.
- Đề xuất bộ giải pháp công nghệ đồng bộ từ khâu lựa chọn chủng loại xi măng, thiết kế cấp phối, kiểm soát nhiệt độ đổ ($T_p \le 23^\circ\text{C}$), đến quy trình phân đợt đổ 3 lớp dày $0,3\text{ m}$ phù hợp cho từng vùng khí hậu đặc thù, đã được ứng dụng thử nghiệm hiệu quả tại đập thủy điện Trung Sơn.
- Kết quả nghiên cứu là nền tảng khoa học vững chắc để chuẩn hóa tiêu chuẩn thiết kế và thi công đập bê tông đầm lăn tại Việt Nam, đóng góp tích cực vào sự phát triển bền vững của cơ sở hạ tầng năng lượng quốc gia.
