考虑剪胀性和应变软化粉细砂双屈服面本构模型.pdfVIP

水 利 学 报 年 月 第 卷 第 期 文章编号 考虑剪胀性和应变软化的粉细砂双屈服面本构模型 黄茂松 扈 萍 张宏博 同济大学岩土及地下工程教育部重点实验室 上海 同济大学地下建筑与工程系 上海 山东大学土建与水利学院 山东济南 摘要 针对上海粉细砂不存在唯一临界状态线的特点 对传统的砂土本构模型进行了改进 提出了一个能合理描述 粉细砂剪胀性和应变软化特性的弹塑性本构模型 该模型采用双屈服面形式 可同时反映剪切变形及压缩变形机 理 模型对传统修正剑桥模型中的剪胀性公式进行了改进 考虑了反映塑性体积应变由正转为负时对应的特征状 态应力比与初始有效围压的相关性 为了描述应变软化特性 提出了一个利用残余状态应力比和峰值应力比的应 变软化公式 可较为合理地反映粉细砂的应变软化特性 通过对上海粉细砂的多组试验结果模拟 验证了本文模 型的合理性和有效性 关键词 粉细砂 剪胀性 应变软化 双屈服面 中图分类号 文献标识码 研究背景 土体的剪胀性一直是本构模型中研究的关键问题 对于砂土而言 剪胀性与砂土的密实度和有效 围压密切相关 自 提出应力剪胀理论以来 其作为流动法则早期曾广泛应用于描述砂土力学 特性的本构模型中 然而 由于这些模型忽略了剪胀性与材料状态之间的关系 将不同密实度的同一种 砂土作为不同的材料 因而使得同一种砂土有多组材料参数 并且不能很好地描述砂土剪切过程中的力 学响应 近年来 一些学者 如 等 蔡正银等 基于 等 提出的状态 参数的概念 针对不同类型砂土在不同初始密度及有效围压下的变形特性进行了大量试验研究 并在此 基础上成功地建立了基于材料状态相关临界状态理论的砂土模型 这些模型将不同初始密实度的砂土 视为同一种材料 针对同一种砂土采用统一的一套材料参数 当然砂土的初始密实度和有效围压也是输 入参数 尽管如此 由于砂土变形特征的复杂性 使研究者在砂土是否存在唯一临界状态线问题 上一直存在较大争议 而且有关粉细砂包括粉煤灰在内的试验结果也显示并不存在唯一的临界状态线 这种基于材料状态相关临界状态理论的砂土模型一般只对洁净的中粗砂比较适合 因此 为了合理地 描述上海粉细砂的力学特性 克服无法确定临界状态线的局限性 本文对上海粉细砂的研究仍然基于传 统的砂土本构模型框架 但传统的砂土本构模型存在着一个明显的缺陷 就是对于同一初始密实度不同 初始有效围压的砂在剪缩向剪胀转换时采用同一特征状态应力比 应力比为偏应力与平均有效应力之 比 而本文所采用的特征状态应力比将与初始有效围压有关 应变软化也是砂土的一个重要变形特性 大量三轴试验 表明 砂土在剪切强度达到峰值应 力比后往往会出现应变软化 特别是对于密砂或中密砂 因此在本构模型中应该合理地反映这种应变软 化现象 在基于材料状态相关临界状态理论的砂土模型中 比如 等 都是利用峰值应力比和当 前状态参数 的关系来隐式地模拟砂土的应变软化现象 问题是需要利用临界状态线来确定当前状态 收稿日期 基金项目 国家自然科学基金重点项目 作者简介 黄茂松 男 浙江玉环人 教授 主要从事岩土工程的研究 参数 为此 本文提出一个利用残余状态应力比与峰值应力比的应变软化公式 此外 现有的砂土本构模型多为单屈服面模型 在模拟砂土变形特别是一维压缩变形时存在明显的 缺陷 因此有必要引进双屈服面加以改进 目前已有一些适用于砂土的双屈服面模型提出较早期的工 作包括 等 以及 等 本文则采用引入形状参数 的修正剑桥模型形式作为压 缩屈服面 建立一个适用于粉细砂的双屈服面本构模型 模型描述 弹性模量 弹性剪切模量 的变化取决于当前孔隙比 与有效平均应力 这里采用 等 给出的经验表达式 即 式中 为孔隙比 为大气压力