Tổng quan nghiên cứu

Thủy điện hiện đóng góp khoảng 20% tổng sản lượng điện năng toàn cầu và chiếm tỷ trọng xấp xỉ 60% tổng công suất lắp máy của hệ thống điện Việt Nam, với nguồn trữ năng kỹ thuật ước tính đạt 90 tỷ kWh tương đương 21 triệu kW. Trong tổng thể cụm công trình đầu mối, nhà máy thủy điện là hạng mục quan trọng hàng đầu, chiếm tỷ lệ vốn đầu tư xây dựng rất lớn. Kết cấu nhà máy được phân chia thành hai phần riêng biệt: khối bê tông dưới nước chứa tuabin và gian máy hở phía trên nước bố trí máy phát cùng hệ thống cầu trục vận hành.

Trong thực tế thiết kế xây dựng, các đơn vị tư vấn thường chỉ áp dụng phương pháp giả tĩnh bằng cách nhân tải trọng tĩnh của cầu trục với một hệ số động lực học cố định. Cách tiếp cận này chỉ phân tích nội lực trong mặt phẳng hai chiều vuông góc với trục nhà máy mà bỏ qua các lực quán tính dọc trục phát sinh khi cầu trục di chuyển, phanh hãm hoặc va chạm vào gối chắn. Điều này dẫn đến nguy cơ sai lệch đáng kể trong việc đánh giá trạng thái ứng xử thực tế của công trình, tiềm ẩn nguy cơ mất ổn định khi cẩu lắp tổ máy nặng hàng chục tấn.

Mục tiêu cốt lõi của nghiên cứu là xây dựng mô hình phân tích không gian ba chiều phản ánh chính xác trạng thái ứng suất và biến dạng của kết cấu phần trên nước dưới tác động của các chế độ tải trọng động cầu trục. Nghiên cứu được áp dụng tính toán cụ thể cho công trình Nhà máy Thủy điện Suối Sập 3 tại huyện Bắc Yên, tỉnh Sơn La với công suất lắp máy 14 MW gồm 2 tổ máy tuabin Francis trục đứng, được hoàn thành năm 2015. Kết quả nghiên cứu cung cấp cơ sở khoa học giúp tối ưu hóa từ 10% đến 15% khối lượng thép kết cấu, bảo đảm hệ số an toàn công trình đạt mức 1,25 và kéo dài tuổi thọ vận hành trên 50 năm.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng dựa trên sự kết hợp giữa lý thuyết phần tử hữu hạn và lý thuyết động lực học công trình dành cho hệ nhiều bậc tự do. Trong cơ học kết cấu đàn hồi tuyến tính, nguyên lý cực tiểu thế năng toàn phần xác định trạng thái cân bằng bền của hệ khi năng lượng biến dạng toàn phần đạt cực trị. Quan hệ ứng suất và biến dạng bé tuân theo định luật Hooke tổng quát thông qua ma trận đàn hồi với mô đun đàn hồi cốt thép đạt 200.000.000 kN/m2 và hệ số Poisson là 0,3.

Theo nguyên lý d'Alembert, phương trình vi phân chuyển động của hệ nhiều bậc tự do chịu tải trọng động cưỡng bức được mô tả qua hệ phương trình ma trận gồm ba thành phần lực cản trở: ma trận khối lượng tập trung, ma trận cản nhớt Rayleigh tỷ lệ thuận với vận tốc và ma trận độ cứng tổng thể của hệ thanh. Thuật toán tích phân từng bước Newmark với các tham số gia tốc trung bình không đổi được áp dụng nhằm giải phương trình vi phân chuyển động phi tuyến theo từng bước thời gian, giúp xác định chính xác đáp ứng động học của kết cấu mà không phụ thuộc vào phương pháp giải tích giản lược.

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

Nguồn dữ liệu đầu vào của luận văn được thu thập từ hồ sơ thiết kế kỹ thuật Nhà máy Thủy điện Suối Sập 3 kết hợp với chuỗi quan trắc khí tượng thủy văn hơn 15 năm tại trạm Phiềng Hiềng trên lưu vực rộng 252,7 km2. Cỡ mẫu khảo sát bao gồm toàn bộ hệ khung không gian gian máy kích thước 27,8 x 21,1 m, chiều cao cột chịu lực 12,5 m từ cao trình 181,0 m đến 191,8 m cùng các hệ dầm thép chữ I và thanh giằng chữ C.

Phương pháp chọn mẫu điển hình toàn diện được áp dụng nhằm đại diện cho mô hình kết cấu nhà máy thủy điện kiểu hở quy mô vừa và nhỏ phổ biến tại khu vực miền núi phía Bắc. Thay vì tính toán phẳng đơn giản, tác giả lựa chọn phương pháp phân tích phần tử hữu hạn ba chiều trên phần mềm SAP2000 và phân tích động học cục bộ trên ADINA. Lý do lựa chọn bộ công cụ này xuất phát từ khả năng chia nhỏ miền liên tục thành các phần tử thanh không gian, áp đặt điều kiện biên đàn hồi chính xác và mô phỏng tải trọng biến thiên theo chu kỳ bước thời gian 0,01 giây. Toàn bộ quá trình nghiên cứu, mô hình hóa số và xử lý số liệu được thực hiện trong thời gian 12 tháng từ đầu năm 2014 đến tháng 03 năm 2015.

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 bốn kịch bản tải trọng động mang lại những phát hiện định lượng rõ nét về sự phân bố lại nội lực của hệ khung thép:

Thứ nhất, ở trường hợp cầu trục di chuyển dọc gian máy, lực cắt dọc trục V3 và mômen uốn ngoài mặt phẳng tăng từ 25% đến 30% so với phương pháp tính toán tĩnh truyền thống, làm tăng chuyển vị ngang cục bộ tại các mắt khung liên kết dầm mái.

Thứ hai, khi xe con mang tải thực hiện phanh hãm đột ngột, lực xô ngang làm ứng suất kéo cục bộ trong tiết diện thân cột I600 tại đoạn cao trình từ 181,0 m đến 188,8 m tăng vọt lên mức xấp xỉ 185 N/mm2, chiếm 54,4% giới hạn cường độ chịu kéo tính toán của thép AIII đạt 340.000 kN/m2.

Thứ ba, trường hợp cầu trục va chạm vào gối chắn cuối đường ray tạo ra xung lực va đập nguy hiểm nhất, làm chuyển vị tổng thể tại đỉnh cột ở cao trình 191,8 m tăng 42% so với tải trọng tĩnh, đạt giá trị biên độ lớn nhất là 18,6 mm.

Thứ tư, khi nâng bốc vật cẩu đột ngột khỏi mặt sàn lắp ráp ở cao trình 181,0 m, hệ số động lực học đo được trong mô hình dao động từ 1,35 đến 1,48, vượt xa mức hệ số an toàn 1,20 thường được áp dụng ước tính trong các quy trình thiết kế thông thường.

Thảo luận kết quả

Hiện tượng gia tăng đột biến ứng suất và chuyển vị bắt nguồn từ lực quán tính xuất hiện tức thời khi cầu trục đổi trạng thái động học. Trong không gian ba chiều, ma trận cản và ma trận khối lượng tạo nên các dạng dao động riêng phức tạp, khiến lực quán tính dọc trục truyền thẳng vào các thanh giằng chữ C tại cao độ 185,5 m và 187,0 m thay vì chỉ phân bố đều lên móng bê tông.

So sánh với các nghiên cứu kết cấu tĩnh truyền thống, phương pháp giả tĩnh đã bỏ qua tác động xoắn và lực xô dọc, dẫn đến việc đánh giá thấp nội lực nguy hiểm tại các cột biên từ 18% đến 22%. Toàn bộ dữ liệu tính toán được minh chứng trực quan qua các biểu đồ bao mômen M3, biểu đồ phân bố lực cắt V2, V3 và đồ thị lịch sử thời gian của điểm nút 62 trong mô hình số. Kết quả thảo luận khẳng định rằng việc bổ sung hệ giằng không gian và kiểm tra tải trọng động đa chiều là yêu cầu bắt buộc nhằm loại bỏ rủi ro nứt mỏi kết cấu kim loại trong suốt chu kỳ làm việc 50 năm của nhà máy.

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

Nhằm nâng cao chất lượng thiết kế và an toàn vận hành cho kết cấu phần trên nước của các nhà máy thủy điện, luận văn đưa ra bốn nhóm giải pháp cụ thể:

Thứ nhất, cải tiến quy trình tính toán thiết kế kết cấu bằng việc chuẩn hóa mô hình phần tử hữu hạn ba chiều trên SAP2000 hoặc ADINA, áp dụng hệ số động lực học nâng tải tối thiểu từ 1,35 đến 1,45 đối với các dòng cầu trục gian máy; mục tiêu hoàn thành tích hợp quy trình trong vòng 6 tháng bởi các viện nghiên cứu và công ty tư vấn thiết kế thủy điện.

Thứ hai, gia cường cục bộ hệ khung thép bằng cách nâng cao khả năng chịu uốn xoắn của hệ giằng C400 tại cao trình 185,5 m, đồng thời lắp đặt hệ thống giảm chấn thủy lực tại gối chắn cuối đường ray nhằm triệt tiêu từ 35% đến 40% xung lực va đập; giao ban quản lý dự án và nhà thầu xây lắp thi công hoàn tất trong giai đoạn xây lắp 3 tháng.

Thứ ba, thiết lập hệ thống quan trắc dao động tự động thông qua việc lắp đặt 4 đến 6 cảm biến gia tốc và đầu đo laser đo độ võng thời gian thực tại các điểm nút cột I600 nhạy cảm, đặt ngưỡng cảnh báo dịch chuyển tức thời không vượt quá 15 mm; đơn vị vận hành chịu trách nhiệm triển khai trong 90 ngày đầu chạy thử nghiệm tổ máy.

Thứ tư, ban hành quy chế vận hành cẩu trục nghiêm ngặt, giới hạn tốc độ di chuyển có tải không vượt quá 0,5 m/s và nghiêm cấm phanh hãm đột ngột khi nâng vật nặng quá 10 tấn; do bộ phận quản lý kỹ thuật nhà máy kiểm tra và đánh giá định kỳ 12 tháng một lần.

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

Nội dung và phương pháp luận của đề tài mang lại giá trị thực tiễn cho nhiều nhóm chuyên gia:

Kỹ sư tư vấn thiết kế kết cấu công trình thủy: Tiếp cận phương pháp lập sơ đồ tính không gian 3D, tối ưu hóa kích thước hình học dầm cột thép và kiểm soát chính xác các tổ hợp tải trọng động cầu trục phức tạp.

Chủ đầu tư và Ban quản lý dự án thủy điện: Có cơ sở khoa học để thẩm tra, thẩm định hồ sơ thiết kế kỹ thuật, từ đó tiết kiệm từ 8% đến 12% chi phí vật liệu thép kết cấu mà vẫn bảo đảm tiêu chuẩn an toàn công trình cấp III.

Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Xây dựng Thủy lợi - Thủy điện: Sử dụng làm tài liệu tham khảo chuyên sâu về ứng dụng phương pháp phần tử hữu hạn, thuật toán Newmark và phân tích động lực học kết cấu công trình ngầm và hở.

Đội ngũ kỹ sư quản lý vận hành và bảo dưỡng nhà máy: Nắm rõ các vị trí tập trung ứng suất nguy hiểm trên hệ dầm ray cầu trục để lập kế hoạch kiểm tra không phá hủy, chống mỏi kim loại và bảo trì định kỳ hiệu quả.

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

Phương pháp phân tích động trong luận văn khác gì so với phương pháp giả tĩnh thông thường? Phương pháp giả tĩnh chỉ nhân tải trọng thẳng đứng với hệ số động cố định và xét khung phẳng 2D. Phương pháp phân tích động của luận văn giải phương trình vi phân chuyển động theo thời gian thực trong không gian 3D, phản ánh chính xác cả lực xô ngang, lực hãm dọc trục và lực quán tính va chạm.

Chế độ làm việc nào của cầu trục gây ra ứng suất và chuyển vị nguy hiểm nhất cho nhà máy? Trường hợp cầu trục va chạm vào gối chắn cuối đường ray tạo ra xung động lớn nhất, khiến chuyển vị đỉnh cột thép tăng khoảng 42% so với tải tĩnh, đạt mức cực đại 18,6 mm. Tiếp theo là chế độ xe con phanh hãm khi mang tải làm ứng suất kéo trong cột I600 đạt 185 N/mm2.

Tại sao cần phải xét đến lực dọc trục nhà máy khi tính toán hệ khung thép gian máy? Khi cầu trục di chuyển và phanh gấp, lực quán tính dọc trục tác động trực tiếp lên hệ dầm ray cao trình 189,5 m, sinh ra mômen xoắn và lực cắt V3 ngoài mặt phẳng. Nếu bỏ qua lực dọc, kết cấu hệ giằng chữ C sẽ bị tính toán thiếu an toàn từ 18% đến 22%.

Luận văn đã sử dụng các phần mềm chuyên dụng nào để mô phỏng kết cấu Suối Sập 3? Tác giả đã kết hợp xây dựng mô hình phần tử hữu hạn không gian 3D trên phần mềm SAP2000 để tính toán trường hợp di chuyển dọc, phanh hãm và nâng tải, đồng thời sử dụng phần mềm ADINA để phân tích chi tiết ứng xử động phi tuyến khi xảy ra va chạm gối chắn.

Kết quả nghiên cứu của công trình Thủy điện Suối Sập 3 có thể mở rộng áp dụng cho dự án khác không? Quy trình mô hình hóa và phương pháp xác định tải trọng động trong luận văn hoàn toàn có thể áp dụng rộng rãi cho hàng trăm nhà máy thủy điện vừa và nhỏ có kết cấu gian máy khung thép đúc sẵn tương tự trên phạm vi toàn quốc.

Kết luận

Luận văn thạc sĩ đã hoàn thành xuất sắc các mục tiêu nghiên cứu với năm đóng góp cốt lõi:

  • Hệ thống hóa hoàn chỉnh cơ sở lý thuyết động lực học công trình và phương pháp tích phân Newmark giải bài toán dao động kết cấu nhà máy thủy điện.
  • Xây dựng thành công mô hình phần tử hữu hạn không gian ba chiều mô phỏng chính xác sự làm việc của khung thép dưới tác động động lực học của cầu trục.
  • Định lượng chi tiết mức độ gia tăng ứng suất và chuyển vị ngang tại bốn chế độ vận hành đặc trưng của công trình Thủy điện Suối Sập 3.
  • Chứng minh rõ nét những thiếu sót của phương pháp giả tĩnh truyền thống khi bỏ qua lực quán tính và lực cắt dọc trục ngoài mặt phẳng.
  • Đề xuất các thông số thiết kế hệ giằng C400, dầm I600 cùng giải pháp giảm chấn thủy lực đạt hiệu quả giảm thiểu xung lực va đập tới 40%.

Trong giai đoạn tiếp theo từ 6 đến 12 tháng tới, các đơn vị thiết kế cần nhanh chóng cập nhật mô hình tính toán động 3D vào tiêu chuẩn cơ sở. Các chủ đầu tư và kỹ sư công trình hãy áp dụng ngay phương pháp luận từ nghiên cứu này để nâng cao chất lượng thẩm định và bảo đảm an toàn bền vững cho các dự án thủy điện.