1、第 26 卷第 2 期2023 年 2 月建筑材料学报JOURNAL OF BUILDING MATERIALSVol.26，No.2Feb.，2023定向钢纤维增强水泥基复合材料断裂细观数值模拟卿龙邦*，杨子美，慕儒，张金鑫，籍晓凯（河北工业大学 土木与交通学院，天津 300401）摘要：采用随机生成算法投放钢纤维，建立了随机乱向、定向钢纤维增强水泥基复合材料（SFRC、ASFRC）三点弯曲梁细观有限元数值模型，计算了不同纤维掺量下 SFRC试件和 ASFRC试件加载断裂的全过程，分析了三点弯曲梁开裂截面处的纤维应力，研究了定向钢纤维的细观增强机理.结果表明：SFRC试件和 ASFRC试件荷
2、载-裂缝张开口位移全曲线的模拟值与试验值符合较好，峰值荷载的误差在 10%以内；SFRC试件和 ASFRC试件的峰值荷载与纤维合力的最大值均随着纤维掺量的增加而增大，当纤维掺量为 0.8%、1.2%、2.0%时，ASFRC 试件的峰值荷载较 SFRC 试件提高了 75%、111%、141%，纤维合力的最大值较 SFRC 试件增大了 202%、144%、127%；定向钢纤维可以有效改善水泥基复合材料的断裂性能，显著提高钢纤维的利用率，延缓裂缝的扩展.关键词：钢纤维；水泥基复合材料；细观数值模拟；等效应力-应变关系；增强机理中图分类号：TB332文献标志码：A doi：10.3969/j.issn
3、.1007-9629.2023.02.001Meso-Numerical Simulation on Fracture of Aligned Steel Fiber Reinforced Cementitious CompositesQING Longbang*，YANG Zimei，MU Ru，ZHANG Jinxin，JI Xiaokai（School of Civil and Transportation Engineering，Hebei University of Technology，Tianjin 300401，China）Abstract:A random generation
4、 algorithm was used to steel fiber，and then the meso-scale finite element numerical mode of three point bending beam of random and aligned steel fiber reinforced cementitious composites（SFRC，ASFRC）were established.The whole process of loading and fracture of SFRC specimens and ASFRC specimens with d
5、ifferent fiber contents were calculated.The fiber stress at the cracked section of the three point bending beam was analyzed，and the meso-scale reinforcement mechanism of aligned steel fiber was studied.The results show that the whole load-crack opening displacement curves of SFRC specimens and ASFR
6、C specimens obtained from test and numerical simulation are in good agreement and the error of the peak load is within 10%.The peak load and the maximum fiber force of the two kinds of specimens both increase with the increase of fiber content.When the fiber content is 0.8%，1.2%and 2.0%，the peak loa
7、d of ASFRC specimens increased by 75%，111%and 141%compared to SFRC specimens，and the maximum fiber force of ASFRC specimens increased by 202%，144%，and 127%compared with SFRC specimens.Aligned steel fiber can effectively improve the fracture performance of cementitious composites，significantly increa
8、se the fiber efficiency，and delay the propagation of cracks.Key words:steel fiber；cementitious composite；meso-numerical simulation；equivalent stress-strain relationship；enhancement mechanism文章编号：1007-9629（2023）02-0111-11收稿日期：2021-12-03；修订日期：2022-02-21基金项目：国家自然科学基金资助项目（52022027，52078180）；河北省自然科学基金资助项
9、目（E2022202141，E2020202151）；河北省高等学校科学技术研究重点项目（ZD2019072，ZD2020190）第一作者（通讯作者）：卿龙邦（1982），男，湖北天门人，河北工业大学教授，博士生导师，博士.E-mail：建筑材料学报第 26卷在水泥基体中加入适量的钢纤维可以有效抑制裂缝的扩展，改善水泥基复合材料的断裂性能1-3.研究表明，钢纤维的增强作用受其分布方向的影响4-5.在传统随机乱向钢纤维增强水泥基复合材料（SFRC）中，纤维呈随机乱向分布，部分纤维分布方向与拉应力方向并不一致，难以充分发挥其增强作用，将钢纤维按特定方向分布则可以明显提高其增强效率6.采用磁场法可以
10、实现钢纤维的定向分布.慕儒等7-9在盛有钢纤维增强水泥砂浆的试模外施加均匀电磁场，同时施加振动辅助，制备了定向钢纤维增强水泥基复合材料（ASFRC）并对其力学性能以及断裂性能进行了研究.Wijffels等10利用钢纤维在磁场中磁化的特性，制备了钢纤维定向分布的水泥基复合材料，研究了定向钢纤维对复合材料弯曲性能的影响.王哲伟等11通过对钢纤维增强地质聚合物施加均匀电磁场，制备了定向钢纤维增强地质聚合物预制缺口梁.汪洋等12采用磁场诱导定向技术，制备了定向碳纤维增强水泥砂浆，研究了碳纤维的取向性对碳纤维增强水泥砂浆力学性能提升效果的影响.细观数值模拟便于分析基体和钢纤维在断裂过程中任意时刻的受力状
11、态，有利于深入研究钢纤维的增强机理.Wang等 13-14 通过分析钢纤维混凝土材料的应力-应变曲线，建立了钢纤维混凝土材料的本构关系，采用最大主应变破坏准则模拟了钢纤维混凝土在动态压缩荷载及爆破荷载下的破坏过程.Soetens等 15 假定纤维与基体无相对滑动，将纤维单元嵌入基体单元，建立了三维钢纤维混凝土有限元模型，模拟了端钩形钢纤维混凝土四点弯曲梁开裂的全过程.Fang等 16 建立了可预测SFRC动态破坏的三维有限元模型，研究了高应变率下SFRC的动态压缩性能.Yu等 17 在预计产生裂缝的截面处插入黏聚单元，将与裂纹张开速率相关的内聚定律作为黏聚单元的失效准则，模拟了端钩形钢纤维混凝
12、土在不同落锤冲击速度下的断裂失效行为.本文采用随机生成算法投放钢纤维，建立了SFRC和 ASFRC三点弯曲梁细观有限元数值模型，计算了不同钢纤维掺量（体积分数）下 SFRC 试件和ASFRC 试件加载断裂的全过程，结合断裂试验验证了数值模型的有效性，分析了开裂截面处的纤维应力，研究了定向钢纤维的细观增强机理.1数值模拟1.1纤维投放在440 mm 100 mm 100 mm区域内投放直径（Df）为 0.5 mm，长度（Lf）为 30 mm 的圆直型钢纤维，钢纤维总数量（N）通过以下方法确定18：N=4VVfLfD2f（1）式中：V为投放区域的体积，mm3；Vf为钢纤维掺量，%.在SFRC试件中
13、，钢纤维呈随机分布.本文参考文献18 ，确定钢纤维掺量为0.8%，试件尺寸为440 mm100 mm100 mm，控制纤维投影与X轴的夹角（）分别为 0，30）、30，45）、45，60）、60，90 ，根据角度占比确定每种角度范围内的纤维数量Ni（i=1，2，3，4），分次投放钢纤维，直到投放的纤维总数量为N，建立不同角度范围内的随机乱向钢纤维细观数值模型，如图1所示.图 1不同角度范围内的随机乱向钢纤维细观数值模型Fig.1Meso-numerical model of random steel fibers in different angle ranges112第 2期卿龙邦，等：定向
14、钢纤维增强水泥基复合材料断裂细观数值模拟在 ASFRC试件中，钢纤维分布方向与主拉应力方向一致，基于随机乱向钢纤维的投放算法，确定钢纤维掺量为 0.8%，试件尺寸为 440 mm100 mm100 mm，控制纤维投影与 X轴夹角为 0，建立定向钢纤维细观数值模型，如图2所示.1.2建立模型1.2.1基体本构关系基体选用三维八节点实体单元建模，钢纤维选用三维二节点桁架单元建模.采用黏聚裂纹模型19计算基体的开裂.在三点弯曲梁有限元模型跨中处插入零厚度的黏聚单元，将牵引分离定律作为黏聚单元的失效准则20-22，以黏聚单元的失效来模拟试件裂缝的产生和扩展.三点弯曲梁模型跨中处钢纤维的分布如图3所示（
15、以钢纤维掺量0.8%的试件为例）.1.2.2钢纤维与水泥界面的本构关系在建立钢纤维增强水泥基复合材料细观有限元数值模型时，通常认为钢纤维和基体完全黏结，无相对滑动，因此将单根钢纤维与水泥砂浆基体拉拔试验所得拉拔力-位移曲线等效为钢纤维的拉伸应力-应变关系23-25，以此来模拟纤维与砂浆基体的黏结滑移作用.拉拔试验使用 P O 42.5普通硅酸盐水泥、细度模数为 2.6 的天然河砂和自来水来制作水泥砂浆基体，钢纤维为弹性模量 210 GPa的圆直型钢纤维.钢纤维埋深（Hf）为 20 mm，埋入角度为 0、30、45、60，每 种角度制备 10 个试件，试件尺寸为 100 mm100 mm100
16、mm.拉拔试验在 SUNS20kN 型万能试验机上进行，加载速度为 0.4 mm/min，当钢纤维被完全拔出或被拔断时结束试验.图 4 为倾斜纤维加载示意图.钢纤维的拉拔力-位移曲线如图 5所示.每种埋入角度取 3组离散性较小的数据作为试验结果.采用理 论 计 算 公 式 对 试 验 拉 拔 力-位 移 曲 线 进 行 计算26-27，结果如图5所示.由图5可见：（1）峰值荷载的理论计算结果与试验结果相差不大，位移存在一定误差.这是钢纤维倾斜角度增大导致的测量误差，总体上理论计算结果与试验结果吻合.（2）当埋入角度小于 45时，钢纤维拉拔力-位移曲线的峰值荷载随着埋入角度的增加逐渐增大且向右偏移；当埋入角度大于 45时，峰值荷载有减小的趋势.这是因为当纤维倾斜角度过大时，纤维埋入端处的基体应力集中程度加重，水泥砂浆破坏速度加快，使钢纤维被拔出所需要的拉拔力变小.为计算简便，在开展有限元模拟分析时，将理论计算的拉拔力-位移曲线等效为钢纤维的拉伸应力-应变关系，如图 6所示.当纤维倾斜角度为 0时，参考文献 18 进行等效计算；当倾斜角度不为 0时，引入几何校正系数（(Sf)）、等效应力（
