Phân tích ổn định phi tuyến của vỏ trụ sandwich cơ tính biến thiên có gân gia cường trực giao chịu tải xoắn trong môi trường nhiệt độ

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - NGUYỄN TIẾN THIỆP PHÂN TÍCH ỔN ĐỊNH PHI TUYẾN CỦA VỎ TRỤ SANDWICH CƠ TÍNH BIẾN THIÊN CĨ GÂN GIA CƢỜNG TRỰC GIAO CHỊU XOẮN TRONG MƠI TRƢỜNG NHIỆT ĐỘ LUẬN VĂN THẠC SĨ KHOA HỌC Hà Nội – 2018 ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - NGUYỄN TIẾN THIỆP PHÂN TÍCH ỔN ĐỊNH PHI TUYẾN CỦA VỎ TRỤ SANDWICH CƠ TÍNH BIẾN THIÊN CĨ GÂN GIA CƢỜNG TRỰC GIAO CHỊU XOẮN TRONG MÔI TRƢỜNG NHIỆT ĐỘ Chuyên ngành: Cơ học vật rắn Mã số: LUẬN VĂN THẠC SĨ KHOA HỌC NGƯỜI HƯỚNG DẪN KHOA HỌC: TS LÊ KHẢ HÕA Hà Nội – Năm 2018 LỜI CẢM ƠN Trước trình bày nội dung luận văn thạc sỹ khoa học này, em xin bày tỏ lòng biết ơn sâu sắc đến q Thầy Cơ khoa tốn – Cơ – Tin học, trường Đại học Khoa học tự nhiên – ĐHQGHN với tri thức tâm huyết để truyền đạt kiến thức quý báu cho chúng em suốt thời gian học tập trường Em xin chân thành cảm ơn TS Lê Khả Hịa tận tâm hướng dẫn để em hoàn thành luận văn Em xin chân thành cảm ơn ban lãnh đạo khoa Toán – Cơ- Tin học, phòng sau Đại học- Trường Đại học Khoa Học Tự Nhiên- ĐHQGHN Qua xin chân thành cảm ơn tới gia đình, bạn bè, đồng nghiệp ln ủng hộ suốt q trình học tập Bước đầu nghiên cứu, kiến thức em nhiều hạn chế Do vậy, em không tránh khỏi thiếu sót mong nhận đóng góp quý giá quý thầy cô bạn đồng môn Sau em xin kính chúc TS Lê Khả Hịa tồn thể q thầy khoa Tốn – Cơ – Tin học thật nhiều sức khỏe, niềm vui để tiếp tục cống hiến cho khoa học mang hết tâm huyết truyền đam mê khoa học cho hệ học trò mai sau Hà Nội, ngày 19 tháng 10 năm 2018 Học viên Nguyễn Tiến Thiệp Mục lục Danh sách hình vẽ i Danh sách bảng vẽ ii MỞ ĐẦU Chương – Vỏ trụ hệ thức .5 1.1 Vỏ trụ sandwich đàn hồi Pasternak Các hệ thức phương trình cân Chương – Phương pháp giải 14 2.1 Điều kiện biên dạng nghiệm 14 Chương – Tính tốn số .21 3.1 So sánh kết 21 3.2 Kết số cho vỏ nón cụt ES – FGM 21 KẾT LUẬN 30 TÀI LIỆU THAM KHẢO .32 Danh sách hình vẽ Hình 1: Vỏ trụ sandwich có gân gia cường Hình Ảnh hưởng hcore/hFG lên tải xoắn tới hạn τupper (k2=k3=N=k=1) 23 Hình Ảnh hưởng hcore/hFG lên tải xoắn tới hạn τlower (k2=k3=N=k=1) 23 Hình 4b Ảnh hưởng hcore/hFG lên đường cong τ – ψ (Case 1: k2=k3=N=k=1) 24 Hình 4a Ảnh hưởng hcore/hFG lên đường cong τ - Wmax/h (Case k2=k3=N=k=1) Hình 5a Ảnh hưởng e0 lên đường cong τ - Wmax/h (Case 1: k2=k3=N=k=1) Hình 5b Ảnh hưởng e0 lên đường cong τ – ψ (Case 1: k2=k3=N=k=1) Hình Ảnh hưởng e0 lên tải xoắn tới h Hình Ảnh hưởng ΔT lên đường cong Hình 8a Ảnh hưởng số tỉ phần thể tích k lên đường cong τ - Wmax/h (Case 1) 26 Hình 8b Ảnh hưởng số tỉ phần thể tích k lên đường cong τ – ψ (Case 1) 26 Hình Ảnh hưởng số tỉ phần thể tích k lên tải xoắn tới hạn τupper Hình 10 Ảnh hưởng R/h lên tải xoắn tới hạn τupper (Case , k=1 Hình 11a Ảnh hưởng tỉ số R/h lên đường cong τ - Wmax/h (Case k2=k3=N=k=1) Hình 11b Ảnh hưởng tỉ số R/h lên đường cong τ – ψ (Case 1: k2=k3=N=k=1) 27 Hình 12a Ảnh hưởng L/R lên đường cong τ - Wmax/h (Case 1: k2=k3=N=k=1) 28 Hình 12b Ảnh hưởng tỉ số L/R lên đường cong τ – ψ (Case 1: k2=k3=N=k=1) 28 Hình 13a Ảnh hưởng gân lên đường cong τ - Wmax/h (Case 1: k2=k3=N=k=1) 29 Hình 13b Ảnh hưởng gân lên đường cong τ – ψ (Case 1: k2=k3=N=k=1) 29 i Danh sách bảng Bảng 1: So sánh tải tới hạn xoắn (psi) với [30] [29] vỏ trụ không gân chịu xoắn 21 Bảng 2: So sánh tải tới hạn xoắn (psi) với [29] [31] vỏ trụ không gân chịu xoắn 21 Bảng :Ảnh hưởng hệ số độ xốp e0 độ dầy lớp lõi hcore đến tải tới hạn τcr, k2=k3=N=k=1, h=0.006m, L/R=1.5, R/h=80, ΔT=300K, hs=hr=0.006m, bs=br=0.006m, ns=20, nr=20, K1=6×107 N/m3, K2=4×105 N/m 22 Bảng 4: Ảnh hưởng nhiệt độ số tỉ phần thể tích đến tải xoắn tới hạn tải tới hạn h=0.006m, L/R=1.5, R/h=80, a/hFG = 2, e0=0.5, hs=hr=0.006m, bs=br=0.006m, ns=20, nr=20, , K1=6×107 N/m3, K2=4×105 N/m 25 Bảng 5: Ảnh hưởng thơng số hình học đến tải xoắn tới hạn tải xoắn tới hạn k2=k3=N=k=1, h=0.006m, a/hFG = 2, e0=0.5, hs=hr=0.006m, bs=br=0.006m, ns=20, nr=20, , K1=6×107 N/m3, K2=4×105 N/m, ΔT=300K 27 Bảng 6: Ảnh hưởng gân đến tải xoắn tới hạn k2=k3=N=k=1, h=0.006m, L/R=1.5, R/h=80, a/hFG = 2, e0=0.5, hs=hr=0.006m, bs=br=0.006m, ns=20, nr=20, K1=6×107N/m3, K2=4×105 N/m, ΔT=300K (Trường hợp 1: gân gia cưởng phía trong) 28 ii MỞ ĐẦU Các kết cấu chế tạo từ vật liệu tính biến thiên (Functionally graded Material-FGM) sử dụng ngày nhiều công nghiệp hàng khơng vũ trụ, lị phản ứng hạt nhân lĩnh vực làm việc môi trường nhiệt độ cao chịu tải phức tạp Do tính chất lý biến đổi trơn liên tục từ mặt đến mặt nên kết cấu FGM tránh tập trung ứng suất bề mặt tiếp xúc lớp, tránh bong tách rạn nứt kết cấu Do nghiên cứu ổn định, dao động độ bền kết cấu FGM thu hút ý đặc biệt cộng đồng nhà khoa học nước Một kết cấu sử dụng nhiều thực tế vỏ trụ Đã có nhiều kết nghiên cứu toán ổn định dao động vỏ trụ làm vật liệu FGM không gân chịu tải xoắn, tác giả Batra [1] nghiên cứu toán xoắn vỏ trụ với modul vật liệu thay đổi dọc theo trục Tác giả Wang cộng [2] đưa nghiệm xác đáp ứng dao động vỏ trụ FGM rỗng hữu hạn Các tác giả Zhang Fu [3] nghiên cứu trạng thái tới hạn vỏ trụ đàn hồi có lớp phủ cứng bề mặt biến dạng lớp lõi bề mặt phân tích thơng qua phương trình Navier mơ hình lớp vỏ mỏng, tương ứng Các tác giả Sofiyev Schnack [4] nghiên cứu ổn định vỏ trụ FGM chịu tải xoắn động tăng tuyến tính Trong đó, họ sử dụng phương trình ổn định động dạng Donnell cải tiến phương pháp Galerkin với liên hệ hình học tuyến tính nghiệm gần lựa chọn số hạng Tác giả Li Wang [5] nghiên cứu ổn định đàn hồi vỏ trụ FGM sandwich với điều kiện biên tựa đơn chịu tải xoắn phương pháp bán giải tích Các phương trình chủ đạo với thành phần chuyển vị toán ổn định kết cấu thiết lập sử dụng lý thuyết vỏ mỏng Flugge quan hệ biến dạng – chuyển vị tuyến tính Về vỏ đàn hồi, Dựa lý thuyết biến dạng trượt bậc cao tác giả Bagherizadeh cộng [6] nghiên cứu ổn định vỏ trụ FGM không gân bao quanh đàn hồi Pasternak Sử dụng lý thuyết vỏ Donnell cải tiến, nghiên cứu ổn định phân tích dao động vỏ trụ FGM khơng gân đàn hồi tác giả Sofiyev cộng nghiên cứu [7] Bằng phương pháp Galerkin, tác giả đưa biểu thức xác định tải xoắn tới hạn tần số dao động vỏ Về tốn ổn định vỏ có gân gia cƣờng: Năm 2009, tác giả Najafizadeh cộng [8] với phương trình ổn định tuyến tính theo chuyển vị, nghiên cứu trạng thái tới hạn vỏ trụ FGM gia cường gân dọc gân vòng chịu nén dọc trục Gân vỏ vật liệu có tính biến thiên Năm 2011, tác giả Đào Huy Bích cộng [9] trình bày cách tiếp cận giải tích để nghiên cứu ổn định tĩnh phi tuyến FGM vỏ thoải FGM có gân gia cường lệch tâm dựa lý thuyết vỏ cổ điển gân giả thiết Dựa vào lý thuyết vỏ cổ điển có tính tới yếu tố phi tuyến hình học von Karman sử dụng kỹ thuật san tác dụng gân, Dung Hoa [10] nghiên cứu ổn định phi tuyến đáp ứng sau tới hạn vỏ trụ có gân FGM gia cường chịu tải xoắn đàn hồi Pasternak môi trường nhiệt độ Vỏ gia cường gân vịng gân dọc vỏ gân làm vật liệu tính biến thiên Với phương pháp tương tự, Thang Trung [11] nghiên cứu ảnh hưởng gân đến ổn định phi tuyến vỏ trụ với lớp phủ FGM chịu tải xoắn Vỏ trụ gia cường gân dọc trục Sử dụng lý thuyết vỏ cổ điển, Thang cộng [12] trình bày đáp ứng vỏ trụ có lớp phủ vật liệu FGM hai lớp phủ làm vật liệu FGM lớp lõi làm vật liệu Vỏ gia cường gân dọc nằm phía ngồi vỏ Với hàm độ võng ba số hạng, Ninh Bich [13] nghiên cứu ổn định vỏ trụ với ba lớp gốm – FGM – kim loại (C-FGM-M) chịu tải xoắn đàn hồi Vỏ gia cường gân dọc gân vịng đặt mơi trường nhiệt độ Về vỏ làm việc môi trƣờng nhiệt độ: Dựa lý thuyết biến dạng trượt bậc cao, tác giả Shen [14] thu kết toán ổn định vỏ trụ FGM chịu xoắn môi trường nhiệt độ Trong báo ông sử dụng phương pháp tham số bé để xác định tải tới hạn đường cong sau tới hạn Bằng cách sử dụng lý thuyết biến dạng lớn phi tuyến phương pháp Ritz, tác giả Huang Han [15] nghiên cứu ổn định phi tuyến vỏ trụ không gân chịu tải xoắn môi trường nhiệt Trong nghiên cứu họ, thành phần độ võng phi tuyến đề xuất dựa vào quan sát từ thực nghiệm Vật liệu xốp (porous material) loại vật liệu vật liệu FGM tính chất xốp vật liệu đặc trưng phân bố lỗ rỗng bên Các vật liệu xốp bọt kim loại (metal foams) loại vật liệu nhẹ quan trọng với khả hấp thụ lượng tuyệt vời Do loại vật liệu nhận ứng dụng rộng rãi thực tế cho cấu trúc chịu tải động Trong báo [16], tác giả Magnucki Stasiewicz nghiên cứu ổn định đàn hồi dầm xốp Magnucka-Blandzi [17, 18, 19] trình bày mơ hình tốn học sandwich hình chữ nhật hình trịn với phần lõi kim loại bọt Các nghiên cứu đưa biểu thức xác định tải tới hạn Phân tích ổn định đàn hồi dầm xốp chịu uốn trình bày tác giả Chen nhóm nghiên cứu [20] Tác giả Kitipornchai nhóm nghiên cứu nghiên cứu dao động tự ổn định đàn hồi dầm xốp báo [21] Sử dụng lý thuyết biến dạng trượt bậc cao, tác giả Jabbari cộng [22, 23] trình bày ổn định phi tuyến tròn làm vật liệu xốp Sử dụng phương pháp Galerkin, tác giả Tú cộng [24] nghiên cứu ổn định phi tuyến chịu nén Bài báo đưa biểu thức biểu thức xác định tải tới hạn mô tả đường cong sau tới hạn Các kết cấu làm vật liệu xốp gần chủ đề quan tâm đặc biệt nhà thiết kế xây dựng Để tăng cường khả làm việc kết cấu người ta thường giảm trọng lượng kết cấu cách sử dụng vật liệu xốp Cách làm có ưu điểm giảm trọng lượng cơng trình, tiết kiệm vật liệu đặc biệt giảm ảnh hưởng thành phần lực quán tính xuất kết cấu có ảnh hưởng tải trọng động Vì vậy, việc nghiên cứu kết cấu làm vật liệu xốp khơng tối ưu vật liệu mà cịn tối ưu giá thành thể tích k tăng Giá trị tải xoắn tới hạn với k=0 (kim loại) nhỏ Tính chất phù hợp với thực tế vật liệu giá trị k lớn tương ứng với hàm lượng gốm lớp phủ tăng lên dẫn tới tính kháng nhiệt vỏ độ cứng kết cấu vỏ trụ tăng lên 25 τ τ (MPa) (MPa) h= 0.006m, R/h=80, L/R=1.5 hcore/hFG = 2, e0=0.5, ΔT=300K ψ (độ) Hình 8b Ảnh hưởng số tỉ phần thể tích k lên đường cong τ – ψ (Case 1) Hình 8a Ảnh hưởng số tỉ phần thể tích k lên đường cong τ - Wmax/h (Case 1) h= 0.006m, hcore/hFG=2 e0=0.5, ΔT=300K τupper (MPa) upper (3) (4) hs=hr=0.0 06m bs=br=0.0 06m (1 ) (2) ns=nr=20 Nền τ hcore/hFG = 2, e0=0.5 (MPa) Gân h= 0.006m, R/h=80, L/R=1.5 Nền K1=6×107N/m3 K2=4×10 k Hình Ảnh hưởng số tỉ phần thể tích k lên tải xoắn tới hạn τupper d) K1=6 ×10 N/m3 K 2= 4×10 N/ m R/h Hình 10 Ảnh hưởng R/h lên tải xoắn tới hạn τupper (Case 1, k=1) Ảnh hưởng thơng số hình học (ảnh hưởng tỉ số R/h L/R) Bảng biểu diễn ảnh hưởng tỉ số R/h L/R lên tải xoắn tới hạn vỏ trụ xốp có gân gia cường trực giao đàn hồi môi trường nhiệt độ (ΔT=300K) Các hình 10, 11 12 ảnh hưởng tỉ số R/h, L/R lên đường cong mô tả quan hệ τupper − R / h , τ −Wmax / h τ −ψ Các ảnh hưởng tỉ số R/h lên đường cong τ −Wmax / h τ −ψ biểu diễn hình 11 Các kết thu khả chịu lực vỏ giảm đáng kể mà tỉ số R/h tăng 26 Bảng 5: với tỉ số R/h, tải xoắn tới hạn giảm mà tỉ số L/R tăng Ví dụ, với R/h=80 (trường hợp 1), tải xoắn tới hạn  lower =1105.4533 MPa tương ứng với L/R=1 tương ứng lớn tải xoắn tới hạn tỉ số L/R=1.5, khoảng 1.24, 1.38 lần Bảng 5: Ảnh hưởng thơng số hình học đến tải xoắn tới hạn tải xoắn tới hạn k2=k3=N=k=1, h=0.006m, a/hFG = 2, e0=0.5, hs=hr=0.006m, bs=br=0.006m, ns=20, nr=20, , K1=6×107 N/m3, K2=4×105 N/m, ΔT=300K τ (MPa) L/R=1 Trường L/R=1.5 hợp L/R=2 L/R=1 Trường L/R=1.5 hợp L/R=2 a Mode vồng (n, λ); m=1 h= 0.006m, L/R=1.5, hcore/hFG =2 e0=0.5, ΔT=300K h= 0.006m, L/R=1.5 τ (MPa) -0.15 -0.24 Gân hs=hr=0.006m bs=br=0.006m ns=nr=20 (1) (2) (3) -0.29 (2) (3) (4) 952 MPa 741 MPa (4) Nền K1=6×107N/m3 K2=4×105N/m -0.53 (1) 1735 MPa τ (MPa) hcore/hFG = 2, e0 = 0.5, ΔT=300K Wmax/h 387 MN Gân Nền K1=6×107N/m3 K2=4×105N/m ψ (độ) Hình 11a Ảnh hưởng tỉ số R/h lên đường cong τ - Wmax/h (Case 1: k2=k3=N=k=1) Hình 11b Ảnh hưởng tỉ số R/h lên đường cong τ – ψ (Case 1: k2=k3=N=k=1) hs=hr=0.006m bs=br=0.006m ns=nr=20 27 h= 0.006m, R/h=80 hcore/hFG = 2, e0=0.5 ΔT=300K h= 0.006m, R/h=80, hcore/hFG = e0=0.5, ΔT=300K τ(MPa) τ (MPa) (1) 1169 MPa (2) (3) 952 MPa 845 MPa (2) (3) Nền Gân hs=hr=0.006m bs=br=0.006m ns=nr=20 K1=6×107N/m3 K2=4×10 N/m Wmax Hì nh 12 a Ả nh hư ởn g củ a L/ R lê n đư ờn g co ng τW ma x/ h (C as e 1: k2 = k3 = N = k= 1) ψ (độ) H ì n h b Ả n h h n g c ủ a t ỉ s ố L / R l ê n đ n g c o n g τ – ψ (Case 1: k2=k 3=N=k=1) Bảng 6: Ảnh hưởng gân đến tải xoắn tới hạn k2=k3=N=k=1, h=0.006m, L/R=1.5, R/h=80, a/hFG = 2, e0=0.5, hs=hr=0.006m, bs=br=0.006m, ns=20, nr=20, K1=6×107N/m3, K2=4×105 N/m, ΔT=300K (Trường hợp 1: gân gia cưởng phía trong) τ (MPa) K1=3×107 N/m3 K1=6×107 N/m3 K1=9×107 N/m3 28 h= 0.006m, R/h=80, L/R=1.5 h= 0.006m, R/h=80, L/R=1.5 hcore/hFG = 2, e0=0.5, ΔT=300K hcore/hFG = 2, e0=0.5, ΔT=300K τ (MPa) τ(MPa) Gân hs=hr=0.006m bs=br=0.006m 952 MPa 857 MPa (1) 714 MPa (2) ns=nr=20 (3) (4) 609 MN ψ (độ) H ì n h c o n g τ a - Ả n h / h h n g W m a x ( C a s e : k c ủ a g â n v n ề n l ê n đ n g = k = N = k = ) Hình 13b Ảnh hưởng gân lên đường cong τ – ψ (Case 1: k2=k3=N=k =1) e) Ảnh hưởng gân Ảnh hưởng gân đến tải xoắn tới hạn tải xoắn tới hạn vỏ trụ sandwich xốp có gân gia cường mơ tả bảng Như quan sát thấy, tải xoắn tới hạn tăng hệ số hền K1 K2 tăng Điều hình 13 Thêm vào đó, bảng với số gân gia cường vỏ xốp gia cường gân vịng có khả chịu xoắn tốt nhất, đứng thứ hai gân trực giao, sau gân dọc Vỏ khơng gia cường (khơng có gân gia cường) khơng có đàn hồi có tải xoắn tới hạn nhỏ 29 KẾT LUẬN Luận văn trình bày phương pháp giải tích để nghiên cứu tốn ổn định vỏ trụ sandwich xốp có gân gia cường lệch tâm Phương trình cân vỏ trụ nhận dựa vào lý thuyết vỏ Donnell kỹ thuật san tác dụng gân; sử dụng phương pháp Galerkin để tìm phương trình tính tải xoắn tới hạn phương trình khảo sát đáp ứng sau tới hạn Các tính tốn số ảnh hưởng thông số độ xốp, độ dầy lớp lõi, số gân, nhiệt độ, số tỉ phần thể tích k tham số hình học đến tải xoắn tới hạn đáp ứng sau tới hạn vỏ trụ sandwich xốp có gân gia cường Luận văn thu số kết sau: i) Đã trình bày chi tiết kết tính tốn giải tích tìm phương trình hiển xác định tải xoắn tới hạn phương trình khảo sát đáp ứng sau tới hạn vỏ trụ Sandwich xốp có gân gia cường có đàn hồi mơi trường nhiệt độ ii) Khảo sát số nghiên cứu ảnh hưởng tham số đến tải xoắn tới hạn đường cong sau tới hạn vỏ trụ Sandwich xốp có gân gia cường từ rút kết luận sau: - Tải xoắn tới hạn τcr vỏ trụ xốp có gân gia cường giảm hệ số độ xốp e0 tăng - Khi hcore/hFG tăng tải xoắn tới hạn hai trường hợp gân gân giảm - Khả chịu tải xoắn vỏ trụ giảm nhiệt độ ΔT tăng - Giá trị tải tới hạn tăng số tỉ phần thể tích k lớp phủ tăng - Khi tỉ số R/h L/R tăng lên khả mang tải vỏ trụ xốp có gân gia cường trực giao đàn hồi giảm 30 - Sự có mặt gân cách bố trí gân thay đổi ảnh hưởng đáng kể đến khả mang tải vỏ Với vỏ có gân gia cường khả mang tải vỏ cao hẳn so với vỏ khơng có gân gia cường Với việc gắn gân vịng khả mang tải trường hợp cao so với trường hợp gắn gân dọc gân trực giao Hƣớng nghiên cứu luận văn: - Giải toán vỏ trụ sandwich xốp có gân gia cường chịu tải động - Giải kết cấu vỏ khác vỏ nón, panel cầu làm vật liệu xốp - Sử dụng lý thuyết biến dạng trượt để nghiên cứu toán 31 TÀI LIỆU THAM KHẢO [1] Batra R.C (2006), “Torsion of a functionally graded cylinder”, AIAA J 44, pp 1363–1365 [2] Wang H.M., Liu C.B., Ding H.J (2009), “Exact solution and transient behavior for torsional vibration of functionally graded finite hollow cylinders”, Acta Mech Sinica 25, pp 555–563 [3] Zhang P., Fu Y (2011), “Torsional buckling of elastic cylinders with hard coatings”, Acta Mech 220, pp 275–287 [4] Sofiyev A.H., Schnack E (2004), “The stability of functionally graded cylindrical shells under linearly increasing dynamic torsional loading”, Eng Struct 26, pp 1321–1331 [5] Li S, Wang S Elastic stability of circular cylindrical shells with FGMS layers under end torsion force Lixue Xuebao/Chinese Journal of Theoretical and Applied Mechanics 2010; 42 (6): 1172-79 [6] Bagherizadeh E., Kiani Y., Eslami M.R (2011), “Mechanical buckling of functionally graded material cylindrical shells surrounded by Pasternak elastic foundation”, Compos Struct 93, pp 3063-3071 [7] Sofiyev A.H., Kuruoglu N (2013), “Torsional vibration and buckling of the cylindrical shell with functionally graded coatings surrounded by an elastic medium”, Compos Part B: Eng 45, pp 1133–1142 [8] Najafizadeh M.M., Hasani A., Khazaeinejad P (2009), “Mechanical stability of functionally graded stiffened cylindrical shells”, Appl Math Modelling 33, pp 1151-1157 [9] Bich D.H., Nam V.H., Phuong N.T (2011), “Nonlinear post-buckling of eccentrically stiffened functionally graded plates and shallow shells” Vietnam J Mech 3, pp.131-147 32 [10] Dung DV, Hoa LK Nonlinear torsional buckling and postbuckling of eccentrically stiffened FGM cylindrical shells in thermal environment Composites: Part B 69 (2015) 378–388 [11] Thang PT, Trung NT Effect of stiffeners on nonlinear buckling of cylindrical shells with functionally graded coatings under torsional load Composite Structures 153 (2016) 654–661 [12] Thang PT, Duc ND, Trung NT Thermomechanical buckling and postbuckling of cylindrical shell with functionally graded coatings and reinforced by stringers Aerospace Science and Technology 66 (2017) 392–401 [13] Ninh DG, Bich DH Nonlinear torsional buckling and postbuckling of eccentrically stiffened ceramic functionally graded material metal layer cylindrical shell surrounded by elastic foundation subjected to thermo-mechanical load Journal of Sandwich Structures and Materials, 18 (6), (2016), pp 712-738 [14] Shen HS Torsional buckling and postbuckling of FGM cylindrical shells in thermal environments International Journal of Non-Linear Mechanics 2009;44:644-57 [15] Huang H, Han Q Nonlinear buckling of torsion-loaded functionally graded cylindrical shells in thermal environment European Journal of Mechanics-A/Solids 2010;29:42-8 [16] Magnucki, K and P Stasiewicz, Elastic buckling of a porous beam Journal of Theoretical and Applied Mechanics, 2004 42(4): p 859868 [17] Magnucka Blandzi, E Mathematical Modeling of a Rectangular Sandwich Plate with a Non - Homogeneous Core in AIP Conference Proceedings 2007 AIP [18] Magnucka-Blandzi, E., Mathematical modelling of a rectangular sandwich plate with a metal foam core Journal of Theoretical and Applied Mechanics, 2011 49(2): p 439-455 33 [19] Magnucka-Blandzi, E., Axi-symmetrical deflection and buckling of circular porous-cellular plate Thin-walled structures, 2008 46(3): p 333-337 [20] Chen, D., J Yang, and S Kitipornchai, Elastic buckling and static bending of shear deformable functionally graded porous beam Composite Structures, 2015 133: p 54-61 [21] Kitipornchai, S., D Chen, and J Yang, Free vibration and elastic buckling of functionally graded porous beams reinforced by graphene platelets Materials & Design, 2017 116: p 656-665 [22] Jabbari, M., et al., Buckling analysis of a functionally graded thin circular plate made of saturated porous materials Journal of Engineering Mechanics, 2013 140(2): p 287-295 [23] Mojahedin, A., et al., Buckling analysis of functionally graded circular plates made of saturated porous materials based on higher order shear deformation theory Thin-Walled Structures, 2016 99: p 83-90 [24] Tu TM, Hoa LK, Hung DX and Hai LT Nonlinear buckling and postbuckling analysis of imperfect porous plates under mechanical loads Journal of Sandwich Structures & Materials https://doi.org/10.1177/1099636218789612 [25] M.R Feyzi, A.R Khorshidvand Axisymmetric post-buckling behavior of saturated porous circular plates Thin-Walled Structures 112 (2017) 149–158 [26] Brush DO, Almroth BO Buckling of bars, plates and shells Mc Graw-Hill, New York; 1975 [27] Reddy JN Mechanics of laminated composite plates and shells: Theory and Analysis Boca Raton: CRS Press 2004 [28] Volmir AS Stability of elastic systems Science Edition Moscow; 1963 (in Russian) 34 [29] Shen H.S (2009), “Torsional buckling and postbuckling of FGM cylindrical shells in thermal environments”, Int J Non-Linear Mech 44, pp 644-657 [30] Nash W.A (1959), “An experimental analysis of the buckling of thin initially imperfect cylindrical shells subject to torsion”, Proceedings of the Society for Experimental stress Analysis 16(2), pp 55-68 [31] Ekstrom R.E (1963), “Buckling of cylindrical shells under combined torsion and hydrostatic pressure”, Experimental Mechanics 3, pp 192-197 35 ... nghiên cứu ổn định phi tuyến đáp ứng sau tới hạn vỏ trụ có gân FGM gia cường chịu tải xoắn đàn hồi Pasternak môi trường nhiệt độ Vỏ gia cường gân vòng gân dọc vỏ gân làm vật liệu tính biến thiên Với... HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - NGUYỄN TIẾN THIỆP PHÂN TÍCH ỔN ĐỊNH PHI TUYẾN CỦA VỎ TRỤ SANDWICH CƠ TÍNH BIẾN THIÊN CĨ GÂN GIA CƢỜNG TRỰC GIAO CHỊU XOẮN TRONG MƠI... dọc gân vòng đặt môi trường nhiệt độ Về vỏ làm việc môi trƣờng nhiệt độ: Dựa lý thuyết biến dạng trượt bậc cao, tác giả Shen [14] thu kết toán ổn định vỏ trụ FGM chịu xoắn môi trường nhiệt độ Trong