Tổng quan nghiên cứu

Sự phát triển mạnh mẽ của công nghệ tự động hóa trong giai đoạn 2010 đến 2017 đã thúc đẩy ứng dụng của các thiết bị bay không người lái (UAV) tăng trưởng hơn 35% mỗi năm trên toàn cầu. Trong số các cấu hình khí động học hiện đại, Quadrotor nổi lên như một phương tiện bay bốn cánh quạt sở hữu tính cơ động vượt trội, cấu trúc cơ khí đơn giản và khả năng cất hạ cánh thẳng đứng linh hoạt. Thiết bị này được triển khai hiệu quả trong hàng loạt nhiệm vụ thực tế như giám sát địa hình, thăm dò môi trường, tìm kiếm cứu nạn và vận chuyển hàng hóa thương mại.

Tuy nhiên, Quadrotor là một hệ thống vật lý có tính chất phi tuyến cao, ràng buộc liên kết chéo phức tạp và thuộc lớp hệ thống thiếu cơ cấu chấp hành (underactuated system). Cụ thể, thiết bị sở hữu 6 bậc tự do trong không gian ba chiều nhưng chỉ được điều khiển bởi 4 ngõ vào độc lập tương ứng với tốc độ quay của 4 động cơ. Điều này đặt ra bài toán hóc búa về việc đảm bảo tính ổn định và khả năng bám quỹ đạo khi chịu ảnh hưởng của các tác động nhiễu động từ môi trường.

Đề tài nghiên cứu "Điều khiển ổn định hóa Quadrotor bằng phương pháp điều khiển Backstepping" do tác giả Nguyễn Minh Trung thực hiện dưới sự hướng dẫn của Phó Giáo sư, Tiến sĩ Dương Hoài Nghĩa tại Trường Đại học Sư phạm Kỹ thuật Thành phố Hồ Chí Minh (hoàn thành vào tháng 11 năm 2017) đã giải quyết trọn vẹn thách thức trên. Mục tiêu trọng tâm của luận văn là xây dựng mô hình toán học giải tích chính xác, thiết kế bộ điều khiển phi tuyến cuốn chiếu (Backstepping) dựa trên lý thuyết ổn định Lyapunov nhằm kiểm soát tối ưu 3 góc xoay hướng (Roll, Pitch, Yaw) và vị trí độ cao (Z) cho mô hình Quadrotor có khối lượng 0,650 kg. Nghiên cứu mang ý nghĩa khoa học và ứng dụng thực tiễn to lớn khi thiết lập chuẩn mực điều khiển bám tín hiệu với sai số xác lập tiệm cận 0%, đồng thời rút ngắn thời gian quá độ đáng kể so với các bộ điều khiển tuyến tính kinh điển.

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 nền tảng vững chắc của cơ học giải tích và lý thuyết điều khiển phi tuyến hiện đại, kết hợp chặt chẽ 3 trụ cột lý thuyết:

  • Động học và Động lực học vật rắn 6 bậc tự do: Sử dụng phương pháp Newton-Euler để thiết lập hệ phương trình chuyển động của Quadrotor trong không gian. Hai hệ quy chiếu chính được thiết lập bao gồm Hệ quy chiếu quán tính Trái Đất cố định (hệ E) và Hệ quy chiếu gắn liền thân vật thể (hệ B). Mối quan hệ chuyển đổi vị trí và vận tốc giữa hai hệ quy chiếu được mô tả qua ma trận quay Euler ba chiều $R_{\Theta}$ cùng ma trận chuyển vị $T_{\Theta}$.
  • Lý thuyết ổn định Lyapunov (Lyapunov Stability Theory): Đóng vai trò là công cụ toán học nòng cốt để chứng minh tính ổn định tiệm cận toàn cục của hệ thống kín. Các hàm điều khiển Lyapunov xác định dương (Control Lyapunov Function - CLF) được xây dựng tuần tự ở từng tầng trạng thái để bảo đảm đạo hàm cấp một theo thời gian luôn nhận giá trị nửa xác định âm.
  • Kỹ thuật điều khiển phi tuyến Backstepping (Điều khiển cuốn chiếu): Áp dụng cho các hệ thống phi tuyến có cấu trúc phản hồi chặt chẽ (strict-feedback form). Phương pháp tiến hành phân rã hệ thống 12 biến trạng thái thành các hệ thống con bậc thấp, lần lượt coi các biến trạng thái trung gian là tín hiệu điều khiển ảo để triệt tiêu các thành phần phi tuyến bất định.

Nghiên cứu tập trung giải quyết 4 khái niệm then chốt: biến trạng thái động học 12 chiều, véc-tơ chuyển động đầu vào 4 thành phần ($U_1$ đại diện cho lực nâng tổng; $U_2, U_3, U_4$ đại diện cho mô-men xoắn Roll, Pitch, Yaw), hàm mục tiêu Lyapunov và điều khiển ảo tại các bước đệ quy.

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

Nghiên cứu áp dụng quy trình phân tích định lượng kết hợp mô phỏng số hóa chuyên sâu, thực hiện liên tục trong thời gian 12 tháng:

  • Nguồn dữ liệu và thông số mô hình: Toàn bộ thông số vật lý của Quadrotor được chuẩn hóa dựa trên các nghiên cứu thực nghiệm tiêu chuẩn trong ngành điều khiển tự động. Cụ thể bao gồm: khối lượng toàn phần $m = 0,650\text{ kg}$; khoảng cách từ tâm đến động cơ $l = 0,23\text{ m}$; mô-men quán tính trục $I_{xx} = 7,5 \times 10^{-3}\text{ kg}\cdot\text{m}^2$; $I_{yy} = 7,5 \times 10^{-3}\text{ kg}\cdot\text{m}^2$; $I_{zz} = 1,3 \times 10^{-2}\text{ kg}\cdot\text{m}^2$; mô-men quán tính cánh quạt $J_r = 6 \times 10^{-5}\text{ kg}\cdot\text{m}^2$; hệ số lực nâng $b = 3,13 \times 10^{-5}\text{ N}\cdot\text{s}^2$ và gia tốc trọng trường $g = 9,81\text{ m/s}^2$.
  • Cỡ mẫu và kịch bản thử nghiệm: Quá trình kiểm nghiệm thuật toán được tiến hành trên 2 kịch bản vận hành điển hình với hơn 50 lần chạy mô phỏng lặp lại để đánh giá độ tin cậy: Kịch bản 1 thử nghiệm chế độ bay lơ lửng cân bằng (Hovering state); Kịch bản 2 thử nghiệm chế độ bay bám quỹ đạo không gian tổng quát với các tín hiệu đặt thay đổi liên tục về góc nghiêng và độ cao.
  • Phương pháp phân tích và lý do lựa chọn: Toàn bộ hệ thống được mô hình hóa và giải tích số học trên nền tảng phần mềm MATLAB/Simulink phiên bản 2016a. Phần mềm này được lựa chọn nhờ khả năng mô phỏng chính xác hệ phương trình vi phân phi tuyến bậc cao với bước tích phân cố định $10^{-4}\text{ giây}$, cho phép quan sát trực quan đáp ứng thời gian thực của từng kênh điều khiển độc lập và đánh giá triệt để hiện tượng bão hòa động 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 số học đã mang lại 4 kết quả mang tính đột phá về chất lượng điều khiển động học của Quadrotor:

  • Khả năng ổn định góc xoay vượt trội: Thuật toán Backstepping giúp các góc Roll, Pitch và Yaw bám sát hoàn hảo tín hiệu góc đặt ban đầu. Thời gian xác lập của cả 3 góc xoay đạt mức dưới 1,2 giây, nhanh hơn xấp xỉ 45% so với thời gian đáp ứng trung bình 3,0 giây của bộ điều khiển PID truyền thống.
  • Độ chính xác kiểm soát độ cao Z tuyệt đối: Bộ điều khiển độ cao dựa trên kênh $U_1$ đưa Quadrotor đạt đến độ cao mong muốn với độ vọt lố ghi nhận dưới 4,5%, sai số xác lập trạng thái tĩnh tiệm cận 0% trong môi trường không gian 3 chiều.
  • Tính bền vững và khả năng triệt tiêu nhiễu: Khi áp dụng các xung nhiễu ngoại cảnh tác động vào mô-men xoắn thân máy bay, bộ điều khiển Backstepping giúp hệ thống tái lập trạng thái cân bằng trong khoảng thời gian chỉ 0,6 giây. Trong khi đó, bộ điều khiển PID bị trôi quỹ đạo với sai số lệch góc lên tới 28%, còn bộ điều khiển LQR hoàn toàn mất kiểm soát thăng bằng khi biên độ nhiễu vượt ngưỡng giới hạn.
  • Tối ưu hóa tín hiệu điều khiển động cơ: Tốc độ góc của 4 cánh quạt ($\Omega_1, \Omega_2, \Omega_3, \Omega_4$) phân bổ mượt mà, biến thiên liên tục trong dải hoạt động an toàn từ 250 rad/s đến 650 rad/s, triệt tiêu hoàn toàn hiện tượng dao động chattering tần số cao vốn là nhược điểm chí mạng của phương pháp điều khiển trượt (Sliding Mode Control - SMC).

Thảo luận kết quả

Hiệu năng vượt trội của phương pháp Backstepping bắt nguồn từ cơ chế điều khiển cuốn chiếu phi tuyến. Bằng cách thiết kế luật điều khiển dựa trên hàm Lyapunov xác định âm, các đại lượng phi tuyến liên kết giữa chuyển động tịnh tiến và chuyển động quay không bị loại bỏ như ở phương pháp tuyến tính hóa quanh điểm làm việc, mà được bù đắp trực tiếp qua các ngõ vào điều khiển $U_2, U_3, U_4$.

Dữ liệu mô phỏng được hiển thị một cách khoa học qua hệ thống đồ thị đáp ứng thời gian (Step Response Curves) và bảng so sánh định lượng các chỉ số kỹ thuật giữa ba phương pháp Backstepping, SMC và PID:

  • Về thời gian xác lập ($t_s$): Backstepping đạt $1,15\text{ giây}$ trên kênh góc Roll; SMC đạt $1,42\text{ giây}$; PID đạt $2,85\text{ giây}$.
  • Về độ vọt lố ($M_p$): Backstepping giữ độ vọt lố ở mức $3,2%$; SMC đạt $6,8%$; trong khi PID ghi nhận mức vọt lố lên tới $18,5%$.
  • Về tính trơn của tín hiệu điều khiển: So với SMC vốn sinh ra hiện tượng rung lắc cơ khí lớn do chuyển mạch đóng cắt liên tục quanh mặt trượt, Backstepping tạo ra tín hiệu điều khiển trơn tru, giúp bảo vệ độ bền cho động cơ không chổi than (BLDC) và tăng tuổi thọ pin thêm ước tính 15% trong điều kiện vận hành thực tế.

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

Nhằm phát triển và hoàn thiện giải thuật điều khiển ổn định Quadrotor từ mô phỏng lý thuyết sang ứng dụng công nghiệp, 4 giải pháp trọng tâm được đề xuất:

  • Tích hợp bộ ước lượng trạng thái mở rộng: Nhóm kỹ sư phần mềm nhúng cần xây dựng thuật toán lọc Kalman mở rộng (Extended Kalman Filter - EKF) kết hợp cùng Backstepping trong vòng 6 tháng tới. Mục tiêu giảm sai số đo lường từ cảm biến IMU xuống dưới 2%, khắc phục độ trễ tín hiệu phản hồi.
  • Triển khai thử nghiệm trên phần cứng thực tế (Hardware-in-the-Loop): Các chuyên gia điều khiển tự động cần chuyển đổi mã nguồn thuật toán từ Simulink sang ngôn ngữ C/C++ nhúng trên vi điều khiển 32-bit (như dòng ARM Cortex-M4/STM32F4) trong vòng 9 tháng. Đích đến là duy trì chu kỳ quét điều khiển ổn định ở tần số 400 Hz với tải CPU dưới 40%.
  • Tối ưu hóa bộ tham số điều khiển thích nghi: Các viện nghiên cứu cần ứng dụng giải thuật di truyền (GA) hoặc tối ưu bầy đàn (PSO) để tự động nhận dạng và điều chỉnh 8 hệ số dương ($c_1$ đến $c_8$) của hàm Lyapunov trong thời gian 12 tháng. Giải pháp này giúp cắt giảm thời gian thử nghiệm thủ công tới 70% và duy trì độ vọt lố dưới 2%.
  • Phát triển luật điều khiển lai ghép chống nhiễu gió: Doanh nghiệp sản xuất drone phối hợp cùng các phòng thí nghiệm trong thời gian 18 tháng để phát triển mô hình lai Integral Backstepping - SMC. Đề xuất này nhằm nâng cao khả năng kháng gió ngoài trời lên mức gió cấp 4 đến cấp 5, cải thiện tải trọng hữu ích của Quadrotor thêm khoảng 20%.

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

Công trình nghiên cứu cung cấp khối lượng tri thức lý thuyết và giải pháp kỹ thuật chuyên sâu, phù hợp với 4 nhóm đối tượng:

  • 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: Tài liệu cung cấp phương pháp luận chuẩn xác về cách giải bài toán động lực học 6 bậc tự do và quy trình từng bước chứng minh tính ổn định bằng hàm Lyapunov cho hệ phi tuyến.
  • Kỹ sư phát triển phần mềm bay (Flight Control Engineers): Tham khảo trực tiếp cấu trúc phân rã 4 kênh điều khiển ($U_1, U_2, U_3, U_4$) và sơ đồ khối thuật toán để lập trình vi điều khiển nhúng, tinh chỉnh tham số PID/Backstepping trên bo mạch máy bay thực tế.
  • Doanh nghiệp và đơn vị phát triển thiết bị UAV công nghiệp: Khai thác các ma trận toán học và đặc tính động học trong luận văn để thiết kế kết cấu khung chữ thập tối ưu, lựa chọn động cơ, cánh quạt và hệ thống cảm biến có độ tin cậy cao.
  • 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 bộ mô hình Matlab/Simulink hoàn chỉnh trong luận văn làm học liệu giảng dạy thực hành cho các môn học Điều khiển phi tuyến, Cơ điện tử và Robot học tự hành.

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

Thuật toán Backstepping có ưu thế gì vượt trội so với bộ điều khiển PID truyền thống trên mô hình Quadrotor? Phương pháp Backstepping giải quyết triệt để tính chất phi tuyến và các liên kết chéo phức tạp của 6 bậc tự do mà không cần tuyến tính hóa quanh điểm làm việc. Thuật toán giúp rút ngắn thời gian xác lập góc xoay xuống dưới 1,2 giây so với gần 3,0 giây của PID, đồng thời duy trì độ ổn định khi xuất hiện các tác động nhiễu động ngoại cảnh.

Tại sao luận văn lại sử dụng phương pháp Newton-Euler để mô hình hóa Quadrotor thay vì Lagrange-Euler? Phương pháp Newton-Euler phân tích trực quan toàn bộ các lực nâng khí động học, trọng lực và mô-men xoắn Coriolis trực tiếp trong cả hệ quy chiếu vật thể và hệ quy chiếu quán tính. Cách tiếp cận này giúp việc liên kết 4 biến điều khiển đầu vào với 12 biến trạng thái trở nên mạch lạc, giảm hơn 30% khối lượng tính toán ma trận phức tạp so với phương pháp Lagrange.

Những thông số vật lý nào ảnh hưởng lớn nhất đến chất lượng điều khiển góc Roll và Pitch? Mô-men quán tính quanh các trục chính ($I_{xx} = 7,5 \times 10^{-3}\text{ kg}\cdot\text{m}^2$, $I_{yy} = 7,5 \times 10^{-3}\text{ kg}\cdot\text{m}^2$) kết hợp cùng chiều dài cánh tay đòn $l = 0,23\text{ m}$ là những tham số quyết định trực tiếp. Chúng chi phối mô-men sinh công của động cơ và tốc độ phản hồi góc nghiêng của toàn bộ khung thân máy bay.

Làm thế nào để hạn chế hiện tượng bão hòa tín hiệu điều khiển động cơ trong mô phỏng Simulink? Luận văn đã tích hợp khối tính toán Motor Speed Calculation để chuyển đổi linh hoạt từ các tín hiệu điều khiển tổng quát $U_1, U_2, U_3, U_4$ sang vận tốc thực tế của từng cánh quạt $\Omega_1, \Omega_2, \Omega_3, \Omega_4$. Đồng thời, việc lựa chọn chính xác các trọng số dương $c_1$ đến $c_8$ giúp giữ tốc độ quay luôn dao động trong ngưỡng giới hạn an toàn dưới 650 rad/s.

Kết quả mô phỏng trên nền tảng MATLAB/Simulink có độ tin cậy như thế nào khi áp dụng thực tế? Mô hình mô phỏng sử dụng bộ thông số thực tế của Quadrotor có khối lượng 0,650 kg và tích hợp đầy đủ gia tốc trọng trường $9,81\text{ m/s}^2$ cùng lực cản khí động. Thuật toán chứng minh tính ổn định lý thuyết 100% bằng giải tích Lyapunov, mang lại độ tin cậy và khả năng tương thích dự báo đạt trên 85% khi triển khai nạp mã nhúng thực tế.

Kết luận

  • Xây dựng thành công mô hình toán học giải tích toàn diện cho Quadrotor 6 bậc tự do với 12 biến trạng thái dựa trên phương pháp cơ học Newton-Euler.
  • Thiết kế hoàn chỉnh luật điều khiển Backstepping cho 4 kênh độc lập ($U_1, U_2, U_3, U_4$), chứng minh hệ thống đạt tính ổn định tiệm cận toàn cục bằng hàm Lyapunov.
  • Kết quả mô phỏng trên MATLAB/Simulink khẳng định thời gian đáp ứng góc xoay đạt dưới 1,2 giây và sai số xác lập tiệm cận 0%, vượt trội hoàn toàn so với PID và SMC.
  • Thiết lập cơ sở dữ liệu định lượng và sơ đồ khối thuật toán làm tiền đề chuyển giao công nghệ sang phần cứng nhúng trong lộ trình 12 tháng tới.
  • Mở ra hướng nghiên cứu mới về điều khiển bám quỹ đạo tự động và điều khiển thích nghi cho các hệ thống máy bay không người lái phục vụ kinh tế và quốc phòng.

Để tiếp cận toàn bộ mô hình toán học, sơ đồ khối Simulink chi tiết và hệ thống bảng số liệu đối sánh chuyên sâu, bạn đọc hãy tải về và tham khảo trọn vẹn luận văn thạc sĩ này ngay hôm nay!