Tổng quan nghiên cứu

Trong các dây chuyền sản xuất công nghiệp hiện đại, hộp giảm tốc và bộ truyền bánh răng đóng vai trò truyền động then chốt. Thống kê kỹ thuật chỉ ra rằng hơn 60% các sự cố cơ khí trong thiết bị quay xuất phát từ hư hỏng bánh răng và ổ lăn, gây ra thiệt hại do dừng máy ước tính lên tới 15.000 USD mỗi giờ trong các ngành công nghiệp nặng. Vấn đề nghiên cứu trọng tâm là việc chẩn đoán sớm các hư hỏng cục bộ như nứt chân răng, tróc rỗ bề mặt hay mẻ răng gặp rất nhiều thách thức do tín hiệu rung động bị che lấp bởi nhiễu nền và hiện tượng điều chế biên độ, pha phức tạp. Các phương pháp giải mã xung hiện đại như giải mã cực đại kurtosis tương quan (MCKD) hay giải mã cực tiểu Entropy tối ưu đa điểm (MOMEDA) tuy cải thiện khả năng trích xuất xung tuần hoàn nhưng bắt buộc phải biết trước chu kỳ lỗi chính xác. Khi vận tốc quay dao động nhỏ dẫn đến chu kỳ lỗi thay đổi khoảng 5%, các phương pháp này lập tức giảm độ chính xác nghiêm trọng và không thể giải quyết bài toán lỗi kết hợp nhiều chu kỳ trong hộp số nhiều cấp.

Mục tiêu cụ thể của luận văn là nghiên cứu, xây dựng và hoàn thiện phương pháp giải mã cực tiểu Entropy tối ưu đã hiệu chỉnh (Modified OMEDA) kết hợp công cụ trung bình miền bao nhằm trích xuất toàn bộ các xung lỗi chu kỳ và ngẫu nhiên mà không cần khai báo trước chu kỳ hư hỏng. Đề tài được thực hiện từ ngày 14 tháng 02 năm 2022 đến ngày 10 tháng 12 năm 2022 tại Khoa Cơ khí, Trường Đại học Bách Khoa, Đại học Quốc gia Thành phố Hồ Chí Minh. Nghiên cứu tập trung giới hạn vào các dạng lỗi cục bộ của bộ truyền bánh răng trên cơ sở phân tích tín hiệu rung động thu thập từ cảm biến gia tốc. Ý nghĩa học thuật và thực tiễn của công trình thể hiện qua việc nâng cao độ chính xác chẩn đoán lỗi lên trên 95%, giảm thiểu hơn 40% thời gian thiết lập thông số ban đầu và cung cấp giải pháp xử lý tín hiệu vượt trội cho các hệ thống giám sát tình trạng máy móc thời gian thực.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Cơ sở lý thuyết của luận văn được xây dựng dựa trên mô hình toán học về tín hiệu rung động ăn khớp bánh răng và lý thuyết giải mã tích chập khử nhiễu. Mô hình rung động tổng quát của bánh răng khi xuất hiện khuyết tật cục bộ được biểu diễn qua sự kết hợp giữa dao động ăn khớp điều chế biên độ, điều chế pha và phản ứng xung cộng hưởng kết cấu khi các bề mặt tổn thương va đập. Tín hiệu rung động đo được là tích chập của chuỗi xung kích thích hư hỏng với đáp ứng xung hữu hạn (FIR) của hệ thống truyền động, cộng thêm thành phần động học và nhiễu trắng Gaussian.

Nghiên cứu ứng dụng 4 khái niệm và chỉ số then chốt:

  1. Chỉ số Kurtosis: Đo lường độ nhọn của phân bố tín hiệu rung động. Tín hiệu chuẩn Gaussian có Kurtosis bằng 3, khi xuất hiện xung hư hỏng chỉ số này tăng lên từ 4 đến 6, và vượt quá 9 khi hư hỏng cơ khí ở mức nghiêm trọng.
  2. Chuẩn K-Norm và D-Norm: Trong khi D-Norm phản ánh tỷ lệ biên độ cực đại bậc một trên chuẩn Euclidean, K-Norm khai thác mô-men thống kê bậc bốn, mang lại độ nhạy vượt trội đối với các xung kích thích thoáng qua.
  3. Véc-tơ chỉ lỗi thích ứng: Cấu trúc véc-tơ mục tiêu nhị phân xác định vị trí các xung mục tiêu được tối ưu hóa động thay vì phân bổ cố định theo chu kỳ giả định.
  4. Trung bình miền bao: Kỹ thuật làm mịn đường bao biên độ nhằm triệt tiêu dao động ngẫu nhiên và làm nổi bật tần số lặp lại của khuyết tật.

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

Nguồn dữ liệu thực nghiệm được huy động từ 3 nguồn chính: mô hình mô phỏng số tín hiệu ăn khớp bánh răng chứa xung đơn và xung kết hợp hai chu kỳ (chu kỳ thứ nhất bằng 210 mẫu, chu kỳ thứ hai bằng 250 mẫu trong bản ghi 10.000 điểm dữ liệu), bộ dữ liệu chuẩn của tổ chức GPMS Inc trên tuabin gió, và bộ dữ liệu thực nghiệm của Đại học Cincinnati (UoC) với 9 chế độ hư hỏng bánh răng từ mức độ nhẹ đến phá hủy hoàn toàn.

Phương pháp chọn mẫu áp dụng kỹ thuật lấy mẫu phân tầng có chủ đích, bao phủ các điều kiện vận hành từ không tải, 10% tải đến 100% đầy tải, đồng thời thực hiện lấy mẫu lại theo góc quay thông qua kỹ thuật theo dõi bậc (Order Tracking) và trung bình tín hiệu đồng bộ (TSA). Lý do lựa chọn phương pháp phân tích Modified OMEDA là khả năng giải trực tiếp nghiệm bộ lọc thông qua nghịch đảo ma trận Toeplitz không lặp, loại bỏ nguy cơ không hội tụ của các giải thuật truyền thống. Chiều dài bộ lọc FIR được xác lập nghiêm ngặt trong khoảng từ 1% đến 4% độ dài bản ghi tín hiệu (tương đương 20 đến 100 mẫu) nhằm triệt tiêu hoàn toàn hiện tượng chồng lấn xung và ngăn ngừa sai số vị trí. Quy trình nghiên cứu được triển khai chặt chẽ qua 4 giai đoạn kéo dài 10 tháng, bao gồm xây dựng thuật toán, lập trình mô phỏng trên MATLAB, kiểm chứng dữ liệu thực nghiệm và đánh giá độ tin cậy.

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

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

Quá trình mô phỏng và thực nghiệm đã đem lại 4 phát hiện quan trọng:

  • Trích xuất chính xác 100% xung lỗi kết hợp không phụ thuộc chu kỳ: Trên tín hiệu mô phỏng chứa hai chu kỳ lỗi đồng thời (chu kỳ 210 mẫu và 250 mẫu), phương pháp Modified OMEDA đã phân lập hoàn hảo toàn bộ các xung kích thích. Trong khi đó, MOMEDA và MCKD hoàn toàn thất bại khi chu kỳ lỗi thực tế dao động lệch khoảng 5% so với giá trị giả định.
  • Tối ưu hóa dải kích thước bộ lọc triệt tiêu xung giả: Khi áp dụng chiều dài bộ lọc trong ngưỡng khuyến nghị từ 1% đến 4% chiều dài bản ghi, thuật toán đã loại bỏ hơn 85% các đỉnh xung giả phát sinh do nhiễu trắng, đồng thời khống chế sai số vị trí xung luôn nhỏ hơn chiều dài một bộ lọc.
  • Khuếch đại độ nhạy chẩn đoán qua trung bình miền bao: Kết hợp Modified OMEDA với trung bình miền bao giúp chỉ số Kurtosis của tín hiệu sau giải mã tăng từ 3,15 ở trạng thái khỏe mạnh lên mức 14,85 ở chế độ nứt chân răng và 12,60 ở chế độ tróc mặt răng, tạo khoảng cách nhận diện an toàn vượt bậc so với tín hiệu gốc.
  • Phân loại thành công 9 chế độ hư hỏng thực nghiệm: Đánh giá trên tập dữ liệu UoC cho thấy phương pháp nhận diện chính xác 100% các mức độ khuyết tật từ mẻ đầu răng nhẹ (giai đoạn TSA0) đến nứt chân răng sâu (giai đoạn TSA5), chứng minh tính ổn định tuyệt đối trước sự thay đổi tải trọng và vận tốc vận hành.

Thảo luận kết quả

Hiệu quả vượt trội của Modified OMEDA bắt nguồn từ việc thay thế chuẩn D-Norm bằng chuẩn K-Norm kết hợp thuật toán tìm kiếm đỉnh cục bộ thích ứng. K-Norm khuếch đại năng lượng xung bậc bốn, tạo ra sự chênh lệch biên độ rõ nét giữa vùng chứa xung thực và vùng nhiễu nền trong khoảng lân cận một chiều dài bộ lọc. Việc hiệu chỉnh vị trí véc-tơ lỗi trong phạm vi cục bộ giúp bộ lọc FIR thích ứng hoàn hảo với sự trượt pha và biến thiên vận tốc quay của trục máy.

So với phương pháp MCKD truyền thống đòi hỏi từ 15 đến 30 vòng lặp với nguy cơ rơi vào cực trị địa phương, Modified OMEDA chỉ thực hiện phép nghịch đảo ma trận giải tích một lần và tối ưu hóa vị trí xung cục bộ, giúp tiết kiệm khoảng 50% đến 60% thời gian tính toán CPU (tốc độ xử lý đạt 0,08 đến 0,15 giây cho một khung tín hiệu 5.000 mẫu so với 0,42 giây của thuật toán lặp).

Dữ liệu kết quả có thể được trực quan hóa tối ưu thông qua biểu đồ dạng sóng thời gian đa kênh so sánh tín hiệu trước và sau lọc, kết hợp bảng thống kê phân bố năng lượng miền bao. Trên biểu đồ miền thời gian, các đỉnh xung lỗi sau khi xử lý bằng Modified OMEDA hiển thị tách biệt với biên độ vượt trội gấp 4 đến 5 lần mức nhiễu nền, giúp kỹ sư quan sát trực tiếp chu kỳ va đập mà không cần qua các bước biến đổi Fourier phức tạp.

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

Nhằm chuyển giao kết quả nghiên cứu vào thực tiễn sản xuất, 4 giải pháp cụ thể được đề xuất:

  1. Tích hợp thuật toán Modified OMEDA vào hệ thống phần mềm giám sát rung động trực tuyến: Nhúng trực tiếp giải thuật vào các thiết bị đo biên (Edge Computing) và hệ sinh thái SCADA nhằm phát hiện sớm hư hỏng hộp số, hướng tới mục tiêu giảm 25% thời gian ngừng máy đột xuất trong vòng 6 tháng triển khai, do Phòng Bảo trì và Kỹ thuật số chủ trì thực hiện.
  2. Chuẩn hóa quy trình tiền xử lý tín hiệu bằng kỹ thuật theo dõi bậc và TSA: Bắt buộc áp dụng lấy mẫu lại đồng bộ theo tín hiệu cảm biến tốc độ (tachometer) trước khi đưa vào bộ giải mã xung nhằm triệt tiêu trên 80% ảnh hưởng của biến thiên vận tốc, hoàn thành tài liệu hướng dẫn trong 3 tháng, do Đội ngũ Kỹ sư Phân tích Rung động đảm nhiệm.
  3. Tự động hóa việc thiết lập chiều dài bộ lọc FIR thích ứng: Phát triển mô-đun tự động nhận diện độ dài bản ghi và giới hạn kích thước bộ lọc trong dải 1% đến 4% chiều dài mẫu, giúp giảm 90% các cảnh báo lỗi giả do thiết lập sai thông số, thực hiện trong 4 tháng bởi Nhóm Nghiên cứu và Phát triển Thuật toán (R&D).
  4. Xây dựng kho cơ sở dữ liệu đối chuẩn cho 9 dạng hư hỏng truyền động cơ khí: Số hóa và lưu trữ các mẫu phổ trung bình miền bao của các chế độ mẻ răng, nứt chân răng và mòn mặt răng để phục vụ đào tạo mô hình học máy (Machine Learning), nâng cao độ chính xác nhận dạng tự động đạt trên 95% trong vòng 12 tháng, do Ban Quản lý Thiết bị phối hợp cùng các viện nghiên cứu chuyên ngành triển khai.

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

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

  • Kỹ sư bảo trì dự đoán (PdM) và chuyên gia phân tích rung động: Nắm vững phương pháp trích xuất xung kích thích không cần giả định trước chu kỳ lỗi, ứng dụng trực tiếp để chẩn đoán chính xác các sự cố cục bộ trên hộp giảm tốc tuabin gió, máy nghiền xi măng và hộp số tàu thủy.
  • Học viên cao học, nghiên cứu sinh ngành Kỹ thuật Cơ khí và Kỹ thuật Điều khiển: Tiếp cận khung lý thuyết toán học chặt chẽ về giải mã tích chập khử nhiễu, mô hình rung động động lực học bánh răng và các phương pháp xử lý tín hiệu nâng cao phục vụ việc phát triển đề tài nghiên cứu liên quan.
  • Nhà phát triển giải pháp IoT công nghiệp (IIoT) và phần mềm chẩn đoán: Khai thác thuật toán giải tích không lặp với tốc độ tính toán nhanh để lập trình các khối chức năng phân tích rung động trên vi xử lý nhúng và nền tảng giám sát đám mây 24/7.
  • Cán bộ quản lý kỹ thuật và giám đốc vận hành nhà máy: Tham khảo cơ sở khoa học để xây dựng chiến lược bảo trì theo tình trạng máy (CBM), tối ưu hóa thời gian đại tu thiết bị và tiết kiệm hàng trăm triệu đồng chi phí thay thế phụ tùng hàng năm.

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

Sự khác biệt cốt lõi giữa Modified OMEDA và MOMEDA truyền thống là gì? MOMEDA bắt buộc phải khai báo chính xác chu kỳ lỗi và chỉ trích xuất một chu kỳ định sẵn. Ngược lại, Modified OMEDA sử dụng chuẩn K-Norm để tự động tìm kiếm các đỉnh xung thực tế, cho phép trích xuất đồng thời nhiều dạng xung lỗi có chu kỳ khác nhau hoặc biến thiên mà không cần biết trước thông số vận hành.

Tại sao chiều dài bộ lọc FIR cần được khống chế trong khoảng từ 1% đến 4% độ dài bản ghi? Nếu bộ lọc quá lớn (vượt quá khoảng cách giữa hai xung), các vùng tăng biên độ sẽ chồng lấn lên nhau gây sai lệch vị trí xung và xuất hiện xung giả. Giới hạn từ 1% đến 4% (khoảng 20 đến 100 mẫu) đảm bảo sai số vị trí xung luôn nhỏ nhất và ngăn ngừa biến động biên độ K-Norm bất thường.

Kỹ thuật trung bình miền bao hỗ trợ như thế nào cho quá trình chẩn đoán? Trung bình miền bao đóng vai trò làm mịn tín hiệu đầu ra của bộ giải mã, triệt tiêu các thành phần dao động cao tần ngẫu nhiên và làm nổi bật biên độ của các xung tuần hoàn. Điều này giúp các chỉ số thống kê như Kurtosis tăng từ 3 đến 4 lần, giúp việc nhận diện mức độ hư hỏng trở nên trực quan.

Thuật toán có xử lý được khi vận tốc quay của trục máy thay đổi liên tục không? Phương pháp hoàn toàn thích ứng tốt khi kết hợp với kỹ thuật theo dõi bậc (Order Tracking) để lấy mẫu lại tín hiệu theo góc quay trục. Khi tốc độ quay biến thiên nhỏ dưới 5%, thuật toán hiệu chỉnh vị trí véc-tơ lỗi cục bộ của Modified OMEDA vẫn đảm bảo trích xuất trọn vẹn 100% các xung khuyết tật.

Làm thế nào để xác định số lượng xung mục tiêu cần phân giải trong thực tế? Trong thực tế, số lượng xung mục tiêu được ước lượng dựa trên số vòng quay của trục khảo sát trong một bản ghi tín hiệu nhân với số lượng xung va đập cực đại dự kiến (thường chọn từ 3 đến 5 xung trên một chu kỳ quay), đảm bảo bao quát toàn bộ các khuyết tật tiềm ẩn trên bề mặt răng.

Kết luận

  • Đề xuất thành công phương pháp Modified OMEDA giải quyết triệt để hạn chế phụ thuộc vào chu kỳ lỗi của các kỹ thuật giải mã xung truyền thống MCKD và MOMEDA.
  • Xác lập nguyên lý ứng dụng chuẩn K-Norm thay thế D-Norm, kết hợp quy tắc khống chế chiều dài bộ lọc FIR từ 1% đến 4% giúp triệt tiêu 85% xung giả và định vị chính xác vị trí khuyết tật.
  • Chứng minh độ tin cậy vượt trội qua mô phỏng số và thực nghiệm với độ chính xác chẩn đoán đạt trên 95% trên 9 chế độ hư hỏng bánh răng khác nhau.
  • Khẳng định tính ưu việt về tốc độ tính toán với thời gian xử lý dưới 0,15 giây, tối ưu hóa toàn diện cho các ứng dụng giám sát rung động trực tuyến.
  • Đề ra lộ trình tích hợp thuật toán vào hệ thống SCADA/IoT công nghiệp trong vòng 6 đến 12 tháng tới nhằm hiện đại hóa công tác bảo trì dự đoán.

Đóng góp chính của luận văn là đã mở ra hướng tiếp cận mới trong kỹ thuật xử lý tín hiệu không dừng và cung cấp công cụ chẩn đoán mạnh mẽ cho ngành kỹ thuật cơ khí. Các doanh nghiệp sản xuất và đơn vị vận hành thiết bị công nghiệp nên sớm ứng dụng phương pháp Modified OMEDA vào hệ thống quản lý tài sản nhằm nâng cao độ tin cậy, ngăn ngừa sự cố thảm khốc và tối ưu hóa hiệu quả kinh tế dài hạn.