第十章-盾构隧道衬砌计算方法

第十章盾构隧道衬砌计算方法10.1国内外的发展动态-常用模型盾构隧道的设计模型，多用荷载一结构模型。但由于其断面为圆形，地层结构法对均一地层中单孔圆形隧道也取得了精确的解析解，但其他情况仍须借助数值逼近。国际隧道协会(ITA)推荐使用日本(修正)惯用设计法、弹性地基梁法及有限元法。目前简化模型有均质圆环法、弹性铰法、多铰圆环模型和梁-弹簧模型。其中最常用的是均质圆环法，即日本的惯用设计法。日本提出并进一步发展了梁-弹簧模型。10.2荷载结构法的假设荷载结构法对于实际情况的假设如下:(I)对于盾构隧道管片简化、对实际情况的简化及基本假设。①假定地层和衬砌处于平面应变状态。②对于浅埋隧道，作用在衬砌上的主动土压力取为上覆地层全部厚度的土柱重。如有地面附加荷载，附加荷载需予考虑。对于深埋隧道，衬砌上的主动土压力用泰沙基土压力理论或普氏理论。③在地层和衬砌之间存在着约束，这种约束可能是径向的或切向的或两种都有。荷载结构法的假设④由于衬砌和地层的共同作用，衬砌的变形将引起地层的抗力。⑤地层和衬砌材料假定是弹性的。⑥衬砌是均质和连续的。由于接缝之间都设有一定数量的螺栓连接，且相邻环之间采用错缝拼装的形式来弥补接缝刚度的不足，因此衬砌可假定为均质的和连续的。(2)对于外荷载的简化。①竖向水压力和土压力按匀布荷载计算。②水平向水压力和土压力按均匀变化荷载考虑。③衬砌在竖向主动土压力作用下发生压扁变形，引起介质的被动抗力。10.3地层结构法认为衬砌与地层一起构成受力变形的整体，并可按连续介质力学原理来计算衬砌和周边地层的计算方法称为地层结构法。地层结构法由于实际情况的复杂性，只对小部分问题取得了解析解，而大部分问题只能依赖数值解。10.4均质圆环法计算管片衬砌内力自由变形圆环模型是将盾构衬砌结构视作在土体中自由变形的弹性均质圆环，不考虑管片接头刚度变化的影响，将地基抗力假定为三角形分布荷载进行计算。自由变形圆环内力求解采用弹性中心法，即根据荷载和结构均对称于竖直轴，取半结构进行分析，根据弹性中心处的相对角位移和相对水平位移为零的条件，列出力法方程，求出多余未知力，再根据多余未知力求出圆环任意截面上的内力。自由变形圆环法即属于荷载结构模型。它同样基于荷载结构法的假设以及简化。均质圆环法的外荷载包括竖向和水平土压力、水压力、自重压力和地基反力。用弹性中心法计算得到各个截面的内力。41141(12si,n:)(22)24(0.045)(1)4(1)4wlNiihihwhhKEIQKhPQgREIKqqqhRqRR1土水土水圆环水平直径处的水平位移地层弹性抗力系数抗力区截面和水平直径的夹角是按照自由变形圆环法求出的公式如下考虑管片接头影响的系数衬砌材料的刚度g----衬P砌的重度212112212((22)tan()2tan()4242)tan()2tan()4242wwgnwiiwwihhPPRPgQhhQQRRCCX1X211112222211100322220010220001111(cos)211cosPPhhhhhhhpPhpPXXRMRdRdEIEIEIRMRdRdEIEIEIRMdEIRMdEI竖向和水平土压力、水压力土压力——计算土压力时的土壤容重，地下水位以上的取土壤湿容重，地下水位以下的取水下容重。若隧道覆盖层厚度小于2D，则拱顶土压力按土柱高度确定(若有地面附加压力，则需加上附加压力)。若隧道覆盖层厚度不小于2D，则按泰沙基公式算折减后的竖向土压力水压力——一般情况下作用在衬砌上的水压力为静水压力。但为了简化计算，也可以将水压力取为拱顶以上和隧底以下其值分别与该处静水压力相等的匀布竖向压力以及由拱顶至隧底均匀变化的水平荷载，其值分别与在拱顶和隧底处的静水压力相等。地基反力衬砌和地层是相互作用的，对地层反力简化的不同将引起衬砌变形及内力的不同。均质圆环法假定作用于两侧的土反力是随着衬砌向地基内位移而产生的，竖向正负45°范围内为脱离区，而水平方向正负45°范围内，按三角形分布的荷载考虑(以水平直径上的点为顶点)，水平直径上点的土反力与衬砌向地基内的水平位移成正比。土壤介质的侧向弹性抗力系数的取值不同，衬砌结构内力的M、N值也随之相应的变化，对工程的经济意义很大，必须谨慎合理，保证结构具有一定的安全度。均质圆环法内力表10.5弹性地基梁法——结构模型弹性地基梁法将衬砌结构看成弹性地基上的圆环，当土体中管片衬砌产生变形时，衬砌周围的土体将阻止管片变形，即产生土体抗力。弹性地基梁法用弹性地基弹簧来模拟衬砌围岩的相互作用。弹性地基梁法从考虑土体抗力的不同可分为两种模型，一种是全周弹簧模型，另一种是局部弹簧模型。局部弹簧模型假定在拱顶90°范围内为脱离区，不产生土体抗力。全周弹簧模型，见图2。弹性地基梁法计算图式见图3。弹性地基梁法——计算方法用有限元法，把衬砌结构离散为有限个梁单元，对于梁单元，其单元刚度矩阵与结构力学中杆单元的单元刚度矩阵相同。取梁轴线为x轴，可写出该梁单元的单元刚度矩阵。将所有梁单元和弹簧单元在局部坐标系下的单元刚度矩阵变换为整体坐标系下的单元刚度矩阵,再把所有整体坐标系下的单元刚度矩阵组成总体刚度矩阵,然后将围岩压力转化为节点荷载,再利用边界条件求出梁单元的内力和位移。假定各节点位移以使地基弹簧受压为正,若求出某节点位移为负(向隧道内位移),即此处弹簧受拉,则将此处的地基弹簧去掉,重新进行计算,在去掉位移为负的节点处的地基弹簧后，经计算若某个已被去掉地基弹簧的节点处位移又为正,则需将此处的地基弹簧加上再重新计算,直到所有的地基弹簧都受压为止。10.6弹性铰法弹性铰法计算图式见图4(所受外荷载同弹性地基梁法,未示出)。弹性铰法考虑了管片接头对结构内力的影响,由于盾构隧道管片衬砌是通过多块圆弧管片用螺栓拼装起来的,各圆弧管片的接头处存在一个能承担一部分弯矩的弹性铰,其承担弯矩的多少,与接头刚度的大小成正比。这种方法是在围岩条件很好的情况下使用的,由于是非静定结构,将可以从围岩得到较大的地基抗力为前提进行计算。如果是软弱围岩,将产生很大的地基变形,使用时应慎重。3种方法总结(1)自由变形圆环法与弹性地基梁法计算的正最大弯矩位置相同，均在圆环顶部，负最大弯矩位置也相同(均位于拱腰)，弹性铰法计算的最大弯矩位于拱底。最大弯矩值则有一定差异，3种方法计算的正最大弯矩依次为:自由变形圆环法最大，弹性地基梁法次之，弹性铰法最小。而正最大弯矩处轴力则是弹性铰法最大，弹性地基梁法次之，自由变形圆环法最小。(2)3种方法计算的弯矩的差异主要是由这3种方法对管片接头和土层反力考虑的不同引起的。自由变形圆环法不考虑接头影响，而且仅在拱腰部分考虑了土层的被动抗力；弹性地基梁法也不考虑接头影响，但是在衬砌圆周合理地考虑了土层的抗力；而弹性铰法不仅考虑了接头影响(接头刚度的削弱)，还考虑了衬砌圆周土层的被动抗力。所以，从模型简化的角度来看，弹性铰法最合理，但是由于其对地层条件的应用限制，以及计算结果的准确性与接头刚度的取值有很大关系，而接头刚度又只能通过实验得到，所以弹性铰法在管片衬砌内力计算中通常起校核作用,提供参考或工程预估，实际工程设计中并不常用。10.7有限元法隧道衬砌结构杆系有限元分析的基木思想是，采用符合“局部变形原理”的弹簧地基来模拟围岩，而不采用连续弹性体，首先将衬砌与围岩所组成的衬砌结构体系离散化为有限个衬砌单元和弹簧单元所组成的组合体；其次确定由主动荷载(如围岩压力、衬砌自重等)直接作用在计算模型上所引起的变形，变形的轮廓就是衬砌与围岩相互作用区。这可以通过逐次渐进的方法求得。即先假定衬砌的某一段周边处在相互作用区，求出衬砌的变形轮廓，然后从没有相互作用的区域拿掉弹簧单元，再进行第二次计算;经过检查又去掉不起作用的弹簧单元，若在上次被拿掉弹簧单元的地方又发现新的相互作用，则必须再加上弹簧单元。如此反复进行，直到弹簧单元都正好放在相互作用的地方为止。杆系有限元分析的前提是结构的理想化，就是将结构看成为有限个单元的组合体，I、J单元之间仅在单元节点处相接，作用在结构上的外荷载和内力都只能通过节点进行传递，以节点力(轴力N、弯矩M、剪力Q)或节点位移(线位移、转角)代表整个结构的受力状态和变形状态。结构的理想化对于隧道衬砌结构来说，其结构理想化包括以下三部分内容。(1)衬砌本身理想化(2)围岩的理想化（3）荷载的理想化——等效结点荷载有限元分析的计算特点及步骤：（1）计算特点（2）计算步骤衬砌本身理想化对于隧道衬砌的内力，轴力和弯矩是主要内力。所以可将衬砌离散化为一些同时承受弯矩、轴力和剪力的轴心受压等直梁所组成的折线形组合体。衬砌单元的力学性质由弹性梁理论确定，即小变形、符合虎克定律。通过衬砌单元可传递弯矩、轴力和剪力。单元数目视计算精度的需要而定。对于一般的单线隧道，整个衬砌不少于16个单元；由于盾构隧道一般为圆形隧道，结构和荷载都对称，计算只需在一半衬砌上进行，可理想化为隧底是固定的，其水平位移、竖向位移和转角均为零。围岩的理想化把衬砌外围处于弹性抗力区范围以内的围岩看作是若干彼此互不关联的矩形岩柱。矩形岩柱底的一个边长是衬砌计算厚度(一般取b=1m)，另一个边长是相邻两离散化了的衬砌单元长度之半的和h。岩柱的深度与传递轴力无关，故不予考虑。岩柱与衬砌外缘紧密结合，在受力及变形中，相互制约、相互传力而共同作用。岩柱具有弹性地基的性质，采用局部变形理论的文克尔(E.Winkler)假定，以弹性抗力系数k来体现。荷载的理想化——等效结点荷载隧道衬砌所承受的外荷载，除了结构自重外，主要是围岩压力。把不同分布图形的分布荷载转换为等效结点荷载，作用在离散化衬砌结构的结点上，严格的说按静力等效的原则进行置换，即结点荷载所做的虚功应等于外荷载所做的虚功。但隧道所受荷载的准确性较差，故可按简单而近似的办法—即按简支分配的原则进行置换，这是衬砌计算中常用的方法。计算特点将整个结构看成是离散化单元的集合体，作用在结构节点上的荷载与结构节点位移之间的关系以结构刚度矩阵来联系。根据变形协调条件，联接于同一节点上的各单元在该节点上的单元节点位移相等，并等于该点的结构节点位移。又根据静力平衡条件，作用于结构节点上的节点荷载必须与该节点上作用的各个单元的单元节点力相平衡。因此，如果找到了单元节点力和单元节点位移之间的关系式，即可利用静力平衡条件建立结构刚度方程式，并利用边界条件求得结构节点的位移。从而再根据变形协调条件求得联结于该节点的各单元节点位移，再根据应力矩阵求得单元节点力(衬砌截面内力)。计算步骤-1(1)衬砌结构离散化。将衬砌结构离散为若干个衬砌单元和支承链杆单元。衬砌单元的数目可视精度要求和计算机容量而定。支承链杆单元视实际情况预先布置。计算各个衬砌单元的长度l及各节点处的x,y坐标值。定出各支承链杆单元的弹性抗力系数Ko。(2)将所有荷载换算成等效节点荷载。将荷载(包括竖向荷载、水平侧向荷载以及自重荷载)按结构坐标系的坐标轴方向化为节点荷载作用在节点上。(3)计算各衬砌单元的刚度矩阵(按局部坐标系)。(4)计算各支承链杆单元和弹性支座单元的刚度矩阵(局部坐标系和结构坐标系一致)。计算步骤-2(5)将按局部坐标系建立的衬砌单元刚度矩阵转换成结构坐标系的衬砌单元刚度矩阵。(6)形成结构刚度矩阵[K](7)按衬砌结构的边界条件(位移边界条件)修改结构刚度方程。(8)求节点位移的第一次近似值。(9)根据节点水平位移的正负值，判别是否去掉某支承链杆，并相应修正结构刚度矩阵。(10)计算节点位移的第二次近似值，并重复9,10步，直至所有支承链杆处的节点位移均符合条件为止。(11)求解最终的结构节点位移值。(12)计算局部坐标系中各衬砌单元的节点力。(13)进行衬砌单元截面强度校核10.8梁-弹簧模型结构模型