08.4 施加荷载和数值模拟结果

 1 施加荷载和数值模拟结果

加载第二步的计算数据，并初始化模型。因上一步沉积地层时将单元连接全部断开，为此在该步中首先使用d.mo.bFilter(:)=true命令连接全部单元，并使用d.deleteConnection('boundary')命令断开样品和边界的连接。因为第二步中为模拟真实世界中开挖隧道时地面上存在的荷载，在模型中施加等效的荷载。为此在该步中需要重新施加压力，以防止模型发生爆炸。使用fs.setPlatenStress命令施加与第二步中相同的压力后再对模型进行标准平衡。

初始设置完成后，筛选出上压力板的部分区域，大致思路是：1.得到上压力板的单元id矩阵tp；2.提取出每个单元所对应的x坐标矩阵tpX；3.通过条件判断得到过滤器矩阵blockFilter，即以上压力板中心为基准左右4m内的颗粒；4.通过tp(blockFilter)获取满足条件的上压力版单元id矩阵blockId。

在此基础上通过addGroup()命令以blockId建立新组接着用blockForceZ=d.mo.mGZ(blockId)提取block单元的Z方向体力。当在topPlaten上施加压力时，实际上都是在其单元上施加了一个体力，即对应mGZ矩阵，体力如果是重力,其方向就往下。在做完这些处理后，程序可以通过for循环模拟体力增加之后产生的作用。在for循环中，会逐渐增加block单元Z方向的体力，通过多次循环会使这些单元的体力不断地增加，这样就实现了在模型中施加递增的压力。