龐書孟，浣 石，陶為俊，譚湘倩

(廣州大學(xué)土木工程防護(hù)研究中心，廣東廣州 510006)

在混凝土中摻鋼纖維能增強(qiáng)其局部韌性，進(jìn)而可以提高混凝土的抗拉強(qiáng)度和主要由主拉應(yīng)力控制的抗彎、抗剪、抗扭強(qiáng)度［1－4］.試驗(yàn)?zāi)軌驕y(cè)量出動(dòng)態(tài)的應(yīng)力－應(yīng)變曲線，但對(duì)于其裂紋開展以及破壞過(guò)程難以觀測(cè)到.但運(yùn)用計(jì)算機(jī)通過(guò)細(xì)觀數(shù)值模擬不僅能夠觀測(cè)到微秒級(jí)的瞬態(tài)過(guò)程，還能夠觀測(cè)到混凝土的形貌演變過(guò)程以及目標(biāo)內(nèi)部應(yīng)力、應(yīng)變等全過(guò)程［5－8］.然而由于鋼纖維的摻入，使得混凝土內(nèi)部結(jié)構(gòu)更具復(fù)雜性，通過(guò)簡(jiǎn)單的數(shù)值模型難以準(zhǔn)確描述鋼纖維對(duì)混凝土局部的影響，因此，本文在試驗(yàn)的基礎(chǔ)上，通過(guò)軟件LS-DYNA建立含有鋼纖維和由骨料、砂漿和界面過(guò)渡區(qū)(ITZ)三相組成的細(xì)觀數(shù)值模型，對(duì)不同鋼纖維體積率下鋼纖維混凝土的沖擊壓縮過(guò)程進(jìn)行對(duì)比分析.

1 研究對(duì)象

本文細(xì)觀數(shù)值研究的對(duì)象:分離式霍普金森壓桿(SHPB)對(duì)鋼纖維混凝土進(jìn)行沖擊壓縮試驗(yàn)的過(guò)程.典型的SHPB試驗(yàn)裝置見圖1，包括發(fā)射裝置、撞擊桿、入射桿、透射桿和緩沖箱，以及位于透射桿和入射桿之間的試件，并在入射桿和透射桿分別貼上靈敏度較高的半導(dǎo)體應(yīng)變片用來(lái)測(cè)量?jī)筛鶋簵U中的應(yīng)力波信號(hào).在試驗(yàn)過(guò)程中撞擊桿、入射桿和透射桿均保持在線彈性變形狀態(tài)，且撞擊桿長(zhǎng)度為300 mm，直徑50 mm;入射桿長(zhǎng)度為2 700 mm，直徑74 mm;透射桿總長(zhǎng)為1 800 mm，直徑74 mm.SHPB試驗(yàn)的基本原理:發(fā)射裝置中的高壓氣體推動(dòng)撞擊桿，在炮管中加速使其獲得一定的速度，然后撞擊入射桿，同時(shí)在入射桿中產(chǎn)生一維壓縮應(yīng)力波，波傳播到應(yīng)變片時(shí)由示波器記錄下應(yīng)變隨時(shí)間的變化曲線，在通過(guò)試件以后，反射波和透射波再次由示波器記錄下信號(hào)，通過(guò)對(duì)入射波、反射波和透射波處理后，即能得到?jīng)_擊荷載下的應(yīng)力－應(yīng)變曲線.

圖1 典型分離式霍普金森壓桿試驗(yàn)裝置Fig.1 Typical SHPB experimental device

2 細(xì)觀數(shù)值模型

根據(jù)試驗(yàn)中裝置的尺寸建立有限元模型，入射桿和透射桿直徑均為74 mm，為了減小計(jì)算量，提高計(jì)算速度，在靠近試件附近的單元?jiǎng)澐州^小，其他位置劃分較大，如圖2所示，且模型中只建立了1/2模型.分別對(duì)0.0%、0.4%、0.6%和1.0%4種體積率的鋼纖維混凝土進(jìn)行建模，骨料的位置按照蒙特卡洛(Monte Carlo)方法生成，并使得4種模型骨料分布一致.試件的模型包括骨料、砂漿、界面過(guò)渡區(qū)(ITZ)以及鋼纖維，圖2所示的是鋼纖維體積率為1.0%時(shí)的細(xì)觀模型.

圖2 混凝土試件的軸對(duì)稱有限元模型Fig.2 Axisymmetric finite element model of concrete specimen

3 材料模型及參數(shù)

撞擊桿、入射桿和透射桿由于均處于線彈性范圍內(nèi)，因此采用線彈性材料模型［9］，其密度ρ=7 800 kg·m－3，彈性模量E=210 GPa，泊松比 ν =0.3.

砂漿和界面過(guò)渡區(qū)(ITZ)均是非均勻、多孔且強(qiáng)度低的材料，因此一般采用多孔材料模型來(lái)描述這兩種材料，多孔材料在壓力p下的體積剛度α(p，e)可以寫為

其中:αp是原始塑性屈服值;p為當(dāng)前壓力;pe為初始?jí)毫?ps為完全壓實(shí)時(shí)的壓力;n=3，為壓縮指數(shù).

骨料由于強(qiáng)度高，相對(duì)變形小，因此將其看作是小變形的脆性材料模型，選擇線性狀態(tài)方程來(lái)描述，即:

其中:p為壓力;τ=(ρ/ρ0)－1，ρ0是材料的初始密度;ρ表示在壓力p作用下的密度;k為材料的彈性模量.

上述材料模型的主要參數(shù)見表1［9－10］.

表1 材料參數(shù)模型Table 1 Material properties

4 損傷失效模型

圖3為模擬中選擇的材料損傷失效模型，具體的破壞過(guò)程主要分為3個(gè)階段:彈性階段;塑性強(qiáng)化階段;材料的軟化階段.在Pt.1以前的彈性階段表示應(yīng)力并未達(dá)到初始屈服面.Pt.1到Pt.2為塑性強(qiáng)化階段，表示材料的應(yīng)力達(dá)到初始屈服面而未達(dá)到極限強(qiáng)度面;在Pt.2到Pt.3為軟化階段，表示材料的應(yīng)力達(dá)到極限強(qiáng)度面而未達(dá)到失效強(qiáng)度破壞面.

圖3 材料的損傷失效模型Fig.3 Failure surface model with damage and failure

如圖3所示，極限強(qiáng)度面與失效強(qiáng)度面均采用三參數(shù)壓力函數(shù)形式，其極限強(qiáng)度面方程為

其中:a0、a1和a2為方程的3個(gè)參數(shù)，p=－(σxx+σyy+σzz)/3為壓力.

失效破壞強(qiáng)度面方程為

其中:a0f、a1f和a2f為方程的3個(gè)參數(shù).

屈服面位于極限強(qiáng)度面和失效強(qiáng)度之間，則

式中，參數(shù)η是λ的函數(shù)，函數(shù)η(λ)從初始值增加到1，表示從初始條件達(dá)到極限強(qiáng)度面，從1減小到0表示軟化階段.函數(shù)λ通過(guò)下面的方程來(lái)描述:

5 數(shù)值模擬結(jié)果與分析

圖4分別為子彈速度5.2 m·s－1時(shí)，鋼纖維體積率在0%、0.4%、0.6%和1.0%時(shí)各個(gè)不同時(shí)刻下SFRC試件的破壞過(guò)程，4種體積率的鋼纖維混凝土分別用A、B、C和D表示.從圖4(A)素混凝土的破壞過(guò)程來(lái)看，素混凝土受到?jīng)_擊加載后最先破壞，在316 μs就已經(jīng)破壞嚴(yán)重，而此時(shí)含有鋼纖維的混凝土并未破壞.當(dāng)增加鋼纖維含量為0.4%時(shí)，試件在325 μs時(shí)才出現(xiàn)裂紋.對(duì)比B、C和D可見，隨著鋼纖維體積率的不斷增加，裂紋開展的時(shí)間越延遲，表現(xiàn)出越好的延性和韌性.從各個(gè)破壞圖來(lái)看，破壞面主要沿著各個(gè)面上的鋼纖維，試件仍保留成塊狀，在裂縫處仍然有大量鋼纖維連接.當(dāng)試件破壞時(shí)，鋼纖維仍然能夠連接，鋼纖維本身并未被拉斷.

圖4 不同體積率下SFRC試件的破壞變形對(duì)比圖Fig.4 Comparison of failure modes of SFRC specimens with different fiber dosage

5 結(jié)論

利用有限元軟件LS-DYNA對(duì)不同體積率下的鋼纖維混凝土SHPB沖擊壓縮試驗(yàn)進(jìn)行了細(xì)觀數(shù)值模擬，模型中包括了骨料、砂漿和界面過(guò)渡區(qū)以及鋼纖維.通過(guò)模擬再現(xiàn)了鋼纖維混凝土在沖擊過(guò)程中的裂紋產(chǎn)生、擴(kuò)展的破壞過(guò)程，隨著鋼纖維體積率的增加，試件裂紋開展的時(shí)間延后，表現(xiàn)出良好的延性和韌性，且大量裂縫出現(xiàn)在鋼纖維與混凝土的連接處，具有明顯的增韌作用.

［1］ ASHOUR S A，HASANAIN G S，WAFA F F.Shear behavior of high-strength fiber reinforced concrete beams［J］.ACI Struct J，1992，89(2):176-184.

［2］ LI Z J，MU B，CHANG T Y P，et al.Prediction of overall tension behavior of short fiber-reinforced composites［J］.In J Solid Struct，1999，36(27):4071-4087.

［3］ 張亞芳，李琳.不同尺寸鋼纖維混雜增強(qiáng)混凝土彎曲性能數(shù)值分析［J］.廣州大學(xué)學(xué)報(bào):自然科學(xué)版，2012，11(1):49-52.ZHANG Y F，LI L.Numerical analysis on reinforced concrete mixed with steel fibers in varying geometry sizes［J］.J Guangzhou Univ:Nat Sci Edi，2012，11(1):49-52.

［4］ RAO T D G，SESHU D R.Torsion of steel fiber reinforced concrete members［J］.Cement Concr Res，2003，33(11):1783-1788.

［5］ 馬懷發(fā)，陳厚群，黎保琨.混凝土試件細(xì)觀結(jié)構(gòu)的數(shù)值模擬［J］.水利學(xué)報(bào)，2004，35(10):27-35.MA H F，CHEN H Q，LI B K.Meso-structure numerical simulation of concrete specimens［J］.J Hydraul Engin，2004，35(10):27-35.

［6］ 唐欣薇，秦川，張楚漢.基于細(xì)觀力學(xué)的混凝土類材料破損分析［M］.北京:中國(guó)建筑工業(yè)出版社，2012.TANG X W，QIN C，ZHANG C H.Analysis on damage of concrete-like materials based on meso mechanics［M］.Beijing:China Architecture＆ Building Press，2012.

［7］ 黨發(fā)寧，田威，鄭婭娜，等.基于細(xì)觀損傷的混凝土破裂過(guò)程的三維數(shù)值模擬［J］.工業(yè)建筑，2006，36(S):887-891.DANG F N，TIAN W，ZHENG Y N，et al.3-D numerical simulation of concrete failure process based on meso-mechanics damage［J］.Indust Constr，2006，36(S):887-891.

［8］ DU X L，JIN L，MA G W.Numerical simulation of dynamic tensile-failure of concrete at meso-scale［J］.Int J Impact Engin，2014，66(4):5-17.

［9］ HAO Y F，HAO H，LI Z X.Influence of end friction confinement on impact tests of concrete material at high strain rate［J］.Int J Impact Engin，2013，60(10):82-106.

［10］HAO Y F，HAO H.Numerical evaluation of the influence of aggregates on concrete compressive strength at high strain rate［J］.Int J Prot Struct，2011，2(2):177-206.