Tổng quan nghiên cứu

Trong công tác quản lý và vận hành công trình xây dựng hiện đại, giám sát sức khỏe kết cấu (Structural Health Monitoring - SHM) đóng vai trò sống còn nhằm đảm bảo an toàn sinh mạng và tối ưu hóa chi phí bảo trì. Thực tế cho thấy, chi phí khắc phục sự cố kết cấu khi phát hiện muộn có thể tăng gấp 3 đến 5 lần so với việc bảo dưỡng định kỳ dựa trên chẩn đoán sớm. Đặc biệt đối với hệ khung không gian ba chiều trong xây dựng dân dụng và công nghiệp, các hư hỏng cục bộ như nứt do mỏi, ăn mòn vật liệu hay lỏng liên kết bu lông thường diễn biến âm thầm bên trong tiết diện.

Đề tài luận văn thạc sĩ chuyên ngành Kỹ thuật Xây dựng công trình Dân dụng và Công nghiệp tại Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh (thực hiện trong năm 2017) tập trung giải quyết bài toán xác định đồng thời vị trí và mức độ hư hỏng trên hệ kết cấu khung không gian phức tạp. Nghiên cứu đề xuất quy trình chẩn đoán hai giai đoạn tối ưu: Giai đoạn thứ nhất nhận diện vị trí thanh và phần tử hư hỏng thông qua các chỉ số năng lượng biến dạng; Giai đoạn thứ hai định lượng chính xác tỷ lệ suy giảm độ cứng bằng thuật toán di truyền.

Nghiên cứu mang lại ý nghĩa học thuật và thực tiễn sâu sắc khi nâng cao độ chính xác nhận diện hư hỏng đa điểm lên trên 95% trong điều kiện tiêu chuẩn và duy trì độ tin cậy ổn định khi chịu tác động của mức nhiễu đo lường 3%. Đây là cơ sở khoa học tin cậy giúp các kỹ sư công trình chuyển đổi từ hình thức kiểm tra trực quan thủ công sang hệ thống đánh giá động học tự động hóa.

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 trên nền tảng cơ học kết cấu động lực học và lý thuyết tối ưu hóa tiến hóa hiện đại, tích hợp hai khung lý thuyết cốt lõi:

  1. Lý thuyết dao động tự do không cản của hệ nhiều bậc tự do: Dựa trên phương trình vi phân chuyển động thiết lập qua phương pháp phần tử hữu hạn, giải bài toán trị riêng để xác định các tần số góc riêng và véc tơ dạng dao động đặc trưng. Phần tử thanh không gian được mô hình hóa với 12 bậc tự do (gồm 6 bậc tự do tại mỗi nút: 3 chuyển vị tịnh tiến và 3 chuyển vị xoay), sử dụng ma trận độ cứng cục bộ và ma trận khối lượng tương thích có kích thước 12x12.
  2. Lý thuyết năng lượng biến dạng dạng dao động (Modal Strain Energy - MSE): Năng lượng biến dạng được xác định bởi tích vô hướng của véc tơ chuyển vị dạng dao động và ma trận độ cứng phần tử. Do ma trận độ cứng rất nhạy với sự suy giảm tiết diện, năng lượng biến dạng phản ánh rõ nét sự thay đổi cục bộ của thanh.

Các khái niệm chủ chốt trong luận văn bao gồm: Ma trận độ cứng tổng thể, Tỷ số thay đổi năng lượng biến dạng (MSECR), Chỉ số năng lượng biến dạng cơ bản (MSEBI), Chỉ số tương đương năng lượng biến dạng (MSEEI), và Véc tơ suy giảm độ cứng đại diện cho mức độ hư hại của kết cấu.

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

Luận văn áp dụng phương pháp mô phỏng số thông qua việc lập trình tự động hóa hoàn toàn trên nền tảng phần mềm MATLAB kết hợp kiểm chứng chéo với phần mềm SAP2000:

  • Cỡ mẫu và đối tượng khảo sát: Nghiên cứu khảo sát mô hình khung không gian ba chiều gồm hệ thống 63 thanh phần tử chịu lực chính, được chia nhỏ thành nhiều phần tử hữu hạn rời rạc để kiểm tra hư hỏng cục bộ. Đồng thời, nghiên cứu phân tích tập dữ liệu gồm từ 6 đến 10 dạng dao động đầu tiên (từ mode 1 đến mode 10).
  • Phương pháp chọn mẫu và kịch bản hư hỏng: Phương pháp chọn mẫu kịch bản có chủ đích được áp dụng để bao quát các tình huống: hư hỏng đơn điểm, hư hỏng đa điểm đồng thời trên nhiều thanh (như thanh 1, thanh 15, thanh 48, thanh 63), và hư hỏng cục bộ tại các vị trí khác nhau trên cùng một phân đoạn thanh với mức độ giảm độ cứng giả định từ 10% đến 50%.
  • Lý do lựa chọn phương pháp: Việc lựa chọn giải thuật di truyền (GA) giúp giải quyết triệt để bài toán tối ưu phi tuyến toàn cục, tránh rơi vào cực trị địa phương khi không gian tìm kiếm có số lượng biến số lớn. Quy trình hai bước giúp giảm số biến tối ưu từ toàn bộ phần tử trong kết cấu xuống chỉ còn vài phần tử nghi vấn, tăng tốc độ hội tụ và độ chính xác.
  • Thời gian nghiên cứu: Toàn bộ quá trình xây dựng thuật toán, chạy mô phỏng và kiểm chứng dữ liệu được hoàn thành trong thời gian 6 tháng (từ tháng 1 năm 2017 đến tháng 6 năm 2017).

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

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

  1. Chỉ số MSEBI là công cụ nhận diện vị trí hư hỏng vượt trội nhất: Khi phân tích trên 8 dạng dao động đầu tiên, chỉ số MSEBI xác định chính xác 100% các thanh và vị trí phần tử bị tổn thương mà không xuất hiện cảnh báo giả. Trong khi đó, các chỉ số MSECR và MSEEI ghi nhận tỷ lệ nhận diện nhầm lẫn ở các phần tử lân cận dao động từ 15% đến 25%.
  2. Hàm mục tiêu dựa trên năng lượng biến dạng (OF2) tối ưu hơn hàm mục tiêu dạng dao động (OF1): Khi áp dụng thuật toán di truyền với hàm OF2, sai số xác định mức độ suy giảm độ cứng chỉ dao động trong khoảng từ 0,2% đến 1,8%. Ngược lại, hàm mục tiêu OF1 cho độ lệch lớn hơn, với sai số trung bình lên đến 8,5% đến 14,2%.
  3. Khả năng phát hiện hư hỏng cục bộ trên từng phân đoạn thanh: Phương pháp mở rộng phân tích cục bộ (subMSEBI) đã định vị chính xác vị trí khuyết tật trên thanh với sai số kích thước dưới 2%, ngay cả khi mức độ tổn thương thực tế chỉ làm suy giảm 10% độ cứng ban đầu của phần tử.
  4. Khả năng kháng nhiễu ấn tượng: Dưới tác động của mức nhiễu đo lường 3% đưa vào véc tơ dạng dao động, chỉ số MSEBI vẫn giữ vững khả năng khoanh vùng chính xác 4 vị trí hư hỏng giả định với độ tin cậy vượt trên 90%.

Thảo luận kết quả

Nguyên nhân chính giúp chỉ số MSEBI và hàm mục tiêu OF2 đạt hiệu năng cao là do năng lượng biến dạng tích hợp đồng thời cả thông tin hình học, ma trận độ cứng và độ cong dao động của kết cấu. Điều này tạo ra độ nhạy cao hơn nhiều so với việc chỉ so sánh chuyển vị đơn thuần tại các nút.

Khi so sánh với các công bố quốc tế tiêu biểu trong cùng lĩnh vực, phương pháp hai bước đề xuất trong luận văn đạt kết quả tương đương về độ chính xác (độ lệch dưới 1,5%) nhưng giảm hơn 60% thời gian tính toán do không gian tìm kiếm được thu hẹp đáng kể từ giai đoạn đầu.

Về mặt trực quan hóa dữ liệu, kết quả chẩn đoán có thể được biểu diễn trực quan qua biểu đồ phân bố chỉ số năng lượng biến dạng dọc theo trục các phần tử (từ phần tử 1 đến 63) kết hợp bảng tổng hợp so sánh mức độ hư hỏng thực tế với mức độ chẩn đoán. Biểu đồ cột thể hiện giá trị MSEBI xuất hiện các đỉnh đột biến vượt trội tại các vị trí khuyết tật, giúp người vận hành dễ dàng nhận diện bằng mắt thường trước khi tiến hành tối ưu hóa bằng thuật toán di truyền.

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

Nhằm chuyển giao và ứng dụng hiệu quả phương pháp chẩn đoán vào thực tiễn kỹ thuật công trình, luận văn đề xuất 4 nhóm giải pháp chiến lược:

  1. Tích hợp thuật toán vào hệ thống giám sát tự động: Các công ty công nghệ và đơn vị tư vấn giám sát cần chủ động lập trình nhúng mô-đun phân tích MSEBI và thuật toán di truyền vào phần mềm giám sát thời gian thực, đặt mục tiêu rút ngắn thời gian chẩn đoán sự cố xuống dưới 15 phút sau khi thu thập dữ liệu đo, hoàn thành triển khai trong vòng 6 đến 9 tháng.
  2. Tối ưu hóa sơ đồ bố trí cảm biến đo dao động: Chủ đầu tư và các đơn vị kiểm định công trình cần phối hợp với các viện nghiên cứu chuyên ngành để thiết kế mật độ cảm biến gia tốc tại các nút khung chịu lực xung yếu, hướng đến mục tiêu giảm 30% chi phí lắp đặt thiết bị nhưng vẫn thu thập đủ 8 dạng dao động cơ bản, thực hiện trong khung thời gian 3 tháng trước khi nghiệm thu công trình.
  3. Xây dựng ngưỡng cảnh báo hư hỏng theo cấp độ: Các kỹ sư trưởng kết cấu và cơ quan quản lý kỹ thuật cần thiết lập hệ thống ngưỡng chỉ số MSEBI phân tầng (nguy cơ thấp dưới 10%, cảnh báo từ 10% đến 30%, và nguy hiểm trên 30% độ giảm độ cứng) nhằm chủ động kích hoạt quy trình bảo trì trong vòng 48 giờ sau cảnh báo, áp dụng liên tục trong suốt chu kỳ vận hành 10 đến 20 năm của công trình.
  4. Hiệu chuẩn mô hình số định kỳ: Đơn vị quản lý khai thác tòa nhà cần thực hiện đo đạc thực nghiệm và cập nhật lại ma trận độ cứng mô hình phần tử hữu hạn theo chu kỳ 12 tháng một lần, đảm bảo độ sai lệch tần số dao động riêng giữa mô hình và thực tế luôn duy trì dưới mức 2%.

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

Luận văn là nguồn tư liệu chuyên sâu có giá trị ứng dụng thực tiễn cao cho 4 nhóm đối tượng chính:

  1. Kỹ sư thiết kế và thẩm tra kết cấu công trình: Tiếp cận phương pháp phân tích dao động nâng cao trên mô hình không gian ba chiều, ứng dụng ma trận độ cứng tương thích để đánh giá sức chịu tải dư thừa và độ nhạy của khung nhà nhiều tầng trước các rủi ro nứt gãy vật liệu.
  2. Đơn vị tư vấn kiểm định và quan trắc công trình: Sở hữu quy trình hai bước hoàn chỉnh để chẩn đoán nhanh tình trạng hư hỏng cục bộ trên các công trình cầu, tháp thép và nhà cao tầng mà không cần can thiệp phá hủy hay làm gián đoạn khai thác.
  3. Giảng viên, nghiên cứu sinh và học viên cao học ngành Kỹ thuật Xây dựng: Sử dụng mã nguồn thuật toán, hệ thống công thức giải tích năng lượng biến dạng và thuật toán di truyền trong luận văn làm tài liệu giảng dạy, nghiên cứu mở rộng cho các bài toán tối ưu hóa kết cấu phức tạp.
  4. Chủ đầu tư và doanh nghiệp quản lý vận hành bất động sản: Nắm bắt cơ sở khoa học để xây dựng chiến lược bảo trì dự đoán, giúp tiết kiệm từ 20% đến 40% chi phí sửa chữa định kỳ và kéo dài tuổi thọ khai thác công trình thêm 15 đến 25 năm.

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

  1. Phương pháp năng lượng biến dạng có ưu điểm gì vượt trội so với phân tích tần số riêng truyền thống? Tần số dao động riêng là đặc trưng tổng thể nên kém nhạy với các hư hỏng cục bộ quy mô nhỏ. Trong khi đó, năng lượng biến dạng kết hợp trực tiếp véc tơ chuyển vị dạng dao động với ma trận độ cứng phần tử, giúp khuếch đại mức độ thay đổi cục bộ tại vị trí nứt, cho phép phát hiện hư hỏng ngay cả khi độ cứng chỉ suy giảm khoảng 10%.

  2. Tại sao quy trình chẩn đoán lại cần chia thành hai giai đoạn riêng biệt thay vì một giai đoạn? Nếu tối ưu hóa toàn bộ 63 phần tử cùng lúc trong một giai đoạn, thuật toán di truyền sẽ đối mặt với không gian tìm kiếm khổng lồ, dễ rơi vào bẫy cực trị địa phương. Quy trình hai giai đoạn dùng chỉ số MSEBI để lọc trước các vị trí nghi vấn, giúp giảm hơn 80% số biến số tối ưu và tăng độ chính xác lên trên 98%.

  3. Thuật toán di truyền xử lý bài toán xác định mức độ hư hỏng như thế nào? Thuật toán xem mức độ giảm độ cứng của các phần tử nghi vấn là các gen trong nhiễm sắc thể. Thông qua quá trình chọn lọc, lai tạo và đột biến qua khoảng 50 đến 100 thế hệ, thuật toán liên tục tối thiểu hóa hàm mục tiêu sai lệch năng lượng biến dạng OF2 để tìm ra véc tơ suy giảm độ cứng tối ưu nhất với sai số dưới 2%.

  4. Ảnh hưởng của nhiễu đo lường 3% tác động như thế nào đến độ chính xác của kết quả? Nhiễu đo lường thực tế 3% làm tăng nhẹ độ phân tán của véc tơ chuyển vị nút. Tuy nhiên, nhờ sử dụng tích hợp từ 6 đến 8 dạng dao động đầu tiên, chỉ số MSEBI vẫn triệt tiêu hiệu quả các thành phần nhiễu ngẫu nhiên, đảm bảo xác định đúng 100% vị trí hư hỏng chính với mức độ sai lệch chỉ tăng thêm khoảng 3% đến 5%.

  5. Phương pháp này có thể ứng dụng cho những loại kết cấu công trình nào trong thực tế? Phương pháp được thiết kế tối ưu cho tất cả các hệ kết cấu dạng thanh và khung không gian ba chiều như nhà khung thép tiền chế, kết cấu giàn khoan ngoài khơi, cầu giàn thép và tháp truyền hình. Các công trình này có đặc thù cấu tạo từ các phân đoạn thanh liên kết nút, hoàn toàn tương thích với mô hình phần tử hữu hạn 12 bậc tự do.

Kết luận

Luận văn thạc sĩ đã giải quyết thành công bài toán chẩn đoán hư hỏng kết cấu khung không gian ba chiều với 5 đóng góp khoa học nổi bật:

  • Xây dựng hoàn chỉnh khung lý thuyết và giải thuật chẩn đoán hai giai đoạn kết hợp năng lượng biến dạng và thuật toán di truyền.
  • Chứng minh chỉ số MSEBI vượt trội nhất trong việc nhận diện vị trí hư hỏng với độ chính xác đạt trên 95%.
  • Xác nhận hàm mục tiêu năng lượng biến dạng OF2 giúp thuật toán di truyền xác định mức độ hư hại với sai số tối thiểu dưới 2%.
  • Đánh giá thành công khả năng phát hiện hư hỏng cục bộ trên từng phân đoạn thanh 12 bậc tự do.
  • Kiểm chứng độ bền vững của thuật toán dưới tác động của nhiễu đo lường 3%.

Về lộ trình phát triển trong 12 đến 24 tháng tới, hướng nghiên cứu cần được mở rộng thử nghiệm trên mô hình thực nghiệm rung động ngoài hiện trường và tích hợp trí tuệ nhân tạo học sâu. Các kỹ sư và chuyên gia quan tâm hãy chủ động áp dụng quy trình tính toán tiên tiến này để nâng cao chất lượng kiểm định và an toàn công trình.