[发明专利]高应力地下洞室流变性软岩稳定性分析的解析计算方法有效
申请号: | 201810917169.1 | 申请日: | 2018-08-13 |
公开(公告)号: | CN109165434B | 公开(公告)日: | 2020-07-28 |
发明(设计)人: | 崔岚;潘雨;盛谦;沈强;崔博 | 申请(专利权)人: | 中国科学院武汉岩土力学研究所 |
主分类号: | G06F30/13 | 分类号: | G06F30/13;G06F119/14 |
代理公司: | 武汉宇晨专利事务所 42001 | 代理人: | 余晓雪;王敏锋 |
地址: | 430071 湖北*** | 国省代码: | 湖北;42 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 应力 地下 流变 性软岩 稳定性 分析 解析 计算方法 | ||
1.一种高应力地下洞室流变性软岩稳定性分析的解析计算方法,其特征在于:所述高应力地下洞室流变性软岩稳定性分析的解析计算方法包括以下步骤:
1)通过现场勘探、测试与分析,获得施工现场软岩的初始地应力、支护压力、隧道半径、软岩流变参量、软岩强度参量、剪胀系数以及体积模量;
2)根据步骤1)中获得的软岩强度参量以及初始地应力计算弹性与塑性区域交界处径向应力;
3)根据步骤1)中获取得到的隧道半径、软岩强度参量、支护压力以及步骤2)中计算得到的弹性与塑性区域交界处径向应力确定围岩塑性区域半径;
4)根据步骤1)中获取得到的初始地应力、支护压力、隧道半径、软岩流变参量、软岩强度参量、剪胀系数、体积模量、步骤2)计算得到的弹性与塑性区域交界处径向应力以及步骤3)确定的围岩塑性区域半径计算得到洞壁处围岩位移;其具体实现方式是:
计算变量A、B、C的具体表达式是:
其中:
u0是洞壁处围岩位移;
R0是隧道半径,单位m;
Kψ为软岩的剪胀系数,单位是无量纲
ψ是软岩的剪胀角,单位是°;
Rp是围岩塑性区域半径;
Y与α分别是两种软岩强度参量,单位是无量纲;
pi是洞壁处软岩受到的支护压力,单位是MPa;
G(t)是软岩流变参量,单位是MPa;
t是时间,单位是天;
K是软岩体积模量,单位MPa;
σ0是软岩所处位置的初始地应力,单位是MPa;
σr,cri为弹性与塑性区域交界处径向应力;
5)针对不同支护压力和软岩流变过程中的不同时刻,根据步骤3)以及步骤4)分别计算围岩塑性区域范围与洞壁处软岩位移,明确高应力洞室流变性软岩稳定性,进行相应工程支护设计。
2.根据权利要求1所述的高应力地下洞室流变性软岩稳定性分析的解析计算方法,其特征在于:所述软岩强度参量的具体获取方式是:
其中:
Y与α分别是两种软岩强度参量,单位是无量纲;
c是软岩粘聚力,单位是MPa;
是软岩摩擦角,单位是°;
所述软岩流变参量的具体获取方式是:
其中:
Gk是软岩Kelvin体剪切模量,单位是MPa;
ηK是软岩Kelvin体黏滞系数,单位是MPa/d;
GM是软岩Maxwell体剪切模量,单位是MPa;
ηM是软岩Maxwell体黏滞系数,单位是MPa/d;
G(t)是软岩流变参量,单位是MPa;
t是时间,单位是天。
3.根据权利要求2所述的高应力地下洞室流变性软岩稳定性分析的解析计算方法,其特征在于:所述步骤2)的具体实现方式是:
其中:
σr,cri为弹性与塑性区域交界处径向应力;
σ0是软岩所处位置的初始地应力,单位是MPa;
Y与α分别是两种软岩强度参量。
4.根据权利要求3所述的高应力地下洞室流变性软岩稳定性分析的解析计算方法,其特征在于:所述步骤3)的具体实现方式是:
其中:
Rp是围岩塑性区域半径,R0是隧道半径,单位是m;
Y与α分别是两种软岩强度参量,单位是无量纲;
pi是洞壁处软岩受到的支护压力,单位是MPa;
σr,cri为弹性与塑性区域交界处径向应力。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国科学院武汉岩土力学研究所,未经中国科学院武汉岩土力学研究所许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201810917169.1/1.html,转载请声明来源钻瓜专利网。