Xem mẫu
Hội nghị Cơ học toàn quốc lần thứ X
Hà Nội, 8-9/12/2017
Mô phỏng số FFT hệ số dẫn vĩ mô vật liệu hai pha dạng
nền-cốt liệu elliptic và các phương pháp xấp xỉ
Nguyễn Văn Luật1, Nguyễn Trung Kiên2
Đại học Công nghiệp Hà Nội
Đại học Giao thông vận tải Hà Nội
1
2
Email: luatnv1980@gmail.com
Tóm tắt
Bài báo trình bày phương pháp biến đổi Fourier (FFT) để tính hệ số dẫn vĩ mô cho vật liệu hai pha dạng nềncốt liệu trong không gian hai chiều, trong đó pha cốt liệu có hình dạng elliptic. Xác định tính chất dẫn vĩ mô (tính chất
dẫn hiệu quả) của vật liệu bằng phương pháp FFT đối với một số mô hình tuần hoàn trong không gian hai chiều, có tỉ
lệ thể tích giữa các pha thay đổi và so sánh với các phương pháp xấp xỉ khác.
Từ khóa: hệ số dẫn, cốt liệu dạng elliptic, phương pháp biến đổi Fourier
1. Mở đầu
Các loại vật liệu tổ hợp ngày nay được áp dụng trong hầu hết các lĩnh vực khoa học, kỹ thuật và đời
sống. Việc nghiên cứu tính chất dẫn vĩ mô hay đồng nhất hóa vật liệu được nhiều nhà khoa học quan tâm
nghiên cứu và đã đưa ra nhiều kết quả xấp xỉ cho các mô hình vật liệu khác nhau. Đối với các mô hình vật
liệu trong tính toán để cho đơn giản có thể được lý tưởng hóa hình học dưới dạng cốt liệu hình cầu hoặc
trong không gian hai chiều là hình tròn. Tuy nhiên trong thực tế cốt liệu có những hình dạng phức tạp hơn
nhiều và cần xấp xỉ dưới dạng hình học không tròn mà có dạng như elliptic. Tính chất vĩ mô của vật liệu tổ
hợp phụ thuộc vào nhiều yếu tố phức tạp như cấu trúc hình học pha, các tính chất của vật liệu thành phần,
tỷ lệ thể tích giữa các pha. Do đó trong các nghiên cứu chủ yếu chỉ tìm được cận trên, dưới và các công thức
xấp xỉ áp dụng cho một số mô hình vật liệu. Hướng tiếp cận để tính xấp xỉ cho các mô hình như của
(Maxwel,1884), (Winner, 1912), (Voight, 1928), (Reuss, 1929), (Bruggeman, 1935), (Hamilton and crosser,
1962), (Lewis and Nielsen, 1970), (Mori and Tanaka, 1973). Một hướng tiếp cận khác là xây dựng biên trên
và biên dưới cho hệ số dẫn vĩ mô như (Hill, 1952), (Hashin and Strikman, 1962), (Pham DC, 1996)…Ngoài
ra các phương pháp số hiện nay cũng là cách tiếp cận hiệu quả trong việc xác định tính chất vĩ mô của vật
liệu như phương pháp phần tử hữu hạn (FEM), phương pháp biến đổi Fourier (FFT). Phương pháp FFT áp
dụng trong lĩnh vực cơ học vật liệu tính mô đun đàn hồi cho vật liệu tổ hợp được đề xuất đầu tiên vào năm
bởi (Moulicec and Subquet, 1994). Trong bài báo này sử dụng phương pháp FFT để tính hệ số dẫn vĩ mô
cho một số mô hình vật liệu hai pha với pha cốt liệu có hình học dạng elliptic được sắp xếp tuần hoàn trong
pha nền, trong đó có so sánh với các phương pháp xấp xỉ khác.
2. Phương pháp biến đổi Fourier (FFT)
Nội dung cơ bản của phương pháp biến đổi Fourier là thiết lập được phương trình LippmanSchwinger đối với bài toán không đồng nhất và sử dụng toán tử Green tuần hoàn. Sau đó sử dụng thuật toán
lặp. Ứng xử của các vật liệu thành phần được mô tả bởi định luật Fourier:
J (x) = -C(x).E(x)
(1)
Nguyễn Văn Luật, Nguyễn Trung Kiên
trong đó E(x) và J (x) lần lượt là trường gradient nhiệt độ và dòng nhiệt địa phương thỏa mãn phương
trình cân bằng
E(x) T (x), . J ( x) 0
Trường gradient E(x) và T (x) có thể tách thành các thành phần sau:
(2)
E(x) E0 e per
T (x) E0·x T per
(3)
per
trong đó E0 là gradient vĩ mô đồng nhất đối với phần tử đặc trưng, e gọi là thành phần nhiễu có tính
chất tuần hoàn. Do tính chất tuần hoàn nên ta có:
e per (x)V 0; E0 E(x)V
với ký hiệu •V là trung bình trên trên thể tích của phần tử đặc trưng V, •V
(4)
1
•dx.
V V
Bài toán trên phần tử đặc trưng có thể quy về tìm các thành phần T per , e per . Đưa vào môi trường làm chuẩn
có hệ số dẫn C0 , phương trình cân bằng trở thành
·J ·(C 0 C )E(x) 0
(5)
với C (x) C (x) C 0
Thay E(x) từ (2) vào (5) viết lại dưới dạng tương đương sau
·C 0T per ·τ (x) 0
trong đó tenxơ τ( x) gọi là tenxơ "cực" được xác định bởi : τ (x) C (x) E0 e per (x)
(6)
, e per và τ (x) được biểu diễn dưới dạng chuỗi Fourier:
1
(7)
F(x) Fˆ ( )ei .x , Fˆ ( ) F(x)ei .x F(x)ei .x dx
VV
per
per
ˆ per , eˆ per và ˆ .
trong đó F chỉ T , e và τ( x) , còn Fˆ là biến đổi Fourier của các đại lượng này, đó là T
Do tính chất chu kỳ của phần tử đặc trưng nên T
per
Ở đây tính chất tuần hoàn được thể hiện bởi
d
ei .( x h ) ei .x , h 2m j a j e j , ( j 1, 2,..., d )
(8)
j 1
với m j là số nguyên bất kỳ, 2a j là kích thước của phần tử đặc trưng song song với trục x j , e j là vectơ cở
sở theo hướng x j , d là số chiều của không gian.
k ek , j
n j
aj
Thay các biểu diễn dạng chuỗi Fourier của T
(n j 0, 1, 2...) , không tổng theo j
per
, e per và τ(x) vào phương trình (6) thu được
0
( mCmj
j )Tˆ per ( )ei .x i mˆ( )ei .x 0
per
per
từ đó các trường Tˆ và eˆ có thể xác định như sau:
(9)
Mô phỏng số FFT hệ số dẫn vĩ mô vật liệu hai pha dạng nền- cốt liệu elliptic và các phương pháp xấp xỉ
i .ˆ( ) per
.ˆ( )
Tˆ per
, eˆ ( ) i Tˆ per ( )
Γˆ 0 ( ).ˆ( )
0
0
.C
.C
0
trong đó Γ ( ) là toán tử Green phụ thuộc môi trường đồng nhất C0 được xác định bởi
Γ 0 ( )
.C 0
(10)
(11)
Từ đó thu được phương trình Lippman-Schwinger
ˆ ( ) E
ˆ 0 ( ) Γ
ˆ 0 ( ).ˆ( ) E
ˆ 0 ( ) Γ
ˆ 0 ( ). (C ( ) C 0 )* Eˆ ( )
E
(12)
Nghiệm của phương trình được tìm bởi sơ đồ lặp sau:
ˆ i 1 ( ) Γˆ 0 ( )·(C ( ) C 0 ) E
ˆ i ( ), 0
E
(13)
i 1
ˆ E0 , 0
E
Chú ý rằng Γ0·C 0 Ei ( ) Ei ( ) với 0 xem Michel(1999-[5]), phương trình (13) được viết lại
dưới dạng sau:
ˆ i 1 ( ) E
ˆ i ( ) Γˆ 0 ( ).Jˆ i ( ), 0
E
(14)
i 1
ˆ E0 , 0
E
trong đó Jˆ i ( ) là biến đổi Fourier của J i ( x) . Liên hệ giữa trường dòng J và trường gradient E trong
không gian Fourier được biểu diễn bằng biểu thức:
ˆ ( )
Jˆ ( ) C ( )* E
(15)
trong đó ký hiệu "*" là tích "convolution". Biến đổi Fourier của tenxơ hệ số dẫn:
C ( ) C (x)ei .x dx C I ( )
(16)
V
với C , I lần lượt là tenxơ hệ số dẫn và hàm dạng của pha , I ( ) được xác định theo Nemat-Nasser
(1999-[7]):
I ( )
1
ei .x dV
V V
(17)
Thay các biểu thức (15), (16) vào (14) thu được
ˆ i 1 ( ) E
ˆ i ( ) Γˆ 0 ( ). C I ( ) * E
ˆ i ( ), 0
E
i
1
0
ˆ E , 0
E
(18)
Để xác định hệ số dẫn vĩ mô của vật liệu composite, cho phần tử đặc trưng chịu tác dụng của gradient vĩ mô
E0 . Khi quá trình lặp theo (18) hội tụ (số hạng đầu tiên E1 E0 ), ta có
J( 0) C eff E0
eff
(19)
trong đó C là hệ số dẫn hiệu quả của vật liệu composite. Từ đó rút ra thuật toán số để xác định hệ số dẫn
của vật liệu nhiều thành phần có cấu trúc tuần hoàn:
Nguyễn Văn Luật, Nguyễn Trung Kiên
ˆ 1 ( ) 0 0; E
ˆ 1 (0) E0
Bước i=1: E
ˆ 1 ( )
Jˆ 1 ( ) C ( )* E
ˆ i ( ) và Jˆ i ( ) đã biết
Bước i: E
Kiểm tra hội tụ
i
ˆ i 1 ( ) E
ˆ i ( ) Γˆ 0 ( ).J ( )
E
ˆ i 1 ( )
Jˆ i 1 ( ) C ( )* E
Kiểm tra điều kiện hội tụ được xác định bằng biểu thức sau:
‖ Jˆ i 1 ( ) Jˆ i ( ) ‖
, với
‖ Jˆ i ( ) ‖
3
là sai số cho trước ( 10 )
3. Một số phương pháp xấp xỉ
Trong mục này giới thiệu một số phương pháp tính xấp xỉ, đánh giá hệ số dẫn vĩ mô (Ceff) của vật liệu
nền-cốt liệu dạng elliptic đẳng hướng trong không gian hai chiều với các ký hiệu: CI, vI là hệ số dẫn và tỉ lệ
thể tích của pha cốt liệu, CM, vM là hệ số dẫn và tỉ lệ thể tích của pha nền.
3.1.1 Xấp xỉ Hamilton (1962-[1])
C eff CM .
với n
3
CI (n 1)CM (n 1)vI (CM CI )
CI (n 1)CM vI (CM CI )
(20)
, là thông số hình học của cốt liệu, trong trường hợp cốt liệu dạng elliptic 0.4 .
3.1.2 Xấp xỉ Maxwell (1892-[1])
C eff (
vI
vM
) 1 (d 1)CM
CI (d 1)CM CM (d 1)CM
(21)
3.1.3 Xấp xỉ Lewis-Nielsen (1970-[1])
C eff
CI
1
1 m
CM
1 A vI
,
, 1
CM .
vI
2
CI
1 vI
m
A
CM
(22)
Trong đó A, m là các hệ số hình học cốt liệu (lấy giá trị A=1.5, m 1 với hình cầu, A=3, m 0.637
với hình dạng khác).
3.1.4 Xấp xỉ Mori-Tanaka (1973-[4])
Trong không gian hai chiều, cốt liệu elliptic xấp xỉ Mori-Tanaka có dạng:
C eff CM vI (CI CM ).
(23)
Mô phỏng số FFT hệ số dẫn vĩ mô vật liệu hai pha dạng nền- cốt liệu elliptic và các phương pháp xấp xỉ
trong đó
CM (CI CM )(1 r )2
, r là tỉ số giữa hai bán trục chính của elliptic.
2(CI rCM )(CM rCI )
3.1.5 Đánh giá Hashin-Strikman (1962-[3])
Hashin-Strikman dựa trên nguyên lý biến phân riêng đưa vào trường khả dĩ phân cực đã xây dựng
được đánh giá trên và dưới cho hệ số dẫn vĩ mô của vật liệu nhiều thành phần đẳng hướng.
PC ((d 1)Cmin ) C eff PC ((d 1)Cmax ),
(24)
1
v
PC (C0 )
C* ,
C C*
C* (d 1)C0 , Cmin min C1 , , Cn , Cmax max C1 , , Cn
4. Kết quả so sánh
Trong mục này sẽ đưa ra kết quả tính toán FFT hệ số dẫn vĩ mô cho một số mô hình vật liệu có cốt
liệu dạng elliptic trong không gian hai chiều và so sánh với các phương pháp xấp xỉ khác. Hàm dạng (17)
cho cốt liệu elliptic trong không gian hai chiều [2]
(25)
I ( ) 2S J1 ( ) / .ei .xc ( )
Trong đó J1 là hàm Bessel loại 1, S là diện tích bề mặt của cốt liệu, xc ( ) là véc tơ xác định vị trí trọng
tâm của pha cốt liệu
( 12 a12 22 a22 )1/2
a1 , a2 là độ dài các bán trục chính của cốt liệu elliptic, 1 , 2 là các thành phần của theo các trục của
elliptic.
- Kết quả tính toán FFT cho mô hình vật liệu có cốt liệu elliptic sắp xếp tuần hoàn trong phần tử đặc
trưng (unit cell) theo hình vuông (square) và so sánh với các mô hình xấp xỉ khác. Do cốt liệu không chồng
lấn nên tỉ lệ thể tích cốt liệu chỉ có thể tăng đến một giới hạn nhất định. Trên hình 2 với pha nền có hệ số
dẫn CM=1, pha cốt liệu có hệ số dẫn CI=10, a2 0.2a1 , trong trường hợp này cả FFT và các xấp xỉ khác
đều nằm trong đánh giá của Hashin-Strickman (HS) và gần nhau khi tỉ lệ thể tích cốt liệu nhỏ. Kết quả trên
hình 3 với CM=10, CI=2 thì có thể thấy chỉ có FFT nằm trong đánh giá của HS, điều này khẳng định độ tin
cậy của phương pháp FFT so với các phương pháp xấp xỉ trước đó.
Hình 1: Mô hình cốt liệu elliptic có cấu trúc tuần hoàn
nguon tai.lieu . vn
Hà Nội, 8-9/12/2017
Mô phỏng số FFT hệ số dẫn vĩ mô vật liệu hai pha dạng
nền-cốt liệu elliptic và các phương pháp xấp xỉ
Nguyễn Văn Luật1, Nguyễn Trung Kiên2
Đại học Công nghiệp Hà Nội
Đại học Giao thông vận tải Hà Nội
1
2
Email: luatnv1980@gmail.com
Tóm tắt
Bài báo trình bày phương pháp biến đổi Fourier (FFT) để tính hệ số dẫn vĩ mô cho vật liệu hai pha dạng nềncốt liệu trong không gian hai chiều, trong đó pha cốt liệu có hình dạng elliptic. Xác định tính chất dẫn vĩ mô (tính chất
dẫn hiệu quả) của vật liệu bằng phương pháp FFT đối với một số mô hình tuần hoàn trong không gian hai chiều, có tỉ
lệ thể tích giữa các pha thay đổi và so sánh với các phương pháp xấp xỉ khác.
Từ khóa: hệ số dẫn, cốt liệu dạng elliptic, phương pháp biến đổi Fourier
1. Mở đầu
Các loại vật liệu tổ hợp ngày nay được áp dụng trong hầu hết các lĩnh vực khoa học, kỹ thuật và đời
sống. Việc nghiên cứu tính chất dẫn vĩ mô hay đồng nhất hóa vật liệu được nhiều nhà khoa học quan tâm
nghiên cứu và đã đưa ra nhiều kết quả xấp xỉ cho các mô hình vật liệu khác nhau. Đối với các mô hình vật
liệu trong tính toán để cho đơn giản có thể được lý tưởng hóa hình học dưới dạng cốt liệu hình cầu hoặc
trong không gian hai chiều là hình tròn. Tuy nhiên trong thực tế cốt liệu có những hình dạng phức tạp hơn
nhiều và cần xấp xỉ dưới dạng hình học không tròn mà có dạng như elliptic. Tính chất vĩ mô của vật liệu tổ
hợp phụ thuộc vào nhiều yếu tố phức tạp như cấu trúc hình học pha, các tính chất của vật liệu thành phần,
tỷ lệ thể tích giữa các pha. Do đó trong các nghiên cứu chủ yếu chỉ tìm được cận trên, dưới và các công thức
xấp xỉ áp dụng cho một số mô hình vật liệu. Hướng tiếp cận để tính xấp xỉ cho các mô hình như của
(Maxwel,1884), (Winner, 1912), (Voight, 1928), (Reuss, 1929), (Bruggeman, 1935), (Hamilton and crosser,
1962), (Lewis and Nielsen, 1970), (Mori and Tanaka, 1973). Một hướng tiếp cận khác là xây dựng biên trên
và biên dưới cho hệ số dẫn vĩ mô như (Hill, 1952), (Hashin and Strikman, 1962), (Pham DC, 1996)…Ngoài
ra các phương pháp số hiện nay cũng là cách tiếp cận hiệu quả trong việc xác định tính chất vĩ mô của vật
liệu như phương pháp phần tử hữu hạn (FEM), phương pháp biến đổi Fourier (FFT). Phương pháp FFT áp
dụng trong lĩnh vực cơ học vật liệu tính mô đun đàn hồi cho vật liệu tổ hợp được đề xuất đầu tiên vào năm
bởi (Moulicec and Subquet, 1994). Trong bài báo này sử dụng phương pháp FFT để tính hệ số dẫn vĩ mô
cho một số mô hình vật liệu hai pha với pha cốt liệu có hình học dạng elliptic được sắp xếp tuần hoàn trong
pha nền, trong đó có so sánh với các phương pháp xấp xỉ khác.
2. Phương pháp biến đổi Fourier (FFT)
Nội dung cơ bản của phương pháp biến đổi Fourier là thiết lập được phương trình LippmanSchwinger đối với bài toán không đồng nhất và sử dụng toán tử Green tuần hoàn. Sau đó sử dụng thuật toán
lặp. Ứng xử của các vật liệu thành phần được mô tả bởi định luật Fourier:
J (x) = -C(x).E(x)
(1)
Nguyễn Văn Luật, Nguyễn Trung Kiên
trong đó E(x) và J (x) lần lượt là trường gradient nhiệt độ và dòng nhiệt địa phương thỏa mãn phương
trình cân bằng
E(x) T (x), . J ( x) 0
Trường gradient E(x) và T (x) có thể tách thành các thành phần sau:
(2)
E(x) E0 e per
T (x) E0·x T per
(3)
per
trong đó E0 là gradient vĩ mô đồng nhất đối với phần tử đặc trưng, e gọi là thành phần nhiễu có tính
chất tuần hoàn. Do tính chất tuần hoàn nên ta có:
e per (x)V 0; E0 E(x)V
với ký hiệu •V là trung bình trên trên thể tích của phần tử đặc trưng V, •V
(4)
1
•dx.
V V
Bài toán trên phần tử đặc trưng có thể quy về tìm các thành phần T per , e per . Đưa vào môi trường làm chuẩn
có hệ số dẫn C0 , phương trình cân bằng trở thành
·J ·(C 0 C )E(x) 0
(5)
với C (x) C (x) C 0
Thay E(x) từ (2) vào (5) viết lại dưới dạng tương đương sau
·C 0T per ·τ (x) 0
trong đó tenxơ τ( x) gọi là tenxơ "cực" được xác định bởi : τ (x) C (x) E0 e per (x)
(6)
, e per và τ (x) được biểu diễn dưới dạng chuỗi Fourier:
1
(7)
F(x) Fˆ ( )ei .x , Fˆ ( ) F(x)ei .x F(x)ei .x dx
VV
per
per
ˆ per , eˆ per và ˆ .
trong đó F chỉ T , e và τ( x) , còn Fˆ là biến đổi Fourier của các đại lượng này, đó là T
Do tính chất chu kỳ của phần tử đặc trưng nên T
per
Ở đây tính chất tuần hoàn được thể hiện bởi
d
ei .( x h ) ei .x , h 2m j a j e j , ( j 1, 2,..., d )
(8)
j 1
với m j là số nguyên bất kỳ, 2a j là kích thước của phần tử đặc trưng song song với trục x j , e j là vectơ cở
sở theo hướng x j , d là số chiều của không gian.
k ek , j
n j
aj
Thay các biểu diễn dạng chuỗi Fourier của T
(n j 0, 1, 2...) , không tổng theo j
per
, e per và τ(x) vào phương trình (6) thu được
0
( mCmj
j )Tˆ per ( )ei .x i mˆ( )ei .x 0
per
per
từ đó các trường Tˆ và eˆ có thể xác định như sau:
(9)
Mô phỏng số FFT hệ số dẫn vĩ mô vật liệu hai pha dạng nền- cốt liệu elliptic và các phương pháp xấp xỉ
i .ˆ( ) per
.ˆ( )
Tˆ per
, eˆ ( ) i Tˆ per ( )
Γˆ 0 ( ).ˆ( )
0
0
.C
.C
0
trong đó Γ ( ) là toán tử Green phụ thuộc môi trường đồng nhất C0 được xác định bởi
Γ 0 ( )
.C 0
(10)
(11)
Từ đó thu được phương trình Lippman-Schwinger
ˆ ( ) E
ˆ 0 ( ) Γ
ˆ 0 ( ).ˆ( ) E
ˆ 0 ( ) Γ
ˆ 0 ( ). (C ( ) C 0 )* Eˆ ( )
E
(12)
Nghiệm của phương trình được tìm bởi sơ đồ lặp sau:
ˆ i 1 ( ) Γˆ 0 ( )·(C ( ) C 0 ) E
ˆ i ( ), 0
E
(13)
i 1
ˆ E0 , 0
E
Chú ý rằng Γ0·C 0 Ei ( ) Ei ( ) với 0 xem Michel(1999-[5]), phương trình (13) được viết lại
dưới dạng sau:
ˆ i 1 ( ) E
ˆ i ( ) Γˆ 0 ( ).Jˆ i ( ), 0
E
(14)
i 1
ˆ E0 , 0
E
trong đó Jˆ i ( ) là biến đổi Fourier của J i ( x) . Liên hệ giữa trường dòng J và trường gradient E trong
không gian Fourier được biểu diễn bằng biểu thức:
ˆ ( )
Jˆ ( ) C ( )* E
(15)
trong đó ký hiệu "*" là tích "convolution". Biến đổi Fourier của tenxơ hệ số dẫn:
C ( ) C (x)ei .x dx C I ( )
(16)
V
với C , I lần lượt là tenxơ hệ số dẫn và hàm dạng của pha , I ( ) được xác định theo Nemat-Nasser
(1999-[7]):
I ( )
1
ei .x dV
V V
(17)
Thay các biểu thức (15), (16) vào (14) thu được
ˆ i 1 ( ) E
ˆ i ( ) Γˆ 0 ( ). C I ( ) * E
ˆ i ( ), 0
E
i
1
0
ˆ E , 0
E
(18)
Để xác định hệ số dẫn vĩ mô của vật liệu composite, cho phần tử đặc trưng chịu tác dụng của gradient vĩ mô
E0 . Khi quá trình lặp theo (18) hội tụ (số hạng đầu tiên E1 E0 ), ta có
J( 0) C eff E0
eff
(19)
trong đó C là hệ số dẫn hiệu quả của vật liệu composite. Từ đó rút ra thuật toán số để xác định hệ số dẫn
của vật liệu nhiều thành phần có cấu trúc tuần hoàn:
Nguyễn Văn Luật, Nguyễn Trung Kiên
ˆ 1 ( ) 0 0; E
ˆ 1 (0) E0
Bước i=1: E
ˆ 1 ( )
Jˆ 1 ( ) C ( )* E
ˆ i ( ) và Jˆ i ( ) đã biết
Bước i: E
Kiểm tra hội tụ
i
ˆ i 1 ( ) E
ˆ i ( ) Γˆ 0 ( ).J ( )
E
ˆ i 1 ( )
Jˆ i 1 ( ) C ( )* E
Kiểm tra điều kiện hội tụ được xác định bằng biểu thức sau:
‖ Jˆ i 1 ( ) Jˆ i ( ) ‖
, với
‖ Jˆ i ( ) ‖
3
là sai số cho trước ( 10 )
3. Một số phương pháp xấp xỉ
Trong mục này giới thiệu một số phương pháp tính xấp xỉ, đánh giá hệ số dẫn vĩ mô (Ceff) của vật liệu
nền-cốt liệu dạng elliptic đẳng hướng trong không gian hai chiều với các ký hiệu: CI, vI là hệ số dẫn và tỉ lệ
thể tích của pha cốt liệu, CM, vM là hệ số dẫn và tỉ lệ thể tích của pha nền.
3.1.1 Xấp xỉ Hamilton (1962-[1])
C eff CM .
với n
3
CI (n 1)CM (n 1)vI (CM CI )
CI (n 1)CM vI (CM CI )
(20)
, là thông số hình học của cốt liệu, trong trường hợp cốt liệu dạng elliptic 0.4 .
3.1.2 Xấp xỉ Maxwell (1892-[1])
C eff (
vI
vM
) 1 (d 1)CM
CI (d 1)CM CM (d 1)CM
(21)
3.1.3 Xấp xỉ Lewis-Nielsen (1970-[1])
C eff
CI
1
1 m
CM
1 A vI
,
, 1
CM .
vI
2
CI
1 vI
m
A
CM
(22)
Trong đó A, m là các hệ số hình học cốt liệu (lấy giá trị A=1.5, m 1 với hình cầu, A=3, m 0.637
với hình dạng khác).
3.1.4 Xấp xỉ Mori-Tanaka (1973-[4])
Trong không gian hai chiều, cốt liệu elliptic xấp xỉ Mori-Tanaka có dạng:
C eff CM vI (CI CM ).
(23)
Mô phỏng số FFT hệ số dẫn vĩ mô vật liệu hai pha dạng nền- cốt liệu elliptic và các phương pháp xấp xỉ
trong đó
CM (CI CM )(1 r )2
, r là tỉ số giữa hai bán trục chính của elliptic.
2(CI rCM )(CM rCI )
3.1.5 Đánh giá Hashin-Strikman (1962-[3])
Hashin-Strikman dựa trên nguyên lý biến phân riêng đưa vào trường khả dĩ phân cực đã xây dựng
được đánh giá trên và dưới cho hệ số dẫn vĩ mô của vật liệu nhiều thành phần đẳng hướng.
PC ((d 1)Cmin ) C eff PC ((d 1)Cmax ),
(24)
1
v
PC (C0 )
C* ,
C C*
C* (d 1)C0 , Cmin min C1 , , Cn , Cmax max C1 , , Cn
4. Kết quả so sánh
Trong mục này sẽ đưa ra kết quả tính toán FFT hệ số dẫn vĩ mô cho một số mô hình vật liệu có cốt
liệu dạng elliptic trong không gian hai chiều và so sánh với các phương pháp xấp xỉ khác. Hàm dạng (17)
cho cốt liệu elliptic trong không gian hai chiều [2]
(25)
I ( ) 2S J1 ( ) / .ei .xc ( )
Trong đó J1 là hàm Bessel loại 1, S là diện tích bề mặt của cốt liệu, xc ( ) là véc tơ xác định vị trí trọng
tâm của pha cốt liệu
( 12 a12 22 a22 )1/2
a1 , a2 là độ dài các bán trục chính của cốt liệu elliptic, 1 , 2 là các thành phần của theo các trục của
elliptic.
- Kết quả tính toán FFT cho mô hình vật liệu có cốt liệu elliptic sắp xếp tuần hoàn trong phần tử đặc
trưng (unit cell) theo hình vuông (square) và so sánh với các mô hình xấp xỉ khác. Do cốt liệu không chồng
lấn nên tỉ lệ thể tích cốt liệu chỉ có thể tăng đến một giới hạn nhất định. Trên hình 2 với pha nền có hệ số
dẫn CM=1, pha cốt liệu có hệ số dẫn CI=10, a2 0.2a1 , trong trường hợp này cả FFT và các xấp xỉ khác
đều nằm trong đánh giá của Hashin-Strickman (HS) và gần nhau khi tỉ lệ thể tích cốt liệu nhỏ. Kết quả trên
hình 3 với CM=10, CI=2 thì có thể thấy chỉ có FFT nằm trong đánh giá của HS, điều này khẳng định độ tin
cậy của phương pháp FFT so với các phương pháp xấp xỉ trước đó.
Hình 1: Mô hình cốt liệu elliptic có cấu trúc tuần hoàn