公路隧道地层结构法.ppt

有限元正分析 1 1 地下结构有限元理论基础地下结构有限元理论基础 1 有限元法的计算步骤 1) 离散和选择单元类型 2) 单元分析 3) 组装单元方程得出总体方程并引进边界条件 4) 解未知自由度(或广义位移) 5) 求解单元应变和应力 6) 解释结果 a. 选择位移函数 b. 定义应变位移和应力应变 关系 c. 推导单元刚度矩阵和方程 2 地层结构有限元计算方法 p 初始地应力的计算 p 本构模型 p 施工过程的模拟 n 构造应力 构造地应力可假设为均布或线性分布，并直接叠加在自 重地应力上得到初始地应力，计算式为 oo初始地应力的计算初始地应力的计算 n 自重应力 采用有限元方法或给定水平侧压力系数法计算。由有 限元法计算时，将自重荷载转化为等效结点荷载计算初始地 应力场。给定水平侧压力系数法是给定水平侧压力系数K0 ，按下式计算初始地应力： oo本构模型本构模型 n n 岩石岩石 oo线弹性模型线弹性模型 对于平面应变问题，横观各向同性体的弹性应力 增量可表示为: oo非线性弹性模型非线性弹性模型 采用邓肯-张模型的假设，应力-应变用双曲线关系近似描述，在主应力3不变时 轴向应变1和侧向应变3之间也假设是双曲线的关系 在不同应力状态下弹性模量的表达式为 不同应力状态下泊松比的表达式为 由Ei和i即可确定该应力状态下的弹性矩阵D。 oo弹塑性模型弹塑性模型 屈服准则 材料进入塑性状态的判断准则采用Drucker-Prager或MohrCoulomb屈服准则， 当采用Drucker-Prager屈服准则时 式中I1为应力张量第一不变量，J2为应力偏量的第二不变量，并有 对MohrCoulomb屈服准则，有 弹塑性矩阵 材料进入塑性状态后，其弹塑性应力应变关系的增量表达式： 弹塑性分析计算过程 增量时步加荷过程中,部分岩土体进入塑性状态后，由材料屈服引起的过 量塑性应变以初应变的形式被转移,并由整个体系中的所有单元共同负担。每 一时步中,各单元与过量塑性应变相应的初应变均以等效结点力的形式起作用 ，并处理为再次计算时的结点附加荷载,据以进行迭代运算,直至时步最终计算 时间,并满足给定的精度要求。 o各向同性材料的弹性应力增量可表示为 n n 梁梁 设梁在局部坐标系下结点位移为 对应的结点力为 则： n n 杆杆 设杆在局部坐标系下结点位移为 对应的结点力为 则： n n 接触面接触面 接触面采用无厚度节理单元，不考虑法向和切向的耦合作用，用增 量形式可表示为： oo施工过程的模拟施工过程的模拟 n n 一般表达式一般表达式 施工步是指一个相对完整的施工阶段，每一个施工步包含若干个增量步， 该施工步产生的开挖释放荷载在所包含的增量步中逐步释放，每一增量步释放 量由释放系数控制。 对各施工阶段的状态，有限元分析的表达式为： 对每个施工步，增量加载过程的有限元分析的表达式为： oo开挖开挖 开挖效应通过在开挖边界上设置释放荷载，并将其转化为等效结点力。表 达式如下： 开挖释放荷载采用单元应力法或Mana法计算。单元应力法先根据初始地 应力或前一步开挖相应的应力场求得预计开挖边界上各结点的应力，并假定 各结点间应力呈线性分布，然后反转开挖边界上各结点应力的方向(改变其符 号)，据以求得释放荷载。 oo填筑填筑 填筑效应包含两个部分，即整体刚度的改变和新增单元自重荷载的增加, 用如下公式描述： oo结构的施作与拆除结构的施作与拆除 3 地下结构有限元计算特点 pp 地下工程锚喷支护计算中有限元法的特点 p 单元类型选择和网格划分 p 计算范围的选取 p 边界条件和初始应力 p 开挖施工步骤的模拟 p 支护与衬砌的模拟 oo地下工程锚喷支护计算中有限元法的特点地下工程锚喷支护计算中有限元法的特点 n根据地下工程的支护结构与其周围岩体共同作用的特点，通常 把支护结构与岩体作为一个统一的组合体来考虑，将支护结构 及其影响范围内的岩体一起进行离散化。 n作用在岩体上的荷载是地应力，由自重地应力和构造地应力两 部分组成。 n通常把支护结构材料视作线弹性的，而岩体及岩体中节理的应 力应变关系作非线性，因而必须采用材料非线性的有限元法进 行分析。 n计算程序中一般应考虑开挖与支护的施工步骤的影响。 n由于地下工程一般轴线很长，因而通常可视作平面应变问题处 理。 oo单元类型选择和网格划分单元类型选择和网格划分 n单元类型的选择影响到计算的精度，贮存量的多少及运算时间 的长短，因而要求尽量选择合适的单元类型。 n通常在洞周附近区域，单元布置密些，而其他区域可疏些，但 也不宜疏密相差过于悬殊。 n单元边界应当划分在材料的分界面上和开挖的分界线上。 n一个单元内的边长不能相差过于悬殊，否则会增大误差。 n单元的结点应布置在荷载的突变点及锚杆的端点，便于锚杆和 载荷布置。 n单元划分应当充分利用对称性，以减少计算量。 n单元结点编号应注意到每个单元的编号序数尽量靠近，以减小 带宽长度。 oo计算范围的选取计算范围的选取 n在岩体中开挖洞室，应力重分布的范围是有限的，因而计算的 范围也可以是有限的。根据圆洞计算的理论解以及计算的经验 ，一般可取计算范围大于3倍洞径，即L 3D，如取3.5 D或更 大一些。对于高边墙洞室，则应取 L 3H（H洞高）。 oo边界条件和初始应力边界条件和初始应力 n所取岩体边界上的位移边界条件通常两侧边界取按水平方向固 定，铅直方向自由。下边界约束情况一般按铅直方向固定，水 平方向自由，在两个角点上两个方向固定。 oo开挖施工步骤的模拟开挖施工步骤的模拟 n计算结果表明，按一步开挖与多步开挖计算其结果是不同的， 因而计算中应当模拟开挖的步骤。同时它能使我们了解到各施 工阶段围岩的应力变化情况，有助于我们对施工方法的选择。 oo支护与衬砌的模拟支护与衬砌的模拟 用于模拟支护结构与衬砌的单元，通常有杆单元、梁单元和实体单元。 杆单元通常用于模拟“柔性”支护，如锚杆(点锚式锚杆或预应力锚杆、 一般砂浆锚杆、及预应力全长粘结式锚杆)，厚度不大的喷层等。 梁单元通常用于模拟“刚性”支护，如二次衬砌等永久结构。 实体单元或梁单元用于模拟衬砌及喷层。在现浇衬砌与岩体之间可能产 生相对滑移，这时本衬砌与岩体之间需要加上接触面单元，允许两者之间 产生相对滑移。 地下工程的有限单元法计算中，应当考虑开挖面的空间效应和支护的施 作时间。 刚开挖时，锚喷支护或衬砌尚未施作，支护不可能发挥作用，只有当支 护施作后，方能逐渐发挥作用，因而支护的施作时间对计算的结果有重大影 响。通常计算中假定洞室开挖后，支护施作之前已释放了一定次数的增量荷 载。 数值计算中如何合理地模拟锚杆的作用，至今还没有很好的方法，一般 认为，围岩加锚以后，锚固区围岩C，值应适当提高，但究竟提高多少，主 要是提高C值还是提高值至今尚无一致的看法，计算中应当注意，锚杆是点 加荷的，作平面应变计算时，应按锚杆间距折算为线分布的，否则会夸大锚 杆的作用。 根据目前的计算结果，一般计算得到的锚杆支护作用没有实际的 锚杆支护作用大。 2 2 曙光软件有限元算例演示曙光软件有限元算例演示 1 ) 工程概况 工程地质剖面图 隧道结构图 2) 隧道施工全过程受力分析 启动曙光软件，选择加载计算模块 3 3 建立模型建立模型 1. 有限元模型(网格)的试生成，以检查模型中不合格的地方。 2. 重新寻找封闭面，模拟施工过程。 3. 生成有限元模型。 4. 有限元参数的填入。 5. 计算。 6. 计算结果解释。 7. 生成计算报告。 围岩：II类，E：2GPa,取100MPa, =0.350.45，取0.4 c取100kPa 自重：