删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

双向非均质黏性阻尼土中管桩纵向振动特性

本站小编 哈尔滨工业大学/2020-12-05

双向非均质黏性阻尼土中管桩纵向振动特性

崔春义,孟坤,梁志孟,赵会杰

(大连海事大学 交通运输工程学院,辽宁 大连 116026)



摘要:

管桩施工会使土体产生扰动,从而引起桩周土体沿桩基径向的不均匀性,为同时考虑桩周土体此种径向非均质效应和纵向成层特性对管桩纵向振动动力阻抗的影响规律,提出基于土体黏性阻尼的多圈层平面应变模型,建立双向非均质土中管桩纵向振动简化分析模型,利用拉普拉斯变换和阻抗函数的传递性,推导得出管桩桩顶动力阻抗解析解答,计算结果表明,纵向软硬夹层对一定低频区间内动力阻抗曲线振幅水平影响显著,夹层越软(硬),共振频率处对应的动力阻抗曲线振幅水平越高(低).施工扰动引起桩周土径向非均匀性对管桩桩顶复刚度曲线影响显著,在对管桩进行纵向振动特性分析时宜采用双向非均质黏性阻尼土-管桩相互作用模型推导所得解析解答更为合理.

关键词:  管桩  径向非均质  纵向成层  复刚度传递  动力阻抗

DOI:10.11918/201904050

分类号:TU443

文献标识码:A

基金项目:国家自然科学基金面上项目(9,1, 51578100);中央高校基本科研业务费专项资金(3132019601);大连理工大学海岸和近海工程国家重点实验室开放课题基金(LP1904); “双一流”建设专项(CXXM2019BS008)



Vertical vibration impedance of pipe pile embedded in radially and vertically inhomogeneous layered soil with viscous damping

CUI Chunyi,MENG Kun,LIANG Zhimeng,ZHAO Huijie

(School of Transportation Engineering, Dalian Maritime University, Dalian 116026, Liaoning, China)

Abstract:

The construction of pipe pile will cause the inhomogeneity of the surrounding soil. By considering the influence of the radial inhomogeneity effect and the vertically layer characteristic on the vertical vibration dynamic impedance of pipe pile, a multi-circle plane strain model based on soil viscous damping was proposed, and a simplified model for vertical vibration analysis of pipe pile in radially and vertically layered soil was established. Based on the Laplace transformation and the transitivity of impedance function, the analytical solution of the dynamic impedance of the pile head was derived. Computational results show that the vertical soft and hard interlayer had a significant effect on the amplitude of the dynamic impedance curve in a certain low frequency range, where the softer (harder) the interlayer was, the higher (lower) the amplitude of the dynamic impedance curve at the resonance frequency became. The radial inhomogeneity caused by construction disturbance had a significant influence on the complex stiffness curve of the pile head. In the analysis of vertical vibration characteristics of pipe piles, the analytical solution derived by the proposed soil-pile interaction model was more reasonable.

Key words:  pipe pile  radially inhomogeneous  vertically layer  complex stiffness transfer  dynamic impedance


崔春义, 孟坤, 梁志孟, 赵会杰. 双向非均质黏性阻尼土中管桩纵向振动特性[J]. 哈尔滨工业大学学报, 2020, 52(11): 113-119. DOI: 10.11918/201904050.
CUI Chunyi, MENG Kun, LIANG Zhimeng, ZHAO Huijie. Vertical vibration impedance of pipe pile embedded in radially and vertically inhomogeneous layered soil with viscous damping[J]. Journal of Harbin Institute of Technology, 2020, 52(11): 113-119. DOI: 10.11918/201904050.
基金项目 国家自然科学基金面上项目(51878109, 51722801, 51578100);中央高校基本科研业务费专项资金(3132019601);大连理工大学海岸和近海工程国家重点实验室开放课题基金(LP1904); “双一流”建设专项(CXXM2019BS008) 作者简介 崔春义(1978—),男,教授,博士生导师 通信作者 崔春义,cuichunyi@dlmu.edu.cn 文章历史 收稿日期: 2019-04-04



Abstract            Full text            Figures/Tables            PDF


双向非均质黏性阻尼土中管桩纵向振动特性
崔春义, 孟坤, 梁志孟, 赵会杰     
大连海事大学 交通运输工程学院,辽宁 大连 116026

收稿日期: 2019-04-04
基金项目: 国家自然科学基金面上项目(51878109, 51722801, 51578100);中央高校基本科研业务费专项资金(3132019601);大连理工大学海岸和近海工程国家重点实验室开放课题基金(LP1904); “双一流”建设专项(CXXM2019BS008)
作者简介: 崔春义(1978—),男,教授,博士生导师
通信作者: 崔春义,cuichunyi@dlmu.edu.cn


摘要: 管桩施工会使土体产生扰动,从而引起桩周土体沿桩基径向的不均匀性,为同时考虑桩周土体此种径向非均质效应和纵向成层特性对管桩纵向振动动力阻抗的影响规律,提出基于土体黏性阻尼的多圈层平面应变模型,建立双向非均质土中管桩纵向振动简化分析模型,利用拉普拉斯变换和阻抗函数的传递性,推导得出管桩桩顶动力阻抗解析解答,计算结果表明, 纵向软硬夹层对一定低频区间内动力阻抗曲线振幅水平影响显著,夹层越软(硬),共振频率处对应的动力阻抗曲线振幅水平越高(低).施工扰动引起桩周土径向非均匀性对管桩桩顶复刚度曲线影响显著,在对管桩进行纵向振动特性分析时宜采用双向非均质黏性阻尼土-管桩相互作用模型推导所得解析解答更为合理.
关键词: 管桩    径向非均质    纵向成层    复刚度传递    动力阻抗    
Vertical vibration impedance of pipe pile embedded in radially and vertically inhomogeneous layered soil with viscous damping
CUI Chunyi, MENG Kun, LIANG Zhimeng, ZHAO Huijie     
School of Transportation Engineering, Dalian Maritime University, Dalian 116026, Liaoning, China



Abstract: The construction of pipe pile will cause the inhomogeneity of the surrounding soil. By considering the influence of the radial inhomogeneity effect and the vertically layer characteristic on the vertical vibration dynamic impedance of pipe pile, a multi-circle plane strain model based on soil viscous damping was proposed, and a simplified model for vertical vibration analysis of pipe pile in radially and vertically layered soil was established. Based on the Laplace transformation and the transitivity of impedance function, the analytical solution of the dynamic impedance of the pile head was derived. Computational results show that the vertical soft and hard interlayer had a significant effect on the amplitude of the dynamic impedance curve in a certain low frequency range, where the softer (harder) the interlayer was, the higher (lower) the amplitude of the dynamic impedance curve at the resonance frequency became. The radial inhomogeneity caused by construction disturbance had a significant influence on the complex stiffness curve of the pile head. In the analysis of vertical vibration characteristics of pipe piles, the analytical solution derived by the proposed soil-pile interaction model was more reasonable.
Keywords: pipe pile    radially inhomogeneous    vertically layer    complex stiffness transfer    dynamic impedance    
纵向成层土与桩基相互作用体系耦合振动特性研究一直都是岩土工程领域的热点问题[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,厚度为l1, l2, …, li, …, lm,顶部埋深为h1, h2, …,hi,…,hm.管桩面积、弹性模量、密度、内径和外径分别为AiPEiPρiPri0ri1,桩底支承阻尼系数和刚度系数分别为δpkp.第i层段桩芯土密度、剪切模量和黏性阻尼系数分别为ρi0Gi0ηi0.第i层段桩周土体厚度为bi的内部扰动区域划分为n个径向圈层,其中第j圈层桩周土的密度、剪切模量和黏性阻尼系数分别为ρijGijη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)

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

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

Fig. 1
图 1 桩土动力相互作用模型 Fig. 1 Dynamic interaction model for pile-soil system


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)

式中wi0wij1分别为桩芯和桩周土体位移.

设第i层段桩身位移为wip,则桩身纵向振动方程为

$\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)

式中fiS0fiS1分别为第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)

式中Wij1(r, s)为wij1(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)

式中:I0(qij1r)、K0(qij1r)分别为零阶第一类、第二类修正Bessel函数,Aij1Bij1为待定系数.

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

$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)

式中:FFij1FFi(n+1)1分别为第i层段j圈层中外边界(r=ri(j+1))和内边界(r=rij)处的剪切刚度.

$\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)

式中Wi0(r, s)为wi0(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)递推所得FFi11和式(20)计算所得FFi0代入可得

$\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} $代表桩顶动刚度,Ki代表桩顶动阻尼.

3 算例分析参照文献[6]取径向圈层数n=20,纵向层段取m=5,如无特殊说明,管桩具体参数取值如下:ri0=0.38 m,ri1=0.50 m,ρiP=2 500 kg/m3EiP=25 GPa,H=6 m,kp=1 000 kN/m3δp=100 kN/m3.桩周各圈层土体相关参数按如下方式取值:$ {\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)}}} $为施工扰动系数,具体取Vi(n+1)=50 m/s,ηi(n+1)=10 kNs/m2qi=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可以发现,本文所得解析解退化后与已有解答吻合良好.

Fig. 2
图 2 本文解析解退化与已有解[18]对比验证 Fig. 2 Comparison of the proposed analytical solution with the solution in Ref.[18]


Fig. 3
图 3 本文解析解退化与已有解[6]对比验证 Fig. 3 Comparison of the proposed analytical solution with Wang's solution in Ref.[6]


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

Fig. 4
图 4 纵向软、硬夹层对桩顶复刚度曲线的影响 Fig. 4 Effect of vertical soft and hard interlayer on complex stiffness curves of pile head


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

表 1
表 1 桩周各层段施工软化程度和硬化程度工况 Tab. 1 Soil softening and hardening degree due to construction disturbance 工况 q1 q2 q3 q4 q5

S1 0.60 0.55 0.50 0.45 0.40

S2 0.80 0.75 0.70 0.65 0.60

S3 1.00 0.95 0.90 0.85 0.80

S4 1.00 1.00 1.00 1.00 1.00

H1 1.00 1.00 1.00 1.00 1.00

H2 1.20 1.15 1.10 1.05 1.00

H3 1.40 1.35 1.30 1.25 1.20

H4 1.60 1.55 1.50 1.45 1.40



表 1 桩周各层段施工软化程度和硬化程度工况 Tab. 1 Soil softening and hardening degree due to construction disturbance


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

Fig. 5
图 5 桩周土软化程度对桩顶复刚度曲线的影响 Fig. 5 Effect of softening degree of pile surrounding soil on complex stiffness curves of pile head


Fig. 6
图 6 桩周土硬化程度对桩顶复刚度曲线的影响 Fig. 6 Effect of hardening degree of pile surrounding soil on complex stiffness curves of pile head


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

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

Fig. 7
图 7 桩周土体软化范围对桩顶复刚度曲线的影响(软化工况S1) Fig. 7 Effect of softening range of pile surrounding soil on complex stiffness curves of pile head (softening condition S1)


Fig. 8
图 8 桩周土体硬化范围对桩顶复刚度曲线的影响(硬化工况H4) Fig. 8 Effect of hardening range of pile surrounding soil on complex stiffness curves of pile head (hardening condition H4)


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

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

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


参考文献
[1] 崔春义, 孟坤, 武亚军, 等. 考虑竖向波动效应的径向非均质黏性阻尼土中管桩纵向振动响应研究[J]. 岩土工程学报, 2018, 40(8): 1433.
CUI Chunyi, MENG Kun, WU Yajun, et al. Dynamic response of vertical vibration of pipe piles in soils with radial inhomogeneousity and viscous damping considering vertical wave effect[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(8): 1433. DOI:10.11779/CJGE201808008


[2] CUI Chunyi, MENG Kun, WU Yajun, et al. Dynamic response of pipe pile embedded in layered visco-elastic media with radial inhomogeneity under vertical excitation[J]. Geomechanics and Engineering, 2018, 16(6): 609. DOI:10.12989/gae.2018.16.6.609


[3] WU Wenbing, JIANG Guosheng, DOU Bin, et al. Vertical dynamic impedance of tapered pile considering compacting effect[J]. Mathematical Problems in Engineering, 2013(2): 1. DOI:10.1155/2013/304856


[4] NOVAK M, HAN Y C. Impedances of soil layer with boundary zone[J]. Journal of Geotechnical Engineering, 1990, 116(6): 1008. DOI:10.1061/(asce)0733-9410(1990)116:6(1008)


[5] EI NAGGAR M H. Vertical and torsional soil reactions for radially inhomogeneous soil layer[J]. Structural Engineering & Mechanics, 2000, 10(4): 299. DOI:10.12989/sem.2000.10.4.299


[6] 王奎华, 杨冬英, 张智卿. 基于复刚度传递多圈层平面应变模型的桩动力响应研究[J]. 岩石力学与工程学报, 2008, 27(4): 825.
WANG Kuihua, YANG Dongying, ZHANG Zhiqing. Study on dynamic response of pile based on complex stiffness transfer model of radial multizone plane strain[J]. Chinese Journal of Rock Mechanics & Engineering, 2008, 27(4): 825. DOI:10.3321/j.issn:1000-6915.2008.04.024


[7] 杨冬英, 王奎华. 任意圈层径向非均质土中桩的纵向振动特性[J]. 力学学报, 2009, 41(2): 243.
YANG Dongying, WANG Kuihua. Vertical vibration of pile in radially inhomogeneous soil layers[J]. Chinese Journal of Theoretical & Applied Mechanics, 2009, 41(2): 243. DOI:10.6052/0459-1879-2009-2-2007-186


[8] WANG Kuihua, YANG Dongying, ZHANG Zhiqing, et al. A new approach for vertical impedance in radially inhomogeneous soil layer[J]. International Journal for Numerical & Analytical Methods in Geomechanics, 2012, 36(6): 697. DOI:10.1002/nag.1024


[9] WANG Kuihua, YANG Dongying, ZHANG Zhiqing, et al. Vertical dynamic response of pile in a radially heterogeneous soil layer[J]. International Journal for Numerical & Analytical Methods in Geomechanics, 2009, 33(8): 1039. DOI:10.1002/nag.755


[10] 刘汉龙, 丁选明. 现浇薄壁管桩在低应变瞬态集中荷载作用下的动力响应解析解[J]. 岩土工程学报, 2007, 29(11): 1611.
LIU Hanlong, DING Xuanming. Analytical solution of dynamic response of cast-in-situ concrete thin-wall pipe piles under transient concentrated load with low strain[J]. Chinese Journal of Geotechnical Engineering, 2007, 29(11): 1611. DOI:10.3321/j.issn:1000-4548.2007.11.003


[11] 丁选明, 刘汉龙. 低应变下变阻抗薄壁管桩动力响应频域解析解[J]. 岩土力学, 2009, 30(6): 1793.
DING Xuanming, LIU Hanlong. Analytical solution in frequency domain of dynamic response of thin-wall pipe piles with variable wave impedance under low strain transient concentrated load[J]. Rock & Soil Mechanics, 2009, 30(6): 1793. DOI:10.3969/j.issn.1000-7598.2009.06.046


[12] 丁选明, 陈磊. 粘弹性地基中管桩的纵向振动特性研究[J]. 长江科学院院报, 2009, 26(3): 32.
DING Xuanming, CHEN Lei. Study on vertical vibration characteristics of pipe piles in viscoelastic soil[J]. Journal of Yangtze River Scientific Research Institute, 2009, 26(3): 32. DOI:10.3969/j.issn.1001-5485.2009.03.009


[13] 郑长杰, 丁选明, 黄旭, 等. 滞回阻尼土中大直径管桩纵向振动响应解析解[J]. 岩石力学与工程学报, 2014(增刊1): 3284.
ZHENG Changjie, DING Xuanming, HUANG Xu, et al. Analytical solution of vertical vibration response of large diameter pipe pile in hysteretic damping soil[J]. Chinese Journal of Rock Mechanics & Engineering, 2014(S1): 3284.


[14] ZHENG Changjie, DING Xuanming, SUN Yifei. Vertical vibration of a pipe pile in viscoelastic soil considering the three-dimensional wave effect of soil[J]. International Journal of Geomechanics, 2015, 16(1): 1. DOI:10.1061/(asce)gm.1943-5622.0000529


[15] LI Zhenya, WANG Kuihua, WU Wenbing, et al. Vertical vibration of a large-diameter pipe pile considering the radial inhomogeneity of soil caused by the construction disturbance effect[J]. Computers & Geotechnics, 2017, 85: 90. DOI:10.1016/j.compgeo.2016.12.016


[16] 胡昌斌, 王奎华, 谢康和. 桩与黏性阻尼土耦合纵向振动时桩顶时域响应研究[J]. 振动工程学报, 2004, 17(1): 72.
HU Changbin, WANG Kuihua, XIE Kanghe. Time domain axial response of dynamically loaded pile in viscous damping soil layer[J]. Journal of Vibration Engineering, 2004, 17(1): 72. DOI:10.3969/j.issn.1004-4523.2004.01.016


[17] 廖振鹏. 工程波动理论导论[M]. 北京: 科学出版社, 2002.
LIAO Zhenpeng. Introduction to engineering wave theory[M]. Beijing: Science Press, 2002.


[18] 崔春义, 孟坤, 武亚军, 等. 径向非均质黏性阻尼土中管桩纵向振动特性[J]. 哈尔滨工业大学学报, 2019, 51(2): 136.
CUI Chunyi, MENG Kun, WU Yajun, et al. Dynamic impedance of vertical vibration of pipe pile embedded in radially inhomogeneous viscoelastic soil[J]. Journal of Harbin Institute of Technology, 2019, 51(2): 136. DOI:10.11918/j.issn.03676234.201801140



相关话题/动力 计算

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 工业纯钛金属织构标准极图的计算及分析
    工业纯钛金属织构标准极图的计算及分析陈亮维,刘状,虞澜,胡劲,易健宏(昆明理工大学材料科学与工程学院,昆明650093)摘要:工业纯钛中的金属织构会引起各向异性,获得织构信息及分析其演变规律对钛材加工与应用非常重要.本文利用单晶钛的晶体结构数据、乌氏网、极图与织构的定义,建立了纯钛的织构与特定晶面极 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 抗癌药物作用预测计算方法的研究现状与展望
    抗癌药物作用预测计算方法的研究现状与展望顾兆伟1,张立忠2,刘晓峰3,谭先4(1.长春中医药大学附属第三临床医院脑病康复科,长春130000;2.长春市朝阳区清和社区卫生服务中心,长春130000;3.空军杭州特勤疗养中心康复理疗科,杭州310000;4.东北师范大学信息科学与技术学院,长春1300 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 变刚度轴承-碰摩转子机动飞行动力学响应
    变刚度轴承-碰摩转子机动飞行动力学响应李小彭,陈仁桢,尚东阳,陈延炜(东北大学机械工程与自动化学院,沈阳110819)摘要:为研究机动飞行状态下高速滚动轴承-转子耦合系统非线性动力学特性,利用有限元方法建立爬升-俯冲机动飞行状态下含滚动轴承时变刚度的轴承-碰摩转子系统数学模型,并利用Newmark- ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 深沟球轴承复合故障动力学特征
    深沟球轴承复合故障动力学特征王凯1,剡昌锋1,王风涛2,常斌全3,吴黎晓1(1.兰州理工大学机电工程学院,兰州730050;2.兰州理工大学能源与动力工程学院,兰州730050;3.舜宇光学科技有限公司,浙江宁波315400)摘要:为提高状态监测水平,诊断和识别轴承的早期故障,降低因轴承故障导致的损 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 甲烷-正庚烷双燃料二阶段着火的动力学特性
    甲烷-正庚烷双燃料二阶段着火的动力学特性李锋,张尊华,梁俊杰,万琦,李格升(武汉理工大学能源与动力工程学院,武汉430063)摘要:基于CHEMKINPRO软件,选取NUI2016机理对甲烷-正庚烷双燃料进行研究.探究了甲烷-正庚烷双燃料在中低温条件下的二阶段着火特性,开展了相应的化学动力学分析.结 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 一种基于拉-剪破坏的边坡动力稳定性分析方法
    一种基于拉-剪破坏的边坡动力稳定性分析方法彭雪峰1,朱亚林1,2,马驰1,檀昆1(1.合肥工业大学土木与水利工程学院,合肥230009;2.土木工程结构与材料安徽省重点实验室,合肥230009)摘要:为量化地震作用下边坡的稳定性,优化矢量和方法并计算出相应动态安全系数,基于边坡在地震作用下的破坏模式 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 桥梁颤振临界风速的概率密度演化计算
    桥梁颤振临界风速的概率密度演化计算姜保宋,周志勇,唐峰(土木工程防灾国家重点实验室(同济大学),上海200092)摘要:针对桥梁结构自身特性以及外部环境的随机性(如刚度、质量、阻尼比、气动导数等因素)所造成的桥梁的颤振临界风速不确定,难以衡量桥梁颤振稳定性问题.将概率密度演化方法与桥梁颤振多模态耦合 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 桥梁伸缩缝跳车冲击荷载计算方法与模型实验
    桥梁伸缩缝跳车冲击荷载计算方法与模型实验丁勇1,2,王佩1,游玖昂1,诸葛萍1(1.宁波大学土木工程系,浙江宁波315211;2.桥梁工程结构动力学国家重点实验室(重庆交通科研设计院),重庆400067)摘要:为实测移动车辆对桥梁伸缩缝的冲击荷载,防止桥梁伸缩缝在这种冲击荷载作用下发生早期损坏,制作 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 小跨高比RC开口梁弹性扭转性能计算方法
    小跨高比RC开口梁弹性扭转性能计算方法陈圣刚1,谢群1,郭全全2,刁波2,叶英华2(1.济南大学土木建筑学院,济南250022;2.北京航空航天大学交通科学与工程学院,北京100191)[HJ1.5mm]摘要:为改善经典Vlasov弹性扭转理论应用于开口深梁计算时引起的极大计算误差,在Vlasov理 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 悬索桥桁架加劲梁动力等效成等截面欧拉梁方法
    悬索桥桁架加劲梁动力等效成等截面欧拉梁方法祝卫亮,葛耀君(土木工程防灾国家重点实验室(同济大学),上海200092)摘要:针对悬索桥桁架加劲梁在数值计算、模型设计中简化问题,基于单独桁架加劲梁动力特性分析结果,以竖弯、侧弯、扭转基频相等为原则,推导并给出了不同边界条件(悬臂和两端固结)对应等截面欧拉 ...
    本站小编 哈尔滨工业大学 2020-12-05