[发明专利]层状横观黏弹性土中部分埋入群桩的竖向振动分析方法

摘要

本发明公开了层状横观黏弹性土中部分埋入群桩的竖向振动分析方法，属于地基工程试验分析技术领域。本方法的de Boer基于连续介质混合物公理和体积分数概念建立的多孔介质理论，利用体积分数的概念，若干的微观性质可以直接通过宏观性质来描述，并且避免了杂交混合物理论中的繁杂公式。de Boer建立的多孔介质理论不仅符合连续介质力学和热动力学原理，并且具有许多Biot理论不具备的优点。此外，已有饱和土‑群桩动力相互作用的研究大多采用Dobry和Gazetas提出的桩‑土‑桩动力相互作用因子。而该桩‑桩动力相互作用因子为单相土的，对于饱和土问题其适用性尚需进一步研究。可为更为广义和复杂的群桩动力特性问题的研究提供理论指导和参考作用。

主权项

1.层状横观黏弹性土中部分埋入群桩的竖向振动分析方法，其特征在于：该方法的实现步骤如下，1.1 基于Boer多孔介质模型的饱和土体动力控制方程组根据横观各向同性饱和黏弹性运动方程描述桩周土的动力学行为；多层土中的第k层土的运动方程表示为式1.1‑1.5中所有变量中带字母k均表示该变量属于第k层土；和分别表示土骨架和孔隙流体的体积分数，且满足饱和条件和分别表示土骨架和孔隙流体的体积密度；u_sk和u_fk分别表示土骨架和孔隙流体的径向位移；w_sk和w_fk分别表示土骨架和孔隙流体的竖向位移；表示孔隙流体压力；液固耦合系数其中g为重力加速度，而和分别表示r和z方向上的Darcy渗透系数；另外，其中和分别表示水平和竖向弹性模量，表示竖向平面的剪切模量，为由水平应力导致的竖向应变的泊松比，为竖向应力引起的水平应变的泊松比，并且为水平应力引起的正交向水平应变的泊松比；τ_εk和τ_σk为表示土体黏性的参数；为α_k阶的Riemann‑Liouville分数导数，且其中为Gamma函数，0＜α_k＜1；1.2 桩‑桩动力相互作用因子半径r₀，长度L的黏弹性圆柱桩1部分埋入一厚度L₂，其中未埋入段L₁，埋入段L₂，且L＝L₁+L₂，下卧基岩的层状横观各向同性饱和黏弹性土中，同时在黏弹性圆柱桩1的桩顶处受到一竖向谐和激振力N(t)＝N₁e^iωt作用，由于黏弹性圆柱桩1受荷振动引起土中产生径向传播的剪切波，进而使得距离S处的黏弹性圆柱桩2受迫振动；黏弹性圆柱桩2和黏弹性圆柱桩1的物理几何性质相同，黏弹性圆柱桩2和黏弹性圆柱桩1的弹性模量为E_p，密度ρ_p，泊松比υ_p，阻尼系数η_p；k_s表示桩周土体的竖向动力阻抗，表示饱和土竖向位移的衰减函数，w_p1和w_p2分别表示黏弹性圆柱桩1即主动桩1和黏弹性圆柱桩2即被动桩2的竖向位移；假设桩土体系的振动为小变形、小位移振动，且桩土在振动过程中保持紧密接触，即在桩土接触面处的位移、应力连续；采用一系列的Voigt模型模拟相邻土层间的动力相互作用，其弹簧和阻尼器沿土体径向均匀分布；引入量分别表示在k层土顶部Voigt模型的刚度系数和阻尼系数，分别表示在k层土底部Voigt模型的刚度系数和阻尼系数；所有量中的上标或者下标k表示第k层土；由于假设了相邻土层间位移和应力连续，故有忽略土体的径向位移，仅考虑其竖向波动效应，即设usk＝ufk＝0，则式(1.1)～(1.5)被简化为在稳态振动问题中设w_sk＝W_ske^iωt，w_fk＝W_fke^iωt，w_pk＝W_pke^iωt，将上述变量代入式(1.6)～(1.8)得为分析方便，引入如下无量纲参数和变量：其中取G＝G_s1，分别表示第一层土的剪切模量和密度；将无量纲参数和变量代入式(1.9)～(1.11)后整理得联立式(1.12)和(1.13)有式中由式(1.14)得令式(1.15)改写为根据变量分离法设并将其代入式(1.17)中得式中Real(q_k)＞0，Real(g_k)＞0；求解式(1.18)和(1.19)后得则在r方向无穷远处土体位移应衰减为零，考虑到修正Bessel函数的特性，有式中Ak＝CkEk，Bk＝DkEk；根据前述变量及的关系，并再次考虑上述边界条件得则相应的剪应力为桩周土体反力f_sk为剪应力沿桩周长的积分，即定义桩周土竖向复刚度为以及k层土的位移衰减函数为式中为两桩的无量纲间距；根据Rayleigh‑Love杆理论，建立主动桩1未嵌入段的振动方程为及其轴力为将式(1.30)无量纲化为式中式(1.29)的无量纲形式为式中Real(λ₁)＞0；式(1.32)的通解为式中a1和a2为待定系数；由式(1.31)得轴力为式中联立式(1.33)和(1.34)，将主动桩1未嵌入段两端的位移及轴力以矩阵形式统一写为式中建立第k层土中主动桩1的振动方程为将其无量纲化为式(1.37)的通解为式中采用矩阵写法，第k层土桩段两端的位移及轴力关系为式中且l_k为第k层土桩段的长度；根据式(1.35)和(1.39)，主动桩1两端的位移及轴力关系为式中[T¹]＝[T_1n]…[T_1k]…[T₁₁][T₀]，元素形式为主动桩1底部满足边界条件则此为主动桩1顶部位移和轴力之间的关系；为单桩桩顶的竖向动力阻抗；被动桩2未嵌入段的振动方程与主动桩1未嵌入段的一致，则被动桩2未嵌入段两端的位移和轴力关系表达为继而建立被动桩2嵌入段的振动方程为将其无量纲化为式中式(1.45)的通解为式中将式(1.38)代入式(1.46)中得桩2的轴力为联立式(1.47)和(1.48)，被动桩2的第k段两端的位移和轴力关系为式中D_42k＝πζλ_3k，D_41k＝πζλ_2k；考虑式(1.43)和(1.49)，被动桩2整体两端的位移和轴力关系为式中[T2]＝[T2n]…[T2k]…[T21][T0]，[T3]＝[T3n]…[T3k]…[T31][T0]；被动桩2满足在桩顶在桩底即将主动桩1顶部的动力阻抗代入式(1.51)中并定义桩‑桩动力相互作用因子为式中a21即表示所研究饱和土‑群桩振动问题的桩‑桩动力相互作用因子；至此，仅剩未知参数gk，根据相邻土层间的平衡条件求得；在k层土顶面满足在该k层土底面满足将式(1.53)及(1.54)无量纲化为式中将式(1.24)代入式(1.55)及(1.56)中得当2≤k≤n‑1时，gk满足式中当k＝1时，由于β2k＝0，则gk满足当k＝n时，由于β1k＝∞，则gk满足1.3 群桩动力阻抗假设一i×j形式的群桩桩顶整体被刚性连接，且在中心处受到竖向激振力N(t)＝N0eiωt作用，采用所得桩‑桩动力相互作用因子，则该群桩满足如下连续性条件：式中W0为群桩顶竖向位移幅值，Wi为桩i的竖向位移幅值；Nj为桩j顶的轴力幅值；aij为桩j对桩i的动力相互作用因子，且当i＝j时aij＝1；为表示方便，上式(1.60)及(1.61)中的时间项eiωt已被省略；将式(1.60)和(1.61)写成矩阵形式为定义群桩竖向动力阻抗为式中KV表示群桩刚度，CV表示群桩阻尼；联立式(1.62)和(1.63)求得KG；为方便分析，将KG标准化为式中m表示单桩数量，表示单桩顶部的静阻抗；基于桩顶频域响应函数，探讨桩基非嵌入长度、土体各向异性及其液固耦合系数因素对群桩动力特性的影响规律。