碾压混凝土细观结构力学性能的数值模拟.doc

《碾压混凝土细观结构力学性能的数值模拟.doc》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、碾压混凝土细观结构力学性能的数值模拟摘要：本文在细观层次上将碾压混凝土看成是由硬化水泥粉煤灰砂浆体、粗骨料颗粒及二者间的粘结带构成的两相非均质复合材料，以随机骨料模型代表碾压混凝土的细观结构，采用有限元技术模拟了劈拉试件的开裂过程，所得碾压混凝土的劈拉强度与试验结果基本吻合。为碾压混凝土破坏机理的研究、碾压混凝土力学性能的数值模拟提供了新的技术途径。关键词：细观复合材料随机骨料模型有限元破坏机理强度　　碾压混凝土是一种可用土石坝施工机械设备运输及铺筑，用振动碾压实的干硬性混凝土，我国从1978年开始对碾压混凝土筑坝技术进行研究，在大量试验研究

2、的基础上，将其广泛应用于实际工程，于1986年建成我国第一座高57m的坑口碾压混凝土重力坝。目前，我国已建、在建碾压混凝土坝约43座[１]。碾压混凝土材料性能的研究是碾压混凝土筑坝的关键技术课题之一，碾压混凝土材料力学行为的研究，大都建立在试验的研究基础上，需要花费大量的人力、物力，所得试验成果还受试验条件的限制。本文拟采用随机骨料模型[２]及数值模拟技术，在细观层次上探索碾压混凝土细观结构与宏观力学性能的关系。1随机骨料模型的形成　　本文假定已知级配的骨料其颗粒形状为球形，按骨料达到最优密实度条件的Fuller三维级配曲线，其级配曲线表达式

3、为[３]：Y＝(D0/Dmax)1/2(1)式中：Y为骨料通过直径为D０筛孔的重量百分比；Dmax为最大骨料颗粒径。WalravenJ.C.基于Fuller公式，将三维级配曲线转化为试件内截平面上任一点具有骨料直径D<D0的概率PC(D<D0)[４]，其表达式为：PC(D<D0)=Pk[1.065(D0/Dmax)0.5－0.053(D0/Dmax)４－0.012(D0/Dmax)69(2)－0.0045(D0/Dmax)8＋0.0025(D0/Dmax)10式中：Pk为骨料体积占试件总体积的百分比，一般取

4、Pk＝0.75；Dmax为骨料最大粒径。　　根据不同的D０值，由式(2)可求得概率分布曲线PC(D<D0)～D0/Dmax，据此可求得在试件内截平面上各种料径的颗粒数。　　由以上结果，将小于5mm的细骨料计入砂浆均质体，按各种粗骨料在截面上不相重叠的条件，由蒙特卡洛法，随机生成各粗骨料颗粒的形心坐标，即可形成随机骨料模型。　　在本文中，以桃林口水库碾压混凝土劈拉试验为算例进行数值模拟，采用该水库碾压混凝土芯样尺寸，模拟试件的随机骨料模型及尺寸如图1所示。2有限元程序、计算网格及材料参数图1碾压混凝土劈拉试件骨料随机分布2.1有限元程序　

5、本文利用我院自编的结构分析有限元程序NLFE1，以碾压混凝土劈拉试件仿真模拟分析为例，进行细观结构力学性能研究。　　　　对混凝土拌和物的两个主要构元粗骨料和硬化水泥粉煤灰砂浆分别进行试验，都会出现弹性破坏，而碾压混凝土本身在破坏前却表现出非弹性性质[５]。因此本文采用碾压混凝土的各相在破坏前符合弹脆性材料假设，应力-应变关系为线弹性关系。当单元体主应力大于该相材料强度时，则单元发生开裂，程序自动用简单超余应力转移法将超余应力转移至相邻未破坏单元。2.2计算网格　　在随机地生成骨料颗粒系统之后，对结构进行有限元分析时，首先要将研究区域划分为有限

6、元网格，这将涉及到成千上万个单元和节点，数据准备工作则变得极其繁重。特别是对于需要随机改变颗粒组合，反复进行网格剖分的情况下，不寻求某种自动或半自动的方法，是难以实现的。(3454个单元，1833个结点)图2碾压混凝土劈拉试件的有限元计算网格9　　本文采用Delaunay三角剖分的原理[６，７]，通过引进程序实现了对随机骨料颗粒分布区域的全自动剖分。当某三角形单元的3个结点均落在骨料颗粒内时，定为骨料单元；当某三角形单元的3个结点均落在硬化水泥粉煤灰砂浆区域内时，定为砂浆单元；当某三角形单元的3个结点分别落在骨料和砂浆区域内时，则定为粘结带单

7、元。对不同单元分配不同的材料特性。为模拟碾压混凝土层面处理的条件，在试件中部设置厚度为1cm的水泥砂浆层面，考虑施工因素影响，可使层面力学性能弱化，故取水泥砂浆层面力学性能与粉煤灰砂浆相同。以上材料类型的判断完全通过编程由计算机实现。对于劈拉试件的有限元计算网格如图2所示。2.3材料参数随机骨料模型材料性能，综合文献［5，7，8，9］及有关试验资料后，取值见表1.表1随机骨料模型材料性能材料抗压强度fc/MPa抗拉强度ft/MPa弹性模量E/GPa泊松比μ水泥砂浆层面25.092.526.00.22水泥粉煤灰砂浆25.02.526.00.2

8、2骨料80.010.035.50.16粘结带22.01.525.00.1693数值模拟结果归纳及分析3.1碾压混凝土试件受劈拉的开裂过程　为了模拟开裂过程，取单宽厚