CFRP平-折-平胶接接头应力模型

2020-06-24刘志明刘文飞

哈尔滨工业大学学报 2020年7期

许昶，刘志明，刘文飞

(1.北京交通大学机械与电子控制工程学院，北京 100044; 2.台州学院航空工程学院，浙江台州 317000)

胶接连接以其较高的连接效率，较好的水密性和气密性及较优异的疲劳和耐腐蚀性能，广泛地应用于飞机、汽车、轨道交通等复合材料主承载结构中. 应力分析是进行结构分析的必要步骤. 接头应力分布可以借助有限元或者闭式模型来分析.

对于包含复杂几何和材料的结构，利用有限元方法更为合适，而对于较为简单的结构，利用闭式模型可以求得应力的精确解，相较于有限元方法避免了较为烦琐的建模过程.

单搭接和双搭接作为胶接连接常见的连接形式，因其连接结构几何简单而受到学者们的关注. 但是，由于搭接区同时存在剪滞效应、弯曲效应和端部效应，使得任意影响因素的微小变化都会导致胶层应力分布的改变，因此获取较为精确的应力分布的解析模型并不容易. Volkersen[1]首次在解析模型中考虑了被粘体的弹性效应，建立了剪滞模型. Goland等[2]在模型中首次引入了弯曲效应，建立了一维梁理论接头应力模型. G-R模型也被很多学者认为是现代胶接理论的起点，之后很多模型都是基于G-R模型进行改进的，对这些应力理论模型的综述可以参考文献[3-4]. 清华大学的赵波等对单搭接胶接接头应力的分析方法进行了较多的研究[5-6]. 在此之后，又有一些新的理论和方法被用于接头应力的分析. Yousefsani等[7]基于整体层化理论对复合材料单搭接胶接接头在拉伸和弯曲载荷下的应力分布进行了研究，该作者同时运用此模型分析了单搭接和双搭接接头沿厚度方向层间应力的分布[8]. Selahi等[9]则基于能量方法和铁木辛柯梁理论，给出了不同边界条件和载荷情况下的胶接接头解析模型. Paroissien等[10]基于线弹性一维梁模型构造了一个四节点BBe单元，用于分析搭接区胶层的应力分布. Paroissien等[11]进一步对比了一维杆和一维梁宏观单元(ME)闭式模型与有限元模型在分析胶层应力时的差异，同时探究了搭接长度的影响.

为了减缓单搭接接头偏心加载产生的弯曲效应，Kishore等[12]在研究中提出了一种平-折-平 (FJF)胶接接头形式，使得接头在拉伸加载过程中承受面内力，接头强度相较于单搭接接头提升了90%，在结构连接中可以用于替代单搭接接头. 接头应力分析是进行接头强度分析的前提，快速、准确地计算接头应力，需要构造具有足够计算精度的半解析模型. 本文以文献[10]中提出的方法为基础，构造了圆弧段BBe单元，并基于有限元方法，分析了碳纤维复合材料FJF胶接连接接头的胶层应力分布, 进一步探究了胶层宽度和厚度对接头胶层应力的影响.

1 FJF胶接连接接头半解析模型

1.1 FJF胶接连接接头

FJF胶接接头的示意图如图1所示. 其中，L′为圆角以后水平搭接段长度，b为试样宽度，e1为层合板厚度，e为胶层厚度，r1为较大的圆角半径，r2为较小的圆角半径，θ为折线段与水平段之间的夹角.

1.2 BBe单元的建立

1.2.1 基本假设

为简化模型，在建立BBe单元时作出如下假设.

1)搭接区域胶层的厚度恒定，胶层应力沿厚度方向没有变化.

2)层合板被认为是欧拉-伯努利层合梁，遵循经典层合板理论.

3)层合板和胶层连接界面不存在缺陷[13].

1.2.2 水平段BBe单元的建立

本文考虑对称铺层的复合材料层合板，从而消除了层合板的耦合效应，即耦合刚度为0. 层合梁的外力和位移关系为

(1)

(2)

式中：Ni、Mi分别表示层合板i(i=1,2)的轴力和弯矩，ui、wi分别表示层合板i(i=1,2)的纵向位移、挠度.Ai、Di分别表示层合板i(i=1,2)的面内刚度和弯曲刚度.

胶层被模拟为经典双参数弹性模型[2]，仅考虑其剥离应力和剪切应力，表达式为

(3)

(4)

式中：E为胶层剥离方向模量，G为胶层剪切模量，e为胶层厚度，ei表示层合板i(i=1,2)的厚度，θi表示层合板i(i=1,2)的转角.

对层合板和胶层进行受力分析，如图2所示. 由此可以得到平衡方程：

(5)

(6)

(7)

图2 平直段段层合板与胶层受力图

Fig.2 Free body diagram of laminated plates and the adhesive in straight section

式中：Vi为层合板i(i=1,2)的剪力，b表示层合板i(i=1,2)的宽度.

对式(1)～(7)进行化简可以分别得出胶层剥离应力和剪应力的通解形式：

T=K1+K2esx+K3e-sx,

(8)

S=K4etxcostx+K5etxsintx+

K6e-txcostx+K7e-txsintx.

(9)

用所求的应力通解反过来可以求解层合板位移，见式(10)～(15)，其中Δ表示单元长度. 进而由本构关系可以求解内力. 位移和内力的表达式中包含K1～K7、J1～J3、J5～J6这12个待定的系数.

(10)

(11)

(12)

(13)

(14)

(15)

式中：

图3给出了BBe单元的节点位移和节点力14]. 分别令x=0和x=Δ，由位移和内力表达式即可求得单元节点位移和节点力关于12个待定系数的表达式：

(16)

(17)

图3 平直段BBe单元节点位移和节点力

Fig.3 Node displacements and node forces of BBe element in straight section

进而水平段BBe单元的刚度矩阵可以由式(18)求得:

KBBe=NM-1.

(18)

1.2.3 圆弧段BBe单元的建立

此外，Warwick Mills还开发了Metal Flex Armor(MFA)柔性装甲和SoftPlateTM防弹衣，这两种层压柔性复合材料结合了紧密织造以及与NASA合作期间开发的涂层整理技术。其中，MFA由Twaron织物与刚性结构层压制成，可防碎冰锥、皮下注射针和刀具刺穿，提供与刚性防刺装甲相同的性能，但具有更高的灵活性和舒适性。SoftPlateTM与MFA结构相似，但主要目的是防弹。

Mohamad[15]在研究层合壳和层合板的振动时引申了对层合曲梁的分析. 层合曲梁的外力和位移关系为

(19)

(20)

对圆弧段层合板和胶层进行受力分析，如图4所示. 平衡方程为

(21)

(22)

(23)

用极坐标系定义，有αi=Riγ，则外力和位移关系可以转化为

(24)

(25)

平衡方程可以转化为

(26)

(27)

(28)

图4 圆弧段层合板与胶层受力图

Fig.4 Free body diagram of laminated plates and adhesive in curved section

(29)

式中，集合Y有如下表达方式：Y={y1,y2,y3,y4,y5,y6,y7,y8,y9,y10,y11,y12}′. 通过求解矩阵[A]的特征根和特征向量，即可求解出层合板位移的通解表达式，见式(30). 式中nr表示特征根为实根的个数，n11=nr+1. 当特征根为复数根时，将其分解为实部和虚部两部分，nim表示特征根为复数根的个数；COij(i=1,...,nim2;j=1,2,3,4)为特征根对应特征向量的元素. 经化简整理得到位移通解包含C1～C12这12个待定的系数.

(30)

用位移通解求解层合板内力，分别令γ=0 和γ=γ0，即可求得单元节点位移和节点力关于12个待定系数的表达式，进而由式(18)求得圆弧段BBe单元的刚度矩阵.

1.3 非搭接区梁单元和总体刚度矩阵

非搭接区的平衡方程为

(31)

(32)

(33)

求解非搭接区梁单元刚度矩阵参照求解BBe单元刚度矩阵的思路，求得的刚度矩阵为

(34)

式中li为非搭接区梁i的长度，Δi=AiDi.

对整体胶接接头进行网格离散，如图5所示(此处每段仅划分1个单元用于示意). 总体刚度矩阵基于经典有限元方法，由各个单元的刚度矩阵求得.

图5 FJF胶接接头网格离散

Fig.5 Mesh discretization of FJF adhesively bonded joint

2 FJF胶接连接接头有限元模型

胶接接头胶层应力分析采用三维有限元模型，如图6所示. 接头基本几何参数如表1所示.

层合板选用T700/5429碳纤维增强双马来酰亚胺树脂基复合材料，单向板的基本材料属性如下[16]：E1=133 GPa，E2=9.1 GPa，ν12=0.31，G12=5.67 GPa，G13=5.67 GPa，G23=3.5 GPa，XT=2 507 MPa，XC=1201 MPa，YT=61.8 MPa，YC=186 MPa，S12=84.8 MPa，S13=84.8 MPa，S23=41.6 MPa. 其中，E1、E2分别为单向板纵、横向的弹性模量；G12、G13、G23为单向板剪切模量；ν12为单向板纵向的泊松比；XT、XC分别为单向板纵向拉伸、压缩强度；YT，YC分别为单向板横向拉伸、压缩强度；S12、S13、S23为单向板剪切强度. 胶层材料采用J299高韧性双马结构胶膜，弹性模量为2.92 GPa，泊松比为 0.47.

模型中碳纤维层合板属性采用等效工程常数，由单层板的面内材料属性和铺层信息计算得出，单元类型为三维实体单元C3D8R. 胶层同样采用三维实体单元C3D8R，胶层单元与层合板单元通过共节点连接. 在接头搭接区两侧端部和圆角过渡区通过加密网格的方式以便更好地捕捉应力信息，考虑到胶层的厚度尺寸，搭接区两侧端部的网格大小设置为0.2 mm，且向远离端部长度和厚度方向采用渐变种子布局，搭接区端部局部网格加密如图7所示.

表1 FJF胶接接头基本几何参数

Tab 1 Basic geometric parameters of FJF adhesively bonded joint

参数数值参数数值参数数值L'/mm42.95b/mm16e/mm0.2e1/mm2r1/mm30θ/mm7.83

图6 CFRP FJF胶接连接接头三维有限元模型

Fig.6 3D Finite Element Model of CFRP FJF adhesively bonded joint

图7 局部网格加密图Fig.7 Local mesh refinement

3 结果与分析

3.1 网格收敛性分析

对在胶层厚度上划分网格数na取值2、4、8、16、32和64，在外载为1 kN的条件下获取胶层端部最大剥离应力和最大剪切应力值，结果如图8和图9所示.

图8 胶层厚度方向不同网格数量下接头端部最大剥离应力

Fig.8 Maximal peel stresses at the joint end with different mesh numbers along the adhesive thickness

图9 胶层厚度方向不同网格数量下接头端部最大剪切应力

Fig.9 Maximal shear stresses at the joint end with different mesh numbers along the adhesive thickness

由图8和图9可知，随着胶层厚度方向网格数量的不断增大，胶层端部最大剥离应力和最大剪切应力值趋于稳定，表明网格逐渐收敛.

3.2 不同模型胶层应力对比分析

采用收敛的网格划分胶层厚度，在后处理中获得胶层剥离应力和剪切应力沿搭接区长度方向的应力分布，并与半解析模型计算得到的胶层应力分布进行对比，如图10和图11所示. 由图10和图11可知，由半解析模型计算得到的胶层剥离应力和剪切应力分布与三维有限元模型得到的对应应力分布基本一致. 胶层最大剥离应力和剪切应力均在搭接区端部，且两端基本对称. 半解析模型胶层端部剥离应力值为6.46 MPa，相较于三维有限元模型胶层端部剥离应力值(6.83 MPa)误差绝对值为5.4%；半解析模型胶层端部剪切应力值为6.98 MPa，相较于三维有限元模型胶层端部剪切应力值(7.25 MPa)误差绝对值为3.7%.

图10 三维有限元模型与半解析模型胶层剥离应力分布对比

Fig.10 Comparison of adhesive peel stress distributions between 3D Finite Element Model and semi-analytical model

图11 三维有限元模型与半解析模型胶层剪切应力分布对比

Fig.11 Comparison of adhesive shear stress distributions between 3D Finite Element Model and semi-analytical model

4 胶层宽度和厚度的影响

研究表明，胶层厚度[17]和胶层宽度[18]对复合材料胶接接头搭接区胶层应力均有影响. 由此，基于半解析模型对不同胶层宽度和厚度组合条件下胶层应力分布进行对比分析，组合如表2所示. 对比第1、2和3组可以得出在胶层厚度一定的情况下，胶层宽度对搭接区胶层应力的影响；对比第3、4和5组可以得出在胶层宽度一定的情况下，胶层厚度对搭接区胶层应力的影响.

表2 不同胶层宽度和厚度组合

Tab.2 Different width and thickness combinations of the adhesive