V型山谷地形中桩基的动力响应计算方法技术

本发明专利技术涉及一种V型山谷地形中桩基的动力响应计算方法，包括如下步骤：(1)获取山谷、桩基及土体的物理力学参数；(2)基于Graf加法定理计算得出由于地震波震动引起的土体位移场；(3)基于动力Winkler地基模型，将地震波震动引起的土体位移场转化为荷载施加到桩基上，建立桩基的受力平衡方程；(4)在桩顶、桩底均处于自由约束的情况下时，求解均布荷载引起的桩基位移及其所受的弯矩，剪力。上述方法通过结合Graf加法定理与Winkler地基梁模型，首次提出了V型山谷上桩基础震动响应闭合解，并准确揭示V型山谷地震放大效应对桩基础震动响应的影响。响。响。

【技术实现步骤摘要】
V型山谷地形中桩基的动力响应计算方法


[0001]本专利技术涉及抗震设计
，更具体地涉及一种V型山谷地形中桩基的动力响应计算方法。

技术介绍

[0002]我国国土面积辽阔，但也是个多山国家，山地、丘陵、高原的面积占了国土面积的2/3以上。因此，使得我国很多地区都存在泥石流、滑坡等风险，尤其是在地震的作用下，会加大这类风险的发生几率，给路桥建设工程带来了很大的难度，影响了V型山谷桩基的稳定性。对此，正确分析地震作用下山区V型山谷桩基动力响应情况，有利于了解山区V型山谷桩基的实际情况，从而制定正确的施工计划，有利于提高桩基质量，保障人们的生命财产安全，实现其应有的效益。
[0003]相比于传统桩基的设计方法来说，山区V型山谷桩基承载设计要更为复杂，施工难度高。我国关于这方面的研究比较少。因此本文提出一种地震波作用下，V型山谷地形中桩基的位移计算方法，从而更加精准的计算出桩基的位移，防范工程事故的发生。现有技术对于V型山谷桩基础震动响应研究存在空白，目前国际上对V型山谷桩震动响应研究为通过极限法模拟震后V型山谷的残余变形，并将其转化为位移荷载计算V型山谷桩的静力受力与变形。这一方法存在着以下问题：
[0004]1.仅能模拟震后V型山谷残余位移对桩基础受力变形影响，而桩基础抗震安全性最低时往往处于地震过程而不是震后，这导致基于该方法的V型山谷桩基抗震设计理论无法计算V型山谷桩基抗震最不利情况。
[0005]2.V型山谷地形效应在地震过程中会显著增幅桩基础震动响应，而现有研究理论无法揭示这种增幅规律，基于现有理论进行V型山谷桩基础抗震设计会导致安全冗余度严重降低。
[0006]3.无法模拟地震过程中桩基础的震动规律，无法获得地震过程中桩基础的动力响应时程曲线。

技术实现思路

[0007]为了解决现有技术中V型山谷桩基础震动响应研究存在空白、极限法模拟震后V型山谷的残余变形存在的以上技术问题。本专利技术公开了一种V型山谷地形中桩基的动力响应计算方法，包含如下步骤：
[0008](1)根据地勘资料获取V型山谷、桩基及土体的物理力学参数；
[0009](2)基于Graf加法定理计算得出由于地震波震动引起的土体位移场；
[0010](3)基于动力Winkler地基模型，将地震波引起的土体位移场转化为荷载，施加到桩基上，建立桩基的受力平衡方程；
[0011](4)在桩顶、桩底均处于自由约束的情况下时，求解均布荷载引起的桩基位移及其所受的弯矩，剪力。
[0012]所述物理力学参数包括：
[0013](1)山谷参数：山谷半径d，地震波无量纲频率为η，入射角角度α，峡谷的中心距离两边的距离b1，b2；
[0014](2)桩基参数：等效直径D，入土深度L，桩基弹性模量E
p
，桩惯性矩I
p
；
[0015](3)土体参数：弹性模量E
s
，土体泊松比μ
s
。
[0016]所述土体外域位移场u
outer
(r,θ)计算方法为：
[0017]u
outer
(r，θ)＝u
f
(r，θ)+u
s1
(r，θ)+u
s2
(r，θ)；
ꢀꢀ
(1)
[0018]其中，外域自由场u
f
(r，θ)，外域散射场1u
s1
(r，θ)，外域散射场2u
s2
(r，θ)的计算方法分别为：
[0019][0020][0021][0022]内域散射场u
inner1
(r，θ)与内域散射场u
inner2
(r，θ)的计算公式分别为：
[0023][0024][0025][0026][0027][0028][0029]其中，α为地震波入射角度；r是以o点为圆心的极半径，θ为极角；k为入射波波数，其中η为入射波无量纲频率；ε
n
为纽曼系数(ε0＝1，n＝0；ε
n
＝2，n>0)；H
(1)
(
·
)是n阶第三类贝塞尔函数，J
n
(
·
)是n阶第一类贝塞尔函数；A
n
，B
n
，C
n
，D
n
是要确定的系数，m和n代表级数，v代表V型山谷参数。
[0030]所述桩基的受力平衡方程如下：
[0031][0032]其中，W(y)为桩基的位移函数；E
p
为桩基弹性模量；I
p
为桩基横截面惯性矩；D为等效宽度；k1为桩侧土的弹簧刚度；E
s
为土体弹性模量；μ
s
为土体泊松比；P(y)为桩基受到的竖向外荷载。
[0033][0034]P(y)＝k1×
w
f
(y)；
ꢀꢀ
(13)
[0035][0036]式中，w
f
(y)为土体自由场位移；x为桩基放置的位置。
[0037]求解集中荷载引起的桩基剪力Q和弯矩M如下：
[0038]可将上述方程写成有限差分的形式：
[0039]αw
i
‑2+βw
i
‑1+γw
i
+βw
i+1
+αw
i+2
＝p
i
；
ꢀꢀ
(15)
[0040]式中，
[0041]结合边界条件可推导出单桩位移的方程：
[0042][0043]此处的K可以看做一个从n+1*n+1的矩阵；
[0044]而刚度矩阵如下：
[0045][0046]采用有限差分法求解上式。将桩基划分为n段，每段的高度为H＝L/n，L为桩长。桩基节点编号从桩顶至桩端依次为0，1，
…
，n
‑
1，n。计算时分别在桩顶和桩端处增加两个虚拟节点
‑
2，
‑
1，n+1，n+2。
[0047]桩身弯矩如下：
[0048][0049]桩身剪力如下：
[0050][0051]本专利技术的技术方案有益效果在于：首次采用动力Winkler地基模型结合波函数求解了V型山谷桩基础震动响应，获得了V型山谷桩基地震过程中得动力时程曲线，并能揭示V
型山谷地形效应对桩震动响应的增幅规律。
附图说明
[0052]以下，通过示例的方式示出本专利技术的示例性实施例的附图，各附图中使用了相同或相似的附图标记来表示相同或相似的元素。附图中：
[0053]图1是本专利技术所述方法的流程图。
[0054]图2是本专利技术所述方法的计算模型图。
[0055]图3是本专利技术实施例中基于地形效应的V型山谷单桩自由场位移曲线本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种V型山谷地形中桩基的动力响应计算方法，其特征在于，包括如下步骤：(1)获取V型山谷、桩基及土体的物理力学参数；所述物理力学参数包括：山谷参数：山谷半径d，地震波无量纲频率为η，入射角角度α，峡谷的中心距离两边的距离b1，b2；桩基参数：等效直径D，入土深度L，桩基弹性模量E
p
，桩惯性矩I
p
；土体参数：弹性模量E
s
，土体泊松比μ
s
；(2)基于Graf加法定理计算得出由于地震波震动引起的土体位移场；(3)基于动力Winkler地基模型，将地震波震动引起的土体位移场转化为荷载，施加到桩基上，建立桩基的受力平衡方程；(4)在桩顶、桩底均处于自由约束的情况下时，求解均布荷载引起的桩基位移及其所受的弯矩，剪力。2.如权利要求1所述的V型山谷地形中桩基的动力响应计算方法，其特征在于，在步骤(2)中，所述由于地震波震动引起的土体位移场包括地震波入射形成的外域位移场，地震波入射形成的外域位移场的计算方法为：u
outer
(r,θ)＝u
f
(r,θ)+u
s1
(r,θ)+u
s2
(r,θ)；
ꢀꢀꢀꢀ
(1)其中，外域自由场u
f
(r,θ)，外域散射场u
s1
(r,θ)，外域散射场u
s2
(r,θ)的计算公式分别为：为：为：式中：α为地震波入射角度；r是以o点为圆心的极半径，θ为极角；k为入射波波数，其中η为入射波无量纲频率；ε
n
为纽曼系数(ε0＝1，n＝0；ε
n
＝2，n＞0)；H
(1)
(
·
)是n阶第三类贝塞尔函数，J
n
(
·
)是n阶第一类贝塞尔函数；A
n
...

【专利技术属性】
技术研发人员：可文海，李源，管凌霄，吴磊，杨文海，丁海滨，岳成，
申请(专利权)人：华东交通大学，
类型：发明
国别省市：