[建筑]盾构隧道推进时原水管沉降数值分析.pdfVIP

第三章 盾构隧道推进时原水管沉降数值分析 3.1 有限元基本理论 3.1.1有限元理论的发展状况 从 1965 年第一次出现“有限元”这个名词，到 70 年代美国第一个有限元结构分析系统 SAP 问世，到今天已有38 年的时间。整整一代人在这一方面做了艰苦卓绝的努力，已经形 成了若干国际上著名的大型计算力学系统，例如: MSC/NASTRAN, ANSYS, FENRIS, SESAM 和 MIDAS 等,解决了成千上万个工程实际课题，为科学技术的发展和工程应用做出 了巨大的贡献。 在国外，有限单元法在 50 年代中起源于航空工程中的飞机结构应力分析。60 年代初， 许多数学和力学工作者参加了有限单元法的研究，搞清了它的理论基础，使有限单元法得到 了很大的发展。在我国，60 年代初，著名数学家冯康教授和他的研究组提出一种以变分原 理为基础的三角形剖分的近似法，为偏微分方程求得了近似解，并在严密的数学基础上证明 了它的收敛性、稳定性和误差估计。当时还编制了程序，并成功的解决了刘家峡水坝的应力 分析和其它一些课题，这个方法就是人们现在所熟知的有限单元法。 60 年代有限单元法在土工中得到应用。由于它能较容易地处理分析域的复杂形状及边 界条件、材料的物理非线性和几何非线性性状，所以其应用发展非常快，在土工固结、稳定 和变形分析等各个领域得到广泛应用。该法被用于分析浅基础、桩基础及各类深基础、挡土 墙、堤坝、基础基坑和隧道等各类岩土工程问题的分析。不仅用于分析静力问题、动力问题、 抗震分析，还用于分析上部结构、基础和地基的共同作用。在有限元分析中可采用总应力分 析法，也可采用有效应力分析法。有限单元法已成为土工数值分析的主要手段之一。 有限单元法的主要优点是它有很强的适用性，应用范围相当广泛。它可以分析非均质 材料，各向异性材料和非线性应力应变关系的结构，也可以处理复杂边界条件等难题，此外， 有限单元法概念浅近易懂，容易掌掘。它可以在不同水平上去理解，既可以用直观的物理途 径来掌握这一方法，也可以建立严格的数学解释。 目前在对于有关地层变形分析研究中，仍然存在很多问题。解析法得到的是理论解， 精度高，计算量小，但解题范围有限，只有少数简单边界条件下可以得到解析解。对于复杂 的实际工程来说，利用解析法来求解必然要做很多简化处理，很多施上因素都不可能考虑进 去，并且这些简化处理会导致计算结果与实际情况相即甚远。因此本文通过有限元的数值模 拟方法对盾构隧道的施工中原水管的位移进行模拟[11] 。 3.1.2有限元的基本原理和基本步骤 有限单元法是一种数值方法，设想把连续体分割成有限个单元，它们在节点上相连接， 即以一个单元集合体来替代连续体，把作用在单元上的力等效地移到节点上每个单元选择一 个位移函数来表示位移分量的分布规律:按变分原理建立单元的节点力-节点位移的关系式； 然后根据节点平衡条件，把所有的单元关系集合形成一组以节点位移为未知量的代数方程 组，从而解得各节点位移。 有限单元法的数学基础为分割近似原理，例如求曲线围成的面积，我们可以采用分段 以直代曲的方法来求其近似值。分割越细，结果越精确。同样地，有限元法把连续体分割成 许多小单元，这些小单元可选用简单的位移分布函数。分割的网格越细，结果越精确。因此 有限单元法可以分析许多复杂的问题，只要分的网格有足够的细度，就能得到满意的精度。 - 22 - 对于不同物理性质和数学模型的问题，有限元求解法的基本步骤是相同的，只是具体 公式推导和运算求解不同。有限元求解问题的基本步骤通常为: 第一步：问题及求解域定义。根据实际问题近似确定求解域的物理性质和几何区域。 第二步：求解域离散化。将求解域近似为具有不同有限大小和形状且彼此相连的有限 个单元组成的离散域，习惯上称为有限元网格划分。显然单元越小，则离散域的近似程度越 好，计算结果也越精确，但计算量及计算误差都将增大，因此求解域的离散化是有限单元法 的核心技术之一。 第三步：确定状态变量及控制方法。一个具体的物理问题通常可以用一组包含问题状态 变量边界条件的微分方程表示，为适合有限元求解，通常将微分方程化为等价的泛函表示。 第四步：单元推导。对单元构造一个适合的近似解，即推导有限元的列式，其中包括选 择合理的单元坐标系，建立单元试函数，以某种方法给出单元各状