开孔材料的断裂力学：基于Voronoi的方法

《Mechanics of Materials》：Fracture mechanics of open-porous materials: A Voronoi-based approach

　　

开孔材料的断裂行为由微观结构构型与裂纹扩展机制之间的复杂相互作用所主导，因此有必要进行系统研究以确保其在结构、生物医学和热绝缘应用中的可靠使用。量化断裂韧性（fracture toughness, K_IC）和应变能释放率（strain energy release rate, G）对于优化应力下的机械稳定性至关重要。研究人员提出了一种基于Laguerre-Voronoi剖分（Laguerre-Voronoi tessellation, LVT）的计算框架，用于模拟开孔固体的复杂细胞结构。在拉伸载荷下，该模型再现了Maiti、Gibson和Ashby提出的经典标度关系。有效弹性模量（effective elastic modulus）表现出对固体分数（solid fraction, φ_s）的近似二次依赖关系，而归一化断裂行为遵循固体分数的平方根函数。模型对归一化断裂韧性（normalized fracture toughness, K_IC）的预测与文献报道的实验数据吻合良好。这种一致性证实了基于Voronoi的微观结构能够捕捉已建立的幂律，同时也揭示了与显著异质性相关的偏差。随后，研究人员系统地改变了关键微观结构参数，包括固体分数、孔隙尺寸和孔壁厚度，以及宏观裂纹几何特征（如裂纹深度-高度比和裂纹方向），以评估它们对断裂阻力的影响。结果表明，裂纹扩展和能量耗散对微观结构特征具有关键依赖性，展示了如何通过调整孔隙率（porosity）和孔隙尺寸分布（pore size distribution, PSD）来定制机械响应。这些发现有助于更深入地理解多孔系统中的断裂力学，并强调了基于微观结构的材料设计以增强机械性能的必要性。

开孔材料（open-porous materials）因其高孔隙率和连通孔隙网络，在结构、生物医学和热绝缘等领域具有广泛的应用前景，但其断裂行为受微观结构构型与裂纹扩展机制之间复杂相互作用的支配，亟需系统研究以优化机械性能。已有经典标度关系（如Maiti、Gibson和Ashby提出的）主要针对均匀孔结构，但实际材料常呈现高度异质性，导致断裂行为偏离理论预测。本研究旨在通过基于Laguerre-Voronoi剖分（Laguerre-Voronoi tessellation, LVT）的计算框架，系统模拟开孔固体的复杂细胞结构，揭示微观结构参数（如固体分数φ_s、孔隙尺寸分布PSD、孔壁几何）和宏观裂纹几何（如裂纹深度-高度比a/h和裂纹角度α）对断裂韧性和应变能释放率的影响，从而为基于微观结构的材料设计提供理论指导。该研究发表在《Mechanics of Materials》。

研究人员采用随机密堆积多分散球体（RCPPS）结合Laguerre-Voronoi剖分（LVT）生成具有可控孔隙尺寸分布（PSD）和固体分数（φ_s）的开孔微观结构。有限元分析（FE）使用Hughes-Liu共旋梁单元离散孔壁，基于线性弹性本构模型，以von Mises等效应力作为断裂准则，通过单元侵蚀模拟脆性断裂。应用线性弹性断裂力学（LEFM）计算应力强度因子K_IC和应变能释放率G。材料属性取自Polymethacrylimide泡沫（Maiti等，1984）。

**3.1. 孔隙尺寸分布（PSD）的影响**
通过改变正态分布的标准差（s）和均值（μ）研究PSD的影响。**3.1.1. s的变化**：不同s下，孔壁长径比（d/l）与φ_s保持二次标度关系，弹性模量（E/Es）与φ_s的幂律关系不变（指数n≈1.8）。归一化断裂韧性（K_IC/σ_fs√(πl)）和G随s增大而降低，但G随s增大而升高，因异质性导致应力集中和能量耗散增强。**3.1.2. μ的变化**：增大μ导致孔壁更长，等效缺陷尺寸增大，断裂应力σ_f降低，但G升高。实验数据（PMI-E和PUR泡沫）与模型预测吻合，Maiti模型（c=0.65, m=1.5）适用于小μ或小s的均匀分布，Green模型（c=0.28, m=1.3）适用于大μ或大s的宽分布。

**3.2. 孔壁长度（l）和直径（d）的影响**
在恒定直径配置（固定d，变化l）和恒定长度配置（固定l，变化d）下，E/Es与φ_s均呈幂律关系。恒定长度配置中K_IC对φ_s的幂律指数更高，更符合经典理论，表明通过调整d改变φ_s更准确地反映断裂行为。较长和较厚的孔壁储存更多弹性应变能，导致更高的K_IC和G。

**3.3. 裂纹几何的影响**
**3.3.1. 裂纹深度-高度比（a/h）**：随a/h增大，有效承载孔壁减少，刚度和断裂强度降低。对于高孔隙率（φ_s=0.01），K_IC独立于a/h；对于低孔隙率，K_IC在a/h≤0.15时随a/h增大而增大，之后趋于稳定。浅裂纹（a/h=0.05）导致曲折路径，增加能量耗散；深裂纹（a/h>0.2）呈现直线脆性断裂。**3.3.2. 裂纹角度（α）**：断裂应力σ_f与α无关，仅随a/h增大而减小。当α≤60°时，K_IC和G基本恒定，保持纯I型断裂；当α>90°时，K_IC降低而G升高，表明混合模式断裂（含II型）贡献增强。所有角度下裂纹路径均呈现曲折。

**结论**：本研究系统揭示了微观结构特征（PSD、φ_s、孔壁几何）和裂纹几何对开孔材料断裂行为的影响。有效弹性模量满足幂律关系E/Es = C₁φ_sⁿ，n=1.80±0.06；归一化断裂韧性满足K_IC/σ_fs√(πl) ∝ φ_s^m，m=1.40±0.08；应变能释放率G ∝ φ_s^m1，m1=1.00±0.15。PSD的异质性显著影响标度指数，均匀分布更接近经典理论（m=1.48）。模型无需参数拟合，可直接使用实验测量的微观结构和力学性能预测宏观响应，为设计高性能开孔材料提供了有力工具。未来研究可引入非高斯PSD、各向异性及非线性孔壁行为，以扩展框架适用性。