纵向成层土与桩基相互作用体系耦合振动特性研究一直都是岩土工程领域的热点问题^[1-2].桩周土体在桩基施工过程中会产生沿桩基径向的不均匀性^[3]，而此种径向非均质效应对桩基纵向振动特性的影响也受到国内外学者的广泛关注.

Novak等^[4]利用滞回阻尼模型考虑土体黏性，并采用平面应变理论建立了存在单层内部扰动区域的桩身-土体动力相互作用模型，初步分析了土体径向非均质效应对桩身纵向振动特性的影响.在此基础上，EI Naggar^[5]为更合理地考虑桩周土的径向不均匀性，将内部区域进一步划分为多个圈层，分析了施工扰动效应对桩基纵向振动特性的影响.王奎华等^[6-7]指出EI Naggar模型与实测曲线存在较大差别，Wang等^[8-9]基于此对EI Naggar模型进行修正，土体材料采用滞回阻尼模型，提出了较为严格的平面应变径向多圈层模型，进行了桩顶振动特性频域解析和时域半解析研究.

以上研究对象均为实心桩，然而管桩由于存在桩芯土，与实心桩振动特性必然存在一定差异性.刘汉龙等^[10]基于桩身三维模型，忽略土体对管桩的约束作用，研究了管桩桩身平截面假定对桩顶动力响应的影响深度.丁选明等^[11]在此基础上将桩周土及桩芯土对管桩的作用简化为弹簧-阻尼器模型，分析了桩顶速度反射曲线的影响因素.此外，丁选明等^[12]利用简化的平面应变模型考虑土体作用，郑长杰等^[13]基于三维轴对称模型，Zheng等^[14]考虑土体真三维波动效应，对均质黏弹性地基中管桩纵向振动频域特性进行理论研究.为考虑管桩在施工过程中引起土体径向非均质效应，Li等^[15]考虑土层复杂非均匀特性，采用滞回阻尼模型，对管桩-土耦合纵向振动特性进行解析研究.

综上，考虑管桩施工扰动引起径向非均质性的研究较少，且在考虑此种效应时土体材料阻尼均采用滞回阻尼模型(与频率及应变速率无关)，其在解决非谐和激振问题时不考虑土体材料阻尼随频率及应变速率的变化会在概念上引起矛盾，此时采用黏性阻尼模型^[16-17]更为合理.崔春义等^[18]基于黏性阻尼土体模型，考虑土层径向非均质效应，对考虑施工扰动效应的管桩-土耦合纵向振动特性进行研究.因土体天然沉积作用使得其沿纵向成层分布，且桩基施工过程中的扰动会引起土体的径向非均质，为考虑土体这种复杂非均质特性，同时避免土体滞回阻尼模型在应用非谐和激振问题时产生的矛盾，提出一种基于土体黏性阻尼的双向非均质土体-管桩动力相互作用模型，通过解析方法求解了管桩桩顶受任意荷载作用下的纵向振动特性，为管桩抗振防振设计提供理论指导.

1 桩土动力相互作用模型 1.1 力学简化模型

根据土体纵向成层特性，将桩周土-管桩-桩芯土相互作用体系分成m个层段，各层段自桩身底部依次编号为1, 2, …, i, …, m，厚度为l₁, l₂, …, l_i, …, l_m，顶部埋深为h₁, h₂, …，h_i，…，h_m.管桩面积、弹性模量、密度、内径和外径分别为A_i^P、E_i^P、ρ_i^P、r_i0和r_i1，桩底支承阻尼系数和刚度系数分别为δ_p和k_p.第i层段桩芯土密度、剪切模量和黏性阻尼系数分别为ρ_i0、G_i0和η_i0.第i层段桩周土体厚度为b_i的内部扰动区域划分为n个径向圈层，其中第j圈层桩周土的密度、剪切模量和黏性阻尼系数分别为ρ_ij、G_ij和η_ij，其值按式(1)和式(2)确定：

$ {G_{ij}}(r) = \left\{ {\begin{array}{*{20}{c}} {{G_{i1}},}&{r = {r_{i1}},}\\ {{G_{i(n + 1)}} \times {f_i}(r),}&{{r_{i1}} < r < {r_{i(n + 1)}},}\\ {{G_{i(n + 1)}},}&{r \ge {r_{i(n + 1)}},} \end{array}} \right. $

(1)

$ {\eta _{ij}}(r) = \left\{ {\begin{array}{*{20}{c}} {{\eta _{i1}},}&{r = {r_{i1}},}\\ {{\eta _{i(n + 1)}} \times {f_i}(r),}&{{r_{i1}} < r < {r_{i(n + 1)}},}\\ {{\eta _{i(n + 1)}},}&{r \ge {r_{i(n + 1)}}.} \end{array}} \right. $

(2)

式中f_i(r)为二次函数形式^[8].

桩土耦合振动体系力学简化模型如图 1所示.本文基本假定参照文献[16]相关描述.

1.2 定解问题

平面应变条件下基于土体黏性阻尼模型的桩周土和桩芯土纵向振动控制方程为

$ \begin{array}{*{20}{c}} {{G_{ij}}\frac{{{\partial ^2}w_{ij}^1(r,t)}}{{\partial {r^2}}} + {\eta _{ij}}\frac{{{\partial ^3}w_{ij}^1(r,t)}}{{\partial t\partial {r^2}}} + \frac{{{G_{ij}}}}{r}\frac{{{\partial ^2}w_{ij}^1(r,t)}}{{\partial {t^2}}} + }\\ {\frac{{{\eta _{ij}}}}{r}\frac{{{\partial ^3}w_{ij}^1(r,t)}}{{\partial t\partial r}} = {\rho _{ij}}\frac{{{\partial ^2}w_{ij}^1(r,t)}}{{\partial {t^2}}},} \end{array} $

(3)

$ \begin{array}{*{20}{c}} {{G_{i0}}\frac{{{\partial ^2}w_i^0(r,t)}}{{\partial {r^2}}} + {\eta _{i0}}\frac{{{\partial ^3}w_i^0(r,t)}}{{\partial t\partial {r^2}}} + \frac{{{G_{i0}}}}{r}\frac{{\partial w_i^0(r,t)}}{{\partial r}} + }\\ {\frac{{{\eta _{i0}}}}{r}\frac{{{\partial ^3}w_i^0(r,t)}}{{\partial t\partial r}} = {\rho _{i0}}\frac{{{\partial ^2}w_i^0(r,t)}}{{\partial {t^2}}},} \end{array} $

(4)

式中w_i⁰、w_ij¹分别为桩芯和桩周土体位移.

设第i层段桩身位移为w_i^p，则桩身纵向振动方程为

$ \frac{{{\partial ^2}w_i^p}}{{\partial {z^2}}} - \frac{{2\pi {r_{i1}}f_i^{{S_1}}}}{{E_i^PA_i^P}} - \frac{{2\pi {r_{i0}}f_i^{{S_0}}}}{{E_i^PA_i^P}} = \frac{{\rho _i^P}}{{E_i^P}}\frac{{{\partial ^2}w_i^p}}{{\partial {t^2}}}, $

(5)

式中f_i^S₀、f_i^S₁分别为第i层段桩芯土和桩周土对桩身产生的切应力.

桩芯土边界条件为

$ {\mathop {{\rm{lim}}}\limits_{r \to 0} w_i^0(r,t) = {\rm{ 有限值 }},} $

(6)

$ {w_i^0({r_{i0}},t) = w_i^p({r_{i0}},t),} $

(7)

$ {f_i^{{S_0}} = \tau _i^{{S_0}}({r_{i0}}).} $

(8)

桩周土边界条件为

$ {\mathop {{\rm{lim}}}\limits_{r \to \infty } w_{i(n + 1)}^1(r,t) = 0,} $

(9)

$ {w_i^1({r_{i1}},t) = w_i^p({r_{i1}},t),} $

(10)

$ f_i^{{S_1}} = \tau _i^{{S_1}}({r_{i1}}). $

(11)

桩身边界条件为

$ {{{\left. {E_m^PA_m^P\frac{{\partial w_m^P}}{{\partial z}}} \right|}_{z = 0}} = - p(t),} $

(12)

$ {E_1^P\frac{{\partial w_1^P}}{{\partial z}} + {{\left. {({k_{\rm{p}}}w_1^P + {\delta _{\rm{p}}}\frac{{\partial w_1^P}}{{\partial t}})} \right|}_{z = H}} = 0.} $

(13)

2 解析推导过程 2.1 桩周土解析解

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

$ \begin{array}{l} {G_{ij}}\frac{{{\partial ^2}}}{{\partial {r^2}}}W_{ij}^1(r,s) + {\eta _{ij}}s\frac{{{\partial ^2}}}{{\partial {r^2}}}W_{ij}^1(r,s) + \frac{{{G_{ij}}}}{r}\frac{\partial }{{\partial r}}W_{ij}^1(r,s) + \\ \frac{{{\eta _{ij}}s}}{r}\frac{\partial }{{\partial r}}W_{ij}^1(r,s) = {\rho _{ij}}{s^2}W_{ij}^1(r,s), \end{array} $

(14)

式中W_ij¹(r, s)为w_ij¹(r, t)的Laplace变换.

进一步对式(14)进行整理可得

$ \frac{{{\partial ^2}}}{{\partial {r^2}}}W_{ij}^1 + \frac{1}{r}\frac{\partial }{{\partial r}}W_{ij}^1 = {(q_{ij}^{{S_1}})^2}W_{ij}^1, $

(15)

式中$q_{i j}^{1}=\left(\frac{\rho_{i j} s^{2}}{G_{i j}+\eta_{i j} s}\right)^{1 / 2} $.

求解可得方程(15)的通解为

$ W_{ij}^1 = A_{ij}^1{K_0}(q_{ij}^1r) + B_{ij}^1{I_0}(q_{ij}^1r). $

(16)

式中：I₀(q_ij¹r)、K₀(q_ij¹r)分别为零阶第一类、第二类修正Bessel函数，A_ij¹、B_ij¹为待定系数.

进一步利用桩周土体边界条件及各圈层的连续条件可得土层剪切刚度递推公式为

$ FF_{ij}^1 = 2\pi {r_{ij}}({G_{ij}} + {\eta _{ij}}s)q_{ij}^1\frac{{C_{ij}^1 + M_{ij}^1FF_{i(j + 1)}^1}}{{D_{ij}^1 + N_{ij}^1FF_{i(j + 1)}^1}}. $

(17)

式中：FF_ij¹、FF_i(n+1)¹分别为第i层段j圈层中外边界(r=r_i(j+1))和内边界(r=r_ij)处的剪切刚度.

$ \begin{array}{*{20}{l}} {C_{ij}^1 = 2\pi {r_{i(j + 1)}}({G_{ij}} + {\eta _{ij}}s)q_{ij}^1[{I_1}(q_{ij}^1{r_{i(j + 1)}}){K_1}(q_{ij}^1{r_{ij}}) - }\\ {{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {K_1}(q_{ij}^1{r_{i(j + 1)}}){I_1}(q_{ij}^1{r_{ij}})],} \end{array} $

$ \begin{array}{*{20}{l}} {D_{ij}^1 = 2\pi {r_{i(j + 1)}}({G_{ij}} + {\eta _{ij}}s)q_{ij}^1[{I_1}(q_{ij}^1{r_{i(j + 1)}}){K_0}(q_{ij}^1{r_{ij}}) + }\\ {{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {K_1}(q_{ij}^1{r_{i(j + 1)}}){I_0}(q_{ij}^1{r_{ij}})],} \end{array} $

$ {M_{ij}^1 = {K_0}(q_{ij}^1{r_{i(j + 1)}}){I_1}(q_{ij}^1{r_{ij}}) + {K_1}(q_{ij}^1{r_{ij}}){I_0}(q_{ij}^1{r_{i(j + 1)}}),} $

$ {N_{ij}^1 = - {K_0}(q_{ij}^1{r_{i(j + 1)}}){I_0}(q_{ij}^1{r_{ij}}) + {K_0}(q_{ij}^1{r_{ij}}){I_0}(q_{ij}^1{r_{i(j + 1)}}).} $

2.2 桩芯土解析解

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

$ \begin{array}{*{20}{l}} {{G_{i0}}\frac{{{\partial ^2}}}{{\partial {r^2}}}W_i^0(r,s) + {\eta _{i0}}s\frac{{{\partial ^2}}}{{\partial {r^2}}}W_i^0(r,s) + \frac{{{G_{i0}}\partial }}{r}\frac{\partial }{{\partial r}}W_i^0(r,s) + }\\ {\frac{{{\eta _{i0}}s}}{r}\frac{\partial }{{\partial r}}W_i^0(r,s) = {\rho _{i0}}{s^2}W_i^0(r,s),} \end{array} $

(18)

式中W_i⁰(r, s)为w_i⁰(r, t)的Laplace变换.

进一步对式(18)进行整理可得

$ \frac{{{\partial ^2}}}{{\partial {r^2}}}W_i^0 + \frac{1}{r}\frac{\partial }{{\partial r}}W_i^0 = {(q_i^0)^2}W_i^0, $

(19)

式中$ q_{i}^{0}=\left(\frac{\rho_{i 0} s^{2}}{G_{i 0}+\eta_{i 0} s}\right)^{1 / 2}$.

对方程(19)进行求解并利用桩芯土体边界条件可得管桩内壁受到桩芯土体的剪切刚度为

$ \begin{array}{*{20}{l}} {FF_i^0 = - \frac{{2\pi {r_{i0}}\tau _i^{{S_0}}({r_0})}}{{W_i^P}} = }\\ {{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} - 2\pi {r_{i0}}q_i^0({G_{i0}} + {\eta _{i0}}s)\frac{{{I_1}(q_i^0{r_{i0}})}}{{{I_0}(q_i^0{r_{i0}})}}.} \end{array} $

(20)

2.3 桩振动方程求解

对式(5)进行Laplace变换，并将由式(17)递推所得FF_i1¹和式(20)计算所得FF_i⁰代入可得

$ \frac{{{\partial ^2}W_i^P}}{{\partial {z^2}}} - \alpha _i^2W_i^P = 0. $

(21)

式中：$ \alpha_{i}^{2}=\frac{\rho_{i}^{P} A_{i}^{\mathrm{p}} s^{2}+F F_{i 1}^{1}-F F_{i}^{0}}{A_{i}^{P} E_{i}^{P}}, W_{i}^{P}(z, s)$为$ w_{i}^{P}(z, t)$的Laplace变换.

求解可得方程(21)的通解为

$ W_i^P = M_i^P{{\rm{e}}^{{{\bar \alpha }_i}z/{l_i}}} + N_i^P{{\rm{e}}^{ - {{\bar \alpha }_i}z/{l_i}}}. $

(22)

式中：$ {\bar \alpha _i} = {\alpha _i}{l_i}, M_i^P、N_i^P$为待定系数.

进一步利用阻抗函数传递性可得管桩桩顶动力阻抗函数$Z_{m}^{P}(z, s) $为

$ Z_m^P{|_{z = {h_m} = 0}} = \frac{{ - E_m^PA_m^P{{\bar \alpha }_m}({\beta _m} - 1)}}{{{l_m}({\beta _m} + 1)}} = \frac{{ - E_m^PA_m^P}}{{{l_m}}}Z_m^{{P^\prime }}. $

(23)

式中：$ Z_m^{{P^\prime }} = \frac{{{{\bar \alpha }_m}\left( {{\beta _m} - 1} \right)}}{{{\beta _m} + 1}}$为无量纲桩顶阻抗，$ {\beta _m}{\rm{ = }}\left( {E_m^PA_m^P{{\bar \alpha }_m} - Z_{m - 1}^P{l_m}} \right){{\rm{e}}^{ - 2{{\bar \alpha }_m}{h_{m - 1}}/{l_m}}}/E_m^PA_m^P{{\bar \alpha }_m} + Z_{m - 1}^P{l_m}$，令$Z_{m}^{P^{\prime}}=K_{r}+\mathrm{i} K_{i}, \quad K_{r} $代表桩顶动刚度，K_i代表桩顶动阻尼.

3 算例分析

参照文献[6]取径向圈层数n=20，纵向层段取m=5，如无特殊说明，管桩具体参数取值如下：r_i0=0.38 m，r_i1=0.50 m，ρ_i^P=2 500 kg/m³，E_i^P=25 GPa，H=6 m，k_p=1 000 kN/m³，δ_p=100 kN/m³.桩周各圈层土体相关参数按如下方式取值：$ {\rho _{ij}} = 2000{\rm{kg}}/{{\rm{m}}^3}, {b_i} = {r_{1i}}, {q_i} = {V_{i1}}/{V_{i(n + 1)}} = \sqrt {{\eta _{i1}}/{\eta _{i(n + 1)}}} $为施工扰动系数，具体取V_i(n+1)=50 m/s，η_i(n+1)=10 kNs/m²，q_i=2.0.

3.1 解析解对比验证

为了验证本文所提模型及相关解析解答的合理性，分别将本文推导所得桩顶动力阻抗解析解答退化到土体纵向均质和实体桩情况，与文献[18]和文献[6]已有解进行对比，如图 2和3所示.图中$\theta {\rm{ = }}\omega H/{V_{\rm{P}}}, {V_{\rm{P}}} = \sqrt {{E^P}/{\rho ^P}} $，综合图 2和3可以发现，本文所得解析解退化后与已有解答吻合良好.

3.2 参数化分析

桩周土纵向成层特性对桩顶动力阻抗影响如图 4所示.假设桩周存在一个参数不同的某纵向夹层，厚度为1 m，距桩顶距离为3.8 m.定义λ为该夹层土体剪切波速与其余土层剪切波速的比值，λ < 1为软弱夹层，λ>1为硬夹层.由图 4可知，软硬夹层的存在仅对一定低频范围内的振动幅值有较明显影响，其对共振频率的影响则相对较小.剪切波速比越小(大)，夹层越软(硬)，共振频率处对应的动力阻抗振幅水平越高(低).当夹层剪切波速比为1.8时，动力阻抗共振幅值减小约25%，当夹层剪切波速比为0.6时，动力阻抗共振幅值增大约19%.

众所周知，不同深度范围施工扰动引起的桩周土径向软化程度亦不同.桩周土纵向第i层段施工扰动系数q_i随施工软化程度、硬化程度工况变化如表 1所示(i=1, 2, …, 5).

图 5，6分别为施工扰动引起管桩桩周土径向软(硬)化程度对桩顶复刚度曲线的影响.综合图 5，6和表 1可知，桩周土软(硬)化程度对管桩桩顶复刚度曲线的共振幅值及共振频率均有显著影响.具体地，桩顶动力阻抗曲线的振动幅值和共振频率均随土体径向软化程度的增强而增大，随土体径向硬化程度的增强而降低.当桩周土软化到未扰动土的60%时，桩顶动力阻抗曲线共振幅值增大约38%；

当桩周土比未扰动土硬60%时，桩顶动力阻抗曲线共振幅值减小约17%.说明在对管桩纵向振动特性进行分析时，不考虑桩周土体的软化和硬化效应计算所得桩顶纵向振动特性存在较大误差，将对管桩抗振防振设计产生不利影响.

施工扰动范围对桩顶动力阻抗的影响如图 7，8所示.综合两图可以发现，施工扰动范围仅对动力阻抗曲线振动幅值有明显影响，对共振频率的影响可忽略.振动幅值随施工扰动引起土体软(硬)化范围的增大而升高(减小).另外，即使桩周土体扰动范围很小(0.1r₁=0.05 m)，管桩施工仅影响桩身附近几厘米厚度范围内土体性质，其对桩顶动力阻抗特性的影响同样不可忽略，这就说明考虑施工扰动效应的影响对管桩纵向振动特性进行研究的必要性.

4 结论

1) 纵向软硬夹层对一定低频区间内动力阻抗曲线振幅水平影响显著.夹层越软(硬)，共振频率处对应的动力阻抗曲线振幅水平越高(低).当夹层剪切波速比为1.8时，动力阻抗共振幅值减小约25%，当夹层剪切波速比为0.6时，动力阻抗共振幅值增大约19%.

2) 桩周土软(硬)化程度对管桩桩顶复刚度曲线的共振幅值及共振频率均有显著影响，当桩周土软化到未扰动土的60%时，桩顶动力阻抗曲线共振幅值增大约38%；当桩周土比未扰动土硬60%时，桩顶动力阻抗曲线共振幅值减小约17%.说明在对管桩纵向振动特性进行分析时，不考虑桩周土体的软化和硬化效应计算所得桩顶纵向振动特性存在较大误差，将对管桩抗振防振设计产生不利影响.

3) 施工扰动土体软(硬)化范围主要影响管桩桩顶复刚度曲线的共振幅值，即使桩周土体扰动范围很小(0.1r₁=0.05 m)，管桩施工仅影响桩身附近几厘米厚度范围内土体性质，其对桩顶动力阻抗特性的影响同样不可忽略，这就说明考虑施工扰动效应的影响对管桩纵向振动特性进行研究的必要性.