Xem mẫu

  1. Trường Đại học Điện Lực Khoa công nghệ tự động ĐỒ ÁN MÔN HỌC: LÝ THUYẾT ĐIỀU KHIỂN TỰ ĐỘNG Giáo viên hướng dẫn: PHẠM THỊ HƯƠNG SEN Sinh viên thực hiện : NGUYỄN TRỌNG DŨNG LỚP : Đ6-CNTĐ HÀ NỘI 2013
  2. MỤC LỤC Trang CHƯƠNG 1: ĐẶC TÍNH THỜI GIAN CỦA MỘT KHÂU……………………. ……4 1.1. Định nghĩa…………………………………………………………………...........4 1.2. Khảo sát các đường đặc tính thời gian của là lò điện trở có hàm truyền đạt………………………………………………………………………………………4 CHƯƠNG 2: QUỸ ĐẠO NGHIỆM…………………………………………………….7 2.1. Định nghĩa………………………………………………………………………..7 2.2. Quy tắc vẽ quỹ đạo nghiệm số……………………………………………...........7 2.3. Ứng dụng vẽ quỹ đạo nghiệm số của hệ kín có phản hồi âm đơn vị ….............9 CHƯƠNG 3: Thiết kế bộ điều khiển cho lò điện trở theo các luật điều khiển P,PI,PID. Xác định tham số bộ điều khiển theo ba phương pháp khác nhau để hệ có chất lượng tốt nhất............................................................................................................11 3.1. Thiết kế bộ điều khiển theo luật P………………………………………………11 3.2. Thiết kế bộ điều khiển theo luật PI………………………………………..........12 3.3. Thiết kế bộ điều khiển theo luật PID…………………………………………...13 3.4. Các phương pháp thiết kế bộ điều khiển………………………………….........16 3.4.1. Thiết kế bộ điều khiển P,PI,PID bằng phương pháp Zeigler – Nichols ……..........16 3.4.2. Thiết kế bộ điều khiển P,PI,PID bằng Simulink……………………………..22
  3. 3.4.3. Thiết kế bộ diều khiển P, PI, PID bằng phương pháp sisotool……………...28 CHƯƠNG 4: TÌM HIỂU VỀ HỆ ĐIỀU KHIỂN LÒ ĐIỆN TRỞ TRONG THỰC TẾ………………………………………………………………………………………..34 4.1 Nguyên lý làm việc của lò điện trở……………………………………………… 34 4.2 Cấu tạo của lò điện trở……………………………………………………….......35 LỜI NÓI ĐẦU Ngày nay ngành tự động hóa đã trở thành một vấn đề thiết yếu trong ngành công nghiệp.Để thiết kế được các mô hình tự động hóa trong nhà máy công nghiệp thì người thiết kế cần nắm được các kiến thức về Lý thuyết điều khiển tự động – bộ môn cơ bản của ngành tự động hóa. Một trong các kỹ năng mà ng ười h ọc ph ải có sau khi học xong bộ môn này là phải nhận dạng và ổn định các mô hình. Trong đồ án này em đã biết cách khảo sát các đường đặc tính thời gian, vẽ được quỹ đạo nghiệm số của hệ kín phản hồi thông qua Maclab t ừ đó xác định các thông số để hệ thống ổn định rồi từ đó thiết kế các bộ điều khiển P,PI,PID đ ể nâng cao chất lượng đầu ra của hệ thống Trong quá trình thực hiện đồ án này em đã nhận được rất nhiều sự khuyến khích và góp ý từ các ban cũng như thầy cô, đặc biệt là cô Phạm Thị Hương Sen – Giáo viên bộ môn của khoa công nghệ tự động trường Đại học Điện Lực. Với những kiến thức và hiểu biết còn hạn chế, em rất mong nhận được nhiều hơn nữa những sự đóng góp , bổ xung ý kiến của cô và các bạn để
  4. cho đồ án này hoàn thiện hơn,giúp em có kiến thức vững chắc để có thể học tập và nghiên cưu sâu hơn trong ngành công nghệ tự động. Em xin chân thành cảm ơn!
  5. Đềbài: Cho đối tượng cần điều khiển là lò điện trở có hàm truyền đạt: W(s) = Yêu cầu: 1. Cho K = 10, khảo sát các đường đặc tính theo thơi gian của lò điện trở. Nêu nhận xét. 2. Vẽ quỹ đạo nghiệm số của hệ kín có phản hồi âm đơn vị. Dựa vào quỹ đạo nghiệm số, tìmKgh để hệ thống ổn định, chỉ rõ giá trị này trên quỹ đạo nghiệm số. - Tìm K để hệ có tần số dao động tự nhiên = 6 - Tìm K để hệ có hệ số tắt dần0.8 - Tìm K để hệ có độ quá điều chỉnh 25% - Tìm K để hệ có thời gian quá độ là 200s. 3. Thiết kế hệ điều khiển cho lò điện trở theo các luật điều khiển: P, PI, PID. Tiến hành xác định tham số bộ điều khiển theo ba phương pháp khác nhau để hệ có chất lượng tốt nhất. 4. Tìm hiểu về hệ điều khiển nhiệt độ lò điện trở trong thực tế.
  6. CHƯƠNG 1: ĐẶC TÍNH THỜI GIAN CỦA MỘT KHÂU 1.1. Định nghĩa Đặc tính thời gian của hệ thống mô tả sự thay đổi tín hiệu đầu ra của hệ thống khi tín hiệu đầu vào là hàm xung đơn vị hay hàm nấc đơn vị. Hàm quá độ của một khâu Hàm quá độ của một khâu là phản ứng của khâu đó với tín hiệu vào 1(t). Kí hiệu : h(t) Biểu thức : h(t)= L-1{} Hàm trọng lượng của một khâu Hàm trọng lượng của một khâu là phản ứng của khâu đối với tín hiệu vào(t). Kí hiệu : (t). Biểu thức : (t)= L-1{W(s)} 1.2. Khảo sát các đường đặc tính thời gian của là lò điện trở có hàm truyền đạt: W(s) = Sử dụng phần mềm maclab để khảo sát các đường đặc tính thời gjan. Khai báo đối tượng khảo sát *) Hàm quá độ h(t): Mở cửa sổ Command window gõ lệnh: >> num=[10]; >> den=[400 1]; >> W=tf(num,den,'inputdelay',20); >> impulse(W); >> step(W);
  7. *) Hàm quá độ h(t): Step Res pons e 10 9 8 7 6 Am plitude 5 4 3 2 1 0 0 500 1000 1500 2000 25 00 Time (s e c onds ) Nhận xét: +) Do bậc của mẫu lớn hơn bậc của tử nên đường h(t) xuất phát từ gốc tạo độ O. +) Do = =10 nên đường h(t) tiến đến 10. +) Thời gian lên: 880 (s). +) Thời gian xác lập: 1593 (s).
  8. *) Hàm trọng lượng (t): Impu ls e Re s po ns e 0 .02 5 0.0 2 0 .01 5 A m plitud e 0.0 1 0 .00 5 0 0 5 00 10 00 15 0 0 20 0 0 25 00 Time ( s e c on ds ) Nhận xét: +) Đường g(t) xuất phát từ gốc tạo độ O và tiến dần về 0 như hình vẽ. +) Thời gian xác lập: 1653 (s).
  9. CHƯƠNG 2: QUỸ ĐẠO NGHIỆM SỐ 2.1. Định nghĩa Quỹ đạo nghiệm số là quỹ đạo tạo ra từ các nghiệm của phương trình đặc tính của hệ thống khi có một thông số nào đó của hệ thay đổi từ 0 đến . Bằng cách quan sát quỹ đạo nghiệm số thì ta có thể nhận thấy quỹ đạo nghiệm số nào ở bên trái trục ảo thì hệ thống sẽ ổn định, còn những QĐNS nằm ở bên phải trục ảo thì hệ thống không ổn định. Từ đó ta có thể xác định được khoảng thông số thay đổi để hệ thống ổn định. Phương pháp này thường dùng cho hệ số biến đổi là hệ số khuếch đại của hệ thống. 2.2. Quy tắc vẽ quỹ đạo nghiệm số Để vẽ QĐNS, trước tiên ta phải biến đổi tương đương phương trình đặc tính về dạng : 1 + =0 Trong đó : K là thông số không thay đổi. Đặt G0(s) = Gọi n là số cực của G0(s) , m là số zero của G0(s). Ta có điều kiện biên độ và điều kiện pha: |G0(s)| = 1 Điều kiện biên độ
  10. - Quy tắc 4: Một điểm trên trục thực thuộc về quỹ đạo nghiệm số nếu tổng số cực và zero của G0(s) bên phải nó là một số lẻ. - Quy tắc 5: Góc tạo bởi các đường tiệm cân của quỹ đạo nghiệm số với trục thực xác định bởi: (với l=0, - Quy tắc 6: Giao điểm giữa các tiệm cận của quỹ đạo ngiệm số với trục thực là điểm A có tọa độ xác định bởi: OA= Trong đó: pi và zj là các cực và zero của G0(s). - Quy tắc 7: Điểm tách nhập( nếu có ) của quỹ đạo nghiệm số nằm trên trục thực và là nghiệm của phương trình : - Quy tắc 8: Giao điểm của quỹ đạo nghiệm số với trục ảo có thể xác định bằng một trong hai cách sau: +) Áp dụng tiêu chuẩn Routh- Hurwitz. +) Thay s= vào phương trình đặc tính, cân bằng phần thực và phần ảo sẽ tìm được giao điểm với trục ảo và giá trị K. - Quy tắc 9: Góc xuất phát của các quỹ đạo nghiệm số tại cực phức pi được xác định bởi - Quy tắc 10: Tổng các nghiệm là hằng số khi K thay đổi từ 0 đến +. - Quy tắc 11: Hệ số khuếch đại dọc theo quỹ đạo nghiệm số có thể xác định điều kiện biên độ
  11. Từ quỹ đạo nghiệm số ta thấy quỹ đạo nghiệm số của hệ thống nằm ở bên trái trục ảo, do đó hệ thống ổn định. 2.3. Ứng dụng vẽ quỹ đạo nghiệm số của hệ kín có phản hồi âm đơn vị Mở cửa sổ Command window gõ lệnh: >> L=20; >> n=3; >> [num,den]=pade(L,n); >> b=tf(num,den); >> Wđt=tf(1,[400 1])*b; >> rlocus(Wđt); Ta được quỹ đạo nghiệm số như hình vẽ: Root Locus 1 0.8 Sy s tem: w Gain: 307 0.6 Pole: 0.000178 + 0.768i Damping: -0.000231 Ov ershoot (% ): 100 Frequenc y (rad/s): 0.768 0.4 -1 ) 0.2 Im aginary Ax is (s ec onds 0 -0.2 -0.4 -0.6 -0.8 -1 -0.5 0 0.5 1 1.5 2 2.5 Real A x is (s ec onds -1) Từ đồ thị cho ta:
  12. 1. Điểm cực: -0.232; -; -1.84 + 0.175j ; -1.84 - 0.175j . 2. Quỹ đạo nghiệm có 4 nhánh. 3. Điểm zero : 0.232; -1.84 + 0.175j ; -1.84 - 0.175j . 4. Giao điểm của quỹ đạo nghiệm với trục ảo : -0.76j; -0.08j; 0.08j; 0.76j. 5. Từ QĐNS =>> Kgh = 307 *) Dựa vào quỹ đạo nghiệm số ta có: - Với tần số dao động tự nhiên = 6 Root Locus 1 0.8 0.6 System: w Gain: 2.93e+03 0.4 Pole: 6 Imaginary Axis (seconds-1) Damping: -1 0.2 Overshoot (%): 0 Frequency (rad/s): 6 0 -0.2 -0.4 -0.6 -0.8 -1 0 1 2 3 4 5 6 -1 Real Axis (seconds ) => K = 2930 - Hệ có hệ số tắt dần 0.8 =>K = 8.815 - Hệ có độ quá điều chỉnh 25%
  13. System: w Root Locus 1 Gain: 91 Pole: -0.225 + 0.51i 0.8 Damping: 0.404 Overshoot (%): 25 0.6 Frequency (rad/s): 0.558 0.4 Imaginary Axis (seconds-1) 0.2 0 -0.2 -0.4 -0.6 -0.8 -1 -0.5 0 0.5 1 1.5 2 2.5 -1 Real Axis (seconds ) => K = 91 - Hệ có thời gian quá độ là 200s => K = 4.7. CHƯƠNG 3: Thiết kế bộ điều khiển cho lò điện trở theo các luật điều khiển P,PI,PID. Xác định tham số bộ điều khiển theo ba phương pháp khác nhau để hệ có chất lượng tốt nhất. • Thiết kế bộ điều khiển cho lò điện trở theo các luật điều khiển • Đối tượng điều khiển là một khâu quán tính bậc nhất và khâu trễ có hàm truyền: 3.1. Thiết kế bộ điều khiển theo luật P - Luật P là luật điều khiển tỉ lệ tạo ra tín hiệu điều khiển u(t) tỉ lệ với tín hiệu sai lệch e(t).
  14. * Phương trình vi phân : u(t) = * Hàm truyền : Trong đó Kp là hệ số khuếch đại quy luật. Theo tính chất của khâu khuếch đại( hay khâu tỷ lệ) ta thấy tín hiệu ra của khâu luôn luôn trùng pha với tín hiệu vào. Điều này nói lên ưu điểm của khâu khuếch đại là độ tác động nhanh. Vì vậy, trong công nghiệp,quy luật tỉ lệ làm việc ổn định với mọi đối tượng. Tuy nhiên, nhược điểm cơ bản của khâu tỉ lệ là khi sử dụng với các đối tượng tĩnh, hệ thồng điều khiển luôn tồn tại sai lệch tĩnh. Để giảm giá trị sai lệch tĩnh thì phải tăng hệ số khuếch đại nhưng khi đó, tính dao động của hệ thống sẽ tăng lên và có thể làm hệ thống mất ổn định. Quy luật tỉ lệ thường được dùng cho những hệ thống cho phép tồn tại sai lệch tĩnh. K càng lớn thì sai số xác lập càng nhỏ. Nếu tăng K thì rõ rang sai lệch tĩnh giảm nhưng lại có biên độ dao động tăng quá, khi đó hệ thống sẽ mất ổn định vì vậy phải lựa chọn thông số cho phù hợp. 3.2. Thiết kế bộ điều khiển theo luật PI Để hệ thống vừa có tác động nhanh, vừa triệt tiêu được sai lệch dư, người ta kết hợp quy luật tỉ lệ với quy luật tích phân để tạo ra quy luật tỉ lệ - tích phân.
  15. Luật điều khiển PI là cấu trúc ghép song song của khâu P và khâu I. Tín hiệu ra của bộ PI là tổng tín hiệu ra của hai khâu thành phần. * Phương trình vi phân * Hàm truyền : Đồ thị bode của khâu PI Về tốc độ tác động thì quy luật PI chậm hơn quy luật tỉ lệ nhưng nhanh hơn quy luật tích phân. Trong thực tế, quy luật điều khiển PI được sử dụng khá rộng rãi và đáp ứng được chất lượng cho hầu hết các quá trình công nghệ, Tuy nhiên, do có thành phần tích phân nê độ tác động của quy luật bị chậm đi. Vì vậy, nếu đối tượng có nhiễu tác động liên tục mà hệ thống lại đòi hỏi độ chính xác cao thì quy luật PI không đáp ứng được. 3.3. Thiết kế bộ điều khiển theo luật PID Để tăng tốc tác động của quy luật PI, trong thành phần của nó người ta ghép thêm thành phần vi phân và nhận được quy luật điều khiển tỉ lệ tích phân. Có thêm thành phần vi phân làm tăng tốc tác động cho hệ thống. Luật điều khiển PID được tạo bằng cách ghép song song ba khâu: P,I và D.
  16. * Phương trình vi phân: Khâu điều chỉnh PID có hàm truyền đạt: Về tốc độ, quy luật PID còn có thể nhanh hơn cả quy luật tỉ lệ. Nói tóm lại, quy luật PID hoàn hảo nhất.Nó đáp ứng nhu cầu chất lượng của hầu hết các quy trình công nghệ nhưng việc hiệu chnhr các tham số của nó phức tạp, đòi hỏi người sử dụng ở những nơi cần thiết, khi quy luật PI không đáp ứng được yêu cầu về chất lượng điều chỉnh. Đồ thị Bode của khâu PID Nhận xét: Đây là quy luật điều khiển hoàn hảo nhất, nhanh và chính xác độ sai số xác lập nhỏ, độ quá hiệu chỉnh có thể điều chỉnh được. Nhưng nhạy cảm vói nhiễu và việc điều chỉnh 3 thông số trên rất phức tạp. Trên thực tế bộ điều khiển PID có thể được tạo ra từ các mạch từ các mạng điện. điện tử… hoặc tạo ra từ các bộ điều khiển mềm trong máy tính. * Sơ đồ khối hệ thống điều khiển như sau:
  17. 1 Gain 1 1 400s+1 1 s Transport Scope Add Transfer Fcn Step Gain1 Integrato r Delay 1 du/dt Gain2 Derivative Ảnh hưởng của các tham số KP, KI, KD, đối với các chỉ tiêu chất lượng được thể hiện qua bảng sau: Thay đổi tham số Chỉ tiêu chất lượng Tăng Tăng Tăng Thời gian đáp Giảm Giảm ít Giảm ít ứng Thời gian quá Thay đổi ít Giảm Giảm độ Độ quá hiệu Tăng Tăng Giảm ít chỉnh Hệ số tắt dần Thay đổi ít Tăng Giảm Sai lệch tĩnh Giảm Triệt tiêu Thay đồi ít Tín hiệu điều Tăng Tăng Tăng khiển Độ dự trữ ổn Giảm Giảm Tăng định Bền với Giảm Thay đổi ít Giảm nhiễu đo
  18. 3.4. Các phương pháp thiết kế bộ điều khiển 3.4.1. Thiết kế bộ điều khiển P,PI,PID bằng phương pháp Zeigler – Nichols Đây là phương pháp thông dụng nhất để chon thông số cho bộ điều khiển PID thương mại hiện nay. Phương pháp này dựa vào thực nghiệm để thiết kế bộ điều khiển P,PI,PID bằng cách chọn thông số bộ điều khiển PID tùy theo đặc điểm của đối tượng. Khi đó ta có thể xác định các thông số của bộ điều khiển P,PI,PID theo bảng sau: Thông số KP TI TD Bộ ĐK P T/(L*K) 0 PI 0.9*T/(L*K) L/0.3 0 PID 1.2*T/(L*K) 2*L L/2 Với T = 400; L = 20; K=1: Thông số KP TI TD Bộ ĐK P 20 0 PI 18 66.67 0 PID 24 40 10  Thiết kế khâu P: Mở cửa sổ Command window gõ lệnh: >> L = 20; >> n = 3;
  19. >> [num den]=pade(L,n); >> Wtre=tf(num,den); >> Wdt=tf([1],[400 1])*Wtre; >> Wpid=20; >> W=feedback(Wdt*Wpid,1); >>step(W); Step Response 1.6 1.4 System: Wdt Peak amplitude: 1.42 1.2 Overshoot (%): 49.1 System: Wdt At time (seconds): 58.9 Settling time (seconds): 249 1 0.8 Amplitude 0.6 0.4 0.2 0 -0.2 0 50 100 150 200 250 300 350 400 Time (seconds) - Từ đồ thị ta thấy độ quá điều chỉnh là 49,1% > 10% và thời gian quá độ là 249s không thể chấp nhận thông số của bộ điều khiển này.
  20. - Ta giảm KP = 10. Step Response 1.2 System W: dt 1 Settling tim (seconds): 117 e System W: dt 0.8 Peak am plitude: 0.953 Overshoot (% 4.8 ): At tim (seconds): 89 e A plitude 0.6 m 0.4 0.2 0 -0.2 0 20 40 60 80 100 120 140 160 180 T e (seconds) im - Từ đồ thị ta thấy độ quá điều chỉnh là 4.8% < 10% và thời gian quá độ là 117s nên có thể chấp nhận thông số của bộ điều khiển này.  Thiết kế khâu PI: Mở cửa sổ Command window gõ lệnh: >> L=20;n=3; >> [num,den]=pade(L,n); >> w1=tf(num,den); >> kp=18; >> ti=66.67; >> w2=tf(1,[400 1])*w1; >> w3=tf(1,[ti 0]); >> w4=kp*(1+w3); >> w5=w2*w4; >> w6=feedback(w5,1); >> step(w6);
nguon tai.lieu . vn