玻璃/环氧复合材料悬臂板中裂纹取向效应的厚度依赖性反转的数值与实验研究

《Results in Engineering》：Numerical and Experimental Investigation on the Thickness-Dependent Reversal of Crack Orientation Effects in Glass/Epoxy Composite Cantilever Plates

　　

本研究针对正交各向异性FR4/玻璃/环氧复合材料悬臂板中板厚与裂纹取向角对动态振动特性的耦合影响进行了调查，尽管对单个裂纹参数已有广泛研究，但这一关系仍未被充分理解。传统研究孤立地分析这些损伤参数，而本研究系统评估了边缘裂纹的特征，包括不同长度、位置和取向角（90°、45°和135°）。在此背景下，研究人员研究了两种板厚（3 mm和5 mm）对前五阶自然频率及其对应模态振型变化的影响。在数值分析中，研究人员采用有限元法（FEM），借助ANSYS软件中的SHELL181单元。在实验分析中，研究人员使用冲击锤激励系统与压电加速度计和数据采集系统（DAQ）。研究发现，数值与实验结果之间存在极佳的相关性，并明确指出边缘裂纹的存在总是由于弯曲刚度的降低而导致自然频率值下降。高阶模态对边缘裂纹的敏感性高于低阶模态。结果表明，垂直于板轴的裂纹方向（90°）对频率降低影响最大，因为它直接阻断了主应力线。当裂纹靠近临界固定支撑区域（1/4）时，这种影响更为显著。关键的是，研究人员识别出裂纹取向存在明显的厚度依赖性行为：在较薄板（3 mm）中，45°裂纹配置导致最小的频率退化，而在较厚板（5 mm）中，135°角对频率的影响最小。这些发现表明，结构厚度改变了系统对取向特定刚度损失的敏感性，为增强先进复合材料构件中高保真结构健康监测（SHM）算法提供了基准数据集。

**论文解读：玻璃/环氧复合材料悬臂板中裂纹取向效应的厚度依赖性反转**

**研究背景与问题**

复合材料因其高强度重量比、抗热降解、耐腐蚀及可定向设计性能，在航空航天、汽车、建筑和新能源等领域广泛应用。其中，玻璃纤维增强聚合物（GFRP）因优异的力学性能、高损伤容限和成本效益而备受青睐。然而，复合材料结构在服役过程中易产生裂纹等缺陷，早期裂纹难以察觉且会显著降低安全性与寿命。基于振动的损伤检测技术和结构健康监测（SHM）利用模态参数（如自然频率和模态振型）的变化来识别损伤，已成为有效的无损检测方法。尽管已有大量研究探讨了裂纹参数（如长度、位置、深度）对振动特性的影响，但裂纹取向与板厚之间的耦合作用仍未被充分理解。传统研究多孤立分析单个参数，忽略了二者交互效应对高阶模态敏感性的影响。因此，本研究旨在揭示裂纹取向角与板厚如何共同影响正交各向异性复合材料悬臂板的动态响应，为提升SHM算法精度提供关键数据。

**研究内容与结论**

研究人员通过数值模拟（有限元法，FEM）与实验测试相结合的方法，系统研究了边缘裂纹长度、位置和取向角（90°、45°、135°）对两种厚度（3 mm和5 mm）FR4/玻璃/环氧层压板前五阶自然频率及模态振型的影响。结果表明：边缘裂纹总是降低自然频率，且高阶模态更敏感；90°裂纹（垂直于板轴）因直接阻断纵向纤维导致最大频率降低，且裂纹越靠近固定端（1/4位置）影响越显著。关键发现是裂纹取向对频率的影响存在厚度依赖性反转：在3 mm薄板中，45°裂纹影响最小；在5 mm厚板中，135°裂纹影响最小。通过分析4 mm中间厚度板，验证了该反转行为存在临界过渡厚度。该研究为考虑厚度-取向耦合的SHM算法提供了基准数据集，论文发表在《Results in Engineering》。

**主要技术方法**

数值分析采用ANSYS APDL软件中的SHELL181单元（基于一阶剪切变形理论，FSDT），对正交各向异性复合材料进行建模，边缘裂纹作为物理几何不连续体直接建模，并在裂纹尖端进行局部网格细化（2 mm单元尺寸经收敛性验证）。实验测试使用从深圳雄毅华绝缘有限公司采购的标准化FR4玻璃/环氧层压板（3 mm和5 mm厚），通过冲击锤激励系统和压电加速度计（12 g）及数据采集系统（DAQ）进行实验模态分析（EMA），获取前五阶自然频率和模态振型。所有实验边界条件通过重型钢制夹具施加固定支撑。

**研究结果**

**1. 整体动态敏感性**
通过数值与实验对比，有限元模型与实验数据高度吻合（最大偏差3 mm板为5.964%，5 mm板为1.542%）。边缘裂纹始终降低自然频率，且高阶模态（第4、5阶）比低阶模态更敏感，裂纹靠近固定端（x_c/L=1/4）时影响最大。

**2. 裂纹取向与厚度的交互作用**
90°裂纹在所有厚度下均导致最大频率降低（例如3 mm板第5模态降低26.01 Hz，5 mm板第5模态降低43.56 Hz）。但45°和135°裂纹的影响随厚度反转：在3 mm板中，45°裂纹影响最小（频率降低0.0158–10.28 Hz）；在5 mm板中，135°裂纹影响最小（频率降低0.052–9.02 Hz）。通过4 mm板数值模拟，观察到过渡行为（45°和135°曲线交替影响），确认h=4 mm为临界厚度。

**3. 应变能分布分析**
应变能分布模式不随板厚改变，但定量值随厚度增加（3 mm板裂纹尖端峰值应变能密度0.462×10⁹ J/m³，5 mm板为0.763×10⁹ J/m³，增加65.15%）。模态分析表明，裂纹取向的影响取决于裂纹与高应变能区域的相对位置：例如，第1模态应变能集中于固定端，远离该区域的裂纹取向影响小；第2模态应变能分布于固定端和跨中，45°和90°裂纹在靠近自由端时影响显著；第3模态（扭转）中，135°裂纹在自由边界影响弱；第4模态（三弯）中，90°裂纹在自由边界影响最大；第5模态（双扭转）中，90°和135°裂纹在边界和中心位置影响较大。

**讨论与结论**

研究讨论部分指出，厚度依赖性反转的物理机制源于弯曲刚度（D₁₁）和扭转刚度（D₆₆）随厚度立方（h³）缩放。3 mm板D₁₁=46.957 N·m，5 mm板D₁₁=217.391 N·m（增加4.63倍），导致厚板中弯曲变形能量壁垒升高，应变能重新分布，使系统动态行为从弯曲主导转向剪切/扭转主导，从而改变裂纹取向的影响顺序。结论部分总结了三点：（1）边缘裂纹降低自然频率，高阶模态及高应变能区域更敏感，固定端边界（x_c/L=1/4）为最危险区域；（2）90°裂纹因纵向纤维断裂导致最大频率降低，而45°与135°的影响依赖于厚度：3 mm板中45°影响最小，5 mm板中135°影响最小，驱动力为横向剪切变形和弯曲-扭转刚度耦合变化；（3）裂纹正交各向异性板的动态响应由裂纹取向、位置、模态应变能与板厚耦合控制，纳入厚度-取向相关性对增强基于振动的SHM框架至关重要。