定向钢纤维增强水泥基复合材料断裂细观数值模拟_卿龙邦.pdf

第 26 卷第 2 期2023 年 2 月建筑材料学报JOURNAL OF BUILDING MATERIALSVol.26，No.2Feb.，2023定向钢纤维增强水泥基复合材料断裂细观数值模拟卿龙邦*，杨子美，慕儒，张金鑫，籍晓凯（河北工业大学 土木与交通学院，天津 300401）摘要：采用随机生成算法投放钢纤维，建立了随机乱向、定向钢纤维增强水泥基复合材料（SFRC、ASFRC）三点弯曲梁细观有限元数值模型，计算了不同纤维掺量下 SFRC试件和 ASFRC试件加载断裂的全过程，分析了三点弯曲梁开裂截面处的纤维应力，研究了定向钢纤维的细观增强机理.结果表明：SFRC试件和 ASFRC试件荷载-裂缝张开口位移全曲线的模拟值与试验值符合较好，峰值荷载的误差在 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.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 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 different 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 ASFRC 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 load 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 increase 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）；河北省自然科学基金资助项目（E2022202141，E2020202151）；河北省高等学校科学技术研究重点项目（ZD2019072，ZD2020190）第一作者（通讯作者）：卿龙邦（1982），男，湖北天门人，河北工业大学教授，博士生导师，博士.E-mail：建筑材料学报第 26卷在水泥基体中加入适量的钢纤维可以有效抑制裂缝的扩展，改善水泥基复合材料的断裂性能1-3.研究表明，钢纤维的增强作用受其分布方向的影响4-5.在传统随机乱向钢纤维增强水泥基复合材料（SFRC）中，纤维呈随机乱向分布，部分纤维分布方向与拉应力方向并不一致，难以充分发挥其增强作用，将钢纤维按特定方向分布则可以明显提高其增强效率6.采用磁场法可以实现钢纤维的定向分布.慕儒等7-9在盛有钢纤维增强水泥砂浆的试模外施加均匀电磁场，同时施加振动辅助，制备了定向钢纤维增强水泥基复合材料（ASFRC）并对其力学性能以及断裂性能进行了研究.Wijffels等10利用钢纤维在磁场中磁化的特性，制备了钢纤维定向分布的水泥基复合材料，研究了定向钢纤维对复合材料弯曲性能的影响.王哲伟等11通过对钢纤维增强地质聚合物施加均匀电磁场，制备了定向钢纤维增强地质聚合物预制缺口梁.汪洋等12采用磁场诱导定向技术，制备了定向碳纤维增强水泥砂浆，研究了碳纤维的取向性对碳纤维增强水泥砂浆力学性能提升效果的影响.细观数值模拟便于分析基体和钢纤维在断裂过程中任意时刻的受力状态，有利于深入研究钢纤维的增强机理.Wang等 13-14 通过分析钢纤维混凝土材料的应力-应变曲线，建立了钢纤维混凝土材料的本构关系，采用最大主应变破坏准则模拟了钢纤维混凝土在动态压缩荷载及爆破荷载下的破坏过程.Soetens等 15 假定纤维与基体无相对滑动，将纤维单元嵌入基体单元，建立了三维钢纤维混凝土有限元模型，模拟了端钩形钢纤维混凝土四点弯曲梁开裂的全过程.Fang等 16 建立了可预测SFRC动态破坏的三维有限元模型，研究了高应变率下SFRC的动态压缩性能.Yu等 17 在预计产生裂缝的截面处插入黏聚单元，将与裂纹张开速率相关的内聚定律作为黏聚单元的失效准则，模拟了端钩形钢纤维混凝土在不同落锤冲击速度下的断裂失效行为.本文采用随机生成算法投放钢纤维，建立了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试件中，钢纤维呈随机分布.本文参考文献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期卿龙邦，等：定向钢纤维增强水泥基复合材料断裂细观数值模拟在 ASFRC试件中，钢纤维分布方向与主拉应力方向一致，基于随机乱向钢纤维的投放算法，确定钢纤维掺量为 0.8%，试件尺寸为 440 mm100 mm100 mm，控制纤维投影与 X轴夹角为 0，建立定向钢纤维细观数值模型，如图2所示.1.2建立模型1.2.1基体本构关系基体选用三维八节点实体单元建模，钢纤维选用三维二节点桁架单元建模.采用黏聚裂纹模型19计算基体的开裂.在三点弯曲梁有限元模型跨中处插入零厚度的黏聚单元，将牵引分离定律作为黏聚单元的失效准则20-22，以黏聚单元的失效来模拟试件裂缝的产生和扩展.三点弯曲梁模型跨中处钢纤维的分布如图3所示（以钢纤维掺量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 mm.拉拔试验在 SUNS20kN 型万能试验机上进行，加载速度为 0.4 mm/min，当钢纤维被完全拔出或被拔断时结束试验.图 4 为倾斜纤维加载示意图.钢纤维的拉拔力-位移曲线如图 5所示.每种埋入角度取 3组离散性较小的数据作为试验结果.采用理 论 计 算 公 式 对 试 验 拉 拔 力-位 移 曲 线 进 行 计算26-27，结果如图5所示.由图5可见：（1）峰值荷载的理论计算结果与试验结果相差不大，位移存在一定误差.这是钢纤维倾斜角度增大导致的测量误差，总体上理论计算结果与试验结果吻合.（2）当埋入角度小于 45时，钢纤维拉拔力-位移曲线的峰值荷载随着埋入角度的增加逐渐增大且向右偏移；当埋入角度大于 45时，峰值荷载有减小的趋势.这是因为当纤维倾斜角度过大时，纤维埋入端处的基体应力集中程度加重，水泥砂浆破坏速度加快，使钢纤维被拔出所需要的拉拔力变小.为计算简便，在开展有限元模拟分析时，将理论计算的拉拔力-位移曲线等效为钢纤维的拉伸应力-应变关系，如图 6所示.当纤维倾斜角度为 0时，参考文献 18 进行等效计算；当倾斜角度不为 0时，引入几何校正系数（(Sf)）、等效应力（f）和等效应变（f），按式（2）（4）进行计算.图 2定向钢纤维细观数值模型Fig.2Meso-numerical model of aligned steel fibers图 3三点弯曲梁模型跨中处钢纤维的分布Fig.3Distribution of steel fibers in the middle of the three point bending beam model图 4倾斜纤维加载示意图Fig.4Schematic diagram of inclined steel fiber loading（size：mm）113建筑材料学报第 26卷f=4FfD2fcos(Sf)（2）f=SfHf1cos()Sf（3）(Sf)=arctan(Hfsin 0-Sfcot 0Hfcos 0+Sf)（4）式中：Ff为钢纤维拉拔荷载，kN；Sf为钢纤维滑移位 移，mm；0为 钢 纤 维 与 拉 拔 力 之 间 的 初 始 夹角，（）.根据纤维投影与 X 轴的夹角范围，分别指定相应角度的等效应力-应变关系.基于扩展有限元法，采 用 位 移 加 载 的 方 式 模 拟 荷 载 作 用 下 SFRC 和ASFRC三点弯曲梁断裂破坏的全过程.2断裂试验采用 P O 42.5普通硅酸盐水泥、细度模数为 2.6的天然河砂、钢纤维、自来水和聚羧酸型高效减水剂来制作钢纤维增强水泥砂浆.调试后确定基体水灰比（质量比）为 0.36，设计纤维掺量为 0.8%、1.2%、2.0%，具体的材料配合比见表 1.钢纤维的特征参数见表2.其中：E为弹性模量，ft为抗拉强度.图 5钢纤维的拉拔力-位移曲线Fig.5Curves of pull-out force-displacement of steel fiber图 6基于拉拔力-位移曲线确定的钢纤维等效应力-应变关系曲线Fig.6Determination of steel fiber stress-strain curves base on pull-out force-displacement curves114第 2期卿龙邦，等：定向钢纤维增强水泥基复合材料断裂细观数值模拟每种钢纤维掺量下均制备SFRC试件和ASFRC试件，试件尺寸为 440 mm100 mm100 mm，跨高比为 4，跨中设置宽 2 mm的预制裂缝，缝高比为 0.4，标准养护 28 d.制备 SFRC 试件时采用常规方式振捣，制备 ASFRC试件时需将盛有钢纤维增强水泥砂浆的试模置于振动台上的通电线圈内进行振捣，如图7所示.断裂试验在伺服万能试验机上进行，采用位移加载方式，加载速率为 0.15 mm/min.选用量程为20 kN 的外接荷载传感器测量荷载（P），在试件裂缝跨中底部设置 YYJ-10/10 型夹式引伸计，用于测量裂缝张开位移（CMOD）.试验加载示意图如图 8 所示.水泥砂浆基体的材料力学参数如表 3所示.其中：为泊松比，GF为断裂能.3结果及分析3.1P-CMOD曲线及峰值荷载SFRC试件和 ASFRC试件的荷载-裂缝张开口位移（P-CMOD）曲线如图9所示.由图9可见：（1）当 钢 纤 维 掺 量 为 0.8%时，SFRC 试 件 的P-CMOD 曲线有明显的陡降，随着钢纤维掺量的增加，陡降程度有所缓解，而ASFRC试件的P-CMOD曲线在达到峰值荷载后均下降缓慢.这是因为钢纤维定向后，其分布方向与试件所受主拉应力方向一致，桥接裂缝的钢纤维得到了充分利用，阻裂增韧效率提高.（2）当钢纤维掺量为 0.8%和 1.2%时，SFRC 试件模拟曲线的下降段略低于试验曲线，但曲线趋势大致相同，当钢纤维掺量为 2.0%时，曲线吻合较好；ASFRC试件的模拟曲线与试验曲线均吻合较好.因此，模拟结果可以较好地反映2种试件断裂的全过程，表明建立的细观有限元数值模型具有一定的有效性.试验与模拟所得 SFRC 试件和 ASFRC 试件的峰值荷载如表 4所示.每组 3块试件，取平均值作为该组的试验峰值荷载.由表 4 可见：SFRC 试件和ASFRC 试件的试验峰值荷载均随着钢纤维掺量的增加而增大，且 ASFRC试件的峰值荷载的增幅大于SFRC 试件；当钢纤维掺量为 0.8%、1.2%、2.0%时，ASFRC 试件的峰值荷载较 SFRC 试件约提高了75%、111%、141%，表明钢纤维定向可以有效改善水泥基复合材料的断裂性能，提高钢纤维的利用率，充分发挥增强作用；3种钢纤维掺量下，SFRC试件峰值荷载模拟值与试验值的误差分别为6.30%、6.45%、4.07%，ASFRC 试件的误差分别为 4.87%、9.01%、表 1砂浆基体的配合比Table 1Mix proportions of mortar matrixVf/%0.81.22.0Mix proportion/(kg m-3)Water235.00235.00235.00Cement651.00651.00651.00Sand1 302.001 302.001 302.00Superplasticizer2.562.643.28表 2钢纤维的特征参数Table 2Characteristic parameters of steel fibersSteel fiber typeRound straightDensity/(kg m-3)7 850Lf/mm30Df/mm0.5Lf/Df60E/GPa210ft/MPa1 150图 7ASFRC试件的制备装置Fig.7Preparation device of ASFRC specimen8图 8断裂试验加载示意图Fig.8Loading schematic diagram of fracture test（size：mm）表 3水泥砂浆基体的材料力学参数Table 3Material mechanical parameters of cement mortar matrixft/MPa3.2E/GPa300.2GF/(N mm-1)0.148115建筑材料学报第 26卷图 9SFRC试件和 ASFRC试件的荷载-裂缝张开口位移曲线Fig.9P-CMOD curves of SFRC and ASFRC specimens表 4试验与模拟所得 SFRC试件和 ASFRC试件的峰值荷载Table 4Tested and simulated peak loads of SFRC and ASFRC specimensSpecimen No.SFRC-0.8-1SFRC-0.8-2SFRC-0.8-3ASFRC-0.8-1ASFRC-0.8-2ASFRC-0.8-3SFRC-1.2-1SFRC-1.2-2SFRC-1.2-3ASFRC-1.2-1ASFRC-1.2-2ASFRC-1.2-3SFRC-2.0-1SFRC-2.0-2SFRC-2.0-3ASFRC-2.0-1ASFRC-2.0-2ASFRC-2.0-3Specimen typeSFRCASFRCSFRCASFRCSFRCASFRCVf/%0.80.81.21.22.02.0Tested peak load/kN3.092.562.444.404.725.042.762.842.766.485.655.502.773.884.438.428.519.76Mean value/kN2.704.722.795.883.698.90Simulated peak load/kN2.534.492.975.353.848.97116第 2期卿龙邦，等：定向钢纤维增强水泥基复合材料断裂细观数值模拟0.79%；2种试件的误差均在10%以内，误差较小.3.2跨中开裂截面处纤维的应力分析图 10为试件开裂破坏图及跨中开裂截面处的钢纤维分布.由图 10可见，SFRC试件中的钢纤维呈乱向分布，难以充分发挥桥接裂缝的作用，而 ASFRC试件中的钢纤维全部与开裂面垂直.不同加载时刻 ASFRC 试件主拉应力方向（x）的应力云图如图 11所示（以钢纤维掺量 1.2%的试件为例）.由图 11可见：在加载初期，裂缝尚未产生，基体和钢纤维共同发挥承载作用（图 11（a）；当 P=3.52 kN时，基体已产生裂缝（图 11（b）；随着荷载的继续增加，裂缝向上扩展，基体逐渐退出工作，钢纤维与基体之间的黏结滑移发挥主要承载作用，消耗大部分能量，增强基体开裂后的延性（图 11（c）；裂缝处钢纤维的应力减小，表明钢纤维开始滑动拔出，继续加载，试件逐渐丧失承载能力，沿着裂缝破坏（图11（d）、（e）.图 12 为试件跨中开裂截面处的纤维应力分布图（以钢纤维掺量 1.2%的试件为例）.以 ASFRC 试件为例进行分析：在加载初期，纤维合力（F）较小，开裂截面处上部钢纤维受压，下部钢纤维受拉（见图12（b）；随着荷载的增加，受拉区域扩大，受拉钢纤维的数量增加，钢纤维合力逐渐增大（见图 12（d）、（f）；当钢纤维合力增至最大值时，只有试件顶部少量钢纤维受压（见图 12（h）；继续加载，钢纤维合力减 小，试 件 逐 渐 丧 失 承 载 能 力，沿 着 裂 缝 破 坏；ASFRC 试件跨中开裂截面上的钢纤维数量明显大于SFRC试件，SFRC试件和ASFRC试件的钢纤维合力最大值分别为 7.082 kN和 17.283 kN（见图 12（g）、（h），ASFRC试件的钢纤维合力最大值较 SFRC试件约提高了 144%，表明钢纤维定向后的利用率有大幅度提升.图 13为试件跨中开裂截面处的钢纤维合力-荷载（F-P）曲线图.由图13可见：（1）当 钢 纤 维 掺 量 为 0.8%、1.2%、2.0%时，SFRC 试 件 的 钢 纤 维 合 力 最 大 值 分 别 为 3.985、7.082、8.970 kN，ASFRC 试件的钢纤维合力最大值分别为 12.047、17.283、20.361 kN，ASFRC 试件钢纤 维 合 力 的 最 大 值 较 SFRC 试 件 分 别 增 大 了 约202%、144%、127%，且 ASFRC 试件的增幅随着钢纤维掺量的增加而减小.（2）2 种试件的钢纤维合力均先增大后减小，且当钢纤维掺量相同时，在整个加载过程中 ASFRC试件的钢纤维合力始终大于 SFRC试件.在加载初期，钢纤维合力近似为线性增长，钢纤维处于弹性阶段，此时的伸长量为自身弹性变形；随着荷载的增加，裂缝开始向上扩展，钢纤维合力呈非线性增长，当达到峰值荷载时，开裂面上钢纤维的端部开始与基体脱黏，由钢纤维侧面与基体的黏结力承载主要荷载，基体开裂后的延性增强；荷载减小，钢纤维合力继续增大，当钢纤维合力增至最大值时，钢纤维与基体完全脱黏；钢纤维合力减小表明纤维被拔出，此时由钢纤维未拔出部分与基体之间的摩擦力共同抵抗荷载作用，继续加载，试件逐渐丧失承载能力，沿裂缝破坏.当钢纤维掺量为 2.0%时，F-P曲线上的 Pmax点和 Fmax点较钢纤维掺量为 0.8%和 1.2%时更接近，这是由于随着钢纤维纤维掺量的增加，峰值荷载逐渐增大，同时钢纤维掺量的增加使钢纤维间距减小，钢纤维与基体之间的黏结作用减小，更易滑动被拔出，钢纤维合力增大的幅度减小，由此Pmax点和Fmax点较为接近.图 10试件开裂破坏图及跨中开裂截面处的钢纤维分布Fig.10Cracking failure diagram of the specimen and distribution of steel fibers at the cracked section in the middle117建筑材料学报第 26卷图 11不同加载时刻 ASFRC试件的主拉应力方向应力云图Fig.11x stress cloud diagram of ASFRC specimen at different loading moments118第 2期卿龙邦，等：定向钢纤维增强水泥基复合材料断裂细观数值模拟图 12试件跨中开裂截面处纤维的应力分布图Fig.12Fiber stress distribution diagram on cracked surface in mid-span of specimens图 13试件跨中开裂截面处的纤维合力-荷载曲线图Fig.13F-P curves of fibers on cracked surface in mid-span of specimens119建筑材料学报第 26卷4结论（1）钢纤维拉拔力的峰值随着钢纤维倾角的增加呈现先增大后减小的趋势.（2）随机乱向、定向钢纤维增强水泥基复合材料（SFRC、ASFRC）试件的峰值荷载均随着钢纤维掺量的增加而增大，相同钢纤维掺量下 ASFRC试件的峰值荷载明显大于 SFRC 试件.当钢纤维掺量为0.8%、1.2%、2.0%时，ASFRC 试件的峰值荷载较SFRC 试件分别提高了 75%、111%、141%，表明钢纤维定向可以有效改善钢纤维增强水泥基复合材料的断裂性能.（3）SFRC 和 ASFRC 三点弯曲梁断裂全过程荷载-裂缝张开口位移全曲线的计算值与试验值符合较好，峰值荷载误差在 10%以内，验证了有限元模型的有效性.（4）当钢纤维掺量相同时，ASFRC 试件跨中开裂面上的钢纤维数量大于SFRC试件；当钢纤维掺量为 0.8%、1.2%、2.0%时，ASFRC试件纤维合力的最大值较 SFRC分别提高了 202%、144%、127%.表明钢纤维定向后可以显著提高钢纤维的利用率，延缓裂缝的扩展.参考文献：1 徐平，郑满奎，王超，等.考虑尺寸及纤维掺量影响的高强混凝土断裂能试验研究 J.硅酸盐通报，2020，39（11）：3488-3495.XU Ping，ZHENG Mankui，WANG Chao，et al.Experimental study on fracture energy of high strength concrete considering the influence of size and fiber content J.Bulletin of the Chinese Ceramic Society，2020，39（11）：3488-3495.（in Chinese）2 ZHANG X X，ELAZIM A，RUIZ G，et al.Fracture behaviour of steel fibre-reinforced concrete at a wide range of loading ratesJ.International Journal of Impact Engineering，2014，71：89-96.3 GAO D Y，ZHANG T Y.Fracture characteristics of steel fiber reinforced high strength concrete under three-point bending J.Journal of the Chinese Ceramic Society，2007，35（12）：1630-1635.4 FU S Y，LAUKE B.Effects of fiber length and fiber orientation distributions on the tensile strength of short-fiber-reinforced polymers J.Composites Science and Technology，1996，56（10）：1179-1190.5 BRANDT A M.On the optimal direction of short metal fibres in brittle matrix composites J.Journal of Materials Science，1985，20（11）：3831-3841.6 ABRISHAMBAF A，BARROS J A，CUNHA V M.Relation between fibre distribution and post-cracking behaviour in steel fibre reinforced self-compacting concrete panelsJ.Cement Concrete Research，2013，51：51-66.7 慕儒，邱欣，赵全明，等.单向分布钢纤维增强水泥基复合材料（I）：钢纤维方向控制 J.建筑材料学报，2015，18（2）：208-213.MU Ru，QIU Xin，ZHAO Quanming，et al.Aligned steel fiber reinforced cement based composites（I）：Orientation control of steel fibersJ.Journal of Building Materials，2015，18（2）：208-213.（in Chinese）8 MU R，LI H，QING L B，et al.Aligning steel fibers in cement mortar using electro-magnetic field J.Construction and Building Materials，2017，131：309-316.9 慕儒，王成，李辉，等.单向分布钢纤维增强水泥基复合材料（III）：断裂性能 J.建筑材料学报，2016，19（1）：78-82.MU Ru，WANG Cheng，LI Hui，et al.Aligned steel fiber reinforced cement based composites（III）：Fracture properties J.Journal of Building Materials，2016，19（1）：78-82.（in Chinese）10 WIJFFELS M，WOLFS R，SUIKER A，et al.Magnetic orientation of steel fibres in self-compacting concrete beams：Effect on failure behavior J.Cement and Concrete Composites，2017，80：342-355.11 王哲伟，周华飞，谢子令.定向钢纤维增强地质聚合物弯曲破坏的声发射特性 J.材料导报，2021，35（12）：12214-12219.WANG Zhewei，ZHOU Huafei，XIE Ziling.Acoustic emission characteristics of aligned steel fiber reinforced geopolymer subjected to bending failure J.Materials Reports，2021，35（12）：12214-12219.（in Chinese）12 汪洋，徐金霞，蒋林华，等.磁场诱导定向碳纤维增强水泥砂浆的力学性能 J.复合材料学报，2019，36（11）：2726-2733.WANG Yang，XU Jinxia，JIANG Linhua，et al.Mechanical properties of cement mortar reinforced with aligned fibers by magnetic field J.Acta Materiae Compositae Sinica，2019，36（11）：2726-2733.（in Chinese）13 WANG Z L，WU L P，WANG J G.A study of constitutive relation and dynamic failure for SFRC in compressionJ.Construction and Building Materials，2010，24（8）：1359-1363.14 WANG Z L，KONIETZKY H，HUANG R Y.Elastic-plastic-hydrodynamic analysis of crater blasting in steel fiber reinforced concrete J.Theoretical and Applied Fracture Mechanics，2009，52（2）：111-116.15 SOETENS T，MATTHYS S.Different methods to model the post-cracking behaviour of hooked-end steel fibre reinforced concrete J.Construction and Building Materials，2014，73（30）：458-471.16 FANG Q，ZHANG J.Three-dimensional modelling of steel fiber reinforced concrete material under intense dynamic loading J Construction and Building Materials，2013，44：118-132.17 YU R C，CIFUENTES H，RIVERO I，et al.Dynamic fracture behavior in fiber-reinforced cementitious composites J.Journal of the Mechanics and Physics of Solids，2016，93：135-152.18 喻渴来，卿龙邦，王苗.定向钢纤维增强水泥基复合材料断裂特性模拟分析 J.硅酸盐通报，2018，37（3）：810-817.YU Kelai，QING Longbang，WANG Miao.Simulation analysis of fracture characteristics for aligned steel fiber reinforced cementitious compositesJ.Bulletin of the Chinese Ceramic 120第 2期卿龙邦，等：定向钢纤维增强水泥基复合材料断裂细观数值模拟Society，2018，37（3）：810-817.（in Chinese）19 SPRING D W，PAULINO G H.A growing library of three-dimensional cohesive elements for use in ABAQUS J.Engineering Fracture Mechanics，2014，126：190-216.20 PARK K，CHOI H，PAULINO G H.Assessment of cohesive traction-separation relationships in ABAQUS：A comparative study J.Mechanics Research Communications，2016，78：71-78.21 POBLETE F R，MONDAL K，MA Y，et al.Direct measurement of rate-dependent mode I and mode II traction-separation laws for cohesive zone modeling of laminated glass J.Composite Structures，2022，279：47-59.22 PARK K，PAULION G H.Cohesive zone models：A critical review of traction-separation relationships across fracture surf