1. Trang chủ
  2. » Kỹ Thuật - Công Nghệ

TÍNH TOÁN CHU TRÌNH CÔNG TÁC TRONG ĐỘNG CƠ ĐỐT TRONG

35 5,3K 41
Tài liệu đã được kiểm tra trùng lặp

Đang tải... (xem toàn văn)

Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống

THÔNG TIN TÀI LIỆU

Thông tin cơ bản

Định dạng
Số trang 35
Dung lượng 3,26 MB

Nội dung

1- Công suất của động cơ Ne Ne =12 (mã lực) =8,83(Kw) 2- Số vòng quay của trục khuỷu n n =2200 (vg/ph) 3- Đường kính xi lanh D D =95 (mm) 4- Hành trình piton S S =115 (mm) 5- Dung tích công tác Vh

Trang 1

PHẦN I :TÍNH TOÁN CHU TRÌNH CÔNG TÁC TRONG ĐỘNG CƠ ĐỐT TRONG

I ) Trình tự tính toán :1.1 )Số liệu ban đầu :

1- Công suất của động cơ Ne Ne =12 (mã lực) =8,83(Kw) 2- Số vòng quay của trục khuỷu n n =2200 (vg/ph)

3- Đường kính xi lanh D D =95 (mm) 4- Hành trình piton S S =115 (mm) 5- Dung tích công tác Vh

Vh = = 0,81515 (dm3)

6- Số xi lanh i i = 1 7- Tỷ số nén ε ε =16 8- Thứ tự làm việc của xi lanh (1-2-4-3)

9- Suất tiêu hao nhiên liệu ge ge =180 (g/ml.h)

10- Góc mở sớm và đóng muộn của xupáp nạp α1 ; α2 α1 =10 (độ) α2 =29 (độ) 11- Góc mở sớm và đóng muộn của xupáp thải  1, 2 1 =32 (độ) 2 =7 (độ) 12- Chiều dài thanh truyền ltt ltt =205 (mm)

13- Khối lượng nhóm pitton mpt mpt =1,15 (kg) 14- Khối lượng nhóm thanh truyền mtt mtt =2,262 (kg)

1.2 )Các thông số cần chọn :

1 )Áp suất môi trường :pk

Áp suất môi trường pk là áp suất khí quyển trước khi nạp vào đông cơ (với đông cơ không tăng áp ta có áp suất khí quyển bằng áp suất trước khi nạp nên ta chọn pk =po Ở nước ta nên chọn pk =po = 0,1 (MPa)

2 )Nhiệt độ môi trường :Tk

Nhiệt độ môi trường được chọn lựa theo nhiệt độ bình quân của cả năm

Vì đây là động cơ không tăng áp nên ta có nhiệt độ môi trường bằng nhiệt độ trước xupáp nạp nên :

Tk =T0 =24ºC =297ºK

3 )Áp suất cuối quá trình nạp :pa

Áp suất Pa phụ thuộc vào rất nhiều thông số như chủng loại đông cơ ,tính năng tốcđộ n ,hệ số cản trên đường nạp ,tiết diện lưu thông… Vì vậy cần xem xét đông cơ đang tính thuộc nhóm nào để lựa chọn Pa

Áp suất cuối quá trình nạp pa có thể chọn trong phạm vi:

Trang 2

5 )Mức độ sấy nóng của môi chất ∆T

Mức độ sấy nóng của môi chất ∆T chủ yếu phụ thuộc vào quá trình hình thành hhkhí ở bên ngoài hay bên trong xy lanh

Với động cơ ddieeezeel : ∆T=20 ºC -40ºC Vì đây là đ/c D12-3 nên chọn ∆T=29,5ºC

Đối với động cơ đang tính là động cơ diesel có α > 1,4 có thể chọn λ=1,10

8 )Hệ số quét buồng cháy λ :

Vì đây là động cơ không tăng áp nên ta chọn λ =1

9 )Hệ số nạp thêm λ

Hệ số nạp thêm λ phụ thuộc chủ yếu vào pha phối khí Thông thường ta có thể chọn λ =1,02÷1,07 ; ta chọn λ =1,0316

10 )Hệ số lợi dụng nhiệt tại điểm z ξ :

Hệ số lợi dụng nhiệt tại điểm z ξ phụ thuộc vào chu trình công tác của động cơ Với các loại đ/c điezen ta thường chọn : ξ=0,70-0,85

Chọn : ξ=0,75

11 )Hệ số lợi dụng nhiệt tại điểm b ξ :

Hệ số lợi dụng nhiệt tại điểm b ξ tùy thuộc vào loại động cơ xăng hay là động cơ điezel ξ bao giờ cũng lớn hơn ξ

Với các loại đ/c điezen ta thường chọn : ξ =0,80-0,90 ta chọn ξ=0,85

12 )Hệ số hiệu chỉnh đồ thị công φ :

Thể hiện sự sai lệch khi tính toán lý thuyết chu trình công tác của động cơ với chutrình công tác thực tế Sự sai lệch giửa chu trình thực tế với chu trình tính toán của động cơ xăng ít hơn của động cơ điezel vì vậy hệ số φ của đ/c xăng thường chọn hệ số lớn.

Có thể chọn φ trong phạm vi: φ =0,92-0,97

Trang 3

Nhưng đây là đ/c điezel nên ta chọn φ =0,97

II )Tính toán các quá trình công tác :2.1 Tính toán quá trình nạp :

= 0,03823

2 )Nhiệt độ cuối quá trình nạp T

Nhiệt độ cuối quá trình nạp T đươc tính theo công thức:

T= ºK

T=

3 )Hệ số nạp η :

η = η= 1 116

 29729729,5 00,088,11

180.0,59059.297 = 0,8191 (kmol/kg nhiên liệu)

5 )Lượng không khí lý thuyết cần để đốt cháy 1kg nhiên liệu M :

Lượng kk lý thuyết cần để đốt cháy 1kg nhiên liệu M được tính theo công thức : M = (kmol/kg) nhiên liệu

Trang 4

Vì đây là đ/c điezel nên ta chọn C=0,87 ; H=0,126 ;O=0,004 M = ( + - ) =0,4946 (kmol/kg) nhiên liệu

2 )Tỉ nhiệt mol đẳng tích trung bình của sản phạm cháy :

Khi hệ số lưu lượng không khí α >1 tính theo công thức sau :

= + 10 T (kJ/kmol.độ) = 

427 10=20,8537 (kJ/kmol.độ)

3 )Tỉ nhiệt mol đẳng tích trung bình của hỗn hợp :

Tỉ nhiệt mol đẳng tích trung bình của hh trong quá trình nén tính theo công thức sau :

= = 19,845382

(kJ/kmol.độ)

4 ) Chỉ số nén đa biến trung bình n:

Chỉ số nén đa biến trung bình phụ thuộc vào thong số kết cấu và thong số vận hành như kích thước xy lanh ,loại buồng cháy,số vòng quay ,phụ tải,trạng thái nhiệt độ của động cơ…Tuy nhiên n tăng hay giảm theo quy luật sau :

Tất cả những nhân tố làm cho môi chất mất nhiệt sẽ khiến cho n tăng.Chỉ số nén đa biến trung bình n được xác bằng cách giải phương trình sau :

n-1 =

Chú ý :thông thường để xác định được n ta chọn n trong khoảng 1,340÷1,390

Rất hiếm trường hợp đạt n trong khoảng 1,400 ÷ 1,410

→ (theo sách Nguyên Lý Động Cơ Đốt Trong - trang 128 )

Vì vậy ta chọn n theo điều kiện bài toán cho đến khi nao thõa mãn điều kiện bài toán :thay n vào VT và VP của phương trình trên và so sánh,nếu sai số giữa 2 vế của phương trình thõa mãn <0,2% thì đạt yêu cầu Sau khi chọn các giá trị của n ta thấy n =1,3685 thõa mãn điều kiện bài toán

5 )Áp suất cuối quá trình nén P :

Áp suất cuối quá trình nén P được xác định theo công thức : P = P ε = 0,088 161,3685 = 3,9037 (MPa)

6 )Nhiệt độ cuối quá trình nén T

Nhiệt độ cuối quá trình nén T được xác định theo công thức

Trang 5

T = T ε = 340,8 161,36851 = 944,9 ( ºK )

7 )Lượng môi chất công tác của quá trình nén M :

Lượng môi chất công tác của quá trình nén M được xác định theo công thức : M = M+ M = M.(1r) = 0,8191.(1+0,03823) = 0,85

2.3 )Tính toán quá trình cháy :

1 )Hệ số thay đổi phân tử lí thuyết β :

Ta có hệ số thay đổi phần tử lý thuyết β được xác định theo công thức : β = = = 1+

Trong đó độ tăng mol ΔM của các loại động cơ được xác định theo công thức sau:M của các loại động cơ được xác định theo công thức sau: ΔM của các loại động cơ được xác định theo công thức sau:M = 0,21.(1-α)M + ( +  )

Đối với động cơ điezel : ΔM của các loại động cơ được xác định theo công thức sau:M = ( + ) Do đó

β = 1 + = 1 +

= 1,0386

2 )Hệ số thay đổi phân tư thưc tế β: ( Do có khí sót )

Ta có hệ số thay đổi phân tử thực tế β được xác đinh theo công thức : β = = 1,0386 0,0382 1, 0372

1 0,0382

3 )Hệ số thay đổi phân tử thực tế tại điểm z β : (Do cháy chưa hết )

Ta có hệ số thay đổi phân tư thực tế tại điểm z β được xác định theo công thức : β = 1 + χ

Trong đó

χ = = 0,75 0,8824

Nên: β =1 + 1,0386 10,8824 1,03281 0,0382

5 )Nhiệt độ tại điểm z T :

* Đối với động cơ điezel,tính nhiệt độ T bằng cách giải pt cháy :

+ T = β T Trong đó :

Q : là nhiệt trị của dầu điezel ,Q =42,5 10 ( kJ/kgn.l ) :là tỉ nhiệt mol đẳng áp trung bình của sản vật cháy tại z là : =8,314+

:là tỉ nhiệt mol đẳng tích trung bình của sản vật cháy tại z được tính theo ct :

Trang 6

= = a'' + b'' T Chỉnh lý lại ta có : = a'' + b'' T Thay (2) vào (1) ta được:

0,75.42500 (19,845 8,314).944,9 1,0382.

Giải phương trình trên với a'' =29,05697 ; b''=0,00264 ta được : T =2032,7 ; T= -6524,06 (loại)

6 )Áp suất tại điểm z p :

Ta có áp suất tại điểm z p được xác định theo công thức : p =λ P ( MPa )

Với λ là hệ số tăng áp λ= β

T :là nhiêt trị tại điểm b và được xác định theo công thức :

T= ( ºK )

Q :là nhiệt trị tính toán

Đối với động cơ điezel Q= Q Q = 42.500 (kJ/kg n.l)

Qua kiệm nghiêm tính toán thì ta chọn đươc n =1,244.Thay n vào 2 vế của pt trên ta so sánh ,ta thấy sai số giữa 2 vế <0,2% nên n chọn là đúng

Trang 7

4 )Nhiệt độ cuối quá trình giãn nở T :

T= = 1,244 1

1088, 4

5 )Áp suất cuối quá trình giãn nở p :

Áp suất cuối quá trình giãn nở P được xác định theo CT : p = = 1244

0,111088, 4.

3 )Suất tiêu hao nhiên liệu chỉ thị g :

Ta có công thức xác định suất tiêu hoa nhiên liệu chỉ thị g: g= =

432.10 0,8139.0,1

196, 27

4 )Hiệu suất chỉ thi η:

Ta có công thức xác định hiệu suất chỉ thị η : η = =

0, 4316

5 )Áp suất tổn thất cơ giới P :

Áp suất tổn thất cơ giới được xác định theo nhiều công thức khác nhau và đươc biểu diễn bằng quan hệ tuyến tính với tốc độ trung bình của động cơ.Ta có tốc độ trung bình của động cơ là :

V = = 115.10 22003 8, 43330

 (m/s) Đối với động cơ diesel cao tốc dung cho ôtô (V >7) :

P= 0,015+0,0156.V= 0,015+0,0156.8,433 = 0,1466 (MPa)

6 )Áp suất có ích trung bình P :

Ta có công thức xđ áp suất có ích trung bình thực tế được xđ theo CT : P = P – P =0,7364 – 0,1466 = 0,5898 (MPa)

Trang 8

Ta có trị số P tính quá trình nạp P (nạp) =0,6768 va P=0,6736 thì không có sự chênh lệch nhiều nên có thể chấp nhận được

7 )Hiệu suất cơ giới η :

Ta có có thức xác định hiệu suất cơ giới: η = e

8 )Suất tiêu hao nhiên liệu g :

Ta có có thức xác định suất tiêu hao nhiên liệu tính toán là:

D = 4.0,815153

3,14.115.10 = 0,9502419(mm) Ta có sai số so với đề bài là :0,045 (mm)

III ) Vẽ và hiệu đính đồ thị công :

Căn cứ vào các số liệu đã tính pr, p , p , p , p ,n, n, ε ta lập bảng tính đường nén và

đường giãn nở theo biến thiên của dung tích công tác V = i.V V : Dung tích buồng cháy

V = = 0,81515

16 1 = 0,054343 ( dm3 )

Các thông số ban đầu: p = 0 ,11 MPa ; p = 0,088MPa; p= 3,9037 MPa

p = 7,027 MPa ; p = 0,2903 MPa

3.1 ) Xây dựng đường cong áp suất trên đường nén :

- Phương trình đường nén đa biến :

3.2 ) Xây dựng đường cong áp suất trên quá trình giãn nở :

- Phương trình của đường giãn nở đa biến :

Trang 9

V = ρ.V Vậy P = P = = = P

n : Chỉ số giãn nở đa biến trung bình n = 1,2438

P : Áp suất tại điểm z : P = 7,027 (MPa)

=1,1654 0,10733 4,0080 (14,2; 102,8) 9,8832 (14,2; 250) 2 0,18418 1,9148 (24,5; 49) 5,0487 (24,5; 129,4) 3 0,27627 1,0997 (36,7; 28) 3,0490 (36,7; 78) 4 0,36836 0,7419 (49; 19,02) 2,1318 (49; 54,6) 5 0,46045 0,5468 (61,1; 14) 1,6152 (61,1; 41,1) 6 0,55254 0,4261 (73,3; 10,9) 1,2875 (73,3; 33) 7 0,64463 0,3451 (85,6; 8,84) 1,0628 ( 85,6; 27,3) 8 0,73672 0,2875 (97,8; 7,37) 0,9002 (97,8; 23) 9 0,82881 0,2447 (110; 6,27) 0,7775 (110; 20) 10 0,92090 0,2119 (122,3; 5,43) 0,6820 (122,3; 17,4) 11 1,01299 0,1860 (134,5; 4,76) 0,6058 (134,5; 15,5) 12 1,10508 0,1651 (146,8; 4,23) 0,5436 (146,8; 14) 13 1,19718 0,1480 ( 159; 3,79) 0,4921 (159;12,6) 14 1,28927 0,1337 (171,2; 3,42) 0,4488 (171,2; 11,5) 15 1,38136 0,1217 (183,4; 3,12) 0,4119 (183,4; 10,6) 16 1,47345 0,1114 (195,6; 2,86) 0,3801 (195,6; 9,74) 17 1,56554 0,1025 (207,9; 2,62) 0,3525 (207,9; 9,03) 18 1,65763 0,0948 (220; 2,43) 0,3288 (220; 8,43)

Trang 10

3.4 ) Vẽ vòng tròn Brick đặt phía trên đồ thị công :

Ta chọn tỉ lệ xích của hành trình piton S là : μ = = = 115

225 15 = 0,546 Thông số kết cấu động cơ là :

λ = = = 115

2.205 = 0,28 Khoảng cách OO’ là :

OO’= = 0, 28.57,5

2 = 8,05 ( mm ) Giá trị biểu diễn của OO’ trên đồ thị :

gtbd = = 0,5468,05 = 14,91( mm ) Ta có nửa hành trình của piton là :

R = =115

2 =57,5 ( mm ) Giá trị biểu diễn của R trên đồ thị :

gtbd = = 0,54657,5 = 105,31 ( mm )

3.5 ) Lần lượt hiệu định các điểm trên đồ thị :

1 ) Hiệu đính điểm bắt đầu quá trình nạp : (điểm a)

Trang 11

Từ điểm O’ trên đồ thị Brick ta xác định góc đóng muộn xupáp thải β , bán kính này cắt đường tròn tại điểm a’ Từ a’ gióng đường thẳng song song với trục tung cắt đường P tại điểm a Nối điểm r trên đường thải ( là giao điểm giữa đường P và trục tung ) với a ta được đường chuyển tiếp từ quá trình thải sang quá trình nạp

2 ) Hiệu định áp suất cuối quá trình nén : ( điểm c’)

Áp suất cuối quá trình nén thực tế do hiện tượng phun sớm (động cơ điezel ) và hiện tượng đánh lửa sớm (động cơ xăng ) nên thường chọn áp suất cuối quá trình nén lý thuyết P đã tính Theo kinh nghiệm , áp suất cuối quá trình nén thực tế P’ được xác định theo công thức sau :

Vì đây là động cơ điezel :

P’ = P + ( P - P ) = 3,9037 + ( 7,027- 3,9037 ) = 4,9448 ( MPa ) Từ đó xác định được tung độ điểm c’trên đồ thị công :

y = = 4,94480, 2778 = 178,0128 (mm )

3 ) Hiệu chỉnh điểm phun sớm : ( điểm c’’ )

Do hiện tương phun sớm nên đường nén trong thực tế tách khỏi đường nén lý thuyết tại điểm c’’ Điểm c’’ được xác định bằng cách Từ điểm O’ trên đồ thị Brick ta xác định được góc phun sớm hoặc góc đánh lửa sớm θ, bán kính này cắt vòng tròn Brick tại 1 điểm Từ điểm gióng này ta gắn song song với trục tung cắt đường nén tại điểm c’’ Dùng một cung thích hợp nối điểm c’’ với điểm c’

4 )Hiệu đính điểm đạt P thực tế

Áp suất p thực tế trong quá trình cháy - giãn nở không duy trì hằng số như động cơ điezel ( đoạn ứng với ρ.V ) nhưng cũng không đạt được trị số lý thuyết như động cơ xăng Theo thực nghiệm ,điểm đạt trị số áp suất cao nhất là điểm thuộc miền vào khoảng 372° ÷ 375° ( tức là 12° ÷ 15° sau điểm chết trên của quá trình cháy và giãn nở )

Hiệu định điểm z của động cơ điezel :

- Xác định điểm z từ góc 15º Từ điểm O΄trên đồ thị Brick ta xác định góc tương ΄trên đồ thị Brick ta xác định góc tương trên đồ thị Brick ta xác định góc tương ứng với 375º góc quay truc khuỷu ,bán kính này cắt vòng tròn tại 1 điểm Từ điểm này ta gióng song song với trục tung cắt đường P tại điểm z

- Dùng cung thích hợp nối c’ với z và lượn sát với đường giãn nở

5 ) Hiệu định điểm bắt đầu quá trình thải thực tế : ( điểm b’ )

Do có hiện tượng mở sớm xupáp thải nên trong thực tế quá trình thải thực sự diễn ra sớm hơn lý thuyết Ta xác định điểm b bằng cách : Từ điểm O’trên đồ thị Brick ta xác định góc mở sớm xupáp thải β,bán kính này cắt đường tron Brick tại 1 điểm.Từ điểm này ta gióng đường song song với trục tung cắt đường giãn nở tại điểm b’.

6 ) Hiệu định điểm kết thúc quá trình giãn nở : ( điểm b’’ )

Áp suất cuối quá trình giãn nở thực tế P thường thấp hơn áp suất cuối quá trình giãn nở lý thuyết do xupáp thải mở sớm Theo công thức kinh nghiệm ta có thể xác định được :

P= P + ( P - P ) = 0,11 + ( 0,2903- 0,11 ) = 0,2003 (MPa) Từ đó xác định tung độ của điểm b’’ là :

Trang 12

y = = 0,027780, 2003 = 7,209 ( mm )

Trang 13

Đồ thị công chỉ thị

PHẦN II : TÍNH TOÁN ĐỘNG HỌC VÀ ĐỘNG LỰC HỌCI ) Vẽ đường biểu diễn các quy luật động học :

Các đường biểu diễn này đều vẽ trên 1 hoành độ thống nhất ứng với hành trình piston S = 2R Vì vậy độ thị đều lấy hoành độ tương ứng với V của độ thị công ( từ điểm 1.V đến ε.V )

Trang 14

1.1 ) Đường biểu diễn hành trình của piston x = ƒ(α)α))

Ta tiến hành vẽ đường biểu diễn hành trình của piston theo trình tự sau : 1 Chọn tỉ xích góc : thường dùng tỉ lệ xích ( 0,6 ÷ 0,7 ) ( mm/độ ) 2 Chọn gốc tọa độ cách gốc cách độ thị công khoảng 15 ÷ 18 cm

3 Từ tâm O’ của đồ thị Brick kẻ các bán kính ứng với 10° ,20° ,…….180° 4 Gióng các điểm đã chia trên cung Brick xuống các điểm 10° ,20° ,…….180° tương ứng trên trục tung của đồ thị của x = ƒ(α) ta được các điểm xác định chuyển vị x tương ứng với các góc 10°,20°,… 180°

5 nối các điểm xác định chuyển vị x ta được đồ thị biểu diễn quan hệ x = f(α).

1.2 ) Đường biểu diễn tốc độ của piston v = f(α)α))

Ta tiến hành vẽ đường biểu diễn tốc độ của píton v = f(α) Theo phương pháp đồ thị vòng Tiến hành theo các bước cụ thể sau:

1.Vẻ nửa vòng tròn tâm O bán kính R ,phía dưới đồ thị x = f(α) Sát mép dưới của bản vẽ

5 Nối tại các điểm a,b,c,… Tạo thành đường cong giới hạn trị số của tốc độ piton thể hiện bằng các đoạn thẳng song song với tung độ từ các điểm cắt vòng tròn bán kínhR tạo với trục hoành góc α đến đường cong a,b,c….

Đồ thị này biểu diễn quan hệ v = f(α) trên tọa độ độc cực :

Hinh 2.1: Dạng đồ thị v = f(α) 1.3 Đường biểu diễn gia tốc của piston j = f(α) x)

Ta tiến hành vẽ đường biểu diễn gia tốc của piston theo phương pháp Tôlê ta vẽ theo các bước sau :

1.Chọn tỉ lệ xích μ phù hợp trong khoảng 30 ÷ 80 (m/s mm ) Ở đây ta chọn μ = 50 (m/s mm )

V= f  

Trang 15

4 57

O O'

j = –R.ω.( 1– λ ) = –57,5.10.230,3835.( 1–0,28 ) = –2,197.10( m/ s) Vậy ta được giá trị biểu diễn của j là :

gtbd = = 2,197.103

 = –43,497 ( mm ) -Xác định vị trí của EF :

EF = –3.R.λ.ω = –3.57,5.10.0,28.230,3835 = –2,563.10 ( m/s ) Vậy giá trị biểu diễn EF là :

x= f( α)

Trang 16

j= f( x)

v= f( α )

II )Tính toán động học :

2.1 )Các khối lượng chuyển động tịnh tiến :

- Khối lượng nhóm piton m = 3,5 Kg

- Khối lượng thanh truyền phân bố về tâm chốt piston

+ ) Khối lương thanh truyền phân bố về tâm chốt piston m có thể tra trong các các sổ tay ,có thể cân các chi tiết của nhóm để lấy số liệu hoặc có thể tính gần đúng theo bản vẽ

+ ) Hoặc có thể tính theo công thức kinh nghiêm sau : Đối với động cơ điezel ta có :

m = (0, 28 0, 29) mtt

Trong đó mttlà khối lượng thanh truyền mà đề bài đã cho Ta chọn m = 0,28 m = 0,28 2,262= 0,63336

Trang 17

Vậy ta xác định đươc khối lượng tịnh tiến mà đề bài cho là : m = m + m = 1,15 + 0,63336 = 1,78336 (Kg)

2.2 ) Các khối lượng chuyển động quay :

m = π .ρ Trong đó ta có :

d : Là đường kính ngoài của chốt khuỷu : 65 δ : Là đường kính trong của chốt khuỷu : 26 l : Là chiều của chốt khuỷu : 47

ρ : Là khối lượng riêng của vật liệu làm chốt khuỷu ρ : 7800 Kg/ m = 7,8.10 Kg/ mm

Lực quán tính chuyển động tịnh tiến :

P = - m.j = -m.R.ω.( cos α + λ.cos 2α ) = -22,9.103( cos α + λ.cos 2α ) Với thông số kết cấu λ ta co bảng tính P :

Ngày đăng: 26/04/2013, 12:07

TỪ KHÓA LIÊN QUAN

TÀI LIỆU CÙNG NGƯỜI DÙNG

TÀI LIỆU LIÊN QUAN

w