不同體積摻量的鋼纖維在混凝土中的分布研究

畢繼紅，趙云，鮑春，霍琳穎，王照耀，喬浩玥

（1.天津大學(xué)建筑工程學(xué)院，天津 300072；2.天津大學(xué)濱海土木工程結(jié)構(gòu)與安全教育部重點(diǎn)實(shí)驗(yàn)室，天津 300072；3.天津市市政工程設(shè)計(jì)研究院，天津 300392）

關(guān)健詞：纖維體積摻量；纖維分布；纖維取向；混凝土的流變性；數(shù)值模擬

眾所周知，混凝土的主要缺點(diǎn)是抗拉強(qiáng)度低，一旦發(fā)生開裂，在拉應(yīng)力的作用下會(huì)急劇失去其承載能力.這種弱點(diǎn)在很大程度上限制了普通混凝土在土木工程中的應(yīng)用.彌補(bǔ)此缺陷的一種有效方法是將適量的鋼纖維摻入混凝土中.纖維通過抵抗裂縫的形成和傳播，提升了水泥基材料的延性、韌性和耗能能力[1-2].由于易于澆筑，鋼纖維混凝土（SFRC）在不同的結(jié)構(gòu)應(yīng)用（例如工業(yè)地板，外部人行道，噴射混凝土和預(yù)制構(gòu)件）中越來越多地替代了鋼筋混凝土[3].如今，SFRC技術(shù)發(fā)展的新趨勢是將自密實(shí)混凝土和SFRC相結(jié)合，后者被稱為鋼纖維自密實(shí)混凝土（SFRSCC）.這種新型纖維水泥復(fù)合材料通過將自密實(shí)混凝土作為介質(zhì)來輸送鋼纖維，為混凝土結(jié)構(gòu)的設(shè)計(jì)和創(chuàng)新提供靈活性和便利性.

在過去的幾十年中，許多學(xué)者致力于研究纖維特性對SFRC力學(xué)性能的影響.纖維特性對SFRC的抗拉性能有著顯著的影響，例如纖維類型、長徑比、體積摻量以及分布取向等，具有不同纖維參數(shù)的SFRC表現(xiàn)出不同的力學(xué)性能[4].在纖維的各項(xiàng)特性中，增加纖維體積摻量是改善混凝土拉伸或彎曲性能最有效的方法.在探究纖維體積摻量對SFRC的力學(xué)性能的影響時(shí)，大多數(shù)學(xué)者通常保持混凝土基體的配合比不變[5]，將纖維體積摻量作為唯一變量，而纖維體積摻量對新拌混凝土的流變性質(zhì)有著重要的影響[6]，流變性質(zhì)直接影響纖維的分布與取向[7].然而，已有研究往往會(huì)忽略因纖維摻量改變而引起的纖維分布與取向變化，但纖維的分布和取向?qū)FRC的力學(xué)性能有著至關(guān)重要的影響[8].因此，通過試驗(yàn)或數(shù)值模擬方法探究纖維分布與取向隨纖維體積摻量的變化規(guī)律具有重要意義.

目前，用于檢測硬化后混凝土基體內(nèi)部纖維分布與取向的方法包括X射線法、電磁感應(yīng)法[9]、圖像處理技術(shù)以及電阻率測量法[10]等.但是，這些方法需要復(fù)雜的數(shù)據(jù)分析過程，并且只能確定試件某個(gè)切割面上的纖維分布形態(tài)，很難應(yīng)用到實(shí)際的工程結(jié)構(gòu)中[11].數(shù)值模擬允許觀察混凝土流動(dòng)過程中纖維取向的變化，有助于深入理解新拌混凝土中的纖維取向與力學(xué)性能之間的關(guān)系，但目前依然缺乏關(guān)于新拌混凝土流變性質(zhì)對纖維分布和取向的影響的研究.

為此，本文配制了4種不同纖維體積摻量的纖維自密實(shí)混凝土，分析了纖維體積摻量與新拌混凝土的流變性質(zhì)之間的關(guān)系.借助軟件ANSYS CFX模擬具有不同流變性質(zhì)的混凝土的澆筑過程，基于模擬得到的混凝土速度場，將纖維視為剛性連接的多個(gè)球形粒子，采用剛體運(yùn)動(dòng)動(dòng)力學(xué)方程計(jì)算了纖維的運(yùn)動(dòng)，分析了纖維體積摻量對纖維分布與取向的影響規(guī)律.本文研究結(jié)果可為優(yōu)化鋼纖維混凝土的配合比提供相關(guān)依據(jù).

1 試驗(yàn)工況

1.1 試件設(shè)計(jì)與制作

試驗(yàn)共設(shè)計(jì)了4組12個(gè)梁試件，共選取4種鋼纖維體積摻量（0.4%、0.6%、0.8%、1%），每組包含3個(gè)試件，分別記為SF-0.4、SF-0.6、SF-0.8、SF-1.0，例如，SF-0.6表示鋼纖維摻量為0.6%的梁試件.試件為400 mm×100 mm×100 mm的棱柱體.

在只改變纖維體積摻量的情況下，保持混凝土基體成分比例不變，混凝土基體的組成如表1所示.制作各試件所用的混凝土水灰比均為0.45，水泥標(biāo)號(hào)為P.O 42.5R，骨料由細(xì)度模數(shù)為2.3的細(xì)砂與粒徑為5～12 mm的碎石組成.鋼纖維采用長度為30 mm、截面直徑為0.75 mm的長直纖維，抗拉強(qiáng)度為1 050 MPa.為保證鋼纖維均勻分布在混凝土中，首先對粗細(xì)骨料進(jìn)行干拌，邊攪拌邊加入水泥及礦物參合料，再分批加入鋼纖維進(jìn)行拌和，加入80%的水和減水劑，攪拌數(shù)分鐘后加入剩余的水，直到所有混合料拌和均勻后停止攪拌.

表1 混凝土基體配合比Tab.1 Proportion of concrete matrix kg·m-3

評價(jià)自密實(shí)混凝土的工作性能最常用的標(biāo)準(zhǔn)試驗(yàn)有V形漏斗試驗(yàn)、L形箱試驗(yàn)、U形箱試驗(yàn)、J環(huán)試驗(yàn)和坍落擴(kuò)展度流動(dòng)試驗(yàn).本文選用以上標(biāo)準(zhǔn)試驗(yàn)中的坍落擴(kuò)展度試驗(yàn)與J環(huán)試驗(yàn)評價(jià)新拌混凝土的工作性能.

從模板的中心位置澆筑已拌好的混凝土，澆筑完成并在實(shí)驗(yàn)室環(huán)境中養(yǎng)護(hù)48 h后，拆除模板并在溫度為20±2 ℃、相對濕度大于90%的環(huán)境下養(yǎng)護(hù)28 d.養(yǎng)護(hù)完成后，對摻有纖維的各組試件進(jìn)行切割，統(tǒng)計(jì)各切割面纖維數(shù)目以備試驗(yàn).

1.2 梁切割面上纖維數(shù)目的計(jì)算

沿每組梁試件的三個(gè)正交方向進(jìn)行切割，纖維斷面顏色與混凝土基體的顏色具有明顯的差別，故通過人工計(jì)數(shù)的方法統(tǒng)計(jì)梁的三個(gè)正交切割面上的纖維根數(shù)，同時(shí)對纖維根數(shù)的計(jì)數(shù)過程重復(fù)數(shù)遍以排除由于測量導(dǎo)致的誤差.梁三個(gè)正交切割面分別為垂直于梁軸向的橫斷面（α）、垂直于澆筑方向的豎斷面（γ）以及平行于澆筑方向的縱斷面（β），如圖1所示.此外，圖2給出了切割面上的鋼纖維分布情況.

圖1 梁的三個(gè)切割面α，β 和γFig.1 The cutting planes α，β and γ of beam specimen

圖2 切割面α 上的鋼纖維分布Fig.2 Fiber distribution in the cutting planes α

采用式（1）計(jì)算每個(gè)切割面上的纖維密度：

式中：nf為切割面上纖維的總根數(shù)；Ac為切割面的面積.各體積摻量下的纖維密度值取每組試件的平均值[12]，以體現(xiàn)試驗(yàn)結(jié)果的統(tǒng)計(jì)特征.

將切割面分為數(shù)個(gè)大小相等的正方形單元，計(jì)算每個(gè)單元中的纖維根數(shù)，采用纖維分散系數(shù)評價(jià)纖維的分散程度[13]，計(jì)算公式如下：

式中：xi代表切割面上第i個(gè)分割單元中纖維的根數(shù).αf值越接近于1，表明纖維在切割面上的分布越均勻.

2 試驗(yàn)結(jié)果

2.1 混凝土的流變性分析

圖3給出了坍落擴(kuò)展度和J環(huán)測試得到的新拌SFRSCC的流動(dòng)結(jié)果.表2給出了各組新拌混凝土坍落擴(kuò)展度、J環(huán)擴(kuò)展度以及T500（坍落至擴(kuò)展直徑為500 mm所需的時(shí)間）.由表2可知，鋼纖維的加入對自密實(shí)混凝土的T500和總流動(dòng)時(shí)間具有顯著的影響.坍落擴(kuò)展度隨著鋼纖維的體積摻量增加呈現(xiàn)逐漸減小的趨勢，這是因?yàn)殇摾w維的加入增加了混凝土的粘稠性，自密實(shí)混凝土的流動(dòng)性能減弱，致使其擴(kuò)展半徑減小.此外，相同鋼纖維體積摻量下的坍落擴(kuò)展度與J環(huán)擴(kuò)展度的差值較小，表明各體積摻量下的鋼纖維自密實(shí)混凝土具有良好的間隙通過能力.總體上，各組SFRSCC均符合自密實(shí)混凝土的流動(dòng)能力和通過能力的基本要求.

圖3 新拌SFRSCC的工作性能測試Fig.3 Workability tests for fresh SFRSCC

流變儀可以用來測量新拌混凝土的塑性黏度和屈服應(yīng)力，但是對于相同的新拌混凝土，不同的流變儀給出不同的屈服應(yīng)力和塑性黏度測量值.因此，本文采用Ghanbari等[14]提出的黏度模型計(jì)算各組新拌混凝土的塑性黏度.

表2 不同纖維體積摻量下SFRSCC的工作性能Tab.2 Workability of fresh concrete for SFRSCC with different fiber volume contents

該模型考慮了鋼纖維對混合物塑性黏度的影響，并已通過大量實(shí)驗(yàn)數(shù)據(jù)進(jìn)行了驗(yàn)證.預(yù)測塑性黏度的公式為：

式中：μ 為沒有鋼纖維的新拌混凝土的塑性黏度；vf為鋼纖維體積摻量；re為鋼纖維的長徑比.

使用Roussel等[15]描述的分析方法計(jì)算了新拌混凝土的屈服應(yīng)力.該方法建立了坍落度試驗(yàn)與屈服應(yīng)力之間的關(guān)系，屈服應(yīng)力的計(jì)算公式為：

式中：ρ 為新拌混凝土的密度；H為新拌混凝土流動(dòng)結(jié)束后的高度；R為新拌混凝土停止流動(dòng)時(shí)半徑.

基于混凝土工作性能試驗(yàn)結(jié)果，各組新拌混凝土的塑性黏度和屈服應(yīng)力均由式（3）和式（4）計(jì)算得到.表3列出了各組新拌混凝土的塑性黏度和屈服應(yīng)力值.新拌混凝土的流變性質(zhì)與纖維的體積摻量之間的關(guān)系如圖4所示.結(jié)果表明，隨著纖維體積摻量的增加，塑性黏度和屈服應(yīng)力均呈線性增加.Alberti等[16]同樣發(fā)現(xiàn)：增加纖維體積摻量會(huì)導(dǎo)致新拌混凝土的屈服應(yīng)力和塑性黏度升高；流變參數(shù)與纖維的體積摻量之間的關(guān)系是線性的.

表3 各組新拌混凝土的流變性質(zhì)Tab.3 Rheological properties for each type of fresh concrete

圖4 新拌混凝土塑性黏度和屈服應(yīng)力與纖維體積摻量之間的關(guān)系Fig.4 Relationship between plastic viscosity，yield stress of fresh concrete and fiber content

2.2 纖維計(jì)數(shù)分析

表4給出了梁試件的三個(gè)正交切割面測得的平均纖維密度.比較三個(gè)正交切割面上的纖維密度，可以發(fā)現(xiàn)平面α 和平面γ 中的纖維密度相似，均高于平面β 中的纖維密度，表明纖維趨向于分布在平面α 和平面γ 上，這是由于在澆筑過程中混凝土主要沿著水平面流動(dòng)，而纖維的取向逐漸變得與流動(dòng)方向一致.此外，混凝土澆筑點(diǎn)位于模板中心，會(huì)導(dǎo)致α 平面中的纖維密度較大.

表4 各切割面上的平均纖維密度Tab.4 Average fiber density in the cut sections cm2

表5給出了不同體積摻量下切割面α 上的分散系數(shù)值.隨著鋼纖維體積摻量的提高，切割面α 上的分散系數(shù)逐漸增大，表明鋼纖維在切割面上的分布越來越均勻，這可歸因于纖維摻量的增加改變了混凝土的流變性，進(jìn)而導(dǎo)致纖維分布趨于均勻.

表5 切割面α 上的纖維分散系數(shù)值Tab.5 Dispersion coefficient in the cut section α

3 基本理論

3.1 新拌混凝土的本構(gòu)方程

自密實(shí)混凝土被視為非牛頓流體，其流變行為通常由Bingham或Herschel-Bulkley黏度模型來表征.本文選用Bingham黏度模型用于模擬混凝土的流動(dòng)，因?yàn)樵撃Ｐ偷膮?shù)具有物理意義，并且比Herschel-Bulkley模型的參數(shù)更容易確定.Bingham的本構(gòu)方程如下：

式中：τ 為剪應(yīng)力；τy為屈服應(yīng)力；μ 為塑性黏度；γ 為剪應(yīng)變率.

3.2 新拌混凝土的流動(dòng)

新拌混凝土的流動(dòng)模型基于均質(zhì)方法，即將新拌混凝土視為一種均質(zhì)的流體.采用CFD軟件CFX模擬新拌混凝土的流動(dòng)，并求解流體質(zhì)量和動(dòng)量守恒的基本控制方程.基本控制方程為：

式中：u為流體速度；ρ 為流體密度；p為壓強(qiáng)；τc為剪切力張量；g為重力加速度.

混凝土和空氣之間的界面隨時(shí)間而變化，因?yàn)樾掳杌炷恋牧鲃?dòng)是一個(gè)瞬態(tài)過程.采用兩相流VOF方法追蹤混凝土的自由表面，VOF方法的基本思想是計(jì)算和更新各網(wǎng)格中混凝土和空氣各占的體積摻量.采用有限體積法（FVM）將整個(gè)計(jì)算域劃分為網(wǎng)格，每個(gè)網(wǎng)格節(jié)點(diǎn)由給定大小的控制體積圍繞，通過在每個(gè)控制體積上積分Navier-Stokes方程，從而獲得一系列離散方程.求解這些離散方程以計(jì)算混凝土的速度場，速度場用于更新空氣和混凝土所占的體積摻量.

3.3 鋼纖維的處理

如果每根單獨(dú)的纖維被視為單個(gè)球形顆粒，通過計(jì)算顆粒軌跡只能監(jiān)測流動(dòng)過程中該纖維的平動(dòng)，而無法捕獲纖維的取向.為了監(jiān)測纖維在流動(dòng)過程中的取向，Kulasegaram等[17]使用一對距離等于纖維長度的球形顆粒代表單根纖維.實(shí)際上，沿纖維長度的每個(gè)部分在流動(dòng)過程中都會(huì)受到流體的作用力，因此，將纖維的運(yùn)動(dòng)簡化為剛性連接的多個(gè)球形粒子的運(yùn)動(dòng)，能夠更精確地考慮纖維的轉(zhuǎn)動(dòng)與平動(dòng).圖5給出了纖維的簡化方法.

圖5 纖維的簡化方法Fig.5 The numerical model of beam

纖維被簡化為多個(gè)剛性連接的球形粒子，纖維質(zhì)量由各粒子平均分配.因此，可將浸入混凝土中的纖維視為剛性細(xì)長桿，在流動(dòng)過程中僅產(chǎn)生剛性運(yùn)動(dòng)，即轉(zhuǎn)動(dòng)和平動(dòng).選取纖維的中心作為基本點(diǎn)，將纖維的運(yùn)動(dòng)分解為纖維相對于基點(diǎn)的旋轉(zhuǎn)和基點(diǎn)的平移，采用剛體的一般動(dòng)力學(xué)方程計(jì)算纖維的運(yùn)動(dòng).纖維的平動(dòng)方程為：

式中：mf為單根纖維的質(zhì)量；vf為纖維中點(diǎn)的速度；FD，i為第i個(gè)球形粒子的拖曳力；Fg為纖維受到的廣義重力（包括浮力）.

計(jì)算球形粒子所受拖曳力的公式為：

式中：CD為粒子的拖曳力系數(shù)；Ap為粒子在流動(dòng)方向的投影面積；up為與粒子相同位置處的新拌混凝土速度；v為粒子的速度.需要注意的是，混凝土的速度是基于CFX軟件中混凝土的速度場計(jì)算得到的，而CD由Haider等[18]提出的非球形粒子的阻力系數(shù)公式計(jì)算得到.

纖維的廣義重力的計(jì)算公式如下：

纖維的轉(zhuǎn)動(dòng)方程為：

式中：I為纖維的慣性張量；ω 為纖維的角速度；T為各粒子拖曳力對纖維中心的力矩之和.

T的計(jì)算公式如下：

式中：di為第i個(gè)球形粒子中心與纖維中心之間的距離.

根據(jù)上述計(jì)算過程，提出一種與CFD相結(jié)合的數(shù)值方法，用于監(jiān)測計(jì)算域中所有纖維的位置.考慮到計(jì)算時(shí)間與計(jì)算精度的平衡性，將每根纖維分為10個(gè)球形顆粒.

4 數(shù)值方法驗(yàn)證

為了驗(yàn)證上述數(shù)值方法模擬纖維在自密實(shí)混凝土中運(yùn)動(dòng)的可行性，以本文試驗(yàn)中新拌混凝土的澆筑過程為對象進(jìn)行驗(yàn)證.

4.1 模型的建立

按照試件的實(shí)際尺寸，建立混凝土梁的數(shù)值模型，如圖6所示.為便于模擬結(jié)果的分析，對模型的三個(gè)軸進(jìn)行統(tǒng)一的方向標(biāo)定，Z軸為垂直于澆筑面的軸，X軸平行于梁模型主軸方向，Y軸為平行于試件短邊的軸.采用結(jié)構(gòu)化網(wǎng)格劃分方法，進(jìn)行靈敏度分析，六面體網(wǎng)格邊長確定為10 mm，網(wǎng)格總數(shù)為5 816個(gè)，基本實(shí)現(xiàn)計(jì)算精度和效率的平衡.

圖6 梁的數(shù)值模型Fig.6 The numerical model of beam

4.2 邊界條件及參數(shù)設(shè)置

各組新拌混凝土的屈服應(yīng)力和塑性黏度采用表3中的數(shù)值，密度為2 255.3 kg/m3，空氣視為密度為1.185 kg/m3、黏度為1.831×10-5的牛頓流體.根據(jù)實(shí)驗(yàn)中混凝土澆筑點(diǎn)的位置，入口的邊界條件采用速度邊條，設(shè)置在模型上表面正中心處，尺寸為40 mm×40 mm，澆筑速度為0.2 m/s，方向垂直于模型上表面.出口邊界條件采用開放壓力邊條（Opening），允許空氣進(jìn)入和離開計(jì)算域，計(jì)算域內(nèi)外壓力都設(shè)置為1個(gè)大氣壓，內(nèi)外壓力的差值為0.由于混凝土可能在壁面附近滑動(dòng)，因此對模板壁面施加有限滑移邊界條件.時(shí)間步長為0.01 s，計(jì)算總步數(shù)為1 250步，采用二階隱式求解器，用迭代的方法計(jì)算流體質(zhì)量和動(dòng)量守恒的基本控制方程，以此來達(dá)到精確模擬的效果.

4.3 纖維的生成

使用Mersenne-Twister程序在澆筑口正上方生成隨機(jī)分布的纖維，該程序能生成均勻的偽隨機(jī)值，從而纖維的位置和取向都是隨機(jī)的.

各組新拌混凝土中纖維的總根數(shù)由式（13）計(jì)算得到：

式中：V為梁試件的體積；Vfib為單根纖維的體積.

在數(shù)值模擬過程中，纖維體積摻量為0.4%、0.6%、0.8%和1.0%對應(yīng)的纖維總根數(shù)分別為1 207、1 811、2 414和3 018.由于纖維摻量的變化，每種模擬工況中計(jì)算域纖維總根數(shù)是不同的，導(dǎo)致難以對比不同混凝土配合比下纖維的取向.為了方便比較不同纖維體積摻量對應(yīng)的纖維取向分布，采用纖維取向系數(shù)[19]分析纖維與梁主軸（X軸）之間的夾角分布.取向系數(shù)計(jì)算公式為：

式中：θi為纖維與梁主軸之間的夾角.

4.4 纖維運(yùn)動(dòng)的求解流程

圖7給出了計(jì)算纖維分布和取向的流程圖.首先，采用CFX獲得新拌混凝土在每個(gè)時(shí)間步的速度場.隨后，在計(jì)算域的入口上方生成隨機(jī)分布的纖維，纖維根據(jù)混凝土澆筑的總步數(shù)分批進(jìn)入計(jì)算域中.需要注意的是，每一時(shí)間步長中計(jì)算的纖維總數(shù)由上一步已經(jīng)進(jìn)入域中的纖維和本時(shí)間步內(nèi)新進(jìn)入域中的纖維組成.基于CFX計(jì)算得到的速度場，計(jì)算每步進(jìn)入域中纖維的運(yùn)動(dòng)，以確定纖維的位置和取向.在下一時(shí)間步中，本步纖維的位置和取向作為已知值，進(jìn)行循環(huán)運(yùn)算.最后，當(dāng)步數(shù)達(dá)到1 250步，也就是混凝土澆筑完成時(shí)，所有纖維進(jìn)入計(jì)算域，循環(huán)結(jié)束.在Matlab中通過編程實(shí)現(xiàn)上述過程，時(shí)間步長同樣為0.01 s.

4.5 驗(yàn)證結(jié)果

模擬計(jì)算得到的纖維密度列于表6.圖8對比了試驗(yàn)和模擬獲得的纖維分散系數(shù)值.結(jié)果表明，通過數(shù)值模擬得到的纖維密度和纖維分散系數(shù)均與試驗(yàn)結(jié)果吻合較好，說明提出的數(shù)值方法能夠模擬新拌混凝土流動(dòng)過程中纖維的運(yùn)動(dòng).

圖7 計(jì)算纖維分布與取向的流程圖Fig.7 Flow chart of calculating the motion of fibers

表6 模擬得到的各切割面纖維密度值Tab.6 Fiber density in the cut sections obtained by simulation

圖8 試驗(yàn)與數(shù)值模擬計(jì)算得到的纖維分散系數(shù)值Fig.8 Fiber dispersion coefficient calculated by experiment and numerical simulation

5 模擬結(jié)果及分析

5.1 纖維體積摻量對纖維分布的影響

首先將梁的數(shù)值模型在豎直方向（Z軸）上分為兩部分，即下半梁和上半梁.計(jì)算每部分的纖維根數(shù)，除以纖維總根數(shù)，通過纖維的相對含量分析纖維的豎向分布.圖9給出了不同纖維體積摻量下豎直方向上纖維的相對含量.結(jié)果表明，梁下半部中的纖維相對含量多于梁上半部中的纖維相對含量，說明在混凝土澆筑過程中纖維發(fā)生了離析.纖維離析主要是由于重力作用引起的.另外，梁的上半部分和下半部分之間的纖維分布差異隨著纖維體積摻量的增加而減小，兩部分之間的纖維分布的差異值從38.4%減少到27.4%，這主要是因?yàn)樵黾永w維體積摻量會(huì)導(dǎo)致新鮮混合物的塑性黏度和屈服應(yīng)力的增大，纖維與混凝土之間的黏聚力隨之增強(qiáng)，進(jìn)而減弱了纖維的離析程度.Hosseinpoor等[20]得出了類似的結(jié)論：增加新拌混凝土的塑性黏度和屈服應(yīng)力會(huì)減少骨料在垂直方向上從水泥漿中的離析.

圖9 不同纖維體積摻量下豎直方向上纖維的相對含量Fig.9 Relative fiber content in vertical direction under different fiber volume content

圖10給出了不同纖維體積摻量下纖維在水平方向上的分布.結(jié)果表明，纖維相對含量在梁跨中分布較多，在梁左右兩端分布對稱.注意到有一個(gè)增加的趨勢：越接近混凝土澆筑點(diǎn)，纖維數(shù)量越多，這種趨勢主要由澆筑口的位置造成的，澆筑口下方纖維發(fā)生沉積，因而該處纖維數(shù)量較多.

5.2 纖維體積摻量對纖維取向的影響

圖10 不同纖維體積摻量下水平方向上纖維的相對含量Fig.10 Relative fiber content in horizontal direction under different fiber volume contents

數(shù)值模擬得到的混凝土瞬態(tài)流動(dòng)和纖維分布形態(tài)見圖11和圖12，此處僅給出了對應(yīng)于0.4%和1.0%的纖維體積摻量的模擬結(jié)果.由圖11和圖12可知，鋼纖維被混凝土完全包裹.纖維體積摻量為1.0%的混凝土界面的分層現(xiàn)象更為明顯，因?yàn)槔w維體積摻量較高的混凝土具有更大的黏度和屈服應(yīng)力.另外，越靠近壁面，纖維越趨于與壁面平行分布，這種取向趨勢主要由混凝土的流速引起.圖13顯示了距離梁底部0.02 m處水平面上混凝土速度的云圖，可以看到混凝土的速度從壁面處最小值0 m/s組件增大到平面中心處的最大值.因此，由于混凝土在水平面上存在速度差，壁面附近纖維逐漸與壁面對齊.

圖11 兩個(gè)不同時(shí)刻時(shí)纖維體積摻量為0.4%的新拌混凝土的流動(dòng)與纖維分布形態(tài)Fig.11 The flow of fresh concrete and the fiber distribution for the fiber volume content of 0.4%at two instants

圖12 兩個(gè)不同時(shí)刻時(shí)纖維體積摻量為1.0%的新拌混凝土的流動(dòng)與纖維分布形態(tài)Fig.12 The flow of fresh concrete and the fiber distribution for the fiber volume content of 1.0%at two instants

圖13 距離梁底部0.02 m處水平面上的混凝土速度云圖Fig.13 Concrete velocity contour of the longitudinal section at 0.02 m from the bottom of the beam

不同纖維體積摻量下纖維取向角的概率密度分布如圖14所示，圖中纖維取向角是基于域中所有纖維根數(shù)計(jì)算得到的，而不是針對某個(gè)切割面.由圖14可知，纖維體積摻量為0.4%對應(yīng)的取向分布明顯偏左，而纖維體積摻量為1.0%對應(yīng)的取向直方圖偏右分布，說明纖維取向角的平均值隨纖維體積摻量的增加而增加，這同樣是因?yàn)樾掳杌炷恋牧髯冃再|(zhì)改變而導(dǎo)致的，即具有較低塑性黏度和屈服應(yīng)力的新拌混凝土具有更好的流動(dòng)性，纖維取向角在流場的作用下變得更小.另一方面，纖維體積摻量為0.6%和0.8%對應(yīng)的取向直方圖均居中，但摻量為0.8%對應(yīng)的取向分布比摻量為0.6%的略微偏右一些.應(yīng)該注意的是，角度分布在70°～90°之間的概率稍大于文獻(xiàn)[21]對應(yīng)的結(jié)果.產(chǎn)生差異主要有兩個(gè)原因：澆筑方法和試件的幾何尺寸.與從梁端部澆筑混凝土的結(jié)果相比，從梁中點(diǎn)位置澆筑混凝土?xí)p少纖維的流動(dòng)距離.另一方面，試件的寬度（Y軸方向）比纖維的長度大許多，纖維能夠自由旋轉(zhuǎn).以上兩點(diǎn)原因都可能增大纖維趨近于大角度的概率.

圖14 不同纖維體積摻量下纖維取向角的概率密度分布Fig.14 Probability density distribution of fiber orientation angle under different fiber volume contents

圖15給出了纖維取向系數(shù)沿梁主軸（X軸）的變化.由圖15可知，跨中位置切割面上的纖維具有較大的取向系數(shù)，意味著纖維趨于與梁主軸平行分布.相反，梁兩端的纖維均勻分布并且取向系數(shù)較小，這主要?dú)w因于梁中混凝土的流場與“壁面”效應(yīng)，梁端部的壁面會(huì)使得端部區(qū)域的混凝土流動(dòng)方向?yàn)樽韵露希率勾藚^(qū)域的纖維與梁主軸之間的平均角度偏大，而梁跨中位置的混凝土水平地流向梁的兩端，導(dǎo)致纖維具有較大的取向系數(shù).另一方面，隨著纖維體積摻量的增加，各切割面上纖維系數(shù)逐漸減小.

圖15 纖維取向系數(shù)沿梁主軸的變化Fig.15 The change of fiber orientation factor along main axis of beam

6 結(jié)論

1）纖維體積摻量是影響新拌混凝土流變性質(zhì)的敏感因素，隨著纖維體積摻量的增加，新拌混凝土的屈服應(yīng)力和塑性黏度呈線性增加.

2）考慮纖維體積摻量對混凝土流變性質(zhì)的影響，采用CFX軟件模擬新拌混凝土澆筑過程.基于混凝土的速度場，利用Matlab編程求解纖維的運(yùn)動(dòng)方程，通過與試驗(yàn)結(jié)果對比證明了本文提出的模擬方法能夠較好地確定纖維的分布與取向.

3）由于重力作用，不同纖維體積摻量下纖維都發(fā)生了離析，并且纖維的離析程度隨著纖維體積摻量的增加而降低.相反，纖維在梁的水平方向上基本均勻分布，只是在澆筑口處發(fā)生輕微沉積.

4）纖維體積摻量是影響纖維取向的關(guān)鍵因素，纖維體積摻量越低，纖維與梁主軸之間傾斜角越小，梁各切割面上的纖維取向系數(shù)變大.