Tổng quan nghiên cứu

Thảm họa động đất luôn là hiểm họa khó lường đối với các công trình xây dựng trên toàn cầu. Các sự kiện lịch sử như trận động đất tại Chi Lê năm 2010 với cường độ 8.8 độ Richter gây tổn thất hơn 30 tỷ USD, hay thảm họa kép tại Nhật Bản năm 2011 mạnh 9.0 độ Richter khiến 15.896 người thiệt mạng và thiệt hại 360 tỷ USD đã gióng lên hồi chuông cảnh báo về an toàn công trình. Tại Việt Nam, dù không nằm trên vành đai địa chấn mạnh nhất, số liệu quan trắc năm 2017 của Viện Vật lý Địa cầu ghi nhận hàng loạt trận động đất có cường độ từ 2.0 đến 4.0 độ Richter tại Quảng Ninh, Điện Biên, Thái Nguyên. Tiêu chuẩn xây dựng Việt Nam TCXDVN 375:2006 đã đặt ra yêu cầu cấp thiết về việc thiết kế kết cấu kháng chấn nhằm bảo vệ tính mạng con người và tài sản.

Về mặt cơ học, tải trọng động đất tác dụng lên công trình thông qua lực quán tính tỷ lệ thuận với khối lượng bản thân. Nếu tăng độ cứng bằng cách phóng to tiết diện cấu kiện, trọng lượng kết cấu tăng theo, kéo theo lực quán tính lớn hơn và làm gia tăng nội lực nguy hiểm. Do đó, giải pháp lắp đặt các thiết bị tiêu tán năng lượng thụ động như hệ cản nhiều khối lượng (Multiple Tuned Mass Dampers - M-TMD) phân bố dọc theo chiều cao nhà trở thành hướng tiếp cận tối ưu.

Mục tiêu cốt lõi của nghiên cứu là đánh giá định lượng ảnh hưởng của thông số tỷ số khối lượng đến hiệu quả giảm chấn của hệ M-TMD phân bố theo dạng dao động riêng của kết cấu khung 20 tầng. Nghiên cứu thực hiện phân tích đáp ứng động lực học của hệ kết cấu dưới tác dụng của tải trọng điều hòa tuần hoàn và dữ liệu gia tốc nền thực tế từ trận động đất El Centro 1940 và Superstition Hills. Kết quả nghiên cứu xác lập cơ sở khoa học để giảm thiểu chuyển vị đỉnh và lực cắt đáy từ 25% đến hơn 40%, hỗ trợ tối ưu hóa thiết kế kết cấu cao tầng an toàn và tiết kiệm chi phí.

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 phát triển trên nền tảng lý thuyết động lực học công trình và lý thuyết điều khiển dao động bị động của Frahm, Den Hartog và Warburton, kết hợp các mở rộng của Roman Lewandowski và Justyna Grzymislawska về hệ cản nhiều khối lượng. Khung lý thuyết tích hợp ba khái niệm nền tảng: hệ cản khối lượng đơn lẻ (S-TMD), hệ cản nhiều khối lượng phân bố (M-TMD), và phân tích phổ tần số bằng biến đổi Fourier.

Phương trình vi phân dao động của hệ kết cấu khung phẳng nhiều bậc tự do có gắn hệ cản M-TMD chịu tác động động lực học được mô tả thông qua hệ ma trận ghép nối tổng thể:

$$M_{total} \ddot{u}(t) + C_{total} \dot{u}(t) + K_{total} u(t) = P(t)$$

Trong đó $M_{total}$, $C_{total}$, $K_{total}$ lần lượt là các ma trận khối lượng, cản và độ cứng tổng thể có kích thước mở rộng tương ứng với số bậc tự do của khung chính và các thiết bị cản gắn thêm. Tỷ số khối lượng $\mu$ được định nghĩa bằng tỷ số giữa tổng khối lượng của hệ cản M-TMD và khối lượng hữu hiệu của dạng dao động chính. Các thông số độ cứng và hệ số cản của từng bộ cản tại mỗi tầng được tính toán dựa trên biên độ của dạng dao động thứ nhất, với điều kiện tổng hệ số tham gia dao động tương đương của các mode được chọn đạt tối thiểu 90% ($\sum \epsilon_i \ge 0.9$). Biến đổi Fourier phức được ứng dụng nhằm chuyển đổi tín hiệu sóng địa chấn từ miền thời gian sang miền tần số, giúp xác định các đỉnh tần số trội để điều chỉnh tần số riêng của hệ cản $\omega_d$ bám sát tần số riêng $\omega$ của công trình.

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

Nghiên cứu áp dụng phương pháp giải tích số mô phỏng trên mô hình phần tử hữu hạn rời rạc hóa. Không gian khảo sát chính là mô hình khung phẳng 20 tầng với 20 bậc tự do động lực học ngang đại diện cho chuyển vị sàn. Cỡ mẫu nghiên cứu bao gồm: mô hình kiểm chứng kết cấu 16 tầng để đối chiếu độ tin cậy với các công bố quốc tế, và mô hình thực nghiệm số khung 20 tầng chịu 3 kịch bản kích động tải trọng (tải điều hòa, băng gia tốc El Centro 1940 với gia tốc đỉnh 0.34g, và băng gia tốc Superstition Hills).

Lý do lựa chọn phương pháp tích phân từng bước Newmark theo sơ đồ gia tốc trung bình với tham số $\gamma = 0.5$ và $\beta = 0.25$ là nhờ tính ổn định vô điều kiện và độ chính xác cao trong việc giải phương trình vi phân cấp 2 phi tuyến theo miền thời gian với bước chia thời gian $\Delta t = 0.01$ giây. Toàn bộ thuật toán ghép nối ma trận và giải lặp Newmark được lập trình hoàn chỉnh trên môi trường MATLAB. Quá trình thu thập dữ liệu và xử lý mô phỏng được thực hiện xuyên suốt trong giai đoạn từ năm 2017 đến năm 2018.

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

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

Phân tích số liệu từ các kịch bản kích động mang lại các phát hiện kỹ thuật quan trọng:

  • Dưới tác động của tải trọng điều hòa cộng hưởng, hệ cản M-TMD với tỷ số khối lượng $\mu = 0.5%$ giúp giảm chuyển vị ngang tại tầng đỉnh từ 35.0% đến 45.2% so với trường hợp kết cấu ban đầu không gắn cản.
  • Khi chịu sóng địa chấn El Centro 1940, biên độ chuyển vị ngang đỉnh của khung 20 tầng giảm rõ rệt từ 28.5% (khi $\mu = 0.5%$) lên mức 38.2% (khi $\mu = 2.0%$). Đồng thời, lực cắt đáy tại tầng 1 giảm từ 18.4% đến 32.1%, giảm thiểu đáng kể ứng suất uốn chân cột.
  • Đối với kích động động đất Superstition Hills, sự phân bố M-TMD dọc theo chiều cao theo dạng dao động bậc 1 đã phân tán năng lượng đều khắp các sàn, giúp độ giảm lực cắt tầng duy trì ổn định trên 25.0% ở toàn bộ 20 tầng.
  • Khi khảo sát dải biến thiên của tỷ số khối lượng $\mu$ từ 0.5% đến 5.0%, kết quả cho thấy hiệu quả giảm chấn tăng phi tuyến. Khi $\mu$ vượt quá ngưỡng 3.0%, mức tăng biên độ giảm chấn có xu hướng bão hòa (chỉ tăng thêm khoảng 2.0% đến 3.5%), trong khi phụ tải tĩnh đè lên sàn tăng lên đáng kể.

Thảo luận kết quả

Hiệu quả triệt tiêu dao động vượt trội của M-TMD phân bố xuất phát từ cơ chế tạo ra các lực cản quán tính lệch pha với chuyển động của kết cấu sàn. Khác với hệ cản đơn S-TMD đặt tại tầng đỉnh dễ bị mất tác dụng khi tần số kích động thay đổi, việc chia nhỏ khối lượng cản và gắn tại tất cả các tầng theo biên độ dạng dao động giúp hệ thống thích ứng tốt với các dải tần số rộng của sóng động đất.

Dữ liệu mô phỏng được thể hiện rõ ràng qua các biểu đồ lịch sử thời gian chuyển vị (time-history displacement), phổ năng lượng phản ứng Fourier và bảng so sánh chuyển vị cực đại. Khi đối chiếu với các công bố quốc tế của Said Elias (2014) về hệ cản phân bố trên nhà 76 tầng (đạt mức giảm chuyển vị 50.0% so với 15.0% của S-TMD đỉnh), kết quả của luận văn thể hiện sự tương đồng cao về quy luật ứng xử. Điều này chứng minh rằng việc bố trí hệ cản bám sát mode shape mang lại hiệu suất tiêu tán cơ năng tối ưu nhất trên từng đơn vị khối lượng cản lắp đặt.

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

Dựa trên kết quả nghiên cứu định lượng, các giải pháp kỹ thuật cụ thể được đề xuất như sau:

  • Chuẩn hóa quy trình tính toán phân bố M-TMD trong thiết kế: Các kỹ sư kết cấu và đơn vị tư vấn thiết kế cần áp dụng quy trình xác định thông số cản $k_{d,i}$ và $c_{d,i}$ dựa trên ma trận dạng dao động bậc 1 của công trình, hướng tới mục tiêu giảm chuyển vị tổng thể tối thiểu 30.0% trong giai đoạn thiết kế kỹ thuật từ 30 đến 45 ngày.
  • Khống chế tỷ số khối lượng tối ưu trong phạm vi 1.0% đến 2.0%: Chủ đầu tư và nhà quản lý dự án nên phê duyệt phương án chọn $\mu$ trong khoảng 1.0% đến 2.0% khối lượng hữu hiệu, giúp tiết kiệm 15.0% đến 20.0% chi phí kết cấu chịu lực phụ trợ mà vẫn đảm bảo 85.0% đến 90.0% hiệu năng giảm chấn cực đại.
  • Tích hợp module giải thuật Newmark vào phần mềm tính toán: Các viện nghiên cứu và bộ phận phát triển kỹ thuật cần tự động hóa thuật toán Newmark trên nền tảng Python hoặc MATLAB liên kết với SAP2000 và ETABS, nhằm rút ngắn 50.0% thời gian phân tích động lực học công trình cao tầng trong vòng 6 tháng triển khai.
  • Bổ sung hướng dẫn lắp đặt thiết bị cản phân bố vào quy chuẩn quốc gia: Bộ Xây dựng và các hội đồng chuyên ngành cần xây dựng khung chỉ dẫn kỹ thuật áp dụng M-TMD phân bố cho các công trình cấp I từ 20 tầng trở lên, với lộ trình ban hành trong vòng 12 đến 24 tháng.

Đố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 chuyên sâu cho bốn nhóm đối tượng chính:

  • Kỹ sư thiết kế kết cấu công trình dân dụng và công nghiệp: Nắm vững phương pháp thiết lập ma trận khối lượng, ma trận độ cứng ghép nối và thuật toán giải số Newmark để ứng dụng trực tiếp vào các dự án cao ốc từ 15 đến 40 tầng chịu tải trọng gió bão và động đất.
  • Nhà nghiên cứu, giảng viên và học viên cao học ngành Kỹ thuật Xây dựng: Sử dụng khung lý thuyết mẫu mực về tối ưu hóa M-TMD và phân tích phổ Fourier làm tài liệu tham khảo để phát triển tiếp các đề tài về hệ cản bán chủ động (Semi-Active TMD) hoặc cản phi tuyến.
  • Chủ đầu tư và các ban quản lý dự án nhà cao tầng: Có cơ sở khoa học để đánh giá hiệu quả kinh tế - kỹ thuật của giải pháp giảm chấn phân bố, giúp giảm tĩnh tải từ 10.0% đến 15.0% so với phương án dùng quả cầu cản đơn nặng hàng trăm tấn đặt tại đỉnh.
  • Cơ quan thẩm tra và quản lý chất lượng công trình xây dựng: Sử dụng các tiêu chí và ngưỡng suy giảm phản ứng động lực học làm căn cứ thẩm định an toàn chịu lực theo tiêu chuẩn TCXDVN 375:2006 cho các dự án trọng điểm đô thị.

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

Hệ cản nhiều khối lượng M-TMD có ưu điểm gì vượt trội so với hệ cản đơn S-TMD? Thay vì tập trung một khối lượng khổng lồ tại đỉnh công trình như quả cầu thép 660 tấn ở tháp Taipei 101, hệ M-TMD chia nhỏ tải trọng ra tất cả các sàn theo biên độ dạng dao động. Cách bố trí này giảm thiểu áp lực tập trung lên hệ dầm cột tầng trên cùng, kiểm soát đa dạng dao động tốt hơn và nâng cao hiệu quả giảm chấn thêm 15.0% đến 35.0%.

Tại sao tỷ số khối lượng $\mu$ nên được lựa chọn trong khoảng 1.0% đến 2.0%? Mô phỏng trên khung 20 tầng cho thấy khi tăng $\mu$ từ 0.5% lên 2.0%, hiệu quả giảm chuyển vị đỉnh tăng mạnh từ 28.5% lên 38.2%. Tuy nhiên, khi tăng $\mu$ vượt mức 3.0%, hiệu quả chỉ tăng thêm từ 2.0% đến 3.0% nhưng làm tăng đáng kể chi phí chế tạo thiết bị và tải trọng tĩnh tác dụng lên khung nhà.

Vai trò của phương pháp Newmark trong việc giải bài toán dao động là gì? Phương pháp tích phân từng bước Newmark với tham số gia tốc trung bình ($\gamma = 0.5, \beta = 0.25$) cho phép giải hệ phương trình vi phân cấp 2 mở rộng với độ chính xác cao. Thuật toán hoạt động ổn định vô điều kiện theo bước thời gian $\Delta t = 0.01$ giây, cung cấp đầy đủ các giá trị chuyển vị, vận tốc và gia tốc tức thời của từng tầng sàn.

Tại sao cần thực hiện phân tích Fourier cho các dữ liệu sóng động đất? Biến đổi Fourier giúp chuyển đổi tín hiệu sóng gia tốc phức tạp sang miền tần số để nhận diện chính xác các đỉnh tần số trội của trận động đất (như sóng El Centro 1940). Dữ liệu này là điều kiện tiên quyết để kỹ sư điều chỉnh dải tần số riêng của hệ cản $\omega_{d,i}$ trùng khớp với tần số kích động chính nhằm triệt tiêu hiện tượng cộng hưởng.

Mô hình M-TMD phân bố phù hợp nhất với loại công trình xây dựng nào? Giải pháp M-TMD phân bố đặc biệt thích hợp cho các công trình khung bê tông cốt thép hoặc khung thép cao từ 15 đến 30 tầng có tính mềm dẻo cao. Thực tế thế giới tại các tháp như John Hancock (240m) hay Chiba Port Tower (125.1m) đã chứng minh các hệ cản khối lượng giúp công trình đứng vững trước gió bão vận tốc 60 m/s và động đất cấp 7.

Kết luận

  • Nghiên cứu đã xây dựng thành công mô hình toán học và thuật toán số hoàn chỉnh để phân tích phản ứng động lực học của khung 20 tầng có gắn hệ cản M-TMD.
  • Xác lập quy luật phân bố khối lượng cản tối ưu theo biên độ dạng dao động bậc 1, giúp phân tán đều năng lượng tiêu tán trên toàn bộ chiều cao công trình.
  • Chứng minh định lượng rằng tỷ số khối lượng $\mu = 1.0% - 2.0%$ mang lại hiệu quả kinh tế và kỹ thuật cao nhất, giúp giảm chuyển vị đỉnh và lực cắt đáy từ 30.0% đến 40.0%.
  • Cung cấp bộ công cụ lập trình MATLAB tin cậy, làm tiền đề cho việc tự động hóa thiết kế kháng chấn trong kỹ thuật xây dựng công trình dân dụng và công nghiệp.
  • Định hướng nghiên cứu tiếp theo trong vòng 6 đến 12 tháng tới cần tập trung vào việc thử nghiệm mô hình thực nghiệm trên bàn rung và mở rộng sang các kết cấu không gian chịu xoắn phức tạp.

Các kỹ sư và đơn vị tư vấn thiết kế hãy áp dụng ngay các công thức tối ưu hóa thông số M-TMD từ luận văn để nâng cao tính an toàn và tối ưu hóa chi phí cho các công trình cao tầng hiện đại.