一种基于随机凹凸坡面的地震边坡永久位移计算方法技术

本发明专利技术专利公开了一种基于随机凹凸坡面的地震边坡永久位移计算方法，具体涉及地震边坡稳定性技术领域。包括输入相关参数；通过随机生成、获取边坡非规则几何坡面；结合极限分析上限法及强度折减法，通过优化求解安全系数获取边坡在极限状态的临界滑动面参数；假定临界滑动面在地震过程中不再发生变化，通过虚功原理构建边坡的动力平衡方程，获取边坡的角加速度；对角加速度在微小时间段内进行二次积分，通过累积叠加，求得整个地震过程中边坡坡脚处的水平永久位移。采用本发明专利技术技术能快速获取真实水平及竖向地震动作用下考虑坡顶张拉破坏的非规则几何坡面边坡永久位移，解决了地震边坡稳定性分析理论未考虑坡体张拉破坏及实际水平地震动与竖向地震动在地震过程中并非呈确定比例关系的问题，可用于更加准确的计算非规则坡面边坡的地震稳定性。算非规则坡面边坡的地震稳定性。算非规则坡面边坡的地震稳定性。

【技术实现步骤摘要】
一种基于随机凹凸坡面的地震边坡永久位移计算方法


[0001]本专利技术涉及地震边坡稳定性的
，特别涉及一种基于随机凹凸坡面的地震边坡永久位移计算方法。

技术介绍

[0002]近年来，随着川藏铁路工程的逐步推进，在脆弱的生态环境、复杂的地质条件及强烈的活动断裂地区开展基础设施的建设将会不可避免的涉及边坡稳定性的问题。在边坡稳定性分析的理论方法中，多数方法仅考虑了直线型、凸型、凹型坡面形态。然而，在实际工程中，边坡坡面多呈现不规则的凹凸形态，且坡面的凹凸特性对边坡的稳定性具有一定的影响，且坡面的凹凸特性对边坡稳定性的影响规律结论尚不统一。坡面凹凸特性除了关系到边坡的稳定性，更涉及到工程中的放坡开挖方案。除此之外，在以往以永久位移为评价指标的地震稳定性研究中，大多仅考虑了水平地震加速度
‑
时程曲线，或假定竖向地震加速度与水平地震加速度成一定的比例关系。然而，实际水平地震动与竖向地震动在地震过程中并非呈确定比例关系。因此，提出一种可以考虑非规则几何坡面以及真实水平
‑
竖向地震动耦合作用的边坡永久位移计算方法具有理论及实际意义。

技术实现思路

[0003]本专利技术意在提供一种基于随机凹凸坡面的地震边坡永久位移计算方法，解决了边坡稳定性分析理论未考虑实际水平地震动与竖向地震动在地震过程中并非呈确定比例关系的问题。
[0004]为了达到上述目的，本专利技术的技术方案如下：一种基于随机凹凸坡面的地震边坡永久位移计算方法，包括如下步骤：
[0005]S1、输入土体重度γ、粘聚力c、内摩擦角φ、坡顶角α、坡高H和岩土体张拉强度的无量纲参数u；
[0006]S2、设定坡面折线个数n、每段折线高度占坡高的比例为a
i
、每段折线与水平线的夹角为β
i
，通过随机生成a
i
、β
i
获取边坡非规则几何坡面；
[0007]S3、结合极限分析上限法及强度折减法，通过优化求解安全系数获取边坡在极限状态的临界滑动面参数；
[0008]S4、假定临界滑动面在地震过程中不再发生变化，通过虚功原理构建边坡的动力平衡方程，进而获取边坡的角加速度；
[0009]S5、对角加速度在微小时间段内进行二次积分，通过累积叠加，进而求得整个地震过程中边坡坡脚处的水平永久位移。
[0010]进一步的，步骤S2的具体方法：
[0011]S21、设定组成坡面折线个数为n(n＝1,2,3,4
······
),设定a
i
和β
i
满足的约束条件；
[0012]S22、根据S21步骤中的约束条件随机生成坡面折线。
[0013]进一步的，步骤S21的约束条件满足如下公式：
[0014][0015]进一步的，步骤S3的具体方法：
[0016]S31、结合极限分析上限法及强度折减法，获取边坡安全系数的计算公式；
[0017]S32、代入t时刻真实的水平及竖向地震加速度，其中t≥0，计算安全系数；
[0018]S33、当安全系数F
s
≤1.0时，进入步骤S34，否则，返回令t＝t+d
t
，并返回步骤S32；
[0019]S34、假定安全系数在时间步长d
t
内为线性变化，而在时刻t
n
优化获得安全系数F
s(n)
(k
v(n)
，k
h(n)
)≤1.0，并记录前一时间步长t
(n
‑
1)
的边坡安全系数F
s(n
‑
1)
(k
v(n
‑
1)
，k
h(n
‑
1)
)>1.0，通过线性内插法获得F
s
＝1.0时的k
h(s)
和k
v(s)
；
[0020]S35、将k
h(s)
和k
v(s)
代入安全系数计算方程获得临界滑动面参数。
[0021]进一步的，步骤S31中边坡安全系数的计算公式如下：
[0022][0023]上式中，F
s
为安全系数，为内能耗散率系数，k
h
为水平地震加速度系数，k
v
为竖向地震加速度系数，f1～f
n+2
，f1'～f
n+2
'，t1～t3，t1'～t3'均为外功率系数；
[0024]所述内能耗散率系数的计算公式为：
[0025][0026]所述外功率系数f1～f
n+2
的计算公式为：
[0027][0028][0029][0030][0031]所述外功率系数f1'～f
n+2
'的计算公式为：
[0032][0033][0034][0035][0036]所述外功率系数t1～t3的计算公式为：
[0037][0038][0039][0040]所述外功率系数t1'～t3'的计算公式为：
[0041][0042][0043][0044]进一步的，步骤S34中k
h(s)
和k
v(s)
的计算公式为：
[0045][0046]进一步的，步骤S4中边坡的动力平衡方程为：
[0047][0048][0049]上式中，W为滑体的重力大小，l为滑体重心到转动重心P点的距离；
[0050]所述滑体重力的计算公式为：
[0051][0052]所述滑体重心到转动重心P点距离的计算公式为：
[0053][0054]进一步的，步骤S5中水平永久位移的计算公式为：
[0055][0056]与现有技术相比，本方案的有益效果：
[0057]本方案考虑实际工程中边坡的非规则坡面形态和真实水平
‑
竖向地震动耦合效应，使得地震作用下边坡永久位移的计算更加精确，为实际工程中边坡永久位移评价方法奠定了基础。
附图说明
[0058]图1是本专利技术一种基于随机凹凸坡面的地震边坡永久位移计算方法的流程图；
[0059]图2是本专利技术一种基于随机凹凸坡面的地震边坡永久位移计算方法的边坡破坏模式图；
[0060]图3是本专利技术实施例中边坡非规则坡面、临界破坏面图；
[0061]图4是本专利技术实施例中永久位移分析图。
具体实施方式
[0062]下面通过具体实施方式对本专利技术作进一步详细的说明：
[0063]实施例
[0064]如附图1和图2所示：一种基于随机凹凸坡面的地震边坡永久位移计算方法，包括如下步骤：
[0065]S1、输入土体重度γ、粘聚力c、内摩擦角φ、坡顶角α、坡高H和考虑岩土体张拉强度的无量纲参数u；
[0066]S2、设定本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于随机凹凸坡面的地震边坡永久位移计算方法，其特征在于：包括如下步骤：S1、输入土体重度γ、粘聚力c、内摩擦角φ、坡顶角α、坡高H和岩土体张拉强度的无量纲参数u；S2、设定坡面折线个数n、每段折线高度占坡高的比例为a
i
、每段折线与水平线的夹角为β
i
，通过随机生成a
i
、β
i
获取边坡非规则几何坡面；S3、结合极限分析上限法及强度折减法，通过优化求解安全系数获取边坡在极限状态的临界滑动面参数；S4、假定临界滑动面在地震过程中不再发生变化，通过虚功原理构建边坡的动力平衡方程，进而获取边坡的角加速度；S5、对角加速度在微小时间段内进行二次积分，通过累积叠加，进而求得整个地震过程中边坡坡脚处的水平永久位移。2.根据权利要求1所述的一种基于随机凹凸坡面的地震边坡永久位移计算方法，其特征在于：步骤S2的具体方法：S21、设定组成坡面折线个数为n(n＝1,2,3,4
······
),设定a
i
和β
i
满足的约束条件；S22、根据S21步骤中的约束条件随机生成坡面折线。3.根据权利要求2所述的一种基于随机凹凸坡面的地震边坡永久位移计算方法，其特征在于：步骤S21的约束条件满足如下公式：4.根据权利要求1所述的一种基于随机凹凸坡面的地震边坡永久位移计算方法，其特征在于：步骤S3的具体方法：S31、结合极限分析上限法及强度折减法，获取边坡安全系数的计算公式；S32、代入t时刻真实的水平及竖向地震加速度，其中t≥0，计算安全系数；S33、当安全系数F
s
≤1.0时，进入步骤S34，否则，返回令t＝t+d
t
，并返回步骤S32；S34、假定安全系数在时间步长d
t
内为线性变化，而在时刻t
n
优化获得安全系数F
s(n)
(k
v(n)
，k
h(n)
)≤1.0，并记录前一时间步长t
(n
‑
1)
的边坡安全系数...

【专利技术属性】
技术研发人员：程肖，李得建，彭新艳，赵炼恒，张迎宾，余鹏程，贾文韬，
申请(专利权)人：四川大学，
类型：发明
国别省市：