1-数值模拟的概念与方法.ppt

一、工程地质数值计算的意义、任务 工程地质数值计算：应用工程地质学原理、现代数学力学方法和计算机技术，研究、解决工程地质（岩土）学中的定量（半定量）化问题的一门学科。 中国的大型和巨型工程建设 300m级高坝 500m级高边坡 10km以上深埋长隧道 西部地形急变带及极端条件下的交通工程 长距离跨流域调水 跨海大桥 跨海峡隧道 。。。。。 工程地质条件的定量化及岩土体赋存环境的模拟 工程动力地质现象及过程的模拟研究 工程岩土体应力场和位移场的分析 岩土体稳定性的模拟 D.边界单元法BEM 70 年代兴起，具有降维作用，精度高 用于非均质岩土体困难较多 暂无商业软件 E.数值流形元法NMM 90 年代兴起，石根华（Discontinuous Deformation Analysis不连续变形分析DDA ） 连续、非连续，结构面 暂无商业软件 三、工程地质数值计算基本程序及其发展 A 有限差分法（FDM） 有限差分法的基本原理与有限单元法类似，只是它们各自的求解方法有所差别。 有限单元法通过刚度矩阵的形式求解每一单元的应力与应变，而在有限差分中，空间离散点处的控制方程组中每一个导数直接由含场变量的代数表达式替换，通过“显式”的方式逐步求解每一单元的应力与应变。 1.3 其它数值模拟方法 软件：FLAC3D、 FLAC2D B 边界单元法（BEM） 边界单元法是20世纪70年代兴起的一种数值方法。 其通过结点之间插值，把边界积分方程转变为线性代数方程组，由此解出各边界单元的结点处待定的边界值，再利用把边界值与域内函数值联系起来的解析公式，求得计算区域内任一点的函数值。 且计算精度、计算效率高，更适用于均质材料和线性性态情况。 软件：Examine2D、Examine3D C 离散单元法（DEM） 岩体往往为众多的节理或结构面所切割，在某些情况下，岩体不能视为连续介质，具有明显的不连续性，很难用连续介质力学方法如有限单元法来处理。 离散单元法是处理非连续介质力学的数值方法，特别适用于节理岩体的应力分析，在土木工程方面应用广泛，尤其在边坡稳定分析方面。 软件：UDEC、3DEC、PFC2D、pfc3D 边界单元法是在有限单元法以后发展起来的一种数值方法。该方法早在20世纪70年代由英国南安普敦大学土木工程系开始使用。该系的C．A．Brebbia在国际上大力倡导边界单元法。现在这个名词已普遍被科学家接受，边界单元法也逐渐被应用到各个领域中。　 由于边界单元法只在研究区域的边界上剖分单元，从而使求解问题的维数降低：三维问题变为二维问题，二维问题变成一维问题。解一个问题所需计算的方程组规模小，有利于节省内存和计算时间。此外，由于边界单元法引入了基本解，具有解析与离散相结合的特点，因而具有较高的精度。 石根华把有限单元方法、非连续变形分析方法(DDA)和解析方法统一起来，创造了流形单元法。流形元法用一层网格套在结构体上，当结构体分成若干块时，块体之间可以产生接触、滑移和张开，而在块体内部，再用网格进行细分，计算出各个节点上的位移和应力，与有限元分析具有同样的效果。因此，流形元方法既能分析结构体内的应力分布状态，同时又能够模拟块体破坏之后的岩块运动规律，具有很强的适应性。 数学覆盖是流形元方法的最基本概念。描述一物体形状函数称为总体函数，总体函数可以通过许多个重叠的覆盖局部区域的覆盖函数来近似地计算出来；物理覆盖的区域是包含在数学覆盖的材料，是数学覆盖与材料的交集。 石根华回顾自己41年的工程师经历，石根华感慨万千：“在这个世界上，主要是靠解决问题的力量，职务、学位、经历等都不太管用。在出现问题时，能解决问题就成功了；失败一次，可能就是永远的失败。成功靠什么来保证？就是数学，在逻辑上靠数学，靠思维的严密，所有的东西，能够用上的，要武装到牙齿。” 　　“做一个真正的工程师，该有胆量时就要有胆量，甚至把自己的生命赌进去。但赌博不是工程师的性格，工程师是要求绝对可靠的，工程师不是赌徒，在任何情况下都要将所有的东西做好。” 　　他对中科院计算数学所的研究生们说：“从采矿、水库大坝到地下隧道工程等，世界各国的工程师面临太多的危险。在这些方面，数学是非常有用的，我们周围的人都需要数学。我希望下一代的数学家们，特别是你们，站在计算数学与工程之间，最重要的是用发明出的一些数学方法和工具，写出很好的教科书，把数学交给工程师，追上这个时代。” 石根华与数学 1、工程地质数值计算基本程序 明确研究目的 工程地质条件的综合分析 地形地貌 地层岩性 地质构造与岩体结构 地下水条件 地应力条件 地震作用 工程地质模型的建立 工程地质体计算范围 工程地质体属性分区与参数确定 岩土体结构类