(LUẬN văn THẠC sĩ) đánh giá trên và dưới hệ số dẫn nhiệt vật liệu tổ hợp đẳng hướng

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC CÔNG NGHỆ VƯƠNG THỊ MỸ HẠNH ĐÁNH GIÁ TRÊN VÀ DƯỚI HỆ SỐ DẪN NHIỆT VẬT LIỆU TỔ HỢP ĐẲNG HƯỚNG LUẬN VĂN THẠC SỸ CƠ HỌC KỸ THUẬT Hà Nội-2012 TIEU LUAN MOI download : skknchat@gmail.com ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC CÔNG NGHỆ VƯƠNG THỊ MỸ HẠNH ĐÁNH GIÁ TRÊN VÀ DƯỚI HỆ SỐ DẪN NHIỆT VẬT LIỆU TỔ HỢP ĐẲNG HƯỚNG Ngành : Cơ học kỹ thuật Chuyên ngành: Cơ học kỹ thuật Mã số: 60 52 02 LUẬN VĂN THẠC SỸ CƠ HỌC KỸ THUẬT Hà Nội-2012 TIEU LUAN MOI download : skknchat@gmail.com MỤC LỤC Trang Trang phụ bìa ……………………………………………………………… LỜI CẢM ƠN …………………………………………………………… LỜI CAM ĐOAN …………………………………………………………… MỤC LỤC ………………………………………………………………… DANH MỤC CÁC BẢNG VÀ HÌNH VẼ …………………………………… Chương Mở đầu ………………………………………………………… Chương Đánh giá cho hệ số dẫn nhiệt ……………………… 11 2.1 Xây dựng đánh giá ……………………………………… 15 2.2 Xây dựng đánh giá ……………………………………… 20 Chương Áp dụng ………………………………………………………… 23 3.1 Vật liệu tựa đối xứng …………………………………………… 23 3.2 Mơ hình cầu lồng ba pha …………………………… 28 KẾT LUẬN ………………………………………………………………… 35 DANH MỤC CƠNG TRÌNH KHOA HỌC ĐÃ CƠNG BỐ …………………… 36 TÀI LIỆU THAM KHẢO ………………………………………………… 37 PHỤ LỤC ………………………………………………………………… 40 TIEU LUAN MOI download : skknchat@gmail.com DANH MỤC CÁC BẢNG VÀ HÌNH VẼ Trang Hình 3.1: Mơ hình vật liệu tựa đối xứng ba pha ……………………………… 23 Bảng 3.1: Đánh giá hệ số dẫn hiệu qủa vật liệu tựa đối xứng ba pha ………………………………………………………….………… 27 Hình 3.2: Đánh giá hệ số dẫn hiệu vật liệu tựa đối xứng ba pha ……………………………………… …………………… 28 Hình 3.3: Mơ hình cầu lồng hai pha …………………………… … 29 Hình 3.4: Mơ hình cầu lồng ba pha ………………………… ….… 32 Bảng 3.2: Đánh giá hệ sơ dẫn hiệu cho mơ hình cầu lồng ba pha …………………………………………….… … 33 Hình 3.5: Đánh giá hệ sô dẫn hiệu cho mơ hình cầu lồng ba pha ……………………………………… ………… 34 TIEU LUAN MOI download : skknchat@gmail.com Chương Mở đầu Một số lớn vật liệu sử dụng tạo từ nhiều thành phần vật liệu khác nhằm phục vụ cho đòi hỏi nhiều lĩnh vực đời sống người Trong luận văn ta quan tâm nghiên cứu hệ số dẫn nhiệt vĩ mô vật liệu tổ hợp đẳng hướng, phụ thuộc vào tính chất thành phần cấu thành tương tác chúng Việc nghiên cứu mối quan hệ cần thiết có ý nghĩa thực tiễn giúp giải thích mối quan hệ tính chất vĩ mơ vật liệu với tính chất thành phần cấu thành hình học vi mơ, giúp thiết kế vật liệu với tính chất vĩ mô theo yêu cầu Vật liệu tổ hợp cấu tạo vi mô từ cách thành phần vật liệu khác mặt vĩ mô coi đồng tính chất hữu hiệu ( mơ đun đàn hồi, hệ số dẫn nhiệt, điện,…) nói chung khác với tính chất thành phần cấu thành Các cấu trúc vi mô coi đủ lớn so với kích thước phan tử để xem môi trường liên tục Các trường (nội lực, chuyển vị, nhiệt, dòng nhiệt…) liên tục mặt ngăn cách pha Khi thành phần cấu thành phân bố hỗn độn không gian ta có vật liệu tổ hợp đẳng hướng vĩ mơ Ở giới hạn với giả thiết vật liệu thành phần đẳng hướng Việc xác định lý tính hữu hiệu vật liệu tổ hợp vấn đề khoa học vật liệu Hệ số dẫn nhiệt vĩ mô vật liệu tổ hợp đẳng hướng phụ thuộc vào tính chất vi mô vật liệu thành phần cấu tạo nên như: tỷ lệ thể tích, dạng hình học pha, hệ số dẫn nhiệt vi mơ… Vì hình học pha vật liệu thường phức tạp cho trước đầy đủ xác được, nên người ta thường tìm cách xây dựng đánh giá tính chất vĩ mơ vật liệu Trong chương này, muốn giới thiệu đến bạn đọc lịch sử vấn đề, tầm quan trọng ý nghĩa thực tiễn nó, đồng thời sơ lược qua phương pháp nghiên cứu kết mà luận văn đạt Việc xác định hệ số dẫn nhiệt hữu hiệu vật liệu tổ hợp đẳng hướng đưa việc giải tốn tìm điểm cực tri phiếm hàm lượng Đó tồn phức tạp tương đương với việc giải phương trình cân TIEU LUAN MOI download : skknchat@gmail.com tương thích trêm miền khơng đồng V với cấu tạo hình học pha phức tạp, tùy tiện mà thực tế người ta khó mà có thơng tin đầy đủ hình học pha vật liệu tổ hợp sử dụng Do đó, đường lối biến phân đề nghị cách tiếp cận khác thực tiễn hơn: sở thơng tin hạn chế có cấu tạo miền thể tích xét ta tìm trường tốt gần với với điểm cực trị phiếm hàm lượng để từ nhận đánh giá hệ số dẫn nhiệt hiệu Khó khăn ta phải giải tốn biến phân miền thể tích V với cấu trúc phức tạp mà ta thường khơng có đầy đủ thơng tin Các thơng tin xác nhất, đơn giản cho mơ tả (và quan trọng nhất) có thường tỷ lệ thể tích vα tính chất dẫn nhiệt Cα thành phần cấu thành Vào năm 1982, Voight [29] đưa công thức trung bình cộng số học để tính xấp xỉ tính chất hữu hiệu vật liệu tổ hợp đẳng hướng: n (1.1) C eff   v C  1 Reuss [26] số trường hợp cơng thức trung bình cộng điều hịa cho kết xấp xỉ tốt hơn: C eff  n     v C1    1  1 (1.2) Xuất phát từ nguyên lý biến phân nói chọn trường số, Hill [13] Paul [22] chứng minh tính chất hữu hiệu vật liệu tổ hợp đẳng hướng dù hình học pha nằm giá trị (1.1) (1.2) Cụ thể, đánh giá Hill- Paul viết sau: 1 n  n 1  eff   v C C v C        1   1  (1.3) Bước tiến quan trọng thực Hashin Strickman [11]; [12] Đưa vào trường phân cực (polarization fiels) có giá trị trung bình khác pha khác nhau, ông xây dựng thành công đánh giá tốt đánh giá Hill- Paul: PC  2Cmin   C eff  PC  2Cmax  TIEU LUAN MOI download : skknchat@gmail.com 1  n    v    C* PC  C*     1  C  C*      (1.4) Cmin  C |   1, , n , Cmax  max C |   1, , n Đánh giá Hashin- Strickman cho vật liệu tổ hợp đẳng hướng (bất kể hình hịc pha nào) với tỷ lệ thể tích vα tính chất dẫn nhiệt thành phần Cα cho trước- coi thành tựu bật học vật liệu tổ hợp Câu hỏi đặt liệu có tồn đánh giá tốt đánh giá HashinStrickman có vα Cα cho trước? Hoặc đánh giá tối ưu phải mơ hình hình học cho giá trị hữu hiệu trùng với đánh giá (dưới) Hashin- Strickman đánh giá họ cho hệ số dẫn nhiệt Ceff tối ưu trường hợp vật liệu hai pha cách xây dựng mơ hình cầu lồng Milton [18] đánh giá Hashin- Strickman cho Ceff vật liệu n pha tối ưu cho lớp giới hạn giá trị vα, Cα Bên cạnh đánh giá nói trên, nghiên cứu nhiều nhà khoa học thời gian qua đề cập tới lớp vật liệu cụ thể với thơng tin bổ sung hình học pha Trong trường hợp vật liệu hai pha với pha cốt liệu gồm hạt hình ellipsoid có tỷ lệ thể tích nhỏ phân bố rời rạc xa pha thứ hai liên tục, ta nhận biểu thức tiệm cận hệ số hữu hiệu dựa kết Eshelby Người ta xem xét mơ hình mà tỷ lệ thể tích pha cốt liệu khơng nhỏ tính tới tương tác hai hạt, ba hạt cốt liệu gần để đưa công thức xấp xỉ (xem [7]) Với hình học pha đặc biệt: tuần hồn đơn giản, áp dụng trực tiếp phương pháp số để giải Một phương pháp tiếp cận xấp xỉ hay sử dụng thời gian gần gọi phương pháp tự tương hợp (self-consistent)- áp dụng [6], [14], [15]- nhiên bị phê bình thiếu sở tốn học chặt chẽ, có cho kết sai lệch có ý nghĩa hạn chế Để có đánh giá hẹp đánh giá Hashin- Strickman, người ta tìm cách đưa thêm vào thơng tin bổ sung hình học pha vật liệu thông qua hàm ngẫu nhiên ([24], [16], [21], [28], [30], [31]) Các hàm ngẫu nhiên bậc n (n-point correlation functions) phụ thuộc vào xác suất n điểm lấy tình cờ (với khoảng cách định chúng) rơi vào pha Việc đo đạc hàm ngẫu nhiên thường phức tạp chưa nói đến việc sử dụng TIEU LUAN MOI download : skknchat@gmail.com chúng để tìm đánh giá tối ưu Với phương tiện đo đạc thực nghiệm máy tính đại người ta nhận hàm ngẫu nhiên bậc hai bậc ba cho trường hợp cụ thể (về mặt lý thuyết người ta có hệ số hữu hiệu xác khơng đánh giá có tất hàm ngẫu nhiên tới bậc  ) Miller [20] xem xét lớp vật liệu tổ hợp đẳng hướng- gọi vật liệu cấu trúc đối xứng (symmetric cell material) với rang buộc bổ sung lên hình học pha vật liệu: pha có cấu trúc hình học vi mơ tỷ lệ thể tích khác Hiển nhiên ràng buộc khơng địi hỏi phải cho trước hình học cụ thể vật liệu Miller xây dựng thành công đánh giá cho Ceff vật liệu cấu trúc đối xứng hai pha- nằm đánh giá HashinStrickman: PC  2Cl   C eff  PC  2Cu   1 Cl   v1 C2  v2 C1  ,  v1 C1  v2 C2  1 , (1.5) Cu  max v1 C2  v2 C1 , v1 C1  v2 C2  Các tính chất loại vật liệu nhiều tác giả quan tâm- đáng ý nghiên cứu gần Bruno ([5]; [7]) Các phương pháp nghiên cứu vật liệu tổ hợp đẳng hướng không gian ba chiều dùng để xem xét vật liệu đẳng hướng mặt phẳng kết tương ứng nhận ([9], [10], [13], [28], [29], [30]) Trong luận văn này, dựa nguyên lý lượng cực tiểu (và nguyên lý lượng bù cực tiểu) phương pháp biến phân, đề xuất cách xây dựng đánh giá cho hệ số dẫn nhiệt hiệu có tính tới thơng tin ngẫu nhiên bậc ba hình học pha vật liệu nằm đánh giá Hashin-Strickman áp dụng cho số vật liệu đối xứng phi đối xứng Sử dụng trường thử tổng quát trường Hashin- Strickman, từ nguyên lý lượng cực tiểu nhận đánh cho hệ số dẫn hiệu vật liệu tổ hợp đẳng hướng nhiều thành phần tốt đánh giá xây dựng [1],[4] [7] Đánh giá chứa đựng thơng tin bậc ba hình học pha vật liệu áp dụng cho lớp vật liệu tựa đối xứng mơ hình cầu lồng ba pha Bố cục luận văn: Luận văn bao gồm chương: TIEU LUAN MOI download : skknchat@gmail.com 10 Chương Mở đầu Chương Đánh giá cho hệ số dẫn nhiệt Chương Áp dụng Ngoài chương nêu trên, luận văn gồm phần: - Kết Luận: nêu kết luận văn đạt được, ứng dụng ý nghĩa luận văn đề xuất thêm hướng nghiên cứu thời gian tới - Danh Mục Cơng Trình Khoa Học Đã Cơng Bố - Tài Liệu Tham Khảo TIEU LUAN MOI download : skknchat@gmail.com 11 Chương Đánh giá cho hệ số dẫn nhiệt Trước tiên, ta cần quan tâm đến khái niệm phần tử đặc trưng (RVE) vật liệu tổ hợp Khái niệm: Một phần tử gọi phần tử đặc trưng vật liệu đủ lớn so với cấu trúc vi mơ để coi thực đại diện cho vật liệu xem xét phải đủ nhỏ so với kích thước vĩ mô vật thể chịu lực (và độ dài bước sóng trường hợp tốn động) để tính chất hữu hiệu thực có ý nghĩa Ta xem xét phần tử đặc trưng V vật liệu hỗn độn đẳng hướng vĩ mô cầu có tâm trùng với tâm hệ tọa độ Decac vng góc {x} V cấu tạo từ n thành phần chiếm vùng V  V tích vα (α = 1,…,n; thể tích V coi bẳng 1) hệ số dẫn Cα Như vậy, hệ số dẫn địa phương biểu diễn sau: n C ( x )   C I  ( x ) , x  V (2.1) 1 , x  V  I ( x )   0 , x  V \ V (2.2)  1 Trong đó: Hệ số dẫn C(x) hệ số tỷ lệ quan hệ tuyến tính: J (x)  C (x)E(x) (2.3) Với: J (x) : dòng nhiệt E(x) : Trường gradient nhiệt Dòng J(x) cần thoả mãn phương trình cân bằng:   J ( x)  (2.4) Trong trường E(x) gradient: TIEU LUAN MOI download : skknchat@gmail.com 24   p p v02e2   p  1v03e4   p 1  e4   v v e2  v  p    v v e2  v e4  A  v v e3  v e4   A  v e1  2v e2  v e3  v2 e4  Trong bất đẳng thức trên, ta lấy thay cho p  nhận xét, ta p chọn pα lớn tùy ý Bước tìm quan hệ ei Từ lý thuyết hàm điều hòa nhớ V cầu, ta có:  ( x)   4  xy 1 dy  V x.x  const , x V Như thế: v v   ,iidx     , ij,ijdx   V V     ,ij ,ij dx   V v   A   v e1  e2  v e2  e3  e4  Vì ei khơng phụ thuộc vα nhận xét trên, từ đẳng thức cuối ta nhận được: e2  e1 , e4  e1  e3 Do vậy, (3.3) đưa (        ):  A    A    A    A  v v v e1  e3   v v v  1e1  v e3   v v 1  v e3  v e1  (3.4)  v 1  v 1  v e1  v e3  Đánh giá (33) (36) trở thành: C eU  C eff  C eL (3.5) Trong đó: TIEU LUAN MOI download : skknchat@gmail.com 25  C    , e   C C   C eU C 1 , v 1 , e1  C UA C 1 , v 1 , A C eL với: n n n , v  A  (3.4), e  n  1 n L A n n n n , v 1n , A 1n    e3  Đánh giá chung cho vật liệu tựa đối xứng là: U L C sym  C eff  C sym   (3.6)   U C sym C 1n , v 1n  0max C eU e1  /   L n n C eL C sym C 1 , v 1  0min e1  /  (3.7) Cho minh hoạ số cụ thể: Ta chọn vật liệu ba thành phần với tính chất sau: C1: C2: C3 = 1: 10: 20, v1= 2v3,v2 = 0.1  0.9 Đánh giá (3.6) giá trị e1 = emax e1 = emin (tại giá trị max tương ứng đạt được) cho Bảng 3.1 Đánh giá rộng theo kết trước [23] đưa vào để so sánh Các đánh giá này, với đánh giá Hashin - Strikman [11] so sánh Hình 3.2 U L , Csym Bảng 3.1: Đánh giá hệ số dẫn hiệu vật liệu tựa đối xứng pha Csym - đánh giá (3.6); CPU , CPL - đánh giá từ [23]; e1max ,e1min -tại giá trị max (3.6) đạt v2 L C HS CPL L Csym U Csym CPU C UHS e1max e1min 0.1 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.6 0.7 0.8 0.9 2.505 2.841 3.028 3.231 3.450 3.689 3.950 4.236 4.898 5.719 6.761 8.127 2.921 3.421 3.704 3.897 4.074 4.268 4.482 4.721 5.283 5.995 6.925 8.188 2.921 3.421 3.704 3.956 4.158 4.359 4.578 4.819 5.383 6.092 7.010 8.245 5.600 5.832 5.953 6.097 6.305 6.55 6.818 7.093 7.652 8.223 8.805 9.397 5.676 5.983 6.139 6.298 6.459 6.623 6.818 7.093 7.652 8.223 8.805 9.397 6.223 6.614 6.812 7.012 7.213 7.417 7.622 7.828 8.247 8.674 9.108 9.550 0 0.48 0.62 2/3 2/3 2/3 2/3 2/3 2/3 2/3 2/3 2/3 2/3 0.43 0.22 0.07 0 0 0 TIEU LUAN MOI download : skknchat@gmail.com 26 10 Đánh giá HS Đánh giá P Đánh giá (31) 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 Hình 3.2: Đánh giá hệ số dẫn hiệu vật liệu tựa đối xứng pha với tính chất C1: C2: C3 = 1: 10: 20, v1=2 v3, v2 = 0.1  0.9; Đánh giá (3.6) so sánh với đánh giá P (Pham, 1994) Hashin - Strickman (Hashin – Strickman, 1962) Nhận xét: Các đường đánh giá cận luận văn nằm đường đánh giá Hashin- Strickman đánh giá P chứng tỏ kết nghiên cứu hoàn toàn phù hợp với nghiên cứu trước đánh giá luận văn vừa xây dựng tối ưu đánh giá kinh điển Hashin- Strickman đánh giá P 3.2 Mơ hình cầu lồng ba pha: Trước hết quan tâm đến mơ hình mà A tính xác Mơ hình ta nghĩ đến mơ hình cổ điển cầu lồng Để dễ dàng hiểu mơ hình cầu lồng ba pha, ta cần xem xét mơ hình cầu lồng hai pha xây dựng Hashin [8]: Phần tử đặc trưng V hợp thành pha liên tục VM (thể tích vM) với tính chất CM pha cốt liệu VI (thể tích vI) với tính chất CI Mỗi cầu S I  VI nằm cầu đồng tâm lớn SIM (vỏ cầu SM = SIM\ S I  VI ) với tỷ lệ thể tích cho trước (thể tích S I ): (thể tích SIM) = vI : (vI + vM) V hồn tồn lấp kín cầu lồng với kích thước khác (tới vơ bé) Đây mơ hình lý tưởng không gặp thực tế TIEU LUAN MOI download : skknchat@gmail.com 27 coi mức độ đại diện cho loại vật liệu hai pha dạng ma trận- cốt liệu với pha cốt liệu hạt hình cầu phân bố rời rạc pha ma trận liên tục Hình 3.3: Mơ hình cầu lồng hai pha Xét cầu lồng S IM  V xây dựng hệ tọa độ Decac địa phương xi với tâm trùng với tâm cầu SI SIM: xi  cij x j  const, cij  const, cij ckj   ik (3.8) Ký hiệu: G ( x, y )   1 xy 4 (3.9) Vì V SI, SIM cầu, từ lý thuyết hàm điều hịa ta có (S ký hiệu SI hay SIM, R bán kính S): 1    xi xi  const , x  S  S G ( x, y)dy   R3  , x V \ S  xixi (3.10)  G( x, y)dy  x x  const, x  V i i V   V khác với SIM Cũng từ lý thuyết hàm thế, cầu lồng S IM x  S IM , ta có:  G( x, y)dy  v  G( x, y)dy, I S I S IM  G( x, y)dy  v  G( x, y)dy M  SM (3.11)  S IM TIEU LUAN MOI download : skknchat@gmail.com 28 Với (3.10) (3.11) ta kiểm tra ( x  S IM )  G( x, y)dy  v  G( x, y)dy I VI \ S I V \ S IM    vI   G ( x, y)dy   G ( x, y )dy   V S IM     G( x, y )dy  V M \ S IM (3.12) vI xi xi  xixi   const, vM xi xi  xixi   const Từ (3.8), (3.11), (3.12) A   ij ij dx V def     , ij  ,ij    , ij  ,ij  dx V    ,ij ,ij dx  V v     ta tính được: Đối với x  S I  VI :  I ( x)   G ( x, y )dy VI   G ( x, y )dy  SI   G( x, y)dy VI \ S I v xixi  I  xi xi  xixi   const , 6  M ( x)   G ( x, y )dy VM   G ( x, y )dy   G ( x, y )dy   G ( x, y )dy VM \ S M  S IM SI vM xi xi  xixi   const; Như thế:  ,Iij   ij , ijI  0,  ijM  0, (3.13) TIEU LUAN MOI download : skknchat@gmail.com 29 AIII    ijI  ijI dx  VI I ij    I ij (3.14) dx  0, S I VI S I AIIM  AIMM  0; Đối với x  S M  VM :  I ( x)   G ( x, y )dy VI   G ( x, y ) dy  SI   G ( x, y )dy VI \ S I RI3 v  I  xi xi  xixi  const, xixi  M ( x)   G ( x, y )dy VM  RI3 v  xixi  M  xi xi  xixi   const; xixi Như thế: RI3 ij RI3 xk xlcki clj    ,  32 xm xm 5 3 xm xm  I ij I , ij AMII   ijI ijI dx  VM I ij I ij     dx S M V M S M R   S M V M dx 8 M dr RI    RI  r S M  xm xm 3 S M VM RI (3.15)  8 RI  RI3 / RM3  vI vM , S IM V   AMMM   AMIM   vI vM Trên sở lý thuyết mơ hình cầu lồng hai pha, ta xét đến mơ hình cầu lồng ba pha sau: Quả cầu thuộc pha lồng vỏ cầu thuộc pha 2, sau lồng tiếp với cầu thuộc pha Các cầu lồng với kích thước khác lấp kín khơng gian vật liệu, với thứ tự lồng tỷ lệ thể tích pha cầu lồng (xem Hình 3.4): TIEU LUAN MOI download : skknchat@gmail.com 30 Hình 3.4: Mơ hình cầu lồng ba pha Các tham số hình học A cho mơ hình tính xác [6]:  33  A3  v3 v1  v2 ,   A11  v 2v v  v 1  3   A312  v1v2v3 v1  v2 1   22 2 1  A3  v2 v3 v1  v2   11  A1  A12  A13  A23    A 23   v v , A13   v v 3  3   A222  A211   A212  v1v2 v1  v2 1   (3.16) Cho minh hoạ số cụ thể: Ta chọn cầu lồng pha với tính chất sau: C1: C2: C3 = 50: 10: v1= v3, v2 = 0.25  0.75 Các kết tính tốn cho đánh giá từ (2.33), (2.36), (3.16) hội tụ tới giá trị xác Ceff, đánh giá cũ từ [17], đánh giá Hashin - Strikman [11] cho Hình 3.5 TIEU LUAN MOI download : skknchat@gmail.com 31 Bảng 3.2: Đánh giá hệ số dẫn hiệu cho mơ hình cầu lồng ba pha v2 L CHS L3 C PA C AL C UA C UPA3 C UHS 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 4.534 4.738 4.956 5.188 5.437 5.703 5.989 6.297 6.630 6.991 7.383 4.584 4.796 5.021 5.259 5.513 5.783 6.072 6.381 6.714 7.072 7.460 4.774 5.026 5.288 5.563 5.850 6.150 6.464 6.793 7.136 7.494 7.869 4.774 5.026 5.288 5.563 5.850 6.150 6.464 6.793 7.136 7.494 7.869 5.172 5.506 5.849 6.200 6.557 6.919 7.283 7.646 8.005 8.356 8.695 12.496 13.009 13.527 14.050 14.577 15.109 15.646 16.189 16.736 17.288 17.846 18 17 16 15 14 13 12 11 10 0.25 Đánh giá HS Đánh giá LP Đánh giá xác 0.35 0.45 0.55 0.65 0.75 Hình 3.5: Đánh giá hệ số dẫn hiệu cho mô hình cầu lồng ba pha Giá trị xác từ (2.33), (2.36), (3.16) so sánh với đánh giá LP (Le and Pham, 1991) HS (Hashin – Strickman, 1962) TIEU LUAN MOI download : skknchat@gmail.com 32 Nhận xét: Đánh giá luận văn tốt nhiều so với đánh giá HashinStrickman đánh giá P Ngoài ra, đánh giá cận luận văn hồn tồn trùng với nhau, chứng tỏ với mơ hình này, luận văn xây dựng đánh giá tối ưu Đồng thời luận văn giải vấn đề nêu chương Mở Đầu: mơ hình hình học cụ thể mà đánh giá hữu hiệu trùng với đánh giá cận cận TIEU LUAN MOI download : skknchat@gmail.com 33 KẾT LUẬN Đường lối biến phân sử dụng để đánh giá hệ số dẫn nhiệt hữu hiệu vật liệu tổ hợp đẳng hướng Ưu đường lối biến phân chỗ có sở tốn học chặt chẽ sử dụng thành cơng chí khơng có thơng tin đầy đủ hình học pha vật liệu tổ hợp- thực tế thường gặp phải nghiên cứu vật liệu nhiều pha Trong tỷ lệ thể tich tính chất vật liệu thành phần cho xác hình học pha cụ thể nhiều vật liệu tổ hợp thường gặp lại mang tính tùy tiện, ngẫu nhiên khó mà nắm Thay giải phương trình cân tương thích hình học pha phức tạp- mà riêng việc mơ tả tốn học ba chiều xác hình học pha mục tiêu không thực tế, phương pháp biến phân cho ta đánh giá tính chất hữu hiệu Càng có thêm thơng tin hình học pha ta xây dựng đánh giá hẹp tìm cách sử dụng thông tin nhiều trường hợp đánh giá đủ hẹp giá trị xấp xỉ tính chất hữu hiệu Đó nói mặt ngun tắc, cịn thực vấn đề khó khăn nhiều Đường lối biến phân khơng cho cách thức tốn học cụ thể để dễ dàng áp dụng, chẳng hạn phương pháp phần tử hữu hạn để giải toán kết cấu đàn hồi cần có cách sáng tạo Chúng ta cần phải tìm cực trị phiếm hàm … không với biến chưa biết mà cịn miền khơng xác định V với số thông tin hạn chế V tỷ lệ thể tích… ta thấy dù V khơng xác định ta có đánh giá- đánh giá Hill-Paul đặc biệt đánh giá Hashin-Shtrickman Các hệ số dẫn nhiệt vật liệu tổ hợp xây dựng luân văn phụ thuộc vào hệ số hình học bậc ba tích phân hàm V việc xác định giá trị cụ thể chúng trường hợp chung vấn đề khó Các trường thử sử dụng rộng trường thử trước sử dụng [1] [17], cho đánh giá tốt Kết luận văn hoàn toàn phù hợp với kết kinh điển công bố trước Voight, Russe, Hashin- Strickman, P … TIEU LUAN MOI download : skknchat@gmail.com 34 Ngồi ra, luận văn xem xét số mơ hình hình học hệ số A tính xác hay đơn giản hóa đánh giá tương ứng- nằm đánh giá Hashin-Shtrickman Vấn đề tìm đánh giá tốt cho vật liệu đẳng hướng nhiều pha đánh giá cho lớp vật liệu đối xứng, đặc biệt số pha n ≥ vấn đề khó tiếp tục mục tiêu nghiên cứu TIEU LUAN MOI download : skknchat@gmail.com 35 DANH MỤC CƠNG TRÌNH KHOA HỌC ĐÃ CÔNG BỐ “Đánh giá bậc ba hệ số dẫn vật liệu hỗn độn nhiều thành phần”- Phạm Đức Chính Vương Thị Mỹ Hạnh, Tuyển tập cơng trình Hội Nghị Khoa Học Tồn Quốc Cơ học vật rắn biến dạng lần thứ X, Thái Nguyên, Tháng 11/2010 TIEU LUAN MOI download : skknchat@gmail.com 36 TÀI LIỆU THAM KHẢO Tiếng Việt: Phạm Đức Chính Luận án Phó tiến sĩ khoa học tốn lý Đánh giá tính chất lý vĩ mơ vật liệu đẳng hướng nhiều pha (1995) Phạm Đức Chính Đánh giá tính chất vĩ mơ vật liệu tổ hợp đẳng hướng Hội nghị Cơ học vật rắn biến dạng toàn quốc lần thứ IV, 10/1994 Phạm Đức Chính Vương Thị Mỹ Hạnh Đánh giá bậc ba hệ số dẫn vật liệu hỗn độn nhiều thành phần Tuyển tập cơng trình Hội nghị Khoa học toàn quốc Cơ học Vật rắn biến dạng lần thứ X, 11/2010 Lê Khánh Châu, Phạm Đức Chính Về hệ số dẫn nhiệt vật liệu tổ hợp Hội thảo trao đổi nhiệt chất, Hà Nội (1988) Tiếng Anh: Bruno O.P>, Taylor expansions and bounds for the effective conductivity and the effective elastic moduli of multicomponent composites and polycrystals Asymptotic Analysis 4,339-369 (1991) Budiansky B., On the elastic moduli of some heterogeneous materials J Mech Phys Solids 13, 223 (1965) Hale D.K The physical properties of composite materials J Mater Sci 11, 2105 (1796) Hashin Z., The elastic moduli of heterogeneous materials, J Appl Mech 29,143-150 (1962) Hashin Z., On elastic behaviour of fiber reinfoeced materials of arbitrary transverse phase geometry J Mech Phys Solids 13, 119 (1965) 10 Hashin Z., Rosen W., The elastic moduli of fiber reinforced materials J Appl Mech 31, 223 (1964) 11 Hashin Z., Strikman S A variational approach to the theory of the effective magnetic peomeability of multiphase materials J Appl TIEU LUAN MOI download : skknchat@gmail.com 37 12 Hashin Z., Strikman S A variational approach to the theory of elastic behaviour of multiphase materials J Mech Phys Solids 11, 127-140 (1963) 13 Hill R Theory of mechenical properties of fiber-strengthened materials: I Elastic behaviour J Mech Phys Solids (1964) 14 Hill R., Self-consistent mechanics of composite materials J Mech Phys 13, 213-222 (1965) 15 Kroner E., Berechnung der elastischen Konstanten des Vielkristalss aus den Konstanten eds Einkristalls Z Phys (1958) 16 Kroner E., Elastic moduli of perfectly disorderd composite materals J Mech Phys Solids (1977) 17 Le R C., Pham D C On bounding the effective conductivity of isotropic composite materials ZAMP 42, 614 – 622 (1991) 18 Milton G.W The theory of composite Cambridge University Press, 2001 19 Miller M N., Bounds for effective electrical, thermal an magnetic properties of heterogeneous materials J Math Phys 10, 1988-2004 (1969) 20 Miller M N., Bounds for the effective elastic bulk moduls of heterogeneous materials J Math Phys 10, 2005 (1969) 21 Milton G.W., Kohn R.V., Variational bounds on the effective moduli of anisotropic composites J Mech Phys Solids (1988) 22 Paul B., Prediction of elastic constants of multiphasw materials Trans ASME (1960) 23 Pham D C., Bounds for the effective conductivity and alastic modul of fullydisordenad multicomponent materials Arch Rational Mech Anal 127, 191 198 (1994) 24 Pham D C., Torquato S, Strong-contrast expansions and approximations for the effective conductivity of isotropic multiphase composites J Appl Phys 94, 6591- 6602 (2003) 25 Pham D C., Estimations for the overall properties of same isotropic locallyordened composite Acta Mech 121, 177 – 190 (1997) 26 Reuss A., Berechnung der Fliebgrenze von Mischkristallen auf Grund der Plastizitatzsbedingung fur Einkristalle ZAMM (1929) 27 Torqueto S.Random heterogeneous media New York, Springer, 2002 TIEU LUAN MOI download : skknchat@gmail.com 38 28 Torquato S., Random heterogenous media: Microstructure and improved bounds on effective properties Appl Mech Rev (1991) 29 Voight W., Lehrbuch der Krystallphysik Teuber, Leipzig (1928) 30 Willis J.R., Bounds and self-consistent estimates for the overall moduli of anisotropic composite J Mech Phy Solids (1977) 31 Willis J R., Variational and related methods for the overall properties of composite materials In Advances in Appl Mech Academic press (1981) TIEU LUAN MOI download : skknchat@gmail.com ... ĐẠI HỌC CÔNG NGHỆ VƯƠNG THỊ MỸ HẠNH ĐÁNH GIÁ TRÊN VÀ DƯỚI HỆ SỐ DẪN NHIỆT VẬT LIỆU TỔ HỢP ĐẲNG HƯỚNG Ngành : Cơ học kỹ thuật Chuyên ngành: Cơ học kỹ thuật Mã số: 60 52 02 LUẬN VĂN THẠC SỸ CƠ HỌC... hướng Việc xác định lý tính hữu hiệu vật liệu tổ hợp vấn đề khoa học vật liệu Hệ số dẫn nhiệt vĩ mô vật liệu tổ hợp đẳng hướng phụ thuộc vào tính chất vi mơ vật liệu thành phần cấu tạo nên như: tỷ... để đánh giá hệ số dẫn nhiệt hữu hiệu vật liệu tổ hợp đẳng hướng Ưu đường lối biến phân chỗ có sở tốn học chặt chẽ sử dụng thành cơng chí khơng có thơng tin đầy đủ hình học pha vật liệu tổ hợp-