数值分析在岩土工程中的应用
合集下载
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
何谓数值分 析
随着计算机技术的发展,科学技术的进步,科学与工程计算 简称科学计 算 的应用范围已扩大到许多学科领域,形成一些边缘学科,例如计算物理、 计算化学、计算力学等，目前,实验、理论、计算已成为人类进行科学活动 的三大方法，
为了解某科学与工程实际问题,首先是依据物理、力学规律建立问题的 数学模型,这些模型一般为代数方程、微分方程等，科学计算的一个重要方 面就是要研究解这些数学问题的数值计算方法 适合计算机计算的计算方 法 ,然后通过计算软件在计算机上计算出实际需要的结果，数值分析内容 包括；函数的插值与逼近方法,微分与积分计算方法,线性方程组与非线性 方程组计算方法,常微分与偏微分数值解等，
岩土工程中最基本的两种分析方法是总应力分析 Total stress Ana1ysis 法和有效应力分析 Effective stress Analysis 法，与此相应,岩土工程中的有 限元法也可分为总应力分析有限元法和有效应力分析有限元法，l 966 年,cloughu等首先将总应力分析有限元法用于土坝的应力和变形分析，1969 年,sandhu和wilson用有限元法分析了Biot二维固结问题,开创了岩土工程有 效应力分析有限元法的先河，在国内,沈珠江 1977 首先将有效应力分析有 限元法应用于软土地基的固结变形分析，
香港大学佑启开展了结构力学有限元法的卓越研究 例如有限条分法等, 同时他在弹性地基上的研究,事实上奠定了边界元法的基本理论，
1、基本 概 念
数值分析的主要求解 方 法
变分法是讨论泛函的极值问题,对上述差分法及有限元法都可起推导基 本公式的作用,而这方法本身,也是数值方法中最古老的方法，
加权余量法,可以引入试函数和权函数的方法,从微分方程中直接求出近 似的数值解，它的优点是可以避免建立能量方程,使一些无法求得能量方程 的课题,也得到了较精确的解答，
Lemos于1983年开发了离散单元法与边界单元法耦合的半平面程序,并 用于计算节理和断裂介质中的应力分布问题，
2、几种常见数值分析方法的主要 特 点
2.2 离散单元法Discrete Element Method DEM
Lorig于1984年开发了包括前处理和后处理的离散单元法与边界单元法 耦合程序， 翌年他到澳大利亚英联邦科学与工业发展组织的岩土力学研究 所,修改了他原先的程序，这个修改后的程序文本称为HYDEBE HYBRID DISCRETE ELEMENT BOUNDARY ELEMENT ,其功能更强,包括一个前 处理程序CREATE,类似于有限单元法程序中的自动划分网格,一个与边界单 元法耦合程序BOUND和一个离散单元法程序BLOCK，
数值分析在岩土工程中的应用
Add the author and the accompanying title
数值分析方法 在岩土工程中的应用简介
A Brief Introduction to the Application of Numerical Analysis Method in the Geotechnical Engineering
差分法由英国学者Southwell提出,是转化课题的常微分方程或偏微分方 程为差分方程,然后结合初始及边界条件,求解线性代数方程组，这方法比 校直观,容易编制程序,所以从20世纪40年代末期以来,迄今盛行不衰，
有限元法由英国学者Zienkiewicz提出,以研究区域的变分形式,离散化所 研究的区域,使其成为有限数目的单元，对于边界复杂以及材料属非线性的 课题,与差分法相比,有特殊的灵活性，
2005年度辽宁省建设厅执业资格注册中心 注册土木工程师 岩土 继续教育培训专用
陈剑 平
吉林大学朝阳校区建设工程学 院
原长春地质学院水文地质与工程地质系 或长春科技大学环境与建设工程学院
2005年2月25日
讲授提纲
1、基本概 念 2、几种常见数值分析方法的 主要特点 3、对数值分析定量结果的理解
1、基本概 念
2、几种常见数值分析方法的主要 特点 2.2 离散单元Discrete Element Method DEM
一般从宏观意义上说,岩石可以视为连续介质,从而可以用弹性力学或塑 性力学的方法来进行分析和计算，但在某些情况下,岩体却不能视做连续介 质,如地下节理岩体中的巷道 见下图 ,这时,就不宜用处理连续介质的力学方 法来进行计算，于是,离散单元法作为一种处理节理岩体的数值方法就应运 而生，
2、几种常见数值分析方法的主要 特 点
2.2 离散单元法Discrete Element Method DEM
至于三维离散单元法的发展则要迟些,共主要原因是数据结构复杂,要 求计算机应具有较大的容量,并且计算结果的图形显示较为困难,如切一剖 面,块体间一般来讲是不接触的,犹如浮在空中,看起来不直观，三维离散单 元法程序3DEC 3—DIMENSIONAL DISTINCT ELEMENT CODE 已由Cundall与ITASCA咨询集团于1986年开发出来，其基本原理同 UDEC一样,只是数据结构作了较大的改进，
1、基本概 念
数值分析的必要 性
应用力学领域内的很多课题,工程师和科学工作者近年来一般用数值方 法来求解,获得显著成功，这主要由于目前力学课题本身的复杂性——非均 质、非线性以及复杂的加荷条件及边界条件,精确解已无能为力，
数值分析的可能性 计算机的迅速发展,也使数值分析得到有效而经济的成果，
1、基本概 念
离散单元法的基本思想 离教单元法也像有限单元法那样,将区域划分成单元 见图1—3 a ，但
是,单元因受节理等不连续面控制,在以后的运动过程中,单元结点可以分离, 即一个单元与其邻近单元可以接触,也可以分开 见图I—3 b ， 单元之间相 互作用的力可以根据力和位移的关系求出,而个别单元的运动则完全根据该 单元所受的不平衡力和不平衡力矩的要求解方法
在量测手段得到改进的今天,少数课题利用实验室的结构模型试验或离 心模型试验也可以得到定性或甚至定量的结果，这样,上述力学求解方法可 列成下表：
力学求解方法
精确解 数值方法
实验手段
差分法 有限元法
边界元法 变分法
加权余量法
1、基本概 念
数值分析的主要求解方法
图中数值方法列出最常用的五种：差分法、有限无法、边界元法、变 分法和加权余量法，
2、几种常见数值分析方法的主要 特点
2.1 有限元Finite Element Method FEM 对于土体,总应力分析法是将其视为固体来分析的，因此岩土工程中的
总应力分析有限元法与一般固体力学中的有限无法是相同的，有效应力分 析法则严格区分土体中分别由土颗粒骨架、孔隙水和孔隙气传递或承受的 应力,并考虑土骨架变形、孔隙水压力消散和孔隙气压力消散三者的耦合作 用,因而比总应力分析法更接近实际,但远较之复杂，对于岩体,因其本身可 视为固体,故与之相应的分析只能是总应力分析，从有效应力原理出发,可将 总应力分析视为有效应力分析的一种特殊形式，因此,总应力有限元分析法 也可视为有效应力分析有限元法的一种持殊形式，
正确的剖分
不正确的剖分
2、几种常见数值分析方法的主要 特点
2.1 有限元Finite Element Method FEM
研究表明,对于大多数岩土工程问题,无论是进行总应力分析还是进行有 效应力分析,均可归结为对Biot动力固结方程的求解，因此,可将Biot动力固 结方程作为岩土工程问题的总控制方程，
边界元法由英国学者Brebbia提出,化微分方程为边界积分方程,使用类似 于有限元法的离散技术来离散边界，离散化所引起的误差仅来源于边界,因 之提高了计算精度,依靠边界节点上算得的量,即可计算区域内的有关物理 量,从而减少了准备工作量及计算量，边界元又有直接法及间接法之分,间 接法需要先求取一个虚设的量,多了一道手续,所以相对麻烦一些，边界元 的缺点是对变系数或非线性问题的造就性不如有限元法，
2、几种常见数值分析方法的主要 特 点
2.2 离散单元法Discrete Element Method DEM
离散单元法在国内外的发展状况 离散单元法一般认为是Cundall于1971年提出来的，该法适用于研究在
准静力或动力条件下的节理系统或块体集合的力学问题,最初用来分析岩 石边坡的运动，到1974年,二维的离散单元法程序趋于成熟,当时已有屏幕 图形输出的交互会话功能，但由于受计算机内存的限制,不少程序是用汇 编语言写成,到1978年才全部翻译成FORTRAN IV的文本,成为离散单元法 的基本程序，与此同时,Cundall和 Strack还开发了二维圆形块体的BALL程 序,用于研究颗粒介质的力学行为,所得结果与Drescher等人用光弹技术的实 验结果极为吻合,使BALL程序在研究颗粒介质的本构方程方面大放异彩，
Cundall于1980年就开始研究块体在受力后变形以及根据破坏堆则允许 断裂的离散单元法, 这显然是将块体视为刚体的离散单元法的一个进步， Cundall称这种方法为UDEC UNIVERSAL DISTINCT ELEMENT CODE ,该 程序最后于1985年完成，UDEC现已广泛用于岩土力学和采矿工程,被公认 为对节理岩体进行数 值模拟的一种行之有效的方法，
2、几种常见数值分析方法的主要 特点
2.1 有限元Finite Element Method FEM
在众多的工程数值计算方法中,有限单元法 Finite Element Method,缩写 为FEM 很早就以其适用性强和处理非均质、非线性、复杂边界诸多问题方 便等突出优点而成为工程数值分析最有效的通用工具，经过近半个世纪的 发展,有限单元法已十分成熟并在各个领域的工程分析中广泛应用，
2、几种常见数值分析方法的主要 特点
2.2 离散单元法Discrete Element Method DEM
离散单元法是一种显式求解的数值方法，该方法与在时域中进行的其他 显式计算相似,例如与解抛物线型偏微分方程的显式差分格式相似， 显式是 针对一个物理系统进行数值计算时所用的代数方程式的性质而言，在用显式 法计算时,所有方程式一侧的量都是已知的,而另一侧的量只要用简单的代入 法就可求得，这与隐式法不同,隐式法必须求解联立方程组，在用显式法时, 假定在每一迭代时步内,每个块体单元仅对其相邻的块体单元产生力的影响, 这样,时步就需要取得足够小,以使显式法稳定，由于用显式法时不需要形成 矩阵,因此可以考虑大的位移和非线性,而不必花费额外的计算时间，
离散单元法与其他数值方法 如有限单元法、边界单元法等 耦合更能发挥 各自方法的优点，例如,用边界单元法考虑远场应力的影响以模拟弹性的性质, 用有限单元法作为中间过渡考虑塑性变形,再用离散单元法考虑近场不连续变 形的情况,从而极大地扩展了数值方法的解题范围，
2、几种常见数值分析方法的主要 特 点
2.2 离散单元法Discrete Element Method DEM
岩土工程中的基本方程包括土体平衡 或运动 方程、物理 或本构 方程、 几何方程、有效应力原理、孔隙流体 水 平衡方程、连续方程等，总控制方 程 即Biot动力固结方程 由这些基本方程组合而成，
总控制方程的推导基于以下假定： 1 土体是完全饱和的横观各向同性弹性体， 2 土体的变形是微小的， 3 土颗粒和孔隙水不可压缩， 4 孔隙水相对于土骨架的渗流运动服从Darcy定律,其惯性力可不计， 5 应力应变的正负号法则与弹性力学相反，
可行的
2、几种常见数值分析方法的主要 特点
2.2 离散单元Discrete Element Method DEM
近几十年来,离散单元法有了长足的发展,已成为解决岩土力学问题的一个 重要的数值方法,越来越受到人们的重视，因为工程中所见到的岩体其形态常 呈非连续结构,所形成的岩石块体运动和受力情况多是几何或材料非线性问题, 所以很难用解决连续介质力学问题的有限单元法或边界单元法等数值方法来 进行求解,而离散单元法正是充分考虑到岩体结构的不连续性,适用于解决节 理岩体力学问题，离散单元法除了用于边坡、采场和巷道的稳定性研究以及 颗粒介质微观结构的分析外,已扩展到用于研究地震、爆炸等动力过程和地下 水渗流、热传导等物理过程，
随着计算机技术的发展,科学技术的进步,科学与工程计算 简称科学计 算 的应用范围已扩大到许多学科领域,形成一些边缘学科,例如计算物理、 计算化学、计算力学等，目前,实验、理论、计算已成为人类进行科学活动 的三大方法，
为了解某科学与工程实际问题,首先是依据物理、力学规律建立问题的 数学模型,这些模型一般为代数方程、微分方程等，科学计算的一个重要方 面就是要研究解这些数学问题的数值计算方法 适合计算机计算的计算方 法 ,然后通过计算软件在计算机上计算出实际需要的结果，数值分析内容 包括；函数的插值与逼近方法,微分与积分计算方法,线性方程组与非线性 方程组计算方法,常微分与偏微分数值解等，
岩土工程中最基本的两种分析方法是总应力分析 Total stress Ana1ysis 法和有效应力分析 Effective stress Analysis 法，与此相应,岩土工程中的有 限元法也可分为总应力分析有限元法和有效应力分析有限元法，l 966 年,cloughu等首先将总应力分析有限元法用于土坝的应力和变形分析，1969 年,sandhu和wilson用有限元法分析了Biot二维固结问题,开创了岩土工程有 效应力分析有限元法的先河，在国内,沈珠江 1977 首先将有效应力分析有 限元法应用于软土地基的固结变形分析，
香港大学佑启开展了结构力学有限元法的卓越研究 例如有限条分法等, 同时他在弹性地基上的研究,事实上奠定了边界元法的基本理论，
1、基本 概 念
数值分析的主要求解 方 法
变分法是讨论泛函的极值问题,对上述差分法及有限元法都可起推导基 本公式的作用,而这方法本身,也是数值方法中最古老的方法，
加权余量法,可以引入试函数和权函数的方法,从微分方程中直接求出近 似的数值解，它的优点是可以避免建立能量方程,使一些无法求得能量方程 的课题,也得到了较精确的解答，
Lemos于1983年开发了离散单元法与边界单元法耦合的半平面程序,并 用于计算节理和断裂介质中的应力分布问题，
2、几种常见数值分析方法的主要 特 点
2.2 离散单元法Discrete Element Method DEM
Lorig于1984年开发了包括前处理和后处理的离散单元法与边界单元法 耦合程序， 翌年他到澳大利亚英联邦科学与工业发展组织的岩土力学研究 所,修改了他原先的程序，这个修改后的程序文本称为HYDEBE HYBRID DISCRETE ELEMENT BOUNDARY ELEMENT ,其功能更强,包括一个前 处理程序CREATE,类似于有限单元法程序中的自动划分网格,一个与边界单 元法耦合程序BOUND和一个离散单元法程序BLOCK，
数值分析在岩土工程中的应用
Add the author and the accompanying title
数值分析方法 在岩土工程中的应用简介
A Brief Introduction to the Application of Numerical Analysis Method in the Geotechnical Engineering
差分法由英国学者Southwell提出,是转化课题的常微分方程或偏微分方 程为差分方程,然后结合初始及边界条件,求解线性代数方程组，这方法比 校直观,容易编制程序,所以从20世纪40年代末期以来,迄今盛行不衰，
有限元法由英国学者Zienkiewicz提出,以研究区域的变分形式,离散化所 研究的区域,使其成为有限数目的单元，对于边界复杂以及材料属非线性的 课题,与差分法相比,有特殊的灵活性，
2005年度辽宁省建设厅执业资格注册中心 注册土木工程师 岩土 继续教育培训专用
陈剑 平
吉林大学朝阳校区建设工程学 院
原长春地质学院水文地质与工程地质系 或长春科技大学环境与建设工程学院
2005年2月25日
讲授提纲
1、基本概 念 2、几种常见数值分析方法的 主要特点 3、对数值分析定量结果的理解
1、基本概 念
2、几种常见数值分析方法的主要 特点 2.2 离散单元Discrete Element Method DEM
一般从宏观意义上说,岩石可以视为连续介质,从而可以用弹性力学或塑 性力学的方法来进行分析和计算，但在某些情况下,岩体却不能视做连续介 质,如地下节理岩体中的巷道 见下图 ,这时,就不宜用处理连续介质的力学方 法来进行计算，于是,离散单元法作为一种处理节理岩体的数值方法就应运 而生，
2、几种常见数值分析方法的主要 特 点
2.2 离散单元法Discrete Element Method DEM
至于三维离散单元法的发展则要迟些,共主要原因是数据结构复杂,要 求计算机应具有较大的容量,并且计算结果的图形显示较为困难,如切一剖 面,块体间一般来讲是不接触的,犹如浮在空中,看起来不直观，三维离散单 元法程序3DEC 3—DIMENSIONAL DISTINCT ELEMENT CODE 已由Cundall与ITASCA咨询集团于1986年开发出来，其基本原理同 UDEC一样,只是数据结构作了较大的改进，
1、基本概 念
数值分析的必要 性
应用力学领域内的很多课题,工程师和科学工作者近年来一般用数值方 法来求解,获得显著成功，这主要由于目前力学课题本身的复杂性——非均 质、非线性以及复杂的加荷条件及边界条件,精确解已无能为力，
数值分析的可能性 计算机的迅速发展,也使数值分析得到有效而经济的成果，
1、基本概 念
离散单元法的基本思想 离教单元法也像有限单元法那样,将区域划分成单元 见图1—3 a ，但
是,单元因受节理等不连续面控制,在以后的运动过程中,单元结点可以分离, 即一个单元与其邻近单元可以接触,也可以分开 见图I—3 b ， 单元之间相 互作用的力可以根据力和位移的关系求出,而个别单元的运动则完全根据该 单元所受的不平衡力和不平衡力矩的要求解方法
在量测手段得到改进的今天,少数课题利用实验室的结构模型试验或离 心模型试验也可以得到定性或甚至定量的结果，这样,上述力学求解方法可 列成下表：
力学求解方法
精确解 数值方法
实验手段
差分法 有限元法
边界元法 变分法
加权余量法
1、基本概 念
数值分析的主要求解方法
图中数值方法列出最常用的五种：差分法、有限无法、边界元法、变 分法和加权余量法，
2、几种常见数值分析方法的主要 特点
2.1 有限元Finite Element Method FEM 对于土体,总应力分析法是将其视为固体来分析的，因此岩土工程中的
总应力分析有限元法与一般固体力学中的有限无法是相同的，有效应力分 析法则严格区分土体中分别由土颗粒骨架、孔隙水和孔隙气传递或承受的 应力,并考虑土骨架变形、孔隙水压力消散和孔隙气压力消散三者的耦合作 用,因而比总应力分析法更接近实际,但远较之复杂，对于岩体,因其本身可 视为固体,故与之相应的分析只能是总应力分析，从有效应力原理出发,可将 总应力分析视为有效应力分析的一种特殊形式，因此,总应力有限元分析法 也可视为有效应力分析有限元法的一种持殊形式，
正确的剖分
不正确的剖分
2、几种常见数值分析方法的主要 特点
2.1 有限元Finite Element Method FEM
研究表明,对于大多数岩土工程问题,无论是进行总应力分析还是进行有 效应力分析,均可归结为对Biot动力固结方程的求解，因此,可将Biot动力固 结方程作为岩土工程问题的总控制方程，
边界元法由英国学者Brebbia提出,化微分方程为边界积分方程,使用类似 于有限元法的离散技术来离散边界，离散化所引起的误差仅来源于边界,因 之提高了计算精度,依靠边界节点上算得的量,即可计算区域内的有关物理 量,从而减少了准备工作量及计算量，边界元又有直接法及间接法之分,间 接法需要先求取一个虚设的量,多了一道手续,所以相对麻烦一些，边界元 的缺点是对变系数或非线性问题的造就性不如有限元法，
2、几种常见数值分析方法的主要 特 点
2.2 离散单元法Discrete Element Method DEM
离散单元法在国内外的发展状况 离散单元法一般认为是Cundall于1971年提出来的，该法适用于研究在
准静力或动力条件下的节理系统或块体集合的力学问题,最初用来分析岩 石边坡的运动，到1974年,二维的离散单元法程序趋于成熟,当时已有屏幕 图形输出的交互会话功能，但由于受计算机内存的限制,不少程序是用汇 编语言写成,到1978年才全部翻译成FORTRAN IV的文本,成为离散单元法 的基本程序，与此同时,Cundall和 Strack还开发了二维圆形块体的BALL程 序,用于研究颗粒介质的力学行为,所得结果与Drescher等人用光弹技术的实 验结果极为吻合,使BALL程序在研究颗粒介质的本构方程方面大放异彩，
Cundall于1980年就开始研究块体在受力后变形以及根据破坏堆则允许 断裂的离散单元法, 这显然是将块体视为刚体的离散单元法的一个进步， Cundall称这种方法为UDEC UNIVERSAL DISTINCT ELEMENT CODE ,该 程序最后于1985年完成，UDEC现已广泛用于岩土力学和采矿工程,被公认 为对节理岩体进行数 值模拟的一种行之有效的方法，
2、几种常见数值分析方法的主要 特点
2.1 有限元Finite Element Method FEM
在众多的工程数值计算方法中,有限单元法 Finite Element Method,缩写 为FEM 很早就以其适用性强和处理非均质、非线性、复杂边界诸多问题方 便等突出优点而成为工程数值分析最有效的通用工具，经过近半个世纪的 发展,有限单元法已十分成熟并在各个领域的工程分析中广泛应用，
2、几种常见数值分析方法的主要 特点
2.2 离散单元法Discrete Element Method DEM
离散单元法是一种显式求解的数值方法，该方法与在时域中进行的其他 显式计算相似,例如与解抛物线型偏微分方程的显式差分格式相似， 显式是 针对一个物理系统进行数值计算时所用的代数方程式的性质而言，在用显式 法计算时,所有方程式一侧的量都是已知的,而另一侧的量只要用简单的代入 法就可求得，这与隐式法不同,隐式法必须求解联立方程组，在用显式法时, 假定在每一迭代时步内,每个块体单元仅对其相邻的块体单元产生力的影响, 这样,时步就需要取得足够小,以使显式法稳定，由于用显式法时不需要形成 矩阵,因此可以考虑大的位移和非线性,而不必花费额外的计算时间，
离散单元法与其他数值方法 如有限单元法、边界单元法等 耦合更能发挥 各自方法的优点，例如,用边界单元法考虑远场应力的影响以模拟弹性的性质, 用有限单元法作为中间过渡考虑塑性变形,再用离散单元法考虑近场不连续变 形的情况,从而极大地扩展了数值方法的解题范围，
2、几种常见数值分析方法的主要 特 点
2.2 离散单元法Discrete Element Method DEM
岩土工程中的基本方程包括土体平衡 或运动 方程、物理 或本构 方程、 几何方程、有效应力原理、孔隙流体 水 平衡方程、连续方程等，总控制方 程 即Biot动力固结方程 由这些基本方程组合而成，
总控制方程的推导基于以下假定： 1 土体是完全饱和的横观各向同性弹性体， 2 土体的变形是微小的， 3 土颗粒和孔隙水不可压缩， 4 孔隙水相对于土骨架的渗流运动服从Darcy定律,其惯性力可不计， 5 应力应变的正负号法则与弹性力学相反，
可行的
2、几种常见数值分析方法的主要 特点
2.2 离散单元Discrete Element Method DEM
近几十年来,离散单元法有了长足的发展,已成为解决岩土力学问题的一个 重要的数值方法,越来越受到人们的重视，因为工程中所见到的岩体其形态常 呈非连续结构,所形成的岩石块体运动和受力情况多是几何或材料非线性问题, 所以很难用解决连续介质力学问题的有限单元法或边界单元法等数值方法来 进行求解,而离散单元法正是充分考虑到岩体结构的不连续性,适用于解决节 理岩体力学问题，离散单元法除了用于边坡、采场和巷道的稳定性研究以及 颗粒介质微观结构的分析外,已扩展到用于研究地震、爆炸等动力过程和地下 水渗流、热传导等物理过程，