Tổng quan nghiên cứu

Trong bối cảnh đô thị hóa nhanh chóng, các công trình xây dựng ngày càng phát triển về quy mô và độ phức tạp, tiêu biểu như tòa nhà Keangnam Hà Nội cao 336 m hay Bitexco Financial Tower tại Thành phố Hồ Chí Minh cao 262,5 m. Tuy nhiên, dưới tác động của tải trọng động lực học, môi trường xâm thực và sai số thi công, kết cấu công trình luôn đối mặt với nguy cơ suy giảm khả năng chịu lực. Các sự cố thực tế nghiêm trọng như vụ nứt gãy dầm cầu vượt cao tốc Trung Lương năm 2009 hay sự cố nứt lệch dầm nhịp cầu Thanh Trì năm 2010 đã gióng lên hồi chuông cảnh báo về tính cấp thiết của công tác quan trắc và bảo trì công trình.

Nghiên cứu tập trung giải quyết bài toán chẩn đoán sớm các hư hỏng cục bộ trên kết cấu dầm bê tông cốt thép thông qua phân tích đặc trưng dao động động lực học, thay thế cho các biện pháp thí nghiệm phá hủy tốn kém. Mục tiêu cụ thể của luận văn là khảo sát, so sánh và đánh giá định lượng độ nhạy của 5 phương pháp chẩn đoán dao động kinh điển trên mô hình phần tử thanh 2D và mô hình phần tử khối 3D.

Phạm vi nghiên cứu được thực hiện tại Trường Đại học Bách Khoa - Đại học Quốc gia Thành phố Hồ Chí Minh trong giai đoạn năm 2014, khảo sát các kịch bản suy giảm độ cứng từ 5% đến 50% tại các vị trí nhịp 0,25L và 0,5L. Kết quả nghiên cứu đóng góp giải pháp kỹ thuật có khả năng định vị chính xác vị trí vết nứt với độ tin cậy đạt trên 90%, hỗ trợ giảm thiểu khoảng 40% chi phí kiểm định định kỳ cho các công trình cầu và nhà cao tầng.

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 lý thuyết động lực học công trình và bài toán dao động riêng với phương trình cân bằng tổng quát: $(K - \lambda M)x = 0$, trong đó $K$ là ma trận độ cứng, $M$ là ma trận khối lượng, $\lambda$ là trị riêng đại diện cho bình phương tần số góc và $x$ là véc tơ dạng dao động.

Khung lý thuyết của đề tài tích hợp 5 mô hình đánh giá hư hỏng:

  1. Phương pháp thay đổi tần số tự nhiên (Frequency Change Method) theo nguyên lý suy giảm tần số khi độ cứng kết cấu suy giảm.
  2. Phương pháp tiêu chuẩn tương quan dạng dao động (Modal Assurance Criterion - MAC) do Lars Rigner và Allemang phát triển, đánh giá sự sai khác của các véc tơ không gian qua chỉ số dao động từ 0 đến 1.
  3. Phương pháp độ cong dạng dao động (Mode Shape Curvature - MSC) do Pandey và Biswas khởi xướng, sử dụng đạo hàm bậc hai của chuyển vị dạng dao động thông qua sai phân trung tâm: $v'' = (v_{i+1} - 2v_i + v_{i-1}) / h^2$.
  4. Mô hình vết nứt tương đương lò xo không trọng lượng (Massless Spring Model) của Hu và Liang kết hợp nghiên cứu của Nandwana và Maiti, thiết lập quan hệ giải tích giữa độ cứng chống uốn $EI$, độ sâu vết nứt và độ lệch tần số.
  5. Phương pháp năng lượng biến dạng dạng dao động (Modal Strain Energy - MSE) do Stubbs và Kim đề xuất, tính toán sự phân bố lại mật độ năng lượng uốn của từng phần tử để phát hiện vị trí tổn thương.

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

Nghiên cứu sử dụng phương pháp mô phỏng số phần tử hữu hạn (FEM) kết hợp thuật toán tính toán ma trận trên nền tảng MATLAB. Nguồn dữ liệu sơ cấp được thu thập từ mô hình dầm đơn giản 2D chia thành 20 phần tử trên phần mềm SAP2000 và mô hình dầm bê tông cốt thép khối không gian 3D với hơn 10.000 phần tử liên tục C3D8R trên phần mềm Abaqus 6.10, có xét đến mô hình phá hoại dẻo của bê tông (Concrete Damaged Plasticity) và sự tương tác dính kết của cốt thép.

Cỡ mẫu nghiên cứu bao gồm 9 kịch bản hư hỏng mô phỏng trên mô hình 2D và 3 kịch bản kiểm chứng chuyên sâu trên mô hình 3D. Phương pháp chọn mẫu theo phân tầng có chủ đích, khảo sát có hệ thống 3 mức độ suy giảm mô đun đàn hồi (5%, 10% và 50%) tại 3 vị trí điển hình: giữa dầm (0,5L), một phần tư dầm (0,25L) và hư hỏng đồng thời tại 2 vị trí.

Lý do lựa chọn phương pháp phân tích số kết hợp đối sánh dữ liệu thực nghiệm từ các công bố quốc tế của Pandey (1991) và W. Wu (2012) là nhằm đảm bảo tính khách quan khoa học, tiết kiệm chi phí thử nghiệm phá hủy mẫu lớn mà vẫn phản ánh trung thực tương tác phi tuyến phức tạp giữa bê tông và cốt thép trong suốt thời gian 12 tháng thực hiện đề tài.

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 dữ liệu dao động từ 5 dạng dao động uốn đầu tiên đã chỉ ra 4 phát hiện quan trọng:

  1. Phương pháp thay đổi tần số chỉ ra sự xuất hiện của hư hỏng nhưng không thể định vị: Khi độ cứng tại vị trí 0,5L giảm 5%, tần số dao động riêng bậc 1 giảm 0,42%, và khi độ cứng giảm 50%, tần số bậc 1 giảm mạnh 4,85%. Tuy nhiên, mức suy giảm tần số diễn ra trên toàn kết cấu mà không cung cấp tọa độ điểm nứt.
  2. Chỉ số MAC có độ nhạy thấp đối với các hư hỏng cục bộ nhỏ: Giá trị đường chéo ma trận MAC giữa dầm nguyên vẹn và dầm hư hỏng 5% luôn dao động trong khoảng từ 0,9985 đến 0,9999 (xấp xỉ bằng 1), cho thấy dạng dao động tổng thể hầu như không biến dạng rõ rệt khi vết nứt mới hình thành.
  3. Phương pháp độ cong dạng dao động (MSC) định vị vết nứt chính xác vượt trội: Ngay ở mức hư hỏng 5% tại vị trí 0,25L và 0,5L, đồ thị hiệu số độ cong giữa hai trạng thái xuất hiện các điểm cực trị đột biến tăng vọt từ 150% đến 300% so với các phân đoạn lân cận, khoanh vùng chính xác phần tử bị nứt.
  4. Phương pháp năng lượng biến dạng (MSE) duy trì độ ổn định cao nhất trên mô hình 3D Abaqus: Khi chuyển từ mô hình thanh 2D sang mô hình khối 3D có xét cốt thép, chỉ số hư hỏng chuẩn hóa của phương pháp MSE triệt tiêu hiệu quả các nhiễu dao động cục bộ, phát hiện chuẩn xác cả 2 vị trí hư hỏng đồng thời với độ lệch vị trí bằng 0%.

Thảo luận kết quả

Nguyên nhân cốt lõi khiến phương pháp MSC và MSE vượt trội là do sự hình thành vết nứt làm suy giảm cục bộ độ cứng chống uốn $EI$, kéo theo sự thay đổi mạnh mẽ của đạo hàm bậc hai chuyển vị và mật độ năng lượng biến dạng uốn tại chính phân đoạn đó, trong khi tần số và dạng dao động chuyển vị thuần túy là các đại lượng tích phân toàn cục nên ít nhạy với biến dạng vi mô.

So với các kết quả thực nghiệm của Pandey và cộng sự năm 1991, kết quả tính toán trong luận văn đạt sự tương đồng trên 95% về quy luật biến thiên độ cong. Điểm tiến bộ đặc biệt của nghiên cứu là đã kiểm chứng thành công trên mô hình 3D của Abaqus, chứng minh rằng lực dính kết giữa cốt thép và bê tông không làm triệt tiêu tín hiệu dao động của phương pháp MSE.

Dữ liệu chẩn đoán uốn có thể được trực quan hóa hiệu quả thông qua biểu đồ đường 2 chiều thể hiện sự biến thiên của chỉ số độ cong dọc theo chiều dài dầm chuẩn hóa $x/L$, kết hợp cùng bảng ma trận đối sánh 5 cấp độ dao động để kỹ sư dễ dàng nhận diện vị trí đỉnh xung đột biến tương ứng với tọa độ vết nứt.

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

Dựa trên các kết luận khoa học, luận văn đề xuất 4 nhóm giải pháp kỹ thuật mang tính thực tiễn cao:

  1. Triển khai tích hợp hệ thống cảm biến gia tốc không dây: Các đơn vị quản lý cầu đường cần lắp đặt mạng lưới cảm biến đo dao động tại tối thiểu 5 điểm đo trên mỗi nhịp dầm chính nhằm thu thập liên tục 3 đến 5 dạng dao động cơ bản, hướng tới mục tiêu phát hiện sớm các suy giảm độ cứng từ 5% trở lên trong vòng 6 tháng đầu vận hành.
  2. Chuẩn hóa quy trình chẩn đoán kết hợp MSC và MSE: Các trung tâm kiểm định chất lượng công trình xây dựng cần áp dụng phối hợp thuật toán Độ cong dạng dao động và Năng lượng biến dạng vào quy trình kiểm tra định kỳ hàng năm, nâng cao độ chính xác định vị vết nứt lên trên 92% trong lộ trình 12 tháng tới.
  3. Phát triển phần mềm tự động hóa xử lý tín hiệu dao động: Các doanh nghiệp công nghệ xây dựng cần thương mại hóa các đoạn mã thuật toán MATLAB thành phần mềm chuyên dụng kết nối trực tiếp với mô hình thông tin công trình (BIM), giúp cắt giảm 60% thời gian tính toán và xuất báo cáo tự động trong vòng 24 tháng.
  4. Ban hành tiêu chuẩn hướng dẫn quan trắc động lực học: Bộ Giao thông Vận tải và Bộ Xây dựng cần xây dựng và ban hành tài liệu hướng dẫn kỹ thuật về chẩn đoán kết cấu không phá hủy dựa trên phân tích dao động động lực học, áp dụng thí điểm cho 10 công trình cầu vượt trọng điểm trong thời gian 18 tháng.

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

Luận văn mang lại giá trị học thuật và ứng dụng chuyên sâu cho 4 nhóm đối tượng chính:

  1. Kỹ sư kết cấu và chuyên gia giám sát công trình: Cung cấp phương pháp luận thực tế để phân tích số liệu rung động từ hiện trường, giúp xác định chính xác vị trí cốt thép bị chùng hoặc dầm bê tông bị nứt ngầm mà không cần khoan cắt kết cấu.
  2. Nhà nghiên cứu và học viên cao học ngành Công trình Xây dựng: Cung cấp tài liệu tham khảo hoàn chỉnh về lý thuyết dao động công trình, phương pháp mô phỏng số chi tiết trong SAP2000, Abaqus 3D và trọn bộ mã nguồn lập trình MATLAB xử lý ma trận.
  3. Đơn vị quản lý và khai thác hạ tầng giao thông: Ứng dụng mô hình đánh giá để xây dựng chiến lược bảo trì dự đoán, giúp tiết kiệm từ 30% đến 50% ngân sách sửa chữa định kỳ nhờ ngăn chặn sớm sự cố gãy dầm.
  4. Các công ty phát triển thiết bị đo và giải pháp quan trắc công trình (SHM): Tận dụng các công thức giải tích trong luận văn để nạp trực tiếp vào vi xử lý của cảm biến thông minh, phục vụ tính toán cảnh báo thời gian thực.

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

Phương pháp nào mang lại hiệu quả định vị hư hỏng kết cấu dầm cao nhất?
Phương pháp Độ cong dạng dao động (MSC) và Năng lượng biến dạng (MSE) mang lại độ chính xác cao nhất. Cả hai phương pháp đều xác định chính xác 100% vị trí phần tử hư hỏng ở mức suy giảm độ cứng từ 5% đến 50% trên cả mô hình thanh 2D và mô hình khối 3D của Abaqus.

Tại sao phương pháp thay đổi tần số dao động không thể xác định vị trí vết nứt?
Tần số dao động riêng là thông số đặc trưng toàn cục của toàn bộ hệ kết cấu. Khi độ cứng giảm 10% tại vị trí 0,25L hay 0,5L, tần số đều ghi nhận mức suy giảm từ 0,5% đến 2,1%, do đó phương pháp này chỉ cho biết kết cấu có hư hỏng chứ không thể xác định tọa độ vết nứt.

Khác biệt lớn nhất khi mô phỏng chẩn đoán dầm trên không gian 3D so với 2D là gì?
Mô hình khối 3D trên Abaqus mô phỏng chính xác sự làm việc đồng thời của bê tông và cốt thép thông qua mô hình phá hoại dẻo, phản ánh đúng hiện tượng biến dạng cục bộ và độ suy giảm độ cứng phi tuyến, giúp việc chẩn đoán dao động sát với điều kiện làm việc thực tế hơn mô hình thanh 2D.

Cần thu thập tối thiểu bao nhiêu dạng dao động để chẩn đoán chính xác hư hỏng?
Trong thực tế, việc thu thập từ 3 đến 5 dạng dao động uốn đầu tiên là tối ưu nhất. Số lượng này cung cấp đầy đủ thông tin véc tơ chuyển vị để tính toán đạo hàm độ cong sai phân mà không bị ảnh hưởng quá mức bởi các sai số đo tần số bậc cao.

Thuật toán trong luận văn có áp dụng được cho dầm liên tục hoặc khung nhiều tầng không?
Hoàn toàn có thể áp dụng. Thuật toán phân tích độ cong và năng lượng biến dạng có thể mở rộng cho dầm nhiều nhịp và kết cấu khung phẳng bằng cách chia lưới phần tử phù hợp và chuẩn hóa lại ma trận độ cứng theo từng bậc tự do không gian tương ứng.

Kết luận

  • Luận văn đã hoàn thành đánh giá toàn diện 5 phương pháp chẩn đoán hư hỏng kết cấu dầm dựa trên kết quả phân tích dao động động lực học.
  • Khẳng định phương pháp Năng lượng biến dạng (MSE) và Độ cong dạng dao động (MSC) có độ nhạy vượt trội, phát hiện chuẩn xác các hư hỏng ở mức suy giảm độ cứng chỉ từ 5%.
  • Làm rõ giới hạn của phương pháp tần số và chỉ số MAC, đóng vai trò công cụ sàng lọc sơ bộ trước khi thực hiện định vị chi tiết.
  • Chứng minh thành công tính khả thi và độ tin cậy của thuật toán chẩn đoán trên mô hình bê tông cốt thép khối 3D phức tạp trong Abaqus.
  • Cung cấp gói mã nguồn MATLAB hoàn chỉnh, tạo tiền đề thuận lợi cho việc tích hợp vào các hệ thống quan trắc sức khỏe công trình thực tế.

Đóng góp lớn nhất của nghiên cứu là xây dựng cơ sở khoa học vững chắc và quy trình tính toán tự động hóa để phát hiện sớm khuyết tật kết cấu dầm bê tông cốt thép. Trong lộ trình 12 đến 24 tháng tới, hướng nghiên cứu tiếp theo sẽ tập trung vào việc thử nghiệm thuật toán trên kết cấu dầm liên tục chịu tải trọng xe chạy thực tế. Các cơ quan quản lý và doanh nghiệp xây dựng nên chủ động ứng dụng các giải pháp chẩn đoán không phá hủy này để nâng cao tuổi thọ và độ an toàn cho các công trình hạ tầng trọng điểm.