CN106021853B - 一种空间变异特性下边坡稳定可靠度分析软件的开发方法 - Google Patents
一种空间变异特性下边坡稳定可靠度分析软件的开发方法 Download PDFInfo
- Publication number
- CN106021853B CN106021853B CN201610297267.0A CN201610297267A CN106021853B CN 106021853 B CN106021853 B CN 106021853B CN 201610297267 A CN201610297267 A CN 201610297267A CN 106021853 B CN106021853 B CN 106021853B
- Authority
- CN
- China
- Prior art keywords
- slope
- analysis
- stochastic
- value
- unstability
- 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
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
Abstract
本发明属于边坡稳定分析软件开发的技术领域,涉及一种空间变异特性下边坡稳定可靠度分析软件的开发方法,包括以下实现过程:在前处理数据输入的基础上,进行边坡离散化与集成分析、得出给定失稳模式的安全系数,进而进行一定次数的最小安全系数搜索并确定主要的失稳区域,最后通过蒙特卡罗抽样方法实现基于多个主要失稳模式的边坡稳定可靠度分析软件开发。本发明的优点在于:(1)能够充分考虑到土体材料空间变异特性带来的***特性。(2)填补了目前商业软件因空间变异特性带来的***特性的计算空白,具有较强的工程应用价值和商业前景。
Description
技术领域
本发明属于边坡稳定分析软件开发的技术领域,涉及一种空间变异特性下边坡稳定可靠度分析软件的开发方法,特别是一种基于极限平衡方法确定空间变异特性下边坡***的响应,并基于***响应进行边坡稳定可靠度分析的软件开发方法。
背景技术
边坡稳定是指自然边坡或人工边坡保持安全稳定的条件和能力。这两类边坡的岩土体在各种内外因素作用下逐渐发生化,坡体应力状态也随之改变,当滑动力或倾覆力达以至超过抗滑力或抗倾覆力而失去平衡时,即出现变破坏,造成灾害或威胁建筑物安全。
边坡稳定可靠度分析是岩土工程设计领域中非常重要的课题之一,边坡设计或者支护时,需要定量地评估边坡可能的失稳概率以及失稳区域。传统的边坡设计时,需要借助于商业软件进行极限平衡分析或者基于强度折减方法计算边坡稳定的抗滑最小安全系数以及相对应的失稳区域,并基于这个最小安全系数和相应的失稳区域进行设计。另一方面,由于岩土工程材料的不确定性,尤其是在长期的地质作用下,土体材料呈现出明显的空间变异特性,这种空间变异特性导致边坡在失稳时呈现出***特性,也就是说,边坡失稳时有不止一个的主要失稳区域,这种***特性在常规的商业软件中很少涉及。目前,市场上还没有一种能够考虑土体材料空间变异特性的边坡稳定可靠度分析软件。
有鉴于上述分析,开发能够考虑土体材料空间变异特性的边坡稳定可靠度分析软件是非常比要的,同时具有较强的工程应用价值和商业前景。
发明内容
根据以上现有技术的不足,本发明提供一种空间变异特性下边坡稳定可靠度分析软件的开发方法,能够充分考虑到土体材料空间变异特性带来的***特性。
本发明所述的一种空间变异特性下边坡稳定可靠度分析软件的开发方法,其特征在于包括以下实现过程:在前处理数据输入的基础上,进行边坡离散化与集成分析、得出给定失稳模式的安全系数,进而进行一定次数的最小安全系数搜索并确定主要的失稳区域,最后通过蒙特卡罗抽样方法实现基于多个主要失稳模式的边坡稳定可靠度分析软件开发。
具体包括以下实现过程:
前处理分析:
(1)读取用于计算边坡稳定的相关数据,包括读取边坡几何剖面数据、土体材料统计数据、浸润线位置数据或者孔压数据、以及外荷载数据;
(2)假定边坡滑动的失稳模式类型;
(3)根据土体材料统计数据中的波动范围值确定每个边坡土层的虚拟土层个数,并进行边坡几何剖面数据的二次更新;
确定性分析:
(4)将给定失稳模式进行离散化,得到一系列的垂直土条,并结合二次更新后的边坡几何模型,得出每一个垂直土条的相关参数;
(5)给定极限平衡方法的安全系数求解,在得到所有离散垂直土条的相关参数基础上,利用简化Bishop法计算给定失稳模式的安全系数;
(6)变换不同的潜在失稳区域,按照步骤(5)进行安全系数的求解,最终经过比较得出具有最小安全系数的失稳区域;
边坡***响应确定:
(7)对所有的随机变量中的一个进行如下迭代:设置该随机变量的取值为其均值减去1倍标准差值,其余随机变量取为其均值,在确定了随机变量的取值之后,进行步骤(4)、(5)、(6)的计算分析,得出该种工况下的最小安全系数及其对应的失稳区域;
(8)对所有的随机变量进行步骤(7)的迭代计算,最终对应每一个随机变量,均得到一个最小安全系数及其对应的失稳区域;
(9)将步骤(7)和(8)得到的失稳区域进行汇总,将其视为边坡***的主要组成部分,边坡***的整体响应可以由这些离散的失稳区域来代表;
边坡可靠度分析与后处理:
(10)生成符合随机变量分布的蒙特卡罗抽样样本值,在每一个蒙特卡罗抽样样本值下,针对步骤(9)中得到的每一个离散的失稳区域,进行步骤(4)和(5)的计算分析,得出每一个离散的失稳区域的安全系数,比较所有离散的失稳区域对应的安全系数,挑选最小的安全系数作为边坡***在该蒙特卡罗抽样样本值下的响应,同时,可以确定该最小安全系数对应的失稳区域;
(11)在所有样本值对应的边坡响应求出后,可以统计边坡响应小于1的蒙特卡罗样本值个数,将该个数与总的蒙特卡罗样本值个数之比定义为边坡的失效概率,该失效概率值可以代表边坡的可靠度程度;此外,每个失稳区域出现最小安全系数的次数及其相应的几何位置可以显示在边坡剖面图上,以便更直观地确定边坡***响应的主要来源。
其中,优选方案如下:
步骤(1)中所述的边坡剖面数据包括坡高、坡角、土层数目以及坡脚与坡顶延伸范围的范围;土体材料统计数据包括容重、粘聚力、内摩擦角的均值与标准差以及表征土层参数空间变异特性程度的波动范围值。
步骤(2)中所述的边坡滑动的失稳模式包括圆弧失稳模式和非圆弧失稳模式。
步骤(4)中所述的垂直土条的相关参数包括每个垂直土条的条底坐标、条底与地面相对应的坐标,以及垂直土条的自重、条底的内摩擦角、粘聚力以及条底孔压数据。
步骤(5)所述的简化Bishop法的公式为:
步骤(10)中利用Cholesky分解方法生成符合随机变量分布的蒙特卡罗抽样样本值。
本发明的优点在于:(1)充分考虑了土体参数的空间变异特性,在此基础上进行了边坡几何剖面数据的二次更新,并基于传统的极限平衡方法进行确定性分析,通过对随机变量的逐个迭代分析,最终获得边坡***的响应源,利用这些响应源进行边坡***可靠度分析,并在后处理中显示边坡***的失效概率和各个响应源对边坡***的贡献程度。(2)填补了目前商业软件因空间变异特性带来的***特性的计算空白,具有较强的工程应用价值和商业前景。
附图说明
图1为本发明的流程框图。
具体实施方式
以下结合附图和实施例对本发明作进一步说明。
实施例1:
一种空间变异特性下边坡稳定可靠度分析软件的开发方法,包括以下实现过程:
前处理分析:
(1)读取用于计算边坡稳定的相关数据,包括读取边坡几何剖面数据、土体材料统计数据、浸润线位置数据或者孔压数据、以及外荷载数据。其中,边坡剖面数据包括坡高、坡角、土层数目以及坡脚与坡顶延伸范围的范围;土体材料统计数据包括容重、粘聚力、内摩擦角的均值与标准差以及表征土层参数空间变异特性程度的波动范围值。此处需要开设公共数组来保存这些数据。
(2)假定边坡滑动的失稳模式类型,譬如圆弧失稳模式或者非圆弧失稳模式。每种失稳模式的数学描述需要的变量不同,当假定为圆弧失稳模式时,需要圆弧的圆心坐标以及半径等3个变量,此处需要判断圆弧是否与边坡地面有交点,若有,则需保存其与边坡地面的两个交点作为失稳区域的滑入点和滑出点;若没有交点,则需要剔除该失稳模式。
(3)根据土体材料统计数据中的波动范围值确定每个边坡土层的虚拟土层个数,并进行边坡几何剖面数据的二次更新,此处需要开设新的数组来存储二次更新后的边坡几何剖面数据。此处需要设定选项由用户来自主决定剖面数据二次更新的程度,即由用户确定一个0.1~0.2之间的系数,此系数与波动范围的乘积作为更新的依据。
确定性分析:
(4)由用户选择垂直土条的间距,自滑出点开始,至滑入点结束,将给定失稳模式进行离散化,得到一系列的垂直土条,并结合二次更新后的边坡几何模型。此处需要开设数组记录每个垂直土条的条底坐标、条底与地面相对应的坐标,以及垂直土条的自重、条底的内摩擦角、粘聚力以及条底孔压数据。
(5)利用极限平衡方法中的简化Bishop法计算离散后的失稳区域的稳定安全系数,此处需要编制子程序进行迭代计算。
(6)按照一定的变化法则,譬如网格法、复合形法、和声搜索算法来变换不同的潜在失稳区域,调用步骤(5)中的子程序计算不同失稳模式的安全系数,并开设数组进行保存,经过比较或者优化确定最小的安全系数值。
边坡***响应确定:
(7)对所有的随机变量中的一个进行如下迭代:设置该随机变量的取值为均值减去1倍标准差值,其余随机变量取为其均值,在确定了随机变量的取值之后,将这些取值视为确定值并进行步骤(4)、(5)、(6)的计算,直至最终确定出最小的安全系数及其对应的失稳区域,此处需要开设数组记录失稳区域的数学描述变量。
(8)对所有的随机变量进行步骤(7)的迭代计算,最终对应每一个随机变量,均得到一个最小安全系数及其对应的失稳区域,并保存在开设的数组中。
(9)将保存在数组中的失稳区域作为边坡***响应的代表,也就是说边坡在不同的土层参数下,其最小安全系数值可以由这些失稳区域的试算来确定。
边坡可靠度分析与后处理:
(10)利用Cholesky分解方法生成随机变量的抽样值,此处需由用户确定蒙特卡罗抽样次数,针对每一个抽样值,将其视为边坡输入的确定值,确定保存在数组中的失稳区域对应的安全系数值,挑选最小值作为边坡***在该抽样值下的响应,并判断该响应值是否小于1,若是,则称该样本为失效样本,该响应值对应的失稳区域称之为显著区域;若否,则为非失效样本。
(11)统计所有抽样样本中失效样本的个数,该个数与总的抽样样本数目的比值称之为边坡***的失效概率。数组中保存的失稳区域作为显著区域的次数可以直观地判断边坡***的响应源,此处需要开设数组记录每个失稳区域作为显著区域的次数,并最终在边坡剖面上显示每个失稳区域的位置和相对应的显著区域次数,以便为用户显示边坡***的响应源。
Claims (4)
1.一种空间变异特性下边坡稳定可靠度分析软件的开发方法,其特征在于包括以下实现过程:
前处理分析:
(1)读取用于计算边坡稳定的相关数据,包括读取边坡几何剖面数据、土体材料统计数据、浸润线位置数据或者孔压数据、以及外荷载数据;其中,边坡剖面数据包括坡高、坡角、土层数目以及坡脚与坡顶延伸范围的范围;土体材料统计数据包括容重、粘聚力、内摩擦角的均值与标准差以及表征土层参数空间变异特性程度的波动范围值;
(2)假定边坡滑动的失稳模式类型;
(3)根据土体材料统计数据中的波动范围值确定每个边坡土层的虚拟土层个数,并进行边坡几何剖面数据的二次更新;
确定性分析:
(4)将给定失稳模式进行离散化,得到一系列的垂直土条,并结合二次更新后的边坡几何模型,得出每一个垂直土条的相关参数;其中,垂直土条的相关参数包括每个垂直土条的条底坐标、条底与地面相对应的坐标,以及垂直土条的自重、条底的内摩擦角、粘聚力以及条底孔压数据;
(5)给定极限平衡方法的安全系数求解,在得到所有离散垂直土条的相关参数基础上,利用简化Bishop法计算给定失稳模式的安全系数;
(6)变换不同的潜在失稳区域,按照步骤(5)进行安全系数的求解,最终经过比较得出具有最小安全系数的失稳区域;
边坡***响应确定:
(7)对所有的随机变量中的一个进行如下迭代:设置该随机变量的取值为其均值减去1倍标准差值,其余随机变量取为其均值,在确定了随机变量的取值之后,进行步骤(4)、(5)、(6)的计算分析,得出该种工况下的最小安全系数及其对应的失稳区域;
(8)对所有的随机变量进行步骤(7)的迭代计算,最终对应每一个随机变量,均得到一个最小安全系数及其对应的失稳区域;
(9)将步骤(7)和(8)得到的失稳区域进行汇总,将其视为边坡***的主要组成部分,边坡***的整体响应由这些离散的失稳区域来代表;
边坡可靠度分析与后处理:
(10)生成符合随机变量分布的蒙特卡罗抽样样本值,在每一个蒙特卡罗抽样样本值下,针对步骤(9)中得到的每一个离散的失稳区域,进行步骤(4)和(5)的计算分析,得出每一个离散的失稳区域的安全系数,比较所有离散的失稳区域对应的安全系数,挑选最小的安全系数作为边坡***在该蒙特卡罗抽样样本值下的响应,同时,确定该最小安全系数对应的失稳区域;
(11)在所有样本值对应的边坡响应求出后,统计边坡响应小于1的蒙特卡罗样本值个数,将该个数与总的蒙特卡罗样本值个数之比定义为边坡的失效概率,该失效概率值代表边坡的可靠度程度;此外,每个失稳区域出现最小安全系数的次数及其相应的几何位置显示在边坡剖面图上,以便更直观地确定边坡***响应的主要来源。
2.根据权利要求1所述的一种空间变异特性下边坡稳定可靠度分析软件的开发方法,其特征在于步骤(2)中所述的边坡滑动的失稳模式包括圆弧失稳模式和非圆弧失稳模式。
3.根据权利要求1所述的一种空间变异特性下边坡稳定可靠度分析软件的开发方法,其特征在于步骤(5)所述的简化Bishop法的公式为:
其中,FS为安全系数,第i垂直土条宽度bi,底面倾角αi,自重Wi,条底孔压ui,内摩擦角粘聚力ci。
4.根据权利要求1所述的一种空间变异特性下边坡稳定可靠度分析软件的开发方法,其特征在于步骤(10)中利用Cholesky分解方法生成符合随机变量分布的蒙特卡罗抽样样本值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610297267.0A CN106021853B (zh) | 2016-05-09 | 2016-05-09 | 一种空间变异特性下边坡稳定可靠度分析软件的开发方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610297267.0A CN106021853B (zh) | 2016-05-09 | 2016-05-09 | 一种空间变异特性下边坡稳定可靠度分析软件的开发方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106021853A CN106021853A (zh) | 2016-10-12 |
CN106021853B true CN106021853B (zh) | 2019-02-12 |
Family
ID=57081821
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610297267.0A Active CN106021853B (zh) | 2016-05-09 | 2016-05-09 | 一种空间变异特性下边坡稳定可靠度分析软件的开发方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106021853B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106709182B (zh) * | 2016-12-26 | 2020-02-14 | 华中科技大学 | 一种地震作用下顺层岩质边坡稳定可靠性安全评价方法 |
CN107239589A (zh) * | 2017-04-07 | 2017-10-10 | 西安理工大学 | 基于mrvm‑afosm的边坡可靠度分析方法 |
CN108491575A (zh) * | 2018-02-11 | 2018-09-04 | 华北水利水电大学 | 一种基于数字地形进行水库涉水边坡安全系数的计算方法 |
CN108629111B (zh) * | 2018-05-02 | 2019-09-10 | 四川大学 | 一种混凝土重力坝材料参数空间变异性的模拟方法 |
CN109063285B (zh) * | 2018-07-18 | 2022-12-02 | 南昌大学 | 一种土坡钻孔布置方案设计方法 |
CN109446616A (zh) * | 2018-10-18 | 2019-03-08 | 东北电力大学 | 一种均质边坡可靠性评价方法 |
CN109614708B (zh) * | 2018-12-12 | 2019-08-09 | 青岛理工大学 | 边坡上建筑物基础安全距离以及基础埋深的设计方法 |
CN109977554B (zh) * | 2019-03-28 | 2020-06-30 | 青岛理工大学 | 一种边坡滑动面积的评估方法 |
CN110245429B (zh) * | 2019-06-18 | 2020-09-04 | 贵州正业工程技术投资有限公司 | 基于简化Bishop法的环形凸坡稳定性评价方法 |
CN111444649B (zh) * | 2020-03-24 | 2022-10-18 | 成都理工大学 | 基于强度折减法的边坡***可靠度分析方法 |
CN113239435B (zh) * | 2021-05-11 | 2022-09-20 | 青岛理工大学 | 一种水库最优放水速度确定的方法 |
CN113449429A (zh) * | 2021-07-09 | 2021-09-28 | 中国电建集团贵阳勘测设计研究院有限公司 | 一种基于局部平均的边坡稳定性评价修正方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103134433A (zh) * | 2013-02-02 | 2013-06-05 | 青岛理工大学 | 一种利用位移监测鉴别边坡失稳致滑因子的方法 |
CN104899380A (zh) * | 2015-06-11 | 2015-09-09 | 武汉大学 | 一种基于蒙特卡洛模拟的边坡稳定可靠度敏感性分析方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07200534A (ja) * | 1993-12-28 | 1995-08-04 | Toshiba Corp | システム仕様決定装置 |
-
2016
- 2016-05-09 CN CN201610297267.0A patent/CN106021853B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103134433A (zh) * | 2013-02-02 | 2013-06-05 | 青岛理工大学 | 一种利用位移监测鉴别边坡失稳致滑因子的方法 |
CN104899380A (zh) * | 2015-06-11 | 2015-09-09 | 武汉大学 | 一种基于蒙特卡洛模拟的边坡稳定可靠度敏感性分析方法 |
Non-Patent Citations (3)
Title |
---|
Rosenblueth法在边坡可靠度分析中的应用;李亮等;《水利水电科技进展》;20120630;第32卷(第3期);第53-55页 |
基于Bishop法假定的边坡临界滑动场方法及应用;李亮等;《中南大学学报(自然科学版)》;20110930;第42卷(第9期);第2848-2852页 |
边坡临界滑动面与临界可靠度滑动面对比分析;李亮等;《华中科技大学学报(自然科学版)》;20120930;第40卷(第9期);第81-84页 |
Also Published As
Publication number | Publication date |
---|---|
CN106021853A (zh) | 2016-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106021853B (zh) | 一种空间变异特性下边坡稳定可靠度分析软件的开发方法 | |
von Haden et al. | Soils' dirty little secret: Depth‐based comparisons can be inadequate for quantifying changes in soil organic carbon and other mineral soil properties | |
Tricarico | The near-Earth asteroid population from two decades of observations | |
Li et al. | Efficient and consistent reliability analysis of soil slope stability using both limit equilibrium analysis and finite element analysis | |
Van Werkhoven et al. | Sensitivity-guided reduction of parametric dimensionality for multi-objective calibration of watershed models | |
Khoi et al. | The responses of hydrological processes and sediment yield to land‐use and climate change in the Be River Catchment, Vietnam | |
US11268353B2 (en) | Energy deposit discovery system and method | |
Cheng | Seismic lateral earth pressure coefficients for c–φ soils by slip line method | |
Reeves et al. | Extent of coterminous US rangelands: quantifying implications of differing agency perspectives | |
Montgomery Jr et al. | Dipole polarizabilities for a hydrogen atom confined in a penetrable sphere | |
Lloyd et al. | Archaeology and geostatistics | |
Paterson et al. | Spatial variability of Australian soil texture: A multiscale analysis | |
US9897714B2 (en) | Use of monte carlo simulations for establishing correlations and their uncertainties | |
CN113360983B (zh) | 一种边坡可靠度分析与风险评估方法 | |
CN109800459A (zh) | 重力式挡土墙设计方法及装置 | |
Maschberger et al. | The relation between accretion rates and the initial mass function in hydrodynamical simulations of star formation | |
Holzmann et al. | Orbital-dependent backflow wave functions for real-space quantum Monte Carlo | |
Ma et al. | Decomposition of small-footprint full waveform LiDAR data based on generalized Gaussian model and grouping LM optimization | |
Mahmoudi et al. | Stochastic field simulation of slope stability problems: Improvement and reduction of computational effort | |
Rochman et al. | From flatness to steepness: Updating TALYS covariances with experimental information | |
Rahman et al. | Applications of random finite element method in bearing capacity problems | |
Makdisi et al. | Improved computational methods for probabilistic liquefaction hazard analysis | |
Naderi et al. | Application of genetic programming for estimation of soil compaction parameters | |
Nalewajski | Kinetic energy as functional of the correlation hole | |
Sun et al. | Study on displacement back analysis of spatial variability of soil slope |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |