差分法举例.ppt

* 差分法举例：桩基础的弹性理论法 POULOS的弹性理论法 － 单桩分析 1： 基本假定 （1）将土看作为均质的、各向同性的弹性半空间体，具有相同的弹性模量和泊松比，它们都不因有桩的存在而发生变化。 （2）将桩看作长度为L，直径为d，底端直径为 的一根圆桩，桩顶与地表面平齐，并作用有轴向外荷载 P；沿桩身圆周作用有均匀分布的剪应力；在桩端作用有均匀的竖向应力。 （3）桩身侧面假定是完全粗糙的。 （4）只考虑桩与其临近土之间的竖向位移协调，忽略它们之外的径向位移协调。 2：土的位移方程 在一般情况下，将桩划分为n个单元，例如取n =10，其精确度可以满足计算要求。 Fig. Poulos单桩分析 考虑单元i，单元j上的剪应力在i处产生的桩周土位移可表示为： 式中： ——单元j上的剪应力 =1时在i处产生的竖向位移系数。用MINDLIN公式计算。 全部n个单元上的剪应力和桩端上的竖向应力在i处产生的土位移 为： 式中： ——桩端竖向应力 =1 时在i处产生的竖向位移系数。用MINDLIN公式计算。 对于其他单元和桩端可以写出类似的表达式，于是桩所有单元的土位移可用矩阵形式表示为： 式中： ——土位移矢量， ——桩侧剪应力和桩端应力矢量 ——土位移系数的方阵，由下式给出： 3：桩的位移方程 假定桩材料的弹性模量和抗截面积均为常数。将面积定义为桩截面积同桩外周边包围的面积之比值。即： 对实心桩， =1。 对于桩单元，考虑圆柱桩单元的竖向平衡条件可得： 式中： ——桩的轴向应力 ——桩侧面的剪应力 分析桩各单元的位移时，忽略径向力的影响，只计轴向力的压缩作用，单元的轴向应变为： 式中： 为桩的位移。 合并上述两式，可得： （i =1,2………n） 将上述方程用差分方程表示，依次应用于计算点i =1,2…….n，可得桩位移方程为： 式中： ——n +1阶剪切矢量 ——n +1阶桩位移矢量 ——桩作用矩阵（n +1）方阵。如下式： 其中： 其中：P——作用于桩顶上的轴向荷载。 下面推倒矩阵 后两行系数的由来： 原理图如下： 基本原理就是将函数的一点在任意一点以泰勒级数展开如下所示： （1） （2） 对于桩侧，由于是受剪应力 作用则应对 进行差分，因 根据上式将Sn-2, Sn-1, Sb在Sn处展开得： （3） （4） （5） 将（3）*x +（4）*y +（5）*z 得到下式： 由于要得到 的值，则可以令 前的系数为1， 前的 系数为0，且忽略h的高次幂 ，因 时， 以h的四 次及更高次幂趋向于零。得方程组如下： 解方程组得：X=0.2 Y=2 Z=3.2 即有： 将上式中的 代入下式得： 而对于桩端，其受到的是竖向正应力 的作用，因 参看下面的原理图： 根据式(1)、(2)将 Sn-1, Sn在Sb处展开得： （6） （7） 将（6）*x+（7）*y 得到下式： 由于要得到 的值，则可以令 前的系数为4（为了组装成矩阵的需要，但也可以为1，只不过最后组装成矩阵时还是要乘以4的）， 前的系数为0，且忽略h的高次幂 ，因 时， 以h的三次及更高次幂趋向于零。得方程组如下： 解得： X=4/3 Y=-12 即有： 将上式中的 带入下式得： 其中： *