结构面抗剪区域的确定方法与流程

本发明属于岩石力学工程领域，涉及一种结构面抗剪区域的确定方法。

背景技术：

大量的工程实践证实，很多工程岩体的失稳问题都与结构面存在有关，研究结构面的力学行为具有十分重要的工程意义。在结构面的剪切试验或原位抗剪试验中，确定其表面上能够提供抗剪力的区域对分析结构面的变形特征、破坏特征以及剪切机理是十分重要的。学者grasselli通过对结构面剪切试验后破坏特征观察发现只有面向剪切方向上的区域才有可能提供抗剪力(grassellig,wirthj,eggerp.quantitativethree-dimensionaldescriptionofaroughsurfaceandparameterevolutionwithshearing[j].internationaljournalofrockmechanicsandminingsciences,2002,39(6):789-800.即grassellig，wirthj，eggerp.粗糙节理面的三维定量描述及其随剪切过程的演化[j]。国际岩石力学与采矿科学杂志，2002,39(6)：789-800.)。然而当结构面表面的形貌特征十分复杂时，通过直观地观察是很难确定其抗剪区域的，这对研究结构面的剪切机理是不利的，因此，准确地确定结构面的抗剪区域具有十分重要的意义。

技术实现要素：

为了克服现有技术中无法准确地确定结构面的抗剪区域的问题，本发明提供一种结构面抗剪区域的确定方法，该方法能够准确地确定结构面表面够提供抗剪力的潜在区域，进而为进一步分析结构面的变形特征、破坏特征以及剪切机理提供基础。

为了实现上述目的，本发明所采用的技术方案是：

一种结构面抗剪区域的确定方法，包括以下步骤：

1)利用三维扫描仪对结构面进行扫描，获得其表面特征的三维数字信息；

2)利用编程软件对节理面的三维数字信息进行等间距插值，并对其进行delaunay三角剖分，将结构面表面划分为若干个三角形，每个三角形称为三角形微凸体，进而可以获得三角形微凸体的集合t，如公式(1)所示，

t＝[t1(v11，v12，v13)t2(v21，v22，v23)…ti(vi1，vi2，vi3)…tn(vn1，vn2，vn3)](1)

式中ti为结构面表面上第i个三角形微凸体，vi1、vi2、vi3为第i个三角形微凸体的三个顶点，n为三角形微凸体的总数；

3)根据三角形微凸体的顶点坐标及其剪切方向的方向向量，利用编程软件计算每个三角形微凸体的视倾角的计算公式如下述公式(2)：

式中为第i个三角形微凸体的视倾角，s为剪切方向的方向方向向量，no剪切平面的法向量，ni为第i个三角形微凸体的外法向量，nio为第i个三角形微凸体的外法向量ni的在剪切平面的投影向量；

4)为区分哪些三角形微凸体具有潜在抗剪力，根据其视倾角值将其分为两类，视倾角不小于0的微凸体为具有潜在抗剪力的三角形微凸体，将其记为反之为不具有潜在抗剪力的三角形微凸体，并记为其公式如下述公式(3)：

5)找出结构面表面上所有的具有潜在抗剪力的三角形微凸体和不具有潜在抗剪力的三角形微凸体，进而可以获得三角形微凸体的分类集合at，其公式如下所诉：

6)计算每个三角形微凸体在剪切平面的投影重心坐标其中与的计算公式如下：

式中和为三角形微凸体的三个顶点坐标；

7)根据三角形微凸体的分类集合at与其在剪切平面投影的重心坐标构造结构面抗剪区域的点云矩阵st，st为n×3的矩阵，并利用编程画出st的点云图，点云图中所有具有潜在抗剪力的三角形微凸体所组成的区域为剪切过程中能够提供抗剪力的区域，其中结构面抗剪区域的点云矩阵st如下述公式所示：

进一步，所述步骤1)中，三维扫描仪为高精度三维扫描仪。

再进一步，所述步骤2)、步骤3)和步骤7)中，所述编程软件为matlab软件。

本发明的有益效果为：由于采用了上述技术方案，可以有效地识别结构面表面上的潜在抗剪区域，能够预测结构面在剪切过程中哪些区域才有可能提供剪切力，这些区域是分析结构面的变形特征和破坏特征时需要重点关注的区域。

附图说明

图1为本发明流程分析图；

图2为结构面的表面特征及剪切方向；

图3为不同剪切方向上结构面抗剪区域的点云图和破坏特征。

具体实施方式

下面结合附图详细说明本发明。

参照图1～图3，一种结构面剪切试验潜在抗剪区域的识别方法，以结构面a为研究对象，包括以下步骤：

1)利用三维扫描仪对结构面进行扫描，获得其表面特征的三维数字信息；

2)对节理面的三维数字信息进行等间距插值，并对其进行delaunay三角剖分，将结构面表面划分为若干个三角形，每个三角形称为三角形微凸体，进而可以获得三角形微凸体的集合t，如公式(1)所示，

t＝[t1(v11，v12，v13)t2(v21，v22，v23)…ti(vi1，vi2，vi3)…tn(vn1，vn2，vn3)](1)

式中ti为结构面表面上第i个三角形微凸体，vi1、vi2、vi3为第i个三角形微凸体的三个顶点，n为三角形微凸体的总数；

3)根据三角形微凸体的顶点坐标及其剪切方向的方向向量，利用编程软件计算每个三角形微凸体的视倾角的计算公式如下述公式(2)：

式中为第i个三角形微凸体的视倾角，s为剪切方向的方向方向向量，no剪切平面的法向量，ni为第i个三角形微凸体的外法向量，nio为第i个三角形微凸体的外法向量ni的在剪切平面的投影向量；

4)为区分哪些三角形微凸体具有潜在抗剪力，根据其视倾角值将其分为两类，视倾角不小于0的微凸体为具有潜在抗剪力的三角形微凸体，将其记为反之为不具有潜在抗剪力的三角形微凸体，并记为其公式如下述公式(3)：

5)找出结构面表面上所有的具有潜在抗剪力的三角形微凸体和不具有潜在抗剪力的三角形微凸体，进而可以获得三角形微凸体的分类集合at，其公式如下：

6)计算每个三角形微凸体在剪切平面的投影重心坐标其中与的计算公式如下：

式中和为三角形微凸体的三个顶点坐标；

7)根据三角形微凸体的分类集合at与其在剪切平面投影的重心坐标构造结构面抗剪区域的点云矩阵st，st为n×3的矩阵，并利用编程画出st的点云图，点云图中所有具有潜在抗剪力的三角形微凸体所组成的区域为剪切过程中能够提供抗剪力的区域，其中结构面抗剪区域的点云矩阵st如下述公式所示：

实例：一种结构面剪切试验潜在抗剪区域的识别方法，包括以下步骤：

1)以结构面a为研究对象，其形貌特征和剪切方向如图2所示，利用三维扫描仪对其进行扫描，获取结构面a的三维数字信息。

2)利用matlab软件编程对结构面a的三维数字信息进行等间距插值，并对其进行delaunay三角剖分；然后计算每个三角形微凸体的视倾角根据值的大小对结构面表面的三角形微凸体进行分类，进而得到其分类集合at。

3)计算每个三角形微凸体在剪切平面的投影重心坐标，进而得到结构面抗剪区域的点云矩阵st(114242×3的矩阵)，其中点云st矩阵部分值如表1所示。

表1

4)结构面抗剪区域的点云矩阵st绘制其点云图，如图3所示，该图中黑色区域为具有潜在抗剪力的区域。图3中还展示了该结构面在0.6mpa法向应力下剪切试验的破坏特征，对比可知结构面的破坏区域基本都落在了具有潜在抗剪力的区域，这说明在分析剪切力学特征时有潜在抗剪力的区域是需要重点关注的区域。