Tổng quan nghiên cứu

Trong xu thế tự động hóa công nghiệp hiện đại, các dòng tay máy công nghiệp đóng vai trò then chốt tại các dây chuyền lắp ráp linh kiện điện tử, đóng gói và phân loại sản phẩm. Robot SCARA (Selective Compliance Articulated Robot Arm) là cấu hình tay máy tiêu biểu được ứng dụng rộng rãi nhờ tốc độ gắp đạt từ 1,0 đến 2,0 m/s và độ lặp lại vị trí chính xác đến 0,1 mm. Tuy nhiên, thách thức lớn nhất trong việc điều khiển bám quỹ đạo chính xác cho robot SCARA nằm ở đặc tính phi tuyến mạnh, sự tương tác ghép kênh giữa các khớp quay và sự biến thiên quán tính khi mang tải.

Các phương pháp chỉnh định bộ điều khiển vi tích phân tỉ lệ (PID) truyền thống như Ziegler-Nichols hoặc thử sai thủ công thường tốn nhiều thời gian, dễ rơi vào cực trị địa phương và không đảm bảo chất lượng động học tối ưu cho các chuyển động bám quỹ đạo cong liên tục. Đề tài luận văn thạc sĩ tập trung giải quyết bài toán ứng dụng giải thuật di truyền (Genetic Algorithm - GA) để tìm kiếm tự động bộ thông số tối ưu gồm hệ số tỉ lệ Kp, hệ số tích phân Ki và hệ số vi phân Kd cho bộ điều khiển PID của robot SCARA hai bậc tự do.

Mục tiêu cụ thể của công trình là xây dựng mô hình toán học động học và động lực học của robot, thiết lập giải thuật tối ưu hóa GA ngoại tuyến trên nền tảng phần mềm MATLAB/Simulink nhằm bám quỹ đạo đường tròn và đường thẳng, sau đó chế tạo mô hình thực nghiệm để kiểm chứng độ tin cậy. Nghiên cứu được thực hiện tại Đại học Sư phạm Kỹ thuật Thành phố Hồ Chí Minh với việc chế tạo hoàn chỉnh phần cứng điều khiển sử dụng vi điều khiển PIC18F4431. Về mặt định lượng, nghiên cứu đặt mục tiêu kiểm soát độ vọt lố của hệ thống dưới ngưỡng 25% và duy trì thời gian xác lập dưới 4,5 giây, tạo tiền đề vững chắc cho việc ứng dụng các thuật toán tối ưu hóa thông minh vào các hệ thống tự động hóa công nghiệp thực tế.

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 ba trụ cột lý thuyết cốt lõi của ngành kỹ thuật điều khiển và robot học: lý thuyết động học và động lực học tay máy, lý thuyết điều khiển vi tích phân tỉ lệ (PID), và lý thuyết tính toán mềm với giải thuật di truyền (GA).

Cấu hình robot SCARA trong nghiên cứu thuộc kiểu hai khâu quay nối tiếp (RR) chuyển động trong mặt phẳng nằm ngang. Mô hình động học thuận được xác lập thông qua phương pháp Denavit-Hartenberg (D-H) để tính toán ma trận biến đổi tọa độ đồng nhất từ hệ tọa độ gốc sang hệ tọa độ khâu chấp hành cuối. Phương pháp động học ngược được phát triển thông qua giải tích hình học và đại số lượng giác nhằm giải nghiệm tường minh các góc quay khớp từ tọa độ không gian làm việc mong muốn.

Mô hình động lực học được thiết lập thông qua phương trình vi phân Lagrange loại 2 dựa trên hàm Lagrange tổng quát L = K - P, trong đó K là tổng động năng của hai khâu và tải trọng, còn P là thế năng của hệ thống. Phương trình vi phân mô tả mô-men xoắn tại từng trục động cơ bao gồm thành phần ma sát nhớt, quán tính tương hỗ và lực hướng tâm - Coriolis.

Bộ điều khiển PID kinh điển nhận tín hiệu sai lệch góc quay giữa quỹ đạo đặt và phản hồi thực tế để tạo ra tín hiệu điều khiển mô-men. Nhằm khắc phục nhược điểm của phương pháp chỉnh định tuyến tính, giải thuật di truyền (GA) được áp dụng như một công cụ tìm kiếm ngẫu nhiên có định hướng trên không gian giả thuyết đa chiều. Giải thuật mô phỏng quy luật tiến hóa tự nhiên thông qua việc mã hóa bộ tham số Kp, Ki, Kd thành các cá thể nhiễm sắc thể, thực hiện liên tục các toán tử chọn lọc dựa trên hàm thích nghi, lai ghép chéo với tỉ lệ xác định và đột biến gen ngẫu nhiên để tìm ra lời giải tối ưu toàn cục.

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 tính toán giải tích, mô hình hóa mô phỏng số và thực nghiệm phần cứng. Nguồn dữ liệu vật lý của robot bao gồm các thông số kết cấu xác định: chiều dài khâu một a1 = 0,13 m, chiều dài khâu hai a2 = 0,13 m, khối lượng khâu một m1 = 0,5 kg, khối lượng khâu hai m2 = 0,3 kg, khối lượng tải m = 0,05 kg, mô-men quán tính trục một I1 = 0,0028 kg.m² và mô-men quán tính trục hai I2 = 0,0017 kg.m², hệ số ma sát nhớt lần lượt là 1,0 kg.m²/s và 0,5 kg.m²/s.

Cỡ mẫu tính toán trong mô hình tối ưu hóa bao gồm quần thể 50 cá thể được khởi tạo ngẫu nhiên, trải qua 100 thế hệ tiến hóa lặp với tần số lấy mẫu mô phỏng đạt 1 kHz trong môi trường MATLAB. Phương pháp phân tích tối ưu được lựa chọn là giải thuật GA thay vì phương pháp đặt cực hay điều khiển mờ (Fuzzy Logic) truyền thống bởi vì hệ thống động lực học SCARA chứa các thành phần phi tuyến bậc cao và ghép kênh chéo phức tạp, vốn là môi trường mà GA thể hiện ưu thế vượt trội về khả năng tìm kiếm không gian rộng mà không đòi hỏi hàm truyền tuyến tính hóa quanh điểm làm việc.

Tiến trình nghiên cứu diễn ra qua ba giai đoạn liên tục trong khoảng thời gian 6 tháng: giai đoạn 1 giải bài toán động học và thiết lập mô hình toán học; giai đoạn 2 lập trình kịch bản M-file kết hợp khối mô phỏng Simulink để tối ưu hóa tham số PID; giai đoạn 3 thiết kế mạch in điều khiển và công suất trên phần mềm OrCAD, gia công cơ khí cánh tay robot, nạp mã nguồn điều khiển C cho vi điều khiển PIC18F4431 với chu kỳ ngắt điều khiển 10 ms và xây dựng giao diện giám sát HMI trên nền tảng Visual 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à thực nghiệm kiểm chứng giải thuật GA-PID trên robot SCARA đã mang lại bốn phát hiện khoa học và kỹ thuật nổi bật với các số liệu định lượng rõ ràng:

Thứ nhất, giải thuật di truyền đã hội tụ thành công để tìm ra bộ thông số tối ưu cho cả hai khớp quay: đối với khớp một thu được Kp1 = 20,4, Ki1 = 0,12, Kd1 = 1,85; đối với khớp hai thu được Kp2 = 18,6, Ki2 = 0,08, Kd2 = 1,42. Bộ thông số này giúp đáp ứng góc của cánh tay thỏa mãn hoàn toàn chỉ tiêu chất lượng với độ vọt lố đo được là 21,4%, thấp hơn 3,6% so với giới hạn trần 25% đặt ra ban đầu.

Thứ hai, khi kiểm tra đáp ứng của hệ thống với các dạng tín hiệu kích thích khác nhau, tín hiệu hàm bước (Step input) cho thấy khả năng đáp ứng vượt trội với thời gian xác lập đạt 3,0 giây và độ vọt lố giảm xuống còn 20,6%. Kết quả này cải thiện 33,3% về thời gian ổn định so với khi kiểm tra với tín hiệu xung vuông có chu kỳ chuyển trạng thái nhanh.

Thứ ba, mô phỏng bám quỹ đạo hình tròn với phương trình tham số bán kính r = 0,05 m dưới tác động của nguồn sóng sin cho thấy đầu công tác tái tạo biên dạng tròn trơn tru và chính xác ở bán kỳ dương, trong khi nguồn xung vuông có chu kỳ 2,0 giây gây ra hiện tượng méo quỹ đạo nghiêm trọng do hệ thống cần 4,5 giây để hoàn tất thời gian đáp ứng quá độ.

Thứ tư, trên mô hình thực nghiệm phần cứng sử dụng kit vi điều khiển PIC18F4431 kết hợp động cơ DC Servo và cơ cấu truyền động đai tỉ số 1:1, cánh tay robot đã vẽ thành công biên dạng quỹ đạo đường tròn và đường thẳng với sai lệch quỹ đạo hình học so với mô phỏng trên máy tính duy trì ở mức dưới 8,5%.

Thảo luận kết quả

Các dữ liệu nghiên cứu được biểu diễn trực quan thông qua đồ thị đáp ứng miền thời gian, giản đồ mô-men điều khiển khớp và biểu đồ quỹ đạo không gian 2D trên mặt phẳng làm việc Oxy. Biểu đồ mô-men khớp một cho thấy lực kích hoạt tăng cao ở mức 6,2 Nm trong 0,1 giây đầu tiên để thắng lực quán tính tĩnh, sau đó suy giảm và ổn định về giá trị 0 khi hệ thống đạt trạng thái cân bằng động học.

Nguyên nhân chính khiến hệ thống không thể bám tốt nguồn xung vuông chu kỳ 2,0 giây và nguồn xung tam giác chu kỳ 1,0 giây là do độ trễ cơ học và hằng số thời gian quán tính của cơ cấu vật lý. Khi tần số kích thích của tín hiệu đặt vượt quá băng thông đáp ứng quá độ của vòng lặp kín (4,5 giây), bộ điều khiển PID dù đã được tối ưu hóa tham số vẫn không thể tích lũy đủ năng lượng điều khiển để triệt tiêu sai số tức thời, dẫn đến hiện tượng dao động trễ pha.

So sánh với các nghiên cứu trước đây về điều khiển cánh tay một khâu mềm sử dụng Logic Mờ (Fuzzy Logic) và Đặt Cực (Pole Placement) có độ vọt lố từ 0 đến 2,1 độ dưới tải trọng 30 g, phương pháp GA-PID trong nghiên cứu này đã chứng minh được tính khả thi vượt trội khi giải quyết bài toán tay máy nhiều khâu có liên kết động học phi tuyến mà không đòi hỏi thiết kế tập luật mờ phức tạp hay giải hệ phương trình ma trận không gian trạng thái cồng kềnh.

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

Nhằm nâng cao hiệu quả vận hành và mở rộng khả năng ứng dụng thực tế của hệ thống robot SCARA điều khiển bằng giải thuật tối ưu hóa, bốn giải pháp kỹ thuật cụ thể được đề xuất:

  1. Nâng cấp bộ vi xử lý điều khiển từ vi điều khiển 8-bit PIC18F4431 lên vi điều khiển 32-bit kiến trúc ARM Cortex-M4 hoặc STM32F4: Tăng tần số lấy mẫu điều khiển từ 100 Hz (chu kỳ 10 ms) lên 1 kHz (chu kỳ 1 ms) nhằm giảm độ trễ vòng lặp điều chế độ rộng xung PWM, thực hiện trong vòng 3 tháng bởi nhóm kỹ sư phần cứng và hệ thống nhúng.
  2. Ứng dụng giải thuật di truyền thích nghi trực tuyến (Online Adaptive GA-PID): Cải thiện khả năng tự động cập nhật hệ số Kp, Ki, Kd theo thời gian thực khi tải trọng của robot biến thiên từ 0,05 kg lên 1,0 kg trong vòng dưới 0,5 giây, hoàn thiện trong lộ trình 6 tháng bởi nhóm nghiên cứu thuật toán điều khiển.
  3. Bổ sung thuật toán quy hoạch quỹ đạo bậc 5 (S-curve Trajectory Planning) và bộ lọc số thông thấp: Triệt tiêu hiện tượng giật động cơ học tại các điểm đảo chiều, hạ thấp độ vọt lố từ mức 21,4% hiện tại xuống dưới 5,0%, triển khai trong thời gian 4 tháng bởi các chuyên gia kỹ thuật cơ điện tử.
  4. Tối ưu hóa phần mềm giao diện giám sát SCADA trên Visual C# và tích hợp chuẩn truyền thông công nghiệp Modbus/RS485: Nâng tần suất cập nhật dữ liệu giám sát góc quay và tọa độ vị trí lên 50 Hz, giảm độ trễ hiển thị xuống dưới 20 ms, thực hiện trong thời hạn 2 tháng bởi nhóm phát triển phần mềm tự động hóa.

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

Công trình nghiên cứu mang giá trị học thuật và ứng dụng thực tiễn cao, đặc biệt phù hợp cho bốn nhóm đối tượng chuyên môn:

Nhóm thứ nhất là sinh viên, học viên cao học chuyên ngành Kỹ thuật Điều khiển và Tự động hóa, Cơ điện tử và Kỹ thuật Robot. Tài liệu cung cấp quy trình hoàn chỉnh từ giải tích động học Denavit-Hartenberg, thiết lập mô hình Lagrange đến mã nguồn mẫu MATLAB và chương trình vi điều khiển PIC.

Nhóm thứ hai là các kỹ sư nghiên cứu và phát triển (R&D) tại các nhà máy sản xuất thiết bị tự động hóa và chế tạo robot công nghiệp. Họ có thể khai thác trực tiếp mô hình toán học và bộ thông số GA-PID để rút ngắn từ 30% đến 40% thời gian hiệu chỉnh các dòng tay máy gắp đặt sản phẩm.

Nhóm thứ ba là giảng viên và nghiên cứu viên tại các trường đại học khối kỹ thuật. Luận văn là tài liệu tham khảo đắc lực cho các học phần Robot công nghiệp, Điều khiển thông minh và Thực tập hệ thống nhúng với đầy đủ bản vẽ 3D và sơ đồ nguyên lý mạch.

Nhóm thứ tư là các doanh nghiệp khởi nghiệp và đơn vị phát triển giải pháp IoT công nghiệp. Nhóm này có thể ứng dụng kiến trúc kết nối giữa phần mềm giám sát Visual C# và mạch điều khiển phần cứng để xây dựng các trạm làm việc tự động hóa chi phí thấp.

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

Ưu điểm chính của việc dùng giải thuật di truyền (GA) để tìm thông số PID trong luận văn là gì? Giải thuật GA cung cấp khả năng tìm kiếm tối ưu toàn cục trên không gian tham số rộng mà không yêu cầu tuyến tính hóa mô hình hay phụ thuộc vào cấu trúc hàm truyền phức tạp. Phương pháp này giúp hệ thống đạt độ vọt lố 21,4% và bám chính xác quỹ đạo mà không cần thử sai thủ công.

Giới hạn lớn nhất của giải thuật di truyền trong nghiên cứu này được xác định là gì? Hạn chế lớn nhất của GA là tốc độ tính toán tiến hóa ngoại tuyến tốn nhiều tài nguyên máy tính với thời gian chạy mô phỏng qua 100 thế hệ. Do đó, thuật toán mới chỉ dừng lại ở mức tìm kiếm tham số offline trên MATLAB trước khi nạp trực tiếp giá trị vào vi điều khiển.

Phần cứng điều khiển thực tế của robot SCARA được xây dựng dựa trên nền tảng nào? Mô hình thực nghiệm sử dụng kit vi điều khiển PIC18F4431 tích hợp khối phát xung PWM chuyên dụng cho điều khiển động cơ, kết hợp mạch công suất điều khiển DC Servo, hệ thống cảm biến góc Encoder quang học và giao tiếp dữ liệu với máy tính qua cổng nối tiếp UART.

Tại sao cánh tay robot không thể bám chính xác theo nguồn tín hiệu xung vuông? Do thời gian đáp ứng quá độ của hệ thống vòng kín là 4,5 giây, trong khi chu kỳ kích thích của xung vuông chỉ có 2,0 giây. Tốc độ thay đổi trạng thái quá nhanh của tín hiệu đặt khiến bộ truyền động cơ khí không kịp tăng tốc, gây méo biên dạng quỹ đạo di chuyển.

Kích thước hình học và khối lượng của mô hình robot SCARA trong luận văn là bao nhiêu? Cánh tay robot hai bậc tự do có chiều dài khâu một và khâu hai bằng nhau là 0,13 m, tổng tầm với tối đa đạt 0,26 m. Khối lượng khâu một là 0,5 kg, khâu hai là 0,3 kg và mang tải trọng định mức tại khâu công tác là 0,05 kg.

Kết luận

  • Đề tài đã xây dựng thành công mô hình toán học giải tích toàn diện bao gồm động học thuận, động học ngược và động lực học Lagrange cho robot SCARA hai khâu quay.
  • Ứng dụng xuất sắc giải thuật di truyền GA trên nền tảng MATLAB/Simulink để tìm ra bộ thông số PID tối ưu (Kp, Ki, Kd), đưa độ vọt lố hệ thống về mức 21,4% và thời gian đáp ứng bước chỉ còn 3,0 giây.
  • Chế tạo hoàn chỉnh mô hình thực nghiệm cơ điện tử tích hợp vi điều khiển PIC18F4431, mạch công suất DC Servo và giao diện giám sát HMI chuyên nghiệp trên Visual C#.
  • Kiểm chứng thực tế sự bám sát giữa mô phỏng và chuyển động thực nghiệm đối với các quỹ đạo chuẩn như đường tròn bán kính 0,05 m và đường thẳng.
  • Đóng góp giải pháp khoa học kết hợp giữa trí tuệ nhân tạo và điều khiển kinh điển, mở ra hướng nghiên cứu thích nghi trực tuyến và quy hoạch quỹ đạo nâng cao trong 6 tháng tới; các nhà nghiên cứu và kỹ sư quan tâm có thể ứng dụng ngay bộ tham số và sơ đồ mạch để phát triển các thế hệ robot công nghiệp tiếp theo.