Luận án: Phương Pháp Phổ Tần Số Nghiên Cứu Dao Động Dầm Đàn Hồi Có Vết Nứt Chịu Tải

Bài toán tải trọng di động xuất hiện từ thế kỷ 19 gắn liền với sự phát triển đường sắt. Ban đầu, các kỹ sư tập trung vào hiện tượng tăng chuyển vị động của cầu khi tàu hỏa chạy qua. Tải trọng di động tạo ra dao động uốn phức tạp trên dầm. Vận tốc xe chạy ảnh hưởng trực tiếp đến biên độ dao động của kết cấu. Khi vận tốc tiến gần vận tốc tới hạn, dầm xảy ra hiện tượng cộng hưởng động lực. Hiện nay bài toán mở rộng cho kết cấu có khuyết tật cục bộ. Khảo sát đáp ứng động học dầm chịu tải mang tính cấp thiết cao. Sự xuất hiện của vết nứt làm thay đổi quy luật dao động uốn. Dao động của dầm không còn mang tính thuần nhất. Nghiên cứu giúp đảm bảo an toàn giao thông và kéo dài tuổi thọ công trình.

Nhiều phương pháp giải tích và số học được phát triển để giải quyết bài toán dao động dầm. Phương pháp Bubnov-Galerkin đưa phương trình vi phân đạo hàm riêng về hệ phương trình vi phân thường. Phương pháp phần tử hữu hạn chia dầm thành nhiều phần tử nhỏ để giải số. Phương pháp độ cứng động mang lại kết quả chính xác cao trong miền tần số. Mỗi phương pháp đều có ưu điểm và hạn chế riêng khi xét tải trọng di động. Khi dầm xuất hiện vết nứt, phương pháp phần tử hữu hạn đòi hỏi lưới phần tử rất mịn quanh đáy nứt. Điều này làm tăng khối lượng tính toán và thời gian xử lý. Phương pháp giải tích kết hợp hàm suy rộng giúp mô tả vết nứt đơn giản hơn. Việc áp dụng biến đổi Fourier giúp giảm thiểu đáng kể chi phí tính toán đối với hệ phức tạp.

Phân tích phổ tần số đóng vai trò trung tâm trong kiểm tra không phá hủy. Tín hiệu dao động đo được tại một số điểm trên dầm chứa đựng toàn bộ thông tin động học. Phép phân tích phổ cho phép trích xuất các đỉnh tần số tự nhiên và tần số kích động. Sự xuất hiện của vết nứt gây dịch chuyển các đỉnh cộng hưởng trên trục tần số. Biên độ phổ tại các tần số riêng bị suy giảm hoặc mở rộng do cản dao động. Việc theo dõi sự biến thiên này giúp chẩn đoán hư hỏng kết cấu qua dao động một cách tin cậy. Dữ liệu phổ tần số cũng ít bị ảnh hưởng bởi điều kiện đo tức thời so với dữ liệu miền thời gian. Đây là cơ sở cốt lõi để xây dựng hệ thống giám sát sức khỏe công trình tự động.

Phép biến đổi Fourier thuận đưa hàm chuyển vị thời gian về hàm phổ tần số phức. Biến đổi Fourier ngược cho phép khôi phục chuyển vị động tức thời. Hàm đáp ứng tần số thể hiện độ nhạy cảm của kết cấu tại từng tần số kích động riêng biệt. Khi tải trọng điều hòa di động qua dầm, hàm phổ xuất hiện hai đỉnh cực đại rõ rệt. Một đỉnh đại diện cho tần số dao động riêng của dầm. Đỉnh còn lại tương ứng với tần số kích động của lực di động. Hàm truyền đạt dao động phản ánh trung thực đặc trưng độ cứng và khối lượng của dầm. Việc tính toán hàm phổ qua tích phân Fourier giải tích giúp triệt tiêu sai số làm tròn số học. Kết quả phân tích phổ cung cấp bức tranh tổng thể về trạng thái động lực học.

Tần số dao động riêng của dầm nứt luôn nhỏ hơn tần số riêng của dầm nguyên vẹn. Mức độ suy giảm tần số phụ thuộc trực tiếp vào độ sâu và vị trí của vết nứt. Vết nứt nằm tại vị trí có mô men uốn lớn sẽ gây suy giảm tần số mạnh nhất. Ngược lại, vết nứt tại điểm nút dao động ít gây ảnh hưởng đến tần số bậc đó. Dạng dao động dầm đàn hồi bị uốn gãy nhẹ ngay tại vị trí có vết nứt. Sự không liên tục của đạo hàm bậc một của dạng riêng phản ánh góc xoay cục bộ. Khi kết cấu xuất hiện nhiều vết nứt, phổ tần số riêng thay đổi theo quy luật đa biến phức tạp. Khảo sát chính xác tần số riêng giúp nhận diện mức độ suy thoái kết cấu.

Bài toán ngược nhận dạng khuyết tật từ dữ liệu dao động thường không ổn định. Hiện tượng bài toán đặt không chỉnh xuất hiện do sai số đo lường và nhiễu môi trường. Ma trận hệ phương trình nhận dạng có chỉ số điều kiện rất lớn. Phương pháp điều chỉnh Tikhonov giải quyết triệt để sự mất ổn định này. Thuật toán bổ sung hàm mục tiêu ổn định hóa có trọng số điều chỉnh. Tham số điều chỉnh được lựa chọn tối ưu theo quy tắc độ lệch hoặc chuẩn hình L. Phương pháp triệt tiêu ảnh hưởng của nhiễu ngẫu nhiên trong tín hiệu phổ đo được. Nhờ đó, nghiệm nhận dạng hội tụ ổn định về vị trí và kích thước thực tế của vết nứt. Độ chính xác của thuật toán được bảo đảm ngay cả khi mức nhiễu lên đến mười phần trăm.

Mô hình dầm Euler-Bernoulli có vết nứt phù hợp cho dầm mảnh có tỷ lệ chiều dài trên chiều cao lớn. Mô hình này bỏ qua biến dạng trượt và quán tính quay của mặt cắt ngang. Khi tỷ lệ chiều dài trên chiều cao giảm, hiệu ứng trượt trở nên đáng kể. Lúc này, mô hình dầm Timoshenko có vết nứt được áp dụng để tăng độ chính xác. Lý thuyết Timoshenko bổ sung góc trượt cắt và quán tính quay vào phương trình vi phân chuyển động. Tần số riêng tính theo lý thuyết Timoshenko thường thấp hơn so với Euler-Bernoulli ở các bậc cao. Cả hai mô hình đều có thể tích hợp vết nứt qua hàm suy rộng Dirac. Việc lựa chọn mô hình phù hợp giúp cân bằng giữa độ chính xác và khối lượng tính toán.

Vết nứt ngang một phía được mô hình hóa bằng mô hình lò xo xoay tương đương không khối lượng. Độ cứng của lò xo được xác định thông qua hàm mật độ năng lượng biến dạng từ cơ học phá hủy. Độ sâu vết nứt càng lớn thì độ cứng của lò xo càng suy giảm. Ma trận độ cứng dầm nứt được suy dẫn từ điều kiện biên và phương trình tương thích tại vị trí lò xo. Chuyển vị võng vẫn liên tục qua vết nứt nhưng đạo hàm độ dốc bị nhảy bậc. Ma trận chuyển trạng thái kết nối các đoạn dầm nguyên vẹn liền kề nhau. Mô hình lò xo xoay cho phép biểu diễn dầm có nhiều vết nứt một cách thuận tiện. Công thức giải tích dạng đóng giúp việc tính toán phổ tần số đạt tốc độ xử lý vượt trội.

Vết nứt hở và vết nứt thở (breathing crack) có đặc tính cơ học hoàn toàn khác biệt. Vết nứt hở luôn mở trong suốt quá trình dao động dưới tải trọng. Độ cứng kết cấu giữ nguyên giá trị suy giảm không đổi theo thời gian. Ngược lại, vết nứt thở đóng và mở tuần hoàn theo chiều uốn của dầm. Khi vết nứt đóng lại, dầm khôi phục độ cứng gần như ban đầu. Hiện tượng đóng mở tạo ra tính phi tuyến mạnh mẽ trong phương trình dao động. Phổ tần số của dầm có vết nứt thở xuất hiện các hài bậc cao và hài phân kỳ. Việc nhận biết vết nứt thở đòi hỏi phân tích phổ phổ động bậc cao hoặc đáp ứng phi tuyến. Việc mô hình hóa vết nứt hở đóng vai trò cơ sở ban đầu trước khi xét đến hiệu ứng phi tuyến phức tạp.

Đáp ứng động học dầm chịu tải trọng điều hòa di động biểu hiện sự giao thoa sóng phức tạp. Lực kích động có biên độ và tần số kích thích dao động cưỡng bức liên tục. Trên phổ biên độ đáp ứng xuất hiện các đỉnh phổ kép. Đỉnh thứ nhất liên quan đến vận tốc di chuyển của lực dọc theo nhịp dầm. Đỉnh thứ hai phụ thuộc vào tần số dao động riêng của dầm nứt. Biên độ đáp ứng tăng đột biến khi tần số lực kích động trùng với tần số riêng kết cấu. Hiện tượng cộng hưởng kép xảy ra khi cả vận tốc chạy xe và tần số kích động đạt giá trị tới hạn. Khảo sát phổ tần số giúp kỹ sư dự báo và phòng ngừa các tình huống nguy hiểm này.

Dạng dao động dầm đàn hồi chịu ảnh hưởng sâu sắc bởi sự phân bố của nhiều vết nứt cùng lúc. Mỗi vết nứt tạo ra một điểm suy giảm độ dốc cục bộ trên đường cong dạng riêng. Dạng dao động bậc cao nhạy cảm hơn với các vết nứt phân tán dọc nhịp dầm. Khi số lượng vết nứt tăng, các nút dao động có xu hướng dịch chuyển về phía các phân đoạn yếu. Phép phân tích phổ tổng quát cho phép xác định đồng thời ảnh hưởng của nhiều vị trí hư hỏng. Ma trận chuyển trạng thái nối kết giúp tính toán dạng riêng của dầm liên tục nhiều nhịp. Độ nhạy của các dạng riêng cung cấp thông tin không gian quý giá cho bài toán chẩn đoán.

Vận tốc di chuyển của tải trọng quyết định thời gian tác dụng của lực lên dầm đàn hồi. Khi vận tốc tăng, tần số kích thích kéo theo cũng tăng tương ứng. Phổ biên độ bị dịch chuyển và kéo giãn dọc theo trục tần số. Hệ số cản dao động Rayleigh đóng vai trò dập tắt các thành phần dao động quá độ. Tỷ số cản cao làm giảm chiều cao đỉnh phổ cộng hưởng nhưng mở rộng bề rộng đỉnh. Độ sâu vết nứt kết hợp với vận tốc lớn gây ra chuyển vị động vượt mức giới hạn an toàn. Đánh giá đồng thời vận tốc và hệ số cản giúp tối ưu hóa giới hạn tốc độ lưu thông trên cầu.

Chẩn đoán hư hỏng kết cấu qua dao động là phương pháp kiểm tra không phá hủy hiện đại. Bằng cách gắn cảm biến gia tốc trên dầm, toàn bộ lịch sử dao động được ghi nhận. Dữ liệu thời gian được chuyển hóa thành phổ tần số nhờ thuật toán biến đổi Fourier nhanh. Các chỉ số nhạy cảm như hàm truyền và mật độ phổ công suất được trích xuất. So sánh phổ trước và sau khi khai thác giúp xác định kết cấu đã suy giảm chất lượng hay chưa. Kỹ thuật này phát hiện được cả những vết nứt ngầm nằm sâu bên trong thân dầm bê tông cốt thép hoặc thép. Quá trình kiểm tra diễn ra liên tục mà không làm gián đoạn giao thông qua cầu.

Nhận dạng vị trí và độ sâu vết nứt là mục tiêu cốt lõi của bài toán ngược động lực học. Thuật toán sử dụng vector sai lệch phổ tần số đo đạc và phổ lý thuyết. Vị trí vết nứt được định vị thông qua điểm kỳ dị hoặc cực trị của hàm chỉ thị hư hỏng. Độ sâu vết nứt được suy ra từ độ suy giảm tương đối của ma trận độ cứng dầm nứt. Thuật toán tối ưu hóa lặp giúp hiệu chỉnh thông số mô hình hội tụ về giá trị thực tế. Kết quả thử nghiệm số chứng minh thuật toán nhận dạng chính xác đồng thời hai đến ba vết nứt trên cùng một nhịp dầm. Sai số xác định vị trí dưới năm phần trăm và sai số độ sâu dưới tám phần trăm.

Nhiễu đo lường từ môi trường và cảm biến luôn hiện diện trong dữ liệu thực tế. Tín hiệu dao động bị pha trộn bởi nhiễu trắng ngẫu nhiên với các mức biên độ khác nhau. Nhờ tích hợp thuật toán điều chỉnh Tikhonov, độ ổn định của lời giải nhận dạng được duy trì vững chắc. Khi mức nhiễu đạt năm phần trăm hoặc mười phần trăm, vị trí vết nứt vẫn được định vị chính xác. Độ sâu vết nứt có độ lệch nhẹ nhưng vẫn nằm trong phạm vi kỹ thuật cho phép. Sự kết hợp giữa bộ lọc tín hiệu phổ và thuật toán chính quy hóa mang lại độ tin cậy vượt trội. Phương pháp hoàn toàn khả thi để ứng dụng vào hệ thống giám sát sức khỏe công trình thực tế.