端锚CFRP-混凝土界面剥离全过程解析模型

董 坤，荆范华，郝爱江，李 鹏，郝建文

(1.中国海洋大学 工程学院，山东 青岛，266100；2.滨州市公路勘察设计院，山东 滨州，256600；3.山东大学 土建与水利学院，济南 250000)

碳纤维增强聚合物(CFRP)是一种轻质高强的材料，外贴CFRP材料加固和修复混凝土结构已成为目前加固工程中一种普遍的选择。CFRP加固技术不但能够提高结构的极限承载力，同时还能提供较好的抗腐蚀效果[1-2]。在实际工程中，CFRP-混凝土黏结界面的端部剥离是CFRP加固混凝土结构的典型破坏模式[3-4]，不但会导致CFRP强度的有效利用率降低，并且将导致加固构件脆性破坏。因此，ACI委员会建议，可采用附加端部锚固的措施来增加极限荷载和避免黏结界面端部的剥离破坏[5]。

针对端部锚固措施，国内外学者开发了一系列不同的锚固装置和方法，试验表明锚固手段可以有效地控制端部剥离[6]。覃银辉等[7]基于一种自锁式的端部锚固装置，进行了CFRP-混凝土黏结界面剪切试验，试验中锚固效果突出，施加端部锚固后的破坏模式为CFRP片材拉断，黏结界面承载力得到显著提升。Barris等[8]采用钢压板加螺栓的方式对CFRP端部进行锚固，试验表明，通过对螺栓施加足够的扭矩也能够实现端部完全锚固的效果。另外，周英武等[9]、卓静等[10]、Zhou等[11]还研究了不同形式的端部锚固装置对CFRP加固梁抗弯性能的提升，发现通过合理的端部锚固，能够极大的提高构件的受弯承载力和CFRP强度利用效率。在理论研究方面，李春良等[12]建立了端部锚固CFRP加固结构的界面黏结应力解析模型，但该模型忽略黏结界面产生的滑移，会对黏结界面的黏结应力水平造成偏高的预测。Zhang等[13]和Sturm等[14]基于双线性界面黏结滑移本构，采用分段积分求解方法分别模拟了端部单个和多个纤维束锚固黏结界面的受力全过程，获得了黏结界面加载端荷载-滑移关系及界面黏结性能分布模型。Chen等[15]则采用三线性分段界面黏结-滑移本构，给出了加固界面极限承载力预测模型。上述理论研究中均采用线性分段的界面黏结-滑移本构形式，无法连续的体现黏结界面的非线性强化及软化行为，且较少涉及对端部锚固下有效黏结长度和最大界面黏结力的探讨。

本文基于端部锚固CFRP-混凝土黏结界面的剪切受力模型，引入了双参数指数型界面黏结-滑移本构，建立了表征端部锚固CFRP-混凝土黏结界面剥离全过程的解析模型，该解析模型得到了试验结果的良好验证。利用解析模型，建立了最大界面黏结力、黏结界面剥离承载力和有效黏结长度的计算方法，并对不同黏结长度界面的剥离全过程进行了分析。

1 黏结界面剥离行为解析模型

1.1 黏结界面微分平衡方程

在黏结界面剥离行为的解析模型建立过程中，做出如下基本假设：1) 黏结界面仅承受切向黏结应力，不承受法向应力；2) 不考虑黏结层厚度，黏结层的性质在界面黏结-滑移本构中予以体现；3) CFRP、混凝土材料均为弹性体，非线性力学特征仅存在于黏结界面；4) CFRP片材所受正应力沿厚度方向均匀分布，不考虑宽度方向的应力变化；5) 锚固端为完全锚固，在受力过程中锚固位置处的CFRP片材不产生滑动。

图1给出了端部锚固CFRP-混凝土黏结界面的单面剪切受力示意图。

图1 端部锚固CFRP-混凝土的单面剪切受力示意Fig.1 Schematic diagram of single-lap shear test of CFRP-to-concrete interface with end anchorage

根据图1(b)中微元体的受力情况，可建立平衡方程：

(1)

式中：σf和σc分别为CFRP和混凝土的轴向应力,tf和tc分别为CFRP和混凝土的厚度，bf和bc分别为CFRP和混凝土的宽度,τ为界面黏结应力。

CFRP、混凝土以及黏结界面的物理关系为

(2)

式中：uf和uc分别为CFRP和混凝土的轴向变形量，εf和εc分别为CFRP和混凝土的轴向应变，Ef和Ec分别为CFRP和混凝土的弹性模量，s为黏结界面的相对滑移量。

联立式(1)、(2)可得

(3)

式中:ρ=Eftfbf/(Ectcbc)对于常见的CFRP加固混凝土构件，值一般小于0.01，对界面黏结行为的影响较小，在本文后续推导过程中暂忽略不计。

联立式(1)～(3)可得黏结界面微分平衡方程：

(4)

1.2 黏结微分平衡方程的求解

界面黏结-滑移本构是分析CFRP-混凝土黏结界面剥离过程的关键。Dai等[16]提出了如式(5)所示的双参数指数型界面黏结-滑移本构，该本构的待定参数少，表达形式为一条光滑的曲线，能够简洁有效地反映CFRP-混凝土的界面黏结-滑移非线性行为。相比于双线性或三线性本构模型，基于指数型本构所推导得到的荷载-滑移关系式和各物理量分布表达式相对统一，不分段的单一表达式便可反映黏结界面剥离全过程的行为。目前，该本构已广泛应用于CFRP-混凝土界面力学行为的解析分析中[17-18]。

τ(s)=EftfA2B(1-e-Bs)e-Bs

(5)

式中A和B为黏结界面参数，可通过黏结界面剪切试验或简化计算方法[19]求得。

联立式(4)～(5)可得

(6)

由于锚固端CFRP应变暂时是未知的，假定其值为ε0，该值的求解方法将在后续给出。结合式(3)对图1(a)中端部锚固CFRP-混凝土黏结界面进行分析，可得边界条件：

(7)

对式(6)进行一次积分可得

(8)

式(8)中C1为常数项。结合边界条件式(7)中的第一项可得

(9)

(10)

对式(10)进行积分可得

(11)

式(11)中C2为常数项，对上式化简可得变量y的表达式

(12)

根据边界条件知y(0)=1，将其代入式(12)可求解得

(13)

(14)

最后，结合y=e-Bs可得界面滑移分布的计算公式：

(15)

联立式(3)和式(15)，可得CFRP应变分布的计算公式：

(16)

同时结合式(2)可得CFRP的正应力分布计算公式：

(17)

联立式(1)和式(17)可得界面黏结应力分布的计算公式：

(18)

将边界条件式(7)中第二项代入式(3)及式(16)可得

(19)

结合φ=ε0/A，式(19)给出了外荷载P与锚固端CFRP应变的一一对应关系。此时，在给定的荷载下可以求得唯一的锚固端CFRP应变ε0，进而可根据式(15)～(19)求得其他物理量表达式。

将x=L代入式(15)可得黏结界面加载端的滑移计算公式：

(20)

联立式(19)、(20)即可求得黏结界面加载端的荷载-滑移响应曲线。

2 解析模型的验证

为验证解析模型的正确性和适用性，设计并进行了端部锚固CFRP-混凝土黏结界面的单剪试验。试验中试件信息见表1，包括3组不同黏结长度的端部锚固CFRP布加固混凝土试件，以及1组纯外贴CFRP布加固混凝土试件，每组3个试件。如图2所示，本文试验使用5 mm厚的45号钢板设计了新型自锁锚固装置，其工作原理是将CFRP布端部反向覆盖开洞钢板，然后将咬合钢板从上部压紧并采用8.8级M8螺栓固定，利用摩擦力和机械咬合力将CFRP布可靠地固定于锚固装置上。该装置满足JGJ 145—2013《混凝土结构后锚固技术规程》第6.1.1～6.1.3条的机械锚固强度验算要求。被加固试件混凝土弹性模量为25.5 GPa，混凝土28 d立方体抗压强度为37.6 MPa；实测涂胶固化后的CFRP布弹性模量为220 GPa，抗拉强度为2 870 MPa。试验获取了CFRP应变数据和界面荷载-滑移曲线，并与解析模型进行了对比。

表1 试件信息表Tab.1 Information of specimens

图2 新型自锁锚固装置及端锚CFRP加固试件示意Fig.2 Schematic diagram of new type self-locking anchorage device and CFRP strengthened specimen with end anchorage

图3给出了解析模型与试验结果的对比，其中将3个EB-200试件的应变数据按照文献[20]中的积分方法换算，可得距离加载端50 mm处的界面黏结应力与界面滑移数据，通过对试验数据拟合可得黏结-滑移本构公式(5)中的界面参数A=0.007 5，B=12，见图3(a)。根据图3荷载-滑移曲线及CFRP应变分布曲线的对比可知，解析模型与试验结果吻合较好。图3(b)部分试件荷载-滑移曲线的水平段数值和最终的CFRP拉断荷载数值均略低于理论值(13.5 kN和23.96 kN)，推测是由于试验过程中碳纤维布宽度的裁剪误差和胶黏剂的涂刷不均匀导致。

图3 解析模型与本文试验结果对比Fig.3 Comparisons between analytical and experimental results

3 界面黏结特征

3.1 界面平衡方程的求解

根据图1及前文推导的解析模型，可以得到整个加固界面的外荷载由端部锚固力和界面黏结力共同来承担，数学表达式为

P=Pa+Pb

(21)

根据锚固点受力平衡关系和式(17)可知，端部锚固力Pa为

Pa=bftfσf(0)=EfbftfAφ

(22)

结合式(19)、(21)、(22)可得到界面黏结力Pb为

(23)

图4以试验试件EA-150的参数为例，给出了加载过程中上述各部分的变化关系。

图4 P、Pa、Pb与滑移对应关系Fig.4 Relation between P, Pa, Pb, and slip

对式(23)求导取极值，可得到给定黏结长度L后，最大界面黏结力的表达式为

(24)

(25)

将φ0依次带入式(19)和式(22)即可求得对应最大界面黏结力Pb,0的外荷载P0和端部锚固力Pa,0。图5给出了不同黏结长度时，Pb,0及其对应的P0、Pa,0的变化曲线。图中数据均作无量纲处理：纵坐标为各荷载(Pb,0、Pa,0、P0)与EfbftfA的比值，横坐标为界面参数A、B与黏结长度L的乘积。

图5 Pb,0、Pa,0、P0与黏结长度对应关系Fig.5 Relation between Pb,0, Pa,0, P0, and bond length

由图5可知，随着黏结长度不断增加，Pb,0与对应的Pa,0分别呈现单调递增和单调递减的趋势，并分别趋近于1和0；Pb,0对应的P0呈现先递减后小幅度递增的趋势，并最终趋近于1。当黏结长度较大时，P0基本都由黏结力来承担，而端锚装置几乎不提供锚固力，此时界面状态与纯外贴情况几乎完全一致；反之，当黏结长度较小时，P0大部分都由锚固力承担，此时的界面状态与有端锚无黏贴情况几乎一致。

3.2 有效黏结长度

根据图5中3条曲线数据，图6给出了最大界面黏结力Pb,0及对应端部锚固力Pa,0所承担荷载的比例与黏结长度的关系。令Pb,0承担荷载的比例为α，对应Pa,0承担荷载的比例为β，通过对图6的曲线进行拟合，给出β的简化表达式

图6 各部分荷载比例变化曲线Fig.6 Load proportion curves of each component

β=1-α=

(26)

类比纯外贴CFRP加固情况，随着黏结长度增加，当最大界面黏结力Pb,0占外荷载的96%以上时，对应的黏结长度即为黏结界面的有效黏结长度[21]。对式(26)第二段求反函数，即可得到端部锚固CFRP-混凝土黏结界面的有效黏结长度计算公式

(27)

文献[19]中给出了考虑自由端滑移时，纯外贴CFRP-混凝土界面的最大黏结力承担荷载的比例α为

(28)

同样，对式(28)第二段求反函数，可得纯外贴CFRP-混凝土黏结界面的有效黏结长度计算公式

(29)

图7给出了纯外贴和端部锚固黏结界面的有效黏结长度的对比。由图7可知，相比于纯外贴黏结界面，在同一承担荷载的比例α下，端部锚固黏结界面所需的有效黏结长度有所增加，两者差值约为1/AB；端部锚固CFRP-混凝土黏结界面的有效黏结长度至少应为7.2/AB。

图7 有效黏结长度与α对应关系Fig.7 Relation between effective bond length and α

3.3 界面剥离荷载

如图8所示，不同于双线性界面黏结-滑移本构，双参数指数型界面黏结-滑移本构仅存在黏结应力强化段和软化段，软化段后不存在剥离点(即τ=0点)，因此无法直接给出对应于加载点黏结应力为0时的剥离荷载，而双线性界面黏结-滑移本构存在剥离点(τ=0点)。借鉴有效黏结长度的取值方法，可在双参数指数型界面黏结-滑移本构曲线的软化段上取极限滑移sf=6/B，见图8(b)，结合式(5)计算此时本构曲线与x轴所围面积可占完整曲线所围面积(即断裂能Gf)的99.5%。因此，可认为当加载端滑移达到sf=6/B时，该点开始剥离，此时外荷载即为剥离荷载。

图8 黏结-滑移本构模型中剥离点的计算方法Fig.8 Calculation method of debonding point in bond-slip model

将s(L)= 6/B代入式(20)，并结合式(19)即可求得黏结界面剥离荷载为

Pdb=EfbftfAφdb+

(30)

(31)

图9给出了剥离荷载Pdb及其对应的界面黏结力Pb,db和端部锚固力Pa,db随黏结长度变化曲线，同时与对应最大界面黏结力的外荷载P0的对比。

图9 剥离荷载及其对应黏结力、锚固力与黏结长度的关系Fig.9 Relation between debonding load and corresponding bond force, anchoring force, and bond length

对于纯外贴界面，其剥离荷载随黏结长度的增加而增大，且趋近于其上限EfbftfA[21]。而由图9可看出，端部锚固黏结界面的剥离荷载随着黏结长度的增加而减小，且趋近于其下限EfbftfA。这是因为在黏结长度较短时，在CFRP发生剥离前，端部锚固发挥了较大的作用。当ABL<6时，Pa,db可承担一半的剥离荷载；而当ABL<2时，Pa,db几乎承担全部剥离荷载。通过与外荷载P0的对比还可以看出，当界面黏结应力积分最大，也即达到最大界面黏结力时，其对应的外荷载始终小于剥离荷载，且随着黏结长度的增加逐渐趋近于剥离荷载。

4 黏结界面剥离全过程行为

以本文试验的材料属性为计算参数，根据3.2节有效黏结长度的计算公式，取α=0.97得到该工况下的有效黏结长度为Leff=83 mm。分别选取50 mm(LLeff)两种黏结长度对端部锚固CFRP-混凝土黏结界面剥离全过程行为进行分析，同时，选取加载过程中达到的4个荷载时刻进行对比，分别为：1)加载端黏结应力最大时的荷载Pτ,max；2)界面黏结力最大时的荷载P0；3)界面剥离荷载Pdb；4)CFRP片材拉断时的荷载Pmax。图10给出了两种黏结长度下的加载端荷载-滑移曲线。

图10 两种黏结长度下荷载-滑移曲线Fig.10 Load-slip curves under different bond lengths

由图10可看出，当L>Leff时，荷载-滑移响应明显呈现3个阶段。加载初期，荷载-滑移响应呈现快速增长的趋势。随着外荷载的增加，荷载-滑移响应进入一个“平台”阶段，外荷载增长缓慢，整个界面处于剥离由加载端向锚固端发展的过程，“平台”长度取决于黏结长度的大小；加载后期，曲线凹凸性发生变化，存在一条渐近线并向其靠拢。而当L

图11给出了按照式(16)和式(18)计算出的两种黏结长度下CFRP应变和黏结界面黏结应力分布曲线。由图11(a)、(b)可知，相比于黏结长度LLeff时剥离荷载较小，黏结界面会更早剥离，且剥离时靠近锚固端位置的应变值处于较低水平；对于L>Leff结点，黏结界面黏结力最大时与界面剥离时应变分布及数值均相差不大，对于L

图11 两种黏结长度下CFRP应变和界面黏结应力分布曲线Fig.11 Distribution curves of CFRP strain and interfacial bond stress under different bond lengths

5 结 论

1) 基于双参数指数型界面黏结-滑移本构，采用解析方法推导了端部锚固下CFRP-混凝土黏结界面滑移、黏结应力、CFRP应力和应变分布以及荷载-滑移响应的表达式，并将所得解析公式与试验结果进行了对比，吻合情况较好。

2) 对端部锚固CFRP-混凝土界面黏结特征进行了分析，建立了端部锚固下最大界面黏结力、有效黏结长度和界面剥离荷载的计算模型。分析表明，端部锚固CFRP-混凝土界面的有效黏结长度至少应为7.2/AB，黏结界面的剥离荷载随黏结长度的增加而降低，且趋近于EfbftfA。

3) 对于黏结长度较大的端锚CFRP-混凝土界面，其界面黏结应力积分最大时的外荷载与剥离荷载相差不大，且剥离荷载值小于黏结长度较短黏结界面的情况；而对于黏结长度较短的CFRP-混凝土界面，端部锚固能够更早地参与承担荷载，阻止整个界面剥离失效，其剥离荷载更接近CFRP拉断荷载。