粗粒土的非共轴性及其离散元数值模拟.pdfVIP

水 利 学 报 2010年2月 SHUILI XUEBAO 第41卷 第2期 文章编号：0559-9350（2010）02-0173-09 粗粒土的非共轴性及其离散元数值模拟 1，2 1，2 1，2 黄茂松 ，孙海忠 ，钱建固 （1. 同济大学 岩土及地下工程教育部重点实验室，上海 200092； 2. 同济大学 地下建筑与工程系，上海 200092） 摘要：以粗粒土为原材料的土石坝在蓄水期或者水位突降的情况下，粗粒土的主应力轴将产生旋转，从而出现土 体单元主应力方向和主应变率方向的非共轴现象，而传统的弹塑性本构模型基本上是建立在大三轴试验结果的基 础上，并不能考虑粗粒土的非共轴性。本文将传统的弹塑性本构模型三维化，引入三维非共轴塑性流动理论，建 立了粗粒土的非共轴本构模型，该模型对粗粒土的单剪试验进行了数值模拟，预测了在单剪过程中主应力方向和 主应变率方向的变化规律。最后，运用离散单元法（二维颗粒流程序PFC2D）对粗粒土的非共轴性进行验证并进行 了细观分析。研究结果表明，本文模型的计算结果和试验结果较为一致，并能较好地描述粗粒土的非共轴现象。 关键词：粗粒土；颗粒流；单剪试验；本构模型 中图分类号：TU435 文献标识码：A 1 研究背景 近20年来是水坝建设的黄金期，其中已建成的水坝绝大多数是以粗粒土为原材料的面板土石 坝。在面板土石坝设计中备受关注的问题是如何控制混凝土面板的应力与变形，然而面板的应力与 变形在很大程度上又取决于其斜依的粗粒土的变形。目前有限元数值分析方法依然是粗粒土变形分 析的最主要手段，而采用的本构模型又是有限元计算的核心。土石坝在施工期和蓄水期其变形特征 有着明显的不同，在施工阶段，就某一固定区域内的应力变化规律是大小应力同步增加，大小主应 力方向变化很小。然而水库在蓄水期或者水位突降的情况下将导致粗粒土主应力比发生明显的变 ［1］ 化，主应力方向也将发生明显的旋转 ，这将导致主应力方向和主应变率方向的不一致，即非共轴 性。由于传统的本构模型几乎都是在室内大三轴试验基础上建立的，加载过程中荷载增量始终限制 在应变的主方向（即比例加载条件），所得的模型因此存在一个共同的局限，即无法客观地描述粗粒 土在真实荷载作用下由主应力旋转所引起的非共轴性。针对这个问题，一些学者纷纷提出了基于多 ［2］ ［3］ ［4］ ［5］ 面概念 、边界面概念 、亚塑性概念 以及中间面概念 等非共轴模型来预测主应力和主应变率的 变化趋势。然而，这些模型的建立没有基于人们熟悉的弹塑性理论，因此，在实际应用中（如有限元 ［6-7］ 计算）受到一定的限制 。 目前，针对粗粒土所开展的试验研究主要有常规大三轴试验和单剪试验，前者属于比例加载条 件，故不能反映粗粒土的非共轴性，后者在加载过程中，会发生主应力旋转，从而引起主应力和主 应变率的非共轴性。值得一提的是，由于试验条件的限制，很难跟踪监测到主应力方向和主应变率 ［8］ 方向的变化规律。此外，Airey等 汇总了前人的单剪仪研究成果发现，单剪中的应力状态不能完全 确定，只能测出水平面上的正应力和剪应力，计算不出整个应力状态来。因此，粗粒土的非共轴特 性研究无法依赖于现有的试验手段。近年来细观土力学的研究方法—离散单元法的出现为开展粗粒 收稿日期：2008-09-22 基金项目：国家自然科学基金重点项目 作者简介：黄茂松（1965-），男，浙江玉环人，教授，博士，主要从事岩土工程研究。E-mail: mshuang@tongji.edu.cn — 173 — 土非共轴特性这一难题的研究开辟了新的途径。离散单元法将土体处理成离散材料，并着眼于颗粒