波纹钢板纵向接缝高强度螺栓连接承载力数值分析

苏明周，赵 凯

(1.西安建筑科技大学 土木工程学院,西安 710055;2.结构工程与抗震教育部重点实验室(西安建筑科技大学),西安 710055)

1787年波纹钢板诞生于英国[1]，随着科学技术的发展，波纹钢板作为一种特殊结构形式被主要应用于矿山、冶金、化工、污染等领域的涵洞、地下管道、储存工程[2].波纹钢板结构具有良好的受力性能和位移补偿特征，单位面积钢波纹板的惯性矩、截面模量、回转半径远大于平板，具有较强的抗弯能力和抗压能力[3].波纹的存在使局部力的传递方向及变形趋势发生改变，减小局部应力集中现象，起到分散作用[4].

目前，国内外学者对波纹钢板结构研究大多针对整体受力性能，对波纹钢板高强度螺栓连接承载力研究较少，缺少连接承载力的计算方法.现行规范中，美国ASTM A796/A796M—2010[5]设计标准以列表形式给出现有拼装波纹板连接件在具体构造下的承载力，不便于工程设计采用；加拿大《Canadian highway bridge design code》[6]规范中提到埋置式波纹管接缝强度验算按加拿大《Steel structures for buildings》[7]中平板连接计算；中国GB/T 34567—2017《冷弯波纹钢管》[8]中对冷弯波纹钢管进行了分类并在设计和连接上提出相应的规定，但并未给出连接承载力的计算方法.中国大多数学者在验算波纹钢管螺栓连接接缝强度时，一直采用平板连接计算公式，未考虑波纹板与平板受力性能的差异[9-12].

本文基于本课题组前期试验研究成果[13]，采用ABAQUS软件建立有限元模型，并以端距、螺栓预拉力、波形为参数，探讨影响波纹板连接件承载力的因素，深入了解波纹板连接件的受力机理，提出波纹板连接承载力计算方法，为工程应用提供设计建议.

1 有限元建模

有限元模型采用文献[13]中平板试件和波纹板试件(见图1、2)，具体试件参数见文献[13].有无螺纹对高强度螺栓的抗剪承载力基本无影响[14]，故高强度螺栓模型忽略螺纹的存在.

1.1 材料本构关系

有限元模型中钢材应力-应变关系采用“双折线”模型(见图3(a))，对应的参数设置依据板材的材性试验结果，有限元分析中将板材极限应变与抗拉强度试验名义值均转换为真实值(见表1).表中CP代表波纹板试件，PP代表平板试件，平板试件与波纹板试件材性试样相同.由于螺栓一般由中碳钢或低合金钢经热处理后制成，强度较高，其应力-应变曲线没有较大的流幅，故螺栓应力-应变关系采用“理想弹塑性”模型(见图3(b))，未考虑其强化能力[2]，fy=664 MPa、E=2.06×105.

图1 平板连接试件简图(mm)

图2 波形板连接试件简图(mm)

图3 材料应力-应变关系曲线

表1 有限元材性取值

1.2 单元类型及网格划分

连接件模型选用八节点减缩积分的一阶三维实体单元C3D8R进行模拟，该单元适用于接触分析，可以考虑材料非线性，计算结果较精确.连接件网格划分采用扫略式方法，此方法对于不规则模型易得到形状规则的单元.由于高强度螺栓及栓孔处易发生应力集中，故连接件整体种子尺寸取4 mm，螺栓及栓孔局部种子尺寸取3 mm.网格划分情况见图4.

1.3 相互作用

模型中的接触均采用面面接触，法向接触定义为硬接触，切向接触定义为库伦摩擦接触[15].抗滑移系数根据试验结果取值，平板试件、波纹板试件的抗滑移系数试验值分别为0.34和0.36.

1.4 边界条件及加载方式

连接板模型左侧端面为固定约束，右侧端面Z向施加位移进行加载，限制其他5个方向自由度.高强度螺栓预拉力通过施加紧固力实现，即在螺栓横截面上施加荷载控制预拉力大小，螺栓长度不发生改变，对于8.8级M20高强度螺栓，模型验证时取预拉力P=125 kN，对应试验试件为PPX-3/CPX-3.整体模型约束施加情况见图5.

图4 模型网格划分

图5 模型边界条件

2 有限元模型验证

2.1 破坏模式对比

因板厚4 mm与6 mm连接件、8 mm与10 mm连接件具有相似的变形特征和破坏形态，此处仅以板厚4、8 mm连接件为例进行有限元与试验结果对比.图6为板厚4 mm试件破坏模式对比图，PP1-3模型孔壁周边发生较大程度的鼓曲并有大范围环状屈服区域，栓孔被拉长，与试验结果相吻合.CP1-3试件模型孔壁四周应力沿连接件横向呈带状分布，屈服区域延伸至板件边缘.提取CP1-3模型端部单元应变值，已达到极限应变，可判断板件端部发生剪断，与试验破坏模式相同.

图6 板厚4 mm试件对比

图7为板厚8 mm试件破坏模式对比图，PP3-3和CP3-3试件破坏模式均为螺杆剪切破坏，栓孔变形较小，螺杆发生较大剪切变形.PP3-3试件模型破坏时板件距端部约55 mm范围内整体翘曲，未发生孔壁周边小范围凸起现象，与试验现象相同.CP3-3试件模型孔壁处仅发生小面积屈服，螺杆因产生较大区域剪切屈服而破坏，与试验现象吻合良好.

图7 板厚8 mm试件对比

2.2 荷载-位移曲线对比

CP1-3～CP4-3、PP1-3～PP4-3试件荷载-位移曲线的有限元与试验结果对比见图8～11.

图8 板厚4 mm试件荷载-位移曲线对比

由图8(a)～11(a)可看出有限元模型屈服后的刚度大于试验值，图8(b)～11(b)可看出波纹板连接件有限元分析荷载-位移曲线与试验结果吻合较好，原因在于波纹板栓孔主要沿横向变形，纵向波纹对栓孔变形存在约束，使其应力分布均匀，而平板试件栓孔向四周均发生变形，塑性变形范围较大.模型的破坏准则为板件的最大应变达到材料的极限应变或螺栓达到全截面屈服.如果荷载-位移曲线有下降段，抗剪承载力取曲线峰值；若无下降段，则偏于安全取模型破坏时的最大值.

图9 板厚6 mm试件荷载-位移曲线对比

图10 板厚8 mm试件荷载-位移曲线对比

图11 板厚10 mm试件荷载-位移曲线对比

2.3 受剪承载力对比

试验与有限元分析的受剪承载力对比见表2，表中有限元计算的承载力取值为：当荷载-位移曲线有下降段时，连接件承载力取曲线的峰值；当曲线无下降段时，承载力偏于安全的取连接件破坏时的最大值.从表2可得出有限元分析与试验结果相对误差在5%以内，表明有限元分析与试验结果吻合良好，存在微小误差的主要原因为：①模型为理想轴心受力状态，试验中存在一定的偏心加载；②建模时未考虑材料本身缺陷.

表2 试验与模拟分析受剪承载力

3 参数分析

根据试验与有限元验证可知板件端距和预拉力是影响波纹板连接件承载力的主要因素，故参数分析时考虑端距、预拉力和波形.同时，对中国GB/T 34567—2017《冷弯波纹钢管》[8]中多种波形尺寸波纹板连接件进行建模分析，对比平板连接件承载力，给出波纹板连接承载力计算方法的建议.

3.1 端距对承载力的影响

通过2.1节有限元分析可知平板连接件与波纹板连接件在受拉时端距方向应力分布存在差异，两种连接件端距限值可能不同，故分析平板连接件与波纹板连接件不同端距时极限承载力与极限位移的变化，端距取值为45、50、55、60、65、70 mm(2d0～3d0).板厚为4、 6、 8、 10 mm连接件不同端距的荷载-位移曲线见图12～15.

由图12(a)～15(a)看出随着端距的增加，平板连接件的荷载-位移曲线变化较小，即连接件极限位移与承载力在板件端距在不小于45 mm时不发生明显变化，符合GB 50017—2017《钢结构设计标准》[16]中端距限值要求(2d0).

由图12(b)～14(b)可看出波纹板连接件在板厚为4、6、8 mm时荷载-位移曲线随端距的增加变化较大，板厚为4、6 mm的连接件破坏模式随端距增加由端部剪断过渡为孔壁承压破坏(见图16).当破坏模式为孔壁承压时，端距为60 mm(2.7d0)，承载力不再增加.

图12 板厚4 mm连接件不同端距下荷载-位移曲线

图13 板厚6 mm连接件不同端距下荷载-位移曲线

图14 板厚8 mm连接件不同端距下荷载-位移曲线

图15 板厚10 mm连接件不同端距下荷载-位移曲线

图16 6 mm波纹板连接件应力云图随端距变化

为使分析结果具有一般性，对中波、大波、深波波形连接件进行不同端距下极限承载力有限元分析.中波尺寸为150 mm×50 mm、200 mm×55 mm、230 mm×64 mm，深波尺寸为300 mm×110 mm、380 mm×140 mm，大波尺寸为400 mm×150 mm，中波波形连接件采用8.8级M20高强度螺栓连接(d0=22 mm)，深波与大波波形连接件均采用8.8级M24高强度螺栓连接(d0=26 mm).计算结果见表3、4，表中F45表示端距取45 mm时连接件承载力，“―”表示无此项.

表3 中波波形连接件不同端距下的承载力

表4 大波、深波波形连接件不同端距下的承载力

由表3、4可知中波连接件在端距取值不小于60 mm(2.7d0)时，承载力趋于稳定；大波、深波连接件端距取值不小于70 mm(2.7d0)时，承载力基本无明显差异.因此，波纹板连接件不发生端距端部剪断的端距可取3d0.

3.2 预紧力对承载力的影响

根据文献[13]可知波纹板连接件孔壁处为三向受力状态，螺栓预拉力影响孔壁承压强度.为深入探讨预拉力对连接件极限承载力的影响，在其他设计参数不变的情况下，板厚为4、6、8、10 mm的波纹板连接件端距取60 mm(2.7d0)时的不同预拉力下荷载-位移曲线见图17.

由图17(a)可看出板厚为4 mm波纹板连接件承载力随预拉力的增加而增加，主要原因为波纹板连接件受拉时孔壁处为三向应力状态，预拉力的增大导致板件承压能力增强.图17(b)～(d)中预拉力对连接件的荷载-位移曲线基本无影响，可知发生螺杆剪切破坏的连接件承载力基本不受预拉力影响.由图17可看出无论试件发生何种破坏，试件的滑移荷载随着预拉力的增加而增加.

图17 不同预拉力时的荷载-位移曲线

3.3 波形对承载力的影响

不同波形波纹板连接件、平板连接件荷载-位移曲线见图18，中波波纹板连接件的端距取值为60 mm，大波、深波波纹板连接件的端距取值为70 mm，平板连接件的端距取值为45 mm；M20高强度螺栓预拉力为125 kN，M24高强度螺栓预拉力为175 kN.图中M20-PB表示采用M20高强度螺栓连接的平板连接件，M24-PB表示采用M24高强度螺栓连接的平板连接件，图中波纹板连接件的波形尺寸单位为mm×mm.

图18 不同波形连接件荷载-位移曲线

由图18可看出3种中波波形连接件的荷载-位移曲线变化趋势较一致，大波与深波波形连接件荷载-位移曲线趋势基本相同.当板厚不小于6 mm时，3种中波波形连接件均发生螺杆剪切破坏，承载力波动较小.大波与深波波形连接件在板厚不小于7 mm时发生螺杆剪切破坏，承载力无明显变化.

由图18可知采用M20和M24高强度螺栓连接的平板连接件分别在板厚不小于8 mm和9 mm时，连接件发生螺杆剪切破坏，承载力变化较小.因波纹板连接件与平板连接件承载力比值应在同等破坏模式下进行，故当波纹板连接件发生螺杆剪切破坏时，平板连接件对应承载力应为同等破坏模式的承载力，即M20-PB试件在板厚不小于6 mm时承载力均取发生螺杆剪切破坏时承载力157.5 kN，M24-PB试件在板厚不小于7 mm时承载力均取发生螺杆剪切破坏时承载力236.5 kN.板厚为5、6、7、8、9和10 mm的中波、大波和深波波形连接件承载力与同等条件下的平板连接件承载力见表5.

由表5、6可看出当波纹板连接件与同等条件平板连接件均发生孔壁承压破坏时，承载力比值范围为1.19～1.41.当两种类型连接件均发生螺杆剪切破坏时，两者承载力基本相同.因此，当波纹板连接件发生孔壁承压破坏时，建议在中国平板连接承压计算公式乘以1.1的增大系数做为波纹板连接接缝承载力计算公式；当波纹板连接件发生螺杆剪切破坏时，抗剪连接承载力按平板连接计算公式计算.

表5 M20高强度螺栓连接件承载力

表6 M24高强度螺栓连接件承载力

4 结 论

1)利用ABAQUS软件建立了波纹板连接件及平板连接件的有限元模型，有限元分析结果与试验结果吻合较好，验证了模型的有效性.

2)发生孔壁承压破坏的波纹板连接件，承载力大于同等条件的平板连接件，并且承载力随螺栓预拉力的增加而增加.发生螺杆剪切破坏的波纹板连接件，承载力与同等条件的平板连接件承载力基本相同，预拉力对其承载力影响较小.

3)为保证波纹板高强度螺栓连接件不发生端部剪断，板件端距应不小于3d0.

4)当波纹板连接件按孔壁承压破坏设计时，建议承载力按中国平板连接承压计算公式乘以1.1的增大系数；当波纹板连接件按螺杆剪切破坏设计时，抗剪连接承载力按平板连接计算公式计算.