Tổng hợp bộ điều khiển trượt đầu cuối nhanh cho hệ truyền động nhiều động cơ có liên hệ khe hở, ma sát, đàn hồi

Nghiên cứu khoa học công nghệ 
Tạp chí Nghiên cứu KH&CN quân sự, Số 66, 4 - 2020 33
TỔNG HỢP BỘ ĐIỀU KHIỂN TRƯỢT ĐẦU CUỐI NHANH 
CHO HỆ TRUYỀN ĐỘNG NHIỀU ĐỘNG CƠ 
CÓ LIÊN HỆ KHE HỞ, MA SÁT, ĐÀN HỒI 
Trần Xuân Tình1*, Phạm Tuấn Thành1, Trần Văn Tuyên1, Đào Sỹ Luật2 
Tóm tắt: Bài báo trình bày kết quả tổng hợp bộ điều khiển trượt đầu cuối nhanh 
trên cơ sở bất đẳng thức ma trận cho hệ cơ điện nhiều động cơ có liên hệ ma sát, 
đàn hồi. Các kết quả được khảo sát đánh giá bằng mô phỏng trên phần mềm 
Matlab-Simulink cho thấy bộ điều khiển này đảm bảo được các yêu cầu chất lượng 
ngay cả khi hệ thống chịu ảnh hưởng của các yếu tố phi tuyến do cấu trúc phần cơ 
gây ra. 
Từ khóa: Nhiều động cơ; Điều khiển trượt đầu cuối; Bất đẳng thức ma trận tuyến tính. 
1. ĐẶT VẤN ĐỀ 
 Hệ truyền động (HTĐ) nhiều động cơ ứng dụng trong công nghiệp và quốc phòng đều 
là hệ động lực học phi tuyến, chứa các liên hệ chéo; các mối liên hệ này làm cho mô hình 
của đối tượng điều khiển trở nên phi tuyến. Trong điều khiển HTĐ, vật liệu đàn hồi nhiều 
động cơ, việc kiểm soát lực căng trên băng vật liệu là một vấn đề khó. Các bộ điều khiển 
(BĐK) HTĐ này liên tục được nghiên cứu phát triển từ đơn giản như BĐK PID [3] đến 
phức tạp như logic mờ [4], mạng nơ ron [5], điều khiển tối ưu [6] và điều khiển bền vững 
[7]. Tuy nhiên, vẫn chưa có BĐK nào đáp ứng tốt chất lượng của hệ cũng như tính hết 
những tác động của yếu tố phi tuyến của phần cơ. BĐK điều khiển trượt đầu cuối nhanh 
trên cơ sở bất đẳng thức ma trận tuyến tính LMI-FTSM (Linear matrix inequalities-Fast 
Terminal Sliding Mode) được tổng hợp cho HTĐ băng vật liệu đàn hồi nhiều động cơ đã 
cho chất lượng tốt ngay cả khi tính đến các yếu tố phi tuyền như: ma sát, khe hở, đàn hồi. 
2. TỔNG HỢP MÔ HÌNH CƠ HỆ 
Xét một mô hình hệ thống điện cơ hai khối lượng đàn hồi dạng 1 kết hợp dạng 2 [1], 
[2] như hình 1. Ở đây, chỉ xét dao động đàn hồi trong liên kết hai chiều giữa động cơ và 
tải; giữa hai tải của hai động cơ, các mối liên kết khác coi như cứng vững hoàn toàn. 
F
1 1,T J
1
1sT
1r
H 
1 1,c b
1 1,L LT J 1r
1v
2r
2v
2 2,L LT J
2 2,T J
2
2sT
2r
H 
2 2,c b
Hình 1. Mô hình cơ hệ đàn hồi hai động cơ. 
Kỹ thuật điều khiển & Điện tử 
T. X. Tình, , Đ. S. Luật, “Tổng hợp bộ điều khiển trượt  khe hở, ma sát, đàn hồi.” 34 
Trong đó: T1, T2, TL1, TL2 là mômen động cơ,tải; Ts là mômen xoắn trên trục; J1, J2, JL1, 
JL2 là mômen quán tính động cơ, tải; ω1, ω2, ωL1, ωL2 là tốc độ động cơ, tốc độ tải; c1, c2, b1, 
b2 là hệ số cứng, hệ số ma sát nhớt của khớp nối; H là góc khe hở bánh răng; F là lực 
căng của dải vật liệu; r1, r2 tương ứng là bán kính rulo 1, rulo 2; v1, v2 là tốc độ dài của dải 
vật liệu chạy qua. Chỉ số 1, 2 tương ứng của HTĐ động cơ 1 và động cơ 2. Theo [1], [2] có 
hệ phương trình trạng thái của HTĐ băng vật liêu đàn hồi hai động cơ: 
 
 
2 2 5
1 12 1 5 12 2 6
2 2 5
2 12 1 5 12 2 6
3 1 5
4 2 6
5 1 1 3 1 1 3 1 1 2
1
6 2 2 4 2 2 4 2 2 1
2
1
. ( ) . ( ) ( )
1
. ( ) . ( ) ( )
c b L
L
c b L
L
r x x
x C r x C r x
l
r x x
x C r x C r x
l
x u x
x u x
x k f x k g x T r x
K
x k f x k g x T r x
K




 
  
(1) 
Trong đó:     61 2 3 4 5 6 21 12 1 2 1 2
T T
L Lx x x x x x F F   x ; 
    1 21 2 1 2
T T
r ru u  
 u ; Các hàm 3 4 3 4( ), ( ), ( ), ( )f x f x g x g x  là hàm phi tuyến 
phụ thuộc vào khe hở bánh răng hộp số. 
3
3 3 3
3 3
0 khi 
( ) khi >
 khi <-
H
H H
H H
x
f x x x
x x
; 3 3 3( ) . ( )g x x f x 
  ; với 33
3
0 khi 
( )
1 khi 
H
H
x
f x
x
 
Đặt:   1 21 1 2
T
x x z ;   1 42 3 4 5 6
T
x x x x z 
 11
0 0
0 0
A ; 
12 1 12 2
12
12 1 12 2
0 0
0 0
C r C r
C r C r
A ; 1
0 0 0 0
0 0 0 0
B ; 
22
0 0 1 0
0 0 0 1
0 0 0 0
0 0 0 0
A ; 2
1 0 0 0
0 1 0 0
0 0 0 0
0 0 0 0
B ; 121
1
2
2
0 0
0 0
0
0
L
L
r
K
r
K
A ;
  
 
1 1 3 1 1 3 12 2
1
2 2 4 2 2 4 2
2
0
0
1
. ( ) . ( ),
1
. ( ) . ( )
c b L
L
c b L
L
k f x k g x Tt
K
k f x k g x T
K
d z 

(2) 
Nghiên cứu khoa học công nghệ 
Tạp chí Nghiên cứu KH&CN quân sự, Số 66, 4 - 2020 35
Hệ thống (1) được viết lại thành: 
1 11 1 12 2
2 21 1 22 2 ( , t)
z A z A z
z A z A z Bu d z
y Cz

 (3) 
3. TỔNG HỢP BỘ ĐIỀU KHIỂN LMI-FTSM 
Đặt ,ref e z z khi đó chọn mặt trượt: 
 12 1
T S e Λe A Pν e (4) 
Trong đó,   Λ H I với H, I là ma trận 4 × 2 và 4 × 4, 11 1

 ν e e . 
Giả sử rằng  và k là các hằng số dương, trong đó (t) ; f k e . Nếu có các số 
0 , 1 0 và 2 0 và các ma trận 0, 0, 0, 0X Q W G và Y với kích thước 
phù hợp để đáp ứng các điều kiện LMI tiếp theo: 
11 11 12 12 0
T T T 
A X + XA - A Y - Y A + Q X
X -W
 (5) 
12 122 0
T 
A A X
X G
; 1 2 2 0 I Q ; 2 2 2 0 I W ;
2
1 22
0
k

  (6) 
thì hệ (3) sẽ ổn định tiệm cận với thời gian hữu hạn và sai lệch 1 t k e . Với quy ước 
các bất đẳng thức ma trận có dấu 0 là xác định dương. 
Chứng minh: 
Để chứng minh điều kiện LMI trên cần sử dụng 2 bổ đề sau: 
Bổ đề 1 [8]. Cho 2 ma trận hằng số A và B có kích thước phù hợp, với bất kỳ ma trận xác 
định dương Q > 0, ta có bất đẳng thức sau: 
1T T T T A B B A A QA B Q B (7) 
Hay 1T T T T A B + B A - A QA - B Q B xác định âm hoặc bằng 0. 
Bổ đề 2 [9]. Giả sử rằng hàm V(t) dương xác định liên tục thỏa mãn bất đẳng thức vi phân sau: 
 0 0 t t , 0V t V t V t V t
    (8) 
Trong đó α và β là hằng số dương và  là tỷ lệ của hai số nguyên dương lẻ với 1 0. 
Khi đó, với bất kỳ giá trị 0t đã cho nào, V t sẽ hội tụ về 0 trong thời gian hữu hạn: 
 1 0
0
1
ln
1
r
V t
t t
 
  
 (9) 
Đặt 1 P = X và 1 H = YX , P là ma trận đối xứng xác định dương. Chọn hàm ứng viên 
Lyapunov như sau: 
 1 1 1 1 1 1 1 1 1 1( )
T T TV V e e Pe e e Pe e Pe   (10) 
Từ (3) có: 
1 11 1 12 2 1 1 11 1 12 2
1 1 11 12 12 12 1
m
T t
z = A z + A z e - z = A z + A z
e e A HA A A Pν e
 

 (11) 
Kỹ thuật điều khiển & Điện tử 
T. X. Tình, , Đ. S. Luật, “Tổng hợp bộ điều khiển trượt  khe hở, ma sát, đàn hồi.” 36 
Trong đó: 11 1 12 2 1m m mt A z A z z 
Thay (11) vào (10) ta có: 
1 1 1 1 11 12 12 12 1
1 11 12 12 12 1 1 1 1 +
T T
TT T T
V
t t 
e e P e A HA A A Pν e
e A HA A A Pν e Pe e P Pe

 (12) 
Dựa trên bổ đề 1, thu được các bất đẳng thức sau: 
1
1 1 1 1
2 1
1 1 x 
T T T T
T
ma
t t t t t t t t
t t
   
 
e P Pe e PQPe Q
e PQPe Q
 (13) 
Trong đó, Q là ma trận xác định dương tùy ý và 1xma Q là giá trị riêng lớn nhất của 
1 Q . Thay thế (13) vào (12) theo đó: 
1 1 1 11 12 11 12 1
2 1
1 12 12 1 1 12 12 1 x 
TT
T T T T
ma
V
Q  
e e P A HA A HA P PQP e
e PA A Pν e ν e PA A Pe

 (14) 
Với 1 11 2 1
T
k
k
e e
e và 
1
1 1

 ν e e , (14) có thể viết lại: 
2
1
1 1 1 11 12 11 12 x 12
1 1
1 12 12 1 1 12 12 1 
TT T
ma
T
T T T
V
r
 


e e PA PA H A P HA P PQP Q I e
e PA A P e e PA A Pe

 (15) 
Giả sử rằng các bất đẳng thức sau đây được thỏa mãn: 
1
12 122 0
T PA A P Φ (16) 
 111 12 11 12
T P A HA A HA P PQP W (17) 
Khi đó, (15) có thể viết: 
2
11 1 1
1 1 1 1 1 x 1 1 12
2
21 1 1
min x 1 min 12
 -
T T T
ma n m
ma
V I
k
k





  
e e W e e Q e e Φ e
W Q e Φ e

 (18) 
Trong đó, min . là giá trị riêng nhỏ nhất. Khi đó, (18) có thể viết lại như sau: 
 11 1 1 1 1 1 1 1V V V
  e e e (19) 
Với 
2
1 1
2min x
1
x
ma
ma
r
 

W Q
P
; 
1
1
1 min
x
1
ma

 

Φ
P
; 1
2

 
Kết hợp với (6) có: 
2 2
1 1
min x min x2 2
0 0ma ma
k k
 
    Q W W Q 
Nghiên cứu khoa học công nghệ 
Tạp chí Nghiên cứu KH&CN quân sự, Số 66, 4 - 2020 37
Suy ra 1 0 . Giả sử, 
1 X P và nhân X vào trước và sau các số hạng của (17) ta được: 
 111 12 11 12
TT A X HA X XA HA X Q XW X (20) 
Đặt Y HX khi đó từ (20) suy ra được (5) và nếu có điều kiện LMI (6) thì (16) thỏa mãn. 
Từ (19) cho thấy, nếu 1 k e thì 1 1 0V e khi đó sai lệch 1 te sẽ tiến tiệm cận về 0, nằm 
trong vùng không gian 1 1 1:t t k e e với thời gian hữu hạn 1rt , theo bổ đề 2 là: 
 111 1 0 1
1
1 1 1
1
ln
1
r
V t
t
 
  
 (21) 
Kết luận: Nếu BĐK thỏa mãn điều kiện LMI đã đặt ra thì sai lệch 1 te sẽ tiệm cận về 
0 trong thời gian hữu hạn. 
Cũng có thể rút ra từ (4) rằng 2 te sẽ tiến tới vùng lân cận 0 
 12 2 2 12: Tt t r k  e e H A P trong thời gian hữu hạn 2rt . Để thực hiện điều 
đó, tiến hành chọn luật điều khiển như sau: 
 12n st t t
 u u B u (22) 
 11
2 12
T
n d
d
t
dt
ν e
u B ΛAz A P Λz (23) 
 2 sgn sgns t d

  u S e S e S e S e (24) 
với ,   là các hằng số dương. 
Chứng minh: Chúng ta đi chứng minh với luật điều khiển đã chỉ ra hệ thống sẽ tiến về 
mặt trượt (4) và nằm trên đó trong thời gian hữu hạn. 
Xét hàm xác định dương sau: 2
1
2
TV S e S e S e 
1
2 1 12 21 1 22 2 2
1
1 2 2 12 2 2 , = ,
T T T
T T
d d d
d
V
dt
d
d t t
dt
g e
S e S e S e S e Hz A P A z A z B u
g e
Hz z z S e ΛAz A P B u Λz d z
 
  
Thay (22) đến (24) vào ta được: 
 22 2sgn sgn ,T T T TV d t

  S e S e S e S e S e S e S e S e S e d z 
Vì rằng 2 2 x, mad t d z nên có: 
 2
2 1
2 min min 2 2 2 2V V V
     
 S e S e S e S e S e 
Chọn 1 /22 min 2 min 22 0, 2 0, 1 / 2 1
        . Theo bổ đề 2, hệ 
thống sẽ tiến về mặt trượt trong thời gian hữu hạn 2rt được xác định như sau: 
 2
1
2 2 0 2
2
2 2 2
1
ln
1
r
V t
t

 
  
S e
Với 2 0 0 01/ 2 TV t t t S e S e S e 
Kỹ thuật điều khiển & Điện tử 
T. X. Tình, , Đ. S. Luật, “Tổng hợp bộ điều khiển trượt  khe hở, ma sát, đàn hồi.” 38 
Kết luận: Từ các kết quả chứng minh trên cho thấy các sai lệch 1 2; t te e sẽ chuyển 
động tiệp cận về vùng lân cận gốc 0 trong khoảng thơi gian hữu hạn 1 2r r rt t t . 
4. MÔ PHỎNG VÀ ĐÁNH GIÁ KẾT QUẢ 
4.1. Tham số mô phỏng 
Trước hết, tính toán các tham số của HTĐ dùng động cơ không đồng bộ ba pha. Xét hai 
động cơ ba pha roto lồng sóc của hãng Siemens có các tham số giống nhau, cụ thể như sau: 
Công suất Pđm = 4kW; Lm = 0,1958(H); Ls = 0,202(H) ; Lr = 0,2065; Rr = 1,275(); Rs 
=1,663() ; p = 2; nđm = 1400 v/p; J1,2 = 7,47.10
-5Kgm2; JL1,2=8.258.10
-5; Bán kính rulo tải 
1, 2 r1, r2 =0,03m; Hệ số cứng của vật liệu C12 =0,4N/m; Hệ số đàn hồi c=360Nm/rad; Hệ 
số ma sát nhớt b=0,02; Chiều dài dải băng vật liệu =0,6m. 
Dựa trên thông số thực tế của HTĐ hai động cơ trong phụ lục tiến hành chọn: 2 5d , 
 20, 200, 0.2, 0.05, 0.6k    . Tính toán các ma trận A11, A12, A21, 
A22 , B1 và B2. 
1 1,T J
1r
1sT
1L
H 
1 1,c b
1 1,L LT J
1r
1v
2r
2v
2 2,L LT J
2 2,T J
2r
2sT
2L
H 
2 2,c b
12 12 1 1 2 2 12
12
1
1r rF C r r F
C l
 

xˆ dk
u
datF
1r
2r
1 2 3 3
2
2 1 1 3 3 3
3
3 12 1 21 2 22 3 3 3
12
3
ˆ ˆ ˆ3 ( )
1ˆ ˆ ˆ ˆ ˆ. ( ) . ( ) (T . ) .( )
1 1ˆ ˆˆ ˆ (1 ) ( )
.
c b L
L
L
k F k G r
K
C r x r x r
C l J



X u X X X
X X X X X X
X X X X
y X

 

1 1
2 121
2 0.6
=
 +5sgn 200sgn 20
T
n s d
dk
d
t tu
dt
u
ν e
u u B ΛAz A P Λz
u
S e S e S e S e

Fˆ
1u 2u
2f1f

1ik 2i
k
21F 12F
Hình 2. HTĐ nhiều động cơ với BĐK LMI-FTSM và bộ quan sát HGO. 
Giải điều kiện (5,6) bằng LMI toolbox của Matlab xác định được: 
 1
2.1758 1.17
1.17 2.1758
P X 
; 1
0 0
0 0
2.1584 2.1584
1.85 1.85
H YX 
Qua đó, xác định được luật điều khiển uđk. Cấu trúc tổng thể của hệ thống với BĐK 
LMI-FTSM và bộ quan sát HGO (High Gain Observer) đã tình bày trong [1, 2] được trình 
Nghiên cứu khoa học công nghệ 
Tạp chí Nghiên cứu KH&CN quân sự, Số 66, 4 - 2020 39
bày trên hình 2. 
4.2. Sơ đồ mô phỏng 
Sơ đồ mô phỏng hình 3 được xây dựng trên công cụ Matlab-Simulink trong đó: động 
cơ, khớp nối, tải được xây dựng theo (1). BĐK mô phỏng theo công thức (22). 
Hình 3. Mô hình mô phỏng hệ thống và BĐK. 
4.3. Kết quả mô phỏng 
Trường hợp 1: Đáp ứng lực căng khi khe hở bánh răng αH = 0,008 rad (0,458
0) 
(a) 
(b) 
Hình 4. Đáp ứng lực căng với khe hở bánh răng αH = 0,008 rad, 
 tín hiệu vào là hàm bước nhảy (a) và hàm xung (b). 
Trường hợp 2: Đáp ứng lực căng khi khe hở bánh răng αH = 0,08 rad (4,58
0). 
(a) 
(b) 
Hình 5. Đáp ứng lực căng với khe hở bánh răng αH 0,08 rad, 
 tín hiệu vào là hằng số (a )và hàm nhảy bậc (b). 
Kỹ thuật điều khiển & Điện tử 
T. X. Tình, , Đ. S. Luật, “Tổng hợp bộ điều khiển trượt  khe hở, ma sát, đàn hồi.” 40 
Trường hợp 3: Đáp ứng lực căng của hệ thống khi xét đến ảnh hưởng của mômen 
quán tính của tải 
(a) 
(b) 
Hình 6. Đáp ứng lực căng với ảnh hưởng của mômen quán tính JL=4J (a) 
và JL1=4J; JL2=J(b). 
Trường hợp 4: Đáp ứng lực căng khi tải thay đổi hàm nhảy bậc TL1 = 100Nm; TL2 = 
50Nm và khi thay đổi độ cứng của băng vật liệu. 
(a) 
(b) 
Hình 7. Đáp ứng lực căng với ảnh hưởng của tải thay đổi (a) và 
khi giảm hệ số cứng băng vật liệu 12 0.3C (N/m)(b). 
Nhận xét: 
TT Trường hợp 
Bộ LMI-FTSM 
ĐQC Tqđ(s) SSXL 
1 αH = 0.008, Fđ =10N 0 1.2 0 
5 αH = 0.008, lực căng đặt là hàm xung 0 0.8 0 
6 αH = 0.08, Fđ =10N 2% 2 0 
7 αH = 0.08, lực căng đặt là hàm Step 8% 2.5 0 
11 Ma sát coulomb 0,1; ma sát nhớt 0.05 0 4 0 
13 Mômen quán tính JL=4J 7% 3.2 0 
14 Mômen quán tính JL1=4J; JL2=J 12% 2.6 0 
15 TL1 = 100Nm; TL2 = 50Nm 8% 2.5 0 
16 Hệ số cứng băng vật liệu 12 0.3C (N/m) 0 1.8 0 
5. KẾT LUẬN 
Bài báo đã trình bày kết quả tổng hợp BĐK LMI-FTSM cho hệ cơ điện nhiều động cơ 
có liên hệ ma sát, đàn hồi. Phần trình bày được bắt đầu từ việc xây dựng mô hình cơ hệ, 
đặt điều kiện LMI, chứng minh điều kiện, xây dựng mô hình mô phỏng, kiểm nghiệm 
Nghiên cứu khoa học công nghệ 
Tạp chí Nghiên cứu KH&CN quân sự, Số 66, 4 - 2020 41
bằng phần mềm Matlab-Simulink. Qua kiểm tra và so sánh với các kết quả của các công 
bố trước đây [3-7] cho thấy, BĐK đã nâng cao được chất lượng của hệ thống truyền động 
nhiều động cơ thông qua các tiêu chí đánh giá, đó là: tính bền vững với nhiễu, đảm bảo 
khả năng quá tải về mô men, đảm bảo tốt đặc tính động và đặc tính tĩnh, đảm bảo tính 
chính xác bám. 
TÀI LIỆU THAM KHẢO 
[1]. Dao Phuong Nam, Pham Tuan Thanh, Tran Xuan Tinh, Tran Thanh Dat, Pham Van 
Tu, “High-Gain Observer based Output feedback Controller for a Two-Motor Drive 
System: A Separation Principle Approach”, Lecture Note in Electrical Engineering 
465, Scopus Q3, Dec-2017. 
[2]. Pham Tam Thanh, Dao Phuong Nam, Tran Xuan Tinh and Luong Cong Nho, “High-
Gain Observer–Based Sliding Mode Control of Multimotor Drive Systems”, book 
Adaptive Robust Control Systems, Published by InTech Janeza Trdine 9, 51000 
Rijeka, Croatia 2018. 
[3]. B. Allaoua , A. Laoufi and B. Gasbaoui, “Multi-Drive Paper System Control Based on 
Multi-Input Multi-Output PID Controller”, Leonardo Journal of Sciences, 2010. 
[4]. Fawzan Salem, E.H.E. Bayoumi, “Robust fuzzy-PID control of three-motor drive system 
using simulated annealing optimization”, Journal of Electrical Engineering, 2011. 
[5]. Li Jinmei, Liu Xingqiao, “Application of an Adaptive Controller with a Single Neuron 
in Control of Multi-motor Synchronous System”, IEEE, 2008. 
[6]. A. Angermann, M. Aicher, and D. Schroder, “Time-optimal tension control 
forprocessing plants with continuous moving webs”, Proc. 35th Annual Meeting- 
IEEE Industry Applications Society, Rome, Oct. 1999. 
[7]. H. Koc, D. Knittel, M. D. Mathelin, “Robust gain-scheduled control of winding 
systems”, IEEE Conf. Decision and Control, Sidney, Australia, Dec. 2000. 
[8]. C. Scherer, S. Weiland, “Linear Matrix Inequalities in Control, Dutch Institute of 
Systems and Control (DISC)”, The Netherlands, 2005. 
[9]. J. Zhang, “Solving Linear Matrix Inequality (LMI) Problems”, Advances in Industrial 
Control, Springer International Publishing Switzerland 2016, PP 203-224. 
ABSTRACT 
DESIGN FAST TERMINAL SLIDING MODE CONTROLLER FOR ELECTRIC 
DRIVETRAIN MULTI-MOTOR HAVE CONTACTS FRICTION, ELASTIC 
The paper presents the results of a fast terminal sliding controller based on 
linear matrix inequality for many motor systems with friction, elasticity. The results 
surveyed and evaluated by simulation on Matlab-Simulink software show that the 
controller ensures the quality requirements in the system under the influence of 
nonlinear factors as causing the mechanical structure. 
Keywords: Multi-Motor, Terminal Sliding Mode Control, linear matrix inequalities. 
Nhận bài ngày 29 tháng 02 năm 2020 
Hoàn thiện ngày 24 tháng 3 năm 2020 
Chấp nhận đăng ngày 10 tháng 4 năm 2020 
Địa chỉ: 1Khoa Kỹ thuật điều khiển, Học viện Kỹ thuật quân sự; 
 2Đại học Đồng Nai, số 4 Khu phố 3 Lê Quý Đôn, Biên Hòa, Đồng Nai. 
 *Email: tinhpk79@gmail.com.