一种基于虚土桩模型的饱和层状土中桩基纵向振动研究方法技术

本发明专利技术提供一种基于虚土桩模型的饱和层状土中桩基纵向振动研究方法。本发明专利技术方法基于饱和虚土桩模型的浮承桩纵向振动动力阻抗算法系统，其采用的饱和虚土桩模型，能考虑桩周及桩底土体的饱和特性，且能合理考虑桩底土厚度等因素对桩基纵向振动特性的影响，适用于饱和土中浮承桩纵向振动特性研究，可为桩基动力检测提供理论指导和参考作用。

【技术实现步骤摘要】
一种基于虚土桩模型的饱和层状土中桩基纵向振动研究方法
本专利技术涉及土建
，具体而言，尤其涉及一种基于虚土桩模型的饱和层状土中桩基纵向振动研究方法。
技术介绍
目前，针对浮承桩情况的问题，已有研究大多将桩底土简化为Winkler模型，其弹簧和阻尼器系数通常按经验取值，无法合理考虑桩底土体波动效应的影响。为解决此类问题，一些学者将桩底土考虑为单相或饱和弹性半空间介质，计算得出桩底复阻抗函数表达式，并对浮承桩纵向振动特性进行了分析。然而，弹性半空间模型虽可考虑桩底土波动效应，但其只适用于基岩埋深较大情况，且无法考虑桩底土厚度及成层特性对桩基纵向振动特性的影响。基于此点考虑，提出了桩与桩底土完全耦合单相介质虚土桩模型。而单相虚土桩模型均假定桩底土体为单相介质，未考虑桩底土饱和两相介质性，这对于饱和土中浮承桩基纵向振动问题并不适用。
技术实现思路
根据上述提出的技术问题，而提供一种基于虚土桩模型的饱和层状土中桩基纵向振动研究方法。本专利技术采用的技术手段如下：一种基于虚土桩模型的饱和层状土中桩基纵向振动研究方法，包括以下步骤：S1、构建基于平面应变法及桩底饱和虚土桩模型的饱和土-桩-饱和虚土桩耦合体系纵向振动力学简化模型；S2、将桩-土耦合振动系统按土体沿纵向分成m层，其中桩底土分成n层，假定桩周和桩底每层土都为均质、各向同性的饱和线黏弹性介质，根据Biot动力波动理论，建立第j层饱和土体和第k层饱和土体的平面应变条件下饱和土层动力控制方程及实体桩与饱和虚土桩界面处的边界条件和桩-土耦合条件；S3、将第j层饱和土体和第k层饱和土体的谐和激振作用下饱和土层质点位移方程代入所述饱和土层动力学控制方程中，求解桩底第j层饱和土体土骨架径向位移和桩周第k层饱和土体土骨架径向位移，基于求解的径向位移求解出桩底土与饱和虚土桩界面处、桩周土与实体桩界面处剪应力，通过饱和虚土桩边界条件和各层桩界面处位移连续、力的平衡条件及阻抗函数传递公式，对实体桩桩顶动力阻抗函数进行求解；S4、由桩顶位移阻抗函数求得桩顶位移频率响应函数，求得桩顶速度频率响应函数、单位脉冲激励作用下桩顶速度时域响应，进而求得桩顶时域速度响应函数，基于求得的所述桩顶速度频率响应函数和桩顶时域速度响应函数对桩身振动特性及桩身完整性进行评判。进一步地，所述步骤S2假定条件还包括：桩底各层土体为渗透性较差的饱和粘土；桩周土和桩底土为一系列相互独立的薄层，不考虑土层间的相互作用；实体桩为均质等截面弹性体，桩身相邻层之间满足力平衡和位移连续条件；饱和虚土桩为等截面饱和两相介质，饱和虚土桩相邻层之间满足力平衡和位移连续条件，其与实体桩界面位移连续、应力平衡；桩-土耦合振动系统满足线弹性和小变形条件，桩土界面完全接触，不存在滑移和脱离。进一步地，所述步骤S2具体为：S21、根据Biot动力波动理论，建立平面应变条件下饱和土层动力学控制方程为：式中，各参数关系如下：其中，将桩-土耦合振动系统按土体沿纵向分成m层，其中桩底土分成n层，由基岩往上依次编号为1、…、j、…、n、n+1、…、k、…、m层，各层厚度分别为l1、…、lj…、ln、ln+1、…、lk、…、lm，各层顶部埋深分别为h1、…、hj…、hn、hn+1、…、hk、…、hm，j＝1、2、…、n，k＝n+1、n+2、…、m，式中，uj、wj为桩底第j层土骨架纵向位移和流体相对于土骨架纵向位移，Nj、ηj、bj、kj、Gj、ξj分别为桩底第j层饱和土体流体密度、土颗粒密度、孔隙率、流体粘滞系数、骨架与孔隙流体的粘性耦合系数、Biot定义的动力渗透系数、土体达西定律渗透系数、土体复剪切模量、土体剪切模量和阻尼比；uk、wk为桩周第k层土骨架纵向位移和流体相对于土骨架纵向位移，Nk、ηk、bk、kk、Gk、ξk分别为桩周第k层饱和土体流体密度、土颗粒密度、孔隙率、流体粘滞系数、骨架与孔隙流体的粘性耦合系数、Biot定义的动力渗透系数、土体达西定律渗透系数、土体复值剪切模量、土体剪切模量和阻尼比；g为重力加速度，S22、基于Biot动力波动理论得出渗透性较差时饱和土体一维纵向振动控制方程为：其中，为第j层饱和虚土桩纵向位移，λj为第j层桩底土拉梅常数，且λj＝2νjGj/(1-2νj)，νj为桩底第j层饱和土体泊松比。αj、Mj为Biot定义的表征桩底第j层饱和土体土颗粒及流体压缩性常数，和分别为桩底第j层饱和土体土颗粒、流体及土骨架的体积压缩模量；S23、将桩底第j层土在饱和虚土桩界面处的剪切应力代入式(4)中，可得桩底第j层饱和虚土桩纵向振动控制方程为：式中，为剪应力幅值，ω为激振圆频率，r0表示实体桩半径，Ap表示实体桩的横截面面积；取第k段实体桩桩身微元体作动力平衡分析可得第k段实体桩纵向振动控制方程为：式中，为第k段实体桩质点纵向振动位移，EP、ρP分别为实体桩弹性模量和密度，为第k层桩周饱和土体与实体桩界面处土体剪应力，为剪应力幅值。进一步地，所述步骤S2中，饱和土-桩-饱和虚土桩体系边界条件如下：饱和土体：径向无穷远处位移为零，即uj(∞,t)＝0(7a)uk(∞,t)＝0(7b)实体桩：桩顶平衡条件：各层界面位移连续条件：各层界面力平衡条件：饱和虚土桩：桩底位移：各层界面位移连续条件：各层界面力平衡条件：桩与饱和虚土桩界面处边界条件：位移连续条件：力平衡条件：桩-土耦合条件：其中，HP为实体桩桩长，表示桩顶作用谐和激振力，表示激振力幅值，H表示基岩上土层总厚度，Hsp表示桩底土层厚，即饱和虚土桩长度。进一步地，所述步骤S3中，谐和激振作用下饱和土层质点位移方程具体为：式中：为桩底第j层饱和土体土骨架径向位移和流体相对于土骨架径向位移响应幅值，为桩周第k层饱和土体土骨架径向位移和流体相对于土骨架径向位移响应幅值，i为虚数单位，所述步骤S3中，求解纵向位移解具体为：S31、将谐和激振作用下饱和土层质点位移方程(12)代入所述饱和土层动力学控制方程(1)、(2)中，可得：将式(13b)、(14b)分别代入式(13a)、(14a)可得：式中：方程(15)的通解为：式中，Aj、Bj、Ak、Bk为待定常数，I0(qjr)、I0(qkr)、K0(qjr)、K0(qkr)为零阶第一类、第二类虚宗量贝塞尔函数。由边界条件式(7)可知Bj＝Bk＝0，则进一步可得：所述步骤S3中，对实体桩桩顶动力阻抗函数进行求解具体为：S32、谐和激振作用下饱和虚土桩和实体桩质点纵向振动位移满足下式：式中，分别为第j层饱和虚土桩和第k层实体桩质点纵向振动位移响应幅值，则式(5)和式(6)可化简为：由边界条件式(11)可得：将式(20)代入式(17)，并由剪应力与位移的关系可得出桩底土与饱和虚土桩界面处、桩周土与实体桩界面处剪应力分别为：将式(21)代入式(19)，并进一步化简可得：式中，则方程(22)的通解为：式中：Cj、Dj、Ck、Dk为待定常数；S33、由饱和虚土桩边界条件式(9a)可得：根据位移阻抗函数的定义可得第1段饱和虚土桩与第2段饱和虚土桩界面处的位移阻抗函数为：综合第1段饱和虚土桩与第2段饱和虚土桩界面处位移连续、力平衡条件可得：由此可得第2段饱和虚土桩与第3段饱和虚土桩界面处的位移阻抗函数为：综合式(24)—(27)可得饱和虚本文档来自技高网...

【技术保护点】
1.一种基于虚土桩模型的饱和层状土中桩基纵向振动研究方法，其特征在于，包括以下步骤：S1、构建基于平面应变法及桩底饱和虚土桩模型的饱和土‑桩‑饱和虚土桩耦合体系纵向振动力学简化模型；S2、将桩‑土耦合振动系统按土体沿纵向分成m层，其中桩底土分成n层，假定桩周和桩底每层土都为均质、各向同性的饱和线黏弹性介质，根据Biot动力波动理论，建立第j层饱和土体和第k层饱和土体的平面应变条件下饱和土层动力控制方程及实体桩与饱和虚土桩界面处的边界条件和桩‑土耦合条件；S3、将第j层饱和土体和第k层饱和土体的谐和激振作用下饱和土层质点位移方程代入所述饱和土层动力学控制方程中，求解桩底第j层饱和土体土骨架径向位移和桩周第k层饱和土体土骨架径向位移，基于求解的径向位移求解出桩底土与饱和虚土桩界面处、桩周土与实体桩界面处剪应力，通过饱和虚土桩边界条件和各层桩界面处位移连续、力的平衡条件及阻抗函数传递公式，对实体桩桩顶动力阻抗函数进行求解；S4、由桩顶位移阻抗函数求得桩顶位移频率响应函数，求得桩顶速度频率响应函数、单位脉冲激励作用下桩顶速度时域响应，进而求得桩顶时域速度响应函数，基于求得的所述桩顶速度频率响应函数和桩顶时域速度响应函数对桩身振动特性及桩身完整性进行评判。...

【技术特征摘要】
1.一种基于虚土桩模型的饱和层状土中桩基纵向振动研究方法，其特征在于，包括以下步骤：S1、构建基于平面应变法及桩底饱和虚土桩模型的饱和土-桩-饱和虚土桩耦合体系纵向振动力学简化模型；S2、将桩-土耦合振动系统按土体沿纵向分成m层，其中桩底土分成n层，假定桩周和桩底每层土都为均质、各向同性的饱和线黏弹性介质，根据Biot动力波动理论，建立第j层饱和土体和第k层饱和土体的平面应变条件下饱和土层动力控制方程及实体桩与饱和虚土桩界面处的边界条件和桩-土耦合条件；S3、将第j层饱和土体和第k层饱和土体的谐和激振作用下饱和土层质点位移方程代入所述饱和土层动力学控制方程中，求解桩底第j层饱和土体土骨架径向位移和桩周第k层饱和土体土骨架径向位移，基于求解的径向位移求解出桩底土与饱和虚土桩界面处、桩周土与实体桩界面处剪应力，通过饱和虚土桩边界条件和各层桩界面处位移连续、力的平衡条件及阻抗函数传递公式，对实体桩桩顶动力阻抗函数进行求解；S4、由桩顶位移阻抗函数求得桩顶位移频率响应函数，求得桩顶速度频率响应函数、单位脉冲激励作用下桩顶速度时域响应，进而求得桩顶时域速度响应函数，基于求得的所述桩顶速度频率响应函数和桩顶时域速度响应函数对桩身振动特性及桩身完整性进行评判。2.根据权利要求1所述的基于虚土桩模型的饱和层状土中桩基纵向振动研究方法，其特征在于，所述步骤S2假定条件还包括：桩底各层土体为渗透性较差的饱和粘土；桩周土和桩底土为一系列相互独立的薄层，不考虑土层间的相互作用；实体桩为均质等截面弹性体，桩身相邻层之间满足力平衡和位移连续条件；饱和虚土桩为等截面饱和两相介质，饱和虚土桩相邻层之间满足力平衡和位移连续条件，其与实体桩界面位移连续、应力平衡；桩-土耦合振动系统满足线弹性和小变形条件，桩土界面完全接触，不存在滑移和脱离。3.根据权利要求1所述的基于虚土桩模型的饱和层状土中桩基纵向振动研究方法，其特征在于，所述步骤S2具体为：S21、根据Biot动力波动理论，建立平面应变条件下饱和土层动力学控制方程为：式中，各参数关系如下：其中，将桩-土耦合振动系统按土体沿纵向分成m层，其中桩底土分成n层，由基岩往上依次编号为1、…、j、…、n、n+1、…、k、…、m层，各层厚度分别为l1、…、lj…、ln、ln+1、…、lk、…、lm，各层顶部埋深分别为h1、…、hj…、hn、hn+1、…、hk、…、hm，j＝1、2、…、n，k＝n+1、n+2、…、m，式中，uj、wj为桩底第j层土骨架纵向位移和流体相对于土骨架纵向位移，Nj、ηj、bj、kj、Gj、ξj分别为桩底第j层饱和土体流体密度、土颗粒密度、孔隙率、流体粘滞系数、骨架与孔隙流体的粘性耦合系数、Biot定义的动力渗透系数、土体达西定律渗透系数、土体复剪切模量、土体剪切模量和阻尼比；uk、wk为桩周第k层土骨架纵向位移和流体相对于土骨架纵向位移，Nk、ηk、bk、kk、Gk、ξk分别为桩周第k层饱和土体流体密度、土颗粒密度、孔隙率、流体粘滞系数、骨架与孔隙流体的粘性耦合系数、Biot定义的动力渗透系数、土体达西定律渗透系数、土体复值剪切模量、土体剪切模量和阻尼比；g为重力加速度，S22、基于Biot动力波动理论得出渗透性较差时饱和土体一维纵向振动控制方程为：其中，为第j层饱和虚土桩纵向位移，λj为第j层桩底土拉梅常数，且λj＝2νjGj/(1-2νj)，νj为桩底第j层饱和土体泊松比。αj、Mj为Biot定义的表征桩底第j层饱和土体土颗粒及流体压缩性常数，和分别为桩底第j层饱和土体土颗粒、流体及土骨架的体积压缩模量；S23、将桩底第j层土在饱和虚土桩界面处的剪切应力代入式(4)中，可得桩底第j层饱和虚土桩纵向振动控制方程为：式中，为剪应力幅值，ω为激振圆频率，r0表示实体桩半径，Ap表示实体桩的横截面面积；取第k段实体桩桩身微元体作动力平衡分析可得第k段实体桩纵向振动控制方程为：式中，为第k段实体桩质...

【专利技术属性】
技术研发人员：崔春义，孟坤，许成顺，梁志孟，刘海龙，
申请(专利权)人：大连海事大学，
类型：发明
国别省市：辽宁,21