Ổ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

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

Thêm vào BST

Báo xấu

19
lượt xem 2
download

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

Để đả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 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 Backstepping được sử dụng bảo đảm tính ổn định quỹ đạo bay. Hiệu chỉnh góc đường bay theo phương pháp tối ưu hoá từng đoạn.

Chủ đề:

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

Lưu

Nội dung Text: Ổ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

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   k1   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)    IyqPZT MIy      θ  1 (LPsinα mgcosθFIz )    mV   1   g1  QSwCLα α  1 (LPsinα 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α Sx 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 Vn  -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   2k1 - , 2k2 -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 2k3 -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  Ax  B(δ  δ - δ )  b Trong đó:  er en e0 0  (26)  Q CT C y  Cx   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