粘弹塑本构模型及用于冻土数值计算的柔度矩阵.pdfVIP

第 29 卷 　第 2 期 冰 　川 　冻 　土 Vol . 29 　No . 2 2 0 0 7 年 4 月 J OU RN AL O F GL A CIOL O GY AND GEOCR YOL O GY Ap ril . 2 0 0 7 文章编号 : (2007) 粘弹塑本构模型及用于冻土数值 计算的柔度矩阵 李栋伟 , 　汪仁和 , 　林 　斌 (安徽理工大学 土木工程系矿山建筑重点实验室 , 安徽 淮南 　23200 1) 摘 　要 : 随着寒区经济发展及冻结法施工技术在工程中的推广应用 , 对冻土蠕变的深入研究显得尤为 重要. 对冻土单轴流变试验结果分析 , 用非线性牛顿体替代线性牛顿体 , 从而对西原模型进行了改 进 , 得到了冻土三向应力状态下的本构方程 ; 并利用 DP 屈服准则推导出应用与数值计算的柔度矩阵. 通过有限元程序的二次开发 , 将本构模型添加到大型非线性程序 AD IN A 中. 用三轴蠕变试验检验了 粘弹塑本构模型 , 获得相应的试验拟合曲线. 结果表明: 三轴蠕变试验与理论计算结果吻合 良好 , 说明 粘弹塑本构模型可以描述冻土在高应力下反映的非衰减性蠕变变形特征. 关键词 : 蠕变 ; 屈服准则 ; 冻土 ; 弹粘塑本构方程 中图分类号 : TU452 ; TD265 . 3 文献标识码 : A 性 , 通过试验可以发现 : 当应力较低时蠕变变形收 0 　引言 敛于某一变形值 , 而当应力水平较大时冻土产生加 　　对冻土的物理力学性能进行研究表明 , 冻土是 速蠕变 , 其中加速蠕变时间较短且加速度随时间的 一种极为特殊的材料 , 由于其内部存在着未冻水和 增大而增大 ; 冻土的瞬时变形也呈现出较强的非线 冰使其应力 、应变强烈地表现出时间效应 , 它们之 性. 因此为了描述冻土应力应变随时间变化的全 ( ) 间的相 互作 用造 成 了冻土 强度 蠕变 的复杂 过程 , 可以采用非线性牛顿体替代线性牛顿体 , 从 性[ 1 - 5 ] . 随着寒区建筑物和构筑物及冻结法凿井的 而对西原模型进行改进[9 ] . 改进后的粘弹塑本构模 日益发展 , 冻土地区的建筑物和结构物的稳定性问 型能很好的描述冻土蠕变的各个阶段 , 使其适用于 题 日益受人们关注 , 其中 , 由于冻土蠕变变形过大 非线性分析 , 从而使拟合结果更接近于真实的试验 导致冻结管断裂 , 高速公路 、铁路和大型建筑物的 结果. 冻土地基由于蠕变变形而导致基础开裂 , 影响建筑 1 冻土蠕变屈服准则及试验论证 物和构筑物的正常使用. 因此 , 研究冻土的流变特 性具有十分重要的理论意义与工程实际意义[ 6 - 7 ] . 　　屈服准则是用来确定冻土开始产生塑性变形时 对冻土的流变问题进行数值分析一般需要用弹粘塑 的应力大小 , 为使试验实的条件和现场复杂应力相 性有限元的方法来求解 , 而用弹粘塑性有限元系统 联系 , 采用应力不变量函数来表示 : 来求解蠕变