波动方程在动力打桩中的工程应用.pptxVIP

田敏 21160933059 建筑与土木工程 波动方程在动力打桩中的工程应用 前言 1 动力打桩一维波动方程的改进 2 土体模型的选择 3 动力打桩一维波动方程的求解 4 工程实例分析 5 波动方程 的工程应用 结语 6 前言 1 问题由来 对桩的可打入性进行准确的分析和预测是打桩施工的关键环节，同时也关系到整个工程的进度与成败。近年来，在土木工程、海洋工程和港口工程中都出现了大直径、超长桩，由于这种桩的重量大对桩的可打入性有显著影响，考虑桩身自重和桩周土阻力，对打桩分析的一维应力波动方程进行改进、推导。 波动方程 的工程应用 波动方程的改进 2 问题分析 本研究在原有的打桩分析一维应力波动方程中， 考虑桩身自重和桩周土阻力的影响，对波动方程进行修正和推导，并基于桩侧土体和桩端土体在竖向荷载作用下变形和破坏机理的不同，采用不同土体模型分别计算桩侧、桩端土体的静阻力，应用有限差分法求解、编制计算程序，应用于分析实际工程的打桩施工。 波动方程 的改进 分析 第一季度 单击编辑副标题 第一季度 单击编辑副标题 D.V. Isaacs(1931 年)将桩周土体阻力R′ 引入古典的一维波动方程得到 式中 x 为桩截面的位置标； u 为x 处桩截面的质点位移；t 为时间； R′ 为桩身土阻力，C 为弹性应力波波速，C = E / ρ ，E, ρ 分别为弹性模量和密度。 公式中存在两方面的不足： （1）仅仅在古典的一维波动方程中引入桩周土体阻力R，并没有真正反映出打桩过程中的波动响应。对方程进行量纲分析可知，桩周土体阻力R′ 的量纲为1/m，而实际力的量纲应为N，因此R′的物理意义不明确； （2）对于大直径（超大直径）、超长桩，其桩身重量很大，对桩的可打入性有显著影响。因此，有必要在一维动力打桩波动方程中考虑桩身自重和桩周土阻力。 波动方程的改进 2 两点问题 第一季度 单击编辑副标题 第一季度 单击编辑副标题 桩单元受力分析 波动方程的改进 2 波动方程的改进 2 桩单元受力分析 假设桩在贯入过程中，受到冲击荷载的撞击后任一截面在变形仍保持平面，则微分单元i 处的应变ε 为： （式中 u 为截面的位移（m），是时间t 及截面位置x的函数。） 桩微分单元 i 在截面ab 处的受力为： 在截面 cd 处的受力为 波动方程的改进 2 桩单元受力分析 微分单元 i 上所受合力为： （式中 R 为微分单元所受的土阻力（N）；G 为微分单元 i 的重量（N）。） 合力产生的加速度a 可按下式计算： 整理得到： 故一维动力打桩的波动方程为: 波动方程的改进 2 桩单元受力分析 右边第2 项为桩的自重对应力波传播的影响；第3 项为桩周土体阻力对应力波传播的影响。显然，该方程与式原波动方程明显不同，且有： 由此可见， R′ 为反映桩周土体阻力的一个参数，其量纲为1/m，而不是桩周土阻力。 土体模型的选择 3 土体模型的选择 3 动力打桩过程中，桩侧土体单元i 在时刻t 所发挥的静阻力和动阻力分别由非线性弹簧和缓冲壶组成的模型来模拟 桩侧土体模型 波动方程的求解 4 波动方程的求解 4 桩的计算模型 采用差分法求解波动方程。将整个打桩系统（桩锤、垫层、桩帽和桩）在空间上离散成若干个由刚 性质量块和无质量的弹簧（模拟弹性）组成的单元，如下图所示： 波动方程的求解 4 单元的长度为Δl ，在计算中将一次锤击历时划分为多个时间间隔Δt ，选取相当短的时间间隔Δt ，使得弹性应力波在Δt 内来不及由一个单元传播到下一个单元，这样所取的Δt 单元的运动可近似作为等速运动。若以 tcr 为临界时间间隔则： 则任一桩单元 i 在时刻t 时的平衡方程式为： 采用向后差分可得： 波动方程的求解 4 可得单元的位移： 因此， t 时刻单元i 的变形为： 在 t 时刻单元i 的受力为： 在 t 时刻单元i 的受力为： V (i,t)可用于计算下一个Δt的位移： 波动方程的求解 4 以锤心撞击垫层的接触瞬间为计算的初始时刻t = 0，由于在初始时刻之前，整个桩、土系统处于静止状态，故在t =0 时桩单元的弹簧力、土阻力、位移、速度及加速度均为零。以锤心的锤击初速度为已知边界条件，进行第一个Δt内应力波在整个系统内传播的计算分析。锤心在时间段Δt内产生的位移即为锤心弹簧的变形量，从而可算出作用在下一个单元上的力。这个力使锤心速度减小，使锤心下面的单元产生加速度并获得新速度。 同理，在每一个Δt时段内逐个单元进行迭代计算，直到 满足下面两个条件： ① 桩单元的位移量不再增加； ② 各单元的速度均已为零或负值。通过上述的运算就可以得到桩在一次锤击中的性状。 采用上述方法编制计算软件用于实际工程的打桩分析。 工程实例分析 5 工程实例分析 5 各个土层的物理力学指标 某海洋石油平