数值方法作业概要1.docVIP

河南理工大学 土木工程数值方法（作业） 题 目： 姓 名： 学 号： 学 院： 专 业： 班 级： 指导教师： 完成日期： 年 月 日 一、工程概况 1 二、建模方法和求解 1 2.1、建立网格 1 2.2、本构模型和参数的选择 2 2.3、模型求解 2 三、数据分析 3 3.1、垂直应力 3 3.2、垂直位移 4 3.2.1、云图 4 3.2.1、巷道中心截面切片 5 3.3、最大主应力 7 3.3.1、云图 7 3.3.2、巷道中心截面切片 8 3.4、塑性区 9 3.5、特征点监测 9 四、结果和讨论 10 一、工程概况 模拟项目为一煤矿巷道，煤层埋深500m。本次模拟的目的是分析在巷道高度恒定条件下，不同宽度条件下的应力与位移的变化规律，以便为巷道的设计和支护方法的选择提供建议。通过建立模型，在初始地应力平衡之后采用null模型来观察和记录特征点的应力和位移，比较不同方案的演化规律。支护采用cable模型。 二、建模方法和求解 2.1、建立网格 建模的思路为采用radcylinder网格和radtunnel网格拼接的方法，模型的尺寸为50m×50m×50m。通过设置size数的匹配来实现拼接过程的节点的自动融合，而后通过镜像的方法建立整个模型。对于不同宽度的模型，网格的基本结构保持不变，只需要设置不同的dimension数组，如表2.1所示。将模型填充并定义组名，这样在开挖的时候就有明确的对象。以B=6m为例，建立好的网格如图2.1所示： 宽度（m） 网格 Dimension 宽度（m） 网格 Dimension 5 radcylinder 2.5，3，2.5，36 radcylinder (3，3，，3 radtunnel (2，，， radtunnel (2，，，7 radcylinder (3.5，3，，38 radcylinder (4，3，，3 radtunnel (2，，， radtunnel (2，，，9 radcylinder (4.5，3，.5，310 radcylinder (5，3，，3 radtunnel (2，，， radtunnel (2，，， 图2.1 网格建立 图2.2 初始地应力生成 2.2、本构模型和参数的选择 首先采用弹性模型来求初始地应力场，然后将所有位移和速度清零之后换做莫尔-库伦模型求解。通过在模型的五个面施加应力的方法来实现深埋隧道的初始地应力场，生成的初始地应力场如图2.2所示。岩体的特性参数要综合考虑各个岩层的属性，最终做出较为合理的选择。 2.3、模型求解 采用莫尔-库伦模型求解，在求解过程中记录三个特征点的位移，分别是巷道中心截面（y=25m处）的顶板下沉、底鼓以及右侧墙体的位移，每50步记录一次。采用固定时步的方法，固定2000步求解。求解完成后显示竖向位移、最大主应力以及塑性区的范围，数据保存后可以导入到Tecplot软件中进行处理。 采用cable模型进行支护，锚杆长3m自由段和锚固段各1.5m。在巷道中心截面的顶板处设置锚杆，如图3所示。 图2.3 锚索结构支护 三、数据分析 3.1、垂直应力 垂直应力随着宽度的增加影响范围逐渐扩大，高应力区域主要集中在巷道两侧，这是由于巷道开挖后上部岩层的重力主要靠巷道两侧的煤体来承担。 图3.1 垂直应力变化云图 3.2、垂直位移 首先观察FLAC3D中显示的云图的变化规律，随后观察经过Tecplot软件处理的隧道中心截面的切片的演化规律。 3.2.1、云图 图3.2 垂直位移变化云图 垂直位移主要发生在顶板和巷道底部，分别为底部的正向位移和顶板处的负向位移，第3.4节将通过监测数据更加形象地说明这个规律。 3.2.1、巷道中心截面切片 图3.3 垂直位移切片图 通过观察切片的等值线图，我们发现随着宽度的增大，位移场的范围逐渐增大，位移峰值也逐渐增大。这表明扁平状截面的巷道应力和位移更大，需要更加充分的支护形式来保证其安全。 3.3、最大主应力 和垂直位移一样，我们观察云图和切片。 3.3.1、云图 图3.4 最大主应力云图 我们发现应力影响的范围逐渐增大，高应力区域的范围和应力峰值也逐渐增大。 3.3.2、巷道中心截面切片 图3.5 最大主应力切片图 通过将数据导入Tecplot软件中进行切片观察，我们发现，高应力区随着宽度的增加逐步扩大，最大主应力的峰值也逐渐增大。