Tải bản đầy đủ (.docx) (37 trang)

tiểu luận môn học robot công nghiệp

Bạn đang xem bản rút gọn của tài liệu. Xem và tải ngay bản đầy đủ của tài liệu tại đây (826.69 KB, 37 trang )

Robot Công Nghiệp GVHD: TS Nguyễn Sĩ Dũng
ĐỀ TÀI TIỂU LUẬN
Cho tay máy như hình sau :
1. Vẽ sơ đồ nguyên lý của tay máy, xây dựng hệ trục và bảng động học.
2. Giải bài toán vận tốc tay gắp.
3. Giải bài toán động lực học.
I. XÂY DỰNG HỆ TRỤC VÀ BẢNG DH
1. Xây dựng hệ trục:
Page 1
Robot Công Nghiệp GVHD: TS Nguyễn Sĩ Dũng
2. Bảng DH :
DL
Khâu
i
γ

i
α

a
i
d
i
1
1
β

90
o

0 h


1
2
1
β
0 h
2
0
3
1
β
90
o
0 0
4 0 0 0 h
3
+h
II. TÍNH VẬN TỐC TAY GẮP:
Ta có công thức tính vận tốc như sau:
.x J q=
& &
Page 2
Robot Công Nghiệp GVHD: TS Nguyễn Sĩ Dũng
Trong đó
1
2
3
q
d
β
β

β
 
 
 
=
 
 
 
 
&
&
&
&
&
và J là ma trận Jacobi nên ta có phương trình vận tốc
tay máy như sau:
[ ]
0
0
0
0
0
0
3
4
1
2
1 2 3

x

y
z
x
y
z
v
v
v
x J J
d
J J
β
β
ω
ω
ω
β
 
 
 
 
 
= =
 
 
 
 
 

 

 
 
 
 
 



&
&
&
&
&

Trong đó :
0 0
4
1
0
z p
J
z
 
×
=
 
 

( ) ( )
( )

1 0 1 0
4
2
1 0
z p
J
z
 
×
=
 
 
 

( ) ( )
( )
2 0 2 0
4
3
2 0
z p
J
z
 
×
=
 
 
 


( )
3 0
4
0
z
J
 
=
 
 
 Tìm
( ) ( ) ( )
1 0 2 0 3 0
0
, , ,z z z z
. Sử dụng kí hiệu s = sin, c = cos, ta có :

0
0
0
1
z
 
 
=
 
 
 

1 1

1 1
0
1
1
0 0
0 0
0 1 0
0 0 0 1
c s
s c
A
h
β β
β β
 
 

 
=
 
 
 
1 1
0
1 1 1
0 s
s 0 c
0 1 0
c
R

β β
β β
 
 
⇒ = −
 
 
 
( )
1
1 0
0
1 1
0
0
1 0
s
z R c
β
β
   
   
⇒ = × = −
   
   
   

Page 3
Robot Công Nghiệp GVHD: TS Nguyễn Sĩ Dũng


1 2 1 2 1 1 2 2
1 2 1 2 1 1 2 2
0 0 1
2 1 2
2 2 1 2 1
. . . .
. . . .
0 .
0 0 0 1
c c c s s c c h
s c s s c s c h
A A A
s c s h h
β β β β β β β
β β β β β β β
β β β

 
 
− −
 
= × =
 
+
 
 
1 2 1 2 1
0
2 1 2 1 2 1
2 2

. .
. .
0
c c c s s
R s c s s c
s c
β β β β β
β β β β β
β β

 
 
⇒ = − −
 
 
 

( )
1
2 0
0
2 1
0
0
1 0
s
z R c
β
β
   

   
= × = −
   
  
   




( ) ( )
( ) ( )
( ) ( )
1 3 1 1 3 1 2 2
0 0 1 2
1 3 1 1 3 1 2 2
3 1 2 3
3 3 1 2
2
1
2 2
2
2 2
. .
. .
0 .
0 0
.
.
0 1
.

.
. .
c s c c c h
s c s c c h
A A A A
c
s c
c
h
s
s h
s
β β β β β β β β β
β β β β β β β β β
β β β β β
+ +
+
 
 

 
= =
 
− +

+

+
+




( ) ( )
( ) ( )
( ) ( )
1 3 1 1 3
1 3 1 1 3
2 2
0
3 2 2
2 23 3
. .
. .
0
c s c
s c s
c
c
s
c s
R s
β β β β β β β
β β β β β β β
β β β β

+ + 
 
⇒ = + +





 
+ +


( )
1 2 3
0
3
2 3
0
3
3
1 2
0 . ( )
0 . ( )
1 ( )
c s
R s sz
c
β β β
β β β
β β
+
   
   
× = +
   
   

− +
   
⇒ =

 Tìm
( ) ( ) ( )
3 0 2 0 1 0
0
4 4 4 4
, , ,p p p p
:
Ta có:

( ) ( )
( )
( )
( )
( )
1 2 3 3
1 2 3 3
2 3
3 0 4 0
3 0 3
4 4 4
3
3 4
3
. .( )
. .( )
.

0
)
0
0
(
c s
r
h h
s s h h
c h h
p R r p
h h
β β β
β β β
β β
+ +
 
 
= + +
 
 
− +
 
 
= ⇒ = × + =
 


+


+




( ) ( )
( )
( )
( )
1 2 3
2 0 3 0
2 0 2
3
3
1 2 3 3
3
4 2 4
2 3
3
. .(0 )
. .( )
(
0
0 . )
r
c s h h
s s h h
c h
p R r p
h

β β β
β β β
β β
+ +
 
 
+ +
 
 
− + +

 
 
= ⇒ = × + =
 



 
Page 4
Robot Công Nghiệp GVHD: TS Nguyễn Sĩ Dũng

( ) ( )
( ) ( )
( ) ( )
( ) ( )
1 2 2 2 3 3
1 2 2 2 3 3
2 2
1 0 2 0

2 2 2 3 3
1 0 1
2 2 2 4 1 2 4
. . .
. . .
. .
.
.
0
c h c s h h
s h c s
h c
r h s p h h
h s c
R p
h
r
h
β β β β
β β β β
β
β
β
β β
 

 
+ +

= ⇒ = ×

+ 
 
 
+ + + 
 
 
 
− + +
 
 
+ =
 
 
 

( )
( ) ( )
( ) ( )
( ) ( )
1 2 2 2 3 3
1 2
1 0
0 0 0
1 4 2 2 3 3
1 2 2 2
1
1
3 3
4
. . .

. . .
. .
0
0
c h c s h h
s h c s h h
h h s c h
r p r p
h
h
β β β β
β β β β
β β β
 
+ + + 
 
 
+ + + 
 

 
 

 
+ − + +
 

= ⇒ =

+ =

 
 
 
 Tính J
i
:

( ) ( )
( ) ( )
1 2 2 2 3 3
1 2 2 2 3 3
1
. . .
. . .
0
0
0
1
s h c s h h
c h c s h h
J
β β β β
β β β β
 
− + + + 
 
 
+ + + 
 
 

 
 
=
 
 
 
 
 


( ) ( )
( ) ( )
( ) ( )
1 2 2 2 3 3
1 2 2 2 3 3
2 2 2 3 3
2
. . .
. . .
. .
0
0
1
c h s c h h
s h s c h h
h c s h h
J
β β β β
β β β β
β β β

 
− − + + 
 
 
− − + + 
 
 
 
+ + +
 
=
 
 
 
 
 


( )
( )
( ) ( )
1 2 3 3
2 3 3
2 3 3
3
. .( )
.( )
.
0
0

1
c c h h
c h h
s h h
J
β β β
β β
β β
 + + 
 
+ +
 
 
+ +
=
 
 
 
 
 
 


( )
( )
( )
1 2 3
1 2 3
2 3
4

.
.
0
0
0
c s
s s
c
J
β β β
β β β
β β
 + 
 
+
 
 
− +
=
 
 
 
 
 
 

Page 5
Robot Công Nghiệp GVHD: TS Nguyễn Sĩ Dũng
Vậy:
( ) ( )

( ) ( )
( ) ( )
( ) ( )
( ) ( )
( ) ( )
( )
0
0
0
0
0
0
1 2 2 2 3 3 1 2 2 2 3 3
1 2 3 3
2 3
1 2 2 2 3 3 1 2 2 2 3 3
2 2 2 3 3

. . . . . .
. .
.
. . . . . .
0 . .
0 0
0 0
1 1
x
y
z
x

y
z
v
v
v
x
s h c s h h c h s c h h
c c h h
c h
c h c s h h s h s c h h
h c s h h
ω
ω
ω
β β β β β β β β
β β β
β β
β β β β β β β β
β β β
 
 
 
 
 
=
 
 
 
 
 

 
− + + + − − + +
   
+ +
   
+
+ + + − − + +   
   
+ + +
=
&
( )
( ) ( )
( )
( )
( )
1 2 3
3 1 2 3
2 3 3 2 3
1
2
3
.
.
.
0 0
0 0
1 0
c s
h s s

s h h c
d
β β β
β β β
β β β β
β
β
β
 
 
 
 
 
 
 
 
 
 
 
 
 
×
 
 
 
 
+
+ +
+ + − +
&

&
&
&
III. PHƯƠNG TRÌNH ĐỘNG LỰC HỌC CHO TAY MÁY:
1. Xác định vector vị trí tâm khối lượng C
i
:
Ta có:

1
0 0 0
0 1
1
0
. 0
c
p R r
h
 
 
= =
 
 
 


( ) ( )
( )
2 2
1 2 2

1 0 2 0
0 1
1 2 1 2 2
2 2
.
. .0
.
c c
c c h
p R r p s s h
s h
β β
β β
β
+
 
 
= = +



=

+




( )
2 2

1 2 2
1 0
0 0
1 1 2 2
1 2 2
.
.
.
c c
c c h
p r p s s h
h s h
β β
β β
β
=
+
 
 
= + +
 
 
+
 
Page 6
Robot Công Nghiệp GVHD: TS Nguyễn Sĩ Dũng

( ) ( )
( )
3 3

2 0 3 0
0 2
2 3
0
0
. 0
0
c c
p R r p
 
 
= + =


 
=




( ) ( )
3 3
1 2 2
1 0 2 0
0 1
1 2 1 2 2
2 2
. .
. . .
.

c c
c c h
p R r p s c h
s h
β β
β β
β
 
 
= =
 
 
 
+


( )
3 3
1 2 2
1 0
0 0 0
0 1 1 2 2
1 2 2
. .
. . .
.
c c
c c h
p R r p s c h
h s h

β β
β β
β
 
 
= + =
 
 
+
 


( ) ( )
( )
( ) ( )
( ) ( )
( ) ( )
4 4
1 2 3 3
3 0 4 0
0 3
3 4 1 2 3 3
2 3 3
. .
. .
.
0 .
c c
c s h h
p R r p s s h h

c h h
β β β
β β β
β β
+ +
 
 
= = + +
 
 
− + +
 
+ =


( ) ( )
( ) ( )
( ) ( )
( ) ( )
4 4
1 2 3 3
2 0 3 0
0 2
2 3 1 2 3 3
2 3 3
. .
. . .
.
c c
c s h h

p R r p s s h h
c h h
β β β
β β β
β β
+ +
 
 
= + = + +
 
 
− + +
 


( ) ( )
( ) ( )
( ) ( )
( ) ( )
4 4
1 2 2 2 3 3
1 0 2 0
0 1
1 2 1 2 2 2 3 3
2 2 2 3 3
.[ . . ]
. .[ . . ]
. .
c c
c c h s h h

p R r p s c h s h h
s h c h h
β β β β
β β β β
β β β
+ + +
 
 
= + = + + +
 
 
− + +
 


( )
( ) ( )
( ) ( )
( ) ( )
4 4
1 2 2 2 3 3
1 0
0 0 0
0 1 1 2 2 2 3 3
2 2 2 3 3 1
.[ . . ]
. .[ . . ]
. .
c c
c c h s h h

p R r p s c h s h h
s h c h h h
β β β β
β β β β
β β β
+ + +
 
 
= + = + + +
 
 
− + + +
 

2.Tính ma trận quán tính:

2
1 1 1
0 1 0 2
1 1 1 1 1 1 1 1
. 0
1
. .( ) . 0
48
0 0 1
T
s s c
I R I R m s c c
β β β
β β β

 

 
= = −
 
 
 

1 1 2 2
0 2 0
2 2 2 2 2
0 0 . .( )
1
. .( ) 0 1 0
75
0 0 1
T
s c c s
I R I R m
β β β β

 
 
= =
 
 
 

Page 7
Robot Công Nghiệp GVHD: TS Nguyễn Sĩ Dũng


33
0 3 0
3 3
(I3 . . )
T
R I RI = =

( ) ( ) ( ) ( )
( ) ( ) ( ) ( )
( ) ( ) ( ) ( ) ( )
2 2 2 2
1 1 2 3 1 1 2 3 1 2 3 2 3
2 2 2 2
3 1 1 2 3 1 1 2 3 1 2 3 2 3
2
2 3 1 2 3 2 3 1 2 3 2 3
. . . . .
3
. . . . .
400
. . . .
s c s s c c c s c
m s c c c s s s s c
c c s c s s c
β β β β β β β β β β β β β
β β β β β β β β β β β β β
β β β β β β β β β β β β
 
+ + − + − + +

 
= − + − + + − + +
 
 
− + + − + + +
 

0 4 0
4 4 44
. .( )
T
R RI I=

( ) ( ) ( ) ( )
( ) ( ) ( ) ( )
( ) ( ) ( ) ( ) ( )
2 2 2
1 2 3 1 1 2 3 1 2 3 2 3
2 2 2
4 1 1 2 3 1 2 3 1 2 3 2 3
2
2 3 1 2 3 2 3 1 2 3 2 3
. . . . .
1
. . . . .
300
. . . .
c s s c s c s c
m s c s s s s s c
c c s c s s c

β β β β β β β β β β β β
β β β β β β β β β β β β
β β β β β β β β β β β β
 
+ + + +
 
= + + + +
 
 
+ + + + +
 
3. Tính vận tốc khối tâm của từng khâu:
3.1.Tính ma trận Jacobi của từng khối tâm:
 Khâu 1: Vì là khớp quay nên ta có:

( )
1
1
1
0 0
0
0
0
c
v c
J z p
 
 
= × =
 

 
 


( )
1
1
0
0
0
1
c
J z
ω
 
 
= =
 
 
 

Suy ra:
1
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
1 0 0 0
c

J
 
 
 
 
=
 
 
 
 
 

 Khâu 2: Vì là khớp quay nên ta có:
Page 8
Robot Công Nghiệp GVHD: TS Nguyễn Sĩ Dũng

( )
2
2
1 2 2
1
0 0
1 2 2
.
.
0
c
v c
s c h
J z p c c h

β β
β β
− −
 
 
= × = +
 
 
 


( )
2
0
1
0
0
1
c
J z
ω
 
 
= =
 
 
 


( ) ( ) ( )

2
2
1 2 2
2 1 0 1 0
1 2 2
2 2 1 1
. .
. .
. .( )
c
v c
c s h
J z p s s h
c h s c
β β
β β
β β β

 
 
= × = −
 
 
+
 


( ) ( )
2
1

2
1
1 0
0
c
s
J z c
ω
β
β
 
 
= = −
 
 
 

Suy ra:
2
1 2 2 1 2 2
1 2 2 1 2 2
2 2 1 1
1
1
. . . 0 0
. . . 0 0
0 . .( ) 0 0
0 0 0
0 0 0
1 0 0 0

c
s c h c s h
c c h s s h
c h s c
J
s
c
β β β β
β β β β
β β β
β
β
− − −
 
 
+ −
 
 
+
=
 
 
 

 
 

 Khâu 3: Vì là khớp quay nên ta có:

( ) ( ) ( )

3
3
3 2 0 2 0
0
0
0
c
v c
J z p
 
 
= × =
 
 
 


( ) ( )
3
1
3
1
2 0
0
c
s
J z c
ω
β
β

 
 
= = −
 
 
 


( ) ( ) ( )
3
3
1 2 2
2 1 0 1 0
1 2 2
2 2
. .
. .
.
c
v c
c s h
J z p s s h
c h
β β
β β
β

 
 
= × = −

 
 
 

Page 9
Robot Công Nghiệp GVHD: TS Nguyễn Sĩ Dũng

( ) ( )
3
1
2 1 0
1
0
c
s
J z c
ω
β
β
 
 
= = −
 
 
 


( )
3
3

1 1 2
1
0 0
1 2 2
. .
. .
0
c
v c
s c h
J z p c c h
β β
β β

 
 
= × = −
 
 
 


( )
3
1
0
0
0
1
c

J z
ω
 
 
= =
 
 
 

Suy ra:
3
1 1 2 1 2 2
1 2 2 1 2 2
2 2
1 1
1 1
. . . . 0 0
. . . . 0 0
0 . 0 0
0 0
0 0
1 0 0 0
c
s c h c s h
c c h s s h
c h
J
s s
c c
β β β β

β β β β
β
β β
β β
− −
 
 

 
 
=
 
 
 
− −
 
 

 Khâu 4: Vì là khớp trượt nên ta có:

( ) ( ) ( )
4
4
3 04 3 0
0
0
0
v c
c
J z p

 
 
= × =
 
 
 


( ) ( )
( )
( )
( )
4
1 2 3
4 3 0
1 2 3
2 3
.
.
c
c s
J z s s
c
ω
β β β
β β β
β β
+
 
 

= = +
 
 
− +
 


( ) ( ) ( )
( ) ( )
( ) ( )
4
4
1 1 2 3 3
3 2 0 2 0
1 1 2 3 3
. . .
. . .
0
c
v c
s c s h h
J z p s c s h h
β β β β
β β β β
+ + 
 
= × = − + +
 
 
 



( ) ( )
4
1
1
3 2 0
0
c
s
J z c
ω
β
β
 
 
= = −
 
 
 

Page 10
Robot Công Nghiệp GVHD: TS Nguyễn Sĩ Dũng

( ) ( ) ( )
( ) ( )
4
4
2 1 0 1 0
2 3 3 1 1 1 1

0
0
. . . .( )
c
v c
J z p
s h h s c c s
β β β β β β
 
 
= × =
 
 
+ + −
 


( ) ( )
4
1
1
2 1 0
0
c
s
J z c
ω
β
β
 

 
= = −
 
 
 


( )
( ) ( )
( ) ( )
4
4
1 2 2 2 3 3
0
1 2 2 2 3 3
1
0
.[ . . ]
.[ . . ]
0
v c
c
s c h s h h
J z p c c h s h h
β β β β
β β β β
− + + +
 
 
= × = + + +

 
 
 


( )
4
1
0
0
0
1
c
J z
ω
 
 
= =
 
 
 

Suy ra:
4
c
J
=
( ) ( ) ( ) ( )
( ) ( ) ( ) ( )
( ) ( ) ( )

( )
( )
( )
1 2 2 2 3 3 1 1 2 3 3
1 2 2 2 3 3 1 1 2 3 3
2 3 3 1 1 1 1
1 1 1 2 3
1 1 1 2 3
2 3
. . 0 . . . 0
. . 0 . . . 0
0 . . . . 0 0
0 .
0 .
1 0 0
s c h s h h s c s h h
c c h s h h s c s h h
s h h s c c s
s s c s
c c s s
c
β β β β β β β β
β β β β β β β β
β β β β β β
β β β β β
β β β β β
β β
 
− + + + + + 
 

 
+ + + − + +
 
 
+ + −
 
 
+
 
− − +
 
 

 

+

 
3.2. Vận tốc khối tâm cuả cả cơ cấu tay máy 4 khâu:
• Khâu 1:
1
0
1
0
1
0
1
1
0
1

0
1
0
1 1
2 2
3 3
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
1 0 0 0
x
y
z
x
y
z
c
c
c
c
c
c
c
V
V
V
J
d d

β β
β β
ω
β β
ω
ω
 
 
 
 
 
 
 
 
×
 
 

 
 
   
 
   
= =
 
   
 
   
 
 


 

 
   
 
 
 



& &
& &
& &
& &
• Khâu 2:
Page 11
Robot Công Nghiệp GVHD: TS Nguyễn Sĩ Dũng
2
0
2
0
2
0
2
2
0
2
0
2

0
1 2 2 1 2 2
1 2 2 1 2 2
2 2 1 1
1
2 2
3
1
1
1
3
. . . 0 0
. . . 0 0
0 . .( ) 0 0
0 0 0
0 0 0
1 0 0 0
x
y
z
x
y
z
c
c
c
c
c
c
c

s c h c s
V
V
V
J
h
c c h s s h
c h s
s
c
d
c
d
β β β β
β β β β
β β
β β
β β
ω
β β
ω
ω
β
β
β
 
 
 
  
 

  
 
 
= =

− − −
 
 
+ −
 
 
+
×


 
 
 
 
 
  
 
 




 

 


 
 
& &
& &
& &
& &


 
 
 
 

• Khâu 3:
3
0
3
0
3
0
3
3
0
3
0
3
0
1 1
2 2

3 3
1 1 2 1 2 2
1 2 2 1 2 2
2 2
1 1
1 1
. . . . 0 0
. . . . 0 0
0 . 0 0
0 0
0 0
1 0 0 0
x
y
z
x
y
z
c
c
c
c
c
c
c
s c h c s h
c c h s s
V
V
V

J
h
c h
s s
c c
d d
β β β β
β β β β
β
β β
β β
β β
β β
ω
β β
ω
ω
 
 
 
   
 
  
 
  
= =
 
  
 
  

 
  
  
 

− −
 
 

 
 
×
 
 
 
− −
 
 

 
 
& &
& &
& &
& &







• Khâu 4:
4
0
4
0
4
0
4
4
0
4
0
4
0
1
2
3
x
y
z
x
y
z
c
c
c
c
c
c

c
V
V
V
J
d
β
β
ω
β
ω
ω
 
 
 
 
 
 
 
 
=
 
 
 
 
 
 
 
 
 

 
 
&
&
&
&
( ) ( ) ( ) ( )
( ) ( ) ( ) ( )
( ) ( ) ( )
( )
( )
( )
1 2 2 2 3 3 1 1 2 3 3
1 2 2 2 3
1
2
3 1 1 2 3 3
2 3 3 1 1 1 1
1 1 1 2 3
1 1 1 2 3
2 3
. . 0 . . . 0
. . 0 . . . 0
0 . . . . 0 0
0 .
0 .
1 0 0
s c h s h h s c s h h
c c h s h h s c s h h
s h h s c c s

s s c s
c c s s
c
β β β β β β β β
β β β β β β β β
β β β β β β
β β β β
β
β
β β β β β
β β
β
β
=
 
− + + + + +
 
 
 
+ + + − + +
 
 
+ + −
 
 
+
 
− − +
 
 

− +
 
×
 
 
&
&
&
3
d
 
 
 
 
 
 
 
&
Page 12
Robot Công Nghiệp GVHD: TS Nguyễn Sĩ Dũng
4.Tính ma trận quán tính cơ cấu chấp hành:
( ) ( )
1
. . . .
i i i i
T T
v v
c i c
n
i

i
c c
M J m J J I J
ω ω
=
 
= +
 
 


( ) ( ) ( ) ( )
( ) ( )
( ) ( )
[ ]
1 1 1 1 2 2 2 2
3 3 3 3 4 4 4 4
1 1 2 2
3 3 4 4
. . . . . . . .
. . . . . . . .
T T T T
v v v v
c c c c c c c c
T T
T T
v v v v
c c c c c c c c
J m J J I J J m J J I J
J m J

M
A B C D
J I J J m J J I J
ω ω ω ω
ω ω ω ω
⇒ = +
+ +
=
+ +
+ +
=
Với:

( ) ( )
1 1 1 1
1
1 1
1
. 0 0 0
48
. . . . 0 0 0 0
0 0 0 0
T T
v v
c c c c
m
J m J J IA J
ω ω
 
 

 
+ =
 
 
 
 
=


( ) ( )
2 2 2 2
2 2
. . . .
T T
v v
c c c c
J m J J I JB
ω ω
+=


2 2 2
1 2 2 1 2 2 2 2 2 2 1 1
2 2 2 2 2 2
2 2 2 2 1 1 2 2 2 2 1 1
2
1 1
1 1 2 2
2
( . ) ( . ) . . . .( ) 0 0

. . . .( ) . . .( ) 0 0
0( ) 0 0 0
0 0 0 0
0 0 0 1
0 0 . .( )
0 0 0 1
1
. . 0 1 0
0 0 0 0
75
0 0 1
0 0 0 0
s c h c c h h c s h c s
h c s h c s s h h c c s
m
c s
c s c s
m
β β β β β β β β
β β β β β β β β
β β
β β β β
 
+ + + −
 
− + +
 
=
 


 
 
 

 
 
+
 
 
 
0 0 0 0
. 0 0 0 0
1 1 0 0
   
   
   
   
   

( ) ( )
3 3 3 3
3 3
. . . .
T T
v v
c c c c
C J m J J I J
ω ω
= +


Page 13
Robot Công Nghiệp GVHD: TS Nguyễn Sĩ Dũng

( )
1 1 2 1 2 2
1 1 2 1 2 2
1 2 2 1 2 2 2 2
3 1 2 2 1 2 2
2 2
1 1
3
1 1
2 2 2
1 1 2 3
. . . . 0
. . . . 0 0
. . . . .
. . . . 0 0
0 0 0
0 . 0 0
0 0 0
0 0 1
0
3
0
400
0 0 0
.
s c h c c h
s c h c s h

c s h s s h c h
m c c h s s h
c h
s c
m
s c
s c s
β β β β
β β β β
β β β β β
β β β β
β
β β
β β
β β β β

 
− −
 
 
− −
 
 
= × × −
 
 
 
 
 
 

 
 

 
+ ×
 

 
 
+ +
×

( ) ( ) ( )
( ) ( ) ( ) ( )
( ) ( ) ( ) ( ) ( )
2
1 1 2 3 1 2 3 2 3
2 2 2 2
1 1 2 3 1 1 2 3 1 2 3 2 3
2
2 3 1 2 3 2 3 1 2 3 2 3
1 1
1 1
. . . .
. . . . . .
. . . .
0 0
0 0
1 0 0 0
s c c c s c

s c c c s s s s c
c c s c s s c
s s
c c
β β β β β β β β β
β β β β β β β β β β β β β
β β β β β β β β β β β β
β β
β β
 
+ − + +
 
− + − + + − + +
 
 
− + + − + + +
 
 
 
× − −
 
 
 
( ) ( )
4 4 4 4
4 4
. . . .
T T
v v
c c c c

J m J J I JD
ω ω
+=

=
( ) ( ) ( ) ( )
( ) ( )
( ) ( ) ( ) ( )
1 2 2 2 3 3 1 2 2 2 3 3
2 3 3 1 1 1 1
1 1 2 3 3 1 1 2 3 3
. . . . . 0
0 0 . . . .( )
. . . . . . 0
0 0 0
.s c h s h h c c h s h h
s h h s c c s
s c s h h s c s h h
β β β β β β β β
β β β β β β
β β β β β β β β
 
− + + + + + + 
 
 
+ + −
 
 
+ + − + +
 

 
 

 
( ) ( ) ( ) ( )
( ) ( ) ( ) ( )
( ) ( )
1 2 2 2 3 3 1 1 2 3 3
4 1 2 2 2 3 3 1 1 2 3 3
2 3 3 1 1 1 1
. . . 0 . . . 0
. . . 0 . . . 0
0 . . . .( ) 0 0
s c h s h h s c s h h
m c c h s h h s c s h h
s h h s c c s
β β β β β β β β
β β β β β β β β
β β β β β β
 
− + + + + +
 
 
 
× + + + − + + 
 
 
 
+ + −
 

 
( ) ( ) ( )
1 1
4
1 1
1 2 3 1 2 3 2 3
0 0 1
0
1
0
300
. .
s c
m
s c
c s s s c
β β
β β
β β β β β β β β
 
 

 
+ ×
 

 
+ + − +
 


Page 14
Robot Công Nghiệp GVHD: TS Nguyễn Sĩ Dũng
( ) ( ) ( ) ( )
( ) ( ) ( ) ( )
( ) ( ) ( ) ( ) ( )
2 2 2
1 2 3 1 1 2 3 1 2 3 2 3
2 2 2
1 1 2 3 1 2 3 1 2 3 2 3
2
2 3 1 2 3 2 3 1 2 3 2 3
. . . . .
. . . . .
. . . .
c s s c s c s c
s c s s s s s c
c c s c s s c
β β β β β β β β β β β β
β β β β β β β β β β β β
β β β β β β β β β β β β
 
+ + − + +
 
× + + − + +
 
 
− + + − + + +
 
5. Tính các vector V
i

và G
i
:
5.1. TínhV
i
:
j
n
j
n
k
i
i
jk
k
ij
i
qq
q
M
q
M
V

= =

∑∑






=
1 1
)
2
1
(
V
1
=
V
2
=
V
3
=
V
4
=
5.2. Tính G
i
:
( )
2
1
n
T i
i j v
j

G m g J
=
= −

Với
0
0
g g
 
 
=
 
 
 
Suy ra :
6.Phương trình chuyển động Lagrang:
Page 15
Robot Công Nghiệp GVHD: TS Nguyễn Sĩ Dũng
ij
1
i
i i
n
j i i i
j
d L L
Q
dt q q
M q V G Q
=

 
∂ ∂
− =
 ÷
∂ ∂
 
⇔ + + =

& &
&
Lực suy rộng :
T
c r
Q J F f
τ
= + −
Bỏ qua lực ma sát tại các khớp :
0
r
f
=
Vậy lực suy rộng sẽ là :
1
1 1
2 2
2
3 3
3
F
Q

Q
F
Q
F
τ
τ
τ
 
   
 
   
= =
 
   
 
   
   
 
 
Thay vào phương trình động lực học Lagrang, ta được:
ij
1
n
i
j i i i i
j
M q V G Q F
τ
=
+ + = = +


&
IV. CHỌN ĐỘNG CƠ CHO ROBOT:
Do yêu cầu điều khiển chính xác điểm tác động cuối của robot ta chọn động cơ cho
các khớp quay là động cơ servo.
Để phù hợp với tay máy ta chọn động cơ servo có momen xoắn 130 oz-in, thời gian
vận chuyển (vận tốc) 0.21s/60
o
.
1. Giới thiệu về động cơ servo:
Động cơ servo đươc thiết kế cho những hệ thống hồi tiếp vòng kín. Tính hiệu
ra của động cơ được nối với một mạch điều khiển. Khi động cơ quay vị trí và vận tốc
sẽ được hồi tiếp về mạch điều khiển này. Nếu có bất cứ lý do gì cản trở chuyển động
quay của động cơ, cơ cấu hồi tiếp sẽ nhận tín hiệu ra chưa đạt được vị trí mong
muốn. Mạch điều khiển tiếp tục chỉnh sai lệch cho động cơ đạt được điểm chính xác.
Page 16
Robot Công Nghiệp GVHD: TS Nguyễn Sĩ Dũng
Các động cơ servo điều khiển bằng liên lạc vô tuyến được gọi là động cơ
servo R/C (radio-controlled). Trong thực tế bản thân động cơ servo không phải được
điều khiển bằng vô tuyến, nó chỉ nối với máy thu vô tuyến trên máy bay, xe hơi.
Động cơ Servo nhận tín hiệu từ máy thu này. Như vậy có nghĩa là ta không cần phải
điều khiển robot bằng tín hiệu vô tuyến bằng cách sử dụng một động cơ servo, trừ
khi ta muốn thế. Ta có thể điều khiển động cơ servo bằng máy tính, một mạch vi xử
lý hay thậm chí một mạch điện tử cơ bản dùng IC555.
Mô hình động cơ Servo
Bên trong của một động cơ servo R/C bao gồm một động cơ, một chuỗi các
bánh răng giảm tốc, một mạch điều khiển và một vôn kế.
Động cơ và vôn kế nối với mạch điều khiển, tạo thành mạch hồi tiếp kín. Cả
mạch điều khiển và động cơ đều được cấp nguồn DC (thường từ 4,8 – 7,2V).
Để quay động cơ, tín hiệu đó được gởi tới mạch điều khiển. Tín hiệu này khởi

động động cơ thông qua chuỗi bánh răng, nối với trục vôn kế. Vị trí của trục vôn kế
cho biết vị trí trục ra của servo.Khi vôn kế đạt được vị trí mong muốn mạch điều
khiển sẽ tắt động cơ.
Động cơ servo được thiết kế quay có giới hạn chứ không phải quay liên tục
như động cơ DC hay động cơ bước. Mặc dù ta có thể chỉnh động cơ servo quay liên
tục nhưng công dụng chính của động cơ servo là quay đạt một góc chính xác trong
khoảng từ 90
o
– 180
o
. Việc điều khiển này được ứng dụng để lái robot, di chuyển các
tay máy lên xuống, quay một cảm biến để quét khắp phòng.
2. Servo và điều biến độ xung:
Page 17
Robot Công Nghiệp GVHD: TS Nguyễn Sĩ Dũng
Trục của động cơ servo R/C được định vị nhờ vào kỹ thuật gọi là điều biến độ
xung (PWM).Trong hệ thống này, Servo là đáp ứng của một dãy các xung ổn
định.Cụ thể hơn mạch điều khiển là đáp ứng cũa một tín hiệu số có các xung biến đổi
từ 1 – 2ms.Các xung này được gửi đi 50 lần/giây.Chú ý là không phải số xung trong
một giây điều khiển servo mà là chiều dài các xung.Servo đòi hỏi khoảng từ 30 – 60
xung/giây.Nếu số này quá thấp, độ chính xác và công suất duy trì servo sẽ giảm.
Với độ dài xung 1ms, servo được điều chỉnh quay theo một chiều:
Với độ dài xung 2ms, servo quay theo chiều ngược lại. Kỹ thuật này còn được
gọi là tỷ lệ số - chuyển động của servo tỉ lệ với tín hiệu số điều khiển.
Công suất cung cấp cho động cơ bên trong servo cũng tỹ lệ với độ lệch giữa vị
trí hiện tại của trục ra với vị trí nó cần tới.Nếu servo ở gần vị trí đích, động cơ được
chuyển động với tốc độ thấp.Điều này đảm bảo động cơ không vượt quá điểm cần
đến.Nhưng nếu servo ở xa vị trí đích nó sẽ được chuyền động với vận tốc tối đa để
đến đích nhanh.Khi trục ra đến vị trí mong muốn, động cơ giảm tốc.Quá trình phức
tạp này diễn ra trong một thời gian ngắn – một servo trung bình có thể quay 60

o
trong
vòng ¼ - ½ giây.
Vì độ dài xung có thể thay đổi tuỳ theo hãng chế tạo nên ta phải trọn servo và
máy thu vô tuyến thuộc cùng một hãng đẻ đảm bảo độ tương thích. Đối với robot ta
phải làm một vài thí nghiệm để xác định xung tối ưu.
3. Vai trò của vôn kế:
Page 18
Robot Công Nghiệp GVHD: TS Nguyễn Sĩ Dũng
Vôn kế trong Servo giữ vai trò chính trong việc định vị trí của các trục ra. Vôn
kế được gắn vào trục ra (trong một vài servo,Vôn kế chính là trục ra).Bằng cách này,
vị trí của vôn kế phản ánh chính xác vị trí của truc ra của servo.Vôn kế hoạt động
nhờ cung cấp một điện áp biến thiên cho mach điều khiển.Khi cần chạy bên trong
Vôn kế điều khiển, điện thế sẽ thay đổi.
Mạch điều khiển trong servo so sánh này với độ dài các xung số đưa vào và
phát “tín hiệu sai số” nếu điện thế không đúng.Tín hiệu sai số này tỷ lệ với độ chênh
lệch giữa vị trí của vôn kế và độ dài tín hiệu vào.Mạch điều khiển sẽ kết hợp tín hiệu
sai số này để quay động cơ.Khi điện thế của vôn kế và độ dài các xung số bằng nhau,
tín hiệu sai số được loại bỏ và động cơ ngừng.
4. Các giới hạn quay:
Các servo khác nhau ở góc qua được với cùng tín hiệu 1 – 2ms (hoặc
bất kỳ) được cung cấp.Các servo chuẩn được thiết kế quay tới và lui từ 90
o
-180
o
khi
được cung cấp toàn bộ chiều dài xung.
Nếu ta cố điều khiển các servo vượt quá giới hạn cơ học của nó, trục ra
của động cơ sẽ đụng vật cản bên trong, dẫn đến các bánh răng bị mài mòn hay bị rơ.
Hiên tượng này kéo dài hơn 10 giây sẽ làm cho bánh răng của động cơ bị phá huỷ.

5. Hệ thống truyền động bánh răng và truyền công suất:
Động cơ bên trong servo R/C quay khoảng vài ngàn vòng /phút. Tốc
độ này quá nhanh để có thể dùng trực tiếp trên mô hình máy bay, xe hơi hay robot.
Tất cả các servo đều có một hệ thống bánh răng để giảm vận tốc ra của động cơ còn
Page 19
Robot Công Nghiệp GVHD: TS Nguyễn Sĩ Dũng
khoảng 50-100 vòng/phút. Các bánh răng của servo có thể làm plastic, nylon hay kim
loại (thường là đồng thau hay nhôm).
Bánh răng kim loại có tuổi thọ cao nhưng giá thành cũng cao. Các bánh
răng thay thế luôn có sẵn.Khi một hoặc vài bánh răng bị hư servo không khớp và ta
phải thay bánh răng. Trong một vài trường hợp ta có thể nâng cấp bánh răng plastic
thành bánh răng kim loại.
Bên cạch các bánh răng dẫn động trục ra của động cơ cũng thường bị
mòn và xước, trong các servo rẻ nhất trục này được đỡ bằng miếng đệm plastic.
Miếng đệm này rất dễ mất tác dụng nếu động cơ chạy nhiều.Servo sử dụng vòng bi
có tuổi thọ cao.
6. Mạch điều khiển servo:
Không giống động cơ DC ta chỉ cần lắp pin vào là chạy, động cơ DC đòi hỏi
một mạch điện tử chính xác để quay trục của nó. Có thể một mạch điện tử sẽ làm
việc sử dụng servo phức tạp hơn ở một mức độ nào đó nhưng thực ra mạch điện tử
này rất đơn giản. Nếu ta muốn điều khiển servo bằng máy tính hoặc bằng bộvi xử lý
thỉ chỉ cần một vài dòng lệnh là đủ.
Một động cơ DC điển hình cần các transitor công suất, mosfet hay replay nếu
muốn kết hợp với máy tính, còn servo có thể gắn trực tiếp với máy tính hoặc bộ vi xử
lý mà không cần một linh kiện điện tử nào cả. Tất cả yếu tố cần thiết để điều khiển
công suất đều được quản lý bởi mạch điều khiển để tránh rắc rối.Đây là lợi ích chủ
yếu khi sử dụng servo cho các robot điều khiển bằng máy tính.
a. Điều khiển servo IC định thì 555:
Ta có thể không cần đến máy tính để điều khiển servo. Một IC555 có thể cung
cấp các xung cần thiết cho servo.

Page 20
Robot Công Nghiệp GVHD: TS Nguyễn Sĩ Dũng
Một phương pháp phổ biến dùng IC555 để điều khiền servo
Khi hoạt động, IC555 sinh ra một tín hiệu xung có chu kỳ nhiệm vụ khác nhau
để điều khiển hoạt động của servo. Chỉnh Vôn kế để định vị servo.Vì IC555 có thể
dễ dàng tạo xung rất dài và rất ngắn nên servo có thể hoạt động ngoài vị trí biên
thông thường. Khi servo gặp vật cản ta phải ngắt nguồn lập tức, nếu không các bánh
răng bên trong động cơ bị trờn.
b. Dùng bộ xử lý chuyên nghiệp:
Các máy thu R/C được thiết kế tối đa với 8 servo. Máy thu nhận xung
số từ máy phát, bắt đầu bằng một xung dài đồng bộ, sau đó là các xung của 8 servo,
mỗi xung dành cho một servo. 8 xung cộng với xung đồng bộ mất khoảng 20 ms.
Điều này có nghĩa là dãy xung có thể lập lại 50 lần/giây.Ta gọi đó là tỷ lệ lặp lại
(refresh rate). Khi tỉ lệ này giảm các servo cập nhật không đủ nhanh và sẽ bị mất vị
trí.
Trừ khi mạch điện tử ta dung có thể cung cấp xung đồng thời cho nhiều servo
(đa nhiệm vụ - multi tasking), mạch điều khiển không thể cung cấp xung lập lại đủ
nhanh. Vì vậy ta có thể dùng bộ xử lý servo chuyên nghiệp. Bộ này có thể điều khiển
5, 8 động cơ hay nhiều hơn một cách độc lập, sẽ làm giảm bớt chương trình tổng
cộng của máy tính hay bộ vi xử lý mà ta đang dùng.
Ưu điểm của bộ vi xử lý servo chuyên nghiệp là ta có thể điều khiển đồng thời
nhiều servo ngay cả khi máy tính, bộ vi xử lý không đa nhiệm vụ.
Ví dụ: giả sử robot cần 24 servo, có thể là một robot dùng nhện 8 chân, mỗi
chân có 3 servo, mỗi servo điều khiể một bậc tự do của chân. Phương pháp ta sử
dụng là phân chia công việc cho 3 bộ xử lý servo, mỗi bộ có thể điểu khiển 8 servo.
Mỗi bộ chịu trách nhiệm xử lý cho một loại bậc tự do: một cho hoạt động quay của 8
chân, một cho “độ linh hoạt” của các chân và một cho độ quay đốt cuối của chân.
Page 21
Robot Công Nghiệp GVHD: TS Nguyễn Sĩ Dũng
Các bộ xử lý servo chuyên nghiệp phải được dung với máy tính hoặc bộ vi xử

lý vì chúng cần được cung cấp dữ liệu từ thời gian thực để điều khiển servo. Dữ liệu
này thường được gởi trong một công thức dữ liệu chuỗi. Một dãy các byte gởi từ máy
tính hay bộ vi xử lý được bộ vi xử lý servo giải mã, mà mỗi byte tương ứng với một
servo. Mỗi bộ xử lý servo điển hình có ghi chú các ứng dụng và các chương trình
mẫu của các máy tính và bộ vi xử lý thông dụng nhưng để đảm đảm bảo ta cần có
kiến thức về lập trình và truyền chuỗi.
c. Sử dụng lớn hơn 7,2 V:
Servo được thiết kế sử dụng với bộ pin R/C recharge, có điện thế từ 4,8 –
7,2V, phụ thuộc vào số pin sử dụng. Các servo cho phép khoảng điện thế vào khá
rộng, và bộ 4 pin AA 6V đã cung cấp đủ điện. Tuy nhiên khi pin hết điện thế giảm,
servo không còn nhanh như lúc đầu. Khi điện thế xuống 4 -4,5V servo thậm chí
không chạy.
Nhiều servo có thể chạy ở mức điện thế cao khoảng 12V mà ít gây hậu quả.
Tuy nhiên đa số các servo đều nóng lên ở mức điện thế 9-10V, chúng sẽ không hoạt
động lâu nếu chúng không đươc nghỉ hoặc làm nguội.
Trừ khi ta cần tăng momen xoắn hay tốc độ, tốt hơn là cung cấp mức điện thế
cho servo dưới 9V, tốt nhất là trong khoảng 4,8-7,2V. Ta cũng cần tham khảo thêm
bảng dữ kiệu của servo để xác định các yêu cầu điện thế đặc biệt khác.
d. Làm việc với dải chết và tránh dải chết:
Tất cả servo đều có dải chết. Dải chết của servo là thời gian sai lệch lớn nhất
giữa các tín hiệu điều khiển ngõ vào và tín hiệu tham chiếu nội sinh ra bởi vị trí của
Vôn kế. Nếu thời gian sai lệch nhỏ hơn dải chết 5-6ms, servo không cần phải điều
chỉnh động cơ để sửa sai lệch.
Nếu không có dải chết servo phải liên tục dò tới lui để tìm điểm tương thích
giữa tín hiệu vào và tín hiệu chiếu nội của nó. Dải chết cho phép servo giảm thiểu sự
dò tìm này và sẽ lấy lân cận điểm cần tìm mặc dù không chính xác lắm.
Page 22
Robot Công Nghiệp GVHD: TS Nguyễn Sĩ Dũng
Dải chết được thay đổi tùy theo servo và được coi như một thông số của
servo. Dải chết điển hình dài 5s. Nếu servo quay 180

o
trong khoảng 1000s thì dải
chết 5 chỉ chiếm 1/200. Ta không cần lưu ý tới ảnh hưởng của dải chết nếu mạch
điều khiển có độ phân giải thấp hơn dải chết.
Tuy nhiên nếu mạch điều khiển có đọ phân giải cao hơn dải chết, một sự thay
đổi nhỏ về giá trị độ rộng xung có thể không ảnh hưởng. Ví dụ nếu bộ vi xử lý có độ
phân giải là 2s và nếu servo có dải chết 5s thì sự thay đổi một trong hai giá trị- Tức là
2-4s trong bề rộng xung sẽ không ảnh hưởng lên servo.
Như vậy ta lên chọn servo có dải chết hẹp nếu ta cần độ chính xác và mạch
điều khiển hay môi trường lập trình có độ phân giải đủ lớn. Ngược lại ta cần lưu ý tới
dải chết.
e. Dải xung lớn hơn 1-2s:
Servo điển hình có tín hiệu đáp ứng từ 1-2s. Trong thực tế nhiều servo có thể
được cung cấp bởi xung ngắn hay dài hơn để tối đa hóa giới hạn quay. Dải 1-2s thật
ra có thể quay servo theo hai hướng nhưng không thể quay toàn bộ theo cả hai
hướng. Tuy nhiên ta không biết rõ giá trị lớn nhất, nhỏ nhất của servo cho đến khi ta
chạy thử.
Nếu ta chỉ cần quay servo đến giá trị max, hãy chọn lực mạch điều khiển. Bắt
đầu bằng cách thay đồi một lượng nhỏ bề rộng xung hơn 1s, có thể là 10s. Sau mỗi
lần thay đổi, dung chương trình điều khiển đẩy servo trở lại vị trí giữa (vị trí trung
hòa). Khi nghe thấy servo gặp vật cản bên trong ( tiếng lạch cạch), lúc đó ta đã tìm
được giá trị biên dưới của servo, lập lại quá trình để tìm giá trị biên trên. Có những
servo có cận dưới là 250s, cận trện là 2200s. Tuy nhiên các servo khác bị hạn chế đến
nỗi thậm chí chúng không thể hoạt động trong dải 1-2ms. Các giá trị biên này cũng
khác nhau đối với từng loại và từng nhãn hiệu khác nhau.
V. HÀM TRUYỀN CHUYỂN ĐỘNG CỦA MỖI KHỚP
Nội dung phần này là trình bày phương pháp xây dựng với trường hợp chuyển
động 1 bậc tự do, mỗi khớp thường được điều khiển bằng một hệ truyền động
riêng.Phổ biến nhất là dòng điện 1 chiều
Page 23

Robot Công Nghiệp GVHD: TS Nguyễn Sĩ Dũng
Xét sơ đồ chuyển động của động cơ điện một chiều với tín hiệu vào là phần
điện áp U
a
đặt vào phần ứng ,tín hiệu ra là góc quay θ
m
của động cơ,động cơ kiểu
kích thích từ độc lập.
Sơ đồ động cơ điện 1 chiều
Trục động cơ được nối với hộp giảm tốc rồi tới trục phụ tải , gọi n là tỉ số truyền, θ
L

là góc phụ tải .Ta có:
θ
L
(t)= n.θ
m
(t)
L
=n
m
(t)
L
=n.
m
(t)
Sơ đồ cơ điện cùng phụ tải
Page 24
Robot Công Nghiệp GVHD: TS Nguyễn Sĩ Dũng
Momen trên trục động cơ bằng tổng Momen để động cơ quay, cộng với momen phụ

Tải quy về trục động cơ.
Khi đó ta có phương trình M(t)=M
m
(t)+ M
L
*
(t)
Kí hiệu: J
m
: Momen quán tính của động cơ
J
L
: Momen quán tính phụ tải
Ta có:
( ) ( ) ( )
( ) ( ) ( )
m m m m m
L L L L L
M t J t f t
M t J t f t
θ θ
θ θ
= +
= +
&& &
&& &

Trong đó f
m
và f

L
hệ số cản của động cơ và của phụ tải
Công do phụ tải sinh ra tính trên trục phụ tải M
L
θ
L
phải quy về trục động cơ M
L
*
. Từ
đó ta tín được :
( )
( ) ( )
( )
( )
( ) ( ) ( )
( )
( )
( ) ( )
( ) ( ) ( )
*
* 2
2 2
L
L L
m
L L m L m
m L m m L m
m m
M t t

M t nM t
t
M t n J t f t
M t J n J t f n f t
M t J t f t
θ
θ
θ θ
θ θ
θ θ
= =
 
⇔ = +
 
⇔ = + +
⇒ = +
&& &
&& &
&& &

Trong đó:
J = J
m
+ n
2
J
L
: Momen quán tính tổng hiệu dụng
f = f
m

+ n
2
f
L
: hệ số ma sát tổng hiệu dụng
Monen trên trục động cơ phụ thuộc vào tuyến tính với cường độ dòng điện phần ứng
và không phụ thuộc vào góc quay và vận tốc góc ,ta có:
( ) ( )
a a
M t K i t=

Với: i
a
: cường độ dòng điện phần ứng
Page 25

×