3D mesoscale modeling method and dynamic mechanical properties investigation of concrete
-
摘要: 根据混凝土材料的细观组成和结构特点,基于三维Voronoi图形提出了一种简单高效的混凝土细观模型生成方法,利用塑性损伤模型对该细观模型进行了单、多轴应力状态下的准静态分析以及SHPB动态有限元分析。结果表明,数值模拟得到的应力应变曲线和破坏模式与实验结果基本吻合,本文中提出的混凝土三维细观模型可较好地模拟混凝土的静、动态力学特性,为进一步从细观力学角度研究混凝土损伤演化规律和破坏机理提供了模型基础。Abstract: According to the mechanical properties and components of concrete from the mesoscale viewpoint, a simple and high-efficiency method based on 3D-Voronoi is proposed to construct realistic meso-structure model of concrete. In this paper, 3D mesoscale numerical simulation is conducted using damaged plasticity modeling under quasi-static and dynamic load. The simulation results resemble the corresponding experimental results both in strain-stress curves and failure modes. The current method can simulate static and dynamic mechanical properties well. In addition, the proved 3D model can provide an effective approach in elucidating fundamental mechanical behavior of concrete-like materials.
-
混凝土材料从细观层次上由骨料、砂浆以及两者之间的黏接界面层(interfacial transitional zone,ITZ)组成,各组分的细观力学性质直接或间接地决定了混凝土材料的宏观力学性能和破坏模式[1-2]。根据混凝土细观结构特点及其力学性能,深入研究细观结构和宏观力学性能之间的关系,是混凝土力学性质的重要研究内容之一。
Wittmann等[3]首先将多尺度的研究方法应用到混凝土材料的研究中,认为混凝土是一种非均质复合材料,提出了用随机分布的多边形来模拟混凝土骨料的细观力学模型。该方法一直被后来的研究者所采用。Zhou等[4-5]建立了圆形随机骨料模型,利用材料损伤本构关系,系统研究了混凝土细观模型在不同受力条件下的拉伸、压缩和爆炸问题,指出骨料和砂浆之间的黏接界面层在压缩时先出现破坏裂纹,在拉伸时对混凝土失效机理和拉伸强度影响最大。Wang等[6]在椭圆形和多边形骨料的基础上引入孔隙,并研究了骨料形状、骨料体积分数和孔隙率对抗拉强度的影响。吴成等[7]对刚性弹丸侵彻细观混凝土靶进行了数值模拟,分析了砂浆种类、粗骨料种类和粗骨料体积分数等对靶板抗侵彻能力的影响,并通过扩展Forrestal阻力方程,建立了细观混凝土侵彻深度模型。相比于二维模型,三维细观模型更接近于实际混凝土的真实形态,因此能较好地反映混凝土材料的实际变形与损伤破坏情况。Wang等[8]提出了“投放算法”用于产生随机分布的三维骨料模型,并在单元之间嵌入零厚度黏接单元,以此来模拟混凝土在拉应力作用下的开裂行为。Zhang等[9]通过对生成的骨料进行平移和旋转来控制骨料的移动,提出了“随机爬行算法”,该算法可以实现更高骨料含量的混凝土建模。邓勇军等[10]基于骨料随机投放的思想建立了混凝土三维细观几何模型,分析了刚性弹正侵彻过程中发生弹道偏转的原因及可能的影响因素,定量研究了混凝土细观因素对弹道偏转的影响。然而,骨料投放往往需要复杂的算法对骨料相互之间是否发生侵入进行判断,同时需要花费较多的时间,建模效率较低。三维Voronoi图形中胞元形状与混凝土骨料具有良好的相似性,其胞元可实现随机的空间形状和空间分布,可作为建立骨料的一种方法。
本文中基于三维Voronoi图形,提出一种简单、高效的混凝土细观模型生成方法。在传统Voronoi图形的基础上,通过控制多面体的随机度和引入缩放因子等参数,得到具有级配的随机骨料模型。然后,对骨料几何体的外表面进行延伸,得到一定厚度的黏接界面层,从而得到完整的混凝土细观有限元模型。采用连续介质损伤模型分析混凝土材料的静态和动态力学特性,拟为该模型的进一步应用提供基础。
1. 混凝土三维细观模型
1.1 细观几何模型
三维Voronoi图形是一组由连接两顶点直线的垂直平分面形成的连续多面体。在一个特定的空间内,随机分布N个种子点,Si代表第i个胞元的核心坐标,相应的胞体由顶点P围成,顶点P到所属胞元的核心Si的距离小于或等于到其他胞元核心Sj的距离:
Vi=∩i≠j{P∈D|d(P,Si)≤d(P,Sj),i≠j} (1) d(P,S)=√(Px−Sx)2+(Py−Sy)2+(Pz−Sz)2 (2) 为了避免产生过于畸形的胞元,通过下式控制两核心之间的最小距离:
δmin (3) 式中:
\delta_ 0 为两核心之间的平均距离,根据{\delta _{{0}}} = \displaystyle\frac{{\sqrt {{6}} }}{{{2}}}{\left( {\displaystyle\frac{V}{{\sqrt {{2}} N}}} \right)}^{1/3} 计算得到;K为胞元的不规则度,当K=1时,生成胞元完全随机,当K=0时,得到完全规则的多面体。图1给出了K 取0.2时生成的随机凸多面体。对于每一个胞元,向量{ v} = \left[ {\begin{array}{*{20}{c}} {{x_{{S_i}}} - {x_{{p_i}}}} \\ {{y_{{S_i}}} - {y_{{p_i}}}} \\ {{\textit{z}_{{S_i}}} - {\textit{z}_{{p_i}}}} \end{array}} \right] 表示第i个胞元顶点到其核心Si的距离,对胞元进行缩放,得到新的顶点坐标{ p}_i' (见图2):{ p}_i' = {{ p}_i} + q { v} (4) 式中:
q\;\left( {{{0}} \text{<} q \text{<} {{1}}} \right) 为缩放因子,通过改变q,可以实现对多面体缩放大小的控制,从而得到满足一定级配的随机骨料模型(见图3)。在生成的骨料模型基础上,利用下式对多面体顶点沿着向量
v 进行延伸,并保留多面体顶点,得到顶点{{ p}_i{''}} (见图4):{ p}{_i''} = { p}_i' - \eta { v} (5) 式中:
\eta \;\left( {{{0}} \text{<} \eta \text{<} {{1}}} \right) 为控制黏接界面层厚度的延伸因子。图5为利用此方法生成的具有不同界面层厚度的骨料示意图。1.2 细观有限元模型
上述生成细观模型的方法具有过程简单、随机性可控制以及生成速度快[11]等特点,但多面体缩放会降低初始生成骨料的体积含量,达不到实际混凝土对体积分数的要求。通过重力下落过程[12]不仅可以提高模型中骨料的体积分数,同时可以改善骨料分布,增强随机性,使模型更接近真实混凝土。一定区域内不同骨料体积分数的模型如图6所示。由于实际混凝土材料中黏接界面层厚度仅为10~50 μm,远小于骨料尺寸,采用四面体划分网格数量巨大,严重耗费计算时间。因此,对于ITZ的网格划分采用楔形单元,在保留模型计算精度的同时,可以最大程度地降低网格数量,提高计算效率。图7为试样尺寸为25 mm×25 mm×25 mm、骨料体积分数为37.38%的有限元模型,其中绿色区域为骨料(四面体单元72 114个),红色区域为黏接界面层(楔形单元25 408个),界面层厚度0.1 mm,蓝色区域为砂浆(四面体单元36 727个)。
ABAQUS有限元软件中塑性损伤模型主要用于模拟砂浆、岩石和混凝土等材料在低静水压力下由损伤引起的不可恢复的材料性能,该模型考虑了拉压性能的差异,可以较好地描述砂浆和黏接界面层的力学特性[13-16]。将图7模型导入ABAQUS有限元软件,采用塑性损伤模型模拟混凝土细观组分中砂浆和黏接界面层的力学性能,砂浆塑形损伤模型曲线如图8所示,界面层相应曲线由砂浆曲线按强度进行等比例折减得到。通常骨料强度远高于砂浆强度,因此骨料采用线弹性模型。具体材料参数见表1。
表 1 细观组分材料参数Table 1. Material parameters of three-phase materials材料 弹性模量/
GPa泊松比 抗压强度/
MPa抗拉强度/
MPa砂浆 25 0.20 35 3.5 ITZ 18 0.20 20 3.0 粗骨料 43 0.23 − − 2. 三维细观模型静力学分析
2.1 单轴静力特性分析
图9为试样单轴加载示意图,试样上下两端为刚性平面,上刚性板采用位移加载,下刚性板完全固定。对于准静态问题,隐式算法不受试样惯性效应以及加载动能的影响,更能反映试样真实受力状态,本文中采用ABAQUS/standard隐式计算模块模拟准静态加载下的压缩和拉伸。单轴压缩和拉伸应力应变曲线以及相应的实验对照如图10所示。由图10可知,数值模拟与实验结果[17]曲线趋势相同,峰值误差在10%以内,吻合较好。对于单轴压缩,刚性板与试样之间的摩擦条件的不同,对混凝土试样的破坏模式有较大影响。不同摩擦条件下混凝土单轴压缩破坏模式的对比如图11所示。在低摩擦条件下,混凝土裂缝主要为平行于加载方向的竖向裂缝;在高摩擦条件下,裂缝转变为八字形裂缝,与通常实验结果[17]相似。图12为单轴拉伸载荷作用下试样的破坏模式。裂缝主要为垂直于加载方向的环向裂缝,符合实际实验结果。
2.2 双轴静力特性分析
尚世明[18]采用定测压加载方式,得到了标准立方体试样(100 mm×100 mm×100 mm的混凝土试件)在不同侧压下的宏观应力应变曲线和相应的破坏模式。图13为数值模拟试样双轴加载示意图,试样上下两端为刚性平面,下端固定,上端位移加载。侧向方向施加均匀的压应力,分别为0、8、16、24 MPa。图14为数值模拟与实验在位移加载方向的应力应变对比图。由于试样尺寸、骨料级配以及应力加载等原因的影响,导致数值模拟与实验结果存在误差,但两者趋势相同,峰值应变也较接近,可以认为误差在合理范围之内。不同侧压应力条件下试样的破坏模式如图15所示。随着侧向压应力不断增大,对试样的约束作用不断增强,受压面裂缝逐渐沿垂直于自由面方向产生,破坏形态转变为层状破坏。由图15可知,数值模拟与实验破坏模式基本吻合,说明该模型能有效反映混凝土试样在双轴受压载荷下的力学特性。
2.3 三轴静力特性分析
为了进一步研究复杂应力状态下混凝土材料的静力特性,基于本文的三维细观模型,对常规三轴条件下混凝土的应力应变关系进行了研究。图16为不同围压下试样位移加载方向应力应变曲线。如图16所示,随着围压不断增大,峰值应力增大明显,软化段逐渐模糊,当围压达到30 MPa时,会出现应力随应变增大而增大的强化段,与文献[19]的实验曲线趋势相似。
3. 三维细观模型动力学分析
混凝土作为一种典型的率相关材料,在动态载荷作用下表现出不同于静态载荷的力学特性[20-22]。运用ABAQUS显式算法,模拟低应变率下混凝土SHPB冲击破坏实验,其中入射杆和透射杆为弹性金属材料,采用弹性模型,弹性模量为210 GPa,泊松比为0.28。细观混凝土模型仍采用与准静态模拟相同的材料模型与材料参数。图17为计算得到的4种应变率下抗压强度动态增强因子。准静态模拟得到的抗压强度为24 MPa(无摩擦边界),即
f_{\rm c} = 24 MPa。从图17可以看出,抗压强度随着应变率的增大而增大,这与常规混凝土SHPB实验所得结论一致,说明模型在动态载荷作用下同样具有较好的有效性。细观模型在不同应变率下的破坏模式如图18所示,应变率较低时,宏观裂纹发展不明显,结构无明显破坏,随着应变率的增大,试样表面逐渐出现宏观裂纹,且裂纹逐渐贯穿,导致试样发生破坏。4. 结 论
在传统Voronoi图形的基础上,提出了一种简单高效的混凝土三维细观模型建立方法:即通过引入缩放因子产生具有一定级配的随机凸多面体作为骨料,随后对多面体的外表面进行延伸,得到具有一定厚度的黏接界面层,最后通过布尔运算得到包裹骨料和界面层的砂浆部分。利用塑性损伤模型对该细观模型进行准静态和动态加载下的有限元分析,得到如下结论。
(1)利用本文的三维细观模型模拟混凝土单轴、双轴以及三轴应力状态下的静力学性能,得到的应力应变曲线和破坏模式与实验结果基本一致,说明该模型能较好地反映实际混凝土在准静态载荷下的力学特性,这为进一步从细观角度研究混凝土损伤演化规律和破坏机理,提供了模型基础。
(2)混凝土在单轴压缩下的破坏模式与接触面摩擦条件有较大关系,在摩擦较低的条件下,主要产生与受载方向平行的竖向裂纹;在高摩擦条件下,破坏形态为八字型裂缝。拉应力作用下,主要产生垂直于载荷平面的环状裂纹。
(3)动态SHPB数值模拟结果表明,混凝土的抗压强度随着应变率的增大而提高。在较低应变率(50 s−1以下)下,试样外表面未发生明显损伤,随着应变率的提高,试样损伤愈发明显,外表面产生明显裂纹,导致混凝土失效破坏。
-
表 1 细观组分材料参数
Table 1. Material parameters of three-phase materials
材料 弹性模量/
GPa泊松比 抗压强度/
MPa抗拉强度/
MPa砂浆 25 0.20 35 3.5 ITZ 18 0.20 20 3.0 粗骨料 43 0.23 − − -
[1] 张楚汉, 唐欣薇, 周元德, 等. 混凝土细观力学研究进展综述 [J]. 水力发电学报, 2015, 34(12): 1–18. doi: 10.11660/slfdxb.20151201ZHANG Chuhan, TANG Xinwei, ZHOU Yuande, et al. State-of-the-art literature review on concrete meso-scale mechanics [J]. Journal of Hydroelectric Engineering, 2015, 34(12): 1–18. doi: 10.11660/slfdxb.20151201 [2] 杜修力, 金浏. 细观分析方法在混凝土物理/力学性质研究方面的应用 [J]. 水利学报, 2016, 47(3): 355–371.DU Xiuli, JIN Liu. Applications of meso-scale analysis methods on the study of the physical/mechanical properties of concrete [J]. Journal of Hydraulic Engineering, 2016, 47(3): 355–371. [3] WITTMANN F H, ROELFSTRA P E, SADOUKI H. Simulation and analysis of composite structures [J]. Materials Science and Engineering, 1985, 68(2): 239–248. [4] ZHOU X Q, HAO H. Mesoscale modelling of concrete tensile failure mechanism at high strain rates [J]. Computers and Structures, 2008, 86(21−22): 2013–2026. [5] ZHOU X Q, HAO H. Modelling of compressive behaviour of concrete-like materials at high strain rate [J]. International Journal of Solids and Structures, 2008, 45(17): 4648–4661. [6] WANG X F, YANG Z J, YATES J R, et al. Monte Carlo simulations of mesoscale fracture modelling of concrete with random aggregates and pores [J]. Construction and Building Materials, 2015, 75: 35–45. [7] 吴成, 沈晓军, 王晓鸣, 等. 细观混凝土靶抗侵彻数值模拟及侵彻深度模型 [J]. 爆炸与冲击, 2018, 38(6): 1364–1371. DOI: 10.11833/bzycj-2017-0123.WU Cheng, SHEN Xiaojun, WANG Xiaoming, et al. Numerical simulation on anti-penetration and penetration depth model of mesoscale concrete target [J]. Explosion and Shock Waves, 2018, 38(6): 1364–1371. DOI: 10.11833/bzycj-2017-0123. [8] WANG X, ZHANG M, JIVKOV A P. Computational technology for analysis of 3D meso-structure effects on damage and failure of concrete [J]. International Journal of Solids and Structures, 2016, 80: 310–333. [9] ZHANG Z, SONG X, LIU Y, et al. Three-dimensional mesoscale modelling of concrete composites by using random walking algorithm [J]. Composites Science and Technology, 2017, 149: 235–245. [10] 邓勇军, 陈小伟, 姚勇, 等. 基于细观混凝土模型的刚性弹体正侵彻弹道偏转分析 [J]. 爆炸与冲击, 2017, 37(3): 377–386. DOI: 10.11833/1001-1455(2017)03-0377-10.DENG Yongjun, CHEN Xiaowei, YAO Yong, et al. On ballistic trajectory of rigid projectile normal penetration based on a meso-scopic concrete [J]. Explosion and Shock Waves, 2017, 37(3): 377–386. DOI: 10.11833/1001-1455(2017)03-0377-10. [11] ZHANG J, WANG Z, YANG H, et al. 3D meso-scale modeling of reinforcement concrete with high volume fraction of randomly distributed aggregates [J]. Construction and Building Materials, 2018, 164: 350–361. doi: 10.1016/j.conbuildmat.2017.12.229 [12] SHENG P, ZHANG J, JI Z. An advanced 3D modeling method for concrete-like particle-reinforced composites with high volume fraction of randomly distributed particles [J]. Composites Science and Technology, 2016, 134: 26–35. [13] HUANG Y, YANG Z, REN W, et al. 3D meso-scale fracture modelling and validation of concrete based on in-situ X-ray computed tomography images using damage plasticity model [J]. International Journal of Solids and Structures, 2015, 67-68: 340–352. [14] DU X, JIN L, MA G. Numerical simulation of dynamic tensile-failure of concrete at meso-scale [J]. International Journal of Impact Engineering, 2014, 66(4): 5–17. [15] 张锦华, 方秦, 龚自明, 等. 基于三维细观模型的全级配混凝土静态力学性能的数值模拟 [J]. 计算力学学报, 2012, 29(6): 927–933.ZHANG Jinhua, FANG Qin, GONG Ziming, et al. Numerical simulation of static mechanical properties based on 3D mesoscale model of fully-graded concete [J]. Chinese Journal of Computational Mechanics, 2012, 29(6): 927–933. [16] 方秦, 还毅, 张亚栋. ABAQUS混凝土损伤塑性模型的静力性能分析 [J]. 解放军理工大学学报(自然科学版), 2007, 8(3): 254–260.FANG Qin, HUAI Yi, ZHANG Yadong, et al. Investigation into static properties of damaged plasticity model for concrete in ABAQUS [J]. Journal of PLA University of Science and Technology (Natural Science Edition), 2007, 8(3): 254–260. [17] VAN VLIET M R A, VAN MIER J G M. Softening behaviour of concrete under uniaxial compression [C] // WITTMANN F H. Fracuture mechanics of concrete structures. Freiburg: Aedificatio Publishers, 1995: 383−396. [18] 尚世明. 普通混凝土多轴动态性能试验研究[D]. 大连: 大连理工大学, 2013. [19] SFER D, CAROL I, GETTU R, et al. Study of the behavior of concrete under triaxial compression [J]. Journal of Engineering Mechanics, 2002, 128(2): 156–163. [20] 宁建国, 商霖, 孙远翔. 混凝土材料动态性能的经验公式、强度理论与唯象本构模型 [J]. 力学进展, 2006, 36(3): 389–405. doi: 10.3321/j.issn:1000-0992.2006.03.006NING Jianguo, SHANG Lin, SUN Yuanxiang. The developments of dynamic constitutive behavior of concrete [J]. Advance in Mechanics, 2006, 36(3): 389–405. doi: 10.3321/j.issn:1000-0992.2006.03.006 [21] 宁建国, 商霖, 孙远翔. 混凝土材料冲击特性的研究 [J]. 力学学报, 2006, 38(2): 199–208. doi: 10.3321/j.issn:0459-1879.2006.02.008NING Jianguo, SHANG Lin, SUN Yuanxiang. Investigation on impact behavior of concrete [J]. Chinese Journal of Theoretical and Applied Mechanics, 2006, 38(2): 199–208. doi: 10.3321/j.issn:0459-1879.2006.02.008 [22] 方秦, 洪建, 张锦华, 等. 混凝土类材料SHPB实验若干问题探讨 [J]. 工程力学, 2014, 31(5): 1–14.FANG Qin, HONG Jian, ZHANG Jinhua, et al. Issues of SHPB test on concrete-like material [J]. Engineering Mechanics, 2014, 31(5): 1–14. 期刊类型引用(10)
1. 田亮,孟俊良,赵健,樊立龙,王宇宁,张诚至. 基于内聚力模型的混凝土断裂性能三维细观模拟. 工业建筑. 2024(09): 170-176 . 百度学术
2. 姜珊,路国运,杨会伟. 侧向冲击载荷下钢管混凝土结构的动力响应及参数分析. 爆炸与冲击. 2023(11): 28-39 . 本站查看
3. 程想,曹俊鑫,王玲玲,赵银霜,孔德文. 基于细观模型的混凝土力学性能研究进展. 混凝土与水泥制品. 2022(01): 27-32+37 . 百度学术
4. 徐磊,姜磊,周昌巧,任青文. 基于多重点云与分级聚合的全级配混凝土三维细观结构高效生成方法. 水利学报. 2022(02): 188-199 . 百度学术
5. 吴宇航,肖映雄,徐亚飞. 基于Python-Abaqus的混凝土三维细观随机模型的建立. 计算力学学报. 2022(05): 566-573 . 百度学术
6. 张杰,王志华,王志勇,树学峰. 骨料对刚性弹正侵彻混凝土过程的影响机理. 中国科学:技术科学. 2021(03): 272-280 . 百度学术
7. 陈青青,张煜航,张杰,王志勇,王志华. 含孔隙混凝土二维细观建模方法研究. 应用数学和力学. 2020(02): 182-194 . 百度学术
8. 方建银,李娜,党发宁,潘优,任劼. 基于破损分区理论和CT数重建混凝土数值模型. 长安大学学报(自然科学版). 2020(02): 66-73 . 百度学术
9. 张超,刘占芳. 低压对变温环境下高聚物黏结炸药界面损伤的抑制. 应用数学和力学. 2020(10): 1057-1071 . 百度学术
10. 于华洋,马涛,王大为,王朝辉,吕松涛,朱兴一,刘鹏飞,李峰,肖月,张久鹏,罗雪,金娇,郑健龙,侯越,徐慧宁,郭猛,蒋玮. 中国路面工程学术研究综述·2020. 中国公路学报. 2020(10): 1-66 . 百度学术
其他类型引用(12)
-