Processing math: 0%

CFRP/Al抽芯铆钉单剪干涉配合的渐进损伤数值模拟研究

王兵兵, 周钊元, 金万军, 何超, 唐正强

王兵兵, 周钊元, 金万军, 等. CFRP/Al抽芯铆钉单剪干涉配合的渐进损伤数值模拟研究[J]. 复合材料学报, 2024, 41(6): 3258-3270. DOI: 10.13801/j.cnki.fhclxb.20231025.001
引用本文: 王兵兵, 周钊元, 金万军, 等. CFRP/Al抽芯铆钉单剪干涉配合的渐进损伤数值模拟研究[J]. 复合材料学报, 2024, 41(6): 3258-3270. DOI: 10.13801/j.cnki.fhclxb.20231025.001
WANG Bingbing, ZHOU Zhaoyuan, JIN Wanjun, et al. Numerical simulation of progressive damage of single-lap CFRP/Al connected by blind rivet under interference condition[J]. Acta Materiae Compositae Sinica, 2024, 41(6): 3258-3270. DOI: 10.13801/j.cnki.fhclxb.20231025.001
Citation: WANG Bingbing, ZHOU Zhaoyuan, JIN Wanjun, et al. Numerical simulation of progressive damage of single-lap CFRP/Al connected by blind rivet under interference condition[J]. Acta Materiae Compositae Sinica, 2024, 41(6): 3258-3270. DOI: 10.13801/j.cnki.fhclxb.20231025.001

CFRP/Al抽芯铆钉单剪干涉配合的渐进损伤数值模拟研究

基金项目: 贵州省科技计划项目(黔科合支撑[2021]一般281)
详细信息
    通讯作者:

    唐正强,博士,教授,博士生导师,研究方向为智能装备、摩擦磨损、防腐耐磨涂层、功能表面涂层 E-mail: zhengqiagntang@126.com

  • 中图分类号: TH131.1;TB333

Numerical simulation of progressive damage of single-lap CFRP/Al connected by blind rivet under interference condition

Funds: Guizhou Science and Technology Plan Project (Guizhou Science and Technology Support [2021] General 281)
  • 摘要: 干涉配合在改善抽芯铆钉连接复合材料的力学性能方面具有较强优势。然而,干涉配合下抽芯铆钉的安装往往会导致碳纤维增强树脂复合材料(CFRP)层压板的损坏,削弱铆接接头的力学性能。本文采用了完整的抽芯铆钉模型,基于连续损伤力学、扩展的三维破坏准则和应变率效应,通过有限元和Abaqus二次开发相结合的方法,研究了CFRP/Al单剪铆接中抽芯铆钉与连接孔的干涉量和铆钉的安装速度对复合材料损伤的影响。从模拟结果来看受安装速度的影响显著,安装阻力呈现出两个典型的阶段,分别由铝合金板、CFRP板与抽芯铆钉的摩擦产生造成。抽芯铆钉安装在干涉配合情况下,安装速度越大对于降低安装阻力十分有利。然而,过高的安装速度会导致孔壁的损伤增加,尤其在高干涉量下更明显。抽芯铆钉铆接接头的失效模式以CFRP损伤为主,并受干涉尺寸的影响显著。

     

    Abstract: Interference fit has advantages in improving the mechanical performance of blind rivet in connections for composite materials. However, the installation of blind rivets under interference fit leads to the damage of carbon fiber reinforced plastic (CFRP) laminates, weakening the mechanical performance of the joints. In this paper, based on continuous damage mechanics, extended three-dimensional failure criteria and strain rate effect, a complete blind rivet model, which combining finite element method and secondary development of Abaqus was developed to study the effect of interference and installation speed on the damage of composite materials. From the simulation results, the installation resistance force exhibits two typical stages, which are generated by the friction between the aluminum alloy and the CFRP laminate with the blind rivet, respectively. Under interference fit conditions, a higher installation speed is highly advantageous in reducing installation resistance. However, excessive installation speed can lead to increased damage to the hole walls, which is especially noticeable in cases of high interference fit. The failure mode of the joint is mainly caused by CFRP damage and is significantly influenced by the interference size.

     

  • 碳纤维增强树脂复合材料(CFRP)具有轻质、高强度、高耐腐蚀和抗疲劳性、低蠕变、高阻尼等优点,被广泛用于航空航天、汽车、高铁等众多工业领域[1-3]。碳纤维复合材料间的连接主要采用胶接、铆接、螺栓连接及其混合连接[4],但由于行业发展,多种不同材料之间的相互连接需求越来越多,而机械连接具有较高的可靠性、稳定性、对湿度和温度不敏感等优点,成为不同材料连接中最优的选择[5]。在机械连接中,螺栓连接强度高、防腐性好,但是在复合材料部件上会产生蠕变现象,会导致紧螺栓松动,从而导致连接故障[6-7]。传统的捶铆是一种通过捶打铆钉使其产生变形,从而固定两个或多个工件的方法[8]。随着科技的发展,相继出现了压铆、抽芯铆、无铆和自冲孔铆接等新型铆接工艺,如带芯杆的抽芯铆钉和带有环槽的环槽铆钉等[9],弥补了传统铆钉中空所导致剪切强度底的缺点[10-11]。这些铆接方法的连接机制都是通过自身的变形来阻止连接件移动,故铆接强度一方面决定于铆钉自身的变形能力和强度,另一方面决定于连接件的材料和几何形状[12-14]。对碳纤维复合材料,与传统金属材料不同,孔的加工导致了复合材料的纤维被切断,这样孔周围极易产生微裂纹,并在开口处形成高应力区域,成为铆接接头的薄弱环节。同时,复合材料的各向异性会在加载过程中在开口区域引起边缘效应,形成应力集中,从而降低复合材料的强度[15-16]。在早期研究中,通常用螺栓连接复合材料并采用间隙配合,通过进一步研究,证明过盈配合可以通过改善螺栓孔接触条件来提高接头的承载能力。此外,适当过盈配合尺寸的接头也显示出较高的抗疲劳性[17-19]。干涉配合利用直径大于连接孔径的紧固件实施接头连接,利用干涉量增加有效接触面积进而提高接头静强度,利用有益残余应力改善接头疲劳性能。Zuo等[20]通过研究CFRP/CFRP过盈配合螺栓接头的螺栓插入损伤,证明了过盈配合在改善复合材料螺栓连接的疲劳行为方面具有优势。Hu等[21]通过摩擦测试和剪切测试,优化了复合材料过盈配合接头的设计和制造工艺,提高接头的可靠性和性能。Li等[22]通过实验研究了过盈配合尺寸对CFRP/Ti合金螺栓连接微动疲劳寿命、孔伸长率和圆周应变的影响,得到了不同应力水平下的过盈配合尺寸与微动疲劳寿命之间的关系及孔伸长率变化。Li等[23]通过对复合材料铆接接头的拉伸响应和失效的实验研究表明,复合材料接头的平均总能量吸收随着加载速率的增加而增加。Zhao等[24]通过采用经过验证的模型进行参数分析,表明在铆钉空间增加到铆钉直径的3.5倍之前,极限载荷就增加了,抗剪承载力随着盲铆钉数量的增加而线性增加。Xing等[25]对CFRP/CFRP铆接接头进行了拉张疲劳试验,估计了疲劳寿命,结果表明在极限抗拉强度为38.3%的情况下,铆接CFRP/CFRP接头可以达到1×106的无限循环寿命,这意味着在循环加载下,接头可以承受至少100万次的循环载荷而不发生破坏。Kelly等[26]对铆钉设置过程中涉及的成形力进行了研究,观察到固定力与固定铆钉的不对称性之间的关系。

    然而,抽芯铆钉连接复合材料的损伤问题是一个重要且复杂的研究课题。在过去的研究中,学者们通常将铆钉简化为哑铃状,这种简化模型只能够提供一些关于连接性能的信息,但无法完全准确地预测抽芯铆钉连接面临的损伤问题。为了更详细地研究抽芯铆钉连接复合材料的损伤问题,本文采用了完整的抽芯铆钉模型。这种模型考虑了铆钉内部的抽芯结构,及与复合材料之间的相互作用。通过使用这个更复杂的模型,可以更准确地描述抽芯铆钉连接复合材料的损伤问题。

    本文首先建立了CFRP/Al抽芯铆钉铆接的完整三维几何模型,并在Abaqus/Explicit中进行了模拟。模型中考虑了材料的非线性本构关系和破坏准则,并利用VUMAT用户子程序实现了渐进损伤模型。通过控制不同的安装速度和干涉量,分析了抽芯铆钉连接复合材料时的安装阻力和接头变形情况。其次,还详细分析了CFRP材料的损伤情况,通过模拟结果,可以观察到在抽芯铆钉安装和锁紧过程中,CFRP材料发生了局部损伤,如纤维破裂和基体开裂。通过对模拟结果和文献中的实验结论进行对比,可以验证模型的准确性和可靠性。

    随着复合材料在航天航空等领域的大量使用,数值算可以对CFRP铆接极端工况的力学性能进行评价,可以精准的预测各向异性材料的损伤,节约大量的实验成本。本文采用渐进失效的分析方法对复合材料层压板铆接结构进行分析。该方法可以预测起始和局部的渐进损伤、破坏模式和铆接接头的极限强度。

    常用的失效判据有5种,分别为最大应力准则、Tsai-Hill失效判据、Hoffman失效准则、Tsai-Wu失效准则和Hashin失效准则[27-28]。在这些准则中,为了更准确地描述复合材料层压板的损坏行为,本研究采用了改进的三维Hashin失效准则,该准则可适用于多种损坏模式,包括纤维拉伸破坏、纤维压缩破坏、基体拉伸破坏、基体压缩破坏、拉伸分层和压缩分层,这些损伤模式在Abaqus软件中用SDV1~SDV6表示,SDV1、SDV2、SDV3、SDV4、SDV5和SDV6分别代表纤维拉伸破坏、纤维压缩破坏、基体拉伸破坏、基体压缩破坏、拉伸分层和压缩分层,失效准则的表达式如下所示:

    (1) 纤维拉伸失效(σ1>0)

    Rft=(σ1C11XT)2+(τ12C44S12)2+(τ13C55S13)2 (1)

    (2) 纤维压缩失效(σ1≤0)

    {R_{{\text{fc}}}} = \sqrt {{{\left( {\frac{{{\sigma _1}{C_{11}}}}{{{X_{\text{C}}}}}} \right)}^2}} \geqslant 1 (2)

    (3) 基体拉伸失效(σ2>0)

    {R_{{\text{mt}}}} = \sqrt {{{\left( {\frac{{{\sigma _2}{C_{22}}}}{{{Y_{\text{T}}}}}} \right)}^2} + {{\left( {\frac{{{\tau _{12}}{C_{44}}}}{{{S_{12}}}}} \right)}^2} + {{\left( {\frac{{{\tau _{13}}{C_{66}}}}{{{S_{23}}}}} \right)}^2}} \geqslant 1 (3)

    (4) 基体压缩失效(σ2≤0)

    {R_{{\text{mc}}}} = \sqrt {{{\left( {\frac{{{\sigma _2}{C_{22}}}}{{{Y_{\text{C}}}}}} \right)}^2} + {{\left( {\frac{{{\tau _{12}}{C_{44}}}}{{{S_{12}}}}} \right)}^2} + {{\left( {\frac{{{\tau _{13}}{C_{66}}}}{{{S_{23}}}}} \right)}^2}} \geqslant 1 (4)

    (5) 拉伸分层失效(σ3>0)

    {R_{{\text{lt}}}} = \sqrt {{{\left( {\frac{{{\sigma _3}{E_3}}}{{{Z_{\text{T}}}}}} \right)}^2} + {{\left( {\frac{{{\tau _{23}}{G_{23}}}}{{{S_{23}}}}} \right)}^2} + {{\left( {\frac{{{\tau _{13}}{G_{13}}}}{{{S_{13}}}}} \right)}^2}} \geqslant 1 (5)

    (6) 压缩分层失效(σ3≤0)

    {R_{{\text{lc}}}} = \sqrt {{{\left( {\frac{{{\sigma _3}{E_3}}}{{{Z_{\text{C}}}}}} \right)}^2} + {{\left( {\frac{{{\tau _{23}}{G_{23}}}}{{{S_{23}}}}} \right)}^2} + {{\left( {\frac{{{\tau _{13}}{G_{13}}}}{{{S_{13}}}}} \right)}^2}} \geqslant 1 (6)

    其中:RftRfcRmtRmcRltRlc为相对应的失效系数;{\sigma _1}{\sigma _2}{\sigma _3}为3个方向的正应力;{\tau _{12}}{\tau _{13}}{\tau _{23}}为3个方向的剪切应力;{X_{\text{T}}} {X_{\text{C}}} 为单层板的纤维方向拉伸强度和压缩强度;{Y_{\text{T}}}{Y_{\text{C}}}为单层板的基体方向拉伸强度和压缩强度;{Z_{\text{T}}}{Z_{\text{C}}}为单层板的面外拉伸强度和压缩强度;{S_{12}}{S_{13}}{S_{23}}为单层板对应方向上的剪切强度;Cij (i, j = 1, 2, 3, 4, 5, 6)为刚度矩阵的系数;G13G23分别为x-z方向、y-z方向的剪切模量;E3为3方向(面外)的弹性模量。

    CFRP是一种各向异性材料,在受到超负荷时,CFRP材料会出现损伤,导致其承载能力下降,进而影响材料的性能。图1展示了CFRP材料力学性能的3种退化模型。突然退化模型是当CFRP材料在达到一定应力或应变阈值后,其力学性能会突然退化;线性逐渐退化模型表示CFRP材料随着应力或应变的增加,材料的承载能力会线性下降,性能退化的速率是恒定的;非线性逐渐退化模型则表示CFRP材料随着应力或应变的增加,材料的承载能力下降速率逐渐增加。这种非线性退化模型更符合实际情况,原因是材料中的损伤和退化过程往往是复杂的,并且退化速率可能会随着过载程度的增加而加快。因此本文选择非线性逐渐退化模型,可以更准确地描述CFRP材料在过载情况下的力学性能变化。

    图  1  力学性能3种退化方式
    Figure  1.  Three degradation modes of mechanical properties

    考虑到抽芯铆钉在高速安装复合材料孔中所显示的应变率效应,需要考虑应变率对材料性能参数的影响。在材料的本构方程中,通常会通过添加修正项来描述应变率效应,该修正项是根据前人对于应变率效应的研究方法进行借鉴的[29],本文在子程序中引入了动态增加因子(Dynamic increase factor,DIF),复合材料单层板各方向上前应变率下的强度矩阵:

    S_{\text{DIF}}^{\text{S}}=S_0\left(1+C\rm{ln}\frac{|\dot{\varepsilon}|}{\dot{\varepsilon}_0}\right) (7)
    \{ {S_{{\text{RF}}}}\} = \{ {S_0}\} S_{{\text{DIF}}}^{\text{S}} (8)

    其中: {S_0} 是参考应变率下的强度矩阵,也就是准静态下的强度矩阵; {\dot \varepsilon _0} 为准静态下的应变率值,根据文献取值10−3 s−1 S_{{\text{DIF}}}^{\text{S}} 为当前应变率下的强度矩阵,对应的 |\dot \varepsilon | 为当前状态下的应变率值;C为应变率硬化修正系数,本文借鉴Naik等[30]的试验研究结果,取C=0.1;{S_{{\text{RF}}}}为当前应变率下的强度。

    由于本文所研究的抽芯铆钉安装距离只有14 mm,抽芯铆钉安装速度的差异而产生的热量差异较小,并且复合材料具有较低的热导率和较低的热膨胀系数。因此,该热量差异对复合材料损伤影响较小,因此本文忽略了热效应的影响。

    CFRP铆接模拟实施流程如图2所示,首先根据改进的Hashin损伤准则,通过Fortran语言编写Abaqus软件用户子程序VUMAT;随后,子程序读取CFRP材料参数,并进行刚度矩阵和应力应变的计算,根据计算损伤指数来判断材料是否失效,将结果反馈状态变量;同时,利用Abaqus软件建立结构的几何模型,包括需要进行材料定义的部分,对模型进行网格划分,并设置合适的边界条件,在模型中引入刚才编写的VUMAT子程序,并将材料参数传递给子程序进行计算。

    图  2  碳纤维增强树脂复合材料(CFRP)铆接模拟实施流程
    Figure  2.  Simulation implementation process of carbon fiber reinforced plastics (CFRP) riveting
    σ—Stress vector; C—Stiffness matrix of undamaged orthotropic material; ε—Strain vector

    该模型的CFRP层压板和铝合金板根据ASTM D5961[31]标准设计,其几何参数如图3所示,长度为110 mm,宽度为30 mm,孔径为4.8 mm。本文使用的碳纤维增强材料名称为T700/M21[32]层压板,该层压板的铺层顺序为[90/0/90/0/90/0/90/0/90/0/90/0/90/0/90/0/90/0/90/0],包含50%的0°铺层和50%的90°铺层,层压板的厚度为3 mm,单层厚度为0.15 mm。T700/M21的具体属性参数和材料强度如表1表2所示。该模型使用的抽芯铆钉为航天用铆钉,抽芯铆钉由钉套、芯杆、锁环和推压衬套组成,各个零部件和铝合金的材料参数如表3所示。

    图  3  CFRP/Al单钉单剪铆接几何参数
    Figure  3.  Geometrical parameters of CFRP/AL single nail single shear riveting
    表  1  T700/M21复合材料属性参数
    Table  1.  Property parameters of T700/M21 composites
    E11/GPa E22/GPa E33/GPa G12/GPa G13/GPa G23/GPa v12 v13 v23
    117 7.7 7.7 4.8 4.8 2.8 0.34 0.34 0.4
    Notes: E11, E22, E33—Elastic modulus in the 1 (axial), 2 (transverse), and 3 (out-of-plane) directions; G12, G13, G23—Shear modulus for the 12, 13, and 23 planes; v12, v13, v23—Principal Poisson's ratios for the 12, 13, and 23 planes.
    下载: 导出CSV 
    | 显示表格
    表  2  T700/M21复合材料层间强度
    Table  2.  Interlaminar strength of T700/M21 composites
    XT/MPa XC/MPa YT/MPa YC/MPa ZT/MPa ZC/MPa S12/MPa S13/MPa S23/MPa
    2200 1700 50 300 50 300 95 95 95
    Notes: XT, XC—Tensile and compressive strengths in the fiber direction; YT, YC—Tensile and compressive strengths in the matrix direction; ZT, ZC—Tensile and compressive strengths in the out-of-plane direction; S12, S13, S23—Shear strengths for the 12, 13, and 23 planes.
    下载: 导出CSV 
    | 显示表格
    表  3  抽芯铆钉各零件和铝合金的属性参数
    Table  3.  Property parameters of each part of blind rivet and aluminum alloy
    Name Type Young's modulus/GPa Poisson's ratio Density/(g·cm−3)
    Sleeve CPTi 100 0.34 4.51
    Stem Ti3Al8V6Cr4Mo4Zr 107.8 0.32 4.82
    Lock collar A286 2.01 0.305 7.49
    Shifting anvil 45A 206 0.269 7.85
    Al alloy Al 72 0.33 2.70
    下载: 导出CSV 
    | 显示表格

    表4为CFRP/Al抽芯铆钉相对干涉量。在抽芯铆钉安装模拟中考虑了3组过盈配合尺寸,以了解干涉配合条件下CFRP层压板的损伤行为。定义相对干涉尺寸I

    I = \frac{{D - d}}{d} \times 100{{\% }} (9)

    其中:d为CFRP与铝合金的孔径;D为抽芯铆钉的钉套直径。

    表  4  CFRP/Al抽芯铆钉相对干涉量I大小
    Table  4.  Relative interference magnitude I of the CFRP/Al blind rivet
    Number Interference amount/% Sleeve diameter/mm Aperture/mm
    I0 0 4.80 4.80
    I1 0.42 4.80 4.78
    I2 0.84 4.80 4.76
    I3 1.30 4.80 4.74
    下载: 导出CSV 
    | 显示表格

    此模型在Abaqus中选取Abaqus/Explicit显式动力学方法,由于抽芯铆钉锁定时钉套会发生大变形,会产生复杂的接触,需要大量的迭代来满足接触条件,因此隐式动力学和静力学收敛困难。在计算过程中为了避免零件之间发生穿透,采用4对接触,其接触属性为切向“罚”接触与法向“硬”接触,分别是通用接触、芯杆头部和CFRP内孔接触、芯杆头部与铝合金板内孔接触和CFRP与铝合金板接触,摩擦系数采用文献[33]中实验所得,分别是0.1、0.2、0.2、0.3。

    三维有限元模型的边界设置如图4所示,显示了网格和边界条件。在抽芯铆钉的安装过程中,CFRP与Al板固定,限制抽芯铆钉XZ轴的自由度。其中1/4的模型被隐藏,以显示内部结构。为了建立一个准确的抽芯铆钉模型,预先使用Solidworks软件建立了模型并将其导入Abaqus。抽芯铆钉的铆接过程如图5所示,可以被描述为以下5个步骤:首先,第一步,将抽芯铆钉插入CFRP孔中,铆钉的尾部放置在需要铆接的两个工件之间;第二和第三步,通过施加力量,芯杆开始被拉动,带动钉套向内收缩,在钉套收缩过程中,会产生径向力,将连接件夹紧在一起,以产生套筒的变形紧固效果,如图5(b)图5(c)所示;第四步是通过向下挤压和推压衬套,使锁环嵌入钉套头部,从而固定铆钉,这一步骤可以在图5(d)中观察到;最后,第五步是继续拉动芯杆,直到芯杆从缺口处拉断,完成了铆接过程,见图5(e)

    为了获得准确的模拟结果,实现更高的模拟精度,进行了网格优化设计分析以确定最优的全局网格尺寸,网格的全局大小选择为0.6。网格使用8节点缩减积分实体单元(C3D8R)实体单元建模,并加强沙漏控制。因为抽芯铆钉为高度对称模型,所以采用中心轴划分网格。并根据网格划分准则,在CFRP与铝合金板孔周围进行了网格细化,如图4所示。

    图  4  网格划分及边界条件
    Figure  4.  Grid division and boundary conditions
    Ux, Uy, Uz—Displacement along the X-axis, Y-axis, Z-axis
    图  5  抽芯铆钉的铆接过程
    Figure  5.  Riveting process of the blind rivet

    模拟过程分为两步,第一步模拟了3种过盈配合尺寸(即I1I2I3)下抽芯铆钉的安装过程;第二步是对抽芯铆钉芯杆的拉伸和钉套锁紧过程进行了研究,可以获得钉套变形和锁环挤压过程的模拟结果。基于上述模拟结果,通过改变边界条件,研究抽芯铆钉以不同的安装速度和干涉量下安装与锁定过程对复合材料的损伤和力学性能,其次是钉套变形形成“鼓包”对最后一层复合材料的损伤的研究。在数值计算过程中,调用VUMAT子程序来评估仿真模型中损伤的起始和演变。

    抽芯铆钉的安装阻力Rf

    {R_{\text{f}}} = {F_1} + {F_2} (10)

    其中:{F_1}是芯杆头部和钉套与孔壁之间的摩擦力;{F_2}是孔变形引起的剪切阻力和弹塑性变形引起的阻力。

    通过模拟结果,得到了不同干涉量下的位移载荷曲线,不同干涉量和安装速度下的安装阻力不同。3种干涉量下的位移-载荷曲线如图6所示,通过对抽芯铆钉安装的过程和安装阻力的分析,整个过程分为6个阶段如图7所示。第1阶段(点O至点A)。由于惯性作用,这一阶段的安装阻力快速上升,安装阻力是由芯杆尾部和铝合金孔壁的摩擦力引起的。第2阶段(点A至点B),由于第一阶段的惯性作用,导致铝合金孔壁附近的材料发生弹塑性变形,因此安装阻力在这一阶段下降。第3阶段(点B至点C),这一阶段安装阻力增大是由于铝合金弹塑性变形产生了回弹,主要阻力是芯杆头部与铝合金孔壁摩擦和钉套与铝合金孔壁摩擦所引起的。第4阶段(点C至点D),安装阻力上升速度减缓,芯杆头部与钉套完全进入铝合金孔壁,由于安装速度过高,导致摩擦力略微减少,从而减缓了安装阻力的上升速度。第5阶段(点D至点E),这一阶段安装阻力上升速度再次增加,是由于芯杆头部开始与CFRP孔壁接触,复合材料与铝合金材料属性的不同导致安装阻力上升速度增加,这一阶段的主要摩擦力来自芯杆头部与CFRP孔壁之间的摩擦和钉套与铝合金之间的摩擦。第6阶段(点E至点F),这一阶段安装阻力开始下降,由于抽芯铆钉与CFRP孔壁的摩擦力过大,导致CFRP孔壁附近出现不可逆转的损伤,材料的损伤会减弱CFRP力学性能。

    图  6  CFRP/Al单钉单剪铆接在干涉量I1I2I3下的位移-载荷曲线
    Figure  6.  Displacement-load curves of CFRP/Al single-pin single-shear riveting with interference quantities I1, I2 and I3
    图  7  过盈配合接头的载荷-位移曲线的典型模型
    Figure  7.  Typical model of load-displacement curve of interference fit joint

    从干涉量来看,在低干涉量I1下,随着速度的增加平均安装阻力落差较小,主要由于此时干涉量较小,安装阻力整体水平较低,导致安装方法的影响效果不能体现。随着干涉量进一步上升到I2I3,安装阻力也随之增加,并且随着速度的增加平均安装阻力落差较大。

    从安装速度来看,安装速度越高,安装阻力越低,其原因有4点:第一是由于干涉量越大使摩擦系数的影响加大,高的安装速度一定程度上降低了摩擦系数,从而降低了安装阻力;第二是应变率效应。应变率是指材料中变形的速率,而在高速安装过程中,铆钉与CFRP孔壁之间会发生摩擦和转动,这会导致局部应变率的快速变化。当铆钉以高速插入过盈的复合材料孔时,应变率会随着速度增加而增加。应变率效应会导致材料的应力集中,但随着速度的增加,应变率的变化更加迅速,这可以减轻应力集中的程度,从而降低安装阻力。换句话说,高速安装可以减少由于应变率效应引起的额外摩擦和转动阻力;第三是由于孔壁的弹性变形。高速安装会导致螺栓和孔壁的弹性变形,使它们之间的接触面积增加,从而减少局部接触应力,降低安装阻力;第四是由于惯性力效应。高速安装时,铆钉和孔壁受到的惯性力会对安装过程产生一定的冲击和震动,可能有助于铆钉的安装。

    此外,安装的速度越大,安装阻力将更早达到峰值并开始下降,这很有可能是由于抽芯铆钉上的应力波通过摩擦传递到铝合金和CFRP孔壁未变形区域,在抽芯铆钉安装过程中,抽芯铆钉的尖端先进入孔中,然后通过在一侧形成断裂,这会产生巨大的应力波。当安装速度较大时,这种应力波的传播速度也会增加。应力波会通过与铝合金和CFRP孔壁接触的摩擦力传递,导致孔壁发生受压变形。当应力波传递到孔壁的未变形区域时,会产生与传递应力成反比的变形应力。随着安装速度的增加,这些变形应力会更早达到峰值,并随后开始下降。这是由于在高速安装过程中,压力被更快地释放,并在较短的时间内使孔壁达到稳定状态。表明抽芯铆钉安装在干涉配合情况下,较大的安装速度对于降低安装阻力十分有利。

    左杨杰[34]通过实验的方法研究了CFRP/Ti干涉配合螺栓连接动静态安装损伤分析,得到了高速安装的方法优于低速安装,该结果与本文的模拟结果一致。

    抽芯铆钉在安装过程中,随着安装速度的增加,虽然安装阻力变小,但是过高的安装速度会导致CFRP的孔壁损伤增加。图8为在干涉量I1I2I3下的损伤云图。可以看出,3种干涉量下纤维损伤、基体损伤、层合板拉伸损伤和层合板压缩损伤都随着安装速度的增加而增加,尤其在高干涉量I3下最为明显。其主要由以下两个因素所导致,第一是冲击力,较高的安装速度会增加抽芯铆钉进入CFRP孔时产生的冲击力。这种冲击力可能会超过复合材料的耐受能力,导致表面或内部出现裂纹;第二是摩擦热,较高的安装速度会产生更多的摩擦热,摩擦热会引起复合材料的热损伤。热损伤导致复合材料的软化、熔化或炭化,降低其强度和耐久性。

    图  8  干涉量I1I2I3下1 m/s、5 m/s和10 m/s所对应的CFRP的纤维拉伸损伤、基体拉伸损伤、层合板拉伸损伤和层合板压缩损伤
    Figure  8.  Fiber tensile damage, matrix tensile damage, laminate tensile damage and laminate compression damage of CFRP corresponding to 1 m/s, 5 m/s and 10 m/s under interference quantities I1, I2 and I3
    SDV1—Fiber tensile damage; SDV3—Matrix tensile damage; SDV5—Tensile delamination damage; SDV6—Compressive delamination damage

    不同干涉量下抽芯铆钉的安装所造成的CFRP的损伤情况各不相同,如图9所示。在I1干涉量下,铝合金发生了轻微弹塑性变形,CFRP层间和层内的损伤值很小,孔壁附近只发生轻微变形,对材料性能几乎没有影响。在I2干涉量下,随着干涉量的增加,CFRP孔壁发生明显的塑性变形,尤其在CFRP层压板的首层,发生了大量的层内损伤裂纹,层间也产生了大量的损伤,这是由于纤维和基体的脱胶引起。在I3干涉量下,孔壁附近产生了明显的变形,CFRP的损伤更加严重,更多的CFRP层压板内部开始萌生基体裂纹,同时基体裂纹由孔壁附近向层压板内部更深处扩展。

    图  9  干涉量I1I2I3下的纤维拉伸损伤、基体拉伸损伤、拉伸分层损伤和压缩分层损伤
    Figure  9.  Fiber tensile damage, matrix tensile damage, tensile delamination damage and compressive delamination damage under interference quantities I1, I2 and I3

    在不同的干涉量下完成抽芯铆钉的安装之后,通过拉伸芯杆和按压衬套对抽芯铆钉进行锁紧,模拟结果的应力云图如图10所示。可以看出,安装的损伤主要集中在孔壁,而锁定的损伤主要位于靠近钉套折叠位置的区域,也就是形成“鼓包”时,会对CFRP底层造成微量损伤。

    图  10  CFRP/Al抽芯铆钉锁定的应力云图
    Figure  10.  Stress nebulogram of a locked CFRP/Al rivet
    S—Von mises stress (MPa)

    在不同干涉量下最后一层的损伤如图11所示。随着干涉量的增加,锁定过程引起的CFRP材料损失范围逐渐增加,尤其基体损伤和层合板压缩分层损伤的范围随着干涉量的增加在增大。在干涉量为I0时,抽芯铆钉安装对孔壁没有损伤,损伤全部来自于“鼓包”的形成对最后一层的挤压,在I0干涉量下,CFRP没有出现基体拉伸和基体压缩损伤,出现了少量的拉伸分层和压缩分层损伤。通过对比锁定过程前后CFRP材料的损伤位置和损伤尺寸,可以看出锁定过程对CFRP材料损伤形成的影响较小,并且受影响区域仅位于接头孔出口附近。

    图  11  CFRP/Al 抽芯铆钉锁紧后CFRP最后一层的损伤
    Figure  11.  Damage of the last layer of CFRP after the CFRP/Al core-pulling rivet is locked

    建立了一套完整的抽芯铆钉安装和锁紧的模型,研究了碳纤维增强树脂复合材料(CFRP)/Al单剪铆接中抽芯铆钉的安装干涉量及安装速度对复合材料损伤的影响,基于连续损伤力学和扩展的三维破坏准则分析了材料的损伤。结论如下:

    (1) 在干涉配合下,抽芯铆钉的安装阻力随着安装速度的增加而减小。CFRP的损伤随着安装速度的增加而增加,在高干涉量I3下最明显;

    (2) 在低干涉量下,随着安装速度的增加,平均安装阻力的差异较小。而在高干涉量下,平均安装阻力的差异较大;

    (3) 抽芯铆钉的安装过程中,损伤主要发生在CFRP层压板的首层和尾层,层间产生了大量的损伤,随着干涉量的增加也发生了大量的层内损伤裂纹,内部开始萌生基体裂纹,同时基体裂纹由孔壁附近向层压板内部更深处扩展。抽芯铆钉的锁定对CFRP损伤全部来自于“鼓包”的形成对最后一层的挤压,在低干涉量下仅产生了对CFRP的拉伸分层损伤。通过对比锁定过程前后CFRP材料的损伤位置和损伤尺寸,可以看出锁定过程对CFRP材料新损伤形成的影响较小,并且受影响区域仅位于接头孔出口附近。

  • 图  1   力学性能3种退化方式

    Figure  1.   Three degradation modes of mechanical properties

    图  2   碳纤维增强树脂复合材料(CFRP)铆接模拟实施流程

    Figure  2.   Simulation implementation process of carbon fiber reinforced plastics (CFRP) riveting

    σ—Stress vector; C—Stiffness matrix of undamaged orthotropic material; ε—Strain vector

    图  3   CFRP/Al单钉单剪铆接几何参数

    Figure  3.   Geometrical parameters of CFRP/AL single nail single shear riveting

    图  4   网格划分及边界条件

    Figure  4.   Grid division and boundary conditions

    Ux, Uy, Uz—Displacement along the X-axis, Y-axis, Z-axis

    图  5   抽芯铆钉的铆接过程

    Figure  5.   Riveting process of the blind rivet

    图  6   CFRP/Al单钉单剪铆接在干涉量I1I2I3下的位移-载荷曲线

    Figure  6.   Displacement-load curves of CFRP/Al single-pin single-shear riveting with interference quantities I1, I2 and I3

    图  7   过盈配合接头的载荷-位移曲线的典型模型

    Figure  7.   Typical model of load-displacement curve of interference fit joint

    图  8   干涉量I1I2I3下1 m/s、5 m/s和10 m/s所对应的CFRP的纤维拉伸损伤、基体拉伸损伤、层合板拉伸损伤和层合板压缩损伤

    Figure  8.   Fiber tensile damage, matrix tensile damage, laminate tensile damage and laminate compression damage of CFRP corresponding to 1 m/s, 5 m/s and 10 m/s under interference quantities I1, I2 and I3

    SDV1—Fiber tensile damage; SDV3—Matrix tensile damage; SDV5—Tensile delamination damage; SDV6—Compressive delamination damage

    图  9   干涉量I1I2I3下的纤维拉伸损伤、基体拉伸损伤、拉伸分层损伤和压缩分层损伤

    Figure  9.   Fiber tensile damage, matrix tensile damage, tensile delamination damage and compressive delamination damage under interference quantities I1, I2 and I3

    图  10   CFRP/Al抽芯铆钉锁定的应力云图

    Figure  10.   Stress nebulogram of a locked CFRP/Al rivet

    S—Von mises stress (MPa)

    图  11   CFRP/Al 抽芯铆钉锁紧后CFRP最后一层的损伤

    Figure  11.   Damage of the last layer of CFRP after the CFRP/Al core-pulling rivet is locked

    表  1   T700/M21复合材料属性参数

    Table  1   Property parameters of T700/M21 composites

    E11/GPa E22/GPa E33/GPa G12/GPa G13/GPa G23/GPa v12 v13 v23
    117 7.7 7.7 4.8 4.8 2.8 0.34 0.34 0.4
    Notes: E11, E22, E33—Elastic modulus in the 1 (axial), 2 (transverse), and 3 (out-of-plane) directions; G12, G13, G23—Shear modulus for the 12, 13, and 23 planes; v12, v13, v23—Principal Poisson's ratios for the 12, 13, and 23 planes.
    下载: 导出CSV

    表  2   T700/M21复合材料层间强度

    Table  2   Interlaminar strength of T700/M21 composites

    XT/MPa XC/MPa YT/MPa YC/MPa ZT/MPa ZC/MPa S12/MPa S13/MPa S23/MPa
    2200 1700 50 300 50 300 95 95 95
    Notes: XT, XC—Tensile and compressive strengths in the fiber direction; YT, YC—Tensile and compressive strengths in the matrix direction; ZT, ZC—Tensile and compressive strengths in the out-of-plane direction; S12, S13, S23—Shear strengths for the 12, 13, and 23 planes.
    下载: 导出CSV

    表  3   抽芯铆钉各零件和铝合金的属性参数

    Table  3   Property parameters of each part of blind rivet and aluminum alloy

    Name Type Young's modulus/GPa Poisson's ratio Density/(g·cm−3)
    Sleeve CPTi 100 0.34 4.51
    Stem Ti3Al8V6Cr4Mo4Zr 107.8 0.32 4.82
    Lock collar A286 2.01 0.305 7.49
    Shifting anvil 45A 206 0.269 7.85
    Al alloy Al 72 0.33 2.70
    下载: 导出CSV

    表  4   CFRP/Al抽芯铆钉相对干涉量I大小

    Table  4   Relative interference magnitude I of the CFRP/Al blind rivet

    Number Interference amount/% Sleeve diameter/mm Aperture/mm
    I0 0 4.80 4.80
    I1 0.42 4.80 4.78
    I2 0.84 4.80 4.76
    I3 1.30 4.80 4.74
    下载: 导出CSV
  • [1]

    LIU Y, XIE J Z, TAFSIROJJAMAN T, et al. CFRP lamella stay-cable and its force measurement based on microwave radar[J]. Case Studies in Construction Materials, 2022, 16: e00824. DOI: 10.1016/j.cscm.2021.e00824

    [2]

    DU Y, MA Y, SUN W B, et al. Effect of hygrothermal aging on moisture diffusion and tensile behavior of CFRP composite laminates[J]. Chinese Journal of Aeronautics, 2023, 36(3): 382-392. DOI: 10.1016/j.cja.2022.11.022

    [3] 王跃全, 童明波, 朱书华. 三维复合材料层合板渐进损伤非线性分析模型[J]. 复合材料学报, 2009, 26(5): 159-166. DOI: 10.3321/j.issn:1000-3851.2009.05.026

    WANG Yuequan, TONG Mingbo, ZHU Shuhua. Progressive damage nonlinear analysis model of three-dimensional composite laminates[J]. Acta Materiae Compositae Sinica, 2009, 26(5): 159-166(in Chinese). DOI: 10.3321/j.issn:1000-3851.2009.05.026

    [4]

    DENG B Y, LI L Z, TAN D, et al. Sustainable and cost-effective ultra-lightweight engineered cementitious composite: Design and material characterization[J]. Cement and Concrete Composites, 2023, 136: 104895. DOI: 10.1016/j.cemconcomp.2022.104895

    [5]

    WANG H, LI H, FAN J, et al. Numerical analysis of dynamic response: Fatigue behaviour analysis of Al alloy-CFRP riveted single-shear lap joints[J]. International Journal of Fatigue, 2023, 170: 107515. DOI: 10.1016/j.ijfatigue.2023.107515

    [6]

    KAPPE F, ZIRNGIBL C, SCHLEICH B, et al. Determining the influence of different process parameters on the versatile self-piercing riveting process using numerical methods[J]. Journal of Manufacturing Processes, 2022, 84: 1438-1448. DOI: 10.1016/j.jmapro.2022.11.019

    [7] 周金秋, 王振军, 杨思远, 等. 连续石墨纤维增强铝基复合材料准静态拉伸损伤演化与断裂力学行为[J]. 复合材料学报, 2020, 37(4): 907-918.

    ZHOU Jinqiu, WANG Zhenjun, YANG Siyuan, et al. Quasi-static tensile damage evolution and fracture mechanics behavior of continuous graphite fiber reinforced aluminum matrix composites[J]. Acta Materiae Compositae Sinica, 2019, 37(4): 907-918(in Chinese).

    [8]

    SCATTINA A, RONCATO D, BELINGARDI G, et al. Investigation of creep phenomenon on composite material for bolt connections[J]. Composite structures, 2015, 134: 378-383. DOI: 10.1016/j.compstruct.2015.08.061

    [9]

    KARATHANASOPOULOS N, MOHR D. Strength and failure of self-piercing riveted aluminum and steel sheet joints: Multi-axial experiments and modeling[J]. Journal of Advanced Joining Processes, 2022, 5: 100107. DOI: 10.1016/j.jajp.2022.100107

    [10]

    JIANG H, GAO S, LI G, et al. Structural design of half hollow rivet for electromagnetic self-piercing riveting process of dissimilar materials[J]. Materials & Design, 2019, 183: 108141.

    [11]

    HUANG Z C, ZHANG Y K, LIN Y C, et al. Physical property and failure mechanism of self-piercing riveting joints between foam metal sandwich composite aluminum plate and aluminum alloy[J]. Journal of Materials Research and Technology, 2022, 17: 139-149. DOI: 10.1016/j.jmrt.2021.12.132

    [12]

    CAO Z Q, ZUO Y J. Electromagnetic riveting technique and its applications[J]. Chinese Journal of Aeronautics, 2020, 33(1): 5-15. DOI: 10.1016/j.cja.2018.12.023

    [13]

    MUCHA J, KAŠČÁK L, SPIŠÁK E. The experimental analysis of forming and strength of clinch riveting sheet metal joint made of different materials[J]. Advances in Mechanical Engineering, 2013, 5: 848973. DOI: 10.1155/2013/848973

    [14]

    RUSIA A, WEIHE S. Development of an end-to-end simulation process chain for prediction of self-piercing riveting joint geometry and strength[J]. Journal of Manufacturing Processes, 2020, 57: 519-532. DOI: 10.1016/j.jmapro.2020.07.004

    [15]

    ABE Y, KATO T, MORI K. Self-piercing riveting of high tensile strength steel and aluminium alloy sheets using conventional rivet and die[J]. Journal of Materials Processing Technology, 2009, 209(8): 3914-3922. DOI: 10.1016/j.jmatprotec.2008.09.007

    [16]

    ZHAO J Q, ZHOU H J, SUN X Y, et al. Evaluation on compressive properties of composite laminates with a hole reinforced by metal plate[J]. Composite Structures, 2021, 258: 113423. DOI: 10.1016/j.compstruct.2020.113423

    [17]

    KHASHABA U A, SEBAEY T A, ALNEFAIE K A. Failure and reliability analysis of pinned-joints composite laminates: Effects of stacking sequences[J]. Composites Part B: Engineering, 2013, 45(1): 1694-1703. DOI: 10.1016/j.compositesb.2012.09.066

    [18]

    LIU Y, LI M, LU X, et al. Pull-out performance and optimization of a novel interference-fit rivet for composite joints[J]. Composite Structures, 2021, 269: 114041. DOI: 10.1016/j.compstruct.2021.114041

    [19]

    CAO Y J, CAO Z Q, ZUO Y J, et al. Numerical and experimental investigation of fitting tolerance effects on damage and failure of CFRP/Ti double-lap single-bolt joints[J]. Aerospace Science and Technology, 2018, 78: 461-470. DOI: 10.1016/j.ast.2018.04.042

    [20]

    ZUO Y J, YUE T, JIANG R, et al. Bolt insertion damage and mechanical behaviors investigation of CFRP/CFRP interference fit bolted joints[J]. Chinese Journal of Aeronautics, 2022, 35(9): 354-365. DOI: 10.1016/j.cja.2022.01.027

    [21]

    HU J, ZHANG K, YANG Q, et al. Fretting behaviors of interface between CFRP and coated titanium alloy in composite interference-fit joints under service condition[J]. Materials & design, 2017, 134: 91-102.

    [22]

    LI J, ZHANG K, LI Y, et al. Influence of interference-fit size on bearing fatigue response of single-lap carbon fiber reinforced polymer/Ti alloy bolted joints[J]. Tribology International, 2016, 93: 151-162. DOI: 10.1016/j.triboint.2015.08.044

    [23]

    LI Q M, MINES R A W, BIRCH R S. Static and dynamic behaviour of composite riveted joints in tension[J]. International Journal of Mechanical Sciences, 2001, 43(7): 1591-1610. DOI: 10.1016/S0020-7403(00)00099-0

    [24]

    ZHAO J, PENG Y, SU W, et al. Finite element analysis of the shear capacity of stainless-steel blind-rivet con nections[J]. Journal of Constructional Steel Research, 2021, 179: 106558. DOI: 10.1016/j.jcsr.2021.106558

    [25]

    XING P, YU H. Tension-tension fatigue behavior and its modelling for life estimation of rivet-bonding CFRP/CFRP single lap joints[J]. Composite Structures, 2021, 275: 114367. DOI: 10.1016/j.compstruct.2021.114367

    [26]

    KELLY B, COSTELLO C. FEA modelling of setting and mechanical testing of aluminum blind rivets[J]. Journal of Materials Processing Technology, 2004, 153: 74-79.

    [27]

    ZUO Y, CAO Z, ZHENG G, et al. Damage behavior investigation of CFRP/Ti bolted joint during interference fit bolt dynamic installation progress[J]. Engineering Failure Analysis, 2020, 111: 104454. DOI: 10.1016/j.engfailanal.2020.104454

    [28]

    NAKAYASU H, MAEKAWA Z. A comparative study of failure criteria in probabilistic fields and stochastic failure envelopes of composite materials[J]. Reliability Engineering & System Safety, 1997, 56(3): 209-220.

    [29] 彭亮, 黄文博, 毛伟, 等. 碳纤维增强树脂基层板应变率相关损伤数值研究[J]. 北京航空航天大学学报, 2017, 43(4): 694-700.

    PENG Liang, HUANG Wenbo, MAO Wei, et al, Numerical study on strain rate-related damage of carbon fiber reinforced resin substrate[J]. Journal of Beijing University of Aeronautics and Astronautics, 2017, 43(4): 694-700(in Chinese).

    [30]

    NAIK N K, SHRIRAO P. Composite structures under ballistic impact[J]. Composite Structures, 2004, 66(1-4): 579-590. DOI: 10.1016/j.compstruct.2004.05.006

    [31]

    ASTM. Standard test method for bearing response of polymer-matrix composite laminates: ASTM D5961/D5961M—05[S]. West Conshohocken: ASTM International, 2005.

    [32]

    BERTHE J, BRIEU M, DELETOMBE E, et al. Temperature effects on the time dependent viscoelastic behaviour of carbon/epoxy composite materials: Application to T700GC/M21[J]. Materials & Design, 2014, 62: 241-246.

    [33]

    LI S, ZHANG S, LI H, et al. Numerical and experimental investigation of fitting tolerance effects on bearing strength of CFRP/Al single-lap blind riveted joints[J]. Composite Structures, 2022, 281: 115022. DOI: 10.1016/j.compstruct.2021.115022

    [34] 左杨杰. 复合材料/钛合金结构干涉配合螺栓连接动静态损伤及失效研究[D]. 西安: 西北工业大学, 2018.

    ZUO Yangjie. Study on dynamic and static damage and failure of bolt connection of composite/titanium alloy structure with interference fit[D]. Xi'an: Northwestern Polytechnical University, 2018(in Chinese).

  • 目的 

    随着科技的进步,复合材料已成为航空航天轻量化首选材料。然而,复合材料构件之间的连接一直是一个具有挑战性的问题。在复合材料构件连接中,抽芯铆钉连接技术备受关注。抽芯铆钉连接是一种利用铆钉中空的多种功能来取代传统铆钉的连接方法,适用于复合材料与金属材料之间的连接。本文采用了完整的抽芯铆钉模型,基于连续损伤力学、扩展的三维破坏准则和应变率效应,通过有限元和Abaqus二次开发相结合的方法,探究CFRP/Al单剪铆接中抽芯铆钉与连接孔的干涉量和铆钉的安装速度对复合材料损伤的影响。

    方法 

    本文采用完整的抽芯铆钉模型,结合连续损伤力学、扩展的三维破坏准则和应变率效应,通过有限元和Abaqus二次开发相结合的方法进行研究。完整的抽芯铆钉模型考虑了材料的复合性质,包括复合材料和金属材料的力学特性和界面效应。该模型可以准确描述抽芯铆钉连接在受力过程中出现的损伤、破坏和失效行为;连续损伤力学是一种用于描述材料损伤和断裂行为的数学模型。在研究中,该模型被应用于抽芯铆钉连接的损伤和失效分析,以预测连接在受力过程中的损伤演化和破坏模式;扩展的三维破坏准则是一种用于描述材料破坏准则的模型。在研究中,该准则被应用于抽芯铆钉连接的破坏分析,以预测连接的破坏模式和承载能力;应变率效应考虑了材料的变形速率对其力学行为的影响。在抽芯铆钉连接研究中,应变率效应被考虑在内,以更准确地描述连接在动态加载下的响应。通过有限元和Abaqus二次开发相结合的方法,建立抽芯铆钉连接的数值模型,并通过模拟加载过程,获得连接的应力分布、变形行为和破坏模式。利用这种方法,可以评估连接的承载能力、耐久性和安全性,并为设计和优化提供依据。

    结果 

    模拟结果显示安装速度对抽芯铆钉连接的影响是显著的。安装阻力表现出两个典型阶段,分别由铝合金板、CFRP板与抽芯铆钉之间的摩擦所产生。在抽芯铆钉安装在干涉配合情况下,较高的安装速度可以有效降低安装阻力。但是,过高的安装速度可能会增加孔壁的损伤,特别是在高干涉量下,这一影响更加明显。此外,抽芯铆钉铆接接头的失效模式主要以CFRP损伤为主,并且受到干涉尺寸的显著影响。这意味着抽芯铆钉的安装过程会导致CFRP板的损伤,并且干涉配合尺寸的改变会对接头的失效模式产生明显影响。

    结论 

    本文建立了一套完整的抽芯铆钉安装和锁紧的模型,研究了CFRP/Al单剪铆接中抽芯铆钉的安装干涉量及安装速度对复合材料损伤的影响,基于连续损伤力学和扩展的三维破坏准则分析了材料的损伤。获得如下

    结论 

    在干涉配下,抽芯铆钉的安装阻力随着安装速度的增加而减小。CFRP的损伤随着安装速度的增加而增加,在高干涉量下最明显。在低干涉量下,随着安装速度的增加,平均安装阻力的差异较小。而在高干涉量下,平均安装阻力的差异较大。抽芯铆钉的安装过程中,损伤主要发生在CFRP层压板的首层和尾层,层间产生了大量的损伤,随着干涉量的增加也发生了大量的层内损伤裂纹,内部开始萌生基体裂纹,同时基体裂纹由孔壁附近向层压板内部更深处扩展。抽芯铆钉的锁定对CFRP损伤全部来自于“鼓包”的形成对最后一层的挤压,在低干涉量下仅产生了对CFRP的拉伸分层损伤。通过对比锁定过程前后CFRP材料的损伤位置和损伤尺寸,可以看出锁定过程对CFRP材料新损伤形成的影响较小,并且受影响区域仅位于接头孔出口附近。

  • 复合材料具有高强度、轻质和良好的耐腐蚀性能,因此成为了航天航空结构的首选材料。然而,复合材料构件之间的连接却一直是一个具有挑战性的问题。其中抽芯铆钉连接是一种备受关注的技术。相比于传统的连接方法,抽芯铆钉连接具有强度高、操作简单和轻质等优点。然而,抽芯铆钉连接仍然存在一些挑战,铆钉的安装会损伤复合材料的纤维结构,并降低连接的强度。

    本文采用了完整的抽芯铆钉模型,如图1(a)所示。并基于连续损伤力学和扩展的三维破坏准则,通过有限元和Abaqus二次开发相结合的方法,研究了CFRP/Al单剪铆接中抽芯铆钉与连接孔的干涉量和铆钉的安装速度对复合材料损伤的影响。从模拟结果图1(b)来看,抽芯铆钉安装在干涉配合情况下,安装速度越大对于降低安装阻力十分有利。并且抽芯铆钉铆接对CFRP最后一层损伤微小可忽略不计,接头的失效模式以CFRP损伤为主,并受干涉尺寸的影响显著。

    抽芯铆钉铆接后应力云图(a)和CFRP/AL单钉单剪铆接在干涉量为I2下的位移载荷曲线(b)

图(11)  /  表(4)
计量
  • 文章访问数:  440
  • HTML全文浏览量:  243
  • PDF下载量:  35
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-07-31
  • 修回日期:  2023-10-10
  • 录用日期:  2023-10-19
  • 网络出版日期:  2023-10-24
  • 刊出日期:  2024-05-31

目录

/

返回文章
返回