专利名称: |
层状横观黏弹性土中部分埋入群桩的竖向振动分析方法 |
摘要: |
本发明公开了层状横观黏弹性土中部分埋入群桩的竖向振动分析方法,属于地基工程试验分析技术领域。本方法的de Boer基于连续介质混合物公理和体积分数概念建立的多孔介质理论,利用体积分数的概念,若干的微观性质可以直接通过宏观性质来描述,并且避免了杂交混合物理论中的繁杂公式。de Boer建立的多孔介质理论不仅符合连续介质力学和热动力学原理,并且具有许多Biot理论不具备的优点。此外,已有饱和土‑群桩动力相互作用的研究大多采用Dobry和Gazetas提出的桩‑土‑桩动力相互作用因子。而该桩‑桩动力相互作用因子为单相土的,对于饱和土问题其适用性尚需进一步研究。可为更为广义和复杂的群桩动力特性问题的研究提供理论指导和参考作用。 |
专利类型: |
发明专利 |
国家地区组织代码: |
北京;11 |
申请人: |
北京工业大学 |
发明人: |
许成顺;崔春义;张石平;梁志孟;孟坤 |
专利状态: |
有效 |
发布日期: |
2019-01-01T00:00:00+0800 |
申请号: |
CN201810536190.7 |
公开号: |
CN108842831A |
代理机构: |
北京思海天达知识产权代理有限公司 11203 |
代理人: |
沈波 |
分类号: |
E02D33/00(2006.01)I;E;E02;E02D;E02D33;E02D33/00 |
申请人地址: |
100124 北京市朝阳区平乐园100号 |
主权项: |
1.层状横观黏弹性土中部分埋入群桩的竖向振动分析方法,其特征在于:该方法的实现步骤如下,1.1 基于Boer多孔介质模型的饱和土体动力控制方程组根据横观各向同性饱和黏弹性运动方程描述桩周土的动力学行为;多层土中的第k层土的运动方程表示为![]() ![]() ![]() ![]() 式1.1‑1.5中所有变量中带字母k均表示该变量属于第k层土; 和 分别表示土骨架和孔隙流体的体积分数,且满足饱和条件![]() 和 分别表示土骨架和孔隙流体的体积密度;usk和ufk分别表示土骨架和孔隙流体的径向位移;wsk和wfk分别表示土骨架和孔隙流体的竖向位移; 表示孔隙流体压力;液固耦合系数 其中g为重力加速度,而 和 分别表示r和z方向上的Darcy渗透系数;另外,![]() ![]() 其中 和 分别表示水平和竖向弹性模量, 表示竖向平面的剪切模量, 为由水平应力导致的竖向应变的泊松比, 为竖向应力引起的水平应变的泊松比,并且![]() 为水平应力引起的正交向水平应变的泊松比;τεk和τσk为表示土体黏性的参数; 为αk阶的Riemann‑Liouville分数导数,且 其中 为Gamma函数,0<αk<1;1.2 桩‑桩动力相互作用因子半径r0,长度L的黏弹性圆柱桩1部分埋入一厚度L2,其中未埋入段L1,埋入段L2,且L=L1+L2,下卧基岩的层状横观各向同性饱和黏弹性土中,同时在黏弹性圆柱桩1的桩顶处受到一竖向谐和激振力N(t)=N1eiωt作用, 由于黏弹性圆柱桩1受荷振动引起土中产生径向传播的剪切波,进而使得距离S处的黏弹性圆柱桩2受迫振动;黏弹性圆柱桩2和黏弹性圆柱桩1的物理几何性质相同,黏弹性圆柱桩2和黏弹性圆柱桩1的弹性模量为Ep,密度ρp,泊松比υp,阻尼系数ηp;ks表示桩周土体的竖向动力阻抗, 表示饱和土竖向位移的衰减函数,wp1和wp2分别表示黏弹性圆柱桩1即主动桩1和黏弹性圆柱桩2即被动桩2的竖向位移;假设桩土体系的振动为小变形、小位移振动,且桩土在振动过程中保持紧密接触,即在桩土接触面处的位移、应力连续;采用一系列的Voigt模型模拟相邻土层间的动力相互作用,其弹簧和阻尼器沿土体径向均匀分布;引入量 分别表示在k层土顶部Voigt模型的刚度系数和阻尼系数, 分别表示在k层土底部Voigt模型的刚度系数和阻尼系数;所有量中的上标或者下标k表示第k层土;由于假设了相邻土层间位移和应力连续,故有![]() 忽略土体的径向位移,仅考虑其竖向波动效应,即设usk=ufk=0,则式(1.1)~(1.5)被简化为![]() ![]() 在稳态振动问题中设wsk=Wskeiωt,wfk=Wfkeiωt,wpk=Wpkeiωt,![]() 将上述变量代入式(1.6)~(1.8)得![]() ![]() 为分析方便,引入如下无量纲参数和变量:![]() ![]() 其中取G=Gs1, 分别表示第一层土的剪切模量和密度;将无量纲参数和变量代入式(1.9)~(1.11)后整理得![]() ![]() 联立式(1.12)和(1.13)有 式中![]() 由式(1.14)得 令 式(1.15)改写为 根据变量分离法设 并将其代入式(1.17)中得![]() 式中 Real(qk)>0,Real(gk)>0;求解式(1.18)和(1.19)后得![]() 则 在r方向无穷远处土体位移应衰减为零,考虑到修正Bessel函数的特性,有 式中Ak=CkEk,Bk=DkEk;根据前述变量 及 的关系,并再次考虑上述边界条件得 则相应的剪应力为 桩周土体反力fsk为剪应力 沿桩周长的积分,即 定义桩周土竖向复刚度为 以及k层土的位移衰减函数为 式中 为两桩的无量纲间距;根据Rayleigh‑Love杆理论,建立主动桩1未嵌入段的振动方程为 及其轴力为 将式(1.30)无量纲化为 式中 式(1.29)的无量纲形式为 式中 Real(λ1)>0;式(1.32)的通解为 式中a1和a2为待定系数;由式(1.31)得轴力为 式中 联立式(1.33)和(1.34),将主动桩1未嵌入段两端的位移及轴力以矩阵形式统一写为 式中 建立第k层土中主动桩1的振动方程为 将其无量纲化为 式(1.37)的通解为 式中 采用矩阵写法,第k层土桩段两端的位移及轴力关系为 式中 且lk为第k层土桩段的长度; 根据式(1.35)和(1.39),主动桩1两端的位移及轴力关系为 式中 [T1]=[T1n]…[T1k]…[T11][T0],元素形式为 主动桩1底部满足边界条件 则![]() 此为主动桩1顶部位移和轴力之间的关系; 为单桩桩顶的竖向动力阻抗;被动桩2未嵌入段的振动方程与主动桩1未嵌入段的一致,则被动桩2未嵌入段两端的位移和轴力关系表达为 继而建立被动桩2嵌入段的振动方程为 将其无量纲化为 式中 式(1.45)的通解为 式中 将式(1.38)代入式(1.46)中得 桩2的轴力为 联立式(1.47)和(1.48),被动桩2的第k段两端的位移和轴力关系为 式中D42k=πζλ3k,D41k=πζλ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表示单桩数量, 表示单桩顶部的静阻抗;基于桩顶频域响应函数,探讨桩基非嵌入长度、土体各向异性及其液固耦合系数因素对群桩动力特性的影响规律。 |
所属类别: |
发明专利 |