隧道工程—第5章2.ppt
- 【下载声明】
1. 本站全部试题类文档,若标题没写含答案,则无答案;标题注明含答案的文档,主观题也可能无答案。请谨慎下单,一旦售出,不予退换。
2. 本站全部PPT文档均不含视频和音频,PPT中出现的音频或视频标识(或文字)仅表示流程,实际无音频或视频文件。请谨慎下单,一旦售出,不予退换。
3. 本页资料《隧道工程—第5章2.ppt》由用户(罗嗣辉)主动上传,其收益全归该用户。163文库仅提供信息存储空间,仅对该用户上传内容的表现方式做保护处理,对上传内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(点击联系客服),我们立即给予删除!
4. 请根据预览情况,自愿下载本文。本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
5. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007及以上版本和PDF阅读器,压缩文件请下载最新的WinRAR软件解压。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 隧道 工程
- 资源描述:
-
1、中南大学隧道与地下工程系中南大学隧道与地下工程系隧 道 工 程12概述解析法数值分析方法剪切滑移法特征曲线法3 岩体力学方法岩体力学方法的出发点是支护结构与围岩相互的出发点是支护结构与围岩相互作用,组成一个共同承载体系,其中围岩为主要的承作用,组成一个共同承载体系,其中围岩为主要的承载结构。载结构。 计算模式为计算模式为地层地层结构模式结构模式,即处于无限或半,即处于无限或半无限介质中的结构和镶嵌在围岩孔洞上的支护结构所无限介质中的结构和镶嵌在围岩孔洞上的支护结构所组成的组成的复合模式。复合模式。 特点特点:能反映出隧道开挖后的围岩应力状态能反映出隧道开挖后的围岩应力状态。概述概述4 概述概述
2、5 解析法解析法 数值数值分析方分析方法法 特征曲线法特征曲线法 剪切滑移破坏法剪切滑移破坏法概述概述6 解析法解析法是根据实际问题列是根据实际问题列出其出其平衡方程、几何方程和物平衡方程、几何方程和物理方程理方程,而后根据所给定的边,而后根据所给定的边界条件,对问题直接进行求解。界条件,对问题直接进行求解。 由于数学上的困难,由于数学上的困难,目前目前解析法还只能给出少数简单问解析法还只能给出少数简单问题的具体解答(如圆形隧道、题的具体解答(如圆形隧道、见第见第9 9章)章)。 5.3.1 5.3.1 解析法解析法7 数值分析方法的含义与类型数值分析方法的含义与类型 由于由于隧道结构隧道结构
3、大多大多几何形状复杂几何形状复杂,围岩围岩介质介质具有具有不均匀不均匀连续、各向异性等连续、各向异性等非线性特性非线性特性,而且,衬砌支护结构的计算,而且,衬砌支护结构的计算还与开挖方法、支护过程等有关。对于这类复杂问题,一般还与开挖方法、支护过程等有关。对于这类复杂问题,一般需要采取数值分析方法加以解决。需要采取数值分析方法加以解决。 数值分析方法:数值分析方法:是将一些可以用常微分方程、偏微分方是将一些可以用常微分方程、偏微分方程、积分方程、线性方程组或非线性方程组来描述、但用程、积分方程、线性方程组或非线性方程组来描述、但用纯纯解析方法又难以求解,或求解的过程非常繁杂、工作量巨大解析方法
4、又难以求解,或求解的过程非常繁杂、工作量巨大的工程问题,的工程问题,运用插值函数、函数逼近与数据拟合等方法运用插值函数、函数逼近与数据拟合等方法,并借助计算机的优势并借助计算机的优势用计算软件来进行用计算软件来进行的方法。的方法。5.3.2 5.3.2 数值分析方法数值分析方法8 数值分析方法的含义与类型数值分析方法的含义与类型 应用岩土工程中数值分析方法主要有以下主要类型:应用岩土工程中数值分析方法主要有以下主要类型: 有限元法有限元法 有限差分法有限差分法 边界元法边界元法 无界元法无界元法 有限元法与上述几种方法的耦合方法等有限元法与上述几种方法的耦合方法等 这里仅介绍有限元法这里仅介绍
5、有限元法5.3.2 5.3.2 数值分析方法数值分析方法9 有限单元法有限单元法(Finite Element MethodFinite Element MethodFEMFEM)目目前已成为隧道工程围岩稳定性分析和支护结构强度计算的前已成为隧道工程围岩稳定性分析和支护结构强度计算的有力工具。有力工具。 它把它把围岩和支护结构都划分为若干单元围岩和支护结构都划分为若干单元,然后根据,然后根据能量原理建立单元刚度矩阵,并形成整个系统的总刚度矩能量原理建立单元刚度矩阵,并形成整个系统的总刚度矩阵。从而阵。从而求出系统中各节点的位移和单元的应力求出系统中各节点的位移和单元的应力。 以下以以下以平面应
6、变问题平面应变问题来说明有限元解法的一般过程。来说明有限元解法的一般过程。 5.3.2 5.3.2 数值分析方法数值分析方法10 使用有限元法进行使用有限元法进行隧道工程问题的分析时,隧道工程问题的分析时,需要将所计算的区域需要将所计算的区域(即(即隧道及其周边一定范围的隧道及其周边一定范围的岩体)岩体)进行网格划分进行网格划分简称离散。简称离散。 (1 1)计算网格划分)计算网格划分离散离散 5.3.2 5.3.2 数值分析方法数值分析方法11 单元划分的疏密,大小和形状都会影响计算精度单元划分的疏密,大小和形状都会影响计算精度。 理论上讲,单元划分得越密越小,形状越规则,计算理论上讲,单元
7、划分得越密越小,形状越规则,计算精度越高。精度越高。 在实际工程中人们总是在实际工程中人们总是对计算范围中的某些区域更感对计算范围中的某些区域更感兴趣兴趣,如洞室或隧道结构物周围区域,地质构造区域等应力,如洞室或隧道结构物周围区域,地质构造区域等应力位移变化梯度大以及荷载有突变的区域。上述部位的单元划位移变化梯度大以及荷载有突变的区域。上述部位的单元划分可加密,而其它区域则可稀疏一些。分可加密,而其它区域则可稀疏一些。疏密区单元大小相差疏密区单元大小相差不宜过大不宜过大,应尽可能均匀过渡。,应尽可能均匀过渡。 (1 1)计算网格划分)计算网格划分离散离散5.3.2 5.3.2 数值分析方法数值
8、分析方法12 此外,根据此外,根据对称的特性对称的特性,在,在对称轴上的各点无垂直于对称轴方对称轴上的各点无垂直于对称轴方向的位移,因此可以从分析区域中向的位移,因此可以从分析区域中取出一半进行研究取出一半进行研究,这将大大减少,这将大大减少计算工作量和计算时间,计算工作量和计算时间,如右如右图。图。5.3.2 5.3.2 数值分析方法数值分析方法 (1 1)计算网格划分)计算网格划分离散离散 13 (2 2)计算网格划分)计算网格划分离散离散 某连拱隧道计算模型图某连拱隧道计算模型图14 (2 2)计算范围计算范围 大多数隧道工程都涉及无限域或半无限域,而有限元大多数隧道工程都涉及无限域或半
9、无限域,而有限元法处理这类问题通常是在有限区域里进行离散化。为了使这法处理这类问题通常是在有限区域里进行离散化。为了使这种处理方法不至于产生过大的误差,计算区域必须有足够的种处理方法不至于产生过大的误差,计算区域必须有足够的范围,并使区域外边界条件尽可能接近实际状态。范围,并使区域外边界条件尽可能接近实际状态。 理论分析表明,在均质弹性无限域中开挖的圆形洞室,理论分析表明,在均质弹性无限域中开挖的圆形洞室,由于荷载释放而引起的洞周介质应力和位移变化,由于荷载释放而引起的洞周介质应力和位移变化,在五倍洞在五倍洞径范围之外将小于,径范围之外将小于,三倍洞径之外约小于三倍洞径之外约小于。5.3.2
10、5.3.2 数值分析方法数值分析方法15 考虑工程的需要和有限考虑工程的需要和有限元离散误差以及计算误差,一元离散误差以及计算误差,一般选计算范围沿洞径般选计算范围沿洞径各方向均各方向均不小于不小于5 5倍洞径倍洞径。 对非圆形洞室或各向异对非圆形洞室或各向异性岩体材料中开挖的洞室,计性岩体材料中开挖的洞室,计算范围应适当扩大算范围应适当扩大 。 (2 2)计算范围计算范围 5.3.2 5.3.2 数值分析方法数值分析方法16 围岩:围岩:多多采用线性应变和二次采用线性应变和二次应变单元。通常应变单元。通常采用四节点或八节点采用四节点或八节点的四边形等参元的四边形等参元,它能适应曲边形的它能适
11、应曲边形的外形,便于进行网格自动划分,也具外形,便于进行网格自动划分,也具有较高的计算精度有较高的计算精度(3 3)单元类型选择)单元类型选择5.3.2 5.3.2 数值分析方法数值分析方法17 喷射混凝土喷射混凝土层或层或模注模注混凝土衬砌混凝土衬砌:多多都采用与围岩相同都采用与围岩相同的单元类型的单元类型 锚杆锚杆:最常用的是采用轴力最常用的是采用轴力杆单元杆单元来来模拟。模拟。(4 4)边界条件)边界条件5.3.2 5.3.2 数值分析方法数值分析方法 岩体边界上的条件通常两侧岩体边界上的条件通常两侧边界边界按水平方向固定按水平方向固定,铅直方向,铅直方向自由,下边界约束情况一般按自由,
12、下边界约束情况一般按铅铅直方向固定直方向固定,水平方向自由,水平方向自由。 无论采取何种边界条件,都无论采取何种边界条件,都可能会产生可能会产生与与实际情况不完全一实际情况不完全一致致的的误差。这种误差在靠近边界误差。这种误差在靠近边界处比远离边界处的误差大。这一处比远离边界处的误差大。这一现象称为现象称为“边界效应边界效应”。18 (1 1)隧道工程数值分析的特点隧道工程数值分析的特点 隧道施工过程主要包括洞室的开挖、喷射混凝土隧道施工过程主要包括洞室的开挖、喷射混凝土和锚杆的设置、二次衬砌混凝土结构的浇筑等等。和锚杆的设置、二次衬砌混凝土结构的浇筑等等。 这些施工过程都相当于在原始地应力场
13、中增加新这些施工过程都相当于在原始地应力场中增加新的荷载或改变地下结构的材料而产生二次、三次的荷载或改变地下结构的材料而产生二次、三次应力应力场。这是场。这是隧道工程数值分析的一个重要特点隧道工程数值分析的一个重要特点。 5.3.2 5.3.2 数值分析方法数值分析方法19 岩体在开挖隧道之前是处于一定的初始应力状态,岩体在开挖隧道之前是处于一定的初始应力状态,开开挖使隧道周边上各点的应力挖使隧道周边上各点的应力“解除解除”,从而引起围岩应力场,从而引起围岩应力场的变化。的变化。 在进行有限元分析时,必须设法模拟这个开挖卸荷的在进行有限元分析时,必须设法模拟这个开挖卸荷的效果。通用方法就是在隧
14、道周边的点上加效果。通用方法就是在隧道周边的点上加“等效释放荷载等效释放荷载”。 在实际模拟计算中,具体的在实际模拟计算中,具体的“等效释放荷载等效释放荷载”一般是由一般是由计算软件自动完成的计算软件自动完成的。 因隧道是采用多次支护,因隧道是采用多次支护,“等效释放荷载等效释放荷载”由围岩、初由围岩、初支、二衬共同承担,计算中则应根据计算经验人为设定它们支、二衬共同承担,计算中则应根据计算经验人为设定它们各自承担释放荷载的比例,即各自承担释放荷载的比例,即“荷载释放系数荷载释放系数”。 (2)(2)开挖过程的荷载释放开挖过程的荷载释放5.3.2 5.3.2 数值分析方法数值分析方法20(3
15、3)施工过程的模拟)施工过程的模拟5.3.2 5.3.2 数值分析方法数值分析方法 1 1)按照施工要求划分好开挖顺序。按照施工要求划分好开挖顺序。 2 2)按照隧道埋深的地质构造特点,进行开挖前的应力按照隧道埋深的地质构造特点,进行开挖前的应力分析,求出围岩中的初始地应力场和位移场,开挖前的应力分析,求出围岩中的初始地应力场和位移场,开挖前的应力状态可作为原始数据直接输入。状态可作为原始数据直接输入。 3 3)根据每次开挖的尺寸,去掉被开挖的单元,根据去根据每次开挖的尺寸,去掉被开挖的单元,根据去掉单元现时的应力值,求出被开挖出的自由表面各节点处由掉单元现时的应力值,求出被开挖出的自由表面各
16、节点处由这些单元作用的节点力。这些单元作用的节点力。将与这些节点力大小相等、方向相将与这些节点力大小相等、方向相反的力作用于自由表面相同的节点上,这些力就是等效开挖反的力作用于自由表面相同的节点上,这些力就是等效开挖释放荷载。释放荷载。21(3 3)施工过程的模拟)施工过程的模拟5.3.2 5.3.2 数值分析方法数值分析方法 ( 4 ( 4)在等效开挖释放荷载作用下进行分析,求出该开挖在等效开挖释放荷载作用下进行分析,求出该开挖步骤后围岩中的位移、应变、应力,并叠加于以前的状态上,步骤后围岩中的位移、应变、应力,并叠加于以前的状态上,重复以上步骤,直至最后一个开挖步骤结束。重复以上步骤,直至
17、最后一个开挖步骤结束。22 完成以上步骤后,就可运用完成以上步骤后,就可运用有限元法有限元法计算软件计算软件求算求算各各种种单元应力。单元应力。 根据数值分析计算结果,如何合理地判断围岩的稳定根据数值分析计算结果,如何合理地判断围岩的稳定性也是当前尚未解决的一个问题。在数值分析方法中,除性也是当前尚未解决的一个问题。在数值分析方法中,除非将支护结构离散为梁单元,否则都非将支护结构离散为梁单元,否则都只能求得支护结构中只能求得支护结构中的应力的应力,而不能直接采用规范中的公式校核支护结构的强,而不能直接采用规范中的公式校核支护结构的强度。度。 目前采用的判断围岩与支护结构稳定性的方法主要有目前采
18、用的判断围岩与支护结构稳定性的方法主要有如下几种:如下几种:5.3.2 5.3.2 数值分析方法数值分析方法23 (1 1)超载系数法超载系数法 将外荷载乘以系数值,并逐渐增大值进行反复计算,直将外荷载乘以系数值,并逐渐增大值进行反复计算,直到计算不能收敛为止,即认为围岩失稳,值为安全系数。到计算不能收敛为止,即认为围岩失稳,值为安全系数。 (2 2)材料安全储备法(强度折减法)材料安全储备法(强度折减法) 将材料的主要强度特征值,如、乘以值,逐渐降低值并将材料的主要强度特征值,如、乘以值,逐渐降低值并反复计算到围岩失稳(即计算不收敛)为止,就是安全度。反复计算到围岩失稳(即计算不收敛)为止,
19、就是安全度。 (3 3)经验类比法经验类比法 将计算所得洞壁变形值或塑性区范围与按经验所得的围将计算所得洞壁变形值或塑性区范围与按经验所得的围岩失稳时的允许位移值岩失稳时的允许位移值(极限位移值)或允许的塑性区(极限位移值)或允许的塑性区大小大小进行比较,由此确定围岩稳定性的安全度。进行比较,由此确定围岩稳定性的安全度。5.3.2 5.3.2 数值分析方法数值分析方法245.3.2 5.3.2 数值分析方法数值分析方法255.3.2 5.3.2 数值分析方法数值分析方法开挖部分支护后竖向位移图 开挖部分支护后竖向位移图 26FLAC3D 3.00Itasca Consulting Group,
20、 Inc.Minneapolis, MN USASettings: Model Perspective08:50:47 Thu Oct 27 2011Center: X: -1.679e+001 Y: 5.000e-001 Z: 1.906e+000Rotation: X: 0.000 Y: 0.000 Z: 0.000Dist: 4.706e+002Mag.: 3.81Ang.: 22.500Block GroupFLAC3D 3.00Itasca Consulting Group, Inc.Minneapolis, MN USASettings: Model Perspective08:5
21、2:27 Thu Oct 27 2011Center: X: -1.679e+001 Y: 5.000e-001 Z: 1.906e+000Rotation: X: 0.000 Y: 0.000 Z: 0.000Dist: 4.706e+002Mag.: 3.81Ang.: 22.500Block Group5.3.2 5.3.2 数值分析方法数值分析方法27 应力分布图应力分布图 位移分布图位移分布图 计算模型图计算模型图5.3.2 5.3.2 数值分析方法数值分析方法28 计算模型图计算模型图 塑性区图塑性区图 弯矩分布图弯矩分布图5.3.2 5.3.2 数值分析方法数值分析方法29 5.
22、3.3 5.3.3 剪切滑移法剪切滑移法305.3.3 5.3.3 剪切滑移法剪切滑移法315.3.3 5.3.3 剪切滑移法剪切滑移法325.3.3 5.3.3 剪切滑移法剪切滑移法33 124cos2tanexpabar5.3.3 5.3.3 剪切滑移法剪切滑移法345.3.3 5.3.3 剪切滑移法剪切滑移法wAstsapppppap35 ssssbdpsin2sssdcs43. 0030s5.3.3 5.3.3 剪切滑移法剪切滑移法 36 ststststbFpsin2stststFst0455.3.3 5.3.3 剪切滑移法剪切滑移法 37 et etFqAAAAFA5.3.3 5.
展开阅读全文