第4章山岭隧道衬砌计算ppt课件

1、第四章山岭隧道衬砌计算第四章山岭隧道衬砌计算 从目前地下结构理论研究成果看，山岭隧道设从目前地下结构理论研究成果看，山岭隧道设计计算方法大致分为计计算方法大致分为: 结构力学方法结构力学方法 岩体力学方法岩体力学方法 监控设计法及经验设计法监控设计法及经验设计法第一节第一节 结构力学方法结构力学方法一、基本原理一、基本原理 结构力学方法是将支护和围岩分开考虑的，支护结构力学方法是将支护和围岩分开考虑的，支护结构是承载主体。当作用在支护结构上的荷载确定后，结构是承载主体。当作用在支护结构上的荷载确定后，可应用普通结构力学的方法求解超静定结构的内力和可应用普通结构力学的方法求解超静定结构的内力和位

2、移。位移。1主动荷载模式主动荷载模式该模式只考虑支护结构在该模式只考虑支护结构在围岩主动荷载作用下的受围岩主动荷载作用下的受力情况，它主要应用于采力情况，它主要应用于采用浅埋暗挖法或明挖法施用浅埋暗挖法或明挖法施工的城市地铁工程及明洞工的城市地铁工程及明洞工程。工程。2主动荷载加被动荷载弹性抗力形式主动荷载加被动荷载弹性抗力形式 该模式认为围岩不仅该模式认为围岩不仅对支护结构施加主动对支护结构施加主动荷载，而且围岩由于荷载，而且围岩由于受支护结构的挤压作受支护结构的挤压作用，对支护结构还施用，对支护结构还施加约束反力弹性抗加约束反力弹性抗力）。力）。 3实际荷载模式实际荷载模式它是采用量测仪器

3、实地量测作它是采用量测仪器实地量测作用在衬砌上的荷载值，这样能用在衬砌上的荷载值，这样能综合反映围岩与支护结构的相综合反映围岩与支护结构的相互作用。当支护结构与围岩牢互作用。当支护结构与围岩牢固接触时，不仅能量测到径向固接触时，不仅能量测到径向荷载而且还能量测到切向荷载。荷载而且还能量测到切向荷载。及支护结构背后回填的质量。及支护结构背后回填的质量。二、隧道围岩的弹性抗力二、隧道围岩的弹性抗力 隧道衬砌变形隧道衬砌变形弹性抗力：弹性抗力：就是由于支护结构发生向围岩方向的变形而引起就是由于支护结构发生向围岩方向的变形而引起的围岩对支护结构的约束力的围岩对支护结构的约束力弹性抗力的大小，目前常用文

4、克尔弹性抗力的大小，目前常用文克尔Winkler假假定为基础的局部变形弹簧模型理论来确定。该理论认定为基础的局部变形弹簧模型理论来确定。该理论认为围岩的弹性抗力与围岩在该点的变形成正比，用公为围岩的弹性抗力与围岩在该点的变形成正比，用公式表示为：式表示为：iiK文氏假定，相当于把文氏假定，相当于把围岩简化为一系列彼此独围岩简化为一系列彼此独立的弹簧，某一弹簧受到立的弹簧，某一弹簧受到压缩时所产生的反作用力压缩时所产生的反作用力只与该弹簧有关，而与其只与该弹簧有关，而与其他弹簧无关。他弹簧无关。 弹簧模型弹簧模型 三、隧道衬砌荷载和荷载组合三、隧道衬砌荷载和荷载组合（1主要荷载主要荷载长期及经常

5、作用的荷载为主要荷载，如：长期及经常作用的荷载为主要荷载，如：围岩压力围岩压力围岩对支护的弹性抗力围岩对支护的弹性抗力支护结构的自重支护结构的自重地下水压力及列车或汽车活载等地下水压力及列车或汽车活载等（2附加荷载附加荷载即偶然的、非经常作用的荷载，如：即偶然的、非经常作用的荷载，如：温差应力温差应力灌浆压力冻胀力灌浆压力冻胀力地震力地震力其中主要的是地震力，其大小可按有关的抗震其中主要的是地震力，其大小可按有关的抗震设计规范规定执行。设计规范规定执行。公路隧道标准将隧道所承受的荷载分为：公路隧道标准将隧道所承受的荷载分为：永久荷载永久荷载可变荷载可变荷载偶然荷载偶然荷载四、隧道衬砌结构按结构

6、力学原理计算的几种方法四、隧道衬砌结构按结构力学原理计算的几种方法 隧道衬砌结构计算的主要内容有：隧道衬砌结构计算的主要内容有：按工程类比方法初步拟定衬砌断面的几何尺寸，按工程类比方法初步拟定衬砌断面的几何尺寸，确定作用在衬砌结构上的荷载，进行力学计算，求出确定作用在衬砌结构上的荷载，进行力学计算，求出衬砌截面的弯距衬砌截面的弯距M和轴向力和轴向力N等内力；验算衬砌截面的等内力；验算衬砌截面的承载安全系数承载安全系数K值。值。 （一假定抗力图形法（一假定抗力图形法1基本原理基本原理根据经验，某些形式衬砌的弹性抗力有一定的分根据经验，某些形式衬砌的弹性抗力有一定的分布规律。通过假定抗力分布区范围

7、及抗力规律，来求布规律。通过假定抗力分布区范围及抗力规律，来求解衬砌的内力大小。解衬砌的内力大小。采用采用“假定抗力图形法这种方法，求得某一特假定抗力图形法这种方法，求得某一特定点的弹性抗力，就可求出其他各点的弹性抗力值。定点的弹性抗力，就可求出其他各点的弹性抗力值。 抗力区的上零点在拱顶两侧约抗力区的上零点在拱顶两侧约45附近的附近的b点点45），下零点在墙脚），下零点在墙脚a点，最大抗力发生在点，最大抗力发生在h点，点，h点的位置约在抗力区点的位置约在抗力区23高度处，即弧高度处，即弧ah（23ah弧，或者在最大跨度附近；弧，或者在最大跨度附近；bh段抗力图形按二次抛物线分布。段抗力图形按

8、二次抛物线分布。2抗力分布假设抗力分布假设假假定定坑坑力力图图形形法法在在bh弧段上有：弧段上有：hhbibi2222coscoscoscos在在ah弧段上有：弧段上有：hiiiyy2213计算方法计算方法h点的衬砌变形：点的衬砌变形： （1求在主动荷载作用下衬砌截面求在主动荷载作用下衬砌截面i的内力：的内力： （2求在单位被动荷载即单位抗力求在单位被动荷载即单位抗力h1作用下作用下所产生的所产生的i截面的内力：截面的内力：（3求最大抗力求最大抗力（4求衬砌截面的总内力：求衬砌截面的总内力：（二弹性地基梁法（二弹性地基梁法1基本原理基本原理这种方法是将衬砌结构看成是置于弹性地基上的这种方法是将

9、衬砌结构看成是置于弹性地基上的曲梁或直梁，弹性地基上抗力按文克尔假定的局部变曲梁或直梁，弹性地基上抗力按文克尔假定的局部变形理论求解。当曲墙的曲率是常数或直墙时，可采用形理论求解。当曲墙的曲率是常数或直墙时，可采用初参数法求解结构内力，一般直墙式衬砌的直边墙采初参数法求解结构内力，一般直墙式衬砌的直边墙采用此方法求解。用此方法求解。2计算方法计算方法 直墙式衬砌将其拱、墙分开计算。拱圈为一个弹直墙式衬砌将其拱、墙分开计算。拱圈为一个弹性固定在边墙顶上的无铰平拱，边墙为一个置于弹性性固定在边墙顶上的无铰平拱，边墙为一个置于弹性地基上的直梁。边墙在拱脚和墙顶连接处应满足力的地基上的直梁。边墙在拱脚

10、和墙顶连接处应满足力的平衡条件和变形协调条件。拱脚处除有转角和竖直位平衡条件和变形协调条件。拱脚处除有转角和竖直位移外，尚有水平位移。在结构和荷载均对称时，拱脚移外，尚有水平位移。在结构和荷载均对称时，拱脚的竖直位移对拱圈内力无影响，可以不计。的竖直位移对拱圈内力无影响，可以不计。 直墙式衬砌计算图式直墙式衬砌计算图式 ddbibi2222coscoscoscosdhdsin（三弹性链杆法（三弹性链杆法1基本原理基本原理 采用符合采用符合“局部变形原理的弹簧地基来模拟围岩；再将衬局部变形原理的弹簧地基来模拟围岩；再将衬砌与围岩二者所组成的衬砌结构体系离散为有限个衬砌单元和砌与围岩二者所组成的衬

11、砌结构体系离散为有限个衬砌单元和弹簧单元所组成的组合体；其次确定由主动荷载如围岩压力、弹簧单元所组成的组合体；其次确定由主动荷载如围岩压力、衬砌自重等直接作用在计算模式上所引起的变形，变形的轮衬砌自重等直接作用在计算模式上所引起的变形，变形的轮廓就是衬砌与围岩相互作用区域。廓就是衬砌与围岩相互作用区域。 （1衬砌本身理想化衬砌本身理想化 对于隧道衬砌结构来说，其理想化包括下对于隧道衬砌结构来说，其理想化包括下面四部分内容面四部分内容:围岩围岩弹性弹性链杆链杆的设的设置置 （2围岩的理想化围岩的理想化左图示左图示等效节点荷载计算等效节点荷载计算 （3荷载理想化荷载理想化等效节点荷载等效节点荷载

12、（4单元和节点编号，边界条件单元和节点编号，边界条件 2隧道衬砌结构杆系有限元分析隧道衬砌结构杆系有限元分析 （1计算图式计算图式 （2单元刚度矩阵。单元刚度矩阵。计算衬砌结构内力共用计算衬砌结构内力共用3种单元：偏心受压的衬砌单元；支承链杆单元；种单元：偏心受压的衬砌单元；支承链杆单元；墙底地层弹性固定的弹性支座单元。墙底地层弹性固定的弹性支座单元。衬砌单元与坐标系衬砌单元与坐标系衬砌单元的刚度矩阵。衬砌单衬砌单元的刚度矩阵。衬砌单元的截面本来是变化的，但为元的截面本来是变化的，但为了便于计算，把它看作是截面了便于计算，把它看作是截面尺寸不变的在单元端点作用有尺寸不变的在单元端点作用有轴力瓦

13、弯矩元和剪力轴力瓦弯矩元和剪力Q的等直的等直梁单元梁单元简写为：简写为：简写为：简写为：同理，可求得单元节点位移在两个坐标系中的关系式：同理，可求得单元节点位移在两个坐标系中的关系式：式子代入得到：式子代入得到：则有：则有：上式用分块矩阵表达：上式用分块矩阵表达：展开得：展开得：（）（）（）（）（）（）（）（） 支承链杆单元刚度矩阵。支承链杆单元支承链杆单元刚度矩阵。支承链杆单元i的弹性抗力三；与其端点压缩位的弹性抗力三；与其端点压缩位移民的关系按文克尔假定可写为移民的关系按文克尔假定可写为 :也即：也即：衬砌单元与支承链杆单元衬砌单元与支承链杆单元 （）（）（）（）支承链杆单元的刚度矩阵为：

14、支承链杆单元的刚度矩阵为：上图所示的衬砌单元与支承链杆单元上图所示的衬砌单元与支承链杆单元 简写为：简写为：墙底弹性支座单元墙底弹性支座单元 （）（）（）（）（3结构刚度方程的形成结构刚度方程的形成 现将式现将式421改写成如下矩阵形式：改写成如下矩阵形式：（）（）（）（）（）（）根据静力平衡条件：在同一节点上作用的各单元节点力的总和应与该节根据静力平衡条件：在同一节点上作用的各单元节点力的总和应与该节点的节点荷载相平衡，因此有：点的节点荷载相平衡，因此有：（）（）（）（）（4衬砌内力计算衬砌内力计算 在求得最终的结构节点位移值后，便可求解衬砌的内力。在求得最终的结构节点位移值后，便可求解衬砌

15、的内力。根据变形协调条件，结构节点位移与汇交于此节点的单元根据变形协调条件，结构节点位移与汇交于此节点的单元节点位移相等，即：节点位移相等，即：（）（）（）（）（5隧道衬砌结构杆有限元分析基本步骤如下：隧道衬砌结构杆有限元分析基本步骤如下：衬砌结构离散化。衬砌结构离散化。将所有荷载换算成等效节点荷载；将所有荷载换算成等效节点荷载；计算各衬砌单元的刚度矩阵按局部坐标系）；计算各衬砌单元的刚度矩阵按局部坐标系）；计算各支承链杆单元的弹性支座单元的刚度矩阵局部坐标系和结计算各支承链杆单元的弹性支座单元的刚度矩阵局部坐标系和结构坐标系相一致）；构坐标系相一致）；把局部坐标系建立的衬砌单元刚度矩阵转换成

16、结构坐标系的衬砌单把局部坐标系建立的衬砌单元刚度矩阵转换成结构坐标系的衬砌单元刚度矩阵；元刚度矩阵；形成结构刚度矩阵形成结构刚度矩阵K；按衬砌结构的边界条件位移边界条件修改结构刚度方程；按衬砌结构的边界条件位移边界条件修改结构刚度方程；求节点位移的第一次近似值，并检查和取消位移求节点位移的第一次近似值，并检查和取消位移ui0的支承链杆的支承链杆单元，修改刚度矩阵，重新计算；单元，修改刚度矩阵，重新计算；据节点水平位移的第二次近似值，反复据节点水平位移的第二次近似值，反复、步骤，直到所有支承步骤，直到所有支承链杆处的节点位移均符合链杆处的节点位移均符合ui0的条件为止；的条件为止；四求解最终的结

17、构节点位移值；四求解最终的结构节点位移值；（计算局部坐标系中个衬砌单元的节点位移值；（计算局部坐标系中个衬砌单元的节点位移值；（计算拱顶截面支座约束反力（计算拱顶截面支座约束反力 （432） （433） 五、隧道衬砌截面强度检算五、隧道衬砌截面强度检算 检算所得的安全系数是否满足相应规范所要求的检算所得的安全系数是否满足相应规范所要求的数值数值 ： 拱形隧道和明洞衬砌属偏心受压构件，其截面极限承载力Nz的计算，按轴力的偏心距eMN的大小有两种情况： （1抗压强度控制当抗压强度控制当 e0.20 d）：）： （2抗拉强度控制当抗拉强度控制当 e0.20 d）：）：第二节第二节 岩体力学方法岩体力

18、学方法整体复合模型整体复合模型.围岩本身由于支护结构提供了一定的支护抗力，而围岩本身由于支护结构提供了一定的支护抗力，而引起它的应力调整，达到新的稳定；引起它的应力调整，达到新的稳定；.由于支护结构阻止围岩变形，也必然受到围岩给予由于支护结构阻止围岩变形，也必然受到围岩给予的反作用力而发生变形。的反作用力而发生变形。目前对于这种模型求解方法有：目前对于这种模型求解方法有：解析法、数值法、特征曲线法、剪切滑移破坏法解析法、数值法、特征曲线法、剪切滑移破坏法一、解析法一、解析法1基本假设基本假设 （深埋圆形隧道，无限长；（深埋圆形隧道，无限长；（原岩应力各向等压；（原岩应力各向等压；（围岩为理想弹

19、塑性体，符合库仑强度准则。（围岩为理想弹塑性体，符合库仑强度准则。2基本方程基本方程 2rBAr弹性区：积分常数待定的弹性应力解弹性区：积分常数待定的弹性应力解 塑性区：平衡方程塑性区：平衡方程 0rddrr强度准则方程强度准则方程库仑准则库仑准则sin1cos2sin1sin1Cr塑性区：平衡方程塑性区：平衡方程 0rddrr3边界条件边界条件在弹性区内，外边界处有在弹性区内，外边界处有 r，rP0，在弹性区的内边界交界面在弹性区的内边界交界面rRP塑性区半径）。塑性区半径）。2PeerRBA在塑性区的外边界交界面）在塑性区的外边界交界面）ePerPrPRr在内边界周边处：在内边界周边处：r

20、 =0不支护），不支护），r =P支护反力）支护反力） 1cotsin1sin20RrCpr1sin1sin1cotsin1sin20RrCp4问题的求解问题的求解式式438）、式）、式439联解并利用式联解并利用式442有有（443） （44） 式式443代人式代人式439），整理得：），整理得：由式由式437及弹性边界条件有及弹性边界条件有20rBPeer（445） 由式由式443）、式）、式445）、塑性内边界条件式）、塑性内边界条件式441）（44） cossin11022220CPrRrRPPPeer（44） 解的解的B，再代人式，再代人式437），整理得弹性区应力：），整理得弹性区

21、应力：由式由式441）、式）、式444）、式）、式447又有又有1sin1sin1cotsin1sin20RRCPcossin1200CPP（44） 得塑性区半径得塑性区半径sin1sin200cotsin1)cot(CPRRP5结果的列出结果的列出 弹性区应力：弹性区应力： 00cosPCPrsin1sin20cotsin1)cot(CP20rR（44） 塑性区应力塑性区应力1cotsin1sin20RrCr（4） 1sin1sin1cotsin1sin20RrC（4） 塑性区半径塑性区半径sin1sin200cotsin1)cot(CPRRP（4） 若式若式438）、式）、式439联解后，

22、改用边界条件联解后，改用边界条件r =P决定积分常数，则可决定积分常数，则可得有支护反力得有支护反力P情况下的弹性区应力为：情况下的弹性区应力为：00cosPCPrsin1sin200cotsin1)cot(CPP20rR（）（）塑性区应力为塑性区应力为 1)cot(sin1sin200RrCPr1sin1sin1)cot(sin1sin200RrCP（）（）（）（）塑性区半径塑性区半径sin1sin2000cotsin1)cot(CPPRRP（）（）将式将式455改写为改写为 cot)sin1)(cot(sin1sin200CRRCPPP（）（）6讨论讨论Rp与与R0成正比，与成正比，与0成正变，与，成反变关系。成正变，与，成反变关系。塑性区应力与塑性区应力与0原始应力无关；原始应力无关；支护反力时，支护反力时， Rp最大。最大。二、数值方法二、数值方法数值计算方法，通常包括：数值计算方法，通常包括：有限差分法离散控制方程）有限差分法离散控制方程）有限元法离散解题全域）有限元法离散解题全域）边界元法只离散边界）边界元法只离散边界）离散元法块体运动）离散元法块体运动） （一有限元法（一有限元法 图图413 有限元分析系统的结构图有限元分析系统的结构图 （二地