Nghiên cứu hệ thống điều khiển vecto động cơ không đồng bộ dùng cảm biến tốc độ - Pdf 43

BỘ GIÁO DỤC VÀ ĐÀO TẠO
TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI
---------------------------------------

NGUYỄN HOÀNG MINH

NGHIÊN CỨU HỆ THỐNG
ĐIỀU KHIỂN VECTO ĐỘNG CƠ KHÔNG
ĐỒNG BỘ KHÔNG DÙNG CẢM BIẾN TỐC ĐỘ

LUẬN VĂN THẠC SĨ KĨ THUẬT

Hà Nội – Năm 2016


BỘ GIÁO DỤC VÀ ĐÀO TẠO
TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI
---------------------------------------

NGUYỄN HOÀNG MINH

NGHIÊN CỨU HỆ THỐNG
ĐIỀU KHIỂN VECTO ĐỘNG CƠ KHÔNG ĐỒNG
BỘ KHÔNG DÙNG CẢM BIẾN TỐC ĐỘ

Chuyên ngành : ĐIỀU KHIỂN VÀ TỰ ĐỘNG HÓA

LUẬN VĂN THẠC SĨ KĨ THUẬT

NGƯỜI HƯỚNG DẪN KHOA HỌC
TS. NGUYỄN MẠNH TIẾN

Một số quy ước
Chỉ số viết trên cao bên phải.
-

f: là đại lượng mô tả trên hệ tọa độ tựa theo từ thông (hệ tọa độ dq quay đồng
bộ với vecto từ thông).
s: đại lượng mô tả trên hệ tọa độ 𝛼𝛽 cố định với stator.
r: đại lượng mô tả trên hệ tọa độ cố định với roto.

Chỉ số phía dưới bên phải.
-

Chữ cái thứ nhất: d, q các thành phần thuộc hệ tọa độ dq.
𝛼𝛽 các thành phần thuộc hệ tọa độ 𝛼𝛽
Chứ cái thứ hai: s đại lượng mạch stator.
r: đại lượng mạch rotor
Lm: Hỗ cảm giữa rotor và stator
L𝛿s: Điện cảm tiêu tán phía cuộn dây stato
L𝛿r: Điện cảm tiêu tán phía cuộn dây roto đã quy đổi về stator
Ls = Lm + L𝛿s: Điện cảm stator
Lr = Lm + L𝛿r: Điện cảm roto
Ts = Ls/Rs : Hằng số thời gian mạch stato
Tr = Lr/Rr : Hằng số thời gian mạch roto
𝛿 =1−

𝐿2 𝑚

𝐿𝑠 .𝐿𝑟

: Hệ số tiêu tán tổng.

Hình 1.6. Mô hình động cơ trên hệ tọa độ dq

14

Hình 1.7. Sự tương tự giữa phương pháp điều khiển động cơ một chiều và điều

15

khiển vecto động cơ không đồng bộ
Hình 1.8. Đồ thị pha của phương pháp điều khiển vecto

17

Hình 1.9. Sơ đồ khối của hệ điều khiển vecto động cơ không đồng bộ

18

Hình 1.10. Xác định từ thông roto trên hệ tọa độ vuông góc (𝛼𝛽)

19

Hình 1.11. Đồ thị góc pha của phương pháp điều khiển vecto gián tiếp

19

Hình 1.12. Sơ đồ cấu trúc tính toán góc quay từ trường

21

Hình 1.13. Sơ đồ khối hệ thống điều khiển vecto động cơ không đồng bộ


37

Hình 3.1. Hệ thống phản hồi dạng chuẩn

39

Hình 3.2. Sơ đồ cấu trúc hệ thống mô hình chuẩn

43

Hình 3.4. Sơ đồ cấu trúc thuật toán nhận dạng tham số động cơ

46

Hình 3.5. Sơ đồ cấu trúc mô hình quan sát

47

Hình 3.6. Mô hình nhận dạng tốc độ với mô hình quan sát từ thông

47

Hình 4.1. Sơ đồ điều khiển vecto với tính toán tốc độ

48

Hình 4.2. Sơ đồ hệ thống mô phỏng điều khiển vecto động cơ không đồng bộ
dùng cảm biến tốc độ


Hình 4.10. Đồ thị sai lệch giữa tốc độ thực của động cơ và tốc độ tính toán

55

Hình 4.11. Đồ thị từ thông quan sát trong hệ thống điều khiển động cơ không

56

đồng bộ
v


Hình 4.12. Đồ thị sai lệch giữa từ thông thực của động cơ và từ thông tính toán

56

Hình 4.13. Đồ thị dòng điện 𝐼𝑑 trong hệ thống điều khiển động cơ không đồng bộ

57

Hình 4.15. Đồ thị dòng điện 𝐼𝑞 trong hệ thống điều khiển động cơ không đồng bộ

58

Hình 4.14. Đồ thị sai lệch giữa dòng 𝐼𝑑 thực của động cơ và dòng 𝐼𝑑 tính toán

57

Hình 4.16. Đồ thị sai lệch giữa dòng 𝐼𝑞 thực của động cơ và dòng 𝐼𝑞 tính toán


64

Hình 4.23. Đồ thị sai lệch dòng điện 𝐼𝑑 giữa hai hệ thống điều khiển vecto
dùng cảm biến tốc độ và không dùng cảm biến tốc độ

Hình 4.25. Đồ thị sai lệch dòng điện 𝐼𝑞 giữa hai hệ thống điều khiển vecto
dùng cảm biến tốc độ và không dùng cảm biến tốc độ

Hình 4.26. Đồ thị mômen động cơ trong hệ thống điều khiển vecto động cơ
không đồng bộ không dùng cảm biến tốc độ

vi

63

64
65


MỤC LỤC
Chương 1. TỔNG QUAN VỀ HỆ THỐNG ĐIỀU KHIỂN ĐỘNG CƠ KHÔNG ĐỒNG BỘ
TRONG HỆ KHÔNG GIAN VECTƠ ................................... Error! Bookmark not defined.
1.1. Tổng quan động cơ không đồng bộ dưới dạng không gian vectơError! Bookmark
not defined.
1.1.1. Vectơ không gian ..........................................Error! Bookmark not defined.
1.1.2. Chuyển hệ tọa độ cho các vectơ không gian Error! Bookmark not defined.
1.2. Mô hình toán học động cơ khồng đồng bộ ......... Error! Bookmark not defined.
1.2.1. Hệ phương trình cơ bản ................................ Error! Bookmark not defined.
1.2.2. Các phương trình và mô hình trạng thái liên tục của động cơ trên hệ tọa độ 𝛼𝛽
................................................................................Error! Bookmark not defined.

defined.
Chương 3. XÂY DỰNG THUẬT TOÁN XÁC ĐỊNH TỐC ĐỘ VÀ TỪ THÔNG ROTO..
.................................................................................................. Error! Bookmark not defined.
3.1. Khái niệm siêu ổn định của Popov...................... Error! Bookmark not defined.
3.1.1. Hệ thống phản hồi dạng chuẩn ..................... Error! Bookmark not defined.
3.1.2. Các định nghĩa ổn định Popov ...................... Error! Bookmark not defined.
3.1.3. Định lý ổn định Popov .................................. Error! Bookmark not defined.
3.2. Mô hình động cơ trong hệ tọa độ dq ................... Error! Bookmark not defined.
3.3. Thuật toán tính toán tốc độ.................................. Error! Bookmark not defined.
3.4. Mô hình quan sát .................................................Error! Bookmark not defined.
Chương 4. HỆ THỐNG ĐIỀU KHIỂN VECTO KHÔNG DÙNG CẢM BIẾN TỐC
ĐỘ..Error! Bookmark not defined.
4.1. Hệ thống điều khiển vecto không dùng cảm biến tốc độ .. Error! Bookmark not
defined.
4.2. Mô phỏng thuật toán tính toán tốc độ ................. Error! Bookmark not defined.
4.2.1. Các sơ đồ mô phỏng ......................................Error! Bookmark not defined.
4.2.2. Kết quả mô phỏng kiểm nghiệm thuật toán ước lượng tốc độError! Bookmark
not defined.
4.3. Kết quả mô phỏng hệ thống điều khiển vecto không dùng cảm biến tốc độError!
Bookmark not defined.
KẾT LUẬN ............................................................................. Error! Bookmark not defined.
TÀI LIỆU THAM KHẢO ...................................................... Error! Bookmark not defined.
PHỤ LỤC A ............................................................................ Error! Bookmark not defined.
viii


PHỤ LỤC B........................................................................... Error! Bookmark not defined.6

ix



𝚤̅𝑐𝑠 = 𝑖𝑐𝑠 𝑒 𝑗240 = 𝑖𝑐𝑠 𝑎2

Thiết lập một hệ toạ độ trên mặt cắt ngang của động cơ một hệ trục toạ độ phức
(trục thực α, trục ảo β), trục thực trùng với trục cuộn dây pha A (hình 1.2)

Hình 1.2. Thiết lập vectơ không gian từ các đại lượng 3 pha
Vậy vectơ dòng điện không gian của stator được định nghĩa như sau:
2

𝚤̅𝑠 = (𝑖𝑎𝑠 + 𝑎𝑖𝑏𝑠 + 𝑎2 𝑖𝑐𝑠 )
3

(1.3)

Trong đó ı̅s là 1 vectơ có modul không thay đổi, quay trên mặt phẳng phức với tốc

độ góc ωe=2πfs và tạo với trục thực 1 góc γ = ωe(t). Dòng điện của từng pha chính là hình
chiếu của vectơ mới thu được lên trục của cuộn dây pha tương ứng. Đặt hệ tọa độ (αβ)
trong hệ không gian động cơ. Hai hình chiếu của vecto dòng điện được đặt tên là hai dòng
điện iαs và iβs. Vectơ dòng điện stato có thể viết ở dạng:
𝚤̅𝑠 = 𝑖𝛼𝑠 + 𝑗𝑖𝛽𝑠

(1.4)
5


Với giả thiết dòng điện 3 pha là đối xứng, các thành phần thứ tự không bằng
không, hai thành phần iαs và iβs là thành phần dòng điện hình sin. Được tính từ các thành
phần dòng điện ở các pha.

(1.5)

Phép biến đổi ngược lại có thể được thực hiện bởi phương trình sau:
𝑖𝑎𝑠
�𝑖𝑏𝑠 � =
𝑖𝑐𝑠

1
0
⎡ 1 √3 ⎤
𝑖𝛼𝑠
⎢− 2
2 ⎥ �
𝑖𝛽𝑠 �
⎢ 1

√3
⎣ 2 − 2⎦

(1.6)

β

Trục pha B

ic
iβs
Trục pha A (trục α)
ib



Hình 1.4. Biểu diễn vecto không gian trên hệ tọa độ dq

� 𝑟 quay
Hình 1.4 biểu diễn cả 2 vecto dòng điện stator 𝚤̅𝑠 và vecto từ thông rotor Ѱ
𝑑𝜃
với tốc độ góc 𝜔𝑠 = 2𝜋𝑓𝑠 = 𝑠�𝑑𝑡.
7


Hệ trục tọa độ dương trong hình 1.4 có trục thực trùng với trục của vecto từ thông
rotor hệ trục này quay xung quanh điểm gốc chung với tốc độ 𝜔𝑠 .
Hình chiếu của 𝚤̅𝑠 lên hệ trục mới này là Ids và Iqs.
Với quy ước như sau:

𝚤̅𝑠𝑠 : Vecto dòng điện stator quan sát trên hệ tọa độ cố định stator (hệ 𝛼𝛽)
𝑓

𝚤̅𝑠 : Vecto dòng điện stator quan sát trên hệ tọa độ từ thông roto (hệ dq)
Quy ước này được dùng cho các đại lượng khác một cách tương tự.
Từ đó có:
𝚤̅𝑠𝑠 = 𝑖𝛼𝑠 + 𝑗𝑖𝛽𝑠

(1.8)

𝑓

𝚤̅𝑠 = 𝐼𝑑𝑠 + 𝑗𝐼𝑞𝑠

Nhìn vào hình 1.4 thấy rằng, nếu biết góc 𝜃𝑠 , có thể tính được:

𝑟
𝑑𝑡

� 𝑘𝑠
+ 𝑗𝜔𝑘 Ѱ

� 𝑘𝑟
+ 𝑗(𝜔𝑘 − 𝜔𝑟 )Ѱ

(1.11)
(1.12)

Một cách tương tự, có thể biểu diễn cho các đại lượng vecto từ thông stato, dòng
điện roto.

8


1.2. Mô hình toán học động cơ khồng đồng bộ
1.2.1. Hệ phương trình cơ bản
Sử dụng khái niệm vecto không gian viết các phương trình cân bằng điện áp stato
và roto của động cơ KĐB:
𝑑Ѱ𝑎𝑠

𝑢𝑎𝑠 = 𝑅𝑠 . 𝑖𝑎𝑠 +
𝑢𝑏𝑠 = 𝑅𝑠 . 𝑖𝑏𝑠 +
𝑢𝑐𝑠 = 𝑅𝑠 . 𝑖𝑐𝑠 +

𝑑𝑡


𝑢�𝑟 = 𝑅𝑟 . 𝚤̅𝑟 +

�𝑠
𝑑Ѱ

(1.15)

�𝑟
𝑑Ѱ

(1.16)

𝑑𝑡

𝑑𝑡

Phương trình từ thông stato và roto:
� 𝑠 = 𝐿𝑠 . 𝚤̅𝑠 + 𝐿𝑚 . 𝚤̅𝑟
Ѱ

(1.17)

� 𝑟 = 𝐿𝑟 . 𝚤̅𝑟 + 𝐿𝑚 . 𝚤̅𝑠
Ѱ

(1.18)

Phương trình mômen quay:
𝑀=


.

𝑑𝜔𝑟

(1.20)

𝑑𝑡

J:

Mômen quán tính

𝜔𝑟 :

Tốc độ góc điện roto động cơ

pc:

Số đôi cực

𝑀, 𝑀𝑐 :

Mômen động cơ, mômen tải.

1.2.2. Các phương trình và mô hình trạng thái liên tục của động cơ trên hệ tọa độ 𝜶𝜷

Hệ phương trình vecto mô tả động cơ KĐB trên hệ tọa độ (𝛼, 𝛽) gắn trên stato của
động cơ nhận được từ (1.15 – 1.18) có dạng:
𝑢�𝑠𝑠 = 𝑅𝑠 𝚤̅𝑠𝑠 +
0 = 𝑅𝑟 𝚤̅𝑟𝑠 +



� 𝑠𝑠 = 𝐿𝑠 𝚤̅𝑠𝑠 + 𝐿𝑚 (Ѱ
� 𝑟𝑠 − 𝐿𝑚 𝚤̅𝑠𝑠 )
Ѱ
𝐿𝑟

Thay (1.22) vào hai phương trình trên của (1.21) được:

10

(1.22)


𝑢�𝑠𝑠 = (𝑅𝑠 + 𝑅𝑟
0=−

𝐿𝑚

𝑑𝑡

= −�

𝐿𝑟

𝑑𝚤̅𝑠𝑠
𝑑𝑡

𝑇𝑟



Ѱ𝑟𝑠

(1.23)

𝑑𝑡

𝐿𝑚

1

� 𝑟𝑠 +
� − 𝑗𝜔𝑟 � Ѱ

𝐿𝑟 𝐿𝑠 𝜎 𝑇𝑟

Đặt biến trạng thái của động cơ là [𝚤̅𝑠𝑠

� 𝑠𝑠 ]𝑇
Ѱ

1

𝐿𝑠 𝜎

𝑢�𝑠𝑠

(1.24)

Có hệ phương trình trạng thái ở dạng tổng quát mô tả động cơ không đồng bộ roto


𝑠
𝐴21
=

𝐿𝑚

𝐿𝑠 𝜎

𝐿𝑟

1

𝐼

� 𝐼 − 𝑗𝜔𝑟 �

𝐿𝑠 𝐿𝑟 𝜎 𝑇𝑟

𝐿𝑚
𝑇𝑟

𝑠
𝐴22
=−

𝐵1 =

2



Cũng có thể viết ( 1.25) ở dạng sau:

11

(1.25)


𝑑𝑖𝑎𝑠
𝑑𝑡

𝑑𝑖𝛽𝑠
𝑑𝑡

= −�

𝑑Ѱ𝛼𝑟
𝑑𝑡

𝑑Ѱ𝛽𝑟
𝑑𝑡

= −�
=

=

𝐿𝑚

1

1

𝑇𝑟

� 𝑖𝑎𝑠 +

� 𝑖𝛽𝑠 −

1−𝜎

𝜎𝑇𝑟 𝐿𝑚

1−𝜎

𝜎𝐿𝑚

Ѱ𝛼𝑟 +

𝜔𝑟 Ѱ𝛼𝑟 +

Ѱ𝛼𝑟 + 𝜔𝑟 Ѱ𝛽𝑟 +

𝑜𝑢�𝑠𝑠

được biểu thức mômen động cơ:
3
2

𝐿𝑚
𝐿𝑟


𝑢𝛼𝑠

𝑢𝛽𝑠

(1.26)

thay vào phương trình (1.18)

(1.27)

Từ (1.26) và (1.27) xây dựng được mô hình của động cơ không đồng bộ ba pha
roto lồng sóc trong hệ (𝛼𝛽) hình 1.5.

Hình 1.5. Mô hình động cơ không đồng bộ roto lồng sóc trên hệ tọa độ (𝛼𝛽)
Trên hình 1.5 có hệ số:

1

𝑇𝜎

=

1

𝜎𝑇𝑠



1−𝜎

𝑑Ѱ
𝑟
𝑑𝑡

+

𝑓

�𝑠
+ 𝑗𝜔𝑠 Ѱ

� 𝑟𝑓
+ 𝑗(𝜔𝑠 − 𝜔𝑟 )Ѱ

(1.28)

𝑓
𝐿𝑚 𝚤̅𝑟

� 𝑟𝑓 = 𝐿𝑚 𝚤̅𝑠𝑓 + 𝐿𝑟 𝚤̅𝑟𝑓
Ѱ

Tương tự trường hợp hệ tọa độ stato, ta tìm cách khử dòng roto cũng như từ thông
stato khỏi hệ phương trình và ta thu được hệ phương trình mô tả động cơ trên hệ tọa độ dq:
𝑓

𝑢�𝑠 = (𝑅𝑠 + 𝑅𝑟
0=−

𝐿𝑚

� 𝑟𝑓 + (𝜔𝑠 − 𝜔𝑟 )𝑗Ѱ
� 𝑟𝑓 + 𝑑Ѱ𝑟
Ѱ

𝑑𝑡

(1.29)

𝑑𝑡

Dạng tổng quát của phương trình trạng thái mô tả điều khiển trong hệ tọa độ dq có
dạng:
𝑓

𝑓

𝐴11

𝑓

𝚤̅𝑠

𝑑

𝐴12

𝑓
𝐵1 𝑓
𝚤̅𝑠


𝐿𝑠 𝜎

𝐿𝑟

1

𝐼 − 𝜔𝑠 𝐽

� 𝐼 − 𝜔𝑟 𝐽�

𝐿𝑠 𝐿𝑟 𝜎 𝑇𝑟

𝐿𝑚
𝑇𝑟

𝐴22 = −
𝐵1 =

2

𝐿
𝑅𝑠 +𝑅𝑟 𝑚
2

1

𝐿𝑠 𝜎

𝐼


𝑑𝑡

𝑑Ѱ𝑞𝑟
𝑑𝑡

1

𝜎𝑇𝑠

+

𝐼
𝑇𝑟 𝑑𝑠



1−𝜎
𝜎𝑇𝑟

=

=

𝐿𝑚

𝐿𝑚
𝑇𝑟

1


3

𝐿𝑚

𝑓

� 𝐼𝑑𝑠 + 𝜔𝑠 𝐼𝑞𝑠 +

1

= −𝜔𝑠 𝐼𝑑𝑠 − �

𝑑Ѱ𝑑𝑟
𝑑𝑡

= −�

1−𝜎

𝜎𝐿𝑚

𝜔𝑟 Ѱ𝑞𝑟 +

1−𝜎

𝜎𝑇𝑟 𝐿𝑚

Ѱ𝑞𝑟 +

1

Ѱ𝑑𝑟 𝐼𝑞𝑠

(1.32)

Hai phương trình (1.31) và (1.32) được kết hợp với nhau được mô hình trọn vẹn
của động cơ không đồng bộ roto lồng sóc trên hệ tọa độ dq (hình 1.7)

𝑢𝑠𝛼

𝑢𝑠𝛽

𝑢𝑠𝑑
𝑒 𝑗𝜃
𝑢𝑠𝑞

1
𝐿𝑠 𝜎

1
𝐿𝑠 𝜎



×
×







×

1
𝑠

Hình 1.6. Mô hình động cơ trên hệ tọa độ dq

14

𝑚𝑤

𝜔

𝜔𝑟

𝜔𝑠+

+


1.3. Tổng quan về hệ thống điều khiển vecto động cơ không đồng bộ
1.3.1. Nguyên lý điều khiển vectơ
Nguyên lý điều khiển vectơ dựa trên ý tưởng điều khiển động cơ không đồng
bộ tương tự như điều khiển động cơ một chiều. Sự tương tự của các phương pháp điều
khiển đó được trình bày trên hình 1.7.
Ở động cơ một chiều, bỏ qua phản ứng phần ứng và giả thiết mạch từ không bão
hoà, momen động cơ một chiều được tính như sau:
M = KΦmIư = K’IktIư
Trong đó:

𝐼𝑑𝑠

𝐼𝑞𝑠

a. Điều khiển một chiều

Mạch ĐK
nghịch lưu

ĐC

b. Điều khiển xoay chiều

Hình 1.7. Sự tương tự giữa phương pháp điều khiển động cơ một chiều
và điều khiển vecto động cơ không đồng bộ

15


Nhận thấy rằng dòng điện phần ứng và dòng điện kích từ (hoặc từ thông) không
phụ thuộc vào nhau. Vậy có thể điều khiển độc lập dòng kích từ và dòng phần ứng để đạt
được đặc tính mômen tối ưu. Duy trì dòng điện kích từ là không đổi (Ikt = const), mômen
được điều chỉnh bằng điều chỉnh dòng điện phần ứng.
Phương pháp điều khiển này có thể áp dụng cho động cơ không đồng bộ nếu sử
dụng khái niệm vectơ không gian. Với ý tưởng định nghĩa vectơ không gian dòng điện
của động cơ và mô tả động cơ ở hệ tọa độ quay với tốc độ đồng bộ ωs (hệ trục dq), các đại
lượng dòng điện, điện áp, từ thông là các đại lượng một chiều. Vectơ dòng điện có thể
tách ra làm hai thành phần trên hai trục vuông góc d và q: Ids và Iqs. Mômen động cơ
không đồng bộ là hàm của từ thông và dòng điện stato (1.32). Bằng cách chọn trục d
trùng với từ thông roto, phương trình từ thông, mômen sẽ viết lại như sau:

𝐼𝑑𝑠 , 𝐼𝑞𝑠 :

Hình chiếu của vecto từ thông roto trên trục d

Số đôi cực của động cơ

𝑝𝑐 :

Hằng số thời gian roto

p:

Toán tử laplace.

𝑇𝑟 :

Phương trình (1.34a) chỉ ra rằng từ thông roto có thể được tăng giảm gián tiếp
thông qua tăng giảm 𝐼𝑑𝑠 . Tuy nhiên, cần lưu ý rằng quan hệ (1.34a) giữa hai đại lượng

Ѱ𝑑𝑟 và 𝐼𝑑𝑠 là quan hệ trễ bậc nhất với hằng số thời gian 𝑇𝑟 . Nếu thành công trong việc áp

đặt nhanh chóng dòng 𝐼𝑑𝑠 , có thể coi 𝐼𝑑𝑠 là đại lượng điều khiển từ thông, nó tương ứng

với 𝐼𝑘 trong động cơ một chiều kích từ độc lập. Mặt khác, tại mỗi thời điểm làm việc của
16


động cơ, từ thông đã giữ ổn định, khi đó quan hệ (1.346b) cho thấy rằng mômen M có thể
được điều khiển thông qua việc điều khiển 𝐼𝑞𝑠 . Do đó 𝐼𝑞𝑠 còn được gọi là dòng tạo
mômen quay và nó tương ứng với dòng phần ứng 𝐼ư trong động cơ một chiều kích từ độc

Iq2

Id1
Is1

Iq

0

q

Is2
𝜔𝑠

d
a. Điều chỉnh thành phần mômen

Id2

𝜃𝑠1
𝜃𝑠2

q

Is1
Is2

d
b. Điều chỉnh thành phần từ thông


Nhờ tải bản gốc

Tài liệu, ebook tham khảo khác

Music ♫

Copyright: Tài liệu đại học © DMCA.com Protection Status