双连续复合材料：线弹性力学中的通用均匀化(Homogenization)方法

《Mechanics of Materials》：Bicontinuous composite materials: general homogenization methods in linear elasticity

　　

提出新的解析与数值线弹性均匀化（Homogenization）方法，用于具有完全互穿（Interpenetrated）微结构的双连续（Bicontinuous）复合材料。在推广插值双夹杂模型（Interpolative Double Inclusion Model, IDI）的基础上，新的均值场（Mean-field）模型考虑两种极端构型：第一种对应于材料2的低浓度夹杂（特定形状，如扁球Oblate Spheroids）嵌入材料1基体相中；第二种对应于材料1的低浓度夹杂（另一形状，如连续纤维Continuous Fibers）嵌入材料2基体相中。对于中间浓度，通过对两种极端构型偏应变集中张量（Partial Strain Concentration Tensor）之逆进行插值，实现两极端间的连续过渡，并提出一种新的正弦幂函数（Power of a Sine）插值函数，含有一个浓度参数使两极端构型权重相等。随后将两种取向平均（Orientation Averaging）策略应用于单向模型：或在刚度（Stiffness）上进行两步均匀化，或在偏应变集中张量上进行平均。提出第二种基于数值切片（Numerical Slicing）的均匀化方法：沿给定方向对三维微结构切片并均匀化各切片，再将由全粘结切片组成的层合板（Laminate）用全场有限元（Finite Element, FE）进行均匀化。将所提两种方法针对文献数据集及各向同性互穿相的直接有限元参考模拟进行验证，预测的有效各向同性弹性模量（如杨氏模量Young's Modulus与剪切模量Shear Modulus）与参考数据吻合良好，证明了新均值场模型的准确性，以及当仅有少量二维数字图像时可应用切片法（Slicing-based Method）的价值。

论文解读：《Bicontinuous composite materials: general homogenization methods in linear elasticity》

一、研究背景与意义

传统复合材料通常由离散纤维或颗粒分散于连续基体中构成，而双连续（Bicontinuous）复合材料的两相均形成三维全连通网络，任一组分移除后剩余相仍可自支撑（Open-celled Structure）。此类材料（如金属-陶瓷互穿体系）可协同发挥各组分的本征性能优势，表现出优于常规复合材料的力学与物理性能。然而，由于互穿微结构的复杂性及实验数据匮乏，传统基于Eshelby解的微观力学模型（如Mori-Tanaka方案、自洽法Self-Consistent Scheme、Maxwell模型等）主要针对"夹杂-基体"拓扑设计，难以直接适用于无孤立相的双连续体系。已有的互穿结构模型（如Tuchinskii的三光纤模型、Schmauder等的自洽Matricity模型、Franciosi等的IPS-IFS方法）或局限于特定周期性微结构，或依赖经验公式拟合体积分数，缺乏普适且无需经验参数的预测手段。为此，研究人员以Lielens插值双夹杂模型（IDI）和Franciosi等的IPS-IFS方法为基础，发展了两套适用于一般双连续复合材料、无经验拟合参数的新型均匀化方法，并通过文献数据与三维全场的有限元（Finite Element, FE）模拟予以验证。该文发表于《Mechanics of Materials》。

二、主要关键技术方法概述

研究人员基于Lielens插值双夹杂模型（IDI）与Franciosi等的各向同性片层/纤维系统（IPS-IFS）思想，发展了新的均值场均匀化模型（New IDI Model）：定义两相互为基体和夹杂的两个极端构型，采用新型正弦幂插值函数在偏应变集中张量逆空间插值过渡，并对单向模型实施两种取向平均策略（刚度空间两步均匀化或偏应变集中张量平均）。同时提出基于数值切片（Numerical Slicing）的全场均匀化方法：沿正交方向切取三维代表性体积单元（Representative Volume Element, RVE）的二维切片分别均匀化，再将其视为完全粘结层合板用有限元均匀化，经多方向平均获得三维等效性能；当仅有二维显微图像时亦可近似估算。模型预测结果与文献实验/数值数据及直接三维有限元均匀化结果对比校验。

三、研究结果

Background micromechanical models（背景微观力学模型）

研究人员综述了基于Eshelby等效夹杂理论的经典微观力学均匀化框架（稀浓度近似、Mori-Tanaka方案、自洽法、增量法及Maxwell模型）及Lielens插值双夹杂模型（IDI），指出这些模型本质上假设一相为连续基体包裹另一相夹杂，无法直接描述双连续互穿微结构中两相均连续的拓扑特征，为后续模型推广提供理论出发点。

A new mean-field homogenization model for bicontinuous composites（双连续复合材料的新均值场均匀化模型）

研究人员推广IDI模型，将双连续复合材料视为两相互为"基体-夹杂"的两个极端构型：极端Ⅰ—低体积分数材料2（取扁球体，纵横比α₁）嵌入材料1基体；极端Ⅱ—低体积分数材料1（取连续长纤维，纵横比α₂→∞）嵌入材料2基体。中间体积分数f₂通过插值函数w(f₂)=sinⁿ(πf₂/2)在两极端偏应变集中张量A⁽¹⁾与A⁽²⁾的逆间过渡，即[A(f₂)]^-1=w(f₂)[A⁽¹⁾]^-1+[1-w(f₂)][A⁽²⁾]^-1，n为控制对称性的幂指数（n=1时两极端等权位于f₂=0.5）。对单向模型进一步施加两种取向平均——在刚度空间对各取向均匀化结果二次均匀化，或直接对各取向偏应变集中张量平均后再求有效刚度。验证表明该均值场模型对各向同性双连续复合材预测的有效杨氏模量E^*与剪切模量G^*与三维有限元参考解及实验数据偏差较小，尤其在中高含量区间优于传统Mori-Tanaka与逆Mori-Tanaka的简单混合。

A new slicing-based full-field homogenization method for bicontinuous composites（基于切片的全场均匀化新方法）

研究人员提出数值切片法：沿x、y、z三正交方向分别切取RVE的系列平行切片（Slice），对每个切片进行二维全场有限元均匀化得到面内/面外等效刚度，再将同方向切片视作串联（Series）或并联（Parallel）粘结层合板二次均匀化得到该方向宏观刚度，最后对三方向结果平均获得三维等效性能。若实际仅能获得若干二维数字图像（如显微断面），可以其统计特征生成代表切片进行近似评估。算例显示，当切片数足够且RVE统计均匀时，该方法预测的各向同性互穿双相复合材料有效模量与直接三维有限元均匀化高度吻合；即便仅有代表性二维切片时亦能提供工程可用的模量估计。

Numerical simulations and assessment of the proposed homogenization methods（数值模拟与新方法评估）

研究人员选取文献中三种典型双连续微结构（包括金属-陶瓷互穿体系及聚合物基双连续相）分别用新IDI均值场模型、切片有限元法和直接三维有限元均匀化（Reference FE）及可用实验数据对比。结果表明：新均值场模型预测的有效杨氏模量E^*与剪切模量G^*在全浓度范围内与参考有限元及实验值吻合良好（相对误差通常在5%以内），显著优于经典单相基体模型的外推结果；切片法预测值与三维参考有限元结果也具较好一致性，且对仅有二维图像的工况具实用价值。验证了两种方法对一般双连续复合材料预测的可靠性与适用性。

四、讨论与结论翻译

研究中提出了两种受IDI模型与IPS-IFS方法启发的双连续复合材料新均匀化方法并进行了评估。

第3节发展的"新IDI（New IDI）"模型基于两种极端构型：第一种对应于材料2极低浓度、具纵横比α₁的球扁体嵌入材料1基体；第二种对应于材料1极低浓度、具纵横比α₂（连续纤维）嵌入材料2基体……（中间浓度通过正弦幂加权偏应变集中张量逆插值连续过渡）。结果表明，该均值场模型可准确预测各向同性双连续复合材料的有效弹性模量。基于微结构切片再层合均匀化的数值方法在仅有二维显微图像时仍能给出合理的三维等效性能估计。两种提出的无经验拟合均匀化方法经文献数据与直接有限元模拟验证，对双连续复合材料力学行为的预测具良好精度与实用性。