Tổng quan nghiên cứu
Trong các hệ thống máy quay tốc độ cao hiện đại như máy nén khí, turbine khí công nghiệp và máy bơm y tế, ma sát cơ học và hao mòn ổ trục truyền thống có thể gây thất thoát từ 15% đến 25% tổng năng lượng vận hành. Sự phát triển của công nghệ ổ từ chủ động (Active Magnetic Bearings - AMB) đã mở ra bước đột phá vượt bậc nhờ khả năng nâng trục quay không tiếp xúc bằng lực hút điện từ điều khiển vòng kín. Công nghệ này giúp triệt tiêu hoàn toàn ma sát cơ học, loại bỏ hệ thống dầu bôi trơn phức tạp và giảm tới 5% tổng khối lượng thiết bị.
Tuy nhiên, đối với hệ rotor 5 bậc tự do (5-DOF), cấu hình AMB hình trụ kinh điển đòi hỏi ít nhất hai cụm ổ hướng tâm và một cụm ổ hướng trục riêng biệt. Kết cấu này làm gia tăng đáng kể chiều dài trục, khối lượng và chi phí chế tạo. Ổ từ chủ động hình nón (Conical Active Magnetic Bearing - CAMB) ra đời như một giải pháp tối ưu hóa không gian, tích hợp đồng thời lực nâng hướng tâm và lực định vị dọc trục trên bề mặt côn vát góc nghiêng 0.98 rad. Mặc dù mang lại hiệu quả cơ khí vượt trội, CAMB lại sở hữu đặc tính phi tuyến mạnh, hiện tượng ghép kênh hình học phức tạp và độ giảm chấn tự nhiên rất nhỏ.
Luận văn thạc sĩ chuyên ngành Kỹ thuật Điều khiển và Tự động hóa của tác giả Tạ Thế Tài tại Đại học Bách Khoa Hà Nội (thực hiện năm 2023 dưới sự hướng dẫn của Tiến sĩ Nguyễn Danh Huy) tập trung mô hình hóa và phát triển giải thuật điều khiển phi tuyến tiên tiến cho hệ thống CAMB 5 bậc tự do. Đề tài xây dựng thành công bộ điều khiển trượt bậc phân số kết hợp bộ quan sát trạng thái mở rộng (FOSMC-ESO) nhằm triệt tiêu hiện tượng rung lắc chattering, rút ngắn thời gian xác lập dưới 0.02 giây và đảm bảo hệ thống vận hành ổn định trước các biến động tải trọng lên tới 100 N ở dải tốc độ quay 9000 vòng/phút.
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 nghiên cứu được xây dựng dựa trên sự kết hợp giữa định luật mạch từ Ampère, phương trình động lực học Newton-Euler cho vật thể quay 5 bậc tự do, lý thuyết điều khiển trượt bậc phân số và nguyên lý quan sát trạng thái mở rộng.
Mô hình tương tác điện từ xác định lực từ tỉ lệ thuận với bình phương dòng điện và tỉ lệ nghịch với bình phương khe hở không khí danh định 0.45 mm. Hệ thống sử dụng 4 cặp nam châm điện có diện tích mặt cực 118 mm2, mỗi cuộn dây gồm 82 vòng và được cấp dòng định thiên ban đầu 2 A.
Bốn khái niệm kỹ thuật trọng tâm bao gồm:
- Hiệu ứng ghép kênh con quay: Xuất hiện do tương tác giữa mô men quán tính đường kính 1.54e-4 kgm2 và mô men cực 31.68e-4 kgm2 khi rotor quay tốc độ cao.
- Nhiễu gộp hệ thống: Đại diện cho tổng các bất định tham số khối lượng rotor 0.755 kg, lực thủy động và rung động do lệch tâm khối lượng.
- Vi tích phân bậc phân số: Mở rộng bậc đạo hàm Riemann-Liouville với bậc phân số alpha được chọn bằng 0.1, tạo thêm bậc tự do linh hoạt cho mặt trượt.
- Hàm làm mượt Sigmoid: Thay thế hàm dấu Signum gián đoạn nhằm loại bỏ triệt để hiện tượng chattering trên tín hiệu điều khiển.
Phương pháp nghiên cứu
Nghiên cứu sử dụng tập dữ liệu mô phỏng gồm 10.000 điểm dữ liệu thời gian thực với chu kỳ trích mẫu 0.0001 giây trên nền tảng MATLAB và Simulink. Phương pháp chọn mẫu là chọn mẫu xác định theo kịch bản vận hành thực tế (Deterministic Scenario-based Sampling) với 4 kịch bản điển hình: khởi động từ độ lệch ban đầu z0 = 0.1 mm ở tốc độ 3000 vòng/phút, biến động bất định tham số ma trận 20%, chịu xung lực nhiễu ngoài 30 N đến 100 N ở tốc độ 9000 vòng/phút, và chịu lực kích động thủy động suy giảm tại thời điểm 0.04 giây.
Lý do lựa chọn phương pháp phân tích phi tuyến FOSMC-ESO: Hệ thống CAMB có tính ghép kênh dòng điện - hình học cao và đặc tính bất định lớn, các phương pháp tuyến tính hóa cổ điển như PID hoặc LQR không thể duy trì độ bền vững khi rotor lệch khỏi điểm làm việc. Thuật toán FOSMC-ESO cho phép ước lượng nhiễu gộp trực tuyến mà không cần cảm biến vận tốc đắt tiền, kết hợp mặt trượt phân số giúp tăng tốc độ hội tụ sai số. Tính ổn định tiệm cận toàn cục của hệ thống kín được chứng minh toán học nghiêm ngặt thông qua tiêu chuẩn hàm năng lượng Lyapunov. Toàn bộ quá trình nghiên cứu và mô phỏng được hoàn thành trong giai đoạn từ tháng 10 năm 2022 đến tháng 04 năm 2023 tại Trường Điện - Điện tử, Đại học Bách Khoa Hà Nội.
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à đánh giá đối chuẩn giữa FOSMC-ESO, SMC truyền thống và ADRC đã mang lại 4 phát hiện khoa học quan trọng:
Thứ nhất, về tốc độ đáp ứng và độ chính xác xác lập, FOSMC-ESO đưa rotor trở về vị trí cân bằng 5 bậc tự do chỉ trong 0.02 giây (20 ms) từ vị trí lệch ban đầu z0 = 0.1 mm, đạt độ quá điều chỉnh hoàn hảo 0%. Thành tích này nhanh hơn 28.5% so với bộ điều khiển ADRC (mất 28 ms và có độ quá điều chỉnh 3.85%).
Thứ hai, khi hệ thống chịu tác động của sai lệch tham số mô hình lên tới 20% đối với ma trận khối lượng và độ cứng lực từ, FOSMC-ESO duy trì sai số vị trí x, y tiệm cận 0 mm mà không có dao động quá độ, trong khi SMC thông thường dao động với biên độ 0.015 mm quanh bề mặt trượt.
Thứ ba, trước tác động của xung ngoại lực đột ngột Fd = 30 N và tăng mạnh lên 100 N ở tốc độ quay 9000 vòng/phút, bộ quan sát ESO đã ước lượng chính xác 100% biên độ nhiễu thực tế. Tín hiệu điều khiển bù trừ tức thời giúp khống chế độ lệch tâm trục tối đa dưới 0.035 mm (chiếm chưa đầy 7.7% khe hở không khí 0.45 mm), triệt tiêu hoàn toàn nguy cơ va chạm cơ học giữa rotor và stator.
Thứ tư, việc kết hợp toán tử bậc phân số alpha = 0.1 với hàm xấp xỉ Sigmoid đã loại bỏ 100% hiện tượng rung lắc chattering trên dòng điện điều khiển, giữ dao động dòng điện ổn định quanh mức 0.05 A tại trạng thái cân bằng.
Thảo luận kết quả
Hiệu năng vượt trội của FOSMC-ESO bắt nguồn từ cơ chế phân tách thông minh giữa bộ ước lượng nhiễu và bộ bù sai số động học. Bộ quan sát ESO ước lượng đồng thời vận tốc và toàn bộ nhiễu gộp mà không phụ thuộc vào độ chính xác tuyệt đối của phương trình vi phân rotor. Bề mặt trượt bậc phân số cung cấp bộ nhớ động học dài hạn, giúp quỹ đạo trạng thái hội tụ về mặt trượt mượt mà hơn cấu trúc số nguyên.
So với các công trình nghiên cứu trước đây sử dụng bộ điều khiển mờ Interval Type-2 Fuzzy hoặc điều khiển thích nghi đòi hỏi tài nguyên tính toán lớn, FOSMC-ESO có cấu trúc thuật toán gọn nhẹ nhưng vẫn duy trì tính bền vững cao trước các lực ghép kênh con quay và lực thủy động phức tạp.
Dữ liệu nghiên cứu được trình bày trực quan thông qua hệ thống 55 đồ thị thời gian thực và 3 bảng đối chuẩn định lượng. Đồ thị quá độ vị trí rotor mô tả rõ nét đường cong phản hồi của FOSMC-ESO với đặc tính trơn tru, không có gợn sóng dư. Bảng so sánh các chỉ số tích phân sai số tuyệt đối (IAE) và chỉ số năng lượng tiêu hao dòng điện khẳng định FOSMC-ESO tiết kiệm năng lượng hơn 15% so với SMC cổ điển.
Đề xuất và khuyến nghị
Nhằm hoàn thiện và đưa kết quả nghiên cứu vào ứng dụng thực tế công nghiệp, 4 nhóm giải pháp kỹ thuật cụ thể được đề xuất:
Thứ nhất, triển khai thực nghiệm giải thuật trên vi xử lý nhúng chuyên dụng. Chủ thể thực hiện: Nhóm kỹ sư nghiên cứu điều khiển tự động tại các viện nghiên cứu và trường đại học. Hành động: Chuyển đổi mã nguồn thuật toán FOSMC-ESO sang chip xử lý tín hiệu số DSP TMS320F28379D hoặc FPGA. Target metric: Đạt chu kỳ tính toán vòng lặp điều khiển dưới 50 microgiây với tần số lấy mẫu tối thiểu 20 kHz. Timeline: Hoàn thành trong 6 tháng đầu năm 2024.
Thứ hai, xây dựng cơ chế tự chỉnh định thông số bậc phân số trực tuyến. Chủ thể thực hiện: Kỹ sư phần mềm điều khiển công nghiệp. Hành động: Tích hợp thuật toán tối ưu hóa bầy đàn (PSO) hoặc mạng nơ-ron thích nghi để tự động cập nhật hệ số trượt c1, c2 và bậc vi phân alpha trong dải từ 0.05 đến 0.3 theo sự thay đổi của tải trọng. Target metric: Giảm sai số bám quỹ đạo động học xuống dưới 0.01 mm khi rotor thay đổi tốc độ từ 0 đến 12.000 vòng/phút. Timeline: Triển khai kiểm thử trong quý 3 và quý 4 năm 2024.
Thứ ba, thiết kế hệ thống vòng bi đỡ phụ an toàn và mạch bảo vệ quá dòng. Chủ thể thực hiện: Đơn vị chế tạo thiết bị quay chân không và máy nén khí. Hành động: Lắp đặt cơ cấu ổ bi đỡ phụ với khe hở an toàn 0.20 mm (chiếm khoảng 44% khe hở danh định) kết hợp cảm biến dòng xoáy eddy current 2 kênh. Target metric: Đảm bảo thời gian kích hoạt bảo vệ rotor dưới 5 mili-giây khi xảy ra sự cố sụt nguồn. Timeline: Thực hiện kiểm định an toàn định kỳ 12 tháng một lần từ năm 2025.
Thứ tư, chuẩn hóa quy trình gia công cơ khí chính xác và chế tạo lõi từ. Chủ thể thực hiện: Nhà máy chế tạo máy công cụ và vật liệu điện từ. Hành động: Kiểm soát dung sai góc nghiêng hình nón 0.98 rad và tăng hệ số điền đầy cuộn dây trên 75% bằng dây đồng phủ men chịu nhiệt cao cấp. Target metric: Giảm tổn hao dòng xoáy Foucault xuống dưới 3%, nâng tuổi thọ làm việc liên tục của stator lên trên 50.000 giờ. Timeline: Hoàn thiện quy trình trong vòng 18 tháng.
Đối tượng nên tham khảo luận văn
Luận văn là tài liệu tham khảo giá trị cao cho 4 nhóm đối tượng chính:
-
Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Kỹ thuật Điều khiển và Cơ điện tử: Tiếp cận phương pháp thiết lập phương trình động lực học phi tuyến 5 bậc tự do cho hệ CAMB và quy trình chứng minh tính ổn định Lyapunov cho hệ điều khiển phân số. Use case: Ứng dụng làm giáo trình giảng dạy chuyên đề điều khiển bền vững hoặc phát triển đề tài luận án tiến sĩ về hệ thống cơ điện tử phi tuyến.
-
Kỹ sư R&D tại các doanh nghiệp chế tạo máy công nghiệp nặng: Nắm bắt giải pháp tích hợp lực từ hướng tâm và dọc trục giúp rút ngắn 20% chiều dài trục máy và tiết kiệm không gian lắp ráp hộp số. Use case: Thiết kế hệ thống ổ từ cho máy nén ly tâm công nghiệp, máy bơm chân không và turbine phát điện công suất trên 500 kW.
-
Chuyên gia thiết kế thiết bị y tế công nghệ cao: Khai thác công nghệ nâng từ không tiếp xúc nhằm triệt tiêu hoàn toàn sự mài mòn vật liệu và sinh nhiệt cục bộ. Use case: Phát triển máy bơm máu nhân tạo hỗ trợ tâm thất (VAD) với kích thước nhỏ gọn dưới 50 mm, yêu cầu độ tin cậy tuyệt đối và thời gian vận hành liên tục không bảo dưỡng.
-
Kỹ sư quản lý vận hành trong ngành lọc hóa dầu và năng lượng hạt nhân: Hiểu rõ cơ chế bù trừ rung động không cân bằng và chẩn đoán trạng thái rotor từ xa thông qua tín hiệu ước lượng của bộ quan sát. Use case: Nâng cấp cụm ổ đỡ cơ khí sang ổ từ CAMB trong môi trường hóa chất ăn mòn và nhiệt độ cao trên 300 độ C, giúp tiết kiệm 40% chi phí bảo trì.
Câu hỏi thường gặp
Ổ từ chủ động hình nón CAMB có ưu điểm gì vượt trội so với ổ từ hình trụ truyền thống? CAMB kết hợp đồng thời lực nâng hướng tâm và lực định vị dọc trục trên bề mặt côn vát góc 0.98 rad. Thiết kế này loại bỏ hoàn toàn cụm 2 cuộn từ dọc trục độc lập, giúp tiết kiệm 20% chiều dài trục rotor, giảm 5% tổng trọng lượng thiết bị và cắt giảm 15% tổn hao nhiệt lượng trên cuộn dây.
Bộ điều khiển trượt bậc phân số FOSMC giải quyết hiện tượng chattering như thế nào? FOSMC kết hợp toán tử vi phân phân số bậc 0.1 với hàm trơn Sigmoid thay cho hàm dấu Signum gián đoạn. Cấu trúc này tạo ra mặt trượt mềm dẻo, triệt tiêu 100% hiện tượng chattering ở tần số cao, duy trì dao động dòng điện dưới 0.05 A và bảo vệ cuộn từ khỏi hiện tượng quá nhiệt.
Bộ quan sát trạng thái mở rộng ESO xử lý các nhiễu loạn của hệ thống ra sao? ESO coi toàn bộ bất định tham số 20%, lực con quay và tải trọng ngoài từ 30 N đến 100 N là một biến trạng thái mở rộng. Bộ quan sát ước lượng chính xác nhiễu trong vòng 0.02 giây, cho phép hệ thống bù trừ tức thời mà không cần đo đạc vận tốc hay phụ thuộc mô hình toán.
Hệ thống CAMB điều khiển bằng FOSMC-ESO có vận hành an toàn ở tốc độ 9000 vòng/phút không? Hệ thống hoàn toàn đảm bảo an toàn. Ở tốc độ 9000 vòng/phút, thuật toán khống chế lực ly tâm và mô men con quay quán tính cực 31.68e-4 kgm2, giữ độ lệch tâm rotor dưới 0.035 mm, nhỏ hơn nhiều so với khe hở không khí 0.45 mm, ngăn chặn triệt để nguy cơ cọ xát cơ học.
Việc cài đặt thuật toán FOSMC-ESO trên vi xử lý nhúng có khả thi về mặt tính toán thời gian thực không? Giải thuật hoàn toàn khả thi trên các chip xử lý tín hiệu số DSP có tốc độ từ 200 MHz trở lên. Phép tính vi tích phân phân số được xấp xỉ hiệu quả bằng bộ lọc Oustaloup bậc 5, giúp thời gian thực thi vòng lặp dưới 50 microgiây, đáp ứng tần số lấy mẫu điều khiển 20 kHz.
Kết luận
Công trình nghiên cứu luận văn thạc sĩ của tác giả Tạ Thế Tài đã giải quyết thành công bài toán điều khiển phi tuyến phức tạp cho hệ thống ổ từ chủ động hình nón 5 bậc tự do với 5 đóng góp cốt lõi:
- Xây dựng mô hình toán học giải tích 5 bậc tự do tích hợp hiệu ứng con quay, lực thủy động và khe hở không khí danh định 0.45 mm.
- Ứng dụng bộ quan sát trạng thái mở rộng ESO ước lượng chính xác 100% nhiễu gộp từ tải trọng 100 N và bất định tham số 20%.
- Phát triển bộ điều khiển FOSMC với bậc phân số 0.1 kết hợp hàm Sigmoid, triệt tiêu 100% hiện tượng chattering gây hại cơ khí.
- Chứng minh toán học chặt chẽ tính ổn định tiệm cận toàn cục của hệ thống kín bằng hàm năng lượng Lyapunov.
- Xác lập hiệu năng vận hành vượt trội qua mô phỏng với thời gian xác lập dưới 0.02 giây và độ quá điều chỉnh 0%.
Kế hoạch phát triển tiếp theo tập trung vào việc hiện thực hóa thuật toán trên phần cứng DSP nhúng trong vòng 12 tháng tới nhằm thử nghiệm ở dải tốc độ 12.000 vòng/phút.
Quý độc giả, nhà nghiên cứu và doanh nghiệp quan tâm có thể kết nối với tác giả cùng nhóm nghiên cứu tại Đại học Bách Khoa Hà Nội để cùng hợp tác chuyển giao và ứng dụng công nghệ ổ từ tiên tiến vào thực tiễn sản xuất công nghiệp.