基于黏性阻尼模型复杂非均质土中管桩纵向振动分析方法与流程

本发明方法涉及土建理论分析领域，更具体地，涉及一种基于黏性阻尼模型的径向非均质、纵向成层土体中管桩纵向振动分析方法。

背景技术：

纵向成层土中的桩基相互作用体系耦合振动特性研究是桩基抗震、防震设计及桩基动力检测等工程技术领域的理论基础，一直以来亦是土动力学、岩土工程和结构-地基相互作用领域的交叉热点问题。

桩-土耦合振动特性研究是桩基抗震、防震设计及桩基动力检测等工程技术领域的理论基础，一直以来亦是岩土工程和固体力学的热点问题。

众所周知，在桩基施工过程中，由于挤土、松弛以及其它扰动因素的影响，使得桩周土体沿桩基径向存在一定不均匀性，即径向非均质效应。为考虑此种径向非均质效应，国内外诸多学者取得了大量成果。这些成果可从不同角度加以分类，从作用的外荷载来看，可分为谐和荷载作用下的频域响应研究和任意荷载下时域、频域响应研究；从土体的材料阻尼来看，可分为滞回材料阻尼和粘性材料阻尼；从求解方法来看，可分为解析法、半解析法及数值方法。

土体的材料阻尼是由土体内部颗粒摩擦所引起的能量耗散，这种内摩擦是由介质颗粒结晶结构的缺损、介质颗粒之间的非弹性连接及其他热弹性过程引起的，是不可避免的，为了考虑这一内摩擦效应，采用考虑阻尼效应的土体线性本构方程，来研究材料阻尼对桩动力响应的影响是非常必要的。

在观测和实验基础上建立的常用线性阻尼本构方程可分为两类：时域本构方程和频域本构方程，前者从宏观物理模型线性粘弹性体出发直接在时域建立；后者则通过与经典的频域分析方法相匹配在频域内建立。

线性粘弹性体的时域本构模型，可以由线性弹簧和线性阻尼元件构成，线性阻尼元件的粘性应力与应变率成正比，由这两种线性单元可以构成各种线性粘弹性本构模型，可以反映真实固体的应力-应变性质。

线性滞回阻尼主要体现在频域本构中的滞回阻尼比，频域本构可以理解为时域本构的逆傅里叶变换，滞回阻尼比通常假设为常数，即假设材料处于弹性工作区域内，滞回阻尼比的变化不大，或无明显趋向性变化。另外，对谐和荷载下的稳态振动问题的频域分析，能够近似地反映土体的材料阻尼特性。然而，对非谐和振动(瞬态振动或随机振动)问题，滞回阻尼模型是不适合的，特别是在研究瞬态激振条件下桩的时域响应时，土阻尼力与振幅有关也与应变速率有关，采用滞回阻尼模型在概念上会引起矛盾，从而产生所谓“动响应的非因果性”，而此时粘性阻尼模型则比较适合，在物理上也更合理。

另外，目前大部分研究均是针对实心桩展开，而对于大直径管桩，由于桩芯土的存在，必然使得其与实心桩的振动特性存在差异。丁选明等和郑长杰等同时考虑桩周土和桩芯土，对径向均质土中管桩振动特性进行求解，并与实心桩结果进行对比，说明了在竖向荷载作用下管桩表现出与实心桩动力特性的不同。

技术实现要素：

本发明的目的在于克服现有技术存在的技术缺陷，考虑桩周土体施工扰动及土体纵向成层特性，土体采用黏性阻尼模型，基于复刚度传递多圈层平面应变模型，对任意激振力作用下径向非均质、纵向成层黏性阻尼土中管桩纵向振动特性进行解析分析。

为实现上述目的，本发明的技术方案如下：

一种基于黏性阻尼模型的径向非均质、纵向成层土体中管桩纵向振动分析方法，包括以下步骤：

s1：管桩的桩周土系统分为桩周土体和桩芯土体，假定桩周土体和桩芯土体为一系列相互独立的薄层，忽略各个薄层间的相互作用；

s2：桩周土体沿纵向分为任意个层段，每个层段的桩周土体沿径向划分为内部扰动区域和外部区域，并将内部扰动区域沿径向划分任意个圈层，每一圈层内的土体各自为均质、各向同性线性粘弹性体，外部区域的土体径向无限延伸，土体材料阻尼采用黏性阻尼，忽略土体的径向位移；

s3：桩周土界面及各圈层土界面两侧位移连续、应力平衡，且桩周土系统振动为小变形；

s4：管桩桩身的混凝土为线弹性，应力波在管桩桩身中的传播满足平截面假定；

s5：根据弹性动力学基本理论，建立平面应变条件下的桩周土体、桩芯土体和管桩桩身纵向振动方程及边界条件；

s6：使用laplace(拉普拉斯)变换，求解s5中所述的三个振动方程，并利用纵向层段间桩身阻抗函数的传递性得到任意激振力作用在桩顶的时域速度响应函数，以对管桩的纵向振动进行分析。

桩周土体振动方程：

桩芯土体振动方程：

符合平截面假定的管桩桩身纵向振动方程为：

其中，将桩周土体沿纵向划分为m个层段，将桩长为h管桩自桩身底部由下往上依次编号为1、2、…、i、…、m层段，各层段厚度分别为l1、l2、…、li、…、lm，各层段顶部埋深分别为h1、h2、…、hi…、hm。管桩内径、管桩外径、管桩桩段截面积、管桩密度和管桩弹性模量分别为ri0、ri1、和桩底黏弹性支承刚度系数为δp、kp。在纵向第i层段中，桩芯土体剪切模量、黏性阻尼系数和密度分别为gi0、ηi0、ρi0。同时，将纵向第i层段桩周土体沿径向划分为内部扰动区域和外部区域，桩周土体内部扰动区域径向厚度为bi，并将内部扰动区域沿径向划分n个圈层，第j圈层土体剪切模量、黏性阻尼系数和密度分别为gij、ηij、ρij，第j-1个圈层与第j圈层的界面处半径为rij。内部区域和外部区域界面处的半径为ri(n+1)，外部区域则为径向半无限均匀黏弹性介质。管桩桩顶作用任意激振力p(t)，第i层段桩芯土和桩周土对桩身产生的切应力分别为和在第i层段中，设桩芯土体位移为桩周土体第j圈层中某点位移为第i层段桩身位移为r为径向位移，t为时间，z为纵向位移，epi为第i层段桩身弹性模量，api为第i层段桩身截面积；

s5中的边界条件包括：

桩芯土体边界条件：

当r＝0时，位移为有限值：

桩芯土体与桩位移及力连续条件：

ri0为桩内半径，为桩芯土体对桩身产生的切应力，为桩芯土体在管桩内壁的竖向剪应力，顺时针为正；

桩周土体边界条件：

当r＝∞时，位移为零：

式中，代表第i层段土体外部区域位移。

桩周土体与桩位移及力连续条件：

其中，ri1为桩外半径，为第1圈层土体位移，为桩周土体对桩身产生的切应力，为桩周土体在管桩外壁的竖向剪应力，顺时针为正；

桩身边界条件：

桩顶作用力为p(t)：

桩端处边界条件：

为桩身弹性模量，为桩身截面积，kp，δp为桩底黏弹性支承常数。

s6包括以下具体步骤：

步骤1：对方程、进行laplace(拉普拉斯)变换，得到纵向第i层段基于黏性阻尼的多圈层平面应变模型的土层剪切刚度递推公式为：

其中

其中，rij为第i层段第j圈层土的内边界，ri(j+1)为第i层段第j圈层土的外边界，为第i层段第j圈层土固有参数，s为复变量，为第i层段第j圈层土的内边界的剪切刚度，为第i层段第j圈层土的外边界的剪切刚度，i0、i1为零阶和一阶第一类修正bessel(贝塞尔)函数，k0、k1零阶和一阶第二类修正bessel(贝塞尔)函数；

步骤2：对方程和进行laplace(拉普拉斯)变换，得到第i层段管桩内壁受到桩芯土体的剪切刚度公式：

其中，为桩芯土固有参数；

步骤3：对方程、和进行laplace(拉普拉斯)变换，根据桩身相邻层段界面处力平衡和位移连续性条件，得到桩身相邻纵向层段阻抗函数递推公式：

式中，为第i层段、第i-1层段桩身阻抗函数，αi、βi为求解化简参数，li为第i层段桩身长度，hi、hi-1分别为桩顶到第i层段、第i-1层段桩身顶部的长度。

步骤4：利用桩身阻抗函数的传递性得到桩顶复动刚度公式：

其中k′d为桩顶复刚度kd的无量纲参数，令k′d＝kr+iki，kr代表桩顶动刚度，ki代表桩顶动阻尼，αm、βm为求解化简参数，lm为第m段桩长；

步骤:5：根据式(17)得到桩顶速度导纳函数hv：

其中，为第m段桩身密度，为第m段桩身弹性波速，hv′为桩顶速度导纳函数hv的无量纲化；ω为纵向振动圆频率。

步骤6：根据(18)得到单位脉冲激励的时域响应为：

式中t′＝t/tc为无量纲时间，θ为无量纲频率；ift为快速傅里叶逆变换符号；

步骤7：根据卷积定理得到任意激振力p(t)作用在桩顶的时域速度响应函数

g(t)＝p(t)*h(t)＝ift[p(iω)·h(iω)](20)

其中，h(t)为单位脉冲激励作用下时域速度响应，h(iω)为桩顶速度频率响应函数。

步骤8：激振力p(t)为半正弦脉冲激励t∈(0,t)时，t为脉冲宽度，桩顶时域速度响应的半解析解答为：

其中，qmax为半正弦脉冲振幅，vv′为时域响应无量纲速度。

从上述技术方案可以看出，本发明通过采用径向非均质、纵向成层黏性阻尼土体模型对大直径管桩的纵向振动进行分析，黏性阻尼土体模型的阻尼力与应变速率相关，能适用于非谐和激振问题，特别是瞬态激振条件下时的桩体时域振动响应问题，同时，径向非均质性能考虑桩周土体施工扰动效应，更接近现实模型，另外，考虑了桩芯土对管桩振动特性的影响，使计算精度更高，可为桩基动力检测提供理论指导和参考作用。

附图说明

图1是本发明的基于黏性阻尼模型的径向非均质、纵向成层土体中管桩纵向振动分析方法的流程图。

图2是本发明的桩土系统纵向耦合振动力学简化模型的示意图。

具体实施方式

下面结合附图，对本发明的具体实施方式作进一步的详细说明。

请参阅图1，图1是本发明的基于黏性阻尼模型的径向非均质、纵向成层土体中管桩纵向振动分析方法的流程图。如图1所示，一种基于黏性阻尼模型的径向非均质、纵向成层土体中管桩纵向振动分析方法，包括以下步骤：

s1：假定桩周土和桩芯土为一系列相互独立的薄层，忽略土层间相互作用；

s2：桩周土体沿纵向分为任意个层段，每个层段桩周土体沿径向划分为内部扰动区域和外部区域，并将内部扰动区域沿径向划分任意个圈层，每一圈层土体各自为均质、各向同性线性粘弹性体，外部区域土体径向无限延伸，土体材料阻尼采用黏性阻尼，忽略土体径向位移；

s3：桩土界面及各圈层土界面两侧位移连续、应力平衡，且桩土系统振动为小变形；

s4：桩身混凝土为线弹性，应力波在桩身中的传播满足平截面假定；

本发明基于平面应变模型，对纵向任意层段、径向任意圈层土中的黏弹性支承管桩的纵向振动特性进行研究，力学简化模型如图2所示。将桩周土体沿纵向划分为m个层段，将桩长为h管桩自桩身底部由下往上依次编号为1、2、…、i、…、m层段，各层段厚度分别为l1、l2、…、li、…、lm，各层段顶部埋深分别为h1、h2、…、hi…、hm。管桩内径、外径、桩段截面积、密度和弹性模量分别为ri0、ri1、ai^p、ρi^p和ei^p，桩底黏弹性支承刚度系数为δp、kp。在纵向第i层段中，桩芯土体剪切模量、黏性阻尼系数和密度分别为gi0、ηi0、ρi0。同时，将纵向第i层段桩周土体沿径向划分为内部扰动区域和外部区域，桩周土体内部扰动区域径向厚度为bi，并将内部扰动区域沿径向划分n个圈层，第j圈层土体剪切模量、黏性阻尼系数和密度分别为gij、ηij、ρij，第j-1个圈层与第j圈层的界面处半径为rij。特别地，内部区域和外部区域界面处的半径为ri(n+1)，外部区域则为径向半无限均匀黏弹性介质。管桩桩顶作用任意激振力p(t)，第i层段桩芯土和桩周土对桩身产生的切应力分别为和

s5：根据弹性动力学基本理论，建立平面应变条件下的桩周土体、桩芯土体和桩身纵向振动方程及边界条件；

s6：使用laplace变换，求解s5中所述的三个振动方程，并利用纵向层段间桩身阻抗函数的传递性得到任意激振力作用在桩顶的时域速度响应函数，以对管桩的纵向振动进行分析。

具体地，包括以下具体步骤：

步骤1：将桩周土体沿纵向划分为m个层段，将桩长为h管桩自桩身底部由下往上依次编号为1、2、…、i、…、m层段，各层段厚度分别为l1、l2、…、li、…、lm，各层段顶部埋深分别为h1、h2、…、hi…、hm。管桩内径、外径、桩段截面积、密度和弹性模量分别为ri0、ri1、和桩底黏弹性支承刚度系数为δp、kp。在纵向第i层段中，桩芯土体剪切模量、黏性阻尼系数和密度分别为gi0、ηi0、ρi0。同时，将纵向第i层段桩周土体沿径向划分为内部扰动区域和外部区域，桩周土体内部扰动区域径向厚度为bi，并将内部扰动区域沿径向划分n个圈层，第j圈层土体剪切模量、黏性阻尼系数和密度分别为gij、ηij、ρij，第j-1个圈层与第j圈层的界面处半径为rij。特别地，内部区域和外部区域界面处的半径为ri(n+1)，外部区域则为径向半无限均匀黏弹性介质。管桩桩顶作用任意激振力p(t)，第i层段桩芯土和桩周土对桩身产生的切应力分别为和在第i层段中，设桩芯土体位移为桩周土体第j圈层中某点位移为第i层段桩身位移为r为径向位移，t为时间，z为纵向位移，epi为第i层段桩身弹性模量，api为第i层段桩身截面积，根据弹性动力学基本理论，建立平面应变条件下的桩周、桩芯土体和桩身纵向振动方程和边界条件分别如下：

桩周土体振动方程：

桩芯土体振动方程：

符合平截面假定的桩身纵向振动方程为：

桩芯土边界条件：

当r＝0时，位移为有限值：

桩芯土与桩位移及力连续条件：

ri0为桩内半径，为桩芯土对桩身产生的切应力，为桩芯土在管桩内壁的竖向剪应力，顺时针为正；

桩周土边界条件：

当r＝∞时，位移为零：

式中，代表第i层段土体外部区域位移。

桩周土与桩位移及力连续条件：

其中，ri1为桩外半径，为第1圈层土体位移，为桩周土对桩身产生的切应力，为桩周土在管桩外壁的竖向剪应力，顺时针为正；

桩身边界条件：

桩顶作用力为p(t)：

桩端处边界条件：

其中，为桩身弹性模量，为桩身截面积，kp，δp为桩底黏弹性支承常数。

进一步的，所述s6包括以下具体步骤：

步骤1：对方程、进行laplace变换，得到纵向第i层段基于黏性阻尼的多圈层平面应变模型的土层剪切刚度递推公式为：

其中

其中，rij为第i层段第j圈层土的内边界，ri(j+1)为第i层段第j圈层土的外边界，为第i层段第j圈层土固有参数，s为复变量，为第i层段第j圈层土的内边界的剪切刚度，为第i层段第j圈层土的外边界的剪切刚度，i0、i1为零阶和一阶第一类修正bessel函数，k0、k1零阶和一阶第二类修正bessel函数；

步骤2：对方程和进行laplace变换，得到第i层段管桩内壁受到桩芯土体的剪切刚度公式：

其中，为桩芯土固有参数；

步骤3：对方程、和进行laplace变换，根据桩身相邻层段界面处力平衡和位移连续性条件，得到桩身相邻纵向层段阻抗函数递推公式：

式中，为第i层段、第i-1层段桩身阻抗函数，αi、βi为求解化简参数，li为第i层段桩身长度，hi、hi-1分别为桩顶到第i层段、第i-1层段桩身顶部的长度。

步骤4：利用桩身阻抗函数的传递性得到桩顶复动刚度公式：

其中k′d为桩顶复刚度kd的无量纲参数，令k′d＝kr+iki，kr代表桩顶动刚度，ki代表桩顶动阻尼，αm、βm为求解化简参数，lm为第m段桩长；

步骤5：根据(17)式得到桩顶速度导纳函数：

其中，为第m段桩身密度，为第m段桩身弹性波速，hv′为桩顶速度导纳函数hv的无量纲化；

步骤6：根据(18)得到单位脉冲激励的时域响应为：

式中t′＝t/tc为无量纲时间，θ为无量纲频率；ift为快速傅里叶逆变换符号；

步骤7：根据卷积定理得到任意激振力p(t)作用在桩顶的时域速度响应函数

g(t)＝p(t)*h(t)＝ift[p(iω)·h(iω)](20)

其中，h(t)为单位脉冲激励作用下时域速度响应，h(iω)为桩顶速度频率响应函数。

步骤8：激振力p(t)为半正弦脉冲激励t∈(0,t)时，t为脉冲宽度时，桩顶时域速度响应的半解析解答为：

其中，qmax为半正弦脉冲振幅，vv′为时域响应无量纲速度。

进一步的，基于桩顶速度导纳函数和桩顶速度时域响应函数，可以对桩身振动特性及桩身完整性进行评价。

综上所述，本发明的基于黏性阻尼模型的径向非均质、纵向成层土体中管桩纵向振动分析方法，其采用的阻尼模型为桩土耦合振动体系提供的阻尼力与应变速率相关，能适用于非谐和激振问题，特别是瞬态激振条件下时，桩体时域振动响应问题，而径向非均质性能考虑桩周土体施工扰动效应，纵向成层可以考虑土体天然沉积引起的成层特性，可为桩基动力检测提供理论指导和参考作用。

以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，根据本发明的技术方案及其发明构思加以等同替换或改变，都应涵盖在本发明的保护范围之内。