孙家崖隧道施工有限元模拟分析报告.doc

孙家崖隧道施工 有限元模拟分析报告 兰州交通大学 中铁七局集团第三工程有限公司 2010年12月  目 录 第一章 有限元分析的理论和计算方法 1 1 材料的弹塑性应力应变关系 1 2 屈服准则 1 3 体系有限元离散化和平衡方程的建立 2 4 荷载列阵的计算 3 5 平衡方程的求解 3 6 衬砌内力计算 4 7 材料强度发挥系数 4 8 衬砌截面强度检算 5 第二章 孙家崖公路隧道右线衬砌和围岩体系的有限元分析 6 1 孙家崖隧道台阶法有限元模拟分析 6 1.1 有限元分析模型 6 1.2 有限元分析结果 7 1.3 小结 16 2 孙家崖隧道单侧壁导坑法有限元模拟分析 18 2.1 有限元分析模型 18 2.2有限元计算结果分析 19 2.3 小结 26 3 两种开挖方式的对比与分析 28 第三章 孙家崖公路隧道左线衬砌和围岩体系的有限元分析 32 1 孙家崖隧道台阶法有限元模拟分析 32 1.1 有限元分析模型 32 1.2 有限元分析结果 33 1.3 小结 40 2 孙家崖公路隧道单侧壁导坑法有限元模拟分析 41 2.1 有限元分析模型 41 2.2有限元计算结果分析 42 2.3 小结 49 3 两种开挖方式的对比与分析 50 第四章 孙家崖公路隧道锚杆施做方式的有限元分析 53 一、隧道径向锚杆法模拟 53 1.1 有限元分析模型 53 1.2 有限元分析结果 54 1.3 小结 62 二、隧道锚杆与水平夹角150方向法模拟 62 2.1 有限元分析模型 62 2.2 有限元分析结果 63 2.3 小结 71 三、对比分析与结论 72 3.1 位移的对比 72 3.2 初衬应力的对比 72 3.3 锚杆轴力的对比 73 3.4 结论 73 第一章 有限元分析的理论和计算方法 在有限元分析中，采用岩土工程二维弹塑性有限元法进行衬砌与围岩的静力分析，其要点简述如下。 1 材料的弹塑性应力应变关系 采用增量弹塑性应力应变关系 (1.1) 式中 应力向量为： (1.2) 应变向量为： (1.3) 弹塑性矩阵为： (1.4) 弹性矩阵为： (平面应变问题) (1.5) ——材料的弹性模量和泊松比 ——材料硬化参数 ——屈服函数 2 屈服准则 采用Mohr-Coulomb屈服准则，即： (1.6) 式中 ——应力的第一不变量； ——偏应力的第二不变量； ——罗地角； ——材料的粘聚力和内摩擦角。 3 体系有限元离散化和平衡方程的建立 采用8节点等参数有限元将围岩与衬砌离散化，单元节点力与节点位移之间的关系为： (1.7) 式中 ——单元刚度矩阵。 (1.8) 其中 分别为节点的水平和竖向位移。 单元刚度矩阵按下式计算： (1.9) 式中 ——单元的几何矩阵，对于8节点等参单元： (1.10) 其中 () (1.11) ——形函数。 将各单元对体系的贡献迭加，即可得到体系的总刚度矩阵： (1.12) 联系体系节点位移与外荷载的平衡方程为： (1.13) 式中 ——体系的节点位移向量； ——荷载列阵。 4 荷载列阵的计算 用释放荷载法进行计算。 在有限元网格上边界和左右边界上作用有根据自重力算出的边界荷载，侧压力系数ξ=0.778。根据作用在单元边界上的分布荷载，按虚功原理可算出等效节点荷载为： (1.14) 式中 、——法向和切向的分布荷载； ——加荷边长度。 围岩和衬砌自重节点力可按下式计算： (1.15) 式中 ——体力，其中分别为体力的水平和竖向分力。通常：。 体系总节点荷载由上述几项构成。 5 平衡方程的求解 采用切线刚度法进行迭代计算，每次迭代中采用高斯消元法求解体系的平衡方程(1.13)。 采用节点位移增量的收敛准则判断求解过程的收敛： (1.16) 式中 ——第次迭代所得的第节点的位移增量； ——第节点的总位移； ——体系的自由度总数； ——收敛容许误差百分数。 在求解的过程中，体系残余不平衡节点力得到消除。 对于超出屈服面的应力进行调整，使之回到屈服面上。 6 衬砌内力计算 求得衬砌各单元中的各个高斯点的应