层合板厚度对桥接和层间断裂韧性影响的研究.docxVIP

? ? 层合板厚度对桥接和层间断裂韧性影响的研究 ? ? 张永贵, 马启文, 刘俭辉 (兰州理工大学 机电工程学院, 甘肃 兰州 730050) 由于复合材料具有高比强度、高比刚度和优秀的抗疲劳性等特点[1-3]，所以被广泛应用于航空航天、机械和建筑等领域[4-5].对于实际工程中常用的复合材料层合板结构，铺层纤维的存在使得其纵向和横向力学性能得到显著增强.然而，因层合板的层间性能较差，使得其极易产生层间分层[6].因此，在对复合材料结构进行设计和强度分析的过程中，需要对复合材料结构的分层损伤起始与扩展过程开展更加深入的研究[7-8]. 复合材料的层间断裂通常伴随着桥接区的发展，在桥接区内存在完整的纤维桥接裂纹面.这种现象在复合材料结构的黏结区建模和损伤容限设计分析中非常重要.此外，在工艺上，层合板厚度过大会增加制造成本，而厚度过小则力学性能达不到相关要求.因此，讨论层合板厚度对力学性能，尤其是对桥接性能的影响，在当前研究中显得尤为重要. 关于碳/环氧复合材料层内断裂过程中R曲线对试样厚度依赖性的探究，已经存在相关的实验研究.Tamuzs等[9]通过光栅实验采集在桥接区末端测得的张开位移和不同裂纹增量所对应的能量释放率，通过能量释放率相对于张开位移的导数求出桥接应力.然而，这种方法容易受到误差率和张开位移实验数据拟合的影响.Frossard等[10]研究发现单层铺层厚度变薄会使层间断裂韧性增大，同时，单层厚度变化对桥接区长度和开口位移也会有所影响，但并未考虑等单层厚度下铺层数变化对层间断裂韧性的影响. 许多学者对纤维桥接和层间断裂韧性进行了探究.Sorensen等[11]和Stutz等[12]提出了运用迭代法来确定层合板的桥接应力法则，通过光纤光栅传感器(FBG)，沿分层方向进行精确应变测量，借助有限元模型，提取分层桥接法则.此方法具有表征复合材料中大规模纤维桥接的强健性与稳定性，但是所需参数较多.Suo等[4]提出了几何效应对R曲线的2个基本特征，即稳态损伤区随梁的厚度增加，稳态应变能释放率与厚度无关，但没有通过实验完全记录下来. 本文基于双线性cohesive内聚力模型，研究单向碳/环氧复合材料I型脱层中大尺度纤维桥接的试样厚度依赖性.结合双线性内聚力模型，借助ABAQUS进行有限元分析，确定不同厚度下层合板试样的载荷变化与层间断裂韧性，并研究相关的分层扩展行为. 1 层合板的数值分析 1.1 修正梁理论(MBT)方法 复合材料I型层间断裂韧性值采用双悬臂梁实验测量，此实验方法在国内已被纳入航空工业HB7402—1996[13]标准.双悬臂梁的应变能释放率G为 (1) 式中：a为分层长度；b为试样宽度；P为载荷；d为载荷施加点产生的位移；F和N′分别为考虑大位移和加载块影响时的修正因子；Δ为分层长度的修正量，该值一般通过最小二乘法拟合确定. 1.2 临界应变能释放率的数值计算 通过文献[14]研究纤维桥接的产生与断开，得到临界应变能释放率GIC与裂纹扩展长度Δα的关系，如图1所示.图中Gi,Init、Gi,Prop分别为初始断裂韧性和临界断裂韧性.由于会产生纤维桥接，所以当裂纹张口位移达到临界值lbz时断裂韧性达到稳定值，即层间裂纹扩展表现出R曲线特征.该曲线能够较好地考虑纤维桥接对应变能释放率的影响. 图1 纤维桥接作用下的R曲线特征图Fig.1 R-curve characteristic diagram under fiber bridging 通过该曲线，对裂纹面和裂尖附近进行J积分计算，得到临界应变能释放率为 (2) 式中：GIC为I型加载模式下的临界应变能释放率，即测量得到的I型层间断裂韧性；Gi,Init通过实验来确定；δ*为桥接区在预制裂纹末端处的开口位移；σb为桥接应力. 基于式(2)，对开口位移进行微分，能够得到桥接力的值. 1.3 双线性内聚力模型 基于损伤和断裂力学，采用牵引-分离的本构关系模拟分层界面处的力学响应，其优点是可以同时模拟分层的萌生与扩展，无需预制分层损伤[15].在双线性本构模型[16]中，当施加能量大于GTC时，损伤开始累积，如图2所示.图中，Kn为最大界面刚度，Nmax为最大界面强度.当内聚力单元上、下界面间的位移达到δ0时，刚度开始下降；达到δf时，内聚力单元完全失效.黏聚力与相对位移本构关系方程为 t=(1-D)KΔ1 (3) 式中：t为黏聚力；Δ1为相对位移；K为界面刚度； 图2 双线性本构模型 D为界面损伤状态变量. 1.4 内聚力单元失效准则 在建立内聚力单元时需要确定相关失效准则，复合材料分层的萌生可通过内聚力分量大小来确定.当界面的内聚力大小等于对应分层模式下的界面强度时，萌生分层.本文选择二次应力准则作为失效判据，其判据为 (4) 式中：σn、σs、σt分