一种岩体随机结构面的三维离散元表征方法技术

本发明专利技术公开了一种岩体随机结构面的三维离散元表征方法，首先基于所需的岩体随机结构面三维网络数值模型范围，确定各随机结构面的等效半径，设置半径阈值，小于阈值的结构面不予以建模；然后对需要建模的结构面，按照等效半径从大到小的顺序在三维离散元软件3DEC中依次生成各结构面，并在生成前隐藏与该结构面不相交的块体；最后将同一个结构面分为真实结构面部分和虚拟结构面部分，对两部分赋予不同的力学参数，以精确界定结构面的真实边界。本发明专利技术保证一个真实结构面只需要切割一次，从而最大可能地减小了虚拟结构面数量，极大地提高了计算效率，克服了虚拟结构面变形参数难取值而导致的计算误差。

【技术实现步骤摘要】
一种岩体随机结构面的三维离散元表征方法
本专利技术涉及岩体结构建模
，具体涉及一种岩体随机结构面的三维离散元表征方法。
技术介绍
室内试验结果和工程实践经验均表明，岩体结构面对岩体力学、水力学性质及行为起决定性作用。其中，与工程尺度相当或略小的岩体随机结构面是应重点关注的问题。为了全面认识工程区岩体随机结构面，需要在随机结构面测量与分析的基础上对某些结构面特性或某项综合指标进行计算，这就要求对随机结构面进行三维网络模拟。目前，Monte-Carlo方法是进行随机结构面三维网络模拟的最有效方法，该方法本质上是统计分析的逆过程。随机结构面三维网络模拟还包含计算机图形可视化技术，由此不但增强工程区岩体结构的感性认识，也可通过布尔运算或其它自编算法方便地求出一些无法测量的几何特性或岩体结构综合指标。目前存在的问题是，随机结构面三维网络模拟的应用仅停留在几何或图形层面，由于要真正解决岩体工程问题无法绕开力学计算，因此研究如何将随机结构面三维网络模拟有效地转换为岩体力学计算模型具有十分重要的意义。当前，将岩体随机结构面在三维离散元软件(3DEC)中进行表征的方法主要是通过构建大量的虚拟结构面进行，即构建一个真实结构面需要增加四个虚拟结构面，造成极低的计算效率，且构建的真实结构面只能是四边形，从而难以与估算的结构面直径相关联。虚拟结构面数量的增加及规模的不准确性会导致以下结果：(1)岩体工程场地通常有成千上万、甚至超十万条结构面，若要构建部分较大规模的结构面，本就有庞大的结构面数量。若以增加四个虚拟结构面的方式来构建一个真实结构面，则大量的切割会造成海量的畸形块体，从而难以进行网格剖分。或者即使网格剖分成功，也会造成海量的单元和节点数据，从而导致极低的计算效率，且由于常应力单元、接触判断算法、接触力学算法等对稍不规则的单元带来的误差，会造成低精度的计算结果。(2)真实结构面规模的不准确以及虚拟结构面变形参数难取值问题，都会导致计算结果的失真，从而难以进行岩体工程应用。
技术实现思路
针对上述现有技术不足，本专利技术提出一种岩体随机结构面的三维离散元表征方法，最大程度地减少虚拟结构面数量及规模，并将真实结构面看作圆形且以能较准确获得的直径进行度量。为实现上述目的，本专利技术的技术方案是：一种岩体随机结构面的三维离散元表征方法，包括以下步骤：步骤1：基于所需的岩体随机结构面三维网络数值模型范围，确定各随机结构面的等效半径，设置半径阈值，小于所述阈值的结构面不予以建模；步骤2：对需要建模的结构面，按照等效半径从大到小的顺序在三维离散元软件3DEC中依次生成各结构面，并在生成前隐藏与该结构面不相交的块体；步骤3：将同一个结构面分为真实结构面部分和虚拟结构面部分，对两部分赋予不同的力学参数，以精确界定结构面的真实边界。进一步地，所述步骤1中各随机结构面的等效半径计算方法如下：基于Monte-Carlo方法生成的各随机结构面数据，根据所需的岩体随机结构面三维网络数值模型范围，首先判断结构面中心点是否在数值模型范围内，判断式如下：其中，JP为判断变量，Pi为中心点P的坐标值，Px、Py和Pz分别为中心点P的x、y、z轴坐标值，和分别为数值模型边界的最大值和最小值；设RP为中心点为P的结构面的半径，则其等效半径的求解方式如下：(1)JP＝0，则中心点P在数值模型内，该结构面等效半径(2)RP>JP>0，则该结构面等效半径计算如下：(3)JP>RP，则该结构面与数值模型不相交，该结构面等效半径进一步地，所述步骤2具体包括：步骤21：采用结构面切割命令，首先生成等效半径最大的结构面；步骤22：提取各块体的中心点坐标和顶点坐标，计算各块体中心点与顶点距离的最大值Rmax，以各个块体中心点为圆心、Rmax为半径确定块体的最大外接球；提取各块体表面的几何特征，计算各块体中心点与表面距离的最小值Rmin，以各个块体中心点为圆心、Rmin为半径确定块体的最小内切球；步骤23：若结构面中心点与某块体中心点的距离Lpc大于结构面等效半径Rpeqv加上块体最大外接球半径Rmax，则对该块体隐藏；若结构面中心点与某块体中心点的距离Lpc小于块体最小内切球半径Rmin减去结构面等效半径Rpeqv，则只显示该块体，隐藏其它所有块体；显示或隐藏块体命令执行后，开始切割生成该结构面；步骤24：用与步骤22、23相同的判断方法，按结构面等效半径从大到小的顺序依次生成其它结构面。进一步地，所述步骤3具体包括：设结构面的中心点为(x0,y0,z0)，半径为r0，遍历该结构面上所有的节点，当该节点满足以下公式时，则赋予真实结构面的力学参数，否则赋予虚拟结构面的力学参数：其中：a、b、c分别为该结构面节点的x、y、z轴坐标值。与现有技术相比，本专利技术具有有益效果：(1)保证一个真实结构面只需要切割一次，无需额外增加任何虚拟结构面，从而最大限度地减少了虚拟结构面数量；(2)按照结构面半径从大到小的顺序，以及采用块体外接球和内切球的判断方式，确保与结构面不相交的块体不进行切割，从而减小了虚拟结构面面积；(3)以力学参数区分真实结构面和虚拟结构面，能确保构建的真实结构面规模与现实相符，且能人为控制虚拟结构面力学参数对计算结果的影响。附图说明图1是本专利技术一种岩体随机结构面的三维离散元表征方法的流程图；图2是结构面等效半径计算模型示意图；图3是结构面与块体外接球、内切球的空间关系示意图；图4是同一个结构面区分为真实结构面部分、虚拟结构面部分的离散元表征图；图5是本专利技术一实施例中岩体随机结构面三维网络模拟图的东北等视图；图6是本专利技术一实施例中岩体随机结构面三维网络模拟图的东南等视图；图7是本专利技术一实施例中岩体随机结构面的离散元表征图。具体实施方式下面结合附图及实施例对本专利技术做进一步说明。如图1所示，一种岩体随机结构面的三维离散元表征方法，具体实施方法如下：(1)采用测线法或测窗法对岩体随机结构面进行采样，基于岩体随机结构面统计分析理论，获得随机结构面产状、直径与体密度等统计特征，确定随机结构面三维网络模拟的生成区与应用区范围，采用Monte-Carlo方法生成各结构面数据，基于所需的数值模型范围，采用下述方式确定各随机结构面的等效半径，并对结构面按照半径从大到小进行排序，设置结构面半径阈值，小于此阈值的结构面不予以建模；基于Monte-Carlo方法生成的各随机结构面数据，根据所需的数值模型范围，首先判断结构面中心点是否在数值模型范围内，其判断式如下：其中，JP为判断变量，Pi为中心点P的x、y、z轴坐标值，和分别为数值模型x、y、z边界的最大值和最小值；如图2所示结构面等效半径计算模型示意图，假定RP为中心点为P的结构面的半径，则其等效半径Rpeqv的求解方式如下：a)JP＝0，则该中心点在数值模型内，取该结构面等效半径Rpeqv＝RP；b)RP>JP>0，则该结构面等效半径的求解方法如下：c)JP>RP，则该结构面不可能与数值模型相交，取该结构面等效半径Rpeqv＝0。(2)对需要建模的结构面，按照等效半径从大到小的顺序在三维离散元软件(3DEC)中依次生成各结构面，并在生成前隐藏与该结构面不本文档来自技高网...

【技术保护点】
一种岩体随机结构面的三维离散元表征方法，其特征在于，包括以下步骤：步骤1：基于所需的岩体随机结构面三维网络数值模型范围，确定各随机结构面的等效半径，设置半径阈值，小于所述阈值的结构面不予以建模；步骤2：对需要建模的结构面，按照等效半径从大到小的顺序在三维离散元软件3DEC中依次生成各结构面，并在生成前隐藏与该结构面不相交的块体；步骤3：将同一个结构面分为真实结构面部分和虚拟结构面部分，对两部分赋予不同的力学参数，以精确界定结构面的真实边界。

【技术特征摘要】
1.一种岩体随机结构面的三维离散元表征方法，其特征在于，包括以下步骤：步骤1：基于所需的岩体随机结构面三维网络数值模型范围，确定各随机结构面的等效半径，设置半径阈值，小于所述阈值的结构面不予以建模；步骤2：对需要建模的结构面，按照等效半径从大到小的顺序在三维离散元软件3DEC中依次生成各结构面，并在生成前隐藏与该结构面不相交的块体；步骤3：将同一个结构面分为真实结构面部分和虚拟结构面部分，对两部分赋予不同的力学参数，以精确界定结构面的真实边界。2.根据权利要求1所述的一种岩体随机结构面的三维离散元表征方法，其特征在于，所述步骤1中各随机结构面的等效半径计算方法如下：基于Monte-Carlo方法生成的各随机结构面数据，根据所需的岩体随机结构面三维网络数值模型范围，首先判断结构面中心点是否在数值模型范围内，判断式如下：其中，JP为判断变量，Pi为中心点P的坐标值，Px、Py和Pz分别为中心点P的x、y、z轴坐标值，和分别为数值模型边界的最大值和最小值；设RP为中心点为P的结构面的半径，则其等效半径的求解方式如下：(1)JP＝0，则中心点P在数值模型内，该结构面等效半径(2)RP>JP>0，则该结构面等效半径计算如下：

【专利技术属性】
技术研发人员：曹洋兵，詹淦基，黄真萍，邱冬冬，陈玉华，曾焕接，陈俊熙，
申请(专利权)人：福州大学，
类型：发明
国别省市：福建,35