• ISSN 1001-1455  CN 51-1148/O3
  • EI Compendex、CA收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊
高级检索 E-mail Alert

混凝土三维细观模型的建模方法与力学特性分析

张煜航 陈青青 张杰 王志勇 李志强 王志华

引用本文:
Citation:

混凝土三维细观模型的建模方法与力学特性分析

    作者简介: 张煜航(1993- ),男,博士研究生,zhangyh951000@163.com;
    通讯作者: 王志华, wangzh077@163.com
  • 中图分类号: O383

3D mesoscale modeling method and dynamic mechanical properties investigation of concrete

    Corresponding author: WANG Zhihua, wangzh077@163.com
  • CLC number: O383

  • 摘要: 根据混凝土材料的细观组成和结构特点,基于三维Voronoi图形提出了一种简单高效的混凝土细观模型生成方法,利用塑性损伤模型对该细观模型进行了单、多轴应力状态下的准静态分析以及SHPB动态有限元分析。结果表明,数值模拟得到的应力应变曲线和破坏模式与实验结果基本吻合,本文中提出的混凝土三维细观模型可较好地模拟混凝土的静、动态力学特性,为进一步从细观力学角度研究混凝土损伤演化规律和破坏机理提供了模型基础。
  • 图 1  K=0.2时三维Voronoi多面体

    Figure 1.  3D Voronoi polyhedron when K=0.2

    图 2  骨料缩放示意图

    Figure 2.  Schematic diagram for single aggregate

    图 3  不同大小骨料模型

    Figure 3.  Aggregate cells with grading sizes

    图 4  ITZ产生示意图

    Figure 4.  Schematic diagram of generating ITZ

    图 5  具有不同ITZ厚度的骨料外轮廓图

    Figure 5.  Geometry outlines of different ITZ layer thicknesses

    图 6  含有不同骨料体积分数的立方体试样

    Figure 6.  Cubic specimens with different volume fractions of aggregate

    图 7  三维混凝土细观模型有限元模型

    Figure 7.  Meshing results for concrete specimen

    图 8  砂浆单轴准静态压缩应力应变曲线

    Figure 8.  Strain-stress curves of mortar phase under uniaxial load

    图 9  三维模型单轴加载示意图

    Figure 9.  3D model for uniaxial quasi-static simulation

    图 10  单轴压缩与拉伸应力应变曲线

    Figure 10.  Comparison of strain-stress curves between experiment and simulation

    图 11  不同摩擦条件下试样破坏模式

    Figure 11.  Failure patterns of different frictional conditions

    图 12  单轴拉伸试样破坏模式

    Figure 12.  Failure patterns under uniaxial tension

    图 13  定测压加载示意图

    Figure 13.  3D model under biaxial simulation

    图 14  双轴压缩应力应变曲线

    Figure 14.  Comparison of strain-stress curves between experiment and simulation

    图 15  不同侧向压应力下试样破坏模式

    Figure 15.  Failure patterns of different frictional conditions

    图 16  三轴压应力作用下应力应变曲线

    Figure 16.  Comparison of stress-strain curves between experiment and simulation

    图 17  不同应变率下的动态增强因子

    Figure 17.  Dynamic increasing factor for concrete at different strain rates

    图 18  不同应变率下混凝土试样破坏模式

    Figure 18.  Failure patternsat differentstrain rates

    表 1  细观组分材料参数

    Table 1.  Material parameters of three-phase materials

    材料弹性模量/
    GPa
    泊松比抗压强度/
    MPa
    抗拉强度/
    MPa
    砂浆250.20353.5
    ITZ180.20203.0
    粗骨料430.23
    下载: 导出CSV
  • [1] 张楚汉, 唐欣薇, 周元德, 等. 混凝土细观力学研究进展综述 [J]. 水力发电学报, 2015, 34(12): 1–18. doi: 10.11660/slfdxb.20151201
    ZHANG 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.006
    NING 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.008
    NING 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.
  • [1] 张德海朱浮声邢纪波 . 梁-颗粒模型在混凝土侵彻问题中的应用. 爆炸与冲击, 2005, 25(1): 85-89. doi: 10.11883/1001-1455(2005)01-0085-05
    [2] 刘海峰宁建国 . 冲击荷载作用下混凝土材料的细观本构模型. 爆炸与冲击, 2009, 29(3): 261-267. doi: 10.11883/1001-1455(2009)03-0261-07
    [3] 邓勇军陈小伟姚勇杨涛 . 基于细观混凝土模型的刚性弹体正侵彻弹道偏转分析. 爆炸与冲击, 2017, 37(3): 377-386. doi: 10.11883/1001-1455(2017)03-0377-10
    [4] 吴成沈晓军王晓鸣姚文进 . 细观混凝土靶抗侵彻数值模拟及侵彻深度模型. 爆炸与冲击, 2018, 38(6): 1364-1371. doi: 10.11883/bzycj-2017-0123
    [5] 宋顺成才鸿年王富耻 . 关于弹体冲击和贯穿混凝土的三维数值模拟. 爆炸与冲击, 2006, 26(1): 1-6. doi: 10.11883/1001-1455(2006)01-0001-06
    [6] 徐松林王鹏飞赵坚胡时胜 . 基于三维Hopkinson杆的混凝土动态力学性能研究. 爆炸与冲击, 2017, 37(2): 180-185. doi: 10.11883/1001-1455(2017)02-0180-06
    [7] 张磊胡时胜陈德兴张守保余泽清刘飞 . 混凝土材料的层裂特性. 爆炸与冲击, 2008, 28(3): 193-199. doi: 10.11883/1001-1455(2008)03-0193-07
    [8] 闫东明林皋 . 单向恒定侧压下混凝土动态抗压特性试验研究. 爆炸与冲击, 2007, 27(2): 121-125. doi: 10.11883/1001-1455(2007)02-0121-05
    [9] 葛涛刘保荣王明洋 . 弹体侵彻与贯穿有限厚度混凝土靶体的力学特性. 爆炸与冲击, 2010, 30(2): 159-163. doi: 10.11883/1001-1455(2010)02-0159-05
    [10] 王永刚王礼立 . 平板撞击下C30混凝土中冲击波的传播特性. 爆炸与冲击, 2010, 30(2): 119-124. doi: 10.11883/1001-1455(2010)02-0119-06
    [11] 胡宏伟冯海云陈朗顾晓辉宋浦 . 非理想炸药在混凝土介质中的爆炸做功特性. 爆炸与冲击, 2018, 38(1): 197-203. doi: 10.11883/bzycj-2016-0123
    [12] 巫绪涛孙善飞李和平 . 用HJC本构模型模拟混凝土SHPB实验. 爆炸与冲击, 2009, 29(2): 137-142. doi: 10.11883/1001-1455(2009)02-0137-06
    [13] 林华令丁育青汤文辉 . 混凝土侵彻数值模拟的影响因素. 爆炸与冲击, 2013, 33(4): 425-429. doi: 10.11883/1001-1455(2013)04-0425-05
    [14] 吴昊方秦龚自明 . 考虑刚性弹弹头形状的混凝土(岩石)靶体侵彻深度半理论分析. 爆炸与冲击, 2012, 32(6): 573-580. doi: 10.11883/1001-1455(2012)06-0573-08
    [15] 李志康黄风雷 . 混凝土材料的动态空腔膨胀理论. 爆炸与冲击, 2009, 29(1): 95-100. doi: 10.11883/1001-1455(2009)01-0095-06
    [16] 杨冬梅王晓鸣 . 混凝土中爆炸数值仿真算法研究. 爆炸与冲击, 2005, 25(6): 569-573. doi: 10.11883/1001-1455(2005)06-0569-05
    [17] 高光发郭扬波 . 高强混凝土动态压缩试验分析. 爆炸与冲击, 2019, 39(3): 033103-1-033103-10. doi: 10.11883/bzycj-2017-0405
    [18] 陶俊林秦李波李奎刘丹贾彬陈小伟陈刚 . 混凝土高温动态压缩力学性能实验. 爆炸与冲击, 2011, 31(1): 101-106. doi: 10.11883/1001-1455(2011)01-0101-06
    [19] 徐伟芳张方举陈裕泽毛勇建任时成陶俊林 . 混凝土薄板的侵彻破坏响应实验. 爆炸与冲击, 2013, 33(2): 169-174. doi: 10.11883/1001-1455(2013)02-0169-06
    [20] 何翔徐翔云孙桂娟沈俊杨建超金栋梁 . 弹体高速侵彻混凝土的效应实验. 爆炸与冲击, 2010, 30(1): 1-6. doi: 10.11883/1001-1455(2010)01-0001-06
  • 加载中
图(18)表(1)
计量
  • 文章访问数:  12
  • HTML全文浏览量:  131
  • PDF下载量:  11
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-10-26
  • 录用日期:  2018-12-10
  • 网络出版日期:  2019-04-25
  • 刊出日期:  2019-05-01

混凝土三维细观模型的建模方法与力学特性分析

    作者简介:张煜航(1993- ),男,博士研究生,zhangyh951000@163.com
    通讯作者: 王志华, wangzh077@163.com
  • 1. 太原理工大学机械与运载工程学院应用力学研究所,山西 太原 030024
  • 2. 太原理工大学材料强度与结构冲击山西省重点实验室,山西 太原 030024

摘要: 根据混凝土材料的细观组成和结构特点,基于三维Voronoi图形提出了一种简单高效的混凝土细观模型生成方法,利用塑性损伤模型对该细观模型进行了单、多轴应力状态下的准静态分析以及SHPB动态有限元分析。结果表明,数值模拟得到的应力应变曲线和破坏模式与实验结果基本吻合,本文中提出的混凝土三维细观模型可较好地模拟混凝土的静、动态力学特性,为进一步从细观力学角度研究混凝土损伤演化规律和破坏机理提供了模型基础。

English Abstract

  • 混凝土材料从细观层次上由骨料、砂浆以及两者之间的黏接界面层(interfacial transitional zone,ITZ)组成,各组分的细观力学性质直接或间接地决定了混凝土材料的宏观力学性能和破坏模式[1-2]。根据混凝土细观结构特点及其力学性能,深入研究细观结构和宏观力学性能之间的关系,是混凝土力学性质的重要研究内容之一。

    Wittmann等[3]首先将多尺度的研究方法应用到混凝土材料的研究中,认为混凝土是一种非均质复合材料,提出了用随机分布的多边形来模拟混凝土骨料的细观力学模型。该方法一直被后来的研究者所采用。Zhou等[4-5]建立了圆形随机骨料模型,利用材料损伤本构关系,系统研究了混凝土细观模型在不同受力条件下的拉伸、压缩和爆炸问题,指出骨料和砂浆之间的黏接界面层在压缩时先出现破坏裂纹,在拉伸时对混凝土失效机理和拉伸强度影响最大。Wang等[6]在椭圆形和多边形骨料的基础上引入孔隙,并研究了骨料形状、骨料体积分数和孔隙率对抗拉强度的影响。吴成等[7]对刚性弹丸侵彻细观混凝土靶进行了数值模拟,分析了砂浆种类、粗骨料种类和粗骨料体积分数等对靶板抗侵彻能力的影响,并通过扩展Forrestal阻力方程,建立了细观混凝土侵彻深度模型。相比于二维模型,三维细观模型更接近于实际混凝土的真实形态,因此能较好地反映混凝土材料的实际变形与损伤破坏情况。Wang等[8]提出了“投放算法”用于产生随机分布的三维骨料模型,并在单元之间嵌入零厚度黏接单元,以此来模拟混凝土在拉应力作用下的开裂行为。Zhang等[9]通过对生成的骨料进行平移和旋转来控制骨料的移动,提出了“随机爬行算法”,该算法可以实现更高骨料含量的混凝土建模。邓勇军等[10]基于骨料随机投放的思想建立了混凝土三维细观几何模型,分析了刚性弹正侵彻过程中发生弹道偏转的原因及可能的影响因素,定量研究了混凝土细观因素对弹道偏转的影响。然而,骨料投放往往需要复杂的算法对骨料相互之间是否发生侵入进行判断,同时需要花费较多的时间,建模效率较低。三维Voronoi图形中胞元形状与混凝土骨料具有良好的相似性,其胞元可实现随机的空间形状和空间分布,可作为建立骨料的一种方法。

    本文中基于三维Voronoi图形,提出一种简单、高效的混凝土细观模型生成方法。在传统Voronoi图形的基础上,通过控制多面体的随机度和引入缩放因子等参数,得到具有级配的随机骨料模型。然后,对骨料几何体的外表面进行延伸,得到一定厚度的黏接界面层,从而得到完整的混凝土细观有限元模型。采用连续介质损伤模型分析混凝土材料的静态和动态力学特性,拟为该模型的进一步应用提供基础。

    • 三维Voronoi图形是一组由连接两顶点直线的垂直平分面形成的连续多面体。在一个特定的空间内,随机分布N个种子点,Si代表第i个胞元的核心坐标,相应的胞体由顶点P围成,顶点P到所属胞元的核心Si的距离小于或等于到其他胞元核心Sj的距离:

      ${V_i} = \underset{i\ne j}{\mathop{\cap }}{} \left\{ {P \in D\;\left| {\;d\left( {P, {S_i}} \right) \text{≤}} \right.\;d\left( {P, {S_j}} \right)\;, \;i \ne j} \right\}$

      $d\left( {P, S} \right) = \sqrt {{{\left( {{P_x} - {S_x}} \right)}^2} + {{\left( {{P_y} - {S_y}} \right)}^2} + {{\left( {{P_z} - {S_z}} \right)}^2}} $

      为了避免产生过于畸形的胞元,通过下式控制两核心之间的最小距离:

      ${\delta _{\min }} = \left( {{{1}} - K} \right){\delta _{{0}}}$

      式中:$\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[ {{array}{*{20}{c}} {{x_{{S_i}}} - {x_{{p_i}}}} \\ {{y_{{S_i}}} - {y_{{p_i}}}} \\ {{\textit{z}_{{S_i}}} - {\textit{z}_{{p_i}}}} {array}} \right]$表示第i个胞元顶点到其核心Si的距离,对胞元进行缩放,得到新的顶点坐标${ p}_i'$(见图2):

      图  1  K=0.2时三维Voronoi多面体

      Figure 1.  3D Voronoi polyhedron when K=0.2

      图  2  骨料缩放示意图

      Figure 2.  Schematic diagram for single aggregate

      $ { p}_i' = {{ p}_i} + q { v} $

      式中:$q\;\left( {{{0}} \text{<} q \text{<} {{1}}} \right)$为缩放因子,通过改变q,可以实现对多面体缩放大小的控制,从而得到满足一定级配的随机骨料模型(见图3)。

      图  3  不同大小骨料模型

      Figure 3.  Aggregate cells with grading sizes

      在生成的骨料模型基础上,利用下式对多面体顶点沿着向量$ v $进行延伸,并保留多面体顶点,得到顶点${{ p}_i{''}}$(见图4):

      图  4  ITZ产生示意图

      Figure 4.  Schematic diagram of generating ITZ

      $ { p}{_i''} = { p}_i' - \eta { v} $

      式中:$\eta \;\left( {{{0}} \text{<} \eta \text{<} {{1}}} \right)$为控制黏接界面层厚度的延伸因子。图5为利用此方法生成的具有不同界面层厚度的骨料示意图。

      图  5  具有不同ITZ厚度的骨料外轮廓图

      Figure 5.  Geometry outlines of different ITZ layer thicknesses

    • 上述生成细观模型的方法具有过程简单、随机性可控制以及生成速度快[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个)。

      图  6  含有不同骨料体积分数的立方体试样

      Figure 6.  Cubic specimens with different volume fractions of aggregate

      图  7  三维混凝土细观模型有限元模型

      Figure 7.  Meshing results for concrete specimen

      ABAQUS有限元软件中塑性损伤模型主要用于模拟砂浆、岩石和混凝土等材料在低静水压力下由损伤引起的不可恢复的材料性能,该模型考虑了拉压性能的差异,可以较好地描述砂浆和黏接界面层的力学特性[13-16]。将图7模型导入ABAQUS有限元软件,采用塑性损伤模型模拟混凝土细观组分中砂浆和黏接界面层的力学性能,砂浆塑形损伤模型曲线如图8所示,界面层相应曲线由砂浆曲线按强度进行等比例折减得到。通常骨料强度远高于砂浆强度,因此骨料采用线弹性模型。具体材料参数见表1

      图  8  砂浆单轴准静态压缩应力应变曲线

      Figure 8.  Strain-stress curves of mortar phase under uniaxial load

      材料弹性模量/
      GPa
      泊松比抗压强度/
      MPa
      抗拉强度/
      MPa
      砂浆250.20353.5
      ITZ180.20203.0
      粗骨料430.23

      表 1  细观组分材料参数

      Table 1.  Material parameters of three-phase materials

    • 图9为试样单轴加载示意图,试样上下两端为刚性平面,上刚性板采用位移加载,下刚性板完全固定。对于准静态问题,隐式算法不受试样惯性效应以及加载动能的影响,更能反映试样真实受力状态,本文中采用ABAQUS/standard隐式计算模块模拟准静态加载下的压缩和拉伸。单轴压缩和拉伸应力应变曲线以及相应的实验对照如图10所示。由图10可知,数值模拟与实验结果[17]曲线趋势相同,峰值误差在10%以内,吻合较好。对于单轴压缩,刚性板与试样之间的摩擦条件的不同,对混凝土试样的破坏模式有较大影响。不同摩擦条件下混凝土单轴压缩破坏模式的对比如图11所示。在低摩擦条件下,混凝土裂缝主要为平行于加载方向的竖向裂缝;在高摩擦条件下,裂缝转变为八字形裂缝,与通常实验结果[17]相似。图12为单轴拉伸载荷作用下试样的破坏模式。裂缝主要为垂直于加载方向的环向裂缝,符合实际实验结果。

      图  9  三维模型单轴加载示意图

      Figure 9.  3D model for uniaxial quasi-static simulation

      图  10  单轴压缩与拉伸应力应变曲线

      Figure 10.  Comparison of strain-stress curves between experiment and simulation

      图  11  不同摩擦条件下试样破坏模式

      Figure 11.  Failure patterns of different frictional conditions

      图  12  单轴拉伸试样破坏模式

      Figure 12.  Failure patterns under uniaxial tension

    • 尚世明[18]采用定测压加载方式,得到了标准立方体试样(100 mm×100 mm×100 mm的混凝土试件)在不同侧压下的宏观应力应变曲线和相应的破坏模式。图13为数值模拟试样双轴加载示意图,试样上下两端为刚性平面,下端固定,上端位移加载。侧向方向施加均匀的压应力,分别为0、8、16、24 MPa。图14为数值模拟与实验在位移加载方向的应力应变对比图。由于试样尺寸、骨料级配以及应力加载等原因的影响,导致数值模拟与实验结果存在误差,但两者趋势相同,峰值应变也较接近,可以认为误差在合理范围之内。不同侧压应力条件下试样的破坏模式如图15所示。随着侧向压应力不断增大,对试样的约束作用不断增强,受压面裂缝逐渐沿垂直于自由面方向产生,破坏形态转变为层状破坏。由图15可知,数值模拟与实验破坏模式基本吻合,说明该模型能有效反映混凝土试样在双轴受压载荷下的力学特性。

      图  13  定测压加载示意图

      Figure 13.  3D model under biaxial simulation

      图  14  双轴压缩应力应变曲线

      Figure 14.  Comparison of strain-stress curves between experiment and simulation

      图  15  不同侧向压应力下试样破坏模式

      Figure 15.  Failure patterns of different frictional conditions

    • 为了进一步研究复杂应力状态下混凝土材料的静力特性,基于本文的三维细观模型,对常规三轴条件下混凝土的应力应变关系进行了研究。图16为不同围压下试样位移加载方向应力应变曲线。如图16所示,随着围压不断增大,峰值应力增大明显,软化段逐渐模糊,当围压达到30 MPa时,会出现应力随应变增大而增大的强化段,与文献[19]的实验曲线趋势相似。

      图  16  三轴压应力作用下应力应变曲线

      Figure 16.  Comparison of stress-strain curves between experiment and simulation

    • 混凝土作为一种典型的率相关材料,在动态载荷作用下表现出不同于静态载荷的力学特性[20-22]。运用ABAQUS显式算法,模拟低应变率下混凝土SHPB冲击破坏实验,其中入射杆和透射杆为弹性金属材料,采用弹性模型,弹性模量为210 GPa,泊松比为0.28。细观混凝土模型仍采用与准静态模拟相同的材料模型与材料参数。图17为计算得到的4种应变率下抗压强度动态增强因子。准静态模拟得到的抗压强度为24 MPa(无摩擦边界),即$f_{\rm c} = $24 MPa。从图17可以看出,抗压强度随着应变率的增大而增大,这与常规混凝土SHPB实验所得结论一致,说明模型在动态载荷作用下同样具有较好的有效性。细观模型在不同应变率下的破坏模式如图18所示,应变率较低时,宏观裂纹发展不明显,结构无明显破坏,随着应变率的增大,试样表面逐渐出现宏观裂纹,且裂纹逐渐贯穿,导致试样发生破坏。

      图  17  不同应变率下的动态增强因子

      Figure 17.  Dynamic increasing factor for concrete at different strain rates

      图  18  不同应变率下混凝土试样破坏模式

      Figure 18.  Failure patternsat differentstrain rates

    • 在传统Voronoi图形的基础上,提出了一种简单高效的混凝土三维细观模型建立方法:即通过引入缩放因子产生具有一定级配的随机凸多面体作为骨料,随后对多面体的外表面进行延伸,得到具有一定厚度的黏接界面层,最后通过布尔运算得到包裹骨料和界面层的砂浆部分。利用塑性损伤模型对该细观模型进行准静态和动态加载下的有限元分析,得到如下结论。

      (1)利用本文的三维细观模型模拟混凝土单轴、双轴以及三轴应力状态下的静力学性能,得到的应力应变曲线和破坏模式与实验结果基本一致,说明该模型能较好地反映实际混凝土在准静态载荷下的力学特性,这为进一步从细观角度研究混凝土损伤演化规律和破坏机理,提供了模型基础。

      (2)混凝土在单轴压缩下的破坏模式与接触面摩擦条件有较大关系,在摩擦较低的条件下,主要产生与受载方向平行的竖向裂纹;在高摩擦条件下,破坏形态为八字型裂缝。拉应力作用下,主要产生垂直于载荷平面的环状裂纹。

      (3)动态SHPB数值模拟结果表明,混凝土的抗压强度随着应变率的增大而提高。在较低应变率(50 s−1以下)下,试样外表面未发生明显损伤,随着应变率的提高,试样损伤愈发明显,外表面产生明显裂纹,导致混凝土失效破坏。

参考文献 (22)

目录

    /

    返回文章
    返回