桩基纵向振动特性研究作为岩土工程和固体力学的热点问题^[1-3]，对桩基抗震、防震设计及动力测桩技术发展具有重要的理论意义和实际应用价值.桩基施工会引起桩周土体产生挤土、松弛现象，即径向非均质效应^[4-5]，基于此，诸多国内外学者做了较多相关研究.Novak等^[6-7]考虑桩周土体单层内部扰动，基于滞回阻尼平面应变模型推导了存在单层软化区域土体的桩基纵向、扭转动力阻抗解析解答，并就桩周土体非均质性对桩-土耦合振动特性的影响规律进行了初步分析.Nogami等^[8-9]基于平面应变模型考虑远场区域土体作用，近场区域土体则采用非线性Winkler模型，并将两者作用相互串联考虑桩周土体径向非均质性，对径向非均质桩周土中桩基纵向振动特性进行分析.进一步地，EI Naggar^[10]将桩周围土体沿径向分为内部扰动区和外部未扰动区，内部区域划分为多个同心圆，首先计算各圈层的复刚度，并将其串联以考虑桩周土的径向非均质特性，分析了径向多圈层土中桩基纵向振动特性.王奎华和杨冬英等^[11-12]研究发现EI Naggar所提出的模型在理论上不严谨，计算结果与实测结果存在较大误差.在此基础上，Wang和Yang等^[13-14]基于径向复刚度传递，土体材料阻尼采用滞回阻尼模型，提出了多圈层平面应变模型，建立了能更为合理地考虑施工扰动引起的桩周土径向非均质效应的桩-土耦合纵向振动简化力学模型，并对桩顶动力阻抗影响因素进行参数化分析.

以上研究均对象均为实体桩，而针对大直径管桩的相关研究较为少见.刘汉龙等^[15]基于径向不变假定，不考虑桩周土体作用探讨了管桩桩顶动力响应及平截面假定的影响深度问题.在此基础上，丁选明等^[16]以Winkler模型模拟桩周土和桩芯土作用，祛除“径向不变假定”，建立了任意段变阻抗薄壁管桩动力响应波动方程，分析了变截面桩和变模量的动力响应特性.进一步地，丁选明等^[17]基于平面应变模型, 土体阻尼采用黏性阻尼，推导出了均质土中管桩动力响应的频域解析表达式.郑长杰等^[18]则基于三维轴对称模型，考虑土体材料的滞回阻尼和桩周土-桩-桩芯土耦合振动，对均质黏弹性地基中现浇大直径管桩纵向振动频域特性进行理论研究.Zheng等^[19]基于三维轴对称模型同时考虑土体纵向位移和径向位移，土体材料阻尼选用滞回阻尼模型，求解了均质地基中管桩纵向振动频域解析解.吴文兵等^[20]基于平面应变模型，土体材料阻尼采用滞回阻尼，桩底支承采用虚土桩模型，以附加质量考虑土塞，研究了均质地基中土塞效应对管桩纵向振动特性的影响.刘林超等^[21]将桩周土、桩芯土和管桩视为一个整体，不考虑径向位移和环向位移，基于三维轴对称模型，利用桩-土接触面的连续条件得到均质土中管桩纵向振动解析解.管桩在施工过程中同样会扰动桩周土体，因此，Li等^[22]考虑土体双向成层特性，土体材料阻尼采用滞回阻尼模型，解析求解了复杂非均质土中管桩纵向振动特性.

综上，考虑桩周土体施工扰动影响的成果多针对实体桩展开，且土体材料阻尼均采用滞回阻尼模型.瞬态激振作用下土阻尼力与应变速率相关，土体滞回阻尼模型会引起概念上的矛盾，此时应用黏性阻尼模型则较合理^[23-28].基于此点考虑，本文基于黏性阻尼模型考虑土体材料阻尼，考虑施工扰动引起的桩周土体径向非均质性，利用多圈层径向复刚度传递，研究考虑施工扰动效应的管桩纵向振动特性.

1 定解问题力学模型建立 1.1 计算简图

基于平面应变模型，对任意激振下黏弹性支承管桩-径向多圈层土体纵向耦合振动特性进行研究，力学简化模型如图 1所示.桩顶作用任意激振力p(t)，桩芯土和桩周土对桩身产生的切应力分别为f₀和f₁，桩长、内半径、外半径、桩身密度、弹性模量和桩底黏弹性支承常数分别为H、r₀、r₁、ρ_p、E_p和δ_p、k_p，桩周第i圈层土的密度、剪切模量、黏性阻尼系数分别为ρ_i、G_i、η_i，桩芯土的密度、剪切模量、黏性阻尼系数分别为ρ₀、G₀、η₀.

基本假定如下：

1) 桩周土和桩芯土为一系列相互独立的薄层，忽略土层间相互作用；

2) 桩周土体分为内部区域和外部区域，内部区域划分任意圈层，各圈层土体均质，土体材料阻尼采用黏性阻尼模型；

3) 各圈层土界面两侧位移连续、应力平衡，且桩土系统振动为小变形，桩土界面完全接触，无脱开和滑移^[29]；

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

1.2 多圈层参数确定

假设内部区域径向厚度为b，内外部区域界面处的半径为r_n+1，内部区域第j圈层土体的剪切模量和黏性阻尼系数G_i、η_i根据该圈层内边界半径r_i，按式(1)和式(2)确定：

$ G\left( r \right) = \left\{ \begin{array}{l} {G_1} \to r = {r_1},\\ {G_{n + 1}}f\left( r \right) \to {r_1} < r < {r_{n + 1}},\\ {G_{n + 1}} \to r \ge {r_{n + 1}}, \end{array} \right. $

(1)

$ \eta \left( r \right) = \left\{ \begin{array}{l} {\eta _1} \to r = {r_1},\\ {\eta _{n + 1}}f\left( r \right) \to {r_1} < r < {r_{n + 1}},\\ {\eta _{n + 1}} \to r \ge {r_{n + 1}}, \end{array} \right. $

(2)

式中：G₁、η₁、G_n+1、η_n+1分别为桩土、内外区域界面处土体剪切模量、黏性阻尼系数，f(r)参照文献[20]采用二次函数描述桩周土内部区域土体性质变化.

1.3 定解问题

设桩周第i圈层土体位移为u_i(r，t)，桩芯土体位移为u₀(r，t)，本文工程背景基于现浇大直径管桩，因此，桩芯土为满布.平面应变条件下桩周、桩芯土体振动方程如下：

桩周土体

$ \begin{array}{*{20}{c}} {{G_i}\frac{{{\partial ^2}{u_i}\left( {r,t} \right)}}{{\partial {r^2}}} + {\eta _i}\frac{{{\partial ^3}{u_i}\left( {r,t} \right)}}{{\partial t\partial {r^2}}} + \frac{{{G_i}}}{r}\frac{{\partial {u_i}\left( {r,t} \right)}}{{\partial r}} + }\\ {\frac{{{\eta _i}}}{r}\frac{{{\partial ^2}{u_i}\left( {r,t} \right)}}{{\partial t\partial r}} = {\rho _i}\frac{{{\partial ^2}{u_i}\left( {r,t} \right)}}{{\partial {t^2}}}.} \end{array} $

(3)

桩芯土体

$ \begin{array}{*{20}{c}} {{G_0}\frac{{{\partial ^2}{u_0}\left( {r,t} \right)}}{{\partial {r^2}}} + {\eta _0}\frac{{{\partial ^3}{u_0}\left( {r,t} \right)}}{{\partial t\partial {r^2}}} + \frac{{{G_0}}}{r}\frac{{\partial {u_0}\left( {r,t} \right)}}{{\partial r}} + }\\ {\frac{{{\eta _0}}}{r}\frac{{{\partial ^2}{u_0}\left( {r,t} \right)}}{{\partial t\partial r}} = {\rho _0}\frac{{{\partial ^2}{u_0}\left( {r,t} \right)}}{{\partial {t^2}}}.} \end{array} $

(4)

设桩身位移为u_p(z，t)，则符合平截面假定的桩身纵向振动基本方程为

$ \frac{{{\partial ^2}{u_{\rm{p}}}\left( {z,t} \right)}}{{\partial {z^2}}} - \frac{{2{\rm{ \mathit{ π} }}{r_0}{f_0}}}{{{E_{\rm{p}}}{A_{\rm{p}}}}} - \frac{{2{\rm{ \mathit{ π} }}{r_1}{f_1}}}{{{E_{\rm{p}}}{A_{\rm{p}}}}} = \frac{{{\rho _{\rm{p}}}}}{{{E_{\rm{p}}}}}\frac{{{\partial ^2}{u_{\rm{p}}}\left( {z,t} \right)}}{{\partial {t^2}}}. $

(5)

边界条件

1) 桩芯土

当r=0时，位移为有限值，即

$ \mathop {\lim }\limits_{r \to 0} {u_0}\left( {r,t} \right) = 有限值. $

(6)

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

$ {u_0}\left( {{r_0},t} \right) = {u_{\rm{p}}}\left( {{r_0},t} \right), $

(7)

$ {f_0} = {\tau _0}\left( r \right)\left| {_{r = {r_0}}} \right.. $

(8)

其中剪应力顺时针为正.

2) 桩周土

当r=∞时，位移为零，即

$ \mathop {\lim }\limits_{r \to \infty } {u_{n + 1}}\left( {r,t} \right) = 0. $

(9)

式中u_n+1(r, t)代表外部区域土体位移.

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

$ {u_1}\left( {{r_1},t} \right) = {u_{\rm{p}}}\left( {{r_1},t} \right), $

(10)

$ {f_1} = - {\tau _1}\left( r \right)\left| {_{r = {r_1}}} \right.. $

(11)

3) 桩身

桩顶作用力为p(t)，即

$ {E_{\rm{p}}}{A_{\rm{p}}}\frac{{\partial {u_{\rm{p}}}\left( {z,t} \right)}}{{\partial z}}\left| {_{z = 0}} \right. = p\left( t \right). $

(12)

桩端处边界条件

$ {E_{\rm{p}}}\frac{{\partial {u_{\rm{p}}}\left( {z,t} \right)}}{{\partial z}}\left| {_{z = H}} \right. = - \left( {{k_{\rm{p}}}{u_{\rm{p}}}\left( {z,t} \right) + {\delta _{\rm{p}}}\frac{{\partial {u_{\rm{p}}}\left( {z,t} \right)}}{{\partial t}}} \right). $

(13)

2 定解问题求解 2.1 桩周土振动方程的求解

对方程(3)进行Laplace变换得

$ \begin{array}{*{20}{c}} {{G_i}\frac{{{\partial ^2}{U_i}\left( {r,s} \right)}}{{\partial {r^2}}} + {\eta _i}s\frac{{{\partial ^2}{U_i}\left( {r,s} \right)}}{{\partial {r^2}}} + \frac{{{G_i}}}{r}\frac{{\partial {U_i}\left( {r,s} \right)}}{{\partial r}} + }\\ {\frac{{{\eta _i}s}}{r}\frac{{\partial {U_i}\left( {r,s} \right)}}{{\partial r}} = {\rho _i}{s^2}{U_i}\left( {r,s} \right).} \end{array} $

(14)

式中:U_i(r, s)为u_i(r, t)的Laplace变换，进一步整理式(14)可得

$ \frac{{{\partial ^2}{U_i}\left( {r,s} \right)}}{{\partial {r^2}}} + \frac{1}{r}\frac{{\partial {U_i}\left( {r,s} \right)}}{{\partial r}} = q_i^2{U_i}\left( {r,s} \right). $

(15)

式中$q_i^2 = \frac{{{\rho _i}{s^2}}}{{{G_i} + {\eta _i}s}}$.

方程(15)的通解为

$ {U_i}\left( {r,s} \right) = {A_i}{K_0}\left( {{q_i}r} \right) + {B_i}{I_0}\left( {{q_i}r} \right). $

(16)

式中：I₀(q, r)、K₀(q_ir)分别为零阶第一类、第二类修正Bessel函数，A_i、B_i为待定系数.

对式(9)进行Laplace变换，代入式(16)可得B_n+1=0，因此，有

$ {U_{n + 1}}\left( {r,s} \right) = {A_{n + 1}}{K_0}\left( {{q_{n + 1}}r} \right). $

(17)

则外部区域土体任意点的竖向剪切应力为

$ \begin{array}{l} {\tau _{n + 1}} = \left( {{G_{n + 1}} + {\eta _{n + 1}}s} \right)\frac{{\partial {U_{n + 1}}\left( {r,s} \right)}}{{\partial r}} = \\ \;\;\;\;\;\;\;\;\;\; - \left( {{G_{n + 1}} + {\eta _{n + 1}}s} \right){q_{n + 1}}{A_{n + 1}}{K_1}\left( {{q_{n + 1}}r} \right). \end{array} $

(18)

内、外部区域分界面处，内部区域受到外部区域的竖向剪切刚度为

$ \begin{array}{l} K{K_{n + 1}} = - \frac{{2{\rm{ \mathit{ π} }}{r_{n + 1}}{\tau _{n + 1}}\left( {{r_{n + 1}}} \right)}}{{{U_{n + 1}}\left( {{r_{n + 1}},s} \right)}} = \\ \;\;\;\;\;\;\;2{\rm{ \mathit{ π} }}{r_{n + 1}}\left( {{G_{n + 1}} + {\eta _{n + 1}}s} \right){q_{n + 1}}\frac{{{K_1}\left( {{q_{n + 1}}{r_{n + 1}}} \right)}}{{{K_0}\left( {{q_{n + 1}}{r_{n + 1}}} \right)}}. \end{array} $

(19)

内部区域圈层i中任意点处的竖向剪切力为

$ \begin{array}{l} {\tau _i} = \left( {{G_i} + {\eta _i}s} \right)\frac{{\partial {U_i}\left( {r,s} \right)}}{{\partial r}} = \\ \;\;\;\;\;\;\;\; - \left( {{G_i} + {\eta _i}s} \right){q_i}\left[ {{A_i}{K_1}\left( {{q_i}r} \right) - {B_i}{I_1}\left( {{q_i}r} \right)} \right]. \end{array} $

(20)

则i圈层的外边界(r=r_i+1)和内边界(r=r_i)处的竖向剪切刚度为

$ \begin{array}{l} K{K_{i + 1}} = - \frac{{2{\rm{ \mathit{ π} }}{r_{i + 1}}{\tau _i}\left( {{r_{i + 1}}} \right)}}{{{U_i}\left( {{r_{i + 1}},s} \right)}} = \\ 2{\rm{ \mathit{ π} }}{r_{i + 1}}{q_i}\left( {{G_i} + {\eta _i}s} \right)\frac{{{A_i}{K_1}\left( {{q_i}{r_{i + 1}}} \right) - {B_i}{I_1}\left( {{q_i}{r_{i + 1}}} \right)}}{{{A_i}{K_0}\left( {{q_i}{r_{i + 1}}} \right) + {B_i}{I_0}\left( {{q_i}{r_{i + 1}}} \right)}}. \end{array} $

(21)

$ \begin{array}{l} K{K_i} = - \frac{{2{\rm{ \mathit{ π} }}{r_i}{\tau _{i\left( {{r_i}} \right)}}}}{{{U_i}\left( {{r_i},s} \right)}} = \\ \;\;\;\;\;\;\;2{\rm{ \mathit{ π} }}{r_i}{q_i}\left( {{G_i} + {\eta _i}s} \right)\frac{{{A_i}{K_1}\left( {{q_i}{r_i}} \right) - {B_i}{I_1}\left( {{q_i}{r_i}} \right)}}{{{A_i}{K_0}\left( {{q_i}{r_i}} \right) + {B_i}{I_0}\left( {{q_i}{r_i}} \right)}}. \end{array} $

(22)

由此可得基于黏性阻尼的多圈层平面应变模型，土层剪切刚度递推公式

$ K{K_i} = 2{\rm{ \mathit{ π} }}{r_i}{q_i}\left( {{G_i} + {\eta _i}s} \right)\frac{{{C_i} + {E_i}K{K_{i + 1}}}}{{{D_i} + {F_i}K{K_{i + 1}}}}. $

(23)

式中：

$ \begin{array}{l} {C_i} = 2{\rm{ \mathit{ π} }}{r_{i + 1}}{q_i}\left( {{G_i} + {\eta _i}s} \right)\left[ {{I_1}\left( {{q_i}{r_{i + 1}}} \right){K_1}\left( {{q_i}{r_i}} \right) - } \right.\\ \left. {{K_1}\left( {{q_i}{r_{i + 1}}} \right){I_1}\left( {{q_i}{r_i}} \right)} \right], \end{array} $

(24)

$ \begin{array}{l} {D_i} = 2{\rm{ \mathit{ π} }}{r_{i + 1}}{q_i}\left( {{G_i} + {\eta _i}s} \right)\left[ {{I_1}\left( {{q_i}{r_{i + 1}}} \right){K_1}\left( {{q_i}{r_i}} \right) + } \right.\\ \left. {{K_1}\left( {{q_i}{r_{i + 1}}} \right){I_0}\left( {{q_i}{r_i}} \right)} \right], \end{array} $

(25)

$ {E_i} = {I_0}\left( {{q_i}{r_{i + 1}}} \right){K_1}\left( {{q_i}{r_i}} \right) + {K_0}\left( {{q_i}{r_{i + 1}}} \right){I_1}\left( {{q_i}{r_i}} \right), $

(26)

$ {F_i} = {I_0}\left( {{q_i}{r_{i + 1}}} \right){K_0}\left( {{q_i}{r_i}} \right) - {K_0}\left( {{q_i}{r_{i + 1}}} \right){I_0}\left( {{q_i}{r_i}} \right). $

(27)

2.2 桩芯土振动方程求解

对方程(4)进行Laplace变换可得

$ \begin{array}{*{20}{c}} {{G_0}\frac{{{\partial ^2}{U_0}\left( {r,s} \right)}}{{\partial {r^2}}} + {\eta _0}s\frac{{{\partial ^2}{U_0}\left( {r,s} \right)}}{{\partial {r^2}}} + \frac{{{G_0}}}{r}\frac{{\partial {U_0}\left( {r,s} \right)}}{{\partial r}} + }\\ {\frac{{{\eta _0}s}}{r}\frac{{\partial {U_0}\left( {r,s} \right)}}{{\partial r}} = {\rho _0}{s^2}{U_0}\left( {r,s} \right).} \end{array} $

(28)

式中，U₀(r, s)为u₀(r, t)的Laplace变换，进一步整理式(28)可得

$ \frac{{{\partial ^2}{U_0}\left( {r,s} \right)}}{{\partial {r^2}}} + \frac{1}{r}\frac{{\partial {U_0}\left( {r,s} \right)}}{{\partial r}} = q_0^2{U_0}\left( {r,s} \right). $

(29)

式中$q_0^2 = \frac{{{\rho _0}{s^2}}}{{{G_0} + {\eta _0}s}}$.

方程(29)的通解为

$ {U_0}\left( {r,s} \right) = {A_0}{K_0}\left( {{q_0}r} \right) + {B_0}{I_0}\left( {{q_0}r} \right). $

(30)

对式(6)进行Laplace变换，代入式(30)可得A₀=0，因此，有

$ {U_0}\left( {r,s} \right) = {B_0}{I_0}\left( {{q_0}r} \right). $

(31)

由式(7)和(31)可得管桩内壁受到桩芯土体的剪切刚度为

$ K{K_0} = - \frac{{2{\rm{ \mathit{ π} }}{r_0}{\tau _0}\left( {{r_0}} \right)}}{{{U_{\rm{p}}}}} = - 2{\rm{ \mathit{ π} }}{r_0}\left( {{G_0} + {\eta _0}s} \right){q_0}\frac{{{I_1}\left( {{q_0}{r_0}} \right)}}{{{I_0}\left( {{q_0}{r_0}} \right)}}. $

(32)

2.3 桩振动方程求解

对式(5)进行Laplace变换，并将由式(23)递推所得桩周最内层土体与桩接触面上剪切刚度KK₁和式(32)代入可得

$ \frac{{{\partial ^2}{U_{\rm{p}}}\left( {z,s} \right)}}{{\partial {z^2}}} - {\alpha ^2}{U_{\rm{p}}}\left( {z,s} \right) = 0. $

(33)

式中：${\alpha ^2} = \frac{{K{K_1}}}{{{E_{\rm{p}}}{A_{\rm{p}}}}} - \frac{{K{K_0}}}{{{E_{\rm{p}}}{A_{\rm{p}}}}} + \frac{{{\rho _{\rm{p}}}}}{{{E_{\rm{p}}}}}{s^2}$，U_p(r, s)为u_p(r, t)的Laplace变换.

方程(33)的通解为

$ {U_{\rm{p}}}\left( {z,s} \right) = C{{\rm{e}}^{\alpha z}} + D{{\rm{e}}^{ - \alpha z}}. $

(34)

式中C、D为由边界条件确定的待定系数.

进一步对式(12)、(13)进行Laplace变换，并将式(34)代入后可得

$ C = \frac{{\xi P\left( s \right)}}{{\alpha \left( {\xi - 1} \right){E_{\rm{p}}}{A_{\rm{p}}}}}, $

(35)

$ D = \frac{{P\left( s \right)}}{{\alpha \left( {\xi - 1} \right){E_{\rm{p}}}{A_{\rm{p}}}}}. $

(36)

式中$\xi = \frac{{\alpha {E_{\rm{p}}} - ({k_{\rm{p}}} + {\delta _{\rm{p}}}s)}}{{\alpha {E_{\rm{p}}} + ({k_{\rm{p}}} + {\delta _{\rm{p}}}s)}}{{\rm{e}}^{ - 2\alpha H}}$.

将式(35)、(36)代入(34)可得

$ {U_{\rm{p}}}\left( {z,s} \right) = \frac{{\left( {\xi {{\rm{e}}^{\alpha z}} + {{\rm{e}}^{ - \alpha z}}} \right)}}{{\alpha \left( {\xi - 1} \right){E_{\rm{p}}}{A_{\rm{p}}}}}P\left( s \right). $

(37)

式中:令s=iω，则Laplace变换相当于单边的Fourier变换，因此, 可以得到管桩位移频域响应函数为U_p(z, iω).由此可得桩顶复动刚度

$ {K_{\rm{d}}} = \frac{{P\left( {{\rm{i}}\omega } \right)}}{{{U_{\rm{p}}}\left( {z,{\rm{i}}\omega } \right)}}\left| {_{z = 0}} \right. = \frac{{\alpha \left( {\xi - 1} \right){E_{\rm{p}}}{A_{\rm{p}}}}}{{\left( {\xi + 1} \right)}} = \frac{{{E_{\rm{p}}}{A_{\rm{p}}}}}{H}{{K'}_{\rm{d}}}, $

(38)

$ {{K'}_{\rm{d}}} = \frac{{\bar \alpha \left( {\xi - 1} \right)}}{{\left( {\xi + 1} \right)}}. $

(39)

式中：α =αH，P(iω)为p(t)的傅里叶变换，K_d^′为无量纲复刚度.令K_d^′=K_r+iK_i，其中K_r代表动刚度，K_i代表动阻尼.

3 算例分析

基于图 1所示径向非均质黏弹性土体中管桩耦合振动力学模型，采用前述推导求解所得管桩纵向振动动力阻抗解析解答，如无特殊说明，具体土层参数取值如下：桩长6 m，外直径1.0 m，壁厚0.12 m，桩身混凝土弹性模量25 GPa，密度2 500 kg/m³；桩周土和桩芯土参数相同，剪切模量为5 MPa，密度为2 000 kg/m³，黏性阻尼系数为η=10 kN · s/m²，桩底刚度系数k_p=1 000 kN/m³，阻尼系数δ_p=100 kN/m³.

当桩周土体径向圈层数大于20时^[11]，径向圈层的数目对计算结果的影响是可以忽略的.因此，随后的分析将桩周围土体的内部区域划分20个圈层.假定外部区域剪切波速${v_{n + 1}}({v_{n + 1}} = \sqrt {{G_{n + 1}}/{\rho _{n + 1}}} )$至第1圈层剪切波速v₁呈线性变化，并定义施工扰动系数β=v₁/v_n+1以描述施工扰动影响.桩周土体的黏性阻尼系数从外区到第1圈呈现二次函数变化，并且$\beta = \sqrt {{\eta _1}/{\eta _{n + 1}}} $.β取值不同代表施工扰动程度不同，β>1代表施工硬化，β < 1代表施工软化，取η_n+1=5 kN · s/m³，v_n+1=50 m/s，b=r₁，β=2.0，施工扰动引起桩周、桩芯土参数变化示意如图 2所示.

3.1 解析解答合理性验证

为了验证本文所推导的径向非均质土中管桩纵向振动动力阻抗解析解答的合理性，桩周围非均质土退化为均质与文献[17]中给出的解析解比较见图 3.此外，管桩退化为实体桩(阻尼系数取0)，与文献[11]已有解析解进行对比验证，如图 3，4所示.图中无量纲频率θ为以特征频率为单位的频率计量，其中特征频率为ω_c=v_P/H，v_P=$\sqrt {{E_{\rm{P}}}/{\rho _{\rm{P}}}} $为桩身压缩波速，则θ=ω/ω_c=ωH/v_P.由图 4可知，本文所退化解答曲线与文献[17, 11]中结果吻合.

3.2 管桩桩身参数对桩顶动力阻抗影响分析

图 5为管桩长径比对桩顶动力阻抗的影响.桩长径比一般控制在50以内，此处取长径比10，15，20，25，30进行对比分析.可以看出，管桩长径比对桩顶动力阻抗曲线影响较大.管桩内外壁侧摩阻力随管桩长径比增加而增大，这就使得桩顶动力阻抗曲线振幅、共振频率均有所减小.桩长径比变化对桩顶动力阻抗影响呈现衰减趋势，当其达到一定幅值后此种影响效应趋于稳定.

图 6为管桩顶动力阻抗随管桩内外径变化情况.可以看出，动力阻抗曲线的振幅、共振频率均随着管桩内外径比减小而增大.且当r₀/r₁=0(实体桩)时，桩顶动力阻抗曲线振幅最大，由此可知桩芯土能对管桩桩体起到减振作用.

图 7为桩弹性模量对桩顶动力阻抗曲线的影响.可以看出，随着管桩弹性模量的增加，管桩桩顶动力阻抗曲线振幅均逐步增大，且管桩弹性模量对管桩动刚度影响较动阻尼更为显著.不同地，管桩弹性模量增大仅引起桩顶动阻尼共振频率小幅减小.管桩弹性模量变化对桩顶动刚度和动阻尼的影响程度逐步衰减，当其达到一定幅值后此种影响趋于稳定.

3.3 桩周土体施工扰动对桩顶动力阻抗影响分析

图 8为桩周土体黏性阻尼系数对桩顶动力阻抗的影响.可以看出，桩顶动刚度和动阻尼振幅和共振频率均随桩周土体黏性阻尼系数增加而显著减小，且桩周土体黏性阻尼系数变化对桩顶动刚度和动阻尼的影响程度随其增加而逐步衰减.

桩身施工扰动会引起原土质较松散或不够密实的砂土、粉土在靠近桩身的内部区域土体硬化.不同地，原土质为密实的砂土、粉土以及高灵敏度黏土、软土时，靠近桩身的内部区域土体在桩身施工扰动影响下则产生软化效应.基于此，对管桩施工扰动引起的桩周土体软化、硬化对桩顶动力阻抗的影响进行分析.

图 9，10分别为桩周土软化、硬化程度对桩顶动力阻抗曲线的影响.可以看出，桩周土的软化和硬化程度对桩顶动力阻抗曲线的幅值有明显的影响，即动力阻抗曲线振幅随桩周土体软化程度升高而增大，随桩周土体硬化程度升高而减小.但桩周土软化、硬化程度对桩顶动力阻抗共振频率影响较小.

图 11，12分别为土体软化和硬化范围对桩顶动力阻抗的影响.可以看出，土体软化和硬化范围对桩顶动力阻抗曲线的幅值有明显影响，振幅随桩周土软化范围增加而增大，随桩周土体硬化范围增加而减小.但与实体桩不同，此种影响由于桩芯土的存在仅在近桩小范围内较为明显，超出该范围即大幅衰减至稳定.

4 结论

1) 桩顶动力阻抗曲线振幅、共振频率均随管桩长径比的增加而减小，且桩长径比变化对桩顶动力阻抗特性影响随其增加而呈现衰减趋势.

2) 桩顶动力阻抗曲线的振幅和共振频率随内外径比的减小而增大，由此可知桩芯土能对管桩桩体起到减振作用.

3) 桩顶动力阻抗曲线振幅随桩身弹性模量的增加而增大，随着桩的弹性模量的增加，这种影响程度逐渐减弱.

4) 管桩桩顶动力阻抗振幅、共振频率随桩周土体黏性阻尼系数的增加均显著减小，且桩周土体黏性阻尼系数变化对桩顶动力阻抗特性影响程度随其增加而逐步衰减.

5) 桩周土体的软(硬)化程度越高，桩顶动力阻抗曲线振幅越大(小).相对而言，桩周土软(硬)化程度对桩顶动力阻抗共振频率影响程度较小.

6) 管桩桩顶动力阻抗振幅水平随桩周土软(硬)化范围增加而增大(减小)，此种影响由于桩芯土的存在仅在近桩小范围内较为明显，超出该范围即大幅衰减至稳定.不同地，桩周土施工扰动范围改变对桩顶动力阻抗共振频率影响很小.