混凝土拉伸时细观数值分析.docVIP

混凝土拉伸时细观数值分析摘要：从细观力学角度,把混凝土看作是由骨料、水泥砂浆及二者之间的界面组成的复合材料，应用TrueGrid软件形成了高质量的分块映射网格，建立了混凝土立方块试件的三维数值模型。模拟了混凝土试件在单向拉伸时裂纹的产生、发展过程。 关键词：混凝土 ;随机骨料模型; 映射网格;细观数值模拟 中图分类号：TU528.01 文献标识码：B 文章编号：1008-0422（2008）08-0163-03 1前言 混凝土作为重要的建筑材料已有百余年历史，人们从不同的结构层次来研究混凝土材料的力学性能。（1）宏观层次。一般认为特征体积相当于3～4 倍的最大骨料体积。当小于特征体积时，当研究对象大于特征体积时，材料假定为均质。（2）细观层次。从分子尺度到宏观尺度之间，其结构单元尺度变化范围在10-4cm至几cm。在这个层次上，混凝土被认为是一种由骨料、水泥砂浆及二者之间的界面组成的三相复合材料[1]。(3) 微观层次。材料的结构单元尺度在原子、分子量级，即从小于10-7cm～10-4cm 。材料由晶体结构及分子结构组成。随着高性能计算机及数值算法的迅速发展，用数值试验来模拟混凝土材料破坏过程的研究已成为热点，相对廉价的细观数值分析有望弥补或替代部分材料试验[2] 2骨料尺寸 选取300mm×300mm×300mm的立方体试件，并将它们取为二级配混凝土，其中大石(直径40mm~80mm)、中石(直径20mm~40mm)进行数值模拟。按富勒理想级配曲线： （1） 其中D为筛孔直径；Dmax为骨料最大粒径；P为能通过筛孔D的骨料的质量百分比。 取Dmax为80mm，求得大石、中石的质量百分比： 由表1可知，粒径为40mm~80mm的骨料质量百分比为100%-70.7%=29.3%，粒径为20mm~40mm的骨料质量百分比为70.7%-50%=20.7%。骨料容重取2800Kg/m3，试件的石料用量为1227Kg/m3。适当选择球体骨料的个数与直径，使之满足级配关系并容易随机填充，本文取大骨料为大石，个数为27,直径为62.58mm。小骨料为中石，个数为189,直径为29.16mm。 3三维有限元网格的划分 随机骨料模型的网格剖分，长期以来是三维细观力学分析的难题之一。分网技术简单地说分为三种：六面体映射或拉伸网分、四面体全自动网分、混合网分。（1）商用软件中，TRUEGRID、ANSYS ICEM CFD、AVL_FAME等可以完成六面体映射网格划分。划分随机骨料模型时，主要思想是：先生成均匀的有限元网格覆盖整个分析区域，根据蒙特卡罗方法和骨料从大到小的投放原则，随机生成一套虚拟的骨料域。根据单元节点是否全部落在、部分落在和不落在骨料域，判别单元是否属于骨料单元、界面单元和砂浆单元。这种映射网格技术能够解决复杂形状骨料的三维模拟，全部是六面体网格，有限元的单元精度高，但界面区域不够光滑。在计算上，考虑到界面的厚度很小，要采用较小的单元而在全域上采用一样大小的单元，导致单元个数太多，计算量较大，在微机上计算难度大。（２）多数的CAE软件都提供四面体网格全自动划分功能，配合局部加密技术划分单元，可以建立随机骨料网格模型，计算量可以接受，但四面体单元的精度较低，空间填充效率不高，且视觉效果稍差。（３）CUBIT、MARC、HARPOON等软件提供混合分网功能，产生的网格主要为六面体，边界复杂的区域用五面体、棱柱体或四面体填充。分网效率高，但边界区域通常有应力集中，需要特别关注的地方却用的是低精度网格。 本文提出先分块映射网分再随机装配的方案，可以兼顾单元精度、局部加密和计算量的要求。首先，将混凝土试件看成由三种块体组合而成：（1）包含大骨料、界面和砂浆的块体，简称大骨料块。如图１所示。（2）包含小骨料、界面和砂浆的块体，简称小骨料块。如图２所示。（3）只包含砂浆的块体，简称砂浆块。如图３所示。分别对这三种块体用映射网分。 以上是球体骨料的情形，类似地，可作出椭球体骨料的映射网格。如图４所示。 注意到界面的材料刚度与骨料的刚度有较大的差别，为了减少应力集中，将骨料的外围网格作了适当的加密。 三种块体都是正方体，使大骨料块的边长为66.7mm，小骨料块和砂浆块边长33.3mm，界面单元厚度为1mm，骨料和界面单元全部位于块体的内部。大骨料块的边长为小骨料块和砂浆块边长的两倍，表面上的单元尺寸相同，显然这三种块体可以无缝、无重叠地填充三维空间，而且将重合的节点合并后，各块之间形成共节点连接的整体有限元网格，如图5所示。随机装配过程利用了MARC软件的PYTHON语言前处理接口，调用随机函数以及单元的复制、移动功能。整体有限元共有252