Tổng quan về luận án

Công nghệ bê tông đầm lăn (Roller Compacted Concrete - RCC) đại diện cho bước chuyển biến mang tính cách mạng trong lĩnh vực xây dựng công trình thủy lợi và thủy điện toàn cầu suốt hơn bốn thập kỷ qua. Theo định nghĩa của Viện Bê tông Hoa Kỳ (ACI), bê tông đầm lăn là "loại bê tông được đầm bằng máy đầm lăn; hỗn hợp bê tông ở dạng chưa đông cứng có khả năng hỗ trợ cho máy đầm khi hoạt động đầm". Tính đến nay, toàn thế giới đã xây dựng trên 300 đập BTĐL tại hơn 44 quốc gia, dẫn đầu là Trung Quốc, Hoa Kỳ, Nhật Bản và Tây Ban Nha. Tại Việt Nam, công nghệ này đã được ứng dụng xây dựng trên 30 đập bê tông trọng lực quy mô lớn (chiều cao trên 60 m), tiêu biểu như đập thủy điện Sơn La (cao 138,1 m, thể tích BTĐL đạt 3,1 triệu $\text{m}^3$ trong tổng số 4,6 triệu $\text{m}^3$ bê tông toàn đập, xếp thứ 10 về chiều cao và thứ 3 về thể tích BTĐL trên thế giới tại thời điểm xây dựng), Bản Vẽ (cao 137 m), Lai Châu (cao 137 m), và Đồng Nai 4 (cao 129 m).

                      ┌────────────────────────────────────────┐
                      │    ĐẬP TRỌNG LỰC BÊ TÔNG ĐẦM LĂN       │
                      │         (RCC GRAVITY DAM)              │
                      └──────────────────┬─────────────────────┘
                                         │
                 ┌───────────────────────┴───────────────────────┐
                 ▼                                               ▼
  ┌───────────────────────────────┐               ┌───────────────────────────────┐
  │     TRƯỜNG ỨNG SUẤT TĨNH      │               │     TẢI TRỌNG ĐỘNG ĐẤT        │
  │ (Trọng lượng + Thủy tĩnh TH1) │               │   (Phổ thiết kế MDE / OBE)    │
  └──────────────┬────────────────┘               └──────────────┬────────────────┘
                 │                                               │
                 └───────────────────────┬───────────────────────┘
                                         │
                                         ▼
                      ┌────────────────────────────────────────┐
                      │     KHỞI SINH VÀ PHÁT TRIỂN NỨT        │
                      │  (Chân thượng lưu, tiếp giáp móng,     │
                      │   khe thi công & thân đập C2, C3)      │
                      └──────────────────┬─────────────────────┘
                                         │
                                         ▼
                      ┌────────────────────────────────────────┐
                      │        MÔ HÌNH HÓA SỐ NÂNG CAO         │
                      │  - AFEM: Lưới tự thích ứng (h-adapt)   │
                      │  - ƯSTĐ-PTHH: Khử kỳ dị ứng suất       │
                      │  - Biên đàn hồi - cản nhớt vô hạn      │
                      │  - Tiêu chuẩn phá hoại Willam-Warnke   │
                      └──────────────────┬─────────────────────┘
                                         │
                                         ▼
                      ┌────────────────────────────────────────┐
                      │    ĐÁNH GIÁ KHẢ NĂNG CHỊU TẢI &        │
                      │     ỔN ĐỊNH CÔNG TRÌNH VẬN HÀNH        │
                      └────────────────────────────────────────┘

Mặc dù mang lại hiệu quả kinh tế vượt trội (giảm 25% - 40% chi phí xây dựng và rút ngắn thời gian thi công), hiện tượng nứt kết cấu đập BTĐL diễn ra phổ biến. Theo thống kê của Hội đập lớn thế giới (ICOLD), phần lớn các đập bê tông xây dựng trên thế giới đều xuất hiện các dạng nứt khác nhau; thậm chí nhà khoa học Zhu Bo Fang đã đưa ra nhận định kinh điển: "Không có đập bê tông sẽ không có nứt". Các nghiên cứu trước đây tại Việt Nam và quốc tế chủ yếu tập trung vào bài toán nhiệt thủy hóa, tối ưu hóa cấp phối tro bay hoặc kiểm soát nứt do co ngót trong giai đoạn thi công (Lê Anh Vân, 2006; Đỗ Văn Lượng, 2008; Đinh Hữu Dụng và Nguyễn Văn Mạo, 2012; Nguyễn Như Oanh, 2013). Khoảng trống học thuật cốt lõi (research gap) chính là sự thiếu vắng một công cụ tính toán mô phỏng hoàn chỉnh quá trình khởi sinh, lan truyền vết nứt và định lượng mức độ suy giảm khả năng chịu tải tổng thể của đập trọng lực BTĐL trong giai đoạn vận hành dưới tác động đồng thời của tải trọng thủy tĩnh và động đất cực đoan.

Luận án của nghiên cứu sinh Hà Thanh Dương (2019) tại Trường Đại học Thủy Lợi, dưới sự hướng dẫn khoa học của PGS. Nguyễn Quang Hùng và GS. Vũ Thanh Te, đã giải quyết tường minh khoảng trống này thông qua hệ thống câu hỏi và giả thuyết nghiên cứu:

  • Câu hỏi nghiên cứu 1 (RQ1): Cơ chế tập trung ứng suất cục bộ tại các điểm kỳ dị hình học (chân đập thượng lưu, vị trí thay đổi độ dốc) và các mặt tiếp giáp khối đổ phát triển như thế nào trong điều kiện vận hành?
  • Câu hỏi nghiên cứu 2 (RQ2): Làm thế nào để mô phỏng chính xác quá trình lan truyền vết nứt phi tuyến trong môi trường BTĐL mà không bị phụ thuộc vào sự phân chia mạng lưới phần tử cố định?
  • Câu hỏi nghiên cứu 3 (RQ3): Sự hiện diện của các vết nứt ban đầu phát sinh trong quá trình thi công (như tại khối C2, C3 của đập thủy điện Sơn La) ảnh hưởng ở mức độ nào đến độ bền và an toàn chịu tải của đập dưới các cấp kích động địa chấn thiết kế?

Hệ thống giả thuyết tương ứng:

  • Giả thuyết 1 (H1): Phương pháp phần tử hữu hạn tự thích ứng (AFEM) kết hợp đánh giá sai số chuẩn năng lượng cho phép hội tụ chính xác quỹ đạo nứt thực tế của đập.
  • Giả thuyết 2 (H2): Phương pháp Ứng suất tương đương (ƯSTĐ-PTHH) giải quyết triệt để sự phi lý của hiện tượng kỳ dị ứng suất tại góc chân đập, phản ánh đúng khả năng chịu lực của vật liệu.
  • Giả thuyết 3 (H3): Các vết nứt thi công cục bộ không phát triển xuyên suốt phá hoại thân đập dưới tổ hợp tải trọng cơ bản và động đất cơ sở vận hành (OBE), công trình vẫn bảo đảm điều kiện an toàn chịu tải.

Nghiên cứu được định vị trên khung lý thuyết cơ học rạn nứt phá hủy (Fracture Mechanics), lý thuyết phá hoại phi tuyến vật liệu bê tông đa trục (Willam-Warnke), và phương pháp phần tử hữu hạn tự thích ứng phát triển bằng ngôn ngữ APDL (ANSYS Parametric Design Language). Phạm vi nghiên cứu bao gồm mô hình lý thuyết tổng quát, kiểm chứng số học trên dữ liệu thực nghiệm đập Koyna (Ấn Độ), đập Guandi (Trung Quốc), và ứng dụng thực chứng định lượng trên đập thủy điện Sơn La.

Literature Review và Positioning

Nghiên cứu về đập bê tông đầm lăn trên thế giới đã hình thành ba trường phái công nghệ và thiết kế chính:

  1. Trường phái Hoa Kỳ (RCC - Roller Compacted Concrete): Khởi xướng từ các nghiên cứu của Raphael (1970), Cannon (1972) và Dunstan (Anh, 1970s), phát triển bởi USACE và USBR. Đặc trưng là sử dụng BTĐL nghèo xi măng (hàm lượng chất kết dính $< 100\text{ kg/m}^3$), sử dụng lớp màng chống thấm PVC hoặc tấm panen bê tông cốt thép mặt thượng lưu (điển hình tại đập Willow Creek hoàn thành năm 1980 cao 52 m, đập Upper Stillwater cao 90 m).
  2. Trường phái Nhật Bản (RCD - Roller Compacted Dam): Khởi đầu từ đập Shimajigawa (1980, cao 89 m), áp dụng kết cấu "vàng bọc bạc" với lớp bê tông thường (CVC) mác cao bao quanh khối lõi BTĐL nhằm đáp ứng yêu cầu chống thấm và chống nứt bề mặt.
  3. Trường phái Trung Quốc (RCCD - Roller Compacted Concrete Dams): Tích hợp từ kinh nghiệm của Mỹ và Nhật, đầm nén trên toàn mặt cắt, sử dụng BTĐL cấp phối II giàu chất kết dính ($180 - 220\text{ kg/m}^3$) kết hợp bê tông biến thái (GEVR) làm lớp chống thấm thượng lưu (áp dụng tại đập Longtan cao 216,5 m, đập Guanyinge cao 82 m, đập Jiangya cao 131 m, đập Yushi cao 50,2 m).
Trường phái / Nghiên cứu Đặc điểm công nghệ & Kết cấu Phương pháp phân tích nứt Giới hạn học thuật
Trường phái Hoa Kỳ (USBR, USACE, Raphael 1970) BTĐL nghèo xi măng ($< 100\text{ kg/m}^3$), chống thấm bằng màng PVC hoặc panen đúc sẵn. Cơ học đàn hồi tuyến tính, Sức bền vật liệu (SBVL). Bỏ qua tập trung ứng suất đáy móng; xuất hiện nhiều vết nứt thấm (đập Galesville, Elk Creek).
Trường phái Nhật Bản (RCD Method, đập Shimajigawa) Kết cấu "vàng bọc bạc", lớp vỏ CVC mác cao bảo vệ lõi BTĐL. Phân tích phần tử hữu hạn (PTHH) lưới cố định. Phức tạp trong thi công phân tầng; lưới cố định gây phụ thuộc nghiệm khi mô phỏng nứt.
Trường phái Trung Quốc (RCCD, Zhu Bo Fang 1999) Đầm nén toàn mặt cắt, BTĐL cấp phối II và bê tông biến thái GEVR. Mô hình hóa ứng suất nhiệt thi công, phụ gia MgO bù co ngót. Chủ yếu tập trung giai đoạn thi công; chưa giải quyết tối ưu tương tác nứt động lực học vận hành.
Nghiên cứu Quốc tế (Mridha S 2014, Gaohui Wang 2017) Mô hình thí nghiệm động lực học và mô phỏng số ABAQUS đập Koyna. Mô hình tương tác 3D đập - nước - nền, phần tử tiếp xúc nứt. Khối lượng tính toán cực lớn; chưa tích hợp thuật toán tái tạo lưới tự thích ứng động.
Luận án (Hà Thanh Dương, 2019) Toàn mặt cắt BTĐL nhiệt đới hóa (công trình đập Sơn La, Đồng Nai 4). Tích hợp AFEM (mô hình h) + ƯSTĐ-PTHH + Lập trình APDL ANSYS. Giải quyết triệt để điểm kỳ dị ứng suất, tự động thích ứng lưới theo quỹ đạo phát triển vết nứt.

Trong y văn tồn tại hai luồng tranh luận học thuật sâu sắc:

  • Tranh luận 1 về bản chất nứt đập: Zhu Bo Fang và các nhà nghiên cứu Trung Quốc chỉ ra rằng nguyên nhân gây nứt do biến dạng (nhiệt, co ngót, lún không đều) chiếm tới 80%, trong khi do tải trọng trực tiếp chỉ chiếm khoảng 20%. Tuy nhiên, các nhà cơ học phá hủy quốc tế (như Malvar, Amir Behshad 2018) chứng minh rằng khi đập bước vào giai đoạn tích nước và chịu tải trọng động đất, chính tổ hợp ứng suất thủy tĩnh và sóng chấn động sẽ kích hoạt các vết nứt biến dạng tiềm ẩn phát triển thành vết nứt phá hoại kết cấu.
  • Tranh luận 2 về công cụ tính toán ứng suất: Phương pháp Sức bền vật liệu (SBVL) dựa trên giả thiết phân bố ứng suất tuyến tính $y = \frac{N}{L} \pm \frac{6M}{L^2}$ bỏ qua hoàn toàn sự tập trung ứng suất tại chân đập, dẫn đến đánh giá quá lạc quan về an toàn. Ngược lại, Phương pháp Phần tử hữu hạn (PTHH) đàn hồi truyền thống lại tạo ra giá trị ứng suất kéo vô hạn (điểm kỳ dị toán học) tại góc tiếp giáp nền - đập, khiến kỹ sư không thể đưa ra tiêu chuẩn giới hạn ứng suất vật lý ($f_t$) để đánh giá an toàn.

Luận án định vị chính xác vào điểm giao thoa này bằng cách so sánh đối chiếu trực tiếp với hai nghiên cứu quốc tế chuẩn mực: mô hình mô phỏng số trên ABAQUS kết hợp thí nghiệm mô hình của Mridha S. (2014) và mô hình tương tác 3D đập - nước - nền dưới động đất Koyna của Gaohui Wang (2017). Luận án đã vượt lên giới hạn của lưới cố định nhờ thuật toán sinh lưới tự thích ứng, tái hiện chính xác đường nứt gãy ngang thân đập Koyna và phân bố nứt đập Guandi với độ tin cậy vượt trội.

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

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

Luận án đã đóng góp vào hệ thống lý thuyết cơ học công trình thủy ba nội dung học thuật quan trọng:

  1. Mở rộng lý thuyết lan truyền vết nứt phi tuyến cho môi trường BTĐL phân lớp: Bằng việc tích hợp tiêu chuẩn bền Willam-Warnke với mô hình vết nứt phân tán (smeared crack model) có xét đến hệ số truyền lực cắt khi vết nứt đóng ($\alpha$) và khi vết nứt mở ($\beta$), nghiên cứu đã mô tả chính xác quy luật suy giảm độ cứng và trượt cắt dọc theo các mặt tiếp giáp lớp đầm nén (lift joints).
  2. Xây dựng cơ sở lý thuyết Phương pháp Ứng suất tương đương (ƯSTĐ-PTHH): Thiết lập nguyên lý trung bình hóa trường ứng suất cục bộ tại các miền kỳ dị hình học thông qua tích phân phân bố lực cắt và mô men uốn trên mặt cắt kiểm tra, kết hợp hài hòa ưu điểm phản ánh biến dạng tổng thể của PTHH với tính thực tế của lý thuyết SBVL.
  3. Chứng minh cơ chế kiềm chế nứt thủy tĩnh: Đưa ra luận cứ lý thuyết về "hiệu ứng tự khống chế nứt" ở chân thượng lưu đập trọng lực, giải thích tại sao vết nứt khởi sinh do tải trọng tĩnh chỉ phát triển đến một độ sâu giới hạn rồi dừng lại do đi vào vùng ứng suất nén chính của thân đập.
                           ┌────────────────────────────────────────┐
                           │      LƯỚI PHẦN TỬ BAN ĐẦU (BƯỚC 0)     │
                           └──────────────────┬─────────────────────┘
                                              │
                                              ▼
                           ┌────────────────────────────────────────┐
                           │   GIẢI HỆ PHƯƠNG TRÌNH ĐỘ CỨNG PTHH    │
                           │           [K]{u} = {P}                 │
                           └──────────────────┬─────────────────────┘
                                              │
                                              ▼
                           ┌────────────────────────────────────────┐
                           │    ĐÁNH GIÁ SAI SỐ CHUẨN NĂNG LƯỢNG    │
                           │   e_u = ∫ (σ - σ*)^T D^-1 (σ - σ*) dΩ  │
                           └──────────────────┬─────────────────────┘
                                              │
                                              ▼
                                   ┌────────────────────┐
                        Không      │   Sai số toàn miền │
                     ┌─────────────┤      η ≤ f ?       │
                     │             └──────────┬─────────┘
                     │                        │ Có
                     ▼                        ▼
      ┌─────────────────────────────┐   ┌─────────────────────────────┐
      │   TỰ ĐỘNG CHIA NHỎ PHẦN TỬ  │   │   XUẤT KẾT QUẢ ỨNG SUẤT &   │
      │  (Kỹ thuật h-refinement tại │   │     VÙNG NỨT HỘI TỤ CHUẨN   │
      │    vùng nứt & chân đập)     │   └─────────────────────────────┘
      └──────────────┬──────────────┘
                     │
                     └────────────────────────┘ (Lặp bước tính)

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

Khung phân tích của luận án được xây dựng dựa trên sự hợp nhất của ba trụ cột phương pháp luận:

$$\text{Khung Phân Tích} = \text{Adaptive FEM (Mô hình } h\text{)} \oplus \text{Lý thuyết ƯSTĐ-PTHH} \oplus \text{Mô hình Tương tác Động Đập - Nền vô hạn}$$

  • Tích hợp đa lý thuyết: Kết hợp giải thuật Zienkiewicz-Zhu đánh giá sai số năng lượng với mô hình liên tục đàn - dẻo phá hủy.
  • Tiếp cận phân tích mới: Thay vì sử dụng phần tử tiếp xúc giao diện định trước (vốn đòi hỏi phải biết trước vị trí nứt như các nghiên cứu của Deeja Alora, 2014 hay Shou-yan Jiang, 2012), chương trình tự động dò tìm vị trí đạt giới hạn ứng suất kéo chính $\sigma_1 \ge f_t$ để khởi tạo vết nứt và tự động tái cấu trúc mạng lưới phần tử hữu hạn tại đỉnh vết nứt.
  • Điều kiện biên xác định nghiêm ngặt: Sử dụng hệ biên đàn hồi - cản nhớt (viscous-spring artificial boundary) đặt tại đáy và hai bên biên nền nhân tạo nhằm triệt tiêu hoàn toàn hiện tượng phản xạ sóng địa chấn ngược trở lại công trình, mô phỏng chính xác bán không gian vô hạn của nền đá.

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 định lượng (Quantitative Positivism) kết hợp phương pháp mô phỏng số học giải tích bậc cao. Thiết kế nghiên cứu được tổ chức đa cấp độ:

  • Cấp độ 1 (Vật liệu): Mô hình hóa đặc tính cơ học phi tuyến của BTĐL, phản ánh sự khác biệt giữa mô đun đàn hồi nén $E_c$ và kéo $E_t$, tỷ lệ Poisson $\nu$, và đặc tính dẻo hóa sau nứt (tension softening).
  • Cấp độ 2 (Phần tử & Kết cấu): Mô hình hóa bài toán biến dạng phẳng 2D mặt cắt đập thực tế chiều cao $H = 138,1\text{ m}$, bao gồm tương tác đập - nền đá có kích thước đủ lớn ($L_{nen} \ge 3H$, $H_{nen} \ge 2H$).
  • Cấp độ 3 (Động lực học hệ thống): Phân tích phản ứng động phi tuyến theo lịch sử thời gian (Time-history analysis) dưới tác động của các băng gia tốc động đất thực và nhân tạo.

Quy trình nghiên cứu rigorous

Quy trình tính toán tự thích ứng được tự động hóa hoàn toàn bằng chương trình chuyên dụng RCCD_FAILURE viết trên nền ngôn ngữ tham số APDL của ANSYS:

                                QUY TRÌNH TÍNH TOÁN CỦA CHƯƠNG TRÌNH RCCD_FAILURE
  ┌──────────────────────────────────────────────────────────────────────────────────────────────────────────┐
  │ 1. Khởi tạo mô hình hình học tham số hóa (Hình học thân đập, nền móng, vị trí nứt ban đầu).             │
  ├──────────────────────────────────────────────────────────────────────────────────────────────────────────┤
  │ 2. Thiết lập điều kiện biên cản nhớt và gán tải trọng: Trọng lượng bản thân, Áp lực nước (MNDBT/MNDGC).   │
  ├──────────────────────────────────────────────────────────────────────────────────────────────────────────┤
  │ 3. Giải bước tĩnh (Tổ hợp cơ bản TH1) ──> Đánh giá sai số năng lượng e_u và chuẩn hóa lưới h-adaptive.   │
  ├──────────────────────────────────────────────────────────────────────────────────────────────────────────┤
  │ 4. Áp đặt gia tốc động đất theo bước thời gian Δt ──> Giải phương trình động học phi tuyến.              │
  ├──────────────────────────────────────────────────────────────────────────────────────────────────────────┤
  │ 5. Kiểm tra điều kiện nứt tại từng phần tử: σ_1 ≥ f_t.                                                   │
  │    ├── NẾU ĐẠT: Giảm độ cứng phần tử, giải phóng liên kết, kích hoạt hệ số truyền lực cắt (α, β).       │
  │    └── TỰ ĐỘNG LÀM MỊN LƯỚI tại lân cận đỉnh vết nứt (Adaptive refinement).                              │
  ├──────────────────────────────────────────────────────────────────────────────────────────────────────────┤
  │ 6. Cập nhật ma trận độ cứng [K_t] ──> Tiếp tục bước thời gian t + Δt cho đến khi kết thúc băng sóng.    │
  └──────────────────────────────────────────────────────────────────────────────────────────────────────────┘

Công thức toán học kiểm soát sai số năng lượng toàn phần $\eta$:

$$e_u = \left( \int_{\Omega} (\sigma^* - \hat{\sigma})^T D^{-1} (\sigma^* - \hat{\sigma}) d\Omega \right)^{1/2}$$

$$\eta = \frac{|e_u|}{\left(|\hat{u}|^2 + |e_u|^2\right)^{1/2}} \le f$$

Trong đó: $\hat{\sigma}, \hat{u}$ là trường ứng suất và chuyển vị giải gần đúng; $\sigma^*$ là trường ứng suất cải tiến sau khi làm trơn; $D$ là ma trận hằng số đàn hồi; $f$ là ngưỡng sai số cho phép định mức ($f = 5%$).

Data và phân tích

Thông số đầu vào của mô hình được trích xuất trực tiếp từ hồ sơ thiết kế kỹ thuật và kết quả thí nghiệm cơ lý đập thủy điện Sơn La và đập Đồng Nai 4:

  • Vật liệu BTĐL thân đập: Khối lượng riêng $\rho = 2,45\text{ T/m}^3$; Mô đun đàn hồi $E = 2,4 \times 10^7\text{ kN/m}^2$ ($24\text{ GPa}$); Hệ số Poisson $\nu = 0,18$; Cường độ chịu nén tính toán $f_c = 18,5\text{ MPa}$; Cường độ chịu kéo tính toán $f_t = 1,45\text{ MPa}$; Hệ số truyền lực cắt vết nứt mở $\beta = 0,2$, vết nứt đóng $\alpha = 0,8$.
  • Vật liệu nền đá (Granit/Gneiss phong hóa nhẹ): Khối lượng riêng $\rho_r = 2,65\text{ T/m}^3$; Mô đun biến dạng $E_r = 1,8 \times 10^7\text{ kN/m}^2$ ($18\text{ GPa}$); Hệ số Poisson $\nu_r = 0,22$.
  • Bộ số liệu kích động động đất: Băng gia tốc số 1a chuẩn hóa theo Phổ động đất thiết kế cực đại (MDE/CME) với gia tốc nền cực đại $a_{max} = 0,228g$, chu kỳ kích động $T = 15\text{ s}$; Băng gia tốc số 4a chuẩn hóa theo Phổ động đất cơ sở vận hành (OBE) với $a_{max} = 0,114g$.

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

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

  1. Quy luật tự giới hạn của vết nứt chân đập dưới tải trọng tĩnh: 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), vết nứt ngang xuất hiện tại chân thượng lưu đập chỉ phát triển cục bộ trong phạm vi chiều dài từ $1,2\text{ m}$ đến $2,5\text{ m}$ (chiếm dưới $2%$ bề rộng đáy đập) rồi dừng hẳn. Nguyên nhân do khi tiến sâu vào thân đập, đỉnh vết nứt đi vào vùng ứng suất nén chính $\sigma_3$ do mô men trọng lực bản thân gây ra, triệt tiêu hoàn toàn ứng suất kéo $\sigma_1$.
  2. Cơ chế lan truyền nứt động dưới động đất cực đoan: Khi chịu tác động của động đất phổ MDE ($a_{max} = 0,228g$), tại thời điểm $t = 3,42\text{ s} - 4,15\text{ s}$, vết nứt từ chân thượng lưu tái phát triển mạnh mẽ và đổi hướng xiên chếch xuống phía hạ lưu với góc nghiêng $20^\circ - 32^\circ$. Đồng thời, xuất hiện vùng nứt thứ hai tại vị trí thay đổi độ dốc thân đập phía hạ lưu (vùng cổ đập) do hiệu ứng tập trung ứng suất động uốn.
  3. Định lượng mức độ an toàn của vết nứt thi công tại đập Sơn La: Phân tích cụ thể hai vết nứt thực tế phát hiện trong quá trình thi công tháng 06/2008 tại khối C2 và C3 (chiều dài vết nứt $31,5\text{ m}$, độ sâu $6\text{ m}$, bề rộng khe hở $1,0\text{ mm}$ do hiệu ứng co khô và ứng suất nhiệt khi bóc lộ qua mùa đông) chứng minh rằng: Dưới tác dụng của tải trọng vận hành bình thường và động đất cấp OBE, chuyển vị ngang tại đỉnh đập chỉ đạt cực đại $12,4\text{ mm}$, độ mở rộng thêm của vết nứt ban đầu là $0,00\text{ mm}$ (không phát triển). Đập hoàn toàn bảo đảm ổn định và khả năng chịu tải tổng thể.
  4. Ngưỡng phá hoại sụp đổ động đất của đập BTĐL: Khi nâng gia tốc đỉnh nền động đất lên gấp 3 đến 4 lần giá trị thiết kế ($0,684g - 0,912g$), vết nứt chân đập và vết nứt thân đập phát triển xuyên suốt cắt qua toàn bộ phân khối đổ, dẫn đến hình thành mặt trượt liên tục làm mất ổn định lật và trượt của phần đỉnh đập.
  5. Hiệu chỉnh sai số ứng suất móng: Phương pháp ƯSTĐ-PTHH giúp loại bỏ đỉnh tập trung ứng suất giả tạo tại chân đập, giảm giá trị ứng suất kéo tính toán từ mức phi thực tế $8,6\text{ MPa}$ (theo PTHH lưới mịn cố định) về mức ứng suất trung bình thực chất $1,28\text{ MPa} < f_t$, giải thích thỏa đáng vì sao nhiều đập thực tế không bị nứt toác đáy dù tính toán PTHH cổ điển cảnh báo vượt ứng suất kéo cho phép.
┌─────────────────────────────────────────────────────────────────────────────────────────────────────────┐
│              SƠ ĐỒ PHÂN VÙNG VÀ CƠ CHẾ PHÁT TRIỂN NỨT TRÊN MẶT CẮT ĐẬP BTĐL SƠN LA                     │
│                                                                                                         │
│                              ▼ Cao trình đỉnh +215.0m                                                   │
│                             ┌──┐                                                                        │
│                             │  │                                                                        │
│                             │  │  Chuyển vị đỉnh cực đại: 12.4 mm (OBE)                                 │
│                             │  │                                                                        │
│   Mực nước dâng             │  │◄──────── Vùng nứt cổ đập hạ lưu                                        │
│   bình thường (MNDBT)       │   \         (Xuất hiện dưới động đất cực đoan MDE)                        │
│   ▼ Cao trình +215.0m       │    \                                                                      │
│  ~~~~~~~~~~~~~~~~           │     \                                                                     │
│  │               │          │      \                                                                    │
│  │ Áp lực nước   │          │  C3   \ ◄── Vết nứt khối C3 (Thi công 2008: L=31.5m, d=6m)                │
│  │ tam giác      │          │        \    (Ổn định, không phát triển dưới tải OBE)                      │
│  │               │          │  C2     \                                                                 │
│  │               │          │          \                                                                │
│  ▼               ▼          │           \                                                               │
│ ────────────────────────────┴────────────\─────────────────────────── Cao trình đáy móng +77.0m         │
│  ▲ Vết nứt ngang chân TL                   \                                                            │
│    (Phát triển 1.2 - 2.5m rồi dừng)         \                                                           │
│ ═══════════════════════════════════════════════════════════════════════════════════════════════          │
│                 NỀN ĐÁ BÁN VÔ HẠN (Mô hình biên cản nhớt Viscous-Spring)                                │
└─────────────────────────────────────────────────────────────────────────────────────────────────────────┘

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp phương pháp luận mới về việc đánh giá an toàn đập dựa trên lý thuyết phá hủy thay thế cho phương pháp trạng thái ứng suất đàn hồi giới hạn thuần túy.
  • Về mặt phương pháp: Module lập trình tự động hóa RCCD_FAILURE trong ANSYS APDL có thể chuyển giao và ứng dụng trực tiếp cho mọi loại kết cấu đập bê tông khối lớn có hình học phức tạp.
  • Về mặt thực tiễn kỹ thuật: Đưa ra kết luận có cơ sở khoa học khẳng định đập thủy điện Sơn La an toàn trong quá trình vận hành dù tồn tại các vết nứt thi công chưa xử lý triệt để, giải tỏa các quan ngại kỹ thuật của Hội đồng Nghiệm thu Nhà nước.
  • Về chính sách và tiêu chuẩn: Cung cấp căn cứ thực tế để kiến nghị sửa đổi tiêu chuẩn thiết kế đập BTĐL Việt Nam (TCVN 9137:2012 và TCVN 4253:2012), cho phép áp dụng phương pháp ƯSTĐ-PTHH để tính toán ổn định và độ bền thay cho phương pháp giải tích SBVL đơn giản hóa.

Limitations và Future Research

Nhằm duy trì tính khách quan khoa học, luận án thẳng thắn thừa nhận các giới hạn nội tại:

  1. Áp lực thủy động trong khe nứt: Mô hình chưa xét đến áp lực nước xâm nhập động tức thời (hydrodynamic crack tip pressure) bên trong lòng vết nứt khi khe nứt đóng mở tuần hoàn theo nhịp sóng động đất.
  2. Giới hạn mô hình không gian: Phân tích chủ yếu thực hiện trên bài toán biến dạng phẳng 2D, chưa phản ánh đầy đủ hiệu ứng kiềm chế không gian ba chiều (3D spatial arching effect) của hẻm vực lòng sông hẹp.
  3. Đặc tính vật liệu theo thời gian: Chưa tích hợp quá trình suy thoái mỏi vật liệu (fatigue damage accumulation) và hiện tượng lão hóa bê tông sau 50 - 100 năm vận hành.

Chương trình nghiên cứu tiếp nối trong giai đoạn 10 năm tới bao gồm:

  • Hướng 1: Phát triển mô hình đa trường ghép Thủy - Cơ - Nhiệt (Thermo-Hydro-Mechanical - THM coupling) 3D cho toàn bộ hệ thống đập - hồ chứa - đới đứt gãy địa chất.
  • Hướng 2: Ứng dụng phương pháp Phần tử hữu hạn mở rộng (XFEM) hoặc Phương pháp Phần tử biên (BEM) để mô phỏng sự phân nhánh ngẫu nhiên của vết nứt dưới tải trọng va chạm động năng lớn.
  • Hướng 3: Tích hợp dữ liệu quan trắc thời gian thực từ hệ thống cảm biến đo nứt (Crack meters) và con lắc đo chuyển vị tại hiện trường vào mô hình phần tử hữu hạn tự thích ứng theo thuật toán Digital Twin (bản sao số).

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

  • Tác động học thuật: Mở ra hướng nghiên cứu chuyên sâu về đánh giá an toàn công trình thủy công theo độ tin cậy rạn nứt phá hủy; là tài liệu tham khảo nền tảng cho các chương trình đào tạo sau đại học chuyên ngành Kỹ thuật xây dựng công trình thủy (Mã số: 9.58.02.02).
  • Tác động công nghiệp: Ứng dụng trực tiếp cho Tập đoàn Điện lực Việt Nam (EVN) và các ban quản lý dự án thủy điện trong công tác kiểm định an toàn đập định kỳ, đánh giá rủi ro sự cố xả lũ và động đất kích thích cho chuỗi bậc thang thủy điện sông Đà (Hòa Bình, Sơn La, Lai Châu) và hệ thống thủy điện miền Trung - Tây Nguyên.
  • Tác động xã hội và an toàn cộng đồng: Bảo đảm an toàn tuyệt đối cho vùng hạ lưu sông Hồng với hàng chục triệu dân sinh sống, hạn chế các quyết định gia cố tốn kém hàng trăm tỷ đồng không cần thiết nhờ kết quả đánh giá vết nứt chính xác.

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

                                    CÁC NHÓM ĐỐI TƯỢNG HƯỞNG LỢI
  ┌───────────────────────────────┬──────────────────────────────────────────────────────────────────────────┐
  │ Nghiên cứu sinh & Nhà khoa học│ Tiếp cận thuật toán AFEM và mã nguồn APDL phục vụ nghiên cứu rạn nứt     │
  │ chuyên ngành Thủy công        │ phi tuyến công trình khối lớn.                                           │
  ├───────────────────────────────┼──────────────────────────────────────────────────────────────────────────┤
  │ Kỹ sư tư vấn & Doanh nghiệp   │ Sở hữu công cụ tính toán ƯSTĐ-PTHH chính xác, tối ưu hóa kích thước mặt  │
  │ thiết kế (PECC1, PECC2, VNK)  │ cắt đập và bố trí khe co giãn hợp lý, tiết kiệm 15-20% chi phí xử lý.    │
  ├───────────────────────────────┼──────────────────────────────────────────────────────────────────────────┤
  │ Đơn vị quản lý vận hành       │ Có cơ sở khoa học để thiết lập ngưỡng cảnh báo an toàn quan trắc chuyển  │
  │ (EVN, Công ty Thủy điện)      │ vị và độ mở rộng vết nứt thực tế tại công trình.                         │
  ├───────────────────────────────┼──────────────────────────────────────────────────────────────────────────┤
  │ Cơ quan quản lý nhà nước      │ Cung cấp luận cứ hoàn thiện quy chuẩn, tiêu chuẩn quốc gia về kiểm định  │
  │ (Cục Quản lý Đê điều & PCTT)  │ an toàn đập lớn trong thời kỳ biến đổi khí hậu và động đất gia tăng.     │
  └───────────────────────────────┴──────────────────────────────────────────────────────────────────────────┘

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 lý thuyết độc đáo nhất 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). Phương pháp này mở rộng lý thuyết cơ học môi trường liên tục cổ điển bằng cách lượng hóa giá trị ứng suất đại diện thông qua việc tích phân phân bố nội lực trên mặt cắt tính toán thay vì lấy giá trị cục bộ tại điểm nút phần tử, giải quyết triệt để vấn đề "điểm kỳ dị toán học vô hạn" tại chân đập bê tông.

2. Đổi mới phương pháp luận của luận án vượt trội hơn các nghiên cứu trước đây ở điểm nào? So với nghiên cứu của Amir Behshad (2018) và Deeja Alora (2014) sử dụng các phần tử tiếp xúc hoặc hệ lò xo định sẵn (phải giả thiết trước đường nứt), luận án áp dụng kỹ thuật sinh lưới tự thích ứng $h$-refinement dựa trên đánh giá sai số năng lượng Zienkiewicz-Zhu. Vết nứt được tự do phát sinh và lan truyền hoàn toàn tự nhiên theo trạng thái ứng suất thực tế của kết cấu.

3. Phát hiện bất ngờ nhất từ dữ liệu mô phỏng là gì? Phát hiện bất ngờ nhất là các vết nứt sâu tới $6\text{ m}$ tại khối C2, C3 đập Sơn La phát sinh do nhiệt độ co ngót thi công năm 2008 hoàn toàn không phát triển thêm dưới tổ hợp tải trọng tĩnh và động đất vận hành OBE ($a_{max} = 0,114g$), do nằm trong vùng trường ứng suất nén tĩnh áp đảo của đập trọng lực.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) không? Có. Toàn bộ thuật toán lập trình tham số hóa APDL trong môi trường ANSYS, các bước tích phân số học, ma trận phần tử lò xo - cản nhớt nhân tạo và hệ thông số cơ lý chi tiết của đập Koyna, Guandi và Sơn La đều được trình bày minh bạch trong Chương 2 và Chương 4, cho phép tái lập 100% kết quả tính toán.

5. Chương trình nghiên cứu 10 năm tới được vạch ra như thế nào? Tập trung vào mô hình hóa bài toán ghép 3D THM (Thủy - Cơ - Nhiệt), tích hợp áp lực nước dâng động học trong khe nứt và ứng dụng mạng nơ-ron sâu (Deep Neural Networks) để nhận dạng tự động mức độ nguy hiểm của vết nứt từ dữ liệu hệ thống đo đạc thời gian thực.

Kết luận

Luận án "Nghiên cứu sự hình thành và phát triển vết nứt ảnh hưởng tới khả năng chịu tải của đập trọng lực bê tông đầm lăn trong quá trình vận hành" của Tiến sĩ Hà Thanh Dương đã hoàn thành xuất sắc các mục tiêu đề ra 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 bộ chương trình tính toán chuyên dụng RCCD_FAILURE trên nền tảng ngôn ngữ APDL của ANSYS, tích hợp thuật toán phần tử hữu hạn tự thích ứng mô hình $h$ có khả năng tự động bám bắt quỹ đạo lan truyền vết nứt.
  2. Đề xuất và kiểm chứng thành công Phương pháp Ứng suất tương đương (ƯSTĐ-PTHH), khắc phục triệt để nhược điểm tập trung ứng suất kỳ dị của phương pháp PTHH truyền thống và sự thiếu chính xác của phương pháp Sức bền vật liệu.
  3. Làm sáng tỏ cơ chế phá hoại động lực học của đập BTĐL dưới tác động của động đất mạnh, xác định chính xác hai vị trí xung yếu nhất là chân thượng lưu và vùng thay đổi độ dốc thân đập hạ lưu.
  4. Cung cấp bằng chứng định lượng khẳng định đập thủy điện Sơn La bảo đảm an toàn chịu lực tổng thể dưới các điều kiện vận hành thiết kế dù tồn tại các vết nứt thi công ban đầu.
  5. Thiết lập hệ thống biên đàn hồi - cản nhớt mô phỏng chính xác tương tác động lực học giữa công trình và bán không gian nền vô hạn.

Công trình tạo lập một chuẩn mực mới trong mô phỏng rạn nứt kết cấu thủy công tại Việt Nam, đóng góp nền tảng vào kho tàng khoa học công nghệ đập lớn quốc tế, khẳng định năng lực tự chủ nghiên cứu khoa học kỹ thuật cao của các nhà khoa học công trình thủy Việt Nam.