高承台锥形嵌岩灌注桩的屈曲稳定分析.doc

高承台锥形嵌岩灌注桩的屈曲稳定分析 摘要： 　　为探讨高承台锥形嵌岩灌注桩的屈曲稳定问题，先基于假定的桩身屈曲稳定分析模型和相应的桩身挠曲位移函数建立出桩土系统总势能的表达式，然后利用最小势能原理和变分法导出了桩身屈曲稳定特征方程，由此获得桩身屈曲临界荷载的能量法解答。对比分析表明，与普通的等截面桩相比，当桩材、桩身几何尺寸、入土深度和桩周土质等条件相同时，变截面桩屈曲稳定性更好。 　　 　　关键词：锥形桩；屈曲分析；能量法；桩身屈曲临界荷载 　　 　　1 前 言 　　 　　桩基础通常由于其承载力高、沉降量小在工程中得到广泛应用。特别是近年来，除了采用等直截面桩（如预制方桩、灌注圆桩等），由于在相同的土质条件下，变截面桩与普通桩相比，其单位体积承载力有明显提高，故变截面桩在工程实践中也得到了越来越多的应用。如柔性桩与桩周土一起组成柔性桩复合地基既提高柔性桩的面积置换率，同时也提高了柔性桩的承载力[1]。而对于桥梁工程中多采用的高承台嵌岩灌注桩，桩顶自由长度一般较大，当上覆土层较软弱时，桩身屈曲与后屈曲稳定问题成为难点[2-5]。 　　 　　因此，可考虑利用施工相对不难且承载和稳定性能更好的锥形变截面桩作为高承台桩基础的基桩。然而，类似锥形桩的变截面桩不仅桩型上不同于等截面桩，其承载机理也于等截面桩不尽相同。对于普通等截面桩的屈曲问题，国内外学者做了不少研究工作，也提出了一些方法[5-10]；而变截面桩的屈曲问题，则鲜有人问津，少有报导。因此，本文拟考虑将锥形截面桩用于桥梁工程中的高承台基桩，并对其进行屈曲稳定分析，将具有一定的工程与理论实际意义。 　　 　　2 锥形桩屈曲稳定分析 　　 　　2.1 计算模式 　　 　　考虑高承台桩基中的锥形嵌岩灌注基桩，其在桩顶竖向集中荷载作用下可简化为顶部自由、底部嵌固的弹性地基梁，如图1 所示。由此可建立桩土体系的总势能方程，其为桩身弯曲应变能、桩侧土的弹性变形能和外荷载势能之和，即根据桩端位移边界条件，假定桩身挠曲位移函数为其中，l 为桩长，h 为桩入土深度， EIz 为桩身抗弯刚度，Iz 为桩身截图1 锥形桩屈曲分析模型面惯性矩，对圆桩有Iz=π 4z d /64，dz 为沿桩轴向线性变化的截面直径，即d= a0z+ a1，a0 为桩径沿桩轴向的变化率，a1 为桩端直径。q 为桩侧土体所受水平抗力，根据m 法假定q=mb1(h-z)x，m 为桩侧土水平抗力的比例系数，b1 为桩身计算宽度，对圆形截面单桩b1=0.9(dz+1)，（dz≥1）。 　　 2.2 桩身屈曲临界荷载的求解 　　 　　由最小势能原理，对（1）式取变分，使δΠ = 0即δΠ/δ C = 0， 即其中，Cj（j=1,2,…,n）为待定参数，n 为所取挠度函数的半波数。取n≥16 可满足精度要求。由（3）式可得一齐次线性方程组，由方程组可以定出个n 可变参数Cj（j=1,2,…n），要使Cj 具有非零解，则方程组的系数行列式即桩屈曲稳定方程式D=0，即上式即为桩身屈曲稳定的特征方程，其根即为桩屈曲临界荷载Pcr。此方程略做变换后即可以采用householder 镜像映射法求解。其中系数( , 1,2,..., ; ) ik ik ki i k n λ λ λ = = ， μi(i=1,2,…,n)的计算式如下： 　　 　　2.3 数值计算结果分析 　　 　　为探讨不同因素对锥形桩桩身屈曲稳定的影响规律，根据上述解答进行编程并进行对比分析。计算时，为方便起见，先将公式无量纲化，即令其中，l 为桩身无量纲长度， p l 为桩的无量纲稳定计算长度， 51 1 α = mb /(EI ) ；取b1=0.9(a1+1)；I1=πa12/64，准稳定计算长度20 1 / p cr l =π E I P ， 1 a 为桩端截面直径，E0 为桩身材料弹性模量，E=0.67E0。 　　 　　各参数取值：a1=1.0 米，E0=2.5×1010 帕，m=5×104 牛/米4。 　　 　　（1）取a0=0。此时锥形桩退化成等截面桩，桩身稳定分析结果如图2 所示。由图中可看出，稳定计算长度随着桩长的增大而增大。桩入土深度与桩长之比值越小，稳定计算长度增大的幅度也越大。此比值取最大值1 时，桩完全埋入土中，即低承台桩其稳定计算长度亦随桩长的增大而增大，但增幅极小。 　　 　　图2 所反映的这一p l - l 变化规律与文[8]、[9]关于等截面桩屈曲分析的结论完全一致。 　　 　　（2）取a0≠0。此时桩为严格意义上的锥形桩。分别h / l =0.5，0.7，1.0 得到的p l - l 关系，其曲线分别如图3、图4、图5 所示。从图3～图5 可以看出， 对于桩周土(m)、桩长（l）、桩端桩径(a1)及桩材(E)给定的锥形桩而言，其桩径变