Tổng quan nghiên cứu

Lý thuyết điều khiển trượt (Sliding Mode Control - SMC) được khởi xướng bởi Emelyanov vào năm 1970 và phát triển hoàn thiện qua các công trình kinh điển của Itkis năm 1976 cùng Utkin năm 1977. Đây là phương pháp điều khiển phi tuyến nổi bật nhờ khả năng thích nghi cao, bảo đảm tính bền vững trước các biến động tham số mô hình và nhiễu loạn ngoại lai. Dẫu vậy, ứng dụng thực tế của SMC gặp trở ngại lớn ở việc xác định tham số mặt trượt lamda ($\lambda$). Tham số $\lambda$ là một hằng số quyết định trực tiếp đến băng thông và tốc độ hội tụ của hệ thống kín. Thông thường, việc xác định $\lambda$ hoàn toàn dựa vào phương pháp thử sai thủ công, đòi hỏi kỹ sư phải thực nghiệm từ 15 đến 20 lần mô phỏng, tiêu tốn từ 4 đến 8 giờ làm việc cho mỗi mô hình đối tượng mà vẫn khó đạt được giá trị tối ưu.

Nhằm giải quyết triệt để nút thắt này, luận văn thạc sĩ chuyên ngành Tự động hóa (mã số 605206) của tác giả Hứa Minh Tú, dưới sự hướng dẫn khoa học của Tiến sĩ Trần Thanh Hùng tại Trường Đại học Bách Khoa - Đại học Quốc gia TP.HCM, đã được triển khai từ ngày 14/02/2011 và bảo vệ thành công vào tháng 01/2012. Mục tiêu trọng tâm của đề tài là xây dựng giải thuật tự động điều chỉnh thông số $\lambda$ trên nền tảng phần mềm MATLAB/Simulink thông qua khối S-Function, cho phép hệ thống tự động xác định giá trị $\lambda$ tối ưu tương ứng với độ dốc đáp ứng mong muốn của từng đối tượng điều khiển cụ thể. Nghiên cứu thực hiện kiểm chứng mô phỏng trên các mô hình động học bậc hai và hệ điều khiển vị trí động cơ kết hợp cấu trúc SMC-PID. Kết quả nghiên cứu giúp cắt giảm hơn 90% thời gian cài đặt tham số ban đầu, triệt tiêu hoàn toàn độ vọt lố và duy trì sai số bám quỹ đạo tiệm cận 0 trong thời gian đáp ứng cực ngắn.

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

Khung lý thuyết áp dụng

Luận văn xây dựng trên nền tảng lý thuyết điều khiển trượt hiện đại cho hệ thống phi tuyến biểu diễn dưới dạng phương trình trạng thái tổng quát $x^{(n)} = f(x) + b(x)u$, trong đó $x$ là biến trạng thái ngõ ra và $u$ là tín hiệu điều khiển. Mặt trượt biến thiên theo thời gian $s(x, t) = 0$ được xác định theo phương trình vi phân $s(x, t) = (\frac{d}{dt} + \lambda)^{n-1}\tilde{x}$, với $\tilde{x} = x - x_d$ là vector sai lệch chuyển động và $\lambda$ là hằng số dương đại diện cho hệ số góc mặt trượt. Theo tiêu chuẩn ổn định Lyapunov, điều kiện trượt bền vững được bảo đảm qua bất đẳng thức $s\dot{s} \le -\eta|s|$ (với $\eta > 0$), bảo đảm quỹ đạo trạng thái chạm tới mặt trượt trong khoảng thời gian $t_{reach} \le |s(0)|/\eta$ và hội tụ tiệm cận về điểm cân bằng theo quy luật hàm mũ với hằng số thời gian $(n-1)/\lambda$.

Bên cạnh đó, tác giả tích hợp mô hình điều khiển phối hợp giữa chế độ trượt và bộ điều khiển PID (SMC-PID). Cấu trúc nối tiếp này giúp mô hình hóa hệ thống với bộ điều khiển PID thành hệ bậc hai tương đương, từ đó tận dụng đặc tính kháng nhiễu của SMC để giảm bậc động học và loại bỏ hiện tượng vọt lố. Ba khái niệm cốt lõi được định hình gồm: mặt trượt bất biến, tín hiệu điều khiển tương đương ($u_{eq}$) và độ dốc đáp ứng động học đo bằng đơn vị phần trăm trên giây (%/s).

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

Phương pháp nghiên cứu kết hợp mô hình hóa giải thuật toán học và mô phỏng số học chuyên sâu trên phần mềm MATLAB/Simulink R2010a. Cỡ mẫu thử nghiệm bao gồm 3 nhóm đối tượng động học điển hình: mô hình hàm truyền bậc hai chuẩn $1/(s^2+2s+1)$, hệ thống điều khiển vị trí motor phối hợp SMC-PID, và hàm truyền bậc hai có hệ số cản $1/(s^2+3s+2)$. Việc lựa chọn 3 đối tượng này theo phương pháp chọn mẫu có chủ đích dựa trên thực tế các hệ thống bậc hai chiếm hơn 80% cấu trúc truyền động cơ điện tử công nghiệp, đại diện đầy đủ cho các hành vi dao động, trễ pha và phi tuyến.

Giải thuật tự chỉnh định được lập trình bằng ngôn ngữ M-code trong khối S-Function adjust_lamda_parameter.m kết hợp với khối trễ Transport Delay có giá trị đệm 1024 mẫu, sử dụng thuật toán giải vi phân ode45 với bước tích phân thích nghi. Lý do lựa chọn phương pháp này là khả năng phản hồi trực tuyến: thuật toán liên tục so sánh sai số trễ trước đó ($e_{delay}$) với sai số hiện thời ($e$) để tự động giảm dần $\lambda$ từ giá trị khởi tạo 100000 về ngưỡng tối ưu phù hợp với độ dốc đặt trước. Toàn bộ quy trình thiết kế, mô phỏng và kiểm chứng được hoàn thành trong khung thời gian 10 tháng.

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

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

Thứ nhất, nghiên cứu chứng minh quy luật tác động phi tuyến của tham số $\lambda$ lên tốc độ đáp ứng hệ thống. Thử nghiệm trên đối tượng $1/(s^2+2s+1)$ cho thấy khi giảm $\lambda$ từ 3000 xuống 300, 100 rồi 30, thời gian quá độ giảm rõ rệt hơn 60% và hệ thống đạt đáp ứng nhanh vượt trội. Tuy nhiên, khi tiếp tục giảm $\lambda$ xuống mức 3, tốc độ đáp ứng của hệ thống bị chậm lại. Điều này khẳng định tồn tại một giá trị $\lambda$ tối ưu cục bộ duy nhất cho từng cơ cấu truyền động.

Thứ hai, quy trình 3 bước ước lượng $\lambda$ tối ưu thông qua độ dốc đáp ứng được xây dựng và xác thực thành công trên hệ điều khiển vị trí motor SMC-PID. Khi gán độ dốc kích thích ban đầu ở mức 10000 %/s, độ dốc cực đại thực tế đo được đạt 78 %/s. Bằng cách áp dụng mức kiểm tra 20% tương đương độ dốc 15.6 %/s, đáp ứng hệ thống đạt 19 %/s, xác định tỷ lệ thích ứng đạt khoảng 0.82 (15.6/19). Từ đó, độ dốc giới hạn của hệ thống được ước lượng chuẩn xác ở mức 64 %/s ($78 \times 0.82$), giúp thuật toán tự động dò tìm và khóa giá trị $\lambda = 360$ chỉ trong 0.05 giây.

Thứ ba, tính thích nghi của giải thuật được tái kiểm chứng chính xác trên đối tượng $1/(s^2+3s+2)$. Với độ dốc kích thích 10000 %/s cho độ dốc đỉnh 7.5 %/s và mẫu thử 20% đạt 1.87 %/s (tỷ lệ 0.802), độ dốc tối ưu ước lượng là 6.7 %/s, tương ứng giá trị $\lambda = 65$. Khi thử nghiệm đặt các độ dốc lớn hơn như 10 %/s, 15 %/s và 20 %/s, giá trị $\lambda$ tự động bão hòa ở mức $58 \times 10^3$, bảo đảm hệ thống luôn duy trì đáp ứng tiệm cận giới hạn vật lý mà không gây mất ổn định vòng kín.

Thảo luận kết quả

Cơ chế tự điều chỉnh giải quyết triệt để mâu thuẫn giữa tốc độ quá độ và tính ổn định. Khi sai số trước $e_{delay}$ lớn hơn sai số sau $e$, thuật toán nhận diện hệ thống đang hội tụ đúng hướng và tăng dần độ dốc cho đến khi bằng giá trị đặt trước thì cố định $\lambda$. Dữ liệu nghiên cứu được biểu diễn trực quan qua các bảng thông số quá độ và đồ thị phản hồi bước nhảy (Step Response vs Time), đồ thị quỹ đạo pha sai số ($e - \dot{e}$) cùng biểu đồ biến thiên tham số $\lambda$ theo thời gian thực.

So với các phương pháp chỉnh định bằng giải thuật di truyền (GA) hay mạng nơ-ron Hopfield vốn đòi hỏi hàng trăm vòng lặp tính toán phức tạp, giải thuật thích nghi theo độ dốc trong luận văn có cấu trúc đại số tinh giản, thực thi tức thời trong từng chu kỳ lấy mẫu. Cấu trúc SMC-PID kết hợp tự chỉnh $\lambda$ đã triệt tiêu hoàn toàn độ vọt lố (overshoot bằng 0%), giảm thiểu hiện tượng rung giật chattering quanh mặt trượt và khắc phục nhược điểm phụ thuộc vào kinh nghiệm chuyên gia khi cài đặt hệ số ban đầu.

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

Dựa trên kết quả đạt được, 4 khuyến nghị và giải pháp hành động cụ thể được đề xuất nhằm ứng dụng giải thuật vào thực tiễn:

  1. Triển khai lập trình nhúng thuật toán S-Function lên các dòng vi xử lý số DSP TMS320F28335 hoặc ARM Cortex-M4 trong vòng 6 tháng tới, nhằm thương mại hóa bộ điều khiển tự chỉnh cho hệ thống servo công nghiệp, giúp giảm 95% thời gian căn chỉnh thông số tại nhà máy. Chủ thể thực hiện: Doanh nghiệp chế tạo thiết bị tự động hóa và kỹ sư R&D truyền động điện.

  2. Mở rộng giải thuật cho các hệ thống phi tuyến nhiều ngõ vào - nhiều ngõ ra (MIMO) và hệ thống có bậc tự do $n \ge 3$ trong lộ trình 12 tháng, bảo đảm sai số bám quỹ đạo dưới 0.05% trong điều kiện tải động biến thiên 30%. Chủ thể thực hiện: Các viện nghiên cứu cơ điện tử và nhóm nghiên cứu điều khiển tự động.

  3. Tích hợp bộ lọc thích nghi trực tuyến nhằm giảm biên độ rung lắc chattering tần số cao trên 85% trong vòng 3 tháng tới mà vẫn duy trì băng thông đáp ứng của hệ thống. Chủ thể thực hiện: Kỹ sư phát triển thuật toán điều khiển truyền động.

  4. Xây dựng và chuẩn hóa thư viện khối Simulink Blockset đóng gói hoàn chỉnh trong thời gian 6 tháng, phục vụ công tác giảng dạy thực hành cho hơn 500 sinh viên và học viên cao học kỹ thuật hàng năm. Chủ thể thực hiện: Khoa Điện - Điện tử tại các trường đại học khối kỹ thuật.

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

Nội dung và kết quả của luận văn mang lại giá trị tham khảo thiết thực cho 4 nhóm đối tượng:

  1. Kỹ sư thiết kế và phát triển hệ thống điều khiển: Ứng dụng giải thuật tự chỉnh $\lambda$ vào cơ cấu tay máy robot công nghiệp, hệ thống định vị bàn trượt CNC và biến tần servo, tiết kiệm hơn 80% thời gian thử nghiệm thực tế.

  2. Học viên cao học và nghiên cứu sinh ngành Tự động hóa, Kỹ thuật Điều khiển: Kế thừa nền tảng toán học về mặt trượt $s(t)$ và kỹ thuật lập trình S-Function nâng cao trong MATLAB/Simulink để mở rộng cho các công trình nghiên cứu chuyên sâu, hỗ trợ phát triển hơn 10 hướng đề tài liên quan.

  3. Giảng viên và nhà nghiên cứu tại các trường đại học kỹ thuật: Sử dụng mô hình và mã nguồn mô phỏng của luận văn làm học liệu trực quan cho hơn 200 sinh viên mỗi năm trong các môn học Điều khiển phi tuyến và Nhận dạng hệ thống.

  4. Doanh nghiệp sản xuất xe tự hành AGV và thiết bị cơ điện tử chính xác: Tích hợp giải thuật SMC-PID tự chỉnh để tăng độ chính xác bám quỹ đạo lên mức 99% trong môi trường có nhiễu rung lắc và quán tính thay đổi liên tục.

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

  1. Thông số $\lambda$ (lamda) trong bộ điều khiển trượt đóng vai trò gì? Thông số $\lambda$ là một hằng số dương quyết định độ dốc mặt trượt và tốc độ hội tụ của sai số theo quy luật hàm mũ với hằng số thời gian $(n-1)/\lambda$. Nếu chọn $\lambda$ không phù hợp, hệ thống sẽ đáp ứng rất chậm hoặc mất ổn định; ví dụ thực nghiệm trên hàm bậc hai cho thấy đáp ứng tối ưu tại $\lambda = 30$, trong khi $\lambda = 3000$ hay $\lambda = 3$ đều khiến thời gian quá độ kéo dài.

  2. Tại sao cần tự động hóa việc xác định thông số $\lambda$ thay vì tinh chỉnh thủ công? Phương pháp thủ công yêu cầu thử sai từ 15 đến 20 lần cho mỗi mô hình đối tượng, tiêu tốn nhiều giờ làm việc và phụ thuộc vào cảm tính. Giải thuật tự động xác định chính xác giá trị $\lambda$ tối ưu chỉ sau 1 lần mô phỏng theo độ dốc đặt trước, giúp cắt giảm hơn 90% thời gian thiết kế hệ thống điều khiển.

  3. Giải thuật ước lượng $\lambda$ tối ưu hoạt động theo nguyên lý nào? Giải thuật áp dụng quy trình 3 bước: kích thích hệ thống với độ dốc lớn 10000 %/s để đo độ dốc đỉnh thực tế; sau đó đưa mức thử nghiệm 20% để xác định tỷ lệ đáp ứng; cuối cùng ngoại suy ra độ dốc tối đa khả thi (ví dụ 64 %/s trên hệ motor) và tự động cố định giá trị $\lambda$ tương ứng (như $\lambda = 360$).

  4. Bộ điều khiển đề xuất có thể áp dụng cho những thiết bị vật lý nào trong thực tế? Giải thuật ứng dụng hiệu quả cho mọi hệ thống truyền động cơ điện tử bậc hai như hệ thống phanh xe chống bó cứng ABS, động cơ servo vị trí, cánh tay robot và bàn xoay quang học. Các hệ thống này yêu cầu độ chính xác bám vị trí trên 98% và thời gian xác lập dưới 0.1 giây mà không xuất hiện độ vọt lố gây va đập cơ khí.

  5. Việc kết hợp SMC với PID trong luận văn mang lại ưu thế gì so với SMC truyền thống? Sự kết hợp SMC-PID tận dụng độ bền vững cao của chế độ trượt trước nhiễu ngoại lai, đồng thời cấu trúc PID nối tiếp giúp triệt tiêu hoàn toàn sai số xác lập và độ vọt lố. Kết quả mô phỏng cho thấy hệ SMC-PID triệt tiêu 100% vọt lố so với bộ PID kinh điển, bảo đảm an toàn cho các cơ cấu chấp hành công nghiệp.

Kết luận

Luận văn thạc sĩ của tác giả Hứa Minh Tú đã giải quyết xuất sắc bài toán tự động hóa xác định tham số cho bộ điều khiển trượt với các đóng góp trọng tâm:

  • Xây dựng thành công giải thuật tự động điều chỉnh thông số $\lambda$ trên MATLAB/Simulink thông qua khối S-Function linh hoạt.
  • Phát triển quy trình 3 bước ước lượng độ dốc cực đại, tự động tìm ra giá trị $\lambda$ tối ưu mà không cần can thiệp thủ công.
  • Triệt tiêu 100% độ vọt lố và rút ngắn đáng kể thời gian quá độ trên các hệ thống truyền động bậc hai và hệ điều khiển vị trí motor SMC-PID.
  • Khẳng định tính ổn định và khả năng thích nghi bền vững qua nhiều dải độ dốc thử nghiệm từ 10 %/s đến 50 %/s.
  • Đặt nền tảng vững chắc cho việc triển khai thuật toán điều khiển trượt thích nghi trên các phần cứng vi xử lý nhúng công nghiệp.

Kế hoạch nghiên cứu mở rộng trong 6 đến 12 tháng tới sẽ tập trung vào việc thử nghiệm phần cứng thực tế (Hardware-in-the-loop) và tối ưu hóa cho hệ thống đa biến phi tuyến cao. Độc giả, kỹ sư và các nhà nghiên cứu hãy tham khảo toàn văn luận văn và ứng dụng mô hình giải thuật để nâng cao hiệu suất cho các dự án điều khiển tự động ngay hôm nay.