砂井地基固结的三维有限元分析及应用.doc

作者简介： 砂井地基固结的三维有限元分析及应用 刘国明，曾连灯 (福州大学土木工程学院，福建 350002) 摘要：本文引入一种新型排水板单元，将其结合到三维Biot固结有限元程序中，大大简化了砂井地基三维有限元模拟和计算工作量。同时，在对称性原理和Biot固结理论基础上，将三维非线性Biot固结有限元法与基于遗传算法的参数反演分析方法相结合，编制了能够进行砂井地基正演和反演计算的数值模拟程序，用于堆载预压下砂井路基固结变形的工程模拟。反演结果与实测吻合良好。 关键词：砂井地基；排水板单元；三维Biot有限元；遗传算法 3D FEM analysis of consolidation for sand drained foundations and its application LIU Guo-ming，ZENG Lian-deng (College of Civil Engineering，Fuzhou University，Fuzhou，Fujian 350002，China) Abstreat: A new element type, named PVD element, is introduced and combined with 3D Biot FEM program in this paper,which makes the calculation work of 3D FEA of sand drain foundation reduced enormously. Moreover, On the basis of symmetry principle and Biot’s consolidation theory, the forward and inversion program is worked out and applied to an engineering project to analyze the subsidence of a soft foundation preloaded by embankment. The results calculated by the program were very close to the measured results. Key words: Sand drain foundation; PVD element; 3D Biot FEM；Genetic Algorithm 1 前言 砂井排水法广泛应用于饱和软基的加固工程。作为连续介质的地基软土内穿插具有空间分布形态的砂井，严格地讲，其固结是典型的三维固结问题。但由于三维应力应变有限单元分析工作量很大且由于砂井数量极多，导致砂井地基三维离散模型不可接受的计算规模。在对砂井地基固结进行三维非线性有限元分析时，本文采用一种新的三维排水板单元法，此方法在不增加有限元结点的情况下，能够考虑塑料排水板或者砂井的作用。 反求地基参数是岩土工程数值模拟的重要问题，也是难题。对于高度非线性问题，常规优化方法搜索的最终结果为平稳点并不能保证收敛到全局最优点。遗传算法是一种别具特色的新的搜索优化方法，不受搜索空间限制性假设的制约，不必求诸如连续性、导数存在以及单峰等假设，因此相对于常规优化方法具有特有的特点。 2 排水板单元的有限元分析 2.1 排水板单元的构造 由于塑料排水板的截面尺寸(一般为l00mmx4mm)与周围土体尺寸(一般为1.2m～3.0m)以及自身的长度(一般为10m～25m)比较非常小，其主要作用是增加地基渗透性，而它的力学性能是次要的。因此考虑用线单元来模拟地基中的塑料排水板的主要作用。将排水板考虑为实体单元的网格划分如图1(a)所示，文中将排水板视为线单元，与实体单元共用结点A、B，如图1(b)所示，不考虑排水板的受剪、受扭性能，而只考虑其轴力以及对总体渗透性能的影响。此种思想类似于混凝土结构中考虑热传递的钢筋单元。图2所示的线单元即为排水板单元，不考虑该单元的自重。由于本文的三维单元采用8结点空间等参单元，考虑到单元之间的协调，排水板单元为2结点轴力杆单元。下面在通常的线单元基础上将其扩展为可参与固结计算。 (a)常规单元 (b)排水板单元 图1 塑料排水板单元的有限元划分 图2 排水板单元 Fig.1 PVD element’s grid of FEM Fig.2 PVD element 2.2 位移模式及力学劲度矩阵 由于该单元是一维线单元，