砂土弹塑性损伤本构关系的数值建模.pdfVIP

维普资讯 第 22卷第 1期 华 中 科 技 大 学 学 报 (城市科学版) VolI22NO．1 2005年 3月 J．ofHUST． (Urban ScienceEdition) M ar．2005 砂土弹塑性一损伤本构关系的数值建模 杨 杰 王靖涛 (1．华中科技大学 土木工程与力学学院．湖北 武汉 430074) 摘 要 ：认为岩土介质变形 中出现的剪胀现象实质上是一个损伤过程．通过 中密砂的三轴压缩试验 ，发现其在 剪胀后 ，弹性模量和剪切模量大幅度下降 ，这进一步证实 了剪胀的确是个损伤演化过程．引进了一个损伤变量 描述这个过程 ，结果表明，损伤程度相对塑性体积应变的变化速率与初始固结压力相关 ，且随压力增大而减 小．同时利用数值建模方法建立了砂土的弹塑性一损伤模型 ，考虑了弹塑性变形与损伤的耦合效应 ．获得 了相 应的应力应变关系及它的三维曲面． 关键词 ：损伤 ； 应力路径； 数值建模 ； 本构关系； 剪胀 中图分类号 ：TU452 文献标识码 ：A 文章编号：1672—7037(2005)01—0063—03 岩石、超固结粘土、密砂及 中密砂在弹塑性变 试验 中，在增加轴压的同时减小围压，进行了初始 形过程 中，体积 由压缩转变为后期的膨胀，即出现 固结压力 P。分别 为 50kPa，100kPa，150kPa， 剪胀．作者认为这种现象可 以用材料临近破坏时 200kPa，250kPa，300kPa，350kPa和 400kPa 内部微裂纹的扩展来解释．最近，孙红和葛修润 共八组试验．此外，还进行了 P。各为 150kPa和 等L1]对上海粘性土的CT试验表明，土样是经历 200kPa的两种应力路径下卸载及再加 载试验 了从压密到微裂纹萌生、扩展、连接，直至形成宏 (图 1)． 观裂纹的过程，这进一步证实了上述解释的正确 性．从微裂纹开裂到宏观裂纹和剪切带的形成 ，实 际上是损伤 的演变过程．同时弹塑性变形继续发 刍 8 展 ，必然出现弹塑性变形与损伤相互作用的问题． = Lemaitre等L2]对金属应用应变等价性原理 ，并将 耗散势分解为塑性耗散和损伤耗散 (两者并无耦 合作用)来描述弹塑性变形与损伤的耦合．对于岩 土问题，还需考虑应力路径的影响等． 作者为了求得弹塑性变形与损伤耦合的本构 图1 应力路径示意图 模型，采用 了文献 [3]的数值建模方法．该方法不 需要寻找塑性势和损伤耗散势的解析表达式，直 2 试验结果与分析 接从试验数据 中提炼本构关 系．对 中密砂进行 了 三轴压缩试验 ，获取了弹塑性一损伤耦合的应力应 按照数值建模方法 ，三轴试验的结果可 以绘 变关系，经可视化，绘出了整个应力场 (P，q)中的 出两簇屈服轨迹 ：以塑性体积应变 为硬化参数 剪应变和体应变的三维曲面． 的体积屈服线和以塑性广义剪应变￡p为硬化参数 的剪切屈服线 (图2)．图3显示了应力路径TC下 1 砂土 的三轴试验 的整个应力场 (户，q)中应力应变关系 的全貌． 在 P。一200kPa的试验 中进行 了三次加载 、 作者的三轴试验 以中密砂为试验对象．加载 卸载及再加载试验．分别在剪胀发生前后 ，广义剪 路径为三轴压缩路径TC，即在等 向固结条件下， 应力 q—1．48，q一2：97和 q一4．24三处卸载至零 平均主应力 P为常数的固结排水三轴压缩试验．