基于柔性边界和Hill接触模型的PFC3D非饱和土动力特性分析.pdfVIP

摘 要 摘 要 在自然界中非饱和土可能受到各种荷载作用，特别是受到动力荷载的作用， 针对非饱和土的动力特性研究尤为重要。以往的研究中大多采用理论模型进行 非饱和土的动力特性分析，但理论模型在复杂荷载情况下特别是在动真三轴试 验的验证中，由于仪器较为昂贵且操作较为复杂，很多实验室不具备验证理论 模型的试验条件，所以复杂荷载动真三轴情况下难以判断理论模型的正确性。 采用PFC3D 进行数值模拟不仅可以从细观层面分析试样加载过程中的微观变化， 同时PFC3D 可以对非饱和动真三轴试验及大型工程问题进行数值模拟分析，进 而验证理论模型的正确性。本文针对江西非饱和红土，通过大量土工试验得出 理论模型计算及PFC3D 细观参数标定所需必要参数，采用GDS 非饱和三轴仪 与温控饱和动三轴仪进行室内静、动荷载三轴试验，并在室内进行了隧道振动 对周围土体动力特性影响的模型试验。在进行PFC3D 非饱和土数值模拟时，以 往的研究采用等量代换或二次开发等方法模拟非饱和土试验，但等量代换法存 在较大误差且与实际非饱和土试验不完全相符，二次开发对编程能力要求较高， 实现起来较为困难，难以推广。本文通过调用可以直接给颗粒间赋予基质吸力 的Hill 接触模型对非饱和土三轴试验进行数值模拟，更加方便、准确得实现了 对非饱和土的数值试验研究。通过将非饱和土动力特性理论模型与PFC3D 中隐 藏算例包Hill 接触模型结合建立了柔性边界的常规静、动三轴与静、动真三轴 数值模型，并以室内模型试验箱为基础建立了1:1数值模型，数值模拟结果表明 Hill 接触模型在PFC3D 中可以很好的模拟非饱和土材料，为非饱和土的研究提 供了新方法。研究内容及结论如下： （1）在PFC3D 中存在很多类似Hill 接触模型的自定义接触模型算例包， Hill 接触模型以液桥理论为基础，通过程序自动判断颗粒间隙是否小于等于2 倍 的最小接触颗粒半径，若符合条件则直接给颗粒间赋予基质吸力，通过将Hill 接触模型与柔性边界相结合建立了与室内试验相同的数值模型，结果显示数值 模拟与室内非饱和土试验数据在合理误差范围内，总体上吻合程度较高，为后 续研究提供了理论支持。 （2）通过室内基本土工试验、压缩回弹试验、抗剪强度指标试验得到理论 模型计算及PFC3D 细观参数标定所需必要参数。数值模型以室内试验结果为细 I 摘 要 观参数标定基础，根据模型特点以杨氏模量、泊松比、摩擦系数和颗粒局部阻 尼为主要细观参数进行标定，将所得细观参数带入数值模型中建立非饱和常规 静三轴数值模型并与室内试验结果对比，结果显示二者应力应变曲线及体应变- - 轴向应变曲线误差在5%以内，整体上较为吻合，说明Hill 接触模型可以用来模 拟非饱和土试验。 3 PFC3D （）以室内试验为原型并结合理论模型采用 建立了不同压实度、 不同基质吸力、不同围压的常规静三轴试验，将数值模拟结果与室内试验和理 5% 论计算结果进行对比，发现曲线结果误差在 以内，基质吸力和围压越大，试 样的偏应力与体应变越大，压实度越大试样的偏应力越大、体应变越小。数值 模拟中试样的位移图及力链图为非饱和土的微观层面研究提供了可靠支持，分 析发现围压与基质吸力越大试样抵抗变形的能力越强。根据室内试验和课题组 成员所提三剪次加载面理论建立了不同基质吸力、不同围压、不同振动幅值的 常规动三轴试验，将数值模拟与试验数据进行对比发现，二者应力-应变及轴向 应变 循环次数曲线存在较小误差，误差在合理范围内，分析结果发现围压和循- 环荷载幅值增大，滞回圈面积变小，动剪切模量增大，基质吸力增大累计塑性 变形和滞回圈面积均减小；在进行动三轴数值模拟中配位数的变化反应试样的 宏观力学行为，增大围压和基质吸力配位数变大，加载过程中配位数的下降趋 势与变化幅度均减小，循环荷载幅值增大，配位数下降速度与振动幅度增大。