考虑参数分布特征隧道围岩稳定失效概率计算.pdfVIP

第27卷第1期 　 计 算 力 学 学 报 　 Vol . 2 7 ,No . 1 20 10 年 2 月 Chinese Journal of Computational Mechanics February 20 10 文章编号 (20 10) 0 10 12007 考虑参数分布特征的隧道围岩稳定失效概率计算 苏永华 , 　李 　翔 , 　赵明华 , 　谢志勇 (湖南大学 土木工程学院 地下工程系 ,长沙 4 10082) 摘 　要 :基于围岩变形破坏原理 ,建立了围岩稳定极限状态方程的普遍形式 ,研究了响应面方法求解该方程的过 程 ,得出经典响应面方法仅仅适用于围岩物理力学参数随机分布曲线偏度系数为零的情况 。针对围岩物理力学 参数统计分布曲线的复杂性 ,利用统计矩参数对随机变量分布曲线形态特征的控制作用 ,提出了抽样样本修正方 法和具体计算公式 ,归纳了修正后响应面算法的运行程序 。以某公路隧道为例 ,分别采用经典算法和修正后算法 计算了围岩稳定可靠度 ,然后与蒙特卡洛模拟结果进行了对照 ,经典方法计算结果与准精确解的绝对误差为 10 17 % ,修正方法的计算结果与准精确解的绝对误差为 0388 % 。 关键词 :计算力学 ;隧道变形 ;抽样方法 ;失效概率 ;响应面法 中图分类号 : TU4 136 ;O3 17 + . 3 　　　文献标识码 :A 1 　引 言 2 　基于围岩收敛的结构功能函数 2 1 　地层与支护特征曲线 目前在各类地下结构的信息化设计 、施工中 , 如图 1 所示, 图中曲线 1 是地层围岩压力随坑 主要考虑地下结构所在环境的围岩类别 、埋深等因 道周边位移的函数, 分为两个阶段 :A B 段属于变 素 ,根据地层特征理论曲线 ,采用最大允许变形值 [ 1 ,2 ] 形压力阶段, 随着坑道周边位移的释放, 围岩压力 作为围岩稳定和支护设计的依据 。但是 , 围岩 逐渐减小; 过 B 点后, 围岩出现松动, 压力增大 。从 的实际变形决定于地应力场和岩体的物理力学性 图中可以看出, 当坑道周边位移达到其极限位移 质 ,而岩体是千百万年历史的地质体 ,其物理力学 umax 时, 围岩产生的压力最小, 此时施作衬砌需要 性质及赋存的地应力场具有随机性 ,导致围岩的实 提供的支护抗力最为合理 。曲线 2 和 3 分别为支护 际变形也具有随机性 。特别是在复杂地层中 ,围岩 结构特征曲线, 不同的刚度具有不同的形状, 分别 物理力学参数随机性程度大 ,其变形无法通过岩体 表示在坑道围岩周边的初始位移达到 u0 1 和 u02 时 物理力学参数的明晰函数表达出来 ,所以在地下工 程围岩稳定可靠度分析中 ,响应面算法成为重要的 求解方法之一[35 ] 。本文首先研究了基于围岩变形 的极限状态方程建立准则及其响应面求解方法 ,然 后针对经典响应面算法中采用的基本变量抽样方 法无法反映岩体物理力学参数分布曲线特征的致 命缺陷 ,探讨了抽样的改进方法及相应的算法 , 以 提高响应