有限元法在土力学中的应用.docVIP

有限元法在土力学中的应用 在很多岩土工程的实际问题中，例如档土墙、板桩、基础梁和板等工程，由于岩土的非均质、非线性的性状以及几何形状的任意性、不连续性等因素，在多数情况下不能获得解析解。最近二十多年来，随着电子计算机的迅速兴起，在岩土工程中，数值分析受到了极大的重视，各种数值方法在岩土工程中都得到了广泛地应用，而岩土工程中的各种复杂问题的解决又深化和丰富了数值分析的内容。 目前．在岩土工程的数值分析中，用的最为普遍的是有限元法和差分法，其他方法如边界元法正在兴起。变分法与加权余量法既可以独立地作为数值方法运用于土工实际问题的求解，又可作为推导前几种数值方法的手段。当数值分析中的差分法首先盛行于工程科学时，土工中的渗流及固结问题在四十年代后期也开始采用差分法成功地解决了某些实际问题，如土坝渗流及浸润线的求法、土坝及地基的固结等。五十年代及六十年代初，弹性地基上的梁与板以及板桩也用 差分法来求解。六十年代，土石坝的静力问题用有限元法来求解。由于有限元解法的灵活性，使差分法在土工中的应用暂时趋丁停滞。进入七十年代之后，土石坝及高楼(包括地基)成功地使用有限无法解决了抗震分析。七十午代后期及八十年代，边界元法异军突起。这方法特别适宜于半无限域课题，这些是土力学及地基工程学科经常遇到的边界情况。近十年来，地基的静力及动力问题，例如桩基及强夯(即动力固结)等，都使用边界元法得到了有效地解决。 有限元法是一种十分有效的数值分析方法。它有几个突出的优点：(1)可以用于解非线性]问题，(2)易于处理非均质材料，各向异性材料，(3)能适应各种复杂的边界条件。岩土材料恰恰存在这几方面的问题，因此很适宜采用有限元法。有限元法刚刚发展起来，就引起了岩土力学界的浓厚兴趣。1966电美国克技夫(C1ough)和伍德沃德(Woodward)首先将有限无法应用于土力学，作了土坝的非线性分析。接着大批发土力学工作者从事这方面的研究取得巨大进展。在国内也是这样，1973年河海大学和南京水科院开始作有限元法用于岩土工程的研究，接着发表了我国最早的有关论文。目前国内大型土石坝的设计已普遍应用有限元法，其它岩土工程问题对有限无法的应用也得到了推广。 有限元法应用于岩土工程，显示了强大的生命力。许多过去无法解决的问题用有限元法得到了很好的解决。举例来说，某大坝的围堰，高90m，一种设计方案是在70m深的水中抛石填砂，然后筑两道1米厚的混凝土防渗墙挡水，断面如图1所示。要通过计算分析回答： （1）两道防渗墙在86．5m高的水头作用下受力和变形情况如何？能否承受得了? （2）防渗墙采用低标号的沤凝上有利，还是高标号有利? （3）水中填土是否要部分压实加固?在什么部位压实有利? (4)两墙之间的水位控制多高为宜? 这些问题是无法用通常的计算方法进行分析的。有限元法却能较好地加以解决。此外，按比奥别论解土体固结，在动力作用下解土体的变形与稳定，土体的流变，以及裂隙岩体的应力和变形等复杂问题，也都可以用有限元法。 图1 有限元法简介 有限无法是将一个连续体结构，如图2所示土坝，离散成有限个单元体，如图3，这些单元体在结点处互相铰结，把荷载简化到结点上，计算在外荷裁作用下各结点的位移，进而计算各单元的应力和应变。用离散体的解答近似地代替原连续体解答。当单元划分得足够密时，它与真实解是接近的。 图2 图3 第二节 非线性分析方法 非线性问题有两种：一是材料非线性，即材料的应力－应变关系是非线性的；一是几何非线性，即存在大变形，应变与位移的关系不呈线性，应变不仅包含位移对坐标的一阶导数，还要包含高阶导数。不管哪种非线性，推出的劲度矩阵都将随位移而变。 这是位移的非线性方程组。直接解这样的方程组是因难的，因此将其他为一系列的线性问题，用近似的方法求解。本文以材料非线性为例介绍几种非线性近似解法。 一、 迭代法 迭代法是将荷裁一次施加于结构，不断地修正劲度或调整荷裁， 来逐步接近真实解，而每次迭代作了一次线性有限元计算。迭代法又可分为割线迭代，余量迭代，初应力迭代等。现分述如下： (一)割线迭代 离散结构上作用有荷载{R}，相应地产生位移仍{}，如。果材料是线性的，{}随{R}成正比变化；如果材料是非线性的，{}－{}关系如图4(a)所示，那么{R}－{}关系也是非线性的，见图4 (b)。图中{}和{}的坐标是抽象的广义的，包含了应力各分量与应变各分量之间的关系。其“斜率”可以抽象地理解为[D]；{R}和{}的坐标也是抽象的，其“斜率”为劲度矩阵[K]。{}和{}关系是在试验基础上由本构关系模型给定的， 图4 (二)余量迭代法 余量迭代法是先将总荷载施加于