Tóm tắt luận án Tiến sĩ Cơ học: Xấp xỉ tương đương cơ tính vật liệu tổ hợp có các cốt liệu phức hợp

Chia sẻ: _ _ | Ngày: | Loại File: PDF | Số trang:27

Thêm vào BST

Báo xấu

17
lượt xem 2
download

Download Vui lòng tải xuống để xem tài liệu đầy đủ

Mục tiêu của đề tài là xây dựng các công thức xấp xỉ sử dụng mô hình cốt tương đương để xác định các giá trị hiệu dụng của hệ số dẫn, hệ số đàn hồi vật liệu Composite với cốt phức hợp và áp dụng phương pháp số sử dụng phần tử hữu hạn (PTHH) tính toán cho một số mô hình vật liệu cụ thể.

Chủ đề:

Bình luận(0) Đăng nhập để gửi bình luận!

Lưu

Nội dung Text: Tóm tắt luận án Tiến sĩ Cơ học: Xấp xỉ tương đương cơ tính vật liệu tổ hợp có các cốt liệu phức hợp

BỘ GIÁO DỤC VÀ VIỆN HÀN LÂM KHOA HỌC VÀ ĐÀO TẠO CÔNG NGHỆ VIỆT NAM HỌC VIỆN KHOA HỌC VÀ CÔNG NGHỆ ----------------------------------- NGUYỄN THỊ HƯƠNG GIANG XẤP XỈ TƯƠNG ĐƯƠNG CƠ TÍNH VẬT LIỆU TỔ HỢP CÓ CÁC CỐT LIỆU PHỨC HỢP Chuyên ngành: Cơ học vật rắn Mã số: 944 01 07 TÓM TẮT LUẬN ÁN TIẾN SĨ CƠ HỌC HÀ NỘI 2020
Công trình được hoàn thành tại: Học viện Khoa học và Công nghệ - Viện Hàn lâm Khoa học và Công nghệ Việt Nam. Người hướng dẫn khoa học 1: GS. TSKH. Phạm Đức Chính Người hướng dẫn khoa học 2: PGS. TS. Trần Bảo Việt Phản biện 1: GS. TS. Hoàng Xuân Lượng Phản biện 2: PGS. TS. Nguyễn Trung Kiên Phản biện 3:PGS. TS. Lã Đức Việt Luận án sẽ được bảo vệ trước Hội đồng đánh giá luận án tiến sĩ cấp Học viện, họp tại Học viện Khoa học và Công nghệ - Viện Hàn lâm Khoa học và Công nghệ Việt Nam vào hồi … giờ ..’, ngày … tháng … năm 2020 Có thể tìm hiểu luận án tại: - Thư viện Học viện Khoa học và Công nghệ - Thư viện Quốc gia Việt Nam
1 MỞ ĐẦU 1. Tính cấp thiết của luận án Trong quá trình sản xuất vật liệu Composite do phản ứng hóa học giữa nền – cốt hoặc do kỹ thuật tráng sợi làm hình thành pha trung gian (lớp vỏ bao quanh cốt – interface), mà trong luận án gọi là vật liệu có cốt phức hợp. Cốt phức hợp được hiểu là cốt cùng lớp vỏ bao quanh cốt. Lớp vỏ này có cơ – lý tính không giống pha nền hay cốt và làm ảnh hưởng đến tính chất hiệu dụng vĩ mô của vật liệu. Khi đó, nếu sử dụng phương pháp đường bao để xác định các giá trị hiệu dụng thì các đường bao thường xa nhau, ít có giá trị thực tế. Các phương pháp xấp xỉ trung bình phân tích theo mô hình trụ tròn hoặc hình cầu nhiều lớp tương đối phức tạp khi dùng cho kỹ sư tính toán. Vì vậy, hướng nghiên cứu trong luận án tập trung vào việc xác định các tính chất cơ – lý tính vĩ mô vật liệu tổ hợp có cốt phức hợp, sử dụng phương pháp cốt tương đương để đưa ra các công thức xấp xỉ đơn giản phù hợp với kỹ sư để bước đầu đánh giá tính chất cơ tính của vật liệu sử dụng . Mô phỏng số bằng phương pháp phần tử hữu hạn cũng đươc áp dụng để kiểm nghiệm tính đúng đắn của công thức xấp xỉ. 2. Mục tiêu nghiên cứu của luận án Xây dựng các công thức xấp xỉ sử dụng mô hình cốt tương đương để xác định các giá trị hiệu dụng của hệ số dẫn, hệ số đàn hồi vật liệu Composite với cốt phức hợp và áp dụng phương pháp số sử dụng phần tử hữu hạn (PTHH) tính toán cho một số mô hình vật liệu cụ thể. 3. Các nội dung nghiên cứu chính của luận án Xây dựng công thức xác định hệ số dẫn hiệu dụng của vật liệu Composite cốt sợi phức hợp đồng phương, lớp vỏ bao quanh cốt đẳng hướng hoặc dị hướng (chương 2) Xây dựng được công thức xác định các mô đun đàn hồi hiệu dụng của vật liệu Composite cốt hạt hình cầu phức hợp và Composite cốt sợi phức hợp đồng phương (chương 3) Áp dụng phương pháp PTHH cho bài toán đồng nhất hóa và tính toán số cho một số dạng hình học tuần hoàn nhiều thành phần (chương 4).
2 CHƯƠNG 1. TỔNG QUAN 1.1. Phân loại vật liệu Composite 1.1.1. Theo hình dạng cốt liệu Theo hình dạng cốt liệu, Composite được chia làm hai loại là Composite cốt sợi và Composite cốt hạt. 1.1.2. Theo bản chất vật liệu nền Theo bản chất của vật liệu nền, Composite có 3 dạng: Composite nền hữu cơ, Composite nền khoáng chất, Composite nền kim loại. 1.2. Hệ số dẫn Hệ số dẫn nhiệt c, hệ số dẫn điện c, hệ số khuếch tán D, hệ số thấm nước k, hệ số điện môi (thấm điện) ϵ, hệ số thấm từ (độ từ thẩm) μ. Các hệ số dẫn có công thức toán học dạng giống nhau và đều thỏa mãn phương trình cân bằng giống nhau. Vì thế trong luận án sử dụng hệ số dẫn nhiệt trong các tính toán và ví dụ minh họa làm đặc trưng cho hệ số dẫn. 1.3. Các mô đun đàn hồi Môđun đàn hồi Young E, mô đun đàn hồi trượt µ, mô đun đàn hồi thể tích k, hệ số Poisson ν. 1.4. Phần tử thể tích đặc trưng Phần tử thể tích đặc trưng phải đủ lớn so với các cấu trúc vi mô để đại diện cho các tính chất của vật liệu thành phần và đồng thời phải đủ nhỏ so với kích thước của vật thể để việc xác định các tính chất vĩ mô có ý nghĩa. Hình 1.3 . Phần tử thể tích đặc trưng RVE
3 Phần tử đặc trưng V được cấu thành bởi N thành phần chiếm các không gian V  V có hệ số dẫn nhiệt c , mô đun đàn hồi k ,  (α = 1,…, n). 1.5. Các xấp xỉ và đánh giá xác định các giá trị hiệu dụng của vật liệu 1.5.1. Phương pháp xấp xỉ trung bình 1.5.1.1. Sơ đồ vi phân n  I cI  c Dc cI , c  dc 1  (1.34) dt 1   I t  1 và hệ phương trình:  dk n  dt  1   t  I k I  k Dk k I ,  I , k ,   1  I  1   d  1     D k ,  , k ,   n I I  I I   dt 1   I t  1 n 0  t  1,  I   I (1.36)  1 Với điều kiện ban đầu c0  cM , k 0  k M ,  0   M Dc , Dk , D - tính theo ten xơ Eshelby tương ứng với cốt α. Giá trị hiệu dụng của các mô đun là lời giải của các phương trình tương ứng khi t=1. 1.5.1.2. Phương pháp tự tương hợp n  I c I  c Dc c I ,c   0 (1.38)  1 và hệ phương trình  n   I k I  k Dk k I ,  I , k ,    0  1 (1.39)  n      D k ,  , k ,    0   I I  I I  1 1.5.1.3. Phương pháp Mori – Tanaka
4 Với vật liệu đẳng hướng hai thành phần, kết quả tính các mô đun hiệu dụng theo phương pháp Mori – Tanaka được biểu diễn: 1    p1  p2  p3   cMTA  cM   I  cI  cM  M 1   cI  cM   I  ,    3cM    1  1 k I  kM  1  kMTA  kM   I  k I  kM  M   Pk   I  , .  d kM  d  1   1  I  M  1   MTA  M   I   I  M  M   P   I  .   2 M  2   (1.40) 1.5.1.4. Một số phương pháp xấp xỉ khác - Xấp xỉ tương tác 3 điểm: chứa đựng các thông tin bậc 3 về hình học pha của vật liệu. - Xấp xỉ phân cực: xuất phát từ nguyên lý năng lượng cực tiểu và trường khả dĩ phân cực dạng Hashin-Strickman. 1.5.2. Đường bao của các giá trị hiệu dụng 1.5.2.1. Đường bao Hill – Paul 1.5.2.2. Đường bao Hashin – Strikman Hashin và Strikman đã xây dựng nguyên lý biến phân riêng và đưa vào trường khả dĩ phân cực với các giá trị trung bình khác nhau trên các pha khác nhau và đưa ra đánh giá mới tốt hơn đánh giá của Hill – Paul.  2  d  1   2  d  1  Pk  max   k eff  Pk  min  (1.45)  d   d  P  *max    eff  P  *min  (1.46) Pc   d  1 cmax   ceff  Pc   d  1 cmin  (1.47) 1.5.2.3. Đánh giá bậc 3 của Phạm Đức Chính
5 Cũng xuất phát từ các nguyên lý năng lượng cực tiểu và xây dựng trường khả dĩ phân cực tương tự trường Hashin, Pham đã tìm được đánh giá hẹp hơn đánh giá của Hashin nhờ thành phần nhiễu chứa thông tin bậc ba về hình học pha của vật liệu. 1.5.3. Phương pháp cốt tương đương Khi tính các giá trị hiệu dụng của vật liệu Composite với lớp cốt phủ, trong một số nghiên cứu đã thay thế cốt cùng lớp vỏ bọc bằng một cốt tương đương có cùng kích thước và có các đặc trưng cơ học tương ứng. Có thể kể đến một số công trình như Hashin tính các tính chất đàn nhiệt của vật liệu composite cốt sợi, với lớp cốt bao bọc bởi một lớp mỏng. Qui và Weng tìm mô đun đàn hồi thể tích của vật liệu Composite cốt hạt với cốt được phủ một lớp mỏng. D. C. Pham và B. V. Tran, sử dụng công thức xấp xỉ Maxwell và mô hình cốt tương đương giải quyết bài toán tìm hệ số dẫn của vật liệu Composite cốt hạt với cốt được phủ,.... 1.6. Phương pháp số Một hướng nghiên cứu rất được quan tâm trong lĩnh vực đồng nhất hóa vật liệu đó là phương pháp số mà kỹ thuật số cổ điển đã xây dựng xấp xỉ các trường khả dĩ động học. Phổ biến là phương pháp số sử dụng phần tử hữu hạn (PTHH). 1.7. Kết luận Tìm hiểu lịch sử phát triển của ngành khoa học vật liệu đối với việc xác định các đặc trưng cơ – lý tính vĩ mô giúp tác giả có cái nhìn tổng quan, xuyên suốt với hướng nghiên cứu đặt ra và có kiến thức cơ sở vững chắc. Với việc lựa chọn nghiên cứu theo phương pháp cốt tương đương tác giả mong muốn có được đánh giá dễ hiểu, phù hợp với việc áp dụng rộng rãi cho kĩ sư tính toán và sẽ được trình bày cụ thể trong các chương sau của luận án. CHƯƠNG 2. XẤP XỈ CỐT TƯƠNG ĐƯƠNG HỆ SỐ DẪN VĨ MÔ VẬT LIỆU CÓ CỐT LIỆU PHỨC HỢP 2.1. Vật liệu cốt sợi phức hợp đồng phương, lớp phủ quanh cốt đẳng hướng 2. 1.1. Mô hình vật liệu
6 Hình 2.1. Vật liệu cốt sợi đồng phương 2.1.2. Các công thức đánh giá hệ số dẫn của vật liệu cốt tròn 2.1.2.1. Đường bao Hashin-Shtrikman (HS) 2.1.2.2. Xấp xỉ vi phân (VP) Hệ số dẫn hiệu dụng c eff là nghiệm của phương trình: 1/2  c1  c2  ceff  eff   1 (2.23) c  c2  c1 2.1.2.3. Xấp xỉ tương tác 3 điểm (TT3Đ) 1   2  c eff  Pc 2 (c0 )   1    c0 (2.24)  c0  c1 c0  c2  c0 là lời giải của phương trình: c0  Pc 2  c0  ,  - thông tin hình học bậc 3. 2.1.3. Xấp xỉ cốt tương đương với cốt tròn được phủ (a) (b) Hình 2.4. (a)- Cốt tròn hai lớp đặt trong pha nền vô tận (b)- Cốt tương đương đặt trong pha nền vô tận
7 Gọi T   ;T   ;T   là nhiệt độ trên các pha (1), (2), (3): 1 2 3  B   B  T     A1r  1  cos ; T     A2 r  2  cos ; T 3  A3r cos 1 2  r   r  Các hằng số A, B được xác định nhờ vào các điều kiện: - Điều kiện liên tục về nhiệt độ giữa các pha: T    r2 ,   T    r2 ,  ; T    r3 ,   T    r3 ,  (2.38) 1 2 2 3 - Điều kiện liên tục về dòng nhiệt giữa các pha: q   n  q    n ; q   n  q    n 1 2 2 3 (2.39) - Và nhiệt độ trên biên r1: T    r1 ,    .r1 cos . 1 (2.40) Với β – gradient nhiệt độ áp lên mẫu, β
8 2 3 với 2'  ; 3'  . (2.65) 2  3 2  3 Với hệ số dẫn của cốt tương đương tìm được theo công thức (2.64), ta thay trở lại (2.62), sẽ tìm được hệ số dẫn hiệu dụng của vật liệu. Có thể viết dưới dạng: 1 23  23   c eff  c1    1   c1 (2.66)   1  c23  c1 2c1  1 c23  c1 2c1 Nếu áp dụng vào công thức (2.23) của phương pháp xấp xỉ vi phân thì khi đó hệ số dẫn hiệu dụng là kết quả của phương trình: 1/2  c1  c23  c eff  eff  c  c  1 (2.67) c  23 1 Khi bổ sung thêm thông tin hình học bậc ba cho cấu trúc vi mô của mô hình nền – cốt tương đương, dựa theo công thức (2.24), (2.25) có thể đề xuất xấp xỉ tương tác 3 điểm khi tính hệ số dẫn hiệu dụng: 1     3  ceff   1  2   c0 (2.68)  c1  c0 c0  c23  trong đó c0 là nghiệm của phương trình: 1   2  c0   1    c0 (2.69)  c1  c0 c0  c23  Tổng quát, các cốt tròn được phủ làm từ các vật liệu khác nhau hoặc có thể cùng một loại vật liệu nhưng khác nhau tỉ lệ thể tích (hình 2.1 c).  1  n  23   c eff    1   c1 (2.70)   1 c  c  2c1   1 23 Lấy ví dụ c1  1, c2  5, c3  20 (W .m1.K 1 );2  3 và
9 c1  20, c2  5, c3  1(W .m1.K 1 );2  3 . Hệ số dẫn hiệu dụng thay đổi theo tỉ lệ thể tích giữa các pha, kết quả thể hiện trên hình 2.5, 2.6. Hình 2.5. Hình 2.6 c1  1, c2  5, c3  20 c1  20, c2  5, c3  1 Với HSU, HSL - đường bao trên và dưới của Hashin- Strikman, VP-CTĐ : theo công thức (2.67), ĐTLN: công thức (2.66) , TT-CTĐ: công thức (2.68) với phân bố ngẫu nhiên . 2.1.4. Kết quả thực nghiệm Hình 2.8. Hệ số dẫn của composite sợi abaca Liu và cộng sự đã đo hệ số dẫn ngang của sợi abaca – nền epoxy bằng kỹ thuật khuếch tán flash (TN) 2.2. Vật liệu cốt sợi dọc trục với cốt phức hợp, lớp phủ quanh cốt dị hướng 2.2.1. Mô hình vật liệu
10 cT cT c3 c3 cn cn c1 c1 (a) (b) (a) (a) (b) Hình 2.9. Vật liệu cốt sợi dọc trục với cốt phức hợp, lớp phủ quanh cốt dị hướng Giả thiết lớp phủ quanh cốt mỏng, có tỉ lệ thể tích Δ  1 .Chia lớp đó thành 2m lớp vô cùng mỏng với tỉ lệ thể tích Δ / 2m , hệ số dẫn luân phiên ở các lớp là c12 , c22 . Hệ số dẫn hiệu dụng của hỗn hợp cốt được phủ (cốt tương đương) c23 :  Δ  c2  c3 c3  c2 1  1   c22  c3 c3  c22     Ο c23  c3  2m  m 2c12  m 2c22   Δ  2   Δ  c23  c3  2.cN   cI 2   c3    Ο Δ 2 2   (2.74) với: 1 1 2 cT  2   2c1 c 2 c2  c2 , cN  1 2 22 , cI  cT cN  c12c22 . (2.75) c2  c2 Khi lớp phủ có tỉ lệ thể tích đáng kể, tương tự phương pháp xấp xỉ vi phân, từ (2.74) ta xây dựng được: 1  cI    c  2 2 dc  (2.77) d 1   2cN với điều kiện biên: c  0  c3 ; c23  c 2  (2.78) Trường hợp cN = const, cT =const từ phương trình (2.77): cI    cI  c3  cI   cI  c3  cI  3' cN 3 c23  , 3'  (2.80) cI 3  2  cI  c3    c3  c    I ' c 3 N
11 Hệ số dẫn hữu hiệu c eff của vật liệu đồng nhất 3 pha nền, cốt và lớp phủ dị hướng. 1     c eff   23  1   c1 (2.81)  c23  c1 2c1  2.3. Kết luận Chương 2 đã xây dựng được các công thức xấp xỉ cho việc xác định giá trị hiệu dụng của hệ số dẫn ngang vật liệu cốt sợi phức hợp đồng phương theo cách tiếp cận cốt tương đương. Kết quả nghiên cứu trong chương này đã được công bố trong các công trình khoa học [2, 4, 7, 8]. CHƯƠNG 3. XẤP XỈ CỐT TƯƠNG ĐƯƠNG HỆ SỐ ĐÀN HỒI VĨ MÔ VẬT LIỆU CÓ CỐT LIỆU PHỨC HỢP 3.1. Mô đun đàn hồi của vật liệu composite cốt hạt hình cầu với cốt phức hợp 3.1.1. Mô hình vật liệu Hình 3.1. (a) - Mô hình vật liệu composite với cốt hình cầu được phủ; (b) - Mô hình cốt tương đương 3.1.2. Mô đun đàn hồi thể tích Áp lên mẫu vật liệu hình 3.1a một ten xơ biến dạng vĩ mô E. Phương trình chuyển vị trên các pha: B B B urI 2  AI 2 r  I22 , urI 1  AI 1r  I21 , uM  AM r  M2 (3.4) r r r Ứng suất trên các pha:
12 BI 2 B  rrI 2  3k I 2 AI 2  4 I 2 3 ,  rrI 1  3k I 1 AI 1  4 I 1 I31 , r r (3.5) BM  rrM  3kM AM  4M 3 . r Với điều kiện biên: Khi r  0 thì ur  0  BI 2  0 ; 𝑟 → ∞ thì  r  E0  AM  E0 . Với E0 là hằng số thỏa mãn E = E01. Các hệ số A, B còn lại được tìm dựa vào điều kiện liên tục về chuyển vị và biến dạng giữa các pha khi : urI 2  urI 1 , rI 2   rI 1 ; urI 1  urIM , rI 1   rM . Tương tự đối với mô hình 3. 1b, chuyển vị trên các pha: B B urEI  AEI r  EI2 , uM  AM r  M2 . (3.13) r r Với BEI  0, AM  E0 còn các hệ số khác được tính theo điều kiện liên tục về chuyển vị và biến dạng: r  REI  RI 1  urEI  urM , rEI   rM . (3.14) Chuyển vị trên biên RI1 ở cả 2 mô hình là như nhau, từ đó ta có: B urI 1  RI 1   urEI  RI 1   AI 1RI 1  I21  AEI RI 1 (3.17) RI 1 4akI 2  I 1  4akI 1 I 1  3ak I 2 k I 1  4ak I 1 I 1 R3 kEI   , a  I32 3akI 2  3akI 1  3kI 2  4I 1 RI 1 (3.19) kM  k*M 4 k  kM  EI  kEI  kM  eff , k*M  M (3.20) kEI  k*M 3 Ngoài ra, dựa theo mô hình quả cầu lồng nhau 2 pha của Hashin, tác giả và cộng sự đề xuất công thức tính mô đun đàn hồi thể tích của cốt tương đương với lớp phủ giống pha nền, như sau: 1   I' 2 I' 1  4 kEI      k*I 1 , k*I 1   I 1 (3.21)   I 2 *I 1 k k k I1  k*I 1  3 k EI nhận được theo hai cách tiếp cận khác nhau theo công thức (3.19) và (3.21) có kết quả trùng khớp nhau. 3.1.2. Mô đun đàn hồi trượt
13 Mô đun đàn hồi trượt của mô hình cốt tương đương (hình 3.1b) theo kết quả phân bố thưa của Eshelby có dạng:     eff  M 1   EI AEId  EI  1  (3.27)   M  Với mô hình cốt được phủ (hình 3.1a), tác giả đề xuất một dạng tương tự cho mô đun đàn hồi trượt giống công thức (3.27):       eff  M 1  I 1 AId1  I 1  1  I 2 AId2  I 2  1  (3.30)   M   M  Đồng nhất (3.27) và (3.30), suy ra:  EI   M M  *M  C d *M  (3.32) M  *M  C d M 9 K M  8 M với *M  M 6 K M  12M 1  d  I1    Cd   I 1 AI 1   1  I 2 AId2  I 2  1  (3.33) EI   M   M  Để xác định mô đun trượt hiệu dụng của cốt tương đương theo công thức (3.32), cần xác định các thành phần của biến dạng d d lệch AI 1 , AI 2 . Xem xét vật liệu ở trạng thái trượt thuần túy. Các thành phần chuyển vị trong hệ tọa độ cầu có dạng : ur  U r  r  sin 2  cos 2 , u  U  r  sin  cos  cos 2 , (3.34) u  U  r  sin  sin 2. 6 3c 5  4 d với U r  ar  br 3  4  . , (3.41) 1  2 r 1  2 r 2 7  4 3 2c 2d U  ar  br  4  2 , U  U . (3.42-3.43) 1  2 r r Với các điều kiện biên: - Khi r = 0 : U r  0  cI 2  d I 2  0 .
14 -Khi r → ∞ : U r  E0 r  aM  E0 , bM  0 . Khi r  RI 1 , r  RI 2 : từ điều kiện liên tục về ứng suất  rr , r , r  và chuyển vị ur , u , u  giữa các pha, ta có các phương trình để xác định 8 hệ số còn lại. Véc tơ chuyển vị ở mỗi pha: u  z   ur er  u e  u e Trung bình biến dạng ở pha I2 được xác định qua biểu thức: ε I 2   u z   nz dS . r  RI 2 mà ε I 2  A : E  AId2 E0  e1  e1  e2  e2  ,  21 4  4  5 I 2  d I 2  1  AId2   aI 2  bI 2 RI22   . (3.67)   5 1  2 I 2  5 1  2 I 2  RI32  E0 Tương tự, trung bình biến dạng trên pha I1+I2:  21 4  4  5 I 1  d I 1  ε I 12   aI 1  bI 1RI21    e1  e1  e2  e2  ,   5 1  2 I 1  5 1  2 I 1  RI31   21 4  4  5 I 1  d I 1  1 AId12   aI 1  bI 1RI21   (3.69)   5 1  2 I 1  5 1  2 I 1  RI31  E0 I 2 d Mà AId1  AId2  A  AId2 I 1 I 12   (3.71) Biểu diễn của mô đun thể tích kEI trong (3.21) gợi ý cho chúng tôi một công thức tương tự xấp xỉ Maxwell cho mô đun đàn hồi trượt của cốt tương đương: 1  I' 1 I' 2  9 K I 1  8 I 1 SEI      *I 1 ; *I 1  I1   I 1  *I 1  I 2  *I 1  6 K I 1  12I 1 (3.72) Với giá trị kEI , EI được đưa ra trong công thức (3.21) và (3.72), sẽ được sử dụng như các mô đun đàn hồi của cốt tương đương dạng đơn giản. 3.1.4. Công thức tổng quát
15 1  n  EI M  4 k eff       k*M ; k*M  M (3.73)   1 kEI  k*M kM  k*M  3 1  n EI M  9k M  8 M  eff      *M ; *M  M .   1  EI  *M M  *M  6kM  12M (3.74) 3.1.5. Kiểm tra và so sánh M  1 GPa, I 1  5 GPa, I 2  25 GPa, I 1  I 2  M  0.3, I 1  2I 2 . Hình 3.2. Mô đun đàn hồi thể Hình 3.3. Mô đun đàn hồi tích hiệu dụng trượt hiệu dụng 3.2. Mô đun đàn hồi của vật liệu composite cốt sợi phức hợp đồng phương 3.2.1. Mô hình vật liệu Hình 3.6. Mô hình vật liệu cốt sợi 3 pha đồng phương 3.2.2. Mô đun đàn hồi diện tích hiệu dụng   I' 1 I' 2  K EI      I1 , (3.79)  K I 1  I 1 K I 2  I 2  trong đó KI1, KI2 là các mô đun đàn hồi diện tích của pha I1, I2. Mô đun đàn hồi diện tích hiệu dụng:
16 1   EI M      M . eff K 23 (3.83)  K EI  M K M  M  3.2.3. Mô đun đàn hồi trượt dọc hiệu dụng 1     eff   12 EI  M   M , (3.95)   EI  M 2M  12  I' 2 với 12 EI   I 1  . (3.96) 1 '  I1  I 2   I 1 2 I 1 3.2.4. Mô đun đàn hồi Young dọc trục và hệ số Poisson 4 I' 2 I' 1  I 2  I 1  2 EEI   ' I 2 EI 2  ' I 1E I 1  ,  2 I' 1 1  I 2  2 I22   2 1  ' I2 I1  2 I21   2 1  I1  EI 2 EI 1 EI 1  EI        ' '   I' 2I' 1  I 2  I 1   EI 2 1  I 1  2 I21  EI 1 1  I 2  2 I22  .    I2 I2 ' I1 I1 2  '  EI 1I 1 1  I 2  2 I 2  EI 2I 2 1  I 1  2 I 1  EI 2 1   I 1  2   4 IM  EI  M  2 E1eff   I EEI  M EM  ,  2M 1  EI  2 EI 2   2 1  I M  2 M2   2 1  M EEI EM EM 12eff  EI EI  M M   EIM  EI  M   EEI 1  M  2 M2  EM 1  EI  2 EI2   .      EM M 1  EI  2 EI  EEIEI 1  M  2 M  EEI 1  M  2 2   3.2.5. Mô đun đàn hồi trượt ngang 1  EI M  K M * M 23 eff     *M ; *M  ,   EI  *M  M  *M  K M  2 M
17 1  I' 1 I' 2  K I1 * I1  EI      *I 1; *I 1  .   I 2  *I 1 M  *I 1  K I 1  2 I 1 3.3. Kết luận Chương 3 đã xây dựng được các công thức xấp xỉ xác định giá trị hiệu dụng mô đun đàn hồi của vật liệu Composite có cốt phức hợp với cốt hình cầu hoặc trụ tròn theo phương pháp tiếp cận cốt tương đương. Kết quả nghiên cứu trong chương này đã được tác giả công bố trong các công trình khoa học [1, 3, 5, 6]. CHƯƠNG 4. MÔ PHỎNG SỐ PHẦN TỬ HỮU HẠN VẬT LIỆU PHỨC HỢP 4.1. Vật liệu tuần hoàn (a) (b) (c) (c) (d) Hình 4.2. Một số nhân tử tuần hoàn, a - lập phương đơn giản, b - lập phương tâm khối, c - lập phương tâm mặt, d – hình vuông , e – hình lục giác. 4.2. Các công thức tính xuất phát 4.2.1. Mô đun đàn hồi Trường ứng suất σ(x), chuyển vị u(x) vi mô có đặc điểm: σ  x   σ  x  n1a1e1  n2 a2e2  n3a3e3  ; n1 , n2 , n3  N . (4.3) u  x  = E.x + u'  x  (4.4)
18 trong đó u’(x) là trường chuyển vị rối loạn (nhiễu), điều kiện tuần hoàn là: u'  x   u'  x  n1a1e1  n2 a2e2  n3a3e3  . (4.5) Các phương trình thỏa mãn trên nhân tuần hoàn   σ  x   0 ; ε  x   u  uT  ; σ  x   C  x  :  E + ε'  x   1 2 σ(x)  n : phản tuần hoàn U 4.2.2. Hệ số dẫn Trường dòng nhiệt vi mô tuần hoàn và có dạng: q  x   q  x  n1a1e1  n2 a2e2  (4.13) Trên nhân tử tuần hoàn phải thỏa mãn các phương trình:   q(x)  0 ; q(x) = -c  x  T  x  ; T  x   T 0  T .x  T *  x  trong đó T0 là nhiệt độ bên trong vật liệu, T *  x  là trường nhiệt nhiễu loạn do các cốt sinh ra và thỏa mãn điều kiện tuần hoàn: T *  x   T *  x  n1a1e1  n2 a2e2  (4.18) q(x)  n : phản tuần hoàn U (4.19) 4.3. Phần mềm Cast3M 4.4. Tính toán cho mô hình vật liệu cụ thể và so sánh kết quả 4.4.1. Vật liệu Composite cốt sợi phức hợp đồng phương 4.4.1.1. Hệ số dẫn ngang khi các pha đồng nhất, đẳng hướng Hình 4.4. (a) - mảng lục giác, (b) - mảng hình vuông