Tổng quan nghiên cứu

Công trình thủy điện và hồ chứa nước đóng vai trò then chốt trong việc phát triển kinh tế, điều tiết nguồn nước và cung cấp năng lượng. Đập bê tông trọng lực là một trong những loại hình kết cấu ngăn dòng chủ đạo tại các vị trí có điều kiện địa chất phù hợp. Tuy nhiên, tính an toàn của công trình luôn phải đối mặt với rủi ro địa chấn. Thống kê toàn cầu cho thấy khoảng 75% các trận động đất tập trung tại vành đai Thái Bình Dương, 23% tại vành đai Địa Trung Hải và Biển Đông, và 2% xảy ra sâu trong nội địa. Tại Việt Nam, lịch sử địa chấn từng ghi nhận các trận động đất mạnh như trận Điện Biên năm 1935 với cường độ 6,75 độ Richter trên đới đứt gãy Sông Mã, trận Tuần Giáo năm 1983 đạt 6,8 độ Richter trên đới đứt gãy Sơn La, và trận động đất ngoài khơi Nam Trung Bộ năm 1923 đạt 6,1 độ Richter gắn liền với hoạt động phun trào núi lửa Hòn Tro.

Trước nguy cơ động đất đe dọa an toàn hồ đập, phương pháp thiết kế truyền thống dựa trên công thức Westergaard thường đơn giản hóa áp lực thủy động thành khối lượng nước bổ sung gắn vào mặt đập thượng lưu, đồng thời giả định nền liên kết ngàm cứng hoàn toàn. Cách tiếp cận này bỏ qua sự tương tác động lực học đa chiều giữa khối nước hồ chứa, thân đập và nền đá, dẫn đến nguy cơ đánh giá sai lệch biến dạng và trạng thái ứng suất thực tế.

Nghiên cứu được thực hiện nhằm phân tích toàn diện chuyển vị và ứng suất của đập bê tông trọng lực dưới tác động của tải trọng động đất, thông qua việc mô phỏng đồng thời hệ thống tương tác nước hồ chứa - thân đập - nền công trình. Mô hình nghiên cứu áp dụng cho một đập bê tông trọng lực điển hình có chiều cao 37 m, bề rộng đỉnh 5 m, bề rộng đáy 27 m, với chiều sâu cột nước tại mực nước dâng bình thường là 35 m. Phạm vi tính toán không gian mở rộng 100 m về cả phía thượng lưu và hạ lưu. Luận văn được triển khai và hoàn thành trong giai đoạn từ tháng 7 năm 2017 đến tháng 12 năm 2017, bảo vệ thành công vào tháng 1 năm 2018 tại Trường Đại học Bách Khoa - Đại học Quốc gia Thành phố Hồ Chí Minh. Kết quả nghiên cứu giúp nâng cao độ chính xác trong dự báo ứng xử động học thêm 20% đến 25%, giảm thiểu nguy cơ nứt gãy cục bộ và tối ưu hóa từ 10% đến 15% khối lượng bê tông và cốt thép gia cường.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Hệ thống lý thuyết trong nghiên cứu xây dựng dựa trên ba trụ cột chính của cơ học công trình và thủy động lực học:

  1. Lý thuyết động lực học kết cấu nhiều bậc tự do (MDOF): Mô tả phản ứng của hệ thống kết cấu dưới tác động của kích động nền biến thiên theo thời gian. Phương trình vi phân dao động tổng quát có dạng ma trận: $[M]{\ddot{u}(t)} + [C]{\dot{u}(t)} + [K]{u(t)} = {P(t)}$, trong đó $[M]$ là ma trận khối lượng, $[C]$ là ma trận cản, $[K]$ là ma trận độ cứng, ${u(t)}$ là vectơ chuyển vị nút và ${P(t)}$ là vectơ tải trọng kích động. Phương pháp phân tích lịch sử thời gian tuyến tính kết hợp nguyên lý chồng chất dao động theo hệ tọa độ mode và tích phân Duhamel được dùng để xác định đáp ứng động học tại từng thời điểm với tỷ số cản nhớt quy chuẩn $\xi = 0,05$ (5%).

  2. Lý thuyết lan truyền sóng âm trong chất lỏng và tương tác chất lỏng - kết cấu (Fluid-Structure Interaction - FSI): Áp lực thủy động của khối nước hồ chứa được mô tả qua phương trình sóng âm Helmholtz $\nabla^2 P - \frac{1}{c^2}\frac{\partial^2 P}{\partial t^2} = 0$, trong đó $c = 1460$ m/s là vận tốc truyền âm trong môi trường nước, mật độ khối lượng riêng của nước $\rho = 1000$ kg/m³ ($1$ T/m³), và hệ số nhớt động lực học $\mu = 0,005$.

  3. Lý thuyết tiếp xúc phi tuyến bề mặt: Áp dụng các thuật toán Pure Penalty và Augmented Lagrange để mô tả tương tác cơ học tại mặt phân giới giữa đáy đập và nền đá. Lực pháp tuyến tiếp xúc $F_{normal} = K_{normal} \times p$ (với $K_{normal}$ là độ cứng tiếp xúc, $p$ là độ xuyên sâu giả định) và lực ma sát trượt tiếp tuyến $F_{tangential} = \mu \times F_{normal}$ cho phép tái hiện cơ chế tách khe, trượt đàn hồi và truyền ứng suất thực tế, khắc phục nhược điểm giả định ngàm cứng tuyệt đối.

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

Nghiên cứu lựa chọn phương pháp mô phỏng số phần tử hữu hạn (FEM) làm công cụ phân tích chủ đạo:

  • Dữ liệu và phân chia lưới mô hình: Bài toán được thiết lập dưới trạng thái biến dạng phẳng (Plane Strain) với bề dày đơn vị 1 m. Vật liệu thân đập bê tông có mô đun đàn hồi $E = 26{,}000{,}000$ kN/m² ($26$ GPa), hệ số Poisson $\nu = 0,2$, khối lượng thể tích $\rho = 24$ kN/m³. Nền đá được phân lớp gồm tầng đá cứng chịu lực trực tiếp ($E_1 = 1{,}244{,}220$ kN/m², $\nu_1 = 0,2$) và tầng đá phong hóa kẹp giữa ($E_2 = 911{,}970$ kN/m², $\nu_2 = 0,25$). Miền chất lỏng thượng lưu mở rộng 100 m được rời rạc hóa bằng phần tử chất lỏng âm học 2D bốn nút Fluid29 trong ANSYS với 3 bậc tự do mỗi nút (chuyển vị $U_x$, $U_y$ và áp suất âm học $P$).

  • Cỡ mẫu và phương pháp chọn chuỗi gia tốc động đất: Nghiên cứu sử dụng chuỗi dữ liệu gia tốc kế thực tế từ trận động đất Altadena với cỡ mẫu gồm 1000 điểm đo thời gian liên tục, bước thời gian tích phân $\Delta t = 0,02$ s, tổng thời gian kích động 20,0 s, và gia tốc cực đại nền (PGA) đạt $1,756$ m/s² (tương đương khoảng $0,18$g). Mẫu dữ liệu này được lựa chọn có chủ đích vì sở hữu dải tần số kích động đặc trưng cho các chấn động gần nguồn đứt gãy, phù hợp để kiểm chứng giới hạn an toàn của đập bê tông trọng lực.

  • Phương pháp phân tích so sánh đối chuẩn: Tiến hành mô phỏng song song trên hai phần mềm chuyên dụng ANSYS và SAP2000. Trong SAP2000, tương tác nước - đập được quy đổi thành lực phân bố tĩnh và lực thủy động Westergaard với liên kết đáy ngàm cứng. Trong ANSYS, hệ thống được giải quyết đồng thời qua phần tử chất lỏng Fluid29 và phần tử tiếp xúc tiếp tuyến. Lý do lựa chọn mô hình so sánh kép này nhằm bộc lộ rõ sự khác biệt giữa phương pháp tính toán gần đúng truyền thống và phương pháp tương tác động học hoàn chỉnh.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình mô phỏng và phân tích số liệu trên hai mô hình tính toán đã mang lại các phát hiện khoa học quan trọng:

  1. Phản ứng chuyển vị dưới tải trọng thủy tĩnh: Dưới áp lực cột nước thượng lưu cao 35 m, chuyển vị ngang $U_x$ của cả hai mô hình đều có xu hướng tăng tịnh tiến từ chân đập lên đỉnh đập. Tuy nhiên, tại vị trí chân đập thượng lưu, mô hình ANSYS sử dụng phần tử tiếp xúc ghi nhận chuyển vị trượt vi mô đạt $U_x = 0,002213$ m ($2,21$ mm), trong khi mô hình SAP2000 do ngàm cứng nên chuyển vị chân đập bị triệt tiêu về $0$ mm. Tại đỉnh đập, chuyển vị ngang trong SAP2000 lớn hơn mô hình ANSYS khoảng 12% đến 15% do không xét đến sự tiêu tán năng lượng qua biến dạng trượt của nền.

  2. Đáp ứng động học trong miền thời gian động đất: Đồ thị lịch sử thời gian chuyển vị của 1000 điểm đo chỉ ra rằng biên độ dao động lớn nhất tập trung trong khoảng thời gian kích động mạnh từ $2,0$ s đến $5,0$ s, trùng khớp với thời điểm gia tốc nền đạt đỉnh $1,756$ m/s². Đỉnh đập là vị trí có biên độ dao động ngang lớn nhất, đạt giá trị cực đại tại thời điểm $t = 3,42$ s.

  3. Phân bố ứng suất pháp đứng $S_y$ và ứng suất pháp ngang $S_x$: Ứng suất nén đứng $S_y$ đạt giá trị lớn nhất tại chân đập phía hạ lưu (vượt mức $-652$ kN/m² trong điều kiện tĩnh và tăng vọt lên trên $-1200$ kN/m² khi xảy ra động đất). Ở mép thượng lưu, ứng suất $S_x$ chịu tác động mạnh của dao động sóng áp lực nước, ghi nhận giá trị tĩnh $-236,62$ kN/m² tại cao độ $+0$ m và $-197,41$ kN/m² tại cao độ $+13$ m.

  4. Hiện tượng tập trung ứng suất tại vị trí thay đổi tiết diện: Tại điểm chuyển tiếp hình học giữa đoạn đỉnh đập thẳng đứng và mái dốc hạ lưu (cao độ từ $+13$ m đến $+26$ m), mô hình SAP2000 xuất hiện nồng độ ứng suất cục bộ rất cao (đạt $-207,68$ kN/m² tại nút số 2), trong khi ANSYS cho kết quả phân phối mượt mà hơn ở mức $-81,32$ kN/m² tại cao độ $+26$ m nhờ cơ chế truyền lực linh hoạt qua phần tử tiếp xúc.

Thảo luận kết quả

Sự khác biệt rõ nét giữa hai mô hình bắt nguồn từ bản chất cơ học của liên kết và việc mô phỏng môi trường nước. Trong SAP2000, giả thiết nút cứng tuyệt đối tại đáy đập tạo ra các điểm kỳ dị toán học (stress singularity), làm gia tăng nhân tạo ứng suất nén tại góc hạ lưu và triệt tiêu khả năng trượt đáy. Ngược lại, mô hình ANSYS sử dụng phần tử tiếp xúc phi tuyến cho phép bề mặt đáy trượt tương đối $2,21$ mm và tự điều chỉnh trạng thái cân bằng lực, phản ánh đúng hơn bản chất địa kỹ thuật của nền đá nứt nẻ có lớp phong hóa ($E_2 = 911{,}970$ kN/m²).

Mặt khác, việc sử dụng phần tử Fluid29 cho khối nước thượng lưu giúp mô phỏng chính xác hiện tượng sóng âm phản xạ và áp lực thủy động biến thiên pha theo từng bước thời gian $0,02$ s. Điều này hoàn toàn phù hợp với các công trình nghiên cứu quốc tế của Manzori và cộng sự (2012) cũng như nghiên cứu của Ghaemian trên đập Pine Flat, khẳng định rằng nếu bỏ qua tính nén được của nước và độ mềm của nền đá, phản ứng chuyển vị đỉnh đập có thể bị ước tính sai lệch từ 10% đến 18%.

Toàn bộ diễn biến số liệu trong nghiên cứu được trực quan hóa thông qua hệ thống 109 hình vẽ và bảng tra cứu, bao gồm biểu đồ đẳng trị ứng suất 2D, đồ thị lịch sử thời gian chuyển vị 1000 bước tích phân, và các bảng so sánh trực tiếp giá trị ứng suất $S_x, S_y$ tại 4 mặt cắt cao độ trọng yếu ($+0$ m, $+13$ m, $+26$ m và $+37$ m). Cách trình bày này giúp các kỹ sư dễ dàng nhận diện vùng nguy hiểm để bố trí cốt thép chống nứt hiệu quả.

Đề xuất và khuyến nghị

Dựa trên các kết quả mô phỏng động lực học, nghiên cứu đưa ra 4 nhóm giải pháp kỹ thuật cụ thể:

  1. Áp dụng quy trình tính toán tương tác đồng thời FSI trong thiết kế kết cấu: Yêu cầu các đơn vị tư vấn thiết kế công trình thủy lợi - thủy điện bắt buộc áp dụng mô hình liên hợp nước - đập - nền cho các công trình cấp I và cấp II (chiều cao đập trên 30 m hoặc dung tích hồ chứa trên 50 triệu m³). Mục tiêu nhằm nâng cao độ tin cậy dự báo ứng suất thêm 20% đến 25%, hoàn thành việc cập nhật quy chuẩn thiết kế trong giai đoạn quý 1 đến quý 4 năm 2026 do Bộ Xây dựng và Bộ Nông nghiệp & Phát triển Nông thôn chủ trì.

  2. Cải tạo hình học tại các vị trí chuyển tiếp tiết diện: Thiết kế vát góc hoặc bo tròn cung lượn với bán kính cong $R \ge 1,5$ m tại các cao độ thay đổi độ dốc thân đập (từ $+13$ m đến $+26$ m). Giải pháp này giúp giảm thiểu từ 30% đến 40% ứng suất kéo cục bộ phát sinh trong 3 giây đầu tiên khi sóng địa chấn truyền qua.

  3. Lắp đặt hệ thống quan trắc địa chấn và chuyển vị tự động: Trang bị cảm biến đo gia tốc nền 3 trục với tần số lấy mẫu $\ge 100$ Hz và hệ thống định vị GNSS đo chuyển vị đỉnh đập thời gian thực với độ chính xác dưới $1$ mm. Mục tiêu phát hiện sớm các dao động bất thường trong vòng 0,5 giây kích hoạt, thực hiện lắp đặt hoàn thiện tại các đập trọng điểm trong giai đoạn 2026-2027 do các công ty cổ phần thủy điện vận hành.

  4. Tăng cường xử lý nền và bố trí khoan giảm áp: Thi công màng chống thấm kết hợp hệ thống khoan phụt vữa gia cố và hành lang thoát nước đáy đập nhằm duy trì hệ số an toàn chống trượt $K_{trượt} \ge 1,30$ dưới tác động của gia tốc nền cực đại vượt $1,5$ m/s². Chủ đầu tư và ban quản lý dự án thực hiện kiểm tra định kỳ 6 tháng một lần đối với hệ thống thoát nước nền.

Đối tượng nên tham khảo luận văn

Tài liệu luận văn mang lại giá trị học thuật và ứng dụng thực tiễn cao cho 4 nhóm đối tượng chuyên môn:

  1. Kỹ sư tư vấn thiết kế công trình thủy lợi và thủy điện: Nắm vững quy trình mô hình hóa số bằng phần tử Fluid29 và phần tử tiếp xúc trong ANSYS để tính toán chuyển vị và ứng suất cho các dự án đập bê tông trọng lực có chiều cao từ 30 m đến 100 m.

  2. Chuyên gia thẩm tra và cán bộ quản lý an toàn hồ đập: Khai thác hơn 10 biểu đồ đối chuẩn chuyển vị và ứng suất giữa ANSYS và SAP2000 làm thước đo đánh giá độ an toàn ổn định trượt lật của công trình trước các kịch bản động đất kích thích từ cấp 6 đến cấp 8 độ Richter.

  3. Giảng viên, nghiên cứu sinh và học viên cao học ngành Công trình Thủy: Sử dụng làm tài liệu tham khảo chuẩn mực về phương pháp số giải bài toán động lực học công trình nhiều bậc tự do, tích phân từng bước Duhamel và mô hình hóa phi tuyến tiếp xúc nền đá.

  4. Đơn vị quản lý vận hành nhà máy thủy điện: Vận dụng các phát hiện về vùng tập trung ứng suất tại chân đập hạ lưu và điểm gãy tiết diện để xây dựng kịch bản ứng phó sự cố, bố trí điểm đặt thiết bị quan trắc vi địa chấn tại các công trình thực tế như Thủy điện Đồng Nai 5, Sông Tranh 2 hoặc Lai Châu.

Câu hỏi thường gặp

  1. Tại sao phương pháp Westergaard truyền thống chưa phản ánh chính xác tác động của động đất lên đập? Phương pháp Westergaard chỉ tính áp lực thủy động dưới dạng một khối lượng nước tĩnh gắn thêm vào thân đập mà bỏ qua tính nén được của nước hồ chứa (vận tốc truyền âm $1460$ m/s) và khả năng hấp thụ sóng chấn động của nền đá, dẫn đến sai lệch biên độ dao động đỉnh đập từ 10% đến 18%.

  2. Phần tử Fluid29 trong ANSYS có ưu thế gì vượt trội so với gán lực trong SAP2000? Phần tử Fluid29 là phần tử âm học 4 nút 2D giải trực tiếp phương trình sóng Helmholtz với 3 bậc tự do mỗi nút ($U_x, U_y, P$). Mô hình này cho phép tái hiện chân thực quá trình tương tác tương hỗ hai chiều giữa dao động thân đập và sóng áp suất chất lỏng trong suốt 1000 bước thời gian tích phân.

  3. Việc sử dụng phần tử tiếp xúc tại đáy đập ảnh hưởng như thế nào đến kết quả phân tích? Phần tử tiếp xúc giải phóng liên kết ngàm cứng cưỡng bức, cho phép phát sinh chuyển vị trượt ngang vi mô $2,21$ mm và xoay nhẹ tại chân đập. Nhờ đó, mô hình loại bỏ được hiện tượng tập trung ứng suất nén giả tạo tại mép góc và phản ánh chính xác hệ số ổn định trượt thực tế.

  4. Vị trí nào trên thân đập bê tông trọng lực có nguy cơ nứt hỏng cao nhất khi xảy ra địa chấn? Nghiên cứu chỉ ra hai vị trí hiểm yếu nhất: chân đập phía hạ lưu chịu ứng suất nén cực đại vượt $-1200$ kN/m², và điểm thay đổi tiết diện hình học tại cao độ $+13$ m đến $+26$ m chịu nồng độ ứng suất cục bộ $-207,68$ kN/m², đòi hỏi phải bo tròn góc lượn và tăng cường cốt thép chịu cắt.

  5. Ý nghĩa của chuỗi gia tốc 1000 bước thời gian với bước tích phân 0,02 giây là gì? Chuỗi gia tốc nền ghi nhận từ trận động đất Altadena có gia tốc cực đại $1,756$ m/s² giúp phân tích chi tiết phản ứng quá độ của kết cấu theo từng khoảng thời gian $0,02$ s, làm rõ thời điểm xuất hiện ứng suất nguy hiểm nhất tập trung trong cửa sổ kích động từ $2,0$ s đến $5,0$ s.

Kết luận

  • Xây dựng thành công mô hình phần tử hữu hạn phân tích động lực học tương tác nước - đập - nền cho công trình đập bê tông trọng lực cao 37 m bằng phần mềm ANSYS và SAP2000.
  • Khẳng định tính ưu việt của phần tử tiếp xúc phi tuyến khi ghi nhận chuyển vị trượt đáy thực tế đạt $2,21$ mm, khắc phục hoàn toàn nhược điểm tập trung ứng suất giả tạo của mô hình ngàm cứng truyền thống.
  • Xác định chính xác hai vị trí xung yếu chịu ứng suất nén và cắt lớn nhất tại chân đập hạ lưu và điểm chuyển tiếp tiết diện ở cao độ $+13$ m đến $+26$ m với ứng suất cục bộ lên tới $-207,68$ kN/m².
  • Phân tích chi tiết miền thời gian kích động 20,0 s với 1000 bước tích phân, chứng minh phản ứng dao động đỉnh đập đạt cực đại trong khoảng thời gian từ $2,0$ s đến $5,0$ s ứng với gia tốc đỉnh $1,756$ m/s².
  • Đóng góp luận cứ khoa học quan trọng phục vụ công tác rà soát tiêu chuẩn thiết kế kháng chấn công trình thủy lợi - thủy điện, định hướng nâng cấp hệ thống quan trắc địa chấn tự động trong lộ trình 2026-2027.

Các chủ đầu tư, kỹ sư thiết kế và cơ quan quản lý hồ đập nên nhanh chóng tham khảo và áp dụng phương pháp mô phỏng số tương tác chất lỏng - kết cấu để nâng cao chất lượng thẩm tra, đảm bảo an toàn tuyệt đối cho các công trình thủy công trước biến động địa chấn phức tạp.