Tổng quan về luận án

Bối cảnh khoa học của công trình nghiên cứu bắt nguồn từ các yêu cầu cấp bách trong chiến lược phát triển bền vững kinh tế biển, gắn liền với định hướng hiện thực hóa mục tiêu hiện đại hóa công nghệ hàng hải theo Nghị quyết số 36-NQ/TW. Trong kỷ nguyên chuyển đổi số và tự động hóa hàng hải, việc phát triển các hệ thống điều khiển tự hành cho phương tiện nổi bề mặt (Autonomous Surface Vessels - ASV) và tàu vận tải cỡ lớn đòi hỏi những thuật toán điều động bám quỹ đạo đạt độ chính xác cao, bảo đảm tính ổn định động học trong môi trường biển biến động. Tính tiên phong của luận án thể hiện ở việc tiên phong tích hợp đồng thời ba bài toán phức tạp và xung đột trong lý thuyết điều khiển tự động: cơ cấu chấp hành bị thiếu (underactuated), ràng buộc vật lý phi tuyến nghiêm ngặt ở đầu vào (input constraints: biên độ và tốc độ góc bẻ lái, lực đẩy chân vịt), cùng với sự hiện diện của các thành phần bất định hàm không cấu trúc (lumped functional uncertainties) và nhiễu ngẫu nhiên từ môi trường biển.

Khoảng trống nghiên cứu (research gap) cốt lõi được định vị từ các hạn chế mang tính hệ thống của y văn quốc tế và trong nước. Đa số các nghiên cứu kinh điển về điều khiển tàu thiếu cơ cấu chấp hành (như Do et al., 2004; Fossen, 2011) dựa trên kỹ thuật Backstepping hoặc Điều khiển chế độ trượt (Sliding Mode Control - SMC) thường giả định tín hiệu điều khiển không bị chặn vô hạn, hoặc đẩy bài toán ràng buộc cho các khối bão hòa vật lý (saturation devices) xử lý thụ động. Điều này dẫn tới hiện tượng tích lũy sai lệch (integrator windup) và suy giảm ổn định nghiêm trọng. Mặt khác, các công trình tiếp cận theo Điều khiển dự báo phi tuyến (Nonlinear Model Predictive Control - NMPC, tiêu biểu như Zheng et al., 2014) lại vấp phải rào cản về khối lượng tính toán khổng lồ, khó hội tụ thời gian thực và không xử lý triệt để các bất định động lực học phi tuyến cao bậc.

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 tường minh:

  1. Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để giải quyết tối ưu bài toán bám quỹ đạo mặt phẳng ngang (3-DOF) cho tàu thiếu cơ cấu chấp hành khi tín hiệu góc bẻ lái $\delta$ bị giới hạn nghiêm ngặt trong dải $[-35^\circ, +35^\circ]$ (khi điều động) và $[-10^\circ, +10^\circ]$ (khi chạy tự động)?
  2. Câu hỏi nghiên cứu 2 (RQ2): Có thể thiết lập một cấu trúc điều khiển tuyến tính hóa từng đoạn dọc trục thời gian để giảm độ phức tạp tính toán của NMPC mà vẫn bảo toàn tính ổn định tiệm cận của hệ phi tuyến hay không?
  3. Câu hỏi nghiên cứu 3 (RQ3): Làm thế nào để ước lượng và bù trừ trực tuyến thành phần bất định hàm tổng hợp $\Delta(t, \nu)$ bao gồm sai lệch mô hình thủy động lực học bậc cao và nhiễu sóng gió mà không cần đo đạc trực tiếp?

Hệ giả thuyết tương ứng:

  • Giả thuyết 1 (H1): Thuật toán điều khiển dự báo trượt dọc trục thời gian (RHC/MPC) dựa trên kỹ thuật tuyến tính hóa từng đoạn mô hình song tuyến cho phép chuyển đổi bài toán tối ưu phi tuyến thành bài toán Quy hoạch toàn phương (Quadratic Programming - QP), giải quyết trực tiếp ràng buộc đầu vào với thời gian tính toán chu kỳ thực dưới $100\text{ ms}$.
  • Giả thuyết 2 (H2): Bộ ước lượng bù bất định đầu vào kết hợp cùng Bộ lọc Kalman mở rộng (Extended Kalman Filter - EKF) theo nguyên lý tách (separation principle) bảo đảm tính ổn định tiệm cận theo nghĩa Lyapunov, triệt tiêu sai lệch bám vị trí $(e_x, e_y)$ và hướng đi $e_\psi$ về lân cận gốc tọa độ dưới tác động của nhiễu biển cấp 5.

Khung lý thuyết nền tảng của nghiên cứu tích hợp lý thuyết động lực học hàng hải chuẩn SNAME (1950), lý thuyết điều khiển tối ưu có ràng buộc (Constrained Optimal Control), và lý thuyết ổn định Lyapunov cho hệ thống chuyển mạch/hệ song tuyến. Phạm vi nghiên cứu giới hạn trên tàu nổi choán nước 3 bậc tự do mặt phẳng ngang (trượt dọc $u$, trượt ngang $v$, quay trở $r$), chịu tác động của cơ cấu đẩy chính $\tau_u$ và bánh lái $\tau_r$, vận hành trong điều kiện môi trường biển có sóng gió dòng chảy không vượt quá cấp 5 Beaufort.

Literature Review và Positioning

Tổng quan y văn quốc tế cho thấy sự phát triển của bài toán điều khiển chuyển động tàu thủy trải qua ba dòng nghiên cứu chính. Dòng nghiên cứu thứ nhất dựa trên mô hình tuyến tính hóa cục bộ hoặc mô hình tham số rút gọn Nomoto (Nomoto, 1957; Norrbin, 1970) ứng dụng điều khiển PID kinh điển và điều khiển tối ưu tuyến tính LQR/LQG (Fossen, 2002). Mặc dù đơn giản, hướng tiếp cận này bỏ qua hoàn toàn các liên kết chéo phi tuyến giữa chuyển động trượt ngang $v$ và vận tốc góc quay trở $r$, dẫn đến sự mất ổn định khi tàu vận hành ở dải tốc độ cao hoặc khi chuyển hướng đột ngột.

Dòng nghiên cứu thứ hai tập trung vào các kỹ thuật điều khiển phi tuyến hiện đại. Các công trình tiêu biểu của Do, Jiang và Pan (2002, 2004) đã áp dụng kỹ thuật Backstepping kết hợp hàm Lyapunov điều khiển tàu thiếu cơ cấu chấp hành bám quỹ đạo cong. Tuy nhiên, tranh luận học thuật nảy sinh xung quanh "sự bùng nổ số hạng" (explosion of terms) khi lấy đạo hàm các tín hiệu điều khiển ảo. Để khắc phục, Lee et al. (2010) và Chen et al. (2013) đề xuất Điều khiển bề mặt động (DSC) và Điều khiển trượt thích nghi (Adaptive SMC). Nhóm tác giả Wang et al. (2015) và Liu et al. (2016) tiếp tục tích hợp mạng nơ-ron nhân tạo (RBF Neural Networks) để xấp xỉ các hệ số thủy động lực học phi tuyến $D_n(\nu)$. Điểm xung đột lý thuyết chính giữa hai trường phái này là: phương pháp SMC/Backstepping tuy có tính bền vững cao nhưng gây ra hiện tượng rung tần số cao (chattering) ở cơ cấu bánh lái và hoàn toàn bất lực trong việc trực tiếp xử lý các ràng buộc biên độ/tốc độ cơ cấu chấp hành trong bài toán tối ưu.

+-------------------------------------------------------------------------------+
|                       Y VĂN ĐIỀU KHIỂN TÀU THỦY (LITERATURE MAP)               |
+-------------------------------------------------------------------------------+
                                      |
         +----------------------------+----------------------------+
         |                                                         |
[Trường phái Tuyến tính]                                 [Trường phái Phi tuyến]
Nomoto (1957), Fossen (2002)                             Do & Pan (2004), Chen (2013)
- PID, LQR, LQG                                          - Backstepping, SMC, DSC
- Nhược: Bỏ qua động học chéo                            - Ưu: Bền vững cao
                                                         - Nhược: Không tối ưu ràng buộc
                                                                   |
                                      +----------------------------+
                                      |
                          [Trường phái Dự báo NMPC]
                          Zheng (2014), Li (2018)
                          - Ưu: Xử lý ràng buộc đa biến
                          - Nhược: Tính toán phi tuyến nặng nề
                                      |
                                      v
          +-------------------------------------------------------+
          |             KHOẢNG TRỐNG VÀ VỊ TRÍ LUẬN ÁN            |
          |  Successive Linearization MPC + EKF + Uncertainty Est |
          |  - Ràng buộc vật lý trực tiếp                         |
          |  - Bù trừ bất định phi tuyến & nhiễu sóng cấp 5       |
          |  - Tối ưu hóa QP thời gian thực trên HIL              |
          +-------------------------------------------------------+

Dòng nghiên cứu thứ ba tập trung vào Điều khiển dự báo mô hình (MPC). Zheng et al. (2014) và Belleter et al. (2016) sử dụng NMPC để xử lý ràng buộc trạng thái và điều khiển. Tuy nhiên, nhược điểm chí tử của NMPC là sự phức tạp tính toán phi tuyến (Non-convex Optimization), tiềm ẩn nguy cơ rơi vào cực tiểu địa phương và không khả thi khi triển khai trên phần cứng nhúng thời gian thực.

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

  1. Nghiên cứu của Li et al. (2018) trên tạp chí IEEE Transactions on Industrial Electronics: Sử dụng NMPC kết hợp thuật toán tối ưu bầy đàn (PSO) để điều khiển tàu ASV. So với Li et al., luận án tạo ra bước đột phá khi không sử dụng thuật toán tiến hóa ngẫu nhiên phức tạp mà áp dụng kỹ thuật tuyến tính hóa từng đoạn mô hình song tuyến dọc trục thời gian, biến đổi thành bài toán Quy hoạch toàn phương lặp (Sequential Quadratic Programming - SQP). Giải pháp này rút ngắn thời gian tính toán hơn 4 lần trên cùng cấu hình vi xử lý.
  2. Nghiên cứu của Sonnenburg và Woolsey (2013): Áp dụng phương pháp điều khiển phân tách với bộ quan sát trạng thái phi tuyến. Tuy nhiên, công trình này giả định các ma trận khối lượng kèm $M_A$ và ma trận giảm chấn $D(\nu)$ hoàn toàn tiền định. Luận án của NCS. Nguyễn Hữu Quyền vượt trội hơn khi xây dựng cấu trúc bù bất định hàm $\Delta(t, \nu)$, cho phép hệ thống duy trì độ chính xác bám quỹ đạo ngay cả khi các hệ số thủy động lực học $X_{ui}, Y_{vr}, N_{vr}$ biến thiên bất định trên $30%$.

Đó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 lý thuyết điều khiển phi tuyến ứng dụng trong hàng hải thông qua ba đóng góp lý thuyết nền tảng:

Thứ nhất, luận án đóng góp vào Lý thuyết Điều khiển Dự báo (Model Predictive Control Theory) bằng cách thiết lập thuật toán tuyến tính hóa từng đoạn mô hình song tuyến (Bilinear Model) dọc theo quỹ đạo dự báo. Bản chất mô hình động lực học tàu 3-DOF chứa các tích chéo giữa vận tốc dài và vận tốc góc (thành phần Coriolis $C(\nu)\nu$). Luận án chứng minh rằng việc rời rạc hóa và đóng băng ma trận hệ thống tại từng điểm làm việc $k$ dọc theo chân trời dự báo $N_p$ cho phép tái cấu trúc toàn bộ mô hình phi tuyến thành hệ tuyến tính biến đổi theo thời gian (LTV). Điều này mang lại sự chuyển dịch mô thức (paradigm shift) từ tối ưu phi tuyến phức tạp sang bài toán tối ưu lồi QP có nghiệm tường minh và duy nhất.

Thứ hai, công trình mở rộng Lý thuyết Quan sát và Lọc tối ưu (State Estimation & Filtering) qua việc dẫn xuất thành công Bộ quan sát trạng thái trực tiếp từ mô hình liên tục (Continuous-time Direct State Observer - QSTT) và kết hợp với Bộ lọc Kalman mở rộng (EKF) để khôi phục toàn vẹn véc-tơ trạng thái $[\hat{u}, \hat{v}, \hat{r}, \hat{x}, \hat{y}, \hat{\psi}]^T$ từ tín hiệu đầu ra đo đạc bị nhiễu bởi hệ thống định vị GPS và La bàn điện (Gyrocompass).

Thứ ba, luận án thiết lập Mô hình bù bất định hàm (Lumped Uncertainty Compensation Formulation). Bằng cách gộp toàn bộ các thành phần phi tuyến khó mô hình hóa bao gồm lực quán tính khối lượng nước kèm không đối xứng, ma trận suy giảm phi tuyến bậc cao $D_n(\nu)$ theo chuẩn SNAME (1950): $$D_n(\nu) = \text{diag}\left(\sum_{i=2}^3 X_{ui}|u|^{i-1}, \sum_{i=2}^3 Y_{vi}|v|^{i-1}, \sum_{i=2}^3 N_{ri}|r|^{i-1}\right)$$ cùng các lực nhiễu môi trường đại dương $\tau_w$ thành một véc-tơ bất định duy nhất $\Delta(t, \nu) = [\Delta_u, \Delta_v, \Delta_r]^T$, tác giả đã thiết kế thuật toán ước lượng trực tuyến $\hat{d}(t)$ phản hồi về đầu vào điều khiển, đưa đối tượng thực về tương đương với một mô hình danh định xác định.

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

Khung phân tích của luận án tích hợp ba trụ cột lý thuyết: (1) Lý thuyết Thủy động lực học chuyển động tàu (Fossen, 2011), (2) Lý thuyết Tối ưu hóa có ràng buộc trượt dọc trục thời gian RHC (Maciejowski, 2002), và (3) Lý thuyết Điều khiển học lặp (Iterative Learning Control - ILC).

+---------------------------------------------------------------------------------------------------+
|                                KHUNG PHÂN TÍCH VÀ ĐIỀU KHIỂN DMPC-O                               |
+---------------------------------------------------------------------------------------------------+
                                                                                                    
   Quỹ đạo đặt                                                                                      
   eta_d = [x_d, y_d, psi_d]^T                                                                      
         |                                                                                          
         v        e(k)      +--------------------+   tau(k)     +--------------------+   tau_real   
       ( + ) -------------> |   BỘ TỐI ƯU QP     | -----------> | BỘ BÙ BẤT ĐỊNH     | ----------+  
         ^ -                |   (Ràng buộc a, b) |              | tau = tau_0 - d_hat|           |  
         |                  +--------------------+              +--------------------+           |  
         |                            ^                                   ^                      v  
         |                            |                                   |             +--------------+
         |                  +--------------------+              +--------------------+  |  TÀU THỦY    |
         |                  |  MÔ HÌNH DỰ BÁO    |              |  BỘ ƯỚC LƯỢNG      |  |  3-DOF       |
         |                  |  TUYẾN TÍNH HÓA    |              |  BẤT ĐỊNH d_hat    |  | (Underact.)  |
         |                  +--------------------+              +--------------------+  +--------------+
         |                            ^                                   ^                      |  
         |                            |                                   |                      | y(t)
         |                            +-----------------+-----------------+                      v  
         |                                              |                               +--------------+
         +----------------------------------------------+------------------------------ |   BỘ LỌC     |
                   Vec-to trang thai uoc luong: x_hat = [u, v, r, x, y, psi]^T           |   EKF/QSTT   |
                                                                                        +--------------+

Đặc tính độc đáo của khung phân tích nằm ở Cấu trúc tách rời hai giai đoạn (Two-stage Decoupled Architecture):

  • Giai đoạn 1 (Fast-loop): Bộ ước lượng bất định thực hiện triệt tiêu tức thời ảnh hưởng của nhiễu động lực học $\Delta(t, \nu)$ tại tầng đầu vào chấp hành.
  • Giai đoạn 2 (Optimization-loop): Bộ điều khiển MPC tính toán lời giải tối ưu cho bài toán quy hoạch toàn phương trên mô hình tuyến tính danh định với hàm mục tiêu trượt: $$J(k) = \sum_{j=1}^{N_p} |\hat{\eta}(k+j|k) - \eta_d(k+j)|Q^2 + \sum{j=0}^{N_c-1} |\Delta \tau(k+j|k)|R^2$$ thỏa mãn điều kiện ràng buộc biên độ lực và góc lái: $\tau{\min} \le \tau(k+j) \le \tau_{\max}$, $\Delta \tau_{\min} \le \Delta \tau(k+j) \le \Delta \tau_{\max}$.

Điều kiện biên (Boundary Conditions): Khung phân tích bảo đảm tính hội tụ khi vận tốc trượt dọc $u(t) \ge u_{\min} > 0$ (loại trừ điểm kỳ dị động học của tàu thiếu cơ cấu chấp hành) và ma trận trọng số $Q \ge 0, R > 0$.

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

Thiết kế nghiên cứu

Nghiên cứu được định vị vững chắc trên hệ hình Thực chứng luận sau hiện đại (Post-positivism / Critical Realism), kết hợp chặt chẽ giữa diễn dịch toán học thuần túy (Deductive Mathematical Derivations) và kiểm chứng thực nghiệm thực tế. Thiết kế nghiên cứu gồm 4 giai đoạn logic chặt chẽ:

  1. Giai đoạn thiết lập mô hình giải tích 6-DOF và rút gọn về 3-DOF phi tuyến có bất định.
  2. Giai đoạn tổng hợp thuật toán điều khiển dự báo phân lớp: MPC-S (State feedback), MPC-O (Output feedback), DMPC-S (Disturbance state feedback) và DMPC-O (Disturbance output feedback).
  3. Giai đoạn chứng minh toán học tính ổn định tiệm cận của hệ thống kín sử dụng tiêu chuẩn ổn định Lyapunov cho hệ điều khiển dự báo có ràng buộc.
  4. Giai đoạn mô phỏng số học (Numerical Simulation) và kiểm chứng thực nghiệm bán vật lý trên mô hình Phần cứng trong vòng lặp (Hardware-In-the-Loop - HIL).

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu áp dụng kỹ thuật tam giác hóa phương pháp luận (Methodological Triangulation):

  • Triangulation lý thuyết: So sánh đối chiếu kết quả thiết kế giữa mô hình giải tích phi tuyến liên tục, mô hình song tuyến rời rạc hóa và mô hình xấp xỉ không gian trạng thái.
  • Triangulation công cụ kiểm chứng: Đối chuẩn thuật toán đồng thời trên môi trường mô phỏng thuần túy MATLAB/Simulink, bộ công cụ mô phỏng chuẩn quốc tế MSS-GNC (Marine Systems Simulator - Guidance Navigation and Control Toolbox) và hệ thống thực nghiệm thời gian thực HIL.
+-------------------------------------------------------------------------------+
|                      CẤU HÌNH THỰC NGHIỆM HIL (HARDWARE-IN-THE-LOOP)          |
+-------------------------------------------------------------------------------+
|                                                                               |
|  +---------------------------+                      +----------------------+  |
|  |     MÁY TÍNH PC1          |                      |     MÁY TÍNH PC2     |  |
|  | (Thuật toán MPC/DMPC-O)   |                      | (Mô phỏng động lực   |  |
|  |                           |                      |  học tàu MSS-GNC)    |  |
|  +---------------------------+                      +----------------------+  |
|          |         ^                                        |         ^       |
|   Lực    |         | Tín hiệu đo                            | Trạng   | Lực   |
|   điều   |         | GPS/Gyro                               | thái    | điều  |
|   khiển  |         | (NMEA-0183)                            | thực    | khiển |
|   tau    |         |                                        | y(t)    | tau   |
|          v         |                                        v         |       |
|  +-----------------------------------------------------------------+  |
|  |         CARD GHÉP NỐI VẬT LÝ ARDUINO DUE R3 (ATmega16U2)        |  |
|  |      - Chuyển đổi giao thức chuẩn Hàng hải NMEA-0183            |  |
|  |      - Xử lý truyền thông UART/Serial thời gian thực            |  |
|  |      - Tạo lập nhiễu đo trắng Gauss (Noise Injection)           |  |
|  +-----------------------------------------------------------------+  |
|                                                                               |
+-------------------------------------------------------------------------------+

Quy trình thiết kế bộ lọc EKF được bảo đảm tính hợp thức thông qua việc xác định chính xác ma trận hiệp phương sai nhiễu quá trình $Q_k$ và ma trận hiệp phương sai nhiễu đo $R_k$, loại bỏ hoàn toàn hiện tượng phân kỳ do sai số làm tròn số học.

Data và phân tích

Môi trường thực nghiệm và thông số đối tượng: Mô hình thực nghiệm sử dụng thông số động lực học chuẩn của tàu chở hàng choán nước trong thư viện MSS-GNC Toolbox, với các thông số vật lý danh định:

  • Chiều dài tàu $L = 32.0\text{ m}$, khối lượng $m = 23.8 \times 10^4\text{ kg}$, mô-men quán tính trục đứng $I_z = 1.96 \times 10^7\text{ kg}\cdot\text{m}^2$.
  • Ma trận quán tính khối lượng kèm: $X_{\dot{u}} = -2.0 \times 10^4\text{ kg}$, $Y_{\dot{v}} = -1.0 \times 10^5\text{ kg}$, $N_{\dot{r}} = -8.0 \times 10^6\text{ kg}\cdot\text{m}^2$.
  • Ràng buộc tín hiệu điều khiển vật lý: Lực đẩy trượt dọc $\tau_u \in [0, 1.5 \times 10^5]\text{ N}$, mô-men bẻ lái $\tau_r \in [-3.0 \times 10^5, +3.0 \times 10^5]\text{ N}\cdot\text{m}$, tương ứng góc bẻ lái $\delta \in [-35^\circ, +35^\circ]$.

Nền tảng thực nghiệm HIL (Hardware-In-The-Loop):

  • Phần cứng: Vi điều khiển Arduino Due R3 (32-bit ARM Cortex-M3 SAM3X8E xung nhịp 84 MHz kết hợp chip truyền thông ATmega16U2).
  • Giao thức: Chuẩn truyền tin điện tử Hàng hải Quốc tế NMEA-0183 mô phỏng câu bản tin $GPGGA (định vị vệ tinh GPS) và $HEHDT (la bàn điện Gyrocompass).
  • Phần mềm: Hệ thống phân tán trên 2 máy trạm (PC1 thực thi bộ điều khiển MPC trong Simulink Real-Time, PC2 mô phỏng môi trường biển 3-DOF và tương tác sóng gió). Chu kỳ trích mẫu điều khiển cố định $T_s = 0.1\text{ s}$, chân trời dự báo $N_p = 10$, chân trời điều khiển $N_c = 3$.

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

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

Công trình đã chứng minh và lượng hóa 4 phát hiện đột phá thông qua thực nghiệm HIL:

+-------------------------------------------------------------------------------+
|            BẢNG ĐỐI CHUẨN HIỆU NĂNG BÁM QUỸ ĐẠO VÀ KHÁNG NHIỄU BIỂN           |
+------------------------------------+------------------+-----------------------+
| Tiêu chí Đánh giá                  | MPC Kinh điển    | Đề xuất (DMPC-O)      |
+------------------------------------+------------------+-----------------------+
| Sai lệch vị trí xác lập (Vòng tròn)| 3.82 m           | 0.41 m (Giảm 89.3%)   |
| Sai lệch vị trí xác lập (Hình Sin) | 4.15 m           | 0.53 m (Giảm 87.2%)   |
| Biên độ góc bẻ lái cực đại delta   | 38.6 độ (Vi phạm)| 28.4 độ (Thỏa mãn <35)|
| Thời gian tính toán 1 chu kỳ lặp   | 340 ms (NMPC)    | 42 ms (Giảm 87.6%)    |
| Khả năng kháng nhiễu biển cấp 5    | Phân kỳ nhẹ      | Ổn định tiệm cận      |
+------------------------------------+------------------+-----------------------+
  1. Triệt tiêu sai lệch bám quỹ đạo vượt trội: Dưới tác động của mô hình bất định và nhiễu sóng gió cấp 5, thuật toán DMPC-O duy trì sai lệch bám vị trí $(e_x, e_y)$ trên quỹ đạo cong hình tròn ($R = 100\text{ m}$) và quỹ đạo hình sin tuần hoàn dưới $0.53\text{ m}$, giảm hơn $87%$ sai số so với các bộ điều khiển MPC tuyến tính tiêu chuẩn không có khối bù bất định (sai số $> 4.1\text{ m}$).
  2. Tuân thủ tuyệt đối ràng buộc cơ cấu chấp hành: Tín hiệu góc bẻ lái $\delta(t)$ và lực đẩy $\tau_u(t)$ luôn nằm hoàn toàn trong miền ràng buộc kỹ thuật $[-35^\circ, +35^\circ]$ mà không xảy ra hiện tượng bão hòa tín hiệu hay sốc cơ học, loại bỏ triệt để hiện tượng rung giật bánh lái thường thấy ở bộ điều khiển trượt SMC.
  3. Hiện tượng tách ghép thích ứng bất định (Counter-intuitive Decoupling): Dữ liệu thực nghiệm chỉ ra rằng, khi tàu chuyển hướng gắt trên quỹ đạo sin, lực cản trượt ngang $v$ tăng vọt phi tuyến. Bộ ước lượng $\hat{d}(t)$ đã phản ứng bù trừ chính xác mô-men bất định $\Delta_r$ trong vòng $0.2\text{ s}$, cho phép tàu giữ vững hướng mũi mà không cần tăng đột biến công suất máy chính.
  4. Hiệu năng tính toán thời gian thực: Kỹ thuật tuyến tính hóa từng đoạn mô hình song tuyến đã giảm thời gian thực thi thuật toán tối ưu từ $340\text{ ms}$ (đối với NMPC giải bằng fmincon) xuống còn trung bình $42\text{ ms}$ trên mỗi chu kỳ lấy mẫu, hoàn toàn đáp ứng yêu cầu khắt khe của hệ thống nhúng điều khiển tàu thời gian thực ($T_s = 100\text{ ms}$).

Implications đa chiều

  • Về lý thuyết: Khẳng định tính khả thi của việc tích hợp nguyên lý trượt thời gian RHC với kỹ thuật bù bất định hàm, cung cấp một khung phương pháp luận mới cho lớp đối tượng phi tuyến thiếu cơ cấu chấp hành có ràng buộc.
  • Về phương pháp luận: Mô hình thực nghiệm HIL phân tán thông qua vi điều khiển Arduino Due và giao thức NMEA-0183 chuẩn hóa quy trình kiểm thử trước khi lắp đặt thực tế trên tàu, giảm thiểu rủi ro và chi phí thử nghiệm đường biển.
  • Về thực tiễn ứng dụng: Cung cấp giải pháp sẵn sàng triển khai cho các nhà máy đóng tàu và viện nghiên cứu hàng hải để nâng cấp hệ thống Lái tự động (Autopilot) thông thường thành hệ thống Tự hành bám quỹ đạo thông minh trên các dòng tàu kiểm ngư, tàu chở hàng tải trọng lớn và phương tiện mặt nước không người lái USV.
  • Về chính sách: Đóng góp cơ sở khoa học và bằng chứng thực nghiệm vững chắc phục vụ việc xây dựng quy chuẩn kỹ thuật quốc gia về thử nghiệm phương tiện hàng hải tự hành trong vùng nước cảng biển Việt Nam.

Limitations và Future Research

Luận án thừa nhận một số giới hạn nghiên cứu khách quan:

  1. Giới hạn bậc tự do: Nghiên cứu mới tập trung giải quyết bài toán chuyển động trong mặt phẳng ngang 3-DOF (surge, sway, yaw), chưa xét đến động học ghép kênh đầy đủ 6-DOF bao gồm lắc ngang (roll), lắc dọc (pitch) và trượt đứng (heave).
  2. Giới hạn cường độ môi trường: Các kiểm chứng thực nghiệm mới đánh giá trong điều kiện sóng gió không vượt quá cấp 5 Beaufort; chưa phân tích trạng thái giới hạn khi gặp sóng biển cấp bão (cấp 8-9) dẫn đến hiện tượng trượt chân vịt hoặc bánh lái nhô khỏi mặt nước.
  3. Đặc tính cơ cấu chấp hành: Mô hình nghiên cứu chưa tích hợp sâu phương trình động lực học trễ của hệ thống thủy lực máy lái (steering gear dynamics) và vùng chết của chân vịt.

Chương trình nghiên cứu 10 năm tiếp theo tập trung vào các hướng mở rộng:

  • Mở rộng thuật toán DMPC-O cho mô hình chuyển động không gian 6-DOF, tích hợp bài toán giảm lắc ngang (Roll Stabilization) đồng thời với bám quỹ đạo.
  • Nghiên cứu cơ chế MPC phân tán (Distributed MPC) ứng dụng cho bài toán điều khiển đội hình (flocking/formation control) của cụm nhiều tàu tự hành USV phối hợp tuần tra.
  • Ứng dụng mạng học sâu tăng cường (Deep Reinforcement Learning) để tối ưu hóa trực tuyến các ma trận trọng số $Q$ và $R$ trong bộ điều khiển MPC khi môi trường biển biến đổi phức tạp.
  • Triển khai thử nghiệm thực địa quy mô đầy đủ (Full-scale sea trials) trên tàu vận tải biển thực tế tại vịnh Bắc Bộ.

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

  • Ảnh hưởng học thuật: Công trình mở ra hướng nghiên cứu liên ngành giữa lý thuyết điều khiển tối ưu hiện đại và kỹ thuật hàng hải; các công bố liên quan đóng góp trực tiếp vào hệ thống tài liệu tham khảo cho các chương trình đào tạo tiến sĩ ngành Kỹ thuật Điều khiển và Tự động hóa.
  • Chuyển đổi công nghiệp hàng hải: Cung cấp lõi công nghệ điều khiển bám quỹ đạo chính xác cho các doanh nghiệp đóng tàu trong nước (như Tổng công ty Công nghiệp Tàu thủy - SBIC), hỗ trợ tăng hàm lượng nội địa hóa công nghệ điều khiển trên các tàu chuyên dụng cảnh sát biển và vận tải hàng rời.
  • Lợi ích kinh tế - xã hội: Tối ưu hóa quỹ đạo chuyển động của tàu giúp rút ngắn chiều dài hành trình thực tế từ $3 - 5%$, giảm thiểu tiêu hao nhiên liệu máy chính, giảm phát thải khí nhà kính ($CO_2, NO_x$) ra môi trường biển, đồng thời hạn chế tối đa nguy cơ đâm va hàng hải do lỗi con người.

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

  • Nghiên cứu sinh và học viên cao học: Tiếp cận phương pháp luận thiết kế hệ thống điều khiển dự báo cho đối tượng phi tuyến thiếu cơ cấu chấp hành và quy trình xây dựng giàn thử nghiệm HIL thời gian thực.
  • Các nhà khoa học đầu ngành: Tài liệu tham khảo giá trị trong việc phát triển các thuật toán điều khiển tối ưu thích nghi, điều khiển bền vững cho các phương tiện không người lái dưới nước (AUV) và trên không (UAV).
  • Bộ phận R&D tại các doanh nghiệp hàng hải: Cơ sở công nghệ để thiết kế, chế tạo các thế hệ máy lái tự động thông minh tích hợp tính năng bám tuyến đường (Track-keeping Autopilot).
  • Cơ quan quản lý nhà nước (Cục Hàng hải Việt Nam): Nguồn dữ liệu khoa học tin cậy phục vụ công tác xây dựng hành lang pháp lý, tiêu chuẩn an toàn cho tàu tự hành hoạt động trong luồng hàng hải nội địa.

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? Trả lời: Đóng góp độc đáo nhất là việc phát triển thành công giải pháp điều khiển dự báo MPC hệ song tuyến trên cơ sở tuyến tính hóa từng đoạn mô hình phi tuyến dọc trục thời gian, tích hợp bộ bù bất định hàm đầu vào cho đối tượng tàu thiếu cơ cấu chấp hành. Công trình đã mở rộng Lý thuyết Điều khiển Dự báo Mô hình (MPC) của Camacho & Bordons (2004) và Maciejowski (2002) sang lớp phương tiện phi tuyến cao bậc có bất định tham số thủy động lực học, giải quyết đồng thời bài toán tối ưu hóa có ràng buộc và bài toán khử bất định trong một cấu trúc điều khiển phân tầng đồng nhất.

  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ế tiền nhiệm? Trả lời: So với công trình của Do et al. (2004) (sử dụng Backstepping thuần túy không xử lý được ràng buộc góc lái) và công trình của Zheng et al. (2014) (sử dụng NMPC phi tuyến nguyên bản với thời gian tính toán chu kỳ vượt quá $300\text{ ms}$), phương pháp luận của luận án mang tính đột phá khi chuyển đổi bài toán tối ưu phi tuyến thành bài toán Quy hoạch toàn phương (QP) liên tục. Kết hợp với việc xây dựng mô hình thực nghiệm HIL ghép nối phần cứng vi điều khiển Arduino Due thời gian thực giao tiếp chuẩn NMEA-0183, nghiên cứu đã lấp đầy khoảng cách giữa mô phỏng số học lý thuyết và khả năng thực thi phần cứng thực tế.

  3. Phát hiện bất ngờ và thú vị nhất thu được từ dữ liệu thực nghiệm là gì? Trả lời: Phát hiện bất ngờ nhất là khi tàu hoạt động trên quỹ đạo cong có bán kính nhỏ, lực cản trượt ngang phi tuyến $v$ tăng vọt dẫn tới mô-men cản quay trở lớn; tuy nhiên, thay vì làm tăng góc bẻ lái $\delta$ đến mức bão hòa biên độ $35^\circ$ như các bộ điều khiển thông thường, bộ điều khiển DMPC-O tự động điều chỉnh tốc độ đặt dọc $u_d$ thông qua cơ chế học lặp, giúp con tàu tự hãm tốc độ tiến một cách nhịp nhàng để hoàn thành vòng quay với góc lái duy trì ổn định dưới $28.4^\circ$, bảo đảm tuyệt đối an toàn hàng hải.

  4. Luận án có cung cấp đầy đủ giao thức để tái lập thực nghiệm (Replication Protocol) hay không? Trả lời: Luận án cung cấp tường minh toàn bộ giao thức tái lập bao gồm: hệ phương trình vi phân mô hình toán 3-DOF với đầy đủ hệ số ma trận khối lượng $M$, ma trận Coriolis $C(\nu)$ và ma trận suy giảm $D(\nu)$; thuật toán chi tiết từng bước của bộ quan sát EKF và bộ tối ưu QP; sơ đồ nguyên lý đấu nối chân phần cứng Arduino Due R3 với máy tính PC qua cổng UART; cấu hình bản tin NMEA-0183 và bộ tham số cài đặt của thư viện MSS-GNC Toolbox trong môi trường MATLAB/Simulink.

  5. Chương trình nghiên cứu 10 năm tới được phác thảo như thế nào? Trả lời: Lộ trình 10 năm được xác lập qua 3 giai đoạn: (1) 2020-2023: Hoàn thiện thuật toán điều khiển thích ứng 6-DOF và tích hợp mô hình cơ cấu chấp hành thủy lực phi tuyến; (2) 2024-2027: Phát triển giải pháp điều khiển phân tán cho cụm tàu tự hành USV phối hợp tránh va tự động theo quy tắc COLREGs; (3) 2028-2030: Chuyển giao công nghệ, chế tạo thử nghiệm bo mạch điều khiển công nghiệp chuyên dụng và thực hiện thử nghiệm đường biển quy mô đầy đủ trên tàu kiểm ngư.

Kết luận

Luận án của NCS. Nguyễn Hữu Quyền đã giải quyết trọn vẹn và xuất sắc mục tiêu nghiên cứu đặt ra, khẳng định những đóng góp khoa học và thực tiễn mang tính then chốt:

  1. Xây dựng và chuẩn hóa hệ thống mô hình toán học động lực học 3 bậc tự do mô tả chuyển động tàu thủy thiếu cơ cấu chấp hành có xét đến đầy đủ các thành phần bất định mô hình và nhiễu sóng gió môi trường biển.
  2. Đề xuất thành công cấu trúc điều khiển dự báo MPC tuyến tính hóa từng đoạn mô hình song tuyến dọc trục thời gian, giải quyết triệt để bài toán tối ưu hóa bám quỹ đạo có ràng buộc tín hiệu điều khiển đầu vào với chi phí tính toán thấp.
  3. Thiết kế thành công Bộ quan sát trạng thái trực tiếp (QSTT) và Bộ lọc Kalman mở rộng (EKF) khôi phục chính xác trạng thái chuyển động của tàu từ các cảm biến hàng hải bị nhiễu.
  4. Sáng tạo thuật toán ước lượng và bù trực tuyến thành phần bất định hàm đầu vào $\Delta(t, \nu)$, nâng cao năng lực kháng nhiễu và độ bền vững của hệ thống điều khiển dưới tác động của môi trường biển lên tới cấp 5 Beaufort.
  5. Thiết lập hoàn chỉnh giàn thử nghiệm bán vật lý Hardware-In-the-Loop (HIL) trên nền tảng vi điều khiển nhúng Arduino Due và chuẩn giao tiếp hàng hải NMEA-0183, chứng minh tính khả thi ứng dụng công nghiệp của thuật toán.
  6. Mở ra ba dòng nghiên cứu học thuật mới đầy triển vọng: Điều khiển dự báo phân tán cho hệ đa phương tiện nổi tự hành, điều khiển bám quỹ đạo tích hợp giảm lắc 6-DOF, và ứng dụng trí tuệ nhân tạo tối ưu hóa thông số điều khiển dự báo thích ứng trong ngành tự động hóa hàng hải Việt Nam và quốc tế.