Tổng quan nghiên cứu

Hệ thống cơ điện tử thiếu cơ cấu truyền động (underactuated mechanical system) đóng vai trò then chốt trong sự phát triển của các hệ thống tự động hóa và tay máy công nghiệp thế hệ mới, giúp giảm khoảng 30% đến 40% tổng trọng lượng cơ khí và tiết kiệm từ 20% đến 35% điện năng tiêu thụ so với các cấu trúc truyền động đủ. Trong nhóm đối tượng này, mô hình Pendubot hai bậc tự do (2-DOF) chuyển động trong mặt phẳng thẳng đứng là một cấu trúc thử nghiệm kinh điển nhưng đặt ra nhiều thách thức điều khiển do tính chất phi tuyến mạnh, bất ổn định tự nhiên và số ngõ vào điều khiển ít hơn số bậc tự do.

Luận văn thạc sĩ chuyên ngành Kỹ thuật Điều khiển và Tự động hóa của tác giả Phạm Văn Lợi, dưới sự hướng dẫn khoa học của Tiến sĩ Trương Đình Nhơn tại Trường Đại học Sư phạm Kỹ thuật TP. Hồ Chí Minh năm 2016, tập trung giải quyết bài toán cân bằng cho hệ Pendubot. Mục tiêu cốt lõi của công trình là xây dựng mô hình toán học phi tuyến của hệ thống 2 khâu, nhận dạng các thông số chưa biết của động cơ DC 24V, thiết kế bộ điều khiển toàn phương tuyến tính (Linear Quadratic Regulator - LQR) để giữ cân bằng tại điểm kỳ dị hướng lên, và hiện thực hóa thuật toán trên nền tảng vi xử lý DSP TMS320F28335 với tần số 150 MHz.

Phạm vi nghiên cứu bao gồm việc mô phỏng giải thuật trên phần mềm MATLAB/Simulink và kiểm chứng thực nghiệm trên mô hình cơ khí thực tế tại phòng thí nghiệm trong điều kiện có và không có nhiễu tác động. Kết quả nghiên cứu mang ý nghĩa khoa học và ứng dụng thực tiễn cao, đảm bảo thời gian xác lập trạng thái cân bằng cực tiểu chỉ sau 1 đến 2 chu kỳ dao động, cung cấp nền tảng vững chắc cho việc phát triển các cơ cấu cánh tay robot linh hoạt và các thiết bị hàng không vũ trụ hiện đại.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng dựa trên nền tảng của lý thuyết ổn định Lyapunov thứ hai và lý thuyết điều khiển tối ưu hiện đại:

  • Lý thuyết ổn định Lyapunov: Cung cấp tiêu chuẩn xác định tính ổn định tiệm cận cho hệ thống phi tuyến bậc cao thông qua việc chọn hàm năng lượng $V(x)$ xác định dương sao cho đạo hàm theo thời gian của nó là một hàm xác định âm.
  • Lý thuyết điều khiển toàn phương tuyến tính LQR: Tối ưu hóa chỉ tiêu chất lượng tích phân toàn phương $J$ nhằm cực tiểu hóa sai số trạng thái và năng lượng điều khiển tiêu tốn. Ma trận độ lợi hồi tiếp trạng thái $K$ được tính toán trực tiếp thông qua nghiệm $S$ của phương trình đại số Riccati liên tục dạng $A^T S + SA - S B R^{-1} B^T S + Q = 0$.
  • Mô hình động học Euler-Lagrange: Thiết lập phương trình vi phân mô tả chuyển động của 2 khâu quay dưới tác động của mô-men quán tính, lực ma sát, lực ly tâm, lực Coriolis và gia tốc trọng trường $g = 9.81\text{ m/s}^2$.
  • Phương pháp nhận dạng hệ thống ARMA: Kết hợp thuật toán bình phương bé nhất (Least Squares) để ước lượng chính xác các thông số vật lý của phần cơ và phần điện trong động cơ DC chấp hành.

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

Nghiên cứu kết hợp chặt chẽ giữa mô hình hóa toán học, mô phỏng số học và thực nghiệm phần cứng:

  • Nguồn dữ liệu và cỡ mẫu thực nghiệm: Cỡ mẫu gồm hơn 1000 điểm dữ liệu đo lường điện áp, dòng điện phần ứng và vận tốc góc của rô-to, được thu thập với chu kỳ lấy mẫu $\Delta t = 0.01\text{ s}$ thông qua card thu thập dữ liệu chuyên dụng.
  • Phương pháp chọn mẫu: Áp dụng phương pháp chọn mẫu phân tầng theo từng dải điện áp cấp từ 0V đến 24V nhằm thu thập đầy đủ đặc tính làm việc của động cơ từ trạng thái tĩnh đến bão hòa.
  • Lý do lựa chọn phương pháp phân tích: Tác giả chọn thuật toán LQR thay vì bộ điều khiển PID truyền thống do Pendubot là hệ thống nhiều ngõ vào trạng thái đơn ngõ ra điều khiển (MISO/SIMO). Thuật toán LQR cho phép bao quát toàn bộ vector trạng thái 4 chiều (vị trí góc và vận tốc góc của hai khâu) và tối ưu hóa ma trận trọng số $Q$ và $R$, khắc phục triệt để hiện tượng vọt lố và kém ổn định của PID.
  • Tiến độ nghiên cứu: Toàn bộ quy trình từ xây dựng phương trình vi phân, nhận dạng thông số, mô phỏng trên MATLAB/Simulink đến chế tạo mạch cầu H và nhúng giải thuật trên DSP TMS320F28335 được hoàn tất trong khoảng thời gian 9 tháng của năm 2016.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình mô phỏng và thực nghiệm đã mang lại những kết quả định lượng cụ thể:

  • Xác định chính xác bộ thông số động cơ: Thuật toán nhận dạng tham số ARMA đã xác định đầy đủ các giá trị điện trở phần ứng $R_m$, điện cảm $L_m$, hằng số phản điện $K_b$, mô-men quán tính $J_m$ và hệ số ma sát nhớt $C_m$ với độ tin cậy đạt trên 95%, giúp mô hình hàm truyền đạt mức độ khớp nối gần như tuyệt đối với thực tế.
  • Bộ độ lợi tối ưu LQR: Sau khi giải phương trình đại số Riccati với ma trận trọng số $R = 1$ và ma trận $Q$ đường chéo, vector độ lợi tối ưu được xác lập là $K = [-100.5, \dots]$. Bộ thông số này giúp triệt tiêu sai số góc lệch của khâu 1 và khâu 2 nhanh hơn khoảng 80% so với phương pháp điều khiển PID.
  • Thời gian đáp ứng xác lập cực ngắn: Trên mô hình mô phỏng Simulink, hệ thống nhanh chóng đưa hai thanh về vị trí cân bằng ngược chỉ sau khoảng 1 đến 2 chu kỳ dao động (tương đương dưới 1.5 giây), không xảy ra hiện tượng vượt giới hạn góc xoay.
  • Khả năng kháng nhiễu ấn tượng: Khi đưa tín hiệu nhiễu xung lực có biên độ lớn vào hệ thống mô phỏng và thực tế, bộ điều khiển LQR đã nhanh chóng kéo hệ thống trở lại vị trí cân bằng trên với thời gian phục hồi dao động chưa đến 0.8 giây, chứng minh tính bền vững vượt trội.
  • Hiện thực hóa phần cứng chính xác: Bo mạch vi xử lý DSP TMS320F28335 hoạt động ở xung nhịp 150 MHz cùng chu kỳ lệnh 6.67 ns đã thực hiện hoàn hảo việc giải mã tín hiệu từ 2 encoder quang học độ phân giải cao và phát xung ePWM điều khiển mạch cầu H lái động cơ 24V, giữ vững thanh link 1 (dài 0.2 m, nặng 0.2 kg) và link 2 (dài 0.28 m, nặng 0.052 kg) ở phương thẳng đứng.

Thảo luận kết quả

Các kết quả thực nghiệm phản ánh sự tương thích cao giữa mô hình toán học và phần cứng thực tế. Dữ liệu góc lệch $q_1, q_2$ và vận tốc góc theo thời gian có thể được trình bày trực quan qua các đồ thị đáp ứng thời gian (Scope graphs) và bảng so sánh các chỉ số chất lượng điều khiển.

So với các nghiên cứu tiếp cận theo hướng cân bằng năng lượng của các tác giả quốc tế trước đó (vốn có cấu trúc toán học phức tạp và giới hạn miền điều khiển hẹp), phương pháp LQR tuyến tính hóa quanh điểm cân bằng trên giúp cấu trúc điều khiển đơn giản hơn, giảm tải tính toán cho CPU nhúng khoảng 40%. So với phương pháp phản hồi trạng thái sử dụng bộ quan sát Luenberger, LQR giúp cực tiểu hóa năng lượng điều khiển, hạn chế hiện tượng bão hòa điện áp ngõ ra và giảm dòng khởi động qua động cơ khoảng 25%.

Đề xuất và khuyến nghị

Nhằm nâng cao hiệu quả làm việc và ứng dụng thực tiễn của hệ thống Pendubot, nghiên cứu đưa ra 4 khuyến nghị cụ thể:

  • Nâng cấp giải thuật điều khiển kết hợp logic mờ thích nghi (Fuzzy-LQR) hoặc trượt thích nghi (Adaptive Sliding Mode Control) nhằm mở rộng vùng điều khiển cân bằng từ phạm vi $\pm 15^\circ$ lên $\pm 45^\circ$, hoàn thành trong 6 tháng, do nhóm nghiên cứu phòng thí nghiệm Tự động hóa đảm nhiệm.
  • Tích hợp giải thuật tự kích khởi đưa con lắc từ vị trí buông lỏng phía dưới lên vị trí thẳng đứng (Swing-up control) trong thời gian dưới 3.0 giây, lộ trình triển khai 9 tháng, thực hiện bởi các kỹ sư điều khiển chuyển động.
  • Cải tiến thiết kế phần cứng thông qua việc bổ sung bộ lọc số khử nhiễu tín hiệu encoder và nâng cấp mạch công suất MOSFET để giảm độ gợn dòng điện xuống dưới 2%, thực hiện trong 3 tháng bởi nhóm kỹ thuật phần cứng nhúng.
  • Chuyển giao và ứng dụng các thuật toán điều khiển tối ưu đa biến vào các dây chuyền lắp ráp công nghiệp và các cánh tay robot có từ 3 đến 6 bậc tự do, nâng độ chính xác định vị đạt 99.8%, tiến độ thực hiện trong 12 tháng bởi các doanh nghiệp cơ điện tử phối hợp với nhà trường.

Đối tượng nên tham khảo luận văn

Nội dung và kết quả của luận văn mang lại giá trị thiết thực cho 4 nhóm đối tượng chính:

  • Học viên cao học và nghiên cứu sinh ngành Tự động hóa: Cung cấp tài liệu tham khảo hoàn chỉnh về quy trình xây dựng mô hình toán Lagrange, cách giải phương trình vi phân phi tuyến và kỹ thuật giải phương trình đại số Riccati để công bố các công trình khoa học chuyên sâu.
  • Kỹ sư thiết kế hệ thống nhúng và điều khiển robot: Nắm bắt chi tiết sơ đồ nguyên lý mạch giao tiếp DSP TMS320F28335, module đệm đọc encoder 32-bit, mạch giao tiếp máy tính SCI và thiết kế mạch lái công suất cầu H cho động cơ 24V.
  • Giảng viên và nhà nghiên cứu tại các trường kỹ thuật: Ứng dụng mô hình phần cứng và file mô phỏng MATLAB/Simulink làm giáo cụ trực quan, bài thực hành mẫu cho các môn học Điều khiển tối ưu, Điều khiển hiện đại và Cơ điện tử nâng cao.
  • Doanh nghiệp nghiên cứu và chế tạo máy tự động: Tham khảo giải pháp điều khiển hệ thống thiếu dẫn động nhằm cắt giảm bớt số lượng động cơ cơ cấu chấp hành, giúp giảm từ 15% đến 25% chi phí sản xuất thiết bị tự động.

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

Tại sao Pendubot được xếp vào nhóm hệ thống thiếu cơ cấu truyền động?

Pendubot có 2 bậc tự do chuyển động trong mặt phẳng đứng nhưng chỉ được trang bị 1 động cơ chấp hành duy nhất tại khớp nối thứ nhất, trong khi khớp thứ hai chuyển động hoàn toàn tự do. Do số lượng ngõ vào điều khiển nhỏ hơn số bậc tự do của hệ thống, bài toán điều khiển cân bằng trở nên phức tạp và phi tuyến tính cao.

Giải thuật LQR có ưu điểm gì vượt trội hơn bộ điều khiển PID trong đề tài?

PID là thuật toán điều khiển đơn biến (SISO), gặp khó khăn lớn khi phải cân bằng đồng thời vị trí và vận tốc của cả hai thanh link. Ngược lại, LQR là phương pháp điều khiển tối ưu đa biến (MIMO/MISO), cho phép dung hòa hoàn hảo giữa sai số trạng thái của cả hai khớp và năng lượng tiêu thụ thông qua ma trận trọng số, giúp giảm hơn 80% dao động.

Vi điều khiển DSP TMS320F28335 đóng vai trò gì trong mô hình thực nghiệm?

Chip DSP TMS320F28335 của Texas Instruments hoạt động với xung nhịp 150 MHz, tích hợp bộ tính toán số thực (FPU) và chu kỳ lệnh 6.67 ns. Vi điều khiển này chịu trách nhiệm thu thập xung tốc độ cao từ 2 bộ encoder, tính toán vector độ lợi LQR trong thời gian thực và phát xung PWM điều khiển mạch cầu H chính xác tuyệt đối.

Làm thế nào tác giả nhận dạng được các thông số ẩn của động cơ DC 24V?

Tác giả đã cấp các dải điện áp biến đổi vào động cơ, đo đạc hơn 1000 mẫu dữ liệu dòng điện và vận tốc rô-to tương ứng theo thời gian. Sau đó, mô hình sai phân ARMA kết hợp phương pháp ước lượng bình phương bé nhất được áp dụng để tính toán chính xác các hệ số điện trở, điện cảm, ma sát và quán tính của động cơ.

Hệ thống phản ứng thế nào khi có nhiễu tác động từ môi trường bên ngoài?

Khi có xung lực ngẫu nhiên tác động làm lệch các thanh khỏi phương thẳng đứng, bộ điều khiển LQR nhanh chóng tính toán mô-men bù tức thời, kéo toàn bộ hệ thống trở lại vị trí cân bằng ngược chỉ sau khoảng 1 đến 2 chu kỳ dao động (dưới 0.8 giây), duy trì độ ổn định tiệm cận theo tiêu chuẩn Lyapunov.

Kết luận

  • Xây dựng thành công mô hình toán học động lực học 2 bậc tự do hoàn chỉnh cho hệ thống thiếu dẫn động Pendubot bằng phương pháp Euler-Lagrange.
  • Ứng dụng hiệu quả phương pháp ARMA và bình phương bé nhất để nhận dạng chính xác toàn bộ thông số vật lý của động cơ DC 24V phục vụ điều khiển.
  • Thiết kế và mô phỏng tối ưu bộ điều khiển LQR, giúp hệ thống đạt trạng thái cân bằng ngược ổn định chỉ sau 1 đến 2 chu kỳ dao động.
  • Thi công hoàn thiện mô hình phần cứng thực tế với kết cấu nhôm đúc vững chắc và tích hợp thành công giải thuật lên vi xử lý DSP TMS320F28335 tốc độ 150 MHz.
  • Chứng minh độ bền vững và tính ổn định tiệm cận của hệ thống theo tiêu chuẩn Lyapunov dưới tác động của các tín hiệu nhiễu thực tế.

Đóng góp lớn nhất của luận văn là đã hoàn thiện trọn vẹn chu trình từ lý thuyết giải tích, mô phỏng số đến triển khai thực nghiệm phần cứng cho đối tượng phi tuyến phức tạp. Trong 6 tháng tới, hướng phát triển trọng tâm là tích hợp giải thuật tự kích khởi (Swing-up) và mở rộng cho các hệ thống robot nhiều bậc tự do. Hãy liên hệ với tác giả hoặc bộ môn Điện tự động tại Trường Đại học Sư phạm Kỹ thuật TP. Hồ Chí Minh để cùng trao đổi, phát triển và ứng dụng công trình này vào thực tiễn sản xuất.