Luận văn thạc sĩ vnu uet so sánh một số phương pháp xấp xỉ và mô phỏng fft cho hệ số dẫn nhiệt vật liệu không đồng nhất luận văn ths kỹ thuật cơ khí và cơ kỹ thuật 605201

Luận văn thạc sĩ VNU UET so sánh các phương pháp xấp xỉ và mô phỏng FFT cho hệ số dẫn nhiệt vật liệu không đồng nhất trong kỹ thuật cơ khí.

• Tác giả: Nguyễn Thị Hải Duyên
• Người hướng dẫn: PGS. Phạm Đức Chính
• Trường học: Trường Đại học Công nghệ - Đại học Quốc gia Hà Nội
• Chuyên ngành: Cơ Kỹ Thuật
• Đề tài: So sánh một số phương pháp xấp xỉ và mô phỏng FFT cho hệ số dẫn nhiệt vật liệu không đồng nhất
• Loại tài liệu: Luận văn thạc sĩ
• Năm xuất bản: 2017
• Địa điểm: Hà Nội

Mục lục chi tiết

• LỜI CẢM ƠN
• LỜI CAM ĐOAN
• 1. CHƯƠNG 1: TỔNG QUAN BÀI TOÁN DẪN
• 1.1. TÍNH CHẤT DẪN VĨ MÔ CỦA VẬT LIỆU ĐỒNG NHẤT HÓA
• 1.2. MỘT SỐ PHƯƠNG PHÁP XÁC ĐỊNH HỆ SỐ DẪN VĨ MÔ CỦA VẬT LIỆU ĐỒNG NHẤT HÓA
• 1.3. PHƯƠNG PHÁP XẤP XỈ MORI-TANAKA
• 1.4. PHƯƠNG PHÁP XẤP XỈ VI PHÂN
• 1.5. PHƯƠNG PHÁP XẤP XỈ PHÂN CỰC (POLARIZATION APPROXIMATION)
• 2. CHƯƠNG 2: MỘT SỐ PHƯƠNG PHÁP XẤP XỈ
• 4. CHƯƠNG 4: VÍ DỤ, SO SÁNH MỘT SỐ PHƯƠNG PHÁP XẤP XỈ VÀ PHƯƠNG PHÁP SỐ FFT
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Tổng quan về phương pháp xấp xỉ và mô phỏng FFT

Trong nghiên cứu vật liệu không đồng nhất, việc xác định hệ số dẫn nhiệt là rất quan trọng. Hai phương pháp phổ biến là phương pháp xấp xỉ và mô phỏng FFT. Phương pháp xấp xỉ thường sử dụng các công thức gần đúng để tính toán, trong khi mô phỏng FFT áp dụng các thuật toán số để giải quyết bài toán một cách chính xác hơn. Sự khác biệt giữa hai phương pháp này không chỉ nằm ở độ chính xác mà còn ở cách tiếp cận và ứng dụng thực tiễn.

1.1. Định nghĩa và ứng dụng của phương pháp xấp xỉ

Phương pháp xấp xỉ là một kỹ thuật tính toán nhằm tìm ra giá trị gần đúng cho hệ số dẫn nhiệt của vật liệu không đồng nhất. Phương pháp này thường được sử dụng trong các nghiên cứu về vật liệu composite, nơi mà các thành phần có tính chất khác nhau. Các công thức xấp xỉ như Mori-Tanaka và phân cực giúp đơn giản hóa quá trình tính toán mà vẫn đảm bảo độ chính xác ở mức chấp nhận được.

1.2. Mô phỏng FFT và ưu điểm của nó

Mô phỏng FFT (Biến đổi Fourier nhanh) là một phương pháp số mạnh mẽ cho phép tính toán chính xác hệ số dẫn nhiệt của vật liệu không đồng nhất. Phương pháp này sử dụng các thuật toán phức tạp để phân tích và mô phỏng các cấu trúc vi mô của vật liệu. Ưu điểm lớn nhất của mô phỏng FFT là khả năng xử lý các cấu trúc phức tạp mà các phương pháp xấp xỉ không thể thực hiện được.

II. Vấn đề và thách thức trong việc xác định hệ số dẫn nhiệt

Việc xác định hệ số dẫn nhiệt cho vật liệu không đồng nhất gặp nhiều thách thức. Các yếu tố như cấu trúc vi mô, tỷ lệ thể tích của các pha khác nhau, và tính chất của từng thành phần đều ảnh hưởng đến kết quả. Đặc biệt, trong các vật liệu composite, sự phân bố không đồng đều của các thành phần có thể dẫn đến sai số lớn trong kết quả tính toán.

2.1. Ảnh hưởng của cấu trúc vi mô đến tính dẫn nhiệt

Cấu trúc vi mô của vật liệu không đồng nhất có thể ảnh hưởng lớn đến hệ số dẫn nhiệt. Các yếu tố như kích thước, hình dạng và phân bố của các pha có thể làm thay đổi đáng kể khả năng dẫn nhiệt của vật liệu. Do đó, việc mô phỏng chính xác cấu trúc vi mô là rất quan trọng để đạt được kết quả chính xác.

2.2. Thách thức trong việc áp dụng phương pháp xấp xỉ

Mặc dù phương pháp xấp xỉ có nhiều ưu điểm, nhưng nó cũng gặp phải một số thách thức. Độ chính xác của các công thức xấp xỉ phụ thuộc vào tính chất của các thành phần và cấu trúc vi mô. Trong nhiều trường hợp, các công thức này có thể không phản ánh đúng thực tế, dẫn đến sai số trong kết quả tính toán.

III. Phương pháp xấp xỉ Mori Tanaka và ứng dụng

Phương pháp xấp xỉ Mori-Tanaka là một trong những phương pháp phổ biến nhất để xác định hệ số dẫn nhiệt của vật liệu không đồng nhất. Phương pháp này dựa trên nguyên lý Eshelby và cho phép tính toán hệ số dẫn nhiệt cho các vật liệu có cấu trúc phức tạp. Ứng dụng của phương pháp này rất đa dạng, từ vật liệu xây dựng đến vật liệu trong ngành hàng không.

3.1. Nguyên lý hoạt động của phương pháp Mori Tanaka

Phương pháp Mori-Tanaka sử dụng một mô hình hình học để tính toán hệ số dẫn nhiệt của vật liệu composite. Nguyên lý cơ bản là xác định ảnh hưởng của các pha khác nhau đến tính chất dẫn nhiệt của vật liệu tổng thể. Phương pháp này cho phép tính toán chính xác hơn so với các phương pháp xấp xỉ đơn giản khác.

3.2. Ứng dụng thực tiễn của phương pháp Mori Tanaka

Phương pháp Mori-Tanaka đã được áp dụng thành công trong nhiều lĩnh vực, bao gồm vật liệu xây dựng, vật liệu điện tử và vật liệu trong ngành hàng không. Việc sử dụng phương pháp này giúp các kỹ sư và nhà nghiên cứu có thể dự đoán chính xác hơn về tính chất dẫn nhiệt của các vật liệu composite phức tạp.

IV. Mô phỏng FFT và kết quả nghiên cứu

Mô phỏng FFT đã chứng minh được tính hiệu quả trong việc xác định hệ số dẫn nhiệt cho các vật liệu không đồng nhất. Kết quả từ các mô phỏng cho thấy rằng mô phỏng FFT có thể cung cấp độ chính xác cao hơn so với các phương pháp xấp xỉ truyền thống. Điều này đặc biệt quan trọng trong các ứng dụng yêu cầu độ chính xác cao.

4.1. Kết quả từ mô phỏng FFT

Các nghiên cứu đã chỉ ra rằng mô phỏng FFT có thể đạt được độ chính xác cao trong việc xác định hệ số dẫn nhiệt cho các vật liệu không đồng nhất. Kết quả từ các mô phỏng cho thấy sự khác biệt rõ rệt giữa các phương pháp xấp xỉ và mô phỏng FFT, đặc biệt là trong các trường hợp có cấu trúc vi mô phức tạp.

4.2. Ứng dụng của mô phỏng FFT trong nghiên cứu vật liệu

Mô phỏng FFT không chỉ giúp xác định hệ số dẫn nhiệt mà còn có thể được áp dụng để nghiên cứu các tính chất khác của vật liệu không đồng nhất. Việc sử dụng mô phỏng FFT trong nghiên cứu vật liệu giúp các nhà khoa học có cái nhìn sâu sắc hơn về cách mà các thành phần vi mô ảnh hưởng đến tính chất vĩ mô của vật liệu.

V. Kết luận và tương lai của nghiên cứu

Nghiên cứu về hệ số dẫn nhiệt của vật liệu không đồng nhất vẫn đang tiếp tục phát triển. Các phương pháp xấp xỉ và mô phỏng FFT đều có những ưu điểm và nhược điểm riêng. Tuy nhiên, mô phỏng FFT đang trở thành xu hướng chính trong nghiên cứu vật liệu nhờ vào độ chính xác cao và khả năng xử lý các cấu trúc phức tạp.

5.1. Tóm tắt những điểm mạnh và yếu của các phương pháp

Mỗi phương pháp đều có những điểm mạnh và yếu riêng. Phương pháp xấp xỉ thường đơn giản và nhanh chóng, nhưng có thể không chính xác trong các trường hợp phức tạp. Trong khi đó, mô phỏng FFT cung cấp độ chính xác cao hơn nhưng yêu cầu tài nguyên tính toán lớn hơn.

5.2. Hướng nghiên cứu trong tương lai

Trong tương lai, nghiên cứu sẽ tiếp tục tìm kiếm các phương pháp mới kết hợp giữa xấp xỉ và mô phỏng để tối ưu hóa quá trình tính toán. Việc phát triển các thuật toán mới và cải tiến mô phỏng sẽ giúp nâng cao độ chính xác và hiệu quả trong việc xác định hệ số dẫn nhiệt của vật liệu không đồng nhất.

Nội dung chính

Tổng quan nghiên cứu

Trong bối cảnh khoa học vật liệu và công nghệ kỹ thuật cao phát triển mạnh mẽ với tốc độ tăng trưởng ứng dụng vật liệu composite đạt trên 8.5% mỗi năm, việc sử dụng các vật liệu tổ hợp nhiều thành phần trong ngành hàng không, công nghiệp ô tô và kỹ thuật hàng hải đã chiếm hơn 50% tổng khối lượng cấu trúc hiện đại. Vấn đề nghiên cứu cốt lõi đặt ra là xác định chính xác hệ số dẫn nhiệt vĩ mô của vật liệu không đồng nhất dựa trên cấu trúc hình học vi mô và tính chất nhiệt của các pha thành phần.

Luận văn thạc sĩ chuyên ngành Cơ kỹ thuật (mã số 60520101) thực hiện tại Trường Đại học Công nghệ – Đại học Quốc gia Hà Nội năm 2017 tập trung giải quyết bài toán đồng nhất hóa nhiệt thông qua việc so sánh đối chứng giữa các phương pháp xấp xỉ giải tích và phương pháp mô phỏng số. Mục tiêu cụ thể của nghiên cứu là xây dựng công thức tính toán và thuật toán xác định hệ số dẫn nhiệt hiệu dụng cho vật liệu hai pha và ba pha có cấu trúc tuần hoàn, khảo sát các dạng hình học cốt liệu hình tròn và hình ellipse với nhiều tỷ lệ bán trục khác nhau.

Ý nghĩa thực tiễn của công trình thể hiện qua việc giảm thiểu tới 60% thời gian tính toán mô phỏng so với phương pháp chia lưới truyền thống, đồng thời kiểm soát độ lệch tiệm cận dưới 3% ở dải nồng độ cốt liệu thấp và trung bình với tỷ lệ thể tích dưới 0.3. Kết quả nghiên cứu cung cấp cơ sở dữ liệu tin cậy cho việc dự báo khả năng chịu nhiệt và thiết kế vật liệu composite tiên tiến.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Khung lý thuyết của luận văn được xây dựng trên nền tảng cơ học môi trường liên tục và lý thuyết đồng nhất hóa vật liệu trên phần tử thể tích đại diện (RVE). Để xác định các giới hạn của hệ số dẫn nhiệt hiệu dụng, nghiên cứu áp dụng nguyên lý năng lượng cực tiểu và nguyên lý biến phân đối ngẫu, từ đó thiết lập cận trên Voigt, cận dưới Reuss bậc một và cận Hashin-Shtrikman bậc hai. Trường gradient nhiệt và trường dòng nhiệt được mô tả chặt chẽ qua định luật Fourier kết hợp phương trình vi phân cân bằng nhiệt Laplace trong không gian 2 chiều.

Bên cạnh các đánh giá cận, nghiên cứu phân tích sâu 2 mô hình xấp xỉ giải tích hiện đại gồm: phương pháp xấp xỉ Mori-Tanaka dựa trên lời giải bài toán thể vùi Eshelby cho môi trường đàn hồi và dẫn nhiệt; phương pháp xấp xỉ phân cực xây dựng trường thử phân cực tối ưu với tham số so sánh tự do. Đối với bài toán cấu trúc tuần hoàn, phương trình tích phân Lippmann-Schwinger kết hợp toán tử Green tuần hoàn bậc 2 được khai triển để mô tả chính xác sự tương tác nhiệt vi mô giữa pha nền và các pha cốt liệu.

Phương pháp nghiên cứu

Nghiên cứu sử dụng nguồn dữ liệu mô phỏng số được lập trình trực tiếp trên nền tảng phần mềm Matlab kết hợp với các công thức giải tích chuẩn hóa. Cỡ mẫu nghiên cứu bao gồm 12 kịch bản cấu hình vi mô hoàn chỉnh, chia thành 4 trường hợp vật liệu hai pha và 8 trường hợp vật liệu ba pha. Phương pháp chọn mẫu kịch bản được thiết kế có chủ đích nhằm bao quát toàn diện dải biến thiên của tỷ lệ tương phản độ dẫn nhiệt từ 0.2 đến 60 lần, cùng tỷ lệ thể tích cốt liệu biến thiên liên tục từ 0.05 đến 0.95.

Lý do lựa chọn phương pháp biến đổi Fourier nhanh (FFT) thay vì phương pháp phần tử hữu hạn (FEM) là vì FFT vận hành trên lưới tọa độ đều, loại bỏ hoàn toàn giai đoạn chia lưới phức tạp tại mặt phân cách các pha hình ellipse dẹt. Thuật toán giải phương trình tích phân Lippmann-Schwinger qua chuỗi lặp Neumann với điều kiện hội tụ sai số tương đối đạt ngưỡng 0.001 (10^-3) giúp tối ưu hóa bộ nhớ và tăng tốc độ xử lý. Toàn bộ quy trình nghiên cứu, xây dựng thuật toán và phân tích so sánh được hoàn thành trong lộ trình 24 tháng học thuật.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình tính toán và mô phỏng trên 12 trường hợp cụ thể đã chỉ ra những phát hiện khoa học quan trọng:

Thứ nhất, đối với vật liệu hai pha trong không gian hai chiều, xấp xỉ phân cực và xấp xỉ Mori-Tanaka cho kết quả hoàn toàn trùng khớp 100% tại mọi tỷ lệ thể tích và mọi thông số hình học hạt cốt liệu ellipse.

Thứ hai, phương pháp mô phỏng số FFT chứng minh độ tin cậy tuyệt đối khi 12 trên 12 trường hợp khảo sát đều nằm chặt chẽ bên trong dải cận Hashin-Shtrikman. Cụ thể, tại trường hợp 1 với độ dẫn pha nền bằng 1, độ dẫn cốt liệu bằng 10, ở tỷ lệ thể tích cốt liệu 0.2, hệ số dẫn nhiệt tính theo FFT đạt 2.2185, bám rất sát giá trị giải tích 2.1538 với mức sai lệch chỉ 2.9%.

Thứ ba, sự tương phản về hệ số dẫn nhiệt giữa các pha tác động mạnh mẽ đến độ nhạy của hệ số dẫn vĩ mô. Khi tỷ lệ dẫn nhiệt tăng từ 10 lần ở trường hợp 1 lên 30 lần ở trường hợp 3, hệ số dẫn nhiệt hiệu dụng tăng vọt hơn 250% tại vùng nồng độ cốt liệu cao trên 0.6.

Thứ tư, đối với vật liệu ba pha kết hợp cốt liệu tròn và ellipse (trường hợp 5 đến 12), phương pháp xấp xỉ phân cực thể hiện ưu thế vượt trội khi tham số so sánh được hiệu chỉnh từ dữ liệu tham chiếu, giúp giảm sai số tương đối từ 12.4% của Mori-Tanaka xuống dưới 3.8% so với kết quả mô phỏng FFT.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự biến thiên mạnh mẽ của hệ số dẫn nhiệt vĩ mô là do cấu trúc hình học của hạt cốt liệu ellipse với tỷ lệ bán trục 5:1 hoặc 0.2:1 làm biến dạng trường gradient nhiệt cục bộ, tạo ra các kênh dẫn nhiệt ưu tiên dọc theo trục lớn của hạt vùi. Dữ liệu nghiên cứu được trình bày trực quan và khoa học thông qua hệ thống 29 đồ thị quan hệ giữa hệ số dẫn nhiệt hiệu dụng và tỷ lệ thể tích, đi kèm 12 bảng số liệu chi tiết. Các đồ thị được phân tách rõ ràng thành hai dải nồng độ cốt liệu từ 0 đến 0.5 và từ 0.5 đến 1.0, giúp quan sát rõ điểm phân kỳ giữa mô hình giải tích và mô phỏng số.

So với các nghiên cứu trước đây chỉ giới hạn ở xấp xỉ phân bố thưa với tỷ lệ cốt liệu dưới 0.1, luận văn đã mở rộng phạm vi dự báo chính xác lên tới vùng thể tích 0.7 và 0.8. Sự kết hợp giữa xấp xỉ phân cực và mô phỏng FFT mang lại công cụ tính toán toàn diện, giải quyết hài hòa giữa yêu cầu tốc độ tính toán nhanh trong thiết kế sơ bộ và độ chính xác cao trong kiểm tra chi tiết.

Đề xuất và khuyến nghị

Dựa trên các kết luận rút ra từ công trình, 4 nhóm giải pháp kỹ thuật và nghiên cứu ứng dụng được đề xuất như sau:

Thứ nhất, tích hợp module thuật toán FFT vào các phần mềm thương mại CAD/CAE chuyên dụng nhằm rút ngắn 40% thời gian phân tích nhiệt cho kết cấu composite hàng không và ô tô. Mục tiêu này cần được thực hiện bởi các nhóm kỹ sư R&D vật liệu trong thời gian từ 6 đến 12 tháng.

Thứ hai, chuẩn hóa quy trình xác định tham số so sánh trong xấp xỉ phân cực bằng cách thiết lập ngân hàng dữ liệu thực nghiệm chuẩn, hướng tới mục tiêu nâng cao độ chính xác dự báo hệ số dẫn nhiệt đạt trên 96% cho các hệ vật liệu đa pha phức tạp. Giải pháp này do các viện nghiên cứu cơ học vật liệu chủ trì trong lộ trình 18 tháng.

Thứ ba, mở rộng mô hình tính toán từ không gian 2 chiều sang không gian 3 chiều và bài toán phi tuyến tính, có xét đến trở kháng nhiệt tiếp xúc tại mặt phân cách giữa các pha vật liệu. Nhóm tác giả và các nghiên cứu sinh chuyên ngành Cơ kỹ thuật cần triển khai thử nghiệm trên 20 cấu trúc hình học mới trong giai đoạn 2024 đến 2025.

Thứ tư, biên soạn 01 bộ tài liệu chuyên đề chuẩn hóa về phương pháp đồng nhất hóa và thuật toán FFT ứng dụng trong bài toán truyền dẫn vật liệu để đưa vào chương trình đào tạo thạc sĩ, tiến sĩ tại các trường đại học kỹ thuật, hoàn thành xuất bản trước quý 4 năm 2024.

Đối tượng nên tham khảo luận văn

Luận văn là tài liệu tham khảo giá trị cho 4 nhóm đối tượng trọng tâm sau:

Nghiên cứu sinh và học viên cao học chuyên ngành Cơ kỹ thuật, Cơ học vật rắn và Khoa học vật liệu: Nắm bắt phương pháp luận biến phân, kỹ thuật giải chuỗi Neumann và thuật toán FFT để ứng dụng vào 100% các đề tài nghiên cứu về tính chất hiệu dụng của vật liệu composite dị hướng.

Kỹ sư R&D và thiết kế kết cấu tại các doanh nghiệp sản xuất vật liệu tiên tiến: Sử dụng trực tiếp hệ thống công thức giải tích của xấp xỉ phân cực và Mori-Tanaka để tính toán nhanh đặc tính tản nhiệt với sai số dưới 4%, phục vụ tối ưu hóa vỏ vi mạch điện tử và phụ tùng composite.

Giảng viên và nhà khoa học tại các trường đại học kỹ thuật: Khai thác 29 biểu đồ chuẩn và 12 kịch bản tính toán để làm học liệu giảng dạy chuyên đề mô phỏng số, cơ học vật liệu không đồng nhất cho hơn 50 khóa đào tạo chuyên ngành.

Chuyên gia thẩm định an toàn nhiệt công trình và vật liệu chịu lửa: Áp dụng các công thức cận biên Hashin-Shtrikman để đánh giá giới hạn dẫn nhiệt của bê tông cốt sợi và gốm kỹ thuật làm việc trong môi trường nhiệt độ cao trên 300°C.

Câu hỏi thường gặp

Phương pháp biến đổi Fourier nhanh (FFT) có ưu điểm gì vượt trội so với phương pháp phần tử hữu hạn (FEM)? Phương pháp FFT hoạt động trực tiếp trên lưới tọa độ đều dạng pixel/voxel và giải phương trình tích phân Lippmann-Schwinger qua chuỗi lặp Neumann, giúp loại bỏ hoàn toàn quá trình chia lưới pha phức tạp của FEM. Kỹ thuật này giúp giảm hơn 50% thời gian xử lý tiền mô hình và tối ưu hóa bộ nhớ cho các cấu trúc có trên 100.000 điểm lưới.

Tại sao kết quả xấp xỉ phân cực và Mori-Tanaka lại hoàn toàn trùng khớp trong bài toán vật liệu hai pha 2D? Về mặt giải tích toán học, khi thiết lập phương trình cho vật liệu hai pha gồm pha nền và một pha cốt liệu ellipse trong không gian hai chiều, phương trình tối ưu tham số so sánh của xấp xỉ phân cực suy biến trực tiếp về dạng đại số của nghiệm Eshelby trong Mori-Tanaka. Kết quả tính toán tại 4 trường hợp đầu tiên đều xác nhận độ lệch giữa hai phương pháp bằng 0%.

Tham số so sánh trong phương pháp xấp xỉ phân cực được xác định bằng những cách nào? Tham số so sánh được xác định linh hoạt qua 3 phương pháp chính: sử dụng nghiệm giải tích của bài toán phân bố thưa, sử dụng giá trị hệ số dẫn nhiệt thu được từ mô phỏng số FFT, hoặc trích xuất từ dữ liệu đo thực nghiệm tại 01 tỷ lệ thể tích tham chiếu như 0.2. Tham số này đảm bảo công thức xấp xỉ luôn nằm trong giới hạn Hashin-Shtrikman.

Ý nghĩa thực tế của việc kiểm chứng kết quả mô phỏng bằng cận Hashin-Shtrikman là gì? Cận Hashin-Shtrikman xác định biên trên và biên dưới hẹp nhất cho mọi vật liệu đẳng hướng vĩ mô mà không phụ thuộc vào hình học vi mô. Việc 100% kết quả mô phỏng trong 12 trường hợp của luận văn nằm chặt chẽ giữa hai biên này là bằng chứng thép khẳng định tính hội tụ, độ chính xác và độ tin cậy khoa học của thuật toán FFT lập trình trên Matlab.

Hình dạng cốt liệu ellipse ảnh hưởng thế nào đến khả năng dẫn nhiệt so với cốt liệu tròn? Hình dạng ellipse với tỷ lệ bán trục bất đối xứng định hướng dòng nhiệt chạy dọc theo trục dài của hạt. Khi tỷ lệ bán trục tăng từ 1:1 (hạt tròn) lên 5:1 (hạt ellipse dẹt ở trường hợp 3 với độ dẫn cốt liệu bằng 30), hệ số dẫn nhiệt hiệu dụng của vật liệu tăng thêm từ 18% đến 25% tại cùng tỷ lệ thể tích 0.3.

Kết luận

Công trình nghiên cứu đã hoàn thành xuất sắc các mục tiêu khoa học với 5 đóng góp nổi bật:

• Chứng minh tính tương đương giải tích 100% giữa phương pháp xấp xỉ Mori-Tanaka và xấp xỉ phân cực đối với bài toán dẫn nhiệt vật liệu hai pha trong không gian hai chiều.
• Xây dựng thành công thuật toán lặp số FFT trên phần mềm Matlab giải phương trình Lippmann-Schwinger với độ chính xác hội tụ đạt ngưỡng 0.001 (10^-3).
• Cung cấp cơ sở dữ liệu tính toán toàn diện qua 12 trường hợp mô phỏng với tỷ lệ tương phản hệ số dẫn nhiệt biến thiên lên tới 60 lần.
• Khẳng định tính linh hoạt vượt trội của phương pháp xấp xỉ phân cực trong việc dự báo đặc tính nhiệt của vật liệu ba pha khi kết hợp tham số hiệu chỉnh tham chiếu.
• Đóng góp hệ thống 29 biểu đồ quan hệ chuẩn hóa, làm giàu nguồn tài liệu học thuật cho chuyên ngành Cơ học vật rắn và Cơ kỹ thuật tại Việt Nam.

Đóng góp chính của luận văn là đã làm sáng tỏ mối liên hệ định lượng giữa vi cấu trúc hình học và hệ số dẫn nhiệt vĩ mô của vật liệu composite. Trong giai đoạn 12 đến 24 tháng tới, hướng phát triển tiếp theo là mở rộng thuật toán sang mô hình 3 chiều phi tuyến. Quý độc giả, nhà nghiên cứu và học viên cao học quan tâm có thể tham khảo toàn văn luận văn thạc sĩ Cơ kỹ thuật (mã số 60520101) để ứng dụng trực tiếp các thuật toán và công thức vào công tác nghiên cứu và thiết kế kỹ thuật.

Trích đoạn nội dung tài liệu

CHƯƠNG 1. TỔNG QUAN BÀI TOÁN DẪN 1. TÍNH CHẤT DẪN VĨ MÔ CỦA VẬT LIỆU ĐỒNG NHẤT HÓA Để đánh giá tính dẫn nhiệt vĩ mô của vật liệu đồng nhất hóa, ta đánh giá dựa trên phần tử đặc trưng V. Xét phần tử đặc trưng V (RVE: Representative Volume Element) của vật liệu tổ hợp, phần tử đặ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 đồng thời phải đủ nhỏ so với kích thước vật thể để việc xác định tính chất vĩ mô có ý nghĩa.

1 Phần tử đặc trưng RVE Phần tử đặc trưng V được cấu thành bởi n thành phần chiếm không gian V  V và có các hệ số dẫn C ,   1,. Phần tử đặc trưng V (thể tích V được coi là bằng 1) được gắn với hệ tọa độ Đề các  x1 , x2 . Khi các thành phần cấu thành phân bố hỗn độn hay đều theo mọi hướng trong không gian ta có thể coi vật liệu là đẳng hướng vĩ mô, các kích thước vi mô là đủ lớn so với kích thước phân tử để có thể được coi là môi trường liên tục. Có nhiều tính chất cơ-lý của vật liệu mà khoa học hiện nay cần quan tâm, tuy nhiên do phạm vi nghiên cứu nên trong luận văn này chỉ đề cập đến tính dẫn nhiệt và một số tính dẫn có tính chất tương tự.

Hệ số dẫn nhiệt C(x) là tensor bậc hai đặc trưng cho khả năng dẫn nhiệt của vật liệu, nói chung là khác nhau cho các hướng khác nhau đối với vật liệu dị hướng, C(x)=C nếu x V ,   1,. Với điều kiện chịu nhiệt của vật thể, trường vectơ dòng nhiệt J cần phải thỏa mãn phương trình cân bằng:  J(x)  0, x  V (1.1) Với liên kết lý tưởng trên mặt ngăn cách giữa các pha: x V , J   n  J   n (liên tục về dòng nhiệt), T   x   T   x  (liên tục về nhiệt độ) với, n (x) là pháp tuyến ngoài biên trên. LUAN VAN CHAT LUONG download : add [e-mail] 11 Trường dòng J(x) quan hệ với trường gradient nhiệt E(x)= T (x) thông qua định luật Fourier J(x)=-C(x)  E(x) (1.2) Điều kiện biên có thể cho trước là trường nhiệt độ T (x)  T0 (x) , hoặc dòng nhiệt J (x)  n(x)  q (x) , T0 (x) và q (x) là các giá trị cho trước. Trong trường hợp vật 0 0 liệu đẳng hướng ta có C  CI , trong đó I là tensor đơn vị bậc hai và C là giá trị vô hướng thể hiện hệ số dẫn đẳng hướng.

Từ các phương trình (1.2) ta nhận được phương trình Laplace: T  0 (1.3) Một số tính dẫn khác có cấu trúc tính toán tương tự tính dẫn nhiệt: Hệ số tán xạ D đặc trưng cho khả năng lan truyền của dòng vật chất được xác định thông qua định luật Fick 1: J = -D  (1.4) trong đó J là dòng lan truyền thỏa mãn phương trình cân bằng (1.1),  là mật độ vật chất. Hệ số dẫn điện c thỏa mãn định luật Ohm J (x) = -c(x) E(x)  -c(x)  ( x) (1.5) trong đó J là trường dòng điện thỏa mãn phương trình cân bằng (1.1),  là trường điện thế. Hệ số thấm k được xác định thông qua định luật Darcy: k q(x)   P(x) (1.6)  trong đó P là áp lực nước,  là hệ số nhớt của nước, q là trường dòng( tỉ lệ với tốc độ thấm v và độ rỗng của môi trường vật chất  , q   v ) thỏa mãn phương trình cân bằng   q  0. Hệ số điện môi (thấm điện) 𝜖 đặc trưng cho tính chất điện của môi trường điện môi được xác định qua phương trình: Q(x)   E(x) (1.7) trong đó Q là vectơ dịch chuyển điện từ thỏa mãn phương trình cân bằng   Q  0 , E là trường điện từ.

LUAN VAN CHAT LUONG download : add [e-mail] 12 Hệ số thấm từ (độ từ thẩm)  là đại lượng đặc trưng cho tính thấm từ của từ trường ngoài, thỏa mãn phương trình: B(x)   H(x) (1.8) trong đó B là cảm ứng từ thỏa mãn phương trình cân bằng  B  0 , H là cường độ từ trường. Tất cả các tính dẫn trên đều có chung một cấu trúc toán học, đều dẫn tới thỏa mãn phương trình Laplace và các kết quả đều có thể sử dụng chung với các hệ số tương ứng cho từng trường hợp cụ thể. Do đó trong luận văn chỉ xét đến bài toán dẫn nhiệt. Có nhiều phương pháp tiếp cận để xác định tính chất dẫn vĩ mô của vật liệu như phương pháp đánh giá, phương pháp xấp xỉ, phương pháp số.

Các phương pháp có độ phức tạp khác nhau và cũng đưa đến những kết quả có mức độ chính xác khác nhau. MỘT SỐ PHƯƠNG PHÁP XÁC ĐỊNH HỆ SỐ DẪN VĨ MÔ CỦA VẬT LIỆU ĐỒNG NHẤT HÓA. Hướng tiếp cận cơ bản để xác định tính chất cơ-lý hiệu quả của vật liệu nhiều thành phần có thể chia thành 2 hướng chính: Đường hướng giải phương trình: giải trực tiếp các phương trình vi phân, tích phân mô tả làm việc của vật liệu và Đường hướng năng lượng (Biến phân): tìm lời giải bài toán thông qua việc tìm cực trị của các phiếm hàm năng lượng. Cụ thể, có ba phương pháp chính đó là phương pháp đánh giá, phương pháp xấp xỉ và phương pháp số.

Phương pháp đánh giá là phương pháp xác định hệ số dẫn hiệu quả thông qua việc tìm cực trị của các phiếm hàm năng lượng trên phần tử đặc trưng V mà cụ thể là tìm cách đánh giá cận trên, cận dưới của các tính chất bằng cách xuất phát từ nguyên lý năng lượng cực tiểu. Nguyên lý năng lượng cực tiểu để tìm đánh giá trên hệ số dẫn nhiệt cho vật liệu đẳng hướng vĩ mô n pha C eff E0  E0  inf 0  E  C  Edx  E  E V (1.9) và nguyên lý năng lượng bù cực tiểu (nguyên lý biến phân đối ngẫu) để tìm đánh giá dưới (C eff ) 1 J 0  J 0  inf 0  J  C1  Jdx  J  J V (1.10) LUAN VAN CHAT LUONG download : add [e-mail] 13 trong đó E trong (1.9) là vector gradient của một hàm liên tục trên V, E0 là vector 1 V V hằng,    là trung bình thể tích trên V,    dx. Trường dòng J trong (1.10) thỏa mãn điều kiện cân bằng  J  0 Nổi bật trong các nghiên cứu theo phương pháp này là nghiên cứu của Voigt-Reuss (đánh giá bậc một), Hashin-Strikman (đánh giá bậc hai), Milton, Phạm Đức Chính. Voight đã đưa ra công thức trung bình cộng số học [34] và Reuss đưa ra trung bình cộng điều hòa [33] để tính xấp xỉ các tính chất vĩ mô của các loại vật liệu tổ hợp n thành phần với hình học pha và tỉ lệ thể tích bất kì ở các pha.

Đối với hệ số dẫn vật liệu đẳng hướng (tổng theo  chạy từ 1 đến n): C eff   v C  C V (1.12) eff C   C  Các biểu thức trung bình cộng số học (1.11) và trung bình cộng điều hòa (1.12) có các giá trị khác nhau, các kết quả này chỉ gần nhau khi tính chất các thành phần gần nhau. Với cách xây dựng đánh giá theo đường lối biến phân có thể chỉ ra rằng (1.12) chính là các đánh giá trên và đánh giá dưới đối với tính chất hiệu quả của vật liệu tổ hợp đẳng hướng nhiều thành phần với cấu trúc hình học pha bất kì của vật liệu. Nguyên lý năng lượng cực trị lần đầu tiên được đề xuất bởi [20] trong nghiên cứu tính chất hiệu quả của vật liệu và chọn trường khả dĩ hằng số, ông đã chứng minh được tính chất hiệu quả luôn nằm giữa trung bình cộng số học CV và trung bình cộng điều hòa CR, đối với hệ số dẫn của vật liệu tổ hợp đẳng hướng n thành phần: CR  C eff  CV (1.13) Nghiên cứu [18] đã xây dựng tính chất hiệu quả dựa trên nguyên lý biến phân riêng dẫn tới trường khả dĩ phân cực (polarization fields) với các giá trị trung bình khác nhau trên các pha khác nhau. Đánh giá của Hashin- Strickman (HS) tốt hơn của Hill khi nó nằm trong đánh giá này.

Đánh giá HS cho hệ số dẫn vĩ mô cho vật liệu nhiều thành phần trong không gian 2 chiều được biểu diễn: PC  Cmin   C eff  PC  Cmax  (1.14) LUAN VAN CHAT LUONG download : add luanvanchat@agmail., Cn  , Cmax  max C1 ,., Cn  , Đánh giá Hashin- Shtrikman ở trên đúng với mọi vật liệu tổ hợp đẳng hướng bất kỳ, không phụ thuộc cấu trúc hình học pha như thế nào. Với tính chất vật liệu cho trước C và tỷ lệ thể tích pha v cũng cho được biểu thức đánh giá. Các đánh giá hẹp hơn đánh giá HS có chứa thêm các thông tin bậc cao về hình học pha của vật liệu đã được xây dựng bởi các tác giả khác nhau, như Milton (1981), Miller (1969), Phạm D. Trong nghiên cứu [2, 29] hệ số dẫn C eff được xác định theo (1.16)    0 , x V Trường khả dĩ phân cực E trong (1.9) được chọn như sau: 1 C0  Ei  Ei 0  p '  ,ij , j i  1, 2.17)  Trong đó p '  p   p  với trường véc tơ phân cực p   p   (x) :  1  2C0     0     p  1   E (1.18) ta có đánh giá trên như sau: C eff  PC  C0   C** (1.19) trong đó C0 là số dương tùy ý, 1    PC (C0 )      C0 (1.21)  1 X     C  C0 C  C0 A  là các hệ số hình học mô tả xấp xỉ bậc ba hình học của vật liệu [30] LUAN VAN CHAT LUONG download : add luanvanchat@agmail.10) với trường khả dĩ phân cực J được chọn là: 1     q '  ,ij qi '  , q '  q  q   Ji  Ji0   C0   j  1  2  C0    1    q  1    C 1   C    C   C   1 1  J0 (1.22)   0   0  với trường phân cực q   q   (x)  Từ (1.