地下洞室群有限元分析地应力场计算方法.pdfVIP

第21 卷 第11 期 岩石力学与工程学报 21(12) ：1639～1644 2002 年11 月 Chinese Journal of Rock Mechanics and Engineering Nov. ，2002 地下洞室群有限元分析的地应力场计算方法* 1 2 3 杨 柯 张立翔 李仲奎 (1 温州大学建筑工程学院 温州 325027) (2 昆明理工大学工程力学系 昆明 650051) (3 清华大学水电工程系 北京 100084) 摘要 提出一种用于地下洞室群有限元计算的三维地应力回归分析的位移函数法。它在得到地应力场合理拟合的同 时，程式化地得到了施加于计算模型边界上满足变形协调条件的位移场。用这个位移场作为有限元计算的边界条件， 可避免附加约束对计算结果的影响，同时对地应力的大小和方向均可合理描述。最后给出了用该方法对二滩地下厂 房洞室群的计算结果。 关键词 地下洞室群，地应力，有限元 分类号 TU 94，TV 731.6 文献标识码 A 文章编号 1000-6915(2002)11-1639-06 因难在于边界上合理地假定剪应力(或平行该面方 1 问题的提出 向的位移)的大小和分布规律。目前经常被采用的一 种简化处理方法是，把边界面分成若干个区域，在 初始应力场是影响地下洞室群围岩应力、位移 每个区域上仅施加与该面垂直的应力或位移。虽然 分布以及破坏形式的重要因素之一。合理地得到建 这种做法对于减小拟合时人工调整边界条件工作量 模区域内的地应力场，并将其合理地施加于有限元 是有效的，但是难以保证平衡条件和协调条件同时 模型上进行计算分析，是所面临的两个基本问题。 满足。更重要的是，这种做法实际上是假设了在整 在实施地下洞室群三维有限元分析时，经常采 个区域内地应力的三个主方向平行于建模区域边界 用的两种地应力处理方法是：在模型的边界上施加 的法线方向。在很多实际情况中这种假设不能被认 外部荷载，或施加已知位移来获得建模区域内部的 为是合理的。例如在二滩工程中，三大洞室的上游 初始应力场。两种方法可以分别称为边界荷载调整 拱肩处均出现破坏，而下游拱肩则基本上没有破坏， 法和边界位移调整法。这两种方法所要解决的问题 这就明显与大主应力方向有关。文[1]也指出这是引 都是：如何调整边界荷载或边界位移，以便根据个 发洞室轴向位移的主要原因。总之，这两种分析方 别点的实测结果获得整个建模区域内的完整的地应 法，都难以避免在计算中只能选用边界的应力或位 力场。 移，它们的分布规律也是未知的，而需要加以某种 在建模区域任何一个边界面上，均存在 3 个方 [2] 假定 的缺陷。 向的应力或位移，那么6 个面上将有18 个待调整的 文[2，3]提出的应力函数法可以避免上述问题。 应力或位移。如果考虑边界面上应力或位移的非均 由于应力函数满足协调方程，所求得的应力，是满 匀分布，要调整的应力或位移分量还会大大增加。 足平衡方程和变形协调方程的。该方法能够很好地 组合所有这些分量来拟合内部的实测应力实在不是 解决 “获得合理的应力场”这一基本问题，利用它 一件轻松的