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

PAGE1 / NUMPAGES2 PAGE1 / NUMPAGES2 碾压混凝土细观结构力学性能的数值模拟 摘要：本文在细观层次上将碾压混凝土看成是由硬化水泥粉煤灰砂浆体、粗骨料颗粒及二者间的粘结带构成的两相非均质复合材料，以随机骨料模型代表碾压混凝土的细观结构，采用有限元技术模拟了劈拉试件的开裂过程，所得碾压混凝土的劈拉强度与试验结果基本吻合。为碾压混凝土破坏机理的研究、碾压混凝土力学性能的数值模拟提供了新的技术途径。 关键词：细观 复合材料 随机骨料模型 有限元 破坏机理 强度 nbsp; 　　碾压混凝土是一种可用土石坝施工机械设备运输及铺筑，用振动碾压实的干硬性混凝土，我国从1978年开始对碾压混凝土筑坝技术进行研究，在大量试验研究的基础上，将其广泛应用于实际工程，于1986年建成我国第一座高57m的坑口碾压混凝土重力坝。目前，我国已建、在建碾压混凝土坝约43座[１]。碾压混凝土材料性能的研究是碾压混凝土筑坝的关键技术课题之一，碾压混凝土材料力学行为的研究，大都建立在试验的研究基础上，需要花费大量的人力、物力，所得试验成果还受试验条件的限制。本文拟采用随机骨料模型[２]及数值模拟技术，在细观层次上探索碾压混凝土细观结构与宏观力学性能的关系。 1 随机骨料模型的形成 　　本文假定已知级配的骨料其颗粒形状为球形，按骨料达到最优密实度条件的Fuller三维级配曲线，其级配曲线表达式为[３]： Y＝(D0/Dmax)1/2 (1) 式中：Y为骨料通过直径为D０筛孔的重量百分比；Dmax为最大骨料颗粒径。Walraven J.C.基于Fuller公式，将三维级配曲线转化为试件内截平面上任一点具有骨料直径Dlt;D0的概率PC(Dlt;D0)[４]，其表达式为： PC(Dlt;D0)= Pk[1.065(D0/Dmax)0.5－0.053(D0/Dmax)４－0.012(D0/Dmax)6 (2) －0.0045(D0/Dmax)8＋0.0025(D0/Dmax)10 式中：Pk为骨料体积占试件总体积的百分比，一般取Pk＝0.75；Dmax为骨料最大粒径。 　　根据不同的D０值，由式(2)可求得概率分布曲线PC(Dlt;D0)～D0/Dmax，据此可求得在试件内截平面上各种料径的颗粒数。 　　由以上结果，将小于5mm的细骨料计入砂浆均质体，按各种粗骨料在截面上不相重叠的条件，由蒙特卡洛法，随机生成各粗骨料颗粒的形心坐标，即可形成随机骨料模型。 　　在本文中，以桃林口水库碾压混凝土劈拉试验为算例进行数值模拟，采用该水库碾压混凝土芯样尺寸，模拟试件的随机骨料模型及尺寸如图1所示。 2 有限元程序、计算网格及材料参数 图1 碾压混凝土劈拉试件骨料随机分布 2.1有限元程序 1　本文利用我院自编的结构分析有限元程序NLFE1，以碾压混凝土劈拉试件仿真模拟分析为例，进行细观结构力学性能研究。　　 　　对混凝土拌和物的两个主要构元 粗骨料和硬化水泥粉煤灰砂浆分别进行试验，都会出现弹性破坏，而碾压混凝土本身在破坏前却表现出非弹性性质[５]。因此本文采用碾压混凝土的各相在破坏前符合弹脆性材料假设，应力-应变关系为线弹性关系。当单元体主应力大于该相材料强度时，则单元发生开裂，程序自动用简单超余应力转移法将超余应力转移至相邻未破坏单元。 2.2 计算网格 　　在随机地生成骨料颗粒系统之后，对结构进行有限元分析时，首先要将研究区域划分为有限元网格，这将涉及到成千上万个单元和节点，数据准备工作则变得极其繁重。特别是对于需要随机改变颗粒组合，反复进行网格剖分的情况下，不寻求某种自动或半自动的方法，是难以实现的。 (3454个单元，1833个结点) 图2 碾压混凝土劈拉试件的 有限元计算网格 　　本文采用Delaunay三角剖分的原理[６，７]，通过引进程序实现了对随机骨料颗粒分布区域的全自动剖分。当某三角形单元的3个结点均落在骨料颗粒内时，定为骨料单元；当某三角形单元的3个结点均落在硬化水泥粉煤灰砂浆区域内时，定为砂浆单元；当某三角形单元的3个结点分别落在骨料和砂浆区域内时，则定为粘结带单元。对不同单元分配不同的材料特性。为模拟碾压混凝土层面处理的条件，在试件中部设置厚度为1cm的水泥砂浆层面，考虑施工因素影响，可使层面力学性能弱化，故取水泥砂浆层面力学性能与粉煤灰砂浆相同。以上材料类型的判断完全通过编程由计算机实现。对于劈拉试件的有限元计算网格如图2所示。 2.3 材料参数 随机骨料模型材料性能，综合文献［5，7，8，9］及有关试验资料后，取值见表1. 表1 随机骨料模型材料性能 材料 抗压强度fc/MPa 抗拉强度ft/MPa 弹性模量E/GPa 泊松比μ 水泥砂浆层面 25.0 2.5 26.0 0.