Nghiên Cứu Điều Khiển Cánh Tay Robot Thiếu Dẫn Động Hai Bậc Tự Do (Pendubot): Từ Mô Hình Toán Học Đến Thực Nghiệm Nhúng


Tóm tắt nghiên cứu

Nghiên cứu tập trung giải quyết câu hỏi trọng tâm: Làm thế nào để mô hình hóa chính xác, nhận dạng hệ thống và thiết kế thuật toán điều khiển ổn định cánh tay robot thiếu dẫn động hai bậc tự do (Pendubot) tại các điểm cân bằng bất ổn định trong điều kiện có nhiễu thực tế?

Hệ Pendubot thuộc lớp đối tượng SIMO (Single-Input Multiple-Output) phi tuyến cao với một ngõ vào điều khiển tại khớp một và khớp hai quay tự do. Phương pháp luận nghiên cứu kết hợp giữa xây dựng mô hình động lực học Lagrange loại 1, khảo sát nhận dạng hệ thống đa cấu trúc (phi tham số: xung, nấc, quét tần số sóng sin; tham số: ARX, ARMAX, OE, BJ), và áp dụng thuật toán điều khiển lai hai giai đoạn: Swing-up theo nguyên lý năng lượng kết hợp điều khiển tối ưu tuyến tính toàn phương LQR (Linear-Quadratic Regulator) với ma trận trọng số tối ưu hóa bằng giải thuật di truyền GA (Genetic Algorithm).

Kết quả thực nghiệm trên nền tảng vi điều khiển Texas Instruments LAUNCHXL-F28379D (C2000) và phần mềm Matlab/Simulink thông qua Embedded Coder chứng minh: mô hình Output Error (OE) cho độ bám nhận dạng cao nhất (best fit đạt 99.4%); bộ điều khiển LQR-GA đưa hệ thống thăng bằng chính xác tại cả hai vị trí kỳ dị TOP ($q_1 = \pi/2, q_2 = 0$) và MID ($q_1 = -\pi/2, q_2 = \pi$).

Nghiên cứu cung cấp giải pháp toàn diện từ lý thuyết giải tích đến hiện thực hóa phần cứng nhúng thời gian thực, mở ra nền tảng quan trọng cho việc làm chủ công nghệ điều khiển các hệ cơ điện tử phức tạp.


Bối cảnh và tầm quan trọng

Trong lý thuyết điều khiển học hiện đại và kỹ thuật robot, hệ thống thiếu dẫn động (underactuated systems) – hệ thống có số lượng cơ cấu chấp hành ít hơn số bậc tự do – luôn là một trong những đối tượng nghiên cứu thách thức nhất. Đặc tính phi tuyến mạnh, sự ràng buộc gia tốc không tích phân được (non-holonomic constraints) cùng tính nhạy cảm cao với các tác động ngoại cảnh khiến các thuật toán điều khiển tuyến tính kinh điển không thể áp dụng trực tiếp nếu thiếu mô hình hóa chính xác và kỹ thuật tuyến tính hóa tối ưu.

       [Khớp 1: Có động cơ DC]
                  |
              (Link 1)
                  |
       [Khớp 2: Tự do - Không động cơ]
                  |
              (Link 2)

Mặc dù có nhiều tài liệu học thuật tiếp cận bài toán điều khiển con lắc ngược và cánh tay robot, khoảng trống nghiên cứu (research gap) vẫn hiện hữu rõ rệt giữa mô phỏng thuần túy trên máy tính và việc triển khai trên phần cứng thực tế. Sai số mô hình do ma sát khớp, điện cảm cuộn dây động cơ DC, độ phân giải bộ đọc xung (Encoder), cùng tác động của các dạng nhiễu ngẫu nhiên (nhiễu DC dạng Mean và nhiễu AC dạng Variance) thường khiến các bộ điều khiển lý thuyết bị mất ổn định trên mô hình thực nghiệm. Hơn nữa, việc so sánh có hệ thống giữa các phương pháp nhận dạng hệ thống tham số và phi tham số trên mô hình Pendubot vẫn chưa được hệ thống hóa đầy đủ.

Trong kỷ nguyên Công nghiệp 4.0, sự phát triển của công nghệ vi xử lý DSP tốc độ cao và quy trình thiết kế dựa trên mô hình (Model-Based Design) đòi hỏi các nghiên cứu phải tiệm cận quy trình kỹ thuật chuẩn công nghiệp. Đề tài này giải quyết trực tiếp thách thức trên, mang lại giá trị học thuật và thực tiễn sâu sắc: vừa xây dựng mô hình giáo cụ trực quan phục vụ đào tạo kỹ sư chất lượng cao, vừa tạo bước đệm hoàn thiện các giải thuật điều khiển phi tuyến ứng dụng cho các hệ thăng bằng động như tên lửa đẩy, robot hai chân, và hệ thống cẩu trục công nghiệp.


Methodology và approach

Nghiên cứu triển khai quy trình kỹ thuật khép kín theo 4 giai đoạn logic chặt chẽ:

[Mô hình hóa toán học (Lagrange, DH)]
[Nhận dạng hệ thống (Tham số & Phi tham số)]
[Thiết kế bộ điều khiển (Energy Swing-up & LQR-GA)]
[Nhúng mã nguồn C (LAUNCHXL-F28379D & Matlab Embedded Coder)]

1. Mô hình hóa toán học và phân tích động lực học

Nghiên cứu ứng dụng quy tắc Denavit-Hartenberg (bảng DH) để thiết lập ma trận biến đổi tọa độ thuần nhất, kết hợp phương trình vi phân Lagrange loại 1 dựa trên hàm Lagrange $L = K - V$ (với $K$ là tổng động năng và $V$ là thế năng hệ thống). Mô hình phi tuyến 2 bậc tự do được ghép nối với phương trình cân bằng điện - cơ của động cơ điện một chiều 24VDC: $$\tau_1 = -J\ddot{q}_1 - \theta_6\dot{q}_1 + \theta_7 u$$ Hệ thống được chứng minh tính điều khiển được (Controllability) thông qua tiêu chuẩn Kalman $Rank(R) = 4 = n$ tại cả hai điểm làm việc bất ổn định TOP và MID.

2. Kỹ thuật nhận dạng hệ thống đa phương pháp

Nghiên cứu tiến hành thu thập và phân tích dữ liệu thực nghiệm qua hai nhóm kỹ thuật:

  • Nhận dạng phi tham số (Non-parametric): Phân tích đáp ứng xung (Impulse response), đáp ứng nấc (Step response), phân tích tương quan (Correlation analysis), và phân tích đáp ứng tần số bằng phương pháp kiểm tra sóng sin (Sine wave test) với 40 tập dữ liệu thực nghiệm. Tất cả được kiểm chứng qua 4 kịch bản nhiễu: không nhiễu, nhiễu DC (Mean), nhiễu AC (Variance), và kết hợp DC + AC.
  • Nhận dạng có tham số (Parametric): Khảo sát 4 cấu trúc ước lượng hàm truyền rời rạc bậc 4: Auto-Regressive with External Input (ARX), Auto-Regressive Moving Average with External Input (ARMAX), Output Error (OE), và Box-Jenkins (BJ).

3. Thiết kế chiến lược điều khiển tối ưu phân tầng

  • Giai đoạn Swing-up: Sử dụng nguyên lý năng lượng điều khiển phi tuyến dựa trên hàm Lyapunov nhằm tích lũy động năng, đưa hai thanh con lắc từ vị trí buông tự do hướng xuống lên lân cận góc thẳng đứng.
  • Giai đoạn Cân bằng (Balancing): Thiết kế bộ điều khiển tối ưu toàn phương tuyến tính LQR ($u = -Kx$) nhằm tối thiểu hóa hàm mục tiêu chi phí $J = \int (x^T Q x + u^T R u) dt$.
  • Tối ưu hóa tham số bằng GA: Tích hợp giải thuật di truyền (Genetic Algorithm) để tìm kiếm tự động cặp ma trận trọng số $Q$ và $R$, loại bỏ phương pháp thử sai truyền thống và khắc phục sai lệch cơ khí thực tế.

4. Nền tảng thực thi phần cứng và kiểm chứng

Hệ thống sử dụng vi điều khiển 32-bit kép DSP TMS320F28379D xung nhịp 200MHz, tích hợp bộ đọc xung tốc độ cao eQEP đọc tín hiệu từ Encoder Omron E6B2-CWZ6C (1000 xung/vòng), mạch công suất cầu H lái MOSFET IRF3205, và giao tiếp dữ liệu UART qua chip CP2102. Mã điều khiển được tự động sinh và nạp trực tiếp từ Simulink qua gói thư viện Embedded Coder Support Package for Texas Instruments C2000 Processors.


Phát hiện chính

Quá trình mô phỏng và thực nghiệm trên đối tượng thực đem lại 4 phát hiện khoa học quan trọng:

+-------------------------------------------------------------+
| Cấu trúc Mô hình       | Độ trùng khớp (Best Fit %)          |
+------------------------+------------------------------------+
| Output Error [oe421]   | 99.4%                              |
| ARX [arx441]           | 96.3%                              |
| ARMAX [amx2221]        | 95.3%                              |
| Box-Jenkins [bj42122]  | 89.4%                              |
+-------------------------------------------------------------+

1. Cấu trúc Output Error (OE) vượt trội tuyệt đối trong nhận dạng tham số

Trong số các mô hình tham số khảo sát, cấu trúc OE bậc 4 [oe421] đạt chỉ số khớp (Best Fit) cao nhất lên đến 99.4% trong điều kiện không nhiễu, vượt trội hơn so với ARX (96.3%), ARMAX (95.3%) và Box-Jenkins (89.4%). Khi xuất hiện nhiễu DC và AC, cấu trúc OE tiếp tục thể hiện sự ổn định với sai số ngõ ra thấp nhất, phản ánh chính xác hàm truyền đạt thực tế của cánh tay robot bậc 4.

2. Ưu thế của phương pháp kiểm tra sóng sin trong phân tích tần số

Đối với các phương pháp phi tham số, phân tích tương quan và đáp ứng xung bộc lộ độ nhạy cảm cao với nhiễu Variance (AC), dẫn đến sai số nhận dạng lớn. Ngược lại, phương pháp quét tần số sóng sin với thời gian lấy mẫu thích hợp (500s đến 30.000s) đã tái tạo giản đồ Bode biên độ và Bode pha tại Link 1 và Link 2 với độ tương quan cao so với hệ thống thực, chứng minh đây là công cụ nhận dạng phi tham số tin cậy nhất cho các cơ cấu khớp nối phi tuyến.

3. Thuật toán di truyền (GA) giải quyết triệt để bài toán chọn ma trận LQR

Việc tối ưu hóa ma trận $Q$ và $R$ thông qua GA giúp tìm ra bộ hệ số hồi tiếp trạng thái $K$ tối ưu. Khi áp dụng lên thực nghiệm, thời gian xác lập trạng thái cân bằng được rút ngắn đáng kể, độ vọt lố góc quay của Link 1 và Link 2 giảm thiểu rõ rệt, đồng thời giữ biên độ điện áp kích từ động cơ luôn nằm trong giới hạn an toàn ($< 24V$).

4. Cơ chế chuyển trạng thái mượt mà giữa Swing-up và Balancing

Sự kết hợp giữa giải thuật năng lượng và LQR tạo nên quá trình chuyển đổi mượt mà khi hệ thống di chuyển từ vùng góc lệch lớn vào lân cận điểm cân bằng ($|q_1 - q_{1_target}| \le \epsilon$). Hiện tượng giật cục cơ khí (chattering) và quá dòng trên cầu H được triệt tiêu hoàn toàn trên cả hai vị trí TOP và MID.


Đóng góp khoa học

Nghiên cứu mang lại những đóng góp rõ nét cả về lý thuyết, phương pháp luận lẫn ứng dụng thực tiễn:

                     ĐÓNG GÓP KHOA HỌC
[Lý thuyết & Mô hình]   [Phương pháp luận]     [Ứng dụng thực tế]
• Chuẩn hóa Lagrange-DH • So sánh nhận dạng     • Bộ phần cứng nhúng
• Phân tích kỳ dị TOP/MID • LQR kết hợp GA      • Model-Based Design
  • Về mặt lý thuyết: Chuẩn hóa toàn bộ hệ phương trình vi phân mô tả chuyển động của cánh tay robot 2 bậc tự do có tính đến đặc tính phi tuyến của động cơ DC. Minh chứng tường minh tính điều khiển được của hệ SIMO tại các điểm kỳ dị không ổn định.
  • Về mặt phương pháp luận: Xây dựng một quy trình đối sánh toàn diện giữa các kỹ thuật nhận dạng hệ thống hiện đại, cung cấp bộ dữ liệu thực nghiệm về đáp ứng của hệ thiếu dẫn động dưới nhiều điều kiện nhiễu khác nhau.
  • Về mặt ứng dụng thực tiễn: Hiện thực hóa thành công sản phẩm phần cứng Pendubot hoàn chỉnh kết nối với vi điều khiển công nghiệp TMS320F28379D. Thiết lập chuẩn quy trình thiết kế Model-Based Design, cho phép lập trình, mô phỏng và nạp mã điều khiển tự động hoàn toàn trên nền tảng Matlab/Simulink.
  • Ý nghĩa đào tạo và xã hội: Mô hình đóng vai trò là giáo cụ trực quan hiện đại cho các phòng thí nghiệm Tự động hóa và Robot học, hỗ trợ sinh viên và học viên cao học kiểm chứng các thuật toán điều khiển tiên tiến (Backstepping, Sliding Mode, Neural Network, Fuzzy Logic).

Đối tượng quan tâm

Nghiên cứu này đặc biệt hữu ích và mang lại giá trị trực tiếp cho các nhóm đối tượng sau:

  • Giảng viên và Nhà nghiên cứu học thuật: Sử dụng làm tài liệu tham khảo chuyên sâu về mô hình hóa hệ thống SIMO, nhận dạng tham số và điều khiển phi tuyến. Cung cấp dữ liệu gốc để đối chiếu khi phát triển các thuật toán điều khiển thông minh mới.
  • Kỹ sư R&D trong ngành Cơ điện tử và Robotics: Tham khảo phương pháp thiết kế phần cứng, mạch công suất bảo vệ opto cách ly, và kỹ thuật đọc xung encoder đa kênh thời gian thực trên dòng chip điều khiển động cơ C2000.
  • Sinh viên, Học viên cao học chuyên ngành Tự động hóa & Robot: Nắm bắt quy trình thực hành từ lý thuyết giải tích đến triển khai nhúng thực tế thông qua công cụ Matlab Embedded Coder, rút ngắn thời gian làm quen với các dòng vi điều khiển DSP chuyên dụng.
  • Lợi ích mang lại: Cung cấp giải pháp tối ưu hóa chi phí chế tạo mô hình thí nghiệm điều khiển tự động mà vẫn đảm bảo độ chính xác và tính năng kỹ thuật tương đương các hệ thống thương mại đắt tiền.

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

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

Phát hiện quan trọng nhất là việc xác định mô hình Output Error (OE) có cấu trúc bậc [oe421] là mô hình tham số tối ưu nhất cho hệ Pendubot (đạt Best Fit 99.4%), kết hợp với giải thuật di truyền (GA) để tự động hóa hoàn toàn việc tìm ma trận trọng số LQR, giúp hệ thống cân bằng ổn định ở cả hai vị trí khó (TOP và MID) trên mô hình thực nghiệm.

2. Phương pháp nghiên cứu có điểm gì đặc biệt so với các đề tài trước?

Nghiên cứu không dừng lại ở mô phỏng lý thuyết mà thực hiện quy trình khép kín: khảo sát nhận dạng hệ thống dưới nhiều kịch bản nhiễu thực tế (DC và AC), đồng thời áp dụng quy trình Model-Based Design nhúng mã C trực tiếp từ Simulink xuống vi điều khiển Texas Instruments F28379D, giúp đảm bảo tính đồng nhất giữa mô phỏng và thực nghiệm.

3. Kết quả điều khiển của đề tài có thể tổng quát hóa (generalize) cho các hệ thống khác không?

Hoàn toàn có thể. Mô hình toán học, phương pháp nhận dạng và thuật toán điều khiển Swing-up kết hợp LQR-GA trong đề tài này có thể tổng quát hóa và áp dụng trực tiếp cho các hệ thống thiếu dẫn động khác như: con lắc ngược trên xe (Inverted Pendulum on Cart), con lắc quay (Furuta Pendulum), hệ bóng và thanh (Ball and Beam), hoặc cánh tay robot công nghiệp nhiều bậc tự do.

4. Hướng phát triển tiếp theo (Next Steps) của nghiên cứu là gì?

Các bước tiếp theo bao gồm: thử nghiệm các thuật toán điều khiển thích nghi (Adaptive Control), điều khiển trượt thích nghi (Adaptive Sliding Mode Control), và ứng dụng mạng nơ-ron sâu (Deep Reinforcement Learning) để điều khiển bám quỹ đạo phức tạp hoặc tự phục hồi khi có ngoại lực va chạm mạnh.

5. Ứng dụng thực tế ngoài phòng thí nghiệm của nghiên cứu này là gì?

Nguyên lý điều khiển hệ thiếu dẫn động trong nghiên cứu được ứng dụng trực tiếp trong việc thiết kế hệ thống thăng bằng cho robot hình người (Humanoid Robot), hệ thống ổn định tên lửa đẩy trong giai đoạn phóng, hệ thống giảm rung lắc cho cẩu trục cảng biển, và hệ thống ổn định góc nghiêng trên tàu thủy chuyên dụng.


Kết luận

Nghiên cứu đã giải quyết trọn vẹn bài toán điều khiển hệ cánh tay robot thiếu dẫn động hai bậc tự do (Pendubot) từ khâu giải tích toán học, nhận dạng hệ thống đa phương pháp đến hiện thực hóa phần cứng nhúng thời gian thực. Việc tích hợp thành công giải thuật Swing-up năng lượng và bộ điều khiển tối ưu LQR-GA trên vi xử lý LAUNCHXL-F28379D khẳng định tính đúng đắn và hiệu quả vượt trội của phương pháp đề xuất.

Trong tương lai, việc mở rộng nghiên cứu sang các thuật toán điều khiển phi tuyến bền vững và trí tuệ nhân tạo sẽ tiếp tục hoàn thiện nền tảng công nghệ này. Các nhóm nghiên cứu và đơn vị đào tạo được khuyến khích ứng dụng mô hình thực nghiệm này nhằm nâng cao chất lượng nghiên cứu khoa học và phát triển các ứng dụng cơ điện tử thông minh.