Giáo trình kỹ thuật điều khiển 7

Chia sẻ: Cindy Cindy | Ngày: | Loại File: PDF | Số trang:10

Thêm vào BST

Báo xấu

170
lượt xem 39
download

Download Vui lòng tải xuống để xem tài liệu đầy đủ

Mô tả hiệu suất trong miền thời gian: Các mô tả định lượng yêu cầu về hiệu suất trong miền thời gian là những chỉ số quan trọng vì các hệ thống điều khiển đều là những hệ thống trong miền thời gian. Điều đó có nghĩa là, hiệu suất nhất thời hay hiệu suất theo thời gian của hệ thống là mối quan tâm chủ yếu đối với các hệ thống điều khiển. Điều cần làm đầu tiên là xác định xem hệ thống có ổn định hay không. ...

Chủ đề:

Bình luận(0) Đăng nhập để gửi bình luận!

Lưu

Nội dung Text: Giáo trình kỹ thuật điều khiển 7

hiện nhiệm vụ được thiết kế cho nó tốt tới mức nào? 5.2. Mô tả hiệu suất trong miền thời gian Các mô tả định lượng yêu cầu về hiệu suất trong miền thời gian là những chỉ số quan trọng vì các hệ thống điều khiển đều là những hệ thống trong miền thời gian. Điều đó có nghĩa là, hiệu suất nhất thời hay hiệu suất theo thời gian của hệ thống là mối quan tâm chủ yếu đối với các hệ thống điều khiển. Điều cần làm đầu tiên là xác định xem hệ thống có ổn định hay không. Chúng ta sẽ nghiên cứu các kỹ thuật phân tích tính ổn định của hệ thống ở các chương sau. Nếu hệ thống ổn định, đáp ứng của hệ thống với một tín hiệu vào nhất định sẽ cung cấp một số số đo của hiệu suất. Tuy nhiên, do tín hiệu vào thực sự của hệ thống thường khó xác định, các tín hiệu vào thử (test input signal) chuẩn thường được sử dụng. Phương pháp này rất hữu ích bởi vì tồn tại một tương quan giữa đáp ứng của hệ thống với một tín hiệu vào thử chuẩn và khả năng của hệ thống khi hoạt động ở những điều kiện bình thường. Hơn nữa, việc sử dụng tín hiệu vào chuẩn cho phép người thiết kế so sánh nhiều thiết kế khác nhau. Nhiều hệ thống điều khiển có các tín hiệu vào cũng tương tự các tín hiệu vào thử chuẩn. Các tín hiệu thử vào chuẩn thường được sử dụng là: (1) tín hiệu nhảy bậc (step), (2) tín hiệu dốc (ramp), (3) tín hiệu parabol và (4) tín hiệu xung đơn vị (unit impulse). Định nghĩa và biến đổi Laplace của các tín hiệu này được trình bày trong bảng sau. Bảng 5.1. Các tín hiệu vào thử chuẩn Tín hiệu r(t) R(s) ⎧ K khi t ≥ 0 K ⎨ Nhảy bậc ⎩ 0 khi t < 0 s ⎧ Kt khi t ≥ 0 K ⎨ Dốc s2 ⎩ 0 khi t < 0 ⎧ Kt 2 khi t ≥ 0 ⎪ 2K ⎨ Parabol s3 ⎪ 0 khi t < 0 ⎩ δ (t ) = lim f ε (t ) , ở đó ε →0 ⎧1 Xung đơn vị 1 khi 0 ≤ t ≤ ε ⎪ f ε (t ) = ⎨ ε ⎪0 khi t < 0 hay t > ε ⎩ ⎧ Kt n khi t ≥ 0 ⎪ n! K ⎨ Dạng tổng quát s n +1 ⎪ 0 khi t < 0 ⎩ 67
Với một hệ thống vòng hở có tín hiệu vào là r(t), tín hiệu ra là c(t) và hàm chuyển là G(s): C(s) = G(s)R(s), nếu tín hiệu vào là hàm xung đơn vị δ(t), chúng ta sẽ có R(s) = 1, nghĩa là C(s) = G(s) hay c(t) = g(t). Mối quan hệ giữa tín hiệu ra và tín hiệu vào của một hệ thống điều khiển phản hồi được biểu diễn bằng phương trình: G ( s) C (s) = (5.1) R(s) 1 + G ( s) Xem xét một hệ thống phản hồi với hàm chuyển G(s) của quá trình có dạng như sau: K G ( s) = (5.2) s( s + p) Thay (5.2) vào (5.1), chúng ta có: K C (s) = (5.3) R( s) 2 s + ps + K Viết lại phương trình (5.3) dưới dạng của tỷ số cản ζ và tần số tự nhiên ωn: 2 ωn C (s) = (5.4) R( s) 2 2 s + 2ζω n s + ω n ở đó, ζ = p (2 K ) và ω n = K . Nếu tín hiệu vào r(t) là một hàm nhảy bậc đơn vị, nghĩa là R(s) = 1/s, phương trình (5.4) trở thành: 2 ωn C (s) = (5.5) s ( s 2 + 2ζω n s + ω n ) 2 Đáp ứng theo thời gian của hệ thống có dạng như sau: 1 e −ζω n t sin(ω n 1 − ζ 2 t + θ ) c(t ) = 1 − (5.6) 1− ζ 2 1−ζ 2 ở đó θ = arctan . Đồ thị của c(t) với tín hiệu vào là hàm nhảy bậc đơn vị ζ cho các giá trị khác nhau của tỷ số cản ζ được biểu diễn trong Hình 5.1. Chúng ta có thể thấy được trên đồ thị, khi ζ càng giảm, biên độ dao động của đáp ứng càng tăng. Nếu tín hiệu vào r(t) là một hàm xung đơn vị, nghĩa là R(s) = 1, phương trình (5.4) trở thành: 2 ωn C (s) = (5.7) s 2 + 2ζω n s + ω n 2 68
Đáp ứng theo thời gian của hệ thống có dạng như sau: ωn e −ζω n t sin(ω n 1 − ζ 2 t ) c(t ) = (5.8) 2 1− ζ Đồ thị của c(t) với tín hiệu vào là hàm xung đơn vị cho các giá trị khác nhau của tỷ số cản ζ được biểu diễn trong Hình 5.2. Chúng ta có thể chọn vài số đo hiệu suất từ đáp ứng nhất thời của hệ thống với các tín hiệu vào là tín hiệu nhảy bậc hay tín hiệu xung đơn vị. c(t) ζ = 0,1 ζ = 0,2 ζ = 0,4 ζ = 0,8 ζ = 1,0 ζ = 2,0 ωnt Hình 5.1. Đáp ứng nhất thời của hệ thống với tín hiệu vào là hàm nhảy bậc đơn vị cho các giá trị khác nhau của tỷ số cản ζ Các số đo hiệu suất chuẩn thường được định nghĩa trên cơ sở đáp ứng của hệ thống với tín hiệu vào nhảy bậc (Hình 5.3). Tốc độ của đáp ứng được đo bằng thời gian lên (rise time) Tr và thời gian tới đỉnh (peak time) Tp. Thời gian lên Tr là khoảng thời gian để đáp ứng của hệ thống tăng từ 10% lên 90% của giá trị cuối cùng, còn thời gian tới đỉnh Tp là khoảng thời gian để đáp ứng của hệ thống đạt tới mức cực đại. Một số đo nữa là phần trăm quá mức (percent overshoot) Po được định nghĩa như sau: M p − fv Po = × 100% (5.9) fv ở đó, Mp là giá trị cực đại của đáp ứng và fv là giá trị cuối cùng của đáp ứng. Với tín hiệu vào nhảy bậc, fv thường có giá trị bằng độ lớn của tín hiệu vào. Thời gian quá độ (settling time) Ts được định nghĩa là khoảng thời gian cần thiết để hệ thống ổn định trong một khoảng δ nhất định của giá trị cuối cùng fv. 69
Với hệ thống được biểu diễn bằng phương trình (5.4), giá trị của Ts với δ bằng 2% của fv là bốn lần giá trị của hệ số thời gian τ của hệ thống: 4 Ts = 4τ = (5.10) ζω n c(t) ζ = 0,1 ζ = 0,2 ζ = 0,4 ζ = 0,8 ζ = 1,0 ζ = 2,0 ωnt Hình 5.2. Đáp ứng nhất thời của hệ thống với tín hiệu vào là hàm xung đơn vị cho các giá trị khác nhau của tỷ số cản ζ Để xác định thời gian tới đỉnh Tp, chúng ta cần giải phương trình sau: dc(t ) =0 (5.11) dt Biến đổi Laplace của đạo hàm khi các điều kiện ban đầu bằng không là: L⎡ dc(t ) ⎤ = sC ( s ) (5.12) ⎢ ⎣ dt ⎥ ⎦ Thay (5.5) vào (5.12), chúng ta có được phương trình: 2 ωn L⎡ dc(t ) ⎤ = (5.13) ⎣ dt ⎥ s 2 + 2ζω n s + ω n ⎢ 2 ⎦ Lấy biến đổi Laplace nghịch của phương trình (5.13): ωn dc(t ) e −ζω n t sin(ω n 1 − ζ 2 t ) = (5.14) dt 2 1−ζ Từ (5.11) và (5.14), chúng ta có được phương trình cần giải để xác định Tp là: 70
ωn e −ζω n t sin(ω n 1 − ζ 2 t ) = 0 (5.15) 2 1− ζ Vế trái của phương trình (5.15) bằng không khi ω n 1 − ζ 2 t = kπ với k là một số nguyên. Chúng ta dễ dàng thấy được, khi k = 1, nghiệm của phương trình chính là Tp, nghĩa là: ωn 1 − ζ 2 T p = π (5.16) hay: π Tp = (5.17) 2 ωn 1 − ζ c(t) Mp fv+δ fv fv−δ t Tp Ts Tr Hình 5.3. Các số đo hiệu suất trên đáp ứng của hệ thống với tín hiệu vào nhảy bậc đơn vị Thay (5.17) vào (5.6), chúng ta sẽ tính được giá trị cực đại Mp của đáp ứng: 1 −ζω nT p sin(ω n 1 − ζ 2 T p + θ ) M p = c(T p ) = 1 − e 2 1− ζ (5.18) 1 1 −ζω nT p −ζωnT p sin( π + θ ) = 1 + sin θ = 1− e e 1− ζ 2 1− ζ 2 71
1−ζ 2 , vì vậy sinθ được tính như sau: ở đó θ = arctan ζ ⎛ ⎞ 1−ζ 2 sin θ = sin ⎜ arctan ⎟ = 1−ζ 2 (5.19) ⎜ ⎟ ζ ⎜ ⎟ ⎝ ⎠ Thay (5.19) vào (5.18), chúng ta có được giá trị của Mp theo ζ và ωn: −ζω n T p 1−ζ 2 = 1 + e −ζπ M p = 1+ e (5.20) Đây là giá trị cực đại của đáp ứng với tín hiệu vào là hàm nhảy bậc đơn vị, trong trường hợp này giá trị cuối cùng của đáp ứng là fv = 1. Giá trị phần trăm quá mức Po được tính như sau: 1−ζ 2 1 + e −ζπ −1 1−ζ 2 × 100% = 100e −ζπ Po = (5.21) 1 Khi thiết kế hệ thống, chúng ta thường muốn đáp ứng của hệ thống có cả Po và Tp càng nhỏ càng tốt. Tuy nhiên, trong khi Po giảm theo chiều tăng của tỷ số cản ζ thì Tp lại tăng, như được thể hiện trên đồ thị của ví dụ đang xét (Hình 5.4). Vì vậy, chúng ta cần điều chỉnh lại các yêu cầu ban đầu của Po và Tp để có thể xác định được giá trị phù hợp cho ζ. Po Tp 5,0 4,8 4,6 4,4 Tp 4,2 Po 4,0 3,8 3,6 3,4 3,2 3,0 ζ Hình 5.4. Đồ thị của các số đo phần trăm quá mức Po và thời gian tới đỉnh Tp của đáp ứng hệ thống khi tỷ số cản ζ thay đổi (đặt ωn = 1) Vị trí các nghiệm của phương trình đặc trưng của hệ thống được biểu diễn 72
bằng phương trình (5.4) thể hiện trong Hình 5.5. Đồ thị này cho chúng ta thấy được sự tương quan giữa đáp ứng trong miền thời gian của hệ thống với vị trí các điểm cực của hàm chuyển trong mặt phẳng s thông qua các hệ số ζ, ωn và góc θ. iω iω n 1 − ζ 2 × ωn θ σ -ζωn 0 − iω n 1 − ζ 2 × Hình 5.5. Đồ thị vị trí các nghiệm của phương trình đặc trưng của hệ thống trong mặt phẳng s Đồ thị trong Hình 5.4 biểu diễn các số đo hiệu suất của một hệ thống bậc hai. Tuy nhiên, trong nhiều trường hợp, các số đo hiệu suất của hệ thống có bậc lớn hơn hai cũng có thể xấp xỉ được bằng đồ thị này. Ví dụ, xem xét một hệ thống bậc ba có hàm chuyển vòng kín như sau: 2 ωn T (s) = (5.22) ( s 2 + 2ζω n s + ω n )( s + γ ) 2 Phương trình đặc trưng của hệ thống nói trên có hai nghiệm là các nghiệm của phương trình đặc trưng của hệ thống bậc hai biểu diễn bởi phương trình (5.4) và nghiệm thứ ba bằng −γ. Các số đo hiệu suất của hệ thống này sẽ có giá trị xấp xỉ các số đo hiệu suất của hệ thống bậc hai nếu như: | γ |≥ 10 | ζω n | (5.23) Khi đó cặp nghiệm của hệ thống bậc hai được gọi là các nghiệm trội (dominant roots) của hệ thống bậc ba. Nói một cách khác, đáp ứng của một hệ thống bậc ba có thể xấp xỉ được bằng đáp ứng của một hệ thống bậc hai nếu độ lớn phần thực của các nghiệm trội nhỏ hơn 1/10 độ lớn phần thực của nghiệm thứ ba của phương trình đặc trưng. Ví dụ 5.1 Xem xét hệ thống được biểu diễn bằng phương trình (5.3). Chúng ta muốn chọn giá trị các tham số K và p sao cho giá trị phần trăm quá mức Po của đáp ứng của hệ thống đối với tín hiệu vào nhảy bậc không vượt quá 5% và thời gian quá độ Ts không quá 4s. Ts được tính theo biểu thức (5.10), do vậy chúng ta có: 4 ≤ 4 hay ζω n ≥ 1 Ts = (5.24) ζω n 73
Sử dụng công thức (5.21), chúng ta tính được điều kiện để Po ≤ 5% là ζ ≥ 1 2 hay θ ≤ 45o. Để thỏa mãn các yêu cầu, đồng thời có được giá trị của Ts và Tp càng nhỏ càng tốt, chọn giá trị nhỏ nhất có thể cho ζ là 1 2 và ω n = 1 ζ = 2 . Khi 2 đó, chúng ta có được các giá trị K = ω n = 2 và p = 2ζωn = 2. 5.3. Chỉ số hiệu suất Chỉ số hiệu suất (performance index) là một số đo định lượng hiệu suất của một hệ thống và được lựa chọn sao cho phù hợp với các yêu cầu đặt ra cho hệ thống. Sau khi hiệu suất được mong muốn cho hệ thống đã được mô tả một cách định lượng, chỉ số hiệu suất sẽ được tính toán hay đo đạc và được dùng để đánh giá hiệu suất của hệ thống. Trong kỹ thuật điều khiển hiện đại, chỉ số hiệu suất là một khái niệm rất quan trọng và cần thiết cho các hệ thống điều khiển thích nghi, cho việc tối ưu hóa các tham số hệ thống một cách tự động và thiết kế các hệ thống điều khiển tối ưu. Một hệ thống được coi là hệ thống điều khiển tối ưu (optimum control system) khi các tham số của hệ thống được điều chỉnh được điều chỉnh sao cho các chỉ số đạt được cực trị, thường là giá trị cực tiểu. Các chỉ số hiệu suất thường được sử dụng bao gồm: 1. Tích phân của sai số bình phương (ISE): T ∫ I1 = e 2 (t )dt (5.25) 0 Giới hạn trên T là một giá trị thời gian hữu hạn được chọn sao cho tích phân tiếp cận giá trị ở trạng thái thường trực. Người ta thường chọn T = Ts, tức là thời gian quá độ của hệ thống. Chỉ số này thường được dùng để phát hiện ra các hệ thống có tỷ số cản quá lớn hay quá nhỏ so với mức phù hợp. 2. Tích phân của sai số tuyệt đối (IAE): T ∫ I 2 = | e(t ) | dt (5.26) 0 Chỉ số này thường được dùng khi chúng ta quan tâm nhiều hơn tới sai số ở đoạn cuối của đáp ứng nhất thời. 3. Tích phân của thời gian nhân sai số tuyệt đối (ITAE): T ∫ I 3 = t | e(t ) | dt (5.27) 0 4. Tích phân của thời gian nhân sai số bình phương (ITSE): 74
T ∫ I 4 = te 2 (t )dt (5.28) 0 Các chỉ số I3 và I4 được dùng vì chúng khuyếch đại sự thay đổi của sai số khi các tham số của hệ thống thay đổi, trong đó chỉ số I3 phân biệt rõ nhất, cũng có nghĩa là điểm cực tiểu của nó được thể hiện rõ ràng nhất. Ngoài các chỉ số hiệu suất nêu trên, chúng ta có thể định nghĩa các chỉ số khác dưới dạng tổng quát: T ∫ f [e(t ), r (t ), c(t ), t ]dt I= (5.29) 0 ở đó f là một hàm của sai số, tín hiệu vào, tín hiệu ra và thời gian. Ví dụ 5.2 Xem xét một hệ thống điều khiển phản hồi có hàm chuyển của toàn hệ thống như sau: 1 T (s) = (5.30) 2 s + 2ζs + 1 Các đồ thị trong Hình 5.6 thể hiện sự biến thiên của các chỉ số hiệu suất của hệ thống khi tỷ số cản ζ biến đổi, với tín hiệu vào là một hàm nhảy bậc đơn vị. Nhìn vào hình vẽ chúng ta thấy được chỉ số I3 có sự thay đổi rõ ràng nhất, và chỉ số này đạt được giá trị cực tiểu khi tỷ số cản ζ có giá trị khoảng 0,7. ζ Hình 5.6. Đồ thị các chỉ số hiệu suất khi tỷ số cản ζ thay đổi 75
5.4. Sai số ở trạng thái thường trực của hệ thống điều khiển phản hồi Một trong những nguyên nhân chủ yếu của việc sử dụng phản hồi, cho dù làm tăng giá thành và độ phức tạp của hệ thống, là khả năng làm giảm sai số ở trạng thái thường trực của hệ thống. Như đã được đề cập ở chương trước, sai số ở trạng thái thường trực của một hệ thống vòng kín thường nhỏ hơn vài lần so với sai số của hệ thống vòng hở. Chúng ta đã dùng ký hiệu Ea(s) để chỉ tín hiệu sai khác được dùng để điều khiển quá trình trong hệ thống vòng kín. Tuy nhiên, sai số thực sự của hệ thống phải là E(s) = R(s) − C(s). Với hệ thống phản hồi như trong Hình 4.1, chúng ta có: 1 + G (s) H ( s) − G ( s) G ( s) E (s) = R(s) − R(s) = (5.31) R(s) 1 + G (s) H (s) 1 + G (s) H (s) Nếu H(s) = 1: 1 E ( s) = (5.32) R( s) 1 + G (s) Khi đó, E(s) = Ea(s). Sai số ở trạng thái thường trực của hệ thống khi H(s) = 1 là: sR ( s ) ess = lim e(t ) = lim (5.33) s →0 1 + G ( s ) t →∞ Sai số ở trạng thái thường trực của hệ thống cho một số tín hiệu vào thử chuẩn khi H(s) = 1 được xác định như sau: − Tín hiệu nhảy bậc: s( A s) A ess = lim = (5.34) s →0 1 + G ( s ) 1 + G ( 0) ở đó A là độ lớn của tín hiệu nhảy bậc. Giả sử hàm chuyển G(s) của quá trình có dạng như sau: M ∏ (s − z ) K i i =1 G ( s) = (5.35) Q ∏ (s − p ) sN j j =1 Giá trị N chính là số lần tích phân, hay còn gọi là số định kiểu (type number) của hệ thống, bởi vì sai số ở trạng thái thường trực phụ thuộc vào giá trị này, cụ thể như sau: o Hệ thống kiểu-0 (type-zero): A A ess = = (5.36) 1+ K p ⎛ ⎞ ⎛ ⎞ Q M ∏ ∏ ⎜ − pj ⎟ 1+ K ⎜ − zi ⎟ ⎜ ⎟ ⎜ ⎟ ⎝ ⎠ ⎝ ⎠ i =1 j =1 76