水平荷载作用下承台桩承载性状有限元分析.docVIP

水平荷载作用下承台桩承载性状有限元分析 　　【摘要】： 结合某具体工程，采用弹性地基反力法理论和杆系有限元模型对承台桩进行模拟。在水平荷载作用下进行有限元计算，验证了m所取值的可行性。并在此基础上，分析并判断承台桩在设计荷载下是否满足使用要求,探讨了m法在基础工程中适用性。 　　【关键词】：双排承台桩；有限元分析； m法； 桩土相互作用 　　中图分类号：F54 文献标识码：B 文章编号：1002-6908(2008)0810037-01 　　 　　引言 　　 　　某地的高速公路栈桥路段发生多次坍塌，路段桩基出现坏损，根据现场观测及分析，认为桩基还可以承受较大的水平力作用。 因此对塌桥过程中产生的水平力，采用有限元法进行数值模拟分析、判断桩的破坏程度[1-3]，为下一步指导施工提供依据。 　　 　　1. 桩地基水平抗力法 　　 　　桩基础中的桩，一般作为埋入地基中的弹性地基梁来进行分析，根据其在土中的受力状况由结构力学导出桩的绕曲基本微分方程式为： 　　 　　式中： E为桩材料的弹性模量，I为桩的截面惯性矩，z为桩在地面下任一点的深度，y为桩身深度z处一点的水平位移，N(z)为作用在桩身的轴向力，p为土的水平抗力。 　　桩在水平荷载及弯矩作用下与桩周地基土相互作用，桩的变形及内力取决于桩身所受到的外荷载类型及大小、桩周土性质及抗力以及桩顶约束条件等多种因素。桩周土的抗力与土体所受的状态、水平位移有关，可表示为： 　　 　　式中k为水平抗力系数；n为指数(一般n不大于1.0)；b为桩的计算宽度。 　　求解问题的关键是选择符合实际的土抗力函数p(z,y)，选择弹性地基反力法，把桩周土离散为一个个单独作用的弹簧。某一弹簧的受力时，仅与该弹簧发生于作用力成比例的压缩而和其他弹簧无关。为了分析的方便，使用弹性抗力法，取n=1。水平抗力系数k随深度的分布形式及大小对桩的内力与变形有直接的影响，假定k(z)=mz，地基反力系数沿深度按线性规律增大。 　　采用杆系有限元法[4]，水平承载桩的有限元分析以弹性地基梁绕曲微分方程为依据，其分析计算实质是个矩阵分析，将桩离散为若干单元，两单元段为一节点，外力作用于节点，可得到外节点力P同单元内节点力F的关系： 　　P=AF （A为比例常数）（3） 　　节点处变形e同外节点位移x之间的关系： 　　e=Bx （B为A矩阵的转置）（4） 　　对于所有单元，单元内力同内变形e的关系： 　　F=Se(S为单元刚度矩阵)（5） 　　有上式可得： 　　 P=ASATX （6） 　　求解得： 　　 X=（ASAT）-1P（7） 　　把x值代入可求得单元内力F。 　　 　　2. 模型的建立 　　 　　如图1所示，双排承台桩采用直径1.2m、长45m的灌注桩，桩中心距为3.6m，弹性模量为30MPa；与桩体相比，承台的整体刚度较大，采用杆系有限元进行计算时，通过增大其弹性模量来考虑这种因素，取承台的弹性模量为3000 MPa。 　　采用有限元计算时，取前排桩桩底为坐标原点，向右为x轴，向上为y轴。桩和承台均用beam3梁单元模拟。对桩体，单元长度取1m，共90个单元；对于承台，单元长度取1.8，共2个单元。桩周土用combin14单元模拟。 　　 　　3. 模型的计算与分析 　　 　　3.1 模型参数的选取 　　有限单元法求解的关键在于正确选择模型的设计参数[5]。在该地区m值在设计中一般取为5000，考虑到桩间土体的扰动等因素，对桩间土体水平向基层系数随深度的比例系数m乘以0.8的折减系数，m值取为4000(单位)。设计作用在承台上的水平荷载为kN, 距承台顶面为6.2m，进行简化，为作用在承台中点E的水平荷载kN及弯矩1277.2，采用弹性地基梁有限元法求解。 　　3.2 双排承台桩节点位移和弯矩计算结果 　　有限元计算的双排承台桩节点位移和弯矩计算结果见图2、图3所示。 　　 　　从图2、图3得出的计算结果可以看出： 　　1） 在水平荷载及弯矩1277.2作用下，节点最大位移为2.22毫米，发生在桩顶处，距离地面8米处，节点位移为零，距离地面8~16米处出现x方向正位移，距离地面16米至桩底，几乎没有位移。 　　2） 在水平荷载及弯矩1277.2作用下，单元最大负弯矩出现在AB桩的桩顶处，为-2.07E5，最大正弯矩也出现在AB桩上，距地面5m的位置，为1.89e5。 　　3） 荷载作用方向上前桩的内力比后桩稍大，但是各桩的桩身弯矩分布较均匀。 　　为了研究水平设计值是不是在桩的水平抗力之内，取m=4000，考虑水平力的不同对节点位移和内力的影响，计算结果如表1所示： 　　 　　由表1可知，当水平力在设计荷载附近变化时