Tài liệu ổn định phi tuyến của tấm có cơ tính biến thiên, không hoàn hảo

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN --------------------- Nguyễn Thị Nga ỔN ĐỊNH PHI TUYẾN CỦA TẤM CÓ CƠ TÍNH BIẾN THIÊN, KHÔNG HOÀN HẢO BẢN TÓM TẮT LUẬN VĂN THẠC SĨ KHOA HỌC Hà Nội – 2011 ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN --------------------- Nguyễn Thị Nga ỔN ĐỊNH PHI TUYẾN CỦA TẤM CÓ CƠ TÍNH BIẾN THIÊN, KHÔNG HOÀN HẢO Chuyên ngành: Cơ học vật thể rắn Mã số: 60 44 21 BẢN TÓM TẮT LUẬN VĂN THẠC SĨ KHOA HỌC NGƢỜI HƢỚNG DẪN KHOA HỌC: PGS. TS. ĐÀO VĂN DŨNG Hà Nội – 2011 Mục lục MỞ ĐẦU ………………………………………………………………………………2 CHƢƠNG 1. ỔN ĐỊNH PHI TUYẾN CỦA TẤM FGM ……………………………..5 1.1. Đặc trƣng cơ học của vật liệu FGM ………………………………………...5 1.2. Tiêu chuẩn ổn định tĩnh …………………………………………………….7 1.3. Các hệ thức cơ bản và hệ phƣơng trình ổn định ……………………………8 1.4. Phƣơng pháp giải…………………………………………………………..13 1.4.1. Phân tích ổn định của tấm chịu tải cơ ……………………………16 1.4.2. Phân tích ổn định của tấm chịu tải nhiệt………………………….17 1.4.3. Phân tích ổn định của tấm chịu đồng thời tải cơ nhiệt…………...21 CHƢƠNG 2. KHẢO SÁT BẰNG SỐ ……………………………………………... 23 2.1. So sánh với kết quả của Shen ……………………………………………..23 2.2. Khảo sát ảnh hƣởng của tỷ phần thể tích k và k1.………………………….25 2.3. Khảo sát ảnh hƣởng của điểu kiện biên……………………………………26 2.4. Khảo sát ảnh hƣởng của độ không hoàn hảo………………………………28 2.5. Khảo sát ảnh hƣởng của nhiệt độ ………………………………………….30 KẾT LUẬN …………………………………………………………………………..31 TÀI LIỆU THAM KHẢO ……………………………………………………………32 PHỤ LỤC ……………………………………………………………………………37 1 Mở đầu Vật liệu có cơ tính biến thiên (FGM) là loại vật liệu mới đang đƣợc các nhà khoa học quan tâm đặc biệt hiện nay và đƣợc sử dụng nhiều trong kỹ thuật, nhất là trong các môi trƣờng nhiệt độ cao nhƣ lò phản ứng hạt nhân và công nghiệp vũ trụ [20]. Do vậy, nghiên cứu về ổn định nhƣ ứng xử vồng và sau vồng của kết cấu FGM là một trong những mối quan tâm chính vì sự an toàn trong sử dụng và tối ƣu trong thiết kế. Javaheri và Eslami [17  19], Shariat và Eslami [34] đã nghiên cứu ổn định của các tấm chữ nhật FGM chịu tải cơ và tải nhiệt dựa trên lý thuyết tấm cổ điển [17, 19] và lý thuyết biến dạng trƣợt bậc cao [18, 34]. Họ đã sử dụng phƣơng pháp năng lƣợng và tìm đƣợc lời giải giải tích về tải tới hạn. Ảnh hƣởng của sự không hoàn hảo hình học ban đầu đến tải tới hạn đã đƣợc Shariat và nhóm nghiên cứu (NNC) [30, 33] nghiên cứu. Nhóm tác giả Najafizadeh và Eslami đã nghiên cứu sự ổn định nhiệt đàn hồi của tấm tròn có cơ tính biến thiên [27, 32]. Ma và Wang đã sử dụng lý thuyết biến dạng trƣợt bậc ba và lý thuyết tấm cổ điển để giải gần đúng bài toán uốn đối xứng trục và vồng của tấm có cơ tính biến thiên [25]. Tác giả Lanhe [21] dùng lý thuyết biến dạng trƣợt bậc nhất để dẫn ra đƣợc liên hệ giải tích của gia số nhiệt tới hạn của tấm hơi dầy, chữ nhật có cơ tính biến thiên, tựa bản lề. Phân tích tới hạn nhiệt ba chiều với cơ tính biến thiên bằng phƣơng pháp phần tử hữu hạn đã đƣợc Na và Kim đề cập đến [26]. Yang và Shen [16] đã đƣa ra phân tích phi tuyến của tấm FGM chịu tải vuông góc và tải nằm trong mặt trung bình. Một số kết quả về ứng xử sau tới hạn của tấm có cơ tính biến thiên với các phƣơng pháp khác nhau cũng đã đƣợc nghiên cứu bởi Shen [35, 36] bằng phƣơng pháp nhiễu; Zhao và Liew [39] bằng phƣơng pháp phần tử tự do kp-Ritz; Liew và NNC của ông [23, 24] bằng phƣơng pháp cầu phƣơng vi phân. Reddy và Chin 2 đã xét bài toán phân tích cơ nhiệt của trụ và tấm FGM [29]. Lý thuyết phi tuyến hình học của tấm composite lớp đẳng hƣớng ngang và sử dụng lý thuyết này vào việc phân tích sau tới hạn cũng đã đƣợc Librescu và Stein nghiên cứu [22]. Trong những năm gần đây, ở nƣớc ta đã có một số kết quả nghiên cứu quan trọng về lĩnh vực này. Các tác giả Đào Huy Bích và Lê Khả Hòa [6] nghiên cứu dao động phi tuyến của vỏ cầu thoải có cơ tính biến thiên. Đáp ứng đối xứng trục phi tuyến của vỏ cầu thoải có cơ tính biến thiên chịu áp lực ngoài có tính đến nhiệt độ đƣợc nghiên cứu bởi các tác giả Đào Huy Bích và Hoàng Văn Tùng [5]. Phân tích động lực phi tuyến của vỏ thoải không hoàn hảo có cơ tính biến thiên đƣợc các tác giả Đào Huy Bích, Vũ Đỗ Long [2] và Đào Văn Dũng, Vũ Hoài Nam [9] nghiên cứu. Các tác giả Đào Văn Dũng, Lê Khả Hòa [10] đã nghiên cứu phân tích phi tuyến sự vồng và sau vồng của panel trụ có cơ tính biến thiên bị nén dọc trục với hệ số Poisson thay đổi theo bề dày. Tuy nhiên, các kết quả tiếp cận giải tích về ổn định phi tuyến của tấm có cơ tính biến thiên chịu tải cơ nhiệt đồng thời vẫn còn ít và do vậy những nghiên cứu về vấn đề này cần đƣợc quan tâm. Gần đây, kết quả về phân tích phi tuyến về sự ổn định của tấm FGM chịu tải cơ và nhiệt đã đƣợc các tác giả Hoàng Văn Tùng và Nguyễn Đình Đức xem xét [38]. Họ đã trình bày cách tiếp cận giải tích để nghiên cứu ổn định của tấm FGM và thu đƣợc biểu thức hiển của tải vồng và đƣờng cong tải - độ võng sau vồng. Tuy nhiên các tác giả đó chỉ xét với hệ số Poisson  không đổi. Với hệ số Poisson  thay đổi theo bề dày z , cũng đã có những nghiên cứu của [4, 8, 28, 37], [11  15], chẳng hạn Fung và Chen [11] đã nghiên cứu ảnh hƣởng của tính không hoàn hảo đến sự dao động phi tuyến của tấm FGM, Navazi và NNC đã đƣa ra lời giải giải tích cho sự uốn phi tuyến của tấm FGM [28], Chen và Tan đã trình bày về phân tích độ nhạy của sự không hoàn hảo đến sự dao động phi tuyến của tấm FGM có ứng suất ban đầu [4]. Darabi và NNC phân tích ổn định động phi tuyến của vỏ trụ FGM chịu tải dọc trục 3 tuần hoàn [8]. Sofiyev và NNC đã xem xét sự ổn định của nón cụt composite ba lớp FGM chịu áp lực phân bố không đều [37]. Đặc biệt các tác giả Huang và Han đã nghiên cứu sự vồng và sau khi vồng của vỏ trụ FGM chịu tải trọng tĩnh và động [12  15] với hệ số    ( z ) theo quy luật lũy thừa của z , nhƣng các hệ số độ cứng Aij đƣợc xác định dƣới dạng tích phân, chƣa đƣợc giải tích hóa. Vì vậy, luận văn này sẽ nghiên cứu giải tích gần đúng về ứng xử vồng phi tuyến của tấm chữ nhật, không hoàn hảo có cơ tính biến thiên, kể cả hệ số Poisson  cũng thay đổi theo quy luật lũy thừa của z , chịu tải nén cơ ở mặt giữa của tấm hoặc chịu tải nhiệt hoặc chịu tải cơ - nhiệt đồng thời. Giả thiết độ chênh lệch nhiệt T không đổi. Sử dụng lý thuyết tấm cổ điển trong đó tính phi tuyến hình học theo nghĩa von Karman và phƣơng pháp Bubnov - Galerkin, để xây dựng hệ thức cho phép tìm lực tới hạn đối với ba trƣờng hợp tải tƣơng ứng ở trên. Dẫn ra phƣơng trình liên hệ tải - độ võng sau tới hạn. Tính toán các hệ số độ cứng Aij dƣới dạng giải tích hiển. Khảo sát ảnh hƣởng của điều kiện biên, tỷ phần thể tích, độ không hoàn hảo và nhiệt độ đến các đặc trƣng ổn định của tấm. 4 Chương 1 Ổn định phi tuyến của tấm FGM 1.1 Đặc trưng cơ học của vật liệu FGM Vật liệu composite là vật liệu cấu thành từ hai hay nhiều vật liệu khác nhau nhằm đạt đƣợc các tính ƣu việt nhƣ khối lƣợng nhẹ và độ bền cao, khả năng chống nhiệt, chống ăn mòn hóa học tốt, … Gần đây, một số vật liệu composite có chức năng thông minh đã ra đời nhằm đáp ứng nhu cầu thực tiễn trong việc chế tạo các kết cấu hiện đại và thỏa mãn các điều kiện làm việc khắc nghiệt nhƣ là các vật liệu gia cƣờng sợi SMA, vật liệu cơ tính biến thiên, … Vật liệu có cơ tính biến thiên là một loại vật liệu composite không thuần nhất đƣợc tạo ra từ hỗn hợp kim loại (metal) và gốm (ceramic) với các tính chất cơ học biến thiên trơn và liên tục từ mặt này đến mặt kia của lớp vật liệu. Bằng sự thay đổi liên tục về tỷ lệ thể tích của các thành phần gốm và kim loại ngƣời ta có thể thu đƣợc loại vật liệu FGM có các đặc trƣng vật liệu biến đổi dần theo bề dày. Do sự liên tục về đặc trƣng vật liệu mà FGM đã tránh đƣợc hiện tƣợng tập trung ứng suất, sự bong tách giữa các lớp nhƣ loại vật liệu composite lớp truyền thống. Vật liệu có cơ tính biến thiên thƣờng đƣợc sử dụng ở kết cấu phẳng nhẹ và chịu nhiệt độ cao. Các vật liệu cơ tính biến thiên cấu thành từ gốm và kim loại trong đó tỷ phần thể tích của mỗi phần biến đổi một cách liên tục từ mặt này sang mặt kia của kết cấu. 5 Chúng khắc phục đƣợc những nhƣợc điểm của các loại vật liệu kim loại truyền thống và composite thông thƣờng về khả năng chịu nhiệt và chịu lực. Do thành phần gốm có mođun đàn hồi 𝐸 cao và có các hệ số truyền nhiệt 𝐾, hệ số dãn nở nhiệt 𝛼 thấp làm cho vật liệu cơ tính biến thiên có độ cứng cao và khả năng chịu nhiệt tốt. Còn thành phần kim loại làm cho vật liệu FGM trở nên mềm dẻo hơn, bền hơn và khắc phục đƣợc sự rạn nứt có thể xảy ra do tính giòn của vật liệu ceramic khi chịu nhiệt độ cao. Sau đây là một số vật liệu thành phần với các tính chất cơ lý đƣợc cho bởi bảng 1.1 và 1.2. Các tính chất Vật liệu Kim loại Ceramic   kg / m3  E (N / m )  ( C ) K (W / mK ) Nhôm ( Al ) 70,0. 109 0,30 23.10−6 204 2707 Ti  6 Al  4V 105,7. 109 0,298 6,9.10−6 18,1 4429 Zirconi (ZrO2 ) 151. 109 0.30 10. 10−6 2,09 3000 Nhôm oxit 320. 109 0.26 7,2. 10−6 10,4 3750 2 o 1 Bảng 1.1. Tính chất của một số vật liệu thành phần của vật liệu có cơ tính biến thiên. Vật liệu Các tính chất E (GPa)  6  (1/ K) SUS304 207,7 0,3177 15,32.106 Si3 N 4 322,2 0,24 7, 4746.106 Bảng 1.2. Các tính chất vật liệu của tấm FGM ( T  300K ) [4, 37]. Một trong những công nghệ chế tạo của vật liệu có cơ tính biến thiên đƣợc minh họa ở mô hình sau đây: Hình 1.1. Một ví dụ về mô hình vật liệu cơ tính biến thiên FGM. Trong mô hình vật liệu FGM (hình 1.1) ta thấy mặt trên của kết cấu hoàn toàn là ceramic, mặt dƣới hoàn toàn là kim loại (metal), phần trong kết cấu là sự pha trộn giữa hai loại vật liệu theo tỷ phần thể tích. 1.2 Tiêu chuẩn ổn định tĩnh Xét trạng thái cân bằng vô cùng gần trạng thái cân bằng xuất phát (hay trạng thái cân bằng cơ bản). Với một giá trị nào đấy của lực ngoài có thể tồn tại dạng cân bằng mới đồng thời với dạng cân bằng cơ bản. Nói một cách khác, với cùng giá trị của tải có thể tồn tại nhiều dạng cân bằng khác nhau, có thể xem các dạng này là dạng 7 chuyển tiếp từ dạng cân bằng ổn định sang dạng mất ổn định. Giá trị lực ngoài nhỏ nhất để tồn tại các dạng cân bằng khác nhau gọi là tải tới hạn [1]. 1.3 Các hệ thức cơ bản và hệ phương trình ổn định Xét tấm chữ nhật có cơ tính biến thiên với các cạnh dài a , rộng b và độ dày h. Chọn hệ tọa độ Đề các vuông góc Oxyz sao cho  x, y  thuộc mặt giữa của tấm và Oz h  h theo hƣớng bề dầy    z   (Hình 1.2). 2  2 Hình 1.2. Giả thiết rằng các tính chất vật liệu của tấm thay đổi theo quy luật lũy thừa nhƣ sau [12, 14]  2z  h  k E  E  z   Em   Ec  Em     Em  Ecm r ,  2h  (1a)  2z  h  k     z    m   c  m      m   cm r ,  2h  (1b)  2z  h  k     z    m   c   m      m   cm r ,  2h  (1c) k k1 1 k trong đó 8 Ecm  Ec  Em , r  2z  h ,  cm   c  m , cm  c   m , k  0, k1  0. 2h (2) Các đại lƣợng Em , Ec ;  m ,  c ;  m ,  c là các mođun Young; hệ số Poisson; hệ số dãn nở nhiệt tƣơng ứng của kim loại và gốm. Đối với tấm không hoàn hảo, các thành phần biến dạng ở mặt giữa với tính phi tuyến hình học theo nghĩa von Karman là [3] 1 2 1  yo  v, y  w,2y  w, y w,*y , 2 o  xy  u, y  v, x  w, x w, y  w, y w,*x  w, x w,*y ,  xo  u, x  w,2x  w, x w,*x , (3) ở đây u  u  x, y  , v  v  x, y  , w  w  x, y  tƣơng ứng là các chuyển dịch theo các trục x , y, z và w*  w*  x, y  là độ không hoàn hảo ban đầu của tấm. Đại lƣợng w* đƣợc giả thiết là nhỏ. Biến dạng của điểm thuộc tấm ở khoảng cách z từ mặt giữa của tấm là  x   xo  z x ,  y   yo  z y ,  xy   xyo  2 z xy ;  x   w, xx ,  y   w, yy ,  xy   w, xy . (4) Liên hệ ứng suất - biến dạng có kể đến nhiệt độ đối với tấm đƣợc xác định bởi quy luật Hooke nhƣ sau E  x   y  (1   )T  , 1  2  E  y   x  (1   )T  , y  1  2  E  xy   xy , 2(1   ) x  (5) 9 ở đây giả thiết tấm chịu nhiệt độ thay đổi đều tức là T không đổi. Lực dãn và mômen đƣợc biểu diễn dƣới dạng N , N x y , N xy   h/2   x ,  y ,  xy dz, h / 2 M x h/2 , M y , M xy     x (6) ,  y ,  xy  zdz. h / 2 Thay các liên hệ (3)  (5) vào (6), sau khi tính toán ta nhận đƣợc  N x   A10 N    y   A20  N xy   0   M x    A11  M y   A21     M xy   0 A20 A10 0 A21 A11 0 0 0 A30 0 0 A31 A11 A21 0 A12 A22 0 A21 A11 0 A22 A12 0 0    xo  1    0    yo  2    A31    xyo   0      T ,  0    x  1  0    y  2      A32   2 xy   0  (7) ở đây các biểu thức của Aij  i  1, 2,3; j  0,1, 2  và của 1 , 2 đƣợc xác định bởi E( z) j z dz,  h / 2 1  2 ( z ) h/2 A1 j   E ( z ) ( z ) j z dz,  h / 2 1  2 ( z ) h/2 A2 j   A3 j  (8) 1 h / 2 E( z) j 1 z dz   A1 j  A2 j  , 2  h / 2 1  ( z ) 2 E ( z ) ( z ) dz,  h / 2 1  ( z ) 1    h/2 E ( z ) ( z ) zdz; j  0, 1, 2.  h / 2 1  ( z ) 2   h/2 10 Biểu thức giải tích hiển của Aij và 1 , 2 đƣợc tính toán và cho ở phần phụ lục. Phƣơng trình cân bằng của tấm không hoàn hảo có dạng [30, 31, 33] N x , x  N xy , y  0, (9a) N xy , x  N y , y  0, (9b) M x, xx  2M xy , xy  M y , yy  N x  w, xx  w,*xx   2 N xy  w, xy  w,*xy   N y  w, yy  w,*yy   0. (9c) Phƣơng trình tƣơng thích biến dạng cho tấm không hoàn hảo đƣợc viết nhƣ sau  xo, yy   yo, xx   xyo , xy  w,2xy  w, xx w, yy  2w,xy w,*xy  w,xx w,*yy  w, yy w,*xx . (10) Đƣa vào hàm ứng suất Airy f  f  x, y  sao cho N x  f, yy , N y  f, xx , N xy   f, xy . (11) Dễ dàng thấy rằng hai phƣơng trình đầu (9a), (9b) tự động thỏa mãn. Thế các liên hệ (11) vào các biểu thức (7) ta thu đƣợc  xo  J 0  A10 f, yy  A20 f, xx  J1w, xx  J 2 w, yy  J 31T  ,  yo  J 0  A10 f, xx  A20 f, yy  J1w, yy  J 2 w, xx  J 31T  , (12)  xyo   2 A31w, xy  f, xy  / A30 , trong đó 2 J 0  1/( A102  A20 ), J1  A10 A11  A20 A21, J 2  A10 A21  A20 A11, J 3  A10  A20. (13) Thế lần nữa biểu thức của  ij ở (12) vào hệ thức của momen M ij ở (7) và sau đó thay M ij vào phƣơng trình (9c) với lƣu ý (11) dẫn đến C34 f  C44 w  f, yy  w, xx  w,*xx   2 f, xy  w, xy  w,*xy   f, xx  w, yy  w,*yy   0, 11 (14) ở đây C3  J 0 J 2 ; C4  J 0  A11J1  A21J 2   A12 ; 4   4 / x4  2 4 / x2y 2   4 / y 4 . Phƣơng trình (14) này có hai hàm cần tìm là w và f , do vậy cần phải tìm thêm phƣơng trình thứ hai liên quan đến hai hàm này, bằng cách sử dụng phƣơng trình tƣơng thích biến dạng (10). Muốn vậy, thế các biểu thức ở trên của  ij trong (12) vào phƣơng trình (10), thu đƣợc 4 f  C14 w  C2  w,2xy  w, xx w, yy 2w, xy w,*xy  w, xx w,*yy  w, yy w,*xx   0, (15) trong đó C1  J 2 / A10 , C2  1/ J 0 A10. Hai phƣơng trình (14) và (15) là các phƣơng trình chủ đạo dùng để nghiên cứu sự ổn định phi tuyến của tấm có cơ tính biến thiên. Trƣờng hợp w*  0 , từ (14) và (15), ta đƣợc các phƣơng trình ổn định cơ bản của tấm hoàn hảo có cơ tính biến thiên. Lƣu ý rằng các phƣơng trình (12), (14) và (15) tƣơng tự nhƣ các phƣơng trình trong bài báo [12], nhƣng các hệ số Aij và 1 , 2 ở đây đã đƣợc giải tích hóa dƣới dạng hiển. 1.4 Phương pháp giải Giả thiết xét ba trƣờng hợp điều kiện biên sau đối với tấm chữ nhật nhƣ sau [22, 38] 12 Trường hợp (1). Tất cả các cạnh của tấm tựa bản lề và có thể di chuyển tự do (FM) trong mặt phẳng của tấm, tức là w  M xx  N xy  0, N x  N x 0 tại x  0, x  a, (16) w  M yy  N xy  0, N y  N y 0 tại y  0, y  b. Trường hợp (2). Tất cả các cạnh của tấm tựa bản lề nhƣng không dịch chuyển đƣợc (IM) trong mặt phẳng của tấm, tức là w  u  M xx  0, N x  N x 0 tại x  0, x  a, w  v  M yy  0, N y  N y 0 tại y  0, y  b. (17) Trường hợp (3). Tất cả các cạnh của tấm tựa bản lề. Tải tác động theo trục x . Hai cạnh x  0, x  a, có thể di chuyển đƣợc trong mặt phẳng của tấm, còn hai cạnh y  0, y  b. thì không di chuyển đƣợc. Điều kiện biên trong trƣờng hợp này là w  M xx  N xy  0, N x  N x 0 tại x  0, x  a, (18) w  v  M yy  0, N y  N y 0 tại y  0, y  b. ở đây N x 0 , N y 0 là các lực trƣớc vồng tƣơng ứng theo hƣớng x và y đối với trƣờng hợp (1) và (3a) và chúng là các phản lực trên các cạnh không đƣợc dịch chuyển của tấm (trƣờng hợp (2) và trƣờng hợp (3b)). Điều kiện biên ở trên sẽ đƣợc thỏa mãn đối với w và theo nghĩa trung bình đối với f nếu nhƣ độ võng w và hàm ứng suất f đƣợc biểu diễn bằng [22, 38] w  W sin m x sin n y (19a) 13 m n , n  , m, n  1,2,... là số bán sóng tƣơng ứng theo hƣớng x và a b y, còn W là biên độ của độ võng. trong đó m  * * Liên quan đến độ không hoàn hảo ban đầu w  w  x, y  , ta giả thiết nó có dạng giống nhƣ độ lệch w của tấm, tức là w*  w*  x, y    h sin m x sin n y; m, n  1, 2, ... (19b) ở đây hệ số  [0,1] biểu thị cỡ của độ không hoàn hảo của tấm. Bằng cách thay các biểu thức (19a) và (19b) vào phƣơng trình (15), ta tìm đƣợc 2 1 4 f  C1W  m2  n2  sin m x sin n y  C2W m2 n2 W  2 h  cos 2m x 2 1  C2W m2 n2 W  2 h  cos 2n y  0 . 2 Nghiệm của phƣơng trình này đƣợc tìm dƣới dạng f  A cos 2m x  B cos 2n y  C sin m x sin n y  1 1 N x0 y 2  N y0 x2. 2 2 (20) trong đó A, B, C là các hệ số cần xác định. Thay (19a, b), (20) vào vế trái phƣơng trình (15) ta đƣợc C2 n2W W  2 h  A , 32m2 C2m2W W  2 h  B , 32 n2 (21) C  C1W . Thay (19a, b), (20) vào vế trái phƣơng trình (14) và sau đó áp dụng phƣơng pháp Bubnov - Galerkin, ta nhận đƣợc phƣơng trình để xác định tải vồng và xây dựng 14 đƣờng cong sau vồng của tấm chữ nhật có cơ tính biến thiên chịu tải cơ, tải nhiệt và tổ hợp tải cơ nhiệt.  C4  C1C3   m2  n2  2 W   N x 0m2  N y 0 n2  W   h   32 C1m n W W   h  3 ab C2  m4  n4  16 C2C3m2 n2  W W  2 h   W W   h W  2 h   0 3 mn  16 (22) ở đây m và n là các số lẻ. Nhận xét:  Nếu hệ số Poisson  là hằng số, thì k1  0 , C1  0, C3  0, C2  E1 , C4   D . Phƣơng trình (22) sẽ trùng với phƣơng trình đã đƣợc cho trong bài báo [38].  Nếu hệ số không hoàn hảo   0 thì phƣơng trình (22) cho ta phƣơng trình dùng để tìm tải tới hạn của tấm hoàn hảo có cơ tính biến thiên. Sau đây, chúng ta sẽ xét một cách chi tiết ba bài toán tƣơng ứng với ba dạng tải trọng đã nêu ở trên. 1.4.1 Phân tích ổn định của tấm chịu tải cơ Xét tấm chữ nhật không hoàn hảo, có cơ tính biến thiên tựa bản lề tại bốn cạnh và các cạnh có thể di chuyển tự do trong mặt phẳng của tấm (trƣờng hợp (1)) và tấm chịu các lực nén trong mặt phẳng Px và Py phân bố đều dọc theo các cạnh x  0, x  a, và y  0, y  b tƣơng ứng. Trong trƣờng hợp này các lực nén trong giai đoạn trƣớc tới hạn cho bởi [3, 17] 15 N x 0   Px h, N y 0   Py h. (23) Thế biểu thức này vào phƣơng trình (22) và đặt W  W / h,   Py / Px , m / n  Ba m / n, Ba  b / a , ta nhận đƣợc Px      C 4  C1 C 3  m 2 Ba2  n 2  B m B   n 2 h 2 2 a 2 2 W W        16 C 2 C 3 Ba2 mn W W  2 3 Bh2  m 2 Ba2   n 2  W     C 2  m 4 Ba4  n 4  C1Ba2 mn 32  W W W  2 3 Bh2  m 2 Ba2   n 2  16 Bh2  m 2 Ba2   n 2      (24) ở đây C1  C C1 C C b , C 2  2 , C 3  3 , C 4  34 , Bh  . 2 h h h h h Phƣơng trình (24) có thể dùng để vẽ đƣờng cong tải - độ võng sau vồng của tấm có cơ tính biến thiên chịu nén bởi tải cơ phẳng. Với   const, tức là k1  0 thì C1  0, C2  E1 , C3  0, C4   D, phƣơng trình (24) dẫn tới phƣơng trình đã đƣợc thiết lập của bài báo [38]. Với tấm hoàn hảo,   0 , từ phƣơng trình (24) ta thu đƣợc phƣơng trình cho phép xác định tải nén tới hạn nhƣ sau Px      C 4  C1 C 3  m2 Ba2  n 2   Bh2  m2 Ba2   n 2    C 2  m4 Ba4  n 4  16 Bh2  m2 Ba2   n 2  2 16 C 2 C 3 B 2 mn  C1Ba2 mn 32 a W   2 2 2 2 2 2 2 2  3 Bh  m Ba   n  3 Bh  m Ba   n   2 W . (25) Trƣờng hợp   const và W  0 , phƣơng trình (25) trở thành phƣơng trình đã đƣợc trình bày bởi Javaheri và Eslami [17] 16 Pxb    D  m2 Ba2  n 2  2 Bh2  m 2 Ba2   n 2  . (26) 1.4.2 Phân tích ổn định của tấm chịu tải nhiệt Giả thiết tấm không hoàn hảo tựa bản lề với bốn cạnh không dịch chuyển đƣợc trong mặt phẳng (trƣờng hợp (2)). Khi đó điều kiện tại các cạnh x  0, x  a và y  0, y  b sẽ đƣợc thỏa mãn theo nghĩa trung bình nhƣ sau [22, 36] u dxdy  0,  x 0 0 b a a b v  y dydx  0. (27) 0 0 Để giải bài toán này, trƣớc hết ta cần phải xác định N x 0 , N y 0 sao cho đảm bảo các cạnh của tấm không dịch chuyển đƣợc trong mặt phẳng. Các đại lƣợng N x 0 , N y 0 này đƣợc tìm từ hệ phƣơng trình cân bằng sau N x , x  N xy , y  0, (28) N xy , x  N y , y  0, và điều kiện (27). Vì đã đƣa vào hàm ứng suất Airy f  x, y  ở (11), nên hệ phƣơng trình (28) đƣợc tự động thỏa mãn. Do vậy chỉ còn điều kiện (27) có thể dùng để tìm N x 0 , N y 0 . Nhằm mục đích này, từ các liên hệ (3) và (12) ta thu đƣợc các phƣơng trình sau 1 u, x  J 0  A10 f, yy  A20 f, xx  J1w, xx  J 2 w, yy  J 31T   w,2x  w, x w,*x , 2 1 v, y  J 0  A10 f, xx  A20 f, yy  J1w, yy  J 2 w, xx  J 31T   w,2y  w, y w,*y . 2 17 (29) Thế biểu thức của w , w* từ (19a, b) và f từ (20) vào (29), và sau đó vào điều kiện (27), cho ta N x 0  1T  N y0 4 A11 1 ( A10m2  A20 n2 )W  ( A10m2  A20 n2 )W W  2 h  , 2 mn A10 8 4 A11 1  1T  ( A10 n2  A20m2 )W  ( A10 n2  A20m2 )W W  2 h  . 2 mn A10 8 (30) Nếu   const thì biểu thức (30) dẫn đến kết quả của [38]. Nếu   const và W  0 thì biểu thức (30) trùng với kết quả của Javaheri và Eslami [3], tức là N x 0  N y 0  1T  m . 1  (31) Thay liên hệ (30) vào phƣơng trình (22), ta có 4 4 4 2 2 2  mnBa   4 A11  A10  m Ba  n   2 A20 m n Ba  32 C1 1T     W 2 2 2  2 2 2 2 A 3 ab ( m B  n ) mnb m B  n   10 a  a  1 2  2 A10  C2   m 4 Ba4  n 4   4 A20 m 2 n 2 Ba2  W W  2 h   2 2 2 2   16 b  m Ba  n   2  m2 Ba2  n2  W W W  2 h  16 C2C3mnBa2    C4  C1C3  . 3 b 2  m2 Ba2  n 2  W h b2 W h (32) Mặt khác tham số nhiệt 1 của (8) đƣợc cho bởi (xem phần phụ lục)  1 1 an   Em cm  Ecm m   an n0 nk1  1 n0 k  nk1  1 1  h  Em m    1 an   hP . n0 2k  nk1  1   Ecm cm  (33) 18