Tổng quan nghiên cứu

Sự phát triển mạnh mẽ của các phương tiện bay không người lái bốn cánh quạt (quadrotor) đang tạo ra bước đột phá lớn trong các nhiệm vụ thám hiểm, giám sát không gian và vận chuyển cứu hộ. Theo ước tính từ các báo cáo kỹ thuật hàng không, phân khúc thiết bị bay cỡ nhỏ chiếm hơn 65% tổng số các nghiên cứu ứng dụng robot bay hiện đại. Tuy nhiên, quadrotor là hệ thống phi tuyến điển hình với 6 bậc tự do nhưng chỉ được điều khiển bởi 4 tín hiệu lực nâng độc lập từ các động cơ không chổi than, tạo nên cấu trúc thiếu cơ cấu chấp hành (underactuated system) phức tạp.

Vấn đề cốt lõi đặt ra trong nghiên cứu là duy trì trạng thái cân bằng ổn định, kiểm soát triệt để các góc nghiêng động học Roll (lăn), Pitch (chúi), Yaw (quay hướng) và khả năng bám quỹ đạo chính xác khi chịu tác động của nhiễu môi trường. Mục tiêu của luận văn là xây dựng mô hình toán học giải tích chi tiết của hệ thống quadrotor, thiết kế và so sánh hiệu năng giữa thuật toán điều khiển cổ điển vi tích phân tỷ lệ - vi phân (PD) với thuật toán điều khiển phi tuyến hiện đại cuốn chiếu (Backstepping), đồng thời hiện thực hóa giải thuật trên nền tảng vi điều khiển nhúng 32-bit hiệu năng cao.

Nghiên cứu được triển khai thực nghiệm tại Trường Đại học Bách Khoa – ĐHQG TP.HCM trong giai đoạn 2011 đến 2013. Luận văn tiến hành khảo sát toàn diện qua 20 giây mô phỏng trên từng kịch bản tải trọng từ 0,5 kg đến 1,0 kg và góc lệch ban đầu lên tới 1,57 rad. Ý nghĩa khoa học và thực tiễn của công trình thể hiện qua việc rút ngắn thời gian đáp ứng cân bằng xuống xấp xỉ 1,0 giây, giảm độ vọt lố dưới 8% và duy trì ổn định góc nghiêng với sai lệch tĩnh nhỏ hơn 0,05 rad, cung cấp giải pháp khả thi cho việc chế tạo thiết bị bay tự hành quy mô 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 tích hợp hai nền tảng lý thuyết điều khiển tự động kinh điển và hiện đại để làm chủ động học quadrotor:

  • Lý thuyết động lực học vật rắn phi tuyến: Mô hình hóa 6 bậc tự do dựa trên phương trình Newton-Euler và biến phân năng lượng Euler-Lagrange. Hệ thống được biểu diễn qua 12 biến trạng thái bao gồm tọa độ vị trí không gian (x, y, z), vận tốc dài, ba góc Euler (Roll $\phi$, Pitch $\theta$, Yaw $\psi$) và vận tốc góc tương ứng.
  • Lý thuyết điều khiển cuốn chiếu (Backstepping Control): Phương pháp điều khiển đệ quy từng bước dựa trên hàm mục tiêu xác định dương Lyapunov. Bằng cách phân rã hệ thống bậc cao thành các chuỗi hệ thống con bậc thấp, thuật toán tổng hợp luật điều khiển phản hồi phi tuyến đảm bảo tính ổn định tiệm cận toàn cục cho toàn bộ hệ thống.
  • Lý thuyết điều khiển vi phân tỷ lệ (PD): Cấu trúc điều khiển vòng kín rút gọn từ bộ PID tiêu chuẩn nhằm triệt tiêu sai lệch tức thời thông qua độ lợi khuếch đại sai số và vận tốc thay đổi sai số.

Các khái niệm then chốt trong nghiên cứu bao gồm: hiệu ứng hồi chuyển thân và cánh quạt ($J_r \Omega_r$), ma trận chuyển đổi tọa độ từ hệ quy chiếu gắn trên thân bay sang hệ quy chiếu cố định mặt đất, và cơ chế bù trừ ngẫu lực xoắn thông qua hai cặp cánh quạt quay nghịch chiều.

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

Nghiên cứu sử dụng phương pháp mô phỏng số kết hợp thực nghiệm nhúng với các thông số cụ thể:

  • Cỡ mẫu và kịch bản thử nghiệm: Thực hiện trên 25 kịch bản mô phỏng độc lập với chu kỳ lấy mẫu 0,01 giây trong thời gian khảo sát 20 giây mỗi lượt. Quá trình quét thông số bao gồm hệ số tỷ lệ $K_p$ biến thiên từ 0,1 đến 1,2; hệ số vi phân $K_d$ biến thiên từ 0,2 đến 0,8; tải trọng thử nghiệm gồm 0,5 kg, 0,8 kg và 1,0 kg; góc lệch ban đầu khảo sát tại các mức 0,2 rad, 0,4 rad, 0,5 rad và $\pi/2$ rad ($1,57$ rad).
  • Phương pháp chọn mẫu: Chọn mẫu thực nghiệm ngẫu nhiên có kiểm soát thông qua việc đưa vào tín hiệu nhiễu trắng nhân tạo tác động trực tiếp lên kênh đo cảm biến góc với biên độ 0,05 rad.
  • Lý do lựa chọn phương pháp phân tích: Sử dụng môi trường Matlab/Simulink kết hợp khối hàm nhúng S-Function giúp phân tách rõ ràng các kênh điều khiển Roll, Pitch, Yaw độc lập. Phương pháp này cho phép kiểm chứng sự hội tụ của hàm Lyapunov trước khi biên dịch mã nguồn thực thi xuống vi điều khiển ARM Cortex-M4 STM32F407VG hoạt động ở xung nhịp 168 MHz.

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 các kết quả nổi bật về đáp ứng động học:

  • Tối ưu hóa thông số bộ điều khiển PD: Bộ thông số tốt nhất đạt được tại $K_p = 1,0$ và $K_d = 0,4$, cho thời gian xác lập cực nhanh khoảng 1,0 giây. Khi giảm $K_p$ xuống 0,1 hoặc 0,4, thời gian đáp ứng bị kéo dài hơn 250%, làm thiết bị mất thăng bằng nghiêm trọng. Khi tăng $K_p$ lên 1,2, hệ thống xuất hiện hiện tượng vọt lố vượt mức 15% đi kèm dao động điều hòa kéo dài.
  • Giới hạn chịu tải của thuật toán: Với tải trọng 0,5 kg và 0,8 kg, bộ điều khiển PD duy trì tốt khả năng bám điểm đặt với độ vọt lố góc Roll ban đầu dưới 8%. Tuy nhiên, khi tăng khối lượng tải lên 1,0 kg (tương đương mức tăng 100% so với cấu hình chuẩn), hệ thống PD hoàn toàn mất kiểm soát và phân kỳ trạng thái.
  • Ngưỡng lệch góc ban đầu: Bộ điều khiển PD chỉ có khả năng ổn định góc quay hướng Yaw khi góc lệch ban đầu nhỏ hơn 0,5 rad. Cụ thể, tại giá trị 0,2 rad và 0,4 rad, hệ thống bám điểm cân bằng trong 1,2 giây, nhưng tại ngưỡng 0,5 rad, thiết bị quay tròn liên tục mất định hướng. Tương tự, đối với góc Roll và Pitch, góc lệch đạt ngưỡng $\pi/2$ rad khiến hệ dao động mạnh và mất thời gian hồi phục rất lâu.
  • Khả năng kháng nhiễu cảm biến: Khi đưa tín hiệu nhiễu trắng vào cảm biến góc ở mức tải 0,8 kg, thuật toán vẫn giữ vững khả năng bám tín hiệu đặt, sai số dao động quanh vùng cân bằng không vượt quá 0,05 rad.

Thảo luận kết quả

Nguyên nhân chính khiến bộ điều khiển PD mất ổn định ở tải trọng 1,0 kg và góc lệch lớn xuất phát từ việc thuật toán PD tuyến tính hóa hệ thống quanh điểm làm việc cân bằng tĩnh. Khi khối lượng hoặc góc nghiêng tăng cao, các thành phần phi tuyến như mô-men quán tính chéo, lực hồi chuyển cánh quạt và ma sát khí động học tăng theo hàm phi tuyến bậc hai, vượt ra ngoài vùng bù sai lệch của các hệ số cố định.

Để minh họa trực quan, dữ liệu đáp ứng động học có thể được trình bày thông qua biểu đồ chuỗi thời gian so sánh góc nghiêng từ 0 đến 20 giây và bảng ma trận đối chuẩn chỉ số: thời gian xác lập ($t_s$), độ vọt lố ($M_p$) và sai số xác lập ($e_{ss}$). Trong khi PD bị giới hạn nghiêm ngặt về biên độ điều khiển, thuật toán Backstepping nhờ vào cơ chế cuốn chiếu tích phân và triệt tiêu thành phần phi tuyến bằng đạo hàm Lyapunov đã kiểm soát xuất sắc góc lệch ngay cả khi chịu tải lớn và góc ban đầu vượt quá 1,57 rad. Kết quả này hoàn toàn tương thích với các nghiên cứu robot bay quốc tế, khẳng định tính tất yếu của việc ứng dụng điều khiển phi tuyến trên quadrotor.

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

Dựa trên kết quả nghiên cứu, 4 khuyến nghị then chốt được đề xuất nhằm nâng cấp chất lượng điều khiển quadrotor:

  • Nâng cấp thuật toán sang cấu trúc lai phi tuyến: Triển khai bộ điều khiển kết hợp giữa Backstepping và chế độ trượt (Sliding Mode) hoặc kỹ thuật thích nghi (Adaptive Backstepping). Mục tiêu giảm thời gian xác lập xuống dưới 0,8 giây, triệt tiêu hiện tượng rung chấn góc dưới 2% và mở rộng dải tải trọng hữu ích lên mức 2,0 kg đến 2,5 kg. Thời gian thực hiện dự kiến 6 tháng, do nhóm nghiên cứu thuật toán điều khiển đảm nhiệm.
  • Tích hợp bộ lọc Kalman mở rộng (EKF) đa tầng: Xây dựng giải thuật ước lượng trạng thái kết hợp dữ liệu từ con quay hồi chuyển (gyroscope), cảm biến gia tốc (accelerometer) và từ kế (magnetometer). Mục tiêu lọc bỏ 95% thành phần nhiễu trắng môi trường và khắc phục hiện tượng trôi điểm không trong vòng 3 tháng, do kỹ sư xử lý tín hiệu thực hiện.
  • Tối ưu hóa phần cứng và tần số điều khiển nhúng: Nâng cấp mạch điều khiển trung tâm trên chip STM32F407VG, thiết lập xung nhịp xử lý vòng lặp điều khiển đạt tần số 400 Hz và giảm độ trễ phản hồi của mạch điều tốc ESC xuống dưới 5 mili giây. Kế hoạch triển khai trong 4 tháng, do nhóm phát triển phần cứng phụ trách.
  • Mở rộng kịch bản thử nghiệm bay ngoài trời: Tiến hành thử nghiệm bay bám quỹ đạo 3D thực tế dưới tác động của gió giật cấp 3 đến cấp 4 (vận tốc gió 3,5 m/s đến 5,5 m/s). Mục tiêu duy trì sai số định vị không gian điểm treo (hovering) nhỏ hơn 0,2 mét. Thời gian kiểm định 6 tháng, thực hiện bởi nhóm thử nghiệm bay hiện trường.

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

Nội dung luận văn mang lại giá trị học thuật và ứng dụng thực tiễn cao cho 4 nhóm đối tượng:

  • Học viên cao học và sinh viên chuyên ngành Tự động hóa, Kỹ thuật Điều khiển và Robotics: Nắm bắt quy trình mô hình hóa giải tích 6 bậc tự do và phương pháp thiết kế bộ điều khiển Backstepping toàn diện từ lý thuyết đến mô phỏng Simulink để ứng dụng trực tiếp vào các đồ án tốt nghiệp và đề tài nghiên cứu chuyên sâu.
  • Kỹ sư phát triển hệ thống điều khiển bay (Flight Control Engineers): Tham khảo kiến trúc mã nguồn nhúng trên dòng vi điều khiển ARM Cortex-M4 STM32F407VG, phương pháp lập trình điều chế xung PWM 4 kênh động cơ và kỹ thuật đọc cảm biến thời gian thực nhằm chế tạo các bo mạch lái tự động thương mại.
  • Giảng viên và nhà nghiên cứu tại các viện, trường đại học: Sử dụng các số liệu mô phỏng, phương trình động học Newton-Euler và kết quả khảo sát tham số làm tài liệu giảng dạy cho các môn học Điều khiển phi tuyến, Robot tự hành và Kỹ thuật điều khiển nhúng.
  • Doanh nghiệp nghiên cứu và phát triển thiết bị bay không người lái (UAV R&D Centers): Đánh giá chính xác ưu nhược điểm giữa cấu hình điều khiển PD đơn giản và thuật toán Backstepping cao cấp để tối ưu hóa chi phí sản xuất, lựa chọn cấu hình phần cứng phù hợp cho từng dòng drone chuyên dụng.

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

Tại sao luận văn lại lựa chọn cấu trúc bộ điều khiển PD thay vì PID đầy đủ cho mô hình quadrotor? Trong các hệ thống bay cơ động nhanh như quadrotor, khâu tích phân (I) dễ gây ra hiện tượng bão hòa tích phân (integrator windup) làm trễ đáp ứng và tăng nguy cơ mất ổn định khi đổi hướng đột ngột. Do đó, việc sử dụng bộ điều khiển PD rút gọn giúp tối ưu hóa tốc độ phản hồi, giảm thiểu gánh nặng tính toán trên bộ xử lý nhúng mà vẫn duy trì được độ chính xác vị trí góc.

Yếu tố nào dẫn đến hiện tượng mất ổn định của bộ điều khiển PD khi tải trọng đạt mức 1,0 kg? Khi tải trọng đạt 1,0 kg, quán tính thân bay thay đổi vượt quá giới hạn bù trừ của các hệ số khuếch đại cố định $K_p = 1,0$ và $K_d = 0,4$. Động cơ không đủ dải điều khiển tuyến tính để triệt tiêu mô-men quán tính phát sinh, dẫn đến việc vòng lặp điều khiển bị trễ pha và hệ thống phân kỳ nhanh chóng chỉ sau khoảng 1,5 giây.

Thuật toán điều khiển Backstepping thể hiện ưu thế vượt trội nào so với giải thuật PD cổ điển? Backstepping giải quyết triệt để các thành phần liên kết phi tuyến chéo giữa 3 trục Roll, Pitch, Yaw bằng cách sử dụng các hàm biến ảo Lyapunov qua từng bước tích phân. Nhờ đó, thuật toán đảm bảo tính ổn định tiệm cận toàn cục, cho phép quadrotor thích nghi với góc lệch lớn trên 1,57 rad và duy trì trạng thái bay ổn định trong điều kiện biến động tải trọng.

Nền tảng vi điều khiển STM32F407VG đáp ứng các yêu cầu điều khiển thời gian thực như thế nào? Vi điều khiển STM32F407VG sở hữu lõi 32-bit ARM Cortex-M4 với bộ xử lý dấu phẩy động FPU phần cứng chạy ở tốc độ 168 MHz. Nền tảng này cho phép thực thi các phép tính ma trận lượng giác phức tạp của thuật toán Backstepping với thời gian trễ dưới 2,5 mili giây, đồng thời xuất tín hiệu PWM tần số cao điều khiển 4 mạch ESC mượt mà.

Làm thế nào để hệ thống duy trì được độ chính xác khi cảm biến đo lường bị tác động bởi nhiễu trắng? Hệ thống tận dụng khâu vi phân có lọc kết hợp thuật toán thích nghi biên độ giúp cô lập thành phần nhiễu tần số cao. Trong các mô phỏng thử nghiệm với biên độ nhiễu 0,05 rad, thuật toán điều khiển vẫn giữ vững quỹ đạo cân bằng, đưa sai số góc trung bình của cả 3 trục về mức dưới 0,02 rad sau khi lọc tín hiệu.

Kết luận

  • Luận văn đã thiết lập thành công mô hình toán học giải tích 6 bậc tự do hoàn chỉnh của máy bay bốn cánh quạt quadrotor dựa trên phương trình Newton-Euler và Lagrange.
  • Khảo sát và tối ưu hóa thành công bộ điều khiển PD tại bộ thông số $K_p = 1,0$ và $K_d = 0,4$, giúp hệ thống đạt thời gian cân bằng lý tưởng xấp xỉ 1,0 giây.
  • Chứng minh rõ nét giới hạn vật lý của bộ điều khiển PD khi khối lượng tải đạt 1,0 kg hoặc góc nghiêng vượt quá 0,5 rad, khẳng định tính tất yếu của việc triển khai thuật toán phi tuyến Backstepping.
  • Hiện thực hóa cấu trúc thuật toán điều khiển nhúng trên vi điều khiển ARM Cortex-M4 STM32F407VG với chu kỳ lấy mẫu cao và độ trễ thấp.
  • Đặt nền móng vững chắc cho lộ trình 12 tháng tiếp theo nhằm tích hợp bộ lọc Kalman mở rộng và thử nghiệm bay tự hành ngoài trời với độ chính xác bám quỹ đạo 100%.

Hãy áp dụng ngay các kết quả nghiên cứu và phương pháp thiết kế bộ điều khiển trong luận văn để nâng cao độ tin cậy và khả năng ổn định cho các dự án phát triển thiết bị bay không người lái của bạn.