成层广义Gibson地基中桩的水平动力响应
上海交通大学 船舶海洋与建筑工程学院,上海 200240
Horizontal Dynamic Response of Pile Based on Layered Generalized Gibson Foundation
School of Navy Architecture, Ocean and Civil Engineering, Shanghai Jiao Tong University, Shanghai 200240, China
Received: 2020-11-19
作者简介 About authors
邱杰凯(1996-),男,四川省南充市人,硕士生,主要从事土与结构的动力相互作用.
为直接求解桩周土切变模量随深度线性变化情形下的单桩水平动力响应问题,基于成层广义Gibson地基,引入Adomian分解法,提出了一种求解非均质地基中桩水平动力响应的方法.相较于初参数法和传递矩阵法,该方法无需对非均质地基进行离散化处理.与数值方法相比,具有计算成本低、精度高和收敛速度快等优点.通过将本文方法和分层方法、均质地基解析方法、数值方法的计算结果进行对比,验证了本文方法的正确性和合理性.分析了桩底边界条件、地基土参数及桩身长细比对桩水平动力响应的影响.结果表明:在非均质地基中,桩土刚度比是影响桩水平动力响应的重要因素,随着地基土切变模量的增大,桩身水平位移幅值减小且分布趋于平缓.此外,相比于其他参数,土层泊松比和阻尼比对桩的水平动力响应影响较小.
关键词:
In order to directly solve the horizontal dynamic response of a single pile under the condition that the shear modulus of the soil around the pile varies linearly with the depth, a method for solving the horizontal dynamic response of piles in nonhomogenous foundation is proposed based on the layered generalized Gibson foundation and the Adomian decomposition method. Compared with the initial parameter method and the transfer matrix method, the proposed method does not need to discretize the nonhomogeneous foundation. Compared with the numerical method, it has the advantages of low computational cost, high accuracy, and fast convergence speed. The correctness and rationality of this method are verified by comparing it with the calculation results of the layered method, the analytical method of homogeneous foundation, and the numerical method. The influence of boundary conditions on the pile bottom, soil parameters, and pile slenderness on the horizontal dynamic response of the pile is investigated. The result show that in nonhomogeneous foundation, the pile-soil elastic modulus ratio is an important factor affecting the horizontal dynamic response of pile. As the shear modulus of soil increases, the amplitude of horizontal displacement of the pile decreases, and the distribution tends to be gentle. In addition, compared with other parameters, Poisson’s ratio and damping ratio of soil have less influence on the horizontal dynamic response of pile.
Keywords:
本文引用格式
邱杰凯, 丁肇伟, 宋春雨, 陈龙珠.
QIU Jiekai, DING Zhaowei, SONG Chunyu, CHEN Longzhu.
在桩的水平动力学理论研究中,关于桩周土非均质性的研究主要集中在土层的竖向自然分层非均质性[5,6,7,8],对由于上覆土压力、土体自重等原因引起自然土层中土体切变模量随深度连续变化的非均质性研究较少.对于自然土层内的非均质性,一种思路是用多层均质地基近似非均质地基,多见于荷载作用下地基的变形和边坡稳定问题研究[9,10,11].值得注意的是,采用有限元模拟非均质地基时,一般需要将不同深度单元的剪切刚度设置为与深度相关,这本质上也是一种分层[12,13].用分层方法考虑地基土的竖向非均质性时,划分层数越多,结果越接近实际情况.文献[14]在群桩的动阻抗研究中基于初参数法和传递矩阵法,采用15层均质地基近似单层广义Gibson地基,得到了较好的近似结果.但传递矩阵等解析方法需进行大量的矩阵拼装和迭代,层数过多会给计算带来极大困难,难以得到问题的封闭解,即使采用有限元方法,在单元网格划分过细时往往也存在内存不足、计算耗时严重和修正误差过大的问题.因此,在非均质地基相关问题的研究中,相较于多层均质地基,采用连续地基形式显得更为简单合理.
Gibson首先应用解析法考虑了一种土体切变模量
上述对各类Gibson地基中桩的特性研究都集中于静力学问题,而未涉及对动力学问题的研究.本文针对成层广义Gibson地基,基于Novak薄层单元法和Adomian微分分解法(ADM),得到了非均质土体内桩在各种边界条件下的水平动力响应近似解.Adomian方法能很好地解决特征方程法难以求解高阶变系数微分方程的问题,无需对非均质地基进行离散化,从而不用修正离散误差,相比于数值解法,具有计算成本低、精度高和收敛速度快等优点[19,20],并且可通过符号运算来实现[21].通过将本文方法与分层方法、均质地基解析方法、数值方法的计算结果进行对比,验证了本文方法的正确性和合理性.最后,结合具体算例分析了桩底边界条件、地基土参数及桩身长细比对桩水平动力响应的影响.
1 计算模型
桩土计算模型如图1所示,构建成层广义Gibson地基中单桩水平振动模型,设桩土连续紧密接触,桩为线弹性材料的圆截面桩,考虑桩周土由于历史沉积等原因引起的自然分层,设土层沿埋深方向自然分为
图1
图1
成层广义Gibson地基中单桩水平振动模型
Fig.1
Model of lateral vibration of a single pile in layered generalized Gibson foundation
式中:Ep、Ip、Ap、ρp分别为桩身弹性模量、截面惯性矩、截面积及密度;
2 基于成层Gibson地基和ADM的计算方法
假设桩身运动为简谐稳态振动,桩身位移可表达为
将式(5)代入式(1),有:
沿埋深各层有:
参考文献[23]中的Adomian分解法,桩水平动力位移复振幅可以分解为无穷级数形式:
式中:
将式(10)代入式(9),有:
式(11)的近似解可用递推形式写出:
式中:
对如图1所示的桩土相互作用模型,设桩顶边界条件为约束转角,桩底边界条件分别考虑自由、铰支和固定3种情况,结合桩身在土层交界面的连续性条件,具体边界条件如下:
(1) 桩顶边界条件.
(2) 桩底边界条件.
(a) 桩底自由:
(b) 桩底铰支:
(c) 桩底固定:
(3) 桩身连续条件.
将式(15)代入式(12)、(13),有:
式中:M0为桩顶弯矩.
将式(19)~(22)代入式(12)、(13),有:
整理递推式(23)~ (26),代入式(16)~(18)中可得:
式中:
3 算例对比与参数分析
通过理论推导和数值计算,得到成层广义Gibson地基中桩的水平动力位移幅值
由于篇幅所限,只针对部分边界条件做详细讨论,后续分析如不作特殊说明,边界条件均为桩头约束转角、桩底自由的情况.首先,基于单层广义Gibson地基,分析了分解级数N对结果的影响,验证了本文方法的收敛性.通过和已有文献的对比,验证本文解的合理性和正确性.最后,基于双层广义Gibson地基,分析了桩底边界条件和桩土参数对桩水平动力响应的影响.
3.1 收敛性验证
基于单层广义Gibson地基验证分解级数N对桩顶和桩底位移幅值影响,桩、土参数取值分别取值如下:dp=0.8 m, l/dp=10, Ep/Es1(0)=5000, ρp/ρs1=1.25, βs1=5%, vs1=0.4, a01=0.5,Es1是第一层土体表面弹性模量.土层切变模量变化率 k1=Gs1(0),计算结果如表1所示.随着分解级数N的增加,桩顶、桩底水平振动位移响应幅值逐渐收敛,不妨取一个极小精度值10-7,当N≥7时,有:
式中:
表1
桩端水平动力位移幅值与
Tab.1
N | ||
---|---|---|
2 | 0.819792340 | 2.708522741 |
3 | 0.575465393 | 0.438747062 |
4 | 0.588914776 | 0.094956192 |
5 | 0.589908739 | 0.103157714 |
6 | 0.589898262 | 0.103125582 |
7 | 0.589898281 | 0.103125562 |
8 | 0.589898281 | 0.103125562 |
9 | 0.589898281 | 0.103125562 |
10 | 0.589898281 | 0.103125562 |
11 | 0.589898281 | 0.103125562 |
12 | 0.589898281 | 0.103125562 |
3.2 和已有文献对比
本节通过和已有文献对比,分别将本文方法的计算结果与分层方法[14,25]、均质地基解析方法[24]、数值方法[12]的计算结果进行对比,从而验证本文解的正确性.首先,对单层广义Gibson地基中的桩-土水平动力问题用本文方法和分层方法进行计算比较,分层方法参考了文献[14]将广义Gibson地基近似为多层均质地基的思路,将单个自然非均质土层划分为ζ层均质土层,分层后的土层在各交界面处的位移、转角、弯矩和剪力均满足连续性条件,每个均质土层的剪切刚度取为层内均值,然后利用文献[25]的初参数法和传递矩阵法,得到单层广义Gibson地基内桩的水平动力位移响应.边界条件为桩头、桩底均自由,桩、土参数为:dp=0.5 m, l/dp=10, Ep/Es1(0)=800, ρs1/ρp=0.7, βs1=5%, vs1=0.35, a01=0.5, 土体模量随深度线性增加,底部模量是顶部的两倍.对比结果如图2所示,可以看到,随着划分层数的不断增加,采用分层方法的结果越来越接近实际的非均质情况,当划分层数足够大,可用多层均质地基近似代替非均质地基,但其计算量极大,耗时甚长,对于多层非均质地基处理起来较为困难,存在一定的局限性.
图2
图2
本文方法与分层方法的结果的对比
Fig.2
Comparisons of results of the method proposed and those of layered method
图3
图3
本文退化解结果与均质地基中结果的对比
Fig.3
Comparisons of results of the degenerate solution proposed and those of homogeneous foundation
图4
图4
本文结果与数值计算结果的对比
Fig.4
Comparisons of results of this paper and results of numerical calculation
3.3 参数分析
采用双层广义Gibson地基研究桩底边界条件、桩身参数
图5
图5
不同λ下双层地基切变模量的变化
Fig.5
Variation of shear modulus of double-layers subgrade at different λ values
为研究在不同非均质条件下不同桩底条件对桩身水平动力位移幅值的影响,给出3种桩底边界条件下的桩水平动力位移幅值曲线,结果如图6所示.桩头均为约束转角边界条件,桩身直径dp=0.8 m,上下土层切变模量比λ分别为0.5、2,m取为0.5、1.5.通常情况下,桩土刚度比都在 1000 甚至 10000 以上,为使结果具有代表性,桩土刚度比取为 5000.此外,如不作特殊说明,后续分析中采用的桩、土参数均为:l/dp=10, l1=l2, ρp/ρs1=ρp/ρs2=1.25, βs1=βs2=5%, vs1=vs2=0.4, a01=0.5.
图6
图6
桩底边界条件对桩身水平位移幅值的影响
Fig.6
Influence of tip conditions on amplitude of horizontal displacement of pile
可见,桩底边界条件对桩身位移曲线的影响主要体现在桩底部的水平动力响应,3种边界条件对应的位移曲线差别明显.相比于下部土层较软(λ>1)的情况,下部土层较硬(λ<1)时位移曲线差别较小.图6反应了桩周土刚度变化对桩水平动力响应的影响.桩周土刚度越大,桩身水平动力位移幅值越小,不同边界条件对应的位移幅值曲线差别也愈小.
为研究上下土层土刚度比在不同非均质条件下对桩的影响,令m取为0.5、1.5,分两种情况进行讨论:① 保持桩土刚度比Ep/Es2(0)=5000 不变,Es1(0)/Es2(0) 分别取为0.5、1、5,探讨上部土层性质对桩水平动力响应的影响,结果如图7所示.②保持桩土刚度比Ep/Es1(0)=5000 不变,Es2(0)/Es1(0) 分别取为0.5、1、5,探讨下部土层性质对桩水平动力响应的影响,结果如图8所示.从图7、8中可以看出,表层土对桩的水平动力响应影响极为显著,随着上部土层刚度的增大,桩的水平动力位移幅值迅速减小,下部土层刚度变化对桩的水平动力响应影响和上部土层类似,但远不如上层土影响显著,图8(a)、8(b)表明,随着桩周土刚度增大,下部土层刚度对桩水平动力响应的影响愈小,表明Ep/Es1(0)是影响桩水平动力响应的重要因素,下层土的桩土刚度比对桩的水平动力响应影响远不如上层土.可见在桩的水平动力响应中表层土起重要作用,桩的水平动力承载能力主要取决于上部土层性质.
图7
图7
桩土刚度比Ep/Es1(0)对桩身水平位移幅值的影响
Fig.7
Influence of Ep/Es1(0) on amplitude of horizontal displacement of pile
图8
图8
桩土刚度比Ep/Es2(0)对桩身水平位移幅值的影响
Fig.8
Influence of Ep/Es2(0) on amplitude of horizontal displacement of pile
为研究不同荷载频率的影响,取无量纲荷载频率a01为0.1、1.0,上、下土层切变模量比λ分别为0.5、2,其桩身水平动力位移幅值如图9所示.从图中可知,在土层软硬不同的情况下,保持其他条件不变,桩身水平位移响应均随无量纲频率a01增大而减小,这是因为桩的惯性力与荷载作用方向相反,频率越高,惯性力越大,抵消的载荷越大,进而导致桩身水平位移响应减小.因此,低频简谐振动对桩的水平位移响应影响更显著,在设计水平动力受荷桩时,应尤其注意低频荷载的影响.
图9
图9
λ不同时无量纲荷载频率对桩身水平位移幅值的影响
Fig.9
Influence of dimensionless load frequency on amplitude of horizontal displacement at different values of λ
为研究桩身长细比的影响,取长细比分别为5、15、40,上下层切变模量比λ分别为0.5、2,结果如图10所示.从图中可知,中短桩的桩身长细比对桩的水平动力位移幅值影响较大,当长细比较小时(l/dp=5),桩的相对刚度较大,可按刚性桩考虑, 桩身位移以平动和转动为主;当桩周土整体剪切刚度较大且下部土层较硬(λ<1)时,桩更容易发生转动位移,这在实际工程中需要引起注意,尤其当土体非均质性较为明显时;当下部土层较软(λ>1)时,随着长细比增加,桩周土土体复动阻抗增大,桩身水平位移幅值逐渐减小.当长细比增加到一定程度,桩的水平动力响应随深度减小,深度超过10倍桩径以后桩的水平动力响应很小且几乎无变化,下部土层表现出一定的嵌固效应.这与均质地基求解结果不同,这是因为均质地基未考虑桩周土刚度随深度增加,在桩的长细比较大时,底部依旧有较明显的动力响应[24],这与实际是不相符的.随着长细比的进一步增大,由图10(e)、10(f)可知当桩长超过有效桩长时,桩深部的水平动力响应变得非常小,几乎可忽略不计,这与胡安峰等[24]的临界长径比结论是一致的.图10(e)、10(f)中λ不同时,其桩身水平动力响应几乎一致,这是因为计算时图10(e)、10(f)中上部土层刚度是一致的,变化的仅是下部土层刚度,可以看到长桩的有效桩长小于上部土层厚度时,底部土层的软硬对长桩动力响应影响极小,这与前面探讨上下土层刚度影响得到的结论是一致的.因此,考虑长桩的水平动力承载能力时应更多关注有效桩长范围内的土层性质.
图10
图10
λ不同时桩长细比对桩身水平位移幅值的影响
Fig.10
Influence of pile slenderness on amplitude of horizontal displacement at different values of λ
为分析土体泊松比的影响,分别取上、下土层的泊松比νs1、νs2为0.1、0.3,上、下层切变模量比λ分别为0.5、2,分析结果如图11所示.可以看出:随着土层泊松比增加,桩身水平动力位移幅值逐渐增大,这是因为泊松比跟桩周土刚度存在一定关系,泊松比增加,桩周土刚度减小,导致桩身水平动力位移幅值增加,这与实际工程是相符合的.值得指出的是,如果不考虑泊松比跟桩周土刚度的关系,直接定义桩周土剪切刚度取值而非采用地基土弹性模量换算,在理论分析中会出现随着泊松比增加,桩的水平动力位移幅值减小的情况.此外,可以看出土层泊松比的变化对桩上部的水平响应影响较大,土层泊松比对桩的位移中性点几乎无影响.
图11
图11
λ不同时土体泊松比对桩身水平位移幅值的影响
Fig.11
Influence of Poisson's ratio of soil on amplitude of horizontal displacement of pile at different values of λ
为分析土体阻尼比的影响,分别取上下土层阻尼比为0.01、0.2,分析结果如图12所示.m较小时,土体阻尼比的增加会导致桩身水平动力位移幅值减小,这是因为阻尼比增加,导致土体复动阻抗增大,桩身水平动力响应减小.与泊松比变化类似,土层阻尼比变化对桩的位移中性点几乎无影响, 阻尼比变化对桩顶部的水平动力位移幅值影响较大,因此,在实际工程中考虑桩的水平振动,需特别关注桩顶部的动力响应.综合来看,相比其他因素,土体泊松比和阻尼比对桩的水平动力响应影响较小.
图12
图12
λ不同时土体阻尼比对桩身水平位移幅值的影响
Fig.12
Influence of damping ratio of soil on amplitude of horizontal displacement of pile at different values of λ
综合分析图4~10,可以看出切变模量随深度连续变化的非均质性对桩的水平动力位移幅值影响很大,即m越大,桩身水平动力位移幅值越小且分布趋于平缓,由物理意义可知,当地基模型退化为均质地基时(m=0),水平位移幅值达到最大.随着m逐渐增大,桩身位移中性点逐渐上移,上部土层承担更多荷载.故可以考虑在实际工程中可通过地基加固的方式来控制桩的水平振动幅值,避免桩产生过大位移.
分析上下土层软硬对桩的水平动力位移幅值影响,保持上部土层刚度不变,当下部土层较硬(λ<1)且m较小时,桩身水平动力位移幅值较下部土层较软(λ>1)时整体偏小,桩身位移中性点下移,下部土层承受更多的动力荷载并表现出一定的嵌固效应.需要指出的是,在上部土层剪切刚度相同情况下,下部土层剪切刚度变化对桩顶部水平动力响应影响程度跟地基土的非均质性有关,m较大时影响很小,当m较小时这种影响较为明显,这与工程实际是相符合的.
4 结论
本文基于Novak薄层单元法和Adomian分解方法,提出了一种计算竖向非均质地基中单桩水平振动的方法.与分层法相比,本方法对切变模量随深度连续变化的地基无需简化为多层均质地基求解,计算过程简单,收敛速度快,适用范围广,避免了数值方法的繁冗计算量,具有良好的收敛性和适用性.利用该计算方法详细分析了土体非均质性、上下土层软硬、桩土参数等对桩的水平振动响应影响,得到结论如下:
(1) 地基土切变模量的变化对桩的水平动力位移幅值影响很大,随着表征地基土非均质性的m值的增大,桩的水平动力位移幅值减小且分布趋于平缓,桩中性点上移,上部土层承受更多动力荷载.
(2) 双层地基中,当上部土层刚度保持一定时,下部土层较硬(λ<1)情况下的水平动力位移幅值小于下部土层较软(λ>1)情况下的位移幅值,桩位移中性点下移,下部土层承受更多的动力荷载并表现出一定的嵌固效应.其他条件保持不变时,随着荷载频率的增加,桩的水平动力响应幅值逐渐减小.在长细比较小时(l/dp=5),桩的相对刚度较大,桩的位移以平动和转动为主.当下部土层较软时,随着长细比增加,桩周土土体复动阻抗增大,桩身水平位移幅值逐渐减小.当长细比增加到一定程度,桩的水平动力响应随深度减小,深度超过10倍桩径以后桩的水平动力响应很小且几乎无变化,下部土层表现出一定的嵌固效应.随着长细比进一步增大,当桩长超过有效桩长时,桩的水平动力响应变得非常小,几乎可忽略不计.当长桩的有效桩长小于上部土层厚度时,底部土层的软硬对长桩动力响应影响极小.因此,考虑长桩的水平动力承载能力时应更多关注有效桩长范围内的土层性质.
(3) 土性参数中,桩土相对刚度比是影响桩的水平动力响应的重要因素,下层土性质对桩水平动力响应的影响远小于表层土.因此,在设计水平动力受荷桩时,要尤其注意表层土的工程性质.
(4) 相比于其他参数,土层泊松比和阻尼比对桩的水平动力响应影响较小,且影响主要集中在桩顶部的动力响应.因此,在进行桩的水平动力分析时,需特别关注桩的顶部响应.
参考文献
Dynamic pile impedances for fixed-tip piles
[J].DOI:10.1016/j.soildyn.2017.03.025 URL [本文引用: 1]
Dynamic pile impedances for laterally-loaded piles using improved Tajimi and Winkler formulations
[J].DOI:10.1016/j.soildyn.2016.09.020 URL [本文引用: 1]
A closed-form solution for kinematic bending of end-bearing piles
[J].DOI:10.1016/j.soildyn.2017.09.004 URL [本文引用: 1]
剪切模量随深度呈指数函数增大地基中锚杆弹塑性分析
[J].
Elastic-plastic analysis for bolts in ground with shear modulus increasing following an exponential function with depth
[J].
Transient analysis of a fixed-head pile group in multi-layered transversely isotropic media due to horizontal loadings
[J].DOI:10.1016/j.compgeo.2020.103772 URL [本文引用: 1]
A coupling rational finite element-boundary element method for a laterally loaded pile in transversely isotropic poroelastic soils
[J].DOI:10.1016/j.compgeo.2019.103227 URL [本文引用: 1]
Time-dependent response of laterally loaded piles and pile groups embedded in transversely isotropic saturated viscoelastic soils
[J].DOI:10.1016/j.compgeo.2020.103815 URL [本文引用: 1]
层状地基中基于Laplace变换的桩基横向振动阻抗计算
[J].
Lateral vibration impedance of piles embedded in layered soil based on Laplace transform
[J].
高速铁路路基动力学研究进展
[J].
The review of high-speed railway track foundation dynamics
[J].
层状路基路面体系温度应力和变形问题的位移函数解法
[J].
Displacement function method on thermal stress and thermal deformation problems for layered pavement structure
[J].
Stability analysis of a layered slope with failure mechanism of a composite slip surface
[J].DOI:10.1061/(ASCE)GM.1943-5622.0001417 URL [本文引用: 1]
Effects of separation on the behavior of soil-pile interaction in liquefiable soils
[J].DOI:10.1061/(ASCE)GM.1943-5622.0000074 URL [本文引用: 5]
Kinematic bending of fixed-head piles in nonhomogeneous soil
[J].DOI:10.1061/(ASCE)GT.1943-5606.0001270 URL [本文引用: 1]
Horizontal impedance of pile groups considering shear behavior of multilayered soils
[J].DOI:10.1016/j.sandf.2014.09.001 URL [本文引用: 3]
Estimation of the dynamic surface modulus of a generalized Gibson soil from the rocking frequency of rectangular foundations
[J].DOI:10.1680/geot.1973.23.1.23 URL [本文引用: 1]
成层Gibson地基中单桩沉降的非线性分析
[J].
Nonlinear analysis of settlement of axially loaded single pile in Gibson soil
[J].
非均质地基中V-T联合受荷桩承载力分析
[J].
Bearing capacity analysis of a single pile under combined V-T Loads in non-homogeneous subsoil
[J].
Investigation of the torsional behaviour of circular piles in double-layered nonhomogeneous soil
[J].DOI:10.1016/j.apor.2020.102110 URL [本文引用: 1]
Solving nonlinear integral equations and nonlinear integro-differential equations using laplace adomian decomposition method through sagemath
[J].
Accelerating the convergence of Adomian decomposition method (ADM)
[J].DOI:10.1016/j.jocs.2018.12.014 URL [本文引用: 1]
An algorithm for Adomian decomposition method
[J].DOI:10.1016/j.amc.2003.10.037 URL [本文引用: 1]
Dynamic soil reactions for plane strain case
[J].DOI:10.1061/JMCEA3.0002392 URL [本文引用: 1]
基于Pasternak地基模型和Adomian分解方法的扩底桩水平动力响应分析
[J].
Horizontal dynamic response of a belled pile based on the Pasternak foundation model and the Adomian decomposition method
[J].
水平荷载下单桩动力反应分析
[J].
Dynamic response analysis for a single pile subjected to lateral load ing
[J].
饱和土-管桩-刚性承台水平振动研究
[J].
A study on the lateral vibration of saturated soil-pipe pile-rigid platform
[J].
Nonlinear dynamic behavior of pile foundations: Effects of separation at the soil-pile interface
[J].DOI:10.3208/sandf.46.437 URL [本文引用: 1]
/
〈 |
|
〉 |
