摘 要：在桩基设计中，地基土水平抗力系数之比例系数m是一个极为重要的参数。目前在工程中，对系数m的取值主要借助于经验公式、查规范表和单桩水平静载试验。本文提出一种新的方法—数值试验法，在一模型中运用该方法成功地求出已知地基土弹性模量E下的水平抗力系数之比例系数m。基于该模型，我们进行单因素数值试验来探究地基土水平抗力系数之比例系数m与其弹性模量E之间关系。通过对试验数据进行回归分析得出在单因素数值试验中，地基土水平抗力系数m与其弹性模量E呈线性正相关的联系。最后，提出正交试验来从整体上研究桩刚度、地质条件等因素对地基土水平抗力系数之比例系数m的影响程度。
　　关键词：地基土水平抗力系数之比例系数；弹性模量；数值试验；回归分析
　　0 引 言
　　桩基础是应用广泛的一种基础类型^[1]，除了承受竖向荷载还往往受水平荷载的作用。鉴于横向受荷桩承载力的重要性，我国建筑、路桥、铁道等部门的有关规范中对横向受荷桩的计算与设计要求均有不同程度的反映。由于地基土横向抗力系数的选取直接影响着横向荷载下桩基的位移和内力计算结果，因此如何选择并确定该系数，使之既安全可靠又合理适用是非常重要的^[2]。
　　在桩基设计中地基土水平抗力系数之比例系数m是一个较为关键的参数，目前工程中对参数m的取值主要由经验公式、查规范表和单桩水平静载试验来确定^[3-4]。在工程实际运用中却发现：经验公式所需参数单桩水平临界荷载H_cr、单桩水平临界荷载对应的位移x_cr均需做单桩水平静载试验并对试验数据进行一定的处理后方可得到。由于载荷试验费用高，在一般工程设计中，多数利用规范查表确定m值的方法^[5]。确定m值的规范法，由于每类土的m值上下限相差很大，如果没有其他土质性能指标可供参考，势必在运用时造成很大的任意性，从而影响计算结果^[6]。同时，根据地质勘查报告等资料，我们易获取土体变形模量、泊松比、内摩擦角、粘聚力、密度等参数，可据此反算出参数m。
　　基于上述问题，本文提出利用“m”法及大型通用有限元分析软件ANSYS做数值试验的方法来确定地基土水平抗力系数的比例系数m值。
　　1 m法简介
　　在线弹性地基抗力法中，地基抗力系数C如何随深度z变化，国内外学者对此有多观点，比较适用的一种^[2]： m法——假定地基系数C随深度呈正比例增加，即C=mz，如图1所示。 　　“m”法的基本前提是Winkel假定，即线弹性地基抗力法，它是用一组独立的离散型弹簧代替支挡结构开挖侧土体对支挡结构的抗力作用。该法克服了张氏法k为常数的缺点，能较好地反映地基系数沿深度分布的情况。此时桩的微分方程式为： 　　m法是迄今为止世界上运用最为广泛的计算方法，在英国、美国、前苏联、中国都把该方法列入规范之中，其比例系数m值可通过侧向荷载试验直接获得。
　　m法仅能反映土的弹性性能，在桩身位移较小时，能很好地反映桩土相互作用。当水平荷载较大且桩在泥面处发生较大位移时，桩侧土体进入塑性工作状态，用m法计算则出现较大误差，并且随着桩顶位移的增大，误差也随之较多的增大。m法在桩身位移较小的情况下比较适宜，而在实际工程中结构的位移均比较小，因此该方法也在建筑工程和公路桥涵的基桩基础设计中被广泛使用。
　　2 研究思路
　　根据《建筑桩基技术规范》^[3]第5.7.5条第2款，地基土水平抗力系数的比例系数m，宜通过单桩水平静载试验来确定，我们可用数值试验来模拟进行。基于ANSYS平台可计算出给定荷载下的桩顶水平位移，利用数值法和“m”法算出的桩顶水平位移相等，通过构造方程进行数值分析可反算出参数m。
　　3 数值实验分析
　　地基土水平抗力系数的比例系数m（MN/m⁴）（为地面以下2(d+1) m深度内各土层的综合值^[3]。通过ANSYS进行数值模拟求得桩顶水平位移，根据《建筑边坡工程技术规范》[7]查表可大致确定土体水平抗力系数的比例系数m的分布区间，借鉴《数值分析》二分法及迭代法^[8]原理，通过“m”法法来计算桩顶水平位移，直至所选土体水平抗力系数的比例系数m达到一定精度。
　　3.1 物理模型
　　选取钢筋混凝土圆截面桩作为研究对象，研究土体弹性模量E与其水平抗力系数的比例系数m的关系。
　　（1）几何模型：桩长23.5m，直径d=1.6m，桩的自由段长为1m，嵌岩段长3m，如图2所示。土体、基岩平面尺寸在桩周各外扩2倍桩长，桩底基岩向下延伸1倍桩长。
　　（2）材料特性及参数
　　本模型是以贵州某一水平受荷嵌岩桩为工程背景，场地各岩土层自上而下为：红粘土，中风化石灰岩。桩为线弹性材料；土体(可塑)为E-µ模型,屈服准则是DP破坏准则；基岩为DP材料。通过查阅大量的工程地质资料及工程地质手册[9]，桩基各参数选取见表1。
　　（3）有限元模型
　　有限元计算采用大型通用有限元软件ANSYS,土体和基岩单元为SOLID45，桩单元为SOLID65。桩身本构关系采用线弹性，而土体本构采用非线性E-µ模型[10]。本算例采用的计算模型如图3所示，该模型的4个侧面及底部边界均被固定约束，在桩顶施加100kN的水平荷载。 　　3.2 确定地基土水平抗力系数的比例系数m
　　这里以E=19MPa时地基土水平抗力系数之比例系数m的确定为例来展示推导过程。在本模型中，我们基于ANSYS软件平台编辑APDL程序来计算，在ANSYS后处理进程中可查询出桩顶水平位移。
　　基于此，我们构造出参数m的方程： 　　令
　　这里: m—土体水平抗力系数的比例系数(MN/m⁴);
　　—“m”法计算得出的桩顶位移（mm）。
　　结合土体物理力学参数等信息，通过查阅《建筑边坡工程技术规范》表G.0.1-2 土质地基系数^[7]可初步确定土体水平抗力系数的比例系数MN/m⁴。
　　由桩基传统解法可计算出
由此说明方程在区间是有解的，我们选取的区间