岩土工程数值分析学习笔记

岩土工程数值分析读书笔记岩土工程数值分析读书笔记 学号 003 姓名 贾红晶 1 20 岩土工程数值分析读书笔记岩土工程数值分析读书笔记 摘要 阅读笔记分为两部分 理论学习和摘要 阅读笔记分为两部分 理论学习和 plaxis 模拟相关问题 模拟相关问题 理论部分理论部分 0 岩土工程数值分析简介岩土工程数值分析简介 岩土工程问题解析分析是以弹塑性力学理论和结构力学作为理论依据 适用于解决连续 介质 各向同性材料 未知量少 边界条件简单的工程问题 存在很大的局限性 岩土工程问题数值分析是借助于计算机的计算能力 适用于解决材料复杂 边界条件复 杂 任意荷载 任意几何形状 适用范围广 岩土工程数值分析发展过程 20 世纪 40 年代 使用差分法解决了土工中的渗流及固结问题 如土坝渗流及浸润线的求法 土坝及地基的固结等 20 世纪 60 年代 使用有限元法成解决了土石坝的静力问题的求解 20 世纪 70 年代 使用有限元法解决了土石坝及高楼 包括地基 的抗震分析 20 世纪 80 年代 边界元法异军突起 解决了半无限域的边界问题 地基的静力及动力问题 都使用边界元法得到了有效地解决 岩土工程数值分析的方法有两类 一类方法是将土视为连续介质 随后又将其离散化 如有限单元法 有限差分法 边界单元法 有限元线法 无单元法以及各种方法的耦合 另 一类计算方法是考虑岩土材料本身的不连续性 如裂缝及不同材料间界面的界面模型和界面 单元的使用 离散元法 不连续变形分析 流形元法 颗粒流等数值计算方法 1 数值分析过程中存在的问题及解决措施数值分析过程中存在的问题及解决措施 问题 1 对岩土工程数值分析方法缺乏系统的知识和深入的理解 出现问题时不知道 在什么情况下属于理论问题或数学模型问题 在什么情况下是属于计算方法问题或本构模型 问题 在什么情况下是参数的确定问题或计算本身的问题等 2 各种本构模型固有的局限性 具有多相性土的物理力学性质太复杂 难以准确地用数学 模型和本构模型描述 例如邓肯一张模型不能反映剪胀性 不能反映压缩与剪切的交叉影响 3 现有的试验手段和设备不能提供适当 合理和精确的参数 靠少数样本点所获得的参数 难以准确地描述整个空间场地的物理力学性能 土的参数因土样扰动难以高质量的获取 其 精度很差 4 数学模型还会给人造成一种错觉 让人觉得其计算结果也一定会更好 更可靠 这样可 能使人们忽略了精确的数学公式也照样会有出错的可能性 只有当输入参数的质量和精度很 高 并能与数学模型的精度相匹配时 才有可能得到较为准确的计算结果 措施 1 加强对土的本构模型的教学与培训 了解和掌握各种土的本构模型的优点和 局限性以及模型参数的离散性 岩土工程数值分析读书笔记岩土工程数值分析读书笔记 学号 003 姓名 贾红晶 2 20 2 在使用数值分析方法的同时 不断地积累使用经验 包括他人的经验 3 在丰富的使用经验基础上 建立相应的使用规范 2 有限元法软件有限元法软件 plaxis 理论学习理论学习 1 模型的选用和分析 Mohr Coulomb 模型 MC 这种模型被推荐用于问题的初步分析 对于每个土层 可 以估计出一个平均刚度常数 由于这个刚度是常数 计算往往会相对较快 可以得到变形的 一个初步印象 五个输入参数 即 表示土体弹性的 E 和 表示土体塑性的 和 c 以及 剪胀角 节理岩石模型 JR 节理模型是一种各向异性的弹塑性模型 特别适用于模拟包括层 理尤其是断层方向在内的岩层行为等 两个参数 和 c Hardening Soil 模型 HS 是一种改进了的模拟岩土行为的模型 适合所有类型的土 不 能用来模拟滞后或者反复循环加载情形 三个参数三轴加载刚度 E50 三轴卸载刚度 Eur和固 结仪加载刚度 Eoed 我们一般取 Eur 3E50和 Eoed E50作为不同土体类型的平均值 软土蠕变模型 SSC 软土蠕变模型是一个新近开发的应用于地基和路基等的沉陷问 题的模型 通常会过高地预计弹性岩土的行为范围 软土模型 SS 软土模型是一种 Cam Clay 类型的模型 特别适用于接近正常固结的粘性 土的主压缩 可以被 Hardening Soil 模型取代 改进的 Cam Clay 模型 MCC 主要用于模拟接近正常固结的粘性土 在实际应用中 是不被推荐的 如果要对所考虑的问题进行一个简单迅速的初步分析 我们建议使用 Mohr Coulomb 模 型 当缺乏好的土工数据时 进一步的高级分析是没有用的 在许多情况下 当你拥有主导土层的好的数据时 可以利用 Hardening Soil 模型 来进行一个额外的分析 毫无疑问 同时拥有三轴试验和固结仪试验结果的可能 性是很小的 但是 原位实验数据的修正值对高质量实验数据来说是一个有益的补充最后 软土蠕变模型可以用于分析蠕变 即 极软土的次压缩 用不同的土工模型来分析同一个 岩土问题显得代价过高 但是它们往往是值得的 首先 用 Mohr Coulomb 模型来分析是相 对较快而且简单的 其次 这一过程通常会减小计算结果的误差 2 各模型的参数 MOHR COULOMB 模型的基本参数 E 杨氏模量 v 泊桑比 内摩擦角 c 内聚力 剪胀角 这些参数可以从土样的基本试验得到 杨氏模量 E 在土力学中 初始斜率用 E0表示 50 强度处的割线模量由 E50表 示 对于具有大范围线弹性行为的材料来 说 使 E0是符合实际的 但是对于土体加 载问题一般使用 E50 如果考虑隧道和开 挖问题中的卸载问题要用 Eur替换 E50 泊桑比 V 在许多情况下得到的 v 值是介于 0 3 和 0 4 之间的 一般地说 除了一维压 Tu 图 1 项目设置 Tu 图 2 强度参数设置 岩土工程数值分析读书笔记岩土工程数值分析读书笔记 学号 003 姓名 贾红晶 3 20 缩 这个范围的值还可以用在加载条件下 在卸载条件下 使用 0 15 和 0 25 之间的值更为普遍 内摩擦角 PLAXIS 可以处理无粘性砂土 c 0 但是有一些选项的执行不太好 为了避免 复杂性 我们建议不熟练的用户至少输入一个较小值 使用 c 0 2kPa 剪胀角 以度的方式指定的 除了严重的超固结土层以外 粘性土通常没有什么剪胀性 0 砂土的剪胀性依赖于密度和摩擦角 对于石英砂土来说 30 的值比 的值小 30 度 然而 剪胀角在多数情况下为零 的小的负值仅仅对 极松的砂土是实际的 HARDENING SOIL 模型的参数 软土模型参数 Tu 图 2 HS 模型强度参数设置 Tu 图 3 软土模型强度参数设置 岩土工程数值分析读书笔记岩土工程数值分析读书笔记 学号 003 姓名 贾红晶 4 20 3 单位和符号规定 符号规定 在平面应变分析里 zz指向平面外 轴对称分析里 x 代表径向坐标 y 代 表轴向坐标 z 代表切向 此时 xx表示径向应力 zz表示环向应力 在所有输出数据里 压应力 包括孔隙压力 和压力设为负值 而拉应力和拉力设为正值 家 建模过程 绘制几何轮廓线 设定边界条 件 添加荷载 添加材料特性 划分网格 设置初始条件 计算工序分层 计算 警警 警告 如果修改已经建好的模型 要意识到必警告 如果修改已经建好的模型 要意识到必 须重新生成网格 有时还包括重新生成初始条件 须重新生成网格 有时还包括重新生成初始条件 以便和修改后的模型一致 以便和修改后的模型一致 Tu 图 4 基本界面 岩土工程数值分析读书笔记岩土工程数值分析读书笔记 学号 003 姓名 贾红晶 5 20 4 荷载问题 荷载子菜单包含的选项 用于在几何模型里引入分布荷载 线荷载或集中荷载以及指定 位移 荷载和指定位移既可用于模型边界 也可用于模型内部 指定位移是为了控制某些点 的位移而强加于模型的特殊条件 大小为零的指定位移即为约束 注意 注意 如果在一条几何线上同时施加指定位移和荷载 那么在计算过程中指定位移优先于 荷载 在全约束线上施加指定位移 约束则优先于指定位移 因而 在全约束线段上没有必 要施加指定位移 集中荷载 分布荷载 如果只有一个位移方向是预加的 而其他方向是自 由的 那么就可以在自由位移的方向上施加集中荷载 分布荷载 5 材料性状种类 材料类型 排水性状 该设置不产生超静水压 很明显 它适用于干土 以及由于高渗透性 砂土 和 或缓慢加载的完全排水条件 该选项也可以用来模拟长期岩土性质 其中不需要模拟不 排水加载和固结的精确历史 不排水性状 该设置用于研究超静水压的完全发展过程 由于低透水性 粘土 和 或 快速加载 孔隙水渗流有时可以忽略 注意 要输入有效的模型参数 非多孔性状 使用该设置的类组 不论是初始孔压还是超静水压 都不予考虑 有关 应用可以在模拟混凝土或结构性状的例子里找到 非多孔性状通常和线 弹性模型组合应用 输入的饱和容重和渗透性参数对非多孔材料是无效 的 6 网格生成 网格的基本单元类型为 15 节点的三角形单元或 6 节点的三角形单元 全局疏密度水平 很粗疏 粗疏 中等 细密和很细密 默认情况下 全局疏密度水平 为粗疏 局部疏密度 在可能出现强烈的应力集中或大变形梯度的区域 最好是用更精确 细密 的有限元元网格来模拟 这通常是在几何模型里出现边缘 角或结构对象的时候用到 这时 候的局部单元尺寸因子就应当等于 0 5 关于网格生成的几点建议 为了优化有限元计算 可以先用相对粗疏的网格作初始分析 这样可以检查几何模型的 大小是否合适 还可以预计发生应力集中和大变形梯度的位置 所得成果主要用来加密有限 元网格 生成加密的有限元模型 为了能够顺利地生成详细的有限元网格 首先应当从网格 子菜单里选择需要的全局疏密度 此外 如果还要进行局部加密 那么应当首先加密类组 然后加密几何线 最后再加密几何点 需要时 可以直接给几何点设局部单元尺寸因子 7 初始条件 岩土工程数值分析读书笔记岩土工程数值分析读书笔记 学号 003 姓名 贾红晶 6 20 在几何模型和有限元网格都建好之后 必须明确初始应力状态和初始构造 初始条件由 两个不同的模式组成 一个模式是生成初始水压 水力条件模式 另一个模式是定义初始 几何构造和生成初始有效应力场 几何构造模式 注意 注意 从初始条件返回到几何图形创建模式是允许的 但是一般不这么做 因为这样做可 能会丢失一些有关初始条件的信息 8 计算 计算分类 分三种基本计算类型 塑性计算 固结分析和 Phi c 折减 安全分析 塑性计算 在弹 塑性变形分析当中 如果不需要考虑超静水压随时间的减小 那么应 当选塑性计算 一般性塑性计算根据未变形的几何图形得出刚度矩阵 这类计算适用于大多 数岩土工程的实际应用情况 固结分析 要分析饱和粘性土的超静水压随时间的变化与消散 应该选择固结计算 一 般情况下 执行完不排水塑性计算之后 再执行不添加另外荷载的固结分析 也可以在固结 分析过程中施加荷载 但是 在接近破坏时要引起注意 因为此时迭代过程可能会不收敛 Phi c 折减 安全分析 可以用减小土的强度参数的方法来执行安全分析 这个过程就叫做 Phi c 折减 是一个单独 的计算类型 如果要计算全局安全系数 就应当选择 Phi c 折减 建议在计算工序列表的末 尾再定义所有安全分析 并且用起始工序号参数注明是在给哪个计算工序求算安全系数 plaxis 例题模拟例题模拟 算例一砂土层上圆形基础的沉降 选用柔性基础 图 1 表示放置在 4m 厚砂土层上半径为 1m 的一个圆形基础 砂土层下是深厚的坚硬岩 石层 计算土体在上部荷载作用下产生的位移和应力 选用轴对称模型 竖直方向取 4m 到 岩石层顶部 水平方向取 5m 5 倍的基础半径 计算过程如下 岩土工程数值分析读书笔记岩土工程数值分析读书笔记 学号 003 姓名 贾红晶 7 20 图 5 新建项目参数设置 图 6 几何模型 图 7 基础模型参数设置 图 8 地基模型的初始应力场 图 9 计算参数设置 放大 图 10 定义计算位置 岩土工程数值分析读书笔记岩土工程数值分析读书笔记 学号 003 姓名 贾红晶 8 20 计算结果分析 基础的沉降值随着距离荷载位置 的增大而逐渐衰减 在 350KN 的均布 荷载下 基础的最大沉降量是 0 129m 最大弯矩是 43 178KN m 最 大剪力是 44 063KN 图 11 变形网格 图 12 变形总位移场 图 13 内力数据 图 14 基础的荷载位移曲线 岩土工程数值分析读书笔记岩土工程数值分析读书笔记 学号 003 姓名 贾红晶 9 20 算例二 不排水条件下河堤的受力分析 如图 1 所示的堤身高 5m 由渗透性差的粘土构成 其地基上部 6m 由软土构成 上层 3m 为一粘土 clay 层 下层 3m 为一泥炭 peat 层 软土层基本上是不透水的 所以河流水位 的短期变化并不影响这一部分的孔隙压力分布 软土层下较大深度范围内为渗透性很好的砂 土 sand 层 其上部 4m 包括在有限元模型之中 这里假设砂土层中的水与河流的水相通 这意味着砂土层的水头紧密跟随河流的水位变化 图 15 新建项目设置 图 16 几何模型 岩土工程数值分析读书笔记岩土工程数值分析读书笔记 学号 003 姓名 贾红晶 10 20 图 17 划分网格 图 18 第一计算工序一般设置 特 特殊位置加密区 图 19 初始应力的设置 图 18 第二计算工序一般设置 岩土工程数值分析读书笔记岩土工程数值分析读书笔记 学号 003 姓名 贾红晶 11 20 图 22 孔隙水压场 图 21 第二计算工序定义孔隙水压力 图 20 第二计算工序一般设置 图 19 初始应力的设置 图 23 因水位变化而产生的位移总量 岩土工程数值分析读书笔记岩土工程数值分析读书笔记 学号 003 姓名 贾红晶 12 20 图 24 水位上升后堤坝的有效应力 应 有效应力几乎为零区域 应 图 25 堤址后有效应力横断面 图 26 水位上升后堤坝的超孔隙水压力 岩土工程数值分析读书笔记岩土工程数值分析读书笔记 学号 003 姓名 贾红晶 13 20 计算结果分析 图 24 可以发现模型右边砂土层顶部的有效应力几乎为零 这是由于砂土层内孔隙压力 的增加造成的 从这一应力图 也可以发现土堤的变形在其后的粘土层中产生了一个被动应 力区 粘土和泥炭层的不排水性能导致了超静水压的发展 图 23 可以看出堤坝在水位变化 情况下 土体的总变形 堤顶和堤址后土体隆起 故在施工过程中需要严格处理堤址后土体 算例三 水下土体开挖施工 鉴于几何模型的对称性 分析时只考虑其中一半 左边 开挖过程模拟为三个不同的 阶段 地下墙使用板单元来模拟 两侧地下墙和土的相互作用用界面单元来模拟 界面单元 可以模拟墙体与土之间产生的滑动摩擦 此摩擦较土之间的摩擦要小些 地下墙之间的横向 支撑用弹簧单元来模拟 岩土工程数值分析读书笔记岩土工程数值分析读书笔记 学号 003 姓名 贾红晶 14 20 材料定义 砂土层 粘土层材料定义如上例 不再叙述 地下连续墙和横向锚杆定义如下 全局网格选用中等级别 在地下连续墙处网格加密 图 28 定义地下连续墙的材料 图 29 定义锚杆的材料 图 30 生成网格 图 31 水压力分布场 图 32 初始土压力分布场 图 33 定义计算工序 图 34 工序 1 激活地下连续 墙 图 27 几何模型 岩土工程数值分析读书笔记岩土工程数值分析读书笔记 学号 003 姓名 贾红晶 15 20 图 30 划分网格 图 31 水压力分布场 图 32 初始土压力分布场 岩土工程数值分析读书笔记岩土工程数值分析读书笔记 学号 003 姓名 贾红晶 16 20 图 33 定义计算工序 图 33 工序 1 激活地下连续墙 岩土工程数值分析读书笔记岩土工程数值分析读书笔记 学号 003 姓名 贾红晶 17 20 图 35 工序 2 第一开挖阶段 图 36 工序 3 激活锚杆 图 37 工序 4 第二开挖阶段 岩土工程数值分析读书笔记岩土工程数值分析读书笔记 学号 003 姓名 贾红晶 18 20 图 38 工序 5 第三开挖阶段 图 39 水下开挖后变形的网格 岩土工程数值分析读书笔记岩土工程数值分析读书笔记 学号 003