Tổng quan về luận án

Công nghệ bê tông đầm lăn (Roller Compacted Concrete - RCC hay BTĐL) là bước đột phá kỹ thuật quan trọng nhất trong lĩnh vực xây dựng đập thủy lợi - thủy điện hơn bốn thập kỷ qua. Với hơn 300 đập BTĐL được xây dựng tại 44 quốc gia và hơn 30 công trình đập lớn đã hoàn thành tại Việt Nam (tiêu biểu như Pleikrông 71m, Định Bình 52,3m, Bản Vẽ 137m, Sơn La 138,1m, Lai Châu 137m), kết cấu đập trọng lực BTĐL đóng vai trò huyết mạch trong an ninh năng lượng và quản trị tài nguyên nước quốc gia. Mặc dù sở hữu các ưu thế vượt trội về tốc độ thi công nhanh, giảm lượng xi măng và hạ giá thành từ 25% đến 40% so với bê tông truyền thống (CVC), đập BTĐL đối mặt với thách thức lớn về an toàn công trình do hiện tượng nứt thân đập. Theo báo cáo của Hội đập lớn thế giới (ICOLD), nứt là thuộc tính phổ biến của các kết cấu bê tông khối lớn; như chuyên gia Zhu Bo Fang từng nhận định: "không có đập bê tông sẽ không có nứt".

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực tế: các công trình khoa học trong nước và quốc tế trước đây chủ yếu tập trung kiểm soát nhiệt thủy hóa, tối ưu hóa cấp phối tro bay/puzolan hoặc phân tích nhiệt - ứng suất trong giai đoạn thi công (Lê Anh Vân & Nguyễn Văn Mao, 2006; Đỗ Văn Lượng, 2008; Nguyễn Như Oanh, 2011, 2013). Các nghiên cứu mô hình hóa toán học về sự hình thành, mở rộng vết nứt và cơ chế suy giảm khả năng chịu tải của đập trọng lực BTĐL trong giai đoạn vận hành dưới tổ hợp tải trọng tĩnh (áp lực nước, áp lực đẩy nổi, bùn cát) và tải trọng động chấn động cực hạn (MCE, OBE) vẫn còn phân tán, thiếu công cụ tính toán tự thích ứng có độ tin cậy cao. Tài liệu thực tiễn chỉ ra rằng "nguyên nhân gây ra nứt kết cấu công trình thuộc về biến dạng là chính và chiếm tới 80%, do tải trọng là chính chiếm khoảng 20%". Tuy nhiên, trong thời kỳ khai thác, sự kết hợp giữa vết nứt khởi sinh sẵn có và tải trọng thủy lực - động đất là mối nguy cơ trực tiếp gây mất ổn định tổng thể.

Luận án xác lập 3 câu hỏi nghiên cứu và hệ giả thuyết khoa học:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để loại bỏ hiện tượng kỳ dị ứng suất tại mũi vết nứt và vùng chân đập trong phương pháp phần tử hữu hạn (PTHH) đàn hồi tuyến tính truyền thống?
    • Giả thuyết 1 (H1): Ứng dụng kỹ thuật chia lưới tự thích ứng (Adaptive Mesh Refinement - mô hình h) kết hợp chuẩn sai số năng lượng Zienkiewicz-Zhu (ZZ) và phương pháp Ứng suất tương đương (ƯSTĐ-PTHH) sẽ chuẩn hóa trường ứng suất cục bộ và tiệm cận chính xác nghiệm giải tích sức bền vật liệu.
  • Câu hỏi nghiên cứu 2 (RQ2): Quy luật phát triển của các dạng vết nứt điển hình (vết nứt ngang chân thượng lưu, vết nứt đứng xuyên nền, vết nứt dọc trục thân đập và đa vết nứt đồng thời) diễn biến như thế nào dưới tác động của động đất?
    • Giả thuyết 2 (H2): Tải trọng động làm tăng đột biến ứng suất kéo chính tại đỉnh vết nứt, kích hoạt cơ chế phá hủy hỗn hợp (Mode I và Mode II), trong đó vết nứt ngang chân thượng lưu có mức độ nguy hiểm cao nhất đối với ổn định trượt và lật của đập.
  • Câu hỏi nghiên cứu 3 (RQ3): Các vết nứt thi công thực tế tại đập thủy điện Sơn La có đe dọa trực tiếp đến an toàn chịu tải của công trình trong điều kiện vận hành cực đoan hay không?
    • Giả thuyết 3 (H3): Vết nứt ban đầu tại khối đổ C2 và C8 đập Sơn La (chiều dài 31,5m, độ sâu 6m, bề rộng ~1mm ghi nhận tháng 06/2008) chỉ gây phá hoại cục bộ hữu hạn và đập vẫn bảo đảm an toàn chịu lực tổng thể dưới động đất thiết kế cực đại.

Khung lý thuyết của nghiên cứu được tích hợp từ Cơ học rạn nứt phá hủy (Fracture Mechanics), Lý thuyết cơ học đàn - dẻo và hư hỏng phi tuyến của bê tông (Concrete Damage Plasticity / Smeared Crack Model), cùng Cơ học tính toán nâng cao (Computational Adaptive FEM). Phạm vi nghiên cứu tập trung vào các mặt cắt đại diện của đập trọng lực BTĐL ở Việt Nam, kiểm chứng qua các bộ dữ liệu thực nghiệm quốc tế (đập Koyna - Ấn Độ, đập Guandi - Trung Quốc) và áp dụng trực tiếp vào công trình đập thủy điện Sơn La (cao 138,1m).

Literature Review và Positioning

Lịch sử phát triển công nghệ BTĐL khởi đầu từ các nghiên cứu tiên phong của Jerome Raphael (1970) với báo cáo "Đập trọng lực tối ưu", tiếp nối bởi Cannon R.W. (1972-1974) và Dunstan M.R.H. (1976-1979) tại Mỹ và Anh, mở đường cho việc xây dựng đập Willow Creek (1980) và Upper Stillwater (1985). Từ đây hình thành ba trường phái công nghệ chính: (1) Trường phái Mỹ (USBR, USACE) sử dụng BTĐL nghèo xi măng (<100 kg/m³) kết hợp màng chống thấm mặt thượng lưu; (2) Trường phái Nhật Bản (RCD - Roller Compacted Dam) sử dụng kết cấu "vàng bọc bạc" với lớp vỏ bê tông thường (CVC) dày bao bọc lõi BTĐL; (3) Trường phái Trung Quốc (RCCD) thi công đầm lăn toàn mặt cắt với bê tông biến thái (GEVR) và BTĐL cấp phối II giàu chất kết dính (130-190 kg/m³) làm khối chống thấm chính.

Tổng quan y văn chỉ ra các cuộc tranh luận học thuật sâu sắc về nguyên nhân và cơ chế phá hoại do nứt trong đập BTĐL:

  1. Tranh luận về kiểm soát nứt vật liệu và thi công: Nhóm nghiên cứu của Zhu Bo Fang (1999) và Nguyễn Minh Việt (2016) nhấn mạnh nứt chủ yếu do ứng suất nhiệt phát sinh từ nhiệt thủy hóa chậm và chênh lệch nhiệt độ môi trường (đặc biệt tại các đập Guanyinge, Jiangya, Yushi ở Trung Quốc). Ngược lại, Nguyễn Chiến (2014) và Nguyễn Tiến Oanh phân tích sự cộng hưởng phức tạp giữa co khô bề mặt phơi lộ, xung nhiệt mùa đông và khoảng cách bố trí khe lún/khe co giãn ngang (khoảng cách 31,5m tại Sơn La bị xem là quá lớn so với khuyến cáo 16-20m).
  2. Tranh luận về phương pháp mô hình hóa cơ học rạn nứt: L.J. Malvar và L.J. Warren (1990) khởi xướng phương pháp PTHH với lưới cố định, nhưng bộc lộ nhược điểm phụ thuộc mạnh vào kích thước phần tử và không mô tả chính xác sự dịch chuyển vết nứt. Các nghiên cứu hiện đại của Mridha S. (2014) trên phần mềm ABAQUS, Gaohui Wang (2017) với mô hình tương tác 3D đập - nước - nền, Amir Behshad (2018), Shou-yan Jiang (2012) và Deeja Alora (2014) áp dụng phần tử tiếp xúc ma sát và hệ lò xo - cản để mô phỏng vết nứt dưới tác động động đất tại đập Koyna. Tuy nhiên, việc tái tạo động lực học mạng lưới phần tử tự thích ứng (adaptive re-meshing) đồng thời xử lý triệt để bài toán sóng phản xạ vô hạn ở biên nền vẫn còn nhiều hạn chế.

Tại Việt Nam, các tác giả Đinh Hữu Dụng & Nguyễn Văn Mao (2012), Lê Quốc Toản, Vũ Thanh Tễ & Vũ Hoàng Hưng (2015) đã thiết lập các mô-đun ANSYS để phân tích trường nhiệt - ứng suất thi công. Luận án định vị khoảng trống học thuật bằng việc chuyển tiếp trọng tâm từ "kiểm soát nhiệt thi công" sang "đánh giá sự lan truyền vết nứt và khả năng chịu tải trong vận hành", thiết lập chương trình tính toán số tích hợp ngôn ngữ APDL trên nền ANSYS với lưới tự thích ứng mô hình h và biên đàn hồi - cản nhớt, đối sánh trực tiếp với 2 ca nghiên cứu kinh điển thế giới là đập Koyna (Ấn Độ) và đập Guandi (Trung Quố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 Cơ học Phá hủy Đàn dẻo Tuyến tính và Phi tuyến (Linear/Nonlinear Fracture Mechanics) áp dụng cho vật liệu bê tông đầm lăn thông qua các đóng góp cụ thể:

  1. Hoàn thiện Tiêu chuẩn phá hoại bê tông nhiều trục Willam-Warnke: Tích hợp phương trình mặt phá hủy 5 tham số của Willam-Warnke trong phần tử SOLID65 của ANSYS, phản ánh chính xác hành vi phi tuyến của BTĐL bao gồm phá hủy do nứt kéo (khi ứng suất chính $\sigma_1 \ge f_t$) và vỡ vụn do nén (crushing khi ứng suất nén vượt $f_c$). Tác giả lượng hóa hệ số truyền lực cắt qua bề mặt vết nứt mở ($\beta_t$) và vết nứt đóng ($\beta_c$), giải thích cơ chế dẻo dai tăng dần theo tuổi của BTĐL (tỷ lệ cường độ kéo/nén đạt 0,10 - 0,14 sau 90 ngày).
  2. Xây dựng Lý thuyết Lưới tự thích ứng dựa trên Chuẩn sai số năng lượng Zienkiewicz-Zhu (ZZ): Luận án chứng minh sự vượt trội của kỹ thuật làm mịn lưới tự động (mô hình h) trong việc xóa bỏ hiện tượng giả tạo ứng suất vô hạn tại các góc đổi hướng hình học và đỉnh vết nứt.
  3. Chuyển dịch hệ hình (Paradigm shift) trong đánh giá an toàn đập: Chuyển đổi từ quan điểm truyền thống "chấp nhận trạng thái ứng suất đàn hồi thuần túy dựa trên ứng suất cho phép" sang hệ hình "đánh giá dung sai vết nứt động lực học" (Dynamic Damage Tolerance Assessment), cho phép định lượng khả năng chịu tải dư của kết cấu bê tông khối lớn khi đã tồn tại khuyết tật hình học.
+-----------------------------------------------------------------------------------+
|                        KHUNG PHÂN TÍCH TỔNG QUAN LUẬN ÁN                          |
+-----------------------------------------------------------------------------------+
                                          |
    +-------------------------------------+-----------------------------------+
    |                                     |                                   |
    v                                     v                                   v
[LÝ THUYẾT NỀN TẢNG]             [PHƯƠNG PHÁP TÍNH TOÁN SỐ]         [ĐỐI TƯỢNG VÀ KIỂM CHỨNG]
- Cơ học phá hủy bê tông         - PTHH tự thích ứng (AFEM)          - Đập Koyna (Ấn Độ - 1967)
- Mặt phá hủy Willam-Warnke      - Chuẩn sai số năng lượng ZZ        - Đập Guandi (Trung Quốc)
- Mô hình Smeared Crack          - Phương pháp ƯSTĐ-PTHH             - Đập thủy điện Sơn La
- Tương tác Đập - Nước - Nền     - Ngôn ngữ lập trình APDL           (Khối C2, C8 - Cao 138,1m)
    |                                     |                                   |
    +-------------------------------------+-----------------------------------+
                                          |
                                          v
                [CHƯƠNG TRÌNH TÍNH TOÁN CHUYÊN DỤNG RCCD_FAILURE]
                - Sinh lưới tự động h-refinement theo sai số năng lượng
                - Xử lý biên nhân tạo đàn hồi - cản nhớt vô hạn
                - Mô phỏng lan truyền nứt dưới tải trọng tĩnh & động chấn
                                          |
                                          v
                [KẾT QUẢ ĐÁNH GIÁ KHẢ NĂNG CHỊU TẢI & ỔN ĐỊNH CÔNG TRÌNH]

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

Khung phân tích của luận án tích hợp 3 trụ cột lý thuyết: (1) Cơ học tính toán tự thích ứng; (2) Cơ học rạn nứt phi tuyến công trình thủy; (3) Động lực học công trình tương tác kết cấu - chất lưu - nền móng vô hạn.

Điểm độc đáo nằm ở việc phát triển Phương pháp Ứng suất tương đương Phần tử hữu hạn tự thích ứng (ƯSTĐ-PTHH). Phương pháp này thực hiện trung bình hóa trường ứng suất phức tạp tại mặt cắt đáy khối đập thu được từ PTHH tự thích ứng theo nguyên lý tích phân nội lực của lý thuyết Sức bền vật liệu (SBVL):

$$\sigma_{Y}(x) = \frac{N}{F} \pm \frac{M}{I_x} \cdot x$$

Trong đó $N$ là tổng lực dọc, $M$ là mômen uốn tổng cộng tại mặt cắt, $F$ là diện tích tiết diện và $I_x$ là mômen quán tính. Giải pháp hỗn hợp này loại bỏ hoàn toàn hiện tượng tập trung ứng suất ảo tại mép biên tiếp giáp nền - bê tông, cung cấp biểu đồ phân bố ứng suất pháp thực chất, phục vụ kiểm tra điều kiện bền theo các quy chuẩn thiết kế đập hiện hành.

Điều kiện biên nghiên cứu được xác lập chặt chẽ: áp dụng điều kiện biên nhân tạo đàn hồi - cản nhớt (Viscoelastic artificial boundary) dựa trên hệ thống lò xo - giảm chấn phân bố liên tục tại đáy và các mặt bên của mô hình nền, triệt tiêu hoàn toàn sự phản xạ giả tạo của sóng địa chấn ngược trở lại thân đập.

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

Thiết kế nghiên cứu

Nghiên cứu theo đuổi quan điểm thực chứng luận (Positivism) kết hợp phương pháp mô phỏng số tất định (Deterministic numerical simulation) và phân tích đối chuẩn thực nghiệm. Thiết kế nghiên cứu đa tầng bao gồm:

  • Cấp độ vật liệu (Micro/Meso level): Mô hình hóa quan hệ ứng suất - biến dạng phi tuyến của BTĐL, mô đun đàn hồi $E$, hệ số Poisson $\nu$, cường độ kháng nén đa trục $f_c$, cường độ kháng kéo uốn $f_t$, hệ số truyền lực cắt $\beta_t$ và $\beta_c$.
  • Cấp độ phần tử (Element level): Ứng dụng phần tử không gian 2D 8-nút (PLANE82 / PLANE183) và phần tử khối 3D SOLID65 chuyên dụng cho bê tông có khả năng tái tạo vết nứt và biến dạng dẻo trong ANSYS.
  • Cấp độ kết cấu (Macro/Structural level): Mô hình hóa tổng thể mặt cắt ngang đập trọng lực BTĐL chiều cao lớn, gắn kết với hệ thống nền đá vô hạn và cột nước thượng lưu thông qua áp lực thủy động học Westergaard cải tiến.

Quy trình nghiên cứu rigorous

Quy trình tính toán số tự thích ứng được tự động hóa hoàn toàn thông qua mô-đun lập trình tham số RCCD_FAILURE bằng ngôn ngữ APDL, tuân thủ các bước:

[BƯỚC 1: Khởi tạo mô hình hình học & Chia lưới phần tử ban đầu (Lưới thô)]
                                |
                                v
[BƯỚC 2: Gán điều kiện biên đàn hồi - cản nhớt & Đặt tổ hợp tải trọng]
                                |
                                v
[BƯỚC 3: Giải bài toán PTHH phi tuyến đàn dẻo / động lực học thời gian thực]
                                |
                                v
[BƯỚC 4: Tính toán sai số chuyển vị e_u, sai số ứng suất e_sigma]
                                |
                                v
[BƯỚC 5: Đánh giá Sai số Năng lượng phần tử ||e||_E,i & Toàn miền ||e||_E]
                                |
                                v
                < Kiểm tra tiêu chuẩn sai số: eta <= eta_cho_phep? >
                               /                  \
                             YES                   NO
                             /                      \
                            v                        v
            [BƯỚC 7: Xuất trường nội lực,    [BƯỚC 6: Tự động mịn hóa lưới cục bộ
             phát hiện nứt & chuyển bước]      (h-refinement) tại vùng sai số cao]
                                                     \
                                                      +---> (Quay lại BƯỚC 3)

Thuật toán đánh giá sai số năng lượng của phần tử thứ $i$ và sai số toàn miền được xác định chuẩn xác theo công thức Zienkiewicz-Zhu:

$$|e|{E, i}^2 = \int{\Omega_i} \left({\sigma^} - {\hat{\sigma}}\right)^T [D]^{-1} \left({\sigma^} - {\hat{\sigma}}\right) d\Omega$$

$$\eta = \frac{|e|_E}{\sqrt{|\hat{e}|_E^2 + |e|_E^2}} \times 100% \le \bar{\eta}$$

Trong đó ${\hat{\sigma}}$ là trường ứng suất giải gần đúng, ${\sigma^*}$ là trường ứng suất làm mịn cải tiến thu được qua phép nội suy liên tục, $[D]$ là ma trận độ cứng đàn hồi của vật liệu, và $\bar{\eta}$ là ngưỡng sai số cho phép định trước. Khi $\eta > \bar{\eta}$, chương trình tự động gia tăng số lượng phần tử $N_{new}$ tại các miền tập trung ứng suất cao trước khi tiến hành bước lặp tải trọng tiếp theo.

Data và phân tích

Bộ dữ liệu kiểm chứng và thông số tính toán được chuẩn hóa từ các công trình thực tế:

  • Bộ thông số vật liệu đập Sơn La: Khối lượng riêng $\rho = 2,4 , \text{T/m}^3$; Mô đun đàn hồi bê tông $E = 2,1 \times 10^7 , \text{kN/m}^2$; Hệ số Poisson $\nu = 0,18$; Cường độ chịu nén tính toán $f_c = 15.000 , \text{kN/m}^2$; Cường độ chịu kéo $f_t = 1.500 , \text{kN/m}^2$; Hệ số truyền lực cắt khi nứt mở $\beta_t = 0,2$, khi nứt đóng $\beta_c = 0,8$. Nền đá bazan: $E_n = 1,8 \times 10^7 , \text{kN/m}^2$, $\nu_n = 0,22$.
  • Dữ liệu kích động địa chấn: Sử dụng băng gia tốc thực tế của trận động đất Koyna (11/1967, gia tốc đỉnh $0,49g$) và các băng gia tốc số nhân tạo tương thích với phổ phản ứng thiết kế cực đại MCE ($a_{max} = 0,158g$) và phổ động đất cơ sở vận hành OBE ($a_{max} = 0,079g$) phê duyệt riêng cho tuyến thủy điện Sơn La.
  • Kiểm chứng độ tin cậy phần mềm (Model Validation): Mô-đun RCCD_FAILURE được đối chuẩn với kết quả mô phỏng ABAQUS và thí nghiệm mô hình vật lý của Mridha S. (2014) trên đập Koyna, cũng như kết quả tính toán 3D của Gaohui Wang (2017) và đập Guandi (Trung Quốc). Hình thức phá hủy nứt gãy uốn ngang tại vị trí thay đổi độ dốc mái hạ lưu đập Koyna và phân bố nứt trượt khớp nối đập Guandi hoàn toàn trùng khớp với dữ liệu thực nghiệm đã công bố, khẳng định độ chính xác tuyệt đối của thuật toán.

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

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

Qua hàng loạt kịch bản tính toán mô phỏng số, luận án xác lập 4 phát hiện khoa học đột phá có giá trị thực chứng cao:

  1. Cơ chế tập trung ứng suất và lan truyền vết nứt ngang chân thượng lưu: Dưới tác dụng của tổ hợp tải trọng cơ bản (trọng lượng bản thân, áp lực nước ở mực nước dâng bình thường MNDBT, áp lực đẩy nổi), ứng suất kéo chính tại góc chân thượng lưu đạt cực đại. Khi tồn tại vết nứt ngang ban đầu, vết nứt chỉ mở rộng cục bộ trong phạm vi nhỏ (<5% chiều rộng đáy đập). Tuy nhiên, khi chịu tổ hợp tải trọng đặc biệt có động đất Koyna hoặc MCE, ứng suất kéo chính vượt ngưỡng $f_t$, vết nứt phát triển xiên chúc xuống nền đá rồi bẻ góc nằm ngang hướng về phía hạ lưu, làm suy giảm nghiêm trọng diện tích chịu nén hiệu dụng của mặt cắt đáy.
  2. Quy luật phân nhánh của vết nứt đứng xuyên từ nền: Dưới tải trọng tĩnh, vết nứt đứng từ nền bị khống chế bởi ứng suất nén thẳng đứng do trọng lượng bản thân thân đập. Khi gia tốc đỉnh động đất tăng gấp 2 lần, 3 lần và 4 lần, vết nứt đứng phát triển vọt lên cao trong thân đập, sau đó xuất hiện hiện tượng phân nhánh (bifurcation), hình thành vết nứt ngang thứ cấp tại chân thượng lưu và làm tăng biên độ chuyển vị ngang đỉnh đập lên gấp 3,2 lần so với trạng thái đập nguyên vẹn.
  3. Hiện tượng che chắn ứng suất (Stress Shielding) khi đa vết nứt đồng thời xuất hiện: Khi trong thân đập đồng thời tồn tại vết nứt ngang chân thượng lưu, vết nứt đứng từ nền và vết nứt dọc trục khối đổ, trường ứng suất bị tái phân bố sâu sắc. Vết nứt ngang thượng lưu đóng vai trò kích hoạt sự giải phóng năng lượng biến dạng, làm chậm tốc độ phát triển của vết nứt đứng nhưng lại gia tăng chuyển vị đỉnh đập và làm tăng nguy cơ phá hủy uốn cục bộ tại các vị trí tiếp giáp hành lang kiểm tra.
  4. Đánh giá an toàn thực tế đập thủy điện Sơn La: Mô phỏng số chi tiết cho khối C2 và C8 đập Sơn La với vết nứt phát hiện tháng 06/2008 (dài 31,5m, sâu 6m tại cao trình +169,0m) khẳng định:
    • Dưới tổ hợp tải trọng cơ bản (MNDBT = +215m, mực nước hạ lưu cực tiểu), vùng nứt không có xu hướng phát triển thêm, trạng thái ứng suất toàn đập hoàn toàn nằm trong giới hạn cho phép của vật liệu ($|\sigma_{nén}| \le 8,5 , \text{MPa} < f_c$, $|\sigma_{kéo}| \le 0,8 , \text{MPa} < f_t$).
    • Dưới tác động của trận động đất MCE cực đoan, tại vị trí vết nứt ban đầu chỉ xuất hiện vùng phá hoại biến dạng dẻo cục bộ, không hình thành vết nứt xuyên cắt đứt thân đập. Chuyển vị đỉnh đập theo phương ngang dao động ổn định trong biên độ cực đại $\le 18,4 , \text{mm}$, khẳng định kết cấu đập BTĐL Sơn La vẫn bảo đảm khả năng chịu tải và độ ổn định tổng thể.

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp minh chứng định lượng cho lý thuyết tương tác nứt - kết cấu đập khối lớn, mở rộng hiểu biết về cơ chế lan truyền nứt động lực học trong môi trường vật liệu dị tính phân lớp của BTĐL.
  • Về mặt phương pháp luận: Thiết lập quy trình chuẩn tích hợp giữa phần mềm thương mại ANSYS, ngôn ngữ APDL và phương pháp làm mịn lưới tự thích ứng (AFEM), cung cấp một công cụ mô phỏng số mạnh mẽ có thể áp dụng cho các dạng kết cấu thủy công phức tạp khác (đập vòm, đập tràn ngầm, cửa lấy nước).
  • Về mặt thực tiễn và chính sách: Cung cấp cơ sở khoa học then chốt để Hội đồng Nghiệm thu Nhà nước và Tập đoàn Điện lực Việt Nam (EVN) tự tin quyết định tích nước và đưa công trình thủy điện Sơn La vào vận hành an toàn vượt tiến độ 3 năm, tiết kiệm hàng chục nghìn tỷ đồng cho ngân sách. Đồng thời, kết quả nghiên cứu đặt nền tảng cho việc biên soạn Tiêu chuẩn thiết kế và Hướng dẫn quan trắc đánh giá vết nứt đập BTĐL tại Việt Nam.

Limitations và Future Research

Nghiên cứu thẳng thắn thừa nhận các giới hạn học thuật nội tại:

  1. Mô hình hóa hình học 2D chiếm ưu thế: Mặc dù đã kiểm chứng 3D trên mô hình đập Koyna và Guandi, các phân tích lan truyền nứt chi tiết tại đập Sơn La chủ yếu thực hiện trên mô hình bài toán biến dạng phẳng (Plane Strain) cho mặt cắt đại diện, chưa phản ánh đầy đủ hiệu ứng không gian 3D tương tác giữa các phân đoạn đập dọc trục thung lũng hẹp.
  2. Đơn giản hóa áp lực thủy động học bên trong vết nứt (Hydro-fracture coupling): Áp lực nước thấm thâm nhập trực tiếp vào bên trong khe nứt khi vết nứt mở ra trong quá trình động đất mới được mô phỏng dưới dạng tải trọng phân bố tĩnh tương đương, chưa giải quyết trọn vẹn bài toán ghép cặp đầy đủ chất lưu - chất rắn (Full Hydro-Mechanical FSI Coupling).
  3. Mô hình vật liệu đồng nhất hóa theo lớp: Chưa xét đến sự suy thoái độ bám dính theo thời gian tại các mặt tiếp giáp thi công (lift joints) do tác động xâm thực hóa học hoặc phản ứng kiềm - cốt liệu (ASR).

Chương trình nghiên cứu 10 năm tiếp theo (Future Research Agenda) cần tập trung:

  • Phát triển mô hình lan truyền nứt 3D tích hợp Phương pháp Phần tử Hữu hạn Mở rộng (XFEM) và Cơ học Hư hỏng Liên tục (Phase-Field Fracture).
  • Nghiên cứu ghép cặp Thủy lực - Nhiệt - Cơ học (THM) mô tả sự thâm nhập áp lực nước áp lực cao vào mũi vết nứt chuyển động nhanh.
  • Tích hợp mạng nơ-ron học máy (Physics-Informed Neural Networks - PINNs) với dữ liệu quan trắc cảm biến sợi quang thời gian thực để chẩn đoán và dự báo sớm nguy cơ mất an toàn đập.

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

  • Tác động học thuật (Academic Impact): Tạo ra bước tiến phương pháp luận trong tính toán rạn nứt kết cấu thủy công tại Việt Nam; cung cấp tài liệu tham khảo nền tảng cho các nghiên cứu sau đại học, có tiềm năng trích dẫn cao trong các tạp chí chuyên ngành như Journal of Structural Engineering, Computers & Structures, và Tạp chí Khoa học & Công nghệ Thủy lợi.
  • Chuyển đổi ngành công nghiệp Thủy điện (Industry Impact): Ứng dụng trực tiếp cho công tác vận hành, bảo trì và kiểm định an toàn đập tại hơn 30 nhà máy thủy điện BTĐL trên toàn quốc, trực tiếp bảo vệ an toàn cho các công trình trọng điểm quốc gia như Sơn La (2.400 MW), Lai Châu (1.200 MW), Bản Vẽ (320 MW).
  • Lợi ích xã hội và môi trường: Bảo đảm an toàn tuyệt đối cho vùng hạ lưu sông Hồng và đồng bằng Bắc Bộ, nơi sinh sống của hơn 20 triệu người dân, loại bỏ nguy cơ thảm họa vỡ đập dây chuyền.

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

  • Nghiên cứu sinh & Nhà khoa học chuyên ngành Thủy công: Tiếp cận quy trình lập trình APDL và phương pháp phần tử hữu hạn tự thích ứng tiên tiến để giải quyết các bài toán phi tuyến phức tạp.
  • Kỹ sư Tư vấn Thiết kế (PECC1, PECC2, HEC): Ứng dụng phương pháp ƯSTĐ-PTHH để chuẩn hóa tính toán kiểm tra ứng suất và thiết kế gia cố khe co giãn, khe dẫn nứt đập trọng lực.
  • Cơ quan Quản lý Nhà nước & Chủ đập (EVN, Bộ Nông nghiệp & PTNT, Bộ Công Thương): Sở hữu cơ sở khoa học xác đáng để ban hành quy trình đánh giá rủi ro, vận hành hồ chứa an toàn và xây dựng kịch bản ứng phó khẩn cấp động đất.

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 Phương pháp Ứng suất tương đương Phần tử hữu hạn tự thích ứng (ƯSTĐ-PTHH) kết hợp tiêu chuẩn phá hủy bê tông Willam-Warnke. Phương pháp này giải quyết triệt để nghịch lý tập trung ứng suất vô hạn tại góc lượn biên và chân đập của phương pháp PTHH đàn hồi truyền thống bằng cách kết hợp ưu điểm của PTHH tự thích ứng (mô hình h theo chuẩn sai số ZZ) với lý thuyết tích phân nội lực Sức bền vật liệu, tạo ra lời giải chuẩn xác phục vụ thiết kế công trình.

2. Đột phá phương pháp luận so với các nghiên cứu trước đây thể hiện thế nào?
So với mô hình lưới cố định của L.J. Malvar & L.J. Warren (1990) vốn gây phụ thuộc kích thước lưới, và mô hình tiếp xúc đơn giản của Shou-yan Jiang (2012) hay Deeja Alora (2014), luận án xây dựng mô-đun RCCD_FAILURE trên ANSYS APDL tự động tái tạo lưới theo năng lượng sai số, đồng thời tích hợp biên đàn hồi - cản nhớt vô hạn loại bỏ phản xạ sóng địa chấn, đạt độ tương thích cao với kết quả thực nghiệm đập Koyna của Mridha S. (2014) và Gaohui Wang (2017).

3. Phát hiện bất ngờ nhất từ dữ liệu nghiên cứu là gì?
Dưới tác động của tải trọng động chấn gia cường (PGA tăng 2 đến 4 lần), vết nứt đứng xuyên từ nền không phát triển thẳng đứng lên mặt tràn đập mà bị uốn cong và phân nhánh (bifurcation) do tương tác với trường ứng suất trọng lực bản thân, kích hoạt vết nứt ngang thứ cấp tại chân thượng lưu trước khi xảy ra phá hủy dẻo tổng thể.

4. Nghiên cứu có cung cấp giao thức tái lập (Replication Protocol) không?
Có. Luận án cung cấp toàn bộ quy trình thuật toán APDL, các phương trình xác định sai số năng lượng Zienkiewicz-Zhu, thông số vật liệu phi tuyến đập và nền, cùng các bước thiết lập phần tử lò xo - cản nhớt trên ANSYS, cho phép các nhà nghiên cứu độc lập tái lập và kiểm chứng hoàn toàn kết quả tính toán.

5. Chương trình nghị sự nghiên cứu 10 năm tới được phác thảo như thế nào?
Mở rộng phương pháp sang mô hình nứt 3D bằng phương pháp Phase-Field / XFEM, giải quyết trọn vẹn bài toán tương tác thủy - nhiệt - cơ học (THM) khi nước áp lực cao thâm nhập vào khe nứt chuyển động, và tích hợp AI/PINNs với hệ thống quan trắc thời gian thực phục vụ bảo trì dự đoán đập thông minh.

Kết luận

Luận án của NCS Hà Thanh Dương đã hoàn thành xuất sắc các mục tiêu nghiên cứu với 5 đóng góp học thuật và thực tiễn cốt lõi:

  1. Xây dựng thành công mô-đun tính toán tự động RCCD_FAILURE bằng ngôn ngữ lập trình tham số APDL trên nền tảng ANSYS, ứng dụng kỹ thuật phần tử hữu hạn tự thích ứng mô hình h kiểm soát bởi chuẩn sai số năng lượng Zienkiewicz-Zhu.
  2. Đề xuất và chứng minh tính ưu việt của Phương pháp Ứng suất tương đương (ƯSTĐ-PTHH), kết nối hoàn hảo giữa phương pháp PTHH hiện đại và lý thuyết Sức bền vật liệu truyền thống, loại bỏ nghiệm kỳ dị ứng suất cục bộ.
  3. Làm sáng tỏ cơ chế phát sinh, lan truyền và phân nhánh của 4 dạng vết nứt điển hình trong thân đập BTĐL dưới tác động của tổ hợp tải trọng tĩnh và động đất gia cường.
  4. Kiểm chứng thành công thuật toán số qua hai chuẩn đối sánh quốc tế kinh điển là đập Koyna (Ấn Độ) và đập Guandi (Trung Quốc), khẳng định độ tin cậy và tính khái quát cao của mô hình.
  5. Định lượng chính xác trạng thái chịu lực và khẳng định độ an toàn ổn định tuyệt đối của công trình đập thủy điện lớn nhất Đông Nam Á - đập BTĐL Sơn La - khi tồn tại các vết nứt thi công ban đầu dưới các kịch bản động đất thiết kế cực đại.

Công trình tạo tiền đề mở ra 3 hướng nghiên cứu mới về cơ học phá hủy đập bê tông khối lớn, để lại giá trị khoa học lâu dài và đóng góp thiết thực cho sự nghiệp phát triển bền vững hạ tầng thủy năng quốc gia.