传统桩基础设计方法通常假定上部荷载全部由桩来承担,不考虑桩间土对基础承载力的贡献。在过去几十年的工程实践中,越来越多的设计人员认识到传统桩基础设计方法在某些情况下过于保守,由于筏板与地基土客观存在的接触作用,筏板能分担相当一部分上部荷载[1-3]。因此,如果在桩筏基础设计中考虑筏板的荷载分担作用,除了能满足建筑物对基础安全性和变形的要求外,还能减少较为可观的桩数,降低基础工程造价。
早在20世纪30年代,工程设计人员根据简单的共同作用原理,采用地基土分担“老8吨”(即80 kPa)的规定对上海高层建筑进行桩基设计,相关建筑至今仍能满足正常使用要求。这种设计方法虽然较为简单,未考虑基础形式对筏板荷载分担比的影响,但其充分体现了桩-筏共同作用的设计理念。基于桩-筏共同作用设计方法的核心之一是合理确定筏板荷载分担比的大小。相关学者基于工程经验、理论分析及数值模拟等手段对桩筏基础筏板荷载分担问题进行了研究,并提出了相应的计算方法[4-11]。但需要指出的是,总体来看,这些计算方法计算参数的确定较为困难,经验性较强,不便于在实际工程设计中广泛应用。
鉴于此,笔者基于刚性板桩筏基础共同作用计算模型,通过分析一系列不同土体参数及基础几何参数下的桩筏荷载分担比变化规律,提出了考虑基础整体安全系数影响的归一化筏板荷载分担比分析模型,并通过10个实际工程案例验证了归一化模型的准确性。
表 1为相关文献中给出的桩筏荷载分担比计算方法。从表中可以看出,大部分桩、筏荷载分担比计算方法都需要较多参数,如方法1、2、4,而这些参数的确定又较为复杂,存在较大的不确定性,因此,在实际工程应用中存在诸多不便。
此外,目前的桩、筏荷载分担比计算方法多为“确定性”计算方法,无法考虑荷载分担比在加载全过程中动态变化的特点。方法1、4考虑了桩基荷载分担比随沉降的变化过程,但由于桩筏基础的沉降预测是一个复杂的课题,在设计阶段往往难以准确预测出基础沉降值的大小。因此,在实际设计中难以准确计算在工作荷载作用下桩、筏荷载分担比的大小。
考虑到基础整体安全系数Fs较基础变形更容易在设计阶段确定,因此,以工程设计中常用的安全系数Fs作为变量,建立起筏板归一化荷载分担比分析模型,更便于实际工程的设计应用。
图 1为桩筏基础共同作用分析模型,为简化计算过程,提高计算效率,筏板假定为刚性筏板。在刚性筏板假定的基础上,筏板下各桩头具有相同的沉降值。对于实际工程而言,由于筏板自身、群桩和上部结构对基础体系刚度的贡献很大,采用刚性筏板基础也能较好地反映桩筏基础的荷载分担特性[12-14]。
刚性板桩筏基础共同作用模型的分析:
1) 土体位移方程。对筏板和桩基进行计算单元剖分,假定桩侧和桩端单元界面的剪应力和压应力为均匀分布,通过对桩侧和桩端面积进行数值积分,可建立反映桩-土-筏共同作用的土反力与土体位移方程
式中:ws和Ps分别为土体竖向位移和土体反力列向量,其展开式为
式中:k和n分别为筏板和桩身节点数。
Fs为土体柔度系数矩阵,矩阵各元素反映桩、土、筏之间的相互作用
式中:Frr为筏-土相互作用柔度系数矩阵,主要包括筏-土柔度系数(主对角线元素)Fi, i和筏-土-筏相互作用柔度系数Fi, j。筏-土柔度系数Fi, i按式(4)计算。
式中:Hsoil为可压缩土层深度,该深度以下为不可压缩土层;requ为筏板单元等效半径(按矩形筏板单元与圆形筏板单元面积相等原则计算得到);Gs为土体剪切模量。
Frr中非对角线上元素筏-土-筏柔度系数Fi, j按式(5)计算。
式中:ri.j为筏板单元i和单元j中心点之间的水平距离;E为土体弹性模量。
Fpp为桩-土相互作用柔度系数矩阵;Frp和Fpr为桩-土-筏相互作用柔度系数矩阵。各柔度系数均可按照Mindlin弹性理论计算,具体计算方法可参考文献[12-14]。
2) 桩身位移方程。考虑桩身单元划分的灵活性,根据杆系有限元理论,将桩看作通过一系列结点连接的杆单元,通过联立各单元刚度矩阵得到基础整体刚度矩阵Kp,进而得到桩身位移与荷载的关系式,即:
式中:wp为节点竖向位移列向量;Qtop为基础外荷载列向量
3) 桩土界面位移协调。假定桩-土界面无相对位移产生,即wp=ws,并对土体柔度系数矩阵Fs进行求逆得到地基刚度矩阵Ks。将式(1)代入式(6),可建立桩筏基础体系外部作用荷载与基础位移的数学关系式,即
式中:Ktot为桩筏整体刚度矩阵,为群桩竖向刚度矩阵Kp和地基土刚度矩阵Ks之和。
4) 在模型中考虑土体塑性影响,令桩、筏各节点极限抗力为
式中:Su为土体不排水强度,可以基于原位试验或室内三轴试验加以确定;η为桩-土摩阻力系数,一般取值为0.35~1.25。在计算过程中,当节点反力大于极限反力时,用极限反力值代替该结点反力。
5) 为考虑桩筏基础的非线性沉降特性,土体模量E采用式(9)所示双曲线型变化函数。
式中:Es为土体的初始弹性模量;q、qu为节点当前节点反力值和极限抗力;Rf为非线性系数, 取值范围在0~1。
需要指出的是,计算过程中土体弹性模量E需要根据节点相应荷载值进行不断调整修正。计算得到各节点应力后,按照式(9)对土体实际弹性模量进行修正,然后按照修正后的弹性模量E计算地基刚度矩阵Ks,并代入整体荷载-位移方程(7)进行求解。
6) 采用增量法对式(7)进行求解。由于假定为刚性筏板,为了保证筏板的各节点的位移相同,需通过乘(或加)大数法对式(7)两侧元素修改后进行求解来实现,其具体处理方法可参考文献[12-14]。对于筏板荷载分担比的计算,根据式(7)可得到不同筏板沉降下桩身、筏板节点的位移wp;考虑桩-土界面位移协调条件wp=ws,基于方程(1)可得到各节点的土反力Ps,通过对各节点土反力进行求和计算,即可分别得到桩基、筏板各自所承担的荷载值,从而得到筏板的荷载分担比大小αr。
以Basile等[15]给出的经典计算案例为例,其桩筏基础布置及相应材料参数如图 2所示。筏板尺寸为6 m×10 m,筏板厚度为1 m,筏下布置9根桩径为0.5 m、桩长为20 m的桩基,单桩基础竖向抗压承载力为873 kN。
由图 3所示的荷载-沉降及荷载分担比变化曲线可以看出,计算结果与Basile等[15]给出的计算结果十分吻合,证明了采用刚性板桩筏共同作用分析模型可以准确地评估基础整体沉降及荷载分担特性。
采用上述介绍的刚性桩筏基础共同作用分析模型进行筏板荷载分担比的参数分析,其计算模型如图 4所示。
刚性桩筏基础桩基均匀布置,桩径为d,桩长为lp,桩间距为S,筏板边缘距离桩基外围尺寸为a=d。分析桩径、桩长、桩间距、桩数及土体参数对筏板荷载分担比的影响,共计算了5大类工况,相应计算工况如表 2所示。
图 5所示为在不同土体参数条件下筏板荷载分担比随基础整体安全系数Fs的变化规律。基础整体安全系数Fs的定义为
式中:N为桩数;Qp为单桩极限承载力;Qur为筏板极限承载力;Q为上部总荷载。
从图 5可以看出,随着基础整体安全系数Fs的减小,筏板荷载分担能力逐渐发挥,其相应的荷载分担比αr逐渐增大。总体来看,当基础整体安全系数Fs>4时,筏板荷载分担比的变化速率显著降低,筏板荷载分担比基本维持不变。这可能是因为上部总荷载Q相对较小,大部分桩侧与筏板节点应力仍处于弹性受荷状态,没有进入屈服阶段,桩基础与筏板同步承担上部荷载。此外,从图 5中还可以看出,土体参数对筏板荷载分担比基本没有影响,这主要与模型假定的均质土体条件有关,桩基及筏板的承载性能同步、成比例发挥。
图 6所示为在不同桩数条件下筏板荷载分担比随基础整体安全系数Fs的变化规律。与直观理解不同的是,在初始加载阶段,随着桩数的不断增加,筏板荷载分担比也逐渐增大。其原因是桩间距保持恒定,因此,单桩所分摊的筏板面积一致,而由于群桩效应的存在,桩数越多,筏板中间部位的桩基荷载发挥能力越弱,导致在基础整体安全系数较大时,筏板荷载分担比随着桩数的增加而不断增大。
图 7所示为在不同桩间距条件下筏板荷载分担比随基础整体安全系数Fs的变化规律。在此只介绍桩数为100时的变化规律,其他工况条件下的变化趋势与之完全相同。筏板的荷载分担比随着桩间距的增大而显著提高,因此,在桩筏基础设计中采用大桩距桩基布置方案不仅能有效发挥筏板荷载的承载潜力,还能减少桩基数量,节约建设成本。
不同桩长条件下筏板荷载分担比随基础整体安全系数Fs的变化规律如图 8所示。从图中可以看出,随着桩长的增加,桩基荷载承载力提高,筏板所承担的荷载比例逐渐降低,但降低的幅度明显减小。
桩径对桩筏荷载分担比的影响如图 9所示,从图中可以看出,随着桩径的增大,筏板荷载分担比也逐渐增大。这是因为在保持桩间距比S/d恒定情况下,桩径d越大,相应的桩间距离S也越大,单桩分摊的筏板面积也更大,单桩承载力的增长水平(基本与桩径d的一次方呈比例)要低于筏板承载力(与桩径d的二次方呈比例)的提高幅度。
经过大量拟合分析,提出了归一化筏板荷载分担比α'r概念。图 10所示为归一化筏板荷载分担比随基础整体安全系数变化曲线,其中归一化筏板荷载分担比α'r定义为
式中:αr为筏板实际荷载分担比;d0为标准桩径,d0=1 m;其余参数物理意义同前。
从上述归一化筏板荷载分担比α′r定义来看,桩间距S对筏板归一化荷载分担比的影响较为显著,桩长lp次之,桩数N的影响程度较弱。
从图 10中可以看出,引进归一化荷载分担比后,不同工况下筏板荷载分担比曲线能较好地得到了统一。相应的拟合函数如图 10所示。
为了验证所提出的归一化模型的有效性,搜集了10个位于不同区域的工程案例,场地包括软、硬黏土及粉土等,桩长12~60 m,桩径0.3~0.9 m,平均桩间距2.6d~9.7d,桩数20~351,基本涵盖了工程上常见的桩筏基础尺寸,具体工程实例资料如表 3所示。其中,若相应文献中没有明确给出建筑物安全系数Fs,则根据场地地质条件、荷载大小等信息,通过式(10)计算得到。
从表 3中可以看出:采用传统桩基础设计方法设计的建筑物整体安全系数大于5;相反,采用减沉桩设计的建筑物,由于本身已考虑了筏板的荷载分担作用,其整体安全系数在2~3之间,实现了基础设计安全高效、经济合理的目的。
图 11为工程实测归一化荷载分担比分布图,从图中可以看出,实测工程案例的归一化荷载分布规律与计算得到的分布模型吻合较好,尤其是当基础整体安全系数Fs≥5时,验证了归一化分析模型的适用性。
对于实际桩筏基础设计,如图 11所示,总体而言,当整体安全系数控制在Fs>4时,归一化筏板荷载分担比随受基础整体安全系数影响基本可以忽略,此时归一化筏板荷载分担比基本为一常数,可以近似按照归一化荷载分担比α′r=10进行筏板荷载分担比计算,即筏板荷载分担比
需要说明的是,该归一化荷载分担比计算模型主要适用于均质地基或吉布森地基,对于场地表面具有较厚硬壳层或者典型端承桩的情况,模型不适用,可能会带来较大误差。
基于桩筏基础共同作用模型分析了影响桩筏基础筏板荷载分担比的各因素,并根据计算结果提出了均质地层中归一化筏板荷载分担比计算模型,根据实际工程案例验证了归一化模型的有效性,得到以下结论:
1) 总体而言,筏板荷载分担比随着基础整体安全系数Fs的增大而逐渐减小,但当基础整体安全系数Fs>4时,筏板荷载分担比基本保持不变。
2) 桩数、桩长、桩间距直接影响归一化筏板荷载分担比的大小,但影响程度存在差别:桩间距影响较为显著,桩长次之,桩数影响程度较弱。
3) 对于实际桩筏基础设计,当整体安全系数基本控制在Fs>4时,可以按照归一化荷载分担比α′r=10进行筏板荷载分担比计算。