CN111651901A - 一种昔格达地层岩质边坡动力稳定性的时程分析方法 - Google Patents

一种昔格达地层岩质边坡动力稳定性的时程分析方法 Download PDF

Info

Publication number
CN111651901A
CN111651901A CN202010610995.9A CN202010610995A CN111651901A CN 111651901 A CN111651901 A CN 111651901A CN 202010610995 A CN202010610995 A CN 202010610995A CN 111651901 A CN111651901 A CN 111651901A
Authority
CN
China
Prior art keywords
xigeda
rock slope
slope
strip
ith
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.)
Granted
Application number
CN202010610995.9A
Other languages
English (en)
Other versions
CN111651901B (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.)
China Nonferrous Metals Industry Kunming Survey And Design Institute Co ltd
Kunming University of Science and Technology
Original Assignee
China Nonferrous Metals Industry Kunming Survey And Design Institute Co ltd
Kunming University of Science and Technology
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 China Nonferrous Metals Industry Kunming Survey And Design Institute Co ltd, Kunming University of Science and Technology filed Critical China Nonferrous Metals Industry Kunming Survey And Design Institute Co ltd
Priority to CN202010610995.9A priority Critical patent/CN111651901B/zh
Publication of CN111651901A publication Critical patent/CN111651901A/zh
Application granted granted Critical
Publication of CN111651901B publication Critical patent/CN111651901B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/12Timing analysis or timing optimisation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/23Dune restoration or creation; Cliff stabilisation

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

一种昔格达地层岩质边坡动力稳定性的时程分析方法
技术领域
本发明涉及一种昔格达地层岩质边坡动力稳定性的时程分析方法,属于岩质边坡稳定性分析的技术领域。
背景技术
昔格达组地层广泛分布于我国攀枝花、西昌地区,其是河湖相沉积地层,形成于晚更新世与上第三系之间。昔格达组地层由不同颜色的粉砂岩和粘土岩互层的一套半成岩。昔格达组地层的粉砂岩和粘土岩均属于极软岩,其抗压强度、抗拉强度和抗剪强度较低、容易软化、极易崩解。攀昌地区碰到了大量昔格达地层岩质边坡,由于昔格达组地层的工程性质极差,昔格达地层岩质边坡的稳定性一直是工程界关注的焦点。
昔格达广泛分布的攀西地区是我国地震高发区,地震非常频繁,该地区的地震烈度达到7-9度。地震作用是昔格达岩质边坡失稳破坏的主要因素之一,地震作用下昔格达岩质边坡的动力稳定性是保证工程安全建设和运行的关键。要准确计算昔格达岩质边坡的动力稳定性需同时考虑地震波传播的时间效应和空间效应。
在当前,岩质边坡动力稳定性分析方法主要包括:拟静力法、拟动力法、基于振型叠加的反应谱法。但这几种方法都有其不足之处,比如:拟静力法可考虑边坡地震拟静力荷载随边坡高度的变化,但其不能考虑地震的“行波效应”且放大了地震荷载效应;拟动力法作为拟静力法的一种改进,可同时考虑地震波在岩体中传播的时间效应,但由于其将地震波假设为间谐波导致计算结果与实际地震的反应有一定差距;反应谱法基于动力频域分析,虽然在结合有限元方法等方面有优势,但不能适用于对时域分析。拟静力法和反应谱法均只能得到单一的边坡安全系数值,而不能获得边坡安全系数的时程曲线。
为了进一步完善对昔格达组地层岩质边坡的动力稳定性的研究,应该增加动力稳定性的时程分析方法,以填补现有方法的不足,并为昔格达组地层岩质边坡的动力稳定性的研究提供依据。
发明内容
本发明提供了一种昔格达地层岩质边坡动力稳定性的时程分析方法,以用于获得动力稳定性的时程曲线。
本发明的技术方案是:一种昔格达地层岩质边坡动力稳定性的时程分析方法,以昔格达地层岩质边坡为研究对象,首先假设昔格达地层岩质边坡的滑动面,并将滑动面以上的滑体划分成竖直条块系统,并进行条块的受力分析;基于岩质边坡坡脚处的地震时程加速度采用动力学原理计算每个条块的地震时程加速度;基于条块的动力平衡方程、条块底滑面和条间接触面的屈服条件,并结合目标函数,建立昔格达地层岩质边坡动力稳定性的非线性数学规划模型,计算昔格达地层岩质边坡动力稳定性安全系数,并绘制昔格达地层岩质边坡动力稳定性安全系数的时程曲线。
所述方法具体步骤包括:
步骤一、拟定昔格达地层岩质边坡动力稳定性计算的参数,用于参与步骤二、步骤三、步骤四;
步骤二、假设昔格达地层岩质边坡的滑动面,将滑动面以上的滑体划分为条块系统,并进行条块的受力分析;
步骤三、计算昔格达地层岩质边坡条块的地震时程加速度;
步骤四、建立昔格达地层岩质边坡动力稳定性的非线性数学规划模型:结合步骤二的条块的质量、步骤三的地震时程加速度、目标函数、条块的动力平衡方程、条块底滑面的屈服条件、条间接触面的屈服条件建立非线性数学规划模型;
步骤五、采用循环迭代方法求解昔格达地层岩质边坡动力稳定性的非线性数学规划模型,获得昔格达地层岩质边坡的动力稳定性的安全系数;
步骤六、绘制昔格达地层岩质边坡动力稳定性安全系数的时程曲线。
所述拟定昔格达地层岩质边坡动力稳定性计算的参数包括:①确定昔格达地层岩质边坡的几何参数;②确定昔格达地层中粉砂岩地层、粘土岩地层的分布情况;③确定昔格达地层岩体的材料参数;④确定昔格达地层岩质边坡的地震动参数;⑤确定昔格达地层岩质边坡的坡脚处的地震时程加速度数据。
所述昔格达地层岩质边坡的几何参数包括:岩质边坡的高度H,边坡几何形状控制点的位置坐标;地层的分布情况包括:粉砂岩和粘土岩地层的数量,粉砂岩和粘土岩地层每一层的厚度,地层的倾角;材料参数包括:昔格达地层岩体的容重γ,岩体的内摩擦角
Figure BDA0002560919470000021
和岩体的凝聚力cs;地震动参数包括:地震振动时长Ts、地震设计烈度;地震时程加速度数据包括:昔格达地层岩质边坡的坡脚处输入的水平方向地震时程加速度αH(tj),昔格达地层岩质边坡的坡脚处输入的竖直方向地震时程加速度αV(tj),其中j=(1,…,Nt),Nt是地震振动时长内的时步数量。
所述受力分析具体为:
(1)昔格达地层岩质边坡的坐标系,以边坡的坡脚为坐标原点,水平轴为坐标系的x轴,x轴向右为正,竖直轴为坐标系的y轴,y轴向上为正;
(2)假设昔格达地层岩质边坡的滑动面为折线滑动面,即为EFG的连线,并确定滑动面关键点的位置坐标;
(3)将滑动面以上的滑体划分成竖直条块,具体为:将EFG折线滑动面以上的滑体区域EFGI划分成竖直的条块i=(1,…,Ns),条块的数量为Ns;其中,条块是由按逆时针排列的A、B、C、D四个关键点形成的四边形/三角形:当第一个条块时,A和D两点重合,条块变成三角形,当最后一个条块时B和C两点重合,条块变成三角形,其余为四边形;第i个条块的形心为oi,第i个条块底滑面的倾角θi
Figure BDA0002560919470000031
是第i个条块底滑面的长度;
Figure BDA0002560919470000032
是第i个条块左侧条间接触面的长度;
Figure BDA0002560919470000033
是第i个右侧条间接触面的长度,当i=1时
Figure BDA0002560919470000034
当i=Ns
Figure BDA0002560919470000035
(4)计算昔格达地层岩质边坡条块的质量:
mi=Siγ/g
式中:mi是第i个条块的质量,Si是第i个条块的面积,γ是昔格达地层岩体的容重;g为重力加速度;
(5)条块的受力分析
条块是由按逆时针排列的A、B、C、D四个关键点形成的四边形/三角形:当i=1时A和D两点重合,条块变成三角形,左侧条间接触面退化成一个点,当i=Ns时B和C两点重合,条块变成三角形,右侧侧条间接触面退化成一个点;其余为四边形,包含底滑面AB、左侧条间接触面AD、右侧条间接触面BC、顶面CD;假设底滑面和条间接触面的法向力以受压为正、剪力以对条块产生逆时针的转动效果为正,条块受力分析如下:
①昔格达地层岩质边坡第i个条块的底滑面AB作用有两个力:底滑面的法向力
Figure BDA0002560919470000036
底滑面的剪力
Figure BDA0002560919470000037
②昔格达地层岩质边坡第i个左侧条间接触面AD作用有两个力:第i个条块左侧条间接触面的法向力
Figure BDA0002560919470000041
第i个条块左侧条间接触面的剪力
Figure BDA0002560919470000042
当i=1时
Figure BDA0002560919470000043
③昔格达地层岩质边坡第i个右侧条间接触面BC作用有两个力:第i个条块右侧条间接触面的法向力
Figure BDA0002560919470000044
第i个条块右侧条间接触面的剪力
Figure BDA0002560919470000045
当i=Ns
Figure BDA0002560919470000046
④昔格达地层岩质边坡第i个条块形心oi处作用有三个力:第i个条块的自重Gi、第j个时步昔格达地层岩质边坡第i个条块的水平向地震时程惯性力
Figure BDA0002560919470000047
第j个时步昔格达地层岩质边坡第i个条块的竖直向地震时程惯性力
Figure BDA0002560919470000048
所述地震时程加速度计算具体为:
(1)地震作用下昔格达地层岩质边坡每个时步对应的计算时间按下式计算:
Figure BDA0002560919470000049
式中,tj是第j个时步对应的计算时间,j=(1,…,Nt),Nt是地震振动时长内的时步数量,Ts为地震的振动时长;
(2)昔格达地层岩质边坡条块的水平方向地震时程加速度按下式计算:
Figure BDA00025609194700000410
式中:
Figure BDA00025609194700000411
是第j个时步昔格达地层岩质边坡第i个条块沿水平方向的地震时程加速度;i=(1,…,Ns),Ns是滑体条块的数量;ξ为地震作用效应的折减系数;αH(tj)是第j个时步昔格达地层岩质边坡的坡脚处输入的水平方向地震时程加速度;λm是地震加速度最大分布系数;Hi是第i个条块的形心至边坡的坡脚处的垂直高度;H是昔格达地层岩质边坡的高度;
(3)昔格达地层岩质边坡条块的竖直方向地震时程加速度按下式计算:
Figure BDA00025609194700000412
式中:
Figure BDA00025609194700000413
是第j个时步昔格达地层岩质边坡第i个条块沿竖直方向的地震时程加速度;αV(tj)是第j个时步昔格达地层岩质边坡的坡脚处输入的竖直方向地震时程加速度。
所述建立昔格达地层岩质边坡动力稳定性的非线性数学规划模型具体为:
(1)建立目标函数
设昔格达地层岩质边坡的安全系数为目标函数,并求解边坡安全系数的最大值,目标函数如下:
Maximize:Kj
式中:Kj是第j个时步昔格达地层岩质边坡动力稳定性的安全系数,j=(1,…,Nt),Nt是地震振动时长内的时步数量;Maximize表示“使最大”;
(2)条块的动力平衡方程:条块受到自重、地震力、底滑面和条间接触面上的法向力、剪力的作用而保持平衡,条块的动力平衡方程按以下流程建立:
①建立第j个时步昔格达地层岩质边坡第i个条块的水平方向的动力平衡方程,具体如下:
Figure BDA0002560919470000051
式中:i=(1,…,Ns),Ns是条块的数量;θi第i个条块底滑面的倾角;
Figure BDA0002560919470000052
是昔格达地层岩质边坡第i个条块底滑面的法向力,
Figure BDA0002560919470000053
以受压为正;
Figure BDA0002560919470000054
是昔格达地层岩质边坡第i个条块底滑面的剪力,
Figure BDA0002560919470000055
以对条块产生逆时针的转动效果为正;
Figure BDA0002560919470000056
是昔格达地层岩质边坡第i个条块左侧条间接触面的法向力,
Figure BDA0002560919470000057
以受压为正,当i=1时
Figure BDA0002560919470000058
Figure BDA0002560919470000059
是昔格达地层岩质边坡第i个条块右侧条间接触面的法向力,
Figure BDA00025609194700000510
以受压为正,当i=Ns
Figure BDA00025609194700000511
是第j个时步昔格达地层岩质边坡第i个条块的水平向地震时程惯性力,
Figure BDA00025609194700000512
mi是第i个条块的质量,
Figure BDA00025609194700000513
是第j个时步昔格达地层岩质边坡第i个条块沿水平方向的地震时程加速度;
②建立第j个时步昔格达地层岩质边坡第i个条块的竖直方向的动力平衡方程,具体如下:
Figure BDA00025609194700000514
式中:
Figure BDA0002560919470000061
是昔格达地层岩质边坡第i个条块左侧条间接触面的剪力,
Figure BDA0002560919470000062
以对条块产生逆时针的转动效果为正,当i=1时
Figure BDA0002560919470000063
是昔格达地层岩质边坡第i个条块右侧条间接触面的剪力,
Figure BDA0002560919470000064
以对条块产生逆时针的转动效果为正,当i=Ns
Figure BDA0002560919470000065
Gi是第i个条块的自重,Gi=miγ,γ是昔格达地层岩体的容重;
Figure BDA0002560919470000066
是第j个时步昔格达地层岩质边坡第i个条块的竖直向地震时程惯性力,
Figure BDA0002560919470000067
是第j个时步昔格达地层岩质边坡第i个条块沿竖直方向的地震时程加速度;
(3)建立条块底滑面的屈服条件,具体如下:
Figure BDA0002560919470000068
式中:
Figure BDA0002560919470000069
是昔格达地层岩质边坡第i个条块底滑面的长度;
Figure BDA00025609194700000610
是岩体的内摩擦角,cs是岩体的凝聚力;
(4)建立左侧条间接触面和右侧条间接触面的屈服条件,具体如下:
Figure BDA00025609194700000611
式中:
Figure BDA00025609194700000612
是昔格达地层岩质边坡第i个条块左侧条间接触面的长度;
Figure BDA00025609194700000613
是昔格达地层岩质边坡第i个条块右侧条间接触面的长度;
(5)建立昔格达地层岩质边坡动力稳定性的非线性数学规划模型,具体为:
Figure BDA00025609194700000614
式中:Si是第i个条块的面积,g为重力加速度;ξ为地震作用效应的折减系数;αH(tj)是第j个时步昔格达地层岩质边坡的坡脚处输入的水平方向地震时程加速度;λm是地震加速度最大分布系数;Hi是第i个条块的形心至边坡的坡脚处的垂直高度;H是昔格达地层岩质边坡的高度;αV(tj)是第j个时步昔格达地层岩质边坡的坡脚处输入的竖直方向地震时程加速度。
本发明的有益效果是:本发明基于动力学理论、极限平衡理论和非线性数学规划理论,研究地震作用下昔格达地层岩质边坡的动力稳定性;本发明基于岩质边坡坡脚处的地震时程加速度采用动力学原理计算条块的地震时程加速度,建立昔格达地层岩质边坡动力稳定性的非线性数学规划模型;使用循环迭代方法计算昔格达地层岩质边坡动力稳定性安全系数,获得边坡动力稳定性安全系数的时程曲线。本发明方法具有概念明确、计算效率和计算精度高等特点,可将其应用于昔格达地层岩质边坡的动力稳定性分析;同时使用动力学理论和非线性规划理论求解昔格达地层岩质边坡的动力稳定性,编制程序简单、计算精度较高;整个方法理论严谨、工程应用简便,可将其应用于昔格达地层岩质边坡的动力稳定性计算领域,可为岩质边坡的抗震设计提供科学依据。
附图说明
图1为昔格达地层岩质边坡的滑动面示意图;
图2为昔格达地层岩质边坡的滑体划分条块示意图;
图3为昔格达地层岩质边坡的条块受力分析示意图;
图4为本发明技术流程图;
图5为实施例昔格达地层岩质边坡的几何形状示意图(单位:m);
图6为实施例昔格达地层岩质边坡的坡脚处输入的水平向地震时程加速度时程曲线;
图7为实施例昔格达地层岩质边坡的坡脚处输入的竖直向地震时程加速度时程曲线;
图8为实施例昔格达地层岩质边坡的条块划分示意图;
图9为实施例昔格达地层岩质边坡的第1个条块沿水平方向的地震时程加速度时程曲线;
图10为实施例昔格达地层岩质边坡的第1个条块沿竖直方向的地震时程加速度时程曲线;
图11为实施例昔格达地层岩质边坡的第2个条块沿水平方向的地震时程加速度时程曲线;
图12为实施例昔格达地层岩质边坡的第2个条块沿竖直方向的地震时程加速度时程曲线;
图13为实施例昔格达地层岩质边坡的第3个条块沿水平方向的地震时程加速度时程曲线;
图14为实施例昔格达地层岩质边坡的第3个条块沿竖直方向的地震时程加速度时程曲线;
图15为实施例昔格达地层岩质边坡的第4个条块沿水平方向的地震时程加速度时程曲线;
图16为实施例昔格达地层岩质边坡的第4个条块沿竖直方向的地震时程加速度时程曲线;
图17为实施例昔格达地层岩质边坡的第5个条块沿水平方向的地震时程加速度时程曲线;
图18为实施例昔格达地层岩质边坡的第5个条块沿竖直方向的地震时程加速度时程曲线;
图19为实施例昔格达地层岩质边坡的第6个条块沿水平方向的地震时程加速度时程曲线;
图20为实施例昔格达地层岩质边坡的第6个条块沿竖直方向的地震时程加速度时程曲线;
图21为实施例昔格达地层岩质边坡的第7个条块沿水平方向的地震时程加速度时程曲线;
图22为实施例昔格达地层岩质边坡的第7个条块沿竖直方向的地震时程加速度时程曲线;
图23为实施例昔格达地层岩质边坡的第8个条块沿水平方向的地震时程加速度时程曲线;
图24为实施例昔格达地层岩质边坡的第8个条块沿竖直方向的地震时程加速度时程曲线;
图25为实施例昔格达地层岩质边坡的安全系数时程曲线。
具体实施方式
下面结合附图和实施例,对本发明作进一步说明,但本发明的内容并不限于所述范围。
实施例1:如图1-25所示,一种昔格达地层岩质边坡动力稳定性的时程分析方法,以昔格达地层岩质边坡为研究对象,首先假设昔格达地层岩质边坡的滑动面,并将滑动面以上的滑体划分成竖直条块系统,并进行条块的受力分析;基于岩质边坡坡脚处的地震时程加速度采用动力学原理计算每个条块的地震时程加速度;基于条块的动力平衡方程、条块底滑面和条间接触面的屈服条件,并结合目标函数,建立昔格达地层岩质边坡动力稳定性的非线性数学规划模型,计算昔格达地层岩质边坡动力稳定性安全系数,并绘制昔格达地层岩质边坡动力稳定性安全系数的时程曲线。流程图,如图4所示。
进一步地,可以设置所述方法具体步骤包括:步骤一、拟定昔格达地层岩质边坡动力稳定性计算的参数,用于参与步骤二、步骤三、步骤四;步骤二、假设昔格达地层岩质边坡的滑动面,将滑动面以上的滑体划分为条块系统,并进行条块的受力分析;步骤三、计算昔格达地层岩质边坡条块的地震时程加速度;步骤四、建立昔格达地层岩质边坡动力稳定性的非线性数学规划模型:结合步骤二的条块的质量、步骤三的地震时程加速度、目标函数、条块的动力平衡方程、条块底滑面的屈服条件、条间接触面的屈服条件建立非线性数学规划模型;步骤五、采用循环迭代方法求解昔格达地层岩质边坡动力稳定性的非线性数学规划模型,获得昔格达地层岩质边坡的动力稳定性的安全系数;步骤六、绘制昔格达地层岩质边坡动力稳定性安全系数的时程曲线。
更进一步地,本发明结合具体步骤,给出如下过程:
步骤一、拟定昔格达地层岩质边坡动力稳定性计算的参数。
根据昔格达地层岩质边坡的实际情况,拟定其动力稳定性计算的参数,具体为:
①确定昔格达地层岩质边坡的几何参数,包括:边坡的高度H=50m,边坡几何形状控制点的位置坐标如图5所示,边坡的4个控制点的位置坐标分别为:O点的位置坐标为(0,0)、P点的位置坐标为(80,0)、Q点的位置坐标为(80,50)、R点的位置坐标为(28.86,50)。
②确定昔格达地层中粉砂岩地层、粘土岩地层的分布情况如图5所示,包括:粉砂岩和粘土岩地层的数量9层,第1层粉砂岩的厚度为4.86m,第2层至第9层的粉砂岩和粘土岩的厚度均为8.0m,地层的倾角为30°。
③确定昔格达地层岩体的材料参数,包括:昔格达地层岩体的容重γ为2500kN/m3,岩体的内摩擦角
Figure BDA0002560919470000108
为25°,岩体的凝聚力cs为160kPa。
④确定昔格达地层岩质边坡的地震动参数,包括:地震振动时长T为29.8s;地震设计烈度为7度;(其中,当地震设计烈度为7度时λm=3.0、地震设计烈度为8度时λm=2.5、地震设计烈度为9度时λm=2.0);
⑤确定昔格达地层岩质边坡的坡脚处的地震时程加速度数据,包括:昔格达地层岩质边坡的坡脚处输入的水平方向地震时程加速度αH(tj)绘制成时程曲线如图6所示,αH(tj)共计1490个点,j=(1,…,1490);昔格达地层岩质边坡的坡脚处输入的竖直方向地震时程加速度αV(tj)绘制成时程曲线如图7所示,其中j=(1,…,1490),地震振动时长内的时步数量Nt为1490。
步骤二、假设昔格达地层岩质边坡的滑动面,将滑动面以上的滑体划分为条块系统,并进行条块的受力分析。
具体为:
(1)昔格达地层岩质边坡的坐标系,以边坡的坡脚为坐标原点,水平轴为坐标系的x轴,x轴向右为正,竖直轴为坐标系的y轴,y轴向上为正,坐标系如图5中所示。
(2)假设昔格达地层岩质边坡失稳的滑动面。具体为:如图8所示,假设昔格达地层岩质边坡的滑动面为折线滑动面,即为EFG的连线,滑动面关键点的位置坐标为:E点的位置坐标为(8,13.86)、F点的位置坐标为(39.3,31.93)、G点的位置坐标为(51.34,50)。
(3)将滑动面以上的滑体划分成竖直条块。具体为:如图8所示,将EFG折线滑动面以上的滑体区域EFGR划分成竖直的8个条块,条块的数量为Ns=8。第i个条块的示意图如图3所示,i=(1,…,Ns),条块是由按逆时针排列的A、B、C、D四个关键点形成四边形,当第一个条块i=1时A和D两点重合,条块变成三角形,当最后一个条块i=Ns时B和C两点重合,条块变成三角形;第i个条块的形心为oi,第i个条块底滑面的倾角θi
Figure BDA0002560919470000101
是第i个条块底滑面的长度;
Figure BDA0002560919470000102
是第i个条块左侧条间接触面的长度;
Figure BDA0002560919470000103
是第i个右侧条间接触面的长度,当i=1时
Figure BDA0002560919470000104
当i=Ns
Figure BDA0002560919470000105
具体的图8中数据如下:
θ1=30°、θ2=30°、θ3=30°、θ4=30°、θ5=30°、θ6=30°、θ7=56°、θ8=56°;
Figure BDA0002560919470000106
Figure BDA0002560919470000107
Figure BDA0002560919470000111
(4)计算昔格达地层岩质边坡条块的质量。昔格达地层岩质边坡第i个条块的质量按下式计算:
mi=Siγ/g(1)
式中:i=(1,…,Ns),Ns是滑体条块的数量,mi是第i个条块的质量,Si是第i个条块的面积,γ是昔格达地层岩体的容重;g为重力加速度,取9.81m/s2
(5)条块的受力分析
条块是由按逆时针排列的A、B、C、D四个关键点形成的四边形/三角形:当i=1时A和D两点重合,条块变成三角形,左侧条间接触面退化成一个点,当i=Ns时B和C两点重合,条块变成三角形,右侧侧条间接触面退化成一个点;其余为四边形,包含底滑面AB、左侧条间接触面AD、右侧条间接触面BC、顶面CD;假设底滑面和条间接触面的法向力以受压为正、剪力以对条块产生逆时针的转动效果为正,条块受力分析如下:
①昔格达地层岩质边坡第i个条块的底滑面AB作用有两个力:底滑面的法向力
Figure BDA0002560919470000112
底滑面的剪力
Figure BDA0002560919470000113
②昔格达地层岩质边坡第i个左侧条间接触面AD作用有两个力:第i个条块左侧条间接触面的法向力
Figure BDA0002560919470000114
第i个条块左侧条间接触面的剪力
Figure BDA0002560919470000115
当i=1时
Figure BDA0002560919470000116
③昔格达地层岩质边坡第i个右侧条间接触面BC作用有两个力:第i个条块右侧条间接触面的法向力
Figure BDA0002560919470000117
第i个条块右侧条间接触面的剪力
Figure BDA0002560919470000118
当i=Ns
Figure BDA0002560919470000119
④昔格达地层岩质边坡第i个条块形心oi处作用有三个力:第i个条块的自重Gi、第j个时步昔格达地层岩质边坡第i个条块的水平向地震时程惯性力
Figure BDA00025609194700001110
第j个时步昔格达地层岩质边坡第i个条块的竖直向地震时程惯性力
Figure BDA00025609194700001111
步骤三、计算昔格达地层岩质边坡条块的地震时程加速度。
具体为:
地震作用下昔格达地层岩质边坡每个时步对应的计算时间按下式计算:
Figure BDA0002560919470000121
式中,tj是第j个时步对应的计算时间,j=(1,…,Nt),Nt是地震振动时长内的时步数量,Nt取1490,Δt是时步的时长,Δt=Ts/Nt=0.02s,Ts为地震的振动时长取29.80s。
昔格达地层岩质边坡条块的水平方向地震时程加速度按下式计算:
Figure BDA0002560919470000122
上式中:i=(1,…,Ns),Ns是滑体条块的数量,Ns取8;j=(1,…,Nt),Nt是地震振动时长内的时步数量,Nt取1490;
Figure BDA0002560919470000123
是第j个时步昔格达地层岩质边坡第i个条块沿水平方向的地震时程加速度;αH(tj)是第j个时步昔格达地层岩质边坡的坡脚处输入的水平方向地震时程加速度;ξ为地震作用效应的折减系数,ξ取0.25;Hi是第i个条块的形心至边坡的坡脚处的垂直高度;H是昔格达地层岩质边坡的高度,H取50m;λm是地震加速度最大分布系数,实施例的地震设计烈度为7度,λm=3.0。
昔格达地层岩质边坡条块的竖直方向地震时程加速度按下式计算:
Figure BDA0002560919470000124
上式中:i=(1,…,Ns),Ns是滑体条块的数量,Ns取8;j=(1,…,Nt),Nt是地震振动时长内的时步数量,Nt取1490;
Figure BDA0002560919470000125
是第j个时步昔格达地层岩质边坡第i个条块沿竖直方向的地震时程加速度;αV(tj)是第j个时步昔格达地层岩质边坡的坡脚处输入的竖直方向地震时程加速度;ξ为地震作用效应的折减系数,ξ取0.25;Hi是第i个条块的形心至边坡的坡脚处的垂直高度;H是昔格达地层岩质边坡的高度,H取50m;λm是地震加速度最大分布系数,实施例的地震设计烈度为7度,λm=3.0。
将第1至第8个条块的沿水平方向地震时程加速度
Figure BDA0002560919470000126
沿竖直方向的地震时程加速度
Figure BDA0002560919470000127
绘制成时程曲线分别如图9~图24所示。
步骤四、建立昔格达地层岩质边坡动力稳定性的非线性数学规划模型。具体为:
(1)建立目标函数
本发明设昔格达地层岩质边坡的安全系数为目标函数,并求解边坡安全系数的最大值。目标函数如下:
Maximize:Kj(5)
上式中:Kj是第j个时步昔格达地层岩质边坡动力稳定性的安全系数,j=(1,…,Nt),Nt是地震振动时长内的时步数量;Maximize表示“使最大”。
(2)条块的动力平衡方程;
条块受到自重、地震力、底滑面和条间接触面上的法向力、剪力的作用而保持平衡,条块的动力平衡方程按以下流程建立:
①建立第j个时步昔格达地层岩质边坡第i个条块的水平方向的动力平衡方程。
具体如下:
Figure BDA0002560919470000131
上式中:i=(1,…,Ns),Ns是条块的数量;j=(1,…,Nt),Nt是地震振动时长内的时步数量;θi第i个条块底滑面的倾角;
Figure BDA0002560919470000132
是昔格达地层岩质边坡第i个条块底滑面的法向力,
Figure BDA0002560919470000133
以受压为正;
Figure BDA0002560919470000134
是昔格达地层岩质边坡第i个条块底滑面的剪力,
Figure BDA0002560919470000135
以对条块产生逆时针的转动效果为正;
Figure BDA0002560919470000136
是昔格达地层岩质边坡第i个条块左侧条间接触面的法向力,
Figure BDA0002560919470000137
以受压为正,当i=1时
Figure BDA0002560919470000138
是昔格达地层岩质边坡第i个条块右侧条间接触面的法向力,
Figure BDA0002560919470000139
以受压为正,当i=Ns
Figure BDA00025609194700001310
是第j个时步昔格达地层岩质边坡第i个条块的水平向地震时程惯性力;
Figure BDA00025609194700001311
mi是第i个条块的质量,
Figure BDA00025609194700001312
是第j个时步昔格达地层岩质边坡第i个条块沿水平方向的地震时程加速度。
②建立第j个时步昔格达地层岩质边坡第i个条块的竖直方向的动力平衡方程。
具体如下:
Figure BDA00025609194700001313
上式中:i=(1,…,Ns),Ns是条块的数量;j=(1,…,Nt),Nt是地震振动时长内的时步数量;θi第i个条块底滑面的倾角;
Figure BDA00025609194700001314
是昔格达地层岩质边坡第i个条块底滑面的法向力,
Figure BDA0002560919470000141
以受压为正;
Figure BDA0002560919470000142
是昔格达地层岩质边坡第i个条块底滑面的剪力,
Figure BDA0002560919470000143
以对条块产生逆时针的转动效果为正;
Figure BDA0002560919470000144
是昔格达地层岩质边坡第i个条块左侧条间接触面的剪力,
Figure BDA0002560919470000145
以对条块产生逆时针的转动效果为正,当i=1时
Figure BDA0002560919470000146
是昔格达地层岩质边坡第i个条块右侧条间接触面的剪力,
Figure BDA0002560919470000147
以对条块产生逆时针的转动效果为正,当i=Ns
Figure BDA0002560919470000148
Gi是第i个条块的自重,Gi=miγ,mi是第i个条块的质量;γ是昔格达地层岩体的容重;
Figure BDA0002560919470000149
是第j个时步昔格达地层岩质边坡第i个条块的竖直向地震时程惯性力,
Figure BDA00025609194700001410
是第j个时步昔格达地层岩质边坡第i个条块沿竖直方向的地震时程加速度。
(3)建立条块底滑面的屈服条件。具体如下:
Figure BDA00025609194700001411
上式中:i=(1,…,Ns),Ns是条块的数量;j=(1,…,Nt),Nt是地震振动时长内的时步数量;
Figure BDA00025609194700001412
是昔格达地层岩质边坡第i个条块底滑面的法向力,
Figure BDA00025609194700001413
以受压为正;
Figure BDA00025609194700001414
是昔格达地层岩质边坡第i个条块底滑面的剪力,
Figure BDA00025609194700001415
以对条块产生逆时针的转动效果为正;
Figure BDA00025609194700001416
是昔格达地层岩质边坡第i个条块底滑面的长度;
Figure BDA00025609194700001417
是岩体的内摩擦角,cs是岩体的凝聚力;Kj是第j个时步昔格达地层岩质边坡动力稳定性的安全系数。
(4)建立左侧条间接触面和右侧条间接触面的屈服条件。具体如下:
Figure BDA00025609194700001418
上式中:i=(1,…,Ns),Ns是条块的数量;j=(1,…,Nt),Nt是地震振动时长内的时步数量;
Figure BDA00025609194700001419
是昔格达地层岩质边坡第i个条块左侧条间接触面的法向力,
Figure BDA00025609194700001420
以受压为正,当i=1时
Figure BDA00025609194700001421
是昔格达地层岩质边坡第i个条块左侧条间接触面的剪力,
Figure BDA00025609194700001422
以对条块产生逆时针的转动效果为正,当i=1时
Figure BDA00025609194700001423
是昔格达地层岩质边坡第i个条块右侧条间接触面的法向力,
Figure BDA00025609194700001424
以受压为正,当i=Ns
Figure BDA00025609194700001425
Figure BDA00025609194700001426
是昔格达地层岩质边坡第i个条块右侧条间接触面的剪力,
Figure BDA00025609194700001427
以对条块产生逆时针的转动效果为正,当i=Ns
Figure BDA0002560919470000151
是昔格达地层岩质边坡第i个条块左侧条间接触面的长度;
Figure BDA0002560919470000152
是昔格达地层岩质边坡第i个条块右侧条间接触面的长度;
Figure BDA0002560919470000153
是岩体的内摩擦角,cs是岩体的凝聚力;Kj是第j个时步昔格达地层岩质边坡动力稳定性的安全系数。
(5)建立昔格达地层岩质边坡动力稳定性的非线性数学规划模型。具体为:
将目标函数、条块动力平衡方程、底滑面和条间接触面的屈服条件集成,可得到昔格达地层岩质边坡动力稳定性的非线性数学规划模型如下:
Figure BDA0002560919470000154
上式中:i=(1,…,Ns),Ns是条块的数量,Ns取8;j=(1,…,Nt),Nt是地震振动时长内的时步数量,Nt取1490。
步骤五、昔格达地层岩质边坡动力稳定性的时程分析循环迭代求解。具体为:
昔格达地层岩质边坡的动力非线性数学规划模型式(10)是一个非线性数学规划模型,其目标函数是步昔格达地层岩质边坡动力稳定性的安全系数Kj,决策变量是
Figure BDA0002560919470000155
本申请基于循环迭代的思想求解昔格达地层岩质边坡的动力稳定性的安全系数Kj,具体数值求解方法如下:
(1)根据式(1),计算昔格达地层岩质边坡滑体的每个条块的质量mi,并根据式(7)中的Gi=miγ计算自重Gi,其中i=(1,…,Ns),Ns是条块的数量;
(2)根据式(2)计算昔格达地层岩质边坡每个时步对应的计算时间tj,其中j=(1,…,Nt),Nt是地震振动时长内的时步数量;
(3)根据式(3)、式(4),结合αH(tj)、αV(tj)、λm、H、Hi、mi分别计算每个条块沿水平方向和沿竖直方向的地震时程加速度
Figure BDA0002560919470000161
进而根据式(6)中的
Figure BDA0002560919470000162
式(7)中的
Figure BDA0002560919470000163
分别计算每个条块沿水平方向和沿竖直方向的地震时程惯性力
Figure BDA0002560919470000164
其中i=(1,…,Ns),Ns是条块的数量,j=(1,…,Nt),j=(1,…,Nt),Nt是地震振动时长内的时步数量。
(4)将Gi
Figure BDA0002560919470000165
cs
Figure BDA0002560919470000166
从j=1到j=Nt循环带入昔格达地层岩质边坡动力稳定性的非线性数学规划模型式(10),并使用“内点算法”求解昔格达地层岩质边坡动力稳定性的非线性数学规划模型,求解得到Nt个昔格达地层岩质边坡动力稳定性的安全系数Kj,并获得与安全系数Kj对应的底滑面、条间接触面的法向力和剪力
Figure BDA0002560919470000167
步骤六、绘制昔格达地层岩质边坡动力稳定性安全系数的时程曲线。具体如下:
(1)以时间为水平轴,安全系数为竖直轴,形成时间-安全系数坐标系如图25所示;
(2)将昔格达地层岩质边坡每个时步j和对应的安全系数Kj结合起来,形成Nt个坐标点(j,Kj),其中j=(1,…,Nt);
(3)从j=1到j=Nt循环,将坐标点(j,Kj)逐次绘制在时间-安全系数的坐标系中,并依次从连线形成昔格达地层岩质边坡动力稳定性的安全系数时程曲线,实施例昔格达地层岩质边坡动力稳定性的安全系数时程曲线图25所示。
上面结合附图对本发明的具体实施方式作了详细说明,但是本发明并不限于上述实施方式,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下作出各种变化。

Claims (7)

1.一种昔格达地层岩质边坡动力稳定性的时程分析方法,其特征在于:以昔格达地层岩质边坡为研究对象,首先假设昔格达地层岩质边坡的滑动面,并将滑动面以上的滑体划分成竖直条块系统,并进行条块的受力分析;基于岩质边坡坡脚处的地震时程加速度采用动力学原理计算每个条块的地震时程加速度;基于条块的动力平衡方程、条块底滑面和条间接触面的屈服条件,并结合目标函数,建立昔格达地层岩质边坡动力稳定性的非线性数学规划模型,计算昔格达地层岩质边坡动力稳定性安全系数,并绘制昔格达地层岩质边坡动力稳定性安全系数的时程曲线。
2.根据权利要求1所述的昔格达地层岩质边坡动力稳定性的时程分析方法,其特征在于:所述方法具体步骤包括:
步骤一、拟定昔格达地层岩质边坡动力稳定性计算的参数,用于参与步骤二、步骤三、步骤四;
步骤二、假设昔格达地层岩质边坡的滑动面,将滑动面以上的滑体划分为条块系统,并进行条块的受力分析;
步骤三、计算昔格达地层岩质边坡条块的地震时程加速度;
步骤四、建立昔格达地层岩质边坡动力稳定性的非线性数学规划模型:结合步骤二的条块的质量、步骤三的地震时程加速度、目标函数、条块的动力平衡方程、条块底滑面的屈服条件、条间接触面的屈服条件建立非线性数学规划模型;
步骤五、采用循环迭代方法求解昔格达地层岩质边坡动力稳定性的非线性数学规划模型,获得昔格达地层岩质边坡的动力稳定性的安全系数;
步骤六、绘制昔格达地层岩质边坡动力稳定性安全系数的时程曲线。
3.根据权利要求2所述的昔格达地层岩质边坡动力稳定性的时程分析方法,其特征在于:所述拟定昔格达地层岩质边坡动力稳定性计算的参数包括:①确定昔格达地层岩质边坡的几何参数;②确定昔格达地层中粉砂岩地层、粘土岩地层的分布情况;③确定昔格达地层岩体的材料参数;④确定昔格达地层岩质边坡的地震动参数;⑤确定昔格达地层岩质边坡的坡脚处的地震时程加速度数据。
4.根据权利要求3所述的昔格达地层岩质边坡动力稳定性的时程分析方法,其特征在于:所述昔格达地层岩质边坡的几何参数包括:岩质边坡的高度H,边坡几何形状控制点的位置坐标;地层的分布情况包括:粉砂岩和粘土岩地层的数量,粉砂岩和粘土岩地层每一层的厚度,地层的倾角;材料参数包括:昔格达地层岩体的容重γ,岩体的内摩擦角
Figure FDA0002560919460000021
和岩体的凝聚力cs;地震动参数包括:地震振动时长Ts、地震设计烈度;地震时程加速度数据包括:昔格达地层岩质边坡的坡脚处输入的水平方向地震时程加速度αH(tj),昔格达地层岩质边坡的坡脚处输入的竖直方向地震时程加速度αV(tj),其中j=(1,…,Nt),Nt是地震振动时长内的时步数量。
5.根据权利要求1或2所述的昔格达地层岩质边坡动力稳定性的时程分析方法,其特征在于:所述受力分析具体为:
(1)昔格达地层岩质边坡的坐标系,以边坡的坡脚为坐标原点,水平轴为坐标系的x轴,x轴向右为正,竖直轴为坐标系的y轴,y轴向上为正;
(2)假设昔格达地层岩质边坡的滑动面为折线滑动面,即为EFG的连线,并确定滑动面关键点的位置坐标;
(3)将滑动面以上的滑体划分成竖直条块,具体为:将EFG折线滑动面以上的滑体区域EFGI划分成竖直的条块i=(1,…,Ns),条块的数量为Ns;其中,条块是由按逆时针排列的A、B、C、D四个关键点形成的四边形/三角形:当第一个条块时,A和D两点重合,条块变成三角形,当最后一个条块时B和C两点重合,条块变成三角形,其余为四边形;第i个条块的形心为oi,第i个条块底滑面的倾角θi
Figure FDA0002560919460000022
是第i个条块底滑面的长度;
Figure FDA0002560919460000023
是第i个条块左侧条间接触面的长度;
Figure FDA0002560919460000024
是第i个右侧条间接触面的长度,当i=1时
Figure FDA0002560919460000025
当i=Ns
Figure FDA0002560919460000026
(4)计算昔格达地层岩质边坡条块的质量:
mi=Siγ/g
式中:mi是第i个条块的质量,Si是第i个条块的面积,γ是昔格达地层岩体的容重;g为重力加速度;
(5)条块的受力分析
条块是由按逆时针排列的A、B、C、D四个关键点形成的四边形/三角形:当i=1时A和D两点重合,条块变成三角形,左侧条间接触面退化成一个点,当i=Ns时B和C两点重合,条块变成三角形,右侧侧条间接触面退化成一个点;其余为四边形,包含底滑面AB、左侧条间接触面AD、右侧条间接触面BC、顶面CD;假设底滑面和条间接触面的法向力以受压为正、剪力以对条块产生逆时针的转动效果为正,条块受力分析如下:
①昔格达地层岩质边坡第i个条块的底滑面AB作用有两个力:底滑面的法向力
Figure FDA0002560919460000031
底滑面的剪力Fi d
②昔格达地层岩质边坡第i个左侧条间接触面AD作用有两个力:第i个条块左侧条间接触面的法向力
Figure FDA0002560919460000032
第i个条块左侧条间接触面的剪力Fi L,当i=1时
Figure FDA0002560919460000033
③昔格达地层岩质边坡第i个右侧条间接触面BC作用有两个力:第i个条块右侧条间接触面的法向力
Figure FDA0002560919460000034
第i个条块右侧条间接触面的剪力Fi R,当i=Ns
Figure FDA0002560919460000035
④昔格达地层岩质边坡第i个条块形心oi处作用有三个力:第i个条块的自重Gi、第j个时步昔格达地层岩质边坡第i个条块的水平向地震时程惯性力
Figure FDA0002560919460000036
第j个时步昔格达地层岩质边坡第i个条块的竖直向地震时程惯性力
Figure FDA0002560919460000037
6.根据权利要求1或2所述的昔格达地层岩质边坡动力稳定性的时程分析方法,其特征在于:所述地震时程加速度计算具体为:
(1)地震作用下昔格达地层岩质边坡每个时步对应的计算时间按下式计算:
Figure FDA0002560919460000038
式中,tj是第j个时步对应的计算时间,j=(1,…,Nt),Nt是地震振动时长内的时步数量,Ts为地震的振动时长;
(2)昔格达地层岩质边坡条块的水平方向地震时程加速度按下式计算:
Figure FDA0002560919460000039
式中:
Figure FDA00025609194600000310
是第j个时步昔格达地层岩质边坡第i个条块沿水平方向的地震时程加速度;i=(1,…,Ns),Ns是滑体条块的数量;ξ为地震作用效应的折减系数;αH(tj)是第j个时步昔格达地层岩质边坡的坡脚处输入的水平方向地震时程加速度;λm是地震加速度最大分布系数;Hi是第i个条块的形心至边坡的坡脚处的垂直高度;H是昔格达地层岩质边坡的高度;
(3)昔格达地层岩质边坡条块的竖直方向地震时程加速度按下式计算:
Figure FDA0002560919460000041
式中:
Figure FDA0002560919460000042
是第j个时步昔格达地层岩质边坡第i个条块沿竖直方向的地震时程加速度;αV(tj)是第j个时步昔格达地层岩质边坡的坡脚处输入的竖直方向地震时程加速度。
7.根据权利要求1或2所述的昔格达地层岩质边坡动力稳定性的时程分析方法,其特征在于:所述建立昔格达地层岩质边坡动力稳定性的非线性数学规划模型具体为:
(1)建立目标函数
设昔格达地层岩质边坡的安全系数为目标函数,并求解边坡安全系数的最大值,目标函数如下:
Maximize:Kj
式中:Kj是第j个时步昔格达地层岩质边坡动力稳定性的安全系数,j=(1,…,Nt),Nt是地震振动时长内的时步数量;Maximize表示“使最大”;
(2)条块的动力平衡方程:条块受到自重、地震力、底滑面和条间接触面上的法向力、剪力的作用而保持平衡,条块的动力平衡方程按以下流程建立:
①建立第j个时步昔格达地层岩质边坡第i个条块的水平方向的动力平衡方程,具体如下:
Figure FDA0002560919460000043
式中:i=(1,…,Ns),Ns是条块的数量;θi第i个条块底滑面的倾角;
Figure FDA0002560919460000044
是昔格达地层岩质边坡第i个条块底滑面的法向力,
Figure FDA0002560919460000045
以受压为正;Fi d是昔格达地层岩质边坡第i个条块底滑面的剪力,Fi d以对条块产生逆时针的转动效果为正;
Figure FDA0002560919460000046
是昔格达地层岩质边坡第i个条块左侧条间接触面的法向力,
Figure FDA0002560919460000047
以受压为正,当i=1时
Figure FDA0002560919460000048
Figure FDA0002560919460000049
是昔格达地层岩质边坡第i个条块右侧条间接触面的法向力,
Figure FDA00025609194600000410
以受压为正,当i=Ns
Figure FDA0002560919460000051
Figure FDA0002560919460000052
是第j个时步昔格达地层岩质边坡第i个条块的水平向地震时程惯性力,
Figure FDA0002560919460000053
mi是第i个条块的质量,
Figure FDA0002560919460000054
是第j个时步昔格达地层岩质边坡第i个条块沿水平方向的地震时程加速度;
②建立第j个时步昔格达地层岩质边坡第i个条块的竖直方向的动力平衡方程,具体如下:
Figure FDA0002560919460000055
式中:Fi L是昔格达地层岩质边坡第i个条块左侧条间接触面的剪力,Fi L以对条块产生逆时针的转动效果为正,当i=1时Fi L=0;Fi R是昔格达地层岩质边坡第i个条块右侧条间接触面的剪力,Fi R以对条块产生逆时针的转动效果为正,当i=Ns时Fi R=0;Gi是第i个条块的自重,Gi=miγ,γ是昔格达地层岩体的容重;
Figure FDA0002560919460000056
是第j个时步昔格达地层岩质边坡第i个条块的竖直向地震时程惯性力,
Figure FDA0002560919460000057
Figure FDA0002560919460000058
是第j个时步昔格达地层岩质边坡第i个条块沿竖直方向的地震时程加速度;
(3)建立条块底滑面的屈服条件,具体如下:
Figure FDA0002560919460000059
式中:
Figure FDA00025609194600000510
是昔格达地层岩质边坡第i个条块底滑面的长度;
Figure FDA00025609194600000511
是岩体的内摩擦角,cs是岩体的凝聚力;
(4)建立左侧条间接触面和右侧条间接触面的屈服条件,具体如下:
Figure FDA00025609194600000512
式中:
Figure FDA00025609194600000513
是昔格达地层岩质边坡第i个条块左侧条间接触面的长度;
Figure FDA00025609194600000514
是昔格达地层岩质边坡第i个条块右侧条间接触面的长度;
(5)建立昔格达地层岩质边坡动力稳定性的非线性数学规划模型,具体为:
Figure FDA0002560919460000061
式中:Si是第i个条块的面积,g为重力加速度;ξ为地震作用效应的折减系数;αH(tj)是第j个时步昔格达地层岩质边坡的坡脚处输入的水平方向地震时程加速度;λm是地震加速度最大分布系数;Hi是第i个条块的形心至边坡的坡脚处的垂直高度;H是昔格达地层岩质边坡的高度;αV(tj)是第j个时步昔格达地层岩质边坡的坡脚处输入的竖直方向地震时程加速度。
CN202010610995.9A 2020-06-29 2020-06-29 一种昔格达地层岩质边坡动力稳定性的时程分析方法 Active CN111651901B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010610995.9A CN111651901B (zh) 2020-06-29 2020-06-29 一种昔格达地层岩质边坡动力稳定性的时程分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010610995.9A CN111651901B (zh) 2020-06-29 2020-06-29 一种昔格达地层岩质边坡动力稳定性的时程分析方法

Publications (2)

Publication Number Publication Date
CN111651901A true CN111651901A (zh) 2020-09-11
CN111651901B CN111651901B (zh) 2023-03-07

Family

ID=72346202

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010610995.9A Active CN111651901B (zh) 2020-06-29 2020-06-29 一种昔格达地层岩质边坡动力稳定性的时程分析方法

Country Status (1)

Country Link
CN (1) CN111651901B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113235671A (zh) * 2021-05-07 2021-08-10 中国建筑西南勘察设计研究院有限公司 一种岩质边坡地基的承载力载荷试验系统及方法

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4404320B1 (ja) * 2008-11-05 2010-01-27 有限会社秋山調査設計 斜面の表土変位量算出方法及び防災情報システム
CN103485353A (zh) * 2013-09-24 2014-01-01 昆明理工大学 基于全局最优化的边坡稳定性分析条分法
CN105160093A (zh) * 2015-08-28 2015-12-16 重庆交通大学 强地震荷载下含裂隙岩质边坡沿深度变化的受力分析方法
CN106709182A (zh) * 2016-12-26 2017-05-24 华中科技大学 一种地震作用下顺层岩质边坡稳定可靠性安全评价方法
CN107341285A (zh) * 2017-05-26 2017-11-10 昆明理工大学 一种地震力作用下节理岩质边坡倾倒破坏极限承载力分析下限法
CN107630474A (zh) * 2017-08-17 2018-01-26 青岛理工大学 一种边坡抗滑桩抗滑耐久性的动力测定方法
CN109214054A (zh) * 2018-08-06 2019-01-15 重庆交通大学 基于位移的边坡动力稳定性分析方法
US20190177932A1 (en) * 2016-12-02 2019-06-13 Dalian University Of Technology Combined energy dissipation scaffolding structure for preventing falling rock for high and steep slope in seismic region
CN109914379A (zh) * 2019-04-03 2019-06-21 昆明理工大学 一种基于条间法向力分布特征的边坡稳定性极限平衡计算方法
US20190250291A1 (en) * 2018-02-09 2019-08-15 China University Of Geosciences, Beijing Method and system for acquiring probability of slope failure and destabilization caused by earthquake
CN110688747A (zh) * 2019-09-19 2020-01-14 昆明理工大学 一种基于安全性和经济性多目标优化的岩质边坡锚固计算方法
CN111324942A (zh) * 2019-12-27 2020-06-23 昆明理工大学 一种考虑滑面动力渐进破坏的地震边坡稳定性分析方法

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4404320B1 (ja) * 2008-11-05 2010-01-27 有限会社秋山調査設計 斜面の表土変位量算出方法及び防災情報システム
CN103485353A (zh) * 2013-09-24 2014-01-01 昆明理工大学 基于全局最优化的边坡稳定性分析条分法
CN105160093A (zh) * 2015-08-28 2015-12-16 重庆交通大学 强地震荷载下含裂隙岩质边坡沿深度变化的受力分析方法
US20190177932A1 (en) * 2016-12-02 2019-06-13 Dalian University Of Technology Combined energy dissipation scaffolding structure for preventing falling rock for high and steep slope in seismic region
CN106709182A (zh) * 2016-12-26 2017-05-24 华中科技大学 一种地震作用下顺层岩质边坡稳定可靠性安全评价方法
CN107341285A (zh) * 2017-05-26 2017-11-10 昆明理工大学 一种地震力作用下节理岩质边坡倾倒破坏极限承载力分析下限法
CN107630474A (zh) * 2017-08-17 2018-01-26 青岛理工大学 一种边坡抗滑桩抗滑耐久性的动力测定方法
US20190250291A1 (en) * 2018-02-09 2019-08-15 China University Of Geosciences, Beijing Method and system for acquiring probability of slope failure and destabilization caused by earthquake
CN109214054A (zh) * 2018-08-06 2019-01-15 重庆交通大学 基于位移的边坡动力稳定性分析方法
CN109914379A (zh) * 2019-04-03 2019-06-21 昆明理工大学 一种基于条间法向力分布特征的边坡稳定性极限平衡计算方法
CN110688747A (zh) * 2019-09-19 2020-01-14 昆明理工大学 一种基于安全性和经济性多目标优化的岩质边坡锚固计算方法
CN111324942A (zh) * 2019-12-27 2020-06-23 昆明理工大学 一种考虑滑面动力渐进破坏的地震边坡稳定性分析方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
CHEN YULONG: "Research on landslide mechanism and stability of the rock slope", 《2011 INTERNATINAL SYMPOSIUM ON WATER RESOURCE AND ENVIRONMENTAL PROTECTION》 *
刘文连等: "昔格达组地方研究中需要注意的若干问题", 《地学前缘》 *
彭勇胜等: "云南某悬崖酒店边坡稳定性分析", 《中国水运(下半月)》 *
彭勇胜等: "云南腾冲县曲石乡某边坡稳定性分析", 《中国水运(下半月)》 *
胡鑫等: "地震作用下滑坡稳定性分析", 《中国优秀硕士学位论文全文数据库 (基础科学辑)》 *
门妮等: "地震作用下边坡动力响应及影响因素研究", 《世界地震工程》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113235671A (zh) * 2021-05-07 2021-08-10 中国建筑西南勘察设计研究院有限公司 一种岩质边坡地基的承载力载荷试验系统及方法
CN113235671B (zh) * 2021-05-07 2022-12-02 中国建筑西南勘察设计研究院有限公司 一种岩质边坡地基的承载力载荷试验系统及方法

Also Published As

Publication number Publication date
CN111651901B (zh) 2023-03-07

Similar Documents

Publication Publication Date Title
Li et al. Microplate tectonics: New insights from micro-blocks in the global oceans, continental margins and deep mantle
Bi et al. The 3D numerical simulation for the propagation process of multiple pre-existing flaws in rock-like materials subjected to biaxial compressive loads
Zhao et al. Compression–shear strength criterion of coal–rock combination model considering interface effect
Pelanti et al. A Roe-type scheme for two-phase shallow granular flows over variable topography
Hopkins et al. On the numerical simulation of the sea ice ridging process
Liu et al. A coupling model of XFEM/peridynamics for 2D dynamic crack propagation and branching problems
CN105138731B (zh) 一种水合物分解引起海底斜坡不稳定性评价系统及方法
Chang et al. Seismic anisotropy of upper mantle in eastern China
Sokoutis et al. Indentation of a continent with a built-in thickness change: experiment and nature
CN108090268B (zh) 一种粘弹性边界下地震时程波的集成化添加方法
Song et al. Multistage structural deformations of a superimposed basin system and its tectonic response to regional geological evolution: A case study from the Late Jurassic-Early Cretaceous Tanan depression, Hailar-Tamtsag basin
CN111651901B (zh) 一种昔格达地层岩质边坡动力稳定性的时程分析方法
Yanpeng et al. Coupled FEM/SBFEM investigation on the characteristic analysis of seismic motions of a trapezoidal canyon in a layered half-space
Dasgupta et al. Controls of subducting slab dip and age on the extensional versus compressional deformation in the overriding plate
Zhang et al. Geometrical characteristic investigation of the Baihetan irregular columnar jointed basalt and corresponding numerical reconstruction method
CN111651902B (zh) 一种昔格达地层岩质边坡稳定性计算的拟动力下限法
CN112100859B (zh) 一种边坡岩土体抗剪强度参数多重反演方法
CN108363115A (zh) 一种致密石英砂岩孔隙度恢复模拟方法
CN104678439A (zh) 一种多波连续介质速度分析方法
WO2023124664A1 (zh) 岩体稳定性极限分析方法
CN109101684B (zh) 一种规则层状半无限体中瑞利波基阶模态频散曲线的解析算法
Nikitin et al. Explicit-implicit schemes for calculating the dynamics of layered media with nonlinear conditions at contact boundaries
Lou et al. Influence of multi-trend major fault reactivation during multiphase rifting: Beier Depression, Hailar Basin, NE China
Sheng et al. Modified Bishop method for stability analysis of weakly sloped subgrade under centrifuge model test
CN111596356A (zh) 一种昔格达组地层岩质边坡的地震惯性力计算方法

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