单轴加载下脆性TPMS裂纹扩展预测：基于测地线和截面分析的双概念方法

《International Journal of Ceramic Engineering & Science》：Predicting Crack Propagation in Brittle TPMS Under Uniaxial Loading: A Dual Concept Approach via Geodesic Lines and Cross-Sectional Analysis

　　

三周期极小曲面（TPMS）结构因其优越的力学性能、几何连续性和适用于增材制造而在工程应用中迅速被采用。本研究探讨了四种代表性TPMS单胞（Primitive、Gyroid、Neovius和IWP）的独特形态对脆性裂纹扩展的影响。提出了两个几何驱动概念：（1）裂纹沿表面测地线（geodesic lines）扩展的预测，以及（2）通过沿加载方向的截面面积分布评估损伤局部化。研究人员开发了一个数学框架来计算曲面上的解析测地线，并将其与数值获得的裂纹路径进行比较。此外，还分析了每个TPMS单胞中不规则截面形状的面积，以评估截面变化如何影响抗压性能。TPMS结构的基体材料是氧化铝陶瓷，在有限元（FE）分析中使用Johnson–Holmquist II本构模型进行建模。结果表明，数值裂纹扩展路径与解析测地线高度一致，且损伤局部化由小截面面积区域控制。该方法为预测复杂晶格拓扑中的断裂提供了稳健的几何基础。

论文解读文章

**研究背景与问题**
人工多孔结构的发展旨在模拟自然系统中固有的优化拓扑，其中孔隙分布决定力学性能。现有三类多孔结构：随机泡沫（孔隙率控制困难）、点阵结构（节点处应力集中）和三周期极小曲面（TPMS）结构。TPMS由数学函数精确定义，具有零平均曲率、高连通性和均匀应力分布特性，且自支撑特性使其适合增材制造，已在航空航天、汽车吸能和生物医学植入物中广泛应用。然而，尽管已有大量研究强调TPMS的力学优势，但针对其独特拓扑如何影响损伤行为（尤其是脆性材料中的裂纹扩展）的研究仍有限。本研究基于前期对表面倾斜角影响的研究，提出两个几何驱动概念：（1）裂纹沿表面测地线（geodesic lines）扩展，遵循最小能量消耗原理；（2）通过沿加载方向截面面积分布评估损伤局部化，依据经典力学中最小截面处应力最大导致裂纹萌生的规律。研究人员旨在回答：结构几何如何影响单轴加载下的裂纹扩展？具体分为两个子问题：裂纹是否沿最小能量路径扩展？截面面积变化如何影响损伤局部化？论文发表在《International Journal of Ceramic Engineering》。

**研究内容与意义**
研究人员选取四种典型TPMS单胞（Primitive、Gyroid、Neovius、IWP），以氧化铝陶瓷为基体材料，采用Johnson–Holmquist II（JH-2）本构模型进行有限元（FE）分析，模拟单轴压缩和拉伸下的脆性断裂。通过求解解析测地线方程（作为初值问题）预测裂纹路径，并计算截面面积沿加载方向的分布来识别易破坏区域。结果表明，数值裂纹路径与解析测地线高度一致，且损伤局部化由小截面面积区域主导。该框架揭示了复杂拓扑与断裂行为的关联，为脆性晶格结构的抗断裂设计提供了几何基础。

**关键技术方法**
（1）TPMS几何生成：利用数学公式（Primitive、Gyroid、Neovius、IWP）在MATLAB中生成曲面，经MeshLab重新网格化后导入FreeCAD生成实体，最后在ABAQUS/Explicit中用四节点四面体单元（C3D4，约20万单元）进行FE建模。（2）有限元模拟：建立单轴压缩模型（顶部施加1%应变，底部完全约束，周期性边界条件保持表面平面性）和单轴拉伸预裂模型（用于测地线分析）。采用JH-2本构模型模拟氧化铝陶瓷的脆性行为，通过删除达到临界损伤阈值的单元实现裂纹扩展可视化。（3）测地线计算：将裂纹路径预测转化为求解测地线方程的初值问题，初始点由预裂纹缺口确定，初始方向由FE裂纹点或基于应力强度因子（K_I、K_II）的裂纹偏转角确定。通过计算路径的测地曲率（接近零表示路径为测地线）评估裂纹路径与解析测地线的一致性。（4）截面面积分析：将单胞沿加载方向分割为64个平行截面，通过数值算法（最近邻点排序、对称性利用、角点添加）提取每个截面的不规则形状并计算面积，进而分析截面面积分布与损伤局部化的相关性。

**研究结果**
**4.1 不同TPMS的损伤模拟结果**
通过JH-2本构模型模拟四种TPMS单胞（体积分数均为22%）的单轴压缩，得到应力-应变曲线。结果显示抗压强度排序为：IWP最高，其次为Neovius、Gyroid，Primitive最低。裂纹萌生最早发生在Primitive和Gyroid，随后是IWP和Neovius。Primitive曲线呈平台状而非突然应力下降，表明其能承受更长的压力。应力-应变曲线呈锯齿状，由间歇性局部开裂和载荷重分布引起。在0.3%、0.5%、0.7%压缩应变下的损伤发展表明，Primitive在0.3%应变时已达到最大压缩应力，损伤更显著；至0.7%应变时，各结构裂纹扩展特征不同，最终失去承载能力。累积能量吸收（0.8%应变）随壁厚增加而增大，且同样遵循IWP > Neovius > Gyroid > Primitive的顺序。

**4.2 截面面积分析**
通过将单胞沿高度方向分割为64个截面并数值积分计算体积分数，结果略低于理论值（误差可通过增加截面数减小）。各TPMS的截面面积分布曲线呈对称模式：Primitive的极值点最少，变化平缓；Gyroid的峰谷数量与Neovius、IWP相同但过渡更平滑；Neovius和IWP的截面面积在最大值和最小值之间变化剧烈。各结构的相对最大值和最小值变化较小，表明最大和最小截面形状不同但面积相近。频率分布显示仅Neovius和IWP的截面面积超过0.3（归一化值），与其高抗压强度一致。Gyroid的最大值较低、最小值较高，导致其应力-应变曲线居中。Primitive截面曲线的微弱振荡可能是其曲线呈平台状的原因。将截面面积分布与von Mises应力分布和最终损伤状态对比，发现高应力区集中在小截面区域，且截面面积突变（如IWP）会加剧应力集中，导致局部裂纹形成。

**4.3 利用测地线预测裂纹路径**
**4.3.1 Gyroid结构**
建立带预裂纹缺口的Gyroid拉伸模型，分别采用粗网格（111551单元）和细网格（5475150单元）进行损伤模拟。提取裂纹节点坐标并拟合三阶多项式曲线，计算其测地曲率。结果发现细网格裂纹路径与解析测地线高度吻合，测地曲率接近零；粗网格路径偏离较大。三种路径的测地曲率均呈初始下降趋势，表明裂纹在早期传播中趋于自校正，逐步接近测地线轨迹。

**4.3.2 IWP结构**
对IWP（体积分数50%）带预裂纹缺口模型，基于应力强度因子K_I、K_II计算裂纹偏转角（通过Franc3D获得K_I=0.0214 MPa·m^1/2，K_II=0.087 MPa·m^1/2，偏转角θ=24.6°），建立解析测地线初值问题。结果解析测地线与有限元裂纹路径高度一致，测地曲率均接近零，证实裂纹倾向于沿测地线扩展。与已有PLA基IWP研究（裂纹沿垂直于加载方向的平面扩展，且与高应变集中区域一致）对比，本研究结果与之相似。

**讨论与结论**
讨论部分指出，裂纹路径与解析测地线的吻合验证了最小能量消耗原理在脆性TPMS中的适用性；截面面积分析揭示了损伤局部化由小截面区域和面积突变主导，这解释了不同TPMS的抗压强度差异。该框架为理解复杂拓扑与断裂行为的关联提供了稳健的几何基础。
**研究结论**：本研究引入了两个基于拓扑的概念来预测脆性断裂中的裂纹扩展。第一种方法通过求解解析测地线方程作为初值问题来预测裂纹路径。初始裂纹方向通过有限元损伤模拟中的第二个裂纹点或基于应力强度因子（stress intensity factors, SIFs）的裂纹偏转角来确定。裂纹偏转准则独立于复杂的损伤模型，并提供了与有限元结果高度一致的裂纹轨迹。第二种方法通过沿压缩加载轴分析截面面积变化来识别易破坏区域。小截面面积和截面面积的突变被证明控制损伤局部化，并显著影响不同TPMS设计之间的相对抗压强度。所提出的框架有助于理解复杂拓扑与断裂行为之间的相关性。