浅析岩土工程数值分析方法.docVIP

浅析岩土工程数值分析方法 　　摘要:岩土工程领域的岩土力学数值分析方法得到了迅速发展,出现了各种各样的数值分析方法。本文浅述了各种方法的研究现状和最新进展,着重介绍了各种方法的优缺点,并提出了解决问题的思路、方法和建议。 　　关键词:数值分析方法;岩土工程;连续变形;非连续变形 　　中图分类号:P62文献标识码:ADOI:10.3969/j.issn.1672-0407.2010.06.024 　　文章编号:1672-0407(2010)06-056-02收稿日期:2010―05―02 　　 　　岩土工程是一门综合应用岩石力学、土力学、工程地质学与基础工程学等基本知识和手段解决实际工程建设中有关岩体与土体变形及稳定问题的学科。近几十年来,随着大型水利水电工程、土木工程、铁路遂道工程等的迅速发展,作为岩土工程基础学科的岩土力学得到了长足的发展,而解决岩土力学问题的方法主要有实验方法、理论分析方法和数值模拟方法三大类。 　　 　　一、连续变形数值分析方法 　　 　　这类方法主要包括有限差分法、有限元法、边界元法、无限元法等,其中以有限元法应用最为广泛,这类方法主要分析岩土介质的连续小变形和小位移特性。 　　有限单元法是近似解法,单元剖分的疏密程度与质量、效益密切相关,在理论上如何把握好这个度且保证收敛是有待研究的课题。边界元法由于所考虑问题的维数降低一维,只须对研究区的边界离散化,因此输入数据少,节省计算机内存,解题花时少。又由于只对边界离散,离散化误差仅来源于边界,区域内有关物理量是用精确的解析公式计算的,故边界元法的计算精度较高,能直接对无限域或半无限域求解。拉格朗日元法运用流体力学中跟踪质点运动的物质描述方法,即拉格朗日拖带坐标系方法。这种方法避免了有限元法进行大型矩阵的复杂计算,但时间步长的选择成了一个新的突出问题,时步过大会导致解答的不稳定,时步太小则会使计算时间过长。 　　 　　二、非连续变形数值分析方法 　　 　　(一)基于块体理论的非连续变形分析方法 　　基于块体理论的非连续变形分析方法主要以离散块体系统为研究对象,针对岩土体的非连续与非均匀的特点,将岩土体视为完全非连续介质,对构成离散系统的各个子块的运动和变形进行数值分析。 　　1988年,石根华博士发表了博士学位论文“Discontinuous Deformation Analysis:ANew Numerical Model for the Static and Dynamics ofB1ock Systems ”,这标志着DDA方法的诞生。该方法引入了运动方程,用最小势能原理把块体之间的接触问题和块体本身的变形问题统一到矩阵的求解中,理论严密,精度较高,而且把静力和动力、正分析和反分析统一起来,不仅可以计算破坏前的小位移,也可以计算破坏后的大变形,对滑动、崩塌、爆炸和贯入等问题也十分有效,是一种不同于DEM的新的数值计算方法。 　　由于DDA方法在运动约束方面做得比较充分,理论上比较严密,且静力和动力分析采用了统一的数值计算格式。因此,在结构、岩体和土体的非连续大变形力学过程模拟方面发挥了较大潜力。 　　(二)界面单元有限元法 　　为了推广有限单元法,使其能够处理简单的非连续问题,人们提出了各种能够反映非连续变形性质的简单力学模型和特殊界面单元,并将其应用于单一性非连续界面力学行为的模拟中。这些模型主要包括联结单元、无厚度接触单元或节理单元、薄层单元和接触摩擦单元等。 　　界面单元有限元法将岩土介质视为准连续介质,仍以连续分析为主,但可以对个别具有控制作用的宏观非连续面的变形与破坏力学行为进行重点分析,在工程中得到了广泛应用。同时,该方法也存在诸多的局限性,如只能对原生的非连续界面进行计算,对次生的非连续界面无法进行处理,界面单元数目不能设置太多,界面弹簧刚度的选取较为困难等。因此,该方法在处理复杂的非连续变形问题时仍显得无能为力。 　　(三) 流形方法 　　有限单元法将岩体理想化为完全连续介质,离散单元法等则将岩体抽象为完全非连续介质即离散介质,这其实是走向了问题的两个极端。将连续变形与非连续变形统一起来的方法将更适于进行实际岩土体变形的分析。因此,近年来出现了另一种更新的数值方法――流形方法。 　　 流形方法是石根华博士应用流形的覆盖技术建立的一种把有限元法、非连续变形分析方法和解析方法包含在内的全新的统一计算方法。由于流形方法可在统一的理论框架下处理连续与非连续性变形问题而引起了许多学者的兴趣,成为目前计算岩土力学的热门课题。 　　流形方法的优点主要表现在它具有相对完善的非连续变形处理功能,可以在统一的数学理论框架下同时处理连续问题与非连续问题。流形方法较有限单元法更适于进行开裂模拟,但由于受网格连