ổn định của dầm liên tục và của dàn

10 790 4
ổn định của dầm liên tục và của dàn

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

Thông tin tài liệu

Trong chương này chúng ta sẽ nghiên cứu các phương pháp xác định lực tới hạn cho các dầm liên tục và các thanh trong dàn phẳng.

Chương 5. Ổn định của dầm liên tục của dàn 5-1 Chương 5. ỔN ĐỊNH CỦA DẦM LIÊN TỤC CỦA DÀN Trong chương này chúng ta sẽ nghiên cứu các phương pháp xác định lực tới hạn cho các dầm liên tục các thanh trong dàn phẳng. Những bài toán này thường được đưa về bài toán ổn định của thanh đơn giản, thanh liên tục trên gối đàn hồi hoặc thanh làm việc trong môi trường đàn hồi. 5.1 Cách tính ổn định của dầm liên tục theo phương pháp lực, phương trình ba mô men Ta vận dụng phương pháp lực đã nghiên cứu trong chương 4 để tính ổn định của dầm liên tục. Giả sử dầm có tiết diện không đổi trong từng nhịp chịu lực dọc trục đặt ở các gối tựa (hình 5-1a). Ta chọn hệ cơ bản như trên (hình 5-1b) thì phương trình chính tắc viết cho gối thứ i bất kỳ sẽ chỉ có ba số hạng: 0MδMδMδ1)(i1)i(iiii1)(i1)i(i=++++−− (5-1) Đó là phương trình ba mô men. Các hệ số của phương trình được xác định theo các trạng thái đơn vị như trên (hình 5-2). Những dầm đơn giản chịu tải trọng như trên (hình 5-2) đã được nghiên cứu trong mục 2, chương 4. Theo công thức (4-5) ta có: kP1i+1i2ki-1kik Pi+1k PPPl lk1P kP2k Pi-1k Pii+1k Pi-1MMiMi+1Pa,b,Hình 5-1. Sơ đồ tính ổn định dầm liên tục. i+1ii-1 M M iM i+1P Pl ll M =1P ii-1 δi(i-1)ll iM ii+1βiαi+1i(i+1)lδi+1M i+1δ =α +βii ii+1 Hình 5-2. Sơ đồ tính theo ba mô men. Chương 5. Ổn định của dầm liên tục của dàn 5-2 ()ii1)i(ivβ6EJδil=− (5-2) () ()1i1i1ii1iiiivα3EJvα3EJδ+++++=+=iillϕϕ (5-3) ()1i1i1i1)i(ivβ6EJδ++++=l (5-4) Trong đó các hàm số α(vi) β(vi) được xác định theo công thức sau: ()⎟⎟⎠⎞⎜⎜⎝⎛−=ii2iitgvv1v3vα (5-5) ()⎟⎟⎠⎞⎜⎜⎝⎛−= 1sinvvv6vβii2ii (5-6) Với iiiEJPkvil= (5-7) Thay (5-2), (5-3) vào (5-1) ta được: () () ( )[ ]( )0MvβλMvαλvαλ2Mvβλ1i1i1ii1i1iii1iii=++++++++− (5-8) Đó là phương trình ba mô men khi tính ổn định của dầm liên tục chịu tác dụng của lực dọc trục. Trong đó: i0iJJil=λ (5-9) Với J0 là đại lượng bất kỳ, thường lấy bằng mô men quán tính của một nhịp nào đó. Trình tự tính toán theo các bước sau: 1- Xác định các chiều dài quy ước λi theo (5-9). 2- Xác định các đại lượng vi theo (5-7). Trong trường hợp vi có các giá trị khác nhau ta cần biểu thị các vi theo một đại lượng vk nào đó theo biểu thức sau: kiikkkii vJEkJEkvkill= (5-10) Trong đó: kkkkk JEPkv l.= (5-11) 3- Thiết lập các phương trình ba mô men theo (5-8). Dầm có n gối tựa trung gian ta sẽ lập được n phương trình. Hệ phương trình này là thuần nhất. 4- Thiết lập phương trình ổn định bằng cách cho định thức các hệ số của hệ phương trình bằng không. Chương 5. Ổn định của dầm liên tục của dàn 5-3 5- Giải phương trình ổn định ta sẽ tìm được nghiệm vk từ (5-11) suy ra lực tới hạn cần tìm. Lực tới hạn tìm được ở trên tương ứng với trường hợp dầm bị mất ổn định với các mô men ở gối tựa khác không. Trong bài toán dầm liên tục, ngoài nghiệm tìm được theo cách trình bày kể trên ta còn phải tìm nghiệm tương ứng với trường hợp dầm bị mất ổn định với các mô men gối tựa M1 = M2 = M3 = . = Mn = 0, nếu về ý nghĩa vật lý, trường hợp này có thể xảy ra. Nếu dầm liên tục có gối bên trái là ngàm thì hệ số δ11 được xác định như sau: () ()22211111vEJlvEJlαγδ+=. Trong đó: ()()v41vtgv2v2vtg.v2tgvvcosv)v(sinvvsinv2cosv2v2ϕ=−−=−−−=γ (5-12) Bảng các hàm số ϕ2(v) cho trong phần phụ lục. Do đó phương trình ba mô men thứ nhất có dạng: () ()[]()0MvMv2v622212211=++ βλγλγλ (5-13) Nếu dầm liên tục có gối bên phải là ngàm thì () ()1n1n1nnnnnnv3EJv3EJ++++= γαδll Do đó, phương trình ba mô men cuối cùng có dạng: ()()()0Mv6v2Mvn1n1nnn1nnn=⎥⎦⎤⎢⎣⎡+++++−γλαλβλ (5-14) Các phương trình viết cho những gối tựa khác vẫn có dạng (5-8). 5.2 Cách tính ổn định của dầm liên tục theo phương pháp chuyển vị Nội dung phương pháp chuyển vị đã trình bày trong chương 4. Hệ cơ bản chọn như trên (hình 5-3a b). Theo (4-27), phương trình chính tắc biểu thị điều kiện phản lực mô men tại liên kết đặt thêm vào thứ k bằng không có dạng: 0ZrZrZr1k1)k(kkkk1k1)k(k=++++−−, (5-15) k+1kl la,lln+11Z1k-1ZkZk+1ZnZb,l11Zk-1Z Zk+1kZ ZnHình 5-3. Sơ đồ tính dầm theo PP. chuyển vị. Chương 5. Ổn định của dầm liên tục của dàn 5-4 trong đó k biến thiên từ 1 đến n. Phương trình (5-15) chỉ gồm ba số hạng nên được gọi là phương trình ba góc xoay. Các hệ số của phương trình này được xác định theo các công thức sau: ()k31)k(kv2EJrϕkkl=− (5-16) () ()1k21k1kk2kkkkv4EJv4EJr++++=ϕϕll (5-17) ()1k31k1k1)k(kv2EJr+++−=ϕl (5-18) Công thức (5-17) chỉ nghiệm đúng với trường hợp k > 1 k < n. Khi k = 1 k = n, biểu thức của hệ số rkk phụ thuộc vào điều kiện liên kết ở các đầu dầm. Nếu các gối biên là ngàm như trên (hình 5-3a) thì các hệ số r11 rnn được xác định theo công thức (5-17). Nếu các gối biên là khớp tựa như trên (hình 5-3b) thì r11 rnn được xác định theo công thức sau: () ()2222111111v4EJv3EJrϕϕll+= (5-19) () ()1n11n1nn2nnnnv3EJv4EJr++++=ϕϕll (5-20) Theo các phương trình thuần nhất (5-15) ta có thể thiết lập phương trình ổn định của dầm theo phương pháp chuyển vị bằng cách cho định thức các hệ số của hệ phương trình đó bằng không: D = 0. (5-21) Sau khi khai triển định thức D giải phương trình (5-21) ta dễ dàng tìm được lực tới hạn cho dầm liên tục. Quá trình thực hiện tương tự như đã trình bày trong mục 5, chương 4. Lực tới hạn tìm được theo điều kiện (5-21) xảy ra tương ứng với trường hợp dầm bị mất ổn định trong đó các chuyển vị xoay Zk khác không. Trong thực tế có thể xảy ra trường hợp dầm bị mất ổn định với các chuyển vị xoay Zk bằng không. Bởi vậy, khi nghiên cứu ổn định của dầm liên tục ta cần phải xét cả 2 trường hợp có thể xảy ra. Chẳng hạn đối với dầm liên tục hai nhịp, phương trình chính tắc có dạng: () ()0Zv3EJv3EJ111=⎥⎦⎤⎢⎣⎡+ .ϕϕll. Phương trình này có thể thoả mãn với Z1 = 0 tương ứng với dạng mất ổn định đối xứng cũng có thể thoả mãn với Z1 ≠ 0 tương ứng với dạng mất ổn định phản đối xứng. Trong trường hợp Z1 ≠ 0, ta có: ()()0vtgv3tgvvv21=−=ϕ. 5.3 Ổn định của các thanh chịu nén trong dàn Dưới tác dụng của tải trọng, các thanh chịu nén của dàn bị có thể bị mất ổn định làm cho toàn dàn bị phá hoại. Những thanh chịu nén trong dàn có thể là: Chương 5. Ổn định của dầm liên tục của dàn 5-5 1. Các thanh đứng, thanh biên hoặc thanh xiên không cắt qua các thanh khác, thí dụ như thanh AB, AC CD trên (hình 5-5). Để kiểm tra ổn định, ta coi thanh là thanh đơn giản có liên kết khớp ở hai đầu, sau đó ta có thể tính theo công thức đã nghiên cứu trong mục 2 hoặc mục 8 của chương 2. (Giả thiết thanh có khớp ở hai đầu chỉ là gần đúng). 2. Những thanh đứng hoặc thanh xiên cắt qua một, hai hoặc nhiều thanh đứng hoặc thanh xiên khác. Trên (hình 5-4) trình bày một số thí dụ về những thanh thuộc loại này. Các thanh chịu nén ACB trên (hình 5-4avà b) là thanh xiên cắt qua một thanh xiên khác ở giữa nhịp. Các thanh ACDB trên (hình 5-4c d) là thanh xiên cắt qua hai thanh đứng hoặc hai thanh xiên khác của dàn. Khi mất ổn định, những thanh này làm việc giống như những thanh đặt trên hai khớp tựa cứng ở hai đầu có một, hai hoặc nhiều gối tựa đàn hồi ở trong nhịp (hình 5-6). Như vậy, bài toán ổn định của những loại thanh này được đưa về bài toán ổn định của thanh liên tục có các gối tựa trung gian là gối tựa đàn hồi. Cách giải quyết bài toán này sẽ nghiên cứu trong mục 4. Khi số lượng các gối tựa đàn hồi trung gian tương đối lớn ta có thể giải quyết bài toán này theo trường hợp thanh làm việc trong môi trường đàn hồi. ABCDHình 5-5. Sơ đồ dàn. BFA1JJCEACBACBDABCe qhfDJJ1a,b,c,d, Hình 5-4. Sơ đồ các dạng dàn . a), b)- các thanh xiên cắt thanh xiên khác ở giữa nhịp; c), d)- thanh xiên cắt qua hai thanh đứng hoặc hai thanh xiên khác của dàn. Chng 5. n nh ca dm liờn tc v ca dn 5-6 3. H thanh biờn trờn ca cu dn h tc l cu dn khụng cú ging giú phớa trờn. H thanh biờn ny chu lc nộn thay i dc theo chiu di thanh. Khi mt n nh, h thanh biờn trờn b cong ra ngoi mt phng dn, cỏc khung ngang trong cu gm dm ngang ca b phn mt cu v thanh ng ngn cn khụng cho phộp h thanh biờn chuyn v t do v lm vic ging nh nhng liờn kt n hi. Kinh nghim cho bit, khi s t ca dn ln hn bn thỡ ta cú th thay cỏc gi n hi bng nn n hi. Nh vy, bi toỏn ny c a v trng hp thanh lm vic trong mụi trng n hi chu lc nộn thay i dc theo chiu di thanh. 5.4 n nh ca thanh liờn tc cú gi ta n hi 5.4.1 n nh ca thanh liờn tc hai nhp cú gi trung gian l gi n hi (hỡnh 5-6). Gi c l cng ca liờn kt n hi. cng c chớnh l phn lc cn tỏc dng ti liờn kt n hi sao cho liờn kt bin dng vi giỏ tr bng n v. gii bi toỏn ny theo phng phỏp chuyn v ta cú h c bn nh trờn (hỡnh 5-7a). Vỡ h ang xột cú tớnh cht i xng nờn ta cú th phõn tớch bi toỏn thnh hai trng hp: Thanh b mt n nh theo dng i xng (ng I trờn hỡnh 5-6) v thanh b mt n nh theo dng phn i xng (ng II trờn hỡnh 5-6). 5.4.1.1 Trng hp thanh b mt n nh theo dng i xng. Trong trng hp ny ta cú Z1 0, Z2 = 0. l/2IACBl/2IIHỡnh 5-6. S thanh ZEJ=const21Z3i (v)11l (v)3i2Z3i (v)1a,b,c,M12M Hỡnh 5-7.H c bnvcỏcbiu mụ men043214 8 12 16P/PƠlecl32EJ Hỡnh 5-8. S bin thiờn ca t s P/Ple vi cng ca gi n hi. Chương 5. Ổn định của dầm liên tục của dàn 5-7 Phương trình chính tắc có dạng: 0Zr111=. Phương trình ổn định: r11 = 0. Từ biểu đồ đơn vị 1M vẽ trên (hình 5-7b) ta dễ dàng xác định được: ()0cv3i2.r12111=+= ηl, trong đó: EJPEJPv121ll ==. Suy ra: ()48EJcv31l−=η. (5-22) Như vậy, nếu biết độ cứng của liên kết đàn hồi thì ta có thể sử dụng bảng 2 trong phần phụ lục để xác định thông số v từ đó suy ra lực tới hạn. 5.4.1.2 Trường hợp thanh bị mất ổn định theo dạng phản đối xứng. Trong trương hợp này ta có Z1 = 0; Z2 ≠ 0. Phương trình chính tắc có dạng: 0Zr222=. Phương trình ổn định: r22 = 0. Từ biểu đồ 2M vẽ trên (hình 5-7c) ta dễ dàng xác định được phản lực đơn vị r22: r22 = 2.3iϕ1(v) = 0. Suy ra: ϕ1(v) = 0. Phương trình này thoả mãn với v = π. Do đó, ta có: 22212thEJ4EJPllππ==. (5-24) Sau khi xác định lực tới hạn tương ứng với hai trường hợp biến dạng nói trên ta sẽ chọn giá trị nhỏ nhất làm giá trị tới hạn. Trên (hình 5-8) là đồ thị biểu diễn luật biến thiên của tỷ số lực tới hạn với lực Ơle theo các độ cứng khác nhau của gối đàn hồi. Ta thấy: Khi c < 16π2EJ/l2 thì thanh sẽ mất ổn định theo dạng đối xứng, còn khi c > 16π2EJ/l2 thì thanh sẽ bị mất ổn định theo dạng phản đối xứng, lúc này lực tới hạn được xác định theo công thức (5-24) không phụ thuộc độ cứng c. Để nghiên cứu sự ổn định của các thanh chịu nén ACB trong dàn (hình 5-4a b) ta có thể áp dụng lời giải vừa tìm được ở trên. Gọi EJ là độ cứng của thanh chịu nén ACB đang khảo sát; EJ1 là độ cứng của thanh bị cắt ECF. Muốn áp dụng kết quả trên ta cần Chương 5. Ổn định của dầm liên tục của dàn 5-8 xác định độ cứng c của liên kết đàn hồi. Trong trường hợp này độ cứng c có giá trị bằng lực cần tác dụng tại giữa nhịp của thanh đơn giản ECF để sao cho điểm C chuyển vị theo phương vuông góc với mặt phẳng dàn có giá trị bằng đơn vị, ta có: 1δ12==48EJcl. Suy ra 31l48EJc = (5-25) Thay (5-25) vào (5-22) ta được phương trình ổn định tương ứng với dạng mất ổn định đối xứng: ()JJv11−=η (5-26) Nếu biết tỷ số J1/J thì ta có thể tìm được thông số v từ đó suy ra lực tới hạn. Có thể biểu thị lực tới hạn của thanh theo công thức quen thuộc như sau: ()2µπlEJP2th= (5-27) Hệ số µ phụ thuộc tỷ số J1/J có các giá trị cho trong bảng 5-1. Bảng 5-1. Bảng các giá trị của hệ số m J1/J 0,1 0,2 0,4 0,5 0,6 0,8 1,0 2,0 3,0 π2/3 µ 0,950 0,912 0,845 0,818 0,793 0,750 0,714 0,586 0,516 0,500 Nếu xét dạng mất ổn định phản đối xứng thì theo (5-24) ta có µ = 0,5. Theo (5-26), trị số này tương ứng với khi J1/J = π2/3 = 3,2898. Như vậy khi J1/J < π2/3 thì thanh sẽ bị mất ổn định theo dạng đối xứng còn khi J1/J > π2/3 thì thanh sẽ bị mất ổn định theo dạng phản đối xứng hệ số µ có giá trị không đổi bằng 0,5. Trên (hình 5-9) là đồ thị biến thiên của hệ số µ theo các tỷ số J1/J khác nhau. 5.4.2 Ổn định của thanh liên tục ba nhịp có gối trung gian là gối đàn hồi Xét hệ cho trên (hình 5-10). Giả sử độ cứng c của hai gối trung gian như nhau. Hệ đang xét có tính chất đối xứng nên ta có thể phân tích bài toán thành hai trường hợp: biến dạng đối xứng biến dạng phản đối xứng. Pl l l000PIIIlHình 5-10. Sơ đồ thanh 00,80,60,40,21 2 3,29J1J1,00,53µ0,7140,5860,5160,5Hình 5-9. Biến thiên hệ số m với J1/J. Chương 5. Ổn định của dầm liên tục của dàn 5-9 5.4.2.1 Trường hợp hệ bị mất ổn định theo dạng đối xứng. Ta có sơ đồ tính có dạng như trên (hình 5-11a). - Hệ cơ bản (hình 5-11b). - Phương trình ổn định: 0rrrrD22211211== - Vẽ biểu đồ đơn vị (hình 5-11c d) - Xác định được phản lực đơn vị : () ()⎥⎦⎤⎢⎣⎡+=+= fv3EJcv3EJr111152ηη3030ll; ()v3EJrr12112ϕ30l==; () ()⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡+=+=2vtg2vv3EJ2vtg2v2EJv3EJr112232000ϕϕlll; trong đó: EJPEJPv30ll ==; (5-28) EJc6EJ5cf3ll.162530=−= (5-29) Sau khi thay phản lực đơn vị vào phương trình ổn định khai triển ra ta được: Z2Z11lϕ (v)3EJa, b, c, 1Z =102d, Z =123EJ2l0ϕ (v)10l2EJ v/2tgv/2ll /200M21MP P P P Hình 5-11. Mất ổn định theo dạng đối xứng. Chương 5. Ổn định của dầm liên tục của dàn 5-10() () ()()2312312511121vtgvvvtgvvvvf+⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡+−=ϕϕηϕ (5-30) Nếu biết kích thước của thanh độ cứng của liên kết đàn hồi, tức là biết giá trị của đại lượng f thì sau khi giải phương trình (5-30) ta sẽ xác định được thông số v từ đó suy ra lực tới hạn. Theo (5-28) ta có: ()222thEJEJ9vPllµπ2==, Trong đó: 3vπµ =. 5.4.2.2 Trường hợp thanh bị mất ổn định theo dạng phản đối xứng. Sơ đồ tính của thanh có dạng như trên (hình 5-12a). Hệ cơ bản vẽ trên (hình 5-12b). Phương trình ổn định: 0rrrrD22211211==. Từ các biểu đồ đơn vị vẽ trên hình 5-12c, d ta dễ dàng xác định được các phản lực đơn vị: () ()⎥⎦⎤⎢⎣⎡+⎟⎠⎞⎜⎝⎛+=+⎟⎠⎞⎜⎝⎛+= f522vη8vη3EJc2vη24EJvη3EJr113013013011lll ; ()⎥⎦⎤⎢⎣⎡⎟⎠⎞⎜⎝⎛+=2vv3EJr1122ϕϕ20l; ()⎥⎦⎤⎢⎣⎡⎟⎠⎞⎜⎝⎛−==2vv3EJrr112112ϕϕ430l. Thay các số liệu vừa tìm được ở trên vào phương trình ổn định sau khi khai triển ta được: () () ()()⎟⎠⎞⎜⎝⎛+⎥⎦⎤⎢⎣⎡⎟⎠⎞⎜⎝⎛+⎥⎦⎤⎢⎣⎡⎟⎠⎞⎜⎝⎛+−⎥⎦⎤⎢⎣⎡⎟⎠⎞⎜⎝⎛−=2v2v2v2v2v8v2v4v25f111111211ϕϕϕϕϕϕηη (5-32) trong đó f v được xác định theo các công thức (5-29) (5-28). Tương tự như trên, nếu biết f độ cứng của gối đàn hồi thì sau khi giải (5-32) ta sẽ xác định được v từ đó suy ra lực tới hạn theo (5-31). Z2Z1b, c, 1Z =1d, Z =123EJ2l0ϕ (v)1ll /200M21MP P P P ϕ (v/2)12EJ2l10ϕ (v)120l3EJϕ (v/2)16EJl20a, Hình 5-12. Mất ổn định theo dạng phản xứng [...]...Chương 5. Ổn định của dầm liên tục của dàn 5-10 () () () () 2 3 1 2 3 1 2 5 1 11 2 1 v tg v v v tg v vvv f + ⎥ ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎢ ⎣ ⎡ +− = ϕ ϕηϕ (5-30) Nếu biết kích thước của thanh độ cứng của liên kết đàn hồi, tức là biết giá trị của đại lượng f thì sau khi giải phương trình (5-30) ta sẽ xác định được thơng số v từ đó suy ra lực tới hạn. Theo (5-28)... trên vào phương trình ổn định sau khi khai triển ta được: () () () () ⎟ ⎠ ⎞ ⎜ ⎝ ⎛ + ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ ⎟ ⎠ ⎞ ⎜ ⎝ ⎛ + ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ ⎟ ⎠ ⎞ ⎜ ⎝ ⎛ +− ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ ⎟ ⎠ ⎞ ⎜ ⎝ ⎛ − = 2 v 2v 2 v 2v 2 v 8v 2 v 4v 2 5 f 11 1111 2 11 ϕϕ ϕϕϕϕ ηη (5-32) trong đó f v được xác định theo các công thức (5-29) (5-28). Tương tự như trên, nếu biết f độ cứng của gối đàn hồi thì sau khi giải (5-32) ta sẽ xác định. .. () 2 22 th EJEJ9v P l l µ π 2 == , Trong đó: 3v π µ = . 5.4.2.2 Trường hợp thanh bị mất ổn định theo dạng phản đối xứng. Sơ đồ tính của thanh có dạng như trên (hình 5-12a). Hệ cơ bản vẽ trên (hình 5- 12b). Phương trình ổn định: 0 rr rr D 2221 1211 == . Từ các biểu đồ đơn vị vẽ trên hình 5- 12c, d ta dễ dàng xác định được các phản lực đơn vị: () () ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ + ⎟ ⎠ ⎞ ⎜ ⎝ ⎛ +=+ ⎟ ⎠ ⎞ ⎜ ⎝ ⎛ +=... như trên, nếu biết f độ cứng của gối đàn hồi thì sau khi giải (5-32) ta sẽ xác định được v từ đó suy ra lực tới hạn theo (5-31). Z 2 Z 1 b, c, 1 Z =1 d, Z =1 2 3EJ 2 l 0 ϕ (v) 1 l l /2 0 0 M 2 1 M P P P P ϕ (v/2) 12EJ 2 l 1 0 ϕ (v) 1 2 0 l 3EJ ϕ (v/2) 1 6EJ l 2 0 a, Hình 5-12. Mất ổn định theo dạng phản xứng . 5. Ổn định của dầm liên tục và của dàn 5-1 Chương 5. ỔN ĐỊNH CỦA DẦM LIÊN TỤC VÀ CỦA DÀN Trong chương này chúng ta sẽ nghiên cứu các phương pháp xác định. hệ số của hệ phương trình bằng không. Chương 5. Ổn định của dầm liên tục và của dàn 5-3 5- Giải phương trình ổn định ta sẽ tìm được nghiệm vk và từ

Ngày đăng: 17/10/2012, 11:30

Từ khóa liên quan

Tài liệu cùng người dùng

Tài liệu liên quan