摘要
考虑纤维在自密实混凝土中的分布,提出一种用于模拟钢纤维自密实混凝土(SFRSCC)抗折性能的数值方法.首先,采用基于流体动力学的方法模拟新拌SFRSCC的浇筑过程,以确定钢纤维的分布;然后,基于钢纤维在SFRSCC试件中的分布,考虑钢纤维取向角度对钢纤维拔出响应的影响,建立由混凝土基体和钢纤维组成的SFRSCC细观有限元模型;最后采用提出的数值方法对SFRSCC梁试件的浇筑过程和抗折试验进行细观模拟.模拟结果与试验结果吻合良好,说明提出的数值方法能够模拟SFRSCC的抗折性能.
关键词
钢纤维混凝土(SFRC)力学性能的增强程度与钢纤维拔出响应有着重要的联系.为考虑钢纤维拔出响应对SFRC的增强作用,研究人员提出了一系列模拟SFRC力学性能的数值模
为此,本文基于钢纤维在自密实混凝土(SCC)中的分布,提出了一种模拟SFRSCC抗折性能的数值方法.首先,采用基于流体动力学(CFD)的数值方法模拟SFRSCC的浇筑过程并确定钢纤维在混凝土基体中的分布;然后,基于钢纤维分布,同时考虑钢纤维的拔出响应,建立模拟SFRSCC抗折性能的细观有限元模型;最后,基于试验,开展了对SFRSCC浇筑过程和抗折模拟的研究.
数值方法包括用于确定钢纤维分布的CFD数值方法和用于模拟SFRSCC抗折性能的有限元方法.
采用文献[
CFD数值方
(1)根据试件形状和新拌混凝土的材料性能,利用软件CFX建立CFD模型.
(2)对CFD模型进行数值求解,得到新拌混凝土的速度场.
(3)利用编程手段在计算域浇筑入口上方生成随机分布的钢纤维,假设其初始速度与混凝土浇筑速度相等,如

图1 钢纤维与计算域的初始状态
Fig.1 Initial state of steel fibers and computational domain
(4)根据新拌混凝土的速度场,计算各时间步中钢纤维的平动和转动.根据所有钢纤维的坐标,确定钢纤维的取向和分布.
完成新拌SFRSCC的浇筑计算后,确定各计算步中的钢纤维取向和分布.将新拌混凝土中钢纤维的最终分布作为有限元模型的输入参数,进而模拟SFRSCC的抗折性能.
采用ABAQUS软件中混凝土损伤塑性(CDP)模型描述混凝土基体的非线性行为.为定义CDP模型,需在模型中输入混凝土单轴拉伸和压缩的力学响应.采用GB 50010—2010《混凝土结构设计规范》中的公式描述单轴受压状态下混凝土的应力-应变关系.通过定义断裂能的形式来表征混凝土的拉伸开裂行为,定义时只需要输入混凝土的断裂能和抗拉强度.然后,在模型中输入其他描述CDP模型的相关参数:膨胀角为31°,偏心率为0.1,双轴受压与单轴受压强度比值为1.16,拉压子午线第二应力不变量之比为0.667,黏性系数为0.000 1.
在有限元模型中,钢纤维基于新拌SFRSCC浇筑过程中最终钢纤维的分布而生成.为考虑钢纤维对混凝土基体的增强作用,参考Yu
将钢纤维视为一种新的材料,其本构关系基于钢纤维的拔出荷载-滑移曲线计算而得.钢纤维的应力和应变按下式计
(1) |
(2) |
式中:为钢纤维直径;为钢纤维的长度.
钢纤维的取向角度对钢纤维的拔出响应有着重要影响,不同取向角度下钢纤维的拔出响应具有明显的差
平直钢纤维拔出行为的典型特征主要包括纤维脱黏、基体剥落、摩擦滑动和纤维拔出.相比之下,端钩型钢纤维除了包含平直钢纤维的所有拔出特征外,基体的剥落程度更显著,端钩的存在会使得钢纤维发生塑性变形.基于端钩型钢纤维从基体中拔出过程的特征,Laranjeira

图2 端钩型钢纤维拔出响应的确定
Fig.2 Determination of pull‑out response of hooked end fibers
为了验证提出的数值方法,采用2种钢纤维体积分数的SFRSCC浇筑和制作了梁试件,采用四点弯曲试验评价梁试件的抗折性能.采用图像处理技术评价了钢纤维的取向和分布.
采用中心浇筑方式制备了6个尺寸为100 mm×100 mm×400 mm的SFRSCC梁试件,其中SCC配合比为:水泥428 kg/
Specimen | Slump flow diameter/mm | T50/s | V‑funnel time/s | Plastic viscosity/(Pa·s) | Yield stress/Pa |
---|---|---|---|---|---|
SFRSCC‑0.6 | 680 | 6.1 | 12.6 | 30.7 | 34.1 |
SFRSCC‑0.9 | 650 | 6.4 | 14.5 | 42.4 | 38.7 |
在完成四点弯曲试验后,采用图像处理技术确定SFRSCC试件截面上的钢纤维分布和取向,试件切割位置如

图3 梁试件的切割位置
Fig.3 Cutting positions of cross section of beam specimens(size: mm)
首先,基于图像处理技术,确定试件截面上的钢纤维根数及钢纤维与梁截面之间的夹角(取向角度).然后,采用2个钢纤维系数来评价截面上钢纤维取向和分布:一是钢纤维密度,即截面上的钢纤维根数除以截面面积;二是钢纤维取向系数,其计算式为,的值越趋近于1,表明截面上的钢纤维越趋近于与截面法向平行分布,此时钢纤维对抗折性能的提升效果越显著.
利用CFD软件CFX建立梁试件的数值模型(

图4 梁试件的数值模型
Fig.4 Numerical model of beam specimens(size: mm)
试件SFRSCC‑0.6和SFRSCC‑0.9对应的钢纤维根数分别为1 494、2 241根.根据软件CFX计算得到的新拌混凝土速度场,计算了各时间步钢纤维的运动,进而得到了各时间步的钢纤维取向和分布结果.详细的新拌混凝土浇筑过程模拟和钢纤维运动计算过程可参考文献[
新拌SFRSCC‑0.6和SFRSCC‑0.9的流动状态与钢纤维分布如

图5 新拌SFRSCC‑0.6和SFRSCC‑0.9的流动状态与钢纤维分布
Fig.5 Transient flow of fresh SFRSCC‑0.6, SFRSCC‑0.9 and fiber distribution

图6 5.0 s时距梁底部0.02 m水平面上新拌混凝土的速度云图
Fig.6 Concrete velocity contour on horizontal plane 0.02 m away from bottom of the beam at 5.0 s
Specimen | Closer to the end of beam | Closer to the crack plane | On the longitudinal section | ||||||
---|---|---|---|---|---|---|---|---|---|
Experimental | Simulated | Error/% | Experimental | Simulated | Error/% | Experimental | Simulated | Error/% | |
SFRSCC‑0.6 | 0.707 | 0.670 | 5.2 | 1.053 | 1.050 | 0.3 | 0.653 | 0.690 | 5.7 |
SFRSCC‑0.9 | 1.057 | 0.970 | 8.2 | 1.510 | 1.593 | 5.5 | 0.977 | 1.010 | 3.4 |
Specimen | Closer to the end of beam | Closer to the crack plane | On the longitudinal section | ||||||
---|---|---|---|---|---|---|---|---|---|
Experimental | Simulated | Error/% | Experimental | Simulated | Error/% | Experimental | Simulated | Error/% | |
SFRSCC‑0.6 | 0.690 | 0.686 | 0.6 | 0.777 | 0.813 | 4.6 | 0.727 | 0.737 | 5.9 |
SFRSCC‑0.9 | 0.681 | 0.656 | 3.7 | 0.753 | 0.784 | 4.1 | 0.704 | 0.723 | 2.7 |

图7 SFRSCC梁试件的有限单元网格
Fig.7 Finite element mesh of SFRSCC beam specimens
Poisson ratio | Elastic modulus/GPa | Compressive strength/MPa | Tensile strength/MPa | Fracture energy/(N·m |
---|---|---|---|---|
0.2 | 37.5 | 31.7 | 1.6 | 0.114 |
采用桁架单元T3D2对钢纤维进行建模,根据新拌SFRSCC浇筑数值模拟得到的钢纤维最终位置来确定钢纤维在混凝土基体中的分布.将钢纤维取向角度设置为12.5°、37.5°和62.5°,其拔出响应曲线分配给取向角度范围为[0°,25°]、[25°,50°]和[50°,75°]的桁架单元,并假设取向角度处于[75°,90°]的钢纤维对混凝土基体不起增强作

图8 3种取向角度下端钩型钢纤维的拔出荷载-滑移曲线
Fig.8 P‑s curves of hooked end fibers with three fiber orientation angles

图9 荷载-挠度曲线的试验结果和模拟结果
Fig.9 Experimental and simulated results of load‑deflection curves

图10 混凝土基体和钢纤维在峰值荷载时的竖向位移云图
Fig.10 Vertical displacement contours of concrete matrix and steel fibers at peak load
在不同钢纤维体积分数下,对比模拟和试验得到的峰值荷载Fp、峰值荷载对应的挠度p以及荷载-挠度曲线与坐标轴围成的面积(韧性K),结果见
Specimen | Fp/kN | Error of Fp /% | /mm | Error of /% | K/(kN·mm) | Error of K/% | |||
---|---|---|---|---|---|---|---|---|---|
Experimental | Simulated | Experimental | Simulated | Experimental | Simulated | ||||
SFRSCC‑0.6 | 16.63 | 16.45 | 1.1 | 0.173 | 0.181 | 4.6 | 20.99 | 20.60 | 1.9 |
SFRSCC‑0.9 | 19.27 | 18.74 | 2.8 | 0.290 | 0.340 | 17.2 | 25.42 | 24.81 | 2.4 |
(1)采用流体动力学(CFD)数值方法模拟新拌钢纤维自密实混凝土SFRSCC浇筑过程中钢纤维的运动.基于流动模拟得到的钢纤维分布,考虑钢纤维取向对钢纤维拔出响应的影响和钢纤维对混凝土基体的增强作用,建立了用于模拟SFRSCC抗折性能的细观有限元模型.将CFD数值方法和有限元方法结合,建立了能够考虑钢纤维分布的SFRSCC抗折性能模拟方法.
(2)钢纤维在梁试件截面上的取向和分布不是完全随机均匀的,而是显著受SFRSCC流变性和壁面效应影响的.越靠近模板壁面,钢纤维越倾向于与模板壁面平行分布,且模板端部附近的钢纤维呈现出自模板底面倾斜向上分布的趋势.
(3)钢纤维密度、钢纤维取向系数和荷载-挠度曲线的模拟结果与试验结果吻合较好,说明提出的钢纤维自密实混凝土数值方法能较为准确地模拟SFRSCC浇筑过程中钢纤维的运动和SFRSCC的抗折性能.
参考文献
BI J H, HUO L Y, ZHAO Y, et al. Modified the smeared crack constitutive model of fiber reinforced concrete under uniaxial loading[J]. Construction and Building Materials, 2020, 250:985‑995. [百度学术]
LUCCIONI B, RUANO G, ISLA F, et al. A simple approach to model SFRC[J]. Construction and Building Materials, 2012, 37:111‑124. [百度学术]
DENNEMAN E, WU R Z, KEARSLEY E P, et al. Discrete fracture in high performance fibre reinforced concrete materials[J]. Engineering Fracture Mechanics, 2011, 78 (10):2235‑2245. [百度学术]
RADTKE F K F, SIMONE A, SLUYS L J. A computational model for failure analysis of fibre reinforced concrete with discrete treatment of fibres[J]. Engineering Fracture Mechanics, 2010, 77 (4):597‑620. [百度学术]
BITENCOURTJR L A G, MANZOLI O L, BITTENCOURT T N, et al. Numerical modeling of steel fiber reinforced concrete with a discrete and explicit representation of steel fibers[J]. International Journal of Solids and Structures, 2019, 159:171‑190. [百度学术]
YU R C, CIFUENTES H, RIVERO I, et al. Dynamic fracture behaviour in fibre‑reinforced cementitious composites[J]. Journal of the Mechanics and Physics of Solids, 2016, 93:135‑152. [百度学术]
卿龙邦,喻渴来,慕儒,等. 定向钢纤维增强水泥基复合材料抗拉细观模拟[J]. 建筑材料学报, 2018, 21(4):561‑567. [百度学术]
QING Longbang, YU Kelai, MU Ru, et al. Meso‑scale simulation of tension properties for aligned steel fiber reinforced cementitious composites[J]. Journal of Building Materials, 2018, 21(4):561‑567. (in Chinese) [百度学术]
赵健,廖霖,张帆,等.钢纤维混凝土弯曲性能和纤维分布试验研究[J]. 建筑材料学报, 2020, 23(4):838‑845. [百度学术]
ZHAO Jian, LIAO Lin, ZHANG Fan, et al. Experimental study on flexural properties and fiber distribution of steel fiber reinforced concrete[J]. Journal of Building Materials, 2020, 23(4):838‑845. (in Chinese) [百度学术]
MARTINIE L, ROUSSEL N. Simple tools for fiber orientation prediction in industrial practice[J]. Cement and Concrete Research, 2011, 41(10):993‑1000. [百度学术]
RAJU R A, LIM S, AKIYAMA M, et al. Effects of concrete flow on the distribution and orientation of fibers and flexural behavior of steel fiber‑reinforced self‑compacting concrete beams[J]. Construction and Building Materials, 2020, 262:119963. [百度学术]
CUNHA V M C F, BARROS J A O, SENA‑CRUZ J M. An integrated approach for modelling the tensile behaviour of steel fibre reinforced self‑compacting concrete[J]. Cement and Concrete Research, 2011, 41(1):64‑76. [百度学术]
MECHTCHERINE V, GRAM A, KRENZER K, et al. Simulation of fresh concrete flow using discrete element method(DEM):Theory and applications[J]. Materials and Structures, 2014, 47(4):615‑630. [百度学术]
毕继红,赵云,鲍春,等. 不同体积掺量的钢纤维在混凝土中的分布研究[J]. 湖南大学学报(自然科学版), 2021, 48(3):99‑108. [百度学术]
BI Jihong, ZHAO Yun, BAO Chun, et al. Research on distribution of steel fiber with different volume content in concrete[J]. Journal of Hunan University(Natural Sciences), 2021, 48(3):99‑108. [百度学术]
LARANJEIRA F, MOLINS C, AGUADO A. Predicting the pullout response of inclined hooked steel fibers[J]. Cement and Concrete Research, 2010, 40 (10):1471‑1487. [百度学术]
CAO Y Y Y, YU Q L. Effect of inclination angle on hooked end steel fiber pullout behavior in ultra‑high performance concrete[J]. Composite Structures, 2018, 201:151‑160. [百度学术]
ISLA F, RUANO G, LUCCIONI B. Analysis of steel fibers pull‑out. Experimental study[J]. Construction and Building Materials, 2015, 100:183‑193. [百度学术]
YOO D Y, JE J, CHOI H J, et al. Influence of embedment length on the pullout behavior of steel fibers from ultra‑high‑performance concrete[J]. Materials Letters, 2020, 276:128233. [百度学术]
ZHANG Y, ZHU Y P, QU S Q, et al. Improvement of flexural and tensile strength of layered‑casting UHPC with aligned steel fibers[J]. Construction and Building Materials, 2020, 251:118893. [百度学术]
ZHOU B, UCHIDA Y. Influence of flowability, casting time and formwork geometry on fiber orientation and mechanical properties of UHPFRC[J]. Cement and Concrete Research, 2017, 95:164‑177. [百度学术]
ZHANG H, HUANG Y J, YANG Z J, et al. A discrete‑continuum coupled finite element modelling approach for fibre reinforced concrete[J]. Cement and Concrete Research, 2018, 106:130‑143. [百度学术]
DENG F Q, DING X X, CHI Y, et al. The pull‑out behavior of straight and hooked‑end steel fiber from hybrid fiber reinforced cementitious composite:Experimental study and analytical modelling[J]. Composite Structure, 2018, 206:693‑712. [百度学术]