连续-非连续性岩体的数值方法.pptVIP

图1-1(b) DDA计算的边坡不稳定的情况 1、DDA：边坡落石分析 图1-2(a) 9号国道A017断面边坡坍塌安全分析评价的DDA模型（Step=0 ） （日本）国道９号美方郡温泉歌長、養父郡関宮町尾崎斜面落石、崩壊検討 图1-2(b) Step=142 1、DDA：边坡稳定分析 图1-3(a) 金居原水电站水库边坡潜在滑动面稳定分析的DDA模型（Step=0 ） （日本）金居原発電所潜在的岩盤地すべり安定率研究 图1-3(b) Step=800 1、DDA：边坡稳定分析 图1-4(a) 殿山水电站水库岩体边坡滑动分析的DDA模型 （Sliding Face: B ） （日本）殿山発電所下露地すべりDDA解析、安定性評価 图1-4(b) 滑动体的位移随滑动面上接触参数变化的情况 1、DDA：洞室围岩稳定分析 图1-5(a) 龙滩水电站地下厂房开挖的DDA模型（开挖前） 仅考虑断层的情况 图1-5(b) 开挖后 图1-6(a) 开挖前 考虑断层和节理的情况：节理平均密度为1条/2.5m和1条/5.8m 图1-6(b) 开挖后 2、NMM：边坡稳定分析 图2-1 NMM概化模型（其中的三角形是数学覆盖） （日本）下小鳥ダム湛水池不連続岩盤斜面の安定性評価 Step=100，稳定情况 图2-2 边坡滑动 图2-3 边坡失稳 图2-4 拱座变形失稳 图2-5 挤压变形破坏 图2-6 裂缝扩展 3、EFM：半无限平面受集中力作用 考虑平面应变情况。分别用有限元法及无单元法进行计算。有限元采用96个节点、25个8结点矩形二次单元，无单元法采用同一组节点，积分域为10×10个1m×1m的矩形子域，每个子域采用2×2高斯积分，基函数选用一次基。 图3-1(a) 半无限平面受集中力作用 半无限平面受集中力作用 图3-1(b) x=2m 处应力精确解与计算值的比较。 采用一次基函数的无单元法计算结果比采用二次元的有限元法计算结果更接近于精确解，尤其是剪应力 3、EFM：边坡应力分析 对四川锦屏水电站普斯罗沟坝址边坡应力进行了弹性分析。普斯罗沟基岩裸露，岩壁耸立，为典型的深切V型河谷，是一优越的拱坝坝址，但地质勘探显示，在左岸某些部位，具有深拉裂区，对其成因需进行探讨。 计算采用1915个结点，50×50m2矩形积分子域，每个子域内采用3×3高斯积分，基函数为一次基，考虑平面应变情况。 图3-2(a) 计算范围及岩体分区 图3-2(b) 拉应力区和第一主应力等值线 在不同岩体的交界处产生了较大片的拉应力区，这些区域分布与地质勘探的裂隙区分布接近。可见，裂隙区的形成主要是由于软硬岩层间错动引起拉应力所致。 边坡应力分析 3、EFM：拉剪裂缝扩展 图3-3 拉剪裂缝扩展 (a)一边长10m的含裂缝正方形板，裂缝两端坐标为(4.5,4.5)和(5.5,5.5)，板受单轴拉应力1kPa的作用。 (b)开裂跟踪计算10步后的裂缝轨迹。(c)开裂过程中有效应力强度因子随开裂长度的变化曲线。 材料的断裂韧度为0。计算中采用一次基函数，积分子域采用由结点为角点的积分子域，在每个子域上采用5×5高斯积分。可见，拉剪情况下，有效应力强度因子随裂缝扩展不断增大，一旦裂缝起裂，将扩展失稳。 3、EFM：压剪裂缝扩展 图3-4 压剪裂缝受力示意图 (a) 侧向拉力 (b) 无侧向力 (c) 侧向压力 材料的断裂韧度为0，裂缝面的摩擦系数为0，凝聚力为0。计算中采用一次基函数，积分子域采用由结点为角点的积分子域，在每个子域上采用5×5高斯积分。 图3-5 计算30步后的裂缝轨迹 (a) F=－1kPa，P=100Pa (b) F=－1kPa，P=0 (c) F=－1kPa，P=－100Pa 谢谢！ 人有了知识，就会具备各种分析能力， 明辨是非的能力。 所以我们要勤恳读书，广泛阅读， 古人说“书中自有黄金屋。 ”通过阅读科技书籍，我们能丰富知识， 培养逻辑思维能力； 通过阅读文学作品，我们能提高文学鉴赏水平， 培养文学情趣； 通过阅读报刊，我们能增长见识，扩大自己的知识面。 有许多书籍还能培养我们的道德情操， 给我们巨大的精神力量， 鼓舞我们前进。 * 6、非连续变形分析(DDA) 石根华和Goodman于1989年提出了一种称为非连续变形分析 (Discontinuous Deformation Analysis — DDA)的数值方法，它是在非连续体位移反分析法的基础上推广而来的一种正分析方法，它可以从块体结构的几何参数、力学参数、外荷载及约束情况计算出块体的位移、变形、