- Trang Chủ
- Tự động hoá
- Ổn định máy bay cánh cụp cánh xòe trên cơ sở điều khiển Backstepping và hiệu chỉnh RHO
Xem mẫu
- KHOA HỌC CÔNG NGHỆ P-ISSN 1859-3585 E-ISSN 2615-9619
ỔN ĐỊNH MÁY BAY CÁNH CỤP CÁNH XÒE TRÊN CƠ SỞ
ĐIỀU KHIỂN BACKSTEPPING VÀ HIỆU CHỈNH RHO
THE COMBINATION OF THE BACKSTEPPING METHOD AND OPTIMAL FILTER TO THE FLIGHT STABILITY
OF MORPHING AIRCRAFT
Nguyễn Ngọc Tuấn1, Trần Xuân Tình2,*,
Nguyễn Tuấn Anh2, Trần Thị Nga2
cản gió, bay cho nhanh. Cấu hình cánh cụp cánh xòe hữu
TÓM TẮT
dụng nhất trong trường hợp máy bay được yêu cầu phải
Để đảm bảo độ ổn định bay của máy bay cánh cụp, cánh xòe, bài báo đề xuất hoạt động tốt ở tốc độ cao lẫn tốc độ thấp, và vì vậy nó chủ
một thiết kế điều khiển - hiệu chỉnh kết hợp. Phương pháp điều khiển trượt yếu được sử dụng trong việc thiết kế máy bay quân sự, như
Backstepping được sử dụng bảo đảm tính ổn định quỹ đạo bay. Hiệu chỉnh góc mẫu máy bay Mirage G8 của Không quân Pháp ở hình 1. Ví
đường bay theo phương pháp tối ưu hoá từng đoạn. Từ kết quả mô phỏng cho dụ với máy bay chiến đấu chủ yếu bay với tốc độ hạ âm,
thấy góc đường bay, quỹ đạo bay đã bám sát tín hiệu đặt trong quá trình biến đổi nhưng bất thình lình tăng tốc độ lên siêu âm khi cần thiết.
cấu hình cánh, hệ thống điều khiển bay đáp ứng tốt các yêu cầu về tính bền vững
và tính tác động nhanh. Tuy nhiên khi cấu hình cánh biến đổi khi bay sẽ ảnh
hưởng xấu đến đặc tính động học của máy bay, dẫn đến
Từ khóa: Điều khiển trượt cuốn chiếu, tối ưu hóa từng đoạn, máy bay, góc những thay đổi phi tuyến của lực khí động và mô-men
đường bay
xoắn, làm tăng độ bất định của mô hình, sự ổn định bay sẽ
ABSTRACT bị ảnh hưởng rất nhiều. Để đáp ứng yêu cầu về chất lượng
bay, hệ thống điều khiển bay được thiết kế cần đảm bảo sự
To ensure the flight stability of morphing aircraft, the paper proposed a
ổn định bay trong quá trình này.
correction system combined with control system. The nominal controller was
designed by the backstepping technique, which was used to provide the basic
flight stability as well as the tracking performance. The correction device was
designed by the receding horizon optimal algorithm based on the command
filter. The retrofit value was calculated to compensate for the nominal controller.
The simulation results show that the flight path angle can track the command
signal regardless of morphing process, and the flight control system satisfies the
requirements of real-time and robustness.
Keywords: Backstepping controller, Receding Horizon Optimal, aircraft,
altitude, flight path angle.
1
Học viện Kỹ thuật Quân sự
2
Học viện Phòng không - Không quân
*
Email: tinhpk79@gmail.com
Ngày nhận bài: 18/10/2021
Ngày nhận bài sửa sau phản biện: 10/12/2021 Hình 1. Máy bay cánh cụp cánh xòe
Ngày chấp nhận đăng: 27/12/2021 Phương pháp H∞ đã được sử dụng thiết kế bộ điều
khiển cho UAV biến đổi cấu hình cánh [1], tuy nhiên độ ổn
định quỹ đạo bay giảm khi tốc độ biến đổi tăng lên. Bộ điều
1. ĐẶT VẤN ĐỀ khiển Backstepping trong [2-4] đã giảm thiểu sai số bám sát
Máy bay cánh cụp cánh xòe là dạng máy bay có hai quỹ đạo, tăng khả năng chống nhiễu. Mạng neural
cánh chính có thể xoay để thay đổi góc giữa cánh và thân, networks đã được sử dụng để ước lượng tham số mô hình
tức có thể "xòe ra" hay "cụp vào" được. Điều này giúp cấu khí động học của máy bay khi biến đổi cấu hình [3, 4], với
hình cánh của máy bay có thể được thay đổi trong quá phương pháp này cần bộ dữ liệu lớn để huấn luyện mạng.
trình bay. Khi bay chậm, để tránh bị tròng trành ở tốc độ Phương pháp điều khiển trượt [5] được đề xuất ổn định góc
thấp, máy bay xòe hai cánh dang rộng ra để giữ thăng chúc ngóc, kết quả mô phỏng cho thấy hệ bền vững và có
bằng. Còn khi bay nhanh, máy bay cụp cánh ra sau để giảm thể đảm bảo sự ổn định bay trong quá trình gấp cánh.
30 Tạp chí KHOA HỌC VÀ CÔNG NGHỆ ● Tập 57 - Số 6 (12/2021) Website: https://jst-haui.vn
- P-ISSN 1859-3585 E-ISSN 2615-9619 SCIENCE - TECHNOLOGY
Trong bài báo này đề xuất bộ điều khiển trượt
θ f1(x) g1α d1(t)
Backstepping kết hợp tối ưu hoá từng đoạn RHO (Receding
α f2 (x) q d2 (t)
Horizon Optimal) để tái thiết cấu trúc bộ điều khiển, để (5)
điều khiển bám sát góc đường bay của máy bay khi biến
q f3 (x) g2 δ e d3 (t)
đổi cấu hình cánh. q
2. MÔ HÌNH ĐỘNG HỌC
Trong đó f1(x), f2(x), f3(x), g1 và g2 là các thông số hệ
Khoang tĩnh S của máy bay có thể được xấp xỉ bằng thống được biểu diễn như sau:
công thức sau [6]:
1
3
f1(x) QSw(CL0 CLδe δe ) Psinα-mgcosθ-FIz
S mkrk 2mr m3r3x i Oj Ok (1)
mV
1 1x
k1
1
f1(x) (-L-Psinα mgcosθ FIz )
Khoang tĩnh máy bay dọc theo trục x của thân là
mV
S x 2m1r1x m3r3x . Hệ thống điều khiển góc đường bay qcA
1 QSwcA (Cm0 Cmαα Cmq )-Sxgcos
theo chiều dọc của máy bay biến hình là:
f3(x) 2V
Iy (6)
IyqPZT MIy
θ 1 (LPsinα mgcosθFIz )
mV
1
g1 QSwCLα
α 1 (LPsinα mgcosθ F )
mV
mV Iz
1
(2)
g2 QSwcACmδe
Iy
q Iy q 1 (S gcos M PZ M )
Iy Iy
x A T Iy
x θ α q là vector trạng thái, d1(t), d2(t), d3(t)
T
q
là những nhiễu loạn bên ngoài. Vì giá trị chính xác của các
Trong đó, P: Lực đẩy của động cơ, có phương, chiều tham số hệ thống và các thông số khí động học không
trùng với trục dọc của máy bay; ϑ: Góc chúc ngóc; α: Góc được biết trong quá trình bay thực tế, phương trình (5) có
tấn; θ: Góc nghiêng quỹ đạo; m: Khối lượng máy bay; g: Gia thể được biểu diễn dưới dạng phi tuyến tính như sau:
tốc trọng trường; V: Tốc độ bay, để nghiên cứu ảnh hưởng
θ f10 (x) g10 α 1(x)
của quá trình biến hình đối với góc đường bay, trong [6]
chọn điều khiển V ở một giá trị cố định; Iy: Quán tính
α f20 (x) q 2 (x)
(7)
chuyển động xung quanh trục y của máy bay; T: Lực đẩy
q f30 (x) g20 δ e 3 (x)
song song với trục x của máy bay; ZT: Vị trí đặt lực đẩy T; FIz
q
và MIy: Lực quán tính và moment quán tính gây ra bởi quá
trình biến đổi cấu hình cánh: Trong đó, f10(x), f20(x), f30(x), g10 và g20 là các tham số hệ
F S (qcosα q2 sinα) 2S xqcosα Sx sinα thống danh định đã biết; ∆1(x), ∆2(x), và ∆3(x) là số hạng
Ix x
(3) không xác định tổng quát không tuyến tính của hệ, được
MIy Sx (Vsinα
Vαcosα Vqcosα)
biểu diễn dưới dạng sau:
L là lực nâng; D là lực cản; MA là momen khí động học 1(x) f1(x) g1 d1(t)
góc chúc ngóc.
2 (x) f2 (x) q d2 (t) (8)
3 (x) f3 (x) g2 e d3 (t)
L QSw (CL0 CLα α CLδe δe )
∆f1(x), ∆f2(x), ∆f3(x), ∆g1, ∆g3 là số hạng không xác định
D QSw (CD0 CDαα CDα2 α )
2
(4)
phi tuyến của tham số hệ thống.
qc
MA QSw cA (Cm0 Cmα α Cmδe δe Cmq A ) Giả thiết 1: Các biến trạng thái cần thiết của bộ điều
2V khiển thu được bằng phép đo. Có tồn tại hằng số đã biết
Trong đó Q = ρV2/2 là mật độ không khí; Sw là diện tích D∆ > 0 thỏa mãn i x D i 1, 2, 3 .
tham chiếu của cánh; cA là độ dài hình học trung bình của
cánh; δe là góc suy giảm của bề mặt nâng hạ; CL0 là hệ số lực Định nghĩa:
nâng của góc tấn tại gốc; CLα và CLδe là hệ số lực nâng với α θ θ - θ
d
và δe; CD0 là hệ số lực cản nâng tại gốc; CDα và CDα2 là lực cản α α - α
(9)
với α và α2. d
q q - qd
3. BỘ ĐIỀU KHIỂN BACKSTEPPING
Xem xét nhiễu bên ngoài, công thức (2) được viết như Trong đó, θd là tín hiệu đặt góc đường bay; αd và qd là tín
sau: hiệu điều khiển ảo; là sai số bám tín hiệu đặt góc đường
Website: https://jst-haui.vn Vol. 57 - No. 6 (Dec 2021) ● Journal of SCIENCE & TECHNOLOGY 31
- KHOA HỌC CÔNG NGHỆ P-ISSN 1859-3585 E-ISSN 2615-9619
bay; và q là các sai số bám tín hiệu điều khiển ảo. Sử Đạo hàm của sai số lọc như sau:
dụng bộ lọc bậc hai [2] có được θd và tín hiệu vi phân bậc
1
ζ α ζ +α
nhất của nó θ d :
k fα α d
(20)
θd ωn2
ζ 1 ζ +q
(10)
q
k fq q d
θ d0 s2 - 2ξ n ωn s ωn2
Theo biểu thức của αd và qd, từ giả thiết 2 có thể thấy
θd0 là đầu vào bộ lọc; ωn là băng thông; n là hệ số tắt
tồn tại các hằng số đã biết Bα0 > 0 và Bq0 > 0 thỏa mãn
dần. α d Bα0 ; q d Bq0 .
Giả thuyết 2: θd và θ d là các hàm bị chặn.
Hàm Lyapunov của thiết kế là:
Bước 1: Xét phương trình vi phân góc đường bay trong 1 2 1 2 1 2 1 2 1 2
hệ (7): Vn θ α q ζ α ζ q (21)
2 2 2 2 2
γ f10 (x) g10 α 1(x) (11)
Dựa trên các bước thiết kế trên, sử dụng giả thiết 1 và
Tín hiệu điều khiển ảo thiết kế là: X2 Y 2
bất đẳng thức XY
1 2 2
αd -k θ f10 (x) θ d (12)
g10 1 Ta có:
Trong đó k1 > 0 là tham số thiết kế. 1 1 2 1 2
Vn -k1 - θ 2 - k2 - 1- α - k3 - 1- q -
Bước 2: Xét phương trình vi phân góc tấn công trong hệ 2 2kfα 2kfq
(7): (22)
1 1 1 1
-
3
α f20 (x) q 2 (x) - ζ2α - - ζ2q D2 B2α0 B2q0
(13)
2kfα 2 2kfq 2 2
Tín hiệu điều khiển ảo thiết kế là:
Đặt:
q -k α g θ - f (x) α (14)
d 2 10 20 d 1 1
2k1 - , 2k2 -1- ,
Trong đó k2 > 0 là tham số thiết kế, α d là xấp xỉ của α d , 2 2kfα
z1 min
(23)
để lấy α sử dụng bộ lọc bậc nhất:
1 1 1
d
2k3 -1- , -1, -1
2kfq kfα kfq
k fα α d α d α d
(15)
α (0) α d (0) 3 2
d z2 D B2α0 B2q0 (24)
2
Trong đó, k fα 0 là hằng số thời gian của bộ lọc.
k1, k2, k3, kf và kfq thiết kế sao cho z1 > 0. Từ phương trình
Bước 3: Xét phương trình vi phân vận tốc góc trong hệ (24) và các bước thiết kế bộ điều khiển trên ta thấy z2 > 0 và
phương trình (7): bị giới hạn. Thế (23) và (24) vào (22) được:
q f30 (x) g20 δe 3 (x) (16) V -z V z (25)
n 1 n 2
Bộ điều khiển thiết kế là: Theo phương trình (25), có thể thấy sai số bám hội tụ tới
1 vùng lân cận gốc với bán kính là 2z2 /z1 . Khi sai lệch mô
δe -k q - a - f30 (x) q d (17)
g20 3
hình và nhiễu loạn bên ngoài nhỏ, bộ điều khiển danh
k3 > 0 là tham số thiết kế. Tương tự, sử dụng bộ lọc bậc nghĩa có thể đảm bảo sự ổn định bay. Khi sai lệch mô hình
và nhiễu loạn bên ngoài lớn, sai số bám tín hiệu đặt góc
một để lấy q d :
đường bay sẽ tăng lên. Đặc biệt trong quá trình biến đổi
hình dạng, sai số mô hình và nhiễu động chưa biết xuất
k fq qd qd qd
(18) hiện sẽ ảnh hưởng lớn hơn đến sự ổn định bay. Dưới đây
q (0) qd (0)
d tiến hành thiết kế bộ hiệu chỉnh, để bù cho bộ điều khiển
danh định, nhằm đảm bảo sự ổn định bay trong quá trình
k fq 0 là hằng số bộ lọc.
biến đổi hình dạng.
Định nghĩa sai số lọc như sau: 4. HIỆU CHỈNH RHO
ζ α - α Phương pháp Jacobian được sử dụng để tuyến tính hoá
α d d
(19) phương trình (7):
ζ q qd - qd
32 Tạp chí KHOA HỌC VÀ CÔNG NGHỆ ● Tập 57 - Số 6 (12/2021) Website: https://jst-haui.vn
- P-ISSN 1859-3585 E-ISSN 2615-9619 SCIENCE - TECHNOLOGY
x Ax B(δ δ - δ ) b Trong đó:
er en e0 0
(26)
Q CT C
y Cx
P r r - QP CrT C 0
Trong đó, x là lượng tăng của vectơ trạng thái đầy đủ x QR -QP C T Cr QP C T C 0
đối với giá trị tham chiếu x0, x 0 0 α0 0 0 ; α0 và ϑ0 0 0 QI
là góc tấn và góc chúc ngóc tương ứng; δen là đầu ra bộ Các tham số biến KR, kR thu được bằng cách giải phương
điều khiển backstepping; δer là đầu ra của bộ hiệu chỉnh; δe0 trình vi phân Riccati thể hiện trong công thức sau:
là góc lệch của cánh lái nâng hạ; b0 là sai lệch mô hình và K -K A - A TK R-1K B BTK - Q
R R R R R
(33) R R R R R
nhiễu bên ngoài; C = [1 0 0 0], A và B thay đổi theo sự thay
đổi hình dạng của cánh. k R -K R bR - ART k R R -1K RBRBRT K R (34)
Dạng không gian trạng thái của bộ lọc chỉ thị (10) như Vì KR, kR có thể được ổn định với bất kỳ giá trị nào tại
sau: biên giới của khoảng thời gian hữu hạn [t0, tf], giả sử KR, kR
x A x B θ tại biên thời gian tf là 0. Trong một thời gian hữu hạn, lấy
r r r r d0
(27) tích phân của (33) và (34) được KR(t0) và kR(t0). Độ lợi bộ hiệu
y r Cr x r
chỉnh được tính toán bằng cách sử dụng công thức sau:
Trong đó: F f -R -1B T K (t ) k (t ) (35)
R R R 0 R 0
θ -2ξ ω 1 Trong đó, FR = [Fr Fx FI], Fr, Fx và FI là độ lợi các bộ hiệu
x r d Ar n 2n
θd - ωn 0 chỉnh của trạng thái sai số tích phân, trạng thái máy bay và
trạng thái bộ lọc chỉ thị. Qua (36) tính được lượng hiệu
Br 0 ωn2
T
Cr 1 0 chỉnh:
Độ lợi bộ hiệu chỉnh được tính theo thời gian thực trong δer FR xR f (36)
một khoảng thời gian hữu hạn dựa trên các công thức (26) Phương pháp RHO không chỉ có những ưu điểm của
và (27) sao cho các chỉ số sau được giảm thiểu: phương pháp toàn phương tuyến tính, mà còn có thể dự
1
tf đoán các đặc tính đáp ứng mong muốn của đối tượng, độ
QP yr - y QI yI2 Rδer2 dt
2 t
2
min JR (28) chính xác của mô hình không đòi hỏi cao, độ lợi các bộ hiệu
δer
0 chỉnh điều chỉnh thời gian thực cho trị sai lệch giữa tín hiệu
Trong đó, tf và t0 là giới hạn trên và dưới của miền thời gian chỉ thị và góc đường bay, để hạn chế tác động của nhiễu
hữu hạn, Qp và QI là các giá trị trọng số của sai số bám và tích bên ngoài và sai lệch mô hình. Cấu trúc bộ hiệu chỉnh thể
phân sai số bám, R là giá trị trọng số của lượng hiệu chỉnh, yI là hiện trong hình 2.
sai số tích phân và có thể thu được từ công thức sau: od xr BỘ ĐIỀU KHIỂN
BỘ LỌC
BACKSTEPPING
x I =y r y
(29) en
yI =xI en
BỘ HIỆU CHỈNH RHO e
Kết hợp phương trình trạng thái tuyến tính của máy bay Fr , Fx , FI , f
với phương trình trạng thái của bộ lọc và phương trình sai số
tích phân, nhận được phương trình trạng thái mở rộng là: x, y, A, B
x r A r 0 0 x r 0 Hình 2. Cấu trúc bộ hiệu chỉnh
x 0 A 0 Δx + B δ er
5. MÔ PHỎNG VÀ NHẬN XÉT
x C -C 0 x 0
I r I
(30)
Br θ d0 0
Góc gập cánh [độ]
+ B(δ en -δ e0 ) + b 0
0 0
Phương trình (30) viết lại như sau:
x R AR xR BR δer bR b (31)
Phương trình (28) có thể được viết lại:
tf
1
min JR xRT QR xr Rδ2er dt (32) Thời gian [s]
δR 2t Hình 3. Đường cong góc gập cánh
0
Website: https://jst-haui.vn Vol. 57 - No. 6 (Dec 2021) ● Journal of SCIENCE & TECHNOLOGY 33
- KHOA HỌC CÔNG NGHỆ P-ISSN 1859-3585 E-ISSN 2615-9619
Tham số mô phỏng: Diện tích cánh: S = 23m2; Trọng Khi có nhiễu loạn khí động học sử dụng tín hiệu nhiễu
lượng: G = 7300kg; Momen quán tính: Jx = 630KgmS2; d3(t) = 5sin(t) (°/s2). Đáp ứng góc đường bay thể hiện trên
Jz = 5250KgmS2 ; Dây cung khí động trung bình: bA = 4m; hình 7, 8, 9.
Kích thước khí động: L = 7,15m; Độ cao bay H = 5000m,
M = 0,5; Tốc độ V = 35m/s. Trong quá trình bay, các thông
số khí động học thay đổi khi điều kiện bay và góc xoè cánh
thay đổi. Sai số mô hình 20% được cộng thêm vào thông số
khí động học. Khi cho máy bay thay đổi góc gập cánh tại
thời điểm 14s, đường cong biến đổi cấu hình cánh được thể
hiện trong hình 3.
Khi không có nhiễu khí động học bộ điều khiển
Backstepping đáp ứng tốt các góc điều khiển. Đồ thị đáp
ứng góc θ, ϑ, α thể hiện trên hình 4, 5, 6.
Hình 7. Đáp ứng góc θ khi có nhiễu
Hình 4. Đáp ứng góc nghiêng quỹ đạo θ
Hình 8. Đáp ứng góc ϑ khi có nhiễu
Hình 5. Đáp ứng góc chúc ngóc ϑ
Hình 9. Đáp ứng góc α khi có nhiễu
Trong khoảng thời gian từ 14 ~ 27s, máy bay thay đổi
hình dạng cánh để chuyển sang cơ động, góc đường bay
biến đổi mạnh. Khi có tác động của nhiễu, bộ hiệu chỉnh
RHO giúp loại bỏ phần lớn tín hiệu nhiễu, góc đường bay
đạt được giá trị xác lập trong thời gian ngắn, hạn chế ảnh
Hình 6. Đáp ứng góc tấn α hưởng bởi quá trình biến hình của cánh.
34 Tạp chí KHOA HỌC VÀ CÔNG NGHỆ ● Tập 57 - Số 6 (12/2021) Website: https://jst-haui.vn
- P-ISSN 1859-3585 E-ISSN 2615-9619 SCIENCE - TECHNOLOGY
6. KẾT LUẬN
Bộ điều khiển Backstepping kết hợp hiệu chỉnh RHO
cho hiệu quả cao khi ổn định đường bay của máy bay biến
đổi cấu hình cánh, một hệ thống phi tuyến và có nhiều
tham số bất định. Trong quá trình biến đổi hình dạng, các
đặc tính động của máy bay sẽ thay đổi rất nhiều, bộ hiệu
chỉnh RHO kết hợp giữa hiệu chỉnh phản hồi và tối ưu hóa
từng đoạn cho mô hình, bù đầu ra của bộ điều khiển
backstepping, đảm bảo sự ổn định đường bay của máy bay
trong quá trình thay đổi hình dạng, nâng cao tính bền
vững, giảm ảnh hưởng của sai lệch mô hình và các nhiễu
loạn bên ngoài.
TÀI LIỆU THAM KHẢO
[1]. Yue T, Wang L X, Ai J Q, 2013. Gain self-scheduled H∞ control for
morphing aircraft in the wing transition process based on an LPV model. J. Chinese
Journal of Aeronautics, Vol.26, No.4, pp.909-917.
[2]. Sonneveldt L, Chu Q P, Mulder J A. 2007. Nonlinear flight control design
using constrained adaptive backstepping. Journal of Guidance, Control, and
Dynamics, 30 (2): 322-335.
[3]. T. Lee, Y. Kim, 2001. Nonlinear adaptive flight control using backstepping
and neural networks controller. Jnl. Guidance, Control, and Dynamics, 24(4):675–
682.
[4]. M. Sharma, D. G. Ward, 2002. Flight-path angle control via neuro-
adaptive backstepping. In Proceedings of the AIAA Guidance, Navigation, and
Control Conference, AIAA-2002-4451, Monterey, Californ.
[5]. Sangchul Lee, Kwangjin Kim, 2012. A Sliding Mode Control with
Optimized Sliding Surface for Aircraft Pitch Axis Control System. The Japan Society
for Aeronautical and Space Sciences, Vol. 55, No. 2, pp. 94-98.
[6]. Sridhar Seshagiri, Ekprasit Promtun, 2018. Sliding Mode Control of F-16
Longitudinal Dynamics. 2018 American Control Conference Westin Seattle Hotel,
Seattle, Washington, USA June 11-13, 2018.
AUTHORS INFORMATION
Nguyen Ngoc Tuan1, Tran Xuan Tinh2,
Nguyen Tuan Anh2, Tran Thi Nga2
1
Military Technical Academy
2
Air Defence - Air Force Academy
Website: https://jst-haui.vn Vol. 57 - No. 6 (Dec 2021) ● Journal of SCIENCE & TECHNOLOGY 35
nguon tai.lieu . vn