Tổng quan nghiên cứu

Trong bối cảnh cuộc cách mạng công nghiệp lần thứ tư, các hệ thống cơ điện tử thiếu cơ cấu chấp hành (under-actuated) với đặc tính phi tuyến cao đóng vai trò then chốt trong việc kiểm chứng và phát triển các thuật toán điều khiển hiện đại. Hệ thống bánh quay trên con lắc ngược tự do (Reaction Wheel Inverted Pendulum - RWIP) là một mô hình thực nghiệm tiêu biểu với khả năng quay tự do 360 độ quanh trục cố định. Đối tượng này chịu tác động phức tạp từ trọng lực, lực ly tâm cùng lực ma sát lăn tại ổ bi KBC 608 RS với hệ số ma sát khoảng 0,12 x 10^-3 và hệ số ma sát động cơ khoảng 0,5. Vấn đề cốt lõi của bài toán là con lắc luôn có xu hướng rơi tự do khỏi vị trí cân bằng chỉ trong khoảng 0,3 giây nếu không có mô-men quán tính phản kháng kịp thời từ bánh đà. Việc duy trì con lắc đứng thẳng trong phạm vi góc hẹp từ -15 độ đến +15 độ đòi hỏi giải thuật phải đạt độ chính xác cao và triệt tiêu nhanh các dao động góc.

Luận văn thạc sĩ chuyên ngành Kỹ thuật Điều khiển và Tự động hóa (mã ngành: 8520216) của tác giả Nguyễn Bình Hậu, dưới sự hướng dẫn khoa học của PGS.TS Nguyễn Minh Tâm tại Trường Đại học Sư phạm Kỹ thuật TP. Hồ Chí Minh (bảo vệ tháng 4/2019), đã giải quyết trọn vẹn thách thức trên. Mục tiêu cụ thể của công trình là thiết lập hệ phương trình toán học phi tuyến toàn diện, xây dựng và so sánh hiệu quả của các bộ điều khiển từ kinh điển (PID, LQR) đến thông minh (Mamdani Fuzzy, Sugeno Fuzzy) kết hợp thuật toán tối ưu hóa tiến hóa (GA, PSO), đồng thời hiện thực hóa thuật toán trên phần cứng nhúng thời gian thực STM32F407 Discovery. Nghiên cứu mang lại ý nghĩa học thuật và ứng dụng to lớn khi rút ngắn thời gian xác lập hệ thống xuống dưới 1,5 giây, giảm sai số góc lắc về mức dưới 0,05 radian và tiết kiệm hơn 25% năng lượng điều khiển động cơ DC.

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 xây dựng trên nền tảng cơ học giải tích Euler-Lagrange với hàm Lagrange tổng quát L = T - V, trong đó T là tổng động năng của thanh lắc cùng bánh quay và V là tổng thế năng dưới trường trọng lực gia tốc 9,5 m/s^2. Phương trình trạng thái phi tuyến của hệ thống được thiết lập thông qua 4 biến trạng thái: góc lệch thanh lắc, vận tốc góc thanh lắc, góc quay bánh đà và vận tốc góc bánh đà. Về mặt điều khiển tối ưu, lý thuyết điều khiển tuyến tính toàn phương (LQR) được áp dụng dựa trên việc giải phương trình vi phân đại số Riccati để tìm ma trận khuếch đại hồi tiếp K nhằm cực tiểu hóa hàm mục tiêu chất lượng J kết hợp giữa sai số trạng thái và năng lượng điều khiển.

Bên cạnh đó, khung lý thuyết điều khiển thông minh ứng dụng logic mờ Mamdani và Sugeno với công cụ mạng nơ-ron mờ thích nghi ANFIS (Adaptive Neuro-Fuzzy Inference System). Cấu trúc ANFIS kết hợp khả năng suy diễn tường minh của luật mờ If-Then với khả năng tự học trọng số của mạng nơ-ron nhiều lớp thông qua thuật toán lan truyền ngược kết hợp bình phương cực tiểu. Để giải quyết bài toán tối ưu hóa đa biến phi tuyến, đề tài tích hợp giải thuật di truyền (GA) với cơ chế chọn lọc tuyến tính, lai ghép BLX-alpha và giải thuật tối ưu hóa bầy đàn (PSO) dựa trên cập nhật vận tốc và vị trí của các hạt trong không gian tìm kiếm đa chiều.

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

Phương pháp nghiên cứu kết hợp chặt chẽ giữa mô hình hóa lý thuyết, mô phỏng số học trên máy tính và kiểm chứng thực nghiệm phần cứng. Nguồn dữ liệu phục vụ nghiên cứu bao gồm bộ thông số vật lý thực đo đạc trực tiếp từ mô hình cơ khí: khối lượng thanh lắc 0,92856 kg, khối lượng bánh đà 1,97856 kg, mô-men quán tính thanh lắc 6,9642 x 10^-3 kg.m^2, mô-men quán tính rotor 2,4732 x 10^-3 kg.m^2, điện trở phần ứng 6,5 Ohm và điện cảm 4,5 mH. Mô hình 3D được thiết kế chính xác trên phần mềm Autodesk Inventor 2017 để trích xuất trọng tâm và mô-men quán tính hình học.

Quá trình mô phỏng thực hiện trên môi trường Matlab/Simulink 2015 với hơn 1.000 mẫu dữ liệu động học phục vụ huấn luyện mạng ANFIS và 50 thế hệ tiến hóa cho các thuật toán GA và PSO nhằm đảm bảo tìm kiếm được nghiệm hội tụ toàn cục. Lý do lựa chọn bộ công cụ Waijung Blockset và vi điều khiển STM32F407ARM Cortex-M4 tốc độ 168 MHz là khả năng biên dịch trực tiếp từ sơ đồ khối Simulink sang mã nguồn C nhúng thời gian thực với tần số lấy mẫu 1 kHz. Quy trình thực nghiệm được triển khai lặp lại 10 lần độc lập trên mô hình vật lý tại phòng thí nghiệm tự động hóa trong giai đoạn từ năm 2017 đến đầu năm 2019 để thu thập dữ liệu dòng điện, điện áp cầu H và phản hồi góc từ 2 bộ encoder quang học độ phân giải cao.

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

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

Kết quả nghiên cứu thực nghiệm và mô phỏng đã chứng minh sự vượt trội rõ rệt của các giải thuật điều khiển thông minh so với các phương pháp truyền thống:

Thứ nhất, đối với bộ điều khiển PID, việc ứng dụng giải thuật di truyền GA để tìm kiếm bộ ba tham số Kp, Ki, Kd giúp giảm độ vọt lố tới 38% và rút ngắn thời gian xác lập góc lắc từ 2,8 giây (theo phương pháp Ziegler-Nichols cổ điển) xuống còn 1,6 giây, đồng thời biên độ dao động điện áp cấp cho động cơ DC giảm khoảng 30%.

Thứ hai, trong hệ thống điều khiển LQR tối ưu, giải thuật bầy đàn PSO cho thấy tốc độ hội tụ hàm mục tiêu J nhanh hơn GA khoảng 25% số thế hệ, giúp hệ thống đạt trạng thái cân bằng ổn định chỉ sau 1,2 giây với điện áp điều khiển duy trì mượt mà trong ngưỡng an toàn từ -12V đến +12V.

Thứ ba, việc cải tiến bộ điều khiển mờ Sugeno từ cấu trúc 81 luật ban đầu xuống bộ luật thỏa hiệp 2x9 luật thông qua mạng ANFIS đã giảm tải hơn 75% khối lượng tính toán của bộ xử lý, trong khi vẫn duy trì sai số bám góc tĩnh cực nhỏ dưới 0,02 radian.

Thứ tư, thử nghiệm thời gian thực trên bo mạch nhúng STM32F407 khẳng định mô hình phần cứng hoạt động ổn định, con lắc giữ vững vị trí thẳng đứng ngay cả khi chịu tác động của lực nhiễu xung cơ học tác động từ bên ngoài.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự cải thiện vượt bậc của các giải thuật thông minh là khả năng thích ứng linh hoạt với các thành phần phi tuyến khó mô hình hóa chính xác như ma sát trượt trong động cơ và hiện tượng bão hòa điện áp phần ứng. Bộ điều khiển LQR kết hợp PSO tối ưu hóa phân bố ma trận trọng số Q và R, tạo ra lực quán tính rotor cân bằng chính xác với mô-men trọng lực tức thời. So với các nghiên cứu trước đây vốn chỉ khảo sát trên mô hình con lắc ngược xe đẩy (Cart-Pole), hệ RWIP phức tạp hơn do cơ cấu chấp hành gắn trực tiếp trên đỉnh thanh lắc tự do.

Dữ liệu thực nghiệm được trực quan hóa qua biểu đồ đáp ứng thời gian của góc thanh lắc và vận tốc góc bánh đà trên màn hình giám sát Terminal. Biểu đồ dạng sóng thể hiện rõ đường cong góc lắc dưới tác động của bộ điều khiển Mamdani-GA triệt tiêu hoàn toàn hiện tượng rung lắc bậc cao (chattering) thường gặp ở các bộ điều khiển tuyến tính kinh điển. Bảng so sánh chỉ tiêu chất lượng động học chứng minh rằng phương pháp Sugeno ANFIS 9 luật đạt mức tiêu hao năng lượng trung bình thấp nhất, mở ra giải pháp tối ưu cho các thiết bị điều khiển hạn chế về dung lượng pin.

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

Dựa trên các kết quả đạt được từ luận văn thạc sĩ, một số giải pháp kỹ thuật cụ thể được đề xuất nhằm hoàn thiện và mở rộng khả năng ứng dụng thực tế của hệ thống:

Thứ nhất, nâng cấp phần cứng cơ cấu chấp hành bằng việc thay thế động cơ DC chổi than sang động cơ DC không chổi than (BLDC) công suất cao tích hợp cảm biến góc tuyệt đối 16-bit và khối cảm biến quán tính IMU 9 trục nhằm giảm thiểu độ trễ thu thập tín hiệu xuống dưới 1 mili-giây, thời gian triển khai trong vòng 6 tháng bởi nhóm kỹ sư cơ điện tử.

Thứ hai, phát triển thuật toán điều khiển trượt thích nghi (Adaptive Sliding Mode Control) kết hợp mạng nơ-ron hàm cơ sở xuyên tâm (RBF) nhằm tăng cường độ bền vững của hệ thống lên 20% đến 30% khi đối mặt với tải trọng thay đổi liên tục, thực hiện trong 12 tháng do nhóm nghiên cứu thuật toán phụ trách.

Thứ ba, chuẩn hóa và tối ưu hóa mã nguồn C/C++ nhúng trực tiếp trên vi điều khiển kiến trúc ARM 32-bit nhằm giảm mức chiếm dụng bộ nhớ RAM xuống dưới 40%, đáp ứng tiêu chuẩn thời gian thực khắt khe trong quý 3 năm tới bởi các chuyên gia phần mềm nhúng.

Thứ tư, mở rộng thuật toán từ chế độ điều khiển cân bằng cục bộ sang bài toán tự vung (Swing-up) toàn dải 360 độ từ vị trí buông thõng ban đầu lên vị trí đỉnh thẳng đứng trong vòng 9 tháng, giao cho các học viên cao học và nghiên cứu sinh thực hiện.

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

Công trình nghiên cứu mang tính hệ thống cao từ lý thuyết đến thực nghiệm, đặc biệt hữu ích cho 4 nhóm đối tượng trọng tâm sau:

Thứ nhất, học viên cao học và nghiên cứu sinh ngành Kỹ thuật Điều khiển và Tự động hóa, Cơ điện tử và Robot: Tài liệu cung cấp quy trình mẫu chuẩn mực về giải tích động lực học Lagrange, thiết kế hệ mờ lai nơ-ron ANFIS và lập trình giải thuật tối ưu hóa GA, PSO.

Thứ hai, kỹ sư phát triển hệ thống nhúng và điều khiển chuyển động: Cung cấp tài liệu thực chiến về cách tích hợp môi trường Matlab/Simulink với phần cứng STM32F4 thông qua thư viện Waijung Blockset để hiện thực hóa các thuật toán phức tạp trên chip 32-bit.

Thứ ba, giảng viên và các nhà nghiên cứu tại các viện, trường đại học kỹ thuật: Có thể sử dụng toàn bộ mô hình toán học và sơ đồ khối điều khiển làm học liệu giảng dạy chuyên đề Điều khiển nâng cao và thí nghiệm điều khiển hệ phi tuyến.

Thứ tư, chuyên gia phát triển công nghệ hàng không vũ trụ và thiết bị tự cân bằng: Tham khảo nguyên lý ổn định hướng bằng bánh đà phản lực (Reaction Wheels) để ứng dụng trên vệ tinh nhỏ CubeSat, kính thiên văn không gian và robot thám hiểm hai bánh.

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

Hệ thống bánh quay trên con lắc ngược (RWIP) có điểm gì khác biệt so với hệ con lắc ngược xe đẩy truyền thống?
Điểm khác biệt căn bản là hệ RWIP sử dụng mô-men xoắn quán tính sinh ra từ gia tốc góc của bánh đà đặt trên đỉnh thanh lắc để tạo lực giữ cân bằng, thay vì dùng lực kéo tuyến tính của xe đẩy di chuyển qua lại trên ray. Điều này làm cho hệ RWIP có tính linh hoạt cao hơn, kích thước nhỏ gọn hơn nhưng mức độ bất định và phi tuyến phức tạp hơn nhiều.

Tại sao giải thuật GA lại tối ưu hóa bộ điều khiển PID hiệu quả hơn phương pháp Ziegler-Nichols?
Phương pháp Ziegler-Nichols dựa trên việc xấp xỉ tuyến tính tại điểm làm việc tĩnh nên dễ gây dao động lớn và bão hòa điện áp. Ngược lại, giải thuật di truyền GA tìm kiếm trên toàn bộ không gian tham số đa chiều thông qua toán tử chọn lọc và lai ghép, giúp giảm 38% độ vọt lố và đáp ứng tối ưu phi tuyến thực tế.

Lợi ích của việc rút gọn bộ mờ Sugeno từ 81 luật xuống bộ 9 luật thông qua mạng ANFIS là gì?
Bộ 81 luật đòi hỏi năng lực xử lý ma trận lớn gây trễ chu kỳ lấy mẫu trên vi điều khiển nhúng. Mạng ANFIS học các đặc tính cốt lõi và tối ưu hóa lại thành 2 khối mờ 9 luật, giúp tiết kiệm hơn 75% thời gian tính toán của CPU STM32F4 mà vẫn đảm bảo độ chính xác bám góc với sai số dưới 0,02 radian.

Thuật toán tối ưu bầy đàn (PSO) hỗ trợ thiết kế bộ điều khiển LQR như thế nào?
Trong bộ điều khiển LQR, việc chọn ma trận trọng số Q và R thường dựa trên thử sai tốn nhiều thời gian. Giải thuật PSO tự động hóa quá trình này bằng cách phân tán bầy hạt tìm kiếm tổ hợp tham số tối ưu, giúp hàm mục tiêu J hội tụ nhanh hơn GA khoảng 25% và đạt chất lượng động học vượt trội.

Thư viện Waijung Blockset trên chip STM32F4 mang lại ưu thế gì cho quá trình thực nghiệm?
Waijung Blockset cho phép người dùng giao tiếp phần cứng trực tiếp từ sơ đồ Simulink và tự động dịch thành mã nhúng nạp vào chip ARM Cortex-M4 168 MHz trong vài giây. Điều này loại bỏ hoàn toàn các lỗi lập trình thủ công bằng tay, rút ngắn thời gian phát triển và đảm bảo tính đồng nhất giữa mô phỏng và thực tế.

Kết luận

Luận văn thạc sĩ của tác giả Nguyễn Bình Hậu đã giải quyết thành công bài toán điều khiển cân bằng cho một đối tượng phi tuyến cao và phức tạp trong kỹ thuật tự động hóa:

  • Thiết lập hệ phương trình toán học phi tuyến Euler-Lagrange hoàn chỉnh cho cơ hệ RWIP tổng khối lượng 2,907 kg có xét đến đầy đủ các thành phần ma sát cơ học và điện từ.
  • Xây dựng và kiểm chứng thành công các bộ điều khiển kinh điển PID, LQR kết hợp kỹ thuật tối ưu hóa toàn cục tự động GA và PSO.
  • Ứng dụng xuất sắc công cụ ANFIS để phát triển bộ điều khiển mờ Sugeno 2x9 luật tinh gọn, giảm 75% chi phí tính toán phần cứng.
  • Hoàn thiện mô hình thực nghiệm thời gian thực trên nền tảng vi điều khiển 32-bit STM32F407 kết hợp mạch công suất cầu H.
  • Đạt chỉ tiêu chất lượng động học xuất sắc với thời gian xác lập dưới 1,5 giây và sai số góc lắc duy trì ổn định dưới 0,05 radian.

Trong giai đoạn 6 đến 12 tháng tới, hướng nghiên cứu tiếp theo sẽ tập trung vào giải thuật tự vung lắc 360 độ và tích hợp bộ lọc Kalman mở rộng. Độc giả quan tâm và các nhà nghiên cứu trong lĩnh vực điều khiển tự động được khuyến khích tham khảo toàn văn luận văn để ứng dụng trực tiếp các mô hình và giải thuật này vào công tác nghiên cứu khoa học cũng như các dự án công nghệ thực tiễn.