摘要：考虑井阻和涂抹区的影响，建立成层未打穿砂井地基固结的数值模型和算法。将未打穿砂井底面下的软土层视为虚拟排水井，同时将未打穿砂井转化为打穿的多段砂井，使得土层边界条件和数值计算过程简单化。对超孔隙水压力用Lagrange插值多项式等结点插值近似，把砂井打设区固结方程转化为结点的超孔隙水压力的简单微分方程组，砂井固结方程转化为一线性代数方程组，以4层土为例，导出详细的计算公式。如增加土层数，仅需改变边界结点系数公式。最后，把该解法与相关近似方法的计算结果进行分析比较，并给出一个多层软土路基预压加固工程案例。结果表明，将下卧层简化为一维固结与考虑其三维固结对计算结果影响很小。

关键词：土力学 固结 成层地基 数值算法 未打穿砂井 

单位：郑州大学土木工程学院; 河南郑州450002