- Trang Chủ
- Tự động hoá
- Tổng hợp hệ thống điều khiển phi tuyến bậc cao có tham số biến đổi dùng điều khiển trượt bậc hai
Xem mẫu
- Nghiên cứu khoa học công nghệ
TỔNG HỢP HỆ THỐNG ĐIỀU KHIỂN PHI TUYẾN BẬC CAO CÓ
THAM SỐ BIẾN ĐỔI DÙNG ĐIỀU KHIỂN TRƯỢT BẬC HAI
Lê Trần Thắng1, Trần Văn Nhân2*, Nguyễn Khắc Tuấn3
Tóm tắt: Điều khiển trượt – Sliding mode control (SMC) là một phương pháp
điều khiển được áp dụng rộng rãi với những ưu điểm nổi bật là sự bền vững chống
lại nhiễu loạn và sự biến đổi của tham số, thiết kế hệ thống giảm bậc và cấu trúc
điều khiển đơn giản. Tuy nhiên, nếu bậc của hệ thống lớn hơn 2 thì phải xử lý vấn
đề siêu bề mặt trượt, điều này đòi hỏi những tính toán phức tạp. Bài báo này trình
bày một phương pháp tổng hợp một lớp hệ thống phi tuyến bậc cao trên cơ sở sử
dụng SMC bậc hai nhằm đơn giản hóa tính toán và tận dụng các ưu điểm của trượt
bậc hai.
Từ khóa: Điều khiển trượt bậc hai, Giảm chattering, Hàm chuyển mạch.
1. ĐẶT VẤN ĐỀ
Trong thực tế, nhiều giải pháp tổng hợp bộ điều khiển theo các phương pháp kinh điển
gặp khó khăn với những hệ thống phi tuyến, hệ thống có các tham số thay đổi... Nếu
không kiểm soát được sự biến động của các tham số, không loại trừ được nhiễu thì chất
lượng động học hệ thống sẽ bị thay đổi, hệ thống sẽ mất ổn định. Có nhiều phương pháp
để giải quyết những khó khăn này một trong số đó là sử dụng bộ điều khiển trượt (SMC)
[3]. Đây là phương pháp điều khiển hiện đại có thế mạnh trong tổng hợp các hệ thống điều
khiển phức tạp. Tuy nhiên, hạn chế cơ bản của SMC là hiện tượng “chattering” [6]. Hiện
nay, SMC đang là một hướng được quan tâm nghiên cứu, phát triển và áp dụng cho nhiều
hệ thống kỹ thuật, nhiều đối tượng, trong đó có các hệ thống tự động bám sát mục tiêu
(HTĐBSMT) [3, 4].
Khi thiết kế bộ điều khiển SMC cho các hệ thống bậc cao thì phương trình mặt trượt có
bậc lớn, hay gọi là siêu mặt trượt, điều đó tăng mức độ phức tạp và khó khăn trong tính
toán. Để giảm đi khó khăn trên bài báo này sẽ đề xuất một phương pháp tổng hợp SMC
cho các hệ thống điều khiển phi tuyến bậc cao có tham số biến đổi trên cơ sở bộ điều khiển
trượt bậc hai.
2. MÔ TẢ VÀ TỔNG HỢP BỘ ĐIỀU KHIỂN CHO HỆ THỐNG
2.1. Mô tả hệ thống
Xét một lớp hệ thống điều khiển phi tuyến SISO bậc n có tham số biến đổi với nhiễu
loạn, các yếu tố không chắc chắn được mô tả như sau:
i =n
(i )
a x
i =1
i (t ) + f ( x) + d (t ) = u (t ) (1)
Trong phương trình (1): x(t) là tín hiệu ra; x(i)(t) là vi phân cấp i của x(t); ai là hệ số
tương ứng của x(i)(t); u(t) tín hiệu điều khiển đầu vào; d(t) là thành phần nhiễu bên ngoài
và phụ thuộc vào thời gian; f(x) là nhiễu phụ thuộc trạng thái; ai là các hệ số thỏa mãn
điều kiện Huzwit, lớn hơn không.
Trong hệ thống (1) không mất tính tổng quát nếu chúng ta giả thiết các thành phần ai,
d(t) và f ( ) bị chặn theo qui luật sau:
ai aiM ; d d M và f f M
Tạp chí Nghiên cứu KH&CN quân sự, Số 49, 06 - 2017 35
- Kỹ thuật điều khiển & Điện tử
Trong điều khiển nếu bậc của hệ thống lớn hơn hai thì phải giải quyết vấn đề siêu bề
mặt trượt trong không gian nhiều chiều, điều này đòi hỏi những kỹ thuật tính toán phức
tạp. Trong bài báo này đưa ra một phương pháp sử dụng bộ điều khiển trượt bậc hai cho hệ
thống (1) được sử dụng bằng cách xem các thành phần bậc ba trở lên là nhiễu bậc cao.
Trong [9] đã bỏ qua các thành phần bậc cao nếu hệ số nhỏ hơn nhiều so với thành phần
bậc hai mà không ảnh hưởng nhiều đến các chỉ tiêu đặc trưng của quá trình điều khiển. Vì
i=n
(i )
thế, trong hệ thống (1) nếu các thành phần a x
i =3
i (t ) bị chặn khi t tiến đến ∞ và các hệ
số ai đủ nhỏ thì (1) được xem như hệ thống bậc hai được biểu diễn như sau:
a2
x(t ) + a1 x (t ) + F ( x) + d (t ) = u (t ) , (2)
trong đó:
i =n
F ( x) = ai x (i ) (t ) + f ( x) , F ( x) FM
i =3
là thành phần nhiễu phụ thuộc trạng thái.
Việc nghiên cứu hệ thống (1) sẽ được thay thế bởi hệ thống (2). Sơ đồ điều khiển (2)
được mô tả như sau:
Hình 1. Sơ đồ cấu trúc hệ thống (2).
Do ai biến đổi theo thời gian và hệ thống có các yếu tố nhiễu loạn, không chắc chắn
nên chúng ta dùng mô hình danh định sau làm mô hình tham chiếu cho hệ thống (2):
a2 n
xn (t ) + a1n xn (t ) = (t ) (3)
Trong đó: a1n , a2n – Các hệ số danh định, thỏa mãn điều kiện Huzwit, lớn hơn không;
xn(t) – Tín hiệu ra danh định; xd(t) – Tín hiệu đặt; (t ) – Tín hiệu điều khiển danh định.
Để điều khiển hệ thống (3) sử dụng bộ điều khiển sau [3]:
a1 n
= a 2 n (
xd + x d h1e h2 e )
a2 n
a1n
với các hệ số h1 = k 2 và h2 = 2k , k 0 ; e = xn xd .
a2 n
2.2. Thiết kế bộ điều khiển cho hệ thống
Với hệ thống (2) trong [2] các tác giả đã đề xuất một bộ điều khiển trượt và một hàm
chuyển mạch mới trong vùng gần bề mặt trượt cho một hệ thống bậc hai có tham số biến
đổi. Bộ điều khiển này đảm bảo tổng hợp hệ thống có các chỉ tiêu chất lượng cao như thời
36 L. T. Thắng, T. V. Nhân, N. K. Tuấn, “Tổng hợp hệ thống điều khiển… trượt bật hai.”
- Nghiên cứu khoa học công nghệ
gian quá độ và sai số xác lập nhỏ, khả năng chống nhiễu tốt. Tuy nhiên, bài báo [2] chưa
chứng minh được tính ổn định tiệm cận của hệ thống khi sử dụng hàm chuyển mạch mới
trong vùng gần bề mặt trượt. Để giải quyết vấn đề đó cần tổng hợp một luật điều khiển cho
hệ thống (2) và chứng minh tính ổn định tiệm cận của hệ thống (2) trong vùng gần bề mặt
trượt theo tiêu chuẩn ổn định Lyapunov.
Ta có bổ đề sau.
Bổ đề 1: Xét hệ thống điều khiển (2) với p là số nguyên không âm; , , , K , h là các
số thực dương và nếu hệ (2) thỏa mãn các điều kiện sau:
1
d M + FM 2 p +1
h 0 khi s
K
1
d + F
h = d + F khi s M M 2 p +1
M M K
h
d + F sin( s ) ; khi s
trong đó, s là mặt trượt, thì bộ điều khiển:
(2 p +1) hsign ( s ) + a ( l x ) + a x khi s
Ks 2
a2 n 1
u= (4)
Ks (2 p +1) h sin( s ) + a ( l x ) + a x khi s
2 1
a2 n
đảm bảo cho hệ thống (2) ổn định tiệm cận.
Sau đây, ta sẽ chứng minh bổ đề 1 cho từng trường hợp.
2.2.1. Chứng minh trường hợp s
Để chứng minh bổ đề trên đặt en = x xn , lựa chọn mặt trượt như sau: s = en + l en
a1n
trong đó l = 0 . Vậy s = en + l en .
a2 n
1
a2 s . Đạo hàm của V là: V = a2 ss , hay :
2
Chọn hàm Lyapunov dạng: V =
2
V = a2 s (en + l en ) = a2 s[(
x
xn ) + l ( x xn )]
a a a
= s[( a
x + a1 x ) a1 x 2
a2 n
xn 2
a1n xn + 2
a1n xn + a l ( x xn )] .
2 a2 n a2 n a2 n 2
a1n
Thay l= vào biểu thức trên, ta có:
a2 n
a
V = s[( a2 x + a1 x ) 2 ( a2 n xn + a1n xn ) a1 x + a2l x ]
a2 n
Tạp chí Nghiên cứu KH&CN quân sự, Số 49, 06 - 2017 37
- Kỹ thuật điều khiển & Điện tử
a2
= s[u d F a1 x + a2l x ] . (5)
a2 n
Kết hợp (4) (trường hợp s ) và (5) được:
V = a2 ss = Ks 2( p +1) h s s (d + F ) Ks 2( p +1) s (h d F ) , vì a2 0
suy ra:
1
ss
( Ks 2( p +1) + s (h d F ) . (6)
a2
Theo lý thuyết hàm Lyapunov để hệ thống (2) ổn định tiệm cận thì đạo hàm
V = a2 ss 0 , tương ứng với vế phải của (6) phải nhỏ hơn không, điều đó đạt được khi bất
đẳng thức sau thỏa mãn:
(2 p +1)
s (K s ( d + F )) + s h 0 . (7)
Do d d M và F FM nên chỉ cần chứng minh:
(2 p +1)
s (K s (d M + FM )) + s h 0 (8)
Từ (7) cho thấy:
1
d + FM 2 p +1
- Nếu s M thì (7) luôn đúng nếu h>0 nhỏ tùy ý.
K
1
d + FM 2 p +1
- Nếu s M thì (8) được viết lại như sau:
K
(2 p +1)
s [( K s + h (d M + FM )] 0 ,
Khi đó do h = d M + FM (9)
Nên (7) được thỏa mãn, hay đạo hàm của hàm Lyapunov V = a2 ss 0 . Điều kiện
trượt được thỏa mãn chặt và hệ thống (2) ổn định tiệm cận. (đpcm)
(2 p +1)
Trong (4), do p là số nguyên không âm nên thành phần Ks sẽ nhanh chóng đưa
hệ thống (1) tiến về mặt trượt khi s lớn và cũng nhanh chóng giảm về 0 khi hệ thống đã ở
gần mặt trượt (s gần bằng không) mà không làm tăng “chattering” trong hệ thống do đây là
một hàm liên tục.
1
d + FM 2 p +1
Nhận xét: Khi s ở xa mặt trượt, s M thì chọn h>0 nhỏ tùy ý, khi s ở
K
1
d + FM 2 p +1
gần bề mặt trượt hơn s M thì phải chọn h thỏa mãn (9).
K
Vế thứ nhất của bộ điều khiển (4) đảm bảo cho hệ thống (2) ổn định tiệm cận [2], lúc
này thành phần không liên tục trong bộ điều khiển: hsign( s ) chỉ phụ thuộc vào biên độ
nhiễu, loại bỏ được sự ảnh hưởng về thay đổi các tham số của hệ thống. Điều này cho
phép làm giảm “chattering” trong hệ thống. Tuy nhiên, hsign( s ) vẫn giao động với tần số
38 L. T. Thắng, T. V. Nhân, N. K. Tuấn, “Tổng hợp hệ thống điều khiển… trượt bật hai.”
- Nghiên cứu khoa học công nghệ
cao trong biên độ h, h . Điều này gây ra “chattering” ảnh hưởng xấu đến các hệ thống
cơ khí và có thể phá hủy phần dẫn động.
2.2.2. Chứng minh trường hợp s
Để giảm ảnh hưởng của “chattering”, trong vùng gần bề mặt trượt khi s tác giả sử dụng
hàm liên tục sin( s ) / [2]. Trong [2] đã hạn định 1, mặt khác ta chỉ cần xét trong
trường hợp s / 2 tức là / 2 vậy:
1 / 2 (10)
Kết hợp (4) (trường hợp s ) và (5) được:
2( p +1) sin( s )
V = Ks hs + s (d + F )
Hệ thống sẽ ổn định theo tiêu chuẩn Lyapunov nếu:
2( p +1) sin( s )
Ks hs + s (d + F )0 (11)
sin( s ) 2( p +1)
Hay hs + s ( d + F ) 0 (do Ks 0 và p > 0)
sin( s )
s h ( d + F ) 0
sin( s )
h ( d + F ) 0 (12)
Do giả thiết d + F sin( s ) (với là hằng số nào đó) khi s , khi đó chọn
h h
hay thì (12) đúng nên (11) đúng. (đpcm)
Thông thường, hệ số được chọn đủ lớn để sin( s ) giảm về 0 từ giá trị biên ,
giá trị này gấp lần giá trị biên chuyển mạch nhằm hạn định giá trị không quá lớn.
3. KẾT QUẢ MÔ PHỎNG
Mô phỏng HTĐBSMT sử dụng động cơ một chiều không chổi than (Torque BLDC
motor) có mô hình toán học đầy đủ được mô tả như sau [1],[6]:
K e [TeTm
x (t ) + Tm
x (t ) + x (t )] + f ( x ) + d (t ) = u (t )
Trong đó, Te là hằng số điện từ, Tm là hằng số điện cơ, Ke là hệ số tỷ lệ.
Vậy: a3=KeTeTm; a2=KeTm; a1=Ke;
Đây là trường hợp riêng của (1) khi n=3. Biểu diễn dưới dạng (2) được:
K eTm x + K e x + F ( x ) + d (t ) = u (t )
F ( x ) = f ( x ) + K eTeTm
x
Thực hiện mô phỏng với động cơ điện S Series Brushless Torque Motor- Aerotech S-
180-44-A có công suất 313W, các tham số được quy đổi theo [6], [7]. Tham số momen
Tạp chí Nghiên cứu KH&CN quân sự, Số 49, 06 - 2017 39
- Kỹ thuật điều khiển & Điện tử
quán tính của hệ thống được tính trên cơ sở tải cân bằng và có khối lượng 30 kg. Các tham
số của đối tượng:
1 2 3
K t = 2, 2 Nm / A; J = 7, 40 10 kgm ; L = 3, 4 10 H
1 3
R = 12,8 ; K e = 2,5650897; Tm = 17 10 ( s ); Te = L R = 0, 27 10 (s)
3
K eTm = 4,3606525; K eTeTm = 1,177376 10
Giả sử các hệ số ai biến thiên theo quy luật sau:
a1 = 2,5650897 + sin(t ); a2 = 4,3606525 + 2sin(2t );
a3 = 1,177376 103 + 103 sin(3t ).
Vậy a1n = 2,5650897; a2 n = 4,3606525 .
f x = 10 [ sin( x ) + 2 sin 2 x + sin x + cos x ]
Đặt :
d t = 20 [ sin(t ) + sin 2t + 3sin 3t + cos t ]
Vậy: F ( x ) = f ( x ) + 1,177376 10 3
x
Tín hiệu đặt: xd t = sin( t ). Chọn = 0.02 . Trong (4), nếu chọn K > 0 quá lớn
sẽ làm biên độ điện áp đầu ra bộ điều khiển quá lớn, có thế vượt qua các giới hạn vật lý
của hệ thống. Trong trường hợp này, ví dụ chọn K = 10.
Giá trị lớn nhất của nhiễu: f M = 50; d M = 120; FM = 50 + 1,177376 103.
Đối với bộ điều khiển của hệ thống (3), do điều kiện k > 0 và k càng lớn thì thời gian
phản ứng của hệ thống càng nhanh, ví dụ chọn k = 8, tính được: h1=k2=64 và
a1n 2,5650897
h2 = 2k = 2* 8 = 15.412 .
a2 n 4,3606525
Hình 2. Sơ đồ mô phỏng hệ thống (2).
1 - Bộ điều khiển cho hệ thống danh định; 2 - Mô hình danh định của hệ thống (2);
3 - Bộ điều khiển đề xuất (4); 4 - Mô hình thực tế của hệ thống (2).
Mô phỏng và so sánh trong trong các trường hợp sau :
- Trường hợp 1: Mô phỏng hệ thống (2) trong trường hợp bộ điều khiển (4) chỉ sử dụng
40 L. T. Thắng, T. V. Nhân, N. K. Tuấn, “Tổng hợp hệ thống điều khiển… trượt bật hai.”
- Nghiên cứu khoa học công nghệ
hàm dấu chuyển mạch sign(s).
- Trường hợp 2: Mô phỏng hệ thống (2) trong trường hợp sử dụng bộ điều khiển (4)
với các giá trị thỏa mãn (10). Kết quả cho thấy sai số xác lập giảm khi tăng. Lấy ba
đại lượng sau làm ví dụ mô phỏng: α = 2; α = 20; α = π/2δ = 78.25975.
Hình 3. Sai số bám khi sử dụng (4) với H=FM+dM ;sign(s).
Hình 4. Sai số bám khi sử dụng (4) với α = 2.
Hình 5. Sai số bám khi sử dụng (4) với α = 20.
Tạp chí Nghiên cứu KH&CN quân sự, Số 49, 06 - 2017 41
- Kỹ thuật điều khiển & Điện tử
Hình 6. Sai số bám khi sử dụng (4) với α = 78.25975.
Bảng 1. Bảng so sánh các kết quả mô phỏng và tính toán.
Trường Bộ điều Giá trị h, tại Thời gian quá độ Sai số xác
hợp khiển = 0.02 ( tqd ) lập ( )
1 (4) H=FM+dM ;sign(s) 0.5669 (s) 1.3 102
2 (4) h=FM+dM ; = 2 0.5539 (s) 102
3 (4) h=FM+dM ; = 20 0.5894 (s) 5 103
h=FM+dM
4 (4) 0.5926 (s) 2 104
= 78.25975
Nhận xét : Khi mô phỏng với bộ điều khiển (4) trong các trường hợp khác nhau thì
thời gian quá độ chênh lệch không đáng kể (nhỏ hơn 0.03 s). Trong các trường hợp, thời
gian quá độ tqd 0.6(s) . Sai số xác lập khi sử dụng bộ điều khiển (4) với hàm chuyển
mạch sign(s) lớn hơn đáng kể khi sử dụng bộ điều khiển (4) với hàm sign(s) kết hợp hàm
đề xuất sin( s ) / . Sai số xác lập giảm nhanh khi tăng, đặc biệt khi = max = / 2
thì sai số xác lập 2 104.
4. KẾT LUẬN
Bài báo đã đề xuất một phương pháp điều khiển hệ phi tuyến bậc cao sử dụng bộ điều
khiển trượt bậc hai và tính toán hạn định một tham số trong hàm chuyển mạch đã đề xuất
trong vùng gần bề mặt trượt. Tính ổn định theo tiêu chuẩn Lyapunov đã được chứng minh
trong vùng xa và gần bề mặt trượt. Hệ thống tương đương (2) đạt được các chỉ tiêu chất
lượng cao, có khả năng chống nhiễu tốt, không bị ảnh hưởng bởi sự thay đổi của các tham
số và giảm thiểu được chattering trong quá trình xác lập.
TÀI LIỆU THAM KHẢO
[1]. Nguyễn Trung Kiên, “Xây dựng phương pháp tổng hợp hệ thống điều khiển các đài
quan sát tự động định vị từ xa các đối tượng di động”, Luận án tiến sỹ kỹ thuật, Viện
Khoa học và Công nghệ quân sự, 2015.
[2]. Trần Văn Nhân, Lê Trần Thắng, Nguyễn Thượng San, “Một phương pháp tổng hợp
hệ thống điều khiển phi tuyến có tham số biến đổi trong vùng gần bề mặt trượt” , Tạp
chí Nghiên cứu KH và CN Quân sự, số đặc san Tên lửa, tháng 9/2016, tr. 77-83.
42 L. T. Thắng, T. V. Nhân, N. K. Tuấn, “Tổng hợp hệ thống điều khiển… trượt bật hai.”
- Nghiên cứu khoa học công nghệ
[3]. Jinkun Liu, Xinhua Wang, “Advanced Sliding Mode Control for Mechanical
Systems”, Tsinghua University Press, Beijing, 2012.
[4]. Xinghuo Yu, Okyay Kaynak,“Sliding-Mode Control With Soft Computing: A Survey”,
IEEE Transactions on Industrial Electronics, Vol. 56, N0. 9, September 2009.
[5]. Lee, H, & Utkin, V. I. “The Chattering Analysis”, Advances in Variable Structure and
Sliding Mode Control, Lecture Notes in Control and Information Sciences (Vol. 334),
2006.
[6]. PITTMAN Motors, “Servo motor application note”.
[7].http://www.aerotech.com/product-catalog/motors/rotary-motors/bm-
series.aspx?p=%2fproduct-catalog%2fmotors.aspx
[8].Смольников П.П, “Синтез квазиоптимальных систем aвтоматического
управления”. Энергия, Ленинград 1967.
ABSTRACT
A SYNTHESIS METHOD OF HIGHER-ORDER NONLINEAR CONTROL SYSTEM
WITH VARIABLE PARAMETERS USING 2ND SLIDING MODE CONTROL.
Sliding mode control (SMC) is a widely used control method with outstanding
advantages of durability against disturbance and variable parameters, reduced-
order system design and simple control structure. However, if the system order is
higher than 2, then the problem of hyperplane which requires complex calculations
has to be treated. In this article, a second-order sliding mode control (2nd-SMC) is
proposed for synthesizing a higher-order nonlinear system class based on the use of
a second-order SMC to simplify calculations and take advantage of 2nd-SMC.
Keywords: Second-order sliding mode control, Chattering reduction, Switching function.
Nhận bài ngày 11 tháng 5 năm 2017
Hoàn thiện ngày 01 tháng 6 năm 2017
Chấp nhận đăng ngày 20 tháng 6 năm 2017
Địa chỉ: 1 Viện Tự động hóa Kỹ thuật quân sự;
2
Học viện Hải quân;
3
Đại học Sư phạm Vinh.
*
Email: tranvannhannt@gmail.com.
Tạp chí Nghiên cứu KH&CN quân sự, Số 49, 06 - 2017 43
nguon tai.lieu . vn
