岩体力学结课报告2.docVIP

岩 体 力 学 结 课 报 告 班级： 姓名： 学号： 摘要: 大温差变化对隧道衬砌受力的影响在隧道计算与设计中是一个不能忽视的问题。以武夷山隧道衬砌为例，考虑温度荷载影响，采用有限元法对隧道衬砌结构的变形和受力进行了计算分析，得到了温度应力影响下隧道衬砌和围岩的变形和受力一些规律，此结论可为考虑温度应力影响下隧道衬砌设计与施工提供一定参考价值。 引言： 由于隧道内二次衬砌表面温度及二次衬砌背后一定深度的围岩体温度差的存在，在混凝土衬砌内部会产生压应力，表面会产生拉应力。而大温度梯度会引起较大的表面拉应力或者收缩应力，可能会在混凝土表面产生表面裂缝或收缩裂缝，对衬砌结构带来严重的危害[1~4]。因此，在隧道衬砌设计与施工中有必要对考虑温度影响下的隧道衬砌受力规律进行分析研究。 本文以新建向塘至莆田铁路项目武夷山隧道为例，考虑温度荷载影响，采用有限元法对隧道衬砌结构的变形和受力进行了计算分析，为考虑温度应力影响下隧道衬砌设计与施工提供一定的参考依据。 基于温度应力的有限元分析法 1、单元特性 在对结构离散化后，对三角形常应变单元进行特性分析，它是以位移为基本未知量，即用结点位移表示单元体的位移、应变和应力，通常要建立单元结点位移与单元内任意一点位移之间的关系，图1所示三角形常应变单元，其中结点按逆时针排列，整个单元有六个结点位移分量。 图1 三角形常应变单元 在三角形常应变有限元法中，坐标变换式和位移模式选择最简单的线性函数形式。利用平面问题的几何方程得出应变分量，用结点位移表示单元应变的矩阵方程为 （1） 式中，为单元应变矩阵，为三角形的面积。按照本构关系，将各向同性材料弹性矩阵代入，应力的计算公式是 （2） 由于和形函数位移系数、、、、、都是常量，因而单元中任意一点的应变分量也都是常量，故通常称为这种单元为常应变单元。 2、单元刚度矩阵 由系统方程，计算三角形单元刚度矩阵为 （3） 式中，是单元厚度，是常数。将应变矩阵公式（1）和平面应力问题的弹性矩阵公式（2）代入公式（3）中，得到常应变三角形单元的刚度矩阵的显式为 （4） 式中 3、自重应力和温度应力 设单元厚度为，单元体积密度为，自重方向为y负方向，则均质等厚度平面三角形单元自重引起的等效结点力为 （5） 考虑温度荷载产生的初应力场在三角形单元是线性分布的，则单元各结点上加上温度荷载引起的等效结点力为 （6） 计算应力的公式相应地修改为 （7） 实例分析 1、工程概况 新建向塘至莆田铁路项目武夷山隧道地处江西省南丰县及福建省西北部建宁县境内，为单洞双线隧道。隧道全长14658m，隧道最大埋深约348m。隧道位于剥蚀中低山区，测区内出露地层较简单，地质构造一般，水文地质条件良好，无不良地质。隧道轨面以上有效净空面积按不小于92m2设计，其中断面设计采用《通隧（2005）0201》通用图。武夷山隧道围岩分级为Ⅴ级的设计长度有935m，占设计全长6.38%。根据相关技术规范，Ⅴ级围岩断面设计详见图2。 图2 Ⅴ级围岩衬砌断面图 2、计算模型 利用平面应变假定、变分法和最小势能原理，分析围岩和衬砌在其自重以及衬砌内外温差作用下的变形和应力分布。分析中采用三角形平面单元。按照图纸设计尺寸，利用前处理工具选择单元形状和分割方案确定单元和结点的数目，首先建立关键点，然后采用六心圆法构造计算模型，并选择单元进行定义属性和网格划分。为保证衬砌以及附近的网格密度大于围岩部分的网格，对整个计算模型进行分区细化处理，以满足区域计算精度和节约计算时间的要求。为简化计算，模型中围岩和衬砌均采用各向同性材料，并且没有考虑锚固支护作用。整个计算模型共分为1953个结点，3604个单元。其中，除上层岩体边界和隧道断面内为自由边界外，其他三个边界均施加法向第一约束，如图3所示,其中深色部分为混凝土衬砌，浅色部分为围岩，左右及下方位置为约束边界。相关计算参数见表1所示。 图3 计算模型示意图 表1 模型相关参数 T模型 弹性模量E0(kPa) 泊松比ν 单元厚度h(m) 密度(kN/m3) 线膨胀系数α 岩体 8.50e5 0.3 1.0 19 1e-5 衬砌 2.81e6 0.2 1.0 25.5 5e-6 3、计算结果分析 本论文假定围岩温度在分析过程中始终不改变，在衬砌的自由边界设定-20～20℃之间等差为5℃温度值，即取9种温度值进行分析，以研究衬砌在隧道内温度变化下的应变和位移变化特征。图4~6为不考虑温度应力的计算结果，图7~9为考虑温度应力的计算结果。 图4 y方向位移云图（不考虑温度）