基于ABAQUS的红黏土地基排水固结沉降分析.pdf

第 3 2卷 第 2期 2 0 1 5年 2 月 长 江科 学 院 院报 J o u r n a l o f Ya n g t z e R i v e r S c i e n t i fi c Re s e a r c h I n s ti t u t e V0 1 32 NO 2 F e b 2 0 1 5 D OI ： 1 0 3 9 6 9 j i s s n 1 0 0 1 5 4 8 5 2 0 1 5 0 2 0 2 1 基于 A B A Q U S的红黏土地基排水固结沉降分析 赵蕾 ， 陈筠 ， 邬忠虎 ， 刘磊磊 ( 1 贵州大学 资源与环境工程学院， 贵阳5 5 0 0 2 5 ； 2 中南大学 地球科学与信息物理学院， 长沙4 1 0 0 8 3 ) 摘要 ： 针对红 黏土地 区地基排水 固结过程 的复杂性 ， 对 红黏土 地基排 水 固结 沉降进 行数值 分析 。采用 大型通 用 有限元分析软件 A B A Q U S S t a n d a r d中的固结计算模块 S o i l s 和 自带的修正剑桥模型， 以某铝厂的改造工程项 目为 算例 ， 对红黏土地基中的沉降随时间的变化规律 、 孑 L 隙水压力消散规律、 孔隙比随深度的非线性分布以及土体的有 效应 力增长进行 了研究 。研究表 明 ： 计算结果与 土力学规律 吻合 ， 能准确反映红黏 土地基 固结沉 降的实际情况 ， 可 为红黏土地基 上的建 ( 构 ) 物改造工程地基处 理和工程设计 提供参考 。 关键词： 红黏土地基 ； A B A Q U S ； 修正剑桥模型 ； 固结沉降； 孔隙水压力 中图分类号 ： T U 4 4 6 文献标志码 ： A 文章编号 ： 1 0 0 1 5 4 8 5 ( 2 0 1 5 ) 0 2 0 1 0 3 0 5 1 研 究背景 2 0世纪修建 的建 ( 构) 筑物部分以红黏土为地 基 ， 随着城镇化 战略 的实施 ， 要不 断通过对 既有 建 ( 构) 筑物 的改造来提高其使用年限 ， 实现可持续发 展。为了使既有建 ( 构 ) 筑物改造工程在保证 安全 的同时 ， 充分利用土体的固结沉降效应 l 1 ， 使设计 更加经济合理 ， 首要工作则是需要对红黏土地基 固 结沉降随时间发展的规律 、 孔隙水压力消散规律 、 孔 隙比沿着深度 分布的规律要有准确的认识 。因此 ， 红黏土地基固结沉降计算和沉降规律的研究对指导 既有建筑物改造工程 中的地基处理和工程设计有着 重要 的指导 ： 意义 。 目前 ， 红黏土地基沉降计算方法主要有 3类 ： 理论计算法 ， 如分层总和法 、 弹性理论公式法等 ， 虽然地基沉降计 算理论有了很大 的改进 ， 但是 由于 土性材料的复杂性和各向异性 ， 计算过程 中必须采 用一些假设作为前提条件 ， 通过理论公式所计算 出 来的沉降和实测值仍然存在一定差距。对地基进 行沉降观测 ， 根据观测结果可 以较准确地预测地基 最终沉降量及工后沉降量以指导现场施工 ， 这类方 法有 ： T a y l o r 法 、 双曲线法 J 、 指数 曲线法 J 、 A s a o k a 法等 J 。由于这类方法需要 现场监测数据 ， 不但操 作起来 比较困难 ， 而且对施工有一定妨碍 ， 加之施工 过程 中的很多不确定性因素可能导致监测数据的不 准确。数值计 算法 ， 主要有差分法l 、 有 限单元 法 、 离散元 法。其 中， 有 限元法 是 目前运用最广 泛的数值计算方法 ， 其特点是适应性强 ， 可以较全面 地考虑土体的变形及其边 界条件 ， 能够方便地处理 非均质 、 非线性 、 复杂边界条件。基于此 ， 本文利用 大型通用有限元软件 A B A Q U S l 9 强大的非线性分析 能力 ， 结合工程实例 ， 对红黏土地基 固结沉降进行数 值分析 “ 2 工程概况 贵州某工 厂车 间为一层排架 结构 ， 高 1 3 5 m， 属二类安全等级建筑物， 对差异沉降中度敏感。基 础埋置深 度 一1 5 1 3 3 ， 基 础形式为条 形基础。地 貌 为溶蚀残积地貌。场地土质的分布情况如下 。 植物层 ： 主要由褐黄 、 灰褐色耕植土组成 ， 含 植物根系， 结构松散 ， 层厚0 51 5 m。红黏土 ： 褐黄色 ， 硬塑状态 ， 稍湿 ， 层厚0 6 4 3 0 1T I ， 属 中等 偏低压缩性土。红黏土 ： 褐黄色 ， 可塑状态， 湿 ， 层 厚0 92 2 i n ， 属 中等偏 高压缩性土 。红黏土 ： 褐黄色 ， 软塑状态 ， 湿 ， 层厚 1 03 5 IT I ， 属 高压缩 性土。在地基正式加载前 ， 土体在 自重和表面均布 荷载 1 0 k P a 作用下完成 固结 ， 随后地基承受来 自于 建筑物施加的均布荷载 1 4 0 k P a ， 总载荷 1 5 0 k P a 。 为了便于数值分析 ， 依据工程地质勘察报告 ， 将 地基影响深度内的主要地层物理力学参数列于表 1 。 收稿 日期 ： 2 0 1 31 01 1 ； 修回 日期 ： 2 0 1 31 11 0 基金项目： 贵州大学 2 0 1 4年研究生创新基金项 目( 研理工 2 0 1 4 0 1 9 ) ； 贵州省地质灾害监测预警与决策支持平 台建设项 目( z l 1 3 2 7 8 ) 作者简介： 赵蕾 ( 1 9 9 0一) ， 女 ， 贵州贵 阳人， 硕士研究 生， 从事岩土力学 性质及其应用方面 的研究工作 ， ( 电话 ) o 8 5 18 3 6 2 1 8 6 9 ( 电子信 箱) 3 6 1 1 9 3 2 5 4 q q c o m 。 通讯作者 ： 陈筠 ( 1 9 7 0一) ， 女 ， 贵州贵阳人， 副教授 ， 硕士 ， 主要从事区域稳定与岩体稳定 、 岩溶工程地 质、 地基与基础工程、 边坡工程 等方面 的研究 ， ( 电话 ) 0 8 5 1 8 3 6 2 1 8 6 9 ( 电子信箱) 4 0 9 2 7 2 2 7 1 q q C O B。 1 0 4 长江科 学院 院报 2 0 1 5年 3 模型建立 地基变形可简化为平面变形 、 空间渗流问题 ， 因 此 ， 模型选取平面应变进行分析 ， 本构模型选择修正 剑桥模型。 3 1 修正剑桥模型 剑桥模型假定 的能量方程式为 d W o =脾 ( 其 中， 为塑性应变能 ) ， 其确定的屈服轨迹在 P 一 9平 面上是子弹头形的 。首先这种屈服面在各 向 等压试验施加应力增量 d p 0及 d q = 0时， 会产生 塑性 剪 应 变 增 量 及 总 剪 应 变 增 量 ， d = = d e M， 这显然是不合理的。另外 ， 许多试验结果也 表明 ， 用剑桥模型计算 的三轴试验 的应力 一应变关 系与试验结果相差较大。因此本 文采用 B u r l a n d提 出的修正剑桥模型 ( 见 图 1 ) ， 该模型在 p - q平 面上 相应的屈服轨迹为一椭圆。修正剑桥模型的屈服方 程 H。 为 ， 2 、 P ( ) P 。 ( ) 其中， ： 曼 。 ( 2 ) ) 一Sl n 式中： P ， q分别 为土体有效平均正应力和有效偏应 力 ； 为内摩擦角 ； p 具有硬化量 的意义 ； M 为剑桥 模型的参 数 之 一 ， 表 示 临 界 状 态 线 ( c r i t i c a l s t a t e l i n e ， C S L ) 在 p - t 平 面上 的斜率 ， 值可以通过 三轴 不排水剪切或排水剪切试验得出。取塑性体积应变 为硬化参数 ， 将 p 表示成 的函数 ： ， 1+ 、 P = Pa e L v ，P 。 ( 3 ) 式中： p 为初始应力 ， 实际计算中近似取为大气压力； e 为初始孔隙比； A， ， c 为剑 桥模 型的 2个参 数。A， ， ( 值可以利用不同的o r 等向 压缩与膨胀试验给出。 A B A Q U S 中 的C l a y P l a s t i c i t y模 型 正 是 基 于 塑 性理论 ， 在原有的剑桥模型 基础上进行了扩展。修正 的剑桥模 型屈服 面如图 1 剑桥模 型屈服 面 Fi g 1 Yi e l d s ur f a c e on P- t pl a ne o f modi fie d Cambr i dg e m o de l 所示 ， 屈服面函数与土力学 中的表达形式虽略有差 异 ， 但实质上都是一样的， 屈服面函数 。 。 为 石 1 z I 2 。 _ 一 1 ) + ( t ) 一 1 = 0。 ( 4 ) 式中： t 为偏应力参数 ； 0为椭圆与 C S L的交点所对 应的P值 ； 为控制屈服面形状的参数 ， 在 t 的 一 侧 =1 ， 在 t 的一侧 ， 可不等于 1 ， 影响该 屈 服面 的形状 。 3 2模型 描述 该地基的深度影响范 围取 1 0 m， 土层底部小排 水 ， 顶面排水 ， 固结时间总长为 1 0 a 。在基础底部埋 设单点沉降计和孔隙水压力计分别对沉降量和孔隙 水压力变化值进行监测， 由于时间跨度太大， 有 的监 + 2 2 22 e - 1 9 62 e - 1 圈 收 集 到 某 几 年 的 实 测 值 ， ； 但这并不会对结果产生 ， + 1 1 8 3 e - 1 ； 影 响 。 +6 6i4 0e-1 土 层 表 面 作 用 的 均 布 目 三 ? m m m一 _ 蒸 篓 用 固 体 与 液 体 耦 合 的 平 面 一 ” 。 一 L i 应变单元 c P E 8 R P n ， 模型 2 初 初始平衡状态如图2所示。 i i a l e q u b r i u m s t a t e 3 3 初始条件和边界条件 3 3 1 初 始条件 在 A B A Q U S的 S t e p模块巾首先建 名为 的 G e o s t a t i c分析步， 通过 G e o s t a t i c分析步施 加 一 个重 力或体积力 ， 使得土体的初始位移为 0或近似为 0 ， 在 G e o s t a t i c分析步中， A B A Q U S自动获得 给定边 界条件相平衡的应力状态 ， 并将这个应力状态作为 后续分析的初始应力场。 3 3 2 边界 条件 模型中地基 的左右两侧水平位移为 0 ， 不限 向位移 ； 地基底部水平位移和竖 向位移均为 0 。除 地基表面外， 气体边界都为不排水边界 4 计算结果及分析 4 1地 基 沉降分 析 红黏土地基 固结分析 中， 通过分别定义初始孔 隙比、 初始有效应力 和初始孑 L 隙水压力得到地翠 第2期 赵 蕾 等 基于 A B A Q U S的红黏土地基排水固结沉降分析 1 0 5 未建造建筑物前的初始应力场 ， 形成地应力平衡 状 态。由图2 可以看出， 初始平衡状态地基土层表面 处 的最大竖 向位移为2 2 21 0 m， 已达到应力平 衡 状 态 。 计算的基础下表面随时间沉降曲线如图 3所示 ， 从 图中可以看出， 有限元计算结果与土力学变化规律 一 样。从图 3可以看出， 经过 1 a 地基的固结沉降量 占总沉降量的5 0 ， 所用时间占总时间的 1 0 ， 当到 第 5年的时候地基固结沉降量达 9 0 ， 占总时间的 5 0 ， 余下 1 0 的沉降大约用了总时间的 5 0 。 对图 3中的基础下表面沉降 一时间曲线进行拟 合， 得出基础下表面沉降随时间的规律符合 E x p D e c 2模型 ， 公式为 Y = Y o+A1 e一 +A2 e 。 ( 5 ) 式中 ： Y 。 ， t ， A ， t ： ， A ： 均为常数 ， 与土的性质有关 。 图 3 地表沉 降 一时间曲线 Fi g 3 Cu r v e s o f s e t t l e me n t i n e a r t h s ur f a c e V S t i me 地基沉降拟合参数如表 2所示。从图 3中的拟 合 曲线和表 2中的相关系数可 以看出 ， 拟合 曲线与 A B A Q U S有限元计算的模 拟曲线相关性很好 ， 基 础 下表面的沉降规律可以用公式( 5 ) 进行描述。 表2 地基沉降拟合参数 Tab l e 2 Fi t t e d pa r a m e t e r s o f f oun dat i o n s e t t l e me nt 4 2孔隙水压 力分析 孔隙水压力是反映固结发展过程的一个重要指 标。图4表明： 孔隙水压力 的衰减可分为 2个阶段 ， 即快速衰减阶段和缓慢衰减阶段。在上部建筑物荷 载加载初期， 土体中的水不能及时排除， 产生超孔隙 水压力 ， 从图 4可以看出， 土层底部超孔 隙水压力迅 速增大到 1 4 0 k P a左右。在快速衰减阶段 ， 孔 隙水压 力消散了约 8 0 的压力值， 而所用时间为总时间的 5 5 。在缓慢衰减阶段， 孔隙水压力的消散时间占整 个时 间的 4 5 左右 ， 消散量 却 只 占孔 隙水压 力 的 2 0 左右。随着土体 固结时间的增长， 超孔隙水压力 慢慢消散， 到第 1 0 a ， 孔隙水压力基本全部消散。 对图4中的孔隙水压力消散曲线进行拟合 ， 得出 孔隙水压力的消散规律符合 G a u s s 模型 ， 公式为 Y ： Y o+ 。一 2 。( 6 ) + 。 。( o ) ，己 式 中： Y 。 ， 。 ， W， A均为常数 ， 与土的性质有关 。 孔压消散拟合参数如表 3所示 。从图 4中的拟 合 曲线和表 3中的相关系数可以看 出， 拟合 曲线与 有限元计算 的模拟 曲线相关性很好 ， 孔 隙水压力的 消散规律可以用公式( 6 ) 进行描述。 图 4孔隙水压 力消散曲线 Fi g 4 Cu r ve s o f wa t e r pr e s s u r e di s s i pa t i on 表 3孔压消散拟合参数 Ta b l e 3 F i t t e d p a r a me t e r s o f p o r e wa t e r p r e s s u r e d i s s i p a t i o n 4 3孔隙比随深度变化 在实际工程 中， 土体的孔隙比并不是一个常数 ， 而是随深度逐渐减小 。在许多数值模拟中给出土体 顶部和底部 的孔 隙比， 而中间部位的数值则 由这 2 个数值线性插值得到 ， 即土体的孔 隙比沿深度线性 递减 ， 但这种线性变化不符合实 际工程情况 。土体 中孔隙比随深度增加而不 断减小 ， 但却不是线性递 减 ， 孔隙比数值 的变化曲线 近似于二次抛物线 。因 此编写了孔隙比随土体深度非线性分布的用户子程 序 V O I D R I ， 在模 型的运算过程 中使用这个子程序 ， 主要是为了使模拟计算更加接近工程实际情况。 从图5中可以看到初始地应力平衡后， 土体的 孔隙比随深度的增加而非线性递减 ， 这 与工程实际 情况一致。因此可 以得到结论 ： 在模拟 的初始分析 步中添加的子程序 已经产生效果 ， 使孔 隙比的数值 在随土体的深度变化上更加接近于工程实际情况 ， 从而为模拟精度的进一步提高奠定了基础。 长 江科 学 院院报 1 - 2 6 1 2 4 1 2 2 1 2 0 1 - 1 8 1 1 6 1 1 4 1 1 2 1 1 0 1 08 1 O6 深 度 m 图 5孔隙比随深度分布 曲线 F i g 5 Cu r v e s o f v o i d r a t i o V S d e p t h 对 图 5中的孔隙比随深度递减曲线进行拟合 ， 得出孔隙比随深度的递减规律符合 E x p D e c l 模 型， 公式为 Y=Y 0+A l e 。 ( 7 ) 式中 ： 。 A 均为常数 ， 与土的性质有关。 拟合参数如表 4所示 。从图 5中的拟合曲线和 表 4中的相关系数可 以看出 ， 拟合 曲线与有限元计 算的模拟曲线相关性很好， 孔隙比递减规律可以用 公式 ( 7 ) 进行描述。 表 4孔隙 比递减拟合参数 Tab l e 4 Fi t t e d par a m e t e r s o f de c r e as e i n v oi d r a t i o ： 笪 ： 丝 ： ： ! 4 4 有效应力随时间的增长规律 图6是初始应力平衡之后的竖向有效应力场 随 时间 的变 化 ， 在 土体底 部 ， 初 始竖 向有效 应 力为 一 9 0 k P a ， 经过 1 0 a固结之后 ， 土体的有效应力达到 一 2 2 7 k P a ， 增长 了 1 5 2 。在第 4年左右时 ， 土体的 有效应力就达到了 一1 8 5 k P a ， 增长了 1 0 5 ， 土体有 效应力增长了 1倍所用的时间约 4 0 ， 到第 7年之 后土体有效应力达到 一2 2 0 k P a ， 增长 了 1 4 4 ， 后面 3 a 增长率越来越小， 有效应力无太大变化。计算结 果较好地反映了固结的本质 ， 即有效应力随时间的 增长而增加 的过程。 对图 6中的有效应力随时间的增长曲线进行拟 时间 0 1 2 3 4 5 6 7 8 9 1 O 图 6 竖 向有效应力 随时 间的变化 Fi g 6 Var i at i o n o f v e r t i c al e ffe c t i ve s t r e s s wi t h t i nl e 合 ， 用多项式拟合的效果较好 ， 拟合公式为 Y = 0 0 2 0 5 x 一 0 5 70 8 x +5 63 0 4 x 一 21 23 4 x + 1 0 8 0 6 x 一 8 9 8 48 R ：0 9 9 9 8。 ( 8 ) 5 结 论 ( 1 )对红黏土地基进行有 限元分析 ， 土体采用 修正剑桥模型 ， 经过 1 0 a固结 ， A B A Q U S计算结果 为 3 2 m m。得出基础下表面随时间沉降的规律符合 E x p D e c 2模型， 经过 1 a 地基的固结沉降量 占总沉降 量的 5 0 ， 所用时间占总时间的 1 0 ， 最后 1 0 的 沉降大约用 了 5 0 的时间 ， 这符合土力学的基本规 律 ， 表明了修正剑桥模型模拟红黏土是合理 的。 ( 2 )地基经过 1 a固结的沉降量 占总沉降量 的 5 0 ， 在第 5 a 左右沉降量达到 9 0 ， 随后的沉降是 个缓慢的过程。 ( 3 )地 基 土体 孔 隙 水 压力 消散 规 律 可 以用 G a u s s 模型来描述， 孔隙比随深度的变化规律可以 用 E x p D e c l 模型来描述 。 ( 4 )地基土体的有效应力随着时间的发展是一 个逐渐增长 的过程 ， 增长率越来越小 ， 在第 4 a左右 时， 有效应力增长能达到 1倍 ， 在第 7 a左右 ， 有效 应力增长率几乎为0 ， 能较好反映固结本质。 ( 5 )本文的工程实例分析表明， 用 A B A Q U S自 带 的修正剑桥模型模 拟红黏土 固结是完全可行的 ， 具有一定的合理性 和实用性 ， 可 以应用于同类工程 实际工程问题的分析。 参考文献 ： 1 李振明地基土随排水固结而强度增长在工程上的应 用 J 工程勘察， 1 9 9 3 ， ( 3 ) ： 1 31 7 ( L I Z h e n m i n g Ap pl i c a t i o n o f S ub s o i l S t r e n g t h Gr o wt h wi t h Dr a i n a g e C o n s o l i d a t i o n i n E n g i n e e ri n g J G e o t e c h n i c a l I n v e s t i g a t i o n ， 1 9 9 3 ， ( 3 ) ： 1 31 7 ( i n C h i n e s e ) ) 2 G B 5 0 0 0 7 -2 0 1 1 ， 建筑地基基础设计规范 M 北京 ： 中国建筑工业 出版社 ， 2 0 1 1 ( G B 5 0 0 0 7 -2 0 1 1 ， C o d e f o r D e s i g n o f B u i l d i n g F o u n d a t i o n M B e i j i n g ：C h i n a A r c h i t e c t u r e a n d B u i l d i n g P r e s s 2 0 1 1 ( i n C h i n e s e ) ) 3 甘友文 ， 王志亮 ， 郑华地基沉降预测中的双曲线 模型修正 J 水文地质工程地质， 2 0 0 4 ， ( 1 ) ： 9 8 1 0 0 ( G A N Y o u w e n ， WA N G Z h i l i a n g ， Z H E N G H u a Hy p e r b o l i c F o u n d a t i o n S e t t l e me n t P r e d i c t i o n Mo d e l Mo d i fi c a t i o n J H y d r o g e o l o g y a n d E n g i n e e r i n g G e o l o g y ， 2 0 0 4， ( 1 ) ： 9 81 0 0 ( i n C h i n e s e ) ) 4 杜海伟，张建军， 李新生， 等地基沉降计算指数曲线 配合法及其应用 J 土工基础 ， 2 0 0 5 ， 1 9 ( 1 ) ： 5 4 第 2期 赵 蕾 等 基于 A B A Q U S的红黏土地基排水 固结沉 降分析 1 0 7 5 6 7 8 9 5 6 ( D U Ha i w e i ， Z H A N G J i a n - j u n ， L I X i n s h e n g ， e t a 1 Th e Ex p o n e nt i a l Cu r v e Fi t t i ng o f t he Ca l c u l a t i o n o f t h e S u b s o i l S e t t l e m e n t Me t h o d a n d I t s A p p l i c a t i o n J S o i l E n g i n e e ri n g a n d F o u n d a t i o n ， 2 0 0 5 ， 1 9 ( 1 ) ： 5 45 6 ( i n C h i n e s e ) ) 王志亮，黄景忠， 李永池沉降预测中的 A s a o k a法应 用研究 J 岩 土 力 学 ， 2 0 0 6， 2 7( 1 1 ) ： 2 0 2 52 0 2 8 ( WA N G Z h i - l i a n g ， H U A N G J i n g - z h o n g ， L I Y o n g c h i S t ud y o n App l i c a t i on o f As a ok as Me t ho d t o S e t tle me n t P r e d i c t i o n J R o c k a n d S o i l Me c h a n i c s ， 2 0 0 6 ， 2 7 ( 1 1 ) ： 2 0 2 52 0 2 8 ( i n C h i n e s e ) ) 李秀珍 ， 许强， 孔纪名， 等九寨黄龙机场高填方地 基沉降的数值模拟分析 J 岩石力学与工程学报， 2 0 0 5， 2 4 ( 1 2 ) ： 2 1 8 82 1 9 3 ( L I X i u z h e n ， X U Q i ang ， KONG J i - mi n g ， e t a 1 Nu me ri c a l Mo d e l i n g An aly s i s o f S e t t l e me n t s o f Hi g h F i l l F o u n d a t i o n f o r J i u z a i Hu ang l o n g A i r p o r t J C h i n e s e J o u r n a l o f R o c k Me c h a n i c s and E n g i n e e r i n g ， 2 0 0 5 ， 2 4 ( 1 2 ) ： 2 1 8 8 2 1 9 3 (i n C h i n e s e ) ) 郭志川 ， 刘宁地基固结沉降随机有限元计算和可 靠度分 析 J 岩土力 学， 2 0 0 1 ， 2 2( 4) ： 4 8 14 8 5 ( G U O Z h i c h u a n ， L I U N i n g T h e S F E M A n al y s i s a n d R e - l i a b i l i t y As s e s s me n t o f C o n s o l i d a t i o n S e t t l e me n t B a s e d o n B i o t C o n s o l i d a t i o n T h e o r y J R o c k and S o i l Me c h a n i c s ， 2 0 0 1 ， 2 2 ( 4 ) ： 4 8 1 4 8 5 ( i n C h i n e s e ) ) 杨坪， 杨军 ， 许德鲜， 等基于 A D I N A的软土坝基 沉降分析 J 工程地质学报， 2 0 0 8 ， ( 4 ) ： 5 3 35 3 8 ( Y A N G P i n g ， Y A N G J u n ， X U D e x i a n ， e t a 1 F i n i t e E l e me n t An aly s i s o f S e t tl e me n t o f S o f t S o i l F o u n d a t i o n u n d e r R e s e r v o i r D a m J J o u rna l o f E n g i n e e r i n g G e o l o g y ， 2 0 0 8 ， ( 4 ) ： 5 3 3 5 3 8 ( i n C h i n e s e ) ) 费康， 张建伟A B A Q U S在岩土工程中的应用 M 北京 ：中 国水利 水 电出版 社，2 0 1 0 ( F E I K a n g ， Z H A N G J i a n w e i A p p l i c a t i o n o f A B A Q U S i n G e o t e c h n i - 1 0 1 1 1 2 1 3 1 4 1 5 c a l E n g i n e e ri n g M B e ij i n g ： C h i n a Wa t e r P o w e r P r e s s ， 2 0 1 0 ( i n C h i n e s e ) ) 王金昌， 陈页开A B A Q U S在土木工程中的应用 M 杭州 ：浙 江 大 学 出 版 社 ，2 0 0 6 ( WA N G J i n c h a n g ， C H E N Y e k a i A p p l i c a t i o n o f A B A Q U S i n C i v i l E n g i - n e e r i n g M Hang z h o u ： Z h e j i a n g U n i v e r s i t y P r e s s ， 2 0 0 6 ( i n C h i n e s e ) ) 李广信高等土力学 M 北京 ：清华大学出版社 ， 2 0 0 4( L I G u a n g x i n Ad v a n c e d S o i l Me c h a n i c s f M B e i j i n g ： T s i n g h u a U n i v e r s i t y P r e s s ， 2 0 0 4 ( i n C h i n e s e ) ) 龚晓南高等土力学 M 杭州：浙江大学出版社 ， 1 9 9 6 ( G O N G X i a o n a n A d v a n c e d S o i l Me c h a n i c s M Ha n g z h o u ： Z h e j i a n g U n i v e r s i t y P r e s s ， 1 9 9 6 ( i n C h i n e s e ) ) 廖红建岩土工程数值分析 M 北京 ： 机械工业出 版 社 ， 2 0 0 9 ( L I A O H o n g - j i a n N u me ri c al A n al y s i s o f G e o t e c h n i c a l E n g i n e e ri n g M B e i j i n g ：C h i n a Ma c h i n e P r e s s ， 2 0 0 9 ( i n C h i n e s e ) ) 周葆春， 王靖涛 ， 杨晓东武汉黏土修正剑桥模型参数 研究 J 武汉理工大学学报， 2 0 0 7 ， 2 9 ( 3 ) ： 1 0 31 0 7 ( Z H O U B a o - c h u n ， WA N G J i n g t a o ， Y A N G X i a o d o n g Re s e a r c h o n P a r a me t e r s o f Mo d i f i e d C a mb rid g e Mo d e l for Wu h a n C l a y l J 1 J o u rna l o f Wu h a n U n i v e r s i t y o f T e c h n o l o g y ， 2 0 0 7 ， 2 9 ( 3 ) ： 1 0 31 0 7 ( i n C h i n e s e ) ) 高晖A B A Q U S在软基固结过程分析中的应用研究 D 武汉： 武汉理工大学， 2 0 0 6 ( G A O H u i S t u d y o n A p p l i c a t i o n o f A B A Q U S i n A n a l y s i s d u ri n g C o n s o l i d a t i o n o f S o ft S o i l F o u n d a t i o n f D 1 Wu h an ： Wu h a n U n i v e r s i t y o f T e c h n o l o gy， 2 0 0 6 ( i n C h i n e s e ) ) An a l y s i s o n Dr a i n a g e Co n s o l i d a t i o n S e t t l e me nt o f ( 编辑 ： 黄玲 ) Re d C l a y F o u n d a t i o n b y Us i n g AB AQUS Z HAO L e i ， C HEN J u n ，W U Z h o n g h u ， L I U L e i l e i ( 1 C o l l e g e o f R e s o u r c e a n d E n v i r o n m e n t a l E n g i n e e ri n g ， G u i z h o u U n i v e r s i t y ， G u i y a n g 5 5 0 0 2 5， C h i n a ； 2 S c h o o l o f G e o s c i e n c e a n d I n f o - P h y s i c s ， C e n t r al S o u t h U n i v e r s i t y ， C h a n g s h a 4 1 0 0 8 3 ， C h i n a ) Abs t r a c t： Du e t o t h e c o mp l e x i t y o f d r a i n a g e c o n s o l i d a t i o n o f r e d c l a y f o u nd a t i o n，we s i mu l a t e d a n d a n aly z e d t h e c o n s o l i da t i o n s e t t l e me n t b y u s i n g t h e c o ns o l i d a t i o n c a l c u l a t i o n mo d u l e S o i l s a n d mo d i f i e d Ca mb r i dg e mo d e l o f A B A Q U S S t a n d a r d T h e r e c o n s t r u c t i o n p r o j e c t o f a n a l u m i n u m f a c t o r y w a s t a k e n a s a n e x a m p l e T h e v a ri a t i o n o f r e d c l a