CN113536420A - 一种考虑节理蠕变的大跨度隧道围岩时效安全度分析方法 - Google Patents
一种考虑节理蠕变的大跨度隧道围岩时效安全度分析方法 Download PDFInfo
- Publication number
- CN113536420A CN113536420A CN202110725694.5A CN202110725694A CN113536420A CN 113536420 A CN113536420 A CN 113536420A CN 202110725694 A CN202110725694 A CN 202110725694A CN 113536420 A CN113536420 A CN 113536420A
- Authority
- CN
- China
- Prior art keywords
- rock mass
- creep
- joint
- stage
- rock
- 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
Links
- 239000011435 rock Substances 0.000 title claims abstract description 227
- 230000032683 aging Effects 0.000 title claims abstract description 13
- 238000004458 analytical method Methods 0.000 title claims description 11
- 238000000034 method Methods 0.000 claims abstract description 39
- 230000008569 process Effects 0.000 claims abstract description 18
- 230000001133 acceleration Effects 0.000 claims abstract description 9
- 239000004927 clay Substances 0.000 claims abstract description 7
- 239000000463 material Substances 0.000 claims description 52
- 230000035882 stress Effects 0.000 claims description 42
- 230000006378 damage Effects 0.000 claims description 31
- 230000007774 longterm Effects 0.000 claims description 31
- 238000009795 derivation Methods 0.000 claims description 6
- 230000014509 gene expression Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims description 3
- 238000009412 basement excavation Methods 0.000 abstract description 22
- 238000010276 construction Methods 0.000 abstract description 12
- 230000011218 segmentation Effects 0.000 abstract 1
- 238000006073 displacement reaction Methods 0.000 description 14
- 230000008859 change Effects 0.000 description 6
- 230000008093 supporting effect Effects 0.000 description 6
- 238000012544 monitoring process Methods 0.000 description 4
- 230000001550 time effect Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 239000002689 soil Substances 0.000 description 3
- 230000007423 decrease Effects 0.000 description 2
- 239000011229 interlayer Substances 0.000 description 2
- 230000002787 reinforcement Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 229920000034 Plastomer Polymers 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005489 elastic deformation Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- -1 is gneiss Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Theoretical Computer Science (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Architecture (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种考虑节理蠕变的大跨度隧道围岩时效安全度分析方法,定义了非线性黏壶元件,提出一种将节理与蠕变相结合的非线性西原遍布节理应变软化蠕变模型,描述节理发生屈服后岩质体强度参数软化程度;推导出该模型ZSIsu分段函数。本发明通过引入的非线性黏壶元件建立了非线性西原遍布节理应变软化蠕变模型,弥补了传统模型中无法模拟节理岩体加速蠕变阶段的特点,且将节理与蠕变相结合,充分反应出节理岩体的特性,同时考虑节理、岩质体属性,能更好地反映岩质体节理特征及工程实际。同时结合单元状态指标,推导出该模型的ZSIsu分段函数,在大跨度暗挖隧道施工过程中,实现隧道开挖各个阶段的量化控制,有效指导施工的安全进行。
Description
技术领域
本发明涉及大跨度暗挖隧道施工领域,尤其涉及一种考虑节理蠕变的大跨度隧道围岩时效安全度分析方法。
背景技术
在大跨度隧道施工过程中,隧道结构两侧岩体的变形不仅仅与空间效应相关,往往与时间也有关系,所谓时间效应就是岩体在外力、水力、温度、地质特征及施工扰动等各种环境、地质与工程因素作用下,岩土材料及岩体结构与时间相隔的力学行为、本构关系、失稳与破坏规律。如隧道主体结构开挖初期变形并不明显,但随着时间的增长,隧道的变形越来越严重,甚至出现了失稳破坏。因此,系统地研究岩土体在特定环境和条件下的流变理论的工程应用、更加有效地评价岩体工程的长期稳定性和运营安全,具有重要的理论价值和工程意义。
现阶段,隧道穿越不良岩质体(节理面、断层、夹层)是不可避免的技术难题,但对于不良岩质体(节理面、断层、夹层)隧道尚未形成成熟的工程经验,目前的研究方法中无法模拟岩质体加速蠕变阶段的特点,不能充分反应节理岩体的特性。
发明内容
本发明提供一种考虑节理蠕变的大跨度隧道围岩时效安全度分析方法,以克服目前的研究方法中无法模拟岩质体加速蠕变阶段的特点,不能充分反应节理岩体特性的技术问题。
为了实现上述目的,本发明的技术方案是:
一种考虑节理蠕变的大跨度隧道围岩时效安全度分析方法,包括如下步骤:
S1:定义具有非线性黏滞系数η2(t)的非线性黏壶元件,所述非线性黏壶元件用以表征岩质体的在加速蠕变阶段的特性;
S2:提出一种非线性西原遍布节理应变软化蠕变模型,将围岩节理与蠕变参数相结合,以描述节理发生屈服后,岩质体的黏聚力c和内摩擦角φ所代表的岩质体软化的程度;
S3:推导出所述非线性西原遍布节理应变软化蠕变模型的改进单元安全度分段函数,以定量评价岩质体材料从弹性到屈服最后破坏的安全程度。
进一步的,所述非线性西原遍布节理应变软化蠕变模型包括串联的Hooke弹簧、黏弹性Kelvin体、黏塑性Binham体和应变软化塑性元件;所述黏塑性Binham体包括并联的非线性黏壶元件和塑性元件;所述黏弹性Kelvin体包括并联的第二弹性元件和第二黏壶元件。
进一步的,所述非线性西原遍布节理应变软化蠕变模型的建立过程为:
假设所述非线性西原遍布节理应变软化蠕变模型中的非线性黏滞系数η2(t)的衰减趋势呈指数型关系,因此:
式中:tc为进入加速蠕变阶段的起始时间;α为调节时间量纲的系数;εps是岩质体进入加速蠕变阶段的对应剪切等效塑性应变;εpt是岩质体进入加速蠕变阶段的拉伸等效塑性应变;是岩质体开始进入加速蠕变阶段的等效塑性剪切应变阈值;是岩质体开始进入加速蠕变阶段的等效拉伸塑性应变阈值;
则加速蠕变阶段非线性黏塑性Binham体蠕变速率为:
将公式(2)两端进行积分,得非线性黏塑性体的蠕变方程为:
则所述非线性西原遍布节理应变软化蠕变模型的方程为:
式中:σs是长期屈服应力;E0是Hooke弹簧的弹性模量;E1黏弹性Kelvin体中弹簧的弹性模量;η1是黏弹性Kelvin体中黏壶的黏滞系数;当σ≤σs时,岩质体材料只发生衰减蠕变,即岩质体处于弹性阶段;当σ>σs,t≤tc时,岩质体材料只发生等速蠕变,即岩质体处于屈服阶段;当σ>σs,t>tc时,岩质体材料发生加速蠕变,即岩质体处于破坏阶段;
由于岩质体的长期强度是衡量岩质体工程耐久性和长期稳定性的重要指标,因此在三维空间中,考虑岩质体的长期强度条件下的非线性西原遍布节理应变软化蠕变模型的方程为:
式中:G0是Hooke弹簧的剪切模量;G1是黏弹性Kelvin体剪切模量;Sij是应力偏张量;Ss是岩质体的长期强度;当Sij≤Ss时,所述考虑岩质体的长期强度条件下的非线性西原遍布节理应变软化蠕变模型中的黏塑性Binham体不发挥作用,岩质体产生衰减蠕变,即岩质体处于弹性阶段;当Sij>Ss,t≤tc时,黏塑性Binham体发挥作用,岩质体产生等速蠕变阶段,即岩质体处于屈服阶段;当Sij>Ss,t>tc时,岩质体进入加速阶段,即破坏阶段。
进一步的,所述非线性西原遍布节理应变软化蠕变模型的改进单元安全度的推导过程如下:
S21:在弹性阶段时,岩质体的塑性剪切应变εps=0,εpt=0;
则推导岩质体主应力表示形式下的岩质体单元安全度ZSImc为:
式中:σ1是最大岩质体主应力;σ3是最小岩质体主应力;c是黏聚力;φ是内摩擦角;
对于节理面,岩质体主应力表示形式下的节理面单元安全度ZSIu推导公式为:
式中:Cj为黏聚力;φj为节理内摩擦角;β为节理面倾角;
得出弹性阶段的岩质体单元安全度ZSImc为:
弹性阶段的节理面单元安全度ZSIu为:
在拉伸屈服阶段时,
则可以得出:
在屈服阶段岩质体单元安全度ZSImc为:
在屈服阶段的节理面单元安全度ZSIu为:
S23:在破坏阶段时,ZSI的计算过程如所述S22;
由以上公式可得,基于所述非线性西原遍布节理应变软化蠕变模型的改进单元安全度为:
ZSIsu=min(ZSImc,ZSIu)。
进一步的,所述评价岩质体材料从弹性到屈服最后破坏的安全程度的分析方法为:
S31:对所述非线性西原遍布节理应变软化蠕变模型中原岩和节理岩质体主应力表达式进行变换,以获得在弹性阶段的节理面单元安全度;
S32:将岩质体材料的塑性应变和岩质体材料的极限塑性应变进行对比,以能够量化表达岩质体材料在屈服阶段距离破坏阶段的接近程度,以获得在屈服阶段的节理面单元安全度;
S33:当所述岩质体材料在屈服阶段距离破坏阶段的岩质体材料的塑性应变和岩质体材料的极限塑性应变相等时,岩质体开始进入破坏阶段。
有益效果:本发明的一种考虑节理蠕变的大跨度隧道围岩时效安全度分析方法,引入了非线性黏壶元件,建立了非线性西原遍布节理应变软化蠕变模型,弥补了传统模型中无法模拟节理岩体加速蠕变阶段的特点,并且将节理与蠕变相结合,能充分反应出节理岩体的特性,同时考虑节理、岩质体属性,能更好地反映岩质体节理特征及工程实际。同时结合单元状态指标,推导出节理蠕变方程的ZSI分段函数,在大跨度暗挖隧道施工过程中,实现隧道开挖各个阶段的量化控制,简化施工控制参数。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明分析方法的步骤图;
图2为本发明非线性西原遍布节理应变软化蠕变模型元件组合图;
图3为本发明所建立的非线性西原遍布节理应变软化蠕变模型图;
图4为本发明一个实施例中施工过程中的监测点布置图;
图5a为本发明的无支护隧道开挖后12天时围岩竖向位移等值线图;
图5b为本发明的无支护隧道开挖后48天时围岩竖向位移等值线图;
图5c为本发明的无支护隧道开挖后60天时围岩竖向位移等值线图;
图5d为本发明的无支护隧道开挖后72天时围岩竖向位移等值线图;
图6a为本发明的无支护隧道开挖后12天时围岩竖向位移等值线图;
图6b为本发明的无支护隧道开挖后48天时围岩竖向位移等值线图;
图6c为本发明的无支护隧道开挖后60天时围岩竖向位移等值线图;
图6d为本发明的无支护隧道开挖后72天时围岩竖向位移等值线图;
图7为本发明的无支护隧道围岩监测点的位移变化曲线图;
图8a为本发明的无支护隧道开挖后12天的ZSI等值线图;
图8b为本发明的无支护隧道开挖后48天的ZSI等值线图;
图8c为本发明的无支护隧道开挖后60天的ZSI等值线图;
图8d为本发明的无支护隧道开挖后72天的ZSI等值线图。
其中:E0是Hooke弹簧的弹性模量;E1是黏弹性Kelvin体中弹簧的弹性模量,η1是黏弹性Kelvin体中黏壶的黏滞系数;σs是基于Mohr-Coulomb准则的长期屈服应力;σS-J是遍布节理应变软化模型短期加载屈服应力;εH是Hooke体弹性应变,εk是黏弹性Kelvin体黏弹性应变,εM是黏塑性Binham体中的黏塑性应变,εp是应变软化塑性元件的塑性应变。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本实施例提供了一种考虑节理蠕变的大跨度隧道围岩时效安全度分析方法,如图1所示,包括如下步骤:
S1:定义具有非线性黏滞系数η2(t)的非线性黏壶元件,所述非线性黏壶元件用以表征岩质体的在加速蠕变阶段的特性;
S2:提出一种非线性西原遍布节理应变软化蠕变模型,将围岩节理与蠕变参数相结合,以描述节理发生屈服后,岩质体的强度参数黏聚力c和内摩擦角φ所代表的岩质体软化的程度;;所述非线性西原遍布节理应变软化蠕变模型将传统的西原模型中的线性黏壶元件替换为非线性黏壶元件。由于传统的西原模型基本元件均为线性元件,而岩质体在蠕变的过程中伴随着应变的增加出现内部损伤,一般认为损伤是岩质体材料在加载条件下其黏聚力呈渐进性减弱,进而导致其体积元裂化和破坏的现象,影响岩质体的有效抗剪强度参数,这种岩质体软化现象,通过非线性黏壶元件的引用便可以获得较好的表达。所述非线性西原遍布节理应变软化蠕变模型的元件组合图如附图1所示。
所述非线性西原遍布节理应变软化蠕变模型包括串联的Hooke弹簧、黏弹性Kelvin体、黏塑性Binham体和应变软化塑性元件;所述黏塑性Binham体包括并联的非线性黏壶元件和塑性元件;所述黏弹性Kelvin体包括并联的第二弹性元件和第二黏壶元件。具体的,在本实施例中,所述第二弹性元件是弹簧结构。其中,模型中的η2(t)是与时间有关的非线性黏性系数;如附图2所示。
当岩质体开始进入加速蠕变时,对应剪切等效塑性应变阈值和拉伸等效塑性应变阈值当模型的剪切等效塑性应变阈值εps或者拉伸等效塑性应变阈值εpt小于相应的阈值时,所述非线性黏壶元件转化为线性黏壶元件,此时η2(t)的值为初始值η2并保持不变;当模型的剪切等效塑性应变阈值εps或者拉伸等效塑性应变阈值εpt大于相应阈值时,所述非线性黏壶元件的非线性特征被触发,η2(t)变成随时间增加而衰减的非线性参数;非线性黏滞系数η2(t)随时间不断减小,初始时,减小速度较慢,但当随着时间的增加,非线性黏滞系数η2(t)减小的速度逐渐加快,最终减小为0,从而使得岩质体由于应变过大而被破坏。因岩质体的损伤变量与时间呈指数关系,因此所述非线性西原遍布节理应变软化蠕变模型中的非线性黏滞系数η2(t)的衰减趋势呈指数型关系,因此:
式中:tc为进入加速蠕变阶段的起始时间;α为调节时间量纲的系数;εps是岩质体进入加速蠕变阶段的对应剪切等效塑性应变;εpt是岩质体进入加速蠕变阶段的拉伸等效塑性应变;是岩质体开始进入加速蠕变阶段的等效塑性剪切应变阈值;是岩质体开始进入加速蠕变阶段的等效拉伸塑性应变阈值;
则加速蠕变阶段非线性黏塑性Binham体体蠕变速率为:
将公式(2)两端进行积分,得非线性黏塑性Binham体的蠕变方程为:
则所述非线性西原遍布节理应变软化蠕变模型的方程为:
式中:σs是长期屈服应力;E0是Hooke弹簧的弹性模量;E1是Kelvin体中弹簧的弹性模量;η1是Kelvin体中黏壶的黏滞系数;当σ≤σs时,岩质体材料只发生衰减蠕变,即岩质体处于弹性阶段;当σ>σs,t≤tc时,岩质体材料只发生等速蠕变,即岩质体处于屈服阶段;当σ>σs,t>tc时,岩质体材料发生加速蠕变,即岩质体处于破坏阶段;
此时,材料的长期强度用σs来代替,以σ≤σs和σ>σs,来表示岩质体发生蠕变的不同阶段;
长期强度是衡量岩质体工程耐久性和长期稳定性的重要指标,确定岩质体的长期强度可根据岩质体的蠕变数据作岩质体破坏应力和破坏前历时时间的长期强度曲线,通过长期强度曲线的渐进线确定长期强度,或者根据岩质体等时应力-应变曲线,做其时间趋近于无穷大时的水平渐近线,该水平渐近线与应力坐标轴交点的应力即为岩质体长期强度。
在三维空间中,三维应力状态下应力张量σij可以分解为应力球张量σm和应力偏张量Sij,应变张量球张量εm和应变偏张量εij,σm应力球张量只改变物体的体积,而不能改变其形状,而应力偏张量Sij只影响形状的变化而不引起体积改变。
因而在考虑岩质体的长期强度条件下的非线性西原遍布节理应变软化蠕变模型的方程建立如下:
式中:G0是Hooke弹簧的剪切模量;G1是黏弹性Kelvin体剪切模量;Sij是应力偏张量;Ss是岩质体的长期强度;
由于三维模型的蠕变方程中存在长期强度Ss,因此,三维情况下的蠕变方程由公式(5)所示分段函数表示,蠕变性质与岩质体所处的应力状态有关。
(1)当Sij≤Ss时,所述考虑岩质体的长期强度条件下的非线性西原遍布节理应变软化蠕变模型中的黏塑性Binham体不发挥作用,产生蠕变第1阶段,即衰减蠕变的弹性阶段。
(2)当Sij>Ss,t≤tc时,应力水平超过长期强度,所述考虑岩质体的长期强度条件下的非线性西原遍布节理应变软化蠕变模型中的黏塑性Binham体发挥作用,但应变未超过加速段应变阈值,变形进入蠕变第2阶段即等速蠕变阶段的屈服阶段。
(3)当Sij>Ss,t>tc时,应力水平超过长期强度,等效塑性应变超过加速段的临界应变,变形进入蠕变第3阶段即加速蠕变阶段,变形迅速增大,材料发生破坏,即为破坏阶段;
此时,岩质体发生蠕变的各个阶段用Sij≤Ss和Sij>Ss来判别,这里的判别条件是6个独立的应力分量,岩质体处于哪一个阶段与其应力状态在三维空间中到屈服面的距离有关。
S3:推导出所述非线性西原遍布节理应变软化蠕变模型的改进单元安全度分段函数,以定量评价岩土材料从弹性到屈服到最后破坏的安全程度;推导过程如下:
岩质体强度破坏理论是判断岩体失稳的重要判据,而Mohr-Coulomb准则是岩土界应用广泛的屈服准则。
S21:在弹性阶段时,此阶段的岩质体的塑性剪切应变εps=0,εpt=0;由Mohr-Coulomb准则岩质体主应力表达式推导岩质体单元安全度ZSImc为:
式中:σ1是最大岩质体主应力;σ3是最小岩质体主应力;c是黏聚力;φ是内摩擦角;
对于节理面,岩质体主应力表示形式下的节理面单元安全度ZSIu推导公式为:
式中:Cj为节理黏聚力;φj为节理内摩擦角;β为节理面倾角;
即得出弹性阶段的岩质体单元安全度ZSImc为:
弹性阶段的节理面单元安全度ZSIu为:
当岩质体材料在拉伸屈服阶段时,
则可以得出:
在屈服阶段岩质体单元安全度ZSImc为:
在屈服阶段的节理面单元安全度ZSIu为:
在屈服阶段的ZSI∈(0,1),其值越小,表示岩质体材料被破坏程度越高。
S23:在破坏阶段时,ZSI的计算过程如所述S22;
综合上公式可得,基于所述非线性西原遍布节理应变软化蠕变模型的改进单元安全度为:
ZSIsu=min(ZSImc,ZSIu)。
所述评价岩质体材料从弹性到屈服最后破坏的安全程度的分析方法为:
S31:对所述非线性西原遍布节理应变软化蠕变模型中原岩和节理岩质体主应力表达式进行变换,以获得在弹性阶段的节理面单元安全度;
S32:将岩质体材料的塑性应变和岩质体材料的极限塑性应变进行对比,以能够量化表达岩质体材料在屈服阶段距离破坏阶段的接近程度,以获得在屈服阶段的节理面单元安全度;
S33:当所述岩质体材料在屈服阶段距离破坏阶段的岩质体材料的塑性应变和岩质体材料的极限塑性应变相等时,岩质体开始进入破坏阶段。
具体的,本实施例中选取某地区隧道施工点为分析对象,根据现场勘察发现,隧道围岩主要为片麻岩,岩体节理裂隙发育,根据现场采样得到的试验数据,根据本发明所公开的内容,建立隧道围岩(片麻岩)的非线性西原遍布节理应变软化蠕变模型。根据现场地质状况,建立三维数值模型,对施工过程进行分析。如附图3和4所示为所建立的三维数值模型及施工过程中的监测点布置图。
为了更好地分析隧道开挖后围岩的变形、应力状态等随时间效应的变化规律,首先对无支护结构隧道进行模拟分析,以充分了解隧洞开挖后时间效应在不受其他支护力的干扰下对隧道结构的影响,同时也可为支护效果的评价提供对比的依据。
由于围岩体是一种流变介质,隧道开挖后初期,新的临空面的形成,使得围岩的应力、应变状态变化速度较快,随后围岩的应力、应变将缓慢变化。对此,根据数值模型的模拟结果,本实施例对隧道开挖完成后不同时间的模拟结果进行分析和研究。如附图5和图6分别为无支护结构隧道开挖后不同时间点围岩竖向位移等值线图及水平位移等值线图,由图可见,随着时间的增加,隧道周边位移逐渐增大。
从图7所示的隧道围岩监测点的位移变化曲线可以看出,在隧道形成后的6天内围岩变形速率较快,6天后围岩变形速率趋于稳定,说明围岩的应力状态已经超过了长期强度,黏塑性体中的非线性黏糊元件发生了作用。随着时间的增加,变形量不断增大,在63天左右围岩的变形量超过了加速段等效塑性应变阈值,拱顶,拱肩和边墙的位移速率急剧增加,从开始进入加速变形阶段后,在11天的时间内,隧道便失稳破坏。从量值来看,隧道围岩的初始的弹性变形量较小,拱顶沉降量在瞬时平衡后仅为6.8mm左右,而时间效应所产生的变形量已远远超过该量值。
围岩的极限等效塑性切应变如图8为蠕变过程中,无支护结构下不同时间点围岩的ZSI等值线变化云图,从图中可以看出,隧道结构的在开挖后最危险部位主要分布在两侧边墙处,最小值为-1。随着时间的增长,围岩的破坏区域区主要分布隧道拱顶和底板两个部位并沿着铅垂方向扩展,最后在加速段出现大面积的破坏区域。因此陈家店隧道开挖后如不及时支护,即使在短期内隧道结构是稳定的,但一旦时间较长,由于软岩的流变性,隧道的结构仍然可能被破坏。
工程实际中,随着时间的发展,围岩的强度和稳定性将呈现出持续降低的趋势。因此,为了保证隧道的正常使用,需要根据围岩的变形和破坏规律制定合理的支护方案,控制围岩的变形量及减小屈服区范围。
根据本实施例的分析结果来看,为了降低围岩流变性及节理存在对隧道结构产生的破坏,在隧道施工时采用超前小管棚,超前帷幕预注浆等加固措施,在隧道周围一定范围内形成加固保护区,能够有效的增强隧道结构的强度。
最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (5)
1.一种考虑节理蠕变的大跨度隧道围岩时效安全度分析方法,其特征在于,包括如下步骤:
S1:定义具有非线性黏滞系数η2(t)的非线性黏壶元件,所述非线性黏壶元件用以表征岩质体的在加速蠕变阶段的特性;
S2:提出一种非线性西原遍布节理应变软化蠕变模型,将围岩节理与蠕变参数相结合,以描述节理发生屈服后,岩质体的黏聚力c和内摩擦角φ所代表的岩质体软化的程度;
S3:推导出所述非线性西原遍布节理应变软化蠕变模型的改进单元安全度分段函数,以定量评价岩质体材料从弹性到屈服最后破坏的安全程度。
2.根据权利要求1所述的一种考虑节理蠕变的大跨度隧道围岩时效安全度分析方法,其特征在于,所述非线性西原遍布节理应变软化蠕变模型包括串联的Hooke弹簧、黏弹性Kelvin体、黏塑性Binham体和应变软化塑性元件;所述黏塑性Binham体包括并联的非线性黏壶元件和塑性元件;所述黏弹性Kelvin体包括并联的第二弹性元件和第二黏壶元件。
3.根据权利要求1所述的一种考虑节理蠕变的大跨度隧道围岩时效安全度分析方法,其特征在于,所述非线性西原遍布节理应变软化蠕变模型的建立过程为:
所述非线性西原遍布节理应变软化蠕变模型中的非线性黏滞系数η2(t)的衰减趋势呈指数型关系,因此:
式中:tc为进入加速蠕变阶段的起始时间;α为调节时间量纲的系数;εps是岩质体进入加速蠕变阶段的对应剪切等效塑性应变;εpt是岩质体进入加速蠕变阶段的拉伸等效塑性应变;是岩质体开始进入加速蠕变阶段的等效塑性剪切应变阈值;是岩质体开始进入加速蠕变阶段的等效拉伸塑性应变阈值;
则加速蠕变阶段非线性黏塑性Binham体蠕变速率为:
将公式(2)两端进行积分,得非线性黏塑性体的蠕变方程为:
则所述非线性西原遍布节理应变软化蠕变模型的方程为:
式中:σs是长期屈服应力;E0是Hooke弹簧的弹性模量;E1黏弹性Kelvin体中弹簧的弹性模量;η1是黏弹性Kelvin体中黏壶的黏滞系数;当σ≤σs时,岩质体材料只发生衰减蠕变,即岩质体处于弹性阶段;当σ>σs,t≤tc时,岩质体材料只发生等速蠕变,即岩质体处于屈服阶段;当σ>σs,t>tc时,岩质体材料发生加速蠕变,即岩质体处于破坏阶段;
由于岩质体的长期强度是衡量岩质体工程耐久性和长期稳定性的重要指标,因此在三维空间中,考虑岩质体的长期强度条件下的非线性西原遍布节理应变软化蠕变模型的方程为:
式中:G0是Hooke弹簧的剪切模量;G1是黏弹性Kelvin体剪切模量;Sij是应力偏张量;Ss是岩质体的长期强度;当Sij≤Ss时,所述考虑岩质体的长期强度条件下的非线性西原遍布节理应变软化蠕变模型中的黏塑性Binham体不发挥作用,岩质体产生衰减蠕变,即岩质体处于弹性阶段;当Sij>Ss,t≤tc时,黏塑性Binham体发挥作用,岩质体产生等速蠕变阶段,即岩质体处于屈服阶段;当Sij>Ss,t>tc时,岩质体进入加速阶段,即破坏阶段。
4.根据权利要求3所述的一种考虑节理蠕变的大跨度隧道围岩时效安全度分析方法,其特征在于,所述非线性西原遍布节理应变软化蠕变模型的改进单元安全度的推导过程如下:
S21:在弹性阶段时,岩质体的塑性剪切应变εps=0,εpt=0;
则推导岩质体主应力表示形式下的岩质体单元安全度ZSImc为:
式中:σ1是最大岩质体主应力;σ3是最小岩质体主应力;c是聚力;φ是内摩擦角;
对于节理面,岩质体主应力表示形式下的节理面单元安全度ZSIu推导公式为:
式中:Cj为节理黏聚力;φj为节理内摩擦角;β为节理面倾角;
得出弹性阶段的岩质体单元安全度ZSImc为:
弹性阶段的节理面单元安全度ZSIu为:
在拉伸屈服阶段时,
则可以得出:
在屈服阶段岩质体单元安全度ZSImc为:
在屈服阶段的节理面单元安全度ZSIu为:
S23:在破坏阶段时,ZSI的计算过程如所述S22;
由以上公式可得,基于所述非线性西原遍布节理应变软化蠕变模型的改进单元安全度为:
ZSIsu=min(ZSImc,ZSIu)。
5.根据权利要求4所述的一种考虑节理蠕变的大跨度隧道围岩时效安全度分析方法,其特征在于,所述评价岩质体材料从弹性到屈服最后破坏的安全程度的分析方法为:
S31:对所述非线性西原遍布节理应变软化蠕变模型中原岩和节理的岩质体主应力表达式进行变换,以获得在弹性阶段的节理面单元安全度;
S32:将岩质体材料的塑性应变和岩质体材料的极限塑性应变进行对比,以能够量化表达岩质体材料在屈服阶段距离破坏阶段的接近程度,以获得在屈服阶段的节理面单元安全度;
S33:当所述岩质体材料在屈服阶段距离破坏阶段的岩质体材料的塑性应变和岩质体材料的极限塑性应变相等时,岩质体开始进入破坏阶段。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110725694.5A CN113536420A (zh) | 2021-06-29 | 2021-06-29 | 一种考虑节理蠕变的大跨度隧道围岩时效安全度分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110725694.5A CN113536420A (zh) | 2021-06-29 | 2021-06-29 | 一种考虑节理蠕变的大跨度隧道围岩时效安全度分析方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113536420A true CN113536420A (zh) | 2021-10-22 |
Family
ID=78126163
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110725694.5A Pending CN113536420A (zh) | 2021-06-29 | 2021-06-29 | 一种考虑节理蠕变的大跨度隧道围岩时效安全度分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113536420A (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012112111A (ja) * | 2010-11-22 | 2012-06-14 | Shimizu Corp | 安定性評価方法および安定性評価装置 |
CN108959758A (zh) * | 2018-06-27 | 2018-12-07 | 大连海大恒信土木工程技术有限公司 | 一种考虑渗流蠕变的隧道局部安全性分析方法 |
CN110274835A (zh) * | 2019-07-12 | 2019-09-24 | 中南大学 | 一种改进Burgers岩石剪切蠕变模型的方法 |
-
2021
- 2021-06-29 CN CN202110725694.5A patent/CN113536420A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012112111A (ja) * | 2010-11-22 | 2012-06-14 | Shimizu Corp | 安定性評価方法および安定性評価装置 |
CN108959758A (zh) * | 2018-06-27 | 2018-12-07 | 大连海大恒信土木工程技术有限公司 | 一种考虑渗流蠕变的隧道局部安全性分析方法 |
CN110274835A (zh) * | 2019-07-12 | 2019-09-24 | 中南大学 | 一种改进Burgers岩石剪切蠕变模型的方法 |
Non-Patent Citations (4)
Title |
---|
卢小雨;王游;: "非线性西原体模型及其应用", 河南城建学院学报, no. 04 * |
张启航: "《沿海山岭大跨度小净距隧道施工稳定性分析》", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》, no. 6, pages 1 - 81 * |
张旭东;: "基于修正西原模型与Lade-Duncan准则的隧道围岩蠕变效应研究", 科学技术与工程, no. 18 * |
张权: "《基于改进单元安全度ZSI的富水区隧道稳定性分析》", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》, no. 1, pages 1 - 71 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Idriss et al. | Nonlinear behavior of soft clays during cyclic loading | |
CN108959758B (zh) | 一种考虑渗流蠕变的隧道局部安全性分析方法 | |
Xiao et al. | A unified strain-hardening and strain-softening elastoplastic constitutive model for intact rocks | |
Chandler | Quantifying long-term strength and rock damage properties from plots of shear strain versus volume strain | |
CN110457853B (zh) | 一种含结构面的拱坝坝肩岩体失稳判断方法 | |
CN113376005B (zh) | 基于颗粒长期强度水平的堆石料最终流变应变预测方法 | |
Gooch et al. | Application of stochastic numerical analyses in the assessment of spatially variable unreinforced masonry walls subjected to in-plane shear loading | |
Zhang et al. | Damage constitutive model of uniaxially compressed coal material considering energy dissipation | |
Baktheer et al. | Modeling of bond fatigue in reinforced concrete based on cumulative measure of slip | |
LI et al. | Effect of loading rate change on the mechanical properties of mudstone under uniaxial compression | |
D'Altri et al. | Environmentally-induced loss of performance in FRP strengthening systems bonded to full-scale masonry structures | |
Kamdem et al. | New description of the mechanical creep response of rocks by fractional derivative theory | |
Mašın et al. | Numerical analyses of a tunnel in London clay using different constitutive models | |
Cai et al. | An improved hardening-damage creep model of lean coal: a theoretical and experimental study | |
Haghgou et al. | Determination of optimal intensity measure for probabilistic seismic demand analysis of intake towers | |
Li et al. | A modified cutting‐plane time integration scheme with adaptive substepping for elasto‐viscoplastic models | |
El Yassari et al. | Numerical simulation of fiber-reinforced concrete under cyclic loading using extended finite element method and concrete damaged plasticity | |
CN113536420A (zh) | 一种考虑节理蠕变的大跨度隧道围岩时效安全度分析方法 | |
Yao et al. | SSUH model: A small-strain extension of the unified hardening model | |
Yuan et al. | Model for predicting and controlling creep settlements with surcharge loading | |
CN104196024B (zh) | 基于锚索预应力损失与岩土体蠕变耦合的预应力损失计算方法 | |
PENG et al. | Research on elasto-plastic viscous fatigue constitutive model of jointed rock | |
He et al. | Analytical and numerical solutions for shear mechanical behaviors of structural plane | |
CN113806966A (zh) | 基于流变模型应用的非线性岩石疲劳本构模型的构建方法 | |
CN112329224B (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 |