混合法计算分析单桩承载力.docVIP

混合法计算分析单桩承载力0 引 言 Poulos（1968）认为桩被插入一个理想均质的弹性半空间内，桩身粗糙而桩底平滑，即桩身与桩周土之间具有位移协调条件，并在实践中得到了验证。Cooke (1974)提出摩擦桩的荷载传递模型，在发生剪切变形的桩周土中，桩周剪应力沿桩侧表面径向向四周土体传递。以此理论为基础，Y.K.chow等对群桩基础进行理论数值分析。下面提出的混合法－维杆单元模拟计算桩身轴力，结合cooke荷载传递模型分析桩周土，具有简单实用的特点。 1、理论基础 1.1 桩身模拟计算 桩身有限元分析中用一维杆单元模拟，这在桩的长径比大于5时计算精度是可以满足的，桩体模型见图1，长度为的桩单元刚度为，单元刚度矩阵为：。 1.2 桩周土分析 据cooke理论假定，桩周土中竖向主应力在上下土层间没有相互作用，只沿径向传递，在传递中衰减，大约20倍桩径外应变忽略不计，这在桩周土小变形下是适用的。结合弹性理论，桩周土体微单元的竖向微分极坐标平衡方程式： ，式中为桩周土距桩轴心距离处的剪应力，据cooke理论联合求解：，是土的剪切模量，为桩周土体距桩轴心距离处的竖向位移，为桩－土接触部位的剪切应力，为桩身半径，若土的不均匀系数为桩入土深度1/2处与桩底土剪切模量的比值，桩周土中剪切应力可以忽略的范围为，其中是桩的入土深度，是土的泊松比，饱和不排水时为0.5。 随桩周土变形的加大，桩周土的应力－应变关系曲线随土的类型与加载方式的不同发生变化，即加工硬化型与加工软化型。图2为土的加工硬化型曲线。 桩与土之间剪切应力与相对位移曲线可以用加工硬化型双曲线模型来拟合：，式中为位移的值，为位移的值，用标识，即桩周土的抗剪强度。 桩周土与桩体之间位移－剪应变之间的关系， 为经验系数，由桩身粗糙度决定，在0.002－0.001之间，土的加工硬化型曲线割线模量表达式：，：初始切线模量。一般情况下，桩土之间剪切破坏时剪应力达不到，两者之间的比值在0.7－0.95之间，通常取0.9，即：。 均质弹性半空间，剪切模量恒为常数，单由于土是一种典型非线性介质，在很小的变形下就表现为塑性特征，所以在数值计算中采用割线模量比较方便实用。此时刚度矩阵中的土刚度：，利用Bussinesq课题解答，桩端沉降：，采用非线性关系割线模量，则桩底土的刚度系数： ，式为桩底土极限荷载，为传递至桩底土的荷载。 2、非线性计算的实现方式 2.1计算矩阵合成 由桩侧土的刚度矩阵，桩的刚度矩阵，桩土体系总刚度矩阵：。有限单元法可以较好的模拟现实土体中的各种边界条件，本方法所需计算参数为常规地质报告所能提供，如由库伦理论得到，：静止土压力系数，：土重度，：计算点深度。 2.2计算迭代过程 计算过程采用增量法，利用土的割线模量分步加载。土体变形是非线性的，反映在土的刚度矩阵元素不再是常数，而是土体应变量的函数，整个刚度矩阵是非线性的。计算误差控制采用力的形式，在一定桩顶荷载、初始切线模量下求出位移 与桩周土各节点剪应力，利用式(7)求出桩周土各节点割线模量，重组刚度矩阵，在割线模量与位移下的反算桩顶力值，达到所要求的控制条件，为一无穷小量，如此反复，直至满足精度要求，图3为迭代过程示意图。 3、算例 场地地质条件：地表至8.0m范围内为均质砂土，8.0m以下为流－软塑粉质粘土，地表水位在地表以下8.0m左右，土的物理力学指标如下表。钻孔桩直径25cm，桩长4.5m，桩身弹性模量23Gpa，桩侧极限侧摩阻力，为深度，计算结果见图4。 深度（m） 土类 含水量 溶重 孔隙比 0.0~8.0 粉砂 25.3 19.1 0.70 内聚力 摩擦角 压缩系数 压缩模量 内聚力 4.0 35.5 0.116 19.0 4.0 3、结论 (1) 以杆单元模拟桩身，以cooke剪切传递模型分析桩周土的变形，Bussinesq用于桩底土变形计算，实测与计算比较结果令人满意。 (2) 本方法可用于分层土的计算，所需计算参数少、简单实用，可由常规地质报告提供的数据选用。 注：文章内所有公式及图表请用PDF形式查看。 1