结构面抗剪区域的确定方法技术

一种结构面抗剪区域的确定方法，利用三维扫描仪对结构面进行扫描，获得其表面特征的三维数字信息；利用编程软件对节理面的三维数字信息进行等间距插值，并对其进行Delaunay三角剖分；根据三角形微凸体的顶点坐标及其剪切方向的方向向量，利用编程软件计算每个三角形微凸体的视倾角；为区分哪些三角形微凸体具有潜在抗剪力，根据其视倾角值将其分为两类；找出结构面表面上所有的具有潜在抗剪力的三角形微凸体和不具有潜在抗剪力的三角形微凸体；利用编程画出S

Determination of shear zone of structural plane

【技术实现步骤摘要】
结构面抗剪区域的确定方法
本专利技术属于岩石力学工程领域，涉及一种结构面抗剪区域的确定方法。
技术介绍
大量的工程实践证实，很多工程岩体的失稳问题都与结构面存在有关，研究结构面的力学行为具有十分重要的工程意义。在结构面的剪切试验或原位抗剪试验中，确定其表面上能够提供抗剪力的区域对分析结构面的变形特征、破坏特征以及剪切机理是十分重要的。学者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。...

【技术保护点】
1.一种结构面抗剪区域的确定方法，其特征在于，所述方法包括以下步骤：/n1)利用三维扫描仪对结构面进行扫描，获得其表面特征的三维数字信息；/n2)利用编程软件对节理面的三维数字信息进行等间距插值，并对其进行Delaunay三角剖分，将结构面表面划分为若干个三角形，每个三角形称为三角形微凸体，进而可以获得三角形微凸体的集合T，如公式(1)所示，/nT＝[T

【技术特征摘要】
1.一种结构面抗剪区域的确定方法，其特征在于，所述方法包括以下步骤：
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的...

【专利技术属性】
技术研发人员：宋磊博，江权，罗战友，雍睿，张贺，杜时贵，
申请(专利权)人：绍兴文理学院，
类型：发明
国别省市：浙江;33