I. Tổng quan bài thực hành cơ sở lý thuyết điều khiển tự động
Bài thực hành cơ sở lý thuyết điều khiển tự động đóng vai trò nền tảng, giúp chuyển hóa kiến thức lý thuyết trừu tượng thành kỹ năng ứng dụng cụ thể. Nội dung thực hành tập trung vào việc sử dụng phần mềm MATLAB và Simulink, hai công cụ mô phỏng mạnh mẽ, để phân tích và thiết kế các hệ thống điều khiển. Mục tiêu chính là khảo sát các đặc tính động học cơ bản, đánh giá tính ổn định, thiết kế bộ điều khiển PID và mô phỏng các hệ thống phức tạp trong môi trường công nghiệp. Thông qua các bài thí nghiệm, người học có thể trực quan hóa các khái niệm như hàm truyền, đáp ứng thời gian, đáp ứng tần số và hiểu rõ ảnh hưởng của các tham số đến chất lượng hệ thống. Báo cáo thực hành này trình bày chi tiết các bước tiến hành, từ việc viết mã lệnh trong MATLAB để khảo sát các khâu động học như khâu tích phân, khâu vi phân, đến việc sử dụng công cụ Quỹ đạo nghiệm số (QĐNS) để xác định giới hạn ổn định. Ngoài ra, một phần quan trọng của bài thực hành là thiết kế và tinh chỉnh bộ điều khiển PID, một bộ điều khiển được ứng dụng rộng rãi nhất trong công nghiệp, nhằm cải thiện các chỉ tiêu chất lượng như thời gian đáp ứng, độ vọt lố và sai lệch xác lập. Cuối cùng, Simulink được sử dụng để xây dựng và mô phỏng một hệ thống điều khiển nhiệt độ lò nung, minh họa cho quy trình từ lý thuyết đến ứng dụng thực tiễn một cách sinh động. Việc thành thạo các kỹ năng này là yêu cầu bắt buộc đối với kỹ sư ngành Kỹ thuật Điện - Điện tử, mở ra nhiều cơ hội trong lĩnh vực tự động hóa và điều khiển quá trình.
1.1. Mục tiêu chính của môn học điều khiển tự động
Mục tiêu cốt lõi của môn học là trang bị cho người học khả năng phân tích, đánh giá và thiết kế các hệ thống điều khiển tự động tuyến tính. Các bài thực hành được xây dựng để hiện thực hóa mục tiêu này, tập trung vào việc: 1) Nhận dạng và khảo sát các đặc tính của các khâu động học cơ bản. 2) Sử dụng các công cụ phân tích trong miền thời gian (hàm quá độ, hàm trọng lượng) và miền tần số (biểu đồ Bode, biểu đồ Nyquist). 3) Áp dụng các tiêu chuẩn ổn định đại số và tần số để đánh giá tính ổn định của hệ thống. 4) Thiết kế các bộ điều khiển hiệu quả, đặc biệt là bộ điều khiển PID, để hệ thống đạt được các chỉ tiêu chất lượng mong muốn. Quá trình này giúp củng cố kiến thức lý thuyết và phát triển tư duy giải quyết vấn đề kỹ thuật.
1.2. Vai trò của phần mềm MATLAB và Simulink trong thực hành
Phần mềm MATLAB và môi trường Simulink là những công cụ không thể thiếu trong lĩnh vực điều khiển tự động hiện đại. MATLAB cung cấp một môi trường tính toán số và lập trình mạnh mẽ, cho phép thực hiện các phân tích phức tạp như tìm nghiệm phương trình đặc trưng, vẽ quỹ đạo nghiệm số, và phân tích tần số một cách nhanh chóng và chính xác. Trong khi đó, Simulink cung cấp một giao diện đồ họa trực quan, cho phép xây dựng mô hình hệ thống dưới dạng sơ đồ khối, giống như trong thực tế. Điều này giúp đơn giản hóa việc mô phỏng các hệ thống phi tuyến hoặc có cấu trúc phức tạp, chẳng hạn như hệ thống điều khiển lò nhiệt. Sự kết hợp giữa MATLAB và Simulink giúp người học kiểm chứng lý thuyết, thử nghiệm các ý tưởng thiết kế và quan sát kết quả một cách trực quan trước khi triển khai trên phần cứng thực tế, qua đó tiết kiệm thời gian và chi phí.
II. Cách phân tích đặc tính động học hệ thống điều khiển
Phân tích đặc tính động học là bước đầu tiên và quan trọng nhất trong việc tìm hiểu và đánh giá một hệ thống điều khiển. Quá trình này giúp xác định cách hệ thống phản ứng với các tín hiệu vào khác nhau, từ đó làm cơ sở để thiết kế bộ điều khiển phù hợp. Việc phân tích thường được thực hiện trên hai miền chính: miền thời gian và miền tần số. Trong miền thời gian, các đặc tính của hệ thống được đánh giá thông qua đáp ứng của nó với các tín hiệu vào chuẩn như hàm nấc (step) và hàm xung (impulse). Hàm quá độ h(t) (đáp ứng nấc) cho thấy tốc độ đáp ứng, độ vọt lố và thời gian xác lập. Hàm trọng lượng w(t) (đáp ứng xung) thể hiện các đặc tính nội tại của hệ thống. Trong khi đó, phân tích trong miền tần số, thông qua biểu đồ Bode và biểu đồ Nyquist, lại cung cấp thông tin về độ dự trữ ổn định, băng thông và khả năng lọc nhiễu của hệ thống. Bài thực hành này sử dụng MATLAB để tự động hóa quá trình khảo sát. Ví dụ, để khảo sát khâu tích phân có hàm truyền W(s) = 5/s, các lệnh step(w) và impulse(w) được dùng để vẽ đáp ứng thời gian, trong khi bode(w) và nyquist(w) được dùng để vẽ các đặc tính tần số. Cách tiếp cận này không chỉ nhanh chóng mà còn giảm thiểu sai sót so với tính toán thủ công, cho phép tập trung vào việc diễn giải kết quả và hiểu sâu hơn về bản chất của hệ thống.
2.1. Phân tích đáp ứng thời gian Hàm quá độ và hàm trọng lượng
Đáp ứng thời gian mô tả hành vi của tín hiệu ra theo thời gian khi có một tín hiệu vào cụ thể. Hàm quá độ h(t), hay đáp ứng với tín hiệu nấc đơn vị, là một trong những chỉ số quan trọng nhất. Nó cho biết hệ thống có đạt được giá trị mong muốn hay không, tốc độ đạt được nhanh hay chậm (thời gian lên), mức độ dao động quanh giá trị xác lập (độ vọt lố), và thời gian để đi vào trạng thái ổn định (thời gian xác lập). Hàm trọng lượng w(t), hay đáp ứng với tín hiệu xung đơn vị, lại mô tả động lực học tự nhiên của hệ thống. Trong báo cáo, việc khảo sát khâu quán tính với các hằng số thời gian T khác nhau cho thấy rõ: khi T tăng, hệ thống đáp ứng chậm hơn, thể hiện qua đồ thị hàm quá độ thoải hơn.
2.2. Đánh giá hệ thống trong miền tần số Biểu đồ Bode Nyquist
Phân tích trong miền tần số cung cấp một góc nhìn khác về hệ thống, đặc biệt hữu ích cho việc đánh giá tính ổn định. Biểu đồ Bode bao gồm hai đồ thị: biên độ theo tần số và pha theo tần số. Nó cho biết độ lợi của hệ thống tại các tần số khác nhau và độ trễ pha mà hệ thống gây ra. Thông qua biểu đồ Bode, có thể xác định các thông số quan trọng như độ dự trữ biên và độ dự trữ pha, là các thước đo định lượng về mức độ ổn định. Biểu đồ Nyquist là một đồ thị phức biểu diễn đặc tính tần biên pha của hệ thống hở. Tiêu chuẩn ổn định Nyquist, dựa trên việc quan sát quỹ đạo Nyquist có bao điểm (-1, j0) hay không, là một công cụ mạnh mẽ để xác định tính ổn định của hệ thống kín. Các lệnh bode và nyquist trong MATLAB giúp vẽ các biểu đồ này một cách dễ dàng.
III. Bí quyết khảo sát tính ổn định hệ thống bằng QĐNS MATLAB
Tính ổn định là yêu cầu tiên quyết của mọi hệ thống điều khiển. Một hệ thống không ổn định sẽ có đáp ứng phân kỳ theo thời gian, dẫn đến phá hủy thiết bị và gây nguy hiểm. Phương pháp Quỹ đạo nghiệm số (QĐNS) là một công cụ đồ họa mạnh mẽ trong MATLAB để phân tích tính ổn định và chất lượng của hệ thống hồi tiếp âm khi một tham số (thường là độ lợi k) thay đổi. QĐNS là quỹ tích của các nghiệm của phương trình đặc trưng của hệ thống kín khi k biến thiên từ 0 đến vô cùng. Vị trí của các nghiệm này trên mặt phẳng phức quyết định trực tiếp đến hành vi của hệ thống: nếu tất cả các nghiệm nằm bên trái trục ảo, hệ thống ổn định. Nếu có bất kỳ nghiệm nào nằm bên phải trục ảo, hệ thống không ổn định. Nếu nghiệm nằm trên trục ảo, hệ thống ở biên giới ổn định. Báo cáo thực hành đã ứng dụng phương pháp này để tìm giá trị hệ số khuếch đại giới hạn (k_gh). Bằng cách vẽ QĐNS và tìm giao điểm của nó với trục ảo, giá trị k_gh được xác định. Từ đồ thị, báo cáo chỉ ra k_gh = 421. Với k < k_gh, hệ thống ổn định; với k > k_gh, hệ thống mất ổn định. Phương pháp này không chỉ cho biết hệ thống ổn định hay không mà còn cho phép thiết kế hệ thống có các chỉ tiêu chất lượng mong muốn bằng cách chọn vị trí nghiệm phù hợp trên quỹ đạo.
3.1. Hướng dẫn vẽ Quỹ đạo nghiệm số Root Locus trong MATLAB
Việc vẽ Quỹ đạo nghiệm số trong MATLAB được thực hiện chủ yếu qua lệnh rlocus(G), trong đó G là hàm truyền vòng hở của hệ thống. Sau khi vẽ, có thể sử dụng lệnh rlocfind(G) để chọn một điểm bất kỳ trên quỹ đạo. MATLAB sẽ trả về giá trị độ lợi k tương ứng với điểm đó cùng với tọa độ của các cực của hệ kín. Lưới tọa độ cực (sgrid) có thể được thêm vào đồ thị để dễ dàng xác định các vùng có cùng hệ số tắt (damping ratio) hoặc tần số dao động tự nhiên, giúp cho việc thiết kế bộ điều khiển trở nên trực quan hơn.
3.2. Cách xác định hệ số khuếch đại giới hạn k_gh từ QĐNS
Hệ số khuếch đại giới hạn (k_gh) là giá trị của độ lợi k mà tại đó hệ thống bắt đầu mất ổn định. Trên đồ thị QĐNS, đây chính là giá trị k tại các điểm mà quỹ đạo cắt trục ảo. Để tìm giá trị này, có thể sử dụng công cụ rlocfind và di con trỏ chuột đến chính xác giao điểm của quỹ đạo và trục ảo. Báo cáo thực hành đã áp dụng cách này cho hệ thống có hàm truyền vòng hở G(s) = k / ((s+3)(s^2+8s+20)) và xác định được k_gh = 421. Đây là thông tin cực kỳ quan trọng trong thực tế, giúp giới hạn phạm vi hoạt động an toàn cho hệ thống.
3.3. Tối ưu chất lượng hệ thống Độ vọt lố và thời gian xác lập
Ngoài việc xác định tính ổn định, QĐNS còn là công cụ hiệu quả để tối ưu chất lượng hệ thống. Các chỉ tiêu như độ vọt lố (overshoot) và thời gian xác lập (settling time) có liên quan trực tiếp đến vị trí của cặp cực phức quyết định. Báo cáo đã minh họa việc tìm giá trị k để hệ thống có độ vọt lố POT = 25% (tương ứng với k = 140) hoặc thời gian xác lập t_s = 4s (tiêu chuẩn 2%, tương ứng với k = 43). Bằng cách lựa chọn các điểm mong muốn trên QĐNS, kỹ sư có thể tìm ra giá trị k phù hợp để hệ thống vừa ổn định, vừa đáp ứng được các yêu cầu về chất lượng động học.
IV. Hướng dẫn thiết kế bộ điều khiển PID tối ưu cho hệ thống
Bộ điều khiển Tỷ lệ - Tích phân - Vi phân (bộ điều khiển PID) là bộ điều khiển được sử dụng phổ biến nhất trong các ứng dụng công nghiệp nhờ cấu trúc đơn giản, hiệu quả và dễ tinh chỉnh. Mục tiêu của việc thiết kế bộ điều khiển PID là cải thiện đáp ứng của hệ thống, bao gồm: giảm thời gian đáp ứng, giảm thiểu độ vọt lố và triệt tiêu sai lệch tĩnh (steady-state error). Bài thực hành số 3 đã khảo sát chi tiết ảnh hưởng của từng thành phần (P, I, D) lên chất lượng hệ thống. Thành phần Tỷ lệ (P - Kp) giúp tăng tốc độ đáp ứng nhưng thường để lại sai lệch tĩnh. Thành phần Tích phân (I - Ki) có tác dụng triệt tiêu sai lệch tĩnh nhưng có thể làm giảm tính ổn định và tăng độ vọt lố. Thành phần Vi phân (D - Kd) giúp cải thiện tính ổn định, giảm độ vọt lố và tăng tốc độ đáp ứng bằng cách dự đoán xu hướng của sai lệch. Báo cáo đã tiến hành mô phỏng hệ thống với các bộ điều khiển P, PI, PD và PID. Kết quả cho thấy, khi chỉ dùng bộ P, hệ thống có sai lệch tĩnh lớn. Khi thêm khâu I (bộ PI), sai lệch tĩnh bị triệt tiêu nhưng độ vọt lố tăng. Khi sử dụng bộ PID được tinh chỉnh hợp lý, hệ thống đạt được chất lượng tốt nhất: đáp ứng nhanh, độ vọt lố rất nhỏ (0.0076) và sai lệch tĩnh gần như bằng không (0.0004). Quá trình này minh họa rõ sức mạnh của việc kết hợp cả ba thành phần để đạt được sự cân bằng tối ưu giữa tốc độ, độ ổn định và độ chính xác.
4.1. Phân tích đáp ứng hệ thống với bộ điều khiển P PI và PD
Việc phân tích riêng lẻ từng bộ điều khiển cho thấy rõ vai trò của mỗi thành phần. Bộ điều khiển P (Proportional) cung cấp một hành động điều khiển tỷ lệ với sai lệch hiện tại. Kết quả mô phỏng cho thấy đáp ứng nhanh nhưng không thể đạt chính xác giá trị đặt, gây ra sai lệch tĩnh. Bộ điều khiển PI (Proportional-Integral) thêm thành phần tích phân của sai lệch, giúp loại bỏ hoàn toàn sai lệch tĩnh nhưng có thể làm hệ thống đáp ứng chậm hơn và dao động nhiều hơn. Ngược lại, bộ điều khiển PD (Proportional-Derivative) sử dụng đạo hàm của sai lệch để chống lại sự thay đổi đột ngột, giúp giảm vọt lố và cải thiện độ ổn định, nhưng không loại bỏ được sai lệch tĩnh.
4.2. Tinh chỉnh thông số Kp Ki Kd để đạt chất lượng mong muốn
Quá trình tinh chỉnh các thông số Kp, Ki, Kd là bước quyết định để bộ điều khiển PID hoạt động hiệu quả. Không có một bộ thông số nào là hoàn hảo cho mọi hệ thống. Việc lựa chọn các giá trị này là một sự đánh đổi. Tăng Kp giúp hệ thống nhanh hơn nhưng dễ gây mất ổn định. Tăng Ki loại bỏ sai lệch tĩnh nhanh hơn nhưng tăng vọt lố. Tăng Kd giúp giảm vọt lố nhưng làm hệ thống nhạy cảm với nhiễu tần số cao. Các phương pháp tinh chỉnh kinh điển như Ziegler-Nichols hoặc các công cụ tự động trong MATLAB (PID Tuner) có thể được sử dụng để tìm ra bộ thông số khởi đầu tốt, sau đó tinh chỉnh thủ công để đạt được hiệu suất tối ưu theo yêu cầu cụ thể của ứng dụng.
4.3. Đánh giá hiệu suất bộ điều khiển PID qua các chỉ tiêu chất lượng
Để đánh giá hiệu suất của bộ điều khiển PID một cách khách quan, cần sử dụng các chỉ tiêu chất lượng định lượng. Lệnh stepinfo trong MATLAB là một công cụ hữu ích, cung cấp một cấu trúc dữ liệu chứa các thông số quan trọng như: RiseTime (thời gian lên), SettlingTime (thời gian xác lập), Overshoot (độ vọt lố), và Peak (giá trị đỉnh). Dựa vào các thông số này, có thể so sánh hiệu quả của các bộ điều khiển khác nhau hoặc các bộ thông số PID khác nhau. Ví dụ, trong báo cáo, bộ điều khiển PID cuối cùng đạt được RiseTime: 0.1657s và Overshoot: 0.0076, những con số này cho thấy một đáp ứng vừa nhanh vừa ổn định.
V. Ứng dụng Simulink mô phỏng hệ thống điều khiển lò nhiệt
Việc ứng dụng lý thuyết điều khiển vào các bài toán thực tế là mục tiêu cuối cùng của môn học. Bài thực hành số 4 sử dụng Simulink, một môi trường mô phỏng đồ họa của MATLAB, để xây dựng và khảo sát một hệ thống điều khiển nhiệt độ lò nung. Ưu điểm của Simulink là cho phép người dùng xây dựng mô hình hệ thống bằng cách kéo và thả các khối chức năng, giúp trực quan hóa cấu trúc và dòng tín hiệu một cách rõ ràng. Mô hình lò nhiệt được xây dựng từ các khối hàm truyền (Transfer Fcn), khối tạo tín hiệu đầu vào (Step) và khối hiển thị kết quả (Scope). Đầu tiên, hệ thống được khảo sát ở chế độ vòng hở để quan sát đặc tính tự nhiên của lò nhiệt. Kết quả cho thấy quá trình tăng nhiệt độ diễn ra rất chậm và không có khả năng tự ổn định tại một giá trị đặt trước. Đây là một thách thức điển hình trong điều khiển quá trình. Sau đó, một cấu trúc điều khiển vòng kín được xây dựng bằng cách thêm vào một bộ điều khiển PID. Việc mô phỏng hệ thống vòng kín cho thấy sự cải thiện vượt bậc: nhiệt độ lò tăng nhanh đến giá trị mong muốn, dao động được kiểm soát và hệ thống giữ ổn định nhiệt độ trước các tác động. Bài thực hành này là một minh chứng xuất sắc về cách Simulink có thể được sử dụng để thiết kế, kiểm tra và tối ưu hóa các hệ thống điều khiển phức tạp trước khi triển khai trong thực tế, giúp tiết kiệm chi phí và đảm bảo an toàn.
5.1. Xây dựng mô hình lò nhiệt vòng hở trên Simulink
Mô hình hệ thống lò nhiệt vòng hở được xây dựng trong Simulink bằng cách kết nối các khối theo đúng sơ đồ. Khối Step được dùng để mô phỏng việc cấp nguồn (công suất 100%). Tín hiệu này đi qua hai khối Transfer Fcn nối tiếp nhau, đại diện cho động học của bộ gia nhiệt (1/(30s+1)) và của buồng lò (300/(120s+1)). Khối Scope được kết nối ở đầu ra để quan sát và vẽ đồ thị quá trình thay đổi nhiệt độ theo thời gian. Việc thiết lập thời gian mô phỏng (Stop time = 600s) là cần thiết để quan sát toàn bộ quá trình quá độ của hệ thống có hằng số thời gian lớn như lò nhiệt.
5.2. Phân tích quá trình quá độ của hệ thống lò nung
Sau khi chạy mô phỏng vòng hở, đồ thị trên Scope cho thấy đáp ứng của hệ thống. Đặc điểm nổi bật là quá trình tăng nhiệt độ diễn ra rất từ từ, thể hiện quán tính nhiệt lớn của lò nung. Hệ thống cần một khoảng thời gian dài để đạt đến trạng thái xác lập. Hơn nữa, ở chế độ vòng hở, không có cơ chế nào để điều chỉnh công suất gia nhiệt khi nhiệt độ đạt đến giá trị mong muốn. Điều này có nghĩa là hệ thống không thể tự ổn định nhiệt độ tại một điểm đặt cụ thể, chứng tỏ sự cần thiết phải có một vòng điều khiển hồi tiếp.
5.3. Tích hợp bộ điều khiển PID để cải thiện hiệu suất
Để cải thiện hiệu suất, một vòng lặp hồi tiếp được thêm vào mô hình Simulink. Vòng lặp này bao gồm khối Sum để tính toán sai lệch giữa tín hiệu đặt và tín hiệu nhiệt độ thực tế. Sai lệch này được đưa vào khối PID Controller. Đầu ra của khối PID sẽ điều khiển công suất cấp cho lò nhiệt. Bằng cách tinh chỉnh các thông số Kp, Ki, Kd của khối PID, hệ thống có thể được tối ưu hóa để đạt được đáp ứng nhanh, không vọt lố và giữ nhiệt độ ổn định chính xác tại giá trị mong muốn. Đây là mô hình hoạt động của hầu hết các hệ thống điều khiển nhiệt độ trong công nghiệp.
