Phương pháp bao hàm tương đương giải quyết bài toán cơ học vi mô nào?

Phương pháp này giải quyết bài toán xác định trường ứng suất và biến dạng cục bộ xung quanh các hạt dị tật đàn hồi khác pha, từ đó cho phép ước lượng các mô đun cơ học vĩ mô hiệu dụng của vật liệu composite.

Điểm cốt lõi trong phát kiến của Eshelby năm 1957 là gì?

Eshelby đã chứng minh rằng nếu một thể bao hàm hình elipsoid nhúng trong môi trường đàn hồi vô hạn chịu một biến dạng tự do đồng nhất, thì trường biến dạng cưỡng bức sinh ra bên trong thể bao hàm đó cũng hoàn toàn đồng nhất.

Mô hình Mori-Tanaka khác gì so với mô hình pha loãng cổ điển?

Mô hình pha loãng giả định các hạt không tương tác trong trường ứng suất từ xa vô hạn, trong khi mô hình Mori-Tanaka nhúng thể bao hàm vào trường ứng suất trung bình hữu hiệu của pha nền để tính đến ảnh hưởng nồng độ hạt.

Phương pháp bao hàm tương đương: Nguyên lý và ứng dụng

Q: Tại sao hình học elipsoid lại giữ vai trò trung tâm trong phương pháp này?

Theo giả thuyết Eshelby, elipsoid là dạng hình học duy nhất bảo đảm trường biến dạng cưỡng bức đồng nhất bên trong hạt đối với mọi biến dạng tự do tùy ý, cho phép thiết lập nghiệm giải tích tường minh.

Phương pháp bao hàm tương đương là một kỹ thuật giải tích nền tảng trong cơ học vi mô, được sử dụng để xác định trường ứng suất và biến dạng cục bộ cũng như dự đoán các đặc tính cơ học vĩ mô hiệu dụng của vật liệu composite và môi trường không đồng nhất, thông qua việc thay thế một thể dị tật hoặc thể hạt gia cường bằng một thể bao hàm giả định có cùng các hằng số đàn hồi với môi trường nền nhưng mang một trường biến dạng tự do tương đương thích hợp.

Bản chất vật lý và cơ sở cơ học vi mô

Trong cơ học vật liệu, một thể bao hàm được định nghĩa là một miền phụ hữu hạn bên trong môi trường vô hạn có các hằng số đàn hồi giống hệt môi trường nền xung quanh, nhưng trải qua một biến dạng tự do do các hiệu ứng vật lý như chuyển pha tinh thể, giãn nở nhiệt hoặc từ biến. Ngược lại, một dị tật hay thể không đồng nhất là một miền phụ có các hằng số đàn hồi hoàn toàn khác biệt so với vật liệu nền, chẳng hạn như hạt gia cường gốm hoặc lỗ xốp vi mô phân bố trong nền kim loại hoặc polymer.

Khi vật liệu chịu tác động của tải trọng cơ học từ xa, sự hiện diện của các thể dị tật sẽ làm nhiễu loạn trường ứng suất và biến dạng đồng nhất trong môi trường. Việc giải trực tiếp bài toán truyền ứng suất qua các mặt phân cách có tính chất cơ học gián đoạn là thách thức toán học phức tạp. Phương pháp bao hàm tương đương của Eshelby giải quyết bài toán này bằng cách đưa ra giả định tương đương cơ học: trường ứng suất sinh ra bên trong và xung quanh thể dị tật thực tế có thể được tái lập chính xác bằng một thể bao hàm tưởng tượng có cùng mô đun đàn hồi với nền, miễn là thể bao hàm này được gán một trường biến dạng tự do tương đương phù hợp.

Lịch sử phát triển và các công trình nền tảng

Năm 1957, nhà vật lý và cơ học John Douglas Eshelby đã công bố công trình trên tạp chí Proceedings of the Royal Society of London Series A, đặt nền móng lý thuyết cho phương pháp bao hàm tương đương. Eshelby đã phát hiện một tính chất toán học quan trọng: nếu một thể bao hàm dạng elipsoid nhúng trong môi trường đàn hồi vô hạn đồng nhất chịu một biến dạng tự do đồng nhất, thì biến dạng cưỡng bức sinh ra bên trong thể bao hàm đó cũng hoàn toàn đồng nhất. Đến năm 1959, Eshelby tiếp tục mở rộng công trình nghiên cứu trên cùng tạp chí nhằm thiết lập các nghiệm giải tích tường minh cho trường đàn hồi bên ngoài thể bao hàm elipsoid.

Dựa trên nghiệm giải tích của Eshelby cho một hạt đơn lẻ trong môi trường vô hạn, các nhà khoa học đã phát triển nhiều mô hình trường trung bình để tính toán tính chất vĩ mô cho các hệ vật liệu chứa nhiều hạt phân bố. Năm 1973, Mori và Tanaka đã công bố trên tạp chí Acta Metallurgica mô hình trường trung bình hữu hiệu. Mô hình Mori-Tanaka sử dụng phương pháp bao hàm tương đương của Eshelby nhưng đặt thể bao hàm vào một trường ứng suất trung bình hữu hiệu của vật liệu nền thay vì trường đặt từ xa, từ đó tính đến tương tác tương đối giữa các hạt khi nồng độ pha gia cường tăng lên. Đến năm 1987, Yakov Benveniste công bố trên tạp chí Mechanics of Materials công trình hệ thống hóa lý thuyết Mori-Tanaka thông qua ten-xơ tập trung ứng biến trực tiếp, củng cố tính nhất quán nhiệt động học của phương pháp đối với composite hai pha.

Nguyên lý toán học của phương pháp bao hàm tương đương

Xét một môi trường đàn hồi vô hạn đồng nhất có ten-xơ độ cứng đàn hồi $C_{ijkl}$ , chứa một thể dị tật elipsoid có ten-xơ độ cứng đàn hồi $C^*_{ijkl}$ . Dưới tác dụng của trường biến dạng đều đặt từ xa $\varepsilon^0_{ij}$ , trường ứng suất thực tế bên trong thể dị tật được xác định theo định luật Hooke tổng quát:

\sigma_{ij} = C^*_{ijkl} (\varepsilon^0_{kl} + \varepsilon_{kl})

Trong phương trình trên, đại lượng $\varepsilon_{kl}$ là trường biến dạng nhiễu loạn sinh ra do sự hiện diện của thể dị tật.

Phương pháp bao hàm tương đương thay thế thể dị tật trên bằng một thể bao hàm có cùng ten-xơ độ cứng đàn hồi $C_{ijkl}$ của vật liệu nền, đồng thời mang một trường biến dạng tự do tương đương giả định $\varepsilon^{**}_{kl}$ . Trường ứng suất tương đương trong thể bao hàm được biểu diễn theo công thức:

\sigma_{ij} = C_{ijkl} (\varepsilon^0_{kl} + \varepsilon_{kl} - \varepsilon^{**}_{kl})

Theo định lý Eshelby, đối với thể bao hàm dạng elipsoid nhúng trong môi trường đàn hồi vô hạn, biến dạng nhiễu loạn bên trong thể bao hàm liên hệ tuyến tính với biến dạng tự do thông qua ten-xơ Eshelby bậc bốn $S_{klmn}$ :

\varepsilon_{kl} = S_{klmn} \varepsilon^{**}_{mn}

Đồng nhất hai biểu thức ứng suất thực tế và ứng suất tương đương, ta thu được phương trình điều kiện tương đương Eshelby:

C^*_{ijkl} (\varepsilon^0_{kl} + S_{klmn} \varepsilon^{**}_{mn}) = C_{ijkl} (\varepsilon^0_{kl} + S_{klmn} \varepsilon^{**}_{mn} - \varepsilon^{**}_{kl})

Từ hệ phương trình đại số tuyến tính này, ta giải ra trường biến dạng tự do tương đương $\varepsilon^{**}_{kl}$ theo biến dạng đặt từ xa $\varepsilon^0_{kl}$ và các đặc trưng đàn hồi của nền và hạt:

\left[ C^*_{ijmn} S_{mnkl} + C_{ijmn} (I_{mnkl} - S_{mnkl}) \right] \varepsilon^{**}_{kl} = (C_{ijkl} - C^*_{ijkl}) \varepsilon^0_{kl}

Trong biểu thức trên, $I_{mnkl}$ là ten-xơ đơn vị bậc bốn đối xứng. Khi đã xác định được biến dạng tự do tương đương, toàn bộ trường ứng suất và biến dạng bên trong cũng như bên ngoài hạt đều được giải quyết hoàn toàn bằng giải tích.

Ten-xơ Eshelby và các dạng hình học đặc biệt

Ten-xơ Eshelby là một ten-xơ bậc bốn có tính đối xứng phụ đối với hai chỉ số đầu và hai chỉ số cuối ( $S_{ijkl} = S_{jikl} = S_{ijlk}$ ), nhưng trong trường hợp tổng quát không thỏa mãn tính đối xứng chính ( $S_{ijkl} \neq S_{klij}$ ). Thành phần của ten-xơ Eshelby chỉ phụ thuộc vào tỷ lệ hình học giữa các bán trục của hình elipsoid và hệ số Poisson của vật liệu nền, hoàn toàn không phụ thuộc vào kích thước tuyệt đối hay độ cứng của pha gia cường. Nhờ tính chất này, nghiệm Eshelby có thể áp dụng cho các cấu trúc kích thước đa dạng từ quy mô nanomet đến quy mô công trình vĩ mô.

Nghiệm elipsoid tổng quát cho phép suy biến thành các dạng hình học vi mô đặc biệt:

Hạt hình cầu: Các bán trục elipsoid bằng nhau, các thành phần của ten-xơ Eshelby chỉ phụ thuộc duy nhất vào hệ số Poisson của vật liệu nền, thường dùng mô hình hóa hạt độn đẳng hướng hoặc lỗ xốp khí hình cầu.
Sợi trụ tròn dài vô hạn: Một bán trục elipsoid tiến tới vô cùng lớn so với hai bán trục còn lại, mô hình hóa ứng xử cơ học của vật liệu composite gia cường bằng sợi liên tục.
Đĩa mỏng hoặc vảy dẹp: Một bán trục elipsoid tiến tới rất nhỏ so với hai bán trục phẳng, mô hình hóa các hạt độn dạng phiến như graphene, đất sét nano hoặc vết nứt vi mô phẳng trong vật liệu giòn.

Mở rộng sang các mô hình đồng hóa trường trung bình

Phương pháp bao hàm tương đương cho một hạt đơn lẻ đóng vai trò viên gạch nền tảng để xây dựng các mô hình đồng hóa vi mô cho vật liệu đa pha:

Mô hình pha loãng: Bỏ qua hoàn toàn tương tác giữa các hạt, giả định mỗi hạt cư xử độc lập trong trường ứng suất từ xa của vật liệu nền.
Mô hình Mori-Tanaka: Giả định mỗi thể bao hàm tương đương được nhúng vào một môi trường có trường ứng suất bằng ứng suất trung bình của pha nền, dự đoán chính xác mô đun hiệu dụng cho composite có pha nền liên tục.
Mô hình tự nhất quán: Giả định từng hạt gia cường được nhúng trực tiếp vào môi trường đồng nhất tương đương có các mô đun cơ học vĩ mô chưa biết, dẫn đến hệ phương trình phi tuyến cần giải lặp.

Bảng dưới đây tóm tắt sự so sánh giữa các mô hình vi cơ học dựa trên nền tảng Eshelby:

Mô hình lý thuyết	Trường nền tham chiếu	Phạm vi nồng độ hạt	Ưu điểm và đặc điểm nổi bật
Mô hình pha loãng	Trường đặt từ xa vô hạn	Nồng độ rất thấp	Công thức giải tích đơn giản, nghiệm chính xác tiệm cận khi nồng độ tiến tới không
Mô hình Mori-Tanaka	Ứng suất trung bình trong pha nền	Nồng độ trung bình khi pha nền còn liên tục	Dạng giải tích tường minh, tính toán thuận tiện, phù hợp cho hệ composite hai pha có pha nền liên tục chiếm ưu thế
Mô hình tự nhất quán	Môi trường đồng nhất vĩ mô hiệu dụng	Nồng độ cao hoặc đa tinh thể	Đối xử bình đẳng giữa các pha vật liệu, phù hợp mô hình hóa đa tinh thể hoặc vật liệu không có pha nền rõ rệt

Nghiên cứu ứng dụng tại Việt Nam và trên thế giới

Tại Việt Nam, các phương pháp cơ học vi mô dựa trên nguyên lý Eshelby và mô hình Mori-Tanaka được nghiên cứu và ứng dụng sâu rộng trong tính toán vật liệu tiên tiến. Điển hình là công trình nghiên cứu của Trần Trung Dũng, Trần Thanh Tùng, Phạm Đức Chính và Nguyễn Đình Đức (2021) công bố trên tạp chí Journal of Mechanical Science and Technology. Các tác giả đã phân tích và so sánh phương pháp tiếp cận xấp xỉ phân cực với xấp xỉ Mori-Tanaka để xác định các đặc tính đàn hồi vĩ mô hiệu dụng của vật liệu đàn hồi đẳng hướng nhiều thành phần chứa các thể bao hàm hình cầu hoặc định hướng ngẫu nhiên, làm sáng tỏ phạm vi ứng dụng và độ chính xác của các sơ đồ trường trung bình.

Trên thế giới, phương pháp bao hàm tương đương Eshelby là công cụ tiêu chuẩn tích hợp trong các phần mềm mô phỏng đa tỉ lệ hiện đại như Digimat hoặc ABAQUS. Kỹ thuật này hỗ trợ các kỹ sư dự đoán chính xác ứng xử nhiệt đàn hồi, độ từ biến và sự tích tụ hư hại của composite hàng không vũ trụ và vật liệu ô tô thế hệ mới.

Giới hạn lý thuyết và ranh giới phương pháp luận

Dù có giá trị lý thuyết sâu sắc, phương pháp bao hàm tương đương Eshelby vẫn tồn tại một số giới hạn phương pháp luận quan trọng:

Giới hạn hình học elipsoid: Theo giả thuyết Eshelby, trường biến dạng cưỡng bức đồng nhất bên trong hạt đối với mọi biến dạng tự do tùy ý chỉ đúng cho hình học elipsoid. Đối với các thể bao hàm có góc nhọn như hình lập phương, hình đa diện hoặc hạt có cạnh sắc nhọn, trường biến dạng bên trong hạt không còn đồng nhất mà tập trung ứng suất mạnh tại các đỉnh nhọn, đòi hỏi phải sử dụng phương pháp phần tử hữu hạn hoặc các kỹ thuật giải tích số.

Giới hạn giả thiết biến dạng bé và đàn hồi tuyến tính: Lý thuyết Eshelby cổ điển xây dựng trên giả thiết đàn hồi tuyến tính và biến dạng bé. Khi vật liệu nền hoặc hạt gia cường chuyển sang trạng thái biến dạng dẻo, phá hủy phi tuyến hoặc chịu biến dạng hữu hạn lớn, phương pháp bao hàm tương đương cần được mở rộng phức tạp thông qua các mô hình tiếp tuyến hoặc mô hình biến phân nâng cao.