ANSY和FLAC3D模拟隧道沉降对建筑物的影响教程详解.doc

本文采用大型有限元软件ANSYS建立分析模型，利用APDL（ansys参数化编程语言）编写命令流，实现房屋、隧道以及土层几何模型的建立，利用mesh模块划分网格，本文采用自由网格划分，自动生成四面体单元，一共10497个节点，30727个单元。 本次模型水平Y方向长度为200米，水平X方向宽为1米，竖直Z方向高度为100米。 房子长40米，宽1米，高40米，其中地下部分15米，地上部分25米。 隧道直径10米，上部隧道顶采用近视半圆处理，下部仰拱设置为弧形。隧道顶距离地面约15米，隧道底距离地面约22.5米。 土层基本参数如下表所示： 0~-3m为素填土 -3~-10m为粉质粘土 -10~-40m为强风化花岗岩，隧道处于该地层-15~-22.5m -40~-100m为弱风化花岗岩 弹性体积模量 Bulk（Pa） 弹性剪切模量 Shear（Pa） 抗拉强度 Tension（Pa） 内聚力 Cohesion（N） 内摩擦角Friction(°) 密度 Density 房屋 1.5e10 1.2856e10 2500 素填土 5.63e7 5.83e6 0 1.9e4 14 1500 粉质粘土 1.5e7 6.92e6 0 2.5e4 25 1700 强风化花岗岩 3.241e9 1.326e9 4.0e6 4e5 32 2300 弱风化花岗岩 3.894e9 2.008e9 4.0e6 7e5 37 2500 通过ANSYS-FLAC3D接口程序，输出单元和节点信息到FLAC3D文件，导入文件到FLAC3D进行计算和后处理，计算结果通过FLAC3D-Tecplot接口程序，进一步导入到Tecplot中做后处理，最终结果如下图所示： 注：本文房屋采用model elastic弹性本构模型分析，可以较好的模拟钢筋混泥土结构的特性，土层及隧道采用model mohr摩尔库伦本构模型，可以较真实的模拟岩土体材料的特性。开挖采用设置null单元方法进行模拟，隧道支护结构采用shell壳型结构单元进行模拟。 未开挖前的状态 图1 未开挖隧道时房屋和地基在自重作用下的变形云图 图2 未开挖隧道时房屋和地基在自重作用下的垂直变形云图 图3 未开挖隧道时房屋和地基在自重作用下的y方向水平变形云图 图4 未开挖隧道时房屋和地基在自重作用下的x方向水平变形云图 图5 未开挖隧道时房屋和地基在自重作用下的垂直应力分布图 图6 未开挖隧道时房屋和地基在自重作用下的SYY水平应力分布图 图7 未开挖隧道时房屋和地基在自重作用下的SXX水平应力分布图 图8 未开挖隧道时房屋和地基在自重作用下的SYZ剪应力分布图 图9 未开挖隧道时房屋和地基在自重作用下的最大主应力分布图 图10 未开挖隧道时房屋和地基在自重作用下的中间主应力分布图 图11 未开挖隧道时房屋和地基在自重作用下的最小主应力分布图 图11 未开挖隧道时房屋和地基在自重作用下的塑性区分布图 图12 未开挖隧道时房屋和地基在自重作用下的变形矢量图 图12 未开挖隧道时房屋在自重作用下的变形矢量图 隧道上半部开挖 图1 开挖隧道上部后的变形云图 由图可知：开挖后，最大变形发生在隧道底部，最大变形值为0.634mm，隧道顶部最大变形值约为0.35mm。受到隧道开挖的影响，临近房屋在不同程度上发生变形，其右上角的变形值较大，约0.1mm，向下变形值逐渐减小，房屋左下角的变形值较小，几乎可以忽略不计。 图2 开挖隧道上半部后的垂直变形图 隧道左半段开挖后，周围土体和房屋垂直方向发生变形，隧道顶部的土体产生沉降较为明显，最大值约0.4mm，变形竖直向上逐渐减小；隧道底部产生隆起变形较为明显，最大值约0.63mm，变形值竖直向下逐渐减小。房屋由于受到开挖的影响，靠近隧道的一侧变形较大，约0.1mm，变形值向远离隧道的方向递减。 图3 开挖隧道上半部后的Y方向变形云图 由以上y方向附件变形等值线图可知：开挖后，隧道左侧产生的附加水平变形方向与y轴方向相反，最大值约0.1mm，右侧产生的附加水平变形与y轴方向相同，最大值发生在隧道右侧壁上，约为0.2mm 图4开挖隧道上半部后的X方向变形云图 图5 隧道开挖上半部后的垂直应力分布图 图6 隧道开挖上半部后的Y方向水平应力分布图 图7 隧道开挖上半部后的X方向水平应力分布图 图8 隧道开挖上半部后的τyz剪应力分布图 图9 隧道开挖上半部后的最大主应力分布图 图10 隧道开挖上半部后的中间主应力分布图 图11 隧道开挖上半部后的最小主应力分布图 图11 隧道开挖上半部后的塑性区分布图 Shear-n 正处于破坏状态或正位于屈服面上 Shear-p 弹性，但之前发生过剪切破坏 Tension-n 正在拉伸破坏 Tension-p 弹性，但之前发生过拉