数值分析在岩土工程中的应用.pptx

数值分析方法 在岩土工程中的应用简介;数值分析方法 在岩土工程中的应用简介;1、基本概 念; 随着计算机技术的发展，科学技术的进步，科学与工程计算(简称科学计算)的应用范围已扩大到许多学科领域，形成一些边缘学科，例如计算物理、计算化学、计算力学等。目前，实验、理论、计算已成为人类进行科学活动的三大方法。 为了解某科学与工程实际问题，首先是依据物理、力学规律建立问题的数学模型，这些模型一般为代数方程、微分方程等。科学计算的一个重要方面就是要研究解这些数学问题的数值计算方法(适合计算机计算的计算方法)，然后通过计算软件在计算机上计算出实际需要的结果。数值分析内容包括；函数的插值与逼近方法，微分与积分计算方法，线性方程组与非线性方程组计算方法，常微分与偏微分数值解等。; “应用力学”领域内的很多课题，工程师和科学工作者近年来一般用数值方法来求解，获得显著成功。这主要由于目前力学课题本身的复杂性——非均质、非线性以及复杂的加荷条件及边界条件，精确解已无能为力。 ;力学求解方法; 图中数值方法列出最常用的五种：差分法、有限无法、边界元法、变分法和加权余量法。 差分法由英国学者Southwell提出，是转化课题的常微分方程或偏微分方程为差分方程，然后结合初始及边界条件，求解线性代数方程组。这方法比校直观，容易编制程序，所以从20世纪40年代末期以来，迄今盛行不衰。 有限元法由英国学者Zienkiewicz提出，以研究区域的变分形式，离散化所研究的区域，使其成为有限数目的单元。对于边界复杂以及材料属非线性的课题，与差分法相比，有特殊的灵活性。 边界元法由英国学者Brebbia提出，化微分方程为边界积分方程，使用类似于有限元法的离散技术来离散边界。离散化所引起的误差仅来源于边界，因之提高了计算精度，依靠边界节点上算得的量，即可计算区域内的有关物理量，从而减少了准备工作量及计算量。边界元又有直接法及间接法之分，间接法需要先求取一个虚设的量，多了一道手续，所以相对麻烦一些。边界元的缺点是对变???数或非线性问题的造就性不如有限元法。 香港大学张佑启开展了结构力学有限元法的卓越研究(例如有限条分法等，同时他在弹性地基上的研究，事实上奠定了边界元法的基本理论。; 变分法是讨论泛函的极值问题，对上述差分法及有限元法都可起推导基本公式的作用，而这方法本身，也是数值方法中最古老的方法。 加权余量法，可以引入试函数和权函数的方法，从微分方程中直接求出近似的数值解。它的优点是可以避免建立能量方程，使一些无法求得能量方程的课题，也得到了较精确的解答。;2.1 有限元Finite Element Method (FEM);2.1 有限元Finite Element Method (FEM); 研究表明，对于大多数岩土工程问题，无论是进行总应力分析还是进行有效应力分析，均可归结为对Biot动力固结方程的求解。因此，可将Biot动力固结方程作为岩土工程问题的总控制方程。 岩土工程中的基本方程包括土体平衡(或运动)方程、物理(或本构)方程、几何方程、有效应力原理、孔隙流体(水)平衡方程、连续方程等。总控制方程(即Biot动力固结方程)由这些基本方程组合而成。 总控制方程的推导基于以下假定： (1) 土体是完全饱和的横观各向同性弹性体。 (2) 土体的变形是微小的。 (3) 土颗粒和孔隙水不可压缩。 (4) 孔隙水相对于土骨架的渗流运动服从Darcy定律，其惯性力可不计。 (5) 应力应变的正负号法则与弹性力学相反。;2.2 离散单元Discrete Element Method (DEM) ; 近几十年来，离散单元法有了长足的发展，已成为解决岩土力学问题的一个重要的数值方法，越来越受到人们的重视。因为工程中所见到的岩体其形态常呈非连续结构，所形成的岩石块体运动和受力情况多是几何或材料非线性问题，所以很难用解决连续介质力学问题的有限单元法或边界单元法等数值方法来进行求解，而离散单元法正是充分考虑到岩体结构的不连续性，适用于解决节理岩体力学问题。离散单元法除了用于边坡、采场和巷道的稳定性研究以及颗粒介质微观结构的分析外，已扩展到用于研究地震、爆炸等动力过程和地下水渗流、热传导等物理过程。 离散单元法与其他数值方法(如有限单元法、边界单元法等)耦合更能发挥各自方法的优点。例如，用边界单元法考虑远场应力的影响以模拟弹性的性质，用有限单元法作为中间过渡考虑塑性变形，再用离散单元法考虑近场不连续变形的