计算机内插在土工分析中的应用.docVIP

计算机内插在土工分析中的应用 周逢君 （深圳市工勘岩土工程有限公司） 摘 要 本文对岩土工程的土工分析中常用的有限元法等进行了介绍，然后再对计算机内插方法，如拉格朗日内插法、不规则三角网（TIN）内插、标准格网内插法进行了算法研究并编写计算机程序，最后根据实际工程项目，应用软地基处理的监测点进行内插计算，对软地基的沉降和回弹问题进行分析。 关键词 土工分析、内插方法、软地基沉降 土工计算机分析的几种常用方法 从上世纪20年代起，数值分析方法开始逐步应用于土工分析，特别是计算机技术的发展，使得数值计算在土工分析中的应用更为方便和快捷，计算机软硬件的发展使得很多理论得以在实际工程中应用，有些数学领域的方法借助计算机技术广泛应用于土工分析中。 土工分析中常用的数值分析方法主要有有限差分法、有限元法、离散单元法、拉格朗日（Lagrange））i = 1,2,···n)，写成多项式为： 对大多数问题是求解该多项式的系数ci，将n个离散点的（1）式写成矩阵方程： 求解的系数Ci的方程式为： 对于二元五次拉格朗日多项式来说（1）式可写为： 2．2 计算机编程内插算法分析 对于（4）式，要求解的系数Cki个数为： (个) 实际工程应用中，往往样本点（m个观测值）多于所要求解的变量Cn(mn） 再求σ对系数Ckj的偏导数，并令其为0： 这样，线性代数方程组有21个方程和21个未知系数Ckj，求出Ckj后就确定了二元五次内插的表达式。 拉格朗日内插方法一般用于曲线和曲面拟合，在岩土工程中广泛应用于边坡稳定性分析中的滑裂面曲线拟合、基坑回弹分析，地下巷道底臌分析及软地基处理的沉降和回弹分析等。 不规则三角网（TIN）及标准格网内插及计算机编程在土工分析中的应用 3．1 TIN的基本理论 不规则三角网（Triangulation Irregular Network，简写为TIN）内插法是近几年，特别是计算机编程广泛应用于工程实践以来最为普遍的内插方法，它最早由俄国数学家Delaunay于1934年提出的三角形最小内角最大化准则，并证明了在没有四点或四点以上共圆条件下的平面离散点存在一种三角化方式，使得连成的三角网中的三角形满足三角形外接圆不包含第四个及以上的点的条件，该三角形最接近等边三角形，所以也通常称这类三角形为Delaunay三角形。 由离散点组成的三角网内插法一般用于空间内插，如表面插值等。其内插方式如图1，内插点P满足共面条件： P1（x1,y1,z1） 稰(x,y,z) P2(x2,y2,z2) P3(x3,y3,z3) 近 图1：三角形内插共面条件 A、B、C三个待求系数由三角形三个顶点的已知值求得，求出该三角形的三个系数后，该三角形内任一点（x,y值已知）内插的z值都可以用（7）式求得。 图1和（7）式为线性内插，它在计算上是最简单的内插方法，线性内插只能保证在边界处是连续的，但不能保证光滑，如图2。 在边界处不光滑 在边界处光滑 `法线方向 法线方向 法线方向 法线方向 图2：线性内插在边界处的非光滑性 图3：高幂次内插在边界处的光滑特性 高幂次多项式内插（以二元五次内插为例）的数学公式见（4）式，它类似于线性内插，认为表面模型是连续和光滑的（其二阶导数不为0），从数学上来说，由于五次多项式在内插时考虑到相邻三角形的几何特性，则垂直于每一个三角形的法线是连续变化的，相邻三角形边线上的法线无断裂，内插的表面是光滑的。 3．2 TIN和标准格网内插在土工分析中的应用 笔者应用以上各类内插理论结合计算机编程应用于岩土工程实际，取得了比较好的效果，在此以深圳软地基处理工程为例作一详细介绍。 我们承担了深圳西部地区一填海项目的软基处理工程，软基场地面积约1.5万平方米，全部为填海软地基。工程完成后按规范要求进行了4次地基沉降监测，监测期限为1