CN109344526A - 一种基于虚土桩模型的饱和层状土中桩基纵向振动研究方法 - Google Patents
一种基于虚土桩模型的饱和层状土中桩基纵向振动研究方法 Download PDFInfo
- Publication number
- CN109344526A CN109344526A CN201811217886.XA CN201811217886A CN109344526A CN 109344526 A CN109344526 A CN 109344526A CN 201811217886 A CN201811217886 A CN 201811217886A CN 109344526 A CN109344526 A CN 109344526A
- Authority
- CN
- China
- Prior art keywords
- pile
- soil
- saturated
- layer
- displacement
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000002689 soil Substances 0.000 title claims abstract description 237
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000011160 research Methods 0.000 title abstract description 5
- 229920006395 saturated elastomer Polymers 0.000 claims abstract description 190
- 238000006073 displacement reaction Methods 0.000 claims description 111
- 239000007787 solid Substances 0.000 claims description 81
- 230000005284 excitation Effects 0.000 claims description 30
- 239000012530 fluid Substances 0.000 claims description 30
- 230000008878 coupling Effects 0.000 claims description 27
- 238000010168 coupling process Methods 0.000 claims description 27
- 238000005859 coupling reaction Methods 0.000 claims description 27
- 238000005316 response function Methods 0.000 claims description 27
- 230000004044 response Effects 0.000 claims description 24
- 230000035699 permeability Effects 0.000 claims description 18
- 239000002245 particle Substances 0.000 claims description 17
- 230000009471 action Effects 0.000 claims description 12
- 238000013016 damping Methods 0.000 claims description 9
- 230000001133 acceleration Effects 0.000 claims description 6
- 210000001124 body fluid Anatomy 0.000 claims description 6
- 239000010839 body fluid Substances 0.000 claims description 6
- 230000005484 gravity Effects 0.000 claims description 6
- 239000011148 porous material Substances 0.000 claims description 6
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 239000004927 clay Substances 0.000 claims description 3
- 230000006835 compression Effects 0.000 claims description 3
- 238000007906 compression Methods 0.000 claims description 3
- 230000002950 deficient Effects 0.000 claims description 3
- 229920001971 elastomer Polymers 0.000 claims description 3
- 239000000806 elastomer Substances 0.000 claims description 3
- 230000003993 interaction Effects 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims description 3
- 238000010008 shearing Methods 0.000 claims description 3
- 238000001514 detection method Methods 0.000 abstract description 2
- 230000000694 effects Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000009933 burial Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Piles And Underground Anchors (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明提供一种基于虚土桩模型的饱和层状土中桩基纵向振动研究方法。本发明方法基于饱和虚土桩模型的浮承桩纵向振动动力阻抗算法系统，其采用的饱和虚土桩模型，能考虑桩周及桩底土体的饱和特性，且能合理考虑桩底土厚度等因素对桩基纵向振动特性的影响，适用于饱和土中浮承桩纵向振动特性研究，可为桩基动力检测提供理论指导和参考作用。
Description
技术领域
本发明涉及土建技术领域，具体而言，尤其涉及一种基于虚土桩模型的饱和层状土中桩基纵向振动研究方法。
背景技术
目前，针对浮承桩情况的问题，已有研究大多将桩底土简化为Winkler模型，其弹簧和阻尼器系数通常按经验取值，无法合理考虑桩底土体波动效应的影响。为解决此类问题，一些学者将桩底土考虑为单相或饱和弹性半空间介质，计算得出桩底复阻抗函数表达式，并对浮承桩纵向振动特性进行了分析。然而，弹性半空间模型虽可考虑桩底土波动效应，但其只适用于基岩埋深较大情况，且无法考虑桩底土厚度及成层特性对桩基纵向振动特性的影响。基于此点考虑，提出了桩与桩底土完全耦合单相介质虚土桩模型。而单相虚土桩模型均假定桩底土体为单相介质，未考虑桩底土饱和两相介质性，这对于饱和土中浮承桩基纵向振动问题并不适用。
发明内容
根据上述提出的技术问题，而提供一种基于虚土桩模型的饱和层状土中桩基纵向振动研究方法。
本发明采用的技术手段如下：
一种基于虚土桩模型的饱和层状土中桩基纵向振动研究方法，包括以下步骤：
S1、构建基于平面应变法及桩底饱和虚土桩模型的饱和土-桩-饱和虚土桩耦合体系纵向振动力学简化模型；
S2、将桩-土耦合振动系统按土体沿纵向分成m层，其中桩底土分成n层，假定桩周和桩底每层土都为均质、各向同性的饱和线黏弹性介质，根据Biot动力波动理论，建立第j层饱和土体和第k层饱和土体的平面应变条件下饱和土层动力控制方程及实体桩与饱和虚土桩界面处的边界条件和桩-土耦合条件；
S3、将第j层饱和土体和第k层饱和土体的谐和激振作用下饱和土层质点位移方程代入所述饱和土层动力学控制方程中，求解桩底第j层饱和土体土骨架径向位移和桩周第k层饱和土体土骨架径向位移，
基于求解的径向位移求解出桩底土与饱和虚土桩界面处、桩周土与实体桩界面处剪应力，通过饱和虚土桩边界条件和各层桩界面处位移连续、力的平衡条件及阻抗函数传递公式，对实体桩桩顶动力阻抗函数进行求解；
S4、由桩顶位移阻抗函数求得桩顶位移频率响应函数，求得桩顶速度频率响应函数、单位脉冲激励作用下桩顶速度时域响应，进而求得桩顶时域速度响应函数，
基于求得的所述桩顶速度频率响应函数和桩顶时域速度响应函数对桩身振动特性及桩身完整性进行评判。
进一步地，所述步骤S2假定条件还包括：
桩底各层土体为渗透性较差的饱和粘土；
桩周土和桩底土为一系列相互独立的薄层，不考虑土层间的相互作用；
实体桩为均质等截面弹性体，桩身相邻层之间满足力平衡和位移连续条件；
饱和虚土桩为等截面饱和两相介质，饱和虚土桩相邻层之间满足力平衡和位移连续条件，其与实体桩界面位移连续、应力平衡；
桩-土耦合振动系统满足线弹性和小变形条件，桩土界面完全接触，不存在滑移和脱离。
进一步地，所述步骤S2具体为：
S21、根据Biot动力波动理论，建立平面应变条件下饱和土层动力学控制方程为：
式中，各参数关系如下：
其中，将桩-土耦合振动系统按土体沿纵向分成m层，其中桩底土分成n层，由基岩往上依次编号为1、…、j、…、n、n+1、…、k、…、m层，各层厚度分别为l1、…、lj…、ln、ln+1、…、lk、…、lm，各层顶部埋深分别为h1、…、hj…、hn、hn+1、…、hk、…、hm，j＝1、2、…、n，k＝n+1、n+2、…、m，
式中，uj、wj为桩底第j层土骨架纵向位移和流体相对于土骨架纵向位移，Nj、ηj、bj、kj、Gj、ξj分别为桩底第j层饱和土体流体密度、土颗粒密度、孔隙率、流体粘滞系数、骨架与孔隙流体的粘性耦合系数、Biot定义的动力渗透系数、土体达西定律渗透系数、土体复剪切模量、土体剪切模量和阻尼比；uk、wk为桩周第k层土骨架纵向位移和流体相对于土骨架纵向位移，Nk、ηk、bk、kk、 Gk、ξk分别为桩周第k层饱和土体流体密度、土颗粒密度、孔隙率、流体粘滞系数、骨架与孔隙流体的粘性耦合系数、Biot定义的动力渗透系数、土体达西定律渗透系数、土体复值剪切模量、土体剪切模量和阻尼比；g为重力加速度，
S22、基于Biot动力波动理论得出渗透性较差时饱和土体一维纵向振动控制方程为：
其中，为第j层饱和虚土桩纵向位移，λj为第j层桩底土拉梅常数，且λj＝2νjGj/(1-2νj)，νj为桩底第j层饱和土体泊松比。αj、Mj为Biot定义的表征桩底第j层饱和土体土颗粒及流体压缩性常数， 和分别为桩底第j层饱和土体土颗粒、流体及土骨架的体积压缩模量；
S23、将桩底第j层土在饱和虚土桩界面处的剪切应力代入式(4)中，可得桩底第j层饱和虚土桩纵向振动控制方程为：
式中，为剪应力幅值，ω为激振圆频率，r0表示实体桩半径，Ap表示实体桩的横截面面积；
取第k段实体桩桩身微元体作动力平衡分析可得第k段实体桩纵向振动控制方程为：
式中，为第k段实体桩质点纵向振动位移，EP、ρP分别为实体桩弹性模量和密度，为第k层桩周饱和土体与实体桩界面处土体剪应力，为剪应力幅值。
进一步地，所述步骤S2中，饱和土-桩-饱和虚土桩体系边界条件如下：
饱和土体：
径向无穷远处位移为零，即
uj(∞,t)＝0 (7a)
uk(∞,t)＝0 (7b)
实体桩：
桩顶平衡条件：
各层界面位移连续条件：
各层界面力平衡条件：
饱和虚土桩：
桩底位移：
各层界面位移连续条件：各层界面力平衡条件：
桩与饱和虚土桩界面处边界条件：
位移连续条件：
力平衡条件：
桩-土耦合条件：
其中，HP为实体桩桩长，表示桩顶作用谐和激振力，表示激振力幅值，H表示基岩上土层总厚度，Hsp表示桩底土层厚，即饱和虚土桩长度。
进一步地，所述步骤S3中，谐和激振作用下饱和土层质点位移方程具体为：
式中：为桩底第j层饱和土体土骨架径向位移和流体相对于土骨架径向位移响应幅值，为桩周第k层饱和土体土骨架径向位移和流体相对于土骨架径向位移响应幅值，i为虚数单位，
所述步骤S3中，求解纵向位移解具体为：
S31、将谐和激振作用下饱和土层质点位移方程(12)代入所述饱和土层动力学控制方程(1)、(2)中，可得：
将式(13b)、(14b)分别代入式(13a)、(14a)可得：
式中：
方程(15)的通解为：
式中，Aj、Bj、Ak、Bk为待定常数，I0(qjr)、I0(qkr)、K0(qjr)、K0(qkr)为零阶第一类、第二类虚宗量贝塞尔函数。
由边界条件式(7)可知Bj＝Bk＝0，则进一步可得：
所述步骤S3中，对实体桩桩顶动力阻抗函数进行求解具体为：
S32、谐和激振作用下饱和虚土桩和实体桩质点纵向振动位移满足下式：
式中，分别为第j层饱和虚土桩和第k层实体桩质点纵向振动位移响应幅值，则式(5)和式(6)可化简为：
由边界条件式(11)可得：
将式(20)代入式(17)，并由剪应力与位移的关系可得出桩底土与饱和虚土桩界面处、桩周土与实体桩界面处剪应力分别为：
将式(21)代入式(19)，并进一步化简可得：
式中，
则方程(22)的通解为：
式中：Cj、Dj、Ck、Dk为待定常数；
S33、由饱和虚土桩边界条件式(9a)可得：
根据位移阻抗函数的定义可得第1段饱和虚土桩与第2段饱和虚土桩界面处的位移阻抗函数为：
综合第1段饱和虚土桩与第2段饱和虚土桩界面处位移连续、力平衡条件可得：
由此可得第2段饱和虚土桩与第3段饱和虚土桩界面处的位移阻抗函数为：
综合式(24)—(27)可得饱和虚土桩阻抗函数传递公式：
式中：
由式(28)递推可得饱和虚土桩与实体桩界面处阻抗函数：
式中：
S34、根据实体桩与饱和虚土桩界面耦合条件式(10)可得第n+1段实体桩与第n+2段实体桩界面处阻抗函数为：
式中：
由此可得实体桩阻抗函数传递公式：
式中：
S35、由式(31)递推可得实体桩桩顶处动力阻抗函数：
式中：
进一步地，所述步骤S4具体求解步骤如下：
S41、根据实体桩桩顶处动力阻抗函数求解实体桩桩顶复刚度，具体为：
式中：Kr代表桩顶动刚度，Ki代表桩顶动阻尼，
由桩顶位移阻抗函数可得桩顶位移频率响应函数为：
S42、桩顶速度频率响应函数为：
根据傅里叶变换的性质，由桩顶速度频率响应函数式(35)可得单位脉冲激励作用下桩顶速度时域响应为：
由卷积定理知，在任意激振力p(t)，桩顶时域速度响应为：
g(t)＝p(t)*h(t)＝IFT[F(iω)×Hv(iω)] (37)
其中，P(iω)为p(t)的傅里叶变换；
S43、半正弦脉冲激振力作用下桩顶速度时域响应半解析解答为：
式中，v(t)为桩顶速度时域响应，桩顶处激励p(t)为半正弦脉冲时，
T为脉冲宽度；
S44、基于求得的所述桩顶速度导纳函数和桩顶速度时域响应函数对桩身振动特性及桩身完整性进行评价。
较现有技术相比，本发明具有以下优点：
本发明基于饱和虚土桩模型的浮承桩纵向振动动力阻抗算法系统，其采用的饱和虚土桩模型能同时考虑桩周、桩底土体饱和特性及桩底土体的波动效应，能适用饱和土中浮承桩纵向振动问题，纵向成层可以考虑土体自然沉积形成的成层特性，通告对桩底土和桩周土的分层处理，同时考虑桩底土层和桩侧土体的纵向成层特性，为桩基动力检测提供理论指导和参考作用。
基于上述理由本发明可在土建技术领域广泛推广。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图做以简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明一种基于虚土桩模型的饱和层状土中桩基纵向振动研究方法流程图。
图2为本发明基于Novak平面应变法及桩底饱和虚土桩模型的纵向振动力学简化模型图。
图3为本发明桩顶动力阻抗函数求解过程具体流程图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分的实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。
需要说明的是，本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
如图1所示，本发明提供了一种基于虚土桩模型的饱和层状土中桩基纵向振动研究方法，包括以下步骤：
S1、构建基于平面应变法及桩底饱和虚土桩模型的饱和土-桩-饱和虚土桩耦合体系纵向振动力学简化模型；
S2、将桩-土耦合振动系统按土体沿纵向分成m层，其中桩底土分成n层，假定桩周和桩底每层土都为均质、各向同性的饱和线黏弹性介质，根据Biot动力波动理论，建立第j层饱和土体和第k层饱和土体的平面应变条件下饱和土层动力控制方程及实体桩与饱和虚土桩界面处的边界条件和桩-土耦合条件；
S3、将第j层饱和土体和第k层饱和土体的谐和激振作用下饱和土层质点位移方程代入所述饱和土层动力学控制方程中，求解桩底第j层饱和土体土骨架径向位移和桩周第k层饱和土体土骨架径向位移，
基于求解的径向位移求解出桩底土与饱和虚土桩界面处、桩周土与实体桩界面处剪应力，通过饱和虚土桩边界条件和各层桩界面处位移连续、力的平衡条件及阻抗函数传递公式，对实体桩桩顶动力阻抗函数进行求解；
S4、由桩顶位移阻抗函数求得桩顶位移频率响应函数，求得桩顶速度频率响应函数、单位脉冲激励作用下桩顶速度时域响应，进而求得桩顶时域速度响应函数，
基于求得的所述桩顶速度频率响应函数和桩顶时域速度响应函数对桩身振动特性及桩身完整性进行评判。
基于Novak平面应变法及桩底饱和虚土桩模型的纵向振动力学简化模型如图2所示，所述步骤S2假定条件还包括：
桩底各层土体为渗透性较差的饱和粘土；
桩周土和桩底土为一系列相互独立的薄层，不考虑土层间的相互作用；
实体桩为均质等截面弹性体，桩身相邻层之间满足力平衡和位移连续条件；
饱和虚土桩为等截面饱和两相介质，饱和虚土桩相邻层之间满足力平衡和位移连续条件，其与实体桩界面位移连续、应力平衡；
桩-土耦合振动系统满足线弹性和小变形条件，桩土界面完全接触，不存在滑移和脱离。
所述步骤S2具体为：
S21、根据Biot动力波动理论，建立平面应变条件下饱和土层动力学控制方程为：
式中，各参数关系如下：
其中，将桩-土耦合振动系统按土体沿纵向分成m层，其中桩底土分成n层，由基岩往上依次编号为1、…、j、…、n、n+1、…、k、…、m层，各层厚度分别为l1、…、lj…、ln、ln+1、…、lk、…、lm，各层顶部埋深分别为h1、…、hj…、hn、hn+1、…、hk、…、hm，j＝1、2、…、n，k＝n+1、n+2、…、m，
式中，uj、wj为桩底第j层土骨架纵向位移和流体相对于土骨架纵向位移，Nj、ηj、bj、kj、Gj、ξj分别为桩底第j层饱和土体流体密度、土颗粒密度、孔隙率、流体粘滞系数、骨架与孔隙流体的粘性耦合系数、Biot定义的动力渗透系数、土体达西定律渗透系数、土体复剪切模量、土体剪切模量和阻尼比；uk、wk为桩周第k层土骨架纵向位移和流体相对于土骨架纵向位移，Nk、ηk、bk、kk、 Gk、ξk分别为桩周第k层饱和土体流体密度、土颗粒密度、孔隙率、流体粘滞系数、骨架与孔隙流体的粘性耦合系数、Biot定义的动力渗透系数、土体达西定律渗透系数、土体复值剪切模量、土体剪切模量和阻尼比；g为重力加速度，
S22、基于Biot动力波动理论得出渗透性较差时饱和土体一维纵向振动控制方程为：
其中，为第j层饱和虚土桩纵向位移，λj为第j层桩底土拉梅常数，且λj＝2νjGj/(1-2νj)，νj为桩底第j层饱和土体泊松比。αj、Mj为Biot定义的表征桩底第j层饱和土体土颗粒及流体压缩性常数， 和分别为桩底第j层饱和土体土颗粒、流体及土骨架的体积压缩模量；
S23、将桩底第j层土在饱和虚土桩界面处的剪切应力代入式(4)中，可得桩底第j层饱和虚土桩纵向振动控制方程为：
式中，为剪应力幅值，ω为激振圆频率，r0表示实体桩半径，Ap表示实体桩的横截面面积；
取第k段实体桩桩身微元体作动力平衡分析可得第k段实体桩纵向振动控制方程为：
式中，为第k段实体桩质点纵向振动位移，EP、ρP分别为实体桩弹性模量和密度，为第k层桩周饱和土体与实体桩界面处土体剪应力，为剪应力幅值。
所述步骤S2中，饱和土-桩-饱和虚土桩体系边界条件如下：
饱和土体：
径向无穷远处位移为零，即
uj(∞,t)＝0 (7a)
uk(∞,t)＝0 (7b)
实体桩：
桩顶平衡条件：
各层界面位移连续条件：
各层界面力平衡条件：
饱和虚土桩：
桩底位移：
各层界面位移连续条件：各层界面力平衡条件：
桩与饱和虚土桩界面处边界条件：
位移连续条件：
力平衡条件：
桩-土耦合条件：
其中，HP为实体桩桩长，表示桩顶作用谐和激振力，表示激振力幅值，H表示基岩上土层总厚度，Hsp表示桩底土层厚，即饱和虚土桩长度。
所述步骤S3中，谐和激振作用下饱和土层质点位移方程具体为：
式中：为桩底第j层饱和土体土骨架径向位移和流体相对于土骨架径向位移响应幅值，为桩周第k层饱和土体土骨架径向位移和流体相对于土骨架径向位移响应幅值，i为虚数单位，
所述步骤S3中，求解纵向位移解具体为：
S31、将谐和激振作用下饱和土层质点位移方程(12)代入所述饱和土层动力学控制方程(1)、(2)中，可得：
将式(13b)、(14b)分别代入式(13a)、(14a)可得：
式中：
方程(15)的通解为：
式中，Aj、Bj、Ak、Bk为待定常数，I0(qjr)、I0(qkr)、K0(qjr)、K0(qkr)为零阶第一类、第二类虚宗量贝塞尔函数。
由边界条件式(7)可知Bj＝Bk＝0，则进一步可得：
如图3所示，所述步骤S3中，对实体桩桩顶动力阻抗函数进行求解具体为：
S32、谐和激振作用下饱和虚土桩和实体桩质点纵向振动位移满足下式：
式中，分别为第j层饱和虚土桩和第k层实体桩质点纵向振动位移响应幅值，则式(5)和式(6)可化简为：
由边界条件式(11)可得：
将式(20)代入式(17)，并由剪应力与位移的关系可得出桩底土与饱和虚土桩界面处、桩周土与实体桩界面处剪应力分别为：
将式(21)代入式(19)，并进一步化简可得：
式中，
则方程(22)的通解为：
式中：Cj、Dj、Ck、Dk为待定常数；
S33、由饱和虚土桩边界条件式(9a)可得：
根据位移阻抗函数的定义可得第1段饱和虚土桩与第2段饱和虚土桩界面处的位移阻抗函数为：
综合第1段饱和虚土桩与第2段饱和虚土桩界面处位移连续、力平衡条件可得：
由此可得第2段饱和虚土桩与第3段饱和虚土桩界面处的位移阻抗函数为：
综合式(24)—(27)可得饱和虚土桩阻抗函数传递公式：
式中：
由式(28)递推可得饱和虚土桩与实体桩界面处阻抗函数：
式中：
S34、根据实体桩与饱和虚土桩界面耦合条件式(10)可得第n+1段实体桩与第n+2段实体桩界面处阻抗函数为：
式中：
由此可得实体桩阻抗函数传递公式：
式中：
S35、由式(31)递推可得实体桩桩顶处动力阻抗函数：
式中：
所述步骤S4具体求解步骤如下：
S41、根据实体桩桩顶处动力阻抗函数求解实体桩桩顶复刚度，具体为：
式中：Kr代表桩顶动刚度，Ki代表桩顶动阻尼，
由桩顶位移阻抗函数可得桩顶位移频率响应函数为：
S42、桩顶速度频率响应函数为：
根据傅里叶变换的性质，由桩顶速度频率响应函数式(35)可得单位脉冲激励作用下桩顶速度时域响应为：
由卷积定理知，在任意激振力p(t)，桩顶时域速度响应为：
g(t)＝p(t)*h(t)＝IFT[F(iω)×Hv(iω)] (37)
其中，P(iω)为p(t)的傅里叶变换；
S43、半正弦脉冲激振力作用下桩顶速度时域响应半解析解答为：
式中，v(t)为桩顶速度时域响应，桩顶处激励p(t)为半正弦脉冲时，
T为脉冲宽度；
S44、基于求得的所述桩顶速度导纳函数和桩顶速度时域响应函数对桩身振动特性及桩身完整性进行评价。
上述本发明实施例序号仅仅为了描述，不代表实施例的优劣。
在本发明的上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述的部分，可以参见其他实施例的相关描述。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (6)
1.一种基于虚土桩模型的饱和层状土中桩基纵向振动研究方法，其特征在于，包括以下步骤：
S1、构建基于平面应变法及桩底饱和虚土桩模型的饱和土-桩-饱和虚土桩耦合体系纵向振动力学简化模型；
S2、将桩-土耦合振动系统按土体沿纵向分成m层，其中桩底土分成n层，假定桩周和桩底每层土都为均质、各向同性的饱和线黏弹性介质，根据Biot动力波动理论，建立第j层饱和土体和第k层饱和土体的平面应变条件下饱和土层动力控制方程及实体桩与饱和虚土桩界面处的边界条件和桩-土耦合条件；
S3、将第j层饱和土体和第k层饱和土体的谐和激振作用下饱和土层质点位移方程代入所述饱和土层动力学控制方程中，求解桩底第j层饱和土体土骨架径向位移和桩周第k层饱和土体土骨架径向位移，
基于求解的径向位移求解出桩底土与饱和虚土桩界面处、桩周土与实体桩界面处剪应力，通过饱和虚土桩边界条件和各层桩界面处位移连续、力的平衡条件及阻抗函数传递公式，对实体桩桩顶动力阻抗函数进行求解；
S4、由桩顶位移阻抗函数求得桩顶位移频率响应函数，求得桩顶速度频率响应函数、单位脉冲激励作用下桩顶速度时域响应，进而求得桩顶时域速度响应函数，
基于求得的所述桩顶速度频率响应函数和桩顶时域速度响应函数对桩身振动特性及桩身完整性进行评判。
2.根据权利要求1所述的基于虚土桩模型的饱和层状土中桩基纵向振动研究方法，其特征在于，所述步骤S2假定条件还包括：
桩底各层土体为渗透性较差的饱和粘土；
桩周土和桩底土为一系列相互独立的薄层，不考虑土层间的相互作用；
实体桩为均质等截面弹性体，桩身相邻层之间满足力平衡和位移连续条件；
饱和虚土桩为等截面饱和两相介质，饱和虚土桩相邻层之间满足力平衡和位移连续条件，其与实体桩界面位移连续、应力平衡；
桩-土耦合振动系统满足线弹性和小变形条件，桩土界面完全接触，不存在滑移和脱离。
3.根据权利要求1所述的基于虚土桩模型的饱和层状土中桩基纵向振动研究方法，其特征在于，所述步骤S2具体为：
S21、根据Biot动力波动理论，建立平面应变条件下饱和土层动力学控制方程为：
式中，各参数关系如下：
其中，将桩-土耦合振动系统按土体沿纵向分成m层，其中桩底土分成n层，由基岩往上依次编号为1、…、j、…、n、n+1、…、k、…、m层，各层厚度分别为l1、…、lj…、ln、ln+1、…、lk、…、lm，各层顶部埋深分别为h1、…、hj…、hn、hn+1、…、hk、…、hm，j＝1、2、…、n，k＝n+1、n+2、…、m，
式中，uj、wj为桩底第j层土骨架纵向位移和流体相对于土骨架纵向位移，Nj、ηj、bj、kj、Gj、ξj分别为桩底第j层饱和土体流体密度、土颗粒密度、孔隙率、流体粘滞系数、骨架与孔隙流体的粘性耦合系数、Biot定义的动力渗透系数、土体达西定律渗透系数、土体复剪切模量、土体剪切模量和阻尼比；uk、wk为桩周第k层土骨架纵向位移和流体相对于土骨架纵向位移，Nk、ηk、bk、kk、 Gk、ξk分别为桩周第k层饱和土体流体密度、土颗粒密度、孔隙率、流体粘滞系数、骨架与孔隙流体的粘性耦合系数、Biot定义的动力渗透系数、土体达西定律渗透系数、土体复值剪切模量、土体剪切模量和阻尼比；g为重力加速度，
S22、基于Biot动力波动理论得出渗透性较差时饱和土体一维纵向振动控制方程为：
其中，为第j层饱和虚土桩纵向位移，λj为第j层桩底土拉梅常数，且λj＝2νjGj/(1-2νj)，νj为桩底第j层饱和土体泊松比。αj、Mj为Biot定义的表征桩底第j层饱和土体土颗粒及流体压缩性常数， 和分别为桩底第j层饱和土体土颗粒、流体及土骨架的体积压缩模量；
S23、将桩底第j层土在饱和虚土桩界面处的剪切应力代入式(4)中，可得桩底第j层饱和虚土桩纵向振动控制方程为：
式中，为剪应力幅值，ω为激振圆频率，r0表示实体桩半径，Ap表示实体桩的横截面面积；
取第k段实体桩桩身微元体作动力平衡分析可得第k段实体桩纵向振动控制方程为：
式中，为第k段实体桩质点纵向振动位移，EP、ρP分别为实体桩弹性模量和密度，为第k层桩周饱和土体与实体桩界面处土体剪应力，为剪应力幅值。
4.根据权利要求3所述的基于虚土桩模型的饱和层状土中桩基纵向振动研究方法，其特征在于，所述步骤S2中，饱和土-桩-饱和虚土桩体系边界条件如下：
饱和土体：
径向无穷远处位移为零，即
uj(∞,t)＝0 (7a)
uk(∞,t)＝0 (7b)
实体桩：
桩顶平衡条件：
各层界面位移连续条件：
各层界面力平衡条件：
饱和虚土桩：
桩底位移：
各层界面位移连续条件：
各层界面力平衡条件：
桩与饱和虚土桩界面处边界条件：
位移连续条件：
力平衡条件：
桩-土耦合条件：
其中，HP为实体桩桩长，表示桩顶作用谐和激振力，表示激振力幅值，H表示基岩上土层总厚度，Hsp表示桩底土层厚，即饱和虚土桩长度。
5.根据权利要求4所述的基于虚土桩模型的饱和层状土中桩基纵向振动研究方法，其特征在于，所述步骤S3中，谐和激振作用下饱和土层质点位移方程具体为：
式中：为桩底第j层饱和土体土骨架径向位移和流体相对于土骨架径向位移响应幅值，为桩周第k层饱和土体土骨架径向位移和流体相对于土骨架径向位移响应幅值，i为虚数单位，
所述步骤S3中，求解纵向位移解具体为：
S31、将谐和激振作用下饱和土层质点位移方程(12)代入所述饱和土层动力学控制方程(1)、(2)中，可得：
将式(13b)、(14b)分别代入式(13a)、(14a)可得：
式中：
方程(15)的通解为：
式中，Aj、Bj、Ak、Bk为待定常数，I0(qjr)、I0(qkr)、K0(qjr)、K0(qkr)为零阶第一类、第二类虚宗量贝塞尔函数。
由边界条件式(7)可知Bj＝Bk＝0，则进一步可得：
所述步骤S3中，对实体桩桩顶动力阻抗函数进行求解具体为：
S32、谐和激振作用下饱和虚土桩和实体桩质点纵向振动位移满足下式：
式中，分别为第j层饱和虚土桩和第k层实体桩质点纵向振动位移响应幅值，则式(5)和式(6)可化简为：
由边界条件式(11)可得：
将式(20)代入式(17)，并由剪应力与位移的关系可得出桩底土与饱和虚土桩界面处、桩周土与实体桩界面处剪应力分别为：
将式(21)代入式(19)，并进一步化简可得：
式中，
则方程(22)的通解为：
式中：Cj、Dj、Ck、Dk为待定常数；
S33、由饱和虚土桩边界条件式(9a)可得：
根据位移阻抗函数的定义可得第1段饱和虚土桩与第2段饱和虚土桩界面处的位移阻抗函数为：
综合第1段饱和虚土桩与第2段饱和虚土桩界面处位移连续、力平衡条件可得：
由此可得第2段饱和虚土桩与第3段饱和虚土桩界面处的位移阻抗函数为：
综合式(24)—(27)可得饱和虚土桩阻抗函数传递公式：
式中：
由式(28)递推可得饱和虚土桩与实体桩界面处阻抗函数：
式中：
S34、根据实体桩与饱和虚土桩界面耦合条件式(10)可得第n+1段实体桩与第n+2段实体桩界面处阻抗函数为：
式中：
由此可得实体桩阻抗函数传递公式：
式中：
S35、由式(31)递推可得实体桩桩顶处动力阻抗函数：
式中：
6.根据权利要求5所述的基于虚土桩模型的饱和层状土中桩基纵向振动研究方法，其特征在于，所述步骤S4具体求解步骤如下：
S41、根据实体桩桩顶处动力阻抗函数求解实体桩桩顶复刚度，具体为：
式中：Kr代表桩顶动刚度，Ki代表桩顶动阻尼，
由桩顶位移阻抗函数可得桩顶位移频率响应函数为：
S42、桩顶速度频率响应函数为：
根据傅里叶变换的性质，由桩顶速度频率响应函数式(35)可得单位脉冲激励作用下桩顶速度时域响应为：
由卷积定理知，在任意激振力p(t)，桩顶时域速度响应为：
g(t)＝p(t)*h(t)＝IFT[F(iω)×Hv(iω)] (37)
其中，P(iω)为p(t)的傅里叶变换；
S43、半正弦脉冲激振力作用下桩顶速度时域响应半解析解答为：
式中，v(t)为桩顶速度时域响应，桩顶处激励p(t)为半正弦脉冲时，
T为脉冲宽度；
S44、基于求得的所述桩顶速度导纳函数和桩顶速度时域响应函数对桩身振动特性及桩身完整性进行评价。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811217886.XA CN109344526B (zh) | 2018-10-18 | 2018-10-18 | 一种基于虚土桩模型的饱和层状土中桩基纵向振动研究方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811217886.XA CN109344526B (zh) | 2018-10-18 | 2018-10-18 | 一种基于虚土桩模型的饱和层状土中桩基纵向振动研究方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109344526A true CN109344526A (zh) | 2019-02-15 |
| CN109344526B CN109344526B (zh) | 2022-09-23 |
Family
ID=65308777
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811217886.XA Active CN109344526B (zh) | 2018-10-18 | 2018-10-18 | 一种基于虚土桩模型的饱和层状土中桩基纵向振动研究方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109344526B (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110147631A (zh) * | 2019-05-29 | 2019-08-20 | 大连海事大学 | 一种径向非均质饱和土中实体桩水平振动分析方法 |
| CN110162908A (zh) * | 2019-05-29 | 2019-08-23 | 大连海事大学 | 一种径向非均质饱和土中管桩水平振动分析方法 |
| CN111962572A (zh) * | 2020-08-07 | 2020-11-20 | 河海大学 | 一种部分埋入式桩基的水平-摇摆振动阻抗确定方法 |
| CN112302061A (zh) * | 2020-09-29 | 2021-02-02 | 安徽省（水利部淮河水利委员会）水利科学研究院（安徽省水利工程质量检测中心站） | 低应变基桩完整性检测信号的智能化快速解释方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1323869A1 (en) * | 2001-12-27 | 2003-07-02 | Mitsubishi Denki Kabushiki Kaisha | Dynamic test method for bearing capacity of piles |
| CN104596852A (zh) * | 2014-11-25 | 2015-05-06 | 宁波大学 | 一种岩土体温控动力特性试验系统及试验方法 |
| CN108107190A (zh) * | 2018-01-29 | 2018-06-01 | 铜陵学院 | 一种饱水土体局部振动液化试验装置和方法 |
| CN108446460A (zh) * | 2018-03-01 | 2018-08-24 | 大连海事大学 | 一种径向非均质、纵向成层土体中桩基纵向振动分析方法 |
-
2018
- 2018-10-18 CN CN201811217886.XA patent/CN109344526B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1323869A1 (en) * | 2001-12-27 | 2003-07-02 | Mitsubishi Denki Kabushiki Kaisha | Dynamic test method for bearing capacity of piles |
| CN104596852A (zh) * | 2014-11-25 | 2015-05-06 | 宁波大学 | 一种岩土体温控动力特性试验系统及试验方法 |
| CN108107190A (zh) * | 2018-01-29 | 2018-06-01 | 铜陵学院 | 一种饱水土体局部振动液化试验装置和方法 |
| CN108446460A (zh) * | 2018-03-01 | 2018-08-24 | 大连海事大学 | 一种径向非均质、纵向成层土体中桩基纵向振动分析方法 |
Non-Patent Citations (1)
| Title |
|---|
| 何伟杰等: "基于扩散虚土桩法的大直径桩纵向振动研究", 《地震工程与工程振动》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110147631A (zh) * | 2019-05-29 | 2019-08-20 | 大连海事大学 | 一种径向非均质饱和土中实体桩水平振动分析方法 |
| CN110162908A (zh) * | 2019-05-29 | 2019-08-23 | 大连海事大学 | 一种径向非均质饱和土中管桩水平振动分析方法 |
| CN110162908B (zh) * | 2019-05-29 | 2022-09-27 | 大连海事大学 | 一种径向非均质饱和土中管桩水平振动分析方法 |
| CN110147631B (zh) * | 2019-05-29 | 2022-09-27 | 大连海事大学 | 一种径向非均质饱和土中实体桩水平振动分析方法 |
| CN111962572A (zh) * | 2020-08-07 | 2020-11-20 | 河海大学 | 一种部分埋入式桩基的水平-摇摆振动阻抗确定方法 |
| CN111962572B (zh) * | 2020-08-07 | 2021-08-24 | 河海大学 | 一种部分埋入式桩基的水平-摇摆振动阻抗确定方法 |
| CN112302061A (zh) * | 2020-09-29 | 2021-02-02 | 安徽省（水利部淮河水利委员会）水利科学研究院（安徽省水利工程质量检测中心站） | 低应变基桩完整性检测信号的智能化快速解释方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109344526B (zh) | 2022-09-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109344526B (zh) | 一种基于虚土桩模型的饱和层状土中桩基纵向振动研究方法 | |
| CN109446637B (zh) | 一种基于层状饱和虚土桩的浮承桩纵向振动分析方法 | |
| Ecemis | Simulation of seismic liquefaction: 1-g model testing system and shaking table tests | |
| Bendoni et al. | On salt marshes retreat: Experiments and modeling toppling failures induced by wind waves | |
| CN109214123B (zh) | 一种基于饱和虚土桩模型的浮承桩纵向振动分析方法 | |
| Sawicki et al. | Developments in modeling liquefaction of granular soils, caused by cyclic loads | |
| CN108446460B (zh) | 一种径向非均质、纵向成层土体中桩基纵向振动分析方法 | |
| Li et al. | Wave-induced dynamic response in a transversely isotropic and multilayered poroelastic seabed | |
| CN109359390B (zh) | 一种轴对称双向非均质黏性阻尼土中桩基扭转振动分析方法 | |
| Khalil et al. | Dynamic behavior of pile foundations under vertical and lateral vibrations | |
| Kumhar et al. | Green’s function and surface waves in a viscoelastic orthotropic FGM enforced by an impulsive point source | |
| CN105869208A (zh) | 三维饱和沉积盆地地震动模拟方法 | |
| Dhar et al. | Different approaches for numerical modeling of seismic soil-structure interaction: impacts on the seismic response of a simplified reinforced concrete integral bridge | |
| Haigh et al. | Importance of time step discretisation for nonlinear dynamic finite element analysis | |
| CN110093951B (zh) | 一种基于虚土桩模型的摩擦桩纵向振动分析方法 | |
| CN108842831B (zh) | 层状横观黏弹性土中部分埋入群桩的竖向振动分析方法 | |
| CN109372035B (zh) | 一种基于饱和虚土桩的浮承桩纵向振动分析方法 | |
| CN110219324B (zh) | 一种成层土中非完全黏结的摩擦桩纵向振动分析方法 | |
| Hu et al. | Dynamic torsional response of pre-strained end bearing pile embedded in pre-strained isotropic saturated soil medium | |
| Kwon et al. | 3D dynamic numerical modeling for soil-pile-structure interaction in centrifuge tests | |
| Lai et al. | Numerical modeling of dynamic centrifuge experiments on a saturated dense sand stratum | |
| Li et al. | Tide-induced airflow in a two-layered coastal land with atmospheric pressure fluctuations | |
| Gupta et al. | Effect of inhomogeneity, initial stress, sandiness, and viscosity on the propagation of torsional wave | |
| Kudarova | Effective models for seismic wave propagation in porous media | |
| Sawada et al. | Applicability of the generalised scaling law to pile-inclined ground system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |