THÔNG TIN TÀI LIỆU
i H c Qu c Gia Tp H Chí Minh TR NG I H C BÁCH KHOA PHAN CƠNG BÌNH NGHIÊN C U S THÔNG S NH H NG C A CÁC K THU T TRONG THI T B C P RUNG N N NG SU T Chuyên ngành: Công ngh ch t o máy LU N V N TH C S TP H CHÍ MINH, tháng 07 n m 2009 CƠNG TRÌNH C HỒN THÀNH T I TR NG I H C BÁCH KHOA I H C QU C GIA TP H CHÍ MINH Cán b h ng d n khoa h c: Cán b ch m nh!n xét 1: Cán b ch m nh!n xét 2: Lu!n v n th c s" # $c b%o v t i H&I LU N V N TH C S TR NG NG CH M B O V' I H C BÁCH KHOA, Ngày…….tháng… n m 2009 TR NG I H C BÁCH KHOA PHÒNG ÀO T O S H C&NG HÒA XÃ H&I CH NGH A VI'T NAM &C L P-T DO-H NH PHÚC ……………………… …………………………… Tp HCM, ngày….tháng….n m 2009 NHI M V LU N V N TH C S H tên h c viên: Ngày, tháng, n m sinh: Chun ngành: Phan Cơng Bình 14/06/1982 Công ngh ch t o máy Phái: Nam N i sinh: Bình D ng MSHV: 00407203 I – TÊN ( TÀI: NGHIÊN C U S NH H NG C A CÁC THÔNG S THU T TRONG THI T B C P RUNG K N N NG SU T II – NHI'M V) VÀ N&I DUNG: Kh o sát nh h ng thông s k thu t thi t b c p rung - Thi t l p toán ng l c h c - Mô ph ng k t qu nh h - Thi t k ch t o mơ hình th nghi m - n n ng su t: ng c a thông s thi t b b ng ph n m m Matlap a nh ng thơng s có nh h ng l n tìm mi n n ng su t cao nh t III – NGÀY GIAO NHI'M V): 05/02/2009 IV – NGÀY HOÀN THÀNH NHI'M V): 03/07/2009 V – CÁN B& H *NG D+N: T.S L u Thanh Tùng CÁN B& H *NG D+N CN B& MÔN QL CHUYÊN NGÀNH T.S L u Thanh Tùng N i dung c ng lu n v n th c s! "c H i #ng chuyên ngành thông qua Ngày ….,tháng … n m 2009 TR NG PHÒNG T-S H TR NG KHOA QL NGÀNH iv L I CAM OAN Tôi xin cam oan ây cơng trình nghiên c u c a riêng Các s li u, k t qu nêu lu n v n trung th c ch a t ng c công b b t k cơng trình khác Ký tên Phan Cơng Bình xiii DANH M C CÁC HÌNH V , S hi u hình v TH Tên hình v Trang 1.1 Dây chuy n cung c p chi ti t 1.2 Thi t b c p rung 1.3 C c u phân lo i 1.4 Mu ng b t gi t 2.1 S chi ti t chuy n 2.2 S phân tích l c hành trình thu n 2.3 2.4 nh h ng chi ti t 16 ng máng th gia t c t i h n J1 S 19 19 21 phân tích l c hành trình ng c 22 2.5 th gia t c t i h n J2 24 2.6 th gia t c t i h n J0 25 2.7 th tr ng thái chuy n ng chi ti t 27 2.8 Mơ hình ba b c t 39 2.9 S phân tích l c c a ph u 40 2.10 S phân tích l c c a 40 2.11 Bi u di n hình h c quy ho ch d ng B không gian 46 3.1 Chuy n 56 3.2 ng chi ti t b ném lên th chuy n ng áp ng chi ti t 3.3 th li ω = 100π giao ng c a máng F=18668N 3.4 th li dao ng giai o n ch a bình n 62 63 64 xiv 3.5 th li 3.6 th biên dao ng theo t n s 65 3.7 th biên dao ng A=0.298 t i t n s f=50hz 66 3.8 th biên dao ng theo t n s giai o n bình n 68 th biên dao ng theo t n s giai o n ch a bình 3.10 th biên A t n s f giai o n bình n 70 3.11 th biên A t n s f giai o n ch a bình n 71 3.12 th biên A t n s f giai o n bình n 72 3.13 th biên A t n s f giai o n ch a bình n 72 3.14 th mi n gi i h n thông s 3.9 dao ng giai o n bình n n 64 69 73 4.1 Ph u 76 4.2 Lò xo 76 4.3 T m 77 4.4 Mơ hình thi t k VBF 78 4.5 Ki m tra mơ hình 78 4.6 Mơ hình th c t VBF 79 5.1 th n ng su t thi t b tính tốn mơ ph ng 89 5.2 th n ng su t thi t b th nghi m 91 5.3 T ng quan mô ph ng th n ng su t th nghi m tính tốn 92 ABSTRACT Vibratory bowl feeders (VBF) are the most versatile devices used in the feeding of small mechanical parts during part assembly processes They are used in high volume, high-speed production This work develops a mathematical model of a bowl feeder to evaluate productivity base on vibration amplitude and vibration frequency From the mathematical method, a computer simulation can predict the velocity of the parts in the bowl or productivity of VBF In addition, experiment’s results were compared with computer simulation to determine the technical parameters, which have a large influence on productivity of VBF The results of this research might be used to optimize the design of VBF or provide data for future work TÓM T T VBF thi t b a n ng nh t c s d ng vi c cung c p nh ng chi ti t c khí nh su t trình l p ráp Chúng c s d ng s n xu t hàng kh i, s n xu t n ng su t cao Nghiên c u phát tri n mơ hình tốn h c c a ph u o n ng su t thi t b d a biên ch t n s dao ng trình mơ ph ng máy tính có th d ng T mơ hình tốn h c, ốn v n t c v n chuy n chi ti t ph u hay n ng su t thi t b Ngoài ra, k t qu th c nghi m ch ng tình mơ ph ng máy tính xác c so sánh v i v i nh thông s k thu t có s n ng su t VBF K t qu nghiên c u có th c s d ng thi t k VBF ho c làm tài li u cho nh ng nghiên c u ti p theo nh h ng n t i u thơng s vii Trang TRANG PH BÌA .i L I CAM OAN .iv I C M N v TÓM T T LU N V N .vi M C L C vii DANH M C CÁC KÝ HI U, CÁC CH VI T T T ix DANH M C CÁC B NG xii DANH M C CÁC HÌNH V , TH xiii xv VIBRATORY BOWL FEEDER i thi u VBF 1.2 M c tiêu ph m vi nghiên c u 16 .17 ! "#$ %&'(# %$ )*% +',- &/ $0#$ 1$ % )*% %2)/# '(#3 17 2.2 Thi t l p mơ hình ng h c c a ph u 39 2.3 Lý thuy t quy ho ch th c nghi m (QHTN) b c d ng B 44 56 54 viii :'(1 ; ?&, 1'(1 %$&/#3 =&* @AB %$C',% 54 ! D )E# %$'A ;&F 1'(1 %$&/#3 =&* A, ω , α 60 99 "#$ %&'(# G'H &/ -$&I#3 J'& ;&,#3 1CI' 1> $), ?'K#3 D'%+'- 63 9L ng su t thi t b VBF 73 THI T K -CH T O MƠ HÌNH TH! NGHI M 75 4.1 Thi t k mơ hình 75 4.2 Ki"m tra mơ hình thi t k 78 4.3 Ch t o mơ hình 79 4.4 Ti n hành thí nghi m 80 K T QU VÀ TH O LU N 87 5.1 Tính tốn n ng su t theo mô ph#ng 87 5.2 Tính tốn n ng su t theo th$ nghi m 89 5.3 Nh%ng v n & th'o lu n 92 K T LU N VÀ KI N NGH D7 N O ML O N$CP +CP1 : 96 $ )*% +',- 1$()ng trình tính tốn k t qu' b*ng Maple N$CP CP1 Q 106 Ch()ng trình tính tốn mơ ph#ng k t qu' b*ng Matlap ix B NG DANH M C CÁC KÝ HI U Ký hi u vi t t t β Ý ngh a Góc nghiêng máng ϖ T n s góc µ H s ma sát chi ti t ρ Kh i l f cb T n s dao ng c fr T n s dao ng riêng σu ng riêng c a ng b c ng su t b n u n µ H s poison λ H s λ1 Chuy n v xoay c a t m λ2 Chuy n v xoay c a ph u A Biên b ng l c dao ng dày c a máng B Chi u r ng c a máng xo n bd H s gi m ch n chuy n v d bh H s gi m ch n xoay cao su bv H s gi m ch n cao su D d ng kính c a ph u Chuy n v ph u > L c kích rung c n thi t cho lò xo > > L c kích rung c n thi t cho thi t b > Gia t c J gia t c t i h n J1 > Gia t c J gia t c ti t ki m n ng l ng Jt > Gia t c J gia t c t i h n J2 > > Gia t c J gia t c t i h n Jo > GVHD:TS L U THANH TÙNG HVTH: PHAN CÔNG BÌNH > V n t c Vo v n t c gi i h n > > Biên A bi n d ng l n nh t c a lị xo > B ng 2A Ch ng trình v th gia t c t i h n b ng Maple > > > > > > > > > > > > GVHD:TS L U THANH TÙNG HVTH: PHAN CƠNG BÌNH > > > > > > > GVHD:TS L U THANH TÙNG HVTH: PHAN CƠNG BÌNH B ng 3A Ch ng trình tính tốn ma tr n b ng Maple: > > > > > > Nhân ma tr n GVHD:TS L U THANH TÙNG HVTH: PHAN CƠNG BÌNH GVHD:TS L U THANH TÙNG HVTH: PHAN CƠNG BÌNH Ghép ma tr n [A11], [A12], [A21], [A22] có ma tr n [A] GVHD:TS L U THANH TÙNG HVTH: PHAN CƠNG BÌNH Ghép ma tr n [B1], [B2] có ma tr n [B]: GVHD:TS L U THANH TÙNG HVTH: PHAN CƠNG BÌNH Sau tìm c ma tr n [A] [B] dùng ph n m m MatLap gi i C n tìm nghi m d, • X = [ A] X + [ B ]Fd y1 λ1 • d Có d = [C ] Y = [C ] y• [C ] = [0 0 0] • λ1 • d GVHD:TS L U THANH TÙNG HVTH: PHAN CƠNG BÌNH PH L C B B ng 1B Ch ng trình mơ ph ng b ng Maplap %Dao ng h b c t c a thi t b VBF %Dao ng c a ph u %PT dao ng d lò xo clear; m1 = 90; m2 = 30; J1 = 2.38; J2 = 1.35; r0 = 2; r1 = 265; r2 = 265; theta = (70*pi)/180; Ks = 2.687*10^6; Kh = 393700; Kv = 1592520; bh = 3563; bv = 4783.5; bd = 0; F = 18668; f=50; w=2*pi*f; ts=0.0001; t=0:ts:10; A=[0 0 0; 0 0 0; 0 0 1; GVHD:TS L U THANH TÙNG HVTH: PHAN CƠNG BÌNH -J2*J1*Kv/(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2)m1*(J1*r2+J2*r1)*Kv/(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2) m1*J1*r2*Kh*r0^2/((m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2)*r1)m1*(J1*r2+J2*r1)*Kh*r0^2/(cos(theta)*(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1 *m2)*r1) m1*(J1*r2+J2*r1)*Ks/(cos(theta)*(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2)) -J2*J1*bv/(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2)m1*(J1*r2+J2*r1)*bv/(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2) m1*J1*r2*bh*r0^2/((m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2)*r1)m1*(J1*r2+J2*r1)*bh*r0^2*sin(theta)/(cos(theta)*(m1^2*J1*r2+m1^2*J2*r1+J2*J1* m1+J2*J1*m2)*r1) m1*(J1*r2+J2*r1)*bd/(cos(theta)*(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2)); J2*m1*r1*Kv/(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2)(m1+m2)*J2*r1*Kv/(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2) m1^2*r2*Kh*r0^2/(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2)(m1+m2)*J2*Kh*r0^2/(cos(theta)*(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2)) (m1+m2)*J2*r1*Ks/(cos(theta)*(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2)) J2*m1*r1*bv/(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2)(m1+m2)*J2*r1*bv/(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2) m1^2*r2*bh*r0^2/(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2)(m1+m2)*J2*bh*r0^2*sin(theta)/(cos(theta)*(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J 2*J1*m2)) (m1+m2)*J2*r1*bd/(cos(theta)*(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2)); m1*(J1*r2+J2*r1)*Kv/(cos(theta)*(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2)) +(m1+m2)*(J1*r2+J2*r1)*Kv/(cos(theta)*(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2* J1*m2)) (m1+m2)*J1*r2*Kh*r0^2/(cos(theta)*(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1* m2)*r1)+(m1+m2)*(J1*r2+J2*r1)*Kh*r0^2/(cos(theta)^2*(m1^2*J1*r2+m1^2*J2*r1 +J2*J1*m1+J2*J1*m2)*r1) (m1+m2)*(J1*r2+J2*r1)*Ks/(cos(theta)^2*(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2* J1*m2)) m1*(J1*r2+J2*r1)*bv/(cos(theta)*(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2)) +(m1+m2)*(J1*r2+J2*r1)*bv/(cos(theta)*(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J 1*m2)) (m1+m2)*J1*r2*bh*r0^2/(cos(theta)*(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m 2)*r1)+(m1+m2)*(J1*r2+J2*r1)*bh*r0^2*sin(theta)/(cos(theta)^2*(m1^2*J1*r2+m1^ 2*J2*r1+J2*J1*m1+J2*J1*m2)*r1) - GVHD:TS L U THANH TÙNG HVTH: PHAN CƠNG BÌNH (m1+m2)*(J1*r2+J2*r1)*bd/(cos(theta)^2*(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2* J1*m2))]; B=[0; 0; 0; m1*(J1*r2+J2*r1)/(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2); (m1+m2)*J2*r1/(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2); (m1+m2)*(J1*r2+J2*r1)/(cos(theta)*(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m 2))]; C=[0 0 0]; D=0; uvao=F*(1-cos(w*t)); u=[uvao]; h=lsim(A,B,C,D,u,t); plot(t,h,'r'); grid; B ng 2B Ch ng trình mơ ph ng b ng ph n m m Matlap %Dao ng h b c t %Dao ng c a ph u %Kh o sát quan h c a thi t b VBF S va w clear; m1 = 90; m2 = 30; J1 = 2.38; GVHD:TS L U THANH TÙNG HVTH: PHAN CƠNG BÌNH J2 = 1.35; r0 = 2; r1 = 265; r2 = 265; theta = (70*pi)/180; Ks = 2.687*10^6; Kh = 393700; Kv = 1592520; bh = 3563; bv = 4783.5; bd = 0; F = 18668; f=50; w=2*pi*f; ts=0.0001; t=0:ts:10; A=[0 0; 0 0 0 0 0 0; 1; -J2*J1*Kv/(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2)m1*(J1*r2+J2*r1)*Kv/(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1* m2) m1*J1*r2*Kh*r0^2/((m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2 )*r1)m1*(J1*r2+J2*r1)*Kh*r0^2/(cos(theta)*(m1^2*J1*r2+m1^2*J2*r 1+J2*J1*m1+J2*J1*m2)*r1) m1*(J1*r2+J2*r1)*Ks/(cos(theta)*(m1^2*J1*r2+m1^2*J2*r1+J2* J1*m1+J2*J1*m2)) J2*J1*bv/(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2)m1*(J1*r2+J2*r1)*bv/(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1* GVHD:TS L U THANH TÙNG HVTH: PHAN CƠNG BÌNH m2) m1*J1*r2*bh*r0^2/((m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2 )*r1)m1*(J1*r2+J2*r1)*bh*r0^2*sin(theta)/(cos(theta)*(m1^2*J1*r 2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2)*r1) m1*(J1*r2+J2*r1)*bd/(cos(theta)*(m1^2*J1*r2+m1^2*J2*r1+J2* J1*m1+J2*J1*m2)); J2*m1*r1*Kv/(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2)(m1+m2)*J2*r1*Kv/(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2) m1^2*r2*Kh*r0^2/(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2)(m1+m2)*J2*Kh*r0^2/(cos(theta)*(m1^2*J1*r2+m1^2*J2*r1+J2*J 1*m1+J2*J1*m2)) (m1+m2)*J2*r1*Ks/(cos(theta)*(m1^2*J1*r2+m1^2*J2*r1+J2*J1* m1+J2*J1*m2)) J2*m1*r1*bv/(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2)(m1+m2)*J2*r1*bv/(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2) m1^2*r2*bh*r0^2/(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2)(m1+m2)*J2*bh*r0^2*sin(theta)/(cos(theta)*(m1^2*J1*r2+m1^2 *J2*r1+J2*J1*m1+J2*J1*m2)) (m1+m2)*J2*r1*bd/(cos(theta)*(m1^2*J1*r2+m1^2*J2*r1+J2*J1* m1+J2*J1*m2)); m1*(J1*r2+J2*r1)*Kv/(cos(theta)*(m1^2*J1*r2+m1^2*J2*r1+J2* J1*m1+J2*J1*m2))+(m1+m2)*(J1*r2+J2*r1)*Kv/(cos(theta)*(m1^ 2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2)) (m1+m2)*J1*r2*Kh*r0^2/(cos(theta)*(m1^2*J1*r2+m1^2*J2*r1+J 2*J1*m1+J2*J1*m2)*r1)+(m1+m2)*(J1*r2+J2*r1)*Kh*r0^2/(cos(t heta)^2*(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2)*r1) (m1+m2)*(J1*r2+J2*r1)*Ks/(cos(theta)^2*(m1^2*J1*r2+m1^2*J2 *r1+J2*J1*m1+J2*J1*m2)) m1*(J1*r2+J2*r1)*bv/(cos(theta)*(m1^2*J1*r2+m1^2*J2*r1+J2* J1*m1+J2*J1*m2))+(m1+m2)*(J1*r2+J2*r1)*bv/(cos(theta)*(m1^ 2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2)) (m1+m2)*J1*r2*bh*r0^2/(cos(theta)*(m1^2*J1*r2+m1^2*J2*r1+J 2*J1*m1+J2*J1*m2)*r1)+(m1+m2)*(J1*r2+J2*r1)*bh*r0^2*sin(th eta)/(cos(theta)^2*(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m 2)*r1) - GVHD:TS L U THANH TÙNG HVTH: PHAN CƠNG BÌNH (m1+m2)*(J1*r2+J2*r1)*bd/(cos(theta)^2*(m1^2*J1*r2+m1^2*J2 *r1+J2*J1*m1+J2*J1*m2))]; B=[0; 0; 0; m1*(J1*r2+J2*r1)/(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2) ; (m1+m2)*J2*r1/(m1^2*J1*r2+m1^2*J2*r1+J2*J1*m1+J2*J1*m2); (m1+m2)*(J1*r2+J2*r1)/(cos(theta)*(m1^2*J1*r2+m1^2*J2*r1+J 2*J1*m1+J2*J1*m2))]; C=[0 0 0]; D=0; f=[]; am=[]; for k=0.1:0.005:100 uvao=1700*(1-cos(k*2*pi*t)); u=[uvao]; h=lsim(A,B,C,D,u,t); tam=max(h); f=[f k]; am=[am tam]; end plot(f,am,'r'); grid; GVHD:TS L U THANH TÙNG HVTH: PHAN CƠNG BÌNH LÝ L CH TRÍCH NGANG LÝ L CH H tên: Phan Cơng Bình Ngày, tháng, n m sinh: 14/06/1982 a ch liên l c: 5/11, Bình Ph N i sinh: Bình D ng c B, xã Bình Chu n, Thu n An, Bình D ng Q TRÌNH ÀO T O T 2007-nay h c ch ng trình t o Th c s t i i h c Bách Khoa Tp H Chí Minh Q TRÍNH CƠNG TÁC N m 2005-2007 làm vi c t i nhà máy P&G Bình D ng ... ng MSHV: 00407203 I – TÊN ( TÀI: NGHIÊN C U S NH H NG C A CÁC THÔNG S THU T TRONG THI T B C P RUNG K N N NG SU T II – NHI'M V) VÀ N&I DUNG: Kh o sát nh h ng thông s k thu t thi t b c p rung -... chuyên ngành thông qua Ngày ….,tháng … n m 2009 TR NG PHÒNG T-S H TR NG KHOA QL NGÀNH iv L I CAM OAN Tôi xin cam oan ây cơng trình nghiên c u c a riêng Các s li u, k t qu nêu lu n v n trung th c... xác c so sánh v i v i nh thông s k thu t có s n ng su t VBF K t qu nghiên c u có th c s d ng thi t k VBF ho c làm tài li u cho nh ng nghiên c u ti p theo nh h ng n t i u thông s vii Trang TRANG
Ngày đăng: 08/03/2021, 18:33
Xem thêm: