混凝土动态力学性能蒙特卡洛法模拟计算研究.pdfVIP

? 建筑材料? 混凝土动态力学性能蒙特卡洛 法模拟计算研究 孙智勇 中铁第五勘察设计院集团有限公司 北京 摘 要 采用蒙特卡洛方法建立随机骨料模型,应用断裂力学原理建立材料破损有限单元刚度矩阵,编制非线形 动态力学方程求解程序,在细观结构层次上进行混凝土动态力学性能模拟计算,得到了混凝土试块中裂缝产生、开 展,直至混凝土试块破坏的全部形态和过程,能较好地与实验结果相吻合,证明了本方法可以成为有效研究混凝土 材料动态力学性能的新途径、新手段。 关键词 混凝土 动态力学性能 断裂力学 模拟计算 中图分类号 . 文献标识码文章编号?? ?, ? , . ,,. ,,,. ,. ;; ; 言 则 模拟混凝土的断裂过程也有了一定进展。 年理论假设 . 提出了网络结构 混凝土是由水泥 、砂子、石子和水按一定比例 模型,把骨料、固化砂浆及胶结面特性随机地分配 配合比 拌和 ,并经过水化形成的人工石材。水泥 在网格模型中,并模拟了混凝土及其他非均质材料 和水在凝结硬化过程中形成水泥胶块 ,将砂子和石 所表现的典型破坏机理和开裂面的贯通过程。 子等骨料粘结在一起 。 砂子和石子等骨料和水泥 年刘光廷、王宗敏用随机骨料模型数值模拟了混凝 胶块中的结晶体组成弹性骨架承受外力 。 弹性骨 簟:。 手毒庭罄 模拟 架使混凝土具有弹性变形的特点,而水泥胶块中未 度及劈裂破坏过程 胶 至真 荏 击嘉 。 , 并引入了断裂力学的强度准 ‘享? 建筑材料. 混凝土结构层次 力,而是先局部开裂,然后随着荷载的增加裂缝逐 微观结构。微观结构指原子、分子量级,其 渐贯通整个单元。此过程中,。当到达某一时 单元尺度~~ ~ ,着眼于水泥水化物的 刻,该单元在垂直 微观结构分析。 裂缝方向上的抗 细观结构。从分子尺度到宏观尺度拉能力全部失去,一一 ,着眼于粗骨料、水泥水化物、孔隙、界面 此时/.,如图 所示。 等细观结构,组成多相复合材料,可按损伤力学或 随机骨料模型进行破坏机理分析。在这个层次上, 图 ,占 // , 为开裂单元 图 垂直裂缝方向刚度变化 混凝土被认为是一种由骨料、硬化水泥砂浆和它们 之间的过渡区 粘结带 组成的三相材料。由于泌 的抗拉强度, 为开裂单元的弹性模量。 水、干缩和温度变化引起骨料和水泥浆体之间产生 为单元完全开裂时的应变。实际上,由于外 粘结裂缝,细观内部裂隙的发展直接影响混凝土的 加荷载的继续增大,开裂单元的应变也会继续增 宏观力学性能。细观层次模型一般在毫米到几个 加,此时的应变可能会大于占 为了能够模拟出破 厘米量级。本文就是在这个层次上进行研究。 坏单元在②阶段时的刚度和内力值,本程序作了如 宏观结构。混凝土实体结构,如梁、板、柱、 下处理: 混凝土坝体等, ,可按连续介质力学分析 线 ①在单元达局部开裂后,弹性矩阵 ? 式中 弹性或非线形有限元、边界元 。 的。 . 与混凝土破坏有关的细观内部结构 ②当荷载继续增加时,利用子程序判断该单元 与混凝土破坏有关的细观内部结构为:硬化水 应变是否达到 ;如果超过占 令/并撤消该 泥胶体、骨料、裂缝、水泥浆体一骨料界面。 单元的反向作用力。此时该单元已经完全开裂且 . 逐渐开裂模型 在垂直裂缝方向上没有任何力的“贡献”。 破坏单元的局部坐标下弹性矩阵: 此时计算所得静荷载抗压强度与试验值接近。 本模型是在不知道细观裂缝一次性开裂尺寸的前提下所建立的细观模型;并且比较这两个模型是为【 了讨论人为的有限元的网格划分 包括有限单元的 绝对尺寸、相对于骨料粒径的相对尺寸、相对于细 式中, 为破坏单元垂直裂缝方向弹性模量保 观裂缝的尺寸等 对细观结构的计算的影响。 留系数; ≥ ≥ ,为剪切模量保留系数,用以反映 混凝土所有组成材料都是线弹性、脆断的,本 骨料的咬合作用;忽略泊松比的影响,即认为垂直 文统一将混凝土的各种材料设定为线弹性一脆性断 裂缝方向与平行于裂缝的方向没有相互作用。 裂的方式。 采用右手螺旋定则,以垂直裂缝方向为局部坐 标的 方向。 蒙特卡洛方法 整体坐标下的弹性矩阵为: 。按照蒙特卡洛方法 在试 件内随机生成各种骨料颗粒的形心坐标,即可生成 随机骨料模型,其步骤如下: ?给定边界范围及坐标系。 一 一产生在 , 区间上的均匀分布随机数 伪 式中, , , 是裂缝面的法线方 随机数 。 向与整体坐标 轴的夹角。 本文采用混合同余法,其递推公式为: 在此基础上提出了逐渐开裂模型:当裂缝单元 在单元主应变达到或超过 时,其垂直裂缝方向上一一 并未立即失去其在垂直裂缝方向上的全部抗拉能 铁道建筑技术建筑材料? 时只是短短几秒,因此本程序描述了在瞬间冲击荷 式 表示以 为模数的同余式,即以 除 载作用下的动力反应。同时加入材料单元的破坏 后得到的余数