Tổng quan nghiên cứu

Trong xu hướng tự động hóa công nghiệp hiện đại, robot và các cơ cấu chấp hành ngày càng giữ vai trò cốt lõi trong các dây chuyền lắp ráp vi mạch chính xác cũng như môi trường độc hại. Tuy nhiên, các hệ thống cơ điện tử truyền thống thường đòi hỏi số lượng động cơ tương ứng với số bậc tự do, dẫn đến việc gia tăng đáng kể trọng lượng và chi phí năng lượng. Hệ thống thiếu cơ cấu truyền động (underactuated mechanical systems) nổi lên như một hướng tiếp cận đột phá với tỷ lệ cơ cấu chấp hành trên số bậc tự do chỉ đạt 1:2, giúp cắt giảm khoảng 30% khối lượng cơ khí và tiết kiệm hơn 25% năng lượng vận hành. Điển hình cho lớp đối tượng này là mô hình Pendubot – cánh tay robot hai bậc tự do hoạt động trong mặt phẳng thẳng đứng nhưng chỉ được trang bị duy nhất một động cơ tại khớp nối thứ nhất, trong khi khớp nối thứ hai hoàn toàn tự do.

Vấn đề nghiên cứu trọng tâm là giải quyết tính phi tuyến cao, khớp nối động lực học phức tạp và trạng thái cân bằng không ổn định của hệ Pendubot. Mục tiêu cụ thể của luận văn là thiết kế thành công bộ điều khiển lai (hybrid control) gồm hai giai đoạn phối hợp: giai đoạn Swing-up đưa hai thanh liên kết từ vị trí cân bằng dưới lên lân cận vị trí thẳng đứng ngược, và giai đoạn cân bằng LQR (Linear Quadratic Regulator) để duy trì ổn định hệ thống tại điểm làm việc mong muốn. Nghiên cứu được triển khai thực nghiệm toàn diện tại Trường Đại học Sư phạm Kỹ thuật Thành phố Hồ Chí Minh trong giai đoạn 2016 - 2017, tích hợp thuật toán trên phần cứng vi xử lý DSP TMS320F28335 với tần số vận hành 150 MHz. Ý nghĩa khoa học và thực tiễn của công trình được chứng minh qua việc tối ưu hóa thời gian đưa hệ thống vào vùng ổn định xuống dưới 1,5 giây, đồng thời mở rộng miền làm việc và duy trì sai lệch góc quay nhỏ hơn 0,05 rad ngay cả khi chịu tác động của xung nhiễu ngoại cảnh.

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 trên nền tảng vững chắc của lý thuyết điều khiển phi tuyến và tối ưu hóa hệ thống động lực học. Đầu tiên, phương pháp tuyến tính hóa hồi tiếp riêng phần (Partial Feedback Linearization) được áp dụng để triệt tiêu các thành phần phi tuyến trong phương trình trạng thái của khớp có cơ cấu truyền động, tạo ra mối liên hệ tuyến tính cục bộ phục vụ quá trình Swing-up. Thứ hai, lý thuyết điều khiển tối ưu toàn phương tuyến tính LQR kết hợp với phương pháp giải phương trình đại số Riccati liên tục (Continuous Algebraic Riccati Equation - CARE) được sử dụng để xác định vector độ lợi hồi tiếp trạng thái K, nhằm cực tiểu hóa hàm chỉ tiêu chất lượng toàn phương J.

Khung nghiên cứu tích hợp năm khái niệm cốt lõi:

  • Hệ thống thiếu cơ cấu truyền động (Underactuated System): Hệ có số ngõ vào điều khiển nhỏ hơn số bậc tự do không gian.
  • Động lực học Euler-Lagrange: Mô hình hóa tương tác giữa thế năng, động năng và moment quán tính của hai thanh liên kết.
  • Tiêu chuẩn ổn định tiệm cận Lyapunov thứ hai: Đảm bảo đạo hàm hàm năng lượng xác định âm trong lân cận điểm cân bằng.
  • Ma trận trọng số trạng thái Q và ma trận trọng số điều khiển R: Thiết lập cân bằng tối ưu giữa tốc độ đáp ứng trạng thái và mức tiêu hao năng lượng điện áp.
  • Mô hình nhận dạng chuỗi thời gian ARMA (AutoRegressive Moving Average): Xác định chính xác các tham số vật lý chưa biết của động cơ chấp hành.

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

Nguồn dữ liệu của luận văn được thu thập từ mô hình phần cứng thực nghiệm thông qua module giao tiếp nối tiếp SCI kết nối với máy tính, kết hợp cảm biến quang encoder độ phân giải cao và mạch đệm đo dòng điện - điện áp chuyên dụng. Cỡ mẫu nghiên cứu bao gồm 1.000 điểm dữ liệu trạng thái được lấy mẫu liên tục trong các chu kỳ thử nghiệm thực tế. Phương pháp chọn mẫu là lấy mẫu tất định theo chu kỳ thời gian thực (deterministic periodic sampling) với chu kỳ lấy mẫu danh định 6,67 ns trên chip DSP TMS320F28335, đảm bảo bắt trọn các biến thiên cơ học nhanh và đáp ứng điện từ của động cơ một chiều.

Lý do lựa chọn phương pháp phân tích ARMA kết hợp thuật toán bình phương tối thiểu (Least Squares) là nhằm loại bỏ sai số mô hình hóa do nhà sản xuất không cung cấp thông số động cơ. Quá trình phân tích thực nghiệm đã xác định chính xác điện trở phần ứng Rm bằng 2,6824 Ohm, hằng số phản điện Kb bằng 0,0924 V.s/rad, hệ số ma sát nhớt Cm bằng 1,0171 N.m.s, moment quán tính rotor Jm bằng 3,4812e-04 kg.m2 và hệ số điện kháng Lm bằng 20,1113 H. Timeline nghiên cứu được thực hiện liên tục trong 12 tháng, bắt đầu từ việc thiết lập phương trình vi phân mô phỏng trên Matlab/Simulink, nhận dạng thông số, cho đến chế tạo cơ khí chính xác và lập trình nhúng thời gian thực.

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à kiểm nghiệm thực nghiệm đã mang lại bốn phát hiện khoa học mang tính đột phá:

  • Quá trình Swing-up diễn ra vượt trội: Thuật toán tuyến tính hóa hồi tiếp riêng phần đã điều khiển hệ Pendubot chuyển từ vị trí cân bằng dưới lên lân cận điểm cân bằng trên chỉ sau 1 đến 2 chu kỳ dao động, tương ứng thời gian xác lập cực ngắn khoảng 1,5 giây.
  • Khả năng ổn định vững chắc của bộ điều khiển LQR: Với ma trận trọng số Q dạng đường chéo có các phần tử tương ứng là 10, 1, 10, 1 và trọng số năng lượng R bằng 10, hệ thống đạt được vector độ lợi tối ưu K với hệ số phản hồi trạng thái đầu tiên lên tới 110,1744. Nhờ đó, sai số góc lệch của cả hai thanh liên kết được duy trì ổn định dưới mức 0,05 rad trong điều kiện làm việc bình thường với nguồn điện 24V DC.
  • Độ bền vững cao trước nhiễu động: Trong môi trường mô phỏng có chèn tín hiệu xung nhiễu biên độ lớn 0,5 rad tác động trực tiếp vào góc quay, hệ thống tự động cân bằng trở lại vị trí thẳng đứng trong thời gian dưới 0,8 giây, giảm thiểu hơn 40% độ vọt lố so với các thuật toán PID đơn vòng cổ điển.
  • Tính chính xác của mô hình nhận dạng tham số: Dữ liệu thực nghiệm chứng minh mô hình động cơ DC nhận dạng qua ARMA phản ánh chính xác 96,5% đặc tính động học thực tế, triệt tiêu hiện tượng trôi dạt điểm cân bằng do sai số quán tính cơ học.

Thảo luận kết quả

Nguyên nhân cốt lõi giúp hệ thống đạt hiệu năng vượt trội là sự phân tách hợp lý giữa bài toán phi tuyến diện rộng (Swing-up) và bài toán tuyến tính hóa quanh điểm làm việc (LQR). Phương pháp tuyến tính hóa hồi tiếp riêng phần đã tận dụng hiệu quả sự ghép kênh quán tính giữa thanh chủ động và thanh thụ động để tích lũy động năng, trong khi thuật toán LQR giải phương trình Riccati giúp tối ưu hóa điện áp cấp vào mạch cầu H, tránh hiện tượng bão hòa tín hiệu điều khiển.

So với các công trình nghiên cứu trước đây như giải thuật điều khiển dựa trên năng lượng của Isabelle Fantoni và Spong (vốn có thời gian dao động kéo dài trên 3,0 giây và miền hút cân bằng hẹp) hay nghiên cứu của Phan Việt Hùng sử dụng bộ điều khiển PID kinh điển (thời gian xác lập trên 2,5 giây với độ vọt lố lớn), bộ điều khiển lai đề xuất trong luận văn này đã rút ngắn thời gian đáp ứng hơn 35% và mở rộng đáng kể góc lệch khởi tạo an toàn. Dữ liệu thực nghiệm của nghiên cứu có thể được trực quan hóa sinh động thông qua biểu đồ đáp ứng thời gian thực của góc lệch q1, q2 từ 0 đến 10 giây, kết hợp bảng so sánh định lượng các chỉ số: thời gian quá độ, điện áp đỉnh tiêu thụ và sai số xác lập giữa các phương pháp điều khiển khác nhau.

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

Nhằm hoàn thiện và mở rộng khả năng ứng dụng thực tiễn của công trình nghiên cứu, bốn giải pháp trọng tâm được khuyến nghị thực hiện:

  • Tích hợp bộ lọc thích nghi Kalman mở rộng: Nhóm nghiên cứu tại các phòng thí nghiệm tự động hóa cần triển khai thuật toán lọc thích nghi để triệt tiêu tới 85% nhiễu đo lường từ cảm biến quang encoder và rung động cơ khí trong vòng 6 tháng tới.
  • Phát triển thuật toán điều khiển trượt mờ thích nghi (Adaptive Fuzzy Sliding Mode Control): Các nhà khoa học cơ điện tử cần kết hợp mạng nơ-ron nhân tạo với luật điều khiển mờ nhằm mở rộng biên độ góc lệch thích nghi lên trên 15 độ, thích ứng hoàn hảo khi thông số tải trọng thay đổi trong lộ trình 9 tháng.
  • Tối ưu hóa gia công cơ khí và giảm ma sát ổ trục: Xưởng cơ khí chính xác và đơn vị chế tạo robot cần nâng cấp chất lượng gia công cơ cấu khớp quay, sử dụng vật liệu hợp kim nhôm siêu nhẹ và ổ bi chất lượng cao để giảm ma sát tĩnh xuống dưới 5% trong thời hạn 3 tháng.
  • Hoàn thiện giao diện giám sát thời gian thực SCADA/HMI: Đội ngũ kỹ sư phần mềm nhúng cần xây dựng ứng dụng giám sát trực quan trên máy tính qua giao thức truyền thông SCI tốc độ cao, đạt tần số cập nhật dữ liệu 100 Hz phục vụ nghiên cứu và giảng dạy chuyên sâu trong 4 tháng tới.

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

Công trình luận văn thạc sĩ này là nguồn tài liệu học thuật và kỹ thuật vô cùng giá trị đối với các nhóm đối tượng sau:

  • Học viên cao học và nghiên cứu sinh chuyên ngành Kỹ thuật Điều khiển và Tự động hóa: Tiếp cận phương pháp luận chuẩn mực về mô hình hóa Euler-Lagrange, giải thuật tuyến tính hóa phản hồi và kỹ thuật tổng hợp bộ điều khiển tối ưu LQR cho hệ phi tuyến bậc cao.
  • Kỹ sư R&D thiết kế hệ thống Robot và Cơ điện tử: Ứng dụng nguyên lý điều khiển hệ thiếu dẫn động để phát triển các cánh tay robot công nghiệp thế hệ mới có trọng lượng nhẹ, giảm số lượng động cơ chấp hành mà vẫn đảm bảo độ chính xác chuyển động.
  • Giảng viên và cán bộ nghiên cứu tại các trường đại học kỹ thuật: Sử dụng toàn bộ cấu trúc bài toán, sơ đồ khối Simulink và kết quả thực nghiệm làm giáo trình giảng dạy chuyên đề điều khiển hiện đại và bài thí nghiệm mô hình thực hành.
  • Kỹ sư lập trình hệ thống nhúng và xử lý tín hiệu số: Khai thác quy trình nhận dạng thông số động cơ bằng ARMA và phương pháp lập trình tự động hóa nhúng mã nguồn từ MATLAB xuống dòng chip chuyên dụng DSP TMS320F28335.

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

  • Tại sao hệ Pendubot được phân loại vào nhóm hệ thống thiếu cơ cấu truyền động? Hệ Pendubot gồm 2 bậc tự do chuyển động trong mặt phẳng đứng nhưng chỉ có 1 động cơ tại khớp thứ nhất, khớp thứ hai quay tự do. Tỷ lệ cơ cấu chấp hành trên bậc tự do là 1:2. Thiết kế này giúp giảm khoảng 30% trọng lượng cơ khí nhưng tạo ra bài toán phi tuyến phức tạp do không thể tác động lực trực tiếp lên khâu thứ hai.

  • Vai trò cốt lõi của phương pháp tuyến tính hóa hồi tiếp riêng phần trong luận văn là gì? Phương pháp tuyến tính hóa hồi tiếp riêng phần được sử dụng trong giai đoạn Swing-up để chuyển đổi hệ phương trình phi tuyến phức tạp thành dạng tuyến tính tương đương đối với biến khớp xác định. Nhờ đó, hệ thống bơm năng lượng chuẩn xác và đưa con lắc lên lân cận điểm cân bằng chỉ trong khoảng 1,5 giây.

  • Bộ điều khiển LQR được cấu hình như thế nào để giữ cân bằng cho mô hình? Bộ điều khiển LQR hoạt động dựa trên việc cực tiểu hóa hàm chỉ tiêu chất lượng toàn phương với ma trận Q dạng đường chéo có các phần tử 10, 1, 10, 1 và trọng số R bằng 10. Thông qua việc giải phương trình đại số Riccati trên ma trận trạng thái tuyến tính hóa, hệ thống đạt được vector độ lợi K tối ưu với hệ số 110,1744.

  • Vi điều khiển DSP TMS320F28335 mang lại những ưu thế gì cho việc thực thi giải thuật? Chip DSP TMS320F28335 sở hữu xung nhịp 150 MHz, kiến trúc xử lý số thực 32-bit FPU và chu kỳ máy 6,67 ns, cho phép thực hiện các phép tính ma trận phức tạp theo thời gian thực. Khả năng liên kết trực tiếp với Simulink giúp nạp mã nhúng tự động, giảm thiểu 50% thời gian thử nghiệm phần cứng.

  • Làm thế nào để xác định các thông số vật lý của động cơ khi thiếu dữ liệu nhà sản xuất? Nghiên cứu áp dụng mô hình chuỗi thời gian ARMA kết hợp thuật toán bình phương tối thiểu để nhận dạng thông số từ 1.000 mẫu đo dòng điện, điện áp và tốc độ. Kết quả đã xác định chính xác điện trở Rm bằng 2,6824 Ohm và hằng số Kb bằng 0,0924 V.s/rad với sai số cực nhỏ.

Kết luận

  • Nghiên cứu đã xây dựng thành công mô hình toán học giải tích hai bậc tự do của hệ Pendubot dựa trên phương trình động lực học Euler-Lagrange.
  • Ứng dụng xuất sắc phương pháp ARMA và bình phương tối thiểu để nhận dạng 5 thông số vật lý then chốt của động cơ chấp hành với độ chính xác trên 96,5%.
  • Đề xuất và phối hợp hoàn hảo giải thuật điều khiển lai, rút ngắn thời gian Swing-up xuống dưới 1,5 giây và giữ ổn định góc lệch dưới 0,05 rad.
  • Chế tạo hoàn chỉnh mô hình phần cứng thực nghiệm, tích hợp mạch driver công suất và nhúng mã điều khiển thời gian thực lên DSP TMS320F28335 tốc độ 150 MHz.
  • Kiểm nghiệm toàn diện tính kháng nhiễu trong môi trường mô phỏng và thực tế, chứng minh hệ thống tự phục hồi trạng thái cân bằng trong thời gian dưới 0,8 giây khi có xung nhiễu 0,5 rad.

Đóng góp lớn nhất của luận văn là cung cấp giải pháp điều khiển lai hoàn chỉnh từ lý thuyết, mô phỏng đến hiện thực hóa phần cứng nhúng cho hệ thống thiếu cơ cấu truyền động. Trong lộ trình 6 đến 12 tháng tới, hướng phát triển trọng tâm là tích hợp mạng nơ-ron thích nghi và bộ lọc Kalman để triệt tiêu hoàn toàn rung động cơ khí. Quý độc giả, kỹ sư và nhà nghiên cứu quan tâm hãy tải toàn văn luận văn để áp dụng ngay giải thuật tối ưu này vào các dự án robot tiên tiến.