Tổng quan về luận án

Nghiên cứu về hệ thống ổ từ chủ động (Active Magnetic Bearing - AMB) đại diện cho một bước chuyển dịch công nghệ quan trọng trong kỹ thuật cơ điện tử hiện đại, cho phép nâng và định vị rotor quay với tốc độ siêu cao mà hoàn toàn không có tiếp xúc cơ học, triệt tiêu ma sát, loại bỏ nhu cầu bôi trơn và nâng cao tuổi thọ vận hành trong các môi trường khắc nghiệt như chân không, nhiệt độ cực hạn, thiết bị bán dẫn và máy quang khắc chính xác cao. Tuy nhiên, rào cản kỹ thuật cốt lõi đặt ra đối với các cơ cấu truyền động ổ từ dọc trục (Axial Magnetic Bearing) nguyên khối là hiện tượng dòng điện xoáy (dòng Foucault) cảm ứng xuất hiện trong cấu trúc sắt từ của stator và rotor khi cuộn dây được cấp dòng điện biến thiên theo thời gian. Định luật cảm ứng điện từ Faraday khẳng định rằng dòng xoáy sinh ra từ trường ngược chiều chống lại sự biến thiên của từ trường kích thích, dẫn đến sự suy giảm nghiêm trọng của lực điện từ động, gây trễ pha lớn và làm biến dạng động lực học phi tuyến của hệ thống.

Khoảng trống nghiên cứu (research gap) trọng tâm được xác định rõ rệt: Các công trình nghiên cứu kinh điển về mô hình hóa dòng xoáy trong ổ từ (Zmood, 1987; Feeley, 1993; Kucera & Ahrens, 1995; Lei Zhu, Carl R. Knospe & Eric H. Maslen, 2005–2011) chủ yếu dừng lại ở việc mô hình hóa toán học cho cấu trúc một cực từ đơn lẻ dạng hình trụ hoặc hình chữ C trên miền tần số. Trong khi đó, các nghiên cứu trong nước (Nguyễn Thị Thanh Bình, 2013; Nguyễn Danh Huy, 2019; Vũ Đình Đạt, 2021) đều giả định stator và rotor ghép từ các lá thép kỹ thuật mỏng để triệt tiêu hoàn toàn tổn hao dòng xoáy. Việc thiếu vắng một mô hình toán học giải tích chính xác trên miền thời gian có dạng vi phân bậc phân số cho cấu trúc cực từ hình chữ E, đồng thời chưa có cấu trúc điều khiển phi tuyến hiện đại hoàn chỉnh cho hệ thống 3 cặp cực từ đặt lệch nhau $120^\circ$ có xét đến tương tác xen kênh (cross-coupling), tốc độ quay của rotor và bất định tham số, chính là nút thắt lý thuyết và thực tiễn mà luận án tiến sĩ kỹ thuật chuyên ngành Kỹ thuật điều khiển và Tự động hóa (Mã số: 9520216) của nghiên cứu sinh Lê Ngọc Hợi dưới sự hướng dẫn của PGS. Nguyễn Quang Địch và PGS. Nguyễn Tùng Lâm tại Đại học Bách Khoa Hà Nội (2024) tập trung giải quyết.

Luận án thiết lập hệ thống câu hỏi nghiên cứu và giả thuyết khoa học chặt chẽ:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để xây dựng mô hình toán học giải tích trên miền thời gian có xét ảnh hưởng của dòng xoáy cho cực từ dạng hình chữ E và hệ thống ổ từ dọc trục 3 cặp cực bố trí không gian đối xứng $120^\circ$?
  • Câu hỏi nghiên cứu 2 (RQ2): Quy luật trễ pha và suy giảm biên độ lực từ biến thiên theo tần số kích thích được lượng hóa chính xác như thế nào thông qua mô hình vi tích phân bậc phân số so với phân tích phần tử hữu hạn (FEM)?
  • Câu hỏi nghiên cứu 3 (RQ3): Cấu trúc điều khiển phi tuyến nào có khả năng bù trừ hoàn toàn ảnh hưởng của dòng xoáy, triệt tiêu hiện tượng rung giật (chattering) và ước lượng tổng nhiễu phi tuyến mà không cần sử dụng cảm biến đo tốc độ quay?

Các giả thuyết nghiên cứu bao gồm:

  • Giả thuyết 1 (H1): Cấu trúc cực từ hình chữ E mở rộng đường dẫn từ thông, làm gia tăng mật độ từ thông chính và lực từ định mức so với cực từ hình chữ C trên cùng kích thước bao hình và điều kiện kích từ.
  • Giả thuyết 2 (H2): Mô hình mạch từ tương đương tích hợp từ trở hiệu dụng bậc phân số ($R_i = R_{io} + k_i \sqrt{s}$) mô tả chính xác đặc tính suy giảm lực từ động học theo tần số kích từ $\omega \in [0, 200\text{ Hz}]$.
  • Giả thuyết 3 (H3): Bộ điều khiển trượt (Sliding Mode Control - SMC) tích hợp bộ quan sát trạng thái mở rộng bậc phân số (Fractional Order Extended State Observer - FO ESO) vượt trội hơn các bộ điều khiển PD và Backstepping (BSC) cổ điển về khả năng kháng nhiễu xen kênh, bù bất định mô hình và triệt tiêu cảm biến tốc độ.

Khung lý thuyết của luận án được xây dựng dựa trên sự giao thoa giữa Lý thuyết điện từ trường Maxwell - Faraday, Giải tích vi phân bậc phân số (Leibniz, Podlubny, Matignon) và Lý thuyết điều khiển phi tuyến kháng nhiễu chủ động (Sliding Mode Control, Extended State Observer). Đóng góp đột phá của công trình được định lượng trên mô hình thực nghiệm với rotor khối lượng $m = 2.25\text{ kg}$, khe hở không khí danh định $g_0 = 0.2\text{ mm}$ ($2 \times 10^{-4}\text{ m}$), dòng tiền từ hóa $i_0 = 0.2\text{ A}$, số vòng dây $N = 1200\text{ vòng}$, hệ số độ cứng dòng điện $K_i = 440\text{ N/A}$, hệ số độ cứng dịch chuyển $K_z = 3.4744 \times 10^5\text{ N/m}$, và hệ số ảnh hưởng dòng xoáy $k = 6.6344 \times 10^4\text{ A/Wb}$.

Literature Review và Positioning

Lịch sử nghiên cứu về ổ từ chủ động ghi nhận bước phát triển quan trọng từ phát minh ban đầu của Jesse Beams tại Đại học Virginia trong Thế chiến II cho đến Hội nghị Quốc tế về Công nghệ Treo từ tính (ISMB) lần thứ nhất tại Thụy Sĩ năm 1988. Tuy nhiên, việc mô hình hóa ảnh hưởng của dòng Foucault trong các cơ cấu chấp hành sắt từ nguyên khối luôn là một bài toán hóc búa, phân nhánh thành hai trường phái tiếp cận chính trong y văn quốc tế:

Trường phái giải tích mạch từ một chiều và hai chiều: Zmood (1987) khởi xướng việc thiết lập mô hình từ trường 1D cho cực từ hình chữ C nguyên khối có tỷ lệ giữa chiều rộng và chiều cao cực từ lớn. Dựa trên phân bố cường độ từ trường bề mặt $H = H_s \frac{\cosh(\eta y)}{\cosh(\eta b)}$ với $\eta = \sqrt{s \sigma \mu}$, Zmood đưa ra mô hình hàm truyền từ thông khe hở dạng rút gọn $\Phi(s) = \frac{\Phi(0)}{1 + sT_1}$ với $T_1 = \frac{4\sigma \mu b^2}{\pi^2}$. Tuy nhiên, mô hình này mất tính chính xác khi tỷ lệ hình học giữa chiều rộng và chiều cao cực xấp xỉ nhau. Ngay sau đó, Feeley (1993) mở rộng sang mô hình từ trường 2D cho thanh chữ nhật và đưa ra hàm xấp xỉ $\Phi(s) = \frac{\mu N I A}{l_{eff}(x) [1 + b \sqrt{2\sigma \mu_r \mu_0 s}]}$. Điểm hạn chế chung của Zmood (1987) và Feeley (1993) được Lei Zhu, Carl R. Knospe và Eric H. Maslen (2005–2011) chỉ ra là: cả hai đều giả định cấu hình mật độ từ thông trong mặt cắt ngang khe hở không khí giống hệt mặt cắt ngang phần tử sắt – một giả thiết chỉ đúng ở trạng thái tĩnh ($\omega = 0$) nhưng sai lệch nghiêm trọng khi từ trường biến thiên theo thời gian ($\omega > 0$). Kucera & Ahrens (1995) giải bài toán Maxwell cho tấm bán vô hạn đối với cực từ hình trụ, song tham số suy giảm $d$ trong mô hình lại phụ thuộc hoàn toàn vào thực nghiệm mà không gắn liền với kích thước hình học giải tích.

Trường phái điều khiển và bù trừ động lực học: Whitlow (2001) đã khảo sát các bộ điều khiển PID kết hợp hồi tiếp tuyến tính hóa (Feedback Linearization - FBL), trong khi Jianpeng Zhong và Lichuan Li (2015) đề xuất bộ điều khiển PID bậc phân số (FO PID) trên miền tần số cho cơ cấu chấp hành một cực từ. Mặc dù vậy, các nghiên cứu này chưa phát triển cấu trúc điều khiển trên không gian trạng thái miền thời gian cho hệ đa cực từ. Tại Việt Nam, các luận án tiến sĩ tiêu biểu như Nguyễn Thị Thanh Bình (2013), Nguyễn Danh Huy (2019) và Vũ Đình Đạt (2021) đều giải bài toán điều khiển trên giả định rotor và stator ghép lá thép kỹ thuật mỏng để bỏ qua dòng xoáy ($f_p = K_i i_p + K_z z$). Khi chuyển sang hệ rotor nguyên khối chịu tải trọng va đập và tốc độ quay cao, mô hình tuyến tính hóa này xuất hiện sai số mô hình cực lớn, dẫn đến suy giảm dự trữ ổn định.

Vị thế khoa học của luận án được xác lập rõ ràng tại giao điểm giữa mô hình hóa vi tích phân bậc phân số và điều khiển phi tuyến bền vững: Luận án không chỉ mở rộng mô hình giải tích cho cực từ hình E (với việc phân rã mạch từ thành 6 phần tử tương đương gồm 3 phần tử sắt và 3 phần tử không khí tích hợp toán tử vi phân bậc nửa $\sqrt{s}$), mà còn lần đầu tiên chuyển đổi toàn bộ hệ thống động lực học 3 cặp cực từ trục dọc $120^\circ$ sang hệ phương trình trạng thái bậc phân số trên miền thời gian, tích hợp bộ điều khiển trượt SMC cùng bộ quan sát trạng thái mở rộng FO ESO.

So sánh đối chuẩn quốc tế:

  • So với mô hình cực từ hình C của Lei Zhu & Carl R. Knospe (2005–2011), mô hình cực từ hình E của luận án gia tăng tiết diện dẫn từ, tăng từ thông cực đại và nâng biên độ lực điện từ hữu dụng trong cùng một thể tích lắp đặt.
  • So với các nghiên cứu tối ưu hóa thông số vật liệu của Henry Spece, Roger Fittro và Carl Knospe (2018), luận án tiến xa hơn bằng việc chế tạo thực nghiệm hoàn chỉnh phần cứng trên card dSPACE DS1104, chứng minh khả năng ổn định động học mà không cần sử dụng cảm biến đo vận tốc quay.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án tạo ra bước phát triển sâu sắc đối với Lý thuyết điều khiển hệ thống bậc phân số (Fractional-Order System Theory) được khởi xướng bởi Leibniz (1695), phát triển bởi Podlubny (1998, 1999), Matignon (1996, 1998), Mehmet Önder Efe (2012) và Takahiro Takamatsu, Hiromitsu Ohmori (2016). Bằng chứng thực nghiệm và mô phỏng số đã chứng minh tính chất vật lý của dòng Foucault trong môi trường sắt từ liên tục vốn dĩ mang bản chất vi tích phân bậc phân số (fractional calculus) với toán tử Laplace $s^{1/2}$ (tương ứng đạo hàm bậc $q = 0.5$).

Khung phân tích động lực học của hệ thống ổ từ dọc trục được thiết lập dựa trên hệ tọa độ không gian trạng thái mở rộng. Đối với mỗi nhánh cực từ $u, v, w$, phương trình trạng thái phi tuyến bậc phân số được mô tả tổng quát dưới dạng:

$$\begin{cases} D^{0.5} z_1 = z_2 \ D^{0.5} z_2 = z_3 \ D^{0.5} z_3 = z_4 \ D^{0.5} z_4 = z_5 \ D^{0.5} z_5 = f(z) + g(z) u(t) + d(t) \end{cases}$$

Trong đó, $z_1$ đại diện cho độ dịch chuyển khe hở $z$, $z_3$ là vận tốc của rotor, $z_5$ là gia tốc của rotor, $z_2$ và $z_4$ là các biến trạng thái trung gian phát sinh từ toán tử phân số của dòng xoáy, $u(t)$ là tín hiệu điều khiển dòng điện, và $d(t)$ là tổng nhiễu phi tuyến tập trung. Bằng cách thiết lập hàm Lyapunov bậc phân số $V(e) = \frac{1}{2} e^T e > 0$ và chứng minh đạo hàm phân số cấp $\alpha$ thỏa mãn $D^\alpha V(e) \le -\kappa e^T e < 0$, luận án cung cấp bằng chứng giải tích chặt chẽ về tính ổn định tiệm cận toàn cục theo chuẩn Mittag-Leffler cho hệ ổ từ dọc trục 3 cặp cực.

+-----------------------------------------------------------------------------+
|                      KHUNG PHÂN TÍCH TỔNG THỂ HỆ THỐNG                     |
+-----------------------------------------------------------------------------+
|                                                                             |
|   +--------------------------+          +-------------------------------+   |
|   |   Lý thuyết Điện từ      |          |  Giải tích Bậc Phân số        |   |
|   |   (Maxwell - Faraday)    |          |  (Leibniz - Podlubny)         |   |
|   +-------------+------------+          +---------------+---------------+   |
|                 |                                       |                   |
|                 +-------------------+-------------------+                   |
|                                     |                                       |
|                                     v                                       |
|                  +--------------------------------------+                   |
|                  | Mô hình Mạch từ Cực từ Hình E 6 Phần |                   |
|                  | tử có xét Tổn hao Dòng xoáy (s^0.5)  |                   |
|                  +------------------+-------------------+                   |
|                                     |                                       |
|                                     v                                       |
|                  +--------------------------------------+                   |
|                  | Hệ Phương trình Trạng thái Bậc Phân  |                   |
|                  | số 3 Cặp cực lệch nhau 120 độ        |                   |
|                  +------------------+-------------------+                   |
|                                     |                                       |
|                 +-------------------+-------------------+                   |
|                 |                                       |                   |
|                 v                                       v                   |
|   +---------------------------+           +-----------------------------+   |
|   | Điều khiển Trượt (SMC)    |           | Bộ quan sát Trạng thái Mở   |   |
|   | Kháng nhiễu cùng kênh     |           | rộng Phân số (FO ESO)       |   |
|   +-------------+-------------+           +--------------+--------------+   |
|                 |                                        |                  |
|                 +-------------------+--------------------+                  |
|                                     |                                       |
|                                     v                                       |
|                  +--------------------------------------+                   |
|                  | Điều khiển Không cảm biến Tốc độ     |                   |
|                  | (Sensorless Control), Triệt Chattering|                  |
|                  +--------------------------------------+                   |
+-----------------------------------------------------------------------------+

Khung phân tích độc đáo

Sự đột phá trong phương pháp phân tích của luận án nằm ở việc tích hợp đồng thời ba cấu trúc lý thuyết:

  1. Lý thuyết mạch từ phi tuyến phân tán: Tách cấu trúc cực từ hình E thành 6 phần tử từ trở phụ thuộc tần số: 3 phần tử sắt ($R_1, R_2, R_3$) và 3 phần tử khe hở không khí ($R_4, R_5, R_6$). Trong đó, nhánh giữa có tiết diện $A = 150 \times 10^{-6}\text{ m}^2$ dẫn gấp đôi từ thông so với hai nhánh biên.
  2. Lý thuyết quan sát trạng thái mở rộng phân số (FO ESO): Mở rộng công trình của Zhe Gao (2015) sang hệ cơ điện từ 3 pha không gian, biến tổng các thành phần bất định bao gồm: nhiễu xen kênh giữa các cặp cực từ $u, v, w$, nhiễu do lực hồi chuyển khi rotor quay quanh trục $x, y$ với mômen quán tính $J_x, J_y$, bất định độ dẫn điện $\sigma$ và độ từ thẩm $\mu_r$, thành một biến trạng thái mở rộng $z_6 = d(z, t)$ để ước lượng trực tuyến và triệt tiêu trong luật điều khiển trượt.
  3. Điều kiện biên xác định: Giả thiết từ thông tản và dòng rò trong khe hở không khí được giới hạn ở mức không đáng kể, vật liệu sắt từ làm việc dưới ngưỡng bão hòa từ ($B < 1.8\text{ T}$), và độ dịch chuyển rotor thỏa mãn điều kiện vi sai $|z| \ll g_0$.

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

Thiết kế nghiên cứu

Nghiên cứu vận hành theo thế giới quan thực chứng (positivism paradigm) với quy trình định lượng đa tầng, kết hợp chặt chẽ giữa suy luận giải tích toán học, mô phỏng trường vật lý 3D bằng phần tử hữu hạn (FEM), mô phỏng hệ thống điều khiển số trên MATLAB/Simulink (sử dụng FOMCON Toolbox) và thực nghiệm phần cứng Hardware-in-the-Loop trên nền tảng vi xử lý DSP dSPACE DS1104.

Cấu trúc ổ từ dọc trục thực nghiệm bao gồm một đĩa rotor quay bằng thép nguyên khối đặt giữa hai cụm stator, mỗi cụm gồm 3 cặp cực hình E bố trí lệch nhau $120^\circ$. Kích thước hình học chính xác của cực từ hình E được thiết lập: chiều cao cực từ $2a = 15\text{ mm}$ ($a = 7.5\text{ mm}$), chiều rộng cực từ $2b = 10\text{ mm}$ ($b = 5\text{ mm}$), chiều dài stator $c = 20\text{ mm}$, chiều dài rotor $h = 30\text{ mm}$, độ dẫn điện của sắt $\sigma = 2.5 \times 10^6\text{ S/m}$, độ từ thẩm tương đối $\mu_r = 1000$, độ từ thẩm chân không $\mu_0 = 4\pi \times 10^{-7}\text{ H/m}$, độ từ thẩm tuyệt đối $\mu = \mu_r \mu_0 = 1.256 \times 10^{-3}\text{ H/m}$.

Quy trình nghiên cứu rigorous

Quy trình giải tích mạch từ cực từ hình chữ E được thực hiện nghiêm ngặt qua 4 bước:

  • Bước 1: Xác định từ trở hiệu dụng các phần tử sắt:

    • Nhánh giữa: $R_1(s) = R_{1o} + k_{i1}\sqrt{s}$, trong đó $R_{1o} = \frac{c}{\mu_r \mu_0 A}$, $k_{i1} = \frac{c}{4(a+b)\mu_r \mu_0}\sqrt{\sigma}$.
    • Hai nhánh bên: $R_2(s) = R_3(s) = R_{2o} + k_{i2}\sqrt{s}$, trong đó $R_{2o} = \frac{4h + 2c}{\mu_r \mu_0 A}$, $k_{i2} = \left(\frac{2h+c}{4a+2b}\right)\sqrt{\frac{\sigma}{\mu_r \mu_0}}$.

  • Bước 2: Xác định từ trở hiệu dụng các phần tử khe hở không khí:

    • Hai nhánh bên: $R_4(s) = R_5(s) = R_{4o} + k_{g4}\sqrt{s}$, trong đó $R_{4o} = \frac{2g_0}{\mu_0 A}$, $k_{g4} = \left[\frac{b}{12a} - \frac{8b^2}{\pi^5 a^2}\tanh\left(\frac{\pi a}{b}\right)\right]\sqrt{\frac{\sigma}{\mu_r \mu_0}}$.
    • Nhánh giữa: $R_6(s) = R_{6o} + k_{g6}\sqrt{s}$, trong đó $R_{6o} = \frac{g_0}{\mu_0 A}$, $k_{g6} = \left[\frac{b}{6a} - \frac{32b^2}{\pi^5 a^2}\tanh\left(\frac{\pi a}{2b}\right)\right]\sqrt{\frac{\sigma}{\mu_r \mu_0}}$.

  • Bước 3: Thiết lập phương trình cân bằng sức từ động (Ampere): Tổng sức từ động tác dụng trong mạch kín tuân theo: $E_1(s) + E_2(s) + E_4(s) + E_6(s) = N I(s)$. Sau khi phân tách thành phần tĩnh $E_{io}$ và thành phần điều khiển biến thiên $E_{ip}(s)$, từ thông điều khiển nhánh giữa $\Phi_p(s)$ được rút ra chính xác:

$$\Phi_p(s) = \frac{N I_p(s) - \frac{N i_0}{g_0} z(s)}{R_0 + k\sqrt{s}}$$

Trong đó, tổng từ trở tĩnh $R_0 = R_{1o} + \frac{1}{2}R_{2o} + \frac{1}{2}R_{4o} + R_{6o} = 2.6526 \times 10^6\text{ A/Wb}$, và hệ số dòng xoáy tổng $k = k_{i1} + \frac{1}{2}k_{i2} + \frac{1}{2}k_{g4} + k_{g6} = 6.6344 \times 10^4\text{ A/Wb}$.

  • Bước 4: Thiết lập phương trình lực điện từ phi tuyến: Lực điện từ tuyến tính hóa quanh điểm làm việc cân bằng có dạng:

$$F_p(s) = \frac{K_i I_p(s) + K_z Z(s)}{1 + \frac{k}{R_0}\sqrt{s}}$$

Data và phân tích

Dữ liệu mô phỏng trường điện từ 3D trên ANSYS Maxwell được thiết lập với dải tần số kích từ biến thiên từ $0\text{ Hz}$ đến $200\text{ Hz}$. Khi đặt dòng điện điều hòa $\Delta I = 0.01\omega\text{ (A)}$ và dòng tĩnh $I_1 = 0.2\omega\text{ (A)}$, kết quả FEM khẳng định: thành phần lực tĩnh không phụ thuộc vào tần số, nhưng biên độ lực điều hòa biến thiên giảm mạnh khi tần số tăng lên. Sai lệch giữa công thức truyền thống (bỏ qua dòng xoáy $F_p(s) = K_i I_p(s)$) và kết quả mô phỏng FEM tăng dần theo tần số $\omega$, chứng minh mô hình bậc phân số của luận án là hoàn toàn khớp với bản chất vật lý thực tế.

Trong giai đoạn thiết kế điều khiển, 4 thuật toán được lập trình và đối chuẩn chi tiết:

  1. Bộ điều khiển PD: Dựa trên cấu trúc tuyến tính cổ điển, điều chỉnh tham số theo mô hình cực.
  2. Bộ điều khiển Backstepping (BSC): Thiết kế cho đối tượng đạo hàm bậc phân số với các biến điều khiển ảo $\alpha_1, \dots, \alpha_5$.
  3. Bộ điều khiển trượt (SMC): Xây dựng mặt trượt $\sigma = c_1 e_1 + c_2 e_2 + c_3 e_3 + c_4 e_4 + e_5$ với luật điều khiển bù trượt phi tuyến $u_{SMC} = -K_s \text{sign}(\sigma) - w \sigma$.
  4. Bộ điều khiển SMC kết hợp FO ESO: Tích hợp bộ quan sát bậc phân số 6 trạng thái với hệ số khuếch đại $\beta_i = C_6^i \omega_o^i$, ước lượng đồng thời biến trạng thái vận tốc rotor $z_3$ và hàm tổng nhiễu $z_6 = d(z, t)$.

Phát hiện đột phá và implications

Những phát hiện then chốt

  1. Quy luật suy giảm lực từ động học theo tần số: Phân tích FEM và mô hình giải tích chỉ ra rằng khi tần số kích thích tăng từ $0\text{ Hz}$ đến $200\text{ Hz}$, biên độ lực từ biến thiên của cực từ nguyên khối bị suy giảm hơn $35%$ so với giá trị tính toán theo lý thuyết truyền thống bỏ qua dòng xoáy.
  2. Ưu thế cấu trúc của cực từ hình chữ E: Với cùng thể tích đóng gói, cùng kích thước dây quấn ($N = 1200\text{ vòng}$) và cùng dòng điện $i_0 = 0.2\text{ A}$, cực từ hình E tăng mật độ từ thông trong khe hở không khí lên $18.4%$ và tăng lực từ cực đại lên $22.6%$ so với cực từ hình C nhờ việc chia đôi đường dẫn từ thông về hai nhánh biên đối xứng.
  3. Khả năng triệt tiêu cảm biến đo tốc độ của FO ESO: Bộ quan sát trạng thái mở rộng bậc phân số FO ESO ước lượng chính xác vận tốc quay và các trạng thái nội tại không đo được với sai số ước lượng hội tụ về lân cận 0 trong thời gian $t < 0.05\text{ s}$, cho phép thực hiện điều khiển không cảm biến (sensorless control) hoàn toàn tin cậy.
  4. Triệt tiêu hiện tượng rung giật (Chattering) và kháng nhiễu xen kênh: Khi kết hợp FO ESO với SMC, tín hiệu bù nhiễu $z_6$ được đưa trực tiếp vào luật điều khiển, cho phép giảm hệ số trượt $K_s$ từ $150$ xuống $25$, loại bỏ hơn $90%$ biên độ chattering trên tín hiệu dòng điện điều khiển, đồng thời duy trì độ lệch rotor trong khoảng ổn định $\Delta z \le \pm 0.005\text{ mm}$ ngay cả khi chịu nhiễu tải va đập $f_z = 50\text{ N}$ và nhiễu xen kênh giữa các trục $u, v, w$.

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp cơ sở toán học tường minh cho việc mô hình hóa và điều khiển các hệ cơ điện từ phi tuyến có tính chất tiêu tán năng lượng phức tạp thông qua giải tích vi phân bậc phân số, mở rộng biên giới của lý thuyết điều khiển tự động.
  • Về mặt phương pháp luận: Thiết lập quy trình chuẩn tắc kết hợp 4 cấp độ: Giải tích mạch từ $\to$ FEM 3D Maxwell $\to$ Mô phỏng trạng thái bậc phân số $\to$ Thực nghiệm DSP dSPACE DS1104, có thể tái sử dụng cho các cơ cấu chấp hành điện từ cao tốc khác.
  • Về mặt ứng dụng công nghiệp: Giải pháp loại bỏ cảm biến tốc độ giúp giảm chi phí chế tạo phần cứng, thu gọn kích thước cụm ổ đỡ, triệt tiêu rung động cơ học do lắp đặt cảm biến trên trục quay, mở đường cho việc thương mại hóa các hệ thống máy nén ly tâm siêu tốc, bánh đà tích năng lượng (Flywheel Energy Storage Systems) và máy phay cao tốc trong công nghiệp gia công chính xác.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn chỉ ra các giới hạn nghiên cứu:

  • Giả định bỏ qua hiện tượng trễ từ (hysteresis) và bão hòa từ sâu: Mô hình toán học hiện tại chỉ đúng trong miền từ hóa tuyến tính ($B < 1.6\text{ T}$). Khi dòng điện kích từ quá lớn hoặc rotor vận hành ở tốc độ siêu cao ($> 60.000\text{ vòng/phút}$), tổn hao trễ từ và hiệu ứng phát nhiệt sẽ làm thay đổi độ dẫn điện $\sigma$ và độ từ thẩm $\mu_r$.
  • Bỏ qua từ thông tản phức tạp tại các góc cực từ: Mạch từ tương đương 6 phần tử sử dụng các hệ số hiệu chỉnh hình học giải tích 2D, chưa phản ánh trọn vẹn sự méo dạng từ trường 3D tại các cạnh biên cực từ.
  • Thử nghiệm thực nghiệm ở dải tốc độ giới hạn: Do điều kiện an toàn phòng thí nghiệm, kiểm chứng thực nghiệm trên card dSPACE DS1104 chủ yếu tập trung vào chế độ nâng tĩnh (levitation) và tốc độ quay danh định dưới tải trọng nhân tạo.

Chương trình nghiên cứu 5–10 năm tiếp theo bao gồm:

  1. Tích hợp mô hình toán học bậc phân số có xét đồng thời hiện tượng trễ từ phi tuyến (Jiles-Atherton hysteresis model) và bão hòa từ vào bộ điều khiển thích nghi.
  2. Thiết kế bộ điều khiển nơ-ron mờ bậc phân số tự chỉnh định thông số (Adaptive Neuro-Fuzzy Fractional-Order SMC) để tự động thích ứng với sự thay đổi nhiệt độ cuộn dây và stator.
  3. Mở rộng điều khiển tích hợp hoàn chỉnh 5 bậc tự do (5-DOF AMB) kết hợp cả ổ từ hướng kính (radial AMB) và ổ từ dọc trục (axial AMB) trên cùng một nền tảng vi xử lý FPGA tốc độ cao.
  4. Triển khai thử nghiệm trên các hệ thống turbine khí công nghiệp và thiết bị lưu trữ năng lượng bánh đà hoạt động trong môi trường chân không sâu.

Tác động và ảnh hưởng

  • Tác động học thuật: Luận án mở ra một hướng nghiên cứu tiên phong trong cộng đồng khoa học điều khiển tại Việt Nam và khu vực, kết nối lĩnh vực toán học thuần túy (vi tích phân phân số) với kỹ thuật cơ điện tử công nghiệp. Các công trình công bố từ luận án trên các tạp chí thuộc danh mục ISI/Scopus và hội nghị quốc tế uy tín khẳng định năng lực hội nhập khoa học đỉnh cao.
  • Chuyển đổi công nghiệp: Cung cấp tài liệu thiết kế và giải pháp công nghệ làm chủ quy trình tính toán, chế tạo ổ từ chủ động không tiếp xúc tại Việt Nam. Các doanh nghiệp sản xuất thiết bị quay công nghiệp, máy nén khí không dầu và máy bơm chân không bán dẫn có thể ứng dụng trực tiếp thuật toán và cấu hình cực từ hình E để nội địa hóa sản phẩm.
  • Lợi ích xã hội và môi trường: Công nghệ ổ từ không tiếp xúc giúp triệt tiêu hoàn toàn nhu cầu sử dụng dầu bôi trơn độc hại, giảm tổn hao cơ học, tiết kiệm năng lượng điện vận hành lên đến $15–20%$ trong các hệ thống truyền động công nghiệp, đóng góp trực tiếp vào mục tiêu trung hòa carbon và phát triển công nghiệp xanh.

Đối tượng hưởng lợi

  • Nghiên cứu sinh và học viên cao học: Tiếp cận phương pháp luận mô hình hóa trạng thái bậc phân số hoàn chỉnh và kỹ thuật thiết kế bộ quan sát FO ESO trên MATLAB/Simulink.
  • Các nhà khoa học và chuyên gia điều khiển: Sở hữu khung lý thuyết đối chuẩn tin cậy về điều khiển trượt kháng nhiễu đa kênh cho các đối tượng cơ điện tử phi tuyến phức tạp.
  • Kỹ sư R&D trong ngành cơ điện tử và năng lượng: Nắm vững bản vẽ, thông số thiết kế mạch từ cực từ hình E và cấu trúc mạch biến đổi công suất điều khiển ổ từ thực tế.
  • Các cơ quan hoạch định chính sách công nghệ: Có cơ sở thực tiễn vững chắc để tài trợ cho các đề tài nghiên cứu làm chủ công nghệ lõi trong chiến lược phát triển công nghiệp phụ trợ công nghệ cao và thiết bị cơ khí chính xác thế hệ mới.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc thiết lập thành công mô hình toán học giải tích không gian trạng thái trên miền thời gian có dạng vi phân bậc phân số cấp $0.5$ ($\sqrt{s}$) cho cơ cấu chấp hành ổ từ dọc trục 3 cặp cực hình E lệch nhau $120^\circ$. Công trình này mở rộng trực tiếp Lý thuyết vi tích phân bậc phân số của Podlubny (1999) và Lý thuyết quan sát trạng thái mở rộng của Han (2009) & Zhe Gao (2015), chuyển hóa bài toán dòng xoáy vật lý thành một bài toán điều khiển bền vững kháng tổng nhiễu phi tuyến tập trung.

2. Đột phá về phương pháp luận so với ít nhất hai nghiên cứu quốc tế tiền nhiệm?

So với nghiên cứu của Lei Zhu & Carl R. Knospe (2005–2011) vốn chỉ dừng lại ở hàm truyền Laplace 1 cực từ hình chữ C trên miền tần số, luận án đã giải quyết trọn vẹn bài toán ghép nối 3 cặp cực hình E trên miền thời gian có xét tương tác xen kênh và động học quay rotor. So với công trình của Jianpeng Zhong & Lichuan Li (2015) chỉ áp dụng bộ điều khiển FO PID kinh điển, luận án đã đề xuất cấu trúc SMC kết hợp FO ESO, giải quyết triệt để vấn đề kháng nhiễu phi tuyến cùng kênh và ước lượng trạng thái không cần cảm biến tốc độ.

3. Phát hiện bất ngờ nhất trong quá trình phân tích dữ liệu là gì?

Phát hiện bất ngờ nhất là sự phi tuyến hóa cực mạnh của thành phần lực từ điều hòa theo tần số kích từ: Khi tần số dòng điện tăng, biên độ lực từ suy giảm phi tuyến theo hàm căn bậc hai của tần số ($\sqrt{\omega}$), khiến bộ điều khiển tuyến tính PD hoàn toàn mất khả năng duy trì vị trí cân bằng khi có tải trọng động tần số cao, trong khi bộ điều khiển SMC + FO ESO vẫn duy trì sai số bám vị trí tiệm cận 0 nhờ thành phần bù nhiễu chủ động.

4. Luận án có cung cấp đầy đủ giao thức tái lập nghiên cứu (replication protocol) không?

Có. Luận án công bố chi tiết toàn bộ thông số hình học và vật lý của cực từ hình E ($a = 7.5\text{ mm}$, $b = 5\text{ mm}$, $c = 20\text{ mm}$, $h = 30\text{ mm}$, $\sigma = 2.5 \times 10^6\text{ S/m}$, $\mu_r = 1000$, $N = 1200\text{ vòng}$), cấu trúc thuật toán điều khiển SMC và FO ESO, sơ đồ kết nối chân tín hiệu ADC/DAC của card dSPACE DS1104, cấu trúc bộ biến đổi công suất PWM và giao diện điều khiển thực nghiệm trên phần mềm ControlDesk.

5. Định hướng nghiên cứu 10 năm tới được phác thảo như thế nào?

Định hướng 10 năm tới tập trung vào: (1) Thương mại hóa hệ thống điều khiển ổ từ tích hợp trên chip nhúng FPGA tốc độ lấy mẫu cỡ micro-giây; (2) Ứng dụng giải thuật học máy thích nghi (Reinforcement Learning / Deep Neuro-Fuzzy) để tự động nhận dạng biến thiên tham số dòng xoáy theo nhiệt độ; (3) Triển khai thực tế trên hệ thống bánh đà lưu trữ năng lượng cao tốc trong lưới điện thông minh.

Kết luận

  1. Luận án đã giải quyết trọn vẹn bài toán mô hình hóa giải tích có xét ảnh hưởng của dòng điện xoáy cho cơ cấu chấp hành cực từ hình chữ E, chỉ ra rõ rệt ưu thế vượt trội của cực từ hình E so với hình C về khả năng mở rộng đường dẫn từ thông và gia tăng lực nâng điện từ.
  2. Thiết lập thành công hệ phương trình trạng thái phi tuyến bậc phân số trên miền thời gian cho hệ thống ổ từ dọc trục 3 cặp cực hình E bố trí lệch nhau $120^\circ$, tích hợp đầy đủ các yếu tố động học: dòng xoáy, xen kênh từ trường, tốc độ quay rotor và bất định tham số.
  3. Thiết kế thành công bộ điều khiển trượt (SMC) kết hợp bộ quan sát trạng thái mở rộng bậc phân số (FO ESO), chứng minh bằng giải tích Lyapunov về tính ổn định tiệm cận toàn cục của hệ kín.
  4. Triệt tiêu hoàn toàn nhu cầu sử dụng cảm biến đo tốc độ quay cơ học (sensorless control), ước lượng chính xác các biến trạng thái nội và triệt tiêu hơn $90%$ hiện tượng rung giật chattering trên cơ cấu chấp hành công suất.
  5. Kiểm chứng tính đúng đắn của mô hình lý thuyết và thuật toán điều khiển thông qua mô phỏng phần tử hữu hạn FEM trên ANSYS Maxwell, mô phỏng số FOMCON/Simulink và thực nghiệm phần cứng trên hệ thống dSPACE DS1104.
  6. Mở ra nhánh nghiên cứu chuyên sâu về ứng dụng vi tích phân bậc phân số trong kỹ thuật điều khiển các hệ truyền động cơ điện từ cao tốc, đóng góp nền tảng công nghệ quan trọng cho ngành chế tạo máy và tự động hóa công nghiệp chính xác cao.