Tổng quan về luận án

Nghiên cứu điều khiển các hệ thống cơ điện tử - khí động học phi tuyến đa biến luôn là một trong những bài toán phức tạp và thách thức nhất của chuyên ngành Kỹ thuật Điều khiển và Tự động hóa. Hệ Twin Rotor MIMO System (TRMS) do Feedback Instruments Ltd chế tạo là thiết bị thực nghiệm kinh điển đại diện cho lớp đối tượng cơ điện tử đa vào - đa ra (MIMO), mang bản chất động học tương tự máy bay trực thăng (Helicopter Dynamics) nhưng có điểm quay cố định với hai bậc tự do: chuyển vị góc trong mặt đứng $\alpha_v$ (Pitch angle, phạm vi hoạt động danh định $1.22\text{ rad}$) và chuyển vị góc trong mặt ngang $\alpha_h$ (Yaw angle, phạm vi $2.82\text{ rad}$). Đặc trưng kỹ thuật nổi bật của TRMS nằm ở tính phi tuyến bậc cao, tham số biến đổi theo thời gian, bất định mô hình, tác động xen kênh (cross-coupling effect) chéo mạnh mẽ giữa hai kênh chuyển động, cùng sự nhạy cảm trước các nhiễu khí động học bên ngoài.

Khoảng trống nghiên cứu (research gap) trọng tâm nằm ở sự thiếu hụt các mô hình toán học giải tích hoàn chỉnh phản ánh đúng cấu trúc hình học thực (như chiều dài chốt quay $h$, hiệu ứng bề mặt $\mu$, động học cơ cấu chấp hành động cơ một chiều) và sự hạn chế của các giải pháp điều khiển kinh điển (PID, LQR tuyến tính hóa tĩnh) khi đối mặt với tính phi tuyến mạnh cùng nhiễu bất định. Các công trình trước đây, tiêu biểu như tài liệu gốc của Feedback Co. (1998) hay nghiên cứu của Akbar Rahideh (2008), thường sử dụng mô hình Newton giản lược hoặc mô hình MPC rời rạc hóa điểm tĩnh vốn gây độ trễ tính toán lớn và sai lệch tĩnh đáng kể khi thực thi trên phần cứng thời gian thực.

Luận án của nghiên cứu sinh Đinh Văn Nghiệp (2018), dưới sự hướng dẫn khoa học của PGS.TS Nguyễn Như Hiển và GS.TS Nguyễn Doãn Phước tại Trường Đại học Kỹ thuật Công nghiệp – Đại học Thái Nguyên, đã thiết lập hệ thống câu hỏi và giả thuyết khoa học tường minh:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để xây dựng một mô hình toán học giải tích hoàn chỉnh cho TRMS trong không gian biến khớp có độ chính xác cao hơn mô hình Newton của nhà sản xuất bằng cách tính đến đầy đủ các yếu tố hình học và hiệu ứng khí động học phụ trợ?
  • Câu hỏi nghiên cứu 2 (RQ2): Cấu trúc thuật toán điều khiển tối ưu - thích nghi nào có khả năng kháng nhiễu, bù trừ sai lệch mô hình và đảm bảo tính bám quỹ đạo bền vững theo thời gian thực cho hệ phi tuyến liên tục TRMS?
  • Giả thuyết khoa học 1 (H1): Việc áp dụng giải tích cơ học Euler-Lagrange trong không gian biến khớp 4 bậc tự do $q = [\alpha_v, \alpha_h, \alpha_m, \alpha_t]^T$ kết hợp hiệu chỉnh hình học chốt quay sẽ giảm thiểu sai số mô hình hóa so với thực nghiệm xuống dưới $35%$ so với mô hình Newton truyền thống.
  • Giả thuyết khoa học 2 (H2): Thuật toán điều khiển tích hợp Chân trời lùi (Receding Horizon Control - RHC) với Bộ điều chỉnh Toàn phương Tuyến tính (LQR) dựa trên khai triển Jacobi cục bộ thời gian thực kết hợp mạch vòng phụ áp đặt mô-men sẽ loại bỏ triệt để hiện tượng sai lệch tĩnh và triệt tiêu ảnh hưởng xen kênh khi hệ chịu nhiễu ngoại cảnh.

Khung lý thuyết của luận án được xây dựng dựa trên sự giao thoa giữa Lý thuyết Cơ học giải tích Euler-Lagrange, Lý thuyết Điều khiển Tối ưu động RHC/MPC, Bộ điều chỉnh tối ưu LQR giải phương trình vi phân Riccati thời gian thực và Lý thuyết ổn định hệ phi tuyến liên tục. Phạm vi thực nghiệm của nghiên cứu được kiểm chứng trực tiếp trên dàn thử nghiệm TRMS Feedback 33-220 kết nối card xử lý tín hiệu số dSPACE DS1103, tạo nên bước đột phá về tính ứng dụng trong công nghệ điều khiển thiết bị bay không người lái (UAV) và các cơ cấu chấp hành cánh quạt phản lực.

Literature Review và Positioning

Các công trình nghiên cứu về TRMS trên thế giới và tại Việt Nam tập trung vào hai nhánh chính: Mô hình hóa động lực học và Thiết kế thuật toán điều khiển thích nghi - bền vững.

Về mô hình hóa, Feedback Instruments Ltd (1998) thiết lập mô hình toán dựa trên định luật Newton 2. Tuy nhiên, mô hình này bỏ qua chiều dài chốt quay $h$, giả định ma sát trơn tuyến tính và không tính đến sự thay đổi của cánh tay đòn do lực nâng tạo bởi cánh quạt chính khi góc nghiêng $\alpha_v$ thay đổi. Akbar Rahideh (2008) đã thực hiện khảo sát toàn diện các phương pháp mô hình hộp trắng (White-box via Newton và Euler-Lagrange), hộp đen (Black-box via Neural Networks) và hộp xám (Grey-box). Nghiên cứu của Rahideh chỉ ra rằng phương pháp Euler-Lagrange cho độ chính xác cao hơn Newton nhưng lại làm tăng độ phức tạp tính toán và chưa tích hợp động học phản ứng điện từ của động cơ chấp hành DC. Tại Việt Nam, nghiên cứu của tác giả Nguyễn Công Định và cộng sự (2016) đã ứng dụng mô hình Euler-Lagrange nhưng bỏ qua động học động cơ, dẫn đến việc thiết kế các vòng điều khiển vận tốc góc cánh quạt bằng PID vẫn tồn tại sai lệch động lớn.

Về phương pháp điều khiển, tồn tại cuộc tranh luận học thuật sâu sắc giữa hai trường phái:

  • Trường phái Điều khiển kinh điển & Trí tuệ nhân tạo: Sumit Kumar Pandey và Vijaya Laxmi (2014) đề xuất bộ điều khiển PID có khâu lọc vi phân; Abdullah Ates cùng các cộng sự ứng dụng hàm phi tuyến Sigmoid để chỉnh định tham số PID liên tục nhằm giảm giật điều khiển; Ahamdal (2014) và các nghiên cứu khác kết hợp PID với giải thuật di truyền (GA), tối ưu bầy đàn (PSO) hoặc Logic mờ (Fuzzy Logic). Nhược điểm cốt lõi của trường phái này là thiếu bảo chứng toán học chặt chẽ về tính ổn định bền vững khi biên độ nhiễu vượt ngưỡng huấn luyện, đồng thời thời gian quá độ còn lớn ($15\text{ s} - 20\text{ s}$).
  • Trường phái Điều khiển bền vững & Tối ưu: Jahed và Farrokhi thiết kế bộ điều khiển trượt (SMC) dựa trên việc phân tách TRMS thành hai hệ con, coi xen kênh là nhiễu bất định. Tuy nhiên, hiện tượng chattering và thời gian quá độ kéo dài (20 giây cho mặt đứng và 17 giây cho mặt ngang) làm suy giảm tuổi thọ phần cứng cơ khí. Bhanu Pratap ứng dụng LQR trên mô hình tuyến tính hóa tĩnh quanh điểm cân bằng, khiến chất lượng bám quỹ đạo phi tuyến suy giảm nghiêm trọng khi góc lệch lớn. Akbar Rahideh (2012) áp dụng Newton MPC phản hồi trạng thái rời rạc kết hợp bộ quan sát, song thuật toán này đòi hỏi chi phí tính toán cao, gây khó khăn cho việc xử lý thời gian thực.

Luận án định vị chính xác vị trí học thuật bằng cách khắc phục đồng thời các hạn chế trên: Tích hợp mô hình Euler-Lagrange hoàn chỉnh trong không gian biến khớp (bổ sung hiệu ứng hình học $h, h_1$ và động học cơ cấu chấp hành) với thuật toán điều khiển tối ưu RHC-LQR phi tuyến thời gian liên tục. Giải pháp này giúp tuyến tính hóa vi phân từng bước bằng ma trận Jacobi thời gian thực thay vì tuyến tính hóa điểm làm việc cố định, kết hợp mạch vòng dòng điện/mô-men giúp nâng cao chất lượng bám quỹ đạo vượt trội so với các công trình quốc tế cùng thời điểm.

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

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

Luận án đóng góp quan trọng vào lý thuyết điều khiển tự động hiện đại bằng việc mở rộng và kết hợp thành công hai trụ cột lý thuyết: Lý thuyết Chân trời lùi (Receding Horizon Control - RHC) và Bộ điều chỉnh Toàn phương Tuyến tính (Linear Quadratic Regulator - LQR) cho lớp đối tượng phi tuyến liên tục có tham số biến đổi theo thời gian (Linear Time-Varying - LTV) chịu nhiễu bất định:

$$\dot{x} = f(x, u) + \Delta_k \approx A_k x + B_k u + d_k + \Delta_k$$

Khác với các phương pháp MPC truyền thống vốn giải bài toán tối ưu hóa phi tuyến tĩnh qua quy hoạch phi tuyến (NLP) phức tạp, mô hình lý thuyết đề xuất thực hiện xấp xỉ chuỗi Taylor bậc nhất xung quanh trạng thái tức thời $(x_k, u_{k-1})$ tại mỗi chu kỳ trích mẫu rất nhỏ $t_k \le t \le t_k + \delta_k$:

$$A_k = \left. \frac{\partial f}{\partial x} \right|{(x_k, u{k-1})}, \quad B_k = \left. \frac{\partial f}{\partial u} \right|{(x_k, u{k-1})}, \quad d_k = f(x_k, u_{k-1}) - A_k x_k - B_k u_{k-1}$$

Luật điều khiển tối ưu phản hồi trạng thái trong cửa sổ dự báo hữu hạn được dẫn xuất thông qua việc giải phương trình vi phân ma trận Riccati thời gian thực:

$$u(t) = -R^{-1} B_k^T P(t) x(t) + u_{\text{comp}}(t)$$

Nghiên cứu chứng minh rằng sự kết hợp này đảm bảo tính ổn định tiệm cận của hệ thống theo tiêu chuẩn Lyapunov, đồng thời hàm mục tiêu toàn phương giúp hạn chế tối đa năng lượng điều khiển đặt lên các động cơ truyền động cánh quạt trong điều kiện có chặn biên độ điện áp phần ứng $u \in [u_{\min}, u_{\max}]$.

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

Khung phân tích của luận án tích hợp liên hoàn 3 mô thức lý thuyết:

  1. Mô thức Động lực học Euler-Lagrange Đa vật thể: Biểu diễn động học 4 bậc tự do trong không gian biến khớp $q = [\alpha_v, \alpha_h, \alpha_m, \alpha_t]^T$, bao gồm đầy đủ ma trận quán tính phi tuyến $M(\alpha_v)$, véc-tơ lực Coriolis và ly tâm $C(q, \dot{q})$, véc-tơ trọng trường $G(\alpha_v)$ và véc-tơ mô-men ngoại lực $M_{\text{ext}}$.
  2. Mô thức Tuyến tính hóa Cục bộ Biến thiên Thời gian (LTV-Jacobian Approximation): Chuyển hóa bài toán điều khiển phi tuyến toàn cục phức tạp thành chuỗi các bài toán LQR biến thiên liên tục, giúp loại bỏ hoàn toàn nhu cầu sử dụng các bộ tối ưu hóa số nặng nề trên máy tính nhúng.
  3. Mô thức Bù đảo Động học Cơ cấu Chấp hành (Inverse Actuator Dynamic Mapping): Bổ sung mạch vòng mô-men nội suy trực tiếp từ dòng điện phần ứng động cơ DC, biến đổi tín hiệu điều khiển mô-men yêu cầu thành điện áp kích thích tức thời, triệt tiêu độ trễ cảm ứng điện từ $L_m, L_t$.

Điều kiện biên (boundary conditions) của mô hình được xác lập chặt chẽ: Dải góc pitch danh định $\alpha_v \in [-1.22\text{ rad}, +1.22\text{ rad}]$, dải góc yaw $\alpha_h \in [-2.82\text{ rad}, +2.82\text{ rad}]$, giới hạn điện áp điều khiển phần ứng $|u_{am}| \le 24\text{ VDC}, |u_{at}| \le 24\text{ VDC}$, và tần số lấy mẫu hệ thống $T_s = 0.01\text{ s}$ trên nền tảng DSP.

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 khoa học thực chứng (Positivism) kết hợp chủ nghĩa hiện thực phản biện (Critical Realism). Thiết kế nghiên cứu là sự kết hợp khép kín giữa mô hình hóa giải tích toán học, mô phỏng số học trên môi trường MATLAB/Simulink và kiểm chứng thực nghiệm trên dàn thử nghiệm vật lý thời gian thực.

       +----------------------------------------------------------------------+
       |                    QUY TRÌNH NGHIÊN CỨU THỰC NGHIỆM                  |
       +----------------------------------+-----------------------------------+
                                          |
       +----------------------------------v-----------------------------------+
       | 1. XÂY DỰNG MÔ HÌNH TOÁN EULER-LAGRANGE (4 BẬC TỰ DO)                |
       |    - Thiết lập hệ tọa độ quy chiếu động O1xyz & cố định OXYZ        |
       |    - Tính toán động năng K, thế năng V và phương trình Lagrange      |
       +----------------------------------+-----------------------------------+
                                          |
       +----------------------------------v-----------------------------------+
       | 2. MÔ PHỎNG & ĐỐI SÁNH TRÊN MATLAB/SIMULINK                          |
       |    - So sánh đáp ứng bước 1(t) & hàm sin(t)                          |
       |    - Đối chứng mô hình đề xuất vs Mô hình Newton Feedback Co.        |
       +----------------------------------+-----------------------------------+
                                          |
       +----------------------------------v-----------------------------------+
       | 3. THIẾT KẾ THUẬT TOÁN ĐIỀU KHIỂN RHC VỚI LQR                        |
       |    - Tuyến tính hóa Jacobi trực tuyến tại điểm (xk, uk-1)            |
       |    - Thiết lập thuật toán giải phương trình Riccati thời gian thực    |
       +----------------------------------+-----------------------------------+
                                          |
       +----------------------------------v-----------------------------------+
       | 4. TRIỂN KHAI THỰC NGHIỆM TRÊN PHẦN CỨNG dSPACE DS1103               |
       |    - Ghép nối Encoder quang, máy phát tốc, cảm biến dòng Hall       |
       |    - Bơm nhiễu khí động học cưỡng bức bằng quạt công suất phụ trợ   |
       +----------------------------------------------------------------------+

Thiết kế đa tầng được phân định rành mạch:

  • Tầng 1 (Physical System Level): Khảo sát các thông số vật lý thực của TRMS Feedback 33-220, bao gồm khối lượng thanh ngang $m_{T1}$, khối lượng đối trọng $m_{cb} = 0.068\text{ kg}$, chiều dài thanh đối trọng $l_b = 0.29\text{ m}$, cánh tay đòn cánh quạt chính $l_m = 0.24\text{ m}$, cánh tay đòn đuôi $l_t = 0.25\text{ m}$.
  • Tầng 2 (Actuator & Sensor Level): Tích hợp cảm biến góc Encoder quang 2 kênh A/B độ phân giải cao, máy phát tốc DC đo vận tốc góc cánh quạt $\omega_m, \omega_t$, và mô-đun cảm biến dòng điện độ nhạy cao đo $i_{am}, i_{at}$.
  • Tầng 3 (Computational & Control Level): Xử lý thuật toán điều khiển RHC-LQR trên vi xử lý trung tâm dSPACE DS1103 PowerPC 750GX chạy ở xung nhịp $1\text{ GHz}$ kết hợp DSP phụ trợ TMS320F240.

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 nghiệm tuân thủ các giao thức chuẩn mực:

  • Giao thức cảm biến và chuyển đổi tín hiệu: Tín hiệu chuyển vị góc được đọc qua bộ đếm xung quang học Encoder đưa vào khối ghép nối số, tín hiệu vận tốc góc từ máy phát tốc và dòng điện phần ứng được đưa qua bộ chuyển đổi tương tự - số (A/D Converter 16-bit) với bộ lọc thông thấp tương tự nhằm triệt tiêu nhiễu tần số cao của chuyển mạch công suất PWM.
  • Triangulation (Tam giác đạc): Dữ liệu phản hồi được kiểm chứng chéo qua 3 kênh độc lập: Dữ liệu mô hình lý thuyết giải tích $\rightarrow$ Dữ liệu mô phỏng S-Function Simulink $\rightarrow$ Dữ liệu đo đạc thực tế từ ControlDesk trên cùng một chuỗi tín hiệu kích thích (Step, Sine, Mixed Sin-Step).
  • Thử nghiệm áp lực và kháng nhiễu (Robustness Test Protocol): Luận án trang bị hệ thống quạt thổi gió cưỡng bức có điều tần độc lập đặt trực diện vào TRMS để tạo nhiễu khí động học phi tuyến ngẫu nhiên, mô phỏng các đợt gió giật đột ngột tác động lên thiết bị bay.

Data và phân tích

Toàn bộ thông số định lượng của hệ thống TRMS được nhận dạng chính xác và cấu hình trong bảng tham số chuẩn:

Tham số vật lý Đại lượng mô tả Giá trị danh định Đơn vị
$g$ Gia tốc trọng trường $9.81$ $\text{m/s}^2$
$m_{cb}$ Khối lượng quả đối trọng $0.068$ $\text{kg}$
$m_{mr}, m_{tr}$ Khối lượng động cơ chính / đuôi $0.22 / 0.18$ $\text{kg}$
$l_m, l_t$ Chiều dài đòn bẩy cánh quạt chính / đuôi $0.24 / 0.25$ $\text{m}$
$l_{cb}$ Khoảng cách từ đối trọng đến chốt quay $0.05 - 0.29$ $\text{m}$
$J_v, J_h$ Tổng mô-men quán tính trục đứng / ngang Phụ thuộc $\alpha_v$ $\text{kg}\cdot\text{m}^2$
$k_g$ Hệ số hiệu ứng Gyroscope $0.05$ $\text{N}\cdot\text{m}\cdot\text{s/rad}$
$B_v, B_h$ Hệ số ma sát nhớt khớp quay đứng / ngang $0.01 / 0.01$ $\text{N}\cdot\text{m}\cdot\text{s/rad}$
$F_v, F_h$ Ma sát trượt Coulomb khớp đứng / ngang $0.005 / 0.005$ $\text{N}\cdot\text{m}$

Hệ thống phần mềm phân tích chuyên sâu bao gồm MATLAB/Simulink R2015b, Real-Time Interface (RTI1103), và dSPACE ControlDesk v5.0. Các bài toán tối ưu hóa ma trận trọng số LQR $Q = \text{diag}(q_1, q_2, \dots, q_n)$ và $R = \text{diag}(r_1, r_2)$ được phân tích độ nhạy (sensitivity analysis) để tìm điểm cân bằng tối ưu giữa tốc độ đáp ứng và biên độ điện áp kích từ.

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

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

Quá trình mô phỏng và thực nghiệm kiểm chứng trên dàn thử nghiệm TRMS đã mang lại 4 phát hiện đột phá mang tính định lượng rõ rệt:

+-----------------------------------------------------------------------------+
|               SO SÁNH SAI LỆCH GÓC YAW αh DƯỚI TÍN HIỆU BƯỚC 1(t)           |
+-----------------------------------------------------------------------------+
| Mô hình Newton (Feedback Co.) : [========================] 0.9 rad          |
| Mô hình Euler-Lagrange đề xuất : [========] 0.3 rad (-66.7% Error)           |
|                                                                             |
|                    THỜI GIAN QUÁ ĐỘ BÁM QUỸ ĐẠO GÓC PITCH                   |
+-----------------------------------------------------------------------------+
| Bộ điều khiển Trượt SMC (Jahed): [===================] 20.0 giây            |
| Thuật toán RHC-LQR của Luận án:  [===] 3.2 giây (-84.0% Settling Time)      |
+-----------------------------------------------------------------------------+
  1. Độ chính xác vượt trội của Mô hình Euler-Lagrange cải tiến: Khi kích thích hệ thống bằng tín hiệu hàm bước nhảy đơn vị $1(t)$, trích xuất dữ liệu thực nghiệm cho thấy: "Sai lệch vị trí góc $\alpha_h$ giữa đối tượng thực với mô hình nhà sản xuất cung cấp là $0.9\text{ rad}$, còn với mô hình nghiên cứu sinh xây dựng là $0.3\text{ rad}$". Mô hình mới đã giảm thiểu sai số mô hình hóa ở kênh mặt phẳng ngang tới $66.7%$. Hiện tượng lệch pha nghiêm trọng từng xuất hiện ở mô hình Newton đã được khắc phục hoàn toàn nhờ việc tính toán chính xác cánh tay đòn chốt quay $h$ và hàm tương quan ma sát phi tuyến $F_h \text{sign}(\dot{\alpha}_h)$.
  2. Khắc phục triệt để hiện tượng trôi sai số tĩnh (Zero Steady-state Error): Trong các thử nghiệm bám quỹ đạo phức tạp kết hợp đồng thời (Sin-Sin, Sin-Step, Step-Step), cấu trúc RHC-LQR tích hợp bù sai số bám ở bước trước $\Delta w_k$ đã triệt tiêu hoàn toàn sai lệch tĩnh, đưa sai số bám xác lập $e_{ss} \to 0$, trong khi các bộ điều khiển PID kinh điển hoặc LQR tĩnh luôn tồn tại độ lệch tĩnh từ $0.08\text{ rad}$ đến $0.15\text{ rad}$ do ảnh hưởng của trọng lực quả đối trọng.
  3. Triệt tiêu hiện tượng xen kênh và hiệu ứng con quay hồi chuyển: Khi góc pitch $\alpha_v$ thay đổi đột ngột với vận tốc lớn, thành phần ngẫu lực Gyroscope $k_g \cdot \omega_m \cdot \dot{\alpha}_v \cos(\alpha_v)$ tác động trực tiếp gây rung lắc dữ dội lên góc yaw $\alpha_h$. Thuật toán RHC-LQR với ma trận Jacobi cập nhật trực tuyến đã dự báo trước tác động chéo này và xuất tín hiệu bù trừ tức thời lên động cơ đuôi, giữ cho độ nhấp nhô góc yaw không vượt quá $\pm 0.02\text{ rad}$.
  4. Tối ưu hóa thời gian quá độ và khả năng kháng nhiễu cưỡng bức: So với công trình của Jahed và Farrokhi (2008) có thời gian quá độ lên tới $20\text{ s}$ (mặt đứng) và $17\text{ s}$ (mặt ngang), bộ điều khiển đề xuất rút ngắn thời gian quá độ xuống dưới $3.5\text{ s}$ đối với cả hai kênh chuyển động. Khi kích hoạt quạt tạo nhiễu khí động học công suất lớn, hệ thống nhanh chóng tái lập trạng thái cân bằng trong vòng $1.2\text{ s}$ mà không xảy ra hiện tượng bão hòa điện áp phần ứng ($u_{am}, u_{at} \le 24\text{ V}$).

Implications đa chiều

  • Về mặt Lý thuyết: Luận án mở ra phương pháp luận thiết kế điều khiển tối ưu LTV-RHC khả thi cho lớp hệ thống phi tuyến MIMO có ma trận quán tính biến thiên, lấp đầy khoảng cách giữa lý thuyết MPC phi tuyến trừu tượng và khả năng hiện thực hóa thuật toán giải Riccati thời gian thực.
  • Về mặt Phương pháp luận: Cung cấp quy trình định chuẩn mô hình hộp trắng đa vật thể có thể chuyển giao cho việc mô hình hóa các thiết bị bay nhiều cánh quạt (Quadrotor, Hexacopter, Tilt-rotor UAVs).
  • Về mặt Thực tiễn Kỹ thuật: Cấu trúc điều khiển phân tầng (RHC-LQR bên ngoài kết hợp mạch vòng dòng điện/mô-men bên trong) giúp bảo vệ phần cứng chấp hành, giảm thiểu hao mòn cơ khí do hiện tượng chattering của cánh quạt và tối ưu hóa hiệu suất tiêu thụ năng lượng của pin trên các thiết bị bay tự hành.
  • Về mặt Chính sách và Công nghệ Quốc gia: Đóng góp giải pháp tự chủ công nghệ điều khiển thiết bị bay không người lái phục vụ công tác giám sát nông - lâm nghiệp, cứu hộ cứu nạn và an ninh quốc phòng tại Việt Nam.

Limitations và Future Research

Dù đạt được những kết quả xuất sắc, luận án vẫn ghi nhận các giới hạn học thuật cần tiếp tục hoàn thiện:

  • Giới hạn 1 (Giả thiết vật rắn tuyệt đối): Mô hình giải tích Euler-Lagrange vẫn coi các thanh nối và cánh quạt là liên kết cứng tuyệt đối, chưa xét đến biến dạng đàn hồi (aeroelastic flexibility) của thân cánh quạt khi vận hành ở dải vận tốc góc siêu cao.
  • Giới hạn 2 (Ước lượng ma sát chưa thích nghi online): Các hệ số ma sát Coulomb và ma sát nhớt ($B_v, B_h, F_v, F_h$) được xác định bằng phương pháp nhận dạng thực nghiệm cố định, chưa có cơ chế ước lượng thích nghi trực tuyến khi nhiệt độ ổ bi khớp quay tăng cao trong quá trình vận hành liên tục.
  • Giới hạn 3 (Ràng buộc biến trạng thái): Thuật toán LTV-RHC giải phương trình Riccati hiện thời xử lý trực tiếp ràng buộc biên độ điện áp điều khiển đầu vào $u$, nhưng chưa tích hợp trực tiếp các ràng buộc bất đẳng thức ngặt trên biến trạng thái trạng thái (State constraints: $\alpha_{\min} \le \alpha(t) \le \alpha_{\max}$) vào hàm mục tiêu tối ưu.

Chương trình nghiên cứu tương lai (Future Research Agenda) được định hình qua 4 hướng chính:

  1. Tích hợp mạng nơ-ron sâu hoặc bộ lọc Kalman mở rộng thích nghi (Adaptive EKF) để nhận dạng trực tuyến các tham số ma sát biến thiên và hệ số suy giảm lực nâng khí động học theo nhiệt độ môi trường.
  2. Mở rộng thuật toán RHC sang không gian 6 bậc tự do (6-DOF) áp dụng trực tiếp cho mô hình máy bay trực thăng không người lái kích thước lớn hoạt động trong điều kiện gió bão cấp độ cao.
  3. Nghiên cứu giải pháp nhúng trực tiếp giải thuật giải Riccati song song trên phần cứng vi mạch chuyên dụng FPGA/SoC nhằm nâng tần số lấy mẫu lên mức $1\text{ kHz} - 10\text{ kHz}$.
  4. Phát triển cơ chế điều khiển chịu lỗi (Fault-Tolerant Control - FTC) tự động bù đắp khi một trong hai động cơ truyền động bị suy giảm công suất đột ngột.

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

  • Ảnh hưởng Học thuật: Luận án đóng góp nguồn tài liệu tham khảo mẫu mực cho các chương trình đào tạo sau đại học ngành Tự động hóa tại Đại học Thái Nguyên, Đại học Bách khoa Hà Nội và các viện nghiên cứu chuyên ngành. Các công bố khoa học từ đề tài trên các tạp chí uy tín đã tạo lập diễn đàn học thuật về điều khiển hệ phi tuyến khí động học.
  • Chuyển dịch Công nghệ Công nghiệp & R&D: Mô hình điều khiển RHC-LQR cung cấp giải pháp lõi cho các đơn vị nghiên cứu phát triển thiết bị bay không người lái (UAV/Drone) công nghiệp, hệ thống camera gimbal chống rung 3 trục, và các cánh tay robot chuyển động tốc độ cao.
  • Giá trị Xã hội và Môi trường: Nâng cao hiệu quả điều khiển chính xác của các phương tiện bay tự hành ứng dụng trong quan trắc rừng, phun thuốc bảo vệ thực vật chính xác trong nông nghiệp công nghệ cao, và tuần tra biên giới, góp phần giảm thiểu rủi ro tai nạn lao động và tiết kiệm tài nguyên năng lượng.

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

+--------------------+--------------------------------------------------------+
| ĐỐI TƯỢNG HƯỞNG LỢI| GIÁ TRỊ VÀ LỢI ÍCH CHUYÊN BIỆT TIẾP NHẬN               |
+--------------------+--------------------------------------------------------+
| Nghiên cứu sinh &  | Tiếp cận khung giải tích Euler-Lagrange 4-DOF chuẩn xác|
| Học viên Cao học   | và thuật toán LTV-RHC giải Riccati thời gian thực.     |
+--------------------+--------------------------------------------------------+
| Giảng viên &       | Khung bài giảng và tài liệu thực nghiệm mẫu trên nền   |
| Nhà khoa học       | tảng phần cứng chuẩn công nghiệp dSPACE DS1103.        |
+--------------------+--------------------------------------------------------+
| Kỹ sư R&D Drone &  | Thuật toán bám quỹ đạo kháng nhiễu và bù trễ động cơ   |
| Thiết bị bay tự hành| giúp nâng cao độ ổn định bay của UAV trong gió mạnh.   |
+--------------------+--------------------------------------------------------+
| Chuyên gia Phân tích| Cơ sở dữ liệu và phương pháp luận đối chuẩn giữa mô    |
| & Mô phỏng Cơ học  | phỏng S-Function và thực nghiệm vật lý thực.          |
+--------------------+--------------------------------------------------------+

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à nó mở rộng lý thuyết nào? Trả lời: Đóng góp độc đáo nhất là việc phát triển thành công thuật toán điều khiển Chân trời lùi tích hợp Bộ điều chỉnh Toàn phương Tuyến tính (RHC với LQR) thời gian liên tục cho hệ phi tuyến có nhiễu và sai lệch mô hình. Luận án đã mở rộng Lý thuyết Điều khiển Tối ưu cổ điển bằng cách áp dụng phép vi phân ma trận Jacobi trực tuyến xung quanh điểm làm việc động $(x_k, u_{k-1})$, biến bài toán phi tuyến phức tạp thành chuỗi bài toán LTV tối ưu giải được qua phương trình đại số/vi phân Riccati trong từng chu kỳ lấy mẫu, loại bỏ hoàn toàn sự phụ thuộc vào các thuật toán quy hoạch phi tuyến nặng nề.

2. Đổi mới về mặt phương pháp nghiên cứu so với các công trình quốc tế tiền nhiệm? Trả lời: So với công trình của Feedback Co. (1998) vốn chỉ dùng mô hình Newton thiếu chốt quay $h$, và công trình của Akbar Rahideh (2012) sử dụng MPC rời rạc hóa điểm tĩnh kèm bộ quan sát trạng thái gây trễ, luận án đã thiết lập mô hình Euler-Lagrange hoàn chỉnh 4 bậc tự do trong không gian biến khớp, đồng thời bổ sung mạch vòng phản hồi dòng điện/mô-men để bù động học điện từ của động cơ chấp hành DC. Toàn bộ thuật toán được biên dịch và chạy thực thời gian thực ở tần số cao trên phần cứng DSP dSPACE DS1103.

3. Phát hiện thực nghiệm nào bất ngờ và có ý nghĩa thống kê rõ nét nhất? Trả lời: Phát hiện bất ngờ nhất là sai lệch mô hình hóa ở kênh mặt ngang $\alpha_h$ của mô hình Newton Feedback Co. lên tới $0.9\text{ rad}$ (tương đương sai số góc trên $51^\circ$) dưới kích thích bước nhảy $1(t)$, trong khi mô hình Euler-Lagrange của luận án khống chế sai lệch chỉ còn $0.3\text{ rad}$. Hơn thế nữa, bộ điều khiển RHC-LQR đã cắt giảm thời gian quá độ từ $20\text{ s}$ (ở các bộ điều khiển trượt SMC trước đây) xuống còn $3.2\text{ s} - 3.5\text{ s}$, triệt tiêu hoàn toàn dao động do hiệu ứng Gyroscope gây ra.

4. Luận án có cung cấp đầy đủ giao thức tái lập thực nghiệm (Replication Protocol) không? Trả lời: Có. Toàn bộ sơ đồ khối ghép nối phần cứng DS1103, thư viện kết nối thời gian thực RTI1103_Simulink, mạch đo dòng điện Hall độ nhạy cao, module phát tốc, module công suất PWM, giao diện giám sát dSPACE ControlDesk, bảng tham số cơ lý tính chi tiết của TRMS và mã nguồn lập trình S-Function đều được đính kèm chi tiết trong phần phụ lục của luận án.

5. Chương trình nghị sự nghiên cứu 10 năm tới mở ra từ luận án là gì? Trả lời: Mở rộng mô hình từ hệ 2 bậc tự do TRMS sang hệ 6 bậc tự do hoàn chỉnh cho UAV tự hành; tích hợp trí tuệ nhân tạo (Deep Reinforcement Learning) phối hợp với RHC-LQR để thích nghi tham số khí động học trong môi trường gió lốc cực đoan; và chế tạo chip điều khiển chuyên dụng (ASIC/FPGA) nhúng thuật toán giải Riccati thời gian thực ở cấp micro-giây.

Kết luận

Luận án Tiến sĩ Kỹ thuật của tác giả Đinh Văn Nghiệp đã hoàn thành xuất sắc các mục tiêu nghiên cứu khoa học đặt ra, tạo nên những đóng góp học thuật và thực tiễn vững chắc:

  1. Xây dựng thành công mô hình toán động lực học giải tích Euler-Lagrange hoàn chỉnh cho hệ Twin Rotor MIMO System trong không gian biến khớp $q = [\alpha_v, \alpha_h, \alpha_m, \alpha_t]^T$, giải quyết triệt để các khiếm khuyết của mô hình Newton kinh điển thông qua việc tính đến chiều dài chốt quay $h$, hiệu ứng bề mặt và ma sát phi tuyến.
  2. Thiết kế thành công thuật toán điều khiển tối ưu thích nghi RHC kết hợp LQR cho hệ phi tuyến liên tục, cho phép tuyến tính hóa vi phân Jacobi từng bước và giải phương trình vi phân Riccati thời gian thực để tạo tín hiệu điều khiển bám quỹ đạo chính xác cao.
  3. Đề xuất cấu trúc điều khiển phân tầng bổ sung mạch vòng áp đặt mô-men giúp bù trừ triệt để quán tính điện từ và động học phi tuyến của cơ cấu chấp hành động cơ DC.
  4. Hiện thực hóa và kiểm chứng thành công toàn diện trên dàn thử nghiệm thực dSPACE DS1103, chứng minh tính vượt trội về tốc độ quá độ (giảm trên $80%$ thời gian xác lập so với SMC), triệt tiêu sai lệch tĩnh ($e_{ss} \to 0$) và duy trì tính bền vững tuyệt đối trước nhiễu khí động học cưỡng bức.
  5. Định hình bước tiến mới về phương pháp luận điều khiển tối ưu hóa thời gian thực cho các hệ thống cơ điện tử - khí động học phức tạp, mở ra các hướng ứng dụng chuyển giao công nghệ quy mô lớn cho công nghiệp hàng không vũ trụ và thiết bị bay không người lái tại Việt Nam.