隧道力学 第5讲(中南大学 施成华)

1、中南大学土木工程学院 隧道工程系第5讲 岩体力学方法本讲主要内容：1、 解析法是根据所给定的边界解析法是根据所给定的边界条件，对问题的平衡方程、几何条件，对问题的平衡方程、几何方程和物理方程直接求解，而后方程和物理方程直接求解，而后根据所给定的边界条件，对问题根据所给定的边界条件，对问题直接进行求解。由于数学上的困直接进行求解。由于数学上的困难，目前解析法还只能给出少数难，目前解析法还只能给出少数简单问题的具体解答。简单问题的具体解答。 土木工程学院隧道工程系土木工程学院隧道工程系施成华1 1、计算范围的确定和离散方法、计算范围的确定和离散方法（1 1）计算范围计算范围( (scope of

2、problem)scope of problem) 大多数隧道工程都涉及无限域或半无限域，而有限元法处大多数隧道工程都涉及无限域或半无限域，而有限元法处理这类问题通常是在有限区域里进行离散化。为了使这种处理理这类问题通常是在有限区域里进行离散化。为了使这种处理方法不至于产生过大的误差，离散区域必须有足够的范围，并方法不至于产生过大的误差，离散区域必须有足够的范围，并使区域外边界条件尽可能接近实际状态。理论分析表明，在均使区域外边界条件尽可能接近实际状态。理论分析表明，在均质弹性无限域中开挖的圆形洞室，由于荷载释放而引起的洞周质弹性无限域中开挖的圆形洞室，由于荷载释放而引起的洞周介质应力和位移变

3、化，在五倍洞径范围之外将小于，介质应力和位移变化，在五倍洞径范围之外将小于，三倍三倍洞径之外约小于洞径之外约小于。土木工程学院隧道工程系土木工程学院隧道工程系施成华 考虑工程的需要和有限元离散误差以及计算误差，一般考虑工程的需要和有限元离散误差以及计算误差，一般选计算范围沿洞径选计算范围沿洞径各方向均不小于倍洞径各方向均不小于倍洞径为好。但计为好。但计算实践表明，对非圆形洞室或各向异性岩体材料中开挖的洞算实践表明，对非圆形洞室或各向异性岩体材料中开挖的洞室，则计算范围应适当扩大或取上限尺寸。如果只考虑自重室，则计算范围应适当扩大或取上限尺寸。如果只考虑自重应力场，则应力场，则可借助于无限域单元

4、可借助于无限域单元，免去计算范围选取的麻烦，免去计算范围选取的麻烦，但是无限元和有限元的交接位置的确定仍要考虑上述原则，但是无限元和有限元的交接位置的确定仍要考虑上述原则，只是范围可略小一些或取下限。只是范围可略小一些或取下限。土木工程学院隧道工程系土木工程学院隧道工程系施成华使用有限元法进行隧道工程分析，在计算范围确定使用有限元法进行隧道工程分析，在计算范围确定之后并非任何一种离散形式都可以得到同样的结果。之后并非任何一种离散形式都可以得到同样的结果。单单元划分的疏密，大小和形状都会影响计算精度元划分的疏密，大小和形状都会影响计算精度。理论上。理论上讲，单元划分得越密越小，形状越规则，计算精

5、度越高。讲，单元划分得越密越小，形状越规则，计算精度越高。据误差分析，据误差分析，应力误差与单元尺寸一次方成正比应力误差与单元尺寸一次方成正比，位移位移误差与单元尺寸二次方成正比误差与单元尺寸二次方成正比。土木工程学院隧道工程系土木工程学院隧道工程系施成华 但在实际工程中人们总是但在实际工程中人们总是对计算范围中的某些区域更对计算范围中的某些区域更感兴趣感兴趣，如洞室或隧道结构物周围区域，地质构造区域等，如洞室或隧道结构物周围区域，地质构造区域等应力位移变化梯度大以及荷载有突变的区域。上述部位的应力位移变化梯度大以及荷载有突变的区域。上述部位的单元划分可加密，而其它区域则可稀疏一些。单元划分可

6、加密，而其它区域则可稀疏一些。疏密区单元疏密区单元大小相差不宜过大大小相差不宜过大，应尽可能均匀过渡。有人研究认为，应尽可能均匀过渡。有人研究认为边边缘区域单元尺寸可为稠密区域单元的倍缘区域单元尺寸可为稠密区域单元的倍。土木工程学院隧道工程系土木工程学院隧道工程系施成华 单元形式可采用单元形式可采用三节点三角形常应变元三节点三角形常应变元，六节点三六节点三角形变应变元角形变应变元，四节点四边形四节点四边形和和八节点四边形八节点四边形等参元等。等参元等。对三维问题则常用对三维问题则常用八至二十节点六面体单元和壳单元八至二十节点六面体单元和壳单元等。等。三角形单元的优点是对复杂的几何形状适应性强。

7、在洞三角形单元的优点是对复杂的几何形状适应性强。在洞室周围附近区域的应力变化较大，采用较室周围附近区域的应力变化较大，采用较密的三角形单密的三角形单元元往往比多节点四边形单元取得更高的计算精度。三角往往比多节点四边形单元取得更高的计算精度。三角形单元形成的整体刚度阵带宽较小。形单元形成的整体刚度阵带宽较小。土木工程学院隧道工程系土木工程学院隧道工程系施成华 在非线性分析中，由于单元个数多，更能反映材料弹在非线性分析中，由于单元个数多，更能反映材料弹模和泊松比的非线性变化。缺点是三角形单元应力波动大，模和泊松比的非线性变化。缺点是三角形单元应力波动大，相邻单元应力往往不连续相邻单元应力往往不连续

8、。四边形单元的优点是能够较好。四边形单元的优点是能够较好地反映应力变化。当节点数相同时精度高于三角形单元，地反映应力变化。当节点数相同时精度高于三角形单元，且在边界较规则时用四边形单元较为简单。如果程序许可，且在边界较规则时用四边形单元较为简单。如果程序许可，也可以混合使用三角形与四边形单元，但公共边上位移必也可以混合使用三角形与四边形单元，但公共边上位移必须协调。在离散计算区域时还需注意以下几方面的问题。须协调。在离散计算区域时还需注意以下几方面的问题。 土木工程学院隧道工程系土木工程学院隧道工程系施成华1 1）一个单元各边长相差不能过大，两边夹角不能过小，一个单元各边长相差不能过大，两边夹

9、角不能过小，各夹角最好尽量相等。各夹角最好尽量相等。2 2）一个单元中不能包含两种或两种以上的材料。）一个单元中不能包含两种或两种以上的材料。3 3）集中荷载作用点或荷载突变处必须布置节点。）集中荷载作用点或荷载突变处必须布置节点。4 4）如隧道结构和岩体结构具有对称性时取部分计算范围）如隧道结构和岩体结构具有对称性时取部分计算范围进行离散。几何形状和材料特性方面都具有对称性进行离散。几何形状和材料特性方面都具有对称性时，可利用该对称性取部分计算范围进行剖分。时，可利用该对称性取部分计算范围进行剖分。土木工程学院隧道工程系土木工程学院隧道工程系施成华5 5）洞室边缘两侧的对应单元，其大小形状尽

10、量一致。）洞室边缘两侧的对应单元，其大小形状尽量一致。6)6)洞室边缘及附近单元的布置应考虑设置锚杆的方向及洞室边缘及附近单元的布置应考虑设置锚杆的方向及深度，以便施加锚固力。深度，以便施加锚固力。7)7)洞室内单元的划分要考虑到分部开挖的分界线和部分洞室内单元的划分要考虑到分部开挖的分界线和部分开挖区域的分界线。开挖区域的分界线。8 8）计算范围内的单元划分还要考虑到地下水位的变化分）计算范围内的单元划分还要考虑到地下水位的变化分界面。界面。 土木工程学院隧道工程系土木工程学院隧道工程系施成华 计算范围的外边界可采取两种方式处理；其一计算范围的外边界可采取两种方式处理；其一为为位移边界条件位

11、移边界条件，即一般假定边界点位移为零（也，即一般假定边界点位移为零（也有假定为弹性支座或给定位移的，但地下工程分析有假定为弹性支座或给定位移的，但地下工程分析中很少用）。其二是假定为中很少用）。其二是假定为力边界条件力边界条件，包括自由，包括自由边界边界( (P=P=) )条件。还可以给定混合边界条件，即节条件。还可以给定混合边界条件，即节点的一个自由度给定位移，另一个自由度给定节点点的一个自由度给定位移，另一个自由度给定节点力（二维问题）。当然无论哪种处理都有一定的误力（二维问题）。当然无论哪种处理都有一定的误差，且随计算范围的减小而增大，差，且随计算范围的减小而增大，靠近边界处误差靠近边界

12、处误差最大最大，这叫做，这叫做“边界效应边界效应”。在动力分析中影响更。在动力分析中影响更为显著，需妥善处理。为显著，需妥善处理。土木工程学院隧道工程系土木工程学院隧道工程系施成华 在生成初始应力场后边界条件应该重新设置，取消边在生成初始应力场后边界条件应该重新设置，取消边界力，改为位移边界条件。界力，改为位移边界条件。 对于采用无界元的情况，则在无限远处边界条件自然对于采用无界元的情况，则在无限远处边界条件自然满足，不用再设边界条件。但是要注意产生刚体位移的可满足，不用再设边界条件。但是要注意产生刚体位移的可能。能。土木工程学院隧道工程系土木工程学院隧道工程系施成华 初始应力场主要由初始应力

13、场主要由岩体自重和地质构造力岩体自重和地质构造力产生。但如何产生。但如何正确地确定这种应力场至今未得到妥善解决，因为构造应力正确地确定这种应力场至今未得到妥善解决，因为构造应力常常分布极不均匀，而费用昂贵的现场地应力测量只能给出常常分布极不均匀，而费用昂贵的现场地应力测量只能给出计算范围中少数几个点的地应力值。计算范围中少数几个点的地应力值。 一种常用的方法是根据自重应力场及构造应力场的特点，一种常用的方法是根据自重应力场及构造应力场的特点，确定较符合计算区域地质特点的力边界条件，并利用部分量确定较符合计算区域地质特点的力边界条件，并利用部分量测数据进行调整和修正。这与简单地硬凑已知点地应力是

14、不测数据进行调整和修正。这与简单地硬凑已知点地应力是不同的。同的。 土木工程学院隧道工程系土木工程学院隧道工程系施成华 自重应力场的特点是垂直及水平方向为主应力方向，主自重应力场的特点是垂直及水平方向为主应力方向，主应力均为压应力，其大小仅与深度有关而与水平位置及时间应力均为压应力，其大小仅与深度有关而与水平位置及时间无关。一般认为由自重产生的垂直及水平地应力为：无关。一般认为由自重产生的垂直及水平地应力为： v v= = H H, , h h= = v vK K= = H H( ( /(1-/(1- ) ) ) 而构造应力场主要与岩性分布和构造形式有关。而构造应力场主要与岩性分布和构造形式有

15、关。坚硬完坚硬完整的岩体中往往构造残余应力较高整的岩体中往往构造残余应力较高，而破碎松软岩体中就较，而破碎松软岩体中就较低，沿河谷附近的岩体由于卸荷作用会使地应力方向和大小低，沿河谷附近的岩体由于卸荷作用会使地应力方向和大小发生改变。一般也不考虑构造应力场与时间的关系。隧道工发生改变。一般也不考虑构造应力场与时间的关系。隧道工程地下洞室往往埋深不会太大，不会受到地热影响，不必考程地下洞室往往埋深不会太大，不会受到地热影响，不必考虑温度应力场。虑温度应力场。土木工程学院隧道工程系土木工程学院隧道工程系施成华 式中式中x,yx,y为坐标为坐标; ;为岩体材料参数为岩体材料参数; ;F Fx x，F

16、 Fy y为构造作用力；为构造作用力； 为边界位移为边界位移, , T T为温度。它们可由各种计算模型给出。先将为温度。它们可由各种计算模型给出。先将计算区域离散化，给定计算区域离散化，给定 ， F Fx x，F Fy y 等的初值，用有限元法计等的初值，用有限元法计算域内应力，求出相应于特定算域内应力，求出相应于特定 , , F Fx, x, F Fy y等因素的基本初始应力等因素的基本初始应力 , , Fx,Fx, FyFy ，将其与相应回归系数相乘并迭加，得到初始，将其与相应回归系数相乘并迭加，得到初始应力场的回归方程：应力场的回归方程：) 1 (.),(TFFEyxfyx) 2 (32

17、1kyFxbbb 土木工程学院隧道工程系土木工程学院隧道工程系施成华式中式中1 1，2 2为回归系数；为回归系数；k k为观测误差。当有个观测为观测误差。当有个观测值时，应有：值时，应有：1 1）观测误差）观测误差kk的数学期望值全为零，即的数学期望值全为零，即) 3(.3 , 2 , 1, 0)(nkEk) 4()()(0,cov2hkhkhk 各测点的现场量测值各测点的现场量测值k k为个独立观测值，为个独立观测值，为个观测值的为个观测值的总体。由各基本因素总体。由各基本因素, , F Fx x，F Fy y 所得的基本初始应力所得的基本初始应力 k,k, Fxk,Fxk, FkyFky,

18、为方程（为方程（2 2）的自变量。根据各实测点提供的组实测值，）的自变量。根据各实测点提供的组实测值，以及由数值方法计算的以及由数值方法计算的“数字观测值数字观测值”给出的各回归系数估计值给出的各回归系数估计值b b1 1，b b2 2，可以算出误差估计值，可以算出误差估计值e ek k和残差平方和和残差平方和： 土木工程学院隧道工程系土木工程学院隧道工程系施成华 形成由个方程组成的法方程组，解出回归系数形成由个方程组成的法方程组，解出回归系数b bi i，再用再用b bi i乘以基本初始应力即得到初始应力场。为使方程组有唯一解，乘以基本初始应力即得到初始应力场。为使方程组有唯一解，地应力量测

19、点数至少应等于地应力量测点数至少应等于m m。回归分析质量可由相关分析，回归分析质量可由相关分析，方差估计，显著性检验等予以检查验证。方差估计，显著性检验等予以检查验证。) 5(.2112112mkFykFxkkkmkkbbbeQ)6().3,2, 1(0mibQi土木工程学院隧道工程系土木工程学院隧道工程系施成华 还有人采用还有人采用模拟地形演化过程来生成地应力场模拟地形演化过程来生成地应力场。即第一。即第一步先模拟一块均匀的地壳，范围要求取得与工程范围相比足步先模拟一块均匀的地壳，范围要求取得与工程范围相比足够大。地层分布与实际工程区域相同。除了自重之外，边界够大。地层分布与实际工程区域相

20、同。除了自重之外，边界上还作用有根据区域地应力量值和方向所确定的分布力。进上还作用有根据区域地应力量值和方向所确定的分布力。进行平衡计算，迭代收敛到两次计算误差小于给定误差后，即行平衡计算，迭代收敛到两次计算误差小于给定误差后，即认为形成了远古应力场。然后采用分期开挖的方式，模拟地认为形成了远古应力场。然后采用分期开挖的方式，模拟地形、地貌形成的过程，一直到模拟的地形与实际地形相一致。形、地貌形成的过程，一直到模拟的地形与实际地形相一致。这时计算域中的应力场分布就可以认为是当前的应力场。这时计算域中的应力场分布就可以认为是当前的应力场。土木工程学院隧道工程系土木工程学院隧道工程系施成华各部分材

21、料的容重采用实际值各部分材料的容重采用实际值将各种材料均变为线弹性材料，不考虑非线性和弹将各种材料均变为线弹性材料，不考虑非线性和弹塑性塑性每计算一步，将各结点位移重新赋为零值每计算一步，将各结点位移重新赋为零值生成满意的初始应力场后，再将材料参数置换成实生成满意的初始应力场后，再将材料参数置换成实际参数，然后进行开挖计算。际参数，然后进行开挖计算。土木工程学院隧道工程系土木工程学院隧道工程系施成华 隧道施工过程主要包括洞室的开挖、喷射混凝土和锚隧道施工过程主要包括洞室的开挖、喷射混凝土和锚杆锚索的设置、二次衬砌混凝土结构的浇筑等等。这些施工杆锚索的设置、二次衬砌混凝土结构的浇筑等等。这些施工

22、过程都相当于在原始地应力场中增加新的荷载或改变地下结过程都相当于在原始地应力场中增加新的荷载或改变地下结构的材料而产生二次、三次构的材料而产生二次、三次应力场。这是隧道工程数值应力场。这是隧道工程数值分析的一个重要特点。分析的一个重要特点。 土木工程学院隧道工程系土木工程学院隧道工程系施成华 岩体在开挖隧道之前是处于一定的初始应力状态，开挖使隧岩体在开挖隧道之前是处于一定的初始应力状态，开挖使隧道周边上各点的应力道周边上各点的应力“解除解除”，从而引起围岩应力场的变化。，从而引起围岩应力场的变化。如果在开挖的同时设置了支护结构并与围岩密贴，则支护结构如果在开挖的同时设置了支护结构并与围岩密贴，

23、则支护结构将约束围岩因应力场变化而产生的位移，支护结构中也将产生将约束围岩因应力场变化而产生的位移，支护结构中也将产生应力和位移。所以，在进行有限元分析时，必须设法模拟这个应力和位移。所以，在进行有限元分析时，必须设法模拟这个开挖卸荷的效果。通用的方法就是在隧道周边的点上加开挖卸荷的效果。通用的方法就是在隧道周边的点上加“等效等效释放荷载释放荷载”。 这些这些“等效释放荷载等效释放荷载”是由于隧道周边各点的应力是由于隧道周边各点的应力“解除解除”而形成的，因此，可以根据沿预定周边上的初始应力来确定。而形成的，因此，可以根据沿预定周边上的初始应力来确定。 土木工程学院隧道工程系土木工程学院隧道工

24、程系施成华 设沿预计开挖边界上各点的初始应力为已知，在离散化的情设沿预计开挖边界上各点的初始应力为已知，在离散化的情况下，可假定沿开挖面上两相邻节点之间的初始应力呈线性变况下，可假定沿开挖面上两相邻节点之间的初始应力呈线性变化，对于任一开挖边界点，开挖引起的化，对于任一开挖边界点，开挖引起的“等效释放荷载等效释放荷载”（等（等效节点力）即可按简支梁分配的原则进行置换。效节点力）即可按简支梁分配的原则进行置换。 在实际模拟计算中，具体的在实际模拟计算中，具体的“等效释放荷载等效释放荷载”一般是由计算一般是由计算软件自动完成的。软件自动完成的。 考虑到隧道施工过程中采用多次支护，考虑到隧道施工过程

25、中采用多次支护，“等效释放荷载等效释放荷载”由由围岩、初支、二衬共同承担，计算中则应根据计算经验人为设围岩、初支、二衬共同承担，计算中则应根据计算经验人为设定它们各自承担释放荷载的比例，即定它们各自承担释放荷载的比例，即“荷载释放系数荷载释放系数”。 土木工程学院隧道工程系土木工程学院隧道工程系施成华隧道工程的开挖和支护过程都是分期进行，相互交替的，隧道工程的开挖和支护过程都是分期进行，相互交替的，因此数值分析过程也要模拟这种施工过程。首先，因此数值分析过程也要模拟这种施工过程。首先，在划分洞室在划分洞室内部单元时就必须考虑整个施工程序内部单元时就必须考虑整个施工程序，所有开挖和不同衬砌部，所

26、有开挖和不同衬砌部分的边线都必须是单元的边线，而不能在单元内部。衬砌施作分的边线都必须是单元的边线，而不能在单元内部。衬砌施作过程的模拟比较简单，即在开挖之后某一规定的分期内，将衬过程的模拟比较简单，即在开挖之后某一规定的分期内，将衬砌部分对应的砌部分对应的“空单元空单元”重新赋予衬砌材料的参数后再进行计重新赋予衬砌材料的参数后再进行计算。算。适当改变开挖和衬砌施作方案，比较围岩应力和变形情况，适当改变开挖和衬砌施作方案，比较围岩应力和变形情况，对确定最优施工程序是非常有效的。对确定最优施工程序是非常有效的。土木工程学院隧道工程系土木工程学院隧道工程系施成华隧道开挖施工过程的具体模拟方法如下：

27、隧道开挖施工过程的具体模拟方法如下：（1 1）按照施工要求划分好开挖顺序。）按照施工要求划分好开挖顺序。（2 2）按照隧道埋深的地质构造特点，进行开挖前的应力分析，）按照隧道埋深的地质构造特点，进行开挖前的应力分析，求出围岩中的初始地应力场和位移场，开挖前的应力状态可作求出围岩中的初始地应力场和位移场，开挖前的应力状态可作为原始数据直接输入。为原始数据直接输入。（3 3）根据每次开挖的尺寸，去掉被开挖的单元，根据去掉单）根据每次开挖的尺寸，去掉被开挖的单元，根据去掉单元现时的应力值，求出被开挖出的自由表面各节点处由这些单元现时的应力值，求出被开挖出的自由表面各节点处由这些单元作用的节点力。将与

28、这些节点力大小相等、方向相反的力作元作用的节点力。将与这些节点力大小相等、方向相反的力作用于自由表面相同的节点上，这些力就是等效开挖释放荷载。用于自由表面相同的节点上，这些力就是等效开挖释放荷载。（4 4）在等效开挖释放荷载作用下进行分析，求出该开挖步骤）在等效开挖释放荷载作用下进行分析，求出该开挖步骤后围岩中的位移、应变、应力，并叠加于以前的状态上，重复后围岩中的位移、应变、应力，并叠加于以前的状态上，重复以上步骤，直至最后一个开挖步骤结束。以上步骤，直至最后一个开挖步骤结束。土木工程学院隧道工程系土木工程学院隧道工程系施成华 对于施工中采用锚喷支护的模拟有如下几种考虑。对于施工中采用锚喷支

29、护的模拟有如下几种考虑。（1) 1) 端部锚固的锚杆的轴向刚度端部锚固的锚杆的轴向刚度 121122241TTuuuuQQEdlrSSKebebbbilcb土木工程学院隧道工程系土木工程学院隧道工程系施成华2112112dddEGdKbga4312142dddIEEIEKbgbs土木工程学院隧道工程系土木工程学院隧道工程系施成华（3 3）L.J. Lorig L.J. Lorig 锚杆最大抗剪力模型锚杆最大抗剪力模型 0, 021maxcoscos18sbysbusigndF土木工程学院隧道工程系土木工程学院隧道工程系施成华 这样，可以用在不连续面附近建立两组弹簧来模拟这这样，可以用在不连续面

30、附近建立两组弹簧来模拟这一活动段，一组沿锚杆轴线方向，一组垂至于轴线。当一活动段，一组沿锚杆轴线方向，一组垂至于轴线。当发生剪切位移时，轴线方向的弹簧发生偏转，但仍保持发生剪切位移时，轴线方向的弹簧发生偏转，但仍保持与活动段轴线方向一致，切向弹簧的方向不出现改变。与活动段轴线方向一致，切向弹簧的方向不出现改变。土木工程学院隧道工程系土木工程学院隧道工程系施成华2122nsafuuur2maxmaxmax)(sbssbssbsssssFFFFFFfFfuKF土木工程学院隧道工程系土木工程学院隧道工程系施成华 锚杆、锚索的设置也有人不考虑对整体刚度的影响，锚杆、锚索的设置也有人不考虑对整体刚度的影

31、响，而作为一种附加荷载施加于相应位置的节点上，尤其是端而作为一种附加荷载施加于相应位置的节点上，尤其是端部锚固的锚杆和锚索是这样处理的。部锚固的锚杆和锚索是这样处理的。通过计算两锚固点之通过计算两锚固点之间的相对位移，可得到锚杆（索）内的拉应力间的相对位移，可得到锚杆（索）内的拉应力，乘以锚杆，乘以锚杆断面面积，则可得到锚固节点力，反加到节点上再进行下断面面积，则可得到锚固节点力，反加到节点上再进行下一期的计算。计算中采用锚杆材料的本构模型进行判断，一期的计算。计算中采用锚杆材料的本构模型进行判断，如锚杆应力。超过屈服强度人则进入塑性阶段，超过强度如锚杆应力。超过屈服强度人则进入塑性阶段，超过

32、强度极限则发生拉断，锚杆拉力释放作用到锚固节点上去。极限则发生拉断，锚杆拉力释放作用到锚固节点上去。土木工程学院隧道工程系土木工程学院隧道工程系施成华 对于对于全长锚固的锚杆全长锚固的锚杆，也可将其分割为若干段，每，也可将其分割为若干段，每段均简化为端部锚固的锚杆来处理。沿锚杆分布的剪应段均简化为端部锚固的锚杆来处理。沿锚杆分布的剪应力，按其分布现律化为等效节点力施加于锚杆通过的各力，按其分布现律化为等效节点力施加于锚杆通过的各个节点，再作用到岩体上个节点，再作用到岩体上土木工程学院隧道工程系土木工程学院隧道工程系施成华 对数值计算结果的分析是非常重要的，而这一点对对数值计算结果的分析是非常重

33、要的，而这一点对于经验不多的人来说又非常困难。有的人完成计算后于经验不多的人来说又非常困难。有的人完成计算后得出了一大堆数据和图形，但不知道如何下手分析它，得出了一大堆数据和图形，但不知道如何下手分析它，看不出这些数据说明什么问题，又能得出什么结论。看不出这些数据说明什么问题，又能得出什么结论。 根据多年来进行隧道工程计算的经验和教训，感到根据多年来进行隧道工程计算的经验和教训，感到可以从如下几个方面对计算结果进行分析和判断：可以从如下几个方面对计算结果进行分析和判断：土木工程学院隧道工程系土木工程学院隧道工程系施成华 位移的时间空间分布形式的比较，包括方向和大小位移的时间空间分布形式的比较，

34、包括方向和大小洞周位移是否连续，有无突变，位移值是洞周位移是否连续，有无突变，位移值是否均匀，离散程度如何。否均匀，离散程度如何。位移随深度的变化情况，位移梯度的大小，位移随深度的变化情况，位移梯度的大小，位移场分布是否合理位移场分布是否合理不同洞室相应部位的位移的比较不同洞室相应部位的位移的比较位移值随开挖或衬砌施作等施工过程的变位移值随开挖或衬砌施作等施工过程的变化化土木工程学院隧道工程系土木工程学院隧道工程系施成华 分析应力等值线、主应力矢量图分析应力等值线、主应力矢量图应力分布与地形地貌的关系是否符合一般应力分布与地形地貌的关系是否符合一般规律？应力随深度的变化情况，应力梯度的规律？应

35、力随深度的变化情况，应力梯度的大小，应力场分布是否合理大小，应力场分布是否合理洞周应力流的分布如何。洞周应力流的分布如何。应力集中出现的部位的比较，应力集中系应力集中出现的部位的比较，应力集中系数的大小数的大小拉应力出现的部位，量值拉应力出现的部位，量值应力值随开挖或衬砌施作等施工过程的变应力值随开挖或衬砌施作等施工过程的变化化土木工程学院隧道工程系土木工程学院隧道工程系施成华屈服区、塑性区分析屈服区、塑性区分析 屈服区分布型式屈服区分布型式屈服区的深度屈服区的深度屈服破坏的类型，与应力的关系屈服破坏的类型，与应力的关系屈服范围的面积或体积量值屈服范围的面积或体积量值屈服区随开挖或衬砌施作等施

36、工过程的发屈服区随开挖或衬砌施作等施工过程的发展展从能量角度分析屈服区从能量角度分析屈服区土木工程学院隧道工程系土木工程学院隧道工程系施成华锚杆、锚索应力分析锚杆、锚索应力分析 锚杆应力沿洞周分布型式锚杆应力沿洞周分布型式锚杆应力占屈服极限和强度极限的比例。锚杆应力占屈服极限和强度极限的比例。有无受压的锚杆有无受压的锚杆最大拉应力锚杆出现的位置最大拉应力锚杆出现的位置锚杆应力随开挖或建造等施工过程的发展锚杆应力随开挖或建造等施工过程的发展增加锚杆后对位移的影响增加锚杆后对位移的影响增加锚杆后对应力分布的影响增加锚杆后对应力分布的影响增加锚杆后对屈服区大小的影响增加锚杆后对屈服区大小的影响土木工

37、程学院隧道工程系土木工程学院隧道工程系施成华 根据数值分析计算结果，如何合理地判断围岩的稳根据数值分析计算结果，如何合理地判断围岩的稳定性也是当前尚未解决的一个问题。在数值分析方法定性也是当前尚未解决的一个问题。在数值分析方法中，除非将支护结构离散为梁单元，否则都只能求得中，除非将支护结构离散为梁单元，否则都只能求得支护结构中的应力，因而不能直接采用规范中的公式支护结构中的应力，因而不能直接采用规范中的公式校核支护结构的强度。目前采用的判断围岩与支护结校核支护结构的强度。目前采用的判断围岩与支护结构稳定性的方法主要有如下几种：构稳定性的方法主要有如下几种：土木工程学院隧道工程系土木工程学院隧道

38、工程系施成华（1 1）超载系数法超载系数法。将外荷载乘以系数值，并逐渐增大值。将外荷载乘以系数值，并逐渐增大值进行反复计算，直到计算不能收敛为止，即认为围岩进行反复计算，直到计算不能收敛为止，即认为围岩失稳，值为安全系数。失稳，值为安全系数。（2 2）材料安全储备法（强度折减法）材料安全储备法（强度折减法）。将材料的主要强。将材料的主要强度特征值，如、乘以值，逐渐降低值并反复计算到围度特征值，如、乘以值，逐渐降低值并反复计算到围岩失稳（即计算不收敛）为止，就是安全度。岩失稳（即计算不收敛）为止，就是安全度。（3 3）经验类比法经验类比法。将计算所得洞壁变形值或塑性区范围。将计算所得洞壁变形值或塑性区范围与按经验所得的围岩