一种深埋硐室围岩压力的计算方法技术

本发明专利技术公开了一种深埋硐室围岩压力的计算方法，步骤如下：构建破坏机制；根据破坏机制建立运动许可的速度场；由内功率与外功率相等推导并得到破坏载荷的上限解析解；在多个破坏荷载中求得最小值，即获得最优上限解。本发明专利技术提出了“向下滑落的楔形塌落体+以顶角为圆心的转动圆弧体+n个平动三角形体”的破坏机制，然后由内功率与外功率相等推导并得到破坏载荷的上限解析解，在多个破坏荷载中采用最优化方法求得最小值，获得最优上限解。该上限解是围岩压力的一个严格的上限值，可为此类工程结构的支护设计提供理论参考。

A method for calculating the surrounding rock pressure of a deep cavern

The invention discloses a deep buried tunnel surrounding rock pressure calculation method comprises the following steps: constructing failure mechanism; according to the failure mechanism of motion velocity field obtained by the license; internal power and external power is equal and the upper bound analytical solution of failure load; obtain the minimum value in more than one load, namely the optimal upper bound solution. The invention proposes a downward slide wedge to collapse + angle for failure mechanism of a center rotating arc +n translational triangle \, and then the internal power and external power equation is derived and obtained the upper bound analytical solution of the failure load, failure load in multiple optimization method was adopted to achieve the minimum value, the optimal upper bound solution get. The upper bound solution is a strict upper limit value of surrounding rock pressure, which can provide a theoretical reference for the support design of this kind of Engineering structures.

【技术实现步骤摘要】
一种深埋硐室围岩压力的计算方法
本专利技术涉及岩石力学围岩压力领域，特别涉及一种深埋硐室围岩压力的计算方法。
技术介绍
在对地下工程进行分析时，一般需要研究土体的应力-应变关系，即岩土体的本构关系。在加载初期，土体将处于弹性阶段，土体的力学特性可用虎克定律描述。当应力强度增至屈服极限时，土体受力特性将由弹性阶段进入塑流阶段，即应力不变而变形将继续增大。当土体的变形增大到某一临界值时，即达到极限荷载，土体将发生滑移或破坏。极限分析法在计算过程中不需要对约束塑流的全过程进行分析，直接关注土体进入流塑状态的极限荷载。与极限平衡法和滑移线法不同，极限分析法考虑了岩土体的应力-应变关系，可以较好地描述由应力与位移表示的边界条件，从而可以得到相比于其他方法更精确的解。在运用极限分析法分析隧道支护力时，AtkinsonandPotts将圆形隧道在无粘性土中施工时需要提供的支护力限定在上限解和下限解组成的区间内，且模型试验证明了该解的有效性。Davisetal应用下限定理计算了隧道在不排水情况下的稳定系数下限解，然后构建了四种有效的浅埋隧道拱顶坍落模式，此破坏机制是AtkinsonandPotts单块体楔形坍落机制的改进，即提出了用不同角度的参数来控制不同几何形状的多个刚性块体的形状和面积，此处Davisetal的结果比AtkinsonandPotts的结果更为精确。Takemuraetal根据隧道围岩变形规律提出了一种通过角度关系表示的5个刚性块体组成的破坏机制，然后分析了土体强度随深度变化和土体各向异性对隧道极限支护力的影响，并通过具体的实验证实了所得结果的有效性。杨峰和阳军生构建了一个由n个块体构成的坍落机制，应用数值分析方法对由虚功率方程构成的支护力目标函数进行优化，可得到浅埋隧道支护力的上限解。FraldiandGuarracino构建了一种以函数f(x)表示的拱顶坍落破坏模式，根据上限定理和变分法则得到了描述速度间断面f(x)的解析表达式和坍落体的规模，通过数值分析验证了所得结论的正确性。HuangandYang引用此破坏机制分析了深埋隧道拱顶在孔隙水压力作用下的稳定性问题，其中将孔隙水压力所做的功率作为外力功率的一部分然后引入到上限定理中。Huang也对深埋矩形和圆形隧道拱顶在孔隙水压力作用下的稳定性问题进行了系统的分析，并计算得到了上限意义下坍落形状的解析解。尤春安通过构建运动许可的围岩速度场和静力许可的应力场计算了极限状态下条带碹所提供的上下限支护力解。但上述利用极限分析法进行分析的过程中，当前存在考虑围岩压力计算假设过多、推理不够严谨、计算结果精度不够的问题。
技术实现思路
为了解决上述技术问题，本专利技术提供一种计算简单、精度高的深埋硐室围岩压力的计算方法。本专利技术解决上述问题的技术方案是：一种深埋硐室围岩压力的计算方法，步骤如下：步骤一：构建破坏机制；步骤二：根据破坏机制建立运动许可的速度场；步骤三：基于极限分析上限定理，由内功率与外功率相等推导并得到破坏载荷的上限解析解；步骤四：在多个破坏荷载中求得最小值，即获得最优上限解。上述深埋硐室围岩压力的计算方法，所述步骤一中，首先构建单圆弧夹层破坏机制、双圆弧夹层破坏机制、圆弧与对数螺旋线夹层破坏机制，然后根据三种破坏机制的优化结果构建最优破坏机制。上述深埋硐室围岩压力的计算方法，所述步骤一中，单圆弧夹层破坏机制由“向下滑落的楔形塌落体+以顶角为圆心的转动圆弧体+三个平动三角形体”构成，双圆弧夹层破坏机制由“向下滑落的楔形塌落体+以顶角为圆心的转动圆弧体+平动三角形体+以顶角为圆心的转动圆弧体+平动三角形体”构成，圆弧与对数螺旋线夹层破坏机制由“向下滑落的楔形塌落体+以顶角为圆心的转动圆弧体+平动三角形体+以顶角为圆心的转动对数螺旋线体+平动三角形体”构成，最优破坏机制由“向下滑落的楔形塌落体+以顶角为圆心的转动圆弧体+n个平动三角形体”构成。上述深埋硐室围岩压力的计算方法，所述步骤三中，外部功率和内部能量耗散率计算如下：①速度大小及速度间断线长度计算：V1＝V0根据圆弧辐射受剪区速度场的分析得出：根据破坏机制的速度矢量图及几何三角关系得出：各速度间断线的长度递推关系为：GB＝GC1②外部功率计算1)重力做功功率2)支护反力功率e＝Kq3)外部总功率外部总功率等于围岩自重功率与支护反力功率之和，即：③内能耗散功率沿AB间断线的能量耗散：圆弧受剪面BC1和受剪区GBC1的能量耗散：沿间断线CkCk+1(k∈Z,1≤k≤n-1)的能量耗散：沿间断线CnF的能量耗散：沿GCk(k∈Z,2≤k≤n)间断线的能量耗散：内部能量总耗散：④支护反力计算根据虚功率原理，得：Wext＝Dint得出围岩压力的表达式：上述公式中，γ表示土体重度，单位为kg/m3；c表示岩土内聚力，单位为kPa；为岩土内摩擦角，单位为度；l为硐室宽度，单位为m；h为硐室高度，单位为m；q为拱顶或顶板单位面积上的支护反力，单位为kPa；e为边墙或两帮单位面积上的支护反力，单位为kPa；K为侧压力系数；V0、V1、V2、…、Vk、Vk+1为速度场中的速度；α1、α2、α3、α4、…、αn为破坏机制中对应三角形的角度，单位为度；W为外部功率，Wexit为外部总功率，Wsoil为围岩自重功率，WT为支护反力功率，单位为瓦；W'、W”、Wk、Wn为变量；D为内能耗散功率，DiDi+1(i＝1,…,n)、分别为对应间断线的能量耗散率，单位为瓦；g'、g”、g”'；f'、f”；g0、g1、g2、…、gk、…gn；f0、f1、f2、…、fk、…fn；l1、l2、…、lk、…ln为变量。本专利技术的有益效果在于：1、本专利技术首先构建一个破坏机制，根据破坏机制建立运动许可的速度场，然后基于极限分析上限定理，由内功率与外功率相等推导并得到破坏载荷的上限解析解，最后在多个破坏荷载中求得最小值，即获得最优上限解，该上限解是围岩压力的一个严格的上限值，可为此类工程结构的支护设计提供理论参考。2、本专利技术推导了深埋硐室围岩压力的计算公式，讨论了硐室断面尺寸、侧压力系数、初始粘聚力与轴向拉应力等因素对硐室拱顶和边墙围岩压力的影响规律，并结合数值模拟与工程实例监测数据验证了公式的合理性，所构建的速度场与数值模拟方法得到的速度场几乎一致，并且本文基于构建深埋硐室的破坏机制，采用极限分析法所得到的围岩压力与普氏理论、数值模拟方法计算所得到的围岩压力以及现场监测结果比较相近，精确度更高。附图说明图1为本专利技术的流程图。图2为本专利技术单圆弧夹层破坏机制示意图。图3为本专利技术双圆弧夹层破坏机制示意图。图4为本专利技术圆弧与对数螺旋夹层破坏机制示意图。图5为转动圆弧与n个平动三角形破坏机制示意图。图6为速度矢量图。图7为围岩压力理论计算值与现场监测值对比图。具体实施方式下面结合附图和实施例对本专利技术作进一步的说明。如图1所示，一种深埋硐室围岩压力的计算方法，步骤如下：步骤一：构建破坏机制。参考众多关于土压力、地基承载力、挡土墙压力、边坡、浅埋硐室、深埋硐室等上限分析破坏机制和研究成果，依据数值分析方法构建三种破坏机制，分别为单圆弧夹层破坏机制、双圆弧夹层破坏机制、对数螺旋与圆弧夹层破坏机制，然后根据三种破坏机制的优化结果构建最优破坏机制。单圆弧夹层破坏本文档来自技高网...

【技术保护点】
一种深埋硐室围岩压力的计算方法，步骤如下：步骤一：构建破坏机制；步骤二：根据破坏机制建立运动许可的速度场；步骤三：基于极限分析上限定理，由内功率与外功率相等推导并得到破坏载荷的上限解析解；步骤四：在多个破坏荷载中求得最小值，即获得最优上限解。

【技术特征摘要】
1.一种深埋硐室围岩压力的计算方法，步骤如下：步骤一：构建破坏机制；步骤二：根据破坏机制建立运动许可的速度场；步骤三：基于极限分析上限定理，由内功率与外功率相等推导并得到破坏载荷的上限解析解；步骤四：在多个破坏荷载中求得最小值，即获得最优上限解。2.根据权利要求1所述的深埋硐室围岩压力的计算方法，其特征在于：所述步骤一中，首先构建单圆弧夹层破坏机制、双圆弧夹层破坏机制、圆弧与对数螺旋线夹层破坏机制，然后根据三种破坏机制的优化结果构建最优破坏机制。3.根据权利要求2所述的深埋硐室围岩压力的计算方法，其特征在于：所述步骤一中，单圆弧夹层破坏机制由“向下滑落的楔形塌落体+以顶角为圆心的转动圆弧体+三个平动三角形体”构成，双圆弧夹层破坏机制由“向下滑落的楔形塌落体+以顶角为圆心的转动圆弧体+平动三角形体+以顶角为圆心的转动圆弧体+平动三角形体”构成，圆弧与对数螺旋线夹层破坏机制由“向下滑落的楔形塌落体+以顶角为圆心的转动圆弧体+平动三角形体+以顶角为圆心的转动对数螺旋线体+平动三角形体”构成，最优破坏机制由“向下滑落的楔形塌落体+以顶角为圆心的转动圆弧体+n个平动三角形体”构成。4.根据权利要求1所述的深埋硐室围岩压力的计算方法，其特征在于：所述步骤三中，外部功率和内部能量耗散率计算如下：①速度大小及速度间断线长度计算：V1＝V0根据圆弧辐射受剪区速度场的分析得出：根据破坏机制的速度矢量图及几何三角关系得出：各速度间断线的长度递推关系为：GB＝GC1②外部功率计算1)重力做功功率

【专利技术属性】
技术研发人员：张道兵，张佳华，
申请(专利权)人：湖南科技大学，
类型：发明
国别省市：湖南,43