1. 引言
国外学者Viest I.M.[1]通过栓钉连接件的推出试验,研究了栓钉的材料属性如何决定其抗剪承载力,并指出了栓钉的极限抗拉强度及混凝土强度等级也对其有重要影响。Ollgaard[2]经过对常规栓钉连接件的测试,得出了荷载–滑移曲线公式及抗剪承载力计算方法,并分析了影响栓钉承载力的各种参数。Spencer和Neillel[3]为通过循环荷载下的推出试验,研究了栓钉连接件的抗剪性能。Oehlers和Johnson[4]发现,栓钉连接件在单调荷载下可能剪断,且在达到极限抗剪强度后荷载会在一定滑移范围内保持稳定。他们还指出,在计算栓钉抗剪承载力时,需要考虑试验与实际状态的差异及连接件对单栓钉的整体影响。Bursi和Gramola[5]通过循环荷载试验,探究了地震荷载下栓钉连接件的强度和延性。Shim[6]研究了大直径栓钉与普通栓钉在抗剪承载力上的差异,发现大直径栓钉虽可用常规公式估算,但对布置有更严格要求。Pallares和Hajjarl[7]分析了三百多个推出试验,为栓钉连接件设计提供了规范参考。Lam D[8]通过考虑材料非线性行为,使用ABAQUS软件模拟得出了预制空心板组合梁中栓钉连接件的抗剪承载力。
国内学者聂建国[9]研究了高强混凝土中钢–混组合梁的栓钉抗剪性能,提出了一种适用于此环境的栓钉抗剪承载力公式。汪炳等[10]利用ABAQUS建立了焊钉连接件的数值模型,探究了焊钉的尺寸和与钢板的连接方式如何影响组合梁的承载力,建议焊钉的长径比应大于4以防止拉拔破坏。蒋丽忠[11]等和余志武[12]等分别对钢–混凝土组合梁在均布和集中荷载下的滑移和变形进行了理论计算,推导出了相关的计算公式。罗子文[13]等和薛伟辰[14]等通过静力推出试验研究了栓钉的各种参数如何影响其抗剪性能,结果表明栓钉直径和混凝土强度的增加可提高抗剪承载力。丁发兴[15]等通过试验和ABAQUS模拟研究了栓钉抗剪承载力,发现随直径和屈服强度的增大而提高。蔺钊飞[16]进行了大直径栓钉的推出试验,分析了破坏模式和抗剪刚度,并与规范计算结果对比,发现规范在大直径栓钉抗剪承载力计算上较为保守。
剪力连接件作为钢–混组合结构的关键部件,其主要功能是将钢结构与混凝土结构组合构成一个整体,它的性能优越与否直接影响钢–混凝土组合结构构件共同受力与协调变形性能。剪力连接件的抗剪性能一般通过推出试验得到,但费时费力;而通过数值分析方法可以模拟推出试验从而预测剪力件性能,再经过少量的试验验证即可,是一种值得推广的方法[17]。
在波纹钢腹板组合箱梁桥结构中,剪力连接件的选择和设计非常关键,它将钢材与混凝土两种材料组合并使其共同工作。其具有如下特点:受力性能良好,各方向受力均匀且抗剪刚度相同,设计布置时不需要考虑受力的方向性,栓钉便于就地取材,对焊接工艺要求质量较高。栓钉剪力连接件是通过在钢翼缘板上焊接栓钉并与波纹钢腹板结合在一起而形成栓钉剪力连接件,其结构形式如图1所示:
Figure1.Schematic diagram of the cutting force connector
图1.栓钉剪力连接件示意图
为研究栓钉连接件在实际工程中受剪性能的表现,以某地区在建波形钢腹板钢–混组合梁桥钢混结合部位为研究对象,采用有限元分析软件ANSYS APDL建立栓钉连接件推出试件的精确三维有限元局部模型,研究栓钉剪力连接件的抗剪滑移性能,研究不同混凝土强度、栓钉强度、对剪力连接件的抗剪承载力和滑移性能的影响。
2. 有限元分析
2.1. 有限元模型
利用有限元软件ANSYS APDL编程语言以某波形钢腹板钢箱组合梁桥为工程背景,局部模型及网格划分如图1所示,材料参数、结构尺寸以实际工程为准,如表1和表2所示。因剪力钉和钢板在受剪过程中,栓钉的材料性质起主导作用,所以在建模时,剪力钉和翼缘板统一采用栓钉的材料性质,取GBT10433-2002《电弧螺柱焊用圆柱头焊钉》中ML15AL的强度。混凝土取C60,钢筋取HPB300。整体采用3D实体实体单元建模,混凝土与钢筋选用SOLD65实体单元进行一体式建模,用实常数定义钢筋的位置以及配筋率,混凝土采用多线性随动强化模型(MKIN)和混凝土专用破坏准则(CONCR)。栓钉和钢板采用SOLID185实体单元进行建模,选用多线性等向硬化模型(MISO)。
Figure 2.Finite element model and boundary conditions
图2.有限元模型及边界条件
通过施加面约束的形式来给模型添加边界条件,对于混凝土,在与Y轴垂直的底面施加面约束限制其在Y方向的位移,在与X轴垂直的混凝土两侧面添加面约束限制混凝土X方向的位移,在垂直于Z轴的混凝土顶面施加面约束限制其在Z方向的位移。对于钢翼缘板,在垂直于X轴的两侧施加面约束限制其在X方向位移,在垂直于Z轴的钢翼缘板底面上施加限制Z方向的面约束。在施加荷载时,将荷载F均分两半,对称施加在垂直于Y轴的两端面上,顶端面受压,底端面受拉,形成合力F。定义单荷载步,多荷载子步逐步加载,直至结构材料屈服而破坏。网格划分后的有限元模型如图2所示。
Table 1.Finite element model material parameters
表1.有限元模型材料参数
材料 |
混凝土 |
钢翼缘板 |
栓钉 |
钢筋 |
型号 |
C60 |
Q390ME |
ML15AL |
HPB300 |
Table2.Finite element model size parameters
表2.有限元模型尺寸参数
尺寸项目 |
混凝土板厚 |
栓钉直径 |
栓钉长度 |
翼缘板厚度 |
参数值(mm) |
250 |
16 |
150 |
20 |
2.2. 本构关系
混凝土的本构模型选用Hongestad公式,如下式所示:
(1)
式中:fc为峰值应力(棱柱体抗压强度),ε0为达到峰值应力时所对应的应变,εu为极限压应变。
剪力钉的本构模型选用三折线模型,如下式所示:
(2)
式中:σi为钢材的等效应力,fy为钢材的屈服强度,fu为钢材的极限强度,栓钉一般屈服强度较高(300~600 MPa之间),弹塑性强化现象较明显,屈强比较高,根据本文和已有试验结果[15],fu/fy在1.2~1.3之间,取中间值1.25,fu= 1.25fy,ES为钢材弹性模量,取ES= 2.06 × 105MPa,εi为钢材的等效应变,εy为钢材屈服时的应变,εy为钢材达极限强度时的应变,εu= 21εy。
钢筋的材料本构模型选用理性弹塑性模型,如下式所示:
(3)
2.3. 接触界面模拟
因为涉及不同性质的材料间的相对滑移,所以需要在交界面上定义接触面,来模拟交界面上的相对滑移。对于接触界面的模拟,根据与Solid65与Solid185实体单元特性,选择与之相匹配的conta173和targe170来模拟。由ANSYS接触面定义的刚–柔原则,混凝土表面用Conta173模拟,剪力钉与钢翼缘板表面用Targe170模拟。为模拟受力作用时交界面的真实情况,接触界面的接触模式定义为标准接触,接触面法向传递压力,允许接触面法向受拉分离,切向允许交界面产生相对滑移,摩擦系数取为0.25,为防止在非线性计算时由于过大的侵入而在材料屈服前收敛,根据下覆面应力大小,在每一荷载子步自动调整接触面的刚度大小。
3. 剪力连接件受力机理及有限元应力与变形分析结果
国内外学者通过研究表明,栓钉剪力连接件底部承受的压力最大沿着高度方向逐渐减小,在其顶部出现拉应力,压应力沿栓钉高度的变化规律如图所示。在实际的组合梁结构中,栓钉承受双向压应力,受力形式如同一根弹性地基梁,在钢–混组合界面处传递剪力能够提高剪力连接件的抗剪强度,弹性地基梁模型如图3所示。
Figure3.Dipping diagram of the stress of the nail
图3.栓钉受力示意图
Figure4.The result of the stress cloud diagram when it is destroyed
图4.破坏时应力云图结果
如图4所示,由栓钉翼缘板应与翼缘板有限元分析力云图分析可得,栓钉与翼缘板上的最大应力产生在栓钉根部,最终因栓钉根部应力达到极限而破坏,符合栓钉破坏的实际情况,综上验证了通过上述方法进行有限元分析的合理性。
3.1. 抗剪承载力与理论值
《钢结构设计标准》中给出单个栓钉连接件抗剪承载力计算式为[18]:
(4)
式中:As为栓钉截面面积,fc为混凝土抗压强度,fu为栓钉抗拉强度。
《钢–混凝土组合桥梁设计规范》中规定单个栓钉连接件抗剪承载力应取如下两式中的较小值[19]:
(5)
式中:η为群钉折减系数,取1.0。
规范Eurocode4规定栓钉连接件抗剪承载力计算式为
(6)
式中:α为栓钉高度影响系数,α= 0.2 (hs/ds+ 1) ≤ 1.0;
为混凝土圆柱体抗压强度;γy为栓钉抗力分项系数,取1.25。
规范AASHTO2007中给出栓钉连接件抗剪承载力计算式为
(7)
式中:φsc为抗力系数,取0.85;Qn为抗剪承载力名义值。
将抗剪承载力公式计算计算值与模拟所得结果进行比较,如表3所示。各国规范所建议值是加了安全系数之后的计算值,取值往往小于栓钉剪力连接件的实际值,由表3可得出,欧洲规范和《钢–混凝土组合桥梁设计规范》的计算值最为保守,设计值也更加偏于安全,而美国规范所建议公式的计算值则相比较更加激进,计算结果更加偏向于实际值,美国规范建议公式在去除抗力系数之后,承载力计算值为80.39 kN,与有限元计算结果很接近,误差为0.74%,验证了有限元模型计算结果的准确性。
Table3.Calculation results of load bearing capacity
表3.承载力计算结果
规范 |
抗剪承载力(kN) |
与本模型结果比值 |
《钢结构设计标准》 |
56.27 |
0.71 |
《钢–混凝土组合桥梁设计规范》 |
51.20 |
0.64 |
规范Eurocode4 |
51.45 |
0.64 |
规范AASHTO2007 |
68.33 |
0.86 |
本模型计算结果 |
79.80 |
1 |
3.2. 荷载位移曲线分析
钢–混凝土组合梁交界面的相对滑移效应是不能忽视的,并有学者以滑移量限制作为栓钉连接件断裂的破坏标准。因而,在研究抗剪连接件承载能力的同时,通常伴随着组合梁交界面相对滑移的研究。
(1) Buttry提出了分式形式的荷载–位移曲线[20]
(8)
式中:Vmax为所能承载最大剪切力,s为位移值
(2) Fisher等提出了指数形式的荷载–位移曲线
(9)
(3) N.Gattesco等在指数的基础上加了修正,提出了修正的指数形式[21]
(10)
Figure5.Comparison of finite element calculated values, equation calculated values, and measured values for load-displacement curves
图5.荷载–位移曲线的有限元计算值、公式计算值以及实测值比较
如图5所示,为验证有限元模型结果的正确性,从ANSYS APDL当中提取模型的荷载–位移曲线图,通过与上述几位学者的荷载–位移理论曲线研究结果以及文献中实测数据[22]进行比对,以横轴为模型栓钉与混凝土的相对滑移值,纵轴为抗剪承载力与最大抗剪承载力的比值。通过分析得出有限元模型的结果曲线整体走势与三个理论曲线吻合程度良好,且与Fisher的公式曲线的吻合程度最好,验证了有限元模型的正确性。
4. 参数化分析
4.1. 不同钢材强度
以混凝土立方体抗压强度为60 MPa、直径为16 mm的剪力连接件模型为例,现将栓钉的极限强度分别调整为360 MPa、380 MPa、420 MPa和440 MPa,计算并绘制不同栓钉屈服强度的荷载相对滑移曲线如图6所示。随着栓钉强度的增大,其极限承载能力也呈现明显的逐级递增趋势。栓钉极限强度360 MPa极限承载能力72.2 kN,栓钉极限强度380 MPa为78.3 kN,栓钉极限强度400 MPa为79.8 kN,栓钉极限强度420 MPa为84.05 kN,栓钉极限强度440 MPa为87.65 kN,各级栓钉强度之间承载能力增加分别为5.7%、4.6%、5.3%、4.3%。
Figure6.Different nail strength loads-bit displacement curve summary
图6.不同栓钉强度荷载–位移曲线汇总
4.2. 混凝土强度
以栓钉屈服强度为360 MPa,直径为16 mm的剪力连接件模型为例,现将混凝土的强度等级分别调整为C40、C60、C70和C80,计算并绘制不同混凝土强度的荷载–相对滑移曲线如图7所示。随着混凝土强度的提升,栓钉剪力连接件的承载能力有较小幅度的提升,从C40到C50,承载力从72.96 kN到78.05 kN,提升了6.98%;从C50到C60,承载力从78.05 kN到79.85 kN,提升了2.4%;从C60到C70,承载力从79.85kN到80.50kN,提升了0.8%;从C70到C80,极限承载力从80.50 kN到81.64 kN,提升了1.4%。从C40提升到C50时,承载力的提升是最大的,之后随着混凝土强的提升承载力并未提升太多,说明此时混凝土强度的不再是影响承载力的主要因素。
Figure7.Different concrete strength lotus load-bit displacement curve summary
图7.不同混凝土强度荷载–位移曲线汇总
4.3. 栓钉直径
Figure8.Different nail diameter load-displacement curve
图8.不同栓钉直径荷载–位移曲线
以栓钉屈服强度为360 MPa,长度为150 mm的剪力连接件模型为例,现将栓钉直径分别调整为11 mm、13 mm、19 mm和22 mm,计算并绘制不同混凝土强度的荷载–相对滑移曲线如图8所示。随着栓钉直径的增大,其极限承载能力也呈现明显的逐级递增趋势,而且增加量也较大。栓钉直径11 mm时,有限元模型极限承载能力51.6 kN;栓钉直径13 mm时型试件为59.8 kN,较11 mm直径栓钉上升了15.89%;栓钉直径16 mm时,有限元模型极限承载力为79.8 kN,较13 mm直径栓钉提升了33.4%;栓钉直径为19 mm时,有限元模型极限承载力为104.05 kN,较16 mm直径栓钉提升了30.39%;栓钉直径为22 mm时,极限承载力为128.05 kN,较19 mm直径栓钉上升了23%。而对于极限位移值,其随着栓钉直径的增加,极限位移值基本维持在同一水平,增加栓钉直径发现极限位移值差别较小,说明栓钉直径的增加对于极限位移值的影响较小。
4.4. 栓钉长度
以栓钉屈服强度为360 MPa,直径为16 mm的剪力连接件模型为例,现将栓钉长度分别调整为80 mm、120 mm和220 mm,计算并绘制不同混凝土强度的荷载–相对滑移曲线如图9所示。当栓钉长度为80 mm时,有限元模型的极限承载力为72.04 kN;当栓钉长度为120 mm时,有限元模型的极限承载力为77.6 kN;栓钉长度为150 mm时,为79.8 kN;栓钉长度为180 mm时,为81.76 kN;栓钉长度为220 mm时,承载力为81.76 kN。栓钉长度从80 mm增至180 mm,每增加一级其承载力增加了7.7%,2.8%,2.5%,而栓钉长度从180 mm到220 mm,极限承载力没有增加。随着栓钉长度的增加,极限承载力的增加在逐渐减少,长度增加至一定程度后对承载力大小无影响。
Figure9.Different armor length load-bit displacement curve summary
图9.不同栓长度荷载–位移曲线汇总
5. 结论
本文主要通过有限元分析软件ANSYS APDL建立栓钉连接件推出试件的三维有限元模型,研究构造参数对栓钉剪力连接件的对抗剪承载力影响,得出以下结论:
(1) 通过有限元分析计算,将有限元模型计算所得承载力值与不同规范所建议的公式计算值比较发现规范计算结果偏保守,相对而言美国规范计算值与有限元计算结果最接近。
(2) 有限元计算的荷载–相对滑移曲线和国外经验公式预测值对比,二者吻合较好,因此,采用非线性有限元法可以较好地模拟栓钉推出试验。
(3) 参数化分析表明:栓钉剪力连接件受剪承载力随混凝土强度提高而提高,但增速明显变慢;受剪承载力随栓钉屈服强度增加而增加且几乎呈线性增加;栓钉直径对抗剪承载力影响很小,长度达到一定程度承载力不再增加;栓钉直径对于抗剪承载力的影响最大,增加栓钉直径明显提高栓钉连接件抗剪承载力。