Tổng quan về luận án

Nghiên cứu của Ngô Mạnh Tùng (2022) tại Trường Đại học Bách khoa Hà Nội, dưới sự hướng dẫn khoa học của TS. Phạm Quang Đăng và PGS. Nguyễn Huy Phương, tập trung vào đề tài: "Nghiên cứu hệ điều khiển động cơ đồng bộ kích thích vĩnh cửu từ trường dọc trục không dùng cảm biến" (Chuyên ngành: Kỹ thuật điều khiển và tự động hóa - Mã số: 9520216). Động cơ đồng bộ kích thích vĩnh cửu từ trường dọc trục (Axial-Flux Permanent Magnet Synchronous Motor - AFPMSM / KTVC TTDT) tích hợp ổ từ dọc trục đại diện cho bước đột phá trong kỹ thuật truyền động điện tốc độ cao, loại bỏ hoàn toàn ma sát cơ học, dầu bôi trơn và nguy cơ cháy nổ trong các môi trường làm việc khắc nghiệt như máy bơm phân tử chân không, máy nén khí hóa lỏng và bánh đà tích trữ năng lượng (UPS).

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực tế: phần lớn các nghiên cứu kinh điển về động cơ từ thông dọc trục (như nghiên cứu của Profumo và cộng sự, 2000; Parviainen, 2005) chỉ tiếp cận đối tượng dưới dạng hệ một bậc tự do (1-DOF - chuyển động quay), bỏ qua sự dịch chuyển dọc trục. Trong khi đó, các công trình khảo sát cấu trúc tích hợp ổ từ - động cơ (Chiba và cộng sự, 2005; Nguyen Huy Phuong, 2012) thường dựa trên mô hình toán học rút gọn, tuyến tính hóa quanh điểm làm việc tĩnh để áp dụng bộ điều khiển PID hoặc trượt cơ bản. Cách tiếp cận này tạo ra sai số tĩnh lớn, làm suy giảm tính bền vững của hệ thống trước biến thiên thông số và không triệt tiêu được hiện tượng tương tác xen kênh (cross-coupling) phi tuyến mạnh mẽ giữa lực nâng dọc trục $F$ và momen quay $T$. Hơn nữa, việc sử dụng cảm biến cơ học (encoder, resolver) ở dải siêu cao tốc làm tăng chi phí phần cứng, gia tăng kích thước trục, suy giảm độ cứng vững cơ khí và phát sinh rung lắc do quán tính.

Để giải quyết triệt để các hạn chế trên, nghiên cứu thiết lập 3 câu hỏi nghiên cứu và 3 giả thuyết tương ứng:

  • RQ1: Làm thế nào để xây dựng mô hình toán học phi tuyến toàn phần mô tả chính xác tương tác 2 bậc tự do (chuyển động quay $\theta_z$ và chuyển động tịnh tiến dọc trục $z$) có tính đến biến thiên khe hở không khí $g(t) = g_0 \pm z(t)$?
    • H1: Mô hình động học phi tuyến trên hệ tọa độ tựa theo từ thông rotor ($dq$) với luật phân bổ dòng điện vi sai sẽ giải trừ được ghép kênh giữa trục $d$ (lực dọc trục) và trục $q$ (momen quay).
  • RQ2: Cấu trúc điều khiển phi tuyến nào có khả năng khử hoàn toàn tác động của thành phần bất định và nhiễu tải mà không gặp hiện tượng bùng nổ toán hạng (explosion of complexity)?
    • H2: Sự kết hợp giữa kỹ thuật Backstepping cải tiến (bổ sung khâu tích phân sai số và bộ lọc bậc nhất từ lý thuyết Dynamic Surface Control) cùng bộ quan sát High-Gain (HGO) với khâu đổi biến trạng thái phụ sẽ đảm bảo tính ổn định tiệm cận toàn cục theo tiêu chuẩn Lyapunov.
  • RQ3: Phương pháp ước lượng không cảm biến nào đạt độ chính xác cao ở dải tốc độ trung bình và cao khi điện cảm stator biến thiên phi tuyến theo khe hở không khí?
    • H3: Bộ quan sát trượt (Sliding Mode Observer - SMO) và High-Gain Observer ước lượng sức điện động cảm ứng (Back-EMF) trên hệ tọa độ cố định $\alpha\beta$ kết hợp khâu hiệu chỉnh góc pha cải tiến sẽ tái cấu trúc chính xác vị trí góc $\theta_e$ và tốc độ $\omega_e$ với sai số quá độ tối thiểu.

Khung lý thuyết của công trình dựa trên lý thuyết điều khiển vector FOC (Field-Oriented Control), lý thuyết ổn định Lyapunov, lý thuyết điều khiển phi tuyến Backstepping thích nghi của Krstic - Kanellakopoulos - Kokotovic (1995), lý thuyết quan sát High-Gain của Hassan K. Khalil (2002) và điều khiển trượt của Vadim Utkin (1992). Nghiên cứu mang lại đóng góp đột phá khi kiểm chứng thành công trên hệ thống mô phỏng thời gian thực HIL (Hardware-in-the-Loop) với thiết bị Typhoon HIL402 và card DSP Texas Instruments TMS320F28379D ở tần số đóng cắt nghịch lưu PWM 20 kHz, xác nhận khả năng triệt tiêu tương tác xen kênh, triệt tiêu sai số vị trí dọc trục $\Delta z \approx 0$ và bám tốc độ chính xác từ $150\text{ rad/s}$ đến $200\text{ rad/s}$.

Literature Review và Positioning

Tổng quan tài liệu trong lĩnh vực truyền động ổ từ tích hợp và động cơ đồng bộ kích thích vĩnh cửu từ trường dọc trục cho thấy ba luồng nghiên cứu chính với các tranh luận học thuật sâu sắc:

                      ┌────────────────────────────────────────────────────────┐
                      │    HỆ TRUYỀN ĐỘNG ĐỘNG CƠ Ổ TỪ TÍCH HỢP TỐC ĐỘ CAO     │
                      └──────────────────────────┬─────────────────────────────┘
                                                 │
         ┌───────────────────────────────────────┼────────────────────────────────────────┐
         ▼                                       ▼                                        ▼
┌─────────────────────────────┐   ┌─────────────────────────────┐   ┌─────────────────────────────┐
│  Luồng 1: Cấu trúc tích hợp │   │ Luồng 2: Chiến lược điều    │   │ Luồng 3: Ước lượng Sensorless│
│   & Mô hình hóa 2-DOF       │   │ khiển phi tuyến & Bù nhiễu  │   │  dựa trên Sức điện động     │
├─────────────────────────────┤   ├─────────────────────────────┤   ├─────────────────────────────┤
│• Chiba et al. (2005)        │   │• PID / PI kinh điển         │   │• Sóng mang cao tần/INFORM:  │
│• Nguyen Huy Phuong (2012)   │   │  (Dixon et al., 2004)       │   │  Schrödler (2001)           │
│• Tranh luận: Mô hình rút gọn│   │• Sliding Mode / Backstepping│   │• Back-EMF Observer:         │
│  vs Mô hình đầy đủ chứa     │   │  (Swaroop et al., 2000;     │   │  Khalil (2002), Utkin (1992)│
│  thành phần $L(g)$ phi tuyến│   │   Krstic et al., 1995)      │   │• Tranh luận: HGO vs SMO     │
└──────────────┬──────────────┘   └──────────────┬──────────────┘   └──────────────┬──────────────┘
               │                                 │                                 │
               └─────────────────────────────────┼─────────────────────────────────┘
                                                 ▼
                      ┌────────────────────────────────────────────────────────┐
                      │        POSITIONING & ĐÓNG GÓP CỦA LUẬN ÁN              │
                      │ • Mô hình hóa phi tuyến đầy đủ 2-DOF cấu trúc Sandwich │
                      │ • Backstepping cải tiến + Bộ lọc DSC + HGO bù $d_1, d_2$│
                      │ • Bộ quan sát Back-EMF SMO/HGO kháng nhiễu đo qua biến $\xi$ │
                      └────────────────────────────────────────────────────────┘

Luồng nghiên cứu thứ nhất tập trung vào cấu trúc động cơ ổ từ và mô hình hóa máy điện. Các công trình kinh điển của Chiba và cộng sự (2005) cùng Schweitzer và Maslen (2009) đã phân tích chi tiết hệ động cơ ổ từ 6 bậc tự do truyền thống với 2 ổ từ ngang trục và 1 ổ từ dọc trục riêng biệt. Nhược điểm của cấu trúc này là kích thước trục quá dài, dẫn đến suy giảm tần số uốn tự nhiên và làm hệ thống cồng kềnh. Để khắc phục, cấu trúc tích hợp ổ từ dọc trục - động cơ dạng đĩa (Sandwich structure: 1 rotor kẹp giữa 2 stator) được đề xuất (Nguyen Huy Phuong, 2012; Deng et al., 2018). Tuy nhiên, một tranh luận học thuật lớn tồn tại: hầu hết các tác giả quốc tế (như Zhang et al., 2016; Amara et al., 2009) chấp nhận giả thiết đơn giản hóa, bỏ qua độ lệch trục $z$ khi tính toán từ thông và điện áp stator, dẫn đến việc coi các điện cảm $L_d, L_q$ là hằng số. Luận án chỉ ra rằng trong quá trình vận hành tốc độ cao, khe hở $g_1 = g_0 - z$ và $g_2 = g_0 + z$ thay đổi liên tục, làm điện cảm pha biến thiên theo quy luật: $$L_{sd} = \frac{3L'{sd0}}{2g} + L{sl}$$ Nếu bỏ qua thành phần phi tuyến này, mô hình sẽ mất độ chính xác nghiêm trọng khi có lực tải tác động.

Luồng nghiên cứu thứ hai liên quan đến thuật toán điều khiển vị trí dọc trục và tốc độ. Nhóm tác giả theo trường phái tuyến tính (Dixon et al., 2004; Zhao et al., 2014) sử dụng các bộ điều khiển PI/PID phân tán cho từng mạch vòng. Ngược lại, nhóm tác giả theo trường phái điều khiển phi tuyến hiện đại (Swaroop et al., 2000; Utkin et al., 2009) chỉ ra rằng PID không thể triệt tiêu sai số động học khi xuất hiện thành phần ghép kênh phi tuyến bậc hai của dòng điện ($i_d^2, i_q^2$). Mặc dù kỹ thuật điều khiển trượt (Sliding Mode Control - SMC) hoặc Backstepping kinh điển được ứng dụng (như trong nghiên cứu của Tarczewski & Grzesiak, 2014), nhưng SMC lại tạo ra hiện tượng dao động cưỡng bức (chattering) làm nóng cuộn dây và gây rung trục, còn Backstepping tiêu chuẩn lại gặp rào cản "bùng nổ toán hạng" do phải tính đạo hàm liên tiếp các tín hiệu điều khiển ảo.

Luồng nghiên cứu thứ ba tập trung vào kỹ thuật điều khiển không cảm biến (Sensorless Control). Ở dải tốc độ thấp và điểm dừng, phương pháp chèn tín hiệu điện áp tần số cao (High-Frequency Signal Injection) hoặc đo điện kháng trực tuyến INFORM (Indirect Flux Detection by Online Reactance Measurement - Schrödler, 2001; Vogelsberger et al., 2010) chiếm ưu thế nhưng đòi hỏi phần cứng phức tạp và phụ thuộc vào tính lồi cực từ. Ở dải tốc độ trung bình và cao, việc ước lượng dựa trên sức điện động cảm ứng (Back-EMF) thông qua bộ lọc Kalman mở rộng (EKF) hoặc bộ quan sát Luenberger (Luenberger, 1971; Bolognani et al., 2003) được sử dụng phổ biến. Tuy nhiên, EKF đòi hỏi khối lượng tính toán ma trận hiệp phương sai cực lớn, gây nghẽn vi xử lý DSP.

So sánh với hai nghiên cứu quốc tế tiêu biểu:

  1. Nghiên cứu của Zhang et al. (IEEE Trans. on Industrial Electronics, 2016) về điều khiển động cơ từ thông dọc trục không cảm biến sử dụng bộ quan sát trượt thông thường, nhưng chỉ áp dụng cho hệ 1-DOF có khe hở cố định, khiến hệ thống mất ổn định khi trục bị dịch chuyển dọc trục.
  2. Công trình của Chiba et al. (IEEE Trans. on Magnetics, 2005) về động cơ ổ từ 2-DOF sử dụng điều khiển tách kênh dựa trên mạng nơ-ron nhưng cấu trúc quá nặng nề, khó thực thi thời gian thực trên vi điều khiển nhúng thông dụng.

Luận án của Ngô Mạnh Tùng đã định vị chính xác điểm giao thoa: giải quyết bài toán điều khiển 2-DOF phi tuyến toàn phần trên cấu trúc Sandwich bằng cách tích hợp thuật toán Backstepping cải tiến với bộ lọc DSC, kết hợp bộ quan sát High-Gain kháng nhiễu và bộ quan sát trượt SMO ước lượng Back-EMF thời gian thực.

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

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

Luận án mở rộng và phát triển các lý thuyết điều khiển phi tuyến kinh điển trong bối cảnh máy điện đặc biệt:

  • Phát triển lý thuyết Backstepping thích nghi của Krstic et al. (1995): Luận án đã tích hợp thành công khâu tích phân sai số $k_{i1} \int e_1 dt$ trực tiếp vào từng bước lặp tổng hợp hàm Lyapunov, biến cấu trúc Backstepping truyền thống thành bộ điều khiển có khả năng triệt tiêu triệt để sai số xác lập (zero steady-state error) trước sự xuất hiện của các nhiễu hằng số và thành phần bất định.
  • Mở rộng lý thuyết Dynamic Surface Control của Swaroop et al. (2000): Đưa tín hiệu điều khiển ảo $x_1^, x_2^, x_3^$ qua các bộ lọc quán tính bậc nhất $\tau_i \dot{x}{id} + x{id} = x_i^$, loại bỏ hoàn toàn việc lấy đạo hàm giải tích các hàm phi tuyến phức tạp, giải quyết triệt để vấn đề "bùng nổ toán hạng" vốn là nhược điểm lớn nhất của Backstepping khi ứng dụng vào hệ thống nhúng.
  • Hiện thực hóa lý thuyết quan sát High-Gain của Hassan K. Khalil (2002): Chứng minh toán học rằng bộ quan sát High-Gain khi ứng dụng cho động cơ AFPMSM 2-DOF có sai số quan sát hội tụ theo hàm mũ về một lân cận $\mathcal{O}(\varepsilon)$ của gốc tọa độ: $$|d_i - \hat{d}_i| \le e^{-(1/\varepsilon_i)t} |d_i(0)| + \mathcal{O}(\varepsilon_i)$$
                                  KHUNG ĐIỀU KHIỂN TÍCH HỢP BÙ NHIỄU PHI TUYẾN
 ┌────────────────────────────────────────────────────────────────────────────────────────────────────────┐
 │                                                                                                        │
 │  ┌──────────────────────┐   Tín hiệu ảo $x_1^*$   ┌───────────────────┐  Tín hiệu đặt $x_{1d}$         │
 │  │ Bước 1: Vị trí trục z│────────────────────────►│ Bộ lọc bậc nhất   │──────────────────────┐         │
 │  │ $e_1 = x_{1d}-x_1    │                         │ $\tau_1\dot{x}_{1d}│                      │         │
 │  │      + k_{i1}\int e_1│                         │   + x_{1d} = x_1^*│                      │         │
 │  └──────────────────────┘                         └───────────────────┘                      │         │
 │             ▲                                                                                ▼         │
 │             │ Bù ước lượng $\hat{d}_1$                                             ┌──────────────────┐│
 │             │                                                                      │ Bước 2: Dòng $i_d││
 │  ┌──────────────────────┐                         ┌───────────────────┐            │ $e_2 = x_{2d}-x_2││
 │  │ Bộ quan sát High-Gain│                         │ Bộ lọc bậc nhất   │◄───────────│      +k_{i2}\inte││
 │  │ Ước lượng $d_1, d_2$ │                         │ $\tau_2\dot{x}_{2d}│ Tín hiệu ảo│ └────────┬───────┘│
 │  │ qua biến phụ $\xi$   │                         │   + x_{2d} = x_2^*│  $x_2^*$   │          │        │
 │  └──────────┬───────────┘                         └───────────────────┘            │          │        │
 │             │                                                                      │          ▼        │
 │             │ Bù ước lượng $\hat{d}_2$                                             │    Luật dòng $i_d$│
 │             ▼                                                                      │    (Lực nâng F)   │
 │  ┌──────────────────────┐   Tín hiệu ảo $x_3^*$   ┌───────────────────┐            │                   │
 │  │ Bước 3: Tốc độ $\omega_e$ │───────────────────►│ Bộ lọc bậc nhất   │────────────┴──────────┼────────┘
 │  │ $e_3 = x_{3d}-x_3    │                         │ $\tau_3\dot{x}_{3d}│   Tín hiệu đặt $x_{3d} │
 │  │      + k_{i3}\int e_3│                         │   + x_{3d} = x_3^*│                       ▼
 │  └──────────────────────┘                         └───────────────────┘               Luật dòng $i_q$
 │                                                                                       (Momen quay T)
 └────────────────────────────────────────────────────────────────────────────────────────────────────────┘

Mô hình toán học hoàn chỉnh thiết lập hệ 3 phương trình trạng thái phi tuyến với các biến $x_1 = z$ (độ lệch vị trí dọc trục), $x_2 = \dot{z}$ (vận tốc tịnh tiến), $x_3 = \omega_e$ (tốc độ góc rotor): $$\begin{cases} \dot{x}1 = x_2 \ \dot{x}2 = g_1 i_d + d_1 \ \dot{x}3 = g_2 i_q + d_2 \end{cases}$$ Trong đó các hằng số điều khiển được xác định cụ thể: $g_1 = \frac{4K{Fd}i_f}{m}$, $g_2 = \frac{2K_T}{J}$. Các thành phần bất định và nhiễu tải $d_1, d_2$ chứa đựng toàn bộ các tương tác ghép kênh và biến thiên thông số: $$d_1 = d_z - \frac{F_L}{m} = \frac{4K{Fd}i{d0}i_d + \frac{4}{g_0}\left[K_{Fd}(i_d^2 + i_{d0}^2 + i_f^2) + 2K_{Fd}i_f i_{d0} + K_{Fq}i_q^2\right]z - F_L}{m}$$ $$d_2 = d_\omega - \frac{T_L}{J} = \frac{2K_R i_q i_{d0} + \frac{2K_R i_d i_q z}{g_0} - T_L}{J}$$

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

Khung phân tích của luận án tích hợp đồng thời ba cấu trúc lý thuyết:

  1. Nguyên lý phân bổ dòng điện vi sai cho cấu trúc 2 Stator: Thiết lập luật điều khiển dòng điện stator 1 và stator 2 nhằm tách biệt hoàn toàn việc tạo lực và tạo momen: $$\begin{cases} i_{d1} = i_{d0} - i_d \ i_{d2} = i_{d0} + i_d \ i_{q1} = i_{q2} = i_q \end{cases}$$ Nguyên lý này đảm bảo: tổng lực dọc trục $F = F_1 - F_2$ tỷ lệ tuyến tính với dòng điều khiển $i_d$, trong khi tổng momen quay $T = T_1 + T_2$ tỷ lệ thuận với $i_q$, đưa hệ thống về dạng điều khiển decoupled quasi-linear.

  2. Kỹ thuật đổi biến trạng thái phụ chống khuếch đại nhiễu đo: Điểm nghẽn cố hữu của bộ quan sát High-Gain truyền thống là việc sử dụng hệ số khuếch đại lớn $\frac{1}{\varepsilon_i}$ sẽ khuếch đại trực tiếp nhiễu đo đạo hàm trạng thái $\dot{x}_2, \dot{x}_3$. Luận án đề xuất cấu trúc đổi biến phụ: $$\xi_1 = \hat{d}_1 - \frac{x_2}{\varepsilon_1}, \quad \xi_2 = \hat{d}_2 - \frac{x_3}{\varepsilon_2}$$ Dẫn đến phương trình ước lượng nhiễu ở miền Laplace có dạng bộ lọc thông thấp tự nhiên: $$\hat{d}_1 = \frac{-\frac{x_2}{\varepsilon_1} + g_1 i_d}{\varepsilon_1 s + 1}, \quad \hat{d}_2 = \frac{-\frac{x_3}{\varepsilon_2} + g_2 i_q}{\varepsilon_2 s + 1}$$ Cấu trúc này triệt tiêu hoàn toàn sự bùng nổ nhiễu đo tần số cao ở đầu ra bộ quan sát.

  3. Điều kiện biên xác định (Boundary Conditions):

    • Giới hạn chuyển động dọc trục: $|z| \ll g_0$ (thực tế $|z| < 0.1\text{ mm}$ so với $g_0 = 1.0\text{ mm}$).
    • Giới hạn dòng phân cực: $i_{d0} \approx 0\text{ A}$ nhằm cực tiểu hóa momen từ trở ký sinh $T_{12}$.
    • Dải tốc độ làm việc của bộ quan sát Back-EMF: $\omega_e \ge 0.1\omega_{rated}$ (đảm bảo biên độ sức điện động vượt ngưỡng tỷ số tín hiệu trên nhiễu SNR).

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ chặt chẽ thế giới quan thực chứng khoa học (Positivism) và chủ nghĩa hiện thực phản biện (Critical Realism), kết hợp giữa mô hình hóa giải tích toán học thuần túy, mô phỏng số học định lượng và kiểm chứng thực nghiệm phần cứng thời gian thực HIL.

                          QUY TRÌNH PHƯƠNG PHÁP NGHIÊN CỨU TOÀN DIỆN
┌──────────────────────────────────────────────────────────────────────────────────────────────────┐
│ GIAI ĐOẠN 1: MÔ HÌNH HÓA TOÁN HỌC GIẢI TÍCH                                                      │
│ • Thiết lập hệ phương trình vi phân phi tuyến 2-DOF trên hệ trục $dq$ và $\alpha\beta$          │
│ • Mô tả định lượng quan hệ điện cảm - khe hở $L_d(g), L_q(g)$ & Khai triển Maclaurin chuỗi Taylor│
└────────────────────────────────┬─────────────────────────────────────────────────────────────────┘
                                 │
                                 ▼
┌──────────────────────────────────────────────────────────────────────────────────────────────────┐
│ GIAI ĐOẠN 2: THIẾT KẾ ĐIỀU KHIỂN & BỘ QUAN SÁT                                                   │
│ • Tổng hợp bộ điều khiển vị trí/tốc độ: Modified Backstepping + DSC Filter                       │
│ • Thiết kế HGO ước lượng nhiễu $d_1, d_2$ & SMO/HGO ước lượng Back-EMF $e_{s\alpha}, e_{s\beta}$ │
│ • Chứng minh ổn định tiệm cận Lyapunov toàn phần: $\dot{V} \le -\sum \alpha_i e_i^2 + \delta$    │
└────────────────────────────────┬─────────────────────────────────────────────────────────────────┘
                                 │
                                 ▼
┌──────────────────────────────────────────────────────────────────────────────────────────────────┐
│ GIAI ĐOẠN 3: KIỂM CHỨNG TRÊN PHẦN CỨNG THỜI GIAN THỰC (HIL SIMULATION)                          │
│ • Mô phỏng Plant phi tuyến trên Typhoon HIL402 (Thời gian trích mẫu phần cứng: $1\ \mu\text{s}$) │
│ • Nạp thuật toán điều khiển C/C++ vào DSP TMS320F28379D (Tần số PWM nghịch lưu: 20 kHz)         │
│ • Đánh giá đáp ứng động học khi quá độ tốc độ ($150 \to 200\text{ rad/s}$) & Đột biến tải $F_L, T_L$│
└──────────────────────────────────────────────────────────────────────────────────────────────────┘

Thiết kế đa tầng (Multi-level design) bao gồm:

  • Cấp 1 (Cơ sở vật lý máy điện): Mô hình hóa cấu trúc Sandwich gồm 2 stator 3 pha và 1 đĩa rotor nam châm vĩnh cửu đất hiếm (NdFeB), số đôi cực $p_c = 2$.
  • Cấp 2 (Mạch vòng điều khiển bên trong - Inner loops): Hai mạch vòng điều khiển dòng điện stator ($i_{d1}, i_{q1}$ và $i_{d2}, i_{q2}$) hoạt động ở tần số lấy mẫu 20 kHz, tích hợp trong Inverter Board.
  • Cấp 3 (Mạch vòng điều khiển bên ngoài - Outer loops): Mạch vòng điều khiển vị trí dọc trục $z$ và mạch vòng điều khiển tốc độ góc $\omega_e$ được tổng hợp bằng thuật toán đề xuất.
  • Cấp 4 (Mạch vòng ước lượng không cảm biến - Sensorless Estimation loops): Bộ quan sát trạng thái tái tạo sức điện động cảm ứng $e_{s\alpha}, e_{s\beta}$ và vòng khóa pha PLL / khâu hiệu chỉnh góc pha tái tạo $\theta_e, \omega_e$.

Quy trình nghiên cứu rigorous

Quy trình thu thập dữ liệu và xử lý tín hiệu thực hiện qua các bước chuẩn hóa nghiêm ngặt:

  1. Thu thập tín hiệu dòng - áp: Cảm biến dòng điện Hall đo dòng 2 pha $i_u, i_v$, suy ra $i_w = -(i_u + i_v)$, chuyển đổi qua hệ tọa độ cố định $\alpha\beta$ bằng phép biến đổi Clarke.
  2. Ước lượng Sức điện động (Back-EMF Estimation):
    • Bộ quan sát High-Gain: Sử dụng mô hình động học dòng điện stator trên hệ $\alpha\beta$: $$\dot{\hat{i}}{s\alpha} = -\frac{R_s}{L_s}\hat{i}{s\alpha} + \frac{1}{L_s}u_{s\alpha} - \frac{1}{L_s}\hat{e}{s\alpha} + \frac{1}{\varepsilon\alpha}(i_{s\alpha} - \hat{i}{s\alpha})$$ $$\dot{\hat{e}}{s\alpha} = \frac{1}{\varepsilon_\alpha^2}(i_{s\alpha} - \hat{i}_{s\alpha})$$
    • Bộ quan sát trượt (SMO): Thay thế khâu hiệu chỉnh tuyến tính bằng hàm dấu (switching function): $$\dot{\hat{i}}{s\alpha} = -\frac{R_s}{L_s}\hat{i}{s\alpha} + \frac{1}{L_s}u_{s\alpha} - \frac{1}{L_s}k_t \text{sgn}(\hat{i}{s\alpha} - i{s\alpha})$$
  3. Triangulation (Tam giác đạc phương pháp): Đánh giá chéo kết quả giữa 3 phương pháp độc lập: (1) Mô phỏng số trên môi trường MATLAB/Simulink; (2) Mô phỏng thời gian thực HIL cấp tín hiệu với Typhoon HIL402; (3) Đối chuẩn với mô hình lý thuyết giải tích giải bằng hàm Lyapunov.

Data và phân tích

Thông số thực nghiệm và cấu hình phần cứng:

  • Động cơ KTVC TTDT cực lồi thử nghiệm: Công suất định mức $P_n = 1.5\text{ kW}$, điện trở stator $R_s = 1.2\ \Omega$, điện cảm trục $d$ định mức $L_{sd0} = 8.5\text{ mH}$, điện cảm trục $q$ định mức $L_{sq0} = 12.4\text{ mH}$, từ thông nam châm vĩnh cửu $\lambda_m = 0.175\text{ Wb}$, khe hở không khí danh định $g_0 = 1.0\text{ mm}$, khối lượng rotor $m = 2.5\text{ kg}$, momen quán tính $J = 0.0048\text{ kg}\cdot\text{m}^2$.
  • Hệ thống mô phỏng thời gian thực: Thiết bị Typhoon HIL402 với độ trễ bước tính (time-step) cực thấp $1\ \mu\text{s}$, kết nối trực tiếp qua cổng giao tiếp quang và card DSP Interface với chip xử lý số tín hiệu 32-bit kép Texas Instruments TMS320F28379D chạy ở xung nhịp 200 MHz.
  • Phần mềm điều khiển: Typhoon HIL Control Center v2022.1, Code Composer Studio (CCS) v10.x biên dịch mã nguồn C tối ưu hóa toán học số thực dấu phẩy động (floating-point unit).

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

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

              BẢNG SO SÁNH HIỆU NĂNG ĐIỀU KHIỂN & ƯỚC LƯỢNG TRONG LUẬN ÁN
┌──────────────────────────────────────────┬────────────────────────────┬────────────────────────────┐
│ Chỉ tiêu đánh giá                        │ Thuật toán kinh điển (PID) │ Đề xuất của Luận án        │
├──────────────────────────────────────────┼────────────────────────────┼────────────────────────────┤
│ Sai số vị trí dọc trục $\Delta z$ tĩnh   │ $0.045\text{ mm}$ (khi có  │ $\approx 0.000\text{ mm}$  │
│                                          │ tải $F_L = 50\text{ N}$)   │ (triệt tiêu hoàn toàn)     │
├──────────────────────────────────────────┼────────────────────────────┼────────────────────────────┤
│ Độ sụt tốc độ khi đóng tải $T_L = 5\text{ Nm}$│ $\Delta\omega = 18.5\text{ rad/s}$ │ $\Delta\omega \le 2.1\text{ rad/s}$ │
├──────────────────────────────────────────┼────────────────────────────┼────────────────────────────┤
│ Thời gian phục hồi tốc độ $t_s$          │ $120\text{ ms}$            │ $15\text{ ms}$             │
├──────────────────────────────────────────┼────────────────────────────┼────────────────────────────┤
│ Độ lệch pha ước lượng vị trí góc $\Delta\theta_e$│ $8.5^\circ$ (HGO cơ bản)   │ $\le 1.2^\circ$ (SMO cải tiến)│
├──────────────────────────────────────────┼────────────────────────────┼────────────────────────────┤
│ Hiện tượng bùng nổ toán hạng             │ Nghiêm trọng (Backstepping)│ Không xảy ra (nhờ DSC)     │
└──────────────────────────────────────────┴────────────────────────────┴────────────────────────────┘
  1. Triệt tiêu hoàn toàn sai số tĩnh vị trí dọc trục: Khi có lực tải dọc trục $F_L$ tác động đột ngột, bộ điều khiển Backstepping cải tiến kết hợp bộ quan sát High-Gain ước lượng chính xác giá trị $d_1$, bù trực tiếp vào kênh điều khiển dòng $i_d$. Kết quả thực nghiệm HIL xác nhận rotor nhanh chóng trở về vị trí cân bằng $z = 0\text{ mm}$ trong thời gian dưới $8\text{ ms}$, trong khi thuật toán không có bù nhiễu duy trì sai lệch tĩnh lên đến $0.045\text{ mm}$.
  2. Kháng nhiễu momen tải vượt trội: Khi áp đặt bước nhảy momen tải $T_L = 5\text{ Nm}$ tại thời điểm $t = 0.5\text{ s}$, độ sụt tốc độ của hệ thống đề xuất chỉ là $2.1\text{ rad/s}$ (tương đương $1.05%$ ở tốc độ $200\text{ rad/s}$), thời gian hồi phục $15\text{ ms}$, vượt trội hoàn toàn so với bộ điều khiển PID truyền thống (độ sụt $18.5\text{ rad/s}$, thời gian hồi phục $120\text{ ms}$).
  3. Phát hiện bất ngờ về so sánh HGO và SMO: Mặc dù bộ quan sát High-Gain có tốc độ hội tụ giải tích nhanh hơn về mặt lý thuyết, nhưng trong môi trường thực nghiệm thời gian thực có chứa sóng hài bậc cao từ biến tần PWM 20 kHz, bộ quan sát trượt SMO kết hợp khâu lọc thông thấp và hiệu chỉnh góc pha mang lại chất lượng ước lượng góc $\theta_e$ mịn hơn, sai số góc $\Delta\theta_e \le 1.2^\circ$, trong khi HGO có độ gợn sai số $\Delta\theta_e \approx 3.5^\circ$.
  4. Khả năng tách kênh hoàn hảo giữa trục $d$ và trục $q$: Quá trình thay đổi tốc độ đặt từ $150\text{ rad/s}$ lên $200\text{ rad/s}$ và ngược lại chỉ gây ra biến động nhất thời trên dòng điện $i_q$, trong khi vị trí dọc trục $z$ và dòng điện $i_d$ hoàn toàn không bị ảnh hưởng (biến thiên vị trí $\Delta z < 0.002\text{ mm}$). Điều này chứng minh hiệu quả tuyệt đối của luật phân bổ dòng điện vi sai và cấu trúc bù nhiễu $d_1, d_2$.

Implications đa chiều

  • Về mặt lý thuyết: Chứng minh tính khả thi của việc tổng hợp hệ thống điều khiển phi tuyến đa biến phức tạp bậc cao bằng cách kết hợp có chọn lọc giữa kỹ thuật Backstepping, Dynamic Surface Control và High-Gain Observer, mở ra hướng tiếp cận chuẩn tắc cho các hệ cơ điện tử phi tuyến 2-DOF khác.
  • Về mặt phương pháp luận: Thiết lập quy trình chuẩn kiểm chứng thuật toán điều khiển máy điện phức tạp trên nền tảng Typhoon HIL402 kết nối DSP thực, rút ngắn khoảng cách từ nghiên cứu mô phỏng thuần túy đến chế tạo sản phẩm công nghiệp.
  • Về mặt ứng dụng thực tiễn: Cho phép loại bỏ hoàn toàn encoder/resolver quang học đắt tiền (thường chiếm $25 - 40%$ giá thành phần cứng hệ thống truyền động tốc độ cao), giảm kích thước chiều dài trục động cơ ít nhất $20%$, nâng cao độ tin cậy và tuổi thọ vận hành trong các môi trường độc hại, nhiệt độ cực hạn.

Limitations và Future Research

Nghiên cứu thẳng thắn thừa nhận các giới hạn học thuật và biên giới ứng dụng:

  1. Hiện tượng mất quan sát ở dải tốc độ cực thấp và điểm dừng ($\omega_e \to 0$): Do sức điện động cảm ứng tỷ lệ thuận với tốc độ ($e_s \propto \omega_e \lambda_m$), khi động cơ dừng hoặc quay rất chậm, biên độ Back-EMF tiệm cận 0 và bị chìm hoàn toàn trong nhiễu đo dòng điện và điện áp rơi trên van bán dẫn. Hệ thống bắt buộc phải sử dụng phương pháp khởi động vòng hở $I-f$ (dòng điện không đổi, tăng dần tần số) trước khi chuyển đổi mượt sang chế độ vòng kín không cảm biến.
  2. Ảnh hưởng của biến thiên nhiệt độ cuộn dây: Điện trở stator $R_s$ thay đổi theo nhiệt độ vận hành thực tế (lên đến $+40%$ khi nhiệt độ cuộn dây tăng từ $25^\circ\text{C}$ lên $120^\circ\text{C}$). Luận án coi $R_s$ là hằng số danh định, phần sai lệch được dồn vào thành phần bất định $d_1, d_2$, điều này có thể làm tăng nhẹ gánh nặng tính toán của bộ quan sát nhiễu ở các chế độ quá tải dài hạn.
  3. Giới hạn bậc tự do cơ học: Nghiên cứu giả định các chuyển động ngang trục ($x, y, \theta_x, \theta_y$) đã được kiểm soát hoàn hảo bởi hai ổ từ ngang trục độc lập, chưa xét đến hiện tượng cộng hưởng cơ học chéo giữa trục ngang và trục dọc khi hệ thống vận hành vượt qua các tần số uốn nguy hiểm (critical bending speeds).

Chương trình nghiên cứu tiếp theo (Future Research Agenda) gồm 4 hướng cụ thể:

  • Tích hợp thuật toán tiêm tín hiệu điện áp tần số cao (HFI) ở dải tốc độ thấp với bộ quan sát Back-EMF ở dải tốc độ cao để tạo ra giải pháp Sensorless toàn dải tốc độ ($0 - 100%$ tốc độ định mức).
  • Bổ sung thuật toán nhận dạng trực tuyến điện trở stator $R_s$ và nhiệt độ cuộn dây dựa trên trí tuệ nhân tạo (AI/Neural Networks) thích nghi.
  • Mở rộng mô hình điều khiển tích hợp đồng thời 5 bậc tự do (5-DOF Active Magnetic Bearing System) để triệt tiêu hoàn toàn rung động không gian.
  • Hiện thực hóa hệ thống trên mô hình động cơ vật lý hoàn chỉnh công suất lớn ($> 10\text{ kW}$) trong các ứng dụng máy nén ly tâm công nghiệp.

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

  • Tác động học thuật (Academic Impact): Dự báo thu hút sự quan tâm và trích dẫn lớn từ cộng đồng nghiên cứu máy điện và điều khiển tự động quốc tế (IEEE Transactions on Industrial Electronics, IEEE Transactions on Power Electronics, Mechatronics). Đóng vai trò tài liệu tham khảo chuẩn mực cho các nghiên cứu sinh tiến sĩ về truyền động động cơ ổ từ tích hợp.
  • Chuyển đổi công nghiệp (Industry Transformation): Ứng dụng trực tiếp trong các ngành công nghiệp công nghệ cao:
    • Sản xuất chất bán dẫn & Quang điện tử: Bơm hút chân không phân tử siêu sạch (Turbomolecular pumps) không dầu, tốc độ quay $30.000 - 60.000\text{ vòng/phút}$.
    • Năng lượng tái tạo & Lưu trữ: Hệ thống bánh đà tích năng (Flywheel Energy Storage) với hiệu suất chu trình $> 90%$, tuổi thọ $> 20\text{ năm}$.
    • Dầu khí & Khí hóa lỏng (LNG): Máy nén khí dạng ống chìm vận hành an toàn tuyệt đối chống cháy nổ.
  • Giá trị kinh tế & Xã hội: Giảm chi phí sản xuất thiết bị truyền động tốc độ cao từ $30 - 45%$ nhờ tinh giản cấu trúc cơ khí và loại bỏ cảm biến đo vị trí đắt tiền; cắt giảm $100%$ chi phí bảo trì thay thế dầu bôi trơn và vòng bi cơ khí định kỳ.

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

  • Nghiên cứu sinh & Nhà khoa học trẻ (Doctoral Researchers): Tiếp cận phương pháp luận toán học mẫu mực về thiết kế điều khiển Backstepping cải tiến kháng nhiễu, kỹ thuật đổi biến High-Gain và quy trình thực nghiệm HIL chuẩn quốc tế.
  • Giảng viên & Chuyên gia học thuật cao cấp (Senior Academics): Nguồn tài liệu chuyên sâu để giảng dạy các học phần Điều khiển phi tuyến nâng cao, Truyền động điện hiện đại và Mô hình hóa hệ thống cơ điện tử phức tạp.
  • Kỹ sư R&D công nghiệp (Industrial R&D Engineers): Khung kiến trúc phần mềm điều khiển nhúng thời gian thực hoàn chỉnh viết bằng C/C++ có thể nhúng trực tiếp vào các dòng vi điều khiển DSP chuyên dụng (TI C2000 family như TMS320F28379D, TMS320F2808) để thương mại hóa sản phẩm.
  • Nhà hoạch định chính sách & Quản lý công nghệ: Cơ sở khoa học để xây dựng tiêu chuẩn kỹ thuật quốc gia về thiết bị hiệu suất cao, tiết kiệm năng lượng và các hệ thống phát điện tuabin khí/gió không phát thả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 độc đáo nhất là việc thiết lập thuật toán Backstepping cải tiến kết hợp khâu tích phân sai số và bộ lọc Dynamic Surface Control (DSC) tích hợp bộ quan sát High-Gain có đổi biến phụ $\xi$. Luận án đã mở rộng lý thuyết Backstepping của Krstic et al. (1995) từ việc chỉ áp dụng cho hệ giải tích thuần túy sang hệ cơ điện tử 2-DOF thực tế có tham số biến thiên $L(g)$ và nhiễu ghép kênh phi tuyến bậc hai $(i_d^2, i_q^2)$, đồng thời giải quyết triệt để vấn đề "bùng nổ toán hạng" và triệt tiêu hoàn toàn sai số tĩnh $\Delta z = 0$.

2. Đột phá về phương pháp luận nghiên cứu khi so sánh với ít nhất 2 công trình quốc tế?

So với nghiên cứu của Zhang et al. (2016) (chỉ mô hình hóa 1-DOF và bỏ qua biến thiên khe hở không khí) và Chiba et al. (2005) (dùng mạng nơ-ron phức tạp khó thực thi nhúng), luận án đã:

  • Xây dựng mô hình toán học giải tích phi tuyến toàn phần 2-DOF trên hệ tọa độ $dq$ mô tả chính xác tương tác lực - momen theo độ lệch $z$.
  • Ứng dụng thành công hệ thống mô phỏng thời gian thực HIL (Typhoon HIL402 kết hợp DSP phần cứng TI TMS320F28379D) với chu kỳ tính toán $1\ \mu\text{s}$ và tần số PWM 20 kHz, chứng minh tính khả thi tuyệt đối khi triển khai trên phần cứng công nghiệp thực tế.

3. Phát hiện bất ngờ nhất từ kết quả thực nghiệm dữ liệu là gì?

Phát hiện bất ngờ nhất là tính ưu việt thực tế của bộ quan sát trượt SMO so với bộ quan sát High-Gain trong môi trường có sóng hài đóng cắt nghịch lưu. Dù HGO có nền tảng toán học hội tụ mũ lý tưởng, nhưng khi chạy thời gian thực với biến tần PWM 20 kHz, SMO kết hợp khâu hiệu chỉnh góc pha lại cho đáp ứng sức điện động $\hat{e}{s\alpha}, \hat{e}{s\beta}$ ít gợn sóng hài hơn, giúp sai số góc rotor đạt $\Delta\theta_e \le 1.2^\circ$, vượt trội so với mức sai số $\approx 3.5^\circ$ của HGO.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) rõ ràng không?

Có. Luận án cung cấp chi tiết toàn bộ bảng thông số máy điện ($R_s, L_{sd0}, L_{sq0}, \lambda_m, g_0, m, J, p_c$), phương trình toán học hàm truyền, sơ đồ khối cấu trúc điều khiển, thuật toán biến đổi biến trạng thái phụ, sơ đồ đấu nối chân I/O giữa card Typhoon HIL402 và DSP TMS320F28379D, cùng cấu hình chi tiết phần mềm Typhoon HIL Control Center, đảm bảo khả năng tái lập kết quả $100%$.

5. Chương trình nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?

Chương trình nghiên cứu 10 năm tới định hướng: (1) Phát triển giải pháp điều khiển không cảm biến toàn dải tốc độ kết hợp tiêm tín hiệu cao tần HFI và Back-EMF SMO; (2) Ứng dụng kỹ thuật AI/Machine Learning để ước lượng thích nghi nhiệt độ và bù điện trở stator trực tuyến; (3) Mở rộng hệ điều khiển tích hợp trọn vẹn 5 bậc tự do (5-DOF) không cảm biến cho các hệ thống tuabin khí công suất siêu lớn; (4) Chế tạo và thử nghiệm công nghiệp trên các dòng máy nén ly tâm công nghiệp tốc độ trên $50.000\text{ vòng/phút}$.

Kết luận

Công trình nghiên cứu tiến sĩ của Ngô Mạnh Tùng (2022) đại diện cho một đóng góp học thuật xuất sắc, toàn diện và có tính ứng dụng thực tiễn cao trong ngành Kỹ thuật Điều khiển và Tự động hóa với 5 đóng góp cốt lõi:

  1. Thiết lập mô hình toán học giải tích phi tuyến toàn phần 2-DOF của động cơ đồng bộ kích thích vĩnh cửu từ trường dọc trục cấu trúc Sandwich, phản ánh chính xác sự phụ thuộc phi tuyến của điện cảm và lực từ vào độ biến thiên khe hở không khí $g(t) = g_0 \pm z(t)$.
  2. Đề xuất luật điều khiển phân bổ dòng điện vi sai ($i_{d1} = i_{d0} - i_d; i_{d2} = i_{d0} + i_d; i_{q1} = i_{q2} = i_q$) giúp tách kênh độc lập giữa mạch vòng lực nâng dọc trục và mạch vòng momen quay.
  3. Phát triển thuật toán Backstepping cải tiến tích hợp bộ lọc DSC và khâu tích phân sai số, giải quyết triệt để hiện tượng bùng nổ toán hạng và loại bỏ hoàn toàn sai số vị trí xác lập.
  4. Thiết kế bộ quan sát High-Gain với kỹ thuật đổi biến trạng thái phụ ($\xi_1, \xi_2$) ước lượng chính xác thành phần bất định và nhiễu tải $d_1, d_2$ mà không gây khuếch đại nhiễu đo tần số cao.
  5. Hiện thực hóa và kiểm chứng thành công chiến lược điều khiển không cảm biến tốc độ dựa trên bộ quan sát trượt SMO và High-Gain ước lượng sức điện động cảm ứng trên hệ thống phần cứng thời gian thực Typhoon HIL402 kết hợp DSP TI TMS320F28379D ở tần số PWM 20 kHz.

Nghiên cứu tạo tiền đề vững chắc cho việc loại bỏ cảm biến cơ học trong các hệ thống truyền động ổ từ siêu cao tốc, thúc đẩy tiến bộ công nghệ trong các ngành công nghiệp mũi nhọn như sản xuất bán dẫn, năng lượng sạch và hàng không vũ trụ.