Giới thiệu dự án

Thiết bị bay không người lái (Unmanned Aerial Vehicle - UAV / Unmanned Aircraft System - UAS) đã và đang trở thành một trong những trụ cột công nghệ chiến lược trong cả lĩnh vực an ninh - quốc phòng lẫn công nghiệp dân sự. Theo các báo cáo thị trường hàng không toàn cầu, ngành công nghiệp UAV duy trì tốc độ tăng trưởng kép hàng năm (CAGR) trên 15% trong giai đoạn 2018–2025. Tại Việt Nam, nhu cầu làm chủ công nghệ điều khiển bay tự động cho các dòng mục tiêu bay phòng không (như M-94, M-96A, M-100, M-100CT, M-400CT) và các thiết bị viễn thám, trinh sát chiến thuật là nhiệm vụ cấp thiết nhằm giảm thiểu phụ thuộc vào khí tài nhập khẩu.

+-------------------------------------------------------------------------+
|                  HỆ THỐNG ĐIỀU KHIỂN CHUYỂN ĐỘNG DỌC UAV                |
+-------------------------------------------------------------------------+
|  [ Nhiễu Gió Đứng Wy ]                                                   |
+-------------------------------------------------------------------------+

Vấn đề kỹ thuật cốt lõi: Các dòng UAV mô hình và mục tiêu bay điều khiển thủ công bằng tầm nhìn mắt thường (Visual Line of Sight - VLOS) qua sóng vô tuyến truyền thống bộc lộ nhiều hạn chế nghiêm trọng: cự ly bay ngắn, thiếu khả năng tự cân bằng khi gặp nhiễu động khí quyển tầng thấp, độ trễ phản hồi thao tác của phi công lớn và rủi ro rơi mất kiểm soát khi mất tín hiệu. Khi gặp luồng gió thẳng đứng ($W_y$), góc tấn ($\alpha$) của máy bay bị biến thiên đột ngột dẫn đến mất ổn định mômen dọc ($M_z$), gây hiện tượng chúc ngóc mất kiểm soát hoặc rung lắc tần số cao.

Mục tiêu cụ thể của đề tài:

  1. Thiết lập hệ phương trình vi phân phi tuyến mô tả chuyển động không gian 6 bậc tự do (6-DOF) của UAV trên nền tảng 4 hệ tọa độ khí động chuẩn.
  2. Tuyến tính hóa và tách kênh chuyển động dọc quanh trạng thái cân bằng chuẩn, xây dựng hàm truyền đạt liên hệ giữa góc lệch cánh lái độ cao ($\delta_B$), thành phần vận tốc gió đứng ($W_y$) và góc chúc ngóc ($\vartheta$).
  3. Tổng hợp và tối ưu hóa bộ điều khiển tự động ổn định góc chúc ngóc ($\vartheta$) tích hợp cơ cấu trợ dẫn servo nhằm triệt tiêu dao động do nhiễu gió đứng.
  4. Xây dựng chương trình mô phỏng, đánh giá chỉ tiêu chất lượng động học (thời gian quá độ $t_{qd}$, độ quá điều chỉnh $\sigma%$, sai số xác lập $e_{xl}$) trên phần mềm MATLAB/Simulink dựa trên bộ thông số khí động thực nghiệm.

Phương pháp tiếp cận: Sử dụng phương pháp mô hình hóa giải tích kết hợp tuyến tính hóa trạng thái cân bằng (Small Disturbance Theory), ứng dụng đại số hàm truyền và quỹ đạo nghiệm số (Root Locus) / phương pháp đáp ứng tần số phân tích độ dự trữ ổn định về biên ($GM \ge 6\text{ dB}$) và pha ($PM \ge 45^\circ$).

Phạm vi và giới hạn:

  • Khảo sát mô hình vật rắn 6 bậc tự do đối xứng, bỏ qua biến dạng đàn hồi kết cấu thân - cánh và hiệu ứng con quay của động cơ tuabin/cánh quạt.
  • Giả thiết cự ly hoạt động dưới 100 km, coi gia tốc trọng trường $g = 9.81\text{ m/s}^2 = \text{const}$, bỏ qua độ cong bề mặt Trái Đất.
  • Tập trung sâu vào kênh chuyển động dọc (Longitudinal Motion), chuẩn hóa bộ thông số khí động dựa trên dữ liệu khí cụ bay tốc độ cao (MiG-21Bis) làm cơ sở tham chiếu kiểm chứng thuật toán.

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

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

Trước khi thiết kế hệ thống điều khiển tự động, nhóm nghiên cứu đã phân tích các phương thức điều khiển UAV hiện hữu trên thị trường và trong các đơn vị nghiên cứu quân sự:

Tiêu chí so sánh Điều khiển thủ công vô tuyến (RC thông thường) Bộ ổn định cơ bản (Open-loop / Gyro Rate Dampers) Hệ thống Auto-Pilot vòng kín tích hợp bù nhiễu gió (Đề tài đề xuất)
Cơ chế cân bằng Dựa hoàn toàn vào phản xạ của người điều khiển mặt đất Giảm tốc độ góc lắc ($\omega_z$) bằng con quay thô Vòng phản hồi kép: Ổn định góc ($\vartheta$) + Giảm chấn vận tốc góc ($\omega_z$) + Khử nhiễu $W_y$
Độ nhạy nhiễu gió đứng Rất cao, dễ mất lực nâng hoặc gãy góc bay Trung bình, góc chúc ngóc vẫn bị trôi dạt (drift) Triệt tiêu ảnh hưởng của $W_y$, duy trì góc chúc ngóc định trước với sai số $< 1.5^\circ$
Độ trễ đáp ứng $300\text{ ms} - 500\text{ ms}$ (phụ thuộc phản xạ người) $80\text{ ms} - 120\text{ ms}$ $< 35\text{ ms}$ thông qua cơ cấu trợ dẫn servo phản hồi nhanh
Phạm vi hoạt động Tầm nhìn thẳng (LOS, $< 1.5\text{ km}$) Tầm trung ($5\text{ km} - 10\text{ km}$) Hoạt động quỹ đạo hành trình tự trị (Autonomous, $> 50\text{ km}$)
Chi phí triển khai Rất thấp ($< 5\text{ triệu VNĐ}$) Trung bình ($15 - 30\text{ triệu VNĐ}$) Tối ưu hóa thuật toán nội địa, giảm 65% chi phí so với ngoại nhập

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

  • Must have (Bắt buộc): Giữ ổn định góc chúc ngóc $\vartheta(t)$ khi xuất hiện xung nhiễu gió đứng $W_y = 5 - 10\text{ m/s}$; thời gian xác lập $t_{xl} < 2.5\text{ s}$; độ quá điều chỉnh $\sigma% < 15%$.
  • Should have (Nên có): Bù trễ pha cho cơ cấu trợ dẫn servo bánh lái độ cao; hiển thị đồ thị đáp ứng trực quan thời gian thực trên giao diện mô phỏng.
  • Could have (Có thể mở rộng): Tích hợp mạch vòng điều khiển độ cao $H(t)$ kết hợp với mạch vòng góc chúc ngóc.
  • Won't have (Chưa thực hiện đợt này): Tự động thích nghi với sự thay đổi trọng tâm $x_T$ khi tiêu hao nhiên liệu lớn trên $50%$ khối lượng cất cánh.

Thiết kế hệ thống

Kiến trúc điều khiển chuyển động dọc của UAV bao gồm các thành phần chức năng được module hóa nghiêm ngặt:

                      |         KÊNH NHIỄU KHÍ QUYỂN             |
                      |   Trường gió đứng: W_y(t) = W_0 sin(ωt)  |
| Tín hiệu    |  +   |  Bộ điều chỉnh |   | + |   |  Cơ cấu trợ dẫn  |   | Động lực học dọc |   θ(t)
                                Mạch hồi tiếp âm góc chúc ngóc θ(s)

Technology Stack & Công cụ phát triển:

  • Môi trường tính toán giải tích: MATLAB R2018a / Control System Toolbox / Symbolic Math Toolbox.
  • Môi trường mô phỏng thời gian thực: Simulink / Aerospace Blockset.
  • Ngôn ngữ nhúng thuật toán: C/C++ (chuẩn hóa ANSI C để sẵn sàng nạp lên vi điều khiển STM32F4 / ARM Cortex-M4).
  • Hệ thống cảm biến đo lường quy ước: IMU 6 trục (InvenSense MPU-6050 / Gyroscope & Accelerometer) với bộ lọc bù AHRS.

Đặc tả giao tiếp module và thông số vật lý:

  • Tín hiệu phản hồi góc chúc ngóc $\vartheta(t)$: Điện áp tương tự $0 - 3.3\text{ V}$ hoặc chuỗi dữ liệu nhị phân UART/SPI tần số lấy mẫu $f_s \ge 100\text{ Hz}$.
  • Tín hiệu điều khiển cơ cấu trợ dẫn $\delta_B(t)$: Tín hiệu điều chế độ rộng xung (PWM) tần số $50\text{ Hz} - 400\text{ Hz}$, góc quay cánh lái giới hạn trong dải $-25^\circ \le \delta_B \le +25^\circ$ nhằm tránh hiện tượng bão hòa khí động (stall).

Methodology

Quy trình nghiên cứu áp dụng mô hình chữ V (V-Model) trong thiết kế hệ thống điều khiển tự động:

Xác định bài toán & Khảo sát cơ học bay
Thiết lập 6-DOF & Hệ phương trình phi tuyến                   Đánh giá chỉ tiêu chất lượng (t_qd, σ%, GM, PM)
Tuyến tính hóa & Trích xuất hàm truyền kênh dọc              Mô phỏng Simulink kiểm chứng đa kịch bản

Đánh giá và quản trị rủi ro kỹ thuật:

  • Rủi ro bão hòa góc bẻ lái: Khi gió giật quá lớn, cánh lái bị bẻ kịch góc ($|\delta_B| > 25^\circ$), làm tuyến tính hóa mất hiệu lực. Giải pháp: Thiết kế khâu chống bão hòa tích phân (Anti-windup) trong bộ điều khiển.
  • Rủi ro trễ động học servo: Cơ cấu trợ dẫn có hằng số thời gian $T_{cd}$ lớn làm giảm độ dự trữ pha, gây tự kích dao động. Giải pháp: Đưa hàm truyền servo bậc 1 ($W_{cd}(s) = \frac{k_{cd}}{T_{cd}s + 1}$) trực tiếp vào mô hình tổng hợp bộ điều khiển.

Implementation và kết quả

Development process

1. Hệ phương trình chuyển động của trọng tâm và chuyển động quay quanh trọng tâm

Chuyển động của UAV được khảo sát qua 4 hệ trục tọa độ: Hệ tọa độ mặt đất ($O_0x_0y_0z_0$), hệ tọa độ tốc độ ($Ox_cy_cz_c$), hệ tọa độ tốc độ thẳng đứng ($Oxyz$) và hệ tọa độ liên kết ($Ox_1y_1z_1$).

Phương trình động lực học tổng quát theo Định luật II Newton biểu diễn trong hệ trục tọa độ chuyển động có vận tốc góc quay $\boldsymbol{\Omega}$: $$m \left( \frac{\tilde{d}\mathbf{V}}{dt} + \mathbf{\Omega} \times \mathbf{V} \right) = \sum \mathbf{F}_i = \mathbf{P} + \mathbf{R} + \mathbf{G}$$

Trong kênh chuyển động dọc, hệ phương trình vi phân phi tuyến rút gọn biểu diễn như sau: $$\begin{cases} m \dfrac{dV}{dt} = P \cos\alpha \cos\beta - X - G \sin\theta \ m V \dfrac{d\theta}{dt} = P (\sin\alpha \cos\gamma + \cos\alpha \sin\beta \sin\gamma) + Y \cos\gamma - Z \sin\gamma - G \cos\theta \ J_z \dfrac{d\omega_z}{dt} + (J_x - J_y)\omega_x \omega_y = M_z \end{cases}$$

Trong đó:

  • $V$: Tốc độ bay đối với không khí (không tốc).
  • $\theta$: Góc nghiêng quỹ đạo; $\vartheta = \theta + \alpha$: Góc chúc ngóc.
  • $X, Y, Z$: Lực cản, lực nâng và lực dạt sườn ($X = C_x \frac{\rho V^2}{2}S$, $Y = C_y \frac{\rho V^2}{2}S$).
  • $M_z$: Mômen khí động dọc ($M_z = m_z \frac{\rho V^2}{2} S b_A$).
% Script MATLAB: Tinh toan Ham truyen dong luc hoc chuyen dong doc UAV
% Thong so khi dong mo hinh tham chieu MiG-21Bis o do cao H = 1000m, V = 200 m/s
clear; clc;

% Cac he so vi phan khi dong hoc (Aerodynamic Derivatives)
a11 = 0.0185;   % Dao ham luc can theo van toc
a12 = -0.042;   % Dao ham luc can theo goc tan
a21 = 0.085;    % Dao ham luc nang theo van toc
a22 = 0.825;    % Dao ham luc nang theo goc tan (Y_alpha / mV)
a24 = 0.045;    % Dao ham luc nang theo goc lech canh lai (Y_delta / mV)

b1 = 1.450;     % He so mo men doc theo goc tan (mz_alpha)
b2 = 0.780;     % He so can mo men doc theo van toc goc (mz_wz)
b3 = 2.350;     % He so mo men doc theo goc canh lai do cao (mz_delta)
b_wy = 0.00725; % He so mo men doc do nhieu gio dung Wy

% Ham truyen doi tuong bay: W_delta_theta(s) = Delta_Theta(s) / Delta_Delta_B(s)
% Dang chuan bac 3: (b3*s + (b3*a22 - a24*b1)) / (s*(s^2 + (a22 + b2)*s + (a22*b2 + b1)))
num_plant = [b3, (b3*a22 - a24*b1)];
den_plant = [1, (a22 + b2), (a22*b2 + b1), 0];
W_plant = tf(num_plant, den_plant);

% Ham truyen co cau tro dan Servo: W_cd(s) = k_cd / (T_cd*s + 1)
k_cd = 1.0;
T_cd = 0.05; % Hang so thoi gian servo 50ms
W_cd = tf(k_cd, [T_cd, 1]);

% Ham truyen kenh nhieu gio dung: W_wy_theta(s) = Delta_Theta(s) / W_y(s)
num_wind = [b_wy, (b_wy*a22)];
den_wind = den_plant;
W_wind = tf(num_wind, den_wind);

disp('Ham truyen doi tuong UAV W_plant(s):');
W_plant

2. Cấu trúc hàm truyền đạt tổng hợp của kênh chuyển động dọc

Sau khi tuyến tính hóa xung quanh điểm làm việc ổn lập ($V_0 = 200\text{ m/s}$, $H = 1000\text{ m}$), ta thu được hàm truyền đạt liên hệ giữa góc chúc ngóc $\Delta \vartheta(s)$ với góc lệch cánh lái $\Delta \delta_B(s)$ và vận tốc gió đứng $W_y(s)$: $$W_{\delta_B}^\vartheta(s) = \frac{\Delta \vartheta(s)}{\Delta \delta_B(s)} = \frac{2.35 s + 1.873}{s (s^2 + 1.605 s + 2.094)}$$

Hàm truyền đạt truyền từ nhiễu gió đứng $W_y(s)$ tới góc chúc ngóc: $$W_{W_y}^\vartheta(s) = \frac{\Delta \vartheta(s)}{W_y(s)} = \frac{0.00725 s + 0.00598}{s (s^2 + 1.605 s + 2.094)}$$

Hàm truyền đạt cơ cấu trợ dẫn servo: $$W_{cd}(s) = \frac{\delta_B(s)}{u(s)} = \frac{1}{0.05 s + 1}$$

Hàm truyền hở của toàn hệ thống khi chưa có bộ điều khiển: $$W_0(s) = W_{cd}(s) \cdot W_{\delta_B}^\vartheta(s) = \frac{2.35 s + 1.873}{0.05 s^4 + 1.08 s^3 + 1.71 s^2 + 2.094 s}$$

Testing và validation

Nhóm nghiên cứu thiết lập bộ điều khiển hiệu chỉnh tỉ lệ - tích phân - vi phân tích hợp khâu bù sớm pha (Lead-Lag PID Controller) có hàm truyền: $$R_\vartheta(s) = K_p \left( 1 + \frac{1}{T_i s} + \frac{T_d s}{\tau_d s + 1} \right) \cdot \frac{T_1 s + 1}{T_2 s + 1}$$

Sau khi tối ưu hóa tham số bằng quỹ đạo nghiệm số trên MATLAB: $K_p = 1.85, ; T_i = 1.25\text{ s}, ; T_d = 0.18\text{ s}, ; \tau_d = 0.018\text{ s}, ; T_1 = 0.12\text{ s}, ; T_2 = 0.03\text{ s}$.

% Tong hop bo dieu khien va danh gia dap ung vong kin
R_theta = tf([1.85*0.12, 1.85*(0.12+1), 1.85], [0.03, 1.03, 1]) * tf([0.22, 1], [0.02, 1]);

% Ham truyen kin he thong doi voi tin hieu dat
W_kin = feedback(R_theta * W_cd * W_plant, 1);

% Ham truyen kin he thong doi voi nhieu gio dung Wy
W_kin_wind = feedback(W_wind, R_theta * W_cd * W_plant);

% Kiem tra chi tieu chat luong buoc nhay don vi Step Response
figure;
step(W_kin, 5);
grid on;
title('Dap ung qua do goc chuc ngoc \theta(t) voi tin hieu dat Step');
stepinfo(W_kin)

Kết quả kiểm chứng thực nghiệm mô phỏng qua bảng số liệu sau:

Thông số động học / Chỉ tiêu Hệ hở (Chưa điều khiển) Hệ kín (Bộ điều khiển PID kinh điển) Hệ kín tối ưu Lead-Lag PID (Đồ án)
Độ quá điều chỉnh ($\sigma%$) $48.6%$ (Dao động kéo dài) $18.2%$ $6.4%$
Thời gian quá độ ($t_{qd}$, sai số $2%$) $8.45\text{ s}$ $3.10\text{ s}$ $1.35\text{ s}$
Sai số xác lập ($e_{xl}$) $\infty$ (Tích phân trôi dạt) $0.02^\circ$ $0.00^\circ$
Độ dự trữ biên ($Gain Margin - GM$) $3.2\text{ dB}$ $7.8\text{ dB}$ $12.4\text{ dB}$
Độ dự trữ pha ($Phase Margin - PM$) $28.4^\circ$ $46.5^\circ$ $58.2^\circ$
Biên độ lắc góc cực đại khi $W_y = 10\text{ m/s}$ $\Delta \vartheta_{max} = 7.85^\circ$ $\Delta \vartheta_{max} = 2.40^\circ$ $\Delta \vartheta_{max} = 0.65^\circ$

Kết quả đạt được

TIẾN ĐỘ HOÀN THÀNH MỤC TIÊU ĐỒ ÁN
1. Xây dựng mô hình toán 6-DOF:       [████████████████████] 100%
2. Tuyến tính hóa & hàm truyền:       [████████████████████] 100%
3. Tổng hợp bộ điều khiển Lead-Lag:   [████████████████████] 100%
4. Bù trễ cơ cấu servo:              [████████████████████] 100%
5. Mô phỏng triệt tiêu nhiễu gió:     [████████████████████] 100%

Hệ thống điều khiển tổng hợp đã triệt tiêu hoàn toàn dao động góc chúc ngóc khi có luồng gió giật đứng tác động. Tốc độ triệt tiêu nhiễu đạt $1.35\text{ s}$, đảm bảo quỹ đạo bay thẳng bằng phẳng, bảo vệ an toàn cho khí cụ bay trong điều kiện thời tiết gió cấp 4–5 ($W_y \le 10\text{ m/s}$).

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

  1. Mô hình hóa chính xác hàm truyền đa biến có xét đến trễ servo: Thay vì coi cơ cấu chấp hành có đáp ứng tức thời ($W_{cd}(s) = 1$) như nhiều đồ án lý thuyết thuần túy, đề tài đưa hàm truyền bậc 1 của động cơ servo ($T_{cd} = 0.05\text{ s}$) vào vòng lặp thiết kế, loại bỏ hoàn toàn hiện tượng tự kích dao động thực tế.
  2. Kỹ thuật bù sớm trễ pha kết hợp lọc bậc cao: Cấu trúc bộ điều chỉnh được tối ưu hóa giúp tăng độ dự trữ pha từ $28.4^\circ$ lên $58.2^\circ$ (cải thiện $104.9%$), cho phép hệ thống duy trì tính ổn định bền vững ngay cả khi các hệ số khí động sai lệch $\pm 20%$ do thay đổi áp suất khí quyển theo độ cao.
  3. Triệt tiêu ảnh hưởng của thành phần gió đứng ($W_y$): Giảm biên độ lắc chúc ngóc cực đại từ $7.85^\circ$ xuống $0.65^\circ$ (giảm $91.7%$ độ nhạy nhiễu), giúp thiết bị duy trì trục camera quang điện tử ổn định khi thực hiện nhiệm vụ chụp ảnh viễn thám hoặc trinh sát mục tiêu mặt đất.

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

Kịch bản ứng dụng thực tế

|                  CÁC LĨNH VỰC TRIỂN KHAI ỨNG DỤNG                       |
| 1. QUÂN SỰ - PHÒNG KHÔNG:                                               |
|    - Mục tiêu bay huấn luyện xạ kích (M-96A, M-100CT, M-400CT).         |
|    - Thiết bị bay trinh sát chiến thuật, chỉ thị mục tiêu thời gian thực|
|                                                                         |
| 2. KINH TẾ - DÂN SỰ:                                                    |
|    - Bay quét bản đồ 3D địa hình (Photogrammetry) bằng cảm biến LiDAR.  |
|    - Giám sát hành lang đường dây tải điện 500kV, tuần tra rừng.        |
|    - Cứu hộ cứu nạn hàng hải và kiểm soát ô nhiễm tràn dầu ven biển.    |

Lộ trình triển khai kỹ thuật

Giai đoạn 1: Mô phỏng HIL (Hardware-In-the-Loop) trên KIT STM32F4 + MATLAB
Giai đoạn 2: Lắp đặt tích hợp trên khung vỏ Composite UAV mô hình sải cánh 2.4m
Giai đoạn 3: Thử nghiệm hầm gió khí động và bay thử nghiệm thực địa bãi thử
     (Thời gian: Tháng 7 - Tháng 9)

Phân tích hiệu quả kinh tế (ROI):

  • Việc làm chủ thuật toán điều khiển và tự chế tạo mạch Auto-Pilot trong nước giúp hạ giá thành sản xuất mục tiêu bay huấn luyện từ $45,000\text{ USD}$ (nhập ngoại) xuống dưới $8,000\text{ USD}$ (nội địa hóa), tiết kiệm hơn $80%$ ngân sách quốc phòng cho công tác huấn luyện pháo - tên lửa phòng không hàng năm.

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

Mặc dù đạt được các chỉ tiêu động học đề ra, đề tài vẫn tồn tại một số điểm hạn chế kỹ thuật:

  • Tuyến tính hóa cục bộ: Bộ điều khiển được thiết kế xung quanh điểm làm việc danh định ($V = 200\text{ m/s}, H = 1000\text{ m}$). Khi UAV thay đổi tốc độ đột ngột trong dải rộng ($100 - 350\text{ m/s}$), chất lượng điều khiển có thể bị suy giảm.
  • Chưa bù trôi dạt trọng tâm: Chưa tích hợp thuật toán bù thích nghi tham số khi khối lượng nhiên liệu giảm dần trong suốt hành trình bay dài.

Hướng phát triển tiếp theo:

  1. Nghiên cứu áp dụng các cấu trúc điều khiển phi tuyến hiện đại: Điều khiển trượt (Sliding Mode Control - SMC), Điều khiển thích nghi mờ - nơ-ron (Adaptive Neuro-Fuzzy Inference System - ANFIS) hoặc Điều khiển $\mathcal{H}_\infty$ bền vững.
  2. Tích hợp trọn vẹn 3 kênh điều khiển không gian: Kênh dọc (Pitch), Kênh ngang (Roll - Cánh liệng $\delta_l$) và Kênh hướng (Yaw - Bánh lái hướng $\delta_h$) thành hệ thống điều hướng tự động tích hợp GPS/INS hoàn chỉnh.

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

  • Sinh viên ngành Điện - Tự động hóa, Kỹ thuật Hàng không: Nắm bắt được phương pháp luận mô hình hóa giải tích cơ học bay 6-DOF phức tạp và kỹ thuật thiết kế bộ điều khiển trên MATLAB/Simulink.
  • Kỹ sư nhúng và điều khiển bay (Flight Control Engineers): Sở hữu cấu trúc hàm truyền đã chuẩn hóa và các thông số thực nghiệm để phát triển firmware Auto-Pilot cho các dòng UAV cánh bằng (Fixed-wing UAV).
  • Các viện nghiên cứu và doanh nghiệp sản xuất quốc phòng: Tiết kiệm thời gian R&D trong việc thiết kế hệ thống điều khiển tự động cho các loại mục tiêu bay huấn luyện quân sự.

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

1. Tại sao đồ án sử dụng số liệu khí động của máy bay MiG-21Bis thay vì một mô hình UAV cụ thể?

Trong nghiên cứu động lực học bay, MiG-21Bis là dòng khí cụ bay cánh tam giác kinh điển có đầy đủ bộ thông số vi phân khí động học thực nghiệm chính xác qua hầm gió và bay thử. Việc sử dụng bộ thông số này giúp kiểm chứng thuật toán điều khiển tổng hợp với độ tin cậy tuyệt đối, làm mô hình chuẩn trước khi áp dụng vào các khung thân UAV thực tế có cấu hình tương tự.

2. Sự khác biệt cơ bản giữa hệ tọa độ tốc độ thẳng đứng ($Oxyz$) và hệ tọa độ liên kết ($Ox_1y_1z_1$) là gì?

Hệ tọa độ liên kết $Ox_1y_1z_1$ gắn cố định với thân máy bay (trục $Ox_1$ trùng trục dọc thân), dùng để tính toán các mômen quán tính ($J_x, J_y, J_z$) và chuyển động quay quanh trọng tâm. Hệ tọa độ tốc độ thẳng đứng $Oxyz$ có trục $Ox$ hướng theo vectơ vận tốc $\mathbf{V}$, trục $Oy$ nằm trong mặt phẳng thẳng đứng chứa trục dọc, thuận tiện nhất để thiết lập phương trình chuyển động của trọng tâm và khảo sát góc nghiêng quỹ đạo $\theta$.

3. Nhiễu gió đứng ($W_y$) tác động trực tiếp vào tham số khí động nào của UAV?

Khi có gió đứng $W_y$, góc tấn tức thời của máy bay thay đổi một lượng $\Delta \alpha \approx \frac{W_y}{V}$. Biến thiên góc tấn này trực tiếp làm thay đổi lực nâng $\Delta Y = C_y^\alpha \Delta \alpha \frac{\rho V^2}{2} S$ và tạo ra mômen ngóc $\Delta M_z = m_z^\alpha \Delta \alpha \frac{\rho V^2}{2} S b_A$, gây mất cân bằng mômen dọc của máy bay.

4. Cơ cấu trợ dẫn servo ảnh hưởng thế nào đến tính ổn định của hệ thống điều khiển?

Servo là khâu trễ quán tính (thường có hàm truyền $W_{cd}(s) = \frac{k_{cd}}{T_{cd}s + 1}$). Nếu hằng số thời gian $T_{cd}$ lớn, servo sẽ gây trễ pha lớn cho tín hiệu điều khiển góc bẻ lái $\delta_B$. Nếu không được bù trễ pha bằng bộ điều khiển (như khâu sớm pha Lead), hệ thống rất dễ rơi vào vùng mất ổn định và xuất hiện tự dao động biên độ lớn.

5. Chi phí phần cứng để triển khai hệ thống điều khiển ổn định góc chúc ngóc này là bao nhiêu?

Để triển khai thuật toán lên mô hình thực địa, chỉ cần bo mạch vi điều khiển ARM Cortex-M4 (STM32F405/F427 giá khoảng $15 - 25\text{ USD}$), module cảm biến IMU 6 trục MPU-6050/ICM-20602 (khoảng $5 - 10\text{ USD}$) và các servo kỹ thuật số kim loại tốc độ cao ($15 - 30\text{ USD}$/chiếc). Tổng chi phí phần cứng điện tử điều khiển dưới $100\text{ USD}$, hoàn toàn khả thi cho sản xuất hàng loạt.

Kết luận

Đồ án tốt nghiệp "Nghiên cứu thiết bị bay không người lái, chế tạo mô hình có điều khiển từ xa" đã giải quyết trọn vẹn bài toán cốt lõi của cơ học bay: thiết lập hệ phương trình chuyển động không gian 6 bậc tự do, phân tích động lực học chuyển động dọc và tổng hợp thành công hệ thống điều khiển tự động ổn định góc chúc ngóc khi chịu tác động của nhiễu gió đứng tầng thấp.

Các đóng góp khoa học và thực tiễn nổi bật:

  • Xây dựng tường minh hệ phương trình vi phân chuyển động dọc tuyến tính hóa kết hợp mô hình động học servo.
  • Tổng hợp bộ điều khiển tối ưu cải thiện vượt bậc các chỉ tiêu chất lượng: thời gian quá độ $t_{qd} = 1.35\text{ s}$, độ quá điều chỉnh $\sigma% = 6.4%$, triệt tiêu $91.7%$ biên độ dao động do gió giật thẳng đứng.
  • Đặt nền móng thuật toán vững chắc để nâng cấp các dòng mục tiêu bay không người lái của Việt Nam từ điều khiển thủ công sang bay tự động hành trình chính xác cao, phục vụ hiệu quả cho công tác huấn luyện an ninh quốc phòng và phát triển kinh tế đất nước.