CN109493380B - 一种岩石节理面抗剪试验中不规则剪切面面积的计算方法 - Google Patents
一种岩石节理面抗剪试验中不规则剪切面面积的计算方法 Download PDFInfo
- Publication number
- CN109493380B CN109493380B CN201710813314.7A CN201710813314A CN109493380B CN 109493380 B CN109493380 B CN 109493380B CN 201710813314 A CN201710813314 A CN 201710813314A CN 109493380 B CN109493380 B CN 109493380B
- Authority
- CN
- China
- Prior art keywords
- point
- joint surface
- rock
- projection
- calculating
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/28—Measuring arrangements characterised by the use of optical techniques for measuring areas
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明属于岩土工程技术领域,涉及岩石材料性能测试技术方向,具体涉及一种岩石节理面抗剪试验中不规则剪切面面积的计算方法。本发明公开了一种岩石节理面抗剪试验中不规则剪切面面积的计算方法,包括以下步骤,步骤一:将待检测的岩石试样放入三维激光扫描仪进行扫描,通过三维激光扫描仪获取岩石试样的三维点云模型并提取岩石节理面点云数据;步骤二:根据所述步骤一得到的岩石节理面点云数据,采用基于散乱点的点云几何信息计算方法拟合切平面,并求出点云中各点在切平面的投影坐标;步骤三:将所述步骤二得到的投影坐标点集采用逐点插入的Delaunay三角剖分算法,由此确定投影坐标点集所构成的凸包(即最小凸多边形),其面积即为岩石节理面面积。
Description
技术领域
本发明属于岩土工程技术领域,涉及岩石材料性能测试技术方向,具体涉及一种岩石节理面抗剪试验中不规则剪切面面积的计算方法。
背景技术
在岩体边坡工程中,研究岩体不连续面的力学性质以及预想滑动剪切面的力学性质是研究岩体边坡稳定性的重要组成部分,这些不连续面和滑动剪切面的力学性质中最重要的一个力学性质就是抗剪强度,因此对岩体的抗剪强度进行研究就显得十分重要。
通过给定不同法向应力N进行多次试验,求出不同的剪应力τ,岩石的抗剪强度σ-τ曲线,从而求得岩石结构面的粘聚力c和内摩擦角由公式(1)和公式(2)可知正应力σ和剪应力τ的计算与剪切面面积A有关,因此剪切面面积A的计算精度对岩石的抗剪强度σ-τ曲线有重要影响。
由于岩石破裂方向的不确定性,以及破裂面的不规则特性,因此在进行岩体抗剪强度计算时,影响其计算结果准确度的一个重要因素便是岩体不规则剪切面面积A的计算。目前人们进行的岩石节理面抗剪试验,在测量和计算剪切面面积时都取其在水平面上的投影面积(即试件的横截面面积)来近似代替,或是对于具有一定倾角的规则剪切面(如椭圆形剪切面等),通常的做法都是用其投影面的面积除以其与水平面夹角的余弦值,而这种计算方法的前提是假设剪切面为平面,没有考虑剪切面为空间曲面这一实际情况。特别是对于岩体裂缝只沿表面分布,或是裂缝未贯穿岩体的情况下,要测试裂缝所在面与水平面的夹角就极为困难,且选择哪一段裂缝的切线作为倾斜剪切面的角度测量标准线也是一个问题。对于目前测量破裂面与水平面的夹角,大多采用量角器进行测量,该方法测量的结果可靠度较低,对于破裂面倾斜度较小的情况,人工测量难度较大。
针对上述问题,我们亟需一种岩石节理面抗剪试验中不规则剪切面面积的计算方法。针对任意破裂面,该方法能够较为精确的计算出其面积,且适用于裂缝未贯穿的情况。
发明内容
本发明的目的是:旨在提供一种岩石节理面抗剪试验中不规则剪切面面积的计算方法,该方法能够对任意破裂面进行较为精确的面积计算,且适用于裂缝未贯穿的情况。
为实现上述技术目的,本发明采用的技术方案如下:
一种岩石节理面抗剪试验中不规则剪切面面积的计算方法,包括以下步骤,
步骤一:将待检测的岩石试样放入三维激光扫描仪进行扫描,通过三维激光扫描仪获取岩石试样的三维点云模型并提取岩石节理面点云数据;
步骤二:根据所述步骤一得到的岩石节理面点云数据,采用基于散乱点的点云几何信息计算方法拟合切平面,并求出点云中各点在切平面的投影坐标;
步骤三:将所述步骤二得到的投影坐标点集采用逐点插入的Delaunay三角剖分算法,由此确定投影坐标点集所构成的凸包(即最小凸多边形),其面积即为岩石节理面面积。
作为本发明的一种优选方案,所述步骤一中获取岩石试样的三维点云模型和岩石节理面点云数据是通过三维激光扫描技术以及点云数据后处理技术实现的。
作为本发明的一种优选方案,所述步骤二中基于散乱点的点云几何信息计算方法拟合切平面并求出点云中各点在切平面的投影坐标,即通过点pi及其近邻k构造岩石节理面点集{pi,Nk(pi)},采用主元分析法(PCA)计算任意点pi的法向量最后由点pi,j∈Nk(pi)和法向量求出岩石节理面点云中各点在切平面Γ的投影坐标{(uj,vj)}。
作为本发明的一种优选方案,所述步骤二中计算投影坐标的具体计算方法如下:
c、计算点pi,j∈{pi,Nk(pi)}的投影点qi,j在O-UV系下的坐标(Ui,j,Vi,j),其计算公式如下:
其中,pi点的投影点qi对应的投影坐标为(Ui,0,Vi,0)=(0,0),即投影坐标系原点;
d、计算投影坐标的归一化坐标{(uj,vj)},其计算公式如下:
作为本发明的一种优选方案,所述步骤三中对投影坐标点集采用逐点插入的Delaunay三角剖分算法,包括以下步骤,
步骤1:构造一个极大三角形,包含待剖分点集,构建三角形链表;
步骤2:将点集中的散点依次插入,同时搜索位于外接圆内部的插入点的三角形,删除这些三角形的公共边,将插入点同三角形的全部顶点连接起来,从而完成一个散点在Delaunay三角形链表中的插入;
步骤3:按照“最小角最大化和空外接圆”优化准则对局部新形成的三角形进行优化并存入Delaunay三角形链表;
步骤4:重复步骤2~步骤3,直到所有散点全部插入完毕。
步骤5:投影坐标点集经过Delaunay三角剖分后便可确定其形成的凸包,凸包的面积即为岩石节理面面积。
与现有技术相比,本发明的有益效果是:
1、提出了一种岩石节理面抗剪试验中不规则剪切面面积的计算方法,该方法能够对任意破裂面进行较为精确的面积计算,且适用于裂缝未贯穿的情况;
2、提出了一种基于散乱点的点云几何信息方法来计算岩石节理面抗剪试验中不规则剪切面面积,并将其计算方法程序化、自动化和结果精确化。本方法使用方便,可操作性强;
3、提出了一种将测绘技术与土木工程相结合来研究岩石节理面性质的新思路,对测绘技术和土木工程相融合的多学科交叉研究发展起到了推动作用。
附图说明
本发明可以通过附图给出的非限定性实施例进一步说明;
图1为本发明一种岩石节理面抗剪试验中不规则剪切面面积的计算方法实施例中岩石试样的三维点云模型;
图2为本发明一种岩石节理面抗剪试验中不规则剪切面面积的计算方法实施例中岩石试样的不规则剪切面点云;
图3为本发明一种岩石节理面抗剪试验中不规则剪切面面积的计算方法实施例中Delaunay三角剖分优化准则“最小角最大和空外接圆”示意图;
图4为本发明一种岩石节理面抗剪试验中不规则剪切面面积的计算方法实施例中Delaunay三角剖分逐点插入示意图;
图5为本发明一种岩石节理面抗剪试验中不规则剪切面面积的计算方法实施例中Delaunay三角剖分算法流程图;
图6为本发明一种岩石节理面抗剪试验中不规则剪切面面积的计算方法实施例的流程示意图;
具体实施方式
为了使本领域的技术人员可以更好地理解本发明,下面结合附图和实施例对本发明技术方案进一步说明。
如图6所示,一种岩石节理面抗剪试验中不规则剪切面面积的计算方法,包括以下步骤,
步骤一:将待检测的岩石试样放入三维激光扫描仪进行扫描,通过三维激光扫描仪获取岩石试样的三维点云模型并提取岩石节理面点云数据,分别如图1和图2所示;
步骤二:根据步骤一得到的岩石节理面点云数据,采用基于散乱点的点云几何信息计算方法拟合切平面,并求出点云中各点在切平面的投影坐标;即通过点pi及其近邻k构造岩石节理面点集{pi,Nk(pi)},采用主元分析法(PCA)计算任意点pi的法向量最后由点pi,j∈Nk(pi)和法向量求出岩石节理面点云中各点在切平面Γ的投影坐标{(uj,vj)}。
计算投影坐标的具体计算方法如下:
c、计算点pi,j∈{pi,Nk(pi)}的投影点qi,j在O-UV系下的坐标(Ui,j,Vi,j),其计算公式如下:
其中,pi点的投影点qi对应的投影坐标为(Ui,0,Vi,0)=(0,0),即投影坐标系原点;
d、计算投影坐标的归一化坐标{(uj,vj)},其计算公式如下:
步骤三:将步骤二得到的投影坐标点集采用逐点插入的Delaunay三角剖分算法,由此确定投影坐标点集所构成的凸包(即最小凸多边形),其面积即为岩石节理面面积。
根据Delaunay三角剖分算法流程图5,借助Matlab编程实现对投影坐标点集{(uj,vj)}的Delaunay三角剖分算法,由此确定投影坐标点集所构成的凸包(即最小凸多边形),其面积即为岩石节理面面积。其中Delaunay三角剖分剖分优化准则“最小角最大和空外接圆”示意图如图3所示;Delaunay三角剖分逐点插入示意图如图4所示。
上述实施例仅示例性说明本发明的原理及其功效,而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下,对上述实施例进行修饰或改变。因此,凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变,仍应由本发明的权利要求所涵盖。
Claims (2)
1.一种岩石节理面抗剪试验中不规则剪切面面积的计算方法,其特征在于:包括以下步骤,
步骤一:将待检测的岩石试样放入三维激光扫描仪进行扫描,通过三维激光扫描仪获取岩石试样的三维点云模型并提取岩石节理面点云数据;
步骤二:根据所述步骤一得到的岩石节理面点云数据,采用基于散乱点的点云几何信息计算方法拟合切平面,并求出点云中各点在切平面的投影坐标;
步骤三:将所述步骤二得到的投影坐标点集采用逐点插入的Delaunay三角剖分算法,由此确定投影坐标点集所构成的凸包,其面积即为岩石节理面面积;
所述步骤一中获取岩石试样的三维点云模型和岩石节理面点云数据是通过三维激光扫描技术以及点云数据后处理技术实现的;
所述步骤二中基于散乱点的点云几何信息计算方法拟合切平面并求出点云中各点在切平面的投影坐标,即通过点pi及其近邻k构造岩石节理面点集{pi,Nk(pi)},采用主元分析法计算任意点pi的法向量最后由点pi,j∈Nk(pi)和法向量求出岩石节理面点云中各点在切平面Γ的投影坐标{(uj,vj)};
所述步骤二中计算投影坐标的具体计算方法如下:
c、计算点pi,j∈Nk(pi)的投影点qi,j在O-UV系下的坐标(Ui,j,Vi,j),其计算公式如下:
其中,pi点的投影点qi对应的投影坐标为(Ui,0,Vi,0)=(0,0),即投影坐标系原点;
d、计算投影坐标的归一化坐标{(uj,vj)},其计算公式如下:
2.根据权利要求1所述的一种岩石节理面抗剪试验中不规则剪切面面积的计算方法,其特征在于:所述步骤三中对投影坐标点集采用逐点插入的Delaunay三角剖分算法,包括以下步骤,
步骤1:构造一个极大三角形,包含待剖分点集,构建三角形链表;
步骤2:将点集中的散点依次插入,同时搜索位于外接圆内部的插入点的三角形,删除这些三角形的公共边,将插入点同三角形的全部顶点连接起来,从而完成一个散点在Delaunay三角形链表中的插入;
步骤3:按照“最小角最大化和空外接圆”优化准则对局部新形成的三角形进行优化并存入Delaunay三角形链表;
步骤4:重复步骤2~步骤3,直到所有散点全部插入完毕;
步骤5:投影坐标点集经过Delaunay三角剖分后便可确定其形成的凸包,凸包的面积即为岩石节理面面积。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710813314.7A CN109493380B (zh) | 2017-09-11 | 2017-09-11 | 一种岩石节理面抗剪试验中不规则剪切面面积的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710813314.7A CN109493380B (zh) | 2017-09-11 | 2017-09-11 | 一种岩石节理面抗剪试验中不规则剪切面面积的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109493380A CN109493380A (zh) | 2019-03-19 |
CN109493380B true CN109493380B (zh) | 2022-05-17 |
Family
ID=65687727
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710813314.7A Active CN109493380B (zh) | 2017-09-11 | 2017-09-11 | 一种岩石节理面抗剪试验中不规则剪切面面积的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109493380B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111127312B (zh) * | 2019-12-25 | 2023-08-22 | 武汉理工大学 | 一种复杂物体点云提取圆的方法及扫描装置 |
CN112347901A (zh) * | 2020-11-04 | 2021-02-09 | 长江岩土工程总公司(武汉) | 一种基于三维激光扫描技术岩体块体分析方法 |
CN113435087B (zh) * | 2021-06-24 | 2022-04-22 | 中铁二院工程集团有限责任公司 | 一种分析洞穴围岩局部稳定性的方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104318622A (zh) * | 2014-10-25 | 2015-01-28 | 大连理工大学 | 一种室内场景非均匀三维点云数据的三角网格建模方法 |
CN106803281A (zh) * | 2017-02-17 | 2017-06-06 | 武汉大学 | 一种开挖边坡岩块信息提取与三维重构方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102930091B (zh) * | 2012-10-24 | 2015-05-20 | 浙江大学 | 一种面向牙齿牙龈的高效网格修补方法 |
CN103985155B (zh) * | 2014-05-14 | 2017-01-25 | 北京理工大学 | 基于映射法的散乱点云Delaunay三角剖分曲面重构方法 |
US20170078593A1 (en) * | 2015-09-16 | 2017-03-16 | Indoor Reality | 3d spherical image system |
CN105628581B (zh) * | 2015-12-24 | 2018-07-03 | 中国石油大学(华东) | 一种基于高光谱技术的致密砂岩储层露头孔隙度表征方法 |
CN106227923B (zh) * | 2016-07-14 | 2018-02-27 | 中国地质大学(武汉) | 一种岩体结构面三维粗糙度评价方法 |
-
2017
- 2017-09-11 CN CN201710813314.7A patent/CN109493380B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104318622A (zh) * | 2014-10-25 | 2015-01-28 | 大连理工大学 | 一种室内场景非均匀三维点云数据的三角网格建模方法 |
CN106803281A (zh) * | 2017-02-17 | 2017-06-06 | 武汉大学 | 一种开挖边坡岩块信息提取与三维重构方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109493380A (zh) | 2019-03-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109493380B (zh) | 一种岩石节理面抗剪试验中不规则剪切面面积的计算方法 | |
Valença et al. | Automatic crack monitoring using photogrammetry and image processing | |
Marani et al. | A modified iterative closest point algorithm for 3D point cloud registration | |
CN102096944A (zh) | 地质体结构面三维激光扫描点云识别方法 | |
CN106844960B (zh) | 一种管状焊接钢结构网架的逆向建模的方法 | |
CN104848838A (zh) | 两种构形条件下岩土试样剪切带倾角演变规律的观测方法 | |
CN103542820A (zh) | 一种检测风洞内表面平整度的方法 | |
Yan et al. | Automated structural modelling of bridges from laser scanning | |
CN104851093A (zh) | 基于距离与角度控制的矢量面数据接边检测方法 | |
Wei et al. | Automated inspection and monitoring of member deformation in grid structures | |
CN115546116A (zh) | 全覆盖式岩体不连续面提取与间距计算方法及系统 | |
Herráez et al. | Optimal modelling of buildings through simultaneous automatic simplifications of point clouds obtained with a laser scanner | |
WO2020258426A1 (zh) | 一种考虑结构面粗糙度的裂隙岩体稳定性分析系统 | |
Uotinen et al. | Photogrammetric prediction of rock fracture properties and validation with metric shear tests | |
Zheng | Crack detection and measurement utilizing image-based reconstruction | |
Sun et al. | Rock mass discontinuity extraction method based on multiresolution supervoxel segmentation of point cloud | |
Lange et al. | Crack detection at concrete construction units from photogrammetric data using image processing procedures | |
WO2023155228A1 (zh) | 矿用三维彩色点云重建系统、方法 | |
KR20200064421A (ko) | 3차원 파이프 라인 검출 자동화 장치 | |
CN102564383B (zh) | 基于轴面角确定三轴磁传感器轴间角的方法 | |
CN111982152B (zh) | 点云地图量化方法、装置、计算机设备和存储介质 | |
CN111721268B (zh) | 一种精确反演建筑物高度的方法与装置 | |
Habib et al. | Error budget of lidar systems and quality control of the derived point cloud | |
CN112785135A (zh) | 工程质量检验方法、装置、计算机设备和存储介质 | |
Chen et al. | Automatic approach for measuring deformations in complex structures using photogrammetry technique |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |