Tổng quan về luận án

Sự phát triển mạnh mẽ của phương tiện bay không người lái (Unmanned Aerial Vehicle - UAV) trong tác chiến hiện đại và ứng dụng dân sự đã đặt ra những yêu cầu khắt khe về tính tự chủ, độ tin cậy và mức độ an toàn vận hành. Trong toàn bộ chu trình thực hiện nhiệm vụ, giai đoạn tiếp cận và hạ cánh là quá trình động học phức tạp nhất, tiềm ẩn nguy cơ tai nạn cao nhất do vận tốc bay giảm dần, hiệu ứng mặt đất gia tăng, và ảnh hưởng mạnh mẽ của nhiễu động khí quyển tầng thấp. 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ố: 9.16 / 9.52.03.01) của nghiên cứu sinh Ngô Văn Toàn, dưới sự hướng dẫn khoa học của Đại tá, TS. Lê Thanh Phong và Đại tá, TS. Nguyễn Xuân Căn tại Học viện Kỹ thuật Quân sự, mang tên "Tổng hợp thuật toán điều khiển hạ cánh theo chương trình cho máy bay không người lái cỡ nhỏ", đã giải quyết trực tiếp bài toán then chốt này.

Khoảng trống nghiên cứu (research gap) cốt lõi trong y văn hiện đại nằm ở sự thiếu hụt các thuật toán tối ưu hóa quỹ đạo hạ cánh phi tuyến có tính đến giới hạn quá tải kết cấu và khả năng kháng nhiễu gió thời gian thực cho UAV cỡ nhỏ cánh cố định. Trong khi các giải pháp thương mại và nghiên cứu trước đây chủ yếu dựa trên quỹ đạo hình học giản đơn (các cung tròn liên tiếp như nghiên cứu MiG-21 [6]), phụ thuộc vào hạ cánh vô tuyến (đài chuẩn hướng $\text{KРМ}$, đài chuẩn tầm $\text{ГРМ}$ [4, 5]), hoặc dẫn đường quang học bằng camera [31, 41] vốn rất nhạy cảm với điều kiện sương mù và mưa bão, việc xây dựng quỹ đạo tối ưu động lực học khép kín theo chương trình có xét đến hạn chế quá tải đứng chưa được giải quyết trọn vẹn. Luận án xác lập hệ thống câu hỏi nghiên cứu và giả thuyết khoa học:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để giải bài toán biên hai điểm phi tuyến tính cao (TPBVP) nhằm xác lập quỹ đạo hạ cánh tối ưu mượt mà cho UAV cỡ nhỏ dưới các ràng buộc nghiêm ngặt về khí động và cơ học tiếp đất?
  • Câu hỏi nghiên cứu 2 (RQ2): Cấu trúc mạch vòng điều khiển chuyển động dọc nào có khả năng bám chính xác quỹ đạo tối ưu danh định đồng thời triệt tiêu độ lệch sinh ra bởi nhiễu động gió dọc và gió đứng?
  • Giả thuyết khoa học 1 (H1): Ứng dụng Nguyên lý cực đại Pontryagin kết hợp phương pháp tham số liên tục (Homotopy Continuation Method) sẽ giải quyết triệt để sự phân kỳ số của bài toán biên, tìm ra quỹ đạo tối ưu toàn cục thỏa mãn điều kiện biên tiếp đất.
  • Giả thuyết khoa học 2 (H2): Cấu trúc điều khiển kết hợp phản hồi tín hiệu quá tải đứng $n_y$ và góc chúc ngóc $\vartheta$ sẽ duy trì góc tấn không tốc $\alpha_r$ dưới ngưỡng thất tốc, triệt tiêu hiện tượng dập đuôi hoặc lún quá tốc khi tiếp đất.

Khung lý thuyết của công trình tích hợp lý thuyết động lực học bay phi tuyến 6 bậc tự do, lý thuyết điều khiển tối ưu biến phân Pontryagin và lý thuyết điều khiển tự động tuyến tính hóa. Đối tượng khảo nghiệm cụ thể là mẫu UAV-70V cánh cố định (khối lượng cất cánh $m = 56.5\text{ kg}$, diện tích cánh hiệu dụng $S = 1.05\text{ m}^2$). Phạm vi khảo sát từ độ cao tiếp cận ban đầu $H_0 = 60\text{ m}$, vận tốc ban đầu $V_0 = 50\text{ m/s}$ tại cự ly $x_0 = 0\text{ m}$, góc nghiêng quỹ đạo $\theta_0 = 0\text{ rad}$, cho đến điểm tiếp đất mong muốn $x_f = 500\text{ m}$, độ cao $y_f = 0.7\text{ m}$ (khoảng cách trọng tâm đến bánh càng sau khi ngóc mũi $12^\circ$), vận tốc tiếp đất danh định $V_f = 31\text{ m/s}$, góc nghiêng quỹ đạo tiếp đất $\theta_f = 0\text{ rad}$, vận tốc lún $V_{yhc} \le 1.0\text{ m/s}$, và sai số tọa độ tiếp đất cực kỳ chặt chẽ: sai số độ cao $0\text{ m} \le \Delta y \le 0.3\text{ m}$ cùng sai số cự ly dọc $\Delta x \le 30\text{ m}$.

Literature Review và Positioning

Khảo cứu y văn quốc tế và trong nước chỉ ra ba luồng nghiên cứu chủ đạo trong lĩnh vực điều khiển hạ cánh thiết bị bay không người lái:

Luồng nghiên cứu thứ nhất tập trung vào hệ thống thông tin dẫn đường quang học và thị giác máy tính. Điển hình là các công trình của Barber et al. [31], Cesetti et al. [41, 43] và Yakimenko et al. [50], sử dụng camera trên khoang để nhận dạng vạch biên đường băng [59], ước lượng tư thế 3D theo thời gian thực và cung cấp tín hiệu điều khiển phản hồi vị trí. Nhược điểm chí tử của luồng tiếp cận này là sự phụ thuộc tuyệt đối vào tầm nhìn khí tượng; khi gặp thời tiết phức tạp (mây mù, ban đêm không đèn chiếu, mưa lớn), thuật toán xử lý ảnh lập tức mất dấu mục tiêu.

Luồng nghiên cứu thứ hai khai thác các giải pháp điều khiển thông minh và thích nghi phi tuyến. Các công trình ứng dụng logic mờ (Fuzzy Logic) của Santos et al. [17], suy diễn mờ kết hợp mạng nơ-ron [22, 23, 33], điều khiển trượt thích nghi phi tuyến [13, 16, 19], thuật toán di truyền cải tiến (GA) kết hợp giải thuật tìm đường A* [40], và lý thuyết ổn định Lyapunov [34]. Mặc dù các bộ điều khiển thích nghi có khả năng bù trừ bất định mô hình, phần lớn các công trình này chỉ giải quyết bài toán bám các điểm mốc (waypoints) được nối bằng đoạn thẳng hoặc cung tròn tĩnh, chưa tối ưu hóa chỉ tiêu năng lượng và quá tải cơ học trong toàn bộ quá trình hạ độ cao, kéo bằng ($L_{kb}$), giữ bằng ($L_{gb}$) và tiếp đất.

Luồng nghiên cứu thứ ba tập trung vào bài toán hạ cánh chuyên biệt như hạ cánh trên tàu sân bay dưới dao động sóng biển của Yang et al. [24, 44], sử dụng khai triển chuỗi Fourier để mô tả động lực học chuyển động sàn tàu. Tuy nhiên, các mô hình này đòi hỏi hệ thống cảm biến phức tạp và công suất tính toán lớn, khó tích hợp trên máy tính nhúng kích thước nhỏ của UAV hạng nhẹ.

Trong bức tranh học thuật đó, tồn tại một cuộc tranh luận sâu sắc giữa hai trường phái:

  • Trường phái tiếp cận hình học/kinematic: Ưu tiên các quỹ đạo đa thức bậc cao hoặc đường cong Dubins nhằm giảm thiểu tải tính toán cho máy tính trên khoang, nhưng chấp nhận việc quá tải đứng $n_y$ và góc tấn $\alpha$ biến thiên tự do, dẫn đến nguy cơ vượt ngưỡng bền kết cấu khi gặp gió chướng.
  • Trường phái tiếp cận tối ưu hóa động lực học (Dynamic Optimal Control): Đòi hỏi giải hệ phương trình vi phân phi tuyến đầy đủ gắn liền với hàm mục tiêu tối ưu. Trở ngại lớn nhất của trường phái này là bài toán biên hai điểm (TPBVP) cực kỳ nhạy cảm với điều kiện ban đầu của các biến đồng trạng thái ($P_V, P_\theta, P_x, P_y$), dẫn đến hiện tượng phân kỳ toán học khi giải bằng các phương pháp bắn (shooting methods) thông thường.

Luận án của Ngô Văn Toàn đã định vị chính xác khoảng trống này bằng cách kết hợp Nguyên lý cực đại Pontryagin với phương pháp giải liên tục theo tham số $s \in [0, 1]$. So sánh với nghiên cứu hạ cánh bám ảnh của Barber [31] và điều khiển mờ của Santos [17], nghiên cứu này vượt trội ở khả năng áp đặt trực tiếp các ràng buộc bất đẳng thức vật lý: quá tải đứng $1.0 \le n_y \le 3.5$, góc chúc ngóc tiếp đất $0^\circ \le \vartheta \le 12^\circ$ (tránh va quệt đuôi), đồng thời thiết lập cơ chế bù góc tấn cảm ứng do gió $\alpha_w \approx W_y/V_k$ trong kênh chuyển động dọ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 lý thuyết điều khiển tối ưu Pontryagin cho bài toán cơ học bay phi tuyến của UAV cánh cố định thông qua việc thiết lập hệ phương trình trạng thái 4 biến trong hệ tọa độ quỹ đạo ($Ox_k y_k z_k$):

  • Phương trình vi phân vận tốc: $\frac{dV_k}{dt} = g(n_x - \sin\theta)$
  • Phương trình vi phân góc nghiêng quỹ đạo: $\frac{d\theta}{dt} = \frac{g}{V_k}(n_y - \cos\theta)$
  • Phương trình tọa độ dọc: $\frac{dx_g}{dt} = V_k \cos\theta$
  • Phương trình độ cao: $\frac{dy_g}{dt} = V_k \sin\theta$

Hàm Hamilton của hệ thống được tổng hợp dưới dạng: $$H = P_V g(n_x - \sin\theta) + P_\theta \frac{g}{V_k}(n_y - \cos\theta) + P_x V_k \cos\theta + P_y V_k \sin\theta - L(u, t)$$

Trong đó, $P_V, P_\theta, P_x, P_y$ là các biến đồng trạng thái (co-state variables) thỏa mãn hệ phương trình vi phân liên hợp Euler-Lagrange: $$\dot{P}V = -\frac{\partial H}{\partial V_k}; \quad \dot{P}\theta = -\frac{\partial H}{\partial \theta}; \quad \dot{P}_x = -\frac{\partial H}{\partial x_g} = 0; \quad \dot{P}_y = -\frac{\partial H}{\partial y_g} = 0$$

Đóng góp lý thuyết mang tính đột phá là việc xây dựng thuật toán nhúng tham số biến dạng liên tục $s$ để chuyển hóa bài toán phi tuyến phức tạp từ một nghiệm giải tích đã biết tại $s=0$ dần về bài toán phi tuyến thực tế tại $s=1$. Giải pháp này triệt tiêu hoàn toàn sự phân kỳ nghiệm số của bài toán biên hai điểm, chứng minh bằng toán học sự tồn tại và duy nhất của quỹ đạo tối ưu hạ cánh có ràng buộc quá tải đứng.

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

Khung phân tích của luận án tích hợp chặt chẽ ba trụ cột lý thuyết: (1) Lý thuyết động lực học và khí động lực học công trình; (2) Lý thuyết điều khiển tối ưu giải tích Pontryagin; và (3) Lý thuyết điều khiển tuyến tính đa biến cho hệ phi tuyến cục bộ.

Mối quan hệ không gian giữa các hệ tọa độ được chuẩn hóa thông qua các ma trận quay Euler: ma trận chuyển từ hệ tọa độ chuẩn ($Ox_g y_g z_g$) sang hệ liên kết ($Ox_1 y_1 z_1$) qua các góc chúc ngóc $\vartheta$, góc nghiêng $\gamma$, góc hướng $\psi$; ma trận chuyển sang hệ tốc độ ($Ox_r y_r z_r$) qua góc tấn không tốc $\alpha_r$ và góc trượt không tốc $\beta_r$; và ma trận chuyển sang hệ quỹ đạo ($Ox_k y_k z_k$) qua góc nghiêng quỹ đạo $\theta$ và góc hướng quỹ đạo $\Psi$.

Điều kiện biên xác định cho bài toán bao gồm:

  • Biên ban đầu ($t=t_0$): $V(0) = 50\text{ m/s}, \theta(0) = 0\text{ rad}, x(0) = 0\text{ m}, y(0) = 60\text{ m}$.
  • Biên kết thúc ($t=t_f$): $V(t_f) = 31\text{ m/s}, \theta(t_f) = 0\text{ rad}, x(t_f) = 500\text{ m}, y(t_f) = 0.7\text{ m}$.
  • Vận tốc tiếp đất danh định được rút ra từ phương trình cân bằng trọng lực và lực nâng tại thời điểm chạm bánh: $$V_{hc} = \sqrt{\frac{2G}{\rho C_{yHC} S}} = \sqrt{\frac{2 \times 56.5 \times 9.81}{1.225 \times 1.05 \times C_{yHC}}} \approx 31\text{ m/s}$$ với gia tốc trọng trường $g = 9.80665\text{ m/s}^2$ và mật độ không khí tầng sát đất $\rho = 1.225\text{ kg/m}^3$.

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

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên lập trường nhận thức luận duy thực hậu thực chứng (Post-positivism / Critical Realism), kết hợp suy diễn toán học hình thức với mô phỏng số kiểm chứng trên phần mềm kỹ thuật tiêu chuẩn MATLAB 2015 / Simulink. Thiết kế nghiên cứu phân tầng hai mức:

  1. Tầng tối ưu hóa cấp cao (High-level trajectory optimization): Sử dụng giải thuật số giải bài toán cực đại Pontryagin để trích xuất các hàm quỹ đạo chương trình thời gian thực: độ cao chuẩn $H_{ct}(t)$, góc chúc ngóc chuẩn $\vartheta_{ct}(t)$, quá tải đứng chuẩn $n_{yct}(t)$, và vận tốc chuẩn $V_{ct}(t)$.
  2. Tầng điều khiển bám cấp thấp (Low-level tracking control): Tổng hợp bộ điều khiển phản hồi trạng thái kênh chuyển động dọc để điều khiển góc lệch cánh lái độ cao $\delta_c$ và tay ga động cơ $\delta_g$, kháng lại tác động của vector gió $\vec{W}$.

Quy trình nghiên cứu rigorous

Quy trình mô hình hóa động học kênh dọc bắt đầu từ hệ phương trình vi phân chuyển động quanh trọng tâm: $$J_z \frac{d\omega_z}{dt} = M_z(\alpha, \omega_z, \delta_c, \dot{\alpha}, V_r)$$ trong đó $J_z$ là mô men quán tính quanh trục ngang $Oz_1$, $\omega_z$ là tốc độ góc chúc ngóc, $M_z$ là mô men khí động học được tuyến tính hóa: $$M_z = m_z^\alpha \alpha + m_z^{\omega_z} \bar{\omega}_z + m_z^{\delta_c} \delta_c + m_z^{\dot{\alpha}} \dot{\bar{\alpha}}$$ với $\bar{\omega}_z = \frac{\omega_z b_a}{V_r}$ và $b_a$ là dây cung khí động trung bình.

Hàm truyền đạt kênh chuyển động dọc từ góc lệch cánh lái độ cao $\delta_c(p)$ đến góc chúc ngóc $\vartheta(p)$ và quá tải đứng $n_y(p)$ được thiết lập qua phép biến đổi Laplace: $$W_\vartheta^{\delta_c}(p) = \frac{k_\vartheta (T_\vartheta p + 1)}{p(T_1^2 p^2 + 2\xi T_1 p + 1)}; \quad W_{n_y}^{\delta_c}(p) = \frac{k_{n_y} (T_{n1}^2 p^2 + 2\xi_n T_{n1} p + 1)}{T_1^2 p^2 + 2\xi T_1 p + 1}$$

Để khử dao động và nâng cao độ ổn định cho mạch vòng bám, nghiên cứu sinh tích hợp hệ số cản dịu kênh dọc $K_z$ vào vòng phản hồi tốc độ góc $\omega_z$, đồng thời thiết kế mạch vòng điều khiển kết hợp: sử dụng tín hiệu sai lệch độ cao $\Delta H = H_{ct} - H$ qua khâu hiệu chỉnh PID tạo ra tín hiệu quá tải yêu cầu $n_{yct}$, từ đó phối hợp với phản hồi góc chúc ngóc $\vartheta$ để sinh lệnh điều khiển cánh lái $\delta_c$.

Mô hình gió tầng sát đất được mô hình hóa theo quy luật hàm mũ thực nghiệm: $$W_0(y_0) = W_{01} \left(\frac{y_0}{y_{01}}\right)^n$$ với số mũ khí quyển $n \in [0.1, 0.4]$, độ cao tham chiếu $y_{01} = 10\text{ m}$, vận tốc gió trung bình $W_{01} = 3 \div 4\text{ m/s}$. Vector gió tổng quát $\vec{W} = [W_x, W_y, W_z]^T$ được phân rã thành thành phần gió dọc $W_x$ làm biến đổi trực tiếp không tốc $V_r = V_k \pm W_x$ và gió đứng $W_y$ làm biến đổi góc tấn không tốc một lượng $\Delta \alpha_w \approx \frac{W_y}{V_k}$.

Data và phân tích

Các tham số hình học và khí động của UAV-70V được số hóa chính xác: trọng lượng cất cánh $G = 56.5 \times 9.80665\text{ N} \approx 554.1\text{ N}$, diện tích cánh $S = 1.05\text{ m}^2$, giới hạn quá tải kết cấu $1.0 \le n_y \le 3.5$.

Kiểm tra tính vững (Robustness checks) được thực hiện qua hàng loạt kịch bản mô phỏng số:

  1. Kịch bản cơ sở: Hạ cánh trong môi trường tĩnh không gió ($W_x = 0, W_y = 0$).
  2. Kịch bản gió dọc: Khảo sát gió ngược $W_x = +5\text{ m/s}$ và gió xuôi $W_x = -5\text{ m/s}$.
  3. Kịch bản gió đứng: Khảo sát luồng khí thổi thẳng đứng từ dưới lên $W_y = +2\text{ m/s}$ và luồng khí thổi từ trên xuống $W_y = -2\text{ m/s}$.
  4. Kịch bản biến thiên điều kiện biên: Thay đổi độ cao ban đầu $H_0 \in [50\text{ m}, 80\text{ m}]$ và vận tốc ban đầu $V_0 \in [45\text{ m/s}, 55\text{ m/s}]$.

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

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

Kết quả mô phỏng số trên môi trường MATLAB/Simulink chứng minh hiệu năng vượt trội của hệ thống điều khiển được tổng hợp:

  1. Tối ưu hóa quá tải và góc chúc ngóc: Quỹ đạo hạ cánh tối ưu Pontryagin kiểm soát quá tải đứng cực đại ở mức $n_{y\max} = 1.82$, nằm hoàn toàn trong giới hạn an toàn kết cấu ($1.0 \le n_y \le 3.5$), và đưa quá tải tiếp đất về chính xác giá trị $n_{yhc} = 1.002 \approx 1.0$. Góc chúc ngóc tại thời điểm tiếp đất đạt $\vartheta(t_f) = 6.8^\circ$, triệt tiêu hoàn toàn khả năng chạm đuôi (ngưỡng giới hạn $\vartheta \le 12^\circ$).
  2. Phát hiện về hiệu ứng gió đứng ($W_y$): Khi xuất hiện gió đứng $W_y = +2\text{ m/s}$, góc tấn không tốc tức thời tăng thêm $\Delta \alpha_w \approx 2.3^\circ$. Nếu chỉ sử dụng bộ điều khiển ổn định góc chúc ngóc thông thường, lực nâng tăng đột ngột làm UAV bị bốc cao ("nhảy cóc" - ballooning), kéo dài điểm tiếp đất ra ngoài đường băng ($\Delta x > 65\text{ m}$). Ngược lại, bộ điều khiển kết hợp phản hồi quá tải $n_y$ đã phát hiện tức thì sự gia tăng gia tốc pháp tuyến, nhanh chóng điều khiển cánh lái độ cao $\delta_c$ chúc mũi xuống, giữ sai số độ cao tiếp đất trong mức $\Delta y = 0.12\text{ m}$ (nhỏ hơn nhiều so với tiêu chuẩn cho phép $0.3\text{ m}$).
  3. Độ êm thuận khi tiếp đất: Vận tốc lún thẳng đứng tại thời điểm chạm bánh đạt $V_{yhc} = 0.38\text{ m/s}$ khi không gió, và đạt $V_{yhc} = 0.62\text{ m/s}$ khi có gió đứng cực đại, hoàn toàn đáp ứng chỉ tiêu kỹ thuật khắt khe $V_{yhc} \le 1.0\text{ m/s}$ để bảo vệ hệ thống giảm chấn càng máy bay.
  4. Bám cự ly hạ cánh: Sai số vị trí tiếp đất dọc theo trục đường băng được kiểm soát chặt chẽ: $\Delta x = 12.4\text{ m}$ trong trường hợp gió dọc $W_x = \pm 5\text{ m/s}$, thỏa mãn hoàn toàn chỉ tiêu $\Delta x \le 30\text{ m}$.

Implications đa chiều

  • Về mặt lý thuyết: Khẳng định tính hiệu quả của phương pháp tham số liên tục trong việc giải bài toán điều khiển tối ưu biến phân không gian trạng thái khí động, cung cấp một khung giải thuật mẫu để xử lý các bài toán điều khiển phi tuyến có ràng buộc bất đẳng thức.
  • Về mặt phương pháp luận: Thiết lập phương pháp tích hợp hai tầng (two-tier integration) giữa tối ưu hóa quỹ đạo off-line/on-line và điều khiển bám phản hồi kín đa vòng có cản dịu, có khả năng chuyển giao áp dụng cho nhiều lớp phương tiện bay khác nhau.
  • Về mặt thực tiễn và chính sách: Cung cấp giải pháp công nghệ tự chủ, giảm thiểu sự phụ thuộc vào các module điều khiển thương mại nước ngoài (như Pixhawk/ArduPilot vốn chỉ bám điểm đơn giản). Kết quả này cho phép các đơn vị công nghiệp quốc phòng như Tập đoàn Viettel, Viện Kỹ thuật Phòng không - Không quân nạp trực tiếp các bảng dữ liệu quỹ đạo tối ưu (Lookup Tables) vào máy tính điều khiển bay (Flight Control Computer - FCC), nâng cao độ tin cậy khi UAV thực thi nhiệm vụ tại các vùng biển đảo có địa hình hiểm trở và gió giật phức tạp.

Limitations và Future Research

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

  1. Giả định phân tách kênh độc lập: Nghiên cứu chỉ tập trung chuyên sâu vào kênh chuyển động dọc, giả định rằng kênh chuyển động cạnh (góc liệng $\gamma$, góc trượt $\beta$, góc hướng $\psi$) đã được giữ ổn định độc lập bởi các mạch vòng riêng biệt. Trong thực tế hạ cánh khi có gió ngang mạnh ($W_z$), hiện tượng ghép kênh khí động phi tuyến (cross-coupling) giữa chuyển động dọc và chuyển động cạnh có thể làm suy giảm chất lượng bám quỹ đạo.
  2. Mô hình gió dừng: Luận án mới khảo sát thành phần gió ổn định và biến thiên theo độ cao theo hàm mũ, chưa tích hợp mô hình nhiễu loạn gió ngẫu nhiên dạng phổ Dryden hoặc von Kármán (Turbulent Wind Gusts).
  3. Thử nghiệm phần cứng: Các kết quả đánh giá mới dừng lại ở mức mô phỏng số chuyên sâu trên phần mềm MATLAB/Simulink, chưa tiến hành thử nghiệm mô phỏng bán tự nhiên (Hardware-in-the-Loop - HIL) hoặc bay thực địa trên bãi thử.

Chương trình nghiên cứu tiếp theo cần tập trung vào:

  • Mở rộng bài toán tối ưu quỹ đạo 6 bậc tự do (6-DOF) ghép kênh không gian 3 chiều đầy đủ, đồng thời hạn chế góc nghiêng cánh liệng $\gamma$ khi tiếp đất.
  • Ứng dụng các thuật toán điều khiển bền vững thích nghi phi nonlinear $H_\infty$ hoặc mạng nơ-ron sâu (Deep Reinforcement Learning) để bù trừ trễ chấp hành của servo cánh lái trong điều kiện gió giật mạnh.
  • Tiến hành thử nghiệm nhúng thuật toán trên vi xử lý ARM Cortex-M7/FPGA và thực hiện các chuyến bay thử nghiệm thực địa trên dòng UAV-70V.

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

Luận án tạo ra tác động sâu rộng trên nhiều bình diện:

  • Tác động học thuật: Đóng góp 02 bài báo khoa học chuyên ngành chất lượng cao trên các tạp chí uy tín thuộc danh mục Scopus/Văn mục Hội đồng Giáo sư Nhà nước; thiết lập nền tảng lý thuyết điều khiển tối ưu hạ cánh cho các chương trình đào tạo sau đại học ngành Tự động hóa và Kỹ thuật Hàng không tại Học viện Kỹ thuật Quân sự.
  • Chuyển đổi công nghiệp quốc phòng: Cung cấp thuật toán cốt lõi cho các dự án phát triển UAV trinh sát chiến dịch - chiến thuật của Quân đội nhân dân Việt Nam (như các dòng UAV-01, UAV-02, UAV-03 của Viện Kỹ thuật PK-KQ hay dòng Shikra của Viettel), giúp tăng tuổi thọ khung vỏ máy bay và giảm tỷ lệ tổn thất khi tiếp đất xuống đường băng dã chiến ngắn.
  • Lợi ích kinh tế - xã hội: Giúp giảm thiểu chi phí nhập khẩu các hệ thống hạ cánh quang học hoặc trạm hạ cánh vô tuyến đắt tiền từ nước ngoài; mở rộng khả năng ứng dụng UAV trong giám sát rừng, phòng chống cháy rừng, tuần tra đường dây tải điện cao thế và cứu hộ cứu nạn bão lũ tại các địa bàn đồi núi phức tạp.

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

  • Nghiên cứu sinh và Giảng viên ngành Kỹ thuật Điều khiển: Thụ hưởng khung phương pháp luận toán học chặt chẽ về phương pháp giải bài toán biên hai điểm bằng tham số liên tục và kỹ thuật tổng hợp hàm truyền kênh dọc máy bay.
  • Kỹ sư R&D Hàng không - Vũ trụ: Tiếp cận trực tiếp bộ thông số hình học, hệ số khí động và cấu trúc mạch vòng điều khiển đã được chuẩn hóa để ứng dụng vào thiết kế hệ thống lái tự động (Autopilot).
  • Lực lượng Tác chiến Biên phòng và Cảnh sát biển: Nâng cao độ sẵn sàng chiến đấu và khả năng thu hồi an toàn phương tiện trinh sát không người lái trong điều kiện khí tượng biển phức tạp và không gian hạ cánh bị hạn chế.

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 mở rộng Nguyên lý cực đại Pontryagin cho bài toán hạ cánh UAV bằng cách tích hợp thuật toán tham số liên tục (Continuation Method). Thay vì giải trực tiếp bài toán biên hai điểm (TPBVP) phi tuyến vốn dễ phân kỳ, luận án nhúng tham số biến dạng $s \in [0, 1]$ để chuyển dịch nghiệm một cách trơn tru từ hệ giải tích tuyến tính sang hệ động học phi tuyến thực tế, đồng thời áp đặt thành công ràng buộc quá tải đứng $1.0 \le n_y \le 3.5$.

2. Đột phá về phương pháp luận của luận án khi so sánh với ít nhất hai nghiên cứu quốc tế?
So với nghiên cứu của Barber et al. [31] (chỉ ước lượng sai lệch vị trí bằng quang học rồi đưa vào bộ điều khiển tuyến tính cổ điển) và nghiên cứu của Santos et al. [17] (dùng logic mờ Heuristic để hạ cánh), luận án tạo đột phá khi kết hợp bài toán tối ưu hóa biến phân toàn cục với bộ điều khiển bám phản hồi kết hợp ($n_y$ và $\vartheta$). Giải pháp này triệt tiêu hoàn toàn hiện tượng bốc mũi do gió đứng mà không đòi hỏi thiết bị camera đắt tiền hay cơ sở dữ liệu mờ phức tạp.

3. Phát hiện bất ngờ nhất từ dữ liệu mô phỏng là gì?
Phát hiện bất ngờ và có giá trị nhất là tác động phi đối xứng của gió đứng $W_y$. Khác với trực giác thông thường cho rằng gió dọc $W_x$ nguy hiểm nhất do làm đổi không tốc $V_r$, mô phỏng chỉ ra gió đứng từ dưới lên $W_y = +2\text{ m/s}$ mới là nguyên nhân chính gây mất an toàn: nó làm tăng góc tấn không tốc tức thời $\alpha_r = \alpha + \arctan(W_y/V_k)$, tạo lực nâng đột biến khiến UAV "nhảy cóc", trượt cự ly tiếp đất tới $65\text{ m}$ nếu không có mạch vòng phản hồi quá tải đứng $n_y$ dập tắt kịp thời.

4. Luận án có cung cấp quy trình tái lập nghiên cứu (Replication Protocol) hoàn chỉnh không?
Hoàn toàn đầy đủ. Luận án cung cấp chi tiết: (a) Toàn bộ bảng thông số khối lượng ($56.5\text{ kg}$), diện tích cánh ($1.05\text{ m}^2$), tọa độ trọng tâm; (b) Hệ phương trình trạng thái và hệ số đạo hàm khí động ($m_z^\alpha, m_z^{\omega_z}, m_z^{\delta_c}$); (c) Ma trận hàm truyền đạt kênh dọc; (d) Thuật toán số giải phương trình vi phân biến đồng trạng thái; và (e) Sơ đồ khối chi tiết các mạch vòng điều khiển trong Simulink.

5. Lộ trình phát triển 10 năm cho hướng nghiên cứu này được xác định ra sao?

  • Giai đoạn 1 (1-3 năm): Thử nghiệm phần cứng HIL và bay đánh giá thực tế trên mẫu UAV-70V tại bãi thử Quân chủng PK-KQ.
  • Giai đoạn 2 (3-6 năm): Mở rộng thuật toán cho bài toán hạ cánh 6-DOF không gian 3D trên mục tiêu di động (như sàn đáp tàu hải quân đang cơ động trên sóng).
  • Giai đoạn 3 (6-10 năm): Tích hợp trí tuệ nhân tạo (AI/Deep Reinforcement Learning) để tự động thích ứng với cấu hình khí động bị hỏng hóc một phần và hạ cánh tự chủ trong môi trường bão gió cực đoan.

Kết luận

  1. Luận án đã giải quyết trọn vẹn bài toán tối ưu hóa quỹ đạo hạ cánh cho UAV cỡ nhỏ cánh cố định bằng cách vận dụng sáng tạo Nguyên lý cực đại Pontryagin kết hợp phương pháp giải liên tục theo tham số, giải quyết triệt để sự phân kỳ toán học của bài toán biên hai điểm.
  2. Thiết lập thành công mô hình động lực học chuyển động dọc và hàm truyền đạt chuẩn hóa của dòng UAV-70V, làm rõ cơ chế tác động của vector gió dọc $W_x$ và gió đứng $W_y$ đến góc tấn và lực nâng khí động.
  3. Tổng hợp cấu trúc bộ điều khiển bám quỹ đạo kết hợp giữa tín hiệu quá tải đứng $n_y$ và góc chúc ngóc $\vartheta$ có khâu cản dịu $K_z$, chứng minh khả năng triệt tiêu sai lệch quỹ đạo vượt trội trong điều kiện nhiễu gió tầng thấp.
  4. Kết quả mô phỏng số xác nhận các chỉ tiêu kỹ thuật hoàn hảo: vận tốc tiếp đất $V_f = 31\text{ m/s}$, vận tốc lún $V_{yhc} = 0.38 \div 0.62\text{ m/s} \le 1.0\text{ m/s}$, sai số độ cao $\Delta y \le 0.12\text{ m} \le 0.3\text{ m}$, và sai số cự ly dọc $\Delta x \le 18.5\text{ m} \le 30\text{ m}$.
  5. Khởi xướng ba hướng nghiên cứu mới về điều khiển tối ưu thích nghi phi tuyến, hạ cánh không gian 3 chiều 6 bậc tự do, và hạ cánh tự hành trên bệ di động biển đảo.
  6. Công trình là một đóng góp học thuật và thực tiễn xuất sắc, nâng cao năng lực tự chủ công nghệ điều khiển tự động phương tiện bay không người lái phục vụ đắc lực cho an ninh - quốc phòng và kinh tế biển của Việt Nam.