Tổng quan nghiên cứu

Thảm họa địa chấn và rung chấn công trình luôn là một trong những thách thức nghiêm trọng nhất đối với sự an toàn của cơ sở hạ tầng nhân loại. Theo thống kê lịch sử từ năm 1755 đến năm 2005, chỉ riêng 10 trận động đất lớn nhất trên thế giới đã cướp đi sinh mạng của hơn 5.227.000 người và gây tổn thất tài sản vượt mức 7,1 tỷ USD. Điển hình là trận siêu địa chấn 9,0 độ Richter tại Nhật Bản năm 2011 kèm sóng thần cao 10 m đã làm hơn 100.000 người thiệt mạng, để lại tổng thiệt hại lên tới 100 tỷ USD, bao gồm 40 tỷ USD phá hủy trực tiếp hạ tầng kỹ thuật. Tại Việt Nam, nguy cơ rủi ro địa chấn ngày càng hiện hữu rõ nét qua các sự kiện rung lắc tại Tuần Giáo năm 1983, ngoài khơi Vũng Tàu năm 2005, động đất kích thích 4,2 độ Richter với độ sâu chấn tiêu 7,3 km tại khu vực thủy điện Sông Tranh 2 năm 2012, cũng như hoạt động của đới đứt gãy sông Sài Gòn có khả năng phát sinh địa chấn mạnh đến 5,5 độ Richter gây rung động cấp 7 tại đô thị trung tâm.

Trước thực trạng đó, việc xác định chính xác các tham số động lực học để đánh giá khả năng chịu lực và chẩn đoán hư hại kết cấu trở thành bài toán sống còn trong quản lý xây dựng. Luận văn tập trung giải quyết vấn đề nhận dạng tần số dao động riêng, dạng dao động và độ cứng từng tầng cho công trình thông qua phương pháp phân giải trong miền tần số (FDD - Frequency Domain Decomposition). Mục tiêu cụ thể là xây dựng mô hình toán học giải tích, thiết lập hệ thống thực nghiệm thu thập dữ liệu gia tốc và thuật toán phân tích tín hiệu trên mô hình kết cấu 2 tầng. Nghiên cứu mang ý nghĩa thực tiễn to lớn trong việc cung cấp công cụ định lượng độ cứng, giúp phát hiện sớm các tầng yếu có nguy cơ sụp đổ, đồng thời tối ưu hóa thông số lắp đặt các thiết bị giảm chấn thụ động TMD hoặc bán chủ động MR Damper nhằm triệt tiêu tới 40% biên độ rung lắc nguy hiểm do động đất gây ra.

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 ứng dụng nền tảng của phương pháp phân tích modal hoạt động (OMA - Operational Modal Analysis), cho phép nhận dạng các đặc tính động học của công trình chỉ dựa trên tín hiệu đáp ứng ngõ ra từ các rung động môi trường xung quanh (ambient vibration) như gió, sóng biển hay hoạt động giao thông mà không cần nguồn kích thích cưỡng bức nhân tạo tốn kém. Phương pháp phân giải trong miền tần số (FDD) do Brincker phát triển được sử dụng làm lý thuyết cốt lõi, mở rộng từ kỹ thuật trích lọc đỉnh cổ điển (Peak Picking). FDD phân giải ma trận mật độ phổ công suất PSD và mật độ phổ chéo CSD của hệ nhiều bậc tự do thành tập hợp các mật độ phổ của hệ một bậc tự do độc lập thông qua giải thuật phân giải giá trị kỳ dị (SVD - Singular Value Decomposition).

Trong mô hình này, tín hiệu kích thích từ môi trường được chuẩn hóa dưới dạng nhiễu trắng Gaussian dừng có kỳ vọng bằng 0, kết hợp giả thiết kết cấu có tỷ số giảm chấn yếu, thông thường dao động quanh mức 5% đối với công trình xây dựng dân dụng chịu động đất. Để xác thực tính chính xác của các mode dao động tách chiết được và loại bỏ các đỉnh phổ nhiễu, chỉ tiêu đảm bảo modal (MAC - Modal Assurance Criterion) do Allemang và Brown đề xuất năm 1982 được áp dụng với ngưỡng tương thích tin cậy trên 80% (tương đương 0,80). Kết cấu tòa nhà được mô hình hóa theo sơ đồ lực cắt dầm (shear beam), trong đó chuyển vị ngang của mỗi tầng phụ thuộc trực tiếp vào các tầng liền kề, làm tiền đề để thiết lập công thức giải tích xác định chính xác ma trận độ cứng K từ các vectơ riêng và tần số dao động riêng nhận dạng được.

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

Nghiên cứu kết hợp chặt chẽ giữa mô phỏng số giải tích và thực nghiệm vật lý đa kênh trên mô hình thực tế. Cỡ mẫu dữ liệu thực nghiệm được xác lập với độ lớn 30.720 điểm đo liên tục trong khoảng thời gian 15 giây cho mỗi kênh cảm biến, đảm bảo đáp ứng đầy đủ định lý lấy mẫu Shannon-Nyquist và ghi nhận trọn vẹn chu kỳ dao động tự do của kết cấu.

Phương pháp chọn mẫu dữ liệu được thực hiện thông qua việc bố trí 3 cảm biến gia tốc áp điện công nghiệp IMI 603C01 tại 3 vị trí chiến lược dọc theo chiều cao mô hình tòa nhà thép carbon 2 tầng: kênh 0 gắn tại nền móng có độ nhạy 94 mV/g, kênh 1 đặt tại sàn tầng 1 có độ nhạy 96 mV/g, và kênh 2 đặt tại sàn tầng 2 có độ nhạy 103 mV/g. Tín hiệu gia tốc được số hóa trực tiếp bằng bộ chuyển đổi tương tự sang số Delta-Sigma 24-bit trên phần cứng NI-9234 thuộc khung máy NI cDAQ-9171 với tốc độ lấy mẫu đồng bộ 2048 mẫu/giây (2048 Hz), phân chia từ tần số thời gian nền chủ 13,1072 MHz.

Lý do lựa chọn phương pháp phân tích FDD kết hợp môi trường LabVIEW và Matlab/Simulink là vì FDD có khả năng tách các dạng mode có tần số gần nhau với độ chính xác cao, tính toán trực tiếp trên miền tần số thông qua phép ước lượng phổ Welch với cửa sổ Hamming di động, giúp triệt tiêu tối đa ảnh hưởng của sai số phân cực dòng common-mode 50 ppm và loại bỏ hoàn toàn yêu cầu đo lực kích thích ngõ vào.

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

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

Quá trình phân tích thực nghiệm và mô phỏng trên mô hình tòa nhà 2 tầng bằng thép carbon có tổng chiều cao 0,753 m, kích thước sàn 0,6 m x 0,4 m và tổng khối lượng thực nghiệm lớn nhất 30,0494 kg (bao gồm 4 đĩa thép 2,233 kg và 2 khối trụ thép 1,0257 kg) đã đem lại những phát hiện quan trọng:

Thứ nhất, giải thuật phân giải giá trị kỳ dị SVD từ ma trận phổ mật độ công suất kích thước 2x2 đã trích xuất rõ nét các đỉnh giá trị kỳ dị đơn, nhận dạng chính xác tần số dao động tự nhiên cơ bản của mode 1 ở mức 4,625 Hz và dạng mode chuyển vị dọc tương ứng. Các đỉnh phổ giả do nhiễu môi trường đã được loại bỏ hoàn toàn khi giá trị MAC chỉ đạt dưới 0,30, trong khi các mode dao động thực tế đều đạt hệ số MAC vượt mức 88%.

Thứ hai, công thức giải tích nghịch đảo ma trận độ cứng theo mô hình lực cắt dầm đã xác định chính xác độ cứng từng tầng với giá trị tầng 1 đạt k1 = 3414,90 N/m và tầng 2 đạt k2 = 3949,38 N/m. Mức độ chênh lệch giữa kết quả nhận dạng FDD từ dữ liệu thực nghiệm so với giá trị tính toán lý thuyết từ mô hình cơ học kết cấu chỉ dao động trong khoảng sai số dưới 4,5%.

Thứ ba, khảo sát phản ứng động học dưới tải trọng xung gió 10 N tác dụng vào sàn tầng 2 trong khoảng thời gian ngắn 5 ms cho thấy biên độ chuyển vị và vận tốc dao động của cả hai tầng nhanh chóng tắt dần và trở về trạng thái ổn định sau khoảng 14 đến 18 giây nhờ tác động của hệ số cản nhớt nội tại.

Thứ tư, khi mô phỏng phản ứng dưới tác động kích thích của gia tốc nền đất dạng hàm sin biên độ ngẫu nhiên đạt cực đại 0,48 m/s2 trong 0,5 giây, hệ thống thể hiện sự khuếch đại dao động mạnh mẽ khi tần số rung động nền đạt ngưỡng 5 Hz, tiệm cận vùng tần số cộng hưởng của công trình.

Thảo luận kết quả

Hiện tượng khuếch đại dao động mất ổn định tại tần số 5 Hz phản ánh đúng bản chất vật lý của sự trùng phùng tần số giữa ngoại lực địa chấn và tần số riêng của hệ kết cấu. Nếu tỷ số cản nhớt không được bổ sung, biên độ dao động chuyển vị tại tầng 2 sẽ tăng vọt gấp hơn 4 lần so với điều kiện làm việc thông thường, đe dọa trực tiếp đến tính bền vững của các liên kết hàn và dầm thép chịu lực.

Dữ liệu nghiên cứu được biểu diễn trực quan và khoa học thông qua hệ thống các đồ thị mật độ phổ công suất tự thân PSD, mật độ phổ chéo CSD và đường cong giá trị kỳ dị đơn trên dải tần số từ 0,01 Hz đến 1000 Hz. Cùng với đó, bảng so sánh đối chuẩn ma trận độ cứng giữa các mode dao động cho phép định vị trực tiếp tầng có độ cứng suy giảm, chứng minh rằng sàn tầng 1 chịu lực cắt lớn hơn nhưng có độ cứng thấp hơn tầng 2 khoảng 13,5%, tạo nên vùng biến dạng uốn tập trung cần được ưu tiên gia cố.

So sánh với các nghiên cứu nhận dạng mạng nơ-ron lan truyền ngược của Olivera Jovanović (1997) hay mô hình thích nghi trực tuyến GADALINE và RARX của Maosheng Gong (2010), phương pháp FDD thể hiện ưu thế vượt bậc về tính tường minh vật lý. Trong khi mạng nơ-ron chỉ dự báo trạng thái dịch chuyển kế tiếp mà không lượng hóa được các tham số kết cấu, thì FDD trích xuất trực tiếp giá trị độ cứng thực tế của từng tầng mà không cần cập nhật mô hình phi tuyến phức tạp.

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

Dựa trên kết quả nghiên cứu lý thuyết và thực nghiệm nhận dạng tham số động học bằng FDD, 4 giải pháp trọng tâm được khuyến nghị thực hiện:

Thứ nhất, tích hợp thiết bị giảm chấn thụ động TMD (Tuned Mass Damper) hoặc bộ giảm chấn từ biến bán chủ động MR Damper cho các công trình cao tầng. Các kỹ sư thiết kế cần sử dụng chính xác tần số dao động riêng nhận dạng được (quanh ngưỡng 4,625 Hz) để hiệu chỉnh tần số làm việc của bộ TMD, hướng tới mục tiêu cắt giảm từ 35% đến 50% biên độ rung lắc đỉnh khi xảy ra dư chấn động đất cấp 6 hoặc cấp 7, hoàn thành việc ứng dụng thử nghiệm trong vòng 18 tháng. Chủ thể thực hiện là các đơn vị tư vấn thiết kế kết cấu và ban quản lý dự án xây dựng.

Thứ hai, triển khai quy trình gia cường kháng chấn cục bộ cho các tầng yếu có độ cứng k đo đạc thực tế thấp hơn ngưỡng thiết kế an toàn 3500 N/m. Sử dụng giải pháp bọc sợi carbon composite (CFRP) hoặc bổ sung hệ giằng thép chữ X để tăng độ cứng tầng lên tối thiểu 25%, tiến hành định kỳ theo chu kỳ kiểm định kỹ thuật 3 năm một lần do các nhà thầu thi công chuyên nghiệp đảm nhiệm.

Thứ ba, xây dựng hệ thống quan trắc sức khỏe kết cấu trực tuyến (SHM) ứng dụng kỹ thuật FDD kết hợp phần cứng thu thập dữ liệu 24-bit với tốc độ quét tối thiểu 2048 Hz. Thiết lập cảnh báo tự động khi phát hiện tần số riêng sụt giảm quá 8% hoặc chỉ số MAC giữa các chu kỳ đo lệch trên 15%, chuyển giao cho các trung tâm kiểm định an toàn công trình thuộc Sở Xây dựng vận hành trong lộ trình 12 tháng tới.

Thứ tư, hoàn thiện quy chuẩn thiết kế kháng chấn bằng cách bắt buộc tính toán kiểm tra kép giữa phương pháp phổ phản ứng nhiều dạng dao động và phương pháp tĩnh lực ngang tương đương theo tiêu chuẩn TCXDVN 375:2006, đặc biệt đối với khu vực nền móng và các tầng từ 50% chiều cao công trình trở xuống, áp dụng ngay vào quy trình thẩm định hồ sơ kỹ thuật của các cơ quan quản lý nhà nước từ quý 4 năm 2024.

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

Nội dung và kết quả của luận văn mang lại giá trị khoa học và ứng dụng thực tiễn cao cho 4 nhóm đối tượng cụ thể:

Nhóm thứ nhất là các kỹ sư kết cấu, chuyên gia tư vấn thiết kế công trình dân dụng và công nghiệp. Luận văn cung cấp phương pháp luận chuẩn xác để trích xuất ma trận độ cứng thực tế, hỗ trợ tối ưu hóa thiết kế kháng chấn cho các tòa nhà nhiều tầng đặt tại khu vực có nguy cơ động đất cấp 7 và cấp 8.

Nhóm thứ hai là các cán bộ kỹ thuật, cơ quan quản lý đô thị và đơn vị kiểm định chất lượng xây dựng. Dữ liệu thực nghiệm và quy trình đo đạc FDD giúp thiết lập phương án đánh giá nhanh chóng mức độ suy thoái, phát hiện tầng yếu của hàng trăm công trình cũ có tuổi thọ trên 20 năm mà không cần phá hủy hay làm gián đoạn công năng sử dụng.

Nhóm thứ ba là các giảng viên, nghiên cứu sinh và học viên cao học thuộc chuyên ngành Kỹ thuật Xây dựng, Cơ kỹ thuật và Kỹ thuật Điều khiển - Tự động hóa. Đề tài đóng vai trò như một tài liệu tham khảo chuyên sâu kết hợp hài hòa giữa cơ học dao động công trình, lý thuyết xử lý tín hiệu số trong LabVIEW và giải thuật ma trận SVD trên Matlab.

Nhóm thứ tư là các doanh nghiệp nghiên cứu, chế tạo thiết bị giảm chấn và giải pháp quan trắc công trình thông minh. Kết quả nhận dạng tham số động học là cơ sở then chốt để lập trình thuật toán điều khiển thích nghi cho các thiết bị giảm chấn TMD và MR Damper, nâng cao hiệu quả triệt tiêu rung động môi trường lên hơn 40%.

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

Phương pháp phân giải trong miền tần số FDD có ưu thế gì vượt trội so với phân tích modal thực nghiệm truyền thống EMA? Phương pháp FDD thuộc nhóm phân tích modal hoạt động OMA, chỉ yêu cầu đo đạc tín hiệu đáp ứng gia tốc ngõ ra từ các rung động môi trường tự nhiên mà không cần trang bị các thiết bị kích thích nhân tạo đắt tiền như búa kích rung công suất lớn hay bàn rung thủy lực. Điều này giúp tiết kiệm hơn 60% chi phí thực nghiệm, bảo vệ an toàn tuyệt đối cho kết cấu công trình hiện hữu và đạt độ phân giải tín hiệu cao nhờ phần cứng thu thập 24-bit.

Chỉ tiêu đảm bảo modal MAC đóng vai trò gì trong việc loại bỏ các tần số dao động giả? Chỉ tiêu MAC lượng hóa mức độ tương đồng giữa hai vectơ dạng dao động thông qua giá trị vô hướng từ 0 đến 1. Trong phân tích FDD, các đỉnh phổ xuất hiện do nhiễu hoặc do nguồn kích thích cục bộ thường có hệ số MAC thấp dưới 0,30. Chỉ các đỉnh phổ đạt hệ số MAC từ 80% (tương đương 0,80) trở lên so với dạng mode chuẩn mới được xác nhận là mode dao động riêng thực tế của tòa nhà.

Vì sao nghiên cứu lại thiết lập tần số lấy mẫu 2048 Hz và thời gian ghi nhận 15 giây? Tốc độ lấy mẫu 2048 mẫu/giây được cấu hình đồng bộ từ tần số xung nhịp chủ 13,1072 MHz của module NI-9234, giúp bao phủ trọn vẹn dải tần số động học rộng và tránh hiện tượng méo tín hiệu do chồng phổ. Khoảng thời gian thu thập 15 giây cung cấp đúng 30.720 điểm dữ liệu, đủ lớn để giải thuật biến đổi Fourier nhanh FFT và phép ước lượng phổ Welch với cửa sổ Hamming đạt độ hội tụ tối ưu mà không làm quá tải bộ nhớ xử lý.

Làm thế nào để xác định độ cứng từng tầng từ các tham số modal nhận dạng được? Dựa trên mô hình lực cắt dầm, phương trình cân bằng động học của hệ dao động được biến đổi giải tích thành mối quan hệ giữa khối lượng tập trung, bình phương tần số góc riêng và chuyển vị tương đối giữa các tầng liền kề. Bằng cách giải hệ phương trình từ tầng mái xuống tầng 1, nghiên cứu xác định chính xác độ cứng tầng 1 là 3414,90 N/m và tầng 2 là 3949,38 N/m với mức sai số đối chuẩn dưới 5%.

Kết quả nhận dạng FDD hỗ trợ việc lắp đặt thiết bị giảm chấn TMD như thế nào? Để thiết bị hấp thụ dao động khối lượng TMD phát huy tối đa hiệu quả triệt tiêu rung chấn, tần số dao động riêng của TMD phải được điều chỉnh trùng khớp với tần số cơ bản của công trình. Nhờ xác định chính xác tần số riêng 4,625 Hz và dạng mode chuyển vị bằng FDD, các kỹ sư có thể tối ưu hóa khối lượng và độ cứng lò xo của TMD, giúp tăng tỷ số cản hệ thống từ 5% lên trên 15%, ngăn chặn hiệu quả hiện tượng cộng hưởng nguy hiểm tại tần số 5 Hz.

Kết luận

  • Luận văn đã hoàn thiện toàn diện quy trình nhận dạng tham số động học kết cấu bằng phương pháp phân giải trong miền tần số FDD dựa trên các rung động ngẫu nhiên xung quanh.
  • Hệ thống thực nghiệm được xây dựng đồng bộ với 3 kênh cảm biến gia tốc công nghiệp và module NI-9234 24-bit, xử lý chính xác tập dữ liệu 30.720 điểm đo ở tần số lấy mẫu 2048 Hz.
  • Nhận dạng thành công tần số dao động riêng cơ bản 4,625 Hz và xác định chính xác ma trận độ cứng từng tầng k1 = 3414,90 N/m, k2 = 3949,38 N/m với sai số dưới 4,5%.
  • Chứng minh rõ nét nguy cơ mất ổn định do cộng hưởng khi công trình chịu tác động rung chấn nền đất tại tần số 5 Hz và khả năng dập tắt dao động trong 18 giây của hệ số cản nhớt.
  • Đề xuất giải pháp định lượng kết cấu phục vụ thiết kế lắp đặt thiết bị giảm chấn TMD, định hình lộ trình mở rộng cho hệ kết cấu nhiều tầng n bậc tự do trong giai đoạn 24 đến 36 tháng tới.

Công trình khẳng định đóng góp khoa học quan trọng trong việc hiện đại hóa công tác kiểm định an toàn, cung cấp công cụ nhận dạng tin cậy, chi phí thấp cho ngành xây dựng và điều khiển công trình. Các đơn vị quản lý, kỹ sư kết cấu và nhà nghiên cứu được khuyến khích khai thác ngay các mô hình giải tích và quy trình đo đạc này để nâng cao năng lực phòng chống rủi ro động đất cho hệ thống công trình hạ tầng đô thị.