CN115238486A - 一种隧道顶板失效概率分析方法 - Google Patents

一种隧道顶板失效概率分析方法 Download PDF

Info

Publication number
CN115238486A
CN115238486A CN202210820495.7A CN202210820495A CN115238486A CN 115238486 A CN115238486 A CN 115238486A CN 202210820495 A CN202210820495 A CN 202210820495A CN 115238486 A CN115238486 A CN 115238486A
Authority
CN
China
Prior art keywords
point
tunnel
failure probability
checking
parameters
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.)
Pending
Application number
CN202210820495.7A
Other languages
English (en)
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.)
Changan University
Original Assignee
Changan University
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 Changan University filed Critical Changan University
Priority to CN202210820495.7A priority Critical patent/CN115238486A/zh
Publication of CN115238486A publication Critical patent/CN115238486A/zh
Pending legal-status Critical Current

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/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
    • 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)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提供了一种隧道顶板失效概率分析方法，属于隧道工程领域，包括如下步骤：建立隧道顶板楔形体失稳的安全系数FS；根据安全系数FS确定岩体力学参数；利用主成分分析方法分析岩体力学参数中对隧道顶板失效概率贡献最大的参数，并将其设为初始验算点；根据安全系数FS确的功能函数及其梯度，在初始验算点处按照Taylor级数展开并取至一次项；采用转换后的一次二阶矩的功能函数对贡献最大的参数进行迭代计算，获得满足迭代终止条件的确定的验算点并计算其对应的灵敏度系数；根据确定的验算点及其对应的灵敏度系数计算可靠度指标；根据可靠度指标计算隧道顶板楔形体的失效概率；本发明分析隧道顶板失效点的精度更高。

Description

一种隧道顶板失效概率分析方法
技术领域
本发明属于隧道工程领域，具体涉及一种隧道顶板失效概率分析方法。
背景技术
随着我国交通科技的发展，隧道在交通运输领域发挥着至关重要的作用，然而由于其结构劣化引发的隧道顶板的失稳坍塌、渗水等事故对国民经济以及人员生命安全产生了重大影响。经研究，隧道围岩地层环境及力学特性复杂，导致决定其稳定性的各种因素存在不可忽略的不确定性，利用可靠性分析方法研究这些不确定性对结构稳定性的影响，计算结构失稳风险，能够更为全面地对隧道安全度进行评估。并且在隧道掘进建设过程中隧道顶板处容易形成楔形体，它是横剖面为V字形剖面,两个结构面相交所构成的势能有减小趋势的封闭结构体。楔形体受围岩应力状态、岩体结构及洞室断面形状和尺寸等多种因素影响，变形破坏过程较为复杂，故其稳定性计算也是一个难题。因此，深入开展隧道顶板楔形体稳定性的研究对于研究隧道稳定性是十分必要的。
目前，考虑岩体力学参数随机性常用的研究方法为蒙特卡罗随机抽样法，根据岩体力学参数的先验概率产生N个样本来计算失效概率，原理简单，但缺点是要进行多次样本抽样，比较繁琐、效率较低；基于中心点法的一次二阶矩法，解决了蒙特卡罗方法需要多次重复抽样的问题，利用岩体力学参数随机变量的前两阶矩来计算失效概率，故中心点法计算简便，若分析精度要求不太高，仍有一定的实用价值，但精度较低、误差相对来说较大。
发明内容
为了克服上述现有技术存在的不足，本发明提供了一种隧道顶板失效概率分析方法。
为了实现上述目的，本发明提供如下技术方案：
一种隧道顶板失效概率分析方法，包括如下步骤：
根据极限平衡法确定隧道顶板楔形体的破坏模式，建立隧道顶板楔形体失稳的安全系数FS；
根据所述安全系数FS确定影响隧道顶板楔形体失稳的岩体力学参数；
利用主成分分析方法分析所述岩体力学参数中对隧道顶板失效概率贡献最大的参数，并将所述贡献最大的参数的均值设为初始验算点；
根据所述安全系数FS确定一次二阶矩的功能函数及其梯度，在所述初始验算点处对所述一次二阶矩的功能函数及其梯度按照Taylor级数展开并取至一次项；
采用转换后的一次二阶矩的功能函数及其梯度对所述贡献最大的参数进行迭代计算，获得满足迭代终止条件的确定的验算点并计算其对应的灵敏度系数；
根据所述确定的验算点及其对应的所述灵敏度系数计算可靠度指标；
根据所述可靠度指标计算隧道顶板楔形体的失效概率。
优选的，所述安全系数FS的计算公式为：
Figure BDA0003744135450000021
所述功能函数的计算公式为：
Z＝gX(X)＝FS-1
即，
Figure BDA0003744135450000022
其中，N为法向力，W为楔形体的重量，α为隧道楔形体半顶角，
Figure BDA0003744135450000025
为节理摩擦角。
优选的，在所述初始验算点处按照Talyor级数展开并取至一次项，将所述一次二阶矩的功能函数Z转换为以下形式：
Figure BDA0003744135450000023
其中，gX(x*)＝0，gX(x*)＝0表示点x*为极限状态面上的一点，x*为验算点
Figure BDA0003744135450000024
Figure BDA0003744135450000031
的集合，
Figure BDA0003744135450000032
为该集合中的一个验算点，Xi为岩体力学参数，gX(x*)为功能函数的梯度。
优选的，在所述贡献最大的参数的空间，ZL＝0为过所述初始验算点处的极限状态平面的平面，利用相互独立正态分布随机变量线性组合的性质，ZL的均值和标准差分别为：
Figure BDA0003744135450000033
Figure BDA0003744135450000034
其中，
Figure BDA0003744135450000035
为岩体力学参数的均值，
Figure BDA0003744135450000036
为岩体力学参数Xi的标准差。
优选的，所述岩体力学参数包括节理摩擦角、楔形体的中心高度、楔形体半顶角、侧压力系数、垂直应力、容重、圆形隧道半径、节理的切向和法向刚度比。
优选的，所述岩体力学参数通过现场或工地实验室测得。
优选的，所述确定的验算点的计算过程如下：
①假定所述初始验算点，初始验算点为贡献最大的参数的均值；
②计算灵敏度系数
Figure BDA0003744135450000037
所述灵敏度系数
Figure BDA0003744135450000038
公式为：
Figure BDA0003744135450000039
③计算结构可靠度指标β；
所述结构可靠度指标β公式为：
Figure BDA00037441354500000310
其中，
Figure BDA0003744135450000041
为岩体力学参数的均值，
Figure BDA0003744135450000042
为岩体力学参数Xi的标准差；
Figure BDA0003744135450000043
为验算点；
④计算新的验算点
Figure BDA0003744135450000044
其公式为：
Figure BDA0003744135450000045
⑤若满足终止迭代条件，终止迭代；否则，以新的验算点
Figure BDA0003744135450000046
重复步骤②～④，直至前后两次验算点之差小于允许误差ε，即可得到确定的验算点值。
优选的，所述隧道顶板楔形体失效概率Pf公式为：
Pf＝Φ(-β)＝1-Φ(β)。
优选的，利用所述主成分分析方法分析所述岩体力学参数中对隧道顶板失效概率贡献最大的参数之前，需进行如下步骤：
统计所述岩体力学参数的均值和标准差；
根据所述岩体力学参数的均值和标准差产生服从正态分布的随机场；
根据所述随机场确定所述岩体力学参数主成分分析的样本量；
再利用所述主成分分析方法分析所述样本量，得到所述贡献最大的参数。
优选的，所述主成分分析方法所研究的指标为主成分贡献率和初始特征值，所述主成分贡献率的计算公式为：
Figure BDA0003744135450000047
其中，p表示主成分，λ表示主成分的方差。
本发明提供的一种隧道顶板失效概率分析方法具有以下有益效果：
本发明通过隧道顶板楔形体的破坏模式可以确定隧道顶板楔形体的安全系数；通过安全系数可以确定影响隧道顶板楔形体失效的岩体力学参数；通过主成分分析法可以计算出不同岩体力学参数中影响隧道顶板楔形体稳定的权重，找到影响最大的参数；通过基于验算点法的一次二阶矩方法，将功能函数的线性化Taylor展开点选在失效面上，并对其进行不断迭代，从而找到更精确的失效点。本发明通过基于验算点法的一次二阶矩方法，将计算机技术引入到岩土领域中，能够更好的应用于计算失效概率上，计算效率和准确性更高，计算误差较小。
附图说明
为了更清楚地说明本发明实施例及其设计方案，下面将对本实施例所需的附图作简单地介绍。下面描述中的附图仅仅是本发明的部分实施例，对于本领域普通技术人员来说，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的隧道顶板失效概率分析方法流程图；
图2为本发明实施例提供的隧道顶板楔形体破坏模式示意图；
图3为本发明实施例提供的隧道顶板失效概率分析方法与蒙特卡罗随机抽样法(MCS)的隧道顶板失效概率方法的对比验证图。
具体实施方式
为了使本领域技术人员更好的理解本发明的技术方案并能予以实施，下面结合附图和具体实施例对本发明进行详细说明。以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本发明的保护范围。
在本发明的描述中，需要理解的是，术语“中心”、“纵向”、“横向”、“长度”、“宽度”、“厚度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”、“轴向”、“径向”、“周向”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明的技术方案和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。
此外，术语“第一”、“第二”等仅用于描述目的，而不能理解为指示或暗示相对重要性。在本发明的描述中，需要说明的是，除非另有明确的规定或限定，术语“相连”、“连接”应作广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体式连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以是通过中间媒介间接相连。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本发明中的具体含义。在本发明的描述中，除非另有说明，“多个”的含义是两个或两个以上，在此不再详述。
实施例
本发明提供了一种隧道顶板失效概率分析方法具体如图1-3所示，包括如下步骤：
步骤1：根据极限平衡法确定隧道顶板楔形体的破坏模式，建立隧道顶板楔形体失稳的安全系数FS；
安全系数FS的计算公式为：
Figure BDA0003744135450000061
安全系数FS由其物理意义FS＝2Scosα/(2Nsinα+W)推导得出。
其中，N为法向力，W为楔形体的重量，α为隧道楔形体半顶角，
Figure BDA0003744135450000066
为节理摩擦角。
安全系数FS中的N和W经由以下公式推导而来：
Figure BDA0003744135450000062
Figure BDA0003744135450000063
Figure BDA0003744135450000064
Figure BDA0003744135450000065
Figure BDA0003744135450000071
Figure BDA0003744135450000072
Figure BDA0003744135450000073
其中，S为切向力，H0为楔形体内部的水平力，h为楔形体的中心高度，K0为侧压力系数，p为垂直应力，γ为容重，ks/kn为节理的切向和法向刚度比，β为结构可靠度指标。
步骤2：根据安全系数FS确定影响隧道顶板楔形体稳定性的岩土力学参数，如图2所示，岩土力学参数包括节理摩擦角
Figure BDA0003744135450000075
楔形体的中心高度h、楔形体半顶角α、侧压力系数K0、垂直应力p、容重γ、节理的切向和法向刚度比ks/kn，再通过现场或工地实验室获取岩体力学参数并统计岩体力学参数的均值和标准差。
步骤3：利用主成分分析方法分析岩体力学参数中对隧道顶板失效概率贡献最大的参数，并将贡献最大的参数的均值设为初始验算点；
在利用主成分分析方法分析岩体力学参数中对隧道顶板失效概率贡献最大的参数之前，需进行如下步骤：
(1)统计岩体力学参数的均值和标准差；
岩土力学参数的取值范围如表1所示：
表1不同岩土力学参数的取值范围
Figure BDA0003744135450000074
(2)根据岩体力学参数的均值和标准差产生服从正态分布的随机场,根据随机场确定岩体力学参数主成分分析的样本量；假设不同的岩土力学参数服从如X(μ,σ)形式的正态分布，借助MATLAB软件生成岩体力学参数的随机场，之后将X～N(μ,σ)随机场的岩体力学参数代入安全系数公式FS，统计FS＜1的样本量n。具体不同样本N下隧道顶板楔形体失效概率见表2，在图3中可以清晰地看出样本量N变化下Pf的变化。当N＝10时，n＝1，Pf＝1/10＝0.1；当N＝100时，n＝16，Pf＝16/100＝0.16；当N＝1000时，n＝155，Pf＝155/1000＝0.155；当N＝10000时，n＝1558，Pf＝1558/10000＝0.1558；当N＝100000时，n＝15590，Pf
＝15590/100000＝0.1559。可以看出当N较小时，Pf会有一定的波动，当N为较大之时，Pf会收敛于固定的值，本发明的实施例中，Pf＝0.155。
表2不同样本量对应的隧道顶板失效概率
Figure BDA0003744135450000081
Figure BDA0003744135450000091
根据表2可知，在本实施例中样本量N为1000时，Pf基本稳定了。
(3)再利用主成分分析方法分析样本量，得到贡献最大的参数。
在进行主成分分析法之前，需要对原始数据资料进行标准化处理，以此来消除单位不同带来的不可公度性，然后以N＝1000时的随机场作为进行主成分分析的原始数据资料。
主成分分析方法所研究的指标为主成分贡献率和初始特征值，主成分贡献率的计算公式为：
Figure BDA0003744135450000092
其中，p表示主成分，λ表示主成分的方差。
然后通过SPSS软件计算主成分贡献率、初始特征值。
最后通过表3的总方差解释对主成分贡献率和初始特征值进行总结和分析，在本实施例中，可以发现节理摩擦角
Figure BDA0003744135450000094
的特征值最大为1.164，因此，将节理摩擦角
Figure BDA0003744135450000095
作为最主要的影响成分，其余为次要因素。
表3总方差解释
Figure BDA0003744135450000093
步骤4：根据安全系数FS确定一次二阶矩的功能函数及其梯度，在初始验算点处对一次二阶矩的功能函数及其梯度按照Taylor级数展开并取至一次项；
功能函数的计算公式为：
Z＝gX(X)＝FS-1，
即，
Figure BDA0003744135450000101
其中，W为楔形体的重量，α为隧道楔形体半顶角，
Figure BDA00037441354500001012
为节理摩擦角；
将功能函数Z转换为以下形式：
Figure BDA0003744135450000102
其中，gX(x*)＝0，gX(x*)＝0表示点x*为极限状态面上的一点，x*为验算点
Figure BDA0003744135450000103
Figure BDA0003744135450000104
的集合，
Figure BDA0003744135450000105
为该集合中的一个验算点，Xi为岩体力学参数，gX(x*)为功能函数的梯度；
在贡献最大的参数的空间，ZL＝0为过初始验算点处的极限状态平面的平面，利用相互独立正态分布随机变量线性组合的性质，ZL的均值和标准差分别为：
Figure BDA0003744135450000106
Figure BDA0003744135450000107
其中，
Figure BDA0003744135450000108
为岩体力学参数的均值，
Figure BDA0003744135450000109
为岩体力学参数Xi的标准差。
步骤5：采用转换后的一次二阶矩的功能函数及其梯度对贡献最大的参数进行迭代计算，获得满足迭代终止条件的确定的验算点并计算其对应的灵敏度系数。
确定的验算点的计算过程如下：
①假定初始验算点，初始验算点为贡献最大的参数的均值；
②计算灵敏度系数
Figure BDA00037441354500001010
灵敏度系数
Figure BDA00037441354500001011
公式为：
Figure BDA0003744135450000111
③计算结构可靠度指标β；
结构可靠度指标β公式为：
Figure BDA0003744135450000112
其中，μ为岩体力学参数Xi的均值，σ为岩体力学参数Xi标准差，
Figure BDA0003744135450000113
为验算点；
④计算新的验算点
Figure BDA0003744135450000114
其公式为：
Figure BDA0003744135450000115
⑤若满足终止迭代条件，终止迭代。否则，以新的验算点
Figure BDA0003744135450000116
重复步骤②～④，直至前后两次验算点之差小于允许误差ε，即可得到确定的验算点值。
步骤6：根据步骤②～④计算确定的验算点对应的灵敏度系数及可靠度指标；
步骤7：根据可靠度指标计算隧道顶板楔形体的失效概率Pf，隧道顶板楔形体失效概率Pf公式为：
Pf＝Φ(-β)＝1-Φ(β)。
以上实施例仅为本发明较佳的具体实施方式，本发明的保护范围不限于此，任何熟悉本领域的技术人员在本发明披露的技术范围内，可显而易见地得到的技术方案的简单变化或等效替换，均属于本发明的保护范围。

Claims (10)

1.一种隧道顶板失效概率分析方法，其特征在于，包括如下步骤：
根据极限平衡法确定隧道顶板楔形体的破坏模式，建立隧道顶板楔形体失稳的安全系数FS；
根据所述安全系数FS确定影响隧道顶板楔形体失稳的岩体力学参数；
利用主成分分析方法分析所述岩体力学参数中对隧道顶板失效概率贡献最大的参数，并将所述贡献最大的参数的均值设为初始验算点；
根据所述安全系数FS确定一次二阶矩的功能函数及其梯度，在所述初始验算点处对所述一次二阶矩的功能函数按照Taylor级数展开并取至一次项；
采用转换后的一次二阶矩的功能函数及其梯度对所述贡献最大的参数进行迭代计算，获得满足迭代终止条件的确定的验算点并计算其对应的灵敏度系数；
根据所述确定的验算点及其对应的所述灵敏度系数计算可靠度指标；
根据所述可靠度指标计算隧道顶板楔形体的失效概率。
2.根据权利要求1所述的隧道顶板失效概率分析方法，其特征在于，所述安全系数FS的计算公式为：
Figure FDA0003744135440000011
所述功能函数的计算公式为：
Z＝gX(X)＝FS-1，
即，
Figure FDA0003744135440000012
其中，N为法向力，W为楔形体的重量，α为隧道楔形体半顶角，
Figure FDA0003744135440000013
为节理摩擦角。
3.根据权利要求2所述的隧道顶板失效概率分析方法，其特征在于，在所述初始验算点处按照Talyor级数展开并取至一次项，将所述一次二阶矩的功能函数Z转换为以下形式：
Figure FDA0003744135440000014
其中，gX(x*)＝0，gX(x*)＝0表示点x*为极限状态面上的一点，x*为验算点
Figure FDA0003744135440000021
Figure FDA0003744135440000022
的集合，
Figure FDA0003744135440000023
为该集合中的一个验算点，Xi为岩体力学参数，gX(x*)为功能函数的梯度。
4.根据权利要求3所述的隧道顶板失效概率分析方法，其特征在于，在所述贡献最大的参数的空间，ZL＝0为过所述初始验算点处的极限状态平面的平面，利用相互独立正态分布随机变量线性组合的性质，ZL的均值和标准差分别为：
Figure FDA0003744135440000024
Figure FDA0003744135440000025
其中，
Figure FDA0003744135440000026
为岩体力学参数的均值，
Figure FDA0003744135440000027
为岩体力学参数Xi的标准差。
5.根据权利要求1所述的隧道顶板失效概率分析方法，其特征在于，所述岩体力学参数包括节理摩擦角、楔形体的中心高度、楔形体半顶角、侧压力系数、垂直应力、容重、圆形隧道半径、节理的切向和法向刚度比。
6.根据权利要求5所述的隧道顶板失效概率分析方法，其特征在于，所述岩体力学参数通过现场或工地实验室测得。
7.根据权利要求4所述的隧道顶板失效概率分析方法，其特征在于，所述确定的验算点的计算过程如下：
①假定所述初始验算点，初始验算点为贡献最大的参数的均值；
②计算灵敏度系数
Figure FDA0003744135440000028
所述灵敏度系数
Figure FDA0003744135440000029
公式为：
Figure FDA00037441354400000210
③计算结构可靠度指标β；
所述结构可靠度指标β公式为：
Figure FDA0003744135440000031
其中，
Figure FDA0003744135440000032
为岩体力学参数的均值，
Figure FDA0003744135440000033
为岩体力学参数Xi的标准差，
Figure FDA0003744135440000034
为验算点；
④计算新的验算点
Figure FDA0003744135440000035
其公式为：
Figure FDA0003744135440000036
⑤若满足终止迭代条件，终止迭代；否则，以新的验算点
Figure FDA0003744135440000037
重复步骤②～④，直至前后两次验算点之差小于允许误差ε，即可得到确定的验算点值。
8.根据权利要求7所述的隧道顶板失效概率分析方法，其特征在于，所述隧道顶板楔形体失效概率Pf公式为：
Pf＝Φ(-β)＝1-Φ(β)。
9.根据权利要求1所述的隧道顶板失效概率分析方法，其特征在于，利用所述主成分分析方法分析所述岩体力学参数中对隧道顶板失效概率贡献最大的参数之前，需进行如下步骤：
统计所述岩体力学参数的均值和标准差；
根据所述岩体力学参数的均值和标准差产生服从正态分布的随机场；
根据所述随机场确定所述岩体力学参数主成分分析的样本量；
再利用所述主成分分析方法分析所述样本量，得到所述贡献最大的参数。
10.根据权利要求9所述的隧道顶板失效概率分析方法，其特征在于，所述主成分分析方法所研究的指标为主成分贡献率和初始特征值，所述主成分贡献率的计算公式为：
Figure FDA0003744135440000038
其中，p表示主成分，λ表示主成分的方差。
CN202210820495.7A 2022-07-13 2022-07-13 一种隧道顶板失效概率分析方法 Pending CN115238486A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210820495.7A CN115238486A (zh) 2022-07-13 2022-07-13 一种隧道顶板失效概率分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210820495.7A CN115238486A (zh) 2022-07-13 2022-07-13 一种隧道顶板失效概率分析方法

Publications (1)

Publication Number Publication Date
CN115238486A true CN115238486A (zh) 2022-10-25

Family

ID=83673761

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210820495.7A Pending CN115238486A (zh) 2022-07-13 2022-07-13 一种隧道顶板失效概率分析方法

Country Status (1)

Country Link
CN (1) CN115238486A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115587417A (zh) * 2022-11-09 2023-01-10 长安大学 一种利用向量投影响应面法确定隧道顶板失效概率的方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107239645A (zh) * 2017-08-10 2017-10-10 西南交通大学 落石冲击下无回填土拱形明洞结构概率可靠度设计方法
CN110362905A (zh) * 2019-07-03 2019-10-22 中国科学院武汉岩土力学研究所 一种盾构隧道开挖面失稳概率确定方法
US20200342346A1 (en) * 2019-04-24 2020-10-29 Cisco Technology, Inc. Adaptive threshold selection for sd-wan tunnel failure prediction
CN113987816A (zh) * 2021-11-01 2022-01-28 长安大学 基于蒙特卡罗随机抽样和主成分分析的隧道顶板监测方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107239645A (zh) * 2017-08-10 2017-10-10 西南交通大学 落石冲击下无回填土拱形明洞结构概率可靠度设计方法
US20200342346A1 (en) * 2019-04-24 2020-10-29 Cisco Technology, Inc. Adaptive threshold selection for sd-wan tunnel failure prediction
CN110362905A (zh) * 2019-07-03 2019-10-22 中国科学院武汉岩土力学研究所 一种盾构隧道开挖面失稳概率确定方法
CN113987816A (zh) * 2021-11-01 2022-01-28 长安大学 基于蒙特卡罗随机抽样和主成分分析的隧道顶板监测方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
YAQIONG WANG ET AL.: ""Countermeasures to treat collapse during the construction of road tunnel in fault zone a case study from the Yezhuping Tunnel in south Qinling, China"", 《ENVIRONMENTAL EARTH SCIENCES》 *
吴雪莉等: ""改进的一次二阶矩法在隧道可靠度中的应用"", 《工业建筑》 *
王亚琼等: ""隧道施工技术风险评估及风险控制措施研究"", 《公路》 *
田荣生等: ""公路隧道上覆溶洞水压影响隧道稳定性的数值模拟研究"", 《筑路机械与施工机械化》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115587417A (zh) * 2022-11-09 2023-01-10 长安大学 一种利用向量投影响应面法确定隧道顶板失效概率的方法

Similar Documents

Publication Publication Date Title
CN106120878B (zh) 一种地基极限承载力及容许承载力的载荷试验测定方法
CN105279325B (zh) 考虑整体空间作用的钢管混凝土锥形柱稳定计算方法
CN115238486A (zh) 一种隧道顶板失效概率分析方法
CN111238931B (zh) 基于能量演化的页岩脆性指数评价方法
WO2023093874A1 (zh) 建筑物坍塌风险评估方法、装置、计算机设备及存储介质
Zhou et al. Flexural strength evaluation for cold-formed steel lip-reinforced built-up I-beams
CN104729988A (zh) 一种测试新老混凝土黏结滑移性能的试验装置及试验方法
Kettler et al. Experimental tests on bolted steel angles in compression with varying end supports
CN115165629A (zh) 一种岩石岩爆类型倾向性的评估方法
CN113536646B (zh) 一种单层球壳地震失效荷载计算方法
CN115266348B (zh) 一种基于动静载叠加试验的岩爆倾向性等级判断方法
CN113987816A (zh) 基于蒙特卡罗随机抽样和主成分分析的隧道顶板监测方法
CN113326547B (zh) 一种基于应力补偿的边坡锚杆锚固方法
Johansson et al. Application of reliability-based design methods to underground excavation in rock
Williams et al. Comparison of the measured and predicted crack propagation behaviour of Zr–2.5 Nb pressure tube material
Karagiozova et al. Stress–strain relationship for metal hollow sphere materials as a function of their relative density
CN111709129A (zh) 一种类岩堆体裂隙围岩中隧道开挖安全系数确定方法
Aoh et al. On the improvement of calibration coefficients for hole-drilling integral method: part II—experimental validation of calibration coefficients
Gan et al. The Model of Simply Supported Beam Force Monitoring Based on Inertial Internet of Things
CN115587417A (zh) 一种利用向量投影响应面法确定隧道顶板失效概率的方法
Somodi et al. Applicable Equivalent Bow Imperfections in GMNIA for Eurocode Buckling Curves–in Case of Box Sections
Fahool et al. A Numerical Investigation for In-Situ Measurement of Rock Mass Mechanical Properties with a CCBO Probe and Evaluation of the Method’s Error in Estimating the In-Situ Stresses with the Overcoring Technique
CN117235839B (zh) 基于自平衡应力场分析的井筒水泥环安全载荷计算方法
Shen et al. Margin of bearing capacity of one-way straight mortise-tenon joints based on inspection of apparent damage and rotation angle
CN112380610B (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
RJ01 Rejection of invention patent application after publication

Application publication date: 20221025

RJ01 Rejection of invention patent application after publication