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

支承链杆单元刚度矩阵。支承链杆单元i的弹性抗力三；与其端点压缩位移民的关系按文克尔假定可写为 : 也即： 衬砌单元与支承链杆单元 （４－２１） （４－２０） * 支承链杆单元的刚度矩阵为： 上图所示的衬砌单元与支承链杆单元 简写为： 墙底弹性支座单元 （４－２４） （４－２２） * （3）结构刚度方程的形成 现将式（4—21）改写成如下矩阵形式： （４－２５） （４－２６） （４－２７） * 　　根据静力平衡条件：在同一节点上作用的各单元节点力的总和应与该节点的节点荷载相平衡，因此有： （４－２８） （４－２９） * （4）衬砌内力计算 在求得最终的结构节点位移值后，便可求解衬砌的内力。 　　根据变形协调条件，结构节点位移与汇交于此节点的单元节点位移相等，即： （４－３０） （４－３１） （5）隧道衬砌结构杆有限元分析基本步骤如下： 　①衬砌结构离散化。 　②将所有荷载换算成等效节点荷载； 　③计算各衬砌单元的刚度矩阵（按局部坐标系）； * ④计算各支承链杆单元的弹性支座单元的刚度矩阵（局部坐标系和结构坐标系相一致）； ⑤把局部坐标系建立的衬砌单元刚度矩阵转换成结构坐标系的衬砌单元刚度矩阵；⑥形成结构刚度矩阵［K］； ⑦按衬砌结构的边界条件（位移边界条件）修改结构刚度方程； ⑧求节点位移的第一次近似值，并检查和取消位移ui＜0的支承链杆单元，修改刚度矩阵，重新计算； ⑨据节点水平位移的第二次近似值，重复⑧、⑨步骤，直到所有支承链杆处的节点位移均符合ui＞0的条件为止； * ⑩四求解最终的结构节点位移值； （１１）计算局部坐标系中个衬砌单元的节点位移值； （１２）计算拱顶截面支座约束反力． （4－32） （4－33） * 五、隧道衬砌截面强度检算 　　检算所得的安全系数是否满足相应规范所要求的数值 ： * 　拱形隧道和明洞衬砌属偏心受压构件，其截面极限承载力Nz的计算，按轴力的偏心距e＝M／N的大小有两种情况： （1）抗压强度控制（当 e＜0.20 d）： （2）抗拉强度控制（当 e＞0.20 d）： * 第二节 岩体力学方法 整体复合模型 １.围岩本身由于支护结构提供了一定的支护抗力，而引起它的应力调整，达到新的稳定； ２.由于支护结构阻止围岩变形，也必然受到围岩给予的反作用力而发生变形。 目前对于这种模型求解方法有： 　　解析法、数值法、特征曲线法、剪切滑移破坏法 * 一、解析法 　1．基本假设 （１）深埋圆形隧道，无限长； （２）原岩应力各向等压； （３）围岩为理想弹塑性体，符合库仑强度准则。 　2．基本方程 * 弹性区：积分常数待定的弹性应力解 塑性区：平衡方程 * 强度准则方程——库仑准则 塑性区：平衡方程 * 3．边界条件 　　在弹性区内，外边界处有 r→∞，σr＝σθ＝P0，在弹性区的内边界（交界面）r＝RP（塑性区半径）。 * 在塑性区的外边界（交界面） 在内边界（周边）处： σr =0（不支护），σr =P（支护反力） * 4．问题的求解 式（4—38）、式（4—39）联解并利用式（4—42）有 （4－43） （4－4４） 式（4—43）代人式（4—39），整理得： * 由（式4—37）及弹性边界条件有 （4－45） 由式（4—43）、式（4－45）、塑性内边界条件式（4－41） （4－4６） * （4－4７） 解的B，再代人式（4—37），整理得弹性区应力： 由式（4—41）、式（4—44）、式（4—47）又有 ＝ * （4－4８） 得塑性区半径 5．结果的列出 弹性区应力： （4－4９） * 塑性区应力 （4－５１） （4－５０） 塑性区半径 （4－５２） * 　　若式（4—38）、式（4—39）联解后，改用边界条件σr =P决定积分常数，则可得有支护反力P情况下的弹性区应力为： （４－５３） 塑性区应力为 （４－５４） （４－５５） * 塑性区半径 （４－５６） 将式（4—55）改写为 （４－５７） * 6．讨论 Rp与R0成正比，与Ｐ0成正变，与Ｃ，Ｐ成反变关系。 塑性区应力与Ｐ0原始应力无关； 支护反力Ｐ＝０时， Rp最大。 * 二、数值方法 数值计算方法，通常包括： 　　　　　　有限差分法（离散控制方程） 　　　　　　有限元法（离散解题全域） 　　　　　　边界元法（只离散边界） 　　　　　　离散元法（块体运动） * （一）有限元法 图4—13 有限元分析系统的结构图 * （二）地下工程有限元特点 1．计算范围的选取（划域） 2．单元划分与选择 单元节点编号顺序 3.　开挖效果的模拟 ４.　求单元应力 * （三）　典型数值计算程序介绍 １.　ＦＬＡＣ３.３软件系统 ２.　ＡＤＩＮ