Báo Cáo Khoa Học: Nghiên Cứu Điều Khiển Trượt Tối Ưu Bám Quỹ Đạo Cho Hệ Pendubot Tại HCMUTE

Tóm tắt nghiên cứu (Research Summary)

  • Câu hỏi nghiên cứu chính (Research Question): Làm thế nào để thiết kế một giải thuật điều khiển phi tuyến bền vững, có khả năng tối ưu hóa các tham số mặt trượt nhằm điều khiển bám quỹ đạo chính xác cho hệ con lắc ngược nhiều khâu (Pendubot 3-link) – một hệ thống cơ điện tử phi tuyến thiếu cơ cấu chấp hành (underactuated SIMO) bậc cao phức tạp?
  • Tổng quan phương pháp tiếp cận (Methodology Snapshot): Đề tài ứng dụng phương pháp cơ học giải tích Euler-Lagrange để xây dựng mô hình toán học giải tích toàn diện cho hệ Pendubot 3 khâu. Trên cơ sở đó, nhóm nghiên cứu phát triển giải thuật Điều khiển trượt tối ưu (Optimal Sliding Mode Control - OSMC) kết hợp với hàm mục tiêu tối ưu ma trận trạng thái sai số (LQR) và tiêu chuẩn ổn định Lyapunov; đồng thời ứng dụng Giải thuật Di truyền (Genetic Algorithm - GA) để tự động hóa quá trình dò tìm bộ thông số điều khiển tối ưu toàn cục.
  • Phát hiện cốt lõi (Key Findings): Giải thuật trượt tối ưu đề xuất đã giải quyết triệt để điểm nghẽn toán học của phương pháp trượt thứ bậc (Hierarchical Sliding Mode Control) truyền thống khi chứng minh được mọi sai số trạng thái đều hội tụ tiệm cận về 0. Kết quả mô phỏng trên nền tảng MATLAB/Simulink cho thấy hệ thống bám sát quỹ đạo sóng sin mong muốn chỉ sau 2 giây, đạt trạng thái xác lập toàn phần trong 5 giây với sai số xác lập xấp xỉ bằng 0 và biên độ dao động của các khâu thụ động được triệt tiêu tối đa.
  • Hàm ý khoa học và ứng dụng (Implications): Kết quả nghiên cứu tạo nền tảng lý thuyết và thuật toán vững chắc để mở rộng điều khiển cho các hệ thống robot công nghiệp nhiều bậc tự do, tay máy cộng tác (cobots), cũng như các hệ thống cân bằng động phức tạp trong công nghiệp 4.0.

Bối cảnh và tầm quan trọng (Context & Significance)

Tổng quan hiện trạng nghiên cứu (Current State of Knowledge)

Trong lý thuyết điều khiển tự động và cơ điện tử, hệ con lắc ngược nhiều khâu (Pendubot) đại diện cho lớp đối tượng kinh điển thuộc nhóm SIMO (Single-Input Multiple-Output) thiếu cơ cấu chấp hành (underactuated). Với đặc tính phi tuyến mạnh, khớp nối động học phức tạp và số lượng biến trạng thái cần điều khiển vượt quá số lượng tín hiệu tác động đầu vào (chỉ có 1 động cơ dẫn động tại khâu đầu tiên), Pendubot luôn là mô hình chuẩn mực để kiểm nghiệm tính hiệu quả của các giải thuật điều khiển phi tuyến hiện đại.

Phần lớn các công trình nghiên cứu trong và ngoài nước trước đây chủ yếu tập trung vào mô hình Pendubot 2 khâu (2-link) và dừng lại ở bài toán điều khiển cân bằng tại chỗ quanh các điểm làm việc tĩnh (vị trí TOP hoặc MID) bằng các bộ điều khiển tuyến tính như PID, LQR, hoặc điều khiển trượt kinh điển.

                  [Khâu 3 (Link 3 - Thụ động)]
                             o (Khớp 3 - Passive)
                            /
                           /
            [Khâu 2 (Link 2 - Thụ động)]
                       o (Khớp 2 - Passive)
                      /
                     /
      [Khâu 1 (Link 1 - Chủ động)]
                 |
                 O (Khớp 1 - Động cơ cấp Momen Tm)

Khoảng trống nghiên cứu (Research Gap)

Mặc dù đã có nhiều nghiên cứu về điều khiển hệ SIMO, hai khoảng trống học thuật lớn vẫn tồn tại:

  1. Hạn chế của điều khiển bám quỹ đạo trên hệ bậc cao: Các bộ điều khiển tuyến tính (LQR, PID) chỉ làm việc quanh vùng lân cận điểm cân bằng và nhanh chóng mất tính ổn định khi biên độ quỹ đạo thay đổi rộng. Trong khi đó, việc mở rộng từ hệ 2 khâu lên hệ 3 khâu (3-link) làm tăng bậc của hệ thống lên 6 trạng thái phi tuyến cao, khiến việc tổng hợp luật điều khiển bám quỹ đạo động học trở nên cực kỳ thách thức.
  2. Nhược điểm của giải thuật trượt thứ bậc (Hierarchical Sliding Mode Control - HSMC): Mặc dù HSMC được thiết kế cho hệ SIMO, phương pháp này chỉ chứng minh được các hàm mặt trượt tổng hội tụ về 0, nhưng không chứng minh được từng thành phần sai số trạng thái $e_i$ triệt tiêu về 0. Hơn nữa, việc lựa chọn vector tham số mặt trượt $c_i$ phần lớn dựa vào phương pháp thử sai (trial-and-error), thiếu cơ sở toán học tường minh.

Tính cấp thiết và tác động kỳ vọng (Timeliness & Potential Impact)

Trong bối cảnh tự động hóa và công nghiệp 4.0, nhu cầu điều khiển chính xác các tay máy robot linh hoạt đòi hỏi những giải thuật phi tuyến vừa đảm bảo tính bền vững (robustness), vừa có cấu trúc toán học chặt chẽ để tiết kiệm tài nguyên tính toán nhúng. Đề tài NCKH (Mã số: T2019-04NCS) do NCS. Huỳnh Xuân Dũng và nhóm nghiên cứu tại Trường Đại học Sư phạm Kỹ thuật TP.HCM (HCMUTE) thực hiện đã giải quyết trực diện bài toán này, đem lại giá trị học thuật cao và bổ sung công cụ mô phỏng chuẩn mực cho công tác đào tạo, nghiên cứu chuyên sâu.

Phương pháp luận và cách tiếp cận (Methodology & Approach)

+-------------------------------------------------------------------------+
|                  1. Mô hình hóa Euler-Lagrange                          |
|  Xây dựng hệ phương trình động lực học phi tuyến 6 trạng thái (3-link)  |
+------------------------------------+------------------------------------+
                                     |
                                     v
+-------------------------------------------------------------------------+
|             2. Thiết kế Mặt trượt Tối ưu (OSMC Design)                  |
|  Thiết lập mặt trượt nối tầng s1..s5 + Đồng nhất tham số qua ma trận K  |
+------------------------------------+------------------------------------+
                                     |
                                     v
+-------------------------------------------------------------------------+
|            3. Chứng minh Tính Ổn định Tiệm cận (Lyapunov)               |
|      Xác định tín hiệu điều khiển U = Ueq + Usw thỏa mãn dV5/dt <= 0    |
+------------------------------------+------------------------------------+
                                     |
                                     v
+-------------------------------------------------------------------------+
|        4. Tối ưu hóa Tham số Toàn cục bằng Giải thuật Di truyền (GA)     |
|   Tìm nghiệm ma trận Q, L và hệ số k, eta triệt tiêu sai số bám J(e)    |
+------------------------------------+------------------------------------+
                                     |
                                     v
+-------------------------------------------------------------------------+
|            5. Kiểm chứng Mô phỏng Động học trên MATLAB/Simulink         |
|         Thời gian lấy mẫu 10ms, quỹ đạo hình sin x1d = sin(2*pi*t/30)   |
+-------------------------------------------------------------------------+

Thiết kế nghiên cứu và mô hình hóa toán học (Research Design & Modeling)

Nghiên cứu áp dụng quy trình kết hợp giữa phân tích giải tích giải thuật và thực nghiệm mô phỏng số. Phương trình động lực học tổng quát của hệ Pendubot 3 khâu được thiết lập thông qua phương pháp Euler-Lagrange:

$$\frac{d}{dt}\left(\frac{\partial L}{\partial \dot{q}_k}\right) - \frac{\partial L}{\partial q_k} = \tau_k \quad (k=1, 2, 3)$$

Trong đó $L(q, \dot{q}) = T(q, \dot{q}) - V(q)$ là toán tử Lagrange (hiệu giữa tổng động năng $T$ và thế năng $V$ của hệ 3 thanh liên kết). Hệ phương trình vi phân phi tuyến bậc hai được đưa về dạng ma trận không gian trạng thái:

$$D(q)\ddot{q} + C(q, \dot{q})\dot{q} + G(q) = F T_m$$

Bằng kỹ thuật đặt lại biến trạng thái quanh vị trí thẳng đứng cân bằng (TOP), hệ thống được hạ bậc thành hệ 6 phương trình vi phân phi tuyến bậc một:

$$\begin{cases} \dot{x}_1 = x_2 \ \dot{x}_2 = f_1(x) + b_1(x)u \ \dot{x}_3 = x_4 \ \dot{x}_4 = f_2(x) + b_2(x)u \ \dot{x}_5 = x_6 \ \dot{x}_6 = f_3(x) + b_3(x)u \end{cases}$$

Kỹ thuật tổng hợp giải thuật điều khiển trượt tối ưu (OSMC Synthesis)

Để giải quyết bài toán bám quỹ đạo $x_{1d}(t)$ cho Link 1 và giữ Link 2, Link 3 ở vị trí cân bằng ($x_{3d} = x_{5d} = 0$), nhóm tác giả đề xuất cấu trúc mặt trượt đa tầng tối ưu:

$$s_1 = c_1 e_1 + c_2 e_2$$ $$s_2 = s_1 + c_3 e_3$$ $$s_3 = s_2 + c_4 e_4$$ $$s_4 = s_3 + c_5 e_5$$ $$s_5 = s_4 + c_6 e_6$$

Điểm đổi mới căn bản là việc đưa vector mặt trượt về dạng quan hệ không gian trạng thái sai số: $S = Z_r - K Z$. Bằng phương pháp đồng nhất hệ số với lý thuyết điều khiển tối ưu toàn phương tuyến tính (LQR), các hệ số mặt trượt $c_1, c_2, c_3, c_4, c_5, c_6$ được tính toán giải tích trực tiếp từ ma trận độ lợi tối ưu $K$, loại bỏ hoàn toàn yếu tố chọn tham số theo cảm tính.

Phân tích tính ổn định Lyapunov và tối ưu hóa bằng GA

  • Chứng minh ổn định Lyapunov: Luật điều khiển tổng hợp $u = u_{eq} + u_{sw}$ được xây dựng từ hàm Lyapunov cấp cao $V_5 = \frac{1}{2}s_5^2$. Đạo hàm $\dot{V}_5 \le -k |s_5| \le 0$ bảo đảm hệ thống hội tụ tiệm cận toàn cục về mặt trượt $s_5 = 0$, kéo theo toàn bộ các mặt trượt thành phần và vector sai số $e(t) \to 0$.
  • Tối ưu hóa đa tham số bằng GA: Để tìm bộ thông số tối ưu tuyệt đối cho các ma trận trọng số $Q$, $L$ và các hệ số chuyển mạch $k, \eta$, nghiên cứu tích hợp Giải thuật Di truyền với cấu hình: kích thước quần thể 100, số thế hệ 100 (tổng cộng 6000 lượt đánh giá thế hệ tích lũy), mã hóa số thập phân, lai ghép đa điểm (multipoint crossover với hệ số 0.2) và hàm mục tiêu cực tiểu hóa tích phân bình phương sai số:

$$J = \sum_{i=1}^{n} \sum_{k=1}^{6} e_k(i) e_k^T(i)$$

Phát hiện chính và kết quả mô phỏng (Major Findings)

Quá trình mô phỏng kiểm chứng được thực hiện trên môi trường MATLAB/Simulink với chu kỳ lấy mẫu khắt khe $T_s = 10\text{ ms}$, thông số hình học và quán tính thực tế ($m_1=0.055\text{ kg}$, $l_1=l_2=0.2\text{ m}$, $l_3=0.22\text{ m}$), quỹ đạo đặt cho Link 1 là tín hiệu sóng sin điều hòa:

$$x_{1d}(t) = \sin\left(\frac{2\pi t}{30}\right)$$

Góc (rad)
 ^
 |             Quỹ đạo đặt x1d vs Đáp ứng góc Link 1
1+                 _--_                _--_
 |   /----------/        \            /        \
 |  /                      \        /            \
0+--------------------------\------/--------------\---> Thời gian (s)
 |                           \    /
-1+                            --
 |  [Quá độ 0-2s] [ Bám hoàn hảo sau 2s, sai số e1 -> 0 ]

1. Khả năng bám quỹ đạo vượt trội và triệt tiêu sai số xác lập

  • Đáp ứng góc Link 1: Link 1 nhanh chóng chuyển từ trạng thái ban đầu và bắt kịp quỹ đạo sin mong muốn chỉ sau 2.0 giây. Trong suốt khoảng thời gian vận hành tiếp theo, đáp ứng góc của Link 1 trùng khớp hoàn toàn với tín hiệu đặt, sai số bám xấp xỉ bằng 0.
  • Ổn định các khâu thụ động (Link 2 & Link 3): Dưới tác động phi tuyến của khớp chủ động Link 1, các góc lệch $q_2$ và $q_3$ chỉ dao động quá độ nhẹ trong khoảng 0.04 rad trong 5 giây đầu. Sau thời gian xác lập 5s, Link 2 và Link 3 duy trì ổn định tuyệt đối quanh vị trí cân bằng 0 rad với dao động điều hòa biên độ vi mô cùng pha tần số của Link 1.

2. Sự hội tụ đồng thời của các mặt trượt và tính ổn định bền vững

  • Hội tụ mặt trượt: Toàn bộ 5 mặt trượt ($S_1, S_2, S_3, S_4, S_5$) đều suy giảm nhanh chóng và triệt tiêu hoàn toàn về 0 trong vòng 5 giây. Điều này khẳng định tính đúng đắn của phép chứng minh toán học Lyapunov trong việc kéo mọi biến trạng thái sai số về gốc tọa độ.
  • Tối ưu hóa năng lượng momen điều khiển: Tín hiệu momen điều khiển $T_m$ sinh ra từ động cơ biến thiên êm thuận theo dạng sóng sin tương ứng với chuyển động của hệ, không xuất hiện hiện tượng dao động tần số cao biên độ lớn (chattering phá hủy cơ khí) nhờ khâu chuyển mạch trượt $u_{sw}$ được tối ưu hóa biên độ bằng GA.

3. Hiệu quả hội tụ của thuật toán di truyền (GA)

Qua 6000 lượt đánh giá thế hệ, hàm thích nghi $J$ đã giảm từ giá trị phân kỳ ban đầu về giá trị cực tiểu toàn cục $J = 25.0$. Bộ ma trận trọng số tối ưu được xác lập:

$$Q = \text{diag}(452, 366, 479, 315, 25, 113)$$ $$L = \text{diag}(1951, 935, 935, 730, 265, 270)$$

Cùng với các hệ số điều khiển $k = 11.6130$ và $\eta = 2.0$, tạo ra vector khuếch đại $K$ bảo đảm hệ số suy giảm sai số nhanh nhất mà không gây quá tải dòng điện động cơ.

Chỉ số đánh giá Trượt thứ bậc cổ điển (HSMC) Điều khiển tối ưu LQR Trượt tối ưu đề xuất (OSMC + GA)
Chứng minh $e_i \to 0$ Không bảo đảm về toán học Chỉ quanh điểm cân bằng Chứng minh chặt chẽ bằng Lyapunov
Thời gian xác lập ($t_s$) $> 8.0\text{ s}$ (hoặc không hội tụ) Dao động lớn khi biên độ cao $< 5.0\text{ s}$
Sai số bám quỹ đạo ($e_1$) Còn tồn tại sai số tĩnh Lệch lớn khi xa điểm làm việc Xấp xỉ 0 (Bám hoàn hảo sau 2s)
Xác định thông số $c_i$ Thử sai (Trial-and-error) Giải Riccati tuyến tính Tối ưu hóa giải tích kết hợp GA

Đóng góp khoa học (Scientific Contributions)

  1. Đóng góp về mặt lý thuyết điều khiển học:
    • Mở rộng thành công lý thuyết điều khiển trượt cho lớp đối tượng SIMO bậc cao từ 2 khâu lên 3 khâu bám quỹ đạo động học.
    • Xây dựng khung toán học hoàn chỉnh kết hợp giữa hàm Lyapunov nhiều tầng và tối ưu hóa ma trận trạng thái sai số LQR, giải quyết triệt để hạn chế "không thể chứng minh triệt tiêu sai số" tồn tại lâu năm của phương pháp trượt thứ bậc Qian.
  2. Đổi mới về mặt phương pháp luận tính toán:
    • Đề xuất quy trình thiết kế khép kín: từ giải tích Euler-Lagrange $\to$ thiết lập quan hệ mặt trượt không gian sai số $\to$ tự động hóa tìm tham số tối ưu bằng thuật toán di truyền GA, thay thế hoàn toàn kỹ thuật dò thông số thủ công.
  3. Giá trị ứng dụng và đào tạo:
    • Bàn giao bộ công cụ mô phỏng toán học chi tiết trên MATLAB/Simulink cho Phòng thí nghiệm Điều khiển Tự động – Khoa Điện-Điện tử, Trường ĐH Sư phạm Kỹ thuật TP.HCM.
    • Công bố công trình trên tạp chí khoa học chuyên ngành quốc tế (International Journal of Robotica & Management), đóng góp trực tiếp vào mục tiêu đào tạo Nghiên cứu sinh chất lượng cao.

Đối tượng quan tâm và lợi ích ứng dụng (Target Audience & Benefits)

  • Các nhà nghiên cứu học thuật (Academic Researchers): Cung cấp tài liệu tham khảo giải tích mẫu mực về phương pháp mô hình hóa Euler-Lagrange cho hệ con lắc N-khâu và kỹ thuật chứng minh ổn định Lyapunov nhiều tầng cho hệ phi tuyến thiếu cơ cấu chấp hành.
  • Kỹ sư R&D và Cơ điện tử (Industry & Robotics Engineers): Cung cấp giải thuật điều khiển chuyển động có độ bền vững cao, khả năng kháng nhiễu tốt để ứng dụng trực tiếp cho các hệ thống tay máy công nghiệp, robot gắp thả tốc độ cao, hệ thống cân bằng giàn khoan trên biển hoặc thiết bị ổn định camera chống rung (gimbal 3 trục).
  • Giảng viên và Sinh viên/NCS ngành Tự động hóa: Là bộ học liệu chuyên sâu phục vụ giảng dạy các môn học: Điều khiển phi tuyến nâng cao, Điều khiển tối ưu, Mô hình hóa và mô phỏng hệ thống cơ điện tử.

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

1. Phát hiện quan trọng nhất của công trình nghiên cứu này là gì?

Trả lời: Đề tài đã xây dựng và chứng minh thành công bằng toán học giải thuật Điều khiển trượt tối ưu (OSMC) cho hệ Pendubot 3 khâu, khắc phục hoàn toàn nhược điểm của phương pháp trượt thứ bậc truyền thống khi bảo đảm tất cả các biến sai số trạng thái đều triệt tiêu về 0, giúp hệ thống bám quỹ đạo hình sin chỉ trong 2 giây.

2. Phương pháp xác định thông số mặt trượt trong đề tài có gì vượt trội?

Trả lời: Thay vì lựa chọn thông số mặt trượt $c_i$ theo phương pháp thử sai cảm tính, nghiên cứu đã thiết lập mối quan hệ giải tích trực tiếp giữa vector hệ số mặt trượt với ma trận độ lợi tối ưu LQR $K$, sau đó sử dụng Giải thuật Di truyền (GA) để tối ưu hóa toàn cục các ma trận trọng số $Q$ và $L$.

3. Kết quả của nghiên cứu này có thể tổng quát hóa cho các hệ thống khác không?

Trả lời: Hoàn toàn có thể. Khung thuật toán OSMC đề xuất có cấu trúc tổng quát, sẵn sàng áp dụng cho toàn bộ lớp đối tượng SIMO phi tuyến bậc cao khác như: Hệ con lắc ngược quay (Rotary Inverted Pendulum), Hệ con lắc trên xe con (Inverted Pendulum on Cart), Hệ bóng-thanh (Ball and Beam), cũng như các cánh tay robot công nghiệp nhiều khớp.

4. Bước tiếp theo của hướng nghiên cứu này là gì?

Trả lời: Nhóm nghiên cứu định hướng triển khai thực nghiệm giải thuật trên mô hình phần cứng thực tế (yêu cầu vi xử lý DSP/FPGA tốc độ cao có khả năng xử lý chu kỳ lấy mẫu $< 10\text{ ms}$ và cảm biến encoder góc phân giải cao), đồng thời mở rộng thuật toán cho hệ Pendubot 4 khâu và 5 khâu.

5. Tại sao nghiên cứu không sử dụng mạng nơ-ron (AI/Neural Networks) thay cho điều khiển trượt?

Trả lời: Mặc dù AI rất mạnh, các bộ điều khiển dựa trên mạng nơ-ron/Fuzzy thuần túy hoạt động như "hộp đen" (black-box), đòi hỏi tập dữ liệu huấn luyện khổng lồ, tiêu tốn tài nguyên phần cứng lớn và khó chứng minh ổn định tiệm cận tuyệt đối. Giải thuật trượt tối ưu dựa trên phương trình toán học giải tích bảo đảm 100% độ tin cậy vận hành thông qua tiêu chuẩn ổn định Lyapunov.

Kết luận (Conclusion)

Công trình nghiên cứu "Điều khiển trượt tối ưu bám quỹ đạo cho hệ Pendubot" (Mã số: T2019-04NCS) tại Trường Đại học Sư phạm Kỹ thuật TP.HCM đã giải quyết xuất sắc bài toán điều khiển bám quỹ đạo cho hệ SIMO phi tuyến bậc cao. Bằng việc kết hợp hài hòa giữa mô hình hóa Euler-Lagrange, lý thuyết điều khiển trượt bền vững, hàm tối ưu Lyapunov và giải thuật di truyền GA, nghiên cứu đã mang lại một bước tiến quan trọng cả về mặt lý thuyết giải tích lẫn kỹ thuật điều khiển thực nghiệm mô phỏng. Đây là tài liệu khoa học giá trị, mở ra nhiều tiềm năng ứng dụng đột phá trong việc phát triển các thế hệ robot tự hành và cơ cấu chấp hành chính xác cao trong tương lai.