高等土力学-有限单元法.pptVIP

3 有限单元法的求解方法 补充解释 由于找不到起始点（即零点），所以只能得到土的切线模量，得不到割线模量 ? F 割线模量 切线模量只对局部有效，所以严格讲，哪一荷载对应的模量都不合适 ?F Fn-1 Fn Kulhawy、Duncan等 对一侧限约束的一维柱体进行计算，并与理论解比较 中点增量法误差约为基本增量法的一半，如果分层较密，中点增量法与理论解误差有望5% 表明： 本构关系的积分 3 有限单元法的求解方法 参考：Desai《岩土工程数值方法》 3 有限单元法的求解方法 2 迭代法 直接迭代法： Newton-Raphson迭代法: 用割线模量直接迭代 ? F 不能适用于土 用切线模量进行迭代 ? F 3 有限单元法的求解方法 修正的Newton-Raphson迭代法: ? F ? F 迭代过程中切线模量不变 N-R迭代法 修正的N-R迭代法 常用方法 虽然收敛慢些，但不必每次计算[K]阵，所以整体速度更快 3 有限单元法的求解方法 3 增量迭代法 增量法和迭代法混合，在每一增量步内用迭代法计算 4 有限元计算中的处理方法 4.1 地基土单元处理 4.2 新填土单元处理 4.3 邓肯-张模型的处理 4.4 界面单元 4.5 不可压缩材料 4 有限元计算中的处理方法 4.1 地基土单元处理 4.2 新填土单元处理 方法1： 方法2: 直接给应力赋值，应变、位移置0 用有限元计算，应变、位移置0，只保留应力 1 新填土单元应变、位移置0，只保留应力， 应力可通过有限元计算或直接赋值 2 为防止新填土单元约束下部土体，可取其 模量为下部土体模量的1/500~1/1000 最大竖向位移一般在中间 4 有限元计算中的处理方法 4.3 邓肯-张模型的处理 1 加卸载判断 参考Duncan的FEDAM84程序 SS ? SS上： 加载，模量取 Et SS SS上： 卸载， 弹性模量 切线模量Et 弹性模量Eur ? ? 4 有限元计算中的处理方法 2 土单元受拉时，即?3 0 本人建议： 对砂土，模量取 1/100~1/1000Et 对粘性土， 模量取 1/3Et 左右 模量取小值 1/100~1/1000Et 未拉坏时， 拉坏时， 4 有限元计算中的处理方法 3 ?3 0 ，但?3较小 切线模量 ?3＝0 时无法计算 通常取 计算Et ?0 取 50kPa 或 100kPa 4 关于泊松比 4 有限元计算中的处理方法 4.4 界面单元 两种材料性质相差较大时，两者之间应设界面单元 如土与混凝土，不设界面单元，混凝土中计算应力可能过大 常用： Goodman单元 Desai接触面单元 1 4 3 2 Goodman单元 殷宗泽、朱伯芳书中介绍较好 推导时注意正负号的定义 4 有限元计算中的处理方法 4.5 不可压缩材料 一般情况下： ?＝0~0.5 ? ? 0.5 称不可压缩材料，会出现自锁现象 如果 ? 0.49， 取?＝0.49 如饱和土在不排水情况下 不可压缩材料在力学上有专门的处理方法 5 文献阅读和编程中注意事项 关于文献阅读 ① 书中难免错误和疏漏，应多拿几本书比较着看， 不可只盯一本书 原因 a 问题本身不易说清 b 作者有的熟悉，有的不熟悉 c 作者以讹传讹 d 打印错误 e 其它原因 轮扁斫轮 Kanali Canal 七月流火 Hammer积分 历史原因 lg ln 图表 软件 关于比奥理论的书 床 井台 马踏飞“隼” “里昂”→“第戎” q = Mp 5 文献阅读和编程中注意事项 ② 关键的、基本的原理要自己推导一下 举例 自己推导能使概念清晰 实在无法推导，应多找几本书互相校对一下 ③ 适当做笔记 繁简适当 太简不起作用，太繁令自己疲惫不堪,望而生畏 比奥固结理论及有限元格式推导 三维接触面刚度矩阵、坐标转换关系推导 渗流自由面，5本书5中写法 总之：应当博采众长，自己把握清楚 四面体单元形函数 顺缝渗流，顺杆渗流 5 文献阅读和编程中注意事项 关于程序编写 ① 先大致学习一下别人程序 ② 编写前写好变量说明和程序框图， 至少概要框图，不要拿起来就编 ③ 基本公式、算法要自己推导 以免书中有错，且使自己概念清晰 ④ 做好记录，有来有回 两个本子，一为日常记录，一为阶段性工作整理 ⑤ 要在精力充沛的时段编程序。合理安排时间，不要勉强应付 ⑥ 修改程序时不要随手。要做好修改记录，并事先做好计划 ⑦ 及时整理好文档材料及最终框图 文档材料 内容建议： 理论原理 变量说明 程序框图 处理技巧 繁简适当 5 文献阅读和编程中注意事项 主要流程 读书 读文献 学习相关程序 程序概要框图 编程 总结报告 5 文献阅读和编程中注意事项 关于程序使用 了解其计算原理、处理方法，