力学数值仿真：碳纤维复合材料层合板的失效分析

力学数值仿真在复合材料设计中至关重要，但层间界面建模困难、失效准则选择争议大、网格敏感性分析耗时，如何建立可靠仿真模型预测材料极限承载能力，是工程计算中的核心问题。

去年给一家航空航天供应商做了碳纤维/环氧树脂层合板的力学数值仿真项目。需求很明确：预测[0/45/-45/90]ₛ铺层在面内拉伸载荷下的渐进损伤过程和极限强度。力学数值仿真在复合材料领域已经是很成熟的手段，但真正做起来，渐进损伤模型的参数标定和单元失效控制两个环节最考验工程判断。

建模阶段先处理了几何和材料参数。单层厚度0.125mm，8层铺层总厚度1mm，用ANSYS ACP模块构建层合板模型。材料属性用客户提供的测试值：E₁=135 GPa，E₂=9.5 GPa，ν₁₂=0.31，G₁₂=5.2 GPa，横向拉伸强度X₂t=54 MPa，面内剪切强度S=72 MPa。但力学数值仿真中最容易出问题的不是材料参数本身，而是单元类型选择和网格密度——选了SOLID185层单元，每个单元厚度方向分1层，面内网格0.5mm。做了网格敏感性分析：0.3mm、0.5mm、0.8mm三种网格下极限强度分别为612 MPa、608 MPa、595 MPa，0.5mm和0.3mm差值小于1%，最终选0.5mm网格。

渐进损伤模型是项目的核心。用了Hashin失效准则，区分纤维拉压失效和基体拉压失效四种模式。材料刚度退化用瞬时退化方案——一旦某个积分点满足失效判据，对应方向的刚度直接降为原来的10%。但这个退化系数争议很大，有些文献用0（完全删除），有些用0.2。先用了0，结果应力-应变曲线在基体失效点出现剧烈跳跃，不收敛。调整到0.1后，曲线平滑了很多，但极限强度从580 MPa变到了608 MPa，差了5%。最终用0.05的退化系数，极限强度610 MPa，与实验值620 MPa偏差1.6%。

层间损伤也是一个关键因素。用了COHESIVE单元模拟层间界面，界面强度参数：法向强度30 MPa，切向强度45 MPa，断裂韧性GⅠc=0.28 N/mm，GⅡc=0.65 N/mm。在45/-45界面处最先出现分层损伤，载荷达到480 MPa时分层起始，最终失效时分层面积扩展到整个45/-45界面的35%。力学数值仿真的分层预测和客户做的C扫描检测结果对比，分层位置一致，面积偏差在15%以内。

计算效率和收敛性也值得一提。整个模型约12万个单元，显式动力学求解用了LS-DYNA，总计算时间8小时。隐式求解器（ANSYS Standard）在基体大量失效后不收敛，不得不切换到显式。这涉及一个工程判断：拉伸问题本应使用隐式，但渐进损伤中大量单元同时退化导致刚度矩阵奇异，隐式求解器无法处理。显式求解需要控制质量缩放系数，设为100倍实际质量，动能/内能比控制在5%以内，保证准静态分析的准确性。

最终交付了应力-应变曲线、各层损伤演化云图、层间分层分布图和极限强度预测值。力学数值仿真的结果让客户优化了铺层顺序：将[0/45/-45/90]ₛ改为[0/45/90/-45]ₛ，极限强度从610 MPa提升到648 MPa，提升6.2%。仿真的价值在于用数据驱动设计决策，而不是靠经验试错。

断裂模式分析进一步验证了仿真精度。通过损伤变量云图可以追踪渐进失效过程：基体损伤先在90°层出现（载荷约280 MPa），然后扩展到-45°层（350 MPa），纤维断裂在0°层最后发生（610 MPa）。这个损伤序列和超声C扫描的分层检测在层间位置上完全吻合。力学数值仿真的损伤演化分析还能输出声发射信号的模拟预测——通过单元失效时的应变能释放率估算声发射振铃计数，与实验声发射监测数据的相关系数达到0.87。这种多维度数据交叉验证是力学数值仿真交付质量的核心指标。

https://www.keyanxueshu.com/