Tổng quan nghiên cứu

Twin Rotor MIMO System (TRMS) do Feedback Instruments Ltd chế tạo từ năm 1998 là mô hình khí động học 2 bậc tự do (2-DOF) gồm 4 biến trạng thái chuyển động cơ bản ($\alpha_v, \alpha_h, \omega_v, \omega_h$). Thiết bị này mô phỏng sát thực tế động học bay của máy bay trực thăng với 2 cánh quạt nâng chính và cân bằng đuôi được dẫn động độc lập bởi 2 động cơ điện một chiều. Vấn đề nghiên cứu trọng tâm xuất phát từ đặc tính phi tuyến cực mạnh, tính bất định tham số ma sát và quán tính, cùng hiện tượng xen kênh chéo nghiêm trọng giữa chuyển vị góc đứng (pitch angle $\alpha_v$) và góc ngang (yaw angle $\alpha_h$). Khi tốc độ động cơ đuôi thay đổi để quay thân máy bay theo phương ngang, lực khí động học lập tức tác động làm thay đổi góc nâng đứng và ngược lại.

Mục tiêu cụ thể của luận văn là thiết kế thành công cấu trúc điều khiển nâng cao kết hợp giữa mạch vòng phản hồi PID tối ưu hóa bằng giải thuật di truyền (GA) và mạch vòng truyền thẳng học thích nghi (LFFC) trên cơ sở hệ thích nghi theo mô hình mẫu (MRAS) theo phương pháp bù tổng. Đề tài giữ nguyên bản chất phi tuyến của đối tượng thay vì tuyến tính hóa cục bộ. Nghiên cứu được triển khai thực nghiệm tại Viện nghiên cứu phát triển công nghệ cao về kỹ thuật công nghiệp và Đại học Kỹ thuật Công nghiệp – Đại học Thái Nguyên, hoàn thành nghiệm thu vào năm 2017. Kết quả nghiên cứu có ý nghĩa lớn trong việc nâng cao độ chính xác bám quỹ đạo, giảm sai số toàn phương trung bình (MSE) xuống dưới mức 0.006 và triệt tiêu tương tác xen kênh thời gian thực cho các hệ thống nhiều đầu vào nhiều đầu ra (MIMO).

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

Khung lý thuyết áp dụng

Khung lý thuyết của công trình được xây dựng trên 3 trụ cột toán học và điều khiển học nâng cao:

  1. Mô hình động học Euler - Lagrange phi tuyến toàn phần: Khác với phương pháp Euler - Newton tuyến tính hóa bậc 1 thường gặp, mô hình Euler - Lagrange thiết lập phương trình vi phân chuyển động thông qua phân tích biến thiên động năng và thế năng của 3 cụm kết cấu: thanh dầm tự do, cần đối trọng và trục quay thẳng đứng. Mô hình này bảo toàn 100% các thành phần mô men ma sát, mô men hồi chuyển $M_{gyro}$ và mô men cáp $M_{cable}$.

  2. Tối ưu hóa tham số PID bằng giải thuật di truyền (GA): Dựa trên lý thuyết tiến hóa của Holland (1975), thuật toán GA mô phỏng quá trình chọn lọc tự nhiên để tìm kiếm bộ 3 tham số ($K_p, K_i, K_d$) tối ưu toàn cục. Quá trình xử lý trải qua 4 bước: mã hóa chuỗi nhị phân, khởi tạo quần thể, lai ghép đột biến và chọn lọc cá thể theo 3 hàm mục tiêu sai số gồm IAE (Integral of Absolute Error), ITAE (Integral of Time Absolute Error) và MSE (Mean Squared Error).

  3. Bộ điều khiển truyền thẳng học thích nghi (LFFC trên cơ sở MRAS): Kế thừa lý thuyết điều khiển thích nghi mô hình mẫu của Whitaker (1958) và nguyên lý ổn định tiệm cận Lyapunov của Park (1966). Khâu truyền thẳng thích nghi đóng vai trò nhận dạng mô hình ngược động học của TRMS để sinh tín hiệu điều khiển truyền thẳng $u_F$, trong khi mạch vòng phản hồi PID xử lý sai lệch tĩnh và cung cấp tín hiệu học trực tuyến.

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

Nghiên cứu sử dụng nguồn dữ liệu thực nghiệm thu thập trực tiếp từ hệ thống phần cứng TRMS kết nối với card xử lý tín hiệu số dSPACE DS1103, vận hành đồng bộ qua phần mềm MATLAB/Simulink và giao diện ControlDesk.

Phương pháp chọn mẫu dữ liệu được thực hiện theo nguyên tắc lấy mẫu ngẫu nhiên phân tầng theo dải góc hoạt động. Cỡ mẫu đạt 10000 điểm dữ liệu trong mỗi chu kỳ đo 10 giây với tần số lấy mẫu cố định 1000 Hz (chu kỳ trích mẫu 1 ms). Dải góc khảo sát bao phủ từ âm 60 độ đến dương 60 độ trên mặt đứng và toàn dải 360 độ trên mặt ngang nhằm kiểm chứng đầy đủ mọi vùng động học phi tuyến.

Lý do lựa chọn phương pháp phân tích miền thời gian kết hợp đánh giá chỉ số sai số MSE là nhằm định lượng chính xác sai lệch bám quỹ đạo và tốc độ dập tắt dao động. Để trích xuất đạo hàm sai lệch an toàn mà không làm khuếch đại nhiễu đo lường, một bộ lọc biến trạng thái (SVF) bậc 2 được tích hợp với dải thông thiết lập gấp 10 lần dải thông của mô hình mẫu tham chiếu. Toàn bộ tiến trình nghiên cứu kéo dài qua 4 giai đoạn logic: thiết lập phương trình vi phân giải tích, mô phỏng kiểm chứng trên Simulink, tối ưu hóa thông số qua GA và nạp mã thực thi thời gian thực trên phần cứng dSPACE DS1103.

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 số và thử nghiệm thời gian thực trên thiết bị TRMS mang lại 3 phát hiện khoa học cốt lõi:

  1. Độ chính xác vượt trội của mô hình Euler - Lagrange: Mô hình động học giải tích phản ánh chính xác đáp ứng vật lý của hệ thống với sai số toàn phương trung bình MSE của góc trệch ngang ($\alpha_h$) đạt 0.0053 và góc chao nghiêng đứng ($\alpha_v$) đạt 0.0041. Mức sai số này thấp hơn đáng kể so với các mô hình tuyến tính hóa thông thường.

  2. Hiệu năng tối ưu hóa tham số của GA: Từ 6 giá trị khởi tạo được xác định sơ bộ theo phương pháp Ziegler - Nichols, giải thuật GA nhanh chóng tìm ra bộ tham số PID tối ưu cho cả 2 kênh rotor. Bộ điều khiển PID-GA giúp giảm độ quá điều chỉnh từ 24.5% xuống dưới 4.6%, đồng thời giảm chỉ số tích phân sai số ITAE hơn 42%.

  3. Năng lực khử xen kênh của cấu trúc LFFC-MRAS: Mạch vòng truyền thẳng học thích nghi tạo ra tín hiệu bù tổng chính xác, triệt tiêu hơn 78% tác động ghép kênh chéo giữa 2 rotor. Thời gian xác lập của góc $\alpha_v$ và $\alpha_h$ được rút ngắn hơn 36% so với bộ điều khiển PID độc lập, đưa sai số xác lập về tiệm cận 0 ngay cả khi chịu tác động của nhiễu ngẫu nhiên.

Thảo luận kết quả

Hiệu quả vượt trội của cấu trúc đề xuất xuất phát từ cơ chế bù tổng động học ngược. Thay vì cố gắng bóc tách từng kênh riêng biệt gây phức tạp cho thuật toán, bộ điều khiển LFFC xem toàn bộ tương tác phi tuyến và xen kênh là một thành phần nhiễu hệ thống tổng hợp. Bằng việc giải phương trình đại số Lyapunov với ma trận xác định dương $P$ và hệ số tăng tốc $\lambda = 10$, các tham số thích nghi ($a_m, b_m, c_m, d_m$) hội tụ nhanh chóng, tạo ra đáp ứng bù tức thời cho điện áp cấp vào 2 động cơ DC.

Dữ liệu thực nghiệm được trình bày trực quan thông qua các đồ thị đáp ứng thời gian thực trên phần mềm ControlDesk. Biểu đồ so sánh đường cong đáp ứng góc $\alpha_v$ và $\alpha_h$ giữa phương án PID đơn thuần và phương án kết hợp PID-GA với LFFC-MRAS cho thấy sự biến mất hoàn toàn của các dao động quá độ kéo dài. Bảng tổng hợp các chỉ số định lượng (độ quá điều chỉnh, thời gian quá độ, chỉ số IAE và MSE) khẳng định rằng cấu trúc bù tổng thích nghi duy trì tính ổn định tiệm cận của hệ thống tốt hơn hẳn so với phương pháp tách kênh truyền thống hay điều khiển deadbeat từng được công bố trong các nghiên cứu trước đây.

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

Dựa trên kết quả nghiên cứu lý thuyết và thực nghiệm, tác giả đưa ra 4 khuyến nghị kỹ thuật cụ thể:

  1. Chuẩn hóa thuật toán ước lượng ma sát và nhiệt độ trực tuyến: Các kỹ sư điều khiển cần lập trình tích hợp khối nhận dạng thích nghi thời gian thực trên bo mạch DSP ở tần số 1000 Hz. Giải pháp này giúp tự động bù trôi dạt thông số ma sát và điện trở cuộn dây khi động cơ DC hoạt động liên tục trên 45 phút, giữ sai số bám góc ổn định dưới 0.8 độ.

  2. Nâng cấp giải thuật GA sang cơ chế tự chỉnh trực tuyến (Online Adaptive GA): Nhóm nghiên cứu phần mềm nhúng nên triển khai thuật toán GA song song trên các dòng vi xử lý đa lõi trong lộ trình 6 tháng tới. Mục tiêu là tái cấu hình bộ 6 tham số PID trong thời gian thực dưới 50 ms khi đối trọng hoặc tải trọng thay đổi đột ngột.

  3. Tối ưu hóa băng thông bộ lọc biến trạng thái SVF: Đơn vị vận hành phòng thí nghiệm cần duy trì dải thông SVF trong khoảng từ 10 đến 12 lần dải thông mô hình mẫu. Việc này đảm bảo khả năng lọc sạch xung nhiễu tần số cao từ cảm biến encoder quang học mà không gây trễ pha tín hiệu điều khiển.

  4. Chuyển giao công nghệ điều khiển bù tổng cho thiết bị bay không người lái (UAV): Các viện nghiên cứu hàng không và doanh nghiệp công nghệ cao nên ứng dụng cấu trúc PID-GA kết hợp LFFC-MRAS lên trực thăng tự hành và drone 4 cánh trong giai đoạn từ 9 đến 12 tháng tới, hướng tới cải thiện 30% độ ổn định bay treo khi gặp gió giật cấp 4.

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

Công trình luận văn mang lại giá trị học thuật và thực tiễn cao cho 4 nhóm đối tượng:

  1. Kỹ sư nghiên cứu và phát triển khí cụ bay (UAV/Drone): Tiếp cận giải pháp xử lý triệt để hiện tượng ghép kênh khí động học 2 trục mà không cần tuyến tính hóa mô hình, giúp tiết kiệm hơn 40% thời gian thử nghiệm thuật toán trên phần cứng bay thực tế.

  2. Giảng viên và học viên cao học ngành Kỹ thuật Điều khiển và Tự động hóa: Sử dụng làm tài liệu tham khảo chuẩn mực về phương pháp mô hình hóa Euler - Lagrange, kỹ thuật tối ưu hóa GA đa mục tiêu và quy trình chứng minh ổn định Lyapunov cho hệ thích nghi MRAS.

  3. Kỹ sư thiết kế robot công nghiệp và cơ điện tử: Vận dụng cấu trúc học truyền thẳng LFFC để khử mô men quán tính biến thiên và lực ma sát phi tuyến trên các cánh tay robot có từ 2 đến 6 bậc tự do.

  4. Chuyên gia phát triển hệ thống nhúng thời gian thực: Khai thác quy trình cấu hình giao tiếp giữa MATLAB/Simulink và card dSPACE DS1103 qua phần mềm ControlDesk với chu kỳ lấy mẫu chuẩn xác 1 ms.

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

  1. Tại sao nghiên cứu lựa chọn mô hình Euler - Lagrange thay vì Euler - Newton? Mô hình Euler - Lagrange giữ nguyên bản chất phi tuyến và tương tác xen kênh thực tế của TRMS, khắc phục triệt để các sai số do giả thiết tuyến tính hóa bậc 1 của Euler - Newton. Kết quả thực nghiệm cho thấy mô hình này đạt độ chính xác cao với sai số toàn phương trung bình MSE của góc trệch ngang là 0.0053 và góc chao đứng là 0.0041.

  2. Giải thuật di truyền GA đóng vai trò gì trong việc thiết kế bộ điều khiển PID? Giải thuật GA tiếp nhận 6 tham số khởi tạo từ phương pháp Ziegler - Nichols và tiến hành tìm kiếm tối ưu toàn cục không gian tham số. GA cực tiểu hóa đồng thời 3 hàm mục tiêu IAE, ITAE và MSE, giúp giảm độ quá điều chỉnh từ 24.5% xuống dưới 4.6% mà không cần can thiệp thủ công vào mô hình toán.

  3. Nguyên lý bù tổng của bộ điều khiển LFFC trên cơ sở MRAS là gì? Thay vì phân tách và bù từng kênh đơn lẻ, cấu trúc LFFC-MRAS coi toàn bộ hiện tượng phi tuyến và ghép kênh là một nguồn nhiễu tổng thể. Thuật toán tự học trực tuyến giải phương trình Lyapunov với hệ số tăng tốc $\lambda = 10$, sinh ra tín hiệu truyền thẳng đối nghịch để triệt tiêu hơn 78% tương tác chéo giữa 2 rotor.

  4. Phần cứng dSPACE DS1103 và bộ lọc SVF có vai trò gì trong nghiên cứu? Card xử lý dSPACE DS1103 đóng vai trò là bộ điều khiển thời gian thực với chu kỳ trích mẫu 1 ms (1000 Hz), kết nối đồng bộ giữa giao diện ControlDesk và mô hình TRMS. Bộ lọc biến trạng thái SVF bậc 2 với dải thông gấp 10 lần mô hình tham chiếu giúp trích xuất đạo hàm sai lệch an toàn và lọc sạch nhiễu cảm biến.

  5. Cấu trúc lai đề xuất có ưu thế gì vượt trội so với bộ điều khiển PID thông thường? Bộ điều khiển PID truyền thống chỉ hoạt động tốt quanh điểm làm việc tuyến tính hóa và dễ mất ổn định khi có hiện tượng xen kênh mạnh. Cấu trúc kết hợp PID-GA và LFFC-MRAS giúp rút ngắn 36% thời gian xác lập, dập tắt dao động nhanh chóng và đảm bảo hệ thống ổn định tiệm cận ngay cả khi có nhiễu ngẫu nhiên tác động.

Kết luận

Những đóng góp khoa học và giá trị ứng dụng chính của luận văn bao gồm:

  • Thiết lập thành công mô hình toán học Euler - Lagrange phi tuyến toàn phần cho hệ Twin Rotor MIMO System với sai số mô hình MSE đạt mức 0.0041 trên trục đứng và 0.0053 trên trục ngang.
  • Tối ưu hóa chuẩn xác bộ 6 tham số điều khiển PID kép bằng giải thuật di truyền GA dựa trên việc cực tiểu hóa đồng thời 3 hàm mục tiêu IAE, ITAE và MSE.
  • Thiết kế hoàn chỉnh mạch vòng truyền thẳng học thích nghi LFFC dựa trên MRAS, chứng minh ổn định tiệm cận theo tiêu chuẩn Lyapunov với hệ số tăng tốc $\lambda = 10$.
  • Thử nghiệm thành công trên card xử lý dSPACE DS1103 ở tần số 1000 Hz, chứng minh khả năng triệt tiêu hơn 78% tác động xen kênh và rút ngắn 36% thời gian xác lập.
  • Cung cấp giải pháp điều khiển bù tổng tin cậy, ứng dụng hiệu quả cho các hệ thống MIMO phi tuyến phức tạp trong công nghiệp tự động hóa.

Trong kế hoạch phát triển từ 6 đến 12 tháng tới, đề tài định hướng mở rộng thử nghiệm trên mô hình khí cụ bay 3 bậc tự do (3-DOF) và tích hợp mạng nơ-ron sâu thích nghi. Các nhà nghiên cứu và kỹ sư quan tâm có thể khai thác trực tiếp sơ đồ cấu trúc điều khiển và thuật toán mô phỏng Simulink trong luận văn để ứng dụng vào các dự án phát triển thiết bị bay thông minh.