三级配混凝土单轴破坏的细观数值模拟
- 格式:pdf
- 大小:322.69 KB
- 文档页数:4
《基于三维细观模型的混凝土损伤力学行为研究》篇一一、引言混凝土作为建筑结构的主要材料,其力学性能的研究对于保障建筑安全具有重要意义。
混凝土损伤力学行为的研究是该领域的重要方向之一,而基于三维细观模型的混凝土损伤力学行为研究更是近年来研究的热点。
本文旨在通过建立三维细观模型,对混凝土损伤力学行为进行深入研究,以期为混凝土结构的设计和施工提供理论支持。
二、混凝土三维细观模型的建立建立混凝土三维细观模型是研究混凝土损伤力学行为的基础。
该模型应包含混凝土的基本组成成分,如骨料、砂浆和气孔等。
在模型中,需要考虑到各种组成成分的形状、大小、分布和排列方式等因素,以及它们之间的相互作用关系。
目前,随着计算机技术的不断发展,通过数值模拟方法建立混凝土三维细观模型已成为可能。
在建立混凝土三维细观模型时,需要考虑到模型的尺度问题。
由于混凝土结构的尺度较大,因此在建模时需要选择合适的尺度范围。
同时,需要考虑模型的精细度问题,即在保证计算效率的同时尽可能地提高模型的精度。
三、混凝土损伤力学行为的研究方法混凝土损伤力学行为的研究方法主要包括实验研究和数值模拟两种方法。
实验研究可以通过对混凝土试件进行加载、卸载等操作,观察其力学性能的变化,从而得出混凝土的损伤规律。
而数值模拟则可以通过建立混凝土三维细观模型,利用有限元等方法对混凝土进行力学分析,从而得出混凝土的损伤情况。
在研究混凝土损伤力学行为时,需要考虑到混凝土的多种因素,如骨料类型、砂浆性质、气孔分布等。
因此,在实验研究和数值模拟中,需要设计多种方案,以全面了解混凝土的损伤规律。
四、基于三维细观模型的混凝土损伤力学行为研究基于三维细观模型的混凝土损伤力学行为研究,可以通过对模型进行力学分析,得出混凝土的损伤情况。
在分析过程中,需要考虑到混凝土的多种因素,如骨料形状、大小和分布,砂浆的强度和弹性模量等。
同时,还需要考虑到混凝土的应力状态、加载速率等因素对损伤的影响。
在研究中,可以采用有限元等方法对模型进行力学分析。
第15卷 第12期2020年12月中国科技论文CHINASCIENCEPAPERVol.15No.12Dec.2020自密实混凝土单轴压缩试验细观模拟王晓琴1,杨名超2,肖 明3,程瑾俊4,孙辅皓5(1.武汉科技大学城市学院,武汉430083;2.中国市政工程中南设计研究总院有限公司,武汉430010;3.武汉绕城高速公路管理处,武汉430000;4.安庆职业技术学院,安徽安庆246003;5.江苏大彭工程项目管理咨询有限公司,江苏徐州221300)摘 要:为了研究自密实混凝土(self compactingconcrete,SCC)力学性能及破坏机理,采用离散元方法,建立了SCC单轴压缩试验离散元模型,进行了不同骨料形状、不同骨料占比和界面过渡区(interfacialtransitionzone,ITZ)不同强度下单轴压缩模拟试验。
结果表明:建立SCC离散元模型可以很好地描述其裂纹的产生、发展和破坏过程;不规则骨料颗粒比圆形骨料颗粒试件更能反映真实骨料分布,碎石骨料抗压强度更接近室内试验;ITZ强度对试件的抗压强度有一定的影响,当ITZ强度大于砂浆强度80%时,试件抗压强度基本不受ITZ强度的影响。
关键词:自密实混凝土;细观模拟;离散元;单轴压缩;界面过渡区中图分类号:TU528.5 文献标志码:A文章编号:20952783(2020)12141007开放科学(资源服务)标识码(OSID):犕犲狊狅 狊犮犪犾犲狊犻犿狌犾犪狋犻狅狀狅狀狌狀犻犪狓犻犪犾犮狅犿狆狉犲狊狊犻狅狀狋犲狊狋狅犳狊犲犾犳 犮狅犿狆犪犮狋犻狀犵犮狅狀犮狉犲狋犲WANGXiaoqin1,YANGMingchao2,XIAOMing3,CHENGJinjun4,SUNFuhao5(1.犆犻狋狔犆狅犾犾犲犵犲,犠狌犺犪狀犝狀犻狏犲狉狊犻狋狔狅犳犛犮犻犲狀犮犲犪狀犱犜犲犮犺狀狅犾狅犵狔,犠狌犺犪狀430083,犆犺犻狀犪;2.犆犲狀狋狉犪犾犪狀犱犛狅狌狋犺犲犿犆犺犻狀犪犕狌狀犻犮犻狆犪犾犈狀犵犻狀犲犲狉犻狀犵犇犲狊犻犵狀牔犚犲狊犲犪狉犮犺犐狀狊狋犻狋狌狋犲犆狅.,犔狋犱.,犠狌犺犪狀430010,犆犺犻狀犪;3.犠狌犺犪狀犚狅狌狀犱 狋犺犲 犆犻狋狔犈狓狆狉犲狊狊狑犪狔犃犱犿犻狀犻狊狋狉犪狋犻狅狀犗犳犳犻犮犲,犠狌犺犪狀430000,犆犺犻狀犪;4.犃狀狇犻狀犵犞狅犮犪狋犻狅狀犪犾犪狀犱犜犲犮犺狀犻犮犪犾犆狅犾犾犲犵犲,犃狀狇犻狀犵,犃狀犺狌犻246003,犆犺犻狀犪;5.犑犻犪狀犵狊狌犇犪狆犲狀犵犈狀犵犻狀犲犲狉犻狀犵犘狉狅犼犲犮狋犕犪狀犪犵犲犿犲狀狋犆狅狀狊狌犾狋犻狀犵犆狅.,犔狋犱.,犡狌狕犺狅狌,犑犻犪狀犵狊狌221300,犆犺犻狀犪)犃犫狊狋狉犪犮狋:Inordertostudythemechanicalpropertiesandfailuremechanismofself compactingconcrete(SCC),theuniaxialcom pressiontestsofSCCweresimulatedwithdifferentaggregateshapes,proportionsandthestrengthofinterfacialtransitionzone(ITZ)accordingtoestablisheddiscreteelementmodel.Thetestresultsshowthattheproposedmodelisfeasibletodescribethegeneration,developmentandfailureprocessofSCCcracks.Theirregularaggregateparticlescanreflecttherealaggregatedistri butionbetterthantheroundones,andthecompressivestrengthisclosertothelaboratorytest.ITZstrengthhasacertaininflu enceonthecompressivestrengthofthespecimen.WhenITZstrengthisgreaterthan80%ofmortarstrength,thecompressivestrengthofthespecimenisnotaffectedbyITZstrength.犓犲狔狑狅狉犱狊:self compactingconcrete(SCC);meso scalesimulation;discreteelement;uniaxialcompression;interfacialtransitionzone(ITZ)收稿日期:2020 03 11基金项目:湖北省教育厅科学研究计划指导性项目(B2017425);湖北省高校省级教学改革研究项目(2018512);安徽高校自然科学研究重点项目(KJ2017A787,KJ2019A1193);安徽省高等学校省级质量工程项目(2016xn2x025)第一作者:王晓琴(1978—),女,副教授,主要研究方向为混凝土材料和混凝土结构,751658171@qq.com 自密实混凝土(self compactingconcrete,SCC)作为高性能混凝土之一,无需振捣就可达到密实效果,硬化后具有抗裂能力强、抗渗透能力高等特点,可有效避免外界环境的影响,因而被广泛应用到各种实际工程中。
混凝土组份的数值模拟及其性能分析混凝土是一种重要的建筑材料,其强度、耐久性和可靠性都对建筑结构的承载能力和使用寿命产生重要影响。
为了更好地了解混凝土的性能,采用数值模拟技术可以预测混凝土的力学性能,从而更好地设计和优化混凝土结构,提高其可靠性和耐久性。
本文将对混凝土组份的数值模拟及其性能分析进行探讨。
一、混凝土组份的数值模拟1.1 有限元模拟方法有限元模拟方法是混凝土组份的数值模拟中最广泛应用的方法之一。
该方法基于物质力学原理和数字计算方法,将混凝土结构划分为离散的网格单元,并在每个网格内求解相应的方程组。
通过有限元分析,可以得到混凝土的应力、应变、变形和裂纹等相关信息,从而预测其力学性能。
1.2 计算流体力学模拟方法计算流体力学模拟方法是近年来发展起来的一种混凝土组份数值模拟方法。
该方法将混凝土视为一种流体,采用流体动力学算法求解其空气、水、固体颗粒等的流动和相互作用过程,在混凝土的力学性能预测中发挥着越来越重要的作用。
1.3 离散元模拟方法离散元模拟方法是一种基于微观颗粒运动学原理的混凝土组份数值模拟方法。
该方法将混凝土视为由许多颗粒组成的离散体系,并在每个颗粒间求解相互作用的力学和动力学方程。
通过离散元分析,可以得到混凝土颗粒的运动学、动力学和受力状态等信息,从而预测其力学性能。
二、混凝土性能分析2.1 压缩强度和抗拉强度压缩强度和抗拉强度是混凝土性能中最基本的两个指标之一,通常通过实验测量来确定。
数值模拟方法可以预测混凝土的压缩强度和抗拉强度,从而优化其配合比、控制其强度和耐久性。
2.2 变形和裂纹混凝土材料的变形和裂纹是混凝土结构在使用过程中非常重要的性能指标之一。
通过数值模拟方法,可以预测混凝土的变形和裂纹行为,从而优化其配合比和结构设计,提高混凝土结构的耐久性和可靠性。
2.3 疲劳性能混凝土材料在长期使用过程中经历的疲劳循环应力作用对其性能有着重要的影响。
通过数值模拟方法,可以预测混凝土在长期循环载荷下的疲劳性能,为结构耐久性的分析和改进提供依据和指导。
2。
3。
2 细观力学与混凝土破坏由混凝土材料组成的工程结构,如高坝、桥梁、海洋平台、核电站、隧道、地基基础及边坡等是基础设施建设中重要的组成部分. 混凝土材料是以水泥为主要胶结材料, 拌合一定比例的砂、石和水,经过搅拌、振捣、养护等工序后,逐渐凝固硬化而成的复合材料. 粗骨料和硬化水泥砂浆两种主要组成材料的成分、性质、配比以及粘结作用均对混凝土的力学特性有不同程度的影响,这使混凝土比其他单一材料具有更为复杂的力学性能.混凝土力学特性(宏观应力–应变关系和破坏机制) 是进行大坝、海洋平台、边坡等混凝土结构抗震设计及静、动力仿真分析的重要基础之一,也是目前研究的薄弱环节.混凝土是由水、水泥和粗细骨料组成的复合材料。
一般从特征尺寸和研究方法的侧重点不同将混凝土内部结构分为三个层次[24](如图2—2):(1)微观层次(Micro—level).材料的结构单元尺度在原子、分子量级,即从小于10-7cm~10—4cm着眼于水泥水化物的微观结构分析。
由晶体结构及分子结构组成,可用电子显微镜观察分析,是材料科学的研究对象;(2)细观层次(Meso-level)。
从分子尺度到宏观尺度,其结构单元尺度变化范围在10—4厘米至几厘米,或更大些,着眼于粗细骨料、水泥水化物、孔隙、界面等细观结构,组成多相复合材料,可按各类计算模型进行数值分析。
在这个层次上,混凝土被认为是一种由粗骨料、硬化水泥砂浆和它们之间的过渡区(粘结带)组成的三相材料。
砂浆中的孔隙很小而量多,且随机分布,水泥砂浆力学性能可以看作细观均质损伤体。
相同配合比、相同条件的砂浆试件,通常其力学性能也比较稳定,可以由试验直接测定。
由泌水、干缩和温度变化引起粗骨料和水泥砂浆之间产生初始粘结裂缝,而这些细观内部裂隙的发展将直接影响混凝土的宏观力学性能;(3)宏观层次(Macro-level)。
特征尺寸大于几厘米,混凝土作为非均质材料存在着一种特征体积,一般认为是相当于3~4倍的最大骨料体积。
混凝土抗压强度尺寸效应的细观数值模拟的开题报告一、选题背景混凝土在现代建筑和工程中具有广泛的应用,其中最基本的力学性质是其抗压强度。
在实际工程中,混凝土的尺寸及其不同的应力状态会对其抗压强度产生影响。
目前,国内外学术界及工程领域已经形成了比较完整的混凝土尺寸效应研究成果,但尺寸效应的细观数值模拟方法研究相对较少。
基于此,本研究拟通过有限元方法,结合理论分析和数值模拟技术,针对混凝土尺寸效应问题,开展深入研究,期望为混凝土工程设计提供有力的理论支持和指导。
二、研究目的本研究旨在通过对混凝土尺寸效应问题的细观数值模拟,探究混凝土尺寸对其抗压强度的影响规律,分析不同尺寸下混凝土结构的力学性能特点,为混凝土结构设计提供实用的依据和参考。
具体研究任务包括:1. 综合国内外尺寸效应研究成果,建立基本的混凝土尺寸效应理论模型,分析其基本原理、规律和特点;2. 建立二维和三维混凝土有限元数值模拟模型,设计合理的数值计算实验方案,通过对不同尺寸混凝土试件进行抗压试验,获取试验数据,验证数值模拟结果的准确性和可靠性;3. 在有限元软件中完成对不同尺寸混凝土试件的细观数值模拟,分析不同尺寸下混凝土的应力分布、破坏形态和抗压性能指标等,比较不同尺寸下混凝土试件力学性质的差异;4. 提出混凝土尺寸效应的设计规范建议,为工程实际应用提供理论基础和指导思路。
三、研究方法本研究主要采用理论分析和数值模拟相结合的方法,旨在通过建立混凝土尺寸效应的理论模型,结合有限元数值模拟技术,对不同尺寸的混凝土试件进行细观数值模拟,并分析其力学特点和性能差异。
具体的研究步骤包括:1. 理论分析:综合相关文献,建立混凝土尺寸效应的理论模型,分析其基本原理和规律,探究混凝土抗压强度与尺寸的关系;2. 数值模拟建模:根据建立的理论模型,利用有限元软件(如ABAQUS等),建立二维和三维混凝土有限元模型,设置材料参数、边界条件,进行实际计算;3. 计算实验:设计合理的数值模拟实验方案,仿真不同尺寸下混凝土试件的抗压试验,获取试验数据,验证数值模拟结果的准确性和可靠性;4. 数值模拟:在有限元软件中运行数值模拟计算,通过对不同尺寸混凝土试件应力分布、破坏形态和抗压性能等方面的分析,探究混凝土尺寸效应的影响规律和特点;5. 结果分析与评价:将数字模拟试验和实验结果进行比较,得出混凝土尺寸效应的结论,分析不同尺寸下混凝土结构的力学特性差异,并提出混凝土尺寸效应的设计规范建议。
超高性能混凝土细观毁伤有限元数值模拟及试验探究关键词:超高性能混凝土;细观有限元模型;力学性能;毁伤特性超高性能混凝土细观毁伤有限元数值模拟及试验探究引言超高性能混凝土由于其高强度、高流淌性、高耐久性等特性,近年来成为了建筑工程领域中备受关注的材料之一。
然而,超高性能混凝土的力学性能和毁伤特性对于其实际应用而言至关重要。
因此,本文旨在探究超高性能混凝土的力学性能及毁伤特性,为建筑材料领域的探究提供参考。
试验材料及细观有限元模型试验材料:本试验选取了一种常用的自然石英粉、高性能水泥、微细粉末等为主要原料,制备了一种超高性能混凝土材料。
细观有限元模型:本文构建了一个三维细观有限元模型,对超高性能混凝土的各向异性进行了分析。
在模型中,接受了非线性弹性理论和各向异性本构干系,思量了混凝土材料的各向异性和非线性变形特性。
试验方法试验包括压缩试验、拉伸试验和剪切试验。
在试验中,实行了标准试样制备方法和设备进行试验,记录了材料在不同荷载下的变形和应力等数据。
通过试验结果,探究了材料的应力-应变特性和毁伤行为。
结果分析通过试验数据,分析了超高性能混凝土的力学性能。
试验结果表明,所制超高性能混凝土的抗压强度、抗拉强度和抗剪强度均高于平凡混凝土,并且具有较好的变形能力。
通过细观有限元模型的分析,得出了混凝土各向异性的分布和变化规律。
此外,试验还反映了超高性能混凝土的初始毁伤阶段和主要毁伤阶段的毁伤特点。
结论通过细观有限元模型和试验,探究了超高性能混凝土的力学性能和毁伤特性。
结果表明,所制超高性能混凝土具有高强度、高韧性、高耐久性等特性,可以广泛应用于各种建筑工程中。
但是,还需要进一步改善其力学性能和毁伤韧性,以满足复杂工程要求。
此外,还需要注意控制超高性能混凝土的生产过程,以确保材料的匀称性和一致性。
在实际施工中,也需要注意超高性能混凝土的施工工艺和保卫措施,以充分发挥其优异的性能。
另外,还有一些探究表明,通过添加纳米材料,如氧化硅、氧化铝等,可以进一步提杰出高性能混凝土的力学性能和耐久性。
《基于三维细观模型的混凝土损伤力学行为研究》一、引言混凝土作为建筑工程中最重要的材料之一,其力学性能和损伤行为的研究一直是学术界和工程界关注的焦点。
传统的混凝土损伤力学研究大多基于宏观尺度上的实验观测和理论分析,然而这些方法无法准确描述混凝土内部复杂的细观结构和损伤机制。
近年来,随着计算机技术的快速发展,基于三维细观模型的研究方法逐渐成为混凝土损伤力学研究的新趋势。
本文旨在通过建立三维细观模型,研究混凝土的损伤力学行为,为混凝土结构的设计和施工提供理论依据。
二、三维细观模型的建立为了准确描述混凝土内部的细观结构,我们建立了基于三维空间的三细观模型。
该模型考虑了混凝土内部的骨料、砂浆和孔隙等组成部分,通过数字图像处理技术和随机骨料模型生成技术,构建了真实的三维细观结构模型。
在此基础上,我们采用了离散元方法对模型进行数值计算和分析。
三、混凝土的损伤力学行为研究在三维细观模型的基础上,我们通过施加外部荷载和内部因素的作用,研究混凝土的损伤力学行为。
具体而言,我们主要关注以下几个方面:1. 骨料与砂浆的相互作用:骨料和砂浆是混凝土的主要组成部分,它们之间的相互作用对混凝土的力学性能具有重要影响。
我们通过数值计算,分析了骨料和砂浆在受力过程中的相互作用机制,揭示了混凝土在荷载作用下的变形和损伤规律。
2. 混凝土内部的裂纹扩展:混凝土的损伤主要表现为裂纹的扩展和贯通。
我们通过观察和分析模型中的裂纹扩展过程,探讨了混凝土损伤的机理和规律。
同时,我们还研究了不同因素对裂纹扩展的影响,如骨料类型、骨料含量、砂浆的强度等。
3. 混凝土在动态荷载下的响应:在实际工程中,混凝土结构经常受到动态荷载的作用,如地震、爆炸等。
我们通过模拟动态荷载作用下的混凝土响应,研究了混凝土的动态力学性能和损伤行为。
这有助于评估混凝土结构在极端荷载下的安全性和稳定性。
四、结果与讨论通过三维细观模型的研究,我们得到了以下结论:1. 骨料与砂浆的相互作用对混凝土的力学性能具有重要影响。
基于离散元法的混凝土三轴压缩破坏试验数值模拟耿纪莹;卜建清;高焱【摘要】混凝土是一种非均质的多相复合材料,其破坏过程具有非常复杂的应力应变关系.利用离散元软件PFC3D,模拟混凝土试件三轴压缩破坏试验,分析了摩擦系数、孔隙率、平行粘结半径、平行粘结刚度、平行粘结强度等细观参数对混凝土试件峰值强度的影响,并对混凝土试件的破坏形式进行分析.【期刊名称】《石家庄铁道大学学报(自然科学版)》【年(卷),期】2015(028)002【总页数】4页(P88-91)【关键词】混凝土;离散元;细观参数;破坏形式【作者】耿纪莹;卜建清;高焱【作者单位】石家庄铁道大学土木工程学院,河北石家庄050043;石家庄铁道大学土木工程学院,河北石家庄050043;西南交通大学土木分院,四川成都610031【正文语种】中文【中图分类】U4440 引言混凝土因其良好的性能而被广泛应用于土木工程领域,它是一种非均质多相复合材料,破坏过程具有非常复杂的应力应变关系。
为了研究混凝土试件破坏的全过程,许多学者通常将混凝土视为由粗骨料、砂浆基质以及两者之间的粘结面组成的三相复合材料。
国内类似工作主要有马怀发等[1]利用有限元法进行了混凝土试件单轴受压细观结构数值模拟;周维恒等[2]研究了混凝土试件微裂缝起裂、扩展、相互贯通的全过程,并由此提出了混凝土材料的断裂准则;肖辉等[3]进行了混凝土单轴压缩破坏试验的二维离散元数值模拟。
实际上,混凝土骨料分布和静荷载作用下的裂纹发展都是三维的,若将混凝土试件破坏过程看成是二维平面问题,则无法真实地模拟混凝土试件复杂的应力应变关系。
本文采用颗粒流软件PFC3D,从细观层面分析了混凝土试件三轴压缩破坏试验中微裂纹产生—贯通—完全破坏的全过程。
1 颗粒流方法的基本思想颗粒流程序(Particle Flow Code)基于离散单元法。
离散单元法最早由Cundall [4]在1971年提出,刚开始适用于岩石力学,1979年Cundall与Strack又联合提出了适于土力学的离散单元法,并推出二维圆盘BALL程序和三维圆球TRUBAL程序,从而形成较系统的模型与方法。
《基于三维细观模型的混凝土损伤力学行为研究》篇一一、引言混凝土作为建筑结构的主要材料,其力学性能的研究对于保障建筑安全具有重要意义。
混凝土损伤力学行为的研究是该领域的重要方向之一,而基于三维细观模型的混凝土损伤力学行为研究更是近年来研究的热点。
本文旨在通过对混凝土三维细观模型的构建与分析,探究混凝土损伤的力学行为,为混凝土结构的设计、施工和维护提供理论支持。
二、三维细观模型的构建1. 模型选取与参数设定本文采用数字图像处理技术,通过对混凝土试样的微观结构进行扫描和图像处理,构建了混凝土的三维细观模型。
在模型构建过程中,需要考虑混凝土的骨料类型、骨料粒径分布、骨料含量、水泥浆与骨料的比例等参数。
2. 模型构建方法采用离散元法构建混凝土三维细观模型。
该方法能够有效地模拟混凝土内部骨料和水泥浆的分布和相互作用,从而更真实地反映混凝土的力学性能。
在模型构建过程中,需要对混凝土试样的尺寸、骨料形状和分布等进行精确控制,以保证模型的准确性。
三、混凝土损伤力学行为的模拟与分析1. 损伤模型的建立基于三维细观模型,建立混凝土损伤模型。
该模型能够考虑混凝土在受力过程中的裂纹扩展、骨料与水泥浆的相互作用等因素,从而更真实地反映混凝土的损伤过程。
2. 模拟与分析通过有限元分析软件,对混凝土损伤模型进行模拟和分析。
在模拟过程中,需要考虑混凝土的应力-应变关系、裂纹扩展规律等因素。
通过对模拟结果的分析,可以了解混凝土在受力过程中的损伤过程和损伤程度,为混凝土结构的设计、施工和维护提供理论支持。
四、结果与讨论1. 结果展示通过模拟和分析,得到混凝土在不同应力下的损伤情况。
结果显示,混凝土在受力过程中会出现裂纹扩展、骨料与水泥浆的相互作用等现象,这些现象会导致混凝土的损伤程度加剧。
此外,还发现混凝土的损伤程度与骨料类型、骨料粒径分布、骨料含量等因素有关。
2. 讨论本文研究的结论为:混凝土在受力过程中会出现明显的损伤现象,且损伤程度与骨料类型、骨料粒径分布、骨料含量等因素有关。
收稿日期:2004Ο07Ο19基金项目:国家自然科学基金资助项目(59739180);江苏省建筑工程管理局资助项目(JS 200348);河海大学学科(院士)基金资助项目(02Ο01)作者简介:任朝军(1980—),男,河南许昌人,硕士,主要从事混凝土的细观数值模拟研究.三级配混凝土单轴破坏的细观数值模拟任朝军1,杜成斌1,戴春霞2(1.河海大学土木工程学院,江苏南京 210098;2.江苏省鸿源招标代理有限公司,江苏南京 210017)摘要:将混凝土看成是由骨料、界面和砂浆组成的非均质复合材料,采用非线性有限元软件Marc 对东江拱坝三级配混凝土试件单轴受拉、单轴受压破坏的全过程进行数值模拟.按照试件的实际配比计算出各种粒径骨料的数目,采用蒙特卡罗法对混凝土骨料进行随机投放,生成试件的随机骨料模型.研究中,单轴受拉时采用最大拉应力准则,并考虑了材料的软化;单轴受压时采用线性莫尔库伦准则.结果表明,混凝土的力学性能与其内部的细观结构组成密切相关,材料的宏观破坏是细观结构破坏累积的结果,无论是受拉时的开裂还是受压时的屈服,几乎都首先发生在界面上.关键词:细观力学;单轴破坏;三级配混凝土;随机骨料模型中图分类号:T V431 文献标识码:A 文章编号:1000Ο1980(2005)02Ο0177Ο04混凝土是一种多相非均质复合材料,其力学性能受材料细观结构的影响较大.为了得到混凝土的各种力学参数,通常要做大量的试验,耗费大量的人力物力,而所得的结果却有一定的局限性.以细观力学为基础的数值模拟方法可以方便准确地研究混凝土的各种力学性能,许多学者在该领域进行了卓有成效的研究,提出了多种数学模型.G urs on [1]于1977年提出了G urs on 模型,用以描述微孔洞损伤对材料塑性变形的影响.T vergaard [2]在G urs on 工作基础上,提出了改进的G urs on 模型,使其更具有一般性,并且考虑了微孔洞之间相互作用的效应.Bazant [3]等提出了随机粒子模型,主要用于模拟颗粒固体材料.本文以计算细观力学为基础,从细观层次上研究混凝土力学破坏的机理,为大体积混凝土静、动力学性能研究提供可能的新途径.图1 非线性本构模型Fig.1 N onlinear constitutive model 1 数值模拟的本构模型单轴受拉试件采用的非线性本构模型见图1,该模型能够很好地模拟材料的脆性破坏[4].模型采用最大拉应力破坏准则,并且考虑了材料的软化现象,当材料中的最大主拉应力超过给定的值时,在垂直于最大主拉应力的方向产生裂纹.采用软化模量对裂纹出现后材料的软化现象进行描述.对于轴压试件,采用线性莫尔库仑本构关系.2 随机骨料模型的生成和网格剖分本文采用随机骨料模型对东江拱坝三级配混凝土试件的单轴破坏进行数值模拟[5].轴拉试件采用450mm ×450mm ×3600mm 的等截面试件,其应力均匀段为450mm ×450mm ×1400mm.平面等效试件尺寸为450mm ×1400mm ,采用Walraven 公式将满足三级配的累积分布概率转化为二维的D <D 0概率,计算出3种代表粒径(60mm ,30mm 和12mm )的骨料个数分别为31,106和825.受压试件采用 450mm ×450mm 的圆柱体,其平面等效尺寸为39818mm ×450mm ,计算出3种代表粒径骨料数目分别为9,30和235,采用蒙特卡罗方法生成随机骨料模型,见图2(a ),2(b ).网格剖分之前,预先将颗粒的半径放大0105倍,将其作为界面层的外边界,内边界为将原半径内缩0105倍.网格剖分分骨料、界面和砂浆3个部分分别进行,然后对3种不同的材料参数分别赋值.计算采用的力学第33卷第2期2005年3月河海大学学报(自然科学版)Journal of H ohai University (Natural Sciences )V ol.33N o.2Mar.2005图2 试件随机骨料模型Fig.2 R andom aggregate model of specimen表1 计算采用的力学参数T able 1 Mech anical p arameters adopted in calculation 类型弹性模量E/G Pa 泊松比软化模量E s /G Pa 抗拉强度σt /MPa内摩擦角φ/(°)黏聚力c/MPa 骨料551501161111610421523153界面2510011651021030104180砂浆2610012251221530105188参数见表1[6~8].本文采用Advancing F ont网格自动剖分方法,单元为四结点等参元.剖分后轴拉试件的单元总数为92260,节点总数为92614,轴压试件的单元总数为26228,节点数目为26428.二者采用的单元边长约为3mm ,界面单元的厚度为颗粒半径的011倍.3 单轴受拉试件的计算结果分析对于轴拉试件,用弧长法计算时,初始弧长取011,最大弧长取013,加载方向为水平方向,得到试件的极限承载力为29316kN ,极限抗拉强度为1145MP a ,与试验得到的极限抗拉强度116MP a 误差为913%,计算的总增量步数为图3 单轴受拉试件最终裂缝Fig.3 Final cracks of specimen subjected to uniaxial tension22步.由计算结果可知:在第16个增量步,荷载为26513kN ,达到了极限荷载的90%,A 点附近(界面,坐标为(18019,3115),坐标原点在左下角)和B 点附近(界面,坐标为(124014,22316))几个界面单元的应力达到了极限拉应力(A ,B位置见图2(a )),最先出现微裂纹,而试件其他部分的骨料和砂浆仍然处于弹性阶段.随着荷载的继续增加,A ,B 之间一些界面单元的应力逐渐达到极限拉应力,试件在这些地方出现微裂纹.荷载继续增加,这些裂纹不断扩张、合并,第22增量步时在A 点附近形成一条贯穿试件的大裂缝,试件破坏,见图3.图5 A 点的应力应变关系Fig.5 Stress 2strain curve at point A图4 荷载位移曲线Fig.4 Load 2displacement curve 荷载位移曲线见图4,曲线明显出现下降段,并且在初始破坏点开裂很长时间后,试件的荷载位移曲线才出现弯曲,即试件宏观裂缝的出现要明显滞后于微裂缝,由此可以得出试件的宏观破坏是细观破坏累积的结果,细观层次的破坏和损伤是宏观破坏的内因和基础的结论.图5为A 点的应力应变关系曲线,其有关数值与计算采用的参数基本一致,说明计算的精度和结果是可靠的.图6(a )给出了纵坐标为170m m 的水平面上(见图2(a )左边局部放大图)各种材料开裂前后(分别为第17步和第18步)水平应力比较,其中1,2为砂浆节点,3,4为界面节点,5,6为骨料节点.由图6(a )可以看出,开裂后界面单元的应力急剧下降到零.由于界面的开裂,相邻骨料和砂浆单元的应力得到较大释放,均比开裂前有较大的减小,其中3,4点附近为最终破坏裂缝的位置.图6(b )为横坐标为180mm 的竖直面上(见图2(a )左边的局部放大图)不同材料开裂前后(分别为第17步和第18步)的应力比较,其中1,2为砂浆节点,3,4为界面节点,5,6为骨料节点.界面开裂后,立即退出工作,原本由界面和与其相邻的砂浆、骨料共同承受的荷载全部由砂浆和骨料来承受,表现在图6(b )上就是界面节点的应力急剧下降,而相邻砂浆和骨料节点的应力明显增大.通过砂浆和骨料的传递,该竖直面上另外一些界面单元逐渐开裂.871河海大学学报(自然科学版)第33卷图6 不同平面上各点开裂前后的水平应力Fig.6 H orizontal stress at eachpoint in different planes before and after cracking图7(a )给出了界面强度分别为118MPa ,210MPa 和212MPa 时B 点(界面)的应力应变关系,对应的破坏荷载分别为16713kN ,29316kN 和31319kN.由图7(a )可以看出,试件的极限承载力随着界面强度的增大而增大,且界面强度的大小只影响试件的极限承载力,对其他力学参数影响很小.图7(b )给出了界面软化模量不同时B 点的应力应变关系比较.随着软化模量的增大,下降段逐渐变陡,但下降段对软化模量的变化不是太敏感.图7 不同界面参数时B 点的应力应变关系Fig.7 Stress 2strain curves for point B with different interface p arameters图8 荷载位移曲线Fig.8 Load 2displacement curve图9 界面黏聚力c 不同时C 点的应力应变关系Fig.9 Stress 2strain curvesof different cohesions4 单轴受压试件的计算结果分析弧长法计算时初始弧长取011,最大弧长取013,加载方向为竖直方向,得到试件的极限承载力为416513kN ,抗压强度为23121MPa ,与试验得到抗压强度2111MPa 的误差为813%,计算总增量步为87步.图8为计算得到的荷载位移曲线.由于采用了线性莫尔库伦准则,认为足够多的单元产生屈服结构即破坏,故曲线未见有下降段.第37增量步的荷载为极限荷载的61%,试件中C 点(32414,18913)(位置见图2(b ))处几个界面单元的等效应力达到了屈服应力,材料开始屈服,荷载位移曲线出现非线性.随着荷载的增大,其他骨料处的界面单元逐渐出现屈服,大约在第60增量步,几乎每个骨料处都有界面单元屈服,此时,试件的屈服开始由逐渐增多变为沿着界面不断扩展、贯通,最终在第87增量步导致试件破坏.同时,随着荷载的加大,屈服应变的变化率也在逐渐增大.图9是界面黏聚力c 分别为414MPa ,418MPa 和514MPa 时C 点的应力应变关系,对应的极限承载力分别为396719kN ,416513kN 和433014kN.结果表明,随着界面黏聚971第2期任朝军,等 三级配混凝土单轴破坏的细观数值模拟081河海大学学报(自然科学版)第33卷力的提高,试件的承载能力逐渐提高,而且界面的黏聚力对应力应变关系曲线的塑性段有一定的影响,黏聚力越大,屈服段的长度越长.5 结 语a.计算结果与试验结果比较一致,表明采用计算细观力学方法研究混凝土结构的破坏过程是可行的.对于极限荷载的求解问题,弧长法是非常有效的计算方法.b.混凝土宏观破坏明显滞后于细观破坏,也就是说,宏观破坏是细观破坏逐渐积累的结果,而细观破坏则是宏观破坏的内因.c.界面是混凝土的薄弱部分.随着界面强度的提高,混凝土的极限承载能力逐渐提高,而随着界面软化模量的增大,下降段逐渐变陡,但下降段对软化模量的变化不是太敏感.影响混凝土结构强度的因素还有很多,如材料的初始损伤、颗粒的任意形状,以及对混凝土强度有较大影响界面参数的选取,这些都有待于以后进一步深入研究.致谢:本文研究过程中得到了中国水利科学研究院陈厚群院士的关心和指导,在此谨向陈院士表示感谢!参考文献:[1]G URS ON A L.C ontinuum theory of ductile rupture by v oid nucleation and growth,I.Y ield criteria and flow rules for porous ductilemedia[J].J Eng Mater T ech,1977,99(1):2—15.[2]T VERG AARD V.In fluence of v oids on shear band instabilities under plane strain condition[J].International Journal of Fracture,1981,17(4):389—407.[3]BAZ ANT Z P,OH B H.Microplane m odel for progressive fracture of concrete and rock[J].Journal of Engineering Mechanics,ASCE,1985,111(4):559—582.[4]宋玉普.多种混凝土材料的本构关系和破坏准则[M].北京:中国水利水电出版社,2002.140—141.[5]刘文彦,叶文瑛,葛辉,等.东江拱坝全级配混凝土力学性能的试验研究[J].水力发电,1986,(5):8—14.[6]尚岩.大体积混凝土静、动力学性能的数值模拟[D].南京:河海大学,2004.[7]刘斌.碾压混凝土细观损伤断裂的仿真模拟[D].北京:北京工业大学,2000.[8]唐春安,傅宇方,朱万成.界面性质对颗粒增强复合材料破坏模式影响的数值模拟分析[J].复合材料学报,1999,(10):113—120.Meso2structure numerical simulation of uniaxial failureof three2graded concreteREN Chao2jun1,DU Cheng2bin1,DAI Chun2xia2(1.College o f Civil Engineering,Hohai Univ.,Nanjing210098,China;2.Hongyuan IF B Surrogate Limited Company o f Jiangsu Province,Nanjing210017,China)Abstract:C oncrete was taken as inhom ogeneous com posites consisting of m ortar,aggregate and their interface. Numerical simulation was carried out by use of nonlinear finite element s oftware Marc,which could simulate the whole failure process of specimens of three2graded concrete subjected to uniaxial com pression and tension for the Dongjiang arch dam.The M onte Carlo method was used to generate the random aggregate m odel,and the number of aggregate of different sizes was determined according to the actual mix proportion.The maximum tensile stress criterion and the linear M ohr2 C oulomb criterion were adopted in the study of specimens subjected to uniaxial tension and uniaxial com pression, respectively,with the s oftening of materials taken into account.The results indicate that the mechanical performance of concrete is closely related to its mes o2structure,and the mes o2failure is the result of accumulation of damages to its mes o2 structure.The results als o show that the interface is the weakest part,where tension2induced cracks or com pression2 induced yield are m ost likely to occur.K ey w ords:mes o2mechanics;uniaxial failure;three2graded concrete;random aggregate m odel。