饱和土中端承桩纵向振动引起土层动力反应阻抗研究.doc

饱和土中端承桩纵向振动引起土层动力反应阻抗研究 摘要：基于Biot动力固结方程，忽略桩土系统竖向振动时土体的径向位移，直接对方程进行求解，得到端承桩在简谐荷载作用下引起饱和土层动力反应的阻抗函数表达式，并讨论了桩的长径比、桩土模量等对土层动力频响的影响规律，并得到若干重要结论。 关键词： 桩土动力相互作用；饱和土；纵向振动；复阻抗 引言 在桩土动力相互作用理论分析中，基于Novak提出的平面应变假定及其演化理论形式的模型[1][2][3]，以其简便实用的特性得到了广泛的应用。平面应变模型将桩周土离散为相互独立的一系列水平薄层，忽略各土层之间的联系，忽略桩土系统竖向振动时土的径向位移，Mylonalis系统的分析了忽略径向运动的含义[4]，指出这以近似的误差对于桩来说较小，可以满足工程要求。但这些研究成果均为在单相土介质当中取得的，和实际有一定出入，尚需在多相土介质中进行进一步研究。 Biot首先建立饱和土中波的传播理论，并提出理论土体三维固结问题的理论及相应的精确表达式[5][6]。由于该方程考虑土骨架与孔隙液体的耦合，虽然更符合实际，但也增加了求解的难度。目前，对于流固耦合的解耦，通常有两种方法，一种是通过积分变换将动力方程转化为常微分方程，再进行求解，另一种是引入势函数，对方程进行解耦，然后求解。以上两种方法均要进行复杂的数学求解。 本文从Biot饱和土动力固结理论出发，忽略土体径向运动，在不引入势函数和不进行积分变化情况下，用解析方法直接对方程进行求解，对端承桩在饱和土层中做纵向振动时引起的土层动力反应特性方面进行了一些研究工作。 基本假定 本文研究简谐荷载竖向振动下，端承桩与饱和土耦合振动问题。计算简图如图1所示。为简化计算，引入下列假定：（1）、桩周土体性质符合Biot模型要求；（2）、土层表面为自由边界，桩底为刚性支承；（3）、桩是垂直设置，具有圆形截面，桩身材料为理想线弹性材料，桩长为H；（4）、竖向振动时，桩周土的径向位移可以忽略；（5）、桩土系统振动为小变形，桩与周围土体完全连续接触；（6）、桩顶外荷载为简谐荷载，其幅值为，系统的振动达到稳态振动。 图1 桩土系统动力学模型简图 定解问题 桩土之间相互作用主要是由桩与土界面的连续条件来体现，可以通过以下4个方面来解决：1土体的振动分析2、桩身的振动分析3、土与桩的接触连续性4、桩土系统的初始条件及边界条件。 2.1土体动力平衡方程 土体的振动符合Biot振动方程。在轴对称条件下，依据Biot波动原理，忽略流体对土骨架的惯性项，在竖向荷载作用下饱和土的控制方程可以表示为： （1） （2） 流体的运动方程为： （3） （4） 式中：，为土骨架径向和轴向位移；，为流体相对土骨架的径向和轴向位移；，，分别为两相介质、土骨架和流体的密度，且有，为土体的孔隙率；，为Lame常数；pf为超静孔隙水压力；kd’为土的动力渗透系数，，g为重力加速度；e为土骨架的体积应变，；为Laplace算子，。 假设土颗粒和流体都不可压缩，土体的渗流连续方程为： （5） 2.2桩动力平衡方程 取桩身纵向振动位移幅值为u(z),取桩身微元体做动力平衡分析，可得桩作纵向振动的方程为： （6） 式中：E p为桩身材料杨氏模量，A为桩截面面积，m为单位桩长的质量，r0为桩半径， f（z）为桩周土对桩身单位面积的侧壁切应力。 2.3桩土连续条件及定解条件 桩土接触连续性表现为接触面的位移协调及力协调。故此，有及。 土体与桩初始条件及外部边界如下： （1）土层表面为自由表面，即； （2）土、桩底部支承为刚性基岩，即； （3）桩顶为幅值为P0的简谐荷载激振，即。 方程求解 对式（1）~(2)做的运算，可有： （7） 式中。 对式（3）~(4)做的运算，可有： (8) 联立（7）、（8）式，可得： （9） 由于是稳态简谐振动，令（以下讨论均省略项），代入（9）式中有： （10） 两边同除以W（z）R（r）有 （11） 上式可以分离成以下两个微分方程 （12） （13） 式中常数h，q