岩石流变力学-数值.pptVIP

5 岩石流变问题数值方法 然后求： ③若 则 5 岩石流变问题数值方法 否则： 9.根据 求 然后求以下积分 求和对所有高斯点 10. 进行收敛性检查 岩石流变力学 1 绪论 2 岩石流变的力学特性 3 岩石流变本构理论 岩石流变问题解析方法 5 岩石流变问题数值方法 6 岩石流变室内试验 7 岩体蠕变现场试验与变形时效监测 8 岩石流变问题的工程应用 5 岩石流变问题数值方法 主要介绍弹——粘塑性问题的有限元解法。 平衡方程： 在tn~tn+1时间间隔内（Δtn） 5 岩石流变问题数值方法 对弹－粘塑性变形，总应变： 5 岩石流变问题数值方法 平面应力： 平面应变： 5 岩石流变问题数值方法 则在Δtn ＝tn+1-tn内，采用Euler时间积分法 令 流动矢量 设tn时刻为 设tn＋1时刻为 其中： 5 岩石流变问题数值方法 将 展开为Taylor级数 其中： 5 岩石流变问题数值方法 则： 在Δtn 内应力增量 其中 由前述 ① 5 岩石流变问题数值方法 则： 式中： ② 将②式代入平衡方程： 上式可简写为： 或： ③ 式中： 5 岩石流变问题数值方法 将 代入②求得： 将 代入①求得： 于是在tn+1时刻的 5 岩石流变问题数值方法 平衡修正： 由于 并不完全满足平衡方程 设 为第n个迭代步结束时的不平衡力，将 该力转移到第n+1个荷载增量上去 5 岩石流变问题数值方法 ⅰ ⅱ 对所有高斯点 时间步长： 5 岩石流变问题数值方法 ⅲ 对Tresca 当 对显式时间积分 时 对Mises 对M-C 5 岩石流变问题数值方法 对8节点四边形等参单元 对每个单元 数值积分： 为16 ×16矩阵 为计算方便，将 分解为8×8个2阶子阵， 即： 为2 ×2矩阵 5 岩石流变问题数值方法 令： 则： 5 岩石流变问题数值方法 对单元内n×n个抽样点： 上述积分： 为加权因子 对体积力，边界力转化为单元节点力的计算类同 5 岩石流变问题数值方法 H矩阵的计算： 5 岩石流变问题数值方法 其中： 5 岩石流变问题数值方法 5 岩石流变问题数值方法 对平面问题： 5 岩石流变问题数值方法 5 岩石流变问题数值方法 等效节点荷载计算： ⅰ 重力荷载 5 岩石流变问题数值方法 ⅱ 边界分布荷载 由微单元平衡得： 由于： 于是： 令 则 5 岩石流变问题数值方法 求解步骤： 求解从t=0开始。t=0时为静弹性解其位移 应力 应变 ，初始屈服 均已知， ，设tn时刻的 已求出 1. 计算n+1步单元节点荷载 5 岩石流变问题数值方法 2. 计算 为时间步进参数 5 岩石流变问题数值方法 3. 组集总刚矩阵 计算位移增量 4.由 计算应变增量 由 计算应力增量 5 岩石流变问题数值方法 5. 6.求高斯点 7.进行屈服判断，计算 5 岩石流变问题数值方法 8. 调整时间步长 ①比较 和总等效应变 是否小于给定的值， 若是，就不计算新的时间步长；否则 ②计算等效粘塑性应变率