考虑围岩峰后行为的隧道初期支护预留变形量确定方法

技术领域

本发明属于隧道支护结构设计技术领域，涉及一种考虑围岩峰后行为的隧道初期支护预留变形量确定方法。

背景技术

人们早在19世纪即开始探索隧道的支护设计理论，早期工程界将衬砌作为主要受力构件，承受上部围岩荷载，将衬砌厚度设计较大。20世纪60年代后，新奥法作为新兴的隧道设计方法得到关注，新奥法的基本指导思想是充分发挥围岩自身承载能力，采用喷混凝土、锚杆等柔性结构作为主要初期支护形式，通过适时加固围岩、控制围岩变形来充分发挥围岩自承能力。过早支护时，支护结构将承受较大变形压力，并不经济；而过迟支护时，围岩会由于变形较大而发生松弛失稳。因此，围岩开挖过程中初次支护存在合理的支护时机，在该条件下，围岩变形与支护受力都不至于过大。合理支护时机的确定，可确保围岩稳定并节约成本造价，具有重要的现实意义。

隧道约束收敛理论是新奥法的基础，该理论认为，可通过初期支护的预留变形量直观反映合理支护时机。洞室开挖所预留的变形量小时，初期支护承受较大荷载，支护可能发生破坏，而变形量较大，则围岩发生失稳坍塌的可能性大。尽管已有众多研究对合理预留变形量的确定展开了讨论，但所提判据众多，很少从定量角度给出某一地质条件下支护时机，难以被工程人员使用。另外，隧道开挖之后，由于应力释放，围岩初始应力场将重新分布，当围岩强度较小或初始应力较大时，围岩发生塑性屈服破坏，破坏区域范围与塑性变形对预留变形量的预测产生直接影响，故准确运用岩体力学模型与塑性流动法尤为重要。应变软化效应与剪胀效应是岩石峰后的典型力学特征，然而，现有预留变形量估算方法对此涉及较少。

发明内容

为了解决背景技术中存在的上述技术问题，本发明提供了一种降低在支护时机设计方面的主观性、可快速计算预留变形量具体范围、具有直观及简捷的考虑围岩峰后行为的隧道初期支护预留变形量确定方法。

为了实现上述目的，本发明采用如下技术方案：

一种考虑围岩峰后行为的隧道初期支护预留变形量确定方法，其特征在于：所述考虑围岩峰后行为的隧道初期支护预留变形量确定方法包括以下步骤：

1)通过施工现场勘探、测试与分析，获得围岩基本力学参数与支护力学参数；所述围岩基本力学参数包括初始地应力、隧道半径、应变软化系数、剪胀系数、强度参数与变形参数；所述支护力学参数包括支护的厚度、弹性模量、泊松比与单轴抗压强度；

2)建立隧道开挖理论模型并基于该理论模型对最大围岩变形、最大围岩塑性范围半径、模型每层环向应变与径向应力一一对应关系、掌子面开挖时围岩初始变形进行求解；

3)根据步骤2)的求解结果计算预留变形量的上限阈值；

4)根据步骤2)的求解结果计算预留变形量的下限阈值；

5)根据预留变形量的上限阈值与预留变形量的下限阈值来明确预留变形量的具体范围，施工时从该范围中选取合适的预留变形量取值，直接为工程施工现场使用。

作为优选，本发明所采用的步骤2)的具体实现方式是：

2.1)基于隧道开挖理论模型，求解弹性与塑性区域边界处径向应力，作为围岩塑性区域分层时径向应力初始值：

式中：

σ

σ

c

σ

2.2)基于隧道开挖理论模型，根据步骤2.1)的结果求解塑性区域每层径向应力增量：

Δσ

式中：

n是塑性区域内部总分层数；

p

Δσ

2.3)基于隧道开挖理论模型，求解应变软化系数初始值、切向应变初始值、径向应变初始值：

当i＝0时，应变软化系数初始值为0：

η

式中：

η

当i＝0时，切向应变与径向应变初始值为：

ε

ε

式中：

ε

ε

2.4)基于隧道开挖理论模型，经过推导，可给出塑性区域每层中切向应变表达式、径向应力表达式、围岩变形表达式：

σ

u

式中：

σ

i是塑性区域内部某一层数，始于塑性与弹性区域交界处，止于隧道边界；

ε

ε

r

r

K

式(5a)中的计算参量r

r

式中：

ε

ε

Δσ

μ是围岩泊松比，单位是无量纲；

K

σ

σ

E为围岩变形模量，单位是MPa；

σ

式(7)中，

式(6)中，H(σ

式(8)、(9)中：

c

c

式(10)、(11)中：

c

c

η

η

2.5)基于隧道开挖理论模型，经过推导，可给出围岩塑性范围半径表达式：

式中：

R

R

Δσ

σ

η

H(σ

2.6)基于隧道开挖理论模型，通过逐步求解算得支护压力为0时最大围岩变形、最大围岩塑性范围半径、每层环向应变与径向应力的一一对应关系：

以上式(5)到式(13)左侧变量均指第i层变量，右侧所含计算变量均指第i-1层变量与已知量，故基于初始变量，根据式(5)到(13)，对式(5c)进行i＝0到i＝n逐步求解，得到当支护压力p

ε

2.7)根据最大围岩变形、最大塑性区域半径获得掌子面开挖时围岩初始变形：

式中：

u

u

R

作为优选，本发明所采用的步骤3)的具体实现方式是：

3.1)根据围岩变形减小比率与围岩最大变形算得最终围岩变形：

u

式中：

f

u

u

3.2)根据切向应变与径向应力的一一对应关系来确定围岩与支护接触后隧道边界的最终支护压力；

假设在第i环处，切向应变根据下式计算得到：

ε

式中：

ε

u

R

将ε

p

故根据式(14)、(17)、(18)获取围岩与支护接触后隧道边界的最终支护压力p

3.3)通过支护弹性模量、支护厚度、支护泊松比算得支护刚度：

式中：

R

t

E

K

μ

3.4)根据步骤3.2)、步骤3.3)获得的最终支护压力与支护刚度算得支护变形：

Δt

式中：

Δt

K

p

3.5)根据步骤3.4)、步骤3.1)以及步骤2.7)获得的支护变形、最终围岩变形与掌子面开挖时围岩初始变形算得预留变形量的上限阈值：

ΔR

式中：

ΔR

u

u

Δt

作为优选，本发明所采用的步骤4)的具体实现方式是：

4.1)根据支护单轴抗压强度、支护厚度、隧道半径求解支护所能承受最大支护压力：

式中：

R

t

σ

p

4.2)根据步骤4.1)与步骤3.3)获得的支护所能承受最大压力与初期支护刚度求解支护变形：

Δt

式中：

K

p

Δt

4.3)对步骤2.1)中围岩变形表达式展开逐步求解，假定支护压力p

4.4)根据步骤4.2)、步骤4.3)、步骤2.7)获得的支护变形、最终围岩变形和掌子面开挖时围岩初始变形得到预留变形量的下限阈值：

ΔR

式中：

ΔR

u

u

Δt

本发明具有以下优点：

本发明提供了一种考虑围岩峰后行为的隧道初期支护预留变形量确定方法，包括：1)通过施工现场勘探、测试与分析，获得围岩基本力学参数与支护力学参数；所述围岩基本力学参数包括初始地应力、隧道半径、应变软化系数、剪胀系数、强度参数与变形参数；所述支护力学参数包括支护的厚度、弹性模量、泊松比与单轴抗压强度；2)建立隧道开挖理论模型并基于该理论模型对最大围岩变形、最大围岩塑性范围半径、模型每层环向应变与径向应力一一对应关系、掌子面开挖时围岩初始变形进行求解；3)根据步骤2)的求解结果计算预留变形量的上限阈值；4)根据步骤2)的求解结果计算预留变形量的下限阈值；5)根据预留变形量的上限阈值与预留变形量的下限阈值来明确预留变形量的具体范围，施工时从该范围中选取合适的预留变形量取值，直接为工程施工现场使用。本发明采用理论方法推导出了预留变形量上限与下限阈值的计算公式，具有严密的逻辑性，故算得的预留变形量范围具有合理性；预留变形量计算公式考虑了隧道开挖扰动后围岩峰后的应变软化与剪胀效应，使得预留变形量具体范围与之前方法相比更符合实际。通过代入相关参数，可快速便捷获取预留变形量上限与下限阈值，方法简洁、结果具体明确，可直接为工程施工现场使用。本发明可充分结合围岩峰后力学行为包括应变软化、剪胀特性对围岩预留变形量范围进行计算，降低在支护时机设计方面的主观性，提高其科学性；同时，可快速计算预留变形量具体范围，具有直观、简捷的优点。

附图说明

图1为预留变形量计算流程图；

图2为隧道开挖理论模型示意图。

具体实施方式

下面结合附图，对发明的技术方案作进一步具体说明。

本发明为一种考虑围岩峰后行为的隧道初期支护预留变形量确定方法，实施流程图如图1所示，具体步骤如下：

1)通过施工现场勘探、测试与分析，获得围岩基本力学参数如初始地应力、隧道半径、应变软化系数、剪胀系数、强度参数、变形参数等，获取初期支护的厚度、弹性模量、泊松比、单轴抗压强度等(复合式初期支护相关参数可进行换算)。围岩基本力学参数与支护力学参数分别列于表1与表2中。

表1围岩基本力学参数

表2支护基本力学参数

2)建立隧道开挖理论模型(如图2所示)，基于该理论模型对最大围岩变形、最大围岩塑性范围半径、模型每层环向应变与径向应力一一对应关系、掌子面开挖时围岩初始变形进行求解，具体包括以下子步骤：

2.1)基于隧道开挖理论模型，求解弹性与塑性区域边界处径向应力，作为围岩塑性区域分层时径向应力初始值：

式中：

σ

σ

c

σ

将c

2.2)基于隧道开挖理论模型，求解塑性区域每层径向应力增量：

Δσ

式中：

n是塑性区域内部总分层数；

p

Δσ

此算例中，将n设置为1000，由于需要计算的围岩变形与塑性区域半径均对应最大值状态，故将支护压力p

2.3)基于隧道开挖理论模型，求解应变软化系数初始值、切向应变初始值、径向应变初始值：

当i＝0时，应变软化系数初始值为0：

η

式中：

η

当i＝0时，切向应变与径向应变初始值为：

ε

ε

式中：

ε

ε

将σ

2.4)基于隧道开挖理论模型，经过推导，可给出塑性区域每层中切向应变表达式、径向应力表达式、围岩变形表达式：

ε

σ

u

式中：

σ

i是塑性区域内部某一层数，始于塑性与弹性区域交界处，止于隧道边界；

ε

ε

r

r

K

式(5a)中的计算参量r

r

式中：

ε

ε

Δσ

μ是围岩泊松比，单位是无量纲；

K

σ

σ

E为围岩变形模量，单位是MPa；

σ

式(7)中，

式(6)中，H(σ

式(8)、(9)中：

c

c

式(10)、(11)中：

c

c

η

η

η

2.5)基于隧道开挖理论模型，经过推导，可给出围岩塑性范围半径表达式：

式中：

R

R

Δσ

σ

η

H(σ

2.6)基于隧道开挖理论模型，通过逐步求解算得支护压力为0时最大围岩变形、最大围岩塑性范围半径、每层环向应变与径向应力的一一对应关系：

以上式(5)到式(13)左侧变量均指第i层变量，右侧所含计算变量均指第i-1层变量与已知量，故基于初始变量，根据式(5)到(13)，对式(5c)进行i＝0到i＝n逐步求解，可得到当支护压力p

ε

2.7)根据最大围岩变形、最大塑性区域半径获得掌子面开挖时围岩初始变形：

式中：

u

u

R

将u

3)根据步骤2)的求解结果计算预留变形量的上限阈值，具体包括以下子步骤：

3.1)根据围岩变形减小比率与围岩最大变形算得最终围岩变形：

u

式中：

f

u

u

假定所需f

3.2)根据切向应变与径向应力的一一对应关系来确定围岩与支护接触后隧道边界的最终支护压力。

假设在第i环处，切向应变可根据下式计算得到：

ε

式中：

ε

u

R

将u

将ε

围岩与支护接触后隧道边界的最终支护压力p

p

因此可得到：p

3.3)通过支护弹性模量、支护厚度、支护泊松比算得支护刚度：

式中：

R

t

E

K

μ

将E

3.4)根据步骤3.2)、步骤3.3)获得的最终支护压力与支护刚度算得支护变形：

Δt

式中：

Δt

K

p

将p

3.5)根据步骤3.4)、步骤3.1)、步骤2.4)获得的支护变形、最终围岩变形与掌子面开挖时围岩初始变形算得预留变形量的上限阈值：

ΔR

式中：

ΔR

u

u

Δt

将u

4)根据步骤2)的求解结果计算预留变形量的下限阈值，具体包括以下子步骤：

4.1)根据支护单轴抗压强度、支护厚度、隧道半径求解支护所能承受最大支护压力：

式中：

R

t

σ

p

将R

4.2)根据步骤4.1)与步骤3.3)获得的支护所能承受最大压力、初期支护刚度求解支护变形：

Δt

式中：

K

p

Δt

K

4.3)对步骤2.1)中围岩变形表达式展开逐步求解，假定支护压力p

假定p

4.4)根据步骤4.2)、步骤4.3)、步骤2.7)获得的支护变形、最终围岩变形和掌子面开挖时围岩初始变形得到预留变形量的下限阈值：

ΔR

式中：

ΔR

u

u

Δt

可根据步骤2.7)(见式(15))算得，将u

5)根据预留变形量的上限阈值与预留变形量的下限阈值来明确预留变形量的具体范围，施工时从该范围中选取合适的预留变形量取值，直接为工程施工现场使用。

ΔR

将ΔR