基于近场动力学三维计算均匀化方法预测纤维增强复合材料有效模量的研究

《Composite Structures》：Three-dimensional computational homogenization of cracked composite materials using state-based peridynamics and MPI parallelization

本文针对纤维增强复合材料微观结构与宏观性能的关联问题，提出了一种基于近场动力学（Peridynamics, PD）的三维计算均匀化框架。研究人员通过建立代表性体积单元（RVE），结合周期性边界条件，精确预测了硼-铝和玻璃-环氧树脂复合材料的全刚度张量。该方法有效克服了传统连续介质力学在模拟材料界面和损伤演化时的局限性，为多尺度力学分析和复合材料设计提供了新工具。

　　

复合材料在现代工程中扮演着至关重要的角色，从航空航天结构到生物医学器械，其应用无处不在。然而，这些材料的宏观性能很大程度上取决于其微观结构的复杂相互作用——纤维如何嵌入基体、界面特性如何、以及损伤如何萌生和扩展。传统连续介质力学理论在描述这类不连续和异质材料行为时面临挑战，特别是在处理裂纹扩展和材料失效等问题时。这促使研究人员寻求新的计算框架，能够更自然地捕捉材料从微观到宏观的跨尺度力学行为。

正是在这样的背景下，近场动力学这一非局部连续介质力学理论应运而生。与经典理论不同，近场动力学通过积分方程描述材料点之间的相互作用，避免了基于偏微分方程的传统方法在位移不连续处出现的奇异性问题。这项发表在《Composite Structures》上的研究，成功开发了一套完整的三维近场动力学计算均匀化框架，用于预测纤维增强复合材料的有效弹性性能。

研究人员采用了基于普通状态型近场动力学公式，结合周期性边界条件，对包含单一纤维的代表性体积单元进行均匀化分析。通过施加六种独立的宏观应变状态，计算相应的宏观应力响应，从而提取完整的有效刚度张量。该方法的核心优势在于其能够自然处理材料界面和潜在的不连续问题，为研究复合材料的损伤和失效行为奠定了基础。

关键技术方法包括：基于MPI的并行计算架构实现大规模三维问题求解；周期性边界条件通过拉格朗日乘子法精确施加；材料界面处采用调和平均法处理属性过渡；基于M?ller-Trumbore算法的射线追踪技术用于识别键与裂纹表面的相互作用。样本为硼-铝和玻璃-环氧树脂两种复合材料体系。

研究结果

有效刚度预测的收敛性分析

通过系统研究不同离散化水平（20^3到90^3）下的计算结果，发现当分辨率达到70^3以上时，大多数刚度分量的相对误差降至1%以下。轴向刚度分量C₁₁₁₁在50^3分辨率时误差已低于0.1%，而剪切分量如C₂₃₂₃收敛相对较慢，表明其对纤维-基体界面离散化更为敏感。

硼-铝复合材料的性能评估

对于纤维体积分数为0.47的硼-铝体系，计算得到的有效刚度张量表现出明显的横观各向同性特征。轴向模量C₁₁₁₁为231.1 GPa，横向模量C₂₂₂₂和C₃₃₃₃为164.5 GPa，剪切模量C₂₃₂₃为47.5 GPa，与理论预测吻合良好。耦合项的值比主项小几个数量级，符合物理预期。

数值方法的稳健性验证

采用PETSc库中的KSPBCGSL迭代求解器，在相对容差10^-10下所有加载情况均稳定收敛。MPI并行实现显示出近线性强扩展性，在130万自由度规模仍保持高效计算。

讨论与结论

本研究发展的三维近场动力学计算均匀化框架，成功实现了对纤维增强复合材料有效弹性性能的精确预测。方法的核心价值在于其能够自然处理材料不连续问题，为后续研究损伤和失效行为奠定了坚实基础。

研究结果表明，近场动力学作为一种非局部理论，在复合材料多尺度建模中具有独特优势：其一，积分形式的控制方程避免了传统连续介质力学在材料界面处的导数计算困难；其二，通过键的断裂直接描述损伤演化，物理意义明确；其三，周期性边界条件的引入确保了均匀化过程的数学严谨性。

该方法在70^3至90^3分辨率范围内表现出良好的收敛性，证实了其数值稳健性。硼-铝复合材料的算例验证了框架预测横观各向同性材料有效刚度的能力，各刚度分量的相对误差在合理范围内，满足工程分析精度要求。

值得关注的是，研究中发展的材料界面处理策略——特别是基于调和平均的参数分配方法，有效解决了异质材料界面处的力学行为描述问题。而射线追踪技术的引入，为后续研究含裂纹复合材料的性能退化提供了技术准备。

这项工作为复合材料的多尺度分析提供了新范式，将传统均匀化方法的应用范围扩展至包含损伤和失效的更一般情形。未来可进一步研究循环加载下的疲劳损伤、动态冲击响应等更复杂的力学问题，推动近场动力学在复合材料工程中的应用深度。

从更广阔的视角看，这种基于非局部理论的均匀化框架有望应用于更具挑战性的材料体系，如功能梯度材料、多孔材料等，为先进材料设计提供可靠的数值工具。随着计算能力的不断提升和算法的进一步优化，近场动力学有望成为复合材料多尺度模拟的重要支柱方法。