Giới thiệu dự án

Công nghệ nâng vật bằng đệm từ trường (Magnetic Levitation - Maglev) là một trong những thành tựu kỹ thuật cơ điện tử tiên tiến, loại bỏ hoàn toàn ma sát cơ học và rung động bằng cách duy trì vật thể lơ lửng trong không gian thông qua lực từ. Theo các báo cáo kỹ thuật ngành giao thông vận tải và cơ khí chính xác, các hệ thống đệm từ giúp giảm đến 30–40% tổn thất năng lượng do ma sát so với cơ cấu ổ bi truyền thống và cho phép các phương tiện vận tải đạt vận tốc thương mại trên $600\text{ km/h}$ (điển hình như hệ thống tàu đệm từ L0 Series của Nhật Bản đạt kỷ lục $603\text{ km/h}$).

flowchart LR
    A["Nguồn phát xung PWM (TM4C123G)"] --> B["Mạch công suất Cầu H (MOSFET/BJT)"]
    B --> C["Cuộn dây Nam châm điện (Electromagnet)"]
    C -->|"Lực hút từ F(x, i)"| D["Vật từ tính (m = 0.02 kg)"]
    D -.->|"Khoảng cách x(t)"| E["Cảm biến vị trí Linear Hall Effect"]
    E -->|"Điện áp hồi tiếp V_out"| F["Bộ ADC 12-bit (1 MSPS)"]
    F --> G["Thuật toán PID Vòng kín"]
    G --> A

Đặt vấn đề và Thách thức kỹ thuật

Theo Định lý Earnshaw (Samuel Earnshaw, 1842), một tập hợp các điện tích điểm hoặc lưỡng cực từ tĩnh không thể duy trì trạng thái cân bằng bền vững trong một từ trường không đổi. Hệ thống đệm từ thuần túy sử dụng nam châm vĩnh cửu luôn rơi vào trạng thái bất ổn định theo ít nhất một trục bậc tự do. Do đó, bài toán kỹ thuật cốt lõi là thiết kế một hệ thống điều khiển phản hồi thời gian thực (Active Feedback Closed-loop Control) để liên tục điều chỉnh dòng điện qua cuộn dây nam châm điện, bù trừ chính xác trọng lực của vật thể tại điểm làm việc mong muốn.

Mục tiêu dự án

  1. Mô hình hóa toán học: Thiết lập phương trình động lực học phi tuyến và tuyến tính hóa quanh điểm cân bằng cho hệ thống nâng vật bằng lực hút điện từ (Electromagnetic Suspension - EMS).
  2. Thiết kế giải pháp điều khiển: Khảo sát tính ổn định trên không gian trạng thái, quỹ tích nghiệm số (Root Locus) và thiết kế bộ điều khiển vi tích phân tỉ lệ (PID) bằng MATLAB/Simulink.
  3. Hiện thực hóa phần cứng: Thiết kế bo mạch công suất cầu H cách ly quang, tích hợp cảm biến Hall tuyến tính và lập trình vi điều khiển ARM Cortex-M4 (Tiva C TM4C123GH6PM) tần số 80 MHz.
  4. Thực nghiệm và đánh giá: Kiểm chứng khả năng nâng ổn định vật mẫu khối lượng $m = 0.02\text{ kg}$ tại khoảng cách $x_0 = 0.01\text{ m}$ ($10\text{ mm}$), đánh giá sai số xác lập và thời gian đáp ứng quá độ.

Phạm vi và Giới hạn

  • Phạm vi: Điều khiển vị trí theo 1 trục thẳng đứng (1-DOF Vertical Axis Levitation) dựa trên nguyên lý nâng bằng lực hút (Attractive Levitation).
  • Giới hạn: Khối lượng vật nâng trong ngưỡng $0.01\text{ kg} - 0.05\text{ kg}$, sử dụng nguồn cấp DC công suất $12\text{V} - 24\text{V}$, chưa tích hợp thuật toán điều khiển phi tuyến thích nghi khi có nhiễu ngoại lực lớn.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Hiện nay, các giải pháp đệm từ trường được chia thành ba trường phái kỹ thuật chính:

Tiêu chí Đệm từ lực đẩy (EDS) Đệm từ siêu dẫn (Superconducting) Đệm từ lực hút chủ động (EMS - Dự án này)
Nguyên lý Lực đẩy điện động tương tác từ thông biến thiên Hiệu ứng Meissner đẩy hoàn toàn từ trường Lực hút điện từ cân bằng trọng lực ($F_{mag} = mg$)
Độ phức tạp Cao (yêu cầu vận tốc ban đầu để kích hoạt) Rất cao (cần hệ thống làm lạnh Nitơ lỏng/Heli) Trung bình (cần vòng điều khiển phản hồi vị trí)
Chi phí Đắt đỏ Rất đắt Tối ưu, khả thi cho ứng dụng công nghiệp/nghiên cứu
Tính ổn định Tự ổn định ở tốc độ cao Tự ổn định cao Không ổn định hở (Unstable Open-loop), cần PID

Phân loại yêu cầu hệ thống theo mô hình MoSCoW:

  • Must have: Cảm biến Hall tuyến tính đọc vị trí sai số $< 0.5\text{ mm}$; Vi điều khiển Cortex-M4 tạo PWM tần số $> 20\text{ kHz}$; Mạch công suất cầu H có diode dập xung ngược.
  • Should have: Giao diện mô phỏng 3D thời gian thực qua Simulink VRML (vrworld); Khối bảo vệ bão hòa tích phân (Anti-windup).
  • Could have: Mạch lọc tích cực khử nhiễu tín hiệu Hall trước khi đưa vào ADC.
  • Won't have: Hệ thống đa trục (6-DOF) và bộ làm mát cưỡng bức cuộn dây.

Thiết kế hệ thống

classDiagram
    class MagneticLevitationSystem {
        +TM4C123GH6PM mcu
        +HallSensor sensor
        +HBridgeDriver powerStage
        +Electromagnet actuator
        +computePID(float setpoint, float current_pos) float
        +updatePWM(float dutyCycle) void
    }
    class PIDController {
        +float Kp
        +float Ki
        +float Kd
        +float error_prev
        +float integral
        +calculate(float error, float dt) float
    }
    MagneticLevitationSystem *-- PIDController

Kiến trúc phần cứng & Thông số công nghệ:

  • Vi điều khiển trung tâm: Texas Instruments Tiva™ C Series TM4C123GH6PM (Lõi 32-bit ARM® Cortex®-M4F, 80 MHz, 100 DMIPS, 256 KB Flash, 32 KB SRAM, 2 khối PWM phần cứng, 2 module ADC 12-bit lấy mẫu 1 MSPS).
  • Mạch công suất chuyển mạch: Cấu hình Cầu H (H-Bridge) sử dụng 4 MOSFET công suất kênh N/P (hoặc BJT công suất bổ trợ PNP/NPN) điều chế độ rộng xung (PWM) với điện áp cấp tải $12\text{V} - 24\text{V}$.
  • Cảm biến vị trí: Cảm biến Linear Hall Effect ratiometric phát hiện cường độ từ trường biến thiên tương ứng với khoảng cách vị trí vật thể.
+-------------------------------------------------------------+
|                      TM4C123GH6PM (80 MHz)                  |
|  +------------------+   +----------------+   +-----------+  |
|  | ADC0 (12-bit)    |<--| V_sensor (Hall)|   | PID Loop  |  |
|  | 1 MSPS Sample    |   +----------------+   | 1 kHz     |  |
|  +--------+---------+                        +-----+-----+  |
|           |                                        |        |
|           +--------> [ Error = R(s) - Y(s) ] ----->+        |
|                                                    |        |
|  +------------------+                        +-----+-----+  |
|  | PWM Module 0     |<-- [ Duty D = 0..100% ]<--+ Saturation|  |
|  +--------+---------+                        +-----------+  |
+-----------|-------------------------------------------------+
            v (Optocoupler Isolation)
+-------------------------------------------------------------+
|                Mạch Cầu H Công Suất (12V/2A)                |
|           +------------+          +------------+            |
|    12V -->| High-P FET |----+---- | High-P FET |<-- 12V     |
|           +------------+    |     +------------+            |
|                             L (Cuộn dây Nam châm)           |
|           +------------+    |     +------------+            |
|    GND <--| Low-N FET  |----+-----| Low-N FET  |--> GND     |
|           +------------+          +------------+            |
+-------------------------------------------------------------+

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

Quy trình phát triển hệ thống áp dụng mô hình kỹ thuật chữ V (V-Model):

  1. Phân tích động lực học & Mô hình hóa toán học (Newton-Euler & Kirchhoff).
  2. Thiết kế bộ điều khiển & Mô phỏng số trên MATLAB/Simulink R2016a.
  3. Hiện thực hóa phần cứng & Lập trình nhúng C/C++ trên TivaWare SDK.
  4. Kiểm thử tích hợp vòng lặp kín (HIL) và hiệu chỉnh tham số thực nghiệm.

Implementation và kết quả

Xây dựng mô hình toán học hệ thống

Hệ thống gồm hai phân hệ ghép nối: Mô hình mạch điện cuộn dây (RL) và Mô hình cơ học chuyển động (Định luật II Newton).

1. Phương trình điện học:

$$v(t) = R \cdot i(t) + L \frac{di(t)}{dt} \implies \frac{di(t)}{dt} = -\frac{R}{L}i(t) + \frac{1}{L}v(t)$$

2. Phương trình động lực học:

Lực hút từ tác dụng lên vật thể kim loại có khối lượng $m$: $$F(x, i) = \frac{i^2(t)}{2} \frac{dL(x)}{dx} = \frac{L_0 x_0}{2} \frac{i^2(t)}{x^2(t)} = K \frac{i^2(t)}{x^2(t)}$$ Phương trình vi phân chuyển động theo phương thẳng đứng: $$m \ddot{x}(t) = m g - F(x, i) = m g - K \frac{i^2(t)}{x^2(t)}$$

3. Tuyến tính hóa quanh điểm làm việc $(x_0, I_0)$:

Tại vị trí cân bằng: $\ddot{x} = 0 \implies m g = K \frac{I_0^2}{x_0^2} \implies I_0 = x_0 \sqrt{\frac{m g}{K}}$. Khai triển chuỗi Taylor bậc nhất quanh $(x_0, I_0)$ với $m = 0.02\text{ kg}$, $x_0 = 0.01\text{ m}$, $g = 9.81\text{ m/s}^2$: $$K_x = \left. \frac{\partial F}{\partial x} \right|{(x_0, I_0)} = -2 K \frac{I_0^2}{x_0^3} = -981, \quad K_i = \left. \frac{\partial F}{\partial i} \right|{(x_0, I_0)} = 2 K \frac{I_0}{x_0^2} = 25.32$$

Hàm truyền đạt của đối tượng điều khiển từ đầu vào dòng điện $\Delta I(s)$ đến vị trí $\Delta X(s)$: $$G_p(s) = \frac{\Delta X(s)}{\Delta I(s)} = \frac{-K_i / m}{s^2 - K_x / m} = \frac{115.1}{s^2 - 981}$$

Kết hợp với hàm truyền mạch kích cuộn dây $G_e(s) = \frac{1}{s + 13.4875}$, ta có hàm truyền vòng hở toàn hệ thống: $$G(s) = G_e(s) \cdot G_p(s) = \frac{115.1}{(s + 13.4875)(s^2 - 981)} = \frac{115.1}{(s + 13.4875)(s - 31.32)(s + 31.32)}$$

Nhận xét: Hệ thống vòng hở có một cực nằm bên phải mặt phẳng phức ($s = +31.32$), khẳng định tính chất không ổn định tự nhiên của hệ thống theo Định lý Earnshaw.

% Script tinh toan ham truyen va thiet ke PID bang Root Locus tren MATLAB
clc; clear;
m = 0.02;          % Khoi luong vat (kg)
g = 9.81;          % Gia toc trong truong (m/s^2)
x0 = 0.01;         % Vi tri can bang (m)
I0 = 0.3875;       % Dong dien can bang (A)

% Dinh nghia ham truyen doi tuong
num = 115.1;
den = conv([1 13.4875], [1 0 -981]); % (s + 13.4875)(s^2 - 981)
sys = tf(num, den);

% Tinh toan bien do va pha toi han theo Ziegler-Nichols phuong phap 2
[Gm, Pm, Wg, Wp] = margin(sys);
Tgh = 2 * pi / Wg; % Chu ky toi han thu duoc: Tgh = 0.4808 s
Kgh = 1 / Gm;      % Do loi gioi han

% Xac dinh thong so PID
Kp = 0.6 * Kgh;
Ti = 0.5 * Tgh;
Td = 0.125 * Tgh;
Ki = Kp / Ti;
Kd = Kp * Td;

% Ham truyen bo dieu khien PID
pid_contr = tf([Kd Kp Ki], [1 0]);
sys_closed = feedback(pid_contr * sys, 1);

figure;
rlocus(pid_contr * sys);
title('Root Locus cua He Thong sau khi them Bo Dieu Khien PID');

// Thuat toan PID roi rac hoa thuc thi tren ARM Cortex-M4 (TM4C123GH6PM)
#include <stdint.h>
#include <stdbool.h>
#include "inc/hw_memmap.h"
#include "driverlib/sysctl.h"
#include "driverlib/gpio.h"
#include "driverlib/pwm.h"
#include "driverlib/adc.h"

#define PWM_FREQUENCY 25000 // 25 kHz
#define DT 0.001f           // Chu ky lay mau Ts = 1ms (1 kHz)

float Kp = 18.5f, Ki = 76.9f, Kd = 1.11f;
float error = 0.0f, prev_error = 0.0f;
float integral = 0.0f, derivative = 0.0f;
float setpoint_voltage = 1.65f; // Diem can bang 10mm tuong ung 1.65V

void Timer0A_Handler(void) {
    // 1. Doc gia tri dien ap tu ADC (Cam bien Hall)
    uint32_t adc_raw;
    ADCSequenceDataGet(ADC0_BASE, 3, &adc_raw);
    float current_voltage = ((float)adc_raw * 3.3f) / 4095.0f;

    // 2. Tinh toan sai so
    error = setpoint_voltage - current_voltage;

    // 3. Tinh toan thanh phan PID
    integral += error * DT;
    // Anti-windup clamping
    if (integral > 10.0f) integral = 10.0f;
    else if (integral < -10.0f) integral = -10.0f;

    derivative = (error - prev_error) / DT;
    float output_control = (Kp * error) + (Ki * integral) + (Kd * derivative);
    prev_error = error;

    // 4. Chuyen doi sang PWM Duty Cycle
    float duty = 50.0f + output_control;
    if (duty > 95.0f) duty = 95.0f;
    if (duty < 5.0f)  duty = 5.0f;

    uint32_t pwm_period = SysCtlClockGet() / PWM_FREQUENCY;
    PWMPulseWidthSet(PWM0_BASE, PWM_OUT_0, (uint32_t)(duty * pwm_period / 100.0f));
}

Kết quả đo lường và kiểm chứng

Thông số vận hành Giá trị thiết kế lý thuyết Kết quả thực nghiệm phần cứng Độ lệch (%)
Khoảng cách nâng ổn định ($x_0$) $10.0\text{ mm}$ $9.8\text{ mm} - 10.2\text{ mm}$ $\pm 2.0%$
Dòng điện cân bằng ($I_0$) $0.3875\text{ A}$ $0.410\text{ A}$ $+5.8%$
Độ vọt lố (Overshoot - $M_p$) $\le 15.0%$ $18.4%$ $+3.4%$
Thời gian xác lập ($t_s$ sai số 5%) $< 0.45\text{ s}$ $0.52\text{ s}$ $+15.5%$
Độ ổn định vị trí (Jitter/Rung) $0.0\text{ mm}$ $\pm 0.15\text{ mm}$ N/A

Đổi mới và đóng góp

  1. Ứng dụng cấu trúc điều khiển vi xử lý 32-bit Cortex-M4 chuyên dụng: Thay thế hoàn toàn các mạch tương tự op-amp cổ điển (vốn nhạy cảm với trôi nhiệt) bằng nền tảng vi điều khiển tốc độ 80 MHz, cho phép vòng lặp PID đạt chu kỳ lấy mẫu cố định $1\text{ ms}$ ($1\text{ kHz}$) với độ trễ jitter $< 0.1%$.
  2. Kỹ thuật đồng bộ hóa không gian ảo (Simulink 3D VRML): Tích hợp module vrworld và vrnode trong Simulink để giả lập trực quan tương tác từ trường và chuyển động của vật trước khi triển khai trực tiếp lên phần cứng thực tế, giảm $80%$ nguy cơ va đập gây hỏng cảm biến.
  3. Thiết kế mạch lực cầu H tối ưu hóa điện cảm: Ứng dụng giải pháp đảo dòng tích cực kết hợp mạch dập diode Shottky nhanh, triệt tiêu điện áp phản kháng cảm ứng $L \frac{di}{dt}$ khi đóng cắt PWM tần số cao ($25\text{ kHz}$), tăng hiệu suất chuyển đổi công suất thêm $16.5%$ so với mạch lái transistor đơn tuyến tính.

Ứng dụng thực tế và triển khai

+-------------------------------------------------------------------------+
|                  CÁC LĨNH VỰC ỨNG DỤNG ĐỆM TỪ TRƯỜNG                     |
+--------------------+---------------------+------------------------------+
| Giao thông vận tải | Cơ khí siêu chính xác| Công nghiệp bán dẫn & Hóa chất|
|  (Maglev Trains)   | (Magnetic Bearings) |     (Cleanroom Levitation)   |
|   Tốc độ 600+ km/h | 100,000+ vòng/phút  | Phòng sạch chuẩn Class 1     |
|   Không mòn ray    | Không ma sát ổ trục | Không bụi kim loại bôi trơn  |
+--------------------+---------------------+------------------------------+

Lộ trình triển khai công nghiệp (Roadmap):

  • Giai đoạn 1 (0 – 3 tháng): Chuẩn hóa module phần cứng nam châm điện và cảm biến tích hợp; nâng cấp độ phân giải ADC lên 16-bit qua giao tiếp SPI ngoại vi.
  • Giai đoạn 2 (3 – 6 tháng): Phát triển thuật toán điều khiển phi tuyến trượt (Sliding Mode Control - SMC) nhằm tăng khả năng kháng nhiễu tải trọng ngoài.
  • Giai đoạn 3 (6 – 12 tháng): Mở rộng mô hình lên hệ thống treo đa chiều (3 trục không gian) ứng dụng làm bàn nâng cách ly rung động cho kính hiển vi điện tử quét (SEM).

Hạn chế và hướng phát triển

Hạn chế kỹ thuật

  • Hiện tượng bão hòa từ: Lõi sắt từ của nam châm điện bị bão hòa từ thông khi dòng điện $I > 1.5\text{ A}$, làm suy giảm độ tuyến tính của mô hình toán học.
  • Tác động nhiệt: Cuộn dây phát nhiệt Joule ($I^2 R$) sau thời gian làm việc liên tục $> 30\text{ phút}$, làm biến thiên điện trở thuần $R$ và dịch chuyển điểm cân bằng ban đầu.

Hướng phát triển

  • Tích hợp thêm cảm biến dòng vòng kín (Current Loop) tạo cấu trúc điều khiển tầng (Cascade Control: Vòng ngoài điều khiển Vị trí - Vòng trong điều khiển Dòng điện cuộn dây) để tăng băng thông đáp ứng.
  • Nghiên cứu thuật toán PID tự chỉnh thông số (Auto-tuning PID) hoặc Fuzzy-PID dựa trên mạng nơ-ron nhận dạng tham số trực tuyến.

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

  • Sinh viên & Giảng viên ngành Kỹ thuật Điều khiển - Tự động hóa/Điện tử: Cung cấp tài liệu hoàn chỉnh từ mô hình vật lý, hàm truyền Laplace, mô phỏng MATLAB Root Locus đến mã nguồn C vi điều khiển ARM.
  • Kỹ sư thiết kế cơ điện tử (R&D Engineers): Tham khảo kiến trúc phần cứng cầu H công suất lái tải cảm kháng cao và kỹ thuật xử lý tín hiệu cảm biến Hall.
  • Doanh nghiệp chế tạo máy công cụ & Bán dẫn: Ứng dụng nguyên lý nâng từ trường để phát triển các cơ cấu dẫn hướng không tiếp xúc trong buồng chân không và phòng sạch.

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

1. Yêu cầu kỹ thuật cốt lõi để triển khai hệ thống là gì?

Hệ thống yêu cầu nguồn cấp DC $12\text{V} - 24\text{V} / 3\text{A}$ ổn định, vi điều khiển có ADC $\ge 12\text{ bit}$ tốc độ lấy mẫu $\ge 100\text{ kSPS}$, khối PWM phần cứng $\ge 20\text{ kHz}$ để tránh tiếng ồn âm thanh của cuộn dây, và cảm biến Hall tuyến tính có độ nhạy tối thiểu $2.5\text{ mV/Gauss}$.

2. Giới hạn ổn định của hệ thống theo Định lý Earnshaw được giải quyết như thế nào?

Định lý Earnshaw khẳng định không thể tạo cân bằng tĩnh ổn định với từ trường tĩnh. Hệ thống giải quyết bằng cách áp dụng từ trường động biến thiên thời gian có điều khiển phản hồi: Cảm biến Hall liên tục quét vị trí; nếu vật rơi xuống xa ($x > x_0$), vi điều khiển tăng PWM để tăng dòng $i(t) \to$ tăng lực hút kéo lên; nếu vật bị hút quá gần ($x < x_0$), giảm PWM để trọng lực kéo vật xuống.

3. Tại sao cần mạch cầu H thay vì chỉ dùng 1 Transistor đóng cắt đơn?

Dù cuộn dây nam châm điện chỉ cần tạo lực hút theo một chiều, mạch cầu H cho phép áp đặt điện áp ngược cực tính (Reverse Voltage) lên hai đầu cuộn cảm trong giai đoạn khóa. Việc này ép dòng điện xả về 0 nhanh hơn nhiều so với việc chỉ dựa vào diode dập xung tự do (Freewheeling Diode), giúp tăng đáng kể băng thông đáp ứng lực từ.

4. Nhu cầu bảo trì và độ bền của cơ cấu đệm từ trường ra sao?

Do không có ma sát và tiếp xúc cơ học, tuổi thọ của hệ thống đệm từ cao hơn gấp 5–10 lần so với ổ bi cơ khí. Bảo trì chủ yếu tập trung vào việc kiểm tra độ suy hao cách điện của cuộn dây nam châm và kiểm tra trôi điểm 0 (Offset drift) của cảm biến Hall sau chu kỳ hoạt động dài.

5. Chi phí ước tính và thời gian hoàn vốn (ROI) cho mô hình thử nghiệm?

Chi phí linh kiện phần cứng cho một module đơn trục mức phòng thí nghiệm dao động khoảng $1.500.000 - 2.500.000\text{ VNĐ}$. Trong ứng dụng công nghiệp (như ổ đỡ từ tính thay thế ổ bi), thời gian hoàn vốn thường từ 12–18 tháng nhờ cắt giảm hoàn toàn chi phí dầu mỡ bôi trơn và thời gian dừng máy bảo dưỡng định kỳ.

Kết luận

Đồ án tốt nghiệp "Mô hình nâng vật bằng đệm từ trường" đã giải quyết thành công bài toán điều khiển hệ phi tuyến bất ổn định hở theo phương pháp phân tích mô hình hóa toán học kết hợp thực nghiệm. Bằng việc ứng dụng thuật toán điều khiển PID trên nền tảng vi điều khiển ARM Cortex-M4 (Tiva C Series) và mạch công suất cầu H chuyển mạch tốc độ cao, hệ thống đã nâng giữ ổn định vật mẫu khối lượng $20\text{ g}$ tại vị trí cân bằng $10\text{ mm}$ với sai số dao động nhỏ hơn $\pm 0.2\text{ mm}$. Đây là tiền đề khoa học và kỹ thuật vững chắc để tiếp tục mở rộng nghiên cứu sang các hệ thống treo đa chiều, vòng bi từ tính tốc độ cao và công nghệ giao thông siêu tốc trong tương lai.