CN111651902B - 一种昔格达地层岩质边坡稳定性计算的拟动力下限法 - Google Patents

一种昔格达地层岩质边坡稳定性计算的拟动力下限法 Download PDF

Info

Publication number
CN111651902B
CN111651902B CN202010611011.9A CN202010611011A CN111651902B CN 111651902 B CN111651902 B CN 111651902B CN 202010611011 A CN202010611011 A CN 202010611011A CN 111651902 B CN111651902 B CN 111651902B
Authority
CN
China
Prior art keywords
xigeda
rock slope
block unit
slope
stratum
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.)
Active
Application number
CN202010611011.9A
Other languages
English (en)
Other versions
CN111651902A (zh
Inventor
刘文连
李泽
眭素刚
唐果
陆得志
宁飞
安家金
樊亚红
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Kunming University of Science and Technology
Kunming Prospecting Design Institute of China Nonferrous Metals Industry Co Ltd
Original Assignee
Kunming University of Science and Technology
Kunming Prospecting Design Institute of China Nonferrous Metals Industry Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Kunming University of Science and Technology, Kunming Prospecting Design Institute of China Nonferrous Metals Industry Co Ltd filed Critical Kunming University of Science and Technology
Priority to CN202010611011.9A priority Critical patent/CN111651902B/zh
Publication of CN111651902A publication Critical patent/CN111651902A/zh
Application granted granted Critical
Publication of CN111651902B publication Critical patent/CN111651902B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/12—Timing analysis or timing optimisation
    • 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/14—Force analysis or force optimisation, e.g. static or dynamic forces

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)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种昔格达地层岩质边坡稳定性计算的拟动力下限法，本发明以昔格达地层岩质边坡为研究对象，采用块体单元离散昔格达地层岩质边坡，将地震波假设为简谐波，并使用拟动力法原理计算块体单元的地震惯性力，根据下限法原理建立地震作用下昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型；使用“序列二次规划算法”循环求解昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型，获得地震作用下岩质边坡稳定性的安全系数与时间的关系曲线。本发明方法具有概念明确、计算效率和计算精度高等特点，可将其应用于地震作用下昔格达地层岩质边坡的稳定性分析，可为岩质边坡的抗震设计提供科学依据。

Description

一种昔格达地层岩质边坡稳定性计算的拟动力下限法
技术领域
本发明涉及一种昔格达地层岩质边坡稳定性计算的拟动力下限法，属于边坡稳定性分析的技术领域。
背景技术
在四川云南交界地区广泛分布着昔格达地层。昔格达地层形成与晚更新世与上第三系之间，属于河湖相沉积地层，其最典型的特征是粉砂岩和粘土岩互层交替出现。昔格达地层岩体是半层岩，属于极软岩，其力学特点是：抗压强度低、容易软化、遇水易崩解。
在四川省、云南省交界地区，地震频发、地震烈度较高，在工程建设中，经常碰到昔格达地层岩质边坡在地震作用下发生失稳破坏的情况。在当前对于昔格达地层岩质边坡地震作用下的稳定计算主要采用拟静力法，拟静力法将地震惯性力简化成静力荷载，并假设地震惯性力是随边坡高度变化的函数，求解地震拟静力载荷以后基于强度折减思想使用刚体极限平衡法求解边坡的动力稳定安全系数。但拟静力法存在一些不足，比如其不仅没有考虑地震的振幅、频率以及地震持续时间等，而且还忽略了边坡岩体材料的阻尼特性，因此拟静力法的计算结果的精度并不高；同时，拟静力法不能获得安全系数与地震时间的关系，因此该方法并不能满足工程的实际需求。
发明内容
本发明提供了一种昔格达地层岩质边坡稳定性计算的拟动力下限法，以用于通过下限法获得昔格达地层岩质边坡稳定性的安全系数分布规律。
本发明的技术方案是：一种昔格达地层岩质边坡稳定性计算的拟动力下限法，以昔格达地层岩质边坡为研究对象，采用块体单元离散昔格达地层岩质边坡，将地震波假设为简谐波正弦波，并使用拟动力法原理计算块体单元的地震惯性力，进行地震作用下块体的受力分析；根据下限法原理建立目标函数、建立地震作用下块体单元的动力平衡方程和结构面的屈服条件，进而建立地震作用下昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型；求解地震作用下昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型。
所述方法具体步骤包括：步骤一、拟定昔格达地层岩质边坡稳定性计算的参数；步骤二、采用块体单元离散昔格达地层岩质边坡；步骤三、进行地震作用下昔格达地层岩质边坡块体单元的受力分析；步骤四、根据拟动力法原理计算受力分析中的昔格达地层岩质边坡块体单元的地震惯性力；步骤五、结合地震惯性力，根据下限法原理建立地震作用下昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型；步骤六、求解地震作用下昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型，获得地震作用下岩质边坡稳定性的安全系数与时间的关系曲线。
所述拟定昔格达地层岩质边坡稳定性计算的参数包括：①确定昔格达地层岩质边坡的几何参数；②确定昔格达地层中粉砂岩地层、粘土岩地层的分布情况；③确定昔格达地层岩体材料的物理力学参数；④确定昔格达地层岩质边坡的地震参数。
所述昔格达地层岩质边坡的几何参数包括：岩质边坡的高度H，岩质边坡的宽度，边坡几何形状控制点的坐标；地层的分布情况包括：粉砂岩和粘土岩地层每一层的厚度，地层的倾角；物理力学参数包括：粉砂岩和粘土岩容重，粉砂岩的内摩擦角和凝聚力，粘土岩的内摩擦角和凝聚力；地震参数包括：地震振动周期T、地震放大系数fs、边坡岩体的剪切波波速Vs、边坡岩体的纵波波速Vp、水平方向的地震加速度系数kh、竖直方向的地震加速度系数kv。
所述采用块体单元离散昔格达地层岩质边坡具体为：每个块体单元为外凸多边形，每条边为一条结构面，结构面上作用有法向力和剪力；块体单元为刚性体，不会发生任何变形和破坏，边坡只会沿着相邻块体单元之间的结构面破坏，假设块体单元的结构面仅发生剪切破坏；粉砂岩和粘土岩按岩层划分块体单元，粉砂岩和粘土岩的交界面为相邻单元的结构面。
所述受力分析具体为：
块体单元i有ni条结构面，块体单元的形心为Oi，块体单元的受力如下：
(1)建立昔格达地层岩质边坡的坐标系，以边坡的坡脚为坐标原点，水平轴为坐标系的x轴，x轴以水平向右为正，竖直轴为坐标系的y轴，y轴以竖直向上为正；昔格达地层岩质边坡第i个块体单元的第k条结构面的长度为
Figure GDA0002586913620000021
第i个块体单元中第k条结构面的倾角为
Figure GDA0002586913620000022
第i个块体单元的第k条结构面上作用有两个力：法向力
Figure GDA0002586913620000023
以受拉为正；剪力
Figure GDA0002586913620000024
以对块体单元产生逆时针的转动效果为正；其中，k＝(1,…,ni)，ni是第i个块体单元中结构面的数量；
(2)第i个块体单元的形心Oi处作用有三个力：块体单元的自重Gi；第j个时步昔格达地层岩质边坡第i个块体单元沿水平方向的地震惯性力
Figure GDA0002586913620000031
第j个时步昔格达地层岩质边坡第i个块体单元沿竖直方向的地震惯性力
Figure GDA0002586913620000032
所述地震惯性力计算具体为：
(1)根据拟动力法原理计算块体单元水平方向的地震惯性力，具体如下：
Figure GDA0002586913620000033
式中：
Figure GDA0002586913620000034
是第j个时步昔格达地层岩质边坡第i个块体单元沿水平方向的地震惯性力；i＝(1,…,Nb)，Nb是块体单元的数量；j＝(1,…,Nt)，Nt是地震振动周期内的时步数量；Ai是第i个块体单元的面积；γ是容重：位于粉砂岩地层时γ取粉砂岩的容重、位于粘土岩地层时γ取粘土岩的容重；kh是水平方向的地震加速度系数；g是重力加速度；Hi是第i个块体单元的形心至边坡的坡脚处的垂直高度；H是昔格达地层岩质边坡的高度；fs为边坡岩体的地震放大系数；π为圆周率；T为地震振动周期；Vs是边坡岩体的剪切波波速；
(2)根据拟动力法原理计算块体单元竖直方向的地震惯性力，具体如下：
Figure GDA0002586913620000035
式中：
Figure GDA0002586913620000036
是第j个时步昔格达地层岩质边坡第i个块体单元沿竖直方向的地震惯性力；kv是竖直方向的地震加速度系数；Vp是边坡岩体的纵波波速。
所述建立地震作用下昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型具体为：
(1)建立目标函数
设昔格达地层岩质边坡的安全系数为目标函数，并求解边坡安全系数的最大值，目标函数如下：
Maximize:Kj
式中：Kj是第j个时步昔格达地层岩质边坡动力稳定性的安全系数，j＝(1,…,Nt)，Nt是地震振动周期内的时步数量；Maximize表示“使最大”；
(2)建立地震作用下块体单元的动力平衡方程
①建立地震作用下第j个时步昔格达地层岩质边坡第i个块体单元的水平方向的动力平衡方程，具体如下：
Figure GDA0002586913620000041
式中：
Figure GDA0002586913620000042
是昔格达地层岩质边坡第i个块体单元中第k条结构面的法向力，
Figure GDA0002586913620000043
以受拉为正；
Figure GDA0002586913620000044
第i个块体单元中第k条结构面的倾角；i＝(1,…,Nb)，Nb是块体单元的数量；j＝(1,…,Nt)，Nt是地震振动周期内的时步数量；k＝(1,…,ni)，ni是第i个块体单元中结构面的数量；
Figure GDA0002586913620000045
是昔格达地层岩质边坡第i个块体单元中第k条结构面的剪力，
Figure GDA0002586913620000046
以对块体单元产生逆时针的转动效果为正；
Figure GDA0002586913620000047
是第j个时步昔格达地层岩质边坡第i个块体单元沿水平方向的地震惯性力；
②建立地震作用下第j个时步昔格达地层岩质边坡第i个块体单元的竖直方向的动力平衡方程，具体如下：
Figure GDA0002586913620000048
式中：
Figure GDA0002586913620000049
是第j个时步昔格达地层岩质边坡第i个块体单元沿竖直方向的地震惯性力；Gi是昔格达地层岩质边坡第i个块体单元的自重，Gi＝Aiγ，Ai是第i个块体单元的面积，γ是容重：位于粉砂岩地层时γ取粉砂岩的容重、位于粘土岩地层时γ取粘土岩的容重；
(3)建立块体单元结构面的屈服条件，具体如下：
Figure GDA00025869136200000410
式中：
Figure GDA00025869136200000411
是内摩擦角：位于粉砂岩地层时
Figure GDA00025869136200000412
取粉砂岩的内摩擦角、位于粘土岩地层时
Figure GDA00025869136200000413
取粘土岩的内摩擦角；
Figure GDA00025869136200000414
是昔格达地层岩质边坡第i个块体单元中第k条结构面的长度；cs是凝聚力：位于粉砂岩地层时cs取粉砂岩的凝聚力、位于粘土岩地层时cs取粘土岩的凝聚力；
(4)建立地震作用下昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型
将目标函数、块体单元的动力平衡方程、结构面的屈服条件、地震惯性力方程集成，可得到地震作用下昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型如下：
Figure GDA0002586913620000051
式中：
Figure GDA0002586913620000052
是第j个时步昔格达地层岩质边坡第i个块体单元沿水平方向的地震惯性力；kh是水平方向的地震加速度系数；g是重力加速度；Hi是第i个块体单元的形心至边坡的坡脚处的垂直高度；H是昔格达地层岩质边坡的高度；fs为边坡岩体的地震放大系数；π为圆周率；T为地震振动周期；Vs是边坡岩体的剪切波波速；
Figure GDA0002586913620000053
是第j个时步昔格达地层岩质边坡第i个块体单元沿竖直方向的地震惯性力；kv是竖直方向的地震加速度系数；Vp是边坡岩体的纵波波速。
所述求解地震作用下昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型，具体为：将已知参数从j＝1到j＝Nt循环带入地震作用下昔格达地层岩质边坡稳定性的拟动力法非线性数学规划模型式，并使用“序列二次规划法”求解地震作用下昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型，求解得到Nt个昔格达地层岩质边坡动力稳定性的安全系数Kj；然后以安全系数为竖直轴、以时间为横轴绘制昔格达地层岩质边坡动力稳定性的安全系数Kj与时间jT/Nt的关系曲线；其中，j＝(1,…,Nt)，Nt是地震振动周期内的时步数量；T为地震振动周期。
本发明的有益效果是：本发明以昔格达地层岩质边坡为研究对象，采用块体单元离散昔格达地层岩质边坡，将地震波假设为简谐波，并使用拟动力法原理计算块体单元的地震惯性力，根据下限法原理建立地震作用下昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型；使用“序列二次规划算法”循环求解昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型，获得地震作用下岩质边坡稳定性的安全系数与时间的关系曲线。本发明方法具有概念明确、计算效率和计算精度高等特点，可将其应用于地震作用下昔格达地层岩质边坡的稳定性分析，可为岩质边坡的抗震设计提供科学依据。
附图说明
图1为本发明技术流程图；
图2为昔格达地层岩质边坡的块体单元受力示意图；
图3为实施例昔格达地层岩质边坡的几何形状示意图(单位：m)；
图4为实施例昔格达地层岩质边坡块体单元离散示意图和单元编号；
图5为实施例昔格达地层岩质边坡的第1个块体单元沿水平方向的地震惯性力与时间的关系曲线；
图6为实施例昔格达地层岩质边坡的第1个块体单元沿竖直方向的地震惯性力与时间的关系曲线；
图7为实施例昔格达地层岩质边坡的第9个块体单元沿水平方向的地震惯性力与时间的关系曲线；
图8为实施例昔格达地层岩质边坡的第9个块体单元沿竖直方向的地震惯性力与时间的关系曲线；
图9为实施例昔格达地层岩质边坡的第20个块体单元沿水平方向的地震惯性力与时间的关系曲线；
图10为实施例昔格达地层岩质边坡的第20个块体单元沿竖直方向的地震惯性力与时间的关系曲线；
图11实施例昔格达地层岩质边坡在地震作用下的安全系数与时间的关系曲线。
具体实施方式
下面结合附图和实施例，对本发明作进一步说明，但本发明的内容并不限于所述范围。
实施例1：如图1-11所示，一种昔格达地层岩质边坡稳定性计算的拟动力下限法，以昔格达地层岩质边坡为研究对象，采用块体单元离散昔格达地层岩质边坡，将地震波假设为简谐波正弦波，并使用拟动力法原理计算块体单元的地震惯性力，进行地震作用下块体的受力分析；根据下限法原理建立目标函数、建立地震作用下块体单元的动力平衡方程和结构面的屈服条件，进而建立地震作用下昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型；求解地震作用下昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型。本发明流程如图1所示。
进一步地，可以设置所述方法具体步骤包括：步骤一、拟定昔格达地层岩质边坡稳定性计算的参数；步骤二、采用块体单元离散昔格达地层岩质边坡；步骤三、进行地震作用下昔格达地层岩质边坡块体单元的受力分析；步骤四、根据拟动力法原理计算受力分析中的昔格达地层岩质边坡块体单元的地震惯性力；步骤五、结合地震惯性力，根据下限法原理建立地震作用下昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型；步骤六、求解地震作用下昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型，获得地震作用下岩质边坡稳定性的安全系数与时间的关系曲线。
更进一步地，本发明结合具体步骤，给出如下过程：
步骤一、拟定昔格达地层岩质边坡稳定性计算的参数。
根据昔格达地层岩质边坡的实际情况，拟定稳定性计算参数，具体为：
①确定昔格达地层岩质边坡的几何参数，实施例边坡的几何形状如图3所示，包括：边坡的高度H＝60m，边坡的宽度为100m；边坡几何形状控制点的坐标分别为：O点的坐标(0,0)、A点的坐标(100,0)、B点的坐标(100,60)，C点的坐标(27.87,60)；
②确定昔格达地层中粉砂岩地层、粘土岩地层的分布情况，包括：粉砂岩和粘土岩地层的分布情况如图3所示，粉砂岩和粘土岩地层的厚度均为10m，粉砂岩和粘土岩的交界面的倾角为30°；
③确定昔格达地层岩体材料的物理力学参数，包括：粉砂岩和粘土岩容重均取2500kN/m3，粉砂岩的内摩擦角取26°，粉砂岩的凝聚力100kPa，粘土岩的内摩擦角取24°，粘土岩凝聚力160kPa。
④确定昔格达地层岩质边坡的地震参数，包括：地震振动周期T＝0.2s；地震放大系数fs＝1.1；边坡岩体的剪切波波速Vs＝2747m/s；边坡岩体的纵波波速Vp＝5139m/s；水平方向的地震加速度系数kh＝0.1，竖直方向的地震加速度系数kv＝0.05。
步骤二、采用块体单元离散昔格达地层岩质边坡。
具体为：使用块体单元离散昔格达地层岩质边坡，块体单元如图2所示，每个块体单元为外凸多边形，每条边为一条结构面，结构面上作用有法向力和剪力；块体单元为刚性体，不会发生任何变形和破坏，边坡只会沿着相邻块体单元之间的结构面破坏，假设块体单元的结构面仅发生剪切破坏；粉砂岩和粘土岩按岩层划分块体单元，粉砂岩和粘土岩的交界面为相邻单元的结构面。实施例昔格达地层岩质边坡块体单元离散的示意和单元编号如图4所示，共划分块体单元86个、结构面183条。
步骤三、进行地震作用下昔格达地层岩质边坡块体单元的受力分析。
昔格达地层岩质边坡的第i个块体单元的受力情况如图2所示，块体单元i有ni条结构面，块体单元的形心为Oi，块体单元的受力如下：
(1)昔格达地层岩质边坡第i个块体单元的第k条结构面的长度为
Figure GDA0002586913620000081
第i个块体单元中第k条结构面的倾角为
Figure GDA0002586913620000082
第i个块体单元的第k条结构面上作用有两个力：法向力
Figure GDA0002586913620000083
以受拉为正；剪力
Figure GDA0002586913620000084
以对块体单元产生逆时针的转动效果为正；其中k＝(1,…,ni)，ni是第i个块体单元中结构面的数量；
(2)第i个块体单元的形心Oi处作用有三个力：块体单元的自重Gi；第j个时步昔格达地层岩质边坡第i个块体单元沿水平方向的地震惯性力
Figure GDA0002586913620000085
第j个时步昔格达地层岩质边坡第i个块体单元沿竖直方向的地震惯性力
Figure GDA0002586913620000086
步骤四、根据拟动力法原理计算昔格达地层岩质边坡块体单元的地震惯性力。
本发明将地震波假设为简谐波，采用拟动力法原理计算昔格达地层岩质边坡块体单元的地震惯性力，具体如下：
(1)根据拟动力法原理计算块体单元水平方向的地震惯性力。
Figure GDA0002586913620000087
式中：i＝(1,…,Nb)，Nb是块体单元的数量，Nb＝86；j＝(1,…,Nt)，Nt是地震振动周期内的时步数量，Nt＝20；
Figure GDA0002586913620000091
是第j个时步昔格达地层岩质边坡第i个块体单元沿水平方向的地震惯性力；Hi是第i个块体单元的形心至边坡的坡脚处的垂直高度；H是昔格达地层岩质边坡的高度，H＝60m；fs为边坡岩体的地震放大系数；Ai是第i个块体单元的面积；γ是块体单元的容重：当块体单元位于粉砂岩地层时γ取粉砂岩的容重、当块体单元位于粘土岩地层时γ取粘土岩的容重；kh是水平方向的地震加速度系数，kh取0.1；g是重力加速度，g取9.81m/s2；T为地震的振动周期，T＝0.2s；地震振动周期内的时步数量Nt取20，Δt＝T/Nt是时步的时间长度；π为圆周率，π取3.14；Vs是边坡岩体的剪切波波速，Vs＝2747m/s。
根据上式计算了实施例昔格达地层岩质边坡所有块体单元的沿水平方向的地震惯性力，第1个、第9个、第20个块体单元的沿水平方向的地震惯性力分别如图5、如图7、图9所示。
(2)根据拟动力法原理计算块体单元竖直方向的地震惯性力。
Figure GDA0002586913620000092
上式中：i＝(1,…,Nb)，Nb是块体单元的数量，Nb＝86；j＝(1,…,Nt)，Nt是地震振动周期内的时步数量，Nt＝20；
Figure GDA0002586913620000093
是第j个时步昔格达地层岩质边坡第i个块体单元沿竖直方向的地震惯性力；Hi是第i个块体单元的形心至边坡的坡脚处的垂直高度；H是昔格达地层岩质边坡的高度，H＝60m；Δt是时步的时间长度，Δt＝0.01s；fs为边坡岩体的地震放大系数，fs＝1.1；Ai是第i个块体单元的面积；γ是块体单元的容重：当块体单元位于粉砂岩地层时γ取粉砂岩的容重、当块体单元位于粘土岩地层时γ取粘土岩的容重；kv是竖直方向的地震加速度系数，kh取0.01；g是重力加速度，g取9.81m/s2；T为地震的振动周期，T＝0.2s；π为圆周率，π取3.14；Vp是边坡岩体的纵波波速，Vp＝5139m/s。
根据上式计算了实施例昔格达地层岩质边坡所有块体单元的沿竖直方向的地震惯性力，第1个、第9个、第20个块体单元的沿竖直方向的地震惯性力分别如图6、如图8、图10所示。
步骤五、根据下限法原理建立地震作用下昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型。具体如下：
(1)建立目标函数。
本发明设昔格达地层岩质边坡的安全系数为目标函数，并求解边坡安全系数的最大值。目标函数如下：
Maximize:Kj(3)
式中：Kj是第j个时步昔格达地层岩质边坡动力稳定性的安全系数，j＝(1,…,Nt)，Nt是地震振动周期内的时步数量；Maximize表示“使最大”。
(2)建立地震作用下块体单元的动力平衡方程。
块体单元受自重、地震惯性力、结构面的法向力和剪力作用而保持平衡。地震作用下块体单元的动力平衡方程按以下流程建立：
①建立地震作用下第j个时步昔格达地层岩质边坡第i个块体单元的水平方向的动力平衡方程。具体如下：
Figure GDA0002586913620000101
式中：i＝(1,…,Nb)，Nb是块体单元的数量，Nb＝86；j＝(1,…,Nt)，Nt是地震振动周期内的时步数量，Nt＝20；
Figure GDA0002586913620000102
第i个块体单元中第k条结构面的倾角；ni是第i个块体单元中结构面的数量；
Figure GDA0002586913620000103
是昔格达地层岩质边坡第i个块体单元中第k条结构面的法向力，
Figure GDA0002586913620000104
以受拉为正；
Figure GDA0002586913620000105
是昔格达地层岩质边坡第i个块体单元中第k条结构面的剪力，
Figure GDA0002586913620000106
以对块体单元产生逆时针的转动效果为正；
Figure GDA0002586913620000107
是第j个时步昔格达地层岩质边坡第i个块体单元沿水平方向的地震惯性力。
②建立地震作用下第j个时步昔格达地层岩质边坡第i个块体单元的竖直方向的动力平衡方程。具体如下：
Figure GDA0002586913620000108
式中：i＝(1,…,Nb)，Nb是块体单元的数量，Nb＝86；j＝(1,…,Nt)，Nt是地震振动周期内的时步数量，Nt＝20；
Figure GDA0002586913620000111
第i个块体单元中第k条结构面的倾角；ni是第i个块体单元中结构面的数量；
Figure GDA0002586913620000112
是昔格达地层岩质边坡第i个块体单元中第k条结构面的法向力，
Figure GDA0002586913620000113
以受拉为正；
Figure GDA0002586913620000114
是昔格达地层岩质边坡第i个块体单元中第k条结构面的剪力，
Figure GDA0002586913620000115
以对块体单元产生逆时针的转动效果为正；Gi是昔格达地层岩质边坡第i个块体单元的自重,Gi＝Aiγ，Ai是第i个块体单元的面积，γ是块体单元的容重，当块体单元位于粉砂岩地层时γ取粉砂岩的容重、当块体单元位于粘土岩地层时γ取粘土岩的容重；
Figure GDA0002586913620000116
是第j个时步昔格达地层岩质边坡第i个块体单元沿竖直方向的地震惯性力。
(3)建立块体单元结构面的屈服条件。具体如下：
Figure GDA0002586913620000117
式中：i＝(1,…,Nb)，Nb是块体单元的数量，Nb＝86；j＝(1,…,Nt)，Nt是地震振动周期内的时步数量，Nt＝20；k＝(1,…,ni)，ni是第i个块体单元中结构面的数量；
Figure GDA0002586913620000118
是块体单元的内摩擦角：粉砂岩的内摩擦角取26°、粘土岩的内摩擦角取24°，cs是块体单元的凝聚力：粉砂岩的凝聚力100kPa、粉砂岩的凝聚力100kPa；
Figure GDA0002586913620000119
是昔格达地层岩质边坡第i个块体单元中第k条结构面的法向力，
Figure GDA00025869136200001110
以受拉为正；
Figure GDA00025869136200001111
是昔格达地层岩质边坡第i个块体单元中第k条结构面的剪力，
Figure GDA00025869136200001112
以对块体单元产生逆时针的转动效果为正；Kj是第j个时步昔格达地层岩质边坡动力稳定性的安全系数。
(4)建立地震作用下昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型。将目标函数、块体单元的动力平衡方程、结构面的屈服条件、地震惯性力集成，可得到地震作用下昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型如下：
Figure GDA0002586913620000121
上式中：i＝(1,…,Nb)，Nb是块体单元的数量；j＝(1,…,Nt)，Nt是地震振动周期内的时步数量。
步骤六、求解地震作用下昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型。
具体如为：将已知参数从j＝1到j＝Nt循环带入地震作用下昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型式(7)，并使用“序列二次规划法”求解地震作用下昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型，求解得到Nt个昔格达地层岩质边坡动力稳定性的安全系数Kj如表1所示，并绘制昔格达地层岩质边坡动力稳定性的安全系数Kj与时间的关系曲线如图11所示。
表1昔格达地层岩质边坡动力稳定性的安全系数
Figure GDA0002586913620000122
上面结合附图对本发明的具体实施方式作了详细说明，但是本发明并不限于上述实施方式，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下作出各种变化。

Claims (8)

1.一种昔格达地层岩质边坡稳定性计算的拟动力下限法，其特征在于：以昔格达地层岩质边坡为研究对象，采用块体单元离散昔格达地层岩质边坡，将地震波假设为简谐波正弦波，并使用拟动力法原理计算块体单元的地震惯性力，进行地震作用下块体的受力分析；根据下限法原理建立目标函数、建立地震作用下块体单元的动力平衡方程和结构面的屈服条件，进而建立地震作用下昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型；求解地震作用下昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型；
所述建立地震作用下昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型具体为：
(1)建立目标函数
设昔格达地层岩质边坡的安全系数为目标函数，并求解边坡安全系数的最大值，目标函数如下：
Maximize:Kj
式中：Kj是第j个时步昔格达地层岩质边坡动力稳定性的安全系数，j＝(1,…,Nt)，Nt是地震振动周期内的时步数量；Maximize表示“使最大”；
(2)建立地震作用下块体单元的动力平衡方程
①建立地震作用下第j个时步昔格达地层岩质边坡第i个块体单元的水平方向的动力平衡方程，具体如下：
Figure FDA0004005230290000011
式中：
Figure FDA0004005230290000012
是昔格达地层岩质边坡第i个块体单元中第k条结构面的法向力，
Figure FDA0004005230290000013
以受拉为正；
Figure FDA0004005230290000014
第i个块体单元中第k条结构面的倾角；i＝(1,…,Nb)，Nb是块体单元的数量；j＝(1,…,Nt)，Nt是地震振动周期内的时步数量；k＝(1,…,ni)，ni是第i个块体单元中结构面的数量；
Figure FDA0004005230290000015
是昔格达地层岩质边坡第i个块体单元中第k条结构面的剪力，
Figure FDA0004005230290000016
以对块体单元产生逆时针的转动效果为正；
Figure FDA0004005230290000017
是第j个时步昔格达地层岩质边坡第i个块体单元沿水平方向的地震惯性力；
②建立地震作用下第j个时步昔格达地层岩质边坡第i个块体单元的竖直方向的动力平衡方程，具体如下：
Figure FDA0004005230290000021
式中：
Figure FDA0004005230290000022
是第j个时步昔格达地层岩质边坡第i个块体单元沿竖直方向的地震惯性力；Gi是昔格达地层岩质边坡第i个块体单元的自重，Gi＝Aiγ，Ai是第i个块体单元的面积，γ是容重：位于粉砂岩地层时γ取粉砂岩的容重、位于粘土岩地层时γ取粘土岩的容重；
(3)建立块体单元结构面的屈服条件，具体如下：
Figure FDA0004005230290000023
式中：
Figure FDA0004005230290000024
是内摩擦角：位于粉砂岩地层时
Figure FDA0004005230290000025
取粉砂岩的内摩擦角、位于粘土岩地层时
Figure FDA0004005230290000026
取粘土岩的内摩擦角；
Figure FDA0004005230290000027
是昔格达地层岩质边坡第i个块体单元中第k条结构面的长度；cs是凝聚力：位于粉砂岩地层时cs取粉砂岩的凝聚力、位于粘土岩地层时cs取粘土岩的凝聚力；
(4)建立地震作用下昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型
将目标函数、块体单元的动力平衡方程、结构面的屈服条件、地震惯性力方程集成，可得到地震作用下昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型如下：
Figure FDA0004005230290000028
式中：
Figure FDA0004005230290000031
是第j个时步昔格达地层岩质边坡第i个块体单元沿水平方向的地震惯性力；kh是水平方向的地震加速度系数；g是重力加速度；Hi是第i个块体单元的形心至边坡的坡脚处的垂直高度；H是昔格达地层岩质边坡的高度；fs为边坡岩体的地震放大系数；π为圆周率；T为地震振动周期；Vs是边坡岩体的剪切波波速；
Figure FDA0004005230290000032
是第j个时步昔格达地层岩质边坡第i个块体单元沿竖直方向的地震惯性力；kv是竖直方向的地震加速度系数；Vp是边坡岩体的纵波波速。
2.根据权利要求1所述的昔格达地层岩质边坡稳定性计算的拟动力下限法，其特征在于：所述方法具体步骤包括：
步骤一、拟定昔格达地层岩质边坡稳定性计算的参数；
步骤二、采用块体单元离散昔格达地层岩质边坡；
步骤三、进行地震作用下昔格达地层岩质边坡块体单元的受力分析；
步骤四、根据拟动力法原理计算受力分析中的昔格达地层岩质边坡块体单元的地震惯性力；
步骤五、结合地震惯性力，根据下限法原理建立地震作用下昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型；
步骤六、求解地震作用下昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型，获得地震作用下岩质边坡稳定性的安全系数与时间的关系曲线。
3.根据权利要求2所述的昔格达地层岩质边坡稳定性计算的拟动力下限法，其特征在于：所述拟定昔格达地层岩质边坡稳定性计算的参数包括：①确定昔格达地层岩质边坡的几何参数；②确定昔格达地层中粉砂岩地层、粘土岩地层的分布情况；③确定昔格达地层岩体材料的物理力学参数；④确定昔格达地层岩质边坡的地震参数。
4.根据权利要求3所述的昔格达地层岩质边坡稳定性计算的拟动力下限法，其特征在于：所述昔格达地层岩质边坡的几何参数包括：岩质边坡的高度H，岩质边坡的宽度，边坡几何形状控制点的坐标；地层的分布情况包括：粉砂岩和粘土岩地层每一层的厚度，地层的倾角；物理力学参数包括：粉砂岩和粘土岩容重，粉砂岩的内摩擦角和凝聚力，粘土岩的内摩擦角和凝聚力；地震参数包括：地震振动周期T、地震放大系数fs、边坡岩体的剪切波波速Vs、边坡岩体的纵波波速Vp、水平方向的地震加速度系数kh、竖直方向的地震加速度系数kv。
5.根据权利要求1或2所述的昔格达地层岩质边坡稳定性计算的拟动力下限法，其特征在于：所述采用块体单元离散昔格达地层岩质边坡具体为：每个块体单元为外凸多边形，每条边为一条结构面，结构面上作用有法向力和剪力；块体单元为刚性体，不会发生任何变形和破坏，边坡只会沿着相邻块体单元之间的结构面破坏，假设块体单元的结构面仅发生剪切破坏；粉砂岩和粘土岩按岩层划分块体单元，粉砂岩和粘土岩的交界面为相邻单元的结构面。
6.根据权利要求1或2所述的昔格达地层岩质边坡稳定性计算的拟动力下限法，其特征在于：所述受力分析具体为：
块体单元i有ni条结构面，块体单元的形心为Oi，块体单元的受力如下：
(1)建立昔格达地层岩质边坡的坐标系，以边坡的坡脚为坐标原点，水平轴为坐标系的x轴，x轴以水平向右为正，竖直轴为坐标系的y轴，y轴以竖直向上为正；昔格达地层岩质边坡第i个块体单元的第k条结构面的长度为
Figure FDA0004005230290000041
第i个块体单元中第k条结构面的倾角为
Figure FDA0004005230290000042
第i个块体单元的第k条结构面上作用有两个力：法向力
Figure FDA0004005230290000043
Figure FDA0004005230290000044
以受拉为正；剪力
Figure FDA0004005230290000045
Figure FDA0004005230290000046
以对块体单元产生逆时针的转动效果为正；其中，k＝(1,…,ni)，ni是第i个块体单元中结构面的数量；
(2)第i个块体单元的形心Oi处作用有三个力：块体单元的自重Gi；第j个时步昔格达地层岩质边坡第i个块体单元沿水平方向的地震惯性力
Figure FDA0004005230290000047
第j个时步昔格达地层岩质边坡第i个块体单元沿竖直方向的地震惯性力
Figure FDA0004005230290000048
7.根据权利要求1或2所述的昔格达地层岩质边坡稳定性计算的拟动力下限法，其特征在于：所述地震惯性力计算具体为：
(1)根据拟动力法原理计算块体单元水平方向的地震惯性力，具体如下：
Figure FDA0004005230290000049
式中：
Figure FDA00040052302900000410
是第j个时步昔格达地层岩质边坡第i个块体单元沿水平方向的地震惯性力；i＝(1,…,Nb)，Nb是块体单元的数量；j＝(1,…,Nt)，Nt是地震振动周期内的时步数量；Ai是第i个块体单元的面积；γ是容重：位于粉砂岩地层时γ取粉砂岩的容重、位于粘土岩地层时γ取粘土岩的容重；kh是水平方向的地震加速度系数；g是重力加速度；Hi是第i个块体单元的形心至边坡的坡脚处的垂直高度；H是昔格达地层岩质边坡的高度；fs为边坡岩体的地震放大系数；π为圆周率；T为地震振动周期；Vs是边坡岩体的剪切波波速；
(2)根据拟动力法原理计算块体单元竖直方向的地震惯性力，具体如下：
Figure FDA0004005230290000051
式中：
Figure FDA0004005230290000052
是第j个时步昔格达地层岩质边坡第i个块体单元沿竖直方向的地震惯性力；kv是竖直方向的地震加速度系数；Vp是边坡岩体的纵波波速。
8.根据权利要求1或2所述的昔格达地层岩质边坡稳定性计算的拟动力下限法，其特征在于：所述求解地震作用下昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型，具体为：将已知参数从j＝1到j＝Nt循环带入地震作用下昔格达地层岩质边坡稳定性的拟动力法非线性数学规划模型式，并使用“序列二次规划法”求解地震作用下昔格达地层岩质边坡稳定性的拟动力下限法非线性数学规划模型，求解得到Nt个昔格达地层岩质边坡动力稳定性的安全系数Kj；然后以安全系数为竖直轴、以时间为横轴绘制昔格达地层岩质边坡动力稳定性的安全系数Kj与时间jT/Nt的关系曲线；其中，j＝(1,…,Nt)，Nt是地震振动周期内的时步数量；T为地震振动周期。
CN202010611011.9A 2020-06-29 2020-06-29 一种昔格达地层岩质边坡稳定性计算的拟动力下限法 Active CN111651902B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010611011.9A CN111651902B (zh) 2020-06-29 2020-06-29 一种昔格达地层岩质边坡稳定性计算的拟动力下限法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010611011.9A CN111651902B (zh) 2020-06-29 2020-06-29 一种昔格达地层岩质边坡稳定性计算的拟动力下限法

Publications (2)

Publication Number Publication Date
CN111651902A CN111651902A (zh) 2020-09-11
CN111651902B true CN111651902B (zh) 2023-03-31

Family

ID=72344352

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010611011.9A Active CN111651902B (zh) 2020-06-29 2020-06-29 一种昔格达地层岩质边坡稳定性计算的拟动力下限法

Country Status (1)

Country Link
CN (1) CN111651902B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113235671B (zh) * 2021-05-07 2022-12-02 中国建筑西南勘察设计研究院有限公司 一种岩质边坡地基的承载力载荷试验系统及方法
CN113701976B (zh) * 2021-08-17 2022-10-25 昆明理工大学 一种振动方向可调的土木工程结构抗震试验装置

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103135128B (zh) * 2013-01-31 2015-07-01 重庆大学 一种地震荷载作用下三维边坡稳定性预测方法
US20180292299A1 (en) * 2014-01-13 2018-10-11 Hubei University Of Technology Method of critical displacement forecast based on the deformation failure mechanism of slope
CN106126796B (zh) * 2016-06-20 2019-01-04 昆明理工大学 一种带岩桥的岩质边坡极限承载力的塑性极限分析下限法
DE102016122467A1 (de) * 2016-11-22 2018-05-24 Geo Explorers AG Verfahren und Vorrichtung zur Einschätzung des Risikos eines natürlich verursachten Fels- oder/und Bergsturzes an einem Hang, insbesondere einem Gebirgshang
CN106709182B (zh) * 2016-12-26 2020-02-14 华中科技大学 一种地震作用下顺层岩质边坡稳定可靠性安全评价方法
CN107341285B (zh) * 2017-05-26 2020-08-25 昆明理工大学 一种地震力作用下节理岩质边坡倾倒破坏极限承载力分析下限法
CN107330145B (zh) * 2017-05-26 2020-07-31 昆明理工大学 一种考虑转动效应的岩质边坡极限承载力分析下限法
CN108170898B (zh) * 2017-12-06 2021-12-17 昆明理工大学 一种节理岩质边坡可靠度分析下限法
CN108254782B (zh) * 2018-02-09 2019-11-05 中国地质大学（北京） 一种边坡地震破坏失稳概率的获取方法及系统

Also Published As

Publication number Publication date
CN111651902A (zh) 2020-09-11

Similar Documents

Publication Publication Date Title
CN111651902B (zh) 一种昔格达地层岩质边坡稳定性计算的拟动力下限法
Yagoda‐Biran et al. Constraining paleo PGA values by numerical analysis of overturned columns
Zhao et al. Three-dimensional distance potential discrete element method for the numerical simulation of landslides
Ning et al. A complex rockslide developed from a deep-seated toppling failure in the upper Lancang River, Southwest China
CN114297864B (zh) 一种受陡缓倾角控制的碎裂松动岩体边坡稳定性分析方法
CN102677792A (zh) 一种高耸结构风致振动的综合控制装置
CN108229056A (zh) 建筑抗震支吊架的抗震可靠度设计方法
Tang et al. Response of a pile group behind quay wall to liquefaction-induced lateral spreading: a shake-table investigation
Sarhosis et al. Three dimensional modelling of ancient colonnade structural systems subjected to harmonic and seismic loading
Kalani Sarokolayi et al. Dynamic analysis of elevated water storage tanks due to ground motions’ rotational and translational components
Yu et al. Centrifuge modeling of a dry sandy slope response to earthquake loading
Lu et al. Numerical and experimental analyses for bearing capacity of rigid strip footing subjected to eccentric load
CN111046480B (zh) 一种基于能量耗散原理的挡土墙地震土压力时程计算方法
CN111596356B (zh) 一种昔格达组地层岩质边坡的地震惯性力计算方法
CN111651900B (zh) 一种昔格达地层岩质边坡稳定性计算的拟动力上限法
CN111651901A (zh) 一种昔格达地层岩质边坡动力稳定性的时程分析方法
Dong et al. Analysis of deep dynamic sliding stability of gravity dam foundation based on DDA method
Hu et al. Response of beams under the impact of freak waves
Bestuzheva et al. Forms of natural vibrations of a gravitational dam on an elastic foundation, taking into account the attached mass of water
Zafeirakos et al. The role of soil and interface nonlinearities on the seismic response of caisson supported bridge piers
El Shamy et al. DEM simulation of the seismic response of gravity retaining walls
Dang et al. Dynamic response analysis of intake tower in hydroelectric power station with high surrounding rock
Klishin et al. Rolling Friction in Loose Media and its Role in Mechanics Problems
Li et al. 3D Valley Site Effect on the Nonlinear Seismic Response of Earth Dam
Markov et al. About effect of dynamic loading caused by concrete pump on bending moments, anchor forces and deformations in single level anchored retaining wall

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