隧道工程第六讲隧道支护结构设计

southwestjIaotongunIversIty?隧道工程?课程第六讲隧道支护结构设计〔2〕1整理ppt内容提要1结构力学设计方法2其他设计方法3常用计算软件介绍隧道支护结构设计2整理ppt一结构力学设计方法3整理ppt隧道支护结构设计一、隧道工程的受力特点〔1〕荷载的模糊性〔2〕围岩物理力学参数难以准确获得〔3〕围岩压力承载体系不确定◆围岩不仅是荷载，同时又是承载体；

◆地层压力由围岩和支护结构共同承受；

◆充分发挥围岩自身承载力的重要性。〔4〕设计受施工方法和施作时机的影响很大〔5〕与地面结构受力的不同—围岩抗力的存在4整理ppt隧道支护结构设计二、隧道结构体系计算模型1.结构力学模型

特点：◆以支护结构作为承载主体；

◆围岩对支护结构的作用间接地表达为两点：

①围岩压力；②围岩弹性抗力。

◆采用结构力学方法计算。

适用于：模筑砼衬砌5整理ppt隧道支护结构设计二、隧道结构体系计算模型2.岩体力学模型

特点：◆支护结构与围岩视为一体，共同承受荷载，且以围岩作为承载主体；

◆支护结构约束围岩的变形；

◆采用岩体力学方法计算；

◆围岩表达为形变压力。

适用于：锚喷支护6整理ppt隧道支护结构设计三、结构力学方法计算原理它将支护结构和围岩分开来考虑:〔1〕支护结构是承载主体，围岩作为荷载的来源和支护结构的弹性支承〔2〕隧道支护与围岩的相互作用是通过弹性支承对支护结构施加约束来表达的〔3〕围岩的承载能力那么在确定围岩压力和弹性支承的约束能力时间接地考虑7整理ppt结构力学设计方法主动荷载模型假定弹性反力模型计算弹性反力模型常用计算模型8整理ppt结构力学设计方法

主动荷载模型适用于围岩与支护结构“刚度比〞较小，或饱和含水或用于初步设计9整理ppt结构力学设计方法主动荷载模型假定弹性反力模型计算弹性反力模型常用计算模型10整理ppt结构力学设计方法

假定弹性反力模型

几乎能适用于所有围岩类型，只不过抗力的大小和范围不同而已11整理ppt结构力学设计方法主动荷载模型假定弹性反力模型计算弹性反力模型常用计算模型12整理ppt结构力学设计方法

计算弹性反力模型

如弹性地基上的闭合框架、弹性支承法等13整理ppt结构力学设计方法1矩形框架结构〔地铁车站〕在?地下铁道?中讲述2装配式衬砌〔地铁区间〕在?水下隧道?中讲述3拱形结构〔山岭隧道〕本节重点讲述!与结构形式相适应的计算方法14整理ppt结构力学设计方法—拱形结构与结构形式相适应的计算方法半拱结构：不考虑弹性反力曲墙式衬砌：假定弹性反力直墙式衬砌：假定弹性反力+弹性地基梁15整理ppt结构力学设计方法计算模型：①荷载-结构模型②只有竖向荷载的作用，无侧向荷载的作用半拱形结构计算适用条件：①地质条件好，不需修边墙的山岭隧道；②大型落地拱结构,如飞机库；半拱形结构的适用条件及计算模型16整理ppt结构力学设计方法用结构轴线代替原衬砌的横断面,纵向长度取1m；拱脚弹性固定在围岩上，相当于弹性固定的无铰拱;因拱脚截面的剪力很小，而与围岩间摩擦较大，故径向位移为0，用径向刚性连杆表示;半拱形结构计算计算图式、根本结构及典型方程的建立17整理ppt结构力学设计方法半拱形结构计算由于半拱圈的拱矢和跨度的比值不大，在竖向荷载的作用下，结构为自由变形，无弹性抗力(脱离区);因此，半拱形结构为拱脚弹性固定的无铰拱〔考虑底部地基变形〕，为三次超静定结构。计算图式、根本结构及典型方程的建立18整理ppt结构力学设计方法半拱形结构计算典型方程的建立正对称的结构，作用有正对称的荷载，利用对称性，从拱顶切开，取根本结构如右图19整理ppt结构力学设计方法半拱形结构计算计算关键：拱顶单位位移和荷载位移的计算；拱脚位移的计算。20整理ppt结构力学设计方法1〕拱顶单位位移和荷载位移的计算：根据结构力学中位移计算方法，可求的某一点在单位力作用下，沿k方向的位移〔忽略剪力作用〕为：21整理ppt结构力学设计方法1〕拱顶单位位移和荷载位移的计算：将X1〔弯矩〕，X2〔轴力〕，X3〔剪力，取零〕以及外荷载作用下结构各截面内力代入可得：22整理ppt结构力学设计方法1〕拱顶单位位移和荷载位移的计算：单位位移和荷载位移的计算分两种情况：〔1〕当半拱为一段规那么圆弧时，可直接积分计算；〔2〕当半拱为多段圆弧时，将圆弧衬砌分块，采用辛普生求和公式；23整理ppt结构力学设计方法典型方程半拱形结构计算24整理ppt结构力学设计方法2〕拱脚位移的计算：根本思路：求出拱脚支承面在单位力〔力矩、竖直力、水平力〕作用下位移〔转角、水平位移〕；求出赘余力及外荷载作用下的拱脚反力；由拱脚反力与各单位力位移计算拱脚位移25整理ppt结构力学设计方法2〕拱脚位移的计算：26整理ppt结构力学设计方法2〕拱脚位移的计算：27整理ppt结构力学设计方法2〕拱脚位移的计算：28整理ppt结构力学设计方法2〕拱脚位移的计算：29整理ppt结构力学设计方法2〕拱脚位移的计算：〔b〕拱脚支承反力的计算30整理ppt结构力学设计方法31整理ppt结构力学设计方法2〕拱脚位移的计算：〔c〕拱脚位移计算32整理ppt结构力学设计方法3〕衬砌截面内力的计算：33整理ppt结构力学设计方法3〕衬砌截面内力的计算：34整理ppt结构力学设计方法曲墙拱形结构计算①计算图示的建立用结构的轴线代替原衬砌断面;在计算中将拱圈和边墙作为一个整体，把它看成是一个支承在弹性围岩上的高拱结构；因底部摩擦力很大，无水平位移；衬砌在以竖向压力为主的主动荷载作用下，拱圈的顶部发生向坑道内的变形不受围岩的约束，形成“脱离区〞。衬砌结构的侧面局部那么压向围岩，形成“弹性反力区〞，引起相应的弹性反力。计算原理35整理ppt结构力学设计方法曲墙拱形结构计算计算原理弹性反力按苏联布加耶娃法假定分布图形为月牙形分布，分布图形用3个特征点控制上零点b〔即脱离区的边界〕与对称轴线间的夹角一般采用下零点a取在墙底，因该处无水平位移；最大弹性反力点h可假定在衬砌最大跨度处。36整理ppt结构力学设计方法曲墙拱形结构计算计算原理最大弹性反力点以下围岩弹性反力对于衬砌的变形还会在围岩与衬砌间产生相应的摩擦力分布规律最大弹性反力点以上37整理ppt结构力学设计方法曲墙拱形结构计算计算原理弹性地基上的高拱荷载：围岩压力弹性反力基底位移〔约束〕衬砌与围岩的摩擦力38整理ppt结构力学设计方法曲墙拱形结构计算计算思路在结构与荷载均为对称的条件下，可以从拱顶切开，以一对悬臂曲梁作为根本结构，切开处赘余力为X1、X2，剪力X3＝0。在主动荷载和弹性反力作用下，根据拱顶相对转角及相对水平位移为零的条件，可以得到2个典型方程式；X1X2基本结构39整理ppt结构力学设计方法利用h点变形协调条件来增加1个方程式。根据叠加原理，h点的最终位移即为:而h点的位移与该点的弹性反力存在下述关系：由以上可以推导出：40整理ppt结构力学设计方法因此，综上两点可得：X1X2基本结构计算关键：拱顶单位位移和荷载位移的计算；墙脚位移的计算；最大跨度处荷载位移的计算41整理ppt结构力学设计方法曲墙拱形结构计算1.求主动荷载作用下的衬砌内力42整理ppt结构力学设计方法曲墙拱形结构计算1.求主动荷载作用下的衬砌内力43整理ppt结构力学设计方法曲墙拱形结构计算2.求单位弹性反力作用下的衬砌内力44整理ppt结构力学设计方法曲墙拱形结构计算3.位移及最大弹性反力值的计算①主动荷载作用下最大抗力点h点位移的计算外荷载h点的位移＝基底转动引起的h点的位移＋受力引起的h点的位移基底转动引起的h点的位移：45整理ppt结构力学设计方法3.位移及最大弹性反力值的计算受力引起的h点的位移：那么外荷载作用下，h点的位移为：46整理ppt结构力学设计方法3.位移及最大弹性反力值的计算③最大弹性反力的计算47整理ppt结构力学设计方法4.衬砌内力计算及校核计算结果的正确性①利用叠加原理可以求出任意截面最终的内力值曲墙拱形结构计算②拱脚截面最终转角48整理ppt结构力学设计方法4.衬砌内力计算及校核计算结果的正确性曲墙拱形结构计算③按变形协调条件，校核整个计算过程：49整理ppt结构力学设计方法直墙拱形结构计算①结构拱圈支承在弹性地基梁上的弹性固定无铰拱；

边墙双向弹性地基梁

计算原理50整理ppt结构力学设计方法直墙拱形结构计算计算原理②弹性反力拱圈：任意截面弹性反力荷载图形假设为二次抛物线，作用方向为径向；计算公式如下；边墙：用弹性地基梁的方法计算51整理ppt结构力学设计方法直墙拱形结构计算计算原理③附加一个方程：墙顶变位④拱圈内力的计算：在原理上与弹性固定的高拱结构完全相同，只是计及墙顶变位⑤边墙：作为弹性地基上的直梁来计算52整理ppt结构力学设计方法直墙拱形结构计算计算原理弹性地基梁，按其换算长度l的不同，可分为3种情况：①长梁l≥2.75②短梁1<l<2.75③刚性梁1≥ll为梁的长度〔即边墙高度〕，为弹性地基梁的弹性特征值边墙53整理ppt结构力学设计方法1.弹性地基梁在梁端荷载作用下的梁端位移计算〔仍然采用叠加原理〕直墙拱形结构计算在弹性地基梁的c端，作用有拱脚传来的外力Mc和Hc，要求c端的位移βc和Uc54整理ppt结构力学设计方法2.拱圈内力计算直墙拱形结构计算55整理ppt结构力学设计方法2.拱圈内力计算拱圈在拱部单位弹性反力图作用下的计算公式与主动荷载的情况相似；从上述结果可求出拱部在主动外荷载和单位弹性反力作用下最后的内力。此时墙顶位移利用叠加原理，并附加方程。可解出墙顶的所有参数。56整理ppt结构力学设计方法3.边墙内力和位移计算墙端初参数：利用弹性地基梁的初参数公式求得边墙各截面的内力和位移:

直墙拱形结构计算57整理ppt结构力学设计方法3.边墙内力和位移计算边墙为短梁时，距墙顶为x的任一截面的内力和位移的计算公式如下：

58整理ppt结构力学设计方法3.边墙内力和位移计算——拱脚〔墙顶〕最终位移值，根据地基局部变形理论求得边墙各截面的抗力为59整理ppt结构力学设计方法3.边墙内力和位移计算边墙为长梁时，距墙顶为x的任一截面的内力和位移的计算公式为：60整理ppt岩体力学〔地层-结构〕设计方法根底：弹性理论、弹塑性理论——解析解这里仅介绍弹性状态下的平面轴对称问题，即初始应力为轴对称分布的圆形隧道问题，围岩视为无重平面，初始应力作用在无穷远处，并假定支护结构与围岩密贴，即外径r0与围岩开挖半径相等，且与开挖同时瞬间完成。下面以均匀内压水工隧洞的计算为例，说明解析法计算的根本思路。61整理ppt

水工隧洞衬砌的材料主要有混凝土、钢筋混凝土和锚喷支护等。厚度一般在20cm以上，故力学分析中可将其视为厚壁圆筒。在均匀内水压力作用下，厚壁圆筒的内力分析是轴对称问题。衬砌应力分析岩体力学〔地层-结构〕设计方法62整理ppt将作用于衬砌内外表的水压力记为pw，地层对衬砌外外表作用的形变压力记为pa。在pw作用下，圆环将向外扩张。设衬砌在半径r处由pw引起的径向位移为u，那么该处的圆周长度必从2πr增加到2π(r+u)，由此可得到衬砌的切向应变为衬砌应力分析岩体力学〔地层-结构〕设计方法63整理ppt

图中示有从衬砌圆环中取出的单元体。因r的增量为u，故单元体一边的长度dr的增量可记为du，该边长度可记为。由此可得衬砌的径向应变为衬砌应力分析〔6.1.2〕岩体力学〔地层-结构〕设计方法64整理ppt衬砌材料的弹性常数为，记，并依据习惯近似按平面应变问题分析衬砌，那么由平面问题极坐标解的物理方程可写为衬砌应力分析可解得:〔6.1.4〕〔6.1.3〕岩体力学〔地层-结构〕设计方法65整理ppt

因作用在单元体上的外荷载为零，且在轴对称情况下单元体内力分量中的剪应力也为零，故根据平面问题极坐标解的静力平衡方程式，得出方程衬砌应力分析或写成：〔6.1.5〕岩体力学〔地层-结构〕设计方法66整理ppt将式〔6.1.4〕代入式〔6.1.5〕，可得衬砌应力分析〔6.1.6〕〔6.1.7〕〔6.1.8〕岩体力学〔地层-结构〕设计方法67整理ppt衬砌应力分析〔6.1.9〕〔6.1.10〕由边界条件可得：岩体力学〔地层-结构〕设计方法68整理ppt将求得的和代入式〔6.1.9〕、〔6.1.10〕，经整理可得衬砌内的应力衬砌应力分析由式〔6.1.8〕可写出衬砌内任一点的径向位移为〔6.1.11〕式中，岩体力学〔地层-结构〕设计方法69整理ppt分析均匀内力圆形水工隧洞围岩的应力仍可采用厚壁圆筒原理。围岩的弹性常数为，并记，那么由式〔6.1.9〕、〔6.1.10〕可写出围岩应力的表达式为

洞室围岩应力分析〔6.1.13〕岩体力学〔地层-结构〕设计方法70整理ppt如下图，洞室围岩的应力边界条件为

洞室围岩应力分析〔6.1.14〕，将代入式〔6.1.13〕，可得将式〔6.1.13〕改为：〔6.1.15〕岩体力学〔地层-结构〕设计方法71整理ppt

洞室围岩应力分析〔6.1.16〕将代入式〔6.1.14〕，可得将上式代入式〔6.1.14〕和〔6.1.15〕，即得围岩内的应力〔6.1.17〕仿照式〔6.1.8〕可写出围岩径向位移的计算式为岩体力学〔地层-结构〕设计方法72整理ppt由式〔6.1.16〕可知，内水压力使围岩产生的切向应力是拉应力。设假设的量值大于围岩中原来存在的压应力，且差值超过岩体的抗拉强度，那么当衬砌抗拉强度缺乏时岩体将与衬砌一起发生开裂。某些有压水工隧洞出现新的、平行于洞轴线且沿圆周均匀间隔分布的裂缝，原因就在于围岩在环向出现了较大的拉应力。

洞室围岩应力分析岩体力学〔地层-结构〕设计方法73整理ppt将式〔6.1.16〕中的r0理解为毛洞半径，pa理解为内压力，那么该式就成为无衬砌圆形水工隧洞围岩应力的计算式。显而易见，环向拉应力的存在必然对无衬砌水工隧洞的适用性起限止作用。故工程设计中常需设置衬砌或锚喷支护，使支护和围岩共同承受内水压力。

洞室围岩应力分析岩体力学〔地层-结构〕设计方法74整理ppt

假设在内压力pw作用下隧洞衬砌对围岩产生的作用力如图，为

衬砌与围岩共同作用的计算那么由式〔6.1.14〕可写出：由此可得：〔6.1.18〕岩体力学〔地层-结构〕设计方法75整理ppt将上式代入式〔6.1.17〕，可得在r＝r0处围岩的径向位移为

衬砌与围岩共同作用的计算〔6.1.19〕由式〔6.1.9〕可知，在衬砌内外表有关系式在衬砌外外表有关系式岩体力学〔地层-结构〕设计方法76整理ppt

将以上2式联立，可解得

衬砌与围岩共同作用的计算由此可写出衬砌外缘的径向位移为岩体力学〔地层-结构〕设计方法77整理ppt因在r＝r0处衬砌与围岩的径向位移应相等，故由式〔6.1.19〕、〔6.1.20〕可得

衬砌与围岩共同作用的计算〔6.1.21〕上式经整理，即得求得了值以后，由式〔6.1.11〕、〔6.1.16〕即可算出衬砌与围岩的应力，并可据以验算围岩的稳定性及进行衬砌截面的设计。岩体力学〔地层-结构〕设计方法78整理ppt数值计算方法1、数值计算的开展主要于70年代中期，随着计算机的广泛应用和计算技术的不断成熟而逐渐开展起来的一种地下结构计算方法。数值解法包括有限差分法、有限元法、边界元法、离散元法、颗粒流法等，其中有限元和边界元建立在连续介质力学的根底上，适合于小变形分析，是开展较早、较成熟的方法，尤以有限元在地下工程中应用更为广泛。有限元法不但可以模拟各种施工过程和各种支护效果，同时可以考虑复杂地层情况〔如断层、节理等地质构造以及地下水等〕，对材料的各种非线性也可考虑。一、概述

79整理ppt2、地下工程有限元法特点支护结构与岩体共同作用，统一的组合体来考虑，共同离散。作用在岩体上的荷载是地应力。深埋，简化为均布垂直地应力和水平地应力。浅埋，垂直应力和侧压系数均按自重应力场确定。支护结构视作线弹性，围岩的应力-应变关系视作非线性可以选用弹塑性、粘塑性、粘弹塑性的力学模型开挖及支护导致应力重分布，能够考虑开挖与支护施工步骤的影响。地下结构工程一般轴线很长，平面应变问题，从而使计算大大简化。一、概述

数值计算方法80整理ppt数值计算方法二、有限元原理

隧道围岩与支护结构的几何模型模型离散为由各种单元组成的计算模型；有限元中分析的结构已不是原有的物体或结构物，而是由众多单元通过节点相连的离散物体。用有限元分析计算所获得的结果只是近似的。如果划分单元数目非常多而又合理，那么所获得的结果就与实际情况相符合。81整理ppt数值计算方法三、单元类型

按单元形状分：线单元、面单元、体单元单元的阶次分为线性单元、二次单元、高次单元。合理选择单元类型对问题的正确解决至关重要，而单元类型的选择又与要解决问题的的本身存在密切关系。在选择单元类型前，首先要对问题本身有非常明确的认识，然后对每一种单元的特性都要详细了解，才能选择适宜的单元。82整理ppt数值计算方法三、单元类型

阶次选择需要在计算精度和计算规模间综合衡量；线性单元的严重扭曲变形可能引起计算精度下降，而高阶次的单元对变形扭曲不敏感；有曲边或曲面存在时，通常推荐使用高阶单元，以获得较高的曲面精度；对于非线性问题，高阶单元并不比线性单元更有效；单元阶次对求解精度影响，相对平面单元和三维实体单元之间简化的差异来说影响要小得多；高阶单元允许忽略局部或所有边的中节点，，此时高阶单元那么转换为线性单元，计算精度也随之降低。83整理ppt数值计算方法三、单元类型

84整理ppt数值计算方法三、单元类型

85整理ppt数值计算方法四、计算步骤

结构力学方法86整理ppt围岩强度判据的理论根底是强度破坏理论，其中尤以剪切屈服破坏为主。如应用比较广泛的摩尔-库伦和D-P准那么。摩尔-库伦准那么：五、屈服准那么数值计算方法87整理ppt摩尔-库伦准那么：五、屈服准那么数值计算方法88整理pptD-P准那么：该准那么假定物体内单位体积的应变能到达某一极限值时所对应的点开始屈服，并考虑到静水压力的作用，对于平面应变状态：

五、屈服准那么数值计算方法89整理ppt六、计算范围的选取三维问题，平面应变问题半无限空间问题，半无限平面问题依据圣维南原理，选取有限计算区域，在满足精度要求的根底上，节省计算费用。计算边界即可确定在3～5倍开挖宽度。根据对称性的特点，分析区域可以取1/2或1/4。数值计算方法90整理ppt七、分析过程问题及求解域定义；求解域离散化确定状态变量及控制方法单元推导总装联立方程组求解和结果解释有限元分析可分成三个阶段，前置处理、计算求解和后置处理。数值计算方法91整理ppt特征曲线设计方法1、围岩特征曲线的根本概念：

说明围岩在洞室周边所需提供的支护阻力与其周边位移的关系：在洞周极限位移范围内，容许围岩的位移增加，所需要的支护阻力减小，而应力重分布的结果大局部由围岩承担。一、围岩特征曲线

92整理ppt2、不同应力状态下的围岩特征曲线：〔1〕、弹性应力状态下的围岩特征曲线〔2〕、弹塑性应力状态下的围岩特征曲线①、不考虑塑性区体积扩容的方程②、考虑体积扩容的方程一、围岩特征曲线

特征曲线设计方法93整理ppt弹性应力状态下的围岩特征曲线一、围岩特征曲线

当pa=0时，uamax=(1+μ)/E*σyr0；当ur0=0时，pamax=σy，即欲坑道周边不产生位移，就需要有相当于初始地应力大小的支护阻力，是坑道处于初始地应力状态，显然这是不可能的。特征曲线设计方法94整理ppt一、围岩特征曲线

由于塑性区c、φ

值是变化的，代以不同的c、φ值就可以得到不同的收敛曲线。通常采用平均的c、φ

值来确定收敛线。特征曲线设计方法95整理ppt一、围岩特征曲线

弹塑性应力状态下围岩特征曲线-考虑塑性区体积扩容式中：M为弹塑性边界上应力差。当μ=0.5时，那么上式成为不考虑塑性区体积扩容的方程。特征曲线设计方法96整理ppt1、支护特征曲线的根本概念：指作用在支护上的荷载与支护变形的关系曲线，支护结构所能提供的支护阻力随着支护结构刚度增大而增大。支护结构的刚度和支护与围岩的接触状态有关。二、支护结构特征曲线

围岩初始位移最大支护阻力特征曲线设计方法97整理ppt2、不同类型支护结构的特征曲线：〔1〕混凝土或喷混凝土的支护特征曲线二、支护结构特征曲线

〔σcc为喷混凝土或混凝土的单轴抗压强度〕上式是按照外压作用下的圆管理论推导得来的，它只能用于圆形衬砌，而且超挖值也受到限制。特征曲线设计方法98整理ppt二、支护结构特征曲线

tc<0.04a时，支护阻力可以按照上式计算〔薄壁圆管计算公式〕，式中σcc为喷混凝土单轴抗压强度tc<0.04a龄期28天喷混凝土厚度对特征曲线的影响特征曲线设计方法99整理ppt2、不同类型支护结构的特征曲线：〔2〕钢支撑的支护特征曲线二、支护结构特征曲线

钢支撑通常是在工厂预制的,因此在工地安装时,必须采用不同厚度的木垫块使支撑与围岩密贴,甚至给予一定的预加荷载。这样在绘制钢支撑的特性曲线时,就不仅要考虑支撑本身的构造,而且要考虑垫块的刚度。楔点钢支撑模式特征曲线设计方法100整理ppt2、不同类型支护结构的特征曲线：〔2〕钢支撑的支护特征曲线木模块刚度对钢支撑特征曲线的影响s-沿隧道极度方向的支撑间距；θ-楔点间夹角之半〔弧度〕；w-支撑翼缘宽度；As-钢支撑截面积；Is-支撑截面惯性矩；Es-钢的弹性模量；tB-垫块厚度；EB-垫块材料的弹性模量;σys-钢的屈服强度;X-钢支撑截面厚度。特征曲线设计方法101整理ppt2、不同类型支护结构的特征曲线：〔3〕灌浆锚杆的支护特征曲线二、支护结构特征曲线

目前，灌浆锚杆的支护机理通过拉拔试验进行一些定性讨论。通过拉拔试验得出的拉拔荷载与位移的关系如右图所示,是非线性的,它说明拉拔荷载与锚杆本身的强度、直径、长度以及使围岩与锚杆胶结在一起的材料〈砂浆或树脂〉的强度等有关。拉拔位移曲线特征曲线设计方法102整理ppt三、特征曲线设计方法

特征曲线法的根本原理是利用围岩特征曲线和支护结构特征曲线交会的方法来决定支护体系的最正确平衡条件。支护结构特征曲线与围岩特征曲线交点处的横坐标为形成平衡体系时洞周发生的位移。交点纵坐标以下的局部为支护结构上承受的荷载，以上局部由围岩来承担。刚度大的支护结构承受较大的围岩作用力〔压力〕；反之，柔性较好的支护结构所承受的围岩压力要小得多。所以，我们在工程中强调采用柔性支护以节约本钱，但它也应有必要的刚度，以便有效地控制围岩变形，而到达稳定。特征曲线设计方法103整理ppt(①—锚杆支护曲线；②、④—喷混凝土支护曲线；③—组合结构支护曲线)〔1〕不同刚度的支护结构与围岩达成平衡时的pa和ur0是不同的图中，锚杆的支护结构特征曲线①没有能和围岩特征曲线相交，说明锚杆的刚度太小，它所能提供的支护阻力满足不了围岩稳定性的需求，这种供不应求的状况最终将导致围岩失稳。当然，增加支护结构的刚度并不意味着要增加支护结构的尺寸和数量，重要的是支护结构及早地形成闭合断面。特征曲线设计方法104整理ppt(①—锚杆支护曲线；②、④—喷混凝土支护曲线；③—组合结构支护曲线)〔2〕同样刚度的支护结构，假设时间不同，最后达成平衡状态也不同如图中曲线②和④。支护结构架设的越早，它所承受的围岩压力就越大。但这不等于说支护结构参与相互作用的时间愈迟愈好，因为初始变形不加控制会导致围岩迅速松弛而崩塌。当然，这个范围的大小视围岩的特性和埋置深度而变。特征曲线设计方法105整理ppt隧道工程结构计算软件FastLagrangianAnalysisofContinua美国Itasca咨询公司开发2D程序(1986)1990年代初引入中国有限差分法(FDM)DOS版→2.0→2.1→3.0→5.0→6.0Itasca其他软件1.FLAC3D106整理ppt隧道工程结构计算软件1.FLAC3D应用：岩土力学分析，例矿体滑坡、煤矿开采沉陷预测、水利枢纽岩体稳定性分析、采矿巷道稳定性研究等岩土工程、采矿工程、水利工程、地质工程特色：大应变模拟完全动态运动方程使得FLAC3D在模拟物理上的不稳定过程不存在数值上的障碍显示求解具有较快的非线性求解速度107整理ppt隧道工程结构计算软件根本原理有限差分法Lagrangian网格空间混合离散技术Lagrangian格式动量平衡方程FLAC3D的求解过程FLAC3D的本构模型108整理ppt隧道工程结构计算软件求解过程109整理ppt隧道工程结构计算软件本构模型开挖模型null3个弹性模型各向同性弹性横观各向同性弹性正交各向同性弹性8个塑性模型〔Drucker-Prager模型、Morh-Coulomb模型、应变硬化/软化模型、遍布节理模型、双线性应变硬化/软化遍布节理模型、修正剑桥模型和胡克布朗模型〕110整理ppt隧道工程结构计算软件111整理ppt隧道工程结构计算软件FLAC3D前后处理命令栏112整理ppt隧道工程结构计算软件FLAC3D前后处理113整理ppt隧道工程结构计算软件一个简单算例genzonbrisize333;建立网格modelelas;材料参数propbulk3e8shear1e8inidens2000;初始条件fixzranz-.1.1;边界条件fixxranx-.1.1fixxranx2.93.1fixyrany-.1.1fixyrany2.93.1setgrav00-10solve;求解appnstr-10e4ranz3x12y12solveRUNFLAC3D114整理ppt隧道工程结构计算软件前后处理优缺点多种zone类型后处理快捷、方便、丰富计算过程中的hist变量动态显示FISH可进行参数化模型设计单元状态的可编程计算暂停时的后处理与可保存复杂模型的建模功能不强可以编程导入其他软件形成的网格(比方：Ansys、Adina、GeoCAD)无等值线的后处理功能(3D)可编程将.sav文件写入TecPlot等其他后处理软件全命令操作，学习困难鼠标功能单一(双击取击点坐标)优点：缺点：115整理ppt隧道工程结构计算软件2.ANSYSANSYS是世界上著名的大型通用有限元计算软件，它包括热、电、磁、流体和结构等诸多模块，具有强大的求解器和前、后处理功能，为我们解决复杂、庞大的工程工程和致力于高水平的科研攻关提供了一个优良的工作环境，更使我们从繁琐、单调的常规有限元编程中解脱出来。ANSYS本身不仅具有较为完善的分析功能，同时也为用户自己进行二次开发提供了友好的开发环境。ANSYS程序自身有着较为强大三维建模能力，仅靠ANSYS的GUI(图形界面)就可建立各种复杂的几何模型；此外，ANSYS还提供较为灵活的图形接口及数据接口。因而，利用这些功能，可以实现不同分析软件之间的模型转换。116整理ppt隧道工程结构计算软件主要功能1.结构分析1〕静力分析-用于静态载荷。可以考虑结构的线性及非线性行为。●线性结构静力分析●非线性结构静力分析♦几何非线性：大变形、大应变、应力强化、旋转软化♦材料非线性：塑性、粘弹性、粘塑性、超弹性、多线性弹性、蠕变、肿胀等♦接触非线性：面面/点面/点点接触、柔体/柔体刚体接触、热接触♦单元非线性：死/活单元、钢筋混凝土单元、非线性阻尼/弹簧元、预紧力单元等117整理ppt隧道工程结构计算软件2〕模态分析-计算线性结构的自振频率及振形.谱分析是模态分析的扩展，用于计算由于随机振动引起的结构应力和应变(也叫作响应谱或PSD).3〕谐响应分析-确定线性结构对随时间按正弦曲线变化的载荷的响应.4〕瞬态动力学分析-确定结构对随时间任意变化的载荷的响应.可以考虑与静力分析相同的结构非线性行为.5〕谱分析6〕随机振动分析等7〕特征屈曲分析-用于计算线性屈曲载荷并确定屈曲模态形状.(结合瞬态动力学分析可以实现非线性屈曲分析.)8〕专项分析: