胡文旭，陳寶春,2，李聰,3，蘇家戰(zhàn),2，黃文金

(1.福州大學(xué)土木工程學(xué)院，福建 福州 350108；2.橋梁技術(shù)創(chuàng)新與風(fēng)險(xiǎn)防治國(guó)家級(jí)國(guó)際聯(lián)合研究中心，福建 福州 350108； 3.廣西大學(xué)土木建筑工程學(xué)院，廣西 南寧 530004；4.福建農(nóng)林大學(xué)交通與土木工程學(xué)院，福建 福州 350108)

0 引言

超高性能混凝土(ultra-high performance concrete，UHPC)是一種超高強(qiáng)、高韌性、高彈性模量的新型水泥基復(fù)合材料[1].利用UHPC的上述優(yōu)良特性，將其作為應(yīng)力過(guò)渡層與鋼橋面結(jié)合形成的鋼-UHPC組合橋面板，是解決瀝青鋪裝層損壞和鋼橋面板疲勞開裂的有效途徑之一[2].

栓釘為鋼-UHPC組合板常用的剪力連接件，能承受和傳遞UHPC與鋼板間的水平剪力(主要)和豎向掀起力，是結(jié)構(gòu)協(xié)同受力的關(guān)鍵部件.為此，已開展許多試驗(yàn)研究[3-6]，結(jié)果表明，UHPC中的栓釘均發(fā)生釘桿剪斷破壞，而UHPC僅在栓釘根部位置出現(xiàn)局部壓潰現(xiàn)象，抗剪承載力主要受栓釘直徑、強(qiáng)度及栓釘焊縫的影響.此外，為了防止橋面板邊緣發(fā)生豎向掀起，需對(duì)板外周一圈的栓釘進(jìn)行加密，因此該處可能存在群釘效應(yīng)，然而，目前尚未見鋼-UHPC組合板有關(guān)群釘效應(yīng)的試驗(yàn).

由于鋼-UHPC組合橋面板栓釘連接件推出試驗(yàn)數(shù)量有限，且試驗(yàn)結(jié)果受試驗(yàn)量測(cè)方法的限制，常通過(guò)開展有限元分析來(lái)深入了解栓釘?shù)目辜魴C(jī)理，探究影響栓釘抗剪性能的因素.但在現(xiàn)有有限元分析中，參數(shù)分析多集中于栓釘?shù)闹睆?、長(zhǎng)度、強(qiáng)度和UHPC強(qiáng)度[4,7-9]，未見對(duì)鋼-UHPC組合板群釘效應(yīng)的分析.

為此，針對(duì)鋼-UHPC組合板中栓釘?shù)娜横斝?yīng)，建立考慮栓釘焊縫影響的抗剪有限元模型，以文獻(xiàn)[5]試驗(yàn)為基礎(chǔ)數(shù)據(jù)進(jìn)行模型驗(yàn)證.應(yīng)用經(jīng)驗(yàn)證的模型，通過(guò)加密栓釘間距，分析群釘效應(yīng)對(duì)鋼-UHPC組合板中栓釘連接抗剪性能的影響，并提出群釘效應(yīng)折減系數(shù).

1 分析對(duì)象

本研究的有限元分析以文獻(xiàn)[5]的試件為原型.文獻(xiàn)[5]中的推出試驗(yàn)試件，分SI(單排雙列)和SII(雙排雙列)組，其中UHPC板厚50 mm，配有Φ10 mm@50 mm鋼筋網(wǎng)，鋼構(gòu)件的翼緣板尺寸為360 mm×550 mm×14 mm，腹板尺寸為272 mm×550 mm×14 mm，并焊接ML15短栓釘，直徑為13 mm，高度為45 mm(焊后高度為40 mm).UHPC的抗壓強(qiáng)度為175.2 MPa，抗拉強(qiáng)度為10.2 MPa，彈性模量為46.2 GPa；栓釘?shù)那?qiáng)度為350.0 MPa，極限強(qiáng)度為440.0 MPa，彈性模量為207.0 GPa；鋼構(gòu)件和鋼筋的屈服強(qiáng)度分別為345.5和500.4 MPa.

2 有限元模型

2.1 單元選取和接觸關(guān)系

模型中的鋼構(gòu)件、UHPC 板和栓釘(包括焊縫)采用三維實(shí)體單元(C3D8R)進(jìn)行模擬；鋼筋網(wǎng)采用兩節(jié)點(diǎn)三維桁架單元(T3D2)進(jìn)行模擬.UHPC板與鋼構(gòu)件之間采用切向無(wú)摩擦、法向硬接觸的面-面接觸方式.栓釘頭部的上、下表面與UHPC的接觸同樣采用切向無(wú)摩擦和法方硬接觸[10].栓釘側(cè)面和焊縫表面與UHPC的接觸均采用法向硬接觸、切向采用罰函數(shù)摩擦，方向摩擦系數(shù)取0.4[11].鋼筋網(wǎng)與UHPC的接觸采用內(nèi)置形式(Embed)進(jìn)行模擬，而栓釘與鋼構(gòu)件間采用綁定的接觸方式.

2.2 邊界條件

考慮到試件具有雙軸對(duì)稱形式，僅建立1/4模型進(jìn)行數(shù)值模擬.在邊界條件上，對(duì)兩個(gè)對(duì)稱面施加對(duì)稱邊界.對(duì)稱面1需對(duì)x方向上平動(dòng)進(jìn)行約束，y軸、z軸限制轉(zhuǎn)動(dòng)自由度；對(duì)稱面2需對(duì)z方向上平動(dòng)進(jìn)行約束，將x軸、y軸限制轉(zhuǎn)動(dòng)自由度.約束鋼構(gòu)件頂面的所有方向上的的位移和轉(zhuǎn)動(dòng).

2.3 材料本構(gòu)2.3.1 鋼材

栓釘(包括焊縫)采用三折線模型[12]，在ABAQUS程序中延性金屬損傷包含延性損傷和剪切損傷兩類，有限元建模中同時(shí)為栓釘定義這兩種損傷[8].

據(jù)文獻(xiàn)[5]試驗(yàn)結(jié)果知，鋼構(gòu)件和鋼筋為碳素鋼材料，且應(yīng)變不大，故本構(gòu)關(guān)系采用理想彈塑性模型.

2.3.2 UHPC材料本構(gòu)關(guān)系

UHPC材料本構(gòu)選用ABAQUS中的混凝土損傷塑性模型.UHPC受壓本構(gòu)關(guān)系參考文獻(xiàn)[13]，如圖1(a)所示.受拉本構(gòu)關(guān)系參考文獻(xiàn)[14]，該本構(gòu)模型包含UHPC的初裂強(qiáng)度且充分體現(xiàn)鋼纖維對(duì)抗拉強(qiáng)度的貢獻(xiàn)，如圖1(b)所示.

圖1中，fc為UHPC抗壓強(qiáng)度；e為自然常數(shù)；ξ為應(yīng)變比，ξ=ε/ε0，ε為UHPC應(yīng)變，ε0為峰值應(yīng)變；n=E/E0，E為UHPC彈性模量，E0為峰值點(diǎn)的割線模量；fct為UHPC抗拉強(qiáng)度；K為鋼纖維增強(qiáng)因子，K=(lf/df)Vf，lf為鋼纖維長(zhǎng)度，df為鋼纖維直徑，Vf為鋼纖維體積率；m=0.85-0.47K+0.12K2,0

圖1 UHPC材料本構(gòu)曲線Fig.1 Material constitutive curves of UHPC

2.4 有限元模型驗(yàn)證2.4.1 試件的有限元模型

栓釘焊縫對(duì)抗剪承載力具有顯著影響，有的有限元模型中給以考慮[4,7 ,11]，有的卻沒(méi)有考慮[6,8-9]，尚未形成共識(shí).因此，分別建立考慮焊縫和不考焊縫的模型進(jìn)行對(duì)比，焊縫尺寸參考《電弧螺柱焊用圓柱頭栓釘(GB/T 10433—2002)》[15]的規(guī)定.文獻(xiàn)[7]對(duì)比了帶倒角和直角邊的焊縫有限元模型，結(jié)果表明采用直角邊的焊縫模擬焊縫效果能滿足要求，且便于建模，因此本研究在焊縫有限元模擬時(shí)也選用直角邊模型.

模型整體網(wǎng)格尺寸為8 mm，在栓釘(焊縫)、與栓釘相接觸的UHPC及鋼構(gòu)件局部區(qū)域內(nèi)采用更精細(xì)化的網(wǎng)格，網(wǎng)格尺寸不超過(guò)為4 mm.栓釘有限模型及網(wǎng)格劃分如圖2所示，邊界條件如圖3所示.整個(gè)模型共計(jì)16 370個(gè)節(jié)點(diǎn)和11 777個(gè)單元.

圖2 有限元模型圖Fig.2 Finite element model

圖3 邊界條件Fig.3 Boundary conditions

2.4.2 有限元分析結(jié)果與試驗(yàn)結(jié)果對(duì)比分析

圖4為有限元分析獲取的荷載-位移曲線與文獻(xiàn)[5]試驗(yàn)所得實(shí)測(cè)曲線對(duì)比結(jié)果.由圖可知，荷載-滑移曲線由彈性段OA、彈塑性段AB和下降段BC等3個(gè)階段組成，計(jì)算結(jié)果與試驗(yàn)實(shí)測(cè)曲線吻合良好，驗(yàn)證了該數(shù)值模型的正確性.

表1為有限元分析提取的極限荷載和位移與文獻(xiàn)[5]試驗(yàn)所測(cè)的對(duì)比.由表1可知，有限元模型中考慮焊縫影響時(shí)，無(wú)論是極限荷載還是極限滑移值均與試驗(yàn)值較為吻合，計(jì)算值與實(shí)測(cè)值的比值均值分別為1.04和0.94；而不考慮焊縫影響時(shí)，計(jì)算值與實(shí)測(cè)值的比值則相差較大，分別為1.37和1.31.這也表明，焊縫對(duì)栓釘?shù)目辜舫休d力有顯著提高，這與文獻(xiàn)[4]和文獻(xiàn)[11]的結(jié)論相符.所以，有限元模型中應(yīng)考慮焊縫的作用，否則計(jì)算結(jié)果過(guò)于保守.

圖4 有限元結(jié)果與試驗(yàn)對(duì)比Fig.4 Comparison between FEA and test results

表1 試驗(yàn)值與有限元結(jié)果對(duì)比Tab.1 Comparison between FEA and test results

3 推出試件受力過(guò)程有限元分析

3.1 栓釘受力過(guò)程分析

通過(guò)ABAQUS中輸出等效塑性應(yīng)變?cè)茍D(PEEQ)來(lái)判斷栓釘?shù)钠茐男螒B(tài).PEEQ>0表示單元開始屈服；PEEQ>0.2表示栓釘單元趨近破壞[11].以SII組試件為例進(jìn)行說(shuō)明.圖5表示栓釘從屈服到峰值再到破壞的過(guò)程(白色區(qū)域表示該區(qū)域趨近破壞，箭頭指示受力方向).由圖5可知，當(dāng)達(dá)到A點(diǎn)的荷載時(shí)，栓釘根部的焊趾區(qū)域最先進(jìn)入屈服狀態(tài)，表征模型將進(jìn)入彈塑性階段；隨著荷載增加，塑性區(qū)域擴(kuò)大，當(dāng)達(dá)到B點(diǎn)的荷載時(shí)，栓釘根部出現(xiàn)明顯變形，該區(qū)域已處于高塑性應(yīng)變區(qū)；當(dāng)達(dá)到C點(diǎn)的荷載時(shí)，白色區(qū)域貫穿栓釘根部，表征栓釘根部被剪斷，右側(cè)栓釘首先被剪斷，左側(cè)栓釘隨之被剪斷.

圖6為栓釘上6個(gè)測(cè)點(diǎn)的荷載-應(yīng)變曲線(紅色箭頭指示栓釘受剪方向)，其應(yīng)變測(cè)點(diǎn)沿栓釘高度方向的布置.

圖5 栓釘?shù)刃苄詰?yīng)變分布云圖Fig.5 Equivalent plastic strain of head studs

圖6 荷載-應(yīng)變曲線Fig.6 Load-strain curves

由圖5、6可知，栓釘在測(cè)點(diǎn)1處由于受到焊縫的約束，在加載全過(guò)程中基本為負(fù)值；加載初期，應(yīng)變?cè)鲩L(zhǎng)緩慢，接近B點(diǎn)后，增速加快.在測(cè)點(diǎn)2和4處，應(yīng)變?cè)诩虞d全過(guò)程中均為正值，處于受拉狀態(tài)；在B點(diǎn)之前，二點(diǎn)處的應(yīng)變變化趨勢(shì)與應(yīng)變值均接近，接近B點(diǎn)時(shí)，測(cè)點(diǎn)2的應(yīng)變繼續(xù)增大，而測(cè)點(diǎn)4的應(yīng)變逐漸變?。粶y(cè)點(diǎn)3和6點(diǎn)處的荷載-應(yīng)變曲線大致關(guān)于y軸對(duì)稱，表明栓釘中部的上側(cè)全程受拉，下側(cè)全程受壓.

3.2 UHPC板受力過(guò)程分析

通過(guò)在ABAQUS中輸出混凝土受壓、受拉損傷云圖(DAMAGEC和DAMAGET)來(lái)判斷UHPC的破壞形態(tài).當(dāng)應(yīng)力下降到峰值應(yīng)力的85%以下時(shí)，可認(rèn)為混凝土材料處于破壞狀態(tài)[16]，將其對(duì)應(yīng)的應(yīng)變定義為UHPC的極限壓應(yīng)變.由計(jì)算可知，DAMAGEC>0.24時(shí)UHPC單元受壓破壞，DAMAGET>0.57時(shí)UHPC單元受拉開裂.

圖7表示SII組試件UHPC板的破壞過(guò)程.由圖7(a)知，當(dāng)達(dá)到A點(diǎn)的荷載時(shí)，UHPC板與焊縫前端相互擠壓，焊趾處首先出現(xiàn)壓潰；當(dāng)達(dá)到B點(diǎn)的荷載時(shí)，損傷區(qū)域高度沿著焊縫前端向外擴(kuò)展；當(dāng)達(dá)到C點(diǎn)的荷載時(shí)，UHPC除局部壓潰外，其余部分基本無(wú)損傷.由圖7(b)知，在荷載作用下，栓釘后端與UHPC板脫空，因此開裂區(qū)域主要集中于栓釘孔后端區(qū)域.在試件受力的全過(guò)程中，UHPC板頂處均無(wú)開裂出現(xiàn).

圖7 UHPC損傷Fig.7 Damage of UHPC

3.3 栓釘與UHPC板相互作用

圖8 荷載-接觸應(yīng)力曲線Fig.8 Load-contact stress curves

圖8為荷載-栓釘與UHPC接觸應(yīng)力曲線.由圖可知，在荷載作用下，栓釘上部與UHPC分離，因此栓釘上部的點(diǎn)3、點(diǎn)6和點(diǎn)9與UHPC的接觸應(yīng)力始終為0.而栓釘下部與UHPC擠壓，在A點(diǎn)之前，栓釘根部變形不大，點(diǎn)1、點(diǎn)4和點(diǎn)7的接觸應(yīng)力隨著荷載增大呈線性增長(zhǎng).當(dāng)接近B點(diǎn)時(shí)，由于栓釘根部變形增大，與其接觸的UHPC區(qū)域進(jìn)入塑性區(qū)域，點(diǎn)1和點(diǎn)4的接觸應(yīng)力逐漸減小，而釘桿中部的點(diǎn)7處的UHPC未進(jìn)入塑性階段，因此接觸應(yīng)力仍逐漸增大.

4 群釘效應(yīng)有限元分析與設(shè)計(jì)建議

4.1 群釘效應(yīng)有限元計(jì)算結(jié)果

圖9為栓釘間距為15.4d(d為栓釘直徑，表示間距為15.4倍的栓釘直徑)時(shí)，直徑分別為13、16 和19 mm的荷載-滑移曲線.直徑16和19 mm的栓釘面積比直徑13 mm的栓釘面積，增大了51.5%和113.6%，承載力也得到了相似比例的提高(58.0%和100.7%).圖10為不同間距下直徑13 mm栓釘?shù)暮奢d-滑移曲線，可知隨著栓釘間距減小，群釘效應(yīng)程度加劇，栓釘抗剪承載力隨著減小.

圖9 不同直徑栓釘荷載-滑移曲線Fig.9 Different diameter stud load-slip curves

圖10 不同間距下直徑13 mm栓釘荷載-滑移曲線Fig 10 Different spacing of 13 mm diameter stud load-slip curves

圖11～13分別對(duì)比在承載力最大時(shí)，不同間距下栓釘直徑為13 、16 和19 mm模型中的UHPC受壓、受拉損傷云圖.對(duì)于直徑13 mm的栓釘，在間距為5.0d時(shí)，受壓高損傷應(yīng)變區(qū)域未發(fā)生重疊，無(wú)明顯的群釘效應(yīng)；在間距為3.5d時(shí)，出現(xiàn)重疊現(xiàn)象，導(dǎo)致混凝土對(duì)栓釘?shù)募s束作用顯著減弱，抗剪承載力下降幅度變大，如圖11所示.對(duì)于直徑16 和19 mm栓釘，在間距5.0d時(shí)就出現(xiàn)重疊，且在板頂出現(xiàn)開裂區(qū)域；當(dāng)間距為3.5d兩栓釘中間的UHPC已完全壓潰，同時(shí)板頂開裂區(qū)域擴(kuò)大，說(shuō)明在推出試驗(yàn)中UHPC可能發(fā)生開裂，如圖12～13所示.

圖11 直徑13 mm栓釘不同間距下UHPC損傷對(duì)比Fig.11 Comparison of damages of UHPC with different spacing of 13 mm diameter stud

圖12 直徑16 mm栓釘不同間距下UHPC損傷對(duì)比Fig.12 Comparison of damages of UHPC with different spacing of 16 mm diameter stud

圖13 直徑19 mm栓釘不同間距下UHPC損傷對(duì)比Fig.13 Comparison of damages of UHPC with different spacing of 19 mm diameter stud

進(jìn)一步，由圖14還可知，直徑19 mm栓釘?shù)乃ㄡ敹嗣鄙喜繀^(qū)域出現(xiàn)大面積開裂，可能發(fā)生UHPC板破壞.為此，同樣定義UHPC的極限拉應(yīng)變?yōu)?5%的峰值拉應(yīng)力對(duì)應(yīng)的應(yīng)變，計(jì)算得DAMAGEC>0.84表示UHPC單元受拉破壞.圖14對(duì)比了在承載力最大時(shí)，不同間距下直徑為19 mm栓釘模型的UHPC受拉破壞云圖.由圖可知，當(dāng)栓釘間距小于3.5d時(shí)，局部位于板頂?shù)腢HPC單元趨近破壞.因此，對(duì)于直徑為19 mm栓釘時(shí)，當(dāng)栓釘間距小于3.5d時(shí)，試件的破壞模式會(huì)轉(zhuǎn)變?yōu)閁HPC板破壞.

圖14 直徑19 mm栓釘不同間距下UHPC受拉破壞對(duì)比Fig.14 Comparison of tensile damage of UHPC with different spacing of 19 mm diameter stud

圖15為不同間距下直徑13、16和19 mm的栓釘抗剪承載力.為便于對(duì)比，以不同間距下栓釘抗剪承載力與15.4d間距下的比值定為群釘效應(yīng)折減系數(shù)，繪制了群釘折減系數(shù)隨著栓釘間距的變化圖，如圖16所示.

圖15 不同間距下栓釘抗剪承載力Fig.15 Stud shear capacity at different spacing

圖16 群釘效應(yīng)折減系數(shù)隨著栓釘間距的變化Fig.16 Reduction coefficient of group studs varies with the spacing

由圖16可知，對(duì)于直徑為13 mm栓釘，當(dāng)栓釘間距小于3.5d后，栓釘承載力折減系數(shù)迅速減小，表明群釘效應(yīng)程度明顯；當(dāng)栓釘間距為2.5d時(shí)，承載力的降幅為5.1%.然而，對(duì)于直徑16和19 mm的栓釘，當(dāng)栓釘間距小于5.0d后，栓釘承載力折減系數(shù)就開始迅速減小，說(shuō)明在相同栓釘間距下，隨著栓釘直徑的增大，群釘效應(yīng)程度也會(huì)加??；當(dāng)栓釘間距為8.0d時(shí)，承載力的降幅不到1.0%.

4.2 群釘效應(yīng)設(shè)計(jì)建議

我國(guó)的《公路鋼混組合橋梁設(shè)計(jì)與施工規(guī)范(JTG/T D64-01—2015)》[17]中規(guī)定，混凝土中的焊釘連接件剪力作用方向上的間距不宜小于焊釘直徑的5倍且不應(yīng)小于100 mm；而在我國(guó)《鋼-混凝土組合橋梁設(shè)計(jì)規(guī)范(GB 50917—2013)》[18]中規(guī)定，當(dāng)6≤β<13時(shí)，需考慮群釘效應(yīng)；當(dāng)β≥13時(shí)，不考慮群釘效應(yīng).這里，β=ld/d為栓釘間距與栓釘直徑的比值，β=ld/d，ld為栓釘間距.

然而，對(duì)于鋼-UHPC組合板，上述的規(guī)定并不一定能直接應(yīng)用.由節(jié)4.1分析結(jié)果可知，鋼-UHPC組合板水平剪切破壞時(shí)，是否出現(xiàn)UHPC破壞的栓釘最小間距，與栓釘直徑有關(guān).以不出現(xiàn)UHPC破壞為原則，根據(jù)節(jié)4.1的分析，建議直徑13 mm的栓釘，其間距不宜小于2.5d；直徑16和19 mm的栓釘，其間距不宜小于3.5d.

對(duì)于栓釘承載力考慮群釘效應(yīng)的折減系數(shù)ηG，其也與栓釘直徑有關(guān).由圖13可知，對(duì)于3種直徑的栓釘，當(dāng)β≥8.0時(shí)，均可不考慮群釘效應(yīng).對(duì)于直徑13 mm的栓釘，當(dāng)栓釘間距為2.5d時(shí)，該折減系數(shù)為0.949；對(duì)于直徑16和19 mm的栓釘，當(dāng)栓釘間距為3.5d時(shí)，該折減系數(shù)為0.958和0.964.考慮到安全性和方便，建議當(dāng)間距大于建議的最小值、小于8.0d時(shí)，ηG均取0.94.由于上述建議的規(guī)定是基于有限元模型分析結(jié)果，其合理性和具體值還有待今后進(jìn)一步的研究.

5 結(jié)論

1) 鋼-UHPC組合板的抗剪破壞以栓釘剪斷為主，栓釘與UHPC受剪過(guò)程的相互作用主要集中于栓釘根部焊縫處，栓釘焊縫能顯著提高栓釘?shù)目辜舫休d力.因此，有限元模型中應(yīng)考慮焊縫的作用，否則計(jì)算結(jié)果過(guò)于保守.

2) 栓釘抗剪承載力受直徑和間距的影響顯著.隨著栓釘直徑的增大和間距的減小，群釘效應(yīng)程度加劇，且UHPC板會(huì)出現(xiàn)開裂甚至可能破壞.現(xiàn)有《公路鋼混組合橋梁設(shè)計(jì)與施工規(guī)范(JTG/T D64-01—2015)》和《鋼-混凝土組合橋梁設(shè)計(jì)規(guī)范(GB 50917—2013)》有關(guān)栓釘最小間距和群釘效應(yīng)的規(guī)定，不能直接用于鋼-UHPC組合板.

3) 對(duì)于鋼-UHPC組合板，建議栓釘?shù)淖钚￠g距為13 mm栓釘不宜小于2.5d；16和19 mm栓釘不宜小于3.5d.對(duì)于栓釘承載力考慮群釘效應(yīng)的折減系數(shù)ηG，建議當(dāng)間距不小于8d時(shí)，可不考慮群釘效應(yīng)；當(dāng)間距大于建議的最小值且小于8.0d時(shí)，取0.94.上述基于有限元模型分析結(jié)果提出的建議，其合理性和具體值還有待今后的深入研究.